第六讲空间杆系分析

1第六讲空间杆系结构有限元分析杆系结构定义：由杆件组成的空间结构称杆系结构。杆件指细长的构件，构件的几何特征是一个方向的尺寸远大于其它两个方向的尺寸。杆系结构根据杆件两端与其它构件连接方式的不同，可以分为桁架和刚架。桁架(杆)——由两端铰接的杆件构成的空间结构。刚架(梁)——亦称为框架，指由两端固接的杆件构成的空间结构。2§6-1空间桁架结构有限元分析（桁架单元）eiueivejuejvejweiweiqejqijYxXZxu由于两端铰接，不能承受力矩，假定每个节点三个自由度，在整体参考系有：jjjiiiewvuwvu为方便研究，在单元上取一个局部参考系，以为原点，由点指向点为轴的正向。iijx在局部参考系中有：ejeieqqq一、桁架单元刚度矩阵3ijxxu在单元中假定：bxaxu)(利用节点处的位移连续性条件可以解出和，代回原式有：eejeieeqNqqLxLxxu1ab其中单元长度为：222ijijijeZZYYXXLeeeeqLxLxxqNxxxu1)(ejeiejeieeqqBqqLL1,1eiqejq4dxBDBAdvBDBkeeLTeTve11111,111eeeeeeeeeeLEALLLLELA其中为单元截面积，为单元材料杨氏模量。eAeE需要指出，以上单元刚度是基于局部参考系得到的，由于不同的单元，其局部参考系的方位各不相同。因而，局部参考系中的单元节点位移向量的方向也各不相同，局部参考系中的刚度矩阵不具有可累加性。这就要求我们将局部参考系下的单元刚度矩阵变换到整体参考系中。5二、坐标变换矩阵局部参考系与整体参考系中单元节点位移向量的变换关系为：ewejeiejeieieieiwnvmulqwnvmulq写成矩阵形式为：jjjiiuejeiwvuwvunmlnmlqq000000简记为：eeq其中：eijLXXleijLYYmeijLZZn6单元弹性应变能在整体参考系中的表达式为：eeTeeK21单元弹性应变能在局部参考系中的表达式为：eeTeeqkq21eeq将：代入上式有：eeTTeeeTeekqkq2121将该式与整体参考系中的单元弹性应变能对比可知：eTekK7三、等效单元体力载荷向量桁架单元的等效体力载荷向量可以简单地表示为：),,(),,(),,(10001000110001000121cccbzcccbycccbxeeebZYXFZYXFZYXFLAPcccZYX,,其中：为单元形心坐标，可以由下式给出：jicXXX21jicYYY21jicZZZ218四、单元应力计算eeeeeBEqBEE9§6-2平面刚架（框架）结构有限元分析一、平面刚架单元受力分析与变形分解XYjiyxeiueieivejejvejueq1eq3eq2eq5eq4eq6该单元可以承受和方向的力以及绕轴方向的力矩。XYZ在整体参考系中单元节点位移向量为：ejejejeieieievuvu为方便研究，在单元上取一个局部参考系，以为原点，由点指向点为轴的正象。iijx10yxeq1eq3eq2eq5eq4eq6)(xv)(xu在局部参考系中单元所受载荷可以分解为：1）轴向拉压，主要引起单元的轴向变形。2）横向载荷及绕平面法向的弯矩，主要引起单元的弯曲。为研究方便，我们在做刚架结构有限元分析时，也将单元变形分解为轴向变形与弯曲变形两部分，分别进行研究。在小变形条件下，可以认为单元的轴向变形与弯曲变形是相互独立的并符合迭加原理。因此，我们可以分别计算两种基本变形的刚度，在进行迭加。刚架结构的有限元分析在研究弯曲变形时采用了材料力学梁的基本假定，因此，有时也将刚架单元称为梁单元。111）单元轴向变形刚度矩阵yxeiqejq)(xu与轴向变形相关的节点位移分量有和，假定在轴向载荷作用下单元内单元内任意点的位移为：eiqejqexejeeeqNqqLxLxxu11)(与轴向变形对应的刚度矩阵为：1111eeeexLEAkeiqeiqejqejqji二、平面刚架单元两种基本变形的的刚度矩阵122）单元弯曲变形刚度yxeieivejvej)(xv与弯曲变形相关的节点位移分量有、、和。eiv假定在轴向载荷作用下单元内单元内任意点的挠度为：中性层上的362543)(xxxxv根据节点处的位移连续性条件有：6,5,4,3,ii解出，代回挠度函数经整理可得：在点：eivv)0(jiji挠度的导数eixxv00x在点：ejevLv)(ejLxexveLx小变性、中心平面、eiejvej13挠度函数为：其中：单元中性层外任意点沿局部参考系方向得位移可以表示为：x与假定坐标方向相隔转角方向xxvyyxu)(),(yx),(yx),(yxuxyexyNxv)(4321NNNNN223432332223233231//32/2/32LLxxxNLLxxxNLxLLxxxNLLLxxxNejejeieiexyvv14xyxxBxxvyxu22单元中性层外任意点沿局部参考系方向得应变可以表示为：xByLxLLxxLLxLyB26612466123222232046266126122646612612LLLLLLLLLLLLLEIDEDdAydxBBEdVBDBkxxALTTVexye弹性矩阵在这里其中：在局部参考系中单元刚度矩阵可以表示为其中：dAyIAzz2——绕轴截面弯曲惯性矩yejejve3eiveive3ejvej15将轴向变形与弯曲变形的刚度迭加，可以得到平面刚架单元的刚度矩阵。LEILEILEILEILEILEILEILEILEALEALEILEILEILEILEILEILEILEILEALEAkzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzze460260612061200000260460612061200000222323222323eiueieivejuejvejeiueieivejuejvej16§6-3空间刚架（框架）结构有限元分析XYjiZeQ1eQ3eQ2eQ6eQ5eQ4eQ7eQ8eQ9eQ10eQ11eQ12一、空间刚架单元受力分析与变形分解该单元可以承受三个方向的力以及绕三个轴的力矩。在整体参考系中单元节点位移向量为：eeeeeQQQQQ12321为方便研究，在单元上取一个局部参考系，以为原点，由点指向点为轴。以截面的两个惯性主轴分别为局部参考系的轴和轴。iijxyz17xyjizeq7eq9eq8eq12eq11eq10eq1eq2eq3eq4eq5eq61)在轴向载荷的作用下，将产生轴向拉压变形，与之相关的节点位移分量为：eeexqqq712)在绕轴向载荷（扭矩）的作用下，将产生扭转变形，与之相关的节点位移分量为：eeetqqq1043)在向载荷及绕轴力矩的作用下，将产生平面内的弯曲变形，与之相关的节点位移分量为：zxyeeeeexyqqqqq12862y4)在向载荷及绕轴力矩的作用下，将产生平面内的弯曲变形，与之相关的节点位移分量为：yxzzeeeeexzqqqqq1195318二、空间刚架单元四种基本变形的刚度矩阵1）单元轴向变形刚度矩阵yxeq1eq7)(xu与轴向变形相关的节点位移分量有和，假定在轴向载荷作用下单元内单元内任意点的位移为：eq1eq7exeeeeqNqqLxLxxu711)(与轴向变形对应的刚度矩阵为：1111eeeexLEAKeq1eq1eq7eq7ji192）绕单元轴向扭转变形刚度矩阵yxeq4eq10)(x与扭转变形相关的节点位移分量有和，假定在轴向载荷作用下单元内单元内任意点的扭转角为：eq4eq10eteeeeqNqqLxLxx1041)(ji在单元中假定：bxax)(在点：eq4)0(ji0x在点：eeqL10)(eLx解出和代回插值函数有：abeeeexqqLrLrdxdr104eeLrLrB对于圆形截面有：其几何矩阵为：20与轴向变形对应的刚度矩阵为：eq4eq4eq10eq10xxGGDdArdxLLLLGdVBDBkALeeeeTVete2011111111eetLGJk根据应力应变关系可以得到则扭转变形刚度矩阵为：dArJA2令：——圆形截面的扭转惯性矩。说明：对于非圆截面前面的应变位移关系及扭转惯性矩算式不成立，但刚度矩阵是普遍适用的，只是扭转惯性矩需要用其它方法获得，如：近似公式计算、实验等。213）单元平面内弯曲变形刚度矩阵yxeq6eq2eq8eq12)(xv与弯曲变形相关的节点位移分量有、、和。eq6eq2eq8eq12假定在横向载荷和弯矩作用下单元内单元内任意点的挠度为：362543)(xxxxv根据节点处的位移连续性条件有：6,5,4,3,ii解出，代回挠度函数经整理可得：在点：eqv2)0(jijiexqxv600x在点：eeqLv8)(eLxqxve12eLxxy22挠度函数为：其中：单元中性层外任意点沿局部参考系方向得位移可以表示为：xxxvyyxuxy)(),(yx),(yx),(yxuxy223432332223232231//32/2/32LLxxxNLLxxxNLxLLxxxNLLLxxxNeeeeexyqqqqq128624321NNNNNexyqNxv)(23xyxyxxqBxxvyxyxu22,单元中性层外任意点沿局部参考系方向得应变可以表示为：xByLxLLxxLLxLyB26612466123其中：在局部参考系中单元刚度矩阵可以表示为其中：dAyIAzz2——绕轴截面惯性矩yeq12eq8eq6eq2eq6eq8eq12eq2244）单元平面内弯曲变形刚度矩阵z