EquationChapter1Section1弹塑性力学结课报告有限元计算结果的.docVIP

弹塑性力学结课报告 有限元计算结果的平滑性处理 指导教师：景作军 专业：机研-07 姓名：李 昊 学号：2007摘 要 有限元法是最重要的工程分析技术之一。它广泛应用于弹塑性力学、断裂力学、流体力学、热传导等领域。有限元法有限元法 应变矩阵B是形函数矩阵付N对坐标进行偏导后得到的矩阵。偏导一次，位移模式的多项式降阶—次。所以通过偏导运算得到的应变ε和应力σ精度较位移u低，这在总体上表示的应力场中会出现缺乏平滑性从而降低大多数算法的强大功能，并可能在结果中引入不期望的波动现象发生。目前在选择导致平滑求解的方法和处理不连续的响应方面还缺乏完全规则化的统一的确定模式。这里讨论一些简单易行而又在一定范围内行之有效的方法，以期达到抛砖引玉的目的。 1.2 等参元的最佳应力点 众所周知，应用最小势能原理求得的位移解具有下限性质。也就是说一般近似解的总势能总是大于真解的总势能，而近似解的应变能总是小于真解的应变能，因此位移解小。 下面我们讨论由有限元求得的应力又将如何。如果位移、应变和应力的真解是 、、。．而近似解表示为、、，他们之间可表示为 按最小势能原理的泛函可表示为式（2.1） （1.1） 由于泛函变分要求=O，上式可表示为 上式表示求的极小值问题可归结为求的极小值问题，按有限元离散形式，可表示成泛函 （1.2） 显然，由（1.1），（1.2）式可改写为 （1.3） 式（1.2）和（1.3）的极小值问题，从数学上可以理解为求解应变值和应力值的加权二乘最小值问题。它显然可以表明，应变近似解和应力近似解必然是在真解左右振荡，即在单元内存在最佳应变点和最佳应力点，这个特点有助于我们对应力计算结果进行处理，从而提高应力解的精度。 由上述可知，如令离散形式的应力泛函式（1.3）记为，则有限元应力分析归结为对泛函求变分 （1.4） 或 （1.5） 假如近似解为p次多项式，是m阶微分算子，令n=p-m，则或是n次多项式。为使以上两式能够精确积分，至少应采用n+1阶高斯积分，被积函数是2n+1次多项式情况仍可达到精确积分。若每一个单元中的高斯积分点上的每个分量的变量是独立的，则必须保证 而就是高斯积分点上的应力近似解、它可以达到n+1次精度。因此选择等参元高斯点上的应力是具有较高精度的应力值，它比选用节点应力更为合理。 1.3 全域应力场的平滑性处理 采用有限元法求得的应力场往往不连续，而实际应力场的真解应该是连续的，这是近似解法带来的问题。为改进计算结果，有必要对应力进行平滑性处理，以便获得全域连续的应力场。 全域应力场的平滑性处理，有时也称总体应力磨平。其磨平原则是平滑性处理后的应力与有限元法求得的应力解应满足加权最小二乘法的规则，即 （1.6） 式中，m为求解域单元总数。 为获得平滑处理后的应力值，设它也可以表示为形函数形式 （1.7） 式中，为待求节点应力；为单元结点数，为形函数矩阵：和可以与求位移解时的单元节点数和形函数相同，也可以不同。 将式(1.7)代人式(1.6)，进行变分，由于的任意性，可得 （=1,2,…,l） 或 （1.8） 其中为进行应力平滑性处理时全部单元的节点总数。 求解式(1.8)的线性代数方程组，可得各节点平滑性处理后的应力，再代人式(1.7)，获得单元内经平滑处理后的应力场，它将能保持求解全域应力场的平滑性。 除此以外，平滑处理方法还有相邻单元的单元平均应力处理方法和单元节点平均应力简单处理、子域局部平滑处理、外推平滑处理、罚函数法平滑处理等。特别是在非线性有限元分析、摩擦、碰撞等问题中，应力平滑性处理显得尤为重要和突出，应用方法也较多。 1.4 外推插值有限元法 根据外推理论和有限元插值理论，对于二维的场问题，我们可以运用式（1.9）、（1.10）求得网格上的外推解，它具有比线性有限元解高一阶的精度。若在网格上，在以外推解进行插值有限元计算，得到