极端加载条件下材料行为的数值仿真技术.pdfVIP

CAE自主创新发展战略＼学术研究水平 极端加载条件下材料行为的数值仿真技术 张 雄 (清华大学航天航空学院，北京100084) 摘 要在极端加载条件下，物质将发生很大的变形，导致网格畸变，给数值模拟造成很大的 困难．近年来发展起来的无网格法直接利用离散点构造近似函数，有效地克服了因网格畸变而 产生的求解困难，成功地应用于求解动态裂纹扩展、金属加工成形、高速碰撞、流体流动等问 题中，是分析极端加载条件下材料行为的有效方法。本文简要介绍了我们课题组近年来利用无 网格法研究极端加载条件下材料行为的一些工作。 关键词极端加载条件；有限元；无网格法；物质点法 一、引言 极端加载条件(如爆炸、超高速碰撞、金属加工成型等)涉及到特大变形、高 温高压和高应变率，实验研究难度大，成本高。数值模拟技术较易模拟各种工况和 材料，因而成为研究极端加载条件下材料行为的重要方法。 在极端加载条件下，物质将发生很大的变形，导致网格畸变，给拉格朗日有限 元法造成很大的困难。网格畸变不但严重降低计算精度，且使得显式积分的稳定时 间步长变得极小，大幅度地增加了计算量。另外网格畸变可能使单元雅可比行列式 变为负值，导致计算失败。欧拉法虽可以避免网格畸变的问题，但在求解方程中引 入了非线性对流项，给数值求解造成了困难，且难于准确描述物质界面。 任意拉格朗日．欧拉方法(ALE)[1l较好地克NT拉格朗日方法和欧拉方法的缺 点，吸收了各自的优点，．成功地应用于二维问题中。在ALE法中网格可以根据需 要独立于物质运动，以便在准确跟踪物质界面的同时能够维持合理的单元形状，但 在求解方程中仍然存在对流项。对于复杂的三维问题，为ALE设计有效的网格运 动算法也是非常困难的。 · 近年来发展起来的无网格法‘2】[31摒弃了网格，直接利用离散点来构造近似函数， ⑩ 学术研究水平／CAE自主创新发展战田各 有效地克服了网格畸变而产生的求解困难。目前已提出了十余种无网格法，成功地 应用于求解动态裂纹扩展、金属加工成形、高速碰撞、流体流动等问题中，是分析 极端加载条件下材料行为的有效方法。 本文简要介绍了我们课题组近年来利用无网格法研究极端加载条件下材料行 为的一些工作。 二、质点积分无单元伽辽金法 Galerkin Method)141采用背景网格覆盖求解 无单元伽辽金法(EFG,Element-Free 域，将对求解域的积分化为对各背景格子的积分之和，然后在每个格子中使用高斯 积分。EFG具有良好的稳定性，但因其近似函数不是多项式，需要采用较高阶的高 斯积分，计算量很大。Beissel等人采用节点积分方案来代替EFG中背景网格的积 ofthe 分，提出了节点积分无单元伽辽金法(NIEFG,Nodal Element—Free Integration Galerkin Method)t51。和有限元法中的单点高斯积分类似，节点积分无单元伽辽金法 也会出现因零能模态不受控制而产生的不稳定现象。Beissel等人建议将平衡方程残 差的平方引入到伽辽金弱形式中来消除节点积分的不稳定性。与背景网格积分相 比，节点积分可以大幅度提高伽辽金无网格方法的计算效率，但是也显著降低了计 算精度。在伽辽金弱形式中引入稳定项虽然可以有效地消除节点积分的不稳定性， 但对于不存在零能模态的问题，稳定项将大幅度地降低解的精度。 Chen等人提出用光滑应变稳定化的方法来消除节点积分的不稳定性【6】。该方法 通过节点邻域内应变的加权平均来计算该节点的应变值，能有效地消除零能模态导 致的不稳定性，显著地提高了节点积分的精度和稳定性。但是这种方法在计算刚度 矩阵的时候需要在各节点子域的边界上进行积分，增加了节点积分的计算量。