有限元基础研究.ppt

* 有限元分析基础理论与方法 等效节点载荷所做的虚功为 　　根据能量等效原则,有 由于虚位移是任意的，可以从上式两边同时约去，故有 或 　　由上两式可见，载荷移置的结果仅与单元形函数有关， (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. 计算结果处理 　　由于在相邻单元的公共边界上应力和应变存在不连续而发生突变，为使计算结果接近实际值,通常需要对单元应力、应变进行平均化处理。处理的方法有两种:一是绕节点平均;另一种是两单元之间平均。绕节点平均就是将环绕节点的所有单元的应力加以平均,并用该平均值表示该节点的应力。两单元平均就是将两个相邻单元的常量应力加以平均,并用该平均值表示公共边中点的应力。 　　8. 结果显示 　　将计算得到的各种物理量以一定方式显示出来,研究结果的合理性和可靠性,评估结构性能和设计方案的优劣,做出相应的改进措施。目前数据显示方式有等值线、变形图、矢量图、二维或三维曲线以及动画等。 二　有限元分析基础理论与方法 * 弹性力学的基本理论(补充) (4) 弹性力学问题的解法 位移法: 取位移分量u,v,w为基本未知量，先利用位移表示的平衡微分方程和边界条件求解位移，再利用几何方程求应变，利用物理方程求应力。 应力法: 取应力分量为基本未知量。 混合法：同时取部分位移分量和应力分量作为基本未知量。 求解上述方程可知，弹性力学问题在数学上是由偏微分方程及其边界条件描述的,微分方程的建立和求解比较复杂, 只有在物体形状和受力较简单的情况下才能获得精确解。为了避免直接建立和求解这些微分方程的困难,提出了弹性力学微分方程的等价表达形式及其建立方法,即能量原理。 * 弹性力学的基本理论(补充) 3. 弹性力学中的能量原理 (1) 应变能 弹性体在受到外载荷作用发生变形的过程中,将把克服内力所做的功作为应变能存储在弹性体内部,当外力去除后,应变能做功,使弹性体恢复原状。 (2) 虚位移原理 将理论力学中刚体的虚功原理之推广到弹性体上,由于外力会使弹性体变形,因此,如果假定不存在热能和动能的改变,根据能量守恒定律,当处于平衡状态的弹性体发生体系所允许的任意微小位移时,外力在虚位移上所做的功等于虚位移发生时引起的弹性体的应变能增量,称为