一种新的基于余能原理的应力恢复方法.pdfVIP

Vbl．33 第33卷第6期 燕山大学学报 No．6 2009年11月 JournalofYanshan Nov．2009 University 文章编号：1007-791X(2009)06-0505-05 一种新的基于余能原理的应力恢复方法 夏怀孝‘ (三峡大学土木水电学院，湖北宜昌443002) 摘 要：提出了一种新的基于余能原理的应力恢复方法：首先将单元小片从弹性体中隔离出来，然后构造一个 高次的待定平衡应力场，应用余能原理通过解一个最小二次规划问题来确定这个应力场，得到恢复应力。本方 法不需要单元中存在超收敛点，直接利用有限元计算结果构造满足平衡条件的应力场，方法简便可靠。最后通 过一个经典的线弹性算例验证了方法的有效性。 关键词：有限元；应力恢复；余能原理；二次规划；小片 中图分类号：034 文献标识码：A 0引言 在超收敛点，在这些点应力具有较高精度这一事 实，用最小二乘法来拟合这些超收敛点的应力，构 由最小势能原理导出的有限元计算格式通常 造一个由相邻单元组成的小片上的光滑应力场。随 得到不连续的应力结果，而且应力的精度一般低于 后一些改进的SPR方法被提出，这些方法将平衡 位移，研究者为了从有限元计算结果中得到较精确 方程和应力边界条件引入，来提高SPR方法的精 的应力，发展了许多被称为应力恢复的方法，显然 度[5-7101997年另外一种具有超收敛性质的应力恢 这样的方法有着重要的工程意义。 复方法recovery by 从20世纪80年代起，自适应有限元成为研究种方法不需要单元存在超收敛点，该方法通过施加 的热点，不管是h型或者p型自适应有限元都包括一 弱的平衡条件来获得恢复应力㈣。 个重要的部分即误差估计“1。目前使用最普遍的误 利用最小余能原理来改进有限元应力精度的 差估计方法是由Zienkiewicz-Zhu提出的z2方法跚， 研究也一直在进行，文献[10-11】是其中较突出 该方法将应力恢复得到的应力代替原问题的精确 的，这些方法的显著特点是不需要单元存在超收敛 解进行误差估计的计算，很显然这种误差估计方法 点，这对于某些有限元单元(三角形单元或者空间 的有效性直接取决于应力恢复的精度，Zienk- 单元)具有特殊意义。这些文献采用的方法都是将 iewicz·Zhu曾证明如果应力恢复得到的应力是超 由相邻单元组成的小片或者某点的邻域从弹性体 收敛的，那么Z2得到的结果总是近似正确的Ⅲ。正 中隔离出来，将有限元得到的位移施加于小片或者 是因为应力恢复方法有着上述重要的理论