弹性力学及有限元法 chapter9.ppt

第九章 杆件系统的有限元法 要 点 了解杆件系统有限元法的解是精确解； 杆件系统中的单元都在单元的局部座标中讨 论，然后再投影到整体座标中去； 杆件系统中的单元是梁单元，其变形以弯曲 和拉伸变形为主，常忽略剪切变形的影响。 第九章 杆件系统的有限元法 所研究的对象是细长的杆件，即轴线方向的尺寸远比其他二个方向的尺寸大的多； 杆件系统可分为平面杆件系统和空间杆件系统，若组成结构物的杆件都在同一个平面内，外力也在这一平面内，则称为平面杆件系统； 若组成结构物的杆件不在同一个平面内，则称为 空间杆件系统； 杆件系统中所用的单元主要为梁单元，如果各杆件只受拉压作用则采用桁架单元。 第九章 杆件系统的有限元法 在梁单元上建立局部座标 oxyz，以梁单元的轴线ij作为x轴， y轴垂直x轴。以右手法则确定座标轴系。对于等截面梁单元有二个节点i, j，每个 节点有三个自由度 ui，vi， θi，它的 节点位移{δ}e 第九章 杆件系统的有限元法 它的节点力{R}e 第九章 杆件系统的有限元法 把上述位移函数写成矩阵形式 令[1 x ]为[h(x)], 令[1 x x2 x3]为[H(x)], 为{a} 为{b} 现根据节点位移解出{a}和{b} 第九章 杆件系统的有限元法 当 x=0 时 ui=a0 当 x=l 时 uj=a0+a1l 当 x=0 时 vi=b0 当 x=l 时 vj=b0+b1l+b2l2+b3l3 当 x=0 时 θi=dv/dx=b1 当 x=l 时 θj=dv/dx=b1+2b2l+3b3l2 写成矩阵形式 第九章 杆件系统的有限元法 其中 第九章 杆件系统的有限元法 又可写成下式: 又令 第九章 杆件系统的有限元法 把[h(x)], [H(x)]扩大，并把[A1]-1, [A2]-1扩大 第九章 杆件系统的有限元法 合并后[A] 如下式: 第九章 杆件系统的有限元法 梁单元的应变{ε} 和应力{σ} 梁单元的应变{ε} 第九章 杆件系统的有限元法 梁单元的的刚度矩阵[K]e 根据虚功原理: 第九章 杆件系统的有限元法 展开后得: 第九章 杆件系统的有限元法 相乘以后 dV=dAdx 其中 为梁截面的主惯性矩 第九章 杆件系统的有限元法 §9-3 等截面空间梁单元的刚度矩阵 第九章 杆件系统的有限元法 第九章 杆件系统的有限元法 分布轴向力p(x)的移置: 由于是轴向拉伸，形态矩阵用[Nu]。 第九章 杆件系统的有限元法 第九章 杆件系统的有限元法 分布扭转力矩Mx(x)的移置: 第九章 杆件系统的有限元法 分布横向力q(x)的移置: 第九章 杆件系统的有限元法 若是均布力 q(x)=q 则 第九章 杆件系统的有限元法 可得等效节点力如下: 第九章 杆件系统的有限元法 §9-5 座标变换 第九章 杆件系统的有限元法 {R’}e和{R}e之间， {δ’}e和 {δ}e之间的转换关系 第九章 杆件系统的有限元法 根据投影规则 第九章 杆件系统的有限元法 写成矩阵形式 第九章 杆件系统的有限元法 [K’]e和 [K]e之间的转换关系 * * §9-1 杆件系统的有限元法的研究对象 §9-2 等截面平面梁单元的刚度矩阵 梁单元的位移模式 我们选择梁单元的位移模式 在梁单元中，规定位移u,v的正方向与座标轴方向一致，转角θ的正向为逆时针转向；规定力N,Q的正方向与座标轴方向一致，力矩M的正方向为逆时针转向。 由此解出{a},{b} 把{a},{b}求出后就可获得梁单元的位移模式，代入上式 可得: 由于在位移函数u中的节点位移为ui，uj，位移函数v中的节点位移为 vi，θi， vj， θj ，而现在要把它改写成{δ}6*1e，所以需要改写。 把[A1]-1和[A2]-1合并成[A] 则位移模式如下: 于是{ε}=[B]{δ}e。其中 梁单元的应力{σ} 则梁单元的刚度矩阵[K]e如下:当l/h 5时忽略剪切的影响 空间梁单元的节点位移{δ}e 空间梁单元的节点力{R}e 同理可得空间梁单元的刚度矩阵 其中