有限元法基础-5等参元与数值积分答题.pptVIP

5.4 数值积分方法 有限元法基础 * 5.4 数值积分方法 有限元法基础 * 5.5 数值积分阶次的选择 有限元法基础 * 积分点个数的选取是对数值积分阶次的选择 计算精度 计算工作量 计算成本 5.5 数值积分阶次的选择 有限元法基础 * 选取积分点个数的原则 1）保证积分精度 2）保证总体刚度矩阵满秩 3）有较好的计算效率 5.5 数值积分阶次的选择 有限元法基础 * 1）保证积分精度 以一维单元刚度矩阵积分为例 积分限标准化，并设Jacobi行列式为常数 5.5 数值积分阶次的选择 有限元法基础 * 对多数弹性力学问题 Ni 插值函数： p 阶多项式 D 微分算子：最高导数 阶次 m 原被积函数为2（p － m）阶多项式 5.5 数值积分阶次的选择 有限元法基础 * 为保证积分精度，Gauss积分点数为n，应有 按此规则选取积分点个数，才能使被积函数达到精度。 5.5 数值积分阶次的选择 有限元法基础 * 对二维和三维单元 按一维的方法选 nxn 或nxnxn 个积分点，可能被积函数达不到精确积分的要求！ 原因1：Jacobi行列式可能不是常数, 这样提高了被积函数的阶次。 当物理坐标中的单元 平行四变形（2D） 平行六面体（3D） 5.5 数值积分阶次的选择 有限元法基础 * 解决办法 1）适当提高积分点数，以适应精度 2）剖分网格时，尽量避免过分扭曲单元 5.5 数值积分阶次的选择 有限元法基础 * 例：不同形状网格剖分的悬臂梁 5.5 数值积分阶次的选择 有限元法基础 * 5.5 数值积分阶次的选择 有限元法基础 * 原因2：B矩阵中包含有高阶非完全项 原插值函数：p阶完备多项式 p’阶非完全项 采用精确积分方案，应以p’为准，即 5.5 数值积分阶次的选择 有限元法基础 * 例：二维4节点单元 优化积分方案： p=1, n = p-m+1 = 1, 一点积分 非完全项含有 , 精确积分方案：积分点 2x2 5.5 数值积分阶次的选择 有限元法基础 * 例：二维8节点单元 优化积分方案： p=2, n = p-m+1 = 2, 2x2积分 精确积分方案： 非完全项含有 ， 积分点 3x3 5.5 数值积分阶次的选择 有限元法基础 * 在实际计算单元刚度矩阵时，还有其他方面的考虑。 实际的Gauss积分点数 精确积分要求的积分点数 5.5 数值积分阶次的选择 有限元法基础 * 精确积分是非完全项所要求，决定有限元的精度是完全多项式，选取较低的Gauss积分已保证完备项的要求，可提高计算效率 位移元得到的解是下限解，往往偏刚，降阶积分可减少单元的刚度 对于畸形单元，积分中含有 本身也不可能精确积分 5.5 数值积分阶次的选择 有限元法基础 * 2）保证总体刚度矩阵满秩 系统总自由度数 N＝系统节点数 X 每节点自由度－约束数 约束数 刚体位移数 施加边界条件后，总刚度矩阵非奇异 存在，方程有解。 5.5 数值积分阶次的选择 有限元法基础 * 矩阵的秩 1）矩阵相乘 2）矩阵相加 5.5 数值积分阶次的选择 有限元法基础 * 单元刚度矩阵的计算公式 C是dXd的方阵，d是应变数量，三维问题为6，平面问题为3，轴对称问题为4。 一般情况下，秩B＝d M个单元的结构 5.5 数值积分阶次的选择 有限元法基础 * 是K非奇异的必要条件。 系统独立的自由度数N超过(或 ) 全部积分点nG提供的独立关系数，则K必然奇异。 5.5 数值积分阶次的选择 有限元法基础 * 例：平面8节点单元 独立DOF数 N＝2x8-3=13 精确积分3x3 减缩积分2x2 减缩积分下K是奇异的。 5.5 数值积分阶次的选择 有限元法基础 * 2x2积分的平面8节点单元的零能变形模式 5.5 数值积分阶次的选择 有限元法基础 * 零能变形模式（Zero Energy Deformation Mode）： 一种由刚度矩阵产生的变形能为零的非刚体位移的变形模式