CN102519801B - 用真应力-真应变曲线建立铝合金动态再结晶模型的方法 - Google Patents
用真应力-真应变曲线建立铝合金动态再结晶模型的方法 Download PDFInfo
- Publication number
- CN102519801B CN102519801B CN201110425352.8A CN201110425352A CN102519801B CN 102519801 B CN102519801 B CN 102519801B CN 201110425352 A CN201110425352 A CN 201110425352A CN 102519801 B CN102519801 B CN 102519801B
- Authority
- CN
- China
- Prior art keywords
- strain
- stress
- true
- peak
- dynamic recrystallization
- 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.)
- Expired - Fee Related
Links
- 238000001953 recrystallisation Methods 0.000 title claims abstract description 29
- 238000000034 method Methods 0.000 title claims abstract description 22
- 229910000838 Al alloy Inorganic materials 0.000 title claims abstract description 15
- 238000013178 mathematical model Methods 0.000 claims abstract description 19
- 238000012360 testing method Methods 0.000 claims abstract description 12
- 239000000463 material Substances 0.000 claims abstract description 11
- 238000009864 tensile test Methods 0.000 claims description 3
- 230000004913 activation Effects 0.000 claims description 2
- 230000003111 delayed effect Effects 0.000 abstract description 3
- 238000011160 research Methods 0.000 abstract description 3
- 238000012113 quantitative test Methods 0.000 abstract description 2
- 238000012417 linear regression Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 229910052751 metal Inorganic materials 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000007493 shaping process Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 239000011149 active material Substances 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公布了一种用真应力-真应变曲线建立铝合金动态再结晶模型的方法。属于材料微观组织分析领域。采用高温拉伸试验得到的真应力-真应变曲线求出材料的硬化率,通过分析硬化率的变化确定不同条件下的峰值应变、峰值应力和稳态应力,再根据经验公式建立峰值应变、临界应变数学模型和动态再结晶体积分数。本发明克服了用金相方法研究材料的再结晶行为的滞后、费时、不全面和不精确的缺点,结合有限元软件后可以定量分析铝合金的动态再结晶体积分数变化情况。
Description
技术领域
本发明属于材料微观组织分析领域,主要涉及一种用真应力-真应变曲线建立铝合金动态再结晶数学模型的方法。
背景技术
利用铝合金在超塑性状态下的优异变形性能发展起来的超塑性成形技术,可以在小吨位设备上成形形状复杂、其他塑性加工工艺难以或不能成形的零件,能够满足轨道车辆铝合金零件的成形需要。铝合金具有较高的强度,良好的塑性、抗腐蚀性及加工性等优点,已成为汽车、轨道交通工业超塑成形热门材料,可大大减轻车辆自重,促进交通产业节能和环保发展。控制成形性能最有效的方法就是控制晶粒大小。铝合金超塑变形中经常伴随晶粒长大,晶粒长大将增加流动应力,降低材料的超塑性能,严重的会造成成形失败。研究晶粒长大及其对零件成形的影响是铝合金超塑成形微观组织模拟的主要目的。晶粒的长大与热成形过程中或热处理过程中的动态再结晶,亚动态再结晶和静态再结晶密切相关。再结晶受一些过程变量,如成形温度,应变速率,原始晶粒大小和变形量的影响。研究材料的再结晶行为常采用金相观察方法,这种方法滞后、费时又不全面和精确。近年来发展了用有限元方法模拟材料在变形过程中的微观组织演变,建立微观组织演变的数学模型是关键。用金相方法来得到数学模型滞后又费时。该发明用真应力-真应变曲线建立铝合金动态再结晶数学模型,弥补了金相方法的不足,可以定量分析铝合金的动态再结晶体积分数变化情况。
发明内容
针对现有用金相方法确定动态再结晶数学模型的方法,本发明提出一种用真应力-真应变曲线建立铝合金动态再结晶数学模型的方法。
本发明基于高温拉伸试验和线性回归方法,采用下述步骤建立铝合金动态再结晶数学模型。
1)在微机控制电子万能热拉伸试验机RG2000-20上进行高温拉伸试验。应变速率范围选择为6.56×10-5~6.56×10-3S-1。试验温度定为623K、673K、723K、773K、793K、808K。
2)用高温拉伸试验得到的真应力-真应变曲线求解硬化率。当应变速率与变形温度一定时,应力随应变的变化率称为硬化率,即θ=dσ/dε。
3)通过不同变形条件下θ-ε和θ-σ曲线求出不同变形条件下的峰值应变εp、峰值应力σp和稳态应力σss。
4)用3)求得的峰值应变、峰值应力建立峰值应变及临界应变数学模型,可求出数学模型中各系数值。
临界应变数学模型:εc=0.8εp (1)
峰值应变数学模型:
其中a1、n1、m1为材料常数,d0为原始晶粒大小,为应变速率,单位为S-1,Q1为再结晶激活能,单位为J/mol,R为气体常数,T为试验温度;
5)用3)峰值应力和稳态应力确定动态再结晶体积分数。
动态再结晶体积分数可表示为Xdrex=(σp-σ)/(σp-σss) (3)
根据Avrami方程,动态再结晶体积分数还可以表示为:
Xdrex(ε≤εc)=0 (4)
将上式进行转换,得式ln[-ln(1-Xdrex)]=lnβd+kd lnt,其中
画出不同变形条件下的ln[-ln(1-Xdrex)]和lnt的线性关系。从曲线的斜率和截距可得常数βd和kd。从而确定数学模型中各系数。
附图说明
图1为高温拉伸试样几何尺寸,板材厚度为2mm。
图2中a,b分别为在623K和673K条件下,不同应变速率时AA5083硬化率θ与应力之间的关系。
图3中a,b分别为在623K和673K条件下,不同应变速率时AA5083硬化率θ与应变之间的关系。
图4为温度623K下不同应变速率下ln[-ln(1-Xdrex)]-lnt。
图5为应变速率3.28E-3下不同温度下ln[-ln(1-Xdrex)]-lnt。
具体实施方式
实施例
试验采用5083板材,其化学成分(质量分数)见表1
表1 试验5083板材化学成分
Si | Fe | Cu | Mn | Mg | Cr | Ni | Zn | Ti | Zr | Be | Na | V | Al |
0.08 | 0.05 | 0.01 | 0.6 | 4.61 | 0.1 | 0.15 | 0.034 | 0.029 | 0.24 | 0.0006 | 0.0008 | 0.048 | 其余 |
高温单拉试验在微机控制电子万能热拉伸试验机RG2000-20上进行,试样的几何尺寸如图1所示,板材厚度为2mm。
采用的应变速率为6.56×10-5~6.56×10-3S-1;温度为623K、673K、723K、773K、793K、808K。拉伸变形前试样需保温5分钟,使试样均匀达到设定目标温度。
用高温拉伸试验得到的真应力-真应变曲线求解硬化率。见附图2和3。通过不同变形条件下θ-ε和θ-σ曲线求出不同变形条件下的峰值应变、峰值应力和稳态应力。
用式(1)和式(2)运用不同变形条件下AA5083的εp,用线性回归确定公式(2)中常数,这里令n1=0,即忽略d0的影响。
将公式(2)两边取对数,转化为线性回归问题
分别画出定温度下和定应变速率下用origin7.0线性拟合求平均值后得a1=0.009728;m1=0.06539;Q1=5114.849J/mol。
可得出AA5083临界应变方程:εc=0.8εp
画出不同变形条件下的ln[-ln(1-Xdrex)]和lnt的线性关系。如附图4,5所示。从曲线的斜率和截距可得常数βd和kd。见表2和表3.
表2 温度623K下不同应变速率的动态再结晶动力学参数
表3 应变速率3.28E-3下不同温度的动态再结晶动力学参数
因此就建立了AA5083的动态再结晶数学模型。
Claims (1)
1.一种用真应力-真应变曲线建立铝合金动态再结晶模型的方法,其特征在于所述方法包括如下步骤:
a、在微机控制电子万能热拉伸试验机RG2000-20上进行高温拉伸试验:应变速率范围选择为6.56×10-5~6.56×10-3S-1;试验温度定为623K、673K、723K、773K、793K、808K;
b、用高温拉伸试验得到的真应力-真应变曲线求解硬化率:当应变速率与变形温度一定时,应力随应变的变化率称为硬化率,即θ=dσ/dε;
c、通过不同变形条件下θ-ε和θ-σ曲线求出不同变形条件下的峰值应变εp、峰值应力σp和稳态应力σss;
d、用c中求得的峰值应变、峰值应力建立峰值应变及临界应变数学模型,可求出数学模型中各系数值;
临界应变数学模型:εc=0.8εp;
峰值应变数学模型:
其中a1、n1、m1为材料常数,d0为原始晶粒大小,为应变速率,单位为S-1,Q1为再结晶激活能,单位为J/mol,R为气体常数,T为试验温度;
e、用c中峰值应力和稳态应力确定动态再结晶体积分数:
动态再结晶体积分数表示为Xdrex=(σp-σ)/(σp-σss);
根据Avrami方程,动态再结晶体积分数还可以表示为:
Xdrex(ε≤εc)=0;
将上式进行转换,得式ln[-ln(1-Xdrex)]=lnβd+kdlnt,其中
画出不同变形条件下的ln[-ln(1-Xdrex)]和lnt的线性关系,从曲线的斜率和截距得常数βd和kd,从而确定数学模型中各系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110425352.8A CN102519801B (zh) | 2011-12-16 | 2011-12-16 | 用真应力-真应变曲线建立铝合金动态再结晶模型的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110425352.8A CN102519801B (zh) | 2011-12-16 | 2011-12-16 | 用真应力-真应变曲线建立铝合金动态再结晶模型的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102519801A CN102519801A (zh) | 2012-06-27 |
CN102519801B true CN102519801B (zh) | 2015-07-29 |
Family
ID=46290797
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110425352.8A Expired - Fee Related CN102519801B (zh) | 2011-12-16 | 2011-12-16 | 用真应力-真应变曲线建立铝合金动态再结晶模型的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102519801B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109255175A (zh) * | 2018-08-31 | 2019-01-22 | 西北工业大学 | 一种基于临界再结晶残余应力的再结晶控制方法 |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102912200A (zh) * | 2012-11-13 | 2013-02-06 | 东北轻合金有限责任公司 | 一种超大规格铝合金扁铸锭的制造方法 |
CN103471932B (zh) * | 2013-09-26 | 2015-09-23 | 北京机电研究所 | 金属材料的应力-应变曲线测量及应用方法 |
CN105740513B (zh) * | 2016-01-22 | 2019-01-08 | 东北大学 | 一种gh4169合金热变形动态再结晶模拟方法 |
CN106053754B (zh) * | 2016-07-06 | 2017-11-14 | 中南大学 | 一种预测时变工况下高合金化材料动态再结晶分数的方法 |
CN107058923B (zh) * | 2017-01-16 | 2018-12-07 | 青岛理工大学 | 一种镁合金热加工工艺优化新方法 |
CN107326314B (zh) * | 2017-07-05 | 2019-04-05 | 中南大学 | 一种预测含铌镍基合金中δ相动态溶解体积分数的方法 |
CN108169019B (zh) * | 2017-12-26 | 2020-06-16 | 中国兵器工业第五九研究所 | 一种准静态塑性压缩应力应变参量的识别方法 |
CN110068507B (zh) * | 2018-01-22 | 2021-07-23 | 中国科学院金属研究所 | 一种对传统再结晶模型进行修正的方法 |
CN108489820B (zh) * | 2018-02-13 | 2019-06-28 | 鞍钢股份有限公司 | 一种高温拉伸试验真应力的获取方法 |
CN110245382B (zh) * | 2019-05-10 | 2023-08-22 | 本钢板材股份有限公司 | 一种确定金属动态再结晶体积分数的Avrami数学模型系数的方法 |
CN110296892B (zh) * | 2019-08-08 | 2020-09-11 | 中国矿业大学(北京) | 基于能量分析的岩石损伤演化过程中特征应力确定方法 |
CN110779797B (zh) * | 2019-11-08 | 2021-10-01 | 南京航空航天大学 | 一种测量金属拉伸试验过程中塑性应变比的方法 |
CN110702711B (zh) * | 2019-11-21 | 2022-02-25 | 山东建筑大学 | 非晶态合金结晶动力学特征的原位exafs表征方法 |
CN111380899A (zh) * | 2019-11-29 | 2020-07-07 | 中国科学院金属研究所 | 一种通过轧制模拟过程温升修正锆合金流变应力的方法 |
CN113702613B (zh) * | 2021-08-19 | 2024-02-13 | 鞍钢股份有限公司 | 一种确定材料发生动态再结晶临界条件的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1641356A (zh) * | 2004-12-13 | 2005-07-20 | 武汉科技大学 | 硬线产品显微组织与力学性能预报*** |
CN101046682A (zh) * | 2007-04-28 | 2007-10-03 | 武汉科技大学 | 一种预测热轧含Nb带钢组织及力学性能的方法 |
CN101319977A (zh) * | 2008-07-18 | 2008-12-10 | 哈尔滨工业大学 | 真应力-真应变计算模型及测试*** |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8515688B2 (en) * | 2009-03-12 | 2013-08-20 | GM Global Technology Operations LLC | Systems and methods to predict fatigue lives of aluminum alloys under multiaxial loading |
-
2011
- 2011-12-16 CN CN201110425352.8A patent/CN102519801B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1641356A (zh) * | 2004-12-13 | 2005-07-20 | 武汉科技大学 | 硬线产品显微组织与力学性能预报*** |
CN101046682A (zh) * | 2007-04-28 | 2007-10-03 | 武汉科技大学 | 一种预测热轧含Nb带钢组织及力学性能的方法 |
CN101319977A (zh) * | 2008-07-18 | 2008-12-10 | 哈尔滨工业大学 | 真应力-真应变计算模型及测试*** |
Non-Patent Citations (3)
Title |
---|
Prediction of hot deformation behaviour of Fe-25Mn-3Si-3Al TWIP steel;Li, Dejun等;《MATERIALS SCIENCE AND ENGINEERING A-STRUCTURAL MATERIALS PROPERTIES MICROSTRUCTURE AND PROCESSING》;20111025;第528卷(第28期);第8084-8089页 * |
基于温度变化的AA5083合金拉伸变形行为描述;徐雪峰,等;《塑性工程学报》;20110228;第18卷(第1期);第96-100页 * |
铝合金6061的热变形力学行为与微观组织演化规律;金泉林;《材料热处理学报》;20110630;第32卷(第6期);第51-57页 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109255175A (zh) * | 2018-08-31 | 2019-01-22 | 西北工业大学 | 一种基于临界再结晶残余应力的再结晶控制方法 |
CN109255175B (zh) * | 2018-08-31 | 2021-04-02 | 西北工业大学 | 一种基于临界再结晶残余应力的再结晶控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102519801A (zh) | 2012-06-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102519801B (zh) | 用真应力-真应变曲线建立铝合金动态再结晶模型的方法 | |
Song et al. | Improving tensile and compressive properties of magnesium alloy plates by pre-cold rolling | |
Fan et al. | Experimental investigation on hot forming–quenching integrated process of 6A02 aluminum alloy sheet | |
Singh et al. | Prediction of mechanical properties of extra deep drawn steel in blue brittle region using Artificial Neural Network | |
Wang et al. | Study of the dynamic recrystallization of Ti–6.5 Al–3.5 Mo–1.5 Zr–0.3 Si alloy in β-forging process via Finite Element Method modeling and microstructure characterization | |
CN107254646B (zh) | 提高6000系铝合金自然时效稳定性的热处理方法及铝合金板材 | |
CN110245382B (zh) | 一种确定金属动态再结晶体积分数的Avrami数学模型系数的方法 | |
CN109855958B (zh) | 一种金属材料拉伸性能的预测方法 | |
Ji et al. | DRX behavior and microstructure evolution of 33Cr23Ni8Mn3N: Experiment and finite element simulation | |
Guo et al. | Microstructure evolution and mechanical responses of Al–Zn–Mg–Cu alloys during hot deformation process | |
CN102224267A (zh) | 铝-铜-锂合金制成的产品 | |
Quan et al. | Correspondence between grain refinements and flow softening behaviors at Nimonic 80A superalloy under different strain rates, temperatures and strains | |
CN107976462A (zh) | 一种优化铝合金热加工工艺的方法 | |
CN107058923B (zh) | 一种镁合金热加工工艺优化新方法 | |
Stachowicz et al. | Warm forming of stainless steel sheet | |
CN108710729A (zh) | 一种基于组织与性能控制的铝合金热成形工艺制定方法 | |
CN101748255B (zh) | 一种提高430铁素体不锈带钢成形性能的方法 | |
Jiang et al. | Hot deformation behavior of Al-Zn-Mg-Cu alloy during compression tests | |
CN103834885A (zh) | 一种提高铝合金板材塑性的热处理方法 | |
Sajadifar et al. | Dynamic recrystallization behavior and hot deformation characteristics in 4340 steel | |
Cai et al. | Study of springback for high strength aluminium alloys under hot stamping | |
CN106086658B (zh) | 一种高强度中碳含铬低锰轻质钢及制备方法 | |
Guo et al. | Hot compressive behavior of Ti–3.0 Al–3.7 Cr–2.0 Fe–0.1 B titanium alloy | |
CN103388116B (zh) | 一种高效轧制Mg-Al-Zn系镁合金的方法 | |
CN114021339A (zh) | 数据驱动的金属薄壁零件性能优化方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150729 Termination date: 20181216 |