CN108153928A - 一种含夹杂物粉末高温合金裂纹萌生寿命预测方法 - Google Patents

一种含夹杂物粉末高温合金裂纹萌生寿命预测方法 Download PDF

Info

Publication number
CN108153928A
CN108153928A CN201711111270.XA CN201711111270A CN108153928A CN 108153928 A CN108153928 A CN 108153928A CN 201711111270 A CN201711111270 A CN 201711111270A CN 108153928 A CN108153928 A CN 108153928A
Authority
CN
China
Prior art keywords
crack initiation
powder metallurgy
life
metallurgy superalloy
cycle
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
CN201711111270.XA
Other languages
English (en)
Other versions
CN108153928B (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.)
AECC Beijing Institute of Aeronautical Materials
Original Assignee
AECC Beijing Institute of Aeronautical Materials
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 AECC Beijing Institute of Aeronautical Materials filed Critical AECC Beijing Institute of Aeronautical Materials
Priority to CN201711111270.XA priority Critical patent/CN108153928B/zh
Publication of CN108153928A publication Critical patent/CN108153928A/zh
Application granted granted Critical
Publication of CN108153928B publication Critical patent/CN108153928B/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
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Investigating And Analyzing Materials By Characteristic Methods (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明属于寿命预测技术,涉及一种含夹杂物粉末高温合金裂纹萌生寿命预测方法。寿命预测的步骤如下:定义描述夹杂物周围材料裂纹萌生进程的损伤变量场D(xi,N);定义表征参量;定义循环特征参数R;定义损伤变量场变化率方程;裂纹萌生寿命预测。本发明提出了一种含夹杂物粉末高温合金裂纹萌生寿命预测方法,建立了一套可以预测含夹杂物粉末高温合金裂纹萌生寿命的方法,能预测粉末高温合金制件的使用寿命,建立了粉末高温合金制件的疲劳失效性能判据,进而为航空发动机部件设计过程中提供了安全性评估依据。

Description

一种含夹杂物粉末高温合金裂纹萌生寿命预测方法
技术领域
本发明属于寿命预测技术,涉及一种含夹杂物粉末高温合金裂纹萌生寿命预测方法。
背景技术
航空发动机的设计和制造是现代航空器的关键技术,其发展水平代表着国家航空工业的整体水平。涡轮盘是航空发动机的关键部件,工作时承受较高的离心载荷、热载荷、气动载荷、振动载荷和环境介质的腐蚀氧化作用,对材料性能要求严苛。粉末高温合金是航空发动机涡轮盘的理想材料,它具有非常优异的性能:粉末细小,冷速极快,合金成分均匀,无宏观偏析,因而制件性能稳定,热加工变形性能较好。然而,由于其工艺特点,颗粒中的夹杂物不可避免。这些夹杂物对粉末高温合金的力学性能有显著影响,易成为裂纹源,大大降低材料的疲劳性能,严重影响和制约合金的发展应用。目前,国内国际对于含夹杂物粉末高温合金裂纹萌生机理和寿命预测方法尚缺乏系统全面的研究。因此,无论是在粉末高温合金材料的失效机理研究还是在其工程应用中,建立描述含夹杂物粉末高温合金的裂纹扩展寿命预测理论模型是非常必要的。“镍基粉末高温合金中夹杂物导致裂纹萌生和扩展行为的研究”,曾燕屏等,《材料工程》,2005(3):10-13,采用扫描电镜原位拉伸和原位疲劳的方法,跟踪观察了人工植入Al2O3夹杂物的镍基粉末高温合金P/M Rene95中夹杂物导致裂纹萌生、扩展乃至断裂的过程。扫描电镜原位拉伸和原位疲劳的方法一方面成本较高,另一方面该方法只能通过断口反推寿命而无法进行预测。
发明内容
本发明的目的是:提出一种含夹杂物粉末高温合金裂纹萌生寿命预测方法,以便建立一套可以预测含夹杂物粉末高温合金裂纹萌生寿命的方法,预测粉末高温合金制件的使用寿命,建立粉末高温合金制件的疲劳失效性能判据,进而为航空发动机部件设计过程中提供安全性评估依据。
本发明的技术方案是:一种含夹杂物粉末高温合金裂纹萌生寿命预测方法,其特征在于,寿命预测的步骤如下:
1、定义描述夹杂物周围材料裂纹萌生进程的损伤变量场D(xi,N):
式[1]中,E为粉末高温合金的弹性模量,ED为局部细观缺陷造成损伤后的宏观等效模量;损伤变量场D的取值范围从0至1,当其为0时表示材料没有开始损伤,当其为1时表示有宏观裂纹萌生;xi表示空间位置坐标,N为循环周次;此变量场为一标量场,是循环周次和空间位置的函数,其在材料不同位置的分布随时间变化;
2、定义表征参量Y为如下应变分量的函数:
式中,W为应变能密度,μ为泊松比,εx,εy,εz,γxy,γyz,γzx为应变分量;
3、定义循环特征参数R:
若在一个循环过程中,表征参量只有一个峰值,
若在一个循环过程中,表征参量有两个峰值,
4、定义损伤变量场变化率方程:
式中,a和m为材料参数,需要利用粉末高温合金标准试样的高温低周疲劳实验数据,由最小二乘法拟合确定;
5、裂纹萌生寿命预测:
将上述损伤变量场变化率方程嵌入有限元软件Ansys或Abaqus中,在有限元软件中安照实际结构的几何尺寸和材料参数建立相应的模型,并施加载荷边界条件,求解时在每一载荷周次步长内逐单元计算损伤场变量并据此更新单元的材料参数,迭代计算直至D在局部单元达到1时跳出计算,统计步长就可以预测出材料的裂纹萌生寿命。
本发明的优点是:提出了一种含夹杂物粉末高温合金裂纹萌生寿命预测方法,建立了一套可以预测含夹杂物粉末高温合金裂纹萌生寿命的方法,能预测粉末高温合金制件的使用寿命,建立了粉末高温合金制件的疲劳失效性能判据,进而为航空发动机部件设计过程中提供了安全性评估依据。
附图说明
图1a是指一个循环过程中,表征参量只有一个峰值,横坐标为时间,纵坐标为表征参量。
图1b是指一个循环过程中,表征参量有两个峰值,横坐标为时间,纵坐标为表征参量。
图2是有限元软件中建立的含圆形缺陷的二维模型。
图3是标准实验数据和在对数坐标下的拟合曲线,横坐标为最大应变载荷对数值,纵坐标为寿命对数值。
具体实施方式
下面对本发明做进一步详细说明。一种含夹杂物粉末高温合金裂纹萌生寿命预测方法,其特征在于,寿命预测的步骤如下:
1、定义描述夹杂物周围材料裂纹萌生进程的损伤变量场D(xi,N):
式[1]中,E为粉末高温合金的弹性模量,ED为局部细观缺陷造成损伤后的宏观等效模量;损伤变量场D的取值范围从0至1,当其为0时表示材料没有开始损伤,当其为1时表示有宏观裂纹萌生;xi表示空间位置坐标,N为循环周次;此变量场为一标量场,是循环周次和空间位置的函数,其在材料不同位置的分布随时间变化;
2、定义表征参量Y为如下应变分量的函数:
式中,W为应变能密度,μ为泊松比,εx,εy,εz,γxy,γyz,γzx为应变分量;
3、定义循环特征参数R:
若在一个循环过程中,表征参量只有一个峰值,
若在一个循环过程中,表征参量有两个峰值,
4、定义损伤变量场变化率方程:
式中,a和m为材料参数,需要利用粉末高温合金标准试样的高温低周疲劳实验数据,由最小二乘法拟合确定;
5、裂纹萌生寿命预测:
将上述损伤变量场变化率方程嵌入有限元通用软件Ansys或Abaqus中,在有限元软件中安照实际结构的几何尺寸和材料参数建立相应的模型,并施加载荷边界条件,求解时在每一载荷周次步长内逐单元计算损伤场变量并据此更新单元的材料参数,迭代计算直至D在局部单元达到1时跳出计算,统计步长就可以预测出材料的裂纹萌生寿命。
本发明的工作原理是:在粉末高温合金承受交变载荷时,由于夹杂物的存在必然导致夹杂物周边产生局部应力应变集中,使得此区域内实际应力应变大于远场的名义应力应变。故夹杂物的存在加速了裂纹萌生过程,缩短了材料使用的疲劳寿命。
实施例
待预测寿命的构件为一个尺寸为长宽均是10毫米的平板形构件,设在其中心存在一个半径R为100微米的圆形氧化铝夹杂物,其几何形式如图2所示。构件基体材料为FGH95粉末高温合金,其弹性模量为E=193GPa,泊松比为0.3;中心夹杂物设为粉末高温合金常见的氧化铝硬球状夹杂物,其弹性模量E2=400GPa,泊松比0.3。设该构件受到远场单轴拉伸交变载荷作用,拉伸最大应变为0.8%,循环特征为脉动载荷,即最大远场拉伸应变为0.8%,而最小应变为0。
1、定义描述夹杂物周围材料裂纹萌生进程的损伤变量场D(xi,N):
2、定义表征参量Y为如下应变分量的函数:
3、定义循环特征参数R:
根据发明内容中循环特征参数R的定义及实际加载,其取值应为R=1;
4、定义损伤变量场变化率方程:
标准实验数据和在对数坐标下的拟合曲线如图3所示。实施例的损伤参数已经由不含夹杂物的标准光滑疲劳试件的实验数据拟合确定为m=5.594;a=5.206×10-8
5、裂纹萌生寿命预测:
基于上述步骤,将实施例所述构件在数值计算软件中建立相应模型并计算其裂纹萌生寿命。实施例使用的是有限元法的通用软件Abaqus。
几何建模:将待计算裂纹萌生寿命的构件按设计几何尺寸与材料性能在计算软件中建立相应的几何模型,如图2所示;
划分网络:在上述方法建立的几何模型基础上划分网络,为保证计算精度,实施例中选择了8节点二次平面单元,网格划分为四边形结构网格。为控制计算规模并保证关键区域计算精度,采取在缺陷处局部加密,然后经过渡区到基体逐渐稀疏的方式划分网格;
边界条件施加:在上述方法划分网络的基础上施加相应边界条件。在实施例中,边界条件为远场施加位移,根据试件尺寸和载荷情况计算为上下边界沿拉伸方向各位移0.04mm;
计算结果提取:在上述各方法实施后计算全场应力应变的结果,并在后处理程序中调取应变结果;
迭代计算萌生寿命结果:上述方法包括在每一计算步长内,提取全场应力应变数据并逐单元按式[2]计算表征参量Y,再根据式[6]公式计算损伤变量场变化率dD/dN,可以得到此刻构件内全场的损伤值;由每个单元的损伤值根据式[1]更新单元的材料参数,进入下一步长计算。迭代计算直至损伤变量场D在裂纹萌生处的值达到1时跳出计算。统计计算迭代次数就可以计算出材料的裂纹萌生寿命了。在实施例中,按此计算的裂纹萌生寿命为N=6500周次。
本实施例,具有如下有益效果:
(1)在本发明中,可以根据已知结构形式和载荷情况以及夹杂物特征,定量计算粉末高温合金的裂纹萌生寿命;
(2)利用本发明给出的计算公式嵌入数值计算软件进行二次开发,可以方便快捷地得到裂纹萌生过程的损伤变量场变化过程图示;
(3)不必实施过多实验,仅通过理论计算即可得到裂纹萌生寿命结果,节省设计成本;
(4)本发明可以为粉末高温合金构件在含夹杂物情况下的安全性评估提供参考。

Claims (1)

1.一种含夹杂物粉末高温合金裂纹萌生寿命预测方法,其寿命预测的步骤如下:
1.1、定义描述夹杂物周围材料裂纹萌生进程的损伤变量场D(xi,N):
式[1]中,E为粉末高温合金的弹性模量,ED为局部细观缺陷造成损伤后的宏观等效模量;损伤变量场D的取值范围从0至1,当其为0时表示材料没有开始损伤,当其为1时表示有宏观裂纹萌生;xi表示空间位置坐标,N为循环周次;此变量场为一标量场,是循环周次和空间位置的函数,其在材料不同位置的分布随时间变化;
1.2、定义表征参量Y为如下应变分量的函数:
式中,W为应变能密度,μ为泊松比,εx,εy,εz,γxy,γyz,γzx为应变分量;
1.3、定义循环特征参数R:
若在一个循环过程中,表征参量只有一个峰值,
若在一个循环过程中,表征参量有两个峰值,
1.4、定义损伤变量场变化率方程:
式中,a和m为材料参数,需要利用粉末高温合金标准试样的高温低周疲劳实验数据,由最小二乘法拟合确定;
1.5、裂纹萌生寿命预测:
将上述损伤变量场变化率方程嵌入有限元软件Ansys或Abaqus中,在有限元软件中安照实际结构的几何尺寸和材料参数建立相应的模型,并施加载荷边界条件,求解时在每一载荷周次步长内逐单元计算损伤场变量并据此更新单元的材料参数,迭代计算直至D在局部单元达到1时跳出计算,统计步长就可以预测出材料的裂纹萌生寿命。
CN201711111270.XA 2017-11-10 2017-11-10 一种含夹杂物粉末高温合金裂纹萌生寿命预测方法 Active CN108153928B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711111270.XA CN108153928B (zh) 2017-11-10 2017-11-10 一种含夹杂物粉末高温合金裂纹萌生寿命预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711111270.XA CN108153928B (zh) 2017-11-10 2017-11-10 一种含夹杂物粉末高温合金裂纹萌生寿命预测方法

Publications (2)

Publication Number Publication Date
CN108153928A true CN108153928A (zh) 2018-06-12
CN108153928B CN108153928B (zh) 2021-02-09

Family

ID=62468924

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711111270.XA Active CN108153928B (zh) 2017-11-10 2017-11-10 一种含夹杂物粉末高温合金裂纹萌生寿命预测方法

Country Status (1)

Country Link
CN (1) CN108153928B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109163990A (zh) * 2018-08-28 2019-01-08 中国科学院金属研究所 一种轴向加载高周疲劳裂纹萌生寿命的测量方法
CN109255202A (zh) * 2018-11-01 2019-01-22 上海工程技术大学 一种用于机械构件疲劳裂纹萌生寿命的预估方法
CN112163352A (zh) * 2020-08-13 2021-01-01 西安建筑科技大学 中碳钢超细晶棒材3d-spd轧制损伤预测方法及模型
CN115034092A (zh) * 2022-08-09 2022-09-09 中国航发北京航空材料研究院 含夹杂物的粉末高温合金低周疲劳寿命的预测方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101126799A (zh) * 2007-09-28 2008-02-20 中国人民解放军装甲兵工程学院 利用铁磁材料表面杂散磁场信号监测疲劳损伤的方法
CN101509855A (zh) * 2009-03-17 2009-08-19 华东理工大学 一种高温管道损伤及寿命的在线预测方法
CN102426137A (zh) * 2011-08-31 2012-04-25 中国船舶重工集团公司第七二五研究所 一种轴向加载低周疲劳裂纹萌生的试验测量方法
CN103778292A (zh) * 2014-01-23 2014-05-07 北京航空航天大学 一种热振联合载荷下bga焊点疲劳寿命预测方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101126799A (zh) * 2007-09-28 2008-02-20 中国人民解放军装甲兵工程学院 利用铁磁材料表面杂散磁场信号监测疲劳损伤的方法
CN101509855A (zh) * 2009-03-17 2009-08-19 华东理工大学 一种高温管道损伤及寿命的在线预测方法
CN102426137A (zh) * 2011-08-31 2012-04-25 中国船舶重工集团公司第七二五研究所 一种轴向加载低周疲劳裂纹萌生的试验测量方法
CN103778292A (zh) * 2014-01-23 2014-05-07 北京航空航天大学 一种热振联合载荷下bga焊点疲劳寿命预测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
许捷 等: "夹杂物对粉末高温合金材料裂纹萌生的影响", 《西北工业大学学报》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109163990A (zh) * 2018-08-28 2019-01-08 中国科学院金属研究所 一种轴向加载高周疲劳裂纹萌生寿命的测量方法
CN109255202A (zh) * 2018-11-01 2019-01-22 上海工程技术大学 一种用于机械构件疲劳裂纹萌生寿命的预估方法
CN109255202B (zh) * 2018-11-01 2023-06-06 上海工程技术大学 一种用于机械构件疲劳裂纹萌生寿命的预估方法
CN112163352A (zh) * 2020-08-13 2021-01-01 西安建筑科技大学 中碳钢超细晶棒材3d-spd轧制损伤预测方法及模型
CN112163352B (zh) * 2020-08-13 2024-04-19 西安建筑科技大学 中碳钢超细晶棒材3d-spd轧制损伤预测方法及模型
CN115034092A (zh) * 2022-08-09 2022-09-09 中国航发北京航空材料研究院 含夹杂物的粉末高温合金低周疲劳寿命的预测方法

Also Published As

Publication number Publication date
CN108153928B (zh) 2021-02-09

Similar Documents

Publication Publication Date Title
CN108153928A (zh) 一种含夹杂物粉末高温合金裂纹萌生寿命预测方法
Sadowski et al. Numerical modelling of two-phase ceramic composite response under uniaxial loading
Neukamm et al. Considering damage history in crashworthiness simulations
Wan et al. Multi-scale damage mechanics method for fatigue life prediction of additive manufacture structures of Ti-6Al-4V
Du et al. Effects of grain size on the high-cycle fatigue behavior of IN792 superalloy
Kiarasi et al. A review on functionally graded porous structures reinforced by graphene platelets
Zhou et al. Effect of skew angle of holes on the tensile behavior of a Ni-base single crystal superalloy
CN108627406B (zh) 一种高强金属材料变幅超高周疲劳寿命预测方法
Flores et al. A fatigue analysis of a hydraulic Francis turbine runner
CN111339703A (zh) 大应变条件下材料硬化行为的虚拟预测方法
CN103914593A (zh) 层状复合材料力学行为的虚拟预测方法
CN113758815A (zh) 一种搅拌摩擦焊构件变幅疲劳寿命预测方法
Wu et al. Micromechanical modeling of damage evolution and fracture behavior in particle reinforced metal matrix composites based on the conventional theory of mechanism-based strain gradient plasticity
Li et al. Effect of phase coarsening on the mechanical properties of alloys: I. Numerical simulations
Wang et al. Investigation on experimental load spectrum for high and low cycle combined fatigue test
CN107219136A (zh) 一种考虑不同减薄率的超塑性成形件疲劳寿命预测方法
Han et al. Research on the effect of micro-voids on the deformation behavior and crack initiation lifetime of titanium alloy under cyclic loading by crystal plasticity finite element method
Wang et al. Numerical simulation and experimental studies on elastic-plastic fatigue crack growth
CN113782104A (zh) 镍基单晶合金模型构建方法、装置、存储介质及设备
BHADAURIA et al. Effect of stress triaxiality on yielding of anisotropic materials under plane stress condition
Sunder A unified model of fatigue kinetics based on crack driving force and material resistance
Nozhnitsky et al. Numerical simulation of spin testing for turbo machine disks using energy-based fracture criteria
Liu et al. Fatigue behavior and dislocation substructures for 6063 aluminum alloy under nonproportional loadings
Archut et al. Estimation of damage value of SAE4340 steel and 5A02 aluminium alloy
Chen et al. The “Inverse Hall-Petch” effect on the impact response of single crystal copper

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