一种基于二元随机退化过程的可靠性评估方法技术

本发明专利技术涉及复杂系统可靠性与安全，具体公开了一种基于二元随机退化过程的可靠性评估方法，通过退化测试或加速退化测试分别获取两种性能特征退化数据；判断两性能特征间相关性；采用边缘分布函数，表征性能特征的随机退化过程；采用退化数据，拟合边缘随机退化模型，采用贝叶斯方法估计参数；设定性能失效阈值，计算单一性能特征可靠度；两性能特征相互独立，按独立系统计算系统可靠度；两性能特征存在相关性，用copula函数建立两特征之间的相关性关系，并用贝叶斯方法对copula函数中的参数进行估计；最后利用copula函数计算特征间存在相关关系的系统可靠度。本发明专利技术能较精准地判断系统部件维护与更换时间，为维护与检修提供指导意见，从而节约成本。从而节约成本。从而节约成本。

【技术实现步骤摘要】
一种基于二元随机退化过程的可靠性评估方法


[0001]本专利技术专利涉及复杂系统可靠性与安全
，具体而言，涉及一种基于二元随机退化过程的可靠性评估方法。

技术介绍

[0002]在工程应用中，复杂的工程系统因其能满足无数的功能需求而被建造，然而随着时间的推移，无论是串联系统还是并联系统中的多个部件，由于磨损，缺口等原因会导致应力集中，从而产生疲劳，其功能会逐渐退化，从而降低系统的可靠性。对许多系统来说，退化是导致系统发生故障的主要原因之一。
[0003]近二十年来，已有许多基于简单系统可靠性的研究，这些研究往往依赖于完备数据。然而，随着科学技术的发展，许多复杂工程系统高度可靠且有更复杂的失效机制，在正常情况下，故障不会在短时间内发生，使得获取完备的寿命数据不仅具有技术挑战性，而且成本较高，因此传统的依赖于完备数据的简单系统可靠性评估方法无法满足要求。与此同时，系统的复杂性和失效机制的复杂性使得系统更容易产生大量故障，故障的产生与服役时间、次数、载荷大小及服役环境密切相关。而大多数产品的故障归因于一些潜在的退化机制，例如机械部件的磨损，电子器件的电阻和电池的容量，且具有时变的特征，当退化达到故障阈值时，最终导致产品故障。因此状态检测和退化分析等方法可以为现代复杂系统的可靠性分析提供一种有效的方法。
[0004]为此，本专利技术旨在通过对退化数据进行分析，建立退化模型，从而计算复杂系统可靠度，提供一种基于二元随机退化过程的可靠性评估方法。

技术实现思路

[0005]本专利技术的目的在于克服上述现有技术的不足，提供一种基于二元随机退化过程的可靠性评估方法，通过本专利技术的方法，能更好地对复杂系统进行故障预测和健康管理。
[0006]本专利技术的目的通过以下方案实现：
[0007]一种基于二元随机退化过程的可靠性评估方法，包括以下步骤：
[0008]S1、通过退化测试(DT)或加速退化测试(ADT)获取两种性能特征退化数据；
[0009]S2、判断两性能特征间的相关性；
[0010]S3、基于边缘分布函数表征两性能特征的随机退化过程；
[0011]S4、用退化数据来拟合随机退化模型并估计模型参数；
[0012]S5、设定两种性能的失效阈值并分别计算单一性能特征可靠度；
[0013]S6、若两性能特征相互独立，则按一般串联系统计算系统可靠度；若两性能特征存在相关性，则用copula函数来建立两特征之间的相关性关系，并对参数进行估计，最后通过copula函数计算考虑特征间相关关系的系统可靠度。
[0014]进一步地，步骤S1中通过退化测试(DT)或加速退化测试(ADT)获取两种性能特征退化数据为：
[0015][0016]式中，Y1，Y2为两种性能特征的退化数据集；y为t时刻，在应力s条件下测得的退化数据，j指示测试单元的次序，N为总的测试单元数；s'
j
和s'
j
'为第j个测试单元的压力变量值，如电流电压等；K
j
为第j个测试单元中的测量总次数，各个测试单元的K值可能有所不同；
[0017]由退化数据可得退化增量为：
[0018]Δy
ij
(t
k
)＝y
ij
(t
k
)
‑
y
ij
(t
k
‑1)；
[0019]其中，j为特征的序号，t
k
，t
k
‑1为在第j个测试单元下相邻的两个测试时刻。
[0020]进一步地，步骤S2中判断两种性能特征间的相关性的方法为：
[0021]采用Kendall秩相关系数(R)判断两特征间是否存在相关性，Kendall相关系数是指设有n个统计对象，每个对象有两个属性的系数；将所有统计对象按属性1取值排列，不失一般性，设此时属性2取值的排列是乱序的；设P为两个属性排列大小关系一致的统计对象的数目；则R的表达式为：
[0022][0023]其中，R＝0时，表示两特征独立；R＝1时，表示两特征呈现正相关；R＝
‑
1时，表示两特征负相关；所述判断两性能特征相关性的方法中，为消除数据随机性带来的影响，当|R|<0.3时，不考虑两特征间的相关性关系，即认为是独立的串联系统；当|R|>0.3时，判定两性能特征有相关性，即在计算系统可靠度时，考虑两特征之间的相关关系。
[0024]进一步地，步骤S3中基于边缘分布函数表征两性能特征的边缘随机退化过程为：
[0025]假设两性能的退化过程为随机退化过程，分别服从于某边缘退化分布函数，即：
[0026][0027]式中，ΔY
ij
(t
k
)为在第i个测试单元，第t
k
时刻，第j个性能的退化过程，MDP为边缘退化过程，为边缘退化过程(MDP)的参数；所述边缘退化过程用维纳(Wiener)分布函数表征，如下式所示：
[0028][0029]式中，α，β为维纳分布函数中待求的两个未知参数，控制随机退化过程中退化量的均值和方差，h(.)为边际模型中协变量的函数；ΔΛ(t
k
；γ
i
)＝Λ(t
k
；γ
i
)
‑
Λ(t
k
‑1；γ
i
)＝
t
k
^γ
i
‑
t
k
‑1^γ
i
，表示参数γ
i
对原始时间进行幂变换后的时间间隔。
[0030]进一步地，步骤S4中用退化数据来拟合边缘随机退化模型并估计模型参数的方法为：
[0031]基于维纳(Wiener)退化过程，通过贝叶斯方法来估计边缘模型参数，边缘退化过程的后验分布表达式为：
[0032][0033]式中，π(θ
i
Mar)为上述维纳(Wiener)分布函数中未知参数的先验概率密度函数。
[0034]进一步地，构建所述后验分布，采用马尔科夫链蒙特卡罗方法(MCMC)来生成参数的后验采样点，MCMC方法采用Metropolis算法，根据所述Metropolis算法，已知采样点θ
(s)
，生成新采点θ
(s+1)
的步骤为：
[0035](1)由θ
*
～J(θ|θ
(s)
)对θ
*
进行采样，其中J为参数为θ
(s)
的条件分布；
[0036](2)通过下式计算接收率r：
[0037][0038]其中，Data表示观测数据点，此处表示退化量；
[0039](3)假设参数u服从[0,1]上的均匀分布，对u进行采样，根据下式来判断是否更新采样点：
[0040][0041]得本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于二元随机退化过程的可靠性评估方法，其特征在于，包括以下步骤：S1、通过退化测试(DT)或加速退化测试(ADT)获取两种性能特征退化数据；S2、判断两性能特征间的相关性；S3、基于边缘分布函数表征两性能特征的随机退化过程；S4、用退化数据来拟合随机退化模型并估计模型参数；S5、设定两种性能的失效阈值并分别计算单一性能特征可靠度；S6、若两性能特征相互独立，则按一般串联系统计算系统可靠度；若两性能特征存在相关性，则用copula函数来建立两特征之间的相关性关系，并对参数进行估计，最后通过copula函数计算考虑特征间相关关系的系统可靠度。2.根据权利要求1所述的一种基于二元随机退化过程的可靠性评估方法，其特征在于，步骤S1中通过退化测试(DT)或加速退化测试(ADT)获取两种性能特征退化数据为：式中，Y1，Y2为两种性能特征的退化数据集；y为t时刻，在应力s条件下测得的退化数据，j指示测试单元的次序，N为总的测试单元数；s'
i
和s”i
为第i个测试单元的压力变量值，如电流电压等；K
j
为第j个测试单元中的测量总次数，各个测试单元的K
i
值可能有所不同；由退化数据可得退化增量为：其中，i为特征的序号，为在第i个测试单元下相邻的两个测试时刻。3.根据权利要求1所述的一种基于二元随机退化过程的可靠性评估方法，其特征在于，步骤S2中判断两种性能特征间的相关性的方法为：采用Kendall秩相关系数(R)判断两特征间是否存在相关性，Kendall相关系数是指设有n个统计对象，每个对象有两个属性的系数；将所有统计对象按属性1取值排列，不失一般性，设此时属性2取值的排列是乱序的；设P为两个属性排列大小关系一致的统计对象的数目；则R的表达式为：其中，R＝0时，表示两特征独立；R＝1时，表示两特征呈现正相关；R＝
‑
1时，表示两特征负相关；所述判断两性能特征相关性的方法中，为消除数据随机性带来的影响，当|R|<0.3时，不考虑两特征间的相关性关系，即认为是独立的串联系统；当|R|>0.3时，判定两性能特
征有相关性，即在计算系统可靠度时，考虑两特征之间的相关关系。4.根据权利要求1所述的一种基于二元随机退化过程的可靠性评估方法，其特征在于，步骤S3中基于边缘分布函数表征两性能特征的边缘随机退化过程为：假设两性能的退化过程为随机退化过程，分别服从于某边缘退化分布函数，即：式中，ΔY
ij
(t
k
)为在第j个测试单元，第t
k
时刻，第i个性能的退化过程，MDP为边缘退化过程，为边缘退化过程(MDP)的参数；所述边缘退化过程用维纳(Wiener)分布函数表征，如下式所示：式中，α，β为维纳分布函数中待求的两个未知参数，控制随机退化过程中退化量的均值和方差，h(
·
)为边际模型中协变量的函数；ΔΛ(t
k
；λ
i
)＝Λ(t
k
；γ
i
)
‑
Λ(t
k
‑1；γ
i
)＝t
k
^γ
i
‑
t
k
‑1^γ
i
，表示参数γ
i
对原始时间进行幂变换后的时间间隔。5.根据权利要求1所述的一种基于二元随机退化过程的可靠性评估方法，其特征在于，步骤S4中用退化数据来拟合边缘随机退化模型并估计模型参数的方法为：基于维纳(Wiener)退化过程，通过贝叶斯方法来估计边缘模型参数，边缘退化过程的后验分布表达式为：式中，为上述维纳(Wiener)分布函数中未知参数的先验概率密度函数。6.根据权利要求1所述的一种基于二元随机退化过程的可靠性评估方法，其特征在于，构建所述后验分布，采用马尔科夫链蒙特卡罗方法(MCMC)来...

【专利技术属性】
技术研发人员：李伟，汪硕文，孙传文，
申请(专利权)人：北京理工大学，
类型：发明
国别省市：