Interaction of Process Parameters, Forming Mechanisms, and Residual Stresses in Single Point Incremental Forming

: The residual stress state of a sheet metal component manufactured by metal forming has a signiﬁcant inﬂuence on the mechanical properties, and thus determines the time until the component fails, especially for dynamic loads. The origin of the resulting residual stress state of incrementally formed parts with regard to the forming mechanisms of shearing, bending, and the normal stress component is still under investigation. The relationship between the process parameters, the forming mechanisms, and the resulting residual stress state for a complex part geometry manufactured by single point incremental forming (SPIF) is presented in this publication. For this purpose, a validated numerical process model is used to analyze the inﬂuence of the step-down increment ∆ z for truncated cones on the characteristics of the forming mechanisms and the resulting residual stress state. For the ﬁrst time the forming mechanisms are evaluated numerically on both sides of the formed component. A relationship between the process parameters, forming mechanisms, residual stresses, and the mechanical properties of an incrementally formed component is shown. Shearing-induced hardening is identiﬁed as a relevant inﬂuence on the residual stress state of cones.


Introduction
The high formability in incremental sheet metal forming is a result of stress superposition during the process in a locally confined area that reduces the process forces and enables the processing of high-strength materials and complexly shaped components in a single setup compared to conventional forming processes. The complex relation of stress formation during the forming process that leads to the increased formability is due to the superposition of the forming mechanisms bending, shearing, and normal components [1]. The appearance of the forming mechanisms is described in various publications. The influence of each forming mechanism on the forming process has been discussed in different ways and is still not clear. Silva et al. [2] indicate the influence of membrane stretching as an important mechanism in single point incremental forming (SPIF) related to the sheet thickness reduction. Sebastiani et al. [3] experimentally demonstrated membrane stretching perpendicular to the tool movement direction. This effect increases by decreased vertical tool step-down increments and forming tool diameters [4]. The influence of the bending mechanism as an essential forming mechanism is proposed by Emmens et al. [1]. Larger strains on the tool side than on the tool-averted side and non-linear stresses on the tool side are shown as results of the bending mechanism during the forming operation [5]. In SPIF, bending always occurs in combination with a stretching component perpendicular to the tool feed [6], therefore additional tensile residual stresses on the tool side can be measured. Sebastiani et al. [3] detected larger plastic strains perpendicular to the tool feed on the tool side than on the tool-averted side. Recently, Maqbool et al. [4] explained the poor geometry accuracy as a result of bending-related springback effects. The importance of the shearing mechanism for this process is related to the comparison of single point incremental forming and the process of metal spinning [7,8]. Three different types of shearing that are relevant for SPIF can be distinguished. In-plane shear [9,10], shearing through the sheet thickness, and out-of-plane shear [3]. Maqbool et al. [4] demonstrated increasing through-thickness shear in the tool direction with decreasing step-down increment, increasing sheet thickness, or increased friction. Additionally, through-thickness shear increases with higher wall angles [11]. The shear due to the tool feed leads to a shear gradient in through-thickness direction [12], but is not an essential forming mechanism for incremental sheet metal forming [3]. Silva et al. [2] explained by membrane theory that stretching, not shear, is the dominant forming mechanism in the SPIF process. Out of plane shear was only proposed by numerical investigations and not yet validated [8]. The contact pressure of the forming tool during the forming operation leads to additional compressive stresses in the forming zone. These compressive stresses stabilize the forming process and lead to an increased formability by increasing the flow stress of the material [13]. In a recent publication [14], the influence of the forming mechanism on the residual stress development for simple linear groove geometries was analyzed. The influence of the step-down increment was shown by texture analysis.
In this work, the dependence of the forming mechanisms shearing, bending, and normal components on the process parameter ∆z/R Tool is shown by forming a complex component geometry. The step-down increment ∆z is varied in a numerical process model while the tool radius R Tool and the other process parameters as well as the final geometry are kept constant. The numerical process model is validated by comparing the vertical tool force, the component geometry, and the sheet thickness with experiments. The determination of the resulting residual stresses and their development with regard to the forming mechanisms, as a result of the process parameter adjustment, is then carried out numerically. It is shown how the residual stress state is influenced by the forming mechanisms. Based on published results linking the relationship between the residual stress state and the fatigue strength of a component, it is shown how the forming mechanisms affect the mechanical properties of the manufactured components.

Materials and Methods
For the validation of the numerical process model truncated cone geometries are manufactured by single point incremental forming. The dimensions of the truncated cone geometry are shown in Figure 1a. A conventional 5-axis DMU 50 milling machine from DMG Mori (Bielefeld, Germany) was used for this step (Figure 2a). The truncated cones were made of aluminum alloy AA5083 with an initial sheet thickness of t 0 = 1 mm. The production is carried out with a bidirectional tool path with a driven forming tool with a spherical end and a radius R Tool = 7.5 mm (Figure 1b). The step-down increment is varied in three steps: ∆z 24 = 1.875 mm, ∆z 12 = 3.75 mm, and ∆z 8 = 5.625 mm. All other process parameters are kept constant according to Table 1. The sheet was fixed with a clamping frame (Figure 2b). A deep-drawing oil (Castrol Iloform PN 226) was applied to the sheet metal surface as lubrication.
Metals 2020, 10, x FOR PEER REVIEW 2 of 12 stresses on the tool side can be measured. Sebastiani et al. [3] detected larger plastic strains perpendicular to the tool feed on the tool side than on the tool-averted side. Recently, Maqbool et al. [4] explained the poor geometry accuracy as a result of bending-related springback effects. The importance of the shearing mechanism for this process is related to the comparison of single point incremental forming and the process of metal spinning [7,8]. Three different types of shearing that are relevant for SPIF can be distinguished. In-plane shear [9,10], shearing through the sheet thickness, and out-of-plane shear [3]. Maqbool et al. [4] demonstrated increasing through-thickness shear in the tool direction with decreasing step-down increment, increasing sheet thickness, or increased friction. Additionally, through-thickness shear increases with higher wall angles [11]. The shear due to the tool feed leads to a shear gradient in through-thickness direction [12], but is not an essential forming mechanism for incremental sheet metal forming [3]. Silva et al. [2] explained by membrane theory that stretching, not shear, is the dominant forming mechanism in the SPIF process. Out of plane shear was only proposed by numerical investigations and not yet validated [8]. The contact pressure of the forming tool during the forming operation leads to additional compressive stresses in the forming zone. These compressive stresses stabilize the forming process and lead to an increased formability by increasing the flow stress of the material [13]. In a recent publication [14], the influence of the forming mechanism on the residual stress development for simple linear groove geometries was analyzed. The influence of the step-down increment was shown by texture analysis. In this work, the dependence of the forming mechanisms shearing, bending, and normal components on the process parameter Δz/RTool is shown by forming a complex component geometry. The step-down increment Δz is varied in a numerical process model while the tool radius RTool and the other process parameters as well as the final geometry are kept constant. The numerical process model is validated by comparing the vertical tool force, the component geometry, and the sheet thickness with experiments. The determination of the resulting residual stresses and their development with regard to the forming mechanisms, as a result of the process parameter adjustment, is then carried out numerically. It is shown how the residual stress state is influenced by the forming mechanisms. Based on published results linking the relationship between the residual stress state and the fatigue strength of a component, it is shown how the forming mechanisms affect the mechanical properties of the manufactured components.

Materials and Methods
For the validation of the numerical process model truncated cone geometries are manufactured by single point incremental forming. The dimensions of the truncated cone geometry are shown in Figure 1a. A conventional 5-axis DMU 50 milling machine from DMG Mori (Bielefeld, Germany) was used for this step (Figure 2a). The truncated cones were made of aluminum alloy AA5083 with an initial sheet thickness of t0 = 1 mm. The production is carried out with a bidirectional tool path with a driven forming tool with a spherical end and a radius RTool = 7.5 mm (Figure 1b). The step-down increment is varied in three steps: Δz24 = 1.875 mm, Δz12 = 3.75 mm, and Δz8 = 5.625 mm. All other process parameters are kept constant according to Table 1. The sheet was fixed with a clamping frame (Figure 2b). A deep-drawing oil (Castrol Iloform PN 226) was applied to the sheet metal surface as lubrication.  The process force of the force component Fz was determined using a 3-component piezoelectric force sensor (Figure 2b). Only the force component in the z-direction is considered here. To determine the geometrical deviation and the distribution of sheet thickness, a three-dimensional measurement is performed using the optical digitizer ATOS Triple Scan form GOM (Braunschweig, Germany). A precise fringe pattern is projected onto the surface of the object and recorded by two cameras. Based on the stereo camera principle, a 3D model of the part is generated. The sheet thickness distribution of the formed component can be determined locally using a scan of the front and back side of the part. Since the residual stress state of a component is highly dependent on the geometry, it must be ensured that a variation of the process parameters and a resulting change of the residual stress is not due to a deviation of the final geometry. The numerical model of the process was built in ABAQUS/Explicit and consists of the forming tool as a rigid body and the elastic-plastic sheet metal built of solid elements (C3D8) with an edge length of 1 mm and five elements over the sheet thickness ( Figure 3). An explicit time integration with a time step of 1 × 10 −06 was used. The computations are done using the high performance cluster (HPC) LiDO3. The total CPU time spent was approximately 32 h. The material model for the sheet includes combined isotropic-kinematic hardening based on Chaboche. The AA5083 base material was characterized using the uniaxial tensile test and in-plane The process force of the force component F z was determined using a 3-component piezoelectric force sensor (Figure 2b). Only the force component in the z-direction is considered here. To determine the geometrical deviation and the distribution of sheet thickness, a three-dimensional measurement is performed using the optical digitizer ATOS Triple Scan form GOM (Braunschweig, Germany). A precise fringe pattern is projected onto the surface of the object and recorded by two cameras. Based on the stereo camera principle, a 3D model of the part is generated. The sheet thickness distribution of the formed component can be determined locally using a scan of the front and back side of the part. Since the residual stress state of a component is highly dependent on the geometry, it must be ensured that a variation of the process parameters and a resulting change of the residual stress is not due to a deviation of the final geometry. The numerical model of the process was built in ABAQUS/Explicit and consists of the forming tool as a rigid body and the elastic-plastic sheet metal built of solid elements (C3D8) with an edge length of 1 mm and five elements over the sheet thickness ( Figure 3). An explicit time integration with a time step of 1 × 10 −6 was used. The computations are done using the high performance cluster (HPC) LiDO3. The total CPU time spent was approximately 32 h. The process force of the force component Fz was determined using a 3-component piezoelectric force sensor (Figure 2b). Only the force component in the z-direction is considered here. To determine the geometrical deviation and the distribution of sheet thickness, a three-dimensional measurement is performed using the optical digitizer ATOS Triple Scan form GOM (Braunschweig, Germany). A precise fringe pattern is projected onto the surface of the object and recorded by two cameras. Based on the stereo camera principle, a 3D model of the part is generated. The sheet thickness distribution of the formed component can be determined locally using a scan of the front and back side of the part. Since the residual stress state of a component is highly dependent on the geometry, it must be ensured that a variation of the process parameters and a resulting change of the residual stress is not due to a deviation of the final geometry. The numerical model of the process was built in ABAQUS/Explicit and consists of the forming tool as a rigid body and the elastic-plastic sheet metal built of solid elements (C3D8) with an edge length of 1 mm and five elements over the sheet thickness ( Figure 3). An explicit time integration with a time step of 1 × 10 −06 was used. The computations are done using the high performance cluster (HPC) LiDO3. The total CPU time spent was approximately 32 h. The material model for the sheet includes combined isotropic-kinematic hardening based on Chaboche. The AA5083 base material was characterized using the uniaxial tensile test and in-plane The material model for the sheet includes combined isotropic-kinematic hardening based on Chaboche. The AA5083 base material was characterized using the uniaxial tensile test and in-plane torsion test ( Figure 4) according to Maaß et al. [14]. The material coefficients in Table 2 are determined inversely using the cyclic loading test results according to Chaboche's procedure [15]. The initial flow stress of the material was determined to be σ f,0 = 165 MPa. A penalty contact is used for forming tool and sheet. The friction coefficient between the well lubricated forming tool and the sheet was set to µ = 0.03 as revealed in strip-tension tests. The numerical simulation starts with a simulation of the incremental single point forming process with ABAQUS/Explicit including mass scaling (mass scaling factor 10 3 ). The explicit simulation to evaluate the forming mechanisms is followed by a subsequent implicit simulation of the unclamping to determine the residual stresses.
Metals 2020, 10, x FOR PEER REVIEW 4 of 12 torsion test ( Figure 4) according to Maaß et al. [14]. The material coefficients in Table 2 are determined inversely using the cyclic loading test results according to Chaboche's procedure [15]. The initial flow stress of the material was determined to be σf,0 = 165 MPa. A penalty contact is used for forming tool and sheet. The friction coefficient between the well lubricated forming tool and the sheet was set to μ = 0.03 as revealed in strip-tension tests. The numerical simulation starts with a simulation of the incremental single point forming process with ABAQUS/Explicit including mass scaling (mass scaling factor 10 3 ). The explicit simulation to evaluate the forming mechanisms is followed by a subsequent implicit simulation of the unclamping to determine the residual stresses.

Results
The process force, the geometry, and sheet thickness distribution from the experiments are used for the validation of the numerical model.

Process Force
A comparison of the experimentally measured in-situ process forces Fz for the three different step-down increments Δz is shown in Figure 5.
An increasing step-down increment Δz results in a reduction in the process time. The maximum process force Fz of the considered process window corresponds to a reduction of 46% comparing 24 increments (Fz,24 = 1.89 ± 0.04 kN) and 8 increments (Fz,8 = 3.49 ± 0.07 kN). There is good agreement between the experimental and the numerical results with a maximum deviation of 8%. The peak forces are a result of the tool step-down. The time span between two peaks decreases due to the decreasing diameter of the tool path circle with increasing depth of the cone.

Results
The process force, the geometry, and sheet thickness distribution from the experiments are used for the validation of the numerical model.

Process Force
A comparison of the experimentally measured in-situ process forces F z for the three different step-down increments ∆z is shown in Figure 5.
An increasing step-down increment ∆z results in a reduction in the process time. The maximum process force F z of the considered process window corresponds to a reduction of 46% comparing 24 increments (F z,24 = 1.89 ± 0.04 kN) and 8 increments (F z,8 = 3.49 ± 0.07 kN). There is good agreement between the experimental and the numerical results with a maximum deviation of 8%. The peak forces are a result of the tool step-down. The time span between two peaks decreases due to the decreasing diameter of the tool path circle with increasing depth of the cone.

Geometry
A comparison of the geometry (cut A-A) of experimental results is shown in Figure 6a. With a decreasing number of step-down increments, there is an increasing waviness of the inner surface due to the increase of the step-down increment amount (Figure 6b). The maximum geometric deviation from the CAD target geometry is 11% due to process-related springback effects. The deviation is defined as the maximum distance between two points on the profiles with the same x-coordinate. The highest deviation is measured in the flange area after unclamping. However, the accuracy of the final geometry increases with a decreasing number of step-down increments due to an increased springback as shown in Ambrogio et al. [16], the maximum deviation between the three final geometries is 4.8%. Therefore, it can be assumed that a variation of the step-down increment, within the examined limits, has no significant effect on the resulting geometry of the component. Regarding the numerical results, a comparison of the final experimental geometry and the numerically determined geometry of the truncated cone is shown in Figure 7.

Geometry
A comparison of the geometry (cut A-A) of experimental results is shown in Figure 6a. With a decreasing number of step-down increments, there is an increasing waviness of the inner surface due to the increase of the step-down increment amount (Figure 6b). The maximum geometric deviation from the CAD target geometry is 11% due to process-related springback effects. The deviation is defined as the maximum distance between two points on the profiles with the same x-coordinate. The highest deviation is measured in the flange area after unclamping. However, the accuracy of the final geometry increases with a decreasing number of step-down increments due to an increased springback as shown in Ambrogio et al. [16], the maximum deviation between the three final geometries is 4.8%. Therefore, it can be assumed that a variation of the step-down increment, within the examined limits, has no significant effect on the resulting geometry of the component.

Geometry
A comparison of the geometry (cut A-A) of experimental results is shown in Figure 6a. With a decreasing number of step-down increments, there is an increasing waviness of the inner surface due to the increase of the step-down increment amount (Figure 6b). The maximum geometric deviation from the CAD target geometry is 11% due to process-related springback effects. The deviation is defined as the maximum distance between two points on the profiles with the same x-coordinate. The highest deviation is measured in the flange area after unclamping. However, the accuracy of the final geometry increases with a decreasing number of step-down increments due to an increased springback as shown in Ambrogio et al. [16], the maximum deviation between the three final geometries is 4.8%. Therefore, it can be assumed that a variation of the step-down increment, within the examined limits, has no significant effect on the resulting geometry of the component. Regarding the numerical results, a comparison of the final experimental geometry and the numerically determined geometry of the truncated cone is shown in Figure 7. Regarding the numerical results, a comparison of the final experimental geometry and the numerically determined geometry of the truncated cone is shown in Figure 7. An increasing number of increments leads to an improved geometrical accuracy. The maximum deviation improves Δzmax,8 = 0.67 ± 0.01 mm for eight increments to just Δzmax,24 = 0.49 ± 0.01 mm for 24 increments used.

Sheet Thickness
The sheet thickness distribution along the cut A-A is shown in Figure 8. The initial sheet thickness of t0 = 0.96 ± 0.01 mm decreases with an increasing number of step-down increments. The decreasing sheet thickness reduction varies between 26.8% for 8 increments (Figure 8a) and 44.7% for 24 increments (Figure 8c). The experimentally determined sheet thickness is smoothened to reduce the high level of noise due to the light reflections during the optical sheet thickness measurement. The maximum deviation between experiment and numerical results is 9% for the sheet thickness. With a decreasing number of increments, the thickness distribution along the cone wall gets more homogeneous. Depending on the incremental depth Δz some part of the material is deformed repeatedly due to an overlap when the forming tool is moving forward along the tool trajectory. If the number of increments is decreased there are less over-laps [17]. Due to the incremental forming of the cone wall, the initial sheet thickness is only reduced in the cone wall area. In the clamped flange area and the unformed bottom area of the part, the initial sheet thickness remains constant.

Hardness
Vickers hardness measurements enable the identification of various mechanical properties. The Vickers hardness measurement corresponding to HV 5.0 is used to determine the average hardness of the initial sheet (Figure 9a). An increasing number of increments leads to an improved geometrical accuracy. The maximum deviation improves ∆z max,8 = 0.67 ± 0.01 mm for eight increments to just ∆z max,24 = 0.49 ± 0.01 mm for 24 increments used.

Sheet Thickness
The sheet thickness distribution along the cut A-A is shown in Figure 8. The initial sheet thickness of t 0 = 0.96 ± 0.01 mm decreases with an increasing number of step-down increments. The decreasing sheet thickness reduction varies between 26.8% for 8 increments (Figure 8a) and 44.7% for 24 increments (Figure 8c). The experimentally determined sheet thickness is smoothened to reduce the high level of noise due to the light reflections during the optical sheet thickness measurement. The maximum deviation between experiment and numerical results is 9% for the sheet thickness. With a decreasing number of increments, the thickness distribution along the cone wall gets more homogeneous. Depending on the incremental depth ∆z some part of the material is deformed repeatedly due to an overlap when the forming tool is moving forward along the tool trajectory. If the number of increments is decreased there are less over-laps [17]. Due to the incremental forming of the cone wall, the initial sheet thickness is only reduced in the cone wall area. In the clamped flange area and the unformed bottom area of the part, the initial sheet thickness remains constant. An increasing number of increments leads to an improved geometrical accuracy. The maximum deviation improves Δzmax,8 = 0.67 ± 0.01 mm for eight increments to just Δzmax,24 = 0.49 ± 0.01 mm for 24 increments used.

Sheet Thickness
The sheet thickness distribution along the cut A-A is shown in Figure 8. The initial sheet thickness of t0 = 0.96 ± 0.01 mm decreases with an increasing number of step-down increments. The decreasing sheet thickness reduction varies between 26.8% for 8 increments (Figure 8a) and 44.7% for 24 increments (Figure 8c). The experimentally determined sheet thickness is smoothened to reduce the high level of noise due to the light reflections during the optical sheet thickness measurement. The maximum deviation between experiment and numerical results is 9% for the sheet thickness. With a decreasing number of increments, the thickness distribution along the cone wall gets more homogeneous. Depending on the incremental depth Δz some part of the material is deformed repeatedly due to an overlap when the forming tool is moving forward along the tool trajectory. If the number of increments is decreased there are less over-laps [17]. Due to the incremental forming of the cone wall, the initial sheet thickness is only reduced in the cone wall area. In the clamped flange area and the unformed bottom area of the part, the initial sheet thickness remains constant.

Hardness
Vickers hardness measurements enable the identification of various mechanical properties. The Vickers hardness measurement corresponding to HV 5.0 is used to determine the average hardness of the initial sheet (Figure 9a).

Hardness
Vickers hardness measurements enable the identification of various mechanical properties. The Vickers hardness measurement corresponding to HV 5.0 is used to determine the average hardness of the initial sheet (Figure 9a). The measurements are carried out using the cross-sections of the component opposite to the step-down mark in the middle of the formed cone wall with three repetitions considering the whole sheet thickness. The initial material hardness of the AA 5083 base material of 68.9 HV 5.0 is significantly increased by the process.
Regarding the number of increments, the strain-induced hardening is stationary with an increasing number of increments (Figure 9a). Figure 9b presents the normalized flow stress σf/σf,0 at the tool-side and the tool-averted side for a different number of increments. The normalized flow stress is increased with an increasing number of increments while the experimental hardness measurements are stationary (Figure 9a). The increase on the tool-averted side is higher than on the tool side for the amount of increments, but there is a turnover with an increasing number of increments.

Residual Stresses
The residual stress prediction of the numerical model was already experimentally validated by residual stress measurements using X-ray diffraction in Maaß et al. [14]. In this publication residual stresses were measured on both sides of linear grooves in tool feed direction and perpendicular to the tool path direction in the same measuring point. The overall variation of the numerical results compared to the experimental results was determined as 7%. For the present study, the residual stress amplitudes are output on both sides of the sheet in the center of the formed component wall at a height of 22.5 mm (middle of the cone wall, Figure 1a). The stress distribution of σ22 residual stresses in the cone wall is almost rotationally symmetric (except for the tool step-down area) and follows the pattern shown as an example for a sliced area around the evaluation point in Figure 10.  The measurements are carried out using the cross-sections of the component opposite to the step-down mark in the middle of the formed cone wall with three repetitions considering the whole sheet thickness. The initial material hardness of the AA 5083 base material of 68.9 HV 5.0 is significantly increased by the process.
Regarding the number of increments, the strain-induced hardening is stationary with an increasing number of increments (Figure 9a). Figure 9b presents the normalized flow stress σ f /σ f,0 at the tool-side and the tool-averted side for a different number of increments. The normalized flow stress is increased with an increasing number of increments while the experimental hardness measurements are stationary (Figure 9a). The increase on the tool-averted side is higher than on the tool side for the amount of increments, but there is a turnover with an increasing number of increments.

Residual Stresses
The residual stress prediction of the numerical model was already experimentally validated by residual stress measurements using X-ray diffraction in Maaß et al. [14]. In this publication residual stresses were measured on both sides of linear grooves in tool feed direction and perpendicular to the tool path direction in the same measuring point. The overall variation of the numerical results compared to the experimental results was determined as 7%. For the present study, the residual stress amplitudes are output on both sides of the sheet in the center of the formed component wall at a height of 22.5 mm (middle of the cone wall, Figure 1a). The stress distribution of σ 22 residual stresses in the cone wall is almost rotationally symmetric (except for the tool step-down area) and follows the pattern shown as an example for a sliced area around the evaluation point in Figure 10. The measurements are carried out using the cross-sections of the component opposite to the step-down mark in the middle of the formed cone wall with three repetitions considering the whole sheet thickness. The initial material hardness of the AA 5083 base material of 68.9 HV 5.0 is significantly increased by the process.
Regarding the number of increments, the strain-induced hardening is stationary with an increasing number of increments (Figure 9a). Figure 9b presents the normalized flow stress σf/σf,0 at the tool-side and the tool-averted side for a different number of increments. The normalized flow stress is increased with an increasing number of increments while the experimental hardness measurements are stationary (Figure 9a). The increase on the tool-averted side is higher than on the tool side for the amount of increments, but there is a turnover with an increasing number of increments.

Residual Stresses
The residual stress prediction of the numerical model was already experimentally validated by residual stress measurements using X-ray diffraction in Maaß et al. [14]. In this publication residual stresses were measured on both sides of linear grooves in tool feed direction and perpendicular to the tool path direction in the same measuring point. The overall variation of the numerical results compared to the experimental results was determined as 7%. For the present study, the residual stress amplitudes are output on both sides of the sheet in the center of the formed component wall at a height of 22.5 mm (middle of the cone wall, Figure 1a). The stress distribution of σ22 residual stresses in the cone wall is almost rotationally symmetric (except for the tool step-down area) and follows the pattern shown as an example for a sliced area around the evaluation point in Figure 10.  The residual stress in the evaluation point are output nearest to the surface along through thickness direction at the tool side and the tool averted. The local stress components in tangential direction σ 11 and meridional direction σ 22 to the cone wall surface for the tool side and tool-averted side are shown in Figure 11.
Metals 2020, 10, x FOR PEER REVIEW 8 of 12 The residual stress in the evaluation point are output nearest to the surface along through thickness direction at the tool side and the tool averted. The local stress components in tangential direction σ11 and meridional direction σ22 to the cone wall surface for the tool side and tool-averted side are shown in Figure 11. Considering the residual stresses σ22, that are perpendicular to the tool feed direction during one cycle transverse to the plane, the tensile residual stresses on the tool side increase monotonously from zero up to σ22 = 144 MPa when increasing the number of increments. The change of the residual stresses on the tool-averted side is less distinct with an increasing number of increments. The residual stresses rise to σ22 = −179 MPa for 24 increments from σ22 = −161 MPa for eight increments. In an opposing trend on the tool side, residual stresses in the tool feed direction σ11 decrease with an increasing number of increments. The tensile residual stresses in the tool feed direction σ11 on the tool-averted side follow those of the tool feed direction σ22 on the tool-averted side and are rather stationary. The residual shear stress component σ12 remains almost stationary with an increasing number of increments in the range of the σ11 component for 24 increments on the tool side. On the tool averted-side, σ12 also remain stationary in the range of σ12 = −20 MPa.

Forming Mechanisms
The forming mechanisms are quantified as described in Maaß et al. [18]. Derived from the output of the stress and strain components of the numerical model, the forming mechanisms are evaluated at the point of interest depicted in Figure 1a according to the locally defined 1-, 2-and 3-directions. The integration points of five linear brick elements in through thickness direction with four integration points in two planes per element are evaluated regarding the relevant stress and strain components σij (Cauchy stress) and dεij (total logarithmic strain increment). The normal stress components of a volume element, which represent the tension compression part, are named σn1, σn2, σn3, and are the average values of, e.g., σ11 over all volume elements through the sheet thickness. This gives four values of e.g., σn1 per element. The shearing components σ12, σ23, σ13 are the direct tensor output values of ABAQUS at a considered integration point (1). The incremental strain components are treated equally (2 As an assumption, pure bending components in terms of both stress and strain cancel out when integrating over the sheet thickness giving the normal components mentioned above. Only the bending components around the 1-and 2-directions are considered. Bending around the 3-direction Considering the residual stresses σ 22 , that are perpendicular to the tool feed direction during one cycle transverse to the plane, the tensile residual stresses on the tool side increase monotonously from zero up to σ 22 = 144 MPa when increasing the number of increments. The change of the residual stresses on the tool-averted side is less distinct with an increasing number of increments. The residual stresses rise to σ 22 = −179 MPa for 24 increments from σ 22 = −161 MPa for eight increments. In an opposing trend on the tool side, residual stresses in the tool feed direction σ 11 decrease with an increasing number of increments. The tensile residual stresses in the tool feed direction σ 11 on the tool-averted side follow those of the tool feed direction σ 22 on the tool-averted side and are rather stationary. The residual shear stress component σ 12 remains almost stationary with an increasing number of increments in the range of the σ 11 component for 24 increments on the tool side. On the tool averted-side, σ 12 also remain stationary in the range of σ 12 = −20 MPa.

Forming Mechanisms
The forming mechanisms are quantified as described in Maaß et al. [18]. Derived from the output of the stress and strain components of the numerical model, the forming mechanisms are evaluated at the point of interest depicted in Figure 1a according to the locally defined 1-, 2-and 3-directions. The integration points of five linear brick elements in through thickness direction with four integration points in two planes per element are evaluated regarding the relevant stress and strain components σ ij (Cauchy stress) and dε ij (total logarithmic strain increment). The normal stress components of a volume element, which represent the tension compression part, are named σ n1 , σ n2 , σ n3 , and are the average values of, e.g., σ 11 over all volume elements through the sheet thickness. This gives four values of e.g., σ n1 per element. The shearing components σ 12 , σ 23 , σ 13 are the direct tensor output values of ABAQUS at a considered integration point (1). The incremental strain components are treated equally (2).
As an assumption, pure bending components in terms of both stress and strain cancel out when integrating over the sheet thickness giving the normal components mentioned above. Only the bending components around the 1-and 2-directions are considered. Bending around the 3-direction is not present in the SPIF process and is therefore neglected. The components for bending (index b) are then defined by the superposition according to Equations (3) and (4).
The determination of the specific dissipated energy is then done by evaluating each forming mechanism for the considered volume element according to Equations (5)- (7) where the integration limits are from t = 0 until the process end. It is noted that these values are the average values of all integration points of an element.
An evaluation of the forming mechanisms on the tool side shows a significant and continuous increase of shear with an increasing number of step-down increments (Figure 12a). is not present in the SPIF process and is therefore neglected. The components for bending (index b) are then defined by the superposition according to Equations (3) and (4).
The determination of the specific dissipated energy is then done by evaluating each forming mechanism for the considered volume element according to Equations (5)- (7) where the integration limits are from t = 0 until the process end. It is noted that these values are the average values of all integration points of an element.
The resulting values are verified by comparing the sum of all energy values with the total plastic internal energy dissipation given by ABAQUS, wtot = wNormal + wShearing + wBending. The forming mechanisms are calculated at the near-surface elements at the tool side and at the opposite toolaverted side.
An evaluation of the forming mechanisms on the tool side shows a significant and continuous increase of shear with an increasing number of step-down increments (Figure 12a). Simultaneously, the bending mechanism decreases continuously with an increasing number of increments, in an opposing trend. The normal component on the tool side varies within a small range, remaining almost constant with an increasing number of step-down increments. At the tool-averted side, the qualitative development of the forming mechanisms is equal to the tool side ( Figure 12b). However, the forming mechanisms shearing and the normal components are at a lower level than on the tool side. The bending mechanism is slightly higher on the tool-averted side than on the tool side. Although the used element type may be prone to shear locking, this was not an obvious problem in the presented study as both the displacements or geometry and the forces were validated experimentally, meaning that a stiffness over prediction was not observed in the simulation results.

Discussion
The accordance of the geometry profiles within the examined number of step-down increments (Figure 6b) suggests that a change in the residual stress state with regard to the number of increments used from σ22.8 = 5 MPa to σ22.24 = 140 MPa on the tool side ( Figure 11) is not primarily caused by a change in geometry. It can be seen that the increase in residual stress amplitudes is due to the process- Simultaneously, the bending mechanism decreases continuously with an increasing number of increments, in an opposing trend. The normal component on the tool side varies within a small range, remaining almost constant with an increasing number of step-down increments. At the tool-averted side, the qualitative development of the forming mechanisms is equal to the tool side ( Figure 12b). However, the forming mechanisms shearing and the normal components are at a lower level than on the tool side. The bending mechanism is slightly higher on the tool-averted side than on the tool side. Although the used element type may be prone to shear locking, this was not an obvious problem in the presented study as both the displacements or geometry and the forces were validated experimentally, meaning that a stiffness over prediction was not observed in the simulation results.

Discussion
The accordance of the geometry profiles within the examined number of step-down increments (Figure 6b) suggests that a change in the residual stress state with regard to the number of increments used from σ 22.8 = 5 MPa to σ 22.24 = 140 MPa on the tool side ( Figure 11) is not primarily caused by a change in geometry. It can be seen that the increase in residual stress amplitudes is due to the process-related decrease in sheet thickness ( Figure 8) and the associated material hardening of the initial AA5083 base material (Figure 9a). The reasons for these effects are identified evaluating the acting forming mechanisms shearing, bending, and the normal components ( Figure 12). In particular, the shearing, which increases solely with an increasing number of increments, seems to have an influence on the decrease in sheet thickness and the hardening of the material during the process. An increasing number of step-down increments strengthens shearing due to the high portion of already plasticized and hardened areas in the forming zone. The normal component, which reaches a higher level on the tool side due to the forming tool pressure, shows no significant dependence on the number of increments. The fact that σ 22 is the only residual stress that increases in the tensile range on the tool side ( Figure 11a) could, thus, be explained by the increasing amount of shearing when the number of increments is increased (Figure 12a) since no other mechanism increases. This increase in shearing, which can be observed on both sides of the sheet, therefore causes the strain hardening ( Figure 9). The hardening-induced increasing normalized flow stress (Figure 9b) enables the material to take a higher absolute amount of residual stress amplitudes [19]. Due to the kinematic process of the SPIF process, residual stresses are more pronounced on the tool side than on the tool-averted side.
On the tool-averted side, the amplitudes of the compressive residual stresses are stationary with an increasing number of increments. Analyzing the forming mechanisms shows an increase in shearing and a decrease in bending as well as a stationary behavior of the normal component on the tool-averted side with an increasing number of step-down increments (Figure 12b). Due to the lower level of shearing compared to the tool side (Figure 12a), the bending gains a relative influence. Although the experimentally determined hardness measurements show a notable increase in material hardness compared to the initial state, the different effects on the tool and tool-averted sides are blurred due to the large test indention of HV 5.0. A clear distinction between tool and tool-averted side is not possible.
Regarding the effects of residual stresses on the mechanical properties, Maaß et al. [20] provide detailed information. Using a cyclic loading test setup of heat-treated and not heat-treated specimens, as explained in Maaß et al. [20], with an increasing number of increments, the relationship between residual stresses on the tool side and the failure of a component is shown ( Figure 13). The failure can be clearly assigned to the residual stresses by comparing stressed specimens and stress-relieved specimens. A relationship between the tensile residual stresses on the tool side and a 42% reduction of the lifetime in cyclic load tests until a component failure is experimentally shown.
Metals 2020, 10, x FOR PEER REVIEW 10 of 12 related decrease in sheet thickness ( Figure 8) and the associated material hardening of the initial AA5083 base material (Figure 9a). The reasons for these effects are identified evaluating the acting forming mechanisms shearing, bending, and the normal components ( Figure 12). In particular, the shearing, which increases solely with an increasing number of increments, seems to have an influence on the decrease in sheet thickness and the hardening of the material during the process. An increasing number of step-down increments strengthens shearing due to the high portion of already plasticized and hardened areas in the forming zone. The normal component, which reaches a higher level on the tool side due to the forming tool pressure, shows no significant dependence on the number of increments. The fact that σ22 is the only residual stress that increases in the tensile range on the tool side ( Figure 11a) could, thus, be explained by the increasing amount of shearing when the number of increments is increased (Figure 12a) since no other mechanism increases. This increase in shearing, which can be observed on both sides of the sheet, therefore causes the strain hardening ( Figure 9). The hardening-induced increasing normalized flow stress (Figure 9b) enables the material to take a higher absolute amount of residual stress amplitudes [19]. Due to the kinematic process of the SPIF process, residual stresses are more pronounced on the tool side than on the tool-averted side.
On the tool-averted side, the amplitudes of the compressive residual stresses are stationary with an increasing number of increments. Analyzing the forming mechanisms shows an increase in shearing and a decrease in bending as well as a stationary behavior of the normal component on the tool-averted side with an increasing number of step-down increments (Figure 12b). Due to the lower level of shearing compared to the tool side (Figure 12a), the bending gains a relative influence. Although the experimentally determined hardness measurements show a notable increase in material hardness compared to the initial state, the different effects on the tool and tool-averted sides are blurred due to the large test indention of HV 5.0. A clear distinction between tool and tool-averted side is not possible.
Regarding the effects of residual stresses on the mechanical properties, Maaß et al. [20] provide detailed information. Using a cyclic loading test setup of heat-treated and not heat-treated specimens, as explained in Maaß et al. [20], with an increasing number of increments, the relationship between residual stresses on the tool side and the failure of a component is shown ( Figure 13). The failure can be clearly assigned to the residual stresses by comparing stressed specimens and stress-relieved specimens. A relationship between the tensile residual stresses on the tool side and a 42% reduction of the lifetime in cyclic load tests until a component failure is experimentally shown.

Conclusion and Outlook
It can be concluded that there is a direct relation between the process parameter step-down increment, the forming mechanisms, and the resulting residual stress state, which significantly influences the component properties. The component properties can be optimized for a specific load

Conclusions and Outlook
It can be concluded that there is a direct relation between the process parameter step-down increment, the forming mechanisms, and the resulting residual stress state, which significantly influences the component properties. The component properties can be optimized for a specific load case by a targeted induction of the residual stresses in the component. According to the load case of the formed component, the product properties can be directly influenced by adjusting the residual stress state. A basic understanding of the relationship between the process parameters and the residual stress formation in the process is necessary for the targeted setting of the residual stress state according to the application. This publication explains the influence of the process parameter relative step-down increment ∆z/R Tool on the forming mechanisms of the incremental sheet metal forming process and the resulting residual stress state in the component by a validated numerical process model. Due to the increasingly dominant shearing mechanism with an increasing number of step-down increments, the material thickness is reduced and hardened, particularly on the tool side. Beside other effects, the material hardening enables a higher residual stress amplitude in the material, the leads to higher residual stress amplitudes. As a result, the increased initial tensile residual stresses on the tool side lead to an early failure of a truncated cone under cyclic loading.
Further investigations should focus on the examination of the transferability of results to materials with significantly different hardening behavior compared to the analyzed material AA5083. A suitable experimental validation of the forming mechanisms could be carried out to confirm the numerically determined portions of the forming mechanisms.