Numerical Simulation on Cutting and Fracturing of Rock Plate with One Side Fixed and Three Sides Free

Conical pick is a rock cutting tool that is commonly used in roadway driving. Pick wear frequently happens in the course of breaking hard rock. The current paper shows a new method to solve the problem of pick wear. The rock is preslit with the saw blade and then broken by the conical pick. In order to study the cutting force and features of rock fragment, the numerical model is built between rock plate and conical pick. And element erosion is added in the code to obtain the fracture result. The rock plate cutting testbed is made to testify the correctness of numerical simulation. The width, height, and thickness of the rock plate, as well as cutting angle and cutting position, which influence cutting force and rock fracture are studied. According to the results, there exist exponential relationships between cutting force and width and thickness of rock plate. In addition, a linear relationship is found between the cutting force and the height of rock plate. Furthermore, both the cutting angle and cutting depth have an influence on cutting force. In particular, the factors of thickness and height or rock plate have the most obvious influence on cutting force. It is proven that what is beneficial to rock fracture is higher height and lower thickness of rock plate.


Introduction
Mechanical excavation is the most widely used rock excavation approach in the mining and civil industries [1]. As a common cutting tool, the selection of conical pick is influenced by many factors, such as the strength of rock, rock cutting depth, and the cost of the cutting tool itself. erefore, the research of the rock fracture mechanism is particularly important for conical pick design.
In order to improve the efficiency of rock breaking, scholars have done comprehensive research on the rock cutting mechanism by various methods. e most classical theoretical study shows the cutting depth of conical pick and compressive and tensile strength of rock as the main parameter influencing cutting force [2]. However, when the tip angle of the conical pick is equal to zero, the cutting force still exists. Friction angle is introduced to overcome the shortcoming of the theoretical model; meanwhile, the experiments indicate that the results obtained from the new theories are closer to reality [3][4][5][6]. Liu and Jin [7] prove that cutting force is also in relevance to the number and layout of the conical pick. e influential factors on cutting force are systematically researched in experimental approach. e factors influencing cutting force are more than uniaxial compressive strength and cutting depth; meanwhile, elastic modulus and tensile strength are also the significantly influencing factors [8]. To reduce the experimental cost, Kang et al. [9] build a small capacity cutting testbed, and the feasibility of the scheme is proved by experiment and simulation. Dewangan and Chattopadhyaya [10] give a cutting angle optimization suggestion based on cutting efficiency and temperature of the conical pick. With the rock cutting machine, the cutting performance at the relieved cutting method is predicted [11]. Huang et al. [12] take deep mining condition into account and research the relationship between confining stress and cutting force. Wang et al. [13] propose a specific energy model with regression methods, which is verified by experiment. Su [14] employs PFC 3D code to build the interaction model between cutter and rock and gives the relationship between rake angle and cutting depth. Meanwhile, he demonstrates that the Nishimatsu theory is only suitable for hard rock. Wang and Su [15] point out that the maximum volume of rock fragment is proportional to cutting depth and that the specific energy decreases exponentially with the increase of cutting depth. Yasar [16,17] points out that there are inherent defects in the traditional theoretical model. Test results of different rocks show that semitheoretical model of rock cutting can better reflect the actual situation. Meanwhile, Yasar demonstrates that it is feasible to use small block sample to predict cutting force. Lu et al. [18] conclude that cutting speed does not influence cutting force or fracture characteristics, after researching the fracture features of a free rock plate with three sides fixed and one side free. Pick wear is an important factor that influences cutting efficiency and production cost. To reveal the mechanism of pick wear and the influencing factors, pick wear at different working conditions is analyzed [19,20]. e rock cutting method assisted by a high-pressure water jet is firstly proposed by Michael [21]. In recent years, how to improve cutting ability and avoid pick wear in the hard rock cutting process has become a research hotspot. Stoxreiter et al. [22] combine high-pressure water jet and mechanical excavation technology, and their research shows that high-pressure jetting is an effective approach to break the hard rock. Poláček and Janurová [23] research with an experiment the relationship between the pressure of water jet and rock breaking state. Wang et al. [24] propose a stress released rock cutting method to improve rock cutability via experimental and simulation method. By releasing high stress concentration to lower or none stress condition, rock produces a damage area, thereby the rock cutting ability is improved. However, due to the limitation of mine and rock production conditions, the methods mentioned before are difficult to apply and popularize. In order to reveal the wear mechanism of conical pick, Liu et al. [25] propose finite element method to research cutting form and the characteristics of cutting force influencing pick wear. It is considered that the interference between the back of the pick and the rock is the main factor that causes wear of the pick. eir research shows that rotation of the pick is also the main reason influencing the wear of pick [26]. Li et al. [27] carry out the linear cutting test of rock by using conical pick, pointing out that the main failure mode of pick is plastic failure. Kim et al. [28] research the skew angle influencing the temperature and wear of the conical pick by bench-scale device. e results show that large tip surface area has a significant influence on temperature.
Although scholars have conducted extensive researches on rock cutting mechanism and wear mechanism of conical pick, pick wear is still ubiquitous in hard rock cutting. A diamond saw blade and conical pick combined rock cutting method is proposed in this paper. And the interaction model between the conical pick and rock plate is established to study the fracture morphology and cutting force of rock plate formed by sawing. e factors influencing rock fracture morphology and cutting force are studied. In order to verify the correctness and feasibility of the numerical method, experiments are conducted to test the numerical results.

Proposal of the New Rock Cutting Method.
Conical pick wear is a serious problem in roadway excavation, and it is an unfavorable factor to increase the cost of production. As is well known, diamond is a very hard substance; therefore, diamond saw blade is a common cutting tool in the stone processing industry, which can slit the rock easily. us, the rock plate with multiple free surfaces can be formed. As a result, the cutability of rock can be effectively improved.
en, the cutting head of roadheader will be used to cut the rock plate, as is shown in Figure 1. e schematic drawing of the new roadheader is shown in Figure 2. Compared with conventional roadheader, the saw blade assembly and its drive unit are added. e rock plate with three sides fixed and one side free is formed by cutting with diamond saw blade, as is shown in Figure 3(a). But in rock cutting, it is most possible that the rock plate with one side fixed and three sides free is formed. erefore, this paper takes rock plate with one side fixed and three sides free as the focus of research.

Numerical Model.
e software chosen is very important for getting reliable and scientific results. LS-DYNA is a typical nonlinear dynamic analysis software, with great advantages in collision analysis. In the paper, LS-DYNA is employed to establish the finite element model which consists of two parts, that is, rock and conical pick. Meanwhile, the rock also contains two parts, base rock and rock plate. It should be pointed out that the rock plate and base bock are an integral structure instead of two independent parts, as shown in Figure 4. e hexahedral grid element is employed for rock division. All the nodes at the bottom of base rock are defined as a node group, and all degrees of freedom of the node group are constrained. In order to eliminate stress wave influencing rock fracture and stress distribution, nonreflection boundary condition is added to all of the elements. e conical pick is used to cut the fixed rock with constant velocity. l x , l y , and l z are the width, thickness, and height of rock plate, respectively. In order to research cutting position influencing rock cutting performance, l xp was employed to express the distance between the conical pick and left boundary of the rock plate. In the xoy plane, the angle between velocity of conical pick and the y-axis is defined as cutting angle, expressed by β in Figure 4. JOHNSON_HOLMQUIST_CONCRETE (JHC) [29] and EROSION was assigned to the rock for obtaining the superior rock fracture result. e contact model 2 Advances in Materials Science and Engineering ERODING_SURFACE_TO_SURFACE is assigned between conical pick and rock. Because the finite element model in this paper contains a small number of elements, the conical pick part is a rigid body. With workstations, each group of data takes about three hours to calculate, which is not long. erefore, mass scaling is not used for numerical simulation in this paper. e hourglass energy has a significant influence on the simulation calculation results. In order to avoid the hourglass, in this paper, full integration is used in model calculation, and there is no hourglass energy generated in the simulation calculation process, so the influence of hourglass does not need to be considered in this paper. Taking granite as an example, its main parameters are shown in Table 1.

Experimental Verification.
Mechanical properties are basic for experimental and numerical simulation, especially the Brazilian tensile strength (BTS), uniaxial compressive strength (UCS), elastic modulus (E), and density (D) of rock. An electrohydraulic servo testing machine is used to obtain the uniaxial compressive strength and elastic modulus of rock specimens, as is shown in Figure 5(a), and the stressstrain curve, as is shown in Figure 5(b). e maximum value of the curve is the UCS. e elastic modulus is calculated at the 50% maximum compressive strength; similarly, the average value of three tests is taken. e tensile strength can be obtained by the Brazilian split test, and the density of rock can be obtained from balance. In order to ensure the accuracy of the test results, three groups of tests are conducted for each data. e mechanical properties of marble, Granite, sandstone-1, and sandstone-2 are shown in Table 2.
Only few studies on the fracture characteristics of rock plate cutting are found. In order to verify the correctness of the numerical simulation in the paper, a rock plate cutting testbed was built, as is shown in Figure 6. Figure 7 shows the cutting results of a rock plate of granite with width of 500 mm, height of 200 mm, and thickness of 20 mm. On the same numerical simulation and experimental cutting conditions, the fracture morphology of rock plate is basically the same. e fracture occurs at the junction of rock plate and base rock and in the middle of rock plate. erefore, the numerical simulation method employed in the paper is accurate and feasible.

Influence of Width on Rock Plate Fracture.
e cutting results of granite with different widths are shown in Figure 8. At width of 60 mm and 80 mm, rock plate integral fractures at the junction between base rock and rock plate are shown in Figures 8(a) and 8(b). Under the action of cutting force, element fails at the contact area between rock plate and conical pick. However, when width of rock plate increases to 200 mm and 400 mm, the rock plate is broken into two parts. e main cracks are located in the middle of rock plate as well as at the junction between rock plate and base rock. e numerical results are the same as the experimental results, which demonstrate that the numerical simulation is feasible. As the conical pick moves, the maximum bending moment is produced at the junction of base rock and rock plate, which leads to rock plate broken. When width of rock plate increases, larger deflection and internal force are produced in the width direction of rock plate, so that the rock plate is broken in the middle. e peak cutting force (PCF) variation with width and the fitting curve between them are shown in Figures 9 and  10. It can be seen that PCF increases when the width of rock plate increases. e PCF at 500 mm increases with small value, compared with that of 400 mm. us, with the increase of width, the increasing amplitude of PFC decreases. In the industrial production process, even if the width of rock plate is very large, the PCF will not be very large. e fitting results between width and PCF with a thickness of 20 mm are shown in Table 3. It can be seen from the fitting equation that the exponential-positive relationship exists between PCF and width of rock plate. All of the correlation coefficients are greater than 0.9, so the fitting results are reliable. e statistical P values are less than 0.01, which show that the results are significant.   Advances in Materials Science and Engineering 3 40 mm are shown in Figure 11. At the height of 20 mm, the middle area of rock plate is broken, but the rock plate on both sides remains on the base rock, as is shown in Figure 11(a). At height of 40 mm, rock plate fracture occurs at the middle part of the rock plate; in the meantime, the rock plate on both sides completely is fractured from the base rock, as can been seen in Figure 11(b). When the heights of rock plate increase to 80 mm and 160 mm, the fracture of rock plate shows better integrity, as shown in Figures 11(c) and 11(d). erefore, the conclusion can be drawn: the higher the rock plate is, the easier it is fractured. e variations of PCF with rock plate height and the fitting relationship between them are shown in Figures 12 and  13. PCF decreases when the height of rock plate increases. From the fitting line, shown in Figure 13, and the fitting results, shown in Table 4, a linear negative correlation exists between PCF and height of rock plate. It is different from the exponential-positive relationship between width of rock plate and PCF. In the process of width increase from 200 mm to 500 mm, the slope of the fitting line increases, which demonstrates that at the higher width, PCF changes more obviously with the height of rock plate. e minimum linear correlation coefficient between the height of rock plate and PCF is 0.95 and the P value is less than 0.01, indicating that the relationship is strong and the fitting results are credible. Figure 14 shows the cutting results of rock plates with widths of 400 mm, height of 160, and thicknesses of 16 mm and 30 mm, respectively. ere exists some similarity in fracture morphology between Figures 11(d) and 14(a) with thickness of 16 mm and 20 mm; that is, the cracks of rock plate are very regular. When the thickness of rock plate increases to 30 mm, the fracture morphology changes significantly, and cracks are generated randomly, no longer along with the edge of the rock plate. So, if the rock plate is thicker, the shape of the broken rock plate is more random. e relationships between PCF and thickness of rock plates with width of 400 mm and height of 160 mm are shown in Figure 15 and Table 4. It can be seen from Figure 15 that the PCF increases significantly with the increase of rock plate thickness. e fitting curve and equation show that exponential-positive relationship exists between PCF and thickness of rock plate. Simultaneously, the fitting results are satisfying, as is shown in Table 5. Figures 11(d) and 16 show the cutting results of rock plate with cutting angles of 0°, 10°, and 30°at width of 400 mm, height of 160 mm, and thickness of 20 mm. It can be seen from the comparison of the three figures that there is no difference in the fracture morphology of rock plate with different cutting angles. Figure 17 shows the PCF variation with the cutting angle. With the increase of cutting angle, PCF increases first, with the cutting force reaching the maximum value at the cutting angle of 5°. Subsequently, PCF decreases with the increase of cutting angle, and when the cutting angle increases to 20°,  PCF keeps stable. erefore, the PCF increases and decreases with the increase of cutting angle and then keeps stable. It can be concluded that the cutting angle had little effect on rock fracture morphology, but it had an obvious effect on PCF. Figure 18 shows that the cutting results of rock plate with cutting position of l xp are equal to 1/8l x and 1/4l x , respectively, at width of 500 mm, height of 200 mm, and thickness of 20 mm. At the cutting position of 1/8l x , rock plate fractures regionally, with only the area around conical pick separated from base rock. At the cutting position of 1/4l x , the broken area of rock plate becomes larger, and, simultaneously, rock plate presents a tendency of overall fracture from the base rock. It can be seen that cutting positions have a significant impact on rock fracture. e PCF variation with cutting position is shown in Figure 19, which shows that PCF increases with the increase of cutting position. At the widths of 500 mm and 400 mm, the difference values between different cutting positions are 0.28 kN and 0.33 kN, respectively. erefore, from the perspective of numerical value, the PCF changes little.

Influence of Rock Property on Rock Plate Fracture.
e cutting results with different rocks are shown in Figures 8(d) and 20. It can be clearly seen that the cutting results of rock plates with different properties are obviously   different. e marble with the lowest compressive strength has a smaller fragment size and the largest number of fragments, as is shown in Figure 20(a). With increasing compressive strength, the number of fragmentations decreases, but the size of significant fragment increases. erefore, the compressive strength of rock has obvious influence on the result of rock cutting. Figure 21 shows the PCF increases obviously with the increase of UCS of rock plate. e peak cutting force of marble with compressive strength of 43.8 MPa is 0.82 kN, while as to sandstone-2 with compressive strength of 139.1 MPa, the peak cutting force is 2.51. According to the fitting curve, a strong exponential relationship exists between peak cutting force and rock compressive strength.

Cutting Process of the Rock Plate.
When the conical pick acts upon the middle of rock plate, the cutting process of sandstone-1 with width of 400 mm, height of 120 mm, and thickness of 20 mm is shown in Figure 22, with the cutting results shown in Figure 20(b). With the contact between conical pick and rock plate, the cutting force increases. e elements of rock contacting with the pick are compressed and fail, thus forming a groove on the rock plate, as is shown in Figure 22(a). e groove can effectively increase the contact area between the pick and the rock plate, which causes the rock plate to bear higher cutting force. Figure 22(b) shows that a main crack is generated in the middle of rock plate when cutting force increases. en, based on the main crack, cracks start to expand outward, as is shown in Figure 22(c). Finally, the cracks are connected with each other, further leading to rock fracture.

Conclusion
In order to study the fracture mechanism of rock plate with three sides fixed and one side free, the interaction model between the conical pick and rock plate is built. e PCF and rock fracture characteristic with different cutting parameters and rock parameters are researched.
(1) In order to generate crack and separate fragment in the process of rock plate cutting, ERO-DING_SURFACE_TO_SURFACE and EROSION are added in the numerical model. By comparing with the experimental results, it is shown that the numerical model proposed in the paper is feasible and accurate.  Thickness of rock plate (mm) Figure 15: e relationship between PCF and thickness of rock plate.  (2) e cutting results obtained from different widths, heights, and thicknesses of rock plates show that PCF increases when the thickness and width of rock plate increase, and that it decreases when the height of rock plate increases. e positive exponential relationship exists between PCF and thickness and width of rock plate, and the linear relationship exists between PCF and height of rock plate. Furthermore, thickness has the most significant influence on PCF. (3) PCF also varies with different cutting angles and cutting positions. With the increasing cutting angle, PCF increases and then decreases, and, finally, it tends to be stable. With the increasing cutting position, PCF increases. However, the change of cutting  force with cutting angle and cutting position is relatively small. (4) From the perspective of rock fragmentation, although the width of the rock plate has an effect on the rock fracture results, rock plates are completely broken. Otherwise, lower thickness and higher height benefit rock fracture integrally. Cutting angle has little influence on rock fracture results. e closer the conical pick is to the middle of the rock plate, the more favorable it is for rock plate fracture.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.