CN107220410A - 参量对焊接残余应力及变形影响敏感度获取方法 - Google Patents

参量对焊接残余应力及变形影响敏感度获取方法 Download PDF

Info

Publication number
CN107220410A
CN107220410A CN201710305152.6A CN201710305152A CN107220410A CN 107220410 A CN107220410 A CN 107220410A CN 201710305152 A CN201710305152 A CN 201710305152A CN 107220410 A CN107220410 A CN 107220410A
Authority
CN
China
Prior art keywords
parameter
welding
residual stress
susceptibility
performance
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201710305152.6A
Other languages
English (en)
Other versions
CN107220410B (zh
Inventor
李晓东
孟立春
孙晓红
孙晓光
王海波
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
CRRC Qingdao Sifang Co Ltd
Original Assignee
CRRC Qingdao Sifang Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by CRRC Qingdao Sifang Co Ltd filed Critical CRRC Qingdao Sifang Co Ltd
Priority to CN201710305152.6A priority Critical patent/CN107220410B/zh
Publication of CN107220410A publication Critical patent/CN107220410A/zh
Application granted granted Critical
Publication of CN107220410B publication Critical patent/CN107220410B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

本发明涉及焊接领域,公开了一种参量对焊接残余应力及变形影响敏感度获取方法,建立焊接接头焊接过程的仿真模型;选定母材材料的多个性能参数,依次实验获取各个性能参数的参量;分别以各个性能参数的参量为基准,取该参量的设定差值范围内的多个不同参量数值依次输入仿真模型;重新实验获取敏感度最高的性能参数的参量,以该性能参数的参量为基准,取该参量的设定差值范围内的多个不同参量数值依次输入仿真模型;对比各个性能参数的参量变化对应的双峰特性曲线,确定敏感度最高的性能参数,形成闭环分析。本发明能够构建核心焊接仿真数据库和提高焊接仿真可靠性,减少测试试验成本。

Description

参量对焊接残余应力及变形影响敏感度获取方法
技术领域
本发明涉及焊接技术领域,特别是涉及一种热物理、热力学参量对焊接残余应力及变形影响敏感度获取方法。
背景技术
通过查阅教科书或参考文献得知,很多外在特征主要受材料单项性能变化的影响:低碳钢焊件纵向焊接残余应力峰值呈现单峰分布特征,而且单峰峰值接近材料屈服极限;铝合金平板中纵向残余应力峰值比其屈服极限低,呈现双峰分布特征;铝合金焊件纵向残余应力峰值约为屈服极限的0.5~0.8倍的主要原因是铝合金材料的热导率高;线膨胀系数是影响铝合金焊接横向收缩的主要原因;焊缝和近缝区金属在力学熔点及以下温度冷却收缩受阻是形成焊接残余应力产生的根本原因。通过对碳钢焊接残余应力场进行测量,发现低碳钢焊接件中纵向焊接残余应力峰值与其屈服极限接近,而铝、钛合金等金属纵向残余应力峰值则低于屈服强度;通过切条释放法和实验室研制的超声波测量系统对薄板铝合金平板焊件中纵向残余应力进行测量,以及参阅其它文献,均证实在一般情况下铝合金焊件纵向残余应力峰值比其屈服极限低。
由于焊接热过程的复杂性,不能通过简单的数学计算和试验测量来获得残余应力和变形的演变过程,更无法获知各材料性能影响的程度。目前工程上焊接变形仿真已应用广泛,但针对残余应力分布等内在关系的模拟仿真因其数学模型复杂、见效慢等问题仍未受到普遍重视。
发明内容
(一)要解决的技术问题
本发明的目的是提供一种参量对焊接残余应力及变形影响敏感度获取方法,以解决针对残余应力分布等内在关系的模拟仿真因其数学模型复杂、见效慢等问题仍未受到普遍重视的问题。
(二)技术方案
为了解决上述技术问题,本发明提供一种参量对焊接残余应力及变形影响敏感度获取方法,其包括如下步骤:
步骤1:利用有限元法建立焊接接头焊接过程的仿真模型;
步骤2:选定母材材料的多个性能参数,依次实验获取各个性能参数的参量;
步骤3:分别以各个性能参数的参量为基准,取该参量的设定差值范围内的多个不同参量数值依次输入仿真模型,获得该性能参数的参量变化对应的双峰特性曲线;
步骤4:对比各个性能参数的参量变化对应的双峰特性曲线,确定敏感度最高的性能参数;
步骤5:重新实验获取敏感度最高的性能参数的参量,以该性能参数的参量为基准,取该参量的设定差值范围内的多个不同参量数值依次输入仿真模型,获得该性能参数的参量变化对应的双峰特性曲线,并替换上一次仿真获取的双峰特性曲线;
步骤6:对比各个性能参数的参量变化对应的双峰特性曲线,确定敏感度最高的性能参数;
若步骤6中确定的敏感度最高的性能参数与上一次确定的敏感度最高的性能参数不相同,则继续执行步骤5和6;
步骤7:将步骤6中确定得到的敏感度最高的性能参数代入步骤1的仿真模型中,形成闭环分析
其中,步骤2中,母材材料的参量包括力学熔点、比热容、热导率、线膨胀量、弹性模量、泊松比和屈服强度中的多项或者全部。
其中,所述双峰特性曲线包括X、Y、Z方向的焊接变形及纵向残余应力分布曲线。
其中,步骤3中,纵向残余应力包括纵向拉应力和纵向压应力。
其中,步骤2中,所述母材材料为铝合金。
其中,步骤5中,实验时通过小孔法或X射线法实测焊接接头的焊接残余应力,通过卷尺或样板实测其焊接变形。
其中,所述仿真模型包括:几何形貌、焊接热输入、夹具约束中的一种或者多种数学物理模型。
(三)有益效果
与现有技术相比,本发明具有以下优点:
本发明提供的一种参量对焊接残余应力及变形影响敏感度获取方法,此方法基于WPS(焊接工艺规范)获取了各参量对焊接残余应力及变形影响敏感度,根据各参量影响特点,进而可仅对影响较大的参量进行重点测试,修正测试大纲后实测材料性能,构建核心焊接仿真数据库和提高焊接仿真可靠性的同时,减少了测试试验成本。
此方法打破了原有通过简单的数学计算和试验测量来获得残余应力和变形的演变过程,由原来的定性分析各参量影响变为定量分析各参量影响,且通过实测和理论验证,形成了有效的分析闭环,可靠性较高。
此方法充分考虑了各影响参量的完整性、可分析性和设置合理性,完善了数学模型的工程应用价值,有效地反映出焊接过程中应力应变的演变历程和规律,更有效的利用数值模拟技术来解决焊接结构应力及焊接变形问题。
此方法通过开源有限元法很容易实现,可在参数敏感性分析、仿真模型评定等场合加以推广。
附图说明
图1为本发明一种参量对焊接残余应力及变形影响敏感度获取方法的获取流程图;
图2为本发明实施例各参量对铝合金焊接残余应力及变形影响敏感度的柱状示意图。
具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明,但不用来限制本发明的范围。
如图1所示,本发明提供了一种热物理、热力学参量对焊接残余应力及变形影响敏感度获取方法,本发明的实施例中,以铝合金作为母材材料,其包括如下步骤:
步骤1:利用有限元法基于WPS(焊接工艺规范)的数据特征建立焊接接头焊接过程的仿真模型,仿真模型具体可以包括几何形貌、焊接热输入、夹具约束中的一种或者多种数学物理模型;
步骤2:选定母材材料的多个性能参数,依次实验获取各个性能参数的参量,理论假设母材材料的其他参量不变,母材材料的一个参量各单项性能由实测材料性能的设定差值范围依次变化,其中,母材材料的参量包括力学熔点、比热容、热导率、线膨胀量、弹性模量、泊松比和屈服强度等,例如,可以先假设母材材料的比热容、热导率、线膨胀量、弹性模量、泊松比和屈服强度的各项性能不变,力学熔点由80%、90%、100%(实测材料性能)、110%、120%依次变化,其他参量以此类推;
步骤3:分别以各个性能参数的参量为基准,取该参量的设定差值范围内的多个不同参量数值依次输入仿真模型,获得该性能参数的参量变化对应的双峰特性曲线,双峰特性曲线包括X(横向)、Y(纵向)、Z(垂向)方向的焊接变形及纵向残余应力分布曲线,其中,纵向残余应力具体包括纵向拉应力和纵向压应力;
步骤4:对比各个性能参数的参量变化对应的双峰特性曲线,确定敏感度最高的性能参数即双峰特性变化曲线变化度最快对应的性能参数,例如,如图2所示,柱形图中依次示出了右侧各参量对铝合金焊接变形、纵向残余应力分布双峰特征的影响程度,铝合金纵向残余拉应力影响较大的性能由大到小依次有:力学熔点(24.29%)、线膨胀系数(10.84%),对铝合金纵向残余压应力影响较大的性能由大到小依次有:泊松比(90.90%),弹性模量(50.54%)、线膨胀系数(47.85%)和力学熔点(38.76%),使铝合金焊接纵向残余应力双峰效应增强的性能由大到小依次有:力学熔点、热导率和泊松比;
步骤5:重新实验获取敏感度最高的性能参数的参量,以该性能参数的参量为基准,取该参量的设定差值范围内的多个不同参量数值依次输入仿真模型,获得该性能参数的参量变化对应的双峰特性曲线,并替换上一次仿真获取的双峰特性曲线,实验时通过小孔法、X射线法等实测焊接结构残余应力,通过卷尺、样板等实测其焊接变形,进行实测验证;
步骤6:对比各个性能参数的参量变化对应的双峰特性曲线,确定敏感度最高的性能参数,修正材料热物理及热力学性能测试大纲,对影响较大的参量重点进行检测,并得到更饱满的数据;
若步骤6中确定的敏感度最高的性能参数与上一次确定的敏感度最高的性能参数不相同,则继续执行步骤5和6;
步骤7:将步骤6中确定得到的敏感度最高的性能参数代入步骤1的仿真模型中,形成闭环分析。
当经过多个闭环分析后,得到稳健的规律性结论,最终得出可靠的影响规律,确定出各参量对焊接残余应力及变形影响敏感度。
由以上实施例可以看出,本发明采用数值模拟技术和实验测量相结合,分析典型焊接接头焊接残余应力和变形,针对焊接残余应力分布内在关联性问题,利用有限元法建立焊接接头的数学物理模型,采用热弹塑性法进行焊接仿真分析,保证接头结构参数和焊接参数一致的情况下,当假设母材材料的其他参量不变,材料性能在实测材料性能前后小范围变化时,分别输出纵向残余应力分布特征,即可总结材料性能的影响规律。
建立精确的材料模型是数值模拟精确建模工作的重要组成部分。本方法预利用铝合金热物理和热力学性能参数的实测和理论假设,通过计算机仿真获知其影响焊接残余应力和变形的共性规律,并实测残余应力及变形、经验公式和理论分析同步验证仿真结果的准确性,经反复的循环分析后,得到稳健的分析结论:热物理性能中的比热容对焊接变形影响较大,热力学性能中的力学熔点、线膨胀系数和泊松比对焊接残余应力影响较大。进而仅对影响较大的参量进行重点测试,构建核心焊接仿真数据库,从而提高焊接仿真可靠性,减少测试试验成本。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (7)

1.一种参量对焊接残余应力及变形影响敏感度获取方法,其特征在于,包括如下步骤:
步骤1:利用有限元法建立焊接接头焊接过程的仿真模型;
步骤2:选定母材材料的多个性能参数,依次实验获取各个性能参数的参量;
步骤3:分别以各个性能参数的参量为基准,取该参量的设定差值范围内的多个不同参量数值依次输入仿真模型,获得该性能参数的参量变化对应的双峰特性曲线;
步骤4:对比各个性能参数的参量变化对应的双峰特性曲线,确定敏感度最高的性能参数;
步骤5:重新实验获取敏感度最高的性能参数的参量,以该性能参数的参量为基准,取该参量的设定差值范围内的多个不同参量数值依次输入仿真模型,获得该性能参数的参量变化对应的双峰特性曲线,并替换上一次仿真获取的双峰特性曲线;
步骤6:对比各个性能参数的参量变化对应的双峰特性曲线,确定敏感度最高的性能参数;
若步骤6中确定的敏感度最高的性能参数与上一次确定的敏感度最高的性能参数不相同,则继续执行步骤5和6;
步骤7:将步骤6中确定得到的敏感度最高的性能参数代入步骤1的仿真模型中,形成闭环分析。
2.根据权利要求1所述的参量对焊接残余应力及变形影响敏感度获取方法,其特征在于,步骤2中,母材材料的参量包括力学熔点、比热容、热导率、线膨胀量、弹性模量、泊松比和屈服强度中的多项或者全部。
3.根据权利要求1所述的参量对焊接残余应力及变形影响敏感度获取方法,其特征在于,所述双峰特性曲线包括X、Y、Z方向的焊接变形及纵向残余应力分布曲线。
4.根据权利要求3所述的参量对焊接残余应力及变形影响敏感度获取方法,其特征在于,步骤3中,纵向残余应力包括纵向拉应力和纵向压应力。
5.根据权利要求1所述的参量对焊接残余应力及变形影响敏感度获取方法,其特征在于,步骤2中,所述母材材料为铝合金。
6.根据权利要求1所述的参量对焊接残余应力及变形影响敏感度获取方法,其特征在于,步骤5中,实验时通过小孔法或X射线法实测焊接接头的焊接残余应力,通过卷尺或样板实测其焊接变形。
7.根据权利要求1所述的参量对焊接残余应力及变形影响敏感度获取方法,其特征在于,所述仿真模型包括:几何形貌、焊接热输入、夹具约束中的一种或者多种数学物理模型。
CN201710305152.6A 2017-05-03 2017-05-03 参量对焊接残余应力及变形影响敏感度获取方法 Active CN107220410B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710305152.6A CN107220410B (zh) 2017-05-03 2017-05-03 参量对焊接残余应力及变形影响敏感度获取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710305152.6A CN107220410B (zh) 2017-05-03 2017-05-03 参量对焊接残余应力及变形影响敏感度获取方法

Publications (2)

Publication Number Publication Date
CN107220410A true CN107220410A (zh) 2017-09-29
CN107220410B CN107220410B (zh) 2020-09-04

Family

ID=59943943

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710305152.6A Active CN107220410B (zh) 2017-05-03 2017-05-03 参量对焊接残余应力及变形影响敏感度获取方法

Country Status (1)

Country Link
CN (1) CN107220410B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110866358A (zh) * 2019-11-12 2020-03-06 桂林电子科技大学 基于实测温度数据的再流焊接工艺仿真模型修正方法
CN111488714A (zh) * 2020-04-10 2020-08-04 桂林电子科技大学 一种热风再流焊喷嘴风速准确计算方法
CN111581862A (zh) * 2020-04-20 2020-08-25 湖南大学 一种焊接接头微区力学性能的等效测试方法
CN111597742A (zh) * 2020-04-17 2020-08-28 扬州大学 一种基于焊接构件联动效应的农机结构残余变形预测方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5992241A (en) * 1995-05-09 1999-11-30 Magyar Allamvasutak Reszvenytarsasag Method and device for determining the neutral temperature of welded tracks
CN103177149A (zh) * 2012-10-29 2013-06-26 中国二十二冶集团有限公司 钢结构焊接变形分析方法
CN103364336A (zh) * 2013-07-15 2013-10-23 江苏科技大学 一种测试环焊缝内部双向残余应力全貌的方法
CN104057210A (zh) * 2013-03-20 2014-09-24 上海通用汽车有限公司 基于强度理论的焊接设计方法
CN105798425A (zh) * 2016-03-28 2016-07-27 南昌航空大学 一种外加磁场装置控制tig焊接残余应力的系统

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5992241A (en) * 1995-05-09 1999-11-30 Magyar Allamvasutak Reszvenytarsasag Method and device for determining the neutral temperature of welded tracks
CN103177149A (zh) * 2012-10-29 2013-06-26 中国二十二冶集团有限公司 钢结构焊接变形分析方法
CN104057210A (zh) * 2013-03-20 2014-09-24 上海通用汽车有限公司 基于强度理论的焊接设计方法
CN103364336A (zh) * 2013-07-15 2013-10-23 江苏科技大学 一种测试环焊缝内部双向残余应力全貌的方法
CN105798425A (zh) * 2016-03-28 2016-07-27 南昌航空大学 一种外加磁场装置控制tig焊接残余应力的系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张艳梅 等: ""一种新的焊接残余应力测量方法"", 《太原理工大学学报》 *
闫忠杰 等: ""TIG熔修对7N01铝合金焊接接头疲劳性能的影响"", 《焊接学报》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110866358A (zh) * 2019-11-12 2020-03-06 桂林电子科技大学 基于实测温度数据的再流焊接工艺仿真模型修正方法
CN110866358B (zh) * 2019-11-12 2023-06-06 桂林电子科技大学 基于实测温度数据的再流焊接工艺仿真模型修正方法
CN111488714A (zh) * 2020-04-10 2020-08-04 桂林电子科技大学 一种热风再流焊喷嘴风速准确计算方法
CN111597742A (zh) * 2020-04-17 2020-08-28 扬州大学 一种基于焊接构件联动效应的农机结构残余变形预测方法
CN111597742B (zh) * 2020-04-17 2023-05-23 扬州大学 一种基于焊接构件联动效应的农机结构残余变形预测方法
CN111581862A (zh) * 2020-04-20 2020-08-25 湖南大学 一种焊接接头微区力学性能的等效测试方法
CN111581862B (zh) * 2020-04-20 2022-04-08 湖南大学 一种焊接接头微区力学性能的等效测试方法

Also Published As

Publication number Publication date
CN107220410B (zh) 2020-09-04

Similar Documents

Publication Publication Date Title
Seleš et al. Numerical simulation of a welding process using a prescribed temperature approach
Ye et al. Influence of groove type on welding-induced residual stress, deformation and width of sensitization region in a SUS304 steel butt welded joint
Jiang et al. A time-variant reliability analysis method based on stochastic process discretization
CN107220410A (zh) 参量对焊接残余应力及变形影响敏感度获取方法
Pu et al. Simulating welding residual stress and deformation in a multi-pass butt-welded joint considering balance between computing time and prediction accuracy
Lee et al. Temperature fields and residual stress distributions in dissimilar steel butt welds between carbon and stainless steels
CN108491652A (zh) 一种基于abaqus的焊接残余应力测量方法
CN109900554B (zh) 一种利用压痕法计算断裂韧性的方法
Barroso et al. Prediction of welding residual stresses and displacements by simplified models. Experimental validation
Xu et al. Temperature and residual stress simulations of the NeT single-bead-on-plate specimen using SYSWELD
Deng et al. Influence of material model on prediction accuracy of welding residual stress in an austenitic stainless steel multi-pass butt-welded joint
Kartal et al. Determination of the profile of the complete residual stress tensor in a VPPA weld using the multi-axial contour method
Carpinteri et al. Fracture mechanics based approach to fatigue analysis of welded joints
CN105975708A (zh) 一种基于数值模拟与数据分析的优化钢管焊接参数方法
Gallo et al. Analysis of creep stresses and strains around sharp and blunt V-notches
Meng et al. An inverse approach to the accurate modelling of 3D-printed sandwich panels with lattice core using beams of variable cross-section
Liang et al. An inverse analysis method to estimate inherent deformations in thin plate welded joints
CN110434443A (zh) 一种电阻点焊仿真方法及系统
Liu et al. Three-dimensional finite element simulation of welding residual stress in RPV with two J-groove welds
CN104008219A (zh) 一种焊接变形分析预测方法
Venkata et al. Numerical simulation of a three-pass TIG welding using finite element method with validation from measurements
Jeyakumar et al. Residual stress evaluation in butt-welded steel plates
Wu et al. Experimental investigation and numerical simulation of welding residual stress in orthotropic steel deck with diaphragm considering solid-state phase transformation
Acevedo et al. Numerical modelling and experimental investigation on welding residual stresses in large‐scale tubular K‐joints
CN104057210B (zh) 基于强度理论的焊接设计方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant