Optional Features

Optional features are specialized modules that can be added to UDEC at an additional cost. For more information or to purchase, please click here.

Barton-Bandis Joint Model

The Barton-Bandis joint model utilizes a series of empirical relations for joint normal behavior and joint shear behavior based on the effects of surface roughness on discontinuity deformation and strength as described by Barton (1982) and Bandis et al. (1985). The Barton-Bandis joint model encompasses the following features.

Joint Shear Behavior

  • Dilation as function of normal stress and shear displacement.
  • Joint damage due to post-peak shear.
  • Reduced secondary peak shear upon post-peak shear reversal.

Joint Normal Behavior

  • Hyperbolic stress-displacement path.
  • Hysteresis due to successive load/unload cycles.
  • Normal stiffness increase due to successive load/unload cycles.
  • Normal stiffness change due to surface mismatch caused by shear displacement.
  • Hydraulic aperture calculation based on joint closure and joint roughness.
Average dilation (mm) versus shear displacement (mm) for a direct shear test UDEC model.


Average shear stress (MPa) versus shear displacement (mm) for a direct shear test UDEC model.


Creep Material Models

The creep option can be used to simulate the behavior of materials that exhibit creep (i.e., time-dependent material behavior). The major difference between creep and other constitutive models is the concept of problem time in the simulation. For creep runs, the problem time and timestep represent real time, while for static analysis (in the other constitutive models), the timestep is an artificial quantity used only as a means of stepping to a steady-state condition. The timestep may be set by the user to a constant value, or controlled by UDEC to change automatically. If the timestep is changed automatically, it can be decreased whenever the maximum unbalanced force exceeds some threshold, and increased whenever it goes below some other level. For some of the creep models available, the creep rate is temperature dependent. Temperatures may be either specified as a model property or calculated during cycling using the thermal mode of the code. For either case, a temperature gradient also may be specified.

Unconfined compression tests are performed with a creep model as a demonstration simulating localization given an appropriate loading rate. Under a low strain rate (top) the response is monotonic with the sample deforming uniformly over 139 hours. However, given a 10x faster strain rate (bottom), shear bands are predicted to form. Although the maximum load is nearly the same for the two tests, the latter one exhibits global softening behavior.


Eight creep models have been implemented in UDEC as described below in order of increasing complexity.

Classical Viscoelastic (Maxwell) Model

Viscoelastic materials exhibit both viscous and elastic behaviors. The classical notion of Newtonian viscosity is that the rate of strain is proportional to stress. Stress-strain relations can be developed for viscous flow in a way similar to the way relationships are developed for elastic deformation.

Burgers Viscoelastic Model

The Burgers model is composed of a Kelvin model and a Maxwell model connected in series.

Two-Component Power Law

The Norton power law (Norton, 1929) is commonly used to model the creep behavior of salt or potash. The standard form of this law calculates the creep rate as a function of two material constants (A and n) and deviatoric stress:

WIPP Model

The WIPP (Waste Isolation Pilot Plant) model is a reference (empirical) viscoelastic creep law commonly used in thermomechanical analyzes associated with studies of underground nuclear-waste isolation in salt. This model describes the time- and temperature-dependent creep of natural rock salt.

Burgers Viscoplastic Model

The Burgers-creep viscoplastic model in UDEC is characterized by a visco-elasto-plastic deviatoric behavior and an elasto-plastic volumetric behavior. The viscoelastic and viscoplastic strain-rate components are assumed to act in series. The viscoelastic constitutive law corresponds to a Burgers model (Kelvin cell in series with a Maxwell component), and the plastic constitutive law corresponds to a Mohr-Coulomb model.

Power Law Viscoplastic Model

The viscoplastic model combines the behavior of the viscoelastic two-component Norton power law and the Mohr-Coulomb elastoplastic models.

WIPP Viscoplastic Model

Viscoplasticity is also modeled by combining the viscoelastic WIPP model with the Drucker-Prager plasticity model. Of the plasticity models currently available in UDEC, the Drucker-Prager model is the most compatible with the WIPP-reference creep law, because both models are formulated in terms of the second invariant of the deviatoric stress tensor.

Crushed Salt Model

The crushed-salt constitutive model is implemented in UDEC to simulate volumetric and deviatoric creep compaction behaviors. The model is a variation of the WIPP-reference creep law and is based on the model described by Sjaardema and Krieg (1987) with an added deviatoric component as proposed by Callahan and DeVries (1991).

User-Defined Models (UDM)

With this option, users may create their own contact or zone constitutive model for use in UDEC. The model must be written in C++ (CPP) and compiled as a DLL (dynamic link library) file using Microsoft Visual Studio 2010. The DLL then can be loaded whenever it is needed. The main function of the zone model is to return new stresses, given strain increments. The main function of the contact constitutive model is to return forces given displacements. However, the model must also provide other information (such as name of the model and material property names) and describe certain additional details about how the model interacts with the code. It is assumed that the user has a working knowledge of the C++ programming language.

UDMs may be used in other Itasca software provided this option is also available for that software. New DLL models can be obtained from an Itasca website devoted specifically to model development and exchange: www.itasca-udm.com.

Instructions for writing a constitutive model in C++ for operation in UDEC and the implementation of a DLL model are provided, including:

  • Descriptions of the base class and member functions;
  • Registration of models;
  • Information passed between the model and UDEC;
  • The model state indicators;
  • Descriptions of the support functions used by the model;
  • The source code for an example model;
  • FISH support for user-written models;
  • The mechanism for creating and loading a DLL; and
  • Microsoft Visual Studio 2010 constitutive model templates are also provided to assist in creating new constitutive models.