Modelling and Optimization of Machined Surface Topography in Ball-End Milling Process

Abstract

In order to optimize machined surface topography, this paper presents a novel algorithm for simulating the surface topography and predicting the surface roughness of a ball-end milling process. First, a discrete workpiece model was developed using the Z-map method, and the swept surface of a cutter edge was represented using triangular approximation. The workpiece surface was updated (i.e., material removal process) using the intersection between the vertical reference line and the triangular facet under a cutting judgement. Second, the proposed algorithm was verified by comparing the simulated 3D surface topography as well as 2D surface profile and average roughness (Sa) with experimental measurements. Then, numerical simulation examples planed by the Box–Behnken design methods were carried out to investigate the Sa in the ball-end milling operation. The correlations of S_a and cutting parameters were represented by a response surface reduced quadratic model based on the ANOVA results. Finally, the feed per tooth, radial depth of cut, and tilt and lead angles were optimized for improving the machining efficiency under the Sa constraints. This study presents an effective method for simulating surface topography and predicting the Sa to optimize the cutting parameters during ball-end milling process.

1. Introduction

Surface topography has a clear effect on the performance of components, such as their contact, friction, wear, and lubrication properties [1,2]. Surface topography is characterized by the surface roughness quantity for evaluating surface quality. Ball-end milling is one of the most widely used machining operations in various industries, such as aerospace, automobile, bioengineering and molding, as a semi-finish and finish method to replace grinding or polishing for improving the machining efficiency [3]. Therefore, the surface finish and dimensional accuracy of the final part are mainly determined by the ball-end milling operation in the finishing stage. According, the surface topography is easily influenced by many factors (cutting parameters) via ball-end milling. Chen et al. [4] investigated the effects of cutter inclination angle on the surface topography geometrical features using cutting experiments. Zhang et al. [5] investigated the ratio and product of the feed per tooth and radial depth of cut using an analytical model for reducing the surface roughness at a constant material removal rate. Peng et al. [6] investigated the effects of the cutter initial phase angle on the surface topography in micro-ball-end milling, and the surface texture was controlled by planning the noncutting tool path. Zhang et al. [7] used a modified cutting edge to develop a surface topography model considering the tool flank wear effect. Zhang et al. [8] considered stochastic tool wear for simulating the surface topography in micro-milling. Lotfi et al. [9] developed a modified cutter/workpiece engagement algorithm for simulating the surface model in five-axis ball-end milling considering the cutter deflection and cutter runout. Hao et al. [10] introduced a correction model of ball-end milling by wear and force deformation for simulating the surface topography. Omar et al. [11] proposed an integral model for simulating the surface topography and predicting the cutting force simultaneously, and investigated the effects of cutter axis tilt, cutter runout, tool deflection, systems dynamics, and face wear for the surface topography in the ball-end milling. Wang et al. [12] investigated the effects of cutting vibration on the surface topography in ball-end milling of thin-walled parts using a numerical method. Li et al. [13] investigate the geometrical errors imposed by the machine tool on the surface topography in five-axis ball-end milling processes, and evaluated quantity by means of modal coefficients. Biondani et al. [14] developed a surface topography model for investigating the effect of cutter edge micro-geometry on the surface topography and roughness based on an actual cutting-edge model and a judgment of minimum chip thickness. Therefore, numerical simulation has more applications compared with milling experiments for investigating the surface topography in ball-end milling operation. The surface topography simulation involves lower costs and is less time-consuming.

The surface topography simulation methods can be classed into three groups: solid modeling, iterative methods, and discretization. In the solid modeling methods [15], the material removal process is simulated by the Boolean operations between the cutter swept volume and workpiece. The iterative method was first proposed by Zhang et al. [16]. Mathematically, it is like solving the intersection of the swept surface of the cutter edge and the reference line at a selected workpiece point by solving the nonlinear equation using iterative methods. In the discretization methods, first, the cutter is dispersed into a series of units. The workpiece is constructed by means of the Z-map, N-buffer, Dexel and Voxel elements, cut plane, points set, etc. Then, machining time (or tool path) is dispersed, and at each discrete time point (or discrete cutter location point), the surface topography value (i.e., coordinate of workpiece surface) can be determined by the cutter edge unit intersecting with the workpiece. According to related research works, these are the main problems with the solid modeling as well as iterative methods and discretization. Indeed, in the solid modeling techniques, the computation of the intersection between geometrical objects is difficult and time-consuming. In iterative methods, a good initial value (also called starting point) must be determined first, while the determination of the initial value is not an easy task due to the nonlinearity of the equation of the swept surface and multiple intersections with the reference lines. In discretization methods, the computation efficiency and accuracy are low because of the plentiful discretization for the cutter and workpiece. However, the discretization method is still widely employed for modeling as a result of its extensibility and robustness with respect to solid modeling and iterative methods. Moreover, some study works have been performed to improve the accuracy and efficiency of the discretization method. Li et al. [17] introduced iterative methods to improve the discrete methods. In other words, they proposed a new initial value selection algorithm using discretization methods for the iterative methods. Xu et al. [18], Yang et al. [19], and Chen et al. [20] improved the discrete method using different interpolation methods. Therefore, the discretization method was adopted to develop the surface topography model. Meanwhile, an efficient algorithm for calculating the intersection point of the swept surface and workpiece was proposed, which is the key technology for simulating surface topography successfully.

The optimization of cutting parameters is simple and effective in machining process for obtaining the desired surface topography and roughness. Dikshit et al. [21,22] optimized the cutting parameters for minimizing the surface roughness using the response surface methodology in the ball-end milling of AL2014-T6. Zhou et al. [23] developed a modified grey relational analysis model for optimizing the surface roughness and residual stress in the ball-end milling of Inconel 718. Buj-Corral et al. [24] developed an artificial neural network for correlating surface roughness with cutting parameters to obtain the optimal machining strategy in the five-axis ball-end milling of W-Nr. 12344. Venkata-Rao et al. [25] considered the surface roughness as a constraint for reducing power consumption in the micro-ball-end milling of D2 steel. Sonawane et al. [26] used a response surface methodology for optimizing the cutting parameters and obtaining the minimum surface roughness in the high-speed ball-end milling of Inconel-718 thin cantilevers. Torres et al. [27] proposed a structured numerical simulation using several friction coefficient approximations, and adapted it to a pre-textured martensitic precipitation hardening stainless-steel. Cross et al. [28] investigated the surface topography of material measures using confocal microscope and atomic force microscope measurements, and it was found that the surface topography depends primarily on the spindle speed and much less on the feed per tooth, whereas differences between the milling tools are observed in the tendency to form substructures. Klauer et al. [29] found that the ploughing of the material instead of a cutting occurs when too-small tilt angles are applied, and investigated the tilt angle’s influence on the resulting roughness and machined geometry using micro-milling sinusoidal freeform surfaces at different tilt angles. Arruda et al. [30] optimized the surface roughness using the Taguchi method in the ball-end milling of AISI P20. Despite all this, there is a lack of research on improving the machining efficiency when surface roughness is applied as a constraint.

The objective of this study is to present a surface topography simulation algorithm for replacing a milling experiment to study the surface roughness during ball-end milling processes. As mentioned earlier, in order to model the surface topography, the intersection point between the swept surface of the cutter edge and the workpiece must be solved first. Hence, in the current study, an efficient numerical algorithm was proposed for calculating the intersection point between the swept surface and workpiece using an improved discretization method. To this end, first, the swept surface of the cutter edge was formulated by means of the homogeneous coordinate transformation and the kinematics of the ball-end milling cutter. Then, a discrete workpiece model was developed using the Z-map method, and the swept surface was represented using a triangular approximation. Finally, the intersection is calculated from the vertical reference lines and triangular facet. Once the intersection is solved, the mesh point of the workpiece surface can be updated by the intersection point with a smaller z-coordinate. The main contribution of the current study pertains to the central part of the proposed model, i.e., an algorithm for calculating the intersection point of the swept surface and the workpiece for updating the workpiece surface.

Henceforth, the paper is organized as follows. In Section 2, the surface topography simulation algorithm is presented in detail. In Section 3, milling experiments are carried out to valid the proposed algorithm. A response surface model has been developed for improving the machining efficiency under the average roughness (Sa) constraint. Finally, the conclusions of the work are summarized in Section 4.

2. Modeling of Machined Surface Topography

In order to obtain the equation of the swept surface of a cutter edge, the coordinate systems involved in the ball-end milling of flat surface processes are established as shown in Figure 1.

Machine tool coordinate system, O_m − X_mY_mZ_m = {O_m; e₁^M, e₂^M, e₃^M}: the global coordinate system, O_m, is placed at the origin of the machine tool, and e₁^M, e₂^M and e₃^M can be determined by the kinematic structure of the machine tool.

Cutter coordinate system, O_c − X_cY_cZ_c = {O_c; e₁^C, e₂^C, e₃^C}: the local coordinate system is fixed on the ball-end milling cutter. O_c is placed at the ball center and defined as the cutter location (CL) point. e₃^C is along the cutter axis. e₁^C is along the tangent of the cutter edge at the cutter tip. e₂^C is the cross product of e₃^C and the e₁^C, e₂^C = e₃^C × e₁^C. The cutter revolves around the spindle axis at an angular speed ω, and the cutter axis and the spindle axis show eccentricity, e.

Spindle coordinate system, O_s − X_sY_sZ_s = {O_s; e₁^S, e₂^S, e₃^S}: local coordinate system attached to the spindle. O_s is placed at the intersection between the axis of the spindle and the X_cY_c plane. e₃^S is placed along the axis of the spindle, which is parallel to the e₃^C. The spindle translates along the cutter path at r_Oc = (x_Oc, y_Oc, z_Oc). When the tilt and lead angles are equal to zero, e₁^S is parallel to e₁^M and e₂^S is parallel to e₂^M, while e₃^S is parallel to e₃^M.

Process coordinate system, O_p − X_pY_pZ_p = {O_P; e₁^P, e₂^P, e₃^P}: the local coordinate system is placed at the cutter path. O_p is placed at the cutter path and consistent with O_s. e₁^P is along the feed direction and assumed to be e₂^W, e₃^P is along the normal direction of the cutter center location surface and is parallel to e₃^M, e₂^P is along the cross-feed direction and is also the cross product of e₃^P and e₁^P, e₂^P = e₃^P × e₁^P.

Workpiece coordinate system, O_w − X_wY_wZ_w = {O_w; e₁^W, e₂^W, e₃^W}: the global coordinate system is attached to the workpiece, in which the cutter path and surface topography are described. O_w is placed at an arbitrarily selected point on the workpiece and the directions of e₁^W, e₂^W, and e₃^W are assigned to be consistent with the O_m − X_mY_mZ_me₁^M, e₂^M and e₃^M, respectively.

The transformation of a selected point from O_c − X_cY_cZ_c to O_s − X_sY_sZ_s (see Figure 2) can be performed by a rotation around e₃^C and a translation along the O_sO_c. The transformation matrices M^CS can be expressed as

(1)${\mathit{M}}^{\mathrm{CS}}=\left[\begin{array}{cccc}\cos \left(-\omega \cdot t+{\mu }_0\right)& -\sin \left(-\omega \cdot t+{\mu }_0\right)& 0& e\cdot \cos \left(-\omega \cdot t+{\upsilon }_0\right)\\ \sin \left(-\omega \cdot t+{\mu }_0\right)& \cos \left(-\omega \cdot t+{\mu }_0\right)& 0& e\cdot \sin \left(-\omega \cdot t+{\upsilon }_0\right)\\ 0& 0& 1& 0\\ 0& 0& 0& 1\end{array}\right],$

where ω = 2πn/60 (rad/s) is the angular speed of the cutter and n (r/min) is the rotational speed of the spindle, t is the machining time, eccentricity e is the distance from O_c and O_s, μ₀ is the initial angle of e₁^C and e₁^S, and υ₀ is the initial phase angle of the cutter runout.

For a type-AB five-axis machining center, we must keep the cutter location point r_Oc in the proper position and the cutter axis in proper orientation. The spindle translates along the cutter path, and rotates around e₁^M at β₂, and rotates around e₂^M at β₁ relative to the body of the machine tool (see Figure 3). Therefore, the transformation matrix M^SW from O_s − X_sY_sZ_s to O_w − X_wY_wZ_w is defined as follows:

(2)${\mathit{M}}^{\mathrm{SW}}=\left[\begin{array}{cccc}\cos {\beta }_1& \sin {\beta }_1\cdot \sin {\beta }_2& -\sin {\beta }_1\cdot \cos {\beta }_2& x_{Oc}\\ 0& \cos {\beta }_2& \sin {\beta }_2& y_{Oc}\\ \sin {\beta }_1& -\cos {\beta }_1\cdot \sin {\beta }_2& \cos {\beta }_1\cdot \cos {\beta }_2& z_{Oc}\\ 0& 0& 0& 1\end{array}\right],$

where (x_Oc, y_Oc, z_Oc) is the coordinate of in O_w − X_wY_wZ_w, β₁ is the tilt angle, and β₂ is the lead angle.

As shown in Figure 4, with the parameter of the axial position angle θ of a selected point p along the cutter edge, in a cutter coordinate system, O_c − X_cY_cZ_c, the spherical cutter edge can be expressed as:

(3)$\begin{array}{cc}\left[\begin{array}{c}{x_{\mathrm{p}}}^{\mathrm{C}}\\ {y_{\mathrm{p}}}^{\mathrm{C}}\\ {z_{\mathrm{p}}}^{\mathrm{C}}\end{array}\right]=\left[\begin{array}{c}{R}_0\cdot \cos \theta \cdot \cos \phi \\ R_0\cdot \cos \theta \cdot \sin \phi \\ -R_0\cdot \cos \theta \end{array}\right],& \left(0\le \theta \le \tan {\beta }_0\right)\end{array},$

where (x_p^C, y_p^C, z_p^C) is the coordinate of the point P in O_c − X_cY_cZ_c. R₀ is the radius of the cutter. For a spherical cutter edge with equal lead, the lag angle φ can be calculated as follows:

(4)$\phi =\tan {\beta }_0\left(1-\cos \theta \right)+\left(j-1\right)2\mathsf{\pi }/N_{\mathrm{c}},$

where β₀ is the nominal helix angle as measured in the connection of the spherical and cylindrical cutter edges, j is the index of the cutter edge, and N_c is the number of cutter edges.

Given the cutter geometry as well as the cutter position and the cutter orientation, the swept surface of the cutter edge can be calculated in terms of homogeneous coordinate transformation. Mathematically, it is expressed as follows:

(5)$\left[\begin{array}{c}{x_{\mathrm{p}}}^{\mathrm{W}}\\ {y_{\mathrm{p}}}^{\mathrm{W}}\\ {z_{\mathrm{p}}}^{\mathrm{W}}\\ 1\end{array}\right]={\mathit{M}}^{\mathrm{SW}}{\mathit{M}}^{\mathrm{CS}}\left[\begin{array}{c}{x_{\mathrm{p}}}^{\mathrm{C}}\\ {y_{\mathrm{p}}}^{\mathrm{C}}\\ {z_{\mathrm{p}}}^{\mathrm{C}}\\ 1\end{array}\right],$

where (x_p^W, y_p^W, z_p^W) is the coordinate of a selected point P on the swept surface in O_w-X_wY_wZ_w.

The workpiece surface is updated (i.e., material removal process) based on a ray-triangle algorithm. First, the workpiece surface is dispersed into a series of discrete points using Z-map methods, as shown in Figure 5. Second, the swept surface of the cutter edge is dispersed into a series of triangle facets, as shown in Figure 6. Then, we update the discrete points of the workpiece surface using the intersection point between the vertical line of the workpiece surface and the triangular facet of the swept surface. Finally, we draw the surface topography as well as the surface profile, and calculate the surface roughness according to the coordinate of the discrete point of the workpiece surface. A sketched map of the surface topography simulation algorithm is shown in Figure 7.

3. Results and Discussion

3.1. Experimental Verification

In this study, the typical hot work die steel, AISI P20/3Cr2Mo (33 ± 3 HRC), combines a very good machinability with mirror polishing performance, and covers a wide variety of applications; it was chosen as the workpiece material here. The nominal chemical compositions and mechanical properties of the workpiece material are listed in Table 1 and Table 2, respectively. The workpiece samples were prepared as rectangular blocks (250 mm × 150 mm × 40 mm). Four samples were prepared. The top surfaces of each workpiece with sizes of 250 mm × 150 mm were divided into eight areas. Before milling experiments, the blocks were face-milled on the top and bottom surfaces to remove the heat treatment-related surface defects and flatness to eliminate errors in the experimental results. The cutter utilized in milling experiments was an integral carbide ball-end mill with two teeth (JH970100-Tribon, Seco Tools Company^®, Troy, MI, USA) with a diameter of 10 mm and a 30° helix angle.

The flat surface milling tests were carried out on a five-axis computer numerical controlled vertical machining center (DMU 60P duo BLOCK, DMG^®, DMP, Frankfurt, Germany) in a dry cutting environment with air-cooling. The milled samples were degreased by ultrasonic cleaning for 15 min in acetone and then rinsed with deionized water. The surface topography, as well as surface profile and Sa, of the machined surface was measured using an interferometer (Wyko NT9300, Bruker Corporation^®, Billerica, MA, USA), and the vertical resolution was 0.1 nm. A rectangular region (1 mm × 1 mm) was chosen in order to compared and analyze the characteristics of surface topography during the ball-end milling process. In order to reduce accidental errors in the measurement, three different regions of the machined surface topography were measured under the same machining parameters. According to the engineering and cutting parameters manual, the value range of feed per tooth is 0~1 mm, the value range of the radial depth of cut is 0~1 mm, the value range of the axial depth of cut is 0~1 mm, the value of the spindle speed is 3000~10,000 r/min, and the lead and tilt angles are 0~45°. In order to deeply understand the influences of changes in cutting parameters on the five-axis ball-end milling process, five values of each factor have been selected according to the actual engineering and cutting parameters manual, as shown in [31]. Three arbitrary experiments under different cutting parameters (see Table 3) were selected to validate the developed surface topography model.

The discrete precision of the workpiece was set to 0.01 mm, which is the less than the minimum value of 0.1f_z and 0.1a_e. The discrete precision of the swept surface is that the angle interval Δθ = 0.2° and Δt is the time it takes the cutter to rotate 2°.

The predicted surface 3D topography and 2D surface profiles were compared with the experimentally measured results from the second trial, as shown in Figure 8 and Figure 9. The results show that the experiments and simulations were very similar to each other, both in the numbers of peaks and valleys and their distribution. There is a certain difference between the simulated and experimental results, such as the fact that the surface topography and profile were smooth in the simulations, while they were rough in the experiments. The main difference between the simulated and experimental results is in the sliding; the abrasion and elastic strain of the material sliding were not considered in the simulation.

According to the simulated and experimental surface topography results, the average roughness, Sa, has been calculated and used to evaluate the developed algorithm quantitatively. The Sa of trials 1#, 2# and 3# are listed in Table 4. The derivation between the simulated and experimental results with a maximum deviation is 13.48%, which shows that the developed model has a high prediction accuracy when predicting the surface roughness. The standard deviations are all less 0.1, which indicates that the measurement results are relatively concentrated, and the values of dispersion are small.

The verified results from Figure 8 and Figure 9 and Table 4 indicate that the surface topography simulation algorithm developed for the ball-end milling of AISI P20 is reliable. Therefore, this surface simulation algorithm can accurately simulate the surface topography and roughness, and thus replace milling experiments.

3.2. Cutting Parameters Optimization

The Sa is a basic requirement of the machining product’s quality, and the surface quality is classified into several grades using the value of Sa. Therefore, the Sa has been taken as a response parameter of the milled surface, which is calculated with the help of simulated surface topography. The feed per tooth (f_z), radial depth of cut (a_e), tilt (β₁) and lead (β₂) angle are selected as the numeric factors because they all have clear effects on the surface topography and can be changed by adjusting the feed rate, tool path and tool orientation. The ranges of the level of each numeric factor were set to feed per tooth—0.1~0.5 mm/tooth, radial depth of cut—0.1~0.5 mm, tilt angle—−8~8° and lead angle—−8~8° for the response analysis using the Box–Behnken design. The designs of the surface topography simulation examples and the simulated Sa values are presented in Table 5.

The results in Table 6 show that the model’s F-value of 145.73 implies the model is significant. There is only a 0.01% chance that a model F-value this large could occur due noise. Values of “Prob>F” less than 0.05 indicate that the model terms are significant. In this case, f_z, a_e, β₁, β₂, f_z × a_e, f_z × β₁, f_z × β₂, f_z² and a_e² are significant model terms. In contrast, values greater than 0.10 indicate the model terms are not significant; in this case, a_e × β₁, a_e × β₂, β₁² and β₂² are not significant model terms, and are reduced to simplify and improve the regression model, as shown in Equation (6). As shown in Figure 10, the comparison of predicted and simulated Sa indicates the high accuracy of the response surface model.

(6)$\begin{array}{l}Sa=-0.2235+0.3955f_z+1.118a_e+0.0078\left({\beta }_1+{\beta }_2\right)-8.2513f_z\cdot a_e\\ -0.3219f_z\left({\beta }_1+{\beta }_2\right)+0.0116{\beta }_1\cdot {\beta }_2+16.681f_z^2+5.936a_e^2\end{array}$

As shown in Figure 10, the predicted Sa using the response surface model was consistent with the simulated Sa using the proposed surface topography model, which indicates the high accuracy of the response surface model. The Studentized residuals are in the range of ±3, which is consistent with the results of the test. In addition, no abnormal data points were found in the residual plot (see Figure 11). The effects of machining parameters on the Sa are shown in Figure 12. The Sa increases with the increase in feed per tooth (f_z) and radial depth of cut (a_e) (see Figure 12a–c). The Sa decreases with the increase in tilt angle (β₁) and lead angle (β₂) (see Figure 12b–d). The feed per tooth (f_z) and radial depth of cut (a_e) have significant effects on material residual height, and the material residual height increases with the feed per tooth (f_z) and radial depth of cut (a_e). The tilt angle (β₁) and lead angle (β₂) can improve the interaction between the tool and workpiece, and reduce the Sa.

The response surface model’s accuracy has here been investigated. Next, the ability of the model to improve the machining efficiency by selecting the cutting parameters is analyzed (see Equation (7)). For this purpose, an optimal cutting parameter with the maximum material removal rate is obtained with an Sa of less than 0.7 by solving the response surface model; the optimal cutting parameters are f_z = 0.2 mm/tooth, a_e = 0.34 mm, β₁ = 7.23° and β₂ = 2.54°. The predicted Sa is 0.6841 μm, in contrast to the simulated Sa, which is 0.7791 μm. The devotion is −12.19%, which denotes the regression model is effective. In contrast to 20#, 21# and 23#, for which the Sa is less than 0.7 μm, the material removal rate can be improved from 0.03 to 0.068 at a growth of 126.67%. The result of the cutting parameters’ optimization is a substantial improvement in the machining efficiency under the surface roughness constraints.

(7)$\left\{\begin{array}{ll}\max \hfill & f_z\times a_e\hfill \\ s.t.\hfill & Sa\le 0.7,0<f_z\le 0.5,0<a_e\le 0.5,0\le {\beta }_1\le 8,0\le {\beta }_2\le 8\hfill \end{array}\right.$

4. Conclusions

In this paper, a novel surface topography simulation algorithm was developed for investigating the surface roughness in the five-axis ball-end milling process. The main conclusions as follows:

• (1). A novel surface topography model was developed using triangular approximation and Z-map methods. The consistency between the simulated and experimental results shows that the model can replace the milling experiment when studying the surface topography and roughness during ball-end milling processes;

• (2). A response surface-reduced quadratic model was developed based on the proposed surface topography simulation algorithm. The model can effectively characterize the correlation of Sa and cutting parameters (i.e., feed per tooth, radial depth of cut, tilt, and lead angles) based on ANOVA results;

• (3). An optimization model was developed for improving the machining efficiency by means of the response surface model. The material removal rate (i.e., product of feed per tooth and radial depth of cut) can be improved effectively under the surface roughness constraints;

• (4). The complex interaction between cutting edge and workpiece is neglected in the proposed model, so the cutting edge trajectory error, cutting edge micro-geometry and workpiece material deformation should be considered in the next study to secure a reliable prediction of surface topography.

Author Contributions

Conceptualization, R.W.; methodology, R.W.; software, R.W.; validation, R.W.; investigation, W.W.; writing—original draft preparation, R.W.; writing—review and editing, B.Z. and D.T.; project administration, R.W. and B.Z. All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

Footnotes

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Figures and Tables

Figure 1. Definition of coordinate systems in ball-end milling of flat surface.

Figure 2. Definition of cutter runout: (a) front view; (b) top view.

Figure 3. Definition of lead and tilt angles.

Figure 4. Definition of cutting edge model: (a) general view; (b) top view.

Figure 5. Z-map model of the workpiece.

Figure 6. Triangular approximation of the swept surface.

Figure 7. Sketched map of surface topography simulation algorithm.

Figure 8. 3D surface topography: (a) experimental; (b) simulated.

Figure 9. 2D surface profile: (a) experimental; (b) simulated.

Figure 10. Comparison of average roughness between surface topography simulation and response surface models.

Figure 11. Residuals of response surface quadratic model: (a) residuals to predict surface roughness; (b) residuals to run number. The numbers represents the amount of similar results, and the colors make it easier to identify corresponding data points in (a,b).

Figure 12. Interaction effect of machining parameters on average roughness: (a) feed per tooth and radial depth of cut; (b) feed per tooth and tilt angle; (c) feed per tooth and lead angle; (d) tilt and lead angles. Colors represent the design points above predicted value.

Table 1

Nominal chemical composition of AISI P20 steel (wt.%).

C	Si	Mn	Cr	Mo	Ni	Fe
0.28~0.40	0.20~0.80	0.60~1.00	1.40~2.00	0.30~0.55	0.05~0.10	Bal

Table 2

Mechanical properties of AISI P20 steel at room temperature.

Density(kg/m3)	Young’s Modulus(GPa)	Hardness(HRC)	Yield Strength(MPa)	Thermal Conductivity(W/m·K)
7800	207	30~36	1140	29.0

Table 3

Machining parameters of ball-end milling experiments.

No.	Spindle Speed, n (r/min)	Feed per Tooth, fz (mm/tooth)	Radial Depth of Cut, ae (mm)	Axial Depth of Cut, ap (mm)	Tilt Angle, β1 (°)	Lead Angle, β2 (°)
1#	5000	0.36	0.25	0.3	−8	4
2#	6500	0.36	0.3	0.15	−12	0
3#	6500	0.56	0.2	0.1	−16	4

Table 4

Experimental and simulated roughness.

No.	Experimental Sa (μm)		Predicted Sa (μm)	Relative Deviation (%)
	Average Value	Standard Deviation
1#	0.8158	0.0656	0.8227	0.85
2#	0.9289	0.0871	0.9087	−2.17
3#	1.8723	0.0985	2.1246	13.48

Table 5

Box–Behnken design and simulated results of average roughness.

Std	Feed per Tooth, fz (mm)	Radial Depth of Cut, ae (mm)	Tilt Angle, β1 (°)	Lead Angle, β2 (°)	Simulated Roughness, Sa (μm)
1	0.1	0.1	4	4	0.0937
2	0.5	0.1	4	4	2.7981
3	0.1	0.5	4	4	1.6139
4	0.5	0.5	4	4	2.9981
5	0.3	0.3	0	0	1.644
6	0.3	0.3	8	0	0.8779
7	0.3	0.3	0	8	0.8779
8	0.3	0.3	8	8	0.8536
9	0.1	0.3	4	0	0.5844
10	0.5	0.3	4	0	3.2895
11	0.1	0.3	4	8	0.5847
12	0.5	0.3	4	8	2.2598
13	0.3	0.1	0	4	0.9095
14	0.3	0.5	0	4	1.6984
15	0.3	0.1	8	4	0.6976
16	0.3	0.5	8	4	1.6937
17	0.1	0.3	0	4	0.5844
18	0.5	0.3	0	4	3.2895
19	0.1	0.3	8	4	0.5847
20	0.5	0.3	8	4	2.2598
21	0.3	0.1	4	0	0.9095
22	0.3	0.5	4	0	1.6984
23	0.3	0.1	4	8	0.6976
24	0.3	0.5	4	8	1.6937
25	0.3	0.3	4	4	0.925
26	0.3	0.3	4	4	0.925
27	0.3	0.3	4	4	0.925
28	0.3	0.3	4	4	0.925
29	0.3	0.3	4	4	0.925

Table 6

ANOVA for response surface quadratic model.

Source	Sum of Squares	df	Mean Square	F Value	p-Value (Prob>F)
Model	21.11	10	2.11	145.73	<0.0001
fz	13.76	1	13.76	949.88	<0.0001
ae	2.33	1	2.33	161.02	<0.0001
β 1	0.35	1	0.35	23.86	0.0001
β 2	0.35	1	0.35	23.86	0.0001
fz × ae	0.44	1	0.44	30.08	<0.0001
fz × β 1	0.27	1	0.27	18.31	0.0005
fz × β 2	0.27	1	0.27	18.31	0.0005
ae × β 1	0.011	1	0.011	0.79	0.3883
ae × β 2	0.011	1	0.011	0.79	0.3883
β1 × β 2	0.14	1	0.14	9.50	0.0064
fz 2	3.07	1	3.07	212.15	<0.0001
ae 2	0.39	1	0.39	26.87	<0.0001
β 1 2	0.03	1	0.03	2.18	0.1615
β 2 2	0.03	1	0.03	2.18	0.1615
Residual	0.26	14	0.014
Lack of fit	0.26	14	0.019
Pure error	0.000	4	0.000
Cor Total	21.37	28

References

1. Zhao, B.; Zhang, S.; Man, J.; Zhang, Q.; Chen, Y. A modified normal contact stiffness model considering effect of surface topography. Proc. Inst. Mech. Eng. J.-J. Eng.; 2015; 229, pp. 677-688. [DOI: https://dx.doi.org/10.1177/1350650114558099]

2. Milan, N.; Sorgato, M.; Parenti, P.; Annoni, M.; Lucchetta, G. Effects of micromilled NiP mold surface topography on the optical characteristics of injection molded prismatic retroreflectors. Precis. Eng.; 2020; 61, pp. 126-135. [DOI: https://dx.doi.org/10.1016/j.precisioneng.2019.10.006]

3. Ghorbani, M.; Movahhedy, M.R. Extraction of surface curvatures from tool path data and prediction of cutting forces in the finish milling of sculptured surfaces. J. Manuf. Process.; 2019; 45, pp. 273-289. [DOI: https://dx.doi.org/10.1016/j.jmapro.2019.07.008]

4. Chen, X.; Zhao, J.; Dong, Y.; Han, S.; Li, A.; Wang, D. Effects of inclination angles on geometrical features of machined surface in five-axis milling. Int. J. Adv. Manuf. Technol.; 2012; 65, pp. 1721-1733. [DOI: https://dx.doi.org/10.1007/s00170-012-4293-y]

5. Zhang, Q.; Zhang, S. Effects of feed per tooth and radial depth of cut on amplitude parameters and power spectral density of a machined surface. Materials; 2020; 13, 1323. [DOI: https://dx.doi.org/10.3390/ma13061323]

6. Peng, F.; Wu, J.; Fang, Z.; Yuan, S.; Yan, R.; Bai, Q. Modeling and controlling of surface micro-topography feature in micro-ball-end milling. Int. J. Adv. Manuf. Technol.; 2013; 67, pp. 2657-2670. [DOI: https://dx.doi.org/10.1007/s00170-012-4681-3]

7. Zhang, C.; Zhang, H.; Li, Y.; Zhou, L. Modeling and on-line simulation of surface topography considering tool wear in multi-axis milling process. Int. J. Adv. Manuf. Technol.; 2015; 77, pp. 735-749. [DOI: https://dx.doi.org/10.1007/s00170-014-6485-0]

8. Zhang, X.; Yu, T.B.; Zhao, J. Surface generation modeling of micro milling process with stochastic tool wear. Precis. Eng.; 2020; 61, pp. 170-181. [DOI: https://dx.doi.org/10.1016/j.precisioneng.2019.10.015]

9. Lotfi, S.; Wassila, B.; Gilles, D. Cutter workpiece engagement region and surface topography prediction in five-axis ball-end milling. Mach. Sci. Technol.; 2017; 22, pp. 181-202. [DOI: https://dx.doi.org/10.1080/10910344.2017.1337131]

10. Hao, Y.; Liu, Y. Analysis of milling surface roughness prediction for thin-walled parts with curved surface. Int. J. Adv. Manuf. Technol.; 2017; 93, pp. 2289-2297. [DOI: https://dx.doi.org/10.1007/s00170-017-0615-4]

11. Omar, O.; El-Wardany, T.; Ng, E.; Elbestawi, M.A. An improved cutting force and surface topography prediction model in end milling. Int. J. Mach. Tools Manuf.; 2007; 47, pp. 1263-1275. [DOI: https://dx.doi.org/10.1016/j.ijmachtools.2006.08.021]

12. Wang, Z.; Wang, B.; Yuan, J. Modeling of surface topography based on cutting vibration in ball-end milling of thin-walled parts in ball-end milling of thin-walled parts. Int. J. Adv. Manuf. Technol.; 2018; 101, pp. 1837-1854. [DOI: https://dx.doi.org/10.1007/s00170-018-3095-2]

13. Li, B.; Cao, Y.; Chen, W.; Pan, J. Geometry simulation and evaluation of the surface topography in five-axis ball-end milling. Int. J. Adv. Manuf. Technol.; 2017; 93, pp. 1651-1667. [DOI: https://dx.doi.org/10.1007/s00170-017-0505-9]

14. Biondani, F.G.; Bissacco, G. Effect of cutting edge micro geometry on surface generation in ball end milling. CIRP Ann.-Manuf. Technol.; 2019; 68, pp. 571-574. [DOI: https://dx.doi.org/10.1016/j.cirp.2019.04.017]

15. Imani, B.M.; Elbestawi, M.A. Geometric simulation of ball-end milling operations. J. Manuf. Sci. Eng.; 2001; 123, pp. 177-184. [DOI: https://dx.doi.org/10.1115/1.1347034]

16. Zhang, W.; Tan, G.; Wan, M.; Gao, T.; Bassir, D.H. A new algorithm for the numerical simulation of machined surface topography in multiaxis ball-end milling. J. Manuf. Sci. Eng.; 2008; 130, 0110031. [DOI: https://dx.doi.org/10.1115/1.2815337]

17. Li, S.J.; Dong, Y.H.; Li, Y.; Li, P.Y.; Yang, Z.C.; Landers, R.G. Geometrical simulation and analysis of ball-end milling surface topography. Int. J. Adv. Manuf. Technol.; 2019; 102, pp. 1885-1900. [DOI: https://dx.doi.org/10.1007/s00170-018-03217-5]

18. Xu, J.T.; Zhang, H.; Sun, Y.W. Swept surface-based approach to simulating surface topography in ball-end CNC milling. Int. J. Adv. Manuf. Technol.; 2018; 98, pp. 107-118. [DOI: https://dx.doi.org/10.1007/s00170-017-0322-1]

19. Yang, L.; Wu, S.; Liu, X.L.; Liu, Z.J.; Zhu, M.W.; Li, Z.H. The effect of characteristics of free-form surface on the machined surface topography in milling of panel mold. Int. J. Adv. Manuf. Technol.; 2018; 98, pp. 151-163. [DOI: https://dx.doi.org/10.1007/s00170-017-0698-y]

20. Chen, H.Q.; Wang, Q.H. Modeling and simulation of the surface topography in ball-end milling based on biharmonic spline interpolation. Int. J. Adv. Manuf. Technol.; 2018; 99, pp. 2451-2466. [DOI: https://dx.doi.org/10.1007/s00170-018-2615-4]

21. Dikshit, M.K.; Puri, A.B.; Maity, A. Optimization of surface roughness in ball-end milling using teaching-learning-based optimization and response surface methodology. Proc. Inst. Mech. Eng. B-J. Eng.; 2017; 231, pp. 2596-2607. [DOI: https://dx.doi.org/10.1177/0954405416634266]

22. Dikshit, M.K.; Puri, A.B.; Maity, A. Modelling and application of response surface optimization to optimize cutting parameters for minimizing cutting forces and surface roughness in high-speed, ball-end milling of Al2014-T6. J. Braz. Soc. Mech. Sci.; 2017; 39, pp. 5117-5133. [DOI: https://dx.doi.org/10.1007/s40430-017-0865-y]

23. Zhou, J.H.; Ren, J.X.; Yao, C.F. Multi-objective optimization of multi-axis ball-end milling Inconel 718 via grey relational analysis coupled with RBF neural network and PSO algorithm. Measurement; 2017; 102, pp. 271-285. [DOI: https://dx.doi.org/10.1016/j.measurement.2017.01.057]

24. Buj-Corral, I.; Ortiz-Marzo, J.A.; Costa-Herrero, L.; Vivancos-Calvet, J.; Luis-Pérez, C. Optimal machining strategy selection in ball-end milling of hardened steels for injection molds. Materials; 2019; 12, 860. [DOI: https://dx.doi.org/10.3390/ma12060860] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30875801]

25. Venkata-Rao, K. Power consumption optimization strategy in micro ball-end milling of D2 steel via TLBO coupled with 3D FEM simulation. Measurement; 2019; 132, pp. 68-78. [DOI: https://dx.doi.org/10.1016/j.measurement.2018.09.044]

26. Sonawane, H.A.; Joshi, S.S. Modeling of machined surface quality in high-speed ball-end milling of Inconel-718 thin cantilevers. Int. J. Adv. Manuf. Technol.; 2015; 78, pp. 1751-1768. [DOI: https://dx.doi.org/10.1007/s00170-014-6759-6]

27. Torres, A.; Amini, C.; Cuadrado, N.; Travieso-Rodriguez, J.A.; Llumà, J.; Vilaseca, M. Experimental validation of ball burnishing numerical simulation on ball-end milled martensitic stainless-steel considering friction and the initial surface topography. J. Mater. Res. Technol.; 2023; 22, pp. 3352-3361. [DOI: https://dx.doi.org/10.1016/j.jmrt.2022.12.100]

28. Cross, K.; Eifler, M.; Klauer, K.; Huttenlochner, K.; Kirsch, B.; Ziegler, C.; Aurich, J.C.; Seewig, J. Determination of the surface topography of ball end micro milled material measures. Eng. Sci. Technol.; 2021; 24, pp. 543-555. [DOI: https://dx.doi.org/10.1016/j.jestch.2020.08.005]

29. Klauer, K.; Eifler, M.; Kirsch, B.; Seewig, J.; Aurich, J.C. Ball end micro milling of areal material measures: Influence of the tilt angle on the resulting surface topography. Prod. Eng.-Res. Dev.; 2020; 14, pp. 239-252. [DOI: https://dx.doi.org/10.1007/s11740-019-00943-x]

30. Arruda, E.M.; Brandão, L.C.; Ribeiro Filho, S.L.M.; De Oliveira, J.A. Integrated optimization using mixture design to confirm the finishing of AISI P20 using different cutting strategies and ball nose end mills. Measurement; 2014; 47, pp. 54-63. [DOI: https://dx.doi.org/10.1016/j.measurement.2013.08.052]

31. Wang, P.; Zhang, S.; Yan, Z.G. Study on surface defects in five-axis ball-end milling of tool steel. Int. J. Adv. Manuf. Technol.; 2016; 89, pp. 599-609. [DOI: https://dx.doi.org/10.1007/s00170-016-9113-3]