CN106570281A - 基于相似产品信息的小子样产品贝叶斯可靠性评估方法 - Google Patents
基于相似产品信息的小子样产品贝叶斯可靠性评估方法 Download PDFInfo
- Publication number
- CN106570281A CN106570281A CN201610982338.0A CN201610982338A CN106570281A CN 106570281 A CN106570281 A CN 106570281A CN 201610982338 A CN201610982338 A CN 201610982338A CN 106570281 A CN106570281 A CN 106570281A
- Authority
- CN
- China
- Prior art keywords
- product
- sample
- parameter
- reliability
- bayesian
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Complex Calculations (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种基于相似产品信息的小子样产品贝叶斯可靠性评估方法,包含:S1、对相似产品进行相似度打分,确定先验样本集;S2、获取先验信息,包括相似产品的失效寿命数据与伪失效寿命数据;S3、确定先验分布形式;S4、相似产品失效机理一致性检验;S5、确定先验分布和后验分布,实现可靠度的贝叶斯估计。本发明有效增加了小子样产品的样本量,先验信息包括失效寿命数据和伪失效寿命数据,并将双参数贝叶斯可靠性估计问题转换成单参数贝叶斯可靠性估计,有效提高了产品贝叶斯可靠性评估的精度,简化了贝叶斯可靠性参数估计的计算方法。
Description
技术领域
本发明涉及一种小子样产品贝叶斯可靠性评估方法,具体是指一种基于相似产品信息的小子样产品贝叶斯可靠性评估方法,属于可靠性建模、可靠性评估的技术领域。
背景技术
卫星在航天领域发挥着非常重要的作用,其可靠性和寿命直接影响卫星运行姿态控制、测量精度等指标,是影响航天任务成败的关键。卫星在轨长期服役,这就要求卫星的有效载荷必须可靠性高、寿命长。因此,对于卫星有效载荷,在全寿命周期的各个阶段,必须开展对应的可靠性工作。然而,由于卫星有效载荷本身具有的造价高昂、系统复杂、小子样的特点,不能实现将大量产品用于寿命试验,从而缺乏相关的寿命数据用于可靠性分析和验证工作。其次,在轨卫星有效载荷处于失重、真空、强辐射、高温差等复杂、恶劣的环境中,地面模拟试验环境很难做到对真实运行环境的全覆盖,传统的可靠性评估技术难以准确反映产品失效的本质。因此,利用有限的失效寿命数据统计出来的失效率模型指导建立产品可靠性模型,对于产品的可靠性评估会存在较大的误差,无法准确判断产品的可靠性水平。
对于小子样失效寿命数据缺乏统计推断的问题,常规方法具有很大的局限性,目前有许多学者提出了新的小子样产品可靠性分析方法。其中,《基于最大熵试验法的小样本产品可靠性验证方法》,通过加大样本载荷或减少样本承载能力,使单个样本的试验熵尽可能大,达到用小样本验证产品高可靠性、高置信度水平的目的。《基于相似产品信息的某系统任务可靠性预计方法》,利用环境试验、可靠性研制试验等不同试验的相似产品试验信息,扩大了统计样本量,提高了小子样系统可靠性评估的可靠性。该方法可用于同一型号不同研制阶段的产品,也适用于具有一定继承性的新旧不同型号的产品。
对于小子样产品的可靠性统计分析问题,贝叶斯方法已得到越来越多的应用。贝叶斯统计推断方法把概率分布中的参数作为未知的服从一定概率分布的随机变量,因此可以充分利用大量的先验信息,结合小样本信息进行产品可靠性评估。贝叶斯方法在提高产品可靠性评估准确性的同时,可以减少试验样本、缩短试验时间从而节约试验经费、缩短研制周期。由于贝叶斯方法的优越性,许多学者在这一领域开展了大量新的研究工作。其中,有学者针对高可靠长寿命产品失效寿命数据越来越难以获得的现状,基于加速退化试验实时获取的大量的产品性能退化数据,结合贝叶斯方法提出《基于多源退化数据的贝叶斯可靠性综合评估方法》。也有学者提出了一种基于混合贝叶斯先验分布的可靠性验证测试方法,通过先验矩方法和最大熵方法分别求出共轭先验分布的两组不同参数,按第二类极大似然方法确定的两种先验分布的权重融合两组参数,进而得到最终的先验分布,比单纯使用其中一种方法得到的先验分布更加准确。还有学者建立了系统动态贝叶斯网络和系统失效模式的动态贝叶斯网络模型,并离散化转化成多个单连通或多连通静态贝叶斯网络,运用静态贝叶斯网络和桶消元法进行双向推导分别求解的方法求解航空发动机涡轮叶盘系统的失效率和各失效模式的故障率。
但目前贝叶斯可靠性评估方法的研究多集中在利用同型号或同类产品寿命数据或性能退化数据作为先验信息评估产品可靠性方面,未将相似产品寿命数据或性能退化数据纳入先验信息范畴。对于单批次小批量生产的产品,可以获取的同型号产品样本较少,缺乏稳定的寿命数据,利用常规的统计分析方法,可靠度模型参数容易存在较大的误差,不能很好的反映产品真实的可靠性。许多单批次小批量产品,如航天、航空类产品,由于自身的特殊性,在设计上往往具有很强的继承性,很多新产品是在老产品的基础上继承发展的,因此可以利用与本产品相似的已有产品的寿命数据来评估本产品的可靠性。通过对相似产品可靠性数据的挖掘获取更多的有效先验信息来构造先验分布,结合现场试验数据,建立产品可靠性参数的后验分布,实现融合相似产品可靠性信息的小子样产品贝叶斯可靠性评估。
发明内容
本发明的目的在于提供一种基于相似产品信息的小子样产品贝叶斯可靠性评估方法,有效增加了小子样产品的样本量,先验信息包括失效寿命数据和伪失效寿命数据,并将双参数贝叶斯可靠性估计问题转换成单参数贝叶斯可靠性估计,有效提高了产品贝叶斯可靠性评估的精度,简化了贝叶斯可靠性参数估计的计算方法。
为了达到上述目的,本发明提供一种基于相似产品信息的小子样产品贝叶斯可靠性评估方法,其中,需要进行可靠性评估的产品为本产品,与本产品相似的其他产品为相似产品,包含以下步骤:
S1、对不同的相似产品进行相似度打分,把具有相同相似度的相似产品归类为一个先验样本,确定先验样本集;
S2、获取先验信息,包括相似产品的失效寿命数据与利用性能退化失效模型得到的伪失效寿命数据;
S3、采用双参数威布尔分布为总体分布函数,构造范·蒙特福特统计量,验证先验样本的样品寿命服从双参数威布尔分布,确定先验分布形式;
S4、利用威布尔分布的形状参数表征产品失效机理,构造巴特利特检验统计量,验证相似产品失效机理的一致性;
S5、利用失效机理一致性原则,结合相似度转化成的先验样本权重,把双参数威布尔分布贝叶斯估计问题转化成单参数贝叶斯估计,利用贝叶斯公式获取后验分布,得到本产品可靠度的贝叶斯估计。
所述的S1中,具体包含以下步骤:
S11、将相似产品的相似程度取为nj个离散值,相似度记为sj,j=1,2,...nj,且
S12、记表示第i个专家zi对相似产品pk的相似度打分,i=1,2,...,ni,k=1,2,...,nk,对相似产品pk的相似度计算加权和,得到
S13、为建立样本集,若则认为相似产品pk的相似度为sj,将具有相同相似度sj的产品归类为一个先验样本Xj,从而形成先验样本集。
所述的S2中,具体包含以下步骤:
S21、记先验样本Xj的样品为 若样品在使用期间未发生失效,利用其性能退化数据及退化失效模型 由最小二乘估计法,计算得到:
估算出之后,确定退化失效模型;
S22、记性能失效阈值为利用方程求解t,记为即得到先验样本Xj中的样品的伪失效寿命数据;
S23、若先验样本Xj中的样品在使用期间已发生失效,失效寿命为记则得到先验样本Xj中所有样品的失效寿命数据分别为
S24、对先验样本Xj的所有样品的失效寿命数据进行大小排序,得到顺序统计量
所述的S3中,具体包含以下步骤:
S31、建立假设:H0表示先验样本X的样品寿命服从威布尔分布W(m,η),m表征失效机理的形状参数,η表征特征寿命的尺度参数;令:
Zg=(Yg-μ)/σ;
μ=lnη;
σ=1/m;
在假设H0成立的条件下,是极值分布FY(y)=1-exp(-exp((y-μ)/σ))的前ng个顺序统计量,是标准极值分布FZ(z)=1-exp(-exp(z))的前ng个顺序统计量;
S32、构造范·蒙特福特统计量为:
取n′g=[ng/2],则统计量为:
在假设H0成立的条件下,渐近服从分布F(2(ng-n′g-1),2n′g);
S33、对于给定的显著性水平α,如果有F<Fα/2(2(ng-n′g-1),2n′g),或者有F>F1-α/2(2(ng-n′g-1),2n′g),则认为假设H0不成立;否则,接收假设H0成立,即先验样本X的样品寿命服从威布尔分布W(m,η)。
所述的S4中,失效机理是指引起故障的物理的、化学的、生物的或其他的过程,是产品失效的内在原因;利用相似产品信息评估本产品的可靠性时,要求两个产品的失效机理一致,即先验样本Xj所服从的威布尔分布的尺度参数mj是相同的。
所述的S4中,具体包含以下步骤:
S41、建立假设:H1表示m1=m2=...=mnj;由于σj=1/mj,此假设等价于极值分布G(σj,μj)中各σj相等,即等价于假设:H1′表示
S42、记是σj的线性无偏估计,的方差为:
S43、按照近似服从自由度有为的χ2分布,构造巴特利特检验统计量:
S44、在假设H′1成立的条件下,B2/C近似服从自由度为nj-1的χ2分布;对于给定的显著性水平α,当时,认为假设H′0不成立;否则,认为假设H′0成立,即假设H1成立。
所述的S5中,具体包含以下步骤:
S51、记先验样本Xj的失效分布函数为:
进行数据变换得到:
μj=lnηj;
x=lnt;
令Bj=mjμj;则得到Aj=mjx-Bj;
利用失效分布函数和最小二乘估计法,根据先验样本Xj的失效寿命数据,获得参数估计值
S52、记为本产品的现场试验信息,其为本产品的现场试验的失效寿命数据或伪失效寿命数据的顺序统计量,按S4所述的方法对其进行威布尔分布假设的验证,并按S51所述的方法获取其形状参数估计值
S53、经过S4的威布尔分布假设的验证,相似产品与本产品失效机理一致,即则把相似产品的相似度sj转化成形状参数mj的权重,融合先验信息和现场试验信息,得到后验分布形状参数m的估计值
其中,为先验样本Xj服从的威布尔分布形状参数mj的估计值,为现场试验信息服从的威布尔分布形状参数m0的估计值;从而将双参数贝叶斯可靠性估计转换成单参数贝叶斯可靠性估计;
S54、此时后验分布形状参数参数m已知,则是以为参数的指数分布的前g0个顺序统计量,把参数η的贝叶斯可靠性估计转化成指数分布进行处理;
S55、根据贝叶斯公式,后验分布为:
采用共轭先验分布,取先验分布为Γ(a,b),则有:
π(λ)∝λa-1exp(-bλ);
其中,符号“∝”表示两边仅相差一个不依赖于λ的常数因子;此时后验分布为:
则参数λ的贝叶斯估计为:
由于可靠度R(t)=exp(-λt),lnR(t)=-λt,其是λ的线性函数,从而求得可靠度R(t)的后验分布为:
进而得到可靠度R(t)的贝叶斯估计,即后验期望为:
τ=(T+b)/t。
综上所述,本发明提供的基于相似产品信息的小子样产品贝叶斯可靠性评估方法,有效弥补了小子样产品在进行可靠性评估时样本量不足的问题,提高了对产品可靠性评估的精度;先验信息不仅包括失效寿命数据,还包括伪失效寿命数据,提高了贝叶斯可靠性评估的精度;将双参数贝叶斯可靠性估计问题转换成单参数贝叶斯可靠性估计,简化了贝叶斯可靠性参数估计的计算方法。
附图说明
图1为本发明中的基于相似产品信息的小子样产品贝叶斯可靠性评估方法的流程图。
具体实施方式
以下结合图1,详细说明本发明的优选实施例。
如图1所示,为本发明所述的基于相似产品信息的小子样产品贝叶斯可靠性评估方法,其中,需要进行可靠性评估的产品称为本产品,与本产品相似的其他产品称为相似产品,包含以下步骤:
S1、对不同的相似产品进行相似度打分,把具有相同相似度的相似产品归类为一个先验样本,确定先验样本集;
S2、获取先验信息,包括相似产品的失效寿命数据与利用性能退化失效模型得到的伪失效寿命数据;
S3、采用双参数威布尔分布为总体分布函数,构造范·蒙特福特统计量,验证先验样本的样品寿命服从双参数威布尔分布,确定先验分布形式;
S4、利用威布尔分布的形状参数表征产品失效机理,构造巴特利特检验统计量,验证相似产品失效机理的一致性;
S5、利用失效机理一致性原则,结合相似度转化成的先验样本权重,把双参数威布尔分布贝叶斯估计问题转化成单参数贝叶斯估计,利用贝叶斯公式获取后验分布,得到本产品可靠度的贝叶斯估计。
所述的S1中,在工程实际中,相关技术领域的专家对产品的设计、制造、功能和使用环境具有丰富的经验,对相似产品进行相似度打分具体是指:在充分考虑相似要素数量和相似要素相似度的基础上,对相似产品与本产品的相似程度给出一个相对可信性较高的评估;具体包含以下步骤:
S11、将相似产品的相似程度取为nj个离散值,相似度记为sj,j=1,2,...nj,且
S12、记表示第i个专家zi对相似产品pk的相似度打分,i=1,2,...,ni,k=1,2,...,nk,对相似产品pk的相似度计算加权和,得到
S13、为建立样本集,若则认为相似产品pk的相似度为sj,将具有相同相似度sj的产品归类为一个先验样本Xj,从而形成先验样本集。
所述的S2中,若相似产品在使用期间已发生失效,则可以直接获取该相似产品的失效寿命数据;若相似产品在使用期间未发生失效,将关键性能退化信息或与之相关的参数作为时间的函数,获取退化轨迹,使用退化失效模型或回归方程进行描述,当性能退化信息达到预先设定的失效阈值所需要的时间,即可获取该相似产品的伪失效寿命数据。
所述的S2中,具体包含以下步骤:
S21、记先验样本Xj的样品为 若样品在使用期间未发生失效,利用其性能退化数据及退化失效模型 由最小二乘估计法,计算得到:
估算出之后,确定退化失效模型;
S22、记性能失效阈值为利用方程求解t,记为即得到先验样本Xj中的样品的伪失效寿命数据;
S23、若先验样本Xj中的样品在使用期间已发生失效,失效寿命为记则得到先验样本Xj中所有样品的失效寿命数据分别为
S24、对先验样本Xj的所有样品的失效寿命数据进行大小排序,得到顺序统计量
所述的S3中,采用双参数威布尔分布作为总体分布函数,并进行各样本分布函数的假设检验。该双参数威布尔分布是可靠性定量分析中常用的一类典型寿命分布,在可靠性工程中得到广泛的应用。
所述的S3中,具体包含以下步骤:
S31、不失一般性,建立假设:H0表示先验样本X的样品寿命服从威布尔分布W(m,η),m表征失效机理的形状参数,η表征特征寿命的尺度参数;令:
Zg=(Yg-μ)/σ;
μ=lnη;
σ=1/m;
在假设H0成立的条件下,是极值分布FY(y)=1-exp(-exp((y-μ)/σ))的前ng个顺序统计量,是标准极值分布FZ(z)=1-exp(-exp(z))的前ng个顺序统计量;
S32、构造范·蒙特福特统计量为:
取n′g=[ng/2],则统计量为:
在假设H0成立的条件下,渐近服从分布F(2(ng-n′g-1),2n′g);
S33、对于给定的显著性水平α,如果有F<Fα/2(2(ng-n′g-1),2n′g),或者有F>F1-α/2(2(ng-n′g-1),2n′g),则认为假设H0不成立;否则,接收假设H0成立,即先验样本X的样品寿命服从威布尔分布W(m,η)。
所述的S4中,失效机理是指引起故障的物理的、化学的、生物的或其他的过程,是产品失效的内在原因。利用相似产品信息评估本产品的可靠性时,要求两个产品的失效机理一致,即先验样本Xj所服从的威布尔分布的尺度参数mj是相同的。
所述的S4中,具体包含以下步骤:
S41、建立假设:H1表示由于σj=1/mj,此假设等价于极值分布G(σj,μj)中各σj相等,即等价于假设:H1′表示
S42、记是σj的线性无偏估计,的方差为:
S43、按照近似服从自由度有为的χ2分布,构造巴特利特检验统计量:
S44、在假设H1′成立的条件下,B2/C近似服从自由度为nj-1的χ2分布;对于给定的显著性水平α,当时,认为假设H′0不成立;否则,认为假设H′0成立,即假设H1成立。
所述的S5中,具体包含以下步骤:
S51、记先验样本Xj的失效分布函数为:
进行数据变换得到:
μj=lnηj;
x=lnt;
令Bj=mjμj;则得到Aj=mjx-Bj;
利用失效分布函数和最小二乘估计法,根据先验样本Xj的失效寿命数据,获得参数估计值
S52、记为本产品的现场试验信息,其为本产品的现场试验的失效寿命数据或伪失效寿命数据的顺序统计量,按S4所述的方法对其进行威布尔分布假设的验证,并按S51所述的方法获取其形状参数估计值
S53、经过S4的威布尔分布假设的验证,相似产品与本产品失效机理一致,即则把相似产品的相似度sj转化成形状参数mj的权重,融合先验信息和现场试验信息,得到后验分布形状参数m的估计值
其中,为先验样本Xj服从的威布尔分布形状参数mj的估计值,为现场试验信息服从的威布尔分布形状参数m0的估计值;从而将双参数贝叶斯可靠性估计转换成单参数贝叶斯可靠性估计;
S54、此时后验分布形状参数参数m已知,则是以为参数的指数分布的前g0个顺序统计量,把参数η的贝叶斯可靠性估计转化成指数分布进行处理;
S55、根据贝叶斯公式,后验分布为:
采用共轭先验分布,取先验分布为Γ(a,b),则有:
π(λ)∝λa-1exp(-bλ);
其中,符号“∝”表示两边仅相差一个不依赖于λ的常数因子;此时后验分布为:
则参数λ的贝叶斯估计为:
由于可靠度R(t)=exp(-λt),lnR(t)=-λt,其是λ的线性函数,从而求得可靠度R(t)的后验分布为:
进而得到可靠度R(t)的贝叶斯估计,即后验期望为:
τ=(T+b)/t。
综上所述,本发明所提供的基于相似产品信息的小子样产品贝叶斯可靠性评估方法,通过确定已有的相似产品与需要进行可靠性评估的本产品的相似度;利用性能退化模型推导相似产品的伪失效寿命数据,拓展了贝叶斯先验信息的获取,结合已有的失效寿命数据,建立先验分布形式;根据失效机理一致性原则,结合由相似度转化的先验样本权重,把双参数贝叶斯可靠性估计问题转换成单参数贝叶斯可靠性估计;最后,根据贝叶斯公式融合先验分布和试验信息获取后验分布,获取可靠度的贝叶斯估计。
因此,本发明所提供的基于相似产品信息的小子样产品贝叶斯可靠性评估方法,具有以下优点和有益效果:
1、充分融合了专家经验、相似产品信息、现场试验信息,弥补了小子样产品在进行可靠性评估时样本量不足的问题,提高了对产品可靠性评估的精度。
2、先验信息不仅包括已有的失效寿命数据,还利用产品性能退化数据建立退化模型推导伪失效寿命数据,从而获取了更多的先验信息,提高了贝叶斯可靠性评估的精度。
3、基于失效机理一致性的原则,结合由相似度转化的先验样本权重,把双参数贝叶斯可靠性估计问题转换成单参数贝叶斯可靠性估计,简化了贝叶斯可靠性参数估计的计算方法。
尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。
Claims (7)
1.一种基于相似产品信息的小子样产品贝叶斯可靠性评估方法,需要进行可靠性评估的产品为本产品,与本产品相似的其他产品为相似产品,其特征在于,包含以下步骤:
S1、对不同的相似产品进行相似度打分,把具有相同相似度的相似产品归类为一个先验样本,确定先验样本集;
S2、获取先验信息,包括相似产品的失效寿命数据与利用性能退化失效模型得到的伪失效寿命数据;
S3、采用双参数威布尔分布为总体分布函数,构造范·蒙特福特统计量,验证先验样本的样品寿命服从双参数威布尔分布,确定先验分布形式;
S4、利用威布尔分布的形状参数表征产品失效机理,构造巴特利特检验统计量,验证相似产品失效机理的一致性;
S5、利用失效机理一致性原则,结合相似度转化成的先验样本权重,把双参数威布尔分布贝叶斯估计问题转化成单参数贝叶斯估计,利用贝叶斯公式获取后验分布,得到本产品可靠度的贝叶斯估计。
2.如权利要求1所述的基于相似产品信息的小子样产品贝叶斯可靠性评估方法,其特征在于,所述的S1中,具体包含以下步骤:
S11、将相似产品的相似程度取为nj个离散值,相似度记为sj,j=1,2,...nj,且
S12、记表示第i个专家zi对相似产品pk的相似度打分,i=1,2,...,ni,k=1,2,...,nk,对相似产品pk的相似度计算加权和,得到
S13、为建立样本集,若则认为相似产品pk的相似度为sj,将具有相同相似度sj的产品归类为一个先验样本Xj,从而形成先验样本集。
3.如权利要求2所述的基于相似产品信息的小子样产品贝叶斯可靠性评估方法,其特征在于,所述的S2中,具体包含以下步骤:
S21、记先验样本Xj的样品为 若样品在使用期间未发生失效,利用其性能退化数据及退化失效模型由最小二乘估计法,计算得到:
估算出之后,确定退化失效模型;
S22、记性能失效阈值为利用方程求解t,记为即得到先验样本Xj中的样品的伪失效寿命数据;
S23、若先验样本Xj中的样品在使用期间已发生失效,失效寿命为记则得到先验样本Xj中所有样品的失效寿命数据分别为
S24、对先验样本Xj的所有样品的失效寿命数据进行大小排序,得到顺序统计量
4.如权利要求3所述的基于相似产品信息的小子样产品贝叶斯可靠性评估方法,其特征在于,所述的S3中,具体包含以下步骤:
S31、建立假设:H0表示先验样本X的样品寿命服从威布尔分布W(m,η),m表征失效机理的形状参数,η表征特征寿命的尺度参数;令:
Zg=(Yg-μ)/σ;
μ=l nη;
σ=1/m;
在假设H0成立的条件下,是极值分布FY(y)=1-exp(-exp((y-μ)/σ))的前ng个顺序统计量,是标准极值分布FZ(z)=1-exp(-exp(z))的前ng个顺序统计量;
S32、构造范·蒙特福特统计量为:
取n′g=[ng/2],则统计量为:
在假设H0成立的条件下,渐近服从分布F(2(ng-n′g-1),2n′g);
S33、对于给定的显著性水平α,如果有F<Fα/2(2(ng-n′g-1),2n′g),或者有F>F1-α/2(2(ng-n′g-1),2n′g),则认为假设H0不成立;否则,接收假设H0成立,即先验样本X的样品寿命服从威布尔分布W(m,η)。
5.如权利要求4所述的基于相似产品信息的小子样产品贝叶斯可靠性评估方法,其特征在于,所述的S4中,失效机理是指引起故障的物理的、化学的、生物的或其他的过程,是产品失效的内在原因;利用相似产品信息评估本产品的可靠性时,要求两个产品的失效机理一致,即先验样本Xj所服从的威布尔分布的尺度参数mj是相同的。
6.如权利要求5所述的基于相似产品信息的小子样产品贝叶斯可靠性评估方法,其特征在于,所述的S4中,具体包含以下步骤:
S41、建立假设:H1表示由于σj=1/mj,此假设等价于极值分布G(σj,μj)中各σj相等,即等价于假设:H′1表示
S42、记是σj的线性无偏估计,的方差为:
S43、按照近似服从自由度有为的χ2分布,构造巴特利特检验统计量:
S44、在假设H′1成立的条件下,B2/C近似服从自由度为nj-1的χ2分布;对于给定的显著性水平α,当时,认为假设H′0不成立;否则,认为假设H′0成立,即假设H1成立。
7.如权利要求6所述的基于相似产品信息的小子样产品贝叶斯可靠性评估方法,其特征在于,所述的S5中,具体包含以下步骤:
S51、记先验样本Xj的失效分布函数为:
进行数据变换得到:
μj=l nηj;
x=l n t;
令Bj=mjμj;则得到Aj=mjx-Bj;
利用失效分布函数和最小二乘估计法,根据先验样本Xj的失效寿命数据,获得参数估计值
S52、记为本产品的现场试验信息,其为本产品的现场试验的失效寿命数据或伪失效寿命数据的顺序统计量,按S4所述的方法对其进行威布尔分布假设的验证,并按S51所述的方法获取其形状参数估计值
S53、经过S4的威布尔分布假设的验证,相似产品与本产品失效机理一致,即则把相似产品的相似度sj转化成形状参数mj的权重,融合先验信息和现场试验信息,得到后验分布形状参数m的估计值
其中,为先验样本Xj服从的威布尔分布形状参数mj的估计值,为现场试验信息服从的威布尔分布形状参数m0的估计值;从而将双参数贝叶斯可靠性估计转换成单参数贝叶斯可靠性估计;
S54、此时后验分布形状参数参数m已知,则是以为参数的指数分布的前g0个顺序统计量,把参数η的贝叶斯可靠性估计转化成指数分布进行处理;
S55、根据贝叶斯公式,后验分布为:
采用共轭先验分布,取先验分布为Γ(a,b),则有:
π(λ)∝λa-1exp(-bλ);
其中,符号“∝”表示两边仅相差一个不依赖于λ的常数因子;此时后验分布为:
则参数λ的贝叶斯估计为:
由于可靠度R(t)=exp(-λt),l n R(t)=-λt,其是λ的线性函数,从而求得可靠度R(t)的后验分布为:
进而得到可靠度R(t)的贝叶斯估计,即后验期望为:
τ=(T+b)/t。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610982338.0A CN106570281B (zh) | 2016-11-08 | 2016-11-08 | 基于相似产品信息的小子样产品贝叶斯可靠性评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610982338.0A CN106570281B (zh) | 2016-11-08 | 2016-11-08 | 基于相似产品信息的小子样产品贝叶斯可靠性评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106570281A true CN106570281A (zh) | 2017-04-19 |
CN106570281B CN106570281B (zh) | 2020-04-17 |
Family
ID=58540259
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610982338.0A Active CN106570281B (zh) | 2016-11-08 | 2016-11-08 | 基于相似产品信息的小子样产品贝叶斯可靠性评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106570281B (zh) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107220500A (zh) * | 2017-05-27 | 2017-09-29 | 上海无线电设备研究所 | 基于逆高斯过程的性能退化试验贝叶斯可靠性评估方法 |
CN107766884A (zh) * | 2017-10-19 | 2018-03-06 | 中国人民解放军国防科技大学 | 一种基于代表点优化的Bayes融合评估方法 |
CN109145460A (zh) * | 2018-08-27 | 2019-01-04 | 枘熠集成电路(上海)有限公司 | 一种半导体可靠性评估方法及装置 |
CN109614586A (zh) * | 2018-12-03 | 2019-04-12 | 中国人民解放军国防科技大学 | 基于贝叶斯理论和线性拟合的可靠度置信下限构建方法 |
CN110889189A (zh) * | 2018-09-10 | 2020-03-17 | 湖南银杏可靠性技术研究所有限公司 | 基于等退化增量时间间隔Spearman系数加速贮存与自然贮存退化数据一致性检验法 |
CN111125918A (zh) * | 2019-12-26 | 2020-05-08 | 中国人民解放军国防科技大学 | 基于相似寿命模型和相似寿命的复杂产品可靠性评定方法 |
CN111178713A (zh) * | 2019-12-19 | 2020-05-19 | 中国航天标准化研究所 | 一种基于贝叶斯网络的Ka相控阵天线在轨可靠性评估方法 |
CN111241743A (zh) * | 2020-01-02 | 2020-06-05 | 中国辐射防护研究院 | 一种基于贝叶斯方法的放射性物品运输容器评价方法 |
CN111811827A (zh) * | 2020-07-17 | 2020-10-23 | 中国人民解放军国防科技大学 | 基于瑞利分布的产品性能一致性检验方法 |
CN112199776A (zh) * | 2020-10-10 | 2021-01-08 | 中车大同电力机车有限公司 | 机车全寿命周期的评估方法、评估系统及优化方法 |
CN112528523A (zh) * | 2020-12-29 | 2021-03-19 | 中国人民解放军国防科技大学 | 失效信息已知下的卫星动量轮表决系统剩余寿命预测方法 |
CN113704988A (zh) * | 2021-08-20 | 2021-11-26 | 中国石油大学(华东) | 一种小样本下化工安全相关设备使用寿命快速判定方法 |
CN114692087A (zh) * | 2022-04-19 | 2022-07-01 | 中国航空发动机研究院 | 一种产品的可靠性评价方法、可读存储介质和电子设备 |
CN111125918B (zh) * | 2019-12-26 | 2024-05-28 | 中国人民解放军国防科技大学 | 基于相似寿命模型和相似寿命的复杂产品可靠性评定方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2003005248A2 (en) * | 2001-07-05 | 2003-01-16 | Motorola Inc | Using a bayesian belief network |
CN102830305A (zh) * | 2012-08-14 | 2012-12-19 | 广东电网公司佛山供电局 | Gis母线接头过热性故障模拟方法 |
CN103279657A (zh) * | 2013-05-21 | 2013-09-04 | 北京航空航天大学 | 一种基于工程经验的产品加速退化试验方案设计方法 |
-
2016
- 2016-11-08 CN CN201610982338.0A patent/CN106570281B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2003005248A2 (en) * | 2001-07-05 | 2003-01-16 | Motorola Inc | Using a bayesian belief network |
CN102830305A (zh) * | 2012-08-14 | 2012-12-19 | 广东电网公司佛山供电局 | Gis母线接头过热性故障模拟方法 |
CN103279657A (zh) * | 2013-05-21 | 2013-09-04 | 北京航空航天大学 | 一种基于工程经验的产品加速退化试验方案设计方法 |
Non-Patent Citations (4)
Title |
---|
程皖民: "基于小子样复杂信息集的可靠性评估方法及其引用研究", 《中国优秀博士学位论文辑·工程科技II辑》 * |
蒋喜 等: "基于伪寿命分布和Bayes法的电主轴可靠性对比", 《机械科学与技术》 * |
邓新文: "基于贝叶斯原理的小子样试验导航误差评定方法", 《探测与控制学报》 * |
马静 等: "基于多源信息融合技术的光纤陀螺可靠性估计", 《中国惯性技术学报》 * |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107220500A (zh) * | 2017-05-27 | 2017-09-29 | 上海无线电设备研究所 | 基于逆高斯过程的性能退化试验贝叶斯可靠性评估方法 |
CN107220500B (zh) * | 2017-05-27 | 2020-07-31 | 上海无线电设备研究所 | 基于逆高斯过程的性能退化试验贝叶斯可靠性评估方法 |
CN107766884A (zh) * | 2017-10-19 | 2018-03-06 | 中国人民解放军国防科技大学 | 一种基于代表点优化的Bayes融合评估方法 |
CN109145460A (zh) * | 2018-08-27 | 2019-01-04 | 枘熠集成电路(上海)有限公司 | 一种半导体可靠性评估方法及装置 |
CN110889189A (zh) * | 2018-09-10 | 2020-03-17 | 湖南银杏可靠性技术研究所有限公司 | 基于等退化增量时间间隔Spearman系数加速贮存与自然贮存退化数据一致性检验法 |
CN110889189B (zh) * | 2018-09-10 | 2020-12-11 | 湖南银杏可靠性技术研究所有限公司 | 一种产品在不同试验应力水平下的退化数据一致性检验法 |
CN109614586A (zh) * | 2018-12-03 | 2019-04-12 | 中国人民解放军国防科技大学 | 基于贝叶斯理论和线性拟合的可靠度置信下限构建方法 |
CN111178713A (zh) * | 2019-12-19 | 2020-05-19 | 中国航天标准化研究所 | 一种基于贝叶斯网络的Ka相控阵天线在轨可靠性评估方法 |
CN111178713B (zh) * | 2019-12-19 | 2022-11-18 | 中国航天标准化研究所 | 一种基于贝叶斯网络的Ka相控阵天线在轨可靠性评估方法 |
CN111125918A (zh) * | 2019-12-26 | 2020-05-08 | 中国人民解放军国防科技大学 | 基于相似寿命模型和相似寿命的复杂产品可靠性评定方法 |
CN111125918B (zh) * | 2019-12-26 | 2024-05-28 | 中国人民解放军国防科技大学 | 基于相似寿命模型和相似寿命的复杂产品可靠性评定方法 |
CN111241743A (zh) * | 2020-01-02 | 2020-06-05 | 中国辐射防护研究院 | 一种基于贝叶斯方法的放射性物品运输容器评价方法 |
CN111241743B (zh) * | 2020-01-02 | 2023-06-23 | 中国辐射防护研究院 | 一种基于贝叶斯方法的放射性物品运输容器评价方法 |
CN111811827B (zh) * | 2020-07-17 | 2022-09-30 | 中国人民解放军国防科技大学 | 基于瑞利分布的产品性能一致性检验方法 |
CN111811827A (zh) * | 2020-07-17 | 2020-10-23 | 中国人民解放军国防科技大学 | 基于瑞利分布的产品性能一致性检验方法 |
CN112199776A (zh) * | 2020-10-10 | 2021-01-08 | 中车大同电力机车有限公司 | 机车全寿命周期的评估方法、评估系统及优化方法 |
CN112528523A (zh) * | 2020-12-29 | 2021-03-19 | 中国人民解放军国防科技大学 | 失效信息已知下的卫星动量轮表决系统剩余寿命预测方法 |
CN113704988A (zh) * | 2021-08-20 | 2021-11-26 | 中国石油大学(华东) | 一种小样本下化工安全相关设备使用寿命快速判定方法 |
CN114692087A (zh) * | 2022-04-19 | 2022-07-01 | 中国航空发动机研究院 | 一种产品的可靠性评价方法、可读存储介质和电子设备 |
Also Published As
Publication number | Publication date |
---|---|
CN106570281B (zh) | 2020-04-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106570281A (zh) | 基于相似产品信息的小子样产品贝叶斯可靠性评估方法 | |
CN102270302B (zh) | 一种基于灰色支持向量机的多应力加速寿命试验预测方法 | |
CN101576443B (zh) | 基于灰色rbf神经网络的加速寿命试验寿命预测方法 | |
CN104750932B (zh) | 一种混合不确定性下基于代理模型的结构可靠性分析方法 | |
López-Vizcaíno et al. | Multinomial-based small area estimation of labour force indicators | |
CN110472846A (zh) | 核电厂热工水力安全分析最佳估算加不确定性方法 | |
CN109934337A (zh) | 一种基于集成lstm的航天器遥测数据异常的检测方法 | |
Abdel-Nasser et al. | A novel smart grid state estimation method based on neural networks | |
CN103942457A (zh) | 基于关联向量机回归的水质参数时间序列预测方法 | |
CN110031872A (zh) | 一种导航星座的系统效能综合评估方法 | |
CN103488869A (zh) | 一种最小二乘支持向量机的风力发电短期负荷预测方法 | |
CN105653450A (zh) | 基于改进遗传算法与Adaboost结合的软件缺陷数据特征选择方法 | |
CN104794735A (zh) | 基于变分贝叶斯期望最大化的扩展目标跟踪方法 | |
CN107423508B (zh) | 卫星导航系统试验验证与测试评估数学模型的建立方法 | |
Wang et al. | Multivalued decision diagram-based common cause failure analysis in phased-mission systems | |
CN106169124A (zh) | 系统级产品可靠性综合评估置信推断方法 | |
CN108879732A (zh) | 电力系统暂态稳定评估方法及装置 | |
CN105320805A (zh) | 一种皮卫星多源可靠性信息融合方法 | |
CN105471647A (zh) | 一种电力通信网故障定位方法 | |
Zhang et al. | Evolving Elman neural networks based state-of-health estimation for satellite lithium-ion batteries | |
CN109961172A (zh) | 一种基于统计模型检验的cps稀有事件概率预测方法 | |
Myren et al. | A comparison of Gaussian processes and neural networks for computer model emulation and calibration | |
Balakrishnan et al. | Step-stress tests for interval-censored data under gamma lifetime distribution | |
Pang et al. | Temporal dependence Mahalanobis distance for anomaly detection in multivariate spacecraft telemetry series | |
CN105426665B (zh) | 基于状态监测的动态可靠度确定方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |