FLAC3D Theory and Background • Constitutive Models

Constitutive Models

The basic mechanical constitutive models provided in FLAC2D, FLAC3D, and 3DEC are arranged into null, elastic, and plastic model groups:

Null model group

  1. Null model

A null material model is used to represent material that is removed or excavated.

Elastic model group

  1. Elastic model

The elastic, isotropic model provides the simplest representation of material behavior. This model is valid for homogeneous, isotropic, continuous materials that exhibit linear stress-strain behavior with no hysteresis on unloading.

  1. Orthotropic model (**)

The elastic, orthotropic model represents material with three mutually perpendicular planes of elastic symmetry. For example, this model may simulate columnar basalt loaded below its strength limit.

  1. Anisotropic model

The elastic, transversely isotropic model gives the ability to simulate layered elastic media in which there are distinctly different elastic moduli in directions normal and parallel to the layers.

Plastic model group

  1. Von-Mises model

A basic elastoplastic model accounting for plastic deviatoric strain but neglecting plastic volumetric strain, typical elastoplastic model for metal-like materials.

  1. Drucker-Prager model

The Drucker-Prager plasticity model may be useful to model soft clays with low friction angles. However, this model is not generally recommended for application to geologic materials. It is included here mainly to permit comparison with other numerical program results.

  1. Mohr-Coulomb model

The Mohr-Coulomb model is the conventional model used to represent shear failure in soils and rocks. Vermeer and de Borst (1984), for example, report laboratory test results for sand and concrete that match well with the Mohr-Coulomb criterion.

  1. Ubiquitous-Joint model

The ubiquitous-joint model is an anisotropic plasticity model that includes weak planes of specific orientation embedded in a Mohr-Coulomb solid.

  1. Ubiquitous-Anisotropic

The ubiquitous-joint model is a model with a plane of weakness and with anisotropic elasticity.

  1. Strain-Softening/Hardening model

The strain-softening/hardening model allows representation of nonlinear material softening and hardening behavior based on prescribed variations of the Mohr-Coulomb model properties (i.e., cohesion, friction, dilation and tensile strength) as functions of the deviatoric plastic strain.

  1. Bilinear-Softening/Hardening Ubiquitous-Joint model

The bilinear strain-softening/hardening ubiquitous-joint model allows representation of material softening and hardening behavior for the matrix and the weak plane based on prescribed variations of the ubiquitous-joint model properties (i.e., cohesion, friction, dilation, and tensile strength) as functions of deviatoric and tensile plastic strain. The variation of material strength properties with mean stress can also be taken into account by using the bilinear option.

  1. Double-Yield model

The double-yield model is intended to represent materials in which there may be significant irreversible compaction in addition to shear yielding, such as hydraulically placed backfill or lightly cemented granular material.

  1. Modified Cam-Clay model

The modified Cam-Clay model may be used to represent materials when the influence of volume change on bulk property and resistance to shear need to be taken into consideration, as in the case of soft clay.

  1. Hoek-Brown model

A Hoek-Brown model provides an alternative to the Hoek-Brown-PAC model with a stress-dependent plastic flow rule, described above. The model characterizes post-failure plastic flow by simple flow rule choices given in terms of a user-specified dilation angle. This model also contains a tensile strength limit similar to that used by the Mohr-Coulomb model. In addition, a factor-of-safety calculation based on the strength reduction method can be run with the Hoek-Brown model.

  1. Hoek-Brown-PAC model

The Hoek-Brown failure criterion characterizes the stress conditions that lead to failure in intact rock and rock masses. The failure surface is nonlinear and is based on the relation between the major and minor principal stresses. The model incorporates a plasticity flow rule that varies as a function of the confining stress level.

  1. Cap-Yield (CYSoil) model (*)

The cap-yield (CY) soil model provides a comprehensive representation of the nonlinear behavior of soils. The model includes frictional strain-hardening and softening shear behavior, an elliptic volumetric cap with strain-hardening behavior, and an elastic modulus function of plastic volumetric strain. The model allows a more realistic representation of the loading/unloading response of soils.

  1. Simplified CYSoil (CHSoil) model (*)

A simplified version of the CYSoil model, called the CHSoil model, offers built-in features including a friction-hardening law that uses hyperbolic model parameters as direct input, and a Mohr-Coulomb failure envelope with two built-in dilation laws.

  1. Plastic-Hardening (PH) model (*)

An elasto-plasticity model with shear and volumetric hardening.

  1. Swell model (*)

A Mohr-Coulomb elasto-plasticity model considering wetting-induced deformation.

  1. Mohr-Coulomb-Tension (MohrT) model (**)

A Mohr-Coulomb elasto-plasticity model with tensile strain crack tracking.

  1. Soft-Soil model (*)

A elasto-plastic model for soft soils with high compressibility.

  1. NorSand model (*)

A critical state model applicable to soils in which particle to particle interactions are controlled by contact forces and slips rather than bonds.

  1. Finn model (*)

A simple Mohr-Coulomb type model accounting for the build-up of excess pore pressure during dynamic analyses.

  1. P2PSand model (*)

A Practical TWO-surface Plastic SAND constitutive model for general 3D geotechnical earthquake engineering application aimed at capturing essential soil dynamic characteristics.

  1. Concrete model

A plastic-damage concrete model.

  1. Columnar-Basalt (COMBA) Model (**)

A model accounting for the presence of up to four arbitrary orientations of weakness (ubiquitous joint) in a non-isotropic matrix.

  1. IMASS Model (**)

An Itasca Model for advanced strain softening

Other mechanical material models are provided with the optional features in FLAC2D and FLAC3D: There are eleven time-dependent creep constitutive models available in the creep model option for FLAC3D and ten for FLAC2D. Nine of the 3D creep constitutive models are included with 3DEC.

In addition, there are several pore-pressure generation material models available in the dynamic analysis option in FLAC2D and FLAC3D. The models are described in Dynamic Analysis.

There is also a modified version of the Drucker-Prager model that is provided to simulate the mechanical behavior associated with thermal hydration, called Hydration-Ducker-Prager model (FLAC2D, FLAC3D only).

The source codes of most of the models are included in the “\pluginfiles\cmodels” sub-directory. Users can modify these models or create their own constitutive models as DLLs by following the procedures given in Writing New Constitutive Models. The list of user-defined models (UDMs) is at Itasca UDM Website.

Fluid and thermal constitutive models are also included in FLAC2D and FLAC3D. Three fluid constitutive models are available. These models are described in Fluid-Mechanical Interaction. Four basic thermal constitutive models are provided. These models are described in Thermal Analysis. A thermal hydration model is also included.

The formulation of the basic mechanical models is addressed in general terms first in the section of Incremental Formulation. Separate discussions of the theoretical background and specific implementation for each model are then presented.

Note

*Not available in 3DEC.

**Not available in FLAC2D.