考虑位错滑移,trip,twip的高锰钢变形机制的晶体塑性建模
文章名称《A crystal plasticity model for twinning- and transformation-induced plasticity》DOI:10.1016/j.actamat.2016.07.032 高锰奥氏体钢有一个很有意思的特点:同样是在拉伸变形,不同温度下,它的微观塑性机制可能完全不同。低温时,奥氏体会发生 ε-马氏体相变,表现出 TRIP 效应;中等温度下,形变孪晶大量生成,带来 TWIP 效应;而高温时,这些机制又逐渐退场,普通位错滑移成为主角。今天推荐的文章来自马普所的一篇经典晶体塑性论文提出了一个统一的晶体塑性模型,同时描述位错滑移、形变孪晶和马氏体相变,用于描述高锰奥氏体钢的变形行为。对于塑性速度梯度,作者扩展了经典的速度梯度公式,纳入了孪生和相变的贡献:通过位错密度演化,模型可以描述加工硬化。文章中位错密度演化包括 edge 位错产生、偶极形成、湮灭和 climb 恢复等机制。孪晶不是凭空出现的。它需要 Shockley partial dislocations 在相邻{111}面上形成一个稳定的三层 stacking fault 结构,然后在外加应力作用下长大。1,层错能 SFE:层错能越低,形成 stacking fault 越容易;2,孪晶核尺度 Ltw:可长大的孪晶核越大,临界应力越低。ε-马氏体相变和孪晶类似,也由 Shockley partials 触发,但多了一个相变能项。文章把 ε-马氏体形核的总能量写成:其中 是奥氏体到 ε-马氏体的 Gibbs 自由能变化项。对应临界相变应力为: 这部分是文章的亮点:TRIP 的发生不仅由应力控制,也由 SFE、界面能和 Gibbs 自由能共同控制。低温:TRIP 更强
在低温下,SFE 降低,cross-slip 更困难,ε-马氏体更容易形成。因此文章模拟和实验都显示:ε-马氏体体积分数随温度降低而增加。(温度越低,奥氏体越“不稳定”,更容易通过 ε-马氏体相变释放变形能,这就是 TRIP 效应增强的原因。)
中温:TWIP 更强
在室温附近或中等温度范围内,孪晶开始发挥重要作用。孪晶既能贡献塑性变形,又会作为位错运动障碍,提高加工硬化。(孪晶像是在晶粒内部不断生成的细小“墙壁”,它们一方面参与变形,另一方面阻碍位错继续滑移,从而提高硬化能力。)
高温:主要是位错滑移
文章中 673 K 和 773 K 的高温曲线主要表现为纯滑移机制,TRIP/TWIP 基本不活跃。原因是高温下 SFE 升高,相变和孪晶的驱动力下降。
作者的模拟效果:
值得注意的是,作者在文章中没有显式描述孪晶/马氏体片层形貌,只用一个唯象的变量去表征滑移和孪晶的体积分数演化,也没有考虑孪晶或者马氏体内部的后续滑移行为.
使用作者完整的本构框架,可以是现在abaqus的umat或者vumat子程序里面,方便后续对高锰钢的复杂变形特征分析研究。
测试了包含200个晶粒的二维拉伸模拟效果:
初始晶粒形貌:
变形后的总位错密度:
变形后的孪生体积分数:
变形后的相变体积分数:
变形后的等效应力云图: