CN106644490A - 一种涡轮榫接结构高低周复合疲劳裂纹扩展寿命预测方法 - Google Patents
一种涡轮榫接结构高低周复合疲劳裂纹扩展寿命预测方法 Download PDFInfo
- Publication number
- CN106644490A CN106644490A CN201611266906.3A CN201611266906A CN106644490A CN 106644490 A CN106644490 A CN 106644490A CN 201611266906 A CN201611266906 A CN 201611266906A CN 106644490 A CN106644490 A CN 106644490A
- Authority
- CN
- China
- Prior art keywords
- max
- ccf
- crack
- lcf
- under
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M15/00—Testing of engines
- G01M15/02—Details or accessories of testing apparatus
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明涉及一种涡轮榫接结构高低周复合疲劳裂纹扩展寿命预测方法,步骤:(1)建立考虑裂纹闭合效应的裂纹扩展模型;(2)确定高周载荷应力强度因子模型;(3)确定低周载荷应力强度因子模型;(4)判断当前循环是否为高周载荷,若是,执行(5),否则执行(6);(5)计算复合疲劳下裂纹增量,判断下一循环是否为低周疲劳载荷,若是,获取当前循环最大裂尖张开位移,执行(7),否则执行(8);(6)计算低周载荷下裂尖张开位移变化量;(7)根据(5)和(6)输入的结果,计算残余裂尖张开位移及裂纹增量;(8)更新裂纹长度,若最大应力强度因子小于断裂韧度,重复步骤(4)到(7),否则计算完成,得到裂纹扩展寿命。
Description
技术领域
本发明是一种针对航空发动机涡轮榫接结构高低周复合疲劳裂纹扩展寿命的预测方法,它是一种能够高低周耦合作用、裂纹闭合效应和小时间尺度的计算方法,属于航空航天发动机技术领域。
背景技术
涡轮盘是航空发动机为数不多的关键件之一,涡轮盘普遍采用枞树形榫槽结构与叶片榫头联接。涡轮榫接同时承受叶片离心力、热载荷等构成的低周疲劳载荷和由气动载荷诱发叶片横向振动及自身振动的小幅值、高频率的高周疲劳载荷,即高低周复合疲劳载荷。由于结构和受载的双重复杂性,不仅多型在役航空发动机长期受到涡轮榫接裂纹故障困扰,某些新型发动机也多次发生涡轮榫接裂纹故障,严重威胁飞机的飞行安全。为此,建立可准确预测涡轮盘榫槽高低周复合疲劳裂纹扩展寿命的分析方法研究具有重要意义。
目前对于高低周复合疲劳裂纹扩展寿命的预测方法主要利用传统方法分别计算高周载荷和低周载荷下的应力强度因子,基于线性累积损伤理论对裂纹扩展进行分析,这具有明显的局限性:(1)没有考虑高周载荷与低周载荷的耦合作用,难以保证预测精度;(2)传统方法将高频载荷处理成叠加在低周载荷上的等幅静态载荷,未考虑结构振动特性对高周应力强度因子的影响,未体现结构振动特性随裂纹长度的增加而发生的改变。
现有文献Hu D,Yang Q,Liu H,et al.Crack closure effect and crack growthbehavior in GH2036superalloy plates under combined high and low cycle fatigue[J].Int J Fatigue.2017,95:90-103从试验角度出发研究了GH2036材料高低周复合疲劳行为,并对其裂纹扩展寿命进行了预测,但其研究成果仅针对实验室情况下标准试件,未考虑涡轮榫接部件的结构特征,未对涡轮榫接结构进行瞬态分析,也未考虑载荷历程对高低周复合疲劳裂纹扩展寿命的影响,其结果针对材料较为单一,缺乏一定的应用性。此外,由于涡轮榫接结构在裂纹扩展过程中,会改变榫接的接触状态,造成应力的重新分布,在计算涡轮榫接结构的裂纹扩展寿命时,必须考虑这方面的影响。
发明内容
本发明技术解决方案:克服现有技术的不足,提供一种涡轮榫接结构高低周复合疲劳裂纹扩展寿命预测方法,考虑裂了纹闭合效应和结构振动特性,保证了预测精度。
本发明技术解决方案:一种涡轮榫接结构高低周复合疲劳裂纹扩展寿命预测方法,概括起来,主要包括:考虑裂纹闭合效应裂纹扩展寿命模型的建立,高低周复合疲劳裂纹扩展寿命模型的建立。
实现步骤如下:
(1)根据不同应力比下低周疲劳裂纹扩展试验,计算涡轮盘榫接材料的张开应力强度因子Kop,有效应力强度因子ΔKeff,并确定残余裂尖张开位移δres与张开应力强度因子Kop的关系,并建立考虑裂纹闭合效应的疲劳裂纹扩展寿命模型,为后续计算疲劳载荷下的裂纹扩展增量提供基础;
(2)对不同裂纹长度下涡轮榫接结构施加高低周复合疲劳载荷,然后进行瞬态分析,根据分析结果拟合确定复合疲劳载荷作用下,高周载荷应力强度因子模型ΔKCCF,max=f(a),和最大裂尖张开位移ΔδCCF,max与高周最大应力强度因子ΔKCCF,max关系模型;
(3)对不同裂纹长度下涡轮榫接结构施加低周疲劳载荷,然后进行瞬态分析,根据分析结果,拟合确定低周载荷应力强度因子模型KLCF,max=f(a),和最大裂尖张开位移δLCF,max与低周最大应力强度因子KLCF,max模型;
(4)设定涡轮榫接结构所受复合疲劳载荷的一个循环起始于卸载,终止于加载到最大载荷,在某一榫齿当前循环裂纹长度an下,判断当前复合疲劳载荷下是否为高周疲劳载荷,若是则执行步骤(5),否则执行步骤(6);
(5)由该榫齿裂纹长度an和步骤(2)确定的高周载荷应力强度因子模型,计算当前循环下的高周复合应力强度因子ΔKCCF,max,计算高低周复合载荷下的裂纹增量Δan,判断下一循环状态是否为低周疲劳载荷,若是低周疲劳载荷,由步骤(2)和(3)确定的最大裂尖张开位移模型,计算当前循环复合疲劳载荷下最大裂尖张开位移δCCF,max,残余裂尖张开位移δCCF,res,张开应力强度因子Kop,CCF,执行步骤(7),否则执行步骤(8);
(6)由该榫齿裂纹长度an和步骤(3)确定的低周应力强度因子模型,计算低循环应力强度因子KLCF及卸载过程中裂尖张开位移的变化量Δδ;
(7)根据步骤(5)计算的复合疲劳载荷下最大裂尖张开位移δCCF,max和步骤(6)计算的裂尖张开位移变化量Δδ,计算低周载荷下残余裂尖张开位移δLCF,res=δCCF,max-Δδ,根据步骤(1)建立的考虑裂纹闭合效应的疲劳裂纹扩展寿命模型,计算低周载荷下的裂纹增量Δan;
(8)更新该榫齿裂纹长度an+1=an+Δan,判断最大应力强度因子是否小于断裂韧度,若小于断裂韧度,再次重复步骤(4)至(7);反之若大于等于断裂韧度,计算完成,得到最终裂纹长度,根据初始和最终裂纹长度即可得出高低周复合疲劳裂纹扩展循环寿命N。
所述步骤(1)中,残余裂尖张开位移δres与张开应力强度因子Kop的关系为:
Kop=(a1×R+b1)δres+c1×R+d1
其中,a1、b1、c1、d1分别为与材料和温度有关的常数,可通过试验数据或数值模拟结果数据拟合获得;R为载荷的应力比;
考虑裂纹闭合效应的裂纹扩展寿命模型为:
其中,da/dN为裂纹扩展速率;C0,n0为与材料有关常数,可通过试验数据或数值模拟结果拟合获得;ΔKeff为有效应力强度因子幅值,且ΔKeff=Kmax-Kop,Kmax为最大应力强度因子。
所述步骤(2)中,高周载荷应力强度因子ΔKCCF,max=KCCF,max-KLCF,max,其中KCCF,max为榫接结构复合疲劳载荷最大应力强度因子,
拟合得到的高周载荷应力强度因子模型ΔKCCF,max=f(a)形式如下:
△KCCF,max=a2a2+b2a+c2
其中,a2、b2、c2分别为与材料和温度有关的常数,可通过试验数据或数值模拟结果数据拟合获得,a为裂纹长度。
高周载荷下最大裂尖张开位移ΔδCCF,max与高周最大应力强度因子ΔKCCF,max的关系模型为:
其中,C1,n1为与材料有关常数,可通过数值模拟结果拟合获得。
所述步骤(3)中,拟合得到的低周载荷应力强度因子模型KLCF,max=f(a)形式如下:
KLCF,max=a3a2+b3a+c3
其中,a3、b3、c3分别为与材料和温度有关的常数,可通过试验数据或数值模拟结果数据拟合获得,a为裂纹长度。
拟合得到的低周载荷下最大裂尖张开位移δLCF,max与低周最大应力强度因子KLCF,max的关系模型为:
其中,C2,n2为与材料有关常数,可通过数值模拟结果拟合获得。
所述步骤(5)中,计算高低周复合疲劳载荷下裂纹增量Δan时,采用以下方法:
其中,ΔKCCF,max根据步骤(2)建立的高周应力强度因子模型确定;C0,n0由步骤(1)确定;
计算当前循环下复合疲劳载荷下最大裂尖张开位移δCCF,max时,应基于步骤(2)和步骤(3)确定的最大裂尖张开位移模型进行计算,
δCCF,max=δLCF,max+△δCCF,max;
计算当前循环下复合疲劳载荷下残余裂尖张开位移δCCF,res,采用以下方法:
其中,KCCF,max,KCCF,min为本循环最大和最小复合疲劳载荷应力强度因子;Kop,pre为上一循环张开应力强度因子;C1,n1由步骤(2)确定;
计算当前循环下复合疲劳载荷下张开应力强度因子Kop,CCF,应基于步骤(1)建立的残余裂尖张开位移δres与张开应力强度因子Kop的关系:
Kop,CCF=(a1×RCCF+b1)δres,CCF+c1×RCCF+d1
其中,RCCF为复合疲劳载荷下应力比,RCCF=KCCF,min/KCCF,max。
所述步骤(6)中,计算低周载荷卸载过程中裂尖张开位移的变化量Δδ,采用以下方法:
其中,KCCF,max,pre为上一循环最大复合疲劳载荷下应力强度因子;KLCF,min为当前循环低周最小应力强度因子;Kop,CCF,pre为上一循环复合疲劳载荷下张开应力强度因子;C1,n1由步骤(2)确定。
所述步骤(7)中,
首先根据低周载荷下残余裂尖张开位移δLCF,res=δCCF,max-Δδ,和步骤(1)建立的残余裂尖张开位移δres与张开应力强度因子Kop的关系,计算低周载荷下张开应力强度因子Kop,LCF:
Kop,LCF=(a1×RLCF+b1)δLCF,res+c1×RLCF+d1
其中,RLCF为低周载荷下应力比,RLCF=KLCF,min/KLCF,max,
然后采用以下方法计算低周载荷下裂纹增量Δan:
其中,ΔKeff=KLCF,max–Kop,LCF,KLCF,max为低周最大应力强度因子;C0,n0由步骤(1)确定。
本发明与现有技术有益效果在于:
(1)从高低周复合疲劳裂纹扩展机理出发:a.高周载荷的叠加,增加了裂纹扩展损伤,导致裂纹扩展加速;b.高周载荷的存在,残余裂尖张开位移增大,使裂纹闭合水平增加,对裂纹扩展抑制。考虑了裂纹闭合效应对裂纹扩展速率的影响,只有当复合疲劳载荷强度因子大于于张开应力强度因子时,裂纹才会扩展,否则不扩展;
(2)传统复合疲劳裂纹扩展计算方法,将高周与低周割裂开来,线性叠加求得复合疲劳裂纹扩展量,本发明引入了小时间尺度概念,对每个循环内裂纹扩展进行分析计算,考虑了时间历程对裂纹扩展的影响,计算结果更加准确;
(3)基于瞬态动力学分析对高周载荷下榫接结构振动特性进行了分析,可准确模拟高周载荷应力强度因子,从而提高裂纹扩展预测的精度;
(4)现有文献Hu D,Yang Q,Liu H,et al.Crack closure effect and crackgrowth behavior in GH2036superalloy plates under combined high and low cyclefatigue[J].Int J Fatigue.2017,95:90-103从试验角度出发研究了GH2036材料高低周复合疲劳行为,并对其裂纹扩展寿命进行了预测,但其研究成果仅针对实验室情况下标准试件,未考虑涡轮榫接部件的结构特征,未对涡轮榫接结构进行瞬态分析,也未考虑载荷历程对高低周复合疲劳裂纹扩展寿命的影响,其结果针对材料较为单一,缺乏一定的应用性。此外,由于涡轮榫接结构在裂纹扩展过程中,会改变榫接的接触状态,造成应力的重新分布,在计算涡轮榫接结构的裂纹扩展寿命时,必须考虑这方面的影响。本发明补充了该文献的不足,发展了文献提出的方法,首先针对涡轮榫接结构材料开展低周不同应力比下裂纹扩展试验,得到考虑裂纹闭合效应的裂纹扩展模型;在此基础上,对不同裂纹长度下的涡轮榫接结构进行瞬态分析或静力分析,确定榫接结构接触状态对裂纹扩展规律的影响,考虑载荷时间历程的影响,计算涡轮榫接结构高低周裂纹扩展寿命。本发明将文献方法发展成一种可应用于涡轮榫接结构,且适用于多种材料的复合疲劳裂纹扩展寿命预测方法。
附图说明
图1为本发明的涡轮榫接结构高低周复合疲劳裂纹扩展寿命预测方法流程图;
图2为本发明的涡轮榫接结构有限元模型示意图;
图3为本发明的涡轮榫接结构高低周复合疲劳载荷示意图;
图4某涡轮榫接结构高低周复合疲劳裂纹扩展寿命预测结果与试验结果对比。
具体实施方式
下面结合附图,对本发明涡轮叶盘结构强度可靠性设计方法的技术方案做进一步说明。
本发明从高低周复合疲劳裂纹扩展机理出发,即高周载荷的叠加,虽然增加了裂纹扩展损伤,导致裂纹扩展加速,但同时使残余裂尖张开位移增大,使裂纹闭合水平增加,对裂纹扩展抑制,复合疲劳载荷下涡轮榫接结构的裂纹扩展受这两方面移速的相互制约。考虑涡轮榫接结构高周载荷和低周载荷的耦合作用,引入小时间尺度概念,即分别计算每一个循环内裂尖张开位移、应力强度因子和裂纹增量,考虑裂纹闭合效应和结构振动特性,提出一种裂纹闭合模型与瞬态分析相结合的高低周复合疲劳裂纹扩展分析方法,同时保证了预测精度。
如图1所示,本发明涡轮榫接结构高低周复合疲劳裂纹扩展寿命预测方法具体如下:
(1)开展涡轮榫接材料带缺口平板试件低周裂纹扩展试验,根据试验结果计算张开应力强度因子Kop、有效应力强度因子Keff,确定残余裂尖张开位移δres与张开应力强度因子Kop的关系:
Kop=(a1×R+b1)δres+c1×R+d1
其中,a1、b1、c1、d1分别为与材料和温度有关的常数,可通过试验数据或数值模拟结果数据拟合获得;R为载荷的应力比;
并建立考虑裂纹闭合效应的疲劳裂纹扩展寿命模型:
其中,C0,n0为与材料有关常数,可通过试验数据或数值模拟结果拟合获得;ΔKeff为有效应力强度因子幅值,且ΔKeff=Kmax-Kop,Kmax为最大应力强度因子。为后续计算疲劳载荷下的裂纹扩展增量提供基础;
(2)如图2所示,为带裂纹涡轮榫接结构有限元数值模型示意图,对不同裂纹长度下涡轮榫接结构施加高低周复合疲劳载荷,然后进行瞬态分析,根据分析结果拟合确定复合疲劳载荷作用下,高周载荷应力强度因子模型:
△KCCF,max=a2a2+b2a+c2
其中,a2、b2、c2分别为与材料和温度有关的常数,可通过试验数据或数值模拟结果数据拟合获得,a为裂纹长度。
最大裂尖张开位移ΔδCCF,max与高周最大应力强度因子ΔKCCF,max关系模型:
其中,C1,n1为与材料有关常数,可通过数值模拟结果拟合获得。
(3)如图2所示,为带裂纹涡轮榫接结构有限元数值模型示意图,对不同裂纹长度下涡轮榫接结构施加低周疲劳载荷,然后进行瞬态分析,根据分析结果,拟合确定低周载荷应力强度因子模型:
KLCF,max=a3a2+b3a+c3
其中,a3、b3、c3分别为与材料和温度有关的常数,可通过试验数据或数值模拟结果数据拟合获得,a为裂纹长度。
最大裂尖张开位移δLCF,max与低周最大应力强度因子KLCF,max模型:
其中,C2,n2为与材料有关常数,可通过数值模拟结果拟合获得。
(4)如图3所示,为本发明的涡轮榫接结构高低周复合疲劳载荷示意图,假设一个循环起始于卸载,终止于加载到最大载荷,在某一榫齿当前循环裂纹长度下an,判断当前复合疲劳载荷下是否为高周疲劳载荷,若是则执行步骤(3),否则执行步骤(4);
(5)由该榫齿裂纹长度an和步骤(2)确定的高周载荷应力强度因子模型,计算当前循环下的高周复合应力强度因子ΔKCCF,max,计算高低周复合载荷下的裂纹增量Δan:
其中,ΔKCCF,max根据步骤(2)建立的高周应力强度因子模型确定,C0,n0由步骤(1)确定;
判断下一循环状态是否为低周疲劳载荷,若是低周疲劳载荷,计算当前循环复合疲劳载荷下最大裂尖张开位移δCCF,max,应基于步骤(2)和步骤(3)确定的最大裂尖张开位移模型进行计算:
δCCF,max=δLCF,max+△δCCF,max;
计算残余裂尖张开位移δCCF,res:
其中,KCCF,max,KCCF,min为本循环最大和最小复合疲劳载荷应力强度因子;Kop,pre为上一循环张开应力强度因子;C1,n1由步骤(2)确定;
计算张开应力强度因子Kop,CCF:
Kop,CCF=(a1×RCCF+b1)δres,CCF+c1×RCCF+d1
其中RCCF为复合疲劳载荷下应力比,RCCF=KCCF,min/KCCF,max。
然后执行步骤(7);
否则,下一循环状态不是低周疲劳载荷,执行步骤(8);
(6)由该榫齿裂纹长度an和步骤(3)确定的低周应力强度因子模型,计算低循环应力强度因子KLCF及卸载过程中裂尖张开位移的变化量Δδ:
其中,KCCF,max,pre为上一循环最大复合疲劳载荷下应力强度因子;KLCF,min为当前循环下低周最小应力强度因子KLCF,min;Kop,CCF,pre为上一循环复合疲劳载荷下张开应力强度因子;C1,n1由步骤(2)确定。
(7)根据步骤(5)计算的复合疲劳载荷下最大裂尖张开位移δCCF,max和步骤(6)计算的裂尖张开位移变化量Δδ,计算低周载荷下残余裂尖张开位移:
δLCF,res=δCCF,max-Δδ
将该值带入步骤(1)建立的残余裂尖张开位移δres与张开应力强度因子Kop的关系,计算低周载荷下张开应力强度因子Kop,LCF:
Kop,LCF=(a1×RLCF+b1)δLCF,res+c1×RLCF+d1
其中,RLCF为低周载荷下应力比,RLCF=KLCF,min/KLCF,max。
然后计算低周载荷下裂纹增量Δan:
其中,ΔKeff=KLCF,max–Kop,LCF,KLCF,max为低周最大应力强度因子,C0,n0由步骤(1)确定;
(8)更新该榫齿裂纹长度an+1=an+Δan,判断最大应力强度因子是否小于断裂韧度,若小于断裂韧度,再次重复步骤(4)至(7);反之若大于等于断裂韧度,计算完成,得到最终裂纹长度,根据初始和最终裂纹长度即可得出高低周复合疲劳裂纹扩展循环寿命N。
根据本发明所述的预测方法进行了某涡轮榫接结构高低周复合疲劳裂纹扩展寿命预测,其预测结果如图4所示,其中横坐标代表通过试验方法获得的涡轮榫接结构裂纹扩展寿命,纵坐标代表通过本发明获得的涡轮榫接结构预测寿命。通过对比可以发现,本发明建立的预测方法预测精度在2.2倍分散带内。本发明首先针对不同裂纹长度涡轮榫接结构进行瞬态分析,然后考虑裂纹闭合效应和载荷历程影响,发展了一种可应用于涡轮榫接结构,且适用于多种材料的复合疲劳裂纹扩展寿命预测方法,具有较强的应用价值。
提供以上实施例仅仅是为了描述本发明的目的,而并非要限制本发明的范围。本发明的范围由所附权利要求限定。不脱离本发明的精神和原理而做出的各种等同替换和修改,均应涵盖在本发明的范围之内。
Claims (7)
1.一种涡轮榫接结构高低周复合疲劳裂纹扩展寿命预测方法,其特征在于实现步骤如下:
(1)根据不同应力比下低周疲劳裂纹扩展试验,计算涡轮盘榫接材料的张开应力强度因子Kop,有效应力强度因子ΔKeff,并确定残余裂尖张开位移δres与张开应力强度因子Kop的关系,并建立考虑裂纹闭合效应的疲劳裂纹扩展寿命模型,为后续计算疲劳载荷下的裂纹扩展增量提供基础;
(2)对不同裂纹长度下涡轮榫接结构施加高低周复合疲劳载荷,然后进行瞬态分析,根据分析结果拟合确定复合疲劳载荷作用下,高周载荷应力强度因子模型ΔKCCF,max=f(a),和最大裂尖张开位移ΔδCCF,max与高周最大应力强度因子ΔKCCF,max关系模型;
(3)对不同裂纹长度下涡轮榫接结构施加低周疲劳载荷,然后进行瞬态分析,根据分析结果,拟合确定低周载荷应力强度因子模型KLCF,max=f(a),和最大裂尖张开位移δLCF,max与低周最大应力强度因子KLCF,max模型;
(4)设定涡轮榫接结构所受复合疲劳载荷的一个循环起始于卸载,终止于加载到最大载荷,在某一榫齿当前循环裂纹长度an下,判断当前复合疲劳载荷下是否为高周疲劳载荷,若是则执行步骤(5),否则执行步骤(6);
(5)由该榫齿裂纹长度an和步骤(2)确定的高周载荷应力强度因子模型,计算当前循环下的高周复合应力强度因子ΔKCCF,max,计算高低周复合载荷下的裂纹增量Δan,判断下一循环状态是否为低周疲劳载荷,若是低周疲劳载荷,由步骤(2)和(3)确定的最大裂尖张开位移模型,计算当前循环复合疲劳载荷下最大裂尖张开位移δCCF,max,残余裂尖张开位移δCCF,res,张开应力强度因子Kop,CCF,执行步骤(7),否则执行步骤(8);
(6)由该榫齿裂纹长度an和步骤(3)确定的低周应力强度因子模型,计算低循环应力强度因子KLCF及卸载过程中裂尖张开位移的变化量Δδ;
(7)根据步骤(5)计算的复合疲劳载荷下最大裂尖张开位移δCCF,max和步骤(6)计算的裂尖张开位移变化量Δδ,计算低周载荷下残余裂尖张开位移δLCF,res=δCCF,max-Δδ,根据步骤(1)建立的考虑裂纹闭合效应的疲劳裂纹扩展寿命模型,计算低周载荷下的裂纹增量Δan;
(8)更新该榫齿裂纹长度an+1=an+Δan,判断最大应力强度因子是否小于断裂韧度,若小于断裂韧度,再次重复步骤(4)至(7);反之若大于等于断裂韧度,计算完成,得到最终裂纹长度,根据初始和最终裂纹长度即可得出高低周复合疲劳裂纹扩展循环寿命N。
2.根据权利要求1所述的一种涡轮榫接结构高低周复合疲劳裂纹扩展寿命预测方法,其特征在于:所述步骤(1)中,残余裂尖张开位移δres与张开应力强度因子Kop的关系为:
Kop=(a1×R+b1)δres+c1×R+d1
其中,a1、b1、c1、d1分别为与材料和温度有关的常数,可通过试验数据或数值模拟结果数据拟合获得;R为载荷的应力比;
考虑裂纹闭合效应的裂纹扩展寿命模型为:
其中,da/dN为裂纹扩展速率;C0,n0为与材料有关常数,可通过试验数据或数值模拟结果拟合获得;ΔKeff为有效应力强度因子幅值,且ΔKeff=Kmax-Kop,Kmax为最大应力强度因子。
3.根据权利要求1所述的一种涡轮榫接结构高低周复合疲劳裂纹扩展寿命预测方法,其特征在于:所述步骤(2)中,高周载荷应力强度因子ΔKCCF,max=KCCF,max-KLCF,max,其中KCCF,max为榫接结构复合疲劳载荷最大应力强度因子,
拟合得到的高周载荷应力强度因子模型ΔKCCF,max=f(a)形式如下:
△KCCF,max=a2a2+b2a+c2
其中,a2、b2、c2分别为与材料和温度有关的常数,可通过试验数据或数值模拟结果数据拟合获得,a为裂纹长度;
高周载荷下最大裂尖张开位移ΔδCCF,max与高周最大应力强度因子ΔKCCF,max的关系模型为:
其中,C1,n1为与材料有关常数,可通过数值模拟结果拟合获得。
4.根据权利要求1所述的一种涡轮榫接结构高低周复合疲劳裂纹扩展寿命预测方法,其特征在于:所述步骤(3)中,拟合得到的低周载荷应力强度因子模型KLCF,max=f(a)形式如下:
KLCF,max=a3a2+b3a+c3
其中,a3、b3、c3分别为与材料和温度有关的常数,可通过试验数据或数值模拟结果数据拟合获得,a为裂纹长度。
拟合得到的低周载荷下最大裂尖张开位移δLCF,max与低周最大应力强度因子KLCF,max的关系模型为:
其中,C2,n2为与材料有关常数,可通过数值模拟结果拟合获得。
5.根据权利要求1所述的一种涡轮榫接结构高低周复合疲劳裂纹扩展寿命预测方法,其特征在于:所述步骤(5)中,计算高低周复合疲劳载荷下裂纹增量Δan时,采用以下方法:
其中,ΔKCCF,max根据步骤(2)建立的高周应力强度因子模型确定;C0,n0由步骤(1)确定;
计算当前循环复合疲劳载荷下最大裂尖张开位移δCCF,max时,应基于步骤(2)和步骤(3)确定的最大裂尖张开位移模型进行计算,
δCCF,max=δLCF,max+△δCCF,max;
计算当前循环下复合疲劳载荷下残余裂尖张开位移δCCF,res,采用以下方法:
其中,KCCF,max,KCCF,min为本循环最大和最小复合疲劳载荷应力强度因子;Kop,pre为上一循环张开应力强度因子;C1,n1由步骤(2)确定;
计算当前循环下复合疲劳载荷下张开应力强度因子Kop,CCF,应基于步骤(1)建立的残余裂尖张开位移δres与张开应力强度因子Kop的关系:
Kop,CCF=(a1×RCCF+b1)δres,CCF+c1×RCCF+d1
其中,RCCF为复合疲劳载荷下应力比,RCCF=KCCF,min/KCCF,max。
6.根据权利要求1所述的一种涡轮榫接结构高低周复合疲劳裂纹扩展寿命预测方法,其特征在于:所述步骤(6)中,计算低周载荷卸载过程中裂尖张开位移的变化量Δδ,采用以下方法:
其中,KCCF,max,pre为上一循环最大复合疲劳载荷下应力强度因子;KLCF,min为当前循环低周最小应力强度因子;Kop,CCF,pre为上一循环复合疲劳载荷下张开应力强度因子;C1,n1由步骤(2)确定。
7.根据权利要求1所述的一种涡轮榫接结构高低周复合疲劳裂纹扩展寿命预测方法,其特征在于:所述步骤(7)中,
首先根据低周载荷下残余裂尖张开位移δLCF,res=δCCF,max-Δδ,和步骤(1)建立的残余裂尖张开位移δres与张开应力强度因子Kop的关系,计算低周载荷下张开应力强度因子Kop,LCF:
Kop,LCF=(a1×RLCF+b1)δLCF,res+c1×RLCF+d1
其中,RLCF为低周载荷下应力比,RLCF=KLCF,min/KLCF,max,
然后采用以下方法计算低周载荷下裂纹增量Δan:
其中,ΔKeff=KLCF,max–Kop,LCF,KLCF,max为低周最大应力强度因子;C0,n0由步骤(1)确定。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611266906.3A CN106644490B (zh) | 2016-12-31 | 2016-12-31 | 一种涡轮榫接结构高低周复合疲劳裂纹扩展寿命预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611266906.3A CN106644490B (zh) | 2016-12-31 | 2016-12-31 | 一种涡轮榫接结构高低周复合疲劳裂纹扩展寿命预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106644490A true CN106644490A (zh) | 2017-05-10 |
CN106644490B CN106644490B (zh) | 2018-08-03 |
Family
ID=58838366
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611266906.3A Active CN106644490B (zh) | 2016-12-31 | 2016-12-31 | 一种涡轮榫接结构高低周复合疲劳裂纹扩展寿命预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106644490B (zh) |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107449894A (zh) * | 2017-08-23 | 2017-12-08 | 中国人民解放军海军工程大学 | 氮化硅结合碳化硅耐火陶瓷的疲劳寿命预测方法 |
CN108444720A (zh) * | 2018-05-29 | 2018-08-24 | 北京航空航天大学 | 一种涡轮小叶片榫接结构的高低周复合疲劳试验夹具 |
CN109142049A (zh) * | 2018-07-24 | 2019-01-04 | 北京工业大学 | 一种基于小裂纹扩展速率模型的疲劳寿命预测方法 |
CN109344553A (zh) * | 2018-11-29 | 2019-02-15 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种高低周复合疲劳载荷作用下结构细节寿命预测方法 |
CN109374449A (zh) * | 2018-09-25 | 2019-02-22 | 南京航空航天大学 | 一种考虑高、低周疲劳的叶片前后缘裂纹型硬物损伤可用极限确定方法 |
CN109376417A (zh) * | 2018-10-15 | 2019-02-22 | 中国海洋大学 | 一种海洋工程结构低周疲劳裂纹扩展速度与方向的预测方法 |
CN109885920A (zh) * | 2019-02-13 | 2019-06-14 | 西安电子科技大学 | 一种航空发动机涡轮叶片的高低周复合疲劳寿命预测方法 |
CN109918701A (zh) * | 2018-12-29 | 2019-06-21 | 北京航空航天大学 | 一种基于分段权函数的涡轮盘裂纹扩展模拟方法 |
CN109918776A (zh) * | 2019-03-05 | 2019-06-21 | 北京航空航天大学 | 基于双步最小二乘法的机械产品的疲劳裂纹预测方法 |
CN110147618A (zh) * | 2019-05-22 | 2019-08-20 | 电子科技大学 | 基于断裂力学的航空发动机涡轮叶片可靠性评估方法 |
CN110765695A (zh) * | 2019-11-22 | 2020-02-07 | 昆明理工大学 | 一种基于高阶有限元法获取混凝土重力坝裂纹扩展路径的模拟计算方法 |
CN110987676A (zh) * | 2019-12-23 | 2020-04-10 | 北京工业大学 | 随机多轴载荷下考虑裂纹闭合效应的全寿命预测方法 |
CN112906281A (zh) * | 2021-03-15 | 2021-06-04 | 中国航发湖南动力机械研究所 | 一种基于拟蒙特卡洛抽样的涡轮盘裂纹扩展可靠性分析方法 |
CN113514351A (zh) * | 2020-12-23 | 2021-10-19 | 中国特种设备检测研究院 | 一种考虑预应力重分布的疲劳裂纹扩展行为预测方法 |
CN114840943A (zh) * | 2022-05-16 | 2022-08-02 | 北京航空航天大学 | 一种基于裂纹扩展路径和应力强度因子一致的疲劳裂纹扩展模拟件设计方法 |
CN115436033A (zh) * | 2022-07-26 | 2022-12-06 | 北京航空航天大学 | 一种考虑变幅载荷作用的榫接多通道传力结构全寿命模拟方法 |
CN115527635A (zh) * | 2022-09-19 | 2022-12-27 | 南京航空航天大学 | 一种适用于谱载荷下的疲劳裂纹扩展寿命的预测方法 |
CN116776586A (zh) * | 2023-06-15 | 2023-09-19 | 上海发电设备成套设计研究院有限责任公司 | 核电汽轮机的转子应力腐蚀与疲劳长寿命监控方法和装置 |
CN118013666A (zh) * | 2024-04-09 | 2024-05-10 | 中国航发四川燃气涡轮研究院 | 用于气冷涡轮叶片高周疲劳裂纹试验的叶片及设计方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020162400A1 (en) * | 2001-03-05 | 2002-11-07 | Ming Xie | Multiaxial high cycle fatigue test system |
JP2013088262A (ja) * | 2011-10-17 | 2013-05-13 | National Institute For Materials Science | 極低温超音波疲労非破壊試験評価装置及び解析・評価方法 |
CN103592018A (zh) * | 2013-11-14 | 2014-02-19 | 北京航空航天大学 | 一种高低周复合疲劳试验高周振幅测量装置及方法 |
CN104535438A (zh) * | 2014-12-31 | 2015-04-22 | 北京航空航天大学 | 一种试验件高温高低周复合疲劳裂纹扩展试验系统及测量方法 |
CN105404756A (zh) * | 2015-12-17 | 2016-03-16 | 北京航空航天大学 | 一种涡轮叶盘结构寿命可靠性设计方法 |
-
2016
- 2016-12-31 CN CN201611266906.3A patent/CN106644490B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020162400A1 (en) * | 2001-03-05 | 2002-11-07 | Ming Xie | Multiaxial high cycle fatigue test system |
JP2013088262A (ja) * | 2011-10-17 | 2013-05-13 | National Institute For Materials Science | 極低温超音波疲労非破壊試験評価装置及び解析・評価方法 |
CN103592018A (zh) * | 2013-11-14 | 2014-02-19 | 北京航空航天大学 | 一种高低周复合疲劳试验高周振幅测量装置及方法 |
CN104535438A (zh) * | 2014-12-31 | 2015-04-22 | 北京航空航天大学 | 一种试验件高温高低周复合疲劳裂纹扩展试验系统及测量方法 |
CN105404756A (zh) * | 2015-12-17 | 2016-03-16 | 北京航空航天大学 | 一种涡轮叶盘结构寿命可靠性设计方法 |
Non-Patent Citations (1)
Title |
---|
田常海 等: "复合型裂纹断裂和扩展速率实验夹具及裂纹长度测量方法", 《实验力学》 * |
Cited By (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107449894A (zh) * | 2017-08-23 | 2017-12-08 | 中国人民解放军海军工程大学 | 氮化硅结合碳化硅耐火陶瓷的疲劳寿命预测方法 |
CN107449894B (zh) * | 2017-08-23 | 2019-08-20 | 中国人民解放军海军工程大学 | 氮化硅结合碳化硅耐火陶瓷的疲劳寿命预测方法 |
CN108444720A (zh) * | 2018-05-29 | 2018-08-24 | 北京航空航天大学 | 一种涡轮小叶片榫接结构的高低周复合疲劳试验夹具 |
CN109142049A (zh) * | 2018-07-24 | 2019-01-04 | 北京工业大学 | 一种基于小裂纹扩展速率模型的疲劳寿命预测方法 |
CN109374449A (zh) * | 2018-09-25 | 2019-02-22 | 南京航空航天大学 | 一种考虑高、低周疲劳的叶片前后缘裂纹型硬物损伤可用极限确定方法 |
CN109376417A (zh) * | 2018-10-15 | 2019-02-22 | 中国海洋大学 | 一种海洋工程结构低周疲劳裂纹扩展速度与方向的预测方法 |
CN109344553A (zh) * | 2018-11-29 | 2019-02-15 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种高低周复合疲劳载荷作用下结构细节寿命预测方法 |
CN109344553B (zh) * | 2018-11-29 | 2023-04-18 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种高低周复合疲劳载荷作用下结构细节寿命预测方法 |
CN109918701A (zh) * | 2018-12-29 | 2019-06-21 | 北京航空航天大学 | 一种基于分段权函数的涡轮盘裂纹扩展模拟方法 |
CN109885920B (zh) * | 2019-02-13 | 2020-09-08 | 西安电子科技大学 | 一种航空发动机涡轮叶片的高低周复合疲劳寿命预测方法 |
CN109885920A (zh) * | 2019-02-13 | 2019-06-14 | 西安电子科技大学 | 一种航空发动机涡轮叶片的高低周复合疲劳寿命预测方法 |
CN109918776B (zh) * | 2019-03-05 | 2020-09-11 | 北京航空航天大学 | 基于双步最小二乘法的机械产品的疲劳裂纹预测方法 |
CN109918776A (zh) * | 2019-03-05 | 2019-06-21 | 北京航空航天大学 | 基于双步最小二乘法的机械产品的疲劳裂纹预测方法 |
CN110147618A (zh) * | 2019-05-22 | 2019-08-20 | 电子科技大学 | 基于断裂力学的航空发动机涡轮叶片可靠性评估方法 |
CN110765695A (zh) * | 2019-11-22 | 2020-02-07 | 昆明理工大学 | 一种基于高阶有限元法获取混凝土重力坝裂纹扩展路径的模拟计算方法 |
CN110765695B (zh) * | 2019-11-22 | 2022-11-15 | 昆明理工大学 | 一种基于高阶有限元法获取混凝土重力坝裂纹扩展路径的模拟计算方法 |
CN110987676A (zh) * | 2019-12-23 | 2020-04-10 | 北京工业大学 | 随机多轴载荷下考虑裂纹闭合效应的全寿命预测方法 |
CN113514351A (zh) * | 2020-12-23 | 2021-10-19 | 中国特种设备检测研究院 | 一种考虑预应力重分布的疲劳裂纹扩展行为预测方法 |
CN113514351B (zh) * | 2020-12-23 | 2024-04-19 | 中国特种设备检测研究院 | 一种考虑预应力重分布的疲劳裂纹扩展行为预测方法 |
CN112906281A (zh) * | 2021-03-15 | 2021-06-04 | 中国航发湖南动力机械研究所 | 一种基于拟蒙特卡洛抽样的涡轮盘裂纹扩展可靠性分析方法 |
CN114840943A (zh) * | 2022-05-16 | 2022-08-02 | 北京航空航天大学 | 一种基于裂纹扩展路径和应力强度因子一致的疲劳裂纹扩展模拟件设计方法 |
CN114840943B (zh) * | 2022-05-16 | 2022-11-25 | 北京航空航天大学 | 一种基于裂纹扩展路径和应力强度因子一致的疲劳裂纹扩展模拟件设计方法 |
CN115436033A (zh) * | 2022-07-26 | 2022-12-06 | 北京航空航天大学 | 一种考虑变幅载荷作用的榫接多通道传力结构全寿命模拟方法 |
CN115436033B (zh) * | 2022-07-26 | 2024-04-30 | 北京航空航天大学 | 一种考虑变幅载荷作用的榫接多通道传力结构全寿命模拟方法 |
CN115527635B (zh) * | 2022-09-19 | 2023-11-10 | 南京航空航天大学 | 一种适用于谱载荷下的疲劳裂纹扩展寿命的预测方法 |
CN115527635A (zh) * | 2022-09-19 | 2022-12-27 | 南京航空航天大学 | 一种适用于谱载荷下的疲劳裂纹扩展寿命的预测方法 |
CN116776586A (zh) * | 2023-06-15 | 2023-09-19 | 上海发电设备成套设计研究院有限责任公司 | 核电汽轮机的转子应力腐蚀与疲劳长寿命监控方法和装置 |
CN118013666A (zh) * | 2024-04-09 | 2024-05-10 | 中国航发四川燃气涡轮研究院 | 用于气冷涡轮叶片高周疲劳裂纹试验的叶片及设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106644490B (zh) | 2018-08-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106644490A (zh) | 一种涡轮榫接结构高低周复合疲劳裂纹扩展寿命预测方法 | |
US10401270B2 (en) | Method of design and manufacturing concrete structures based on the verification of concrete fatigue strength by test | |
CN106840877B (zh) | 一种基于应力的多轴小裂纹全寿命预测方法 | |
CN107145641B (zh) | 叶片振动疲劳概率寿命预测方法 | |
Citarella et al. | FEM-DBEM approach for crack propagation in a low pressure aeroengine turbine vane segment | |
CN109918789B (zh) | 一种多轴变幅加载下基于短裂纹扩展的全寿命预测方法 | |
CN106644784A (zh) | 一种考虑多部位及多失效模式的涡轮盘损伤容限评估方法 | |
CN115640666B (zh) | 一种基于损伤等效的航空发动机加速任务试车谱编制方法 | |
CN109255202A (zh) | 一种用于机械构件疲劳裂纹萌生寿命的预估方法 | |
Schijve | Fatigue crack growth under variable-amplitude loading | |
CN110083984A (zh) | 一种涡轮叶片区间损伤容限分析方法 | |
CN111507038B (zh) | 一种陶瓷基复合材料结构疲劳寿命预测方法 | |
CN107122531A (zh) | 一种基于加速寿命试验的高速叶轮寿命评估方法 | |
CN107449894A (zh) | 氮化硅结合碳化硅耐火陶瓷的疲劳寿命预测方法 | |
Xiannian et al. | Fatigue crack propagation analysis in an aero-engine turbine disc using computational methods and spin test | |
Zhang et al. | A multiaxial probabilistic fatigue life prediction method for nickel‐based single crystal turbine blade considering mean stress correction | |
Qiu et al. | A dual-threshold modelling approach for fatigue life prediction under combined high and low cycle fatigue | |
Měšťánek | Low cycle fatigue analysis of a last stage steam turbine blade | |
CN108334666A (zh) | 一种车轴抗疲劳断裂分析方法 | |
CN116050202B (zh) | 适应多种平均应力表现的Chaboche叶片疲劳寿命预测模型及方法 | |
Qi et al. | Fatigue Reliability Analysis System for Key Components of Aero‐Engine | |
Yu et al. | A LCF life assessment method for steam turbine long blade based on Elastoplastic analysis and local strain approach | |
Colombo et al. | Simple modelling of the constitutive behaviour of a 1% CrMoV rotor steel in service-like thermo-mechanical fatigue tests | |
Inozemtsev et al. | Life validation strategy | |
CN106442298A (zh) | 一种自然暴晒试验中金属构件应力腐蚀模拟装置及方法 |
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 |