中国有色金属学报 2014,24(09),2295-2301 DOI:10.19476/j.ysxb.1004.0609.2014.09.014
冷却速率对Cu70Zr30合金初生相及包晶转变的影响
杨伟 张燕龙 余欢 陈寿辉 卢百平 汪志太
1. 南昌航空大学轻合金加工科学与技术国防重点学科实验室
3. 西北工业大学凝固技术国家重点实验室
摘 要:
采用差热分析与阶梯铜模喷铸技术, 研究不同冷却速率作用下Cu70Zr30合金初生相尺寸、形貌以及包晶层厚度等的演化规律, 讨论相选择与包晶转变的形成机理。结果表明:近平衡凝固条件下, 差热分析试样中初生相形成与包晶转变的发生温度分别为1042℃和957℃, 均滞后于平衡相图中相应的温度。由于包晶相两侧成分区间相差不大, 且固态中原子扩散系数较小, 因此, 平均包晶层厚度仅为5μm, 并存在残余的Cu51Zr14初生相。铜模喷铸条件下, 非平衡凝固组织中初生相结构未发生改变, 但包晶转变得到有效抑制。随冷却速率提高, 过冷度的增加有利于形核发生, 初生相形貌从定向束状向等轴晶发生转变。
关键词:
Cu-Zr合金;冷却速率;非平衡凝固;相选择;包晶转变;
中图分类号: TG146.11
作者简介:杨伟, 博士;电话:0791-86453167;传真:0791-86453167, E-mail:weiyang@mail.nwpu.edu.cn;
收稿日期:2013-12-02
基金:国家自然科学基金资助项目 (51164028);国家高技术研究发展计划资助项目 (2013AA064001);江西省教育厅资助项目 (GJJ14504);轻合金加工科学与技术国防重点学科实验室开放课题 (GF201201003);
Effects of cooling rate on primary phase and peritectic transformation in Cu70Zr30 alloy
Abstract:
The effects of cooling rate on the size and morphology of primary phase in Cu70Zr30 alloy, as well as the thickness of peritectic layer were investigated by differential thermal analysis and step copper mould spray casting techniques. The formation mechanisms of phase selection and peritectic transformation were discussed. The results show that the measured temperatures for primary phase and peritectic transformation are respectively 1042 ℃ and 957 ℃ in near-equilibrium solidification, which are smaller than the values from equilibrium phase diagram. Due to slight composition interval between the adjacent regions of peritectic phase and sluggish solute diffusion rate in solid, the average thickness of peritectic layer is merely 5 μm with the presence of residual primary phase Cu51Zr14. The structure of primary phase is not changed by copper mould spray casting, while peritectic transformation is suppressed effectively. With increasing cooling rate, undercooling is enhanced to activate more nucleation events, which generates the morphology transition of primary phase from directionally bundle to equiaxed-grain microstructure.
Keyword:
Cu-Zr alloy; cooling rate; non-equilibrium solidification; phase selection; peritectic transformation;
Received: 2013-12-02
Cu-Zr合金作为非晶形成能力较强的二元合金, 具有强度高、硬度高、耐腐蚀以及成本低廉等诸多优点[1,2]。通过添加合适的元素还可以进一步提高其非晶形成能力与综合性能, 因此, 成为极具应用前景的非晶结构材料[3,4,5]。XU等[6]采用铜模吸铸方式, 研究了Cu100-xZrx (x=34~40) 非晶合金的临界尺寸, 指出Cu64Zr36的非晶形成能力最强。WANG等[7]采用Cu50Zr50合金制备出直径为2 mm的非晶棒。ABE等[8]采用热力学计算手段, 讨论了初生相种类对相变驱动力的影响。然而, 现有研究大多集中于非晶形成能力的预测、热力学数据的优化以及晶化动力学的描述方面, 有关冷却速率对该合金系晶化组织影响规律的研究仍不系统。
由相图[2]可知, 某些成分的Cu100-xZrx (x=10~38) 合金凝固阶段存在着诸多复杂相变过程, 其中包括初生相的形成、包晶及共晶等多相凝固[2]。由于相组成和微观形貌对材料性能具有重要影响, 因此, 为最大限度地开发合金应用潜力, 准确理解并有效控制凝固过程中各相及组织的形成尤为重要。
非平衡凝固作为材料科学领域中的研究热点, 在改善传统材料性能及合成新型亚稳相材料方面已发挥出巨大潜力[9,10,11,12], 其中所形成的复杂组织结构既涉及热力学和动力学耦合作用下的相选择[13], 又涉及形核和生长过程间的相互竞争[14]。本文作者针对Cu70Zr30合金, 采用差热分析与铜模喷铸技术, 对比研究不同冷却速率作用下合金的相组成及晶粒尺寸、形貌、取向等微观组织的演化规律, 从而为理解并有效控制该类合金凝固过程提供理论参考。
1 实验
采用高真空电弧熔炼炉制备出Cu70Zr30合金, 其中Cu和Zr的纯度分别为99.99%和99.9%。将称量好的合金原料 (约80 g) 置于水冷铜坩埚内, 抽真空至炉内压力达6.67×10-4Pa后, 返充99.999%高纯氩气至0.03 MPa。为保证合金成分均匀性, 熔炼过程中在坩埚下方进行电磁搅拌并翻转, 反复熔炼3次。
差热分析在德国NETZSCH STA449F3同步热分析仪上进行。首先采用线切割方式从合金锭中取出直径为4 mm、厚度为5 mm的圆片, 经砂纸打磨及超声波清洗去除表面油污后, 制备出质量约120 mg薄片试样。实验在动态高纯氩气下进行, 其中温度范围为25~1200℃, 升温速率为20 K/min。为获得近平衡凝固组织, 冷却速率选择为5 K/min。为消除热滞现象, 仪器所测温度与In、Sn、Zn、Al、Ag和Au等高纯金属的熔点进行校正, 确保温度精度为±1 K。
非平衡凝固实验在高真空感应熔炼炉内进行。首先将约20 g合金小块装入石英坩埚中 (尺寸d15 mm×150 mm, 下方开有直径为0.5 mm的小孔) 。采用机械泵+分子泵二级泵抽气系统实现6.67×10-4Pa极限真空, 并返充高纯氩气至0.03 MPa。通过改变线圈输出功率, 最高熔炼温度控制在1300℃, 并采用高精度红外测温仪进行实时测温。合金熔化并保温10min确保均匀混合后, 在坩埚上方通入0.02 MPa氩气, 金属液垂直喷射入不同内径的阶梯铜模内 (外径80mm, 内径由下向上依次为4、6、8和10 mm) , 制备出不同冷却速率下的亚快速凝固试样。
试样经镶嵌、打磨、抛光及腐蚀 (40%氢氟酸+60%乙醇腐蚀液, 质量分数) 等工序后, 在VHX-600E三维超景深显微镜上进行低倍光学显微组织分析。采用Quanta 200扫描电镜进行高倍组织分析, 利用INCA250 X-Max 50型能谱仪对微区成分进行测定。采用D8-X射线衍射仪进行物相结构分析。
2 实验结果
2.1 差热分析曲线
根据相图可知[2], 平衡条件下Cu70Zr30合金相变过程是:当熔体温度1061℃时, 首先从液相中析出Cu51Zr14初生相;当温度降到包晶温度TP (975℃) 时, 发生包晶反应:L+Cu51Zr14→Cu8Zr3;当温度继续降到885℃时, 发生共晶反应L→Cu8Zr3+Cu10Zr7。凝固结束后合金平衡组织由包晶相 (Cu8Zr3) 和共晶相 (Cu8Zr3+Cu10Zr7) 组成。然而, 由于包晶相的形成需要通过固相原子扩散来实现, 因此, 常规铸态组织中包晶转变并不彻底。
图1所示为Cu70Zr30合金的差热分析曲线。由图1可知, 其升温过程中出现3个主要吸热峰, 峰值温度分别为912, 941和1111℃。这与实际相图偏差较大, 可解释为在电弧熔炼过程中水冷铜坩埚强烈的激冷效果导致合金内部形成复杂的相。在5 K/min的冷却阶段, 凝固放热峰减少为两个, 分别为1042℃时的初生相凝固与957℃时的包晶转变。由于存在过冷度及包晶反应滞后现象, 因此, 所测凝固发生温度与包晶转变温度均低于平衡相图中相应的温度[15]。此外, 由于包晶转变结束后残余液相体积分数较小, 平衡相图中的共晶反应未在差热分析曲线中出现。
![]()
图1 加热速率为20 K/min、冷却速度为5 K/min时Cu70Zr30合金的差热分析曲线Fig.1 Differential thermal analysis curves of Cu70Zr30alloy at heating rate of 20 K/min and cooling rate of 5 K/min
2.2 冷却速率对Cu70Zr30合金组织形貌的影响
图2所示为Cu70Zr30合金在冷却速率为5 K/min时的近平衡凝固以及铜模内径为10、6和4 mm时的亚快速凝固光学组织。由图2 (a) 可知, 差热分析试样中存在着明显的包晶转变层, 厚度大约为5μm。即使在5 K/min缓冷条件下, 包晶转变进行得仍不彻底, 可清楚地看到包晶层内部所包裹的残余Cu51Zr14初生相, 并且晶界区域存在部分共晶组织。相比图2 (a) 近平衡凝固组织, 该合金的非平衡凝固组织呈现出较为复杂的变化 (见图2 (b) ~ (d) ) 。铜模喷铸条件下, 包晶转变得到有效抑制, 即使在d 10 mm试样中包晶层也不明显, 并且在靠近铜模内壁处, 初生相以定向束状生长为主 (见图2 (b) ) 。随铜模内径的减小, 合金冷却速率不断提高, 初生相形貌呈现出束状和等轴晶混合组织 (见图2 (c) ) , 并且体积分数不断增加。当铜模内径为4 mm时, 非平衡凝固组织由完全非定向等轴晶组成, 尤其在试样边缘部位, 初生相晶粒数目显著增多 (见图2 (d) ) 。
为确定各相成分, 图3给出了差热分析试样与d 4 mm试样的扫描电镜组织与EDS能谱分析结果。冷却速率5 K/min时, 包晶相呈现非小平面特征 (见图3 (a) ) 。铜模喷铸条件下, 由于包晶转变受到抑制, 初生相的原始小平面特征得以保存, 同时该相周围存在一定量的共晶组织 (见图3 (b) 中E点) , 说明冷却速率的提高不能抑制共晶反应的发生。EDS能谱分析 (见图3 (c) ~ (e) ) 表明:近平衡条件下, 沿包晶层外侧向中心移动, Zr元素含量不断降低, 由19.25% (C点) 下降到15.81% (A点) , 其中包晶层Zr含量为16.28% (B点) , 从而说明包晶转变中的固态原子扩散特征。非平衡条件下, 初生相中Zr含量为15.56% (D点) , 与A点的基本接近, 由此可判断初生相结构未发生改变。
![]()
图2 不同冷却条件下Cu70Zr30合金光学显微组织Fig.2 Optical microstructures of Cu70Zr30alloy under different cooling conditions: (a) Cooling rate of 5 K/min; (b) Copper mould spray casting with inner diameter of 10 mm; (c) Copper mould spray casting with inner diameter of 6 mm; (d) Copper mould spray casting with inner diameter of 4 mm
![]()
图3 不同冷却条件下Cu70Zr30合金SEM像及EDS能谱分析结果Fig.3 SEM images and EDS spectra of Cu70Zr30alloy under different cooling conditions: (a) Cooling rate of 5 K/min; (b) Copper mould spray casting with inner diameter of 4 mm; (c) - (g) EDS spectra
![]()
图3 不同冷却条件下Cu70Zr30合金SEM像及EDS能谱分析结果Fig.3 SEM images and EDS spectra of Cu70Zr30alloy under different cooling conditions: (a) Cooling rate of 5 K/min; (b) Copper mould spray casting with inner diameter of 4 mm; (c) - (g) EDS spectra
2.3 冷却速率对Cu70Zr30合金相结构影响
图4所示为差热分析试样与d 4 mm Cu70Zr30的XRD谱。其中, 差热分析试样中主要相组成为Cu51Zr14、Cu8Zr3和Cu10Zr7, 这与图2 (a) 和3 (a) 组织分析结果一致。此外, 铜模喷铸组织中没有新相形成, 主要相仍为Cu51Zr14相, 这与图3 (c) 中能谱分析结果一致。然而由于包晶转变的抑制, 亚快速凝固试样中包晶相Cu8Zr3不再存在, 因此, 初生相与共晶相衍射峰数目增多。
![]()
图4 不同冷却条件下Cu70Zr30合金的XRD谱Fig.4 XRD patterns of Cu70Zr30alloy under different cooling conditions
3 分析与讨论
3.1 包晶转变层分析
基于溶质扩散的包晶转变过程, St JOHN等[16]提出, 当低于包晶温度时, 包晶层厚度Δx随时间t的变化关系为
![]()
式中:γ是与合金相图有关的因子;C1和C2分别为包晶层两侧溶质含量;D为包晶相中的平均互扩散系数。
由式 (1) 可知, 包晶相两侧成分区间越大, 包晶转变进行得越彻底。然而, 由于本研究中包晶相Cu8Zr3为固定化学计量比, 并且固相中溶质扩散系数非常小, 因此, 包晶层厚度Δx随时间增加极为缓慢。5 K/min缓冷条件下, Δx仅为5μm左右。随冷却速率提高, 包晶转变易被抑制。
3.2 初生相结构选择
假设液态金属冷却过程中出现一个稳定晶核时即开始凝固, 该过程可表示为
![]()
式中:TL为液相线温度;TN为形核温度;I为形核速率;V为熔体体积;Rc为冷却速率。
铜模喷铸过程中Rc的计算式为[17,18]
![]()
式中:T0为铜模温度;ρ为合金密度;cp为合金的比热容;d为铜模内径;α为换热系数。
由式 (2) ~ (3) 可知, 当d减小时, 冷却速率Rc增加。在I、Vm及TL不变的条件下, 形核温度TN下降, 即冷却速率越大, 过冷度ΔT越大。因此, 为研究相选择及形核过程, 需要定量分析不同相形成时的热力学驱动力变化。
Cu-Zr合金系具有较大的负混合焓, 相图中存在6种化合物 (Cu9Zr2, Cu51Zr14, Cu8Zr3, Cu10Zr7, Cu Zr和Cu Zr2) 。根据热力学自协模型, 并利用Cu-Zr热力学数据[8], 图5给出了Cu51Zr14、Cu8Zr3和Cu10Zr7析出时所对应的驱动力变化ΔG。由图5可知, 随过冷度增加, 形核驱动力不断增加。根据最大驱动力原理, 当初生相为Cu51Zr14时, 系统自由能差最大, 因此从热力学上讲Cu51Zr14相最易析出, 这与本研究结果一致, 即铜模喷铸高冷却速率条件下初生相结构未发生变化。类似结果在Cu50Zr50合金静电悬浮实验中也有报道[19], 当熔体过冷到310 K时, X射线同步辐射检测出的初生相B2结构没有改变。
![]()
图5 Cu70Zr30合金不同相形成时热力学驱动力随过冷度的变化Fig.5 Driving forces for formation of different compounds from Cu70Zr30alloy as function of undercooling degree
3.3 形核过程分析
通过上述分析可知, 冷却速率对熔体的形核及形貌取向也有影响。对于d 10 mm的试样, 形核一旦在型壁处发生, 将在熔体内部快速生长, 从而形成束状的定向组织, 此时较大的冷却速率提高了初生相的择优生长取向, 但不影响形核数目 (见图2 (b) ) 。然而, 随冷却速率提高, d 4 mm试样的非平衡凝固组织以细小等轴晶为主, 越靠近铜模处, 晶粒数目越多, 说明该冷却速率条件下形核率得到显著提高 (见图2 (d) ) 。
根据经典形核理论, 衡量形核难易程度的临界形核功 (ΔG*) 和临界形核半径 (r*) 表达式分别为[20]
![]()
式中:ΔGV为体积自由能差, 可表示为ΔGV=ΔG/Vm, Vm为摩尔体积;σ为液-固自由能, 可表示为[21]
![]()
式中:γ为晶体结构因子;ΔSf为熔化熵;NA为阿伏伽德罗常数。
合金凝固过程中非均质形核率为[20]
![]()
式中:AV为单位体积内有效形核质点面积;Ns0为与形核质点接触的液态原子数目;ν为原子振动频率;f (θ) 为接触因子。
利用合金热物性参数[2], 由式 (4) ~ (6) 可计算出不同过冷度条件下Cu51Zr14初生相形成时的ΔG*和r* (见图6 (a) ) 。分析可知, 随过冷度增加, 二者均出现不同程度的下降, 说明有利于形核。由式 (7) 可计算出合金形核率随过冷度的变化 (见图6 (b) ) 。可以看出, 当过冷度达到一定程度时, 形核率急剧增加。这与图2 (d) 中结果一致。
4 结论
1) 冷却速率为5 K/min时, Cu70Zr30合金差热分析实验中初生相凝固及包晶转变温度分别为1042和957℃, 均滞后于平衡相图, 凝固组织由Cu51Zr14、Cu8Zr3和Cu10Zr7组成。由于包晶相两侧成分相差不大, 且固态原子扩散系数较小, 因此, 包晶层厚度仅为5μm, 存在残余初生相Cu51Zr14。
2) 铜模喷铸条件下, Cu70Zr30合金非平衡凝固组织由Cu51Zr14和Cu10Zr7组成, 包晶转变得到有效抑制。相变驱动力计算结果表明, Cu51Zr14相析出时热力学驱动力最大。在所研究冷却速率范围内, 初生相结构未发生改变。
![]()
图6 Cu70Zr30合金临界形核半径与临界形核功随过冷度变化以及形核速率随过冷度变化Fig.6 Critical nucleation radius, critical nucleation work (a) and nucleation rate (b) for Cu70Zr30alloy as function of undercooling degree
3) 冷却速率的提高导致形核过冷度增加, 形核率显著提高, 初生相形貌从定向束状向等轴晶发生转变, 越靠近铜模位置处, 晶粒数目越多。
参考文献
[1]张程煜, 姚可夫.Cu64Zr36非晶合金的晶化动力学[J].稀有金属材料与工程, 2006, 35 (1) :158-160.ZHANG Cheng-yu, YAO Ke-fu.Crystallization kinetics of Cu64Zr36 amorphous alloy[J].Rare Metal Materials and Engineering, 2006, 35 (1) :158-160.
[2] YANG W, LIU F, LIU H, WANG H F, CHEN Z, YANG G C.Glass forming ability in Cu-Zr binary alloy:Effect of nucleation mode[J].Journal of Alloys and Compounds, 2009, 484:702-707.
[3]寇生中, 陈延, 刘广桥, 李娜, 郑宝超, 黄文军, 索红莉.半固态处理对Zr60Cu17.5Al7.5Ni10Ti5非晶形成能力和力学性能的影响[J].中国有色金属学报, 2011, 21 (8) :1875-1880.KOU Sheng-zhong, CHEN Yan, LIU Guang-qiao, LI Na, ZHENG Bao-chao, HUANG Wen-jun, SUO Hong-li.Effect of semi-solid processing on glass forming ability and mechanical properties of Zr60Cu17.5Al7.5Ni10Ti5 bulk amorphous alloys[J].The Chinese Journal of Nonferrous Metals, 2011, 21 (8) :1875-1880.
[4]潘冶, 纪秀林, 皮锦红, 张露.微量元素对Cu-Zr-Al块体金属玻璃形成能力及力学性能的影响[J].中国有色金属学报, 2011, 21 (3) :583-587.PAN Ye, JI Xiu-lin, PI Jing-hong, ZHANG Lu.Effects of minor alloying additions on glass forming ability and mechanical properties of Cu-Zr-Al bulk metallic glasses[J].The Chinese Journal of Nonferrous Metals, 2011, 21 (3) :583-587.
[5] DEO L P, MENDES M A B, COSTA A M S, CAMPOS NETO N D, de OLIVEIRA M F.Applying a new criterion to predict glass forming alloys in the Zr-Ni-Cu ternary system[J].Journal of Alloys and Compounds, 2013, 553:212–215.
[6] XU D H, LOHWONGWATANA B, DUAN G, JOHNSON W L, GARLAND C.Bulk metallic glass formation in binary Cu-rich alloy series-Cu100-xZrx (x=34, 36, 38.2, 40 at%) and mechanical properties of bulk Cu64Zr36 glass[J].Acta Materialia, 2004, 52 (9) :2621-2624.
[7]WANG D, LI Y, SUN B B, SUI M L, LU K, MA E.Bulk metallic glass formation in the binary Cu-Zr system[J].Applied Physics Letters, 2004, 84 (20) :4029-4031.
[8] ABE T, SHIMONO M, ODE M, ONODERA H.Thermodynamic modeling of the undercooled liquid in the Cu-Zr system[J].Acta Materialia, 2006, 54 (4) :909-915.
[9] YANG W, YU H, WANG J H, CAI C C, XU Z F, LI S, LIU F, YANG G C.Application of dendrite fragmentation to fabricate the homogeneous dispersed structure in undercooled Cu-Co immiscible alloy[J].Journal of Alloys and Compounds, 2011, 509:9675-9678.
[10]赵国际, 盛光敏, 罗军.快速凝固对Sn-6.5Zn钎料合金特性及钎料/Cu焊点力学性能的影响[J].中国有色金属学报, 2012, 22 (10) :2805-2810.ZHAO Guo-ji, SHENG Guang-min, LUO Jun.Influence of rapid solidification on characteristics of Sn-6.5Zn solder alloy and mechanical properties of solder/Cu joints[J].The Chinese Journal of Nonferrous Metals, 2012, 22 (10) :2805-2810.
[11]穆丹宁, 杨长林, 魏晓伟, 刘峰.深过冷铁钴基块体非晶合金的细晶化研究[J].金属学报, 2012, 48 (12) :1409-1414.MU Dan-ning, YANG Chang-lin, WEI Xiao-wei, LIU Feng.Research on grain refinement in bulk undercooled Fe-Co base alloys[J].Acta Metallurgica Sinica, 2012, 48 (12) :1409-1414.
[12]李圣, 刘峰, 杨根仓, 杨伟.溶质含量对快速凝固Co-Cu不混溶合金晶粒细化组织的影响[J].中国有色金属学报, 2013, 23 (4) :1012-1018.LI Sheng, LIU Feng, YANG Gen-cang, YANG Wei.Effect of solute content on grain refinement microstructure of rapid solidified Co-Cu immiscible alloy[J].The Chinese Journal of Nonferrous Metals, 2013, 23 (4) :1012-1018.
[13]HERLACH D M, GILLESSEN F, VOLKMANN T, WOLLGARTEN M, URBAN K.Phase selection in undercooled quasicrystal forming Al-Mn alloy melts[J].Physical Review B, 1992, 46:5203-5210.
[14]LIU F, CHEN Y Z, YANG G C, LU Y P, CHEN Z, ZHOU Y H.Competitions incorporated in rapid solidification of the bulk undercooled eutectic Ni78.6Si21.4 alloy[J].Journal of Materials Research, 2007, 22:2953-2963.
[15]胡小武, 李双明, 高斯峰, 刘林, 傅恒志.Pb-Bi包晶合金凝固过程的DSC差热分析[J].特种铸造及有色合金, 2010, 30 (7) :589-592.HU Xiao-wu, LI Shuang-ming, GAO Si-feng, LIU Lin, FU Heng-zhi.DSC analysis of Pb-Bi peritectic alloys during solidification process[J].Special Casting&Nonferrous Alloys, 2010, 30 (7) :589-592.
[16]STJOHN D H, HOGAN L M.A simple prediction of the rate of the peritectic transformation[J].Acta Metallurgica, 1987, 35 (1) :171-174.
[17] KONG J, YE Z T, LüF.Non-equilibrium solidification character of Zr56.2Ti13.8Nb5.0Cu6.9Ni5.6Be12.5 bulk metallic glass composites containing ductile dendrite phase[J].Journal of Alloys and Compounds, 2009, 478:202-205.
[18]杨伟, 陈寿辉, 张守银, 余欢, 严青松, 蔡长春.冷却速率对AZ91D镁合金非平衡凝固组织的影响[J].中国有色金属学报, 2014, 24 (3) :593-599.YANG Wei, CHEN Shou-hui, ZHANG Shou-yin, YU Huan, YAN Qing-song, CAI Chang-chun.Effect of cooling rate on non-equilibrium solidified microstructure of AZ91D magnesium alloy[J].The Chinese Journal of Nonferrous Metals, 2014, 24 (3) :593-599.
[19] GEGNER J, SHULESHOVA O, KOBOLD R, HOLLANDMORITTZ D, YANG F, HORNFECK W, BEDNARCIK J, HERLACH D M.In situ observation of the phase selection from the undercooled melt in Cu-Zr[J].Journal of Alloys and Compounds, 2013, 576:232-235.
[20]李建国, 胡侨丹.凝固原理[M].北京:高等教育出版社, 2010:17.LI Jian-guo, HU Qiao-dan.Principle of solidification[M].Beijing:Higher Education Press, 2010:17.
[21]SPAEPEN F.A structural model for the solid-liquid interface monatomic system[J].Acta Metallurgica, 1975, 23 (6) :729-743.