Springback Prediction of Aluminum Alloy Sheet under Changing Loading Paths with Consideration of the Influence of Kinematic Hardening and Ductile Damage

Springback prediction of sheet metal forming is always an important issue in the industry, because it greatly affects the final shape of the product. The accuracy of simulation prediction depends on not only the forming condition but also the chosen material model, which determines the stress and strain redistributions in the formed parts. In this paper, a newly proposed elastoplastic constitutive model is used, in which the initial and induced anisotropies, combined nonlinear isotropic and kinematic hardenings, as well as isotropic ductile damage, are taken into account. The aluminum alloy sheet metal AA7055 was chosen as the studied material. In order to investigate springback under non-proportional strain paths, three-point bending tests were conducted with pre-strained specimens, and five different pre-straining states were considered. The comparisons between numerical and experimental results highlighted the hard effect of both kinematic hardening and ductile damage on the springback prediction, especially for a changed loading path case.


Introduction
In recent years, with the rapid development of the need for lightweight materials, more high-strength aluminum alloy sheets have been widely used in the products of automotive, aerospace, and medical health fields. However, some undesirable problems during the forming of high strength sheet metal have been observed. Springback is one of these important issues, and it will affect the final appearance and accuracy of the workpiece. Sheet metal forming usually induces important and complicated plastic strain state, particularly under complex loading paths and large strain. Understanding these behaviors of sheet metal forming becomes more important, which can give great help during the tool adjustments in the metal forming process, especially for the new lightweight and high-strength materials.
where e e denotes the deviatoric part of the elastic strain, · + and · − denote positive and negative parts. µ e and K e are the classical shear and bulk elastic moduli. h is the microcracks closure parameter 0 ≤ h ≤ 1, C and Q are the hardening moduli for kinematic and isotropic hardening, respectively, γ is the parameter of damage effect on isotropic hardening. Y e , Y α , and Y r are the density energy release rates corresponding respectively to elastoplastic flow, kinematic hardening, and isotropic hardening. Details about the state variables and evolution equations are given in Appendix A.
The given model will be implemented in both ABAQUS/Standard and ABAQUS/Explicit finite element codes through the user subroutines UMAT and VUMAT. The implementation of the developed model is based on purely local implicit integration scheme used with the elastic predictor-plastic corrector approach.

Experimental Investigations
The AA7055 aluminum alloy sheet with 1.60 mm thickness was used for the experimental study. In order to investigate the influence of kinematic hardening and ductile damage on springback under the loading path changing, a series of tests were conducted. Firstly, the simple uniaxial tensile test, cyclic shear test, and cyclic loading-unloading tensile tests were conducted separately to determine the material parameters. In order to investigate the springback of AA7055 alloy sheet the tensile specimens were subjected to five different pre-strain states, and then bent with 25 mm depth. With the neutral layer as the boundary, the inner sheet material experiences the cyclic tension-compression loading process during tensile-bending tests. During the bending tests, the internal damage appearing in the sheet increases with plastic strains and causes the decrease of the rigidity of the material. The increase of micro-cracks and the micro-voids in the sheet reduces the stress levels in the sheet, which will affect the final formed shape.

Uniaxial Tensile Tests
The geometry of the uniaxial tensile test specimen is given in Figure 1. The specimens were cut along the orientations (0 • , 45 • , 90 • ) according to the rolling direction. The tensile velocity was fixed to 2.0 mm/min. To ensure the repeatability of the results, each test was conducted at least three times. The uniaxial tensile test is used to identify the material properties, including the yield function, anisotropy, and hardening laws. The obtained plastic flow curves 0 • , 45 • , and 90 • according to the rolling direction are given in Figure 2. kinematic hardening, isotropic hardening, and isotropic ductile damage, respectively. Detailed description of the fully coupled relationship can be found in the published research works [22,23].
where e e denotes the deviatoric part of the elastic strain, .  and .  denote positive and negative parts. e μ and e K are the classical shear and bulk elastic moduli. h is the microcracks closure parameter 0 ≤ h ≤ 1, C and Q are the hardening moduli for kinematic and isotropic hardening, respectively, γ is the parameter of damage effect on isotropic hardening. Y e , Y α , and Y r are the density energy release rates corresponding respectively to elastoplastic flow, kinematic hardening, and isotropic hardening. Details about the state variables and evolution equations are given in Appendix A. The given model will be implemented in both ABAQUS/Standard and ABAQUS/Explicit finite element codes through the user subroutines UMAT and VUMAT. The implementation of the developed model is based on purely local implicit integration scheme used with the elastic predictorplastic corrector approach.

Experimental Investigations
The AA7055 aluminum alloy sheet with 1.60 mm thickness was used for the experimental study. In order to investigate the influence of kinematic hardening and ductile damage on springback under the loading path changing, a series of tests were conducted. Firstly, the simple uniaxial tensile test, cyclic shear test, and cyclic loading-unloading tensile tests were conducted separately to determine the material parameters. In order to investigate the springback of AA7055 alloy sheet the tensile specimens were subjected to five different pre-strain states, and then bent with 25 mm depth. With the neutral layer as the boundary, the inner sheet material experiences the cyclic tension-compression loading process during tensile-bending tests. During the bending tests, the internal damage appearing in the sheet increases with plastic strains and causes the decrease of the rigidity of the material. The increase of micro-cracks and the micro-voids in the sheet reduces the stress levels in the sheet, which will affect the final formed shape.

Uniaxial Tensile Tests
The geometry of the uniaxial tensile test specimen is given in Figure 1. The specimens were cut along the orientations (0°, 45°, 90°) according to the rolling direction. The tensile velocity was fixed to 2.0 mm/min. To ensure the repeatability of the results, each test was conducted at least three times. anisotropy, and hardening laws. The obtained plastic flow curves 0°, 45°, and 90° according to the rolling direction are given in Figure 2.

Cyclic Shear Loading Tests
The geometry of the cyclic shear loading test specimen is given in Figure 1. In the proposed model, the kinematic hardening parameters should be determined through the cyclic loading tests. Through the literature study, the cyclic shear test can be one valid method to investigate the Bauschinger phenomenon of sheet metal [24]. In this study, the cyclic shear loading test was conducted with a butterfly clamp on the universal tensile machine. The Digital image correlation (DIC) system was used to capture the vertical displacement between point A and point B shown in Figure 3. The engineering strain and stress can be calculated by: where Δh denotes the vertical displacement; h0 denotes the initial height of the connect zone; t0 denotes the initial thickness of the sheet; F denotes the loading force during the tests.

Cyclic Shear Loading Tests
The geometry of the cyclic shear loading test specimen is given in Figure 1. In the proposed model, the kinematic hardening parameters should be determined through the cyclic loading tests. Through the literature study, the cyclic shear test can be one valid method to investigate the Bauschinger phenomenon of sheet metal [24]. In this study, the cyclic shear loading test was conducted with a butterfly clamp on the universal tensile machine. The Digital image correlation (DIC) system was used to capture the vertical displacement between point A and point B shown in Figure 3. The engineering strain and stress can be calculated by: where ∆h denotes the vertical displacement; h 0 denotes the initial height of the connect zone; t 0 denotes the initial thickness of the sheet; F denotes the loading force during the tests. anisotropy, and hardening laws. The obtained plastic flow curves 0°, 45°, and 90° according to the rolling direction are given in Figure 2.

Cyclic Shear Loading Tests
The geometry of the cyclic shear loading test specimen is given in Figure 1. In the proposed model, the kinematic hardening parameters should be determined through the cyclic loading tests. Through the literature study, the cyclic shear test can be one valid method to investigate the Bauschinger phenomenon of sheet metal [24]. In this study, the cyclic shear loading test was conducted with a butterfly clamp on the universal tensile machine. The Digital image correlation (DIC) system was used to capture the vertical displacement between point A and point B shown in Figure 3. The engineering strain and stress can be calculated by: where Δh denotes the vertical displacement; h0 denotes the initial height of the connect zone; t0 denotes the initial thickness of the sheet; F denotes the loading force during the tests.  With the obtained stress-strain curves in uniaxial and cyclic shear tests, the material parameters suited to aluminum alloy sheet AA7055 can be determined through the use of an appropriate inverse hybrid numerical-experimental methodology [25]. The relevant numerical and experimental result are the displacement-load responses. The inverse method here is to search the minimum error value between simulation and experiment responses. The Trust Region reflective method is used, which suits itself to the nonlinear least squares optimization problem. The optimization process involves the approximate solution of a large linear system based on the method of preconditioned conjugate gradients. The optimization algorithm is coded within MATLAB and Python script in conjunction with the FE commercial software ABAQUS.
Note that the use of only these two tests cannot allow the accurate determination of all model parameters such as the determination of micro-cracks closure parameter h and distortional parameters (X c l1 , X c l2 , X p l1 ). In this study, the micro-cracks closure parameter is assumed to be h = 0.2. The distortional effect in this study will be ignored by considering X c l1 = X c l2 = X p l1 = ∞. The set of material parameters for the AA7055 is summarized and given in Table 1.
In this study, the decomposition of kinematic hardening and isotropic hardening was determined through cyclic shear tests ( Figure 4). But they also can be determined with the help of the springback angle after the bending test with different pre-stretch strain [13]. The initial yield criteria (Hill 48) in the yield function f and potential equation F are assigned the same values (F = F*, G = G*, H = H*, L = L*, M = M*). In the coming studies, this decomposition ratio will be changed in order to study the influence of kinematic hardening on the prediction of the final springback values.  With the obtained stress-strain curves in uniaxial and cyclic shear tests, the material parameters suited to aluminum alloy sheet AA7055 can be determined through the use of an appropriate inverse hybrid numerical-experimental methodology [25]. The relevant numerical and experimental result are the displacement-load responses. The inverse method here is to search the minimum error value between simulation and experiment responses. The Trust Region reflective method is used, which suits itself to the nonlinear least squares optimization problem. The optimization process involves the approximate solution of a large linear system based on the method of preconditioned conjugate gradients. The optimization algorithm is coded within MATLAB and Python script in conjunction with the FE commercial software ABAQUS.
Note that the use of only these two tests cannot allow the accurate determination of all model parameters such as the determination of micro-cracks closure parameter h and distortional parameters ( l1 c X , l2 c X , l1 p X ). In this study, the micro-cracks closure parameter is assumed to be h = 0.2. The distortional effect in this study will be ignored by considering l1 l2 l1 of material parameters for the AA7055 is summarized and given in Table 1.

Cyclic Loading-Unloading Uniaxial Tension Tests
In order to better investigate the damage initiation and growth during the plastic flow, uni-axial tensile test with loading-unloading cycles was conducted, and the evolution of elastic modulus was recorded during the unloading phases ( Figure 5). The specimen is prepared following the rolling direction. The strain levels of unloading are chosen (1.6%, 2.7%, 3.7%, 4.7%, 5.7%, 6.7%). The average value of elastic moduli during the unloading process and the decrease ratios (E i /E 0 ) are given in Figure 6. In this figure, we can observe clearly the decreases of elastic modulus according to the tensile strain which can be related to the damage evolution during the tensile loading. The decrease ratio of E can be regarded as the damage value in a way. During the simulation on the prediction of the springback phenomenon, the E decrease can be taken into account through user subroutine, and this result can compared with ones obtained with ductile damage model (marked as Simu-E in Figure 13). springback phenomenon, the E decrease can be taken into account through user subroutine, and this result can compared with ones obtained with ductile damage model (marked as Simu-E in Figure 13).

Bending Test with Pre-Stretch Strain
Three-point bending tests were conducted to investigate the springback phase of AA7055 specimens. The test design and device are given in Figure 7, and also the specimen is prepared following the rolling direction. Considering the maximum strain of 0.12 (where the local necking happens after this point) under uniaxial tensile loading path, five pre-straining states are used of 4.6% (2.4 mm), 6.4% (3.3 mm), 8.0% (4.2 mm), 9.7% (5.1 mm), and 11.3% (6.0 mm) were selected to assign on the bending specimens. The bending depth is about 25 mm, and the fixed velocity of the upper die is about 5.0 mm/min. The geometry sizes of the bending dies are shown in Figure 7. The springback angle is measured for the different specimens including the pre-straining states. springback phenomenon, the E decrease can be taken into account through user subroutine, and this result can compared with ones obtained with ductile damage model (marked as Simu-E in Figure 13).

Bending Test with Pre-Stretch Strain
Three-point bending tests were conducted to investigate the springback phase of AA7055 specimens. The test design and device are given in Figure 7, and also the specimen is prepared following the rolling direction. Considering the maximum strain of 0.12 (where the local necking happens after this point) under uniaxial tensile loading path, five pre-straining states are used of 4.6% (2.4 mm), 6.4% (3.3 mm), 8.0% (4.2 mm), 9.7% (5.1 mm), and 11.3% (6.0 mm) were selected to assign on the bending specimens. The bending depth is about 25 mm, and the fixed velocity of the upper die is about 5.0 mm/min. The geometry sizes of the bending dies are shown in Figure 7. The springback angle is measured for the different specimens including the pre-straining states.

Bending Test with Pre-Stretch Strain
Three-point bending tests were conducted to investigate the springback phase of AA7055 specimens. The test design and device are given in Figure 7, and also the specimen is prepared following the rolling direction. Considering the maximum strain of 0.12 (where the local necking happens after this point) under uniaxial tensile loading path, five pre-straining states are used of 4.6% (2.4 mm), 6.4% (3.3 mm), 8.0% (4.2 mm), 9.7% (5.1 mm), and 11.3% (6.0 mm) were selected to assign on the bending specimens. The bending depth is about 25 mm, and the fixed velocity of the upper die is about 5.0 mm/min. The geometry sizes of the bending dies are shown in Figure 7. The springback angle is measured for the different specimens including the pre-straining states.  The test design and part of final obtained samples are showed in Figure 8. In Figure 9 the obtained average springback angles after bending tools unloading are shown. From the comparison of the results it can be seen that plastic pre-straining affects greatly the resulting springback angles. We observe that the increase of tensile pre-strain induces an increase of springback angle. On the other hand, when the plastic pre-strain increases, the increases of hardening stress and plastic strain will be obtained on one side accompanying the increase of ductile damage. This damage evolution allows a significant reduction of rigidity of the material so that when the bending tools are unloaded the critical parts are submitted to the high level of stress with the low rigidity which will contribute to maximizing the springback phase.
(a) Test design (b) Test samples   The test design and part of final obtained samples are showed in Figure 8. In Figure 9 the obtained average springback angles after bending tools unloading are shown. From the comparison of the results it can be seen that plastic pre-straining affects greatly the resulting springback angles. We observe that the increase of tensile pre-strain induces an increase of springback angle. On the other hand, when the plastic pre-strain increases, the increases of hardening stress and plastic strain will be obtained on one side accompanying the increase of ductile damage. This damage evolution allows a significant reduction of rigidity of the material so that when the bending tools are unloaded the critical parts are submitted to the high level of stress with the low rigidity which will contribute to maximizing the springback phase.  The test design and part of final obtained samples are showed in Figure 8. In Figure 9 the obtained average springback angles after bending tools unloading are shown. From the comparison of the results it can be seen that plastic pre-straining affects greatly the resulting springback angles. We observe that the increase of tensile pre-strain induces an increase of springback angle. On the other hand, when the plastic pre-strain increases, the increases of hardening stress and plastic strain will be obtained on one side accompanying the increase of ductile damage. This damage evolution allows a significant reduction of rigidity of the material so that when the bending tools are unloaded the critical parts are submitted to the high level of stress with the low rigidity which will contribute to maximizing the springback phase.
(a) Test design (b) Test samples    The test design and part of final obtained samples are showed in Figure 8. In Figure 9 the obtained average springback angles after bending tools unloading are shown. From the comparison of the results it can be seen that plastic pre-straining affects greatly the resulting springback angles. We observe that the increase of tensile pre-strain induces an increase of springback angle. On the other hand, when the plastic pre-strain increases, the increases of hardening stress and plastic strain will be obtained on one side accompanying the increase of ductile damage. This damage evolution allows a significant reduction of rigidity of the material so that when the bending tools are unloaded the critical parts are submitted to the high level of stress with the low rigidity which will contribute to maximizing the springback phase.
(a) Test design (b) Test samples  . Springback angles after bending tools unloading versus pre-straining states. Figure 9. Springback angles after bending tools unloading versus pre-straining states.
Note that when the bending depth is the same, under different pre-strain levels, the stress and strain distributions across the critical deformed zone situated under the punch radius shall be quite different. Without considering pre-strains, the normally neutral layer is situated in the central side of the sheet thickness. Tensile pre-strains shall be added to the bending strain state and causes the sliding of the neutral layer in the direction of the punch. The movement of the neutral layer will reduce the volume ratio of the compression stress state in the inner thickness section. This fact will contribute also to spring back increase.

Simulation Results
In the simulation process, the influences of the kinematic hardening and the ductile damage on the springback were studied. The loading processes include pre-straining and bending phases that were simulated using ABAQUS/Explicit ® with user subroutine (VUMAT), which takes into account the dynamic process and avoid the non-convergence problem. The unloading process was simulated using ABAQUS/Standard ® with user subroutine (UMAT) in order to reflect the final static shape after unloading. The decrease of Young's modulus during the deformation process can be defined using two approaches: the direct assignment according to the decrease ratios of E given in Figure 6. The second one is obtained automatically through the full coupling with ductile damage. In the coming sections, the comparison of these two approaches is analyzed. In order to better investigate the effect of the stress gradient along thickness, the 2D plane strain geometrical model was chosen in the simulation process. Constant mesh size 0.1 mm is used in the critical deformation zone, as shown in Figure 10. Note that when the bending depth is the same, under different pre-strain levels, the stress and strain distributions across the critical deformed zone situated under the punch radius shall be quite different. Without considering pre-strains, the normally neutral layer is situated in the central side of the sheet thickness. Tensile pre-strains shall be added to the bending strain state and causes the sliding of the neutral layer in the direction of the punch. The movement of the neutral layer will reduce the volume ratio of the compression stress state in the inner thickness section. This fact will contribute also to spring back increase.

Simulation Results
In the simulation process, the influences of the kinematic hardening and the ductile damage on the springback were studied. The loading processes include pre-straining and bending phases that were simulated using ABAQUS/Explicit ® with user subroutine (VUMAT), which takes into account the dynamic process and avoid the non-convergence problem. The unloading process was simulated using ABAQUS/Standard ® with user subroutine (UMAT) in order to reflect the final static shape after unloading. The decrease of Young's modulus during the deformation process can be defined using two approaches: the direct assignment according to the decrease ratios of E given in Figure 6. The second one is obtained automatically through the full coupling with ductile damage. In the coming sections, the comparison of these two approaches is analyzed. In order to better investigate the effect of the stress gradient along thickness, the 2D plane strain geometrical model was chosen in the simulation process. Constant mesh size 0.1 mm is used in the critical deformation zone, as shown in Figure 10.

Effect of Kinematic Hardening Fraction on Springback Prediction
The hardening stresses obtained through the uniaxial tensile test combines isotropic and kinematic hardenings. Without the cyclic loading tests, it is difficult to separate these two hardenings and determine the exact fraction of every part. In this study, the kinematic hardening fraction is considered through a parametric study. The fraction ratio of kinematic hardening with respect to total hardening is defined by: where Xsat is the kinematic hardening saturation value and Rsat is the isotropic hardening saturation value. In the parametric study, the kinematic hardening fraction values chosen are t = (0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0). For all these kinematic hardening fractions, the values of hardening modulus Q and C are re-identified in order to obtain the same fit with experimental tensile test curves. In Table 2 are given the obtained sets of hardening parameters corresponding to each considered kinematic hardening fraction.

Effect of Kinematic Hardening Fraction on Springback Prediction
The hardening stresses obtained through the uniaxial tensile test combines isotropic and kinematic hardenings. Without the cyclic loading tests, it is difficult to separate these two hardenings and determine the exact fraction of every part. In this study, the kinematic hardening fraction is considered through a parametric study. The fraction ratio of kinematic hardening with respect to total hardening is defined by: where X sat is the kinematic hardening saturation value and R sat is the isotropic hardening saturation value. In the parametric study, the kinematic hardening fraction values chosen are t = (0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0). For all these kinematic hardening fractions, the values of hardening modulus Q and C are re-identified in order to obtain the same fit with experimental tensile test curves. In Table 2 are given the obtained sets of hardening parameters corresponding to each considered kinematic hardening fraction. In order to enhance the stress gradient description through the thickness, element size of 0.1 mm was considered along with the thickness direction. The simulation was consistent with the experimental setup, and five different tensile pre-strains 4.6%, 6.4%, 8.0%, 9.7%, and 11.3% were applied during a step before the bending process.
The friction coefficient between the specimen and the dies was fixed to 0.05. The bending depth of the upper die was about 25.0 mm. The pre-tension and three-point bending processes were simulated through the dynamic explicit methodology, while the unloading process after bending was conducted by the static analysis methodology. The state variables after the bending analysis were imported as the initial state of the unloading process. After the unloading process, the resulting values obtained of springback angles are given in Figure 11 with respect to the pre-strain states.  In order to enhance the stress gradient description through the thickness, element size of 0.1 mm was considered along with the thickness direction. The simulation was consistent with the experimental setup, and five different tensile pre-strains 4.6%, 6.4%, 8.0%, 9.7%, and 11.3% were applied during a step before the bending process.
The friction coefficient between the specimen and the dies was fixed to 0.05. The bending depth of the upper die was about 25.0 mm. The pre-tension and three-point bending processes were simulated through the dynamic explicit methodology, while the unloading process after bending was conducted by the static analysis methodology. The state variables after the bending analysis were imported as the initial state of the unloading process. After the unloading process, the resulting values obtained of springback angles are given in Figure 11 with respect to the pre-strain states. In Figure 11 are given the springback angles obtained after the bending tests applied for different considered pre-strained specimens without considering ductile damage. When pre-tensile plastic strain is zero, there is still a gap of 3°, which can be caused by the friction and the thickness variation of the sheet, and also the anisotropy following the thickness direction is ignored during the simulation. In this figure, the results of the simulations using the various kinematic hardening fraction cases (t = 0.0~1.0) are compared with the experimental results. The case of t = 0.0 value fraction corresponds to a purely isotropic hardening model; however, the case of t = 1.0 value fraction corresponds to purely kinematic hardening model. As mentioned before, the springback angle In Figure 11 are given the springback angles obtained after the bending tests applied for different considered pre-strained specimens without considering ductile damage. When pre-tensile plastic strain is zero, there is still a gap of 3 • , which can be caused by the friction and the thickness variation of the sheet, and also the anisotropy following the thickness direction is ignored during the simulation. In this figure, the results of the simulations using the various kinematic hardening fraction cases (t = 0.0~1.0) are compared with the experimental results. The case of t = 0.0 value fraction corresponds to a purely isotropic hardening model; however, the case of t = 1.0 value fraction corresponds to purely kinematic hardening model. As mentioned before, the springback angle increases when the pre-strain increases. In this figure, we observe an important effect of kinematic hardening fraction on the springback predicted angle. It needs to be highlighted that for high kinematic hardening fraction (t > 0.8) the increase of pre-strain state induces a decrease of springback angle. For these cases, the difference among numerical results is important. The stress contours before springback are given in Figure 12, which clearly present the big effect of kinematic hardening fraction on the final obtained stress contour on the thickness section during the bending process. increases when the pre-strain increases. In this figure, we observe an important effect of kinematic hardening fraction on the springback predicted angle. It needs to be highlighted that for high kinematic hardening fraction (t > 0.8) the increase of pre-strain state induces a decrease of springback angle. For these cases, the difference among numerical results is important. The stress contours before springback are given in Figure 12, which clearly present the big effect of kinematic hardening fraction on the final obtained stress contour on the thickness section during the bending process. Through comparing the experimental and numerical obtained springback angles, it is recommended that the kinematic hardening fraction shall be defined between t = 0.1 and t = 0.5 for Al7050. In the coming section about the ductile damage, the ratio of kinematic hardening in total hardening was assigned to be 0.5.

Effect of Damage Coupling on Springback Prediction
The Young's modulus plays an important role in the springback prediction. In this section, the decrease of Young's modulus during the deformation process is considered through the two approaches of including Young's moduli variation according to loading strain as discusses in the last section. The direct evolution of Young's moduli with respect to accumulating plastic strain is defined with tabular form according to Figure 6. Linear interpolation is considered between given points. For the ductile damage model, the material parameters were defined according to Table 2. The corresponding obtained numerically predicted springback angles are shown in Figure 13. These two methods can both well predict the tendency of the influence of pre-stretch on the final springback angles. Together with the increase of the pre-stretch displacement, the springback angles almost increase linearly with the same slope. Also, it can be found that the obtained springback angle with the ductile damage model can be closer to the experimental results. Through comparing the experimental and numerical obtained springback angles, it is recommended that the kinematic hardening fraction shall be defined between t = 0.1 and t = 0.5 for Al7050. In the coming section about the ductile damage, the ratio of kinematic hardening in total hardening was assigned to be 0.5.

Effect of Damage Coupling on Springback Prediction
The Young's modulus plays an important role in the springback prediction. In this section, the decrease of Young's modulus during the deformation process is considered through the two approaches of including Young's moduli variation according to loading strain as discusses in the last section. The direct evolution of Young's moduli with respect to accumulating plastic strain is defined with tabular form according to Figure 6. Linear interpolation is considered between given points. For the ductile damage model, the material parameters were defined according to Table 2. The corresponding obtained numerically predicted springback angles are shown in Figure 13. These two methods can both well predict the tendency of the influence of pre-stretch on the final springback angles. Together with the increase of the pre-stretch displacement, the springback angles almost increase linearly with the same slope. Also, it can be found that the obtained springback angle with the ductile damage model can be closer to the experimental results. increases when the pre-strain increases. In this figure, we observe an important effect of kinematic hardening fraction on the springback predicted angle. It needs to be highlighted that for high kinematic hardening fraction (t > 0.8) the increase of pre-strain state induces a decrease of springback angle. For these cases, the difference among numerical results is important. The stress contours before springback are given in Figure 12, which clearly present the big effect of kinematic hardening fraction on the final obtained stress contour on the thickness section during the bending process. Through comparing the experimental and numerical obtained springback angles, it is recommended that the kinematic hardening fraction shall be defined between t = 0.1 and t = 0.5 for Al7050. In the coming section about the ductile damage, the ratio of kinematic hardening in total hardening was assigned to be 0.5.

Effect of Damage Coupling on Springback Prediction
The Young's modulus plays an important role in the springback prediction. In this section, the decrease of Young's modulus during the deformation process is considered through the two approaches of including Young's moduli variation according to loading strain as discusses in the last section. The direct evolution of Young's moduli with respect to accumulating plastic strain is defined with tabular form according to Figure 6. Linear interpolation is considered between given points. For the ductile damage model, the material parameters were defined according to Table 2. The corresponding obtained numerically predicted springback angles are shown in Figure 13. These two methods can both well predict the tendency of the influence of pre-stretch on the final springback angles. Together with the increase of the pre-stretch displacement, the springback angles almost increase linearly with the same slope. Also, it can be found that the obtained springback angle with the ductile damage model can be closer to the experimental results. The damage parameters S will be assigned different values of (0.5, 0.75, 1.0, 1.7, 3.4). In Figure 14, the effect of the damage parameters S on the strain-stress curves is compared. It is found that with the decrease of the damage parameters S, the ductility of the material highly increased. In Figure 15, the ductile damage contour across the thickness zone after the bending process with pre-strains of 4 mm was plotted and compared. When S is smaller, more ductile damage appears on the section, and the springback angle is smaller due to the smaller stress caused by damage. The comparisons with experimental and other numerical obtained results are also given in Figure 15. The damage parameters S will be assigned different values of (0.5, 0.75, 1.0, 1.7, 3.4). In Figure  14, the effect of the damage parameters S on the strain-stress curves is compared. It is found that with the decrease of the damage parameters S, the ductility of the material highly increased. In Figure 15, the ductile damage contour across the thickness zone after the bending process with pre-strains of 4 mm was plotted and compared. When S is smaller, more ductile damage appears on the section, and the springback angle is smaller due to the smaller stress caused by damage. The comparisons with experimental and other numerical obtained results are also given in Figure 15.

Conclusions
In this study, the springback phenomena of AA7055 alloy sheet metal under three points bending test after pre-stretch were investigated. Under the same bending depth, five different prestretch displacements were conducted on the samples, and the obtained springback angles were discussed by comparing the experimental and numerical obtained results. Meanwhile, in the numerical aspect, the influence of the kinematic hardening and ductile damage on the springback was compared and discussed. The conclusions in this study are given as follows:  Through the experimental observation, it was found that under the same bending depth the springback angle increases with the increase of the pre-stretch displacement. The predeformation plays an important role in the final sample shapes.


The influence of the kinematic hardening on the springback prediction of sheet metal is relatively large. The numerically obtained springback angle decreases with the increase of the ratio of kinematic hardening in total hardening. With the increase of pre-stretch displacement, the The damage parameters S will be assigned different values of (0.5, 0.75, 1.0, 1.7, 3.4). In Figure  14, the effect of the damage parameters S on the strain-stress curves is compared. It is found that with the decrease of the damage parameters S, the ductility of the material highly increased. In Figure 15, the ductile damage contour across the thickness zone after the bending process with pre-strains of 4 mm was plotted and compared. When S is smaller, more ductile damage appears on the section, and the springback angle is smaller due to the smaller stress caused by damage. The comparisons with experimental and other numerical obtained results are also given in Figure 15.

Conclusions
In this study, the springback phenomena of AA7055 alloy sheet metal under three points bending test after pre-stretch were investigated. Under the same bending depth, five different prestretch displacements were conducted on the samples, and the obtained springback angles were discussed by comparing the experimental and numerical obtained results. Meanwhile, in the numerical aspect, the influence of the kinematic hardening and ductile damage on the springback was compared and discussed. The conclusions in this study are given as follows:  Through the experimental observation, it was found that under the same bending depth the springback angle increases with the increase of the pre-stretch displacement. The predeformation plays an important role in the final sample shapes.


The influence of the kinematic hardening on the springback prediction of sheet metal is relatively large. The numerically obtained springback angle decreases with the increase of the ratio of kinematic hardening in total hardening. With the increase of pre-stretch displacement, the

Conclusions
In this study, the springback phenomena of AA7055 alloy sheet metal under three points bending test after pre-stretch were investigated. Under the same bending depth, five different pre-stretch displacements were conducted on the samples, and the obtained springback angles were discussed by comparing the experimental and numerical obtained results. Meanwhile, in the numerical aspect, the influence of the kinematic hardening and ductile damage on the springback was compared and discussed. The conclusions in this study are given as follows: • Through the experimental observation, it was found that under the same bending depth the springback angle increases with the increase of the pre-stretch displacement. The pre-deformation plays an important role in the final sample shapes.

•
The influence of the kinematic hardening on the springback prediction of sheet metal is relatively large. The numerically obtained springback angle decreases with the increase of the ratio of kinematic hardening in total hardening. With the increase of pre-stretch displacement, the influence of kinematic hardening will be enlarged. Meanwhile, though this methodology, the ratio of kinematic hardening in total hardening can be coarsely determined.

•
The ductile damage has a great influence on the springback prediction of sheet metal. It not only affects the Young's modulus but the mechanical state of the formed samples. The numerical simulation with consideration of the ductile damage gives more accurately predicted springback angles.