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

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