CN103488826A - Experience acceleration model based degradation amount distribution parameter modeling and extrapolating method - Google Patents

Experience acceleration model based degradation amount distribution parameter modeling and extrapolating method Download PDF

Info

Publication number
CN103488826A
CN103488826A CN201310415072.8A CN201310415072A CN103488826A CN 103488826 A CN103488826 A CN 103488826A CN 201310415072 A CN201310415072 A CN 201310415072A CN 103488826 A CN103488826 A CN 103488826A
Authority
CN
China
Prior art keywords
formula
degradation
stress
distribution parameter
amount
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
Application number
CN201310415072.8A
Other languages
Chinese (zh)
Other versions
CN103488826B (en
Inventor
党香俊
姜同敏
封雷
孙富强
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Original Assignee
Beihang University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beihang University filed Critical Beihang University
Priority to CN201310415072.8A priority Critical patent/CN103488826B/en
Publication of CN103488826A publication Critical patent/CN103488826A/en
Application granted granted Critical
Publication of CN103488826B publication Critical patent/CN103488826B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses an experience acceleration model based degradation amount distribution parameter modeling and extrapolating method. The method includes the following steps: firstly, selecting an experience acceleration model corresponding to an accelerated life test; secondly, establishing a degradation amount distribution parameter acceleration model of a relation of degradation amount distribution parameters and stress; thirdly, extrapolating a degradation amount distribution parameter value under target stress by the aid of the degradation amount distribution parameter acceleration model. The experience acceleration model which already undergoes a large amount of research and is proved by application in the field of accelerated life tests is adopted for modeling, relevant test information in the product research process and test information of relevant products can be utilized effectively, and credibility of test results is improved. On the basis of model extrapolation, the problem that accuracy is difficult to guarantee when a regression method is purely adopted to describe the relation of the distribution parameter and the stress under the circumstance of small samples is avoided.

Description

Amount of degradation distribution parameter modeling Extrapolation method based on the experience acceleration model
Technical field
The present invention relates to a kind of amount of degradation distribution parameter modeling Extrapolation method based on the experience acceleration model, belong to the accelerated degradation test technical field.
Background technology
When suffered stress, for normal operation during stress, the performance degradation of high reliability long life product is process very slowly.Under this condition, expectation obtains enough degraded datas, certainly will will experience the longer time and pay corresponding manpower and materials cost.The data volume needed in order to obtain analyzing and processing within a short period of time, in laboratory environment, adopt the method that improves the harsh level of the suffered stress of product, the degenerative process of expedite product, be one of focus direction of current solution high reliability long life prediction and evaluation problem.Develop into acceleration stress degradation experiment from the degradation experiment of normal applied stress, test period and relevant support expense cost have obtained decrease.And, when this changes into and shortens product development cycle and facilitate, corresponding data are processed and are also become more difficult.The degraded data will speed up under stress is transformed under regular service condition, is one of Major Difficulties accelerated the degraded data processing.List of references [1] (Zhao Jianyin, Sun Quan, Peng Baohua, Zhou Jinglun. the fail-safe analysis based on the accelerated degradation test data, electron mass, 2005 (7), acceleration degeneration equation concept has been proposed 30-33), suppose that there is derivative in the amount of degradation function, the Arrhenius of take in the reaction rate model, as basis, has provided acceleration degeneration equation and logarithmic form thereof.Yet, in existing acceleration model, except the reaction rate model, also have application very extensively against the experience acceleration models such as power rate model (list of references [2]: yellow graceful, Jiang Tongmin. statistics acceleration model summary in accelerated life test, equipment Environmental Engineering, 2010,7 (4), 57-62).Because the experience acceleration models such as contrary power rate are not that to take reaction rate be basis, the modeling process that accelerates the degeneration equation in list of references [1] also is not suitable for the experience acceleration model.
For this problem, the very ripe experience acceleration model in the accelerated life test field be take in this paper is basis, sets up the amount of degradation distribution parameter acceleration model of reflection amount of degradation distribution parameter and stress relation.In the model inference process, the widely used contrary power rate model of take is example, sets up corresponding numerical relationship model, and the amount of degradation distribution parameter under the extrapolation normal stress, be convenient to the practical application of accelerated degradation test.
Summary of the invention
The objective of the invention is, in order to solve the relationship modeling technical matters of amount of degradation distribution parameter and stress in accelerated degradation test, to propose a kind of amount of degradation distribution parameter modeling Extrapolation method based on the experience acceleration model.
Described method comprises the steps:
(1) select the experience acceleration model of corresponding accelerated life test;
(2) set up the amount of degradation distribution parameter acceleration model of amount of degradation distribution parameter and stress relation;
(3) utilize the amount of degradation distribution parameter value under gained amount of degradation distribution parameter acceleration model extrapolation target stress.
The invention has the advantages that:
(1) adopt and carried out modeling in the accelerated life test field by the experience acceleration model of large quantity research and application attestation, can effectively utilize correlation test information in the product research process and the Test Information of Related product, improved the credibility of test findings;
(2) avoided under Small Sample Size, while adopting merely homing method to describe distribution parameter and stress relation, accuracy is difficult to the problem guaranteed.
The accompanying drawing explanation
Fig. 1 is method flow diagram of the present invention.
Embodiment
Below in conjunction with drawings and Examples, the present invention is described in further detail.
The present invention is a kind of amount of degradation distribution parameter modeling Extrapolation method based on the experience acceleration model, flow process as shown in Figure 1, and specific implementation method is as follows:
Step 1, select the experience acceleration model of corresponding accelerated life test;
For the experience acceleration model, be based on the slip-stick artist and the properties of product long-term observation summed up and proposed, and obtained in actual applications the confirmation of great number tested data.For example contrary power rate model description the relation between voltage or pressure and life of product, the Coffin-Manson model description relation between temperature cycles stress and life of product.In the hypothesis test process, the acceleration stress that product is applied is voltage stress, and corresponding acceleration model is selected contrary power rate model.
Contrary power rate model representation for accelerated life test is:
ξ=Av -C (1)
In formula: ξ is certain life characteristics; A is a normal number; C is a normal number relevant with activation energy; V is stress, and normal power taking is pressed.
Step 2, set up the amount of degradation distribution parameter acceleration model of amount of degradation distribution parameter and stress relation;
For degenerating the life-span, be commonly defined as product performance index from the degeneration initial time to the time span that reaches defined threshold and experience, this defined threshold is the degradation failure threshold value of product, usually uses D fmean.
ξ=f(D f) -1 (2)
In formula, the universal expression formula that f (.) is the performance degradation function, adopt the form of implicit function herein, so that the expression in formulation process.
In fact, the determined value on the degradation failure threshold value is not strict, relevant with actual demand.For example, for the degeneration life-span of critical product or parts, in order to reduce the risk of product failure, improve margin of safety, it is more conservative that the degradation failure threshold value can be determined; Otherwise, can be decided to be more optimistic index.Degradation failure is also referred to as " soft failure ", with demblee form, lost efficacy, " hard failure " compare, certain difference is arranged on the determinacy in " life-span ".During the degradation failure changes of threshold, " life-span " of degradation failure also changes thereupon.Therefore, under theoretical case, the degradation failure threshold value can get arbitrarily on the occasion of, correspondingly have:
t=f(d f) -1,d f≥0 (3)
In formula, variable d fmean degradation failure threshold value D fcan in the degeneration value is positive scope, change continuously, the variable t life-span ξ that means to degenerate accordingly is also continually varying time value.
Contrary power rate model description the relation of life of product and suffered stress, with the variation of life of product, do not change.So same life-span constantly, during t=ξ, formula (1) and formula (2) can be combined into:
f(d f) -1=Av -C (4)
No longer include time dependent amount in formula.And d fin fact just be selected as the degeneration value of " degradation failure threshold value ", with general amount of degradation, there is identical field of definition, might as well replace d with general amount of degradation M f, so that the derivation of equation is more directly perceived, obtain:
f(M) -1=Av -C,M=d f (5)
Formula (5) is taken the logarithm:
ln(f(M) -1)=lnA-Clnv (6)
Can find out the functional transformation ln of general amount of degradation M (f (M) from formula (6) -1) and the natural logarithm ln ν of stress between be linear relationship, coefficient C wherein and lnA(or other constant terms) can obtain by the linear regression to data.Because the distribution parameters such as average μ, scale parameter η are the index that embodies the data general characteristic, its meaning has embodied the amount of degradation of " totally ", the category that belongs to general amount of degradation, therefore itself and the same coincidence formula of relation (5) of stress and the model of formula (6), be the relation that formula (5) and formula (6) have been described amount of degradation distribution parameter and stress, this conclusion is applicable to formula (9) and formula (10) hereinafter too.
It is example that the present invention be take against power rate model, has provided the amount of degradation distribution parameter acceleration model derivation based on the experience acceleration model.For other experience acceleration models, suppose that its life-span acceleration model has following common version:
Figure BDA0000381098750000042
Wherein
Figure BDA0000381098750000043
for the function of stress S, be the generic representation form on experience acceleration model equation right side.Identical with the front derivation, when t=ξ, in conjunction with formula (7) and formula (3), formula (4) and the corresponding universal expression form of formula (5) are arranged:
Figure BDA0000381098750000044
Figure BDA0000381098750000045
Model parameter to be estimated in formula (9) can be carried out proper transformation with reference to formula (6) in conjunction with concrete equation form, by the simple regression method, just can obtain, and does not repeat them here simultaneously.
Step 3, utilize the amount of degradation distribution parameter value under gained amount of degradation distribution parameter acceleration model extrapolation target stress;
Take against power rate model is example, according to Degradation path, obtains performance degenrate function f (.), and under application of formula (6) employing acceleration stress, the data of amount of degradation distribution parameter under different stress, return the parameters C obtained in i moment acceleration model formula iand lnA i(or other constant terms).Under normal stress to any i moment extrapolation, can be by normal stress value ν 0substitution extrapolation formula now:
ln ( f ( M i 0 ) - 1 ) = ln A i - C i ln v 0 - - - ( 10 )
Solve and obtain.Herein, the extrapolation of other experience acceleration models solves identical therewith.
Embodiment 1: in the accelerated degradation test of certain product, the acceleration stress applied is voltage stress, and the Mean Parameters that under 90h measured each stress of the moment, amount of degradation distributes is: (7.5V, 0.14750), (7.0V, 0.0652), (6.7V, 0.0316).In the situation that do not affect method validation, the magnitude of Mean Parameters has been done relevant treatment, and unit omits.The degeneration initial value of this product is 0, and Degradation path function f (.) has the form of f (t)=at.Below adopt the average under the extracting method extrapolation 6.5V of institute of the present invention normal stress.
Step 1, select the experience acceleration model of corresponding accelerated life test;
Because the stress that product is applied is voltage stress, the contrary power rate model in the experience acceleration model is applicable to describing the relation of this life of product and stress, and selecting contrary power rate model is accelerated life model, and its expression formula is as shown in formula (1).
Step 2, set up the amount of degradation distribution parameter acceleration model of amount of degradation distribution parameter and stress relation;
The modeling process of amount of degradation distribution parameter acceleration model, referring to the content in step 2 above, no longer repeats herein.
Step 3, utilize the amount of degradation distribution parameter value under gained amount of degradation distribution parameter acceleration model extrapolation target stress;
For the 90h data in the moment, by Degradation path function f (.) substitution formula (5), have
a - 1 μ 90 = A 90 v - C 90 - - - ( 11 )
After taking the logarithm, conversion obtains
lnμ 90=ln(aA 90)-C 90lnv (12)
Through linear regression, try to achieve
C 90=-13.5003,ln(aA 90)=-29.0834
By ν 0=6.5 substitution formula (12) obtain the average under the 6.5V normal stress .

Claims (1)

1. the amount of degradation distribution parameter modeling Extrapolation method based on the experience acceleration model, is characterized in that, comprises following step:
Step 1, select the experience acceleration model of corresponding accelerated life test;
If the acceleration stress that product is applied is voltage stress, select contrary power rate model;
Now the contrary power rate model representation for accelerated life test is:
ξ=Av -C (1)
In formula: ξ is certain life characteristics; A is a normal number; C is a normal number relevant with activation energy; V is stress; Step 2, set up the amount of degradation distribution parameter acceleration model of amount of degradation distribution parameter and stress relation;
Use D fmean the degradation failure threshold value of product, the life characteristics of product is:
ξ=f(D f) -1 (2)
The universal expression formula that in formula, f (.) is the performance degradation function, adopt the form of implicit function herein;
Under theoretical case, the degradation failure threshold value get arbitrarily on the occasion of, correspondingly have:
t=f(d f) -1,d f≥0 (3)
In formula, variable d fmean degradation failure threshold value D fin the degeneration value is positive scope, change continuously, the variable t life-span ξ that means to degenerate accordingly is also continually varying time value;
Contrary power rate model description the relation of life of product with suffered stress, with the variation of life of product, do not change, so, in the same moment in life-span, during t=ξ, formula (1) and formula (2) are combined into:
f(d f) -1=Av -C (4)
Replace d with general amount of degradation M f, have:
f(M) -1=Av -C,M=d f (5)
Formula (5) is taken the logarithm:
ln(f(M) -1)=lnA-Clnv (6)
Formula (5) and formula (6) have been described the amount of degradation distribution parameter and relation stress;
Suppose that the life-span acceleration model has following common version:
Figure FDA0000381098740000011
Wherein
Figure FDA0000381098740000012
for the function of stress S, when t=ξ, in conjunction with formula (7) and formula (3), formula (4) and the corresponding universal expression form of formula (5) are arranged:
Figure FDA0000381098740000022
Figure FDA0000381098740000023
Step 3, utilize the amount of degradation distribution parameter value under gained amount of degradation distribution parameter acceleration model extrapolation target stress;
Obtain degenrate function f (.) according to Degradation path; the data of amount of degradation distribution parameter under different stress under application of formula (6) employing acceleration stress; recurrence obtains the parameter in i moment acceleration model formula, to amount of degradation under the normal stress in any i moment
Figure FDA0000381098740000024
extrapolation, by normal stress value v 0substitution extrapolation formula now:
ln ( f ( M i 0 ) - 1 ) = ln A i - C i ln v 0 - - - ( 10 )
Solve and obtain.
CN201310415072.8A 2013-09-12 2013-09-12 Amount of degradation distributed constant modeling Extrapolation method based on experience acceleration model Expired - Fee Related CN103488826B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310415072.8A CN103488826B (en) 2013-09-12 2013-09-12 Amount of degradation distributed constant modeling Extrapolation method based on experience acceleration model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310415072.8A CN103488826B (en) 2013-09-12 2013-09-12 Amount of degradation distributed constant modeling Extrapolation method based on experience acceleration model

Publications (2)

Publication Number Publication Date
CN103488826A true CN103488826A (en) 2014-01-01
CN103488826B CN103488826B (en) 2016-05-04

Family

ID=49829046

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310415072.8A Expired - Fee Related CN103488826B (en) 2013-09-12 2013-09-12 Amount of degradation distributed constant modeling Extrapolation method based on experience acceleration model

Country Status (1)

Country Link
CN (1) CN103488826B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104181457A (en) * 2014-08-15 2014-12-03 中国电子科技集团公司第二十四研究所 Method for selecting optimal semiconductor device temperature and humidity combined stress acceleration model
CN107515965A (en) * 2017-07-27 2017-12-26 北京航空航天大学 A kind of acceleration degeneration modelling evaluation method based on uncertain course
CN107704663A (en) * 2017-09-14 2018-02-16 中国电子科技集团公司第二十四研究所 A kind of semiconductor device temperature pulsating stress acceleration model method for optimizing

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101793927A (en) * 2010-01-12 2010-08-04 北京航空航天大学 Optimization design method of step-stress accelerated degradation test
US20120235701A1 (en) * 2011-03-14 2012-09-20 Universal Display Corporation Method for Accelerated Lifetesting of Large Area OLED Lighting Panels

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101793927A (en) * 2010-01-12 2010-08-04 北京航空航天大学 Optimization design method of step-stress accelerated degradation test
US20120235701A1 (en) * 2011-03-14 2012-09-20 Universal Display Corporation Method for Accelerated Lifetesting of Large Area OLED Lighting Panels

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张恒: "基于性能退化数据的可靠性建模与评估理论研究", 《万方学术论文库》, 31 December 2010 (2010-12-31), pages 1 - 75 *
茆诗松: "加速寿命试验的加速模型", 《质量与可靠性》, no. 02, 30 April 2003 (2003-04-30), pages 15 - 17 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104181457A (en) * 2014-08-15 2014-12-03 中国电子科技集团公司第二十四研究所 Method for selecting optimal semiconductor device temperature and humidity combined stress acceleration model
CN104181457B (en) * 2014-08-15 2017-01-18 中国电子科技集团公司第二十四研究所 Method for selecting optimal semiconductor device temperature and humidity combined stress acceleration model
CN107515965A (en) * 2017-07-27 2017-12-26 北京航空航天大学 A kind of acceleration degeneration modelling evaluation method based on uncertain course
CN107515965B (en) * 2017-07-27 2019-07-23 北京航空航天大学 A kind of acceleration degeneration modelling evaluation method based on uncertain course
CN107704663A (en) * 2017-09-14 2018-02-16 中国电子科技集团公司第二十四研究所 A kind of semiconductor device temperature pulsating stress acceleration model method for optimizing

Also Published As

Publication number Publication date
CN103488826B (en) 2016-05-04

Similar Documents

Publication Publication Date Title
Yang et al. A novel Gaussian process regression model for state-of-health estimation of lithium-ion battery using charging curve
CN101620045B (en) Method for evaluating reliability of stepping stress quickened degradation experiment based on time sequence
CN102542155B (en) Particle filter residual life forecasting method based on accelerated degradation data
CN102622473B (en) Optimization design method for step stress accelerated degradation test based on Bayesian theory
CN107436983B (en) O-shaped rubber sealing ring service life prediction method based on multivariate sample difference
CN111859658B (en) Product storage life and reliability assessment method
CN103310051B (en) Board information terminal Failure Rate Forecasting Method in a kind of life cycle management
CN104156500A (en) Method for predicting material fatigue life
Sun et al. Remaining useful life prediction of aviation circular electrical connectors using vibration-induced physical model and particle filtering method
Mohanty et al. Bayesian statistic based multivariate Gaussian process approach for offline/online fatigue crack growth prediction
CN110686915B (en) Method, system, medium and equipment for determining multi-stress acceleration test profile
CN108804806B (en) Simplified MLE (Multi-level engineering) method for parameters in Weibull distributed comprehensive stress constant addition test
CN108959745B (en) Product performance change and service life high confidence small sample statistics inference and confirmation method
CN104316872A (en) Pressure relay storage life prediction method based on step-down degradation acceleration test
Kang et al. Development of Hartigan’s dip statistic with bimodality coefficient to assess multimodality of distributions
Zhang et al. A novel nonlinear adaptive Mooney-viscosity model based on DRPLS-GP algorithm for rubber mixing process
CN105737881A (en) Acceleration reliability growth test method for electronic equipment
CN117471227B (en) Automobile wire harness parameter performance test method and test system
CN110728038A (en) Dam monitoring method based on regression analysis
CN103488826A (en) Experience acceleration model based degradation amount distribution parameter modeling and extrapolating method
He et al. Valve stiction quantification method based on a semiphysical valve stiction model
CN106056305B (en) Power generation system reliability rapid evaluation method based on state clustering
CN107301285A (en) Non-electronic product Sequential Compliance Method based on predicting residual useful life
Zhao et al. Rolling bearing remaining useful life prediction based on wiener process
CN103630588A (en) Rapid response method of electrochemical gas sensor

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160504

Termination date: 20200912

CF01 Termination of patent right due to non-payment of annual fee