一种基于分区的涡轮盘疲劳寿命及失效风险评估方法

摘要

本发明涉及一种基于分区的涡轮盘疲劳寿命及失效风险评估方法，步骤：(1)获取涡轮盘几何、载荷、材料随机变量概率分布。(2)进行涡轮盘有限元分析，获取应力和应变分布，并基于温度分布和应力分布形成初始分区。(3)确定每一个区域失效相应的疲劳寿命模型，代入应力分布和应变分布计算每一个区域的失效风险。(4)选定某一疲劳寿命，基于步骤(3)中各区域的失效风险，进行分区优化得到优化后的分区方案。(5)基于优化后分区方案，重复步骤(3)得到每一个区域输入变量和疲劳寿命的第一级响应面；同时建立每一个分区失效寿命和涡轮盘失效寿命的第二级响应面。基于两级响应面模型，利用协同的方法得到涡轮盘概率寿命概率寿命分布和累积概率分布，即不同寿命下涡轮盘失效风险。

说明书

技术领域

本发明是一种针对航空发动机涡轮盘结构疲劳概率寿命及失效风险预测方法，它是一种基于分区和协同响应面的涡轮盘结构疲劳概率寿命及失效风险预测方法，属于航空航天发动机技术领域。

背景技术

航空发动机是一种极限产品，工作在高温、高压、高转速等的复杂载荷/环境下；发动机性能及安全性指标的提高，要求发动机重量轻、长寿命、高可靠性(如，安全飞行对发动机结构件则要求低的破坏概率，可达10^-5-10^-7次/飞行小时)。传统的可靠性分析方法只考虑危险点处单一失效模式对疲劳寿命的影响，预测精度较差难以指导工程实际应用。

涡轮盘疲劳寿命及失效风险预测时，未考虑涡轮盘不同部位失效模式的差异及其对疲劳寿命的影响，而本发明考虑涡轮盘不同部位对涡轮盘疲劳寿命和失效风险的影响，可以得到具有一定可靠度的涡轮盘疲劳寿命和失效风险，用于涡轮盘可靠性评估。

现有文献Enright M P,Moody J P,Sobotka J C.Optimal automated fracture risk assessment of 3d gas turbine engine components[C]//ASME Turbo Expo 2016:Turbomachinery Technical Conference and Exposition.American Society of Mechanical Engineers,2016:V07BT29A005-V07BT29A005.中提出了针对发动机轮盘的分区失效风险评估方法，文中从轮盘中缺陷出发进行分区研究，但在分区计算疲劳失效时并未考虑各区域失效模式的差异性，采用统一的疲劳失效模型进行研究；同时在得到各分区失效寿命后，没有考虑各失效之间的相关性对疲劳失效寿命的影响。本发明补充了该发明的不足，在对涡轮盘进行分区后，考虑各区域失效模式的差异性进行疲劳寿命计算，并利用协同响应面方法考虑不同区域失效模式相关性对疲劳寿命的影响。

发明内容

本发明技术解决方案：克服现有技术无法考虑涡轮盘不同部位失效模式的差异及其对疲劳寿命影响的不足，提供了一种考虑涡轮盘不同部位对涡轮盘疲劳寿命和失效风险影响的疲劳寿命预测方法，可以得到具有一定可靠度的涡轮盘疲劳寿命和失效风险，用于涡轮盘可靠性评估。

本发明一种基于分区的涡轮盘疲劳寿命及失效风险评估方法，实现步骤如下：

从涡轮盘几何、载荷、材料等随机变量的分散性出发，结合粉末涡轮盘应力、温度分布及有限元网格进行区域划分，通过引入协同响应面方法量化各区域失效风险相关性对整盘失效的影响，本发明提出基于分区的涡轮盘失效风险研究方法，其流程见图1。

(1)获取随机变量分布：针对涡轮盘可靠性分析中几何、载荷、材料三类随机变量，获取随机变量概率分布。

(2)建立初始分区：对涡轮盘进行有限元网格划分，结合几何、载荷、材料三类随机变量分布施加边界条件进行整盘有限元分析，获取应力分布和应变分布。以有限元网格为基本单元，结合随机变量中温度分布和有限元分析得到的应力分布，对应力和温度分布相近的相邻有限元网格单元进行合并，形成涡轮盘初始分区。

(3)分区失效风险分析：针对各分区中温度、载荷特征，判断每一个区域失效模式特征，选用相应失效的疲劳寿命模型结合应力/应变分布确定该区域疲劳失效寿命Y_i，并基于试验设计DOE方法建立随机输入变量和疲劳寿命的多项式代理模型。在此基础上利用蒙特卡罗算法对几何、载荷、材料三类随机变量进行抽样，获取初始分区中每一个区域的疲劳失效寿命的概率分布，从而得到相应的失效风险。

(4)分区方案优化：选定某一疲劳寿命，基于步骤(3)中初始分区中每一个区域对应的失效风险，对于高风险区域，基于应力和温度分布重新细化分区，对于低风险区域，合并应力和温度相邻区域，得到优化后的分区方案，共n分区。

(5)涡轮盘概率寿命分析：基于优化后分区方案，重复步骤(3)得到各分区输入随机变量和疲劳寿命的多项式代理模型，为第一级响应面；考虑到每一个区域中有一处发生失效则涡轮盘发生失效，建立每一个区域失效和涡轮盘失效的响应面关系为第二级响应面。基于两种响应面，利用协同的方法，得到整盘概率寿命概率寿命分布，从而得到不同寿命下涡轮盘失效风险，从而实现对涡轮盘的概率寿命评估。

所述步骤(1)中随机变量概率分布获取方法为：①通过灵敏度分析筛选涡轮盘关键几何尺寸，针对关键几何尺寸通过真实涡轮盘测量得到几何随机变量概率分布；②基于外场统计数据采用载荷谱转换方法获取载荷随机变量概率分布；③针对优化后分区方案中每一个区域对应涡轮盘的相应区域进行取样，开展取样件标准件试验，利用贝叶斯推理方法获取材料随机变量概率分布。

所述步骤(2)中对应力和温度分布相近的相邻有限元网格单元进行合并时以温度50℃-100℃和应力50Mpa-100Mpa为间隔合并相邻有限元网格。

所述步骤(3)中每一个区域失效特征包括低循环疲劳、高低周复合疲劳、蠕变-疲劳。

所述步骤(3)中，具体实施过程如下：针对初始分区中每个分区确定失效模式，利用有限元分析得到的应力分布和应变分布确定每个分区的疲劳失效寿命；然后利用试验设计DOE方法建立各分区随机变量和疲劳寿命间多项式代理模型，并利用蒙特卡罗算法对随机变量进行抽样，得到每一个区域的疲劳失效寿命的概率分布。

所述步骤(4)中分区优化方法为对于失效风险大于整体失效风险水平2/3的区域，基于应力分布和温度分布，缩小应力和温度分布的间隔，对区域进行重新细化分区，对于失效风险小于整体失效风险水平1/3的区域，扩大应力和温度分布的间隔，合并应力和温度相邻区域。

所述步骤(5)中协同方法为首先对第一级响应面中几何、载荷和材料三类随机因素抽样获取概率寿命分布，在此基础上，对第一级响应面的输出疲劳失效寿命进行抽样，作为第二级响应面的输入变量，最终得到涡轮盘概率寿命分布。

本发明与现有技术有益效果在于：现有文献Enright M P,Moody J P,Sobotka J C.Optimal automated fracture risk assessment of 3d gas turbine engine components[C]//ASME Turbo Expo 2016:Turbomachinery Technical Conference and Exposition.American Society of Mechanical Engineers,2016:V07BT29A005-V07BT29A005.中提出了针对发动机轮盘的分区失效风险评估方法，文中从轮盘中缺陷出发进行分区研究，但在分区计算疲劳失效时并未考虑各区域失效模式的差异性，采用统一的疲劳失效模型进行研究；同时在得到各分区失效寿命后，没有考虑各失效之间的相关性对疲劳失效寿命的影响。本发明补充了该发明的不足，在对涡轮盘进行分区后，分区各分区的失效模式，基于该区域失效模式所对应的概率寿命模型，进行概率疲劳寿命分析，在此基础上利用协同响应面方法考虑不同区域失效模式相关性对疲劳寿命的影响，得到准确度更高的概率疲劳寿命及失效风险。

附图说明

图1为本发明的方法流程图。

具体实施方式

下面结合附图，对本发明涡轮盘失效风险评估研究方法做进一步说明。

从涡轮盘几何、载荷、材料等随机变量的分散性出发，结合粉末涡轮盘应力、温度分布及有限元网格进行区域划分，通过引入协同响应面方法量化各区域失效风险相关性对整盘失效的影响，本发明提出基于分区的涡轮盘失效风险研究方法，其流程见图1。

(1)获取随机变量分布：针对涡轮盘可靠性分析中几何、载荷、材料等n个随机变量{X_n}，获取随机变量概率分布f(X_i)，具体方法为：①以涡轮盘最大等效应力为目标函数进行灵敏度分析，筛选出涡轮盘关键几何尺寸，针对关键几何尺寸通过真实涡轮盘测量得到尺寸统计数据，考虑到如果加工工艺过程稳定，而且严格按照尺寸公差中心加工，机械零件尺寸为随机变量且服从正态分布，选用K-S检验方法确定几何尺寸概率分布；②基于外场统计数据，采用循环矩阵法进行载荷谱转换方法获取温度、转速等载荷参量分布，具体方法为：首先去除伪读数，然后提取载荷峰谷值，筛选有效峰谷值，通过雨流计数法法进行载荷谱转化；③开展涡轮盘不同部位取样件标准件试验，基于确定部位的失效概率模型，利用贝叶斯推理方法结合已有材料手册上材料参数分布和试验数据，确定材料参数分布。

(2)建立初始分区：对涡轮盘进行有限元网格划分，结合几何、载荷、材料等随机变量分布施加边界条件进行涡轮盘有限元分析，获取应力和应变分布。以有限元网格为基本单元，结合随机变量中温度分布和有限元分析得到的应力分布，对应力和温度分布相近的相邻有限元网格单元进行合并，如以温度50℃-100℃和应力50Mpa-100Mpa为间隔合并相邻有限元网格，形成涡轮盘初始分区。

(3)分区失效风险分析：针对(2)得到的涡轮盘初始分区中温度分布范围，判断初始分区中每一个区域对应的失效模式，选用相应失效的疲劳寿命模型，如盘心区域为低循环疲劳选择Manson-Coffin公式、盘缘榫接部位为高低周复合疲劳选择Miner损伤模型、盘缘部位为蠕变-疲劳选择基于机械功的寿命模型，在寿命模型中代入步骤(2)有限元获取的应力和应变确定该区域疲劳失效寿命Y_i；然后采用试验设计DOE方法，如拉丁超立方抽样，对几何、载荷、材料等随机变量进行抽样，计算相应的该区域疲劳失效寿命分布，基于随机变量和疲劳失效寿命间对应关系，建立多项式代理模型即响应面Y＝F₀(X₁,X₂,…,X_n)；最后利用蒙特卡罗算法对响应面Y＝F(X₁,X₂,…,X_n)中几何、载荷、材料等随机变量X_i进行蒙特卡罗抽样，获取各区域疲劳失效的概率寿命分布，从而得到各区域疲劳寿命对应的失效风险。

(4)分区方案优化：选定某一疲劳寿命，如N＝10000循环，基于步骤(3)初始分区中每一个区域对应的失效风险，对于失效风险大于等于整体失效风险水平2/3的区域，基于应力分布和温度分布，缩小应力和温度分布的间隔，对区域进行重新细化分区，对于失效风险小于整体失效风险水平1/3的区域，扩大应力和温度分布的间隔，合并应力和温度相邻区域，得到优化后的分区方案，共n分区。

(5)涡轮盘概率寿命分析：基于优化后分区方案，重复步骤(3)得到各分区输入随机变量和疲劳失效寿命之间的多项式代理模型，为第一级响应面，即Y＝F(X₁,X₂,…,X_n)；考虑到优化后的分区方案中每一个区域中有一处发生失效则涡轮盘发生失效，建立各区域失效和涡轮盘失效的响应面关系为第二级响应面，即，Z＝G(Y₁,Y₂,…,Y_n)＝MIN(Y₁,Y₂,…,Y_n)，其中MIN表示最小值，即涡轮盘总寿命等于每一个区域失效寿命最小值。基于两级响应面模型，利用协同的方法，即首先对第一级响应面中几何、载荷、材料等随机变量抽样得到每一个区域疲劳失效寿命概率分布，然后对第一级响应面抽样得到的每一个区域失效寿命进行第二次抽样作为第二级响应面输入变量，从而得到涡轮盘疲劳失效寿命概率分布和疲劳失效寿命累积概率分布即为涡轮盘在。

提供以上实施例仅仅是为了描述本发明的目的，而并非要限制本发明的范围。本发明的范围由所附权利要求限定。不脱离本发明的精神和原理而做出的各种等同替换和修改，均应涵盖在本发明的范围之内。