# Numerical Simulation of Rock Cutting with a Diamond Sawblade Based on LS-DYNA

Authors: Zhiwen Wang; Qingliang Zeng; Zhenguo Lu; Lirong Wan; Xin Zhang
Journal: Mathematical Problems in Engineering
 (2019)

---

## Abstract

Aiming at the complex nonlinear dynamic problem of cutting hard rock with diamond sawblades, the process of hard rock cutting with a diamond sawblade was studied based on ANSYS/LS-DYNA. A numerical simulation model of a diamond sawblade cutting hard rock with different rotation speeds, feed speeds, and cutting depths was established to study the effects of the cutting parameters on the cutting force and specific cutting energy consumption and on the damage to the rock. The numerical simulation results demonstrated that the feed speed and cutting depth of the diamond saw are quadratically correlated with the cutting force, but the rotation speed is negatively linearly correlated. The damage region of the rock is positively correlated with the feed speed and cutting depth of the diamond sawblade and has a negative correlation with the rotation velocity. The cutting parameters have a great influence on specific cutting energy consumption. Analysis of the relationship among the cutting parameters and the specific cutting energy with multivariate linear fitting indicated that the cutting speed and cutting depth have a great influence on the cutting energy.


---

## Body

## 1. Introduction

Natural stone is widely used in the construction industry, and diamond sawblades are mostly used in the process of cutting natural stone. To meet the basic technological demand, improve rock quality, and reduce the cost, the influence of cutting parameters on the cutting performance in the process of rock cutting with diamond sawblades and optimizing the cutting parameters of the sawblade to improve rock quality and reduce the cost were studied. However, research on diamond sawblades is mainly based on research on specific cutting energy consumption, cutting force, and sawblade wear. Yurdakul and Akdas used statistical methods for the prediction of specific cutting energy values [1] and investigated the effects of the cutting mode, cutting depth, and feed speed on the level of consumed power [2]. Konstanty explained the sawing mechanism and structural properties of diamond impregnation and presented a theoretical model of natural stone sawing by means of diamond-impregnated tools [3]. Yu et al. investigated the characteristics of specific energy in the grinding of granite using diamond abrasives [4]. Sun et al. studied the influence of rock properties on segment wear during the cutting of granite using an innovative diamond frame saw prototype and developed a mathematical model to predict segment wear [5]. Karakurt [6] and Sengun and Altindag [7] investigated the cutting force and specific energy during the diamond sawblade rock-cutting process. Hellstrom et al. performed a theoretical and experimental study of the circular sawing process [8]. Yasitli et al. established a discrete-element model of the sawing process and studied the different peripheral speeds and advanced speed in compliance with the actual cutting operation [9]. Asche et al. showed empirical results about the influence of process parameters on tool wear [10] and used an empirical model to predict the variation in the cutting energy and to guide the selection of cutting conditions [11]. Xu et al. studied the force and energy of grinding granite with circular sawing [12] and researched the force ratio by experimental methods [13]. Turchetta established empirical models to investigate the influence of cutting conditions on cutting force, and cutting energy was related to the shape of the idealized chip thickness [14]. Zeng et al. researched the influences of diamond sawblades with various cutting speeds and rotation speeds on the damage to the rock [15]. Meanwhile, studies of rock fragmentation have been conducted. Lu et al. used the numerical code LS-DYNA(3D) [16] to numerically simulate fragment separation [17]. Evans proposed the classic and the most authoritative theory of cutting rock, which has been accepted by scholars [18], and obtained the optimum line spacing for cutting picks [19]. Tang studied the statistical rock damage constitutive model to investigate rock damage conditions [20].Previous studies have obtained fruitful advancements in diamond sawblades, and studies of diamond sawblades for cutting rocks have mostly studied their wear, cutting force, and specific energy consumption. However, there are fewer studies on rock damage with diamond sawblade cutting. Moreover, most studies were carried out by experimental and theoretical methods, and thus, a more efficient numerical simulation method was used to study sawblade rock cutting in this paper. Through the simulation method to research the process of cutting rock with diamond sawblades, the influences of cutting parameters on the performance of sawblades were studied. The simulation results could be used to guide the production of natural stone.


## 2. Theoretical and Simulation Studies

### 2.1. Theoretical Study

During the cutting of hard rock with a diamond sawblade, the cutting force of the sawblade, which is composed of the normal force, tangential force, and axial force, is very complicated, as shown in Figure1. The normal force is attributable to the impact of the diamond sawblade on the rock, placing the sawblade and stone on the portions of the contact force. The tangential force is the friction against the sawblade when the rock is cut by the sawblade. The axial force forms as a result of the rock fragments compressing the diamond sawblade. The random distribution of rock fragments on both sides of the diamond sawblade, which is the interaction medium between the sawblade and rock, leads to a large fluctuation in the axial force of the diamond sawblade. Meanwhile, the sawblade rapidly rotates, which will remove the rock debris from sawing, causing the curve of the axial force to fluctuate around the axial force equalling zero.Figure 1
The model for cutting hard rock with a diamond sawblade:Fv = vertical force, Fh = horizontal force, Ft = tangential force, Fn = normal force, Vf = feed speed, Vr = rotation speed, and dc = cutting depth.


### 2.2. Mathematical Expression of Cutting Force

A diamond sawblade cuts hard rock with a constant cutting depthdc, constant feed speed Vf, and cutting speed Vs. The rock is cut by a diamond sawblade with the down-cut style. When cutting rock with a diamond sawblade, the tangential and normal forces are not easy to obtain directly, but the horizontal and vertical forces are relatively easy to obtain. The process of solving for the normal and tangential force by horizontal and vertical forces is as follows:The relation of forces can be expressed as follows:(1)Fc=Fh2+Fv2,Fh2+Fv2=Fn2+Ft2,where Fc is the cutting force, Fh is the horizontal force, Fv is the vertical force, Ft is the tangential force, Fn is the normal force, and d is the sawblade diameter. Therefore, the relationship of force is shown in Figure 1, which can be obtained from the geometric relation diagram, as follows:(2)β=cos−11−2dcd,δ=tan−1FhFv.It can be concluded from the graph that the tangential forcesFt and normal force Fn in the rock cutting process can be expressed as follows:(3)Ft=Fcsinα−δ,Fn=Fccosα−δ,Ft=Fvsinα−Fhcosα,Fn=Fvcosα−Fhsinα,α=ε+δ,α=kβ,k=ε+δβ.Among them, the value ofk is determined by the contact position between the diamond sawblade and the rock. Some scholars believe that k = 0.5 [21], but after further study, the calculation of the tangential and normal forces is found to be inaccurate. Through the study of Guo, the research was relatively accurate when k ≈ 0.67 [22], and thus, k=0.7 in this paper.In summary, the calculation formulas of the tangential and normal forces in the process of cutting hard rock with a diamond sawblade can be obtained as follows:(4)Ft=Fvsinkcos−11−2dcd−Fhcoskcos−11−2dcd,(5)Fn=Fvcoskcos−11−2dcd+Fhsinkcos−11−2dcd.


### 2.3. Establishing the Numerical Simulation Model

The diamond sawblade-rock model was established with ANSYS/LS-DYNA(3D). The diameter of the diamond sawblade was 380 mm, and the diamond sawblade had 24 blocks 40 mm in length, 15 mm in height, and 3.4 mm in thickness. The rock was a cuboid block 500 mm in length, 115 mm in width, and 150 mm in height. The 3D 8-node hexahedron element was used to mesh the rock. The meshing of the diamond sawblade and rock was carried out by the meshing tool meshwork of ANSYS/LS-DYNA. The eight-node hexahedron element SOLID164 was adopted to improve the accuracy of the calculation and to accurately simulate the rock-cutting force and crushing process. The method of local thinning for the rock model was used with a smaller mesh in the cutting area and a larger mesh in the noncutting area. The diamond sawblade was set as a rigid body, and the middle part of the noncutting part was removed to reduce the number of grids and improve the calculation efficiency.Constraints were added to the rock model. Full constraints were added to the rock bottom surface, y-axial constraints were added to the front and back surfaces, andx-direction constraints were added to the left and right surfaces of the rock. Nonreflecting boundaries were implemented for all of the rock surfaces except the upper surface, which can prevent stress reflection from affecting the numerical simulation results. Displacement constraints along the y- and z-axes and rotation constraints around the x- and z-axes were added to the diamond sawblade. The sawblade moved uniformly and linearly along the forward direction of the x-axis and rotated around the central axis as shown in Figure 2.Figure 2
Rock-cutting model of the rock with diamond sawblade.The RHT constitutive model, which contains the damage constitutive relation and erosion pattern, was chosen as the hard-rock material model, and the main material parameters of the rock are shown in Table1. The contact form of the diamond sawblade and the rock was defined as the eroding surface and surface (ERODING SURFACE AND SURFACE) contact and automatic universal contact (AUTOMATIC_GENRAL). In the process of cutting rock with a diamond sawblade, the deformity element may appear. It is necessary to control the hourglass energy to prevent the influence of negative volume without the method of mass scaling to ensure the accuracy of the calculation. Based on ensuring the calculation accuracy, the element quantity of the rock model is reduced. The rock was meshed by local refinement, and the cutting area of the diamond sawblade was meshed with smaller elements. The feed speed of the diamond sawblade was set as 0.2, 0.225, 0.25, 0.3, and 0.35 m/min; the rotation speed was set as 1000, 1500, 2000, 2500, and 3500 r/min; and the cutting depth of the diamond sawblade was set as 1, 3, 5, 7, and 9 cm. We set the simulation time as 2.5 minutes, and each step lasts 1 s. The solver in LS-DYNA has no limits for CPU cores, and the results of each step are automatically exported as a d3plot file; the diamond sawblade rock-cutting simulation model is efficiently performed on a workstation for 21 h.Table 1
The main material parameters of rock.
ρ
0
Shear
B0
B1
T1
FC
2760 kg/m3
1.68e10
0.9
0.9
4.539e10
4.8e7
EC
ET
FS∗
FT∗
PEL
PCO
3.0e − 8
3.0e22
0.18
0.1
2.33e7
6e8
E0C
E0T
GC∗
GT∗
D1
D2
3e − 9
3e22
0.53
0.7
0.04
1


## 2.1. Theoretical Study

During the cutting of hard rock with a diamond sawblade, the cutting force of the sawblade, which is composed of the normal force, tangential force, and axial force, is very complicated, as shown in Figure1. The normal force is attributable to the impact of the diamond sawblade on the rock, placing the sawblade and stone on the portions of the contact force. The tangential force is the friction against the sawblade when the rock is cut by the sawblade. The axial force forms as a result of the rock fragments compressing the diamond sawblade. The random distribution of rock fragments on both sides of the diamond sawblade, which is the interaction medium between the sawblade and rock, leads to a large fluctuation in the axial force of the diamond sawblade. Meanwhile, the sawblade rapidly rotates, which will remove the rock debris from sawing, causing the curve of the axial force to fluctuate around the axial force equalling zero.Figure 1
The model for cutting hard rock with a diamond sawblade:Fv = vertical force, Fh = horizontal force, Ft = tangential force, Fn = normal force, Vf = feed speed, Vr = rotation speed, and dc = cutting depth.


## 2.2. Mathematical Expression of Cutting Force

A diamond sawblade cuts hard rock with a constant cutting depthdc, constant feed speed Vf, and cutting speed Vs. The rock is cut by a diamond sawblade with the down-cut style. When cutting rock with a diamond sawblade, the tangential and normal forces are not easy to obtain directly, but the horizontal and vertical forces are relatively easy to obtain. The process of solving for the normal and tangential force by horizontal and vertical forces is as follows:The relation of forces can be expressed as follows:(1)Fc=Fh2+Fv2,Fh2+Fv2=Fn2+Ft2,where Fc is the cutting force, Fh is the horizontal force, Fv is the vertical force, Ft is the tangential force, Fn is the normal force, and d is the sawblade diameter. Therefore, the relationship of force is shown in Figure 1, which can be obtained from the geometric relation diagram, as follows:(2)β=cos−11−2dcd,δ=tan−1FhFv.It can be concluded from the graph that the tangential forcesFt and normal force Fn in the rock cutting process can be expressed as follows:(3)Ft=Fcsinα−δ,Fn=Fccosα−δ,Ft=Fvsinα−Fhcosα,Fn=Fvcosα−Fhsinα,α=ε+δ,α=kβ,k=ε+δβ.Among them, the value ofk is determined by the contact position between the diamond sawblade and the rock. Some scholars believe that k = 0.5 [21], but after further study, the calculation of the tangential and normal forces is found to be inaccurate. Through the study of Guo, the research was relatively accurate when k ≈ 0.67 [22], and thus, k=0.7 in this paper.In summary, the calculation formulas of the tangential and normal forces in the process of cutting hard rock with a diamond sawblade can be obtained as follows:(4)Ft=Fvsinkcos−11−2dcd−Fhcoskcos−11−2dcd,(5)Fn=Fvcoskcos−11−2dcd+Fhsinkcos−11−2dcd.


## 2.3. Establishing the Numerical Simulation Model

The diamond sawblade-rock model was established with ANSYS/LS-DYNA(3D). The diameter of the diamond sawblade was 380 mm, and the diamond sawblade had 24 blocks 40 mm in length, 15 mm in height, and 3.4 mm in thickness. The rock was a cuboid block 500 mm in length, 115 mm in width, and 150 mm in height. The 3D 8-node hexahedron element was used to mesh the rock. The meshing of the diamond sawblade and rock was carried out by the meshing tool meshwork of ANSYS/LS-DYNA. The eight-node hexahedron element SOLID164 was adopted to improve the accuracy of the calculation and to accurately simulate the rock-cutting force and crushing process. The method of local thinning for the rock model was used with a smaller mesh in the cutting area and a larger mesh in the noncutting area. The diamond sawblade was set as a rigid body, and the middle part of the noncutting part was removed to reduce the number of grids and improve the calculation efficiency.Constraints were added to the rock model. Full constraints were added to the rock bottom surface, y-axial constraints were added to the front and back surfaces, andx-direction constraints were added to the left and right surfaces of the rock. Nonreflecting boundaries were implemented for all of the rock surfaces except the upper surface, which can prevent stress reflection from affecting the numerical simulation results. Displacement constraints along the y- and z-axes and rotation constraints around the x- and z-axes were added to the diamond sawblade. The sawblade moved uniformly and linearly along the forward direction of the x-axis and rotated around the central axis as shown in Figure 2.Figure 2
Rock-cutting model of the rock with diamond sawblade.The RHT constitutive model, which contains the damage constitutive relation and erosion pattern, was chosen as the hard-rock material model, and the main material parameters of the rock are shown in Table1. The contact form of the diamond sawblade and the rock was defined as the eroding surface and surface (ERODING SURFACE AND SURFACE) contact and automatic universal contact (AUTOMATIC_GENRAL). In the process of cutting rock with a diamond sawblade, the deformity element may appear. It is necessary to control the hourglass energy to prevent the influence of negative volume without the method of mass scaling to ensure the accuracy of the calculation. Based on ensuring the calculation accuracy, the element quantity of the rock model is reduced. The rock was meshed by local refinement, and the cutting area of the diamond sawblade was meshed with smaller elements. The feed speed of the diamond sawblade was set as 0.2, 0.225, 0.25, 0.3, and 0.35 m/min; the rotation speed was set as 1000, 1500, 2000, 2500, and 3500 r/min; and the cutting depth of the diamond sawblade was set as 1, 3, 5, 7, and 9 cm. We set the simulation time as 2.5 minutes, and each step lasts 1 s. The solver in LS-DYNA has no limits for CPU cores, and the results of each step are automatically exported as a d3plot file; the diamond sawblade rock-cutting simulation model is efficiently performed on a workstation for 21 h.Table 1
The main material parameters of rock.
ρ
0
Shear
B0
B1
T1
FC
2760 kg/m3
1.68e10
0.9
0.9
4.539e10
4.8e7
EC
ET
FS∗
FT∗
PEL
PCO
3.0e − 8
3.0e22
0.18
0.1
2.33e7
6e8
E0C
E0T
GC∗
GT∗
D1
D2
3e − 9
3e22
0.53
0.7
0.04
1


## 3. Numerical Results and Discussion

### 3.1. Grid Independence Validation

To fully illustrate the reliability of the numerical simulation results and the scientific basis of mesh generation, the relation rates of the vertical force and horizontal force of the numerical simulation in different numerical grids were used to judge the grid independence. When the relative change rate is within the set range, it can be concluded that the results are accurate and that the mesh is scientifically based. The number of rock grids has a great influence on the numerical simulation results, so it is necessary to test the independence of the rock grids. Finally, a reasonable mesh density is selected to test the grids of the rock mesh. The relative variation in horizontal and vertical forces under different mesh numbers of the rock is plotted in Figure3.Figure 3
Independence test of the rock mesh density.When the mesh number reached 231340, the relative change rate of the vertical and horizontal forces was less than 0.2%, and the relative change rate was controlled to within 0.2%. Considering the accuracy of the results and computer performance, the number of grids should be more than 231340.


### 3.2. Calibration and Verification

The accuracy of the numerical simulation is verified by comparing the numerical simulation with experimental data. The results of the numerical simulation with a peripheral speed of the diamond sawblade of 30 m/s, a feed speed of 3 m/min and a cutting depth of 15 mm are presented in Figure4(a), and the results of the experiment [13] under the same cutting conditions are shown in Figure 4(b). The average values of the vertical force and horizontal force curves at the stable stage in the experiments are 190.06 N and 137.12 N, respectively, and those of the numerical simulation are 195.26 N and 131.35 N, respectively. Owing to the low frequency of the data acquisition in the experiment, the sudden increase in the force curve at the moment of contact between the sawblade and the rock was not collected. This result indicates that the numerical simulation model accurately simulates the diamond sawblade cutting rock.Figure 4
The results of the numerical simulation and experiment.
(a)
(b)


### 3.3. Failure Mode of Rock

To research the rock failure mode, two test point elements at different positions were selected, and the frequency of the extracting force value improved. A diamond sawblade cutting rock with a high rotation speed and cutting speed causes diverse rock failure modes. The diamond sawblade interacts with the rock, and the cutting force acts on the rock, causing rock failure and fragmentation. To understand the failure forms clearly, the pressure and shear force must be extracted. The pressure is positive, which means that the force is a compressive force; negative pressure denotes a tensile force.Two points in the broken zone were selected to investigate the failure mode of cutting rock with a diamond sawblade, and the results are presented in Figure5. For test element 1 in the broken zone, the pressure force, shear force, and damage value versus time are shown in Figure 6(a). For test point element 1, the damage value increased to 1 within 0.02 ms, the compressive force reached 18.8 MPa, and the maximum shear force reached 12.47 MPa. Therefore, test point 1 is compressive failure. As shown in the Figure 6(b), the tensile force of the test point 2 reached 16.6 MPa, while the shear force reached 11.2 MPa, and the damage value reached 1, which means the element is in tensile failure.Figure 5
The test point elements in the failure mode.Figure 6
The variations in the stress and corresponding damage and shear with time.
(a)
(b)


### 3.4. Influences of the Cutting Parameters on Cutting Force

The cutting force of the diamond sawblade is mainly composed of the tangential force, normal force, and axial force. Diamond sawblades cut hard rock with a feed speed of 0.25 m/min, a rotation speed of 1500 r/min, and a cutting depth of 5 cm. The cutting force, tangential force, normal force, and axial force curves of a diamond sawblade in the cutting process are shown in Figure7. The cutting force, normal force, tangential force, and axial force appear simultaneously when the diamond sawblade is in contact with the hard rock. The normal force and tangential force are the main components of the cutting force of the diamond sawblade. The change tendency of the normal force is similar to that of the cutting force. While the diamond sawblade is touching the rock, the normal force begins to increase rapidly and beyond the fluctuation of stable cutting. The fluctuation ranges of the initial normal and the tangential force are larger than those during steady cutting. The cutting force and normal force of the diamond sawblade decrease rapidly to a stable fluctuation range as the cutting distance increases and fluctuates within a certain range. The axial force appears when the diamond sawblade touches the rock. The fluctuation amplitude of the axial force increases with the cutting distance at the beginning. The axial force fluctuation is stable in a certain range after a certain cutting distance.Figure 7
Force curves of a diamond sawblade cutting hard rock.When the diamond sawblade was cutting rock with a uniform feed speed, rotation speed, and cutting depth, the cutting force, normal force, and tangential force appeared together as the diamond sawblade contacted the hard rock. The main cause of the phenomenon is that the rotation speed and feed speed are higher. The sawblade has a large impact on the rock as the sawblade contacts the rock. In the instant of contact between the sawblade and rock, the peak of the normal force and the peak of the cutting force increase beyond the maximum values of the fluctuation ranges during stable cutting. In actual production, the sawblade contacts the rock with a small feed speed and then increases the feed speed slowly to meet the requirements of the actual work. At the initial stage of cutting, the axial force of the sawblade was smaller, and it increased until it reached a stable fluctuation range as the cutting distance increased. With the distance of the sawblade cutting into the rock increasing, the overlapping area between the sawblade and the rock increased and reached the maximum value, and the axial action increased until reaching the maximum, which caused the axial force to increase and reach the maximum stability fluctuation range as the force changed.The curves of the cutting forces, normal forces, tangential forces, and axial forces when the diamond sawblade cut the hard rock with a rotation speed of 2000 r/min, a cutting depth of 5 cm, and feed speeds of 0.2, 0.225, 0.25, 0.3, and 0.35 m/min are plotted in Figure8. It is obvious that the force of the diamond sawblade increases with increasing feed speed. Additionally, the fluctuation amplitude of force increases with increasing feed speed.Figure 8
Various forces versus distance for different feed speeds: (a) cutting forces, (b) normal forces, (c) tangential forces, and (d) axial forces.
(a)
(b)
(c)
(d)The following is a discussion of the influence of the feed speed on the cutting force of the diamond sawblade cutting rock with a constant rotation speed and cutting depth. The average cutting forces, average normal forces, average tangential forces, and average axial forces were fitted with the feed speeds, as shown Figure9. The average cutting force, average normal force, average tangential force, and average axial force are called the average force. It is clear that the average force increased with increasing feed speed. The relational expression among the average force and the feed speed was y=Ax2+Bx+C. The functional expressions of the feed speed and forces are shown in Figure 9.Figure 9
Relationships between the mean forces and the feed speed: (a) relationship between the mean cutting force and feed speed, (b) relationship between the mean normal force and feed speed, (c) relationship between the mean tangential force and feed speed, and (d) relationship between the mean axial force and feed speed.
(a)
(b)
(c)
(d)To research the influence of the rotation speed on the force of the diamond sawblade cutting rock, the diamond sawblade cut the rock with a feed speed of 0.25 m/min, a cutting depth of 5 cm, and rotation speeds of 1000, 1500, 2000, 2500, and 3500 r/min in the simulations. The force curves of the diamond sawblade versus cutting distance are shown in Figure10, and the cutting force is closely related to the rotation speed. At the instant that the diamond sawblade comes in contact with the rock, the impact leads to an instant increase in the peak force. The main reason is that the rotation speed of the diamond sawblade increases, the cutting speed increases, and the sawblade cutting the hard rock at a constant speed and depth causes the impact force to increase by a larger margin. However, it is obvious that the cutting force, normal force, tangential force, and axial force of the diamond sawblade decreased with increasing rotation speed when the diamond sawblade was cutting stably.Figure 10
Variations in forces with distance: (a) variation in the cutting force with distance for different rotation speeds, (b) variation in the normal force with distance for different rotation speeds, (c) variation in the tangential force with distance for different rotation speeds, and (d) variation in the axial force with distance for different rotation speeds.
(a)
(b)
(c)
(d)With different cutting depths, the rotation speed was fit with the mean of the cutting force, the normal force, the tangential force, and the axial force to obtain a fitting curve and the functional expressions, which are shown in Figure11. The relationship expression between the rotation speed and the force is y=Ax+b. It can be concluded that the rotation speed is linearly related to various forces, and the functional expression and correlation coefficient of each curve are shown in Figure 11.Figure 11
Relationships between the mean force and rotation speed: (a) relationship between the mean cutting force and rotation speed, (b) relationship between the mean normal force and rotation speed, (c) relationship between the mean tangential force and rotation speed, and (d) relationship between the mean axial force and rotation speed.
(a)
(b)
(c)
(d)When the diamond sawblade cuts the rock with cutting depths of 1, 3, 5, 7, and 9 cm, a rotation speed of 2000 r/min and a feed speed of 0.25 m/min, the cutting forces, normal forces, tangential forces, and axial forces on the diamond sawblades versus the cutting distance are shown in Figure12.Figure 12
Variations in forces with distance for different cutting depths: (a) variation in the cutting force with distance for different cutting depths, (b) variation in the normal force with distance for different cutting depths, and (c) variation in the tangential force with distance for different cutting depths.
(a)
(b)
(c)
(d)It is obvious that the force increases with increasing cutting depth. The fluctuation amplitude of the cutting force curve had an increasing trend by comparing the cutting force with different cutting depths. The axial force fluctuation increased noticeably with increasing cutting depth, as seen from the axial force curves of the diamond sawblade. The depth of sawing by the sawblade increases the difficulty of discharging rock fragments and increases the interaction between the diamond sawblade and the rock, causing the fluctuations of the axial force to increase. The cutting depth of the diamond sawblade increased as the cutting force increased, and the fluctuation of the cutting force curve increased. The change trend of the normal force curve is similar to that of the cutting force curve. However, the fluctuation amplitude of the tangential force curve is basically the same, and the fluctuation range is not closely related to the changes in the cutting depth. It can be concluded that the cutting depth has a greater effect on the cutting force. The impact on the cutting force is mainly embodied in the normal force and the tangential force. The increasing normal force and tangential force causes the cutting force to increase. Furthermore, the fluctuation amplitude of the cutting force is mainly influenced by the amplitude of the normal force.The fitting curves between the cutting depth and the cutting force, normal force, tangential force, and axial force are shown in Figure13. In the process of fixed deep cutting, the fitting functional expression between the cutting distance and the cutting force, normal force, tangential force, and axial force is y=Ax2+Bx+C. The relationship expressions and correlation coefficients of the cutting depth and the average force are shown in Figure 13.Figure 13
Relationship between the average force and the cutting depth with different feed speeds: (a) relationship between the average cutting force and cutting depth with different feed speeds, (b) relationship between the average normal force and cutting depth with different feed speeds, (c) relationship between the average tangential force and cutting depth with different feed speeds, and (d) relationship between the average axial force and cutting depth with different feed speeds.
(a)
(b)
(c)
(d)


### 3.5. Influences of the Cutting Parameters on Rock Damage

In the dynamic simulation, the diamond sawblade cut rock with different rotation speeds, cutting depths, and feed speeds. The interaction force of the diamond sawblade and the rock has a certain influence on the hard rock on both sides of the saw seam. After the initial crack of the rock is compressed by the diamond sawblade, a certain amount of damage is produced by transient elastic compression and irreversible plastic deformation. The degree of damage to the rock is closely related to the stress of the rock, so different cutting parameters have a great impact on the degree of rock damage. It can be seen that with the increase in feed speed, the damage zone is increased, and the degree of damage degree is also increased, as shown in the cloud figure of different feeding speeds in Figure14.Figure 14
Damage nephogram for different feed speeds. (a) Feed speed 0.200 m/min. (b) Feed speed 0.225 m/min. (c) Feed speed 0.250 m/min. (d) Feed speed 0.300 m/min. (e) Feed speed 0.350 m/min.
(a)
(b)
(c)
(d)
(e)Rock was cut with the diamond sawblade with a feed speed of 0.25 m/min, a cutting depth of 5 cm and rotation speeds of 1000, 1500, 2000, 2500, and 3500 r/min to study the influence of the rotation speed on the rock damage. The damage zone of the hard rock and the degree of damage were smaller with increasing rotation speed, which is distributed with respect to the sawblade. Damage nephograms with different rotation speeds are presented in Figure15.Figure 15
Damage nephograms for different rotation speeds. (a) Rotation speed 1000 r/min. (b) Rotation speed 1500 r/min. (c) Rotation speed 2000 r/min. (d) Rotation speed 2500 r/min. (e) Rotation speed 3500 r/min.
(a)
(b)
(c)
(d)
(e)Figure16 shows damage nephograms of rock cutting with diamond sawblades at different cutting depths. As the cutting depth increases, the damage zone and the degree of damage increase. When the cutting depth is decreased, an uneven distribution of the damage cloud on both sides of the saw results from the cutting axial forces easily expanding outward.Figure 16
Damage nephograms for different cutting depths. (a) Cutting depth 1 cm. (b) Cutting depth 3 cm. (c) Cutting depth 5 cm. (d) Cutting depth 7 cm. (e) Cutting depth 9 cm.
(a)
(b)
(c)
(d)
(e)


### 3.6. Influences of the Cutting Parameters on Specific Cutting Energy

The product of the average cutting force and cutting distance of a diamond sawblade is defined as the cutting energy consumption of a diamond sawblade, and the specific energy consumption is defined as the ratio of the cutting energy consumption and the volume of rock removed. In this paper, the relational expression was obtained by fitting the specific energy consumption and the feed speed and the rotation speed, as shown in Table2. The influence of the diamond sawblade feed speed and rotation speed on the specific energy consumption is as follows: with increasing feed speed in a quadratic relationship, the diamond sawblade rotation speed linearly decreases, and the specific energy consumption decreases with the increase in cutting depth, as shown in the curves in Figures 17 and 18.Table 2
Regression results to predict the diamond sawblade performance.
Variables
Cutting parameter
Regression equation
Correlation coefficient
F-value
P value
SE-FS
Cutting depth of 1 cm
y = 11.42824x2 − 2.02953x + 0.61463
0.98417
125.30464
5.7e−4
SE-FS
Cutting depth of 3 cm
y
=
11.42824
x
2
−
4.43866
x
+
0.93703
0.99322
293.83917
1.8e−4
SE-FS
Cutting depth of 5 cm
y
=
10.60784
x
2
−
4.34054
x
+
0.38059
0.92829
26.88862
0.00134
SE-FS
Cutting depth of 7 cm
y
=
5.95759
x
2
−
1.98856
x
+
0.585675
0.97888
93.69311
2.7e−4
SE-FS
Cutting depth of 9 cm
y
=
8.99425
x
2
−
3.48346
x
+
0.70192
0.99464
372.24357
9.1e−5
SE-RS
Cutting depth of 1 cm
y
=
−
1.11000
e
−
3
x
+
0.82366
0.92196
48.25547
0.00282
SE-RS
Cutting depth of 3 cm
y
=
−
9.76846
e
−
4
x
+
0.74836
0.92214
48.37552
0.00219
SE-RS
Cutting depth of 5 cm
y
=
−
9.50260
e
−
4
x
+
0.71639
0.96328
105.91909
9.5e−4
SE-RS
Cutting depth of 7 cm
y
=
−
9.35802
e
−
4
x
+
0.65677
0.87846
29.91088
0.00325
SE-RS
Cutting depth of 9 cm
y
=
−
8.95074
e
−
4
x
+
0.58997
0.92743
52.12186
0.00171
SE = specific energy; FS = feed speed; RS = rotation speed.Figure 17
Relationship between the cutting energy consumption and feed speed at different cutting depths.Figure 18
Relationship between the cutting energy consumption and rotation speed at different cutting depths.Not only do the feed speed and the rotation speed affect the cutting energy consumption significantly, but the cutting depth also has a large effect on the cutting energy consumption. Multiple linear regression analysis was conducted to determine the effects of multiple variables on the dependent variables simultaneously. For cutting rock with a diamond sawblade, the specific cutting energy consumption is set as the dependent variable, and the feed speed, rotation speed, and cutting depth are independent variables. Multivariate linear analysis was used to identify the relationship among the specific cutting energy and the 3 independent variables.For convenience of analysis, the specific cutting energy is defined asY, and the independent variables, namely, the feed speed, rotation speed and the cutting depth, are defined as X1, X2, and X3, respectively. The expression of the multivariable linear function is Y = φ + AX1 − BX2 − CX3, where φ is the intercept and A, B, and C are the coefficients of the feed speed, rotation speed, and cutting depth, respectively. The multiple linear function fits the specific cutting energy consumption with respect to the feed speed, the rotation speed, and the cutting depth, and the expression was obtained as follows: Y = 0.41 + 1.679X1 − 0.097X2 − 0.24X3.


## 3.1. Grid Independence Validation

To fully illustrate the reliability of the numerical simulation results and the scientific basis of mesh generation, the relation rates of the vertical force and horizontal force of the numerical simulation in different numerical grids were used to judge the grid independence. When the relative change rate is within the set range, it can be concluded that the results are accurate and that the mesh is scientifically based. The number of rock grids has a great influence on the numerical simulation results, so it is necessary to test the independence of the rock grids. Finally, a reasonable mesh density is selected to test the grids of the rock mesh. The relative variation in horizontal and vertical forces under different mesh numbers of the rock is plotted in Figure3.Figure 3
Independence test of the rock mesh density.When the mesh number reached 231340, the relative change rate of the vertical and horizontal forces was less than 0.2%, and the relative change rate was controlled to within 0.2%. Considering the accuracy of the results and computer performance, the number of grids should be more than 231340.


## 3.2. Calibration and Verification

The accuracy of the numerical simulation is verified by comparing the numerical simulation with experimental data. The results of the numerical simulation with a peripheral speed of the diamond sawblade of 30 m/s, a feed speed of 3 m/min and a cutting depth of 15 mm are presented in Figure4(a), and the results of the experiment [13] under the same cutting conditions are shown in Figure 4(b). The average values of the vertical force and horizontal force curves at the stable stage in the experiments are 190.06 N and 137.12 N, respectively, and those of the numerical simulation are 195.26 N and 131.35 N, respectively. Owing to the low frequency of the data acquisition in the experiment, the sudden increase in the force curve at the moment of contact between the sawblade and the rock was not collected. This result indicates that the numerical simulation model accurately simulates the diamond sawblade cutting rock.Figure 4
The results of the numerical simulation and experiment.
(a)
(b)


## 3.3. Failure Mode of Rock

To research the rock failure mode, two test point elements at different positions were selected, and the frequency of the extracting force value improved. A diamond sawblade cutting rock with a high rotation speed and cutting speed causes diverse rock failure modes. The diamond sawblade interacts with the rock, and the cutting force acts on the rock, causing rock failure and fragmentation. To understand the failure forms clearly, the pressure and shear force must be extracted. The pressure is positive, which means that the force is a compressive force; negative pressure denotes a tensile force.Two points in the broken zone were selected to investigate the failure mode of cutting rock with a diamond sawblade, and the results are presented in Figure5. For test element 1 in the broken zone, the pressure force, shear force, and damage value versus time are shown in Figure 6(a). For test point element 1, the damage value increased to 1 within 0.02 ms, the compressive force reached 18.8 MPa, and the maximum shear force reached 12.47 MPa. Therefore, test point 1 is compressive failure. As shown in the Figure 6(b), the tensile force of the test point 2 reached 16.6 MPa, while the shear force reached 11.2 MPa, and the damage value reached 1, which means the element is in tensile failure.Figure 5
The test point elements in the failure mode.Figure 6
The variations in the stress and corresponding damage and shear with time.
(a)
(b)


## 3.4. Influences of the Cutting Parameters on Cutting Force

The cutting force of the diamond sawblade is mainly composed of the tangential force, normal force, and axial force. Diamond sawblades cut hard rock with a feed speed of 0.25 m/min, a rotation speed of 1500 r/min, and a cutting depth of 5 cm. The cutting force, tangential force, normal force, and axial force curves of a diamond sawblade in the cutting process are shown in Figure7. The cutting force, normal force, tangential force, and axial force appear simultaneously when the diamond sawblade is in contact with the hard rock. The normal force and tangential force are the main components of the cutting force of the diamond sawblade. The change tendency of the normal force is similar to that of the cutting force. While the diamond sawblade is touching the rock, the normal force begins to increase rapidly and beyond the fluctuation of stable cutting. The fluctuation ranges of the initial normal and the tangential force are larger than those during steady cutting. The cutting force and normal force of the diamond sawblade decrease rapidly to a stable fluctuation range as the cutting distance increases and fluctuates within a certain range. The axial force appears when the diamond sawblade touches the rock. The fluctuation amplitude of the axial force increases with the cutting distance at the beginning. The axial force fluctuation is stable in a certain range after a certain cutting distance.Figure 7
Force curves of a diamond sawblade cutting hard rock.When the diamond sawblade was cutting rock with a uniform feed speed, rotation speed, and cutting depth, the cutting force, normal force, and tangential force appeared together as the diamond sawblade contacted the hard rock. The main cause of the phenomenon is that the rotation speed and feed speed are higher. The sawblade has a large impact on the rock as the sawblade contacts the rock. In the instant of contact between the sawblade and rock, the peak of the normal force and the peak of the cutting force increase beyond the maximum values of the fluctuation ranges during stable cutting. In actual production, the sawblade contacts the rock with a small feed speed and then increases the feed speed slowly to meet the requirements of the actual work. At the initial stage of cutting, the axial force of the sawblade was smaller, and it increased until it reached a stable fluctuation range as the cutting distance increased. With the distance of the sawblade cutting into the rock increasing, the overlapping area between the sawblade and the rock increased and reached the maximum value, and the axial action increased until reaching the maximum, which caused the axial force to increase and reach the maximum stability fluctuation range as the force changed.The curves of the cutting forces, normal forces, tangential forces, and axial forces when the diamond sawblade cut the hard rock with a rotation speed of 2000 r/min, a cutting depth of 5 cm, and feed speeds of 0.2, 0.225, 0.25, 0.3, and 0.35 m/min are plotted in Figure8. It is obvious that the force of the diamond sawblade increases with increasing feed speed. Additionally, the fluctuation amplitude of force increases with increasing feed speed.Figure 8
Various forces versus distance for different feed speeds: (a) cutting forces, (b) normal forces, (c) tangential forces, and (d) axial forces.
(a)
(b)
(c)
(d)The following is a discussion of the influence of the feed speed on the cutting force of the diamond sawblade cutting rock with a constant rotation speed and cutting depth. The average cutting forces, average normal forces, average tangential forces, and average axial forces were fitted with the feed speeds, as shown Figure9. The average cutting force, average normal force, average tangential force, and average axial force are called the average force. It is clear that the average force increased with increasing feed speed. The relational expression among the average force and the feed speed was y=Ax2+Bx+C. The functional expressions of the feed speed and forces are shown in Figure 9.Figure 9
Relationships between the mean forces and the feed speed: (a) relationship between the mean cutting force and feed speed, (b) relationship between the mean normal force and feed speed, (c) relationship between the mean tangential force and feed speed, and (d) relationship between the mean axial force and feed speed.
(a)
(b)
(c)
(d)To research the influence of the rotation speed on the force of the diamond sawblade cutting rock, the diamond sawblade cut the rock with a feed speed of 0.25 m/min, a cutting depth of 5 cm, and rotation speeds of 1000, 1500, 2000, 2500, and 3500 r/min in the simulations. The force curves of the diamond sawblade versus cutting distance are shown in Figure10, and the cutting force is closely related to the rotation speed. At the instant that the diamond sawblade comes in contact with the rock, the impact leads to an instant increase in the peak force. The main reason is that the rotation speed of the diamond sawblade increases, the cutting speed increases, and the sawblade cutting the hard rock at a constant speed and depth causes the impact force to increase by a larger margin. However, it is obvious that the cutting force, normal force, tangential force, and axial force of the diamond sawblade decreased with increasing rotation speed when the diamond sawblade was cutting stably.Figure 10
Variations in forces with distance: (a) variation in the cutting force with distance for different rotation speeds, (b) variation in the normal force with distance for different rotation speeds, (c) variation in the tangential force with distance for different rotation speeds, and (d) variation in the axial force with distance for different rotation speeds.
(a)
(b)
(c)
(d)With different cutting depths, the rotation speed was fit with the mean of the cutting force, the normal force, the tangential force, and the axial force to obtain a fitting curve and the functional expressions, which are shown in Figure11. The relationship expression between the rotation speed and the force is y=Ax+b. It can be concluded that the rotation speed is linearly related to various forces, and the functional expression and correlation coefficient of each curve are shown in Figure 11.Figure 11
Relationships between the mean force and rotation speed: (a) relationship between the mean cutting force and rotation speed, (b) relationship between the mean normal force and rotation speed, (c) relationship between the mean tangential force and rotation speed, and (d) relationship between the mean axial force and rotation speed.
(a)
(b)
(c)
(d)When the diamond sawblade cuts the rock with cutting depths of 1, 3, 5, 7, and 9 cm, a rotation speed of 2000 r/min and a feed speed of 0.25 m/min, the cutting forces, normal forces, tangential forces, and axial forces on the diamond sawblades versus the cutting distance are shown in Figure12.Figure 12
Variations in forces with distance for different cutting depths: (a) variation in the cutting force with distance for different cutting depths, (b) variation in the normal force with distance for different cutting depths, and (c) variation in the tangential force with distance for different cutting depths.
(a)
(b)
(c)
(d)It is obvious that the force increases with increasing cutting depth. The fluctuation amplitude of the cutting force curve had an increasing trend by comparing the cutting force with different cutting depths. The axial force fluctuation increased noticeably with increasing cutting depth, as seen from the axial force curves of the diamond sawblade. The depth of sawing by the sawblade increases the difficulty of discharging rock fragments and increases the interaction between the diamond sawblade and the rock, causing the fluctuations of the axial force to increase. The cutting depth of the diamond sawblade increased as the cutting force increased, and the fluctuation of the cutting force curve increased. The change trend of the normal force curve is similar to that of the cutting force curve. However, the fluctuation amplitude of the tangential force curve is basically the same, and the fluctuation range is not closely related to the changes in the cutting depth. It can be concluded that the cutting depth has a greater effect on the cutting force. The impact on the cutting force is mainly embodied in the normal force and the tangential force. The increasing normal force and tangential force causes the cutting force to increase. Furthermore, the fluctuation amplitude of the cutting force is mainly influenced by the amplitude of the normal force.The fitting curves between the cutting depth and the cutting force, normal force, tangential force, and axial force are shown in Figure13. In the process of fixed deep cutting, the fitting functional expression between the cutting distance and the cutting force, normal force, tangential force, and axial force is y=Ax2+Bx+C. The relationship expressions and correlation coefficients of the cutting depth and the average force are shown in Figure 13.Figure 13
Relationship between the average force and the cutting depth with different feed speeds: (a) relationship between the average cutting force and cutting depth with different feed speeds, (b) relationship between the average normal force and cutting depth with different feed speeds, (c) relationship between the average tangential force and cutting depth with different feed speeds, and (d) relationship between the average axial force and cutting depth with different feed speeds.
(a)
(b)
(c)
(d)


## 3.5. Influences of the Cutting Parameters on Rock Damage

In the dynamic simulation, the diamond sawblade cut rock with different rotation speeds, cutting depths, and feed speeds. The interaction force of the diamond sawblade and the rock has a certain influence on the hard rock on both sides of the saw seam. After the initial crack of the rock is compressed by the diamond sawblade, a certain amount of damage is produced by transient elastic compression and irreversible plastic deformation. The degree of damage to the rock is closely related to the stress of the rock, so different cutting parameters have a great impact on the degree of rock damage. It can be seen that with the increase in feed speed, the damage zone is increased, and the degree of damage degree is also increased, as shown in the cloud figure of different feeding speeds in Figure14.Figure 14
Damage nephogram for different feed speeds. (a) Feed speed 0.200 m/min. (b) Feed speed 0.225 m/min. (c) Feed speed 0.250 m/min. (d) Feed speed 0.300 m/min. (e) Feed speed 0.350 m/min.
(a)
(b)
(c)
(d)
(e)Rock was cut with the diamond sawblade with a feed speed of 0.25 m/min, a cutting depth of 5 cm and rotation speeds of 1000, 1500, 2000, 2500, and 3500 r/min to study the influence of the rotation speed on the rock damage. The damage zone of the hard rock and the degree of damage were smaller with increasing rotation speed, which is distributed with respect to the sawblade. Damage nephograms with different rotation speeds are presented in Figure15.Figure 15
Damage nephograms for different rotation speeds. (a) Rotation speed 1000 r/min. (b) Rotation speed 1500 r/min. (c) Rotation speed 2000 r/min. (d) Rotation speed 2500 r/min. (e) Rotation speed 3500 r/min.
(a)
(b)
(c)
(d)
(e)Figure16 shows damage nephograms of rock cutting with diamond sawblades at different cutting depths. As the cutting depth increases, the damage zone and the degree of damage increase. When the cutting depth is decreased, an uneven distribution of the damage cloud on both sides of the saw results from the cutting axial forces easily expanding outward.Figure 16
Damage nephograms for different cutting depths. (a) Cutting depth 1 cm. (b) Cutting depth 3 cm. (c) Cutting depth 5 cm. (d) Cutting depth 7 cm. (e) Cutting depth 9 cm.
(a)
(b)
(c)
(d)
(e)


## 3.6. Influences of the Cutting Parameters on Specific Cutting Energy

The product of the average cutting force and cutting distance of a diamond sawblade is defined as the cutting energy consumption of a diamond sawblade, and the specific energy consumption is defined as the ratio of the cutting energy consumption and the volume of rock removed. In this paper, the relational expression was obtained by fitting the specific energy consumption and the feed speed and the rotation speed, as shown in Table2. The influence of the diamond sawblade feed speed and rotation speed on the specific energy consumption is as follows: with increasing feed speed in a quadratic relationship, the diamond sawblade rotation speed linearly decreases, and the specific energy consumption decreases with the increase in cutting depth, as shown in the curves in Figures 17 and 18.Table 2
Regression results to predict the diamond sawblade performance.
Variables
Cutting parameter
Regression equation
Correlation coefficient
F-value
P value
SE-FS
Cutting depth of 1 cm
y = 11.42824x2 − 2.02953x + 0.61463
0.98417
125.30464
5.7e−4
SE-FS
Cutting depth of 3 cm
y
=
11.42824
x
2
−
4.43866
x
+
0.93703
0.99322
293.83917
1.8e−4
SE-FS
Cutting depth of 5 cm
y
=
10.60784
x
2
−
4.34054
x
+
0.38059
0.92829
26.88862
0.00134
SE-FS
Cutting depth of 7 cm
y
=
5.95759
x
2
−
1.98856
x
+
0.585675
0.97888
93.69311
2.7e−4
SE-FS
Cutting depth of 9 cm
y
=
8.99425
x
2
−
3.48346
x
+
0.70192
0.99464
372.24357
9.1e−5
SE-RS
Cutting depth of 1 cm
y
=
−
1.11000
e
−
3
x
+
0.82366
0.92196
48.25547
0.00282
SE-RS
Cutting depth of 3 cm
y
=
−
9.76846
e
−
4
x
+
0.74836
0.92214
48.37552
0.00219
SE-RS
Cutting depth of 5 cm
y
=
−
9.50260
e
−
4
x
+
0.71639
0.96328
105.91909
9.5e−4
SE-RS
Cutting depth of 7 cm
y
=
−
9.35802
e
−
4
x
+
0.65677
0.87846
29.91088
0.00325
SE-RS
Cutting depth of 9 cm
y
=
−
8.95074
e
−
4
x
+
0.58997
0.92743
52.12186
0.00171
SE = specific energy; FS = feed speed; RS = rotation speed.Figure 17
Relationship between the cutting energy consumption and feed speed at different cutting depths.Figure 18
Relationship between the cutting energy consumption and rotation speed at different cutting depths.Not only do the feed speed and the rotation speed affect the cutting energy consumption significantly, but the cutting depth also has a large effect on the cutting energy consumption. Multiple linear regression analysis was conducted to determine the effects of multiple variables on the dependent variables simultaneously. For cutting rock with a diamond sawblade, the specific cutting energy consumption is set as the dependent variable, and the feed speed, rotation speed, and cutting depth are independent variables. Multivariate linear analysis was used to identify the relationship among the specific cutting energy and the 3 independent variables.For convenience of analysis, the specific cutting energy is defined asY, and the independent variables, namely, the feed speed, rotation speed and the cutting depth, are defined as X1, X2, and X3, respectively. The expression of the multivariable linear function is Y = φ + AX1 − BX2 − CX3, where φ is the intercept and A, B, and C are the coefficients of the feed speed, rotation speed, and cutting depth, respectively. The multiple linear function fits the specific cutting energy consumption with respect to the feed speed, the rotation speed, and the cutting depth, and the expression was obtained as follows: Y = 0.41 + 1.679X1 − 0.097X2 − 0.24X3.


## 4. Conclusions

The finite element software LS-DYNA was used to study the process of diamond sawblades cutting rock. The rock in the numerical simulation model determined the material model, which included the damage failure criterion, to research the force and rock damage in the process of cutting rock with a diamond sawblade. The influences of feed speed, rotation speed, and cutting depth on forces, rock damage, and specific cutting energy were studied.(1)
The force of the diamond sawblade increases with the increase in feed speed and cutting depth and decreases with the increase in rotation speed during the cutting of hard rock with a diamond sawblade. The fluctuation of the forces increases with increasing feed speed and cutting depth and decreases with increasing rotation speed.(2)
The degree of rock damage increases with increasing cutting depth and feed speed of the diamond sawblade. The increase in rotation speed of the diamond sawblade causes the damage zone of the rock to decrease.(3)
The specific cutting energy consumption increases to the second power with the feed speed, and the rock specific cutting energy consumption decreases linearly with the rotation speed. However, the cutting energy consumption decreases to a certain extent with increasing cutting depth of the diamond sawblade.



---

*Source: 6462909-2019-09-25.xml*