基于光滑位移场的应力恢复方法.pdfVIP

第24卷第5期 计算力学学报 VOI．24。No．5 Chinese of Journal Mechanics October2007 2007年lO月 Computational 基于光滑位移场的应力恢复方法 李 涛。， 左正兴， 廖日东， 蒲大宇， 王舒妍 (北京理工大学机械与车辆工程学院．北京100081) 摘要：为了提高位移元的应力精度．提出了一种应力求解的新方涪．首先在一取定的单元小片内对基本解进 行多项式拟台光滑处理，使得基本解在单元交界处获得更高的光滑性，进而通过微分运算获得节点的应力值}最 后按此方法遍历所有的节点得到全求解域应力场．数值算例表明该方珐是可行的，并且计算结果的精度较高． 该方法计算量小t昔适性强，具有较殛的理论意义与工程应用价值． 关鼙词：应力恢复f超收敛}单元片I光滑位移场 中国分类号：0242+21文献标识码：A 是残值法与SPR法的结合．数值计算发现。sPR 1引言 法和SPRE法在给定边界条件的边界及其临近区 域误差较大，因此文献[9]提出SPREB方法，使得 位移元的应力精度损失是由微分运算引起的。 恢复应力满足边界条件。由于某些类型的单元并 文献[1]将该问题列为有限元法的五个疑难问题之 不确切地存在这类超收敛点，REP[6]方法随之被提 一(作者提出的解决方法为应用Hellinger-Reiss— 出．等等口]．对诸如上述的应力恢复方法一一研究 tier混合变分原理，用宏观的方法代替微分运算)。 后发现，它们有一个共同的特点，那就是按照如下 Lagrange元基本解在全域的G连续性使得基本 途径进行； 解在单元内无穷光滑，而在单元交界处经典导数不 存在，应力在单元之间不连续．对此，世界范围内 有限元基本解j墅}旦不连续应力解苎!坐生较好的应力群 的学者进行了大量的研究．研究了许多应力恢复方 它们的不同之处只是光滑应力的方法不同．在本 法．单元应力或节点应力平均是最原始、最简单的 文中，将交换“微分运算”与“光滑处理”的先后次 一种方法。计算量也比较小．高斯点应力外插平均 序，使有限元后处理过程按如下步骤进行： 能够获得较节点平均更高的应力精度．最小二乘 有限元基本_』塑坚光精位移场!堡墨兰．较好的应力解 意义下的总体应力磨平或在限制在单个单元内的 本文按上述方法进行应力恢复研究，然后结合数值 单元应力磨平能够获得较高精度的解，但计算量比 算例验证该方法的可行性，获得了较高的计算精度。 较大。1992年，基于某些类型的单元存在应力超 2通过光滑位移场获取应力场 收敛点，0’C．Zienklewiez，J．z．zhu[3∞提出了 Patch SPR(SuperconvergenceRecovery)应力恢 复法。该法将相邻单元组成的小片作为研究对象， 每一个连续函数都可以表示成某一多项式序列的 用特定阶效的多项式在单元小片内对应力场进行 一致收敛极限。根据函数逼近理论对其进行推论．