CN108804806A - Weibull分布综合应力恒加试验中参数的简化MLE方法 - Google Patents
Weibull分布综合应力恒加试验中参数的简化MLE方法 Download PDFInfo
- Publication number
- CN108804806A CN108804806A CN201810577377.1A CN201810577377A CN108804806A CN 108804806 A CN108804806 A CN 108804806A CN 201810577377 A CN201810577377 A CN 201810577377A CN 108804806 A CN108804806 A CN 108804806A
- Authority
- CN
- China
- Prior art keywords
- stress
- weibull
- parameter
- maximum
- likelihood
- 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
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000009826 distribution Methods 0.000 claims abstract description 56
- 238000012360 testing method Methods 0.000 claims abstract description 29
- 230000035882 stress Effects 0.000 claims description 69
- 238000007476 Maximum Likelihood Methods 0.000 claims description 22
- 230000007246 mechanism Effects 0.000 claims description 11
- 230000001133 acceleration Effects 0.000 claims description 9
- 238000005457 optimization Methods 0.000 claims description 9
- 238000002474 experimental method Methods 0.000 claims description 8
- 241001269238 Data Species 0.000 claims description 6
- 238000006467 substitution reaction Methods 0.000 claims description 5
- 230000001186 cumulative effect Effects 0.000 claims description 4
- 238000005315 distribution function Methods 0.000 claims description 4
- 230000008646 thermal stress Effects 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 2
- 230000004913 activation Effects 0.000 claims description 2
- 238000007619 statistical method Methods 0.000 abstract description 9
- 238000012795 verification Methods 0.000 abstract 1
- 239000011159 matrix material Substances 0.000 description 13
- 238000013461 design Methods 0.000 description 8
- 238000011160 research Methods 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 5
- 238000007689 inspection Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000000750 progressive effect Effects 0.000 description 4
- 230000002146 bilateral effect Effects 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000000846 Bartlett's test Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- ODINCKMPIJJUCX-UHFFFAOYSA-N Calcium oxide Chemical compound [Ca]=O ODINCKMPIJJUCX-UHFFFAOYSA-N 0.000 description 1
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 238000013398 bayesian method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 235000014347 soups Nutrition 0.000 description 1
- 238000013179 statistical model Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
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
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/17—Function evaluation by approximation methods, e.g. inter- or extrapolation, smoothing, least mean square method
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Tests Of Electronic Circuits (AREA)
- Complex Calculations (AREA)
Abstract
本发明属于加速寿命试验技术领域,具体是涉及Weibull分布综合应力恒加试验中参数的简化MLE方法。本发明针对单应力ALT的试验周期较长和Weibull分布参数的最大似然估计难以求解问题,给出了Weibull分布综合应力恒加试验数据的简化MLE方法,并通过对航天电连接器综合应力恒加试验的Monte Carlo仿真失效时间进行统计分析,验证所给Weibull分布综合应力恒加试验数据的简化MLE方法与可靠性评估的可行性和有效性。
Description
技术领域
本发明属于加速寿命试验技术领域,具体是涉及Weibull分布综合应力恒加试验中参数的简化MLE方法。
背景技术
加速寿命试验(Accelerated life test,ALT)是在保持失效机理不变的条件下,通过加大试验应力来加速产品失效的一种寿命试验方法;其目的是运用加速模型外推加速应力下的试验数据,对产品在正常应力水平下的各种可靠性指标进行统计推断。ALT可缩短试验周期、提高试验效率、降低试验成本,ALT技术的引入为解决高可靠度长寿命产品的寿命与可靠性评估提供了新途径。ALT分为恒定应力ALT(简称恒加试验,CSALT),步进应力ALT(简称步加试验,SSALT)和序进应力ALT(简称序加试验);ALT的研究主要集中在加速模型、ALT的优化设计和ALT数据的统计分析三个方面。ALT数据的统计分析方法主要有经典统计推断方法和贝叶斯(Bayesian)统计推断方法,其中经典统计推断方法中MLE方法应用最为广泛。目前,不少学者在ALT数据的统计分析方面方面已作了大量研究工作。
NELSON(1990)最先对ALT进行了全面性研究,提出了Weibull分布、正态分布和对数正态分布下完全样本和截尾样本情形ALT数据的MLE方法和最小二乘估计方法,且给出模型参数的置信区间;WATKINS(1994)研究了Weibull分布下恒加试验的精简MLE理论,并采用所提出的精简MLE方法刻画了Weibull分布的尺度参数与试验应力间的对数线性关系,算例分析结果表明所提方法是有效可行的;WANG等(2000)采用MLE方法研究了恒加试验Weibull分布对数线性模型的参数估计问题,提出了恒加试验Weibull分布对数线性模型参数的MLE方法,真实试验数据分析结果表明所给方法是有效可行的;王炳兴等(2002)提出了定数截尾恒加试验Weibull分布参数的渐进无偏估计及渐进置信区间估计的估计理论;WATKINS等(2008)提出了Weibull分布下定数截尾恒加试验的精简MLE理论,给出了与似然函数相关的Fisher信息矩阵的精确计算方法,算例分析表明所提精简MLE方法较传统MLE提高了计算效率且降低了参数MLE值的标准偏差;汤银才等(2009)提出了三参数Weibull分布Bayesian估计的Laplace数值积分方法和Gibbs抽样模拟方法,模拟算例说明了所提方法的可行性和有效性;XU等(2012)研究了相依指数分布竞争失效情形下恒加试验的MLE方法,真实试验数据分析结果表明所提MLE方法是可行的;陈文华等(2012)在Weibull分布失效下,利用MLE理论以产品中位寿命MLE值渐进方差的均值和标准离差分别作为估计精度和稳定性的考核指标,提出了航天电连接器步加试验的模拟评价理论与方法,算例分析表明基于MLE理论的模拟评价方法具有可行性;武东等(2013)研究了Weibull分布下步加试验的Bayesian估计方法,Monte Carlo仿真算例说明了所提Bayesian估计方法的有效性;ZHENG等(2013)采用EM算法和最小二乘法对两参数广义指数分布逐步定数混合截尾恒加试验数据进行了统计分析研究,获得了参数的渐进无偏估计和用于构造参数置信区间的Fisher信息矩阵;张详坡(2014)研究了三参数Weibull分布竞争失效场合序加试验的MLE方法,实例分析结果表明MLE方法是正确可行的且具有很好的估计效果。这些研究基本上都集中在单应力ALT方面。
然而,在单应力ALT中,为使产品在试验过程中失效机理保持不变,试验应力水平不宜过高,这样势必会造成试验周期较长;另一方面,运用传统的MLE法对Weibull分布的ALT数据进行统计分析时,参数估计往往解析困难、计算冗繁,通常需要借助数值计算方法。因此,有必要在单应力ALT数据统计分析方法的基础上,研究综合应力ALT数据的统计分析方法。
发明内容
本发明的目的,就是针对上述问题,以航天电连接器为研究对象,针对单应力ALT的试验周期较长和Weibull分布的最大似然估计难以求解问题,提出Weibull分布下综合应力恒加试验的优化MLE方法,并通过对航天电连接器综合应力恒加试验的Monte Carlo仿真失效时间进行统计分析,验证所提Weibull分布下综合应力恒加试验优化MLE方法的可行性和有效性。
本发明的技术方案为:
Weibull分布综合应力恒加试验中参数的简化MLE方法,该方法用于加速寿命试验中,Weibull分布参数的极大似然估计;其特征在于,包括以下步骤:
a、采用热应力(温度)和非热应力(如电压、电流、湿度、振动等)同时作为加速应力进行试验,则建立产品寿命特征与试验应力间关系的广义艾林(Eyring)模型为:
其中,η为特征寿命;T为热应力(绝对温度),S为非热应力;A、B均为待定常数,E为激活能(eV),K=8.617×10-5eV/℃为波耳兹曼常数;
并设定:
(1)在各应力水平组合(Ti,Sj)下,产品的寿命tij服从双参数Weibull分布W(mij,ηij),即tij~W(mij,ηij),i=0,1,2,…,k,j=0,1,2,…,l,其累积分布函数为:
其中,mij>0为形状参数,ηij>0为特征寿命或尺度参数;
设定寿命tij~W(mij,ηij),则对数服从极值分布G(μij,σij),即yij~G(μij,σij),其累积分布函数为:
其中,μij=lnηij为极值分布G(μij,σij)的位置参数,σij=1/mij为极值分布G(μij,σij)的尺度参数;
(2)在各应力水平组合(Ti,Sj)下,产品的失效机理保持不变,即所有应力水平组合下Weibull分布的形状参数mij相同:这等价于极值分布G(μij,σij)的尺度参数σij相同:i=0,1,2,…,k,j=0,1,2,…,l;
(3)在各应力水平组合(Ti,Sj)下,Weibull分布的特征寿命ηij与应力水平组合(Ti,Sj)间满足对数线性广义艾林加速模型:
其中,γ0=lnA,γ1=E/103K和γ2=-B均为待估参数;φ(Sj)=lnSj均为转化应力水平,记
b、根据步骤a设定的条件,对产品进行定时截尾综合应力恒加试验:
设定试验总样本量为n,在应力水平组合(Ti,Sj)下,投入nij个样品进行定时截尾寿命试验,试验截止时间为τij,在[0,τij]内观测到rij个样品发生失效,失效时间依次为:
同时设定其余nij-rij个样品将在(τij,∞)内发生失效;
c、获取样本数据的对数似然函数:
根据产品的失效机理保持不变,和Weibull分布的概率密度函数:
可得应力水平组合(Ti,Sj)下样本数据的似然函数为:
从而应力水平组合(Ti,Sj)下样本数据的对数似然函数为:
可得所有应力水平组合下全部样本数据的对数似然函数为:
其中:
d、获得全部样本数据的优化对数似然函数:
定义非负整数a、b和c,即a,b,c>0,由于
利用Qa,b,c可求得对数似然函数lnL的一阶偏导数令:
可得:
将(9)式代入(8)可得全部样本数据的优化对数似然函数为:
lnL*=D0,0lnm+(m-1)Qe-m(γ1D1,0+γ2D0,1)-D0,0lnQ0,0,0 (10)
e、获得Weibull分布尺度参数的极大似然估计:
根据步骤d,获得公式(10)的一阶偏导数分别为:
lnL*的二阶偏导数分别为:
lnL*的二阶混合偏导数分别为:
通过令对数似然函数lnL*的一阶偏导数为零,然后再解方程组,即可获得m、γ1和γ2的极大似然估计值和然后再将和代入式(9)即可得γ0的极大似然估计值
根据极大似然估计的不变性,即可得公式(1)的极大似然估计为:
在应力水平组合(T0,S0)下,Weibull分布尺度参数的极大似然估计为:
其中,
可靠度为R的可靠寿命tR,00的极大似然估计为
可靠度R00(t)的极大似然估计值为
失效率λ00(t)的极大似然估计值为
平均寿命tE的极大似然估计值为
上述方案中,根据MLE理论,参数MLE值和的协方差矩阵为
式(12)中:协方差矩阵∑用Fisher信息矩阵F来估计,而Fisher信息矩阵F是对数似然函数lnL的负二阶偏导数的数学期望,即
当对数似然函数二阶偏导数的数学期望的精确分布难以确定时,Fisher信息阵F的元素可用对数似然函数lnL的负二阶偏导数值来近似,其计算公式为
可靠性统计模型参数的渐进置信区间:
(1)特征寿命η00的置信区间
根据MLE的渐进正态性,统计量的分布可以采用正态分布近似,即
其中,
方差与协方差可通过Fisher信息矩阵F的逆矩阵,即协方差矩阵∑得到。
由此得到μ00的置信度为1-α的双侧置信线为
从而由μ00=lnη00得到η00的置信度为1-α的双侧置信线为
(2)形状参数m的置信区间
根据MLE理论,统计量的分布亦可用正态分布近似,即
因此,在置信水平为1-α时,形状参数m的双侧置信线为
其中,的方差是协方差矩阵∑的元素,u1-α2是标准正态分布的1-α/2分位数。
(3)加速模型系数γi的置信区间
根据MLE理论,统计量的分布亦可用正态分布近似,即
给定置信水平1-α,对数线性加速模型系数γi的双侧置信线为
与
其中,的方差由协方差矩阵∑给出。
本发明的方案中,与现有技术最直接的差异为,现有技术求解参数是通过令对数似然函数lnL的4个一阶偏导数为零,然后采用Newton法解对数似然方程组求得参数的MLE值 和这种直接求解对数似然方程组的方法,计算效率十分低下。
本发明的有益效果为,通过对对数似然函数的优化,极大的提高了参数估计的效率,因此提高了加速寿命试验的效率,有利于产品的快速测试。
附图说明
图1为实施例在各应力水平组合下的Weibull分布直线;
图2为实施例在正常应力水平组合下产品的可靠度曲线;
具体实施方式
下面结合附图和实施例进一步详细描述本发明的技术方案:
实施例
本例为仿真Y11X系列航天电连接器的失效过程,以环境温度和振动为加速寿命试验的加速应力,并采用广义艾林(Eyring)模型为加速模型。有关Y11X系列航天电连接器加速试验的研究表明,在环境温度和振动应力综合作用下,Y11X系列航天电连接器失效机理的改变点为(158℃,1.0g2/Hz),此即Y11X系列航天电连接器的最高应力水平组合点。另一方面,根据国标GJB101A—1997,Y11X系列航天电连接器的正常应力水平组合可选取为(85℃,0.06g2/Hz)。
通过大量摸底试验数据的统计分析,以及参考现有参考文献在Y11X系列航天电连接器加速试验与可靠性方面的研究,选取
作为广义艾林(Eyring)—Weibull可靠性统计模型参数的真值。
当加速应力水平数k=l=4时,采用均匀设计表U4(42)来进行试验,当给定最高试验应力水平组合后,试验方案的设计变量为和φ2,即设计变量水平组合为若温度与振动应力的转化应力水平取为等间隔水平,即
则温度与振动应力均为4水平数时的均匀设计方案如表1所示:
表1综合应力恒加试验均匀设计方案
若各应力水平组合下试验的截尾时间均为50h,则均匀设计方案下航天电连接器的Monte Carlo仿真失效时间如表2所示:
表2 均匀设计方案的Monte Carlo仿真失效时间
可靠性统计模型的检验
表2列举的仿真失效时间是否满足可靠性统计模型的基本假设,需要进行检验。首先,Weibull分布的拟合优度检验和失效机理一致性的检验,可采用Weibull概率纸法:
采用中位秩公式:
F(tijh)=(h-0.3)/(nij+0.4),h=1,2,…,rij
估计产品的失效概率F(tijh),并将各应力水平组合下的试验数据[tijh,F(tijh)]画在同一张Weibull概率纸上,如图1所示。由图1容易发现,各应力水平组合下试验数据点的趋势均为直线,且数据点的趋势大致平行。因此可以认为在各温度与振动综合应力作用下,航天电连接器的Monte Carlo仿真失效时间均服从Weibull分布,且形状参数m保持不变,即失效机理具有一致性。
Weibull分布的拟合优度检验亦可采用范·蒙特福特检验法,其检验结果如表3所示:
表3 Weibull分布的拟合优度检验结果
由表3的结果可知,范·蒙特福特统计量F值均介于F分布的0.05分位点F0.05(f1,f2)与0.95分位点F0.95(f1,f2)之间,这表明表2列举的航天电连接器的仿真失效时间均服从Weibull分布。
其次,失效机理一致性的检验可采用巴特利特检验法;在显著性水平α=0.05下,巴特利特检验统计量的值B2/C=2.5508,它小于卡方分布临界值χ2(0.95,3)=7.8147。因此可以接受失效机理一致性的原假设,即可以认为在各温度与振动应力水平组合下Weibull分布的形状参数相等。
模型参数的MLE和渐进置信区间:
采用本发明的优化MLE方法,并运用Newton算法来估计可靠性统计模型的未知参数,γ1、γ2和m的初值为和Newton方法迭代运算过程中,参数γ1、γ2和m,对数似然函数lnL*,以及lnL*的所有一、二阶偏导数值的收敛过程如表4所示:
表4 Newton方法迭代运算过程
由表4的结果可知,经过7步迭代运算后,lnL*的一阶偏导数和均收敛到零,且参数γ1、γ2和m收敛到稳定值,由此可得γ1、γ2和m的MLE值为
再将MLE值和代入式(9)可得γ0的MLE值为
从而可得加速模型的估计为
而采用传统MLE方法,Netwon算法经过39步迭代后,参数γ0、γ1、γ2和m才能收敛到稳定值和
由此可得本发明的方法相对于传统方法极大的提高了效率。
已知MLE值和后,根据式(13)可获得与似然函数相关的Fisher信息矩阵
再由式(12)可得与的协方差矩阵
估计出协方差矩阵∑后,根据第4部分介绍的渐进置信区间理论,可容易获得参数的渐进置信区间。在置信水平1-α=90%下,参数γ0、γ1、γ2、m和正常应力水平组合下Weibull分布特征寿命η0x的MLE值,渐进置信区间如表5所示:
表5模型参数的MLE和90%置信区间
正常应力水平组合下可靠性指标的MLE:
利用加速模型(15)式可得正常应力水平组合下特征寿命η00的MLE值为
可靠度R00(t)的MLE值为
失效率λ00(t)的MLE值为
可靠度为R的可靠寿命tR,00的MLE为
平均寿命tE的MLE值为
tE=89.173h
可靠度估计的可靠度曲线如图2所示。
Claims (1)
1.Weibull分布综合应力恒加试验中参数的简化MLE方法,该方法用于加速寿命试验中,Weibull分布参数的极大似然估计;其特征在于,包括以下步骤:
a、采用热应力和非热应力同时作为加速应力进行试验,则建立产品寿命特征与试验应力间关系的模型为:
其中,η为特征寿命;T为热应力,S为非热应力;A、B均为待定常数,E为激活能,K=8.617×10-5eV/℃为波耳兹曼常数;
并设定:
(1)在各应力水平组合(Ti,Sj)下,产品的寿命tij服从双参数Weibull分布W(mij,ηij),即tij~W(mij,ηij),i=0,1,2,…,k,j=0,1,2,…,l,其累积分布函数为:
其中,mij>0为形状参数,ηij>0为特征寿命或尺度参数;
设定寿命tij~W(mij,ηij),则对数服从极值分布G(μij,σij),即yij~G(μij,σij),其累积分布函数为:
其中,μij=lnηij为极值分布G(μij,σij)的位置参数,σij=1/mij为极值分布G(μij,σij)的尺度参数;
(2)在各应力水平组合(Ti,Sj)下,产品的失效机理保持不变,即所有应力水平组合下Weibull分布的形状参数mij相同:这等价于极值分布G(μij,σij)的尺度参数σij相同:i=0,1,2,…,k,j=0,1,2,…,l;
(3)在各应力水平组合(Ti,Sj)下,Weibull分布的特征寿命ηij与应力水平组合(Ti,Sj)间满足对数线性广义艾林加速模型:
其中,γ0=lnA,γ1=E/103K和γ2=-B均为待估参数;φ(Sj)=lnSj均为转化应力水平,记
b、根据步骤a设定的条件,对产品进行定时截尾综合应力恒加试验:
设定试验总样本量为n,在应力水平组合(Ti,Sj)下,投入nij个样品进行定时截尾寿命试验,试验截止时间为τij,在[0,τij]内观测到rij个样品发生失效,失效时间依次为:
同时设定其余nij-rij个样品将在(τij,∞)内发生失效;
c、获取样本数据的对数似然函数:
根据产品的失效机理保持不变,和Weibull分布的概率密度函数:
可得应力水平组合(Ti,Sj)下样本数据的似然函数为:
从而应力水平组合(Ti,Sj)下样本数据的对数似然函数为:
可得所有应力水平组合下全部样本数据的对数似然函数为:
其中:
d、获得全部样本数据的优化对数似然函数:
定义非负整数a、b和c,即a,b,c>0,由于
利用Qa,b,c可求得对数似然函数lnL的一阶偏导数令:
可得:
将(9)式代入(8)可得全部样本数据的优化对数似然函数为:
lnL*=D0,0lnm+(m-1)Qe-m(γ1D1,0+γ2D0,1)-D0,0lnQ0,0,0 (10)
e、获得Weibull分布尺度参数的极大似然估计:
根据步骤d,获得公式(10)的一阶偏导数分别为:
lnL*的二阶偏导数分别为:
lnL*的二阶混合偏导数分别为:
通过令对数似然函数lnL*的一阶偏导数为零,然后再解方程组,即可获得m、γ1和γ2的极大似然估计值和然后再将和代入式(9)即可得γ0的极大似然估计值根据极大似然估计的不变性,即可得公式(1)的极大似然估计为:
在应力水平组合(T0,S0)下,Weibull分布尺度参数的极大似然估计为:
其中,
可靠度为R的可靠寿命tR,00的极大似然估计为
可靠度R00(t)的极大似然估计值为
失效率λ00(t)的极大似然估计值为
平均寿命tE的极大似然估计值为
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810577377.1A CN108804806B (zh) | 2018-06-05 | 2018-06-05 | Weibull分布综合应力恒加试验中参数的简化MLE方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810577377.1A CN108804806B (zh) | 2018-06-05 | 2018-06-05 | Weibull分布综合应力恒加试验中参数的简化MLE方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108804806A true CN108804806A (zh) | 2018-11-13 |
CN108804806B CN108804806B (zh) | 2022-03-08 |
Family
ID=64087444
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810577377.1A Expired - Fee Related CN108804806B (zh) | 2018-06-05 | 2018-06-05 | Weibull分布综合应力恒加试验中参数的简化MLE方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108804806B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110263472A (zh) * | 2019-06-27 | 2019-09-20 | 中国人民解放军国防科技大学 | 基于回归法综合寿命试验数据的机电产品可靠度评估方法 |
CN110426168A (zh) * | 2019-07-09 | 2019-11-08 | 中国航天标准化研究所 | 一种多维应力耦合的星上电子产品加速寿命试验方法 |
CN112446139A (zh) * | 2020-11-05 | 2021-03-05 | 西安电子科技大学 | 加速试验剖面优化方法、***、机电产品、介质及终端 |
CN112528505A (zh) * | 2020-12-14 | 2021-03-19 | 西南交通大学 | 一种指数分布型产品可靠性评估方法 |
CN113688492A (zh) * | 2020-05-19 | 2021-11-23 | 中车株洲电力机车研究所有限公司 | 基于步降综合应力的加速试验可靠性评估方法及*** |
CN113704973A (zh) * | 2021-08-03 | 2021-11-26 | 人本股份有限公司 | 基于vba的轴承寿命数据自动化处理装置 |
CN116736781A (zh) * | 2023-08-15 | 2023-09-12 | 国网浙江省电力有限公司杭州供电公司 | 一种工业自动化控制设备安全状态监测方法及装置 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100138267A1 (en) * | 2009-08-31 | 2010-06-03 | Sameer Vittal | System and method for wind turbine health management |
EP2224256A1 (en) * | 2009-02-25 | 2010-09-01 | C.R.F. Società Consortile per Azioni | Reliability estimate of an apparatus |
CN104077445A (zh) * | 2014-07-01 | 2014-10-01 | 北京航空航天大学 | 基于模糊理论的加速寿命试验统计分析方法 |
FR3020681A1 (fr) * | 2014-04-30 | 2015-11-06 | Commissariat Energie Atomique | Procede de determination de la dispersion de la tenacite et de la temperature de transition fragile-ductile d'un produit en acier soumis a des variations thermiques |
CN106021928A (zh) * | 2016-05-20 | 2016-10-12 | 中国人民解放军总参谋部第六十研究所 | 一种综合应力加速试验方法 |
CN107491694A (zh) * | 2017-08-29 | 2017-12-19 | 西南交通大学 | 用于量化评估scada***信息安全脆弱性的方法 |
-
2018
- 2018-06-05 CN CN201810577377.1A patent/CN108804806B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2224256A1 (en) * | 2009-02-25 | 2010-09-01 | C.R.F. Società Consortile per Azioni | Reliability estimate of an apparatus |
US20100138267A1 (en) * | 2009-08-31 | 2010-06-03 | Sameer Vittal | System and method for wind turbine health management |
FR3020681A1 (fr) * | 2014-04-30 | 2015-11-06 | Commissariat Energie Atomique | Procede de determination de la dispersion de la tenacite et de la temperature de transition fragile-ductile d'un produit en acier soumis a des variations thermiques |
CN104077445A (zh) * | 2014-07-01 | 2014-10-01 | 北京航空航天大学 | 基于模糊理论的加速寿命试验统计分析方法 |
CN106021928A (zh) * | 2016-05-20 | 2016-10-12 | 中国人民解放军总参谋部第六十研究所 | 一种综合应力加速试验方法 |
CN107491694A (zh) * | 2017-08-29 | 2017-12-19 | 西南交通大学 | 用于量化评估scada***信息安全脆弱性的方法 |
Non-Patent Citations (1)
Title |
---|
冯雪峰: "基于Weibull 分布的变压器恒加试验可靠性统计分析", 《重庆工商大学学报》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110263472B (zh) * | 2019-06-27 | 2021-01-29 | 中国人民解放军国防科技大学 | 基于回归法综合寿命试验数据的机电产品可靠度评估方法 |
CN110263472A (zh) * | 2019-06-27 | 2019-09-20 | 中国人民解放军国防科技大学 | 基于回归法综合寿命试验数据的机电产品可靠度评估方法 |
CN110426168A (zh) * | 2019-07-09 | 2019-11-08 | 中国航天标准化研究所 | 一种多维应力耦合的星上电子产品加速寿命试验方法 |
CN113688492B (zh) * | 2020-05-19 | 2024-04-02 | 中车株洲电力机车研究所有限公司 | 基于步降综合应力的加速试验可靠性评估方法及*** |
CN113688492A (zh) * | 2020-05-19 | 2021-11-23 | 中车株洲电力机车研究所有限公司 | 基于步降综合应力的加速试验可靠性评估方法及*** |
CN112446139A (zh) * | 2020-11-05 | 2021-03-05 | 西安电子科技大学 | 加速试验剖面优化方法、***、机电产品、介质及终端 |
CN112446139B (zh) * | 2020-11-05 | 2024-05-10 | 西安电子科技大学 | 加速试验剖面优化方法、***、机电产品、介质及终端 |
CN112528505A (zh) * | 2020-12-14 | 2021-03-19 | 西南交通大学 | 一种指数分布型产品可靠性评估方法 |
CN112528505B (zh) * | 2020-12-14 | 2022-03-25 | 西南交通大学 | 一种指数分布型产品可靠性评估方法 |
CN113704973A (zh) * | 2021-08-03 | 2021-11-26 | 人本股份有限公司 | 基于vba的轴承寿命数据自动化处理装置 |
CN113704973B (zh) * | 2021-08-03 | 2024-03-15 | 人本股份有限公司 | 基于vba的轴承寿命数据自动化处理装置 |
CN116736781B (zh) * | 2023-08-15 | 2023-11-03 | 国网浙江省电力有限公司杭州供电公司 | 一种工业自动化控制设备安全状态监测方法及装置 |
CN116736781A (zh) * | 2023-08-15 | 2023-09-12 | 国网浙江省电力有限公司杭州供电公司 | 一种工业自动化控制设备安全状态监测方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN108804806B (zh) | 2022-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108804806A (zh) | Weibull分布综合应力恒加试验中参数的简化MLE方法 | |
Zhang et al. | A novel approach of battery pack state of health estimation using artificial intelligence optimization algorithm | |
CN107238765A (zh) | 基于加速性能退化参数的led集成驱动电源可靠性分析方法 | |
CN113805064B (zh) | 基于深度学习的锂离子电池组健康状态预测方法 | |
Mathieu et al. | State estimation and control of heterogeneous thermostatically controlled loads for load following | |
CN108535656A (zh) | 基于pca-narx神经网络的锂离子电池剩余使用寿命预测方法和*** | |
CN108959178A (zh) | 无失效数据的可靠性置信限统计评估方法 | |
CN108170991A (zh) | 基于Weibull分布的综合应力恒加试验方案可靠性评价方法 | |
CN110082682A (zh) | 一种锂电池荷电状态估计方法 | |
CN113109717B (zh) | 一种基于特征曲线优化的锂电池荷电状态估算方法 | |
CN107479000A (zh) | 一种基于Box‑Cox变换与蒙特卡罗仿真的锂离子动力电池RUL预测方法 | |
CN106991074B (zh) | 基于加速因子不变原则的加速退化试验优化设计方法 | |
Ali et al. | Estimation of P (Y< X) in a four-parameter generalized gamma distribution | |
CN109507527B (zh) | 一种基于突触可塑性的snp***含分布式电源配电网故障定位方法 | |
CN103472373B (zh) | 一种评价自然污秽与人工污秽绝缘子外绝缘等价性的方法 | |
Gökpınar et al. | A new computational approach for testing equality of inverse Gaussian means under heterogeneity | |
Cao et al. | A flexible battery capacity estimation method based on partial voltage curves and polynomial fitting | |
Kong et al. | New short-long-term degradation model for precise battery health prognostics | |
CN106126798A (zh) | 磷酸铁锂蓄电池soc算法 | |
Zhu et al. | An improved reliability assessment method for lithium-ion battery system considering imbalanced current and uneven cooling | |
CN114280490A (zh) | 一种锂离子电池荷电状态估计方法及*** | |
CN109214112B (zh) | 一种燃料电池***的可靠性仿真分析方法及*** | |
CN116757874A (zh) | 一种基于群体智能分析的用户用能态势推演方法 | |
CN111898829A (zh) | 一种基于稀疏异常感知的边缘设备下的电量预测方法 | |
Dvorak et al. | Arrhenius-equation based approach for modelling lithium-ion battery aging effects |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220308 |
|
CF01 | Termination of patent right due to non-payment of annual fee |