Next Article in Journal
Effect of Intermetallic Compounds on the Thermal and Mechanical Properties of Al–Cu Composite Materials Fabricated by Spark Plasma Sintering
Next Article in Special Issue
PSO-BP Neural Network-Based Strain Prediction of Wind Turbine Blades
Previous Article in Journal
Prediction and Sensitivity Analysis of Bubble Dissolution Time in 3D Selective Laser Sintering Using Ensemble Decision Trees
Previous Article in Special Issue
Application of a New, Energy-Based ΔS* Crack Driving Force for Fatigue Crack Growth Rate Description
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Reliability-Based Low Fatigue Life Analysis of Turbine Blisk with Generalized Regression Extreme Neural Network Method

1
School of Mechanical and Power Engineering, Harbin University of Science and Technology, Harbin 150080, China
2
Department of Aeronautics and Astronautics, Fudan University, Shanghai 200433, China
3
School of Computer Science and Technology, Beihang University, Beijing 10191, China
*
Author to whom correspondence should be addressed.
Materials 2019, 12(9), 1545; https://doi.org/10.3390/ma12091545
Submission received: 23 March 2019 / Revised: 6 May 2019 / Accepted: 8 May 2019 / Published: 10 May 2019

Abstract

:
Turbine blisk low cycle fatigue (LCF) is affected by various factors such as heat load, structural load, operation parameters and material parameters; it seriously influences the reliability and performance of the blisk and aeroengine. To study the influence of thermal-structural coupling on the reliability of blisk LCF life, the generalized regression extreme neural network (GRENN) method was proposed by integrating the basic thoughts of generalized regression neural network (GRNN) and the extreme response surface method (ERSM). The mathematical model of the developed GRENN method was first established in respect of the LCF life model and the ERSM model. The method and procedure for reliability and sensitivity analysis based on the GRENN model were discussed. Next, the reliability and sensitivity analyses of blisk LCF life were performed utilizing the GRENN method under a thermal-structural interaction by regarding the randomness of gas temperature, rotation speed, material parameters, LCF performance parameters and the minimum fatigue life point of the objective of study. The analytical results reveal that the reliability degree was 0.99848 and the fatigue life is 9419 cycles for blisk LCF life when the allowable value is 6000 cycles so that the blisk has some life margin relative to 4500 cycles in the deterministic analysis. In comparison with ERSM, the computing time and precision of the proposed GRENN under 10,000 simulations is 1.311 s and 99.95%. This is improved by 15.18% in computational efficiency and 1.39% in accuracy, respectively. Moreover, high efficiency and high precision of the developed GRENN become more obvious with the increasing number of simulations. In light of the sensitivity analysis, the fatigue ductility index and temperature are the key factors of determining blisk LCF life because their effect probabilities reach 41% and 26%, respectively. Material density, rotor speed, the fatigue ductility coefficient, the fatigue strength coefficient and the fatigue ductility index are also significant parameters for LCF life. Poisson’s ratio and elastic modulus of materials have little effect. The efforts of this paper validate the feasibility and validity of GRENN in the reliability analysis of blisk LCF life and give the influence degrees of various random parameters on blisk LCF life, which are promising to provide useful insights for the probabilistic optimization of turbine blisk LCF life.

1. Introduction

As a heat-end core component of an aeroengine, a turbine blisk endures complex alternating loads due to operation in a severe environment with high temperatures and high rotation speeds. In this case, it is easy to produce large plastic deformation for blisk and to induce the low cycle fatigue (LCF) failure of blisk [1,2]. Most of the parameters that significantly effect blisk LCF failure have some randomness [3]. To improve the safety and reliability of a turbine blisk to ensure the high performance of an aeroengine, it is important to study blisk LCF life reliability from a probabilistic perspective [4,5,6,7,8].
The LCF life of structures has been widely investigated. Sun et al. established a nonlinear model for LCF life of a steam turbine rotor under a temperature-stress coupling field by considering the relationship between cyclic stress and strain and validated the model to be accurate and reasonable in describing damage accumulation [9]. Letcher et al. proposed an energy-based critical fatigue life prediction approach, which derived the approximate failure cycle index from the ratio of the total accumulation of energy in the fracture process to the one-cycle strain energy [10]. Bargmann et al. discussed the full-probability quick integral algorithm based on the Coffin-Manson-Neuber local strain-fatigue theory [11]. Zhu et al. discussed the probabilistic LCF life prediction of a turbine disk under uncertainties [12,13]. Viadro et al. studied the reliability of stiffened bending plates [14]. Repetto et al. discussed the role of parameter uncertainty in the damage prediction of the alongwind-induced fatigue and long term simulation of wind-induced fatigue loadings [15,16]. Most of the above work was conducted based on numerical simulation methods (or-called direct simulation methods) with Monte Carlo (MC) simulation [15,16,17,18,19,20]. Generally, the direct simulation methods are powerful for the deterministic analyses of component LCF life. However, for the probabilistic analyses of component LCF life with thousands of iterations and MC simulations, it is unbelievable to efficiently perform blisk fatigue life analysis owing to excess computational burden (loads) and unacceptable computational efficiency; although this method has satisfactory computing precision against engineering practice. Therefore, it is urgent to seek an alternative effective method for direct methods to address this issue.
In respect of the in-depth investigation of structural fatigue probabilistic analyses, the response surface method (RSM, also called surrogate model method) is indeed an alternative method to direct simulation methods [21,22,23,24,25,26]. With the development of structural reliability theory and methods, various surrogate methods have emerged [27,28,29]. To improve the computational efficiency and accuracy of RSM for complex structural reliability analysis, Bai et al. proposed a distributed collaborative response surface method for the mechanical dynamic assemble reliability analysis of aeroengine high pressure turbine blade-tip clearance [30]. Hurtado et al. proposed a highly efficient surrogate method, a support vector machine, for structural reliability analysis with small samples [31]. Zhang et al. developed an extremum response surface method (ERSM) in respect of the extreme thought, to address the transient problem in the dynamic reliability analysis of a flexible mechanism and validated the ERSM to be precise for the reliability analysis of a flexible manipulator [32]. Lu et al. developed an improved Kriging method by integrating the Kriging algorithm and ERSM for the reliability and sensitivity analyses of a compressor blisk regarding multiple failure (deformation failure, stress failure and strain failure) modes [33]. From the efforts in References [32,33], it can be seen that ERSM has the potential to handle the transient problem in structural dynamic probabilistic analyses with a high simulation accuracy and efficiency, which provides useful insight into the process of the reliability analysis of blisk LCF life with the consideration of aeroengine operating conditions. For another, the developed ERSM does not satisfy the requirement of engineering in computing precision, derived from the weakness in processing the involved nonlinear probabilistic analyses.
With the development of neural network technology recently, the nonlinear problem was skillfully addressed by developing a generalized regression neural network (GRNN) due to the strong nonlinear mapping capability and robustness [34]. Zhao et al. established an all-purpose regression neural network model based on a freight volume condition and validated the effectiveness of this model in freight volume prediction by modeling adaptive training and extrapolation evaluation in terms of historical statistical data of freight volume and related samples and economic indicators [35]. Li et al. fused the drosophila optimization algorithm and GRNN to build the prediction model of power loads for power load prediction and this model had a strong nonlinear fitting ability [36]. Sun et al. compared the GRNN model with the back propagation neural network (BPNN) model based on air quality prediction and the GRNN method needed less training time and had better stability, a higher fitting precision as well, compared to the BPNN model [37]. Wang et al. validated the strengths of the GRNN method again by waveguide orientation [38]. Therefore, the GRNN method has been comprehensively verified to be highly computationally precise and efficient.
To effectively perform the reliability and sensitivity analyses of a turbine blsik LCF life, the generalized regression extremum neural network (GRENN) method is proposed in this paper; by integrating the transient procession ability of ERSM and the nonlinear mapping and small samples of GRNN, to collectively ensure and improve computing precision and efficiency. The reliability analysis of a turbine blisk LCF life was implemented based on the the developed GRENN, by considering random input variables of temperature, rotation speed, material parameters (density, Poisson’s ratio and elastic modulus) and fatigue performance parameters (fatigue ductility coefficient, fatigue strength coefficient and fatigue ductility index and fatigue strength index) as well as the output response of the minimum fatigue life. The developed GRENN method was validated by comparison with the MC method and ERSM.

2. Basic Theory

2.1. Mathematical Model of Low Cycle Fatigue Life

The Mason-Coffin equation indicates the strain-fatigue life equation, which expresses the relationship between strain and the fatigue life of materials [39], i.e.,
Δ ε 2 = σ f E ( 2 N f ) b + ε f ( 2 N f ) c
where Δ ε is the total strain of specific structure; E indicates the elasticity modulus; σ f denotes the fatigue strength coefficient; ε f is the fatigue ductility coefficient; b indicates the fatigue strength exponent; c stands for the fatigue ductility exponent; Nf expresses the LCF life.
Considering the mean stress, σm, inducted by complex loads during aeroengine operation, Equation (1) can be rewritten by the Morrow correction equation [40], i.e.,
Δ ε 2 = σ f σ m E ( 2 N f ) b + ε f ( 2 N f ) c .
Considering many cyclic loads, the LCF life can be classically computed by the line damage accumulation (Miner) law, i.e.,
D = i = 1 r n i N i
in which D indicates the fatigue damage; r is the number of loading levels; ni denotes the cyclic number under the ith loading level; Ni is the fatigue life corresponding to the ith loading level.

2.2. Mathematical Model of Extremum Response Surface Method

To effectively process the transient problem in the dynamic reliability analysis of blisk LCF life involving nonlinear and transient features of numerous parameters, i.e., gas temperature, rotation speed, material parameters (density, Posion’s ratio, elasticity modulus) and fatigue performance parameters (fatigue strength coefficient, fatigue ductility coefficient, fatigue strength exponent and fatigue ductility exponent), the ERSM proposed in Reference [32] was adapted by simplifying the response process of the LCF life as an extreme value (maximum value or minimum value) in an analytical time domain. When X and ye were used to indicate the input parameters set and the output extremum response, the ERSM model ye(X) of the dynamic system [32] can be written as
y e ( X ) = f ( X ) = { y e ( j ) ( X ( j ) ) }
where X ( j ) is the jth group of the input samples; y e ( j ) ( X ( j ) ) indicates the output extremum response during a time domain.
In previous studies, most of the ERSM models were built based on polynomials [33] and these models are usually inefficient in model fitting because the polynomials are unworkable for highly nonlinear problems and the large computing burden (requiring a large number of samples for modeling) required for blisk fatigue life probabilistic analysis. Thus, the GRNN method, with a strong nonlinear mapping ability and robustness, was applied by combining ERSM in this paper to address the issues of modeling precision and efficiency, which result from nonlinearity and large samples.

2.3. Mathematical Model of Generated Regression Extremum Neural Network Method

GRNN is a feedforward neural network model based on nonlinear regression theory, including input layer (first layer), hidden layer (middle layer) and output layer (last layer), as shown in Figure 1.
By inputting train samples into the input layer, the input matrix X and output matrix (denoted by T) are expressed by
X = [ x 11 x 12 x 1 Q x 21 x 22 x 2 Q x R 1 x R 2 x R Q ] ,   T = [ t 11 t 12 t 1 Q t 21 t 22 t 2 Q t S 1 t S 2 t S Q ]
where xji (i = 1,2, …, Q; j = 1,2, …, R) is the jth input sample in the ith group of training samples; tji indicates the jth output sample in the ith group of training samples; R is the number of input variables; S is the number of output variables; Q is the number of samples in the training set.
The number of neurons in the hidden layer was equal to the number of samples in the training set, the layer weight function was the Euclidean distance function (expressed in ||dist||), and the implicit layer weight matrix is calculated as follows:
L W 1 , 1 = X T
The threshold for Q hidden layer neural units is b:
b = [ b 1 , b 2 , b Q ] T
in which b 1 = b 1 = = b Q = 0.8326 σ , σ is the smooth factor of the Gauss function.
The transfer function of the hidden layer is usually based on the Gaussian radial basis function. The number of neurons in the hidden layer Q is equal to the number of training samples and each neuron corresponds to one training sample. When the weight matrix and threshold value of the hidden layer neural unit were determined, the output a i j of the ith hidden layer neuron is
a i j = e x p ( 0.8326 σ L W 1 . i x j 2 ) ,   j = 1 , 2 , , Q   ; i = 1 , 2 , , Q
in which L W 1 . i = [ x h 1 , x h 2 , , x h R ] T (h=1,2,   , Q ) is the vector of the ith implicit layer weight matrix LW1,1; x j = [ x j 1 , x j 2 , , x j R ] T is the vector of the jth training samples. let a j = [ a 1 j , a 2 j , , a i j , a Q j ] , which is the output vector of Q nerve cells corresponding to the jth group of input samples.
Regarding the connection weight LW2,1 between the hidden layer and output layer as the output matrix of the training set of samples, which is denoted by T, i.e.,
L W 2 , 1 = T .
The output layer is the third layer of GRNN. Based on GRNN Equations (8) and (9), vector nj can be computed by
n j = L W 2 , 1 [ a j ] T i = 1 Q a i j .
With regard to the line transfer function   y j = p u r e l i n ( n j ) of nj, the mathematical model of GRNN for the response of the jth group of training samples is expressed by
y j = p u r e l i n ( n j ) = L W 2 , 1 [ a j ] T i = 1 Q a i j
where exp is a natural exponential function.
With respect to the format in Equation (4), the mathematical model of GRENN is
y m i n j = M i n { L W 2 , 1 [ a j ] T i = 1 Q a i j } .

2.4. Reliability Sensitivity Analyses Approaches with GRENN Model

Assuming that y * is the allowable LCF life and y m i n j is the performance function of the structural fatigue life, the limit state function of LCF life is derived as [24]
Z = y m i n j y * .
In Equation (13), Z > 0 indicates that the blisk structure is secure, while Z < 0 reveals a failure. When random input variables are independently mutual, their means and variance are denoted by μ = [ μ 1 , μ 2 μ n ] and D = [ D 1 , D 2 D n ] , respectively, we can gain [18]
{ E ( Z ) = μ Z ( μ 1 , μ 2 , , μ n ; D 1 , D 2 , , D n ) D ( Z ) = D Z ( μ 1 , μ 2 , , μ n ; D 1 , D 2 , , D n )
in which E ( Z ) is mean function and D ( Z ) is variance function.
When the limit state function of the structural LCF life (Equation (13)) obeys a normal distribution, the reliability degree Pr is expressed as [25]
P r = Φ ( μ Z D Z )
where µz is the mean matrix of a limit state function Z; Dz is the variance matrix of a limit state function.
The sensitivity reflects the level of sensitivity of the input random variables on the failure probability of a structural system response, which is promising to determine the extent to which these parameters effect the response and then provide a useful guide for structural design and optimization [41].
With the proposed GRENN method, the sensitivity degree can be determined by the mean matrix µ and variance D of input random variables [42], i.e.,
P r μ T = P r ( μ Z D Z ) ( ( μ Z D Z ) μ Z μ Z μ T + ( μ Z D Z ) D Z μ Z μ T ) ; P r D T = P r ( μ Z D Z ) ( ( μ Z D Z ) μ Z μ Z D T + ( μ Z D Z ) D Z μ Z D T )
in which
{ P r ( μ Z / D Z ) = P r , ( μ Z / D Z ) μ Z = 1 D Z , ( μ Z / D Z ) D Z = μ Z 2 D Z 3 2 μ Z μ T = [ μ Z μ 1 , μ Z μ 2 , , μ Z μ n ] T μ Z D T = [ μ Z D 1 , μ Z D 2 , , μ Z D n ] T D Z μ T = [ D Z μ 1 , D Z μ 1 , , D Z μ 1 ] T D Z D T = [ D Z D 1 , D Z D 1 , , D Z D 1 ] T .
In respect of the GRENN method and thermal-structure coupling, the flowchart of the blisk LCF life reliability analysis is drawn in Figure 2 and its basic procedure is described below.
Step 1: 
Build the finite element (FE) model of blisk in a workbench environment;
Step 2: 
Consider the means of the input random variables (i.e., gas temperature, rotation speed, material parameters and fatigue performance parameters) and set boundary conditions to conduct the blisk FE analysis under the interaction of heat load, centrifugal load and then gain the minimum fatigue point as the design point of the blisk reliability design.
Step 3: 
Extract small samples of the input random variables using the Latin hypercube sampling (LHS) method and perform FE analyses for each group of samples to gain the output responses (blisk LCF life) and extract the minimum values of the responses as a training sample set by combining the input samples.
Step 4: 
Training the GRENN model by computing the optimal smooth factors, radial basis function and connection weights with the cross validation method [26], through the normalization of training samples.
Step 5: 
Structure of the limit state function of blisk LCF life with the established GRENN model.
Step 6: 
Check the precision of the GRENN model. If unacceptable, return to Step 4; if acceptable, conduct Step 7.
Step 7: 
Calculate the reliability degree and sensitivity degree of the fatigue life and input variables, by conducting the reliability and sensitivity analyses of blisk LCF life with thermal-structure coupling, through a large number of samples extracted by the MC method.

3. Reliability and Sensitivity Analyses of Blisk Low Cycle Fatigue Life

3.1. Random Variables Selection

In this study, we selected the high-pressure turbine blisk of an aeroengine with the high-temperature GH4133 as the object of study. In fact, the uncertainty and randomness of some parameters are the basic nature in blisk LCF life design and prediction [15]. By comprehensively regarding the engineering practice, the exiting data and the basic properties of parameter uncertainty studied by Repetto, et al. [15], the probabilistic analysis of the blisk LCF life was performed by the randomness of numerous reasonably-selected parameters, such as rotation speed ω, gas temperature T, material density ρ, heat conductivity λ, elasticity modulus E, fatigue strength efficient σ′f, fatigue ductility coefficient ε′f, fatigue strength index b and fatigue ductility index c. To simplify the calculation by combining engineering practices [43] and the present data, the selected variables were summed to be independent mutually and obey normal distributions. The distributions of the variables are listed in Table 1.

3.2. Deterministic Analysis of Blisk Low Cycle Fatigue Life

For the static analysis of the blisk, the blisk stress inducted by aerodynamic loads can be ignored because it is far less than that caused by the centrifugal load and heat load [43]. The deterministic analysis of the blisk was completed by regarding the interaction of temperature and centrifugal loads, under a workbench 16.0 environment in the computer with a central processing unit (CPU) mode of Xeon E5-2630V3 (Intel Corporation, Santa Clara, CA, USA) and RAM (Intel Corporation) of 64 GB. Due to the symmetry of the blisk, we selected 1/40 of the whole blisk for analysis to reduce the computational burden [44]. The FE models are shown in Figure 3, with 31,380 nodes and 17,111 elements. The thermodynamic analysis of the blisk was implemented in which the heat energy of a high temperature gas is transferred to the surface of the blisk according to the heat conduction law and heat convection. In light of thermodynamic theory, the temperature distribution on the blisk surface can be calculated by the empirical formula, i.e.,
T = T a + ( T a T b ) ( R m R a m R b m R a m )
in which T a is the temperature at blisk-root; R a is the radius of blisk-root edge;   T b is the temperature at blisk-tip;   R b is the radius of blisk-tip; R is the radius of blisk in a different position; m = 2 was determined for the high temperature alloy GH4133B [39].
By the displacement constraint of the blisk’s inner diameter to restrict the degrees of freedom in the directions x, y and z, the deterministic analysis of the blisk was finished based on the means of the input variables in Table 1. The distributions of temperature, equivalent stress and equivalent strain are shown in Figure 4a–c. As seen in Figure 4a–c, the maximum stress of the blisk was 1 057.7 Mpa on the blade-root and the minimum strain was 8.142 7 × 10−3 m/m. Therefore, the node of the maximum strain on the blade-root was selected as the object of study for the blisk LCF life analysis. In terms of the Mason-Coffin formula in Equation (2) and the Miner line accumulative damage rule in Equation (3), the fatigue life values at the max-strain point of the blisk areshown in Figure 4d. It is illustrated in Figure 4d that the minimum fatigue life was 8900.6 cycles. In respect of the double safety coefficients in engineering, the LCF of the blisk should be about 4450 cycles based on the deterministic analysis.

3.3. Low Cycle Fatigue Life Models of Blisk with GRENN Method

With regard to the distribution of the input random variables in Table 1, 150 samples (a small batch) were extracted by LHS technology. Based on these samples and FE analyses, the corresponding output responses (minimum LCF lives) were computed as the samples together with the extracted input samples. One hundred and twenty groups of samples were selected from the pool of training samples as training samples and the remain 30 groups of samples were selected as the test samples for the GRENN model.
Regarding the Gauss function as a transfer function in the hidden layer, the implicit layer weight LW1.1 of the hidden layer was computed using the Euclidean distance method. The outputs of GRENN training were taken as the connection weights LW2,1 between the hidden layer and the output layer. The original samples data should be normalized for each parameter. The normalized data were adopted to train the GRENN and then to gain the parameters of GRENN (the implicit layer weight LW1.1 of the hidden layer, the connection weights LW2,1 and the smooth factor σ) by the cross validation method [35], in which b and LW1.1 and LW1,2 (computed by Equations (6), (7) and (9)) are summarized in Equation(19). By inputting the values of these parameters into Equation (12), the GRENN model can be gained. The remaining 30 groups of samples were employed to test the established GRENN model. The prediction results are shown in Figure 5. From Figure 5, it can be seen that the predicted data were almost consistent with the true sample data, which indicates a high prediction precision for the developed GRENN model.
{ L W 1 , 1 = [ 0.4189 0.0946 0.9054 0.9459 0.7973 0.8912 0.5646 0.7415 0.2517 0.1020 1.0000 0.5646 0.7415 0.7959 0.5510 0.7852 0.7315 0.8792 0.1678 0.4765 0.2245 0.8639 0.2789 0.6190 0.3878 0.5839 0.3020 0.2483 1.0000 0.7315 0.6510 0.2752 0.6644 0.1678 0.5973 0.6892 0.8514 0.6216 0.9865 0.4865 0.0470 0.5302 0.7584 0.6107 0.0336 0.1757 0.8108 0.7432 0.2838 0.3919 ] 10 × 120 T L W 2 , 1 = [   0.9061 0.8701   0.9627 0.9968 0.9400 ] 1 × 120 b = [ 2.8710   2.8710     2.8710   2.8710   2.8710 ] 1 × 120 T

3.4. Reliability Analysis of Blisk Low Cycle Fatigue Life with GRENN Model

In this subsection, the reliability analysis of the blisk LCF life with the GRENN model was performed by 10,000 simulations with the MC method. The historical simulation diagram and histogram of the blisk LCF life are shown in Figure 6. As shown in Figure 6, the blisk minimum fatigue life followed a normal distribution with a mean of 9419 cycles and a standard deviation 967 cycles. As the allowable fatigue life y* = 6 000 cycles, the reliability degree Pr of the blisk LCF life was 0.99848 in line with Equations (13) and (15). The gained reliability degree basically catered for blisk design in engineering. In this case, the obtained fatigue life of a blisk was 6000 cycles in respect to the reliability analysis. However, the minimum LCF life of the deterministic analysis was 8900.6 cycles, as shown Figure 4d. In respect to the double safe coefficients, the safe fatigue life of blisk design was about 4450 cycles in engineering in line with the deterministic analysis. Therefore, it is revealed that the deterministic analysis method is backward-looking relative to ~4450 cycles of the probabilistic analysis method for blisk LCF life prediction at 6000 cycles, because 4450 cycles was far less than 6000 cycles.

3.5. Sensitivity Analysis of Blisk Low Cycle Fatigue Life with GRENN Method

Sensitivity reflects the level of sensitivity of the input random variables on blisk reliability, which is helpful to find the major impact factors and then guide structural design. Sensitivity involves the sensitivity degree and the effect probability. The sensitivity degree is defined by the effect of the input parameters on the output response with positive and negative signs. A positive sign indicates the input parameter was positively correlated with the output response and vice versa for a negative sign. Effect probability is defined as the ratio of the sensitivity degree of one input parameter to the total sensitivity degree of all input parameters. In terms of Equations (12)–(17), the sensitivity results are listed in Table 2 and Figure 7.
As demonstrated in Table 2 and Figure 7, the fatigue ductility index c and gas temperature T were two major influencing parameters because their effect probabilities and sensitivity degree were 41.3% and 0.27929, as well as 26.16% and −0.176122, respectively. Other parameters play small effects on blisk reliability. Therefore, T and c should be considered and controlled in blisk design with priority.

3.6. Validation of GRENN

To validate the effectiveness and validation of the GRENN method, the MC method and ERSM were used in the reliability analyses of blisk LCF life under different simulations based on the same computation conditions and random variables. The computing time and reliability degrees are presented in Table 3 and Table 4. In Table 3 and Table 4, the precision for the method Dp was computed under 10,000 simulations, by
D p = 1 | γ a γ m | γ a × 100 %
in which γa is the reliability degree of the MC method; γm indicates the reliability degree of ERSM or the GRENN method.
As shown in Table 3, with the increasing number of simulations, the computing time increases for the MC method, ERSM and the GRENN method. For the MC method, the simulations larger than 10,000 require an excessive computational burden so that the MC method is unworkable for such large simulations. However, ERSM and the GRENN method only take a few seconds and thus can breezily implement simulations from 100 to 1,000,000. Relative to ERSM, the developed GRENN method spends less time and is highly computationally efficient, and that the strength of the GRENN method becomes more obvious with an increase in the number of simulations. For instance, under 10,000 simulations, the GRENN method reduces the computing time by 0.048 s and improves the computational efficiency by 3.843% relative to the ERSM, while the simulation time is reduced by 2.138 s and the efficiency is improved by 50.42%. Therefore, it is revealed that the proposed GRENN method has a strong computing power and is highly computationally efficient in probabilistic simulations. Meanwhile, the potential of high efficiency becomes stronger and more simulations are required.
In Table 4, the reliability degree computed by the MC method under 10,000 simulations is regarded as a reference. In this case, we find that the reliability degrees of ERSM and the GRENN method were 0.9824 and 0.9923 under 10,000 simulations, and their computational precisions were 98.56% and 99.95% so that the GRENN method improves the precision by 1.39%. Additionally, with an increasing number of simulations, the reliability degree of a turbine blisk increases and the developed GRENN method is more accurate than ERSM.
Therefore, the developed GRENN method is highly computationally precise and efficient and the strengths become more obvious for more simulations.

4. Conclusions

The aim of this paper was to propose a new reliability analysis method, i.e., the generalized regression extreme neural network (GRENN) method, for the reliability analysis of blisk LCF life, to improve the life and performance of turbine blisks. The developed GRENN absorbed the strengths of a generalized regression neural network (GRNN) in nonlinear mapping and small sample-based modeling, and the extremum response surface method (ERSM) for handling the transient problem of the dynamic reliability analysis of blisk LCF life. Through this study, some conclusions are summarized as follows:
(1)
The reliability degree of blisk LCF life was 0.99848 when the life allowable value was 6000 cycles. Relative to 4450 cycles acquired from the deterministic analysis after considering the double coefficient of a safe life, the LCF (6000 cycles to ensure a reliability degree of 0.99848) of the blisk obtained from the reliability design had enough life margin (about 1550 cycles) to ensure the operation of the blisk structure.
(2)
From the sensitivity analysis of a blisk, the fatigue ductility index c and gas temperature T played key roles in blisk LCF life evaluation and design. T and c were positively and negatively correlated with blisk life, respectively. The conclusions can significantly guide the optimization and design of blisk LCF life.
(3)
Through the comparison of the methods, it is demonstrated that the developed GRENN method is far better than ERSM in modeling precision and computing efficiency and is basically consistent with the MC method. Moreover, the strengths of the GRENN method become more obvious with the increasing number of simulations. It is fully supported that the proposed GRENN method is a high-accuracy and high-efficiency method to address the key questions of nonlinearity, transients and large sample-based modeling.
In summary, the efforts of this paper provide a promising method (GRENN method) for the nonlinear dynamic reliability analysis of complex structures and enrich and develop mechanical reliability theory.

Author Contributions

Conceptualization, C.Z. and C.F.; methodology, C.Z.; software, J.W.; validation, J.W., H.J. and C.F.; formal analysis, J.W.; investigation, H.J.; resources, C.Z.; data curation, C.F.; writing—original draft preparation, C.Z. and W.T.; writing—review and editing, C.F.; visualization, H.J.; supervision, C.F. and W.T.; project administration, C.Z.; funding acquisition, C.Z. and C.F.

Funding

This research was funded by National Natural Science Foundation of China, grant numbers 51275138, 51605016 and 51475025, Start-up Research Funding of Fudan University, grant number FDU38341, Shanghai Pujiang Program, grant number PJ20191334 and The APC was funded by 51605016.

Conflicts of Interest

The authors declare that there are no conflict of interest regarding the publication of this article.

References

  1. Bin, B.; Bai, G.C. Dynamic probabilistic analysis of stress and deformation for bladed disk assemblies of aeroengine. J. Cent. South Univ. 2014, 21, 3722–3735. [Google Scholar]
  2. Gao, H.F.; Wang, A.; Bai, G.C.; Wei, C.M.; Fei, C.W. Substructure-based distributed collaborative probabilistic analysis method for low-cycle fatigue damage assessment of turbine blade-disk. Aerosp. Sci. Technol. 2018, 79, 636–646. [Google Scholar] [CrossRef]
  3. Scott-Emuakpor, O.; George, T.; Cross, C.; Shen, M. Multi-axial fatigue-life prediction via a strain-energy method. AIAA J. 2010, 48, 63–72. [Google Scholar] [CrossRef]
  4. Hou, N.X.; Wen, Z.X.; Yu, Q.M.; Yue, Z.Y. Application of a combined high and low cycle fatigue life model on life prediction of SC blade. Int. J. Fatigue 2009, 31, 616–619. [Google Scholar] [CrossRef]
  5. Zhu, S.P.; Liu, Q.; Peng, W.; Zhang, X.C. Computational-experimental approaches for fatigue reliability assessment of turbine bladed disks. Int. J. Mech. Sci. 2018, 142–143, 502–517. [Google Scholar] [CrossRef]
  6. Zhu, S.P.; Liu, Q.; Zhou, J.; Yu, Z.Y. Fatigue reliability assessment of turbine discs under multi-source uncertainties. Fatigue Fract. Eng. Mater. Struct. 2018, 41, 1291–1305. [Google Scholar] [CrossRef]
  7. Zhu, S.P.; Liu, Q.; Lei, Q.; Wang, Q.Y. Probabilistic fatigue life prediction and reliability assessment of a high pressure turbine disc considering load variations. Int. J. Damage Mech. 2018, 27, 1569–1588. [Google Scholar] [CrossRef]
  8. Zhu, S.P.; Huang, H.Z.; Smith, R.; Ontiveros, V.; He, L.P.; Modarres, M. Bayesian framework for probabilistic low cycle fatigue life prediction and uncertainty modeling of aircraft turbine disk alloys. Probab. Eng. Mech. 2013, 34, 114–122. [Google Scholar] [CrossRef]
  9. Sun, Y.; Hu, L.S. Low cycle fatigue life prediction of a 300MW steam turbine rotor using a new nonlinear accumulation approach. In Proceedings of the 24th Chinese Control and Decision Conference (CCDC), Taiyuan, China, 23–25 May 2012; pp. 3913–3918. [Google Scholar]
  10. Letcher, T.; Shen, M.H.H.; Scott-Emuakpor, O.; George, T.; Cross, C. An energy-based critical fatigue life prediction method for AL6061-T6. Fatigue Fract. Eng. Mater. Struct. 2012, 35, 861–870. [Google Scholar] [CrossRef]
  11. Bargmann, H.; Rüstenberg, I.; Devlukia, J. Reliability of metal components in fatigue: A simple algorithm for the exact solution. Fatigue Fract. Eng. Mater. Struct. 2010, 17, 1445–1457. [Google Scholar] [CrossRef]
  12. Zhu, S.P.; Foletti, S.; Beretta, S. Probabilistic framework for multiaxial LCF assessment under material variability. Int. J. Fatigue 2017, 103, 371–385. [Google Scholar] [CrossRef] [Green Version]
  13. Zhu, S.P.; Huang, H.Z.; Peng, W.; Wang, H.K.; Mahadevan, S. Probabilistic Physics of Failure-based framework for fatigue life prediction of aircraft gas turbine discs under uncertainty. Reliab. Eng. Syst. Saf. 2016, 146, 1–12. [Google Scholar] [CrossRef]
  14. Viadero, F.; Bueno, J.I.; de Lacalle, L.L.; Sancibrian, R. Reliability computation on stiffened bending plates. Adv. Eng. Softw. 1995, 20, 43–48. [Google Scholar] [CrossRef]
  15. Pagnini, L.; Repetto, M.P. The role of parameter uncertainties in the damage prediction of the alongwind-induced fatigue. J. Wind Eng. Ind. Aerodyn. 2012, 104–106, 227–238. [Google Scholar] [CrossRef]
  16. Repetto, M.P.; Torrielli, A. Long term simulation of wind-induced fatigue loadings. Eng. Struct. 2017, 132, 551–561. [Google Scholar] [CrossRef]
  17. Marseguerra, M.; Zio, E. The cell-to-boundary method in Monte Carlo-based dynamic PSA. Reliab. Eng. Syst. Saf. 1995, 48, 199–204. [Google Scholar] [CrossRef]
  18. Melchers, R.E.; Ahammed, M. A fast-approximate method for parameter sensitivity estimation in Monte Carlo structural reliability. Comput. Struct. 2004, 82, 55–61. [Google Scholar] [CrossRef]
  19. Puatatsananon, W.; Saouma, V.E. Reliability analysis in fracture mechanics using the first-order reliability method and Monte Carlo simulation. Fatigue Fract. Eng. Mater. Struct. 2010, 29, 959–975. [Google Scholar] [CrossRef]
  20. Pan, Q.; Dias, D. An efficient reliability method combining adaptive support vector machine and Monte Carlo simulation. Struct. Saf. 2017, 67, 85–95. [Google Scholar] [CrossRef]
  21. Fei, C.W.; Bai, G.C.; Tian, C. Extremum response surface method for casing radial deformation probabilistic analysis. AIAA J. Aerosp. Inf. Syst. 2013, 10, 47–52. [Google Scholar]
  22. Tvedt, L. Distribution of quadratic forms in normal space—Application to structural reliability. J. Eng. Mech. 1990, 116, 1183–1197. [Google Scholar] [CrossRef]
  23. Fei, C.W.; Bai, G.C. Distributed collaborative extremum response surface method for mechanical dynamic assembly reliability analysis. J. Cent. South Univ. 2013, 20, 2414–2422. [Google Scholar] [CrossRef]
  24. Kaymaz, I.; Marti, K. Reliability-based design optimization for elastoplastic mechanical structures. Comput. Struct. 2007, 85, 615–625. [Google Scholar] [CrossRef]
  25. Fei, C.W.; Tang, W.Z.; Bai, G.C.; Shuang, M. Dynamic probabilistic design for blade deformation with SVM-ERSM. Aircr. Eng. Aerosp. Technol. 2015, 87, 312–321. [Google Scholar] [CrossRef]
  26. Wei, Z.; Feng, F.; Wei, W. Non-linear partial least squares response surface method for structural reliability analysis. Reliab. Eng. Syst. Saf. 2017, 161, 69–77. [Google Scholar]
  27. Fei, C.W.; Bai, G.C. Distributed collaborative probabilistic design for turbine blade-tip radial running clearance using support vector machine of regression. Mech. Syst. Signal Process. 2014, 49, 196–208. [Google Scholar] [CrossRef]
  28. Fei, C.W.; Choy, Y.S.; Hu, D.Y.; Bai, G.C.; Tang, W.Z. Dynamic probabilistic design approach of high-pressure turbine blade-tip radial running clearance. Nonlinear Dyn. 2016, 86, 205–223. [Google Scholar] [CrossRef]
  29. Goswami, S.; Ghosh, S.; Chakraborty, S. Reliability analysis of structures by iterative improved response surface method. Struct. Saf. 2016, 60, 56–66. [Google Scholar] [CrossRef]
  30. Bai, G.C.; Fei, C.W. Distributed collaborative response surface method for mechanical dynamic assembly reliability design. Chin. J. Mech. Eng. 2013, 26, 1160–1168. [Google Scholar] [CrossRef]
  31. Hurtado, J.E.; Alvarez, D.A. An optimization method for learning statistical classifiers in structural reliability. Probab. Eng. Mech. 2010, 25, 26–34. [Google Scholar] [CrossRef]
  32. Zhang, C.Y.; Bai, G.C. Extremum response surface method of reliability analysis on two-link flexible robot manipulator. J. Cent. South Univ. 2012, 19, 101–107. [Google Scholar] [CrossRef]
  33. Lu, C.; Feng, Y.W.; Liem, R.P.; Fei, C.W. Improved kriging with extremum response surface method for structural dynamic reliability and sensitivity analyses. Aerosp. Sci. Technol. 2018, 76, 164–175. [Google Scholar] [CrossRef]
  34. Nose-Filho, K.; Lotufo, A.D.P.; Minussi, C.R. Short-term multinodal load forecasting using a modified general regression neural network. IEEE Trans. Power Deliv. 2011, 26, 2862–2869. [Google Scholar] [CrossRef]
  35. Zhao, C.; Liu, K.; Li, D.S. Freight volume forecast based on GRNN. J. China Railw. Soc. 2004, 26, 12–15. [Google Scholar]
  36. Li, H.Z.; Guo, S.; Li, C.J.; Sun, J.Q. A hybrid annual power load forecasting model based on generalized regression neural network with fruit fly optimization algorithm. Knowl. Based Syst. 2013, 37, 378–387. [Google Scholar] [CrossRef]
  37. Sun, G.; Hoff, S.J.; Zelle, B.C.; Smith, M.A. Development and comparison of backpropagation and generalized regression neural network models to predict diurnal and seasonal gas and PM 10 concentrations and emissions from swine buildings. Trans. Asabe 2008, 51, 685–694. [Google Scholar] [CrossRef]
  38. Wang, Y.; Peng, H. Underwater acoustic source localization using generalized regression neural network. J. Acoust. Soc. Am. 2018, 143, 2321–2331. [Google Scholar] [CrossRef]
  39. Gao, H.; Fei, C.; Bai, G.; Ding, L. Reliability-based low-cycle fatigue damage analysis for turbine blade with thermo-structural interaction. Aerosp. Sci. Technol. 2016, 49, 289–300. [Google Scholar] [CrossRef]
  40. Liu, C.L.; Lu, Z.Z.; Xu, Y.L.; Yue, Z.F. Reliability analysis for low cycle fatigue life of the aeronautical engine turbine disc structure under random environment. Mater. Sci. Eng. A 2005, 395, 218–225. [Google Scholar] [CrossRef]
  41. Vubac, N.; Lahmer, T.; Keitel, H.; Zhao, J.; Zhuang, X.; Rabczuk, T. Stochastic predictions of bulk properties of amorphous polyethylene based on molecular dynamics simulations. Mech. Mater. 2014, 68, 70–84. [Google Scholar] [CrossRef]
  42. Zhai, X.; Fei, C.W.; Zhai, Q.G.; Bai, G.C.; Ding, L. Reliability and sensitivity analyses of HPT blade-tip radial running clearance using multiply response surface model. J. Cent. South Univ. 2014, 21, 4368–4377. [Google Scholar] [CrossRef]
  43. Shu, M.A.; Wang, K.M.; Hui, M.; Shuai, Z. Research on turbine blade vibration characteristic under steady state temperature field. J. Shenyang Aerosp. Univ. 2011, 28, 18–21. [Google Scholar]
  44. Zhang, C.Y.; Lu, C.; Fei, C.W.; Liu, L.J. Multiobject reliability analysis of turbine blisk with multidiscipline under multiphysical field interaction. Adv. Mater. Sci. Eng. 2015, 2015, 519–520. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of generalized regression extreme neural network (GRENN) method.
Figure 1. Schematic diagram of generalized regression extreme neural network (GRENN) method.
Materials 12 01545 g001
Figure 2. Flow chart of reliability analysis with GRENN method.
Figure 2. Flow chart of reliability analysis with GRENN method.
Materials 12 01545 g002
Figure 3. FE model and gridding of a turbine blisk.
Figure 3. FE model and gridding of a turbine blisk.
Materials 12 01545 g003
Figure 4. Nephgrams of the responses of blisk stress and fatigue life.
Figure 4. Nephgrams of the responses of blisk stress and fatigue life.
Materials 12 01545 g004
Figure 5. Predicted results of the GRENN model with 30 groups of samples.
Figure 5. Predicted results of the GRENN model with 30 groups of samples.
Materials 12 01545 g005
Figure 6. Reliability analysis results of blisk fatigue life with the GRENN method.
Figure 6. Reliability analysis results of blisk fatigue life with the GRENN method.
Materials 12 01545 g006
Figure 7. Sensitivity degree distributions of random parameters on blisk LCF life.
Figure 7. Sensitivity degree distributions of random parameters on blisk LCF life.
Materials 12 01545 g007
Table 1. Distributions of random variables.
Table 1. Distributions of random variables.
Random VariablesMean μStandard Deviation δDistribution
Density ρ, kg.m−38210328.4Normal
Rotate speed ω, rad·s1116835Normal
Heat conductivity λ, W·m−1·°C−1230.005Normal
Modulus of elasticity, E, MPa1630004890Normal
Blade-root temperature Ta , k1173.1535.2Normal
Blade-tip temperature Tb, k1473.1547Normal
Fatigue strength efficient σ′f141942.5Normal
Fatigue ductility coefficient ε′f50.51.53Normal
Fatigue strength index b−0.10.005Normal
Fatigue ductility index c−0.840.042Normal
Table 2. Sensitivity degree and impact probability of the random input parameters.
Table 2. Sensitivity degree and impact probability of the random input parameters.
Random ParametersSensitivity Degree, ×10−3Effect Probability, %
ρ−0.415866.18
ω−0.525657.81
λ +0.01320.20
E+0.169482.52
T−1.7602226.16
σ′f+0.416156.18
ε′f+0.213113.17
b+0.435856.48
c+2.792941.30
Table 3. Computing time of the MC method, ERSM and GRENN.
Table 3. Computing time of the MC method, ERSM and GRENN.
Number of SamplesComputing Time under Different Simulations, sReduced Time, sImproved Efficiency, %
MC MethodERSMGRENN
10254001.2491.2010.0483.843
103144001.2661.2010.0655.134
1044320001.6811.3110.37015.18
1052.4371.3421.09544.93
1064.3122.1382.17450.42
Table 4. Computational precision of the reliability analysis methods under different simulations.
Table 4. Computational precision of the reliability analysis methods under different simulations.
SamplesReliability DegreePrecision/%Improved Precision/%
MC MethodERSMGRENNERSMGRENN
1020.850.760.7976.2479.253.01
1030.9760.9470.96895.0097.112.11
1040.99680.98240.997398.5699.951.39
1050.981810.9984898.4999.831.34
1060.982620.9958798.5899.911.33

Share and Cite

MDPI and ACS Style

Zhang, C.; Wei, J.; Jing, H.; Fei, C.; Tang, W. Reliability-Based Low Fatigue Life Analysis of Turbine Blisk with Generalized Regression Extreme Neural Network Method. Materials 2019, 12, 1545. https://doi.org/10.3390/ma12091545

AMA Style

Zhang C, Wei J, Jing H, Fei C, Tang W. Reliability-Based Low Fatigue Life Analysis of Turbine Blisk with Generalized Regression Extreme Neural Network Method. Materials. 2019; 12(9):1545. https://doi.org/10.3390/ma12091545

Chicago/Turabian Style

Zhang, Chunyi, Jingshan Wei, Huizhe Jing, Chengwei Fei, and Wenzhong Tang. 2019. "Reliability-Based Low Fatigue Life Analysis of Turbine Blisk with Generalized Regression Extreme Neural Network Method" Materials 12, no. 9: 1545. https://doi.org/10.3390/ma12091545

APA Style

Zhang, C., Wei, J., Jing, H., Fei, C., & Tang, W. (2019). Reliability-Based Low Fatigue Life Analysis of Turbine Blisk with Generalized Regression Extreme Neural Network Method. Materials, 12(9), 1545. https://doi.org/10.3390/ma12091545

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop