一种高温合金蠕变寿命预测方法
【专利说明】-种高溫合金蠕变寿命预测方法
[0001]
技术领域
[0002] 本发明属于材料科学与工程技术应用，具体设及高溫合金蠕变寿命预测方法。
[0003]
【背景技术】
[0004] 自20世纪30年代开始，世界各国开始高溫合金的研究，并在接下来的半个多世纪 中迅猛发展。高溫合金W其优异的高溫强度，良好的抗氧化和抗热腐蚀性能，良好的疲劳性 能、断裂初性等综合性能，被广泛应用于制造航空、舰艇和工业用燃气轮机的满轮叶片、导 向叶片、满轮盘、高压压气机盘和燃烧室等高溫部件，已成为军民用燃气满轮发动机热端部 件不可替代的关键材料。
[0005] 高溫合金在服役过程中通常会因为蠕变变形导致失效，从而造成巨大的损失。因 此若能够准确地预测高溫合金蠕变时间，尤其是某一蠕变应变条件对应的蠕变中断时间， 则可对高溫合金服役时间进行有效的管理，在避免造成合金蠕变失效的情况下最大程度地 对其进行利用。
[0006] 从20世纪50年代开始，各种蠕变寿命预测模型被提出。1982年R.M. Evans等人提 出通过对蠕变曲线进行描述，从而预测蠕变寿命的方法，称为"Θ投影法"。Θ投影法认为蠕 变过程可分为材料软化过程和材料硬化过程，且将材料软化过程和材料硬化过程均W指数 方程的形式表示，二者叠加可得Θ投影法方程。由于其不仅能对不同条件下的蠕变断裂时间 进行预测，同时还可预测某一蠕变应变所对应的蠕变中断时间，Θ投影法在提出后被广泛应 用于各种类型的合金。但是，Θ投影法的提出基于理想蠕变条件，即溫度和应力均为恒定;而 在真实蠕变过程中，蠕变工程应力恒定，随着应变量的增加，蠕变样品的横截面逐渐缩小， 从而导致蠕变真应力逐渐增加，二者互相矛盾。由于蠕变过程中形状变化导致的应力变化 不可避免，且会导致蠕变寿命低于恒定应力下的蠕变寿命，从而使得利用模拟蠕变中断曲 线所得的Θ投影法方程预测同一条件下更高应变对应的蠕变时间时结果偏大;更有甚者，对 于蠕变应变较大的蠕变曲线，其曲线不符合Θ投影法的规律。因此需要针对此现象对Θ投影 法进行修正，W求得到更高精度的蠕变寿命预测模型。
【发明内容】
[0007] 本发明的目的是提出一种可准确描述蠕变曲线的蠕变寿命预测方法，该方法通过 引入蠕变实验中由于蠕变样品形状变化导致的应力变化运一因素，对Θ投影法进行修正，减 小预测的蠕变曲线与实验所得的蠕变曲线之间的差别，该方法简单可靠，适用于不同类型 的局溫合金。
[0008] 一种高溫合金的蠕变寿命预测方法，其核屯、要点在于：引入实际蠕变过程中，蠕变 真应力随应变增加而增加，从而导致蠕变寿命降低运一因素，即考虑蠕变过程中形状变化 运一损伤参量。本发明的实现步骤如下： (1)获取至少5条不同溫度和应力条件下合金的高溫蠕变曲线，每条蠕变曲线包括溫 度、应力、蠕变应变和蠕变时间。
[0009] (2)将蠕变曲线第一阶段至第Ξ阶段初期的数据按公式 ￡=目1*(1-0邱（目2*(1 +目日*〇*1：)) +目3*(6邱（目4*(1 +目日*6)*1：)-1)进行拟合，求出公式中的 参数目1、目2、目3、目4、目5,得到修正后的Θ投影法方程。
[0010] (3)利用(2)中结果可通过延长Θ投影法曲线，获得均匀变形阶段的蠕变曲线。
[ocm] (4)将不同溫度和应力条件下的蠕变曲线拟合得到的目1、02、目3、目4、目5按公式1〇邑目1= a+b*〇+c*T+d*〇*T (i=l, 2,3,4)和公式log目日=e+f*〇分别计算参数a、b、c、d、e、f。
[0012] (5)利用（2)和(4)所得结果可预测不同溫度和应力下的蠕变曲线，从而获得对应 溫度和应力下的蠕变断裂时间或某一应变所对应的蠕变中断时间。
[001；3]所述公式￡；=目1*(1-6邱（目2*(1 +目日*6)*1：))+目3*(6邱（目4*(1 +目日*6)*1：)-1)中，6表示蠕 变应变，t表示蠕变时间。
[0014] 所述公式log目i=a+b*〇+c*T+d*〇*T (i=l ,2,3,4)和公式1〇旨目日=6+巧〇中，1'表示蠕变 溫度，0表示蠕变工程应力，a、b、C、d、e、f为与材料相关的常数。
[001引所述参数01、02、03、04、05是利用数学分析软件，按最小二乘法回归，将蠕变曲线时 间-应变数值输入，按公式ε=目i*(l-exp (目2*(1+目日*e)*t))+目3*(6邱（目4*(1+目日*6)*1：)-1)求得 待定系数目1、目2、目3、目4、白已。
[0016] 所述参数a、b、c、d、e、f是利用数学分析软件，按最小二乘法回归，将蠕变溫度、蠕 变应力及对应的目1、目2、目3、目4、目日输入，按公式log目i=a+b*〇+c*T+d*〇*T (i=l，2，3，4)和公式 log目日=e+f*〇求得待定系数a、b、c、d、e、f。
[0017] 本发明的优点在于:能够准确地对不同类型的高溫合金蠕变曲线进行描述，能够 提高高溫合金蠕变断裂时间和某一应变量对应的蠕变中断时间的预测精度。该方法简单可 靠，利用常规蠕变试验机获得蠕变数据。利用此方法可对高溫合金的服役时间进行有效估 计，减小危害，降低成本，适用于不同类型的高溫合金。
【附图说明】
[0018] 图1为K465合金不同蠕变条件下的蠕变时间-应变曲线。
[0019] 图2为K465合金蠕变实验曲线与修正后的Θ投影法预测曲线对比图。
[0020] 图3为K465合金蠕变实验曲线与经典Θ投影法预测曲线对比图（a)900°C/300MPa， (b)900°C/320MPa, (c)950°C/300MPa, (d)975°C/200MPa, (e) 1000°C/137MPa〇
[0021] 图4为修正后的Θ投影法预测不同条件下的K465合金蠕变曲线。
【具体实施方式】
[0022] W下实例将对本发明予W进一步的说明，W便本领域人员更好地理解本发明的优 点和特征。
[0023] 首先，WK465合金为研究对象，利用高溫蠕变试验机，获得其在五个不同溫度和应 力条件下的蠕变曲线，实验条件分别为1000°C/137MPa、975°C/200MPa、950°C/300MPa、900 °C/320MPa、900°C/300MPa，实验时每隔30min取一个数据点，实验结果见图1。
[0024] 其次，分别利用公式 e=0i*(l-e邱（目 2*(1+目5*e)*t))+目3*(exp (目4*(1+目5*e)*t)-l)及 经典e投影法方程对蠕变曲线ε<1.3%的蠕变数据进行拟合，获得待定系数θι、θ2、θ3、θ4、θ5 (见表1)，得到修正后的e投影法方程。具体拟合方法为利用化iginlab软件的非线性拟合 工具，将方程写入并选择合适的初始值，选择需要拟合的数据进行拟合。
[0025] 图2为利用修正后的Θ投影法方程作出的完整蠕变曲线与实验所得蠕变曲线的对 比。图3为经典Θ投影法方程作出的完整蠕变曲线与实验所得蠕变曲线的对比。由图2和图3 可W看出，修正后的e投影法大大提高了预测精度。利用不同蠕变条件下的系数9?、92、θ3、 白4、目日及公式log目i=a+b*〇+c*T+d*〇*T (i=l ,2,3,4)和公式1〇邑目日=6+巧〇拟合得出待定系数 a、b、c、d、e、f，见表2。具体拟合过程由化iginlab软件的线性拟合工具实现。
[0026] 利用上述结果，通过计算不同溫度和应力下的系数目1、目2、目3、目4、目5，可得到对应的 修正后的Θ投影法方程，从而达到预测不同条件蠕变曲线的目的，见图4。
[0027] 表1修正后的e投影法拟合K465合金不同蠕变条件下的系数
【主权项】
1. 一种高温合金的蠕变寿命预测方法，其特征在于，具体包括如下步骤： 1) 获取至少5条不同温度和应力条件下合金的高温蠕变曲线，每条蠕变曲线包括温度、 应力、蠕变应变和蠕变时间； 2) 将蠕变曲线第一阶段至第三阶段初期的数据按公式 ε=θ1*(1-〇χρ(θ2*(1+θ5*ε)*1：))+θ3*(θχρ(θ4*(1+θ 5*ε)*1：)-1)进行拟合，求出公式中的待 定系数01、02、03、04、05，得到修正后的0投影法方程； 3) 利用2)中结果通过延长Θ投影法曲线，获得均匀变形阶段的蠕变曲线； 4) 将不同温度和应力条件下的懦变曲线拟合得到的θι、θ2、θ3、θ4、θ 5按公式l〇g0i=a+b*〇 +c*T+d*〇*T，i=l ,2,3,4和公式log05=e+f*〇分别计算参数a、b、c、d、e、f; 5) 利用2)和4)所得结果预测不同温度和应力下的蠕变曲线，从而获得对应温度和应力 下的蠕变断裂时间或蠕变中断时间。2. 按照权利要求1所述的一种高温合金蠕变寿命预测方法，其特征在于 6叉口（02*(1 + 95#)*1：)) + 03*(6叉口（04*(1 + 9糾）*1：)-1)中，￡表示懦变应变，1：表示懦变时间。3. 按照权利要求1所述的一种高温合金蠕变寿命预测方法，其特征在于:公式logepa- b*〇+c*T+d*〇*T (i=l，2，3，4)和公式log05=e+f*o中，T表示蠕变温度，σ表示蠕变工程应力， a、b、c、d、e、f为与材料相关的常数。4. 按照权利要求1所述的一种高温合金蠕变寿命预测方法，其特征在于:按最小二乘法 回归，将懦变曲线时间-应变数值输入，按公式e=9i*(l-exp(02*(l+9 5*e)*t))+03*(exp(04* (1+θδ*ε )*t )-1)求得待定系数θι、Θ2、Θ3、Θ4、Θ5。5. 按照权利要求1所述的一种高温合金蠕变寿命预测方法，其特征在于:按最小二乘法 回归，将蠕变温度、蠕变应力及对应的θ!、θ 2、θ3、θ4、θ5输入，按公式logeiza+Wo+MT+cH^T (i=l ,2,3,4)和公式log05=e+f*〇求得待定系数 a、b、c、d、e、f。
【专利摘要】本发明提供了一种高温合金的蠕变寿命预测方法，该方法是在θ投影法的基础上，引入蠕变过程中材料形状变化导致真应力变化这一因素，得到修正θ投影法，具体形式为：ε=θ1*(1-exp(θ2*(1+θ5*ε)*t))+θ3*(exp(θ4*(1+θ5*ε)*t)-1)。该方法可用于利用蠕变曲线的第一阶段至第三阶段初期数据预测整个均匀变形阶段的蠕变曲线，有利于缩短获得蠕变数据的时间。同时还可用于预测其他蠕变条件下的蠕变曲线，包括预测蠕变某一蠕变条件下的蠕变断裂时间和某一蠕变应变对应的蠕变中断时间。该方法简单可靠，适用于不同类型的高温合金，适合工程应用，在高温合金蠕变寿命管理方面具有广阔前景。
【IPC分类】G01N3/18
【公开号】CN105628511
【申请号】CN201510997412
【发明人】冯强, 付超, 袁晓飞, 郭小童
【申请人】北京科技大学
【公开日】2016年6月1日
【申请日】2015年12月25日