CN111222280B - 一种高合金超高强度钢焊接升温过程中相份数的计算方法 - Google Patents
一种高合金超高强度钢焊接升温过程中相份数的计算方法 Download PDFInfo
- Publication number
- CN111222280B CN111222280B CN202010051351.0A CN202010051351A CN111222280B CN 111222280 B CN111222280 B CN 111222280B CN 202010051351 A CN202010051351 A CN 202010051351A CN 111222280 B CN111222280 B CN 111222280B
- Authority
- CN
- China
- Prior art keywords
- phase
- curve
- welding
- thermal expansion
- temperature
- 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.)
- Active
Links
Landscapes
- Investigating And Analyzing Materials By Characteristic Methods (AREA)
Abstract
本发明公开了一种高合金超高强度钢焊接升温过程中相份数的计算方法,属于高合金超高强度钢焊接技术领域。该方法包括:(1)获得不同加热速率下的热膨胀曲线;(2)利用杠杆法求解奥氏体转变的相份数曲线;(3)基于FORTRAN语言开发组织计算模型,通过优化模型参数,使计算得到的组织转变曲线与实验曲线吻合;(4)基于有限元软件计算焊接升温过程相份数的大小及分布。本发明可以获取高合金超高强度钢在不同焊接工艺参数下焊接升温过程中奥氏体相份数的大小及分布,为准确预测冷却过程马氏体和残余奥氏体的含量提供了技术支持。
Description
技术领域
本发明涉及高合金超高强度钢焊接技术领域,具体涉及一种高合金超高强度钢焊接升温过程中相份数的计算方法。
背景技术
高合金超高强度钢以其优异的塑韧性、良好的抗应力腐蚀性能及焊接性能,被广泛应用在关键承力构件的制造中。高合金超高强度钢的屈服强度σ0.2≥1350MPa,合金元素含量大于10%,其主要应用钢种为马氏体时效钢和低碳、高Co-Ni二次硬化钢,其组织主要由马氏体和残余奥氏体构成。在经过焊接及热处理过程后,使超高强度钢产生马氏体相变和二次硬化作用,以满足产品的强度和设计要求,但在焊接及热处理过程后,马氏体相变产生较大的残余应力往往使结构件变形严重。在实际生产过程中,结构件形状复杂,构件变形对结构件的生产制造有较大影响。采用数值模拟的方式,可以节省大量的人力物力,减少实验周期,为实际生产提供技术指导。
准确地预测结构材料在不同的焊接热循环作用下的组织含量演变规律,对于分析组织转变对应力变形的影响尤为重要,而只有在焊接升温过程中生成的奥氏体才会发生马氏体相变。因此,在高合金超高强度钢的焊接升温过程中,准确预测生成的奥氏体含量的大小及分布具有重要的意义。
针对升温过程组织转变预测,主要手段分析认为AC1、AC3为定值,而实际中升温速率(尤其高速率)对AC1、AC3的变化影响较大。因此,现有高合金超高强度钢焊接升温过程中没有考虑升温速率对AC1、AC3影响,从而影响预测结果的准确性。
发明内容
本发明的目的在于提供一种高合金超高强度钢焊接升温过程中相份数的计算方法,该方法考虑了升温速率对AC1、AC3的影响,为准确地预测冷却过程中马氏体和残余奥氏体的含量提供了技术支持。
为实现上述目的,本发明所采用的技术方案如下:
一种高合金超高强度钢焊接升温过程中相份数的计算方法,包括以下步骤:
(1)获得不同加热速率下的热膨胀曲线;
(2)根据步骤(1)所得热膨胀曲线,利用杠杆法求解所得奥氏体转变的相份数曲线,称为实验曲线Ⅰ;
(3)根据组织计算模型,计算不同升温速率下奥氏体转变的相份数曲线,通过调整组织计算模型中的参数,使计算出的奥氏体转变的相份数曲线与步骤(1)所得实验曲线相吻合时,则进行步骤(4);与实验曲线吻合时的计算出的奥氏体转变的相份数曲线称为计算曲线;
(4)基于有限元软件计算焊接升温过程中相份数的大小及分布;
建立与实际工况相同的焊接有限元模型,设置热源参数及热边界条件,并利用焊接温度场控制方程及步骤(3)中的组织计算模型求解,从而获得焊接温度场和升温过程相份数的大小及分布。
上述步骤(1)的具体过程为:根据实验需要选定不同加热速率,并利用热膨胀仪记录样品在不同加热速率下的热膨胀曲线(横坐标为温度,纵坐标为样品尺寸变化(即膨胀量))。由于相变引起的体积改变破坏了膨胀量与温度之间的线性关系,所以利用切线法标定不同升温速率下的AC1、AC3点,绘制CHT曲线。
上述步骤(2)中,求解奥氏体转变的相份数曲线的过程为:在热膨胀曲线中,发生相变前的直线段a的末端点(相变点)为A,绘制该直线段a的延长线,相变完全之后的直线段e的前端点(相变点)为E,绘制该直线段e的延长线;由于相变的发生,热膨胀曲线从A到E的中间段的斜率是变化的,温度T时的膨胀量为C,过C点作横坐标的垂线,与直线段a的延长线和直线段e的延长线分别相交于B点和D点;按照杠杆法,根据公式(1)计算在温度T时新相的份数W;并绘制以温度为横坐标、新相份数W为纵坐标的奥氏体转变的相份数曲线;新相是指升温过程中新产生的奥氏体相;
上述步骤(3)中,所述组织计算模型如公式(2)所示;
公式(2)中:P为相份数的计算值,t为时间,为加热速率,T为温度,Peq为平衡时的相份数,τ为延迟时间。
上述步骤(4)中,所述焊接温度场控制方程如公式(3)所示;
公式(3)中,T为温度,t为时间,λ为导热系数,ρ为材料的密度,c为材料的比热容,Q为外部热源输入的能量。
本发明的有益效果:
本发明提供一种高合金超高强度钢焊接升温过程中相份数的计算方法,针对焊接升温过程组织计算,通过不断调整组织计算模型中的参数,使计算得到的组织转变曲线与实验曲线吻合。得到的相份数计算结果与实验得到的结果吻合良好,表明本发明适合高合金超高强度钢焊接升温过程中相份数的计算。
附图说明
图1为杠杆法计算组织含量示意图。
图2为AF1410钢CHT图与模拟结果对比。
图3为奥氏体相的转变曲线与模拟结果对比。
图4为计算得到焊接完成时奥氏体分布云图。
具体实施方式
下面以AF1410钢焊接过程为例对焊接升温过程相份数进行计算,验证本发明的有益效果。
步骤一:获得不同加热速率下的热膨胀曲线;
将AF1410钢加工若干Φ3×10mm的试样,使用L78 RITA快速加热热膨胀仪进行实验。针对焊接升温过程,由于升温速率较高,设定20℃/s、50℃/s、100℃/s、200℃/s的温度梯度。由于相变引起的体积改变破坏了膨胀量与温度之间的线性关系,所以利用切线法标定不同升温速率下的AC1、AC3点,绘制CHT曲线。根据AC1、AC3随温度转变的趋势,线性外推至1000℃/s;
步骤二:利用杠杆法求解奥氏体转变的相份数曲线;
图1为杠杆法计算组织含量示意图。在加热过程中,如果不发生相变,热膨胀曲线按照AB延长。相变完全完成后,热膨胀曲线按照DE延长。由于相变的发生,曲线按照ACE段变化,所以BC段即为相变带来的影响。按照杠杆法,某一温度下新相的份数W为:
求解各温度下的新相转变份数,即可得到奥氏体转变的相份数曲线;
步骤三:基于FORTRAN语言开发组织计算模型公式(2),通过优化模型参数,使计算得到的组织转变曲线与步骤二中的实验曲线吻合;
其中,P为相份数的计算值,t为时间,为加热速率,T为温度,Peq为平衡时的相份数,τ为延迟时间。
模拟后的AF1410钢CHT图及奥氏体相的转变曲线与实验结果对比分别如图2和图3所示。
步骤四:基于有限元软件计算焊接升温过程相份数及分布;
试板尺寸为140mm×120mm×6mm,采用钨极惰性气体保护焊自熔焊焊接方法,焊接工艺参数如表1所示。建立与实际工况相同的焊接有限元模型,划分网格,共23808个单元,29059个节点。设置热源参数及热边界条件,采用ABAQUS有限元软件计算求解温度场控制方程(公式3)及开发的组织计算模型(公式2),获得焊接温度场和升温过程相份数的大小及分布。计算得到焊接完成时奥氏体分布云图如图4所示。
焊接温度场控制方程为:
其中,T为温度,t为时间,λ为导热系数,ρ为材料的密度,c为材料的比热容,4为外部热源输入的能量。
表1AF1410钢自熔焊焊接工艺参数
为了验证参数拟合结果的准确性,将计算结果与实验进行比较。通过比较可以看出,本发明所提供的相份数计算方法可以准确地反映不同升温速率下奥氏体相的转变含量与趋势。
Claims (2)
1.一种高合金超高强度钢焊接升温过程中相份数的计算方法,其特征在于:该方法包括以下步骤:
(1)获得不同加热速率下的热膨胀曲线,具体过程为:根据实验需要选定不同加热速率,并利用热膨胀仪记录样品在某一加热速率下的热膨胀曲线;绘制热膨胀曲线时,由于相变引起的体积改变破坏了膨胀量与温度之间的线性关系,所以利用切线法标定不同升温速率下的AC1、AC3点,绘制CHT曲线;
(2)根据步骤(1)所得热膨胀曲线,利用杠杆法求解奥氏体转变的相份数曲线,称为实验曲线Ⅰ;求解奥氏体转变的相份数曲线的过程为:在热膨胀曲线中,发生相变前的直线段a的末端点为A,绘制该直线段a的延长线,相变完全之后的直线段e的前端点为E,绘制该直线段e的延长线;由于相变的发生,热膨胀曲线从A到E的中间段的斜率是变化的,温度T时的膨胀量为C,过C点作横坐标的垂线,与直线段a的延长线和直线段e的延长线分别相交于B点和D点;按照杠杆法,根据公式(1)计算在温度T时新相的份数W;并绘制以温度为横坐标、新相份数W为纵坐标的奥氏体转变的相份数曲线;新相是指升温过程中新产生的奥氏体相;
(3)根据组织计算模型,计算不同升温速率下奥氏体转变的相份数曲线,通过调整组织计算模型中的参数,使计算出的奥氏体转变的相份数曲线与步骤(1)所得实验曲线相吻合时,则进行步骤(4);与实验曲线吻合时的计算出的奥氏体转变的相份数曲线称为计算曲线;所述组织计算模型如公式(2)所示;
公式(2)中,P为相份数的计算值,t为时间,为加热速率,T为温度,Peq为平衡时的相份数,τ为延迟时间;
(4)基于有限元软件计算焊接升温过程中相份数的大小及分布;
建立与实际工况相同的焊接有限元模型,设置热源参数及热边界条件,并利用焊接温度场控制方程及步骤(3)中的组织计算模型求解,从而获得焊接温度场和升温过程相份数的大小及分布。
2.根据权利要求1所述的高合金超高强度钢焊接升温过程中相份数的计算方法,其特征在于:步骤(4)中,所述焊接温度场控制方程如公式(3)所示;
公式(3)中,T为温度,t为时间,λ为导热系数,ρ为材料的密度,c为材料的比热容,Q为外部热源输入的能量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010051351.0A CN111222280B (zh) | 2020-01-17 | 2020-01-17 | 一种高合金超高强度钢焊接升温过程中相份数的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010051351.0A CN111222280B (zh) | 2020-01-17 | 2020-01-17 | 一种高合金超高强度钢焊接升温过程中相份数的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111222280A CN111222280A (zh) | 2020-06-02 |
CN111222280B true CN111222280B (zh) | 2023-10-20 |
Family
ID=70831156
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010051351.0A Active CN111222280B (zh) | 2020-01-17 | 2020-01-17 | 一种高合金超高强度钢焊接升温过程中相份数的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111222280B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113702420B (zh) * | 2021-08-19 | 2024-01-09 | 鞍钢股份有限公司 | 一种焊接热循环冷却过程中相变温度的测定方法 |
CN114664393A (zh) * | 2022-01-05 | 2022-06-24 | 北京理工大学 | 基于热膨胀标定高瞬时性固态相变jmak方程参数的方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101871899A (zh) * | 2010-06-22 | 2010-10-27 | 天津钢铁集团有限公司 | 钢铁材料控冷cct曲线测定方法 |
CN103993138A (zh) * | 2014-05-29 | 2014-08-20 | 东北大学 | 预测高强塑积钢淬火中马氏体组织演变的方法 |
CN109858085A (zh) * | 2018-12-26 | 2019-06-07 | 钢铁研究总院 | 一种金属材料热处理过程中的奥氏体化测定方法 |
WO2019173681A1 (en) * | 2018-03-08 | 2019-09-12 | Northwestern University | Carbide-free bainite and retained austenite steels, producing method and applications of same |
-
2020
- 2020-01-17 CN CN202010051351.0A patent/CN111222280B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101871899A (zh) * | 2010-06-22 | 2010-10-27 | 天津钢铁集团有限公司 | 钢铁材料控冷cct曲线测定方法 |
CN103993138A (zh) * | 2014-05-29 | 2014-08-20 | 东北大学 | 预测高强塑积钢淬火中马氏体组织演变的方法 |
WO2019173681A1 (en) * | 2018-03-08 | 2019-09-12 | Northwestern University | Carbide-free bainite and retained austenite steels, producing method and applications of same |
CN109858085A (zh) * | 2018-12-26 | 2019-06-07 | 钢铁研究总院 | 一种金属材料热处理过程中的奥氏体化测定方法 |
Non-Patent Citations (1)
Title |
---|
迟露鑫 ; .0Cr18Ni9不锈钢焊接温度场模拟与验证.焊接技术.2015,(11),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN111222280A (zh) | 2020-06-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109858085B (zh) | 一种金属材料热处理过程中的奥氏体化测定方法 | |
Yang et al. | FEM simulation of quenching process in A357 aluminum alloy cylindrical bars and reduction of quench residual stress through cold stretching process | |
CN111222280B (zh) | 一种高合金超高强度钢焊接升温过程中相份数的计算方法 | |
CN103886125B (zh) | 一种钛合金热复合成形数值模拟方法 | |
CN105373683A (zh) | 一种20CrMnTiH钢热变形过程微观组织演化规律的预测方法 | |
Li et al. | Experimental study and numerical simulation of dynamic recrystallization behavior of a micro-alloyed plastic mold steel | |
CN105586476A (zh) | 基于数值模拟技术的高强钢淬火工艺优化方法 | |
US20110098989A1 (en) | Systems and methods for predicting heat transfer coefficients during quenching | |
Cai et al. | Constitutive equation and model validation for 33Cr23Ni8Mn3N heat-resistant steel during hot compression | |
CN102339344B (zh) | 一种动态再结晶模型参数反分析识别方法 | |
CN112926213B (zh) | 一种热损伤边界元测定方法、***、介质、设备、终端 | |
CN103745114A (zh) | 一种钛合金应力松弛数值及回弹计算方法 | |
Yang et al. | Flow behavior and processing maps of a low-carbon steel during hot deformation | |
CN107220410B (zh) | 参量对焊接残余应力及变形影响敏感度获取方法 | |
Xue et al. | Mechanical behaviour and strain rate sensitivity analysis of TA2 by the small punch test | |
CN115795933A (zh) | 一种金属嵌件注塑成型优化方法、装置及存储介质 | |
Taddesse et al. | Combined notch and size effect modeling of metallic materials for LCF using plasticity reformulated critical distance theory | |
Chen et al. | Multiaxial thermo-mechanical fatigue life prediction based on notch local stress-strain estimation considering temperature change | |
Liu et al. | Thermal fatigue life prediction method for die casting mold steel based on the cooling cycle | |
CN102750453B (zh) | 一种用于轧机设计的材料变形抗力统一模型的构造方法 | |
CN108536968A (zh) | 一种焊接过程中的应力与变形的调控方法 | |
CN116756870A (zh) | 基于数值模拟仿真的电弧堆焊工艺参数优化方法 | |
Peng et al. | Dynamic recrystallization behavior under inhomogeneous thermomechanical deformation state | |
Wen et al. | Development and application of cast steel numerical simulation system for heat treatment | |
JP2017138101A (ja) | 硬さ解析方法および硬さ解析プログラム |
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 |