理论教育 杆系结构有限元:结构坐标系中的单元刚度矩阵

杆系结构有限元:结构坐标系中的单元刚度矩阵

时间:2023-08-23 理论教育 版权反馈
【摘要】:单元坐标系与结构坐标系间的关系较规整时,也可按概念直接推导结构坐标系统下的单元刚度矩阵元素。图2.8单元刚度矩阵元素对于单元①:α=90°,若按式计算,执行矩阵乘法计算量较大,故根据单元刚度元素kij的性质进行元素推导,如图2.9所示。图2.9结构整体坐标系下的杆端自由度kij的意义为:在结构坐标系下,杆端自由度j发生单位位移,在i自由度方向上需作用力的大小。

杆系结构有限元:结构坐标系中的单元刚度矩阵

1)截面力向量的坐标变换

平面杆系结构任一截面力向量,需要以3个分量描述其构成元素。此时存在两个坐标系统:单元坐标系统和结构坐标系统(X-Y-θ)。规定由X轴到x轴的夹角α以顺时针旋转为正,如图2.6所示。

图2.6 截面力在不同坐标系统下的分量表示

由式(1.9),并结合考虑平面体系中坐标系旋转变换与转动自由度无关,故有

写成矩阵形式,为

即在不同直角坐标系之间,截面力向量使用坐标变换矩阵T0进行坐标变换。

显然,T0仍为一个正交矩阵,此时存在的线性变换为:

2)杆端力向量的坐标变换

杆单元杆端力向量由杆件的始端杆端力和末端杆端力构成,共包含6个分量,如图2.7所示。

图2.7 单元坐标变换

始端杆端力和末端杆端力在坐标系之间,分别满足以上变换式(2.5)。因此,单元坐标系下的杆端力,结构坐标系下的杆端力Fe若分别表示为:

则根据式(2.5),应有以下坐标变换成立:

或展开写为

或简写成

式(2.6)中

称为一般平面杆单元的坐标转换矩阵。T仍为正交矩阵,因而有

3)杆端位移向量的坐标变换

杆端力向量的坐标变换方式,也同样适用于杆端位移向量之间的坐标变换,即

4)单元刚度矩阵的坐标变换(www.daowen.com)

将式(2.9)代入,有:

式(a)左右同时左乘TT,可得:

即结构坐标系下的杆端位移向量与杆端力向量之间的变换关系,由3个矩阵的乘积决定,式(b)可写为:

此即为结构坐标系中的单元刚度方程。其中

式中,Ke为结构坐标系中的单元刚度矩阵

根据矩阵知识可知,为相似变换,实对称矩阵经此变换后仍应为对称性不变,矩阵特征值也不会改变。因此,结构坐标系统下的单元刚度矩阵Ke仍然是对称矩阵;Ke仍然是奇异矩阵。由于执行了坐标变换,Ke除与单元本身的属性有关外,还与结构坐标系与单元坐标系之间的夹角α有关。

Ke中任一元素表示结构坐标系下杆端位移δe中第m个分量等于1(其余位移分量均等于零)时,所引起的杆端力Fe的第l个分量的值。单元坐标系与结构坐标系间的关系较规整时(如α=0°或90°),也可按概念直接推导结构坐标系统下的单元刚度矩阵元素。

【例2.1】试以单元刚度系数的概念,直接写出如图2.8所示结构各单元在结构坐标系统下的单元刚度矩阵。各杆EI、EA为常数,结构坐标系和单元坐标系如图2.8所示。

图2.8 单元刚度矩阵元素

【解】(1)对于单元①:

α=90°,若按式(2.12)计算,执行矩阵乘法计算量较大,故根据单元刚度元素kij的性质进行元素推导,如图2.9所示。

图2.9 结构整体坐标系下的杆端自由度

kij的意义为:在结构坐标系下,杆端自由度j发生单位位移,在i自由度方向上需作用力的大小。

令杆AB在沿1方向发生水平单位位移,确定对应杆端力,由形常数即有:

令杆AB在沿2方向发生竖向单位位移:

令杆AB在沿3方向发生单位转角:

令杆AB在沿4方向发生水平单位位移:

令杆AB在沿5方向发生竖向单位位移:

令杆AB在沿6方向发生单位转角:

可得单元①的单元刚度矩阵为:

(2)对于单元②:显然,单元②的单元刚度矩阵与单元①的单元刚度矩阵相同,有:

(3)对于单元③:α=0°,单元刚度矩阵在两个坐标系统中是相同的,即式(2.4)的形式,此处不再列写。

免责声明:以上内容源自网络,版权归原作者所有,如有侵犯您的原创版权请告知,我们将尽快删除相关内容。

我要反馈