
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

39237729
71431
10.1038/s41598-024-71431-1
Article
Machine learning-based predictions and analyses of the creep rupture life of the Ni-based single crystal superalloy
Zou Fan 12
Liu Pengjie pengjieliu_1112@163.com

1
Chen Yanzhan 1
Zhao Yaohua 1
1 https://ror.org/00f1zfq44 grid.216417.7 0000 0001 0379 7164 School of Traffic & Transportation Engineering, Central South University, Changsha, 410083 China
2 https://ror.org/04n3k2k71 grid.464340.1 0000 0004 1757 596X School of Intelligent Manufacturing and Mechanical Engineering, Hunan Institute of Technology, Hengyang, 421002 China
5 9 2024
5 9 2024
2024
14 2071617 9 2023
28 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
The evaluation of creep rupture life is complex due to its variable formation mechanism. In this paper, machine learning algorithms are applied to explore the creep rupture life span as a function of 27 physical properties to address this issue. By training several classical machine learning models and comparing their prediction performance, XGBoost is finally selected as the predictive model for creep rupture life. Moreover, we introduce an interpretable method, Shapley additive explanations (SHAP), to explain the creep rupture life predicted by the XGBoost model. The SHAP values are then calculated, and the feature importance of the creep rupture life yielded by the XGBoost model is discussed. Finally, the creep fracture life is optimized by using the chaotic sparrow optimization algorithm. We then show that our proposed method can accurately predict and optimize creep properties in a cheaper and faster way than other approaches in the experiments. The proposed method can also be used to optimize the material design across various engineering domains.

Keywords

Creep property prediction
XGBoost
Shapley additive explanations
Sparrow optimization algorithm
Subject terms

Computational methods
Mechanical engineering
Fundamental Research Funds for the Central Universities of Central South Universityissue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

The nickel-based single-crystal superalloy is the primary material that is used in aero-engine and gas turbine blades with excellent high-temperature performance1. To meet the design requirements of high-performance aero-engines, research and development activities with respect to nickel-based single-crystal superalloys are of great significance. In the current study, the creep property of the given material greatly limits the life cycle of the product. This is because the creep deformation and fracture of the material under high-temperature stress cause the product to fail2. Therefore, creep life is defined as the core research object of this paper because of its high scientific value. Even though the creep rupture life can be determined experimentally, the required measurements are often time-consuming and laborious due to lengthy creep tests and the high cost of alloy manufacturing, especially in cases when the alloy under study consists of noble metal elements such as Ru, Ti, and Ta3.

The creep property of an alloy is related to many variables that can be grouped into four categories: chemical composition factors, processing parameters, test conditions, and microstructural factors. The main chemical compositions of nickel-based superalloys are highly varied4, including nickel (Ni), rhenium (Re), cobalt (Co), chromium (Cr), and other alloy components. Additionally, the processing parameters, including the solution treatment time, heating temperature, test stress, and other parameters, are extremely complicated. The microstructural factors can be calculated from the chemical composition factors and processing parameters.

Since the experimental method is unable to meet the requirements of alloy design, several theoretical methods for accelerating the prediction of the creep fracture life spans of alloys have been proposed. Dang et al.5 performed creep fracture tests on three Mg–Al–Ca diecast alloys and used the Larson–Miller parameter (LMP) to predict the creep life. The results showed that the method performed well for specific temperature and stress levels. Li et al.6 investigated the creep properties of the K417 alloy at different high-temperature conditions and developed a physical analysis model to describe the influence of microstructural factors. As an example of using the time–temperature parameter (TTP) method, Larson and Miller7 developed a predictive model for creep life at low temperatures. Similarly, an isothermal method8 was proposed to describe the relationship between the creep property and stress under the conditions of constant temperature and changing stress for inferring the creep life of an alloy indirectly. In conclusion, the above traditional theoretical methods can achieve rapid creep life predictions, but they cannot reveal the complicated functional relationships between the creep life and its many influencing factors. Therefore, the accuracy of such models is limited.

With the development of big data and artificial intelligence, machine learning has been successfully applied to material performance prediction9–11, the discovery of new materials12, and other applications13–16. Data-driven approaches tend to be superior to other methods in terms of time efficiency and predictive performance. Since machine learning can learn the complex functional relationships between material properties and many factors, it is suitable for the task of predicting creep rupture life. Venkatesh and Rach17 developed and demonstrated a new method for life prediction using a BP (backward propagation) neural network. The results showed that the prediction accuracy of this method was significantly improved over that of linear models. By adopting the Bayesian neural network and Markov-chain Monte Carlo method, Yoo et al.18 predicted creep fracture life with an accuracy of 93.2% and ranked the importance of the constituent elements of the tested alloy.

This paper not only predicts creep life but also explains the results of the predicted creep life and finally optimizes the chemical compositions and processing parameters by an optimization algorithm. The remainder of the paper is structured as follows. In “Machine learning methods” section, we use machine learning algorithms to build a mapping function between creep rupture life and 27 features, including chemical composition factors, processing parameters and microstructural factors. In “Explanatory methods for interpreting creep rupture life” section, a novel method, Shapley additive explanations (SHAP), is introduced to interpret the predictive results regarding creep rupture life. In “Optimization based on improved sparrow search algorithm” section, the chaotic sparrow optimization algorithm is applied to improve the creep rupture life and achieve the best solution. Finally, conclusions are provided in “Conclusions” section. The overall framework of this paper is shown in Fig. 1.Fig. 1 Overall framework of this study.

Machine learning methods

Acquisition of the data structure

The creep rupture life dataset used in this study is obtained from previous publications19–21 and consists of 264 instances of superalloys. Each sample in the dataset contains 27 features describing its physical properties, which can be classified into the following four categories: chemical composition factors, processing parameters, test conditions and microstructural factors. A detailed description of the 27 features in the creep dataset is shown in Table 1. The distribution histograms of the 27 input features, as depicted in Fig. 2, reveal that certain features such as (1satt, 2satt, and 2satT) exhibit a high degree of variability and pronounced heterogeneity. Such characteristics can introduce several challenges to the predictive tasks in machine learning. The inherent complexity of these features may necessitate the employment of more sophisticated modeling techniques that can accommodate the diverse data patterns. Additionally, the variability could lead to increased model sensitivity to feature scaling and parameter tuning, requiring meticulous calibration to achieve optimal performance. Table 1 Description of the data features.

Index	Symbol	Description	Minimum content	Maximum content	
1	Ni	Nickel	57.4135	77.26	
2	Re	Rhenium	0	7.5	
3	Co	Cobalt	0	10	
4	Al	Aluminum	4.15	7.5	
5	Ti	Titanium	0	5	
6	W	Masurium	0	18.6	
7	Mo	Molybdenum	0	4.21	
8	Cr	Chromium	0.6	11.73	
9	Ta	Tantalum	0	12	
10	C	Carbon	0	0.1	
11	B	Boron	0	0.05	
12	Y	Yttrium	0	0.05	
13	Nb	Peloponium	0	2.59	
14	Hf	Hafnium	0	2	
15	Stt	Solution treatment time (h)	2	4	
16	1satt	The first-stage aging treatment time (h)	4	5	
17	2satt	The second-stage aging treatment time (h)	16	32	
18	StT	Solution treatment temperature (°C)	1180	1348	
19	1satT	The first-stage aging treatment temperature (°C)	982	1145	
20	2satT	The second-stage aging treatment temperature (°C)	704	899	
21	T	Test temperature (°C)	204	1800	
22	S	Test stress (MPa)	70	759	
23	Γ	Stacking fault energy (mJ/m2)	34.00666	119.4323	
24	DL	Diffusion coefficient (m2/s)	1.41E−25	0.00131	
25	G	Shear modulus (GPa)	52.20007	71.05828	
26	L	Lattice parameter (nm)	0.313652	0.362381	
27	Xγ′	Mole fraction of the γ′ phase	0.50187	0.90197	

Fig. 2 The histogram distribution for the 27 input features.

Data standardization

Data standardization is often the first step in machine learning tasks. As shown in Table 1, a total of 27 features that affect creep life vary widely in terms of their value ranges. For example, the value of the diffusion coefficient DL ranges from 1.41E−25 to 0.00131 m2/s, while the solution treatment temperature ranges from 1180 to 1348 °C. Therefore, this may become an obstacle during the model training process. Here, we use the min–max scaling method22 to achieve isometric scaling of the original data. The min–max scaling method is employed to preprocess the 27 features, hence, all descriptors are mapped to [0, 1]. The method is defined as follows:1 X∗=X-XminXmax-Xmin,

where Xmin and Xmax represent the lower and upper limits of a feature X, respectively. Furthermore, from the original dataset, it can be seen that the creep life spawns of the samples are distributed in the interval from 30 to 5180 h. We use logarithmic scaling for the creep life, as shown in Eq. (2). This enables us to efficiently convert the predicted creep life. Finally, the value of the feature after data normalization is used as an input parameter for the machine learning model mentioned in Table 2.2 Y∗=logY.

XGBoost algorithm

Before the machine learning task is completed, we perform a Pearson correlation analysis of the creep dataset. Figure 3 shows a heat map that is an indication of the relationships between the 27 features and the creep rupture life. This is, however, only a simple linear analysis and may not be able to explore the more complex nonlinear relationships between the variables. Therefore, several machine learning models with strong representational capabilities are discussed in the following sections.Fig. 3 Heat map based on the Pearson correlation coefficient.

XGBoost (eXtreme Gradient Boosting) is a successful machine learning library based on the gradient boosting algorithm proposed by Chen23. It has been extensively used in industry due to the portability of the distributed version and its powerful predictive capabilities. As an improved GBDT (Gradient Boosting Decision Tree) algorithm24, XGBoost consists of a large number of decision trees strung together to handle classification and regression tasks accurately and efficiently. Compared with GBDT, XGBoost is mainly optimized in the following aspects. First, XGBoost adds a second-order Taylor expansion to the cost function, which improves the accuracy of the algorithm. Second, the regular penalty term25 is utilized in the cost function to reduce the model complexity and prevent overfitting. The Fig. 4 visually illustrates the procedural workflow of the XGBoost algorithm, beginning with the initial step of “Bootstrap aggregation”. This phase involves the utilization of resampling techniques to generate multiple subsets from the original dataset, encompassing a diverse range of classes. These subsets serve as the foundation for the subsequent stage known as ‘CART Trees’, which includes classification and regression trees, are binary models that are constructed and intertwined with each of the data subsets. The subsequent phase, termed ‘Weighting increase’, embodies the iterative nature of gradient boosting. This process dynamically adjusts the decision boundaries by focusing on and augmenting the weights assigned to misclassified points from previous iterations, thereby enhancing the overall accuracy of the model. The final stage cleverly combines the fitting results obtained from all the CART Trees generated throughout the workflow. By employing a weighted average of these results, the XGBoost model is capable of formulating an optimized and robust prediction.Fig. 4 The schematic diagram of XGBoost.

Empirical analysis and comparisons with other models

In this study, the dataset is split into training and test sets using a 10-fold cross-validation method. This method randomly divides the dataset into 10 subsets. Each time one of these subsets is considered as the test set, the remaining 9 subsets are used as the training sets. Using this approach, the machine learning model can be trained and tested 10 times to obtain a more comprehensive evaluation index. To reflect the accuracy of the predictive models, three evaluation parameters are adopted in this paper26, including the coefficient of determination (R2), mean absolute error (MAE), and root mean squared error (RMSE). Suppose that y1,y2,⋯yn are the actual values, y~1,y~2,⋯y~n are the predicted values, and y¯ is the mean of yi; then,3 R2=1-∑i=1ny~i-yi2∑i=1ny¯-yi2,

4 MAE=1n∑i=1ny~i-yi,

5 RMSE=1n∑i=1ny~i-yi2.

Without the loss of generality, we establish other classic predictive models for a performance comparison, including a decision tree27, a support vector machine28, a neural network29, a random forest30 and LightGBM31 and so on. Table 2 provides a brief description of each model. Table 2 Introduction to different machine learning models.

Models	Notes	
LR	Linear regressor	
KNN	K nearest neighbor	
SVM	Support vector machine	
ETR	Extra tree regressor	
DT	Decision tree	
GB	Gradient boosting regressor	
NN	Neural network	
RF	Random Forest	
LightGBM	Light gradient boosting machine	
XGBoost	eXtreme gradient boosting	

Here, we implement the methods in Table 2 using the Scikit-learn library32 with version 3.7 of Python, which is available at https://www.python.org/. In addition, XGBoost and LightGBM are implemented in Python using the XGBoost and LightGBM libraries, respectively. To ensure optimal performance of these models, we utilized the Optuna library33 for hyperparameter tuning. Optuna is a state-of-the-art framework for automated hyperparameter optimization, which allowed us to systematically search for the best hyperparameters for each model. The test platform includes a device equipped with an AMD Ryzen-7 4800H CPU and 16 GB of RAM. The above 10 machine learning models are trained and tested separately by employing 10-fold cross-validation.

The line chart in Fig. 5 shows the R2, RMSE and MAE distributions corresponding to the 10 machine learning models used in this paper. Figure 5 shows that complex machine learning models (i.e., NN, GB, RF, LightGBM and XGBoost) achieve high prediction accuracy, and their corresponding R2 values are larger than 0.95. The results also show that the prediction performance of LightGBM and XGBoost is better than that of the other models. In particular, for XGBoost, R2=0.9759, RMSE=0.1428, and MAE=0.1058; these results are better than those of the other models. It is also seen that the performance of the three relatively simple machine learning models (LR, KNN, and SVM) is poor, and their R2 values are low. Due to the obvious nonlinear relationships between the creep rupture life and the 27 features examined in this study, these three simple models are simply unable to reveal the complex mapping function relationships.Fig. 5 Results of the creep rupture life prediction models, including their R2, RMSE and MAE values. A model with a larger R2 value and lower of RMSE and MAE values has higher fitting accuracy.

Finally, we utilize scatter plots to illustrate the fitting performance of the top three best-performing models and the bottom three underperforming models. As shown in Fig. 6, the experimental value of the creep rupture life is represented on the X-axis, and the predictions of each model are shown on the Y-axis. Each graph contains a 1:1 straight line representing a perfect fit. The closer the scatter points to the straight line, the better the fit of the model. It can be seen that the prediction effect of the XGBoost model is the best. Therefore, in this paper, we use XGBoost as the final model for the prediction of the creep life.Fig. 6 Predictive performance of the six machine learning models on the creep dataset. Note: The X-axis and Y-axis represent the actual and predicted values, respectively. The color bar represents the absolute error between predicted and actual values.

Explanatory methods for interpreting creep rupture life

In machine learning, accuracy, complexity, and interpretability are the three most important concepts. In general, complexity and interpretability are in opposition to each other34. The XGBoost model used in this research has strong fitting capabilities due to its complex structure. However, due to its large number of parameters and complex working mechanisms, XGBoost is often regarded as a black-box model. To reveal the internal mechanisms of the 27 features that affect the creep rupture life, this section attempts to explain the XGBoost model using model interpretation methods.

Principles of the SHAP method

Model interpretability is divided into two categories: global interpretability and local interpretability. Global interpretability is based on the entire studied dataset to help people understand the internal working mechanism and complex logic of the model, while local interpretability analyzes the specific decision-making process of the model for a single sample. In this paper, SHAP is used to improve the interpretability of the machine learning model. SHAP is a method proposed by Lundberg et al.35 that provides global and local interpretations and further possesses a solid theoretical foundation relative to those of other model interpretation methods. SHAP assigns an importance value to each feature of the test learning model. This facilitates the provision of a reliable explanation for a complex predictive model. The key idea of SHAP is to use a surrogate model g to represent the black-box model f that is hard to directly explain. The explanatory model g is defined as a simple linear function, as expressed in Eq. (6):6 fx≈gz′=ϕ0+∑i=1Mϕizi′,

where f(x) is the original complex machine learning model’s prediction function, which takes an input x and outputs a prediction, g(z′) is an approximation of the function f(x), which is constructed to be more interpretable. It is a linear function of the features z′, which are some transformed or perturbed versions of the original features x. zi′∈0,1M denotes a feature being observed (zi′=1) and the unknown feature (zi′=0); M represents the total number of features, which is equal to 27 in this paper; ϕi represents the contribution value of the i-th feature to the model prediction result, and the sum of each feature’s contribution is similarly equal to the output f(x) of the original model.

In SHAP, the Shapley value explanation is represented as an additive feature attribution method, i.e., a linear model. This view also connects the LIME (Local Interpretable Model-Agnostic Explanations)36 and Shapley values37. In previous studies, it was shown that Shapley values represent a unique solution that can satisfy three properties at the same time: symmetry, the dummy property, and additivity. To calculate the Shapley value of each feature, the mathematical expectation of the function conditioned on a subset of the input features is defined as E[f(x)|xS]. The contribution value of the i-th feature (ϕi) is a weighting of all mathematical expectations, as shown in Eq. (7):7 ϕi=∑S⊆N\{i}S!(N-S-1)!N![fx(S⋃{i})-fx(S)],

where ϕi represents the average contribution of feature i to the prediction function, the summation is taken over all subsets S of the set N excluding feature i. The term S denotes the cardinality of subset S, and N is the total number of features. The function fx(S) is the prediction value of the model when only the features in subset S are considered. The coefficient S!(N-S-1)!N! is a weight that reflects the proportion of subsets that include S in all possible subsets of N, ensuring that each feature's contribution is fairly distributed.

Local interpretation based on individual samples

The most prominent advantage of SHAP is that it can correspond to the influence of the features in each sample, and it can also show its positive and negative effects of these features. First, we show how SHAP values are locally attributed to individual samples in the creep dataset.

We take the first sample in the dataset as an example to explore the impacts of these features on the model output. Figure 7 shows how each feature contributes to the predictive result by pushing it from the base value, which is defined as the average value of all sample predicted values in the dataset. In this example, this base value is equal to 5.93, while the predictive result of the XGBoost model is equal to 4.30. As shown in Fig. 7, the red and blue arrows indicate the positive and negative SHAP values of the features, respectively. On the right side of the output value (f(x)), the feature represented in blue in the graph pushes the predicted value toward lower values, while on the left side of the output value, the feature represented in red pushes the predicted value toward larger values. The longer the length of the arrow in the figure, the greater the contribution of the corresponding feature is. We see that 2satt (17.0) and 1satt (4.0) are related to the negative SHAP values, and a consequent reduction in the creep rupture life is expected. On the other hand, Re (0.16) and especially T (256.0) are related to the positive SHAP values, and the creep rupture life is expected to increase.Fig. 7 Explanation for the first sample obtained using the SHAP method.

Global interpretation based on all the samples

The global interpretation is the aggregation of the local interpretation information of all samples so that it can fully reflect how each feature affects the creep life. For each sample, a set of SHAP values can be generated, and Fig. 8 shows the distribution of SHAP values for each feature based on the entire dataset. The Y-axis in the graph represents all features are ranked by their global contributions ∑j=1Nϕi(j), and the X-axis represents the SHAP values ϕi(j). Each point represents a sample, from low (blue) to high (red), with different colors representing the values of the features. In addition, the importance of each feature is ranked from top to bottom. In Fig. 8, it can be seen that XGBoost generates the following four important factors for creep rupture life prediction: the test temperature (T); second-stage aging treatment time (2satt); first-stage aging treatment time (1satt); and mass percent of Cr (Cr). Furthermore, the test temperature is an extremely important feature that is negatively correlated with the creep rupture life. For a high test temperature (such dots in the graph are shown in red), the SHAP values are negative, and for a low test temperature (such dots in the graph are shown in blue), the SHAP values are positive. Shortening the experimental time can help to increase the creep rupture life. According to the same analysis, the contribution of each feature to the creep life can be discovered. It is worth mentioning that although the second-stage aging treatment time is negatively correlated with the creep rupture life, it is seen in Fig. 8 that a short second-stage aging time results in a reduction in the creep life. In terms of chemical composition, Cr has the most obvious effects on the creep rupture life. Figure 8 shows that the creep fracture life is negatively correlated with Cr content because an increase in Cr content compromises the phase stability of the alloy and thus affects the creep resistance of the alloy. In addition, we can also find that increasing Ta, Re, and W, three refractory elements, results in increasing the creep rupture life. This is because these elements can increase the dissolution temperature of the γ′ phase and improve the temperature bearing capacity. To quantitatively analyze the contribution of each feature, we obtain the average of the absolute values of the SHAP values for each feature, and the feature importance rankings are shown in Fig. 9. These feature importance rankings not only enhance the transparency and interpretability of the black-box machine learning models but also guide the design of further experiments for optimizing creep rupture life.Fig. 8 SHAP summary plot of the XGBoost model.

Fig. 9 Feature importance rankings for predicting the creep rupture life.

Optimization based on improved sparrow search algorithm

In this section, the XGBoost model is used as the mapping function between the creep life and the main variables, which is also the objective function. We then establish a single-objective optimization model that maximizes the creep rupture life. To enhance the global search ability and convergence speed of the model, a chaotic sparrow search optimization algorithm is proposed by importing the Tent chaotic sequence and Gaussian mutation.

Model establishment

Here, we choose all 27 features as decision variables and XGBoost as the mapping function. Let X=[X1,X2,X3,⋯Xn](n=27) denote 27 features which can be classified into four categories, that is, chemical composition factors, processing parameters, test conditions and microstructural factors. Furthermore, fX demonstrates the creep rupture life predicted by the XGBoost model. As a result, the optimization algorithm can be constructed as the following mathematical model.maxfX

8 s.tX1min≤X1≤X1max,X2min≤X2≤X2max,X3min≤X3≤X3max,⋮Xnmin≤Xn≤Xnmaxn=27,∑i=114Xi=100.

where Xnmin and Xnmax are the lower and upper bounds according to the raw dataset, respectively, and their specific values are shown in Table 1. Notably, an equality constraint is established to ensure that the total sum of the 14 chemical element contents equals 100%. To solve this model, we employ an improved heuristic optimization algorithm which will described in the next section.

Solution of the model

How to find an optimal set of parameter combinations in a high-dimensional search space is the problem to be solved by combinatorial optimization, and this section uses an optimization algorithm to determine the optimal combination of chemical composition factors and processing parameters. Currently, heuristic optimization algorithms are widely applied in the engineering field because of their simplicity and practicality. The main idea of such algorithms is to search for the optimal solution within a certain range by simulating the collective behavior of decentralized systems, e.g., in nature. Many swarm intelligence optimization algorithms have been proposed based on the swarm behaviors of ants, birds, bees, wolves, whales, sparrows, and other intelligent creatures. The sparrow search algorithm (SSA) proposed by Xue et al.38 in 2020 is a new type of swarm intelligence optimization algorithm. Compared with other algorithms, such as the particle swarm optimization (PSO)39, ant colony optimization (ACO)40, and gray wolf optimization (GWO) algorithms41, it provides higher search accuracy, faster convergence, and higher levels of stability and robustness. Therefore, here, we use an improved SSA algorithm to improve the creep rupture life.

Standard SSA algorithm

The sparrow search algorithm (SSA) is inspired by the foraging and antipredation behaviors of sparrows. Its core idea is to abstract the foraging process of sparrows as a producer-scrounger model with reconnaissance and early warning mechanisms.

The producers are highly adaptable and have a wide search range, guiding the population to search and forage. To obtain better fitness, scroungers follow the producers for food. At the same time, to improve the success rate of predation, some scroungers monitor the producers or forage around them. In cases where the entire population faces predators or realizes danger, the group immediately carries out collective antipredation behavior.

The core of the SSA algorithm is the producer-scrounger model, and its modeling process is as follows. Assuming that there are N sparrows in a D-dimensional search space, the position of the i-th sparrow in the D-dimensional search space is Xi=[Xi1,Xi2,⋯XiD](i=1,2,⋯,N), where Xid represents the position of the i-th sparrow in the d-th dimension.

The producers often account for 10% to 20% of the population, and the equation for position updating is:9 xidt+1=xidt·exp(-iα·T),R2<STxidt+Q·L,R2≥ST,

where t represents the current number of iterations, T is the maximum number of iterations, α is a uniform random number between 0 and 1, Q is a Gaussian random number, L is a 1 × d matrix with all entities equal to 1, and R2∈[0,1] and ST∈[0.5,1] represent the warning and safety values, respectively. In cases where R2<ST, the population does not detect the presence of predators or other dangers, and the search environment is safe. Producers can extensively search to guide the population to obtain a higher degree of fitness. When R2≥ST, a sparrow detects a predator and releases the danger signal, and the population immediately performs antipredation behavior, adjusts its collective search strategy, and quickly moves closer to a safe area.

Aside from the producers, the remaining sparrows act as scroungers and update their positions according to the following formula:10 xidt+1=Q·expxwdt-xidti2,i>n2xbdt+1+1D∑d=1Drand-1,1·xidt-xbdt+1,i≤n2,

where xwdt represents the worst position of the given sparrow in the d-th dimension during the t-th iteration of the population, and xidt+1 represents the optimal position of the sparrow during the t+1-th iteration of the population. For i>n/2, the i-th scrounger does not obtain food and is in a state of starvation and low fitness. To obtain higher energy, it needs to fly to other places for foraging. In cases where i≤n/2, the i-th scrounger randomly finds a location near the current optimal foraging location.

Chaotic sparrow search optimization algorithm

The CSSA (chaotic sparrow search optimization algorithm) algorithm is improved by using Tent chaos and Gaussian mutation. Chaos, as a nonlinear phenomenon in nature, is widely used to optimize search problems due to the randomness, ergodicity, and regularity of the underlying chaotic variables42. This effectively maintains the diversity of the population, further helps the algorithm jump out of local optima and improves the global search performance. Through rigorous mathematical derivation, it is shown that the traversal uniformity and convergence speed of the Tent map are better than those of other forms of chaotic maps. Therefore, the CSSA algorithm uses Tent mapping to generate the chaotic sequence of the optimized algorithm. The Tent mapping expression is shown in Eq. (11). It randomly generates an initial value zi between 0 and 1 and then uses the iterative process in Eq. (11) to generate the Tent sequence.11 zi+1=2zi,0≤z≤122(1-zi),12<z≤1.

The Gaussian variation is obtained through a Gaussian distribution, which specifically refers to replacing the original parameter value with a random number conforming to a normal distribution with mean u and variance σ2 during the mutation process.12 mutationx=x1+N0,1,

where x is the original parameter value, N(0,1) represents a normal random variable with a mean of 0 and a standard deviation of 1, and mutation(x) is the value obtained after Gaussian mutation. The normal distribution characteristics show that the key search area of Gaussian mutation is a local area near the original individual. This helps the algorithm find the global minimum point with high efficiency and high precision and improves the robustness of the algorithm43. In summary, the CSSA algorithm introduces the Tent chaotic search and Gaussian mutation operations. The pseudocode of this algorithm is presented below.

Results and discussion

To verify the optimized performance of the CSSA algorithm, we conduct several additional sets of optimization experiments based on the sparrow search algorithm (SSA), gray wolf optimization (GWO) algorithm, and particle swarm optimization (PSO) algorithm. To compare the convergence speeds and optimization performances of the four algorithms, we obtain the curve of the fitness value (ln(creep rupture life)) versus the number of iterations for each algorithm. As shown in Fig. 10, the CSSA has the fastest convergence rate, and its optimal solution is the best. It is also seen that the PSO and GWO algorithm evolve extremely slowly and tend to quickly reach a local optimum. Despite SSA converging faster than GWO and PSO, this premature convergence leads to it obtaining a suboptimal solution. In addition, the stability of the selected optimization algorithm is also an important criterion. Therefore, we perform 30 separate runs for each algorithm and set the number of iterations to 1000. The results in Table 3 show that the CSSA algorithm has the highest convergence rate and the lowest standard deviation. In addition, the best solution and average solution of the CSSA algorithm are the highest, and this also confirms the optimization efficiency of the CSSA. Figure 11 shows the distributions of the results for the considered optimization algorithms. The stability of the CSSA algorithm is higher than that of the other three algorithms.Fig. 10 Iteration curves of the four optimization algorithms.

Table 3 Statistical results obtained by four different optimization algorithms.

Optimization algorithm	Convergence rate	The best solution	The average solution	Standard deviation	
PSO	25/30	8.8733	8.6878	0.16928	
GWO	26/30	8.8121	8.6429	0.17557	
SSA	25/30	8.5424	8.3536	0.18396	
CSSA	27/30	8.9014	8.7559	0.12881	

Fig. 11 Distributions of the optimized results yielded by four algorithms.

Table 4 shows the best combination of these 27 features obtained by the CSSA algorithm. Initially, we observe that the sum of the chemical element contents (i.e., the numerical values of the first fourteen features) comprising Ni-based single crystal superalloy amounts to 100%, fulfilling the equality constraint specified in model (8). Aside from Ni, Ta has the largest proportion of approximately 11.872%. This is also consistent with our conclusion in “Explanatory methods for interpreting creep rupture life” section that the creep life is positively correlated with the Ta content. Therefore, in industrial production processes, increasing the Ta content is essential for improving the creep rupture life. Moreover, the processing parameters are equally important to the creep rupture life, such as the second-stage aging treatment time (2satt) and the first-stage aging treatment time (1satt). This further suggests that adjustment of 1satt and 2satt enhances the creep resistance of the tested alloy. It should be noted that when all the features are fully considered, a satisfactory creep rupture life can be achieved. Table 4 One of the optimization results of the CSSA.

Ni (%)	Re (%)	Co (%)	Al (%)	Ti (%)	W (%)	Mo (%)	Cr (%)	Ta (%)	
60.199	2.163	0.101	4.150	3.227	10.485	0.674	4.399	11.872	
C (%)	B (%)	Y (%)	Nb (%)	Hf (%)	Stt (h)	1satt (h)	2satt (h)	StT (°C)	
0.025	0.048	0.025	2.444	0.188	3.139	4.138	22.524	1299.178	
1satT (°C)	2satT (°C)	T (°C)	S (MPa)	Γ (mJ/m2)	DL (m2/s)	G (GPa)	L (nm)	Xγ′ (–)	
1104.011	885.346	893.655	119.791	41.332	0.001	67.710	0.343	0.596	

Conclusions

In the field of developing new materials and improving material compositions, the data-driven approach offers irreplaceable advantages. Many of the previous studies on creep rupture life focused on performing large numbers of experiments, but this is not cost-effective and requires considerable time. With the accumulation of industrial big data, machine learning can provide a simpler and more efficient solution for creep life prediction and optimization. In this paper, we establish 10 different machine learning models to predict creep fracture life and finally select XGBoost by comparing the accuracy rates of the tested models. To gain insights into black-box ML models, we then propose using the SHAP model to explain the working mechanism of creep rupture life. The process of interpretation is divided into two steps: (1) visualizing the positive and negative feature contributions of each feature using the corresponding SHAP values and (2) ranking the importance levels of features based on global interpretation methods. Ultimately, we obtain an optimal combination of 27 features with the chaotic sparrow search optimization algorithm. In summary, the proposed method can provide theoretical guidance for improving creep life. In future works, we will extend the approach in this paper by incorporating other parameters, such as activation energy and lattice misfit.

Acknowledgements

The authors wish to thank the supported by the Fundamental Research Funds for the Central Universities of Central South University.

Author contributions

Zou Fan initiated the study and participated in data collection, statistical analysis, and manuscript writing. Pengjie Liu was responsible for data collection, program writing, and manuscript writing. Yanzhan Chen were responsible for the design and implementation of optimization algorithms. Yaohua Zhao initiated the study, supervised it, and revised the draft manuscripts. All authors approved the final manuscript. Information is available upon request from the corresponding author.

Data availability

The datasets used during the current study are available from the corresponding author on reasonable request.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Zhang Y Xu X Lattice misfit predictions via the Gaussian process regression for Ni-based single crystal superalloys Met. Mater. Int. 2020 10.1007/s12540-020-00883-7
Zhang, Y. & Xu, X. Lattice misfit predictions via the Gaussian process regression for Ni-based single crystal superalloys. Met. Mater. Int.10.1007/s12540-020-00883-7 (2020).10.1007/s12540-020-00883-7
2. Choi BG Kim IS Hong HU Do J Jung JE Jo CY Effect of Ti content on creep properties of Ni-base single crystal superalloys Met. Mater. Int. 2017 23 877 883 10.1007/s12540-017-7089-7
Choi, B. G. et al. Effect of Ti content on creep properties of Ni-base single crystal superalloys. Met. Mater. Int. 23, 877–883. 10.1007/s12540-017-7089-7 (2017).10.1007/s12540-017-7089-7
3. Rajan K Materials informatics Mater. Today 2005 8 38 45 10.1016/S1369-7021(05)71123-8
Rajan, K. Materials informatics. Mater. Today 8, 38–45. 10.1016/S1369-7021(05)71123-8 (2005).10.1016/S1369-7021(05)71123-8
4. Reed RC Tao T Warnken N Alloys-by-design: Application to nickel-based single crystal superalloys Acta Mater. 2009 57 5898 5913 10.1016/j.actamat.2009.08.018
Reed, R. C., Tao, T. & Warnken, N. Alloys-by-design: Application to nickel-based single crystal superalloys. Acta Mater. 57, 5898–5913. 10.1016/j.actamat.2009.08.018 (2009).10.1016/j.actamat.2009.08.018
5. Terada Y Sato T Assessment of creep rupture life of heat resistant Mg–Al–Ca alloys J. Alloys Compd. 2010 504 261 264 10.1016/j.jallcom.2010.05.108
Terada, Y. & Sato, T. Assessment of creep rupture life of heat resistant Mg–Al–Ca alloys. J. Alloys Compd. 504, 261–264. 10.1016/j.jallcom.2010.05.108 (2010).10.1016/j.jallcom.2010.05.108
6. Li S Wang B Shi D Yang X Qi H A physically based model for correlating the microstructural degradation and residual creep lifetime of a polycrystalline Ni-based superalloy J. Alloys Compd. 2019 783 565 573 10.1016/j.jallcom.2018.11.417
Li, S., Wang, B., Shi, D., Yang, X. & Qi, H. A physically based model for correlating the microstructural degradation and residual creep lifetime of a polycrystalline Ni-based superalloy. J. Alloys Compd. 783, 565–573. 10.1016/j.jallcom.2018.11.417 (2019).10.1016/j.jallcom.2018.11.417
7. Dang YY Zhao XB Yuan Y Predicting long-term creep-rupture property of Inconel 740 and 740H Mater. High Temp. 2016 33 1 1 5 10.1179/1878641315Y.0000000010
Dang, Y. Y. et al. Predicting long-term creep-rupture property of Inconel 740 and 740H. Mater. High Temp. 33(1), 1–5. 10.1179/1878641315Y.0000000010 (2016).10.1179/1878641315Y.0000000010
8. Bolton J Reliable analysis and extrapolation of creep rupture data Int. J. Press. Vessels Pip. 2017 157 1 19 10.1016/j.ijpvp.2017.08.001
Bolton, J. Reliable analysis and extrapolation of creep rupture data. Int. J. Press. Vessels Pip. 157, 1–19. 10.1016/j.ijpvp.2017.08.001 (2017).10.1016/j.ijpvp.2017.08.001
9. Hong D Kwon S Yim C Exploration of machine learning to predict hot ductility of cast steel from chemical composition and thermal conditions Met. Mater. Int. 2020 10.1007/s12540-020-00713-w
Hong, D., Kwon, S. & Yim, C. Exploration of machine learning to predict hot ductility of cast steel from chemical composition and thermal conditions. Met. Mater. Int.10.1007/s12540-020-00713-w (2020).10.1007/s12540-020-00713-w
10. Thankachan T Soorya Prakash K Kavimani V Silambarasan SR Machine learning and statistical approach to predict and analyze wear rates in copper surface composites Metals Mater. Int. 2020 27 220 10.1007/s12540-020-00809-3
Thankachan, T., Soorya Prakash, K., Kavimani, V. & Silambarasan, S. R. Machine learning and statistical approach to predict and analyze wear rates in copper surface composites. Metals Mater. Int. 27, 220. 10.1007/s12540-020-00809-3 (2020).10.1007/s12540-020-00809-3
11. Lin YC Yang H Chen D-D He D-G Stacked auto-encoder network to predict tensile deformation behavior of a typical nickel-based superalloy considering Portevin–Le Chatelier effects Met. Mater. Int. 2019 10.1007/s12540-019-00435-8
Lin, Y. C., Yang, H., Chen, D.-D. & He, D.-G. Stacked auto-encoder network to predict tensile deformation behavior of a typical nickel-based superalloy considering Portevin–Le Chatelier effects. Met. Mater. Int.10.1007/s12540-019-00435-8 (2019).10.1007/s12540-019-00435-8
12. Liu Y Zhao T Ju W Shi S Materials discovery and design using machine learning J. Materiom. 2017 3 159 177 10.1016/j.jmat.2017.08.002
Liu, Y., Zhao, T., Ju, W. & Shi, S. Materials discovery and design using machine learning. J. Materiom. 3, 159–177. 10.1016/j.jmat.2017.08.002 (2017).10.1016/j.jmat.2017.08.002
13. Jain DK Jain R Lan X Upadhyay Y Thareja A Driver distraction detection using capsule network Neural Comput. Appl. 2020 1 1 14 10.1007/s00521-020-05390-9
Jain, D. K., Jain, R., Lan, X., Upadhyay, Y. & Thareja, A. Driver distraction detection using capsule network. Neural Comput. Appl. 1, 1–14. 10.1007/s00521-020-05390-9 (2020).10.1007/s00521-020-05390-9
14. Jain DK Lan X Manikandan R Fusion of iris and sclera using phase intensive rubbersheet mutual exclusion for periocular recognition Image Vis. Comput. 2020 103 104024 10.1016/j.imavis.2020.104024
Jain, D. K., Lan, X. & Manikandan, R. Fusion of iris and sclera using phase intensive rubbersheet mutual exclusion for periocular recognition. Image Vis. Comput. 103, 104024. 10.1016/j.imavis.2020.104024 (2020).10.1016/j.imavis.2020.104024
15. Jain D Jain R Upadhyay Y Kathuria A Lan X Deep refinement: Capsule network with attention mechanism-based system for text classification Neural Comput. Appl. 2020 32 1 10.1007/s00521-019-04620-z
Jain, D., Jain, R., Upadhyay, Y., Kathuria, A. & Lan, X. Deep refinement: Capsule network with attention mechanism-based system for text classification. Neural Comput. Appl. 32, 1. 10.1007/s00521-019-04620-z (2020).10.1007/s00521-019-04620-z
16. Jain D Kumar A Garg G Sarcasm detection in mash-up language using soft-attention based bi-directional LSTM and feature-rich CNN Appl. Soft Comput. 2020 91 106198 10.1016/j.asoc.2020.106198
Jain, D., Kumar, A. & Garg, G. Sarcasm detection in mash-up language using soft-attention based bi-directional LSTM and feature-rich CNN. Appl. Soft Comput. 91, 106198. 10.1016/j.asoc.2020.106198 (2020).10.1016/j.asoc.2020.106198
17. Venkatesh V Rack HJ A neural network approach to elevated temperature creep–fatigue life prediction Int. J. Fatigue 1999 21 225 234 10.1016/S0142-1123(98)00071-1
Venkatesh, V. & Rack, H. J. A neural network approach to elevated temperature creep–fatigue life prediction. Int. J. Fatigue 21, 225–234. 10.1016/S0142-1123(98)00071-1 (1999).10.1016/S0142-1123(98)00071-1
18. Yoo YS Jo CY Jones CN Compositional prediction of creep rupture life of single crystal Ni base superalloy by Bayesian neural network Mater. Sci. Eng. A 2002 336 22 29 10.1016/S0921-5093(01)01965-7
Yoo, Y. S., Jo, C. Y. & Jones, C. N. Compositional prediction of creep rupture life of single crystal Ni base superalloy by Bayesian neural network. Mater. Sci. Eng. A 336, 22–29. 10.1016/S0921-5093(01)01965-7 (2002).10.1016/S0921-5093(01)01965-7
19. Liu Y Wu J Wang Z Lu X-G Avdeev M Shi S Wang C Yu T Predicting creep rupture life of Ni-based single crystal superalloys using divide-and-conquer approach based machine learning Acta Mater. 2020 195 454 467 10.1016/j.actamat.2020.05.001
Liu, Y. et al. Predicting creep rupture life of Ni-based single crystal superalloys using divide-and-conquer approach based machine learning. Acta Mater. 195, 454–467. 10.1016/j.actamat.2020.05.001 (2020).10.1016/j.actamat.2020.05.001
20. Conduit BD Jones NG Stone HJ Conduit GJ Design of a nickel-base superalloy using a neural network Mater. Des. 2017 131 358 365 10.1016/j.matdes.2017.06.007
Conduit, B. D., Jones, N. G., Stone, H. J. & Conduit, G. J. Design of a nickel-base superalloy using a neural network. Mater. Des. 131, 358–365. 10.1016/j.matdes.2017.06.007 (2017).10.1016/j.matdes.2017.06.007
21. Yamazaki, M., Yamagata, T. & Harada, H. Nickel-Base Single Crystal Superalloy and Process for Production Thereof, US Patent, 4707192 1–12 (1987).
22. Milligan GW Cooper MC A study of standardization of variables in cluster analysis J. Classif. 1988 5 181 204 10.1007/BF01897163
Milligan, G. W. & Cooper, M. C. A study of standardization of variables in cluster analysis. J. Classif. 5, 181–204. 10.1007/BF01897163 (1988).10.1007/BF01897163
23. Chen, T. & Guestrin, C. Xgboost: A scalable tree boosting system. In Proc. 22nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining 785–794. 10.1145/2939672.2939785 (ACM, 2016).
24. Rao H Shi X Rodrigue AK Feng J Xia Y Elhoseny M Yuan X Gu L Feature selection based on artificial bee colony and gradient boosting decision tree Appl. Soft Comput. 2019 74 634 642 10.1016/j.asoc.2018.10.036
Rao, H. et al. Feature selection based on artificial bee colony and gradient boosting decision tree. Appl. Soft Comput. 74, 634–642. 10.1016/j.asoc.2018.10.036 (2019).10.1016/j.asoc.2018.10.036
25. Sutskever I Hinton GE Krizhevsky A Imagenet classification with deep convolutional neural networks Adv. Neural. Inf. Process. Syst. 2012 10.1145/3065386
Sutskever, I., Hinton, G. E. & Krizhevsky, A. Imagenet classification with deep convolutional neural networks. Adv. Neural. Inf. Process. Syst.10.1145/3065386 (2012).10.1145/3065386
26. Weber G Pinz M Ghosh S Machine learning-aided parametrically homogenized crystal plasticity model (PHCPM) for single crystal Ni-based superalloys JOM 2020 72 4404 4419 10.1007/s11837-020-04344-9
Weber, G., Pinz, M. & Ghosh, S. Machine learning-aided parametrically homogenized crystal plasticity model (PHCPM) for single crystal Ni-based superalloys. JOM 72, 4404–4419. 10.1007/s11837-020-04344-9 (2020).10.1007/s11837-020-04344-9
27. Rutkowski L Jaworski M Pietruczuk L Duda P The CART decision tree for mining data streams Inf. Sci. 2014 266 1 15 10.1016/j.ins.2013.12.060
Rutkowski, L., Jaworski, M., Pietruczuk, L. & Duda, P. The CART decision tree for mining data streams. Inf. Sci. 266, 1–15. 10.1016/j.ins.2013.12.060 (2014).10.1016/j.ins.2013.12.060
28. Abe S Abe S Variants of support vector machines Support Vector Machines for Pattern Classification 2010 Springer 163 226
Abe, S. Variants of support vector machines. In Support Vector Machines for Pattern Classification (ed. Abe, S.) 163–226 (Springer, 2010).
29. Liu C-L Fink GA Govindaraju V Jin L Special issue on deep learning for document analysis and recognition Int. J. Doc. Anal. Recogn. 2018 21 3 159 160 10.1007/s10032-018-0310-5
Liu, C.-L., Fink, G. A., Govindaraju, V. & Jin, L. Special issue on deep learning for document analysis and recognition. Int. J. Doc. Anal. Recogn. 21(3), 159–160. 10.1007/s10032-018-0310-5 (2018).10.1007/s10032-018-0310-5
30. Breiman L Random Forests Mach. Learn. 2001 45 5 32 10.1023/A:1010933404324
Breiman, L. Random Forests. Mach. Learn. 45, 5–32. 10.1023/A:1010933404324 (2001).10.1023/A:1010933404324
31. Sun X Liu M Sima Z A novel cryptocurrency price trend forecasting model based on LightGBM Financ. Res. Lett. 2020 32 101084 10.1016/j.frl.2018.12.032
Sun, X., Liu, M. & Sima, Z. A novel cryptocurrency price trend forecasting model based on LightGBM. Financ. Res. Lett. 32, 101084. 10.1016/j.frl.2018.12.032 (2020).10.1016/j.frl.2018.12.032
32. Bisong, E. & Bisong, E. Introduction to Scikit-Learn. Building Machine Learning and Deep Learning Models on Google Cloud Platform: A Comprehensive Guide for Beginners 215–229 (2019).
33. Akiba, T. et al. Optuna: A next-generation hyperparameter optimization framework. In Proc. 25th ACM SIGKDD International Conference on Knowledge Discovery & Data Mining 2623–2631 (2019).
34. Lisboa PJG Masulli F Pasi G Yager R Interpretability in machine learning—Principles and practice Fuzzy Logic and Applications 2013 Springer 15 21
Lisboa, P. J. G. Interpretability in machine learning—Principles and practice. In Fuzzy Logic and Applications (eds Masulli, F. et al.) 15–21 (Springer, 2013).
35. Lundberg, S. M. & Lee, S.-I. A unified approach to interpreting model predictions. In Proc. 31st International Conference on Neural Information Processing Systems 4768–4777. http://arxiv.org/abs/1705.07874 (Curran Associates Inc., 2017).
36. de Sousa IP Local interpretable model-agnostic explanations for classification of lymph node metastases Sensors 2019 19 969 10.3390/s19132969 30823554
de Sousa, I. P. et al. Local interpretable model-agnostic explanations for classification of lymph node metastases. Sensors 19, 969. 10.3390/s19132969 (2019).30823554 10.3390/s19132969
37. Nowak AS Radzik T The Shapley value for n-Person games in generalized characteristic function form Games Econom. Behav. 1994 6 150 161 10.1006/game.1994.1008
Nowak, A. S. & Radzik, T. The Shapley value for n-Person games in generalized characteristic function form. Games Econom. Behav. 6, 150–161. 10.1006/game.1994.1008 (1994).10.1006/game.1994.1008
38. Xue J Shen B A novel swarm intelligence optimization approach: Sparrow search algorithm Syst. Sci. Control Eng. 2020 8 22 34 10.1080/21642583.2019.1708830
Xue, J. & Shen, B. A novel swarm intelligence optimization approach: Sparrow search algorithm. Syst. Sci. Control Eng. 8, 22–34. 10.1080/21642583.2019.1708830 (2020).10.1080/21642583.2019.1708830
39. Sibalija TV Particle swarm optimisation in designing parameters of manufacturing processes: A review (2008–2018) Appl. Soft Comput. 2019 84 105743 10.1016/j.asoc.2019.105743
Sibalija, T. V. Particle swarm optimisation in designing parameters of manufacturing processes: A review (2008–2018). Appl. Soft Comput. 84, 105743. 10.1016/j.asoc.2019.105743 (2019).10.1016/j.asoc.2019.105743
40. Dorigo M Blum C Ant colony optimization theory: A survey Theoret. Comput. Sci. 2005 344 243 278 10.1016/j.tcs.2005.05.020
Dorigo, M. & Blum, C. Ant colony optimization theory: A survey. Theoret. Comput. Sci. 344, 243–278. 10.1016/j.tcs.2005.05.020 (2005).10.1016/j.tcs.2005.05.020
41. Mirjalili S Mirjalili SM Lewis A Grey wolf optimizer Adv. Eng. Softw. 2014 69 46 61 10.1016/j.advengsoft.2013.12.007
Mirjalili, S., Mirjalili, S. M. & Lewis, A. Grey wolf optimizer. Adv. Eng. Softw. 69, 46–61. 10.1016/j.advengsoft.2013.12.007 (2014).10.1016/j.advengsoft.2013.12.007
42. Liu L Sun SZ Yu H Yue X Zhang D A modified fuzzy C-means (FCM) clustering algorithm and its application on carbonate fluid identification J. Appl. Geophys. 2016 129 28 35 10.1016/j.jappgeo.2016.03.027
Liu, L., Sun, S. Z., Yu, H., Yue, X. & Zhang, D. A modified fuzzy C-means (FCM) clustering algorithm and its application on carbonate fluid identification. J. Appl. Geophys. 129, 28–35. 10.1016/j.jappgeo.2016.03.027 (2016).10.1016/j.jappgeo.2016.03.027
43. Rudolph G Local convergence rates of simple evolutionary algorithms with Cauchy mutations IEEE Trans. Evol. Comput. 1997 1 249 258 10.1109/4235.687885
Rudolph, G. Local convergence rates of simple evolutionary algorithms with Cauchy mutations. IEEE Trans. Evol. Comput. 1, 249–258. 10.1109/4235.687885 (1997).10.1109/4235.687885
