Simulation and experimental research on friction stir welding of large thickness aluminum alloys
-
摘要: 针对10 mm大厚度7050 铝合金板搅拌摩擦焊接需求,建立了任意拉格朗日–欧拉模型,结合仿真分析和实验研究,分析了转速、焊接速度及焊接时间下接头温度场分布规律。实测温度和特征点的热循环曲线吻合较好,说明数值模拟可以用来预测温度场的分布规律。结果表明:不同焊接参数下焊缝横截面温度场分布均呈碗状。各位置点的峰值温度随着转速增大及焊接速度减小而升高,其中转速对峰值温度影响更显著。最后,实验验证3组工艺参数下接头的力学性能,实验表明,转速为800 r/min、焊接速度为1.6 mm/s时,接头的力学性能较好。本研究方法为大厚度轻质合金搅拌摩擦焊接工艺的确定提供了参考和借鉴。Abstract: In order to meet the requirement of friction stir welding joint for 7050 aluminum alloy plate with 10 mm thickness, an arbitrary Lagrange-Olah model was established, which was combined with simulation analysis and experimental research, the distribution law of joint temperature field under rotating speed, welding speed and welding time was analyzed. The thermal cycle curves of the measured temperature and the characteristic points are in good agreement, indicating that the numerical simulation can be used to predict the distribution pattern of the temperature field. The results show that the temperature field distribution of the weld cross-section under different welding parameters are bowl-shaped. The peak temperature at each location point increases with increasing rotational speed and decreasing welding speed, where the rotational speed has a more significant effect on the peak temperature. Finally, the experimental verification of the mechanical properties of the joint under three groups of process parameters, found that 800 r/min, 1.6 mm/s group of parameters under the better tensile properties, the study method for the determination of large thickness of light alloy stir friction welding process provides a reference.
-
激光选区熔化(selective laser melting, SLM)是一种使用高能聚束激光逐点扫描熔化金属粉末,逐层凝固累计成形的粉末床成型技术,因具有高自由度、高材料利用率和高生产效率等优点,所增材金属部件成为机械制造行业应用发展的新趋势[1−2]。Inconel 718是高温镍基合金中最常用的材料,因其良好的高温力学性能、抗蠕变性和耐腐蚀性,被誉为航空航天、石油化工、核工业等高端领域最具潜力的材料之一,制造应用于航空航天涡轮发动机、燃气轮机、工作叶片和燃烧室等重要零部件[3−4]。随着航空技术的快速发展,激光选区熔化Inconel 718(SLM Inconel 718)复杂结构零部件的需求日益增长[5]。现阶段,由于SLM Inconel 718成形过程存在粉末黏附、球化飞溅和台阶效应,其表面成形质量[6]、力学性能[7]及显微组织结构[8]等方面与传统Inconel 718合金存在较大差异,且激光成形过程中因热变形而出现尺寸精度低、形状精度低及表面质量差的问题,通常难以满足精密零部件的精度要求,导致无法直接应用[9−11]。
近年来,国内外学者为了解决上述问题,对SLM Inconel 718材料的后处理工艺进行了深入研究。Zhang B C等[12]对SLM Inconel 718零件进行不同持续时间电化学抛光(ECP)的表面形貌和粗糙度演变研究,发现在ECP 5 min的条件下,表面粗糙度从6.05 μm降低到3.66 μm。Kaynak Y等[13]研究干冷空气条件下的精加工车削(FM)、拖动精加工(DF)和振动精加工(VSF)对SLM Inconel 718合金表面和亚表面性能的影响,发现直接成型件经过FM加工后表面粗糙度降低96%,DF和VSF加工后表面粗糙度降低到73%~88%。Ćwikła Michał等[14]研究了在空气和氩气的分别作用下激光抛光(LP)对SLM Inconel 718的表面质量影响,表明在空气和氩气不同条件下,LP使得零件平均粗糙度分别降低82.8%和87.9%。Ji H S等[15]研究SLM Inconel 718+微铣削显微组织和可加工性,发现微铣削使SLM Inconel 718组织出现轻微晶体粗化现象,且残余应力随主轴转速的增加而减小,随进给速度的增加而增大,刀具磨损随主轴转速和进给速度的增大而增大。
上述研究采用不同后处理工艺分析,改善了SLM Inconel 718的表面质量和力学性能,但目前对SLM Inconel 718成形件切削后处理机理和工艺优化的研究鲜见。因此开展相应高效切削加工研究,优化切削工艺,降低刀具磨损,提高表面质量和加工效率对SLM Inconel 718成形件的广泛应用具有重要推动作用。
本研究选择铣削减材后处理工艺,设计四因素四水平正交试验,分析不同条件下的切削力、表面粗糙度的变化规律,同时基于灰色关联分析、逐步回归法和粒子群算法实现以切削力、表面粗糙度和材料去除率为性能指标的多目标参数优化,建立灰色关联度与工艺参数的二阶回归预测模型,获得使得切削力最小、表面粗糙度最小和材料去除率最大的最优参数组合,为增减材Inconel 718加工提供一定的理论参考。
1. 试验
1.1 SLM增材试验
试验设备选用图1所示的EP-M150金属3D打印机,试验材料选用由气雾化法制备的Inconel 718高温合金球形金属粉末,粒径范围为15~53 μm,化学成分见表1。SLM成型工艺参数[16]:激光功率295 W,扫描速度950 mm/min,扫描间距0.1 mm,层厚0.04 mm,制备得到50 mm×30 mm×20 mm的块状样品,如图2所示。
表 1 Inconel 718粉末化学成分(%) Ni Cr Nb Mo Ti Al Co Cu C B Fe 51.08 19.24 5.28 3.13 0.87 0.6 0.19 0.1 0.05 0.006 余 1.2 铣削减材试验
影响增材件铣削加工结果的因素较多,综合考虑后采用 $ {L}_{16}{(4}^{5}) $正交表,设计四因素四水平正交试验,选择数值见表2。试验设备选用DMU50五轴联动加工中心,铣削加工刀具采用厦门鸿鹭联创公司的TiAlN涂层硬质合金立铣刀,所属系列S10-R4-08010,参数见表3。试验采用干式顺铣方式加工,使用瑞士奇石乐Kistler 9527B压力传感器和Kistler 5070A多通道电荷放大器以及数据采集系统5697A1型对三向切削力实时监测,Dyno ware软件记录变化曲线。使用TR200表面粗糙度测量仪测量工件铣削前后的表面粗糙度,使用基恩士VHX-7000超景深显微镜观察铣削前后表面形貌、铣削加工试验与测量数据现场如图3所示。
表 2 试验因素水平表水平 铣削速度
$ {v}_{c} $/(m/min)铣削深度
$ {a}_{p} $/mm铣削宽度
$ {a}_{e} $/mm每齿进给量
$ {f}_{{\textit{z}}} $/(mm/z)1 30 0.25 2.0 0.025 2 60 0.35 2.5 0.050 3 90 0.45 3.0 0.075 4 120 0.55 3.5 0.100 表 3 刀具参数L/mm d/mm $ {L}_{c} $/mm $ {d}_{1} $/mm $ {\gamma }_{0} $/(°) $ {\alpha }_{0} $/(°) $ \beta $/(°) 60 8 20 8 9.37 15.36 35 注:L为总长;d为刃径;$ {L}_{c} $为刃长;$ {d}_{1} $为柄径;$ {\gamma }_{0} $为前角;$ {\alpha }_{0} $为后角;$ \beta $为螺旋角。 2. 结果与分析
为了获得科学观测值,在测力仪软件数据中取实时值相对稳定位置的数据均值作为X、Y、Z方向的铣削力数值,并根据式(1)求其铣削合力F。为提高加工效率引入材料去除率MRR,MRR由式(2)计算得出。表面粗糙度$ {R\mathrm{a}}_{} $选用五点测量法,数据中去掉最大值和最小值后再取平均进行记录。正交试验设计结果见表4,铣削前测得表面粗糙度为14.804 μm。
表 4 正交试验设计与结果编号 $ {v}_{c} $/
(m/min)$ {a}_{p} $/mm $ {a}_{e} $/mm $ {f}_{{\textit{z}}} $/
(mm/z)F/N Ra/μm MRR/
(cm3/min)1 30 0.25 2.0 0.025 117.63 0.325 0.06 2 30 0.35 2.5 0.050 212.77 0.705 0.21 3 30 0.45 3.0 0.075 273.24 1.337 0.48 4 30 0.55 3.5 0.100 340.79 1.823 0.92 5 60 0.25 2.5 0.075 207.45 0.514 0.45 6 60 0.35 2.0 0.100 299.89 0.763 0.67 7 60 0.45 3.5 0.025 225.20 0.541 0.38 8 60 0.55 3.0 0.050 234.71 0.996 0.79 9 90 0.25 3.0 0.100 282.25 1.113 1.07 10 90 0.35 3.5 0.075 294.84 0.705 1.32 11 90 0.45 2.0 0.050 220.77 0.729 0.64 12 90 0.55 2.5 0.025 208.49 0.621 0.49 13 120 0.25 3.5 0.050 271.71 0.931 0.84 14 120 0.35 3.0 0.025 252.98 0.529 0.50 15 120 0.45 2.5 0.100 346.93 1.628 2.15 16 120 0.55 2.0 0.075 295.47 0.890 1.58 $$ F=\sqrt{F_x^2+F_y^2+F_{\text{z}}^2} $$ (1) $$ MRR=\frac{\mathit{{v}_c}a_p{\textit{z}}a_ef_{\textit{z}}}{\text{π}{d}} $$ (2) 式中:$ {F}_{x} $、$ {F}_{y} $、$ {F}_{{\textit{z}}} $分别为X、Y、Z方向的铣削分力;z为铣刀齿数,本试验z取4。
2.1 铣削参数对响应变量的影响分析
根据表4数据进行极差法分析,结果见表5和表6,绘制不同水平下响应变量与铣削参数的影响变化图,结果如图4~图7所示。
表 5 铣削合力F极差表因素 $ {k}_{1} $ $ {k}_{2} $ $ {k}_{3} $ $ {k}_{4} $ 极差R 主次顺序 $ {v}_{c} $ 236.1 241.8 251.6 291.8 55.7 2 $ {a}_{p} $ 219.8 265.1 266.5 269.9 50.1 3 $ {a}_{e} $ 233.4 243.9 260.8 283.1 49.7 4 $ {f}_{{\textit{z}}} $ 263.3 286.2 384.8 395.3 132.0 1 表 6 表面粗糙度Ra极差表因素 $ {k}_{1} $ $ {k}_{2} $ $ {k}_{3} $ $ {k}_{4} $ 极差R 主次顺序 $ {v}_{c} $ 1.048 0.704 0.792 1.002 0.344 3 $ {a}_{p} $ 0.721 0.676 1.066 1.083 0.407 2 $ {a}_{e} $ 0.677 0.875 0.994 1.000 0.323 4 $ {f}_{{\textit{z}}} $ 0.504 0.840 0.862 1.339 0.835 1 注:ki为该因素下第i水平所对应试验结果求和平均值,极差R为ki中最大值与最小值之差。 由表5和表6可知,铣削参数对铣削合力F的影响主次顺序依次是$ {f}_{{\textit{z}}}、{v}_{c} $、$ {a}_{p} $、$ {a}_{e} $。铣削参数对表面粗糙度Ra的影响主次顺序依次是$ {f}_{{\textit{z}}}、{a}_{p}、{v}_{c}{、a}_{e} $。
由图4a可知,在试验四水平内,随着铣削速度的增加,铣削速度与铣削合力F呈正相关,但增加速率不同,增加速率在30~60 m/min过程最小,增加速率在90~100 m/min过程最大。这是因为铣削速度的增加伴随着单位时间内切削厚度的增大,铣削力呈上升趋势。由4b可知,表面粗糙度值随着铣削速度的增大呈现快速下降后递增的趋势,铣削速度为60 m/min时,表面粗糙度值降到最小。这是因为持续高温环境对材料有一定程度的软化,工件与刀具之间的摩擦系数下降,材料表面粗糙度降低,而当速度持续增加时,刀具的磨损和与工件之间的颤振现象会导致表面粗糙度值增大,所以增材件的铣削速度在一定范围内增大可以改善表面质量,但不宜过大。
由图5a可知,随着铣削深度的增加,铣削力也随之增加,铣削深度增加至0.45 mm后上升速率变得缓慢,这主要因为开始加工时,增材件表面硬化,切削层横截面增大,材料去除所需能量增大,铣削力增加,而随铣削深度的增加,加工过程可以避开表面硬化层,且铣削温度上升影响材料软化程度,铣削力增加速率变慢。由图5b可知,表面粗糙度的大小随着铣削深度的增加呈现缓慢下降后急剧上升的趋势,最后速率变慢。这主要因为开始的铣削深度较小,导致刀刃与工件表面加工震动较大,表面质量较差;当铣削深度增至0.35 mm时,加工趋于平稳,表面质量良好,表面粗糙度值降低;当深度继续增加时,刀具和工件之间挤压程度增大,工件塑性变形程度变大,表面粗糙度值增加。所以,增材件的表面质量优化应尽量选择较小的铣削深度。
由图6可知,随着铣削宽度的增大,铣削合力整体呈上升趋势。而当铣削宽度小于3.0 mm时,表面粗糙度值与其呈正相关,当铣削宽度大于3.0 mm时,表面粗糙度值上升趋势明显变缓。这是因为宽度增加使刀刃与工件平面的接触面积增大,单位时间内材料去除率增大,摩擦力增大,铣削力和表面粗糙度值也随之上升,当宽度增大到一定范围时,表面粗糙度受其影响明显变小。
由图7可知,随着每齿进给量的增加,铣削力与表面粗糙度值都呈现上升趋势。这主要因为刀具每齿切削间距增加使得瞬时材料去除量增加,铣削力变大,刀具与工件之间的振动变大,导致表面质量变差。因此,在对增材件的切削过程中应选择较小的每齿进给量,以获得更好的表面质量从而达到使用标准。
2.2 铣削工艺参数对表面形貌的影响分析
通过基恩士VHX-7000超景深显微镜所拍摄的铣削前后表面形貌如图8所示。在图8a中,SLM成形件存在扫描熔道痕迹和球化物,导致表面产生较大波峰波谷,粗糙度测量值过大。而从图8b中表面形貌图可见,微铣削后表面存在铣削刀纹,且去除了SLM技术形成的表层材料使得孔隙减少,同表4表面粗糙度结果综合分析可得,第1组工艺参数下加工效果更为理想,此时表面粗糙度值最低,其表面粗糙度Ra最大下降97.80%,整体平均下降94.03%。
3. 工艺参数优化
提高零件的加工质量和效率是增材件应用领域内的重点关注问题,解决这一问题的最主要方法是选择合适的加工参数。本文为了探究“SLM Inconel 718+铣削”参数的最优组合,基于灰色关联分析(GRA)、逐步回归法(SR)和粒子群算法(PSO)求解。
3.1 GRA分析
本文以铣削合力F、表面粗糙度Ra和材料去除率MRR为响应变量,对其进行多目标优化,而灰色关联分析法[17−18] 可以解决优化过程的复杂性问题,将多目标优化转化为单目标优化。方法步骤如下。
(1)评价指标矩阵
假设评价系统有n个试验方案和p个评价指标,则评价指标矩阵X为
$$ {\boldsymbol{X}} = \left[ \begin{array}{*{20}{c}} {x_{11}} & {x_{12}} & \cdots & {x_{1p}} \\ {x_{21}} & {x_{22}} & \cdots & {x_{2p}} \\ \vdots & \vdots & \ddots & \vdots \\ {x_{n1}} & {x_{n2}} & \cdots & {x_{np}} \\ \end{array} \right] $$ (2)数据归一化
对每个响应变量的试验数据进行归一化处理,使其介于0和1之间,减少可变性。公式如下:
$$ x_i^*(k) = \frac{{\max x_i^0(k) - x_i^0(k)}}{{\max x_i^0(k) - \min x_i^0(k)}} $$ (3) $$ x_i^*(k) = \frac{{x_i^0(k) - \min x_i^0(k)}}{{\max x_i^0(k) - \min x_i^0(k)}} $$ (4) 式中:$ x_i^*(k) $为归一化数值;$ x_i^0(k) $为原始试验数据;$ \max x_i^0(k) $为第i个数列中所有试验的最大水平;$ \min x_i^0(k) $为第i个数列中所有试验的最小水平。
因为铣削合力和表面粗糙度越小越好,所以选择式(3);而对于材料去除率来讲,越大越好,因此选择式(4)。
(3)灰色关联系数(GRC)
将数据归一化后,代入公式计算灰色关联系数。公式如下:
$$ {\zeta _{{i}}}(k) = \frac{{\mathop {\min }\limits_{\forall j \in i} \left\| {\mathop {\min }\limits_{\forall k} x_0^*({{k}}) - x_i^*({{k}})} \right\| + \rho \cdot \mathop {\max }\limits_{\forall j \in i} \left\| {\mathop {\max }\limits_{\forall k} x_0^*({{k}}) - x_i^*({{k}})} \right\|}}{{\left\| {x_0^*({{k}}) - x_i^*({{k}})} \right\| + \rho \cdot \mathop {\max }\limits_{\forall j \in i} \left\| {\mathop {\max }\limits_{\forall k} x_0^*({{k}}) - x_i^*({{k}})} \right\|}} $$ (5) 式中:$ {\zeta _{{i}}}(k) $为灰色关联系数;$ \rho $为分辨系数,主要作用是提高关联系数之间的差异显著性,$ \rho \in \left[ {0,1} \right] $,本次计算取0.5;$ {x}_{0}^{*}(k)和{x}_{{i}}^{*}(k) $分别为参考序列(原始)和比较序列(归一化处理)。
(4)计算灰色关联度(GRG)
取每个性能特征GRC的平均值来计算GRG,GRG值越大,对应的加工组合越接近响应变量的最优值。公式如下:
$$ {\zeta _{{i}}} = \frac{1}{n}\sum\limits_{k = 1}^m {{\zeta _{{i}}}} (k) $$ (6) 式中:$ {\zeta _{{i}}} $为第i次试验的灰色关联度;m为响应变量的数量。通过式(3)~式(6)计算,结果见表7。
表 7 灰色关联分析数据结果序号 归一化处理 GRC GRG F Ra MRR F Ra MRR 1 1 1 0 1 1 0.333 3 0.777 8 2 0.582 1 0.746 3 0.071 4 0.546 5 0.663 4 0.350 0 0.520 0 3 0.321 4 0.324 4 0.202 9 0.424 2 0.425 3 0.385 5 0.411 7 4 0.026 8 0 0.411 5 0.339 4 0.333 3 0.459 4 0.377 4 5 0.608 3 0.873 8 0.185 7 0.560 7 0.798 5 0.380 4 0.579 9 6 0.205 1 0.707 6 0.291 4 0.386 1 0.631 0 0.413 7 0.476 9 7 0.530 9 0.855 8 0.151 4 0.515 9 0.776 2 0.370 8 0.554 3 8 0.489 4 0.552 1 0.348 5 0.494 8 0.527 5 0.434 2 0.485 5 9 0.282 1 0.474 0 0.485 7 0.410 5 0.487 3 0.492 9 0.463 6 10 0.227 2 0.746 3 0.601 4 0.392 8 0.663 4 0.556 4 0.537 6 11 0.550 2 0.730 3 0.280 0 0.526 4 0.649 6 0.409 8 0.528 6 12 0.603 8 0.802 4 0.207 1 0.557 9 0.716 7 0.386 7 0.583 8 13 0.328 0 0.595 5 0.371 4 0.426 6 0.552 8 0.443 0 0.474 1 14 0.409 7 0.863 8 0.211 4 0.458 6 0.785 9 0.388 0 0.544 2 15 0 0.110 1 1 0.333 3 0.359 8 1 0.564 4 16 0.224 4 0.622 8 0.725 7 0.392 0 0.570 0 0.645 8 0.535 9 3.2 SR法建立预测模型
为进行后续优化问题,需建立工艺参数与灰色关联度值之间的映射关系。采用SR法[19]建立数学模型,表达式如下:
$$ \hat y = {\beta _0} + \sum\limits_{i = 1}^4 {{\beta _0}} {x_i} + \sum\limits_{i = 1}^4 {\sum\limits_{j = i + 1}^4 {{\beta _{ij}}} {x_i}} {x_j} + \sum\limits_{i = 1}^4 {{\beta _{ij}}} x_i^2 $$ (7) 式中:$ \hat y $为GRG的预测值;$ \beta $为数学模型的系数;$\displaystyle\sum\limits_{i = 1}^4 {{\beta _0}} {x_i} $为线性效应;$\displaystyle\sum\limits_{i = 1}^4 {\displaystyle\sum\limits_{j = i + 1}^4 {{\beta _{ij}}} {x_i}} {x_j} $为交互性效应;$\displaystyle\sum\limits_{i = 1}^4 {{\beta _{ij}}} x_i^2 $为二次效应。
运用Matlab软件编写逐步回归法程序计算,得到其二阶多项式回归预测模型结果如下:
$$ \begin{split} {GRG}= & 1.0802-2.3989a_p+0.2954a_e-12.3823f_{{\textit{z}}}+ \\ & 0.0146v_ca_p+0.0176v_cf_{{\textit{z}}}+17.0944a_pf_{{\textit{z}}}- \\ & 0.0001v_c^2-0.0506a_e^2+23.9053f_{{\textit{z}}}^2 \end{split} $$ (8) 其所得回归模型P值小于0.05,模型泛化性较高,整体影响显著。为进一步验证模型的准确性,计算得到其残差个案次序,结果如图9所示,观察发现无异常点,模型残差符合正态分布规律,不存在过拟合现象。其预测值与实测值的关系如图10所示,误差率均小于5%,说明其预测形态良好,可作为灰色关联度分析的预测模型。
3.3 粒子群算法优化
基于GRG预测模型进行参数优化,是非线性函数寻找极值的问题,传统优化算法容易陷入局部最优解,而PSO[20]是基于群智能理论的优化问题,全局和局部搜索能力强,故选择PSO。它主要包括对随机粒子群的初始化操作和逐次迭代逼近最优解,其中把每一次迭代过程所计算得到的各粒子目标函数值作为适应度值。设置粒子群群体个数为500,最大迭代次数2 500,粒子加速常数均为2。以试验参数水平范围为该优化问题的约束条件,即
$$ \mathrm{s.t.}\left\{\begin{array}{l}30\leqslant v_c\leqslant120 \\ 0.25\leqslant a_p\leqslant0.55 \\ 2\leqslant a_e\leqslant3.5 \\ 0.025\leqslant f_{{\textit{z}}}\leqslant0.1\end{array}\right. $$ 通过Matlab工具实现程序编写与运行,算法优化结果如图11所示。经过2 500次迭代,目标函数适应度值变化趋于稳定,此时GRG的最优解为0.817 3,高于表5中GRG的最大值0.777 8,说明在此值下的响应变量整体最优,其所对应的工艺参数组合为$ {v}_{c} $=46.96 m/min,$ {a}_{p} $=0.25 mm,$ {a}_{e} $=2.90 mm,$ {f}_{{\textit{z}}} $=0.025 mm/z。
4. 结语
(1)采用SLM增材工艺和铣削减材工艺,设计四因素四水平正交试验,得到16组不同工艺参数下的铣削合力、表面粗糙度值,通过极差分析得到:对F影响的重要性依次为$ {f}_{{\textit{z}}} > {v}_{c} > {a}_{p} > {a}_{e} $;对Ra影响的重要性依次为$ {f}_{{\textit{z}}} > {a}_{p} > {v}_{c}{ > a}_{e} $。同时表明了F、Ra在试验四水平下随$ {v}_{c}、{a}_{p}、{a}_{e}、{f}_{{\textit{z}}} $的影响规律,并通过对比铣削前后表面粗糙度值和表面形貌分析,验证了微铣削工艺可以改善增材件表面完整性,其Ra最大下降97.80%,平均下降94.03%。
(2)在响应变量中引入材料去除率MRR,通过灰色关联分析法将多目标优化问题转化为单目标优化问题,利用逐步回归法建立灰色关联度值与工艺参数的二阶多项式回归预测模型,并通过残差图分析和实测-预测值对比图观察结果,说明了该模型具有较高的预测精度,其最大误差不超过5%。
(3)基于灰色关联度与工艺参数的二阶非线性回归模型采用PSO算法优化,结果表明,GRG的最优解为0.817 3,所对应的工艺参数组合为$ {v}_{c} $=46.96 m/min,$ {a}_{p} $=0.25 mm,$ {a}_{e} $=2.90 mm,$ {f}_{{\textit{z}}} $=0.025 mm/z。
-
表 1 7050铝合金的热物性参数
温度/ ℃ 25 40 100 200 300 400 500 热导率/
(W/(m·K))164.46 166.45 172.63 178.75 181.69 182.54 比热容/
(J/(kg·℃))894.55 925.4 968.85 1 009.95 1 052.05 1 104.1 -
[1] 管仁国,娄花芬,黄晖,等. 铝合金材料发展现状, 趋势及展望[J]. 中国工程科学,2020,22(5):68-75. [2] Xu W F,Luo Y X,Zhang W,et al. Comparative study on local and global mechanical properties of bobbin tool and conventional friction stir welded 7085-T7452 aluminum thick plate[J]. Journal of Materials Science & Technology,2018,34(1):173-184.
[3] 庞嘉尧,余乐,杨宏. 搅拌摩擦焊中铝合金板材温度场的研究现状[J]. 机械管理开发,2021,36(6):277-280. [4] 王磊,谢里阳,张丹,等. 搅拌摩擦焊接过程温度场动态仿真[J]. 东北大学学报:自然科学版,2008(7):1025-1028. [5] 顾乃建. 搅拌摩擦焊温度场及流场数值模拟[D]. 大连: 大连交通大学, 2018 [6] 武凯,贾贺鹏,孙宇,等. 搅拌摩擦焊技术的研究进展[J]. 机械制造与自动化,2020,49(6):1-9. [7] 刘占军,姜秀萍,李亮,等. 焊接工艺参数对LF5铝合金搅拌摩擦焊温度场的影响[J]. 热加工工艺,2012,41(23):145-147. [8] 崔超,吴志生,柴斐,等. 5A06铝合金搅拌摩擦焊温度场数值模拟[J]. 焊接技术,2017,46(3):12-15,97-98. [9] 朱海,孙朝伟,孙金睿,等. 2024铝合金搅拌摩擦焊工艺参数对焊接质量的影响研究[J]. 热加工工艺,2019,48(23):159-162. [10] Mohan R,Rajesh N R,Kumar S S. Finite element modeling for maximum temperature in friction stir welding of AA 1100 and optimization of process parameter by taguchi method[J]. IJRET:International Journal of Research in Engineering and Technology,2014,3(5):728-733.
[11] Lotfi A H,Nourouzi S. Effect of welding parameters on microstructure, thermal, and mechanical properties of friction-stir welded joints of AA7075-T6 aluminum alloy[J]. Metallurgical and Materials Transactions A,2014,45(6):2792-2807. DOI: 10.1007/s11661-014-2235-z
[12] Ma F J, Fu D Y, Liu Y, et al. The numerical simulation of temperature field in friction stir welding of 7075 aluminium alloy[C].IOP Conference Series: Materials Science and Engineering. IOP Publishing, 2020, 751(1): 012081.
[13] 冯莹莹,赵双,刘照松,等. 7075铝合金搅拌摩擦焊模拟与实验研究[J]. 东北大学学报:自然科学版,2021,42(3):340-346. [14] 刘杰. 搅拌摩擦焊工艺参数对厚板铝合金接头组织和性能的影响[C].中国机械工程学会焊接学会第十八次全国焊接学术会议论文集——S01压力焊, 2013: 112-116. [15] 汪洪峰,汪建利,宋娓娓,等. 7022铝合金搅拌摩擦焊工艺参数的优化[J]. 热加工工艺,2013,42(9):156-159. [16] Kumar D A,Biswas P,Tikader S,et al. A study on friction stir welding of 12mm thick aluminum alloy plates[J]. Journal of Marine Science and Application,2013,12(4):493-499. DOI: 10.1007/s11804-013-1221-y
[17] 江旭东,黄俊,周琦,等. 铝-铜异种材料对接搅拌摩擦焊温度场数值模拟[J]. 焊接学报,2018,39(3):16-20,129. [18] 杨明鄂,李治国. 车身搅拌摩擦焊温度场仿真分析及验验证[J]. 机械设计与制造,2021(10):165-168. [19] 夏弘扬. 6082-T6铝合金高转速搅拌摩擦焊工艺研究及温度场模拟[D]. 镇江: 江苏科技大学, 2020. [20] 卢晓红, 乔金辉, 周宇, 等. 搅拌摩擦焊温度场研究进展[J/OL]. 吉林大学学报: 工学版, 2022: 1-17. [21] Dong X W,Li X F,Zuo D W,et al. Numerical simulation of temperature field in the full friction stir welding process for 7022 aluminum alloy[J]. Materials for Mechanical Engineering,2012,36(10):92-96.
[22] 刘其鹏,顾乃建,刘泽,等. AA6061-T6板材搅拌摩擦焊温度场仿真[J]. 大连交通大学学报,2018,39(3):80-85. DOI: 10.13291/j.cnki.djdxac.2018.03.016 [23] 陈彦君. 基于ABAQUS的搅拌摩擦焊仿真模拟分析[D]. 大连: 大连交通大学, 2019.