劉小元，杜立飛，張 蓉，張利民

(西北工業(yè)大學(xué) 理學(xué)院 空間應(yīng)用物理與化學(xué)教育部重點實驗室，西安 710129)

材料的微觀組織結(jié)構(gòu)對材料的性能具有決定行動作用[1?3]。在金屬的凝固過程中，樹枝晶生長是一種最常見的方式，而枝晶間距又是表征凝固組織形態(tài)的重要特征參量，且材料的強(qiáng)度和韌性普遍受到枝晶間距因素的影響，枝晶間距的大小與組織中顯微偏析、夾雜物的形成、微裂紋與縮松的產(chǎn)生等都有密切的關(guān)系，所以材料枝晶間距的研究一直是國內(nèi)外學(xué)者普遍關(guān)注的課題。1997年，黃衛(wèi)東等[4?5]通過定向凝固臺階變速試驗驗證了Al-Cu合金一次枝晶間距選擇的歷史相關(guān)性。翟慎秋等[6]研究了金屬枝晶間距與枝晶生長速度、溫度梯度的關(guān)系，發(fā)現(xiàn)各次枝晶間距與生長速度、溫度梯度基本符合指數(shù)關(guān)系。1999年，常國威等[7]基于M-S理論，提出了電流作用下凝固界面的穩(wěn)定性判據(jù)，并發(fā)現(xiàn)隨著電流增大，枝晶間距以先增大后減小的規(guī)律變化。對于枝晶間距的實驗研究，大多是在定向凝固條件下進(jìn)行的，且研究對象多為溫度梯度、電場強(qiáng)度、枝晶生長速度等與枝晶間距的關(guān)系[8?11]。但是，在材料的實際凝固過程中，往往要經(jīng)歷許多非規(guī)則的橫截面，例如由定向凝固獲得的渦輪葉片就是有不規(guī)則界面的影響下得到的復(fù)雜微觀結(jié)構(gòu)成品。這種不規(guī)則界面的存在對凝固過程有兩個重要的作用：首先，熱量的傳導(dǎo)因界面的存在受到明顯影響，不同導(dǎo)熱性質(zhì)的橫向界面有著不同的導(dǎo)熱規(guī)律；其次，界面的存在會明顯影響的固液界面的推進(jìn)，且不同形狀和性質(zhì)的界面對微觀組織的形成存在不同程度的影響。另外，在合金凝固過程中存在橫向界面時，溶質(zhì)場的分布也會顯著變化。FABIETTI等[12]通過研究橫向限制存在時定向凝固條件下對丁二腈凝固過程，發(fā)現(xiàn)橫向限制對枝晶尖端生長速度及枝晶間距存在重要影響。王華明等[13]研究了側(cè)向約束條件下單晶 Ni基高溫合金的凝固行為及溶質(zhì)再分配規(guī)律，發(fā)現(xiàn)側(cè)向約束條件對單晶Ni基高溫合金的凝固組織、溶質(zhì)偏析及溶質(zhì)再分配均具有極為強(qiáng)烈的影響，提出側(cè)向約束條件即試樣斷面變化對單晶 Ni基高溫合金凝固過程的影響實質(zhì)上是局部凝固工藝參數(shù)的影響。由此可見，研究橫向限制條件下的凝固行為對理解界面存在對凝固過程中溶質(zhì)及溫度傳導(dǎo)規(guī)律具有重要意義：溶質(zhì)及熱量傳導(dǎo)的改變將直接影響材料微觀結(jié)構(gòu)的形成，并最終導(dǎo)致材料具有特異的性能。但是由于實驗條件的限制，目前的實驗手段很難精確地測量材料凝固過程的諸多細(xì)節(jié)，而采用計算模擬的方法不僅可以有效地節(jié)約成本，還能夠有效地反映完整的凝固過程，并且獲得更加豐富的過程參量及數(shù)據(jù)，有助于在理論上研究各種實驗條件對凝固過程的影響，但以往的研究大都集中于各種復(fù)雜條件下的凝固過程[14?19]，而對于橫向限制存在條件下的模擬研究卻鮮有報道，所以本文將以純金屬Ni為研究對象，利用相場方法，模擬不同性質(zhì)橫向擋板作用下的枝晶微觀形貌過程，并著重討論橫向擋板限制對枝晶間距調(diào)整的影響。

1 相場模型及數(shù)值求解

純金屬凝固的非等溫相場方程包括兩個與空間r及時間t相關(guān)的場量，即相場量φ(r,t)和溫度場量T(r,t)，序參量φ(r,t)表征體系的有序程度，定義液相內(nèi)φ=0，在固相內(nèi)φ=1，在固液界面內(nèi)φ∈ ( 0,1)。T(r,t)為溫度場。

相場控制方程是基于 Ginzburg-Landau自由能泛函的變分法推導(dǎo)得到的，純金屬的 Ginzburg-Landau自由能形式如下：

式中：ε為界面自由能梯度系數(shù)；(,)fmφ為雙阱勢能函數(shù)，其具體形式取為

可以看出，式(2)滿足：f′(0,m)=f′(1,m)≡0，即在固液兩相達(dá)到極小。在此基礎(chǔ)上，引入各向異性參數(shù)以及噪聲項，最終得到相場控制方程：

式中：ε=ε01+δcos[j(θ?θ0)]）為各向異性參數(shù)，其中δ表示各向異性強(qiáng)度，j為各向異性模數(shù)，表征枝晶最優(yōu)生長取向的方向數(shù)，θ為最優(yōu)生長方向與參考方向的夾角，θ0為參考角度；τ(弛豫時間)是一個很小的負(fù)常數(shù)；aφ(1?φ)χ為噪聲項，χ∈［? 1/2，1/2］為隨機(jī)數(shù)，a為噪聲強(qiáng)度；m為控制自由能勢阱的參數(shù)，m是溫度相關(guān)的函數(shù)，定義為

式中：Te表示平衡溫度；α為負(fù)常數(shù)值，取α＜1，保證＜1/2；γ為相變驅(qū)動力控制參數(shù)。

式(3)的詳細(xì)推導(dǎo)可參見文獻(xiàn)[20]。

由熱擴(kuò)散方程可得溫度場方程：

式中：K為無量綱潛熱。

以純金屬Ni為研究對象，利用有限差分法對無量綱化后的式(3)和(5)進(jìn)行離散求解，計算域的大小為800 ?x× 8 00?y，相關(guān)物性參數(shù)及模擬參數(shù)取值見表1。

表1 模擬中用到的參數(shù)取值Table 1 Parameters used in simulations

模擬初始時刻計算域的無量綱溫度T=0，擋板的初始無量綱溫度設(shè)置為T=?0.1，在計算區(qū)域底部引入等距分布的晶核，在晶核范圍內(nèi)場量初始值1=φ，在晶核范圍之外，相場量初始值0=φ，且不隨時間發(fā)生改變，相場與溫度場的邊界條件均為周期性邊界條件。在滿足上述設(shè)定的情況下，本研究模擬了不同形狀橫向擋板限制下，不同初始枝晶間距條件下的枝晶生長過程，模擬中擋板的形狀及尺寸設(shè)置見圖1。

2 模擬結(jié)果與討論

2.1 橫向擋板限制下的枝晶演化過程

圖2(a)～(f)所示為上三角型擋板限制存在下的枝晶微觀形貌演化過程。初始枝晶間距d0=100?x，擋板直角邊長度d2=300?x，擋板高度w1=100?x。在枝晶未接觸到擋板時，主枝晶臂均沿豎直方向自由生長，二次枝晶臂發(fā)育被主枝晶軸的快速發(fā)育限制。主枝晶尖端接觸到擋板時，其生長直接被限制，而通過擋板間隙的枝晶繼續(xù)向上發(fā)育，如圖2(b)所示。計算中擋板設(shè)置為完全潤濕，即擋板溫度低于熔體溫度，故固相接觸擋板后將沿?fù)醢逍边吙焖偕L，如圖2(c)所示，固相在擋板 45°角處發(fā)育出新枝晶中心，與通過擋板間隙的枝晶形成競爭生長。而且，固相在沿?fù)醢逍边吙焖偕L的同時，還發(fā)育出了新的縱向枝晶向上生長，如圖2(d)所示，且新發(fā)育縱向枝晶的間距明顯減小，并且?guī)缀鯖]有二次枝晶形成。最終，通過擋板間隙枝晶，擋板 45°角處新發(fā)育枝晶及擋板斜邊新發(fā)育的枝晶臂競爭生長，形成與初始枝晶微觀結(jié)構(gòu)完全不同的不規(guī)則的組織，如圖2(d)所示，此時自擋板斜邊處新發(fā)育枝晶臂受到穿過擋板間隙的枝晶生長的影響，生長方向相對于界面能各向異性決定的最優(yōu)生長發(fā)現(xiàn)發(fā)生了一定的偏轉(zhuǎn)?？梢姡瑱M向擋板的存在使擋板上下兩個區(qū)域的枝晶形貌發(fā)生很大變化：上部新發(fā)育的枝晶形貌變得不規(guī)則，且新發(fā)育枝晶的平均枝晶間距明顯小于初始枝晶間距，同時在擋板尖端出發(fā)育的新枝晶臂具有發(fā)達(dá)的二次枝晶。圖 2(g)所示為時間為10 000t?時刻的溫度場分布情況，可見擋板附近存在明顯的溫度梯度，使得擋板附近的金屬以平面方式生長。但是由于模擬中隨機(jī)噪聲的引入，能夠誘發(fā)快速生長金屬固相表面的枝晶發(fā)育，形成新的枝晶臂向上快速生長，與穿過擋板的枝晶競爭生長。這些與FABIETTI等[12]得到實驗結(jié)果基本一致，如圖 2(h)所示。由此可以看出，橫向限制擋板的存在可以顯著改變晶體微觀結(jié)構(gòu)的發(fā)展過程，且對枝晶間距的調(diào)整具有明顯的作用。為此將設(shè)置不同初始間距的枝晶，來研究橫向擋板對其枝晶間距的調(diào)整作用。計算結(jié)果見圖3。

從圖3看出，初始枝晶間距較大時，枝晶自擋板間隙通過，二次枝晶臂接觸擋板后沿?fù)醢逍泵婵焖偕L，且發(fā)育出新的枝晶臂，與初始枝晶的二次枝晶臂競爭生長，明顯改變了微觀形貌演化，且經(jīng)過擋板后的枝晶間距明顯減小。初始枝晶間距與擋板間隙相差不大時，初始枝晶的發(fā)展受到限制，但枝晶繞過擋板并在其 45°角處發(fā)育出較為發(fā)達(dá)的新的主枝晶軸，且在擋板斜面處仍然會發(fā)育出豎直方向生長的枝晶臂，從而改變枝晶間距，如圖3(b)和(c)所示。當(dāng)擋板間距小于枝晶間距時，部分初始枝晶通過擋板間隙繼續(xù)向上生長，但繞過擋板新發(fā)育的枝晶間距明顯小于初始枝晶間距，且自擋板 45°角處發(fā)育的新枝晶具有相對發(fā)達(dá)的二次枝晶臂。可見，三角形擋板的存在對不同初始間距的枝晶生長都有一定的作用。

圖1 4種不同形狀的擋板尺寸圖Fig. 1 Sketch maps of four different constrains: (a) Rectangle; (b) Triangle with sharp corner on top; (c) Triangle with sharp corner at bottom; (d) Trapezoid

圖2 上三角擋板存在時枝晶微觀結(jié)構(gòu)的演化過程Fig. 2 Microstructure evolution of dendrite growth with lateral constrains of triangle: (a) 2 000?t; (b) 4 000?t; (c) 6 000?t; (d) 8 000?t; (e) 10 000?t; (f) 12 000?t; (g) Temperature distribution at 10 000?t; (h) Experiment result[12]

圖3 不同初始間距枝晶在上三角形擋板存在時的枝晶微觀結(jié)構(gòu)Fig. 3 Microstructures of dendrites with different primary spacings and lateral constrains of triangle with sharp corner on top:(a) d0=800?x; (b) d0=400?x; (c) d0=200?x; (d) d0=100?x

2.2 不同形狀橫向擋板對枝晶間距的影響

圖4 所示為下三角型擋板限制存在時，不同初始枝晶間距條件下的枝晶形貌，擋板尺寸設(shè)置為d3=300?x。對于d0=800?x的情況，主枝晶軸從擋板通過，其二次枝晶臂被被擋板限制，主枝晶通過擋板后其二次枝晶重新發(fā)育，基本保持原有形貌繼續(xù)發(fā)展，此時擋板的存在并不能夠改變枝晶生長過程中二次枝晶的發(fā)育。隨著初始枝晶間距減小，主枝晶軸在接觸擋板后沿其斜邊快速生長，且在擋板 45°角處形成新的晶核，沿左上右3個方向發(fā)育出新的枝晶主軸，發(fā)育形成具有發(fā)達(dá)二次枝晶臂的新的枝晶結(jié)構(gòu)，且由于枝晶間隙相對較小，縱向發(fā)育的枝晶軸完全限制了初始枝晶的生長，新發(fā)育枝晶結(jié)構(gòu)具有統(tǒng)一的生長模式，即 45°角處新發(fā)育的具有發(fā)展二次枝晶臂枝晶與兩側(cè)柱狀晶的競爭生長。

通過以上分析可見，晶體生長在橫向限制擋板的拐角處會重新形成新的枝晶發(fā)展中心，導(dǎo)致枝晶生長模式的改變，從而改變枝晶間距，故有必要進(jìn)一步研究具有多個擋板拐角的橫向限制擋板的存在對枝晶發(fā)育過程的影響。

圖5所示為受矩形擋板限制，不同初始枝晶間距條件下的枝晶形貌，擋板尺寸為d1=300?x，w1=100?x。矩形擋板具有相對規(guī)則的形狀，當(dāng)固相接觸擋板后會沿著擋板邊緣快速生長，且在擋板上拐角處形成二次枝晶臂較發(fā)達(dá)的新主枝晶軸。在初始枝晶間距較大時，繞過擋板間隙的枝晶間距明顯減??；當(dāng)初始枝晶間距等于擋板間隙時(見圖5(c))，通過擋板間隙的枝晶臂仍然有比較充足的發(fā)展空間，二次枝晶臂得到一定的發(fā)育，擋板上方枝晶間距仍然會減小；當(dāng)初始枝晶間距小于擋板間距時，中心枝晶臂通過擋板間隙按原軌跡繼續(xù)生長，并與兩側(cè)新生成的枝晶臂形成競爭生長，最后發(fā)育成不規(guī)則的結(jié)構(gòu)。

圖4 不同初始間距枝晶在下三角形擋板存在時的枝晶微觀結(jié)構(gòu)Fig. 4 Microstructures of dendrites with different primary spacings and lateral constrains of triangle with sharp corner at bottom:(a) d0=800?x; (b) d0=400?x; (c) d0=200?x; (d) d0=100?x

圖5 不同初始間距枝晶在矩形擋板存在時的枝晶微觀結(jié)構(gòu)Fig. 5 Microstructures of dendrites with different primary spacings and lateral constrains of rectangle: (a) d0=800?x; (b) d0=400?x;(c) d0=200?x; (d) d0=100?x

圖6 和7所示分別為兩種不同尺寸梯形擋板限制存在時，不同初始枝晶間距條件下的枝晶形貌，擋板尺寸設(shè)置分別為d4=300?x、w2=100?x和d4=300?x、w2=200?x。相比于矩形擋板，梯形擋板間隙上的兩個邊存在一傾斜度，因此，枝晶在沿?fù)醢蹇焖偕L時，可以在其上下兩個角處均發(fā)育出新的枝晶主軸。對于所有的情況，在梯形的上135°處均發(fā)育出了新的枝晶軸。對于枝晶主軸被阻隔的生長情況(排除d0=800?x情況)，當(dāng)擋板下間隙小于或等于枝晶初始間距時，在擋板下 45°角處也能夠發(fā)育出新枝晶軸，且新的枝晶間距小于初始枝晶間距；當(dāng)擋板下間隙大于枝晶初始間距時，在擋板下 45°角處，由于競爭生長的作用，原主枝晶軸側(cè)向二次枝晶臂阻礙了出擋板下 45°角處新枝晶軸的生長，導(dǎo)致?lián)醢彘g隙附近的枝晶間距變大，但由于擋板上側(cè)柱狀晶的形成，繞過擋板后的平均枝晶間距仍然減小，如圖6(d)所示。當(dāng)擋板斜邊相對較大時，如圖7所示，擋板斜邊上也會有新的柱狀晶形成，明顯改變繞過擋板后的枝晶間距。由此可見，不同尺寸梯形擋板對于枝晶間距調(diào)整具有不同的作用。

2.3 不同形狀擋板限制對枝晶間距調(diào)整的作用

由以上的計算分析可知，對于不同初始間距枝晶，不同形狀橫向限制對枝晶間距改變具有不同的效應(yīng)。盡管橫向擋板的存在會使初始的規(guī)則枝晶在繞過擋板后形成不規(guī)則的微觀組織，但是新發(fā)育的枝晶間距明顯減小，這使得經(jīng)過擋板后的枝晶的平均間距發(fā)生明顯改變。假設(shè)繞過橫向擋板后新發(fā)育枝晶的平均間距為λ1，其大小隨初始枝晶間距的改變見圖8(a)；圖8(b)所示為新的枝晶間距與初始枝晶間距的比值。可見，相對于初始枝晶間距λ0，上三角形擋板對枝晶間距的改變最明顯，同時下三角形及矩形的橫向限制對枝晶間距具有一定的控制作用，即對于不同初始間距枝晶繞過這兩類擋板后具有基本相同的重新發(fā)育模式，枝晶間距只與擋板尺寸有關(guān)，與初始枝晶間距關(guān)系不大。不同尺寸梯形擋板對于枝晶間距具有完全不同的調(diào)整作用，即通過調(diào)整梯形擋板尺寸比例可以有效改變枝晶間距。

圖6 不同初始間距枝晶在梯形擋板存在時的枝晶微觀結(jié)構(gòu)Fig. 6 Microstructures of dendrites with different primary spacings and lateral constrains of trapezoid at d4=300?x and w2=100?x:(a) d0=800?x; (b) d0=400?x; (c) d0=200?x; (d) d0=100?x

圖7 不同初始間距枝晶在梯形擋板存在時的枝晶微觀結(jié)構(gòu)Fig. 7 Microstructures of dendrites with different primary spacings and lateral constrains of trapezoid at d4=300?x and w2=200?x:(a) d0=800?x; (b) d0=400?x; (c) d0=200?x; (d) d0=100?x

圖8 不同形狀擋板對枝晶間距調(diào)整的作用Fig. 8 Functions of lateral constrains with different shapes on dendrite arm spacing changes: (a) λ1[?x]—λ0[?x]; (b) λ1/λ0—λ0[?x];A—Rectangle; B—Triangle with sharp corner on top; C—Triangle with sharp corner at bottom; D—Trapezoid, shown in Fig. 6; E—Trapezoid, shown in Fig. 7

3 結(jié)論

1) 凝固過程中橫向擋板的存在能夠明顯改變枝晶發(fā)育過程：枝晶生長過程中碰到橫向擋板時, 枝晶的生長會直接被擋板下邊緣限制, 但固相繞過擋板后在擋板拐角處形成新的枝晶發(fā)育中心，且固相繞過擋板后會沿其邊緣快速發(fā)展, 并發(fā)育出新的枝晶臂, 與穿過擋板間隙的枝晶形成競爭生長，最終形成不規(guī)則的枝晶微觀組織，與相關(guān)的實驗結(jié)構(gòu)基本一致。

2) 不同形狀的橫向擋板對枝晶間距調(diào)整具有不同的作用：上三角形擋板對枝晶間距的改變最明顯；下三角形及矩形的橫向限制對枝晶間距具有一定的控制作用，即對于不同初始間距枝晶繞過這兩類擋板后具有基本相同的重新發(fā)育模式，枝晶間距只與擋板尺寸有關(guān)，與初始枝晶間距關(guān)系不大。不同尺寸梯形擋板對于枝晶間距具有完全不同的調(diào)整作用，通過調(diào)整梯形擋板尺寸比例可以有效實現(xiàn)改變枝晶間距的作用。

[1]ASTA M, BECKERMANN C, KARMA A, KURZ W,NAPOLITANO A, PLAPP M, PURDY G, RAPPAZ M,TRIVEDI R, Solidification microstructures and solid-state parallels: Recent developments, future directions[J]. Acta Materialia, 2009, 57: 941?971.

[2]YAN Z M, LI X T, CAO Z Q, ZHANG X A, LI T J. Grain refinement of horizontal continuous casting of the CuNi10Fe1Mn alloy hollow billets by rotating magnetic field(RMF)[J]. Materials Letters, 2008, 62: 4389?4392.

[3]YASUDA H, TOH T, IWAI K, MORITA K. Recent progress of EPM in steelmaking, casting, and solidification processing[J].ISIJ International, 2007, 47(4): 619?626.

[4]林 鑫, 黃衛(wèi)東, 潘清躍, 丁國陸, 薛玉芳, 周堯和. Al-4.5Cu單晶定向凝固一次枝晶間距研究[J]. 金屬學(xué)報, 1997, 33(11):1140?1146.LIN Xin, HUANG Wei-dong, PAN Qing-yue, DING Guo-lu,XUE Yu-fang, ZHOU Yao-he. Research on primary dendritic spacing in unidirectional solidification of Al-4.5Cu single crysta[J]. Acta Metallurgical Sinica, 1997, 33(11): 1140?1146.

[5]李林蓄, 林 鑫, 王 猛, 黃衛(wèi)東. Al-3.6%Cu合金定向凝固一次枝晶間距選擇[J]. 鑄造技術(shù), 2008, 193(7): 60?64.LI Kin-xu, LIN Xin, WANG Meng, HUANG Wei-dong. Selection of primary dendritic spacing during directional solidification of Al-3.6%Cu[J]. Foundry Technology, 2008, 29(7): 60?64.

[6]翟慎秋, 魏亞杰. Sn-Pb合金的三次樹枝晶生長[J]. 淄博學(xué)院學(xué)報: 自然科學(xué)與工程版, 2000, 2(3): 38?46.ZHAI Shen-qiu, WEI Ya-jie, Tertiary dendrite growth in Sn-Pb alloy[J]. Journal of Zibo University: Natural Science and Engineering Edition, 2000, 2(3): 38?46.

[7]常國威, 袁軍平, 王自東, 胡漢起, 薛慶國. 電流改變定向凝固單向合金枝晶間距機(jī)理[J]. 北京科技大學(xué)學(xué)報, 1999, 21(2):175?178.CHANG Guo-wei, YUAN Jun-ping, WANG Zi-dong, HU Han-qi, XUE Qing-guo. Mechanism of dendrite armapace of unidirectional solidification single-phase alloy with influence of electric current[J]. Journal of University of Science and Technology of Beijing, 1999, 21(2): 175?178.

[8]EASTON M, DAVISON C, JOHN D Jr. Effect of alloy composition on the dendrite arm spacing of multicomponent aluminum alloys[J].Metallurgical Transactions A, 2010, 41(6): 1528?1538.

[9]楊初斌, 劉 林, 趙新寶, 劉 剛, 張 軍, 傅恒志. 〈001〉和〈011〉取向 DD407單晶高溫合金枝晶間距和微觀偏析[J]. 金屬學(xué)報, 2011, 47(10): 1246?1250.YANG Chu-bin, LIU Lin, ZHAO Xin-bao, LIU Gang, ZHANG Jun, FU Heng-zhi. dendrite arm spacings and microsegregations in 〉〈001 and 〉〈011 orientated single crystal superalloys DD407[J]. Acta Metallurgical Sinica, 2011, 47(10): 1246?1250.

[10]HU Xiao-wu, LI Shuang-ming, GAO Si-feng, LIU Lin, FU Heng-zhi. Effect of melt convection on primary dendrite arm spacing in directionally solidified Pb-26%Bi hypo-peritectic alloys[J]. Transactions of Nonferrous Metals Society of China,2011, 21(1): 65?71.

[11]AN G Y, LIU L X. Dendrite spacing in unidirectional solidified Al-Cu alloy[J]. Journal of Crystal Growth, 1987, 80(2): 383?392.

[12]FABIETTI L M, SEETHARAMAN V, TRIVEDI R. The development of solidification microstructures in the presence of lateral constraints[J]. Metallurgical Transactions A, 1990, 21:1299?1310.

[13]王華明, 唐亞俊, 張靜華, 李英敖, 張志亞, 于 洋, 胡壯麒.單晶Ni基高溫合金側(cè)向約束下的凝固行為[J]. 材料科學(xué)進(jìn)展,1993, 7(2): 99?104.WANG Hua-ming, TANG Ya-jun, ZHANG Jing-hua, LI Ying-ao,ZHANG Zhi-ya, YU Yang, HU Zhuang-qi. Solidification behavior of a single-crystal Ni-base supper alloy under lateral constrains[J]. Materal Science Progress, 1993, 7(2): 99?104.

[14]李俊杰, 王錦程, 許 泉, 楊根倉. 外來夾雜物顆粒對枝晶生長形態(tài)影響的相場法研究[J]. 物理學(xué)報, 2007, 56(3): 1514?1519.LI Jun-jie, WANG Jin-cheng, XU Quan, YANG Gen-cang.Effect of foreign particles on the dendritic growth in phase-field theory[J]. Acta Physica Sinica, 2007, 56(3): 1514?1519.

[15]石玉峰, 許慶彥, 柳百成. 定向凝固共晶生長的元胞自動機(jī)數(shù)值模擬[J]. 金屬學(xué)報, 2012, 48(1): 41?48.SHI Yu-feng, XU Qing-yan, LIU Bai-cheng. Simulation of eutectic growth in directional solidification by cellular automaton method[J]. Acta Metallurgical Sinica, 2012, 48(1): 41?48.

[16]盧艷麗, 陳 錚, 來慶波, 張 靜. Ni75Cr25?xAlx合金中L12相和D022相形核孕育期的微觀相場模擬[J]. 物理學(xué)報, 2009,58(13): 319?326.LU Yan-li, CHEN Zheng, LAI Qing-bo, ZHANG Jing.Microscopic phase field simulation on the nucleation incubation period ofL12andD022phases in Ni75Cr25?xAlxalloy[J]. Acta Physica Sinica, 2009, 58(13): 319?326.

[17]杜立飛, 張 蓉, 刑 輝, 張利民, 張 洋, 劉 林. 橫向限制下凝固微觀組織演化的相場法模擬[J]. 物理學(xué)報, 2013,62(10): 106401-1?106401-6.DU Li-fei, ZHANG Rong, XING Hui, ZHANG Li-min, ZHANG Yang, LIU Lin. Phase-field simulation of solidified microstructure evolution in the presence of lateral constraint[J].Acta Physica Sinica, 2013, 62(10): 106401-1?106401-6.

[18]LU Shi-qiang, LI Xin, WANG Ke-lu, DONG Xian-juan, FU M W. High temperature deformation behavior and optimization of hot compression process parameters in TC11 titanium alloy with coarse lamellar original microstructure[J]. Transactions of Nonferrous Metals Society of China, 2013, 23(2): 353?360.

[19]WANG Gang, ZENG De-chang, LIU Zhong-wu. Phase field calculation of interface mobility in a ternary alloy[J]. Transactions of Nonferrous Metals Society of China, 2012, 22(6): 1711?1716.

[20]KOBAYASHI R. Modeling and numerical simulations of dendritic crystal growth[J]. Physica D, 1993, 63: 410?423.