CN109992853B - 一种表面强化金属部件的残余应力场数值重建方法 - Google Patents

一种表面强化金属部件的残余应力场数值重建方法 Download PDF

Info

Publication number
CN109992853B
CN109992853B CN201910206511.1A CN201910206511A CN109992853B CN 109992853 B CN109992853 B CN 109992853B CN 201910206511 A CN201910206511 A CN 201910206511A CN 109992853 B CN109992853 B CN 109992853B
Authority
CN
China
Prior art keywords
temperature
stress
value
residual
residual stress
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
Application number
CN201910206511.1A
Other languages
English (en)
Other versions
CN109992853A (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.)
Southwest Jiaotong University
Original Assignee
Southwest Jiaotong University
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 Southwest Jiaotong University filed Critical Southwest Jiaotong University
Priority to CN201910206511.1A priority Critical patent/CN109992853B/zh
Publication of CN109992853A publication Critical patent/CN109992853A/zh
Application granted granted Critical
Publication of CN109992853B publication Critical patent/CN109992853B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

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

Abstract

本发明公开一种表面强化金属部件的残余应力场数值重建方法，基于商业CAE软件，利用温度差使强化层区域产生压缩应力，对材料或部件强化层和其他区域分别设置不同的温度值，从而用于模拟各种不同表面强化工艺所产生的残余压缩应力场，并且根据实测部件残余应力数据进行调整，使其与实测数据的残余应力梯度分布相吻合，利于后续的数值仿真计算。本发明可得到较为精确的结果，通用性强，效率高，经济性好，大大简化了仿真和设计过程，便于工程应用中快速、经济地推证出符合实测数据的残余压缩应力场。

Description

一种表面强化金属部件的残余应力场数值重建方法
技术领域
本发明涉及机械零件表面强化技术领域，具体为一种表面强化金属部件的残余应力场数值重建方法。
背景技术
表面强化能够有效改善机械零件和构件表面性能，提高疲劳强度和耐磨性能。承受载荷的零件表面通常处于最大应力状态，并在不同的介质环境中工作。因此，机械零部件的疲劳失效和破坏也大多发生在表面或从表面起始。表面强化能够在零件表层引入一定的残余压应力，抑制疲劳裂纹的萌生及扩展，从而有效提高部件的抗疲劳开裂性能，大幅延长工作寿命。但表面残余压缩应力的存在也给结构强度评价和断裂力学仿真分析带来了巨大困难。
工程实际中部件表面强化方式不尽相同，根据强化层残余应力的实测和经验数据，将残余应力梯度分布在有限元软件中准确重建，是断裂力学、动力学及结构强度数值仿真的关键与前提。
进行有限元仿真计算时，首先要把因表面强化形成的残余压缩应力植入有限元模型中。而工程实际中部件表面强化方式不尽相同，如喷丸、滚压、激光冲击、表面感应淬火等，如果按照相应强化工艺进行模拟，模拟过程就会十分复杂，有效性也很难保证。
发明内容
针对上述问题，本发明的目的在于提供一种能够简化表面残余压缩应力仿真过程，便于工程应用中快速、经济地推证出符合实测数据的残余压缩应力场的表面强化金属部件的残余应力场数值重建方法。技术方案如下：
一种表面强化金属部件的残余应力场数值重建方法，包括以下步骤：
步骤1：使用商业CAD三维软件，根据实际部件的结构尺寸进行三维几何建模，并以STEP格式导出模型文件；
步骤2：将所述模型文件导入商业CAE软件，采用完全六面体单元C3D8R进行网格剖分，并对强化层区域进行网格细分，以inp格式输出网格有限元模型；
步骤3：在网格有限元模型中设置有限元分析所需的基本材料属性，添加材料的线胀系数属性；
步骤4：先将整个模型的初始温度设置为室温，再将强化层区域的温度升高到温度值T high，其余区域的温度仍保持室温不变，从而用温度差异在高温区域产生的压应力来模拟的残余压应力；
步骤5：获取输入的温度值T high与输出的应力值的对应关系，通过调整输入的温度值T high来调整计算得到的应力值，使其与试验测试得到的残余应力值相符；当计算得到的应力分布与试验测试的应力分布相近时，即为合适的T high；保持这一温度场不变，获得部件表面强化层的应力应变场。
进一步的，步骤5之后还包括：将强化层根据深度分成不同区域，在每个区域分别设置一个温度值T high1,T high2,T high3…，从而得到接近真实强化效果的压应力梯度分布。
本发明的有益效果是：本发明不使用复杂的数学模型或迭代公式，比较简便，通用性强，效率高，经济性好，大大简化了仿真和设计过程，适用于工程应用；不必根据实际表面强化工艺进行复杂有限元模拟，而是基于温度场导致的变形机制，简化了表面残余压缩应力仿真过程，便于工程应用中快速、经济地推证出符合实测数据的残余压缩应力场；可用于模拟各种不同工艺所产生的残余应力，通过细化网格和调整强化层不同区域的温度值，可得到较为精确的结果；同时，不受部件形状的制约，方板、圆柱、T型接头等其他形状均可有效模拟。
附图说明
图1为板状表面强化部件的三维建模及网格剖分。
图2为温度场设置示意图。
图3为板状试样表面残余压缩应力应变场及其梯度分布。
图4为残余压缩应力场有限元模拟与实测数据对比。
图5为圆柱表面强化部件有限元模型及其温度场设置示意图。
图6为圆柱试样表面残余压缩应力应变场及其梯度分布。
图7为圆柱试样残余压缩应力场有限元模拟与实测数据对比。
图8为T型焊缝有限元模型。
图9为相同载荷及边界条件下焊缝强化前后应力云图对比。
具体实施方式
下面结合附图和具体实施例对本发明做进一步详细说明。本发明提出一种方便实用和普适性佳的强化层残余压缩应力场数值重建方法。具体步骤如下:
步骤1：使用商业CAD三维软件，根据实际部件的结构尺寸进行详细的三维几何建模，并以STEP格式导出模型文件。
步骤2：如图1所示，使用商业CAE软件，采用完全六面体单元C3D8R进行网格剖分，并对强化层区域进行网格细分，具体网格细化规模依据试验数据和所要求的仿真精度而定，以inp格式输出网格有限元模型。
步骤3：在网格有限元模型中设置有限元分析所需的基本材料属性后，添加材料的线胀系数属性。
步骤4：为模型设置两个分析步，第一步中设置整个模型的初始温度为室温（通常为20°C）；在第二步中将强化层区域的温度升高到某个温度值T high（例如120°C），其余区域的温度仍保持室温（20°C）不变。这种温度差异会使高温区域产生压应力，相当于所要模拟的残余压应力。
步骤5：经过几次试算，将仿真模拟值与经验或者实测值进行比较，找到合适的T high，保持这一温度场不变，可获得部件表面强化层的应力应变场。
获取输入的温度值T high与输出的应力值的对应关系，通过调整输入的温度值T high来调整计算得到的应力值，使其与试验测试得到的残余应力值相符；当计算得到的应力分布与试验测试的应力分布相近时，即为合适的T high。
步骤6：通常为了精确模拟表面强化部件的残余应力场，可进一步将强化层根据深度分成不同区域，在每个区域分别设置一个温度值（ T high1,T high2,T high3…），如图2所示；从而得到更准确的接近真实强化效果的压应力梯度分布，如图3所示。
实施例1：一板状试样其形状如图1所示，试样上下表面经过表面强化处理，划分网格时在其强化层区域网格划分较为细密。将网格模型导入有限元分析软件，在试样强化层和非强化层设置不同的温度值，利用温度差使强化层区域产生压缩应力。为使结果更精确，可将强化层不同区域设置一个温度梯度，如图2所示。图3为有限元模拟的残余压缩应力梯度分布。在设置温度时，要根据试验测试值进行试算以找到合适的T high，最后得到与实测残余应力分布相吻合的模拟残余应力场，如图4所示。
实施例2：具体步骤与例1相同，只是将板状试样改为圆柱形试样，圆柱试样外表面为表面强化层。此例证明发明中所述方法不受部件形状的限制。
实施例3：为一T型焊接接头，其焊缝强化区域如图8所示。焊接结构在疲劳载荷作用下，多在焊缝处萌生疲劳裂纹致部件断裂失效，因此工程中会在焊缝部位进行强化处理，以提高其疲劳性能。图9所示，可在焊缝区域和部件其他区域设置温度差，模拟焊缝区域强化后的残余压缩应力，发明所示方法对于焊接结构这样较为复杂的结构也同样适用。

Claims (2)

1.一种表面强化金属部件的残余应力场数值重建方法，其特征在于，包括以下步骤：
步骤1：使用商业CAD三维软件，根据实际部件的结构尺寸进行三维几何建模，并以STEP格式导出模型文件；
步骤2：将所述模型文件导入商业CAE软件，采用完全六面体单元C3D8R进行网格剖分，并对强化层区域进行网格细分，以inp格式输出网格有限元模型；
步骤3：在网格有限元模型中设置有限元分析所需的基本材料属性，添加材料的线胀系数属性；
步骤4：先将整个模型的初始温度设置为室温，再将强化层区域的温度升高到温度值T high，其余区域的温度仍保持室温不变，从而用温度差异在高温区域产生的压应力来模拟的残余压应力；
步骤5：获取输入的温度值T high与输出的应力值的对应关系，通过调整输入的温度值T high来调整计算得到的应力值，使其与试验测试得到的残余应力值相符；当计算得到的应力分布与试验测试的应力分布相近时，即为合适的T high；保持这一温度场不变，获得部件表面强化层的应力应变场。
2.根据权利要求1所述的表面强化金属部件的残余应力场数值重建方法，其特征在于，步骤5之后还包括：将强化层根据深度分成不同区域，在每个区域分别设置一个温度值T high1,T high2,T high3…，从而得到接近真实强化效果的压应力梯度分布。
CN201910206511.1A 2019-03-19 2019-03-19 一种表面强化金属部件的残余应力场数值重建方法 Active CN109992853B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910206511.1A CN109992853B (zh) 2019-03-19 2019-03-19 一种表面强化金属部件的残余应力场数值重建方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910206511.1A CN109992853B (zh) 2019-03-19 2019-03-19 一种表面强化金属部件的残余应力场数值重建方法

Publications (2)

Publication Number Publication Date
CN109992853A CN109992853A (zh) 2019-07-09
CN109992853B true CN109992853B (zh) 2022-04-29

Family

ID=67130470

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910206511.1A Active CN109992853B (zh) 2019-03-19 2019-03-19 一种表面强化金属部件的残余应力场数值重建方法

Country Status (1)

Country Link
CN (1) CN109992853B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112464528B (zh) * 2020-11-10 2022-02-11 中国核动力研究设计院 一种测量焊接管道内表面及内部残余应力的方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102601167A (zh) * 2012-03-23 2012-07-25 天津大学 一种矫正薄板焊接失稳变形的超声喷丸法及其应用
CN104899388A (zh) * 2015-06-17 2015-09-09 卢伟 一种空间钢结构环境荷载作用下的结构安全评估方法
CN106048143A (zh) * 2016-07-12 2016-10-26 广东工业大学 一种航空发动机叶片边缘预变形激光喷丸强化的方法
CN107354409A (zh) * 2017-07-07 2017-11-17 中南大学 一种温度场均匀的铝合金构件时效成形方法
CN107545110A (zh) * 2017-09-01 2018-01-05 北京强度环境研究所 一种动态应力加速寿命试验剖面编制方法
CN108491652A (zh) * 2018-03-29 2018-09-04 中国计量大学 一种基于abaqus的焊接残余应力测量方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8414563B2 (en) * 2007-12-31 2013-04-09 Deka Products Limited Partnership Pump assembly with switch

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102601167A (zh) * 2012-03-23 2012-07-25 天津大学 一种矫正薄板焊接失稳变形的超声喷丸法及其应用
CN104899388A (zh) * 2015-06-17 2015-09-09 卢伟 一种空间钢结构环境荷载作用下的结构安全评估方法
CN106048143A (zh) * 2016-07-12 2016-10-26 广东工业大学 一种航空发动机叶片边缘预变形激光喷丸强化的方法
CN107354409A (zh) * 2017-07-07 2017-11-17 中南大学 一种温度场均匀的铝合金构件时效成形方法
CN107545110A (zh) * 2017-09-01 2018-01-05 北京强度环境研究所 一种动态应力加速寿命试验剖面编制方法
CN108491652A (zh) * 2018-03-29 2018-09-04 中国计量大学 一种基于abaqus的焊接残余应力测量方法

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
A Guide to Finite Element Simulations of Thermal Barrier Coatings;Martin Baker;《https://link.springer.com/article/10.1007/s11666-017-0592-z》;20170706;1146–1160 *
Finite-Element Analysis and Experimental Test for a Capped-Die Flip-Chip Package Design;Yuci Shen;《 IEEE Transactions on Components, Packaging and Manufacturing Technology ( Volume: 6, Issue: 9, Sept. 2016)》;20160826;1308 – 1316 *
斜轧管材生产中内部组织有限元模拟和预测;双远华;《中国有色金属学报》;20040108;238-242 *
滚压强化的残余应力的数值仿真及工艺分析;王生武等;《计算力学学报》;20081215;113-118 *
激光温喷丸强化航空轻质合金的振动模态及疲劳延寿机理;孟宪凯;《中国博士学位论文全文数据库 (工程科技Ⅱ辑)》;20180115;C031-23 *
箱形结构的焊接残余应力仿真研究;高宗凯等;《机械工程师》;20181110(第11期);40-42,50 *
考虑淬火残余压应力的空心车轴无损探伤周期研究;刘宇轩;《中国力学学会 中国材料科学学会》;20161101;118 *
钢管混凝土拱计算合龙温度试验研究;林春姣等;《广西大学学报(自然科学版)》;20100820(第04期);601-609 *

Also Published As

Publication number Publication date
CN109992853A (zh) 2019-07-09

Similar Documents

Publication Publication Date Title
CN111189699B (zh) 一种基于纳米压痕实验的晶体塑性材料参数反演识别的实现方法
Garcia et al. Fatigue crack growth in residual stress fields
Sun et al. Numerical simulations of the fatigue damage evolution at a fastener hole treated by cold expansion or with interference fit pin
Correia et al. Local unified probabilistic model for fatigue crack initiation and propagation: application to a notched geometry
Sermage et al. Multiaxial creep–fatigue under anisothermal conditions
Herz et al. Numerical simulation of plasticity induced fatigue crack opening and closure for autofrettaged intersecting holes
Xu et al. Modeling of anchor bolt pullout in concrete based on a heterogeneous assumption
CN106407614A (zh) 一种结合纳米压痕与有限元获取焊缝组织力学参数的方法
Yin et al. Development of a new 3D model for the prediction of residual stress and fracture behaviour in Ti-6Al-4V after ultrasonic peening treatment
CN109992853B (zh) 一种表面强化金属部件的残余应力场数值重建方法
Machniewicz Fatigue crack growth prediction models for metallic materials Part II: Strip yield model–choices and decisions
Houghton et al. Multiaxial fretting fatigue testing and prediction for splined couplings
Abdollahi et al. Numerical and experimental study of ratcheting in cold expanded plate of Al‐alloy 2024‐T 3 in double shear lap joints
Zhuang et al. Multipass low-plasticity burnishing induced residual stresses: three-dimensional elastic-plastic finite element modelling
Alshoaibi et al. Finite element simulation of crack growth path and stress intensity factors evaluation in linear elastic materials
Castagne et al. Application of a damage model to an aluminum alloy
Abdollahi et al. Ratcheting behaviour of stainless steel 316L with interference fitted holes in low-cycle fatigue region
Liu et al. Linking process, structure, and property in additive manufacturing applications through advanced materials modelling
Scorza et al. Tension failure assessment at lug hole edges
Mastrone et al. Ductile damage model of an alluminum alloy: Experimental and numerical validation on a punch test
Li et al. Modification of the Johnson–Cook model for metal at a wide range of strain rates and application in the dynamic response of honeycomb panels
Manoucherifar et al. 3D FE analysis of shotpeening process for simulation and research on shotpeening process parameters
Xiao et al. High cycle fatigue life assessment of notched components with induced compressive residual stress
Zhao et al. Effect of parametric variations on the local compression deformation of aluminum foam sandwich panels
Canto et al. Finite element simulation of ceramic powder isostatic pressing process using material parameters for uniaxial compaction

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