镁合金孪晶建模经典文章推荐 原创

文章名称：《A constitutive model for hcp materials deforming by slip and twinning: application to magnesium alloy AZ31B》

DOI：10.1016/S0749-6419(03)00039-1

在镁合金晶体塑性建模的发展脉络里，2003 年 Staroselsky 这篇文章并不是最新的，也不是形式最复杂的。但如果真正从“建立一个能跑、能解释现象、还能继续扩展的模型”这个角度去看，它仍然是一篇非常值得反复阅读的经典工作。

这篇文章研究的对象，是 AZ31B 镁合金在室温条件下的塑性变形行为。作者关注的问题非常明确：为什么这种 HCP 结构材料在不同加载方向下，会表现出强烈的不对称性、明显的织构演化，以及非常突出的孪晶效应？换句话说，这篇文章不是简单去拟合一条应力—应变曲线，而是试图回答：镁合金在室温下究竟是靠哪些机制在变形，这些机制又如何共同决定宏观响应。

这篇工作的建模核心思想：

第一，它非常强调孪晶不是附属机制，而是主导机制之一。在 HCP 镁合金里，单靠常规滑移并不能解释很多室温下的实验现象，尤其是压缩拉伸不对称和织构快速变化。Staroselsky 这篇文章抓住了这一点，把孪晶明确地放到了与滑移同等重要的位置。

第二，它采用了一种非常有启发性的处理方式：把孪晶作为一种具有方向性的 pseudo-slip 来处理。这样做的好处是，一方面保留了孪晶“有极性、不可正反完全对称”的物理特征，另一方面又能把它自然嵌入有限变形晶体塑性框架中。这个思想到今天看仍然非常高明，因为它在“物理真实性”和“程序可实现性”之间找到了很好的平衡。

第三，这篇文章并没有急着把硬化写得非常复杂。恰恰相反，作者采用了一个相对简洁的 non-hardening 假设，认为材料宏观上表现出来的许多“硬化感”，很大程度上其实来自取向变化和孪晶导致的晶格重排，而不只是每个滑移系的阻力不断增加。这个判断非常重要，因为它提醒我们：在 HCP 镁合金里，单纯盯着“硬化参数”往往是不够的，织构和重取向本身就是塑性响应的重要组成部分。

第四，作者还加入了一个各向同性 accommodation 项，用来描述晶界附近非晶体学协调变形的作用。这个处理非常值得重视。很多时候我们做晶体塑性，只把目光放在晶内滑移和孪晶上，却容易忽略多晶材料中晶粒之间并不是天然完全协调的。Staroselsky 这篇文章清楚地认识到：如果不考虑这部分效应，数值计算中的应力水平会偏高，甚至难以合理匹配实验。因此，这个附加项虽然形式上不复杂，但在建模思想上非常成熟。

主要考虑的滑移和孪晶如下：

拉伸变形的实验于模拟结果对比：

压缩变形的模拟和实验结果对比：

从结果上看，这篇文章得到的结论也非常有代表性。作者指出，AZ31B 镁合金室温变形的主要载荷承担机制，不是单一滑移，而是 basal、prismatic、pyramidal <a> 滑移与 {10-12}<10-11> 孪晶共同作用的结果。同时，他们认为材料在变形过程中表现出的显著各向异性和应变硬化特征，很大一部分来源于变形诱导的织构演化，而孪晶又是其中最关键的驱动力之一。这个结论今天看依然不过时，因为它直接抓住了镁合金室温塑性的本质。

推荐文章的最主要原因是：

2003 的这个文章的价值，不在于它把所有机制都做全了，而在于它先把最重要的几个问题讲清楚了——谁在变形、谁在重取向、谁在影响应力水平。对于刚开始做 HCP 晶体塑性的人来说，这种建模路径非常值得学习。

孪晶不能只作为“后处理现象”看待，而应该进入本构主框架。 如果模型里没有真正把孪晶放进去，那么很多镁合金室温响应就只能停留在表面拟合，难以解释深层原因。

一个好模型不仅要能解释，还要能扩展。Staroselsky 这篇文章的框架虽然简洁，但非常适合作为后续更复杂模型的基础。无论是进一步加入显式硬化、温度效应，还是发展到更复杂的 twin-detwin、parent-child 思路，它都可以作为一个非常扎实的起点。

后续很多孪晶模型基于此进行二次开发，因此实现该文章的数值模型对于孪晶的研究非常有帮助：

使用文章的公式，讲整体算法集成到abaqus的vumat子程序相对容易，因为不需要推导一致性雅可比。但是率无关模型通常数值稳定性较差。

使用该本构模型模拟效果如下：

初始RVE模型：

沿着X方式施加拉伸变形，变形结束后应力分布：

变形后的孪晶体积分数：

阈值孪晶体积分数（文章中推荐0.3-1.0的随机值）：

变形结束是否发生孪生变形：