一种低温疲劳性能表征与寿命估算的方法，该方法有三大步骤：步骤一、考虑应力比的影响，在三参数幂函数表达式基础上，采用Goodman等寿命曲线，推导出表征疲劳性能的S‑N‑R曲面，并借助线性回归理论进行数据处理；步骤二、考虑谱载下载荷间的交互作用，在Willenborg/Chang模型和裂尖塑性区理论基础上，提出修正的谱载疲劳寿命估算模型，给出疲劳损伤增量的表示方法；步骤三、采用累积损伤理论计算材料的低温谱载疲劳寿命。本发明专利技术简单实用，仅需要低温环境下材料的恒载疲劳性能曲面和实测飞行载荷谱，便可构建低温疲劳性能表征模型，并估算谱载疲劳寿命，具有重要学术意义和工程应用价值。
【技术实现步骤摘要】
一种低温疲劳性能表征与寿命估算的方法
本专利技术提供一种低温疲劳性能表征与寿命估算的方法，属于金属结构疲劳可靠性
技术介绍
实际工程中，材料常会受到交变载荷的作用而产生损伤，出现疲劳裂纹；并在一定的循环应力作用后发生断裂失效，从而对结构的安全性造成威胁。低温是材料使用过程中不可避免的环境因素，例如：由于飞行高度和气候的影响，航空器在某些区域的工作温度达到-60℃；液氮储存设备和低温超导材料的工作温度更是低于-250℃；低温环境下，材料常会发生疲劳韧脆转变和循环硬化，其宏观性能也通常产生变化，疲劳行为也有所不同，因此，研究材料的低温疲劳行为有重要的实际意义。目前，尚缺乏更为精确而实用的低温疲劳性能表征模型和寿命估算方法，为此，专利技术了一种简单实用的低温疲劳性能表征与寿命估算的方法，该方法仅需要低温环境下材料的恒载疲劳性能曲面和实测飞行载荷谱，便可构建低温疲劳性能表征模型，并估算材料的谱载疲劳寿命，本专利技术具有重要学术意义和工程应用价值。
技术实现思路
1、目的：本专利技术目的是提供了一种低温疲劳性能表征与寿命估算的方法，该方法具有所需计算参数少、计算简便和精度较高等特点，对于低温环境下材料的疲劳性能和谱载寿命评估有重要价值。2、技术方案：一种低温疲劳性能表征与寿命估算的方法，该方法具体步骤如下，步骤一、低温疲劳性能表征模型图1为疲劳中心缺口试件的加载示意图，按照图1的加载形式和标准《金属材料恒幅轴向疲劳试验方法》(HB5287-96)，在低温环境下进行疲劳试验。在单一应力比R下，采用成组试验法和升降试验法，选取不同的最大疲劳应力Smax加载，得到不同寿命范围的疲劳失效循环数N。记录每组疲劳试验结果Smax和N，并采用三参数幂函数经验公式表征单一应力比下的疲劳性能：(Smax-S0)m·N＝C(1)式中，S0是拟合得到的疲劳极限，C和m为材料常数。值得指出的是，三参数幂函数表达式(1)仅能表征单一应力比下的疲劳行为，而实测载荷历程中包含大量不同应力比的载荷循环。因此，有必要借助等寿命曲线经验公式对三参数幂函数表达式(1)进行修正，以表征不同应力比加载下的疲劳性能，表征应力比效应的Goodman等寿命曲线经验公式为：式中，Sa和Sm分别为疲劳应力幅值和应力均值，S-1代表对称循环加载下的疲劳极限，σb为试验测得的材料拉伸强度极限。根据应力比的定义可知：式中Smin是循环加载的最小疲劳应力。对式(3)做恒等变换可以得到：将式(4)和式(5)代入Goodman方程(2)，可以得到疲劳极限S-1的表示方法为：当疲劳试验在对称循环加载(R＝-1)下进行时，最大疲劳应力Smax即为疲劳极限S-1，从而，三参数幂函数式(1)可表示为：(S-1-S0)m·N＝C(7)将式(6)代入式(7)，可以得到表征不同应力比和疲劳极限的S-N-R曲面模型：在S-N-R曲面模型表达式(8)中，C、m和S0为待定参数，需要在试验数据[Smax,R,N]的基础上，借助线性回归方法拟合得到。对式(8)两边取对数，得到：Y1＝a1+b1X1(9)式中Y1＝lgN，a1＝lgC，b1＝-m，根据线性回归理论，待定参数a1和b1的确定方法如下：其中从式(12)至式(15)可以看出，中间变量L11和L12均与待定参数S0有关，也就是说，它们均是S0的函数。因此，待定参数a1和b1也同样是S0的函数，从而，首先需要确定参数S0的取值。根据残差平方和(RSS)理论，可以得到Q(S0)的数值解：具体的求解过程为：(i)确定S0的取值范围为S0∈[0,min(S1,S2,…,Sn)]，其中Si(i＝1,2,…,n)根据试验数据[(Smax)i,Ri]进行确定，(ii)给定S0的初始值(如S0＝0)和计算步长，由式(10)至(15)可以得到参数a1和b1的值，再由式(16)计算函数Q(S0)。如此进行迭代差值计算，可绘制变量S0和函数Q(S0)的关系曲线；(iii)根据S0-Q(S0)关系曲线，确定Q(S0)的最小值及其对应的S0的解。在此基础上，由式(10)和式(11)以及解得的S0值，可以得到参数C和m的拟合结果为：从而，根据式(9)至式(18)并结合测得的试验数据，由曲面模型式(8)可以拟合低温环境下材料的疲劳性能S-N-R曲面(如图2所示)，曲面拟合结果能够有效地表征材料的疲劳性能，并且更直观地反映了低温环境对疲劳行为的影响。步骤二、谱载疲劳寿命估算模型谱载试验采用实测载荷谱加载，图3示出了实测载荷系数谱，再乘以应力水平即为试验加载的实测载荷谱。评估材料的谱载疲劳剩余寿命的线性累积损伤Miner理论为：式中，D是载荷谱中单一应力循环[Smax,R]造成的疲劳损伤，n是单一应力循环在载荷谱中出现的次数，N是疲劳失效循环数，dD/dn是疲劳损伤速率(每载荷循环的疲劳损伤)，ΔDi是随机载荷谱中第i个应力循环造成的疲劳损伤增量，k是一个随机载荷谱块中的应力循环数，T是预测的材料谱载疲劳寿命。但Miner理论需要借助雨流技术处理分离出实测载荷谱中所用不同的应力循环，而忽略了载荷历程中的顺序效应和载荷间的交互作用。实际上，载荷交互作用对疲劳行为有显著影响，包括高载后残余压应力引起的高载迟滞效应，以及低载后残余拉应力引起的抵消迟滞效应等。因此，需要考虑载荷交互作用的影响，构建新的评估谱载疲劳寿命的计算方法。Willenborg/Chang模型考虑了载荷间的交互作用，在塑性区理论基础上，引入有效应力比和有效应力强度因子以表征谱载加载下的裂纹扩展速率。在Willenborg/Chang模型和连续损伤力学理论的基础上，将表征裂纹尖端塑性修正的有效应力比Reff和有效最大疲劳应力Smax,eff引入式(8)和式(19)，考虑载荷间的交互作用表征谱载下的疲劳损伤速率dD/dn，其表示方法为：其中Reff和Smax,eff可由下式确定式中，(Smax)OL是过载应力循环的谱载最大疲劳应力，r是超载截止比，ΔD′是过载后的疲劳损伤增量，zOL是过载迟滞区尺寸参数。对式(21)进行积分变换，可以得到谱载加载下对于第i个循环的疲劳损伤增量：步骤三、谱载疲劳寿命累积损伤算法借助式(25)和式(20)，材料的谱载疲劳剩余寿命可以通过损伤累积计算得到，具体计算方法为：(i)对于随机载荷谱中的第一个应力循环，可由式(22)至(24)得到对应的有效最大疲劳应力(Smax,eff)1和有效应力比(Reff)1，代入式(25),得到第一个循环造成的疲劳损伤增量ΔD1和当前损伤值D1；(ii)同样的，在D1基础上，计算第二个应力循环的有效最大疲劳应力(Smax,eff)2和有效应力比(Reff)2，得到第二个循环造成的疲劳损伤增量ΔD2和当前损伤值D2；(iii)通过如此循环接循环的累积，计算载荷谱中每一后续应力循环造成的疲劳损伤增量直至载荷谱结束，此时，对应的累积损伤即为一个随机载荷谱块所造成的疲劳损伤；(iv)当累积损伤值达到或超过式(20)允许的损伤容限时，疲劳损伤累积计算停止，此时对应的最终失效循环数即为材料的低温谱载疲劳寿命。3、优点及功效：本专利技术提供了一种低温疲劳性能表征与寿命估算的方法，其特点是简单实用，考虑应力比的影响，在三参数幂函数表达式基础上，采用Goodman等寿命本文档来自技高网...
【技术保护点】
一种低温疲劳性能表征与寿命估算的方法，具有所需计算参数少、计算简便和精度较高等特点，该方法具体步骤如下：步骤一、低温疲劳性能表征模型图1为疲劳中心缺口试件的加载示意图，按照图1的加载形式和标准《金属材料恒幅轴向疲劳试验方法》(HB5287‑96)，在低温环境下进行疲劳试验。在单一应力比R下，采用成组试验法和升降试验法，选取不同的最大疲劳应力Smax加载，得到不同寿命范围的疲劳失效循环数N。记录每组疲劳试验结果Smax和N，并采用三参数幂函数经验公式表征单一应力比下的疲劳性能：(Smax‑S0)m·N＝C (1)式中，S0是拟合得到的疲劳极限，C和m为材料常数。值得指出的是，三参数幂函数表达式(1)仅能表征单一应力比下的疲劳行为，而实测载荷历程中包含大量不同应力比的载荷循环。因此，有必要借助等寿命曲线经验公式对三参数幂函数表达式(1)进行修正，以表征不同应力比加载下的疲劳性能，表征应力比效应的Goodman等寿命曲线经验公式为：式中，Sa和Sm分别为疲劳应力幅值和应力均值，S‑1代表对称循环加载下的疲劳极限，σb为试验测得的材料拉伸强度极限。根据应力比的定义可知：式中Smin是循环加载的最小疲劳应力。对式(3)做恒等变换可以得到：将式(4)和式(5)代入Goodman方程(2)，可以得到疲劳极限S‑1的表示方法为：当疲劳试验在对称循环加载(R＝‑1)下进行时，最大疲劳应力Smax即为疲劳极限S‑1，从而，三参数幂函数式(1)可表示为：(S‑1‑S0)m·N＝C (7)将式(6)代入式(7)，可以得到表征不同应力比和疲劳极限的S‑N‑R曲面模型：在S‑N‑R曲面模型表达式(8)中，C、m和S0为待定参数，需要在试验数据[Smax,R,N]的基础上，借助线性回归方法拟合得到。对式(8)两边取对数，得到：Y1＝a1+b1X1 (9) 式中Y1＝lgN，a1＝lgC，b1＝‑m，根据线性回归理论，待定参数a1和b1的确定方法如下：其中从式(12)至式(15)可以看出，中间变量L11和L12均与待定参数S0有关，也就是说，它们均是S0的函数。因此，待定参数a1和b1也同样是S0的函数，从而，首先需要确定参数S0的取值。根据残差平方和(RSS)理论，可以得到Q(S0)的数值解：具体的求解过程为：(i)确定S0的取值范围S0∈[0，min(S1，S2，…,Sn),，其中(ii)给定S0的初始值(如S0＝0)和计算步长，由式(10)至(15)可以得到参数a1和b1的值，再由式(16)计算函数Q(S0)。如此进行迭代差值计算，可绘制变量S0和函数Q(S0)的关系曲线；(iii)根据S0‑Q(S0)关系曲线，确定Q(S0)的最小值及其对应的S0的解。在此基础上，由式(10)和式(11)以及解得的S0值，可以得到参数C和m的拟合结果为：从而，根据式(9)至式(18)并结合测得的试验数据，由曲面模型式(8)可以拟合低温环境下材料的疲劳性能S‑N‑R曲面(如图2所示)，曲面拟合结果能够有效地表征材料的疲劳性能，并且更直观地反映了低温环境对疲劳行为的影响。步骤二、谱载疲劳寿命估算模型谱载试验采用实测载荷谱加载，图3示出了实测载荷系数谱，再乘以应力水平即为试验加载的实测载荷谱。评估材料的谱载疲劳剩余寿命的线性累积损伤Miner理论为：式中，D是载荷谱中单一应力循环[Smax,R]造成的疲劳损伤，n是单一应力循环在载荷谱中出现的次数，N是疲劳失效循环数，dD/dn是疲劳损伤速率(每载荷循环的疲劳损伤)，ΔDi 是随机载荷谱中第i个应力循环造成的疲劳损伤增量，k是一个随机载荷谱块中的应力循环数，T是预测的材料谱载疲劳寿命。但Miner理论需要借助雨流技术处理分离出实测载荷谱中所用不同的应力循环，而忽略了载荷历程中的顺序效应和载荷间的交互作用。实际上，载荷交互作用对疲劳行为有显著影响，包括高载后残余压应力引起的高载迟滞效应，以及低载后残余拉应力引起的抵消迟滞效应等。因此，需要考虑载荷交互作用的影响，构建新的评估谱载疲劳寿命的计算方法。Willenborg/Chang模型考虑了载荷间的交互作用，在塑性区理论基础上，引入有效应力比和有效应力强度因子以表征谱载加载下的裂纹扩展速率。在Willenborg/Chang模型和连续损伤力学理论的基础上，将表征裂纹尖端塑性修正的有效应力比Reff和有效最大疲劳应力Smax,eff引入式(8)和式(19)，考虑载荷间的交互作用表征谱载下的疲劳损伤速率dD/dn，其表示方法为：其中Reff和Smax,eff可由下式确定式中，(Smax)OL是过载应力循环的谱载最大疲劳应力，r是超载截止比，ΔD′是过载后的疲劳损伤增量，zOL是过载迟滞区尺寸参数。对式(21)进行积分变换，可以得到谱载加载下对于第i个循环的疲劳损伤增量：步...
【技术特征摘要】
1.一种低温疲劳性能表征与寿命估算的方法，该方法具体步骤如下：步骤一、低温疲劳性能表征模型按照疲劳中心缺口试件的加载形式和标准《金属材料轴向加载疲劳试验方法》(HB5287-96)，在低温环境下进行疲劳试验；在单一应力比R下，采用成组试验法和升降试验法，选取不同的最大疲劳应力Smax加载，得到不同寿命范围的疲劳失效循环数N；记录每组疲劳试验结果Smax和N，并采用三参数幂函数经验公式表征单一应力比下的疲劳性能：(Smax-S0)m·N＝C(1)式中，S0是拟合得到的疲劳极限，C和m为材料常数；借助等寿命曲线经验公式对三参数幂函数表达式(1)进行修正，以表征不同应力比加载下的疲劳性能，表征应力比效应的Goodman等寿命曲线经验公式为：式中，Sa和Sm分别为疲劳应力幅值和应力均值，S-1代表对称循环加载下的疲劳极限，σb为试验测得的材料拉伸强度极限；根据应力比的定义可知：式中Smin是循环加载的最小疲劳应力；对式(3)做恒等变换可以得到：将式(4)和式(5)代入Goodman方程(2)，可以得到疲劳极限S-1的表示方法为：当疲劳试验在对称循环加载，即R＝-1下进行时，最大疲劳应力Smax即为疲劳极限S-1，从而，三参数幂函数式(1)可表示为：(S-1-S0)m·N＝C(7)将式(6)代入式(7)，可以得到表征不同应力比和疲劳极限的S-N-R曲面模型：在S-N-R曲面模型表达式(8)中，C、m和S0为待定参数，需要在试验数据[Smax,R,N]的基础上，借助线性回归方法拟合得到；对式(8)两边取对数，得到：Y1＝a1+b1X1(9)式中Y1＝lgN，a1＝lgC，b1＝-m，根据线性回归理论，待定参数a1和b1的确定方法如下：其中根据残差平方和(RSS)理论，可以得到Q(S0)的数值解：具体的求解过程为：(i)确定S0的取值范围为S0∈[0,min(S1,S2,…,Sn)]，其中Si(i＝1,2,…,n)根据试验数据[(Smax)i,Ri]进行确定，(ii)给定S0的初始值和计算步长，由式(10)至(15)可以得到参数a1和b1的值，再由式(16)计算函数Q(S0)；如此进行迭代差值计算，可绘制变量S0和函数Q(S0)的关系曲线；(iii)根据S0-Q(S0)关系曲线，确定...
【专利技术属性】
技术研发人员：熊峻江，刘牧东，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：北京;11
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。