陳　瑞，許慶彥，柳百成（清華大學 材料學院 先進成形制造教育部重點實驗室，北京 100084）

Al-7Si-0.36Mg合金定向凝固一次枝晶臂間距實驗和模擬

陳瑞，許慶彥，柳百成
（清華大學 材料學院 先進成形制造教育部重點實驗室，北京 100084）

通過Al-7Si-0.36Mg合金定向凝固實驗和元胞自動機模型，開展定向凝固枝晶形貌演化和一次枝晶臂間距選擇過程的實驗和模擬。結果表明：在給定的凝固條件下，一次枝晶臂間距范圍是一個連續(xù)的變化區(qū)間。在恒定溫度梯度和不同凝固速度條件下，測得Al-7Si-0.36Mg合金一次枝晶臂間距上限值（λmax）、下限值（λmin）和平均值（λave）以及生長速率之間的關系，且上限值和下限值的比值接近3。模擬結果與實驗結果的吻合程度明顯優(yōu)于Hunt-Lu等解析模型的預測結果，表明CA模型在枝晶定向凝固過程枝晶形貌演化模擬和枝晶臂間距預測等方面的準確性。結合模擬研究和文獻調研分析影響定向凝固一次枝晶臂選擇的因素，包括抽拉速度v、溫度梯度G、界面能大小、溶質擴散系數DL、枝晶生長取向與熱流方向的偏離角度θ等。

定向凝固；枝晶臂間距；生長速度；元胞自動機；數值模擬

定向凝固鑄件中枝晶形貌及特征直接影響著鑄件最終力學性能，因此，通過深入了解和掌握定向凝固過程中枝晶的演化規(guī)律，有助于控制枝晶的形成過程，進而達到對鑄件性能進行有效預測和控制的目的[1]。枝晶臂間距（尤其是一次枝晶臂間距），是用于表征定向凝固枝晶細化程度的重要參數，枝晶臂間距的大小直接影響著溶質成分分布、中間相、共晶相、雜質相的尺寸和分布以及裂紋和縮松等缺陷的形成[2-4]，因此，研究定向凝固不同工藝條件下枝晶臂間距的演化過程并進行準確的定量預測，對控制凝固組織、提高鑄件性能具有十分重要的意義。

過去的幾十年中，許多學者基于一些假設條件建立了枝晶臂間距與凝固參數及合金特性相關聯的解析模型[5-7]。HUNT[5]、KURZ等[6]、TRIVEDI[7]分別通過假設不同的枝晶尖端形狀（球形或是橢圓形）和完全一致的枝晶列排布，并利用臨界穩(wěn)定性判據（Marginal stability criterion）建立了穩(wěn)態(tài)條件下定向凝固枝晶一次枝晶臂間距λ1與固液界面溫度梯度G和生長速度v的單值函數關系。這些模型可以用式（1）統(tǒng)一表示：

式中：A表示與材料特性參數相關的函數；C0是合金成分。從式（1）可以看到，在材料和成分固定的情況下，即A和C0是定值，一次枝晶臂間距λ1只由固液界面的溫度梯度G和生長速度v決定，并且在同一種條件下，λ1具有唯一值。隨著枝晶生長界面穩(wěn)定性、枝晶側向分枝機制等因素對枝晶臂間距影響研究的深化[8-10]以及一些實驗數據的驗證[11]，研究人員逐漸意識到枝晶的演化不僅受實驗條件參數的影響，而且與演化歷史密切相關，所以在同一凝固參數下，一次枝晶臂間距λ1并不是一個固定值，而是一個連續(xù)變化區(qū)間。HUNT等[12]通過分析凝固界面前沿溶質橫向擴散所引起的交互作用影響以及基于枝晶陣列穩(wěn)定性分析，推導出了一次枝晶臂間距與溫度梯度、凝固速度等因素的關系，確定了穩(wěn)態(tài)和非穩(wěn)態(tài)條件下一次枝晶臂間距的下限值[12]：

為了能更準確描述枝晶臂間距的演化過程，模型中必須要考慮枝晶臂的淹沒和側向分枝的作用，而這兩種機制主要受枝晶凝固界面穩(wěn)定性（界面張力和界面各向異性）、溶質交互作用以及枝晶復雜的空間擴展結構等因素的影響。顯微組織數值模擬技術作為一種能準確描述這些因素變化的方法，已經成功應用在定向凝固過程枝晶形貌和枝晶臂演化等方面的預測[13-18]，并且證明了數值模擬預測的結果要比上述解析模型預測的結果更為準確[13, 15]。目前，顯微組織數值模擬的主要方法有元胞自動機（CA）法和相場（PF）法，其中，CA方法具有較高的計算效率，通過耦合熱擴散、溶質擴散以及枝晶生長動力學等，建立較復雜的描述枝晶固/液界面形態(tài)以及溶質分布的數學方程，從而在更大尺度上實現對枝晶形貌演化過程的模擬，近些年來發(fā)展較快，在實際中得到了廣泛的應用。WANG等[13]建立了溶質擴散控制的CA-FD模型，該模型考慮了成分過冷和曲率過冷，并改進了GANDIN 等[19]提出的用于模擬晶粒生長的偏心算法（Decentered Square Algorithm），實現了任意取向枝晶生長的模擬，他們采用該模型詳細研究二元鎳基高溫合金定向凝固過程中初始枝晶臂數目對穩(wěn)態(tài)枝晶整列一次枝晶臂間距的影響，同時也分析了凝固速度對枝晶生長的影響，結果表明初始枝晶臂數目對穩(wěn)態(tài)時的枝晶臂間距影響不大，而凝固速度對枝晶臂間距的影響很大，同時也給出了枝晶臂間距上限值與下限值呈接近3倍的關系。隨后，單博煒等[15]也建立了相應的CA模型模擬分析了定向凝固一次枝晶臂間距的選擇過程，分析了初始晶核數目和變抽拉速度對于一次枝晶臂間距允許范圍的影響，結果表明一次枝晶臂間距與這兩種數值實驗方法無關，同時給出了一次枝晶臂間距的上線值和下限值相對于凝固速度的冪函數關系式，并與實驗結果進行了良好對比。ZHANG等[16]采用3D-CA模型模擬了Al-11.6Cu-0.85Mg（質量分數，％）合金定向凝固枝晶演變過程，指出三次枝晶臂能否發(fā)展成為穩(wěn)態(tài)一次枝晶臂主要由固液界面穩(wěn)定性以及相鄰一次枝晶臂間溶質交互作用程度共同決定，同時，模擬得到的平均一次枝晶臂間距與實驗結果相吻合。目前，大多數的研究局限于分析抽拉速度和溫度梯度等工藝參數的影響，對于影響一次枝晶臂間距的其他因素研究的較少。

本文作者首先通過實驗獲得不同凝固速度下Al-7Si-0.36Mg合金的定向凝固組織和枝晶臂間距的統(tǒng)計數據，然后基于之前建立的CA模型基礎之上[20-21]，預測了一次枝晶臂間距的上限值和下限值，并模擬了該合金定向凝固過程中不同凝固速度下枝晶臂間距演化過程，同時與實驗結果進行對比。此外，還對影響定向凝固一次枝晶臂選擇的諸多因素進行了詳細分析。

1　實驗

1.1實驗方法

本實驗用于定向凝固實驗合金的成分為Al-7Si-0.36Mg（質量分數，％）。由于該合金不同組元之間的密度差異小，自然對流對枝晶生長的影響相對較小。定向凝固實驗在Bridgman 定向凝固爐中完成。首先，把直徑為10mm，長度為120mm的試棒放入到高純石墨坩堝中，在含有氬氣保護的氣氛中把試樣加熱到840℃，保溫0.5 h，然后把試樣按設定的速度往下抽拉到裝有Ga-In-Sn液態(tài)金屬的冷卻罐中，使合金從下往上長出柱狀晶組織。實驗過程中溫度梯度為15 K/mm[22]，完成了4種不同抽拉速度下的定向凝固試樣，分別為50、100、150、200μm/s。截取高度在75～85mm區(qū)間的試棒用于金相試樣的制取并在偏振光顯微鏡下觀察縱截面和橫截面的枝晶組織。通過在橫截面上采用截線法來統(tǒng)計各種凝固條件下柱狀晶的一次枝晶臂間距λ1，每種凝固條件都統(tǒng)計200個數據點，從而可以通過統(tǒng)計的方式得到一次枝晶臂間距λ1的概率分布圖以及一次枝晶臂間距最大值λmax、平均值λave、最小值λmin等特征參數。

1.2實驗結果

圖1所示為4種不同抽拉速度下凝固試樣橫截面的定向凝固組織。從圖1中可以看出，在一次枝晶臂上發(fā)展出了明顯的二次枝晶臂，在二次枝晶臂上出現了三次枝晶臂。對于Al-7Si-0.36Mg合金來說，該合金的初生相凝固溫度區(qū)間相對較窄，在共晶形成之前，熔體中的枝晶組織沒有足夠的時間進行枝晶臂的粗化、合并，所以該合金最終的初生相呈現為比較發(fā)達的樹枝晶形貌。比較不同抽拉速度下的枝晶組織可以明顯看到，隨著抽拉速度從50增加到150μm/s，枝晶發(fā)生明顯的細化（見圖1（a）～（c）），而當抽拉速度從150增大到200μm/s時，枝晶的細化效果相對不是很明顯。此外，從圖1中也可以看到在同一凝固條件下，相鄰一次枝晶臂間距也并不是固定的。為了定量表征不同凝固速度下枝晶間距的變化，圖2所示為一次枝晶臂間距的概率分布圖。由圖2可見，一次枝晶臂間距存在較大的穩(wěn)態(tài)區(qū)間，且隨著抽拉速度的增大，一次枝晶臂穩(wěn)態(tài)區(qū)間減小，概率分布更加趨于正態(tài)分布。表1所列為不同凝固條件下一次枝晶臂間距最大值λmax、平均值λave、最小值λmin以及最大值與最小值的比值λmax/λmin。由表1可以看出，λmax/λmin的比值大致分布在3.0左右。許多文獻從實驗和模擬的角度對一次枝晶臂的穩(wěn)態(tài)區(qū)間進行了相關研究。對于高溫合金，QUESTED等[23]、黃太文等[24]、WANG等[13]的研究表明穩(wěn)態(tài)時一次枝晶臂上限值與下限值的比值λmax/λmin接近3。對于丁二腈-2.5％乙醇（質量分數）透明合金，HUANG等[11]和單博煒等[15]研究結果顯示，在不同凝固速度下，λmax/λmin的比值在2左右波動，而對于丁二腈-0.055.％丙酮溶液，HAN等[25]研究發(fā)現，λmax/λmin的比值約為1.3。由此可見，一次枝晶臂間距上下限的比值與合金體系密切相關。

圖1　Al-7Si-0.36Mg合金定向凝固不同凝固速度下柱狀晶橫截面的枝晶形貌Fig.1　Dendrite morphologies on transverse sections of Al-7Si-0.36Mg alloy directional solidified at different growth velocities：（a）50μm/s；（b）100μm/s；（c）150μm/s；（d）200μm/s

圖2　不同凝固速度下一次枝晶臂間距λ1的概率分布Fig.2　Size distribution of primary dendrite spacing λ1solidified at different growth velocities：（a）50μm/s；（b）100μm/s；（c）150μm/s；（d）200μm/s

表1　不同凝固速度下Al-7Si-0.36Mg合金一次枝晶臂間距最大值λmax、平均值λave、最小值λmin以及λmax/λmin比值Table 1　λmax, λave, λminand λmax/λminof characteristic primary dendrite spacing at different growth velocities in directional solidification of Al-7Si-0.36Mg alloy

2　模擬結果與分析

利用CA模型對Al-7Si-0.36Mg合金定向凝固條件下枝晶生長過程開展相應數值模擬研究，模擬采用的算法在文獻[20-21]中進行了詳細描述，合金的熱物性參數請參見文獻[22]。本文模擬內容主要為：首先通過設置不同初始一次枝晶臂間距來考察不同凝固速度下穩(wěn)態(tài)一次枝晶臂間距的上下限值，然后模擬了圖1中不同凝固條件下的枝晶選擇過程，預測平均一次枝晶臂間距，并與實驗結果進行對比，最后利用CA模型對影響枝晶臂選擇過程的因素進行了分析。

2.1一次枝晶臂間距λmax和λmin預測

由于一次枝晶臂間距的分布區(qū)間決定了定向凝固組織的均勻程度，而一次枝晶臂間距上限值λmax和下限值λmin是用于體現該分布區(qū)間的重要參數，因此準確預測不同凝固條件下一次枝晶臂間距的上下限值顯得尤為重要。有研究表明[13, 26]，定向凝固一次枝晶臂間距的選擇是一個自組織過程，在固定的溫度梯度G和凝固速度v作用下，柱狀晶一次枝晶臂間距存在一個穩(wěn)態(tài)區(qū)間λmin～λmax（見表1）。當枝晶臂間距λ＜λmin，由于相鄰枝晶臂之間強烈的溶質交互作用，會出現競爭生長，通過淹沒機制來調整枝晶臂間距至穩(wěn)態(tài)區(qū)間；當枝晶臂間距λ＞λmax，枝晶間的溶質成分富集較低，過冷度大，有利于通過側向分叉機制產生新的枝晶臂來降低枝晶臂間距。為了預測穩(wěn)態(tài)枝晶臂間距的上下限值，在固定的計算區(qū)域（2.6mm×5mm，網格尺寸5μm）中通過設置不同的初始晶粒密度來調整初始一次枝晶臂間距，同時采用與圖1中實驗一致的凝固條件進行模擬。

圖3所示為初始枝晶臂間距在λmax～λmin區(qū)間的枝晶生長情況。以圖3（b1）～（b4）為例，當初始枝晶臂間距λini=170μm時（見圖3（b1）），初始一次枝晶臂間距發(fā)生淹沒現象，淘汰了部分的一次枝晶臂，使得一次枝晶臂間距增大。而當λini超過180μm時（見圖3（b2）），穩(wěn)態(tài)一次枝晶臂間距等于初始值λini，表明在凝固速度v=100μm/s時，最小一次枝晶臂間距λmin約為180μm。當λini=575μm時（見圖3（b3）），凝固過程中二次枝晶臂上形成大量的三次枝晶臂，但是由于相鄰枝晶臂空間上的限制和自身二次枝晶臂的阻礙，這些三次枝晶臂很難發(fā)展成為穩(wěn)定的一次枝晶臂。當λini增大到590μm 時（見圖3（b4）），部分三次枝晶臂發(fā)展成了一次枝晶臂，從而降低了相鄰枝晶臂間距，這表明上限值λmax處于575～590μm。

穩(wěn)態(tài)一次枝晶臂間距的下限值λmin也可以通過Hunt-Lu模型來預測，即式（2）。該模型并沒有給出上限值的函數表達式，而是認為上限值λmax約為下限值λmin的2倍。圖4所示為模擬得到的一次枝晶臂間距上下限值、表1中實驗值以及Hunt-Lu模型預測結果（λmax=2λmin）之間的比較，從圖4中可以看到，CA模型模擬得到的結果與實驗結果非常接近，吻合程度明顯優(yōu)于Hunt-Lu模型的預測結果。DING等[26]基于透明合金丁二腈-乙醇溶液定向凝固實驗結果總結出一次枝晶臂上限值λmax和下限值λmin相對于凝固速度v呈冪函數關系式：λmax=avb、λmin=cvd。將CA模擬結果進行數據擬合，得到的方程為：λmax=8837v-0.61、λmin=3495v-0.65，可以看到上限值λmax和下限值λmin基

本呈3倍的關系，與實驗結果吻合良好。

圖3　不同凝固速度和不同初始一次枝晶臂間距條件下定向生長枝晶陣列Fig.3　Steady state dendrite arrays in directional solidification at different initial primary dendrite arm spacings （λini）and different growth velocities：（a1）-（a4）50μm/s；（b1）-（b4）100μm/s；（c1）-（c4）150μm/s；（d1）-（d4）200μm/s

圖4　不同凝固速度下一次枝晶臂間距上下限λmax、λmin預測結果和實驗結果比較Fig.4　Comparison of allowable range of primary dendrite spacing between CA simulations, analytical predictions and experimental results （λmaxpredicted by Hunt-Lu model is twice of λmin；λmax,simand λmin, simpredicted by CA model）

2.2不同凝固速度下枝晶定向生長過程模擬

圖5　定向凝固條件下模擬得到的不同凝固速度下枝晶穩(wěn)態(tài)生長時的枝晶陣列Fig.5　Simulated steady state dendrite arrays at different growth velocities in directional solidification：（a）200μm/s；（b）150μm/s；（c）100μm/s；（d）50μm/s

文獻[26]指出，一次枝晶臂間距的選擇不僅取決于穩(wěn)態(tài)凝固參數，而且還與枝晶生長的歷史過程相關。本研究中采用的是單一凝固條件，即枝晶在整個生長過程中采用恒定的凝固速度和溫度梯度，因此，可以認為在本研究中枝晶生長歷史對一次枝晶臂間距的影響較小。為了模擬圖1中4種凝固條件下的枝晶生長過程，選取計算域大小為6.25mm×15mm，CA單元尺寸5μm。初始在計算域底部放置一定數量的晶核。圖5所示為模擬得到的4種凝固條件下的枝晶定向凝固演化過程。從圖5中可以看到，在凝固初期，從晶核的二次枝晶臂上發(fā)展出了大量的三次枝晶臂，并且部分三次枝晶臂發(fā)展成了穩(wěn)定的一次枝晶臂。這主要是因為初始的晶核間距λini遠大于一次枝晶臂間距上限值λmax，從而出現側向分枝而形成新的一次枝晶臂，但由于大部分枝晶臂相鄰間距小于一次枝晶臂間距下限值λmin，從而出現淹沒，只有少數能成為穩(wěn)定的一次枝晶臂。此外，從圖5中也可以看到，凝固不久后系統(tǒng)就可以達到穩(wěn)定狀態(tài)，若沒有外界的擾動，則后續(xù)凝固過程中枝晶可以一直穩(wěn)定生長，枝晶臂間距不變。圖6所示為CA模擬結果、解析模型預測結果以及實驗結果的比較。從圖6中可以看到，Kurz-Fisher模型預測的結果普遍高于實驗結果，而Hunt模型預測的結果則小于實驗結果，并且兩個模型的預測結果與實驗值偏離的程度都較大。Trivedi模型相對較為準確，但是與CA模擬的結果相比，誤差還是較大，這表明CA模型在枝晶定向凝固過程枝晶形貌演化模擬和枝晶臂間距預測等方面的準確性。對于上述的解析模型，界面能各向異性的影響是沒有進行考慮的，而事實證明，界面能各向異性對于枝晶生長具有較大影響，而CA模型可以很好地把界面能各向異性的影響耦合進去，這可能是CA模型預測準確性較高的原因之一。此外，Hunt等模型認為一次枝晶臂均勻分布，顯然與實際不相符，這也可能導致預測結果不準確。類似的可以對CA模擬的結果進行數據擬合，得到的關系式為λave=4788v-0.59。

圖6　不同凝固速度下平均一次枝晶臂間距λave預測結果和實驗結果比較Fig.6　Comparison of analytical predictions and experimental results of average primary dendrite arm spacing λavebetween CA simulations

2.3影響一次枝晶臂間距因素分析

對于固定成分的合金，定向凝固枝晶臂間距的選擇過程受外因和內因的共同作用。外因主要是抽拉速度v、溫度梯度G等可控因素，而內因包括界面能大小、溶質擴散系數DL、枝晶生長方向與熱流的偏離角度θ等。外應和內因的共同作用會改變溫度場和溶質場，使枝晶凝固前沿的過冷度ΔT發(fā)生變化，改變枝晶尖端生長速度 （νn∝ΔT2）和界面穩(wěn)定性，從而影響枝晶的生長和一次枝晶臂的選擇過程。

2.3.1界面能和擴散系數DL的影響

界面能大小是影響枝晶尖端生長動力學的一個重要因素，體現了枝晶的分枝能力。LANGER等[27]通過枝晶尖端的線性分析，提出了尖端界面穩(wěn)定性理論，即LMK理論。該理論認為枝晶尖端半徑r和尖端生長速度v之間滿足關系：vr2=2ΓDL/（σ*ΔT0）。式中：Γ為表征界面能的Gibbs-Thomson系數；DL為液相擴散系數；ΔT0為結晶溫度區(qū)間；σ*是尖端穩(wěn)定性參數，約為1/（4π2）。從該模型中可以看到，Γ越小，會促使尖端半徑減小，當尖端半徑減小到不能保證界面穩(wěn)定時，就會通過界面分叉的形式來保證生長過程中的界面穩(wěn)定性，從而減小一次枝晶臂間距。枝晶凝固過程中，起到穩(wěn)定固液界面作用的界面能效應通過Γζ來體現，其中ζ為表征界面能各向異性的函數，表示為ζ=1-15εcos[4（φ-θ）]，ε為界面能各向異性系數。隨著界面各向異性系數ε的增大，ζ將減小，導致界面能穩(wěn)定作用減弱，造成枝晶臂容易出現側向分枝，減小一次枝晶臂間距[28-29]。此外，由文獻[30-31]中可以看出，隨著各向異性系數ε增大，尖端穩(wěn)定性參數σ*呈冪指數增大，這也有利于枝晶臂的分枝形成。YIN 等[32]通過CA方法研究了Gibbs-Thomson系數Γ和界面能各向異性系數ε對于枝晶形貌的影響，模擬得到的結果與上述的解釋相一致。從LMK理論還可以看到，溶質擴散系數DL減小也會降低枝晶界面穩(wěn)定性，從而影響枝晶生長。DL越小，溶質擴散能力減弱，相鄰枝晶間的溶質場交疊作用降低，促進固液界面前沿的過冷度（ΔT=Tliq（CL）-T*，式中：Tliq為成分為CL的液相線溫度；T*為實際界面溫度）增大，降低枝晶界面穩(wěn)定性，有利于枝晶臂的分枝，降低枝晶臂間距以保證穩(wěn)定生長。表2所列為兩組不同液相溶質擴散系數條件下的平均一次枝晶臂間距的模擬結果。綜上所述，界面能和溶質擴散系數的減小都有利于定向凝固過程一次枝晶臂間距的減小。

表2　液相溶質擴散系數DL對一次枝晶臂間距的影響Table 2　Effects of diffusion coefficients on primary dendrite spacing during directional solidification

2.3.2凝固參數（v,G）的影響

定向凝固一次枝晶臂間距主要取決于溫度梯度G和生長速度v這兩個可控的外部因素，許多研究者[24,26]針對這兩個參數對定向凝固一次枝晶臂間距的影響進行了大量實驗研究，表明一次枝晶臂間距隨著溫度梯度和生長速度的增大而減小，這一結論與式（1）相一致。溫度梯度和生長速度對枝晶臂間距的影響機制類似，兩者的變化會改變固液界面前沿的散熱效果和溶質場分布，從而影響界面前沿的過冷度ΔT，導致枝晶臂間距發(fā)生變化。為了說明這一點，圖7所示為初始枝晶臂間距λini=320μm、溫度梯度G=15 K/mm、凝固速度v為100和200μm/s兩種條件下凝固前沿的過冷度分布。從圖7中可以看出，由于相鄰枝晶臂間的“V”形區(qū)域落后于枝晶尖端，使得“V”形區(qū)域的溫度值T*較低，熱過冷度較大，從而使得該處的局部冷度比枝晶尖端大，因此，該處容易發(fā)展出穩(wěn)定的一次枝晶臂來填充該區(qū)域（見圖3）。比較圖7中兩組不同凝固速度下的模擬結果可以看到，當凝固速度從100增加至200μm/s時，由于冷卻速度增大（CR=G·v），所以在凝固至同一高度時，在v=200μm/s條件下的溫度值T*要比在v=100μm/s條件下的低，從而枝晶凝固前沿的過冷度ΔT在v=200μm/s條件下明顯比v=100μm/s時的大。過冷度的增大促進枝晶生長速度增大，從而導致一次枝晶臂的形成和一次枝晶臂間距的減小。

圖7　初始枝晶臂間距λini=320μm、凝固速度v=100μm/s 和v=200μm/s下凝固前沿的過冷度分布Fig.7　Undercooling distributions of solidification front with same initial dendrite spacing （λini）of 320μm at different growth velocities：（a）v=100μm/s；（b）v=200μm/s

2.3.3生長方向偏離熱流方向角度θ的影響

實際定向凝固過程中，枝晶生長方向并不完全平行于熱流方向，兩者往往存在一定角度θ，如圖8所示。從圖8中可以看到，晶粒A和晶粒B的〈100〉方向并不完全沿著溫度梯度G的方向，而是存在一定的夾角，且晶粒B的夾角θB（15°）要大于θA（2°）。通過測定A和B晶粒的一次枝晶臂間距發(fā)現，B晶粒的平均一次枝晶臂間距（234μm）要大于A晶粒的平均一次枝晶間距（206μm），這與文獻[33-34]中的實驗結果是相一致的，即枝晶臂間距會隨著枝晶〈100〉生長方向與溫度梯度方向之間夾角θ的增大而增大。從式（1）和（2）可以看到，HUNT等[5]提出的預測枝晶臂間距的理論模型中并沒有考慮枝晶生長取向對枝晶間距的影響，因此，需要進行修正才能解釋晶體取向對一次枝晶臂間距的影響。若一次枝晶臂的生長方向與熱流方向的夾角為θ，則枝晶沿〈100〉方向的生長速度和溫度梯度分別為v/cosθ、Gcosθ，代入式（1）可得

根據式（3）可以看到，對于固定成分C0，在恒定溫度梯度G和生長速度v作用下，當θ增大時，cosθ減小，從而使得一次臂間距λ1增大。GANDIN等[35]在研究晶體取向對枝晶生長的影響時，指出一次枝晶臂間距對θ的依賴關系主要是由三次枝晶臂和二次枝晶臂相互間的競爭機制決定的，并基于分枝機制建立了一次枝晶臂間距預測模型：

式中：λ（C0,v,G）、e、d均為正數，從而也可以說明枝晶臂間距隨著擇優(yōu)取向偏離角度θ的增大而增大。

圖8　生長方向偏離熱流方向的角度θ大小對枝晶生長和枝晶臂間距的影響（G=15 K/mm, v=200μm/s）Fig.8　Effects of deviation angle θ between dendrite growth direction and heat flux direction on dendrite growth and dendrite spacing （G=15 K/mm, v=200μm/s）

通過模擬研究了θ對一次枝晶臂間距的影響。初始在計算域底部設置不同生長取向的晶核（5°～30°），且G=15 K/mm、v=200μm/s，并測定穩(wěn)態(tài)生長時的枝晶間距。圖9所示為平均枝晶臂間距隨θ的變化曲線，可以看到λave隨著θ的增大而增大，這與式（3）和（4）的解釋相一致。

圖9　平均一次枝晶臂間距隨生長取向角度θ的變化曲線（G=15 K/mm, v=200μm/s）Fig.9　Predicted dendrite arm spacing varying with deviation angle θ （G=15 K/mm, v=200μm/s）

3　結論

1）通過開展Al-7Si-0.36Mg合金在不同凝固速度下的定向凝固實驗，測得了一次枝晶臂間距最大值λmax、平均值λave、最小值λmin以及λmax/λmin比值。

2）采用CA模型預測了Al-7Si-0.36Mg合金不同凝固速度下穩(wěn)態(tài)一次枝晶臂間距的上限值、下限值及平均值，分別可以表示為：λmax=8837v-0.61、λmin=3495v-0.65、λave=4788v-0.59。模擬結果與實驗結果的吻合程度明顯優(yōu)于Hunt-Lu等模型的預測結果，表明本模型在枝晶定向凝固過程枝晶形貌演化模擬和枝晶臂間距預測等方面的準確性。

3）結合模擬研究和文獻調研分析了影響定向凝固一次枝晶臂選擇的因素，包括抽拉速度v、溫度梯度G、界面能大小、溶質擴散系數DL、枝晶生長方向與熱流的偏離角度θ等。結果表明，界面能和溶質擴散系數的減小都有利于定向凝固過程一次枝晶臂間距的減小；溫度梯度G和生長速度v對一次枝晶臂間距的影響主要通過影響凝固前沿的過冷度來體現；枝晶生長取向偏離熱流方向越大，一次枝晶臂間距越大。

REFERENCES

[1] BOETTINGER W J, CORELL S R, GREERA L, KARAM A, KURA W, RAPPAZ M, TRIVEDI R.Solidification microstructures：Recent developments, further directions[J].Acta Materialia, 2000, 48：43-70.

[2] MELO M L N M, RIZZO E M S, SANTOS R G.Predicting dendrite arm spacing and their effect on microporosity formation in directionally solidified Al-Cu alloy[J].Journal of Materials Science, 2005, 40：1599-1609.

[3] ROCHA O L, SIQUEIRA C A, GARCIA A.Heat flow parameters affecting dendrite spacing during unsteady-state solidification of Sn-Pb and Al-Cu alloy[J].Metallurgical and Materials Transactions A, 2003, 34：995-1006.

[4] QUAREMA J M V, SANTOS C A, GARCIA A.Correlation between unsteady-state solidification condition, dendrite spacing and mechanical properties of Al-Cu alloys[J].Metallurgical and Materials Transactions A, 2000, 31：3167-3178.

[5] HUNT J D.Solidification and casting of metals[M].London：The Metal Society, 1979.

[6] KURZ W, FISHER D J.Dendrite growth at the limit of stability：tip radius and spacing[J].Acta Metallurgica, 1981, 29：11-20.

[7] TRIVEDI R.Interdendritic spacing：partⅡ.A comparison of theory and experiment[J].Metallurgical Transactions A, 1984, 15：977-982.

[8] WARREN J A, LANGER J S.Stability of dendritic spacings[J].Physics Review A, 1990, 42：3518-3525.

[9] WARREN J A, LANGER J S.Prediction of dendritic spacings in a directional solidification experiment[J].Physics Review E, 1993, 47：2702-2712.

[10] MA D.Modeling of primary spacing selection in dendrite arrays during directional solidification[J].Metallurgical and Materials Transactions B, 2002, 33：223-233.

[11] HUANG W D, GENG X G, ZHOU Y H.Primary spacing selection of constrained dendritic growth[J].Journal of Crystal Growth, 1993, 134：105-115.

[12] HUNT J D, LU S Z.Modeling of cellular-dendritic array growth：spacing and structure predictions[J].Metallurgical and Materials Transactions A, 1996, 27：611-623.

[13] WANG W, LEE P D, MELEAN M.A model of solidification microstructures in nickel based superalloys：Predicting primary dendrite spacing selection[J].Acta Materialia, 2003, 51：2971-2987.

[14] DIEPERS H J, MA D, STEINBACH I.History effects during the selection of primary dendrite spacing：Comparison of phase-field simulation with experimental observations[J].Journal of Crystal Growth, 2002, 237/239：149-153.

[15] 單博煒, 黃衛(wèi)東, 林鑫, 魏雷.元胞自動機模型模擬枝晶一次間距的選擇[J].金屬學報, 2008, 44：1042-1050.SHAN Bo-wei, HUANG Wei-dong, LIN Xin, WEI Lei.Dendrite primary spacing selection simulation by the cellular automaton model[J].Acta Metallurgica Sinica, 2008, 44：1042-1050.

[16] ZHANG Xian-fei, ZHAO Jiu-zhou.Dendritic microstructure formation in a directionally solidified Al-11.6Cu-0.85Mgalloy[J].Journal of Crystal Growth, 2014, 391：52-58.

[17] 張航, 許慶彥, 史振學, 柳百成.DD6高溫合金定向凝固枝晶生長的數值模擬研[J].金屬學報, 2014, 50：345-354.ZHANG Hang, XU Qing-yan, SHI Zhen-xue, LIU Bai-cheng.Numerical simulation of dendrite grain growth of DD6 superalloy during directional solidification process[J].Acta Metallurgica Sinica, 2014, 50：345-354.

[18] 石玉峰, 許慶彥, 龔銘, 柳白成.定向凝固過程中NH4Cl-H2O 枝晶生長的數值模[J].金屬學報, 2011, 47：620-627.SHI Yu-feng, XU Qing-yan, GONG Ming, LIU Bai-cheng.Simulation of NH4Cl-H2O dendritic growth in directional solidification[J].Acta Metallurgica Sinica, 2011, 47：620-627.

[19] GANDIN C A, RAPPAZ M.A 3D cellular automaton algorithm for the prediction of dendritic grain growth[J].Acta Materials, 1997, 45（5）：2187-2195.

[20] CHEN Rui, XU Qing-yan, LIU Bai-cheng.Cellular automaton simulation of three-dimensional dendrite growth in Al-7Si-Mg ternary aluminum alloys[J].Computational Materials Science, 2015, 105：90-100.

[21] CHEN Rui, XU Qing-yan, LIU Bai-cheng.A modified cellular automaton model for the quantitative prediction of equiaxed and columnar dendritic growth[J].Journal of Materials Science and Technology, 2014, 30：1311-1320.

[22] 石玉峰.鋁合金砂型鑄件微觀組織及力學性能的數值模擬[D].北京：清華大學, 2013.SHI Yu-feng.Numerical simulation of microstructure and mechanical properties of aluminum alloy sand casting[D].Beijing：Tsinghua University, 2013.

[23] QUESTED P N, MCLEAN M.Solidification morphologies in directionally solidified superalloys[J].Material Science and Engineering, 1984, 65：171-180.

[24] 黃太文, 劉林, 張衛(wèi)國, 張軍, 傅恒志.抽拉速度躍遷對定向凝固單晶高溫合金 DD3一次枝晶間距和微觀偏析的影響[J].金屬學報, 2009, 45（10）：1225-1231.HUANG Tai-wen, LIU Lin, ZHANG Wei-guo, ZHANG Jun, FU Heng-zhi.The primary dendrite arm spacing and microsegregation of directionally solidified single crystal superalloy DD3[J].Acta Metallurgica Sinica, 2009, 45（10）：1225-1231.

[25] HAN S H, TRIVEDI R.Primary spacing selection in directionally solidified alloys[J].Acta Metallurgica and Materialia, 1994, 42：25-41.

[26] DING G L, HUANG W D, HUANG X, LIN X, ZHOU Y H.On primary dendritic spacing during unidirectional solidification[J].Acta Metallurgica, 1996, 44：3705-3709.

[27] LANGER J S, MüLLER-KRUMBHAAR H.Theory of dendritic growth-Ⅰ.Elements of a stability analysis[J].Acta Metallurgica, 1978, 26：1681-1687.

[28] 張顯飛, 趙九州.合金定向凝固一次枝晶間距模擬[J].中國有色金屬學報, 2012, 22（10）：2868-2874.ZHANG Xian-fei, ZHAO Jiu-zhou.Simulation of primary dendrite spacing in unidirectionally solidified alloy[J].The Chinese Journal of Nonferrous Metals, 2012, 22（10）：2868-2874.

[29] 張云鵬, 林鑫, 魏雷, 王猛, 彭東劍, 黃衛(wèi)東.用CA方法模擬界面能各向異性對胞晶生長形態(tài)的影響[J].物理學報, 2012, 61：228106.ZHANG Yun-peng, LIN Xin, WEI Lei, WANG Meng, PENG Dong-jian, HUANG Wei-dong.Effect of surface tension anisotropy on the growth patterns of cellulars in directional solidification[J].Acta Physics Sinica, 2012, 61：228106.

[30] KESSLER D A, LEVINE H.Velocity selection in dendritic growth[J].Physics Review B, 1986, 33：7867-7870.

[31] XU J J.Generalized needle solutions, interfacial instabilities and pattern formations[J].Physics Review E, 1996, 53：5051-5062.

[32] YIN H, FELICELLI S D, WANG L.Simulation of a dendrite microstructure with the lattice Boltzmann and cellular automaton methods[J].Acta Materialia, 2011, 59：3124-3136.

[33] 王賢斌, 林鑫, 王理林, 白貝貝, 王猛, 黃衛(wèi)東.晶體取向對定向凝固枝晶生長的影響[J].物理學報, 2013, 62（10）：108103.WANG Xian-bin, LIN Xin, WANG Li-lin, BAI Bei-bei, WANG-Meng, HUANG Wei-dong.Effect of crystallographic orientation on dendrite growth in directional solidification[J].Acta Physics Sinica, 2013, 62（10）：108103

[34] 楊初斌, 劉林, 趙新寶, 劉剛, 張軍, 傅恒志.〈001〉和〈011〉取向DD407單晶高溫合金枝晶間距和微觀偏析[J].金屬學報, 2011, 47：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 Metallurgica Sinica, 2011, 47：1246-1250.

[35] GANDIN C H A, ESHELMAN M, TRIVEDI R.Orientation dependence of primary dendrite spacing[J].Metallurgical and Materials Transactions A, 1996, 27：2727-2739.

（編輯王超）

Experimental and simulation of primary dendrite spacing in directional solidification of Al-7Si-0.36Mg alloy

CHEN Rui, XU Qing-yan, LIU Bai-cheng
（Key Laboratory for Advanced Materials Processing Technology, Ministry of Education, School of Materials Science and Engineering, Tsinghua University, Beijing 100084, China）

The directional solidification dendritic morphology evolution and primary dendrite arm spacing selection process of Al-7Si-0.36Mg alloy were studied by directional solidification experiments and cellular automaton model.The results indicate that there is a continuous range of primary dendrite arm spacing under the given solidification condition.Under the condition of the given temperature gradient and various solidification velocities, the relationships between the primary dendrite arm spacing parameters （λmax, λave, λmin）of Al-7Si-0.36Mg alloy and growth velocity were expressed, and λmax/λmin≈3.The simulated results show a quite good agreement with the experimental results, which is better than predicted results of Hunt-Lu model.The comparisons reveal that the present CA model has a high accuracy in simulating the evolution of dendrite morphology and predicting primary dendrite arm spacing in directional solidification.Based on the predictions and related literatures, the factors influencing the selection of primary dendrite arm spacing, such as growth velocity v, temperature gradient G, interfacial energy effect, solute diffusion coefficient DL, as well as the deviation angle θ between the dendrite growth direction, and heat flux direction were analyzed.

directional solidification；primary dendrite spacing；growth velocity；cellular automaton；numerical simulation

TG290

A

1004-0609（2015）10-2613-10

國家重點基礎研究發(fā)展計劃項目（2011CB706801）；國家自然科學基金資助項目（51374137, 51171089）

2015-01-12；

2015-05-16

許慶彥，教授，博士；電話：010-62795482；E-mail：scjxqy@mail.tsinghua.edu.cn