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