多尺度晶体塑性模拟文章推荐 原创

文章名称《Concurrent multi-scale crush simulations with a crystal plasticity model》

DOI：10.1016/j.tws.2011.12.019

在汽车防撞梁、吸能盒和薄壁管结构中，压溃吸能能力直接影响结构安全性。传统有限元分析通常采用各向同性塑性模型，通过宏观应力–应变曲线描述材料响应。但实际金属材料并不是“均匀黑箱”：晶粒取向、滑移系激活、织构演化都会影响局部塑性变形，尤其在薄壁管压溃这类大变形、强局部化问题中，微观结构可能对吸能行为产生重要影响。

压溃模型示意图：

Najafi、Marin 和 Rais-Rohani 的文章《Concurrent multi-scale crush simulations with a crystal plasticity model》正是针对这一问题展开研究。作者关注的是：在方形薄壁管轴向压溃过程中，材料初始织构及其演化是否会影响整体压溃力、平均吸能能力和局部折叠模式。

为解决这一问题，作者提出了一种并发多尺度建模方法：宏观结构层面采用显式有限元模拟方管压溃；每个积分点内部嵌入一个由多个 FCC 晶粒组成的多晶聚集体；晶粒层面采用 Marin 晶体塑性模型描述滑移、硬化和晶格旋转；最后通过 Taylor 型均匀化获得积分点平均应力。这样，宏观有限元计算不再只依赖经验塑性曲线，而是能够实时考虑晶粒取向和织构演化对结构响应的影响。

文章中，作者首先通过单元模型分别施加拉伸、压缩和简单剪切，生成不同初始织构；随后将这些织构赋予方管模型，并进行轴向压溃模拟。

结果表明，虽然不同织构对整体折叠形貌的影响并不总是非常显著，但对压溃力–位移曲线、平均压溃力和能量吸收能力具有明显影响。尤其是在角部、水平铰线和锥面等局部大塑性区域，晶粒取向会持续演化，形成不同的局部织构模式。文章还指出，拉伸织构和压缩织构在不同压溃模式下表现出不同的吸能优势，这说明“材料制造历史”并不是可以忽略的背景信息，而是可能影响结构服役性能的重要因素。

这篇文章对我们的启发在于：晶体塑性并不只能用于单晶拉伸、RVE 或微观变形分析，也可以嵌入显式动力学框架，用于研究真实工程结构中的局部变形、吸能和织构演化。对于高温合金、铝合金薄壁件、微尺度构件等问题，如果材料存在明显织构或晶粒尺度效应，将晶体塑性与结构有限元耦合，能够提供比传统本构更丰富的物理信息。

我们可以将我之前推文提到的umat-taylor模型转化为vumat子程序，进一步使用晶体塑性模型模拟大变形结构尺度材料变形行为。案例展示如下：

初始模型参考文章的设置（上下两层钢板，中间为薄壳结构）：

使用通用接触，摩擦系数设置为0.5，共4000个单元，每个单元包含50个具有不同初始取向晶粒。共20万晶粒。

边界条件设置为下端钢板固定，上端下压。

模拟结果如下：

应力分布结果：

晶粒1的剪切滑移：

晶粒2的剪切滑移：

晶粒50的剪切滑移：

单元标号5变形结束后的50个欧拉角分布：