Paper contribution

Resolution-induced uncertainty in DEM-predicted brittle damage thresholds in sandstone

Proceedings
EUROCK 2026
Pages
645-650
Submitted
30 April 2026
Accepted
15 June 2026
Published
12 August 2026

Summary

Abstract

Brittle damage thresholds predicted by discrete element method (DEM) modelling are informative for interpreting rock failure processes and supporting engineering design. Nonetheless, these predictions are constrained by the model particle resolution. This study investigates whether variations in model particle resolution affect the interpretation of damage thresholds in sandstone subjected to uniaxial compression. A recently developed Subspring Network Breakable Voronoi model was used to simulate uniaxial compression at four different particle resolutions. The results show a systematic dependence of the predicted response on model particle resolution. With increasing resolution, the predicted uniaxial compressive strength, crack initiation, and crack damage thresholds generally increase, while Young’s modulus decreases slightly and Poisson’s ratio shows limited sensitivity. Finer resolutions also produce significantly higher crack counts and denser fracture patterns. Despite these variations, the dominant macroscopic failure mode remains tensile dominated; axial splitting is observed across all resolutions. These findings indicate that model particle resolution influences quantitative interpretations of brittle damage more strongly than qualitative identification of failure modes. Consequently, model particle resolution should not be treated merely as a numerical setting but rather as a source of modelling uncertainty in DEM based rock engineering analyses. An intermediate model particle resolution is considered to provide a practical compromise between computational efficiency and representative failure behaviour.

Contributors

Authors

  1. Jinci Chen Institute of GeoEnergy Engineering, School of Energy, Geoscience, Infrastructure and Society, Heriot-Watt University, Edinburgh, UK; British Geological Survey, Keyworth, Nottingham, UK ORCID 0009-0002-0920-9778
  2. Elli-Maria Charalampidou Institute of GeoEnergy Engineering, School of Energy, Geoscience, Infrastructure and Society, Heriot-Watt University, Edinburgh, UK
  3. Ali Ozel School of Engineering and Physical Sciences, Heriot-Watt University, Edinburgh, UK
  4. James Woodman Geo-Design Consulting Engineers Ltd, Ringwood, UK
  5. Audrey Ougier-Simonin British Geological Survey, Keyworth, Nottingham, UK

Index terms

Keywords

Full text

Attachments

Works cited

References

  1. Itasca Consulting Group, Inc. (2025) Subspring network model. PFC Documentation. Available at: https://docs.itascacg.com/itasca940/common/contactmodel/rbsn/doc/manual/cmspringnetwork.html?node5092. Accessed 25 May 2026.
  2. Li W, Yu L, Zhang T, Su H, Mi X, Fan D, Jin B (2024) Quantitative analysis of grain size effect on tensile mechanical behavior of granite based on multi-level force chain networks. Computational Particle Mechanics, 11(5): 2245–2266.DOI:10.1007/s40571-024-00790-3
  3. Nicksiar M, Martin C (2012)Evaluation of methods for determining crack initiation in compression tests on low-porosity rocks.Rock Mechanics and Rock Engineering, 45(4): 607–617.DOI: 10.1007/s00603-012-0221-6
  4. Paterson MS, Wong TF (2005)Experimental Rock Deformation – The Brittle Field.2nd edn.Springer,Berlin Heidelberg.
  5. Peng J, Wong LNY, Teh CI (2021) Influence of grain size on strength of polymineralic crystalline rock:Newinsights from DEM grain-based modeling. Journal of Rock Mechanics and Geotechnical Engineering, 13(4): 755–766.DOI:10.1016/j.jrmge.2021.01.011
  6. Potyondy DO, Cundall PA (2004) A bonded-particle model for rock. International Journal of Rock Mechanicsand Mining Sciences, 41(8): 1329–1364. DOI:10.1016/j.ijrmms.2004.09.011
  7. Potyondy DO (2012)A flat-jointed bonded-particle material for hard rock.46th US Rock Mechanics/Geomechanics Symposium,Chicago, Illinois, USA,24–27 June 2012,Paper No. ARMA-2012-501.
  8. Potyondy DO, Fu W (2024)A 3D subspring network breakable Voronoi model for rock: laboratory-scale behavior.58th US Rock Mechanics/Geomechanics Symposium,Golden, Colorado, USA,23–26 June 2024,Paper No. D031S029R005.
  9. Potyondy DO, Fu W, Purvance M, Mas Ivars D (2025) A 3D subspring network breakable Voronoi model for rock. Rock Mechanics and Rock Engineering. DOI: 10.1007/s00603-025-04655-w
  10. Rasmussen LL (2021) Hybrid lattice/discrete element method for bonded block modeling of rocks. Computers and Geotechnics, 130: 103907. DOI:10.1016/j.compgeo.2020.103907
  11. Rasmussen LL, de Farias MM, de Assis AP (2018) Extended Rigid Body Spring Network method for the simulation of brittle rocks. Computers and Geotechnics, 99: 31–41.DOI:10.1016/j.compgeo.2018.02.021
  12. Sun Y, Kwok CY, Duan K (2024) Size effects on crystalline rock masses: Insights from grain-based DEM modeling. Computers and Geotechnics, 171: 106376. DOI:10.1016/j.compgeo.2024.106376
  13. Woodman J (2020) Thermo-mechanical loading of intact rock and discontinuities. PhD thesis, University of Leeds. Available at: https://etheses.whiterose.ac.uk/id/eprint/26306/. Accessed 25 May 2026.
  14. Woodman J, Ougier-Simonin A, Stavrou A, Vazaios I, Murphy W, Thomas ME, Reeves HJ (2021) Laboratory experiments and grain based discrete element numerical simulations investigating the thermo-mechanical behaviour of sandstone. Geotechnical and Geological Engineering, 39(7): 4795–4815.DOI:10.1007/s10706-021-01794-z
  15. Zhang T, Yu L, Peng Y, Ju M, Yin Q, Wei J, Jia S (2022) Influence of grain size and basic element size on rock mechanical characteristics: insights from grain-based numerical analysis. Bulletin of Engineering Geology and the Environment, 81(9): 347.DOI:10.1007/s10064-022-02856-z