增强假设变形梯度有限元方法稳定性计算.pdfVIP

VoL27 No.5 第 27 卷第5 期 北京科技大学学报 2005 年 10 月 Journal of University of Science and Technology Beijing Oct. 2005 增强假设变形梯度有限元方法稳定性计算 陈章华 刘洪波马文江 北京科技大学应用科学学院数力系，北京 100083 摘 要 讨论了非协调增强假设变形梯度奋限元的基本理论、所采用的本构关系和有限元 列式.为了克服原先由Simo建议的方法中所固有的缺陷，即在模拟压缩变形时容易出现奇异 能量模式而导致计算失稳甚至崩溃，对原方法中有关内变量和应力更新算法进行了改进.计 算结果表明，改进的方法在模拟压缩和拉伸模型时，具有较好的计算稳定性. 关键词 大变形:增强假设应变有限元:有限单元法:应力光滑 分类号 024 1.82; TB 115 在塑性成型分析中，有限元法可以有效地处 如颈缩(轴对称问题)或剪切带(平面应变问题) , 理塑性成型中材料、几何、边界三重非线性问题， 从而不会产生符合实际需要的极限载荷.而 是当前的发展趋势[!)然而，常规低阶协调单元存 EASM 单元可以通过采用一附加的应变场，即增 在着计算精度差、过于刚硬、容易自锁等缺点四) 强应变场，来增加单元的自由度数，从而提高单 为解决此问题，很多学者开展了深入研究.其中， 元柔性和适应变形局部化的特点.EASM 单元将 经典的Wilson 非协调元[3) 应用在小变形计算时具 应变场划分为两部分，协调应变场Vu和增强应 变场豆，即: 有很好的性能，但不能直接用于几何非线性问 e=Vu+主(1)， 题.目前又发展了许多低阶高精度非协调单元， 如增强假设应变有限元(Enhanced Assumed Strain 然而，增强应变场主不可任意选取.为满足具 Modes)，简称 EASM. 增强应变列式最初由 Simo 有三个独立场的Hu-Washizu广义变分原理，必须 于 1990 年因无限小问题提出I叫， 1992 年在广义变 在单元内部先对应力场σ施加相对于增强假设应 分原理的基础上给出了双线性四边形及三角形 变场主的L正交变化条件: 模式[旬，在模拟几何非线性问题时，具有较好的特 (σ，趴=fσ.edV=O (2) 性. 1993 年， Simo 进一步发展了一个用于三维问 进一步，为了不使作为独立变量的假设应变出现 题的无自锁单元QMl尼1 ，并且采用修正的积分 在总刚度方程中，独立变量需要在单元内部消 准则，但在压缩变形时仍存在Hourglass 现象 [6) 去，即将独立的增强应变在单元一级用静态凝聚 1996 年Aremero 等提出增强变形梯度场[咽，其优