Next Article in Journal
Design of a Half-Bridge Current-Source Inverter Topology for Avionic Systems
Next Article in Special Issue
Speech GAU: A Single Head Attention for Mandarin Speech Recognition for Air Traffic Control
Previous Article in Journal
Incremental Backstepping Sliding-Mode Trajectory Control for Tailless Aircraft with Stability Enhancer
Previous Article in Special Issue
Semi-Physical Simulation of Fan Rotor Assembly Process Optimization for Unbalance Based on Reinforcement Learning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adaptive Local Maximum-Entropy Surrogate Model and Its Application to Turbine Disk Reliability Analysis

1
School of Energy and Power Engineering, Beihang University, Beijing 100191, China
2
Beijing Key Laboratory of Aero-Engine Structure and Strength, Beijing 100191, China
3
Collaborative Innovation Center of Advanced Aero-Engine, Beijing 100191, China
4
Aero Engine Academy of China, Aero Engine (Group) Corporation of China, Beijing 101304, China
5
School of Mechanical Engineering, Purdue University, West Lafayette, IN 47907, USA
*
Author to whom correspondence should be addressed.
Aerospace 2022, 9(7), 353; https://doi.org/10.3390/aerospace9070353
Submission received: 16 May 2022 / Revised: 28 June 2022 / Accepted: 29 June 2022 / Published: 30 June 2022

Abstract

:
The emerging Local Maximum-Entropy (LME) approximation, which combines the advantages of global and local approximations, has an unsolved issue wherein it cannot adaptively change the morphology of the basis function according to the local characteristics of the sample, which greatly limits its highly nonlinear approximation ability. In this research, a novel Adaptive Local Maximum-Entropy Surrogate Model (ALMESM) is proposed by constructing an algorithm that adaptively changes the LME basis function and introduces Particle Swarm Optimization to ensure the optimality of the adaptively changed basis function. The performance of the ALMESM is systematically investigated by comparison with the LME approximation, a Radial basis function, and the Kriging model in two explicit highly nonlinear mathematical functions. The results show that the ALMESM has the highest accuracy and stability of all the compared models. The ALMESM is further validated by a highly nonlinear engineering case, consisting of a turbine disk reliability analysis under geometrical uncertainty, and achieves a desirable result. Compared with the direct Monte Carlo method, the relative error of the ALMESM is less than 1%, which indicates that the ALMESM has considerable potential for highly nonlinear problems and structural reliability analysis.

1. Introduction

The performance of an aero-engine is determined by numerous factors, such as its geometric dimensions, material properties, and environments. In the process of the equipment design and manufacturing, these factors will inevitably experience certain random changes, which lead to security risks in the deterministic design [1,2,3,4]. The traditional method is to set a safety factor to consider these random uncertainties, which often results in an increased structural weight. This is particularly unfavorable for aero-engines that pursue high thrust-to-weight ratios. Therefore, a reliability analysis that considers the influence of various random factors came into being [5,6,7,8,9,10,11,12]. The Monte Carlo (MC) method [13] depends on numerous sample points to obtain statistical data. It can be performed efficiently when the performance function y(x) is expressed as an explicit form in terms of the random variables x. When the performance functions are implicit, the performance must be obtained through experiments or numerical simulations. Obviously, the MC method is impractical for numerous costly experiments, or time-consuming numerical simulations, such as 3-D Finite Element (FE) simulations [14,15], which are abundant in the design and manufacturing of aero-engines. Therefore, a surrogate model is proposed to replace these expensive experiments or simulations, which utilizes approximate methods to establish an approximate function that can represent the real function. It only needs a limited number of experiments or numerical simulations to establish a surrogate model. After that, The MC method is performed on the established surrogate model to greatly improve the calculation efficiency while meeting the accuracy. This approach has been used in many engineering applications, including predicting the Low Cycle Fatigue (LCF) life of turbine disks [16,17,18,19,20] and structure optimization [21,22,23].
Due to the extremely complex structure and its harsh operating environment, the performance of an aero-engine exhibits high nonlinearity in practice. The original and widely used Response Surface Method (RSM) [24,25,26,27,28] underperforms in dealing with this type of issue. Although it is possible to increase its nonlinear ability by increasing the polynomial order, overfitting will occur. A large amount of work was performed to improve the traditional RSM. Bai and Fei [29] utilized distributed collaborative RSM to design the reliability of the high-pressure turbine blade tip clearance of an aero-engine, and verified the feasibility and effectiveness of the method; Li and Lu [30] and Fei and Bai [31] combined the support vector machine (SVM) with RSM to analyze the reliability of structures; Yan et al. [32], Lu et al. [33] and Ren and Bai [34] combined the neural network (NN) algorithm with RSM to establish the neural network RSM, and used it for the regression analysis of the limit state surface of complex structures to improve the calculation accuracy. However, these methods still have relatively large limitations for highly nonlinear problems.
Generally, for some highly nonlinear problems in engineering, not only global optimal solutions but also local optimal solutions are required. Thus, the surrogate model needs to have the ability to balance the global and local approximations. A radial Basis Function (RBF) uses all the sample points in the design space to approximate the response function at a given point, and has strong global fitting capabilities [14,35,36]. It is inefficient to use sample points around prediction points. The Kriging method (KG) is a local approximation scheme, which largely depends on the sample points around the prediction point [37,38,39]. The global representation ability of KG decreases as the dimension increases. Arroyo and Ortiz [40] proposed a Local Maximum-Entropy (LME) approximation for the data interpolation of mesh-free methods, which combined the global maximum-entropy approximation method and the minimum width basis function based on Delaunay triangulation. Global or local capabilities can be adjusted by a dimensionless parameter γ . Based on this, Fan et al. [18] proposed an LME-based Surrogate Model (LMESM) considering both global and local approximations and found that the LME approximation has a great potential for dealing with highly nonlinear problems through a systematic comparison with RSM, RBF, and KG. In highly nonlinear problems, the degree of nonlinearity varies with space. The number of sample points used for fitting or interpolating different prediction points should be different. For example, in low-order nonlinear or linear regions, only a small number of sample points can achieve high accuracy, whereas highly nonlinear regions need more sample points to ensure their accuracy. Therefore, when constructing a surrogate model, the influence of a different number of sample points on the prediction points within different nonlinearity regions should be considered. In the original LME approximation and LMESM, the number of samples near the prediction points in different regions can be controlled by changing the basis function parameter γ . However, a globally unified γ was adopted in both, instead of updating γ according to the different number of sample points around the different prediction points, which greatly limits its highly nonlinear approximation ability. Moreover, γ was determined manually, which requires numerous trial and error approaches. A large amount of time is spent adjusting γ in the actual analysis. Even so, it is difficult to ensure the optimality of γ .
In this study, a new model for automatically updating basis functions, Adaptive Local Maximum-Entropy Surrogate Model (ALMESM), is proposed. A Particle Swarm Optimization (PSO) algorithm is adopted to optimize an objective function, which is the sum of the squared errors between the real and predicted values of all the sample points, to solve the optimal γ * in different regions automatically. The Centroidal Voronoi Tessellations (CVT) sampling method is adopted to ensure the uniformity of the sample point distribution. The performance of ALMESM is systematically studied by comparing the original LME approximation with a globally unified γ , RBF, and KG, with two explicit highly nonlinear mathematical functions. Finally, we apply this model to the reliability analysis of an aero-engine turbine disk with LCF life considering the uncertainty of the geometric parameters. The results show that ALMESM maintains higher accuracy based on the refinement of local characterization by comparison with the original LME approximation, RBF, and KG. Finally, the combination of ALMESM and the MC method is used to calculate the LCF life under different reliability and compared with the direct MC method. The relative error is within 1%, which proves that ALMESM can be used for the life reliability analysis of turbine disks.
The remainder of this study is organized as follows. In Section 2, the ALMESM approximation scheme is introduced. In Section 3, the performance of ALMESM is presented through two explicit highly nonlinear mathematical functions. The results from applying the proposed surrogate model to an industry application are presented in Section 4. Some conclusions are summarized in Section 5.

2. ALMESM

The original LME approximation method can be found in detail in Arroyo and Ortiz [40]. For the sake of the completeness of this study, here is a brief review of the process of establishing this method. In the original LME approximation, an implicit performance function y(x), at a given prediction point x, is approximated by:
y ^ ( x ) = i = 1 N y ( x ( i ) ) p i ( x )
where x ( i ) is the n dimensional vector x ( i ) = ( x 1 ( i ) , x 2 ( i ) , , x n ( i ) ) T , x ( i ) X . X = { x ( i ) , i = 1 , 2 , , N } is the set of N sample points, and y ( x ( i ) ) is the corresponding response at x ( i ) . p i ( x ) is the basis function of the sample point x ( i ) at a given prediction point x.
Let Ω M a x E n t be the convex hull of the sample point set X, and the basis function of the sample point x ( i ) at a given prediction point x is defined as the solution of problem (2):
( Max - Ent ) β Minimize f β ( x , p ) β i = 1 N p i | x x ( i ) | 2 + i = 1 N p i log p i
given x Ω M a x E n t ,
p i ( x ) 0 , i [ 1 , N ] p i ( x ) = 1 p i ( x ) x ( i ) = x
where β R + is a Pareto optimal parameter used to ponder the locality term against the entropy term in the objective function.
The constraints in the minimization problem ensure that the resulting basis functions satisfy the zeroth and first order consistencies. Consequently, this model can exactly reproduce constant and linear functions. The solution of the minimization problem is unique; therefore, the basis function can be directly calculated as:
p i ( x ) = 1 Z ( x , λ * ) exp [ β | x x ( i ) | 2 + λ * ( x x ( i ) ) ]
where
λ * = arg min λ n log Z ( x , λ )
Z ( x , λ ) = i = 1 N exp [ β | x x ( i ) | 2 + λ ( x x ( i ) ) ]
The calculation of the basis function in (3) requires finding β and λ * . β was determined by a dimensionless parameter γ and the local characteristic length h of the prediction point x :
γ = β h 2 , γ [ 0 , + )
where the characteristic length h of a given prediction point x is calculated by finding the nearest sample points x ( i ) around x , and then finding the nearest sample point x ( j ) to x ( i ) . Thus, the Euclidean distance between x ( i ) and x ( j ) is defined as h of x . In the previous report [18], γ usually ranges from 0.1 to 6.8. It determines the support width of the LME basis functions, which introduces a seamless bridge between local approximation and global approximation as shown in Figure 1. Obviously, the LME basis functions are allowed to depend on the location of the prediction point and have the potential to adjust themselves adaptively by varying γ to achieve optimal degrees of locality, which provides a great advantage for the construction of high-order approximations. The basis function itself covers all global sample points, but the locality of the basis function determines the contribution of each sample point to the predicted point. By setting a truncation error p t o l (e.g., p t o l = 10 5 ), the sample points with minimal contributions can be ignored. Therefore, by controlling the locality of the basis function, the global and local capabilities of the surrogate model can be adjusted. As shown in Figure 1, as γ increases, the basis function support width gradually narrows, and with fewer sample points with a significant contribution included, the local ability will be stronger. On the contrary, the global capability will be strengthened.
Nevertheless, in the previous research work [18,40], the value of γ was manually set to a globally uniform value, which is difficult for determining an optimal value and requires a lot of trial and experience. In this study, by introducing an objective function for optimization, the optimal γ i * of each sample point is obtained.
Optimal γ i * should minimize the fitting error of the surrogate model. That is, the predicted value at the sample point should be as close as possible to the real value. In this study, when calculating the optimal γ * of a sample point, the surrounding sample points with a significant contribution to this sample point are selected as a local neighborhood X l o c a l . Then, we take the sum of the squared errors of the real and predicted values of the sample points in X l o c a l as the objective function for optimization, i.e.,
γ i * = arg min i = 1 k ( y ( x ( i ) ) y ^ ( x ( i ) ) ) 2 , γ i * [ 0 , + )
where y ( x ( i ) ) is the real response value of the sample points, y ^ ( x ( i ) ) is the predicted value of the sample points, and k is the number of sample points in the local set X l o c a l of the prediction sample point by selecting p i ( x ( i ) ) > p t o l .
It can be observed that γ i * and λ * can be obtained by iteratively solving (4) and (7). In the original LME approximation, λ * was optimized by a Newton–Raphson algorithm and a Nelder–Mead simplex algorithm. However, the Newton–Raphson algorithm based on the gradient is usually limited by the initial value and derivative operation, resulting in convergence to the local minimum solution or divergence. The emergence of an intelligent optimization algorithm effectively avoids the abovementioned problems. The PSO algorithm, proposed by Kennedy J. and Eberhart R. [41], has the advantages of a fast convergence and strong global optimization ability [42]. Therefore, the PSO algorithm is adopted in this study to solve (4) and (7) to obtain γ i * and λ * .
Suppose there is a set { x 1 ,   x 2 ,   ,   x m } of m particles in a D-dimensional objective search space. The i-th particle is represented as a D-dimensional solution vector x i = ( x i 1 , x i 2 , , x i D ) . It describes the position of the particle in the search space. The velocity of the individual particles is recorded as V i = ( v i 1 , v i 2 , , v i D ) . The best position that the individual particle has experienced, i.e., the individual extreme value, is recorded as P i = ( p i 1 , p i 2 , , p i D ) . The best position that the entire group has experienced, i.e., the extreme value of the group, is recorded as P g = ( p g 1 , p g 2 , , p g D ) . The PSO algorithm updates the particle velocity and position in the following manner based on P i and P g :
v i = ω v i - 1 + c 1 r 1 ( p i x i ) + c 2 r 2 ( p g x i )
x i + 1 = x i + α v i ,   i = 1 ,   2 ,   , m
where ω 0 is the inertial factor, c 1 0 and c 2 0 are the acceleration factors, r 1 and r 2 are the random number between [0, 1], and α is the constraint factor which is the weight that controls the velocity. The convergence condition adopted in this study is that the values of the objective function (4) and (7) remain unchanged for five consecutive iterations. Finally, the ALMESM algorithm resulting from the preceding scheme is listed in Table 1.
Through the ALMESM algorithm, an optimized γ i * for the location can be automatically solved according to the nonlinear degree of the location of each sample point, which is a better fit to the real model. Figure 2 shows the characteristics of the ALMESM algorithm. The gray solid line represents the original function, and the dots on the line represent the sample points. x a and x b are prediction points, which are in two regions with different nonlinearities. X l o c a l ( x a ) and X l o c a l ( x b ) are the neighborhoods of the two prediction points. It can be observed from Figure 2 that these two neighborhoods contain different numbers of sample points, and the value of γ i * at each sample point is different, leading to a different basis function.

3. Performance

To verify that the ALMESM can obtain the optimal γ * of each sample point, as well as the improvement of the performance of the original LME approximation, we compare it with the original LME approximation, an RBF surrogate model with a Gaussian basis function, and a KG surrogate model with a Gaussian correlation function. The detailed introduction of RBF and KG can be found in Simpson et al. [43], Wild et al. [44], and Forrester and Keane [45]. Two explicit mathematical functions are selected as the test cases from Hock and Schittkowski [46] for testing the highly nonlinear optimization algorithms.
(a)
A 3-D function with high-order nonlinearity
y ( x ) = 10 + i = 1 n ( x i 2 + 5 cos ( 2 π x i ) ) n = 2 ,   1.0 x i 1.0
(b)
A 6-D function with high-order nonlinearity
y ( x ) = i = 1 n ( ( ln ( x i 2 ) ) 2 + ( ln ( 10 x i ) ) 2 ) i = 1 n x i 2 n = 5 , 2.1 x i 9.9
R-squared (R2) is a widely used assessment metric for surrogate models, which is used in the study to make a quantitative assessment of the fitting performance of a surrogate model
R 2 = 1 i = 1 n ( y ( x i ) y ^ ( x i ) ) 2 i = 1 n ( y ( x i ) y ¯ ( x i ) ) 2 = 1 MSE v a r i a n c e
where y ^ ( x i ) is the prediction value obtained by the surrogate model, and y ¯ ( x i ) is the mean of the response values of all the sample points. While Mean Square Error (MSE) represents the departure of the surrogate model from the real function, the variance captures how irregular the problem is. R 2 embodies the fitting degree between the surrogate model and the real model. The larger the value is, the higher the fitting accuracy of the surrogate model will be.
When constructing a surrogate model, it is important to select the sample points scientifically for the accuracy of the surrogate model. Generally, it is required that the selected sample points can evenly cover the entire design space, and the number of sample points should be as few as possible. Considering these two requirements, this article uses the Centroidal Voronoi Tessellations Sampling (CVT) method for sampling. Let { x i } ( i = 1 , 2 , , n ) be n points in M-dimensional hyper cubic space. The Voronoi cell V i corresponding to x i is defined as
V i = { x Ω | | x x i | < | x x j | , j = 1 , 2 , , n , j i }
where | x x i | represents the Euclidian distance between x and x i ; { V i } is called the Voronoi diagram. For each Voronoi unit, the center of gravity x i * is defined as
x i * = V i x d x V i d x
If and only if a Voronoi diagram meets the following conditions (15), it is called the gravity center Voronoi diagram.
x i = x i *
That is, in the gravity center Voronoi diagram, the generating point of each Voronoi unit happens to be the gravity center of the unit. Currently, { x i } ( i = 1 , 2 , , n ) forms a set of CVT points. Figure 3 shows a set of CVT points in a two-dimensional space with 100 sample points and 50 prediction points. It can be observed that the sample points and prediction points evenly cover the entire design space.
We first explored the influence of the randomness of the sampling on the results of the original LME approximation and ALMESM. Assuming that the sampling scheme is represented by (A, B), where A is the number of sample points and B is the number of prediction points, we selected four sampling schemes, i.e., (60,30), (80,40), (100,50), (120,60), and conducted five CVT random samplings for each scheme. Taking the 3-D function (10) as the example, we compared the R 2 of the original LME approximation with the globally unified γ = 0.8 (randomly selected) and ALMESM. It can be observed from Figure 4 that the results of the two algorithms are less dispersive, and as the number of sample points increases, the dispersion tends to further decrease. Moreover, we performed statistics on the dispersion results (see Table 2). R max 2 and R min 2 represents the maximum and minimum value of R 2 in the five CVT samplings, respectively. It can be found that the maximum absolute error in R 2 of the original LME approximation is 0.08, resulting in a relative error of 8.70%. As a contrast, ALMESM gives the absolute error of 0.04 and the relative error of 4.21%. The above results show that the results of the two algorithms are less affected by the randomness of sampling. Moreover, the stability (the degree of variation in the prediction accuracy of a surrogate model that is affected by the number of sampling points and the randomness of sampling) of ALMESM relative to the original LME approximation is further improved.
According to the above mathematical functions, eight sets of sample points with the number in [60, 200] were extracted by the CVT method. Table 3 is the selection scheme of the sample points and prediction points. We compared the R 2 trend graph of ALMESM, the original LME approximation with a globally unified γ , RBF, and KG.
In the 3-D function (Figure 5), the accuracy of the original LME approximation changes significantly with γ when the sample size is small. As the sample size increases, the accuracy gradually stabilizes. However, when the dimension increases (Figure 6), the accuracy of the original LME approximation is greatly influenced by γ and the number of sample points. The accuracy varies widely, and even negative values appear. This is because when the number of sample points increases, the global and local nonlinearity of the original function will be more prominent, making the basis function that was originally applicable to a smaller number of sample points not applicable at this time, resulting in reduced accuracy. The original LME approximation cannot automatically adjust γ to change the locality of the basis function. Thus, it leads to low accuracy. In contrast, it can be observed from the two examples that ALMESM always maintains high accuracy. ALMESM can automatically adjust γ so that the basis function can be adaptively adjusted in various situations. Therefore, under different sample points, it can be stabilized with high accuracy. Figure 7 shows the optimal γ * of each sample point obtained by ALMESM in the two examples. It can be observed that the increase in the dimensionality and the degree of nonlinearity causes the γ of each sample point to vary more drastically than the 3-D ones. The reason is that the nonlinearity of the 6-D function is more severe than that of the 3-D function, and the basis functions with higher γ are required to capture local features. Therefore, the large amount of the optimal γ * obtained by ALMESM for the 6-D function are quite large (>10.0) and it is difficult to achieve higher accuracy with a globally unified γ . In summary, compared to the original LME approximation, ALMESM improves the accuracy and stability of the calculations.
In Figure 8, the accuracy of KG is undesirable when the sample size is small, and high accuracy is not obtained until the number of sample points reaches 140. Whereas RBF and ALMESM can reach high accuracy in the range of 60–200 sample points. With the increase in the dimension and nonlinear order (Figure 9), the fitting accuracy of the three models gradually increases as the number of samples increases, but ALMESM has the highest accuracy. It can be observed from the above two examples that ALMESM realizes the automatic optimization of γ for each sample point and can achieve higher accuracy dealing with highly nonlinear problems. In the following, ALMESM is used for the LCF and reliability analysis of an aero-engine turbine disk considering geometric uncertainty.

4. Application

As a life-limiting part of an aero-engine, the reliability of the turbine disk is a weak link in the aero-engine design and has become a bottleneck restricting engine research and development [47]. The LCF is the most important failure mode of turbine disks. In practice, the structure of the turbine disk is uncertain due to the dispersion of the dimensional tolerances, so that the actual life of the turbine disk is statistically dispersed. Meanwhile, the operation environment of the turbine disk is complex. There are multiple load conditions, which leads to a strong nonlinearity in the reliability prediction of the turbine disk.
Structural reliability is the probability that a structure completes a specified function under specified conditions within a specified time t . Assuming that the performance function of the structure is y ( x ) = y ( x 1 , x 2 , , x n ) , the limit state equation y ( x 1 , x 2 , , x n ) = 0 divides the basic variable space of the structure into a failure domain and a reliability domain. The failure probability P f can be obtained by integrating the joint probability density function of the basic random variable x . For general multi-dimensional problems and complex integration domain or implicit integration domain problems, the integration of the failure probability has no analytical solution. In such a case, the MC method can be used to solve this problem. Based on the MC method, the reliability function R ( t ) can be achieved as:
R ( t ) = 1 P f = 1 n f ( t ) n 0 = n 0 n f ( t ) n 0
where n 0 represents the number of total sample points, and n f ( t ) represents the cumulative number of sample points in the failure domain during the working time from 0 to t .
In this study, the ALMESM was used to predict the LCF life and perform the reliability analysis of a turbine disk considering geometric uncertainty, which has been divided into the following four steps (also as outlined in Figure 10):
(1)
Filter key geometric variables. As shown in Figure 11, there are numerous geometric parameters of the turbine disk, and it is necessary to perform a sensitivity analysis on all the geometric parameters to screen out the geometric variables that have the greatest influence on the LCF. Fan et al. [17] have detailed the process of sensitivity analysis, pointing out that the inner diameter (d_R1), the outer diameter (d_R2), and the rim thickness (d_W6) are the three key variables for LCF analysis.
(2)
Determine the probability distributions of the key variables. For the reliability analysis, the fidelity of the analysis strongly depends on how well we know the underlying distributions of the random variables. In TDRA, the Kolmogorov–Smirnov (K-S) test approach [48,49] is adopted to determine the probability distributions of the three key geometric variables. Fan et al. [17] have found that the three geometric variables obey the normal distribution.
(3)
Construct the ALMESM. For fast and high-fidelity MC simulations, the ALMESM needs to be constructed. First, perform CVT sampling on the three key geometric variables. Then, carry on the FE analysis to the sample points, obtain the LCF life of each sample point according to the LCF life formula, and finally train out ALMESM.
(4)
Predict probability for fatigue life. Perform the MC simulation on the constructed ALMESM to obtain the reliability analysis results.
Figure 10. Flow chart for the reliability analysis of a turbine disk considering geometric uncertainties in the Turbine Disk Reliability Analysis (TDRA).
Figure 10. Flow chart for the reliability analysis of a turbine disk considering geometric uncertainties in the Turbine Disk Reliability Analysis (TDRA).
Aerospace 09 00353 g010
Figure 11. Variable geometric parameter of the GH720Li turbine disk, where the three key parameters are marked with blue circles.
Figure 11. Variable geometric parameter of the GH720Li turbine disk, where the three key parameters are marked with blue circles.
Aerospace 09 00353 g011
The material properties of GH720Li at 650 °C and the boundary conditions applied to the disk in the FE simulations of step (3) are listed in Table 4. According to the LCF test data of GH720Li and stress–strain model and strain–life model in Tang and Lu [50] and Gao et al. [51], respectively, the Mason–Coffin model of GH720Li material at 650 °C was obtained, as shown in (17):
Δ σ 2 × 10 5.2671 0.0003 μ 1 + 10 0.2102 μ 1 ( Δ σ 2 × 3.4326 × 10 4 ) 6.2723 = ( 10 24.1541 + 3.7741 μ 2 12.0456 + 1.7532 μ 2 σ m E ) ( 2 N f ) 1 12.0456 + 1.7532 μ 2 + ( 10 2.4733 0.0834 μ 2 1.9361 + 0.0494 μ 2 ) ( 2 N f ) 1 1.9361 + 0.0494 μ 2
where N f is the LCF life, Δ σ is the stress range, σ m is the mean stress, E is the elastic modulus, and μ 1 and μ 2 are two standard normal random variables.
A total of 10 sets of sample points (Table 5) were selected to test the accuracy of the surrogate model by using CVT sampling. To verify the performance of the model proposed in this article, we compared it with the original LME approximation with a globally unified γ , RBF, and KG. From Figure 12, the accuracy of the original LME approximation is greatly influenced by γ and the number of sample points. An improper value of γ can even lead to a negative value of R 2 ( γ = 6.8, 120 sample points). In contrast, under different sample points, the accuracy of ALMESM can always be maintained at a high level, and even the lowest accuracy is higher than 0.80. Overall, the accuracy tends to increase as the number of sample points increases. When the number of sample points is 240, the accuracy of ALMESM reaches 0.96, and gradually tends to converge. In summary, compared to the original LME approximation, ALMESM realizes the automatic optimization of γ for each sample point, and improves the accuracy and stability of calculations. From Figure 13, KG has a certain prediction accuracy when the sample is small. However, as the number of sample points increases (>80), the predictive ability of KG is lost in this highly nonlinear engineering case. This phenomenon was also reported in literature [18]. The KG model is a local approximation scheme; the correlation matrix can become singular if multiple sample points are spaced close to one another. Under a different number of sample points, the accuracy of ALMESM and RBF can be maintained at a high level. However, the ALMESM has higher accuracy overall. When the number of sample points is 100, the difference of R 2 between ALMESM and RBF is the largest, with a value of 0.11. In summary, ALMESM can achieve higher accuracy for approximating the LCF life of the GH720Li turbine disk considering geometric uncertainty.
Subsequently, the LCF life reliability of the turbine disk was calculated using ALMESM, and the LCF life value under different reliability was obtained. Since the accuracy of ALMESM gradually converges when the number of samples is 240 (Figure 13), the surrogate model established at this time was used for TDRA. A total of 100,000 sample points were created by using the MC random sampling method. By conducting MC simulations 100,000 times on ALMESM the frequency histogram of the LCF life was obtained, as shown in Figure 14. It can be observed that the LCF life approximately obeys the lognormal distribution within the scope of the sampling. The minimum LCF life is 148,392 cycles and the maximum LCF life is 154,896 cycles. When the LCF life equals 150,892 cycles, it reaches the highest frequency of 0.057. In addition, the reliability–life curve was obtained as shown in Figure 15. It is evident that although the variation of the dimensional tolerance is small, it still has a significant contribution to the LCF life of the turbine disk.
To verify the reliability results of the ALMESM, direct MC simulations (FE + MC) were performed on the turbine disk model based on the previous 100,000 samples. Through a statistical analysis, Table 6 displays the results of the same samples to the LCF life of a certain reliability level obtained by the two methods and the total computational cost of the analysis. The results obtained by the two methods are relatively close to each other. The relative error is within 1%. The ALMESM is more efficient than the direct MC simulation. The comparison results fully demonstrate that it is feasible to use ALMESM for a reliability analysis.

5. Conclusions

In this study, the ALMESM, a novel surrogate model that can adaptively adjust and calculate the optimal basis function according to the degree of nonlinearity, is proposed. The performance of the ALMESM is systematically investigated by comparison to the original LME approximation, RBF, and KG model in two explicit highly nonlinear mathematical functions. The performance of the ALMESM is further demonstrated via the application of a highly nonlinear turbine disk reliability analysis. Compared with the original LME approximation, RBF, and KG model, the model in all the test cases maintained high accuracy and robustness. Moreover, there is no need to manually adjust the parameters according to the number of sample points. These attractive properties are derived from the calculation of the basis functions. The basis functions are allowed to depend on the location of the prediction point and can adjust themselves adaptively by varying γ to achieve optimal degrees of locality. These results and capabilities demonstrate that ALMESM is a feasible surrogate model for highly nonlinear problems.

Author Contributions

Conceptualization, J.F.; methodology, Q.Y.; validation, H.X.; formal analysis, Q.M.; investigation, H.X.; resources, F.J.; writing—original draft preparation, Q.Y.; writing—review and editing, J.F., Q.Y. and H.W.; visualization, Q.Y.; supervision, J.F.; project administration, F.J. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the National Natural Science Foundation of China (Grant Nos. 51875023 and 51775015).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data generated or analyzed during this study are included in this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Keshtegar, B. A hybrid conjugate finite-step length method for robust and efficient reliability analysis. Appl. Math. Model. 2017, 45, 226–237. [Google Scholar] [CrossRef]
  2. Song, L.K.; Wen, J.; Fei, C.W.; Bai, G.C. Distributed collaborative probabilistic design of multi-failure structure with fluid-structure interaction using fuzzy neural network of regression. Mech. Syst. Signal. Pr. 2018, 104, 72–86. [Google Scholar] [CrossRef]
  3. 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]
  4. Lu, C.; Feng, Y.W.; Fei, C.W.; Bu, S.Q. Decomposed-coordinated framework with enhanced extremum kriging for multicomponent dynamic probabilistic failure analyses. IEEE Access 2019, 7, 163287–163300. [Google Scholar] [CrossRef]
  5. Grigoriu, M. Stochastic mechanics. Int. J. Solids Struct. 2000, 37, 197–214. [Google Scholar] [CrossRef]
  6. Rahman, S.; Rao, B.N. An element-free Galerkin method for probabilistic mechanics and reliability. Int. J. Solids Struct. 2001, 38, 9313–9330. [Google Scholar] [CrossRef]
  7. Ding, Z.H.; Zhou, Y.; Pu, G.G.; Zhou, M.C. Online failure prediction for railway transportation systems based on fuzzy rules and data analysis. IEEE T. Reliab. 2018, 67, 1143–1158. [Google Scholar] [CrossRef]
  8. Duan, Y.; Li, W.F.; Fu, X.W.; Luo, Y.; Yang, L. A methodology for reliability of WSN based on software defined network in adaptive industrial environment. IEEE-CAA J. Autom. 2018, 5, 74–82. [Google Scholar] [CrossRef]
  9. Zhu, S.P.; Liu, Q.; Zhou, J.; Yu, Z.Y. Fatigue reliability assessment of turbine discs under multi-source uncertainties. Fatigue Fract. Eng. M 2018, 41, 1291–1305. [Google Scholar] [CrossRef]
  10. Jiang, C.; Qiu, H.B.; Yang, Z.; Chen, L.M.; Gao, L.; Li, P.G. A general failure-pursuing sampling framework for surrogate-based reliability analysis. Reliab. Eng. Syst. Saf. 2019, 183, 47–59. [Google Scholar] [CrossRef]
  11. Song, L.K.; Bai, G.C.; Fei, C.W. Probabilistic LCF life assessment for turbine discs with DC strategy-based wavelet neural network regression. Int. J. Fatigue 2019, 119, 204–219. [Google Scholar] [CrossRef]
  12. Wei, Y.X.; Bai, G.C.; Song, L.K.; Bai, B. The estimation of reliability probability of structures based on improved iterative response surface methods. Ksce J. Civ. Eng. 2019, 23, 4063–4074. [Google Scholar] [CrossRef]
  13. Bjerager, P. Probability integration by directional simulation. J. Eng. Mech. 1988, 114, 1285–1302. [Google Scholar] [CrossRef]
  14. Deng, J. Structural reliability analysis for implicit performance function using radial basis function network. Int. J. Solids Struct. 2006, 43, 3255–3291. [Google Scholar] [CrossRef] [Green Version]
  15. Wei, Y.X.; Bai, G.C.; Song, L.K. A novel reliability analysis approach with collaborative active learning strategy-based augmented RBF metamodel. IEEE Access 2020, 8, 199603–199617. [Google Scholar] [CrossRef]
  16. Fan, J.; Wang, H. Probabilistic assessment of turbine disk considering geometry distribution based on surrogate models. In Proceedings of the ASME 2014 International Mechanical Engineering Congress and Exposition, Montreal, QC, Canada, 14–20 November 2014; ASME: New York, NY, USA, 2014. [Google Scholar]
  17. Fan, J.; Liao, H.; Li, D.; Wang, R.Q.; Hu, D.Y. Probabilistic analysis of turbine disk fatigue life considering geometric uncertainties. J. Aerosp. Eng. 2017, 32, 66–74. (In Chinese) [Google Scholar]
  18. Fan, J.; Liao, H.M.; Wang, H.; Hu, J.H.; Chen, Z.Y.; Lu, J.; Li, B. Local maximum-entropy based surrogate model and its application to structural reliability analysis. Struct. Multidiscip. Optim. 2018, 57, 373–392. [Google Scholar] [CrossRef]
  19. Fan, J.; Chen, G.X.; Shen, X.L.; Fan, T.Y. Evaluation method of uncertainty of reliability calculation for turbine disk life. In Proceedings of the ASME Turbo Expo: Turbomachinery Technical Conference and Exposition, Oslo, Norway, 11–15 June 2018; ASME: New York, NY, USA, 2018. [Google Scholar]
  20. Chen, G.X.; Fan, J.; Xu, H.B.; Li, B. Calculation of hybrid reliability of turbine disk based on self-evolutionary game model with few shot learning. Struct. Multidiscip. Optim. 2021, 63, 807–819. [Google Scholar] [CrossRef]
  21. Sun, G.Y.; Song, X.G.; Baek, S.; Li, Q. Robust optimization of foam-filled thin-walled structure based on sequential Kriging metamodel. Struct. Multidiscip. Optim. 2014, 49, 897–913. [Google Scholar] [CrossRef]
  22. Gu, X.G.; Lu, J.W.; Wang, H.Z. Reliability-based design optimization for vehicle occupant protection system based on ensemble of metamodels. Struct. Multidiscip. Optim. 2015, 51, 533–546. [Google Scholar] [CrossRef]
  23. Shang, X.B.; Ma, P.; Yang, M.; Chao, T. An efficient polynomial chaos-enhanced radial basis function approach for reliability-based design optimization. Struct. Multidiscip. Optim. 2021, 63, 789–805. [Google Scholar] [CrossRef]
  24. Guan, X.L.; Melchers, R.E. Effect of response surface parameter variation on structure reliability estimates. Struct. Saf. 2001, 23, 429–444. [Google Scholar] [CrossRef]
  25. Gavin, H.P.; Yau, S.C. High-order limit state functions in the response surface method for structural reliability analysis. Struct. Saf. 2008, 30, 162–179. [Google Scholar] [CrossRef]
  26. Sun, S.; Liu, S.; Jia, D.; Zhang, H.; Liu, F. Multiple nonlinear regression model weld bead shape high nitrogen steal by laser-arc hybrid welding. J. Mech. Eng. 2015, 51, 67–75. [Google Scholar] [CrossRef]
  27. Teng, F.; Zhang, W.X.; Liang, J.C.; Gao, S. Springback prediction and optimization of variable stretch force trajectory in three-dimensional stretch bending process. Chin. J. Mech. Eng. 2015, 28, 1132–1140. [Google Scholar] [CrossRef]
  28. Zhang, J.; Guo, F. Statistical modification analysis of helical planetary gears based on response surface method and Monte Carlo simulation. Chin. J. Mech. Eng. 2015, 28, 1194–1203. [Google Scholar] [CrossRef]
  29. 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]
  30. Li, H.; Lu, Z. A support vector machine response surface method for structural reliability analysis. Chin. J. Comput. Mech. 2009, 2, 199–203. (In Chinese) [Google Scholar]
  31. Fei, C.; Bai, G. Nonlinear dynamic probabilistic analysis for turbine casing radical deformation using extremum response surface method based on support vector machine. J. Comput. Nonlin. Dyn. 2013, 8, 041004. [Google Scholar] [CrossRef]
  32. Yan, H.; Hu, Y.; Niu, Y. Approximate structure reliability analysis based on NN-RSM. Ocean. Eng. 2002, 20, 1–6. [Google Scholar]
  33. Lu, Z.; Yang, Z.; Zhao, J. An artificial neural network method for reliability analysis based on weighted linear response surface. Chin. J. Aeronaut. 2006, 6, 1063–1067. (In Chinese) [Google Scholar]
  34. Ren, Y.; Bai, G. New neural network response surface methods for reliability analysis. Chin. J. Aeronaut. 2011, 24, 25–31. [Google Scholar] [CrossRef] [Green Version]
  35. Li, X.; Gong, C.L.; Gu, L.X.; Gao, W.K.; Jing, Z.; Su, H. A sequential surrogate method for reliability analysis based on radial basis function. Struct. Saf. 2018, 73, 42–53. [Google Scholar] [CrossRef] [Green Version]
  36. Wang, Q.; Fang, H.B. Reliability analysis of tunnels using an adaptive RBF and a first-order reliability method. Comput. Geotech. 2018, 98, 144–152. [Google Scholar] [CrossRef]
  37. Panda, S.S.; Manohar, C.S. Applications of meta-models in finite element based reliability analysis of engineering structures. Comput. Model Eng. Sci. 2008, 28, 161–184. [Google Scholar]
  38. 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]
  39. Song, L.K.; Bai, G.C.; Fei, C.W. Dynamic surrogate modeling approach for probabilistic creep-fatigue life evaluation of turbine disks. Aerosp. Sci. Technol. 2019, 95, 105439. [Google Scholar] [CrossRef]
  40. Arroyo, M.; Ortiz, M. Local maximum-entropy approximation schemes: A seamless bridge between finite elements and meshfree methods. Int. J. Numer. Methods Eng. 2006, 65, 2167–2202. [Google Scholar] [CrossRef] [Green Version]
  41. Kennedy, J.; Eberhart, R. Particle Swarm Optimization. In Proceedings of the 1995 IEEE International Conference on Neural Networks, Perth, Australia, 27 November–1 December 1995; IEEE: New York, NY, USA, 1995. [Google Scholar]
  42. Parsopoulos, K.E.; Vrahatis, M.N. Recent approaches to global optimization problems through particle swarm optimization. Nat. Comput. 2002, 1, 235–306. [Google Scholar] [CrossRef]
  43. Simpson, T.W.; Lin, D.K.; Chen, W. Sampling strategies for computer experiments: Design and analysis. Int. J. Reliab. Appl. 2001, 2, 209–240. [Google Scholar]
  44. Wild, S.M.; Regis, R.G.; Shoemaker, C.A. Orbit: Optimization by radial basis function interpolation in trust-regions. SIAM J. Sci. Comput. 2008, 30, 3197–3219. [Google Scholar] [CrossRef]
  45. Forrester, A.I.J.; Keane, A.J. Recent advances in surrogate-based optimization. Prog. Aerosp. Sci. 2009, 45, 50–79. [Google Scholar] [CrossRef]
  46. Hock, W.; Schittkowski, K. Test examples for nonlinear programming codes. J. Optim. Theory Appl. 1980, 30, 127–129. [Google Scholar] [CrossRef]
  47. Wang, R.; Hu, D.; Shen, X.; Fan, J. Probabilistic design technology on typical structure of aeroengine. Aeronaut. Manuf. Technol. 2014, 7, 26–30. [Google Scholar]
  48. Birnbaum, Z.W. Numerical tabulation of the distribution of Kolmogorov’s statistic for finite sample size. J. Am. Stat. Assoc. 1952, 47, 425–441. [Google Scholar] [CrossRef]
  49. Lilliefors, H.W. On the Kolmogorov-Smirnov test for normality with mean and variance unknown. J. Am. Stat. Assoc. 1967, 62, 399–402. [Google Scholar] [CrossRef]
  50. Tang, J.; Lu, S. Reliability analysis method for the strain fatigue life of a disk. J. Propul. Technol. 2005, 26, 344–347. (In Chinese) [Google Scholar]
  51. Gao, Y.; Yang, C.J.; Lin, K.; Gao, Q. Reliability analysis of LCF life for a turbine disk. Adv. Mater. Res. 2011, 146–147, 1379–1385. [Google Scholar] [CrossRef]
Figure 1. The basis functions with different γ in a two−dimensional convex domain.
Figure 1. The basis functions with different γ in a two−dimensional convex domain.
Aerospace 09 00353 g001
Figure 2. Schematic diagram of the ALMESM basis function in a case with one variable.
Figure 2. Schematic diagram of the ALMESM basis function in a case with one variable.
Aerospace 09 00353 g002
Figure 3. The CVT sampling scheme of a two−dimensional case.
Figure 3. The CVT sampling scheme of a two−dimensional case.
Aerospace 09 00353 g003
Figure 4. Dispersion analysis of 5 CVT samplings of the original LME approximation and ALMESM. (a) 60 sample points; (b) 80 sample points; (c) 100 sample points; (d) 120 sample points.
Figure 4. Dispersion analysis of 5 CVT samplings of the original LME approximation and ALMESM. (a) 60 sample points; (b) 80 sample points; (c) 100 sample points; (d) 120 sample points.
Aerospace 09 00353 g004
Figure 5. R-squared comparison of ALMESM, and the original LME approximation with different globally unified γ under different number of sample points in the 3-D function.
Figure 5. R-squared comparison of ALMESM, and the original LME approximation with different globally unified γ under different number of sample points in the 3-D function.
Aerospace 09 00353 g005
Figure 6. R-squared comparison of ALMESM, and the original LME approximation with different globally unified γ under different number of sample points in the 6-D function.
Figure 6. R-squared comparison of ALMESM, and the original LME approximation with different globally unified γ under different number of sample points in the 6-D function.
Aerospace 09 00353 g006
Figure 7. The γ * at different sample points of the two test functions in the case of ALMESM with 100 sample points.
Figure 7. The γ * at different sample points of the two test functions in the case of ALMESM with 100 sample points.
Aerospace 09 00353 g007
Figure 8. R-squared comparison of ALMESM, RBF, and KG under different number of sample points in the 3-D function.
Figure 8. R-squared comparison of ALMESM, RBF, and KG under different number of sample points in the 3-D function.
Aerospace 09 00353 g008
Figure 9. R-squared comparison of ALMESM, RBF, and KG under different number of sample points in the 6-D function.
Figure 9. R-squared comparison of ALMESM, RBF, and KG under different number of sample points in the 6-D function.
Aerospace 09 00353 g009
Figure 12. R-squared comparison of ALMESM and the original LME approximation with different globally unified γ under different numbers of sample points for approximating the LCF life of the GH720Li turbine disk.
Figure 12. R-squared comparison of ALMESM and the original LME approximation with different globally unified γ under different numbers of sample points for approximating the LCF life of the GH720Li turbine disk.
Aerospace 09 00353 g012
Figure 13. R-squared comparison of ALMESM, RBF, and KG under different numbers of sample points for approximating the LCF life of the GH720Li turbine disk.
Figure 13. R-squared comparison of ALMESM, RBF, and KG under different numbers of sample points for approximating the LCF life of the GH720Li turbine disk.
Aerospace 09 00353 g013
Figure 14. The frequency histogram of the LCF life.
Figure 14. The frequency histogram of the LCF life.
Aerospace 09 00353 g014
Figure 15. The reliability–life curve of the GH720Li turbine disk.
Figure 15. The reliability–life curve of the GH720Li turbine disk.
Aerospace 09 00353 g015
Table 1. The ALMESM algorithm flow.
Table 1. The ALMESM algorithm flow.
Algorithm Calculate   y ^ ( x )   for   Prediction   Point   x
Step 1 Define   sample   point   set   X = { x ( i ) , i = 1 , 2 , , N } ,   corresponding   response   y ( x ( i ) ) ,   and   the   given   prediction   point   x .
Step 2 Calculate   h i = min | x ( i ) x ( j ) | 2 .
Step 3 Initial   γ t   and   λ t ,   t t is the number of iterations.
Step 4 Calculate   β i t = γ i t / h i 2 .
Step 5 Calculate   λ * = arg min λ n log Z ( x ( i ) , λ )   by   PSO ,   and   Z ( x ( i ) , λ ) = i = 1 N exp [ β i t | x ( i ) x ( j ) | 2 + λ ( x ( i ) x ( j ) ) ] .
Step 6 Calculate   p i ( x ( i ) ) = exp [ β i t | x ( i ) x ( j ) | 2 + λ * ( x ( i ) x ( j ) ) ] Z ( x ( i ) , λ * )   at   each   sample   point   x ( i ) .
Step 7 Define   the   tolerance   p t o l .
Step 8 Determine   the   local   neighborhoods   X l o c a l   of   the   sample   point   by   selecting   p i ( x ( i ) ) > p t o l .
Step 9 Calculate   y ^ ( x ( i ) )   for   sample   point   x ( i ) X l o c a l ,   i = 1 , 2 , , k ,   y ^ ( x ( i ) ) = i = 1 k y ( x ( i ) ) p i ( x ( i ) ) .
Step 10 Calculate   i = 1 k ( y ( x ( i ) ) y ^ ( x ( i ) ) ) 2 ,   x ( i ) X l o c a l .
Step 11 Calculate   γ t + 1   and   λ t + 1   by   PSO ,   repeat   step   4 10 ,   until   γ i * = arg min i = 1 k ( y ( x ( i ) ) y ^ ( x ( i ) ) ) 2 ,   γ i * [ 0 , + ) ,   x ( i ) X l o c a l .
Step 12 Calculate   p i ( x )   at   each   prediction   point   x   in   local   sample   points   set   X l o c a l .
Step 13 Calculate   y ^ ( x )   for   prediction   point   x ,   y ^ ( x ) = i = 1 k y ( x ( i ) ) p i ( x ) ,   x ( i ) X l o c a l .
Table 2. Statistical analysis of sampling randomness of the LME and ALMESM.
Table 2. Statistical analysis of sampling randomness of the LME and ALMESM.
Sample Points R max 2 R min 2 R max 2 R min 2 ( R max 2 R min 2 ) / R max 2
LMEALMESMLMEALMESMLMEALMESMLME (%)ALMESM (%)
600.880.950.830.910.050.045.684.21
800.920.980.840.960.080.028.702.04
1000.940.990.900.970.040.024.262.02
1200.960.990.920.980.040.014.171.01
Table 3. Experimental designs for all test problems.
Table 3. Experimental designs for all test problems.
Sample Points SizePrediction Points Size
40 + 20i, i = 1, 2, …, 820 + 10i
Table 4. The material properties of GH720Li at 650 °C and the boundary conditions.
Table 4. The material properties of GH720Li at 650 °C and the boundary conditions.
Density ρ (kg/m3)Elastic Modulus E (Pa)Poisson’s Ratio υRotating Speed n (Rad/s)
81402.07 × 10110.373979.35
Table 5. The distribution and sampling scheme for 3 random variables.
Table 5. The distribution and sampling scheme for 3 random variables.
NameDistribution TypeMean (mm)Standard Deviation (mm)Sample PointsPrediction Points
d_R1Normal23.02750.00376540 + 20i,
i = 1, 2, …, 10
20 + 10i
d_R2Normal81.3850.01075
d_W6Normal10.0150.00246
Table 6. The typical life–reliability results.
Table 6. The typical life–reliability results.
MethodReliabilityComputational Cost (h)
50% (Cycles)90% (Cycles)99% (Cycles)99.87% (Cycles)
ALMESM151,025150,246149,686149,3433.5
Direct MC150,395149,025148,556148,11770.6
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fan, J.; Yuan, Q.; Jing, F.; Xu, H.; Wang, H.; Meng, Q. Adaptive Local Maximum-Entropy Surrogate Model and Its Application to Turbine Disk Reliability Analysis. Aerospace 2022, 9, 353. https://doi.org/10.3390/aerospace9070353

AMA Style

Fan J, Yuan Q, Jing F, Xu H, Wang H, Meng Q. Adaptive Local Maximum-Entropy Surrogate Model and Its Application to Turbine Disk Reliability Analysis. Aerospace. 2022; 9(7):353. https://doi.org/10.3390/aerospace9070353

Chicago/Turabian Style

Fan, Jiang, Qinghao Yuan, Fulei Jing, Hongbin Xu, Hao Wang, and Qingze Meng. 2022. "Adaptive Local Maximum-Entropy Surrogate Model and Its Application to Turbine Disk Reliability Analysis" Aerospace 9, no. 7: 353. https://doi.org/10.3390/aerospace9070353

APA Style

Fan, J., Yuan, Q., Jing, F., Xu, H., Wang, H., & Meng, Q. (2022). Adaptive Local Maximum-Entropy Surrogate Model and Its Application to Turbine Disk Reliability Analysis. Aerospace, 9(7), 353. https://doi.org/10.3390/aerospace9070353

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