• CN:11-2187/TH
  • ISSN:0577-6686

机械工程学报 ›› 2019, Vol. 55 ›› Issue (20): 168-177.doi: 10.3901/JME.2019.20.168

• 材料科学与工程 • 上一篇    

高强低合金钢焊接过程多物理场耦合数值模拟

孙玉杰1, 史清宇2, 臧勇3, 张锁怀4, 崔青春1   

  1. 1. 西北机电工程研究所 咸阳 712099;
    2. 清华大学机械工程系 北京 100084;
    3. 北京科技大学机械工程学院 北京 100083;
    4. 上海应用技术大学机械工程学院 上海 201418
  • 收稿日期:2019-01-27 修回日期:2019-07-10 发布日期:2020-01-07
  • 通讯作者: 孙玉杰(通信作者),男,1981年出生,博士,高级工程师。主要研究方向为焊接过程数值模拟。E-mail:sunyujiegood@eyou.com
  • 作者简介:史清宇,男,1971年出生,博士,教授,博士研究生导师。主要研究方向为搅拌摩擦焊接与加工,焊接力学及模拟仿真。E-mail:shqy@tsinghua.edu.cn;臧勇,男,1963年出生,教授,博士研究生导师。主要研究方向为塑性加工工艺及装备。E-mail:yzang@ustb.edu.cn;张锁怀,男,1962年出生,博士,教授。主要研究方向为机械系统动力学,数字化设计。E-mail:zhangsuohuai@263.net

Numerical Simulation of Multi-physical Coupling of Welding Process for High Strength Low Alloy Steel

SUN Yujie1, SHI Qingyu2, ZANG Yong3, ZHANG Suohuai4, CUI Qingchun1   

  1. 1. Northwest Institute of Mechanical & Electrical Engineering, Xianyang 712099;
    2. Department of Mechanical Engineering, Tsinghua University, Beijing 100084;
    3. School of Mechanical Engineering, University of Science and Technology Beijing, Beijing 100083;
    4. School of Mechanical Engineering, Shanghai Institute of Technology, Shanghai 201418
  • Received:2019-01-27 Revised:2019-07-10 Published:2020-01-07

摘要: 目前高强低合金钢焊接数值模拟中,采用热-力耦合分析时,忽略固态相变效应,残余应力模拟值与试验测量值误差较大。为提高焊接数值模拟精度,根据多场耦合关系,基于传热学、固态相变理论和连续介质力学,建立焊接过程多物理场耦合本构方程,并通过子程序将其嵌入到通用隐式有限元程序中。采用数值模拟与试验分析的方法研究高强低合金钢小试样的自由膨胀试验、相变塑性试验及平板焊接试验应力及各应变分量的演变。研究结果表明:固态相变体积变化引起的相变应变对残余应力有显著影响,不但改变了残余应力的大小,甚至改变了残余应力的符号,考虑相变塑性应变时会降低应力的水平。残余应力改变程度与相变程度有关系:完全相变区影响最大,部分相变区次之,未发生相变区最小。相变应变和相变塑性应变最终大小相当。研究方法为深入了解高强低合金钢焊接过程和焊接工艺优化提供了理论基础。

关键词: 多物理场耦合本构方程, 完全隐式积分算法, 数值模拟, 残余应力, 应变分量

Abstract: At present, the errors between welding residual stress simulation results and experimental measurements are large for high strength low alloy steel when adopting coupled temperature-displacement analysis and omitting solid-state phase transformation effect. Based on the multi-physical coupling relationships, a multi-physical coupling constitute equation is established by using heat transfer, solid-state phase transformation theory and continuum mechanics to improve welding numerical simulation accuracy and then the equation is implemented into the general purpose implicit finite element program via user material subroutine. The evolution of stress and strain components of free dilation test, transformation induced plasticity test and flat plate welding test are investigated by means of both numerical simulations and experiments. The results suggest that transformation strain caused by volume change has significant effect on welding residual stress, which not only changes the magnitude of residual stress and even changes the sign of residual stress. Transformation induced plasticity strain lowers the magnitude of residual stress. The extent of welding residual stress change has a close relationship with the degree of phase transformation, fully-transformed region is larger than partially-transformed region, and partially-transformed region is larger than untransformed region. The final transformation strain is approximate to transformation induced plasticity strain. The presented research method provides a theoretical basis for deeper understanding of welding process and optimization of welding technology.

Key words: multi-physical coupling constitute equation, fully implicit integration algorithm, numerical simulation, residual stress, strain component

中图分类号: