1三角形单元有限元解读.ppt

有限元法 目前在工程技术领域内常用的数值模拟方法有：有限单元法（Finite Element Method，FEM）、边界元法（Boundary Element Method，BEM）和有限差分法（Finite Difference Method，FDM）等，但就其实用性和应用的广泛性而言，主要还是有限单元法。作为一种离散化的数值解法，有限单元法首先在结构分析，然后又在其他领域中得到广泛应用。 发展与现状 离散化的思想可以追溯到20世纪40年代。1941年A.Hrennikoff首次提出用离散元素法求解弹性力学问题，当时仅限于用杆系结构来构造离散模型，但能很好地说明有限元的思想。如果原结构是杆系，这种方法的解是精确的，发展到现在就是大家熟知的矩阵分析法。究其实质这还不能说就是有限单元法的思想，但结合以后的有限元理论，统称为广义有限单元法。1943年R．Courant在求解扭转问题时为了表征翘曲函数而将截面分成若干三角形区域，在各三角形区域设定一个线性的翘曲函数，这实质上就是有限单元法的基本思想（对里兹法的推广），这一思想真正用于工程中是在电子计算机出现后。 20世纪50年代因航空工业的需要，美国波音公司的专家首次采用三节点三角形单元，将矩阵位移法用到平面问题上。同时，联邦德国斯图加特大学的J. H. Argyris教授发表了一组能量原理与矩阵分析的论文，为这一方法的理论基础作出了杰出贡献。1960年美国的R. W. Clough教授在一篇题为“平面应力分析的有限单元法”的论文中首先使用“有限单元法（the Finite Element Method）”一词，此后这一名称得到广泛承认。 20世纪70年代以来，有限单元法进一步得到蓬勃发展，其应用范围扩展到所有工程领域，成为连续介质问题数值解法中最活跃的分支。由变分法有限元扩展到加权残数法与能量平衡法有限元，由弹性力学平面问题扩展到空间问题、板壳问题，由静力平衡问题扩展到稳定性问题、动力问题和波动问题，由线性问题扩展到非线性问题，分析的对象从弹性材料扩展到塑性、粘弹性、粘塑性和复合材料等，由结构分析扩展到结构优化乃至于设计自动化，从固体力学扩展到流体力学、传热学、电磁学等领域。它使许多复杂的工程分析问题迎刃而解。 有限单元法的基本思想是将物体（即连续的求解域）离散成有限个且按一定方式相互联结在一起的单元的组合，来模拟或逼近原来的物体，从而将一个连续的无限自由度问题简化为离散的有限自由度问题求解的一种数值分析法。物体被离散后，通过对其中各个单元进行单元分析，最终得到对整个物体的分析。网格划分中每一个小的块体称为单元。确定单元形状、单元之间相互联结的点称为节点。单元上节点处的结构内力为节点力，外力（有集中力、分布力等）为节点荷载。 (1-61） 由式(1-55）即可解出全部未知的节点位移： 1.7.7 节点位移和单元力的解答 1、结构节点位移 2、支座反力 把解出的??f?代入(1-56），即得支座反力?Pr?： 关于方程 (1-61）的解算方法，当[Kff]采用本章中上述方法组集时，可直接采用结构力学中的高斯（Gauss）法求解。 (1-56） 至此，我们可以看出：系统分析的主要任务是： （1）组集引入约束后的结构刚度矩阵[Kff]； （2）求解式(1-55）给出的线性代数方程组。算出全部未知的节点位移。 至于支座反力的计算，实际计算时，根本不去组集式(1-56）中的矩阵[Krf]，不用式(1-62）而是直接对支承点使用节点平衡方程计算。 3、单元力 (1-62） （1）全解公式 全解是指结构在未经简化的实际荷载作用下（全解系统）的解答。 余解是结构在节点荷载作用下（余解系统）的解答。节点荷载包括原来给定的节点荷载和由分布荷载简化而来的等价节点荷载。 特解是单元从实际分布荷载到等价节点力过程中（特解系统）导致的单元节点力。 全解 = 余解 + 特解 (1-63） 通过式(1-55）解出获得全部未知节点位移???，据此找出任何一个单元的单元位移（单元节点位移），然后可以求得单元力。 图(1-15） i j Mi Mj Ri Rj = + (e) i j （a）全解系统 （b）特解系统 （c）余解系统 （2）余解 -Mi -Mj -Ri -Rj i j （3）特解 把单元等价节点力反向作用在单元的相应节点，即为单元力的特解。 4、结论 （1）全解的节点位移等于余解的节点位移 因为，特解给出的节点位移总是为零。 （2）对不承受分布荷载的单元，内力由余