CN109657937A - 一种基于退化数据的产品可靠性评估与寿命预测方法 - Google Patents
一种基于退化数据的产品可靠性评估与寿命预测方法 Download PDFInfo
- Publication number
- CN109657937A CN109657937A CN201811455684.9A CN201811455684A CN109657937A CN 109657937 A CN109657937 A CN 109657937A CN 201811455684 A CN201811455684 A CN 201811455684A CN 109657937 A CN109657937 A CN 109657937A
- Authority
- CN
- China
- Prior art keywords
- model
- degradation
- stress
- reliability assessment
- data
- 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 59
- 230000015556 catabolic process Effects 0.000 claims abstract description 60
- 238000006731 degradation reaction Methods 0.000 claims abstract description 60
- 230000001133 acceleration Effects 0.000 claims abstract description 29
- 238000007476 Maximum Likelihood Methods 0.000 claims abstract description 18
- 238000010276 construction Methods 0.000 claims abstract description 11
- 238000009826 distribution Methods 0.000 claims abstract description 7
- 230000008080 stochastic effect Effects 0.000 claims abstract description 3
- 230000006870 function Effects 0.000 claims description 17
- 230000005653 Brownian motion process Effects 0.000 claims description 14
- 238000012360 testing method Methods 0.000 claims description 14
- 238000003860 storage Methods 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 5
- 238000009795 derivation Methods 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 3
- 230000003993 interaction Effects 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 claims description 2
- 238000013461 design Methods 0.000 abstract description 3
- 230000035882 stress Effects 0.000 description 53
- 230000008569 process Effects 0.000 description 6
- 238000011160 research Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 4
- 230000014759 maintenance of location Effects 0.000 description 4
- ZVQXQPNJHRNGID-UHFFFAOYSA-N tetramethylsuccinonitrile Chemical compound N#CC(C)(C)C(C)(C)C#N ZVQXQPNJHRNGID-UHFFFAOYSA-N 0.000 description 4
- 238000009792 diffusion process Methods 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 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 description 2
- 230000032683 aging Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000007850 degeneration Effects 0.000 description 2
- 230000003412 degenerative effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 150000001336 alkenes Chemical class 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Development Economics (AREA)
- Strategic Management (AREA)
- Educational Administration (AREA)
- Entrepreneurship & Innovation (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Marketing (AREA)
- Game Theory and Decision Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明属于工程设计与数据建模技术领域，公开了一种基于退化数据的产品可靠性评估方法、参数估计方法；建立基于随机效应的Wiener加速退化模型，其中漂移参数为服从正态分布的随机变量；根据加速退化模型构建剩余寿命预测模型，且推导出剩余寿命分布及剩余寿命的可靠度函数；根据应力状态下的性能退化数据，通过两步极大似然估计法，获得了模型未知参数估计值。本发明通过经典的Wiener模型构建了一种新的加速退化模型、剩余寿命预测模型，求得剩余寿命的估计值；采用两步极大似然估计法，估计了退化模型中的未知参数，弥补传统极大似然估计法的限制性，能更为精确地得到未知参数的最优值，模型能更准确的预测产品的剩余寿命，适用范围较广。
Description
技术领域
本发明属于工程设计于数据建模技术领域，尤其涉及一种基于退化数据的产品可靠性评估方法、模型参数估计方法。
背景技术
目前，业内常用的现有技术是这样的：随着科学技术的不断进步和产品加工技艺的不断提高，越来越多的产品具有高可靠、长寿命的特点。但在相对短期内高长寿命产品几乎不可能失效，甚至会出现“零失效”现象，由此基于失效数据很难对产品可靠性进行评估，则预测其剩余寿命也增加了难度。加速退化试验(Accelerated Degradation Testing，ADT)是通过提高试验应力水平来使产品的性能加速退化，并利用搜集到的产品在高应力平下的性能退化信息来预测产品在正常工作应力水平下产品的可靠性。因此，加速退化试验成为了评估产品可靠性的重要技术之一。加速退化数据是指在加速退退化实验的过程中，通过人为施加应力造成产品性能退化时测量的数据，或者是产品在实际工作中受到外部环境因素(温度、湿度等)的影响导致产品的加速退化获得的退化数据。有很多研究工作是针对加速退化数据的可靠性评估。在这些研究工作中，基于Wiener过程的加速退化模型在很多领域有较为广泛的应用，电动机加速寿命研究、光伏组件可靠性评估、集成逻辑系统的可靠性评估等。Wiener过程的加速退化模型应用于“感应电动机加速退化试验”中时，模型的漂移系数与增加的载荷(可以看成一种应力)有关系，而其扩散系数与增加的载荷无关系。Wiener过程的加速退化模型应用于研究“随机参数逆高斯过程的加速退化建模方法”中时，漂移系数和扩散系数都受到综合应力水平的影响。Wiener过程的加速退化模型应用于研究“光伏组件的可靠性的加速寿命试验”时，通过引入应力并在达到阈值时确定寿命分布。Wiener过程的加速退化模型应用于研究“集成逻辑系统可靠性”时，最大退化量和应力之间的关系用Arrhenius关系描述。由此可以看到，Wiener过程的加速退化模型用于实际系统的可靠性或寿命分布研究时。模型中的漂移系数和扩散参数的随机性、相互关系等是建立模型的关键之一。另外，基于数据的模型参数的估计方法也是可靠性评估和寿命预测关键技术之一。
综上所述，现有技术存在的问题是：
(1)基于Wiener过程的退化模型中，模型中的两个重要参数的假设有些理想化，与实际数据情况不相符合。
(2)基于Wiener过程的加速退化模型中，性能退化量与应力关系的描述较为单一，应力与模型参数关系如何建模仍然是技术难题。
(3)基于Wiener过程的性能退化建模，需要设计更为可行的参数估计方法和算法，参数估计量的优良性需要进一步分析。
解决上述技术问题的难度和意义：难度是如何确立基于Wiener过程加速退化建模中参数的随机性及其参数之间关系的描述；如何建立模型中加速应力和参数的约束关系；如何给出更为精确的模型参数估计算法。解决上述问题的意义在于能提供一种确立模型参数的数据处理方法；建立基于退化数据的更加准确和方便系统可靠性评估方法；提供有效的模型参数估计方法，为基于信息处理的系统可靠性评估提供一种可行的技术方案。
发明内容
针对现有技术存在的问题，本发明提供了一种基于性能退化数据的产品可靠性评估方法、模型参数估计方法。
本发明是这样实现的，一种产品可靠性评估方法，所述产品可靠性评估方法建立了一种基于随机效应的Wiener加速退化模型，其中漂移参数为服从正态分布的随机变量；根据加速退化模型构建剩余寿命预测模型，且推导出剩余寿命分布及剩余寿命的可靠度函数；根据应力状态下的性能退化数据，通过两步极大似然估计法，获得了模型未知参数估计值。
进一步，所述产品可靠性评估方法包括以下步骤：
步骤一，收集性能加速退化数据，产品试样个数为m个，施加的温度应力是Sk(k＝1,2,…,m)；在这m个温度应力水平下，分别对产品试样进行n次性能退化数据监测、记录；表示第i次对产品性能的测量时刻，产品所增加的温度应力为Sk，测得的产品退化量为i＝1,2,…,n；k＝1,2,…,m；
步骤二，构建经典的加速退化模型，加速退化试验中，加速关系常指的是退化率与压力之间的相互作用，加速模型为：
其中a、b为系数，表示以应力S为自变量函数，采用Arrhenius模型表示，即温度应力下加速模型为：
μk＝aexp(-b/Sk)；
式中Sk为第k个应力；μk为漂移系数，
基于Wiener过程的应力加速退化数据建模：
X(t)＝μkt+σB(t)；
其中μk～N(μa exp(-b/Sk),
步骤三，特定产品在时刻tl和t+tl的性能退化量为X(tl)和X(t+tl)，有：
X(t+tl)＝X(tl)+μkt+σB(t)；
其中μk～N(μa exp(-b/Sk),
当前的退化量X(tk)和失效阈值ξ，产品在时刻tk的剩余寿命t表示为：
则产品的剩余寿命t的PDF表示为：
根据剩余寿命t的定义，得到产品可靠度函数为：
步骤四，采用两步极大似然估计方法估计对数极大似然估计函数如下：
式中：是未知参数的集合；lnL(Θ)表示利用fminsearch函数所求解得到的最大值；
ak、的偏导数为零：
则：
根据可得到的数据组及本发明所建模型，模型中的未知的参数可利用两步极大似然估计法来求解，先估计再估计
步骤五：得到加速应力下基于Wiener模型的剩余寿命预测的估计值；采用两步极大似然估计方法得到
本发明的另一目的在于提供一种实施所述产品可靠性评估方法的产品可靠性评估系统，所述产品可靠性评估系统包括：
数据获取模块，用于获取加速退化数据；
模型构建模块，用于构建典型的加速退化模型；
构建及推导模块，用于基于Wiener模型构建加速应力下的剩余寿命预测模型以及推导出函数表达式；
参数估计模块，用于结合性能加速退化数据的特征给出模型未知参数的估计方法；
预测模块，用于得到加速应力下基于Wiener模型的剩余寿命预测的估计值。
本发明的另一目的在于提供一种实现所述产品可靠性评估方法的计算机程序。
本发明的另一目的在于提供一种实现所述产品可靠性评估方法的信息数据处理终端。
本发明的另一目的在于提供一种计算机可读存储介质，包括指令，当其在计算机上运行时，使得计算机执行所述的产品可靠性评估方法。
综上所述，本发明的优点及积极效果为：本发明通过构建经典的加速退化模型，以及构建基于Wiener模型加速应力下的剩余寿命预测模型，并求得剩余寿命的估计值。并采用两步极大似然估计法，对Wiener过程退化模型中未知参数进行估计，可以更好地弥补传统极大似然估计法的限制性，更加精确地求解未知参数的最优值，说明该模型适用范围较广，且更加准确的预测了产品的剩余寿命。
表2是在四种温度应力下分别求出的模型未知参数的结果。从图3中可以知道，和都随着135℃、150℃、165℃、180℃温度的增加依次增大，相对和来说变化不大，说明Wiener退化模型的和都受温度应力的影响，而所受影响较小。
表3是依照建立的Wiener退化模型和已有的两种模型，分别求出的模型未知参数μ、σ2以及产品寿命置信度为95％的区间估计CI值，三种模型的AIC值，三种模型的TMSN值。由图可以看出，对于三种模型的各种参数值的对比，本发明中的模型的CI值，AIC值，TMSN值都是最小的。说明本发明中的模型比比其它两种模型的拟合性要好，也说明本发明的模型精确度高，更加的符合实际情况。
表4是分别取1000h、1500h和2000h的时间点时对电缆剩余寿命预测，得到的剩余寿命均值。本发明的模型在每个时刻预测的剩余寿命结果都是最接近真实值的，反映了本发明的模型比较好，更贴近样本的实际退化情况。另外两种模型与真实值相差比较大，这是因为没有考虑应力以及个体差异性对退化过程的影响。剩余寿命预测的结果也就与实际剩余寿命有偏差。
附图说明
图1是本发明实施例提供的产品可靠性评估系统结构示意图；
图中：1、数据获取模块；2、模型构建模块；3、构建及推导模块；4、参数估计模块；5、预测模块。
图2是本发明实施例提供的产品可靠性评估方法流程图。
图3是本发明实施例提供的产品可靠性评估方法实现流程图。
图4是本发明实施例提供的四种应力下伸长率保留率与时间的关系曲线图；
图中：(a)135℃时裂伸长率保留率与时间的关系曲线；(b)150℃时裂伸长率保留率与时间的关系曲线；(c)165℃时裂伸长率保留率与时间的关系曲线；(d)180℃时裂伸长率保留率与时间的关系曲线。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
本发明通过构建经典的加速退化模型，以及构建基于Wiener模型加速应力下的剩余寿命预测模型，并求得剩余寿命的估计值。
下面结合附图对本发明的应用原理作详细的描述。
如图1所示，本发明实施例提供的产品可靠性评估系统包括：
数据获取模块1，用于获取加速退化数据；
模型构建模块2，用于构建典型的加速退化模型；
构建及推导模块3，用于基于Wiener模型构建加速应力下的剩余寿命预测模型以及推导出函数表达式；
参数估计模块4，用于结合性能加速退化数据的特征给出模型未知参数的估计方法。
预测模块5，用于得到加速应力下基于Wiener模型的剩余寿命预测的估计值。
如图2所示，本发明实施例提供的产品可靠性评估方法包括以下步骤：
S201：获取加速退化数据；
S202：构建典型的加速退化模型；
S203：基于Wiener模型构建加速应力下的剩余寿命预测模型以及推导出函数表达式；
S204：结合性能加速退化数据的特征给出模型未知参数的估计方法；
S205：得到加速应力下基于Wiener模型的剩余寿命预测的估计值。
下面结合附图对本发明的应用原理作进一步的描述。
如图3所示，本发明实施例提供的产品可靠性评估方法具体包括以下步骤：
步骤一，收集性能加速退化数据。试验时产品试样个数为m个，其施加的温度应力是Sk(k＝1,2,…,m)。在这m个温度应力水平下，分别对产品试样进行n次性能退化数据监测、记录。设表示第i次对产品性能的测量时刻，这时产品所增加的温度应力为Sk，测得的产品退化量为
步骤二，构建经典的加速退化模型，加速退化试验中，加速关系常指的是退化率与压力之间的相互作用，一般用来表示这种关系的加速模型形式为：
其中a、b为系数，表示一种以应力S为自变量函数，现有文献中用来描述温度应力的模型，一般采用Arrhenius模型来表示，即
则温度应力下加速模型为：
μk＝aexp(-b/Sk) (2)
式中Sk为第k个应力；μk为漂移系数，
进一步，基于Wiener过程的应力加速退化数据建模如下：
X(t)＝μkt+σB(t) (3)
其中μk～N(μa exp(-b/Sk),
步骤三，基于Wiener模型构建加速应力下的剩余寿命预测模型以及推导出函数表达式：为预测应力下产品的剩余寿命，需要建立此应力下产品剩余寿命的预测模型。假定特定产品在时刻tl和t+tl的性能退化量为X(tl)和X(t+tl)，有：
X(t+tl)＝X(tl)+μkt+σB(t) (4)
其中μk～N(μa exp(-b/Sk),
已知当前的退化量X(tk)和失效阈值ξ，产品在时刻tk的剩余寿命t可表示为：
则产品的剩余寿命t的PDF表示为：
根据剩余寿命t的定义，得到产品可靠度函数为：
步骤四，结合性能加速退化数据的特征给出模型未知参数的估计方法。采用两步极大似然估计方法估计对数极大似然估计函数如下：
式中：是未知参数的集合。lnL(Θ)表示利用fminsearch函数所求解得到的最大值。
分别令式(8)于ak、的偏导数为零，即：
则：
根据可得到的数据组及本发明所建模型，可以计算式(9)、式(10)求出式(9)、式(10)的求解，与b的赋值有关。因此，模型中的未知的参数可利用两步极大似然估计法来求解，先估计再估计
步骤五：得到加速应力下基于Wiener模型的剩余寿命预测的估计值。采用两步极大似然估计方法得到
下面结合仿真对本发明的应用效果作详细的描述。
1.仿真条件
本发明通过对温度应力下的电缆断裂伸长率的试验数据的实验仿真来说明该算法的有效性。仿真实验是在一台4G内存，ENVISION双核3.20GHz，32位win7操作系统下，，使用matlab2016b进行的。
2.仿真内容
采用温度应力下的电缆断裂伸长率的试验数据进行仿真。在仿真过程中，取4个电缆样本，分别在135℃、150℃、165℃、180℃的温度应力下，进行加速退化试验。每个应力下采集了10组数据，且每组数据的采集时间是不确定的。其中，580％是断裂伸长率的原始值，电缆正常工作时的温度是90℃。从而根据这些试验数据计算得出的估计值以及剩余寿命的预测值。
表1和图3是在四种温度应力下电缆断裂伸长率以及电缆断裂伸长率保留率与时间的关系曲线图。
表2是在四种温度应力下分别求出的模型未知参数的结果。从图3中可以知道，和都随着135℃、150℃、165℃、180℃温度的增加依次增大，相对和来说变化不大，说明Wiener退化模型的和都受温度应力的影响，而所受影响较小。
表3是依照建立的Wiener退化模型和已有的两种模型，分别求出的模型未知参数μ、σ2以及产品寿命置信度为95％的区间估计CI值，三种模型的AIC值，三种模型的TMSN值。由图可以看出，对于三种模型的各种参数值的对比，本发明中的模型的CI值，AIC值，TMSN值都是最小的。说明本发明中的模型比比其它两种模型的拟合性要好，也说明本发明的模型精确度高，更加的符合实际情况。
表4是分别取1000h、1500h和2000h的时间点时对电缆剩余寿命预测，得到的剩余寿命均值。本发明的模型在每个时刻预测的剩余寿命结果都是最接近真实值的，反映了本发明的模型比较好，更贴近样本的实际退化情况。另外两种模型与真实值相差比较大，这是因为没有考虑应力以及个体差异性对退化过程的影响。剩余寿命预测的结果也就与实际剩余寿命有偏差。
表1烯烃绝缘电缆在不同老化温度下的测试结果表
表2不同温度应力下的模型参数估计结果表
表3三种方法的估计结果表
表4三种方法的剩余寿命预测对比表
在上述实施例中，可以全部或部分地通过软件、硬件、固件或者其任意组合来实现。当使用全部或部分地以计算机程序产品的形式实现，所述计算机程序产品包括一个或多个计算机指令。在计算机上加载或执行所述计算机程序指令时，全部或部分地产生按照本发明实施例所述的流程或功能。所述计算机可以是通用计算机、专用计算机、计算机网络、或者其他可编程装置。所述计算机指令可以存储在计算机可读存储介质中，或者从一个计算机可读存储介质向另一个计算机可读存储介质传输，例如，所述计算机指令可以从一个网站站点、计算机、服务器或数据中心通过有线(例如同轴电缆、光纤、数字用户线(DSL)或无线(例如红外、无线、微波等)方式向另一个网站站点、计算机、服务器或数据中心进行传输)。所述计算机可读取存储介质可以是计算机能够存取的任何可用介质或者是包含一个或多个可用介质集成的服务器、数据中心等数据存储设备。所述可用介质可以是磁性介质，(例如，软盘、硬盘、磁带)、光介质(例如，DVD)、或者半导体介质(例如固态硬盘SolidState Disk(SSD))等。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (6)
1.一种产品可靠性评估方法，其特征在于，所述产品可靠性评估方法建立在一种基于随机效应的Wiener加速退化模型，其中漂移参数为服从正态分布的随机变量；根据加速退化模型构建剩余寿命预测模型，且推导出剩余寿命分布及剩余寿命的可靠度函数；根据应力状态下的性能退化数据，通过两步极大似然估计法，获得了模型未知参数估计值。
2.如权利要求1所述的产品可靠性评估方法，其特征在于，所述产品可靠性评估方法包括以下步骤：
步骤一，收集性能加速退化数据，产品试样个数为m个，施加的应力如温度是Sk(k＝1,2,…,m)；在这m个温度应力水平下，分别对产品试样进行n次性能退化数据监测、记录；表示第i次对产品性能的测量时刻，产品所增加的温度应力为Sk，测得的产品退化量为i＝1,2,…,n；k＝1,2,…,m；
步骤二，构建经典的加速退化模型，加速退化试验中，加速关系常指的是退化率与压力之间的相互作用，加速模型为：
其中a、b为系数，表示以应力S为自变量函数，采用Arrhenius模型表示，即温度应力下加速模型为：
μk＝aexp(-b/Sk)；
式中Sk为第k个应力；μk为漂移系数，
基于Wiener过程的应力加速退化数据建模：
X(t)＝μkt+σB(t)；
其中
步骤三，特定产品在时刻tl和t+tl的性能退化量为X(tl)和X(t+tl)，有：
X(t+tl)＝X(tl)+μkt+σB(t)；
其中
当前的退化量X(tk)和失效阈值ξ，产品在时刻tk的剩余寿命t表示为：
则产品的剩余寿命t的PDF表示为：
根据剩余寿命t的定义，得到产品可靠度函数为：
步骤四，采用两步极大似然估计方法估计对数极大似然估计函数如下：
式中：是未知参数的集合；lnL(Θ)表示利用fminsearch函数所求解得到的最大值；
ak、的偏导数为零：
则：
根据可得到的数据组及本发明所建模型，模型中的未知的参数可利用两步极大似然估计法来求解，先估计再估计步骤五：得到加速应力下基于Wiener模型的剩余寿命预测的估计值；采用两步极大似然估计方法得到
3.一种实施权利要求1所述产品可靠性评估方法的产品可靠性评估系统，其特征在于，所述产品可靠性评估系统包括：
数据获取模块，用于获取加速退化数据；
模型构建模块，用于构建典型的加速退化模型；
构建及推导模块，用于基于Wiener模型构建加速应力下的剩余寿命预测模型以及推导出函数表达式；
参数估计模块，用于结合性能加速退化数据的特征给出模型未知参数的估计方法；
预测模块，用于得到加速应力下基于Wiener模型的剩余寿命预测的估计值。
4.一种实现权利要求1～2任意一项所述产品可靠性评估方法的计算机程序。
5.一种实现权利要求1～2任意一项所述产品可靠性评估方法的信息数据处理终端。
6.一种计算机可读存储介质，包括指令，当其在计算机上运行时，使得计算机执行如权利要求1-2任意一项所述的产品可靠性评估方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811455684.9A CN109657937B (zh) | 2018-11-30 | 2018-11-30 | 一种基于退化数据的产品可靠性评估与寿命预测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811455684.9A CN109657937B (zh) | 2018-11-30 | 2018-11-30 | 一种基于退化数据的产品可靠性评估与寿命预测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109657937A true CN109657937A (zh) | 2019-04-19 |
| CN109657937B CN109657937B (zh) | 2022-12-30 |
Family
ID=66112139
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811455684.9A Active CN109657937B (zh) | 2018-11-30 | 2018-11-30 | 一种基于退化数据的产品可靠性评估与寿命预测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109657937B (zh) |
Cited By (21)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110245386A (zh) * | 2019-05-17 | 2019-09-17 | 桂林电子科技大学 | 一种功率器件可靠性的评估方法、装置以及存储介质 |
| CN110595742A (zh) * | 2019-09-18 | 2019-12-20 | 广东产品质量监督检验研究院（国家质量技术监督局广州电气安全检验所、广东省试验认证研究院、华安实验室） | 一种机械载荷对光伏组件性能长期潜在影响的检测方法 |
| CN111046564A (zh) * | 2019-12-16 | 2020-04-21 | 重庆大学 | 两阶段退化产品的剩余寿命预测方法 |
| CN111177863A (zh) * | 2019-12-16 | 2020-05-19 | 北京动力机械研究所 | 一种针对机械磨损过程随机特征的退化寿命预估方法 |
| CN111523251A (zh) * | 2020-06-09 | 2020-08-11 | 江苏科技大学 | 一种随机环境应力下的产品寿命快速评估方法 |
| CN112084688A (zh) * | 2020-08-18 | 2020-12-15 | 中国电子产品可靠性与环境试验研究所（（工业和信息化部电子第五研究所）（中国赛宝实验室）） | 阴极寿命预测方法 |
| CN112395805A (zh) * | 2020-10-27 | 2021-02-23 | 中国电子产品可靠性与环境试验研究所（（工业和信息化部电子第五研究所）（中国赛宝实验室）） | Ebsm阴极寿命评估方法 |
| CN112685912A (zh) * | 2021-01-12 | 2021-04-20 | 河南科技大学 | 一种多元广义Wiener过程性能退化可靠性分析方法 |
| CN112800616A (zh) * | 2021-02-05 | 2021-05-14 | 中国人民解放军空军工程大学 | 基于比例加速退化建模的设备剩余寿命自适应预测方法 |
| CN113010971A (zh) * | 2019-12-20 | 2021-06-22 | 新疆金风科技股份有限公司 | 模块化齿形带、性能评估方法及变桨驱动装置 |
| CN113032998A (zh) * | 2021-03-26 | 2021-06-25 | 北京中关村水木医疗科技有限公司 | 医疗器械寿命评估方法和装置 |
| CN113032985A (zh) * | 2021-03-11 | 2021-06-25 | 北京必创科技股份有限公司 | 一种无线感知设备寿命智能评估方法和装置 |
| CN114417534A (zh) * | 2022-02-21 | 2022-04-29 | 北京科技大学 | 一种基于Wiener过程和P-EMD的机械结构剩余寿命预测方法 |
| CN114692542A (zh) * | 2022-04-25 | 2022-07-01 | 西安电子科技大学 | 基于动态步长的模拟电路可靠性仿真方法 |
| CN115841046A (zh) * | 2023-02-10 | 2023-03-24 | 中国电子产品可靠性与环境试验研究所（（工业和信息化部电子第五研究所）（中国赛宝实验室）） | 基于维纳过程的加速退化试验数据处理方法和装置 |
| CN116227240A (zh) * | 2023-05-08 | 2023-06-06 | 中国电子产品可靠性与环境试验研究所（（工业和信息化部电子第五研究所）（中国赛宝实验室）） | 基于综合应力加速试验的产品寿命评价方法、装置及设备 |
| CN116228045A (zh) * | 2023-05-09 | 2023-06-06 | 中国电子产品可靠性与环境试验研究所（（工业和信息化部电子第五研究所）（中国赛宝实验室）） | 基于性能退化的产品可靠性薄弱环节评估方法和装置 |
| CN116502391A (zh) * | 2022-12-12 | 2023-07-28 | 哈尔滨工业大学 | 基于误差函数模型的寿命预测方法 |
| CN116756505A (zh) * | 2023-06-07 | 2023-09-15 | 上海正泰电源系统有限公司 | 一种基于大数据的光伏设备智能管理系统及方法 |
| CN117094169A (zh) * | 2023-09-05 | 2023-11-21 | 西南科技大学 | 基于halt试验的afss吸波器可靠性评估方法 |
| US12079550B1 (en) | 2024-02-06 | 2024-09-03 | China Electronic Product Reliability and Environmental Testing Research Institute (The Fifth Electronic Research Institute of Ministry of Industry and Information Technology) (CEPREI) | Wiener process-based method and device for processing accelerated degradation testing data |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106484949A (zh) * | 2016-09-12 | 2017-03-08 | 西安理工大学 | 基于退化数据的动量轮可靠性分析与剩余寿命预测方法 |
| CN107862134A (zh) * | 2017-11-06 | 2018-03-30 | 河南科技大学 | 一种考虑自相关测量误差的Wiener过程可靠性分析方法 |
| CN107885928A (zh) * | 2017-11-06 | 2018-04-06 | 河南科技大学 | 考虑测量误差的步进应力加速性能退化可靠性分析方法 |
| WO2018076475A1 (zh) * | 2016-10-26 | 2018-05-03 | 广东产品质量监督检验研究院 | 一种基于深度学习法构建的光伏组件加速退化模型及光伏组件寿命预测方法 |
-
2018
- 2018-11-30 CN CN201811455684.9A patent/CN109657937B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106484949A (zh) * | 2016-09-12 | 2017-03-08 | 西安理工大学 | 基于退化数据的动量轮可靠性分析与剩余寿命预测方法 |
| WO2018076475A1 (zh) * | 2016-10-26 | 2018-05-03 | 广东产品质量监督检验研究院 | 一种基于深度学习法构建的光伏组件加速退化模型及光伏组件寿命预测方法 |
| CN107862134A (zh) * | 2017-11-06 | 2018-03-30 | 河南科技大学 | 一种考虑自相关测量误差的Wiener过程可靠性分析方法 |
| CN107885928A (zh) * | 2017-11-06 | 2018-04-06 | 河南科技大学 | 考虑测量误差的步进应力加速性能退化可靠性分析方法 |
Cited By (34)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110245386A (zh) * | 2019-05-17 | 2019-09-17 | 桂林电子科技大学 | 一种功率器件可靠性的评估方法、装置以及存储介质 |
| CN110245386B (zh) * | 2019-05-17 | 2022-05-31 | 桂林电子科技大学 | 一种功率器件可靠性的评估方法、装置以及存储介质 |
| CN110595742A (zh) * | 2019-09-18 | 2019-12-20 | 广东产品质量监督检验研究院（国家质量技术监督局广州电气安全检验所、广东省试验认证研究院、华安实验室） | 一种机械载荷对光伏组件性能长期潜在影响的检测方法 |
| CN111046564A (zh) * | 2019-12-16 | 2020-04-21 | 重庆大学 | 两阶段退化产品的剩余寿命预测方法 |
| CN111177863A (zh) * | 2019-12-16 | 2020-05-19 | 北京动力机械研究所 | 一种针对机械磨损过程随机特征的退化寿命预估方法 |
| CN111177863B (zh) * | 2019-12-16 | 2023-11-07 | 北京动力机械研究所 | 一种针对机械磨损过程随机特征的退化寿命预估方法 |
| CN111046564B (zh) * | 2019-12-16 | 2024-04-05 | 重庆大学 | 两阶段退化产品的剩余寿命预测方法 |
| CN113010971A (zh) * | 2019-12-20 | 2021-06-22 | 新疆金风科技股份有限公司 | 模块化齿形带、性能评估方法及变桨驱动装置 |
| CN111523251A (zh) * | 2020-06-09 | 2020-08-11 | 江苏科技大学 | 一种随机环境应力下的产品寿命快速评估方法 |
| CN111523251B (zh) * | 2020-06-09 | 2023-04-21 | 江苏科技大学 | 一种随机环境应力下的产品寿命快速评估方法 |
| CN112084688A (zh) * | 2020-08-18 | 2020-12-15 | 中国电子产品可靠性与环境试验研究所（（工业和信息化部电子第五研究所）（中国赛宝实验室）） | 阴极寿命预测方法 |
| CN112084688B (zh) * | 2020-08-18 | 2023-10-03 | 中国电子产品可靠性与环境试验研究所（（工业和信息化部电子第五研究所）（中国赛宝实验室）） | 阴极寿命预测方法 |
| CN112395805A (zh) * | 2020-10-27 | 2021-02-23 | 中国电子产品可靠性与环境试验研究所（（工业和信息化部电子第五研究所）（中国赛宝实验室）） | Ebsm阴极寿命评估方法 |
| CN112395805B (zh) * | 2020-10-27 | 2023-10-03 | 中国电子产品可靠性与环境试验研究所（（工业和信息化部电子第五研究所）（中国赛宝实验室）） | Ebsm阴极寿命评估方法 |
| CN112685912A (zh) * | 2021-01-12 | 2021-04-20 | 河南科技大学 | 一种多元广义Wiener过程性能退化可靠性分析方法 |
| CN112685912B (zh) * | 2021-01-12 | 2023-02-28 | 河南科技大学 | 一种多元广义Wiener过程性能退化可靠性分析方法 |
| CN112800616A (zh) * | 2021-02-05 | 2021-05-14 | 中国人民解放军空军工程大学 | 基于比例加速退化建模的设备剩余寿命自适应预测方法 |
| CN113032985A (zh) * | 2021-03-11 | 2021-06-25 | 北京必创科技股份有限公司 | 一种无线感知设备寿命智能评估方法和装置 |
| CN113032985B (zh) * | 2021-03-11 | 2024-04-26 | 北京必创科技股份有限公司 | 一种无线感知设备寿命智能评估方法和装置 |
| CN113032998A (zh) * | 2021-03-26 | 2021-06-25 | 北京中关村水木医疗科技有限公司 | 医疗器械寿命评估方法和装置 |
| CN114417534A (zh) * | 2022-02-21 | 2022-04-29 | 北京科技大学 | 一种基于Wiener过程和P-EMD的机械结构剩余寿命预测方法 |
| CN114692542B (zh) * | 2022-04-25 | 2024-05-14 | 西安电子科技大学 | 基于动态步长的模拟电路可靠性仿真方法 |
| CN114692542A (zh) * | 2022-04-25 | 2022-07-01 | 西安电子科技大学 | 基于动态步长的模拟电路可靠性仿真方法 |
| CN116502391A (zh) * | 2022-12-12 | 2023-07-28 | 哈尔滨工业大学 | 基于误差函数模型的寿命预测方法 |
| CN115841046A (zh) * | 2023-02-10 | 2023-03-24 | 中国电子产品可靠性与环境试验研究所（（工业和信息化部电子第五研究所）（中国赛宝实验室）） | 基于维纳过程的加速退化试验数据处理方法和装置 |
| CN116227240A (zh) * | 2023-05-08 | 2023-06-06 | 中国电子产品可靠性与环境试验研究所（（工业和信息化部电子第五研究所）（中国赛宝实验室）） | 基于综合应力加速试验的产品寿命评价方法、装置及设备 |
| CN116227240B (zh) * | 2023-05-08 | 2023-08-04 | 中国电子产品可靠性与环境试验研究所（（工业和信息化部电子第五研究所）（中国赛宝实验室）） | 基于综合应力加速试验的产品寿命评价方法、装置及设备 |
| CN116228045B (zh) * | 2023-05-09 | 2023-09-01 | 中国电子产品可靠性与环境试验研究所（（工业和信息化部电子第五研究所）（中国赛宝实验室）） | 基于性能退化的产品可靠性薄弱环节评估方法和装置 |
| CN116228045A (zh) * | 2023-05-09 | 2023-06-06 | 中国电子产品可靠性与环境试验研究所（（工业和信息化部电子第五研究所）（中国赛宝实验室）） | 基于性能退化的产品可靠性薄弱环节评估方法和装置 |
| CN116756505A (zh) * | 2023-06-07 | 2023-09-15 | 上海正泰电源系统有限公司 | 一种基于大数据的光伏设备智能管理系统及方法 |
| CN116756505B (zh) * | 2023-06-07 | 2024-05-28 | 上海正泰电源系统有限公司 | 一种基于大数据的光伏设备智能管理系统及方法 |
| CN117094169A (zh) * | 2023-09-05 | 2023-11-21 | 西南科技大学 | 基于halt试验的afss吸波器可靠性评估方法 |
| CN117094169B (zh) * | 2023-09-05 | 2024-05-24 | 西南科技大学 | 基于halt试验的afss吸波器可靠性评估方法 |
| US12079550B1 (en) | 2024-02-06 | 2024-09-03 | China Electronic Product Reliability and Environmental Testing Research Institute (The Fifth Electronic Research Institute of Ministry of Industry and Information Technology) (CEPREI) | Wiener process-based method and device for processing accelerated degradation testing data |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109657937B (zh) | 2022-12-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109657937A (zh) | 一种基于退化数据的产品可靠性评估与寿命预测方法 | |
| Liao et al. | Optimal design for step-stress accelerated degradation tests | |
| CN107862134B (zh) | 一种考虑自相关测量误差的Wiener过程可靠性分析方法 | |
| CN111859658B (zh) | 一种产品贮存寿命与可靠性评估方法 | |
| CN109829137B (zh) | 一种周期应力下非线性退化设备的寿命预测方法及系统 | |
| CN107436983B (zh) | 一种基于多元样本差异的o型橡胶密封圈寿命预测方法 | |
| CN107885928B (zh) | 考虑测量误差的步进应力加速性能退化可靠性分析方法 | |
| CN110081923B (zh) | 野外基线环境参数自动采集系统故障检测方法与装置 | |
| Yi et al. | Acceleration model considering multi‐stress coupling effect and reliability modeling method based on nonlinear Wiener process | |
| CN110059894A (zh) | 设备状态评估方法、装置、系统及存储介质 | |
| EP2838022A1 (en) | User sensory quality estimation device, terminal bottleneck determination device, similar operation extraction device, and methods and programs therefor | |
| CN109345041B (zh) | 一种利用威布尔分布与arma结合的设备故障率预测方法 | |
| CN118657057B (zh) | 基于Tweedie指数性能指标相关加速退化试验方法 | |
| CN115630490A (zh) | 基于多性能参数退化的电机可靠性建模与寿命预测方法 | |
| Sun et al. | Accelerated degradation process analysis based on the nonlinear Wiener process with covariates and random effects | |
| Biswas et al. | A model-agnostic method for PMU data recovery using optimal singular value thresholding | |
| CN119310401B (zh) | 一种配电网中防雷智能监测方法及系统 | |
| Hua et al. | Efficient probabilistic contingency analysis through a stability measure considering wind perturbation | |
| CN116070302A (zh) | 一种电缆绝缘状态的预测方法及装置 | |
| Švenda et al. | Flexible hybrid state estimation for power systems with communication irregularities | |
| CN117561518A (zh) | 用于检测和解释异常的方法和设备 | |
| CN116756505B (zh) | 一种基于大数据的光伏设备智能管理系统及方法 | |
| CN115225455B (zh) | 异常设备检测方法及装置、电子设备、存储介质 | |
| CN114792053B (zh) | 一种基于初值-速率相关退化模型的可靠性评估方法 | |
| CN114021302B (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 |