1. Introduction
The current machine part technology is becoming more integrated (one machine, a few tools, short production time and cost reduction). Research papers and industrial practice focus on the combination of machining types [
1,
2,
3,
4,
5,
6,
7,
8,
9], including turning [
10], milling [
11,
12,
13,
14,
15], drilling, and threading, using machining centers with forming processes like burnishing [
1,
2,
3,
4,
5,
6,
7,
8,
9,
10,
11,
12,
13,
14,
15,
16,
17,
18]. The main aim of the new technology is to meet the smoothness and strength requirements of the surfaces of the part [
1,
2,
3,
4,
5,
6,
7,
8,
9,
10,
11,
12,
13,
14,
15,
16,
17,
18,
19,
20,
21]. The optimization of the technological parameters of turning [
10] or milling [
5,
6] makes it possible to obtain a satisfactory surface geometry (with low roughness and a useful Abbott–Firestone curve). Given the predefined accuracy, the efficiency and production economics are the main concerns for the technologist. The practically unlimited freedom to control the tooth path trajectory on multi-axis milling machines enables the technologist to use a range of smart tools to obtain the required properties for useful surface shapes [
2,
3,
4,
5,
6,
18].
In addition to the stress responsible for the fatigue strength, the performance properties of the surface [
6] also include a new range of additional requirements defined by constructors and R&D specialists. A more common new requirement is a given level of surface isotropy [
7,
8]. It is beneficial to prepare a surface with an isotropic geometric structure when the operation conditions have not been completely defined. An isotropic surface results in the even wear of components during their use. Given a low roughness for the surface geometry, an isotropic surface provides a uniform and high surface reflectance in all the regions, which is particularly important in modern industrial design [
7,
8,
19]. The surface texture obtained with hybrid production techniques depends to a large degree on kinematic–geometric conditions and the technological parameters defined by the technologist. The surface of a ready product often bears traces of the treatment applied prior to surface finishing [
7,
9,
10,
11,
12]. Working through the remains of the last treatment poses a significant limitation to synergy in hybrid technological processes [
7].
In the production of expensive and complex tools, including molds, masters, and dies, there is no room to search for the optimum treatment parameters that would produce the synergy effect. This paper presents a way to develop a finite element method (FEM) model and research methodology for the assessment of the height parameters of a surface geometry. The model data were validated with experimental data. An analysis covered selected parameters of the surface geometry.
2. Experimental Research
Complex spatial surfaces are commonly machined with ball nose end mills and torus heads. Milling cutters leave traces related to the feed rate, interval, and diameter of the cutting insert. The roughness resulting from the feed per revolution is usually negligibly small. In such cases, the standard procedure is to burnish perpendicular to the direction of the crossfeed cutting after milling [
4,
12,
13,
14,
15,
16]. The burnishing itself is accompanied by phenomena occurring in the contact zone of the burnishing tool. The effect of burnishing to a large degree depends on the hardness of the material and roughness resulting from the previous operations. The basic technological parameters of burnishing include the burnishing force, feed, and number of burnishing tool passes. The surface condition is the least affected by the burnishing velocity.
An experiment was planned in line with the literature guidelines [
9,
10,
11,
12,
13,
14,
15,
16] and based on our own experience in conducting experimental research [
5,
6,
7,
8]. First, 100 × 100 × 20 mm samples made of 42CrMo4 steel and thermally improved to 35 ± 2 HRC were milled using a WNT R1000G.42.6.M16.IK torus (Ceratizit Group, Poland) with six inserts to a diameter of
dp = 10 mm (RD.X1003 MOT–WTN1205). The treatment was conducted on a DMG DMU 60MONOBLOCK machining (Deckel Maho Pfronten GmbH, Pfronten, Germany) center, with a spindle angle of 15°. The milling speed was
vc = 1000 mm/min, with a feed rate
fz = 0.1 mm and a cross feed rate
fwm = 0.5 mm. An approximately 11 mm wide sample was left unburnished to measure the roughness obtained in the forming milling (
Figure 1a).
Measurements of the surface geometry (SG) were made using an AltiSurf A520 multisensor (Altimet, Thonon Les Bains, Frence). It was fitted with a chromatic confocal CL1 sensor, with a range of up to 130 µm and a resolution of 8 nm on the Z optical axis. Measurements were conducted on 2 × 2 mm areas. A scanning resolution of 2 µm along the X and Y axes was experimentally selected. As a result, a cloud of 1 million points was obtained. The mapping of this point cloud was conducted using an AltiMap PREMIUM 6.2 (Digitalsurf, Besancon, France), and then 3D SG parameters were determined:
for each of the measured areas, a threshold value of 0.01–99.9% was determined to delete unreliable surface point data (deleted points were set up as unmeasured values);
then, the surface was leveled (with mean area using the method of least squares);
finally, selected stereometric roughness parameters were determined in compliance with EN-ISO 25178.
Selected SG parameters after milling are listed in
Table 1. The cloud of points for a real surface prepared in this way was exported using the TXT file standard. A recorded image of the milled surface was used for a numerical analysis conducted with the software used for modeling the plastic deformation with the FEM.
The experimental research on burnishing was conducted on a 3-axis MIKRON VCE 500 machining center (Haas Automation INC, Oxnard, USA). A prototype hydrostatic tool with a bellows actuator (West Pomeranian University of Technology Szczecin, Szczecin, Poland) (Patent PL 220528 B1 issued on Nov. 30, 2015) with a ZrO
2 ceramic tip with a diameter of
db = 10 mm was used. The burnishing speed
vb was 8000 mm/min, and cross feed
fwb was 0.12 mm. Burnishing was conducted with a force
Fb = 500 N. The surface after burnishing is shown in
Figure 1b.
After burnishing, selected 3D SG parameters were again determined (
Table 1).
3. Computer Model of Burnishing Process
A randomly selected milled profile from an earlier generated point cloud for a real surface was used to analyze the SG changes following hybrid treatment in the FEM environment (
Figure 2).
A randomly selected milled profile was set as a benchmark for the FEM environment results (
Figure 3). After burnishing, the profile was obtained in the same way as the profile of the milled surface. Selected 2D SG parameters (
Table 2) were determined in compliance with the EN-ISO 4287 standard for both profiles (
Figure 2 and
Figure 3).
Milling with ball nose cutters leaves marks, with the feed per revolution and ball nose cutter diameter as the dominant factors in the roughness. Burnishing smooths the marks.
An earlier numerical simulation [
6] demonstrated in a sensitivity analysis that a mesh gap size below 50 µm did not significantly affect the plastic deformation and residual stress of milled and burnished surfaces. Therefore, bearing in mind the hardware limitations of the physical model developed with Nastran FX, a mesh gap size of 20 µm was used for an individual finite element. Roughness shape changes were observed in a selected milled surface profile in the physical model. The changes visible in
Figure 4 are a result of changing the measurement mesh (2 × 2 µm) into a (20 × 20 µm) FEM mesh.
The samples used in the experiments were 0.4 mm wide and 2 mm long. The modeled physical object was similar to a pyramid with the top cut off. It was 10 mm in height, with a base that was 10 mm wide and 25 mm long (
Figure 5).
The material properties of the thermally improved milled sample were taken into consideration. Nonlinearity was attributed to the sample material to account for the residual stress resulting from prior milling. The characteristics of the material used in the physical model were determined in tensile testing (
Table 3) [
6].
Young’s modulus of X42CrMo4 was determined experimentally; it was E = 210 GPa. Poisson’s number for the model was ν = 0.28 [
6]. In the burnishing model, a ball was pressed into the sample surface with a force of 500 N. The ball was then rolled back perpendicularly to marks left by the milling cutter over 1.5 mm (
Figure 6).
The model of a sample shaped like a pyramid with a cut off top consisted of 470,400 elements connected by 247,222 nodes. The base of the sample was fixed so that displacement was possible only in the model’s plane of symmetry. A total of 2100 rigid shell elements connected by 2227 nodes were used to model the ceramic tip of the burnishing tool. The mechanical properties of ZrO2 were taken from research conducted previously: Young’s modulus E = 220 GPa and Poisson’s number ν = 0.3.
The modeling of the burnishing process (
Figure 6) comprised the following steps:
modeling of how the ball is pressed into the sample surface with a burnishing force of 500 N;
then, the burnishing ball is rolled back over the surface over a 1.5 mm section;
the burnishing force acting on the ball is reduced.
When the ball is rolled back, it rotates by approximately 18°. It took 48 h for Nastran FX to calculate the data for one numerical experiment. The results included the plastic deformation of the total model and coexisting residual stress.
4. Experimental Research Results
Because of hardware limitations and the high requirements for computing power in numerical research, the selected segment of a real milled surface was only 1.5 mm long. Therefore, only unfiltered profile P-parameters were calculated in compliance with EN-ISO 4287 (over one measured section) out of the vertical displacement data recorded in the direction of the burnishing ball rolling back (
Figure 7). The value of the Pa-parameter can be easily interpreted by analogy with the Ra-parameter (the arithmetic mean deviation of the unfiltered profile points away from the midline).
To determine the selected 2D SG parameters (
Table 4) after the burnishing simulation, the set of points was exported again (
Figure 8), and the SG parameters were calculated using AltiMapa 6.2.
The value of Pa = 0.55 was obtained in the burnishing simulation for a selected milled profile from the initial roughness of Pa = 1.56. The differences in parameter Pa were convergent to the average parameter values obtained in experimental research (
Table 2 and
Table 4). The experimental and simulated Pa-parameters differed by less than 5% of their values (
Table 4).
The profile changes in the workpiece surface during the numerical research are shown in
Figure 9. The biggest deformations were observed close to “sharp” summits. Summit deformations caused by the force exerted by the tool in the normal and tangential directions to the surface resulted in valleys being evenly filled in the whole area. In this particular burnishing case, lowering the surface peaks rather than elevating the valleys resulted in the plastic “filling” of the profile. As can be seen, there was a larger difference in the values of the Pp-parameters obtained after burnishing compared to the lowering of the Pv value of the burnished surface in relation to the original value obtained after milling. The trend observed in the numerical research was also confirmed in the experiment.
Small differences (of less than 10%) between the experimental and numerical results were observed only for the parameters most commonly used in engineering: Pa and Pq (
Figure 10). Similarly, a small difference was observed in the filtered profile for roughness parameter Ra (
Figure 11). The values of the 2D SG parameters depended on all the points in a set generated by a profile whose changeability was constant (for Pa) or the root mean square (for Pq) relative to a determined value of the arithmetic mean.
The goodness of fit between the experimental model and experimental results for unfiltered profile amplitude parameters Pt and Pz was worse. In this case, the coordinates of individual points were crucial in terms of their values. The extreme displacement of individual points on a surface decreased when the FEM mesh gap size increased (resolution of the FEM model decreased). Similarly, great differences were observed for the 2D SG parameters that defined the profile ordinate distribution, including the Psk skewness and Pku kurtosis. The observed asymmetry of the unfiltered profile point distribution and concentration of points around the average were greater in the experimental research. This was due to the very high resolution of the profilometer used in the study. The observed direction of changes in the surface character following the use of the hybrid technology was consistent. Some similarity can also be seen in the description of the surface changes expressed with the roughness amplitude parameters (
Figure 11).
As before, an excellent goodness of fit was obtained for the Ra and Rq-parameters (differences smaller than 10%). The numerical research could be used to satisfactorily predict the mechanism and magnitude of the deformation of the roughness summits after milling (the Rp differences were around 21%). Unfortunately, the characteristic of the valley filling in the roughness profile defined by the Rv parameter yielded a difference of almost 10%. Just like before, the reasons for such a great magnitude of errors can be found in the FEM model discretization, which was ten times smaller than that used in the experimental research.
The point distributions of the burnished profile were characterized by different directions of asymmetry (parameter Rsk). The high resolution of the measuring equipment also played a role in this respect. High resolution could be used to better scan and map the narrow and deep valleys of the real burnished surface.
The developed FEM model could satisfactorily determine the Rmr and Rdc parameters of the core ratio of the material. Values based on the profile ordinate distribution given by the Abbott–Firestone curve are often used by technologists to define the required tribological properties of a product surface. To determine the material core parameters, the extreme points at a profile’s summits and valleys are rejected.
5. Summary and Conclusions
The following conclusions can be drawn based on a comparison of the numerical research data and empirical experiment data.
Because of hardware and software limitations, only one pass of burnishing was modeled in the initial research. This is why the surface residual deformation results from the numerical experiment and experimental research were somewhat divergent. The origin of the observed differences lay in the propagation of the residual deformation and stress outside the direct contact zone between the workpiece and burnishing tool. This practically means that each successive burnishing pass occurred on a surface that had already been burnished.
In its present state, the developed numerical burnishing model is a good prognostic tool for predicting the condition of a surface layer. It is particularly useful for residual stress forms and values, which are difficult to measure.
In the numerical research, we assumed the real condition of a milled surface scanned with 0.002 mm steps. The first simulation attempts showed that 0.02 mm steps would be good enough for computer modeling. This simplification did not significantly affect the surface height parameters obtained in the study.
The differences in the parameters typically used to define SG were in most cases below 20%. The assumed simplification made it possible to reduce the computing time and did not significantly affect the ability to predict the SG conditions.
The developed model of the treatment process could also be used to predict the fundamental tribological properties of a surface. The profile point distributions generated with the model were consistent with the experimental data.
In special cases, amplitude parameters, defined by individual points with extreme height values, could not yet be estimated. However, our research is ongoing. First, we plan to increase the resolution for the mesh discretization in the FEM model. Then, we want to introduce additional burnishing passes so that stereometric SG parameters (defined by EN-ISO 25178) can be determined.
Only 3D simulations will provide a reliable tool for the optimization of technological parameters in a hybrid treatment that includes milling and burnishing. A 3D model with improved resolution will enable the prediction of other functional surface properties such as isotropy, wetting, or the ability to maintain lubricants and lacquer coatings.
Author Contributions
D.G.: conceptualization, formal analysis, Investigation, methodology, project administration, resources, supervision, writing—review & editing. S.B.: formal analysis, funding acquisition, validation. Z.G.: conceptualization, formal analysis, funding acquisition, investigation, methodology, project administration, writing—review & editing
Funding
This research received no external funding.
Conflicts of Interest
The authors declare no conflict of interest.
References
- Tersing, H.; Lorentzon, J.; Francois, A.; Lundbäck, A.; Babu, B.; Barboza, J.; Bäcker, V.; Lindgren, L.E. Simulation of manufacturing chain of a titanium aerospace component with experimental validation. Finite Elem. Anal. Des. 2012, 51, 10–21. [Google Scholar] [CrossRef]
- Chen, C.H.; Shiou, F.J. Determination of optimal ball-burnishing parameters for plastic injection moulding steel. Int. J. Adv. Manuf. Technol. 2003, 21, 177–185. [Google Scholar] [CrossRef]
- He, D.; Wang, B.; Zhang, J.; Liao, S.; Deng, W.J. Investigation of interference effects on the burnishing process. Int. J. Adv. Manuf. Technol. 2018, 95, 1–10. [Google Scholar] [CrossRef]
- Bounouara, A.; Hamadache, H.; Amirat, A. Investigation on the effect of ball burnishing on fracture toughness in spiral API X70 pipeline steel. Int. J. Adv. Manuf. Technol. 2018, 94, 4543–4551. [Google Scholar] [CrossRef]
- Berczyński, S.; Grochała, D.; Grządziel, Z. Modeling of Surface geometric structure state after integrated formed milling and finisch burnischinl. Manag. Syst. Prod. Eng. 2017, 2, 131–136. [Google Scholar] [CrossRef]
- Grochała, D.; Berczyński, S.; Grządziel, Z. Modeling of burnishing thermally toughened X42CrMo4 steel with a ceramic ZrO2 ball. Arch. Civ. Mech. Eng. 2017, 17, 1011–1018. [Google Scholar] [CrossRef]
- Dudzińska, S.; Szydłowski, M.; Grochała, D.; Bachtiak-Radka, E. Application of correlation function for analysis of surface structure shaping by hybrid manufacturing technology. Adv. Manuf. 2018, 1, 651–659. [Google Scholar] [CrossRef]
- Bachtiak-Radka, E.; Dudzińska, S.; Grochała, D.; Berczyński, S.; Olszak, W. The influence of CNC milling and ball burnishing on shaping complex 3D surfaces. Surf. Topogr. Metrol. Prop. 2017, 5, 1. [Google Scholar] [CrossRef]
- Kalisz, J.; Żak, K.; Grzesik, W.; Czechowski, K. Characteristics of surface topography after rolling burnishing of EM AW-AlCu4MgSi(A) aluminum alloy. J. Mach. Eng. 2015, 15, 71–80. [Google Scholar]
- Żak, K.; Grzesik, W. Investigation of technological effects of ball burnishing after cryogenic turning of hard steel. Adv. Manuf. Sci. Technol. 2014, 38, 37–52. [Google Scholar] [CrossRef]
- López de Lacalle, L.N.; Lamikiz, A.; Muñoa, J.; Sánchez, J.A. Quality improvement of ball-end milled sculptured surfaces by ball burnishing. Int. J. Mach. Tools Manuf. 2005, 45, 1659–1668. [Google Scholar] [CrossRef]
- López de Lacalle, L.N.; Lamikiz, A.; Sánchez, J.A.; Arana, J.L. The effect of ball burnishing on heat-treated steel and Inconel 718 milled surfaces. Int. J. Adv. Manuf. Technol. 2007, 32, 958–968. [Google Scholar] [CrossRef]
- Rodríguez, A.; López de Lacalle, L.N.; Celaya, A.; Lamikiz, A.; Albizuri, J. Surface improvement of shafts by the deep ball-burnishing technique. Surf. Coat. Technol. 2012, 206, 2817–2824. [Google Scholar] [CrossRef]
- Shiou, F.J.; Chen, C.H. Ultra-precision surface finish of NAK80 mould tool steel using sequential ball burnishing and ball polishing processes. J. Mater. Process. Technol. 2008, 201, 554–559. [Google Scholar] [CrossRef]
- Shiou, F.J.; Chuang, C.H. Precision surface finish of the mold steel PDS5 using an innovative ball burnishing tool embedded with a load cell. Precis. Eng. 2010, 34, 76–84. [Google Scholar] [CrossRef]
- Shiou, F.-J.; Banh, Q.-N. Development of an innovative small ball-burnishing tool embedded with a load cell. Int. J. Adv. Manuf. Technol. 2016, 87, 31–41. [Google Scholar] [CrossRef]
- Randjelovic, S.; Tadic, B.; Todorovic, P.M.; Vukelic, D.; Miloradovic, D.; Radenkovic, M.; Tsiafis, C. Modelling of the ball burnishing process with a high-stiffness tool. Int. J. Adv. Manuf. Technol. 2015, 81, 1509–1518. [Google Scholar] [CrossRef]
- Tadic, B.; Randjelovic, S.; Todorovic, P.; Zivkovic, J.; Kocovic, V.; Budak, I.; Vukelic, D. Using a high-stiffness burnishing tool for increased dimensional and geometrical accuracies of openings. Precis. Eng. 2016, 43, 335–344. [Google Scholar] [CrossRef]
- Tamadon, A.; Pons, D.J.; Sued, K.; Clucas, D. Formation Mechanisms for Entry and Exit Defects in Bobbin Friction Stir Welding. Metals 2018, 8, 33. [Google Scholar] [CrossRef]
- Kaczyński, P.; Ptak, M.; Fernandes, F.A.O.; Chybowski, L.; Wilhelm, J.; Alves de Sousa, R.J. Development and Testing of Advanced Cork Composite Sandwiches for Energy-Absorbing Structures. Materials 2019, 12, 697. [Google Scholar] [CrossRef] [PubMed]
- Ratajczak, M.; Ptak, M.; Chybowski, L.; Gawdzińska, K.; Będziński, R. Material and Structural Modeling Aspects of Brain Tissue Deformation under Dynamic Loads. Materials 2019, 12, 271. [Google Scholar] [CrossRef]
© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).