《有限元理论与应用教案》辽宁工程技术大学YXY6杆件系统的有限单元法.pptVIP

* 第六章 ? 杆件系统的有限单元法 返回 第一节 等截面梁单元的刚度矩阵 第二节 等效节点力 第三节 梁单元刚度矩阵的坐标变换 在机械工程领域中，结构件的应用非常广泛。例如，起重运输机械的动臂、汽车车架、输变电塔架、油田抽油机和港口输油臂、大型广告牌等，都是由金属杆件所组成。对于这类复杂的杆系结构，由于结构和受力情况复杂，超静定次数很高，难以作精确分析，通常须采用有限元法来进行分析。应该注意的是，尽管杆件系统的有限单元法与结构力学中的经典位移法并没有很大的区别，但在有限单元法中，其“基本结构”的选取有所不同。在有限单元法中，凡是杆系的交叉点、边界点、集中力作用点都应取为节点，而节点之间的杆件均可作为单元。即用单元取代了经典位 移法中的“基本结构”。 返回 第一节 等截面梁单元的刚度矩阵 在同一平面内的若干杆件以焊接或铆接等方式连接起来的结构，若其所承受的载荷也在该平面内，称此结构为平面杆件系统，如图6-1所示。 图 6-2 平面梁单元 P P P P P 1 2 3 4 5 6 8 7 9 10 11 12 13 14 i i j y x j ui uj Mj Mi Qi Qi Qj Qj Nj Ni Vi Vj 6.1 平面杆件系统 图 6-1 平面杆件系统 返回 在平面杆件系统中，是取节点为i和j之间的梁为梁单元，如图6-2所示。采用右手坐标系，使x轴与梁轴重合，而y轴和z轴为梁截面的主惯性轴方向。因外载荷都在同一平面内，所以梁单元总是处于轴向拉压和平面弯曲的组合变形状态。在节点i和j上所受到的节点力为轴力、剪力和弯矩，即Ni、Qi、Mi和Nj、Qj、Mj；与之相对应的节点位移分别为ui、vi、θi和uj、vj、θj。 6.1.1单元的位移模式 对于轴向位移u，其位移模式可取x的线性函数，而挠度v则可用x的三次多项式来表示，即 （6-1） 返回 记 (6-2) 则 (6-3) 返回 由节点位移可以求得位移模式中的全部参数{a}和{b}。于是，梁单元的位移模式便可用其节点的节点位移来表示，其矩阵形式为 其中 若将节点i和j的位移记为 (6-4) 和 返回 及 则 (6-5) 式中 返回 6.1.2单元的应变和应力 由于梁单元受到拉压和弯曲变形后，其线应变可分为两部分：拉压应变εo弯曲应变εb。一般情况下，剪切应变对梁挠度的影响是微小的，可以忽略不计，于是 式中 (6-6) 由虎克定律，得 (6-7) 返回 6.1.3 单元的刚度矩阵 设单元内各点的虚位移为{f *}，则 , 根据弹性理论，梁单元内的应力由于虚应变所作的虚功为 而单元所受外力（集中力{F}e和沿轴线作用的分布载荷{q}）在虚位移上所作的功为 由虚位移原理，得 返回 令 则 式中 就是由分布载荷移置到节点的等效节点力。 返回 经积分运算，单元刚度矩阵的显式为 （6-8） 式中 是梁截面对主轴的惯性矩，A为梁截面面积。 返回 应该注意，当梁截面的高度大于梁长度的1/5时，剪切应变对挠度的影响就必须予以考虑，尤其是在薄壁截面的情形，剪切对挠度的影响将是巨大的。考虑剪切影响时，只需对梁单元的刚度矩阵作如下修正 式中 为剪切影响系数，As为有效抗剪面积。 （6-9） 返回 图6-3 空间梁单元 i j Myi(Qyi) Qyi(vi) y Mxi(Qxi) Ni(ui) Qxi (Wi) Mxi(Qxi) z Qyj(vj) Nj(uj) Mxi(Qxi) Qxj (Wj) Mxj(Qxj) x 6. 2 空间杆件系统 若杆件系统、截面主 轴或作用载荷不在同 一平面内，则这类情 况属于空间杆件系统 问题。一般情况下， 空间梁单元的每个节 点的位移具有六个自 由度，对应于六个节 点力。如图6-3所示。 返回 与平面杆件系统类似，若记 且 且 式中Ni和Nj表示作用于节点i和j的轴向力；Qyi、Qzi、Qyj、Qzj表示y和z方向的剪力；Mxi、Mxj表示扭矩；Myi、Mzi、Myj、Mzj表示绕y和z轴的弯矩。 同样地，用{Q}e表示由分布载荷移置的等效节点力，则: 返回 (6-10) 式中Iy、Iz是对y和z轴的主惯性矩； 、 是对y和 z轴方向的剪切影响系数；Jk是对x轴的扭转惯性矩；Ay、Az是 梁截面沿y和z轴方向的有效抗剪面积。 返回 第二节