软化材料的弹塑性有限元分析.PDF

第３１卷　第３期 四 川 兵 工 学 报 ２０１０年３月  【其他研究】  软化材料的弹塑性有限元分析 夏怀孝 （三峡大学 土木水电学院，湖北 宜昌　４４３００２） 摘要：给出了软化材料弹塑性有限元分析的基本算法：选用显式Ｅｕｌｅｒ法作为本构积分方案；在结构未达到极限 荷载前采用初应力法配合Ｔｈｏｍａｓ加速方案作为平衡迭代方法；而当结构接近极限荷载时，程序转换为弧长法控 制迭代，使结构越过极值点进入软化区直至破坏。最后，讨论了该领域面临的重要问题———由于材料软化而导 致的应变局部化现象及其有限元模拟。 关键词：有限元分析；弹塑性；软化；应变局部化 中图分类号：Ｏ３４４．３ 文献标识码：Ａ 文章编号：１００６－０７０７（２０１０）０３－０１３０－０３ 　　近３０年来，弹塑性有限元分析一直是国际数值计算领 题，即如何求得极限荷载，以及如何越过该极值点（ｌｉｍｉｔ 域研究的重点。随着工程规模的扩大和相应研究的深入， ｐｏｉｎｔ）使受力—位移过程得以继续发展。这２个问题一直 结构的极限荷载及其越过极限荷载后的破坏情况等逐渐 是非线性有限元分析的难点，研究者也提出了许多方法进 引起了工程界的高度重视。因此，含软化阶段的弹塑性有 行解决，如虚拟弹簧法、位移控制法、强制迭代法等，关于 ［４］ 限元分析也随之被引入。 这个问题Ｂｅｒｇａｎ 也有详细评述。众多研究表明，Ｒｉｓｋｓ和 本构积分是弹塑性分析最重要的环节之一，它对计算 Ｃｒｉｓｆｉｅｌｄ等［５－６］提出的弧长法是较为有效的方法。该法在 结果的精度和计算的收敛性均有着重要影响。弹塑性本 每个增量步中，能通过弧长约束条件使迭代在位移—荷载 构积分方案大致有显式法和隐式法２种。隐式法由于能给 空间沿一弧形路径收敛至平衡位置，且该平衡位置位于上 出一致切线刚度矩阵（ｃｏｎｓｉｓｔｅｎｔｔａｎｇｅｎｔｓｔｉｆｆｎｅｓｓｍａｔｒｉｘ），且 一增量步平衡位置附近，这样就可稳定地越过极值点，从 对于ＮｅｗｔｏｎＲａｐｈｓｏｎ法可达到二阶收敛速度，因而该方法 而避免计算发散的危险。 ［１］ 得到了工程界的广泛应用 ；相比较而言，由于显示法较 针对上述问题，本文中根据计算经验选用了显式Ｅｕｌｅｒ 直接，也易于实现，并且还可通过误差控制的自适应本构 法作为本构积分方案。在结构未达到极限荷载前，采用初 积分来提高计算精度和效率，故近年也得到了新的 应力法配合Ｔｈｏｍａｓ加速方案作为平衡迭代的方法；而当结 发展［２］。 构接近极限荷载时，当迭代转换为弧长法控制，越过极值 平衡迭代是弹塑性非线性有限元分析的又一重要问 点进入软化区直至破坏。 题。平衡迭代就是将不平衡力施加于系统得出位移增量， 重复这个过程直至系统离平衡状态的差距可以忽略为止。 １　弹塑性本构积分 根据所采用线性方程组刚度矩阵的不同，平衡迭代法可分 为 ＮｅｗｔｏｎＲａｐｈｓｏｎ法、拟Ｎｅｗｔｏｎ－Ｒａｐｈｓｏｎ法和初应力法。 塑性屈服面在应力空间中可记为 Ｆ（，ｋ）＝０ （１） 关于各种方法及其比较有很多文献可供参考，这里不再一 σ 一列举。由于初应力法具有极强的稳定性和广泛的适用 　　随着塑性变形的发生