文章名稱《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個歐拉角分布：