CN111046564A - Method for predicting residual life of two-stage degraded product - Google Patents

Method for predicting residual life of two-stage degraded product Download PDF

Info

Publication number
CN111046564A
CN111046564A CN201911296295.0A CN201911296295A CN111046564A CN 111046564 A CN111046564 A CN 111046564A CN 201911296295 A CN201911296295 A CN 201911296295A CN 111046564 A CN111046564 A CN 111046564A
Authority
CN
China
Prior art keywords
stage
degradation
model
life
probability density
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
CN201911296295.0A
Other languages
Chinese (zh)
Other versions
CN111046564B (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.)
Chongqing University
Original Assignee
Chongqing 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 Chongqing University filed Critical Chongqing University
Priority to CN201911296295.0A priority Critical patent/CN111046564B/en
Publication of CN111046564A publication Critical patent/CN111046564A/en
Application granted granted Critical
Publication of CN111046564B publication Critical patent/CN111046564B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30Computing systems specially adapted for manufacturing

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Marketing (AREA)
  • Game Theory and Decision Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Development Economics (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

The invention discloses a residual life prediction method of a two-stage degradation product, belonging to the field of residual life prediction in prediction and health management, and mainly comprising the following steps: modeling a degradation process, estimating model parameters and predicting residual life; the modeling of the degradation process is to establish a two-stage degradation model by utilizing a nonlinear Wiener process; the model parameter estimation comprises: collecting historical degradation data; estimating model parameters based on a unit maximum likelihood estimation method; statistical analysis of random parameter distribution; the remaining life prediction includes: acquiring first arrival time distribution; deducing a state transition probability function; and predicting the residual life of the product based on the estimated parameters. The invention can effectively predict the residual service life of the two-stage degraded product, is favorable for ensuring the operational reliability of the product, reduces the maintenance cost and avoids safety accidents.

Description

Method for predicting residual life of two-stage degraded product
Technical Field
The invention belongs to the field of prediction and health management, and relates to a method for predicting the residual life of a two-stage degraded product.
Background
The residual life prediction is the core content of the prediction and health management technology, is the effective time interval left by the product at the current time when the product loses the specified function, is an important index for reflecting the reliability of the product, and has important significance for practically guaranteeing the operation safety, reliability and economy of the product.
The prediction of remaining life has been of great interest and intensive research over the last decade. Among them, the Wiener process has been developed because it can describe non-monotonic degradation trajectory and has good mathematical characteristics. It is worth noting that in most life prediction methods based on Wiener processes, the degradation rate is generally fixed and does not change over time. However, in practical engineering, the degradation trace of many products usually has significant degradation rate variation due to the variation of external operating conditions and internal mechanisms, and exhibits two-stage degradation phenomena such as high-performance capacitors, LCDs, liquid coupling devices, light emitting diodes, batteries, bearings, etc. as shown in FIG. 2. For such two-stage degraded products, the conventional life prediction method has low prediction accuracy.
In order to improve the prediction accuracy, many researchers have proposed a two-stage degradation model, but the existing method mainly focuses on whether the change points are random or whether individual differences are considered, and neglects the derivation of the analytical solution of degradation nonlinearity and lifetime of each stage. Each stage of the existing degradation model is a linear Wiener process, and as can be seen from fig. 2, in practice, the degradation process of each stage presents a nonlinear characteristic due to changes of load, internal state and external environment. Therefore, it is more reasonable to study the residual life prediction method of the two-stage degradation product based on the nonlinear Wiener process.
Disclosure of Invention
In view of the above, the present invention provides a new method for predicting remaining life of a two-stage degraded product, which makes up for the deficiencies of the prior art and can effectively predict the remaining life of the two-stage degraded product.
In order to achieve the purpose, the invention provides the following technical scheme:
a method for predicting the residual life of a two-stage degradation product comprises three contents of degradation process modeling, model parameter estimation and residual life prediction, and specifically comprises the following steps: the method comprises the following steps: modeling a degradation process, namely segmenting from a degradation variable point, and establishing a two-stage degradation model by utilizing a nonlinear Wiener process model; step two: estimating model parameters, namely estimating unknown parameters of the model by collecting historical degradation data and utilizing a two-stage unit MLE method; step three: predicting the residual life, namely deducing an approximate analytic solution of first-arrival time distribution by utilizing space-time transformation based on the two-stage degradation model obtained in the step one; step four: on the basis of the first arrival time distribution result obtained in the step three, considering the degradation rate of each stage as a random variable, and based on a total probability law, obtaining a first arrival time distribution function considering individual differences; step five: under the condition that the variable point does not appear, according to the characteristics of a Wiener process, a transition probability density function for transferring from an initial degradation state to a variable point degradation state can be deduced, and further, a state transition density function considering individual difference is obtained based on a total probability law; step six: based on the results obtained in the fourth step and the fifth step, based on the general probability law and the properties of Gaussian distribution, a life probability density function under the first arrival time concept of the two-stage degradation model can be deduced; and substituting the estimated value of the model parameter obtained in the second step into a life probability density function, thereby realizing the life prediction of the two-stage degraded product.
Further, in the first step, segmentation is carried out from a change point, each stage carries out degradation modeling by utilizing a nonlinear Wiener process, and the constructed two-stage degradation model is as follows:
Figure BDA0002320641730000021
wherein: x (t) represents the amount of degradation of the product, X (0) and X (τ) represent the state of degradation at the initial time and at the point of change, respectively,
Figure BDA0002320641730000022
and
Figure BDA0002320641730000023
is drift of each stageFunction, σ1And σ2Is the diffusion coefficient of each stage, B (t) is the standard Brownian motion;
further, in the second step, model unknown parameter estimation is performed by the proposed two-stage unit MLE method, which specifically comprises the following processes: firstly, historical degradation data of N degradation products are collected, and a degradation increment training data set delta X is establishednWherein m isnRepresenting the measured number of the nth product, Δ tn,j-1=tn,j-tn,j-1(j=1,..,mn) Indicating the detection interval.
Figure BDA0002320641730000024
In the first stage of parameter estimation, the degradation model parameters are estimated by using the degradation data of each test sample
Figure BDA0002320641730000025
According to the nature of the standard Brownian motion, Δ XnFollowing a normal distribution, the log-likelihood function for the parameter Θ is expressed as:
Figure BDA0002320641730000026
wherein: m isaAnd mbRepresenting the number of measurements in the first and second phases, satisfying ma+mb=mn
Figure BDA0002320641730000027
Figure BDA0002320641730000028
Figure BDA0002320641730000029
Correlating α to log likelihood function1,n,α2,n(n=1,...,N),
Figure BDA0002320641730000031
And
Figure BDA0002320641730000032
so that it is equal to 0, the corresponding parameter estimate is obtained:
Figure BDA0002320641730000033
Figure BDA0002320641730000034
wherein, firstly, the 'FMINSEARCH' in MATLAB is utilized to search the parameter β1And β2The estimated value is obtained by substituting the formula. Then, the maximum log-likelihood function is obtained, and the value of the variable point can be obtained.
Figure BDA0002320641730000035
In a second stage of parameter estimation, using each sample estimated in the first stage
Figure BDA0002320641730000036
And
Figure BDA0002320641730000037
and
Figure BDA0002320641730000038
and carrying out statistical analysis to obtain the hyper-parameters corresponding to the random parameter distribution.
Further, in the third step, based on the space-time transformation, referring to the lifetime derivation process of the nonlinear Wiener process, the first-arrival time distribution function is expressed as:
Figure BDA0002320641730000039
wherein: f. ofT(t | τ) is the lifetime probability density function at time t, xτIs the amount of degradation at the point of change;
further, in step four, based on the result of step three, considering individual differences, the degradation rate of each stage is regarded as a random variable, i.e. set
Figure BDA00023206417300000310
And
Figure BDA00023206417300000311
according to the law of general probability, the first-arrival time distribution function considering individual variability is expressed as:
(1)0<t≤τ
Figure BDA00023206417300000312
(2)t>τ
Figure BDA00023206417300000313
further, in the fifth step, if the change point does not appear, it is difficult to accurately acquire the degradation state at the change point. The occurrence of the second stage degradation process means that the first stage degradation has not reached the failure threshold, i.e., the lifetime is greater than the transition point time. Therefore, the state transition probability entering the second stage is defined as:
mτ(xτ)=Pr{X(τ)=xτ|X(0)=x0,T>τ}Pr{T>τ}
wherein: m isτ(xτ) Representing slave state x0Transition to State xτThe transition probability density of (2).
According to Wiener process characteristics and the law of total probability, the transition probability density function considering the influence of random effects can be expressed as:
Figure BDA0002320641730000041
further, in step six, based on the results of step four and step five, the lifetime probability density function under the first arrival time concept can be further expressed as:
Figure BDA0002320641730000042
if t is greater than 0 and less than or equal to tau, fTThe expression of (t | τ) is the same as the result of step four, if t > τ, based on the nature of gaussian distribution, the integral in the above equation can be solved, and further the lifetime probability density function can be obtained as follows:
Figure BDA0002320641730000043
wherein: a ═ A1-A2,B=B1-B2
Figure BDA0002320641730000044
Figure BDA0002320641730000045
Figure BDA0002320641730000046
Figure BDA0002320641730000047
Wherein: phi (-) and phi (-) represent the cumulative distribution function and probability density function of a standard normal distribution, respectively,
Figure BDA0002320641730000048
Figure BDA0002320641730000051
Figure BDA0002320641730000052
Figure BDA0002320641730000053
Figure BDA0002320641730000054
and substituting the parameters estimated in the second step into the life probability density function obtained in the sixth step to realize the residual life prediction of the two-stage product.
The core idea and principle of the invention are as follows:
the invention provides a two-stage degradation model based on a nonlinear Wiener process, which carries out model unknown parameter estimation by a two-stage unit MLE method, considers the difference of individuals and the uncertainty of a variable point degradation state, and deduces an approximate analytical solution of the residual life under the first arrival time concept.
The invention has the beneficial effects that:
1) compared with the existing two-stage degradation model, the two-stage degradation model is constructed by adopting a nonlinear Wiener process, so that the two-stage degradation process is more reasonably described; 2) the model is subjected to parameter estimation by a two-stage unit MLE method, and the method can independently and effectively estimate the parameters of each sample and is not limited by random parameter distribution; 3) through space-time transformation, an approximate analytic solution of first-arrival time distribution of the two-stage degradation model is obtained, and the service life of the two-stage degradation model can be effectively represented; 4) the uncertainty of the degradation state at the variable point and the difference between individuals are comprehensively considered, and the accuracy of life prediction is improved.
Drawings
In order to make the purpose, technical scheme and beneficial effect of the invention more clear, the invention provides the following drawings for explanation:
FIG. 1 is a flow chart of the residual life prediction for a two-stage degraded product of the present invention;
FIG. 2 is a schematic diagram of the degradation curve of the two-stage degraded product of the present invention;
FIG. 3 is a graph of simulated degradation trajectories;
FIG. 4 is a schematic illustration of remaining life prediction;
FIG. 5 is a schematic diagram of a predicted expected and actual remaining life comparison.
Detailed Description
Preferred embodiments of the present invention will be described in detail below with reference to the accompanying drawings. It should be understood that the preferred embodiments are illustrative of the invention only and are not limiting upon the scope of the invention.
The implementation process of this embodiment specifically includes: the model is established in the degradation process, the model parameter estimation and the residual life prediction are carried out, a numerical simulation is used for explanation based on an MATLAB tool, and the effect of the invention is shown by combining with the attached drawings. Fig. 1 is a flow chart of life prediction of a two-stage degraded product, and as shown in the figure, the method specifically includes the following steps:
the method comprises the following steps: and modeling a degradation process. Fig. 2 is a schematic diagram of the degradation of a two-stage degradation product, and it can be seen from the diagram that the degradation process of the product shows a two-stage degradation phenomenon. Thus, segmentation is performed at the degradation transition point, utilizing two non-linear Wi with different degradation rateseneAnd (3) carrying out degradation modeling on the r process, wherein the constructed two-stage degradation model is as follows:
Figure BDA0002320641730000061
wherein: x (t) represents the amount of degradation of the product, X (0) and X (τ) represent the state of degradation at the initial time and at the point of change, respectively,
Figure BDA0002320641730000062
and
Figure BDA0002320641730000063
is a drift function of each stage, σ1And σ2Is the diffusion coefficient of each stage, and B (t) is the standard Brownian motion. To express the difference between individuals, let
Figure BDA0002320641730000064
And
Figure BDA0002320641730000065
but these values are fixed values for each sample.
Then, the lifetime under the first arrival time concept can be expressed as:
T=inf{t:X(t)≥w|X(0)≤w}
wherein: w is a failure threshold, the value of which is determined by relevant industry standards or expert experience, and is typically a constant.
Step two: and estimating model parameters. A degradation trace for a number of sample samples was first generated using MATLAB, as shown in fig. 3. Constructing a degradation incremental training set for each group of degradation data
Figure BDA0002320641730000066
In the first stage of parameter estimation, the degradation model parameters are estimated by using the degradation data of each test sample
Figure BDA0002320641730000067
According to the nature of the standard Brownian motion, Δ XnAnd obtaining a log-likelihood function of the parameter theta according to normal distribution:
Figure BDA0002320641730000068
wherein: m isaAnd mbRepresenting the number of measurements in the first and second phases, satisfying ma+mb=mn
Figure BDA0002320641730000069
Figure BDA00023206417300000610
Figure BDA00023206417300000611
Correlating α to log likelihood function1,n,α2,n(n=1,...,N),
Figure BDA00023206417300000612
And
Figure BDA00023206417300000613
so that it is equal to 0, the corresponding parameter estimate is obtained:
Figure BDA0002320641730000071
Figure BDA0002320641730000072
wherein, firstly, the 'FMINSEARCH' in MATLAB is utilized to search the parameter β1And β2The estimated value is obtained by substituting the formula. Then, the maximum log-likelihood function is obtained, and the value of the variable point can be obtained.
Figure BDA0002320641730000073
This results in each sample transition point. In a second stage of parameter estimation, using each sample estimated in the first stage
Figure BDA0002320641730000074
And
Figure BDA0002320641730000075
and
Figure BDA0002320641730000076
statistical analysis was performed, and assuming that the change point obeyed a normal distribution, the mean parameter of the normal distribution obeyed the change point was 49.93, and the variance was 0.98.
Step three: and predicting the residual life. Referring to the lifetime derivation process of the nonlinear Wiener process, if the degradation process is:
Figure BDA00023206417300000714
based on the space-time transformation, the following can be obtained:
Figure BDA0002320641730000077
wherein:
Figure BDA0002320641730000078
similarly, it can be obtained that the first-arrival time distribution function of the two-stage degradation model is expressed as:
Figure BDA0002320641730000079
wherein: f. ofT(t | τ) is the lifetime probability density function at time t, xτIs the amount of degradation at the point of change.
Step four: considering individual differences, the degradation rate of each stage is regarded as a random variable, i.e. set
Figure BDA00023206417300000710
And
Figure BDA00023206417300000711
according to the following theorem:
Figure BDA00023206417300000712
based on the law of total probability, the first-arrival time distribution function considering individual differences is expressed as:
(1)0<t≤τ
Figure BDA00023206417300000713
(2)t>τ
Figure BDA0002320641730000081
step five: if the change point does not appear, it is difficult to accurately acquire the degradation state at the change point. The occurrence of the second stage degradation process means that the first stage degradation has not reached the failure threshold, i.e., the lifetime is greater than the transition point time. Therefore, the state transition probability entering the second stage is defined as:
mτ(xτ)=Pr{X(τ)=xτ|X(0)=x0,T>τ}Pr{T>τ}
wherein: m isτ(xτ) Representing slave state x0Transition to State xτThe transition probability density of (2).
According to Wiener process characteristics and the law of total probability, the transition probability density function considering the influence of random effects can be expressed as:
Figure BDA0002320641730000082
step six: based on the results of step four and step five, the lifetime probability density function under the first arrival time concept can be further expressed as:
Figure BDA0002320641730000083
if t is greater than 0 and less than or equal to tau, fTThe expression of (t | τ) is the same as the result of step four, if t > τ, based on the nature of gaussian distribution, the integral in the above equation can be solved, and further the lifetime probability density function can be obtained as follows:
Figure BDA0002320641730000084
wherein: a ═ A1-A2,B=B1-B2
Figure BDA0002320641730000085
Figure BDA0002320641730000086
Figure BDA0002320641730000087
Figure BDA0002320641730000091
Figure BDA0002320641730000092
Wherein: phi (-) and phi (-) represent the cumulative distribution function and probability density function of a standard normal distribution, respectively,
Figure BDA0002320641730000093
Figure BDA0002320641730000094
Figure BDA0002320641730000095
Figure BDA0002320641730000096
Figure BDA0002320641730000097
substituting the parameters estimated in the second step into the life probability density function obtained in the sixth step to obtain a remaining life probability density function of the two-stage product, as shown in fig. 4, and fig. 5 is a comparison of the predicted expectation of the remaining life obtained by the method of the present invention with the actual remaining life, and results of the single-stage nonlinear degradation model and the two-stage linear degradation model. As can be seen from the life prediction result, the method can effectively predict the residual life of the two-stage degraded product.
Finally, the above embodiments are only intended to illustrate the technical solutions of the present invention and not to limit the present invention, and although the present invention has been described in detail with reference to the preferred embodiments, it will be understood by those skilled in the art that modifications or equivalent substitutions may be made on the technical solutions of the present invention without departing from the spirit and scope of the technical solutions, and all of them should be covered by the claims of the present invention.

Claims (7)

1. The method for predicting the residual life of the two-stage degraded product is characterized by comprising the following steps: the method comprises the following steps:
the method comprises the following steps: modeling a degradation process, namely segmenting from a degradation variable point, and establishing a two-stage degradation model by utilizing a nonlinear Wiener process model;
step two: estimating model parameters, namely estimating unknown parameters of the model by collecting historical degradation data and utilizing a two-stage unit MLE method;
step three: predicting the residual life, namely deducing an approximate analytic solution of first-arrival time distribution by utilizing space-time transformation based on the two-stage degradation model obtained in the step one;
step four: on the basis of the first arrival time distribution result obtained in the step three, considering the degradation rate of each stage as a random variable, and obtaining a first arrival time distribution function considering the individual difference on the basis of a total probability law;
step five: under the condition that the change point does not appear, deducing a transition probability density function for transferring from the initial degradation state to the change point degradation state according to the characteristics of a Wiener process, and further obtaining a state transition density function considering individual difference based on a total probability law;
step six: deducing a life probability density function under the first arrival time concept of the two-stage degradation model based on the total probability law and the properties of Gaussian distribution based on the results obtained in the fourth step and the fifth step; and substituting the estimated value of the model parameter obtained in the second step into a life probability density function, thereby realizing the life prediction of the two-stage degraded product.
2. The method of predicting remaining life of a two-stage degraded product according to claim 1, characterized in that: in the first step, segmentation is carried out from a change point, each stage carries out degradation modeling by utilizing a nonlinear Wiener process, and the constructed two-stage degradation model is as follows:
Figure FDA0002320641720000011
wherein: x (t) represents the amount of degradation of the product, X (0) and X (τ) represent the state of degradation at the initial time and at the point of change, respectively,
Figure FDA0002320641720000012
and
Figure FDA0002320641720000013
is a drift function of each stage, σ1And σ2Is the diffusion coefficient of each stage, and B (t) is the standard Brownian motion.
3. The method of predicting remaining life of a two-stage degraded product according to claim 1, characterized in that: in the second step, model unknown parameter estimation is carried out by the proposed two-stage unit MLE method, and the method specifically comprises the following processes:
firstly, historical degradation data of N degradation products are collected, and a degradation increment training data set delta X is establishednWherein m isnRepresenting the measured number of the nth product, Δ tn,j-1=tn,j-tn,j-1(j=1,...,mn) Indicating a detection interval;
Figure FDA0002320641720000014
in the first stage of parameter estimation, the degradation model parameters are estimated by using the degradation data of each test sample
Figure FDA0002320641720000021
According to the nature of the standard Brownian motion, Δ XnFollowing a normal distribution, the log-likelihood function for the parameter Θ is expressed as:
Figure FDA0002320641720000022
wherein: m isaAnd mbRepresenting the number of measurements in the first and second phases, satisfying ma+mb=mn
Figure FDA0002320641720000023
Figure FDA0002320641720000024
Figure FDA0002320641720000025
Correlating α to log likelihood function1,n,α2,n(n=1,...,N),
Figure FDA0002320641720000026
And
Figure FDA0002320641720000027
so that it is equal to 0, the corresponding parameter estimate is obtained:
Figure FDA0002320641720000028
Figure FDA0002320641720000029
wherein, firstly, the 'FMINSEARCH' in MATLAB is utilized to search the parameter β1And β2Substituting the formula into the formula to obtain an estimated value; then, the maximum log-likelihood function is solved, and the value of the variable point can be obtained;
Figure FDA00023206417200000210
in a second stage of parameter estimation, using each sample estimated in the first stage
Figure FDA00023206417200000211
And
Figure FDA00023206417200000212
and
Figure FDA00023206417200000213
and carrying out statistical analysis to obtain the hyper-parameters corresponding to the random parameter distribution.
4. The method of predicting remaining life of a two-stage degraded product according to claim 1, characterized in that: in the third step, based on space-time transformation, referring to the lifetime derivation process of the nonlinear Wiener process, the first-arrival time distribution function is expressed as:
Figure FDA00023206417200000214
wherein: f. ofT(t | τ) is the lifetime probability density function at time t, xτIs the amount of degradation at the point of change.
5. The method of predicting remaining life of a two-stage degraded product according to claim 1, characterized in that: in step four, based on the result of step three, the individual differences are considered, and the degradation rate of each stage is taken as a random variable, i.e. the degradation rate is set
Figure FDA00023206417200000215
And
Figure FDA00023206417200000216
according to the law of general probability, the first-arrival time distribution function considering individual variability is expressed as:
when t is more than 0 and less than or equal to tau:
Figure FDA0002320641720000031
when t > τ:
Figure FDA0002320641720000032
6. the method of predicting remaining life of a two-stage degraded product according to claim 1, characterized in that: in step five, the state transition probability entering the second stage is defined as:
mτ(xτ)=Pr{X(τ)=xτ|X(0)=x0,T>τ}Pr{T>τ}
wherein: m isτ(xτ) Representing slave state x0Transition to State xτA transition probability density of (a);
according to Wiener process characteristics and the law of total probability, the transition probability density function considering the influence of random effects can be expressed as:
Figure FDA0002320641720000033
7. the method of predicting remaining life of a two-stage degraded product according to claim 1, characterized in that: in step six, based on the results of step four and step five, the lifetime probability density function under the first arrival time concept can be further expressed as:
Figure FDA0002320641720000034
if t is greater than 0 and less than or equal to tau, fTThe expression of (t | τ) is the same as the result of step four, if t > τ, based on the nature of gaussian distribution, the integral in the above equation can be solved, and further the lifetime probability density function can be obtained as follows:
Figure FDA0002320641720000035
wherein: a ═ A1-A2,B=B1-B2
Figure FDA0002320641720000036
Figure FDA0002320641720000041
Figure FDA0002320641720000042
Figure FDA0002320641720000043
Wherein phi (-) and phi (-) respectively represent a cumulative distribution function and a probability density function of a standard normal distribution,
Figure FDA0002320641720000044
Figure FDA0002320641720000045
Figure FDA0002320641720000046
Figure FDA0002320641720000047
Figure FDA0002320641720000048
and substituting the parameters estimated in the second step into the life probability density function obtained in the sixth step to realize the residual life prediction of the two-stage product.
CN201911296295.0A 2019-12-16 2019-12-16 Residual life prediction method for two-stage degraded product Active CN111046564B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911296295.0A CN111046564B (en) 2019-12-16 2019-12-16 Residual life prediction method for two-stage degraded product

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911296295.0A CN111046564B (en) 2019-12-16 2019-12-16 Residual life prediction method for two-stage degraded product

Publications (2)

Publication Number Publication Date
CN111046564A true CN111046564A (en) 2020-04-21
CN111046564B CN111046564B (en) 2024-04-05

Family

ID=70236823

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911296295.0A Active CN111046564B (en) 2019-12-16 2019-12-16 Residual life prediction method for two-stage degraded product

Country Status (1)

Country Link
CN (1) CN111046564B (en)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111523251A (en) * 2020-06-09 2020-08-11 江苏科技大学 Method for rapidly evaluating service life of product under random environmental stress
CN111753416A (en) * 2020-06-17 2020-10-09 重庆大学 Lithium ion battery RUL prediction method based on two-stage Wiener process
CN111814331A (en) * 2020-07-08 2020-10-23 重庆大学 Method for predicting residual service life of equipment with multiple variable points under complex conditions
CN111967640A (en) * 2020-07-10 2020-11-20 西北工业大学 Tool residual life prediction method considering tool wear amount and surface roughness
CN111967133A (en) * 2020-07-10 2020-11-20 西北工业大学 Method for predicting residual life of cutter in multiple cutting stages
CN112396215A (en) * 2020-10-22 2021-02-23 国网浙江省电力有限公司嘉兴供电公司 Intelligent prediction method for self-adaptive interval of residual life of equipment
CN112505568A (en) * 2020-11-11 2021-03-16 电子科技大学 Multi-stack solid oxide fuel cell system service life prediction method
CN112765813A (en) * 2021-01-19 2021-05-07 中国人民解放军火箭军工程大学 Method for predicting residual life of equipment under sequential Bayesian framework
CN112949058A (en) * 2021-02-27 2021-06-11 中国人民解放军火箭军工程大学 Heuristic residual life prediction method for lithium battery based on implicit nonlinear wiener process
CN113033015A (en) * 2021-04-09 2021-06-25 中国人民解放军火箭军工程大学 Degraded equipment residual life prediction method considering two-stage self-adaptive Wiener process
CN113378309A (en) * 2021-01-28 2021-09-10 河南科技大学 Rolling bearing health state online monitoring and residual life prediction method
CN114091790A (en) * 2022-01-20 2022-02-25 浙江大学 Life prediction method fusing field data and two-stage accelerated degradation data
CN114169107A (en) * 2021-12-14 2022-03-11 大连理工大学 Service life prediction algorithm based on random stagnation nonlinear degradation model
CN116227366A (en) * 2023-05-08 2023-06-06 浙江大学 Two-stage motor insulation life prediction method
CN116756505A (en) * 2023-06-07 2023-09-15 上海正泰电源系统有限公司 Photovoltaic equipment intelligent management system and method based on big data

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107145720A (en) * 2017-04-19 2017-09-08 浙江大学 It is continuous to degenerate and the unknown equipment method for predicting residual useful life impacted under collective effect
CN109387779A (en) * 2018-09-28 2019-02-26 河北工业大学 A kind of omnipotent breaker operation attachment method for predicting residual useful life based on statistical data driving
CN109657937A (en) * 2018-11-30 2019-04-19 西安电子科技大学 A kind of Reliability Assessment and life-span prediction method based on degraded data
CN110197288A (en) * 2019-05-30 2019-09-03 重庆大学 The remaining life prediction technique of equipment under the influence of failure

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107145720A (en) * 2017-04-19 2017-09-08 浙江大学 It is continuous to degenerate and the unknown equipment method for predicting residual useful life impacted under collective effect
CN109387779A (en) * 2018-09-28 2019-02-26 河北工业大学 A kind of omnipotent breaker operation attachment method for predicting residual useful life based on statistical data driving
CN109657937A (en) * 2018-11-30 2019-04-19 西安电子科技大学 A kind of Reliability Assessment and life-span prediction method based on degraded data
CN110197288A (en) * 2019-05-30 2019-09-03 重庆大学 The remaining life prediction technique of equipment under the influence of failure

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
张正新;胡昌华;高迎彬;陈墨;: "多阶段随机退化设备剩余寿命预测方法", 系统工程学报, no. 01, pages 3 - 9 *
张鹏 等: "考虑随机效应的两阶段退化系统剩余寿命预测方法", 中国测试, vol. 45, no. 1, pages 1 - 7 *
牛一凡;邵景峰;: "基于非线性数据融合的设备多阶段寿命预测", 信息与控制, no. 06, pages 93 - 101 *
黄亮;刘君强;贡英杰;: "基于Wiener过程的发动机多阶段剩余寿命预测", 北京航空航天大学学报, no. 05, pages 190 - 196 *

Cited By (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111523251A (en) * 2020-06-09 2020-08-11 江苏科技大学 Method for rapidly evaluating service life of product under random environmental stress
CN111523251B (en) * 2020-06-09 2023-04-21 江苏科技大学 Method for rapidly evaluating service life of product under random environmental stress
CN111753416A (en) * 2020-06-17 2020-10-09 重庆大学 Lithium ion battery RUL prediction method based on two-stage Wiener process
CN111814331A (en) * 2020-07-08 2020-10-23 重庆大学 Method for predicting residual service life of equipment with multiple variable points under complex conditions
CN111814331B (en) * 2020-07-08 2023-10-20 重庆大学 Multi-point equipment residual service life prediction method under complex condition
CN111967640A (en) * 2020-07-10 2020-11-20 西北工业大学 Tool residual life prediction method considering tool wear amount and surface roughness
CN111967133A (en) * 2020-07-10 2020-11-20 西北工业大学 Method for predicting residual life of cutter in multiple cutting stages
CN111967640B (en) * 2020-07-10 2023-05-23 西北工业大学 Cutter residual life prediction method considering cutter abrasion loss and surface roughness
CN112396215B (en) * 2020-10-22 2022-06-17 国网浙江省电力有限公司嘉兴供电公司 Intelligent prediction method for self-adaptive interval of residual life of equipment
CN112396215A (en) * 2020-10-22 2021-02-23 国网浙江省电力有限公司嘉兴供电公司 Intelligent prediction method for self-adaptive interval of residual life of equipment
CN112505568A (en) * 2020-11-11 2021-03-16 电子科技大学 Multi-stack solid oxide fuel cell system service life prediction method
CN112505568B (en) * 2020-11-11 2022-03-15 电子科技大学 Multi-stack solid oxide fuel cell system service life prediction method
CN112765813A (en) * 2021-01-19 2021-05-07 中国人民解放军火箭军工程大学 Method for predicting residual life of equipment under sequential Bayesian framework
CN113378309A (en) * 2021-01-28 2021-09-10 河南科技大学 Rolling bearing health state online monitoring and residual life prediction method
CN113378309B (en) * 2021-01-28 2023-04-14 河南科技大学 Rolling bearing health state online monitoring and residual life prediction method
CN112949058A (en) * 2021-02-27 2021-06-11 中国人民解放军火箭军工程大学 Heuristic residual life prediction method for lithium battery based on implicit nonlinear wiener process
CN113033015B (en) * 2021-04-09 2024-05-14 中国人民解放军火箭军工程大学 Degradation equipment residual life prediction method considering two-stage self-adaptive Wiener process
CN113033015A (en) * 2021-04-09 2021-06-25 中国人民解放军火箭军工程大学 Degraded equipment residual life prediction method considering two-stage self-adaptive Wiener process
CN114169107A (en) * 2021-12-14 2022-03-11 大连理工大学 Service life prediction algorithm based on random stagnation nonlinear degradation model
CN114169107B (en) * 2021-12-14 2024-04-16 大连理工大学 Life prediction algorithm based on random stagnation nonlinear degradation model
CN114091790B (en) * 2022-01-20 2022-05-03 浙江大学 Life prediction method fusing field data and two-stage accelerated degradation data
CN114091790A (en) * 2022-01-20 2022-02-25 浙江大学 Life prediction method fusing field data and two-stage accelerated degradation data
CN116227366B (en) * 2023-05-08 2023-08-11 浙江大学 Two-stage motor insulation life prediction method
CN116227366A (en) * 2023-05-08 2023-06-06 浙江大学 Two-stage motor insulation life prediction method
CN116756505A (en) * 2023-06-07 2023-09-15 上海正泰电源系统有限公司 Photovoltaic equipment intelligent management system and method based on big data
CN116756505B (en) * 2023-06-07 2024-05-28 上海正泰电源系统有限公司 Photovoltaic equipment intelligent management system and method based on big data

Also Published As

Publication number Publication date
CN111046564B (en) 2024-04-05

Similar Documents

Publication Publication Date Title
CN111046564A (en) Method for predicting residual life of two-stage degraded product
WO2019174142A1 (en) Multi-mode degradation process modelling and remaining service life prediction method
CN102789545B (en) Based on the Forecasting Methodology of the turbine engine residual life of degradation model coupling
CN102042848B (en) Prediction method of multi-functional parameter accelerated degradation testing product life based on multivariate hybrid time sequence analysis
CN110851980A (en) Method and system for predicting residual life of equipment
CN102252843B (en) Assessment method for rolling bearing performance variation
CN110398368B (en) FBM-based bearing inner ring fault residual life prediction method based on long correlation model
CN109813544B (en) A kind of rotating machinery Incipient Fault Diagnosis method and system based on on-line monitoring
CN102629300A (en) Step stress accelerated degradation data assessment method based on gray prediction models
Sreenuch et al. Probabilistic Monte-Carlo method for modelling and prediction of electronics component life
CN106529832A (en) Relay protection system risk assessment method based on Markov reliability correction model
Zhou et al. A control chart based on likelihood ratio test for detecting patterned mean and variance shifts
CN110852509A (en) Fault prediction method and device of IGBT module and storage medium
CN112835769A (en) Service data abnormity diagnosis method, device, equipment and storage medium
CN104463331A (en) Accelerated degradation experiment modeling method based on fuzzy theory
CN112966441A (en) Equipment residual life evaluation method based on continuous Weiner process damage
CN111523727B (en) Method for predicting remaining life of battery by considering recovery effect based on uncertain process
CN114036647A (en) Power battery safety risk assessment method based on real vehicle data
CN117471346A (en) Method and system for determining remaining life and health status of retired battery module
Liao et al. Nonparametric and semi-parametric sensor recovery in multichannel condition monitoring systems
CN113326744B (en) Method and system for detecting on-orbit state abnormity of spacecraft
CN110532629A (en) A kind of method for predicting residual useful life of multi-modal degenerative process
CN105741184B (en) Transformer state evaluation method and device
CN111079270A (en) Bearing residual life prediction method based on binary mixed random process
CN114492507B (en) Bearing residual life prediction method under digital-analog cooperative driving

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