一种考虑多部位失效相关的轮盘裂纹扩展可靠性分析方法

1.本发明属于航空航天发动机技术领域，具体涉及一种考虑多部位失效相关的轮盘裂纹扩展可靠性分析方法。


背景技术：

2.轮盘在载荷、几何、材料等随机因素作用下，盘心、偏心孔、辐板等典型应力集中部位将随机萌生裂纹，引发轮盘多部位失效故障。由于相同/相似载荷、几何随机变量，轮盘中多失效部位具有失效相关性，属于共因失效。现有的可靠性分析方法大多假设各失效完全独立或完全相关，前者忽略各失效部位的相关性，而后者利用最危险失效部位近似代表整体结构。对于多部位串联共因失效的轮盘，完全独立和完全相关方法仅能提供轮盘失效概率的上下界，均无法获取准确的失效概率。有必要研究轮盘盘心、偏心孔、辐板等多失效部位相关性表征方法，建立考虑多部位相关的寿命可靠性分析模型，从而更为精确地预测轮盘裂纹扩展寿命可靠性。
3.现有专利“一种基于拟蒙特卡洛抽样的涡轮盘裂纹扩展可靠性分析方法”，以最危险点为考核点，考虑材料和几何分散性实现涡轮盘概率寿命及可靠性评估，但是该方法无法考虑多危险部位及其相关性对轮盘可靠性的影响。现有文献“基于混合copula函数的风电功率相关性分析[j].季峰,蔡兴国,王俊.电力系统自动化,2014,38(02):1-5+32.”通过混合耦合(copula)函数量化输入变量相关性，但是该方法适用于相关变量维度低的情况，当相关变量维度高时，耦合参数预测困难，不具有工程实用价值。


技术实现要素：

[0004]
为了克服现有技术的不足，本发明提供一种考虑多部位失效相关的轮盘裂纹扩展可靠性分析方法，服务与支撑航空发动机寿命预测技术。
[0005]
本发明解决上述技术问题采用的技术方案为：进行轮盘多失效部位分组，利用分组耦合函数表征多部位相关性，基于试验数据确定耦合函数中耦合参数，实现考虑多部位失效相关的裂纹扩展寿命可靠性分析。具体实现步骤如下：
[0006]
(1)轮盘多失效部位分组:轮盘中存在盘心、偏心孔、辐板多处危险部位，通过分组方法将所述危险部位进行分组；
[0007]
(2)多失效部位失效相关性表征：为描述所述分组的相关性，通过重复建立相邻部位的二元耦合函数h和耦合参数θ，得到多失效部位失效的联合分布函数；利用极大似然估计法估计得到最终的二元耦合函数h和耦合参数θ；
[0008]
(3)多失效部位相关的裂纹扩展寿命可靠性分析：假设轮盘中共存在k个危险部位，各所述危险部位组成串联系统，基于步骤(2)所获取的最终的二元耦合函数h和耦合参数θ，计算轮盘可靠度r为：
[0009][0010]
其中，pr表示概率；i为危险部位编号；k表示危险部位数量；fi表示第i个危险部位的累积概率分布；xi表示第i个危险部位的裂纹扩展长度，xc表示临界裂纹扩展长度。
[0011]
进一步地，所述步骤(1)中，所述分组的具体步骤如下：
[0012]
①
开展轮盘有限元分析，确定轮盘最大应力，选择大于最大应力50％且应力集中系数大于2的部位作为轮盘的所述危险部位；
[0013]
②
根据所述危险部位的载荷和几何特征进行分类，具体规则为：按照盘缘到盘心的顺序将具有相同几何特征的危险部位分为一类，用a、b、c依次编号，其中a、b、c每类包含多个危险部位；
[0014]
③
a、b、c各类组内的单个危险部位从同一起点处以顺时针顺序，用1、2、3依次编号；
[0015]
④
根据
②‑③
中编号，对a、b、c各类组内的危险部位分别建立“藤树”结构；其中各危险部位为节点，相邻节点分成一组，通过边连接，每条边与一个二元耦合函数h关联。
[0016]
进一步地，所述步骤(2)中，所述最终的二元耦合函数h和耦合参数θ的确定步骤如下：
[0017]
①
获取多部位裂纹扩展寿命数据；
[0018]
②
采用边缘函数推断估计法估计所述危险部位极限状态函数的边缘累积分布函数；
[0019]
③
利用极大似然估计法进行二元耦合函数h和耦合参数θ估计，首先估计a、b、c各类组内的危险部位的二元耦合函数h和耦合参数θ，每次均通过极大似然估计法估计相邻所述边所关联的二元耦合函数h及耦合参数θ，最终汇聚到每个节点xa、xb、xc；
[0020]
④
在确定a、b、c系列的内部危险部位的二元耦合函数h和耦合参数θ后，a、b、c各类之间通过节点xa、xb、xc节点相连，利用极大似然估计法依次估计最终的二元耦合函数h和耦合参数θ。
[0021]
进一步地，所述步骤(2)中的第
②
步中的所述危险部位的极限状态函数为临界裂纹扩展长度与实际裂纹长度差。
[0022]
本发明与现有技术对比的优点在于：
[0023]
(1)本发明中通过耦合函数分组方法量化多部位失效相关性，可提高计算效率；
[0024]
(2)本发明可以量化多失效相关性对轮盘可靠性的影响，可提高可靠性预测精度。
附图说明
[0025]
图1为一种考虑多部位失效相关的轮盘裂纹扩展可靠性分析方法流程图；
[0026]
图2为多部位失效相关性分析步骤：(a)分类，(b)编号，(c)“藤树”结构；
[0027]
图3为分组估计耦合参数流程图；
[0028]
图4为分组估计耦合参数的分组示意图。
具体实施方式
[0029]
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
[0030]
下面结合附图1-4，对本发明的一种考虑多部位失效相关的轮盘裂纹扩展可靠性分析方法的技术方案做进一步说明。
[0031]
第一步：轮盘多失效部位分组。轮盘中存在盘心、偏心孔、辐板等多处危险部位，通过分组方法将危险部位分组。以图2所示轮盘结构为例，具体步骤为：
[0032]
①
开展轮盘有限元分析，确定轮盘最大应力，选择大于最大应力50％且应力集中系数大于2的部位作为轮盘的危险部位；
[0033]
②
根据危险部位的载荷和几何特征进行分类，具体规则为：按照盘缘到盘心的顺序将具有相同几何特征的危险部位分为一类，用a、b、c依次编号，其中每类可包含多个危险部位。图2轮盘结构的危险部位可分为三类，如图2(a)所示，其中辐板偏心孔为a类、盘心偏心孔为b类和盘心为c类；
[0034]
③
各类a、b、c组内的单个危险部位从同一起点处以顺时针顺序，用1、2、3依次编号。以辐板偏心孔b为例，编号为b1、b2、b3、
…b10
，如图2(b)所示；
[0035]
④
根据
②‑③
中编号，对a、b、c等各类分别建立“藤树”结构。以b类为例，具有10个危险部位，构建的藤树结构如图2(c)所示。其中各危险部位为“节点”，相邻节点分成一组，通过“边”连接，每条边与一个二元耦合(copula)函数h关联，相邻节点通过一个二元耦合函数h关联后汇聚到新的节，假设节点为b1和b2，则二元耦合函数h标记为h
b1,b2
，二元耦合函数h中的耦合参数θ标记为θ
b1,b2
点，汇聚后的新节点标记为b1,b2。
[0036]
第二步：多失效部位失效相关性表征。为描述分组相关性，通过重复建立相邻部位的二元耦合函数h，得到多失效部位失效的联合分布函数。具体步骤为：
[0037]
①
获取多部位裂纹扩展寿命数据t
ij
,其中i＝a、b、c
…
表示具有相同几何特征的危险部位的一类节点，j＝1、2、3
…
表示每类中危险部位数量。
[0038]
②
采用边缘函数推断估计法估计第i类危险节点第j个危险部位极限状态函数g
ij
的边缘累积分布函数(cdf)f
ij
，
[0039]fij
＝f(g
ij
)＝f(x
c-x
ij
)
ꢀꢀꢀꢀꢀꢀꢀ
(3)
[0040]
其中，xc表示临界裂纹长度，x
ij
表示第i类危险节点第j个危险部位的裂纹长度；
[0041]
③
利用极大似然估计法(mle)进行耦合参数θ估计，首先估计a、b、c等各类组内的危险部位的耦合参数θ，以b类为例，二元耦合函数h和耦合参数θ计算流程如图3虚框所示，第bj和b(j+1)个节点为一组，通过第j条边连接，通过极大似然估计法估计第j条边所关联的二元耦合函数h
bj,b(j+1)
及耦合参数θ
bj,b(j+1)
，得到考虑第j和第j+1个节点的联合累积分布函数为f bj,b(j+1)
，
[0042]fbj,b(j+1)
＝h
bj,b(j+1)
(f
bj
,f
b(j+1)
,θ
bj,b(j+1)
)
ꢀꢀꢀꢀꢀꢀꢀ
(4)
[0043]
其中二元耦合函数h包括高斯耦合函数、clayton耦合函数、gumbel耦合函数和frank耦合函数，通过aic准则筛选最佳的二元耦合函数h。每条边所关联的二元耦合函h及耦合参数θ确定后，汇聚成一个新的节点，相邻新节点继续分组后计算二元耦合函数h及耦
合参数θ。例如第bj和b(j+1)个节点汇聚到新节点bj,b(j+1)，第b(j+1)和b(j+2)个节点汇聚到新节点b(j+1),b(j+2)，新节点bj,b(j+1)和b(j+1),b(j+2)相邻为1组，联合累积分布函数为f bj,b(j+2)|b(j+1)
，
[0044][0045]
重复上述步骤遍历所有的相邻节点得到每个边所对应的最佳二元耦合函数h及耦合参数θ，最终将该类内部所有危险部位汇聚到一个节点xb，如图4所示，包含该类内部危险部位所有边的二元耦合函数h和耦合参数θ信息，节点xb的联合连接分布函数fb表示为：
[0046]
fb＝h
b1,b10|b2,b3,...,b9
(f
b1,b9|b2,b3,...,b8
,f
b2,b10|b3,b4,...,b9
,θ
b1,b10|b2,b3,...,b9
)
ꢀꢀ
(5)
[0047]
④
在确定a、b、c等各类内部危险部位的所有边的二元耦合函数h和耦合参数θ后，a、b、c各类之间通过节点xa、xb、xc等节点相连，利用极大似然估计法依次估计xa、xb、xc等节点相邻边对应的二元耦合函数h和耦合参数θ，最终将所有类汇聚到最终节点，包含所有边的二元耦合函数h和耦合参数θ信息。其中xa、xb、xc分别表示a、b、c等各类内部危险部位通过步骤
③
最终汇聚的节点，其中二元耦合函数h包括高斯耦合函数、clayton耦合函数、gumbel耦合函数和frank耦合函数，通过aic准则筛选最佳的二元耦合函数h。
[0048]
第三步：多部位相关的裂纹扩展寿命可靠性分析。假设轮盘中共存在k个危险部位，各危险部位组成串联系统，基于第二步第
④
步所获取的最终的二元耦合函数h和耦合参数θ，计算轮盘可靠度r为：
[0049][0050]
其中，pr表示概率；i为危险部位编号；k表示危险部位数量；fi表示第i个危险部位的累积概率分布；xi表示第i个危险部位的裂纹扩展长度，xc表示临界裂纹扩展长度。
[0051]
提供以上实施例仅仅是为了描述本发明的目的，而并非要限制本发明的范围。本发明的范围由所附权利要求限定。不脱离本发明的精神和原理而做出的各种等同替换和修改，均应涵盖在本发明的范围之内。