Research Article  Open Access
Jiuqun Zou, Jihuan Han, Weihao Yang, "Investigating the Influences of Indentation Hardness and Brittleness of RockLike Material on Its Mechanical Crushing Behaviors", Mathematical Problems in Engineering, vol. 2020, Article ID 4713532, 16 pages, 2020. https://doi.org/10.1155/2020/4713532
Investigating the Influences of Indentation Hardness and Brittleness of RockLike Material on Its Mechanical Crushing Behaviors
Abstract
Indentation hardness and brittleness are the important factors to be considered in the study of rocklike materials’ mechanical crushing behaviors. The brittleness of rocklike materials is defined as the ratio of uniaxial compressive strength to tensile strength in this paper. In order to investigate the influences of hardness and brittleness on rocklike materials’ crushing behaviors, quartz sand and high strength αhemihydrated gypsum were used to prepare rocklike materials with different hardness and brittleness through different mass ratios. The artificial rocklike materials can eliminate the effects of natural rock’s weak structure plane on experimental results. The indentation test, Uniaxial compressive test, and Brazilian tensile test were conducted for characterizing the indentation hardness and brittleness of this artificial rocklike materials. The experimental results showed that brittleness increased with the increase of indentation hardness with high correlation coefficient. The confining stress presented a positive impact on the indentation hardness of the rocklike materials. Based on those mechanical properties, the numerical rock models were calculated to study rock crushing behaviors using discrete element method (DEM). A series of rock crushing tests were conducted to investigate the influences of hardness and brittleness of rocklike materials on rock crushing behaviors using a conical pick cutter. The numerical results showed that the normalized specific energy was negatively correlated with indentation hardness index (IHI). The normalized specific energy decreased with the increase of brittleness index (BI) with a high correlation coefficient. This study is beneficial in utilizing IHI and BI to evaluate the mechanical properties, failure patterns, and mechanical crushing efficiency of rocklike materials.
1. Introduction
Understanding of rock hardness and brittleness has a crucial importance for practical underground construction engineering. The correlation between hardness and brittleness of rock materials and their mechanical properties and failure patterns has been the focus of research hotspot in related fields. Szwedzicki [1] employed the indentation test for measuring rock hardness and stated that this test could be used for characterization of rocks’ mechanical properties. There exists a relationship between the value of the indentation hardness and the uniaxial compressive strength (UCS). Yagiz [2] stated that the indentation hardness could be directly used for measuring rock toughness and brittleness, and the measured index had a close correlation with UCS and density of rocks. Shalabi et al. [3] suggested that there is a good relationship between the engineering properties (such as strength, modulus of elasticity, and Poisson’s ratio) of the intact rock and its hardness (Schmidt, shore scleroscope, abrasion, and total hardness). Xuefeng et al. [4] investigated the influences of rock brittleness on rock crushing behaviors and suggested that rock brittleness had a dominant effect on rock failure patterns. Fan et al. [5] investigated the mechanical characterization of shale through instrumented indentation test and stated that the indentation hardness could be used to obtain the elastic modulus, hardness, and brittleness. Therefore, hardness and brittleness can be used as the important indexes to characterize the mechanical properties and failure patterns of rock materials.
In addition, many numerical studies had been conducted on the failure of brittle rock materials utilizing discrete element method. Li et al. [6] studied the brittle rock indentation process and rocktool interactions by particle flow code (PFC). The results revealed that the tensile cracks are resulted from the chipping process and the shear cracks are induced by the crushing force. Li et al. [7] investigated the crack initiation and propagation mechanism in brittle rocks using grainbased finitediscrete element method. The compression testing under different end friction, slenderness, loading rate, and confining stress are systematically performed to investigate the external influences on rock fracturing. It can be seen that the fragmentation patterns of brittle rocks are related to many mechanical properties of rocks and external factors. The DEM has advantages in the numerical investigation of brittle rocks’ mechanics fragmentation.
In practical underground constructions, mechanical rock crushing is a continuous, efficient and safe technique exhibiting great superiority in geotechnical engineering [8]. Many researchers [8–10] have investigated rock crushing behaviors through different indentation tests. Specific energy (SE) is always used to evaluate the rock crushing efficiency, which is defined as the amount of work required to break per unit volume of rocks. Yilmaz et al. [11] stated that the hybrid dynamic hardness had the potential to be used as an index rock material property that can be employed in preliminary estimations of specific energy (SE) and cutting force for a wide range of rocks under different operational conditions. Ji et al. [12] investigated the effects of ductilitybrittleness change on rock fragmentation process and its failure mode. As the brittleness of the model increased, the normalized specific energy decreased and the efficiency of the rock crushing increased. Consequently, this paper investigates the correlation between mechanical properties and failure patterns of rock materials and their brittleness as well as indentation hardness based on rock crushing behaviors.
Although using only one mechanical parameter always cannot accurately predict the rock crushing behaviors, a large number of studies have shown that there exists a close correlation between those rock mechanical parameters (such as strength, modulus of elasticity, toughness, Poisson’s ratio, and density) and hardness as well as brittleness. Balci et al. [13] tested 23 different rock specimens and investigated the predictability of SE from physical and mechanical properties. The results showed good correlations between SE and rock properties, such as UCS, Brazilian tensile strength (BTS), static and dynamic elastic moduli, and the Schmidt hammer value. Singh [14] observed that the in situ SE strongly depended upon the brittleness, UCS, Protodyakonov strength index, coefficient of internal friction, and the initial critical penetration stress. Therefore, it is worth investigating the influence of rock indentation hardness and brittleness on rock crushing behaviors. Because rock failure process is not visible, the mechanism of influence of hardness and brittleness on rock crushing behaviors remains unclear. Very few studies have investigated the hardness and brittleness of rock materials on their crushing behaviors utilizing pick cutters. Moreover, the IHI and BI indexes can be used to quantitatively characterize rocks’ mechanical crushing behaviors.
Indentation test is the most common method to obtain IHI of rock materials. In generally, the IHI is calculated by equation (1) [1, 2]. Here, is the maximum thrusting force of indenter at the specific penetration, is the penetration depth. Rock brittleness is defined as the ability of rocks to deform continuously without permanent deformations along with the application of stress surpassing the stress limit corresponding to rock microcracking [4, 15]. Liu [16] studied the relations between rock brittleness and rock fragmentation efficiency by disc cutter and summarized the existing 35 different brittleness indexes, which can be classified into seven categories with respect to strength, strain, strain energy, hardness, Mohr envelope, etc. Although the rock brittleness has not yet been standardized [16], most previous studies stated that BI could be calculated by the equations in Table 1 and stated that it exhibited a good correlation with rock crushing performance [16, 17].

As shown in Figure 1, natural rocks contain many joints and fissures [18], which are the important factors affecting rock mechanical properties. Gypsum has proven to be a good rocklike material for rock mechanics experiments and rock fragmentation behaviors [19, 20], while the artificial rocklike materials may be able to eliminate the testing errors caused by the inhomogeneity and discontinues of natural rocks. Therefore, quartz sand and highstrength αhemihydrated gypsum were used to prepared rocklike materials with different hardness and brittleness through different mass ratios. These rocklike materials can eliminate test errors caused by the inhomogeneity and discontinues of natural rocks. A series of indentation tests were conducted using a new indentation apparatus under confining stress conditions. In addition, the rock crushing models were calculated using a discrete element method code, i.e., Particle Flow Code in two dimensions (PFC2D). PFC2D is widely used in the simulation of rock crushing behaviors [4, 6, 21]. In this paper, the crushing models are calculated based on the UCS and BTS of the rocklike specimens. The cutting force and SE during fragmentation and failure patterns were investigated. The contribution of this study is that for the first time, the influence mechanism of rocklike materials’ hardness and brittleness on their crushing behaviors and failure patterns is investigated by artificial rocklike materials and DEM. The artificial rocklike materials can eliminate the effects of natural rock’s weak structure plane on study results. This investigation is of great significance for understanding the correlation between hardness and brittleness of rock materials and their fragmentation as well as failure patterns in rock crushing behaviors of underground engineering.
2. Indentation Test with Different Confining Stress Levels
2.1. Testing Apparatus and Specimen Preparation
In this paper, the indentation hardness testing is performed on a newly designed apparatus, which enables the indentation test of cylindrical specimen can be conducted under confining stress conditions. A schematic diagram of the newly designed apparatus is shown in Figure 2. The movement of indenter is controlled by a stepping motor. The penetration displacement is monitored by a linear variable differential transformer (LVDT), and the loading force is monitored by a 20 kN load sensor. A volumepressurecontroller (VPC) was used to load confining pressure. The conical indenter has a 90° cone and 0.25 mm radius spherical tip. The bearing platform is a hollow cylinder, and the indenter can pass through the bearing platform. A free rock surface is prepared for indentationtesting on the top of rock specimen.
The detail of the mass ratios of rocklike materials is shown in Table 2. The particle size of quartz sand is 0.18 mm (80 mesh sieve). As shown in Figure 3, the preparation process of specimens is as follows:(1)According to the designed mass ratio, the ingredients were weighed, mixed, and stirred at a high speed for 3 minutes.(2)The mixture was loaded into moulds and placed on the vibrostand to shake for 5 minutes.(3)Specimens were removed from moulds after 2 hours and maintained at room temperature for 28 days. All the end faces of specimens were corrected according to ISRM recommendations. The ratio of retarder to gypsum is equal to 0.1%; the ratio of defoamer to water is equal to 0.5%.

As shown in Figure 3, the dimensions of specimens for uniaxial compressive testing were designed as 50 mm and 100 mm in diameter and height. Moreover, the dimensions of specimens for Brazilian splitting testing were designed as 50 mm and 25 mm in diameter and height. In the indentation tests, the size of the specimens affects the indentation test results and, as such, guidelines as observed from the following studies may be used for eliminating the size effects. According to the cavity expansion model (CEM) theory [9, 22, 23], the ratio of the depth of the plastic zone induced by an axisymmetric indenter to the diameter of the rock specimen should be less than 0.167. Based on CEM theory, the ratio of the plastic zone depthtopenetration depth could be calculated with the equations (2)–(7). In the equations, is the critical radius of the plastic zone, n is a number characterizing the indenter geometry (for spherical and conical indenters, n = 2.0), is the passive coefficient, is the dilatancy coefficient (here, dilatancy angle is equal to friction angle , = 20°), is the Young’s modulus of rocks ( = 3.63 GPa), is Young’s modulus of indenter ( = 450 GPa). is Poisson’s ratio of rocks ( = 0.16), is Poisson’s ratio of indenter ( = 0.27), β is the angle between the indenter and the free surface of the rock (for this indenter, β = 45°), is the UCS of this rocklike material ( = 18.53 MPa), and α is the tip contact radius (α = 0.25 mm). By substituting the parameters to equations (2)–(7), the plastic zone radiustosample diameter ratio is equal to 0.045, which is less than 0.167. Therefore, the dimension of rock sample with 36 mm in diameter is reasonable.
2.2. Uniaxial Compressive Strength and Brazilian Tensile Strength Tests
In order to investigate the effects of hardness and brittleness on rocks fragmentation mechanism, the rocklike specimens with four different mass ratios were prepared in this study. A series of uniaxial compressive tests and Brazilian splitting tests were conducted to calculate the indexes of BI. The stressstrain curves of uniaxial compressive tests are shown in Figure 4, the displacementstress curves of Brazilian tensile tests are shown in Figure 5. It can be seen that the failure patterns all exhibited obvious plunge after peak stress, which suggested that these rocklike materials could characterize the brittle breaking characteristics of rocks. As shown in Table 3, the values range from 4.08 to 6.3, and the values range from 0.62 to 0.73. Figure 6 shows that the specimens in the uniaxial compression test mainly exhibited splitting failure induced by the propagation of vertical cracks. It also can be seen in the Brazilian tensile strength tests, as shown in Figure 7. Consequently, the rock numerical models used for investigating rock crushing behaviors are calculated based on the UCS, BTS, Young’s modulus, and the failure patterns.

2.3. Indentation Tests under the Unconfined Condition and Confining Stress Condition
A series of indentation tests were conducted under unconfined stress and 5 MPa confining stress conditions. The experimental results are shown in Figures 8 and 9. The velocity of the pick cutter is set to 1.0 mm/min according to Wang’s studies [8]. The confining stress restrained the propagation of vertical cracks. Since the slope of the curve is unstable at the initial stage of the indentation test, the slope of the fitted line near the maximum load is used as the indentation hardness, as shown in Figure 8. Moreover, the slope of the fitted line at the penetration displacement of 6 mm is used as the IHI under 5 MPa confining stress conditions. The IHIs under 5 MPa confining stress conditions are always less than that under unconfined stress conditions. According to the results in Table 3 and failure patterns in Figure 10, it can be seen that propagating crack length has negative correlation with IHI. As shown in Figures 10 and 11, there existed a compaction stage at the beginning of indentation test, and then the expansion of the crushed rocks induced the initiation and propagation of tension cracks. It can be seen in Figure 12 that the index [16] increased with the increase of IHI with a high correlation coefficient of 0.999. As shown in Figure 13, exhibited the same change law with IHI while and did not show a prominent correlation with IHI. The Figure 14 shows that the IHI has evident positive correlation with UCS. Therefore, the index is more suitable for the brittleness study in this paper.
3. Discrete Element Analysis of Crushing
3.1. Microparameters Calculating Based on UCS and BTS
In order to investigate the relation between hardness, brittleness, and rock crushing behaviors, the numerical models of rocklike materials were built through the flatjoint contact model (FJM) [24] in PFC2D [25]. As shown in Figure 15, the flatjoint contact simulates the behavior of an interface between two notional surfaces, each of which is connected rigidly to a piece of body. The particles of a flatjointed material are considered as faced grains, each of which is depicted as a core and a number of skirted faces. The interface is evenly discretized into several elements. According to previous researchers’ work, the numerical models of FJM is found to closely match real rocks’ macroproperties [26]. Therefore, the numerical models of rocklike materials were all calculated utilizing FJM in this study.
According to the mechanical parameters in Table 3 and the calibration procedure in Figure 16, the calculating results are shown in Figure 17. The microparameters can be calculated as shown in Table 4. The ratio of maximum to minimum diameter is set to 1.66 [27]. Table 5 shows the comparison of results obtained from numerical modeling and experimental tests. It suggested that the modeling results agree well with the results obtained from physical experiments. It can be seen that the failure patterns of numerical results exhibited the same characteristics with that of experimental results as shown in Figure 6. The different errors of all the parameters can be controlled within ±3.5%. In Figure 17, the blue cracks represent the shear cracks and the green cracks represent the tension cracks. Figures 17(a)–17(d) are the stressstrain curves and failure patterns of UCS tests, and Figure 17(e)–17(h) are the stressstrain curves and failure patterns of BTS tests.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)


3.2. Effects of Hardness and Brittleness on Rock Crushing Behaviors
In the numerical crushing tests, the numerical models of rocklike materials were established according to the experimental results (i.e., UCS, BTS, and Young’s modulus). More microscopic crushing characteristics can be observed in numerical analysis. The SE can be calculated by equation (8), where F_{i} is the recorded crushing force data and V_{chippings} is the total volume of crushing chippings. Since the SE is positively related to UCS, the normalized is proposed by equation (9). The failure patterns of crushing chippings of rocklike materials are shown in Figures 18 and 19.
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
Figure 20 shows that is negatively correlated with IHI. The normalized specific energy decreased exponentially with the increase of with a high correlation coefficient of 0.9217, as shown in Figure 21. Figure 18 shows that the main failure patterns of crushing behaviors keep in line with the failure patterns of rock tests in Figure 17. According to the numerical results of crushing force in Table 6, it can be seen that the mean of crushing forces under unconfined stress conditions is larger than that under 5 MPa confining stress conditions in the crushing tests of material A and C, while the standard deviation values of crushing forces are opposite. Comparing Figures 18 and 19, it can be seen that the propagation depth of cracks was restrained, but the index increased. Under 5 MPa confining stress conditions, is negatively correlated with IHI; the normalized specific energy decreased with the increase of BI with a high correlation coefficient of 0.98256, as shown in Figures 22 and 23.

The numerical results show that the average size of chippings of materials B and D are larger than that of materials A and C. Moreover, the propagation length of tensile cracks is larger than that of shear cracks. As shown in Tables 7 and 8, the number of shear cracks of materials A and C is greater than that of materials B and D, but the number of tension cracks of materials A and C is less than that of materials B and D. Therefore, the failure patterns of rocklike materials with higher brittleness are more prone to tensile failure in the crushing testing.


In equation (10), is the equivalent diameter of chipping, is the number of particles of chipping, and is the mean diameter of particle. The pore ratio of the numerical model is 0.12. Comparing A and C, it is found that the larger the indentation hardness, the larger the proportion of the particles (equivalent diameter at 0–0.008 m range) in the cumulative gradation curve, as shown in Figures 24 and 25. The same is true between B and D.
4. Conclusions
In this paper, quartz sand and highstrength αhemihydrated gypsum were prepared to make rocklike materials with different hardness and brittleness through different mass ratios. Gypsum has proven to be a good rocklike material for rock mechanics experiments and rock fragmentation behaviors. The numerical models were calculated based on the mechanical parameters of those rocklike materials. The influence mechanism of indentation hardness and brittleness of rocklike materials on their crushing behaviors was investigated using the DEM, i.e., PFC2D. The main conclusions can be drawn as follows.
The index increased exponentially with the increase of IHI with highcorrelation coefficient. The index IHI has evident positive correlation with UCS. Hardness and brittleness can be used as important indexes to characterize the mechanical properties and failure patterns of rock materials. According to the numerical results of the crushing tests, the normalized specific energy is negatively correlated with IHI. decreased exponentially with the increase of . The failure patterns of rocklike materials with higher brittleness are more prone to tensile failure in the crushing testing.
This investigation is of great significance for understanding the correlation between hardness and brittleness of rock materials and their fragmentation as well as failure patterns in underground constructions. The FJM in PFC is more suitable for the numerical study of rocklike materials than finite element method, and it can obtain more coincident mechanical parameters and more microscopic failure characteristics.
Data Availability
All data used to support the findings of this study are available from the corresponding author upon request (list items: design drawings of the testing appartuas; the experimental test data; all the PFC codes).
Conflicts of Interest
The authors declare no conflicts of interest.
Acknowledgments
The research was supported by the National Key Research and Development Program of China (Grant no. 2016YFC0600904).
References
 T. Szwedzicki, “Indentation hardness testing of rock,” International Journal of Rock Mechanics and Mining Sciences, vol. 35, no. 6, pp. 825–829, 1998. View at: Publisher Site  Google Scholar
 S. Yagiz, “Assessment of brittleness using rock strength and density with punch penetration test,” Tunnelling and Underground Space Technology, vol. 24, no. 1, pp. 66–74, 2009. View at: Publisher Site  Google Scholar
 F. I. Shalabi, E. J. Cording, and O. H. AlHattamleh, “Estimation of rock engineering properties using hardness tests,” Engineering Geology, vol. 90, no. 34, pp. 138–147, 2007. View at: Publisher Site  Google Scholar
 L. Xuefeng, W. Shibo, G. Shirong, R. Malekian, and L. Zhixiong, “Investigation on the influence mechanism of rock brittleness on rock fragmentation and cutting performance by discrete element method,” Measurement, vol. 113, pp. 120–130, 2018. View at: Publisher Site  Google Scholar
 M. Fan, Y. Jin, M. Chen, and Z. Geng, “Mechanical characterization of shale through instrumented indentation test,” Journal of Petroleum Science and Engineering, vol. 174, pp. 607–616, 2019. View at: Publisher Site  Google Scholar
 X. F. Li, H. B. Li, Y. Q. Liu, Q. C. Zhou, and X. Xia, “Numerical simulation of rock fragmentation mechanisms subject to wedge penetration for TBMs,” Tunnelling and Underground Space Technology, vol. 53, pp. 96–108, 2016. View at: Publisher Site  Google Scholar
 X. F. Li, H. B. Li, L. W. Liu, Y. Q. Liu, M. H. Ju, and J. Zhao, “Investigating the crack initiation and propagation mechanism in brittle rocks using grainbased finitediscrete element method,” International Journal of Rock Mechanics and Mining Sciences, vol. 127, Article ID 104219, 2020. View at: Publisher Site  Google Scholar
 S. Wang, X. Li, K. Du, and S. Wang, “Experimental investigation of hard rock fragmentation using a conical pick on true triaxial test apparatus,” Tunnelling and Underground Space Technology, vol. 79, no. 1, pp. 210–223, 2018. View at: Publisher Site  Google Scholar
 L. H. Chen and J. F. Labuz, “Indentation of rock by wedgeshaped tools,” International Journal of Rock Mechanics and Mining Sciences, vol. 43, no. 7, pp. 1023–1033, 2006. View at: Publisher Site  Google Scholar
 L. J. Yin, Q. M. Gong, H. S. Ma, J. Zhao, and X. B. Zhao, “Use of indentation tests to study the influence of confining stress on rock fragmentation by a TBM cutter,” International Journal of Rock Mechanics and Mining Sciences, vol. 72, pp. 261–276, 2014. View at: Publisher Site  Google Scholar
 N. G. Yilmaz, D. Tumac, and R. M. Goktan, “Rock cuttability assessment using the concept of hybrid dynamic hardness (HDH),” Bulletin of Engineering Geology and the Environment, vol. 74, no. 4, pp. 1363–1374, 2015. View at: Publisher Site  Google Scholar
 P. Q. Ji, X. P. Zhang, and Q. Zhang, “Partide flow code analysis of effect of ductilitybrittleness change on TBM cutters rock fragmentation process and its failure mode,” Yantu Lixue/Rock and Soil Mechanics, vol. 37, no. S2, pp. 724–734, 2016. View at: Google Scholar
 C. Balci, M. A. Demircin, H. Copur, and H. Tuncdemir, “Estimation of optimum specific energy based on rock properties for assessment of roadheader performance,” SIAMM—Journal of the Southern African Institute of Mining and Metallurgy, vol. 104, no. 11, pp. 633–642, 2004. View at: Google Scholar
 S. P. Singh, “Criterion for the assessment of the cuttability of coal,” in Advances In Mining Science And Technology, pp. 225–239, Elsevier, Amsterdam, Netherlands, 1987. View at: Publisher Site  Google Scholar
 G. E. Andreev, Brittle Failure of Rock Materials, CRC Press, Boca Raton, FL, USA, 1995.
 Q. Liu, “Evaluation of rock brittleness indexes on rock fragmentation efficiency by disc cutter,” Yanshilixue Yu Gongcheng Xuebao/Chinese Journal of Rock Mechanics and Engineering, vol. 35, no. 3, pp. 498–510, 2016. View at: Google Scholar
 R. Altindag, “Correlation of specific energy with rock brittleness concepts on rock cutting,” Journal of the Southern African Institute of Mining and Metallurgy, vol. 103, no. 3, pp. 163–171, 2003. View at: Google Scholar
 X.X. Yang and W.G. Qiao, “Numerical investigation of the shear behavior of granite materials containing discontinuous joints by utilizing the flatjoint model,” Computers and Geotechnics, vol. 104, pp. 69–80, 2018. View at: Publisher Site  Google Scholar
 M. Chen, S. Yang, R. P. Gamage et al., “Fracture processes of rocklike specimens containing nonpersistent fissures under uniaxial compression,” Energies, vol. 12, no. 1, p. 79, 2019. View at: Publisher Site  Google Scholar
 L. N. Y. Wong and H. H. Einstein, “Crack coalescence in molded gypsum and carrara marble: part 1. Macroscopic observations and interpretation,” Rock Mechanics and Rock Engineering, vol. 42, no. 3, pp. 475–511, 2009. View at: Publisher Site  Google Scholar
 X. Li, S. Wang, R. Malekian, S. Hao, and Z. Li, “Numerical simulation of rock breakage modes under confining pressures in deep mining: an experimental investigation,” IEEE Access, vol. 4, pp. 5710–5720, 2016. View at: Publisher Site  Google Scholar
 L. J. Yin, Q. M. Gong, and J. Zhao, “Study on rock mass boreability by TBM penetration test under different in situ stress conditions,” Tunnelling and Underground Space Technology, vol. 43, pp. 413–425, 2014. View at: Publisher Site  Google Scholar
 H. Alehossein, E. Detournay, and H. Huang, “An analytical model for the indentation of rocks by blunt tools,” Rock Mechanics and Rock Engineering, vol. 33, no. 4, pp. 267–284, 2000. View at: Publisher Site  Google Scholar
 D. O. Potyondy, “PFC2D FlatJoint Contact Model,” in Technical Memorandum ICG7138L, Itasca Consulting Group, Inc., Minneapolis, MN, USA, 2012. View at: Google Scholar
 Itasca Consulting Group Inc, PFC Manual for Version 5.0, Itasca Consulting Group Inc., Minneapolis, MN, USA, 2018.
 S. Wu and X. Xu, “A study of three intrinsic problems of the classic discrete element method using flatjoint model,” Rock Mechanics and Rock Engineering, vol. 49, no. 5, pp. 1813–1830, 2016. View at: Publisher Site  Google Scholar
 X. Ding and L. Zhang, “A new contact model to improve the simulated ratio of unconfined compressive strength to tensile strength in bonded particle models,” International Journal of Rock Mechanics and Mining Sciences, vol. 69, pp. 111–119, 2014. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Jiuqun Zou et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.