数值模拟与仿真基础——6 等参数单元.pptVIP

数值模拟与仿真技术基础 麻向军 华南理工大学 聚合物新型成型装备国家工程研究中心 聚合物加工工程教育部重点实验室 等参数单元 6.1 等参数单元 6.2 平面问题的等参数单元 6.3 数值积分 6.1 等参数单元 1、等参数单元的概念 2、单元矩阵变换 1、等参数单元的概念 (1) 提高复杂几何计算精度的方法 采用较小的单元或高精度单元 但是这些单元几何形状规则，难以适应复杂边界。不能期望用较少的形状规则单元来离散复杂几何形状。 构造本身形状任意单元存在的问题： 难以构造出满足连续性条件的单元插值函数。 单元分析中出现的积分难以确定积分限。 解决策略： 通过坐标变换，将一个坐标系下形状复杂的几何边界映射到另一个坐标系下，生成形状简单的几何边界。 将满足收敛条件的形状规则的高精度单元作为基本单元，定义于局部坐标系（取自然坐标系），通过坐标变换映射到总体坐标系（取笛卡儿坐标系）中生成几何边界任意的单元，作为实际单元。 只有满足坐标变换的相容性——实际单元和基本单元之间的点一一对应，实际单元同样满足收敛条件。 构造的单元具有双重特性： 实际单元——几何特性、受力情况、力学性能来自真实结构，充分反映了它的属性。 基本单元——形状规则，便于计算和分析。 (2) 等参数变换 定义：坐标变换和单元内的场函数用相同的节点参数及相同的插值函数。 特性：满足坐标变换的相容性。 等参数单元：采用等参数变换的单元。 2、单元矩阵变换 (1) 为什么要进行单元矩阵变换 有限元分析中，为建立求解方程，需要进行各个单元体内和面积上的积分。 (2) 导数之间的变换关系 用复合求导法则： (3) 体积微元、面积微元的变换关系 (3) 体积微元、面积微元的变换关系 (4) 二维情况的蜕变关系 (5) 面积坐标与笛卡儿坐标之间的变换关系 6.2 平面问题等参数单元 1、四节点斜四边形单元 2、坐标变换式及变换行列式 1、四节点斜四边形单元 (1)、矩形单元 在矩形单元分析中，曾通过坐标变换，将其变换为一边长为2的正方形： 矩形单元的位移函数为： 这种双线性插值，对于任何双线性函数都是精确成立的，把它用来描述几何位置的坐标变量也必然成立。即： (2)、4节点斜四边形单元 在 4 边形四个顶点布置节点1 , 2 , 3 , 4 ，节点坐标为(xi , yi ) ( i = 1 , 2 , 3 , 4) 。 如果形函数仍然取矩形单元的双线性函数，则在单元的边界上，一般不能满足相容性条件。 不平行于 x( 或 y ) 轴的任一边上 ( 如 43 边 ) 的直线方程为： y = ax + b（a?0） 代入到位移函数式中，有： u = Ax2 + Bx + C 位移为x的二次函数。这表明在这条边上 u 不再是线性变化的。因此，这条边上的位移就不能由两个节点值的插值函数唯一确定，从而在相邻两个单元的公共边上将不能保证 位移是连续的。即相容性条件不满足。 因此，4 节点斜四边形单元不能象矩形单元那样，直接采用原直角坐标 ( x , y ) 表示的双线性函数为形函数。 为此，引入一个局部坐标系，进行等参数变换，使斜四边形单元能够变换到相应的矩形单元上去。 我们通过总体坐标 ( x , y ) 与局部坐标 ( ? , ? ) 之间的变换（或称为几何映射），使在总体坐标下的斜四边形单元变换为在局部坐标下边长为 2、坐标原点位于单元中心的正方形单元。 ( x , y) 1 , 2 , 3 , 4 ? (? ,?) 1，2，3，4 各边长的等分点 ? 各边长的等分点 从而保证一一对应的关系，即确定相应的坐标变换。 这种变换是否存在？一一对应关系需要满足什么条件？ 假设变换存在：则在局部坐标 ( ? , ? ) 下，单元是一个正方形，其位移函数为： 上式表示的位移函数为局部坐标 (? ,?) 的表达式，在有限元分析时，需要计算位移 u , v 对总体坐标的偏导数（ 等）。因此，采用上述做法时就必须找出总体坐标 (x , y) 对局部坐标 (? , ?) 之间的变换式。 如果采用与 位移函数 相同的形式对总体坐标变量插值：即 这样，我们可以看到位移函数式（4-25