第六讲空间杆系分析

第六讲空间杆系结构有限元分析

杆系结构定义：

由杆件组成的空间结构称杆系结构。杆件指细长的构件，构件的几何特征是一个方向的尺寸远大于其它两个方向的尺寸。杆系结构根据杆件两端与其它构件连接方式的不同，可以分为桁架和刚架。桁架(杆)——由两端铰接的杆件构成的空间结构。刚架(梁)——亦称为框架，指由两端固接的杆件构成的空间结构。1§6-1空间桁架结构有限元分析（桁架单元）

由于两端铰接，不能承受力矩，假定每个节点三个自由度，在整体参考系有：为方便研究，在单元上取一个局部参考系，以为原点，由点指向点为轴的正向。在局部参考系中有：一、桁架单元刚度矩阵2在单元中假定：利用节点处的位移连续性条件可以解出和，代回原式有：其中单元长度为：3其中为单元截面积，为单元材料杨氏模量。

需要指出，以上单元刚度是基于局部参考系得到的，由于不同的单元，其局部参考系的方位各不相同。因而，局部参考系中的单元节点位移向量的方向也各不相同，局部参考系中的刚度矩阵不具有可累加性。

这就要求我们将局部参考系下的单元刚度矩阵变换到整体参考系中。4二、坐标变换矩阵局部参考系与整体参考系中单元节点位移向量的变换关系为：写成矩阵形式为：简记为：其中：5单元弹性应变能在整体参考系中的表达式为：单元弹性应变能在局部参考系中的表达式为：将：代入上式有：将该式与整体参考系中的单元弹性应变能对比可知：6三、等效单元体力载荷向量桁架单元的等效体力载荷向量可以简单地表示为：其中：为单元形心坐标，可以由下式给出：7四、单元应力计算8§6-2平面刚架（框架）结构有限元分析一、平面刚架单元受力分析与变形分解该单元可以承受和方向的力以及绕轴方向的力矩。在整体参考系中单元节点位移向量为：为方便研究，在单元上取一个局部参考系，以为原点，由点指向点为轴的正象。9在局部参考系中单元所受载荷可以分解为：1）轴向拉压，主要引起单元的轴向变形。2）横向载荷及绕平面法向的弯矩，主要引起单元的弯曲。为研究方便，我们在做刚架结构有限元分析时，也将单元变形分解为轴向变形与弯曲变形两部分，分别进行研究。在小变形条件下，可以认为单元的轴向变形与弯曲变形是相互独立的并符合迭加原理。因此，我们可以分别计算两种基本变形的刚度，在进行迭加。刚架结构的有限元分析在研究弯曲变形时采用了材料力学梁的基本假定，因此，有时也将刚架单元称为梁单元。101）单元轴向变形刚度矩阵与轴向变形相关的节点位移分量有和，假定在轴向载荷作用下单元内单元内任意点的位移为：与轴向变形对应的刚度矩阵为：

二、平面刚架单元两种基本变形的的刚度矩阵112）单元弯曲变形刚度与弯曲变形相关的节点位移分量有、、和。假定在轴向载荷作用下单元内单元内任意点的挠度为：中性层上的根据节点处的位移连续性条件有：解出，代回挠度函数经整理可得：在点：在点：小变性、中心平面、12挠度函数为：其中：单元中性层外任意点沿局部参考系方向得位移可以表示为：13单元中性层外任意点沿局部参考系方向得应变可以表示为：其中：在局部参考系中单元刚度矩阵可以表示为其中：——绕轴截面弯曲惯性矩14将轴向变形与弯曲变形的刚度迭加，可以得到平面刚架单元的刚度矩阵。15§6-3空间刚架（框架）结构有限元分析一、空间刚架单元受力分析与变形分解该单元可以承受三个方向的力以及绕三个轴的力矩。在整体参考系中单元节点位移向量为：为方便研究，在单元上取一个局部参考系，以为原点，由点指向点为轴。以截面的两个惯性主轴分别为局部参考系的轴和轴。161)在轴向载荷的作用下，将产生轴向拉压变形，与之相关的节点位移分量为：2)在绕轴向载荷（扭矩）的作用下，将产生扭转变形，与之相关的节点位移分量为：3)在向载荷及绕轴力矩的作用下，将产生平面内的弯曲变形，与之相关的节点位移分量为：4)在向载荷及绕轴力矩的作用下，将产生平面内的弯曲变形，与之相关的节点位移分量为：17二、空间刚架单元四种基本变形的刚度矩阵1）单元轴向变形刚度矩阵与轴向变形相关的节点位移分量有和，假定在轴向载荷作用下单元内单元内任意点的位移为：与轴向变形对应的刚度矩阵为：

182）绕单元轴向扭转变形刚度矩阵与扭转变形相关的节点位移分量有和，假定在轴向载荷作用下单元内单元内任意点的扭转角为：在单元中假定：在点：在点：解出和代回插值函数有：

对于圆形截面有：其几何矩阵为：19与轴向变形对应的刚度矩阵为：

根据应力应变关系可以得到则扭转变形刚度矩阵为：令：——圆形截面的扭转惯性矩。说明：对于非圆截面前面的应变位移关系及扭转惯性矩算式不成立，但刚度矩阵是普遍适用的，只是扭转惯性矩需要用其它方法获得，如：近似公式计算、实验等。203）单元平面内弯曲变形刚度矩阵与弯曲变形相关的节点位移分量有、、和。假定在横向载荷和弯矩作用下单元内单元内任意点的挠度为：根据节点处的位移连续性条件有：解出，代回挠度函数经整理可得：在点：在点：21挠度函数为：其中：单元中性层外任意点沿局部参考系方向得位移可以表示为：22单元中性层外任意点沿局部参考系方向得应变可以表示为：其中：在局部参考系中单元刚度矩阵可以表示为其中：——绕轴截面惯性矩234）单元平面内弯曲变形刚度矩阵与弯曲变形相关的节点位移分量有、、和。假定在单元内单元内任意点的挠度为：根据节点处的位移连续性条件有：在点：在点：其中：——绕轴截面惯性矩对应的单元刚度矩阵为：其他过程与平面内的弯曲相同，不再24对称空间刚架单元在局部参考系中的单元刚度矩阵为：25三、空间刚架单元刚度矩阵的变换关系整体参考系1）空间刚架单元的节点位移向量局部参考系2）空间刚架单元的节点位移向量的变换关系——坐标变换矩阵，具体形式如下：26其中：——局部参考系轴关于整体参考系的方向余弦——局部参考系轴关于整体参考系的方向余弦——局部参考系轴关于整体参考系的方向余弦因为空间刚架单元的局部参考系有三个坐标轴，根据单元的两个节点只能确定其中一个坐标轴的方向，为确定其它两个坐标轴的方向，还需要因入一个辅助量，可以是一个节点或一个参考方向。27四、空间刚架单元局部参考系确定方法（参考点法）所谓参考点法，是在单元截面主轴平面内找一个辅助点（称参考点），用以帮助确定单元局部参考系。首先由点到点作一个向量

其中：取的方向为局部参考系轴的方向，轴与整体参考系的三个方向余弦为：28再由点到点作一个向量

其中：显然的方向垂直于和，以局部参考系轴的方向。由向量差乘规则可知：令：或写成行列式形式：29这样，局部参考系轴的三个方向余弦

可以表示为：显然的方向垂直于和，以局部参考系轴的方向。