HOME > Browse Articles > Archive

Archive

Journal of the Korean Society for Precision Engineering - Vol. 36 , No. 9

[ REGULAR ]
Journal of the Korean Society for Precision Engineering - Vol. 36, No. 9, pp. 883-890
Abbreviation: J. Korean Soc. Precis. Eng.
ISSN: 1225-9071 (Print) 2287-8769 (Online)
Print publication date 01 Sep 2019
Received 26 Mar 2019 Revised 21 May 2019 Accepted 29 May 2019
DOI: https://doi.org/10.7736/KSPE.2019.36.9.883

Application/Comparison Study of a Graphical Method of Forming Limit Curve Estimation for DP590 Steel Sheets
The-Thanh Luyen1 ; Quoc-Tuan Pham2 ; Young-Suk Kim2 ; Duc-Toan Nguyen3, #
1Faculty of Mechanical Engineering, Hungyen University of Technology and Education, Hungyen, Vietnam
2School of Mechanical Engineering, Kyungpook National University, 80, Daehak-ro, Buk-gu, Daegu, 41566, Republic of Korea
3School of Mechanical Engineering, Hanoi University of Science and Technology, 1 Dai Co Viet Road, Hanoi, Vietnam

Correspondence to : #E-mail: toan.nguyenduc@hust.edu.vn, TEL: +84-988-693-047


Copyright © The Korean Society for Precision Engineering
This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
Funding Information ▼

Abstract

Since its introduction, the Modified Maximum Force Criterion proposed by Hora et al., has been widely used to theoretically estimate the forming limit curve of metal sheets. On the basis of this criterion, a graphical method was presented in our previous study to simplify the evaluation of the forming limit curve (FLC) of metal sheets. This paper presents an application of the graphical method to estimate the FLCs of an advanced high-strength steel sheet material, DP590. The material is frequently used in the automotive industry. To verify the ability of the graphical method, various hardening laws and yield functions were used to estimate the forming limiting curves for the examined material. The calculated forming limiting curves are then adopted for the finite element method (FEM) to predict the fracture heights of different notch specimens desired by the Hecker’s punch stretching tests. The results of the finite element method simulations agree well with the values of the fracture heights, in comparison to the experimental data. This verifies the ability and potential of the graphical method in industrial engineering.


Keywords: Forming limit curve (FLC), Hecker’s punch stretching tests, Finite element method (FEM), DP590 steel sheet, Modified maximum force criterion

1. Introduction

Today, sheet metal forming methods based on the deformation of materials play an important role in mechanical production and metallurgy. The growing applications of numerical simulations in the field of sheet metal forming have helped engineers solve various problems in improving the formability and reducing the cost and time of products. Accurate simulation results are necessary for mold and product design. Many factors affect the final simulation results, but the most important input data for the ductile fracture prediction of a product is the forming limit curve (FLC) of the sheet material. Several studies have been carried out to predict and evaluate the FLC by using experimental and theoretical methods.1-3 In addition, this concept has been widely applied in various commercial finite-element software packages for technical studies.

According to the experimental approach, Hecker4 and Nakazima5 tests are popular methods that have been widely used to clarify the levels of FLCs for sheet metals. However, time-consuming and high-cost computing is the main drawback of this testing method. Therefore, considerable effort has been made to obtain FLCs theoretically. Swift6 can be recognized as a pioneering study on predicting FLC. Hill7 then proposed a way to improve the accuracy of FLC prediction by adopting necking point criteria. Stören and Rice8 developed a solution for FLC prediction by applying a force equilibrium between necking and uniform deformed regions. Banabic9 observed and developed a pre-defect in the material and developed a theory of limited deformation based on imperfections of material thickness.

Hora et al.10 upgraded the Swift diffuse necking criteria and set a modified maximum force criterion (MMFC) by effectively examining the instant deformation state changes until the forming force achieved a maximum value. Some new MMFC models proposed to improve the accuracy of FLC prediction1,9 based on theoretical models by solving systems of equations. To simplify the determination of a theoretical FLC based on the theory of MMFC for sheet metal, Tuan and Kim1 discussed a method to predict FLCs based on three special points of the graphical method.

Currently, the finite element (FE) simulation is an indispensable tool to research, evaluate, and discern physical phenomena modeled by various theoretical and experimental equations. To improve the fitting of experimental data with numerical data, Hochholdinger et al.11 upgraded a new FE simulation model by modifying the Hockett–Sherby equation and determining the Zener–Hollomon parameters. To verify the accuracy of FLC prediction compared with corresponding experiments, various previous researchers12-14 imported forming limit diagram (FLD) curves to FEM software in order to predict fractures and compare them with experimental results.

The current study investigates the press formability of DP590 sheet steel. Various hardening equations were adopted to fit the stress-strain curves of DP590 sheet steel. In addition, FLDs were calculated by a graphical method. The predicted FLCs was then imported to an FEM simulation to observe the fracture occurrence of punch stretching with various notch specimens. The fracture heights of specimens measured by FEM prediction were finally compared with corresponding experimental ones to indicate the predicted accuracy of the hardening equations and FLC calculations by the graphical method.


2. Modified Maximum Force Criterion and Graphical Method
2.1. Review of MMFC

To describe the FLC, Hill7 proposed the following Eq. (1):

σ1ε1dε1=σ11+βdε1(1) 

where β=dε2dε1=σ¯/σ2σ¯/σ1=const

However, this equation is only used for the left-hand side of the FLC. Then, Swift6 expressed the right-hand side of the FLC as

σ1ε1dε1=σ1dε1(2) 

However, the above criteria are insufficient when compared with corresponding experimental data. Therefore, Hora10 improved the abovementioned equation by observing diffuse necking based on the strain path ratio (β). Then, a Modified Maximum Fore Criterion (MMFC) could be presented as Eq. (3).

σ1ε1dε1+σ1ββε1dε1σ1dε1(3) 

where β=dε2dε1=σ¯/σ2σ¯/σ1=const

To confirm the accuracy of the FLC, an MMFC model was applied to various sheet metals15,16 with the following assumptions:

α=σ2σ1fα=σ¯σ1gβ=ε¯ε1(4) 

where σ¯ and ε¯ are the equivalent stress and equivalent strain, respectively.

As shown by plasticity flow criterion dεij=dλσ¯/σij and β=dε2dε1=σ¯/σ2σ¯/σ1, it can be proven that β = β(α) and thus g(β(α)) = γ(α). By applying the above explanation to the principle of work equivalence of Eqs. (5)-(7) can be derived.

σ¯Δε¯=σ1 Δε1+σ2Δε2(5) 
gα=1+αβα/fα;(6) 

By solving and inserting all equations into Eq. (3), we can obtain an MMFC model [Eq. (8)].

H'H1gα1-f'αfαββ'α1ε1(7) 

where H=Hε¯ and H' denote the hardening law and the slope of its curve, respectively.

2.2. Application of Graphical Method1 for DP590

As shown in Eq. (7), the calculation can be obtained by the hardening and yield laws. For special strain modes such as plane strain (β = 0), uniaxial tension (β = -1/2), and equibiaxial tension (β = 1), all terms of Eq. (7) are simply calculated.

To verify FLC prediction by a graphical method, DP590 sheet steel was adopted with a thickness of 1.181 mm. Here, tensile tests were performed at a nominal speed of 20 mm/min. The mechanical properties and chemical composition of DP590 steel are listed in Table 1.

Table 1 
Chemical composition and mechanical properties of DP590 material
Mass (%)
C 0.098 S 0.003 Cu 0.025 Ni 0.02
Mn 1.59 Si 0.087 Sn 0.013 Fe Bal
Direction 45° 90°
Young’s modulus (GPa) 202 206 215
Yield stress (MPa) 368 367.67 384.12
Elongation (%) 26.3 27.5 26.8
R value 0.754 1.193 0.995

The experimental stress-strain data of a DP590 steel sheet in the 0 direction are used to fit the hardening parameters (Table 2 and Fig. 1) of the constitutive functions of the Swift,6 Voce,17 and Kim-Tuan18 models as in Eqs. (8)-(10), respectively.

Table 2 
Parameters of Swift, Voce and Kim-Tuan equations of DP590
Swift Voce
C (MPa) ε0 n σY P Q
1027.8 0.0028 0.184 368.51 368.24 20.966
Kim-Tuan
σ0 K(MPa) t h
368 675.8 81.0 0.342


Fig. 1 
Flow stress curves of DP340 with various hardening models

σ¯=Cε0+ε¯n(8) 
σ¯=σY+P1-exp-Qε¯(9) 
σ¯=σ0+K1-exp-tε¯ε¯+0.002h(10) 

To present the anisotropy of DP590 sheet steel owing to the deference R-value (Table 1), the Hill 1948 yield criterion is used as in Eq. (11) and shows good agreement with the experimental data (Fig. 2).


Fig. 2 
Hill anisotropy of DP590 steel sheet prediction for (a) R-value and (b) normalized yield stress

σ¯=Gσ12+Fσ22+Hσ1-σ22(11) 

where G, F, and H are constants calculated from the R-value as follows:

F=12R222+12R332-12R112G=12R112+12R332-12R222H=12R112+12R222-12R332(11-1) 

and R11, R22, and R33 are anisotropic yield stress ratios. For the plane stress case, only four yield stress ratios are needed. Here, we assumed that the rolling direction (RD) involves a user-defined reference yield stress. Then, R11 = 1, and R22 and R33 are defined as 1.077 and 0.999, respectively, as in the following Eq. (11-2):

R22=r90r0+1r0r90+1R33=r90r0+1r0+r90(11-2) 

To review the simple calculation of FLC by a graphical method for DP590 steel, the plane strain stage (β = 0) was first selected and substituted into Eq. (7): H'H1gα at (β = 0). At any strain stage, Eq. (7) can be rewritten as Eq. (12).

H'HA-Bε(12) 

where A = 1/g(α) and B = (f'/f)*(β/β') are constants.

Thus, A and B can be expressed as a function of α (Eqs. (13) and (14)) when a Von Mises yield function is adopted.

Aα=2-α21-α+α2(13) 
Bα=2-α2α-1261-α+α2(14) 

The calculation for special strain modes such as plane strain (P. S.) (β = 0), uniaxial tension (U. T.) (β = -1/2), and equibiaxial tension (B. T.) (β = 1) can be obtained as listed in Table 3 for A and B in Eqs. (13) and (14). Based on the calculated values and the abovementioned equations, Eq. (12) can be plotted as shown in Fig. 3(a) by adopting the Swift, Voce, and Kim-Tuan models.

Table 3 
Calculated constants at plane strain, uniaxial tension, and equibiaxial tension stages
Forming mode (U. T.) (P. S.) (B. T.)
Mises α 0 0.5 1
β -0.5 0 1
A(α) 1 0.866 0.5
B(α) 1/3 0 0.1667
Hill48 α 0.17 0.5 1
β -0.5 0 1
A(α) 1.024 0.885 0.499
B(α) 0.203 0 0.182


Fig. 3 
Evaluation of limit strains based on different hardening laws and yield functions of (a) Von Mises and (b) Hill 48

A(α) and B(α) according to the Hill 48 yield function can be obtained by applying Eq. (11) to Eqs. (4)-(6). The results are also presented in Eqs. (15)-(21), Table 3 and Fig. 3(b).

fα=G+H-2Hα+F+Hα2(15) 
f'α=-H+F+HαG+H-2Hα+F+Hα2(16) 
βα=F+Hα-H-Hα+G+H(17) 
β'α=F+HG+H-H2-Hα+G+H2(18) 
gα=G+H-2Hα+F+Hα2-Hα+G+H(19) 

Then:

Aα=-Hα+G+HG+H-2Hα+F+Hα2(20) 
Bα=-H+F+HαG+H-2Hα+F+Hα2×F+Hα-H×-Hα+G+HF+HG+H-H2(21) 

According to the graphical method, the equivalent plastic strains can be determined. To calculate the major and minor strains, Von Mises (Eq. (22)) and Hill 1948 (Eq. (23)) can be adopted. Then, FLC is obtained by drawing three special strain points as shown in Figs. 4(a) and 4(b).


Fig. 4 
FLC prediction by graphical method for DP590 steel sheet based on various hardening laws in cases of (a) Von Mises and (b) Hill 48

ε1=ε¯/431+β+β2(22) 
ε1=ε¯/gα(23) 

where ε1, ε2 and β = ε2/ε1 are the major and minor strains and the strain ratio, respectively.


3. Experimental and Finite Element Method

In order to verify the accuracy of the FLD predictions, Hecker punch stretching tests4 (Fig. 5(a)) at 15 t of blank holder force and 300 mm/min of punch speed were performed for various specimens (Fig. 5(b)) until the fractures occurred (Fig. 5(c)) in order to measure the fracture heights of the specimens (Fig. 5(d)).


Fig. 5 
Hecker punch stretching (a) for various specimens, (b) experimental results, and (c) fracture heights of DP590 sheet

To analyze the deformation of Hecker’s punch stretching tests, the commercial software package ABAQUS (6.12)19 was used. The FE models are depicted in Fig. 6. Specifically, reduced-integration four-node linear shell elements (S4R) are utilized for the blank. The die/tool and clamping parts are assigned as the rigid body. The size of the element for all parts correspond to approximately 1 × 1 mm (length × width).


Fig. 6 
FEM model for Hecker punch stretching tests

The mesh density is sufficiently fine to ensure result convergence. The mesh size is also refined by comparing the initial simulations and corresponding experiments. The ENCASTRE condition was assigned for the die. The tool and clamming part are permitted to run in the vertical direction. The blank specimen is held by the concentration force acting on the clamping part. The Coulomb friction law is used in the friction condition.


4. Results and Discussion

When the punch stretching tests are performed by using DP590 steel sheets, fractures tend to occur at high tensile stresses and plastic strains owing to the excessive clamping force. Based on the ductile fracture criterion of FLCs, the equivalent plastic strains/stresses change when the specimen shapes are modified. As a follow-up, the values of FLC prediction were altered. In order to verify the effect of the hardening laws and the FLC prediction on the fracture heights of various specimens, FLC input data from the Voce, Swift, and Kim-Tuan equations following the Hill 48 yield function were implemented in order to predict fracture occurrences in various specimens. This is shown in Fig. 7.


Fig. 7 
Deformed shapes of (a) Voce, (b) Swift, and (c) Kim-Tuan’s FLD predictions

In the FE simulation of the forming process, fractures occur when the failure criterion value reaches 1.0, and the element is then deleted. Here, the fracture heights when the fracture elements of deformed shapes occur are verified and compared with the corresponding experiments, as listed in Table 4. To investigate the accuracy of fracture height prediction based on FLC as derived from the Voce, Swift, and Kim-Tuan models, the differences between the fracture heights of the FE method and experimental measurement were estimated as an error of the fracture height Δh(%) = (|hsimulationhexperiment|) / hexperiment.

Table 4 
Comparison of fracture height between experiments and simulation
Fracture height Exp. Simulation Error
Voce Swift Kim-Tuan Voce Swift Kim-Tuan
Specimens No. he (mm) hv (mm) hs (mm) hkt (mm) Δhv (%) Δhs (%) Δhkt (%)
1 36.76 25.75 35.00 36.50 29.95 4.78 0.70
2 37.83 22.50 33.00 35.25 40.52 12.76 6.81
3 31.59 20.25 29.50 31.25 35.89 6.61 1.07
4 25.04 19.75 25.75 26.00 21.12 2.83 3.83
5 23.28 19.00 23.75 23.75 18.38 2.01 2.01
6 24.57 21.00 24.50 24.50 14.52 0.28 0.28
7 21.44 20.00 21.50 21.50 6.71 0.27 0.27
Average error (%) 23.87 4.23 2.15

The results show that the Swift model exhibits reasonably good predictions for almost specimens (Fig. 7(a) and Table 4). Otherwise, the Voce model shows incorrect predictions of fracture heights for DP590 sheet steel. Fig. 7(c) and Table 4 also show the results of finite elements from the simulation of punch stretching tests with Kim-Tuan models. The comparison results exhibit additional good agreement between the simulations and experiments with an average error of 2.15(%). The fracture height error of each specimen of the Kim-Tuan model is more uniform than the others.

To determine the accuracy of the FLCs presented in Table 4, the strain path evolution of each specimen is shown in Fig. 8. Both the Swift and Kim-Tuan models exhibit good predictions on the left-hand side for specimens 5-7, while the Kim-Tuan model shows more accurate predictions on the right-hand side of FLC for the DB590 steel sheet (specimens 1-4).


Fig. 8 
Strain paths by FEM for each specimen


5. Conclusion

This study applied the graphical method to estimate and predict specific FLDs of DP590 steel sheets. Ductile fracture criteria based on the obtained FLDs were implemented by using an FEM simulation and were compared with the corresponding experimental results of punch stretching tests. The prediction of fracture heights using various FLDs derived by difference hardening laws was compared with that of corresponding experiments. The finite element simulation method results indicated that the Kim-Tuan hardening law predicted more accurate fracture heights of various specimens. Meanwhile, the Voce hardening law showed underestimated predictions. The compared results also indicated that the predicted FLCs can be used to predict the fracture position in a simulation with a smallest estimated average error about 2.15% for the Kim-Tuan model and the largest error of 23.87% for the Voce model.


NOMENCLATURE
σ1, σ2 : Two principal stress, σ1 > σ2
ε1, ε2 : Two principal strain, ε1 > ε2
b : Strain path, β=dε2dε1
a : Stress path, a=σ2σ1
σ¯,ε¯ : Equivalent stress and equivalent strain
H(ε¯), H'(ε¯) : Hardening law and its slope curve
C, ε0, n, P, Q, K, t, h : Coefficients in the hardening models
G, F, H, R11, R22, R33 : The anisotropic constants of Hill 1948 yield criterion

Acknowledgments

This paper was presented at PRESM 2019

This research is funded by the Vietnam National Foundation for Science and Technology Development (NAFOSTED) under Grant No. “107.02-2016.01”.


REFERENCES
1. Pham, Q. T., Nguyen, D. T., Kim, J. J., and Kim, Y. S., “A Graphical Method to Estimate Forming Limit Curve of Sheet Metals,” Key Engineering Materials, Vol. 794, pp. 55-62, 2019.
2. Mac, T.-B., Do, V.-C., and Nguyen, D.-T., “A Study of Combined Finite Element Method Simulation/Experiment to Predict Forming Limit Curves of Steel DP350 Sheets,” Advances in Mechanical Engineering, Vol. 10, No. 4, pp. 1-9, 2018.
3. Slota, J. and Spisak, E., “Comparison of the Forming-Limit Diagram (FLD) Models for Drawing Quality (DQ) Steel Sheets,” Metalurgija, Vol. 44, No. 4, pp. 249-253, 2005.
4. Hecker, S. S., “Simple Technique for Determining Forming Limit Curves,” Sheet Metal Industries, Vol. 52, No. 11, pp. 671-676, 1975.
5. Yue, Z., Badreddine, H., Dang, T., Saanouni, K., and Tekkaya, A., “Formability Prediction of Al7020 with Experimental and Numerical Failure Criteria,” Journal of Materials Processing Technology, Vol. 218, pp. 80-88, 2015.
6. Swift, H. W., “Plastic Instability under Plane Stress,” Journal of the Mechanics and Physics of Solids, Vol. 1, No. 1, pp. 1-18, 1952.
7. Hill, R., “On Discontinuous Plastic States, with Special Reference to Localized Necking in Thin Sheets,” Journal of the Mechanics and Physics of Solids, Vol. 1, No. 1, pp. 19-30, 1952.
8. Stören, S. and Rice, J., “Localized Necking in Thin Sheets,” Journal of the Mechanics and Physics of Solids, Vol. 23, No. 6, pp. 421-441, 1975.
9. Banabic, D., Aretz, H., Paraianu, L., and Jurco, P., “Application of Various FLD Modelling Approaches,” Modelling and Simulation in Materials Science and Engineering, Vol. 13, No. 5, pp. 759-769, 2005.
10. Hora, P., Tong, L., and Reissner, J., “A Prediction Method for Ductile Sheet Metal Failure in Fe-Simulation,” Proc. of Numisheet, pp. 252-256, 1996.
11. Hochholdinger, B., Hora, P., Grass, H., and Lipp, A., “Simulation of the Press Hardening Process and Prediction of the Final Mechanical Material Properties,” Proc. of AIP Conference, pp. 618-625, 2011.
12. Nguyen, D.-T., Dinh, D.-K., Nguyen, H.-M. T., Banh, T.-L., and Kim, Y.-S., “Formability Improvement and Blank Shape Definition for Deep Drawing of Cylindrical Cup with Complex Curve Profile from SPCC Sheets Using FEM,” Journal of Central South University, Vol. 21, No. 1, pp. 27-34, 2014.
13. Nguyen, D.-T. and Kim, Y.-S., “A Numerical Study on Establishing the Forming Limit Curve and Indicating the Formability of Complex Shape in Incremental Sheet Forming Process,” International Journal of Precision Engineering and Manufacturing, Vol. 14, No. 12, pp. 2087-2093, 2013.
14. Nguyen, D.-T., Kim, Y.-S., and Jung, D.-W., “Coupled Thermomechanical Finite Element Analysis to Improve Press Formability for Camera Shape Using AZ31B Magnesium Alloy Sheet,” Metals and Materials International, Vol. 18, No. 4, pp. 583-595, 2012.
15. Paraianu, L., Dragos, G., Bichis, I., Comsa, D. S., and Banabic, D., “A New Formulation of the Modified Maximum Force Criterion (MMFC),” International Journal of Material Forming, Vol. 3, No. 1, pp. 243-246, 2010.
16. Pham, Q.-T., Lee, B.-H., Park, K.-C., and Kim, Y.-S., “Influence of the Post-Necking Prediction of Hardening Law on the Theoretical Forming Limit Curve of Aluminium Sheets,” International Journal of Mechanical Sciences, Vol. 140, pp. 521-536, 2018.
17. Voce, E., “The Relationship between Stress and Strain for Homogeneous Deformation,” Journal of the Institute of Metals, Vol. 74, pp. 537-562, 1948.
18. Pham, Q. T. and Kim, Y. S., “Identification of the Plastic Deformation Characteristics of Al5052-O Sheet based on the Non-Associated Flow Rule,” Metals and Materials International, Vol. 23, No. 2, pp. 254-263, 2017.
19. Hibbit, D., Karlsson, B., and Sorensen, P., “ABAQUS User’s Manual,” http://dsk.ippt.pan.pl/docs/abaqus/v6.13/books/usb/default.htm (Accessed 28 AUG 2019)

Luyen The-Thanh

Ph.D. candidate in the Faculty of Mechanical Engineering, Hungyen University of Technology and Education, Vietnam. His research interests are CAD/CAM/CAE and sheet metal forming.

E-mail: luyenthethanh@gmail.com

Pham Quoc-Tuan

Ph.D. candidate in the Department of Mechanical Engineering, National University in Daegu, Korea. His research interests are plasticity and sheet metal forming.

E-mail: tuanphamquoc@knu.ac.kr

Kim Young-Suk

Professor in the Department of Mechanical Engineering, National University in Daegu, Korea. His research interests are sustainable material processing, nano/micro mechanics, FEM and biomechanics.

E-mail: caekim@hnu.ac.kr

Duc-Toan Nguyen

Associate Professor in the School of Mechanical Engineering, Hanoi University of Science and Technology, Vietnam. His research interests are Plasticity, Machining Process, CAD/CAM/CAE.

E-mail: toan.nguyenduc@hust.edu.vn