Developing highly efficient and accurate methodologies is a central challenge in the field of time-dependent reliability. In this paper, a time-dependent-reliability analysis method that couples series expansion with the equivalent-plane approach is proposed. The first-order reliability method (FORM) is applied to evaluate the reliability on a discretized time instant, and expansion optimal linear estimation (EOLE) represents the stochastic process by series expansion and constructs its correlated representation. The equivalent-plane approach (EPA) is used to aggregate the discretized failure events and to estimate the time-dependent failure probability over the considered period. Three numerical examples demonstrate that the proposed method achieves close agreement with Monte Carlo simulation benchmarks while requiring far fewer limit state evaluations.
Introduction
Uncertainty is ubiquitous in engineering systems due to inherent variability, limited knowledge, and modeling errors . To quantify stochastic structural safety under uncertainty, probability-based structural reliability commonly measures safety in terms of failure probability. The first-order reliability method (FORM) is widely adopted because it estimates small failure probabilities efficiently by linearizing the limit state function at the design point in a transformed standard normal space . When the limit state function is strongly nonlinear, the first-order approximation may lose accuracy, and higher-order approaches such as second- and even third-order reliability methods (SORM and TORM) can be used to improve accuracy at increased computational cost . Beyond these design-point-based approximations, methods for reliability analysis include cross-entropy-based adaptive importance sampling , Kriging-based active learning , sparse polynomial chaos expansions , and entropy-based measures . These developments are predominantly formulated for time-invariant reliability evaluation at fixed time instants. However, in time-dependent reliability problems, temporal dependence introduces strong correlation across the service period and increases the effective problem dimension. To address this challenge, numerous methods for time-dependent reliability analysis have been developed. These methods are commonly grouped into three classes: outcrossing-rate methods, sampling-based methods, and equivalent methods.
Outcrossing-rate methods, which are derived from the Rice formula , are widely used. However, their accuracy may deteriorate when multiple dependent upcrossings occur within the considered period. Among the outcrossing methods, a widely used refinement is the PHI2 method developed by and . By introducing bivariate correlation corrections for upcrossing events, PHI2 preserves the low computational cost of outcrossing-rate methods and typically reduces approximation error relative to classical formulations.
Sampling-based methods are conceptually straightforward. However, crude Monte Carlo simulation (MCS) requires very large sample sizes to accurately estimate small failure probabilities, especially with expensive limit state evaluations . As a result, the computational cost can be prohibitive. To alleviate this burden, mitigation strategies include an adaptive importance sampling approach , a first-order sampling approach , random field discretization , and a surrogate model approach . Within random field discretization, series expansion methods including the Karhunen–Loève (KL) expansion , expansion optimal linear estimation (EOLE) , and orthogonal series expansion (OSE) provide effective low-dimensional approximations. A comparative assessment of their accuracy and efficiency is provided by .
Equivalent methods are a key class in time-dependent reliability analysis. The main idea is to transform time-dependent problems into static ones. Equivalent methods include extreme-value methods and envelope function methods . In practice, efficient global optimization (EGO) and mixed EGO are often employed as search strategies to locate critical time instants on surrogate models and thereby predict time-dependent reliability . Envelope function methods estimate the time-dependent reliability by constructing response envelopes, which avoids per-instant analyses . Compared with outcrossing-rate methods, equivalent methods typically achieve higher accuracy when the response exhibits strongly correlated or multiple crossings while remaining far more efficient than direct sampling in many practical settings.
The equivalent-plane approach (EPA) was initially developed to assess the reliability of series and parallel systems . proposed a sequential compounding method and derived analytical expressions for combining series and parallel systems. Gong et al. transformed variables in the standard normal space into a high-dimensional space with the correlation of stochastic processes and obtained the time-dependent reliability by combining the reliability at each time instant.
This work develops a time-dependent reliability method that couples series expansion with the equivalent-plane approach. In Sect. , the problem formulation and notation for time-dependent reliability are reviewed. In Sect. , the reliability at each time instant is evaluated by FORM, EOLE provides a correlated representation of the stochastic process in the standard normal space, and EPA is then used to estimate the time-dependent failure probability. Numerical examples and the conclusions are presented in Sects. and , respectively.
Problem formulation of time-dependent reliability
In the limit state function G(t,X,Y(t)), t denotes time, X=[X1,…,Xq] are time-invariant random variables, and Y(t) is the stochastic process. Failure occurs when G(⋅)≤0. The reliability of a structure with the limit state function is considered within the period of ts,te.
At a time instant ti∈ts,te, the failure event can be expressed as
Ei=Gti,X,Y(ti)≤0,
and the corresponding failure probability is
pfti=PrEi=PrGti,X,Y(ti)≤0.
Over the period ts,te, the failure event is
E=∃t∈ts,te:Gt,X,Y(t)≤0.
The failure probability over ts,te is
pfts,te=PrE=Pr∃t∈ts,te:Gt,X,Y(t)≤0.
When the stochastic process changes slowly or when resistance decays monotonically relative to the load, the failure probability over ts,te can be approximately expressed by pfte. However, when material properties or external loads of a structure are time-dependent, the structural reliability over the period requires a time-dependent reliability method.
The series expansion and equivalent-plane approach
At a discretized time instant ti, the limit state function is written as Gti,X,Y(ti). For reliability evaluation at the fixed time instant ti, the process value Y(ti) is treated as an equivalent random variable so that the problem involves (q+1) random variables (X,Y(ti)). By applying the Nataf transformation to (X,Y(ti)), an independent standard normal vector Vi=[V1(i),…,Vq+1(i)]⊤∈Rq+1 is obtained, and the transformed limit state function is denoted by GV(ti,Vi). The reliability at ti is quantified by the FORM reliability index
βi=minv∈Rq+1‖v‖:GV(ti,v)=0.
The minimizer of Eq. (), denoted by Vi∗, is the most probable point (MPP) at ti, such that βi=‖Vi∗‖.
To account for the dependence of the stochastic process across different time instants, the stochastic process Y(t) can be approximately represented as an infinite linear combination of orthogonal functions . In this paper, Y(t) is approximated using EOLE, a truncated series expansion. Let m(t) and σ(t) denote the mean and standard deviation of Y(t). Select n time instants t1,…,tn in [ts,te], and define the correlation matrix C∈Rn×n with entries Cjk=ρY(tj,tk). Let ρY(t)=[ρY(t,t1),…,ρY(t,tn)]⊤. The truncated series expansion reads as
Y(t)≈Y^(t)=m(t)+σ(t)∑k=1Mξkλkϕk⊤ρY(t),
where ξk∼N(0,1) denotes independent standard normal variables, {(λk,ϕk)}k=1M denotes the M largest eigenvalues and the corresponding eigenvectors of C, and M is the truncation order. With n time instants fixed, the process Y(t) is parameterized by ξ=[ξ1,…,ξM]⊤.
In this work, the truncation order M is selected such that the maximum pointwise approximation error over [ts,te] does not exceed 1 %. The pointwise error is estimated by
ε(t)=1-∑k=1Mϕk⊤ρY(t)2λk.
Based on Eq. (), we define the standardized process as
η(t)=Y^(t)-m(t)σ(t)=∑k=1Mξkbk(t),bk(t)=ϕk⊤ρY(t)λk.
At ti, η(ti)=b(ti)⊤ξ, with b(ti)=[b1(ti),…,bM(ti)]⊤. Here, the last component Vq+1(i) corresponds to the standardized process value at ti. Accordingly, the process-related standard normal variable in the above FORM analysis is linked to the EOLE variables by
Vq+1(i)≈η(ti)=b(ti)⊤ξ.
With the EOLE relations in Eqs. ()–(), the process component at ti is represented in terms of the common EOLE variables ξ. Consequently, its contribution can be mapped onto the ξ subspace along b(ti) when constructing the (q+M)-dimensional direction vector. Therefore, the uncertainty at different time instants can be described in a unified (q+M)-dimensional independent standard normal space. Let U∼N(0,Iq+M) denote the corresponding standard normal vector. The equivalent design point associated with the FORM result at ti is defined as
Ui∗=V1(i)∗,…,Vq(i)∗,Vq+1(i)∗b(ti)⊤⊤∈Rq+M,
and the corresponding unit direction vector is defined as
αi=Ui∗‖Ui∗‖∈Rq+M.
For notational convenience, define β=[β1,…,βn]⊤ and A=[α1,…,αn]∈R(q+M)×n. The first-order linearization of the limit state at ti is
Zi=βi-αi⊤U,i=1,…,n,
such that failure at ti corresponds to {Zi≤0}.
According to first-order reliability theory, the failure probability at ti is
pf(ti)=1-Φ(βi),
where Φ(⋅) denotes the standard normal cumulative distribution function.
For any two time instants ti and tj, the correlation coefficient between the linearized limit state function Zi and Zj is
ρij=αi⊤αj.
From the definitions above, the correlation matrix of Z={Zi}i=1n is R:
R=A⊤A.
Equation () implies that the diagonal entries equal 1 and that the off-diagonal entries coincide with the pairwise correlations. In particular, (R)ij=ρij for i,j=1,…,n.
The time-dependent failure over ts,te can be modeled as a series system across the discretized time instants. Accordingly, the linearized responses at the time instants form the components of a multivariate normal distribution. The system failure probability is then given by
pfts,te=Pr⋃i=1n{Zi≤0}=1-∫-∞βφn(w;R)dw,
where φn(⋅;R) denotes the probability density function of the zero-mean multivariate normal distribution with correlation matrix R∈Rn×n, and the integration is with respect to w.
Multidimensional integrals in Eq. () suffer from the “curse of dimensionality” when the dimensionality increases. Some approximate methods, such as the Monte Carlo and quasi-Monte-Carlo methods, also have difficulty obtaining an acceptable accuracy in low-failure-probability problems. Thus, the EPA is adopted as an effective approximation of multidimensional integrals and is effective for series and parallel systems. As mentioned above, the time-dependent reliability problem is converted into a sequence of static reliability problems. The n linearized limit state responses Z={Zi}i=1n represent the failure events Ei={Zi≤0} at the n time instants. Selecting two instants ti and tj, the joint failure probability is
pf,ij=PrEi∪Ej=1-Φ2βi,βj;ρij,
where Φ2(⋅,⋅;ρ) denotes the binormal cumulative distribution function with mean zero, unit variances, and correlation ρij.
The equivalent reliability index of the joint failure event Ei∪Ej is
βij=-Φ-1pf,ij.
Then, the equivalent limit state function for the union event Ei∪Ej is defined as
Zij=βij-αij⊤U,
where αij∈Rq+M is the unit vector in the direction of ∇Uβij and is defined by
αij=∂βij∂U/∂βij∂U.
The contribution of the kth component Uk of U (k=1,…,q+M) to βij is
∂βij∂Uk=∂βij∂pf,ij∂pf,ij∂βi∂βi∂Uk+∂βij∂pf,ij∂pf,ij∂βj∂βj∂Uk.
Each factor in the first term on the right side of Eq. () reads as
∂βij∂pf,ij=-∂Φ-1pf,ij∂pf,ij=-1φβij∂pf,ij∂βi=-φ(βi)Φβj-ρijβi1-ρij2∂βi∂Uk=αk(i).
The second term in Eq. () is obtained analogously from Eq. () by exchanging i and j.
Equivalent-plane aggregation with correlation-based pair selection.
Failure events are combined pairwise in sequence until the n failure events E1,…,En with Ei={Zi≤0} are reduced to a single equivalent limit state Zend, as shown in Fig. , with the associated failure event Eend={Zend≤0}. The time-dependent failure probability over [ts,te] is approximated by
pf(ts,te)≈Pr(Eend)=1-Φ(βend).
Since each compounding step introduces an approximation, the final time-dependent failure probability may depend on the compounding order, and the approximation error can accumulate through subsequent steps . To reduce this order sensitivity, we repeatedly select the remaining pair with the largest correlation ρij, compound the selected pair into a new equivalent event, update the correlations involving the new event, and continue until the final equivalent event Eend is obtained.
Numerical examples and discussionSupported steel beamProblem definition and solutions
The supported steel beam shown in Fig. is the most classic model of time-dependent reliability . The length of the beam is L, and the width and height of the rectangular section are b0 and h0, respectively. The unit weight of the steel beam is ρs, and the force of the weight p=ρsb0h0 can be regarded to be a fixed uniform load acting on the beam. The time-dependent load is a concentrated force F(t) applied at the midspan. The steel beam is assumed to be complete at initial time t=0. The corrosion of the rectangular section is recognized as uniform, and the corrosion rate is denoted as k. The cross-section width and height with time are b(t)=b0-2kt and h(t)=h0-2kt. The yield strength of the material is σs.
The model of the supported steel beam.
The autocorrelation function of the Gaussian process F(t) is defined as
ρ(ti,tj)=exp-|tj-ti|ℓc2,
where the correlation coefficient depends on the time lag |tj-ti|. The correlation length is set to ℓc=1 year in this example.
The statistical properties of all inputs are summarized in Table . The stochastic process over t∈[0,20] years is discretized into n=100 time instants for the proposed and PHI2 methods. In the proposed method, the EOLE truncation order is set to M=28 to ensure that the maximum pointwise approximation error over t∈[0,20] years does not exceed 1 %. For the MCS method, the stochastic process is discretized into nMCS=500 time instants, and 108 sample paths are simulated. The result of the MCS method is regarded to be the benchmark, and the results obtained by the proposed method and PHI2 are compared with the MCS benchmark in Fig. .
Failure probability of the steel beam with different methods.
Here, the error of the proposed and PHI2 methods is defined as
ε(ts,te)=pf(ts,te)-pf,MCS(ts,te)pf,MCS(ts,te)×100%,
where pf(ts,te) is the time-dependent failure probability over [ts,te] estimated by the proposed method or PHI2, and pf,MCS(ts,te) is the corresponding estimate obtained by the MCS method.
Failure probability of the steel beam over different periods.
For te=[5,10,15,20] years, the comparison is provided in Table . The result of the proposed method shows good agreement with the MCS benchmark, with a maximum relative error of 12.82 % over the considered time periods. In contrast, the relative error of the PHI2 method increases rapidly with time, reaching a maximum value of 46.9 %. In terms of the function calls, the MCS benchmark corresponds to 5×1010 calls. The proposed and PHI2 methods require 2.45×104 and 4.328×104 calls, respectively. The proposed method is more effective than the PHI2 and MCS methods.
Convergence and sensitivity analysesi. Sensitivity to the time discretization level <inline-formula><mml:math id="M177" display="inline"><mml:mi>n</mml:mi></mml:math></inline-formula>
As the number of time discretization points n increases, the failure event over a continuous time period can be approximated on a finer temporal grid, thereby reducing the time discretization error. On the other hand, the EPA approximates the high-dimensional joint failure probability through repeated compounding of two events. It has been reported that, when the number of component events is large, the approximation error and the sensitivity to the compounding order may accumulate across multiple compounding levels so that the resulting system failure probability tends to be biased on the conservative side . Therefore, the choice of n should balance the time discretization error against the system-compounding error.
In the proposed method, a correlation-based compounding rule is adopted (see Fig. ) to reduce order sensitivity and alleviate error accumulation. To investigate the influence of the discretization level on numerical stability and conservatism, the time period is discretized with n ranging from 30 to 400. For each n, the EOLE truncation order M is determined adaptively following the same criterion as in Sect. . The resulting time-dependent failure probability pf(0,20) as a function of n is shown in Fig. . The MCS benchmark pf,MCS(0,20)=1.328×10-4 is also plotted for comparison, together with a ±10 % band of [1.195×10-4,1.461×10-4]. This band corresponds to an interval of reliability indices [3.6341,3.6655], with Δβ≈0.03, indicating that a ±10 % variation in pf at this probability level translates into only a minor change in β.
Convergence of the time-dependent failure probability with respect to the number of time discretization points n.
When n is small, the coarse temporal grid provides insufficient coverage of potential failure events within the considered time period, leading to an underestimated time-dependent failure probability. As n increases, more potential failure instants are captured, and time-dependent failure probability increases and gradually stabilizes. Once n≥70, the results enter the ±10 % band around the MCS benchmark. With further refinement, the variation of the time-dependent failure probability becomes limited; meanwhile, as the number of discretized events increases, the influence of system-compounding error on the final estimate becomes more pronounced, which may still cause an upward drift in failure probability. Throughout this study, the adaptively selected EOLE truncation order varies only within a narrow range of M≈28–31, suggesting that, under the current error threshold, M is relatively insensitive to the choice of n.
Based on the behavior of time-dependent failure probability in Fig. , a discretization level of approximately n=[90,150] is found to provide a good overall performance for this example. The corresponding time step is Δt≈0.225–0.134 years, and the correlation between two adjacent time instants is ρ(Δt)≈0.95–0.98 according to Eq. (). Accordingly, this correlation range is adopted as the discretization guideline in the subsequent analyses.
ii. Sensitivity to the EOLE truncation order <inline-formula><mml:math id="M204" display="inline"><mml:mi>M</mml:mi></mml:math></inline-formula>
The EOLE truncation order M determines the fidelity with which the stochastic process Y(t) is represented in the finite-dimensional standard normal space. To examine its influence on the time-dependent failure probability, the time interval of [ts,te]=[0,20] years is discretized with a fixed level n=100, and M is varied from 5 to 80. The resulting failure probability over [ts,te] is shown in Fig. . In this figure, the dashed line denotes the MCS benchmark pf,MCS(0,20)=1.328×10-4, and the shaded region indicates a ±10 % band around this benchmark.
Convergence of the time-dependent failure probability with respect to the EOLE truncation order M.
As shown in Fig. , the estimate exhibits a pronounced sensitivity to M at low truncation orders. For M<20, the time-dependent failure probability is notably lower than the benchmark and increases rapidly with M, indicating that the truncated expansion is insufficient to preserve the process fluctuations that govern the occurrence of failure within the considered period. When M increases to approximately 20–40, the estimate approaches the benchmark and falls within the ±10 % band. For M≥50, the curve reaches a clear plateau, and further increasing M produces negligible changes in the time-dependent failure probability, demonstrating diminishing returns from higher truncation orders.
To relate this behavior to the truncation quality, the maximum pointwise error is reported as εmax=maxj=1,…,nε(tj). The corresponding trend of εmax with M is shown in Fig. . The results show an orders-of-magnitude decay of εmax with increasing M. εmax=1.21×10-1 at M=20, which is consistent with the marked underestimation observed in Fig. . At M=28, εmax decreases to 9.65×10-3, and the corresponding time-dependent failure probability estimate is already close to the benchmark. Further increasing M reduces εmax to the order of 10-7 or smaller, whereas the change in failure probability remains negligible.
The maximum pointwise EOLE truncation error with respect to M.
Taken together, Figs. and suggest that controlling the maximum pointwise error at the 1 % level is sufficient to obtain a practically stable estimate of time-dependent failure probability for this example. Tightening the threshold to 10-3 or 10-4 would substantially increase the process dimension while yielding only marginal changes in failure probability. Accordingly, M is chosen as the smallest truncation order satisfying εmax≤1 % in the subsequent analyses.
iii. Failure probability versus uncertainty level
To investigate the influence of the uncertainty level of the stochastic load process on the time-dependent failure probability, an uncertainty scaling factor κ is introduced to proportionally scale the standard deviation of the load process F(t). The mean function and correlation structure of F(t) remain the same as those in Sect. . The considered time period t∈[0,20] years is discretized with n=100 time instants, and the EOLE truncation order is selected adaptively using the same 1 % pointwise error criterion. A parametric scan is conducted over κ∈[0.5,2.0], and the proposed results are compared against an MCS benchmark with 108 sample paths.
As shown in Fig. , the time-dependent failure probability pf(0,20) increases markedly with the uncertainty scaling factor κ. For instance, pf(0,20) rises from 8.074×10-6 at κ=0.5 to 1.347×10-4 at κ=1.0 and reaches 4.915×10-3 at κ=2.0, indicating that amplifying the load fluctuations markedly increases the likelihood of failure within the considered period. The proposed method closely follows the MCS benchmark across the entire range of κ, reproducing both the monotonic trend and the magnitude. The largest relative discrepancy is observed at κ=2.0, where the maximum relative error is about 4.15 %, corresponding to a small reliability index difference of |Δβ|<0.015, which confirms the robustness of the proposed estimate under varying process uncertainty.
Failure probability versus uncertainty level.
The 23-bar truss
An elastic 23-bar truss is depicted in Fig. . The upper and lower chords have Young's modulus E1 and cross-sectional area A1, whereas the web members have Young's modulus E2 and cross-sectional area A2. The overall length of the truss is 24 m, and the height is 2 m. The external loads Pi(t)(i=1,…,6) are assumed to share the same stationary Gaussian process P(t). The autocorrelation function is defined as
ρti,tj=exp-|tj-ti|2.
The model of 23-bar truss structure.
The serviceability requirement is that the maximum displacement does not exceed the allowable value dmax, and the limit state function is expressed as
G(X,P(t))=dmax-d(X,P(t)),
where d(X,P(t)) is the maximum displacement of the truss, and dmax=13 cm is the maximum allowable displacement. The probabilistic characteristics of all of the parameters of the truss are provided in Table .
The time-dependent failure probability of the truss is considered over t∈[0,5] years. The stochastic load process is discretized into n=100 time instants for the proposed method. Because the load process is stationary and the model does not include any time degradation variable, the time-independent reliability problem is identical across all discretized instants. Hence, only one representative FORM analysis is required for the proposed method. In the proposed method, the FORM result at a representative time instant is used for all discretized time instants, and the EOLE truncation order is set to M=52 so that the pointwise approximation error of the load process is controlled within 1 % over the considered period. For PHI2, only a few FORM analyses are needed to estimate the outcrossing rate. For the MCS method, the stochastic process is discretized into nMCS=200 time instants, and 1×105 sample paths are simulated. Figure and Table compare the results of the proposed, MCS, and PHI2 methods. The proposed method agrees well with the MCS benchmark, and the maximum relative error over the considered periods is only 6.24 %. Meanwhile, the PHI2 method is unsuitable for this problem and has a maximum relative error of about 295 %. In terms of limit state function evaluations, the MCS benchmark needs 2×107 calls, PHI2 requires 563 calls, and the proposed method requires only 168 calls.
Failure probabilities of the 23-bar truss with different methods.
Corroded pipeline
The buried natural-gas pipeline is subjected to both environmental corrosion and time-varying internal pressure. The internal pressure of the pipeline should be considered to be a stochastic process because it changes with time. The case study considers an in-service natural-gas pipeline located in eastern China. The pipeline is made of X65 steel with an outer diameter of D=813 mm and a wall thickness of w=14.2 mm, and the operating pressure is 6.3 MPa. The most recent in-line inspection (ILI) reported a dominant corrosion metal loss defect with a maximum depth of d0=7.66 mm and an axial length of l0=57 mm (see Fig. ). In Fig. , the metal loss defect is idealized by its axial extent l0 and the maximum wall thickness loss d0 reported by the ILI. Based on the historical inspection records, an average corrosion growth rate of kp=0.226 mm yr−1 is adopted. The parameters of the pipeline and the defect are provided in Table .
The limit state function of a pipeline is defined as
G(X,t)=Pb(t)-P(t),
where P(t) is the operating pressure, and Pb(t) is the burst pressure.
The operating pressure is modeled as a stationary Gaussian process with the autocorrelation function given in Eq. (), where ℓc=0.2 years for this case. The burst pressure is related to the structure, material, and defect state of the pipeline, which can be obtained by using the semi-empirical formula of ASME B31G ,
Pb(t)=2σy+68.95wD1-d(t)w1-d(t)wM(t),
and Folias factor M(t) is
M(t)=1+0.6275l(t)2Dw-0.003375l(t)4(Dw)2,l(t)2Dw≤500.032l(t)2Dw+3.3,l(t)2Dw>50,
where the σy, w, and D are the yield strength, wall thickness, and outer diameter of the pipeline, respectively, and d(t) and l(t) are the defect depth and length at t, respectively.
Corrosion growth is modeled by a linear law and is predicted under a widely used engineering assumption that the depth-to-length ratio is preserved during growth . Accordingly,
d(t)=d0+kpt,l(t)=l01+kpd0t,
where kp is the average corrosion growth rate inferred from historical inspection records.
Failure probabilities of the corroded pipeline over different periods.
The time-dependent failure probability of the corroded pipeline is evaluated over t∈[0,10] years after the most recent ILI. The evaluation period is discretized into n=270 time instants for the proposed method, and the EOLE truncation order is chosen as M=69 so that the maximum pointwise approximation error over the period does not exceed 1 %. The PHI2 method is evaluated on the same time grid to obtain the time-dependent failure probability. The failure probabilities of the corroded pipeline using the three methods are compared in Fig. . The error comparison at different periods is provided in Table . In the MCS simulation, the pressure process is discretized at nMCS=500 time instants, and 1×108 sample paths are simulated. The proposed method agrees well with the MCS benchmark, and the maximum relative error over the considered period is 10.82 %. Meanwhile, PHI2 yields larger discrepancies and has a maximum relative error of 68.03 %. In terms of limit state function calls, the proposed method requires 3.59×104 calls, whereas PHI2 requires 1.08×106 calls. The MCS benchmark corresponds to 5×1010 calls.
Failure probability of corroded pipeline with different methods.
Conclusions
A new time-dependent reliability analysis method is presented in this work. This method employs FORM, EOLE, and EPA. FORM is used to analyze the reliability at each time instant. EOLE is used to represent the stochastic process by series expansion and to construct its correlated representation in the standard normal space, which supports the subsequent system aggregation over discretized time instants. EPA is employed to estimate the time-dependent failure probability of the series system over the considered period. Three numerical examples demonstrate the efficiency and accuracy of the proposed method. The results of the proposed method achieve close agreement with MCS benchmarks while requiring orders-of-magnitude fewer limit state evaluations. Compared with the PHI2 method, the proposed approach provides consistently improved accuracy in the studied examples.
Nevertheless, we also noted some limitations of the proposed method during this study. EPA inevitably introduces approximation errors when progressively compounding multiple failure events. When the stochastic process has a short correlation length and the response exhibits high-frequency fluctuations, a smaller time step is required to maintain adequate temporal resolution, which increases the number of compounding levels and may aggravate error accumulation and the sensitivity to the compounding order. For non-Gaussian stochastic processes, the dependence among time instants cannot be fully characterized by a correlation coefficient alone, and Gaussian approximations may lead to noticeable bias in tail probability estimates. Moreover, when the limit state function is strongly nonlinear, the first order linearization error of FORM may become non-negligible and can be further amplified through multi-aggregation. Future work will focus on mitigating these issues and improving the robustness of the method.
Code and data availability
The code used in this study contains internal core solver routines and proprietary algorithmic modules developed by our research team, which are confidential under our team's internal confidentiality policy and therefore cannot be publicly released. The raw in-line inspection data and corresponding reports used in the pipeline case study are owned by a third party and are not publicly available due to confidentiality restrictions. Reasonable requests for access to the code or processed data may be considered by the corresponding author, subject to confidentiality agreements and appropriate data use arrangements. All model formulations and parameter values necessary to reproduce the results are fully provided within this article.
The contact author has declared that neither of the authors has any competing interests.
Disclaimer
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. The authors bear the ultimate responsibility for providing appropriate place names. Views expressed in the text are those of the authors and do not necessarily reflect the views of the publisher.
Financial support
This work was supported by the Talent Program of Chengdu Technological University (grant no. 2023RC042).
Review statement
This paper was edited by Zhiwei Zhu and reviewed by two anonymous referees.
ReferencesAndrieu-Renaud, C., Sudret, B., and Lemaire, M.: The PHI2 method: a way to compute time-variant reliability, Reliab. Eng. Syst. Safe., 84, 75–86, 10.1016/j.ress.2003.10.005, 2004. ASME: Manual for Determining the Remaining Strength of Corroded Pipelines, American Society of Mechanical Engineers, New York, NY, USA, standard, ISBN 9780791873847, 2023.Balesdent, M., Morio, J., and Marzat, J.: Kriging-based adaptive importance sampling algorithms for rare event estimation, Struct. Saf., 44, 1–10, 10.1016/j.strusafe.2013.04.001, 2013. Baroth, J., Breysse, D., and Schoefs, F.: Construction reliability: Safety, Variability and Sustainability, Wiley-ISTE, ISBN 9781848212305, 2011.Bredow, R. and Kamiński, M.: Structural safety of the steel hall under dynamic excitation using the relative probabilistic entropy concept, Materials, 15, 3587, 10.3390/ma15103587, 2022.Breitung, K.: Asymptotic approximations for multinormal integrals, J. Eng. Mech., 110, 357–366, 10.1061/(ASCE)0733-9399(1984)110:3(357), 1984.Carr, P.: Riser and pipeline corrosion risk assessment, in: Offshore Technology Conference Asia, OTC–24946, Kuala Lumpur, Malaysia, 10.4043/24946-MS, 2014.Chun, J., Song, J., and Paulino, G. H.: Parameter sensitivity of system reliability using sequential compounding method, Struct. Saf., 55, 26–36, 10.1016/j.strusafe.2015.02.001, 2015.Ditlevsen, O. and Bjerager, P.: Methods of structural systems reliability, Struct. Saf., 3, 195–229, 10.1016/0167-4730(86)90004-4, 1986.Du, X.: Time-dependent mechanism reliability analysis with envelope functions and first-order approximation, J. Mech. Design, 136, 081010, 10.1115/1.4027636, 2014.Gollwitzer, S. and Rackwitz, R.: Equivalent components in first-order system reliability, Reliability Engineering, 5, 99–115, 10.1016/0143-8174(83)90024-0, 1983.Gong, C. and Frangopol, D. M.: An efficient time-dependent reliability method, Struct. Saf., 81, 101864, 10.1016/j.strusafe.2019.05.001, 2019.Gong, C. and Zhou, W.: Improvement of equivalent component approach for reliability analyses of series systems, Struct. Saf., 68, 65–72, 10.1016/j.strusafe.2017.06.001, 2017.Hasofer, A.: Reliability index and failure probability, J. Struct. Mech., 3, 25–27, 10.1080/03601217408907254, 1974.Hawchar, L., El Soueidy, C.-P., and Schoefs, F.: Principal component analysis and polynomial chaos expansion for time-variant reliability problems, Reliab. Eng. Syst. Safe., 167, 406–416, 10.1016/j.ress.2017.06.024, 2017.Hu, Z. and Du, X.: A sampling approach to extreme value distribution for time-dependent reliability analysis, J. Mech. Design, 135, 071003, 10.1115/1.4023925, 2013.Hu, Z. and Du, X.: First order reliability method for time-variant problems using series expansions, Struct. Multidiscip. O., 51, 1–21, 10.1007/s00158-014-1132-9, 2015a.Hu, Z. and Du, X.: Mixed efficient global optimization for time-dependent reliability analysis, J. Mech. Design, 137, 051401, 10.1115/1.4029520, 2015b.Hu, Z. and Mahadevan, S.: Time-dependent system reliability analysis using random field discretization, J. Mech. Design, 137, 101404, 10.1115/1.4031337, 2015.Hu, Z. and Mahadevan, S.: A single-loop kriging surrogate modeling for time-dependent reliability analysis, J. Mech. Design, 138, 061406, 10.1115/1.4033428, 2016.Hu, Z., Mansour, R., Olsson, M., and Du, X.: Second-order reliability methods: a review and comparative study, Struct. Multidiscip. O., 64, 3233–3263, 10.1007/s00158-021-03013-y, 2021.Jones, D. R., Schonlau, M., and Welch, W. J.: Efficient global optimization of expensive black-box functions, J. Global Optim., 13, 455–492, 10.1023/A:1008306431147, 1998.Kamiński, M. and Bredow, R.: On application of the relative entropy concept in reliability assessment of some engineering cable structures, Comput. Struct., 305, 107560, 10.1016/j.compstruc.2024.107560, 2024.Kamiński, M. and Strąkowski, M.: An application of relative entropy in structural safety analysis of elastoplastic beam under fire conditions, Energies, 16, 207, 10.3390/en16010207, 2022.Kang, W.-H. and Song, J.: Evaluation of multivariate normal integrals for general systems by sequential compounding, Struct. Saf., 32, 35–41, 10.1016/j.strusafe.2009.06.001, 2010.Kim, J. and Song, J.: Probability-Adaptive Kriging in n-Ball (PAK-Bn) for reliability analysis, Struct. Saf., 85, 101924, 10.1016/j.strusafe.2020.101924, 2020.Kurtz, N. and Song, J.: Cross-entropy-based adaptive importance sampling using Gaussian mixture, Struct. Saf., 42, 35–44, 10.1016/j.strusafe.2013.01.006, 2013.Li, C.-C. and Der Kiureghian, A.: Optimal discretization of random fields, J. Eng. Mech., 119, 1136–1154, 10.1061/(ASCE)0733-9399(1993)119:6(1136), 1993.Li, J., Chen, J.-B., and Fan, W.-L.: The equivalent extreme-value event and evaluation of the structural system reliability, Struct. Saf., 29, 112–131, 10.1016/j.strusafe.2006.03.002, 2007.Liu, Z., Liu, Z., and Peng, Y.: Dimension reduction of Karhunen-Loeve expansion for simulation of stochastic processes, J. Sound Vib., 408, 168–189, 10.1016/j.jsv.2017.07.016, 2017.Marelli, S. and Sudret, B.: An active-learning algorithm that combines sparse polynomial chaos expansions and bootstrap for structural reliability analysis, Struct. Saf., 75, 67–74, 10.1016/j.strusafe.2018.06.003, 2018. Nikolaidis, E., Ghiocel, D. M., and Singhal, S.: Engineering design reliability handbook, CRC Press, ISBN 9780203483930, 2004.Papadrakakis, M., Papadopoulos, V., and Lagaros, N. D.: Structural reliability analyis of elastic-plastic structures using neural networks and Monte Carlo simulation, Comput. Method. Appl. M., 136, 145–163, 10.1016/0045-7825(96)01011-0, 1996. Pillai, S. U.: Probability, random variables, and stochastic processes, McGraw-Hill, ISBN 9780071122566, 2002.Ping, M., Han, X., Jiang, C., and Xiao, X.: A time-variant extreme-value event evolution method for time-variant reliability analysis, Mech. Syst. Signal Pr., 130, 333–348, 10.1016/j.ymssp.2019.05.009, 2019.Rice, S. O.: Mathematical analysis of random noise, Bell Syst. Tech. J., 23, 282–332, 10.1002/j.1538-7305.1944.tb00874.x, 1944.Rice, S. O.: Mathematical analysis of random noise, Bell Syst. Tech. J., 24, 46–156, 10.1002/j.1538-7305.1945.tb00453.x, 1945.Roscoe, K., Diermanse, F., and Vrouwenvelder, T.: System reliability with correlated components: Accuracy of the Equivalent Planes method, Struct. Saf., 57, 53–64, 10.1016/j.strusafe.2015.07.006, 2015.Sahraoui, Y. and Chateauneuf, A.: The effects of spatial variability of the aggressiveness of soil on system reliability of corroding underground pipelines, Int. J. Pres. Ves. Pip., 146, 188–197, 10.1016/j.ijpvp.2016.09.004, 2016.Son, Y. K. and Savage, G. J.: Set theoretic formulation of performance reliability of multiple response time-variant systems due to degradations in system components, Qual. Reliab. Eng. Int., 23, 171–188, 10.1002/qre.783, 2007.Sudret, B.: Analytical derivation of the outcrossing rate in time-variant reliability problems, Struct. Infrastruct. E., 4, 353–362, 10.1080/15732470701270058, 2008.Sudret, B. and Der Kiureghian, A.: Stochastic Finite Element Methods and Reliability: A State-of-the-Art Report, Tech. Rep. UCB/SEMM-2000/08, University of California, Berkeley, CA, USA, https://escholarship.org/uc/item/3d606904 (last access: 25 March 2026), 2000.Zhang, J. and Du, X.: Time-dependent reliability analysis for function generation mechanisms with random joint clearances, Mech. Mach. Theor., 92, 184–199, 10.1016/j.mechmachtheory.2015.04.020, 2015.Zhang, J. and Ellingwood, B.: Orthogonal series expansions of random fields in reliability analysis, J. Eng. Mech., 120, 2660–2677, 10.1061/(ASCE)0733-9399(1994)120:12(2660), 1994.Zhou, W., Gong, C., and Hong, H.: New perspective on application of first-order reliability method for estimating system reliability, J. Eng. Mech., 143, 04017074, 10.1061/(ASCE)EM.1943-7889.0001280, 2017.