CN113378368B - 一种基于非线性退化轨迹模型的加速因子评估方法 - Google Patents
一种基于非线性退化轨迹模型的加速因子评估方法 Download PDFInfo
- Publication number
- CN113378368B CN113378368B CN202110621941.7A CN202110621941A CN113378368B CN 113378368 B CN113378368 B CN 113378368B CN 202110621941 A CN202110621941 A CN 202110621941A CN 113378368 B CN113378368 B CN 113378368B
- Authority
- CN
- China
- Prior art keywords
- performance
- model
- degradation
- temperature
- rate constant
- 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
- 230000001133 acceleration Effects 0.000 title claims abstract description 58
- 238000006731 degradation reaction Methods 0.000 title claims abstract description 53
- 230000015556 catabolic process Effects 0.000 title claims abstract description 50
- 238000011156 evaluation Methods 0.000 title claims abstract description 15
- 238000012360 testing method Methods 0.000 claims abstract description 44
- 238000000034 method Methods 0.000 claims abstract description 23
- 230000004913 activation Effects 0.000 claims abstract description 13
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000000137 annealing Methods 0.000 claims 1
- 238000006243 chemical reaction Methods 0.000 description 10
- 238000004458 analytical method Methods 0.000 description 3
- 238000010200 validation analysis Methods 0.000 description 3
- 238000012795 verification Methods 0.000 description 3
- 230000007850 degeneration Effects 0.000 description 2
- 238000012417 linear regression Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- SYHGEUNFJIGTRX-UHFFFAOYSA-N methylenedioxypyrovalerone Chemical compound C=1C=C2OCOC2=CC=1C(=O)C(CCC)N1CCCC1 SYHGEUNFJIGTRX-UHFFFAOYSA-N 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 239000003380 propellant Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000010998 test method 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/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/08—Thermal analysis or thermal optimisation
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 Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)
Abstract
本发明涉及一种基于非线性退化轨迹模型的加速因子评估方法,其包括如下步骤:(1)收集加速试验数据;(2)建立应力与性能变化速率常数的关系;(3)表观活化能及性能退化曲线拟合;(4)计算加速因子。本发明提出产品的活化能作为产品的固有属性,会随性能退化过程而改变。采用点斜法开展加速寿命试验时,需根据产品性能退化量选取合适的加速因子,以验证产品是否满足贮存使用要求。
Description
技术领域
本发明属于导弹寿命试验与评估技术领域,具体涉及一种基于非线性退化轨迹模型的加速因子评估方法。
背景技术
针对导弹产品贮存寿命加速试验,由于产品价格昂贵,依靠常规的不同类型的加速试验方法,试验件满足不了要求,成本太高,因此,一般采用贮存加速试验是验证性试验。其基本思想是通过确定产品加速因子,利用加速因子将产品的实际贮存时间进行折合得到加速试验时间,以验证产品是否满足贮存使用要求。可以看出,上述加速贮存验证试验中,通过加速因子进行时间折算,求取整机加速因子是关键。
然而,上述方法有一个基本假设,就是产品的加速因子不变,然而实际贮存时,产品的性能退化曲线往往是非线性的,导致其加速因子可能会随着贮存时间而变化,应用上述方法会造成一些理解上的误导,不能指导加速退化试验加速因子的评估工作。
发明内容
本发明的目的是提供一种基于非线性退化轨迹模型的加速因子评估方法。
本发明采用如下技术方案:
一种基于非线性退化轨迹模型的加速因子评估方法,其包括如下步骤:
(1)收集加速试验数据;
(2)建立应力与性能变化速率常数的关系;
(3)表观活化能及性能退化曲线拟合;
(4)计算加速因子。
步骤(1)中,根据加速试验的结果,建立性能变化速率常数K(T)与时间t的关系式;对每个加速试验温度可获得一组加速时间t与P的数据:
步骤(2)中,参照指数模型或对数模型计算每个加速试验温度的性能变化速率常数K;根据计算结果,建立性能变化速率常数K与温度T的关系。
其中,性能变化速率常数K通过如下方法确定:
(a)将指数模型表示为Y=a+bX;用最小二乘法求得系数a、b和相关系数r;
(b)查相关系数表以置信度为99%;自由度f=n-2的r值计算的r值比较,若r计算值的绝对值大于r的查表值,则X与Y线性关系成立,可用Y=a+bX表示,方程的斜率b为相应加速试验温度下的性能变化速率常数K;若r计算值小于r的查表值,则X与Y线性关系不成立;
(c)若X与Y线性关系不成立,则将指数模型换用对数模型重复步骤(a)和(b)。
其中,通过如下方法建立性能变化速率常数K与温度T的关系:
Arrhenius方程式用Y1=a1+b1X1表示,用最小二乘法求系数a1、b1和相关系数r1;
查相关系数表以置信度为95%,自由度f=m-2的r值与r1计算值比较,如果r1计算值大于r查表值,则线性关系成立,可用Y1=a1+b1X1表示;反之,则线性关系不成立。如不成立,则本发明不适用。
其中,Arrhenius方程式用Y1=a1+b1X1表示时,其中X1=1/T,Y1=lnK,a1=ln(Z/Pk),b1=-Ea/R。
步骤(3)具体包括:
若指数模型表示为Y=a+bX后,X与Y线性关系成立,那么对于给定的性能退化量Pk通过下式进行计算:
Pk=P0+K1lnt1=P0+Kslnts
其中,K1-实际贮存温度为t1℃的性能变化速率常数;
Ks-加速贮存温度为ts℃的性能变化速率常数。
步骤(4)具体包括:
若性能参数退化轨迹服从指数模型,对于给定的加速温度应力T1,,加速因子AF可表述为:
若性能参数退化轨迹服从对数模型,对于给定的加速温度应力T1,,加速因子AF可表述为:
本发明的有益效果在于:
(1)加速因子定义为产品在加速应力水平与基准贮存应力水平作用下某一贮存寿命特征的比值,对于非线性退化轨迹模型,其化学反应速率v(T)随性能退化量的改变而改变,不能简单的用性能变化速率常数K(T)来代替化学反应速率v(T)。
根据Arrhenius方程的定义,推导提出用化学化学反应速率v(T)来表征加速模型,可以避免现有方法中,将性能变化速率常数K(T)用Arrhenius方程来表述,缺少明确的物理意义的问题。
(2)提出产品的活化能作为产品的固有属性,会随性能退化过程而改变。采用点斜法开展加速寿命试验时,需根据产品性能退化量选取合适的加速因子。
若性能参数退化轨迹服从指数模型,其加速因子跟线性模型一样,不随性能退化量的改变而改变;若性能参数退化轨迹服从对数模型,其加速因子随性能变化量的改变而改变,即随着加速贮存试验时间t1延长而改变,在开展贮存期验证试验时需特别注意这一点。
附图说明
图1为整机Ⅲ舵输出电流(-)性能参数退化曲线图。
图2为90℃下,ψ波道输入电压0.9950V时,Ⅲ舵输出电流(-)数据拟合图。
图3为100℃下,ψ波道输入电压0.9950V时,Ⅲ舵输出电流(-)数据拟合图。
图4为110℃下,ψ波道输入电压0.9950V时,Ⅲ舵输出电流(-)数据拟合图。
图5为120℃下,ψ波道输入电压0.9950V时,Ⅲ舵输出电流(-)数据拟合图。
图6为加速模型拟合图。
具体实施方式
一、公知技术
对于非线性退化轨迹模型的加速寿命试验模型,主要采用QJ2328A-2005标准规定的方法,这些方法建立在两个假设之上:
假设1:对于大多数具有软故障(即性能退化)的产品来说,尤其是材料/结构类产品(例推进剂等),产品的性能参数指标与试验时间可用以下几种数学模型来进行有效的拟合:
线性模型:P=P0+K·t [1]
指数模型:P=P0·e-Kt [2]
对数模型:P=P0+K·lnt [3]
式[1]~[3]中:P为产品性能参数指标;P0为常数,性能参数初始值;K为与温度有关的性能变化速率常数,d-1;t为试验时间,d。
假设2:产品的性能变化速率常数K与温度T的关系服从Arrhenius方程:
式[4]中:T为绝对温度,K;Ea为表观活化能,J·mol-1;Z为频率因子,d-1;R为气体常数,J·K-1·mol-1。
二、公知技术的缺陷分析
当前贮存期验证试验中,往往认为在整个试验过程中加速因子保持不变,即对于给定的加速应力水平T1,其相对于实际贮存使用温度应力Ts的加速因子是常数,即:
式[5]中:AF为加速因子;K(T1)为加速试验条件下性能变化速率常数,d-1;K(Ts)为实际贮存使用条件下性能变化速率常数,d-1;Ts为实际贮存使用温度应力,K;T1为加速温度应力,K。
但是,这与实际情况可能不符,事实上加速因子一般的定义为产品在加速应力水平与基准贮存应力水平作用下某一贮存寿命特征的比值,即时间之比或速率之比:
针对式[1]~[3]对应的退化轨迹模型,其化学反应速率分别为:
可以看出,对于线性退化轨迹模型,产品的化学反应速率v(T)即为性能变化速率常数K(T)。但是,对于非线性退化轨迹模型,产品的化学反应速率v(T)不等同于性能变化速率常数K(T)。
实际上,对于非线性退化轨迹模型,不能简单的用性能变化速率常数K(T)来代替化学反应速率v(T),产品的化学反应速率v(T)随性能退化量的改变而改变,如果简单的套用线性退化轨迹模型的试验数据处理方法,会造成一些理解上的误导,不能指导加速退化试验加速因子的评估工作。
三、本发明提出的技术方案
本发明严格按照Arrhenius方程的定义,通过严格推导化学反应速率v(T)的表达式,给出基于非线性退化轨迹模型的加速寿命试验评估方法。
1.收集加速试验数据
根据加速试验的结果,建立性能变化速率常数K(T)与时间t的关系式。对每个加速试验温度可获得一组(n个)加速时间t与P的数据。
2.建立应力与性能变化速率的关系
2.1指数模型
参照式[2]计算每个加速试验温度的性能变化速率常数K。令X=t;Y=lnP;a=lnP0;b=-K。则[2]式可用Y=a+bX表示。用最小二乘法求得系数a、b和相关系数r。
其中:
查QJ2328A-2005标准中的相关系数表,以置信度为99%;自由度f=n-2的r值计算的r值比较,若r计算值的绝对值大于r的查表值,则X与Y线性关系成立,可用Y=a+bX表示,方程的斜率b为相应加速试验温度下的性能变化速率常数。若r计算值小于r的查表值,则X与Y线性关系不成立,对于非线性退化轨迹模型,可尝试使用2.2对数模型建立应力与性能变化速率的关系。
2.2对数模型
参照式[3]计算每个加速试验温度的性能变化速率常数K。令X=lnt;Y=P;a=P0;b=K。则式[3]可用Y=a+bX表示。用最小二乘法求得系数a、b和相关系数r,方法见式[11]~式[18]。
查QJ2328A-2005标准中的相关系数表,以置信度为99%;自由度f=n-2的r值计算的r值比较,若r计算值的绝对值大于r的查表值,则X与Y线性关系成立,可用Y=a+bX表示,方程的斜率b为相应加速试验温度下的性能变化速率常数。若r计算值小于r的查表值,则X与Y线性关系不成立。
若指数模型和对数模型均无法使X与Y线性关系成立,则不符合假设1中的非线性退化轨迹模型,考虑使用其他模型进行分析。
2.3应力与性能变化速率的关系
根据以上计算结果,确定非线性退化轨迹模型,建立性能变化速率常数K与温度T的关系。
令X1=1/T,Y1=lnK,a1=ln(Z/Pk),b1=-Ea/R。则式[4]可以用Y1=a1+b1X1表示,用最小二乘法求系数a1、b1和相关系数r1,方法见式[11]~式[18]。
查QJ2328A-2005标准中的相关系数表,以置信度为95%,自由度f=m-2的r值与r1计算值比较,如果r1计算值大于r查表值,则线性关系成立,可用Y1=a1+b1X1表示。反之,则线性关系不成立。
若线性关系不成立,则产品的性能变化速率常数K与温度T的关系不服从Arrhenius方程(不符合假设2)。
3.激活能及性能退化曲线拟合
3.1指数模型
按照线性方程Y1=a1+b1X1计算,表观活化能(激活能)为见式[19],对于给定的性能退化量Pk有式[20]成立。
Ea=-b1·R [19]
其中:
Ks-实际贮存温度为ts℃的性能变化速率常数;
K1-加速贮存温度为t1℃的性能变化速率常数。
3.2对数模型
按照线性方程Y1=a1+b1X1计算,表观活化能(激活能)为见式[19],对于给定的性能退化量Pk有式[21]成立。
Pk=P0+K1lnt1=P0+Kslnts [21]
其中:
Ks-实际贮存温度为ts℃的性能变化速率常数;
K1-加速贮存温度为t1℃的性能变化速率常数。
4.计算加速因子
4.1指数模型
若性能参数退化轨迹服从指数模型,对于给定的加速温度应力T1,,加速因子AF可表述为:
可以看出,若性能参数退化轨迹服从指数模型,其加速因子评估不随性能退化过程而改变。
4.2对数模型
若性能参数退化轨迹服从对数模型,对于给定的加速温度应力T1,,加速因子AF可表述为:
可以看出,若性能参数退化轨迹服从对数模型,其加速因子评估随性能退化过程而改变,即随着加速贮存试验时间t1延长而改变,在开展贮存期验证试验时需特别注意这一点。
四、具体计算实例
某整机产品实际贮存温度为20℃,本次基于4个温度应力水平开展了12个整机产品的加速贮存试验,温度应力水平分别为90℃、100℃、110℃和120℃,得到了加速贮存试验数据,其中以“ψ波道输入电压0.9950V时,Ⅲ舵输出电流(-)”作为典型参数进行加速因子评估案例分析。所得整机性能参数退化曲线如图1所示。
对性能参数退化数据进行线性相关性判断,如表1所示。
表1不同温度下不同产品加速试验数据的相关性
产品编号 | 产品1 | 产品2 | 产品3 | 结果 |
90℃相关系数判定 | 0.9185 | 0.9029 | 0.9334 | 均具有相关性 |
100℃相关系数判定 | 0.8498 | 0.8574 | 0.8770 | 均具有相关性 |
110℃相关系数判定 | 0.9539 | 0.9173 | 0.8196 | 均具有相关性 |
120℃相关系数判定 | 0.8560 | 0.8781 | 0.9348 | 均具有相关性 |
由上表可知,“ψ波道输入电压0.9950V时,Ⅲ舵输出电流(-)”参数退化轨迹服从指数退化轨迹模型。曲线拟合情况如图2~图5所示。
不同温度点下参数估计值如表2所示。
表2参数估计结果
利用相关系数r检验来验证Y1=a1X1+b1方程的线性相关性。通过相关系数临界值表可知,显著性水平为0.1,自由度f=4-2=2的相关系数临界值为r0=0.900。而通过计算得到方程Y1=a1X1+b1的线性相关系数r=0.9160,高于相关系数临界值,因此可以得出线性回归效果显著(如图6所示)。
通过计算得到a1=-3.4481×103,进而计算获得整机产品激活能Ea=-a1*k=0.2971(eV)。
表3线性化方程参数
序号 | 温度T(K) | Y<sub>1</sub>=ln|a| | X<sub>1</sub>=1/T |
1 | 90+273.16 | -3.7781 | 0.0028 |
2 | 100+273.16 | -3.8134 | 0.0027 |
3 | 110+273.16 | -3.2846 | 0.0026 |
4 | 120+273.16 | -3.1458 | 0.0025 |
则不同温度对应的加速因子如表4所示。
表4加速因子
Claims (4)
1.一种基于非线性退化轨迹模型的加速因子评估方法,其特征在于,其包括如下步骤:
(1)收集加速试验数据;根据加速试验的结果,建立性能变化速率常数K(T)与时间t的关系式;对每个加速试验温度可获得一组加速时间t与P的数据:
P为产品性能参数指标;
(2)建立应力与性能变化速率常数的关系;
参照指数模型或对数模型计算每个加速试验温度的性能变化速率常数K;根据计算结果,建立性能变化速率常数K与温度T的关系;
所述性能变化速率常数K通过如下方法确定:
(a)将指数模型表示为Y=a+bX;用最小二乘法求得系数a、b和相关系数r;
(b)查相关系数表以置信度为99%;自由度f=n-2的r值计算的r值比较,若r计算值的绝对值大于r的查表值,则X与Y线性关系成立,可用Y=a+bX表示,方程的斜率b为相应加速试验温度下的性能变化速率常数K;若r计算值小于r的查表值,则X与Y线性关系不成立;
(c)若X与Y线性关系不成立,则将指数模型换用对数模型重复步骤(a)和(b);
(3)表观活化能及性能退化曲线拟合;
(4)计算加速因子;
若性能参数退化轨迹服从指数模型,对于给定的加速温度应力T1,Ts为实际贮存使用温度应力,加速因子AF可表述为:
若性能参数退化轨迹服从对数模型,对于给定的加速温度应力T1,加速因子AF可表述为:
3.根据权利要求2所述的基于非线性退化轨迹模型的加速因子评估方法,其特征在于,Arrhenius方程式用Y1=a1+b1X1表示时,其中X1=1/T,Y1=lnK,a1=ln(Z/Pk),b1=-Ea/R;T为绝对温度;Ea为表观活化能;Z为频率因子;R为气体常数;Pk为性能退化量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110621941.7A CN113378368B (zh) | 2021-06-03 | 2021-06-03 | 一种基于非线性退化轨迹模型的加速因子评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110621941.7A CN113378368B (zh) | 2021-06-03 | 2021-06-03 | 一种基于非线性退化轨迹模型的加速因子评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113378368A CN113378368A (zh) | 2021-09-10 |
CN113378368B true CN113378368B (zh) | 2022-07-22 |
Family
ID=77575794
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110621941.7A Active CN113378368B (zh) | 2021-06-03 | 2021-06-03 | 一种基于非线性退化轨迹模型的加速因子评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113378368B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115659642B (zh) * | 2022-10-26 | 2023-07-18 | 兰州空间技术物理研究所 | 一种离子电推进寿命试验加速因子确定方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103777098A (zh) * | 2014-01-23 | 2014-05-07 | 北京航空航天大学 | 一种基于低频噪声分类的光耦贮存寿命评价方法 |
CN104236889A (zh) * | 2014-09-11 | 2014-12-24 | 航天科工防御技术研究试验中心 | 橡胶密封圈的贮存寿命优化设计方法 |
CN105468907A (zh) * | 2015-11-23 | 2016-04-06 | 北京航空航天大学 | 一种加速退化数据有效性检验及模型选择方法 |
CN108399278A (zh) * | 2018-01-24 | 2018-08-14 | 航天科工防御技术研究试验中心 | 一种电子装备多因素加速因子计算方法 |
CN111678535A (zh) * | 2020-04-27 | 2020-09-18 | 中国电子产品可靠性与环境试验研究所((工业和信息化部电子第五研究所)(中国赛宝实验室)) | 光纤陀螺可靠度加速试验方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130138419A1 (en) * | 2011-11-30 | 2013-05-30 | Oracle International Corporation | Method and system for the assessment of computer system reliability using quantitative cumulative stress metrics |
JP6374447B2 (ja) * | 2016-07-08 | 2018-08-15 | ファナック株式会社 | 温度を含む駆動条件を考慮した実効的駆動時間と残存寿命が算出可能なレーザ装置 |
CN107015875B (zh) * | 2017-03-31 | 2020-09-04 | 北京强度环境研究所 | 一种电子整机贮存寿命评估方法及装置 |
US11204391B2 (en) * | 2019-09-13 | 2021-12-21 | GM Global Technology Operations LLC | Method and apparatus for monitoring a battery state estimator |
CN110851980B (zh) * | 2019-11-11 | 2021-01-29 | 中国人民解放***箭军工程大学 | 一种设备剩余寿命预测方法及*** |
CN111859658B (zh) * | 2020-07-15 | 2023-06-27 | 北京强度环境研究所 | 一种产品贮存寿命与可靠性评估方法 |
-
2021
- 2021-06-03 CN CN202110621941.7A patent/CN113378368B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103777098A (zh) * | 2014-01-23 | 2014-05-07 | 北京航空航天大学 | 一种基于低频噪声分类的光耦贮存寿命评价方法 |
CN104236889A (zh) * | 2014-09-11 | 2014-12-24 | 航天科工防御技术研究试验中心 | 橡胶密封圈的贮存寿命优化设计方法 |
CN105468907A (zh) * | 2015-11-23 | 2016-04-06 | 北京航空航天大学 | 一种加速退化数据有效性检验及模型选择方法 |
CN108399278A (zh) * | 2018-01-24 | 2018-08-14 | 航天科工防御技术研究试验中心 | 一种电子装备多因素加速因子计算方法 |
CN111678535A (zh) * | 2020-04-27 | 2020-09-18 | 中国电子产品可靠性与环境试验研究所((工业和信息化部电子第五研究所)(中国赛宝实验室)) | 光纤陀螺可靠度加速试验方法 |
Non-Patent Citations (3)
Title |
---|
Acceleration factor for ageing measurement of dye solar cells;Laura Ciammaruchi,et al;《Microelectronics Reliability》;20130916;正文第279-281页 * |
Analysis On The Non- Arrhenius Life Prediction Method Of Rubber;Wang You-shan,et al;《Advanced Materials Research》;20130424;正文第366-371页 * |
机电产品加速贮存试验与寿命评估方法研究;张文广 等;《机电工程》;20210531;正文第528-535页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113378368A (zh) | 2021-09-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109977491B (zh) | 一种冲击损伤可恢复条件下的退化建模与寿命预测方法 | |
CN110851980B (zh) | 一种设备剩余寿命预测方法及*** | |
CN101976311B (zh) | 基于漂移布朗运动模型的加速退化试验贝叶斯评估方法 | |
CN105740988A (zh) | 基于灰色关联分析和多元线性回归模型对煤炭热值的预测方法 | |
CN113378368B (zh) | 一种基于非线性退化轨迹模型的加速因子评估方法 | |
CN110110476B (zh) | 基于威布尔分布及右截尾数据估计可靠度的样本取舍方法 | |
CN108918261B (zh) | 一种用少量试验测定材料构件疲劳寿命规律的方法 | |
CN111859658A (zh) | 一种产品贮存寿命与可靠性评估方法 | |
CN107220500B (zh) | 基于逆高斯过程的性能退化试验贝叶斯可靠性评估方法 | |
CN114970157B (zh) | 电子产品在电压应力作用下的小样本试验寿命预测方法 | |
CN103258115A (zh) | 基于多项式拟合和寿命分布的锂离子电池动态可靠性评估方法 | |
CN113011012A (zh) | 一种基于Box-Cox变化的储能电池剩余寿命预测方法 | |
CN110414086B (zh) | 一种基于灵敏度的综合应力加速因子计算方法 | |
CN114295491B (zh) | 一种蠕变损伤及变形随时间演化行为的预测方法 | |
CN113791351B (zh) | 基于迁移学习和差值概率分布的锂电池寿命预测方法 | |
Mehdi et al. | Determination of plotting position formula for the normal, log-normal, Pearson (III), log-Pearson (III) and Gumble distributional hypotheses using the probability plot correlation coefficient test | |
CN112765769B (zh) | 一种多阶段变工况固体推进剂剩余贮存寿命预测方法 | |
CN114169128A (zh) | 一种基于Bayes分析的可靠性强化试验定量评估方法 | |
CN112070030B (zh) | 一种巴克豪森信号随机性度量及转换方法 | |
CN110135002B (zh) | 一种衡量新能源汽车电池充电累积容量损耗可靠性的方法 | |
CN110263472B (zh) | 基于回归法综合寿命试验数据的机电产品可靠度评估方法 | |
CN109359375B (zh) | 一种数据产品贮存寿命预测方法 | |
Tanis et al. | Comparisons of six different estimation methods for log-Kumaraswamy distribution | |
CN114136538B (zh) | 一种基于随机变分贝叶斯学习的压力传感器标定装置温度建模方法 | |
CN114792053A (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 |