平面问题有限元分析 等参数单元.pptVIP

* * §3-3 等参数单元 为降低离散化带来的误差，更好地反映弹性体中的位移状态和应力状态，可以采用较精密的单元(具有较高次位移函数的单元) 常用的较精密单元： 六节点三角形单元，矩形单元 在矩形单元上，取四个角点作为节点，矩形沿x，y轴的边长分别为2a，2b 1.位移函数 4节点，8自由度，列出8个方程，求解?1～ ?8 §3-3 等参数单元 位移函数是否满足有限元解收敛于精确解？ ①满足刚体位移条件 ②满足常应变条件 ③保证相邻单元公共边的位移协调 如12边，y=const 23边，x=const 边界位移变化是线性变化，两节点确定，位移变化惟一确定，满足变形协调条件 §3-3 等参数单元 形函数表示 §3-3 等参数单元 矩形单元的应力矩阵不是常应力矩阵，即应力分量不是常量。在弹性体中采用相同数目的节点时，矩形单元的精度高于三节点三角形单元。虽然相邻的矩形单元在公共边界处的应力有差异，但差异较小 矩形单元缺点： 不能适应斜交边界和曲线边界；不便在不同部位采用大小不同的单元 弥补： 矩形单元和三角形单元混合使用——单元边界上的位移线性变化，可以保证收敛于精确解的充分条件。但程序编制、单元信息填写比较复杂 §3-3 等参数单元 改用斜四边形单元，4节点8自由度，仍采用矩形单元的位移函数，在公共边界上，如12边的直线方程为 位移为x的二次函数，位移在12边上不是线性变化，不能由两个节点的坐标值惟一确定，不能保证在相邻单元的公共边上位移连续，即收敛于精确解的充分条件不能满足 如何解决矛盾，利用坐标变换，就是引入局部坐标系，把斜四边形单元变换成相应的矩形单元 一、等参数单元概念 提高有限元解精度途径： ①单元形态不变，增加单元数，减小单元尺寸 ②在单元数目不变的条件下，增加单元节点数目，增加形函数的阶次 对三角形单元来说，增加单元节点数目 a. 增加边上高阶节点构造二次三角形单元 b. 增加角节点，形成四边形单元 一、等参数单元概念 1．矩形单元 把坐标原点取在矩形单元中心,且边长为2;对在整体坐标中中心在p0(x0,y0)的矩形单元，边长为2a ，2b 通过坐标变换化成边长为2矩形单元 斜四边形单元有四个节点，仍采用有8个待定系数的位移模式。但和矩形单元不同，两个任意四边形单元的公共边上位移不是线性变化，位移的连续性得不到保证 一、等参数单元概念 一、等参数单元概念 解决办法——坐标变换 图中斜四边形单元，作连接对边中点的直线，作为?轴和?轴，取其交点为原点，令四边上的?值、 ?值分别为? 1，使整体坐标(x，y)下的斜四边形，变换为局部坐标(? ，?)下的矩形。对于局部坐标下的矩形单元，有 一、等参数单元概念 根据等参元性质，坐标插值 一、等参数单元概念 求应力、应变需要局部坐标和整体坐标之间的变换关系 已知 对任意四边形的四个节点成立，证明在边上，也正确 一、等参数单元概念 以?=c为例， ?是常数，必有 在整体坐标下与之对应一条直线，两个端点对应，对于四条边也正确 (xi,yi)是整体坐标，Ni是局部坐标(?,?)的函数 坐标变换式 一、等参数单元概念 利用坐标变换式将单元内任一点的整体坐标换算成局部坐标的函数，或是将局部坐标换算成整体坐标 整体坐标中任意四边形单元就可以用局部坐标中正方形单元的位移模式及形函数进行计算 一、等参数单元概念 应力分析需要整体坐标下的形函数，怎样求解？ 已有 若 求z对(?,?)偏导数 一、等参数单元概念 整体坐标和局部坐标变换关系 令 一、等参数单元概念 得到 [J]——雅可比(Jacobi)矩阵 条件： 雅可比矩阵的行列式值不能为零 位移插值和坐标插值采用同样的插值函数——等参数单元或等参元 在等参数单元的特性分析中被积函数比较复杂，通常需要采用数值积分方法进行计算。一般的有限元程序中都已编入这类单元，计算时只要选用即可 二、二维8节点等参数单元 在二维曲边单元中，8节点单元应用最广，即每边增加一个中间节点 图示边长为2的8节点正方形单元——母单元(基本单元)。单元各节点坐标(?i,?i)分别为?1或0 二、二维8节点等参数单元 二维8节点单元，单元有16个自由度，位移模式也有16个待定系数 局部坐标中的位移函数 当?为常数时，位移是?的二次函数；当?为常数时，位移是?的二次函数——双二次函数 若在母单元中选取双二次位移函