3DEC Options

Optional Features


The Itasca Constitutive Model for Advanced Strain Softening (IMASS) has been developed to represent the rock mass response to excavation induced stress changes. IMASS represents the damage around an excavation, slope, or caving process by accounting for the progressive failure and disintegration of the rock mass from intact, jointed, and/or veined rock to a disaggregated, bulked material. IMASS is based on empirical relationships and uses strain and zone-size dependent properties that reflect the impacts of dilation and bulking as a rock mass undergoes plastic deformation.

IMASS uniquely contains two softening (or residual) yield envelopes to represent the two-stage softening behavior for a rock mass that distinguishes between damage (caused by fracturing and the associated loss of cohesion and tensile strength) and the subsequent disturbance (due to bulking) in rock mass behavior. This two-stage softening/weakening behavior in IMASS is critical to accurately represent the rock mass post-peak behavior for underground and surface mining applications.

IMASS is available as a built-in, optional constitutive model for 3DEC (version 7.0 or later) and is sold as a separate, monthly or annual lease, license.


Conceptual model of the rock mass response to caving captured by IMASS

Dynamic Analysis

3DEC simulates the nonlinear response of a system (soil, rock, and structures) to excitation from an external (e.g., seismic) source or internal (e.g. vibration or blasting) sources. It can reproduce the evolution of permanent movements due to yield. This option models the full dynamic response of a system in the time domain. Capabilities include specification of velocity or stress-wave input, quiet (i.e., viscous) boundaries, free-field conditions (ideal for earthquake simulation), and damping.

3DEC model of an arch dam for earthquake analysis. The free-field boundary conditions are displayed along the model boundaries (Lemos, 2012).

Problems such as seismic loading, explosive loading, seismic release of energy, and flow of particles may be modeled.

Comparison between tunnel pillars in jointed rock after an accidental explosion in one tunnel. Shown are deformations with rock reinforcement (top) and large blocks deformations without ground support (bottom). (Brandshaug and Rosengren, 2008).

Two types of damping are available in 3DEC: mass-proportional and stiffness-proportional. Mass-proportional damping applies a force which is proportional to absolute velocity and mass, but in the direction opposite to the velocity. Stiffness proportional damping applies a force, which is proportional to the incremental stiffness matrix multiplied by relative velocities or strain rates, to contacts or stresses in zones. Either form of damping may be used separately or in combination (i.e., Rayleigh damping).

In 3DEC, the dynamic input can be applied as either a prescribed velocity history or as a stress history. An acceleration history needs to be integrated numerically first to produce a velocity history for 3DEC.

Thermal Analysis

The thermal option in 3DEC allows the simulation of transient heat conduction. There are two separate formulations of the thermal logic. The first is a numerical formulation using the explicit or implicit finite difference method. This method is more accurate for short times, and includes thermal-mechanical fluid coupling. The second is an analytical formulation that uses superposition of point heat sources* in an infinite medium. This method is suitable for long thermal times, and is very fast.

Comparison of features available with numerical vs analytical formulations.

Thermal mechanical couplingYesYes
Thermal fluid couplingYesNo
Thermal convection in fluidYesNo
Temperature boundaryYesNo
Thermal flux boundaryYesNo
Thermal convection boundaryYesNo
Adiabatic boundaryYesYes
Homogeneous thermal block propertiesYesYes
Anisotropic thermal block propertiesYesNo
Heat Sources* (e..g, nuclear waste)YesYes

*Point heat sources may be placed individually, in lines, or in grids, to represent point, line, or plane sources of heating. This formulation yields rapid calculations, correct application of mechanical boundary conditions, incorporation of the infinite thermal boundary,and the ability to use inhomogeneous and anisotropic mechanical properties.

3DEC model showing temperature contours and exaggerated mechanical displacement, due to thermal effects, along an underground drift above a nuclear waste repository.

User-Defined Constitutive Models

User-defined constitutive models can be written in C++ for both zoned block materials and joint materials. These are compiled as DLL files that can be loaded whenever needed with this option. Microsoft Visual Studio 2017 or 2019 is used to compile the DLL files. The main function of the constitutive model is to return new stresses, given strain increments. However, the model must also provide other information (such as name of the model and material property names) and describe certain details about how the model interacts with the code. Itasca maintains an online library of UDM C++ models where users can submit and download novel and useful constitutive models.


Brandshaug, T. and L. Rosengren (2008). 3D Numerisk Analys av Explosionslaster I Bergtunnlar (3D numerical analyses of accidental explosions in rock tunnels). Swedish Rock Engineering Research, SveBeFo Rapport 89, Stockholm, ISSN 1104 – 1773, 225 pages (in Swedish).

Lemos, J. (2012). Modelling the failure modes of dams' foundations. In MIR 2012 - Nuovi metodi diindagine, monitoraggio e modellazione degli amassi rocciosi (Eds. G. Barla, M. Barla, A.M. Ferrero, T. Rotonda), Politecnico di Torino, Italy, 2012, pp. 259-272.

Latest News
  • Itasca Celebrates 40 Years Itasca is celebrating 40 years of solving geomechanical and hydrogeological challenges through engineering and computer...
    Read More
  • Stability and Stress-Deformation Analyses of Reinforced Slope Failure at Yeager Airport This paper describes the material properties along with the inverse limit-equilibrium and permanent deformation analyses...
    Read More
  • Computers and Geotechnics: Scott Sloan Best Paper Award for 2019 Itasca is pleased to congratulate Dr. Branko Damjanac and Dr. Peter A, Cundall for their...
    Read More