DOI:10.19476/j.ysxb.1004.0609.2002.03.011
激光熔凝的数值模拟及其在激光定向凝固中的应用
刘振侠 黄卫东 杨森
西北工业大学凝固技术国家重点实验室
西北工业大学凝固技术国家重点实验室 西安710072
摘 要:
建立了一套激光熔凝的数学模型和计算方法 , 给出了两种条件下熔池深度与扫描速度的关系曲线。计算结果与实验值吻合良好 , 验证了模型的有效性。根据计算分析了激光工艺参数与激光定向凝固组织的关系。据此通过工艺参数的匹配 , 在熔池顶部得到了与扫描方向完全一致的凝固组织
关键词:
激光熔凝 ;数值模拟 ;定向凝固 ;
中图分类号: TG665
收稿日期: 2001-07-16
基金: 国家重点基础研究发展规划项目 (G2 0 0 0 0 6 72 0 5 -3); 航空基金重点资助项目 ( 97z5 30 0 2 );
Numerical simulation of laser surface re-melting and its use in laser directional solidification
Abstract:
A model and a numerical method of laser surface re-melting were established. The relationship between the depth of laser melted pools and the laser scanning velocity was calculated for two base materials. The calculations agree with the experiments well, and these validate the model. The relationship between the laser processing parameters and laser directional solidification morphology was analyzed by the model. Based on the model, the morphology with the growth direction same with scanning direction was obtained at the top of the melting pool through the match among the processing parameters.
Keyword:
laser re-melting; numerical simulation; directional solidification;
Received: 2001-07-16
高能密度的激光束以一定的速度扫过工件表面时, 由于吸收了高强度的能量, 在激光光斑的辐照区及其附近, 金属工件表面会形成激光熔池。 在熔池的前端固态金属不断熔化, 而在熔池的后端熔化了的金属又不断凝固, 由此熔池不断向前移动。 熔池内部是熔化了的液态金属, 在浮升力以及表面张力等力的驱动下液态金属会流动。 激光功率及其分布、 光斑形状及尺寸、 扫描速度、 工件材料物性等参数的不同会导致熔池中流场分布和温度分布的变化, 而不同的流场和温度场分布会产生不同的熔池形状, 进而造成不同的熔池固-液界面上的温度梯度分布。 在熔池的后端固液界面上的温度梯度对于晶体的生长有直接的影响
[1 ,2 ,3 ]
。
晶体的生长方向与生长界面上的温度梯度方向有着重要的联系。 尤其对于胞晶的生长, 因其主要受热流方向的影响, 所以我们可以通过控制熔池后端固-液边界的温度梯度的方向, 并使合金生长处于胞晶区来达到合金的定向凝固。
实际上, 熔池的尺寸是非常小的, 根据不同的工艺参数其量级一般在1 mm2 左右
[4 ,5 ,6 ]
, 而且熔池内部的温度非常高, 一般在1 000 ℃左右, 甚至更高, 加之熔池本身位置在不断变化, 所以通过实验的方法测量熔池内的流场和温度场是非常困难的, 甚至是不可能的。 因此激光熔凝的数值模拟在激光定向凝固的研究中具有特别重要的意义。
本文作者采用焓法对激光熔池内的液态金属和熔池外的固态金属建立了统一的控制方程, 用SIMPLE方法进行求解, 得到了熔池形状的有关参数。 通过与实验测得的熔池尺寸进行对比验证了模型的有效性, 在此基础上对与此相对应算例条件下的定向凝固问题进行了初步的研究。 数值模拟与实验进行了对比, 结果比较令人满意。
1 物理模型
1.1 基本假设
如图1所示, 功率为定值P 的高强度激光光束以固定的速度U 沿x 轴方向扫过金属表面, 产生的高温使材料熔化, 在金属表面形成熔池。 在本文计算模型中我们作如下假设。
图1 物理模型示意图
Fig.1 Sketch map of physical model
1) 熔池内为牛顿流体;
2) 为不可压缩流体;
3) 材料的物性参数值不随温度变化且固液态条件下取值相同, 只在考虑浮升力的作用时认为密度随温度变化, 亦即采用Boussinesq假设;
4) 熔池内流动为层流;
5) 熔池顶部为平面;
6) 流动的驱动力只考虑浮升力和熔池自由表面的表面张力。
1.2 控制方程和边界条件
稳态时, 描述工件上的流场和温度场的控制方程可以写成如下形式
[7 ]
:
连续方程
?·V =0 (1)
动量方程
ρ (V ·?) V =-?p -[?·τ ]-ρβg (T -T 0 ) (2)
能量方程
V ·? (ρH -p ) =?· (k ?T ) (3)
式中 V 是速度, ρ 为密度, p 为压力, τ 为剪切力, β 为材料的体积热膨胀系数, g 为重力加速度, T 为绝对温度, T 0 为参考温度, H 为单位体积的热焓, k 为导热系数。
边界条件如下:
1
)
?
件
上
表
面
(
z
=
0
)
w
=
0
μ
?
u
?
z
=
-
?
Τ
?
x
?
γ
?
Τ
μ
?
v
?
z
=
-
?
Τ
?
y
?
γ
?
Τ
在光斑区
-
κ
?
Τ
?
z
=
ε
q
˙
?
?
(
r
≤
a
)
在其他区域
-
κ
?
Τ
?
z
=
α
(
Τ
-
Τ
a
)
+
σ
e
(
Τ
4
-
Τ
a
4
)
?
(
r
>
a
)
2
)
?
件
底
面
(
z
=
δ
)
u
=
U
v
=
w
=
0
-
κ
?
Τ
?
z
=
α
(
Τ
-
Τ
a
)
+
σ
e
(
Τ
4
-
Τ
a
4
)
3
)
扫
描
轨
道
中
心
面
(
y
=
0
)
v
=
0
?
u
?
y
=
?
w
?
y
=
?
Τ
?
y
=
0
4
)
?
件
侧
面
(
y
=
±
l
)
u
=
U
v
=
w
=
0
-
κ
?
Τ
?
y
=
α
(
Τ
-
Τ
a
)
5) 光斑的远前方 (x →-∞)
T =T 0
u =U
v =w =0
6
)
光
斑
的
远
后
方
(
x
→
∞
)
?
Τ
?
x
=
0
u
=
U
v
=
w
=
0
上述式中 u , v , w 分别是速度V 在x , y , z 方向上的投影; μ 是液态金属的粘性系数; γ 是表面张力;
q
˙
是激光功率; ε为材料表面对激光的吸收系数; a是光斑直径; α为表面换热系数; σ为斯蒂芬-波尔滋曼常数; e为材料表面总辐射系数; Ta 为环境温度; l为工件宽度之半, 在本计算中假设工件y尺度很大, 故可认为关于A—A面成对称, 且在y=±l因距扫描中心线较远, 温度较低, 其辐射可以忽略; U为激光束扫描速度。
2 求解方法
本文所求解的控制方程具有两大特点: 其一, 方程组是按稳态条件写出来的, 激光光斑的移动或工件的反向移动通过边界条件反映出来, 即其移动速度隐含在V中; 其二, 固态和液态的方程形式是统一的。 因此求解的过程不需要追踪熔池界面的位置。 另一方面, 固液两相存在着两大区别, 一是液相存在流动, 而固相没有; 二是液固两相的焓不同。 因此, 还需要引入焓—温度, 粘性—温度两个关系。
2.1 焓—温度关系
焓与温度的关系主要考虑液相有凝固潜焓, 而固相没有。
对于式 (3) 中的焓H , 有
[8 ]
H =h +ΔH
式中 h =c p T 为显焓, ΔH =f (T ) 为潜焓。 当考虑糊状区时, 潜焓与糊状区的液相分数有关,
f
(
Τ
)
=
{
L
(
Τ
≥
Τ
l
)
L
f
l
(
Τ
l
>
Τ
≥
Τ
s
)
0
(
Τ
≤
Τ
s
)
式中 T l 是液相温度, T s 是固相温度, f l 为液相分数, L 为相变焓。 当不考虑糊状区, 即采用等温线凝固时 (本文方法) ,
f
(
Τ
)
=
{
L
(
Τ
≥
Τ
m
)
0
(
Τ
<
Τ
m
)
式中 T m 为等温熔点。
这样处理之后, 以焓的形式表达的能量方程式 (3) , 在液、 固两相都能使用。
2.2 粘性—温度关系
方程 (2) 是由N-S方程推导出来的, 在液相中可以直接应用, 而在固相中实际上没有流动存在, 故不能直接应用。 我们通过引入有效粘性的概念, 来统一液固两相的运动方程。
有效粘性是由粘性—温度关系得到的, 这一关系类似于焓—温度关系:
μ
e
f
f
=
{
μ
?
?
?
?
(
Τ
≥
Τ
l
)
C
(
1
-
f
l
)
+
μ
(
Τ
s
<
Τ
<
Τ
l
)
C
?
?
?
?
(
Τ
≤
Τ
s
)
式中 C 为一个很大的数, 其取值视具体计算情形而定, 既不能太大不易收敛, 也不能太小产生物理上的不真实解
[9 ]
, 本文中C =108 。 当采用等温线凝固时,
μ
e
f
f
=
{
μ
(
Τ
≥
Τ
m
)
C
(
Τ
<
Τ
m
)
可以看出, 在方程 (2) 中采用有效粘性后, 液相中方程不变 (μ eff =μ ) , 固相中由于μ eff =C 为一个很大的值, 就使得流动速度 (V -U ) 衰减为非常小 (近似为零) , 这样方程 (2) 既可以用来描述液相也可以描述固相, 达到了两者的统一。
2.3 数值求解
基于本文采用的控制方程的特点, 以及上述焓—温度、 粘性—温度关系的引入, 使得计算过程不需要人为追踪固-液相界面的位置, 从而使求解更加容易。 本文采用SIMPLE算法
[10 ]
求解速度场。 为了适应流场及温度梯度的变化, 采用变间距网格, 靠近激光光斑的位置用细网格, 远离光斑的位置用粗网格。 这样, 既可以提高计算精度也可以加快收敛速度。
3 结论和分析
为了验证方法的可靠性, 计算了两种条件下激光熔凝的熔池深度, 并通过实验对其进行了测量。 实验一: 基体材料为Al-1.2%Mn, 激光功率为1.5 kW, 光斑直径为0.33 mm; 实验二: 基体材料为Al-3.2%Mn, 激光功率为1.5 kW, 光斑直径为0.3 mm。
图2和图3所示是两种不同材料熔池深度与扫描速度之间的关系。 可以看出计算与实验的吻合程度还是令人满意的。 由此可以验证本文计算的可靠性。
图2 熔池深度与扫描速度的关系
Fig.2 Relationship between depth of melting pool and scanning velocity (Al-1.2%Mn)
图3 熔池深度与扫描速度的关系
Fig.3 Relationship between depth of melting pool and scanning velocity (Al-3.2%Mn)
图4是计算得到的Al-1.2%Mn合金在激光输出功率为3.25 kW, 扫描速度为24 mm/s时熔池纵截面 (A —A 面) 固液界面上的温度梯度分布图, 其x 轴是沿着激光扫描方向离开光斑中心的距离, z 轴是工件内某一点离开激光辐照表面的距离 (见图1坐标) 。 通过该图可以预见在熔池底部胞晶组织将垂直向上生长, 随着离开熔池底部距离的增加, 其生长方向将逐渐发生偏转, 到熔池顶部则变成与激光扫描方向基本一致。
在一定条件下熔池顶部的温度梯度方向变得与激光扫描速度方向一致, 这一事实使得有可能发展出一种材料的激光定向凝固技术, 其特点是具有常规定向凝固技术不具备的超高温度梯度和很快的凝固速度, 这对发展新材料及凝固理论研究都具有重要意义。
图4 熔池固液界面上温度梯度的方向
Fig.4 Distribution of temperature gradient on solid/liquid interface of melting pool (Al-1.2%Mn, P =3.25 kW, U =24 mm/s)
图5所示为在同样工艺条件下得到的Al-1.2%Mn合金的熔池纵截面金相组织照片。 可见其生长方向与计算分析结果一致。
图5 熔池纵截面金相组织照片
Fig.5 Microstructure photograph of longitudinal section of melting pool (Al-1.2%Mn, P =3.25 kW, U =24 mm/s)
为了进一步研究激光熔凝各工艺参数与定向凝固组织的关系, 计算了Cu-31.4%Mn合金在相同的激光功率不同的扫描速度下和相同的扫描速度不同的激光功率下, 纵截面上熔池固液界面上温度梯度的分布。 图6所示为计算得到的Cu-31.4%Mn合金在激光输出功率为1.5 kW, 扫描速度分别为5, 19, 30 mm/s时, 熔池中固液界面处温度梯度分布。 图7所示是计算得到的Cu-31.4%Mn合金在扫描速度为24 mm/s, 激光输出功率为1.5, 2.0, 2.75 kW时, 熔池中固液界面处温度梯度分布。 由图可见, 功率一定时, 在扫描速度较低时熔池表面的温度梯度方向与扫描速度方向一致, 随着扫描速度的增加, 熔池顶部温度梯度的方向明显偏离扫描速度方向, 因而当扫描速度增大后, 在熔池顶部将不再能得到与扫描速度方向一致的定向凝固组织。 要想改变温度梯度的方向, 需要适当增大功率。 同样, 在一定的扫描速度下, 功率太小时, 熔池顶部的温度梯度方向不与扫描速度的方向平行, 因此也不能得到与扫描速度方向一致的定向凝固组织。 通过图6和图7我们可以看出, 要想在熔池顶部获得与扫描速度方向一致的温度梯度, 必须要求激光功率与
图6 相同的激光功率不同的扫描速度下纵截面上熔池界面温度梯度分布
Fig.6 Distribution of temperature gradient on solid/liquid interface at various scanning velocity (P =1 500 W, Cu-31.4%Mn)
图7 相同的扫描速度不同的激光功率下纵截面上熔池界面温度梯度分布
Fig.7 Distribution of temperature gradient on solid/liuqid interface at different laser power (U =24 mm/s, Cu-31.4%Mn)
扫描速度相互匹配。
据此对激光功率与扫描速度进行调整匹配, 并在基体Cu-31.4%Mn进行了实验。 图8是基体为Cu-31.4%Mn, 激光功率P =1.0 kW, 扫描速度U =5.1 mm/s时熔池纵截面组织。 由图可看出, 重熔区组织生长非常规则, 从熔池底部到顶部胞晶组织生长方向受熔池中温度梯度场的影响逐渐偏转, 最后在顶部变得完全与激光束的扫描方向一致。 可见, 通过激光功率与扫描速度的匹配, 在熔池顶部可以得到与扫描方向完全一致的定向组织。
图8 熔池纵截面组织
Fig.8 SEM photos of longitudinal section (Cu-31.4%Mn, P =1.0kW, U =5.1mm/s) (a) —Bottom; (b) —Center; (c) —Top
4 结论
1) 建立了一种激光表面熔凝的数学模型及计算方法, 与实验测得的熔池深度的比较证明了数值模拟的有效性。
2) 初步分析了激光定向凝固与工艺参数的关系, 通过激光功率与扫描速度的匹配, 在熔池顶部中心处得到了与扫描速度方向一致的凝固组织。
参考文献
[1] KurzW , FisherDJ.FundamentalsofSolidification[M ].AedermannsdorfTech, 1989.5-8.
[2] ChongE , WilcoxWR .ControlofinterfaceshapeintheverticalBridgamanStockbargertechnique[J].JCrystalGrowth, 1974, 21:135.
[3] FlemingsFC .SolidificationProcessing[M ].NewYork:McgrawHill, 1974.
[4] ChanCL , MazumderJ, ChenMM .Three dimensionalaxisymmetricmodelforconvectioninlaser meltedpools[J].MaterialsScienceandTechnology, 1987, 4:306-311.
[5] HoadleyAFA , RappazM , ZimmermannM .Heat flowsimulationoflaserremeltingwithexperimentalvalidation[J].MetallurgicalTransactionsB , 1991, 22B :101-109.
[6] PaulA , DebroyT .Ferrsurfaceflowandheattransferinconductionmodelaserwelding[J].MetallurgicalTrans actionsB , 1988, 19B :851-858.
[7] KouS , WangYH .Three dimensionalconvectioninlasermeltedpools[J].MetallurgicalTransactionsA , 1986, 17A :2265-2270.
[8] VollerVR , Parkash.Afixedgridnumericalmodelingmethodologyforconvection diffusionmushyregionphase changeproblems[J].IntJHeatMassTransfer, 1987, 30 (8) :1709-1718.
[9] 张寅平, 胡汉平, 孔祥冬, 等.相变储能[M].合肥:中国科技大学出版社, 1996. ZHANGYin ping, HUHan ping, KONGXiang dong, etal.EnergyStoragebyPhaseChange[M ].Hefei:U niversityofScience&TechnologyofChinaPress, 1996.
[10] 陶文铨.数值传热学[M].西安:西安交通大学出版社, 1988. TAOWen quan.NumericalHeatTransfer[M ].Xi′an:Xi′anJiaotongUniversityPress, 1988.