An Efficient Reliability-Based Optimization Method Utilizing High-Dimensional Model Representation and Weight-Point Estimation Method

Xiaoyi Wang,Xinyue Chang,Wenxuan Wang,⋆,Zijie Qiao and Feng Zhang,⋆

1School of Mechatronic Engineering,Lanzhou Jiaotong University,Lanzhou,730070,China

2School of Electronic and Information Engineering,Lanzhou Jiaotong University,Lanzhou,730070,China

3School of Mechanics,Civil Engineering and Architecture,Northwestern Polytechnical University,Xi’an,China

ABSTRACT The objective of reliability-based design optimization (RBDO) is to minimize the optimization objective while satisfying the corresponding reliability requirements.However,the nested loop characteristic reduces the efficiency of RBDO algorithm, which hinders their application to high-dimensional engineering problems.To address these issues, this paper proposes an efficient decoupled RBDO method combining high dimensional model representation(HDMR)and the weight-point estimation method(WPEM).First,we decouple the RBDO model using HDMR and WPEM.Second,Lagrange interpolation is used to approximate a univariate function.Finally,based on the results of the first two steps, the original nested loop reliability optimization model is completely transformed into a deterministic design optimization model that can be solved by a series of mature constrained optimization methods without any additional calculations.Two numerical examples of a planar 10-bar structure and an aviation hydraulic piping system with 28 design variables are analyzed to illustrate the performance and practicability of the proposed method.

KEYWORDS Reliability-based design optimization; high-dimensional model decomposition; point estimation method;Lagrange interpolation;aviation hydraulic piping system

1 Introduction

Many uncertainty factors inevitably exist in practical engineering problems and affect the performance of structural systems to some extent[1–3].Therefore,it is critical to consider the uncertain influence of structural systems on reliability-based design optimization(RBDO)[4,5].RBDO produce more reliable design results compared to those obtained using traditional deterministic optimization methods [6,7].Existing RBDO methods can be roughly divided into three categories: the nested double-loop method(NDLM),single loop method(SLM),and decoupled-loop method(DLM).

NDLM is the original RBDO method and is executed with alternating inner and outer layer calculations.The outer layer performs deterministic optimization analysis for the design variables and the inner layer performs probabilistic analysis for the reliability constraints [8].Because the innerlayer probability analysis requires a large number of computations and each outer-layer optimization iteration must invoke inner-layer probability analysis, the efficiency of RBDO based on NDLM is relatively low.Regardless, this method can be considered as a milestone in RBDO and has received extensive attention in the field of structural system design[9].

Efficient algorithms such as SLM and DLM have emerged to improve the computational efficiency of RBDO.For SLM,the processes of searching for optimal design variables and calculating the most probable failure points occur simultaneously.This process overcomes the drawbacks of NDLM to some extent and has been widely applied in practical engineering problems.Wang et al.[10]proposed a non-probabilistic reliability-based optimization method based on SLM and applied it to the optimization of a supersonic wing.Their method uses first-order interval Taylor expansion, the interval vertex theorem,and direct optimization methods to analyze the uncertainty propagation problem of a system.Additionally, the extended non-probabilistic set theory-stress strength interference model and volume ratio theory are introduced to evaluate the non-probabilistic index reasonably.Jiang et al.[11] proposed an adaptive hybrid single-loop method (AH-SLM).An iterative control strategy(ICS)with two iterative control criteria is proposed.This strategy converges rapidly,regardless of the degree of non-linearity of the performance function, and improves the accuracy of the most probable point (MPP) in highly nonlinear problems.Keshtegar et al.[12] proposed a novel RBDO based on SLM and applied it to the design of aircraft stiffened plates with complex bucking constraints.The results demonstrated that their method has enhanced computational efficiency and robustness.Additional studies on SLM include references[13–15].One can see that SLM is not only suitable for probabilistic RBDO,but also provides a good performance for non-probabilistic RBDO.Currently,the research on this method mainly focuses on handling highly nonlinear probabilistic constraints and rapid convergence.To reduce the computational burden further,the RBDO of complex structure systems through the combination of surrogate models and SLM is an emerging trend[16].

DLM first transforms reliability constraints into general deterministic constraints.Converted RBDO problems can then be solved utilizing traditional deterministic optimization methods.In[17,18],it is replaced the MPP with a quantile-based decoupling method to overcome the drawbacks of the original method that requires solving the MPP.This is particularly important for highly nonlinear problems, multiple-MPPs and non-normally distributed variables.In this approach, quantitative metrics are obtained by sampling a surrogate model and a new sample is selected using a sampleupdating strategy in the significant region [19,20].Overall, this method yields good performance for complex structures.Shi et al.[21] utilized a two-step method to implement a novel form of decoupling.Because this method involves only a small amount of time-varying reliability analysis and deterministic optimization,the efficiency of time-varying reliability optimization is improved to some extent.Zhang et al.[22] used a Kriging surrogate model [23,24] to handle RBDO problems with independent constraint functions and then proposed a novel quantile-based sequential RBDO approach.In this method, an adaptive Kriging scheme with error control is integrated to derive the accuracy information of the surrogate model.The proposed independent training strategy reduces the number of function evaluations while maintaining the accuracy of reliability estimation.Other studies on DLM include references[25,26].Most DLM implementations are similar to SLM.However,there is no method for completely decoupling the reliability of the inner layer and optimization analysis of the outer layer because it is difficult to determine the functional relationships between reliability constraints and design variables directly.

In this study, a fully decoupled RBDO method is developed by combining high dimensional model representation(HDMR)[27]and the weight point estimation method(WPEM)[28].First,the reliability constraint function is decomposed using HDMR and then combined with WPEM and the first-order reliability method(FORM)[29]to decompose the reliability index into a combination of a series of univariate functions.Second,the univariate function is fitted using Lagrange interpolation[30].Finally, a traditional deterministic optimization method is utilized to obtain optimal design variables.

The remainder of this paper is organized as follows.Section 2 briefly introduces the basic principles of RBDO and HDMR.Section 3 presents the derivation of the reliability index based on HDMR and WPEM.The proposed decoupled RBDO model is detailed in Section 4.In Section 5,two numerical examples of a planar 10-bar structure and an aviation hydraulic piping system(AHPS)with 28 design variables involving finite element analysis are presented to illustrate the performance and practicality of the proposed method.Finally,our conclusions are summarized in Section 6.

2 Brief Review of RBDO and HDMR

2.1 Reliability-Based Design Optimization

The goal of RBDO is to derive optimal design variables that meet the reliability requirements of an engineering system.Generally,RBDO is mathematically expressed as follows:minimizef(d)

wherex=(x1,x2,...,xn)Tandd=(d1,d2,...,dm)Tare the random parameter and design variable,respectively,andnandmare the numbers of random parameters and design variables,respectively.fis the objective function,andRk,gk,εkandare the reliability,performance function,security threshold and target reliability levels of thekth reliability constraint,respectively.Ncis the number of reliability constraints,andare the lower and upper bounds of theith design variabledi,respectively,andNdis the number of design variables.

2.2 High Dimensional Model Representation

Two types of HDMR are commonly used in the reliability field, namely the multiplication dimension reduction method(MDRM)[31]and conventional dimension reduction method(CDRM)[32].In this study,MDRM is utilized to handle reliability constraints and it is formulated as follows:

wherexc=represents a specific reference vector of a random vector.For a random variable that follows a normal distribution,its reference point is generally regarded as its mean value.denotes the reference points other thanxiandgiis the univariate function ofxi.

3 Reliability Index Based on MDRM and WPEM

Because the relationship between the reliability level and reliability index isR= Φ(β), the reliability constraint in Eq.(1)can be transformed into

where Φ is the cumulative distribution function, andβkandβtakare the reliability indices of theith eventgk(x,d)≤εkandkth target reliability index,respectively.Then,Eq.(1)can be rewritten as

Many methods can be used to calculateβ,including FORM[27]and the second-order reliability method [33].Because FORM is most commonly used in actual engineering structural reliability analysis,we used FORM to calculateβ,which is expressed as

whereμg,σgandεdenote the mean and standard deviations of the performance function and security threshold,respectively.

According to MDRM,as shown in Eq.(2),a performance function containing both random and design variables can be decomposed into the following form in a similar manner:

wheregjdenotes the univariate function of thejth design variabledjandrepresents the reference vector of the design vector,which can be considered as the most valuable point among the design variables.In this study,the central point of the design variables is considered as the reference point.denotes the reference vector, excludingdj.The performance function is decomposed into two parts using Eq.(6),where one is the product of the random univariate functionsgiand the other is the product of the design univariate functiongj

Next,the relationship between the reliability indexβand design variablesdis deduced based on Eq.(6).The deduction process for the meanμgand varianceis as follows:

whereξandζare expressed as follows:

whereλ1xi,λ2xi,λ3xiandλ4xidenote the first four moments ofxi.The first four moments are the mean,standard deviation,kurtosis and skewness in order.The corresponding standard normally distributed variables are zero,one,zero and three in order.

By combining Eq.(5)with Eq.(7)to Eq.(10),the reliability indexβcan be rewritten as

As shown in Eq.(19),βis only a function of the design variables and is independent of the random variables.

4 Decoupled RBDO Model

According to the basic principles described in the previous section, two types of decoupled RBDO models can be formulated.In the first RBDO model, the reliability level is given byR=Pr(g(x,d)≤ε)and the corresponding reliability constraint can be deduced as follows:

whereC=(ζβta+ξ)is a constant.Furthermore, the RBDO model can be decoupled into the following deterministic optimization model:

In the second RBDO model,the reliability level is expressed asR= Pr(g(x,d)≥ε).Therefore,the reliability index can be rewritten as follows:

Then,the reliability constraint can be deduced in the following form:

whereC′=(ξ-ζβta)is a constant.Furthermore, the corresponding decoupled RBDO model is expressed as follows:

The optimization models defined in Eqs.(21)and(24)are deterministic optimization models with respect to the design variables that can be solved using any traditional optimization algorithm.In this study,Monte Carlo simulation(MCS),sequential quadratic programming(SQP)[35]and particle swarm optimization(PSO)[36]we adopted to verify the performance of the proposed method.

Eqs.(21)or(24)can be used directly to solve the RBDO model when the constraint is an explicit function,whereas for an RBDO model with an implicit constraint function,Lagrange interpolation is adopted to fit the univariate functionas follows:

In Eqs.(25)and(26),Nldenotes the number of samples,xjdenotes thejth input sample,andykdenotes thekth output sample.In this study,the input samples are evenly distributed within the design interval.A flowchart of the proposed RBDO algorithm is presented In Fig.1.

As shown in Fig.1,the calculation steps of the proposed method can be divided into two main parts.The first is calculating two constantsξandζfor the random variables and the second is fitting the univariate functionusing Lagrange interpolation.Therefore, the computational cost of the proposed approach depends on these two operations.Assuming that the numbers of random variables,design variables,and constraint functions aren,mandk,respectively,and that there arepinput samples in each design variable,the number of function evaluations in the proposed method is 3 ×n+ (p×m) ×k.For example, whenn,m,kandpare 2, 2, 2 and 4, respectively, the number of function evaluations for RBDO analysis is 22,which demonstrates that the proposed approach is relatively efficient.

Figure 1:Flowchart of the proposed RBDO algorithm

Notably,the application of the proposed method must focus on the following two points.

(1) The premise of applying the proposed method is assumption that the performance function of all constraints is greater than 0.Otherwise,the performance functions in the reliability constraint must be transformed to ensure that they are greater than 0.For example,suppose that the performance function inR=Pr(g(x,d)≤ε)is less than 0,i.e.,g(x,d)≤0.Then,there is a non-negative constantγthat makes Eq.(27)true.

In Eq.(27),the value ofγshould be sufficiently large to ensure thatg′(x,d)=g(x,d)+γ >0,g′(x,d)is the corresponding new performance function,which must be greater than 0.Assume thatε′=ε+γis the corresponding new threshold.Then,R′=Pr(g′(x,d)≤ε′)is the corresponding new reliability level andR′=R.Finally,the reliability constraint is transformed as follows:

(2)According to Wang et al.[10],there are three combination types between the design variabledand random variablex.The first type is direct mixing of the design variable and random variable(i.e.,(d,x)).The second type is that the design variable is a parameter of the random variable (e.g., such asx∼N(d,σ)).The third type contains both the first type and second type.However,any of these types can be transformed into the first type by the following formula:

whereU∼N(0,1).Therefore,only the first combination type is considered in this study.

5 Illustrative Examples

In this section,two numerical examples of a planar 10-bar structure and AHPS that requires finite element software calculations are presented to illustrate the performance of the proposed method.The proposed method is set up to use SQP and PSO(abbreviated as PM-SQP and PM-PSO,respectively)to perform the optimization analysis,and the results of direct MCS(MCS-D)and a combination of MCS and PSO (MCS-PSO) are used as references.It is noteworthy that MCS-D is a direct MCS method that has the greatest reference value.However, because the optimization module is still MCS, there may be some errors.Therefore,using PSO can further ensure the correctness of the reference results(MCS-PSO).PM-SQP and PM-PSO are implemented based on the proposed method with Lagrange interpolation proxy models after complete decoupling.The difference between these models lies in the optimization methods,which are only used to verify the accuracy of the models.

5.1 Numerical Example 1

A numeric al example with two design variables(d1,d2)and two random variablesx2∼N(0,1)is defined in Eq.(30).The target reliability ofRtais set to 0.9987 withd1∈[2,5],d2∈[2,5],andε=1.

The samples required to approximate the input-output relationship of the univariate function using Eqs.(25) and (26) are defined as=(2,3,4,5).Given 105samples for MCS,the optimization results are presented in Table 1.

Table 1: Results calculated by different RBDO methods

The reliability levelRin Table 1 is calculated by applying MCS 10 times and the design variable values are the corresponding median values.The reliability level obtained by MCS-D is 99.86%,which is lower than the target reliability level of 99.87%.This is because the optimal design variables obtained by MCS-D are the optimal solutions based on 105samples,so the median value of 10 calculations is generally less than the target reliability level and it is reasonable to use the optimal objective function valuef(d)obtained by MCS-D as a reference.For MCS-PRO,the inner layer uses MCS for reliability analysis,and the outer layer uses PRO for optimization analysis.The results obtained by MCS-PRO can be used as an auxiliary reference,because the error of MCS-D is relatively large when the number of design variables is large or the sample size is limited.

In Table 1,one can see that thef(d)values calculated by PM-SQP and PM-PRO are the same at 3.40,which is slightly larger than the reference value of 3.22 calculated by MCS-PRO.The error of the proposed method is approximately 5.92%.The reliability level corresponding to the optimal design variables obtained by PM-SQP and PM-PRO is the same and equal to 99.99%,which is slightly higher than the target reliability level of 99.87%.

The relationship between the number of samples and the minimum value off(d)is presented in Fig.2.One can see that the value off(d) remains largely unchanged after 3 samples and 4 samples can meet the accuracy requirement.The number of function callsNcallis 14 for the proposed method,demonstrating that the computational efficiency of the proposed method is relatively high.

Figure 2:Objective function value changes with the number of samples

5.2 Numerical Example 2

A numerical example with four constraints that has been widely used for benchmarking purposes in previous studies [37,38] is defined in Eq.(31).Two random variables are normally distributed asandx2∼andd1andd2are two design variables with the same design range of[2,5].The target reliabilityRtais set to 0.9987 and four thresholds are assigned asε1=ε2=ε3=ε4= 1.Given an the input sample=(2,3,4,5), the optimization results are presented in Table 2.

Table 2: Results calculated by different RBDO methods

The results in Table 2 are obtained by executing each method once.The MCS results from [38]are presented to verify the correctness of the reference method.The results obtained by MCS and MCS-D are similar, indicating that it is reasonable to use the results of MCS-D as a reference.The minimum values obtained by PM-SQP and PM-PRO are very similar and the error is approximately 5.34%relative to the results obtained by MCS-D.

As shown in Fig.3,because the objective function value is unchanged after 3 samples,4 samples can meet the fitting accuracy requirement.The number of function evaluations for the proposed method is 38,indicating that the efficiency of the proposed method is relatively high.

Figure 3:Objective function value changes with the number of samples

The two numerical examples demonstrate that the accuracy and efficiency of the proposed method are within acceptable ranges.Its accuracy depends on the accuracy of the formula used to calculate the reliability index.In this study, the first-order reliability method, which is the most widely used in practical engineering,resulting in a conservative calculation error of approximately 5%.Additionally,its computational efficiency depends entirely on the number of feature points in WPEM and the number of sample points fitting the univariate function, which are unrelated to subsequent optimization.Next,we apply the proposed method to two practical engineering problems using finite element software to demonstrate the feasibility of the proposed method for practical engineering problems.

5.3 Planar 10-Bar Structure

The plane 10-bar structure is presented in Fig.4.The length of the rods isL,and the corresponding cross-sectional area isAi(i=1,2,...,10), and the elasticity modulus and three external loads areEandPi(i=1,2,3), respectively.LandAi(i=1,2,...,10)are the 11 design variables, andEandPi(i=1,2,3)are 4 random variables that follow normal distribution.An optimization model aimed at minimizing the volume of the 10-bar structure is established and expressed by Eq.(32).

wherex= [x1,x2,x3,x4] = [E,P1,P2,P3], and the meanμand standard deviationσare listed in Table 3.Additionally,d=[d1,d2,...,d11]=[L,A1,A2,...,A10],g(x,d)is the maximum displacement response,andrepresents the lower and upper bounds ofdi,respectively,and the corresponding values are listed in Table 3.

Table 3: Specific parameter values of design and random variables in a 10-bar structure

Figure 4:Planar 10-bar structure

Given the thresholdε= 0.0035(m)and target reliability levelRta= 0.999, the optimization process and results are presented in Fig.5 and Table 4,respectively.

Table 4: Optimization results of a planar 10-bar structure

Because PRO is a heuristic global optimization algorithm that is more suitable for the optimization of complex problems,only utilized PM-PRO to optimize the planar 10-bar structure.The finite element analysis model of the planar 10-bar structure is presented in Fig.6,where the design variables are fixed atd1=1 anddi(i=2,...,11)=0.001,and the maximum displacement is 0.0027(m).It is noteworthy that this example required the use of the ANSYS commercial software and it would be a laborious process to calculate the reliability level using MCS.Therefore,the reliability levels listed in Table 4 are obtained by substituting the optimal design variables into Eq.(21).

Figure 5:Optimization process of a planar 10-bar structure

Figure 6:Displacement contour plot before optimization

For this engineering optimization problem with 11 design variables,Fig.7 reveals that the objective function value remained essentially unchanged after 6 sample points.Therefore,7 samples are taken for each design variable to satisfy the basic accuracy requirements and only 89 function calls are required to solve the optimization problem.For an optimization model with 11 design variables and 4 random variables,the proposed method is relatively efficient.The FE model corresponding to the optimal solution of PM-SQP proposed in Table 4 is presented in Fig.8.

Figure 7:Objective function value changes with the number of samples

5.4 Reliability Optimization of the Clamp Support Position in an AHPS

An AHPS is required for most aircraft operations and is one of the most critical aircraft accessories.Considering the long span of an AHPS, a large number of clamp supports are typically required to fix it and reduce the instability caused by engine vibration and fuselage flutter.A series of AHPS clamp support positions are presented in Fig.9, while Fig.10 presents the corresponding right-side view.

Figure 9:AHPS clamp support positions

Figure 10:Right-side view of an AHPS

In Fig.9,Ais the pump source, andBandCare the two oil outlets.The numbers in Fig.9 corresponding to nodes in the finite-element model,which are also the clamp support positions.The corresponding coordinates are listed in Table 5.The bold coordinates represent the initial values of the design variables and the superscripts are the corresponding serial numbers of the design variables(total of 28 design variables).As shown in Table 5,the design variables are the coordinates of the clamp position nodes and most of the design variables areXcoordinates(e.g.,nodes 14,17,and 51).Some of the design variables are theXandYcoordinates of the corresponding notes(e.g.,nodes 45,62 and 69)and there are three design variables for nodes 82 representingX,Y,andZcoordinates.

Table 5: Coordinates of the clamp support locations(mm)

It is noteworthy that the coordinates in Table 5 are global coordinates,and the pipeline is modeled using the PIPE 16 unite in ANSYS.Considering the characteristics of the PIPE 16 unit,the pipeline is constructed using local coordinate,where the last node of each segment corresponded to the start of the next segment.Fig.11 illustrates this type of construction for node 14 in Fig.10.

Figure 11:Illustration of PIPE16 unit modeling

In Fig.11, one can see that the global coordinate in theXdirection of node 14 is -509, which corresponds to the value for node 14 listed in Table 5.The establishment of this section of the pipeline between nodes 13 and 14 begins at node 13 and the local coordinate inXdirection is -43.For convenience, only local coordinates are used to describe the initial valuesd0, lower boundsd land upper boundsd uof the design variables,as shown in Table 6.The choice of upper and lower bounds for the design variables fo llowed the princip le that the shape o f the pipeline cannot change when the design variables vary between the upper and lower bounds.Furthermore,the shape of the pipeline does changes whenever if a variable exceeds the upper or lower bounds.Therefore, the range o f variable changes in Table 6 is optimal.pipeline wall thicknesst,f luid densityρf,ambient tem peratureT,oil pressureP,elastic modulusE,

Table 6: Parameter information regarding the design variables

The random variables for the AHPS are the pipeline material densityρp,pipeline outer radiusD,and Poisson’s ratioμ.For sim p licity,we letThe specifi c values of these variables are listed in Table 7 and additional values can be found in[39].

Table 7: D istribution types and parameter values of random variables

The random excitation applied to all clam p supports is an acceleration power spectral density function, as shown in Fig.12 [39].The scaling factor of excitation applied at the three endpoints of the AHPS(i.e.,A,B,C)is 1,and 0.025 is applied at the other clam p supports.Because the excitation exerted on an AHPS is mainly caused by engine vibration and fuselage f luttering and mainly occurs in the direction of gravity,only the excitation in theYdirection is considered in the AHPS model.

Figure 12:Acceleration power spectral density function of random excitation

In the current state of the AHPS(i.e.,when the design variable is at the initial value and the random variable value is at parameter 1),the maximum value of the stress response and displacement response are 2.12×107(Pa)and 3.89×10-4(m),as shown in Figs.13 and 14,respectively.

Figure 13:Stress response before optimization

Figure 14: Displacement response before optimization

Both stress and the displacement responses affect the safety of an AHPS.To improve further improve the reliability of the AHPS further,the following RBDO model is established:

In Eq.(33),Smaxdenotes the maximum stress response,d∗represent the optimal design variable,xcis the mean value of the random variables(i.e.,parameter 1 shown in Table 7)andDmaxis the maximum displacement response.It is noteworthy that unlike the previous three examples in which the objective function is only a function of design variables,Smax(d,x)contains both design variables and random variables.Therefore, Eq.(33) cannot be directly solved by the proposed method.Because random variables only affect the reliability constraints,it is reasonable to replaceSmax(d,x)withSmaxas the objective function.Given the thresholdεand target reliabilityRtaof 4 × 10-4(m) and 0.999,respectively,based on PRO,the optimization process and results are presented in Fig.15 and Table 8,respectively.

Table 8: Optimal values of 28 design variables(×10-3m)

Figure 15:Optimization process of the AHPS

For practical problems, only PM-PRO is used to perform optimization analysis because it has better global stability than SQP.When comparing Figs.13 and 16, one can see that the result is reduced from 2.12×107(Pa)to 1.48×107(Pa)through the RBDO design,representing a reduction of approximately 31.19%.Simultaneously,as a reliability constraint,the displacement response is reduced from 3.89×10-4(m)to 2.73×10-4(m),representing a reduction of approximately 27.51%which can be found in the comparison between Figs.14 and 17.When substituting the initial value in Table 6 and the optimal values in Table 8 into Eq.(33),the resulting reliability levels are 0.9980 and 0.9991,respectively, indicating that the optimized reliability level satisfies the requirements.Because there are only one reliability constraint in this reliability optimization problem, the number of function evaluations used to approximate the performance function in the reliability constraint is 3 × 8 +28 × 5 = 164.It is noteworthy that the sample size used to fit each univariate function is 5.This value is selected based on the empirical results obtained from several previous examples (i.e., when the sample size is greater than 4,engineering accuracy can be guaranteed).It should be noted that the RBDO of this AHPS is different from previous optimization problems,where the objective functions are explicit functions that does not need require additional computation.The optimization objective in the AHPS is an implicit function and additional calculations are required for each iteration.To reduce the computational cost, the implicit optimization target is approximated utilizing Lagrange interpolation and the concrete process is similar to that of the optimization constraint.The number of function calls required to approximate the implicit objective function is 28×5=140 and the error is approximately 1.22%.Therefore,the total number of function calls is 304 for the RBDO of the AHPS.

Figure 16:Stress response after optimization

Figure 17:Displacement response after optimization

6 Conclusions

In this study, the RBDO model is completely decoupled into a general constrained optimization problem based on MDRM and WPEM, and Lagrange interpolation is used to approximate a univariate function.Therefore, the RBDO problem is transformed into a general constrained optimization problem.Because the transformed model is a proxy model,no additional computation is required for optimization analysis following transformation.As a result,the computational cost of the proposed method mainly depends on decoupling and the approximation of reliability constraints.In the first example with a single reliability constraint,the number of function calls is 14 and the error is 5.92%.The number of function calls is only 38 for the second example with four reliability constraints.The function calls are 89 and 164 for two finite element examples with 15 and 36 variables,respectively.These examples fully demonstrated the excellent practicability,efficiency,and accuracy of the proposed method.

Acknowledgement:The authors disclosed receipt of the following financial support for the research,authorship and/or publication of this article:This work was supported by the Innovation Fund Project of the Gansu Education Department.

Funding Statement:The authors disclosed receipt of the following financial support for the research,authorship and/or publication of this article:This work was supported by the Innovation Fund Project of the Gansu Education Department(Grant No.2021B-099).

Author Contributions:Wang Wenxuan and Zhang Feng are mainly responsible for theoretical exploration and subsequent syntax modification; Wang Xiaoyi is responsible for theoretical derivation,algorithm programming and paper writing;Chang Xinyue and Qiao Zijie are responsible for example verification and analysis.

Availability of Data and Materials:Algorithm code and related data are available via email:wxwangd@163.com.

Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.