Have a personal or library account? Click to login
Strength of industrial sandstones modelled with the Discrete Element Method Cover

Strength of industrial sandstones modelled with the Discrete Element Method

Open Access
|Oct 2021

Figures & Tables

Figure 1

The scheme of DEM simulation (after O’Sullivan 2011).
The scheme of DEM simulation (after O’Sullivan 2011).

Figure 2

The interaction of the two particles by the force–displacement law. The figure is redrawn from Ferdowsi (2014). Particles a and b of radii Ra and Rb are at a distance D with overlap C. Pa and Pb are the points of intersection of the line connecting the disc centres with the boundaries of the discs; d and t are unit vectors. Derivatives of r and θ stand linear damping and rotational damping, respectively.
The interaction of the two particles by the force–displacement law. The figure is redrawn from Ferdowsi (2014). Particles a and b of radii Ra and Rb are at a distance D with overlap C. Pa and Pb are the points of intersection of the line connecting the disc centres with the boundaries of the discs; d and t are unit vectors. Derivatives of r and θ stand linear damping and rotational damping, respectively.

Figure 3

Forces and moments between the particles bonded through rotational elastic-brittle bonds (Abe et al. 2014).
Forces and moments between the particles bonded through rotational elastic-brittle bonds (Abe et al. 2014).

Figure 4

Location of the active quarries on the map of Poland from which the tested sandstones originated.
Location of the active quarries on the map of Poland from which the tested sandstones originated.

Figure 5

Schematic representation of the uniaxial compression test. Sample is compressed from the top and bottom until failure. On the left – simplified scheme of the test, on the right – laboratory experiment.
Schematic representation of the uniaxial compression test. Sample is compressed from the top and bottom until failure. On the left – simplified scheme of the test, on the right – laboratory experiment.

Figure 6

Cracked samples after uniaxial compression test.
Cracked samples after uniaxial compression test.

Figure 7

Schematic representation of different possible failure patterns of the sample under uniaxial compression extracted and proposed by Basu et al. (2013).
Schematic representation of different possible failure patterns of the sample under uniaxial compression extracted and proposed by Basu et al. (2013).

Figure 8

Comparison between numerical and laboratory stress-strain curves for selected samples: (A) Brenna BR2, (B) Mucharz MU1, (C) Radków RAD1 and (D) Tumlin TU2.
Comparison between numerical and laboratory stress-strain curves for selected samples: (A) Brenna BR2, (B) Mucharz MU1, (C) Radków RAD1 and (D) Tumlin TU2.

Figure 9

Scaled number of bonds Nbonds/Nbonds max (number of bonds divided by maximum number of bonds) inside the sample as a function of axial displacement for all four types of sandstones.
Scaled number of bonds Nbonds/Nbonds max (number of bonds divided by maximum number of bonds) inside the sample as a function of axial displacement for all four types of sandstones.

Figure 10

Scaled number of fractures (Nfrac/Nfrac max-number of fractures divided by maximum number of fractures) as a function of axial strain with the stress–strain curve in the background: (A) Brenna, (B) Mucharz, (C) Radków, (D) Tumlin.
Scaled number of fractures (Nfrac/Nfrac max-number of fractures divided by maximum number of fractures) as a function of axial strain with the stress–strain curve in the background: (A) Brenna, (B) Mucharz, (C) Radków, (D) Tumlin.

Figure 11

(A) Scaled potential energy of bonds Epot/Epot max and (B) scaled kinetic energy Ekin/Ekin max for all four types of sandstones.
(A) Scaled potential energy of bonds Epot/Epot max and (B) scaled kinetic energy Ekin/Ekin max for all four types of sandstones.

Figure 12

Comparison between two scaled different components of kinetic energy (Eckin/Eckin max – kinetic energy divided by maximum kinetic energy) plotted as a function of strain: (A) Brenna, (B) Mucharz, (C) Radków, (D) Tumlin.
Comparison between two scaled different components of kinetic energy (Eckin/Eckin max – kinetic energy divided by maximum kinetic energy) plotted as a function of strain: (A) Brenna, (B) Mucharz, (C) Radków, (D) Tumlin.

Figure 13

Four different components of scaled potential energy of bonds (Ecpot/Ecpot max – potential energy of bonds divided by maximum potential energy of bonds) plotted together: (A) Brenna, (B) Mucharz, (C) Radków, (D) Tumlin.
Four different components of scaled potential energy of bonds (Ecpot/Ecpot max – potential energy of bonds divided by maximum potential energy of bonds) plotted together: (A) Brenna, (B) Mucharz, (C) Radków, (D) Tumlin.

Figure 14

Exemplary presentation of the numerical slab after two different tests: (A) failure resistance test, (B) impact resistance test. View of bonds.
Exemplary presentation of the numerical slab after two different tests: (A) failure resistance test, (B) impact resistance test. View of bonds.

Figure 15

Schematic representation of all four tests: (A) failure resistance test, (B) impact resistance test, (C) vibration resistance test and (D) abrasion resistance test. View of the cross-section of the numerical slab.
Schematic representation of all four tests: (A) failure resistance test, (B) impact resistance test, (C) vibration resistance test and (D) abrasion resistance test. View of the cross-section of the numerical slab.

Figure 16

Results of all four tests: (A) failure resistance test, dependence between scaled maximum force F/Fmax (force divided by maximum force) and sample thickness; (B) impact resistance test, dependence between scaled number of fractures Nfrac/Nfrac max (number of fractures divided by maximum number of fractures) which appeared in the sample after the test and sample thickness; (C) vibration resistance test, number of fragments into which the sample has fallen apart as a function of sample thickness; (D) abrasion resistance test, bulk friction coefficient of the sandstone.
Results of all four tests: (A) failure resistance test, dependence between scaled maximum force F/Fmax (force divided by maximum force) and sample thickness; (B) impact resistance test, dependence between scaled number of fractures Nfrac/Nfrac max (number of fractures divided by maximum number of fractures) which appeared in the sample after the test and sample thickness; (C) vibration resistance test, number of fragments into which the sample has fallen apart as a function of sample thickness; (D) abrasion resistance test, bulk friction coefficient of the sandstone.

Figure 17

Network of fractures which appeared in the sample after impact resistance test. View of bonds.
Network of fractures which appeared in the sample after impact resistance test. View of bonds.

Number of grains after failure_

Total number of grainsNumber of the smallest grainsNumber of middle grainsNumber of the biggest grains
Brenna19,88819,87981
Mucharz23,38423,37743
Radków3379033,7641412
Tumlin152,953152,94481

Comparison of failure patterns in DEM modeling (initial model, and model just after the failure - view of bonds between particles) and in laboratory (laboratory failure)_

TypeInitial numerical modelNumerical failureLaboratory failure
Brenna sample BR2
Mucharz sample MU1
Radków sample RAD1
Tumlin sample TU2

Range of results of geomechanical laboratory tests (mean values in brackets)_

Lithological typeρ (kg/m3)UCS (MPa)E (GPa)ν (−)
Gogula sandstones from Brenna2391 – 2444 (2429)99 – 107 (102)20.2 – 21.8 (20.8)0.21 – 0.25 (0.22)
Krosno sandstones from Mucharz2627 – 2668 (2655)138 – 153 (143)26.6 – 31.8 (29.2)0.18 – 0.23 (0.19)
Joint sandstones from Radków2119 – 2172 (2142)55 – 60 (58)22.5 – 24.9 (23.4)0.20 – 0.29 (0.24)
Sandstones from Tumlin2355 – 2480 (2448)115 – 128 (123)26.2 – 29.4 (27.4)0.22 – 0.28 (0.25)

Detailed parameters of numerical models of sandstones_

Sandstone 1 BrennaSandstone 2 MucharzSandstone 3 RadkówSandstone 4 Tumlin
Sample geometryCylinder:height – 75 mm;radius – 18.75 mmCylinder:height – 75 mm;radius – 18.75 mmCylinder:height – 75 mm;radius – 18.75 mmCylinder:height – 75 mm;radius – 18.75 mm
Particles radii0.25 – 2.0 mm0.25 – 0.5 mm0.25 – 1.0 mm0.125 – 0.5 mm
Particles density2429 kg/m32655 kg/m32142 kg/m32448 kg/m3
Initial number of particles82,266339,324150,6531,155,851
Initial number of bonds297,847821,488500,0243,831,657
Time step1.23×10−5s8.32×10−6s1.06×10−5s5.14×10−6s
Parameters of bonds :- Young's modulus Eb (MPa)- Poisson's ratio vb (dimensionless)- cohesion cb (MPa)- tangent of internal friction angle φb (dimensionless)Eb = 5 284 MPavb = 0.25cb = 13.59 MPatan φb = 1.0Eb= 50 258 MPavb= 0.25cb = 205.33 MPatan φb = 1.0Eb = 12 492 MPavb = 0.25cb = 21.19 MPatan φb = 1.0Eb = 15 144 MPavb = 0.25cb = 44.7 MPatan φb = 1.0
Stroke rate0.2 mm/s0.2 mm/s0.2 mm/s0.2 mm/s
Computational time287 min667 min335 min1440 min
Number of processors9222893320

Comparison between crucial macroscopic parameters used for model calibration_ ‘Sim’ stands for simulation result and ‘lab’ stands for result from the laboratory experiment_

Selected testE (GPa)UCS (MPa)ez cr (%)ν (−)
Brenna sample BR2Sim: 20.8Lab: 20.5Sim: 98Lab: 99Sim: 0.51Lab: 0.52Sim: 0.21Lab: 0.25
Mucharz sample MU1Sim: 29.2Lab: 31.8Sim: 136Lab: 138Sim: 0.47Lab: 0.49Sim: 0.19Lab: 0.23
Radków sample RAD1Sim: 23.4Lab: 22.8Sim: 61Lab: 60Sim: 0.28Lab: 0.29Sim: 0.21Lab: 0.20
Tumlin sample TU2Sim: 27.4Lab: 26.2Sim: 118Lab: 115Sim: 0.46Lab: 0.51Sim: 0.26Lab: 0.22
DOI: https://doi.org/10.2478/sgem-2021-0020 | Journal eISSN: 2083-831X | Journal ISSN: 0137-6365
Language: English
Page range: 346 - 365
Submitted on: Dec 31, 2020
Accepted on: Jul 8, 2021
Published on: Oct 9, 2021
Published by: Wroclaw University of Science and Technology
In partnership with: Paradigm Publishing Services
Publication frequency: 4 issues per year

© 2021 Piotr Klejment, Robert Dziedziczak, Paweł Łukaszewski, published by Wroclaw University of Science and Technology
This work is licensed under the Creative Commons Attribution 4.0 License.