Next Article in Journal
New Concept for the Facile Fabrication of Core–Shell CuO@CuFe2O4 Photocathodes for PEC Application
Previous Article in Journal
Dimensionality of the Superconductivity in the Transition Metal Pnictide WP
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigation of the Robotized Incremental Metal-Sheet Forming Process with Ultrasonic Excitation

by
Vytautas Ostasevicius
1,*,
Agne Paulauskaite-Taraseviciene
2,
Ieva Paleviciute
1,
Vytautas Jurenas
1,
Paulius Griskevicius
3,
Darius Eidukynas
1 and
Laura Kizauskiene
4
1
Institute of Mechatronics, Kaunas University of Technology, Studentu Street 56, LT-51424 Kaunas, Lithuania
2
Department of Applied Informatics, Kaunas University of Technology, Studentu Street 50, LT-51368 Kaunas, Lithuania
3
Department of Mechanical Engineering, Kaunas University of Technology, Studentu Street 56, LT-51424 Kaunas, Lithuania
4
Department of Computer Sciences, Kaunas University of Technology, Studentu Street 50, LT-51368 Kaunas, Lithuania
*
Author to whom correspondence should be addressed.
Materials 2022, 15(3), 1024; https://doi.org/10.3390/ma15031024
Submission received: 16 December 2021 / Revised: 21 January 2022 / Accepted: 24 January 2022 / Published: 28 January 2022

Abstract

:
During the single-point incremental forming (SPIF) process, a sheet is formed by a locally acting stress field on the surface consisting of a normal and shear component that is strongly affected by friction of the dragging forming tool. SPIF is usually performed under well-lubricated conditions in order to reduce friction. Instead of lubricating the contact surface of the sheet metal, we propose an innovative, environmentally friendly method to reduce the coefficient of friction by ultrasonic excitation of the metal sheet. By evaluating the tool-workpiece interaction process as non-linear due to large deformations in the metal sheet, the finite element method (FEM) allows for a virtual evaluation of the deformation and piercing parameters of the SPIF process in order to determine destructive loads.

1. Introduction

Robotized single-point incremental forming (SPIF), whereby a small-sized tool deforms a sheet of metal, is ideal for replacing expensive stamping processes, assuring flexibility and cost-effectiveness, as well as increasing the competitiveness of companies. This technology opens up new possibilities compared to traditional small-series production methods wherein added value is determined by skill and qualification. The forces resulting from friction between the tool and the workpiece during machining play an important role in product quality. The influence of petroleum and vegetable-oil-based greases on the coefficient of friction, abrasion, forming forces and surface roughness of metal sheets produced by SPIF was investigated in [1]. Lubricating oil was found to produce a surface roughness close to Ra = 1.45 μm, which exceeds the surface roughness of an undeformed sheet and degrades the quality of the final product. As lubricants are not environmentally friendly, new methods for reducing and predicting forming forces need to be explored.
Various research methods have been used to evaluate the mechanical features of the robotized SPIF process [2,3,4]. The common objective of research to date [2] has been to analyse the results of SPIF processes employing industrial robots with respect to the potential increase in technological capabilities (five-axis production) and process versatility (ability to produce additional flanges or complete work) compared to traditional three-axis machining. In [3,4] Belchior et al. present an approach to link finite element method (FEM) analysis with a model of the elastic robot structure to improve the geometric accuracy of the formed metal sheet.
The complexity of SPIF mechanics and special conditions, such as bending under stress and cyclic bending/unbending, as well as the shear deformation process, which contribute to the improvement of overall formability, require innovative research efforts and techniques. Major parameters affecting formability are tool size, step down (drawing angle) and thickness and material properties of the metal sheet [5,6,7]. Formability of the metal sheet increases with decreasing tool size, step down per revolution and thickness of the sheet, while the feed rate does not have a significant effect on formability [5]. The SPIF process has much higher forming limits than those achievable by conventional sheet-forming processes, such as deep drawing, which is affected by complicated stress state and strain path. The strain-based forming limit curve criterion is widely used in the sheet-metal forming industry to predict fracture. This criterion is only valid when the strain path is linear [6]. In contrast, during the SPIF process, the strain path is rather nonlinear, and the practice of using a strain-based forming limit criterion often leads to an erroneous assessment of formability and failure prediction. Stress-based forming limits are not sensitive to strain-path changes and are used to model the necking limit, combined with the fracture limit, based on the criterion of maximum shear stress [7].
The behaviour of the ultrasonically excited forming tool during the SPIF process was investigated in [8,9]. An FEM model was employed to evaluate the influence of ultrasonic vibration frequencies and amplitudes on the accuracy of the technological process.
In order to reduce computational time in SPIF simulations, the FEM mesh is composed of several non-overlapping parts subjected to plastic and elastic deformations [10]. The plastic deformations are localized, including two substructures: plastic-non-linear and elastic-pseudo-linear. The group of plastic structures creates a finite-element mesh that is in contact with the tool used to form the material. Studies based on Marciniak–Kuczynski (MK) theory have shown that the MK damage criterion can be used to predict SPIF limits for arbitrary loading paths. An analytical process is required to investigate the extent to which such monoform strain paths can be induced by highly non-monotonic deformation pathways [11]. The MK model can be used to predict the onset of necking when strain paths have been obtained by FEM modelling. The effect of stress rate on the forming parameters of different types of steels was examined in [12]. The results revealed that the impacts of stress rate on the forces and energies used in the forming process and their limits are insignificant and highly material-dependent.
Materials such as aluminium alloys are widely used to manufacture automotive, aerospace and other components due to their good strength, low weight and formability. Therefore, the peculiarities of SPIF processes using these alloys have been extensively studied by researchers [13,14,15,16]. Based on the simplified proportional-load assumption, the fracture-boundary model, in contrast to a conventional neck, can represent the failure mode according to the AA1050-H111 failure mechanism [13]. The SPIF process applies elastic fracture criteria to the generalized tensile stress state based on the degree of triaxiality of the stresses. The research presented in [14] explains the systematic progress of SPIF formation and damage, and the extended Gurson–Tvergaard–Needleman model is used to predict surface damage. AlMgSc alloy, well known in the field of aeronautics, was investigated in [15].
Multiple experiments have been carried out wherein material laws and corresponding material variables were selected and adjusted to correctly define the behaviour of the material. In order to define the deformation behaviour of sheet metal in real production, it is essential to accurately determine the ultimate deformations of the SPIF process compared to traditional material-formation limits and the effect of the process parameters on these deformations [16]. In addition, the force of the SPIF process is very important, especially when designing special equipment or ensuring that all safety regulations are observed when using adapted machinery. Ambrogio et al. [17] examined SPIF failures by experimentally influencing process parameters, such as tool diameter, spindle rotation speed and step down on the SPIF (spifability) of AISI 304 metal sheets and analysed the results according to the circle-grid method.
A review of the scientific literature shows that no attempt has been made so far to reduce friction between the tool and sheet metal by removing environmentally unfriendly petroleum products, which also affect the chemical composition of the formed sheets. Another important point is that the lubricated surfaces need to be cleaned after the forming process, which also increases the development time and cost of the product.
SPIF is a process in which many independent parameters operate simultaneously, and it is an important task for engineers to anticipate these factors. A review of published research results suggests that alternative methods to reduce the coefficient of friction between the tool and the metal-sheet contact surfaces have not yet been proposed; the dynamics of the process have not been sufficiently investigated, and there is a lack of physical evidence suggesting that virtual models are appropriate.

2. Materials and Methods

2.1. Methodology for Evaluation of Forming Process

Metallic materials used in deforming processes, such as single-point incremental forming (SPIF), visibly deform on contact with the tool area. The key challenge of the SPIF process is to effectively assess the shaping forces. The forces generated during SPIF may be controlled by varying the coefficient of friction between the tool and the workpiece. Lubrication is an important factor in the SPIF process, reducing friction in the tool-workpiece contact area, but the use of grease is associated with environmental problems. Therefore, it is necessary to find other ways to reduce the shaping forces associated with the process dynamics. This requires, in particular, the development of a physically adequate mathematical model of the process, the mechanical parameters of which can be adjusted experimentally using the 3D scanning device shown in Figure 1.
A 300 × 300 × 0.5 mm aluminum alloy AW1050 sheet embedded in a 30 × 30 mm L angle profile welded steel frame was formed. The frequencies of the resonant modes were determined using the COMSOL Multiphysics FEM. In order to change the lubrication of the contact surface between the forming tool and the metal sheet, an attempt was made to excite the metal sheet by vibrations. For this purpose, equipment, the scheme of which is presented in Figure 2, was developed to excite 3D ultrasonic frequency vibrations in the metal sheet.
This equipment consists of an angled-profile metal frame (2) to which piezoceramic elements (3 and 4) are glued on different planes in order to excite bending oscillations of the higher harmonics in the frame. Since the metal sheet is rigidly attached to the frame by its outer contour, the vibrations of the frame also excite the metal sheet. This principle of ultrasonic 3D vibration excitation in the metal sheet makes it possible to reduce the friction between the metal sheet and the tool during SPIF in order to improve the lubrication conditions at the tool-workpiece contact pair and, in some cases, to completely eliminate the need for lubrication. The number of piezoceramic elements in the rectangular frame used to excite the 3D vibrations in the metal sheet can range from 2 to 8, according to the given scheme (Figure 2), by alimenting them with harmonic ultrasonic-frequency electrical signals. Experimental studies on a 0.5 mm thick aluminum plate showed that the maximum reduction in friction forces at the tool-workpiece contact pair was observed when the piezoelectric actuators were excited in the 25–35 kHz frequency range.

2.2. Mechanical Testing of Aluminum Alloy Sheets

A 0.5 mm thick sheet of aluminum alloy, EN AW1050A H24, was used for experimental and numerical analysis. Material properties, taken from the supplier datasheet, are presented in Table 1.
As failure strain is a parameter dependent on stress state [18] and the value obtained in the uniaxial test may differ from the failure strain in the SPIF process, simple cupping was selected for calibration of the material model. During the SPIF process, the strain path is nonlinear [19]; therefore, not only the uniaxial failure strain but also the conventional forming-limit diagram criterion can cause inaccurate failure prediction.
The stress-strain relationship needed for numerical analysis was obtained by a reverse-engineering approach using data from the physical cupping test. The physical cupping test was chosen to describe the material behaviour according to the power-law material model in order to reduce the number of material-model constants needed to calibrate. Using the range of material properties from the material datasheet (yield stresses, ultimate strength range and failure strain), the power-law parameters, such as strength coefficient (K), hardening exponent (n) and failure strain, were calibrated using the LS-Opt system. The input range of the stress-strain curves for the power-law models used in LS-Opt is presented in Figure 3, where the solid line represents the calibrated power-law material curve. The power-law equation is provided below.
σ = K ε n = 143 ε 0.097
where ε is uniaxial strain.

2.3. Cupping Test for Material Model Calibration

In order to develop a validated numerical model of the SPIF process, several sets of experimental tests were performed. Firstly, the material model was calibrated by reverse engineering using experimental results of the cupping test. Subsequently, the same material model was used for the SPIF simulation. In order to characterize the plastic behaviour and ductile fracture of the aluminum alloy sheet, the Erichsen cupping test was performed on square aluminum plates, in accordance with ISO 20482: 2013 (metallic materials; sheets and strip; Erichsen tensile test). The Erichsen standard provides information on fracture under the equi-biaxial state of stress. The results are influenced by sheet thickness and the friction between the sheet and the tool surface. For the Erichsen compression test, the specimen was attached to a 55 mm die with a heavy flange and bolts tightened to achieve a force of approximately 10 kN on the workpiece holder (Figure 4). The radius of the hemi-spherical punch was 10 mm, and the test speed was 5 mm/min.
The experiment was carried out by applying a hemispherical punch to a sheet of metal until a crack appeared.

2.4. Experimental Investigation of SPIF Aluminum Alloy Sheet

In order to find a more efficient method to control the forming force, an attempt was made to excite the aluminum-alloy sample with ultrasonic vibrations. The experimental setup created for this purpose is shown in Figure 5. An investigated sheet (2) was clamped to the frame (1), and it was excited with two piezoelectric transducers (Ferroperm Piezoceramics A/S, Kvistgaard, Denmark) (8). A power amplifier (7) was used to generate the vibrations in the piezoelectric transducer. The robotic arm (4) incrementally formed the sheet by moving the tool (3) with the sphere attached to it. Formation force was measured by a pressing-force sensor (5) located under the frame, and values were obtained in the display controller (6).

3. Results and Discussion

3.1. Analysis of Ultrasonic Sheet Vibrations

The metal-sheet vibrational analysis, performed with the Polytec scanning vibrometer, PSV-500-3D-HV, presented in Figure 1, shows that in this frequency band, the workpiece is dominated by planar (XY) higher harmonic vibrations, which are significantly less suppressed in the tool-sheet pair friction contact than in direction Z, perpendicular to the metal sheet (Figure 6). These results show that when two piezoelectric actuators in different planes are excited, in-plane (XY) oscillations in the 30–33 kHz band are greatly enhanced.
In order to assess the effect of vibrations on the SPIF of an aluminum-alloy sheet, of surface-roughness measurements of the sheet were taken with and without ultrasonic vibrations. An advanced surface-roughness tester, TR200 (Beijing TIME High Technology Ltd., China), was used to measure roughness. The measured surface roughness of the sheet formed without ultrasonic vibrations varied in the range of Ra = 0.30–0.33 μm, and with ultrasonic vibrations, surface roughness was in the range of Ra = 0.18–0.25 μm.
In order to determine the frictional force between the tool and the aluminum-alloy sheet (Figure 7a), calculations were carried out in accordance with the diagram in Figure 7b.
Since the roughness momentum gauge was placed with an angle of 30° degrees, both pressure and friction torque could be obtained. Friction force can be calculated as follows:
F T = k T F N
F T O = F T sin α
M O = F O L
F N = F O sin α
M O + M T = F O + F T O L
where L is the length of the tool link connecting the sphere and moment gauge (150 mm); α is an angle between the sample sheet surface and the tool link (30°); MO and MT are tool-pressing-force- and tool-friction-force-generated angular momentums, respectively; FO and FTO are tool pressure force and the friction force projection to tool link, respectively; FN and FT are the tool‘s sphere pressure and friction force, respectively; kT is the friction force coefficient between the sample sheet and the steel-tool sphere; and vs is the steel tool‘s sphere speed on the surface of the sample sheet (mm/s).
The mechanical friction between the steel tool and aluminum-alloy sheet was measured using a torque sensor, STJ100, connected with a BGI series digital force gauge (Mark-10 Corp., New York, NY, USA) and a PC (Figure 7). The torque sensor‘s sensitivity factor was 6 Nm/V, and speed rate of the tool was −1200 mm/min. The measurement results of the tool friction coefficient and friction force on dry, lubricated and ultrasonically excited surfaces in the 30–33 kHz frequency range, which reduces the friction force at the tool-sheet metal contact pair, are shown in Table 2.
Table 2 shows that the coefficient of friction between the steel tool and the aluminum-alloy sheet was close to that of the lubricated surfaces when subjected to ultrasonic vibrations. This makes it possible to solve environmental problems, and the surface of the manufactured part does not need to be cleaned, while at the same time, its roughness is reduced by half. This revealed phenomenon has been patented by the authors [20].

3.2. Results of Theoretical and Experimental Investigation of the Aluminium-Alloy Sheet

Strain path, stress state and prediction of the formability of the SPIF were evaluated by FE analysis. In all subsequent simulations, the explicit FE code LS-Dyna was used. The aluminum sheet was modelled using fully integrated shell elements with thickness stretch allowed. Through the thickness, linear variation of strain was evaluated by five integration points, and a shear correction factor of 0.833 was used. Contact between the punch/tool and the aluminum sheet was described using the keyword: * CONTACT_FORMING _ONE_WAY _SURFACE_TO_SURFACE. There were no damage models involved; only a simple failure model was used, based on maximum failure effective strain criteria.
To simulate the cupping test and calibrate the material model. The FE model consisted of four parts: a deformable aluminum sheet, a rigid spherical punch, a bottom holder and a top holder. The holders were compressed with a force of 10 kN. The size of the shell elements varied between 0.5 at the centre/contact zone and 1 mm on the external side of the specimen. The FEM model (Figure 8a) consisted of approximately 8500 shell elements. The speed of the punch in the numerical explicit model was 1 m/s, while in physical experiments, the punch speed was 5 mm/min. The elastoplastic properties of aluminum were described by the power-law plasticity material model, and three parameters were chosen for calibration: strength coefficient (K), strain-hardening coefficient (n) and failure strain (ε_u). The simulation results are presented in Figure 8b.
After curve fitting, the experimental force-displacement relation and the FE simulation prediction (Figure 8b) correlated well. The calibrated material model curve is presented in Figure 5. The spherical punch penetration into an aluminum-alloy sheet is illustrated in Figure 9. Comparing the fracture paths also revealed a common behaviour (Figure 6b and Figure 9). The most highly stressed point is located at the apex of the dome, where the components of radial and circumferential stress and strain are equal to one another. The fracture lines are diagonal and almost symmetrically located between the centre and the corners. Three diagonal fracture lines dominate in the physical test experiments, while four diagonal fracture lines were obtained in the simulations. One of the reasons for this could be the anisotropy of the material, which is not evaluated in the FE model.

3.3. Numerical Simulations of SPIF Process

The SPIF model consisted of two parts: an aluminum sheet and a forming tool (Figure 10). The aluminum sheet was modelled in the same way as for the cupping test. The size of the shell elements varied between 1 mm at the centre/contact zone and 5 mm on the external side of the specimen. The FE model consisted of about 20,000 shell elements. The main features of the FE model were copied from the calibrated cupping-test model. The elastoplastic properties of aluminum were described by a power-law plasticity-material model calibrated with a cupping test. The forming tool was modelled as a rigid body using shell elements. The aluminum-sheet holder was simulated with fixed nodes (25 mm wide around all four edges). Major diameter of the helix was 140 mm.
The load was applied by three functions of tool displacements: f(x), f(y) and f(z). All functions were continuous and varied so that the horizontal displacement per revolution was 0.5 mm and the vertical displacement per revolution was 0.5 mm. In an explicit numerical model, the speed of the forming tool was about 2 m/s; in physical experiments, it was 1.2 m/min. The trajectories of the forming-tool path are shown in Figure 11.
The experimentally obtained ultrasonic excitation effect on friction forces was evaluated numerically by applying different friction coefficients. Various simulations were carried out, with coefficients of friction ranging from 0 to 0.5. The effect of friction on the horizontal-force component is shown in Figure 12.
It can be seen that increasing the coefficient of friction from 0 to 0.5 increases the amplitudes of the X-force component by a factor of approximately two. In contrast, the effect on the vertical Z-force component is negligible (Figure 13a). The effect of friction on the resultant force is shown in Figure 13b.
SPIF is characterized by a reduction in the wall thickness of the final finished part compared to the initial thickness of the sheet metal. The excessive-thickness reduction ratio in the deformation zones when the sheet metal is formed separately has a significant effect on the forming limit. Prediction of the thickness of the deformation zone is an important approach to control the thinning ratio. With regard to the object of study, aluminum alloy, the principle of thickness deformation in the SPIF process is presented in Figure 14. The relationship between wall thickness and drawing angle (α) can be expressed by the sine law used in SPIF:
t f = t 0 sin 90 α
where t0 is the initial thickness; tf is the final thickness; and α is the drawing angle between the initial flat surface and the deformed surface.
The drawing angle obtained in the simulation varied from 47°; therefore, according to the sine law, this would lead a reduced thickness by up to 0.34 mm. In contrast, the minimum thickness obtained in the simulation was 0.335 mm. Three elements (A, B, C) with different drawing angles of 47°, 30° and 17°, respectively, after a forming depth of 25 mm show the development of the thickness reduction. The results obtained from the SPIF simulation correlate well with the sine law. Strain distribution after a forming depth of 25 mm is presented in Figure 15a. Shell element 6285 of the aluminum sheet comes into contact with the forming tool at approximately 35 revolutions. Figure 15b shows that a significant increase in effective strains starts after contact with the tool and increases with further revolutions of the tool.
Figure 14 and Figure 15 show that the decrease in sheet thickness during SPIF correlates with effective plastic strain changes. The thinning of the forming sheet is one of the main failure modes of SPIF and is related to the drawing angle (α), which is one of the geometrical limits of the SPIF process. Figure 15 shows that the effective strains achieved during incremental forming (maximum effective strain, 60%) are much higher than the values given in the material datasheet since tensile failure strain constitutes 3–8%, and the values of failure strain calibrated by the cupping test constitute 38%. This confirms that failure prediction in the incremental forming process is not so straightforward and that other failure criteria should be applied. Comparing the evolution of stress triaxialities during the cupping test and SPIF (Figure 16), it can be observed that during the cupping test, the stress triaxiality is stable and equal to the 2/3, which corresponds to the equi-biaxial tension stress state, while the stress triaxiality during the SPIF process varies from −0.6 to 0.6. This confirms that the stress state during the SPIF process varies from plane-strain compression up to plane-strain tension. During the SPIF process can be achieved higher formability because at the plane-strain tension state, the failure strain of the material has a lower value compared with failure strains at the other stress states.
Deformation phenomena during incremental formation are still under discussion. In the absence of direct experimental evidence, some researchers have claimed that deformation is caused by stretching, while the others state that it is caused by through-thickness shear [21,22]. If out-of-plane shear dominates, the principal strains are not in the plane [23]; therefore, the use of shell elements can be limited, as shell elements typically do not allow for capture of through-thickness shear properties in simulations. Fully integrated shell elements with a thickness stretch in LS-Dyna allow for the capture of through-thickness shear stresses and shear strains. Three elements (S2681, S2683 and S2685; see Figure 15a) along the x-axis were selected to analyse stresses and strains during simulation of the SPIF process. The strain distribution through the thickness is shown in Figure 17.
The results show that the through-thickness shear strains have significant values compared to the in-plane (x,y) shear strains. In particuar, the transverse through-thickness deformation of element No. 6285 accounts for about 50% of the in-plane shear deformation.
As vibrational excitation of the sheet-metal sample was performed by piezo actuators attached to the steel frame, the calculations were also performed by exciting the frame in which the aluminum-alloy sheet was mounted. The results of experimental and theoretical studies for the first embedded sheet in steel-frame mode, presented in Figure 18, show a good correlation between the dynamic properties of the COMSOL Multiphysics (81.707 Hz) and experimental (78.8 Hz) models.
The force-displacement response was obtained by plotting the reaction force on the punch vs. the displacement of the punch (Figure 19).
The maximum load for all five cupping tests was Fmax = 1160 N, Erichsen index IE = 5.36 mm. Force vs. displacement curve was chosen as the validation criterion for the numerical model.

3.4. Experimental Validation of SPIF Simulation Results

The SPIF experiment was performed under dry, lubricated and vibration-excited friction on the contact surfaces of the tool and the sheet metal. The vertical force dependences on the process conditions are given in Figure 20.
Accordingly, to verify modelling results (presented in Figure 8), the puncture forces of the sheet metal during experimental investigation are shown and Figure 21. It can be observed that at a depth of 0, the vertical force is equal to 224.04 N; at a depth of 10, the vertical force is 778.04 N; at a depth of 16, the vertical force is 1064 N; however, when the depth reaches a value of 15, the vertical force starts to decrease significantly.
Based on the obtained results, curve and magnitude of puncture force correlated well with modelling results; therefore, the theoretical model is verified.
In order to validate the adequacy of the simulated reduction in sheet thickness (Figure 14) for the experimentally obtained SPIF profile, the wall thickness of the formed aluminum specimen was measured at several points of the cross-sectional cut (Figure 22) using a Topex 31C629 micrometer with an accuracy of 0.01 mm.
The wall thickness of the incrementally formed aluminum sheet at measurement points is presented in Table 3.

3.5. Discussion

The SPIF experiment was performed under dry, lubricated and vibration-excited friction on the contact surfaces of the tool and the sheet metal. As depicted in Figure 20, depending on the forming depth, the forming force is the lowest under ultrasonic excitation. The results in Figure 6 show that when two piezoelectric actuators excite the aluminum-alloy sheet in different planes, the in-plane (XY) oscillations of the sheet in the 30–33 kHz ultrasonic frequency band are greatly amplified. Comparing the ultrasonic excitation frequencies obtained above with the simulated COMSOL Multiphysics finite elements and the experimentally obtained eigenfrequencies of the aluminum sheet in the direction perpendicular to its plane, which are 81,707 Hz and 78.8 Hz, respectively (Figure 18), these ultrasonic frequencies are several tens of times higher. The coefficient of friction of the applied concentrated contacts between the forming tool and the surface of the sheet is almost identical to the coefficient of friction between the vibrating sheet surface and the tool (Table 2). It was found that the surface roughness of the formed sheet without ultrasonic excitation varied in the range of Ra = 0.30–0.33 μm, whereas with ultrasonic excitation, it varied in the range of Ra = 0.18–0.25 μm. These SPIF properties of ultrasonic excitation could offer a better solution in terms of environmental protection, product cost and time to market compared to lubrication, which is widely used today. Measurements of the cross section of the formed sheet showed that the most deformed sheet was the thinnest at cross-sections points no 2 and 4 (Figure 22).

4. Conclusions

In this paper, numerical and experimental methods were used to analyse the deformation forces of single-point robotized incremental forming of an aluminum alloy sheet. We proposed an innovative method for reducing the frictional force between the forming tool and the sheet surface by excitation of the workpiece with high-frequency oscillations in two orthogonal directions in the sheet plane. The coefficient of friction between the dry sheet surface subjected to ultrasonic vibrations and the tool was found to be close to the coefficient of friction between the lubricated surfaces, which can reduce the product’s time to market and make the process more environmentally friendly. Thanks to the numerical FEM, the decreasing tendencies of the sheet cross section, depending on the degree of plastic deformation and the processes of sphere penetration into the aluminum-alloy sheet, were investigated and validated. The effective strains achieved during the incremental forming (maximum effective deformation, 60%) were found to be much higher than the values given in the material datasheet since tensile failure strain constitute 3–8%, and the values of failure strain calibrated by the cupping test constitute 38%. This confirms that failure prediction in the incremental forming process is not straightforward.

Author Contributions

V.O.: supervision, idea generation, investigation, funding acquisition, writing; A.P.-T.: experimental data exploration, review and editing; I.P.: investigation, data acquisition, formal analysis; V.J.: experimental research, methodology, resources, investigation; P.G.: modelling, theoretical research, visualization; D.E.: experimental research, data acquisition, visualization, review and editing; L.K.: writing, review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the European Regional Development Fund (project No 01.2.2-LMT-K-718-05-0076) under grant agreement with the Research Council of Lithuania (LMTLT). Funded as European Union’s measure in response to the COVID-19 pandemic.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zavala, J.M.D.; Martínez-Romero, O.; Elías-Zúñiga, A.; Gutiérrez, H.M.L.; Vega, A.E.; Taha-Tijerina, J. Study of friction and wear effects in aluminum parts anufactured via single point incremental forming process using petroleum and vegetable oil-based lubricants. Materials 2021, 14, 3973. [Google Scholar] [CrossRef] [PubMed]
  2. Callegari, M.; Amodio, D.; Ceretti, E.; Giardini, C. Sheet incremental forming: Advantages of robotized cells vs. CNC machines. In Industrial Robotics: Programming, Simulation and Applications; IntechOpen: London, UK, 2006; pp. 493–514. [Google Scholar]
  3. Belchior, J.; Leotoing, L.; Guines, D.; Courteille, E.; Maurine, P. A process/machine coupling approach: Application to robotized incremental sheet forming. J. Mater. Process. Technol. 2014, 214, 1605–1616. [Google Scholar] [CrossRef] [Green Version]
  4. Belchior, J.; Guillo, M.; Courteille, E.; Maurine, P.; Leotoing, L.; Guines, D. Off-line compensation of the tool path deviations on robotic machining: Application to incremental sheet forming. Robot. Comput. Integr. Manuf. 2013, 29, 58–69. [Google Scholar] [CrossRef] [Green Version]
  5. Bhattacharya, A.; Maneesh, K.; Venkata Reddy, N.; Cao, J. Formability and surface finish studies in single point incremental forming. J. Manuf. Sci. Eng. 2011, 133, 061020. [Google Scholar] [CrossRef]
  6. Stoughton, T.B.; Yoon, J.W. A new approach for failure criterion for sheet metals. Int. J. Plast. 2011, 27, 440–459. [Google Scholar] [CrossRef]
  7. Jawale, K.; Duarte, J.F.; Reis, A.; Silva, M.B. Microstructural investigation and lubrication study for single point incremental forming of copper. Int. J. Solids Struct. 2018, 151, 145–151. [Google Scholar] [CrossRef]
  8. Li, P.; He, J.; Liu, Q.; Yang, M.; Wang, Q.; Yuan, Q.; Li, Y. Evaluation of forming forces in ultrasonic incremental sheet metal forming. Aerosp. Sci. Technol. 2017, 63, 132–139. [Google Scholar] [CrossRef] [Green Version]
  9. Li, P.Y.; Liu, Q.; An, W.R.; Li, S.J. Finite element simulation of ultrasonic incremental forming. Mater. Sci. Forum 2016, 836, 452–461. [Google Scholar] [CrossRef]
  10. Hadoush, A.; Van den Boogaard, A.H. Substructuring in the implicit simulation of single point incremental sheet forming. Int. J. Mater. Form. 2009, 2, 181–189. [Google Scholar] [CrossRef] [Green Version]
  11. Mulay, A.; Ben, S.; Syed, I.; Ben, A. Artificial neural network modeling of quality prediction of a single point incremental sheet forming process. Adv. Sci. Technol. Lett. 2017, 147, 224–250. [Google Scholar]
  12. Verleysen, P.; Peirs, J.; Van Slycken, J.; Faes, K.; Duchene, L. Effect of strain rate on the forming behaviour of sheet metals. J. Mater. Process. Technol. 2011, 211, 1457–1464. [Google Scholar] [CrossRef]
  13. Wu, R.; Hu, Q.; Li, M.; Cai, S.; Chen, J. Evaluation of the forming limit of incremental sheet forming based on ductile damage. J. Mater. Process. Technol. 2021, 287, 116497. [Google Scholar] [CrossRef]
  14. Guzmán, C.F.; Yuan, S.; Duchêne, L.; Flores, E.I.S.; Habraken, A.M. Damage prediction in single point incremental forming using an extended Gurson model. Int. J. Solids Struct. 2018, 151, 45–56. [Google Scholar] [CrossRef] [Green Version]
  15. Bouffioux, C.; Lequesne, C.; Vanhove, H.; Duflou, J.R.; Pouteau, P.; Duchêne, L.; Habraken, A.M. Experimental and numerical study of an AlMgSc sheet formed by an incremental process. J. Mater. Process. Technol. 2011, 211, 1684–1693. [Google Scholar] [CrossRef]
  16. Ambrogio, G.; Filice, L.; Guerriero, F.; Guido, R.; Umbrello, D. Prediction of incremental sheet forming process performance by using a neural network approach. Int. J. Adv. Manuf. Technol. 2011, 54, 921–930. [Google Scholar] [CrossRef]
  17. Centeno, G.; Bagudanch, I.; Martínez-Donaire, A.J.; Garcia-Romeu, M.L.; Vallellano, C. Critical analysis of necking and fracture limit strains and forming forces in single-point incremental forming. Mater. Des. 2014, 63, 20–29. [Google Scholar] [CrossRef]
  18. Gorji, M.; Berisha, B.; Manopulo, N.; Hora, P. Effect of through thickness strain distribution on shear fracture hazard and its mitigation by using multilayer aluminum sheets. J. Mater. Process. Technol. 2016, 232, 19–33. [Google Scholar] [CrossRef]
  19. Mirnia, M.J.; Shamsari, M. Numerical prediction of failure in single point incremental forming using a phenomenological ductile fracture criterion. J. Mater. Process. Technol. 2017, 244, 17–43. [Google Scholar] [CrossRef]
  20. Ostasevicius, V.; Jurenas, V.; Grigaliunas, V.; Eidukynas, D.; Bubulis, A.; Paleviciute, I. Incremental Forming Machine for Sheet Metal Parts. LT Patent LT2020 516, 31 March 2020. [Google Scholar]
  21. Silva, M.B.; Skjødt, M.; Atkins, A.G.; Bay, N.; Martins, P.A.F. Single-point incremental forming and formability-failure diagrams. J. Strain Anal. Eng. Des. 2008, 43, 15–35. [Google Scholar] [CrossRef]
  22. Allwood, J.M.; Shouler, D.R.; Tekkaya, A.E. The increased forming limits of incremental sheet forming processes. In Key Engineering Materials; Trans Tech Publications Ltd.: Bäch, Switzerland, 2007; Volume 344, pp. 621–628. [Google Scholar]
  23. Emmens, W.C.; Van den Boogaard, A.H. Strain in shear, and material behaviour in incremental forming. In Key Engineering Materials; Trans Tech Publications Ltd.: Bäch, Switzerland, 2007; Volume 344, pp. 519–526. [Google Scholar]
Figure 1. 3D scanning experimental setup for the investigation. 1: an experimental body—aluminium alloy sheet; 2: steel base frame; 3: piezoelectric actuator; 4: liner amplifier P200 (FLC Electronics AB, Sweden); 5: 3D scanning vibrometer PSV-500-3D-HV (Polytec GmbH, Germany).
Figure 1. 3D scanning experimental setup for the investigation. 1: an experimental body—aluminium alloy sheet; 2: steel base frame; 3: piezoelectric actuator; 4: liner amplifier P200 (FLC Electronics AB, Sweden); 5: 3D scanning vibrometer PSV-500-3D-HV (Polytec GmbH, Germany).
Materials 15 01024 g001
Figure 2. Scheme of excitation of ultrasonic 3D vibrations in the sheet. 1: metal sheet; 2: frame; 3,4: bimorph-type piezoelectric actuators.
Figure 2. Scheme of excitation of ultrasonic 3D vibrations in the sheet. 1: metal sheet; 2: frame; 3,4: bimorph-type piezoelectric actuators.
Materials 15 01024 g002
Figure 3. Stress-strain curves of the power-law model (dashed curves represent min and max range for calibration of the material model, and the solid line represents the calibrated power-law function).
Figure 3. Stress-strain curves of the power-law model (dashed curves represent min and max range for calibration of the material model, and the solid line represents the calibrated power-law function).
Materials 15 01024 g003
Figure 4. Erichsen cupping test: (a) scheme, (b) test die.
Figure 4. Erichsen cupping test: (a) scheme, (b) test die.
Materials 15 01024 g004
Figure 5. Experimental setup of incremental aluminum-alloy sheet forming. 1: frame for aluminum-sheet fixation; 2: aluminum-alloy sheet; 3: forming tool; 4: robot ABB IRB1200; 5: pressing-force sensor, FCS-4035-150; 6: controller of pressing-force sensor; 7: power amplifier; 8: piezoelectric actuators for ultrasonic excitation.
Figure 5. Experimental setup of incremental aluminum-alloy sheet forming. 1: frame for aluminum-sheet fixation; 2: aluminum-alloy sheet; 3: forming tool; 4: robot ABB IRB1200; 5: pressing-force sensor, FCS-4035-150; 6: controller of pressing-force sensor; 7: power amplifier; 8: piezoelectric actuators for ultrasonic excitation.
Materials 15 01024 g005
Figure 6. Amplitude-frequency characteristics of the metal sheet measured with a Polytec 3D scanning vibrometer in the frequency range from 2 kHz to 60 kHz after excitation of two differently arranged piezoelectric actuators (according to Figure 2): (a) actuator 4 is excited; (b) actuator 3 is excited; (c) both actuators are excited, where X (red) and Y (green) represent lateral vibrations and Z (blue) represents vibrations in normal direction.
Figure 6. Amplitude-frequency characteristics of the metal sheet measured with a Polytec 3D scanning vibrometer in the frequency range from 2 kHz to 60 kHz after excitation of two differently arranged piezoelectric actuators (according to Figure 2): (a) actuator 4 is excited; (b) actuator 3 is excited; (c) both actuators are excited, where X (red) and Y (green) represent lateral vibrations and Z (blue) represents vibrations in normal direction.
Materials 15 01024 g006
Figure 7. Experimental setup for investigation of friction force: (a) interaction and (b) calculation scheme of single-point tool in contact with aluminum alloy sheet. 1: frame for sheet; 2: aluminum sheet; 3: forming tool; 4; ABB IRB1200 robot; 5: mechanical torque sensor, STJ100; 6: BGI series digital force gauge.
Figure 7. Experimental setup for investigation of friction force: (a) interaction and (b) calculation scheme of single-point tool in contact with aluminum alloy sheet. 1: frame for sheet; 2: aluminum sheet; 3: forming tool; 4; ABB IRB1200 robot; 5: mechanical torque sensor, STJ100; 6: BGI series digital force gauge.
Materials 15 01024 g007
Figure 8. Modelling of the cupping test: (a) numerical model of the cupping test, (b) force-displacement plot of the cupping experiment and sensitivity of the simulation results to the strength coefficient (K) of the calibrated power-law material model (strain-hardening coefficient, n = 0.097; failure strain, εu = 0.38).
Figure 8. Modelling of the cupping test: (a) numerical model of the cupping test, (b) force-displacement plot of the cupping experiment and sensitivity of the simulation results to the strength coefficient (K) of the calibrated power-law material model (strain-hardening coefficient, n = 0.097; failure strain, εu = 0.38).
Materials 15 01024 g008
Figure 9. Simulated images of a spherical punch piercing an aluminum-alloy sheet during the Erichsen test: (a) isometric view, (b) at the bottom of the sheet.
Figure 9. Simulated images of a spherical punch piercing an aluminum-alloy sheet during the Erichsen test: (a) isometric view, (b) at the bottom of the sheet.
Materials 15 01024 g009
Figure 10. SPIF simulation model.
Figure 10. SPIF simulation model.
Materials 15 01024 g010
Figure 11. Loading-curve displacements in XYZ directions.
Figure 11. Loading-curve displacements in XYZ directions.
Materials 15 01024 g011
Figure 12. X-contact force vs. rotations at the different contact friction coefficients.
Figure 12. X-contact force vs. rotations at the different contact friction coefficients.
Materials 15 01024 g012
Figure 13. Contact force vs. rotations for different coefficients of contact friction: (a) z-force (vertical force), (b) resultant force.
Figure 13. Contact force vs. rotations for different coefficients of contact friction: (a) z-force (vertical force), (b) resultant force.
Materials 15 01024 g013
Figure 14. Reduction in shell thickness: (a) variation at elements, which, after a forming depth of 25 mm, has drawing angles of 47°, 30° and 17°; (b) distribution of shell thickness after a forming depth of 25 mm.
Figure 14. Reduction in shell thickness: (a) variation at elements, which, after a forming depth of 25 mm, has drawing angles of 47°, 30° and 17°; (b) distribution of shell thickness after a forming depth of 25 mm.
Materials 15 01024 g014
Figure 15. Effective strain: (a) distribution after a forming depth of 25 mm; (b) evolution in shell element 6285 during the SPIF simulation process.
Figure 15. Effective strain: (a) distribution after a forming depth of 25 mm; (b) evolution in shell element 6285 during the SPIF simulation process.
Materials 15 01024 g015
Figure 16. Stress-triaxiality evolutions in particular elements that were in contact with tools during the forming process: (a) cupping test; (b) SPIF process.
Figure 16. Stress-triaxiality evolutions in particular elements that were in contact with tools during the forming process: (a) cupping test; (b) SPIF process.
Materials 15 01024 g016
Figure 17. Variation of in-plane and out-of-plane shear strains during the SPIF process.
Figure 17. Variation of in-plane and out-of-plane shear strains during the SPIF process.
Materials 15 01024 g017
Figure 18. The first modes of transverse oscillations of an aluminum-alloy sheet fixed in a steel angle: (a) measured at a frequency of 78.8 Hz and (b) calculated at a frequency of 81.707 Hz.
Figure 18. The first modes of transverse oscillations of an aluminum-alloy sheet fixed in a steel angle: (a) measured at a frequency of 78.8 Hz and (b) calculated at a frequency of 81.707 Hz.
Materials 15 01024 g018
Figure 19. Erichsen cupping test results: (a) cupping curves for specimens 1–5; (b) nature of failure of specimens 1–5 after the cupping test.
Figure 19. Erichsen cupping test results: (a) cupping curves for specimens 1–5; (b) nature of failure of specimens 1–5 after the cupping test.
Materials 15 01024 g019
Figure 20. Graphical representation of vertical force dependency from different frictions when sheet thickness is 0.5 mm.
Figure 20. Graphical representation of vertical force dependency from different frictions when sheet thickness is 0.5 mm.
Materials 15 01024 g020
Figure 21. Graphical representation of the puncture force of an aluminum alloy sheet.
Figure 21. Graphical representation of the puncture force of an aluminum alloy sheet.
Materials 15 01024 g021
Figure 22. Section of incrementally formed aluminum specimen with schematic representation of measurement points.
Figure 22. Section of incrementally formed aluminum specimen with schematic representation of measurement points.
Materials 15 01024 g022
Table 1. Mechanical properties of aluminum alloy EN AW1050A H24 from datasheet.
Table 1. Mechanical properties of aluminum alloy EN AW1050A H24 from datasheet.
Modulus of Elasticity, GPaProof Stress, MPaTensile Strength, MPaElongation A, %Strength Coefficient K, MPaStrain-Hardening Coefficient n
Datasheet Values71>85105–1453–8
Range for curve fitting 2 ÷ 40140 ÷ 2000.05 ÷ 0.2
Calibrated values LS-Opt & LS-Dyna 37.91430.097
Table 2. Friction coefficient and friction-force measurement results.
Table 2. Friction coefficient and friction-force measurement results.
ParameterWithout Lubrication and VibrationWith LubricationWith Vibration
Friction Force, N3.151.631.87
Friction coefficient0.50.10.12
Table 3. Thickness of incrementally formed aluminum sheet at measurement points.
Table 3. Thickness of incrementally formed aluminum sheet at measurement points.
Measurement PointWall Thickness, mmΔ from Initial Thickness, mm
10.510.01
20.360.14
30.490.01
40.370.13
50.510.01
where Δ is the deviation of the deformed sheet cross section. The results of measurement and modelling correlate with one another.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ostasevicius, V.; Paulauskaite-Taraseviciene, A.; Paleviciute, I.; Jurenas, V.; Griskevicius, P.; Eidukynas, D.; Kizauskiene, L. Investigation of the Robotized Incremental Metal-Sheet Forming Process with Ultrasonic Excitation. Materials 2022, 15, 1024. https://doi.org/10.3390/ma15031024

AMA Style

Ostasevicius V, Paulauskaite-Taraseviciene A, Paleviciute I, Jurenas V, Griskevicius P, Eidukynas D, Kizauskiene L. Investigation of the Robotized Incremental Metal-Sheet Forming Process with Ultrasonic Excitation. Materials. 2022; 15(3):1024. https://doi.org/10.3390/ma15031024

Chicago/Turabian Style

Ostasevicius, Vytautas, Agne Paulauskaite-Taraseviciene, Ieva Paleviciute, Vytautas Jurenas, Paulius Griskevicius, Darius Eidukynas, and Laura Kizauskiene. 2022. "Investigation of the Robotized Incremental Metal-Sheet Forming Process with Ultrasonic Excitation" Materials 15, no. 3: 1024. https://doi.org/10.3390/ma15031024

APA Style

Ostasevicius, V., Paulauskaite-Taraseviciene, A., Paleviciute, I., Jurenas, V., Griskevicius, P., Eidukynas, D., & Kizauskiene, L. (2022). Investigation of the Robotized Incremental Metal-Sheet Forming Process with Ultrasonic Excitation. Materials, 15(3), 1024. https://doi.org/10.3390/ma15031024

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