織構模擬的案例
abaqus調用damask實現FCC，BCC，HCP多晶織構演化和應力應變場分布模擬
FCC------以鋁為代表，參數使用原始abaqus提供的參數
織構演化模擬模型使用包含1000個單元的1*1*1mm立方體，其中每個單元表示一個特定取向的單晶，初始織構使用軟件生成1000組隨機取向，并分配給不同的單元，模型和初始織構如下圖所示，
利用周期性邊界條件分別模擬多晶沿著ND方向拉伸，壓縮，以及沿著ND方向進行平面應變壓縮時的織構
RD拉伸織構：
RD壓縮織構：
ND平面應變壓縮織構：
BCC------以鐵素體為代表，參數使用原始abaqus提供的參數
織構演化模擬模型使用包含1000個單元的1*1*1mm立方體，其中每個單元表示一個特定取向的單晶，初始織構使用軟件生成1000組隨機取向，并分配給不同的單元，模型和初始織構如下圖所示，
、
利用周期性邊界條件分別模擬多晶沿著ND方向拉伸，壓縮，以及沿著ND方向進行平面應變壓縮時的織構
拉伸織構：
壓縮織構：
平面應變壓縮織構：
HCP------以鎂為代表，參數使用原始abaqus提供的參數
織構演化模擬模型使用包含1000個單元的1*1*1mm立方體，其中每個單元表示一個特定取向的單晶，初始織構使用軟件生成1000組隨機取向，并分配給不同的單元，模型和初始織構如下圖所示，
利用周期性邊界條件分別模擬多晶沿著ND方向拉伸，壓縮，以及沿著ND方向進行平面應變壓縮時的織構
拉伸織構：
壓縮織構：
平面應變壓縮織構：
多晶局部應力應變場分布模擬與宏觀應力應變響應。以FCC－鋁為例子。BCC與HCP同理。
展開 abaqus調用damask實現FCC，BCC，HCP多晶織構演化和應力應變場分布模擬
FCC------以鋁為代表，參數使用原始abaqus提供的參數
織構演化模擬模型使用包含1000個單元的1*1*1mm立方體，其中每個單元表示一個特定取向的單晶，初始織構使用軟件生成1000組隨機取向，并分配給不同的單元，模型和初始織構如下圖所示，
利用周期性邊界條件分別模擬多晶沿著ND方向拉伸，壓縮，以及沿著ND方向進行平面應變壓縮時的織構
RD拉伸織構：
RD壓縮織構：
ND平面應變壓縮織構：
BCC------以鐵素體為代表，參數使用原始abaqus提供的參數
織構演化模擬模型使用包含1000個單元的1*1*1mm立方體，其中每個單元表示一個特定取向的單晶，初始織構使用軟件生成1000組隨機取向，并分配給不同的單元，模型和初始織構如下圖所示，
、
利用周期性邊界條件分別模擬多晶沿著ND方向拉伸，壓縮，以及沿著ND方向進行平面應變壓縮時的織構
拉伸織構：
壓縮織構：
平面應變壓縮織構：
HCP------以鎂為代表，參數使用原始abaqus提供的參數
織構演化模擬模型使用包含1000個單元的1*1*1mm立方體，其中每個單元表示一個特定取向的單晶，初始織構使用軟件生成1000組隨機取向，并分配給不同的單元，模型和初始織構如下圖所示，
利用周期性邊界條件分別模擬多晶沿著ND方向拉伸，壓縮，以及沿著ND方向進行平面應變壓縮時的織構
拉伸織構：
壓縮織構：
平面應變壓縮織構：
多晶局部應力應變場分布模擬與宏觀應力應變響應。以FCC－鋁為例子。BCC與HCP同理。
展開 基于vpsc7.0的FCC不同工況下織構演變模擬------案例十二
? 基于vpsc7.0的FCC不同工況下織構演變模擬
案例實操
1，建立包含1000個晶粒隨機取向的初始晶粒
2，采用Voce硬化模型，獲得材料的拉伸曲線
3，分別采用單向拉伸，壓縮，平面應變變形（100%）
4，后處理，織構演化
材料的初始織構
材料的應力應變曲線
拉伸100%后的取向分布
壓縮100%后的取向分布
軋制下壓100%后的取向分布
軋制過程中織構分數的演化
晶體塑性模擬中特征織構分量離散為由高斯分布組成的等效織構
參考文獻：《Effects of texture on shear band formation in
plane strain tension/compression and bending》
doi：10.1016/j.ijplas.2006.03.014
三類典型織構與三種離散程度(針對FCC結構)
Cube：
Miller indices{0 0 1} <1 0 0>
euler angles（°） 0 0 0
初始極圖：
距離理想織構2°偏差的高斯分布（2000個取向）極圖分布：
距離理想織構5°偏差的高斯分布（2000個取向）極圖分布：
距離理想織構15°偏差的高斯分布（2000個取向）極圖分布：
Brass：
Miller indices{0 1 1} <2 1 1>
euler angles（°） 35 45 0
初始極圖：
距離理想織構2°偏差的高斯分布（2000個取向）極圖分布：
距離理想織構5°偏差的高斯分布（2000個取向）極圖分布：
距離理想織構15°偏差的高斯分布（2000個取向）極圖分布：
Goss：
Miller indices{0 1 1} <1 0 0>
euler angles（°） 0 45 0
初始極圖：
距離理想織構2°偏差的高斯分布（2000個取向）極圖分布：
距離理想織構5°偏差的高斯分布（2000個取向）極圖分布：
距離理想織構15°偏差的高斯分布（2000個取向）極圖分布：
三類典型織構等比例初始取向離散：
初始極圖：
距離理想織構2°偏差的高斯分布（2000個取向）極圖分布：
距離理想織構5°偏差的高斯分布
展開 
基于粘塑性自洽模型（VPSC）FCC結構金屬拉伸壓縮過程中織構的演化模擬------案例二十五
VPSC是由加州大學洛斯阿拉莫斯國家實驗室（大學）的董事會與美國能源部聯合開發全稱VISCO-PLASTIC
SELF-CONSISTENT (VPSC)
目前使用的是更新于2012年5月1日的vpsc7d，
VPSC 是用 FORTRAN 77 編寫的計算機代碼，用于模擬多晶聚集體的塑性變形。 VPSC 全稱 Visco Plastic Self Consistent，指的是特定的機械狀態 (VP) 和使用的方法 (SC)。 VPSC 是為應用于低對稱材料（六邊形、三角形、正交、三角形）而開發的，盡管它在立方材料上也表現良好。
VPSC 說明了單晶和聚集體的特性和響應的完全各向異性。它模擬了骨料在外部應變和應力作用下的塑性變形。
VPSC 基于滑移和孿晶的物理剪切機制，并考慮了晶粒相互作用效應。除了提供宏觀應力-應變響應外，它還解釋了單個晶粒的硬化、重新定向和形狀變化。
因此，它預測了與塑性成形相關的硬化和織構的演變。模擬程序可應用于金屬、金屬間化合物和地質聚集體的變形。
展開 Damask和abaqus晶體塑性聯合仿真培訓通知
（7）Ubuntu上abaqus的簡單使用以及注意事項
五、相關案例
案例一：FCC------以鋁為代表，參數使用原始abaqus提供的參數
織構演化模擬模型使用包含1000個單元的1*1*1mm立方體，其中每個單元表示一個特定取向的單晶，初始織構使用軟件生成1000組隨機取向，并分配給不同的單元，模型和初始織構如下圖所示，
利用周期性邊界條件分別模擬多晶沿著ND方向拉伸，壓縮，以及沿著ND方向進行平面應變壓縮時的織構
RD拉伸織構：
RD壓縮織構：
ND平面應變壓縮織構：
案例二：BCC------以鐵素體為代表，參數使用原始abaqus提供的參數
織構演化模擬模型使用包含1000個單元的1*1*1mm立方體，其中每個單元表示一個特定取向的單晶，初始織構使用軟件生成1000組隨機取向，并分配給不同的單元，模型和初始織構如下圖所示，
利用周期性邊界條件分別模擬多晶沿著ND方向拉伸，壓縮，以及沿著ND方向進行平面應變壓縮時的織構
拉伸織構：
壓縮織構：
平面應變壓縮織構：
案例三：HCP------以鎂為代表，參數使用原始abaqus提供的參數
織構演化模擬模型使用包含1000個單元的1*1*1mm立方體，其中每個單元表示一個特定取向的單晶，初始織構使用軟件生成1000組隨機取向，并分配給不同的單元，模型和初始織構如下圖所示，
利用周期性邊界條件分別模擬多晶沿著ND方向拉伸，壓縮，以及沿著ND方向進行平面應變壓縮時的織構
拉伸織構：
壓縮織構：
平面應變壓縮織構：
案例四：多晶局部應力應變場分布模擬與宏觀應力應變響應。以FCC－鋁為例子。BCC與HCP同理。
展開 Damask和abaqus晶體塑性聯合仿真培訓通知
（7）Ubuntu上abaqus的簡單使用以及注意事項
五、相關案例
案例一：FCC------以鋁為代表，參數使用原始abaqus提供的參數
織構演化模擬模型使用包含1000個單元的1*1*1mm立方體，其中每個單元表示一個特定取向的單晶，初始織構使用軟件生成1000組隨機取向，并分配給不同的單元，模型和初始織構如下圖所示，
利用周期性邊界條件分別模擬多晶沿著ND方向拉伸，壓縮，以及沿著ND方向進行平面應變壓縮時的織構
RD拉伸織構：
RD壓縮織構：
ND平面應變壓縮織構：
案例二：BCC------以鐵素體為代表，參數使用原始abaqus提供的參數
織構演化模擬模型使用包含1000個單元的1*1*1mm立方體，其中每個單元表示一個特定取向的單晶，初始織構使用軟件生成1000組隨機取向，并分配給不同的單元，模型和初始織構如下圖所示，
利用周期性邊界條件分別模擬多晶沿著ND方向拉伸，壓縮，以及沿著ND方向進行平面應變壓縮時的織構
拉伸織構：
壓縮織構：
平面應變壓縮織構：
案例三：HCP------以鎂為代表，參數使用原始abaqus提供的參數
織構演化模擬模型使用包含1000個單元的1*1*1mm立方體，其中每個單元表示一個特定取向的單晶，初始織構使用軟件生成1000組隨機取向，并分配給不同的單元，模型和初始織構如下圖所示，
利用周期性邊界條件分別模擬多晶沿著ND方向拉伸，壓縮，以及沿著ND方向進行平面應變壓縮時的織構
拉伸織構：
壓縮織構：
平面應變壓縮織構：
案例四：多晶局部應力應變場分布模擬與宏觀應力應變響應。以FCC－鋁為例子。
展開 HCP多晶變形與織構演化模擬
模擬通過編程ｕｍａｔ子程序實現，正確性通過ｄａｍａｓｋ程序進行對比驗證。驗證正確性通過織構演化確定。
模擬考慮三組滑移+一組孿晶
滑移為BASAL <a>{0 0 0 1}<1 1 -2 0>,3組
PYRAMIDAL<c+a>{1 1 -2 2}<-1 -1 2 3> 6組 （二階）
PRISMATIC {1 0 -1 0}<1 1 -2 0> 3組
孿晶為TENSILE TWIN {1 0 -1 2}<-1 0 1 1> 6組
１，沿著RD壓縮50%：
編寫子程序結果
ｄａｍａｓｋ模擬結果
２，沿著RD拉伸50%：
編寫子程序結果
ｄａｍａｓｋ模擬結果
３，沿著ND平面應變壓縮結果：
編寫子程序結果
ｄａｍａｓｋ模擬結果
多晶拉伸變形模擬（參數與ｄａｍａｓｋ保持一致）
包含500個晶粒100000個單元的多晶體模擬，沿著X方向施加50％的工程應變
模擬結果如下：
應力分布
孿晶分布
累計剪切應變分布
BASAL累計剪切分布
PYRAMIDAL<c+a>累計剪切分布
PRISMATIC累計剪切分布
TENSILE TWIN累計剪切分布
模擬結果在織構演化方面與ｄａｍａｓｋ程序具有良好的一致性。并且很好二的表現了孿晶的出現，以及對于塑性變形的相對低的貢獻。
展開 基于粘塑性自恰模型（VPSC）的鈦合金拉伸壓縮織構演變模擬
作者：辭殤
關鍵詞：VPSC；鈦合金；拉伸壓縮；織構演變
粘塑性自恰（VPSC）模型，區別與宏觀本構模型，VPSC模型不僅能夠模擬變形過程中材料宏觀力學性能的演化過程，還可以同時模擬材料內部由于變形引起的織構演化過程，實現宏觀與細觀結合，從而使我們更加深刻地理解材料的變形過程。
本文使用VPSC計算HCP金屬鈦合金的單軸拉伸和單軸壓縮變形過程，實現鈦合金拉伸壓縮過程中的應力應變、織構演變以及滑移孿晶變形機制啟動情況的預測，VPSC程序模擬過程如圖1所示。
圖1 VPSC程序模擬過程圖
VPSC模擬的材料初始極圖由程序隨機生成，其極圖如圖2，可見初始狀態表現為隨機織構，極密度最大值為1.4。在經過25%的拉伸以及壓縮變形后，材料內部織構發生明顯變化，表現出織構特征。圖3所示為單軸拉伸后的織構極圖，圖4所示為單軸壓縮后的織構極圖。
圖2 初始隨機織構極圖
圖3 單軸拉伸織構極圖
圖4 單軸壓縮織構極圖
圖5所示為單軸拉伸過程中的應力應變曲線，圖6所示為單軸壓縮過程中的應力應變曲線。可以看到，由于HCP金屬鈦合金的各向異性導致兩種變形模式下材料的流動應力演變過程以及變形過程中織構的演變有很大差異。
圖5 單軸拉伸應力應變曲線 圖6 單軸壓縮應力應變曲線
圖7所示為VPSC預測的單軸拉伸過程中變形機制相對活性。可以看出，柱面滑移的活性急劇下降至最低點并且之后幾乎為零，基面滑移和錐面滑移占據主導地位，二者的活性隨著變形量的增加持續上升至最大值，拉伸孿晶和壓縮孿晶的活性呈現先緩慢上升又緩慢下降的變化趨勢。
圖8所示為VPSC預測的單軸壓縮過程中變形機制相對活性。
展開 基于VPSC模擬FCC金屬等通道轉角擠壓(ECAE)工藝
在之前的推文中我們使用粘塑性自洽多晶體塑性模型（Visco-plasitic Self Consistant，VPSC）計算了面心立方（fcc）、體心立方（bcc）金屬材料變形過程，實現了織構演變的模擬，應力預測等。本文將介紹VPSC模擬FCC金屬等通道轉角擠壓(ECAE)工藝。等通道轉角擠壓是將多晶試樣壓入一個特別設計的模具中以實現大變形量的剪切變形工藝，主要通過變形過程中的近乎純剪切作用,使材料的晶粒得到細化, 從而材料的機械和物理性能得到顯著改善。等通道轉角擠壓是一種有效的制備超細晶材料的方法。
本處粘塑性自洽多晶體塑性模擬的材料初始取向由程序隨機生成，其(100)、(110)和(111)極圖見圖1，可見初始狀態表現為隨機取向，極密度最大值為1.5。變形過程強加100%的剪切應變，步長為0.2，共50步，用4個過程來描述整個等通道轉角擠壓的變形工藝流程，如圖2，在VPSC模擬中，擠出、擠入、模具的流動軸分別為設置為軸1、2、3。
圖1. 初始隨機織構極圖
ECAE通過90o模反復擠壓樣品,在每道工序中，大約100%的剪切應變被施加，其優點是試樣的截面保持不變，這一過程旨在大幅度減小晶粒尺寸，在保證塑性同時提高屈服應力，模擬結果如下：
(a) ECAE1
(b) 90°CW
(c) 90°CW
(d) ECAE2
圖2 等通道轉角擠壓過程織構模擬結果
從模擬結果可以看到，經過等通道轉角擠壓后的FCC金屬產生了明顯的擇優取向-變形織構，其最大強度為5.5。
最后,有VPSC培訓等相關需求歡迎聯系我們.
VPSC培訓
公眾號:320科技工作室
展開 基于粘塑性自洽模型（VPSC）HCP(AZ31B)結構金屬拉伸壓縮過程中織構的演化與應力應變響應模擬
運行分析時輸入包含文件如下
1，vpsc.in（主文件）
2，TENSIN.3（邊界條件）
3，rand1000.tex（初始取向）
4，AZ31b.sx（單晶屬性）
運行時輸出文件
1，PCYS.OUT（屈服面信息）
2，STR_STR.OUT（應力應變信息）
3, TEX_PHn.OUT（變形結束后取向信息）
輸入
1，具有
1944 個方向的基底織構特征取向文件（文件 RAND1944.TEX）
2，考慮IJP文章選擇合理的材料參數，本案例共有3組滑移系（分別是Prismatic，Basal，Pyramidal 〈c + a〉，）+一個拉伸孿晶系統（拉伸孿晶系統）
材料參數如圖所示：
一，拉伸變形模擬（50%Z方向拉伸）（affine方法）
晶粒初始取向分布
變形結束后晶體取向分布
拉伸過程中滑移系開動情況
二，壓縮變形模擬（50%Z方向壓縮）
晶粒初始取向分布
變形結束后晶體取向分布
壓縮過程中滑移系開動情況
三，平面應變壓縮模擬（50%Z方向軋制）
晶粒初始取向分布
變形結束后晶體取向分布
平面應變壓縮過程中滑移系開動情況
變形過程中的等效應力應變響應
展開 
粘塑性自洽模型(VPSC)在復合工藝中的應用
關鍵詞：晶體塑性 VPSC織構模擬 復合工藝
本期將繼續介紹粘塑性自洽模型(VPSC)在金屬變形過程的應用。VPSC適用于各種金屬材料(如鋁合金、鋼材、鎂合金)，各種加載方式(如單向拉伸、單向壓縮、剪切、平面應變、雙向拉伸等)下的宏觀力學性能和微觀結構演化模擬，也可以針對多相金屬(如雙相鋼等)。在結合有限元軟件后，可擴展VPSC模型的模擬范圍，如扭轉、等通道擠壓及壓剪工藝等。本期將VPSC與宏觀有限元結合，以BCC材料作為研究材料，利用有限元獲取了復合工藝下的邊界條件，分別研究了單向壓縮工藝和復合工藝下的織構演化，對比了不同工藝下產生的織構區別。
圖1為建立的有限元模型及VPSC計算過程，有限元采用abaqus軟件構建，施加以壓縮及復合工藝的邊界條件后進行模擬。
(a) 有限元計算 (b) VPSC計算
圖1 模擬過程
圖2為初始材料的取向，可以看到取向呈現明顯的隨機分布。當在壓縮條件下時，材料中逐漸出現取向聚集，在應變為0.5時出現明顯的<100>//X和<111>//X的絲織構，如圖3所示。在復合工藝下，合金中的織構較為復雜，呈現弱的絲織構和明顯的剪切織構，且隨應變的增加，剪切織構越為顯著，壓縮織構明顯減弱，如圖4所示。
圖1 初始材料的織構
(a) 應變為0.3 (b) 應變為0.5
圖3 壓縮工藝下的織構
(a) 應變為0.3 (b) 應變為0.5
圖4 壓縮+剪切工藝下的織構
從圖5中可以看到，不同工藝下的相對滑移激活完全不同，在單相壓縮工藝下，(101)[1-1-1]處于有利激活位置，而復合工藝下的(101)[11-1]處于最大概率的相對激活位置。
展開 晶體塑性有限元仿真入門(4)--織構演變文獻復現
圖11 織構演變模擬常見的邊界條件
織構演變結果
完成Abaqus構建有限元模型所有關鍵步驟后，輸出inp文件并提交Job，查看織構演變結果如下(由于計算資源的限制，僅計算了simple compression和plane strain compression)：
simple compression
plane strain compression
以多晶體中一號節點為例，在塑性變形過程中它的織構演變如下：
1號節點織構取向演變
參考資料
Polycrystalline Plasticity and the Evolution of Crystallographic Texture in FCC Metals
Texture evolution and mechanical behaviour of irradiated face-centred cubic metals
A User-Material Subroutine Incorporating Single Crystal Plasticity in the ABAQUS Finite Element Program
附件
[取向參數提取程序、歐拉角計算程序、織構程序]
網上參考數據測試：
http://muchong.com/bbs/viewthread.php?
展開 密西西比大學多晶塑性模型建模示例
本案例采用該代碼，模擬相同參數情況下FEPX結果與<a href="/major/abaqus">abaqus結果的比較，以及一些多晶變形案例，織構演化等。
模型運行包含
1，umat_xtal 晶體塑性主程序
2，texture文件，用于織構的輸出
3，test文件，用于確定采取的迭代方案，收斂判據，織構輸出頻率，
4，params_xtal文件，一些參量的類型
5，numbers文件，變量類型，公共區
6，fcc文件，晶體的彈塑性屬性，以及滑移系統
需要注意的是：模型使用的是Voce硬化類似于vpsc模型，具體可以參
參考作者發表的相關文獻以及FEPX使用說明書
使用的材料參數如下圖，對應FEPX模型的案例一
利用FEPX得到的結果與開源代碼得到應力應變響應結果如下圖（存在
較小差異，FEPX使用了高階單元，而在abaqus計算時采用了線性四面
體單元，并且取向存在差異）兩者具有良好的一致性
同時比較了使用taylor模型和多晶模型模擬織構演化的差異
初始取向
FCC拉伸織構
taylor模型結果
多晶模型模擬結果
FCC壓縮織構
taylor模型結果
多晶模型模擬結果
FCC軋制織構
taylor模型結果
多晶模型模擬結果
多晶拉伸變形結果如下
50個晶粒5%應力應變分布
500個晶粒10%變形應力應變響應
在相同初始取向下，其變形過程中應力分布情況與damask的計算結果也保持了良好的一致性。因此感興趣的同學可以嘗試使用該開源代碼。鏈接為Category:Mesoscale - EVOCD (msstate.edu)
展開 干貨：基于MEDC模型的鈦合金熱處理相變模擬
DEFORM MEDC（Microstructure Evolution During Cooling）模型可應用于鈦合金的熱處理冷卻過程相變分析，其目的是模擬雙相α和β 在Ti6-4合金冷卻過程中的顯微組織演化，以預測轉變織構。根據變異選擇規則，MEDC模型預測了初生α相晶粒尺寸、初生α、β相和次生α相的體積分數，并根據變量選擇規則從β轉化為次生α相的織構。
在DEFORM中實現了兩個變量選擇規則，即隨機變體選擇和首選變量選擇規則，以最佳匹配初生α相的極圖。織構可以在羅德里格斯空間或歐拉空間中表示。MEDC模型的輸入數據包括MEDC模擬控制數據、材料定義、從β相到初生α和β到次生α的相變數據、初始織構、α和β相的初始體積分數以及初始的初生α相晶粒尺寸。對于DEFORM數據庫中保存的每個步驟，輸出數據包括上述微觀結構特征。此分析技術提供了織構分析工具，如極圖、反極圖和HCP晶體的Kearns數，特別是對于總α相的統計。
該模型可以在完成變形織構模型后運行，也可以使用已知或假定的變形織構作為獨立的變形織構模型運行。在變形紋理模型完成計算后，可以運行MEDC模型，從而將預測的變形織構作為初始織構輸入到MEDC模型中。如果MEDC被用作一個獨立的模型來預測變形織構，而沒有任何先前的變形織構模型結果，那么用戶可以使用一個典型的DEFORM熱處理數據庫以及測量的變形織構(EBSD)或假定的變形織構數據作為MEDC模型的初始織構。
DEFORM能夠進行α+β相及亞穩β相的熱處理冷卻過程（淬火及空冷或爐冷）β向初生α相、β向次生α相的轉變，從而有利于控制熱處理相成分的最終轉變量及相尺寸。
展開 