Calibration of a discrete element model for intact rock up to its peak strength

When three dimensional, bonded discrete element models (DEMs) are deployed to model intact rock, a basic question is how to determine the micro parameters that control macro properties of the modeled rock. After briefly describing the authors' DEM code, this paper describes algorithms to calibrate the model's micro parameters against standard laboratory tests, such as uniaxial and triaxial tests. Sensitivity analysis is used to identify the deformability micro parameters by obtaining relationships between microscopic and macroscopic deformability properties. The strength model parameters are identified by a global optimization process aimed at minimizing the difference between computed and experimental failure envelopes. When applied to the experimental results of Lac du Bonnet granite, this calibration process produced a good agreement between simulated and experimental results for both deformability and strength properties. Copyright © 2009 John Wiley & Sons, Ltd.

[1]  T. Ng Input Parameters of Discrete Element Methods , 2006 .

[2]  Eugenio Oñate,et al.  Combination of discrete element and finite element methods for dynamic analysis of geomechanics problems , 2004 .

[3]  Martin Berggren,et al.  Hybrid differentiation strategies for simulation and analysis of applications in C++ , 2008, TOMS.

[4]  A. Fakhimi,et al.  Application of Dimensional Analysis in Calibration of a Discrete Element Model for Rock Deformation and Fracture , 2007 .

[5]  V. Hajiabdolmajid,et al.  Modelling brittle failure of rock , 2002 .

[6]  C. Martin,et al.  The strength of massive Lac du Bonnet granite around underground openings , 1993 .

[7]  F. Donze,et al.  Identification and Validation of a Discrete Element Model for Concrete , 2004 .

[8]  N Oreskes,et al.  Verification, Validation, and Confirmation of Numerical Models in the Earth Sciences , 1994, Science.

[9]  C. Martin,et al.  A clumped particle model for rock , 2007 .

[10]  P. Cundall A computer model for simulating progressive, large-scale movements in blocky rock systems , 1971 .

[11]  John R. Williams,et al.  Discrete Element Modeling Applied to Laboratory Simulation of Near-Wellbore Mechanics , 2004 .

[12]  Arnold Neumaier,et al.  SNOBFIT -- Stable Noisy Optimization by Branch and Fit , 2008, TOMS.

[13]  Jeoung Seok Yoon,et al.  Application of experimental design and optimization to PFC model calibration in uniaxial compression simulation , 2007 .

[14]  Ching S. Chang,et al.  Stress-strain relationship for granular materials based on the hypothesis of best fit , 1997 .

[15]  P. A. Cundall,et al.  NUMERICAL MODELLING OF DISCONTINUA , 1992 .

[16]  F. Donze,et al.  Modeling fractures in rock blasting , 1997 .

[17]  P. Cundall,et al.  A bonded-particle model for rock , 2004 .

[18]  Dimitri Debruyne,et al.  Elasto-plastic material parameter identification by inverse methods: Calculation of the sensitivity matrix , 2007 .

[19]  P. Cundall,et al.  A discrete numerical model for granular assemblies , 1979 .