CN105718722A - 基于定时截尾寿命试验数据的产品可靠度估计方法 - Google Patents
基于定时截尾寿命试验数据的产品可靠度估计方法 Download PDFInfo
- Publication number
- CN105718722A CN105718722A CN201610030872.1A CN201610030872A CN105718722A CN 105718722 A CN105718722 A CN 105718722A CN 201610030872 A CN201610030872 A CN 201610030872A CN 105718722 A CN105718722 A CN 105718722A
- Authority
- CN
- China
- Prior art keywords
- eta
- tau
- reliability
- var
- solving
- 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
- 238000012360 testing method Methods 0.000 title claims abstract description 48
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000007476 Maximum Likelihood Methods 0.000 claims abstract description 22
- 239000011159 matrix material Substances 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 abstract description 2
- 238000002474 experimental method Methods 0.000 abstract description 2
- 238000011161 development Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Complex Calculations (AREA)
Abstract
本发明属于可靠性评估领域,具体涉及基于定时截尾寿命试验数据的产品可靠度估计方法,包括以下步骤:(S1)获取产品试验数据,并根据试验数据求解产品可靠度的点估计;(S2)在产品可靠度的点估计的基础上,结合极大似然估计的性质和增量方法,求解产品的可靠度的置信区间估计;本发明提出的方法无需生成大量的自助样本,因此比现有的bootstrap方法时间消耗要少。另外,经过大量的实验验证,采用本发明提出的方法计算得到的结果比bootstrap方法更加准确。
Description
技术领域
本发明属于可靠性评估领域,特指一种针对收集到的定时截尾寿命试验数据,对产品的可靠度进行估计的方法。
背景技术
可靠性是指产品在规定的条件下和规定的时间内完成规定功能的能力(具体见参考文献1:郭波,武小悦.系统可靠性分析[M].长沙:国防科技大学出版社,2002:5-6.),它是产品的固有属性,是衡量产品质量好坏的重要指标。可靠性的概率度量称为可靠度,有时也常用产品的寿命这一指标进行衡量。相应地,产品在规定的条件下丧失规定的功能,则称之为故障。随着现代科学技术的发展,构成产品的元器件越来越多,产品的规模越来越庞大,研制和生产费用越来越高,这使得产品的可靠性问题变得越来越重要。对产品的可靠度进行准确的估计,有助于我们及时了解产品的运行情况,做出正确的决策。
对可靠度的估计,通常是将产品的寿命视为随机变量,并认为产品的寿命服从某个特定分布,然后借助于数理统计理论进行分析。例如在理论分析和工程上,因为威布尔分布的良好特性,常用威布尔分布来拟合产品的寿命分布。威布尔分布的概率密度函数为:
其中t为产品的寿命,m、η为威布尔分布参数,具体为m为威布尔分布的形状参数,η为威布尔分布的尺度参数;exp表示以自然对数e为底的指数函数;在时刻τ处的可靠度为:
由此可知,假如需要估计产品在τ时刻处的可靠度,只要知道分布参数m和η的估计值和借助于式(2)即可求得可靠度R(τ)的估计值因此对可靠度的估计,关键在于对分布参数m和η的估计。
在实际中,往往需要首先利用一批试验样品进行可靠性寿命试验,收集试验样品的试验数据,然后借助于统计分析理论,对分布参数和可靠度进行估计。假如试验中针对所有样品收集得到的试验数据全部为故障数据,则称这种试验数据为完全样本,否则称之为截尾样本。随着科学技术的发展,当前工业上生产出的往往是高可靠、长寿命产品,且产品的成本造价越来越昂贵。因为成本和时间的限制,对一批产品进行可靠性寿命试验时,试图收集到完全样本是不现实的,在实际中,经过寿命试验收集到的一般是截尾样本。根据试验方式的不同,截尾样本又可分为两种,一种是定时截尾样本,另一种是定数截尾样本。定时截尾样本是指在寿命试验中,当试验达到预先设定的终止时刻后立即终止试验而收集到的试验数据。定数截尾样本是指在寿命试验中,当试验样品的故障数目达到预先设定值后立即终止试验而收集到的试验数据。本发明考虑的主要是定时截尾样本,即针对寿命试验预先设定一个终止时刻。一旦试验达到该时刻,就立即终止所有试验样品的试验。当样品在终止时刻前发生故障,则针对该样品收集到的试验时间为故障数据,反之收集到的试验时间就是截尾数据。所有样品的试验时间混合起来就构成定时截尾样本。
工程上对产品的可靠度估计的要求,往往既包括可靠度的点估计,又包括可靠度的置信区间估计。对可靠度的点估计进行求解是相对容易的,相关研究也很充分,因此可直接利用已有的方法对可靠度的点估计进行求解,比如可靠度的极大似然估计。但是对可靠度的置信区间估计进行求解往往比较困难,相关研究相对较少。当产品的寿命服从威布尔分布时,鉴于威布尔分布的复杂性,对置信区间估计进行求解时更为困难,尤其是针对定时截尾样本。Joarder等(具体见参考文献2:Joarder,A.,H.Krishna,andD.Kundu.InferencesonWeibullparameterswithconventionaltype-Icensoring[J].ComputationalStatistics&DataAnalysis,2011,55(1):1-11.)利用自助法(bootstrap)对可靠度的置信区间进行求解,这也是求解置信区间估计的一种常见思路,但该方法从本质上而言是近似方法,且需要生成大量的自助样本,所以相对耗时。本发明要解决的关键技术问题在于当部件的寿命服从威布尔分布时,借鉴已有方法得到产品可靠度的点估计后,如何进一步对产品可靠度的置信区间估计进行求解。
发明内容
为了解决上述技术问题,本发明基于极大似然估计的相关性质及增量方法,通过对现有技术进行整合和改进,提出了一种更快速更准确的求解可靠度置信区间的方法,能够用于当产品寿命服从威布尔分布时,得到产品的可靠度的点估计,并求解可靠度的置信区间估计,具体技术方案如下。
基于定时截尾寿命试验数据的产品可靠度估计方法,包括以下步骤:
(S1)获取产品试验数据,并根据试验数据求解产品可靠度的点估计;
(S11)获取产品试验数据,并计算该试验数据构成的截尾样本的似然函数;
假定针对n个产品进行可靠性寿命试验,并在T时刻终止,收集到的截尾样本中有r个故障数据,记为t1,t2,…,tr,r、n为整数且r≥1,n大于r,则截尾样本中剩余的(n-r)个样本值都为T;t1,t2,…,tr与(n-r)个T混合起来构成截尾样本,得到该样本的似然函数为:
其中f(ti;m,η)和R(T;m,η)分别为 i为自然数,i=1,2,…,r;m、η为威布尔分布参数,具体为m为威布尔分布的形状参数,η为威布尔分布的尺度参数;
(S12)根据极大似然估计方法,求解分布参数的极大似然估计值;
当分布参数(m,η)的取值令式(3)的自然对数值最大时,对应的取值即为(m,η)的极大似然估计值,按下式求解(m,η)的极大似然估计值:
其中分别为式(3)的自然对数关于m和η的一阶偏导数;化简式(4),得形状参数m的极大似然估计值是下式的根:
其中r是样本中的故障数据个数;针对式(5),借助于牛顿迭代公式求解,得到关于式(5)的求解,具体步骤如下:
1.利用下式m0值作为迭代的初值:
其中π为圆周率, 令
2.根据牛顿迭代公式,利用下式获取迭代一步之后的迭代值
其中g'(m0)是函数(7)的一阶导数。
3.继续迭代,并获取每步的迭代值j为自然数,直到满足条件|mj+1-mj|≤ε后终止迭代,其中ε是预先设定的迭代终止的误差;最终,可得分布参数m的极大似然估计值是
对于尺度参数η的极大似然估计值为
(S13)求解产品在时刻τ处的可靠度R(τ)的点估计为:
(S2)求解产品的可靠度的置信区间估计;
(S21)Fisher信息矩阵是由样本内所蕴含的分布参数的信息量所构成的;在威布尔分布下,按照信息量的严格定义,信息矩阵内的元素是很难求解的。
用分布参数的极大似然估计值近似信息矩阵中的元素,将信息矩阵表示为:
其中,
(S22)求解和的协方差矩阵C2×2,即信息矩阵FI2×2的逆
其中表示的方差,表示的方差;表示和的协方差;
(S23)运用增量方法求解得到的方差
其中
统计量的分布服从正态分布,即
其中
(S24)在置信水平(1-α)下,求解得到可靠度的双侧置信区间为:
单侧置信区间的置信下限为:
其中Uα/2,Uα和U1-α/2分别是标准正态分布U(0,1)的α/2,α和(1-α/2)分位数;可以通过查正态分布分位数表获得;α为取置信水平时设置的常量。
采用本发明获得的有益效果:本发明提出的方法无需生成大量的自助样本,因此比参考文献2中运用的bootstrap方法时间消耗要少。另外,经过大量的实验验证,用本发明提出的方法计算得到的结果比bootstrap方法更准确。本发明通过上述步骤很好地解决了当产品寿命服从威布尔分布时,依据定时截尾样本,如何对产品的可靠度的点估计和置信区间估计进行求解的问题。
附图说明
图1是本发明的流程图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明。
如图1所示,为本发明的流程图,本发明基于试验数据样本的似然函数,得到了可靠度的极大似然估计,并作为可靠度的点估计。进一步基于可靠度的点估计,结合极大似然估计的性质和增量方法,求得了可靠度的置信区间估计。
本实施例中假定对8个样品进行可靠性寿命试验,在T=36天处终止试验,则收集到7个故障数据,具体是4,5,8,11,20,29,35,数据的单位是天。剩余1个截尾数据为36天。7个故障数据和1个截尾数据构成的定时截尾样本即为本实施例中运用的数据。然后,按照本发明提出的技术方案,对该产品在时刻为1天处的可靠度进行估计。
根据本发明方法,求解产品可靠度的点估计。依据式(5),求解分布参数m的极大似然估计值并求得接着按照式(9)求解分布参数η的极大似然估计值并求得进一步按照式(10)求得产品在τ=1时的可靠度点估计为
取置信水平为0.9,并求解置信水平为0.9下的可靠度单侧置信区间的置信下限估计。按照式(12)求得和的协方差矩阵为 进一步可依据式(13)得到的方差为最终可按照式(16)求得可靠度的置信下限为RL=0.9512。若按照参考文献2中提出的bootstrap法,则算得R'L=0.9429。类似地,若取置信水平为0.95,可求得置信水平为0.95时可靠度的置信下限为RL=0.9426。若按照参考文献2中提出的bootstrap法,则算得R'L=0.9180。可看出,根据本发明算得的置信下限结果与bootstrap法算得的结果相比,更接近于点估计值;由此可知,本发明算得的结果更准确。
尽管上面是对本发明具体实施方案的完整描述,但是可以采取各种修改、变体和替换方案。这些等同方案和替换方案被包括在本发明的范围内。因此,本发明的范围不应该被限于所描述的实施方案,而是应该由所附权利要求书限定。
Claims (1)
1.基于定时截尾寿命试验数据的产品可靠度估计方法,其特征在于,包括以下步骤:
(S1)获取产品试验数据,并根据试验数据求解产品可靠度的点估计;
(S11)获取产品试验数据,并计算该试验数据构成的截尾样本的似然函数;
假定针对n个产品进行可靠性寿命试验,并在T时刻终止,收集到的截尾样本中有r个故障数据,记为t1,t2,…,tr,r、n为整数且r≥1,n大于r,则截尾样本中剩余的(n-r)个样本值都为T;t1,t2,…,tr与(n-r)个T混合起来构成截尾样本,得到该样本的似然函数为:
其中f(ti;m,η)和R(T;m,η)分别为ti≥0,m>0,η>0,i=1,2,…,r;m、η为威布尔分布参数,具体为m为威布尔分布的形状参数,η为威布尔分布的尺度参数;
(S12)根据极大似然估计方法,求解分布参数的极大似然估计值;
当分布参数(m,η)的取值令式(1)的自然对数值最大时,对应的取值即为(m,η)的极大似然估计值,按下式求解(m,η)的极大似然估计值:
其中分别为式(1)的自然对数关于m和η的一阶偏导数;化简式(2),得形状参数m的极大似然估计值是下式的根:
其中r是样本中的故障数据个数;针对式(3),借助于牛顿迭代公式求解得到
对于尺度参数η的极大似然估计值为
(S13)求解产品在时刻τ处的可靠度R(τ)的点估计为:
(S2)求解产品的可靠度的置信区间估计;
(S21)用分布参数的极大似然估计值近似信息矩阵中的元素,将信息矩阵表示为:
其中,
(S22)求解和的协方差矩阵C2×2,即信息矩阵FI2×2的逆
其中表示的方差,表示的方差;表示和的协方差;
(S23)运用增量方法求解得到的方差
其中
统计量的分布服从正态分布,即
其中
(S24)在置信水平(1-α)下,求解得到可靠度的双侧置信区间为:
单侧置信区间的置信下限为:
其中Uα/2,Uα和U1-α/2分别是标准正态分布U(0,1)的α/2,α和(1-α/2)分位数,α为取置信水平时设置的常量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610030872.1A CN105718722B (zh) | 2016-01-18 | 2016-01-18 | 基于定时截尾寿命试验数据的产品可靠度估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610030872.1A CN105718722B (zh) | 2016-01-18 | 2016-01-18 | 基于定时截尾寿命试验数据的产品可靠度估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105718722A true CN105718722A (zh) | 2016-06-29 |
CN105718722B CN105718722B (zh) | 2018-07-06 |
Family
ID=56147224
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610030872.1A Active CN105718722B (zh) | 2016-01-18 | 2016-01-18 | 基于定时截尾寿命试验数据的产品可靠度估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105718722B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106251044A (zh) * | 2016-07-21 | 2016-12-21 | 中国科学院数学与系统科学研究院 | 多批次成败型试验下产品贮存期评估的Buehler 方法 |
CN106547265A (zh) * | 2016-11-01 | 2017-03-29 | 中车株洲电力机车研究所有限公司 | 一种轨道交通电子控制装置的现场可靠性评估方法及系统 |
CN107632275A (zh) * | 2017-08-28 | 2018-01-26 | 厦门市产品质量监督检验院 | 一种led照明产品寿命和可靠性快速评价方法 |
CN108241790A (zh) * | 2018-01-24 | 2018-07-03 | 西南交通大学 | 基于恒加试验数据的Weibull型产品可靠性估计方法 |
CN108763041A (zh) * | 2018-05-03 | 2018-11-06 | 温州大学 | 一种基于热流量统计和概率分布的服务器热可靠度计算和监测方法 |
CN109325289A (zh) * | 2018-09-17 | 2019-02-12 | 中国人民解放军海军工程大学 | 一种估计电子件可靠性参数的方法 |
CN110260907A (zh) * | 2019-07-11 | 2019-09-20 | 中国电子科技集团公司第四十九研究所 | 一种用于传感器的温度应力无失效加速寿命试验方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103065052A (zh) * | 2013-01-07 | 2013-04-24 | 河南科技大学 | 一种测定机械产品理论寿命的方法 |
CN103218534A (zh) * | 2013-04-22 | 2013-07-24 | 北京航空航天大学 | 一种右截尾型寿命数据分布选择方法 |
CN103218533A (zh) * | 2013-04-22 | 2013-07-24 | 北京航空航天大学 | 一种指数分布的串联系统可靠度置信下限估计方法 |
CN104750982A (zh) * | 2015-03-19 | 2015-07-01 | 北京航空航天大学 | 一种基于比例分组再抽样的可靠度区间估计方法 |
-
2016
- 2016-01-18 CN CN201610030872.1A patent/CN105718722B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103065052A (zh) * | 2013-01-07 | 2013-04-24 | 河南科技大学 | 一种测定机械产品理论寿命的方法 |
CN103218534A (zh) * | 2013-04-22 | 2013-07-24 | 北京航空航天大学 | 一种右截尾型寿命数据分布选择方法 |
CN103218533A (zh) * | 2013-04-22 | 2013-07-24 | 北京航空航天大学 | 一种指数分布的串联系统可靠度置信下限估计方法 |
CN104750982A (zh) * | 2015-03-19 | 2015-07-01 | 北京航空航天大学 | 一种基于比例分组再抽样的可靠度区间估计方法 |
Non-Patent Citations (5)
Title |
---|
JIAN ZHANG ETAL: ""Statistical inference of component lifetimes with location-scale distributions from censored system failure data with known signature"", 《IEEE TRANSACTIONS ON RELIABILITY》 * |
K. KRISHNAMOORTHY ETAL: ""Confidence limits and prediction limits for a Weibull distribution based on the generalized variable approach"", 《JOURNAL OF STATISTICAL PLANNING AND INFERENCE》 * |
K.S. SULTAN ETAL: ""bayesian and maximum likelihood estimations of the inverse weibull parameters under progressive"", 《JOURNAL OF STATISTICAL COMPUTATION & SIMULATION》 * |
王智明: ""数控机床的可靠性评估与不完全预防维修及其应用"", 《中国博士论文全文数据库 工程科技Ⅰ辑》 * |
顾嘉麟 等: ""截尾数据下威布尔分布的参数估计问题"", 《哈尔滨理工大学学报》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106251044A (zh) * | 2016-07-21 | 2016-12-21 | 中国科学院数学与系统科学研究院 | 多批次成败型试验下产品贮存期评估的Buehler 方法 |
CN106251044B (zh) * | 2016-07-21 | 2021-09-10 | 中国科学院数学与系统科学研究院 | 多批次成败型试验下产品贮存期评估的Buehler方法 |
CN106547265A (zh) * | 2016-11-01 | 2017-03-29 | 中车株洲电力机车研究所有限公司 | 一种轨道交通电子控制装置的现场可靠性评估方法及系统 |
CN106547265B (zh) * | 2016-11-01 | 2019-01-29 | 中车株洲电力机车研究所有限公司 | 一种轨道交通电子控制装置的现场可靠性评估方法及系统 |
CN107632275A (zh) * | 2017-08-28 | 2018-01-26 | 厦门市产品质量监督检验院 | 一种led照明产品寿命和可靠性快速评价方法 |
CN108241790A (zh) * | 2018-01-24 | 2018-07-03 | 西南交通大学 | 基于恒加试验数据的Weibull型产品可靠性估计方法 |
CN108241790B (zh) * | 2018-01-24 | 2021-04-02 | 西南交通大学 | 基于恒加试验数据的Weibull型产品可靠性估计方法 |
CN108763041A (zh) * | 2018-05-03 | 2018-11-06 | 温州大学 | 一种基于热流量统计和概率分布的服务器热可靠度计算和监测方法 |
CN109325289A (zh) * | 2018-09-17 | 2019-02-12 | 中国人民解放军海军工程大学 | 一种估计电子件可靠性参数的方法 |
CN109325289B (zh) * | 2018-09-17 | 2023-03-10 | 中国人民解放军海军工程大学 | 一种估计电子件可靠性参数的方法 |
CN110260907A (zh) * | 2019-07-11 | 2019-09-20 | 中国电子科技集团公司第四十九研究所 | 一种用于传感器的温度应力无失效加速寿命试验方法 |
CN110260907B (zh) * | 2019-07-11 | 2021-07-16 | 中国电子科技集团公司第四十九研究所 | 一种用于传感器的温度应力无失效加速寿命试验方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105718722B (zh) | 2018-07-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105718722A (zh) | 基于定时截尾寿命试验数据的产品可靠度估计方法 | |
CN102262191B (zh) | 一种固体钽电解电容器寿命预测方法 | |
CN102252898A (zh) | 基于“寿命-应力”模型的电子产品加速寿命试验方法 | |
CN103279657B (zh) | 一种基于工程经验的产品加速退化试验方案设计方法 | |
CN107357275A (zh) | 非高斯工业过程故障检测方法及系统 | |
CN110260907B (zh) | 一种用于传感器的温度应力无失效加速寿命试验方法 | |
CN102779208A (zh) | 基于相对熵的序贯加速退化试验优化设计方法 | |
CN111458661A (zh) | 一种配电网线变关系诊断方法、装置及系统 | |
CN105069532A (zh) | 一种多应力多退化量步进加速退化试验方案优化设计方法 | |
CN108491637B (zh) | 基于无失效数据超小样本的产品寿命分布评估方法 | |
CN107220500B (zh) | 基于逆高斯过程的性能退化试验贝叶斯可靠性评估方法 | |
CN109253826B (zh) | 一种基于多退化样本数据融合的热量表剩余寿命预测方法 | |
CN102303377B (zh) | 一种橡胶门尼粘度的在线自动测量方法 | |
CN110414086A (zh) | 一种基于灵敏度的综合应力加速因子计算方法 | |
CN105426647A (zh) | 基于可靠度先验信息融合的冷备系统可靠度估计方法 | |
CN102661848A (zh) | 智能电表液晶器件可靠性关键故障特征的测定方法 | |
CN106407592A (zh) | 一种测定疲劳性能的小子样数据方法 | |
CN109857581A (zh) | 一种优化选择软件可靠性增长模型的方法 | |
CN103488826B (zh) | 基于经验加速模型的退化量分布参数建模外推方法 | |
CN112213681A (zh) | 一种电能表寿命预测方法及装置 | |
CN115951292A (zh) | 一种电子式电压互感器误差状态在线评估方法及装置 | |
CN105510745B (zh) | 一种故障录波数据故障起点检测方法 | |
CN102435816B (zh) | 一种moa阻性电流测试仪的工作方法 | |
CN111859296B (zh) | 一种基于装备使用期间的测试性指标评估方法及系统 | |
Sergici et al. | Measurement and verification principles for behavior-based efficiency programs |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |