Multidisciplinary systems with transient behavior under time-varying inputs and coupling variables pose significant computational challenges in reliability analysis. Surrogate models of individual disciplinary analyses could be used to mitigate the computational effort; however, the accuracy of the surrogate models is of concern, since the errors introduced by the surrogate models accumulate at each time-step of the simulation. This paper develops a framework for adaptive surrogate-based multidisciplinary analysis (MDA) of reliability over time (A-SMART). The proposed framework consists of three modules, namely, initialization, uncertainty propagation, and three-level global sensitivity analysis (GSA). The first two modules check the quality of the surrogate models and determine when and where we should refine the surrogate models from the reliability analysis perspective. Approaches are proposed to estimate the potential error of the failure probability estimate and to determine the locations of new training points. The three-level GSA method identifies the individual surrogate model for refinement. The combination of the three modules facilitates adaptive and efficient allocation of computational resources, and enables high accuracy in the reliability analysis result. The proposed framework is illustrated with two numerical examples.
Multidisciplinary analysis (MDA) has been intensively studied during the past two decades [1,2]. Current MDA methods can be roughly classified into three groups, namely, the field elimination method, the monolithic method, and the partitioned method [3,4]. Based on these MDA methods, approaches have been developed for the reliability analysis of multidisciplinary systems by considering the uncertainty sources [5–7]. While most of the current multidisciplinary reliability analysis (MDRA) methods have only considered aleatory uncertainty, a couple of studies have also included model uncertainty [2,8,9]. The above MDRA methods have only considered time-independent reliability analysis, whereas time-dependent behavior and uncertainty sources are of interest in many practical multidisciplinary systems .
In the past decades, reliability analysis methods have been investigated for engineering systems with time-dependent uncertainty sources [11–14]. These methods, however, are developed for single disciplinary systems and cannot be directly applied to multidisciplinary systems, due to the complicated couplings between different disciplinary simulations. Motivated by solving this problem, Zhu et al.  recently developed a time-dependent reliability analysis method for a multidisciplinary system using the first-order reliability method (FORM). Since FORM could have large errors for problems with nonlinear limit-state functions, the method presented in Ref.  can only be applied to multidisciplinary systems where FORM is accurate .
This work aims to develop a general time-dependent MDRA method by exploring the surrogate modeling approach, where expensive physics simulation models are replaced by surrogate models. Monte Carlo simulation (MCS) is then performed on the surrogate models to estimate the failure probability. To guarantee the accuracy of reliability analysis, two types of strategies have been pursued for surrogate model refinement, namely, uncertainty minimization and region of interest (RoI). In the former, the surrogate model is refined by minimizing either the bias or variance in the prediction . In the latter, the surrogate model is refined mainly in the RoI (i.e., the limit state), which is important for reliability analysis . The limit-state surrogate modeling methods have been developed for both time-independent [18–20] and time-dependent reliability analysis [11,21], and have shown excellent performance. However, current limit-state surrogate methods focus only on single disciplinary systems and are inapplicable to multidisciplinary systems.
In this paper, we are interested in pursuing the RoI strategy to improve the efficiency and accuracy of the surrogate model in time-dependent MDRA. Several challenges as detailed in Sec. 2.4 are encountered. The main challenge comes from the fact that the surrogate models built for individual disciplinary systems are not directly connected with the quantities of interest (QoIs) with respect to reliability analysis, which brings challenges in determining the RoI. Some questions that need to be answered include when to refine the surrogate models, which disciplinary surrogate models to refine, and where to refine the individual disciplinary surrogate models. This work solves these challenges through the development of a systematic adaptive surrogate modeling framework. In this framework, approaches are proposed to determine whether the surrogate models need to be refined and at what time instant to refine. Along with that, a three-level global sensitivity analysis (GSA) method is developed to adaptively allocate the computational resources to refine the individual disciplinary surrogate models. The developed framework is generic and applicable to both time-dependent and time-independent MDRA. Two numerical examples are used to demonstrate the effectiveness of the proposed framework.
The remainder of the paper is organized as follows: Section 2 provides background concepts on MDRA under time-dependent uncertainty. Section 3 describes the proposed adaptive surrogate modeling framework. Section 4 illustrates the application of the proposed methodology with two numerical examples. Concluding remarks are provided in Sec. 5.
In this section, we first introduce motivation problems of the proposed research. We then briefly review the differences between time-dependent and time-independent multidisciplinary analyses, and the general procedure of RoI-based surrogate modeling methods for reliability analysis. Based on that, we summarize the challenges in the adaptive surrogate modeling for time-dependent MDRA.
Coupled multidisciplinary analysis is very common in the analysis of many engineering systems. The coupling of individual disciplinary analysis is represented as the exchanging physical quantities between different computational models. The required computational effort is very large when sophisticated simulation models are used in the analysis, even for deterministic MDA . One example is the aero-elastic analysis of an aircraft wing, which involves coupled finite element analysis (FEA) and computational fluid dynamics simulation as shown in Fig. 1 . Since time-dependent MDRA requires repeated MDAs under different realizations of uncertainty sources (i.e., random variables and stochastic processes), it is computationally unaffordable when the original simulation models are used.
Another motivation problem is the reliability analysis of a panel structure on a conceptual hypersonic aircraft vehicle . As the vehicle is subjected to hypersonic flow, an attached oblique shock is created at the forebody leading edge . The resulting aerodynamic pressure causes elastic deformation of the panel, which feeds back to alter the aerodynamic pressure on the panel. The panel is also subjected to aerothermal effects (aerodynamic heating and heat transfer). The aerothermal effects are coupled with the aero-elastic effects . A detailed description of this panel structure is available in Refs.  and . Figure 2 shows the four disciplinary analysis models of the panel (aerodynamics, aerodynamic heating, heat transfer, and structural deformation) .
In the four-disciplinary analysis, finite element models are used in each individual disciplinary analysis to predict the response of coupling variables and response variables. A partitioned approach is employed to solve the four-disciplinary analysis over time since the partitioned method enables us to use relatively independent solution procedures for each discipline . The partition sequence in the analysis is: (1) aerodynamic pressure, (2) aerodynamic heating, (3) thermal, and (4) structural . Since the simulation models in each disciplinary analysis are computationally expensive, directly performing time-dependent MDRA is computationally prohibitive. New strategies are required to reduce the required computational effort in time-dependent MDRA without sacrificing the accuracy of reliability analysis.
The nodal responses of FEA or computational fluid dynamics analyses are usually used as coupling variables and response variables in MDA and MDRA. As a result, the coupling and response variables may have both spatial and temporal variability. In this paper, we only focus on coupled systems with temporal variability alone. The developed method in this paper could be integrated with the method presented in Ref.  to handle spatial and temporal variability simultaneously. In what follows, we first provide a generalized model of multidisciplinary analysis under time-dependent uncertainty. Based on the generalized model, we present the proposed new method.
Generalized Multidisciplinary Analysis Under Time-Dependent Uncertainty.
Figure 3 shows a generalized multidisciplinary system with three individual disciplines . The notations of different variables are also given in the figure. In Secs. 2–4, we use to represent input random variables and to represent input stochastic processes. For the ith discipline, the inputs include random variables (), stochastic processes (), and time (). The subscript s refers to shared inputs across multiple disciplines. There are also coupling variables () between the ith and jth disciplines. The response variables () are functions of the coupling variables and input variables.
In time-independent MDRA, MDA is a deterministic analysis for a given realization of uncertainty sources, and will converge after several iterations. This situation also holds for time-dependent MDRA if the time-dependent inputs only affect the system response variables but not the coupling variables . However, when the time-dependent inputs affect the coupling variables, the coupling variables will change over time; thus we do not have distributions of converged values of these coupling variables and will not be able to decouple the disciplinary analyses. As a result, the computational effort becomes prohibitive in time-dependent MDRA, due to the need to run the fully coupled MDA over the entire duration of interest for each realization of the inputs. In this situation, the partitioned method of MDA needs to be employed. As mentioned previously, a partitioned approach to MDA will enable the use of separate solution procedures for each disciplinary analysis and also handle different time intervals in coupled multidisciplinary analysis, i.e., the same time-step or identical discretization for individual disciplines is not required. In addition, note that the partition sequence may affect the result of MDA; the proposed MDRA method is based on the assumption of a predetermined partition sequence from deterministic MDA.
As shown in Fig. 4, at any time instant , the MDA given in Fig. 3 can be reformulated as Fig. 4. Based on the formulation given in Fig. 4, the MDA becomes a multilevel analysis problem at a given time instant and the analysis is coupled together between time instants.
where is a vector of models of coupling variables , is a vector of models of response variables , is a vector of all coupling variables from the kth discipline, is a vector of coupling variables from the jth discipline to the kth discipline with j > k, and are coupling variables from the jth discipline to the kth discipline with j < k.
Next, we briefly review the procedure of adaptive surrogate-based reliability analysis. After that, we discuss the challenges in adaptive surrogate modeling for time-dependent MDRA.
Adaptive Surrogate Modeling for Reliability Analysis.
In adaptive surrogate-based reliability analysis, the original computer simulation models are substituted with surrogate models. Training points are added adaptively in the RoI (i.e., around the limit state) to guarantee the efficiency and accuracy of reliability analysis. As discussed in Sec. 1, adaptive surrogate modeling approaches were first developed for time-independent reliability analysis and have been extended to time-dependent reliability analysis of single disciplinary systems . The general procedure of adaptive surrogate modeling for time-dependent reliability analysis can be summarized as follows:
Step 1: Build the initial surrogate model and generate random realizations of random variables and stochastic processes.
Step 2: Perform predictions using the surrogate model and samples generated in step 1.
Step 3: Check convergence of the failure probability estimate. If the convergence criterion is satisfied, obtain the failure probability estimate. Otherwise, go to the next step.
Step 4: Identify a new training point in the region of interest, update the surrogate model by adding the new training point into current training points, and then go to step 2.
In step 4, the identification of new training points in the RoI can be achieved based on different learning functions (e.g., the expected feasibility function defined in Ref.  or the U function defined in Ref. ) or performing global sensitivity analysis . Even if the detailed implementation procedures may vary with methods, most current adaptive surrogate modeling approaches can be generalized into the above summarized procedure.
Adaptive Surrogate Modeling for MDRA.
where is the safe state of the kth disciplinary response variables, “” means “for all,” is the time interval of interest, and and are the reliability values of series and parallel system, respectively.
For systems with other configurations, the time-dependent reliability needs to be defined according to the reliability block diagram. The objective of this paper is to develop a framework to efficiently and accurately estimate the time-dependent multidisciplinary system reliability by pursuing the adaptive surrogate modeling strategy. However, the procedure summarized in Sec. 2.3 cannot be directly applied to MDRA due to the difference between multidisciplinary system and single disciplinary system. A major difference is that the QoIs in a single disciplinary system are only affected by the uncertainty sources of that discipline, whereas the QoIs in a multidisciplinary system are affected by not only the uncertainty sources in their own discipline, but also by those in the other disciplines. The uncertainty sources are present in each discipline and the uncertainty accumulates as the analysis proceeds from one model to the other and over time. This is also why current available surrogate modeling method for single disciplinary system reliability analysis  cannot be applied to solve the MDRA problem. To efficiently and accurately estimate the time-dependent failure probability, a question that needs to be answered is how to effectively allocate the computational resources to train the individual disciplinary surrogates. To answer this question, the following challenges need to be addressed:
Whether: After we replace the original simulation models with surrogate models, the first question that needs to be answered is whether the surrogate models are good enough to be used to estimate the failure probability. This question is answered by quantifying the uncertainty in the failure probability estimate and assessing whether the uncertainty is larger than a threshold.
Where: At what realizations of the inputs and uncertainty sources to refine the surrogate models? If the surrogate models of the individual disciplinary analyses are not directly connected with the failure probability estimate, currently available learning functions cannot be used to determine where to refine the surrogate models.
When: As shown in Fig. 4, the multidisciplinary system simulation under time-dependent uncertainty is performed iteratively over time. The uncertainty due to various sources will accumulate over time. At what time instant to refine the surrogate models needs to be addressed?
Which: After we identified when and at what realization of the uncertainty sources to refine the surrogate models, the next challenge is to determine which disciplinary surrogates to refine and which surrogate model within the discipline to refine (since each disciplinary analysis is substituted with multiple surrogate models, one for each output of the disciplinary analysis).
In Sec. 3, we develop a framework called adaptive surrogate-based multidisciplinary analysis of reliability over time (A-SMART) to address the above questions.
In this section, we first provide an overview of the proposed A-SMART framework. Following that, we explain the framework in detail.
Overview of the Proposed Framework.
In the proposed A-SMART framework, for the purpose of generalization, we assume that all the simulation models are expensive. We substitute individual disciplinary simulation models with surrogates. In other words, we build surrogate models for all and and determine how to allocate computational resources to refine these surrogate models. Figure 5 shows the flowchart of the main steps of the proposed framework.
The proposed framework (as given in Fig. 5) consists of three main modules, namely, initialization, uncertainty propagation, and three-level GSA. The initialization and uncertainty propagation modules inherit the basic idea of the single loop Kriging surrogate modeling approach presented by Hu and Mahadevan . In these two modules, initial surrogate models are built for the response and coupling variables. Then the surrogate model uncertainty sources are propagated through the time-dependent MDA to quantify the effects of surrogate model uncertainty on the time-dependent failure probability estimate. Based on the results of uncertainty propagation, we can quantify the uncertainty in the failure probability estimate from the current set of surrogate models and thus determine whether the surrogate models need to be refined, where to refine, and when to refine. In the three-level GSA module, we overcome the challenge that the surrogate model uncertainty sources are not directly connected to the time-dependent failure probability by performing GSA at three levels. The first two modules address the challenges of whether, where, and when summarized in Sec. 2.4. The last module solves challenge of which. A brief summary of the three modules corresponding to the procedure given in Sec. 2.3 is given as follows:
Module 1—Initialization: In this module, random realizations of and , are generated. Initial surrogate models are built for and . This is similar to step 1 of the procedure summarized in Sec. 2.3.
Module 2—Uncertainty propagation: This module includes steps 2–4 shown in Fig. 5. The basic idea is similar to steps 2–3 discussed in Sec. 2.3. Due to the multilevel simulation at each time instant and the coupling between different time instants, the implementation details are quite different from that of single disciplinary analysis even if the basic idea is similar.
Module 3—Three-level GSA: This module covers steps 5–8 of the flowchart presented in Fig. 5. In this module, we first identify the response variable that makes the highest contribution to the uncertainty in the failure probability estimate (this is level 1 GSA) since the response variables are directly connected to the failure probability. After the response variable is identified, we identify which surrogate model (i.e., response surrogate model or coupling variable surrogate model) to refine to reduce the uncertainty in the response variable through level 2 and 3 GSAs.
Module 1: Initial Surrogate Modeling.
According to the description of multidisciplinary system given in Sec. 2.2, at a given time instant , a two disciplinary system is simulated as shown in Fig. 6(a). For each discipline, there are models of coupling variables and response variables. For example, as depicted in Fig. 6(b), for the kth discipline, there are models of and , .
In module 1, we build initial surrogate models for and , . In this paper, the surrogate models are built using the Kriging surrogate method. For the sake of illustration, in Sec. 3.2, we use to represent all the inputs of .
When we are building the surrogate models, we do not know the domains of the coupling variables. To generate training points, we need to get a rough guess of the domains by performing several deterministic MDAs over the time duration of interest with fixed random realizations of random variables and stochastic processes. Once we have a rough guess of the domains, we generate training points of , , and t. We then evaluate the models of individual disciplines separately at these training points. Based on the training points, we build the initial surrogate models for and , . Note that the surrogate models are built for one time-step MDA as indicated in Fig. 4. Here, the generated training points are therefore for partitioned MDA at one time-step. In addition to the generated training points, the MDA iterations used to get the initial guess of coupling variables ranges can also be added into the pool of initial training points for surrogate modeling. Even if the initial training points may not fully cover the ranges of , , and t, more training points will be added in important regions later in the proposed method through adaptive sampling approach to improve the accuracy of reliability analysis.
where and represent the surrogate model approximations of and , respectively.
These surrogate models are then used to substitute the original simulation models in time-dependent MDA. Along with the surrogate models, in the initialization module, we also generate random samples for , , where is the number of time instants in the time interval . We denote the generated random samples as , , where and represent the qth sample of and , and are the qth trajectory of and at time instant , and is the number of MCS samples.
where is the number of response variables from the kth discipline, and and are the mean and variances of the prediction of at time instant . The mean and variance are obtained from Kriging prediction.
The above equations imply that after substituting the original simulation models with surrogate models, we introduce surrogate model uncertainty into MDA. Next, we will discuss the propagation of the surrogate model uncertainty and investigate how to determine the quality of the surrogate models based on the uncertainty propagation.
Module 2: Uncertainty Propagation in Surrogate Model-Based MDA
Uncertainty Analysis of Surrogate Model-Based MDA.
In this work, to investigate the effect of surrogate model uncertainty on the failure probability estimate, we first quantify the uncertainty in the response variables due to the uncertainty of the coupling and response surrogates for any given realization of the input variables.
Taking the multidisciplinary system shown in Fig. 6(a) as an example, for the qth sample of the input variables (i.e., , ) generated in Sec. 3.2, the MDA under uncertainty is a multimodel uncertainty quantification (UQ) problem as shown in Fig. 7. In the network, square nodes represent deterministic nodes and elliptical nodes represent random nodes. At any time instant , a multimodel UQ is performed and the uncertainty quantified at will be propagated to due to the coupling between time instants. This process continues until the maximum number of iterations is reached.
where and are the conditional mean and variance predicted from the Kriging surrogate model by using as inputs.
where is the wth sample of the pth response variable of the kth discipline at .
Error Analysis of the Multidisciplinary System Failure Probability Prediction.
where is the probability of making an error on the sign of the pth response of the kth discipline at time instant i for given , and stands for probability.
where is the safety state indicator of the pth response of the kth discipline for the qth realization of random variables and stochastic processes.
In addition, if there exists such a that and , we can directly assign . This is due to the fact that once the response is classified as failed at a time instant and the probability of making an error on the classification is low, the corresponding trajectory can be classified as failed.
The above equation indicates that if the kth disciplinary response needs to be refined, it should be refined at the input setting . From Eqs. (20)–(25), the questions of whether, where, and when (summarized in Sec. 2.4) have been answered. Next, we investigate which surrogate to refine through the development of a three-level GSA scheme.
Module 3: Three-Level GSA.
As shown in Fig. 5, there are three levels of analysis: coupling variables→response variables→system safety indicator. To minimize the maximum probability of making an error by reducing the surrogate model uncertainty, a three-level GSA scheme is developed.
Level 1 GSA: Which Discipline to Refine?.
where is the ith input variable, is the vector of variables excluding variable , is the variance of quantity of interest (B), and is the expectation by freezing .
where is computed using Eq. (19), is computed using Eq. (16), and is computed by combining the RBD with Eq. (16).
Level 2 GSA: Response Surrogate or Coupling Surrogate?.
where is the wth sample of generated in Sec. 3.3.1, is the joint probability density function of computed using kernel smoothing density function, and is the conditional variance.
From Eq. (38), we have the new training point of as . Based on the new training point, surrogate model is updated. After that, the level 2 GSA is performed again until we have .
If we get in the first round level 2 GSA, it implies that the uncertainty in mainly comes from the coupling surrogate models. In that case, we will go to the level 3 GSA to decide, which coupling surrogate to refine.
Level 3 GSA: Which Coupling Surrogate?.
where are the auxiliary variables corresponding to the coupling variables from the jth discipline to the th discipline at , is a new function defined through a network defined similar to Fig. 8, is the auxiliary variable of the th response of the th discipline at , and includes all other deterministic variables. For the sake of illustration, we do not list all of the involved deterministic variables, which include the realizations of random variables, stochastic processes, and parameters of the Kriging models.
where , represent the other auxiliary variables, is computed by performing uncertainty propagation using the network given in Fig. 8, and is computed in Eq. (37).
After the coupling surrogate model is identified, the new training point for the surrogate model is identified using the weighted mean square error criterion given in Eq. (38) and the samples of the coupling variables at time instant . After the coupling surrogate models are updated, the distributions of the coupling variables are updated. The updated samples of the coupling variables are then used to perform level 2 GSA again. If , we go to the next step. Otherwise, level 3 GSA is performed again to further refine the coupling surrogate.
Until now, we have discussed all the implementation procedure of the proposed A-SMART framework. Next, we will use numerical examples to demonstrate the effectiveness of the proposed method.
In this section, a mathematical example and an engineering example are used to illustrate the proposed method. In each example, the proposed method is compared with two other methods, namely, variance-minimization method and MCS. For MCS, 5 × 105 realizations are generated for random variables and stochastic processes and the Karhunen–Loève (KL) expansion method is used to generate samples for stochastic processes . In the variance-minimization method, the surrogate models are refined by minimizing the maximum variance over iterations. In the Kriging surrogate, the zeroth order trend function is used.
A Mathematical Example.
where and . The time interval is discretized into 20 time instants in MDA.
where is the results estimated from MCS.
It shows that the proposed method is able to estimate the time-dependent system failure probability efficiently and accurately. For the 651 NOF of the proposed method, the NOFs are allocated as , , , , , , , , and , where means the NOF allocated to . It implies that more computational resources are allocated to models , , and than the other models. Figure 10 shows the comparison of the convergence history of from the proposed method and the variance minimization-based method with respect to NOF. It indicates that the proposed method converges much faster than the variance minimization-based method.
A Compound Cylinder.
where and , representing the effect of corrosion on the thickness of cylinders over time.
where is the correlation length of pressure.
where and .
We perform time-dependent system reliability analysis for the compound cylinder. In the initial surrogate modeling, 13 training points are generated for each model. Table 2 gives the results comparison between the proposed method and MCS. Table 3 lists the allocation of NOFs to different simulation models. Following that, Fig. 13 shows the comparison of convergence history of different methods. Similar conclusions are obtained as that from the first example.
The above two numerical examples demonstrate that the proposed method can effectively allocate the computational resources to efficiently and accurately perform time-dependent MDRA. In the proposed method, only temporal variability is considered in MDRA.
As discussed in our motivation problems (i.e., Sec. 2.1), the coupling variables and response variables may have both spatial and temporal variabilities. When the spatial variability is also considered, the coupling and response variables will be high-dimensional, which brings significant challenge to the surrogate modeling. In that situation, dimension reduction methods need to be employed to reduce the dimension and thus make the surrogate modeling possible. For instance, the singular value decomposition or proper orthogonal decomposition methods can be adopted first to map the high-dimensional coupling and response variables into low-dimensional variables in latent space. Surrogate models are then built in the low-dimensional latent space. Such a method has been presented in Ref. .
Based on the dimension reduction methods, the proposed adaptive surrogate modeling framework can be employed to reduce the required computational effort for surrogate modeling to achieve an accurate prediction of the system failure probability. The integration of the proposed adaptive surrogate modeling method with the dimension reduction method, however, is not straightforward. The convergence criterion and the three-level GSA methods need to be integrated with dimension reduction methods, and need to be further extended to account for the spatial variability. Integration of the proposed adaptive surrogate modeling method with the method presented in Ref.  will be studied in our future work.
This paper developed a new A-SMART framework for time-dependent multidisciplinary system reliability analysis. Surrogate models are built separately for each coupling variable and response of each disciplinary analysis. Computational resources are then adaptively allocated to the surrogate models to improve the efficiency and accuracy of reliability analysis, through the integration of three modules: initialization, uncertainty propagation, and three-level GSA. A mathematical example and an engineering application example demonstrated that the proposed method can effectively allocate the computational resources to efficiently and accurately perform time-dependent MDRA.
In the proposed framework, sampling-based approaches are employed to perform uncertainty propagation and global sensitivity analysis. For time-dependent MDRA over a very long time period, the sampling-based approaches might be computationally expensive even if it is much cheaper than using the original simulation models. In addition, the correlation between the predictions of individual samples as discussed in Ref.  is not considered in the developed framework. Accounting for the correlation effect will further improve the accuracy and efficiency of the proposed method.
Air Force Office of Scientific Research (Grant No. FA9550-15-1-0018).