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

第五章杆系结构的有限元法 5.1引言 杆系结构是工程中应用较为广泛的结构体系，包括平面或空间形式的梁、桁架、刚架、拱等。其 组成形式虽然复杂多样，但用计算机进行分析时却较为简单。杆系结构中的每个杆件都是一个明显的 单元。杆件的两个端点自然形成有限元法的节点，杆件与杆件之间则用节点相连接。显然，只要建立 起杆件两端位移与杆端力之间的关系，则整体平衡方程的建立与前几章完全相同。 杆端位移与杆端力之间的关系，可用多种方法建立，包括前面几章一直采用的虚功原理，但是采 用材料力学、结构力学的某些结论，不仅物理概念清晰、直观，而且推导过程简单明了。因此，本章 将采用这种方法进行单元分析。至于整体平衡方程的建立，则和前面几章所讲的方法一样，即借助于 单位定位向量，利用单元集成法进行。 5.2平面桁架的有限元分析 平面桁架在计算上有以下几个特点： 杆件的每个节点仅有两个线位移； 杆件之间的连接为理想铰，即在节点处各杆件可相对自由转动，且杆件轴线交于一点。 外载荷均为作用于节点的集中力。 由于以上特点，所以在理论上各杆件只产生轴向拉、压力，截面应力分布均匀，材料可得到充分 利用，因此桁架结构往往用于大跨结构。 5.2.1局部坐标系下的单元刚度矩阵 从平面桁架中任取一根杆件作为单元，称作桁架单元，单元长为L，横截面面积为A，图5.1。 两端节点分别用i和j表示，规定从i到j的连线方向为局部坐标x轴，垂直于三的方向为亍轴。 由于桁架中各杆只产生轴向力和轴向变形，所以节点i 和 j只发生沿小向的位移，用%和七表 示，相应的杆端轴力分别用匚和Fj表示。由虎克定律可推得 F --^ (u - u ) (u - u )F ^ (u - u ) xj 将这两个式子写成矩阵形式，就是 r \F r \ F e _ EA EA 1 r—? u xi = L L U i EA EA 1 F 1 xj」 _ L L \ u j e (5.1) 显然，在局部坐标系下，i、j两节点沿亍轴方向的位移亏=二=0，在亍轴方向的节点力 F = T 0。因此，可以把61)r——) FeEA0EA — F = T 0。因此，可以把61) r——) F e EA 0 EA — L L F 0 0 0 4^ ]= EA EA F 0 F L L yj 0 0 0 扩大为下面的四阶的形式 0 0 0 0 L〕 u. V ―^ u 广 I jJ 可以简写为 =伙]e {8 }e 其中 称作桁架单元的单元杆端力向量。 称作桁架单元的杆端位移向量 而 F }=卜 xi v i EA EA EA L 0 EA 00 0 0 0 0 (5.2) (5.3) (5.4) (5.5) (5.6) 称作桁架的单元刚度矩阵，式(5.2)或式(5.3)就是桁架的单元刚度方程，它反映了单元杆端力与杆 端位移之间的关系。 5.2.2整体坐标系下的单元刚度矩阵 在一个复杂的结构中，各个杆件的杆轴方向不尽相同，因而各自的局部坐标系也不尽相同。为了 建立结构的整体平衡方程，必须选用一个统一的公共坐标系，称为整体坐标系，用x,y表示。 首先分析单元杆端力在不同坐标系中的关系。图5.2所示任一单元e，其局部坐标系为0X,，整 体坐标系为oxy，由x轴到X轴的夹角a以顺时针转向为正。局部坐标系中的杆端力用厂、匚。表 示。整体坐标系中的杆端力则用F：、Fy表示，如图5.2所示，显然。二者有下列关系。 图5.2F = F . cos以 + F sin 以F = -F sin 以 + F cos 图5.2 F = F . cos以 + F sin 以 F = -F sin 以 + F cos以 F = F cos以 + F sin 以 F = -F sin 以 + F cos 以 (5.7) 将式(5.7)写成矩阵: (——1 F 妇 Fy _y_ = F xj F 1 yj cos a sin a 0 0 -sin a cosa 0 0 0 0 cosa sin a 0 0 -sin a cosa \ FT F y. F F l yj) (5.8) 式简写为 {F }e = [T ]{F }e (5.9) 式中卜］称为单元坐标转换矩阵 -sin a cosa 0 0 0 0 cosa sin a 0 0 -sin a cosa cos a sin a 0 0 [T ]= (5.10) [T ]-i = [T ]t (5.11) 或 [T ][T ]t = [T ]-t [T ] = [ I ] (5.12) 式中］为与T］同阶的单位阵。 结合式(5.12)，由式(5.9)得 {F }e = [T ]{F }e (5.13) 同理，可以求出单元杆端位移在两种坐标系中的转换关系。设局部坐标中单元杆端位移向量为 肉