FLAC3D New 3D Liquefaction Constitutive Model


Fully nonlinear methods can be formulated in terms of effective stresses and the fluid-mechanical interaction behaviors can be simulated before, during, and following dynamic shaking. Multi-dimensional nonlinear constitutive models are being used increasingly to simulate the phenomenon of liquefaction for sand-like geotechnical materials in earthquake engineering practice, e.g., from the relatively simple uncoupled models like the cyclic-counting Dames & Moore model (Dawson et al. 2001) to the relatively complex effective stress-based models like the PDMY02 model (Elgamal et al. 2003), DM04 (Sanisand) model (Dafalias & Manzari 2004) model, UBCSand model (Beaty & Byrne 2011), and PM4Sand model (Boulanger & Ziotopoulou 2015), to name some. While all of the aforementioned models have merit they also have some notable limitations and many of them are not formulated for 3D modeling.

P2P-Sand: An Introduction

Itasca is pleased to announce P2P-Sand (Practical Two-Surface Plasticity model for Sands), our new practice-friendly 3D liquefaction constitutive model for FLAC3D based on the DM04 (Sanisand) model (Dafalias & Manzari, 2004). By revising some formula without destroying the conciseness of the original DM04 (Sanisand) model, the modifications improve the comparative results between model simulations and various laboratory and field observations with various initial and loading conditions. In addition, Itasca's new liquefaction model embraces an easier and more practical calibration procedure in terms of in-situ data, instead of experimental data, greatly reducing the laboriousness of calibration.

Itasca's enhanced model:

  • retains the general 3D formulation, considering Lode angle effect, and thus can be used for general 3D boundary value problems
  • switches all void-ratio related formula, including the elastic modulus, plastic hardening modulus, and critical state, to be relative density related
  • requires only the initial soil relative density (Dr) and initial stresses to match the semi-empirical relation of CRR versus some in-situ parameter such as (N1)60 or qc1N when using the default parameters*
  • preserves the feature that one set of model constants is used to simulate different responses with different initial relative densities and initial stress states

The following table compares a number of factors for selected liquefaction constitutive models.


Single-zone FLAC3D models were used to test and evaluate this new liquefaction constitutive model considering the following:

  • CRR versus (N1)60 and number of cycles
  • Modulus reduction and damping
  • Undrained Cyclic Loading
  • Empirical Kσ curves
Values of CRR predicted by the model and compare very well to semi-empirical curve (left); as do curves of CSR versus number of cycles in undrained DSS cyclic simulations (right).
Model results fall within a reasonable range of G/Gmax ratios from constant mean pressure drained triaxial simulations for sands with and Dr = 35%, 55%, and 75%.
Models are well behaved given damping ratios from constant mean pressure drained triaxial simulations for sands with and Dr = 35%, 55%, and 75% and Po = 50 kPa, 300 kPa, and 1500 kPa.
Undrained cyclic DSS simulations for initial relative densities of 25%, 45%, 65%, and 80%.
Model results compare very well with Empirical Kσ curves (Boulanger and Idriss, 2004) for a range of overburden vertical effective stresses.

Some Examples

The following engineering and laboratory scale simulations using the new liquefaction model are performing very well and behaving as expected.

FLAC3D dynamic simulation of a wharf composed of structural elements resting on a slope subjected to shaking. Itasca's new liquefaction constitutive model was used for Dr = 50% and 65%.

More Information

For more details about Itasca's new liquefaction model, P2P-Sand, please see the following paper:

Cheng, Zhao (2018). "A practical 3D bounding surface plastic sand model for geotechnical earthquake engineering application," Geotechnical Earthquake Engineering and Soil Dynamics (GEESD) Congress V, June 10-13, 2018, Austin, Texas, United States.


Latest News
  • Workshop: The “Big Five”: Numerical Modeling of Cave Mining This MassMin2020 workshop will review caveability, ground subsidence, infrastructure stability, fragmentation, and gravity flow associated...
    Read More
  • 2020 H. Bolton Seed Lecture: Open Issues about Soil Liquefaction Dr. Bruce L. Kutter delivered the 2020 H. Bolton Seed Lecture at Geo-Congress 2020 in...
    Read More
  • Webinar: Using Rigid Blocks in PFC This webinar will provide a brief overview of PFC modeling using convex rigid blocks, including...
    Read More

Upcoming Events
1 Dec
WEBINAR: What’s New in 3DEC Version 7
This webinar is for people who have used 3DEC before and are interested in the latest developments.... Read More
4 Dec
Analyses of Embankment Dams and Slopes Using FLAC
This hands-on, virtual training course is 16 hours total, spread over four days in a 1.5-week period, and covers the analysis of embank... Read More
8 Mar
FLAC3D Remote Introductory Training - Online!
Three days of general feature training addressing basic concepts and recommended procedures for geotechnical numerical analysis.... Read More