最新录用 | 当期目录 | 过刊浏览 | 热点文章 | 虚拟专辑 | 阅读排行 | 下载排行 | 引用排行 | 期刊订阅

金属学报, 2022, 58(12): 1633-1644 DOI: 10.11900/0412.1961.2022.00099

研究论文

GH4169合金拘束相关的疲劳裂纹萌生寿命

郭昊函1, 杨杰,1, 刘芳2, 卢荣生3

1.上海理工大学 能源与动力工程学院 上海 200093

2.上海理工大学 机械工程学院 上海 200093

3.华东理工大学 承压系统与安全教育部重点实验室 上海 200237

Constraint Related Fatigue Crack Initiation Life of GH4169 Superalloy

GUO Haohan1, YANG Jie,1, LIU Fang2, LU Rongsheng3

1.School of Energy and Power Engineering, University of Shanghai for Science and Technology, Shanghai 200093, China

2.School of Mechanical Engineering, University of Shanghai for Science and Technology, Shanghai 200093, China

3.Key Laboratory of Pressure Systems and Safety, Ministry of Education, East China University of Science and Technology, Shanghai 200237, China

通讯作者: 杨 杰，yangjie@usst.edu.cn，主要从事金属材料损伤与断裂等方面的研究

责任编辑: 李海兰

收稿日期: 2022-03-07   修回日期: 2022-04-15  

基金资助: 国家自然科学基金项目(51975378)
上海市浦江人才计划项目(21PJD047)

Corresponding authors: YANG Jie, associate professor, Tel: (021)55272320, E-mail:yangjie@usst.edu.cn

Received: 2022-03-07   Revised: 2022-04-15  

Fund supported: National Natural Science Foundation of China(51975378)
Shanghai Pujiang Program(21PJD047)

作者简介 About authors

郭昊函，男，1996年生，硕士生

摘要

以GH4169镍基高温合金为研究对象，基于低周疲劳率相关的晶体塑性本构，引入累积能量耗散和累积塑性滑移2种疲劳指示因子作为疲劳裂纹萌生判据，对不同微缺口深度和长度下的疲劳裂纹萌生寿命进行研究。并基于统一拘束参数Ap，进一步考察拘束与疲劳裂纹萌生寿命的关联。结果表明：2种疲劳指示因子均可较好地预测疲劳裂纹萌生寿命。随着微缺口深度的增加，疲劳裂纹萌生寿命逐渐减少；随着微缺口长度的增加，疲劳裂纹萌生寿命逐渐增加。在不同的微缺口深度和长度下，疲劳裂纹萌生寿命均与Ap存在线性关系。可以根据该线性关系，确定拘束相关的疲劳裂纹萌生寿命。

关键词： 拘束; 疲劳裂纹; 萌生寿命; 晶体塑性; 累积能量耗散; 累积塑性滑移

Abstract

Nickel-based GH4169 superalloy is used as turbine disc material in aeroengines because of its good oxidation resistance, good formability, weldability, and high strength. However, turbine disc fatigue failure will inevitably occur in onerous service environments and after a long operation time. To ensure the safety and reliability of aeroengines, the fatigue damage behavior and fatigue life of GH4169 superalloy need to be studied. Constraint is an important factor affecting the fatigue fracture behavior of materials, because changing it will impact the fatigue behavior. To achieve a long service life and high reliability of aeroengines, fatigue and constraint effects must be researched. However, there are only limited studies on the effect of constraint on fatigue crack initiation time. In this study, a crystal plasticity constitutive model based on low cycle fatigue rate correlation was applied to the GH4169 superalloy. Two fatigue indicators, namely the cumulated energy dissipation and cumulated plastic slip, were introduced as fatigue crack initiation criteria to study the fatigue crack initiation time for different micro-notch depths and lengths. In addition, the relationship between constraint and fatigue crack initiation life was further investigated using the unified constraint parameter Ap. The results showed that both cumulated energy dissipation and cumulated plastic slip can accurately predict the fatigue crack initiation time. With the increase in micro-notch depth, the fatigue crack initiation time decreased, while it increased with the increase in micro-notch length. A linear relationship between the fatigue crack initiation time and Ap under different micro-notch depths and lengths was observed. Based on this relationship, the constraint related to the fatigue crack initiation time can be determined.

Keywords： constraint; fatigue crack; initiation life; crystal plasticity; cumulated energy dissipation; cumulated plastic slip

PDF (3733KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

郭昊函, 杨杰, 刘芳, 卢荣生. GH4169合金拘束相关的疲劳裂纹萌生寿命[J]. 金属学报, 2022, 58(12): 1633-1644 DOI:10.11900/0412.1961.2022.00099

GUO Haohan, YANG Jie, LIU Fang, LU Rongsheng. Constraint Related Fatigue Crack Initiation Life of GH4169 Superalloy[J]. Acta Metallurgica Sinica, 2022, 58(12): 1633-1644 DOI:10.11900/0412.1961.2022.00099

GH4169镍基高温合金具有良好的抗氧化性、稳定的热加工性、优良的可焊接性和高强度等特点，在航空航天领域常被用作制造涡轮盘的材料[1,2]。然而，在恶劣的服役环境和长时间工作下，涡轮盘将不可避免地发生疲劳失效事故[3~6]。为了保障航空发动机的安全性和可靠性，需对GH4169合金的疲劳损伤行为和疲劳寿命展开深入研究。

在疲劳损伤行为方面，尚德广等[7,8]对GH4169合金在多轴载荷下的力学响应和疲劳损伤进行了研究，指出疲劳寿命和循环硬化或软化特性不仅取决于加载路径，还取决于加载顺序和加载参数。Sun等[9,10]提出了基于临界平面法的多轴疲劳损伤参数用于计算恒定高温下的纯疲劳损伤，并将循环稳定时最大等效应力的一半当量为蠕变应力，以评估多轴蠕变疲劳损伤。Prasad等[11]研究表明，加载波形对GH4169疲劳载荷下的损伤程度和滞后回线形态有显著影响，且在较低应变速率下GH4169合金存在应变时效现象。在疲劳寿命方面，Walz和Riesch-Oppermann[12]基于断裂力学理论和一阶可靠性分析方法，建立了考虑涡轮盘质量和寿命的随机模型。Wei和Yang[13]针对GH4169材料，发展了基于疲劳损伤参数的寿命预测模型，可以预测不同应变幅下的低周疲劳寿命，预测结果在± 2倍误差带范围内。Ayyappan等[14]给出了剩余寿命预测流程，并通过± 3σ (σ代表标准差)应力-应变曲线的统计分析得到了涡轮盘的概率剩余寿命。

然而，上述工作主要在宏观尺度下展开。随着研究的深入，人们发现仅凭宏观尺度下的研究还不足以完全理解疲劳行为，微观尺度下疲劳行为的研究工作相继出现。在实验方面，Krupp[15]研究表明第二相粒子的弹性模量和硬化模量与基体存在较大差异，在γ″相强化的GH4169合金中差异性尤为明显。An等[16]通过对GH4169合金进行不同时间的时效处理，研究了材料微观组织演变与低周疲劳寿命的相关性。研究表明，当时效时间增加时，γ″相转变为δ相，材料的抗疲劳能力显著降低。Alexandre等[17]就晶粒尺寸对裂纹萌生机制影响进行了研究，发现晶粒处于5~10 μm时疲劳裂纹更容易萌生于第二相粒子，晶粒大于150 μm时则会通过穿晶方式滑移萌生。Huang等[18]也得到了相同的研究结果。在有限元数值计算方面，Hu等[19]基于晶体塑性有限元法对GH4169合金在蠕变疲劳载荷下的裂纹扩展行为展开研究，证实晶粒细化和晶界处δ相的存在使GH4169合金在650℃下弱界面的数量快速增长，从而导致沿晶裂纹的萌生和扩展。Guerchais等[20]基于晶体塑性有限元法分析了缺陷对高周疲劳下裂纹萌生的影响，研究表明缺陷尺寸的增加会导致扭转、拉伸疲劳极限间比值的增加。Zhang等[21,22]采用扩展有限元(XFEM)和晶体塑性有限元相结合的方法对镍基高温合金的短裂纹扩展现象展开研究，通过采用各个滑移系上的独立剪切应变作为驱动力成功模拟了单晶fcc结构镍基高温合金中沿滑移方向的裂纹扩展路径。

其中，在对微观疲劳裂纹萌生进行研究时，如何定义裂纹的萌生判据成为研究的焦点。Manonukul等[23]和Dunne等[24]使用塑性滑移作为裂纹萌生的判据，准确预测了多晶镍基合金的疲劳裂纹萌生点、疲劳裂纹扩展路径及塑性滑移带的位置和方向，并与实验得到的结果高度吻合。Guo等[25,26]以能量耗散作为阈值，发展了基于能量耗散的双线法评估疲劳极限。Yuan等[27]基于晶体塑性有限元法对GH4169合金的低周疲劳循环塑性行为进行了研究，结果表明累积塑性滑移和累积能量耗散均能作为疲劳裂纹萌生的判据。

无论在宏观尺度或微观尺度下，拘束均是影响材料疲劳行为的重要因素。拘束可以看作结构对裂尖区材料塑性变形的阻碍，与三轴应力状态密切相关。影响拘束的因素包括裂纹尺寸、试样或结构几何、材料性能失配等方面。通常将裂纹尺寸和试样或结构几何所引起的拘束效应称为几何拘束，将材料性能失配所引起的拘束效应称为材料拘束。拘束的变化会引起疲劳裂纹扩展行为的改变。在现有标准中，疲劳裂纹扩展速率往往采用标准试样在实验室测得，但是在相同应力水平下，工程结构实际工况与实验室条件往往并不一致，与此同时，不同拘束试样和结构的裂纹扩展速率也会存在差别。一些学者观察到高拘束会降低疲劳裂纹扩展速率[28~30]，也有学者观察到了相反的现象[31]。一些学者发现拘束效应对Paris区和近门槛值区的疲劳裂纹扩展速率都有影响[29,32]，也有研究[30]表明拘束效应的影响只存在于近门槛值区。Yang等[33]在最近研究中表明拘束对疲劳裂纹扩展速率的影响与加载方式有关。当用较快的疲劳裂纹扩展速率评定具有较慢疲劳裂纹扩展速率结构的完整性时会产生过于保守的评定结果，从而导致不必要的维修或报废，产生较大的经济损失；反之，当用其评定具有更快疲劳裂纹扩展速率结构的完整性时，可能得到非保守的结果，从而产生危险。

为了准确理解材料的疲劳行为，需要在研究中兼顾疲劳与拘束效应，但目前尚缺乏有关拘束效应对GH4169合金疲劳裂纹萌生寿命的影响研究。鉴于此，本工作基于晶体塑性有限元法，分别选用累积塑性滑移和累积能量耗散2种疲劳指示因子作为疲劳裂纹萌生判据，对不同微缺口深度和长度下(即不同拘束下)的GH4169合金疲劳裂纹萌生寿命进行研究。并基于统一拘束参数Ap，进一步探究疲劳裂纹萌生中的拘束效应，以建立拘束与疲劳裂纹萌生寿命的关联。

1 晶体塑性理论

采用Asaro等[34,35]提出的率相关晶体塑性本构，该本构通过大变形理论分析晶粒之间的局部相互作用。当多晶材料受到外加载荷时，每一个单晶体的总变形梯度张量可用 F 表示[34]：

F=Fe⋅Fp

式中， Fe为晶格畸变和刚性转动所产生的变形梯度，表征弹性部分； Fp为晶体沿着滑移方向的均匀剪切所产生的变形梯度，表征塑性部分。

对于由剪切产生的塑性变形，塑性速率梯度 Lp可以表示为与塑性滑移率γ˙α相关的方程[35]：

Lp=F˙pFp-1=∑αΝγ˙αmα⊗nα

式中，mα和nα分别代表α滑移系上的滑移方向以及滑移面的法向，符号⊗代表张量的并矢积，N代表总滑移系的个数。

在经典唯象本构模型中，通常使用临界剪切应力τcα作为滑移系的状态变量，γ˙α是关于剪切应力τα与τcα的方程[35]：

γ˙α=fτα , τcα , ⋯

其中，剪切应力与当前应力状态相关[35]：

τα=S∶mα⊗nα

式中，S为第二类Piola-Kirchhoff应力。

在第α个滑移系中，τcα与γ˙α之间的关系可以采用如下方程描述[36]：

γ˙α=γ˙0sgnτα-χατα-χαgαn

式中，γ˙0为参考应变率，n为材料属性相关的率敏感系数，本工作中，γ˙0取0.001，n取100[27]；χα 为背应力，用来描述材料疲劳载荷下的循环变形行为；gα 为α滑移系的当前强度参量，定义为[36]：

g˙α=∑βhαβγ˙β

式中，γ˙β为β滑移系的塑性滑移率；hαβ为塑性硬化矩阵，与自硬化矩阵hαα相关[36]：

hαβ=qhαα
hαα=h0sech2h0γαg∞-τ0

式中，q为潜硬化系数，通常取1~1.4，本工作中q取1[36]；参数h0、τ0和g∞分别为初始硬化模量、临界分解剪切应力和饱和应力。

本模型通过引入背应力参量χ˙α进而描述疲劳的随动硬化现象，其非线性的演化方程为[27]：

χ˙α=Cγ˙α-Dχαγ˙α

式中，C代表直接硬化模量，D代表动态回复模量。

晶体塑性理论框架通过用户材料自定义子程序UMAT嵌入ABAQUS软件中，实现对材料疲劳行为的描述。此外，为了表征材料的损伤程度，并确定疲劳裂纹萌生寿命，引入2种疲劳指示因子——累积能量耗散W、累积塑性滑移P，作为微观疲劳裂纹萌生判据[37]。

W可通过τα和γ˙α定义[37]：

W=∑α∫0tταγ˙αdt

P可通过 Lp定义[37]：

P˙=23Lp∶Lp1/2
P=∫0tP˙dt

W和P均与循环周次呈线性关系[37]，根据测定的临界累积能量耗散Wcrit、临界累积塑性滑移Pcrit和单周循环能量耗散Wcyc、单周循环塑性滑移Pcyc，通过 式(13)和(14)计算得到疲劳裂纹萌生寿命Ni。

Ni=WcritWcyc
Ni=PcritPcyc

2 方案设计和有限元建模

2.1 材料

选用航空发动机涡轮盘材料GH4169镍基高温合金，其化学成分(质量分数，%)为[27]：C 0.82，Al 0.44，Si 0.18，Ti 1.16，Cr 19.55，Fe 18.93，Nb 5.19，Mb 2.74，Ni余量。其晶体塑性本构参数为[27]：弹性常数C11 = 266 GPa，C12 = 114 GPa，C44 = 76 GPa；n = 100；γ˙0 = 0.001；h0 = 200 MPa；g∞ = 650 MPa；τ0 = 315 MPa；C = 39500 MPa，D = 475 MPa。Wcrit = 89159 MJ / m3，Pcrit = 118.2[37]。

2.2 代表体积单元建立

通过Voronoi tessellation (VT)法[38]生成GH4169镍基高温合金的代表体积单元(RVE)，如图1所示。图中RVE整体尺寸为200 μm × 200 μm，共包含178个晶粒。对于镍基合金材料，当晶粒数量超过150个，模拟结果会有较好的稳定性[39]。晶粒采用随机分布策略，平均晶粒尺寸为(15 ± 1) μm，与实际晶粒尺寸相符。所有晶粒均为包含12个滑移系的fcc结构，并根据GH4169电子背散射衍射(EBSD)结果，赋予所有晶粒随机取向。

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1   GH4169镍基高温合金的代表体积单元(RVE)

Fig.1   A representative volume element (RVE) of GH4169 nickel base superalloy


2.3 方案设计

为了分析拘束对疲劳裂纹萌生寿命的影响，在图1模型上边界处设置不同尺寸的微缺口，如图2所示。其中，a和b分别代表微缺口的深度和长度，分别通过改变a (a = 15、45、75和105 μm)和b (b = 15、45、75和105 μm)，即1倍、3倍、5倍和7倍晶粒尺寸得到不同拘束度。当改变a时，b固定为15 μm；当改变b时，a固定为15 μm。为进一步研究应变幅对疲劳裂纹萌生寿命的影响，对每种不同的微缺口模型，分别在3种应变幅(Δεt = 0.6%、0.8%和1.0%)下进行有限元数值计算。

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2   模型加载与尺寸示意图

Fig.2   Schematic of model loading and geometry (a—micro-notch depth, b—micro-notch length, RP—reference point)


2.4 有限元建模

在RVE建立的基础上，通过Python语言读取其拓扑信息，并在ABAQUS软件中生成基于二维平面应变的多晶模型。分别读取每个晶粒的取向，随材料参数一起赋给每个晶粒。在模型左右两侧建立参考点，如图2所示，参考点分别与左右两边进行耦合，通过在参考点处施加位移载荷的方式控制加载，加载波形为三角波。为了消除刚体位移的影响，对底边沿y方向的位移进行约束，且允许其在x方向上自由移动。模型网格采用四节点平面应变单元(CPE4)，为了避免网格尺寸对计算结果产生影响，分别对模型划分不同尺寸(2、1和0.5 μm)的网格进行网格无关性验证，如图3所示。由图可见，当网格尺寸降到1 μm时，计算结果不再随网格尺寸的变化而发生改变，故本工作所用的网格尺寸为1 μm。

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3   基于3种不同网格尺寸所得到的应力-应变曲线和塑性滑移云图

Fig.3   Stress-strain curves and plastic slips based on three different mesh sizes


3 拘束对疲劳裂纹萌生的影响

3.1 不同微缺口深度

3.1.1 累积能量耗散

应变幅Δεt = 0.6%时，4个具有不同微缺口深度(a = 15、45、75和105 μm)的模型在循环10 cyc后的累积能量耗散如图4所示(仅截取了微缺口根部处)。可见，在相同的循环周次下，随着微缺口深度的不断增加，累积能量耗散不断增加。

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4   应变幅为0.6%时不同微缺口深度模型在循环10 cyc后的累积能量耗散

Fig.4   Accumulated energy dissipation illustrations of models with different micro-notch depths when the strain amplitude Δεt = 0.6% and the cycle number is 10 cyc (unit: MJ/m3)

(a) a = 15 μm (b) a = 45 μm

(c) a = 75 μm (d) a = 105 μm


在3种应变幅(Δεt = 0.6%、0.8%、1.0%)下对4个不同微缺口深度模型在各循环周次下的累积能量耗散进行计算，可得累积能量耗散与循环周次间的关系，如图5所示。可见，随着循环周次的增加，累积能量耗散线性增加，即对于具有特定微缺口尺寸的模型，在特定的应变幅加载下，能量耗散的累积速率是恒定的，这与文献[37]的结果一致。随着微缺口深度的增加，能量耗散的累积速率逐渐加快(图5中直线的斜率逐渐增加)，说明高拘束对于能量耗散的累积有明显的促进作用。但随着拘束的增加，能量耗散的累积速率变化逐渐减小，有趋于稳定的趋势。这说明拘束虽然对能量耗散的累积有促进作用，但终归有限。应变幅对累积速率的影响更大。随着应变幅的增加，不仅累积能量耗散和能量耗散累积速率在逐渐增加(以a = 105 μm微缺口深度为例，在循环40 cyc后，应变幅为0.6%时，累积能量耗散为2750 MJ/m3，能量耗散累积速率为68 MJ/(m3·cyc)；应变幅为0.8%时，累积能量耗散为5194 MJ/m3，能量耗散累积速率为129 MJ/(m3·cyc)；应变幅为1.0%时，累积能量耗散为10470 MJ/m3，能量耗散累积速率为261 MJ/(m3·cyc))，且最大累积速率需要更高的拘束度才能达到(应变幅为1.0%时，a = 105 μm时能量耗散累积速率未趋于稳定)。

图5

新窗口打开| 下载原图ZIP| 生成PPT

图5   累积能量耗散与循环周次间的关系

Fig.5   Relationships between accumulated energy dissipation and cycle number

(a) Δεt = 0.6% (b) Δεt = 0.8% (c) Δεt = 1.0%


以上结果可以从滞后回线的变化中得到很好的解释。选取微缺口左下角点，对其在应变幅0.6%、不同微缺口深度、循环第10 cyc时的滞后回线进行考察，如图6所示。由图可见，随着微缺口深度的增加，滞后回线的应力变化范围逐渐增大、应变变化范围基本一致，这与微缺口尺寸变化方向和加载方向有关。此时微缺口变化方向与加载方向垂直，随着微缺口深度增加，微缺口应变基本保持不变，但应力随拘束的增加越来越大。

图6

新窗口打开| 下载原图ZIP| 生成PPT

图6   应变幅为0.6%时微缺口左下角点在循环第10 cyc时的滞后回线

Fig.6   Hysteresis loops of the point at the lower left corner of the micro-notch at Δεt = 0.6% and the 10th cycle


3.1.2 累积塑性滑移

应变幅为0.6%时，4个微缺口深度模型在循环10 cyc后的累积塑性滑移如图7所示。由图可见，与累积能量耗散规律一致，在相同的循环周次下，随着微缺口深度的不断增加，累积塑性滑移不断增加。

图7

新窗口打开| 下载原图ZIP| 生成PPT

图7   应变幅为0.6%时不同微缺口深度模型在循环10 cyc后的累积塑性滑移

Fig.7   Accumulated plastic slip illustrations of models with different micro-notch depths when Δεt = 0.6% and the cycle number is 10 cyc

(a) a = 15 μm (b) a = 45 μm (c) a = 75 μm (d) a = 105 μm


不同应变幅下，不同微缺口深度模型累积塑性滑移与循环周次间的关系如图8所示。由图可见，与累积能量耗散一致，累积塑性滑移与循环周次也呈线性关系，塑性滑移的累积速率也是恒定的。且随着微缺口深度的增加，塑性滑移的累积速率逐渐加快。与能量耗散相比，塑性滑移累积速率随着拘束的增加而趋于稳定的趋势不明显。应变幅的增加会显著提高累积塑性滑移和塑性滑移累积速率(以a = 105 μm微缺口深度为例，在循环40 cyc后，应变幅为0.6%时，累积塑性滑移为2.48，塑性滑移累积速率为0.062 cyc-1；应变幅为0.8%时，累积塑性滑移为4.67，塑性滑移累积速率为0.12 cyc-1；应变幅为1.0%时，累积塑性滑移为8.27，塑性滑移累积速率为0.21 cyc-1)。

图8

新窗口打开| 下载原图ZIP| 生成PPT

图8   累积塑性滑移与循环周次间的关系

Fig.8   Relationships between accumulated plastic slip and cycle number

(a) Δεt = 0.6% (b) Δεt = 0.8% (c) Δεt = 1.0%


3.1.3 寿命对比

分别基于累积能量耗散和累积塑性滑移，通过 式(13)和(14)计算得到不同微缺口深度、不同应变幅下的疲劳裂纹萌生寿命，如图9所示。需要说明的是，因为累积能量耗散和累积塑性滑移均与循环周次呈线性关系， 式(13)和(14)中的Wcyc和Pcyc分别与能量耗散和塑性滑移的累积速率相等，在寿命计算时，为方便起见，分别用能量耗散和塑性滑移的累积速率代替Wcyc和Pcyc。由图可见，随着微缺口深度的增加，疲劳裂纹萌生寿命逐渐降低。应变幅对疲劳裂纹萌生寿命的影响尤为明显，随着应变幅的增加，各微缺口深度下的疲劳裂纹萌生寿命均显著降低。基于累积能量耗散进行预测所得寿命比基于累积塑性滑移预测所得寿命略低。

图9

新窗口打开| 下载原图ZIP| 生成PPT

图9   基于累积能量耗散和累积塑性滑移疲劳裂纹萌生寿命预测对比

Fig.9   Comparisons of predicted fatigue crack initiation life based on accumulated energy dissipation and accumulated plastic slip


3.2 不同微缺口长度

3.2.1 累积能量耗散

应变幅为0.6%时，不同微缺口长度的模型(b = 15、45、75和105 μm)在循环10 cyc后的累积能量耗散如图10所示。由图可见，在相同的循环周次下，随着微缺口长度的不断增加，累积能量耗散不断减少，与微缺口深度变化所得趋势相反。

图10

新窗口打开| 下载原图ZIP| 生成PPT

图10   应变幅为0.6%时不同微缺口长度模型在循环10 cyc后的累积能量耗散

Fig.10   Accumulated energy dissipation illustrations of models with different micro-notch lengths when Δεt = 0.6% and the cycle number is 10 cyc (unit: MJ/m3)

(a) b = 15 μm (b) b = 45 μm (c) b = 75 μm (d) b = 105 μm


3种应变幅下4个微缺口长度模型累积能量耗散与循环周次间的关系如图11所示。可见，与不同微缺口深度相同，对于具有特定微缺口尺寸的模型，在特定的应变幅加载下，随着循环周次的增加，累积能量耗散线性增加，能量耗散的累积速率恒定；应变幅对累积能量耗散和能量耗散累积速率的影响更大，随着应变幅增加，二者均逐渐增加。与不同微缺口深度时的情况不同，随着微缺口长度增加，能量耗散累积速率逐渐降低，且其趋于稳定的趋势不明显。如上结果可以从滞后回线的变化中得到解释。

图11

新窗口打开| 下载原图ZIP| 生成PPT

图11   累积能量耗散与循环周次间的关系

Fig.11   Relationships between accumulated energy dissipation and cycle number

(a) Δεt = 0.6% (b) Δεt = 0.8% (c) Δεt = 1.0%


微缺口左下角点在应变幅0.6%、不同微缺口长度、循环第10 cyc时的滞后回线如图12所示。可见，随着微缺口长度的增加，滞后回线的应变变化范围逐渐减小、应力变化范围基本一致(在b = 105 μm时有所降低)，这与微缺口尺寸变化方向和加载方向一致有关。随着微缺口长度增加，沿加载方向应变被平均，应变变化范围均匀减小，应力变化范围基本不变。

图12

新窗口打开| 下载原图ZIP| 生成PPT

图12   应变幅0.6%时微缺口左下角点在循环第10 cyc后的滞后回线

Fig.12   Hysteresis loops of the point at the lower left corner of the micro-notch at Δεt = 0.6% and the 10th cycle


3.2.2 累积塑性滑移

应变幅为0.6%时，4个微缺口长度模型在循环10 cyc后的累积塑性滑移如图13所示。可见，与累积能量耗散规律一致，在相同的循环周次下，随着微缺口长度的不断增加，累积塑性滑移不断减少。

图13

新窗口打开| 下载原图ZIP| 生成PPT

图13   应变幅为0.6%时不同微缺口长度模型在循环10 cyc后的累积塑性滑移

Fig.13   Accumulated plastic slip illustrations of models with different micro-notch lengths when Δεt = 0.6% and the cycle number is 10 cyc

(a) a = 15 μm (b) a = 45 μm (c) a = 75 μm (d) a = 105 μm


不同应变幅下，不同微缺口长度模型累积塑性滑移与循环周次间的关系如图14所示。可见，累积塑性滑移与循环周次同样呈线性关系，塑性滑移的累积速率同样恒定。且随着微缺口长度的增加，塑性滑移的累积速率逐渐降低。应变幅的增加同样会显著提高累积塑性滑移和塑性滑移累积速率。

图14

新窗口打开| 下载原图ZIP| 生成PPT

图14   累积塑性滑移与循环周次间的关系

Fig.14   Relationships between accumulated plastic slip and cycle number

(a) Δεt = 0.6% (b) Δεt = 0.8% (c) Δεt = 1.0%


3.2.3 寿命对比

分别基于累积能量耗散和累积塑性滑移计算得到不同微缺口长度、不同应变幅下的疲劳裂纹萌生寿命，如图15所示。与3.1.3节相同，在寿命计算时，为方便起见，分别用能量耗散和塑性滑移的累积速率代替Wcyc和Pcyc。由图可见，随着微缺口长度的增加，疲劳裂纹萌生寿命逐渐增加；随着应变幅的增加，疲劳裂纹萌生寿命显著降低。基于累积塑性滑移和累积能量耗散进行预测所得寿命相差不大。

图15

新窗口打开| 下载原图ZIP| 生成PPT

图15   基于累积能量耗散和累积塑性滑移疲劳裂纹萌生寿命预测对比

Fig.15   Comparisons of predicted fatigue crack initiation life based on accumulated energy dissipation and accumulated plastic slip


4 疲劳裂纹萌生寿命与拘束间的关联

可以发现，在不同的微缺口深度和长度(拘束)下，疲劳裂纹萌生寿命有很大的不同。为了对含缺口试样或结构进行准确的结构完整性评价，需建立疲劳裂纹萌生寿命与拘束间的关联，在评价中纳入拘束效应的影响。

在之前的研究中，Yang等[40~42]曾基于等效塑性应变(εp)等值线所围绕区域的面积定义了统一拘束参数Ap：

Αp=ΑPEEQ / Aref

式中，APEEQ为不同拘束试样或结构中εp等值线所围绕区域的面积； Aref为标准试样中εp等值线所围绕区域的面积。该参数已在准静态加载下得到了广泛应用[40~42]，并进一步用于疲劳加载，建立了宏观疲劳裂纹扩展速率与拘束间的统一关联[33]。

为了进一步探究疲劳裂纹萌生中的拘束效应，选择Δεt = 1.0%时，对不同拘束模型在10 cyc循环后微缺口处εp = 0.3等值线(在微缺口εp等值线与加载处εp等值线不相连的情况下，可选取任意εp值)所围绕的面积APEEQ进行计算；并以a = b = 15 μm时微缺口处εp = 0.3等值线所围绕的面积作为参考(Aref)，根据 式(14)计算不同拘束下的Ap。并以累积塑性滑移所得的疲劳裂纹萌生寿命为准，建立疲劳裂纹萌生寿命与拘束间的关联，如图16所示。由图可见，在不同的微缺口深度和长度下，疲劳裂纹萌生寿命均与Ap呈线性关系。且无论是改变微缺口深度或是微缺口长度，疲劳裂纹萌生寿命与Ap存在唯一线性关系，如图17所示。

图16

新窗口打开| 下载原图ZIP| 生成PPT

图16   疲劳裂纹萌生寿命与拘束间的关联

Fig.16   Relationships between fatigue crack initiation life and constraint

(a) different micro-notch depth

(b) different micro-notch length


图17

新窗口打开| 下载原图ZIP| 生成PPT

图17   所有模型疲劳裂纹萌生寿命与拘束间的关联

Fig.17   Relationship between fatigue crack initiation life and constraint of all models


这进一步表明Ap可以适用于疲劳领域，可以将其应用范围从准静态加载推广到疲劳加载。另一方面，结果表明疲劳裂纹萌生寿命在很大程度上取决于等效塑性应变等值线所围绕区域的面积。在之前的研究中，Chikh等[43]试图将塑性区尺寸作为一个力学参数引入疲劳裂纹失效准则，然而所使用的塑性区尺寸多指塑性区的长度，并非面积。从本工作结果来看，等效塑性应变等值线所围绕区域的面积或许可以作为一个合适的参数引入疲劳裂纹失效准则。此外，疲劳裂纹萌生寿命与Ap良好的关联性对于不同拘束下准确的结构完整性评定具有意义。可以根据本关联，确定拘束相关的疲劳裂纹萌生寿命。

5 结论

(1) 随着微缺口深度的增加，能量耗散和塑性滑移的累积速率逐渐加快，累积能量耗散和累积塑性滑移逐渐增加，疲劳裂纹萌生寿命减少；随着微缺口长度的增加，变化趋势相反。这与微缺口尺寸变化方向和加载方向有关。

(2) 与微缺口深度和长度的改变相比，应变幅改变所带来的影响更大。随着应变幅的增加，累积能量耗散和累积塑性滑移均显著增加，疲劳裂纹萌生寿命显著减少。

(3) 在不同微缺口深度下，基于累积能量耗散进行预测所得寿命比基于累积塑性滑移进行预测所得寿命略低；在不同微缺口长度下，基于累积能量耗散和基于累积塑性滑移进行预测所得寿命相差不大。

(4) 在不同的微缺口深度和长度下，疲劳裂纹萌生寿命均与Ap存在线性关系。Ap可以适用于疲劳领域，可尝试作为一个合适的参数引入疲劳裂纹失效准则。

参考文献

Azarbarmas M, Aghaie-Khafri M, Cabrera J M, et al.

Dynamic recrystallization mechanisms and twining evolution during hot deformation of Inconel 718

[J]. Mater. Sci. Eng., 2016, A678: 137

[本文引用: 1]

Chen F, Liu J, Ou H G, et al.

Flow characteristics and intrinsic workability of IN718 superalloy

[J]. Mater. Sci. Eng., 2015, A642: 279

[本文引用: 1]

McDowell D, Dunne F P E.

Microstructure-sensitive computational modeling of fatigue crack formation

[J]. Int. J. Fatigue, 2010, 32: 1521

DOI      URL     [本文引用: 1]

Deng G J, Tu S T, Zhang X C, et al.

Small fatigue crack initiation and growth mechanisms of nickel-based superalloy GH4169 at 650℃ in air

[J]. Eng. Fract. Mech., 2016, 153: 35

DOI      URL    

Cini A, Irving P E.

Development of fatigue cracks from mechanically machined scratches on 2024-T351 aluminium alloy—Part I: Experimentation and fractographic analysis

[J]. Fatigue Fract. Eng. Mater. Struct., 2017, 40: 776

DOI      URL    

Cini A, Irving P E.

Development of fatigue cracks from mechanically machined scratches on 2024-T351 aluminium alloy—Part II: Finite element analysis and prediction method

[J]. Fatigue Fract. Eng. Mater. Struct., 2017, 40: 853

DOI      URL     [本文引用: 1]

Shang D G, Chen J H, Sun G Q, et al.

Fatigue characteristics for GH4169 superalloy under uniaxial/multiaxial loading at high temperature

[J]. Acta Metall. Sin., 2005, 41: 785

[本文引用: 1]

尚德广, 陈建华, 孙国芹 等.

单、多轴混合加载下GH4169合金的高温疲劳特性

[J]. 金属学报, 2005, 41: 785

[本文引用: 1]

利用拉扭薄壁管疲劳试样, 在应变控制拉扭循环加载下对高温合金材料GH4169的多轴循环特性进行了实验研究. 高温疲劳实验过程中, 采用单、多轴混合加载路径对薄壁管疲劳试件进行加载, 通过连续记录拉与扭的应力响应值研究了变幅高温多轴疲劳特性. 结果表明, 在高温低周单、多轴混合疲劳块载荷加载下, 试件的疲劳寿命与应力响应特性不但取决于应变加载路径, 而且与加载路径的排列顺序和加载参数的大小有关.

Shang D G, Sun G Q, Chen J H, et al.

Multiaxial fatigue behavior of Ni-based superalloy GH4169 at 650℃

[J]. Mater. Sci. Eng., 2006, A432: 231

[本文引用: 1]

Sun G Q, Shang D G.

Prediction of fatigue lifetime under multiaxial cyclic loading using finite element analysis

[J]. Mater. Des., 2010, 31: 126

DOI      URL     [本文引用: 1]

Sun G Q, Shang D G, Bao M.

Multiaxial fatigue damage parameter and life prediction under low cycle loading for GH4169 alloy and other structural materials

[J]. Int. J. Fatigue, 2010, 32: 1108

DOI      URL     [本文引用: 1]

Prasad K, Sarkar R, Ghosal P, et al.

High temperature low cycle fatigue deformation behaviour of forged IN718 superalloy turbine disc

[J]. Mater. Sci. Eng., 2013, A568: 239

[本文引用: 1]

Walz G, Riesch-Oppermann H.

Probabilistic fracture mechanics assessment of flaws in turbine disks including quality assurance procedures

[J]. Struct. Saf., 2006, 28: 273

DOI      URL     [本文引用: 1]

Wei D S, Yang X G.

Investigation and modeling of low cycle fatigue behaviors of two Ni-based superalloys under dwell conditions

[J]. Int. J. Press. Vessel. Pip., 2009, 86: 616

DOI      URL     [本文引用: 1]

Ayyappan C, Kumar R, Ramesh P, et al.

Experimental and numerical study to predict residual growth in an aeroengine compressor disc after overspeed

[J]. Procedia Eng., 2013, 55: 625

DOI      URL     [本文引用: 1]

Krupp U. Fatigue Crack Propagation in Metals and Alloys: Microstructural Aspects and Modelling Concepts [M]. Weinheim: Wiley-VCH, 2007: 1

[本文引用: 1]

An J L, Wang L, Liu Y, et al.

The role of δ phase for fatigue crack propagation behavior in a Ni base superalloy at room temperature

[J]. Mater. Sci. Eng., 2017, A684: 312

[本文引用: 1]

Alexandre F, Deyber S, Pineau A.

Modelling the optimum grain size on the low cycle fatigue life of a Ni based superalloy in the presence of two possible crack initiation sites

[J]. Scr. Mater., 2004, 50: 25

DOI      URL     [本文引用: 1]

Huang X Y, Yu H C, Xu M Q, et al.

Experimental investigation on microcrack initiation process in nickel-based superalloy DAGH4169

[J]. Int. J. Fatigue, 2012, 42: 153

DOI      URL     [本文引用: 1]

Hu D Y, Wang X Y, Mao J X, et al.

Creep-fatigue crack growth behavior in GH4169 superalloy

[J]. Front. Mech. Eng., 2019, 14: 369

DOI      URL     [本文引用: 1]

Guerchais R, Saintier N, Morel F, et al.

Micromechanical investigation of the influence of defects in high cycle fatigue

[J]. Int. J. Fatigue, 2014, 67: 159

DOI      URL     [本文引用: 1]

Zhang P, Zhang L, Baxevanakis K P, et al.

Modelling short crack propagation in a single crystal nickel-based superalloy using crystal plasticity and XFEM

[J]. Int. J. Fatigue, 2020, 136: 105594

[本文引用: 1]

Zhang P, Baxevanakis K P, Zhao L G.

An investigation of short crack propagation in a single crystal Ni-based superalloy using crystal plasticity and the extended finite element method

[J]. Procedia Struct. Integr., 2020, 28: 1176

DOI      URL     [本文引用: 1]

Manonukul A, Dunne F P E.

High- and low-cycle fatigue crack initiation using polycrystal plasticity

[J]. Proc. Roy. Soc., 2004, 460A: 1881

[本文引用: 1]

Dunne F P E, Wilkinson A J, Allen R.

Experimental and computational studies of low cycle fatigue crack nucleation in a polycrystal

[J]. Int. J. Plast., 2007, 23: 273

DOI      URL     [本文引用: 1]

Guo Q, Guo X L, Fan J L, et al.

An energy method for rapid evaluation of high-cycle fatigue parameters based on intrinsic dissipation

[J]. Int. J. Fatigue, 2015, 80: 136

DOI      URL     [本文引用: 1]

Guo Q, Guo X L.

Research on high-cycle fatigue behavior of FV520B stainless steel based on intrinsic dissipation

[J]. Mater. Des., 2016, 90: 248

DOI      URL     [本文引用: 1]

Yuan G J, Zhang X C, Chen B, et al.

Low-cycle fatigue life prediction of a polycrystalline nickel-base superalloy using crystal plasticity modelling approach

[J]. J. Mater. Sci. Technol., 2020, 38, 28

DOI      URL     [本文引用: 5]

Varfolomeev I, Luke M, Moroz S.

Experimental and numerical investigations of fatigue crack growth in various specimen geometries

[J]. Procedia Eng., 2010, 2: 1829

DOI      URL     [本文引用: 1]

Hutař P, Seitl S, Knésl Z.

Effect of constraint on fatigue crack propagation near threshold in medium carbon steel

[J]. Comput. Mater. Sci., 2006, 37: 51

[本文引用: 1]

Hutař P, Seitl S, Kruml T.

Effect of specimen geometry on fatigue crack propagation in threshold region

[A]. Proceedings of the 12th International Conference on Fracture [C]. icf12 Ottawa, 2009: 4: 2914

[本文引用: 2]

Tong J.

T-stress and its implications for crack growth

[J]. Eng. Fract. Mech., 2002, 69: 1325

DOI      URL     [本文引用: 1]

Seitl S, Hutař P.

Fatigue-crack propagation near a threshold region in the framework of two-parameter fracture mechanics

[J]. Mater. Technol., 2007, 41: 135

[本文引用: 1]

Yang J, Guo G L, Lu R S, et al.

Numerical modelling of a new FCP model and a correlation of the FCP rate with the constraint

[J]. Int. J. Fatigue, 2022, 163: 107036

[本文引用: 2]

Asaro R J.

Micromechanics of crystals and polycrystals

[J]. Adv. Appl. Mech., 1983, 23: 115

[本文引用: 2]

Asaro R J, Rice J R.

Strain localization in ductile single crystals

[J]. J. Mech. Phys. Solids, 1977, 25: 309

DOI      URL     [本文引用: 4]

Peirce D, Asaro R J, Needleman A.

Material rate dependence and localized deformation in crystalline solids

[J]. Acta Metall., 1983, 31: 1951

DOI      URL     [本文引用: 4]

Zhong F.

Researches on tensile properties and fatigue behavior of Ni-based alloy GH4169 based on crystal plasticity finite element

[D]. Shanghai: East China University of Science and Technology, 2017

[本文引用: 6]

钟飞.

基于晶体塑性有限元的镍基合金GH4169拉伸性能及疲劳行为研究

[D]. 上海: 华东理工大学, 2017

[本文引用: 6]

Barbe F, Decker L, Jeulin D, et al.

Intergranular and intragranular behavior of polycrystalline aggregates. Part 1: F.E. model

[J]. Int. J. Plast., 2001, 17: 513

DOI      URL     [本文引用: 1]

Lin B, Zhao L G, Tong J.

A crystal plasticity study of cyclic constitutive behaviour, crack-tip deformation and crack-growth path for a polycrystalline nickel-based superalloy

[J]. Eng. Fract. Mech., 2011, 78: 2174

DOI      URL     [本文引用: 1]

Yang J, Wang G Z, Xuan F Z, et al.

Unified characterisation of in-plane and out-of-plane constraint based on crack-tip equivalent plastic strain

[J]. Fatigue Fract. Eng. Mater. Struct., 2013, 36: 504

DOI      URL     [本文引用: 2]

Yang J, Wang G Z, Xuan F Z, et al.

Unified correlation of in-plane and out-of-plane constraints with fracture toughness

[J]. Fatigue Fract. Eng. Mater. Struct., 2014, 37: 132

DOI      URL    

Yang J, Wang G Z, Xuan F Z, et al.

Unified correlation of in-plane and out-of-plane constraint with fracture resistance of a dissimilar metal welded joint

[J]. Eng. Fract. Mech., 2014, 115: 296

DOI      URL     [本文引用: 2]

Chikh B O, Imad A, Benguediab M.

Influence of the cyclic plastic zone size on the propagation of the fatigue crack in case of 12NC6 steel

[J]. Comput. Mater. Sci., 2008, 43: 1010

DOI      URL     [本文引用: 1]

/

〈 〉


版权所有 © 《金属学报》编辑部
地址: 沈阳市文化路72号, 中国科学院金属研究所(110016)
电话: +86- 024-23971286 　E-mail: jsxb@imr.ac.cn
本系统由北京玛格泰克科技发展有限公司设计开发 技术支持：support@magtech.com.cn