1、微觀至介觀尺度的模擬微觀至介觀尺度的模擬微觀至介觀尺度的模擬微觀至介觀尺度的模擬主要研究內容： 微結構演化 （動力學控制） 微結構與其性質之間關系 結構演化的方向熱力學控制 微結構變化路徑動力學控制 結構演化的這種非平衡特性導致了各種各樣的晶格缺陷結構及其相互作用機制。尺度/m特性、現象或缺陷10-10 10-7點缺陷，原子團簇，短程有序，在玻璃態和界面中的結構單元，位錯芯，裂紋尖端，原子核10-9 10-5失穩分解，涂層，薄膜，表面腐蝕10-9 10-4二嵌段共聚物，三嵌段共聚物，星形共聚物，大質量的非熱變化，界面網格，位錯源，堆積效應10-9 10-3粒子、沉積物，枝晶，共晶，共析10-8

2、 10-5微裂紋，裂紋，粉末，磁疇，內應力10-8 10-4堆垛層錯，微帶，微孿晶，位錯通道10-8 10-3聚合物中的球晶，存在于金屬、陶瓷、玻璃及聚合物中的結構疇或晶粒團簇（對于多晶或非晶的情況）10-8 10-2聚合物中的構象缺陷團簇10-8 10-1位錯，位錯壁，旋錯，磁壁，亞晶粒，大角晶界，界面10-7 10-1晶粒，剪切帶，復合材料的第二相10-7 100擴散，對流，熱傳遞，電流傳輸10-7 100微結構逾滲路徑(斷裂，再結晶，界面潤濕，擴散，腐蝕，電流，布洛赫壁)10-6 101表面，樣品斷面收縮，斷面微結構的實物空間和時間尺度微觀至介觀尺度的模擬微觀至介觀尺度的模擬l非平衡因素

3、 材料性質的多樣性應用性質材料微結構機制微觀至介觀尺度的模擬微觀至介觀尺度的模擬介觀尺度模擬的特點： 處理的原子數目巨大(1023個/cm3)。 排除了（1）嚴格求解薛定諤方程（2）由唯象原子論方法(如與經驗勢相聯系的分子動力學)來完成。 必須建立能覆蓋較寬尺度范圍的恰當的介觀尺度模擬方法，以便給出遠遠超過原子尺度的預測。微觀至介觀尺度的模擬微觀至介觀尺度的模擬l連續體模型 原子運動方程的嚴格解或近似解 （薛定諤方程或分子動力學） 替換為平均本征結構關系式 介觀尺度機理和本構定律的復雜性和多樣性，導致建立介觀尺度模型的方法的不唯一性。空間及時間離散化介觀尺度模擬方法空間及時間離散化位錯動力學(

4、晶體塑性，復原，織構，斷裂)相場動力學或廣義Ginzburg-Landau模型(超導電性，擴散，相變，晶粒生長)確定性或概率性元胞自動機(擴散，熱傳遞，相變，再結晶，晶粒生長)多態動力學波茨(Potts)模型(相變，再結晶，晶粒生長)幾何拓撲和組分模型(相變，再結晶，晶粒生長)拓撲網格和頂點模型(晶界動力學，網格動力學，成核，復原，晶粒生長)典型應用領域中的主要介觀尺度模擬方法微觀至介觀尺度的模擬微觀至介觀尺度的模擬l介觀模擬方法的共同特點：不明顯地包含原子尺度動力學，而是理想化地把材料作為連續體。l由均勻性基體將晶格缺陷之間的相互作用耦合在一起。l控制方程中通常不顯含內秉空間或時間標度。l含

5、有單個晶格缺陷的連續體介觀尺度模型通常由一系列唯象的偏微分速率和本征結構方程組表述。采用有限差分法、有限元法或蒙特卡羅方法可以對這些微分方程近行求解。微觀至介觀尺度的模擬微觀至介觀尺度的模擬l時空標度參數和離散度的確定由微分方程及其系數、變量所擁有的特點和性質決定。l作為態變量(例如原子濃度，位鍺密度，結構參數，位移或品格取向)，通常被并進空間格柵坐標；l控制微分方程被用于局域或整體情況，這取決于相互作用的性質(短程或長程)。l能夠利用連續體近似方法對介觀尺度的結構演化進行預測、意義重大，因為唯象態方程和結構演化定律已在介觀尺度進行很好地研究，其實驗數據的獲得比在微觀尺度更容易，而且數據信息比

6、在宏觀尺度更詳細。第第6章章 元胞自動機元胞自動機l6.1基本原理l元胞自動機是描述和處理復雜系統在離散空間-時間上演化規律的算法，通常采用對晶格格座的局域或整體的確定性和概率性變換規則進行具體操作。l空間變量可以代表實空間、動量空間或波矢空間。l其晶格定義為具有固定數目的點，一般是規則晶格，但其維數及大小可以是任意的。它表述了系統由基礎實體形成的構象。l這些“基礎實體” 可以是任意大小的連續體型體積單元、原子顆粒、晶格缺陷或生物界中的動物等等。元胞自動機的原理應用于城市規劃 6.1 6.1 基本原理基本原理l基本實體，由廣義態變量(諸如無量綱數、粒子密度、晶格缺陷密度、粒子速度、顏色、血壓或

7、動物種類等)進行量化表述。l在每一個獨立的格座，這些態變量的實際取值都是確定的。并且認為，每一個結點代表有限個可能的離散狀態中的一個態。l通過將某些變換規則應用于每個結點狀態，就會發生自動機的演化。這些規則決定著晶格格座的狀態；對于局域規則，格座狀態是其前一狀態及近鄰格點(座)狀態的函數，而在整體變換規則下，則為所有格座狀態的函數。傳統元胞自動機大多采用局域變換規則。6.16.1基本原理基本原理l對于在規則晶格結構方面的應用是比較容易。l對于非均勻介質，在討論的晶格區域采用較小的晶格間距比較妥當；而且，還必須考慮對變換速率進行合理修正和重正化。l元胞自動機以離散時間步發展演化。經過一個時間間隔

8、，要對所有結點的態變量值同時更新。l廣義微結構元胞自動機可以采用元胞或格座的離散空間格柵，在空間上通常被認為是均勻的，所有格座都是等價的，并被排布在規則晶格上，其中的變換規則在各處都是一樣的。假定它們是有限個可能狀態中的一個，并對所有元胞狀態同步更新。此外，它們與常規自動機不同的是，格座變換既可以按照確定性定律，也可以按照概率性定律。6.16.1基本原理基本原理l元胞自動機方法為模擬動力學系統的演化提供了一種直接的手段，這些動力學系統包含有大量基于短程相互作用或長程相互作用的相似組元。 l 對于簡單的物理系統，時間是其惟一個獨立變量(自變量)。這種直接方法，就相當于利用有限差分近似法給出偏微分

9、方程組的離散解。l元胞自動機方法對“基礎實體”類型和選用的變換規則沒有任何限制。它們可以描述：簡單有限差分模擬中態變量值的分布，混合算法的色問題，“教室里的兒童健康情況”，在任何變換條件下的模糊集合元素，以及元胞的初級生長與衰減過程等。6.16.1基本原理基本原理l例如，用于計算高次多項式系數或裴波那契數的帕斯卡三角形，可以作為一維元胞自動機。其中規則三角晶格各個格座對應的值，可通過在其上方的兩個數之和給出。在這種情況下、自動機的“基礎實體”是一些無量綱的整數，其變換定律是求和法則。一維元胞自動機模型帕斯卡三角形示意圖6.16.1基本原理基本原理l另一種自動機是由立方晶格組成的，這時每個點具有

10、一種顏色，并能按照下述簡單的變換規則進行轉換：“如果某點有超過50的近鄰格點(座)是藍色，則該點就由原色變成紅色”；或者“當有超過75的近鄰格點是紅色時，那么所考察格點的顏色也轉換為紅色”。若要描述學校里孩子們之間的相互傳染問題，我們可以通過一個規則，亦即“如果一個教室里有50的孩子得病則該教室里其他所有孩子就被感染”，定義一個元胞自動機。為了使上述簡單唯象模型變得更加合理、真實、可信，應該增加更多的變換規則。上面的例子可補充這樣的規則：“經過一定數目的時間步之后，受感染的孩子已康復”或“每個孩子只能被感染一次等等。6.16.1基本原理基本原理l元胞自動機并不簡單地等同于普通模擬方法，例如各種

11、有限差分法、有限元法、伊辛(Ising)法、波茨(Potts)方法等。l元胞自動機具有廣泛的適用性和多功能的特點，是離散計算方法的普遍化推廣。l這種靈活適用性是基于這樣一個事實：除了采用簡明的數學表達式作為變量和變換規則之外，如果需要的話，自動機能夠包括任何元素或規則。6.16.1基本原理基本原理l在材料科學中，有時對常規有限差分計算方法補充一些“如果就”規則可以為處理“數學上的奇點(即非光滑函數表述中的臨界或自發效應)問題提供了一種簡單有效的途徑。事實上，這些規則經常出現在微結構模擬中。例如， 在離散位錯動力學模擬中，“如果兩個反平行螺位錯相互靠近到其間距小于5個伯格斯矢量時，它們就會自發湮

12、沒”；在斷裂力學或彈簧模型中，會經常包含這樣的規則：“如果裂紋速度達到某一個值，試驗樣品將自發損壞”；在重結晶模擬中，會經常遇到這樣的規則：“如果晶體局城取向誤差達到某一個值，格座將滿足成核的動力學非穩定性臨界條件。”或“局域儲存的彈性能達到某個臨界值，格座將滿足成核的熱力學非穩定性條件”。6.16.1基本原理基本原理l如果對主微分方程補充上述所說的“如果就”的變換規則，我們就可以對復雜系統的動力學行為特性進行模擬。通常而言，所考察粒子之間的局域相互作用是這一問題的根本基礎。l盡管元胞自動機模擬一般是在基本尺度層次(例如原子、原子團簇、位錯段、亞晶粒)上完成的，但是作為對連續體空間進行離散化和

13、映射處理的派生方法，本身不存在物理特征線度或時間刻度的內秉標定問題。l對連續體系統的元胞自動機模擬，需要定義相應的基本單元和對應的變換規則，以便恰當地展現系統在給定層次上的行為特性。6.16.1基本原理基本原理l從物理角度看，分子動力學表示的是真正的微觀模型，而在使用元胞自動機方法時，并不局限任何特定體系，可適用于任何系統。與蒙特卡羅方法相比，由元胞自動機方法得到的平衡系綜的熱力學量，在物理上更缺少依據和基礎。由于這個原因，在進行元胞自動機計算機實驗之前，一個重要工作就是，檢驗基本模擬單元是否切實體現了“基礎物理實體”的特性。由于元胞自動機的應用并不局限于微觀體系，所以它為在微結構模擬中實現不

14、同空間及時間尺度的方法之間的跨越，提供了一個非常方便的數值工具。6.2 CA在材料中的多面性在材料中的多面性l由于在考慮大量可能的空間態變量及變換規則時所展現的廣泛適用性和靈活性。元胞自動機方法在對由再結晶、晶粒生長及相變現象等形成的微結構進行模擬時，表現出特有的多面性。l例如，對于再結晶和晶粒生長，元胞自動機可以離散化方式同時描述局域結晶結構及其形成過程。為了實現對這些特性的描述，一般是將局域晶體取向g、儲存的彈性能(即某種近似可測量，諸如位錯密度或局域泰勒因子M)以及溫度T作為態變量。這些變量都是因變量，也就是它們依賴于自變量，諸如空間坐標(x1，x2，x3)和時間t等。6.2 CA在材料

15、中的多面性在材料中的多面性l就特定的研究對象，狀態參量應包含在所使用的各種局域結構演化定律之中。根據局域的信息、數據及變換規律，可以對諸如復原、成核及其生長等現象的機理結出相應合理的唯象解釋。l通過確定二維或三維空間格柵所對應的態變量，元胞自動機模擬可以應用于對微結構的非均勻性質的研究，其中包括諸如第二相、微帶、剪切帶、過渡帶、異相界面、晶界和孿晶等。這些局域性缺陷結構，可以借助其態變量的相應值或梯度值進行表述；用高位錯密度和大的局域晶格曲率表征剪切帶的特性。對于一個給定雜質含量的晶界遷移率m，可以采用相鄰晶粒之間的取向偏差g和晶界法線的空間取向n來表征。6.3 6.3 元胞自動機的一般表述元

16、胞自動機的一般表述l在元胞自動機中，鄰接格座的局域相互作用，是通過一套確定性或概率件變換規則具體確定的。在時間(t+t)時，對應于某特定格座的態變量值將由目前狀態(t0) (或最接近的幾個態t0，t0-t等)及其鄰近格點的狀態決定。若只考慮最鄰近的兩個時間步，則對于一維元胞自動機的演化來說，可以用公式寫成下式形式：l 表示在時間t0時對應于結點j的態變量值；j+1和j-1表示格點j的兩個最近鄰結點。f具體指定了描述變換規則的函數。00000001111,tjtjtjttjttjttjttjf(6.1)0tj6.3 6.3 元胞自動機的一般表述元胞自動機的一般表述l幾種鄰接狀態l馮諾伊曼鄰接 結

17、點狀態僅取決于最鄰近結點l摩爾鄰接 結點狀態取決于最鄰近結點和次鄰近結點l擴展摩爾鄰接 考慮兩層鄰近的元胞l馬哥勒斯鄰接 每次考慮一個22的元胞塊l鄰接類型影響系統的轉換速率和演化形態。 1,11,21,32,12,22,33,13,23,31,11,21,32,12,22,33,13,23,3圖6.1 馮諾依曼鄰接和摩爾鄰接6.3 6.3 元胞自動機的一般表述元胞自動機的一般表述l對于擴展配置，一維情況下，考慮兩個鄰近時間步時的轉換規則可以寫為：l其中n表示單位晶格元胞變換規則的作用范圍。00000000000,111111tnjtnjtjtjtjttjttjttjttnjttnjttjf

18、(6.2)6.3 6.3 元胞自動機的一般表述元胞自動機的一般表述l元胞自動機存在眾多可行的變換規則。 l馮諾伊曼鄰接的一維二進制元胞自動機，每一個時間步，即j=0或 j=1，轉換規則采取 的形式，其轉換規則有28個。其中之一 000011,tjtjtjttjf0000000000000000000000001111111111111,1,101,1,011,0,101,0,010,1,110,1,0tttttjjjjtttttjjjjtttttjjjjtttttjjjjtttttjjjjttttjjjj00000000111100,0,110,0,00ttttttjjjjtttttjjjj該

19、轉換規則可以以(01011010)2的編碼形式表示。在元胞自動機方法中，一般采用數字編碼方式簡化表述相關變換規則。 6.3 6.3 元胞自動機的一般表述元胞自動機的一般表述l變換規則的數目可以由k (kn)計算得到，其中k為元胞的狀態數，n為包含芯元胞在內的鄰近元胞的數目。對于具有摩爾鄰接的二維元胞自動機（n=9），假設每個元胞具有兩個可能的狀態，則該系統將具有229=262144個不同的轉變規則。6.4 6.4 元胞自動機的分類元胞自動機的分類l沃爾弗拉姆（Wolfram）基于動力學行為的差異將元胞自動機分為四類。l(1)平穩型： 自任何初始狀態開始，經過一定時間運行后，元胞空間趨于一個空間

20、平穩的獨一無二的構形，這里空間平穩即指每一個元胞處于固定狀態，不隨時間變化而變化。l(2)周期型： 產生周期性重復的短周期結構，或者產生穩定結構。在這種元胞自動機中同時呈現出局部和整體的排列次序。這種自動機可以看作是一個濾波器（Filter），這來自于給定轉換規則的離散數據的本質。在向空間中這種系統形成閉環。6.4 6.4 元胞自動機的分類元胞自動機的分類l(3)混沌型： 自任何初始狀態開始，經過一定時間運行后，元胞自動機形成非周期的混沌結構。至少在經過一定的時間后，這種結構的統計特征與初始結構的統計特征大致相同。由第三種元胞自動機生成的結構通常為自相似的分形排列。對于任意的初始配置，經過大量

21、的時間步后，這些結構具有相同的統計特征。人們對這類自動機在幾何方面的應用具有很大的興趣。這種自動機是最常用的一種元胞自動機。6.4 6.4 元胞自動機的分類元胞自動機的分類l(4)復雜型： 第四類元胞自動機產生穩定的，周期性的，可以維持任意長時間的傳播結構。一些元胞自動機在經歷一定時間步以后衰退，即所有元胞的狀態變為零。一些第四類元胞自動機可以形成穩定的周期性結構。通過對這些恰當的傳播結構進行設置，可以得到具有任意循環長度的最終狀態。在演化過程中，第四類元胞自動機表現出高度的不可逆性。這種元胞自動機可以呈現出重要的局域排列。6.5 6.5 概率性元胞自動機方法概率性元胞自動機方法 l將確定性元

22、胞自動機變為非確定性的基本方法有兩種:l第一種方法就是隨機地選擇所研究的晶格格點，而不是系統化地按順序選擇，但是要使用確定性變換規則；l第二種方法就是用概率性變換代替確定性變換，但要系統地研究所有格點。l主要討論第二種方法，并將之歸為概率性或隨機性元胞自動機。 6.5 6.5 概率性元胞自動機方法概率性元胞自動機方法l概率性元胞自動機，就其基本過程和要素方面而言，非常相似于普通的元胞自動機，只不過轉變規則由確定性的換成了隨機性的。l設有N個格點組成一個一維鏈，其中每個格點有k個可能的狀態Sv=0, 1, 2, ., k-1。從而整個鏈共有kN個不同的排列方式。由（S1, S2, , SN）描述

23、的某給定晶格狀態用下式整數標記：(6.3)NvvvkSi116.5 6.5 概率性元胞自動機方法概率性元胞自動機方法l假設每個狀態i的存在概率為Pi。作為時間的函數的Pi(t)，按照其轉變概率以離散時間步t=0, 1, 2, 的方式變換發展。如果只考慮鄰近的時間步(t-1)，這一規則可用下式給出：l轉移概率Tij就表示由前一時刻的j狀態轉變到i狀態的概率。l對于離散型元胞自動機方法，轉移矩陣Tij 是由局部規則決定：(6.4)10) 1()(NkjjijitPTtPNvjvjvjvivijSSSSpT111),|(6.5)6.5 6.5 概率性元胞自動機方法概率性元胞自動機方法l 和 分別表示

24、狀態j和i的格點變量；因而，變量的轉換只有其最近鄰及其自己的狀態有關。 l雖然概率性元胞自動機與Metropolis蒙特卡洛算法之間具有一定的相似性，但二者之間還是有差別的。這種差別主要表現在兩個方面：l第一，蒙特卡洛方法每個時間步只更新一個格點，而概率元胞自動機像大多數自動機一樣，每次要全部一起更新；l第二，元胞自動機沒有本征的長度或時間標度。l元胞自動機的標定參數主要是由構成物理模型的基礎來決定，而不是由所采用的元胞自動機算法來決定。jvSivS6.5 6.5 概率性元胞自動機方法概率性元胞自動機方法l元胞自動機與波茨蒙特卡羅自旋模型的區別：l(1)元胞自動機對微觀體系不存在內稟標度，如果

25、選擇合適的基礎單元，并且建立與場變量相匹配的代數、微分或積分方程，那么元胞自動機可以用于任意空間和時間尺度上的問題處理。蒙特卡羅方法對于微觀體系是有內稟標度的。l(2)在蒙特卡羅方法中，廣義自旋格座是隨機抽樣順序考察的，而元胞自動機則是同步一起更新。l(3)元胞自動機比多態波茨模型使用了更多的確定性或概率性變換規則。6.6 6.6 非平衡現象的模擬非平衡現象的模擬l6.6.1熱力學模擬 l在金屬的熱變形過程中，會發生諸如再結晶、連續與非連續型晶粒生長和不連續沉淀等非平衡轉變現象和微結構瞬態問題。按照微結構的觀點，這些轉變現象都是由于高角晶界的運動引起的。由于吉布斯自由焓存在梯度，原子或原子團將

26、從一個晶粒越遷轉移到其鄰近晶粒。對同相界面，其凈驅動壓強(6.6)dVdGp 6.6.16.6.1熱力學模擬熱力學模擬l在實際材料中，各種貢獻都將影響到局域自由焓的值。l(1)在冷加工金屬中，位錯密度的增加對所儲存的彈性能的貢獻，在對驅動壓強中占最大的份額。l是界面兩邊的位錯密度差；為各向同性極限下的體剪切模量；b表示伯格矢量的大小。l如果將存在于元胞壁的位錯(w)和元胞內的位錯(i)的貢獻分別表述。(6.7) 221bp(6.8)Dbpsubiwi221,lw只能用亞晶粒尺寸D和亞晶粒壁的界面能sub表述6.6.16.6.1熱力學模擬熱力學模擬l(2)作用于各晶粒上的拉普拉斯壓強或毛細壓強的

27、貢獻。l對于常見的晶粒粒度分布和球形晶粒l為23的常數；為界面能；1/R為曲率。l對于薄膜，還有來自表面能梯度的貢獻：l式中，B表示薄膜寬度；h為膜厚；代表表面能變化量。(6.9) Rp(6.10) hBhBdxBdxps22216.6.16.6.1熱力學模擬熱力學模擬l(3)在過飽和態，對驅動壓強還有一項化學貢獻。其對應的轉變稱為非連續沉淀。lkB為波耳茲曼常數；l為原子體積；lT1為（數值）實驗中的實際溫度；lT0為相應于T1時過飽和濃度的平衡溫度；c0為濃度。 (6.11) 0010lnccTTkcpB6.6.16.6.1熱力學模擬熱力學模擬l(4)其他因素對總壓強的貢獻l在總的驅動壓強

28、中，還要考慮冷加工或硬化金屬間化合物中由于損失長程有序而產生的貢獻。更進一步還應當考慮來自于磁性、彈性及溫度場等梯度的貢獻，但是這類貢獻在實際應用中意義不大。l(5)可能的反驅動力l雜質阻力以及在有序化合金中高角晶界運動在遠處產生疇的結構。6.6.2 6.6.2 動力學模擬動力學模擬l為使原級再結晶能夠啟動，要在熱力學、力學和動力學方面有一定的不穩性。l第一類不穩定性就是成核，l第二類就是有凈驅力，l第三類就是高角晶界的運動。l在再結晶過程中，主要是非均勻成核。l可能成核的格點所處的區域應該具有非常高的位錯密度和較小的子晶粒尺寸，以及具有較大的局域晶格取向偏差。l例如：剪切帶、微帶、遷移帶、存

29、在高角晶界、在沉淀周圍的形變區等。 6.6.2 6.6.2 動力學模擬動力學模擬l采用垂直通過均勻晶界的各向同性單原子擴散過程，則用于描述界面運動的對稱速率方程可以寫為：l 表示界面速度；vD是德拜頻率；gb表示通過界面時的跳變寬度；c表示平面內自擴散截體缺陷的固有濃度（如晶界空位或源的重組）；n表示晶界片的法向矢量；Gt是與轉變有關的吉布斯焓; kB為波耳茲曼常數;T為絕對溫度。TkGGTkGGncvxBtBtgbD2exp2exp(6.12)x6.6.2 6.6.2 動力學模擬動力學模擬l將焓、熵及驅動壓強帶入式（6.12），則有：lp為驅動力（如儲存的彈性能或界面曲率）；為原子體積；Sf

30、表示形成熵；Hf表示形成焓；Sm表示運動熵；Hm表示運動焓。 (6.13)TkpSTHTkpSTHTkHkSnvxBmmBmmBfBfgbD2exp2expexpexp6.6.2 6.6.2 動力學模擬動力學模擬lSf主要是振動熵，而Sm包含有組態和振動兩者的貢獻，則式（6.13）變為:l考慮到雙曲函數中的 是個小量(6.14)TkHHTkpkSSbnvxBmfBBmfDexpsinhexpTkpB(6.15)TkHHTkpkSSbnvxBmfBBmfDexpexp6.6.2 6.6.2 動力學模擬動力學模擬l晶界遷移率實驗數據阿倫烏斯分析的著名唯象表達式：l式中，m表示遷移率，Qgb表示晶界

31、運動的激活能。比較式（6.15）和式（6.16）兩式中的系數，則有：l晶界運動的經典動力學圖像。 pTkQnmpmnxBgbexp0(6.16)(6.17)mfgbBmfBHHQkSSTkbvmexp006.6.2 6.6.2 動力學模擬動力學模擬l在原級再結晶的初級階段，局域復原過程促進了晶核的形成。在其最后階段，位錯湮滅及重新排列將引起所儲存能量的不斷降低，從而使局域驅動力明顯減小，最終導致再結晶速度減慢。l(t)是作為時間t的函數的位錯密度，0表示形變后的位錯密度，為弛豫時間。 (6.18) ttexp06.6.3 6.6.3 確定性元胞自動機解法確定性元胞自動機解法l冷加工金屬中原級再

32、結晶模擬的確定性元胞自動機方法。l假定成核和新結晶晶粒長大所需驅動力均來源于局域位錯密度的梯度；并且當有碰撞時生長終止。起始數據應包括格柵幾何參數和態變量取值等信息，例如，溫度、成核概率、晶界遷移率、位錯密度和晶體取向。這些數據必須能夠描述作為空間函數的初始微結構的主要特征。 l原級再結晶的物理過程：復原、成核和晶核生長。6.6.3 6.6.3 確定性元胞自動機解法確定性元胞自動機解法l(1)復原階段l位錯密度與驅動力相聯系并對成核速度有潛在的影響。l在簡單的有限差分公式中，因子f與弛豫時間、溫度T和時間ti有關。 l(2)在成核階段l各元胞或元胞團簇由變形態轉變為再結晶狀態。(6.19) 2

33、1032121321,TtxxxtTfTtxxxii6.6.3 6.6.3 確定性元胞自動機解法確定性元胞自動機解法l根據所構造模型的物理基礎，相對周圍的形變基體，其晶核取向可有三種情況：l即相同取向l相似取向（g15o）l不同取向（g15o）l相同或相似取向的晶核只能存在于高角晶界，而當取向偏差明顯不同時晶粒會生長進入近鄰晶粒。 l把這些反映晶粒取向特性的臨界條件，以及產生的晶核應與基體有相似取向的規則結合起來，就相當于給出一個取向成核的假說。 6.6.3 6.6.3 確定性元胞自動機解法確定性元胞自動機解法l(3)生長階段l在生長階段，對于每個晶粒可執行一個循環，這個循環遍及所有屬于目標晶

34、粒表面的元胞，可以確定表面元胞與其非再結晶近鄰元胞兩者結晶取向偏差g和溫度T的函數。l在原級晶界的情況下，局域驅動力取決于非再結晶元胞的實際位錯密度。驅動力和遷移率決定著晶界運動的速度；晶界速度即是指在單個時間增量內的生長量。6.6.4 6.6.4 概率性元胞自動機解法概率性元胞自動機解法l在概率性元胞自動機方法中，通過采用權重隨機抽樣方案把確定性積分用統計積分代替，為此，必須把式（6.15）或式（6.16）分解成確定性部分和概率性部分w，即有：l其中 為晶界速度，并且;(6.20)TkQTkpTmknwxxBgbBBexp00(6.21)x00TmknxBTkQTkpwBgbBexp6.6.

35、4 6.6.4 概率性元胞自動機解法概率性元胞自動機解法l模擬應是在空間網格上進行，其給定標度m大于原子尺度。如果有轉變現象發生，則晶粒將按生長（或收縮），而不是b3，為了校正標定尺度l根據m及時間標度（1/v）可知，對統計積分施加這樣一個頻率是不合適的。利用沖擊頻率v0把上述方程歸一化(6.22)mBmTmkvwvnwxx0wxwvvxwvvvnwxxm000000(6.23)6.6.4 6.6.4 概率性元胞自動機解法概率性元胞自動機解法l其中l 由網格大小及選取的沖擊頻率決定； 由溫度及實驗輸入數據決定。例如，晶界特征性質依賴于取向偏差和平面傾角，驅動力取決于所存儲的彈性能和局域曲率。 00vnxmTkQvpmTkQTkpvvwBgbmBgbBexpexp000(6.24)0 xw6.7 CA6.7 CA方法在材料科學中的應用方法在材料科學中的應用l6.7.1再結晶的模擬l(1)動態再結晶過程HY-100鋼的顯微