Research on springback characteristic in flexible multi-points 3D stretch-bending process

Abstract

The problem of springback is one of the most significant factors affecting the forming accuracy for aluminum 3D stretch-bending parts. In order to achieve high-efficiency and high-quality forming of such kind of structural components, the springback behaviors of the AA6082 aluminum profiles are investigated based on the flexible multi-points 3D stretch-bending process (3D FSB). Firstly, a finite element simulation model for the 3D FSB process was developed to analyze the forming procedure and the springback procedure. The forming experiments were carried out for the rectangle-section profile to verify the effectiveness of the simulation model. Secondly, the influence of tension on springback was studied, which include the pre-stretching and the post-stretching. Furthermore, the influences of the bending radius and bending sequence are revealed. The results show that: (1) The numerical model can be used to evaluate the effects of bending radius and process parameters on springback in the 3D FSB process effectively. (2) The pre-stretching has little effect on the horizontal springback reduction, but it plays a prominent role in reducing the springback in the vertical direction. (3) The increase of bending deformation in any direction will lead to an increase of springback in its direction and reduce the springback in the other direction. Besides, it reduces the relative error in both directions simultaneously. This research established a foundation to achieve the precise forming of the 3D stretch-bending parts with closed symmetrical cross-section.

Keywords

Stretch-bending springback flexible forming multi-points forming aluminum alloy

Introduction

The aluminum profile stretch-bending process is one of the most widely used techniques for forming spatial structural parts in the manufacturing industry because of its high forming efficiency and excellent repeatability.¹ Generally, a solid die is designed for the manufacturing part with a specified shape in the conventional two-dimensional (2D) stretch-bending process. However, considering the improvement of aerodynamics, lightweight, and safety properties, modern vehicles’ structural design tends to be more complicated. More and more three-dimensional (3D) curved structural parts are designed and applied in the high-end equipment manufacturing industry.² It can provide a lighter, stiffer, and more space-saving skeleton structure than the 2D bending part.³ Therefore, the study of how to achieve accurate stretch-bending for this kind of high-performance component always be a hot issue in the metal forming field.

Due to the restriction of the solid die structure, the conventional stretch-bending process could not realize 3D stretch-bending for the aluminum profile with a complex cross-section. It needs to be improved for its low flexibility and reform its structure to achieve multi-dimensional stretch-bending. Since metal forming’s flexibility plays an essential role in customized mass and individualized production in the manufacturing industry, various flexible forming processes are proposed in recent years.⁴ A typical one is the 3D free bending process, which can be used to produce complex tubes and profiles with asymmetrical cross-sections and continuously-varying radii.^5–7 To realize the flexible forming of the conventional stretch-bending process, a new type of flexible multi-points 3D stretch-bending technology (3D FSB) has been developed.⁸ As shown in Figure 1(a) and (b), based on the idea of multi-points forming technology,⁹ the solid die in the traditional 2D stretch-bending process is replaced by the multi-point dies (MPD). The 3D forming was achieved by the superimposition of the horizontal bending and the vertical bending. The reconfigurable die surface makes the 3D FSB process efficient and flexible for 3D curved parts. Moreover, it is cost-effective for small lot production. The forming equipment is shown in Figure 1(c). Based on the 3D FSB process, a typical 3D curved structure part was precisely formed and employed in the Chinese CRH380 series high-speed train.¹⁰ However, due to the coupling effects of multi-direction force and non-uniform deformation characteristics, it is found that the springback problem of the 3D stretch-bending parts is particularly prominent in the experiment.¹¹ Especially when forming a new part, it needs several times to adjust the envelope shape of the MPD to achieve precise forming, which caused such a waste of time, workforce, and material resources.

Figure 1.

The forming principle of 3D FSB process: (a) the horizontal bending procedure, (b) the vertical bending procedure, and (c) the 3D FSB forming equipment.

The springback problem always is the main factor affecting the forming accuracy in the conventional stretch-bending process.^12–14 The research on the springback mainly focuses on the two directions: (1) the prediction of springback, which is primarily based on the influence of the material constitutive and the bending parameters¹⁵; (2) The method of springback control and compensation, that is, how to realize the precision forming.^16,17 For 2D stretch-bending, Clausen et al.¹⁸ proposed that the most critical factor affecting springback is the material’s strain hardening property and tensile force. The effects of tension, internal pressure, and friction on section distortion, springback, and net extension rate were systematically studied by Miller et al.¹⁹ The parameters design strategy for the 2D stretch-bending process was obtained. Besides, a theoretical model to predict the 2D springback and section distortion by using the incremental strain theory was established.²⁰ Based on this, Corona²¹ presented a formulation that can predict springback with arbitrary thin-wall section shape. The models used to calculate the springback and residual stress of 2D stretch-bending with U-section and T-section were introduced by El-Domiaty and Elsharkawy.^22,23 Yu and Li²⁴ developed the springback prediction model for 2D stretch-bending of L-section aluminum profile. Ruixue et al.²⁵ divided the 2D stretch-bending into 12 cases and established the analytical model to predict the springback for rectangular and U-section profiles. Paulsen and Welo^26,27 investigated the 2D deformation characteristic for rectangular single- and double-chamber AA6060 aluminum profiles with three different wall thicknesses, furthermore, a new, flexible rotary stretch bending machine, and analytical models to predict the springback are developed.²⁸ In addition, the 2D springback behaviors of Al-Li alloys profiles with Z-section and T-section were investigated by Liu et al.²⁹

The precision forming of the 2D stretch-bending parts is well evaluated in the previous research. However, there are relatively few systematic studies on the springback characteristics for 3D parts. With the increase of the forming dimension, the springback shows different characteristics in two different directions. The influence of the 3D bending parameters on springback becomes more complex and needs to be revealed urgently. Therefore, based on the 3D FSB process, in this paper, the springback characteristic of the 3D stretch-bending part is comprehensively and systematically studied by numerical simulation and experiment. The influence law of process parameters and geometric parameters on the 3D springback is discussed, which lays a foundation for accurate forming of the 3D parts with closed symmetrical cross-section.

Numerical simulation models

Based on the 3D FSB technology principle, the numerical simulation models for the forming procedure and the springback procedure were established using commercial finite element software ABAQUS. In order to compare the springback between the horizontal and the vertical direction, the typical symmetrical rectangle-section profile was adopted. The cross-section size of the profile is 40 mm × 40 mm, the thickness of the wall is 4 mm, and the length is 2.2 m. The shape of the target part is a symmetrical circular arc with a radius of 3.82 m and a central angle of 30° both in the horizontal and vertical directions. Considering that the forming process is a highly nonlinear quasi-static process, the dynamic explicit algorithm was used to simulate the forming procedure, while the static implicit algorithm was adopted in the analysis of the springback procedure.

Material model

The widely used AA6082 aluminum is applied in the analysis of the springback characteristics for 3D parts. Table 1 gives the parameters of this material, which are acquired by the uniaxial tensile test. According to the standard of GB/T 228.1-2010, the gauge size of a standard sample is 75.73 mm × 20 mm with a material thickness of 2 mm. The simulation model is established based on the following hypothesis: (1) the material acts agree the Von Mises yield criterion, and its plastic properties follow the Prandtl–Reuss flow rule; (2) The material performs isotropic in the strain-hardening model; (3) For small elastic deformation, the stress satisfies Hooke’s law. The AA6082 constitutive relationship can be express as:

\begin{array}{l} \begin{matrix} σ = E ε & for & σ \end{matrix} < σ_{s} \\ \begin{matrix} σ = K ε^{n} & for & σ \end{matrix} \geq σ_{s} \\ \begin{matrix} σ = - K {| ε |}^{n} & for & σ \end{matrix} \leq - σ_{s} \end{array}

(1)

Where $K$ is the strength coefficient, $n$ is the exponential strain-hardening parameter, $σ_{s}$ is the yield strength. They can be calculated according to the tensile test, which is shown in Table 1.

Table 1.

The AA6082 mechanical properties.

Parameter	Value
Elastic modulus $E$ (GPa)	71.32
Yield stress $σ_{s}$ (MPa)	155.44
Poisson’s ratio $υ$	0.35
Strength factor $K$ (MPa)	398.87
Strain-hardening $n$	0.154

FE simulation model

Figure 2 exhibits the developed simulation model. According to the forming equipment’s actual geometrical dimension, there are nine flexible fundamental units (FFU) employed in the forming process. Since the target part is symmetrical, a 1/2 geometrical model is used to improve calculation efficiency. Assume that only the workpiece is deformed in the forming process, the shell element S4R is adopted to mesh the workpiece, and the rigid element R3D4 is used to mesh the MPD. There are 1920 and 160 elements meshed on the workpiece and each MPD. The elements’ size is 18.33 mm × 5 mm and 6.25 mm × 5 mm for the workpiece and the MPD. The friction coefficient is 0.13 between profile and tools, which is tested by experiment.

Figure 2.

Finite element simulation model for the 3D FSB process.

The 3D FSB forming is a result of multi-body motion and force coupling. In order to reduce the control complexity of the boundary conditions, the FFU and clamps are reasonably simplified. Under the displacement loading control mode, the clamp needs to have 4 degrees of freedom (DOF), as shown in Figure 3(a). They are the movement along the x-axis and z-axis, the rotation around the y-axis and the z-axis. The clamp’s motion trajectory determines the final geometric shape of the formed part. Therefore, as shown in Figure 2, the clamp is simplified as a reference point. It connects with the workpiece by the coupling connection unit. By controlling the reference point’s displacement, the workpiece can be stretched and bent as the actual clamp motion trajectory.

Figure 3.

Structure diagram of the 3D stretch-bending manipulator and the FFU: (a) The 3D stretch-bending manipulator and (b) the flexible fundamental unit.

In the 3D FSB process, the MPD needs to move and rotate with the workpiece, as shown in Figure 3(b), which possesses 4 DOF. They are the movement and rotation on the y-axis and z-axis. Once the workpiece contacts the MPDs, it will be wrapped by the MPDs. Since the MPD could not rotate around the axis of the profile (x-axis), the torsion of the profile will be constrained. Therefore, the torsion problem will not occur during the forming process. The envelope shape of the MPDs in the horizontal direction is adjusted before the bending procedure, while the height control screw constrains the shape in the vertical direction. The axial connect unit is used to express the height control screw to simplify the FFU. Since the MPDs’ envelope is the equidistant curve to the characteristic line of the formed part, the adjustment position parameters of each FFU can be calculated according to the geometric shape of the target part.

Simulation results of 3D FSB process

Generally, the forming process for a 3D curved part can be divided into three steps. Firstly, the workpiece is stretched to the yield point $σ_{s}$ in the pre-stretching step. Secondly, the clamps drive the workpiece to be bent in the horizontal direction. Then, the third step is bending in the vertical direction. The simulation is consistent with the actual forming process. The distribution of stress in this simulation is shown in Figure 4. The maximum stress is 151.02 MPa after pre-stretching, the corresponding maximum strain is 0.00197; after horizontal bending, the maximum stress reached 218.43 MPa, the maximum strain reached 0.03937; then the final maximum stress is 221.61 MPa, the final maximum strain is 0.05425 after vertical bending. It can be seen in Figure 4, the geometric shape of the 3D bending part is smooth and regular, which is close-fitted with the multi-point dies. The variation of stress and strain is continuous.

Figure 4.

Simulation results of FSB forming process: (a) stress state after pre-stretching, (b) stress state after horizontal bending, and (c) stress state after vertical bending.

Springback prediction

Figure 5(a) shows the springback of a 3D stretch-bending part. The springback $Δ α_{i}$ is the dimensional error measured at each MPD, the subscript $i (i = 1, 2, 3 . . . . . N)$ is the MPD’s serial number. The relative error $δ_{i}$ is calculated to analyze the relationship between the formed part and the target part. Besides, the contour error $Δ ϕ$ is defined as the part’s overall forming accuracy. They are expressed mathematically as:

Δ α_{i} = s_{i} - d_{i}

(2)

δ_{i} = \frac{Δ α_{i}}{d_{i}}

(3)

Δ ϕ = \frac{1}{N} \sqrt{\sum_{i = 1}^{N} (\frac{Δ {α_{yi}}^{2} + Δ {α_{zi}}^{2}}{{d_{yi}}^{2} + {d_{zi}}^{2}})}

(4)

Where $s_{i}$ and $d_{i}$ are the coordinates, representing the shape of the 3D stretch-bending part after springback and the shape of the target part.

Figure 5.

Springback prediction results: (a) the definition of springback, (b) the change of the maximum stress, and (c) the distribution of stress before and after springback.

Based on the simulation results, a springback prediction model is established using the ABAQUS/Standard module. In this model, all boundaries defined in the forming process are removed. Only one end of the workpiece is fixed constrained, which is the symmetrical center of the whole actual profile. The formed part’s stress state, mesh, and geometry shape are imported into the model using the predefined field function in ABAQUS. The profile that is constrained only at one end will move freely, resulting in springback deformation. Figure 5 shows the prediction results of springback. The shape deviations occurred both in these two directions. It can be seen that the vertical bending hardly changed the maximum stress in Figure 5(b), and most of the stress is released during the springback procedure. The maximum stress reduced to 48.42 MPa after springback. The springback at the edge of the bent part in the y- and z-axis is 4.04 and 2.48 mm. As shown in Figure 7(c), the stress concentration phenomenon appears after forming, the high-level stress is mainly concentrated on the upper and two side walls of the rectangular section. By comparing with the stress distribution after springback, the upper wall’s residual stress is better released than that on the two side walls.

Experimental verification

The 3D stretch-bending experiments for the rectangle-section profile are conducted to verify the effectiveness of the simulation. In the experiment, the critical process parameters, such as the material, the shape of the target part, the MPD’s adjustment parameters, and the clamps’ trajectories, are in agreement with the FE simulation model. The picture of the experimental part is compared with the simulated part, which is shown in Figure 6.

Figure 6.

Comparison of the forming parts between simulation and experiment.

The geometric shape of the simulation result is consistent with that of the experiment. The dimensional error $Δ α_{i}$ at each MPD is measured by a multi-point springback detector, which can acquire the deviations’ data to the target shape. As shown in Figure 7(a), the simulation and experimental results are compared by extracting the displacement field data from the simulation output. Since the target part is symmetrical, the x-axis is the distance from each measured point to the center of the profile, the vertical axial on the left and right are the coordinates for the horizontal and the vertical springback, respectively.

Figure 7.

The springback simulation results: (a) comparison between simulation and experiment and (b) the ratio error of springback prediction.

By comparing the experimental data with the simulation data, it can be seen that the numerical model has successfully predicted the trend of springback. The maximum prediction error of springback is 0.52 and 0.36 mm in the horizontal and vertical direction, which occurred all at the outermost of the No.5 MPD. The relative error of predicted springback at each MPD is shown in Figure 7(b). The maximum relative error is 14.77%, while the average value is 11.71%. The numerical simulation results meet the springback prediction requirements for the 3D FSB process, which lays a foundation for the study of the following forming characteristics.

Springback characteristic for 3D stretch-bending parts

Effects of tension on springback

A load of axial tension plays an essential role in the stretch bending process. The axial tension is used to stretch the workpiece into the plastic state, which is an effective manner to reduce the springback. Also, it can help to prevent wrinkling caused by the compression. However, the excessive axial tensile force will lead to severe deformation of the cross-section. In this section, the effects of tension on springback in the 3D FSB process are studied. According to the different applied stages of tension, it can be divided into pre-stretching and post-stretching. They are discussed separately.

Effects of pre-stretching

Pre-stretching is the axial tension applied to the profile before the bending process. In order to acquire the influence law of pre-stretching on springback, the pre-stretching takes 0%, 20%, 40%, 60%, 80%, 100%, 120%, and 140% $ε_{s}$ , respectively, in which $ε_{s}$ is the yield strain. The simulation is conducted under the condition that no post-stretching is applied with the same target forming shape.

Figure 8(a) and (b) demonstrate the pre-stretching effects on the horizontal and vertical springback, respectively. With the pre-stretching increase, the horizontal springback does not change obviously, while the vertical springback has a downward trend and varies significantly. Without the post-stretching, the vertical bending is performed on the curved surface after horizontal bending, which results in the uneven distribution of the stress on the formed part. After springback, its residual stress causes the non-monotonic increase of the vertical springback. Figure 8(c) shows the bilinear fitting curves of the maximum springback in the horizontal and vertical directions. Their trends are analyzed separately for below and above 100% $ε_{s}$ pre-stretching. It can be seen that the value of the maximum horizontal springback keeps in a stable interval from 4.43 to 4.95 mm, and it does not change dramatically with the yield state of the material. As the pre-stretching increased from 0% to 100% $ε_{s}$ , the maximum vertical springback was reduced from 5.22 to 2.15 mm. In addition, when the pre-stretching reaches 100% $ε_{s}$ (yield point), the vertical springback tends to be increasing. It can be concluded that the pre-stretching has little effect on the horizontal springback, but it plays a distinct role in reducing the vertical springback. The reason is that the vertical bending has the function of reinforcing and stretching for the horizontal forming in the process of 3D bending. Therefore, when the vertical bending occurs, the horizontal springback will stabilize in a fixed interval, showing no noticeable change. Since the vertical bending is the last step of the whole forming process, its springback mainly depends on the pre-stretching step. That is why pre-stretching has an apparent influence on the springback in the vertical direction.

Figure 8.

The effects of the pre-stretching on the springback: (a) the horizontal direction, (b) the vertical direction, and (c) the fitting curves of the maximum springback.

Effects of post-stretching

The post-stretching is the additional axial load after the bending procedure. In order to clarify its effect on springback, the 0%, 20%, 40%, 60%, 80%, and 100% $ε_{s}$ post-stretching are adopted to analysis the springback characteristics under the 0% and 100% $ε s$ pre-stretching condition, respectively.

The simulation results are shown in Figure 9, in which Figure 9(a) and (b) are the horizontal and vertical springback with no pre-stretching. With the increase of the post-stretching value, there is almost no change for the horizontal springback, while the vertical springback is visibly declined, but it does not continue to decrease with the further increase of the post-stretching. The application of the post-stretching helps to retain the plastic deformation. However, according to the simulation results, when the post-stretching reaches 60% $ε_{s}$ or more, its parameters cause severe cross-section deformation.

Figure 9.

The effects of the post-stretching on the springback: (a) the horizontal direction (without pre-stretching), (b) the vertical direction (without pre-stretching), (c) the horizontal direction (100% $ε_{s}$ pre-stretching), (d) the vertical direction (100% $ε_{s}$ pre-stretching), and (e) the fitting curves with post-stretching.

Figure 9(c) and (d) are the curves of the horizontal and vertical springback under the condition of 100% $ε s$ pre-stretching. In this situation, the horizontal springback tends to decrease, and the maximum springback reduces from 4.75 to 4.19 mm, but the overall change is not apparent. The vertical springback is similar to the horizontal springback. However, as shown in Figure 9(d), the residual stress in the profile’s middle position could not be released without the post-stretching, which helps to minimize the vertical springback. When the post-stretching is applied, the residual stress is released, which promotes the increase of the vertical springback. Similarly, when the post-stretching is more than 60% $ε_{s}$ , the section deformation becomes more serious.

In the condition with post-stretching, the linear fitting curves of the maximum springback with the increase of post-stretching ratio are shown in Figure 9(e). Except for the horizontal fitting curve with no pre-stretching, the other three curves show a downward trend. Whether pre-stretching is applied or not, the maximum horizontal springback does not show noticeable change, while the maximum vertical springback demonstrates apparent decline. Based on the above simulation results and analysis, we can conclude that the post-stretching tension application can effectively release the residual stress and reduce the springback. However, excessive post-stretching makes the section deformation serious. So in the actual production process, <40% $ε_{s}$ of the post-stretching should be selected as the most suitable parameter.

Effects of bending radius

Since the springback is closely related to the bending deformation, in this section, the tension parameters of 100% $ε_{s}$ pre-stretching and 40% $ε_{s}$ post-stretching are adopted to investigate the effects of the bending radius on the springback.

The effect of horizontal bending radius

Under the same load condition, the simulation was conducted with the same vertical bending radius 3.82 m and different horizontal bending radius 3.82, 2.86, 2.29, and 1.91 m. The simulation results are shown in Figure 10, in which Rh is the horizontal bending radius and Rv is the vertical bending radius. As shown in Figure 10(a), the springback in the horizontal direction increases with the decrease of the horizontal bending radius, and the maximum horizontal springback rises from 4.5 to 6.34 mm. On the contrary, the vertical springback tends to decline, which is shown in Figure 10(b), and the maximum vertical springback deduces from 4.49 to 3.32 mm. The decrease of bending radius leads to the increase of horizontal deformation, so the increased plastic deformation accumulates more elastic deformation to be released in the horizontal direction, resulting in the raising of the horizontal springback. However, since the shape of the deformation in the vertical direction is the same, the plastic deformation in the horizontal direction strengthens the pre-stretching effect before the vertical bending, which makes the vertical forming is more sufficient, so the springback is reduced.

Figure 10.

The effects of the horizontal bending radius Rh on the springback: (a) the horizontal springback with different Rh, (b) the vertical springback with different Rh, (c) the horizontal springback ratio error, and (d) the vertical direction springback ratio error.

According to equation (3), the relative error can be calculated at each MPD. So the influence of the deformation degree on the springback can be eliminated. It can be seen from Figure 10(c) that the horizontal relative error at different MPD only has slight variations at each horizontal bending radius, which shows that the deformation degree of each position on the profile is stable. In addition, although the horizontal springback is on the rise, the horizontal relative error shows a decreasing trend with the increase of the deformation. It reduced from 2.82% with Rh = 3.82 m to 1.98% with Rh = 1.91 m at the No. 5 MPD, which indicates that the increase of the horizontal bending plastic deformation suppresses the tendency of the springback to increase in the same proportion. As can be seen from Figure 10(d), the vertical relative error tends to decrease as the bending radius is reduced in the horizontal direction, it declined from 2.92% with Rh = 3.82 m to 2.05% with Rh = 1.91 m at the No. 5 MPD. Under each deformation condition, the vertical relative error at each MPD is the same.

The effect of the vertical bending radius

In this section, the horizontal bending radius was kept at 3.82 m, and the tensile process parameters unchanged. The test is conducted with the vertical bending radius is 3.82, 2.86, 2.29, and 1.91 m, respectively. The simulation results are shown in Figure 11. Since the vertical bending has a post-stretching effect on the horizontal forming, the more deformation occurred in the vertical direction, the more horizontal springback decreased. As shown in Figure 11(a), the maximum horizontal springback reduces from 4.50 to 2.49 mm. Although the target shape of the vertical direction changed considerably, the vertical springback does not vary significantly, and the maximum vertical springback only rises from 4.49 to 5.55 mm.

Figure 11.

The effects of the vertical bending radius Rv on the springback: (a) the horizontal springback with different Rv, (b) the vertical springback with different Rv, (c) the horizontal relative error, and (d) the vertical relative error.

As shown in Figure 11(c) and (d), the relative error in the horizontal and vertical directions decreases with the increase of the vertical deformation. The horizontal and vertical relative error of the workpiece is stable under each deformation condition. Besides, take the relative error at the No. 5 MPD as an example, the relative error in the horizontal and vertical directions are similar in any forming conditions. Its values are related to the deformation degree of the whole workpiece, which is illustrated in Table 2.

Table 2.

The relative error at the No. 5 MPD.

Bending condition		Horizontal relative error (%)	Vertical relative error (%)
Rv = 3.82m	Rh = 3.82 m	2.82%	2.92
	Rh = 2.86 m	2.28	2.27
	Rh = 2.29 m	2.11	2.17
	Rh = 1.91 m	1.98	2.05
Rh = 3.82 m	Rv = 2.86 m	2.24	2.15
	Rv = 2.29 m	1.71	2.05
	Rv = 1.91 m	1.57	1.66

Discussion

In order to achieve the high-quality forming of the 3D stretch-bending parts, the parameters design strategy of the 3D FSB process should be considered from two aspects: the tensile parameters and the geometric parameters. In terms of tensile parameters, the method of pre-stretching and post-stretching should be used together, the process windows for pre-stretching and post-stretching are shown in Figures 8(c) and 9(e). The pre-stretching range from 80% to 120% $ε_{s}$ , and the post-stretching range from 20% to 40% $ε_{s}$ is a reasonable scale for the profile with a closed symmetrical cross-section. Figure 12 shows the influence of tension on the stress distribution. According to the comparison of stress distribution between Case 1 (no pre-stretching and no post-stretching (Figure 12(a)) and Case 2 (100% $ε_{s}$ pre-stretching and no post-stretching (Figure 12(b)), the application of pre-stretching reinforces the plastic deformation of the upper and lower walls of the cross-section, which helps to reduce the vertical springback and the maximum residual stress. In the loading condition without post-stretching, considerable residual stress is generated on the two side walls of the cross-section, both in these two cases. Once the post-stretching is applied, as shown in Figure 12(c) (Case 3: 100% $ε_{s}$ pre-stretching and 40% $ε_{s}$ post-stretching), the residual stress distribution is more uniform, which improves the service performance of the formed part. However, due to the residual stress on the two side walls is effectively released, the springback in the vertical direction increases.

Figure 12.

The effect of the tension on the distribution of stress: (a) case 1: no pre-stretching and no post-stretching, (b) case 2: 100% $ε_{s}$ pre-stretching and no post-stretching, and (c) case 2: 100% $ε_{s}$ pre-stretching and 40% $ε_{s}$ post-stretching

For the geometric parameters, the effect of the bending radius on the max horizontal and vertical springback is shown in Figure 13. It can be concluded that: if the deformation is aggravated in any one direction, the springback in this direction will be increased, while the springback in the other direction will be reduced. Besides, if we convert the springback into the rate-related parameters, we can see that the relative error is always <3%. The relative error does not increase with the rise of the center distance. Its value is maintained as a relatively stable parameter at each MPD. The increase of bending deformation in any one direction will reduce the relative error both in the horizontal and vertical directions at the same time.

Figure 13.

The effects of the bending radius.

Since the workpiece is a center symmetry profile with a rectangular cross-section, the final shape is the same for the forming condition at Case 1: Rh = 3.82 m, Rv = 1.91 m, and Case 2: Rh = 1.91 m, Rv = 3.82 m. The stress distribution of these two cases is shown in Figure 14. The simulation results are displayed in mirror symmetry for comparison. Due to the different bending sequences, the formed parts show different springback under the same tensile parameters. The maximum stress in Case 1 is greater than that in Case 2, both before and after springback. Since the residual stress of Case 1 is not released well, the contour error of the formed part is less than that of Case 2. There is no obvious forming defect occurred in Case 1 and Case 2. According to equation (4), the contour error is calculated in Table 3. Without considering the section deformation, the contour error of the Rh < Rv group is always larger than that of the Rh > Rv group. Therefore, if the cross-section of the formed workpiece is symmetrical, it can be concluded that the bending sequence should be arranged as small deformation in the horizontal direction and the large deformation in the vertical direction to further reduce the springback.

Figure 14.

The effect of the bending sequence on the distribution of stress: (a) case 1: horizontal 3.82 m, vertical 1.91 m and (b) case 2: horizontal 1.91 m, vertical 3.82 m.

Table 3.

The comparison of the contour error with the same forming shape.

Forming shape (radius at two directions)	Rh < Rv (%)	Rh > Rv (%)
3.82 m, 2.86 m	0.91	0.84
3.82 m, 2.29 m	0.87	0.85
3.82 m, 1.91 m	0.84	0.76

Conclusion

The 3D curved structural parts are more and more designed and used in the high-end equipment manufacturing industry. In this paper, extensive numerical simulations for the 3D FSB process have been performed to determine the springback characteristics.

According to the principle of the 3D FSB technology, the numerical simulation models were established. The simulation results of the rectangle-section profile are in good agreement with the experimental results. The prediction error of springback is <0.52 mm in the horizontal direction and 0.36 mm in the vertical direction. The maximum predicted relative error is 14.77%, while the average value is 11.71%.

The pre-stretching has little effect on the horizontal springback reduction, but it plays a prominent role in reducing the springback in the vertical direction. Besides, the application of the post-stretching can effectively release the residual stress and reduce the springback. However, excessive post-stretching makes the section deform severely. In the actual production process, 100% $ε s$ pre-stretching and less than 40% $ε s$ post-stretching should be selected as the most suitable parameters for the profile with a closed symmetrical cross-section.

For the effect of bending radius, the increase of bending deformation in any direction will lead to an increase of springback in its direction and reduce the springback in the other direction. Moreover, it will reduce the relative error both in these two directions simultaneously. The relative error is maintained as a relatively stable parameter which is <3% at each MPD. The bending sequence also affects the forming precision. The large deformation should be arranged in the vertical direction further to reduce the springback in the 3D FSB process.

Footnotes

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the National Natural Science Foundation of China (No.51805045) and Scientific and Technological Developing Scheme of Ji Lin Province (Nos. 2019 0302100GX, 20200401115GX).

ORCID iD

Song Gao

Availability of data and materials

The data during the current study are available from the corresponding author on reasonable request.

References

Vollertsen

Sprenger

Kraus

, et al. Extrusion, channel, and profile bending: a review. J Mater Process Technol 1999; 87(1–3): 1–27.

Guo

Chen

, et al. Simulation and experimental research of the free bending process of a spatial tube. J Mater Process Technol 2018; 255: 137–149.

Chatti

Hermes

Tekkaya

, et al. The new TSS bending process: 3D bending of profiles with arbitrary cross-sections. CIRP Ann Manuf Technol 2010; 59(1): 315–318.

Yang

Bambach

Cao

, et al. Flexibility in metal forming. CIRP Ann Manuf Technol 2018; 67(2): 743–765.

Cheng

Guo

Tao

, et al. Investigation of the effect of relative thickness (t/d) on the formability of the aa6061 tubes during free bending process. Int J Mech Sci 2019; 160: 103–113.

Guo

Chen

, et al. Simulation and experimental research of the free bending process of a spatial tube. J Mater Process Technol 2018; 255: 137–149.

Groth

Frohn

Engel

Product planning system for manufacturing-oriented modeling of freeform bend tubes produced by three-roll-push-bending. Procedia Manuf 2019; 34: 10–18.

Liang

J-C

Gao

Teng

, et al. Flexible 3D stretch-bending technology for aluminum profile. Int J Adv Manuf Technol 2014; 71(9–12): 1939–1947.

Cai

Sui

, et al. Multi-point forming technology for sheet metal. J Mater Process Technol 2002; 129(1–3): 333–338.

10.

Gao

Liang

, et al. Precision forming of the 3D curved structure parts in flexible multi-points 3D stretch-bending process. Int J Adv Manuf Tech 2018; 95(1–4): 1205–1213.

11.

Liang

Lin

, et al. Springback study on profile flexible 3d stretch-bending process using the neural network. Adv Mater Sci Eng 2019; 2019: 6465196.

12.

Liu

Wang

, et al. Minimum bending radius of Al-Li alloy extrusions in stretch bending. Proc Inst Mech Eng Part B J Eng Manuf 2016; 232(2): 281–295.

13.

Choudhury

Ghomi

Springback reduction of aluminum sheet in V-bending dies. Proc Inst Mech Eng Part B J Eng Manuf 2014; 228: 917–926.

14.

Liu

Wang

, et al. Springback behaviors of extruded 6063 aluminum profile in subsequent multi-stage manufacturing processes. Int J Adv Manuf Technol 2020; 109(1): 1–13.

15.

Wagoner

Lim

Lee

M-G

. Advanced issues in springback. Int J Plast 2013; 45: 3–20.

16.

Zhao

Liu

Zhang

, et al. Optimal design of process parameters of rotary-draw bending process for thin-walled rectangular tube of aluminum alloy. Proc Inst Mech Eng Part B J Eng Manuf 2014; 228(11): 1442–1448.

17.

Liu

Zhang

, et al. Optimization of post-stretching elongation in stretch bending of aluminum hollow profile. Int J Adv Manuf Technol 2016; 82(9–12): 1737–1746.

18.

Clausen

Hopperstad

Langseth

. Sensitivity of model parameters in stretch bending of aluminium extrusions. Int J Mech Sci 2001; 43(2): 427–453.

19.

Miller

Kyriakides

Bastard

. On bend-stretch forming of aluminum extruded tubes––I: experiments. Int J Mech Sci 2001; 43(5): 1283–1317.

20.

Miller

Kyriakides

Corona

. On bend-stretch forming of aluminum extruded tubes––II: analysis. Int J Mech Sci 2001; 43(5): 1319–1338.

21.

Corona

. Simple analysis for bend-stretch forming of aluminum extrusions. Int J Mech Sci 2004; 46(3): 433–448.

22.

El-Domiaty

Elsharkawy

. Stretch-bending analysis of U-section beams. Int J Mach Tools Manuf 1998; 38(1–2): 75–95.

23.

Elsharkawy

El-Domiaty

. Determination of stretch-bendability limits and springback for T-section beams. J Mater Process Technol 2001; 110(3): 265–276.

24.

. Theoretical analysis on springback of L-section extrusion in rotary stretch bending process. Trans Nonferrous Met Soc China 2011; 21(12): 2705–2710.

25.

Ruixue

Xuehui

, et al. Stretch bending and springback of profile in the loading method of prebending and tension. Int J Mech Sci 2018; 144: 746–764.

26.

Paulsen

Welo

. Cross-sectional deformations of rectangular hollow sections in bending: Part I––experiments. Int J Mech Sci 2001; 43(1): 109–129.

27.

Paulsen

Welo

. Cross-sectional deformations of rectangular hollow sections in bending: Part II—analytical models. Int J Mech Sci 2001; 43(1): 131–152.

28.

Welo

. Analytical springback assessment in flexible stretch bending of complex shapes. Int J Mach Tools Manuf 2020; 160: 1–19.

29.

Liu

Wang

, et al. Springback analysis of Z & T-section 2196-T8511 and 2099-T83 Al-Li alloys extrusions in displacement controlled cold stretch bending. J Mater Process Technol 2015; 225: 295–309.