Next Article in Journal
A Self-Healing Hierarchical Fiber Hydrogel That Mimics ECM Structure
Next Article in Special Issue
Tensile Deformation Behavior of Typical Porous Laminate Structure at Different Temperatures
Previous Article in Journal
The Effect of Different Plant Oil Impregnation and Hardening Temperatures on Physical-Mechanical Properties of Modified Biocomposite Boards Made of Hemp Shives and Corn Starch
Previous Article in Special Issue
Deformation of Bioinspired MXene-Based Polymer Composites with Brick and Mortar Structures: A Computational Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Researches and Simulation of Elastic Recovery Phenomena during Roller Burnishing Process of Macro-Asperities of Surface

by
Agnieszka Kułakowska
* and
Łukasz Bohdal
Department of Mechanical Engineering, Koszalin University of Technology, Racławicka 15-17 Street, 75-620 Koszalin, Poland
*
Author to whom correspondence should be addressed.
Materials 2020, 13(22), 5276; https://doi.org/10.3390/ma13225276
Submission received: 12 October 2020 / Revised: 9 November 2020 / Accepted: 18 November 2020 / Published: 21 November 2020
(This article belongs to the Special Issue Advances in Mechanical Testing of Engineering Materials)

Abstract

:
The paper presents preliminary studies of a new innovative surface treatment method—the process of roller burnishing of macro-irregularities of the surface. As part of the work, the possibility of plastic shaping of the surface macrostructure with indentations (plateau structure), which will show anti-wear properties through appropriate surface shaping and the compressive stress state in the product’s top layer, was investigated. The essence of the paper is the analysis of one of the aspects of the application of this processing method, i.e., the influence of the elastic recovery of the product on its technological quality measured by dimensional deviation. The main objective of the work is to develop adequate methods and mathematical models to enable the design of the macro-asperities of the surface burnishing process to maintain the dimensional tolerance of the shaped parts. The results of dependencies of elastic recovery of the asperities and the deviation of height, Δht, for sample depths of burnishing were presented. The model tests of the elastic recovery of the model material using the visioplasticity method show that with the increase of the value of the vertical surface asperities, the value of the elastic recovery of the material decreases. The increase of the deviation of the asperities’ height causes a decrease in the value of elastic recovery. With the increase of the value of the vertical angle of the surface roughness, the value of the elastic recovery of the material is smaller.

1. Introduction

Roller burnishing is used increasingly as a finishing operation [1,2,3,4,5,6]. It is known that the surface state after previous treatment has significant influence on its quality after roller burnishing [7,8].
The effect of burnishing is to reduce the roughness of the profile while hardening the surface layer of the product and the formation of compressive stress in it [1,2,3,9]. In the process of burnishing several interacting processes occur simultaneously, phenomena present in the material deformation areas affected on the conditions in the contact zone, and those on some of the burnishing conditions, which also affects its complexity [5,10].
Review of the literature shows that work on the burnishing process has been conducted by many researchers. Several works have investigated the effect of the burnishing process to improve the properties of the parts’ mechanical properties, e.g., increased hardness [2], higher wear resistance [6], surface quality [9,11], and increased maximum residual stress in compression. The parameters affecting the surface finish are burnishing force, feed rate, ball material, number of passes, workpiece material, and lubrication [6]. References [1,2,3,4] showed that the burnishing force and the number of tool passes are the most predominant of the parameters that have an effect on the surface roughness of the workpiece during the burnishing process.
Determining roller burnishing parameters also causes difficulty in achieving the required quality, because the measurement of process parameters (determining the technological quality such as temperature field, stress field, structural changes, the change in the dimensions under the influence of elastic deformation, etc.) during the burnishing process in the current state technique is impossible.
The process parameters can only be inferred from the properties of the product after the process. As a result of burnishing, there are dimensional changes caused by the following factors:
  • Plastic deformation and the decrease in surface roughness,
  • The concentration and structure squeeze,
  • An axial elongation of the product due to plastic deformation and residual stress generated in the surface layer [6],
  • Moving the material after moving away the burnishing element and increasing the asperities height.
Despite numerous works on burnishing treatment, the problem of springback of the material after unloading (or moving away) the burnishing tool, has been considered in a few studies [12,13]. The problem of elastic recovery is most often occupied by researchers working on pressing an indenter into the material [14,15,16,17,18,19], and less frequently, rolling [20] or the process of sheet metal forming [21] or cutting [22].
Stilwel and Tabor [18] found that when a spherical indenter is pressed into a plastic material, it creates a cavity, then after it is pulled back, there is an elastic return of the material of both the indenter and the cavity. Shape changes can be described in relation to the Hertz equation for the elastic deformation of a spherical surface. Chodór, in his work [22], using the results from numerical analyses, shows the elastic recovery of the material during the process of micro-cutting and sliding burnishing. Skalski [23] noticed that during the propagation of the plastic area (in the case of pressing a rigid cylinder into a body with a flat surface), a narrow elastic area appears under the cylinder on a plastic substrate. This area does not plasticize even during increasing load.
Despite the large number of works on burnishing, elastic recovery has not been considered. There are works concerning elastic recovery during the scratch test [14,24] but not during roller burnishing. Also, in many papers, modeling and numerical analysis of the process were presented [25,26,27,28,29,30,31] but without considering the outline of the surface after previous treatment and material elastic recovery. Reference [32] describe the effects of ball-end hard milling to a variation of the lubrication situation. At the reversal, roller burnishing process parameters on the surface point to micro-contacts at the maximum rolling speed elasto-roughness of the raceways of constant velocity joint shafts. The authors of Reference [33] developed a three-dimensional (3D) numerical finite element model of the ball burnishing process including real surface integrity descriptors resulting from a ball-end-milled AISI 1038 surface in the target workpiece. Specifically, its periodical topological features were used to generate the surface geometry. Different models varying the effect of the coefficient of friction and the direction of application of burnishing passes with regards to the original milling direction were calculated. Reference [34] reports the results of experimental studies on the impact of ball burnishing parameters on the roughness, microstructure, and microhardness of the surface layer of laser-cut C45 steel parts. The distribution of residual stresses generated in the surface layer of these parts were analyzed. The superficial effect of hardening caused after vibration-assisted ball burnishing and its consequences in the tensile behavior of a carbon steel material were studied in Reference [35]. Experiments using different amplitudes and new forces were encouraged to obtain more information about how the material can be modified optimally through vibration-assisted ball burnishing.
Burnishing of the macro-irregularities of a surface is a variation of the typical burnishing process, however it differs in the scope of influence on the surface layer of the product. The research on this new process was undertaken due to the conducted preliminary research, which gave very favorable operational effects of the surface layer shaped in this way. It can be seen that the depth of strengthening the surface layer can be even a few mm, while for a typical burnishing process, only a part of a mm or a few µm. Additionally, the literature shows how to solve the problem of the kinematics of the roughness of the material for the triangular asperities of the surface crushed (burnished) with a flat punch under various friction conditions [36,37,38]. Kukiełka [10,39,40] stated that the kinematics of the smoothing of the asperities depends only on the apex angle of the rough surface and determined the angle boundaries for which there is a change in the flow mechanism.
The author of Reference [39] qualitatively distinguished three different cases of flow of the surface layer of material. The considerations were carried out with the assumption that the surface asperity after the previous treatments is triangular, symmetrical, and regular. The research shows that in the burnishing process, these cases depend only on the vertical angle θ of the asperity (Figure 1):
  • For vertical angles θ ≤ 80°, the material core remains undistorted. The indentations of the asperity do not rise, and the deformation of the material occurs only within the unevenness. A gap (discontinuity planes) with the depth 0.5 Rt that separates deep inequalities can be seen. As a result of the material flowing sideways, the surface flattens out.
  • For the vertical angles 80° < θ < 145°, the zone of plastic deformation includes the material core. Gaps are still visible at the junction of adjacent burrs but of less depth. The asperity valleys rise and total strain in the contact zone of neighboring material overflows.
  • For angles θ ≥ 145°, as a result of deformation of asperity and the material core, the surface is evened out. This is not at the expense of the material outflows towards the sides of the asperity. There are no planes of material discontinuity in the surface layer. The rise of the indentation is equal to the decrease in the top of the asperity.
A significant problem is to ensure the regularity of the outline and the appropriate vertical angle θ of the shaped irregularities in the treatment preceding burnishing. Due to the type of the machined part and its intended use, the value of the feed, ft, in the preceding machining and the value of the vertical angle θ of the asperity should be selected each time.
In previous paper [12] on the initial stages of the experiment related to elastic recovery issues, we analyzed the influence of the deviations of the height and distance of the asperities after turning on the deviation of the height of the asperities of the surface as well as elastic recovery after burnishing. Research was performed for different burnishing depths. The calculations were performed for apex angles θ = 90° and θ = 120°. The received results showed that the elastic recovery increases as the burnishing depth increases. Each of the asperities were analyzed separately. The analysis of the state of knowledge shows that none of the researchers analyzed the problem of elastic recovery, which is a novelty and originality of the presented article. In the paper, FEM (Finite Element Method) analysis with using the dynamic explicit method for the roller burnishing process was established. The surface state after turning (as previous treatment) was taken into account. The results of dependencies of elastic recovery of the asperities and the deviation of height, Δht, for exemplary burnishing depth is presented. The influence of the surface after previous treatment on the burnished product quality can be analyzed by developed numerical algorithms. It allows for better understanding of the phenomena which occur in the tool–workpiece contact zones. The basis for the development of guidelines for the selection of the conditions of rolling and burnishing processes considering the required technological quality of the product can be formulated. This paper consists of numerical simulations and analysis, model test, and the research on the steel C45, with the aim to understand elastic recovery as well as analysis of elastic recovery of the asperities with different vertical angles.

2. Idea of Elastic Recovery

It is assumed that during the burnishing process, deformations of the material occur in terms of elastic and visco-plastic strains [41]. It is known that visco-plastic strains will remain in the material, while elastic strains change [41,42,43]. After the burnishing tool has passed, the phenomenon of elastic recovery occurs and visco-plastic deformations remain in the processed material, while the elastic deformations change and settle on the level corresponding to the new state of equilibrium of the object. Thus, the states of displacements, stresses, and strains in the tested object during and after the end of the process differ significantly, because after the process is completed, the nodes are displaced in the elastic range. The value of elastic deformations will decrease by the value Δϕ(E), while the elastic deformations will remain at the value ϕ(E) in the surface layer of the product. This is shown in Figure 2.
During the process, the burnishing tool is plunged to the depth a. It can be noticed that the contact point of the tool and the workpiece material A, as a result of this depression, moves to the point A′. Then, after the rotation of the workpiece and the burnishing tool, the load is unloaded and point A′ is moved to position A″. The radial distance between points A′ and A″ is marked as Δui and represents the elastic recovery of the material after burnishing, and φ is the actual deformation of the material.
The process of roller burnishing is treated in this study as a process with unknown boundary conditions in the tool–workpiece contact area. It is a geometrical and physical boundary and initial value problem. During experimental research and simulations of the process, special attention was drawn to the elastic strain of the material, which takes place after moving away the burnishing element, called elastic recovery, Δui. The state of the surface layer after previous treatment and after roller burnishing are considered together. They take into account two asperities, A1 and A2, after the turning process, which are burnished. Asperity A1 after turning was characterized by constant height (h) and distance (s). Asperity A2 possesses different values of the height or distance. The difference between the asperities was called height deviation after turning (Δht) and between distances: distance deviation after turning (Δst). Similarly, deviations were defined after burnishing and marked as Δhb and Δsb. The schematic diagram of elastic recovery in joint analysis of the surface after turning and burnishing is presented in Figure 3.

3. Researches of the Influence of Height and Distance Deviations on the Elastic Recovery

3.1. Experimental Researches—Materials and Methods, Discussion

The experiment was conducted for the following input factors: height deviation (Δht) and distance deviation (Δst) of the surface asperities after rolling, determined in compliance with Table 1, and the vertical angle of asperities was constant, θ = 60°. The five-level rotatable experiment plan was used to set up the values of the deviations.
The values of the input factors were changed in the following ranges: Δht = 0.01–0.04 mm, Δst = 0.05–0.2 mm. The tests were repeated three times for each point of the experiment plan. Tests were performed at room temperature (T = 19 °C). The external surfaces of the samples made of steel C45 were prepared in the turning process on a NEF 520 numerical lathe in “TEPRO” Vacuous Technology Plant in Koszalin, with a VCMT 160402 machining plate. The following machining parameters were used: velocity vt = 200 m/min and feed ft = 0.15 mm/rev in rough machining and vt = 200 m/min, ft = 0.1 mm/rev in finishing. Six samples were made on the shaft (Figure 4). Each sample included two asperities, which differed in the value of the height and distance deviations, according to the experiment plan. The accuracy of the shaping of the outline was verified on a Werth optoelectronic microscope. Figure 5 shows an example of a shaft with samples prepared in accordance with the experiment plan, Table 1 (first 6 implementations). The values of the height and distance deviations obtained in measurements were not significantly different (at the significance level α = 0.05) from those planned on the basis of single factor analysis of variance.
Next, the samples were burnished using a Fette T27 tangent burnishing head with high stiffness to reduce the influence of dimensional changes of the head on results. The burnishing depth for each sample was calculated separately, while assuming that it is a = 0.5 Rt, where Rt is the asperity height without a deviation (constant one).
After the burnishing process, the values of parameters Δu0 were calculated and profiles with profilometer T8000 were obtained. Exemplary profile of burnished asperities is presented in Figure 6.
The results of the measurements were developed statistically in the Experiment Planner program developed by Kukiełka L. and Kukiełka S. [44,45,46,47,48,49,50]. The following regression equation is obtained of the dependence of the material elastic recovery after burnishing (Δu0) from the deviations of height (Δht) and distance (Δst) of asperities after turning, for correlation coefficient R = 0.985:
Δ u 0 = 0.0007 Δ s t 0.013 Δ h t 2.423 × 10 14 Δ s t Δ h t + + 0.0048 Δ s t 2 + 0.12 Δ h t 2
In order to determine the form of Equation (1), appropriate analytical and experimental tests were carried out, in the following stages: determination of the sets of tested, constants, disturbing and resulting factors, determining the range of variability (research area) of the examined factors, adoption of the class of the mathematical model of the research object, coding of researched factors, implementation of the actual research, experiment plan, the results of the experiment, elimination of results with a gross error, calculating the inter-row variance, checking the homogeneity of variance in the sample, calculating the coefficients as a function of the regression, statistical analysis of the regression function, study of the significance of the multidimensional correlation coefficient, checking the adequacy of the mathematical model, and decoding the regression function. The detailed procedure is given in References [39,41,44]. The graph of the function according to Equation (1) with the correlation coefficient R = 0.86 is presented in Figure 7.
It is visible that an increase of the deviation of the asperities’ height causes a decrease in the value of elastic recovery, while the opposite situation occurs when the deviation of the distance is analyzed. Knowledge about influence of determined profile parameters on the physical phenomena that occur in the burnishing process is unknown. It is important because they impact on the quality of the formed product’s surface layer, tribological properties, and fatigue issues. Many physical phenomena whose observation or measurement is very difficult or impossible take place during processes of burnishing, such as the pressure in contact zones, the friction forces, strain and stresses, elastic recovery, the slip and contact zones, and states of displacement. After passing the burnishing tool, visco-plastic deformations remain in the workpiece, while the elastic deformations change and settle at the level corresponding to the new equilibrium state of the object. Thus, the states of displacements, stresses, and strains in the tested object during and after the end of the process differ significantly, because after the process is completed, the nodes are displaced in the elastic range. In the case of higher asperities, the displacement is greater than in the case of those with a height deviation. The opposite is the case for a distance deviation. These issues are very complex and require further analysis.

3.2. FEM Modeling and Numerical Results

The computer models (2D and 3D) of the process in the real scale were characterized by the possibility of changing the value of the height deviation (Δht) and the distance deviation (Δst). Exemplary calculations were performed for the apex angle θ = 60° and the asperities feed f = 2.7 mm, and the deviations were changed according to Table 1. Numerical analyses were conducted with the aid of Ansys Structural, LS-Dyna and LS-Prepost programs. In the analyzed case, the object was treated as an elastic body (in terms of reversible deformations) and visco-plastic (in terms of irreversible deformations), while the tool (roller) was treated as a perfectly rigid body. Based on the results obtained, various analyses can be carried out. In the discussed case, only some of the results were presented in order to determine the impact of the burnishing process of the rough surface (with regular triangular asperities) on the value of elastic material return after the burnishing process. Moreover, in incremental material models, it was assumed that the object was made of C45 steel, for which the cumulative value of the yield stress in 2D models is described by the regression equation according to Equation (2):
σ p = 976.95   ( 0.0667 + φ i ( VP ) ) 0.2052 [ 2.33 × 10 6 ( φ ˙ i ( VP ) ) 2 + 8.83 × 10 5 φ ˙ i ( VP ) + 1 ] 1.337
where: yield stress is σp, intensity of the real strain is φi, and strain rate is φ ˙ i .
The necessary parameters for this material are Young’s modulus E = 210 GPa, Poisson ratio ʋ = 0.29, yield stress Re = 425 MPa, and hardening modulus ET = 1024 MPa. In 3D models, object discretization was performed using finite elements of the object with a linear shape function, shell elements (tool), and contact elements (contact of the tool with the object). The nonlinear, anisotropic Coulomb friction model and the coefficient of static friction μs = 0.1 and dynamic friction μd = 0.05 were adopted. Then, adequate boundary conditions were set—mainly concerning the displacements: tool cavity (up to the value of 0.5 of the roughness height) and the degrees of freedom for the workpiece, as well as the set rotational speed ensuring the required burnishing speed vb = 0.52 m/s. In this case, the shaft with surface asperities was rotated at n = 477 rpm, then the burnishing element (roll) was moved to the burnished surface. The contact of the roller with the machined surface caused the tool to rotate. In the next step, methods of solving the discrete equation of motion were introduced. For this purpose, in accordance with the adopted algorithm of explicit solution, at each analysis step, an approximation of the column vector of acceleration and velocity increment as a function of the column vector of displacement increment was introduced, all with the following dimension integration and were given the required number of steps, as well as the total calculation time t = 3.5 s. The geometrical model and mesh grid are shown in Figure 8. The 3D FEM model was developed according to Hill plasticity condition and low strain rates. Thermal strains for this model were approximately 2%.
The temporary yield stress in the 3D model was described with the aid of the Cowper–Symonds model. A Cowper–Symonds model allows for linear isotropic (β = 1), kinematic (β = 0), or mixed (used in this study) (0 < β < 1) plastic strain hardening, and the effect of the plastic strain velocity is given by the following power relation:
σ p = [ 1 + ( φ ˙ i ( p ) / C ) ] m   ( R e + β E p φ i ( p ) )
where β is the plastic strain hardening parameter, R e (MPa) is the initial static yield point, φ ˙ i ( p ) (s-1) is the plastic strain velocity, C (s−1) is the material parameter defining the effects of the plastic strain velocity, m = 1/P is the material constant defining the sensitivity of the material to the plastic strain velocity, φ i ( p ) is the plastic strain intensity, and E P = E T E E E T is a material parameter dependent on both the plastic strain hardening modulus, E T = σ p / φ i ( p ) , and Young’s modulus, E.
To find the value of material elastic recovery in performed computer simulations, one had to select areas when the tool affected the subject, and after it moved away. Two kinds of analysis were performed. The first, when Δht (height deviation) values after the turning process were Δht = 0–0.04 mm and Δst (distance deviation) values after the turning process were Δst = 0, and the second, when Δst = 0–0.2 mm but Δht = 0. On the basis of computer simulations, the graphs of dependencies of elastic recovery of the first Δu1 and the second Δu2 asperity and the deviation of height, Δht, for the sample after burnishing were elaborated. The graph of dependencies of elastic recovery Δu1 and Δu2 in the asperities A1 and A2 from distance deviation Δst is presented in Figure 9.
It has been observed that together with the increase of the deviation of surface asperities’ height, the elastic recovery decreases. Comparing the values of elastic recovery for an asperity with constant geometry A1 and changeable geometry A2, it can be seen that the value of elastic recovery was higher for asperity A1.
Computer simulations were also performed to validate the experimental studies according to the experiment plan presented in Table 1. An example of the results obtained from the simulation is shown in Figure 10, including a geometric model with a mesh of finite elements (a), mesh deformation (b), elastic recovery of selected asperity nodes with variable geometry (A) and constant (B).
On the basis of the obtained results, an Equation (4) was statistically developed, describing the influence of height deviation and distance deviation on the elastic recovery of the material. Multidimensional correlation coefficient R = 0.8332.
Δ u ¯ ^ 0 = 0.0048 0.0032 Δ s t 0.026 Δ h t 3.526 × 10 14 Δ s t Δ h t + 0.024 Δ s t 2 + 0.358 Δ h t 2
As in the case of experimental studies, it was found that an increase of the asperities’ height deviation causes a decrease in the value of elastic recovery, while the opposite situation occurs when the distance deviation is analyzed. The results of experimental studies and computer simulations differ slightly from each other. This may result from simplifications and assumptions made when developing the computer model (e.g., refinement of the finite element mesh, shape of elements). Based on Equation (4), a spatial diagram (Figure 11) was prepared.

4. Researches of the Phenomenon of Elastic Recovery Depending on the Vertical Angle of the Asperities

4.1. Modelling Researches of the Phenomenon of Elastic Recovery

Experimental studies with the use of a model material were carried out in order to assess the effect of the flow blocking phenomenon through adjacent irregularities on the elastic return of the surface after burnishing. The tests were carried out in accordance with the rheological, geometric, and time scales [4,36,44]. Modeling consists in replacing a real research object with a model object. The real and model objects form a system of equivalent objects whose similarity is determined by mathematical relationships. Then, the modeling results can be translated into a real research object. Physical modeling on substitute material is used in experimental qualitative and quantitative analysis. The most frequently used model material replacing metal is a non-metallic material, such as plasticine, wax, putty, etc. These are materials with much lower plastic flow resistance. Physical modeling on plasticine is used for experimental qualitative and quantitative analysis of the plastic working processes. Qualitative analysis mainly concerns the kinematics of the material flow process, while the physical analysis covers the study of force and energy parameters (force, unit pressure on the contact surface of the tool—model material). In recent years, the technique of visualizing the flow trajectory of material particles of the shaped material, the so-called method of visioplasticity, has been used. It is based on the study of plastic flow kinematics based on the observation and measurements of the coordination grid. Performing experimental tests of visualization on real metal materials is laborious and costly and requires a station with very high-pressure forces. Therefore, often in experiments, the real material is replaced by the so-called model material or numerical analyses are performed [39,42,45]. In order to ensure the condition of rheological similarity, the characteristics for the model material—plasticine—were determined. This is due to its availability, low price, easy workability, plasticity (plasticizing stresses are 100–1000 times lower than the appropriate metal stresses), the possibility of its modification and regeneration, and therefore, multiple use. The material model was developed based on the static compression test. The results of this test are approximate data due to the presence of disturbances in the measurement of the compression force. The factor disturbing the measurement of the compressive force is the friction occurring on the faces of the cylindrical sample.
Under ideal conditions, during a static compression test, the diameter of a cylindrical specimen would change uniformly over its entire height. However, for real conditions, the effect of friction causes the sample material to flow unevenly (a barrel is formed), which is a result of blocking the flow of material on the front contact surfaces of the sample with the substrate and the punch. The static tensile test does not have this disadvantage, however, due to the nature of the material, it is impossible to carry out.
The compression test was carried out on the stand shown in Figure 12a. A cylindrical sample with the initial height h0 = 50 mm and the initial diameter d0 = 2r0 = 40 mm was used for the tests. On the basis of the compression test, the dependence of the force as a function of the change of sample height P = fh) was developed (Figure 12b). Using the results of the compression test and appropriate dependencies from the works of Dawidenkow and Spiridonowa [39,43,46], a material model was developed in the form of an actual compression diagram (Figure 12c). The obtained material characteristics for the model material are similar to the material characteristics of C 45 steel.
The aim of the model research was to determine the relationship between the vertical angle of triangular asperities and the elastic recovery of the model material. The burnishing element in the test was a flat punch mounted on a hydraulic press. The tests allowed for visual assessment of material displacement, numerical determination of the size of deformation, and elastic recovery of the material. It is also possible to determine the magnitude of stresses and the direction of material flow. The samples made of model material (plasticine) were characterized by the vertical angles of asperities from θ = 60° to θ = 150°, every 15°. Each sample had three asperities with the same vertical angles and was cut along the longer sides to apply a coordination grid. Photos of samples before, during, and after the burnishing process are presented in Table 2.
Then, the samples were placed in the matrix, which task was to properly position the sample in relation to the punch and to ensure appropriate boundary and initial conditions. The stamp was displaced, causing the deformation of the vertices of the triangular asperities and the deformation of the coordination grid. The sample and stamp were covered with a layer of talcum powder to reduce friction between them. Using a coordination grid applied to the samples, it was possible to perform appropriate measurements of the elastic recovery of the material.
The obtained results are presented in the form of a graph in Figure 13. The model tests of the elastic recovery of the model material using the visioplasticity method show that with the increase of the value of the vertical angle of the surface asperities, the value of the elastic recovery of the material decreases. The results were confirmed by numerical analysis. In some areas, the results are very close. The differences may be due to simplifications and assumptions made during creating a computer model. The accuracy of FEM calculations is influenced, among others, by the ratio of the side length of the finite elements used to discretize the object (object and tool), mesh refinement (mainly in areas of strong geometric or physical nonlinearity), and the finite element shape function. Although, in order to determine the significance of the influence of these factors and to determine an effective discrete model for the case of roller burnishing of a rough surface, sensitivity analysis was performed. As a result, acceptable gaps between the results were obtained.

4.2. Computer Simulation of the Phenomenon of Elastic Recovery Depending on the Vertical Angle of the Asperities

Computer simulations were developed for the same values of the asperities’ peak angles as in the research on the model material. The calculations were made in two stages. In the first stage, triangular irregularities were burnished, and in the second stage, the burnishing tool was moved away. It was burnished to a depth of a = 0.5 h, where h is the height of the triangular asperity. This approach to the problem made it possible to calculate the value of the material elastic recovery by determining the displacement of the nodal point of the peaks of the asperity. The difference between the node displacement modules at the moment of contact with the burnishing tool and after the tool is fully retracted is the desired result (Figure 14). Table 3 presents selected results from 2D computer simulations before, during, and after the burnishing process, showing displacement maps.
Figure 14 and Figure 15 show diagrams of the displacement of the node lying on the contact surface (the highest point of asperity) along the Y axis.
In the initial phase of the process, the node was slowly displaced due to the impact of the tool on the burnished material. The asperity was deforming. At the moment when the burnishing tool shifted to half of the roughness, the observed node reached the highest value of displacement. The first stage of numerical analysis has been completed. In the second step, when the time exceeded 1 s, the crushing tool was moved back.
The phenomenon of elastic recovery of the material is visible and then the stabilization of the position of the node (A) (the peak of the roughness), which has stopped moving along the Y axis. Similar results were obtained for the cases of crushing the remaining asperities with different vertical angles. When analyzing the node located between the asperities (B) (on the bottom), no elastic recovery was found for the case of asperities with a constant vertical angle and height.
The difference between the individual simulation results concerns the displacement value, which is related to the burnishing depth. The burnishing depth resulted from the geometry of the outline, i.e., the vertical angle of the asperities, and thus its height. The distance of the asperity was the same in each case. The elastic recovery also changed its value—the greater the vertical angle of the roughness, the smaller the value of the elastic recovery.
Computer simulations made it possible to observe and quantify the amount of material elastic recovery during the burnishing process. It is possible to check the results at any time during the process. It was confirmed that with the increase of the value of the vertical angle of the surface roughness, the value of the elastic recovery of the material is smaller. Moreover, using the developed numerical model, it is possible to qualitatively and quantitatively determine how much the bottoms between surface asperities rise in the process of their crushing. In cases where the vertical angles of the asperities were equal to θ = 60° and θ = 75°, the indentations of the asperities (bottom) did not rise, and the material core remained undeformed, and with complete deformation, deformed asperities are visible, separated from each other by gaps (planes of discontinuity) at depth a = 0.5 Rt. In the range of vertical angle θ = 90–135°, the plastic deformation zone increased to encompass the material core. The indentations between the asperities (bottom) rise (Figure 15), and with complete deformation, gaps are visible, but of less depth. For the angle θ = 150°, deformations of the roughness and the core of the material occur, and the cavities between the asperities (bottoms) rise completely.

5. Conclusions

The phenomenon of elastic recovery of the material during burnishing is a very complicated issue. The presented results of model tests, simulations, and analyses do not reflect the complexity of the problem. They confirm that this phenomenon occurs and there are methods by which it is possible to try to quantify the dependence of the elastic recovery on, for example, the apex angle of the inequality or deviations of the inequality outline. In the conducted analyses, only the elastic recovery of the burnished material was taken into account. The spring-return occurring in the tool material was omitted, treating the tool as perfectly rigid. It should also be noted that the phenomenon of elastic recovery probably also takes place during the treatments preceding burnishing. It is not possible to continuously measure the elastic recovery during modeling or experimental studies. This possibility is provided by computer simulations using the finite element method. However, the history of the material was not taken into account in the simulations performed. Model tests and computer simulations are a cheaper way than carrying out experimental tests each time. They give a guideline to consider the material elastic recovery when designing the burnishing process, especially in the case of fine machining. Then, when calculating the burnishing depth, its value should be increased by the elastic recovery of the material. Based on the research, the following conclusions can be drawn:
  • Knowledge about influence of determined profile parameters on the physical phenomena that occur in the burnishing process is unknown. It is important because they impact on the quality of the formed product’s surface layer, tribological properties, and fatigue issues. Many physical phenomena whose observation or measurement is very difficult or impossible take place during processes of burnishing, such as the pressure in contact zones, the friction forces, strain and stresses, elastic recovery, the slip and contact zones, and states of displacement.
  • After passing the burnishing tool, visco-plastic deformations remain in the workpiece, while the elastic deformations change and settle at the level corresponding to the new equilibrium state of the object. Thus, the states of displacements, stresses, and strains in the tested object during and after the end of the process differ significantly, because after the process is completed, the nodes are displaced in the elastic range. In the case of higher asperities, the displacement is greater than in the case of those with a height deviation. The opposite is the case for a distance deviation. These issues are very complex and require further analysis.
  • The obtained regression equations of the dependence of the material elastic recovery after burnishing (Δu0) from the deviations of height (Δht) and distance (Δst) of asperities after turning are reasonable for specific angles of asperities (θ = 60°) for C45 steel.
  • The increase of the deviation of the asperities’ height causes a decrease in the value of elastic recovery, while the opposite situation occurs when the deviation of the distance was analyzed.
  • Comparing the values of elastic recovery for an asperity with constant geometry, A1, and changeable geometry, A2, it can be seen that the value of elastic recovery was higher for asperity A1.
  • It was confirmed that with the increase of the value of the vertical angle of the surface roughness, the value of the elastic recovery of the material is smaller. Moreover, using the developed numerical model, it is possible to qualitatively and quantitatively determine how much the bottoms between surface asperities rise in the process of their crushing.

Author Contributions

The individual contributions of the authors: conceptualization, A.K.; methodology, A.K. and Ł.B.; investigation, A.K. and Ł.B.; software, A.K.; data curation, Ł.B.; validation, A.K.; writing—original draft preparation, A.K. and Ł.B.; writing—review and editing, A.K. and Ł.B.; visualization, Ł.B.; supervision, A.K. All authors have read and agreed to the published version of the manuscript.

Funding

The authors gratefully acknowledge the financial support from the Koszalin University of Technology, funding No. 524.02.23, entitled “Analysis of the deformation state of burnished surface with taking into account surface layer state and nonlinearity of the process”.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Δϕ(E)increment of logarithmic strain in elastic stage
ϕ(E)logarithmic strain in elastic stage
Δuielastic recovery
θvertical angle
σpyield stress

References

  1. Luca, L.; Neagu-Ventzel, S.; Marinescu, J. Effects of working parameters on surface finish in ball-burnishing of hardened steels. Precis. Eng. 2005, 29, 253–256. [Google Scholar] [CrossRef]
  2. El-Axir, M.H. An investigation into roller burnishing. Int. J. Mach. Tools Manuf. 2000, 40, 1603–1617. [Google Scholar] [CrossRef]
  3. Hassan, A.M. The effects of ball- and roller-burnishing on the surface roughness and hardness of some non-ferrous metals. J. Mater. Process. Technol. 1997, 72, 385–391. [Google Scholar] [CrossRef]
  4. Kułakowska, A. Problems of surface preparation under burnishing rolling in aspect of product quality. Steel Res. Int. 2010, 81, 218–221. [Google Scholar]
  5. Patyk, R. Theoretical and experimental basis of regular asperities about triangular outline embossing technology. Steel Res. Int. 2010, 81, 190–193. [Google Scholar]
  6. Przybylski, W. Technologia Obróbki Nagniataniem; WNT: Warszawa, Poland, 1987. (In Polish) [Google Scholar]
  7. Kulakowska, A.; Kukielka, L. Numerical analysis and experimental researches of burnishing rolling process with taking into account deviations in the surface asperities outline after previous treatment. Steel Res. Int. 2008, 2, 42–48. [Google Scholar]
  8. Patyk, R.; Kukielka, L. Optimization of geometrical parameters of regular triangular asperities of surface put to smooth burnishing. Steel Res. Int. 2008, 2, 642–647. [Google Scholar]
  9. Kulakowska, A. Experimental researches of burnishing rolling process of regular surface asperities prepared in turning process. In Proceedings of the 14th International Conference on Metal Forming, Krakow, Poland, 16–19 September 2012; pp. 127–131. [Google Scholar]
  10. Kukielka, L.; Kukielka, K. Numerical analysis of the physical phenomena in the working zone in the rolling process of the round thread. WIT Trans. Eng. Sci. 2007, 55, 125–134. [Google Scholar]
  11. Patyk, R. New Method of Technological Surface Layer Shaping of Machine Parts during Duplex Burnishing Rolling Process. In Proceedings of the 14th International Conference on Metal Forming, Krakow, Poland, 16–19 September 2012; pp. 115–118. [Google Scholar]
  12. Kułakowska, A.; Kukiełka, L. Analiza Numeryczna Wpływu Odchyłek Zarysu Nierówności Powierzchni po Toczeniu na Odchyłki Zarysu Nierówności Powierzchni po Nagniataniu Powierzchniowym; Projektowanie Procesów Technologicznych; Komisja Budowy Maszyn PAN: Poznań, Poland, 2006; pp. 171–178. ISBN 978-83-903808-7. (In Polish) [Google Scholar]
  13. Li, F.L.; Xia, W.; Zhou, Z.Y.; Zhao, J.; Tang, Z.Q. Analytical prediction and experimental verification of surface roughness during the burnishing process. Int. J. Mach. Tools Manuf. 2012, 62, 67–75. [Google Scholar] [CrossRef]
  14. Bucaille, J.L.; Felder, E.; Hochstetter, G. Mechanical analysis of the scratch test on elastic and perfectly plastic materials with the three-dimensional finite element modeling. Wear 2001, 249, 422–432. [Google Scholar] [CrossRef]
  15. Jardret, V.; Zahouani, H.; Loubet, J.L.; Mathia, T.G. Understanding and quantification of elastic and plastic deformation during a scratch test. Wear 1998, 218, 8–14. [Google Scholar] [CrossRef]
  16. Meza, J.H.; Moré Farias, M.-C.; Martinz de Souza, R.; Riaño, L.J.C. Using the ratio, maximum load over unload stiffness squared, Pm/Su2, on the evaluation of machine stiffness and area function of blunt indenters on depth-sensing indentation equipment. Mater. Res. 2007, 10, 437–447. [Google Scholar] [CrossRef] [Green Version]
  17. Pelletier, H. Predictive model to estimate the stress-strain curves of bulk metals using nanoindentation. Tribol. Int. 2006, 39, 593–606. [Google Scholar] [CrossRef]
  18. Stilwell, A.; Tabor, D. Elastic Recovery of Conical Indentations. Proc. Phys. Soc. 1961, 78. [Google Scholar] [CrossRef]
  19. Zajkowska, A.; Derpeński, Ł.; Seweryn, A. Metody wyznaczania rzeczywistej krzywej umocnienia materiału, Cz. II. Test wciskania kulistego wgłębnika. Acta Mech. Autom. 2011, 5, 143–147. (In Polish) [Google Scholar]
  20. Lee, J.D. A large-strain elastic-plastic finite element analysis of rolling process. Comput. Met. Appl. Mech. Eng. 1998, 161, 315–347. [Google Scholar] [CrossRef]
  21. Zang, S.L.; Liang, J.; Guo, C. A constitutive model for spring-back prediction in which the change of Young’s modulus with plastic deformation is considered. Int. J. Mach. Tools Manuf. 2007, 47, 1791–1797. [Google Scholar] [CrossRef]
  22. Chodór, J. Modelowanie i Analiza Przemieszczania Materiału Obrabianego Podczas Mikroskrawania i Nagniatania Ślizgowego z Uwzględnieniem Nieliniowości. Ph.D. Thesis, Politechnika Koszalińska, Koszalin, Poland, 2011. (In Polish). [Google Scholar]
  23. Skalski, K. Analiza Zagadnienia Kontaktowego Ciała Sprężysto-Plastycznego (na Podstawie Zasad Wariacyjnych i Metod Elementów Skończonych); Wydawnictwa Politechniki Warszawskiej: Warszawa, Poland, 1979. (In Polish) [Google Scholar]
  24. Gauthier, C.; Lafaye, S.; Schirrer, R. Elastic recovery of a scratch in a polymeric surface: Experiments and analysis. Tribol. Int. 2001, 34, 469–479. [Google Scholar] [CrossRef]
  25. Yang, H.; Guo, L.; Zhan, M.; Sun, Z. Research on the influence of material properties on cold ring rolling processes by 3D-FE numerical simulation. J. Mater. Process. Technol. 2006, 177, 634–638. [Google Scholar] [CrossRef]
  26. Schulze, V.; Bleicher, F.; Groche, P.; Guo, Y.B.; Pyun, Y.S. Surface modification by machine hammer peening and burnishing. CIRP Ann. 2016, 65, 809–832. [Google Scholar] [CrossRef]
  27. Janczewski, Ł.; Toboła, D.; Brostow, W.; Czechowski, K.; Kot, M.; Zagórski, K. Effects of ball burnishing on surface properties of low density polyethylene. Tribol. Int. 2016, 93, 36–42. [Google Scholar] [CrossRef]
  28. Nestler, A.; Schubert, A. Roller burnishing of particle reinforced aluminium matrix composites. Metals 2018, 8, 95. [Google Scholar] [CrossRef] [Green Version]
  29. Sánchez Egea, A.J.; Rodríguez, A.; Celentano, D.; Calleja, A.; López de Lacalle, L.N. Joining metrics enhancement when combining FSW and ball-burnishing in a 2050 aluminium alloy. Surf. Coat. Technol. 2019, 367, 327–335. [Google Scholar] [CrossRef] [Green Version]
  30. Juanjuan, H.; Wei, Z.; Guangchun, W. Investigation of influence factors on surface roughness of micro-scale features. Precis. Eng. 2019, 56, 524–552. [Google Scholar]
  31. Kuznetsov, P.; Tarasov, S.Y.; Dmitriev, A.I. Nanostructuring burnishing and subsurface shear instability. J. Mater. Process. Technol. 2015, 217, 327–335. [Google Scholar] [CrossRef]
  32. Denkena, B.; Grove, T.; Maiss, O. Surface texturing of rolling elements by hard ball-end milling and burnishing. Int. J. Adv. Manuf. Technol. 2017, 93, 3713–3721. [Google Scholar] [CrossRef]
  33. Amini, C.; Jerez-Mesa, R.; Travieso-Rodriguez, J.A.; Llumà, J.; Estevez-Urra, A. Finite element analysis of ball burnishing on ball-end milled surfaces considering their original topology and residual stress. Metals 2020, 10, 638. [Google Scholar] [CrossRef]
  34. Skoczylas, A.; Zaleski, K. Selected properties of the surface layer of C45 steel parts subjected to laser cutting and ball burnishing. Materials 2020, 13, 3429. [Google Scholar] [CrossRef]
  35. Llumà, J.; Gómez-Gras, G.; Jerez-Mesa, R.; Rue-Mascarell, J.; Travieso-Rodriguez, J.A. Mechanical strengthening in S235JR steel sheets through vibration-assisted ball burnishing. Metals 2020, 10, 1010. [Google Scholar] [CrossRef]
  36. Hill, R. The Mathematical Theory of Plasticity; Oxford University Press: Oxford, UK, 1986. [Google Scholar]
  37. Szczepiński, W. Mechanika Plastycznego Płynięcia; Małe Monografie; PWN: Warszawa, Poland, 1978. (In Polish) [Google Scholar]
  38. Johnson, W. Impact Strength of Materials; Edwards Arnolds: London, UK, 1972. [Google Scholar]
  39. Kukiełka, L. Teoretyczne i Doświadczalne Podstawy Powierzchniowego Nagniatania Tocznego z Elektrokontaktowym Nagrzewaniem; Monografia Wydziału Mechanicznego Wyższej Szkoły Inżynierskiej w Koszalinie no. 47; Wyższa Szkoła Inżynierska w Koszalinie: Koszalin, Poland, 1994. (In Polish) [Google Scholar]
  40. Kukiełka, K.; Kukiełka, L.; Olender, J.; Kuchta, R. Numerical Analysis of State of Deformation and Stress in Process of Trapezoidal Rolling; Zeszyty Naukowe Wydziału Mechanicznego no. 33; Politechnika Koszalińska: Koszalin, Poland, 2004; pp. 167–178. [Google Scholar]
  41. Kulakowska, A. Modelowanie, Analiza i Prognozowanie Cech i Wyników Procesu Nagniatania Naporowego Tocznego Warstwy Wierzchniej o Zdeterminowanej Strukturze Geometrycznej Powierzchni; Monografia no. 307; Politechnika Koszalińska: Koszalin, Poland, 2016. (In Polish) [Google Scholar]
  42. Kowalczyk, L. Modelowanie Fizykalne Procesów Obróbki Plastycznej; Instytutu Technologii Eksploatacji: Radom, Poland, 1995. (In Polish) [Google Scholar]
  43. Dawidenkow, N.N. Izbrannyje Trudy; Naukowa Dumka: Kijew, Ukraine, 1981; Volume 1. (In Ukrainian) [Google Scholar]
  44. Kukiełka, L. Podstawy Badań Inżynierskich; Politechnika Koszalińska: Koszalin, Poland; PWN: Warszawa, Poland, 2002. (In Polish) [Google Scholar]
  45. Patyk, R. Kształtowanie Regularnych Symetrycznych Nierówności Trójkątnych w Procesie Nagniatania Naporowego tocznego Wałków Stalowych. Ph.D. Thesis, Koszalin University of Technology, Koszalin, Poland, 2006. [Google Scholar]
  46. Kukiełka, L.; Geleta, K.; Kukiełka, K. Modelling of initial and boundary problems with geometrical and physical nonlinearity and its application in burnishing processes. In Proceedings of the 14th International Conference Metal Forming, Krakow, Poland, 16–19 September 2012; Wiley-VCH Verlag GmbH & Co.KGaA: Weinheim, Germany, 2012; pp. 1375–1378, ISBN 978-3-514-00797-0. [Google Scholar]
  47. Kukiełka, L. New damping models of metallic materials and its application in non-linear dynamical cold processes of metal forming. In Proceedings of the 13th International Conference Metal Forming 2010, Toyohashi, Japan, 19–22 September 2010; Verlag Stahleisen GmbH: Düsseldorf, Germany, 2010; pp. 1482–1485, ISBN 978-3-514-00774-1. [Google Scholar]
  48. Kukiełka, L.; Geleta, K.; Kukiełka, K. Modelling and analysis of nonlinear physical phenomena in the burnishing rolling operation with electrical current. In Proceedings of the 14th International Conference Metal Forming, Krakow, Poland, 16–19 September 2012; Wiley-VCH Verlag GmbH & Co.KGaA: Weinheim, Germany, 2012; pp. 1379–1382, ISBN 978-3-514-00797-0. [Google Scholar]
  49. Kukiełka, L.; Kukiełka, S. Zastosowanie Programu Experiment Planner 1.0 do Wspomagania Eksperymentu w Technologii Nagniatania, Obróbka Kształtująca i Powierzchniowa Nagniataniem; Zeszyty Naukowe Wydziału Mechanicznego nr 34; Wydawnictwo Uczelniane Politechniki Koszalińskiej: Koszalin, Poland, 2004; pp. 131–140. (In Polish) [Google Scholar]
  50. Kukielka, L. Mathematical modelling and numerical simulation of non-linear deformation of the asperity in the burnishing cold rolling operation. In Computational Methods in Contact Mechanics V, Proceedings of the 5th International Conference on Computational Methods in Contact Mechanics, Seville, Spain, 7 June 2001; Book Series: Computational and Experimental Methods; WIT Press: Ashurst, UK, 2001; pp. 317–326. [Google Scholar]
Figure 1. The influence of the vertical angle of a symmetrical triangle asperity on the depth of strain zone and surface roughness profile after burnishing: θ ≤ 80° (a), 80° < θ < 145° (b), θ ≥ 145° (c) [39].
Figure 1. The influence of the vertical angle of a symmetrical triangle asperity on the depth of strain zone and surface roughness profile after burnishing: θ ≤ 80° (a), 80° < θ < 145° (b), θ ≥ 145° (c) [39].
Materials 13 05276 g001aMaterials 13 05276 g001b
Figure 2. Elastic recovery during the roller burnishing process.
Figure 2. Elastic recovery during the roller burnishing process.
Materials 13 05276 g002
Figure 3. Schematic diagram of the elastic recovery of the material after burnishing: A1—first asperity, A2—second asperity after previous machining (turning), h—asperity’s height, Δht—height deviation of asperity after turning, Δhb—asperity’s height deviation after burnishing, ab—burnishing depth, Δu01—elastic recovery of the material of first asperity A1, Δu02—elastic recovery of the material of second asperity A2.
Figure 3. Schematic diagram of the elastic recovery of the material after burnishing: A1—first asperity, A2—second asperity after previous machining (turning), h—asperity’s height, Δht—height deviation of asperity after turning, Δhb—asperity’s height deviation after burnishing, ab—burnishing depth, Δu01—elastic recovery of the material of first asperity A1, Δu02—elastic recovery of the material of second asperity A2.
Materials 13 05276 g003
Figure 4. General view of the shaft with the samples after turning process and zoom of the sample (nb 3) after turning process, where Rt is the asperity height.
Figure 4. General view of the shaft with the samples after turning process and zoom of the sample (nb 3) after turning process, where Rt is the asperity height.
Materials 13 05276 g004
Figure 5. Burnishing position (a), general view of the shaft with the samples after burnishing (b), zoom of the asperities after roller burnishing (c).
Figure 5. Burnishing position (a), general view of the shaft with the samples after burnishing (b), zoom of the asperities after roller burnishing (c).
Materials 13 05276 g005
Figure 6. Exemplary profile of the asperities after roller burnishing.
Figure 6. Exemplary profile of the asperities after roller burnishing.
Materials 13 05276 g006
Figure 7. Dependence of the material elastic recovery after burnishing (Δu0) from the deviations of height (Δht) and distance (Δst).
Figure 7. Dependence of the material elastic recovery after burnishing (Δu0) from the deviations of height (Δht) and distance (Δst).
Materials 13 05276 g007
Figure 8. Geometry of the object, tool, and finite element mesh.
Figure 8. Geometry of the object, tool, and finite element mesh.
Materials 13 05276 g008
Figure 9. The graphs of dependencies of elastic recovery Δu1 and Δu2 accordingly in the asperities A1 and A2, from height deviation, Δht (a), and dependencies of elastic recovery Δu1 and Δu2 accordingly in the asperity A1 and A2 from distance deviation, Δst (b).
Figure 9. The graphs of dependencies of elastic recovery Δu1 and Δu2 accordingly in the asperities A1 and A2, from height deviation, Δht (a), and dependencies of elastic recovery Δu1 and Δu2 accordingly in the asperity A1 and A2 from distance deviation, Δst (b).
Materials 13 05276 g009
Figure 10. Geometric model with a mesh of finite elements and mesh deformation (a), elastic return of selected inequality nodes with variable geometry—A and constant—B (b).
Figure 10. Geometric model with a mesh of finite elements and mesh deformation (a), elastic return of selected inequality nodes with variable geometry—A and constant—B (b).
Materials 13 05276 g010
Figure 11. Relation of elastic material recovery after burnishing—Δu0 = fht, Δst).
Figure 11. Relation of elastic material recovery after burnishing—Δu0 = fht, Δst).
Materials 13 05276 g011
Figure 12. Stand for static compression test on the model material: 1—bow dynamometer, 2—dial displacement sensor, 3—compressed sample, and 4—dial gauge calibrated in force units (a), change of the compressive force P as a function of the height change Δh of the cylindrical sample (b), the actual compression diagram for the model material (c).
Figure 12. Stand for static compression test on the model material: 1—bow dynamometer, 2—dial displacement sensor, 3—compressed sample, and 4—dial gauge calibrated in force units (a), change of the compressive force P as a function of the height change Δh of the cylindrical sample (b), the actual compression diagram for the model material (c).
Materials 13 05276 g012
Figure 13. Influence of the vertical angle of surface asperities, θ, on the material elastic recovery Δui0.
Figure 13. Influence of the vertical angle of surface asperities, θ, on the material elastic recovery Δui0.
Materials 13 05276 g013
Figure 14. Diagram of displacement along the Y axis of selected nodes of asperity with a vertical angle in the process of crushing with a flat punch: A—top of asperity, B—bottom between asperities.
Figure 14. Diagram of displacement along the Y axis of selected nodes of asperity with a vertical angle in the process of crushing with a flat punch: A—top of asperity, B—bottom between asperities.
Materials 13 05276 g014
Figure 15. Diagram of displacement along the Y axis of selected asperity nodes with the vertical angle θ = 75° (a), θ = 90° (b), θ = 105° (c), θ = 120° (d), θ = 135° (e), θ = 150° (f), crushed with a flat punch: A—peak of asperity, B—bottom between asperities.
Figure 15. Diagram of displacement along the Y axis of selected asperity nodes with the vertical angle θ = 75° (a), θ = 90° (b), θ = 105° (c), θ = 120° (d), θ = 135° (e), θ = 150° (f), crushed with a flat punch: A—peak of asperity, B—bottom between asperities.
Materials 13 05276 g015
Table 1. Experimental plan with real and coded values.
Table 1. Experimental plan with real and coded values.
12345678910111213
Δht0.010.010.030.030.020.020.0400.020.020.020.020.02
Δst0.050.150.050.150.200.10.10.10.10.10.10.1
Δht−1−1+1+100+1.414−1.41400000
Δst−1+1−1+1+1.414−1.4140000000
Table 2. View of the exemplary samples.
Table 2. View of the exemplary samples.
Before BurnishingDuring ProcessAfter Burnishing
θ = 60° Materials 13 05276 i001 Materials 13 05276 i002 Materials 13 05276 i003
θ = 150° Materials 13 05276 i004 Materials 13 05276 i005 Materials 13 05276 i006
Table 3. Results of sample computer simulations.
Table 3. Results of sample computer simulations.
Before BurnishingDuring ProcessAfter Burnishing
θ = 60° Materials 13 05276 i007 Materials 13 05276 i008 Materials 13 05276 i009
θ = 150° Materials 13 05276 i010 Materials 13 05276 i011 Materials 13 05276 i012
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kułakowska, A.; Bohdal, Ł. Researches and Simulation of Elastic Recovery Phenomena during Roller Burnishing Process of Macro-Asperities of Surface. Materials 2020, 13, 5276. https://doi.org/10.3390/ma13225276

AMA Style

Kułakowska A, Bohdal Ł. Researches and Simulation of Elastic Recovery Phenomena during Roller Burnishing Process of Macro-Asperities of Surface. Materials. 2020; 13(22):5276. https://doi.org/10.3390/ma13225276

Chicago/Turabian Style

Kułakowska, Agnieszka, and Łukasz Bohdal. 2020. "Researches and Simulation of Elastic Recovery Phenomena during Roller Burnishing Process of Macro-Asperities of Surface" Materials 13, no. 22: 5276. https://doi.org/10.3390/ma13225276

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop