顯微組織演化模擬的案例
鍛造工藝對Ti-30Zr-5Al-3V合金顯微組織及性能的影響研究
實驗結果及討論
在鈦合金特別是α+β雙相鈦合金中，可以觀察到各式各樣的組織，鈦合金的力學性能在很大程度上取決于這兩個相的比例、形態、尺寸和分布。目前，最典型且應用最廣泛的TC4鈦合金的四種典型組織對應的性能如表3所示。
表3 TC4典型組織的室溫力學性能
由表3可以看出，不同的顯微組織對TC4的室溫性能有很大的影響。隨著α相尺寸的減小，棒材的強度有所降低，延伸率顯著提高。許多研究者認為，復雜合金化鈦合金的強化水平是由加入的各種元素產生的強化疊加起來的，由此，通過計算可得TC4合金中Al強度當量為9.8，Mo強度當量為2.4，計算強度約為940MPa，這與其真實值基本相當。對于Ti-30Zr-5Al-3V合金，按名義成分進行計算，Al強度當量為15，Mo強度當量為1.76，進而計算出該合金的強度值約為1220MPa。相關學者研究發現：鈦鋯合金相對于純鋯和純鈦，合金熔點降低了，耐腐蝕性變強了，并且由于鈦鋯合金具有更細小的顯微組織，其強度和塑性都比純鋯和純鈦高。
通過表2設計的鍛造工藝，獲得了Ti-30Zr-5Al-3V合金的四類典型組織：魏氏組織、網籃組織、雙態組織和等軸組織。通過分析不同組織形貌對應的力學性能，積累Ti-30Zr-5Al-3V合金的性能特點。
不同鍛造工藝對顯微組織的影響
由圖1可以看出：添加了30%Zr的Ti-30Zr-5Al-3V合金屬于典型的α+β雙相鈦合金，其組織形貌與兩相鈦合金的組織類型一致，這與Zr元素是中性元素的結論是一致的。
1#工藝圖1(a）、1(b）為全β區變形組織，因此顯微組織可以清晰地看到原始β晶粒及三叉晶界，經700℃/1h·AC退火處理后晶粒內的α束集寬化，每個晶粒內的α束取向趨于一致且平直并列，有少量編織狀組織。
展開 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累計剪切分布
模擬結果在織構演化方面與ｄａｍａｓｋ程序具有良好的一致性。并且很好二的表現了孿晶的出現，以及對于塑性變形的相對低的貢獻。
展開 激波作用下顆粒層動態演化的雙流體模擬
圖2（c）給出了顆粒層固含率隨時間的演化。初始時刻顆粒層固含率設置為0.36，由于顆粒層的膨脹，局部最大固含率在4ms時刻降低到0.1。為了更好地展示固含率沿沖擊方向的變化，圖2（c）中的右側圖例最大固含率設置為0.1。隨著時間的演化，顆粒層寬度逐漸增大，固含率逐漸降低。這是因為高速氣體的夾帶作用使得固體顆粒層形成沿沖擊方向的運動、膨脹的趨勢。
(2)模擬與實驗結果對比
圖3給出了顆粒床層兩側監測點處壓力隨時間的變化，圖中的黑色曲線為實驗測量得到的壓力信號。為了檢驗計算網格對模擬結果定量準確性的影響，數值模擬中考察了三個網格尺寸：2mm、4mm、8mm，在數據處理上將沖擊波前沿抵達顆粒層的時刻定義為t=0時刻。
圖3表明，模擬得到的壓力信號能夠在定量上與實驗結果較好吻合。如圖3（a）所示，對于位于顆粒層左側的P1監測點，在沖擊波前沿抵達顆粒物料層表面后反彈至P1監測點的時間間隔，實驗測量結果約為2.7ms，數值模擬得到的結果為2.5ms；P1監測點處的壓力隨后急劇增大至一極大值，實驗測量結果為6.5bar，模擬結果為6.8bar；P1監測點處的壓力隨后逐漸降低，實驗測得的壓力值呈現較明顯的波動，數值模擬中因為將顆粒床層做了擬流體處理，得到的壓力信號非常平緩，但壓力值整體上都處于實驗測得壓力數據的波動范圍內，參見網格尺寸為2mm和4mm的模擬結果。
圖3（b）對比了P2處壓力信號的實驗結果和數值模擬結果。
展開 退火銅晶粒生長模型（熱力耦合），用于TSV、TGV填充晶粒演化（相場模擬） ￥99
結合電子背散射衍射（EBSD）實驗與耦合熱–力的多晶相場模擬，揭示電鍍 TXV-Cu 在退火過程中的晶粒演化行為及其對可靠性的影響;基于相場方法的退火晶粒演化模型，將溫度依賴的界面遷移率、界面能及熱膨脹效應納入描述框架，從而在數值模擬中再現 TXV-Cu 的微觀組織演變過程。該模型不僅能夠為實驗觀察提供理論支撐，還可進一步用于預測不同工藝參數下 TXV-Cu 的組織演化規律，為優化工藝與提升器件可靠性提供指導。
3D生物打印模擬血管組織
2018年8月28日，南極熊從外媒獲悉，布里格姆婦女醫院的研究人員開發了一種3D生物打印管狀結構的方法，可以更好地模擬人體內的天然血管和導管。 3D生物打印技術允許微調打印組織的特性，例如層數和運輸營養素的能力。 這些更復雜的組織為受損組織提供了潛在可行的替代品。
“體內的血管不均勻，”該研究的資深作者，BWH醫學系的生物工程師Yu Shrike Zhang博士說。 “這種生物打印方法可以生成復雜的管狀結構，模仿人體系統中的結構，比以前的技術具有更高的保真度。”
為了制作生物3D打印機的“墨水”，研究人員將人體細胞與水凝膠混合，水凝膠是一種由親水聚合物組成的柔性結構。然后，他們優化了水凝膠的化學性質，使人體細胞在整個混合物中增殖。
接下來，他們用這種生物墨水填充3D生物打印機的墨盒。他們還開發了一種定制噴嘴，可以連續打印最多三層的管狀結構。研究人員解釋說：“這些可灌注的管狀結構可以在生物打印管的長度上以規則的間隔從單層連續調整到三層。”
許多疾病損害管狀組織：動脈炎，動脈粥樣硬化和血栓形成損傷血管，而尿路上皮組織可能遭受炎性病變和有害的先天性異常。
研究人員發現，他們可以打印出模仿血管組織和尿路上皮組織的組織。他們將人尿路上皮和膀胱平滑肌細胞與水凝膠混合，形成尿路上皮組織。為了打印血管組織，他們使用人內皮細胞，平滑肌細胞和水凝膠的混合物。
打印管具有不同的尺寸，厚度和性質。 Zhang表示，生物打印組織的結構復雜性對其作為天然組織替代品的可行性至關重要。那是因為天然組織很復雜。例如，血管由多層組成，而多層又由各種細胞類型組成。
該團隊計劃繼續進行臨床前研究，以在測試安全性和有效性之前優化生物墨水成分和3D打印參數。
“我們目前正在進一步優化參數和生物材料，”Zhang說。
展開 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同理。
展開 deform 組織模擬方面的一些論文
借花獻佛，都是別人的資料，和大家分享一下啊
Grain+growth,recrystallization+model.pdf
4Cr5MoSiV1熱作工具鋼熱形變組織模擬.pdf
Simulation of microstructures for Alloy 718 blade forging using 3D FEM simulator.pdf
Finite-element analysis of microstructure evolution.pdf
DEFORM二次開發模擬組織的學習資料及代碼 ￥99
在熱成型工藝中，通過組織預測調控產品性能是一個重要的手段。借助DEFORM體積成形有限元模擬相應的熱加工工藝，獲取基礎的熱力數據，此外通過自帶GRAIN模塊或DEF_SIM自編程組織本構方程實現組織預測，資料為模擬熱壓縮過程中鋁合金組織變化代碼，僅供初學者學習和參考
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同理。
展開 基于有限元-元胞自動機法（CAFE）的增材制造過程組織模擬
<p>關鍵詞：增材制造；有限元，元胞自動機，凝固組織，晶體塑性</p><p class="ql-align-justify">增材制造技術是一種先進的數字化制造技術，其采用熱源熔融離散材料（如粉末），并逐層逐道沉積成3維實體構建。這與傳統減材制造 (切削、磨削等) 和等材制造 (鑄造、鍛壓等) 加工材料方式的本質不同。增材制造過程伴隨著快速的熔化和凝固循環，材料經歷復雜的熱歷程。這導致熔池內部及相鄰層、道之間形成獨特的微觀結構，包括精細的枝晶結構、晶粒尺寸、晶粒取向（織構）以及由微觀偏析引起的潛在析出相。這些凝固組織特征直接決定了制件最終的力學性能（如強度、韌性）和物理性能。因此，精準預測和控制凝固組織演變對于增材制造的工業化應用至關重要。</p><p>有限元-元胞自動機（CAFE）法是一種強大的跨尺度模擬方法，為研究增材制造凝固組織形成提供了有力工具。其采用有限元法或有限體積法建立起制造過程的宏觀熔池模型，模擬激光/電子束等熱源移動產生的瞬態溫度場（包括熔池形狀、溫度梯度G、冷卻速率R）、熱應力及潛在的熔池流動。</p><div contenteditable="false" width="100%" class="ql-align-justify">
<img src="https://p3-sign.toutiaoimg.com/tos-cn-i-axegupay5k/6d18f544077e4f7891aafa2bda90eca2~tplv-tt-origin-web:gif.jpeg?
展開 基于粘塑性自洽模型（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 基于滑移和孿晶的物理剪切機制，并考慮了晶粒相互作用效應。除了提供宏觀應力-應變響應外，它還解釋了單個晶粒的硬化、重新定向和形狀變化。
因此，它預測了與塑性成形相關的硬化和織構的演變。模擬程序可應用于金屬、金屬間化合物和地質聚集體的變形。
展開 
精沖鋼微觀組織對其力學性能和精沖性能影響的多尺度模擬研究
結束語
通過對精沖用碳鋼不同微觀組織的拉伸、剪切、精沖性能進行宏微觀模擬，模擬結果顯示：球化退火后，當基體中存在較小體積分數的小直徑碳化物顆粒時，材料的拉伸、剪切強度提高，裂紋萌生時間延遲；而碳化物帶的存在使得帶中產生應力集中，加快損傷演化而導致材料過早斷裂，并且碳化物帶中滲碳體顆粒越多或越密集，這種應力集中越嚴重；精沖時裂紋沿著碳化物帶擴展，降低斷面質量。模擬預測的單元失效情況與掃描電鏡觀察的精沖試樣裂紋擴展結果一致。
作者簡介 莊新村，博士，副研究員，碩士研究生導師，現任上海交通大學材料學院塑性成形技術與裝備研究院副院長，主要研究方向為精沖成形工藝及模具設計優化、板料鍛造新工藝和材料韌性損傷演化建模。
《精沖鋼微觀組織對其力學性能和精沖性能影響的多尺度模擬研究》摘自《鍛造與沖壓》雜志第20期
展開 基于粘塑性自洽模型（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方向軋制）
晶粒初始取向分布
變形結束后晶體取向分布
平面應變壓縮過程中滑移系開動情況
變形過程中的等效應力應變響應
展開 基于晶體塑性有限元方法模擬不同取向單晶鋁簡單拉伸過程中的響應情況和取向演化情況------案例二十四
案例說明
1，建立柱狀單晶鋁模型（直徑10um，高度25um）如下：
2，賦予單晶鋁對應的的單晶材料材料參數，（本案例主要考慮在立方金屬軋板中常見的典型取向）見下表（研究選取了前七種情況+taylor取向）
典型取向
3，進行網格劃分，采用C3D8R單元，共包含網格為5004個單元,網格模型如下：
4，X0面所有自由度均為0，X1面施加X正方向20%工程應變的拉伸位移邊界條件
5，后處理與結果展示（默認圖片中單晶取向與表順序相同）
不同取向單晶拉伸的應力分布云圖
不同取向單晶的累計塑性應變分布云圖
不同取向單晶拉伸過程中應力應變響應
同取向單晶拉伸過程中取向演化（紅色為初始取向，藍色為變形后的取向）
立方取向單晶織構演化
旋轉立方取向單晶織構演化
銅型取向單晶織構演化
黃銅取向單晶織構演化
戈斯取向單晶織構演化
S取向單晶織構演化
R取向單晶織構演化
Taylor取向單晶織構演化
展開 