一类计算弹塑性问题的二阶锥互补法.pdfVIP

一类计算弹塑性问题的二阶锥互补法1) *,+, 2) + 李建宇 ，张洪武 * + （天津科技大学 机械工程学院，天津 300222 ） （大连理工大学 工业装备结构分析国家重点实验室，大连 116023) 摘要：基于最大塑性耗散原理和二阶锥规划对偶理论，建立线性等向和运动强化J2 弹塑性问题的二阶锥线性互 补模型；利用参变量变分原理和有限元方法，将混合线性强化J2 弹塑性问题转化为标准的二阶锥线性互补问题， 并利用半光滑Newton 法求解。新的二阶锥线性互补法是以往弹塑性分析线性互补法的推进，它不仅保持了线 性互补法数学模型的简洁性，而且避免了线性互补法对非线性屈服函数的显式线性化需求，因而问题的规模更 小，而解的精度更高。数学规划领域关于锥规划问题内点法、非内点法的研究，为二阶锥互补问题的高效求解 提供了强大算法支撑。。 关键词：J2 弹塑性; 参变量变分原理; 二阶锥; 互补; 半光滑Newton 法 引 言 随着塑性力学 （尤其是经典弹塑性力学）的基本理论日臻完善，在工程分析和设计中更多的考虑材 料的塑性性能成为可能。由于弹塑性力学控制方程的高度非线性，研究适用于高效求解大规模工程实际 [1] 问题的弹塑性数值方法成为必要。Simo和Taylor于1985年提出的回映 （Return-Mapping）算法 被公 [2] [3] 认为计算塑性领域里程碑式的工作，他们在Wilkins 和Krieg等 工作的基础上，引入了一致切线模 量 （或称算法模量），从而获得了具有超线性收敛的求解算法。回映算法的主体一般由两部分构成：全 局的平衡迭代和局部的本构迭代。其中局部本构迭代实质是一种弹性预估/塑性校正算法。所谓“回映” 是指若弹性预估应力状态在塑性屈服面以外，则需要回映到屈服面上。回映算法在本构方程迭代阶段， 首先需要识别出哪些区域的应力状态进入了塑性流动状态，然后利用Newton法求解由进入流动区域的屈 服条件所构成的非线性方程组。在平衡迭代阶段，基于回映后的应力增量求得的切线模量称为一致切线 模量，或算法模量。一致切线模量的获得，极大提高了平衡迭代Newton法的收敛速度。截至目前，回映 [4，5] 算法已经是求解包括金属弹塑性、岩土塑性等工程塑性问题的一种标准算法而被广泛应用 ，多数具 有非线性分析功能的商用CAE软件，如Ansys、Abaqus等均装配了该算法。 事实上，回映类算法并非求解弹塑性问题的唯一有效方法。规划法是几乎与回映算法同时发展起来 的另一类弹塑性计算方法。与回映类算法需要“人工”识别哪些区域进入塑性流动状态不同，规划法在应 力更新阶段利用数学规划领域提出的各类处理不等式约束的模型和算法“自动”识别塑性流动区域。例如 [6] [7] [8] [9] [10] [11] 增广Lagrangian法 ，线性互补法 ，非线性互补法 ，非光滑Newton法 ，SQP 法 以及双势方法 1) 国家自然科学青年基金资助项目 2) E-mail ：lijianyu@ 1 等。与回映类算法比较，规划法的优势不仅在于其简洁、优美的数学模型，而且数学规划领域在算法方 面的不断进展，使其在求解大规模工程实际问题方面更具潜力。特别是在 1994年到2004年间，人们在 将线性规划原对偶内点法推广到锥规划时意外地获得了多项式时间方法，进而在数学规划界出现了一种