等参单元和数值积分.ppt

数值积分方法 二维和三维高斯积分 在 ?，?，? 不同坐标方向上，也可以选不同的积分点数，即根据具体情况采用不同阶的积分方案。 二维 三维 等参元计算中数值积分阶次的选择 保证积分精度 数值积分阶次的选择影响计算的精度和工作量 一维问题刚度矩阵的积分 选择不当，会导致计算失败。 积分阶次的选择原则 插值函数 N 中多项式的阶次 p 微分算子 L 中导数的阶次 m 有限元中被积函数的阶次 2(p-m) 对于等参元假设 为常数 高斯积分的阶次 n=p-m+1 精确积分的多项式阶次 n=2(p-m)+1 满足刚度矩阵精确积分要求 * * 第 4 章等参单元和数值积分 等参变换的概念和实现单元特性矩阵变换的内容和方法。 本章要点 实现等参变换的条件和等参元满足有限元收敛准则的条件。 数值积分的基本思想及以高斯积分为代表的几种常用数值积分方法的特点。 第 4 章 单元和插值函数的构造 刚度矩阵数值积分阶次选择的原则，以及保证这些原则实现的具体方案。 等参单元 采用等参变换的单元 等参变换的概念和单元矩阵的变换 局部坐标单元 几何形状规则 几何形状规则的单元离散几何形状复杂的求解域比较困难 规则形状的单元 边界为曲线或曲面的单元 转化 等参变换 等参变换 单元几何形状的变换和单元内的场函数采用相同数目的结点参数及相同的插值函数的变换 整体坐标单元 几何形状扭曲 坐标变换 等参变换 等参变换的概念和单元矩阵的变换 等参变换 等参变换的概念和单元矩阵的变换 等参变换 等参变换的概念和单元矩阵的变换 坐标变换表示成插值函数的形式 m 用以进行坐标变换的单元结点数 形状函数，也是局部坐标表示的插值函数 这些结点在总体坐标的坐标值 函数插值 形式相同 坐标变换和函数插值采用相同的结点，并且采用相同的插值函数，即 等参变换 超参变换 坐标变换结点数多于函数插值的结点数 亚参变换 坐标变换结点数少于函数插值的结点数 单元矩阵变换 等参变换的概念和单元矩阵的变换 有限元中的单元体积内和面积内的积分 希望在自然（局部）坐标内按规格化进行数值积分 导数之间的变换 G 和 g 中常包含场函数对于整体坐标 x,y,z 的导数 单元矩阵变换 等参变换的概念和单元矩阵的变换 导数之间的变换 J 称为 Jacobi 矩阵，可表示为 单元矩阵变换 等参变换的概念和单元矩阵的变换 导数之间的变换 是 J 的逆矩阵，可按下式计算 单元矩阵变换 等参变换的概念和单元矩阵的变换 体积微元、面积微元的变换 d?,d?,d? 在笛卡儿坐标系内所形成的体积微元是 i, j, k 在是笛卡儿坐标 x, y, z 方向的单位向量 等参变换的概念和单元矩阵的变换 体积微元、面积微元的变换 在?=常数 (c) 的面上 面积微元 在面上的 dA 可以通过轮换 ?，?，? 得到 单元矩阵变换 单元矩阵变换 等参变换的概念和单元矩阵的变换 有限元中的单元体积内和面积内的积分变换到自然坐标系的规则化域内，分别表示成 （在 的面上， ） 其中 体积微元、面积微元的变换 单元矩阵变换 等参变换的概念和单元矩阵的变换 Jacobi 矩阵 体积微元、面积微元的变换 二维问题 单元矩阵变换 等参变换的概念和单元矩阵的变换 二个坐标之间的偏导数关系 体积微元、面积微元的变换 二维问题 d? 和 d? 在笛卡儿坐标系内所形成的面积微元 在?=常数 (c) 的曲线上，在笛卡儿坐标内的线段微元的长度 单元矩阵变换 等参变换的概念和单元矩阵的变换 面积或体积坐标不完全独立 面积（体积）坐标与笛卡儿坐标之间的变换 重新定义自然坐标 三维问题 面积或体积坐标不完全独立 单元矩阵变换 等参变换的概念和单元矩阵的变换 面积（体积）坐标与笛卡儿坐标之间的变换 二维问题 重新定义自然坐标 单元矩阵变换 等参变换的概念和单元矩阵的变换 面积（体积）坐标与笛卡儿坐标之间的变换 积分限的改变 L3=1 L3=0 2 3 4 在 L1=0 的表面上 在 34 边上， L2=0 在 23 边上， L1 +L2+L3+ L4 =1，L1=0, L4=0，即 L2=1-L3 类似地，三维四面体1234，在 L4=0 的表面上有L1=1- L2 -L3 等参变换的条件 等参变换的条件和等参单元的收敛性 等参变换为一种坐标变换，其一一对应的条件是 Jacobi 行列式 不得为0。 体积微元 三维问题 面积微元 二维问题 自然坐标的微元 不是一一对应 笛卡儿坐标的一个点 不成立 等参变换的条件 等参变换的条件和等参单元的收敛性 笛卡儿坐标中 笛卡儿坐