激光焊接传热过程的数值计算讲解.docx
- 文档编号:14414148
- 上传时间:2023-06-23
- 格式:DOCX
- 页数:9
- 大小:284.52KB
激光焊接传热过程的数值计算讲解.docx
《激光焊接传热过程的数值计算讲解.docx》由会员分享,可在线阅读,更多相关《激光焊接传热过程的数值计算讲解.docx(9页珍藏版)》请在冰点文库上搜索。
激光焊接传热过程的数值计算讲解
第29卷增刊
1999年11月
东 南 大 学 学 报JOURNALOFSOUTHEASTUNIVERSITY Vol129Sup.Nov.1999
激光焊接传热过程的数值计算3
徐九华 罗玉梅(,摘 要 ,建立了运动,提出采用位置预置、液相交界面位置进行准确捕捉,并对这一小孔焊接,揭示了材料热物理性能、小孔直径、焊接速度等因素对焊接热过程的影响.
关键词 传热;温度场;数值模拟;激光焊接
分类号 O551.3;O241.8
Ξ国家自然科学基金资助项目(59575069.
收稿日期:
1999-06-14. 第一作者:
男,1964年生,博士,副教授.
当功率密度超过109W/m2的激光束流照射到金属材料表面时,材料将瞬时汽化并在束流压力和蒸汽压力的共同作用下形成一个细长的柱形小孔,小孔中的汽化金属被电离并将射入的能量完全吸收,然后将热量传递给周围材料使其熔化,在小孔周围形成熔池.这种现象常出现在激光束焊过程中.这类焊接能获得极小的熔化区和热影响区并能净化焊缝,从而减少焊缝中存在的诸如内应力、裂纹、气孔等缺陷.因此,激光焊接在高精度、高质量焊接领域中引起了人们的广泛关注.
大量的实验研究表明,在激光焊接中,所形成的小孔深度大于其直径一个数量级,并且熔化金属液主要沿小孔周向流动.因此,在早期的理论研究中,这类问题被模化为具有移动热源的二维导热问题,仅考虑了热扩散作用[1,2];随后相变和流动效应被引入到一系列研究之中[3~5].本文是在已有的研究基础上,借助于曲线坐标系,运用有限控制容积积分法对激光小孔焊接热过程进行数值计算,并对这一热过程进行全面的参数化分析,深入剖析了材料热物理性能、焊接主要过程参数对焊接热过程的综合影响,揭示了以往研究中未曾涉及到的激光小孔焊接模式的一些新特征.
1 数学模型的建立
图1给出了激光小孔焊接过程的物理模型,采用随热源移动的运动坐标系.考虑到对称性,仅取一半区域进行计算.为使问题简化,特作如下假设:
1高能束流以恒速(uw移动,在运动坐标系下,传热过程视为稳态;
2小孔的形状可视为圆形(半径为ra,小孔表面的温度为材料沸点温度;
3熔化金属为牛顿型流体,不可压缩;固相区域和液相区域物性均匀,各向同性,且热物
1性为温度的函数;
4在计及焊件与周围环境换热损失中,忽略辐射换
热损失,仅考虑自然对流情况,且假定对流换热系数为常
量.
采用随热源移动的运动坐标系,可将传热过程视为
稳态.图中ra为小孔半径,rm为固液交界面半径,δLS
分别表示固相层和液相层厚度.可模化为运动坐标系下的稳态、,用控制方程如下:
(v<Γ+yΓ+S(x,y(1式中,uy;和S(x,y为变量<所对应的扩散系数和源项.对于固相区,TS,u=uw,v=0,式(1转化为扩散型导热方程.
相应的边界条件为
y=0:
y=0,v=0;y=0,y=0(2 r=ra:
ur=0,rur
=0,T=Tb(Tb为材料沸点温度(3 r→r∞:
T∞=Ta(Ta为环境温度
(4 r=rm:
u=uw;v=0;TL=TS=Tm(Tm为材料熔点温度
(5在固、液交界面上,还存在一反映局部能量平衡的关系式:
-λsn=-λLn
+ρuwhslcosθn(6式中,hsl为材料熔化潜热;θn为x轴正向和固液交界面法线之间的夹角.
2 计算过程
由于焊接熔池中固、液相交界面位置的不确定性,计算过程分为两部分:
首先是在假定的相界面位置上,对固、液两相区域内的变量各自进行独立求解;其次是确定满足包括界面能量平衡式在内的所有控制方程和边界条件的固、液相界面位置.
2.1 坐标变换
考虑到计算域的不规则性,采用坐标变换法将计算过程转化到规则的计算平面上进行.引
入曲线坐标系x=x(ζ,η
y=y(ζ,η,通用控制方程式(1可以变换为 J
(ρU<+J(ρV<=J+JJ-JJ+J+S(ζ,η(7式中,U=uyη-vxη;V=vxζ-uyζ;α=x2η+y2η;β=xζxη+yζyη;γ=x2ζ+y2ζ;
J=xζyη-yζxη.
网格生成采用代数法.取r∞=Crm,C为定值,其取值范围为200~500.坐标变换采用如
36增刊徐九华等:
激光焊接传热过程的数值计算
下的变换关系式:
r=ra+δLζ
0≤ζ≤1ra+δL+δS(ζ-1 1≤ζ≤2(8
θ= 0≤η≤π
式中,δL和δS分别表示固相层和液相层厚度.
2.2 固、液交界面位置捕捉
固、液界面位置捕捉采用位置预置在曲线坐标系下, =-(9、能量守恒条件,有
^=^FS,^EL=^ES(10其中,上标“^”表示交界面的参数,h为热焓.
对于任一假定的固、液交界面位置,式(10一般难以满足,故需对固、
液交界面位置进行修正.记固、液交界面位置修正系数为ωm,则有rNm=rm(1+ωm,rN∞=r∞(1+ω∞,ω∞=ωm.
上标“N”表示对应于修正后的变量.
将上述表达式代入式(10并对固、
液交界面的有限控制体进行积分,经推导可得到位置修正式如下:
ωm=^E-^EPm,Pm=rδL(^E-^F^hL[6](11
在计算中,为保证ωm的绝对值小于1,提高计算稳定性,对式(11加以改造,令
Pm=|^EL|+|^ES|
固、液交界面位置修正,应能满足:
|ωm|≤0.001.
3 计算结果与分析
利用本文所建立的模型分别对不锈钢(AISI304(算例1和铝基复合材料(20%SiCp+Al6061(算例2的激光小孔焊接热过程进行数值计算.对于不锈钢,密度ρ=7200kg/m3,熔点Tm=1427℃,沸点Tb=2807℃,粘度μ=6142×10-3kg/(m・s,潜热hsl=26512kJ/kg.对于铝基复合材料,密度ρ=2858kg/m3,熔点Tm=627℃,沸点Tb
=2457℃,熔化潜热hsl=395kJ/kg,分子粘度μ=110×10-3kg/(m・s.2个算例的导热系数和比热都是温度的函数,在每一个计算步骤上都要根据当时的温度分布重新确定其取值.这些参数随温度变化的函数关系及其它计算参数取值见文献[6].
图2和图3为在不同焊接速度下计算出的不锈钢板小孔焊接的温度分布.显然,对于本算例,温度分布的等值线具有椭圆形特征,并且在热源前端的等温线要比后端的密集,这与具有移动热源导热的理论分析是一致的.图示结果还表明:
随着焊接速度的提高,熔池形状和焊接温度场的等温线的椭圆形特征越显著.
图4和图5为本文模型对算例2的计算结果(铝基复合材料,ra=0105mm,uw=210
46东南大学学报第29卷
21(uw=0.5mm/s
图3 算例1的焊接温度场(uw=2.5mm/s
mm/s.图4为理论计算所得到的焊接熔池表面形状.图5为理论计算所得到的焊接温度场.图4中熔池宽度的理论计算值约为0155mm,与实验值0153mm十分接近,如图6所示.
图4和图5表明:
在本算例所给定的材料物性和焊接工艺参数下,所获得的熔池和等温线形状不同于算例1,它们的椭圆形特征并不明显,而更接近于圆形.
图4 本文模型对算例2计算出的熔池边界
由上述算例可以得出:
焊接熔池的几何形状
和接温度场分布特征是受材料热物理性能、主要
焊接过程参数(在本文中主要分析了焊接速度和小孔直径两个参数综合作用的.对于同样的焊
接材料,焊接熔池和热影响区的大小将随着小孔
直径的减小而减小,随焊接速度的增大而减小.对
不同的焊接材料,它们的热扩散率不同,扩散传热
的能力也不同,所获得的焊接熔池的形状也不同.因此,可以考虑引入一个能全面反映焊接速度、小孔直径、材料热扩散率三者对焊接热过程的综合作用的无量纲化参数Pe数作为一个基本变量,并定义Pe=uraL
aL为材料热扩散率(m2/s.对于算例1,Pe=
015237(当uw=115mm/s,ra=2mm时,而对于算例2,Pe=010025(当uw=210mm/s,ra=0105mm.对照2个算例所获得的熔池和等温线形状,可以说明当Pe数越大,熔池和等温线越接近于椭圆形并更为扁平;Pe数越小,熔池和等温线的椭圆形特征越不明显,更接近于圆形.关于这一现象,可以由Pe数的定义来解释:
其分母是材料的热扩
5
6增刊徐九华等:
激光焊接传热过程的数值计算
2计算出的焊接温度场
图6 焊接速度对熔池宽度的影响散率a;其分子是焊接速度和
.分子越大,
这种流动的能力越强,对流传热所占的比例越大,熔池和等温
线的流线形特征越明显,此为大Pe数的情形.相反地,分母越
大,焊接熔池中液态金属的扩散传热的能量就越强,导热所占
的比例就大,对流则不太明显,故熔池和等温线接近于圆形,此为小Pe数的情况.当Pe数降至0,则对应于固定热源作用下的焊接热过程,此时,熔池应呈标准的圆形.
4 结 论
1建立了运动热源作用下二维激光小孔焊接熔池中流体流动及传热过程的数值分析模型,并在曲线坐标系下,采用位置预置—修正的方法成功地实现了焊接熔池界面位置的精确捕捉.
2利用本文所建立的数值分析模型对两种不同材料、不同焊接过程参数的焊接热过程进行了数值计算.计算结果表明,本文的计算方法是合理可行的.
3对激光小孔焊接热过程进行了较为全面的参数化分析,探讨了焊接材料热物理性能、主要焊接过程参数等对焊接熔池形状、焊接温度分布趋势的综合影响.并提出了一个可以综合反映各参数之间关系的无量纲参数Pe数.
参考文献
1 Swift2HookDT.Penetrationweldingwithlasers.WeldingJournal,1973
52:
492~499
2 MazumderJ,SteenWM.HeattransfermodelforCWlasermaterialprocessing.JournalofAppliedPhysics,1980,51:
941~947
3 DavisM,KapadiaP,Dowden.Modellingthefluidflowinlaserbeamwelding.WeldingJournal(ResearchSupple2ment,1986,167~174
4 HsuYF,RubinskyB.Two2dimensionalheattransferstudyonthekeyholeplasmaarcweldingprocess.IntJournalofHeatandMassTransfer,1988,31:
1409~142166东南大学学报第29卷
增刊徐九华等:
激光焊接传热过程的数值计算67tionisemployedtoaccuratelyidentifythesolid2liquidinterfacelocationduringlaserwelding.Aparamet2ricanalysisoftheexistingkeyholemodelisperformed,andtheinfluenceofparameterssuchasmaterialcess.Keywords:
heattransfer;temperaturedistribution;numericalsimulation;taserweldingproperties,weldingspeed,keyholeradiusontheheattransferisanalyzedduringkeyholeweldingpro25 Charn2JungKim.Parametricstudyoftwo2dimensionalkeyholemodelforhighpowerdensityweldingprocesses.JournalofHeatTransfer,1994,116:
209~2146 罗玉梅.激光焊接传热过程数值研究:
[学位论文].南京:
南京航空航天大学图书馆,1999NumericalSimulationandParametricStudyfortheHeatTransferinKeyholeHighPowerDensityWeldingProcessXuJiuhua LuoYumei ZhangJingzhou(NanjingUniversityofAeronauticsandAstronautics,Nanjing210016theheattransferinkeyholelaserweldingprocesses.Arecentlydevelopedmethodoftheposition2correc2Abstract:
two2dimensional,quasi2stationarynumericalmodelispresentedtostudythefluidflowandA
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 激光 焊接 传热 过程 数值 计算 讲解
![提示](https://static.bingdoc.com/images/bang_tan.gif)