文章名稱：《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的隨機值）：

變形結束是否發生孿生變形：