第5章有限单元法应用中的若干实际考虑.ppt

* CHONGQING UNIVERSITY CHONGQING UNIVERSITY 工程力学系 第5章 有限单元法应用中的若干实际考虑 主讲教师：严 波 CHONGQING UNIVERSITY CHONGQING UNIVERSITY 工程力学系 5.1引 言 1 计算应力结果的处理和改善 应力和应变精度低于位移（位移元） 2 减少模型自由度，降低计算量 子结构法、结构对称性和周期性的利用 3 Wilson非协调元 改进等参元的计算精度和效率 CHONGQING UNIVERSITY CHONGQING UNIVERSITY 工程力学系 5.2 有限元模型的建立 1 单元类型和形状的选择 1）单元类型 物理属性： 固体力学 流体力学 热传导 压电 电磁场 固体力学： 二维实体 三维实体 梁 板 壳 无限元 2）单元形式 不同形状 不同阶次 CHONGQING UNIVERSITY CHONGQING UNIVERSITY 工程力学系 2 网格划分 1）网格疏密的合理布置 应力集中区 （a）密集网格 （b）局部分析方法 （c）自适应分析方法 2）疏密网格的过渡 （a）过渡单元（变结点单元） （b）不规则形状单元过渡（精度受影响） （c）三角形单元（处理不便时实用） （d）多点约束方法 罚函数法：近似处理 在方程中引入约束方程 CHONGQING UNIVERSITY CHONGQING UNIVERSITY 工程力学系 5.3 应力计算结果的性质和处理 1 位移元应力解的误差 （1）单元内部不精确满足平衡方程 （2）单元交界面上应力不连续 （3）力的边界上一般不精确满足边界条件 在建立平衡方程时，泛函的位移变分并不是完全任意的。 2 近似解的性质 CHONGQING UNIVERSITY CHONGQING UNIVERSITY 工程力学系 求 的极值归结为求 的极小值。 对线弹性体，上式还可以表示为 求 的极值 与 的差值 满足加权的最小二乘 CHONGQING UNIVERSITY CHONGQING UNIVERSITY 工程力学系 最小二乘的概念 测量数据 寻找函数 使得 最小，得到离散数据得最佳逼近。 x1 xn x y 注意与插值函数的区别。 最小二乘曲线拟合 结论： 有限元的应变近似和应力近似解必然在真解的上下震荡 单元内存在（等于真解的）最佳应力点。 CHONGQING UNIVERSITY CHONGQING UNIVERSITY 工程力学系 3 等参元的最佳应力点 为 p 次多项式，L为 m 阶微分 为 n=p-m 次多项式 上述积分表达式的积分阶次：2n 取n+1阶Gauss积分，精度：2(n+1)-1=2n+1次 如果 且 各分量独立，应力和应变在Gauss积分点上 是精确的。 结论：在等参元中，单元中（n+1）阶（n=p-m）Gauss积分点上的应力或应变近似解比其它部位具有较高的精度，称为最佳应力点。 CHONGQING UNIVERSITY CHONGQING UNIVERSITY 工程力学系 图5.4 注意：结点上的应力和应变值的精度较差。 图5.5 CHONGQING UNIVERSITY CHONGQING UNIVERSITY 工程力学系 4 单元平均或结点平均 1）单元平均 ① ② 常用于3结点三角形单元 相邻两单元构成的四边形形心处： （平均） 或 （面积加权平均） 2）结点平均 m是围绕单元 i 的全部单元数。 CHONGQING UNIVERSITY CHONGQING UNIVERSITY 工程力学系 5 总体应力磨平 构造一改进的应力解 此改进解在全域上连续 建立泛函： 插值函数 待求改进后的结点应力 代入 应力磨平的所有单元结点总数 CHONGQING UNIVERSITY CHONGQING UNIVERSITY 工程力学系 CHONGQING UNIVERSITY CHONGQING UNIVERSITY 工程力学系 6 单元应力磨平 令C＝I，未加权的最小二乘 对于等参单元，可利用精度较高的Gauss积分点的应力值来改变结点应力的近似性质。 例子： 4结点四边形单元 1 2 3 4 Ⅰ Ⅱ Ⅲ Ⅳ CHONGQING UNIVERSITY CHONGQING UNIVERSITY 工程力学系 可以推广到三维问题。 CHONGQING UNIVERSITY CHONGQING UNIVERSITY 工程力学系 7 分片应力磨平 总