基于SPH方法的搅拌摩擦焊数值模拟研究

韩泽平; 曹丽杰

doi:10.12299/jsues.23-0060

基于SPH方法的搅拌摩擦焊数值模拟研究

doi: 10.12299/jsues.23-0060

韩泽平,
曹丽杰^,

上海工程技术大学机械与汽车工程学院, 上海 201620

详细信息

作者简介:
韩泽平（1999−），男，硕士生，研究方向为搅拌摩擦焊数值模拟。E-mail：2942151857@qq.com

通讯作者:
曹丽杰（1968−），女，副教授，硕士，研究方向为轻合金激光焊接技术、摩擦学行为。E-mail：clj@sues.edu.cn

中图分类号: TG453.9
计量
- 文章访问数: 215
- HTML全文浏览量: 82
- PDF下载量: 45
- 被引次数: 0
出版历程
- 收稿日期: 2023-03-14
- 刊出日期: 2024-03-30

Numerical simulation research of friction stir welding based on SPH method

HAN Zeping,
CAO Lijie^,

School of Mechanical and Automotive Engineering, Shanghai University of Engineering Science, Shanghai 201620, China

摘要

摘要: 光滑粒子流体动力学(smoothed particle hydrodynamics, SPH)法是基于拉格朗日的一种无网格法，基于SPH方法使用ABAQUS软件进行AZ31B镁合金的搅拌摩擦焊数值模拟，通过改变搅拌头转速，研究残余应力和温度的分布。结果显示：焊缝两侧残余应力分布呈不对称，后退侧应力大于前进侧；应力曲线具有双峰特征，在焊接区域表现为拉应力；随着转速的增加，工件应力分布更加均匀；焊缝两侧的温度呈不对称分布，且前进侧高于后退侧，随着转速的增加，温度呈上升趋势；纵向温度曲线呈M形分布。模拟结果对实际焊接过程中焊接参数的选择具有一定参考意义。
- 搅拌摩擦焊 /
- 光滑粒子流体动力学 /
- 残余应力 /
- 温度场 /
- AZ31B镁合金
Abstract: Smoothed particle hydrodynamics (SPH) is a meshless method based on Lagrangian. ABAQUS software was used to simulate the friction stir welding of AZ31B magnesium alloy based on the SPH method. By changing the rotation speed of the stirring too, the distributions of residual stress and temperature were studied. The results show that the distribution of residual stress on both sides of the weld is asymmetrical, and the stress on the retreating side is greater than that on the advancing side; the stress curve has bimodal characteristics, showing tensile stress in the welding area; as the rotational speed increases, the stress distribution of the workpiece becomes more uniform; the temperature distribution on both sides of the weld is asymmetrical, and the advancing side is higher than the retreating side, and the temperature rises with the increase of the rotational speed; the longitudinal temperature curve presents an M-shaped distribution. The simulation results have certain reference significance for the selection of welding parameters in the actual welding process.
- friction stir welding /
- smoothed particle hydrodynamics (SPH) /
- residual stress /
- temperature field /
- AZ31B magnesium alloy

HTML全文

图 1 SPH原理图^[5]

Figure 1. SPH schematic diagram^[5]

下载: 全尺寸图片幻灯片

图 2 有限元网格和SPH粒子分布

Figure 2. Finite element mesh and SPH particle distribution

下载: 全尺寸图片幻灯片

图 3 网格划分

Figure 3. Meshing

下载: 全尺寸图片幻灯片

图 4 网格单元粒子转换

Figure 4. Grid unit particle conversion

下载: 全尺寸图片幻灯片

图 5 焊接速度40 mm/min时应力云图

Figure 5. Stress cloud diagram with welding speed 40mm/min

下载: 全尺寸图片幻灯片

图 6 焊接速度40 mm/min时焊接截面云图

Figure 6. Welding cross-section nephogram with welding speed 40 mm/min

下载: 全尺寸图片幻灯片

图 7 路径选取

Figure 7. Path selection

下载: 全尺寸图片幻灯片

图 8 工件上表面残余应力分布

Figure 8. Residual stress distribution on the upper surface of workpiece

下载: 全尺寸图片幻灯片

图 9 工件底面残余应力分布

Figure 9. Residual stress distribution on the bottom surface of workpiece

下载: 全尺寸图片幻灯片

图 10 温度分布云图

Figure 10. Temperature distribution cloud map

下载: 全尺寸图片幻灯片

图 11 工件上表面焊缝中心温度

Figure 11. Center temperature of weld seam on the upper surface of workpiece

下载: 全尺寸图片幻灯片

图 12 工件中部上表面横截面温度分布

Figure 12. Temperature distribution of the cross-section on the upper surface of the middle part of workpiece

下载: 全尺寸图片幻灯片

表 1 AZ31B的材料参数^[9]

Table 1. Material parameters of AZ31B^[9]

温度/ ℃	热传导系数/ (W·(m℃)⁻¹)	比热容/ (J·(kg·℃)⁻¹)	热膨胀系数/ ℃⁻¹	杨氏模量/ GPa
20	96.4	1050		40.2
75				37.3
100	101	1130	2.64E-05	34.3
125				30.9
150				30.4
200	105	1170	2.70E-05	29.4
250				27.5
300	109	1210	2.79E-05
350	113	1260

下载: 导出CSV

表 2 AZ31B的Johnson-Cook本构模型参数^[9]

Table 2. Johnson-Cook constitutive model parameters of AZ31B^[9]

材料	A /MPa	B/MPa	C	n	m	T_melt /℃	T_ref /℃
AZ31B	279.827	159	0.013	0.327	1.573	650	22

下载: 导出CSV

参考文献(15)

[1]	FRASER K, ST-GEORGES L, KISS L I. A mesh-free solid-mechanics approach for simulating the friction stir-welding process. ISHAK M. Joining Technologies[M]. London: IntechOpen, 2016. https://api.intechopen.com/chapter/pdf-download/51396.pdf.
[2]	董添文, 柳和生, 黄益宾, 等. 光滑粒子流体动力学的研究进展及应用[J] . 上饶师范学院学报,2014,34(6):28 − 38.
[3]	LIBERSKY L D, PETSCHEK A G. Smooth particle hydrodynamics with strength of materials[C]//Advances in the Free-Lagrange Method Including Contributions on Adaptive Gridding and the Smooth Particle Hydrodynamics Method: Lecture notes in physics. Berlin: Springer, 1991. DOI: 10.1007/3-540-54960-9_58.
[4]	CERVERA M, CHIUMENTI M, VALVERDE Q, et al. Mixed linear/linear simplicial elements for incompressible elasticity and plasticity[J] . Computer Methods in Applied Mechanics and Engineering,2003,192(49/50):5249 − 5263. doi: 10.1016/j.cma.2003.07.007
[5]	FRASER K, St-GEORGES L, KISS L I. Numerical simulation of bobbin tool friction stir welding[C]//Proceedings of the 10th International Friction Stir Welding Symposium. Beijing: the Joining and Welding Research Institute, 2014.
[6]	TIMESLI A, MOUFKI A, ZAHROUNI H, et al. Numerical model based on SPH method to simulate Friction Stir Welding[C]//10e colloque national en calcul des structures. 2011: Clé USB.
[7]	LIU Z, XIU L, WU J, et al. Numerical simulation on residual stress eliminated by shot peening using SPH method[J] . Fusion Engineering and Design,2019,147:111231. doi: 10.1016/j.fusengdes.2019.06.004
[8]	ZAMANI S M M, BEHDINAN K, MOHAMMADPOUR A, et al. Friction stir welding of Al-SiC composite sheets: a numerical simulation of residual stresses[J] . The International Journal of Advanced Manufacturing Technology,2021,116(11):3717 − 3729.
[9]	侯增. 2024铝合金/AZ31B镁合金异种金属搅拌摩擦焊温度场的数值模拟[D]. 北京: 中国石油大学, 2017.
[10]	马核, 田志杰, 熊林玉, 等. 2A14-T6铝合金搅拌摩擦焊温度场及黏流层数值模拟分析[J] . 航空制造技术,2018,61(8):55 − 61.
[11]	武凯, 贾贺鹏, 孙宇, 等. 搅拌摩擦焊技术的研究进展[J] . 机械制造与自动化,2020,49(6):1 − 9.
[12]	殷鹏飞, 张蓉, 熊江涛, 等. 搅拌摩擦焊准稳态温度场数值模拟[J] . 西北工业大学学报,2012,30(4):622 − 627.
[13]	陈彦君. 基于ABAQUS的搅拌摩擦焊仿真模拟分析[D]. 大连: 大连交通大学, 2019.
[14]	KISHTA E E, FARID A, DARRAS B. Nonlinear finite element simulation of friction stir processing of marine grade 5083 aluminum alloy[J] . Engineering Transactions,2014,62(4):313 − 328.
[15]	HASAN A F. CFD modelling of friction stir welding (FSW) process of AZ31 magnesium alloy using volume of fluid method[J] . Journal of Materials Research and Technology,2019,8(2):1819 − 1827. doi: 10.1016/j.jmrt.2018.11.016