西南交大《有限元方法》课件.ppt

Ⅲ、面积坐标的微积分运算 由(2.50)， 一般复合函数求导法则： (2.52) 面积坐标在三角形全面积上的积分公式， (2.53) 在任一条边上如的积分公式 ， (2.54) 三角形单元： Ⅳ、单元形函数 由面积坐标与形函数比较： 直接用面积坐标表示形（插值）函数，而不必先求解待定系数与节点位移之关系，不必求解节点坐标构成的逆矩阵。 对三节点单元， (2.55) 线性单元的三个形函数就是三角形单元的三个面积坐标 对六节点单元： 节点1、2、3坐标， 由单元任意点p移动到1-2边中点，既构成节点4。 在1-2边上， 节点4坐标， 同理有节点5、6坐标， 各直线段方程： 1-4-2： 2-5-3： 3-6-1： 4—6： 4—5： 5—6： 按节点分块形式建立形函数， 对于每个节点，可选择几条直线，使得直线通过除节点本身以外其余全部节点；再利用直线方程的左部为函数因式来构造形函数。 如节点1：可选二条直线，4—6、2—5—3，即包揽了剩余全部节点。 分别取二直线方程的左部为形函数的因式，得到节点1的形函数， (2.56a) 显然，对于2~6节点，代入 中，均有 ，仅在节点1， 。 对节点4，可选二条直线，3—6—1、2—5—3， ，由节点4的坐标和形函数值定待定系数C (2.56b) 同理， (2.56c) 校核： 通过形函数即可建立单元位移模式： (2.57) 2.7计算流程 以平面问题为例，详细介绍有限元求解的步骤。 Step1： 结构离散化 A、选择正确的物理模型 B、选择合理的单元、节点、位移模式 C、选择高精度单元 D、同一个结构，单元类型尽可能少 E、载荷突变点应选为节点 F、单元划分数量 G、节点编号 Step2： 单元特性分析计算 针对不同单元建立其刚度矩阵，并按节点分块； 同时计算单元上外载荷的等效节点载荷 Step3： 形成总刚，总载荷列阵 按总节点编号，把分块后的单元刚度子矩阵叠加存入相应的总刚中对应区域； 由总刚的对称性，一般只需存入下三角元素即可； 将单元的等效节点载荷按总结点编号叠加形成载荷列阵。 Step4： 边界条件 引入边界约束条件，消除物体的刚体位移，消除总体刚度矩阵的奇异性。 位移条件有两类：一是给定边界位移的约束条件，排除了刚体运动的可能性；另一类是给定边界一个有限位移值，类似于位移加载。 引入已知边界条件，未知量减少，有限元方程数目减少，重新组合将导致数据交换复杂化。有二种处理方式，可以维持方程维数、排序不变：局部修正法和大数法 A、局部修正法 将节点位移的指定值置入方程组中，替换相应位移变量所在序号的方程，仍保持方程组维数不变，对相关的刚度矩阵和载荷列阵作修正。 如节点 在y方向有给定位移 ： 令 中的元素 ，且第 行和第 列的其余元素均为零值； 令 中的第 个元素的值等于给定位移值 ， 中的第 j 元素都需减去该节点位移值 和原来 中第 行相应的第j列元素的乘积； 例： 设 ，引入上式 这样 ， 而 均为原方程的解答。 这种方式保持了总刚维数不变化，求解方程组数目不减少，避免有限元方程组的重新组合、节点重新排序等烦恼。 B、大数法 如节点 在y方向有给定位移 ： 将 中的对角线元素 乘以一个大数如 ；同时将对应 中的第 项，替换为给定的位移与大数之积，即 同样地，以上例为例，第三个方程变为： 由于 远大于其它项，这些项的影响可忽略不计， Step5： 求解有限元方程 有限元方程实质上是一组联立的线性代数方程组，在静力平衡中，为一组以节点位移为基本未知量的系统平衡方程； 有限元求解的效率很大程度上取决于这组方程的解法； 目前多采用直接解法求解：如高斯消去法、三角分解法、以及波前等；迭代法有高斯－赛德尔超松弛法等。 Step6： 应力计算（后处理） 从有限元方程，可直接求解出所有单元节点的位移，拟合出结构的位移曲线、变形规律。 由求出的节点位移，从 ，求解各单元内部的应变，应力。 对于二，三维弹性体，我们引入的位移模式只考虑了其位移连续性，未涉及其导数的连续性，所以结果不能保证应变应力在边界上连续。 一般计算各单元的应变应力后