Zhonglai Wang1 School of Mechatronics Engineering, University of Electronic Science and A Simulation Method to Estimate Technology of China, Chengdu, Sichuan 611731, China; The State Key Laboratory of Advanced Design Two Types of Time-Varying and Manufacturing for Vehicle Body, Failure Rate of Dynamic Changsha 410082, China e-mail: [email protected] Systems Xiaoqiang Zhang School of Mechatronics Engineering, The failure rate of dynamic systems with random parameters is time-varying even for lin- University of Electronic Science and ear systems excited by a stationary random input. In this paper, we propose a simulation- Technology of China, based method to estimate two types (type I and type II) of time-varying failure rate of Chengdu, Sichuan 611731, China dynamic systems. The input stochastic processes are discretized in time and the trajecto- e-mail: [email protected] ries of the output stochastic process are calculated. The time of interest is partitioned into a series of time intervals and the saddlepoint approximation (SPA) is employed to Hong-Zhong Huang estimate the probability of failure in each interval. Type I follows the commonly used def- School of Mechatronics Engineering, inition of failure rate. It is estimated at discrete time intervals using SPA and the correla- University of Electronic Science and tion information from a properly selected time-dependent copula function. Type II is a Technology of China, proposed new concept of time-varying failure rate. It provides a way to predict the failure Chengdu, Sichuan 611731, China rate considering a virtual “good-as-old” repair action of repairable dynamic systems. e-mail: [email protected] The effectiveness of the proposed method is illustrated with a vehicle vibration example. [DOI: 10.1115/1.4034300] Zissimos P. Mourelatos Mechanical Engineering Department, Keywords: time-varying failure rate, dynamic systems, simulation-based method, Oakland University, repaired samples, saddlepoint approximation, time-dependent copula function Rochester, MI 48309 e-mail: [email protected]
1 Introduction further studied for time-dependent reliability analysis based on total probability theorem [8]. Also, a time-dependent reliability The failure rate quantifies the probability of failure over a small analysis is performed in Ref. [9] using a nested kriging based pre- time interval dt after time t under the condition that the system did diction model where time-independent reliability methods such as not fail before time t. In reliability engineering, it is common to the first order reliability method (FORM) [10] and the saddlepoint derive the failure rate from the cumulative distribution function approximation method (SPA) [11] can be used. The out-crossing (CDF) [1–4]. Sufficient data samples are needed to estimate the rate method has been widely used for time-dependent reliability CDF accurately with statistical methods, especially for a small problems and many advanced methods have been proposed probability of failure. For example, about 234,000 data samples [12–15] under the assumption of independent out-crossings fol- are needed to estimate a 0.00171 probability of failure with 10% lowing the Poisson distribution, which may result in a poor accu- error under 95% confidence. Estimation of the CDF using physics racy. A parallel system reliability formulation, the so-called PHI2 of failure (POF) has attracted more attention for structural systems method, using the out-crossing rate is proposed in Ref. [16]. A by considering the time-dependent properties. The time horizon is Monte Carlo based set theory method, similar with the PHI2 partitioned into many time intervals for discrete CDF estimation method, is reported in Ref. [17]. Studies reveal that the PHI2 and then a continuous CDF is fitted. Accordingly, discrete and based method shows relatively poor accuracy when dealing with continuous time-varying failure rates could be derived from the nonmonotonic problems [6,18]. In order to avoid the calculation discrete and continuous CDFs, respectively. of out-crossing rate, a time-dependent reliability analysis method To estimate a time-varying failure rate, time-dependent reliabil- based on stochastic process discretization is presented in ity analysis must be conducted over the planning horizon. Many Ref. [19]. methods have been proposed for time-dependent reliability analy- Most of time-dependent reliability analysis methods provide sis of structural systems. Because of the independent non-negative only the reliability over a specific time interval. Design under increments of the Gamma process, a combination of Gamma pro- uncertainty considering maintenance and lifecycle cost requires cess and statistics of live load models is used in Ref. [5]. A distri- however, the reliability at all times within the planning horizon. bution of extreme values approach is employed in Ref. [6]to Singh et al. [20] proposes such a method by estimating the failure transform the time-dependent reliability problem to a time- rate of random dynamic systems using importance sampling. independent reliability problem. An equivalent time-invariant Their method is suitable for low dimensional limit state functions. composite limit state is defined and used in Ref. [7] for time- Wang et al. [21] also presents a method to estimate the failure rate dependent reliability analysis. The composite limit state method is at all times within the planning horizon using an improved subset simulation with splitting by partitioning the original high dimen- 1Corresponding author. sional random process into a series of correlated, short duration, Contributed by the Design Automation Committee of ASME for publication in low dimensional random processes. However in estimating the the JOURNAL OF MECHANICAL DESIGN. Manuscript received January 22, 2016; final manuscript received July 5, 2016; published online September 14, 2016. Assoc. time-varying failure rate, the high computational effort is still a Editor: Xiaoping Du. challenge.
Journal of Mechanical Design Copyright VC 2016 by ASME DECEMBER 2016, Vol. 138 / 121404-1
Downloaded From: http://mechanicaldesign.asmedigitalcollection.asme.org/ on 10/05/2016 Terms of Use: http://www.asme.org/about-asme/terms-of-use Copula functions have been widely used for modeling the cor- provides a new way to predict the failure rate by essentially build- relation of random variables and time series in economic field ing a relationship between the system failure rate and the capacity [22,23]. Because of their great capacity to represent correlation, for maintenance of repairable dynamic systems. Copulas have been recently introduced in the structural safety The contributions of this article are: (1) a new time-varying field. Noh et al. [24] uses Gaussian Copula to describe the correla- failure rate (type I) estimation method is proposed for dynamic tion of input random variables for reliability-based design optimi- systems; (2) a time-dependent copula function approach is pro- zation and further identifies the joint CDFs from provided data posed to establish the failure correlation as time progresses with- using a copula function [25]. In order to express the multidimen- out additional computational effort; and (3) a new concept of sional correlation, the vine-copula has been also introduced in time-varying failure rate (type II) is presented for repaired sam- structural reliability analysis [26]. Tang et al. [27] studies the ples and its expression is derived. effect of different copula models on the structural system reliabil- The reminder of the paper is organized as follows. Section 2 ity results. reviews the SPA method and the method to construct a time- The above studies have played an important role in fostering dependent copula function. We use both methods in the proposed applications of Copula functions for time-independent problems time-varying failure rate approach. Section 3 provides details on in the structural safety field. For dynamic systems, the correlation the proposed approach and Sec. 4 uses a vehicle vibration exam- of performance between adjacent time intervals can be established ple to illustrate its effectiveness. Finally, Sec. 5 summarizes and accounting therefore, for time. In order to address the time- concludes. dependent probability of failure or failure rate, of dynamic systems, a simulation-based method is proposed in this article to 2 Saddlepoint Approximation and Copula Function estimate two types of time-varying failure rate. A simulation approach is used because the failed samples still 2.1 Saddlepoint Approximation. The SPA is an attractive remain in the operating queue, resulting in correlated failures in method to perform structural reliability analysis because of its time. The input stochastic processes are discretized using a small high computational accuracy in approximating the tails of a distri- time step into a set of high-dimensional random variables and the bution and its comparable efficiency with FORM [11,28]. If Y is a trajectories of the output stochastic process are calculated accu- random variable with probability density function (PDF) fY ðyÞ, the rately with simulation. The planning horizon is then partitioned moment generating function (MGF) of Y is given by ð into a series of time intervals and the saddlepoint approximations 1 nY ny (SPA) method is used to calculate the probability of failure in MðnÞ¼E½e ¼ e fY ðyÞdy (1) each time interval. With the same trajectories used in SPA, a 1 time-dependent copula is built to provide the correlation between the maximum response in each time interval and the maximum where E½ denotes expectation. The cumulative generating func- responses up to that time interval. A discrete mean time-varying tion (CGF) KðnÞ is failure rate in each time interval is calculated using an estimated probability of failure from the SPA and the correlation informa- KðnÞ¼log ½MðnÞ (2) tion from the estimated time-dependent copula. A continuous mean time-varying failure rate is also derived from the fitted CDF where log ½ is the natural logarithm. The inverse Fourier trans- from the discrete cumulative probabilities of failure based on a form is used to obtain fY ðyÞ from KðnÞ as validated Weibull distribution. This time-varying failure rate, ð i1 called Type I failure rate, is calculated using the classical 1 fY ðÞy ¼ exp½ KðÞn ny dn (3) definition. 2pi i1 During simulation, we assume that after failure, performance is restored to the state just before failure occurred (good-as-old Analytical forms of the CGF are only available for certain dis- repair assumption). In order to compare the failure evolution of tributions [29]. If there is no analytical form, an approximation the repaired and nonrepaired samples, we propose the concept of can be obtained using samples. It has been shown that an accurate time-varying failure rate for repaired samples, called type II fail- approximation of the CGF can be obtained using the following ure rate, derive its expression and compare it with type I. Type II four cumulants [30]
8 > s1 > k1 ¼ > N > s > 2 > Nss2 s1 <> k2 ¼ NsðÞNs 1 3 2 (4) > 2s 3Nss1s2 þ N s3 > k ¼ 1 s > 3 N ðÞN 1 ðÞN 2 > s s s > 4 2 2 2 > 6s1 þ 12Nss1s2 3NsðÞNs 1 s2 4NsðÞNs þ 1 s1s3 þ Ns ðÞNs þ 1 s4 : k4 ¼ NsðÞNs 1 ðÞNs 2 ðÞNs 3
P Ns j r j where sr ¼ j¼1ðy Þ ðr ¼ 1; 2; 3; 4Þ and y ðj ¼ 1; …; NsÞ is the If the derivative of KðnÞ with respect to n is set equal to y as output of the performance function. NsðNs 4Þ is the sample size. Using the above four cumulants, the CGF is approximated as K0ðnÞ¼y (6)
1 2 1 3 1 4 the solution to Eq. (6) is the saddlepoint n , which is also the KðÞn ¼ k1n þ k2n þ k3n þ k4n (5) sp 2 3! 4! extreme point of KðnÞ ny in Eq. (3).
121404-2 / Vol. 138, DECEMBER 2016 Transactions of the ASME
Downloaded From: http://mechanicaldesign.asmedigitalcollection.asme.org/ on 10/05/2016 Terms of Use: http://www.asme.org/about-asme/terms-of-use Daniels [31] uses a power series expansion to evaluate the inte- gral in Eq. (3) resulting in the following expression of the PDF of Y
"#1 1 2 ½ KðÞnsp nspy fY ðÞy ¼ } e (7) 2pK nsp
The saddlepoint approximation approach can also be used to Fig. 1 Schematic of the time-varying failure rate estimation process approximate the CDF FY ðyÞ. A reasonable and widely used for- mula is provided by Lugannini and Rice [32] 1 1 motion are discretized in time and expressed in a state-space FY ðÞy ¼ PrfgY y ¼ UðÞw þ /ðÞw (8) w v form. The discretized equations are time integrated using for example, a Runge–Kutta method or Newmark-beta method [41] where Uð Þ and /ð Þ are the CDF and PDF of a standard normal and are expressed as distribution, and w and v can be expressed by Xðt þ DtÞ¼hðXðtÞ; Z; MðtÞ; tÞ (13) 1=2 w ¼ sgnðnspÞf2½nspy KðnspÞ g (9) c where Xðt þ DtÞ2< is the vector of uncertain states xsðt þ DtÞ, s ¼ 1; 2; …; c at time t þ Dt, Dt is the integration time step and and hð Þ indicates a function of the arguments. Z 2