CN112926144B - 一种多应力加速寿命试验耦合效应分析及寿命预测方法 - Google Patents
一种多应力加速寿命试验耦合效应分析及寿命预测方法 Download PDFInfo
- Publication number
- CN112926144B CN112926144B CN202110085635.6A CN202110085635A CN112926144B CN 112926144 B CN112926144 B CN 112926144B CN 202110085635 A CN202110085635 A CN 202110085635A CN 112926144 B CN112926144 B CN 112926144B
- Authority
- CN
- China
- Prior art keywords
- model
- stress
- acceleration
- life
- function
- 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.)
- Active
Links
- 238000012360 testing method Methods 0.000 title claims abstract description 71
- 238000000034 method Methods 0.000 title claims abstract description 57
- 230000001808 coupling effect Effects 0.000 title claims abstract description 39
- 230000001133 acceleration Effects 0.000 claims abstract description 102
- 238000009826 distribution Methods 0.000 claims abstract description 60
- 230000007246 mechanism Effects 0.000 claims abstract description 20
- 238000007476 Maximum Likelihood Methods 0.000 claims abstract description 15
- 239000011159 matrix material Substances 0.000 claims description 18
- 230000004913 activation Effects 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 7
- 230000001186 cumulative effect Effects 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 230000009466 transformation Effects 0.000 claims description 2
- 230000035882 stress Effects 0.000 description 68
- 238000010276 construction Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000006353 environmental stress Effects 0.000 description 1
- 238000011156 evaluation 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/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- 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
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明提供一种多应力加速寿命试验耦合效应分析及寿命预测方法,步骤如下:一:利用广义对数线性模型构建多应力加速方程;二:利用整体极大似然理论对模型参数进行估计;三:通过中位寿命的相对误差判断参数估计结果的合理性;四:基于失效机理一致性原则,开展加速模型和寿命分布模型选择;五:外推正常应力下产品的可靠寿命;本发明利用广义对数线性模型构建的多应力加速方程考虑到应力之间复杂的耦合效应,使模型更贴合实际;它采取整体极大似然估计方法,对参数的初值要求低,方法迭代快速简单,可操作性强;它基于失效机理一致性原则,选择最优的加速寿命试验模型,提高了寿命预测精度;方法科学,工艺性好,具有推广应用价值。
Description
技术领域
本发明涉及一种多应力加速寿命试验耦合效应分析及寿命预测方法,它是一种基于广义对数线性模型和失效机理一致性原理的多应力加速寿命试验预测方法。它针对多应力加速寿命试验数据,利用广义对数线性模型构建包含耦合效应的加速方程,运用整体极大似然估计法对模型参数进行估计,并基于失效机理一致性原理进行模型选择,最终对产品的寿命进行预测。适用于多应力加速寿命试验的加速模型确定和寿命分布确定等领域。
背景技术
通常情况下,对于高可靠长寿命的产品,为了在尽可能短的时间内实现对其寿命的评估,加速寿命试验是常用方法。加速寿命试验就是要在保证试验样品失效机理不变的条件下,将其放在更加严苛的环境中试验,使其在更短时间内出现失效,然后再将失效数据反推回正常的工作环境中,以实现提高试验效率的目的。通常的加速寿命试验仅选取对产品寿命影响最大的单个应力,例如温度、湿度、电应力、机械应力等进行试验。
然而在工程实际中,某些产品存在复杂的自然贮存环境导致的多种环境因素综合作用,由于传统的加速模型未考虑到多种环境应力之间存在的耦合效应,会显著影响寿命预测精度。另一方面基于物理化学的失效机理分析手段并不完善,亟需根据统计数据建立多应力耦合效应加速寿命模型,并通过相应模型选择准则选出最符合加速寿命试验数据的模型,给出该模型下产品在正常应力水平下的可靠寿命。
基于此本发明提出一种考虑耦合效应模型构建和失效机理一致性模型选择的寿命预测方法,综合考虑加速试验的加速模型和寿命分布模型的选择,外推出产品在正常应力下的寿命。
发明内容
(1)本发明的目的:针对应力之间存在耦合效应的多应力加速寿命试验数据,提供一种多应力加速寿命试验耦合效应分析及寿命预测方法,它是一种包含耦合效应多应力加速模型构建、参数估计、模型选择、寿命预测的完整产品寿命评估方法。通过分析应力之间的耦合效应,选择最优加速寿命试验模型,提高寿命预测的准确程度。
(2)技术方案:
本发明需建立如下基本设置:
设置1在加速寿命试验中,产品的失效机理恒定,寿命分布的累积失效函数形式相同,加速应力的不同仅改变模型参数的大小;
设置2产品寿命t服从威布尔分布、对数正态分布、伽马分布中的一种,各分布的概率密度函数分别为:
①威布尔分布
其中,η和β分别为威布尔分布的尺度参数和形状参数;
②对数正态分布
其中,μ和σ分别为对数正态分布的位置参数和尺度参数;
③伽马分布
其中,α和λ分别为伽马分布的形状参数和尺度参数;
本发明提出的方法主要针对含有耦合效应的多应力加速寿命试验数据,根据广义对数线性模型构建加速方程,并运用极大似然估计的方法,对模型的参数进行估计,考虑失效机理一致性原则,将对数加速因子变异系数进行比较,确定加速寿命试验数据符合的加速模型并外推其正常应力下的寿命;
基于上述假设,本发明一种多应力加速寿命试验耦合效应分析及寿命预测方法,通过如下步骤实现:
步骤一:利用广义对数线性模型构建多应力加速方程;
基于广义对数线性模型,构建如下多应力加速方程:
广义对数线性模型包含了温度、湿度、电应力;机械应力等常见加速应力对应的阿伦尼乌斯模型与逆幂律模型;通过这些单应力加速模型即可构建上述广义对数线性模型;
具体步骤为:
I.将单应力加速模型转化为广义对数线性模型中的加速应力函数;对于阿伦尼乌斯模型,加速应力Sk的函数在广义对数线性模型中表示为xk=1/Sk;对于逆幂律模型,在广义对数线性模型中表示为xk=ln(Sk);
II.通过单应力加速应力函数构建耦合效应函数;当加速应力Sp与Sq均服从阿伦尼乌斯模型时,两应力的耦合效应函数为ypq=xpxq=1/(Sp·Sq);当加速应力Sp与Sq均服从逆幂律模型时,两应力的耦合效应函数为ypq=xpxq=lnSp·lnSq;当加速应力Sp服从阿伦尼乌斯模型Sq服从逆幂律模型时,两应力的耦合效应函数为ypq=xpxq=(lnSq)/Sp;
III.按照广义对数线性模型的形式,将加速应力函数与耦合效应函数进行组合,构建备选的加速方程;
步骤二:利用整体极大似然理论对模型参数进行估计;
将步骤一中的备选加速方程和设置二中的三种分布进行组合,构成加速寿命试验模型;通过整体极大似然方法分别求解每个模型的参数估计值;
具体步骤为:
I.将广义对数线性模型与寿命分布模型进行组合,构成加速寿命试验模型;当寿命分布模型为威布尔分布时,其尺度参数η不随应力变化,形状参数β与特征寿命相等,是加速应力的函数;当寿命分布模型为对数正态分布时,其尺度参数σ不随应力变化,位置参数μ与特征寿命相等,是加速应力的函数;当寿命分布模型为伽马分布时,其尺度参数λ不随应力变化,形状参数α与特征寿命相等,是加速应力的函数;
II.假设Θ为加速寿命试验模型待估参数向量;加速寿命试验中应力水平组数为m,每组试验包含ni个样本;加速寿命试验数据为完全数据、定数截尾和定时截尾的一种,则对数似然函数ln(Θ)可表示为:
①完全数据:
②定数截尾数据:
③定时截尾数据:
其中,ti,j为第i组应力组合下第j个样本的失效时间;f(t)为产品寿命的概率密度函数,F(t)为产品寿命的累积失效函数;对于定数截尾数据,ri为第i组试验产品失效数;为第ri个产品的失效时间;对于定时截尾数据,t0为定时截尾的截尾时间;
步骤三:通过中位寿命的相对误差判断参数估计结果的合理性;
基于试验数据得到第i组应力水平下的中位寿命基于参数估计结果,计算第i组应力水平下的中位寿命的估计值ti=F-1(0.5,Θ);接着计算每个模型中位寿命的相对误差最后,根据加速寿命试验情况设定相对误差阈值范围,剔除不符合要求的参数估计结果;
步骤四:基于失效机理一致性原则,开展加速模型和寿命分布模型选择;
计算加速寿命试验中各加速应力下限组合Smin相对于各加速应力上限组合Smax的对数加速因子均值和方差;根据机理一致性原则,计算该对数加速因子的变异系数,并选择变异系数最小的模型组合作为该加速寿命试验的最优模型;
具体步骤为:
I.计算参数估计值的协方差;将参数Θ的极大似然估计值Θ′带入Hessian矩阵,并对得到的Hessian矩阵求逆得到参数的协方差矩阵;Hessian矩阵形式为:
则协方差计算公式为:
cov(Θ)=[G(Θ)]-1 (9)
其中,[·]-1为矩阵的逆;
II.计算对数加速因子均值和方差:
其中,μ(AFm)表示对数加速因子均值;分别表示加速应力下限组合Smin和加速应力上限组合Smax的特征寿命;σ2(AFm)表示对数加速因子的方差;为参数估计值和的协方差;xmin,h和xmax,h分别为加速应力Sh取下限和上限时的加速应力函数值;
III.计算对数加速因子的变异系数cv(AFm):
选择变异系数最小的加速寿命试验模型作为最优模型;
步骤五:外推正常应力下产品的可靠寿命;
基于选择的最优寿命分布模型与加速模型,代入正常应力水平,得到其累积失效函数F(t,Θ),给定可靠度R,计算对应的可靠寿命tR:
其中,F-1(·)为总体累积失效函数的反函数。
在步骤一中所述的阿伦尼乌斯模型,是指:
阿伦尼乌斯通过大量数据分析得到的反应速率与温度、激活能之间的关系;阿伦尼乌斯模型的表达形式如下:
对此模型两边作对数变换,可得:
在步骤一中所述的逆幂律模型,是指:
逆幂律模型通常用于描述电应力对产品特征寿命的影响,其表达式为:
将上述表达式两边同时取对数,将逆幂律模型转换为线性形式:
其中,a0=lnA,a1=-c,均为待定参数;因此,特征寿命的对数与电应力的对数函数呈线性关系。
在步骤四中所述的Hessian矩阵,是指:
f的Hessian矩阵表示为:
(3)优点和功效:本发明是一种多应力加速寿命试验耦合效应分析及寿命预测方法,即一种考虑耦合效应模型构建和失效机理一致性模型选择的寿命预测方法,其优点是:
①本发明针对多应力加速寿命试验数据,利用广义对数线性模型构建的多应力加速方程考虑到了应力之间复杂的耦合效应,使模型更加贴合实际;
②本发明采取整体极大似然估计方法,在方法中对参数的初值要求较低,方法迭代快速简单,可操作性强;
③本发明基于失效机理一致性原则,选择最优的加速寿命试验模型,可以提高寿命预测精度;方法科学,工艺性好,具有广阔推广应用价值。
附图说明
图1是本发明所述方法流程图。
图2是本发明所述加速应力进行加速寿命试验得到寿命试验数据,即“失效数据”示意图。
具体实施方式
下面将结合实施例对本发明做进一步详细说明。
某智能电表产品正常工作环境下温度、湿度、电流应力分别为:296K、45%RH、10A;对其按照表1所示加速应力进行加速寿命试验得到寿命试验数据,即“失效数据”见图2所示;
表1某智能电表产品加速寿命试验应力水平
序号 | 温度(K) | 湿度(RH) | 电流(A) |
试验1 | 353 | 80 | 60 |
试验2 | 328 | 80 | 40 |
试验3 | 328 | 95 | 20 |
试验4 | 343 | 95 | 60 |
试验5 | 343 | 95 | 40 |
根据本说明提出的多应力耦合效应分析与寿命预测方法,利用构建多应力加速模型,运用整体极大似然估计方法进行参数估计,利用对数加速因子变异系数选择最优加速模型,并对可靠寿命进行预测。
本发明一种多应力加速寿命试验耦合效应分析及寿命预测方法,见图1所示,其实施步骤如下:
步骤一:利用广义对数线性模型构建多应力加速方程;
首先根据应力的类型确定温度应力和湿度应力符合阿伦尼乌斯模型,电流应力符合逆幂律模型;再将两类加速模型的对数形式按照广义对数线性模型进行组合,构建加速模型如表2所示:
表2广义对数线性模型构建情况
步骤二:利用整体极大似然理论对模型参数进行估计;
首先假设产品的寿命t服从威布尔分布、对数正态分布、伽马分布,且该试验为完全数据试验则通过极大化对数似然函数,可以分别得到威布尔分布、对数正态分布、伽马分布的参数估计结果,如表3、表4、表5所示:
表3威布尔分布参数估计结果
表4对数正态分布参数估计结果
表5伽马分布参数估计结果
步骤三:通过中位寿命的相对误差判断参数估计结果的合理性;
利用参数估计结果计算中位寿命的相对误差,结果如表6所示:
表6中位寿命最大误差结果
在表6的中位寿命最大误差计算结果中,模型一和模型二的中位寿命误差最大值出现在第三组试验中,此时误差最小值达到了25%以上,超出了阈值范围,因此模型一和模型二不适用于该产品。对于其他模型,中位寿命的误差都在25%以下,需要进行进一步判断该产品适用的模型;
步骤四:基于失效机理一致性原则,开展最优加速模型和寿命分布模型选择;
根据步骤三的中位寿命误差结果,计算模型三、四、五、六、七、八的参数估计结果的协方差矩阵,和对数加速因子变异系数;计算结果如表7所示:
表7对数加速因子变异系数计算结果
根据表7的计算结果,模型三对应的三种寿命分布的对数加速因子变异系数都是比较小的,且当寿命服从对数正态分布时,变异系数达到最小值。此时温度与电流存在耦合效应,符合在温度升高时,电流对产品影响增加对物理意义;
步骤五:外推正常应力下产品的可靠寿命;
由步骤四可知最优模型为寿命服从对数正态分布且加速模型为模型三;故选取模型三作为加速模型、对数正态分布作为该产品寿命分布函数;当可靠度R=0.95时,由对数正态函数可解得,产品的可靠寿命为tR=6.8172年;
所以基于前面所述计算评估得产品在正常应力下的95%的可靠寿命为6.8172年;
结果表明,采用本发明方法可以精确外推产品的可靠寿命,达到预期的目的。
综上所述,本发明涉及一种多应力加速寿命试验寿命预测方法,它是一种考虑耦合效应模型构建和失效机理一致性模型选择的寿命预测方法。它针对多应力加速寿命试验数据,利用广义对数线性模型构建包含耦合效应的加速方程,运用极大似然估计法对模型参数进行估计,并基于失效机理一致性原理进行模型选择,最终对产品的寿命进行预测。该方法的具体步骤是:首先用广义对数线性模型构建多种含有耦合效应的多应力加速模型,然后将加速模型和寿命分布模型斤西瓜组合进行整体极大似然估计,之后根据加速寿命试验数据筛选可行的加速模型、进而基于失效机理一致性原理,计算筛选后加速模型的最大最小应力间的对数加速因子变异系数来确定模型选择结果,计算在正常应力下的可靠寿命。本适用于多应力加速寿命试验的加速模型确定和寿命分布确定等领域,具有较强的操作性。
Claims (4)
1.一种多应力加速寿命试验耦合效应分析及寿命预测方法,需要设置如下:
设置1:在加速寿命试验中,产品的失效机理恒定,寿命分布的累积失效函数形式相同,加速应力的不同仅改变模型参数的大小;
设置2:产品寿命t服从威布尔分布、对数正态分布、伽马分布中的一种,各分布的概率密度函数分别为:
①威布尔分布
其中,η和β分别为威布尔分布的尺度参数和形状参数;
②对数正态分布
其中,μ和σ分别为对数正态分布的位置参数和尺度参数;
③伽马分布
其中,α和λ分别为伽马分布的形状参数和尺度参数;
其特征在于,通过如下步骤实现:
步骤一:利用广义对数线性模型构建多应力加速方程;
基于广义对数线性模型,构建如下多应力加速方程:
广义对数线性模型包含了温度、湿度、电应力、机械应力的加速应力对应的阿伦尼乌斯模型与逆幂律模型;通过这些单应力加速模型即构建上述广义对数线性模型;
具体步骤为:
1.1将单应力加速模型转化为广义对数线性模型中的加速应力函数;对于阿伦尼乌斯模型,加速应力Sk的函数在广义对数线性模型中表示为xk=1/Sk;对于逆幂律模型,在广义对数线性模型中表示为xk=ln(Sk);
1.2通过单应力加速应力函数构建耦合效应函数;当加速应力Sp与Sq均服从阿伦尼乌斯模型时,两应力的耦合效应函数为ypq=xpxq=1/(Sp·Sq);当加速应力Sp与Sq均服从逆幂律模型时,两应力的耦合效应函数为ypq=xpxq=lnSp·lnSq;当加速应力Sp服从阿伦尼乌斯模型Sq服从逆幂律模型时,两应力的耦合效应函数为ypq=xpxq=(lnSq)/Sp;
1.3按照广义对数线性模型的形式,将加速应力函数与耦合效应函数进行组合,构建备选的加速方程;
步骤二:利用整体极大似然理论对模型参数进行估计;
将步骤一中的备选加速方程和设置二中的三种分布进行组合,构成加速寿命试验模型;通过整体极大似然方法分别求解每个模型的参数估计值;
具体步骤为:
2.1将广义对数线性模型与寿命分布模型进行组合,构成加速寿命试验模型;当寿命分布模型为威布尔分布时,其尺度参数η不随应力变化,形状参数β与特征寿命相等,是加速应力的函数;当寿命分布模型为对数正态分布时,其尺度参数σ不随应力变化,位置参数μ与特征寿命相等,是加速应力的函数;当寿命分布模型为伽马分布时,其尺度参数λ不随应力变化,形状参数α与特征寿命相等,是加速应力的函数;
2.2设Θ为加速寿命试验模型待估参数向量;加速寿命试验中应力水平组数为m,每组试验包含ni个样本;加速寿命试验数据为完全数据、定数截尾和定时截尾的一种,则对数似然函数ln(Θ)表示为:
①完全数据:
②定数截尾数据:
③定时截尾数据:
其中,ti,j为第i组应力组合下第j个样本的失效时间;f(t)为产品寿命的概率密度函数,F(t)为产品寿命的累积失效函数;对于定数截尾数据,ri为第i组试验产品失效数;为第ri个产品的失效时间;对于定时截尾数据,t0为定时截尾的截尾时间;
2.3通过最大化对数似然函数lnL(Θ)得到未知模型参数的估计结果,即
步骤三:通过中位寿命的相对误差判断参数估计结果的合理性;
基于试验数据得到第i组应力水平下的中位寿命基于参数估计结果,计算第i组应力水平下的中位寿命的估计值ti=F-1(0.5,Θ);接着计算每个模型中位寿命的相对误差最后,根据加速寿命试验情况设定相对误差阈值范围,剔除不符合要求的参数估计结果;
步骤四:基于失效机理一致性原则,开展加速模型和寿命分布模型选择;
计算加速寿命试验中各加速应力下限组合Smin相对于各加速应力上限组合Smax的对数加速因子均值和方差;根据机理一致性原则,计算该对数加速因子的变异系数,并选择变异系数最小的模型组合作为该加速寿命试验的最优模型;
具体步骤为:
4.1计算参数估计值的协方差;将参数Θ的极大似然估计值Θ′带入Hessian矩阵,并对得到的Hessian矩阵求逆得到参数的协方差矩阵;Hessian矩阵形式为:
则协方差计算公式为:
cov(Θ)=[G(Θ)]-1 (9)
其中,[·]-1为矩阵的逆;
4.2计算对数加速因子均值和方差:
其中,μ(AFm)表示对数加速因子均值;分别表示加速应力下限组合Smin和加速应力上限组合Smax的特征寿命;σ2(AFm)表示对数加速因子的方差;为参数估计值和的协方差;xmin,h和xmax,h分别为加速应力Sh取下限和上限时的加速应力函数值;
4.3计算对数加速因子的变异系数cv(AFm):
选择变异系数最小的加速寿命试验模型作为最优模型;
步骤五:外推正常应力下产品的可靠寿命;
基于选择的最优寿命分布模型与加速模型,代入正常应力水平,得到其累积失效函数F(t,Θ),给定可靠度R,计算对应的可靠寿命tR:
其中,F-1(·)为总体累积失效函数的反函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110085635.6A CN112926144B (zh) | 2021-01-22 | 2021-01-22 | 一种多应力加速寿命试验耦合效应分析及寿命预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110085635.6A CN112926144B (zh) | 2021-01-22 | 2021-01-22 | 一种多应力加速寿命试验耦合效应分析及寿命预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112926144A CN112926144A (zh) | 2021-06-08 |
CN112926144B true CN112926144B (zh) | 2022-08-16 |
Family
ID=76164391
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110085635.6A Active CN112926144B (zh) | 2021-01-22 | 2021-01-22 | 一种多应力加速寿命试验耦合效应分析及寿命预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112926144B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113887009B (zh) * | 2021-06-24 | 2024-06-21 | 北京电子工程总体研究所 | 一种定数截尾试验加速失效机理一致性检验方法 |
CN114088117B (zh) * | 2021-11-30 | 2024-06-25 | 中国兵器工业集团第二一四研究所苏州研发中心 | 一种复杂工况下mems惯性器件可靠性评测方法 |
CN114626248B (zh) * | 2022-03-30 | 2023-04-18 | 北京航空航天大学 | 一种基于多应力加速退化数据的螺旋弹簧可靠性评估方法 |
CN114792053B (zh) * | 2022-04-24 | 2022-10-11 | 哈尔滨工业大学 | 一种基于初值-速率相关退化模型的可靠性评估方法 |
CN114925510B (zh) * | 2022-05-06 | 2022-11-11 | 哈尔滨工业大学 | 一种带有自适应交互作用项的多应力加速模型构建方法 |
CN114818348B (zh) * | 2022-05-06 | 2022-10-11 | 哈尔滨工业大学 | 考虑多应力耦合作用对产品退化影响的可靠性评估方法 |
CN114676602B (zh) * | 2022-05-27 | 2022-08-23 | 浙江大学 | 一种改进的电机绝缘多应力可靠性评估方法及系统 |
CN115452358B (zh) * | 2022-09-22 | 2024-01-26 | 中国核动力研究设计院 | 一种电动阀门动作寿命加速试验方法 |
CN115879610B (zh) * | 2022-11-30 | 2024-03-29 | 中国南方电网有限责任公司超高压输电公司广州局 | 接触器使用期限预测方法、装置、设备和存储介质 |
CN115795928B (zh) * | 2023-02-10 | 2023-05-23 | 中国电子产品可靠性与环境试验研究所((工业和信息化部电子第五研究所)(中国赛宝实验室)) | 基于伽马过程的加速退化试验数据处理方法和装置 |
CN117390767B (zh) * | 2023-10-19 | 2024-04-19 | 中国民航大学 | 航空机载电子设备可靠性评估方法、装置及系统和介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110795887A (zh) * | 2019-10-28 | 2020-02-14 | 中国人民解放军国防科技大学 | 多应力加速寿命试验分析方法和装置 |
CN112069697A (zh) * | 2020-09-23 | 2020-12-11 | 中国人民解放军国防科技大学 | 多应力多失效模式相依竞争条件下加速退化试验统计分析方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103336901B (zh) * | 2013-06-27 | 2016-07-06 | 中国人民解放军国防科学技术大学 | 一种竞争失效相关加速寿命试验统计分析方法 |
-
2021
- 2021-01-22 CN CN202110085635.6A patent/CN112926144B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110795887A (zh) * | 2019-10-28 | 2020-02-14 | 中国人民解放军国防科技大学 | 多应力加速寿命试验分析方法和装置 |
CN112069697A (zh) * | 2020-09-23 | 2020-12-11 | 中国人民解放军国防科技大学 | 多应力多失效模式相依竞争条件下加速退化试验统计分析方法 |
Non-Patent Citations (2)
Title |
---|
An Extended Model for Fatigue Life Prediction and Acceleration Considering Load Frequency Effect;Yikun Cai 等;《IEEE Access》;20180404;21064 - 21074 * |
三参数Weibull分布参数的Bayes估计;汤银才等;《系统科学与数学》;20090115(第01期);111-117页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112926144A (zh) | 2021-06-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112926144B (zh) | 一种多应力加速寿命试验耦合效应分析及寿命预测方法 | |
CN113572206B (zh) | 一种风电出力区间预测方法 | |
CN112784415B (zh) | 一种定数截尾加速寿命试验机理等同性检验与寿命预测方法 | |
CN105930976B (zh) | 基于加权理想点法的节点电压暂降严重程度综合评估方法 | |
CN109752656B (zh) | 一种在多时间尺度下用电动汽车电池soh评估方法 | |
CN107436963B (zh) | 一种基于Copula函数多元退化的O型橡胶密封圈寿命预测方法 | |
CN111414991B (zh) | 一种基于多元回归的气象锋面自动识别方法 | |
CN111859658B (zh) | 一种产品贮存寿命与可靠性评估方法 | |
CN107885928B (zh) | 考虑测量误差的步进应力加速性能退化可靠性分析方法 | |
CN107238765A (zh) | 基于加速性能退化参数的led集成驱动电源可靠性分析方法 | |
CN107463993B (zh) | 基于互信息-核主成分分析-Elman网络的中长期径流预报方法 | |
CN110543932A (zh) | 基于神经网络的空调性能预测方法和装置 | |
CN116593896B (zh) | 一种电池储能系统的状态检测方法、系统及电子设备 | |
CN113805138A (zh) | 一种基于参数有向遍历的智能电表误差估计方法及装置 | |
CN113255900A (zh) | 一种考虑改进谱聚类与Bi-LSTM神经网络的冲击性负荷预测方法 | |
CN115686124B (zh) | 一种基于安全保护的储能电池输出功率自调节系统及方法 | |
CN113484813A (zh) | 多环境应力下的智能电表故障率预估方法及系统 | |
CN114839538A (zh) | 一种提取锂离子电池退化特征用于估计剩余寿命的方法 | |
CN115545333A (zh) | 一种多负荷日类型配电网负荷曲线预测方法 | |
CN112305441A (zh) | 一种集成式聚类下的动力电池健康状态评估方法 | |
CN113991711B (zh) | 一种光伏电站储能系统容量配置方法 | |
CN114970157A (zh) | 电子产品在电压应力作用下的小样本试验寿命预测方法 | |
CN113761755A (zh) | 考虑认知不确定性的温湿度双应力下的加速寿命分析方法 | |
CN117200223A (zh) | 日前电力负荷预测方法和装置 | |
CN110163437B (zh) | 基于DPK-means的日前光伏发电功率预测方法 |
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 |