fem-06_平面问题有限元(三角形).pdfVIP

有限元方法与应用 平面问题有限元 有限单元法求解的基本步骤 • 结构的离散化处理 把任意形状的结构或连续体划分成有限个基本单元的组合，并在单元体的指 定点设置节点，把相邻的单元体在节点处连接起来，组成单元的集合体，代 替原有的结构或连续体。单元上所有载荷向节点简化，使弹性体离散后，假 定力可以通过节点由一个单元传递到另外一个单元。 • 单元的特性分析 在分析连续体问题时，必须对单元中的位移分布进行假设，假定位移是坐标 的某种简单函数，称之为位移函数。位移法有限元采用节点位移作为未知量， 在每一个单元内部用该单元的节点位移插值多项式表示单元内的近似位移函 数，{f}=[N]{d}e 根据所选用的单元位移函数，对单元进行力学特性分析： e e e e e {e}=[B]{d} ；{s}=[D][B]{d} ；{F} =[k] {d} • 整体分析求得总体刚度矩阵 包含两个方面的内容，一是将各单元的刚度矩阵集合成整体刚度矩阵，二是 将作用于各单元节点上的等效节点力列阵集合成总的载荷列向量。 [K]{d}= {F} 有限单元法求解的基本步骤(cont.) • 约束处理 总体刚度矩阵此时还是奇异的，不能直接用来求解，需要引入边界条件，修 正后的总刚仍还是对称的，且成为非奇异矩阵。 • 求解修正后的总刚方程，得到节点位移 求解方法一般是根据总刚的特性，选取高斯消去法（高斯循序消去法，三角 分解法，高斯—约当消去法）、迭代法（雅可比迭代法，高斯—赛德尔迭代 法，超松弛迭代法，共轭梯度法）等。 讨论： 直接解法的优点在于对于给定的方程组，一种选定了的直接解法可以在规定 步骤之内完成计算，操作数可以预先计算，方法简单。不足之处在于需要保 存系数矩阵中夹杂在非零元素之中的零元素，因为计算中要使用，因此，增 加了对计算机存储的要求，影响计算效率。同时，不能对解的误差进行检查 和控制。 迭代法的优点在于不保存系数矩阵中夹杂在非零元素之中的零元素，并且不 对它们进行运算，另外，迭代过程中可以检查解的误差，并可以通过增加迭 代次数来降低误差，直至满足要求，不足在于收敛性。 • 求解应变和应力。 有限元法的基本思想 将单元内任一节点 对弹性区域离散化 位移通过函数表达 （位移函数） 进行单元集成， 在节点上加外载荷力 建立单元方程 引入位移边界条件 求解得到节点位移 进行求解 根据弹性力学公式得到单元应变、应力 4 结构的离散化 三节点三角形元较为常见。 边界单元：至少带有两个边界节点的单元； 内部单元：只有一个或者没有边界节点的单元。 结构的离散化(cont.) 1. 划分单元时，必须沿整个求解区域进行。单元的大小和节点的疏密需要根据 求解精度和计算机容量等进行综合考虑。 2. 节点的安排需要考虑物体所承受的载荷和支撑情况。一般集中力作用点应处 理为节点，物体的铰支点应为节点。面力的起始位置和数值突变的部位应设 置为节点。另外，试验测试部位需布置节点。 3. 三角形单元划分时需注意三个边长不能相差悬殊，根据误差分析，应力和位 移的误差都和单元的最小内角的正弦值成反比，采用等边三角形或者直角等 腰三角形误差较小。 4. 采用三角形单元，应使一个单元的顶点必须同时又是相邻三角形单元的顶点， 而不是处于某一边上的内点。 5. 尽可能使每个内节点同时为六个三角形单元的顶点，且六个角度应相差不大。 6. 厚度有突变或者材料性质