有限单元法基本思想-原理-数值计算过程.docVIP

有限单元法学习报告 在对力学问题分析求解过程中，方法可以概括为两种方法，一种为解析法，对具体问题具体分析，通过一定的推导用具体的表达式获得解答，由于实际工程中结构物的复杂性，此方法在处理工程问题是十分困难的；另一种是数值法，有限元法是其中一种方法，其数学逻辑严谨，物理概念清晰，又采用矩阵形式表达基本公式，便于计算机编程，因此在工程问题中获得广泛的应用。 有限元法基本原理是，将复杂的连续体划分为简单的单元体；将无限自由度问题化为有限自由度问题，因为单元体个数是有限的；将偏微分方程求解问题化为有限个代数方程组的求解问题。通常以位移为基本未知量，通过虚功原理和最小势能原理来求解。 基本思想是先化整为零，即离散化整体结构，把整体结构看作是由若干个通过结点相连的单元体组成的整体；再积零为整，通过结点的平衡来建立代数方程组，最后计算出结果。 我将采用最简单的三结点三角形为基本单元体，解决弹性力学中的平面问题为例，解释有限单元法的基本原理、演示数值计算过程和一般性应用结论。 一、离散化 解决平面问题时，主要单元类型包括三角形单元（三结点、六结点）和四边形单元（四结点矩形、四结点四边形、八结点四边形）等。选用不同的单元会有不同的精度，划分的单元数越多，精度越高，但计算量也会越大。因此在边界曲折，应力集中处单元的尺寸要小些，但最大与最小单元的尺寸倍数不宜过大。 在集中力作用点及分布力突变的点宜选为结点，不同厚度，不同材料不能划分在同一单元中。三角形单元以内角接近60°为最好。充分利用对称性与反对称性。 单元分析 将一个单元上的所有未知量用结点位移表示，并将分布在单元上的外力等效到结点上。 1、位移函数选取： 根据有限元法的基本思路，将连续体离散为有限的单元集合后，此时单元体满足连续性、均匀性、各向同性、完全线弹性假设。单元与单元之间通过结点连接并传递力，位移法（应用最广）以结点位移δi=(ui vi)T为基本未知量，以离散位移场代替连续位移场。单元体内的位移变化可以用位移函数（位移模式）来表示，因为有限元分析所得结果是近似结果，为了保证计算精度和收敛性，x位移函数应尽可能反应物体中的真实位移，即满足完备性和连续性的要求： (位移模式必须能反映单元的刚体位移。 (位移模式必须能反映单元的常量应变。 (位移模式应尽可能反应位移的连续性。 设三角形单元三个结点编号为i、j、m。平面三角形单元位移函数选取为 u=α1+α2x+α3y v=α4+α5x+α6y 可以写成的形式，反映了单元的刚体平动，反映了单元的刚体转动，满足完备性和连续性的要求(。 采用插值法由单元结点位移列阵δe=T计算α1、α2、α3、α4、α5、α6.，求出位移d=[u（x，y), v（x，y）]。6个未知量,6个代数方程,得de=Nδe de==T 式中Ni=(ai+bix+ciy)/2A，ai=bi= -ci=（i、j、m轮换）A为三角形面积，为避免A0，i、j、m按逆时针排列。N为形函数矩阵，形函数Ni的性质有： (Ni（xi，yi）=1 Ni（xj，yj）=0 Ni（xm，ym）=0 (Ni（x，y）+Nj（x，y）+Nm（x，y）=1可推出三个形函数中，两个是独立的，反映了刚体平移。 令z=Ni，在直接坐标系中画出Ni、Nj、Nm的函数图形是以Ni（xi，yi）=1为高的四面体，所以结点位移影响单元的位移场，单元的位移场是线性分布的，相邻单元在公共边上的位移是连续的，单元相邻边的位移只取决于单元相邻公共边上的结点而与其他结点无关，无论以哪个单元计算相邻边的位移，结果一定相同。 形函数Nie决定了单元内的位移模式，反映了i结点位移对单元内任意点位移的贡献率。 2、根据几何方程用单元结点位移表示单元应变： T B为几何矩阵 B可写为分块矩阵B=（Bi Bj Bm）T,Bi=，B内所有元素与x，y无关，所以该单元内应变是常量，反映单元的常量应变，满足完备性和连续性的要求(，这是一种常应变单元。 3、根据物理方程用单元结点位移表示单元应力： D为弹性矩阵 S为应力矩阵 S=DB中，每一个元素都是常数，所以的每一个分量与单元内x，y位置无关，这是一种常应力单元。 因为在三结点三角形单元中，位移函数中含有坐标的一次项，其误差为，而应力、 应变是常量，其误差为，比位移精度低。 4、根据虚功原理用单元结点位移表示单元结点力 单元在结点处受力，单元会发生变形，因此单元在结点处所受到的力与单元结点位移肯定有关系。单元间通过结点的相互作用成为整体，因此每一单元的受力——