LEAP-UCD-2017 Simulation Team Fugro

Fugro participated in the Liquefaction Experiment and Analysis Projects (LEAP) by performing numerical simulations using two different constitutive models implemented in the software FLAC. Fugro developed a calibration framework based on soil-specific laboratory test data considering multiple elements of dynamic response such as liquefaction triggering criteria (i.e., number of cycles to a specified strain and pore pressure ratio threshold) as well as post-triggering strain accumulation rate. The calibrated model parameters were subsequently used to obtain “Type B (blind)” predictions of the centrifuge experiments with the opportunity to refine after the centrifuge results were provided (Type C simulations). Overall, Fugro’s blind predictions captured the centrifuge test responses well with small refinements needed during Type C simulations. Estimated deformations were within a factor of about two compared to the observed. The overall good comparison provides confidence in the proposed calibration framework, which can be implemented and used for different sand types and project conditions.


Introduction
The LEAP 2017 simulation exercise consisted of the following four stages: (1) constitutive model calibration, (2) Type-B (blind) predictions, (3) Type-C simulations, and (4) sensitivity analyses. A series of stress-controlled cyclic triaxial laboratory tests were available (El Ghoraiby et al. 2017, 2019 and used in the calibration process. Next, the model parameters obtained during model calibration were used to analyze the centrifuge experiments, providing the base excitation and the soil density without any knowledge of the actual results (Type-B predictions). The simulation teams were subsequently provided with the centrifuge experiment results and simulations were refined, if necessary, to obtain better agreement (Type-C simulations).
Last, a series of sensitivity analyses were performed to assess the influence of soil relative density (D R ), motion intensity, and motion high frequency content on the simulation results.

Introduction
Fugro calibrated two constitutive models, PM4SAND and UBCSAND, to the available cyclic stress-controlled triaxial laboratory tests on Ottawa F65 sand (El Ghoraiby et al. 2017, 2019. The cyclic triaxial tests were performed for three different soil densities. The three groups of specimens tested had void ratios of approximately 0.585, 0.542, and 0.515 corresponding to relative densities of about 65%, 80%, and 90%, respectively, based on a maximum void ratio (e max ) of 0.74 and a minimum void ratio (e min ) of 0.49 per Vasko as tabulated in Table 6 of the El Ghoraiby et al. (2017) report. At the time of this calibration process, no other test data was available.
Constitutive model calibration was performed considering multiple elements of response such as number of cycles to specified strain and pore pressure ratio thresholds as well as post-triggering strain accumulation rate.

PM4SAND Constitutive Model
PM4SAND Version 3 constitutive model developed by Boulanger and Ziotopoulou (2015) was initially calibrated. Stress-controlled plane strain compression (PSC) test simulations were performed using PM4SAND and calibrated against the available stress-controlled triaxial test data. The soil elements were initially confined with vertical and horizontal stresses of one atmosphere, considering coefficient of earth pressures at rest of one (K o ¼ 1). Plane strain compression loading conditions were stress controlled by imposing a velocity until a specific stress ratio was reached and then the sign of velocity was changed.
The model parameters are grouped into two categories: a primary and a secondary set of parameters that may be modified from their default values in special circumstances. The three primary input parameters are the sand's apparent relative density D R , the shear modulus coefficient G o , and the contraction rate parameter h po . The contraction rate parameter h po was used to calibrate the liquefaction triggering resistance, and the shear modulus parameter G o was mainly used to calibrated the strain rate accumulation. PM4SAND calibrated model parameters are tabulated in Table 27.1. A detailed description of the primary and secondary input parameters and default values is provided in Tables 4.1 and 4.2 in Boulanger and Ziotopoulou (2015).
Plane strain compression tests were simulated for the relative densities of 65%, 80%, and 90% for various cyclic stress ratios (CSR) in order to compare the plane strain simulation results with the triaxial test results. Each plane strain compression simulation was compared to a triaxial test of a soil sample with similar relative density and stress ratio (CSR) in terms of liquefaction triggering and strain accumulation rate. Liquefaction strength curves were plotted for both the actual triaxial tests and the simulated plane strain compression tests in terms of numbers of cycles to 98% excess pore pressure ratio, and 2 and 5% double amplitude axial strain. Figure 27.1 presents the observed (from the tests-red) and simulated (with calibrated soil properties-blue) liquefaction strength curve for the case of sand with D R of 65%. Figure 27.2 presents the observed and simulated response in terms of (a) deviatoric stress versus mean effective stress, (b) deviatoric stress versus axial strain, (c) mean effective stress versus axial strain, and (d) pore pressure ratio versus number of cycles for all the plane strain compression simulations and superimposed with the respective triaxial tests for comparison purposes for the case of D R of 65% and CSR of 0.2.
The constitutive model was also calibrated based on the average axial strain accumulation rate (SAR). SAR is defined in units of percent of strain per cycle and indicates the half double-amplitude strain accumulation per cycle (average of triaxial extension and compression) after liquefaction triggering. Figure 27.3 presents a comparison of SAR observed during the triaxial tests and that from the simulations using the calibrated parameters for PM4SAND for the case of D R of 65%. For the PSC simulations, the strain accumulation rate (SAR) increases consistently with larger CSR. However, this trend is not consistent in the actual triaxial tests. Given this discrepancy, our goal was to approximate the average rate of SAR for CSR values between 0.14 and 0.2, which would be more relevant for the subsequent centrifuge tests. The observed versus simulated average rate of strain accumulation over these CSR values are 0.39% versus 0.42% strain per cycle, which are relatively in good agreement.

UBCSAND Constitutive Model
The UBCSAND constitutive model, which has been modified from the commonly available 904ar version by Fugro and Peter Burn, was also used in this study. This model has been modified, subsequently calibrated and validated by Fugro (Giannakou et al. 2011), and used in major projects such as the BART Offshore Transbay Tube Seismic Retrofit (Travasarou et al. 2011). The primary modification consists of introducing one additional model parameter (hfac4), which controls the rate of shear strain accumulation after liquefaction triggering. Detailed description of the model can be found in Beaty and Byrne (1998) and Byrne et al. (2004).
The model is fully defined by means of nine elastic and plastic parameters that can generally capture the observed response from monotonic and cyclic tests. Apart from N 1,60 which is a physical property of sand, four model parameters (hfac1 to hfac4; plastic hardeners) which are not related to physical properties are available to allow model calibration to observed behavior of different sands. Among these four parameters, the first three control the triggering of liquefaction (primarily the first two) while the fourth controls the post-trigger response. Table 27.2 summarizes the UBCSAND calibrated model parameters.
For calibration of the UBCSAND constitutive model, undrained cyclic direct simple shear (DSS) single element tests were simulated and compared against the DSS simulations using the PM4SAND calibrated properties. The UBCSAND model was not directly calibrated against the triaxial tests due to time limitations.
Liquefaction strength curves in terms of number of cycles to 98% pore pressure ratio, 1% shear strain, and 3% shear strain were compared against the cyclic simple shear simulations with the calibrated PM4SAND parameters. Figure 27.4 presents liquefaction strength curves for sand with D R of 65%. Plots of CSR versus normalized vertical effective stress, CSR versus shear strain, normalized vertical effective stress versus shear strain, and pore pressure ratio versus number of cycles are also provided in Fig. 27.5. For comparison purposes, the results from the cyclic DSS simulations using PM4SAND with the calibrated parameters are superimposed with the results from the DSS simulations using UBCSAND with the calibrated parameters.
Overall, the liquefaction strength curves implied by UBCSAND are generally steeper than those implied by PM4SAND and the laboratory tests. Hence, the UBCSAND model may need to be calibrated with a tighter target CSR range, corresponding to that induced by the input motions.

Introduction
During this phase of the project, the model parameters obtained as part of the calibration process were used to simulate the centrifuge experiments. The numerical simulations were conducted without any knowledge of the actual results. The simulation teams were only given information about the main characteristics of the centrifuge experiments, including the geometry of the centrifuge model, the Fig. 27.5 CSR versus normalized vertical effective stress, CSR versus shear strain, normalized vertical effective stress versus shear strain and excess pore pressure ratio versus number of cycles for D R~6 5% and CSR ¼ 0.14, DSS simulations using PM4SAND and UBCSAND centrifugal acceleration, the viscosity of the pore fluid (μ), the base excitation (horizontal and vertical components) and the achieved soil density.
Nine centrifuge experiments were selected by the organizing committee for LEAP-2017 Type-B simulations. The main characteristics of these experiments such as the centrifugal acceleration in g, the ratio of the pore fluid viscosity over the gravity (μ/g), and the achieved dry soil density (ρ d ) are summarized in Table 27.3 This table also contains the relative density (D R ) which is calculated for each test according to Eq. 27.1: where ρ d,min and ρ d,min are considered equal to 1485 and 1773 kg/m 3 , respectively based on the average values of LEAP-2017 laboratory data on Ottawa F65 sand (El Ghoraiby et al. 2017). Six centrifuge experiments were selected by Fugro to be simulated at this phase due to time limitations: Tests Ehime-2, NCU-3, UCD-1, and UCD-3 with a relative density (D R ) of approximately 65%, test KAIST-1 with a D R of approximately 80%, and test ZJU-2 with a D R of approximately 50%.

Analysis Platform
The analysis platform used for the numerical simulation of the centrifuge experiments is FLAC2D, Version 7.0 for PM4Sand and Version 8.0 for UBCSAND, respectively (Itasca 2011(Itasca , 2017. FLAC is a two-dimensional explicit finite difference program for engineering mechanics computation, following a "mixed discretization" scheme (Marti and Cundall 1982).

Numerical Modeling
Numerical modeling in FLAC was performed at the prototype scale. The model is 20 m long and 4 m high at the centerline. The soil surface forms a 5-degree slope. A horizontal ground water table with coordinates (À10, 0.875) and (10, 0.875) was assumed so that the slope is submerged. A baseline schematic of the centrifuge in prototype scale, illustrating the locations of pore pressure transducers, accelerometers and prescribed locations of displacement measurements, is presented in Fig. 27.6a. The grid and its dimensions of the numerical model are shown in Fig. 27.6b. The grid was built so that the grid points correspond to the locations of the accelerometers and the pore water pressure sensors illustrated in Fig. 27.6a.
A triangular pressure was applied at the soil surface simulating the load due to the water mass. Initially, during the gravity loading, Mohr-Coulomb properties were assigned to soil. A gravitational field with a vertical acceleration equal to 9.81 m/s 2 was applied and once equilibrium was reached, PM4SAND or UBCSAND constitutive models were assigned to soil elements and equilibrium was reached again. In the next step, the groundwater flow mode was set on and equilibrium was reached once more. During this process of gravity loading, the error tolerance was on the order of 10 À4 . The pore pressure was fixed at the soil surface. The horizontal displacement was also fixed at the vertical boundaries of the model. Last, the base of the model was constrained to move neither horizontally nor vertically. After the end of the gravity loading phase, horizontal and vertical mean input motions, obtained from recorders AH11-12 and AV1-2, respectively, were applied at the base and the vertical boundaries of the model. Rayleigh damping (stiffnessand mass-proportional) with a minimum value of 0.5% at 3.3 Hz was assigned to the soil elements. After the end of shaking, the analyses continued until the excess pore water pressures dissipated fully. This procedure was followed for each of the simulated centrifuge tests. For each simulation, the corresponding input motions were baseline corrected and then applied to the model.

Material Properties and Constitutive Model Parameters
The material properties assumed in the numerical simulations, such as permeability, the water bulk modulus, the soil relative density and the dry soil density, are given in Table 27.4. These model properties were adopted following sensitivity analyses evaluating the effect of changes in water bulk modulus, permeability anisotropy, and Raleigh damping.
After shaking, the model was allowed to reconsolidate until excess pore pressures dissipated and initial effective stresses were re-established. After end of shaking, permeability was scaled by ten to reduce calculation time, and the time was also scaled in the post-processing phase to ensure that the results are not affected compared to the solutions that would have been obtained by not scaling the permeability.
As aforementioned, the gravity loading was applied in two stages. In the first stage, Mohr-Coulomb properties were assigned to the soil elements. Basically, the gravity loading was initially applied elastically since no failure occurred.
In the second stage of gravity loading, as well as the dynamic part of the numerical analyses, advanced constitutive models were assigned to the soil elements, such as PM4SAND and UBCSAND models. Constitutive models PM4SAND and UBCSAND were calibrated for soil relative densities of 65%, 80%, and 90%. For the Type-B simulations of centrifuge tests with D R of about 65% (Ehime-2, NCU-3, UCD-1, and UCD-3) and 80% (KAIST-1), the initial calibrated parameters for a D R of 65% and 80% were used, respectively. There

Simulation Results
Figure 27.7 illustrates a comparison of the horizontal displacement and excess pore pressure ratio time histories at the middle surface of the model for Type-B simulation results using constitutive models UBCSAND, and PM4SAND versus Ehime 2 centrifuge experiment results. Both constitutive models and the centrifuge test results developed excess pore pressure ratio in excess of 80% while the rate of excess pore pressure dissipation recorded in the centrifuge is bounded by the two constitutive model predictions. Figure 27.8 presents a comparison of permanent horizontal displacement at the middle surface of the numerical models for Type-B simulation results using PM4SAND and UBCSAND versus the centrifuge experiment results for all cases simulated. Simulations with PM4SAND result in displacements closer to those measured and within a factor of 2 from the actual measurements with the exception of NCU3. Simulations with UBCSAND resulted in generally larger displacements than the actual measurements by a factor of about 2. This is likely attributed to the model calibration which was initially targeted at larger CSR ranges compared to those induced on average by the input ground motions. Overall, the blind predictions captured the centrifuge tests' behavior well and bounded the centrifuge response in most cases, indicating the benefits of using more than one constitutive model when performing numerical modeling.

Introduction and Constitutive Model Parameters
At this stage, the centrifuge experiment results were released to the teams who were provided with the opportunity to refine the simulations, if necessary. In the Type-C simulations, the UBCSAND constitutive model parameters hfac1 and hfac2 were increased by 40% (hfac1 ¼ 0.45 Â 1.4 ¼ 0.63 & hfac2 ¼ 0.5 Â 1.4 ¼ 0.7) and the rest of the parameters are the same with the ones presented in Table 27.2. The PM4SAND constitutive model parameters were not modified and are the same with the ones presented in Table 27.1. As discussed before, the liquefaction strength curves implied by the UBCSAND model are steeper than the ones implied by PM4SAND and the laboratory tests. Hence, the UBCSAND model had to be recalibrated for a specific CSR range. Figure 27.9 illustrates the liquefaction strength curves from DSS simulations at a vertical effective stress of 0.2 atm for constitutive models PM4SAND (red) and UBCSAND using Type-B (blue) and Type-C (cyan) calibration parameters.

Simulation Results
Recalibrating UBCSAND constitutive model parameters resulted in permanent displacements closer to the actual measurements. Figure 27.10 illustrates a

Introduction
A sensitivity analysis study was performed to illustrate the sensitivity of the simulations to key input parameters such as relative density, ground motion intensity and the influence of superimposed high frequencies in the ground motions. The properties of the seven sensitivity analyses (NS-1 to NS-7) are tabulated in Table 27.5. All sensitivity analyses were performed using the constitutive model parameters that were used in Type-C simulations discussed in the previous section. Simulations NS-1, NS-4, NS-5, NS-6, and NS-7 with soil relative density of 65% were performed using the Type-C calibrated parameters for the case of D R of 65%. Simulation NS-2 was performed using the Type-C calibrated parameters for the case of D R of 65% but using as input a D R of 50%. Simulation NS-3 was performed using the Type-C calibrated parameters for the case of D R of 80% but using as input a D R of 75%.

Sensitivity Analyses Results
Sensitivity analyses NS-1 to NS-3 were used to explore the sensitivity of simulations to relative density. Sensitivity analyses NS-1, 4, and 5 were used to explore the sensitivity of simulations to motion intensity. Sensitivity analyses NS-6 and NS-7 were used to assess the influence of superimposed ground motion high frequencies on simulation results. Table 27.5 presents the resulting horizontal displacement at the middle of the model surface and the duration of liquefaction at point P4 (P4 location: middle at 1-m depth) after the end of shaking. Figure 27.11 presents a comparison of sensitivity analyses using PM4SAND in terms of excess pore pressure ratio, horizontal displacement time history and amplification ratio at the middle surface of the model for NS1, NS6, and NS7. Input ground motions for analyses NS6 and NS7 have superimposed high frequencies, and input motion NS7 is equal to input motion NS6 multiplied by 1.43. This figure indicates that superimposed high frequencies in the input ground motion as well as motion intensity can have a significant effect on the model response.
The predicted displacements were evaluated against various intensity parameters, and meaningful trends were discovered between permanent displacements and input ground motion Arias intensity. Figure 27.12 presents the predicted horizontal displacements at the middle surface of the models versus Arias intensity for sensitivity analyses with D R of 65%. As expected, for higher Arias intensity motions, the predicted permanent displacements are higher. It appears that motions with high frequency content (NS-6 and NS-7) can result in greater displacements than motions with lower high frequency content and similar Arias intensity.
The sensitivity analyses also show that soil relative density has a significant impact on the model behavior.

Conclusions
We used two different constitutive models (PM4SAND and UBCSAND), implemented in the software FLAC, to simulate liquefaction-induced slope instability from six centrifuge tests. The methodology adopted for the constitutive model calibration process considered both liquefaction triggering and post-liquefaction accumulation of shear strains. Firstly, constitutive models, PM4SAND and UBCSAND, were calibrated against laboratory tests to comply with liquefaction triggering criteria such as the number of cycles required to reach specific strain and pore pressure ratio thresholds. Secondly, the strain accumulation rate (SAR) following liquefaction triggering was used as a target for model parameter calibration. Overall, the blind predictions captured the centrifuge tests behavior well and bounded the centrifuge response in most cases, indicating the benefits of using more than one constitutive model when performing numerical modeling. The predicted permanent displacements based on both constitutive models were generally within a factor of 2 from the actual measurements. The liquefaction strength curves implied by the UBCSAND model are steeper than the ones implied by PM4SAND and the laboratory test data. Hence, during the Type-C simulations, the UBCSAND model calibration was adjusted considering a CSR range corresponding to the CSR induced by the input motions. This should be considered when model calibrations are performed in practice, especially when multiple hazard levels are involved.
Differences between the predicted and measured responses are partially associated with limitations in the laboratory test data available, which affect the results of the numerical model calibration. For example, direct simple shear tests rather than triaxial tests, tests with static bias in addition to the no-bias conditions as well as tests at different CSR levels, representative of the CSR induced by the input motions, would be more relevant and are recommended. Differences between observed and computed responses are also attributed to simplifications in the prototype-scale numerical model (e.g., simulating possibly heterogeneous actual conditions with an idealized uniform relative density), accuracy of measurements and conversions from model to prototype scale, and differences between the input motions in the simulations ("average" of recorded at the centrifuge base) and the actual recorded input motions in the centrifuge which may vary across its base.
The sensitivity analyses suggest that horizontal displacements increase when high frequencies are introduced in the ground motion and for smaller soil relative density and higher ground motion intensity. insights from the study of an immersed railway tunnel. In 5th International Conference on Earthquake Geotechnical Engineering, Paper No. LITUR, January.
Open Access This chapter is licensed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license and indicate if changes were made. The images or other third party material in this chapter are included in the chapter's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the chapter's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.