Next Article in Journal
Optimizing Slag Content to Control Ds-Type Inclusions in 10B21 Cold Heading Steel
Next Article in Special Issue
Industry Scale Optimization: Hammer Crusher and DEM Simulations
Previous Article in Journal
Development of a Physical Separation Route for the Concentration of Base Metals from Old Wasted Printed Circuit Boards
Previous Article in Special Issue
Use of Discrete Element Modelling to Evaluate the Parameters of the Sampling Theory in the Feed Grade Sampler of a Sulphide Gold Plant
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Textural and Mineralogical Controls on Rock Strength Elucidated Using a Discrete Element Method Numerical Laboratory

1
Centre for Minerals Research, Department of Chemical Engineering, University of Cape Town, Cape Town 7701, South Africa
2
Julius Kruttschnitt Mineral Research Centre, Sustainable Minerals Institute, The University of Queensland, Brisbane, QLD 4072, Australia
*
Author to whom correspondence should be addressed.
Minerals 2021, 11(9), 1015; https://doi.org/10.3390/min11091015
Submission received: 19 August 2021 / Revised: 6 September 2021 / Accepted: 9 September 2021 / Published: 18 September 2021

Abstract

:
Numerical modelling techniques such as the discrete element method are now well established and extensively used in many applications including solid earth geoscience, materials science, geotechnical engineering and rock mechanics. The potential for this technique in understanding comminution mechanisms has been identified as highly promising. This work utilizes the discrete element method as a numerical laboratory to conduct investigations relevant to comminution that would otherwise be costly or time-consuming to perform in the field or laboratory. A benchmark numerical model for impact breakage of rock specimens is first established and validated against results of controlled laboratory experiments. Thereafter, the model is utilized to systematically investigate the potential dependency of ore breakage properties upon the prevalence of pre-existing fractures, as well as the mineralogical composition of the ore. These numerical experiments serve to highlight the potential for quantitatively relating the mechanical response of ore to its textural and mineralogical characteristics. Tandem utilization of numerical and laboratory experimentation to formulate and test hypotheses is a promising avenue to illuminate such relationships.

1. Introduction

The Discrete Element method (DEM) is a well-established numerical technique for simulating and studying the flow of granular materials as well as the fracture of rocks and rock-like materials [1,2,3]. A notable advantage is the ability to simulate at a desired scale for different applications [4,5]. Its application to serve as a “numerical laboratory” for predicting rock response or fracture, whilst varying micro-structural parameters and loading conditions, has been demonstrated and comprehensively tested in rock engineering and fracture mechanics [6,7,8,9]. The concept of the DEM as a numerical laboratory entails prescribing the boundary conditions and specifying material properties in the model, and then systemically varying these to investigate the variability of macroscopic measurable quantities [10]. This approach yields quantitative hypotheses that are amenable to subsequent experimental verification [11]. Areas of application where this has been demonstrated include materials science [12], rock mechanics [11,13,14,15], seismology [9] and mining sciences [16].
In minerals processing, the DEM has been demonstrated as a viable tool for simulating breakage of ore and equipment operation at known conditions [17,18,19,20]. The DEM has proven useful for gaining understanding of the underlying causes for operational challenges and to assist in the design of comminution equipment [21]. However, to date, the DEM has not been extensively applied as a numerical laboratory in comminution research, specifically to investigate the phenomenology of rock breakage and its dependence upon the textural characteristics of ore.
Research towards this end typically relies upon laboratory impact breakage experiments studies involving individual rock specimens [22]. The Short Impact Load Cell (SILC) has arisen as the preferred device for such studies [22]. The SILC is a hybrid of the drop weight tester and split Hopkinson bar [23] and is a shortened version of the ultra-fast load cell [23]. Figure 1a shows a schematic of the SILC apparatus, while Figure 1b displays typical force-time profiles obtained from laboratory experiments. It consists of a pneumatic gripper that clamps and releases a steel ball at a prescribed height above a vertically oriented cylindrical steel rod on which the rock specimen rests. The steel rod is instrumented with strain gauges to measure the changing forces acting on rock specimens during breakage events. Analysis of the force time-histories permits estimation of the rock stiffness, fracture force or tensile strength, and fracture energy [24,25,26].
SILC investigations have shown that factors such as rock shape [29] and size [30], pre-weakening and mineralogical composition [22] cause variations in the measured mechanical response during impact breakage experiments [22,27,30]. Rock shape markedly contributes to the variability of mechanical properties measured in SILC experiments, potentially masking the influence of intrinsic textural characteristics of the ore. This has prompted research employing specimens of controlled shape and size [31,32], as well as synthetic rock specimens of homogeneous mineral composition [30,33]. Such studies have confirmed that the mechanical properties measured by the SILC apparatus depend systematically upon the textural characteristics of the rock specimens. Furthermore, these experiments show that the peak fracture force as measured by the SILC is proportional to the tensile strength of the specimens [30,31,32], offering the possibility to harmonize geotechnical and comminution mechanical parameters.
Having demonstrated that mechanical properties and textural characteristics are inter-related [31], the next logical step is to illuminate the quantitative nature of this relationship. A purely empirical approach would involve extensive laboratory investigations in which the textural characteristics of the rock specimens are carefully measured, or synthetic specimens are prepared with strictly controlled internal structure. Such an experimental campaign is likely to be laborious and expensive, with the ever-present concern that any resulting relationships may apply only to the specimens examined.
The primary objective of the presented research is to demonstrate that the bonded particle DEM may be employed as a numerical laboratory to investigate textural controls on comminution that are difficult to access via purely laboratory-based investigations. Numerical simulations are conducted to identify potential relationships between mechanical properties and textural characteristics. Such numerically obtained relationships are treated as hypotheses, to be tested via subsequent controlled laboratory experiments. In the following, two such numerical hypotheses are presented, governing the relationship between peak fracture force and the density of microcracks or mineralogical composition within rock specimens. For a binary phase mineralogical texture, it is shown that fracture force decreases proportionately with an increase in the concentration of the weaker component. Fracture force similarly decreases linearly with an increase in the surface area of initial fractures within pre-weakened rock specimens.
The following section describes the numerical methodology. This entails the description of the numerical SILC apparatus, model validation as well as the construction of numerical rock specimens with prescribed fracture networks or mineralogical composition. The subsequent section presents results of numerical experiments involving these damaged or heterogeneous numerical specimens. The paper concludes with a discussion of the findings, guidelines for future research and experimental validation and the primary conclusions.

2. Materials and Methods

A bonded particle DEM model [34] is employed for the numerical experiments presented herein. Bonded particle DEM models represent rock specimens as an assembly of DEM spheres. Adjacent spheres are bound together via cylindrical linear-elastic beams. Deformation of this bonded assembly induces forces and moments between the bound DEM spheres, which are used to incrementally update the positions and orientations of the spheres. When the shear and normal stresses within an individual beam exceed a prescribed failure criterion, the beam interaction is replaced with a frictional elastic-repulsive interaction. This method mimics the spontaneous formation of frictional fracture surfaces during crack propagation within brittle materials.
Past research has demonstrated that such a formulation reproduces key phenomenology in fracture mechanics such as Mode I tensile crack propagation [35] and wing-crack formation [2,36]. This research employs ESyS-Particle [37], an open-source software package that implements a bonded particle DEM. A detailed mathematical description of the bonded particle formulation of ESyS-Particle can be found elsewhere [38,39]. The software has been successfully utilized to study the behaviour of rock and rock-like materials in the past [15,40,41,42,43,44].

2.1. Numerical Reproduction of the SILC Apparatus

A numerical reproduction of the SILC apparatus has been constructed, as illustrated in Figure 2. The steel ball is represented as a single, large DEM sphere, initially located a small distance above the DEM rock specimen and prescribed with a constant velocity commensurate with having fallen under gravitational acceleration from a user-specified height above the specimen.
The specimen itself is an assembly of non-overlapping DEM spheres filling a cylindrical volume, constructed using a volume-filling algorithm [45], as implemented in the GenGeo software that accompanies ESyS-Particle. The method involves firstly inserting non-overlapping spheres at random locations within the cylindrical volume, with radii in a specified range. Thereafter, four existing objects (spheres and/or volume boundaries) are selected and a new sphere is inserted to touch each of these, so long as the radius of the candidate sphere resides within the specified range. This sphere insertion procedure is repeated until the algorithm fails repeatedly to find a hole of sufficient size to insert new spheres. The resultant assembly consists of non-overlapping spheres that fill the prescribed volume with a final porosity that decreases as the range of sphere radii increases. For the simulations considered herein, the smallest spheres are one quarter the size of the largest spheres, resulting in a porosity of approximately 35%. The porosity of the DEM sphere packings is a major factor influencing the relationship between microscopic parameters governing interactions and the macroscopic mechanical response of bonded assemblies of DEM spheres. Consequently, the porosity of all numerical specimens is held fixed in this investigation.
The lower-most spheres comprising the specimen undergo elastic interactions with a planar substrate that remains stationary throughout simulations. Simple linear, elastic repulsive forces act between such spheres and the substrate. The net force acting on this substrate is recorded throughout simulations, producing force time-histories akin to that measured during laboratory SILC experiments.
At commencement of each simulation, all pairs of adjacent DEM spheres are identified and bound together via linear-elastic beam interactions. Four model parameters govern the mechanical properties of these beams and, by extension, the macroscopic mechanical properties of the specimen. These four parameters are: the elastic beam modulus (Yb), beam Poisson’s ratio (νb), cohesive strength (Cb) and tangent of the internal angle of friction of the beams ( tan θ ).
In fracture mechanics investigations using the bonded particle DEM, it is typical to conduct simulations of uniaxial compression tests and modify these four parameters by trial-and-error, to match the measured geotechnical properties (Young’s modulus, Yc, and unconfined compressive strength, σ c ) of a particular rock type. In this study, an alternative, parametric approach has been employed to calibrate the mechanical properties of numerical rock specimens. Each of the four model parameters were systematically varied to construct empirical calibration relationships between these parameters and measured macroscopic geotechnical properties. The details of this calibration procedure are provided in Appendix A and the so-obtained calibration relationships are as follows:
Y c = K 1 Y b ( 1 + K 2 v b )
σ c = K 3 C b K 4 + tan θ
These calibration relationships are valid only for DEM sphere assemblies constructed using the packing algorithm described above, with a ratio of sphere radii satisfying Rmin:Rmax = 1:4 and the diameter of the largest DEM sphere less than one tenth the length of the smallest linear dimension of the numerical rock specimen. The length of the cylindrical specimen is twice its diameter, ensuring that the well-known aspect-ratio dependency of mechanical properties is avoided. Adherence to these constraints ensures that all numerical specimens constructed and calibrated according to these relationships have statistically identical macroscopic geotechnical properties, regardless of the precise locations and sizes of individual spheres comprising the specimens.

2.2. Model Validation

The purpose of a numerical laboratory is to quantitatively predict measurable responses of a modelled system, under changing initial and boundary conditions. To validate such a numerical laboratory, it is insufficient to simply demonstrate that the model reproduces the results of an individual laboratory test case; the typical method employed to validate DEM models in the literature. The numerical laboratory must demonstrably predict the results of laboratory experiments across a broad range of experimental conditions.
The herein proposed SILC numerical laboratory is validated with reference to the experimental results of Barbosa et al. [30], who conducted laboratory SILC experiments using 3D-printed synthetic rock specimens. These authors manufactured synthetic rocks using silica-sand and jet binder and processed in a three-dimensional (3D) printer programmed with the desired geometry (cylinder) and size. Specimens were produced with 20% binder saturation and dried at 80 °C for about 24 h [33]. The length to diameter ratio (L:D) of all specimens varied between 1:1 to 2:1. Figure 3 shows a picture of different sizes of the 3D printed cylinders that were manufactured for the laboratory studies. A subset of the manufactured cylinders, with diameters of 10.2 mm, were subjected to uniaxial compression tests, yielding measurements of the Young’s modulus and uniaxial compressive strength of the synthetic specimens. Thereafter, the authors [30] conducted SILC experiments and reported the variation in peak fracture force, as a function of the diameter of the cylindrical specimens. These results comprise the target for the validation exercise.
Validation commences with calibration of the numerical model parameters, to match the measured geotechnical properties of the 3D-printed. The above calibration relationships (Equations (1) and (2)) are employed for this purpose and the obtained model parameters are summarized in Table 1. The geotechnical properties measured in the experiments, as well as the obtained values for the calibrated numerical specimens, are summarized in Table 2. According to this table, a comparison of the Young’s modulus and the uniaxial compressive strength from both experiment and numerical simulation shows that there is no statistically significant difference based on the calculated standard deviation. As such, the numerically calibrated parameters are considered a close match to the experiment.
The standard deviation of mechanical properties for the numerical simulations is significantly lower than that of the experiments, demonstrating that the numerical model is an “idealized” representation of rocks that have heterogeneous structure yet homogeneous mechanical properties. This is a distinct advantage when employing the bonded particle DEM as a numerical laboratory. Specimens may be constructed with statistically identical mechanical properties, avoiding the inherent uncertainties arising from the greater variability of rock specimens employed in laboratory studies.
Having calibrated the numerical model parameters, a suite of numerical SILC simulations was conducted, to measure the variation in peak fracture force for cylindrical specimens of differing diameters, in accordance with the experimental conditions employed by Barbosa et al. [30]. For each specimen size, five different numerical samples were fractured using the numerical SILC apparatus. Each numerical sample contained approximately 40,000 DEM spheres. The measured peak fracture force for each of these was then averaged and the standard deviation calculated. The results are illustrated in Figure 4, in comparison with the experimental results of Barbosa et al. [30]. A statistical T-test was conducted to test whether the experimental and numerical results reflect samples from statistically different distributions. For the 2-tailed T-test, a p-value of 0.035 confirms that the two sets of results are statistically identical.
The numerical laboratory matches not only the measured fracture force at the specimen size used for geotechnical characterization (10.2 mm) but also correctly predicts the observed change in peak fracture force as the specimen size is varied. Furthermore, visual inspection of the fragmentation obtained in the numerical simulations also closely matches that reported for the laboratory studies, as shown in Figure 5a,b. For both the laboratory experiments and the numerical simulations, the specimens are observed to cleave into two large fragments, with some small debris between the primary facture surfaces.
It is emphasized that no additional calibration of the numerical model was conducted in order to obtain agreement between the laboratory and numerical SILC results; calibration was only based on the geotechnical properties as reported by Barbosa et al. [30] for specimens with a diameter of 10.2 mm. This is considered a robust validation of the capacity of the numerical SILC model to quantitatively predict the mechanical response of brittle-elastic materials under a broad range of geometrical, initial and boundary conditions, a requisite feature for the model to be considered a numerical laboratory.

2.3. Construction of Pre-Weakened Numerical Rock Specimens

Pre-weakened numerical rock specimens are constructed using a geometrical method wherein pre-existing cracks are represented as a set of triangulated planes that are randomly distributed within the specimen volume. For this study, cracks are constrained to be square surfaces whose side length ranges between 0.8 mm and 2.0 mm. The orientations of these crack surfaces are randomized in all three coordinate directions, to represent pre-weakening due a variety of sequential, independent comminution and materials handling stages, as is typical of ore samples within a minerals processing circuit.
These crack surfaces are imposed as additional internal boundaries during construction of the DEM specimens. The sphere packing algorithm of GenGeo ensures that inserted spheres do not cross these surfaces, forming geometrically flat contacts between spheres either side of the crack surfaces. These DEM spheres are initially unbonded, undergoing only linear elastic repulsion and frictional interactions should they come into contact. This approach captures the essential features of pre-existing cracks: frictional surfaces that constrain deformation in their vicinity to be parallel and perpendicular to the fracture surface.
Illustrations of homogeneous and pre-weakened DEM rock specimens are shown in Figure 6. Figure 6a,b show the homogeneous specimen without micro-cracks and its bond networks, respectively, while Figure 6c shows a thin section from the bond network around the middle of the specimen in the Z-direction. Figure 6d and e show the specimen with micro-cracks and its bond networks, respectively, while Figure 6f shows a thin section from the bond network. The presence of some of the cracks denoted by holes within the bond network are highlighted with red ellipses.
For the numerical SILC experiments, pre-weakened specimens with differing levels of damage were constructed. This is achieved by prescribing the total number of pre-existing cracks, Nc = 0, 30, 90, 180 and 300. For each of these damage levels, 30 DEM rock specimens with different random crack distributions were constructed.

2.4. Construction of Numerical Rock Specimens with Mineralogical Composition

The bonded particle DEM approach is further adapted to investigate the effect of mineralogical composition of the numerical rock specimens. In this study, binary specimens comprised of two distinct mineral phases are considered; one phase representing a valuable mineral and the other gangue. This simplifies the construction and reduces the complexity of analyzing results obtained for more complex, multi-mineral textures.
The numerical rock specimens are constructed such that the mechanical properties of one of the two mineral phases may be weakened with respect to the other mineral phase. For the purposes of discussion, it is assumed that the softer component is the “valuable mineral” whilst the other phase is the “gangue”. The steps involved to build the numerical rock specimen with the desired mineralogy are as follows:
  • First generate a homogeneous rock specimen.
  • Specify random “seed points” within the homogeneous specimen which represent the centroids of individual mineral grains.
  • Group together spheres of the homogeneous rock specimen closest to each of these seed points; a Voronoi tessellation procedure.
  • Randomly assign each group of spheres to be either “valuable” or “gangue” in accordance with a prescribed probability, denoting the mineral grade of the specimen.
These four steps are illustrated in Figure 7. Figure 7a shows the initial homogeneous cylindrical specimen, while Figure 7b shows the insertion of seed points (coloured spheres). Grouping of DEM spheres into individual mineral grains, based upon their proximity to the nearest seed point, is shown in Figure 7c. Figure 7d shows the final two-phase specimen in which a percentage of the mineral grains are selected as representing the valuable mineral phase (red spheres) and gangue (light grey spheres). A 20% grade is illustrated in this case.
For this study, six different mineral grades are examined: 0% (homogeneous gangue), 10%, 25%, 50%, 75% and 100% (homogenous valuable). Figure 8a–f depicts the typical distribution of valuable minerals within the numerical rock specimens at different mineral grades. These were randomly distributed for each specimen.
With respect to mechanical properties of these DEM mineralogical rock specimens, two key assumptions are made:
  • One phase has mechanical properties twice the value of the other, i.e., the hard component is double the strength of the soft component. The weak component represents the valuable mineral whilst the strong component represents the gangue.
  • The mechanical properties of the interconnecting interactions between adjacent grains of differing phases are the same as those of the weaker component phase.

3. Results

This section describes the results of two numerical experiments, one which systematically explores the influence of pre-weakening in the DEM rock specimen and a second which investigates the mechanical response of the DEM rock specimen as the mineralogical composition changes.

3.1. Breakage Behaviour under the Influence of Pre-Existing Cracks

Breakage investigations were conducted on three random specimens at each pre-existing crack level to predict the fracture behaviours. Predictions of the force-time histories and fracture patterns are presented from Figure 9, Figure 10, Figure 11, Figure 12, Figure 13 and Figure 14.
An examination of the force-time histories in these Figures indicate that homogeneous rock specimens with 0 cracks (Figure 9a) and specimens with the least amount pre-existing cracks, i.e., 30 cracks (Figure 10a) exhibit the highest consistency of force-time histories. Histories are fairly consistent in the presence of 90 and 180 cracks as shown in Figure 11a and Figure 12a, respectively, while in the presence of 300 cracks (Figure 13a), the force-time history shows more variability.
The fracture patterns for three numerical specimens at each damage level are shown in Figure 9b–d to Figure 13b–d. The key observation here is also the consistency of the fracture pattern obtained at 0 crack level with specimens showing a tendency of splitting into almost two equal fragments (Figure 9b–d). A similar split of rock specimens into two fragments was observed in the presence of 30 cracks into the rock specimen, where some inconsistencies in the symmetry of the cleavage begin to appear. This can be seen with the fracture line of specimen 1 (Figure 10b), moving slightly more diagonally to the left-hand side whereas the fracture lines of specimens 2 and 3 tend to pass through the center of the specimen (Figure 10c,d). This behaviour can be attributed to the presence of micro-cracks driving the failure of crack propagation in variable directions within the rock specimen.
As the crack density intensifies, the irregularity of the fracture patterns within rock specimens become more apparent. These are shown with the fracture patterns at 90 cracks (Figure 11b–d), 180 cracks (Figure 12b–d) and 300 cracks (Figure 13b–d). Among these, fracture patterns at 180 cracks and 300 cracks show the highest irregularities. The irregularity of fracturing is further emphasized by examining the relative sizes of fragments produced. Figure 14 illustrates fragment masses for differing levels of damage. Relative fragment masses are calculated by identifying all DEM spheres forming a connected cluster and accumulating the masses of these spheres. It is the fragment mass relative to the original specimen. As the damage level increases, the number and range of masses of fragments increases markedly, as compared with pristine specimens (0% cracks) or mildly damaged specimens.
These results may also be examined with reference to the initial fracture surface area of pre-existing cracks, as illustrated in Figure 15. The mathematical expression for the calculation of the initial fracture surface area is given by Equation (3).
F S A = ( L m i n + L m a x 2 ) 2 N c
where, F S A = initial fracture surface area (mm2)
L m i n = minimum length of crack (mm)
L m a x = maximum length of crack (mm)
N c = number of cracks.
It is evident that the peak fracture force for the DEM SILC simulations decreases linearly as the initial fracture surface area increases. This represents a numerically derived hypothesis that is amenable to experimental verification.
Figure 15. Effect of fracture surface area on the mean fracture force. Error bar (SD) represents the standard deviation.
Figure 15. Effect of fracture surface area on the mean fracture force. Error bar (SD) represents the standard deviation.
Minerals 11 01015 g015

3.2. Mechanical Response under the Influence of Varying Mineralogical Composition

In the results that follow, the mineral grade is systematically varied and the resultant change in fracture force measured. Six different mineral grades are examined: 0% (homogeneous gangue), 10%, 25%, 50%, 75% and 100% (homogenous valuable). Thirty SILC simulations are conducted for randomly differing specimens for each level of mineral grade, to provide representative statistics on fracture force versus mineral grade.
Figure 16 shows a typical force-time history for the different grades. While the peak force decreases as the grade increases, it is further observed that there tends to be increasing oscillation of the force time-history in the presence of the valuable mineral. These oscillations are noticeable at even the lowest grade considered (10%). These oscillations are likely a consequence of increased incidence of elastic wave reflections at internal boundaries between the two mineral phases. The disappearance of such oscillations for specimens comprising 100% grade supports this conjecture.
Figure 17 shows the variation of the predicted fracture force as the grade increases. Homogeneous specimens with zero grade, i.e., without any inherent weak-phase mineral, have the highest fracture force of approximately 1.2 kN. The force decreases as the composition of the grade increases, reducing to 0.1 kN for specimens which are entirely composed of the weaker phase (i.e., 100% grade).
Preferential exposure (liberation) of mineral phases is also examined. Figure 18a–r summarize the locations and distributions of the valuable mineral, the fracture paths and locations within the rock specimen as well as the fragmentation obtained for different mineral compositions. At 0% grade (Figure 18a), the fracture path travels almost vertically downward (Figure 18b) with the specimen spitting into almost two halves (Figure 18c). For 10% grade (Figure 18d), fracture paths tend to flow around the valuable minerals (Figure 18e) which results in an irregular fragmentation (Figure 18f). At 25% and 50% grade composition (Figure 18g,j), the fracture paths and directions are controlled and directed along the densely populated regions of valuable minerals which give rise to tributaries within the fracture patterns (Figure 18h,k) and associated irregular fragmentation of the specimens (Figure 18i,l). However, at 75% and 100% grade (Figure 18m,p), when the rock specimen is predominantly composed of valuable minerals, there is no evidence for preferential fracture propagation. Instead, the fracture locations are randomly distributed (Figure 18n,q). This anomalous behaviour is interpreted as a deviation of the mechanical properties of the specimens from distinctly brittle-elastic to a more ductile nature. A lack of distinct through-going fractures, as well as localized indentation and minor fragmentation at the contact points (Figure 18o,r) are further evidence of the ductility of specimens comprised of predominantly weaker minerals.

4. Discussion

The primary objective of the research described herein is to demonstrate that the bonded particle DEM may be employed as a numerical laboratory to investigate textural controls on comminution. The approach involves first validating that the DEM model for a given experimental apparatus, the SILC in this instance, can quantitatively predict the results of laboratory control experiments. The validation exercise conducted herein demonstrates that the proposed DEM SILC model predicts the variation in peak fracture force as both the impact energy (drop height) and physical dimensions of specimens are varied. Qualitative analysis further demonstrates that the DEM SILC model predicts fracture patterns and fragmentation, in reasonable agreement with previously reported laboratory observations.
Having established the validity of the DEM SILC model, it is employed to conduct numerical control experiments that are difficult or tedious to conduct in the laboratory. This work focusses upon two such numerical control experiments, examining the variation in peak fracture force of cylindrical specimens with prescribed pre-existing cracks or differing binary mineral composition. Whilst simplified, these numerical control experiments yield quantitative predictions on the relationship between peak fracture force and textural characteristics of rock.

4.1. Breakage Behaviour under the Influence of Pre-Existing Cracks

The effect of pre-damage was investigated by specifying five levels of pre-damage within the numerical simulations. These pre-weakened rock specimens were simplified representations of rock specimens with inherent flaws or pre-weakening, with the added benefit of being consistent replications for control and statistical reliability through numerical experiments. Results indicated that the fracture patterns became increasingly irregular as the level of pre-damage increased, as more opportunities for the propagation of cracks occurred within specimens. From the perspective of fracture mechanics, this result is not unexpected. Stress concentrates at the tips of pre-existing cracks, forming preferential pathways for crack propagation that would be otherwise absent from a pristine rock specimen. These stress concentrations promote reactivation of pre-existing cracks rather than absorbing surface energy to form new fracture surfaces. Consequently, the required fracture force decreased as the initial fracture surface area of the specimens increased by virtue of the pre-damage. Results of the simulations highlighted this to be a linear relationship (Figure 15), which was not unexpected given that the elasticity of all specimens remained unchanged [34].
The reported linear relationship between peak fracture force and initial fracture surface area is a numerically derived hypothesis, amenable to laboratory verification. Such verification might be achieved through a variety of means. One approach might be to prepare 3D-printed specimens with prescribed internal fractures matching the fracture networks employed in the numerical experiments. Breakage of these synthetic, damaged specimens in the SILC will provide measurements of peak fracture force that may be directly compared with the numerical predictions. An alternative procedure for verifying this numerical hypothesis might involve controlled pre-weakening of rock specimens via compression or low-energy impact, non-destructive imaging of damaged specimens to estimate internal fracture surface area, followed by breakage in a SILC device. Such an experimental campaign should be sufficient to verify whether the reported linear relationship for peak fracture force versus fracture surface area is a reasonable hypothesis.

4.2. Breakage Behaviour under the Influence of Varying Mineralogical Composition

The second examination of the potential application of the proposed methodology was the consideration of mineralogy. For simplification, the construction of specimens specified two mineral phases, a representation of typical ore bodies known to exhibit such distinct behaviour between the mineral and gangue. The grade was systematically varied using mechanical properties to create specimens over five varieties spanning homogenous gangue to homogenous mineral.
The variability in force time histories and nature of fracture were found to be most consistent when closer to homogeneity of either form. Greater variability was observed in the fracture force as the grade increased from completely gangue to become binary, decreasing subsequently as the grade increased to 100% mineral. The fracture force decreased with an increase in ore grade, with fracture occurring predominately within regions in which the weaker mineral phase was present. This continued to occur until the numerical experiments with 100% grade, at which the change in mechanical properties from brittle-elastic to more ductile behaviour meant that specimens yielded readily under low force with no discernible fracture pattern.
The fracture of the rock specimen along the path of the mineral suggests that crack propagation is enhanced under the presence of a weaker phase. This leads to a more exposed surface area (liberation) of the weak phase. These findings demonstrate notable application for investigations of mineral liberation, in which the goal of exposing the highest surface area of the valuable material is the desired outcome. Prior work has examined the manner of preferential breakage through laboratory testing on common variabilities of ore [46]. It is worth highlighting that the main difficulty of these studies lies in the impracticability of being able to facilitate closely controlled experiments. This is mainly due to the inherent variability of specimens in internal characteristics such as texture and shape, as well as controlling the external environment they have been subject to and establishing a repeatable mode of testing. Such issues have been evident in the inconsistencies found in the breakage of ores even considered to be homogenous [47]. The methodology exhibited in this study can be utilized as a means of examining a wide range of effects associated with simple to complex mineralogical structures, with an added advantage of being able to consistently replicate both the specimen characteristics and the manner of impact. This would subsequently enable for the identification of specific scenarios to be closely studied in greater detail via experiments in which consistent specimens can be generated, e.g., by 3-d printing [33] or selective preparation of rock specimens whose mineralogical composition is measured via non-destructive imaging techniques.

5. Conclusions

In this work, DEM was utilized as a numerical laboratory to examine the effect of specimen pre-weakening and varying mineralogical composition on fracture behaviour. The methodology to calibrate BPM parameters to meet macroscopic mechanical properties was demonstrated to be capable of generating results consistent with standard geotechnical tests. It was further shown that empirical calibration relationships could be used to tune these across a range of properties characteristic of common rock types.
For the first application, the effect of pre-damage on the fracture force, nature and extent of fracture was investigated. An inverse linear relationship was exhibited between the fracture force and fracture surface area, highlighting the dominant role of pre-weakening in decreasing the amount of force required to cause fracture while increasing the resulting extent. This observation offers a hypothesis which warrants further scrutiny through further simulations coupled with closely controlled experiments which can investigate the optimal form of pre-damage which maximizes the efficiency of breakage in terms of energy expended to produce a given fracture surface area.
Secondly, this paper also assessed the potential of the numerical approach for studying the effect of ore texture on fracture behaviour. Ore texture is increasingly being considered as the next frontier in the design of flexible circuits of the future, the goal being to advantageously utilize information regarding mineralogy as a proxy to identify the most efficient manner of processing. Obtaining statistically representative information on the patterns of different textures is a crucial step toward this goal. Numerical experiments with a binary mineralogical composition indicated that fracture of specimens occurred mainly along the path of the weaker phase. This highlighted a promising opportunity to facilitate further studies on mineralogical characteristics that can be exploited to enhance liberation. While the texture considered was a simplified scenario, future simulations can consider more complex architectures and a range of parameters consistent with measurements from prior work. This can be supplemented by closely controlled laboratory tests to derive meaningful relationships that can be put to use toward quantitatively characterizing fracture using texture.

Author Contributions

T.O.: conceptualization, data curation, formal analysis, investigation, methodology, writing original draft; L.B.: Conceptualization; funding acquisition, project administration, resources, supervision, writing original draft, reviewing and editing manuscript; D.W.: conceptualization, data curation, methodology, supervision, reviewing and editing manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

Simulations were performed using High performance computers provided by the University of Cape Town’s ICTS Computing team (http://hpc.uct.ac.za). Details of the computational resources that were allocated are as follows: cluster type Dell C6420, 40 cores, 294GB RAM. For the research presented herein, each simulation employed 5 processor cores, with typical execution times of between four and five hours. Multiple simulations were conducted simultaneously for parametric studies, thus reducing the real time required to conduct the numerical experiments.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Calibration of BPM Parameters

In the rock mechanics literature, it is typical to calibrate the macro-mechanical properties of BPM models utilising numerical reproductions of standard geotechnical tests such as the uniaxial compression test (UCT) and either direct or indirect tension tests [48,49]. These model parameters are then applied to the DEM environment specific to the experimental apparatus under investigation. Following this, numerical results are compared against similar laboratory tests to verify that the selected model parameters produce results reasonably consistent with experimentation without the need for further tuning and refinement [49,50].
A range of calibration methodologies have been identified in literature, including trial and error, design of experiments and machine learning or neural networks approaches [51,52,53]. For this study, an approach called the Systematic Parametric Study (SPS) was used to calibrate BPM parameters against standard measures of macro-mechanical properties. The SPS approach, which takes a form similar to design of experiments, entails systematically varying model parameters against macroscopic mechanical properties. These are then regressed against a fitting function to describe an empirical relationship between the model parameters and macroscopic mechanical properties.
Four BPM parameters were calibrated: elastic beam modulus (Yb), beam Poisson’s ratio (νb), cohesive strength (Cb) and tangent of the internal angle of friction ( tan θ ). The first two parameters govern the elastic response of the DEM specimen, whilst the latter two govern the onset and nature of the brittle failure. UCT simulations were conducted to establish the relationship between these four parameters and DEM specimen properties such as the macroscopic Young’s modulus and unconfined compressive strength. The macroscopic elastic properties of a given specimen could subsequently be tuned to a desired value using the established calibration relationships.

Appendix A.1. UCS Simulations

A set of Uniaxial Compressive Strength (UCS) simulations were conducted on a numerical cylindrical-drill core specimen as shown in Figure A1, performed according to International Society of Rock Mechanics (ISRM) and American Society for Testing and Materials (ASTM) geotechnical standards [54,55,56,57]. A cylindrical geometry with an L:D ratio of 2:1 (20 mm:10 mm) was packed with bonded DEM-spheres with a radius ranging between 0.08 mm to 0.32 mm using the volume filling algorithm [45]. Planar walls were placed on at the top and bottom of the cylinder as shown in Figure 3 and prescribed with a downward speed from the top and an upward speed from the bottom of 0.125 m/s to impart compressive stress. The total instantaneous force on the DEM-spheres was recorded and the axial stress was calculated by dividing the instantaneous total force by the initial cross-sectional area of the specimen, in accordance with the standard laboratory procedure. The corresponding strain was determined by monitoring the mean vertical position at 25 and 75% of the height of the specimen as indicated in the Figure. This was calculated by dividing the relative vertical distance by the initial distance between these layers of DEM-spheres.
Figure A1. A typical DEM specimen employed for uniaxial compression calibration simulations.
Figure A1. A typical DEM specimen employed for uniaxial compression calibration simulations.
Minerals 11 01015 g0a1
Figure A2 illustrates a typical stress-strain curve obtained from the numerical UCT, which is observed to follow conventional experimental results. The Young’s modulus of the rock was determined by calculating the slope of the linear portion of the stress-strain curve within zone 1. The peak point between zones 2 and 3 of the graph determines the unconfined compressive strength of the specimen.
Figure A2. Axial stress vs. axial strain and curves for a typical uniaxial compression simulation with an illustration of the calculated Young’s modulus and unconfined compressive strength.
Figure A2. Axial stress vs. axial strain and curves for a typical uniaxial compression simulation with an illustration of the calculated Young’s modulus and unconfined compressive strength.
Minerals 11 01015 g0a2

Appendix A.2. Development of Empirical Calibration Relationships

During the UCT simulations, BPM parameters were paired and varied, whilst measuring the corresponding macroscopic mechanical properties. Herein, the beam elastic modulus and beam Poisson’s ratio were paired and varied against the measured Young’s modulus. The cohesive strength and tangent of the internal angle of friction were paired and varied against the measured unconfined compressive strength. Table A1 summarizes the ranges of the model parameters that were varied, the interval of variation for a successive increment, as well as the total number of levels of the varied individual model parameter.
Table A1. Range of values of model parameters utilised for the SPS.
Table A1. Range of values of model parameters utilised for the SPS.
Model ParameterRange of ValuesIntervalNumber of Levels
Elastic beam modulus (Yb)25–125 (GPa)25 GPa5
Beam Poisson’s ratio (νb)0.1–0.90.25
Cohesive strength (Cb)25–325 (MPa)50 MPa7
Tangent int. friction angle ( tan θ )0.1–2.10.55
Figure A3 shows a graph of the macroscopic Young’s modulus against elastic beam modulus at different Poisson’s ratios. It is observed that there is a linear increase in Yc as the Yb increases and an incremental vertical spread of the data points, i.e., different slopes. This suggests that the Yc is influenced by both Yb and vb. Similar trends can also be observed in the graph of UCS versus the Cb at different values of θ as depicted in Figure A4. The elastic beam parameters (Yb and vb) were found to have no influence on the UCS. Likewise, the beam strength properties (Cb and θ) did not influence the macroscopic Young’s modulus of the DEM specimens.
Figure A3. Macroscopic Young’s modulus under compressive loading (Yc) versus elastic beam modulus (Yb) of bonded DEM-spheres, for a range of beam Poisson’s ratios.
Figure A3. Macroscopic Young’s modulus under compressive loading (Yc) versus elastic beam modulus (Yb) of bonded DEM-spheres, for a range of beam Poisson’s ratios.
Minerals 11 01015 g0a3
Figure A4. Macroscopic UCS versus microscopic cohesive strength (Cb) of bonded DEM-spheres, for a range of beam bond internal angle of frictions.
Figure A4. Macroscopic UCS versus microscopic cohesive strength (Cb) of bonded DEM-spheres, for a range of beam bond internal angle of frictions.
Minerals 11 01015 g0a4
In order to isolate and study the effect of vb independently, Yc was divided by Yb and plotted against vb in Figure A3. This collapsed the data points into an almost linear (monotonic) trend. A similar technique was followed by plotting UCS/Cb against tan θ which resulted in a similar trend.
From inspecting the correlation between the BPM and macroscopic parameters, fitting functions in Equations (A1) and (A2) were proposed.
Y c = K 1 Y b ( 1 + K 2 v b )
σ c = K 3 C b K 4 + tan θ
Ki (i = 1, 2, 3 and 4) are the constants determined from the empirical fit. These reflect different elements of the spatial relationship of connecting adjacent DEM-spheres within the bonded rock specimens [35].
Non-linear regressions of the simulation data points were conducted using the SciPy optimize curve fit function [58] which returned the calculated parameters (K-values) and the coefficient of determination (R2) as summarized in Table A2. Observance of the fitted model and high R-squared suggested the correlations were reasonable to adopt.
Table A2. Summary of K values obtained from models fit.
Table A2. Summary of K values obtained from models fit.
K-ValueK-ValueR2
K1 = 0.568 ± 0.00187K2 = 0.236 ± 0.006120.9996
K3 = 1.850 ± 0.0470K4 = 1.46 ± 0.04910.9956
These constants were substituted into Equations (1) and (2) to establish the relationships between the BPM and macroscopic mechanical properties and selected as a basis for BPM calibration of rock specimens. Figure A5 and Figure A6 show the fitted empirical relationships against collapsed data points for the two scenarios, while Figure A7 and Figure A8 summarize the fitted relationships across the range of tests on a surface dot-plot.
Figure A5. Plot of Yc/Yb against vb fitted with Equation (A1).
Figure A5. Plot of Yc/Yb against vb fitted with Equation (A1).
Minerals 11 01015 g0a5
Figure A6. Plot of σ c /Cb against tan θ fitted with Equation (A2).
Figure A6. Plot of σ c /Cb against tan θ fitted with Equation (A2).
Minerals 11 01015 g0a6
Figure A7. Surface dot-plot and model fit of calibration relationship between measured Young’s modulus at different elastic beam modulus and beam Poisson’s ratio.
Figure A7. Surface dot-plot and model fit of calibration relationship between measured Young’s modulus at different elastic beam modulus and beam Poisson’s ratio.
Minerals 11 01015 g0a7
Figure A8. Surface dot-plot and model fit of calibration relationship between measured unconfined compressive strength at different beam cohesive strengths and internal angle of frictions.
Figure A8. Surface dot-plot and model fit of calibration relationship between measured unconfined compressive strength at different beam cohesive strengths and internal angle of frictions.
Minerals 11 01015 g0a8

References

  1. Cundall, P.A.; Strack, O.D. A discrete numerical model for granular assemblies. Geotechnique 1979, 24, 47–65. [Google Scholar] [CrossRef]
  2. Wang, Y.; Yin, X.; Ke, F.; Xia, M.; Peng, K. Numerical simulation of rock failure and earthquake process on mesoscopic scale. Pure Appl. Geophys. 2000, 157, 1905–1928. [Google Scholar] [CrossRef]
  3. Hentz, S.; Donzé, F.V.; Daudeville, L. Discrete element modelling of concrete submitted to dynamic loading at high strain rates. Comput. Struct. 2004, 82, 2509–2524. [Google Scholar] [CrossRef]
  4. Feng, Y.T.; Owen, D. Discrete element modelling of large scale particle systems—I: Exact scaling laws. Comput. Part. Mech. 2014, 1, 159–168. [Google Scholar] [CrossRef] [Green Version]
  5. Pöschel, T.; Saluena, C.; Schwager, T. Scaling properties of granular materials. In Continuous and Discontinuous Modelling of Cohesive-Frictional Materials; Springer: Berlin/Heidelberg, Germany, 2001; pp. 173–184. [Google Scholar]
  6. Procházka, P.P. Application of discrete element methods to fracture mechanics of rock bursts. Eng. Fract. Mech. 2004, 71, 601–618. [Google Scholar] [CrossRef]
  7. Invernizzi, S.; Lacidogna, G.; Carpinteri, A. Particle-based numerical modeling of AE statistics in disordered materials. Meccanica 2013, 48, 211–220. [Google Scholar] [CrossRef]
  8. Klejment, P.; Dębski, W. Crack nucleation in solid materials under external load-simulations with the Discrete Element Method. MATEC Web Conf. EDP Sci. 2018, 165, 22019. [Google Scholar] [CrossRef]
  9. Mair, K.; Abe, S. Breaking up: Comminution mechanisms in sheared simulated fault gouge. Pure Appl. Geophys. 2011, 168, 2277–2288. [Google Scholar] [CrossRef] [Green Version]
  10. Trent, B.C.; Margolin, L.G. A numerical laboratory for granular solids. Eng. Comput. 1992, 9, 191–197. [Google Scholar] [CrossRef] [Green Version]
  11. Han, Z.; Weatherley, D.; Puscasu, R. Projected area-based strength estimation for jointed rock masses in triaxial compression. Comput. Geotech. 2018, 104, 216–225. [Google Scholar] [CrossRef]
  12. Herrmann, H.; Eik, M.; Berg, V.; Puttonen, J. Phenomenological and numerical modelling of short fibre reinforced cementitious composites. Meccanica 2014, 49, 1985–2000. [Google Scholar] [CrossRef]
  13. Xu, Y.; Dai, F.; Xu, N.; Zhao, T.; Zhou, C. Discrete element simulation of dynamic semi-circular bend flexure tests of rocks using split Hopkinson pressure bar. Arab. J. Geosci. 2016, 9, 1–13. [Google Scholar] [CrossRef]
  14. Feng, P.; Dai, F.; Liu, Y.; Xu, N.; Zhao, T. Effects of strain rate on the mechanical and fracturing behaviors of rock-like specimens containing two unparallel fissures under uniaxial compression. Soil Dyn. Earthq. Eng. 2018, 110, 195–211. [Google Scholar] [CrossRef]
  15. Liu, Y.I.; Dai, F.; Xu, N.; Zhao, T.; Feng, P. Experimental and numerical investigation on the tensile fatigue properties of rocks using the cyclic flattened Brazilian disc method. Soil Dyn. Earthq. Eng. 2018, 105, 68–82. [Google Scholar] [CrossRef]
  16. Abe, S. Comparison of discrete element simulations to theoretical predictions of the elastic moduli of damaged rocks. Int. J. Rock Mech. Min. Sci. 2016, 88, 265–272. [Google Scholar] [CrossRef]
  17. Cleary, P.W. DEM simulation of industrial particle flows: Case studies of dragline excavators, mixing in tumblers and centrifugal mills. Powder Technol. 2000, 109, 83–104. [Google Scholar] [CrossRef]
  18. Cleary, P. Modelling comminution devices using DEM. Int. J. Numer. Anal. Methods Geomech. 2001, 25, 83–105. [Google Scholar] [CrossRef]
  19. Jiménez-Herrera, N.; Barrios, G.K.; Tavares, L.M. Comparison of breakage models in DEM in simulating impact on particle beds. Adv. Powder Technol. 2018, 24, 692–706. [Google Scholar] [CrossRef]
  20. Morrison, R.D.; Shi, F.; Whyte, R. Modelling of incremental rock breakage by impact–for use in DEM models. Miner. Eng. 2007, 20, 303–309. [Google Scholar] [CrossRef]
  21. Weerasekara, N.S.; Powell, M.S.; Cleary, P.W.; Tavares, L.M.; Evertsson, M.; Morrison, R.D.; Quist, J.; Carvalho, R.M. The contribution of DEM to the science of comminution. Powder Technol. 2013, 248, 3–24. [Google Scholar] [CrossRef]
  22. Lois-Morales, P.; Evans, C.; Bonfils, B.; Weatherley, D. The impact load cell as a tool to link comminution properties to geomechanical properties of rocks. Miner. Eng. 2020, 148, 106210. [Google Scholar] [CrossRef]
  23. Bourgeois, F.S.; Banini, G.A. A portable load cell for in-situ ore impact breakage testing. Int. J. Miner. Process. 2002, 65, 31–54. [Google Scholar] [CrossRef] [Green Version]
  24. King, R.P.; Bourgeois, F. Measurement of fracture energy during single-particle fracture. Miner. Eng. 1993, 6, 353–367. [Google Scholar] [CrossRef] [Green Version]
  25. Tavares, L.M.; King, R.P. Measurement of the load-deformation response from impact-breakage of particles. Int. J. Miner. Process. 2004, 74, S267–S277. [Google Scholar] [CrossRef]
  26. Tavares, L.M. Breakage of single particles: Quasi-static. Handb. Powder Technol. 2007, 12, 3–68. [Google Scholar]
  27. Bonfils, B. Quantifying of impact breakage of cylindrical rock particles on an impact load cell. Int. J. Miner. Process. 2017, 161, 1–6. [Google Scholar] [CrossRef] [Green Version]
  28. Tavares, L.M.; King, R.P. Single-particle fracture under impact loading. Int. J. Miner. Process. 1998, 54, 1–28. [Google Scholar] [CrossRef]
  29. Chandramohan, R.; Holtham, P.; Powell, M. The influence of particle shape in rock fracture. In XXV International Mineral Processing Congress; Australasian Institute of Mining and Metallurgy (AusIMM): Brisbane, Australia, 2010; pp. 6–10. [Google Scholar]
  30. Barbosa, K.; Chalaturnyk, R.; Bonfils, B.; Esterle, J.; Chen, Z. Testing Impact Load Cell Calculations of Material Fracture Toughness and Strength Using 3D-Printed Sandstone. Geotech. Geol. Eng. 2019, 38, 1065–1096. [Google Scholar] [CrossRef]
  31. Lois-Morales, P.; Bonfils, B.; Evans, C.; Weatherley, D. Quantifying the effect of Texture and Mineralogy over the physical properties of rocks using the Short Impact Load Cell. In Proceedings of the 14th International Mineral Processing Conference and 5th International Seminar on Geometallurgy, GECAMIN, Santiago, Chile, 28–30 November 2018; pp. 100–115. [Google Scholar]
  32. Barbosa, K.; Esterle, J.; Bonfils, B.; Chen, Z. The use of short impact load cell to derive geomechanical properties of sub-bituminous coal and mudstone. J. Nat. Gas Sci. Eng. 2019, 72, 103018. [Google Scholar] [CrossRef]
  33. Hodder, K. Fabrication, Characterization and Performance of 3D-Printed Sandstone Models. Ph.D. Thesis, University of Alberta, Edmonton, AB, Canada, 2017. [Google Scholar]
  34. Potyondy, D.O.; Cundall, P.A. A bonded-particle model for rock. Int. J. Rock Mech. Min. Sci. 2004, 41, 1329–1364. [Google Scholar] [CrossRef]
  35. Han, Z.; Weatherley, D.; Puscasu, R. A relationship between tensile strength and loading stress governing the onset of mode I crack propagation obtained via numerical investigations using a bonded particle model. Int. J. Numer. Anal. Methods Geomech. 2017, 41, 1979–1991. [Google Scholar] [CrossRef]
  36. Donze, F.V.; Bouchez, J.; Magnier, S.A. Modeling fractures in rock blasting. Int. J. Rock Mech. Min. Sci. 1997, 34, 1153–1163. [Google Scholar] [CrossRef]
  37. Esysparticle. Available online: https://launchpad.net/esys-particle (accessed on 5 August 2019).
  38. Wang, Y.; Abe, S.; Latham, S.; Mora, P. Implementation of particle-scale rotation in the 3-D lattice solid model. Pure Appl. Geophys. 2006, 163, 1769–1785. [Google Scholar] [CrossRef]
  39. Wang, Y.; Mora, P. The ESYS_particle: A new 3-D discrete element model with single particle rotation. In Advances in Geocomputing; Springer: Berlin/Heidelberg, Germany, 2009; pp. 183–228. [Google Scholar]
  40. Schöpfer, M.P.; Abe, S.; Childs, C.; Walsh, J.J. The impact of porosity and crack density on the elasticity, strength and friction of cohesive granular materials: Insights from DEM modelling. Int. J. Rock Mech. Min. Sci. 2009, 46, 250–261. [Google Scholar] [CrossRef] [Green Version]
  41. Abe, S.; Urai, J.L.; Kettermann, M. Fracture patterns in nonplane strain boudinage—Insights from 3-D discrete element models. J. Geophys. Res. Solid Earth 2013, 118, 1304–1315. [Google Scholar] [CrossRef] [Green Version]
  42. Zhao, T.; Crosta, G.B.; Dattola, G.; Utili, S. Dynamic fragmentation of jointed rock blocks during rockslide-avalanches: Insights from discrete element analyses. J. Geophys. Res. Solid Earth 2018, 123, 3250–3269. [Google Scholar] [CrossRef]
  43. Estay, D.; Chacana, F.; Ibarra, J.; Pérez, L.; Lascano, S. Bond calibration method for Young’s modulus determination in the discrete element method framework. Granul. Matter 2017, 19, 60. [Google Scholar] [CrossRef]
  44. Weatherley, D.; Boros, V.; Hancock, W.; Abe, S. Scaling benchmark of esys-particle for elastic wave propagation simulations. In Proceedings of the 2010 IEEE Sixth International Conference on e-Science, Brisbane, Australia, 7–10 December 2010. [Google Scholar]
  45. Gengeo. Available online: https://code.launchpad.net/esys-particle/gengeo (accessed on 5 August 2019).
  46. Little, L.; Mainza, A.N.; Becker, M.; Wiese, J.G. Using mineralogical and particle shape analysis to investigate enhanced mineral liberation through phase boundary fracture. Powder Technol. 2016, 301, 794–804. [Google Scholar] [CrossRef]
  47. Bbosa, L.; Powell, M.S.; Cloete, T.J. An investigation of impact breakage of rocks using the split Hopkinson pressure bar. J. S. Afr. Inst. Min. Metall. 2006, 106, 291–296. [Google Scholar]
  48. Cho, N.A.; Martin, C.D.; Sego, D.C. A clumped particle model for rock. Int. J. Rock Mech. Min. Sci. 2007, 44, 997–1010. [Google Scholar] [CrossRef]
  49. Yoon, J. Application of experimental design and optimization to PFC model calibration in uniaxial compression simulation. Int. J. Rock Mech. Min. Sci. 2007, 44, 871–889. [Google Scholar] [CrossRef]
  50. Kazerani, T.; Zhao, J. Micromechanical parameters in bonded particle method for modelling of brittle material failure. Int. J. Numer. Anal. Methods Geomech. 2010, 34, 1877–1895. [Google Scholar] [CrossRef]
  51. Zhang, Y. Probabilistic Calibration of a Discrete Particle Model. Ph.D. Thesis, Texas A&M University, College Station, TX, USA, 2010. [Google Scholar]
  52. Chehreghani, S.; Noaparast, M.; Rezai, B.; Shafaei, S.Z. Bonded-particle model calibration using response surface methodology. Particuology 2017, 32, 141–152. [Google Scholar] [CrossRef]
  53. Zhai, S.; Zhan, J.; Ba, Y.; Chen, J.; Li, Y.; Li, Z. PFC model parameter calibration using uniform experimental design and a deep learning network. In IOP Conference Series: Earth and Environmental Science; IOP Publishing: Bristol, UK, 2019; p. 032062. [Google Scholar]
  54. Bieniawski, Z.T.; Bernede, M.J. Suggested methods for determining the uniaxial compressive strength and deformability of rock materials: Part 1. Suggested method for determining deformability of rock materials in uniaxial compression. In International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts; Elsevier: Pergamon, Turkey, 1979; Volume 16, pp. 138–140. [Google Scholar]
  55. Hatheway, A.W. The complete ISRM suggested methods for rock characterization, testing and monitoring; 1974–2006. Environ. Eng. Geosci. 2009, 15, 47–48. [Google Scholar] [CrossRef]
  56. American Society for Testing and Materials. ASTM D2438-86/95—Standard Test Method for Unconfined Compressive Strength of Intact Rock Core Samples; American Society for Testing and Materials: West Conshohocken, PA, USA, 1991. [Google Scholar]
  57. American Society for Testing and Materials. ASTM D7012-13—Standard Test Methods for Compressive Strength and Elastic Moduli of Intact Rock Core Samples under Varying States of Stress and Temperatures; American Society for Testing and Materials: West Conshohocken, PA, USA, 2013. [Google Scholar]
  58. Virtanen, P.; Gommers, R.; Oliphant, T.E.; Haberland, M.; Reddy, T.; Cournapeau, D.; Burovski, E.; Peterson, P.; Weckesser, W.; Bright, J.; et al. SciPy 1.0: Fundamental algorithms for scientific computing in Python. Nat. Methods 2020, 17, 261–272. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. (a) A schematic of a SILC device [26]. (b) Typical force-time profiles obtained in SILC experiments [24,27,28].
Figure 1. (a) A schematic of a SILC device [26]. (b) Typical force-time profiles obtained in SILC experiments [24,27,28].
Minerals 11 01015 g001
Figure 2. A depiction of the DEM environment used for SILC tests.
Figure 2. A depiction of the DEM environment used for SILC tests.
Minerals 11 01015 g002
Figure 3. Pictures of 3D-printed cylinders [30].
Figure 3. Pictures of 3D-printed cylinders [30].
Minerals 11 01015 g003
Figure 4. DEM fracture force predicted at different specimen sizes compared against experiment.
Figure 4. DEM fracture force predicted at different specimen sizes compared against experiment.
Minerals 11 01015 g004
Figure 5. (a) Fragments generated from the SILC experiment [30]. (b) Fragments generated from the DEM simulation.
Figure 5. (a) Fragments generated from the SILC experiment [30]. (b) Fragments generated from the DEM simulation.
Minerals 11 01015 g005
Figure 6. (a) Homogeneous specimen without micro-cracks. (b) the bond network of the homogeneous specimen (c) a slice-cut of the bond network of the homogeneous specimen (d) specimen with micro-cracks. (e) the bond network of the specimen with micro-cracks (f) a slice-cut of the bond network of the specimen with micro-cracks.
Figure 6. (a) Homogeneous specimen without micro-cracks. (b) the bond network of the homogeneous specimen (c) a slice-cut of the bond network of the homogeneous specimen (d) specimen with micro-cracks. (e) the bond network of the specimen with micro-cracks (f) a slice-cut of the bond network of the specimen with micro-cracks.
Minerals 11 01015 g006
Figure 7. The four stages of creating specimens with a binary mineralogical texture. (a) homogeneous cylindrical specimen (b) homogenous specimen with seed points (coloured spheres) (c) DEM-spheres close to each seed points are grouped together (d) grouped spheres assigned to be either “valuable” (red spheres) or “gangue” (grey spheres).
Figure 7. The four stages of creating specimens with a binary mineralogical texture. (a) homogeneous cylindrical specimen (b) homogenous specimen with seed points (coloured spheres) (c) DEM-spheres close to each seed points are grouped together (d) grouped spheres assigned to be either “valuable” (red spheres) or “gangue” (grey spheres).
Minerals 11 01015 g007
Figure 8. Typical mineral grades of rock specimen. (a) 0%, (b) 10%, (c) 25%, (d) 50%, (e) 75% and (f) 100% grade composition.
Figure 8. Typical mineral grades of rock specimen. (a) 0%, (b) 10%, (c) 25%, (d) 50%, (e) 75% and (f) 100% grade composition.
Minerals 11 01015 g008
Figure 9. (a) Force-time histories of three rock specimens at 0-crack level and (bd) their respective post-process analyses showing the fracture pattern.
Figure 9. (a) Force-time histories of three rock specimens at 0-crack level and (bd) their respective post-process analyses showing the fracture pattern.
Minerals 11 01015 g009
Figure 10. (a) Force-time histories of three rock specimens at 30-crack level and (bd) their respective post-process analyses showing the fracture pattern.
Figure 10. (a) Force-time histories of three rock specimens at 30-crack level and (bd) their respective post-process analyses showing the fracture pattern.
Minerals 11 01015 g010
Figure 11. (a) Force-time histories of three rock specimens at 90-crack level and (bd) their respective post-process analyses showing the fracture pattern.
Figure 11. (a) Force-time histories of three rock specimens at 90-crack level and (bd) their respective post-process analyses showing the fracture pattern.
Minerals 11 01015 g011
Figure 12. (a) Force-time histories of three rock specimens at 180-crack level and (bd) their respective post-process analyses showing the fracture pattern.
Figure 12. (a) Force-time histories of three rock specimens at 180-crack level and (bd) their respective post-process analyses showing the fracture pattern.
Minerals 11 01015 g012
Figure 13. (a) Force-time histories of three rock specimens at 300-crack level and (bd) their respective post-process analyses showing the fracture pattern.
Figure 13. (a) Force-time histories of three rock specimens at 300-crack level and (bd) their respective post-process analyses showing the fracture pattern.
Minerals 11 01015 g013
Figure 14. Ranges of fragments generated at the different crack levels. (a) fine and big fragments at 0 crack, (b) fine and big fragments at 30 cracks, (c) fine, medium size and big fragments at 90 cracks, (d) fine, medium size and big fragments at 180 cracks, (e) fine, medium size and big fragments at 300 cracks. In each case, the relative fragment mass is the of fragment mass relative to the original specimen.
Figure 14. Ranges of fragments generated at the different crack levels. (a) fine and big fragments at 0 crack, (b) fine and big fragments at 30 cracks, (c) fine, medium size and big fragments at 90 cracks, (d) fine, medium size and big fragments at 180 cracks, (e) fine, medium size and big fragments at 300 cracks. In each case, the relative fragment mass is the of fragment mass relative to the original specimen.
Minerals 11 01015 g014aMinerals 11 01015 g014b
Figure 16. Force-time history of the representative specimen at different percentage grade composition.
Figure 16. Force-time history of the representative specimen at different percentage grade composition.
Minerals 11 01015 g016
Figure 17. Variation of peak fracture force versus mineral grade.
Figure 17. Variation of peak fracture force versus mineral grade.
Minerals 11 01015 g017
Figure 18. Locations and distributions of valuables within rock specimens, fracture locations and paths within rock specimens, fragmentation of specimens. (a) rock specimen with 0% mineral grade, (b) Fracture locations and patterns at 0% mineral grade, (c) Fragment mass at 0% mineral grade, (d) rock specimen with 10% mineral grade, (e) Fracture locations and patterns at 10% mineral grade, (f) Fragment mass at 10% mineral grade, (g) rock specimen with 25% mineral grade, (h) Fracture locations and patterns at 25% mineral grade, (i) Fragment mass at 25% mineral grade, (j) rock specimen with 50% mineral grade, (k) Fracture locations and patterns at 50% mineral grade, (l) Fragment mass at 50% mineral grade, (m) rock specimen with 75% mineral grade, (n) Fracture locations and patterns at 75% mineral grade, (o) Fragment mass at 75% mineral grade, (p) rock specimen with 100% mineral grade, (q) Fracture locations and patterns at 100% mineral grade and (r) Fragment mass at 100% mineral grade.
Figure 18. Locations and distributions of valuables within rock specimens, fracture locations and paths within rock specimens, fragmentation of specimens. (a) rock specimen with 0% mineral grade, (b) Fracture locations and patterns at 0% mineral grade, (c) Fragment mass at 0% mineral grade, (d) rock specimen with 10% mineral grade, (e) Fracture locations and patterns at 10% mineral grade, (f) Fragment mass at 10% mineral grade, (g) rock specimen with 25% mineral grade, (h) Fracture locations and patterns at 25% mineral grade, (i) Fragment mass at 25% mineral grade, (j) rock specimen with 50% mineral grade, (k) Fracture locations and patterns at 50% mineral grade, (l) Fragment mass at 50% mineral grade, (m) rock specimen with 75% mineral grade, (n) Fracture locations and patterns at 75% mineral grade, (o) Fragment mass at 75% mineral grade, (p) rock specimen with 100% mineral grade, (q) Fracture locations and patterns at 100% mineral grade and (r) Fragment mass at 100% mineral grade.
Minerals 11 01015 g018
Table 1. Model parameters used to configure the DEM environment for simulated impact breakage of a SILC.
Table 1. Model parameters used to configure the DEM environment for simulated impact breakage of a SILC.
ParameterSymbolValue(s)
Minimum DEM-sphere radiusRmin0.08 mm
Maximum DEM-sphere radiusRmax0.32 mm
Steel ball radiusRs25.44 mm
Elastic beam modulus Y b 3.35 GPa
Cohesive strength C b 37.2 MPa
Internal angle of friction θ 45°
Coefficient of sliding frictionµf0.6
Beam Poisson’s ratio v b 0.25
Steel ball elastic modulusE360 GPa
Table 2. Geotechnical parameters measured for the laboratory experiments and the calibrated numerical model.
Table 2. Geotechnical parameters measured for the laboratory experiments and the calibrated numerical model.
ExperimentDEM
ParameterValueSDValueSD
Y c 1.8 GPa0.231.8 GPa0.10
σ c 27.9 MPa5.2922.35 MPa1.83
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Oladele, T.; Bbosa, L.; Weatherley, D. Textural and Mineralogical Controls on Rock Strength Elucidated Using a Discrete Element Method Numerical Laboratory. Minerals 2021, 11, 1015. https://doi.org/10.3390/min11091015

AMA Style

Oladele T, Bbosa L, Weatherley D. Textural and Mineralogical Controls on Rock Strength Elucidated Using a Discrete Element Method Numerical Laboratory. Minerals. 2021; 11(9):1015. https://doi.org/10.3390/min11091015

Chicago/Turabian Style

Oladele, Temitope, Lawrence Bbosa, and Dion Weatherley. 2021. "Textural and Mineralogical Controls on Rock Strength Elucidated Using a Discrete Element Method Numerical Laboratory" Minerals 11, no. 9: 1015. https://doi.org/10.3390/min11091015

APA Style

Oladele, T., Bbosa, L., & Weatherley, D. (2021). Textural and Mineralogical Controls on Rock Strength Elucidated Using a Discrete Element Method Numerical Laboratory. Minerals, 11(9), 1015. https://doi.org/10.3390/min11091015

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop