Application of the Improved POA-RF Model in Predicting the Strength and Energy Absorption Property of a Novel Aseismic Rubber-Concrete Material

The application of aseismic materials in foundation engineering structures is an inevitable trend and research hotspot of earthquake resistance, especially in tunnel engineering. In this study, the pelican optimization algorithm (POA) is improved using the Latin hypercube sampling (LHS) method and the Chaotic mapping (CM) method to optimize the random forest (RF) model for predicting the aseismic performance of a novel aseismic rubber-concrete material. Seventy uniaxial compression tests and seventy impact tests were conducted to quantify this aseismic material performance, i.e., strength and energy absorption properties and four other artificial intelligence models were generated to compare the predictive performance with the proposed hybrid RF models. The performance evaluation results showed that the LHSPOA-RF model has the best prediction performance among all the models for predicting the strength and energy absorption property of this novel aseismic concrete material in both the training and testing phases (R2: 0.9800 and 0.9108, VAF: 98.0005% and 91.0880%, RMSE: 0.7057 and 1.9128, MAE: 0.4461 and 0.7364; R2: 0.9857 and 0.9065, VAF: 98.5909% and 91.3652%, RMSE: 0.5781 and 1.8814, MAE: 0.4233 and 0.9913). In addition, the sensitive analysis results indicated that the rubber and cement are the most important parameters for predicting the strength and energy absorption properties, respectively. Accordingly, the improved POA-RF model not only is proven as an effective method to predict the strength and energy absorption properties of aseismic materials, but also this hybrid model provides a new idea for assessing other aseismic performances in the field of tunnel engineering.


Introduction
The research and development of aseismic concrete materials is of great significance to reducing the damage caused by earthquakes on tunnel structures. To determine the material's performance and decide its application in the field of tunnel engineering, the strength and energy absorption properties are widely used as the main performance evaluation indices [1][2][3][4]. Among numerous novel concrete materials and various component combinations, rubber-concrete material has received much attention for its outstanding capabilities in two performance indices mentioned above [5,6]. Meanwhile, the combination of rubber and traditional concrete material is an effective approach to solving waste rubber pollution [7].
Reviewing the published research, the experiment tests are direct methods to calculate the strength and energy absorption properties of the rubber-concrete material [8]. The strength can be calculated by using the universal testing machine and the Split Hopkinson Pressure Bar (SHPB) device is usually used to measure the energy absorption property [9,10].
Nevertheless, elaborate test sample production and complex test programs are extremely time-consuming [11]. The design and construction of short duration tunnels may be delayed due to time loss, resulting in an inability to resist earthquake shocks.
In recent years, soft computing represented by machine learning (ML) has made remarkable achievements in the prediction of rubber-concrete material performance, e.g., artificial neural networks (ANN) [12,13], support vector machine (SVM) [14,15], backpropagation neural network (BPNN) [16,17], extreme learning machine (ELM) [18,19], multi-layer perceptron (MLP) [20,21] and trees-based models [22,23]. Among the ML models, the random forest (RF) model has an excellent resistance to overfitting and a fitting ability in solving prediction problems [24]. Farooq et al. [25] proposed an optimized RF model using the gene expression programming (GEP) to forecast the strength of high strength concretes. This model has obtained a better performance than the decision tree model and artificial neural networks. Mai et al. [26] utilized the RF model to determine the strength of a novel concrete material, which contains ground granulated blast furnace slag. The prediction results showed that the RF model is a more suitable predictor than time-consuming experiments for engineering. However, the RF model is rarely used to predict the aseismic performance of rubber-concrete materials. Sun et al. [27] used the RF model to predict rubber concrete's strength, which obtained a good predictive performance, the correlation coefficient of which (R) was 0.9596 and the root mean square error (RMSE) was 3.9032. In addition, the energy absorption property of the rubber-concrete material has not been investigated using soft techniques, especially the ML models. Furthermore, numerous metaheuristic optimization algorithms based on swarm intelligence (SI) are used to further improve the predictive ability of ML models for predicting the aseismic concrete material's performance, e.g., particle swarm optimization [28], grey wolf optimizers [29], sparrow search algorithms [11], artificial bee colonies [30] and firefly algorithms [31]. In addition, building more reasonable initial populations of SI algorithms is an effective means to strengthen their optimization abilities [32].
Therefore, an optimized RF model using the pelican optimization algorithm (POA) is proposed to predict the strength and energy absorption properties of a novel aseismic rubber-concrete material in this study. Two functions named Latin hypercube sampling (LHS) and Chaotic mapping (CM) are adopted to improve the initial population of POA. In total, 140 concrete specimens are generated to test their strength (70 specimens) and energy absorption property (70 specimens) for training the prediction models. Four statistical indices are utilized to evaluate the predictive performance of the improved POA-based RF model and other comparison models. Finally, sensitive analysis is responsible for calculating the parameters' importance in the strength and energy absorption property prediction.

Research Significance
As a new functional material with special uses, rubber-concrete material is mainly aimed at weakening the dynamic impact of engineering structures under strong earthquake disturbance. The incorporation of rubber means more impact energy is absorbed and converted into other forms of energy removal, but also leads to the reduction of the rubber concrete's strength. A traditional uniaxial compression test and an impact test based on the SHPB device are the main technical approaches to measure the strength and energy absorption properties of aseismic materials. Although the test results are accurate, sample preparation and equipment operation are extremely time consuming. In the face of urgent engineering needs, the performance estimation of seismic materials without timely estimation is bound to delay the construction period and bring certain economic losses. Moreover, ordinary experimental results and regression analysis cannot accurately grasp the nonlinear relationship between material properties, material proportions, geometric properties and material performance. Therefore, the development of a novel performance estimation method is of great significance for the close fit between novel seismic materials and practical engineering.
Soft computing represented by ML models may be used to learn this nonlinear relationship and make fast and accurate predictions of a material's performance. This study contributes (a) to the development of a novel aseismic concrete material consisting of rubber, sand and concrete and (b) to the accurate prediction of the strength and energy absorption properties of aseismic material, which is significant to get rid of the miscellaneous material performance tests based on large-scale equipment and to measure them in real time in field engineering. The improved POA-RF model is provided to estimate the strength and energy absorption properties of a novel aseismic rubber-concrete material using an actual database collected from a series of performance tests in this study.

Experiment of a Novel Aseismic Concrete Material
To enhance tunnel stability under earthquake resistance, a novel aseismic concrete material doped with rubber and river sand is developed and its performance is assessed. The rubber is responsible for increasing the elasticity and damping of the novel aseismic material, and the river sand is utilized to provide enough strength for it. Although some novel fine aggregates were put forward instead of river sand, e.g., blast furnace slag [33], waste glass [34] and plastic box waste [35], low workability and some substandard physical properties (e.g., voids and abrasion) make them difficult to be adopted in practical projects [36]. Therefore, river sand is still considered the ideal fine aggregate for composing the novel concrete material in this work.
Before testing the strength and energy absorption properties of a novel aseismic rubber-concrete material, the concrete specimens should be prepared according to the test standards, and the production procedure is shown in Table 1. As illustrated in Table 2, there are 27 mix proportions designed for the rubber-concrete material to test the strength (cylindrical shape) and energy absorption properties (disc shape). The rubber mass ratios (RM = M rubber /M RSM ), where M rubber is the rubber mass and M RSM is the rubber-sand mixtures (RSM) mass, are set to 0, 10, 30, 50, 70 and 100% and four different ranges of rubber particle size (RPS) are selected. Thus, the sand mass ratios (SM = M sand /M RSM ), where M sand is the sand mass, are set to 0, 30, 50, 70, 90 and 100%. The cement mass ratios (CM = M cement /M RC ), where M cement is the cement mass and M RC is the specimen mass, are selected as 30, 40, 50 and 60%. Three specimens are prepared for each mix proportion. Because of the influence of test operation, mechanical failure and other factors, the accuracy of some test data is suspicious. As a result, a total of 140 concrete specimens are concreted to carry out experimental tests. Seventy specimens are used to calculate the material strength using the universal testing machine and the rest of specimens are tested using SHPB device to obtain the energy absorption property as shown in Figure 1. To accurately describe the aseismic performance of this novel rubber-concrete material, the uniaxial compressive strength (UCS) and the energy transmission rate (ETR) are adopted to quantify the strength and the energy absorption property, respectively [37]. The UCS and the ETR can be defined using Equations (1) and (2), respectively.
where δ specimen and S specimen represent the UCS and the compression area of the tested specimen. P max is the maximum collapse load. E I and ε I represent the incident energy and the measured strains on the incident bar, respectively. E T and ε T indicate the transmitted energy and the measured strains on the transmission bar, respectively. For the SHPB bar, EM bar is the elasticity modulus, ν is the stress wave velocity and S bar is the crosssectional area. Table 1. A brief description of the preparation procedure for novel concrete specimen.

Procedures Description
Step 1-Dosing Mixing the rubber and river sand with the cement in pre-designed proportions and thoroughly stirring.
Step 2-Concreting Stirring the mixture for five minutes and quickly pouring it into the mold.
Step 3-Demolding After 24 h, separating the specimen and polishing it to the specified specification.
Step 5-Testing After 28 days, testing 140 specimens in the laboratory. It should be noted that each specimen is independent and non-repetitive, such as composition content (rubber particle with different size, river sand and cement), physical properties (mass and density) and size (diameter and length). The detailed information of the considered parameters is listed in Tables 3 and 4. Nevertheless, not all parameters are favorable for predicting the UCS and ETR of the novel concrete material. The results of parameter correlation analysis are shown in Figure 2. As illustrated in Figure 2a, the correlation between rubber, river sand, cement and the specimen mass is significantly higher than that between other parameters and is close to 1, respectively. Thus, the specimen mass (M) should be removed to reduce the calculation time for predicting the UCS. On the other hand, the correlation between the M and specimen density (r) and most other parameters is not low, especially the rubber (R), river sand (S) and cement (C) for predicting the ETR (see Figure 2b). In addition, the R and S are negatively correlated, and the correlation values are high in both the UCS prediction and the ETR prediction data sets, but the reason for this result is that the rubber and river sand are taken as important components of this novel material in the form of combination. For example, concrete with low rubber content may lead to insufficient elasticity and require more river sand to increase its strength. Therefore, neither the rubber nor river sand can be deleted in this study.
Step 1-Dosing Mixing the rubber and river sand with the cement in pre-designed proportions and thoroughly stirring.
Step 2-Concreting Stirring the mixture for five minutes and quickly pouring it into the mold.
Step 3-Demolding After 24 h, separating the specimen and polishing it to the specified specification.
Step 5-Testing After 28 days, testing 140 specimens in the laboratory.

Random Forest
Random forest (RF) is generally regarded as an excellent integrated machine learning model for solving classification and prediction problems. The significant advantage of this model is the data utilization and performance evaluation mechanism, i.e., bootstrap resampling technology is responsible for randomly selecting most data from the original database (the ratio is two to three) to form the decision trees (DTs) and the rest of data is utilized to establish the test set for evaluating the DTs performance; the final performance of the RF model is determined by using the average value of the predicted values of all DTs.
As illustrated in Figure 3, the number of trees (Nt) and the random features (Maxdepth) are the main hyperparameters that affect its prediction performance. Although the increase of the Nt will not cause the overfitting of the model, it is difficult to obtain satisfactory predictive performance through the time-consuming manual debugging of hyperparameters' combinations [38,39].

Random Forest
Random forest (RF) is generally regarded as an excellent integrated machine learning model for solving classification and prediction problems. The significant advantage of this model is the data utilization and performance evaluation mechanism, i.e., bootstrap resampling technology is responsible for randomly selecting most data from the original database (the ratio is two to three) to form the decision trees (DTs) and the rest of data is utilized to establish the test set for evaluating the DTs performance; the final performance of the RF model is determined by using the average value of the predicted values of all DTs.
As illustrated in Figure 3, the number of trees (N t ) and the random features (Maxdepth) are the main hyperparameters that affect its prediction performance. Although the increase of the N t will not cause the overfitting of the model, it is difficult to obtain satisfactory predictive performance through the time-consuming manual debugging of hyperparameters' combinations [38,39].

Pelican Optimization Algorithm
Trojovský and Dehghani [40] proposed a novel metaheuristic optimization algorithmsbased swarm intelligence to solve the optimization problem, that is, the pelican optimization algorithm (POA). The development of the algorithm was inspired by the fact that pelicans often hunt prey such as fish in groups. Pelicans' hunting behavior is full of wisdom. For example, pelicans determine the position of prey in advance, and then approach the prey quickly and complete the hunting behavior with a swotting posture when the distance between them is close to 10-20 m [41]. The detailed framework of the POA can be represented using following steps: (a) Population initialization: Pelicans generally search for prey within a certain range of search space, so the initial position of each pelican in the population is random and can be expressed as: where P i represents the initial position of the i-th pelican. I indicates the maximum number of pelicans. S min and S max are the minimum and the maximum boundaries of the search space, respectively. rand · () represents a random number in (0, 1). (b) Exploration phase: At this stage, the main goal of pelicans is to find and determine the prey position and change their positions in preparation for an attack. Therefore, the position of each pelican is updated using Equation (4): where P 1 i and P p represent the updated position of the pelican and the prey position in the exploration phase, respectively. f p and f i indicate the objective functions of the prey and the pelican, respectively. U is also a random number (1 or 2).
(c) Exploitation phase: Once the pelican is in a good position and launches an attack, i.e., flying across the water and pushing the fish up into its throat pouch. This strategy can be mathematically expressed as: where P 2 i represents the updated position of the pelican in the exploitation phase. z is a random number (0 or 2). T is the maximum number of iterations, and t represents the current iteration.

Optimization Methods
For the metaheuristic optimization algorithms-based SI, the incompatibility of population initialization results in low precision and an easy to fall into local minimum problem [42]. Numerous researchers have utilized various strategies to optimize the population initialization of SI algorithms, e.g., opposition-based learning [43,44], Gaussian mutation [45,46], normal distribution [47] and multi-subgroup [48]. In addition, the chaotic mapping (CM) and the Latin hypercube sampling (LHS) have also widely been used to improve the performance of the SI algorithms by adjusting the population initialization [49,50]. The former is characterized by traversal and random, while the latter has significant space-filling impact and convergence features for solving the population initialization problem [51,52]. Therefore, CM and the LHS are adopted to improve the POA for predicting the UCS and the ETR of a novel aseismic rubber-concrete material in this study.

Chaotic Mapping Method (CM)
The aim of the CM method is to use mapping functions to generate a more diverse population of the POA. The logistic mapping function is a popular method to achieve this goal [53,54], which can be expressed using Equation (6).
where Log i and Log i+1 represent the i-th and i+1-th chaotic sequence of the logistic mapping function, respectively. λ is a constant. Thus, Equation (3) can be rewritten as: where P log i represents the new initial position of the i-th pelican. η indicates an iteration factor. f log is the logistic mapping function.

Latin Hypercube Sampling Method (LHS)
The function of the LHS is to stratify the initial positions in three-dimensional space, and then conduct random sampling and disrupt the order of each sample to obtain a more stable and diverse population. Therefore, the new population can be described as follows: where P LHS i represents the new initial position of the i-th pelican. k and K indicate the current dimension of pelican and the maximum dimension of the search space, respectively.

Pelican Optimization Algorithm
Trojovský and Dehghani [40] proposed a novel metaheuristic optimization algorithms-based swarm intelligence to solve the optimization problem, that is, the pelican optimization algorithm (POA). The development of the algorithm was inspired by the fact that pelicans often hunt prey such as fish in groups. Pelicans' hunting behavior is full of wisdom. For example, pelicans determine the position of prey in advance, and then approach the prey quickly and complete the hunting behavior with a swotting posture when the distance between them is close to 10-20 m [41]. The detailed framework of the POA can be represented using following steps: (a) Population initialization: Pelicans generally search for prey within a certain range of search space, so the initial position of each pelican in the population is random and can be expressed as: () rand ⋅ represents a random number in (0, 1).
(b) Exploration phase: At this stage, the main goal of pelicans is to find and determine the prey position and change their positions in preparation for an attack. Therefore, the position of each pelican is updated using Equation (4):

A Novel Combination of the IPOA and RF Model
Two improved optimization algorithms consisting of Latin hypercube sampling-POA (LHSPOA) and chaotic mapping-POA (CMPOA) are used to optimize the RF model for predicting the UCS and the ETR of a novel aseismic rubber-concrete material in this study. As illustrated in Figure 4, the framework of using the proposed models to predict the aseismic performance of the novel concrete material can be organed in three parts: (a) Building the database-as mentioned in Section 3, 70 samples were used to predict the UCS of this novel aseismic concrete material, and another 70 samples were used to predict the ETR. Whether the UCS or the ETR prediction, the ratio of the train set to the test set is four to one. It should be noted that all data is normalized to [−1, 1] to eliminate the data discrepancies; (b) Developing prediction models-the POA and the improved POA were combined with an RF model to generate different prediction models, i.e., POA-RF, LHSPOA-RF and CMPOA-RF. The population size is a key factor affecting the optimization algorithm's performance [55,56], which is set as 20, 40, 60, 80 and 100 to explore the potential of the model. The iteration is 200 and the RMSE is considered as a fitness function to determine the optimal population size. For the RF model, the ranges of N t and Maxdepth are (1, 100) and (1, 10), respectively; (c) Performance evaluation-four statistical indices were used to evaluate the predictive performance of the proposed hybrid RF models. model to generate different prediction models, i.e., POA-RF, LHSPOA-RF and CMPOA-RF. The population size is a key factor affecting the optimization algorithm's performance [55,56], which is set as 20, 40, 60, 80 and 100 to explore the potential of the model. The iteration is 200 and the RMSE is considered as a fitness function to determine the optimal population size. For the RF model, the ranges of Nt and Maxdepth are (1, 100) and (1, 10), respectively; (c) Performance evaluation-four statistical indices were used to evaluate the predictive performance of the proposed hybrid RF models.

Performance Evaluation
In this study, four indices were used to evaluate the performance of the proposed models for predicting the UCS and the ETR of the novel aseismic concrete material. The determination coefficient (R 2 ) is also known as goodness of fit, which reflects the interpretation of the predicted value to the measured value from the fitting perspective. If R 2 is equal to 1, that means the prediction model is perfect. The variance accounted for (VAF) is often used to evaluate the degree to which the prediction model can explain the variance of the considered data. The function of the root mean square error (RMSE) is to evaluate the model performance by measuring the error between the predicted value and the measured value. In addition, the mean absolute error (MAE) can further reflect the real situation of error. These statistical indices have been considered to verify the performance of different prediction models for solving the regression problem [57][58][59][60][61][62].
where T indicates the number of used samples. m t and p t represent the t-th measured and predicted values, respectively. m is the average of the measured values.  Figure 5b,c. However, it should be noted that the fitness values of the best LHSPOA-RF model and the best CMPOA-RF model are lower than the POA-RF model (see Table 5). As illustrated in Figure 5d, the POA-RF model with a population size of 40 achieved a lower fitness value than other models for predicting the ETR of the novel aseismic concrete material. Compared with the best POA-RF model, the most satisfying population sizes of the LHSPOA-RF model and the CMPOA-RF model are 80 and 40, respectively (see Figure 5e,f). As a result, the optimal hyperparameter combinations of RF models based on the three optimization algorithms for predicting the UCS and the ETR are listed in Table 5.

Performance Evaluation for the UCS and the ETR Prediction
After establishing the optimal hybrid RF models, each model was first utilized to predict the UCS and the ETR of this novel aseismic concrete material based on the train set. The performance indices of all the models are listed in Table 6. As can be observed in this table, the LHSPOA-RF model is not only the best model with the highest accuracy (R 2 : 0.9800, RMSE: 0.7057, MAE: 0.4461 and VAF: 98.0005%) for predicting the UCS, but also obtained better performance indices (R 2 : 0.9108, RMSE: 1.9128, MAE: 0.7364 and VAF: 91.0880%) than other models in the ETR prediction. After this model, the performance indices of the CMPOA-RF model are also better than the POA-RF for predicting the UCS and ETR. In addition, the rank score is a useful tool to evaluate the model performance by calculating the score of all the performance indices. For example, the score values corresponding to the highest values of R 2 and RMSE in the same set are 4 and 1, respectively. As demonstrated in Table 6, the LHSPOA-RF model has achieved the highest values of ranking scores (12 and 12) for predicting the UCS and the ETR in the training phase. Figure 6 illustrates the regression diagrams of the proposed hybrid RF models using the train set. In each diagram, the blue line represents the perfect prediction function y = x, i.e., the predicted value is equal to the measured value. Therefore, the more points that locate on the blue line or close to it, the more superior the model's performance compared with others. Of course, the regression lines with 10% have similar functions for evaluating the models' performance. As can be seen in Figure 6b,e, the LHSPOA-RF model not only obtained the most data points close to or located at the blue line in the UCS prediction, but also showed the same excellent prediction performance in the ETR prediction resulting in the best performance indices. After the LHSPOA-RF, the CMPOA-RF model has achieved a more satisfying predictive performance than the POA-RF model resulting in more data points within the 10 % lines for the UCS prediction as shown in Figure 6a,c, and obtained the same performance for predicting the ETR (see Figure 6d,f).   However, the performance evaluation of the model in the training phase is not enough to prove its real prediction performance, especially as it is applied to predict unknown data. Therefore, the test set is responsible for verifying the performance of the proposed models for predicting the UCS and the ETR. Table 7 shows the performance evaluation results of all the models in the testing phase. Compared with the performance results obtained using the ETR dataset in the training phase, the performance indices of each model are somewhat worse, but not by much. The result indicates that the adverse overfitting phenomenon does not appear in the proposed models, which suggests that the proposed model is accepted and may be applied to practical engineering. Furthermore, the LHSPOA-RF model is still the best prediction model for forecasting the UCS and the ETR resulting in the highest values of R 2 (0.9857 and 0.9065) and VAF (98.5909% and 91.3652%), and the lowest values of RMSE (0.5781 and 1.8814) and MAE (0.4233 and 0.9913). After the LHSPOA-RF model, the CMPOA-RF model has better performance indices than the POA-RF model for predicting the UCS and the ETR. The ranking score results of the proposed models for predicting the UCS and ETR in the testing phase are also shown Table 7. Based on these results, the order of the three hybrid RF models is the LHSPOA-RF model, the CMPOA-RF model and the POA-RF model in the UCS prediction. The same ranking score results of each model is obtained by using the ETR test data.
The regression diagrams of the proposed hybrid RF models using the test set are shown in Figure 7. For the UCS prediction, the LHSPOA-RF model has achieved the most satisfactory performance by means of the most points located on the perfect line or closer to it than other models as illustrated in Figure 7b. As shown in Figure 7a,c, the predictive performance of CMPOA-RF is still better than the POA-RF by means of fewer points out of the 10% line. On the other hand, the predictive performance of the proposed models for forecasting the ETR is obviously worse than that of the UCS as shown in Figure 7d,e. However, the regression results showed that the LHSPOA-RF model and the CMPOA-RF model still have better performance than the POA-RF model resulting in the superior performance indices and more effective points close to the perfect line, especially the LHSPOA-RF model. As a results, two improved POA-RF models (i.e., the LHSPOA-RF model and the CMPOA-RF) are considered to be the superior predictors, rather than the initial POA-RF model, in predicting the UCS and the ETR of the novel aseismic rubber-concrete material in this study.
Moreover, four common ML models named the BPNN, support vector regression (SVR), ELM and kernel-extreme learning machine (KELM) are developed to predict the UCS and the ETR of the novel aseismic rubber-concrete material and compare the predictive performances with the proposed models. For the BPNN model, the hyperparameters are the number of hidden layers (N h ) and the number of neurons in each hidden layer (N e ). The penalty parameter (P c ) and the RBF kernel parameter (k 1 ) are the main parameters affecting the SVR model performance. The ELM model is a special ANN model with a single hidden layer, thus the number of neurons in this layer (N n ) can control the performance of the ELM model. The KELM model has similar hyperparameters with the SVR model, i.e., regularization coefficient (R c ) and the RBF kernel parameter (k 2 ). After developing these four ML models, the performance indices of each model with the best hyperparameter combination are shown in Table 8. The performance comparison results of the proposed models and four common ML models in the UCS and the ETR prediction are shown in Figure 8. As can be seen in these pictures, the performance indices of the LHSPOA-RF model and the CMPOA-RF model are obviously better than that of the four ML models, i.e., higher values of R 2 and VAF, and lower values of RMSE and MAE. Meanwhile, it is obvious that the LHSPOA-RF model is the best model to predict the UCS and the ETR of the novel aseismic rubber-concrete material in this study.
performance of CMPOA-RF is still better than the POA-RF by means of fewer points out of the 10% line. On the other hand, the predictive performance of the proposed models for forecasting the ETR is obviously worse than that of the UCS as shown in Figure 7d,e. However, the regression results showed that the LHSPOA-RF model and the CMPOA-RF model still have better performance than the POA-RF model resulting in the superior performance indices and more effective points close to the perfect line, especially the LHSPOA-RF model. As a results, two improved POA-RF models (i.e., the LHSPOA-RF model and the CMPOA-RF) are considered to be the superior predictors, rather than the initial POA-RF model, in predicting the UCS and the ETR of the novel aseismic rubberconcrete material in this study. As illustrated in Figure 9, the Taylor diagram is used to make a clear comparison between the proposed models and the four common ML models. This diagram consists of the RMSE, the St. D and the R calculated from the predicted values. In particular, the measured value is fixed on the horizontal axis since its R is equal to 1 and the RMSE is equal to 0. Assuming that the prediction model achieves the best performance, its position on the Taylor diagram is the closest to the measured value. Based on this criterion, the predictive performance of the proposed models is obviously superior to the other four ML models. The LHSPOA-RF is the best prediction model for predicting the UCS and the ETR of this novel aseismic material.  As illustrated in Figure 9, the Taylor diagram is used to make a clear comparison between the proposed models and the four common ML models. This diagram consists of the RMSE, the St. D and the R calculated from the predicted values. In particular, the measured value is fixed on the horizontal axis since its R is equal to 1 and the RMSE is equal to 0. Assuming that the prediction model achieves the best performance, its position on the Taylor diagram is the closest to the measured value. Based on this criterion, the predictive performance of the proposed models is obviously superior to the other four ML models. The LHSPOA-RF is the best prediction model for predicting the UCS and the ETR of this novel aseismic material.

Sensitively Analysis
Although the best prediction model has been determined through performance evaluation, the input parameters' importance for predicting the UCS and the ETR is unknown. Therefore, the PAWN method proposed by Pianosi and Wagener [63,64] is utilized to conduct sensitive analysis based on the LHSPOA-RF model. The parameter importance results are shown in Figure 10. As demonstrated in Figure 10a, rubber is the most important input parameter with a score of 0.6143 for predicting the UCS of the novel aseismic rubberconcrete material. After this parameter, cement (0.5714), rubber particle size (0.4857) and river sand (0.4009) have higher importance scores than specimen diameter (0.3571), specimen length (0.3250) and specimen density (0.2571). For the ETR prediction (see Figure 10b), the cement has the highest importance score of 0.5841 among all the input parameters. The values of rubber and specimen diameter were significantly higher than the other parameters, i.e., rubber (0.4256), specimen diameter (0.3122), specimen length (0.1521), river sand (0.1365) and rubber particle size (0.1123). Figure 8. The performance comparison of the proposed models and four common ML models in the UCS and the ETR prediction.
As illustrated in Figure 9, the Taylor diagram is used to make a clear comparison between the proposed models and the four common ML models. This diagram consists of the RMSE, the St. D and the R calculated from the predicted values. In particular, the measured value is fixed on the horizontal axis since its R is equal to 1 and the RMSE is equal to 0. Assuming that the prediction model achieves the best performance, its position on the Taylor diagram is the closest to the measured value. Based on this criterion, the predictive performance of the proposed models is obviously superior to the other four ML models. The LHSPOA-RF is the best prediction model for predicting the UCS and the ETR of this novel aseismic material.

Sensitively Analysis
Although the best prediction model has been determined through performance evaluation, the input parameters' importance for predicting the UCS and the ETR is unknown. Therefore, the PAWN method proposed by Pianosi and Wagener [63,64] is utilized to conduct sensitive analysis based on the LHSPOA-RF model. The parameter importance results are shown in Figure 10. As demonstrated in Figure 10a, rubber is the most important input parameter with a score of 0.6143 for predicting the UCS of the novel aseismic rubberconcrete material. After this parameter, cement (0.5714), rubber particle size (0.4857) and river sand (0.4009) have higher importance scores than specimen diameter (0.3571), specimen length (0.3250) and specimen density (0.2571). For the ETR prediction (see Figure  10b), the cement has the highest importance score of 0.5841 among all the input parameters. The values of rubber and specimen diameter were significantly higher than the other parameters, i.e., rubber (0.4256), specimen diameter (0.3122), specimen length (0.1521), river sand (0.1365) and rubber particle size (0.1123).

Conclusions
The strength and energy absorption properties of aseismic concrete materials have a significance effect on tunnel stability when a tunnel is shocked by an earthquake and other dynamic impacts. To clearly investigate the above two aseismic performance indices of a novel aseismic concrete material, experimental tests and soft computation techniques are combined in this study. Based on the uniaxial compression test and impact tests, an improved pelican optimization algorithm (POA) using the Latin hypercube sampling (LHS) method and the Chaotic mapping (CM) method was proposed to optimize the random forest (RF) model for predicting the aseismic performance of this novel material. The prediction results demonstrated that the LHSPOA-RF model and CMPOA-RF model have better predictive performances than the general POA-RF model for predicting the UCS and the ETR, especially the LHAPOA-RF model with a more satisfying performance accuracy of R 2 , RMSE, MAE and VAF (UCS: 0.9857, 0.5781, 0.4233 and 98.5909 %; ETR: 0.9065, 1.8814, 0.9913 and 91.3652 %). Furthermore, the comparison results of four common ML models and the proposed hybrid RF models showed that the LHSPOA-RF model is the best prediction performance among all models.
In addition, rubber plays an important role in strength prediction and cement is the most important parameter for predicting the energy absorption property.
This work indicates that the application of the intelligence models to the aseismic material performance is an effective and more simple method that does not lose accuracy compared with traditional experiment tests in the laboratory. However, the limitation of

Conclusions
The strength and energy absorption properties of aseismic concrete materials have a significance effect on tunnel stability when a tunnel is shocked by an earthquake and other dynamic impacts. To clearly investigate the above two aseismic performance indices of a novel aseismic concrete material, experimental tests and soft computation techniques are combined in this study. Based on the uniaxial compression test and impact tests, an improved pelican optimization algorithm (POA) using the Latin hypercube sampling (LHS) method and the Chaotic mapping (CM) method was proposed to optimize the random forest (RF) model for predicting the aseismic performance of this novel material. The prediction results demonstrated that the LHSPOA-RF model and CMPOA-RF model have better predictive performances than the general POA-RF model for predicting the UCS and the ETR, especially the LHAPOA-RF model with a more satisfying performance accuracy of R 2 , RMSE, MAE and VAF (UCS: 0.9857, 0.5781, 0.4233 and 98.5909%; ETR: 0.9065, 1.8814, 0.9913 and 91.3652%). Furthermore, the comparison results of four common ML models and the proposed hybrid RF models showed that the LHSPOA-RF model is the best prediction performance among all models.
In addition, rubber plays an important role in strength prediction and cement is the most important parameter for predicting the energy absorption property.
This work indicates that the application of the intelligence models to the aseismic material performance is an effective and more simple method that does not lose accuracy compared with traditional experiment tests in the laboratory. However, the limitation of this combined method is the number of samples in the used database. More and more tested or monitored data should be added into the intelligence models' development to improve the prediction accuracy in the further research.