J Matdes 2011 03 069
J Matdes 2011 03 069
J Matdes 2011 03 069
a r t i c l e i n f o a b s t r a c t
Article history: Bent parts of complex shapes with high precision are increasingly required. To achieve a high precision of
Received 15 December 2010 parts, especially the required bending angle, a suitable design of process parameters is strictly consid-
Accepted 28 March 2011 ered. In this study, process parameters of bending angle, material thickness and punch radius were inves-
Available online 31 March 2011
tigated. The finite element method (FEM), in association with the Taguchi and the analysis of variance
(ANOVA) techniques, was carried out to investigate the degree of importance of process parameters in
Keywords: V-bending process. The results revealed that the degree of importance of process parameters in V-bend-
Forming
ing process depended on the spring-back and spring-go. The material thickness has a major influence on
Elastic behaviour
Plastic behaviour
the spring-back. In contrast, in the case of spring-go, the bending angle has a major influence and closely
followed by the material thickness. In addition to predicting the degree of importance of process param-
eters by the combination of the FEM simulation, the Taguchi technique, and the ANOVA technique, by
facilitating an improvement in the quality of the required bending angle was strictly considered by opti-
mization of these process parameters corresponding with the spring-back and spring-go.
Ó 2011 Elsevier Ltd. All rights reserved.
0261-3069/$ - see front matter Ó 2011 Elsevier Ltd. All rights reserved.
doi:10.1016/j.matdes.2011.03.069
S. Thipprakmas, W. Phanitwong / Materials and Design 32 (2011) 4430–4436 4431
rial thickness, and punch radius, in relation to the spring-back and Simulation model Plane strain model
spring-go. The commercial software DEFORM-2D was used as the fi- Object type Workpiece:Elasto-plastic, Punch/Die:rigid
nite-element simulation (FE-simulation) code to carry out the pre- Workpiece material A1100-O (rB = 92.5 MPa, k = 46%)
diction of spring-back and spring-go. The experiments were Flow curve equation r ¼ 160:4e0:22
performed to validate the FEM simulation results. The ANOVA tech- Angular punch radius (Rp) 2.5 and 3.5 mm
Angular die radius (Rd) 5 mm
nique is able to specify the process parameters that markedly influ-
Bending angle (h) 30°, 45°, 120°, and 135°
ence the spring-back and spring-go and yield information about the Workpiece thickness (t) 2.0, 3.0, and 4.0 mm
degree of importance of each process parameter on the spring-back Friction coefficient (l) 0.1
and spring-go. The results elucidated that the degree of importance
of each process parameter in bending angle was different in the case
2.2. Experimental procedure
of spring-back and spring-go. The material thickness has a major
influence on the spring-back. The bending angle has a major influ-
Referring to the author’s past researches [13,14], V-bending dies
ence on the spring-go closely followed by the material thickness.
and a 5-ton universal tensile testing machine (Lloyd Instruments
This difference in the process parameter having a major influence
Ltd) were used for the laboratory V-bending experiments with
on the spring-back and spring-go in the V-bending process could
the same FEM simulation conditions, the laboratory experiments
be clearly identified based on the stress distribution analysis. There-
were carried out to validate the FEM simulation results. Aluminum
fore, the achievement of the required bending angle depended on
(A1100-O) of 30 mm width and 3 mm thickness was used as a
the application of the optimal value of the process parameters and
workpiece material. The bending angle was examined using the
their favorable combination relating to the spring-back and
profile projector. The bending force was recorded and compared
spring-go.
with the bending force analyzed by FEM.
2. The FEM simulation, experimental procedures, and Taguchi 2.3. Taguchi method
method
In this study, two levels of the three process parameters, includ-
2.1. FEM simulation procedure ing bending angle, material thickness, and punch radius were ap-
plied. The experimental design of FEM simulation using the two
Fig. 1 shows the FEM simulation model. A half V-bending sim- level for three-parameter full factorial design (23) is shown in Table
ulation model with a die radius (Rd) of 5 mm was used to reduce 2. The spring-back and spring-go were the process responses. In
the calculated time. In this study, the model was a two-dimen- this study, to understand the effects of the process parameters,
sional plane strain 70 mm in length and 3 mm in thickness. The the Taguchi technique employing a generic signal-to-noise ratio
commercial analytical code for two-dimensional implicit quasi- for the ‘‘Smaller is better (SNS)’’ characteristics was considered
static finite element method (DEFORM-2D) was used as the FEM for the amounts of spring-back and spring-go as depicted in the
simulation tool. Automatic remeshing was set at every five steps following equation:
to prevent a divergence calculation due to excessive deformation. !
With reference to earlier research [12–14], Table 1 shows the de- 1X n
SNS ¼ 10 log y2 ð1Þ
tails of FEM simulation conditions and material properties. A n i¼1 i
punch radius (Rp) of 2.5 and 3.5 mm, a material thickness (t) of 2,
3, and 4 mm, and bending angles (h) of 30°, 45°, 120°, and 135° where n represents the number of experiments and y represents the
were investigated. The workpiece material was set as an elasto- process response.
plastic type with a rectangular element of approximately 3500 ele- In addition to the FEM simulation results of the spring-back and
ments. The punch and die were set as a rigid type. The material spring-go, the ANOVA technique was also applied to illustrate the
used was aluminum A1100-O (JIS) and its properties included ten- degree of importance of each process parameter that markedly
sile strength of 92.5 MPa and elongation of 46%. The constitutive influenced the spring-back and spring-go in the V-bending process.
equation was determined from the SS (Stress Strain) curve ob-
tained by the tensile testing experiment [12–14]. 3. Results and discussion
Table 3
Process parameters and their levels.
1X 8
Table 4 S=N ¼ ðS=NÞk ð2Þ
The amounts of spring-back and spring-go analyzed by FEM simulation. 8 k¼1
Experiment Bending Material Punch The amount
no. angle thickness radius of SB and SG/ 8
X 2
(h)/(°) (t)/(mm) (Rp)/(mm) (°) SS ¼ ðS=NÞij S=N ð3Þ
i¼1
(a) Spring-back (SB)
1 30 2.0 2.5 0.58
2 45 2.0 2.5 0.45 3
X 2
3 30 3.0 2.5 0.11 SSi ¼ ðS=NÞij S=N ð4Þ
4 45 3.0 2.5 0.09 j¼1
5 30 2.0 3.5 0.70
6 45 2.0 3.5 0.60
SSi
7 30 3.0 3.5 0.27 % Contributioni ¼ 100 ð5Þ
8 45 3.0 3.5 0.20 SS
(b) Spring-go (SG) The S=N; SS, and percentage contributions calculated for each
1 120 2.0 2.5 0.40
process parameters in the case of spring-back and spring-go are
2 135 2.0 2.5 0.48
3 120 3.0 2.5 0.46 illustrated in Table 5. In the case of spring-back, the percentage
4 135 3.0 2.5 0.54 contributions of the material thickness, the punch radius, and the
5 120 2.0 3.5 0.33 bending angle were 83.6365, 13.9263, and 2.4372, respectively.
6 135 2.0 3.5 0.47 On the other hand, in the case of spring-go, the percentage contri-
7 120 3.0 3.5 0.44
8 135 3.0 3.5 0.52
butions of the material thickness, the punch radius, and the bend-
ing angle were 33.4050, 6.9883, and 59.6067, respectively. Based
on these statistical analysis results, it was found that the degree
of importance of process parameters in V-bending process de-
pended on the spring-back and spring-go. As shown in the results,
good agreement with the bending theory in which the spring-back
in the case of spring-back the process parameters of the material
and spring-go generated in small and large bending angles, respec-
thickness had the most influence, respectively followed by the
tively. Specifically, spring-back was generated in 30° and 45° bend-
small influence of the punch radius and the bending angle. In the
ing angles while the spring-go was generated in 120° and 135°
case of spring-go it was vice versa, with the process parameters
bending angles. To investigate the degree of importance of the pro-
of bending angle having the highest influence, closely followed
cess parameters, therefore, the ANOVA technique was carried out
by the material thickness.
in case by case of spring-back and spring-go. The ‘‘Smaller is better
(SNs)’’ characteristics, as shown in Eq. (1), were considered for the
3.2. Effects of process parameters on the spring-back and spring-go
spring-back and spring-go. The mean value of the overall S/N (S=N)
is expressed as Eq. (2), where, k is the number of experiments. The
Many past researches have been studied on the effects of pro-
sums of squares owing to the variations of the overall mean (SS)
cess parameters in V-bending process, including bending angle
and the mean of the process parameters (SSi) are expressed as
[4,11], material thickness [3,11,15], punch radius [12,15], and
Eqs. (3) and (4), respectively. The percentage contributions as
punch height [13,14]. However, most of these were studied on
shown in Eq. (5) were calculated which was generally applied for
the spring-back. Therefore, the effects of process parameters in
considering the degree of importance of each process parameter
V-bending process on the spring-go were rarely investigated. In
as reported in the literature [23,24].
this study, the effects of bending angle, material thickness, and
Table 5
ANOVA results.
Process parameters Bending angle (h)/(°) Material thickness (t)/(mm) Punch radius (Rp)/(mm)
30 45 2.0 3.0 2.5 3.5
(a) Spring-back
(S/N)ij 9.59 11.57 4.80 16.36 12.94 8.22
Sum of squares (SSi) 1.95 66.81 11.12
% Contributions 2.4372 83.6365 13.9263
(b) Spring-go
120 135 2.0 3.0 2.5 3.5
(S/N)ij 7.87 5.99 7.63 6.23 6.61 7.25
Sum of squares (SSi) 1.76 0.98 0.21
% Contributions 59.6067 33.4050 6.9883
S. Thipprakmas, W. Phanitwong / Materials and Design 32 (2011) 4430–4436 4433
0.20
/ (°)
0.00 Bending angle (θ ) / (°)
-0.20 30 45 120 135
-0.40
-0.60
-0.80
(a) Rp 2.5 mm
Spring-go value Spring-back value
0.80
t = 2 mm
0.60
t = 3 mm
0.40
0.20
/ (°)
-0.40
-0.60
-0.80
(b) Rp 3.5 mm
Fig. 2. Comparison of the amounts of spring-back and spring-go analyzed by FEM simulation with respect to the process parameters.
3.13 mm
4.75 mm
3.30 mm
10.52 mm
(a-1) θ = 30°, Rp = 3.5 mm, t = 2 mm (b-1) θ = 120°, Rp = 3.5 mm, t = 2 mm
4.72 mm
6.02 mm
4.36 mm
11.33 mm
(a-2) θ = 30°, Rp = 3.5 mm, t = 3 mm (b-2) θ = 120°, Rp = 3.5 mm, t = 3 mm
5.91 mm
6.90 mm
4.80 mm
12.87 mm
(a-3) θ = 30°, Rp = 3.5 mm, t = 4 mm (b-3) θ = 120°, Rp = 3.5 mm, t = 4 mm
Fig. 3. Comparison of mean stress with respect to the material thickness (Rd = 5 mm).
4434 S. Thipprakmas, W. Phanitwong / Materials and Design 32 (2011) 4430–4436
punch radius on spring-back and spring-go were clearly investi- material thickness [11,15], and decrease in punch radius [15].
gated and identified based on the stress distribution analysis. However, the spring-go was rarely reported. Only the author’s past
Fig. 2 shows the amounts of spring-back and spring-go analyzed research reported that the amount of spring-go decreased as the
by FEM with respect to bending angle, material thickness, and punch radius increased [12]. However, the mechanism of spring-
punch radius. The results showed the same tendency of spring- back/spring-go related to these process parameters and the degree
back and spring-go in both cases of punch radius of 2.5 and of important of each process parameter have not been clearly iden-
3.5 mm, as shown in Fig. 2a and b, respectively. Specifically, as tified yet. Base on the stress distribution analysis, the mechanism
the bending angle increased, the amount of spring-back decreased of spring-back and spring-go related to process parameters was
and the amount of spring-go increased. In addition, as the material clearly identified. Fig. 3 shows a comparison of mean stress with
thickness increased, the amount of spring-back decreased whereas respect to the material thickness. The results illustrated that as
the amount of spring-go increased. As the punch radius increased, the material thickness increased, the bending allowance zone
the amount of spring-back increased whereas the amount of and reversed bending zone increased. It was observed that the in-
spring-go decreased. The results obtained in this research are gen- crease in reserved bending zone was larger than the increase in
erally in agreement with those reported in the literature. These re- bending allowance zone. In addition, as the material thickness in-
sults were reported by FEM simulation and experiments that creased, the plastic zone across the thickness increased whereas
spring-back decreased with increase in bending angle [4,11] and the elastic zone over the neutral plane was approximately
3.13 mm
3.35 mm
10.52 mm 10.27 mm
(a) θ = 30°, Rp = 3.5 mm, t = 2 mm (b) θ = 45°, Rp = 3.5 mm, t = 2 mm
5.03 mm
5.00 mm
10.75 mm 11.08 mm
(c) θ = 45°, Rp = 2.5 mm, t = 3 mm (d) θ = 45°, Rp = 3.5 mm, t = 3 mm
Fig. 4. Comparison of mean stress with respect to the process parameters in case of spring-back (Rd = 5 mm).
4.75 mm 4.95 mm
3.30 mm 2.90 mm
5.27 mm 5.18 mm
2.87mm 3.00 mm
(c) θ = 135°, Rp = 2.5 mm, t = 3 mm (d) θ = 135°, Rp = 3.5 mm, t = 3 mm
Fig. 5. Comparison of mean stress with respect to the process parameters in case of spring-go (Rd = 5 mm).
S. Thipprakmas, W. Phanitwong / Materials and Design 32 (2011) 4430–4436 4435
Table 6
Required bending angle
Comparison of the amounts of spring-back and spring-go with reference to the 120.0
process parameters. FEM EXP
Fig. 6. Comparison of the bending angle between FEM and experimental results
(h = 120°).
constant. Therefore, with these characteristics, as the material
thickness increased, the amounts of spring-back and spring-go de-
creased and increased, respectively. Fig. 4 shows a comparison of 600
mean stress in case of spring-back, with respect to the bending an-
spring-go. These ANOVA results and the effects of process param- [6] Kim HS, Koc M. Numerical investigations on springback characteristics of
aluminum sheet metal alloys in warm forming conditions. J Mater Process
eters could be clearly identified based on the stress distribution
Technol 2008;204:370–83.
analysis. In the case of spring-back, the material thickness had a [7] Yu HY. Variation of elastic modulus during plastic deformation and its
great effect on the bending allowance and reversed bending zones. influence on springback. J Mater Des 2009;30:846–50.
It was vice versa in the case of spring-go; the bending angle and [8] Kazan R, Firat M, Tiryaki AE. Prediction of springback in wipe-bending process
of sheet metal using neural network. J Mater Des 2009;30:418–23.
subsequent material thickness had a major effect on the bending [9] Ozturk F, Toros S, Kilic S. Tensile and spring-back Behavior of DP600 advanced
allowance and reserved bending zones. Therefore, the process high strength steel at warm temperatures. J Iron Steel Res Int 2009;16:41–6.
parameter design in V-bending process depended on the spring- [10] Parsa MH, Nasher al ahkami S, Ettehad M. Experimental and finite element
study on the spring-back of double curved aluminum/polypropylene/
back and spring-go. Specifically, the material thickness was aluminum sandwich sheet. J Mater Des 2010;31:4174–83.
strongly considered in the case of spring-back. In contrast, the [11] Panthi SK, Ramakrishnan N, Ahmed M, Singh SS, Goel MD. Finite element
bending angle and subsequent material thickness were strongly analysis of sheet metal bending process to predict the springback. J Mater Des
2010;31:657–62.
considered in the case of spring-go. The combination of the FEM [12] Thipprakmas S, Rojananan S. Investigation of spring-go phenomenon using
simulation, the Taguchi method, and the ANOVA technique was finite element method. J Mater Des 2008;29:1526–32.
an effective tool to predict the degree of importance of the process [13] Thipprakmas S. Finite element analysis of punch height effect on V-bending
angle. J Mater Des 2010;31:1593–8.
parameters in the V-bending process, in addition to aiding in [14] Thipprakmas S. Effect of ratio of punch height to workpiece length in partial V-
improvement of the quality of the required bending angle by opti- bending process. Steel Res Int 2010;81:769–72.
mization of the process parameters relating to the spring-back and [15] Bakhshi-Jooybari M, Rahmani B, Daeezadeh V, Gorji A. The study of springback
of CK67 steel sheet in V-die and U-die bending processes. J Mater Des
spring-go.
2009;30:2410–9.
[16] Yanagimoto J, Oyamada K. Mechanism of springback-free bending of
Acknowledgements highstrength steel sheets under warm forming conditions. Ann CIRP
2007;56:265–8.
[17] Kim H, Nargundkar N, Altan T. Prediction of bend allowance and springback in
The presented research is partially supported by a grant from air bending. ASME 2007;129:342–51.
the National University Research Project of Thailand (NRU Project) [18] Garcia-Romeu ML, Ciurana J, Ferrer I. Springback determination of sheet
under Grant No. 54000298, the Thai Research Fund (TRF) under metals in an air bending process based on an experimental work. J Mater
Process Technol 2007;191:174–7.
Grant No. TRF-MRG5380117, and the Center of Excellence in Sheet [19] Panthi SK, Ramakrishnan N, Pathak KK, Chouhan JS. An analysis of springback
and Bulk Metal Forming Technology (SBMFT). in sheet metal bending using Finite Element Method (FEM). J Mater Process
Technol 2007;186:120–4.
[20] Song JH, Huh H, Kim SH. Stress based springback reduction of a channel shaped
References auto-body part with high-strength steel using response surface methodology. J
Eng Mater Technol 2007;129:397–406.
[1] Zang SL, Liang J, Guo C. A constitutive model for spring-back prediction in [21] Dongjuan Z, Zhenshan C, Xueyu R, Yuqiang L. An analytical model for
which the change of young’s modulus with plastic deformation is considered. predicting springback and side wall curl of sheet after U-bending. Comput
Int J Mach Tool Manuf 2007;47:1791–7. Mater Sci 2007;38:707–15.
[2] Leu DK, Hsieh CM. The influence of coining force on spring-back reduction in [22] Chen C, Ko M, Wenner ML. Experimental investigation of springback variation
V-die bending process. J Mater Process Technol 2008;196:230–5. in forming of high strength steels. J Manuf Sci Eng 2008;130:1–9.
[3] Wang J, Verma S, Alexander R, Gauc J. Spring-back control of sheet metal air [23] Thipprakmas S. Application of Taguchi technique to investigation of geometry
bending process. J Manuf Process 2008;10:21–7. and position of V-ring indenter in fine-blanking process. J Mater Des
[4] Tekaslan O, Gerger N, Seker U. Determination of spring-back of stainless steel 2010;31:2496–500.
sheet metal in ‘‘V’’: bending dies. J Mater Des 2008;29:1043–50. [24] Thipprakmas S. parameter design for the counter-blanking process using the
[5] Meinders T, Burchitz IA, Bonte MHA, Lingbeek RA. Numerical product design: finite element method and the Taguchi approach. Key Eng Mater
springback prediction, compensation and optimization. Int J Mach Tool Manuf 2010;443:134–9.
2008;48:499–514.