Abstract

This article presents an approach that combines the active global Kriging method and multidisciplinary strategy to investigate the problem of evidence-based multidisciplinary design optimization. The global Kriging model is constructed by introducing a so-called learning function and using actively selected samples in the entire optimization space. With the Kriging model, the plausibility, Pl, of failure is obtained with evidence theory. The multidisciplinary feasible and collaborative optimization strategies of multidisciplinary design optimization are combined with the evidence-based reliability analysis. Numerical examples are provided to illustrate the efficiency and accuracy of the proposed method. The numerical results show that the proposed algorithm is effective and valuable, which is valuable in engineering application.

1. Introduction

The performance of structures in engineering system is often influenced by uncertainties. The optimal solution obtained by conventional deterministic design frequently falls on the boundary of the feasible region [1]. The deterministic design may fall into the failure region if the system is disturbed by uncertainties. In the multidisciplinary coupling systems, the propagation of uncertainty factors between different subjects makes the search for an optimal design more complicated as well as difficult than the search for an optimal individual discipline design due to the coupling effect [2, 3]. The influence of uncertainty factors must be considered to ensure the reliability of the optimized design results, and reliability-based multidisciplinary design optimization (RBMDO) needs to be performed [4, 5].

Uncertainties can be categorized as aleatory and epistemic [6, 7]. Aleatory or objective uncertainties arise from the inherent randomness of a system. Probability theory can be adopted to handle random uncertainties. However, sufficient information is required to construct the distribution function of uncertain variables. Epistemic or subjective uncertainties stem from the lack of sufficient information during modeling and optimization. Several theories, including probability box models, fuzzy sets, Bayesian approaches, and evidence theory, have been developed to deal with epistemic uncertainties [815]. Evidence theory provides a general modeling of epistemic uncertainty and can be reduced to other theories. When the interval of evidence theory is infinite, it is equivalent to traditional probability theory. When no conflict exists between uncertain pieces of information, the theory is equivalent to possibility theory. When the subinterval of the evidence variable is unique, the theory degenerates into convex model theory.

Several methods have been proposed for reliability analyses based on evidence theory, and these approaches include Cartesian product method (CPM) [16], Monte Carlo simulation (MCS) [1720], polynomial chaos expansion [21], and surrogate-based models [2224]. Design optimization based on evidence theory has also been studied [2528]. Zhang et al. [29] adopted two-stage framework to handle the evidence-based design optimization problem, where the evidence variables were transformed into random variables and the sequential optimization and reliability assessment were used. He and Qu [30] gave an overview of possibility and evidence theory-based design optimization. However, limited research has been conducted on evidence-based multidisciplinary design optimization (EBMDO) [3133]. Conventional EBMDO has a three-level loop architecture. The outer loop is multidisciplinary optimization, the middle loop is reliability analysis based on evidence theory, and the inner loop is discipline analysis. The computational cost of EBMDO is high, especially in the inner discipline analysis. Practical engineering problems require repeatedly calling the simulation model to evaluate limit state and objective functions, and this process is time consuming.

To address this issue, this study proposes a surrogate model-based method with an active global learning strategy, which constructs surrogate models in entire design space with intelligent sampling approach, to reduce the amount of computation [34]. The active learning Kriging [35] model is utilized as a surrogate model to evaluate the limit state function in evidence-based reliability analysis. Since the original active Kriging model is utilized for reliability analysis, a model in local region, here the Kriging model is constructed in the entire design space called active global Kriging model. Then, two multidisciplinary integration frameworks, namely, Multidisciplinary Feasible (MDF) and Collaborative Optimization (CO), are adopted and combined with the Kriging model. Using the global Kriging model, the computational cost of EBMDO is reduced to the approximate equivalent in deterministic multidisciplinary design optimization (MDO).

The rest of the article is organized as follows. MDF and CO strategies for EBMDO are briefly introduced in Section 2. An evidence-based reliability analysis using the global kriging model for MDO is presented in detail in Section 3. Numerical examples are provided and discussed in Section 4, and the conclusions are presented in Section 5.

2. Methodology

2.1. EBMDO

The general EBMDO model is expressed aswhere d pertains to the deterministic variables, ds pertains to shared design variables that are common among all disciplines, di pertains to the design variables of the ith discipline, X pertains to the uncertain variables, denotes the mean of the X, which is treated as the uncertain design variables, Xs pertains to uncertain shared evidence variables that are the common input variables for all disciplines, and Xi pertains to the local evidence variables of the ith discipline. The shared variable Xs and the input variable X are independent variables, K encompasses coupling variables, and f is the objective function. is the upper bound of failure probability which represents the worst-case and is called plausibility. is the failure probability constraint of the ith discipline, and is the limit state of the ith discipline. The failure model is defined by . is the reliability index, and denotes the cumulative distribution function (CDF) of the standard normal distribution. The coupling variables between design and uncertain variables are expressed aswhere the first subscript j and the second subscript i of denote the values from the jth discipline to the ith discipline. The meaning of can be explained in detail that the coupling variables are inner shared variables in the jth discipline and the ith discipline, and is the output values of the jth discipline and the input values for the ith discipline. For example, in fluid structure thermal coupling analysis, the pressure values from fluid analysis should be transferred to the structure discipline for finite element analysis. The conventional structure of EBMDO contains a triple loop, as shown in Figure 1 [36, 37].

As shown in Figure 1, due to the nested architecture, failure plausibility must be calculated at every design point during optimization. Therefore, the multidisciplinary analysis is repeatedly called to evaluate the limit state function. Furthermore, according to evidence theory, the optimization problem should be implemented to determine the upper bound of the limit state. This nested structure makes EBMDO costly.

2.2. Multidisciplinary Feasible Method

The framework of MDO aims to manage coupled variables and subsystems and improve calculation efficiency. The MDO framework has two main forms, namely, single-level and multilevel structures.

MDF is a traditional single-level optimization method. A systematic analysis is performed for each iteration during optimization, and a multidisciplinary analysis is performed repeatedly to obtain consistent solutions of the coupling variables among disciplines. This optimization framework includes a system of integrated and subsystem analysis modules. Only one interface exists for the input and output between the system module and subsystem modules, as shown in Figure 2. The system module is an optimizer to search for the optimal solution. The subsystem modules are for multidisciplinary and reliability analyses. In each optimization iteration, the coupling variables should be consistent and the constraints should be also satisfied. Therefore, when the optimal solution is obtained, the consistency constraints of coupling variables and reliability constraints are satisfied.

2.3. Collaborative Optimization Method

CO is a bilevel framework of MDO. The framework of MDO is designed to be top level (system level) and low level (discipline level) by CO. The top level is a system optimizer responsible for assigning target values for system-level state variables to disciplines. The low level is a parallel, distributed, multiple-discipline subsystem. In addition to satisfying the constraints of the subsystems, the objective function is the smallest difference between the coupling state variables of the top-level and low-level systems. After low-level system optimization, the objective function is fed back to the top-level system, which constitutes the top-level consistency constraint. Optimization of the top-level system is performed to address the uncoordinated state variables among the subsystems. The framework for CO is defined as follows:where f is the objective function of system-level problems, ds_sl denotes the shared deterministic design variables of system-level problems, dsl denotes the deterministic design variables of discipline-level problems, Xs_sl denotes the shared uncertain design variables of system-level problems, Xsl denotes the uncertain design variables of the discipline-level problems, Psl denotes the global uncertain parameters, and is the constraint of the discipline-level problems.

The ith discipline subproblem iswhere Ji is the objective of the discipline-level problems, dsl denotes the shared design variables, di denotes the design variables of the ith discipline, si denotes the shared uncertain design variables, and Xi denotes the uncertain design variables of the ith discipline.

The diagram of CO for EBMDO is shown in Figure 3.

3. Reliability Analysis

3.1. Reliability Analysis Based on Evidence Theory

The limit state function is denoted as G(x). Failure domain F is defined as , where the evidence variables are . The evidence variable can be described by evidence space , where is the set of all focal elements [21] and , where is the ith focal element. In the evidence analysis, the focal element is usually an interval. is the basic probability assignment (BPA) of each focal element. The joint focal element (JFE) is defined by CPM:

The BPA associated with JFE is calculated bywhere is the number of JFEs. In the evidence space, failure probability is an interval. The lower and upper bounds of the interval are called belief (Bel) and plausibility (Pl), respectively. Bel means that the system is in the absolute failure state, and Pl means that the system is in the possible failure state. The belief and plausibility measures of failure probability are calculated bywhere is the indicator function.

The equivalent form of equations (7) and (8) are

This method involves a combinatorial explosion problem because it needs to explore the extrema of in JFEs. According to equations (7) and (8), if variables and focal elements exist for each variable, then the total number of JFEs will be . The Monte Carlo method is adopted to randomly sample the focal elements and overcome the obstacle of dimensions. JFEs are formed by these focal elements. An extremum analysis is performed based on JFEs, and Bel and Pl can be estimated using statistical methods, shown aswhere N is the number of samples for the Monte Carlo simulation (MCS).

3.2. Evidence Analysis with the Kriging Model

Although the dimensionality challenge is addressed, the computational cost remains large for MCS. To further improve the computational efficiency, the Kriging model is proposed to evaluate the limit state function in a multidisciplinary coupling system. The Kriging model is briefly introduced as follows.

Suppose that a set of N samples called design of experiments (DoE) has the form (x, G), where x is an n-dimensional vector and G is the corresponding output . The predicted value and variance are given aswhere 1 is an n-dimensional unit vector, r is the vector of the correlation function of x and points in DoE, R is the matrix of the correlation function of points in DoE, and are the parameters to be determined. The methods of the fitting parameters are adopted as modified DIRECT algorithm [38]. Further details about the Kriging model can be found in Ref. [39].

For evidence-based design optimization, the Kriging model is used to approximate the limit state function in the entire design space. The idea here is that if the sign given by the Kriging model is the same to that of the true limit state function, then the result of MCS combined with the Kriging model can be considered trustworthy. The following property can ensure that the sign of the extremum value by the Kriging model is consistent with that of the true limit state function.

Property. In the entire design space , if , then in the local design space , and where () is the sign function.

Proof. Consider the case of the minimum. Suppose is the minimum point of , and is the minimum point of , namely, and . We will discuss in two cases.

Case 1. If because , then because , so .

Case 2. If because , then because , so .
Therefore, when the limit state function has a minimum value, the sign by the Kriging model and the true limit state function are the same. The case of the maximum value can be proven in a similar way. Suppose that the Kriging model has already been constructed to correctly estimate the sign of the limit state function, equations (11) and (12) can be written asTo reduce the sample size, we use the expected risk function (ERF) [23, 39] as learning function to adaptively select the samples.where and are the probability density function and CDF of standard normal distribution, respectively. The convergence criterion for the Kriging model is written aswhere is the value for maximum of ERF, is the mean of the X, and is .

3.3. Validation of Reliability Analysis with the Evidence Variables

The following numerical example is adopted to validate the proposed reliability analysis method. The limit state function is written as [13]where is a deterministic design variable and and are evidence variables. When , the system is in the failure region. The BPA structures of and are listed in Table 1.

There are 13 points selected to construct the Kriging model. Then, Monte Carlo simulation is performed to calculate the Pl of the failure probability. The results of the reliability analysis are shown in Table 2. As can be seen, the result is close to the result from literature [13], which means the proposed evidence analysis method has a high accuracy.

3.4. Summary of EBMDO with Kriging-Based Reliability Analysis

First, a large population S that contains points is generated. Points are randomly selected from population S. The limit state function is evaluated based on these points, thus forming the initial DoE, SDoE, for the Kriging model. The Kriging model is then constructed, and the values of ERF in population S are calculated. The point with the largest value of ERF is selected and added to SDoE. The Kriging model is updated by SDoE until convergence is achieved.

Second, a reliability constraint is applied during optimization by the Kriging model. The framework of MDF and CO is adopted to achieve an optimal design.

The procedure of EBMDO is summarized as follows:Step 1. Generate points S according to the range of the variablesStep 2. Randomly select points from S, evaluate the limit state function, and form the initial SDoEStep 3. Construct the Kriging model based on SDoEStep 4. Calculate the ERF on each point in S, and pick out point with the maximum value of ERFStep 5. Evaluate the limit state function on the point , and add the sample in SDoEStep 6. If the Kriging model is converged, go to next step; otherwise, go to Step (3)Step 7. Set the initial value of the optimizationStep 8. At each optimal point during optimization, evaluate the reliability of the system by the Kriging modelStep 9. If the reliability constraint is satisfied and the objective function is converged, stop; otherwise, go to Step (8)

During Kriging modeling, MDA should be performed to evaluate the limit state function. During optimization, MDA is not required for reliability analysis because the Kriging model is constructed to approximate the limit state function in the entire design space. Therefore, computational efficiency is considerably improved.

4. Application Examples

The proposed methodology is validated with two examples of EBMDO. For the inequality constraints of failure probability, evidence theory is used to estimate the epistemic uncertainty in terms of failure plausibility.

4.1. Case 1

The first EMBDO problem is a simple mathematical example. The problem is defined aswhere are the design variables and are the evidence variables, and . The system utilizes two design variables and outputs two states. They are not shared variables, thus the systems are uncoupled. The objective function is nonlinear with two constraints. is the constraint of failure probability. The BPA structure of the uncertain variables is shown in Table 3, and the corresponding bar chart of epistemic variables is illustrated in Figure 4.

The CO architecture of this problem can be written as one system optimization and two discipline optimizations. The system optimization is written aswhere are the system-level design variables, denotes the ith variables in the jth discipline, and are two auxiliary variables and the objective function of the discipline optimization problem. The auxiliary variables aim to reduce the difference of the system and discipline variables to . The optimization problem of discipline 1 is defined as

The optimization problem of discipline 2 is defined as

The initial value of the design variables are set to [1.6, 0.8], which is the result of deterministic optimization (DO). To verify the accuracy of the proposed approach, MCS combined with the MDF architecture is used to directly solve EBMDO. The nonlinear sequential quadratic programming (NLPQL) is adopted as the optimization algorithm for both system and discipline optimization.

Firstly, the candidate points are generated by Latin hypercube sampling. Then, according to experience, 12 points are randomly selected and evaluated by executing the disciplinary analysis. The Kriging model is constructed and updated by the proposed active global learning method. Subsequently, the MDF and CO framework are performed based on the Kriging model to search the optimal solutions.

The results are summarized and compared with those of the deterministic design and MCS in Table 4. The second row shows the deterministic results. EBMDO is performed with three failure probability constraints, namely, . The results show that the objective function of MDO is smaller than that of EBMDO. The constraints of failure probability are not satisfied. For a small failure probability constraint, the design variables and the objective function must be large to meet the requirement of failure plausibility. The obtained results are close to those of MCS, and only the second constraint is active. The outer optimization loop has 100 iterations, and for each design point, times of MCS are performed in the evidence analysis. For each evidence analysis, the limit state function is called for about 10 times to obtain the minimum value. Hence, the limit state function is called for approximately times. The number of limit state function calls by the proposed method is 26 times for this example, which is only of MCS.

4.2. Case 2

The second example contains two disciplines with three design and two coupled variables. The constraint function is nonlinear. The EBMDO problem is defined bywhere are design variables, are coupled variables, are uncertainty variables, and is the constraint of failure probability, . In this example, the coupled variables make the problem highly complex. Probability theory is a special case of evidence theory. To compare the EBMDO design with the RBMDO design, the BPA for each interval of the uncertain variables is considered to be similar to the area under the probability density function used in RBMDO, as shown in Figure 5.

In RBMDO, pertains to normally distributed random parameters with . The BPA structure of the uncertain variables for EBMDO is shown in Table 5.

The samples are generated in the entire design space, namely, . Subsequently, 12 points are randomly selected to construct the initial DoE. Then, the learning strategy mentioned in Section 3 is utilized to select the points for DoE. In each selected point, the multidisciplinary analysis must be performed to obtain the response of the limit state function. The training history curve of the Kriging model is shown in Figure 6.

The Kriging model converges in 9 and 13 iterations for the first and second limit state function, respectively. To verify the accuracy of the Kriging model, the signs of samples that are randomly generated in the design space are compared with that of the true limit state function. The histogram of the sign is illustrated in Figure 7. The predicted signs have a good match with that of the true function. The number of signs that are incorrectly predicted is only four, while the relative error is .

Then, the problem is solved using the CO strategy. The EBMDO problem can be decomposed into one system problem and two corresponding discipline problems because of the presence of two disciplines. The system optimization problem is defined aswhere contains the system design variables and denotes the design variables of the ith discipline. The coupled variables, , are considered the design variables to decouple the disciplines.

The optimization problem of the first discipline is defined by

The optimization problem of the second discipline is defined by

The data flowchart is presented in Figure 8 to show the solving architecture of cooperative optimization. The blue and red lines indicate the input and output, respectively. The system variables are passed to the discipline analysis, and the objective functions of disciplines are returned to system optimization. The objective functions of disciplines are the constraints of system optimization. The reliability constraints are evaluated by evidence analysis with the global Kriging models.

Table 6 compares the DO, EBMDO by MDF, EBMDO by CO, RBMDO, and MCS results. The conclusions are similar to those in the previous example. The evaluation number of the limit state function by the proposed approach is 46, which is about that by MCS. This result indicates that the proposed method can ensure accuracy while having a great advantage over MCS.

4.3. Case 3

The third case is multidisciplinary design optimization for cooling blade which considers the aerodynamic, heat transfer, and strength analysis, revised from [40]. The geometry model of the cooling turbine blade and the design parameters are shown in Figure 9, where is the thickness of ribs and is the thickness of blade wall at different blade profile.

The EBMDO problem is defined bywhere x is the uncertainty design variables, f is the objective function, and D is the damage of cooling turbine blade, which is the function of the design variables. In this case, the damage D is the whole damage, which is constituted by two parts, namely, the creep damage and the fatigue damage. The creep damage is predicted by the Larson–Miller equation [41] and the fatigue damage is calculated by the nominal stress method [42]. The fluid-thermal-solid interaction analysis is carried out to obtain the stress and the temperature in MDF framework. Also, in the CO framework, the three disciplines are decoupled and performed independently. is the average temperature, and are the weight factor, is the maximum temperature of cooling turbine blade, is the maximum stress, P is the reliability, and x is the uncertainty design variables including . The design variables are considered as epistemic parameters. The lower and upper limit of variables are shown in Table 7. The thickness difference obeys three-parameter Weibull distribution, with the location parameter, the shape parameter, and the scale parameter as −0.1540, 2.2979, and 0.2923, respectively. According to the Weibull distribution and Monte Carlo sampling, the BPA of the design variables are shown in Table 8.

The EBMDO results compared with the RBMDO and MDO by [40] are listed in Table 9. The multidisciplinary feasible approach is adopted in this case. As can be seen, the result obtained by EBMDO is very close to RBMDO, which indicates the proposed method is effective.

5. Conclusions

This article presents a novel approach for EBMDO that combines the active global Kriging model, MDF, and CO. By introducing a learning function, the global Kriging model is adaptively constructed to replace the limit state function. Then, the evidence analysis is performed by MCS and the Kriging model. The computational cost of EBMDO is reduced to approximately that of DO. Three examples are provided to illustrate the proposed approach. Compared with MCS, the proposed method can obtain accurate results and hold a significant advantage in terms of computational efficiency. Therefore, the proposed method is expected to be of great value in engineering applications.

As part of further work, some other architectures of multidisciplinary optimization can be taken into account in EBMDO, such as concurrent subspace optimization (CSSO) and bi-level integrated system synthesis (BLISS).

Data Availability

All relevant data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work has been supported by the Natural Science Foundation of Jiangsu Province (BK20190424), the Natural Science Foundation of Shaanxi Province (2019JQ-470, 2019JQ-609), the National Key R&D Program of China (grant no. 2017YFB0202702), and the Center for High Performance Computing and System Simulation of Pilot National Laboratory for Marine Science and Technology (Qingdao).