有限元第五章杆系结构有限元法.docVIP

杆系结构的有限元法 5.1 引言 杆系结构是工程中应用较为广泛的结构体系，包括平面或空间形式的梁、桁架、刚架、拱等。其组成形式虽然复杂多样，但用计算机进行分析时却较为简单。杆系结构中的每个杆件都是一个明显的单元。杆件的两个端点自然形成有限元法的节点，杆件与杆件之间则用节点相连接。显然，只要建立起杆件两端位移与杆端力之间的关系，则整体平衡方程的建立与前几章完全相同。 杆端位移与杆端力之间的关系，可用多种方法建立，包括前面几章一直采用的虚功原理，但是采用材料力学、结构力学的某些结论，不仅物理概念清晰、直观，而且推导过程简单明了。因此，本章将采用这种方法进行单元分析。至于整体平衡方程的建立，则和前面几章所讲的方法一样，即借助于单位定位向量，利用单元集成法进行。 5.2 平面桁架的有限元分析 平面桁架在计算上有以下几个特点： 杆件的每个节点仅有两个线位移； 杆件之间的连接为理想铰，即在节点处各杆件可相对自由转动，且杆件轴线交于一点。 外载荷均为作用于节点的集中力。 由于以上特点，所以在理论上各杆件只产生轴向拉、压力，截面应力分布均匀，材料可得到充分利用，因此桁架结构往往用于大跨结构。 5.2.1 局部坐标系下的单元刚度矩阵 从平面桁架中任取一根杆件作为单元，称作桁架单元，单元长为L，横截面面积为A，图5.1。两端节点分别用i和j表示，规定从i到j的连线方向为局部坐标轴，垂直于的方向为轴。 图5.1 由于桁架中各杆只产生轴向力和轴向变形，所以节点i和j只发生沿方向的位移，用和表示，相应的杆端轴力分别用和表示。由虎克定律可推得 将这两个式子写成矩阵形式，就是 (5.1) 显然，在局部坐标系下，i、j两节点沿轴方向的位移，在轴方向的节点力。因此，可以把（5.1）扩大为下面的四阶的形式 (5.2) 可以简写为 (5.3) 其中 (5.4) 称作桁架单元的单元杆端力向量。 (5.5) 称作桁架单元的杆端位移向量 而 (5.6) 称作桁架的单元刚度矩阵，式(5.2)或式(5.3)就是桁架的单元刚度方程，它反映了单元杆端力与杆端位移之间的关系。 5.2.2 整体坐标系下的单元刚度矩阵 在一个复杂的结构中，各个杆件的杆轴方向不尽相同，因而各自的局部坐标系也不尽相同。为了建立结构的整体平衡方程，必须选用一个统一的公共坐标系，称为整体坐标系，用x,y表示。 首先分析单元杆端力在不同坐标系中的关系。图5.2 所示任一单元e，其局部坐标系为0，整体坐标系为oxy，由x轴到轴的夹角α以顺时针转向为正。局部坐标系中的杆端力用、表示。整体坐标系中的杆端力则用、表示，如图5.2所示，显然。二者有下列关系。 图5.2 (5.7) 将式(5.7)写成矩阵： (5.8) 式简写为 (5.9) 式中称为单元坐标转换矩阵 (5.10) 容易证明，单元坐标转换矩阵 是一个正交矩阵。因此有 (5.11) 或 (5.12) 式中为与 同阶的单位阵。 结合式(5.12)，由式(5.9)得 (5.13) 同理，可以求出单元杆端位移在两种坐标系中的转换关系。设局部坐标中单元杆端位移向量为，整体坐标系中单元杆端位移向量为，则 (5.14) (5.15) 式中 现在来推导单元刚度矩阵在两种坐标系中的转换关系。 单元杆端力与杆端位移在整体坐标系中的关系式可写为 (5.16) 式中称为在整体坐标系中的单元刚度矩阵。 将式中(5.9)和(5.14)代入(5.3)，得 将此式两边各前乘，并利用式(5.12)得 再将上式与式(5.16)比较，可知 (5.17) 这就是单元刚度矩阵在两种坐标系中的转换关系。 5.2.3 整体平衡方程和单元杆端力的计算 整体平衡方程由单元集成法建立，引入约束条件后，求解该方程可