CN112784415A - 一种定数截尾加速寿命试验机理等同性检验与寿命预测方法 - Google Patents
一种定数截尾加速寿命试验机理等同性检验与寿命预测方法 Download PDFInfo
- Publication number
- CN112784415A CN112784415A CN202110085621.4A CN202110085621A CN112784415A CN 112784415 A CN112784415 A CN 112784415A CN 202110085621 A CN202110085621 A CN 202110085621A CN 112784415 A CN112784415 A CN 112784415A
- Authority
- CN
- China
- Prior art keywords
- distribution
- life
- product
- stress
- model
- 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
- 230000001133 acceleration Effects 0.000 title claims abstract description 60
- 238000000034 method Methods 0.000 title claims abstract description 58
- 238000012360 testing method Methods 0.000 title claims abstract description 57
- 230000007246 mechanism Effects 0.000 title claims abstract description 46
- 238000005520 cutting process Methods 0.000 title claims description 5
- 238000009826 distribution Methods 0.000 claims abstract description 110
- 238000013213 extrapolation Methods 0.000 claims abstract description 6
- 238000005315 distribution function Methods 0.000 claims description 10
- 230000001186 cumulative effect Effects 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 4
- 241000022852 Letis Species 0.000 claims description 2
- 230000001419 dependent effect Effects 0.000 claims description 2
- 238000012886 linear function Methods 0.000 claims description 2
- 238000004519 manufacturing process Methods 0.000 claims description 2
- 239000011159 matrix material Substances 0.000 claims description 2
- 239000002994 raw material Substances 0.000 claims description 2
- 238000001514 detection method Methods 0.000 abstract description 7
- 238000010998 test method Methods 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 3
- 238000007689 inspection Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000000528 statistical test Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/04—Ageing analysis or optimisation against ageing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
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)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Complex Calculations (AREA)
Abstract
本发明提供一种定数截尾加速寿命试验机理等同性检验与寿命预测方法,通过如下步骤实现:步骤一:计算相邻加速应力之间的经验对数似然比;步骤二:检验产品在各应力水平下失效机理是否等同;步骤三:用两步估计法对产品的寿命分布模型进行参数估计;步骤四:选择最优寿命分布模型;步骤五:确定模型选择结果,外推计算常规应力下产品的可靠寿命。本发明实现了对产品定数截尾加速寿命数据的失效机理等同性检验,保证了外推预测常规应力下产品失效寿命的准确性和有效性;它适用范围广,检测精度较高,它无需求解超越方程,在保证精度的前提下计算量大大降低,可操作性较强;本发明所述方法科学,工艺性好,具有广阔推广应用价值。
Description
技术领域
本发明涉及一种定数截尾加速寿命试验机理等同性检验与寿命预测方法,它是一种基于经验似然函数与估计方程的机理等同性检验方法和基于加速模型的寿命预测方法。它针对产品加速寿命试验得到的定数截尾失效数据,通过建立经验似然函数与估计方程,对数据的变异系数进行统计检验,判断在加速应力的作用下产品的失效机理是否发生了改变,并在失效机理等同的基础上,基于加速模型对产品的寿命分布模型进行参数估计,最终预测产品在常规应力下的可靠寿命。适用于加速试验失效机理检验与产品寿命评估等领域。
背景技术
加速试验技术用于评估高可靠、长寿命的产品,可以在节省时间和成本的情况下得到足够多的产品可靠性信息。根据失效数据类型的不同,通常将加速试验分为加速寿命试验(ALT)和加速退化试验(ADT)。合理的加速试验需要保证产品在不同应力水平下的失效机理不变,否则外推得到的结论无法真实反映产品在常规应力下的可靠性水平。因此,基于加速数据进行可靠性评估与寿命预测时,检验加速试验中产品在各应力水平下失效机理是否等同是十分重要的步骤,这会直接影响外推预测的有效性与准确性。
目前,针对ADT的失效机理等同性检验方法已经发展出了较为成熟的理论体系,但对于ALT的失效机理等同性检验缺少完整的方法。而基于经验似然函数与估计方程的非参数估计方法的不断完善,为ALT的失效机理等同性检验提供了一种思路。
基于此本发明提出了一种基于经验似然函数与估计方程的定数截尾加速寿命试验机理等同性检验方法,并在失效机理等同的基础上,基于加速模型给出产品寿命分布的参数估计以及外推预测常规应力下的可靠寿命。
发明内容
(1)本发明的目的:针对定数截尾加速寿命试验数据,提供一种失效机理等同性检验方法,即提供一种定数截尾加速寿命试验机理等同性检验与寿命预测方法,在失效机理等同的基础上,基于加速模型对产品的寿命分布模型进行参数估计,最终预测产品在常规应力下的可靠寿命。
(2)技术方案:
本发明需建立如下基本设置(即基本假设):
设置1产品加速寿命t服从威布尔分布、对数正态分布中的一种,各分布的累积失效函数分别为:
①威布尔分布:
其中,η和m分别为威布尔分布的尺度参数和形状参数;
②对数正态分布:
其中,θ和σ分别为对数正态分布的对数均值和对数标准差;
设置2若随机变量X的累积分布函数为:
F(x)=exp(-e-x) (3)
则称Exα为其分位点;其中,P{X≤Exα}=α;
基于上述假设,本发明一种定数截尾加速寿命试验机理等同性检验与寿命预测方法,通过如下步骤实现:
步骤一:计算相邻加速应力之间的经验对数似然比;
加速寿命试验共有k组应力水平,即s1<s2<…<sk,记第i组应力水平下观测得到的失效数据顺序统计量为:
其中,ni为样本数,ri为停止试验时的失效样品数;
定义第i组应力水平下产品寿命分布的变异系数为:
其中,μi和Si分别为第i组应力水平下产品寿命分布的均值和标准差;
将第h组应力水平下产品的变异系数记为Ch=C,第h+1组应力水平下产品的变异系数记为Ch+1=C+δ,则零假设为:
H0:δ=0 (6)
备选假设为:
H1:δ≠0 (7)
对于第h组与第h+1组失效数据,分别构建估计方程:
其中,μh为第h组应力水平下产品寿命分布的均值,Xh,i=th,i,i=1,2,…,rh为第h组应力水平下样本失效数据的顺序统计量;μh+1为第h+1组应力水平下产品寿命分布的均值,Yh+1,j=th+1,j,j=1,2,…,rh+1为第h+1组应力水平下样本失效数据的顺序统计量;
令pi=P(X=Xi),i=1,2,…,rh与qj=P(Y=Yj),j=1,2,…,rh+1为经验似然概率值,则经验对数似然比为:
其中,h=1,2,…,k-1;
步骤二:检验产品在各应力水平下失效机理是否等同;
由步骤一得到相邻加速应力之间的经验对数似然比之后,令:
当k→∞时,有:
则失效机理等同性的检验统计量定义为:
对于给定的置信水平1-α,若Tk<Ex1-α,则认为在各组加速应力水平下产品的失效机理没有发生改变;否则,失效机理发生了变化,并且,失效机理发生改变的位置h为:
步骤三:用两步估计法对产品的寿命分布模型进行参数估计;
Ⅰ.由最佳线性无偏估计法计算各组应力下产品寿命分布模型的参数估计值;
对于威布尔分布,有:
对于对数正态分布,有:
其中,C′(ni,ri,j)和D′(ni,ri,j)为对数正态分布的最佳线性无偏估计系数,可通过查表获得;
Ⅱ.确定产品寿命分布模型的参数估计结果;
根据失效机理等同性条件,在各组应力水平下,对于威布尔分布,形状参数的估计值应统计相等,对于对数正态分布,对数标准差的估计值应统计相等;由于样本原材料、制造过程中存在多种不确定性因素,上述参数的估计值并非完全相等,故采用加权平均的方法确定其估计值;
对于威布尔分布,有:
对于对数正态分布,有:
Ⅲ.结合加速模型,由最小二乘估计法计算加速模型中参数的估计值;
威布尔分布和对数正态分布的加速模型分别为:
①威布尔分布:
lnηi=c0+cρi (24)
其中,c0,c为待估的模型参数,ρi=ρ(si)为加速应力的函数;
②对数正态分布:
θi=b0+bρi (25)
其中,b0,b为待估的模型参数,ρi=ρ(si)为加速应力的函数;
确定加速模型后,通过最小二乘估计法计算加速模型中参数的估计值;
步骤四:选择最优寿命分布模型;
由步骤三得到产品不同寿命分布模型的参数估计值,分别计算其对数似然函数值与AIC值,选择AIC值最小的模型作为该加速寿命试验的最优寿命分布模型;
定数截尾数据的对数似然函数为:
其中,ni为样本数,ri为停止试验时的失效样品数,ti,j为第i组应力水平下失效数据的第j个顺序统计量,i=1,2,…,k,f(·)为产品寿命分布模型的概率密度函数,F(·)为累积分布函数;
AIC值的计算公式为:
AIC=-2lnL+2N (27)
其中,N为模型中未知参数的个数;
步骤五:确定模型选择结果,外推计算常规应力下产品的可靠寿命;
由步骤四确定最优的产品寿命分布模型后,对于给定的可靠度R,将常规应力条件代入加速模型,外推计算常规应力下产品的可靠寿命tR;
威布尔分布和对数正态分布的可靠寿命分别为:
①威布尔分布:
其中,η0为由加速模型得出的常规应力下产品的特征寿命;
②对数正态分布:
tR=exp(σZ1-R+θ0) (29)
其中,θ0为由加速模型得出的常规应力下产品寿命的对数均值,Z1-R为标准正态分布的下分位点。
其中,在步骤三中所述的“最佳线性无偏估计法”,是指:
(3)对θ的任意一个线性无偏估计量θ*,有:
设X为位置-尺度分布族随机变量,其分布函数可以表示为:
C(n,r,j)、D(n,r,j)可查表获得,对于不同的标准分布,将有不同的数值;
①威布尔分布的最佳线性无偏估计:
其中,gr,n为修偏系数,可通过查表获得;
②对数正态分布的最佳线性无偏估计:
其中,在步骤三中所述的“最小二乘估计法”,是指:
设因变量y与自变量x=(1,x1,x2,…,xp-1)T具有线性关系,即:
y=xTθ
其中,θ=(θ0,θ1,…,θp-1)T;
设一共进行n次试验,在第i组自变量xi=(1,xi1,xi2,…,xi,p-1)T下的观测值为yi,i=1,2,…,n;所有自变量xi构成矩阵X=(x1,x2,…,xn),所有yi构成向量Y=(y1,y2,…,yn)T;由于受随机干扰因素的影响,Y与X不严格呈线性关系,Y-Xθ代表误差;
(3)优点和功效:本发明一种定数截尾加速寿命试验机理等同性检验与寿命预测方法,其优点是:
①本发明基于经验似然函数与估计方程实现了对产品定数截尾加速寿命数据的失效机理等同性检验,保证了外推预测常规应力下产品失效寿命的准确性和有效性;
②本发明的失效机理等同性检验方法并不针对特定的寿命分布模型,适用范围广,检测精度较高;
③本发明的参数估计方法无需求解超越方程,在保证精度的前提下计算量大大降低,可操作性较强;
④本发明所述方法科学,工艺性好,具有广阔推广应用价值。
附图说明
图1本发明所述方法流程图。
具体实施方式
下面将结合实施例对本发明做进一步详细说明。
某种电路板产品,在正常环境条件下工作时的温度和相对湿度分别为298K和45%。共取88个样品分为4组,每组22个,分别对其按表1所示的加速应力进行加速寿命试验。
表1某电路板加速寿命试验应力水平
序号 | 温度/K | 相对湿度/% |
试验1 | 348 | 70 |
试验2 | 358 | 85 |
试验3 | 348 | 90 |
试验4 | 368 | 90 |
收集的定数截尾数据如表2所示:
表2电路板加速试验寿命数据
注:失效寿命数据的单位为小时。
本发明一种定数截尾加速寿命试验机理等同性检验与寿命预测方法,见图1所示,通过如下步骤实现:
步骤一:计算相邻加速应力之间的经验对数似然比;
根据式(13)计算各组相邻加速应力之间的经验对数似然比-2lnΛh,h=1,2,3,计算结果如表3所示:
表3相邻加速应力之间的经验对数似然比
-2lnΛ<sub>1</sub> | -2lnΛ<sub>2</sub> | -2lnΛ<sub>3</sub> |
0.0594 | 0.6835 | 0.6159 |
步骤二:检验产品在各应力水平下失效机理是否等同;
取显著性水平α=0.1,根据式(16)可得,Tk=-0.1469<Ex1-0.1=2.2504。因此,可判定在四组加速应力水平下该电路板的失效机理没有发生变化。
步骤三:用两步估计法对产品的寿命分布模型进行参数估计;
常见的加速应力中温度应力(T)和湿度应力(RH)符合阿伦尼斯模型的情况;
对于威布尔分布,有:
其中,c0,c1,c2为待估参数,Ti和RHi为第i组加速寿命试验中的加速应力,i=1,2,3,4;对于对数正态分布,有:
其中,b0,b1,b2为待估参数,Tj和RHj为第j组加速寿命试验中的加速应力,j=1,2,3,4;
①威布尔分布的参数估计结果
各应力水平下参数的最佳线性无偏估计结果如表4所示:
表4威布尔分布参数线性估计结果
②对数正态分布的参数估计结果
各应力水平下参数的最佳线性无偏估计结果如表5所示:
表5对数正态分布参数线性估计结果
步骤四:选择最优寿命分布模型;
根据式(26),分别计算威布尔分布和对数正态分布的对数似然函数值,再根据式(27),计算其AIC值,结果如下;
①威布尔分布:
对数似然函数值:log L=-238.73;AIC值:AIC=485.46;
②对数正态分布:
对数似然函数值:log L=-234.72;AIC值:AIC=477.43;
由以上结果可见,对数正态分布的AIC值小于威布尔分布的AIC值;
步骤五:确定模型选择结果,外推计算常规应力下产品的可靠寿命;
根据步骤四,选择对数正态分布作为产品的寿命分布模型;当给定可靠度R=0.95时,根据式(29)可得,产品的可靠寿命为tR=18378.26h=2.10y;
所以基于前面所述计算评估得产品在正常应力下的95%的可靠寿命为2.10年。
综上所述,本发明涉及一种定数截尾加速寿命试验机理等同性检验与寿命预测方法,它是一种基于经验似然函数与估计方程的机理等同性检验方法和基于加速模型的寿命预测方法。该方法的具体步骤是:一、计算各组相邻加速应力之间的经验对数似然比;二、检验产品在各应力水平下失效机理是否等同;三、用两步估计法对产品的寿命分布模型进行参数估计;四、选择最优寿命分布模型;五、确定模型选择结果,外推计算常规应力下产品的可靠寿命。本发明适用于加速试验失效机理检验与产品寿命评估等领域,具有较强的可操作性。
Claims (3)
1.一种定数截尾加速寿命试验机理等同性检验与寿命预测方法,需要设置如下:
设置1:产品加速寿命t服从威布尔分布、对数正态分布中的一种,各分布的累积失效函数分别为:
①威布尔分布:
其中,η和m分别为威布尔分布的尺度参数和形状参数;
②对数正态分布:
其中,θ和σ分别为对数正态分布的对数均值和对数标准差;
设置2:若随机变量X的累积分布函数为:
F(x)=exp(-e-x) (3)
则称Exα为其分位点;其中,P{X≤Exα}=α;
其特征在于,通过如下步骤实现:
步骤一:计算相邻加速应力之间的经验对数似然比;
加速寿命试验共有k组应力水平,即s1<s2<…<sk,记第i组应力水平下观测得到的失效数据顺序统计量为:
其中,ni为样本数,ri为停止试验时的失效样品数;
定义第i组应力水平下产品寿命分布的变异系数为:
其中,μi和Si分别为第i组应力水平下产品寿命分布的均值和标准差;
将第h组应力水平下产品的变异系数记为Ch=C,第h+1组应力水平下产品的变异系数记为Ch+1=C+δ,则为:
H0:δ=0 (6)
备选为:
H1:δ≠0 (7)
对于第h组与第h+1组失效数据,分别构建估计方程:
其中,μh为第h组应力水平下产品寿命分布的均值,Xh,i=th,i,i=1,2,…,rh为第h组应力水平下样本失效数据的顺序统计量;μh+1为第h+1组应力水平下产品寿命分布的均值,Yh+1,j=th+1,j,j=1,2,…,rh+1为第h+1组应力水平下样本失效数据的顺序统计量;
令pi=P(X=Xi),i=1,2,…,rh与qj=P(Y=Yj),j=1,2,…,rh+1为经验似然概率值,则经验对数似然比为:
其中,h=1,2,…,k-1;
步骤二:检验产品在各应力水平下失效机理是否等同;
由步骤一得到相邻加速应力之间的经验对数似然比之后,令:
当k→∞时,有:
则失效机理等同性的检验统计量定义为:
对于给定的置信水平1-α,若Tk<Ex1-α,则认为在各组加速应力水平下产品的失效机理没有发生改变;否则,失效机理发生了变化,并且,失效机理发生改变的位置h为:
步骤三:用两步估计法对产品的寿命分布模型进行参数估计;
3.1由最佳线性无偏估计法计算各组应力下产品寿命分布模型的参数估计值;
对于威布尔分布,有:
对于对数正态分布,有:
其中,C′(ni,ri,j)和D′(ni,ri,j)为对数正态分布的最佳线性无偏估计系数,可通过查表获得;
3.2确定产品寿命分布模型的参数估计结果;
根据失效机理等同性条件,在各组应力水平下,对于威布尔分布,形状参数的估计值应统计相等,对于对数正态分布,对数标准差的估计值应统计相等;由于样本原材料、制造过程中存在多种不确定性因素,上述参数的估计值并非完全相等,故采用加权平均的方法确定其估计值;
对于威布尔分布,有:
对于对数正态分布,有:
3.3结合加速模型,由最小二乘估计法计算加速模型中参数的估计值;
威布尔分布和对数正态分布的加速模型分别为:
①威布尔分布:
lnηi=c0+cρi (24)
其中,c0,c为待估的模型参数,ρi=ρ(si)为加速应力的函数;
②对数正态分布:
θi=b0+bρi (25)
其中,b0,b为待估的模型参数,ρi=ρ(si)为加速应力的函数;
确定加速模型后,通过最小二乘估计法计算加速模型中参数的估计值;
步骤四:选择最优寿命分布模型;
由步骤三得到产品不同寿命分布模型的参数估计值,分别计算其对数似然函数值与AIC值,选择AIC值最小的模型作为该加速寿命试验的最优寿命分布模型;
定数截尾数据的对数似然函数为:
其中,ni为样本数,ri为停止试验时的失效样品数,ti,j为第i组应力水平下失效数据的第j个顺序统计量,i=1,2,…,k,f(·)为产品寿命分布模型的概率密度函数,F(·)为累积分布函数;
AIC值的计算公式为:
AIC=-2lnL+2N (27)
其中,N为模型中未知参数的个数;
步骤五:确定模型选择结果,外推计算常规应力下产品的可靠寿命;
由步骤四确定最优的产品寿命分布模型后,对于给定的可靠度R,将应力条件代入加速模型,外推计算常规应力下产品的可靠寿命tR;
威布尔分布和对数正态分布的可靠寿命分别为:
①威布尔分布:
其中,η0为由加速模型得出的常规应力下产品的特征寿命;
②对数正态分布:
tR=exp(σZ1-R+θ0) (29)
其中,θ0为由加速模型得出的常规应力下产品寿命的对数均值,Z1-R为标准正态分布的下分位点。
2.根据权利要求1所述的一种定数截尾加速寿命试验机理等同性检验与寿命预测方法,其特征在于:在步骤三中所述的“最佳线性无偏估计法”,是指:
(3)对θ的任意一个线性无偏估计量θ*,有:
设X为位置-尺度分布族随机变量,其分布函数可以表示为:
①威布尔分布的最佳线性无偏估计:
其中,gr,n为修偏系数,通过查表获得;
②对数正态分布的最佳线性无偏估计:
3.根据权利要求1所述的一种定数截尾加速寿命试验机理等同性检验与寿命预测方法,其特征在于:在步骤三中所述的“最小二乘估计法”,是指:
设因变量y与自变量x=(1,x1,x2,…,xp-1)T具有线性关系,即:
y=xTθ
其中,θ=(θ0,θ1,…,θp-1)T;
设一共进行n次试验,在第i组自变量xi=(1,xi1,xi2,…,xi,p-1)T下的观测值为yi,i=1,2,…,n;所有自变量xi构成矩阵X=(x1,x2,…,xn),所有yi构成向量Y=(y1,y2,…,yn)T;由于受随机干扰因素的影响,Y与X不严格呈线性关系,Y-Xθ代表误差;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110085621.4A CN112784415B (zh) | 2021-01-22 | 2021-01-22 | 一种定数截尾加速寿命试验机理等同性检验与寿命预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110085621.4A CN112784415B (zh) | 2021-01-22 | 2021-01-22 | 一种定数截尾加速寿命试验机理等同性检验与寿命预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112784415A true CN112784415A (zh) | 2021-05-11 |
CN112784415B CN112784415B (zh) | 2022-08-16 |
Family
ID=75758500
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110085621.4A Active CN112784415B (zh) | 2021-01-22 | 2021-01-22 | 一种定数截尾加速寿命试验机理等同性检验与寿命预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112784415B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113221286A (zh) * | 2021-05-24 | 2021-08-06 | 北京航空航天大学 | 一种基于加速寿命试验的盘轴螺纹连接器振动可靠性评估方法 |
CN113761755A (zh) * | 2021-09-24 | 2021-12-07 | 北京航空航天大学 | 考虑认知不确定性的温湿度双应力下的加速寿命分析方法 |
CN113792266A (zh) * | 2021-09-16 | 2021-12-14 | 西安太乙电子有限公司 | 恒定应力定时截尾加速寿命试验寿命评估方法及系统 |
CN113887009A (zh) * | 2021-06-24 | 2022-01-04 | 北京电子工程总体研究所 | 一种定数截尾试验加速失效机理一致性检验方法 |
CN114580118A (zh) * | 2022-03-25 | 2022-06-03 | 兰州空间技术物理研究所 | 一种离子推力器寿命及可靠性定量评估方法 |
CN114970157A (zh) * | 2022-05-26 | 2022-08-30 | 北京航空航天大学 | 电子产品在电压应力作用下的小样本试验寿命预测方法 |
CN114996928A (zh) * | 2022-05-26 | 2022-09-02 | 北京航空航天大学 | 电子产品的温压双应力小样本加速试验寿命预测方法 |
CN115544441A (zh) * | 2022-10-11 | 2022-12-30 | 成都海光微电子技术有限公司 | 一种高温工作寿命测试中回测时间点的确定方法及装置 |
CN116956638A (zh) * | 2023-09-12 | 2023-10-27 | 中国人民解放军63921部队 | 一种设备寿命数据联合拟合检验分析方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104680005A (zh) * | 2015-02-11 | 2015-06-03 | 北京航空航天大学 | 基于加速因子可行域选择的非平行贮存寿命试验评估方法 |
CN108241790A (zh) * | 2018-01-24 | 2018-07-03 | 西南交通大学 | 基于恒加试验数据的Weibull型产品可靠性估计方法 |
-
2021
- 2021-01-22 CN CN202110085621.4A patent/CN112784415B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104680005A (zh) * | 2015-02-11 | 2015-06-03 | 北京航空航天大学 | 基于加速因子可行域选择的非平行贮存寿命试验评估方法 |
CN108241790A (zh) * | 2018-01-24 | 2018-07-03 | 西南交通大学 | 基于恒加试验数据的Weibull型产品可靠性估计方法 |
Non-Patent Citations (5)
Title |
---|
A.REGATTIERI 等: "Reliability prediction of a mechanical component through accelerated life testing", 《DESTECH TRANSACTIONS ON ENGINEERING AND TECHNOLOGY RESEARCH》, 31 January 2018 (2018-01-31), pages 726 - 731 * |
王浩伟等: "融合加速退化和现场实测退化数据的剩余寿命预测方法", 《航空学报》, no. 12, 25 December 2014 (2014-12-25), pages 3350 - 3357 * |
王炳兴: "指数分布场合基于逐次定数截尾步加寿命试验数据的统计分析", 《系统工程理论与实践》, no. 11, 25 November 2005 (2005-11-25), pages 105 - 111 * |
葛峰等: "基于通用对数线性模型的多应力加速寿命试验预测方法", 《兵工自动化》, no. 10, 15 October 2020 (2020-10-15), pages 10 - 14 * |
邱雷颦: "基于广义指数分布与区间数据的加速寿命试验的研究", 《中国优秀博硕士学位论文全文数据库(硕士)基础科学辑》, 15 December 2017 (2017-12-15), pages 002 - 421 * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113221286A (zh) * | 2021-05-24 | 2021-08-06 | 北京航空航天大学 | 一种基于加速寿命试验的盘轴螺纹连接器振动可靠性评估方法 |
CN113221286B (zh) * | 2021-05-24 | 2022-08-19 | 北京航空航天大学 | 一种基于加速寿命试验的盘轴螺纹连接器振动可靠性评估方法 |
CN113887009A (zh) * | 2021-06-24 | 2022-01-04 | 北京电子工程总体研究所 | 一种定数截尾试验加速失效机理一致性检验方法 |
CN113792266A (zh) * | 2021-09-16 | 2021-12-14 | 西安太乙电子有限公司 | 恒定应力定时截尾加速寿命试验寿命评估方法及系统 |
CN113761755B (zh) * | 2021-09-24 | 2023-05-23 | 北京航空航天大学 | 考虑认知不确定性的温湿度双应力下的加速寿命分析方法 |
CN113761755A (zh) * | 2021-09-24 | 2021-12-07 | 北京航空航天大学 | 考虑认知不确定性的温湿度双应力下的加速寿命分析方法 |
CN114580118B (zh) * | 2022-03-25 | 2023-12-15 | 兰州空间技术物理研究所 | 一种离子推力器寿命及可靠性定量评估方法 |
CN114580118A (zh) * | 2022-03-25 | 2022-06-03 | 兰州空间技术物理研究所 | 一种离子推力器寿命及可靠性定量评估方法 |
CN114996928A (zh) * | 2022-05-26 | 2022-09-02 | 北京航空航天大学 | 电子产品的温压双应力小样本加速试验寿命预测方法 |
CN114970157A (zh) * | 2022-05-26 | 2022-08-30 | 北京航空航天大学 | 电子产品在电压应力作用下的小样本试验寿命预测方法 |
CN114970157B (zh) * | 2022-05-26 | 2024-05-28 | 北京航空航天大学 | 电子产品在电压应力作用下的小样本试验寿命预测方法 |
CN114996928B (zh) * | 2022-05-26 | 2024-05-28 | 北京航空航天大学 | 电子产品的温压双应力小样本加速试验寿命预测方法 |
CN115544441A (zh) * | 2022-10-11 | 2022-12-30 | 成都海光微电子技术有限公司 | 一种高温工作寿命测试中回测时间点的确定方法及装置 |
CN115544441B (zh) * | 2022-10-11 | 2023-09-08 | 成都海光微电子技术有限公司 | 一种高温工作寿命测试中回测时间点的确定方法及装置 |
CN116956638A (zh) * | 2023-09-12 | 2023-10-27 | 中国人民解放军63921部队 | 一种设备寿命数据联合拟合检验分析方法 |
CN116956638B (zh) * | 2023-09-12 | 2023-12-19 | 中国人民解放军63921部队 | 一种设备寿命数据联合拟合检验分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112784415B (zh) | 2022-08-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112784415B (zh) | 一种定数截尾加速寿命试验机理等同性检验与寿命预测方法 | |
CN112926144B (zh) | 一种多应力加速寿命试验耦合效应分析及寿命预测方法 | |
CN108388694B (zh) | 一种塑封光电耦合器贮存寿命预测方法 | |
CN105930976B (zh) | 基于加权理想点法的节点电压暂降严重程度综合评估方法 | |
CN111859658B (zh) | 一种产品贮存寿命与可靠性评估方法 | |
CN110197020B (zh) | 一种环境变化对水文干旱影响的分析方法 | |
CN105954667A (zh) | 一种基于加速退化试验的pcb绝缘寿命快速预测方法及系统 | |
CN108959745A (zh) | 产品性能变化及寿命高置信度小样本统计推断和确认方法 | |
CN110956014A (zh) | 一种pcb在不同导线间距和电压作用下的寿命预测方法 | |
CN109522517B (zh) | 一种电表失效率计算方法和系统 | |
CN109375143B (zh) | 一种确定智能电能表剩余寿命的方法 | |
CN111914424B (zh) | 一种基于短期测风资料的设计风速取值方法及系统 | |
CN113032999B (zh) | 医疗设备使用寿命的预测方法及装置 | |
CN112329272B (zh) | 综合光伏组件退化和寿命数据的剩余寿命拟合预测方法 | |
CN117131977B (zh) | 一种基于误判风险最小准则的径流预报样本集划分方法 | |
CN116773239A (zh) | 一种智能燃气表控制器可靠性寿命预测方法 | |
CN116957534A (zh) | 一种预测智能电表更换数量的方法 | |
CN112567241A (zh) | 环境传感器协同校准方法 | |
CN115935675A (zh) | 基于失效物理分析的铝电解电容器寿命快速预测方法 | |
CN114611325A (zh) | 一种基于威布尔分布曲线的智能电能表寿命预测方法 | |
CN109978349B (zh) | 一种武器装备指数型和联系统的可靠度置信下限计算方法 | |
CN115455647A (zh) | 一种基于随机过程的小卫星寿命评估方法及系统 | |
CN112989629B (zh) | 一种基于钽电容多性能退化的可靠性评估方法 | |
CN111814360B (zh) | 一种机载机电设备使用可靠性评估方法 | |
CN114970311A (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 |