Calculator
D5
Time-Dependent Deformation: Creep, Relaxation & Squeezing Ground
Rocks can slowly deform over time under constant load—like clay under pressure—causing tunnels or slopes to gradually close, crack, or shift.
Typical Scale
Convergence rates range from 0.1 mm/day (mild) to >5 mm/day (severe)
Key Standards
ISRM Suggested Methods, ASTM D7012 (creep testing), USBR EM-1110-2-2001
Monitoring Threshold
Support systems designed for ≤1 mm/day sustained convergence are considered serviceable
⚠️ Why It Matters
1
High in-situ stress in deep excavations
↓
2
Activation of time-dependent rheology in soft rocks or clay-rich fault zones
↓
3
Unanticipated convergence exceeding support capacity
↓
4
Premature yielding of primary supports (steel sets, shotcrete)
↓
5
Catastrophic serviceability loss or collapse
↓
6
Costly re-excavation, redesign, and extended project delays
📘 Definition
Time-dependent deformation in rock masses encompasses creep (progressive strain under constant stress), stress relaxation (gradual stress reduction under constant strain), and squeezing ground (convergent deformation in weak, ductile, or highly jointed rock subjected to high in-situ stress). These phenomena arise from viscoelastic, viscoplastic, and damage-accumulation mechanisms within the rock matrix and discontinuity networks, and are critically sensitive to stress state, temperature, pore fluid pressure, and time scale.
🎨 Concept Diagram
AI-generated illustration for visual understanding
💡 Engineering Insight
Squeezing is rarely a 'rock property' issue—it’s a *system response* to misaligned support stiffness. Yielding steel sets often outperform rigid concrete linings not because they're stronger, but because their controlled deformation absorbs energy without triggering brittle failure in the surrounding rock mass. Always design support to match the *rate* of convergence—not just its magnitude.
📖 Detailed Explanation
Time-dependent deformation begins with microscopic processes: dislocation glide in quartz/feldspar, interlayer slip in phyllosilicates (e.g., illite), and pressure solution along grain boundaries—all accelerated by pore fluid activity and elevated temperature. In engineering terms, this manifests as slow, non-recoverable strain even below peak strength.
At the rock mass scale, discontinuities dominate behavior: sheared faults with clay gouge exhibit power-law creep, while intact blocks between joints may behave elastically. The transition from stable to squeezing ground occurs when the excavation-induced stress redistribution exceeds the time-dependent yield envelope—often only apparent after days or weeks, not hours.
Advanced analysis requires coupling rheology with damage mechanics: microcrack coalescence under sustained load reduces effective stiffness, accelerating secondary creep and potentially triggering tertiary (accelerated) creep—the precursor to macroscopic instability. Modern practice integrates acoustic emission monitoring with inverse modeling to update Burgers parameters in real time during tunneling.
At the rock mass scale, discontinuities dominate behavior: sheared faults with clay gouge exhibit power-law creep, while intact blocks between joints may behave elastically. The transition from stable to squeezing ground occurs when the excavation-induced stress redistribution exceeds the time-dependent yield envelope—often only apparent after days or weeks, not hours.
Advanced analysis requires coupling rheology with damage mechanics: microcrack coalescence under sustained load reduces effective stiffness, accelerating secondary creep and potentially triggering tertiary (accelerated) creep—the precursor to macroscopic instability. Modern practice integrates acoustic emission monitoring with inverse modeling to update Burgers parameters in real time during tunneling.
🔄 Engineering Workflow
Step 1
Step 1: In-situ stress mapping (hydraulic fracturing + overcoring)
→
Step 2
Step 2: Core sampling from critical horizons (depth-specific, oriented, low-damage recovery)
→
Step 3
Step 3: Time-dependent triaxial creep testing (3–7 stress levels, up to 10⁴ h duration)
→
Step 4
Step 4: Squeezing index assessment & RMR/Q-system recalibration with time-effects weighting
→
Step 5
Step 5: Numerical modeling using Burgers or Nishihara viscoplastic constitutive laws
→
Step 6
Step 6: Support system design iteration (yielding capacity vs. closure rate match)
→
Step 7
Step 7: Field validation via convergence monitoring and adaptive re-design
📋 Decision Guide
| Rock/Field Condition | Recommended Design Action |
|---|---|
| SI > 0.6 + high horizontal stress ratio (K₀ > 2.5) in argillaceous rock | Adopt full-face SEM with 0.8–1.2 m advance, steel arches with yielding connections, and immediate shotcrete with fiber reinforcement. |
| Measured ε̇ₛ > 5×10⁻⁸ s⁻¹ at 70% σ_c in laboratory tests on core samples | Install real-time convergence monitoring (total stations + extensometers) with automated alarm at 0.5 mm/day threshold. |
| RMR < 30 + presence of active shear zones with clay gouge (>30% smectite) | Prevent squeezing via grouted forepoling + radial bolts before face advance; avoid smooth blasting to minimize disturbance. |
📊 Key Properties & Parameters
Creep Compliance (J(t))
10⁻⁴ to 10⁻² MPa⁻¹ (at t = 1–1000 h, depending on rock type)Measure of time-dependent strain per unit applied stress, typically expressed as J(t) = ε(t)/σ₀ for uniaxial loading.
⚡ Engineering Impact:
Directly governs long-term convergence rate predictions for tunnel linings and pillar stability.
Secondary Creep Rate (ε̇ₛ)
10⁻⁹ to 10⁻⁶ s⁻¹ (for shale, phyllite, and mylonite at 50–80% σ_c)Steady-state strain rate during the linear phase of creep, reflecting dominant viscoplastic flow.
⚡ Engineering Impact:
Used to calibrate numerical models for support design lifetime and allowable closure thresholds.
Squeezing Index (SI)
0.2–0.9 (SI > 0.4 indicates severe squeezing risk)Empirical parameter defined as SI = (σₕₘₐₓ / UCS) × (1 − RMR/100), quantifying propensity for squeezing behavior.
⚡ Engineering Impact:
Triggers mandatory use of sequential excavation methods (SEM), forepoling, and yieldable supports.
Burgers Model Parameters (η₁, η₂, E₁, E₂)
E₁: 2–15 GPa; E₂: 0.5–5 GPa; η₁: 10¹²–10¹⁴ Pa·s; η₂: 10¹⁰–10¹² Pa·sFour-parameter viscoelastic model combining instantaneous elasticity, delayed elasticity, and two viscous dashpots to simulate multi-stage creep.
⚡ Engineering Impact:
Required input for FLAC2D/FLAC2D and Phase2 time-dependent simulations of deep tunnel advance rates.
📐 Key Formulas
Squeezing Index (SI)
SI = (σₕₘₐₓ / UCS) × (1 − RMR/100)Empirical indicator of squeezing potential in deep tunnels.
Variables:
| Symbol | Name | Unit | Description |
|---|---|---|---|
| SI | Squeezing Index | Empirical indicator of squeezing potential in deep tunnels | |
| σₕₘₐₓ | Maximum Horizontal Stress | MPa | Maximum principal horizontal stress in the rock mass |
| UCS | Uniaxial Compressive Strength | MPa | Unconfined compressive strength of the rock material |
| RMR | Rock Mass Rating | Empirical rock mass classification index ranging from 0 to 100 |
Typical Ranges:
Moderate squeezing
0.4–0.6
Severe squeezing
0.6–0.9
⚠️ SI < 0.3 indicates negligible squeezing risk
Norton Power-Law Creep
ε̇ = A × σⁿ × exp(−Q/RT)Predicts steady-state creep rate as function of stress, temperature, and activation energy.
Variables:
| Symbol | Name | Unit | Description |
|---|---|---|---|
| ε̇ | creep strain rate | s⁻¹ | rate of plastic deformation under constant stress and temperature |
| A | material constant | s⁻¹·Pa⁻ⁿ | pre-exponential factor dependent on material microstructure |
| σ | applied stress | Pa | mechanical stress causing creep deformation |
| n | stress exponent | dimensionless exponent reflecting stress sensitivity of creep mechanism | |
| Q | activation energy for creep | J/mol | energy barrier for thermally activated creep process |
| R | universal gas constant | J/(mol·K) | fundamental physical constant |
| T | absolute temperature | K | thermodynamic temperature |
Typical Ranges:
Argillite at 25°C
A = 10⁻²⁵–10⁻²³ s⁻¹·MPa⁻ⁿ; n = 3–5
Welded tuff at 25°C
A = 10⁻²⁷–10⁻²⁵ s⁻¹·MPa⁻ⁿ; n = 4–6
⚠️ ε̇ₛ > 10⁻⁸ s⁻¹ warrants time-dependent support design
🏭 Engineering Example
Yucca Mountain Exploratory Studies Facility (ESF), Nevada, USA
Tuff (welded Topopah Spring Tuff)SI
0.68
RMR
42
UCS
75 MPa
η₂
3.8×10¹¹ Pa·s
ε̇ₛ
1.2×10⁻⁷ s⁻¹
Convergence Rate (measured)
1.8 mm/day at 30 m from face
🏗️ Applications
- Deep mining access tunnels (e.g., Kidd Creek, Canada)
- Geological repository drifts (e.g., Onkalo, Finland)
- Hydropower headrace tunnels in metamorphic belts (e.g., Chamera III, India)
🔧 Try It: Interactive Calculator
📋 Real Project Case
Deep-Level Gold Mine Rockburst Mitigation
Mponeng Mine, South Africa — 4.2 km depth expansion
Challenge: Frequent high-energy rockbursts causing fatalities and equipment damage
Read full case study →
❓ Frequently Asked Questions
What is the difference between creep, stress relaxation, and squeezing ground in rock masses?
Creep is time-dependent strain accumulation under constant stress (e.g., tunnel convergence under steady overburden). Stress relaxation is the gradual reduction in stress under constant strain (e.g., decreasing support load in a deforming liner held at fixed displacement). Squeezing ground is a field-scale manifestation—convergent, often rapid, plastic deformation of weak, ductile, or highly fractured rock around excavations under high in-situ stress—driven by combined creep, relaxation, and brittle damage processes.
Why do some rocks exhibit significant time-dependent deformation while others do not?
Time-dependent behavior depends on mineral composition (e.g., phyllosilicates like illite or chlorite enhance viscoplasticity), microstructure (e.g., grain size, clay content, pre-existing fractures), and environmental conditions: elevated temperature accelerates diffusion-controlled mechanisms (e.g., pressure solution); high pore fluid pressure reduces effective stress and lubricates slip surfaces; and high differential stress promotes dislocation glide and intergranular flow. Competent, dry, crystalline rocks (e.g., granite) typically show minimal creep at shallow depths and ambient temperatures.
How does pore fluid pressure influence creep and squeezing behavior?
Elevated pore fluid pressure reduces effective stress (σ' = σ − u), weakening the rock matrix and discontinuities. It enhances pressure solution, facilitates clay mineral hydration and interlayer slip, and promotes hydrolytic weakening of silicate bonds. In squeezing ground, high pore pressure—especially in low-permeability, overpressured formations—can trigger delayed instability, accelerate strain rates, and reduce the time to failure, making deformation more sensitive to excavation-induced stress redistribution.
Can time-dependent deformation be predicted or modeled for underground engineering projects?
Yes—using constitutive models such as Burger’s (viscoelastic), Nishihara’s (viscoelastoplastic), or generalized Maxwell-type models calibrated with laboratory creep tests (e.g., triaxial or uniaxial creep under controlled T, P, and u). Field monitoring (convergence, extensometer, and piezometer data) combined with back-analysis improves model reliability. However, prediction remains challenging due to scale effects, heterogeneity, and evolving discontinuity behavior—requiring conservative design allowances (e.g., staged support, time-dependent yield criteria) in high-risk squeezing zones.
What are practical mitigation strategies for squeezing ground during tunneling?
Key strategies include: (1) preemptive ground improvement (e.g., forepoling, grouting, or ground freezing) to increase stiffness and reduce permeability; (2) flexible, deformable support systems (e.g., yielding steel arches, shotcrete with fiber reinforcement) that accommodate controlled convergence; (3) sequential excavation with minimal disturbance (e.g., NATM principles); (4) real-time monitoring and adaptive design; and (5) drainage where feasible to lower pore pressure—though caution is needed to avoid destabilizing clay-rich or swelling strata.
🎨 Technical Diagrams
📚 References
[1]
ISRM Suggested Methods for the Quantification of Time-Dependent Deformation of Rock
— International Society for Rock Mechanics (ISRM)
[2]
Rock Slope Engineering (4th ed.)
— ERM Ltd. & CRC Press
[3]
Guidelines for the Design of Tunnel Support Systems in Squeezing Ground
— U.S. Bureau of Reclamation (USBR)