🎓 Lesson 2
D2
Fundamentals of Dynamic Process Modeling
Dynamic process modeling is like building a digital copy of a real mining or blasting operation that changes and responds over time—just like the real thing.
🎯 Learning Objectives
- ✓ Calculate blast-induced vibration peak particle velocity (PPV) using scaled-distance laws
- ✓ Design a dynamic charge initiation sequence to control burden movement timing
- ✓ Analyze time-resolved fragment size distribution evolution using kinetic energy partitioning models
- ✓ Explain how damping coefficients and rock mass stiffness influence model fidelity in blast simulation
- ✓ Apply dynamic boundary conditions in ANSYS AUTODYN or LS-DYNA to replicate field-scale blast response
📖 Why This Matters
In modern mining, a single mispredicted blast can cost $500k+ in rework, equipment damage, or regulatory penalties—and delay production by days. Dynamic process modeling lets engineers simulate *how* energy moves through rock *over milliseconds*, not just predict final fragment size. It’s the foundation for digital twins that update in real time using IoT sensor data from blast holes and accelerometers—turning reactive operations into predictive, adaptive systems.
📘 Core Principles
Dynamic modeling begins with conservation laws: mass, momentum, and energy applied to continuum or discrete media. In blasting, this translates to coupled thermomechanical modeling—where detonation chemistry (e.g., JWL equation of state) drives shock wave propagation, which induces elastic-plastic deformation, microcrack coalescence, and eventual macro-fracture. Key theoretical layers include: (1) wave propagation theory (impedance matching, reflection/transmission at interfaces), (2) rate-dependent rock strength models (e.g., Johnson-Holmquist), (3) time-domain coupling between explosive gas expansion and rock motion, and (4) validation via high-speed imaging and seismic array measurements. Fidelity increases with spatial resolution and temporal sampling—but must balance computational cost against decision-critical accuracy.
📐 Scaled-Distance Law for Peak Particle Velocity (PPV)
The scaled-distance law predicts ground vibration intensity based on charge weight and distance—a cornerstone for dynamic compliance and near-field modeling. It’s empirically derived but grounded in spherical wave attenuation theory and widely accepted for regulatory assessment and pre-blast modeling.
Dowding’s Scaled-Distance Equation
PPV = k × (R / W^{0.5})^{-b}Empirical relationship predicting peak particle velocity (mm/s) at distance R (m) from a blast of total charge weight W (kg), calibrated by site-specific constants k and b.
Variables:
| Symbol | Name | Unit | Description |
|---|---|---|---|
| PPV | Peak Particle Velocity | mm/s | Maximum ground vibration velocity measured orthogonal to wave propagation; primary metric for structural risk assessment. |
| R | Distance from Blast Source | m | Shortest horizontal distance from monitoring point to nearest charged hole center. |
| W | Maximum Weight per Delay | kg | Largest charge detonated within a single electronic delay interval—not total blast weight. |
| k, b | Site-Specific Attenuation Coefficients | dimensionless | k reflects geotechnical stiffness and damping; b reflects wave scattering and absorption (typically 1.3–2.0). |
Typical Ranges:
Hard, competent rock (e.g., quartzite): k = 150–300, b = 1.8–2.0
Weathered granite or schist: k = 400–600, b = 1.4–1.7
Soft sedimentary rock (e.g., shale): k = 700–1100, b = 1.2–1.5
💡 Worked Example
Problem: A surface blast uses 85 kg of ANFO in a single delay. A monitoring station is located 120 m from the nearest charged hole. Calculate PPV using k = 500 and b = 1.6 (typical for weathered granite).
1.
Step 1: Compute scaled distance SD = R / W^0.5 = 120 / √85 ≈ 120 / 9.22 ≈ 13.01 m/kg⁰·⁵
2.
Step 2: Apply Dowding’s formula: PPV = k × SD^(-b) = 500 × (13.01)^(-1.6)
3.
Step 3: Calculate exponent: 13.01^1.6 ≈ 42.7 → PPV ≈ 500 / 42.7 ≈ 11.7 mm/s
Answer:
The predicted PPV is 11.7 mm/s, which falls within the safe range of <25 mm/s for residential structures per USBM RI 8507.
🏗️ Real-World Application
At Newmont’s Boddington Mine (Western Australia), dynamic process modeling integrated microseismic event timing, high-speed video of bench face throw, and downhole pressure transducer data to calibrate a 3D AUTODYN model of a 12-row cast blast. The model revealed a 42-ms delay misalignment between rows 5 and 6 caused premature confinement loss—leading to oversize and increased flyrock. Revised initiation timing (adjusted from 25 ms to 38 ms inter-row delay) improved fragmentation uniformity by 27% and reduced oversize (>76 cm) by 41%, validated by post-blast LiDAR and crusher throughput analysis.