董衛(wèi)平

DOI：10.16660/j.cnki.1674-098X.2017.14.151

摘 要：該文基于微觀相場模型，推導出L10結構的長程序參數(shù)方程，并研究了L10結構有序度變化。結果表明，L10結構有序度先增大后減小，且L10結構屬于過渡結構，有序度較低，即L10結構不穩(wěn)定，有序度低時轉變成穩(wěn)定的L12結構。

關鍵詞：相場法 L10結構 鎳基高溫合金 有序度

中圖分類號：TG132 文獻標識碼：A 文章編號：1674-098X（2017）05（b）-0151-02

鎳基高溫合金由于其高溫穩(wěn)定性、耐磨性及耐腐蝕性等特性，在各個領域得到廣泛應用，因此許多研究者對其性能、強化過程等進行研究。在鎳基高溫合金研究中隨著溫度濃度不同，可以得到許多不同相，其中最受關注的是強化相γ-L12結構，其實在得到穩(wěn)定的γ-L12結構前，往往有L10過渡結構出現(xiàn)。研究L10過渡結構對進一步研究γ-L12結構有非常重要的意義。微觀相場方程由于其尺度可以達到原子級，且可以動態(tài)觀察不同結構的沉淀變化過程，非常適合該文用于研究L10過渡結構的有序度。

1 相場方程

微觀相場方程對于三元系統(tǒng)，占位幾率公式為，其中A、B和C表示3種不同合金元素。關于微觀相場法的更詳細的公式和過程，可以參考文獻。

2 結果及討論

L10結構模擬得到的原子演化圖通過晶格位置Ni原子占位幾率反映合金的沉淀過程，晶格位置的亮度越高，Ni原子占位幾率越高；亮度越低，占位幾率越低。L10結構有兩種取向：Type-I型和Type-II型，它們的晶體結構沒有區(qū)別，只是取向不同，Type-I型沿著平行于面心立方結構[100]方向旋轉90°即為Type-II型結構。以下公式推導中，占位幾率都是指溶質原子的占位幾率。

合金原子長程序參數(shù)ηM表達式為：

ηu（i，j） （1）

Pu（i，j）為原子在（i，j）晶格位置的占位幾率，當η=1時為完全有序狀態(tài)，即u原子只占據(jù)有序結構中兩種晶格位置其中一種；而η≠1時則u原子占兩種晶格位置；當η=0時u原子處于無序狀態(tài)。用ηu（i，j）和cu（i，j）分別表示u原子（i，j）晶格位置的長程序參數(shù)和原子在晶格中的平均占位幾率（也是u原子的濃度），即表示有序度和原子簇聚程度。根據(jù)Khachaturyan的占位幾率PX（r）函數(shù)，消去cX，可以得出溶質原子的長程序參數(shù)公式：

（2）

應用上述公式對L10結構進行有序度計算，計算結果如圖1，Ni75Al14Cr11合金溫度為923 K時沉淀過程中取1 000步長時演化圖及不同時間步長時L10結構和L12結構長程序參數(shù)變化曲線。從圖1（a）可以看出該鎳基高溫合金時效過程中已經(jīng)得到L10結構，只是L10結構有序度較低，可見L10結構是模擬過程中過渡有序化階段的預析出相，從圖像的顏色可以看出，L10結構中亮度較低，這就表示有序度較低，這也可以從圖1（b）中第三條原點曲線有序度只有0.5得到進一步印證。在1 500步長之后L10結構漸漸消失，轉變成L12結構。且L10結構最大有序度只能達到0.62，之后就降低，而穩(wěn)定的L12結構可以達到更高的有序度。

3 結語

文章有以下結論：推導的長程有序參數(shù)公式可以用于L10結構有序度的計算；模擬的鎳基高溫合金沉淀過程中可以得到預析出相L10結構；L10結構有序度先增大后減小，且最高有序度較低；L10結構不穩(wěn)定，有序度低時就轉變成穩(wěn)定的L12結構。

參考文獻

[1] 王會陽，安云岐，李承宇，等.鎳基高溫合金材料的研究進展[J].材料導報，2011，25（18）：482-486.

[2] Pareige C，Blavette D.Simulation of the FCC→FCC+L12+D022 kinetic reaction[J].Scripta Mater，2001（44）：243-247.

[3] Poduri R，Chen L Q.Computer simulation of atomic ordering and compositional clustering in the pseudobinary Ni3Al-Ni3V system[J].Acta Mater，1998，46（5）：1719-1729.

[4] Zapolsky H，Pareige C，Marteau L，et al.Atom probe analyses and numerical calculation of ternary phase diagram in Ni-Al-V system[J].Calphad，2001（25）：125-134.

[5] Dong W P，Wang Y X，Yang K，et al.Phase field simulation of interatomic potentials for double phase competition during early stage precipitation[J].Chinese Sci.Bull，2011（56）：2055-2060.

[6] Dong W P，Chen Z，Wang Y X.Energy studies of precipitation sequence in Ni75Al10Cr15 alloy based on the phase field theory[J].Sci.China Phys.Mech.Astron，2011，54（5）：821-826.

[7] Dong W P，Wang Y X，Chen Z，et al.Ageing process of pre-precipitation phase in Ni0.75Al0.05Fe0.2 alloy based on phase field method[J].Trans.Nonferrous Met.Soc.China，2011（21）：1105-1111.

[8] Khachaturyan A G.Theory of Structural Transformation in Solids[M].New York：Wiley，1983.

[9] Chen L Q，Khachaturyan A G.Computer simulation of structural transformations during precipitation of an ordered intermetallic phase[J].Acta Metall. Mater，1991（39）：2533-2551.

[10] Wang Y，Chen L Q，Khachaturyan A G.Kinetics of the strain-induced morphological transformation in cubic alloys with miscibility gap[J].Acta Metall. Mater，1993（41）：279-296.

[11] Chen L Q，Simmons J A.Microscopic master equation approach to diffusional transformations in inhomogeneous systems—single-site approximation and direct exchange mechanism[J].Acta Metall.Mater，1994，42（9）：2943-2954.