1. Introduction
The manufacturing industry is facing higher and higher requirements for the quality of machine parts produced by machining. For economic reasons, it has become essential to reduce the number of finishing operations. Much of the research in this area aims to find some convenient and efficient solutions to these problems. Shoulder milling, which combines peripheral milling with face milling, may provide an answer to these expectations. These days, it is one of the most popular milling methods as it is suitable for both roughing and finishing cuts. Another approach to surface quality enhancement in machining is to use new materials for cutting tools or apply special coatings in order to improve the cutting conditions and, consequently, the process efficiency, and, ultimately, the surface quality. Changes to the tools affect the other factors responsible for the cutting process, especially the relative displacements, i.e., the displacements in the tool and workpiece system, the magnitude and distribution of cutting forces, the tool and workpiece temperatures, the friction in the cutting zone, or the minimum chip thickness. Since some of these factors largely influence the surface quality of the finished product, many studies devoted to milling or turning use them to model the surface roughness parameters, particularly Ra and Rt.
Research into the surface texture formation takes into account the contribution of various factors. For instance, a permanent change in the chip cross-section may be due to some axial and/or radial runout of the inserts employed, with this possibly resulting from insufficient geometrical or dimensional accuracy of the cutter body. The effect of the axial runout of inserts is considered in [
1]; therein, the theoretical values of the surface roughness are modeled and compared with the experimental ones. A similar problem is described in [
2]; the analysis concerns the influence of the tool runout on the surface quality in dynamic milling. The modeling and simulation of the surface topography in peripheral milling based on the radial and axial runout of inserts are discussed in [
3]. Differences in the performance of the inserts employed resulting from their axial and radial runouts are responsible for the occurrence of displacements in the tool–workpiece system. Vibration is not a desirable phenomenon in machining because it has a negative effect on the dimensional and geometrical accuracy of the workpiece and the durability of tools and machine tool systems. The displacements in the tool–workpiece system are responsible for lower machine performance and process efficiency. It is, thus, essential to use lower values of the cutting parameters (especially the cutting speed and feed rate), which makes the process take longer. There have been many studies on the influence of vibration on the cutting process. The effects of vibration on the surface roughness parameter Ra are discussed, for example, in [
4]; the research focused on the relationship between the corner radius and the relative displacements. In a further study [
5], a model is proposed to predict surface topography in milling when affected by tool vibration. The method can be used to develop a formula to optimize the surface roughness parameters. Wu et al. [
6] propose the use of the cutting parameters and signal features of vibration measured in milling to predict the surface roughness of S45C steel. Many studies involved using artificial neural networks or machine learning to predict the roughness parameters. For example, in [
7], a model based on an artificial neural network was employed to determine surface roughness parameters. The use of machine learning methods to predict surface roughness parameters in milling was discussed in [
8]. Another approach is proposed in [
9], where the Long Short-Term Memory (LSTM) method was applied to predict surface roughness in milling for S45C steel. Many studies deal with the prediction of surface roughness in machining. Yan et al. [
5] developed a method for predicting the surface roughness in milling based on the information about the tool vibration. The analysis of the influence of the process parameters shows that the feed per tooth has the greatest influence on surface roughness, whereas the effect of the axial depth of cut is the lowest. Many studies have aimed to determine the influence of the key machining process parameters, i.e., the feed rate, the cutting speed, and the depth of cut, on the roughness of milled surfaces [
10]. Kao et al. [
11] indicate that in milling, surface roughness can be predicted from the cutting forces. Their research focused on determining the relationship between the cutting force and the surface roughness of the workpiece machined on a small-size machine tool. The analysis involved employing machine learning algorithms, multi-dimensional linear regression, and a generalized neural regression network to determine the relationship between the cutting forces and surface roughness. In the machining of free-form surfaces, the major parameter responsible for surface roughness is the shape of the machining path. These days, CAM systems are used to design the shape of the tool path. Research in this area [
12] reveals that tool path strategies greatly affect the actual milling time and the workpiece surface roughness. Uzun et al. [
13] analyze the influence of the four most popular machining methods available in CAM systems. In micromachining, the minimum chip thickness is one of the most important and frequently considered process parameters affecting the surface roughness parameters. There have been many studies aimed to predict the minimum chip thickness by using the stagnation point [
14], or to determine it through experiments [
15]. Another factor responsible for the surface roughness of the workpiece is the cutting tool. Surface roughness is dependent on the tool material as well as its micro- and macrogeometry. Different materials are used for cutting tools, and the choice depends on the tool application. Recently, with the rapid development of 3D technology, additively manufactured tools are also available.
Analysis of the literature shows that much attention is given to the prediction of surface roughness parameters, particularly Ra and Rt, as these are the most popular in industry. It can be concluded that the models developed to predict surface roughness are generally used for one type of material and specific cutting conditions.
The novelty of the current article is the development of a model for predicting the surface roughness parameter Ra in shoulder milling. The parameter is determined on the basis of the following factors: the feed per tooth fz, the corner radius, rε, the actual number of inserts used for cutting, the minimum chip thickness, hmin, and the standard deviation of the relative displacements in the tool–workpiece system, D(ξ). The additional novelty of the present work is the description of how to use the surface profile to determine the minimum chip thickness and how to recognize the selected areas on the profile.
2. Materials and Methods
The material selected for the experiments was X37CrMoV5 hot work tool steel, characterized by good ductility, thermal conductivity, fracture toughness under high temperature and water cooling conditions, hardenability, and resistance to tempering. The steel is suitable for die-casting molds and cores, hot scissors and guillotine cutters, press elements, light-metal forming dies, and plastic molding cores [
16]. The material has been purchased from a commercial source that certifies the chemical composition of the material. The chemical composition of X37CrMoV51 steel is provided in
Table 1.
The cutting was performed using a 490-050Q22-08M milling cutter (Sandvik Coromant, Sandviken, Sweden) with five 8 mm × 3.97 mm 490–08T308M–PL–1030 inserts. The inserts, with a corner radius, re, of 0.8 mm, are designed for light roughing to finishing operations. Although the maximum depth of cut for this type of insert is 5.5 mm, the producer recommends that it should not exceed 4 mm. The tool was mounted in the machine tool spindle using an A1B05–4022035 arbor (Sandvik Coromant, Sandviken, Sweden), which is an ISO 40 taper face mill arbor commonly used in machining centers.
The shoulder milling tests were carried out on a VMC 800 vertical machining center (FOP AVIA, Warsaw, Poland). The VMC 800 is a robust and stable machining center composed of four iron cast elements bolted together. The design of the VMC800 vertical machining center is based on a cross table moving in the X and Y axes and a horizontal spindle box moving in the Z axis along the vertical column. During the test, other factors in addition to the feed rates were controlled, namely, cutting speed, depth of cut, and workpiece position, and the relative displacement between tool and workpiece was measured.
A test rig was fitted on the machining center to measure the relative displacements between the tool and the workpiece; the major measuring device was an XL-80 laser interferometer (Renishaw plc, New Mills Woton-under-Edge Gloucestershire, Kingswood, UK). The measurement was performed for 20 s at a frequency of 500 Hz. A total of 10,000 measurement points were used in each measurement. The laser was mounted outside the machine to avoid vibration affecting its position. Before cutting, the tool was run through to measure any external factors that could affect the results. The displacement was then measured during cutting. The forward trend was then subtracted from the measured signal without removing any material.
Before the cutting tests, the insert mounting accuracy was measured using a Heidenhain TT120 tool touch probe (Dr. Johannes Heidenhain GmbH, Traunreut, Germany).
The surface topography, including roughness, was analyzed using a Taylor Hobson: Talysurf CCI—Lite Non-Contact 3D profiler (Taylor-Hobson Ltd., Leicester, UK) with a 20× magnifying lens. A measurement area of 0.8 mm × 0.8 mm provided a 1024 × 1024 pixel array. The 2D surface roughness parameters were determined using a Gaussian filter with a 0.8 cut-off wavelength.
The primary surface profile measurements for use in determining the minimum chip thickness, h
min, were performed using a TOPO 01P profilometer (Instytut Zaawansowanych Technologii Wytwarzania, Cracow, Poland), which is a contact-type instrument equipped with a non-sliding stylus for measuring surface roughness. The profiles were analyzed using the Topografia software running on the profilometer. The profile showed two clear peaks, which corresponded to the scratches made to act as a frame of reference (a system of coordinates) for a specimen. The distance between the scratches was 3.528 mm. In
Figure 1, the plot is divided into two parts, one representing the average surface profile in milling and the other the average surface profile in grinding. The angle between the lines was measured (
α = 0°23′39″). This angle corresponded to the angle of the workpiece surface inclination to the tool during a cutting test. The profile obtained for the milled surface was smoothened in relation to the profile reported for the ground surface so that the former corresponded to zero on the axis of ordinates. The value of h
min (1.34 µm) was measured directly on the profile.
Because of the principle of operation of the test rig used to measure the displacements, the machining was performed in the presence of a cutting fluid. A selective static program was used to determine the effect of each factor separately. One series of the shoulder milling tests was conducted at a constant depth of cut (ap = 0.2 mm), a constant cutting speed (vc = 300 m/min), and a variable feed per tooth ranging from 0.02 to 0.22, which changed every 0.02 mm/tooth. The other series of shoulder milling tests were carried out at ap = 0.2 mm, fz = 0.1 mm/tooth, and vc = 200–400 m/min changed every 20 m/min.
3. Results
The axial and radial runouts of the tool head were calculated for each insert separately on the basis of 50 measurements. This required determining the average axial insert mounting accuracy. For simplicity and greater transparency, it was assumed that the tool length for the most axially protruding insert was 0 mm. Assumptions about tool engagement and zero tool length were used for calculations only to assess the radial and axial runout of the tool. A model of the cutting tool was developed to perform simulations and assess how many inserts were engaged in the surface texture generation at a given value of the feed per tooth. The model can be used to calculate the load each insert is subjected to, the actual feed per tooth, the instantaneous depth of cut, and the cross-sectional area of the material removed from each insert. Some of the calculation results are given in
Table 2.
From
Table 2, it can be concluded that the insert mounting accuracy and the feed per tooth have a substantial influence on the material removal process. As can be seen from
Table 2, only three out of five inserts were engaged in the cutting at a feed of 0.04 mm/tooth. When, however, the feed exceeded 0.12 mm/tooth, the material was removed by all the five inserts. The third insert at a feed rate of 0.12 mm/tool removes less material than the value of the h
min parameter, and at a feed rate of 0.14 mm/tool, the insert exceeds the value of the h
min parameter. During the cutting tests, the relative displacements in the tool–workpiece system were measured using an XL-80 laser interferometer (Renishaw plc, New Mills Woton-under-Edge Gloucestershire, Kingswood, UK). Some relative displacements in the tool–workpiece system are illustrated in
Figure 2.
Figure 2 shows a signal of the relative displacements registered in the tool–workpiece system during a test measurement under the following cutting conditions:
vc = 300 m/min,
fz = 0.1 mm/tooth, and
ap = 0.2 mm. As can be seen, the signal varies over time. During the first phase, the displacement is stable because this represents the time when the tool approaches the workpiece and the tool is not engaged in cutting (0–2.7 s). In the second phase (2.7–6.8 s), when the tool performs shoulder milling, the displacement increases. The last phase (6.8–11 s) corresponds to the final stage of cutting, during which the tool is no longer in contact with the workpiece. From
Figure 2, it can be concluded that the relative displacements increase with the depth of cut. Thus, it is clear that the displacements occurring in the first phase increase with increasing width of cut until the maximum engagement of the cutting tool is achieved, i.e., the whole diameter of the tool is used. It was then vital to determine how the cutting speed, v
c, and the feed per tooth, f
z, affected the standard deviation of the relative displacements (D(ξ)) for the selected material.
Table 3 shows the standard deviation of the relative displacements measured in the tool–workpiece system, D(ξ), given in μm, for the cutting performed at
vc = 300 m/min,
ap = 0.2 mm, and a variable feed per tooth.
Table 3 shows the values of the standard deviation of the relative displacement in the tool–workpiece system, i.e., the parameter D(ξ) expressed in μm, for the selected workpiece material machined at a constant cutting speed, v
c, of 300 m/min and a variable feed per tooth, and also at a constant feed of 0.1 mm/tooth and a variable cutting speed. From the data in
Table 3, it is clear that increasing the feed per tooth had a negative effect on the D(ξ) parameter. The higher the feed per tooth, the higher the standard deviation of the relative displacement. In the entire range of feeds per tooth, the displacement rose by 0.63 µm. When the feed varied between 0.02 mm/tooth and 0.08 mm/tooth, the parameter D(ξ) first increased from 1.07 µm to 1.53 µm and then dropped slightly to 1.49 µm. At feeds of 0.1–0.16 mm/tooth, D(ξ) improved from 1.49 µm to 1.61 µm to fall slightly at a feed of 0.16 mm/tooth. When the feed was 0.2 mm/tooth, there was a rapid decline in the standard deviation of the displacement to 1.48 µm, followed by a significant increase of 0.23 µm.
At feed rates of 0.10 mm/tooth and 0.16 mm/tooth, D(ξ) decreases because the third and fourth inserts were involved in surface texture generation and the minimum feed rate was higher than the hmin parameter.
Two characteristic maximum values of the standard deviation of the relative displacements (2.6 μm and 3.45 µm) can be observed in
Table 3. They were registered at cutting speeds of 260 m/min and 360 m/min, respectively. The points correspond to the two most efficient ranges of operation of the CoroMill 490 shoulder mill taking into consideration the displacements in the tool–workpiece system. The smallest displacements were reported at cutting speeds of 200–220 m/min, 280–340 m/min, and 400 m/min. The average standard deviation of the relative displacements for the whole range of cutting speeds was 2 μm.
The experiment consisted of machining 22 samples, 11 at constant feed and 11 at constant cutting speed.
Table 3 shows the values of the standard deviation of the relative displacement in the tool–workpiece system for the whole range of the cutting path. Discrepancies in results for the same machining parameters may be the result of material heterogeneity. The difference in the standard deviation of relative displacements for the same machining parameters is 0.39 µm.
The next series of tests aimed to determine the influence of the cutting conditions on the minimum chip thickness (h
min) in shoulder milling.
Table 4 shows the values of the minimum cut thickness measured at different feeds per tooth and cutting speeds.
Table 4 shows the influence of cutting speed and feed rate on the minimum chip thickness, h
min. From the data, it can be seen that in the initial phase, an increase in cutting speed (
vc = 200–260 m/min) contributed to a decrease in h
min. The parameter h
min increased when the cutting speed exceeded 280 m/min. The downward trend continued until the cutting speed reached 380 m/min. The parameter remained stable until v
c was 260 m/min. When the cutting speed exceeded this value, h
min increased and then decreased. The analysis of the influence of the feed per tooth, f
z, on the minimum chip thickness, h
min, for X37CrMoV51 steel shows that the parameter h
min fluctuates up and down. From the relationship between f
z and h
min, it can be concluded that the parameter h
min increases as the feed per tooth increases.
Another objective of this study was to look at the influence of the cutting process parameters in shoulder milling on the surface texture of the workpiece. The shoulder milling tests were carried out on an AVIA VMC800 vertical machining center (FOP AVIA, Warsaw, Poland). First, the tests were conducted under dry cutting conditions, i.e., without the use of cutting fluid, under specific cutting conditions assumed for a given series. The measurement results obtained for different cutting parameters were analyzed thoroughly. The data in the form of tables, plots, and maps were used to determine the 2D surface roughness parameters. Some isometric views of the surfaces of the X37CrMoV51 steel specimens are shown in
Figure 3. The surface topography was analyzed using a Taylor Hobson: Talysurf CCI—Lite Non-Contact 3D profiler measurement interferometer (Taylor-Hobson Ltd., Leicester, UK).
As can be seen from
Figure 3, at a feed of 0.02 mm/tooth, random roughness was observed. The surface roughness increased with increasing feed per tooth; the peak-to-valley distance corresponded to the distance the tool traveled during a single spindle rotation (feed per revolution). At
fz = 0.06 mm/tooth, there were no sharp peaks; burrs were observed instead. A further increase in the feed per tooth caused an increase in both the surface roughness and the surface waviness. From the isometric views, it is evident that at
fz = 0.1–0.16 mm/tooth and 0.22 mm/tooth, a two-directional cut was reported because the milling cutter milled forward and backward.
Table 5 shows the 2D surface roughness parameters for shoulder milling. The lowest value of the maximum profile height, Rz, was recorded at a feed rate of 0.1 mm/tooth, while the highest value was recorded at a feed rate of 0.22 mm/tooth. The profile characteristics, Rt and Rv, were similar to Rz. The arithmetic mean height, Ra, (
Figure 4) increased steadily with increasing feed from 0.2 μm at 0.02 mm/tooth to 0.321 μm at 0.2 mm/tooth. At a feed of 0.22 mm/tooth, Ra decreased slightly to 0.304 μm. The root mean square deviation, Rq, and Rp, the maximum profile peak height, were reported to be identical to Ra over the whole range of feeds. The mean height of the profile elements, Rc, decreased with increasing feed. After the feed exceeded 0.06 mm/tooth, Rc fluctuated around 0.55 μm and then increased when the feed was 0.1–0.18 mm/tooth. At feed rates higher than 0.18 mm/tooth, the parameter Rc decreased. The mean width of the profile elements, RSm, first increased and then decreased with increasing feed per tooth. At
fz = 0.06 mm/tooth, RSm decreased to 0.03 mm. Then, at feed rates of 0.08–0.12 mm/tooth, RSm increased to 0.048 mm and then decreased to 0.034 mm. A further increase in feed to 0.2 mm/tooth caused the parameter RSm to reach a maximum of 0.058 mm and then drop to 0.04 mm.
The next stage of the study was to develop a model for predicting the surface roughness parameter, Ra, in shoulder milling. First, relationship (1) was formulated as a linear function to link the minimum chip thickness with the feed per tooth. The coefficients used in the formula were selected in such a way that the feed could be expressed in mm. As there were substantial fluctuations in h
min, the model took into consideration only the minimum chip thickness, which increased with increasing feed per tooth.
Then, relationship (2) was formulated in the form of a power function to calculate the parameter D(ξ) using the value of the feed per tooth. The coefficient of correlation between the theoretical results determined from relationship (2) and the experimental data obtained through tests was 0.89.
The last stage of the study was to formulate relationship (3), which is the model for predicting the surface roughness parameter, Ra, on the basis of the number of inserts engaged, the corner radius, the minimum chip thickness, and the relative displacements. It was assumed that the relative displacements combined with the geometry of the insert (wiper for finishing) contribute to lower roughness as they help remove some of the machining debris.
where f
za is the actual value of the feed per tooth dependent on the actual number of inserts taking part in the cutting process, z
a.
k is the material coefficient describing the contribution of the minimum chip thickness to the surface roughness formation. The material coefficient, k, is the ratio of the average peak-to-valley distance on a surface profile to the minimum chip thickness.
The analysis to determine the effect of the selected factors required finding the value of the material coefficient, k. The average value of this coefficient (k = 0.4) was calculated using the data from measurements of 11 X37CrMoV51 steel specimens.
The model for predicting the surface roughness parameter, Ra, in shoulder milling was verified using relationship (3) by substituting the actual feed per tooth with relationship (4), the material coefficient, k, with relationship (5), and the parameters hmin and D(ξ) with relationships (1) and (2), respectively. The predicted surface roughness Ra was thus dependent on the feed per tooth, fz, the corner radius, rε, the actual number of inserts taking place in the material removal process, and the parameters hmin and D(ξ).
Table 6 shows the data used to verify the proposed model. As can be seen, the difference between the measured values and the calculated values is very small.
Figure 5 compares the plots obtained from the experiments with those from the calculations.
The red line in
Figure 5 represents the predicted values of the surface roughness parameter Ra in shoulder milling, which were calculated using the proposed model based on relationship (3). The blue line shows the relationship between the actual surface roughness and the feed per tooth in shoulder milling at a cutting speed of 300 m/min and a depth of cut of 0.2 mm. The X37CrMoV51 steel specimens were measured with a Taylor Hobson Talysurf CCI–Lite Non-Contact 3D Profiler (Taylor-Hobson Ltd., Leicester, UK).
From
Table 6 and
Figure 5, it is evident that for feed rates ranging from 0.02 to 0.18 mm/tooth, the predicted surface roughness Ra determined using Equation (3) is similar to the actual surface roughness. At a feed rate of 0.2 mm/tooth, however, the predicted values are slightly different from the actual ones. The coefficient of correlation between the theoretical and experimental data was 0.96. This value confirms that the developed model is correct.