Probabilistic Fatigue/Creep Optimization of Turbine Bladed Disk with Fuzzy Multi-Extremum Response Surface Method.

To effectively perform the probabilistic fatigue/creep coupling optimization of a turbine bladed disk, this paper develops the fuzzy multi-extremum response surface method (FMERSM) for the comprehensive probabilistic optimization of multi-failure/multi-component structures, which absorbs the ideas of the extremum response surface method, hierarchical strategy, and fuzzy theory. We studied the approaches of FMERSM modeling and fatigue/creep damage evaluation of turbine bladed disks, and gave the procedure for the fuzzy probabilistic fatigue/creep optimization of a multi-component structure with FMERSM. The probabilistic fatigue/creep coupling optimization of turbine bladed disks was implemented by regarding the rotor speed, temperature, and density as optimization parameters; the creep stress, creep strain, fatigue damage, and creep damage as optimization objectives; and the reliability and GH4133B fatigue/creep damages as constraint functions. The results show that gas temperature T and rotor speed ω are the key parameters that should be controlled in bladed disk optimization, and respectively reduce by 85 K and 113 rad/s after optimization, which is promising to extend bladed disk life and decrease failure damages. The simulation results show that this method has a higher modeling accuracy and computational efficiency than the Monte Carlo method (MCM). The efforts of this study provide a new useful method for overall probabilistic multi-failure optimization and enrich mechanical reliability theory.

1. Introduction

Mechanical structures are usually assembled by a several components; for example, the rotor system of an aero engine is assembled by a spindle, disk, blade, and other components [1]. If we directly establish the reliability optimization design model of an overall structure involving multi-material, multi-disciplinary, and multi-physics structures, the computational burden will become very large in analysis, so that computational efficiency is unacceptable [2]. Therefore, it is significant to propose an efficient method for an overall reliability optimization design of multi-component and multi-failure modes, to make computational precision and efficiency satisfy engineering requirements.

Recently, numerous methods on structural reliability optimization design have emerged [3,4,5]. The response surface method (RSM) is widely used in reliability optimal design for high efficiency and precision. Zhang et al. [6] firstly proposed an extremum response surface method to complete the reliability optimization of a two-link flexible manipulator; Fei et al. [7,8,9] studied an importance degree model with the extremum response surface method for the dynamic reliability optimization design of a mechanical assembly relationship such as turbine blade-tip radial clearance. However, the traditional RSM can’t meet the reliability optimization design of complex mechanical structures in modeling accuracy and simulation efficiency. To solve this issue, advanced response surface methods were developed recently. Song et al. [10] established a multiple response surface model by using back propagation-artificial neural network to construct a limit state function and establish a multi-objective reliability-based optimization model with a dynamic multi-objective particle swarm optimization algorithm for a reliability optimization design of an aero-engine blisk under fluid–thermal–structure coupling. Hamzaoui et al. [11] proposed an integrated method for estimating the resonance stress of blades with super high strength by combining the inverse of artificial neural network inverse (ANNI) with the Nelder–Mead optimization method. Rodríguez et al. [12] applied a probabilistic design procedure to a group of 10 blades of a low pressure (LP) stage steam turbine of 110 MW, in order to compute the stress changes and reliability due to variations in: damping, natural frequencies, vibration magnitude, and density. The computed vibration stresses were analyzed by applying probability distributions and statistical parameters of input and output to compute the useful life. Wang et al. [13] introduced evidence variables and fuzzy variables to describe cognitive uncertainty parameters and presented a novel dual-stage reliability analysis framework where the first stage incorporates the evidence information by the belief and plausibility measures and the second stage incorporates the fuzzy information by a membership function-like formula. Gao et al. [14] proposed an accurate and efficient fatigue prognosis based on a distributed collaborative response surface method, a substructure-based distributed collaborative probabilistic analysis method (SDCPAM), and a substructure analysis method. Ai et al. [15] discussed a probabilistic framework for fatigue reliability analysis. These works implement reliability-based optimization for many analytical objectives, through analyzing the submodels and then processing the response of submodels to carry out the overall design and analysis. The basic thought in the above works for handling multi-objective design problems provides an enlightened insight to reveal the overall reliability-based optimization design of turbine bladed disks with many failure modes, such as stress failure, strain failure, fatigue damage, creep damage, and so forth. However, since scientific research has its own development laws, the reliability optimization design was carried out in one failure mode at that time, without considering the correlation between the failure modes, and the fuzziness of the constraint boundary conditions.

Most works on aero-engine turbine blades regard the randomness of variables (parameters). Alongside the randomness, actually, some parameters in blade models, such as density, temperature, elastic modulus, boundary conditions, and so forth, possess obvious fuzziness centering on a certain value [16]. In fact, the probabilistic fatigue/creep optimization design of turbine bladed disks involves an obvious fuzziness for design parameters and constraint conditions as well as the coupling among many failure modes such as stress failure, strain failure, creep damage, fatigue damage, and so on [15,16,17,18]. Meanwhile, the fuzziness and coupling seriously negatively influence the design precision and efficiency of multi-object optimization when the above methods are directly applied. Therefore, it is urgent to propose an effective method for multi-object reliability-based optimization, in which the fuzziness for design parameters and constraint conditions as well as the coupling among many failure modes are fully considered in order to improve the modeling accuracy and simulation efficiency.

The objective of this paper is to attempt to propose a fuzzy multi-extremum response surface method (FMERSM) regarding failure correlation and parameter fuzziness, to improve the accuracy and efficiency of the overall dynamic reliability optimization design for a multi-component structure with multi-failure mode, by reasonably handling the transients. Then, the probabilistic fatigue/creep optimization design of an aero-engine bladed disk was effectively implemented with respect to this method, and the developed FMERSM is validated by a comparison of methods.

The remainder of this paper is organized as follows. The fuzzy multi-extremum response surface method (FMERSM) is studied in Section 2, comprising the FMERSM modeling approach, fatigue/creep theory, and the basic thought of the comprehensive probabilistic optimization of a bladed disk with FMERSM. Section 3 implements the fuzzy reliability-based optimization of bladed disk fatigue/creep damage including a parameters selection, finite element (FE) modeling, surrogate modeling, probabilistic fatigue/creep analysis, and method validation. In Section 4, some main conclusions are summarized.

2. Methods and Models

The extremum response surface method (ERSM) was firstly developed to simplify the modeling complexity for the transient probabilistic design of mechanical structures by considering the extreme values of the response process in sample extraction [6]. ERSM has been validated to have high-computational efficiency and acceptable accuracy relative to RSM, in the probabilistic design and optimization of aerospace structures/components [7,8,19,20,21,22]. The multi-extremum response surface method was proposed to handle the multi-model problem in the transient probabilistic analysis of multi-component structures, multi-discipline, and multi-failure modes by assimilating ERSM [19,23,24,25]. In most of the structural probabilistic designs, in fact, influential parameters and constraint conditions hold obvious fuzziness and seriously influence design precision. Therefore, it is reasonable to consider the fuzziness of design parameters and constraint conditions to improve the probabilistic design of structures, especially with multi-failure modes or multi-component structures. In respect of the heuristic thought of MERSM, this paper develops FMERSM with the consideration of fuzzy parameters and constraints to implement the fuzzy reliability-based optimization of bladed disk fatigue/creep damage.

2.1. FMERSM Modeling

Assuming that a structure system includes m components and one component has n failure modes (m,n∈Z) (the sample number of failure models is assumed in this study), as well as indicating the input random variables of the jth failure mode in the ith component (for instance, the creep failure of a blade in a bladed disk system) and y((t,() is the corresponding output response, enough of a data set {ymax(t, (): j∈} consisting of the maximum output responses of y((t,() in the time domain is employed to fit the extremum output response y [24]:

When the quadratic polynomials are considered, Equation (1) is rewritten as:

Regarding the fuzziness and randomness of data in Equation (2), the model comprising numerous sub-models (), the FMERSM model, for multi-failure structure, can be structured as: in which is the fuzzy random input variable vector of the jth failure mode in the ith component, and is the corresponding extremum output response. and are the constant term, linear term, and quadratic term of the jth failure mode in the ith component, respectively. and are denoted by: where (m, n = 1, 2, …, k) are elements (or components) in , and respectively.

The modeling process of Equation (3) regards the randomness and fuzziness of design parameters and constraints based on FMERSM. Therefore, this model (Equation (3)) is called a FMERSM model in this paper.

2.2. Fatigue/Creep Modeling for Probabilistic Optimization of Bladed Disks

Under fatigue/creep coupling failure mode, this paper adopts FMERSM to complete the fuzzy probabilistic fatigue/creep optimization of bladed disks. For a structure system with m components, indicates the fuzzy optimization parameters of the ith component. The main plan is to minimize the objective function subject to the overall reliability performance R(x, w, Dc, D) and coupling critical damage D, which is a single-objective constrained optimization problem. The sub-plan is to maximize the reliability R () of the ith component subject to mechanical load and constraints, which is a multi-objective constrained optimization problem. By introducing pseudo-variables [26], the cyclic optimization between the main plan and sub-plans is done until the convergence condition is satisfied. The fuzzy probabilistic optimization model is shown in Equation (7). where is the ith design variable; and w is the random parameters of mechanical load and material property. represent the lower and upper limit of the ith fuzzy design variables; D is the total amount of creep damages; D is the total amount of fatigue damages; D is fatigue–creep coupled critical damage; denotes the stress and deformation of a component; and is the allowable range of . By the λ level-cut method, the fuzzy subset is decomposed into the common set G(λ*), as explained in Equation (8); then, the problem of fuzzy probabilistic constrained optimization can be transformed into the conventional probabilistic optimization design problem [27]. where is allowable constraint of the jth component stress and deformation; and λ* is optimal horizontal cut set.

2.3. Miner Linear Accumulation Damage Law

Under the interaction between fatigue and creep, the overall damage of the structure is equal to the sum of fatigue damage and creep damage, which is the Miner linear cumulative damage law [28] as follows: in which n is the number of stresses acting on a component; n is the number of cycles acted by the jth stress; N is the fatigue life under the jth acting stress; n is the number of stress levels; t is the hold time of the ith stress; and T is the creep failure time of the ith stress.

When the structure is destroyed (D = 1), the relationship between D and D [29] is: where θ1 and θ2 are fatigue–creep characteristic parameters.

The strain fatigue life prediction model is used to predict the low-cycle fatigue life. in which is the amplitude of total strain; is the fatigue life; is the fatigue strength coefficient; is the fatigue ductility coefficient; is the fatigue strength index; and is the fatigue ductility index.

The creep life prediction equations commonly used in material manuals include creep life prediction equations and thermal strength parameter synthesis equations. The persistence equation is expressed in the form of the thermal intensity parameter synthesis equation. in which is durable strength; is the undetermined coefficient in which r indicates the subscript of the rth coefficient in Equation (12); is the thermal intensity parameter; t is the hold time of the ith stress; and c and C are the constants related to fatigue ductility and temperature, respectively, which were generally gained by experiments.

2.4. Basic Thought of Probabilistic Fatigue/Creep Optimization with FMERSM

The basic thought of probabilistic fatigue/creep optimization with FMERSM is illustrated below. (1) Regard material density, gas temperature, pneumatic pressure, elastic modulus, and thermal expansion coefficient as input variables, and the maximum creep stress, maximum creep strain, maximum fatigue damage, and maximum creep damage as output responses. (2) Carry out the deterministic analysis of a bladed disk based on FE models with the consideration of design parameters. (3) Obtain the fatigue damage and creep damage of a bladed disk under each load by the fatigue–creep damage equation of GH4133B discussed in Section 2.2 and Miner linear damage accumulation law introduced in Section 2.3. (4) Considering the randomness and fuzziness of input variables, enough samples of input random variables are extracted by the Latin hypercube sampling technique [30]. (5) Calculate the dynamic responses of bladed disk creep stress, creep strain, creep damage, and fatigue damage in the time domain for all input samples by FE models, and extract the maximum values of dynamic output responses as new output responses to establish the FMERSM function. (6) Perform a probabilistic analysis of a bladed disk based on the FMERSM function. (7) Complete the probabilistic fatigue/creep optimization of a turbine bladed disk by the fuzzy probabilistic optimization model with FMERSM and decoupling coordination iterative solution. The flowchart is shown in Figure 1.

Figure 1

Flow chart of reliability optimization based on the fuzzy multi-extremum response surface method (FMERSM) method.

3. Fuzzy Probabilistic Fatigue/Creep Optimization of Turbine Bladed Disk

In this section, the fuzzy probabilistic fatigue/creep optimization of a turbine bladed disk is performed with respect to the proposed FMERSM and established probabilistic optimization model in Section 2.

3.1. Parameters Preparation

With respect to the material test, the fatigue/creep material parameters θ1 and θ2 of an aero-engine turbine bladed disk with a GH4133B superalloy at the temperature of 600 °C and experimental load of 18 KN is 0.36 and 6.5, respectively. The fatigue–creep damage curve of GH4133B superalloy (Ni–Cr-based precipitation hardening-type deformation high-temperature alloy) is shown in Figure 2. In this study, we selected a 1/40 turbine bladed disk of an aero-engine as the object of study, and a GH4133B superalloy as the material of the bladed disk [31]. Density ρ, rotational speed ω, temperature T, pneumatic pressure p, elastic modulus E, and thermal expansion coefficient α are considered as fuzzy variables. Moreover, in respect of engineering practice, the length of the fuzzy region is defined as 0.05 times the mean value, as shown in Table 1. The parameters in Table 1 are assumed to obey normal distribution, and are mutually independent.

Figure 2

Curves of GH4133B fatigue–creep damage.

Table 1

Distribution characteristics of input random variables.

Random Variables	Mean	Length of Fuzzy Zone	Distribution
Density, ρ, kg·m−3	8210	410.5	Normal
Rotor speed, ω, rad·s−1	1168	58.4	Normal
Temperature T, K	873.15	43.658	Normal
Pneumatic pressure p, MPa	0.1	0.005	Normal
Elastic modulus E, MPa	163,000	8150	Normal
Thermal expansion coefficient α, ×10−6 °C−1	9.4	0.47	Normal

3.2. Deterministic Analysis of Bladed Disk

The finite element (FE) models of the blade and disk are shown in Figure 3 and Figure 4. The FE model of the blade consists of tetrahedrons with 39,547 elements, and the FE model of the disk consists of tetrahedrons with 58,271 elements. The number of cells (meshes) is required by the convergence analyses of the responses [15,32,33]. The FE basic equations of the bladed disk comprising a shape function of the tetrahedron in Equation (14) [34], geometric equation in Equation (15) [35], physical equation in Equation (16) [36], and Norton implicit creep equation in Equation (17) [37] are analyzed with regard to the means of the parameters in Table 1. From this analysis, the distributions of the creep stress and creep strain of the bladed disk are drawn in Figure 5, Figure 6, Figure 7 and Figure 8. As seen in Figure 5, Figure 6, Figure 7 and Figure 8, the maximum creep stress and maximum creep strain of the bladed disk are at the blade-root and disk tenon groove tip, respectively. where v is the volume of the tetrahedron; a, b, c and d are the related coefficients of node geometry; ε, ε and γ, γ, γ are the elastic line strains and shear strains along the x, y and z directions, respectively; [] = [σ, σ, σ, τ, τ, τ] (σ-structural stress) and [ε] = [ε, ε, ε, γ, γ, γ] are the stress and strain of the components; [] is the elastic matrix; c1, c2, and c3 stand for experimental coefficients.

Figure 3

Finite element model of the blade.

Figure 4

Finite element model of the disk.

Figure 5

Distribution of blade creep stress.

Figure 6

Distribution of disk creep stress.

Figure 7

Distribution of blade creep strain.

Figure 8

Distribution of disk creep strain.

Based on the low-cycle load spectrum described in [34], the Miner linear cumulative damage law in Equation (9) and fatigue-creep damage relation in Equation (10) of GH4133B [38] were resolved by programming in MATLAB (R2017a) simulation environment. When the number of cyclic loads is 5530, the fatigue damage D and creep damage D of bladed disk under a cyclic load are shown in Table 2.

Table 2

Results of bladed disk fatigue-creep damage.

	Fatigue Damage Df	Creep Damage Dc
Blade	0.36363	0.0039
Disk	0.40859	0.0041

3.3. FMERSM Modeling

By the Latin hypercube sampling technique [30], 100 samples on fuzzy input random variables were extracted in respect of the max creep stress, max creep strain, max fatigue damage, and max creep damage of the bladed disk, to acquire model parameters and establish the FMERSM model in Equations (18) and (19).

The FMERSM models in Equations (18) and (19) are used to perform the probabilistic fatigue/creep optimization of a bladed disk involving sensitivity analysis and reliability analysis in the following subsection.

3.4. Probabilistic Fatigue/Creep Optimization of Bladed Disk

Regarding the FMERSM models in Equations (18) and (19), the reliability sensitivity index of input random variables for a bladed disk was obtained in Table 3 and Figure 9 by sensitivity analysis with MC simulation.

Table 3

Sensitivity index of a bladed disk.

Blade			Disk
Variables	Sensitivity	Effect Probability %	Variable	Sensitivity	Effect Probability %
ρ	0.09855	10.55	ρ	0.22067	24.04
ω	0.477722	51.14	ω	0.386637	40.99
T	0.241222	25.82	T	0.252026	26.72
p	0.073895	7.91	p	–0.0121	1.28
E	0.032676	3		0.02191	2.32
α	−0.01	1.07	α	0.043787	4.64

Figure 9

Sensitivity indexes of parameters on bladed disk coupling failure.

As shown in the sensitivity analysis of a bladed disk, rotor speed w and temperature T are the main factors and greatly influence the coupling failure of a bladed disk as the two largest sensitivity degrees and effect probabilities, while the other parameters have little impact on the coupling failure of a bladed disk as smaller sensitivity degrees and effect probabilities. We extract high-sensitivity input random variables as design variables to conduct the fuzzy probabilistic fatigue/creep optimization of a bladed disk. In the main planning model, the reliability product of R1·R2 and the coupling critical damage D for the bladed disk are taken as the constraints. In the sub-planning model, the parameters (ω,T) with a high-sensitivity index are regarded as the design variables. The creep stress σ, creep strain ε, fatigue damage D, creep damage D, maximum blade reliability R1, and maximum disk reliability R2 were evaluated. The allowable comprehensive reliability of the bladed disk is R0 = 0.99. The optimal level λ* is solved by the fuzzy comprehensive evaluation method [13,23,39], and the substitution of λ* into the asymmetric fuzzy optimized conversion condition (Equation (8)). The allowable means of a bladed disk with the corresponding failure modes are shown in Table 4. The fuzzy probabilistic fatigue/creep optimization model of a bladed disk was established as illustrated in Figure 10, and the optimization models were solved by the MATLAB program and iteratively solved for all levels. The optimization results are listed in Table 5.

Table 4

Optimal level threshold and allowable mean of a bladed disk.

Blade				Disk
Optimal Level Threshold		Allowable Mean		Optimal Level Threshold		Allowable Mean
λσ_c_1*	0.3558	σ˜_c_10	677.92	λσ_c_2*	0.6008	σ˜_c_20	654.94
λεc_1*	0.8051	ε˜_c_10	2.010824	λεc_2*	0.8516	ε˜_c_20	1.0076991
λDc_1*	0.6720	D˜_c_10	0.2039	λDc_2*	0.1608	D˜_c_20	0.2041
λDf_1*	0.0260	D˜_f_10	0.96363	λDf_2*	0.8392	D˜_f_20	0.90895

Figure 10

Fuzzy probabilistic fatigue/creep optimization model of a bladed disk.

Table 5

Optimized results of a bladed disk.

Design Variable	Original Data	Optimization Results
ω, rad·s−1	1168	1055.1
T, K	873.15	788.15

3.5. FMERSM Validation

To verify the effectiveness of FMERSM, the reliability-based optimization of a bladed disk was completed with MCM and FMERSM, based on the same variables in Table 1 and computing conditions. For dynamic probabilistic analyses under different simulations (102, 103, 104, and 105), the computing time and reliability degrees of a bladed disk are listed in Table 6. The optimization results of object functions under different simulations are listed in Table 7.

Table 6

Dynamic probabilistic computational results with different methods. FMERSM: fuzzy multi-extremum response surface method, MC: Monte Carlo.

Number of Samples	Computational Time, s		Reliability Degree %		Precision of FMERSM
	MC Method	FMERSM	MC Method	FMERSM
102	32400	0.203	99	98.6	0.996
103	72000	0.279	99.7	99.5	0.998
104	432000	0.437	99.83	99.60	0.9977
105	-	4.43	-	99.962	-

Table 7

Results of bladed disk optimization design with different methods.

Objective Function	Before Optimization	MCM		FMERSM
		After Optimization	Reduction	After Optimization	Reduction
σ_c_1, MPa	607.92	548.66	9.8%	487.08	19.9%
σ_c_2, MPa	454.94	423.03	7%	368.8	18.93%
ε_c_1, m/m	0.010824	0.001298	88%	0.0073988	31.64%
ε_c_2, m/m	0.0076991	0.0076629	0.47%	0.0065652	14.77%
Df1	0.36363	0.30452	16.25%	0.24822	31.74%
Df2	0.40859	0.39375	4.52%	0.29315	28.3%
D_c_1	0.0039	0.0036	7.69%	0.0025	35.9%
D_c_2	0.0041	0.00409	0.24%	0.0032	21.95%
R	95	99.515	-	99.635	-

As shown in Table 6, the following conclusions were obtained from the probabilistic failure analysis of a bladed disk. (1) The MC method does not have computing time at 105 simulations, because the MC method cannot perform the calculation for a too-large computational burden for a probabilistic analysis of bladed disk FE models. Thus, it is inefficient for the MC method to conduct the design analysis of a complex structure with large-scale simulations. (2) The time-cost for the probabilistic analysis of a bladed disk increases with the increase of MC simulations. (3) The time consumption of the FMERSM is far less than that of the MC method for the same number of simulations. For instance, the FMERSM only spends 0.437 s for 10,000 simulations, which is only about 1/106 that of the MC method. Meanwhile, the strength of the FMERSM in time computation is more obvious with increasing simulations. Thus, it is demonstrated that the efficiency of the FMERSM is far higher than that of the MC method in calculation, and the FMERSM is an efficient approach replacing FE models for the probabilistic analysis of a complex structure with many components or multi-failure modes. (4) For the same simulations, the reliability degrees of bladed disk coupling dynamic failure probability with FMERSM are almost consistent with those of the MC method. Moreover, the reliability degree of the bladed disk increases and becomes higher with the rise of simulations. It is illustrated that more precise results such as the reliability degree can be gained by increasing the number of MC simulations against the response surface models, for structure design analysis from a probabilistic perspective.

As revealed in Table 7, summarized from the probabilistic fatigue/creep optimization of a bladed disk, the creep stress, creep strain, fatigue damage, and creep damage of the blade in respect of FMERSM are reduced by 19.9%, 18.93%, 31.64%, and 14.77%, respectively. Meanwhile, the MC method reduces the creep stress, creep strain, fatigue damage, and creep damage of the disk by 9.8%, 7%, 88%, and 0.47%, respectively. The comprehensive reliability index of the bladed disk was increased from 99.515% to 99.635%. It is verified that the FMERSM is workable for the fuzzy probabilistic fatigue/creep optimization of complex structures, similar to a turbine bladed disk.

In summary, the developed FMERSM has high modeling precision and simulation efficiency for the comprehensive reliability optimization design for multi-component structures with multi-failure modes.

4. Conclusions

The objective of this study is to develop a high-efficient reliability-based optimization method, called the fuzzy multi-extremum response surface method (FMERSM), for the probabilistic fatigue/creep coupling optimization of a turbine bladed disk. This paper has investigated the theory and modeling of FMERSM, and gives the procedure of probabilistic optimization of a multi-component structure with multi-failure modes for the fuzzy probabilistic fatigue/creep optimization of a turbine bladed disk with the considerations of the correlation of the failure modes and the fuzziness of the constraint boundary conditions. Through the works in this study, some conclusions are summarized as follows:

(1) With regard to the probabilistic failure analysis of a bladed disk, we find that the FMERSM costs less analytical time (0.437 s for 10,000 simulations), and thus has high computational efficiency relative to the Monte Carlo (MC) method (432,000 s for 10,000 simulations), but has an acceptable computational precision (99.77%) of the reliability degree, which is almost consistent with the FE method based on MC simulation with a reliability degree of 0.9983. Moreover, the strengths of the proposed FMERSM in modeling and simulation become more obvious with the increase of simulations.

(2) In terms of the probabilistic fatigue/creep optimization of a bladed disk, it is illustrated that the developed FMERSM is more workable than the MC method. The reason is that the optimal parameters, including design parameters and optimization objects, are preferable by larger reductions (19.9%, 18.93%, 31.64%, and 14.77% for the creep stress, creep strain, fatigue damage, and creep damage of the blade, respectively), and a higher reliability degree of 99.635%.

The efforts of this paper provide a useful way for high-precise modeling and high-efficient simulation for the fuzzy comprehensive probabilistic optimization of multi-failure/multi-component structures, because the accuracy of the model is close to that of the MC method, while the calculation time is only 1/106. Meanwhile, this work enriches the theory of mechanical reliability.