CN109870257B - 一种板材厚度方向淬火残余应力分布预测方法 - Google Patents
一种板材厚度方向淬火残余应力分布预测方法 Download PDFInfo
- Publication number
- CN109870257B CN109870257B CN201711265746.5A CN201711265746A CN109870257B CN 109870257 B CN109870257 B CN 109870257B CN 201711265746 A CN201711265746 A CN 201711265746A CN 109870257 B CN109870257 B CN 109870257B
- Authority
- CN
- China
- Prior art keywords
- residual stress
- distribution
- thickness direction
- groups
- dispersion
- 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
Images
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
一种板材厚度方向淬火残余应力分布预测方法,包括如下步骤:(1)对n组材料试样厚度方向残余应力分布进行数值模拟或实际测试分析,通过离差标准化处理得到残余应力离差标准化分布结果;(2)对分布结果进行平均,获得残余应力离差标准化统计分布,拟合得到相对于标准化厚度z的残余应力分布函数J(z);(3)统计分析n组材料表面残余应力绝对值与内部残余应力绝对值,拟合得到内部残余应力与表面残余应力的对应关系;(4)对表面残余应力进行测试,得到表面残余应力数值,并计算内部残余应力;(5)利用表面残余应力、内部残余应力数值以及分布函数J(z)预测其他规格样品厚度方向残余应力分布。本发明能够实现材料淬火残余应力的快速预测评价。
Description
技术领域
本发明涉及一种板材厚度方向淬火残余应力分布预测方法,属于材料分析检测领域。
背景技术
钢铁、铝合金、镁合金等金属板材由于其各自的性能优势,在工业领域广泛应用于结构件制备。而目前结构件加工倾向于向大型化整体化发展,以减少零件数量,提高装配精度与装配效率。金属板材在加工过程中由于热处理、机加工等因素,不可避免地会引入残余应力,而材料中存在过高的残余应力可能会导致板材在加工成为零件的过程中发生翘曲,严重的甚至发生开裂,严重影响了制造效率与成本。通过残余应力的表征与预测,获得材料残余应力分布状态,可为材料制备加工工艺优化提供依据,也作为材料出厂状态评价,避免加工变形导致材料报废。但如何实现对板材残余应力准确评价与预测仍然是行业内的共性难题,因此如何快速、准确地获得板材残余应力分布是行业内各企业的迫切需求。
现有残余应力测试方法中,钻孔法、X射线衍射法等方法可应用于工厂规格材料应力测试,但两种方法都只能测试表面残余应力,无法对材料厚度方向残余应力分布进行表征;而可进行内部残余应力测试的裂纹柔度法、应变云图法、中子衍射法、同步辐射X射线衍射法均无法应用于较大规格的样品测试,且裂纹柔度法、应变云图法对材料破坏性大,中子衍射法、同步辐射X射线衍射法只能在特定实验室条件下进行测试。
发明内容
针对上述问题,本发明的目的在于提供一种板材厚度方向淬火残余应力分布预测方法,高效,快速地得到板材厚度方向淬火残余应力分布。
为实现上述目的,本发明采用以下技术方案:
一种板材厚度方向淬火残余应力分布预测方法,包括如下步骤:
(1)对选定材料的n组不同厚度规格、相同状态的试样厚度方向残余应力分布进行数值模拟或实际测试分析,得到n组试样应力分布数据,对n组不同厚度及n组试样应力分布数据分别进行离差标准化处理,离差标准化的转换函数如下式:
式中,max为n组不同厚度数据或n组试样应力分布数据最大值,min为n组不同厚度数据或n组试样应力分布数据最小值,x’代表离差标准化处理后的数据,x代表离差标准化处理前的数据;
(2)对n组离差标准化处理结果进行平均,得到离差标准化统计分布,对离差标准化统计分布进行多项式拟合,得到相对于标准化厚度z的残余应力分布函数J(z);
(3)统计分析n组试样表面残余应力绝对值与厚度方向中心位置残余应力绝对值,拟合得到厚度方向中心位置残余应力绝对值与表面残余应力绝对值的对应关系;
(4)对表面残余应力进行测试,测试方法选择钻孔法或X射线衍射法,得到表面残余应力绝对值,并依据步骤(3)得到的厚度方向中心位置残余应力绝对值与表面残余应力绝对值对应关系,计算厚度方向中心位置残余应力;
(5)利用表面残余应力绝对值、厚度方向中心位置残余应力绝对值以及分布函数J(z)预测其他规格样品厚度方向残余应力分布。
所述步骤(1)中残余应力分布实际测试分析采用裂纹柔度法或中子衍射法进行,对不同厚度的试样测试位置选择在离差标准化厚度相同的位置进行测试,厚度方向实测点数量不少于9个;所述步骤(1)中n应大于等于7个。
本发明的优点在于:
本发明在实验测试确定相同淬火条件下小规格试样淬火残余应力分布和大规格厚板表面中心点应力水平基础上,实现材料淬火残余应力的快速预测评价。本发明通过计算厚度方向残余应力分布函数J(z),结合表面点残余应力测试,高效,快速地得到板材厚度方向淬火残余应力分布,可为产品残余应力状态评价提供参考。
附图说明
图1是本发明的板材厚度方向淬火残余应力分布预测方法的流程图。
图2是板材淬火残余应力标准化分布函数及其拟合图形。
具体实施方式
图1是本发明的板材厚度方向淬火残余应力分布预测方法的流程图。以下参照图1对本发明的具体实施方式进行说明。
以不同规格7055铝合金板材为例,选取2000mm×1000mm×30/50/70/90/110/130/150/170/200mm共9个模型的淬火残余应力分布模拟结果,对厚度和厚度方向应力分布分别进行离差标准化处理。令表面最大压应力为最小值、心部最大拉应力为最大值,处理公式如公式(1)。
由此得到的残余应力离差标准化分布结果如表1所示。
表1不同厚度下淬火残余应力0-1标准化分布
随后对9组不同的应力数据进行平均,绘制淬火残余应力离差标准化统计分布曲线,如图2所示,并对其进行多项式拟合,得到相对于0-1标准化厚度z的淬火残余应力分布函数J(z):
J(z)=3.327e-8z9-65.599z8+262.39z7-556.75z6+751.87z5-606.45z4
+265.91z3-59.113z2+7.7352z-0.0059(0≤z≤1) (4)
对铝合金板材而言,长厚比和宽厚比一般较大,影响残余应力水平的主要是厚度,厚板淬火心部最大拉应力σ淬火内部可表示为:
其中C、A1、A2为待定系数,H为板材厚度。对上式两边取自然对数。
lnσ淬火内部=ln C+A1·ln H+A2·ln|σ淬火表面|(6)
从各厚度的厚向模拟应力分布中确定各厚度的表面最大压应力和心部最大拉应力水平,如表2所示。
表2不同厚度试样表面与内部淬火残余应力值
根据表2中数据代入公式(4),计算得到各项系数,得到板材淬火心部最大拉应力为:
σ淬火内部=0.0021·H-0.7188·|σ淬火表面|1.8308 (7)
利用钻孔法或X射线衍射法可测得σ淬火表面,代入公式(5),可计算得到σ淬火内部。两个数据可分别对应离差标准化中应力0位置与1位置。据此将应力分布函数还原,可得:
σ(z)=(σ淬火内部-σ淬火表面)·J(z)+σ淬火表面 (8)
其中:
z=h/H (9)
h为板材需预测的点在厚度方向的位置,H为所研究的板材总厚度。即可根据公式(6)、(7)计算板材厚度方向淬火残余应力分布。
Claims (3)
1.一种板材厚度方向淬火残余应力分布预测方法,其特征在于,包括如下步骤:
(1)对选定材料的n组不同厚度规格、相同状态的试样厚度方向残余应力分布进行数值模拟或实际测试分析,得到n组试样应力分布数据,对n组不同厚度及n组试样应力分布数据分别进行离差标准化处理,离差标准化的转换函数如下式:
式中,max为n组不同厚度数据或n组试样应力分布数据最大值,min为n组不同厚度数据或n组试样应力分布数据最小值,x’代表离差标准化处理后的数据,x代表离差标准化处理前的数据;
(2)对n组离差标准化处理结果进行平均,得到离差标准化统计分布,对离差标准化统计分布进行多项式拟合,得到相对于标准化厚度z的残余应力分布函数J(z);
(3)统计分析n组试样表面残余应力绝对值与厚度方向中心位置残余应力绝对值,拟合得到厚度方向中心位置残余应力绝对值与表面残余应力绝对值的对应关系;
(4)对表面残余应力进行测试,测试方法选择钻孔法或X射线衍射法,得到表面残余应力绝对值,并依据步骤(3)得到的厚度方向中心位置残余应力绝对值与表面残余应力绝对值对应关系,计算厚度方向中心位置残余应力;
(5)利用表面残余应力绝对值、厚度方向中心位置残余应力绝对值以及分布函数J(z)预测其他规格样品厚度方向残余应力分布。
2.如权利要求1所述的板材厚度方向淬火残余应力分布预测方法,其特征在于,所述步骤(1)中残余应力分布实际测试分析采用裂纹柔度法或中子衍射法进行,对不同厚度的试样测试位置选择在离差标准化厚度相同的位置进行测试,厚度方向实测点数量不少于9个。
3.如权利要求1所述的板材厚度方向淬火残余应力分布预测方法,其特征在于,所述步骤(1)中n大于等于7个。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711265746.5A CN109870257B (zh) | 2017-12-04 | 2017-12-04 | 一种板材厚度方向淬火残余应力分布预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711265746.5A CN109870257B (zh) | 2017-12-04 | 2017-12-04 | 一种板材厚度方向淬火残余应力分布预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109870257A CN109870257A (zh) | 2019-06-11 |
CN109870257B true CN109870257B (zh) | 2020-12-18 |
Family
ID=66916226
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711265746.5A Active CN109870257B (zh) | 2017-12-04 | 2017-12-04 | 一种板材厚度方向淬火残余应力分布预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109870257B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110907082B (zh) * | 2019-12-10 | 2021-03-16 | 湖北汽车工业学院 | 一种残余应力检测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1793872A (zh) * | 2005-12-29 | 2006-06-28 | 哈尔滨工业大学 | 微小区域残余应力的无损检测方法 |
CN101329204A (zh) * | 2008-07-18 | 2008-12-24 | 清华大学 | 薄膜非均匀应力在线测量的方法及装置 |
CN106370336A (zh) * | 2016-11-08 | 2017-02-01 | 桂林电子科技大学 | 一种柴油机曲轴残余应力检测方法 |
JP2017156263A (ja) * | 2016-03-03 | 2017-09-07 | 株式会社神戸製鋼所 | 残留応力算出方法 |
Family Cites Families (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100465615C (zh) * | 2005-11-25 | 2009-03-04 | 中国科学院金属研究所 | 一种测试薄膜残余应力及其沿层深分布的方法 |
CN101477031A (zh) * | 2008-01-02 | 2009-07-08 | 中国人民解放军装甲兵工程学院 | 一种测定基体残余应力的方法 |
CN101782445B (zh) * | 2010-02-04 | 2011-11-30 | 南京航空航天大学 | 基于斜楔结构的内应力分布快速测量及其误差修正方法 |
CN102169033B (zh) * | 2010-12-08 | 2013-04-03 | 北京科技大学 | 一种铝合金板材内部残余应力定点无损检测方法 |
CN102162755B (zh) * | 2010-12-17 | 2014-08-20 | 罗键 | 一种基于反问题的残余应力检测方法 |
CN102564661B (zh) * | 2011-12-27 | 2013-08-21 | 北京工业大学 | 钛合金表面最大应力和表层应力梯度的计算机测算方法 |
JP6180211B2 (ja) * | 2013-07-12 | 2017-08-16 | 富士フイルム株式会社 | ダイアフラム型共振memsデバイス用基板、ダイアフラム型共振memsデバイス及びその製造方法 |
CN103364336B (zh) * | 2013-07-15 | 2015-05-13 | 江苏科技大学 | 一种测试环焊缝内部双向残余应力全貌的方法 |
CN103926025B (zh) * | 2014-04-13 | 2016-03-09 | 北京工业大学 | 一种用于测量涂层残余应力的试验装置及方法 |
CN104792808B (zh) * | 2014-12-01 | 2017-06-06 | 北京理工大学 | 一种金属基复合材料微观残余应力的检测方法 |
JP6439452B2 (ja) * | 2015-01-08 | 2018-12-19 | いすゞ自動車株式会社 | 高周波焼入れ部材の内部残留応力測定方法 |
CN106768505A (zh) * | 2016-12-02 | 2017-05-31 | 山东省特种设备检验研究院 | 一种q245r材料无损检测应力的方法 |
CN108225636B (zh) * | 2016-12-22 | 2020-11-06 | 有研工程技术研究院有限公司 | 一种铝合金预拉伸板材残余应力的表征方法 |
CN107273630B (zh) * | 2017-06-28 | 2020-05-19 | 华中科技大学 | 一种通过参数反演计算调控残余应力加工参数的方法 |
-
2017
- 2017-12-04 CN CN201711265746.5A patent/CN109870257B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1793872A (zh) * | 2005-12-29 | 2006-06-28 | 哈尔滨工业大学 | 微小区域残余应力的无损检测方法 |
CN101329204A (zh) * | 2008-07-18 | 2008-12-24 | 清华大学 | 薄膜非均匀应力在线测量的方法及装置 |
JP2017156263A (ja) * | 2016-03-03 | 2017-09-07 | 株式会社神戸製鋼所 | 残留応力算出方法 |
CN106370336A (zh) * | 2016-11-08 | 2017-02-01 | 桂林电子科技大学 | 一种柴油机曲轴残余应力检测方法 |
Non-Patent Citations (1)
Title |
---|
《铝合金预拉伸板的残余应力评估与变形预测》;吕海洋;《中国优秀硕士学位论文全文数据库(电子期刊)》;20170831;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109870257A (zh) | 2019-06-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhan et al. | Experimental studies and constitutive modelling of the hardening of aluminium alloy 7055 under creep age forming conditions | |
Xu et al. | Geometry and grain size effects on the forming limit of sheet metals in micro-scaled plastic deformation | |
CN108645704B (zh) | 基于纳米压痕和有限元模拟反演标定金属材料微观本构参数的方法 | |
Newman Jr et al. | Fatigue and crack-growth analyses on 7075-T651 aluminum alloy coupons under constant-and variable-amplitude loading | |
CN108225636B (zh) | 一种铝合金预拉伸板材残余应力的表征方法 | |
CN107908917B (zh) | 一种高强板冲压成形回弹预测方法 | |
CN109933925B (zh) | 一种金属板材的冲压成形性能预测方法 | |
CN107220485B (zh) | 一种适用于多道次压缩的本构模型的建立方法 | |
CN103886125A (zh) | 一种钛合金热复合成形数值模拟方法 | |
McMurtrey et al. | The effect of pit size and density on the fatigue behaviour of a pre‐corroded martensitic stainless steel | |
CN109870257B (zh) | 一种板材厚度方向淬火残余应力分布预测方法 | |
Liu et al. | Machine learning for extending capability of mechanical characterization to improve springback prediction of a quenching and partitioning steel | |
CN112304478A (zh) | 一种基于蠕变轮廓法的残余应力测试方法 | |
CN107609223B (zh) | 一种抗拉强度1200MPa的冷轧双相钢动态变形本构模型建立方法 | |
CN114295491A (zh) | 一种蠕变损伤及变形随时间演化行为的预测方法 | |
Zhou et al. | Multi-objective optimization of stamping forming process of head using Pareto-based genetic algorithm | |
CN108169019B (zh) | 一种准静态塑性压缩应力应变参量的识别方法 | |
Li et al. | Identification of material parameters from punch stretch test | |
CN113053471B (zh) | 一种用于风机主轴无损在线检测布氏硬度的方法 | |
Halilovič et al. | Prediction of elastic strain recovery of a formed steel sheet considering stiffness degradation | |
Chen et al. | Hot deformation behaviour and Hansel-Spittel constitutive model of Cr5 alloy for heavy backup roll | |
CN113946975A (zh) | 一种基于Copula的PNET桥梁体系可靠性评估方法 | |
JP7281708B2 (ja) | 不良品の発生に寄与する装置特定のための製造条件計算装置、製造条件計算方法及び製造条件計算プログラム | |
CN114018832A (zh) | 一种钢表面涂层防护等级的评估方法 | |
CN108520167B (zh) | 一种g102钢受热面的高温寿命快速评估方法和系统 |
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 | ||
TA01 | Transfer of patent application right | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20190628 Address after: 101407 No. 11 Xingke East Street, Yanqi Economic Development Zone, Huairou District, Beijing Applicant after: Research Institute of engineering and Technology Co., Ltd. Address before: No. 2, Xinjie street, Xicheng District, Beijing, Beijing Applicant before: General Research Institute for Nonferrous Metals |
|
GR01 | Patent grant | ||
GR01 | Patent grant |