镁合金高速拉压为何不同？从滑移与孪生构建动态本构
DOI：10.1016/j.ijmecsci.2016.10.012
推荐理由
同样提高加载速率，镁合金的拉伸和压缩为何反应不同？Xie等以挤压AZ31B为对象，把速率修正分配给具体变形机制，提供了一条比单纯修正宏观曲线更便于讨论物理含义的建模路线。
作者研究思路
室温下沿垂直挤压方向，分别开展高速拉伸和压缩试验。模型区分孪生与滑移，分别计算母晶、孪晶区滑移并按体积分数加权；继承孪生分阶段阻力图景，借鉴PTW速率控制思想修正滑移支路，再用β规则连接单晶和多晶。
图 4：原文沿用Wu等[27]的孪生机制示意。本文新增重点在动态速率描述，图中没有显式求解空间孪晶前沿。
核心公式
式（1—2）：孪生剪切率由特征剪切gtwin与孪晶分数速率相连，再通过Schmid张量汇总为塑性应变率。
式（7）：孪生阈值先降低，再随f接近ffinal而增强。表2设τ0=80MPa、τg=20MPa、fg=0.06、ffinal=0.18；这些是支路参数，不是实测总孪晶分数。
式（12）：以误差函数erf的对数速率因子修正滑移阻力；r为系数，ε̇Lᵍ为晶粒加载方向塑性应变率，m=2为幂指数。速率修正主要放在滑移支路，不宜据此断言孪生普遍率无关。
式（14—16）：β变量参与晶间协调；局部应力与宏观应力相联系，宏观塑性应变按晶粒体积分数平均。这是β规则，并非EVPSC自洽计算。
核心结果图
图 2：1250、1650、2150s⁻¹拉伸，增速后流动应力明显提高。
图 3：1648—3544s⁻¹压缩，曲线差异较弱，不能扩大成压缩在所有条件下都没有速率效应。
图 5节选：1250和2150s⁻¹下，红线较好描述主要流动段，未跟随实验末段下降；这些动态曲线参与了参数标定。
图 6节选：1648和3544s⁻¹压缩，整体硬化趋势得到描述，但实验屈服峰未重现，低速端后段应力也有偏差。
图 8：作者另用0.01s⁻¹的MTS拉压试验比较预测，提供跨速率对照。模型未含失效判据，输出到0.15是人为终点，不是断裂应变预测。
作者核心结论
将滑移与孪生区分处理，可描述本组AZ31B的拉压速率差异。动态曲线拟合与低速对照支持这一建模路线，但压缩屈服峰仍是明确不足。
对我们的启示
按机制设置速率依赖，还应同时检查额外速率、瞬态峰值和后段偏差。孪生弱速率敏感、忽略温升均属模型近似；仅凭宏观曲线吻合，尚不足以独证各机制的贡献。
来源：Xie等，International Journal of Mechanical Sciences 119（2016），107—113。公式、图像为原文节选与重排；图4示意源自Wu等（2015），版权归原权利人。
