Research Article  Open Access
Discrete Element Simulation of Conical Pick’s Coal Cutting Process under Different Cutting Parameters
Abstract
To explore how the cutting parameters affect the conical pick’s cutting force in coal cutting process, some simulation tests were carried with EDEM based on analysis of coal cutting process of conical pick and micro mechanical characteristic of coal particles. In this paper, traction speed, drum angular velocity, and installation angle of conical pick were considered as the main coal cutting parameters. The discrete element simulations were conducted under 7 sets of various cutting parameters, which were determined in the range of empirical parameters. From the results of the simulation tests, it can be concluded that the fragmentation of coal particles appears in a compacted, squeezed, and crumbling state in coal cutting process; the traction speed has the greatest impact on the conical pick’s cutting force and cutting force fluctuation and the influence of drum angular velocity and installation angle reduce in turn; the force and force fluctuation tend to decrease and obviously increase, respectively, with the increase of drum angular velocity and traction speed; and they will reach the minimum when the pick’s installation angle is 45° under certain traction speed and drum angular velocity.
1. Introduction
Various coal cutting works are carried out in mining engineering by means of different machines and cutting tools [1]. The shearer and conical pick, as the core device and main cutting tool, are widely used in mining engineering, and the basic physical phenomenon occurring during coal cutting process is fracturing and fragmentation of the coal under mechanical action of conical picks, of which the cutting parameters not only affect the cutting efficiency of shearer, but also influence the cutting power and abrasion of the conical picks. Therefore, the research on the cutting parameters of conical pick in coal cutting process has become a hotspot in this field.
From the cutting force model of conical pick proposed by Evans on the basis of tensile stress theory [2], the cutting force was considered to be related to tensile stress of coal, semitip angle of conical pick, cutting angle, and cutting thickness and was proportional to the square of the cutting thickness. According to analyzing the friction between conical picks and coal, the theoretical model of the conical picks has been obtained based on Evans theory by Roxborough and Liu [3]. When the cutting friction angle and front angle of pick were introduced into Evans’ model, the value of cutting force was modified under semitip angle at 0 degree, and then a semiempirical cutting model was proposed by Göktan [4] and Göktan and Gunes [5] with regression method on the basis of a large number of experiments. A theoretical model to estimate the cutting force of conical pick based on fracture mechanics was presented by Bao et al. [6], in which the cutting force work was considered to be proportional to the square of the depth of cut, and the cutting peak force and the depth of cut have a power law relationship with the component equal to 4/3. A mathematical model on the lateral force of conical pick under different cutting conditions has been developed by Liu and Li [7] based on cutting mechanism and mechanics principle; then spectrum of the lateral force load at different slope angles and cutting thicknesses was decomposed in frequency domain and time domain.
With the development of computer technology and modern engineering technology, a number of numerical tests of rock cutting utilized the finite element method have been explored [8–12]. However, the finite element method based on the continuum mechanics theory of material modeling has some drawbacks in representing properly discontinuities of the material occurring during rock cutting [13]. The discrete element method was pioneered by Cundall [14], which has a wide range of application in rock mechanics [15, 16] and soil mechanics [17, 18], as well as various penetration conditions in granular media [19, 20].
The discrete element method that treats coal or rock as being composed of separate rigid particles has extensive prevalence in coal cutting process. Particle Flow Code (PFC) was used to study the cutting force of the picks in linear [21–23] and rotary [1, 24] cutting process when the cutting thickness was constant and to develop the crack propagation [22, 23, 25] of coal and rock in the linear cutting process. Extended Discrete Element Method (EDEM) was used to establish linear conical pick cutting coal model to study the cutting process with different cutting thickness and different installation angles [26], which showed that the cutting force of conical pick increases linearly with the cutting thickness and varied nonlinearly with the installation angle.
In the coal mining operation, conical picks, installed on the shearer drum at different installation angles, are driven by traction speed of the shearer and the angular velocity of the drum to cut the coal. However, the cutting condition, like the prolate cycloid cutting trajectory and the varying cutting thickness, is not distinctly considered in the above scholars’ research [21–26] based on simplified assumptions. Based on the coal mining operation, this paper simulated the coal cutting process of conical pick at various traction speed, drum angular velocity, and installation angle in EDEM, to analyze the phenomenon of coal fragmentation and the cutting force of conical pick, as well as the relation between cutting force and the cutting parameters mentioned above.
2. Cutting Parameters of Conical Pick
2.1. Determination of Cutting Parameters
Figure 1 presents the coal cutting process of conical pick. The pick is mounted on shearer drum with the installation angle and it is driven by the shearer traction speed and the drum angular velocity to cut the coal, where is cutting radius, the points and are coordinate origins located at the drum center axis at the beginning and ending of the coal cutting process, trajectories 1 and 2 are formed by the first and the second cutting cycle, respectively, in coal cutting process, between which is the cutting thickness, the formula of cutting trajectory is given in (1), and is the maximum cutting thickness in cutting process.where ,.
In coal cutting process, the cutting thickness varies with the traction speed of shearer, drum angular velocity, and cutting radius. If (1) represents trajectory 1, the formula of trajectory 2 can be deduced as follows:
where is the number of picks on the same transversal. Then the change of cutting thickness in coal cutting process can be obtained as
The detailed formula of cutting thickness can be calculated according to (1), (2), and (3).
Under certain working condition of shearer, this paper selected the semitip angle of conical pick equal to the optimum value of 38 degrees [2], and the drum diameter was set as 1 m (cutting radius is 0.5 m) for the tests. Then the traction speed of shearer , angular velocity of drum , and installation angle were selected as the parameters to be studied in the coal cutting process.
2.2. The Range of Cutting Parameters
At present, the above three parameters have formed the empirical range, respectively, [27]. In this paper, 7 sets of various parameters were selected to simulate the cutting process on the basis of empirical ranges, which were classified as different angular velocity, traction speed, and installation angle, as shown in Table 1.
 
Notes. The data in the Table 1 are the cutting parameters range of single pick. 
3. Modeling of Coal Cutting Process in EDEM
3.1. Materials and Contact Model of Coal and Conical Pick
Conical pick and coal materials were selected, respectively, as steel and coal from EDEM material database, and the contact model between coal and conical pick was chosen as HertzMindlin (with no slip) contact model in this paper. The parameters of contact property [28] are given in Table 2.

3.2. Building of Coal Model
3.2.1. Parameters of Coal Particle
The structural and mechanical properties of coal are complex, which present certain regularity bedding, joints, and random waste, pores, cracks, and so on. In order to simplify the modeling and simulation experiments, the circular particles with equal radius were simulated as coal particle in EDEM. The parameters of coal particles are shown in Table 3.

3.2.2. Mechanical Parameters of Coal Particles
Since the coal is massive structure state before cutting and appears crumbling during the cutting process, certain bonding force must be existent between the coal particles. The HertzMindlin with bonding was selected as the contact model between coal particles in this paper. The bonding model was depicted in Figure 2 [29–31].
Where , are the radius of particles 1 and 2, respectively, , are the contact radius of particles 1 and 2, is the normal overlap between particles, , are normal and tangential damping, and , are normal contact stiffness and tangential contact stiffness between particles. Their relations [32] are obtained:
where is equivalent elastic modulus, , and are macroscopic elastic modulus of coal particles 1 and 2, respectively, and are macroscopic Poisson ratio of coal particles 1 and 2, is equivalent elastic modulus, , and , are macroscopic shear modulus of coal particles 1 and 2. is equivalent radius of coal particles, . is equivalent mass of coal particles, . is the restitution coefficient of coal particles in Table 2. And the macro mechanical characteristic parameters of coal particles are shown in Table 4 [33].

The normal contact stiffness and the tangential contact stiffness were calculated by the (5) and (6) and the relevant parameters from Tables 3 and 4. Then the Micro mechanical parameters of coal particles were set up in EDEM, which were shown in Table 5, where the critical normal stress and the critical shear stress are equal to the compressive strength and shear strength of the coal, respectively.

3.3. Building of the Geometric Model
In this paper, a simplified threedimensional drum geometric model with a single conical pick was established using UG, and the model was imported into the EDEM geometry module to simulate the coal cutting process of conical pick.
In order to simulate the coal cutting process showed in Figure 1, every simulation test was carried under different coal cutting surface, which was determined by the cutting trajectory (1) based on the cutting parameters shown in Table 1. In EDEM, the length, width, and height of coal model were determined as 0.68 m, 0.15 m, and 1.30 m according to the drum diameter, physical dimension of conical pick, and the related cutting parameters, and its geometry model is consisted by 90000 coal particles, which were dynamically deposited by Rain Drop Method. When the micro mechanical parameters of coal particles were set according to Table 5, the bonding between coal particles was generated. The model of coal cutting process of conical pick in EDEM was given in Figure 3.
3.4. Time Step of Simulation
If the time step is not appropriate, it may cause particle divergence and affect the accuracy of simulation results. The Rayleigh time should be set reasonably for the simulation tests to go smoothly in EDEM. According to the literature [32], the Rayleigh time step is the maximum of the time step of particle set, determined by
where is the Rayleigh time step, is coal particle radius, is particle density, is shear modulus, and is Poisson ratio.
Rayleigh time step was calculated as according to the related parameters in Tables 3 and 4. In the actual simulation test, the time step is usually set from the range of 10–40% . In this paper the time step was 20% .
4. Discussion of Simulation Results
4.1. Analysis of Coal Cutting Process
The simulation of coal cutting process of conical pick under the parameters ① from Table 1 is shown in Figure 4, where (a)~(i) are cutting state, particle velocity vector, and sectional of bonding diagrams when conical pick is at the beginning, middle, and end position of the cutting process. In the particle velocity vector diagrams, the direction of the arrows indicates the instantaneous velocity direction of the coal particles, and the color of the arrows represents the instantaneous velocity of particles, which increases gradually with blue, green, and red in turn. And in the section of bonding diagrams, the strip represents the bonding between the coal particles, which constitutes bond networks; the red and blue color represent large and little cohesive force, respectively.
(a) Cutting state at the beginning position
(b) Particle velocity vector diagram at the beginning position
(c) Section of bonding diagram at the beginning position
(d) Cutting state at the middle position
(e) Particle velocity vector diagram at the middle position
(f) Section of bonding diagram at the middle position
(g) Cutting state at the end position
(h) Particle velocity vector diagram at the end position
(i) Section of bonding diagram at the end position
As indicated in Figure 4, the cutting thickness is small at the beginning of coal cutting process, and the pick causes small clumps of coal particles, when the internal velocity pattern of coal tends to be more stable than that at the tip of the conical pick, and the falling massive coal particles appear irregular motion state; meanwhile, the bonding between coal particles has few fractures that belong to compaction state of coal. With the continuous cutting of the conical pick, the cutting thickness increases, and the cutting thickness reaches the maximum at the middle position of cutting process, and near this position, coal particles are gradually piled up at the joint between the conical pick and pick holder, where the bonding is extruded and deformed, and then a large degree of fracture occurs, showed by blue colored bonding, which results in coal gliding fracture. And the velocity of the coal particles is high at the tip and the end of the conical pick, which causes coals to collapse along the movement direction of the conical pick. When the pick is close to the end of the coal cutting process, cutting thickness decreases gradually, and velocity of coal particles is still high at the tip of conical pick. The falling coal moves with cutting direction of conical pick. The particle velocity vector and bonding diagram proposed in this paper shows fracture mechanism and movement direction of coal particles, which strengthen and prove the rationality of results on coal cutting state [34, 35].
4.2. Analysis of Cutting Force under Different Cutting Parameters
Figures 5, 6, 7, 8, 9, 10, and 11 are cutting force curve of conical pick in the cutting coal process with the parameters ①, ②, ③, ④, ⑤, ⑥, and ⑦, respectively. It can be seen from the diagram that the traction speed is the primary factor affecting the cutting force of the conical pick, and the drum angular velocity is the secondary and the installation angle is the least. In coal cutting process, the cutting force of the conical pick, which tends to increase first and then decrease with the change of cutting thickness, shows irregular fluctuations on the whole because of randomness of the lump coal fragmentation. This trend is consistent with the experimental results [24] and theoretical calculation results [27] in the literature.
The general trend of cutting force is similar to the cutting thickness obtained from (4) and the larger the cutting thickness, the greater the cutting force. At the initial stage of coal cutting process, the coal was effected by cutting force of conical pick, and the coal fragmentation occurs when the cutting force is accumulated to a certain extent, then the cutting force is obviously reduced until the end of cutting process.
4.2.1. Force Analysis at Different Angular Velocity
According to the comparison of the cutting force curves in Figures 5, 6, and 7, the cutting force of conical pick decreases with the increase of the angular velocity , and the smaller the angular velocity of the drum, the greater the variation of the cutting force. When the drum angular velocity is 5.23 rad/s, the maximum fluctuation of the cutting force is about 8000 N. To analyze this phenomenon, the coal particles are easily piled up when the angular velocity of the drum is small, and great pressure is formed before caving; the coal fragmentation is obviously larger than the cutting thickness, which results in sudden increase in cutting force. There will be a temporary noload phenomenon that occurs after the coal was cut.
4.2.2. Force Analysis at Different Traction Speed
Comparing the cutting force curves in Figures 8, 6, and 9, it can be seen that the cutting force of conical pick increases obviously with the increase of the traction speed. Also, the fluctuation of cutting force increases with the increase of traction speed; the maximum fluctuation of the cutting force is about 9000 N when the traction speed is 0.09 m/s.
4.2.3. Force Analysis at Different Installation Angle
The installation angle has little effect on the cutting force, and the average cutting forces of the conical pick are about 6500 N, 6000 N, and 6300 N when the installation angle is set as 40°, 45°, and 50°, respectively, by comparing the cutting force curves in Figures 10, 6, and 11. The interference of extrusion between coal and conical pick and its holder is not prone to happen because the conical pick has an advantage position to cut the coal when the installation angle is 45°. Therefore, the value and fluctuation of the cutting force are smaller when the installation angle is 45°.
5. Conclusions
The simulation tests are further studied based on theoretical study [27] about conical pick cutting coal, which considered the influence of changed cutting thickness in actual working condition. The conclusions in this paper are as follows:
On the basis of analyzing the coal cutting process of conical pick and the micro mechanics of coal particles, EDEM was used to establish coal and conical pick model to simulate the coal cutting process at different parameters. The fragmentation direction and velocity field of coal particle and massive coal, as well as the bonding fracture between coal particles, were revealed in the coal cutting process; the fragmentation of coal particles appears in a compacted, squeezed, and crumbling state in general.
The cutting force and fluctuation of cutting force are greatly impacted by traction speed, affected by drum angular velocity and installation angle decrease in turn. They tend to decrease and obviously increase, respectively, with the increase of drum angular velocity and traction speed. And the cutting force and its fluctuation are least affected by the installation angle.
Both the cutting force and its fluctuation will reach their minimum values when the installation angle is 45° under certain traction speed and drum angular velocity.
The structural characteristics of coal, such as fracture, joint, and pore, are not taken into account because of the software limitations and the particularities of the generated particles method (Rain Drop Method), as well as the vibration, abrasion of conical pick, and the smoothness of cutting surface formed by cutting trajectory. These factors will impact the cutting force of conical pick, but the impact will be far less than the impact of the parameters discussed in this paper.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This research was supported by the National Nature Science Fund (Grant no. 51674155) and High Educational Institution of Shandong Science and Technology Project (Grant no. YB06).
References
 J. Rojek, E. Oñate, C. Labra, and H. Kargl, “Discrete element simulation of rock cutting,” International Journal of Rock Mechanics & Mining Sciences, vol. 48, no. 6, pp. 996–1010, 2011. View at: Publisher Site  Google Scholar
 I. Evans, “A theory of the cutting force for pointattack picks,” International Journal of Mining Engineering, vol. 2, no. 1, pp. 63–71, 1984. View at: Publisher Site  Google Scholar
 F. F. Roxborough and Z. C. Liu, “Theoretical considerations on pick shape in rock and coal cutting in,” in Proceedings of the sixth underground operators conference, pp. 189–193, 1995. View at: Google Scholar
 R. M. Göktan, “A suggested improvement on Evans cutting theory for conical bits,” in Fourth International Symposium on Mine Mechanization and Automation, pp. 57–61, Queensland, Australia, 1997. View at: Google Scholar
 R. M. Göktan and N. Gunes, “A semiempirical approach to cutting force prediction for pointattack picks,” The Journal of The South African Institute of Mining and Metallurgy, vol. 105, no. 4, pp. 257–263, 2005. View at: Google Scholar
 R. H. Bao, L. C. Zhang, Q. Y. Yao, and J. Lunn, “Estimating the peak indentation force of the edge chipping of rocks using single pointattack pick,” Rock Mechanics and Rock Engineering, vol. 44, no. 3, pp. 339–347, 2011. View at: Publisher Site  Google Scholar
 C. Liu and D. Li, “Experimental research and theoretical model on lateral force of conical pick under different cutting conditions,” Journal of China Coal Society (Chinese), vol. 41, pp. 2359–2366, 2016. View at: Google Scholar
 P. Li, H. Zhang, S. Jiang, and W. Zhang, “Analysis and Testing of Load Characteristics for RotaryPercussive Drilling of Lunar Rock Simulant with a Lunar Regolith Coring Bit,” Shock and Vibration, vol. 2017, Article ID 3012749, 2017. View at: Publisher Site  Google Scholar
 S. Q. Kou, P.A. Lindqvist, C. A. Tang, and X. H. Xu, “Numerical simulation of the cutting of inhomogeneous rocks,” International Journal of Rock Mechanics and Mining Sciences, vol. 36, no. 5, pp. 711–717, 1999. View at: Publisher Site  Google Scholar
 J. Jonak and J. Podgórski, “Mathematical model and results of rock cutting modeling,” Journal of Mining Science, vol. 37, no. 6, pp. 615–618, 2001. View at: Publisher Site  Google Scholar
 B. Yu and A. W. Khair, “Numerical modeling of rock ridge breakage in rotary cutting,” in ARMA general meeting, British Columbia, Canada, 2007. View at: Google Scholar
 L. Monforte, M. Arroyo, J. M. Carbonell, and A. Gens, “Numerical simulation of undrained insertion problems in geotechnical engineering with the Particle Finite Element Method (PFEM),” Computers & Geosciences, vol. 82, pp. 144–156, 2017. View at: Publisher Site  Google Scholar
 Z. Z. Liang, H. Xing, S. Y. Wang, D. J. Williams, and C. A. Tang, “A threedimensional numerical investigation of the fracture of rock specimens containing a preexisting surface flaw,” Computers & Geosciences, vol. 45, pp. 19–33, 2012. View at: Publisher Site  Google Scholar
 P. A. Cundall, “A computer model for simulating progressive large scale movements in blocky rock system in blocky rock system,” in third international rock mechanics symposium, vol. 1, p. 2, 1971. View at: Google Scholar
 D. O. Potyondy and P. A. Cundall, “A bondedparticle model for rock,” International Journal of Rock Mechanics and Mining Sciences, vol. 41, no. 8, pp. 1329–1364, 2004. View at: Publisher Site  Google Scholar
 X. Liu, S. Liu, and P. Tang, “Coal fragment size model in cutting process,” Powder Technology, vol. 272, pp. 282–289, 2015. View at: Publisher Site  Google Scholar
 P. A. Cundall and O. D. L. Strack, “The distinct numerical model for granular assemblies,” Géotechnique, vol. 29, no. 1, pp. 47–65, 1979. View at: Publisher Site  Google Scholar
 M. B. Cil and K. A. Alshibli, “3D assessment of fracture of sand particles using discrete element method,” Geotechnique Letters, vol. 2, no. 79, pp. 161–166, 2012. View at: Publisher Site  Google Scholar
 J. K. Holmen, L. Olovsson, and T. Børvik, “Discrete modeling of lowvelocity penetration in sand,” Computers & Geosciences, vol. 86, pp. 21–32, 2017. View at: Publisher Site  Google Scholar
 M. O. Ciantia, M. Arroyo, J. Butlanska, and A. Gens, “DEM modelling of cone penetration tests in a doubleporosity crushable granular material,” Computers & Geosciences, vol. 73, pp. 109–127, 2016. View at: Publisher Site  Google Scholar
 O. Su and N. A. Akcin, “Numerical simulation of rock cutting using the discrete element method,” International Journal of Rock Mechanics and Mining Sciences, vol. 48, no. 3, pp. 434–442, 2011. 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. Yanxin, L. Haibo, Z. Xiaohua, and T. Liping, “Bondedcluster simulation of rockcutting using PFC2D,” Cluster Computing, vol. 20, no. 2, pp. 1289–1301, 2017. View at: Publisher Site  Google Scholar
 Y. Dai, L. S. Chen, X. Zhu, and H. Liu, “Modelling and simulation of a mining machine excavating seabed massive sulfide deposits,” International Journal of Simulation Modelling, vol. 15, no. 2, pp. 377–387, 2016. View at: Publisher Site  Google Scholar
 S.O. Choi and S.J. Lee, “Numerical study to estimate the cutting power on a disc cutter in jointed rock mass,” KSCE Journal of Civil Engineering, 2015. View at: Publisher Site  Google Scholar
 X. Chen, S. A. Miedema, and C. Van Rhee, “Numerical modeling of excavation process in dredging engineering,” in Proceedings of the 7th World Congress on Particle Technology, WCPT 2014, pp. 804–814, China, May 2014. View at: Publisher Site  Google Scholar
 J.X. Liu, C. Ma, Q.L. Zeng, R. Gao, and Z.Y. Jiao, “Research on cutting force of conical pick in cutting coal process,” Meitan Xuebao/Journal of the China Coal Society, vol. 42, no. 5, pp. 1325–1330, 2017. View at: Publisher Site  Google Scholar
 Z. Xiaolong, W. Yiliang, and Y. Zhaojian, “Simulation and analysis on coal cutting process based on EDEM,” Coal Technology (Chinese), vol. 34, no. 8, pp. 243244, 2015. View at: Google Scholar
 A. B. Yu, N. Standish, and L. Lu, “Coal agglomeration and its effect on bulk density,” Powder Technology, vol. 82, no. 2, pp. 177–189, 1995. View at: Publisher Site  Google Scholar
 J.W. Zhou, Y. Liu, C.L. Du, and S.Y. Liu, “Effect of the particle shape and swirling intensity on the breakage of lump coal particle in pneumatic conveying,” Powder Technology, vol. 317, pp. 438–448, 2017. View at: Publisher Site  Google Scholar
 K. Iwashita and M. Oda, “Microdeformation mechanism of shear banding process based on modified distinct element method,” Powder Technology, vol. 109, no. 13, pp. 192–205, 2000. View at: Publisher Site  Google Scholar
 J. Quist and C. M. Evertsson, Cone Crusher Modeling and Simulation, Department of Product and Production Development, Chalmers University of Technology, 2012. View at: Publisher Site
 W. Shen and B. Zhang, “Testing Study on Mechanical Parameters of coal,” Chinese Journal of Rock Mechanics and Engineering (Chinese), vol. 19, pp. 860–862, 2000. View at: Publisher Site  Google Scholar
 P. L. Menezes, M. R. Lovell, I. V. Avdeev, and C. F. Higgs, “Studies on the formation of discontinuous rock fragments during cutting operation,” International Journal of Rock Mechanics and Mining Sciences, vol. 71, pp. 131–142, 2014. View at: Publisher Site  Google Scholar
 L. Scholtes and de ricVictor Donze. F, “Modelling progressive failure in fractured rock masses using a 3D discrete element method,” nternational Journal of Rock Mechanics & amp; Mining Sciences, vol. 52, pp. 18–30, 2012. View at: Google Scholar
Copyright
Copyright © 2018 Jinxia Liu 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.