Options for 3DEC

3DEC simulates the nonlinear response of a system (soil, rock, and structures) to excitation from an external source or internal 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.

Dynamic analysis is required for the following types of geomechanical problems:

  • seismic (i.e., earthquake) loading;
  • explosive or impulsive loading;
  • mining problems involving seismic release of energy (i.e., rockbursts); and
  • flow of particles (angular or rounded) under gravity.

The fundamental assumption in all cases is that time is relevant. In almost all cases, the time of interest is less than one minute; in the case of explosive or impulsive loading, it is less than one second. 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).
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).

Natural dynamic systems contain some degree of damping of the vibrational energy within the system (due to friction, air or fluid resistance); otherwise, the system would oscillate indefinitely when subjected to driving forces. Two types of simulated damping (mass-proportional and stiffness-proportional) are available in 3DEC. 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.

In 3DEC, either form of damping may be used separately or in combination (i.e., Rayleigh damping). Maxwell damping provides an alternative to Rayleigh damping for time-domain seismic deformation analyses. The numerical calculation time using Maxwell damping can be a small fraction of that required for Rayleigh damping. Maxwell damping makes large 3D site response and soil-structure interaction analyses practical.

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. A dynamic input wizard is available for filtering and baseline correction of ground motions and to provide data necessary for dynamic analysis in 3DEC.

LEARN MORE

References

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.

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.

LEARN MORE

Comparison of features available with numerical vs analytical formulations.

Feature
Numerical
Analytical
Thermal mechanical coupling
Yes
Yes
Thermal fluid coupling
Yes
No
Thermal convection in fluid
Yes
No
Temperature boundary
Yes
No
Thermal flux boundary
Yes
No
Thermal convection boundary
Yes
No
Adiabatic boundary
Yes
Yes
Homogeneous thermal block properties
Yes
Yes
Anisotropic thermal block properties
Yes
No
Heat Sources* (e..g, nuclear waste)
Yes
Yes
3DEC model showing temperature contours and exaggerated mechanical displacement, due to thermal effects, along an underground drift above a nuclear waste repository.

*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.

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.

LEARN MORE

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

  • Itasca has announced the release of FLAC2D v9 Itasca has announced the release of FLAC2D v9, revolutionizing the way we analyze and predict...
  • 6th Itasca Symposium on Applied Numerical Modeling The next Itasca Symposium will take place June 3 - 6, 2024, in Toronto, Canada....