An Analytical Model for Rock Cutting with a Chisel Pick of the Cutter Suction Dredger

Cutter suction dredgers are important pieces of rock excavation equipment in port and waterway construction. It is valuable but difficult to properly estimate the cutting force on the chisel pick of the cutter suction dredger. In this paper, an analytical model, called the crushed zone expansion induced tensile failure model (CEIT model), is proposed for rock cutting with a chisel pick in order to predict the peak cutting force (Fc) more accurately. First, a review of the existing models for rock cutting with a chisel pick is presented. Next, based on the tensile breakage theory, cavity expansion theory and some hypotheses, the mathematical formula of the CEIT model is obtained. Different from that in the previous models, the effect of the rock on both sides of the chisel pick on Fc, defined as the sidewall effect is considered in the CEIT model. Then, the predicted Fc by the CEIT model is compared with the predicted Fc by existing theoretical models and experimental results to check the validity of the CEIT model. The results show that the CEIT model can well capture the relationships of Fc to the cutting parameters, including cutting width, cutting depth, and rake angle, and can predict the experimental results much better than the existing models. Finally, the sidewall effect and its influence factors according to the CEIT model are discussed.


Introduction
Rock excavation with tools is widely used for mining, tunneling, drilling, and dredging engineering. The rock excavation devices, such as the polycrystalline diamond compact cutter (PDC cutter) and the cutter head of the cutter suction dredger (CSD), break the rock mainly through the cutting action. The cutter head of the CSD and the rock cutting with a chisel pick (tooth of the cutter head) are shown in Figure 1. With the development of the CSDs, they are widely adopted for rock excavation in port and waterway construction. From a large-scale perspective, rock excavation is the interaction between the rock masses and the excavation devices. However, in practical applications, the interaction zone between every single chisel pick and the rock masses is very small, and the chisel pick is frequently required to break the intact rock, rather than just separate the rock masses along the discontinuities, especially when the rock masses contain very few discontinuities. Considering the fact that the strength of the intact rock is larger than strength of the rock masses, cutting the intact rock requires a larger cutting force. Thus, rock cutting usually refers to the interactions between the cutters and the intact cutters and the intact rock rather than the rock masses, and the rock strength adopted in rock cutting refers to the strength of the intact rock.
In order to enhance cutting performance and machine life, it is very important to reveal the fracture mechanism of rocks and estimate the cutting force in the rock cutting process. The mechanism of rock cutting with PDC cutters has been investigated in numerous pieces of scientific literature [1][2][3][4][5][6][7][8][9][10][11]. Since the cutting depth is shallow, typically less than 1 mm for a medium strength sandstone [12], and the rake angle is negative (or the back rake angle is positive), the ductile fracture regime usually happens in rock cutting with PDC cutters [12][13][14]. In the ductile fracture regime, the cutting force is proportional to the sectional area of the groove [6,12,13,15]. The ratio of the cutting force to the sectional area of the groove is the intrinsic energy, which is found well correlated with the uniaxial compressive strength (UCS) [4,12,16]. Therefore, the cutting force can be estimated by multiplying the sectional area of the groove and the UCS of the rock. However, this method is not suitable for predicting the cutting force on the teeth (chisel picks or conical picks) mounted on the cutter head of the CSD. For rock cutting with the CSD, positive rake angles and large cutting depths, mostly larger than 10 mm to obtain high dredging production, are used [17], as presented in Figure  1b. In this situation, the rock will fracture in the brittle regime rather than the ductile regime [14,18,19]. In the brittle fracture regime, the cutting force is affected by the rake angle, the cutting depth, the shape of the cutter, and the rock strength [20][21][22][23][24][25][26][27]. Evans [20][21][22] believed that rock cutting with a large positive rake angle is governed by tensile failure and the cutting force is related to the tensile strength of the rock. While Nishimatsu [26] believed that rock cutting in the brittle fracture regime is governed by shear failure and the cutting force is related to the shear strength of the rock. Both Evans and Nishimatsu developed theoretical models for rock cutting problems. Their models are widely cited for comparisons [28][29][30][31][32][33][34] and for developing new models [23][24][25]35]. However, Evans' model [20] and Nishimatsu's model [26] for rock cutting with a chisel pick are developed by two-dimensional analysis. They are only suitable for practical rock cutting in shallow cutting depths (the ratio of the cutting width to cutting depth is very large). When the cutting depth is close to or larger than the cutting width, the rock on both sides of the chisel pick has a great effect on the cutting force, as studied in [15,36,37]. This effect is defined as the sidewall effect in this work. The cutting force may increase in a parabolic form rather than a linear form with the increasing cutting depth because of the sidewall effect, both in the ductile and brittle fracture regimes, according to the experimental results, when the ratio of the cutting width to cutting depth is small enough [15,27,36,38]. Therefore, the sidewall effect cannot be ignored during estimating the cutting force on the chisel pick of the CSD because the cutting depth is usually larger than the cutting width.
Even though dredging operations with a CSD are performed under water pressure, a limited water depth does not cause hyperbaric conditions [39]. In addition, the results in [40,41] show that However, Evans' model [20] and Nishimatsu's model [26] for rock cutting with a chisel pick are developed by two-dimensional analysis. They are only suitable for practical rock cutting in shallow cutting depths (the ratio of the cutting width to cutting depth is very large). When the cutting depth is close to or larger than the cutting width, the rock on both sides of the chisel pick has a great effect on the cutting force, as studied in [15,36,37]. This effect is defined as the sidewall effect in this work. The cutting force may increase in a parabolic form rather than a linear form with the increasing cutting depth because of the sidewall effect, both in the ductile and brittle fracture regimes, according to the experimental results, when the ratio of the cutting width to cutting depth is small enough [15,27,36,38]. Therefore, the sidewall effect cannot be ignored during estimating the cutting force on the chisel pick of the CSD because the cutting depth is usually larger than the cutting width.
Even though dredging operations with a CSD are performed under water pressure, a limited water depth does not cause hyperbaric conditions [39]. In addition, the results in [40,41] show that even in deep water conditions, the water pressure has a very small effect on the cutting force in the steady stage of the dredging operation. Therefore, it is acceptable to ignore the effect of the water pressure on the cutting force during rock dredging with a CSD. Then, the cutting action of the chisel pick or conical pick mounted on the CSD is similar to the cutting actions with a continuous miner [38] and a road-header [42]. The rock cutting process with a chisel pick or conical pick of the CSD can be characterized by an increase in the cutting force as the picks penetrate into the rock. At the same time, a crushed zone is generated in front of the pick tip [31,[43][44][45][46]. When the cutting force reaches a critical value, a rock chip is produced with an instantaneous reduction in cutting force. As the pick cuts forward, it will penetrate into the rock again and repeat the chip formation process. The shape of the diagram of the cutting force versus cutting distance is zigzag [31,38,[47][48][49][50][51]. The crushed zone generated near the pick tip, as shown in Figure 1b, is similar to the core region and plastic region of the cavity expansion model [52][53][54][55][56][57]. Therefore, the cavity expansion theory was used for analyzing the fracture mechanism of rock indentation by tools in [53,54,57]. As the pick penetrates into the rock, the size of the crushed zone increases, which will lead to the initiation and propagation of tensile cracks and finally the rock fragmentation in rock indentation. Rock cutting with picks at a positive rake angle is actually a series of indentations on the rock edge. In this situation, the tensile cracks, initiating near the elasto-plastic interface, can easier propagate to the free surface of the rock because the rock edge is unconfined. Macro tensile cracks can be visualized in the rock edge indentation tests when the cutting depth is large [58]. A few models were proposed for predicting the peak cutting force on the picks based on tensile breakage theory [20,22,59,60]. The first indentation cycle of the rock cutting process was considered to obtain these analytical models. As the shape of the chip produced in the first indentation cycle is relatively regular, the chips produced at different cutting depths are similar [58,61].
Based on the above context, an analytical model is proposed in this paper for estimating the peak cutting force on the chisel pick in cutting the intact rock. The new proposed model is based on the tensile breakage theory, similar to Evans' model and the assumption that the tensile failure is induced by the expansion of the crushed zone. The new model is called crushed zone expansion induced tensile failure model (CEIT model). The CEIT model will also count the sidewall effect to improve its accuracy in predicting the cutting force in rock cutting with a smaller ratio of the cutting width to cutting depth. First, a brief review of the existing models for rock cutting with a chisel pick is presented. Next, the mechanical model and mathematical formula of the CEIT model are developed based on the tensile fracture theory and cavity expansion theory. Then, the CEIT model is compared with the existing theoretical models and is verified by published experimental data. Finally, some discussions and conclusions are presented about the CEIT model.

A Brief Review of the Models for Rock Cutting with Chisel Picks
In this paper, we focus on predicting models for rock cutting with chisel picks under the brittle fracture regime. Since Evans [20] introduced the theory of plowing of coal by a sharp wedge into rock cutting, the theory for rock cutting with chisel picks has been developed by many researchers. Numerous models for estimating the cutting forces on a chisel pick have also been proposed, including the theoretical models [20,21,26], semi-empirical models [23,33], and regression models [32,62]. In this section, some of these models are presented and discussed.

Evans' Tensile Breakage Model
Evans [20] assumed that the rock fails along a circular arc due to tensile fracture when a symmetrical chisel pick penetrates normally into the rock edge, as shown in Figure 2a. Then, the tensile force, T, on the tensile failure path is calculated as follows: where σ t is the tensile strength of the intact rock, w is the cutting width, Rc is the radius of the arcuate failure path, and ψ is the breakout angle. Evans also assumed that the size of the interaction zone between the chisel pick and rock is very small when the rock chip is raptured from the rock mass compared with the cutting depth, and there is a concentrated interaction force, S, between the rock chip and rock mass. Then, the relationship between P and T can be obtained by taking moments about point B, as shown in Figure 2a.
where P is the force acting on the rock chip by the chisel pick, as shown in Figure 2a. Considering the geometric relationship in Figure 2a, Rc can be expressed by the cutting depth (d) as: Rc = 0.5d/sin 2 ψ. Thus, the peak cutting force, Fc, can be calculated as follows: The value of ψ can be obtained with the minimum work hypothesis. In the mathematical expression, it is obtained by solving equation ∂Fc/∂ψ = 0. Finally, the peak cutting force on the chisel pick penetrating normally into the rock is calculated as If the friction between the chisel pick and rock is considered, Fc is calculated as where φ is the friction angle between the chisel pick and rock. zone between the chisel pick and rock is very small when the rock chip is raptured from the rock mass compared with the cutting depth, and there is a concentrated interaction force, S, between the rock chip and rock mass. Then, the relationship between P and T can be obtained by taking moments about point B, as shown in Figure 2a.
where P is the force acting on the rock chip by the chisel pick, as shown in Figure  Thus, the peak cutting force, Fc, can be calculated as follows: The value of ψ can be obtained with the minimum work hypothesis. In the mathematical expression, it is obtained by solving equation If the friction between the chisel pick and rock is considered, Fc is calculated as where φ is the friction angle between the chisel pick and rock. In practical applications, continuous rock cutting is more universal than rock indentation. For rock indentation, the clearance angle, β , should be positive, as shown in Figure 2b. In this situation, the cutting action is asymmetrical. Thus, another form of Equation (5) for rock indentation is used for predicting the peak cutting force in continuous rock cutting. The modified form is as follows [63]: In practical applications, continuous rock cutting is more universal than rock indentation. For rock indentation, the clearance angle, β, should be positive, as shown in Figure 2b. In this situation, the cutting action is asymmetrical. Thus, another form of Equation (5) for rock indentation is used for predicting the peak cutting force in continuous rock cutting. The modified form is as follows [63]: where α is the rake angle in rock cutting.

Nishimatsu's Shear Breakage Model
Some researchers believe the rock fails due to shear breakage rather than tensile breakage in rock cutting. Based on the shear fracture theory, Nishimatsu [26] proposed a theoretical model for estimating the cutting force on the chisel pick. Nishimatsu assumed that the shear failure path during rock cutting is a straight line (or a plane), as shown in Figure 3, and that the shear stress (τ) and normal stress (σ) on the shear failure path distribute as follows: where τ A and σ A are the shear stress and normal stress on the shear failure path near point A in Figure 3, λ is the distance between the calculated point and point A along the AB, n is the stress distribution factor. Then, the resultant shear force (Fs) and resultant normal force (Fp) acting on the rock chip by the rock mass can be obtained by integrating τ and σ from point A to B.
where w is the cutting width and d is the cutting depth in rock cutting.

Nishimatsu's Shear Breakage Model
Some researchers believe the rock fails due to shear breakage rather than tensile breakage in rock cutting. Based on the shear fracture theory, Nishimatsu [26] proposed a theoretical model for estimating the cutting force on the chisel pick. Nishimatsu assumed that the shear failure path during rock cutting is a straight line (or a plane), as shown in Figure 3, and that the shear stress ( τ ) and normal stress ( σ ) on the shear failure path distribute as follows: where A τ and A σ are the shear stress and normal stress on the shear failure path near point A in where w is the cutting width and d is the cutting depth in rock cutting. On the other hand, the relationships of Fs and Fp to the resultant force on the chisel pick, F, can be obtained according to the forces equilibrium equations.

cos(
) and sin( ) In addition, Nishimatsu assumed that rock failure takes place when the maximum stresses on line AB in Figure 3 satisfy the Mohr-Coulomb failure criterion. On the other hand, the relationships of Fs and Fp to the resultant force on the chisel pick, F, can be obtained according to the forces equilibrium equations.
In addition, Nishimatsu assumed that rock failure takes place when the maximum stresses on line AB in Figure 3 satisfy the Mohr-Coulomb failure criterion.
where τ 0 is the shear strength of rock and ϕ is the internal friction angle. Combining Equations (9) to (11), and the minimum work hypothesis, the peak cutting force on the chisel pick, Fc, can be obtained.
where n = 11.3 − 0.18 α and the unit of the rake angle α is degree.

Semi-Empirical Models
Both Evans' model and Nishimatsu's model are developed based on a single failure type, tensile failure or shear failure. Thus, the predicted peak cutting force is only affected by the tensile strength or shear strength of rock. However, Goktan [23] believed that the cutting force is affected by both the tensile strength and compressive strength of the rock. He developed a semi-empirical model based on the experimental results of the rock cutting tests of Bilgin [28].
where σ c is the uniaxial compressive strength of the rock. This model is only valid for cutting rock with the value of σ c between 80 and 180 MPa, and the value of the rake angle α between −20 degrees to 10 degrees. Ouyang [33] believed that the influence of the rock strength on rock cutting changes with increasing the rake angle. He claimed that the tensile strength plays a dominant role in the cutting force when the rake angle is very large and that the influence of the compressive strength on the cutting force will increase with decreasing the rake angle. A semi-empirical model was developed in [33] based on these hypotheses and the experimental results in [28].
where g(α, φ) = 0.32 + 0.075 ln(α − φ + 1.05). The unit of the cutting depth d is millimeter in Equation (14). Both Goktan's model, Equation (13), and Ouyang's model, Equation (14), can well predict the experimental results of rock cutting with a chisel pick in [28], as discussed in [33]. However, their validity for other rock cutting cases needs further verification. Comparisons of the models mentioned above for estimating the peak cutting force in rock cutting with a chisel pick are listed in Table 1.
In conclusion, the theoretical models mentioned above were developed through two-dimensional analysis and ignored the sidewall effect. Thus, they usually underestimate the peak cutting force on the chisel pick in practical applications. The semi-empirical models are suitable only for quite limited cutting conditions and they are not suitable for rock cutting with the chisel pick of the CSD. In rock cutting with the CSD, the cutting depth is usually larger than 10 mm, and the rake angle is larger than 20 degrees, and the uniaxial compressive strength of excavating rock is usually less than 80 MPa. Therefore, it is of great value to develop a model to predict properly the cutting force on the chisel pick of the CSD. where: σ t is the tensile strength of rock, φ is the friction angle between the chisel pick and rock, w is the cutting width, d is the cutting depth, α is the rake angle, τ 0 is the shear strength of rock, ϕ is the internal friction angle, n is the stress distribution factor, σ c is the uniaxial compressive strength of rock.

Conception of the CEIT Model
Many studies show that the contact zone between the pick and rock is very small, compared with the size of the rock chip produced, in rock cutting under the brittle fracture regime [20,22,58,61]. It is widely accepted that the small contact zone will lead to plastic failure of rock near the tip of the pick and form a crushed zone. The pressure within the crushed zone increases as the pick cuts forward and causes the initiation of the cracks on the boundary of the crushed zone similar to that in cavity expansion theory [52][53][54][55][56][57]. Thus, we propose an analytical model based on Evans' tensile breakage model.
As illustrated in Figure 4, we first make the following assumptions, similar to those in Evans' model. (1) The rock chip fails along a circular arc on the vertical symmetrical plane of the chisel pick due to tensile fracture. (2) The tensile failure path has a horizontal tangent at point A. (3) The rock chip is in the moment equilibrium condition about point B at the time when it is cut from the rock mass. In addition, a few assumptions are made in our CEIT model, listed as follows: (a) a crushed zone appears near the tip of the chisel pick and its size increases gradually as the chisel pick penetrates into the rock and reaches a maximum value when a rock chip fails. (b) The crushed zone has a fanlike cross-section with a radius of r, which is very small compared with the cutting depth and the cutting width (or the width of the chisel pick), and a center on the tip of the chisel pick. (c) The cutting action of the chisel pick is transmitted to the rock chip by the crushed zone where the rock material is fully crushed and the pressure is uniform. (d) The rock chip is torn off the rock mass by the resultant force acting on it from the crushed zone. (e) The rock chip formed in the first indentation cycle is composed of three parts, as shown in Figure 5. The central part is part of a cylinder and the two side parts are a part of a ball. (f) The chisel pick has a wear flat that is very small compared with the cutting depth but larger than the radius of the cross-section of the crushed zone, as shown in Figure 4. (g) The resultant force on the chisel pick, F, is the sum of the normal force, Fa and friction force Fb generated by the crushed zone, as shown in Figure 4, where Fb = Fa tan φ, and φ is the friction angle between the chisel pick and rock.
The CEIT model is different from Evans' tensile fracture model in the following two aspects: the CEIT model considers the sidewall effect and is a model for three-dimensional rock cutting, the force acting on the chisel pick is from a crushed zone rather than directly from the rock chip. The rock chip in the CEIT has a more complex shape, which makes it more difficult to obtain an analytical solution to the cutting force. In addition, the chisel pick has a wear flat in the CEIT model. This is actually closer to reality because the severe wear of the tip will produce a wear flat on the chisel pick even though the original chisel pick is ideally sharp.
indentation cycle is composed of three parts, as shown in Figure 5. The central part is part of a cylinder and the two side parts are a part of a ball. (f) The chisel pick has a wear flat that is very small compared with the cutting depth but larger than the radius of the cross-section of the crushed zone, as shown in Figure 4. (g) The resultant force on the chisel pick, F, is the sum of the normal force, Fa and friction force Fb generated by the crushed zone, as shown in Figure 4, where , and φ is the friction angle between the chisel pick and rock. The CEIT model is different from Evans' tensile fracture model in the following two aspects: the CEIT model considers the sidewall effect and is a model for three-dimensional rock cutting, the force acting on the chisel pick is from a crushed zone rather than directly from the rock chip. The rock chip in the CEIT has a more complex shape, which makes it more difficult to obtain an analytical solution to the cutting force. In addition, the chisel pick has a wear flat in the CEIT model. This is actually closer to reality because the severe wear of the tip will produce a wear flat on the chisel pick even though the original chisel pick is ideally sharp.

Mechanical Model and Mathematical Formula of the CEIT Model
The indentation force on an indenter is obtained by integrating the pressure acting on the contact surface between the cavity and the indenter [53,54]. This method is adopted to calculate the resultant force, F, acting on the chisel pick from the crushed zone. In addition, the friction between the chisel pick and rock particles in the crushed zone is considered. This friction will hinder the chisel pick from cutting forward. Thus, in the CEIT model, both the normal force, Fa, and friction force, Fb, will act on the chisel pick from the crushed zone, as shown in Figure 4. The radius of the cross-section of the crushed zone is very small compared with the cutting depth and cutting width. Thus, both ends of the crushed zone have negligible effects on the chisel pick compared with the middle part. Therefore, Fa and Fb can be estimated as: where p is the pressure in the crushed zone, w is the cutting width (or the width of the chisel pick), r is the radius of the crushed zone, φ is the friction angle between the chisel pick and rock.
where α is the rake angle, and sgn is the signum function, which is introduced here to change the direction of the friction force along the front face of the chisel pick. In Equations (16) and (17), p and r need to be determined. Essentially, we only need to know the product of p and r. The value of pr can be obtained by analyzing the fracture of the rock chip. As presented in Figure 4, the resultant force, P, acting on the rock chip from the crushed zone is: As mentioned above, we assume that at the time when the rock chip fails, the moment generated by all forces acting on the rock chip about axis B is equilibrium, similar to that described in Evans' model [20]. Axis B is defined as the axis through point B in Figure 4 and perpendicular to the side

Mechanical Model and Mathematical Formula of the CEIT Model
The indentation force on an indenter is obtained by integrating the pressure acting on the contact surface between the cavity and the indenter [53,54]. This method is adopted to calculate the resultant force, F, acting on the chisel pick from the crushed zone. In addition, the friction between the chisel pick and rock particles in the crushed zone is considered. This friction will hinder the chisel pick from cutting forward. Thus, in the CEIT model, both the normal force, Fa, and friction force, Fb, will act on the chisel pick from the crushed zone, as shown in Figure 4. The radius of the cross-section of the crushed zone is very small compared with the cutting depth and cutting width. Thus, both ends of the crushed zone have negligible effects on the chisel pick compared with the middle part. Therefore, Fa and Fb can be estimated as: where p is the pressure in the crushed zone, w is the cutting width (or the width of the chisel pick), r is the radius of the crushed zone, φ is the friction angle between the chisel pick and rock. The peak cutting force, Fc, and the peak normal force, Fn, on the chisel pick, can be estimated as follows: where α is the rake angle, and sgn is the signum function, which is introduced here to change the direction of the friction force along the front face of the chisel pick.
In Equations (16) and (17), p and r need to be determined. Essentially, we only need to know the product of p and r. The value of pr can be obtained by analyzing the fracture of the rock chip. As presented in Figure 4, the resultant force, P, acting on the rock chip from the crushed zone is: As mentioned above, we assume that at the time when the rock chip fails, the moment generated by all forces acting on the rock chip about axis B is equilibrium, similar to that described in Evans' model [20]. Axis B is defined as the axis through point B in Figure 4 and perpendicular to the side view of the rock chip in Figure 5. Based on the assumption that the tensile failure path has a horizontal tangent at point A in Figure 4, the radius of the circular arc, Rc, can be expressed as where d is the cutting depth and ψ is the breakout angle. Therefore, the moment about axis B generated by P, M 1 is calculated as follows: In addition, based on the assumption that the rock chip fails due to tensile fracture and only normal stress with a magnitude equal to the tensile strength of rock σ t acts on the fracture surface, the moment about axis B generated by the force T acting on the rock chip from the rock mass, M 2 can be calculated as follows: where S is the fracture surface of a rock chip, n is a unit normal vector of the fracture surface at the calculated point, r is the vector pointing to the calculated point from point B, and b is a unit vector of axis B. Based on the simplified shape of the rock chip, M 2 can be obtained where f (ψ) = [sin(4ψ) cos(2ψ) − 4ψ cos(2ψ) + π sin 3 (2ψ)]/(8 sin 4 ψ). As the rock chip is in moment equilibrium at the time it starts to fail, taking moment about point B, we can get where the effect of gravity is ignored because it is very small. Combining Equations (16), (17), (20), (22), and (23), the analytical solutions of Fc and Fn on the chisel pick are obtained.
The value of ψ can be obtained with the minimum work hypothesis, similar to that in [20,26], i.e., the value of ψ is obtained by solving the equation

An Approximate Value of the Breakout Angle in the CEIT Model
When the ratio of the cutting width to the cutting depth w/d tends to infinity (d/w tends to zero), Equation (26) can easily be solved, and the breakout angle ψ can be obtained by ψ = (π/2 + α)/4. Otherwise, it is difficult, even impossible, to find an analytical solution to the breakout angle due to the complexity of the expression of the peak cutting force Fc, Equation (24), and the peak normal force Fn, Equation (25). In this situation, an approximate solution to ψ is required to facilitate the application of the CEIT model. It is easy to figure out from Equations (24) and (26) that the value of ψ is only affected by the value of the rake angle α and the ratio of the cutting width to the cutting depth (w/d). By solving Equation (26)  of w/d. The relationships of k and b to w/d may be approximated by hyperbolas, as presented in Figure   6b. Hence, an approximate value of ψ can be obtained by using a regression formula as follows: where the unit of ψ and α is radian. The relationship between the values of ψ is obtained by directly solving Equation (26) and those calculated by the regression formula, Equation (27), are shown in Figure 7, where α is between -1 rad and 1 rad, and w/d is between 0.05 and 20. The results in Figure 7 indicate that the values of ψ calculated by Equation (27) agree well with that obtained by solving Equation (26). In addition, the ranges of α and w/d used during the above analysis are large enough for practical applications. When the value of w/d tends to infinity, Equation (27)     The results, presented in Figure 6a, indicate that there is a strong linear relationship between the breakout angle and the rake angle when w/d is a constant value between 0.05 and 20. The slopes (k) and intercepts (b) of the best fitting lines in Figure 6a between ψ and α depend only on the values of w/d. The relationships of k and b to w/d may be approximated by hyperbolas, as presented in Figure 6b. Hence, an approximate value of ψ can be obtained by using a regression formula as follows: where the unit of ψ and α is radian. The relationship between the values of ψ is obtained by directly solving Equation (26) and those calculated by the regression formula, Equation (27), are shown in Figure 7, where α is between −1 rad and 1 rad, and w/d is between 0.05 and 20. The results in Figure 7 indicate that the values of ψ calculated by Equation (27) agree well with that obtained by solving Equation (26). In addition, the ranges of α and w/d used during the above analysis are large enough for practical applications. When the value of w/d tends to infinity, Equation (27) tends to ψ = 0.246α + 0.404, which is very close to the analytical solution of Equation (26), ψ = (π/2 + α)/4. Thus, Equation (27) can be regarded as an approximate solution to Equation (26), and used for estimating the value of the breakout angle ψ in practical applications.   Substituting the value of the breakout angle estimated by Equation (27) into Equations (24) and (25), the peak cutting force and the peak normal force in rock cutting with a chisel pick can be obtained. The values of cutting width w, cutting depth d, rake angle α, tensile strength of rock σ t and friction angle between the chisel pick and rock φ are required for predicting the peak cutting forces in the CEIT model, similar to that in Evans' model.

Validation of the CEIT Model
As the CEIT model is established for rock cutting under the brittle fracture regime, very few experimental results are available and suitable for validating the CEIT model, even though numerous rock cutting experiments have been performed in the past few decades. In this paper, the rock cutting experiments performed in unrelieved conditions by Bilgin [28] and Allington [27] with chisel picks are used for validation. In addition, the CEIT model is compared with the existing theoretical models, including Evans' tensile breakage model and Nishimatsu's shear breakage model, to examine its performance in predicting the peak cutting forces on the chisel pick. Evans' model, Equation (6), Nishimatsu's model, Equation (12), the CEIT model, and Equations (24) and (27), are used for calculating the peak cutting force. As the CEIT model is an analytical model, the semi-empirical models and regression models, proposed based on a certain set of rock cutting experiments, are not used for comparison in this paper. The mechanical properties of the rock samples are presented in Table 2. The contact area between the chisel pick and rock has been assumed to be quite small in the CEIT model. Therefore, the friction angle between the tungsten carbide and rock, φ c , is used as φ in the CEIT model. While, in Nishimatsu's model, the friction angle between the hardened steel and rock, φ s , is used as φ because the contact area between the chisel pick and rock is larger. Since the friction angle between the chisel pick and rock samples were not tested in [27], the corresponding friction angle in [28] with the same rock type is used in the calculation, as presented in Table 2.

Validity of the CEIT Model in Capturing the Effects of the Cutting Parameters
In this part, the peak cutting force Fc at different values of the cutting parameters, including cutting width w, cutting depth d, and rake angle α, are predicted with the CEIT model. Then, the predicted Fc is compared with Fc calculated by Evans' model and Nishimatsu's model, and Fc measured from experiments. The values of Fc measured in [28] and calculated by the three models at different cutting widths are presented in Figure 8 It can be noticed that the best fitting lines of the experimental data in Figure 8 have positive intercepts for all the tested rock types. This indicates that when w tends to zero, Fc will not tend to zero, but tends to a positive value. This phenomenon is caused by the sidewall effect. The best-fitting lines for w and Fc predicted by the CEIT model also have positive intercepts, because the CEIT model has counted the sidewall effect during predicting Fc. From the phenomenological perspective, the side parts of the rock chip presented in Figure 5 causes the sidewall effect. From the mathematical perspective, the term ( ) is used for considering the sidewall effect. However, neither Evans' model nor Nishimatsu's model has considered the sidewall effect. Thus, Fc calculated by both of them is proportional to w and is usually smaller than Fc measured from experiments. The values of Fc measured in [27,28] and calculated by the three models at different cutting depths are presented in Figure 9, where the values of α and w remain unchanged for each rock type.
The results show that Fc, measured from the experiments in cutting anhydrite, limestone, and granite, increases linear with increasing d, while in cutting greywacke, anhydrite 2 and limestone 2, Fc It can be noticed that the best fitting lines of the experimental data in Figure 8 have positive intercepts for all the tested rock types. This indicates that when w tends to zero, Fc will not tend to zero, but tends to a positive value. This phenomenon is caused by the sidewall effect. The best-fitting lines for w and Fc predicted by the CEIT model also have positive intercepts, because the CEIT model has counted the sidewall effect during predicting Fc. From the phenomenological perspective, the side parts of the rock chip presented in Figure 5 causes the sidewall effect. From the mathematical perspective, the term f (ψ)d/w is used for considering the sidewall effect. However, neither Evans' model nor Nishimatsu's model has considered the sidewall effect. Thus, Fc calculated by both of them is proportional to w and is usually smaller than Fc measured from experiments.
The values of Fc measured in [27,28] and calculated by the three models at different cutting depths are presented in Figure 9, where the values of α and w remain unchanged for each rock type. The results show that Fc, measured from the experiments in cutting anhydrite, limestone, and granite, increases linear with increasing d, while in cutting greywacke, anhydrite 2 and limestone 2, Fc increases almost in a parabolic form with increasing d. These phenomena are perfectly represented by the CEIT model, as shown in Figure 9. However, Fc calculated by both Evans' model and Nishimatsu's model is proportional to d for all the tested rock types. It should be noticed that when the ratio of w to d is very small, the sidewall effect will influence the rock cutting greatly. Since the sidewall effect is proportional to the square of d [15,36,37], there will be a parabolic relationship between Fc and d. This is well captured by the CEIT model with Equation (24). In addition, the results in Figure 9 also indicate that the CEIT model is better than Evan's model and Nishimatsu's model in predicting Fc of the related experiments because the values of Fc predicted by the CEIT model are much closer to the measured values than those calculated by the other two theoretical models.
The values of the peak cutting force Fc measured in [28] and calculated by the three models at different rake angles are presented in Figure 10, where the values of the cutting width w and the cutting depth d remain unchanged for each rock type. The experimental results show that the values of Fc decrease with the increasing rake angle α in most of the cutting situations. Similar to Evans' model and Nishimatsu's model, the CEIT model can represent a decreasing trend of Fc with increasing α. As shown in Figure 10, the values of Fc calculated by the CEIT model show almost the same decrease trend as that of Fc measured from experiments with increasing α. It is obvious that the values of Fc, calculated by the CEIT model, are much closer to the measured Fc than those of Fc calculated by Evans' model and Nishimatsu's model.
In conclusion, the results in Figures 8-10 indicate that the CEIT model can well capture the relationships of Fc to w, d, and α, and the CEIT model can predict Fc on the chisel pick better than Evans' model and Nishimatsu's model.

Comparison of the CEIT Model and the Existing Theoretical Models
In this part, the CEIT model and the existing theoretical models, including Evans' model and Nishimatsu's model, are used to predict the values of the peak cutting force Fc on the chisel pick of the orthogonal experiments presented in [28]. The predicted results are compared with the experimental results to examine the validity of the CEIT model.
The comparisons of Fc calculated by formulas with that measured from rock cutting experiments are presented in Figures 11-14, respectively, for cutting anhydrite, limestone, granite, and greywacke. The best-fitting lines between the predicted Fc and measured Fc are also presented in these figures. Observing Figures 11-14, it is obvious that the values of Fc predicted by the CEIT model agree much better with the measured values than those predicted by Evans' model and Nishimatsu's model. The values of the peak cutting force Fc measured in [28] and calculated by the three models at different rake angles are presented in Figure 10  In conclusion, the results in Figures 8-10 indicate that the CEIT model can well capture the relationships of Fc to w, d, and α , and the CEIT model can predict Fc on the chisel pick better than Evans' model and Nishimatsu's model. Figure 10. Comparison of the calculated and measured peak cutting forces at different rake angles.

Comparison of the CEIT Model and the Existing Theoretical Models
In this part, the CEIT model and the existing theoretical models, including Evans' model and Nishimatsu's model, are used to predict the values of the peak cutting force Fc on the chisel pick of the orthogonal experiments presented in [28]. The predicted results are compared with the experimental results to examine the validity of the CEIT model.
The comparisons of Fc calculated by formulas with that measured from rock cutting experiments are presented in Figures 11-14, respectively, for cutting anhydrite, limestone, granite, and greywacke. The best-fitting lines between the predicted Fc and measured Fc are also presented in these figures. Observing Figures 11-14, it is obvious that the values of Fc predicted by the CEIT model agree much better with the measured values than those predicted by Evans' model and Nishimatsu's model. Figure 15 demonstrates the relationships between the measured Fc and the calculated Fc by the CEIT model, Evans' model, and Nishimatsu's model for all the tested rock types. To quantify the discrepancies between the predicted results and the measured results, root mean square error (RMSE), Equation (28), and bias, Equation (29), are introduced here. Figure 10. Comparison of the calculated and measured peak cutting forces at different rake angles. Figure 15 demonstrates the relationships between the measured Fc and the calculated Fc by the CEIT model, Evans' model, and Nishimatsu's model for all the tested rock types. To quantify the discrepancies between the predicted results and the measured results, root mean square error (RMSE), Equation (28), and bias, Equation (29), are introduced here.
where Fc m i is the i th measured Fc, Fc p i is the i th predicted Fc, and N is the number of the measured Fc. The RMSE provides a quantitative measure of the prediction model. A smaller value of the RMSE indicates the better performance of the prediction model. The bias indicates the correspondence between the mean values of predicted Fc and measured Fc. The closer to zero the bias is, the better the prediction model performs. The values of the RMSE, bias, and R 2 for each model in predicting the experimental results in [28] are listed in Table 3. The results in Table 3 indicate that the CEIT model is much better than Evans' model and Nishimatsu's model in predicting the values of Fc on the chisel pick in rock cutting. To conclude, the CEIT model is quite suitable for predicting the value of Fc in rock cutting with a chisel pick under the brittle fracture regime.

Sidewall Effect in the CEIT Model
One of the great differences between the CEIT model and the existing theoretical models is that the CEIT model considers the sidewall effect and is a three-dimensional model, as illustrated in Figure 5. The sidewall effect on the peak cutting force is counted by using the term f (ψ)d/w in Equation (24), where f (ψ) is a non-dimensional function of the breakout angle ψ, d is the cutting depth, w is the cutting width. When the ratio of w to d tends to infinity, as is the case in shallow cutting or two-dimensional rock cutting, the sidewall effect will tend to zero. In this situation, by assuming f (ψ)d/w = 0 and considering the minimum work hypothesis, Equation (24) can be simplified to the following equation: where the peak cutting force Fc is proportional to w and d, same with that in Evan's model and Nishimatsu's model. On the contrary, when the ratio of w to d tends to zero, Fc is totally caused by the sidewall effect. In this situation, Equation (24) can be rearranged into the following form: where Fc is proportional to the square of d and proportional to the tensile strength of rock σ t . Equation (31) indicates that the influence of the sidewall effect on Fc is proportional to the square of d, and is independent of w, which agrees with the findings in [36,37]. Equation (31) is quite similar to the theoretical model proposed by Evans [22] for predicting Fc on the conical pick. The prediction performance of Equation (31) for Fc on the conical pick is not discussed in this paper. In fact, since the influence of the sidewall effect on Fc is independent of w, it can be calculated by Equation (31) at any value of w/d. To quantify the sidewall effect, the sidewall effect coefficient (η) is introduced here. The value of η represents the proportion of the cutting force caused by the sidewall effect to the total cutting force. As presented in Equation (27), the value of ψ is only affected by the rake angle α and the ratio of w to d. In addition, f (ψ) only depends on ψ. Therefore, the value of η depends on the value of α and the ratio of w to d. For the CEIT model, η is calculated as follows: The relationship between η and w/d predicted by the CEIT model at different values of α is presented in Figure 16. The results show that the value of η decreases with the increasing value of w/d when α remains unchanged. The value of η tends to 100% when the value of w/d tends to zero. The value of η tends to zero when the value of w/d tends to infinity. This indicates that ignoring the sidewall effect has great influence on the models' predicting accuracy when the w/d ratio is very small. Thus, the CEIT model is more suitable for predicting the peak cutting force on the chisel pick than other existing theoretical models particularly when the value of w/d is very small. The value of w/d needed to obtain η = 50% decreases with the increasing value of α. In addition, the value of η will decrease with the increasing value of α and with w/d remaining unchanged. The validity of Equation (32) in calculating the value of η is not discussed in this paper due to the lack of suitable experimental data.

Advantages and Limitations of the CEIT Model
As the sidewall effect is considered in the CEIT model, the physical model adopted in this paper in deriving the formula for predicting the peak cutting force is more similar to the reality than the physical models adopted in Evans' model and Nishimatsu's model. Therefore, the CEIT model is a strict three-dimensional model. By counting the sidewall effect, the CEIT model can capture the parabolic relationship between the peak cutting and the cutting depth obtained from some rock cutting experiments. This relationship cannot be explained by existing theoretical models, including Evans' model and Nishimatsu's model. Additionally, the CEIT model can capture the relationships of the peak cutting force to the cutting width and the rake angle very well, same as the existing models. In fact, the validation of the CEIT model and comparisons of the CEIT model and other theoretical models show that the CEIT model can predict the peak cutting force of the discussed experiments much better than Evans' model and Nishimatsu's model. In addition, the CEIT model is convenient for practical application because only two mechanical parameters of rocks are needed to predict the peak cutting force, same with Evans' model. However, the available experimental results for verifying the CEIT model are inadequate, and further validations are still required through new experimental results and simulation results.
Based on the assumptions made in this paper, the CEIT model has some limitations during application. First, the CEIT model is intended for rock cutting under the brittle fracture regime. It is not suitable for rock cutting under the ductile fracture regime. Next, in the CEIT model, the symmetry plane of the chisel pick is vertical and is parallel to the cutting direction. The effects of the inclination and deflection of this symmetry plane are not considered. In addition, the size of the wear flat is assumed to be much smaller than the cutting depth and cutting width. For a badly worn chisel pick, the CEIT model needs improvement because the friction between the chisel pick and rock will increase with the increasing size of the wear flat [28,64]. Finally, the current form of the CEIT model is intended for unrelieved cutting conditions and is not suitable for relieved cutting conditions. In relieved cutting conditions, the influence of the neighboring picks is complicated, which depends on the ratio of cutting depth to picks spacing and rock properties. Therefore, many works are still required to improve the CEIT model to make it suitable for relieved cutting conditions.

Advantages and Limitations of the CEIT Model
As the sidewall effect is considered in the CEIT model, the physical model adopted in this paper in deriving the formula for predicting the peak cutting force is more similar to the reality than the physical models adopted in Evans' model and Nishimatsu's model. Therefore, the CEIT model is a strict three-dimensional model. By counting the sidewall effect, the CEIT model can capture the parabolic relationship between the peak cutting and the cutting depth obtained from some rock cutting experiments. This relationship cannot be explained by existing theoretical models, including Evans' model and Nishimatsu's model. Additionally, the CEIT model can capture the relationships of the peak cutting force to the cutting width and the rake angle very well, same as the existing models. In fact, the validation of the CEIT model and comparisons of the CEIT model and other theoretical models show that the CEIT model can predict the peak cutting force of the discussed experiments much better than Evans' model and Nishimatsu's model. In addition, the CEIT model is convenient for practical application because only two mechanical parameters of rocks are needed to predict the peak cutting force, same with Evans' model. However, the available experimental results for verifying the CEIT model are inadequate, and further validations are still required through new experimental results and simulation results.
Based on the assumptions made in this paper, the CEIT model has some limitations during application. First, the CEIT model is intended for rock cutting under the brittle fracture regime. It is not suitable for rock cutting under the ductile fracture regime. Next, in the CEIT model, the symmetry plane of the chisel pick is vertical and is parallel to the cutting direction. The effects of the inclination and deflection of this symmetry plane are not considered. In addition, the size of the wear flat is assumed to be much smaller than the cutting depth and cutting width. For a badly worn chisel pick, the CEIT model needs improvement because the friction between the chisel pick and rock will increase with the increasing size of the wear flat [28,64]. Finally, the current form of the CEIT model Figure 16. Relationship between sidewall effect coefficient and the ratio of cutting width to cutting depth at different rake angles.

Conclusions
This paper aims to deepen the understanding of rock cutting with the chisel pick of a cutter suction dredger during rock excavation. An analytical model, the CEIT model, was proposed for rock cutting with a chisel pick under the brittle fracture regime in order to predict the peak cutting force more accurately. The CEIT model is developed based on the tensile breakage theory and cavity expansion theory. Different from the existing theoretical models for a chisel pick, the sidewall effect is taken into account in this three-dimensional CEIT model. After the deduction of the mathematical formula of the CEIT model, the relationships of the breakout angle to the rake angle and the ratio of cutting width to cutting depth were investigated, and a simplified CEIT model was established for the calculation in practical applications.
Validation was made by comparing the predicted results with the available experimental data and the existing theoretical results. The results show that the CEIT model can well reproduce some phenomena observed from experimental studies on rock cutting with a chisel pick, which are as follows: (1) the peak cutting force increases linearly with the increasing cutting width, and decreases with the increasing rake angle; (2) the peak cutting force increases in a parabolic form with the increasing cutting depth when the ratio of cutting width to cutting depth is very small; (3) the peak cutting force is larger than zero when the cutting width tends to zero, which cannot been reproduced by existing theoretical models. In addition, the CEIT model can predict the peak cutting force on the chisel pick more accurately then the existing theoretical models, particularly for rock cutting at a small ratio of cutting width to cutting depth.
Finally, according to the CEIT model in this paper, it was found that the influence of the sidewall effect on the peak cutting force depends on the cutting depth and rake angle but is independent of the cutting width. To quantify this influence, the sidewall effect coefficient was defined, which decreases with the increasing cutting width to cutting depth ratio and the rake angle.