UDEC - Version 7.0 - Universal Distinct Element Code Software
The Universal Distinct Element Code (UDEC) is a two-dimensional numerical program that simulates the quasi-static or dynamic response to loading of media containing multiple intersecting joint structures. The discontinuous medium is represented as an assemblage of discrete blocks while the discontinuities are treated as boundary conditions between blocks. Large displacements along discontinuities and rotations of blocks can occur. UDEC utilizes an explicit solution scheme that can model complex, nonlinear behaviors.
Models may contain a mix of rigid or deformable blocks. Deformable blocks are defined by a continuum mesh of finite-difference zones, with each zone behaving according to a prescribed linear or nonlinear stress-strain law. The relative motion of the discontinuities is also governed by linear or nonlinear force-displacement relations for movement in both the normal and shear directions. Joint models and properties can be assigned separately to individual, or sets of, discontinuities.
Because UDEC is not limited to a particular type of problem or initial condition, it may be applied to a wide variety of physical behaviors or any case where an understanding of the two-dimensional response of such structures is needed. UDEC is capable of simulating a wide range of engineering and scientific analyzes including: stability analysis of jointed rock slopes or underground excavations; fluid or gas flow through jointed rock; stability of masonry structures, dams, and foundations; blasting, earthquakes, and microseismicity; among many more applications.
Optional features are specialized modules that can be added to UDEC at an additional cost for additional simulation tools.
- Barton-Bandis Joint Analysis: 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).
- Creep Material Analysis: The creep option can be used to simulate the behavior of materials that time-dependent material behavior.
- User-Defined Models (UDM): With this option, users may create their own contact or zone constitutive model for use in UDEC using C++ scripting.
- Thermal Analysis: The thermal model simulates the transient flux of heat in materials and the subsequent development of thermally induced stresses. The heat flux is modeled by either isotropic or anisotropic conduction. Heat sources can be added and can be made to decay exponentially with time.
- Designed to accommodate any kind of analyzes of jointed, blocky systems.
- Simulates and tracks large displacements, including slip and opening, along a discontinuous medium (e.g., jointed rock mass, damage models of intact materials).
- Models may contain a mixture of rigid or deformable blocks.
- Deformable blocks consist of a finite-difference mesh, permitting simulation of soils, intact rock, engineered materials and rock masses as continuums).
- Welded construction cracks can be specified to aid model building and permit the simulation of joints that terminate within blocks (e.g., discrete fracture networks).
- Excavation, backfill simulation and “infinite domain” problems can be all be modeled.
- Includes automatic joint set and Voronoi joint generators.
- Includes a built-in scripting language (FISH) and text-editor to customize or automate virtually all aspects of program operation, including user-defined properties and other model variables.
- Libraries of materials and constitutive models for deformable blocks and joints are included.
Special Features come Standard
The following features are included as part of the standard UDEC license.
Structural elements - Structural Elements provide the ability to model support in UDEC. Reinforcement is modeled using cables that may be part of a local or global reinforcement model. Surface support is provided through beam elements (which may represent concrete lining, shotcrete, steel sets, etc.) and one-dimensional support members (which may be hydraulic or wooden props or packs).
Fluid and Gas Flow - UDEC can be used to model fluid or gas flow through the fractures of a system of impermeable blocks. A fully coupled mechanical-hydraulic analysis is performed in which fracture conductivity is dependent on mechanical deformation, and in return joint water pressures affect the mechanical computations. Confined flow, transient flow, two-phase flow and flow with a free surface can be modeled using UDEC.
Thermal - UDEC's thermal capability simulates transient heat conduction in materials and the development of thermally induced displacements and stresses. Heat transfer is modeled as conduction. Several different thermal boundary conditions may be imposed. Any mechanical model may be used with the thermal model. Heat sources may be inserted into the material and may be made to decay with time. Implicit and explicit calculation schemes are available and may be interchanged. The thermal analysis may be coupled to the mechanical and the fluid calculations.
Dynamic - UDEC can be used to simulate dynamic models for seismic, blasting and earthquake analyzes. Absorbing (quiet) and free-field boundaries and seismic wave input are all available. The dynamic formulation can be coupled to structural element, fluid flow and thermal models.
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.
The thermal model simulates the transient flux of heat in materials and the subsequent development of thermally induced stresses. The heat flux is modeled by either isotropic or anisotropic conduction. Heat sources can be added and can be made to decay exponentially with time.
UDEC allows simulation of transient heat conduction in materials and the development of thermally induced displacements and stresses. This includes the following specific features.
- Heat transfer is modeled as conduction, either isotropic or anisotropic, depending on the user’s choice of material properties.
- Several different thermal boundary conditions may be imposed.
- Any of the mechanical block models may be used with the thermal model.
- Heat sources may be inserted into the material as volume sources. The sources may be made to decay exponentially with time.
- Both implicit and explicit calculations schemes are available, and the user can switch from one to the other at any time during a run.
- The thermal analysis provides one-way coupling to the mechanical stress calculation through the thermal expansion coefficient.
- The thermal analysis provides one-way coupling to the calculation for fluid flow in joints through the temperature dependency of fluid density and joint permeability.
The following example shows a transient thermal-mechanical simulation of the behavior of a nuclear-waste emplacement drift in which heat-producing waste is placed vertically beneath the panel floor (Christianson, 1989). The emplacement drift under study is in the center of an emplacement panel. Waste canisters are placed in the floor of the drift. Using model symmetry, only one-half of the disposal room-and-pillar needs to be included in the analysis. The thermal boundary conditions are considered to be adiabatic.
The tributary heating area for the emplacement panel is reported to be 8194.5 m2. The average thermal loading is considered to be 14.1 W/m2. For the panel geometry, this results in an initial heat-generating power per meter of room length of 713.5 W. The initial power of a waste container at the time of emplacement varies from 0.42 kW to 3.2 kW. The following UDEC model shows the rock mass around the panel with the temperate distribution and thermal-induced rock displacements after 50 years since waste emplacement.
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:
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: User-Defined Constitutive Models (UDM).
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.