Open Access

ARTICLE

# Seismic Liquefaction Resistance Based on Strain Energy Concept Considering Fine Content Value Effect and Performance Parametric Sensitivity Analysis

1 School of Civil Engineering and Geomatics, Southwest Petroleum University, Chengdu, 610500, China

2 Key Laboratory of Geological Hazards on Three Gorges Reservoir Area, Ministry of Education, China Three Gorges University, Yichang, 443002, China

3 College of Civil Engineering & Architecture, China Three Gorges University, Yichang, 443002, China

4 Department of Civil Engineering, Faculty of Engineering, International Islamic University Malaysia, Jalan Gombak, Selangor, 50728, Malaysia

5 Department of Civil Engineering, University of Engineering and Technology Peshawar (Bannu Campus), Bannu, 28100, Pakistan

6 Central Queensland University, Queensland, 4740, Australia

* Corresponding Author: Jilei Hu. Email:

(This article belongs to this Special Issue: Computational Intelligent Systems for Solving Complex Engineering Problems: Principles and Applications)

*Computer Modeling in Engineering & Sciences* **2023**, *135*(1), 733-754. https://doi.org/10.32604/cmes.2022.022207

**Received** 26 February 2022; **Accepted** 20 May 2022; **Issue published** 29 September 2022

## Abstract

Liquefaction is one of the most destructive phenomena caused by earthquakes, which has been studied in the issues of potential, triggering and hazard analysis. The strain energy approach is a common method to investigate liquefaction potential. In this study, two Artificial Neural Network (ANN) models were developed to estimate the liquefaction resistance of sandy soil based on the capacity strain energy concept (*W*) by using laboratory test data. A large database was collected from the literature. One group of the dataset was utilized for validating the process in order to prevent overtraining the presented model. To investigate the complex influence of fine content (

*FC*) on liquefaction resistance, according to previous studies, the second database was arranged by samples with

*FC*of less than 28% and was used to train the second ANN model. Then, two presented ANN models in this study, in addition to four extra available models, were applied to an additional 20 new samples for comparing their results to show the capability and accuracy of the presented models herein. Furthermore, a parametric sensitivity analysis was performed through Monte Carlo Simulation (MCS) to evaluate the effects of parameters and their uncertainties on the liquefaction resistance of soils. According to the results, the developed models provide a higher accuracy prediction performance than the previously published models. The sensitivity analysis illustrated that the uncertainties of grading parameters significantly affect the liquefaction resistance of soils.

## Keywords

Nomenclature

uw | Water pressure |

σ’c | Effective confining pressure |

FC% | Fine content in percent |

Cu | Coefficient of uniformity |

D50 | Mean grain size (mm) |

Cc | coefficient of curvature |

W | liquefaction resistance of sandy soil based on capacity strain energy concept |

ANN | Artificial neural network |

GP | Genetic programming |

LGP | Linear genetic programming |

MEP | Multi expression programming |

ANFIS | Neuro-fuzzy interface system |

MARS | multivariate adaptive regression splines |

MCS | Monte Carlo simulation |

FOSM | first order second moment |

PEM | point estimation |

CSR | Cyclic stress ratio |

CRR | Cyclic strength ratio |

E | Unit energy |

γ | Shear strain amplitude |

ν | Standard deviation |

δW | Increment of energy/volume |

NME | Normalized maximum energy |

Δu | Excess pore water pressure |

δ3 | Lateral stress |

τ | Shear stress |

Γ | Shear strain amplitude |

R2 | Coefficient of determination |

ydj | Target output |

yj | Predicted output |

di | Individual sample points indexed with i |

x0 | Mean sample size |

RMSE | Root mean square error |

MAE | Mean absolute error |

COV | Coefficient of variation |

When saturated sand is subjected to an earthquake, because of the rapid vibrations, drainage is prevented and the tendency towards volume reduction, causes the transfer of the effective overburden stress to the pore water until excess pore water pressure becomes equal to the initial effective overburden stress; after which liquefaction happens. The most commonly reported liquefaction manifests in saturated loose or medium sandy soil have been observed during the most massive earthquakes worldwide. The 1964 magnitude 9.2 earthquake in Alaska and the magnitude 7.6 earthquake in Niigata of the same year, prompted extensive research on this phenomenon. Soil liquefaction has also been observed in recent earthquakes in China [1], Japan [2], Indonesia [3] and the USA [4,5].

Three main methods have been employed in relevant studies. The first one is the stress-based method, which was introduced by Seed et al. [6]. And other researchers performed research to develop models using in-situ tests [7–9], laboratory tests [10,11] and numerical simulation [12–19]. Additionally, the strain-based method was first introduced by Dobry et al. [12]. They assumed under cyclic loading approximately 0.01% for threshold strain and initial water pressure (uw) increasing. After that, the shear strain was compared with 0.01 according to Serikawa et al. [2] or 0.02 [3]. Next, water pressure was estimated through experimental graphs. In the end, this water pressure was compared to confining stress to predict the triggering of liquefaction. In other words, depending on the following condition, liquefaction may or may not occur:

uw > σv0 liquefaction occurs.

uw < σv0 liquefaction does not occur.

where σv0 is initial vertical effective stress.

The coupled numerical models, since 1975, have been presented [13–17] based on Biot’s theory [18–20] which was the clarification of effective stress concept and coupled phases interaction between solid porous materials and fluid. Recently, some numerical simulation was also performed by researchers [21–28]. The fourth method includes strain energy-based methods developed by applying seismic energy dissipated in the soil [21–30]. This method has been applied in three main procedures by researchers which are using histories of site exploration liquefied [29–31], and laboratory test results [23,27,29–40] and Arias intensity-based models [32,41]. To evaluate the potential of liquefaction in energy concept method, the capacity strain energy (W) value of the soil is required to be estimated to compare with the energy transferred to the soil by the earthquake loads. Since the energy dissipated by mechanisms (e.g., cohesion and frictional mechanisms) cannot be easily discerned from laboratory and field data, the energy dissipated by frictional mechanisms is estimated by the total dissipated energy as the frictional mechanisms are expected to be dominant growing interest in earthquake engineering. In addition, the total amount of dissipated energy to the liquefaction point should be relatively independent of the sequence for increasing the load. On the contrary, the viscous mechanisms of energy dissipation can be considered for the low increase of strains where the rate of energy dissipated by this mechanism is directly related to the sequence used. In order to enhance the conventional load, the dissipated energy is expected to be greater than the liquefaction point, which is either independent or increases due to the load sequence used.

Based on laboratory test results six input parameters including effective confining pressure (σ’c) kPa, initial relative density (Dr)%, FC%, coefficient of uniformity (Cu), mean grain size (D50) (mm) and coefficient of curvature (Cc), have been identified and confirmed as the most influential factors in modeling liquefaction to estimate liquefaction resistance of sandy soil based on capacity strain energy concept [23,27,29–40]. Clearly, permeability of the soil is considered implicitly in soil properties parameters of Dr, Cu, D50 and Cc.

These studies, except Cabalar et al. [37], extracted Cc in their final correlation due to its limited range values in the datasets and hence, its limited effect. Furthermore, they included all ranges of the parameters in their models without special consideration to their value.

Further, Maurer et al. [42] analyzed 7,000 case histories from Canterbury Earthquakes in 2010–2011 and concluded when soils contain a high value of FC, assessment of liquefaction is less reliable. Zhang et al. via laboratory tests showed that liquefaction potential is closely related to FC [38]. Liu et al. [43] performed some experimental tests on marine sediments and proposed a critical value for FC to evaluate liquefaction resistance. Additionally, Tao [44] defined the limit value of 28% for estimating the liquefaction resistance, and through laboratory test results showed liquefaction resistance becomes more dependent on Dr when FC is higher than 28%, rather than FC value. While, in all presented models there has not been consideration to different effect of FC in different range.

Regarding the evaluation of liquefaction strength by applying the strain energy approach, several models have been developed using artificial neural network (ANN) [35], genetic programming (GP) [45,46], multi expression programming (MEP) [46], neuro-fuzzy Interface system (ANFIS) [37], and multivariate adaptive regression splines (MARS) [38]. Although the importance of the validating phase has been indicated by many researchers [47–49], in all studies, data division was performed randomly in two groups of testing and training phases, without considering the statistical characteristics of the data. Also, no validating phase has been performed in order to avoid overtraining of the models.

Due to the uncertainty of geotechnical problems, particularly, liquefaction phenomena, some studies, such as Bayesian methods, have been performed to develop probabilistic forms and reliability analysis to evaluate the potential of liquefaction [50–56]. Furthermore, artificial intelligence [57–63] and Monte Carlo simulation (MCS) which is a classic approach to assess risk in quantitative analysis, has recently been used in engineering and sciences [64–71] and also in the evaluation of liquefaction potential [72–74]. Jha et al. [73] presented the probability of liquefaction due to factor of safety using FOSM method, an advanced first-order second-moment (FOSM), Hasofer–Lind reliability method, a point estimation (PEM), and an MCS method. They presented a new combined method using both FOSM and PEM to find the cyclic stress ratio (CSR) and cyclic strength ratio (CRR) statistically. They showed that the factor of safety measured by the combined method is similar to the PEM and MCS methods. They also indicated FOSM, PEM, and MCS methods present nearly the same probabilities of liquefaction by considering input variability. Using the jointly distributed random variables method and using the data from triaxial test results, Johari et al. [74] presented a reliability assessment of liquefaction and compared the results with the Monte Carlo simulation. The results exhibited close probability density functions of the safety factor applying both methods.

In this study, to investigate the complex effect of FC on liquefaction resistance of soil in terms of the unit energy, two datasets were arranged. The first dataset was collected from the literature as the largest and likely most complete dataset employed by researchers covering a large range of parameters. Due to the complicated influence of FC on W and the spares attention to this parameter in developing earlier models, in the second database, according to Tao [44] only samples with FC values less than 28% were selected. Two new ANN models were developed based on these two datasets. A multilayer perceptron network with a backpropagation algorithm was constructed and the samples were divided into three groups, including a validation set to avoid overtraining. These sample groups were formed with similar statistics certificates, and avoided random division, according to Tables 1–4 and Tables 6–9 to enhance the accuracy and capability of trained networks.

This study investigates the effects of all parameters, including CC while also paying special consideration to the influence of FC in different range according to its critical value in seismic soil liquefaction assessment. To achieve this goal, two different ANN models, one using the entire dataset and the other using the samples with FC value of less than critical value, were developed to compare their predictions for validating and choosing the best one. In development of the models the validation phase was also performed to eliminate overtraining of the models. The data division was performed considering statistics characteristics of the variables instead of performing randomly to increase the accuracy of the trained models. Furthermore, to the best of the author’s knowledge, there has been no previous due attention to the uncertainties of parameters to predict W, which was the motivation behind performing the sensitivity analysis via MCS simulation to investigate the effect of uncertainty and the mean value of parameters on liquefaction resistance. Because of the numerous samples required by MCS and the relative data scarcity, in this study, the MC simulation was performed based on an ANN model developed in this study.

2 Methods Based on Laboratory test Results

Figueroa et al. [75] developed two equations to evaluate unit energy (E) in a cyclic triaxial test. Alkhatib [76] introduced ER to measure liquefaction resistance. ER is a ratio of the energy computed by area under the stress-strain hysteresis loop to the initial effective confining stress, and through conducting laboratory cyclic triaxial tests. The presented model is as below.

Extensive research has been performed at Case Western Reserve University on energy-based evaluation of liquefaction [31,38–41,75,77–79]. In all procedures, Wu which is the area of the stress-strain hysteresis loops up to the initial liquefaction point was used to define liquefaction resistance. Parameter of δW was introduced first time by Figueroa et al. [31,75,77]. They conducted 27 torsional shear tests on a Reid Bedford sand sample at different shear strain amplitudes and confining pressures and developed a model.

Liang et al. [13,41] conducted 74 liquefaction torsional shear tests on Reid Bedford sand, Lower San Fernando Dam (LSFD) silty sand, and Lapis Luster Dried sand (LSI-30) through random loading. From the test results, they performed a regression analysis and presented an equation to estimate δW.

Kusky [78] developed two equations according to 27 strain-controlled torsional triaxial tests, which were conducted on samples of Reid Bedford.

Rokoff [79] conducted some cyclic torsional shear tests on Nevada sand to investigate the influence of particle size distribution on δW. The regression was limited to special soil properties and geology related to the samples of Nevada region, which contains Cu and Cc as below:

In addition, Figueroa et al. [40] confirmed that Cu and Cc affect δW more than

Baziar et al. [35] collected a large dataset from performed shear, cyclic triaxial, and torsional shear laboratory test results, which contained 284 samples from the literature. They developed two Artificial neural network (ANN) models to obtain a correlation between input parameters and Log (W). The first developed ANN model contained six input parameters (σ’c, Dr%, FC%, Cu, D50, Cc) while the second model was developed by eliminating the parameter Cc. They subsequently demonstrated that FC has the highest effect on W by carrying out a sensitivity analysis. By adding a new dataset to Baziar et al. [35] with the same parameters and applying multigene Genetic Programming, Baziar et al. [45] developed an equation to measure W and then used case histories earthquake data plus laboratory test data to validate and present the accuracy of their model. They concluded that the value of W has a complicated relationship with FC. Alavi et al. [46] presented three equations to estimate Log (W) through applying MEP, GP, and MEP and with the same database and parameters as Baziar et al. [35], as mentioned in the appandix and Table A1. In addition, they conducted sensitivity analysis and confirmed that W is more affected by Dr and σ’cthan other parameters. Zhang et al. [38] collected 302 samples for their database, which contained six cyclic simple shear, 18 centrifuges, six cyclic simple shear, and 217 cyclic tests. They developed a MARS model with the same five input parameters with Cavallaro et al. [10] to evaluate Log (W). They validated their model using 22 centrifuge test results conducted by Dief [81].

Artificial Neural Network (ANN) is defined as brain model systems, which are collections of mathematical models containing cells (here called neurons) interconnected by links. The goal of ANN is to utilize a training process to learn a nonlinear multiplex relation between parameters to approximate a target (output). Training is the process of calculating weights, which indicates the strength of the links between neurons. There are several neural network types proposed, but feed forward neural networks are the most capable and commonly applied type. Among all classes of neural network topologies, Hornik et al. [82] demonstrated multilayer perceptrons (MLP), which are supervised networks, with the best capacity and ability to approximate any function with high accuracy. These include three types of layers: an input layer which distributes the input data and contains one neuron for each input variable, one or more hidden layers which perform non-linear transformations, additions, and multiplications; and an output layer for estimated final results, which contains a number of neurons equal number of targets meant to be approximated by the ANN model. The backpropagation algorithm is one of the most commonly used algorithms for training ANN’s. Here, connection weights are updated by estimating error and distributing it through the layers of neurons. It contains two steps that are iterated to obtain a pre-specified tolerance range of the output. In the first step, the network generates an output, and in the second step, the estimated error at the output layer is distributed to the hidden layers and then to the input layer to modify the weights. Each neuron’s error is calculated by:

Further, the overall neurons’ output is estimated by:

The correlation coefficient (R) is the most common and capable tool to test the performance on networks given by:

Network samples are commonly separated into two subsets randomly; the first one is the training set to train the network by adjusting the weights of the network, and the second one is the testing set. Testing samples are not used in the training step and are applied to assess the performance of the trained network. A new sample set, called validation set, should be selected to prevent overtraining the network which occurs when the accuracy and the correlation coefficient increase, but the accuracy and the correlation coefficient of the validation samples set decreases. When overfitting starts, training should be stopped.

3.2 Monte Carlo Simulation and Uncertainties

The Monte Carlo method was introduced first during research on the atomic bomb in the beginning of the 1940s. The main idea is using random samples of inputs or parameters to discover the response of a complex process or system through observing the fraction of numbers. It involves three main steps:

1. Generating random input samples called scenarios.

2. Simulating each scenario to explore the response.

3. Evaluating outputs of all simulations to estimate statistics certificates and properties such as minimum and maximum values, mean value, and distribution function for each variable.

Conservative values of loads and soil properties cannot be reliably assumed due to inaccuracy in measurements and models, as well as the inherent variability in the systems under consideration. In geotechnical soil properties, these uncertainties are categorized into two groups: aleatory and epistemic [83]. Aleatory uncertainties are defined as natural randomness such as spatial variability of soil properties and are related to inherent randomness, which cannot be reduced by adding new data and information. Epistemic uncertainties on the other hand, are caused by a shortage of data, information, and measurement procedures or model error as well as non-standard equipment, laboratory instruments, and random testing effects [84]. Reliability approaches provide a formal way to deal with uncertainties and quantify them. Monte Carlo Simulation (MCS) conducts risk assessment by providing a probability distribution for any variables due to their uncertainties to estimate the possible outcomes of an uncertain phenomenon.

3.3 MCS Based ANNs Response Surface for Sensitivity Analysis

The main idea of the response surface method is a computational calculation reduction. In the classic form, the surface was approximated through an equivalent function such as polynomial form [85], which is not capable of modeling high nonlinear phenomena [86] such as liquefaction. Thus, in this study, the response surface which belongs to the ANN trained model is applied. The procedure of MCS-based ANNs response surface for sensitivity analysis is described in the flowchart of Fig. 1.

4 Models Presented in this Study

In this study, two different databases were arranged to train two ANN models to investigate the complex influence of FC on liquefaction resistance. According to previous research [23,27,29–34,36–39,41], six parameters of

Tao [44] studied the effect of FC value on the liquefaction resistance by considering the void ratio. He demonstrated that Dr becomes more effective when the FC grows above 28%. He declared that there is no clear correlation between the entire range of FC and liquefaction resistance. Therefore, in this study, the second database was arranged by collecting just samples with FC value lower than 28% to train an ANN model and hence samples with FC higher than 28% were eliminated from the dataset. Consequently, the second dataset contains 309 samples, which were divided into three parts, considering to have similar statistics certificates to achieve a more capable and accurate model. Around 15% of samples (44 samples) were selected for testing, equal portion and numbers for validation, and 221 samples for training the model. The characteristics and statistical factors of the second database are summarized in Tables 6–9. For example, the mean value of FC in all datasets (including training group, validation group, and testing group) were 7.8, 7.9 and 8.1, respectively. Note that because of deleting 94 samples included FC value of larger than 28%, the second dataset and its subsets would provide different statistical factors than the first dataset. The characteristics of the second ANN model are presented in Table 10. It can be seen that the value of R for all groups in dataset division (i.e., all data, training data, validating data, and testing data) is greater than 90%, which indicates a high power fitting model.

4.2 Comparison of the Predicted Value of W Using the ANN Models and Available Models

In this section, the results of the ANN models are compared to the other four well-established models i.e., GP, LGP, MEP and MARS [38,46] to evaluate their capability. For more details about these models, refer to the Appendix.

To achieve this goal, 20 laboratory test results from Dief [81], performed on Nevada sand and Reid Bedford sand, considering the range of applied database are selected. Note that these 20 samples were not used in the database to construct the two ANN models developed in this study. The results predicted by these four models and two presented ANN models are presented in Table 11.

To compare the capability and accuracy of all six models, three criteria of root mean square error (RMSE), mean absolute error (MAE), and R2 are estimated and summarized in Table 12. As can be seen, two ANN models show higher agreement and less error between predicted and measured results in comparison with other available models.

Fig. 2 illustrates that all predicted values through ANN models are close to the measured values for log (W). Two ANN models developed in this study predict log (W) with high accuracy, as presented in Table 12. The first and second ANN models predicted log (W) with R2 of 0.77 and 0.83, respectively that are higher than the value of extra four models. In addition, the first ANN with RMSE and MAE values of 0.13 and 0.11, respectively, and the second ANN (referred to as ANN28 herein) with RMSE and MAE values of 0.1 and 0.09, respectively, demonstrate the highest precision.

Based on the illustrated figures and Table 12, the two presented ANN models are the most accurate and capable models for predicting log (W) and between them, the second model, which contains a dataset with a limited FC value of less than 28%, indicates more accuracy. Note that the ANN28 model was developed based on fewer samples due to eliminating samples with FC values larger than 28%.

As mentioned in Section 4, most geotechnical parameters, soil properties, and applied loads are uncertain. To deal with these uncertainties, reliability methods have been used to quantify and capture these uncertainties. In this study, MC simulation was applied to perform sensitivity analysis and investigate the influence of parameters and their uncertainties by changing their mean values and coefficient of variations (COV) or standard deviation (ν). Monte Carlo simulation requires a large number of samples to present a reliable response. Providing such a large number of samples is costly and time-consuming. Therefore, to overcome this shortage, the second ANN model was applied to provide a response surface for MCS to be able to conduct sensitivity analysis. Phoon et al. [84] suggested a mean COV of 19% for sand with Dr ranging from 11% to 36%; Therefore, in this study to evaluate this parameter’s effect on log (W), it was supposed to have a mean COV equal to 20% with minimum and maximum value of 10% and 30%, respectively. Subsequently,

Furthermore, given the fact that with a small value of ν, the distribution function supposition error is insignificant, normal distribution was assigned to all variables [90,91]. All statistical properties of parameters are summarized in Tables 13 and 14. It should be mentioned that during parametric sensitivity analysis of each variable, the other five variables fixed in their mean value and mean COV value, without changing, then analysis was conducted.

Additionally, in order to conduct a sensitivity analysis through MC simulation, a definition of correlation coefficient (ρ) is required. By considering the independency of all six input parameters, ρ among all parameters is supposed to be 0. The value of 2.9 was chosen for reliability analysis to assess the cumulative probability density function. As can be observed in Fig. 3, the probability of log (W) larger than 2.9 is illustrated as a function of the parameters and their uncertainties.

By considering Fig. 3, which plots parameters vs. probability of if log (W) be higher than 2.9, can be seen, there is a slight increase (i.e., 15%) in log (W) > 2.9 is observed for σ’c from 44 to 250 and then, it grows dramatically to 75% at σ’c beyond 250. Upon increasing COV from 5% to 10% and then 15%, probability grows two times by 1.5%. The probability rises slightly from 4% to 60% in the range of Dr from 5.44 onwards. It experiences an impressive rise to 60% as Dr increases to 71.5% also, by growing uncertainty as COV from 0.1 to 0.2, and then 0.3 in the critical range of Dr between 35% to 70%, the probability shows two increases of 3%. During the growth of the FC value until 28%, the probability shows slight growth from 21% to approximately 24.5% and it experiences a negligible increase while the COV changes. Furthermore, the probability of log (W) > 2.9 illustrates a falling range from 100% to 0% during the range of Cu in this study. Next to that, by increasing COV from 0.1 to 0.2, and subsequently 0.3 in the critical values between 13 to 16, the probability augments around 7% every time. There was a steady climb of around 19% in the probability in the range of the D50 in this study. In addition, by increasing any 0.1 in COV, from 0.1 to 0.2 and then 0.3, the probability rises negligibly less than 1% for log (W) > 2.9. Whereas, by any 10% increase in COV of Cc results show around 2.5% growth in the probability in a sense that the probability goes up around 58% in the range of Cc from 0.74 to 10.89.

In this study, ANN was used to develop models to estimate the liquefaction resistance of sandy soil based on the capacity strain energy concept and using laboratory test data. The validating phase was performed, in addition to the testing and training phase, to avoid overtraining the model and increasing the model’s capability. An extensive database was collected from literature, including triaxial, simple shear, torsional, and centrifuge test results. ANNs are powerful tools for developing models that can take into account the complexity and non-linearity of the liquefaction issue. To inspect the complicated influence of FC on liquefaction resistance of soil, according to research results presented by Tao [44], two ANN models were developed. The first model was developed using a complete dataset, while the second one was based on the samples by FC less than 28%. The accuracy and capability of the presented models were demonstrated by comparing their predicted values for log (W) with four other available well-known equations. To conduct this comparison, 20 liquefaction test results from Nevada sand and Reid Bedford sand [41], which were independent of the two applied datasets for training the models, were considered. Finally, to investigate the effect of uncertainty in geotechnical parameters, a sensitivity analysis was performed using MCS based on the response surface provided by the second presented ANN model, which showed higher accuracy. The results of sensitivity analysis were illustrated through some graphs to indicate the correlation between the variables and their uncertainties with the liquefaction resistance of the soil in order to capacity energy. The limitation of the present study includes its application in the issue of strain energy, not in the other methods such as stress-based or numerical methods.

In conclusion, this study has demonstrated:

1. Artificial neural network (ANN) is a powerful tool to assess liquefaction in soil with high non-linearity. Adding validation phase and performing data division by considering the statistical aspects, instead of random division, provides significant precision on the model.

2. The second ANN model (considering samples with FC less than 28%) is able to predict log (W) with higher accuracy. As it includes a smaller number of samples in the dataset in comparison with the first ANN model, it is evident that different FC values provide a different effect on liquefaction resistance.

3. The parameter of Cc significantly affected W and should be considered to predict the W value.

4. The uncertainty of parameters had a considerable impact on liquefaction resistance. As a result, performing probabilistic frameworks and models are suggested by the authors instead of deterministic models to consider and quantitate these uncertainties’ effects.

Funding Statement: This work is supported by the Scientific Innovation Group for Youths of Sichuan Province under Grant No. 2019JDTD0017.

Conflicts of Interest: The authors declare that they have no conflicts of interest to report regarding the present study.

## References

- Shao, Z. F., Zhong, J. H., Howell, J., Hao, B., & Luan, X. W. (2020). Liquefaction structures induced by the M5.7 earthquake on May 28, 2018 in Songyuan, Jilin province, NE China and research implication.
*Journal of Palaeogeography*,*9(1)*, 3. [Google Scholar] [CrossRef] - Serikawa, Y., Miyajima, M., Yoshida, M., & Matsuno, K. (2019). Inclination of houses induced by liquefaction in the 2018 Hokkaido Iburi-Tobu earthquake, Japan.
*Geoenvironmental Disasters*,*6(1)*, 14. [Google Scholar] [CrossRef] - Jalil, A., Fathani, T. F., Satyarno, I., & Wilopo, W. (2021). Liquefaction in palu: The cause of massive mudflows.
*Geoenvironmental Disasters*,*8(1)*, 21. [Google Scholar] [CrossRef] - Zimmaro, P., Nweke, C., Hernandez, J. L., Hudson, K. S., & Hudson, M. B. (2020). Liquefaction and related ground failure from July 2019 ridgecrest earthquake sequence.
*Bulletin of the Seismological Society of America*,*110(4)*, 1549-1566. [Google Scholar] [CrossRef] - Mavroulis, S., Lekkas, E., & Carydis, P. (2021). Liquefaction phenomena induced by the 26 November 2019, Mw = 6.4 Durrës (Albania) earthquake and liquefaction susceptibility assessment in the affected area.
*Geosciences*,*11(5)*, 215. [Google Scholar] [CrossRef] - Seed, H. B., Idriss, I. M. (1971). Simplified procedure for evaluating soil liquefaction potential. Journal of the Soil Mechanics and Foundations Division, 97(9), 1249–1273. DOI 10.1061/JSFEAQ.0001662. [CrossRef]
- Idriss, I. M., Boulanger, R. W. (2004). Semi-empirical procedures for evaluating liquefaction potential during earthquakes. Proceedings of the 11th International Conference on Soil Dynamics and Earthquake Engineering, and 3rd International Conference on Earthquake Geotechnical Engineering, pp. 32–56. Berkeley, California, USA.
- Idriss, I. M., & Boulanger, R. W. (2006). Semi-empirical procedures for evaluating liquefaction potential during earthquakes.
*Soil Dynamics and Earthquake Engineering*,*26(2–4)*, 115-130. [Google Scholar] [CrossRef] - Idriss, M. (1999). An update to the Seed-Idriss simplified procedure for evaluating liquefaction potential. TRB Workshop on New Approaches to Liquefaction, Washington DC: Federal Highway Administration.
- Cavallaro, A., Capilleri, P. P., & Grasso, S. (2018). Site characterization by dynamic in situ and laboratory tests for liquefaction potential evaluation during Emilia Romagna earthquake.
*Geosciences*,*8(7)*, 242. [Google Scholar] [CrossRef] - Castelli, F., Cavallaro, A., Grasso, S., & Lentini, V. (2019). Undrained cyclic laboratory behavior of sandy soils geosciences.
*Geosciences*,*9(12)*, 512. [Google Scholar] [CrossRef] - Dobry, R., Ladd, R. S., Yokel, F. Y., Chang, R., Powell, D. J. (1982). Prediction of pore water pressure buildup and liquefaction of sands during earthquakes by the cyclic strain method. National Bureau of Standards Building Science Series, 138–168. DOI 10.6028/NBS.BSS.138. [CrossRef]
- Liang, L. (1995). Development of an energy method for evaluating the liquefaction potential of a soil deposit (Ph.D. Thesis), Case Western Reserv University, Ohio. http://rave.ohiolink.edu/etdc/view?acc_num=case1058541489.
- Elgamal, A. W., Parra, E., Yang, Z., Dobry, R., Zeghal, M. (1999). Liquefaction constitutive model. In: Physics and mechanics of soil liquefaction (1st Edition). London: Taylor & Francis Group. DOI 10.1201/9780203743317. [CrossRef]
- Elgamal, A., Yang, Z., & Parra, E. (2002). Computational modeling of cyclic mobility and post-liquefaction site response.
*Soil Dynamics and Earthquake Engineering*,*22(4)*, 259-271. [Google Scholar] [CrossRef] - Dafalias, Y. F., & Manzari, M. T. (2004). Simple plasticity sand model accounting for fabric change effects.
*Journal of Engineering Mechanics*,*130(6)*, 622-634. [Google Scholar] [CrossRef] - Gao, Z., & Zhao, J. (2015). Constitutive modeling of anisotropic sand behavior in monotonic and cyclic loading.
*Journal of Engineering Mechanics*,*141(8)*, 04015017. [Google Scholar] [CrossRef] - Biot, M. A. (1956). Theory of propagation of elastic waves in a fluid-Saturated porous solid. II. higher frequency range.
*The Journal of the Acoustical Society of America*,*28(2)*, 179. [Google Scholar] [CrossRef] - Biot, M. A. (1962). Mechanics of deformation and acoustic propagation in porous media.
*Journal of Applied Physics*,*33(4)*, 1482-1498. [Google Scholar] [CrossRef] - Biot, M. A. (1941). General theory of three-dimensional consolidation.
*Journal of Applied Physics*,*12(2)*, 155-164. [Google Scholar] [CrossRef] - Dashti, S., & Bray, J. D. (2013). Numerical simulation of building response on liquefiable sand.
*Journal of Geotechnical and Geoenvironmental Engineering*,*139(8)*, 1235-1249. [Google Scholar] [CrossRef] - Dashti, S., & Karimi, Z. (2017). Ground motion intensity measures to evaluate I: The liquefaction hazard in the vicinity of shallow-founded structures.
*Earthquake Spectra*,*33(1)*, 241-276. [Google Scholar] [CrossRef] - Fotopoulou, S., Karafagka, S., & Pitilakis, K. (2018). Vulnerability assessment of low-code reinforced concrete frame buildings subjected to liquefaction-induced differential displacements.
*Soil Dynamics and Earthquake Engineering*,*110*, 173-184. [Google Scholar] [CrossRef] - Forcellini, D. (2019). Numerical simulations of liquefaction on an ordinary building during Italian (20 May 2012) earthquake.
*Bulletin of Earthquake Engineering*,*17(9)*, 4797-4823. [Google Scholar] [CrossRef] - Forcellini, D. (2020). Soil-structure interaction analyses of shallow-founded structures on a potential-liquefiable soil deposit.
*Soil Dynamics and Earthquake Engineering*,*133*, 106108. [Google Scholar] [CrossRef] - Petridis, C., & Pitilakis, D. (2020). Fragility curve modifiers for reinforced concrete dual buildings, including nonlinear site effects and soil–structure interaction.
*Earthquake Spectra*,*36(4)*, 1930-1951. [Google Scholar] [CrossRef] - Forcellini, D. (2021). Seismic fragility for a masonry-infilled RC (MIRC) building subjected to liquefaction.
*Applied Sciences*,*11(13)*, 6117. [Google Scholar] [CrossRef] - Karafagka, S., Fotopoulou, S., & Pitilakis, D. (2021). Fragility assessment of non-ductile RC frame buildings exposed to combined ground shaking and soil liquefaction considering SSI.
*Engineering Structure*,*229*, 111629. [Google Scholar] [CrossRef] - Davis, R. O. (1982). Energy dissipation and seismic liquefaction in sands.
*Earthquake Engineering & Structural Dynamics*,*10(1)*, 59-68. [Google Scholar] [CrossRef] - Law, K. T., Cao, Y. L., & He, G. N. (1990). An energy approach for assessing seismic liquefaction potential.
*Canadian Geotechnical Journal*,*27(3)*, 320-329. [Google Scholar] [CrossRef] - Figueroa, J. L., Saada, A. S., Liang, L., & Dahisaria, N. M. (1994). Evaluation of soil liquefaction by energy principles.
*Journal of Geotechnical Engineering*,*120(9)*, 1554-1569. [Google Scholar] [CrossRef] - Kayen, R. E., & Mitchell, J. K. (1997). Assessment of liquefaction potential during earthquakes by arias intensity.
*Journal of Geotechnical and Geoenvironmental Engineering*,*123(12)*, 1162-1174. [Google Scholar] [CrossRef] - Chen, Y. R., Hsieh, S. C., Chen, J. W., & Shih, C. C. (2005). Energy-based probabilistic evaluation of soil liquefaction.
*Soil Dynamics and Earthquake Engineering*,*25(1)*, 55-68. [Google Scholar] [CrossRef] - Green, R. A., & Terri, G. A. (2005). Number of equivalent cycles concept for liquefaction evaluations.
*Journal of Geotechnical and Geoenvironmental Engineering*,*131(4)*, 477-488. [Google Scholar] [CrossRef] - Baziar, M. H., & Jafarian, Y. (2007). Assessment of liquefaction triggering using strain energy concept and ANN model: Capacity energy.
*Soil Dynamics and Earthquake Engineering*,*27(12)*, 1056-1072. [Google Scholar] [CrossRef] - Okur, D. V., & Ansal, A. (2007). Stiffness degradation of natural fine grained soils during cyclic loading.
*Soil Dynamics and Earthquake Engineering*,*27(9)*, 843-854. [Google Scholar] [CrossRef] - Cabalar, A. F., Cevik, A., Gokceoglu, C. (2012). Some applications of adaptive neuro-fuzzy inference system (ANFIS) in geotechnical engineering. Computers and Geotechnics, 40, 14–33. DOI 10.1016/j.compgeo.2011.09.008. [CrossRef]
- Zhang, W., Goh, A. T. C., Zhang, Y., Chen, Y., & Xiao, Y. (2015). Assessment of soil liquefaction based on capacity energy concept and multivariate adaptive regression splines.
*Engineering Geology*,*188*, 29-37. [Google Scholar] [CrossRef] - Ulmer, K. J. (2019). Development of an energy-based liquefaction evaluation procedure (Ph.D. Thesis). Virginia Polytechnic Institute and State University, Blacksburg, Virginia.
- Figueroa, J. L., Saada, A. S., Rokoff, M. D., Liang, L. (1998). Influence of grain-size characteristics in determining the liquefaction potential o f a soil deposit by the energy method. Proceedings of the International Workshop on the Physics and Mechanics of Soil Liquefaction, pp. 237–245. Baltimore, Maryland, USA.
- Liang, L., Figueroa, J. L., & Saada, A. S. (1995). Liquefaction under random loading, unit energy approach.
*Journal of Geotechnical Engineering*,*121(11)*, 776-781. [Google Scholar] [CrossRef] - Maurer, B. W., Green, R. A., Cubrinovski, M., & Bradley, B. A. (2015). Fines-content effects on liquefaction hazard evaluation for infrastructure in Christchurch, New Zealand.
*Soil Dynamics and Earthquake Engineering*,*76*, 58-68. [Google Scholar] [CrossRef] - Liu, B., & Jeng, D. S. (2016). Laboratory study for influence of clay content (CC) on wave-induced liquefaction in marine sediments.
*Marine Georesources & Geotechnology*,*34(3)*, 280-292. [Google Scholar] [CrossRef] - Tao, M. (2003). Case history verification of the energy method to determine the liquefaction potential of soil deposits (Ph.D. Thesis). Case Western Reserve University, Cleveland, Ohio.
- Baziar, M. H., Jafarian, Y., Shahnazari, H., Movahed, V., & Tutunchian, M. A. (2011). Prediction of strain energy-based liquefaction resistance of sand–silt mixtures, an evolutionary approach.
*Computers & Geosciences*,*37(11)*, 1883-1893. [Google Scholar] [CrossRef] - Alavi, A. H., & Gandomi, A. H. (2012). Energy-based numerical models for assessment of soil liquefaction.
*Geoscience Frontiers*,*3(4)*, 541-555. [Google Scholar] [CrossRef] - Kohavi, R. A. (1995). Study of cross-validation and bootstrap for accuracy estimation and model selection. Proceedings of the 14th International Joint Conference on Artificial Intelligence, pp. 1137–1143. Montreal, Quebec, Canada.
- Zeng, X., & Martinez, T. R. (2000). Distribution-balanced stratified cross-validation for accuracy estimation.
*Journal of Experimental & Theoretical Artificial Intelligence*,*12(1)*, 1-12. [Google Scholar] [CrossRef] - Karaci, A. (2019). Estimating the properties of ground-waste-brick mortars using DNN and ANN.
*Computer Modeling in Engineering & Sciences*,*118(1)*, 207-228. [Google Scholar] [CrossRef] - Hwang, J. H., Yang, C. W., & Juang, D. S. A. (2004). Practical reliability-based method for assessing soil liquefaction potential.
*Soil Dynamics and Earthquake Engineering*,*24(9–10)*, 761-770. [Google Scholar] [CrossRef] - Juang, C. H., Fang, S. Y., & Khor, E. H. (2006). First-order reliability method for probabilistic liquefaction triggering analysis using CPT.
*Journal of Geotechnical and Geoenvironmental Engineering*,*132(3)*, 337-350. [Google Scholar] [CrossRef] - Juang, C. H., Ching, J., Lou, Z. (2013). Assessing SPT-based probabilistic models for liquefaction potential evaluation: A 10-year update. Georisk: Assessment and Management of Risk for Engineered Systems and Geohazards, 7(3), 137–150. DOI 10.1080/17499518.2013.778117. [CrossRef]
- Jha, S. K., & Suzuki, K. (2009). Liquefaction potential index considering parameter uncertainties.
*Engineering Geology*,*107(1–2)*, 55-60. [Google Scholar] [CrossRef] - Huang, H. W., Zhang, J., & Zhang, L. M. (2012). Bayesian network for characterizing model uncertainty of liquefaction potential evaluation models.
*KSCE Journal of Civil Engineering*,*16(5)*, 714-722. [Google Scholar] [CrossRef] - Hu, J., & Liu, H. (2019). Bayesian network models for probabilistic evaluation of earthquake-induced liquefaction based on CPT and Vs databases.
*Engineering Geology*,*254*, 76-88. [Google Scholar] [CrossRef] - Hu, J., & Liu, H. (2019). Identification of ground motion intensity measure and its application for predicting soil liquefaction potential based on the Bayesian network method.
*Engineering Geology*,*248*, 34-49. [Google Scholar] [CrossRef] - Gorodnichev, E. E., Kondratiev, K. A., Kuzovlev, A. I., & Rogozhin, D. B. (2020). Propagation and depolarization of a short pulse of light in sea water.
*Journal of Marine Science and Engineering*,*8(5)*, 371. [Google Scholar] [CrossRef] - Alyousef, R., Ali, B., Mohammed, A., Kurda, R., & Alabduljabbar, H. (2021). Evaluation of mechanical and permeability characteristics of microfiber-reinforced recycled aggregate concrete with different potential waste mineral admixtures.
*Materials*,*14(20)*, 5933. [Google Scholar] [CrossRef] - Piro, N. S., Salih, A., Hamad, S. M., & Kurda, R. (2021). Comprehensive multiscale techniques to estimate the compressive strength of concrete incorporated with carbon nanotubes at various curing times and mix proportions.
*Journal of Materials Research and Technology*,*15*, 6506-6527. [Google Scholar] [CrossRef] - Asteris, P. G., Lourenço, P. B., Roussis, P. C., Adami, C. E., & Armaghani, D. J. (2022). Revealing the nature of metakaolin-based concrete materials using artificial intelligence techniques.
*Construction and Building Materials*,*322*, 126500. [Google Scholar] [CrossRef] - Asteris, P. G., Mamou, A., Hajihassani, M., Hasanipanah, M., & Koopialipoor, M. (2021). Soft computing based closed form equations correlating L and N-type Schmidt hammer rebound numbers of rocks.
*Transportation Geotechnics*,*29*, 100588. [Google Scholar] [CrossRef] - Jahed Armaghani, D., Harandizadeh, H., Momeni, E. (2021). Load carrying capacity assessment of thin-walled foundations, an ANFIS–PNN model optimized by genetic algorithm. Engineering with Computers, 1–23. DOI 10.1007/s00366-021-01380-0. [CrossRef]
- Parsajoo, M., Armaghani, D. J., Mohammed, A. S., Khari, M., & Jahandari, S. (2021). Tensile strength prediction of rock material using non-destructive tests, a comparative intelligent study.
*Transportation Geotechnics*,*31*, 100652. [Google Scholar] [CrossRef] - Armaghani, D. J., Harandizadeh, ., Momeni, E., Maiziri, H., & Zhou, J. (2022). An optimized system of GMDH-ANFIS predictive model by ICA for estimating pile bearing capacity.
*Artificial Intelligence Review*,*55(3)*, 2313-2350. [Google Scholar] [CrossRef] - Asteris, P. G., Lemonis, M. E., Le, T. T., & Tsavdaridis, K. D. (2021). Evaluation of the ultimate eccentric load of rectangular CFSTs using advanced neural network modeling.
*Engineering Structures*,*248*, 113297. [Google Scholar] [CrossRef] - Li, C., Zhou, J., Armaghani, D. J., & Li, X. (2021). Stability analysis of underground mine hard rock pillars via combination of finite difference methods, neural networks, and Monte Carlo simulation techniques.
*Underground Space*,*6(4)*, 379-395. [Google Scholar] [CrossRef] - Lin, W., Su, C., & Tang, Y. (2020). Explicit time-domain approach for random vibration analysis of jacket platforms subjected to wave loads.
*Journal of Marine Science and Engineering*,*8(12)*, 1001. [Google Scholar] [CrossRef] - Xiong, M., & Huang, Y. (2020). Static and dynamic reliability analysis of laterally loaded pile using probability density function method.
*Journal of Marine Science and Engineering*,*8(12)*, 994. [Google Scholar] [CrossRef] - Dou, D., Zeng, Z., Yu, W., Zeng, M., & Men, W. (2021). In-situ seawater gamma spectrometry with LaBr3 detector at a nuclear power plant outlet.
*Journal of Marine Science and Engineering*,*9(7)*, 721. [Google Scholar] [CrossRef] - Lin, W., & Su, C. (2021). An efficient Monte-Carlo simulation for the dynamic reliability analysis of jacket platforms subjected to random wave loads.
*Journal of Marine Science and Engineering*,*9(4)*, 380. [Google Scholar] [CrossRef] - Shojaei Barjouei, A., & Naseri, M. A. (2021). Comparative study of statistical techniques for prediction of meteorological and oceanographic conditions, an application in sea spray icing.
*Journal of Marine Science and Engineering*,*9(5)*, 539. [Google Scholar] [CrossRef] - Hwang, J. H., Chen, C. H., Juang, C. H. (2005). Liquefaction hazard analysis, a fully probabilistic method. Geo-Frontiers Congress 2005, pp. 1–15. Austin, Texas, USA. DOI 10.1061/40779(158)22. [CrossRef]
- Jha, S. K., & Suzuki, K. (2008). Reliability analysis of soil liquefaction based on standard penetration test.
*Computers and Geotechnics*,*36(4)*, 589-596. [Google Scholar] [CrossRef] - Johari, A., Javadi, A., Makiabadi, M. H., Khodaparast, A. R. (2012). Reliability assessment of liquefaction potential using the jointly distributed random variables method. Soil Dynamics and Earthquake Engineering, 38, 81–87. DOI 10.1016/j.soildyn.2012.01.017. [CrossRef]
- Figueroa, J. L., Saada, A. S., Liang, L., & Dahisaria, N. M. (1996). Closure to Evaluation of liquefaction by energy principles by J. Ludwig Figueroa, Adel S. Saada, Liqun Liang, and Nitin M. Dahisaria.
*Journal of Geotechnical Engineering*,*122(3)*, 244-244. [Google Scholar] [CrossRef] - Alkahatib, M. (1994). Liquefaction assessment by strain energy aprroach (Ph.D. Thesis). Wayne State University, Detroit, Michigan.
- Figuerao, J. L., Dahisaria, N. (1991). An energy approach in defining soil liquefaction. International Conferences on Recent Advances in Geotechnical Earthquake Engineering and Soil Dynamics, pp. 407–410. Rolla, UK, University of Missouri.
- Kusky, P. J. (1996). Influence of loading rate on the unit energy required for liquefaction (M.S.Thesis). Case Western Reserve University, Cleveland, Ohio.
- Rokoff, M. D. (1999). The influence of grain-size characteristics in determining the liquefaction potential of a soil deposit by the energy method (M.S.Thesis). Case Western Reserve University, Cleveland, Ohio.
- Wallin, M. S. (2000). Evaluation of normalized pore water pressure vs. accumulated unit energy relationships for determining liquefaction potential in soils (Ph.D. Thesis). Department of Civil Engineering, Case Western Reserve University, Cleveland, USA.
- Dief, H. M. (2000). Evaluation of soil liquefaction by energy principles using centrifuge modeling (Ph.D. Thesis). Virginia Polytechnic Institute and State University, Blacksburg.
- Hornik, K., Stinchcombe, M., & White, H. (1989). Multilayer feedforward networks are universal approximators.
*Neural Networks*,*2(5)*, 359-366. [Google Scholar] [CrossRef] - Lacasse, S., & Nadim, F. (1997). Uncertainties in characterising soil properties.
*Publikasjon-Norges Geotekniske Institute*,*201*, 49-75. [Google Scholar] - Phoon, K. K., & Kulhawy, F. H. (1999). Evaluation of geotechnical property variability.
*Canadian Geotechnical Journal*,*36(4)*, 625-639. [Google Scholar] [CrossRef] - Box, G. E. P., & Behnken, D. W. (1960). Some New three level designs for the study of quantitative variables.
*Technometrics*,*2(4)*, 455-475. [Google Scholar] [CrossRef] - Azeiteiro, R. J. N., Coelho, P. A. L. F., Taborda, D. M. G., & Grazina, J. C. D. (2017). Energy-based evaluation of liquefaction potential under non-uniform cyclic loading.
*Soil Dynamics and Earthquake Engineering*,*92*, 650-665. [Google Scholar] [CrossRef] - Green, R. A. (2001). Energy-based evaluation and remediation of liquefiable soils (Ph.D. Thesis). Virginia Polytechnic Institute and State University, USA.
- Arulanandan, K., & Scott, R. F. (1993). Project VELACS 2014, control test results.
*Journal of Geotechnical Engineering*,*119(8)*, 1276-1392. [Google Scholar] [CrossRef] - Kanagalingam, T. (2006). Liquefaction resistance of granular mixes based on contact densityand energy considerations and environmental engineering (Ph.D. Thesis). The State University of New York at Buffalo, USA.
- Juang, C. H., Rosowsky, D. V., & Tang, W. H. (1999). Reliability-based method for assessing liquefaction potential of soils.
*Journal of Geotechnical and Geoenvironmental Engineering*,*125(8)*, 684-689. [Google Scholar] [CrossRef] - Lumb, P. (1966). The variability of natural soils.
*Canadian Geotechnical Journal*,*3(2)*, 74-97. [Google Scholar] [CrossRef]

Alavi et al. [46] developed three equations using genetic programming (GP), linear genetic programming (LGP), and multi expression programming (MEP) to evaluate the strength of soil liquefaction according to the capacity energy as below:

GP model:

LGP model:

MEP model:

The normalized variables used in these three equations are defined as below:

Zhang et al. [38] developed an equation using MARS as below:

Table A1 presents all the coefficients required in Eq. (A5).

## Cite This Article

Pirhadi, N., Wan, X., Lu, J., Hu, J., Ahmad, M. et al. (2023). Seismic Liquefaction Resistance Based on Strain Energy Concept Considering Fine Content Value Effect and Performance Parametric Sensitivity Analysis.*CMES-Computer Modeling in Engineering & Sciences, 135(1)*, 733–754.