第六讲空间杆系分析

1、1第六讲第六讲 空间杆系结构有限元分析空间杆系结构有限元分析 杆系结构定义：杆系结构定义： 由杆件组成的空间结构称杆系结构。 杆件指细长的构件，构件的几何特征是一个方向的尺寸远大于其它两个方向的尺寸。 杆系结构根据杆件两端与其它构件连接方式的不同，可以分为桁架和刚架。桁架(杆)由两端铰接的杆件构成的空间结构。刚架(梁)亦称为框架，指由两端固接的杆件构成的空间结构。26-1 空间桁架结构有限元分析（桁架单元）空间桁架结构有限元分析（桁架单元）eiueivejuejvejweiweiqejqijYxXZ xu 由于两端铰接，不能承受力矩，假定每个节点三个自由度，在整体参考系有： jjjiiiewv

2、uwvu 为方便研究，在单元上取一个局部参考系，以 为原点，由 点指向 点为 轴的正向。iijx在局部参考系中有： ejeieqqq一、桁架单元刚度矩阵一、桁架单元刚度矩阵3ijx xu在单元中假定：bxaxu)( 利用节点处的位移连续性条件可以解出 和 ，代回原式有： eejeieeqNqqLxLxxu 1ab其中单元长度为：222ijijijeZZYYXXLeeeeqLxLxxqNxxxu1)(ejeiejeieeqqBqqLL1,1eiqejq4 dxBDBAdvBDBkeeLTeTve11111,111eeeeeeeeeeLEALLLLELA其中 为单元截面积， 为单元材料杨氏模量。e

3、AeE 需要指出，以上单元刚度是基于局部参考系得到的，由于不同的单元，其局部参考系的方位各不相同。因而，局部参考系中的单元节点位移向量的方向也各不相同，局部参考系中的刚度矩阵不具有可累加性。 这就要求我们将局部参考系下的单元刚度矩阵变换到整体参考系中。5二、坐标变换矩阵二、坐标变换矩阵局部参考系与整体参考系中单元节点位移向量的变换关系为：ewejeiejeieieieiwnvmulqwnvmulq写成矩阵形式为：jjjiiuejeiwvuwvunmlnmlqq000000简记为： eeq其中：eijLXXleijLYYmeijLZZn6单元弹性应变能在整体参考系中的表达式为： eeTeeK21

4、单元弹性应变能在局部参考系中的表达式为： eeTeeqkq21 eeq将： 代入上式有： eeTTeeeTeekqkq2121将该式与整体参考系中的单元弹性应变能对比可知： eTekK7三、等效单元体力载荷向量三、等效单元体力载荷向量桁架单元的等效体力载荷向量可以简单地表示为： ),(),(),(10001000110001000121cccbzcccbycccbxeeebZYXFZYXFZYXFLAPcccZYX,其中： 为单元形心坐标 ，可以由下式给出： jicXXX21jicYYY21jicZZZ218四、单元应力计算四、单元应力计算 eeeeeBEqBEE96-2 平面刚架（框架）结构

5、有限元分析平面刚架（框架）结构有限元分析一、平面刚架单元受力分析与变形分解一、平面刚架单元受力分析与变形分解XYjiyxeiueieivejejvejueq1eq3eq2eq5eq4eq6 该单元可以承受 和 方向的力以及绕 轴方向的力矩。XYZ 在整体参考系中单元节点位移向量为： ejejejeieieievuvu 为方便研究，在单元上取一个局部参考系，以 为原点，由 点指向 点为 轴的正象。iijx10yxeq1eq3eq2eq5eq4eq6)(xv)(xu 在局部参考系中单元所受载荷可以分解为：1）轴向拉压，主要引起单元的轴向变形。2）横向载荷及绕平面法向的弯矩，主要引起单元的弯曲。 为

6、研究方便，我们在做刚架结构有限元分析时，也将单元变形分解为轴向变形与弯曲变形两部分，分别进行研究。 在小变形条件下，可以认为单元的轴向变形与弯曲变形是相互独立的并符合迭加原理迭加原理。因此，我们可以分别计算两种基本变形的刚度，在进行迭加。 刚架结构的有限元分析在研究弯曲变形时采用了材料力学梁的基本假定，因此，有时也将刚架单元称为梁单元。111）单元轴向变形刚度矩阵）单元轴向变形刚度矩阵yxeiqejq)(xu 与轴向变形相关的节点位移分量有 和 ，假定在轴向载荷作用下单元内单元内任意点的位移为：eiqejqexejeeeqNqqLxLxxu11)(与轴向变形对应的刚度矩阵为： 1111eeee

7、xLEAkeiqeiqejqejqji二、平面刚架单元两种基本变形的的刚度矩阵二、平面刚架单元两种基本变形的的刚度矩阵122）单元弯曲变形刚度）单元弯曲变形刚度yxeieivejvej)(xv 与弯曲变形相关的节点位移分量有 、 、 和 。eiv假定在轴向载荷作用下单元内单元内任意点的挠度为：中性层上的362543)(xxxxv根据节点处的位移连续性条件有：6 , 5 , 4 , 3,ii解出 ， 代回挠度函数经整理可得：在 点：eivv)0(jiji挠度的导数eixxv00 x在 点：ejevLv)(ejLxexveLx 小变性、中心平面、eiejvej13挠度函数为：其中：单元中性层外任意

8、点沿局部参考系 方向得位移可以表示为：x与假定坐标方向相隔转角方向xxvyyxu)(),(yx),(yx),(yxuxy exyNxv)( 4321NNNNN 223432332223233231/32/2/32LLxxxNLLxxxNLxLLxxxNLLLxxxN ejejeieiexyvv14 xyxxBxxvyxu22单元中性层外任意点沿局部参考系 方向得应变可以表示为：x ByLxLLxxLLxLyB26612466123 222232046266126122646612612LLLLLLLLLLLLLEIDEDdAydxBBEdVBDBkxxALTTVexye弹性矩阵在这里其中：在局

9、部参考系中单元刚度矩阵可以表示为其中：dAyIAzz2绕 轴截面弯曲惯性矩yejejve3eiveive3ejvej15将轴向变形与弯曲变形的刚度迭加，可以得到平面刚架单元的刚度矩阵。 LEILEILEILEILEILEILEILEILEALEALEILEILEILEILEILEILEILEILEALEAkzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzze460260612061200000260460612061200000222323222323eiueieivejuejvejeiueieivejuejvej166-3 空间刚架（框架）结构有限元分析空间刚架（框架）结构有限

10、元分析XYjiZeQ1eQ3eQ2eQ6eQ5eQ4eQ7eQ8eQ9eQ10eQ11eQ12一、空间刚架单元受力分析与变形分解一、空间刚架单元受力分析与变形分解 该单元可以承受三个方向的力以及绕三个轴的力矩。 在整体参考系中单元节点位移向量为：eeeeeQQQQQ12321 为方便研究，在单元上取一个局部参考系，以 为原点，由 点指向 点为 轴。以截面的两个惯性主轴分别为局部参考系的 轴和 轴。iijxyz17xyjizeq7eq9eq8eq12eq11eq10eq1eq2eq3eq4eq5eq61) 在轴向载荷的作用下，将产生轴向拉压变形，与之相关的节点位移分量为： eeexqqq712)

11、 在绕轴向载荷（扭矩）的作用下，将产生扭转变形，与之相关的节点位移分量为： eeetqqq1043) 在 向载荷及绕轴 力矩的作用下，将产生 平面内的弯曲变形，与之相关的节点位移分量为：zxy eeeeexyqqqqq12862y4) 在 向载荷及绕轴 力矩的作用下，将产生 平面内的弯曲变形，与之相关的节点位移分量为：yxzz eeeeexzqqqqq1195318二、空间刚架单元四种基本变形的刚度矩阵二、空间刚架单元四种基本变形的刚度矩阵1）单元轴向变形刚度矩阵）单元轴向变形刚度矩阵yxeq1eq7)(xu 与轴向变形相关的节点位移分量有 和 ，假定在轴向载荷作用下单元内单元内任意点的位移为

12、：eq1eq7 exeeeeqNqqLxLxxu711)(与轴向变形对应的刚度矩阵为： 1111eeeexLEAKeq1eq1eq7eq7ji192）绕单元轴向扭转变形刚度矩阵）绕单元轴向扭转变形刚度矩阵yxeq4eq10)(x 与扭转变形相关的节点位移分量有 和 ，假定在轴向载荷作用下单元内单元内任意点的扭转角为：eq4eq10 eteeeeqNqqLxLxx1041)(ji在单元中假定：bxax)(在 点：eq4)0(ji0 x在 点：eeqL10)(eLx 解出 和 代回插值函数有：ab eeeexqqLrLrdxdr104 eeLrLrB对于圆形截面有：其几何矩阵为：20与轴向变形对应

13、的刚度矩阵为： eq4eq4eq10eq10 xxG GD dArdxLLLLGdVBDBkALeeeeTVete201111 1111eetLGJk根据应力应变关系可以得到则扭转变形刚度矩阵为：dArJA2令：圆形截面的扭转惯性矩。说明：说明： 对于非圆截面前面的应变位移关系及扭转惯性矩算式 不成立，但刚度矩阵是普遍适用的，只是扭转惯性矩需要用其它方法获得，如：近似公式计算、实验等。213）单元）单元 平面内弯曲变形刚度矩阵平面内弯曲变形刚度矩阵yxeq6eq2eq8eq12)(xv 与弯曲变形相关的节点位移分量有 、 、 和 。eq6eq2eq8eq12假定在横向载荷和弯矩作用下单元内单元

14、内任意点的挠度为：362543)(xxxxv根据节点处的位移连续性条件有：6 , 5 , 4 , 3,ii解出 ， 代回挠度函数经整理可得：在 点：eqv2)0(jijiexqxv600 x在 点：eeqLv8)(eLxqxve12eLx xy22挠度函数为：其中：单元中性层外任意点沿局部参考系 方向得位移可以表示为：xxxvyyxuxy)(),(yx),(yx),(yxuxy 223432332223232231/32/2/32LLxxxNLLxxxNLxLLxxxNLLLxxxN eeeeexyqqqqq12862 4321NNNNN exyqNxv)(23 xyxyxxqBxxvyxyx

15、u22,单元中性层外任意点沿局部参考系 方向得应变可以表示为：x ByLxLLxxLLxLyB26612466123 222232046266126122646612612LLLLLLLLLLLLLEIdAydxBBEdVBDBkxxALTTVexye其中：在局部参考系中单元刚度矩阵可以表示为其中：dAyIAzz2绕 轴截面惯性矩yeq12eq8eq6eq2eq6eq8eq12eq2244）单元）单元 平面内弯曲变形刚度矩阵平面内弯曲变形刚度矩阵zxeq5eq3eq9eq11)(xw 与弯曲变形相关的节点位移分量有 、 、 和 。eq5eq3eq9eq11假定在单元内单元内任意点的挠度为：31

16、02987)(xxxxw根据节点处的位移连续性条件有：在 点：eqv3)0(jijiexqxv500 x在 点：eeqLv9)(eLxqxve11eLx xz 2222346266126122646612612LLLLLLLLLLLLLEIkyyexz其中：dAzIAyy2绕 轴截面惯性矩y对应的单元刚度矩阵为：eq5eq3eq9eq11eq5eq3eq9eq11其他过程与 平面内的弯曲相同，不再xy25 LEILEILEILEILEILEILEILEILGJLGJLEILEILEILEILEALEALEILEILEILEILGJLEILEILEAkzzzzzzzzyyyyyyyyyyyyzz

17、zzzzzzyyyyyyzze400060200060406000206000000000012000001200120000012000000400060406000001200120033333333333333对称空间刚架单元在局部参考系中的单元刚度矩阵为：26三、空间刚架单元刚度矩阵的变换关系三、空间刚架单元刚度矩阵的变换关系整体参考系1）空间刚架单元的节点位移向量局部参考系eeeeeQQQQQ12321 eeeeeqqqqq123212）空间刚架单元的节点位移向量的变换关系 eeQq 坐标变换矩阵，具体形式如下：27 1212000000000000 333222111nmlnmln

18、ml222nml111nml333nml其中：局部参考系 轴关于整体参考系的方向余弦局部参考系 轴关于整体参考系的方向余弦局部参考系 轴关于整体参考系的方向余弦xyz 因为空间刚架单元的局部参考系有三个坐标轴，根据单元的两个节点只能确定其中一个坐标轴的方向，为确定其它两个坐标轴的方向，还需要因入一个辅助量，可以是一个节点或一个参考方向。28四、空间刚架单元局部参考系确定方法（参考点法）四、空间刚架单元局部参考系确定方法（参考点法）jzxyikz 所谓参考点法，是在单元截面主轴平面内找一个辅助点 （称参考点），用以帮助确定单元局部参考系。k首先由 点到 点作一个向量 AijTzyxAAAA,ijzijy