土石混合体细观数值模拟的弹塑性自适应有限元分析.pdfVIP

第４２卷第２期应用力学学报Ｖｏｌ.４２Ｎｏ.２

２０２５年４月ＣｈｉｎｅｓｅＪｏｕｒｎａｌｏｆＡｐｐｌｉｅｄＭｅｃｈａｎｉｃｓＡｐｒ.２０２５

ＤＯＩ:１０.１１７７６/ｊ.ｉｓｓｎ.１０００￣４９３９.２０２５.０２.０１５

土石混合体细观数值模拟的弹塑性

自适应有限元分析

余科ꎬ肖映雄ꎬ徐亚飞

(湘潭大学机械工程与力学学院ꎬ４１１１０５湘潭)

摘要:在土石混合体的细观弹塑性有限元分析中ꎬ为了能有效克服不规则块石附近的应力集中以

及更好地捕捉塑性区大的塑性应变、提高塑性区的计算精度ꎬ需要采用自适应网格技术ꎬ以避免加

密网格时的盲目性ꎮ采用基于后验误差估计的自适应网格重划分算法ꎬ并结合Ａｂａｑｕｓ二次开发ꎬ

编写了相应的弹塑性自适应有限元Ｐｙｔｈｏｎ脚本ꎬ实现了土石混合体细观数值模拟的自适应求解全

过程ꎮ通过将Ｐｙｔｈｏｎ脚本应用于几类典型的二维土石混合体细观结构模型的数值模拟ꎬ验证了基

于网格重划分算法的自适应有限元方法在克服应力集中及更好捕捉塑性区大的塑性应变、提高塑

性区计算精度上的有效性和适应性ꎮ

关键词:细观随机模型ꎻ自适应网格重划分ꎻ等效塑性应变ꎻ应力集中ꎻ数值模拟

中图分类号:Ｏ３４４.３ꎻＯ３０２文献标志码:Ａ文章编号:１０００￣４９３９(２０２５)０２￣０３９５￣１０

Ｅｌａｓｔｉｃ￣ｐｌａｓｔｉｃａｄａｐｔｉｖｅｆｉｎｉｔｅｅｌｅｍｅｎｔａｎａｌｙｓｉｓｉｎｍｅｓｏｓｔｒｕｃｔｕｒａｌ

ｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｏｆｓｏｉｌ￣ｒｏｃｋｍｉｘｔｕｒｅ

ＹＵＫｅꎬＸＩＡＯＹｉｎｇｘｉｏｎｇꎬＸＵＹａｆｅｉ

(ＳｃｈｏｏｌｏｆＭｅｃｈａｎｉｃａｌＥｎｇｉｎｅｅｒｉｎｇａｎｄＭｅｃｈａｎｉｃｓꎬＸｉａｎｇｔａｎＵｎｉｖｅｒｓｉｔｙꎬ４１１１０５ＸｉａｎｇｔａｎꎬＣｈｉｎａ)

Ａｂｓｔｒａｃｔ:Ｉｎｔｈｅｍｅｓｏｓｃｏｐｉｃｅｌａｓｔｉｃ￣ｐｌａｓｔｉｃｆｉｎｉｔｅｅｌｅｍｅｎｔａｎａｌｙｓｉｓｏｆｓｏｉｌ￣ｒｏｃｋｍｉｘｔｕｒｅꎬｉｔｉｓｎｅｃｅｓｓａｒｙｔｏａ￣

ｄｏｐｔａｄａｐｔｉｖｅｍｅｓｈｅｓｔｏａｖｏｉｄｔｈｅｂｌｉｎｄｎｅｓｓｏｆｍｅｓｈｒｅｆｉｎｉｎｇｏｒｒｅｍｅｓｈｉｎｇｉｎｏｒｄｅｒｔｏｏｖｅｒｃｏｍｅｔｈｅｓｔｒｅｓｓ

ｃｏｎｃｅｎｔｒａｔｉｏｎｎｅａｒｔｈｅｉｒｒｅｇｕｌａｒｂｌｏｃｋｂｌｏｃｋｓａｎｄｂｅｔｔｅｒｃａｐｔｕｒｅｔｈｅｌａｒｇｅｐｌａｓｔｉｃｓｔｒａｉｎｉｎｔｈｅｐｌａｓｔｉｃｚｏｎｅｓ

ａｎｄｆｕｒｔｈｅｒｉｍｐｒｏｖｅｔｈｅｃａｌｃｕｌａｔｉｏｎａｃｃｕｒａｃｙｏｆｔｈｅｐｌａｓｔｉｃｚｏｎｅｓ.Ｉｎｔｈｉｓｐａｐｅｒꎬｔｈｅａｄａｐｔｉｖｅｒｅｍｅｓｈｉｎｇａｌ￣

ｇｏｒｉｔｈｍｂａｓｅｄｏｎａｐｏｓｔｅｒｉｏｒｉｅｒｒｏｒｅｓｔｉｍａｔｉｏｎｉｓａｄｏｐｔｅｄａｎｄｔｈｅｒｅｓｕｌｔｉｎｇＰｙｔｈｏｎｓｃｒｉｐｔｏｆｔｈｅｅｌａｓｔｉｃ￣ｐｌａｓ￣

ｔｉｃａｄａｐｔｉｖｅｆｉｎｉｔｅｅｌｅｍｅｎｔｐｒｏｇｒａｍｉｓｗｒｉｔｔｅｎｃｏｍｂｉｎｉｎｇｗｉｔｈＡｂａｑｕｓｓｏｆｔｗａｒｅ.Ｔｈｅｗｈｏｌｅｐｒｏｃｅｓｓｏｆａｄａｐ￣

ｔｉｖｅｓｏｌｕｔｉｏｎｉｓｒｅａｌｉｚｅｄｆｏｒｍｅｓｏｓｃｏｐｉｃｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｓｏｆｓｏｉｌ￣ｒｏｃｋｍｉｘｔｕｒｅ.Ｂｙａｐｐｌｙｉｎｇｔｈｅｐｙｔｈｏｎ

ｓｃｒｉｐｔｔｏｔｈｅｍｅｓｏｓｔｒｕｃｔｕｒｅｎｕｍｅｒｉｃａｌ