浙江国检检测

首页 检测百科

分享:铝-镁异种合金搅拌摩擦搭接焊残余应力数值模拟

2022-09-26 14:05:25 

吴晓虎1,曹丽杰1,苗臣怀1,王 仪2

(1.上海工程技术大学机械与汽车工程学院,上海 201620;2.上海飞机制造有限公司,上海201324)

摘 要:采用顺序热力耦合的方法建模,通过正交试验,利用 ABAQUS-合金搅拌摩擦搭接焊中焊接速度搭接量搅拌头转速对残余应力的影,的工,过试验验证模拟结果的准确性,并研究焊接速度搭接量对残余应力的影响规律结果表明:搅拌 摩擦搭接焊的最优参数为焊接速度60mm·min-1搭接量60mm搅拌1400r·min-1;不同焊接工艺下,热循环曲线和残余应力的模拟结果与试验结果相吻合,7.5% 8.4%,验证了数值模拟结果的准;,最优焊接艺下的最大残余应力为137.7MPa;,的影响较 ,随着焊接速度的增大,纵向压应力峰值升高,压应力作用范围变窄,而随着搭接量的增加, 纵向残余压应力峰值降低,压应力作用范围变宽

关键词:异种合金搅拌摩擦搭接焊;残余应力;数值模拟

中图分类号:TG457 文献标志码:A 文章编号:1000-3738(2022)07-0095-08

0 引 言

搅拌摩擦焊是一种新型而有效的先进固态连接 工艺,在焊接过程中板材无需熔化即可连接传统熔焊技术难以焊接的高强度铝合金以及异种材料,与传统的熔焊相同,搅拌摩擦焊接头中也存在热 影响区和焊后残余应力[1]在搅拌摩擦焊过程中, 焊接接头内部产生的残余应力通过改变焊接工件的 内应力和力学性能来破坏工件的整体性能,导致工 件的 结 构 和 性 能 严 重 受 损陈 彦 君 等[2] 基 于 ABAQUS,模型,6, 且试验结果与仿真结果吻合较好[3]拌摩擦焊产热模型,研究了焊接参数和搅拌头偏置 对铝-钢异种金属焊接过程的影响,发现是残余应力的主要影响因素周文静[4]合金板搅拌摩擦焊的稳态与瞬态过程,发现瞬态与 稳态计算结果基本一致,随着主轴转速与轴肩下压 力的增大,应力与温度的非对称性逐渐增强史学 海等[5]采用热力耦合方法建立了冷热源辅助的搅拌 摩擦焊模型,发现与常规工艺相比,冷热源辅助可以 有效降低接头焊后残余应力峰值目前主要针对搅 拌摩擦对接焊过程进行数值模拟,而搭接也是搅拌 摩擦焊接结构中一种常见接头形式张昭等[6]建立 铝合金搅拌摩擦搭接焊温度-再结晶-析出相-力学性 能的一体化计算模型,实现了基于晶体塑性理论的 力学性能研究,通过试验和数值模拟的对比验证了 模型的有效性胡云瑞等[7]研究发现,2024-T4 铝合金搅拌摩擦搭接焊过程中,搅拌头轴肩边缘处 的材料流动速率较大,并且流速随着距搅拌头距离 的增加而减小但是,关于搅拌摩擦搭接焊的研究 主要集中在同种金属方面,有关异种金属搅拌摩擦 搭接焊的数值模拟研究较少因此,作者以6082合金与 AZ31B镁合金为研究对象,通过正交试验, 采用 ABAQUS行数 值模拟,后残余应力场的影响,获得最佳的工艺参数,并进行 了试验验证;研究了焊接速度搭接量对残余应力的 影响规律,以期为优化焊接参数提供理论指导

1 数值模型的建立

1.1 热源模型

在搅拌摩擦焊接过程中,热量主要来源于搅拌 头与工件接触面上的摩擦产热以及工件剪切层内的 少量塑性变形功[8]产热总功率Qt 的计算公式为 Qt =ηMzω=n6η0Mz (1) 式中:η为热传输效率,0.95[9];Mz 为搅拌头的转;ω 为转动角速度;n 为转速

搅拌头采用平轴肩圆柱状的热源模型,轴肩半 径为5 mm,搅拌针半径为 1.5 mm,轴肩部分的摩擦产热约占75%,变形 ,产热总 功 率 约 等 于 轴 肩 与 搅 拌 针 产 热 功 率 之 [9]轴肩的产热功率Qs 可表示为 Qs =0.75Qt (2) qs(r)= 3Qsr (R31 -R30) (3) 式中:qs 为轴肩热流密度,由于轴肩与焊件之间为 面接,,qs ;r 中心;R1 ;R0 为搅拌针半径搅拌针的产热功率Qp 可表示为 Qp =0.25Qt (4) qp(r)= Qp πR20h (5)式中:h 为搅拌针长度;qp 为搅拌针热流密度,搅拌 针插入焊件内部高速转动并向前移动,属于体热源, 因此qp 流密度

1.2

1.2.1 材料参数

6082[10]AZ31B [11]在不同温 度下12



1.2.2 网格划分

采用有限元软件 ABAQUS 热力学模型模型采用尺200mm×100mm× 3mm 2块异种材料板搭接的形式,6082合金板在下侧,AZ31B镁合金板在上侧为了提高 计算效率并兼顾模拟精度,对焊缝中心处及附近的 网格进行局部加密处理,模型网格划分如图1所示, 网格划分为73200个单元,91254个节点,单元类 型选择为热-力耦合单元 C3D8R

1.2.3 边界条件

在模型中,温度场的模拟计算采用等效热力学 边界条件代替夹具与垫板;铝合金板与下表面的垫 板直接接触,散热较快,因此铝合金板与垫板之间采 用较高的 对 流 换 热 系 数[12],设 为 1000 W·m-2· K-1;铝合金板侧面和镁合金板完全暴露在空气中, 采用较低的对流换热系数,设为30W·m-2·K-1;件初始温度设为室温(20 )。与对接方式相比,拌摩擦搭接焊在接头处具有更复杂的边界条件搅拌头下方网格固定约束,其余搭接面的接触属性 定义为2 个 光 滑 板 材 之 间 的 接 触,摩 擦 因 0.3焊接结束后,当温度趋于稳定时,除垫 板和夹具的相关约束在焊接开始时选择少量节点 约束各个方向的自由度,在焊缝末端选取少量节点 只约束其垂直方向和水平方向的自由度,使工件能 够沿焊缝方向变形以避免刚体位移[13]

2 参数优化与试验验证

2.1 参数优化

焊接速度搅拌头转速和搭接量等焊接参数对 搅拌摩擦搭接焊接头残余应力的大小与分布有重要 影响正交试验设计在多因素多水平的试验中具有 独特的优点,能够在不影响试验效果的同时,尽可能 减少试验次数[14]采用3因素4水平[L16(43)]形式,设计搅拌摩擦搭接焊参数的因素水平表,具体 如表3


残余压应力对搅拌摩擦焊接接头的疲劳寿命有 ,而残余拉应力则可以增加裂纹尖端的应力分布, 严重影响搅拌摩擦焊接头的疲劳寿命[15]搅拌摩 擦焊接头的残余拉应力主要以纵向残余应力为主, 因此模拟过程中主要研究纵向残余拉应力在正交 试验中,信 噪 比 作 为 特 征 数 来 衡 量 焊 接 接 头 的 质 [14]由信噪比的望小特性可知,纵向残余拉应力 小越,即信噪比越小越好信噪比RSN 的计算 RSN =-10lg 1aai=1σ2yi (6) :a ;σyi i 次数值模拟得 到的

由表4,按从大到小,最优的设计方案为 C1A4B3,即焊接速度60mm· min-1,搭接量60mm,搅拌头转速1400r·min-1优方模拟,得到的大值137.7MPa

2.2 试验验证

为了验证数值模拟的准确性,在专用搅拌摩擦 焊接机上进行搅拌摩擦搭接焊试验,焊件材料和尺 寸与模型保持一致,搅拌头材料为 H13,拌头参数选择与热源模型参数相同结果,搅拌头转速对纵向残余拉应力的影响较小,速不,, 接试如表5使K 型热电偶记录焊接中间阶段垂直于焊接方向 距焊缝中心线约10mm 处的瞬态,直径分别为2mm 1 mm,线 10mm



由表4,按从大到小,最优的设计方案为 C1A4B3,即焊接速度60mm· min-1,搭接量60mm,搅拌头转速1400r·min-1优方模拟,得到的大值137.7MPa2.2 试验验证 为了验证数值模拟的准确性,在专用搅拌摩擦 焊接机上进行搅拌摩擦搭接焊试验,焊件材料和尺 寸与模型保持一致,搅拌头材料为 H13,拌头参数选择与热源模型参数相同结果,搅拌头转速对纵向残余拉应力的影响较小,速不,, 接试如表5使K 型热电偶记录焊接中间阶段垂直于焊接方向 距焊缝中心线约10mm 处的瞬态,直径分别为2mm 1 mm,线 10mm


将模拟得到线结果进行对比2,,得到测量点的热循环曲线与试验得到的热循环曲线 吻合较好,相对误差小于7.5%,表明该热源模型可 用来预测搅拌摩擦搭接焊的温度分布规律

采用盲孔法[16]测量焊缝中心处的表面纵向残 余拉应力,并与模拟结果进行对比,结果如表6由表6可知,纵向结果相差较小,大相8.4%,差为3.5%,验证了数值模型的准确性,型及模拟结果对实际焊接有一定的指导意义


3 模拟结果与讨论

3.1 残余应力的分布规律

由图3可知,焊缝中心附近的残余应力集中且 明显高于其他区域,纵向残余应力呈典型的不对称 M 型分布,最大残余应力为102.5MPa,位于前进侧 (AS),退(RS),小最终变为压应力横向残余应力整体上余应力分布趋势相同,但焊缝中心处的横向力主要为压应力,最大值为18.1 MPa,进侧的轴肩边缘附近在焊接过程中,向上的温度梯度变化比其他方向上的变化快[17],此纵向残余应力明显高于横向残余应力


由图4可知,在轴肩范围以内,接头上表面与搭 接面处 的 温 度 分 布 曲 线 几 乎 相 同,最 高 温 度 为 443.3 ,明显高于下表面的最高温度,而在轴肩围以外,上表面搭接面和下表面的温度分布曲线几 乎相同在进行异种金属搅拌摩擦搭接焊时,上下 板之间的搭接面存在接触热阻,热量在搭接面处传 递会发生突变,此时温度梯度大于均匀传递时的温 度梯度,从而导致下板温度过低,容易造成接头的弱 连接[18]在数值模拟中发现最大残余应力基本出现 在焊接末端的搭接面处,主要原因是搭接面的温度梯 度较大随着焊接的进行,热量不断积累,焊接末端 温度快速升高,温度梯度进一步扩大,因此最大残余 应力主要出现在焊接末端的搭接面处

3.2 焊接速度对残余应力的

在正交试验结果分析中,现焊接速度和搭接,响较,表面 残余应力以纵向残余应力为主,因此只分析焊接速 度和搭接量的变化对纵向残余应力的影响保持搭 接量60mm搅拌头转速1400r·min-1,变焊接速度为60,80,100,120 mm·min-1,接速度对温度分布曲线和纵向残余应力的影响


由图5,,度降低当焊60mm·min-1,443.5 ;度从 80 mm·min-1 100mm·min-1 ,温 度 424.6 ℃ 降 低 至 397.5 ;速度继续增大至 120 mm·min-1 ,321.2 度焊缝的热输入,随着,缝相同位置的作用时间减少,单位长度焊缝的热输 变小,,接过程中,匀程,收缩不同 程度,[19]


由图6可知,焊接速度的变化对纵向残余应力 有明显影响随着焊接速度的增大,纵向残余应力 峰值增大,当焊接速度较低时,残余应力的 M 布特征较明显当焊接速度从60 mm·min-1 80 mm·min-1 ,,,;着焊 接速度的增大,残余压应力的变化较小残余拉应 力的作用区域对焊接速度较敏感,主要原因是在较 焊接焊缝中心的材,,轴肩用范材料塑性流动更充分并且熔合更均匀,残余拉应力 的作用区域相对于较高焊接速度分布更广泛[20]


3.3 搭接量对残余应力的影响

保持 焊 接 速 度 60 mm·min-11400r·min-1 不变,仅改变搭接量,量对 接头纵向残余应力的影响由图7可以看出,的即时峰值温度随着搭接量的增加而降低不同搭 接量接头的即时峰值温度随时间的变化趋势相同, 即随着焊接时间的延长,即时峰值温度在初期迅速 上升,25s左右到达焊接稳定,变缓,由于在整个焊接过程中,热略,此即时峰值温度在焊接稳定期一直平稳且缓慢上 ;当热源移动至焊接末端时,由于热量的累积,接末期的 即 时 峰 值 温 度 迅 速 升 高,最 高 温 度 达 到 600 以上,高于焊件的熔点过高的温度易产生 焊接缺陷,在实际的焊接过程中应避免出现末端的 高温现象[13]焊 接 结 束 后,因 没 有 持 续 的 热 量 输 ,所以焊件的温度迅速降低搭接量的变化主要 影响焊件之间的相互作用与传热效率随着搭接量 的增加,两板的接触面积增大,热量传递更加均匀, 焊接过程中的温度梯度减小,因此焊缝表面的残余 拉应力降低此外,垫板对焊接过程中的温度传递 和散热有重要影响,随着搭接量的增加,通过垫板散 失的热 量 增 加,因 此 焊 件 即 时 峰 值 温 度 略 微 降 [18]


由图8可知,随着搭接量的增大,接头的最大残 余拉应力减小,残余拉应力作用的范围扩大当搭 接量从30mm 增大至40mm ,残余拉应力峰值 133.8MPa降至118.0MPa,下降趋势较明显,余压应力变化不明显;当搭接量从 40 mm 增大50 mm ,残 余 压 应 力 变 化 较 明 显,峰 值 从 57.3MPa24.9 MPa,118.0MPa112.3 MPa;,102.5 MPa,,热应 ,的作,,够提高板材 之 间 的 接 触 质 量,使 得 焊 接 接 头 受 到 的机械应力降低,焊后变形量减小,从而降低焊缝 表面残余拉应力[21]

4

(1)用顺序热力耦合的建模方法,得到镁-异种合金搅拌摩擦搭接焊的最优参数为焊接速度 60mm·min-1搭接量60mm搅拌头转速1400r· min-1;在不同焊接工艺下,热循线和残的模拟结果与试验结果相吻合,7.5%8.4%,值模(2)-表面 垂直于焊缝方向的纵向残余应力呈典型的 M 形分 ,前进侧的残余应力峰值高于后退侧;最大残余应 力出现在焊缝末端的搭接面处,最优焊接工艺下的 最大残余拉应力为137.7MPa(3)随着焊接速度的增大,接头上表面纵向残 余拉应力峰值明显增大,作用范围变窄,而残余压应 力变化不明显;随着搭接量的增加,接头上表面残 余拉应力峰值降低,作用范围变宽,残余压应力峰值 在搭接量从40 mm 增至50 mm 时呈明显下降趋

参考文献:

[1] HATTEL J H,SONNE M R,TUTUM C C.Modelling residualstressesinfrictionstirweldingofAlalloysAreview ofpossibilitiesandfuturetrends[J].TheInternationalJournal ofAdvanced Manufacturing Technology,2015,76(9/10/11/ 12):1793-1805. [2] 陈彦君,朱平,张秀娟,.6系铝合金板搅拌摩擦焊残余应力 有限元分析[J].大连交通大学学报,2019,40(6):93-97. CHEN YJ,ZHUP,ZHANGXJ,etal.Finiteelementanalysis ofresidualstressinfrictionstirweldingof6seriesaluminum alloyplate[J].JournalofDalianJiaotongUniversity,2019,40 (6):93-97. [3] 卢翔,邵良臣,韩善灵./铝异种金属搅拌摩擦焊残余应力场 分析[J].兵器材料科学与工程,2019,42(5):85-89. LU X,SHAOLC,HANSL.Analysisofresidualstressfield offrictionstir weldingonsteel/aluminum dissimilar metals [J].OrdnanceMaterialScienceandEngineering,2019,42(5): 85-89. [4] 周文静,许振波,杜柏松.铝合金搅拌摩擦焊温度场及残余应力 场研究[J].兵器材料科学与工程,2020,43(3):74-79. ZHOU WJ,XUZB,DU BS.Temperaturefieldandresidual stressfieldoffrictionstir weldingofaluminum alloy[J]. OrdnanceMaterialScienceandEngineering,2020,43(3):74- 79. [5] 史学海,吕赞,胡云瑞,.-热源辅助对 FSW 残余应力影响 的数值模拟[J].精密成形工程,2020,12(2):55-60. SHIX H,LU Z,HU Y R,etal.Numericalsimulationfor influenceofcoldandthermalsourcesonFSW residualstress [J].JournalofNetshapeFormingEngineering,2020,12(2): 55-60. [6] 张昭,谭治军,李健宇,.铝 合 金 搅 拌 摩 擦 搭 接 焊 温 度-再 结 -析出相-力 学 性 能 一 体 化 数 值 模 拟 [J].机 械 工 程 学 报, 2020,56(6):213-220. ZHANGZ,TANGZJ,LIJY.Integratednumericalsimulation oftemperature-recrystallization-precipitatemechanicalproperty infrictionstirlap weldingofaluminum alloy[J].Journalof MechanicalEngineering,2020,56(6):213-220. [7] ,,,.[J].,2018,47(13):199- 202. HU Y R,HE Z K,ZHANG L G,etal.Effectofshoulder geometryonmaterialflowandjointformationoffrictionstir lapwelding[J].HotWorkingTechnology,2018,47(13):199- 202. [8] ,宿,.数值模拟[J].金属学报,2018,54(2):265-277. WU C S,SU H,SHI L.Numerical simulation of heat generation,heattransferand materialflow infriction stir welding[J].ActaMetallurgicaSinica,2018,54(2):265-277. [9] ,赵双,,.7075研究[J].东北大学学报(自然科学版),2021,42(3):340-346. FENGY Y,ZHAO S,LIU ZS,etal.Experimentstudyand simulationforfrictionstirweldingprocessof7075aluminum alloy[J].Journal of Northeastern University (Natural Science),2021,42(3):340-346. [10] ,,,.度场[J].,2018,47(6):73-76. LIUCY,YUZK,ZHANG Y H,etal.Effectofpinshapes ontemperaturefieldinfriction stir welding[J].Welding Technology,2018,47(6):73-76. [11] 冯可.AZ31BTIG 熔钎焊的数值模 [D].:,2014. FENG K.NumericalsimulationsofdissimilarTIG welding- brazingofTitanium and magnesium alloy AZ31Blapjoints [D].Chongqing:ChongqingUniversity,2019. [12] RIAHIM,NAZARIH.Analysisoftransienttemperatureand residualthermalstressesinfrictionstirweldingofaluminum alloy6061-T6vianumericalsimulation[J].TheInternational JournalofAdvancedManufacturingTechnology,2011,55(1/ 2/3/4):143-152. [13] 才宇.轨道车辆铝镁合金异种材料 FSW 仿真研究[D].大连: 大连交通大学,2019. CAIY.SimulationondissimilarmaterialFSW with Aland Mg alloy for rail vehicles[D].Dalian:Dalian Jiaotong University,2019. [14] .设计与分析[M].北京:清华大学出版社,2005. CHEN K.Experimentaldesingabdanalysis[M].Beijing: TsinghuaUniversityPress,2005. [15] HEXC,GUFS,BALLA.Areviewofnumericalanalysisof frictionstirwelding[J].Progressin MaterialsScience,2014, 65:1-66. [16] MOHARAMI R, SATTARI-FAR I. Experimental and numericalstudyofmeasuringhighweldingresidualstresses byusingtheblind-hole-drillingtechnique[J].TheJournalof StrainAnalysisforEngineeringDesign,2008,43(3):141- 148.