Analytical prediction of cutting forces in cylindrical turning of 304 stainless steel using unequal division shear zone theory

In order to establish a more accurate prediction model of turning forces, this paper proposed an analytical model for cylindrical turning with the consideration of the effect of the main cutting edge angle and the nose radius. Meanwhile, the unequal division shear zone theory in orthogonal free cutting is extended and applied to the oblique non-free cutting in the interaction between the chip units. To take into account the real tool nose geometry, the tool nose involved in the cutting is discretized into a series of cutting edge units. The geometrical parameters associated with the cutting edge units are analysed by using the coordinate transformation approach. Then, the improved oblique cutting model is applied to each cutting edge unit to acquire the component forces along the tool rake face. Finally, the resultant cutting forces in the turning process are calculated by the numerical integration method. To verify the effectiveness of the proposed model, the turning force experiment of 304 stainless steel was carried out by changing the cutting conditions, the main cutting edge angle, and the nose radius. Through the comparative analysis between the measured results and the calculation values of the proposed model, it was found that the analytical prediction of cutting force is in good agreement with the experiment.


Introduction
Through the accurate prediction of cutting force, the tool and machining parameters can be reasonably selected to improve cutting performance. Many researchers around the world have also done a lot of research on oblique cutting, because the majority of metal machining is non-free oblique cutting. At present, although most of the research is focused on the cutting force prediction of orthogonal free cutting, in which the cutting direction is performed along one straight cutting edge, non-free cutting, in which both primary and minor cutting edges are involved, is most widely used such as cylindrical turning in practice. The cutting conditions, the nose radius, and the main cutting edge angle also affect the cutting force in turning; therefore, comprehensive consideration of these factors is essential to accurately predict cutting forces. Besides, the dead metal zone also has major effect on cutting force and cutting performance in micromachining [1,2].
Cylindrical turning, one of the commonest processes in manufacturing components, is widely used in automotive, aerospace, and so on. It is a typical non-free cutting operation whose chip deformation process is very complex. The calculation of cutting forces in turning is roughly divided into three methods: empirical, analytical, and numerical models. Parakkal et al. [3] developed a relatively simple mechanistic model to predict cutting forces for grooved tool in oblique cutting, which considered the effect of the tool edge radius. This model was used to investigate the effects of cutting conditions and grooved parameters on the cutting forces. Kone et al. [4]. identified empirical cutting force equations; the chip flow angle of groove coated tool cutting 304L austenitic steel under dry cutting condition is deduced. Venkatarao [5,6] proposed a power consumption optimization-base strategy couple with 3D FEM simulation to reduce power consumption in micro-machining, the optimal value of process parameters is obtained for minimum cutting force, and a power consumption model was developed in terms of cutting force directions in ultrasonic vibration helical milling. Colwell [7] proposed the equivalent cutting edge that is used to replace the real cutting edge with an imaginary straight line in the tool rake face and assumed that the equivalent cutting edge is perpendicular to the chip flow direction. However, the chip flow direction defined above is not suited for cylindrical turning due to the nose radius. Hu et al. [8] extended Colwell's concept of equivalent cutting edge to the case of non-zero inclination angle. Kushima and Minato [9] divided the cutting edge into many small units and supposed that the component chip flow direction of each unit is perpendicular to the component cutting edge, and then the resultant chip flow direction is obtained by vectorial sum of the component chip flow directions. Young et al. [10] regarded the chip as many independent chip units and assumed that the thickness and direction of undeformed chip part of each unit are variable. The equivalent cutting edge is a straight line perpendicular to the direction of friction and can be obtained by calculating the resultant force per unit friction. Similarly, Wang et al. [11]. extended Young's method to the tool nose with non-zero rake and inclination angle and evaluated cutting force through numerical integration. Endres et al. [12]. introduced a nonlinear relationship among the unit friction, cutting force, and the chip thickness by considering the size effect of rounded edge on the tool nose. Khlifi et al. [13]. applied an equivalent tool geometry method to induce the same cutting force components as the real tool while taking into account the radius of the tool nose and edge. In addition to estimating the tool flank and crater wear, this method can also predict cutting force components in turning.
In the primary shear zone, Merchant [14] proposed a single shear plane analytical model. The shear angle was obtained by the minimum energy theory; many analytical models were then proposed for the modeling of orthogonal cutting process, such as the classical parallel-sided shear zone theory proposed by Oxley [15], which supposed that the primary shear zone is equally divided by the primary shear plane. However, Li et al. [16] established the unequal division shear zone theory by observing the strain rate distribution in the primary shear zone is inconsistent, the analytical results showed that the unequal division shear zone theory has acceptable accuracy. In subsequent research, this model is widely used to predict the cutting force and chip flow direction and to identify material constitutive equation. In addition, the existence of the rounded nose and the main cutting edge angle makes the cutting mechanism more complicated. Simultaneously, they have a significant influence on the geometry of the undeformed chip area, cutting forces, chip geometry, heat generation, and tool wear during turning. Redetzky et al. [17] divided the cutting edge into many trapezoidal units that are perpendicular to the component edge and calculated the cutting force of each discrete unit through unit cutting force, unit edge force, and unit lateral force. Hagiwara et al. [18] developed the hybrid model by Redetzky et al. [17] to predict the chip side-flow angle for complex grooved tool inserts in contour turning operations by applying the effective cutting conditions and tool geometry. Armarego et al. [19]. put forward to the generalized cutting edge concept to replace the rounded edge of tool nose. Storch et al. [20] investigated the unit force distribution on the rounded edge for free cutting and found that both the tool nose and the main cutting edge affect the machined surface quality and chip formation process. Fang [21] developed an improved oblique cutting model different from traditional model by considering the two cases where the tool main cutting edge angle is not equal to 90° and the influence of the tool feed velocity on the resultant cutting velocity is not ignored. Molinari et al. [22] established a new thermomechanical model of oblique cutting to take into account the real cutting edge geometry; the rounded nose is discretized into a set of cutting edge units to obtain the component cutting force and temperature distribution along the rake face, but their oblique model is based on the inertia item, mainly suitable for adiabatic shear process in high-speech machining. Budak et al. [23] proposed a thermo-mechanical model that considered both the sticking and sliding contact on rake face, and applied this model to obtained the local forces acting on each element by dividing the nose of the turning into small elements. Abdellaoui et al. [24] also analysed the effects of tool nose radius and established a thermomechanical model, which allows predicting cutting force components and thermo-mechanical parameters.
In conclusion, the existing turning force prediction model mainly apply the single shear plane theory and the parallel shear zone theory in the primary shear zone. In such case, these models cannot accurately describe the real plastic deformation state and stress-strain relationship of the primary shear zone. In this paper, a novel turning force prediction model based on unequal shear theory is presented for cylindrical turning, considering the influence of the tool nose radius and the main cutting edge angle on cutting force.
In this present work, the oblique cutting theory based on the unequal division shear zone theory investigated by Li et al. [16,25], and combined with the tool discretization approach of Molinari and Moufki [22], is developed to predict the cutting forces in turning of 304 stainless steel. This paper is organized as follows. Section 2 presents the analytical modeling of non-free oblique cutting with geometrical characterization and governing equations. The unequal shear theory and material constitutive equation are described in Sect. 3. In Sect. 4, the discretization of the tool cutting edge and the component geometry of each cutting edge unit are described. The component and resultant cutting force calculation are detailed in Sect. 5. The turning experiments for 1 3 304 stainless steel were conducted, and the effect of cutting parameters on cutting force is also discussed in Sect. 6. Finally, some conclusions are provided in Sect. 7.

Characteristics of the improved oblique cutting
To study the turning process, an analytical model of oblique cutting of the unequal division shear zone theory is developed for the case that the tool cutting edge angle relative to the main cutting edge is not equal to 90°, namely, the feed rate f r is not equal to the undeformed chip thickness t. In view of this, Fang [21] also proposed an improved method for the oblique cutting model and applied it to chip control experiment. Figure 1 illustrates the process of improved oblique cutting. In this model, where t, a w , and φ r are undeformed chip thickness, cutting width, and lead angle (complementary angle of the main cutting edge angle κ r ), respectively. The cutting depth a p , feed rate f r , cutting velocity V, normal rake angle γ n , and inclination angle λ s are known, and the geometric relations of these parameters are shown in Fig. 2. Some planes that can be used to determine the location of the cutting edges, rake face, and flank face are defined.
• The reference plane P r • The working plane P f • The back plane P p • The cutting plane P s • The orthogonal plane P o • The normal plane P n • The shear plane P sh • The equivalent plane P e As shown in Figs. 1 and 2, to calculate the relationships of geometry, motion, and force in the improved oblique cutting model, a series of coordinate systems are defined.
• Frame (x a , y a , z a ), x a is parallel to the cutting velocity direction, z a is parallel to the feed direction, (x a , y a ) is the back plane P p , (y a , z a ) is the reference plane P r , (x a , z a ) is the working plane P f .
The eight frames can be related by rotation transformation, as shown in Fig. 3. Take for example, (x o , y o , z o ) is obtained by the rotation of angle -φ r of the frame (x a , y a , z a ) around the x a axis. It can be expressed as the following form: The rest of the frames can also be given by the above method; Φ n is the normal shear angle, η c is the chip flow The improved oblique cutting process angle, η s is the shear flow angle, and η e is the equivalent plane angle. To locate the equivalent plane of each discrete edge unit, of which k is used to represent the number, the shear flow angle and the equivalent plane angle are taken by the following expression [22].
Since the chip flow angle is affected by the non-free chip in turning, the chip flow angle η c formula of free oblique cutting which was proposed by Moufki and Molinari [26] or others no longer suitable for no-free oblique cutting. The chip velocity can be calculated as: The normal shear angle can be calculated with the Merchants formula [14].
The mean friction angle β k can be defined by mean friction coefficient, and the mean friction coefficient may be a power function of the chip velocity, as shown in Fig. 4. The friction coefficient is given by: where f 0 is constant, f is mean friction coefficient, while p < 0. The coefficients f 0 and p can be identified by the experimental data in orthogonal cutting. Fig. 2 The geometrical relationship of oblique cutting Fig. 3 The rotation transformation relationships between the coordinate systems

Governing equations of the primary shear zone
Johnson-Cook model is used as flow stress model of workpiece material [27].
where k , k , ̇k , and T k are the shear stress, the shear strain, the shear strain rate, and the instantaneous temperature of the workpiece, respectively. ̇0 is reference plastic strain rate, T m is melting temperature of the workpiece material, and T r is room temperature or initial workpiece temperature. A, B, C, m, and n is the material constants which can be identified through split Hopkinson bar (SPHB), the material coefficients of Johnson-Cook model are listed in Table 1.
The unequal division shear zone theory is used to predict the turning forces, as shown in Fig. 5. In this model, the strain, strain rate, and the temperature of discrete chip are calculated as follows [25]: where ̇k m is the maximum strain rate, k k u is unequal division ratio, q is a parameter to characterize the non-uniform distribution of the tangential velocity in the primary shear zone, and h is the thickness of the primary shear zone; 0.025 mm is an acceptable approximation of h in reference by observing Shaw's research [28]. Besides, the boundary of the primary shear zone can be considered adiabatic so that the heat transfer equation can be expressed in the following form, and the chip temperature at the entry boundary is equal to the initial workpiece temperature.  Where ρ, c, and μ represent material density, heat capacity, and Taylor-Quinney coefficient, respectively. These values of 304 stainless steel are ρ=7900kg/m 3 , c=400J/ (kg·K) and μ=0.85, respectively. The material behavior and thermo-mechanical evolution in the primary shear zone are governed by Eqs. (6), (7), (8) and (11). Using the boundary conditions (12), the shear stress k s on the main shear plane can be calculated.

Cylindrical turning geometry and the cutting edge discretization
Considering the complexity of chip formation in turning, the majority of current research is based on the equivalent cutting edge method, i.e., a fictitious straight cutting edge is used to replace the real cutting edge, as shown in Fig. 6. When some studies need to know the local information about stress and temperature on the rake face, this method is very disadvantageous to the study of tool design and tool wear. Therefore, in order to obtain the component state of real cutting edge during chip formation, the cutting edge needs to be discretized. In cylindrical turning, the modeling process is more complicated than the traditional oblique cutting with a single cutting edge on account for the existence of the tool rounded (12) T| z k e =0 = T w nose. In this paper, the discretization method of the cutting edge proposed by Molinnari and Moufki [22] is used to model the cylindrical turning process. The rounded edge is divided into finite straight cutting edge units. For each component edge unit, the expressions of component normal rake angle and inclination angle are derived by using geometric relationships. For the case of cylindrical turning, as shown in Fig. 7. The cutting conditions to be given are cutting velocity V, feed rate f r , cutting depth a p , and tool geometric parameters (the main cutting edge angle κ r , inclination angle λ s , normal rake angle γ n and nose radius r n ); the tool geometry is defined by six planes (P r , P s , P n , P o , P f , and P p ) to determine the spatial position of the cutting edge. Figure 8 illustrates the discrete process of tool nose in cylindrical turning.
On the reference plane P r , the cutting edge can be discretized into K + 2 parts of the straight cutting edge section by the index k ( 0 ≤ k ≤ K + 1 ). The tool nose is discretized into K + 1 straight cutting edge units, and the chip thickness corresponding to each component edge unit is variable as shown in Fig. 8. In the turning process, the rounded edge in contact with the workpiece can be divided into two parts by the angles ϕ s and ϕ a respectively. ϕ st and ϕ ex are the entry and exit angles in the first part (see Fig. 8). The first part (for the interval angle of ϕ s ) can be subdivided into K cutting edge units which are identified as k (1 ≤ k ≤ K), and Δϕ is taken as the increment angle. The second part is also considered as a single unit with the index k = K + 1. Figure 8 shows an Fig. 6 The equivalent cutting edge method example of discretization for K = 3, which corresponds to the projection of a discrete cutting edge on the reference plane. According to the geometric relationship, the following equations of these angles can be obtained: For the rounded edge, the geometric description of the discrete cutting edge is mainly defined by the planes P r ,P 0 s ,P 0 n ,P k s , and P k n , as shown in Fig. 8. P 0 s and P 0 n correspond to the cutting plane and normal plane of the main cutting edge (k = 0). P k s is defined as containing the kth component edge and being perpendicular to P r . P k n is normal to the kth component edge. k r is the angle between the planes P 0 s and P k s measured in the plane P r , k a is the angle between the main cutting edge (k = 0) and the kth component edge (1 ≤ k ≤ K) measured in the tool rake face (see Fig. 8). As shown in Fig. 9, the component parameters k , k r , k r , k r ,a k w , and t k in the reference plane P r can be obtained as below: where k , k r , k r , a k w , and t k represent, respectively, the angular position, the component cutting edge angle, the component lead angle, the component cutting width, and the component cutting thickness corresponding to the kth component cutting edge. The following relations can be obtained for the main cutting edge (k = 0): The chip load A k of the kth (0 ≤ k ≤ K) component edge is defined as the undeformed chip area. For the chip load A K+1 of the k = K + 1 component edge, the following relationship can be obtained.
The relationship between the discrete cutting edge on the rake face and its projection on P r is processed by coordinate transformation. The component cutting conditions corresponding to each component cutting edge, that is, the component cutting angles (cutting edge angle, inclination angle, normal rake angle) and the component undeformed chip area, Fig. 9 Projection of discrete cutting edge on reference plane are also determined. In the rake face. The projection of y k r on the reference plane P r is oriented by the unit vector y k o ; y k r is the vector along the kth component edge and is contained in the P k s as shown in Fig. 9. To calculate the angles k a , k s , and k n , a set of frames is introduced. These frames can be related by rotation transformation, as shown in Fig. 10. Based on the relationship of the coordinate systems shown in Fig. 10, the following expression can be obtained: Therefore, by considering Eq. (27) and after some algebraic manipulation, k r corresponding to the projection of k a in the reference plane P r is obtained: The component inclination angle k s and the component normal rake angle k n for the kth component edge are measured in P k s and P k n , respectively. For the main cutting edge (k = 0), the cutting angles are given by 0 s = s and 0 n = n . The angles k s and k n , for 1 ≤ k ≤ K + 1, are determined from the following relationships by Eq. (27).
The component cutting edge geometry can be characterized by these angles k r , k n , k s , k r , and k a .

Component equilibrium of the discrete chip unit
For the improved oblique cutting model, the shear force F k s and normal force F k ns exerted by the workpiece on the kth chip unit can be expressed as:  (x k , y k , z k ). The following expressions were deduced by Moufki and Molinari [26]. where k s is the shear stress on the main shear plane of the primary shear zone; it can be obtained from the unequal division shear zone theory.

Resultant chip equilibrium and cutting forces.
In order to determine the value of the main chip flow angle In the resultant chip balance, the sum of F k mc is equal to zero for the interaction force F k mc between discrete chip units.
Because the variables k e , k s , k n , k , and F k mc all depend on the main chip flow angle 0 c , numerical iteration is applied to  Eq. (36), it considers the influence of material properties of workpiece, tool geometry, and cutting parameters. Adibi-Sedeh et al. [29] considered the effects of the nose-rounded radius and the cutting angle on the chip flow direction when the main cutting edge angle is not equal to zero for the traditionally defined equivalent cutting edge.
For the improved oblique cutting model of the kth component edge, in Fig. 3, the resultant force F k is expressed as These force components can then be transformed to the frame (x a , y a , z a ) by: where F k x is the tangential force, F k y is radial force, and F k z is the longitudinal force as shown in Fig. 2. The algorithm to calculate cutting forces for cylindrical turning is as follows in Fig. 11. At the beginning, the component cutting angles can be calculated by giving an initial value of 0 c , These output parameters such as strain, strain rate, stress, and temperature distributions in the primary shear zone for each chip unit can be computed by numerical iterative method. If K+1 ∑ k=0 F k mc is not equal to zero, a new estimated value of 0 c will be obtained (39) Flow chart for the improved oblique cutting model applied to cylindrical turning by adding an increment Δ . The computation will stop until the iteration requirement is satisfied. Conclusively, the component cutting force ( F k x , F k y , F k z ) of the kth cutting edge unit is deduced by Eqs. (38) and (39); F mc can be determined by Eq. (36). All steps of calculations were performed programmatically by MATLAB software.

Experimental validation and analysis
In order to validate the proposed method above, a series of turning experiments for 304 stainless steel were conducted without cutting fluid in Mazak lathe. The cutting forces were measured by Kisler Dynamometer (6275A); the force signals that pass through the charge amplifier (5070A) was addressed by the Dynoware software, as shown in Fig. 12.
The experiment was carried out in three cases: • The influence of cutting depth a p , feed rate f r , and cutting velocity V on turning forces was analysed with different cutting parameters. • The influence of the tool nose radius r n on turning forces was analysed with different tool insert. • The influence of the main cutting edge angle κ r on turning forces was analysed with different toolholder.

Effects of cutting conditions
The experiment used SNMG120408-MR2025 insert and DSBNR2525M12 holder of Sandvik to study the influence of cutting parameters on turning forces of 304 stainless steel. The geometry parameters of the insert and the cutting conditions are shown in Tables 2 and 3, respectively. Specifically, there were three procedures in the experiment to obtain the reliable results, namely, the experiment with the same cutting parameters was repeated twice; the measurement interval time of the cutting process is more than 10 s, and the insert was replaced by a new one after each cutting test with the same cutting parameters.
In Fig. 13, it can be seen that the cutting force increases approximately linearly with the increase of cutting depth. When the cutting depth is very small, the change of cutting force shows obvious nonlinearity. The reason is that the straight main cutting edge plays a leading role in the cutting process when the cutting width a w (a w = a p /sin(κ r )) increases in a positive proportion with the cutting depth increases. Furthermore, when the cutting depth is close to the nose radius, the influence of the nose radius on the cutting force must be considered. At this time, the nonlinear variation can only be measured in a more stable environment with precise equipment.
As shown in Fig. 14, the cutting force also shows a linear trend with increasing feed rate, which is different from the effect of cutting depth on cutting force. This is because the chip thickness t (t = f r sin(κ r )) has been increasing in a positive proportion with the increase of feed rate. However, with the increase of feed speed, the increase of temperature will lead to the decrease of friction coefficient and cutting force in the turning process. Considering the positive and negative effects, the influence of feed rate on cutting force should be nonlinear. At the same time, it can be found from Figs. 14a  and b that the cutting force increase is slow with the increase in feed rate when the cutting velocity is very high. For the 304 stainless steels, the nonlinear effect of cutting velocity on turning force is shown in Fig. 15. It can be seen that the cutting force decreases with the increase of cutting velocity in Fig. 15; the variation trend of cutting forces is consistent with the experimental results of Young et al. [10]. and Arsecularate et al. [30]. This is due to the fact that the deformation coefficient and friction coefficient decrease with the cutting velocity increases. It should be noted that the tool-chip friction is mainly caused by high temperature, high pressure, and high chip velocity. In cylindrical turning, the friction coefficient of each discrete chip unit will change, the average tool-chip temperature decreases with the cutting velocity increases, which will lead to the decrease of average tool-chip friction coefficient and cutting force. Meanwhile, it can be seen in Fig. 15 that the changing ratio of cutting force is obviously slow when cutting velocity is very low.

Effects of the main cutting edge angle
The cutting experiments of 304 stainless steel with three kinds of the main cutting edge angle were carried out. Table 4 illustrates the geometric parameters of the tool, and normal rake angle γ n is approximately regarded as 0°, Fig. 13 The effect of cutting depth a p on turning forces  inclination angle λ s is − 6°, and nose radius r n is 0.8 mm.
The predicted values of the proposed model and the experimental results are shown in Fig. 16; the main cutting edge angle κ r is investigated in the range of 30° to 90°. It can be found from Fig. 16 that the influence of the main cutting edge angle is approximately linear. As the main cutting edge angle increases, F x and F z are increasing, F y is decreasing, and this trend is consistent with the experimental results of Young et al. [10]. However, F x basically remains constant with the increase of the main cutting edge angle. Although the incensement of the main cutting edge angle will result in a decrease in the engaged length in cutting, the area A (A = f r a p ) of the undeformed chip which mainly determines the cutting force not depend on the tool main cutting edge angle. It can be seen in Fig. 16 that F x is increased significantly when the main cutting edge angle increases to the range between 60 and 90°. It can generally be considered that the main cutting edge angle mainly affects the cutting force through the cutting thickness. With the increase of κ r , the contact length of the tool rounded edge will also increase (at this time, the non-free cutting effect will be more significant). When the cutting area is constant, the cutting thickness increases with the increase of the main cutting edge angle. As a result, the deformation of chip formation will become large, which will lead to the increase of F x . Nevertheless, the cutting thicknesses of the discrete chip units for tool nose part decreases in order as shown in Fig. 16, and they are smaller than the cutting thickness of the main cutting edge. When the cutting area is constant and the main cutting edge angle is increased, the cutting thickness will increase. The influence of the main cutting edge angle on F z and F y is mainly through the projection components of the resultant cutting force in the Z and Y directions. Therefore, F z will decrease and F y will increase as the main cutting edge angle increases.

Effects of the tool nose radius
The effect of four kinds of the tool nose radius on turning force of 304 stainless steel was studied. The rounded edge parameters are presented in Table 5. Figure 17 illustrates the predicted variation of turning force with the nose radius and the experimental results. In the experiments, the main cutting edge is not involved in cutting process (a p < r n (1-cos(κ r )). It can be seen in Fig. 17a that F x is basically a constant, which may be because the total undeformed chip area A is nearly independent of the nose radius and only related to feed rate and cutting depth. The predicted curve in Fig. 17 shows an increase of F y and a decrease of F z with the increase of the nose radius. The ig. 16 The effect of the main cutting edge angle κ r on turning force experimental results of 304 stainless steel are also basically consistent with the predicted law, which is also consistent with the experimental results of Young [10] and Arsecularatne et al. [30]. According to the previous discussion, F z will increase with the increase of the main cutting edge angle and decrease with the increase of 0 c ; the change law of F y is opposite. However, when the main cutting edge is involved in cutting (a p > r n (1cos(κ r )), as shown in Fig. 17b and c. When the cutting depth a p and feed rate f r are constant, the tool nose radius length and cutting width increase with the increase of the main cutting edge angle. Conversely, the cutting thickness and the tool component cutting edge angle of the discrete edge units will decrease; these lead to the increase of cutting deformation and cutting force. It has a greater influence on F y than on F x as shown in Fig. 17.

Conclusion
This paper extended the previous work for the traditional oblique cutting model to non-free turning with a rounded edge. The main conclusions can be summarized as follows: • The improved oblique cutting model considers the condition that the main cutting edge angle is unequal to be 90° and the tool nose has a rounded edge. • Through applying the coordinate transformation method, the geometry relations for the non-free turning are analysed. The unequal division shear zone theory in orthogonal cutting is used in oblique cutting modeling. • The tool rounded edge was discretized into small oblique cutting edges characterized by component cutting angles.
For every engaged cutting chip unit, the cutting forces in different directions are obtained during oblique cutting. The effects of tool geometrical parameters and workpiece material properties on the cutting force are considered as their coupling interaction.
A MATLAB-based program has been given to facilitate the application of the analytical model for cutting forces. The experiments with three types of cutting conditions have been conducted to validate the prediction model and the deviation between the predicted results and experiment data is within an acceptable range. The proposed models have also been used to investigate the effects of cutting conditions and tool angles on the cutting forces. This quantitative information is helpful in selecting the appropriate cutting parameters and optimization of cutting edge geometry. It should be noted that the above description law of the influence of the machining parameters on the turning forces has appeared in many literatures, but most are obtained by curve fitting from the experimental results. The significance of  Fig. 17 The effect of tool nose radius r n on turning force