第5章New(beam)han教程分析.docVIP

第五章 梁单元 本章介绍利用梁单元进行结构静力学的有限元分析原理。首先介绍了杆单元的分析方法，以杆单元为例给出了有限元的完整思路。然后导出了平面梁单元刚度矩阵，利用最小势能原理推导系统刚度矩阵。最后给出了空间梁单元模式。 5.1杆单元 20世纪50年代，Boeing公司的Turner, Clough，Martin, Topp在分析飞机结构时首先研究了离散杆、梁单元刚度的表达式，这种将复杂结构进行离散化的做法就是有限元法的开端。可以认为，有限元法的基本实质是最小势能原理的应用，其核心在于使用离散的方式组合表达全几何场上的形函数而不是直接寻找全场上的形函数。 有限元的总体思路是用标准化分段的单元逼近原结构（离散化）; 寻找能满足位移边界条件的许可位移场; 用基于位移场的最小势能原理求解位移场。采用的基本变量是节点位移，进而求出位移场、应变场、应力场等。 对于弹性变形体的三类变量（位移、应变、应力）和三大方程，采用试算函数的求解方法可以大大降低难度，且有较好的规范性和可操作性。主要包括基于试函数的加权残值法和基于试函数的最小势能法。在有限元中，采用后者作为理论基础。 5.1 结构力学中的杆件系统分析方法-位移法 下面以结构力学中的杆件系统求解为例，描述有限元的基本过程。对于一个杆系结构，通常是由许多结构杆件（单元）组成，这些单元仅在有限个节点上彼此相连。对每个单元而言，力与位移之间的关系等结构特性可以用节点上所确认的自由度来唯一地加以规定，而整体杆系的结构特性则可以通过组集这些单元特性加以描述。 图5-1所示的是一个由二根杆件组成的铰接桁架。杆件的截面积为A，弹性模量为E，长度分别为l1和l2 。该桁架在各铰接点处有外力X1、Y1、X2、Y2、X3、Y3 。因杆件在节点处是铰接，不承受弯矩，只能承受轴向力，所以，每个节点的力和位移各有两个分量，即每个节点均具有两个自由度，因此每个单元则有四个自由度。为此，用四个方程来描述每个单元的力与位移之间的关系。对于单元①，有 （5.1） (a) (b) (c) 图5-1铰接桁架 式中 U11、V11、U21、V21 分别为节点1和节点2施于单元①的节点力沿坐标方向的分量；u11、v11、u21、v21 分别为节点1和节点2的位移沿相应坐标方向的分量。上标是单元编号，下标为节点编号。 可将式（5-1）写成矩阵形式 （5.2） 或简记为 （5.3） 式中被称为单元①的节点位移向量；称为单元①的节点力向量；[k]1 为单元①的刚度矩阵，其中k11，k12，…，k44 为刚度系数。 若令 u11 =1，v11 = u21 = v21 =0 则得 U11 = k11 ，V11 = k21 ，U21 = k31 ，V21 = k41 这表明，当节点1沿x方向产生一单位位移（u11 =1），而其余节点、其余自由度方向上的位移为零时，刚度系数就等于施加于该单元各自由度方向上的力。 图5-2 杆单元(虚实线对调，加力) 如图5-2所示，在此状态下，单元的长度将缩短 (l1 = cos( 。根据材料力学知识，在节点1处需要向单元①施加的轴向压力F为 由此可以求得作用于单元①在节点1处两个自由度方向上的力，即该轴向压力在x和y 方向上的分量 而作用于节点2的两个自由度方向上的力，则与之大小相等、方向相反。即 再继续作类似的分析，便可得到其余的刚度系数，有 (5.4) 同理，可以求出单元②的刚度矩阵，这时式(5.4)中有。 (5.5) 为了获得整体结构的力与位移之间的关系，需要引入整体结构的节点位移分量u1、v1、u2、v2、u3、v3 和单元节点位移分量、、、 及、 （上标为单元编号）之间的协调关系，即 u2 = u11 v2 = v11 u2 = u21 = u22 v2 = v21 = v22 u3 = u32 v3 = v32 另外，根据力的平衡条件，作用在节点上的外力应该等于与该节点相连的各单元所受到的节点力之和。因此可得到结构的力与位移的关系 (5.6) 写成矩阵形式，有 (5.7) 或简记为 (5.8) 这就是有限单元法所要建立的基本方程组。上式中是由作用在节点上的外载荷所组成的向量，称为载荷向量（或载荷列阵）；是由基本未知