金屬凝固元胞自動機建模的案例
基于有限元-元胞自動機法（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?
展開 元胞自動機模擬動態再結晶
有沒有大佬會用元胞自動機模擬動態再結晶晶粒長大的，有償代做微ddw1679
解析DEFORM軟件中的元胞自動機法
元胞自動機法最早提出用于模擬生命系統所具有的自我復制功能，其數學模型是時間、空間、狀態都離散，空間相互作用和時間因果關系為局部的網格動力學模型，能夠模擬復雜系統時空演化過程，廣泛應用于數學、物理學、生物學、化學、地理學和經濟學等各個學科的非線性現象和分形結構的研究。
Hesselbarth和Gobel最早將元胞自動機法應用到再結晶方面，他們的模型研究了再結晶形核和晶核長大的動力學以及其不同的參數和算法對再結晶行為的影響，結果成功得描述了已被公認的再結晶動力學理論JMAK方程。隨著國內外大量研究人員進一步發展完善模型，將元胞自動機法應用于不同金屬材料再結晶過程，與實驗測試得到的再結晶結果吻合。
DEFORM軟件以模擬金屬變形和熱處理過程為主要目的，在不斷深入研究發展中，加入了金屬微觀組織演變模擬，能夠從宏觀和介觀兩個尺度下模擬金屬材料變形行為和組織演變過程，不但具有經典的JMAK法用于金屬再結晶模擬，而且包含了當前流行的元胞自動機法和蒙特卡洛法，能夠直觀的分析觀察晶粒演變過程。
DEFORM中的CA法介紹
目前CA法在再結晶模擬方面的大部分研究與應用，都是針對具體的材料和特定變形條件下，研究人員通過Fortran、MATLAB等編譯軟件編程定義轉變規則和圖形可視化，無法直接輸入實際復雜的工藝加工過程，適用普遍性不強，難以推廣。
DEFORM軟件做為成熟的商業化軟件，使用向導式界面設置界面，流程化操作，簡單易用，元胞中的轉化規則采用位錯模擬模型，與軟件中的JMAK方法可形成對照，互相印證，模擬結果可靠。
▲ CA模型設置界面
工件研究位置的選擇
元胞自動機晶粒組織演變模擬都是在介觀尺度下的，不可能同時對一個實際工件的所有位置模擬計算，否則計算工作量太大無法實現，因此首先需要確認的是分析哪個位置點。
展開 解析DEFORM軟件中的元胞自動機法
元胞自動機法最早提出用于模擬生命系統所具有的自我復制功能，其數學模型是時間、空間、狀態都離散，空間相互作用和時間因果關系為局部的網格動力學模型，能夠模擬復雜系統時空演化過程，廣泛應用于數學、物理學、生物學、化學、地理學和經濟學等各個學科的非線性現象和分形結構的研究。Hesselbarth和Gobel最早將元胞自動機法應用到再結晶方面，他們的模型研究了再結晶形核和晶核長大的動力學以及其不同的參數和算法對再結晶行為的影響，結果成功的描述了已被公認的再結晶動力學理論JMAK方程。隨著國內外大量研究人員進一步發展完善模型，將元胞自動機法應用于不同金屬材料再結晶過程，與實驗測試得到的再結晶結果吻合。
DEFORM軟件以模擬金屬變形和熱處理過程為主要目的，在不斷深入研究發展中，加入了金屬微觀組織演變模擬，能夠從宏觀和介觀兩個尺度下模擬金屬材料變形行為和組織演變過程，不但具有經典的JMAK法用于金屬再結晶模擬，而且包含了當前流行的元胞自動機法和蒙特卡洛法，能夠直觀的分析觀察晶粒演變過程。
DEFORM中的CA法介紹
目前CA法在再結晶模擬方面的大部分研究與應用，都是針對具體的材料和特定變形條件下，研究人員通過Fortran、MATLAB等編譯軟件編程定義轉變規則和圖形可視化，無法直接輸入實際復雜的工藝加工過程，適用普遍性不強，難以推廣。DEFORM軟件做為成熟的商業化軟件，使用向導式界面設置界面，流程化操作，簡單易用，元胞中的轉化規則采用位錯模擬模型，與軟件中的JMAK方法可形成對照，互相印證，模擬結果可靠。
CA模型設置界面
1、工件研究位置的選擇。元胞自動機晶粒組織演變模擬都是在介觀尺度下的，不可能同時對一個實際工件的所有位置模擬計算，否則計算工作量太大無法實現，因此首先需要確認的是分析哪個位置點。
展開 
python自動元胞機方法實現晶粒生長模擬，二維 ￥39
python模擬晶粒生長
晶體塑性耦合元胞自動機模擬熱壓縮過程中的再結晶行為
本構理論分成晶體塑性和再結晶兩部分，其中晶體塑性部分公式如下：
流動方程（經典的唯象流動）：
硬化方程使用的taylor位錯模型
位錯密度的演化使用經典的KM方程：
再結晶部分公式包含形核和晶界遷移兩部分，其中形核的理論公式是
晶界遷移速度為：
整體數值實現框架示意圖如下：
作者以OFHC銅為研究對象，對775K和875K的熱壓縮進行了研究，分析了溫度對再結晶的影響，以及定向形核和生長選擇兩類機制的差異，同時模擬了順序耦合的 DRX→SRX（退火）過程及異常晶粒長大（AGG），模擬效果如下：
根據作者提供的思路（相對簡單清晰），可以編寫對應的子程序，完整晶體塑性和元胞自動機的完全耦合，同樣使用隱式umat實現。數值案例如下：
建立一個包含20個晶粒8000個單元的RVE模型，如下所示
給定對應的初始形核臨界位錯密度和初始的形核率計算公式以及晶界遷移率公式，通過施加周期性邊界PBC沿著X方向壓縮45%（使用鎳基高溫合金的材料參數）。
根據FCC的取向差計算公式，得到初始的晶界分布：
初始的IPF圖如下：
初始的晶粒尺寸分布（mm）：
變形45%后的IPF圖如下：
變形45%后的晶界分布情況：
變形45%后的應力分布情況：
變形45%后的位錯密度分布情況：
變形45%后的晶粒尺寸分布情況：
感興趣的歡迎加入知識星球交流討論，當前效果是初步的建模分析結果：
展開 