Calibrating the stress-time curve of a combined finite-discrete element method to a Split Hopkinson Pressure Bar experiment

Abstract We present a generic method for automatically calibrating a computer code to an experiment, with uncertainty, for a given “training” set of computer code runs. The calibration technique is general and probabilistic, meaning the calibration uncertainty is represented in the form of a probability distribution. We demonstrate the calibration method by calibrating a combined Finite-Discrete Element Method (FDEM) to a Split Hopkinson Pressure Bar (SHPB) experiment with a granite sample. The probabilistic calibration method combines runs of a FDEM computer simulation for a range of “training” settings and experimental uncertainty to develop a statistical emulator. The process allows for calibration of input parameters and produces output quantities with uncertainty estimates for settings where simulation results are desired. Input calibration and FDEM fitted results are presented. We find that the maximum shear strength σ t max and to a lesser extent maximum tensile strength σ n max govern the behavior of the stress-time curve before and around the peak, while the specific energy in Mode II (shear) E t largely governs the post-peak behavior of the stress-time curve. Good agreement is found between the calibrated FDEM and the SHPB experiment. Interestingly, we find the SHPB experiment to be rather uninformative for calibrating the softening-curve shape parameters (a, b, and c). This work stands as a successful demonstration of how a general probabilistic calibration framework can automatically calibrate FDEM parameters to an experiment.

[1]  M. Paterson Experimental Rock Deformation: The Brittle Field , 1978 .

[2]  Kazushi Sato,et al.  Fracture Toughness Evaluation Based on Tension-softening Model and its Application to Hydraulic Fracturing , 2006 .

[3]  O. K. Mahabadi,et al.  Hybrid Finite-Discrete Element Simulation of the EDZ Formation and Mechanical Sealing Process Around a Microtunnel in Opalinus Clay , 2016, Rock Mechanics and Rock Engineering.

[4]  Kazushi Sato,et al.  Fracture Toughness Evaluation Based on Tension-softening Model and its Application to Hydraulic Fracturing , 2006 .

[5]  David Higdon,et al.  THE COYOTE UNIVERSE. II. COSMOLOGICAL MODELS AND PRECISION EMULATION OF THE NONLINEAR MATTER POWER SPECTRUM , 2009, 0902.0429.

[6]  Juergen H. Schön,et al.  Physical Properties of Rocks: A Workbook , 2011 .

[7]  H. Kolsky An Investigation of the Mechanical Properties of Materials at very High Rates of Loading , 1949 .

[8]  Donald Geman,et al.  Stochastic Relaxation, Gibbs Distributions, and the Bayesian Restoration of Images , 1984, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[9]  Richard J. Beckman,et al.  A Comparison of Three Methods for Selecting Values of Input Variables in the Analysis of Output From a Computer Code , 2000, Technometrics.

[10]  Stefan M. Wild,et al.  A Bayesian approach for parameter estimation and prediction using a computationally intensive model , 2014, Journal of Physics G: Nuclear and Particle Physics.

[11]  M. D. McKay,et al.  A comparison of three methods for selecting values of input variables in the analysis of output from a computer code , 2000 .

[12]  Q. Qian,et al.  Dynamic strength of rocks and physical nature of rock strength , 2009 .

[13]  Giovanni Grasselli,et al.  Continuum–discontinuum analysis of failure mechanisms around unsupported circular excavations in anisotropic clay shales , 2014 .

[14]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[15]  M. J. Bayarri,et al.  Predicting Vehicle Crashworthiness: Validation of Computer Models for Functional and Hierarchical Data , 2009 .

[16]  Brian Williams,et al.  A Bayesian calibration approach to the thermal problem , 2008 .

[17]  J. Gattiker,et al.  A novel approach for determining source–receptor relationships in model simulations: a case study of black carbon transport in northern hemisphere winter , 2013 .

[18]  Elastoplastic damage modeling the mechanical behavior of rock-like materials considering confining pressure dependency , 2013 .

[19]  Davide Elmo,et al.  An Integrated Numerical Modelling–Discrete Fracture Network Approach Applied to the Characterisation of Rock Mass Strength of Naturally Fractured Pillars , 2010 .

[20]  Emanuele Reccia,et al.  FEM-DEM Modeling for Out-of-plane Loaded Masonry Panels: A Limit Analysis Approach , 2012 .

[21]  J. Fredrich,et al.  Effect of grain size on brittle and semibrittle strength: Implications for micromechanical modelling of failure in compression , 1990 .

[22]  A. Munjiza The Combined Finite-Discrete Element Method , 2004 .

[23]  P. Sammonds,et al.  Damage of Saturated Rocks Undergoing Triaxial Deformation Using Complex Electrical Conductivity Measurements: Experimental Results , 1997 .

[24]  B. S. A. Tatone,et al.  A calibration procedure for two-dimensional laboratory-scale hybrid finite–discrete element simulations , 2015 .

[25]  T. Backers Fracture Toughness Determination and Micromechanics of Rock under Mode I and Mode II Loading , 2005 .

[26]  Thomas J. Santner,et al.  Design and analysis of computer experiments , 1998 .

[27]  H. Viswanathan,et al.  A mechanisms-based model for dynamic behavior and fracture of geomaterials , 2014 .

[28]  Tetrahedron to Points Contact Interaction in 3D , 2011 .

[29]  Pathegama Gamage Ranjith,et al.  The Effect of Specimen Size on Strength and Other Properties in Laboratory Testing of Rock and Rock-Like Cementitious Brittle Materials , 2011 .

[30]  R. H. Evans,et al.  Microcracking and stress-strain curves for concrete in tension , 1968 .

[31]  Sheng-Qi Yang,et al.  Experimental Investigation on Strength and Failure Behavior of Pre-cracked Marble Under Conventional Triaxial Compression , 2008 .

[32]  F. Donath,et al.  Dependence of Strain-Rate Effects on Deformation Mechanism and Rock Type , 1971, The Journal of Geology.

[33]  D. Higdon,et al.  Computer Model Calibration Using High-Dimensional Output , 2008 .

[34]  A. O'Hagan,et al.  Curve Fitting and Optimal Design for Prediction , 1978 .

[35]  A. O'Hagan,et al.  Bayesian calibration of computer models , 2001 .

[36]  I W Farmer,et al.  Principles of engineering geology , 1975 .

[37]  Jiansheng Xiang,et al.  Granular packing: numerical simulation and the characterisation of the effect of particle shape , 2009 .

[38]  A. Munjiza,et al.  Large Strain Finite Element Method: A Practical Course , 2015 .

[39]  Y. Obara,et al.  Fracture Toughness of Granite Measured Using Micro to Macro Scale Specimens , 2017 .

[40]  K. Xia,et al.  Dynamic rock tensile strengths of Laurentian granite: Experimental observation and micromechanical model , 2017 .

[41]  Antonio Munjiza,et al.  Validation of a three-dimensional Finite-Discrete Element Method using experimental results of the Split Hopkinson Pressure Bar test , 2014 .

[42]  Antonio Munjiza,et al.  Combined single and smeared crack model in combined finite-discrete element analysis , 1999 .

[43]  Christopher C. Pain,et al.  Modelling of massive particulates for breakwater engineering using coupled FEMDEM and CFD , 2008 .