有限元基础理论与方法预案.ppt

* 有限元分析基础理论与方法 　　对于结构中的所有节点，则有 (21) (22) 　　式中　n－节点总数。 　　将式(22)记为 式中，{q}=(q1， q2 ， … ，qn)T，是由所有节点的位移分量组成的列阵； {R}=(R1，R2，…，Rn)T，是由所有作用在节点上的载荷组成的列阵；[K]是要求的总刚度矩阵，表达式为 (23) * 有限元分析基础理论与方法 　　式(23a)为整个结构的平衡方程，称为有限元方程或刚度方程。 (23a) 　　3. 载荷转移 　　在式(23)中载荷列阵{R}的元素为作用在节点上的外载 荷，是集中力。但结构上的载荷除了集中力之外，往往还有分布的面力和体力，不可能只作用在节点上，另外，即使是集中力也不一定刚好作用在节点上，因此，需要将这些载荷转换为节点载荷,即将载荷向节点移置。 　　载荷移置遵循能量等效原则，即原载荷与移置载荷产生的节点载荷在虚位移上所做的功相等。对于给定位移的函数， * 有限元分析基础理论与方法 这种移置的结果是唯一的。 在线形位移函数情况下，也可以按静力等效原则进行移置。 　　载荷移置是在结构的局部区域内进行的。根据圣维南原理,这种移置可能在局部产生误差,但不会影响整个结构的力学特性。 　　(1) 集中力的移置 　　集中力的移置是面力和体力的基础。如图2所示,设平面单元e中某一点(x,y)作用集中力Pc 　　设Pc移置后产生的等效节点载荷为 　　 * 有限元分析基础理论与方法 如果节点发生虚位移?qe,则单元内任一点位移为 集中力Pc所做的虚功为 * 有限元分析基础理论与方法 等效节点载荷所做的虚功为 　　根据能量等效原则,有 由于虚位移是任意的，可以从上式两边同时约去，故有 或 　　由上两式可见，载荷移置的结果仅与单元形函数有关， (24) (24a) * 有限元分析基础理论与方法 当形函数确定后，移置的结果是唯一的。 　　(2) 面力的移置 　　设厚度为t的平面单元单位面积上作用的面力为Ps ＝(psx, psy)T,如图3所示。并将微元面积dA=tdl上的面力PsdA视为集中力，利用式(24)并积分，可得与面力等效的移置节点载荷为 或 　 (25) (25a) * 有限元分析基础理论与方法 根据形函数的特点，在ij边上有Nk=0，所以 。因此在ij边上作用的面力只能移置到该边的两个节点上。 　　(3) 体力的移置 　　设单元体积内作用的体力为Pv ＝(pvx, pvy)，若将微元体tdxdy上的体力Pv tdxdy视为集中力，则利用式(24)并积分，可得与体力等效的移置节点载荷为 或 　 (26) (26a) * 有限元分析基础理论与方法 　　一般情况下,由叠加原理,一个单元上总的节点载荷应为上述三种载荷移置结果之和,即 　　若是一个节点与多个单元相连,则节点载荷应为所有相关单元向该节点移置的载荷的叠加。因此,对于整个结构而言,有 式中，ne－结构的单元数量。 　 (28) (27) * 有限元分析基础理论与方法 　　4. 约束处理 　　在没有计入约束条件之前,结构在外力作用下,会产生刚体运动,总刚度矩阵Κ是奇异的,因而,无法求解平衡方程获得节点位移值。故必须对结构施加足够的几何约束,消除结构刚体运动,使得Κ为非奇异阵。 　　如果约束是位移边界条件为零,则将总刚度矩阵中的零位移分量所对应行和列的主对角元素置1,其他元素置为0。在节点载荷列阵中，将零位移分量所对应的节点载荷也变为0。 　　如果节点位移是大于零的已知值,则将该位移分量所对应的主对角元素置为大数,再将载荷列阵中对应的分量置为大数乘以已知的节点位移,其余各元素保持不变。 * 有限元分析基础理论与方法 　　5. 求解线性方程组 　　对矩阵K、R进行约束处理之后，式(23)变为 式中　　。　　经约束处理后的总刚度矩阵和载荷列阵。 　　式(29)是一个关于节点位移分量的线性方程组，利用适当方法，如高斯消元法、高斯－赛德尔迭代法或波前法等就可求出所有节点位移。 　　 　　6. 计算其他物理量 　　在解式(29)求出基本未知量－节点位移分量后，便可利用式(4)、式(7)式(8)求出单元上任一点的位移及单元应变和应力。 (29) * 有限元分析基础理论与方法 　　7. 计算结果处理 　　由于在相邻单元的公共边界上应力和应变存在不连续而发生突变，为使计算结果接近实际值,通常需要对单元应力、应变进行平均化处理。处理的方法有两种:一是绕节点平均;另一种是两单元之间平均。绕节点平均就是将环绕节点的所有单元的应力加以平均,并用该平均值表示该节点的应力。两单元平均就是将两个相邻单元的常量应力加以平均,并用该