Numerical Analysis of Stress State of Cemented Rockfill in Underground Backfilled Stopes Using Coupled FD/DE Methods

Chong Wei

Doctoral Candidate | Mining Engineering

University of Alberta (Canada)

Advisor: Derek Apel

Itasca Mentor: Tatyana Katsaga

Introduction

In underground stopes, energy alternation and stress transmission happens along the interfaces of backfilled stopes and the surrounding rockmass due to ubiquitous mechanical interactions between them. In this study, the stress state variation law in cemented rockfill and underground backfilled stopes, the strength characteristics of rockfill material, and the stability of backfilled stopes will be simulated using coupled PFC3D/FLAC3D software to develop linked discrete/continuum models that incorporate in-situ and laboratory experiments.

Five aspects of research contents, mainly about the strength characteristics of cemented rockfill and mechanical interactions between cemented rockfill and surrounding rocks in underground backfilled stopes, will be involved in this study. The research materials and data will be collected from a mine, and the laboratory experiments will be carried out in Rock Mechanics Lab at the University of Alberta.

Numerical Modeling of Cemented Rockfill Material under Uniaxial and Triaxial Compression Based on Discrete Element Method

In this part, material-modeling support package imbedded in software PFC3D was adopted, and numerical uniaxial and triaxial compression tests were carried out on both bonded-particle models (BPMs) with the simplified particle size distribution and the mean weighted particle diameter of the cemented rockfill (CRF) material, of which the microstructural parameters were both well calibrated to ensure the macro-mechanical properties of these BPMs perfectly correspond to the laboratory

testing data of the real CRF core specimens. The strength characteristics and failure pattern of the CRF BPMs under numerical triaxial compression were calibrated and the modeling results highly agreed with the laboratory testing results according to Mohr-Coulomb strength criterion. The particle size distribution and the model-particle diameter ratio of these BPMs also significantly influence the macro-mechanical properties and fracture evolution of the BPMs.

Numerical Modeling of Direct Shear Test on Cemented Rockfill Material Based on the Coupled Finite Difference and Discrete Element Method

This part provides an effective determination process of the interface microstructural parameters in the numerical modelling of the mechanical interaction between the backfill column and the surrounding rocks by using the coupled finite-difference and discrete-element (FD-DE) method. Three-dimensional numerical direct shear tests were conducted on the coupled FD-DE model of the combined rock-cemented rockfill (R-CRF) specimens with the frictional-bonded interface. With the microstructural parameters determined, the macro-mechanical properties of the interface of the coupled model were calibrated to well fit with the theoretical values of the combined specimen under the direct shear condition assumed according to the Mohr-Coulomb strength criterion. Evident elastoplastic shear behaviour was observed on the interfaces, and excellent agreement was attained between the shapes of the theoretical and simulated shear stress-displacement curves of the coupled models at different servo-controlled normal stresses. The test results validated the feasibility and rationality of the coupled FD-DE method for modelling the combined specimens with frictional-bonded interfaces or joints.

Numerical Modeling of Shear Behavior at Interface between Backfill Column and Surrounding Rock Based on Coupled Finite Difference and Discrete Element Methods

In this part, a 1000-meter-deep underground stope was excavated in the numerical model created with FLAC3D by finite difference method (FDM) and then backfilled with cemented rockfill (CRF) which was created with PFC3D by discrete element method (DEM). The shear behavior at the interface of the CRF column and surrounding rock was simulated by wall-zone coupling method (WZCM) through Socket 0/1 connections. The modeling results showed that during the settlement and consolidation of the CRF particles, weak binding effect formed at the interface and arching phenomenon occurred inside the CRF column, which significantly influenced the stress and displacement distributions inside the columns, meanwhile, backfilling with CRF efficiently improved the stress variation at the pillar and reduced the rockburst risk at the stress concentration zones.

Numerical direct shear test on a coupled rock/cemented-rockfill sample: a) model, b) shear stress distribution at the interface (normal stress of 0.5 MPa), and c) x-displ. distribution on the central cross-section at different shear displacement (ds).
Coupled FD-DE model of a single stope backfilled with the cemented rockfill: a) the view along the plane-xz, b) the view along the plane-yz, and c) backfilling process with the cemented rockfill.
Displacement distribution at the backfilled stope caused by the settlement process of the cemented rockfill column: a) vertical displacement and b) horizontal displacement.
Stress distribution at the backfilled stope influenced by the frictional bonded interfaces between the cemented rockfill column and the surrounding walls: a) vertical stress and b) horizontal stress.

Latest News
  • Itasca has announced the release of FLAC2D v9 Itasca has announced the release of FLAC2D v9, revolutionizing the way we analyze and predict...
    Read More
  • 6th Itasca Symposium on Applied Numerical Modeling The next Itasca Symposium will take place June 3 - 6, 2024, in Toronto, Canada....
    Read More
  • Itasca International Inc. announces the Selection of its New CEO Itasca International Inc. announces the Selection of its New CEO ...
    Read More

Upcoming Events
5 Nov
Python in Itasca Software
This course provides an overview of the Python programming language in Itasca software.The course covers major applications of Python t... Read More
19 Nov
Getting Started with FLAC2D/FLAC3D
This training is an introduction to continuous modeling with FLAC2D and FLAC3D. At the end of the course, participants will master the ... Read More