Next Article in Journal
A Quantitative and Optimization Model for Microstructure Uniformity of Sinter Based on Multiple Regression-NSGA2
Next Article in Special Issue
Ultrasonic-Assisted Granular Medium Forming of Aluminum Alloy 6063-T5: Simulations and Experiments
Previous Article in Journal
Structure and Properties of WC-Fe-Ni-Co Nanopowder Composites for Use in Additive Manufacturing Technologies
Previous Article in Special Issue
The Effect of Powder Temperature on Semi-Solid Powder Rolling AA2024 Based on Experiments and Numerical Simulation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Phenomenological Model to Predict Forming Limit Curves from Tensile Properties for Hot-Rolled Steel Sheets

1
Shi-Changxu Innovation Center for Advanced Materials, Institute of Metal Research, Chinese Academy of Sciences, Shenyang 110016, China
2
School of Materials Science and Engineering, University of Science and Technology of China, Shenyang 110016, China
3
Research Institute, Baoshan Iron and Steel Co., Ltd., Shanghai 201999, China
*
Author to whom correspondence should be addressed.
Metals 2024, 14(2), 168; https://doi.org/10.3390/met14020168
Submission received: 2 December 2023 / Revised: 2 January 2024 / Accepted: 10 January 2024 / Published: 29 January 2024
(This article belongs to the Special Issue Advanced Forming Process of Light Alloy)

Abstract

:
A phenomenological model for the prediction of the forming limit curve (FLC) based on basic mechanical properties through a uniaxial tensile test can tremendously shorten the design time of the forming process and reduce the measuring costs. In this paper, a novel phenomenological model named the IMR-Baosteel model (abbreviated as the IB model) is proposed for efficient and accurate FLC prediction of hot-rolled steel sheets featuring distinct variations in thickness and mechanical properties. With a systematic test of the plane strain forming limit (FLC0), it was found that a higher regression correlation exists between the FLC0 and the total elongation under different sheet thicknesses. For accurate assessment of the FLC0 from tensile properties, compared using experiments, the error of FLC0 calculated with the proposed model is within 10%. In the IB model, the left side of FLC can be calculated using a line with a slope of −1 while the right side of the FLC is obtained via a modified Keeler model with the exponent (p) determined as 0.45 for hot-rolled steels. Complete experimental FLCs of hot-rolled steels from measurements and the literature were used to validate the reliability of the proposed model. Resultantly, the prediction of FLCs with the proposed IB model is greatly improved, and agrees much better with the experimental FLCs than the predictions of the well-known Keeler model, Arcelor model and Tata Steel model.

1. Introduction

Hot-rolled high-strength steels are employed widely in chassis parts of passenger cars and commercial vehicles to reduce automotive weight for energy saving and carbon dioxide emission reduction [1,2,3]. In the automotive industry, excessive thinning or necking are both unacceptable. When sheet thickness is reduced and material strength increases, the formability of the material is always decreased. It is an especially great challenge to produce complex components with high-strength steels, which always need trial and error iterations for the design of the component profile and the forming process [4]. Thus, accurate evaluation of hot-rolled sheet formability is essential for improving design efficiency and the quality of complex components.
The concept of the forming limit curve (FLC) was initially proposed in 1963 [5], and was usually determined from standard experiments using the Marciniak [6] test or the Nakazima test [7]. Nowadays, the FLC is popularly chosen to provide efficient prediction of the failure risk in sheet metal forming processes [8,9,10]. In the FLC, the FLC0 is the forming limit under plane strain conditions, which is near the vertical axis and is usually the lowest point of the FLC. The strain values of the FLC are usually determined through strain measurement procedures on gridded specimens. The main problem is that the experimental determination of the FLC is always costly, time-consuming and inconsistent, and the measurement results highly depend on the mechanical properties of the specimen being tested at that time. Another problem is that friction affects the location, deformation behavior and strain path of the necking point during the experimental tests [11]. In order to improve the quality and efficiency of FLC experiments, many test methods have been developed. A new procedure based on the hydraulic bulging of a double specimen was proposed for the experimental determination of the FLCs to reduce the frictional effects [12]. A test procedure for determining the complete left-hand side of the FLC via tensile tests without friction was outlined [13]. A method was developed to determine the width of specimens to obtain the FLC with the minimum number of tests [14]. However, shortcomings in FLC testing such as laborious measurement and data discretization still cannot be overcome.
Then, much work has been performed to determine the FLC more cheaply, efficiently and accurately during the last few years. One research direction is to develop theoretical models to estimate the FLC. The Hill–Swift model is based on the Swift diffuse instability theory [15] and the Hill localized instability theory [16], and has been expanded to many modified models. However, it has been reported that the Hill–Swift model delivers too-small FLD0 values [17]. The Marciniak–Kuczyński (MK) model is the most widely used model to estimate the necking limit strain for sheet metals [18]. However, it is not popular to estimate the formability of hot-rolled steels with the MK model. This is because the characteristics of hot-rolled steels are weak anisotropy and a low strain hardening exponent, which are different from cold-rolled steel, stainless steel and aluminum alloy [19,20]. Furthermore, it is typical to calculate the forming limit under plane stress conditions with these theoretical models, and the relationship between the through-thickness stress and sheet thickness was ignored. The nonconstant through-thickness normal stress was presented as a critical factor in the FLC prediction for medium plate [21]. For hot-rolled steels with a larger thickness range, the complicated calculation process of the FLC cannot meet the need for efficient formability evaluation in the automotive industry [22].
Another research direction is to develop empirical methods from simple and low-cost experiments to predict the FLC. Empirical methods based upon basic mechanical properties through tensile tests to predict the FLC have been popular for many decades [23]. Keeler and Brazier [24] proposed a regression equation to predict the FLC0 with a strain hardening exponent and thickness. Raghavan [25] described an equation to predict the FLC0 with total elongation and sheet thickness. Paul [26,27] proposed a nonlinear regression equation to predict the FLC0 with ultimate tensile strength, total elongation, strain hardening exponent and sheet thickness. Cayssials [28] developed a predictive method with the strain rate sensitivity, the strain hardening and the sheet thickness based on plastic instability and damage theories. Furthermore, the model was extended by Cayssials [29] to predict the FLC of ultrahigh-strength steels with ultimate tensile strength, uniform elongation, the anisotropy coefficient and the sheet thickness; this is called the Arcelor model. Abspoel [30,31] developed a model to predict the FLC with four strain points including the uniaxial tensile necking point, the plane strain point, the intermediate biaxial stretching point and the equi-biaxial stretching point. The four representative points are calculated using the Lankford coefficient, total elongation and sheet thickness. This model is also called the Tata Steel model. Gerlach [32,33] provided equations to calculate the three characteristic points of the FLC based on three parameters, including ultimate tensile strength, total elongation and sheet thickness. Among these empirical methods, the Keeler model, Arcelor model and Tata Steel model have been integrated into commercial finite element simulation software AutoForm R5.2 [34]. However, since these predictive models are mainly developed from cold-rolling steel sheets or aluminum alloy sheets, the FLC prediction results for hot-rolled steel sheets have a large deviation from the experimental results.
Consequently, in order to shorten the time for the forming process design and reduce the test costs, a reliable phenomenological model named the IMR-Baosteel model (IB model) is established based upon tensile properties in this work, and can effectively predict the formability of hot-rolled steel sheets.

2. Prediction of FLD0 with Classic Empirical Models for Hot-Rolled Steel Sheets

2.1. Data Collection from Experimental Tests

The tensile test is the most widely used method to determine the mechanical properties of materials. Therefore, empirical models are always derived from the statistical relations between the experimental FLC points and the mechanical properties. The details of the test data measured from the Baosteel laboratory are tabulated in Appendix A (Table A1). There are almost eighty hot-rolled steel sheets in various thickness and strength ranges for this investigation obtained from regular steel production in the Baosteel workshop. Figure 1a shows the range of the mechanical properties. The ultimate tensile strength varies between 200 MPa and 1100 MPa, and the total elongation varies between 10% and 50%. The thickness of the sheets varies between 1.5 mm and 6.0 mm, as shown in Figure 1b.
The mechanical properties were obtained from uniaxial tensile tests, performed according to ISO 6892-1:2019 [35], on Instron testing equipment. The collected mechanical properties are yield strength Rp, ultimate tensile strength Rm, total elongation At, uniform elongation Ag, the plastic strain ratio r-value and strain hardening exponent n-value. The r-values and the n-values were determined between 2% and 20% strain or between 2% and Ag when the Ag was lower than 20%, according to ISO 10113:2020 [36] and ISO 10275:2020 [37], respectively. The gauge length L0 to measure At is correlated with sheet thickness, which can be determined as
L 0 = 5.65 t b 0
where b 0 is the specimen width of the measurement area in the tensile test.
The experimental strains of FLC0 and complete FLCs were obtained from Nakazima tests according to ISO 12004-2: 2008 [38] using a 750 KN Interlaken sheet metal testing machine with the Vialux photogrammetric measurement system. A pattern of 2 mm square grids was applied to the surface of the specimens using the electrochemical method. Then the specimens were deformed until fracture using a hemispherical punch with a diameter of 100 mm. Ten specimens with the same length of 196 mm but different widths (varying from 20 mm to 180 mm) were selected to obtain limit strains under different loading paths. Finally, the complete experimental FLCs were obtained with these limit strains under different strain states. The specimens were measured transverse to the rolling direction. Additionally, two specimens of the same size were tested to take the average value. The FLC0 is usually determined using the widths of 90 mm and 100 mm, which are nearest to the plane strain state.
FLC0 is the forming limit for plane strain conditions, and denotes the lowest point of the FLC. So the accurate determination of the FLC0 is primarily important for predicting the FLC. Experimental FLC0 values were plotted with mechanical properties, and the influence of these mechanical properties on the characteristics of the FLC0 was studied.
Figure 2 shows the correlation between the FLC0 and tensile properties. Rp (Figure 2a), Rm (Figure 2b), Ag (Figure 2c), At (Figure 2d) and the n-value (Figure 2e) show an approximately linear trend with the FLC0, while there is no significant correlation between the FLC0 and the r-value (Figure 2f) or thickness t (Figure 2g). Furthermore, among these mechanical properties obtained from regular tensile tests, the coefficients of determination of three properties including Ag, At and the n-value are more than 0.8, which shows a stronger correlation with the FLC0.

2.2. Prediction Results for FLC0 with Classic Empirical Models

The Keeler model is the most popular method for predicting the FLC, especially in the automotive industry [39]. However, the comparison showed that the Keeler model was only reliable for classic forming-grade steels [40]. In this section, these classic empirical models including the Keeler model, Raghavan model, Paul model, Tata Steel model and Arcelor model are employed to verify the prediction reliability for hot-rolled steel sheets. The collected FLC0 data of hot-rolled steel sheets are grouped according to sheet thickness and tensile strength.
The Keeler model is shown as
F L C 0 = l n 1 + 23.3 + 14.13 t 21 n
The Raghavan model is shown as
F L C 0 = 2.78 + 3.24 t + 0.892 A t
The Paul model is expressed as
F L C 0 = 7.702 e x p 0.0122 R m 0.1124 r 0.6908 e x p 12.4187 A t + 0.1149 n + 0.0823 t + 0.3011
The Tata Steel model is expressed as
F L C 0 = 0.0084 A t + 0.0017 A t t 1
Furthermore, the equation for the Arcelor model was not provided in papers and the FLC prediction can be obtained from AutoForm R7 [41].
First, the calculation results from the well-known empirical models were compared with the experimental FLC0 according to the stratification of sheet thicknesses. For sheets of a thickness less than 3 mm, as shown in Figure 3a, the Keeler model and Arcelor model can predict well for sheets with high formability, and the predicted deviation is even lower than 10%. However, the Keeler model and Arcelor model severely underestimate the FLC0 of sheets with low formability. The Paul model and Raghavan model slightly overestimate the FLC0, with the predicted deviation between 10% and 30%, while the Tata Steel model underestimates the formability, with the deviation exceeding 10%. For sheets of a thickness greater than 3 mm, as shown in Figure 3b, the prediction results with the above models have a large scatter. The predicted deviations of the Keeler model, Raghavan model and Tata Steel model are barely less than 30%.
Second, the calculation results using the well-known empirical models were compared with the experimental FLC0 according to the stratification of tensile strength. For the sheets with a tensile strength lower than 550 MPa, as shown in Figure 3c, the Paul model and Raghavan model slightly overestimate the FLC0 while the Tata Steel model slightly underestimates the FLC0 with the predicted deviation almost between 10% and 30%. In comparison, the Keeler model has the best prediction accuracy with a deviation less than 10%. For the sheets with a tensile strength higher than 550 MPa, as shown in Figure 3d, the prediction accuracy of these empirical models is much more unreliable. The predicted deviation of the Raghavan model is just near 30%, while the predicted deviation of the other models actually even exceeds 30%.
In summary, for hot-rolled steel sheets, when the sheet thickness is less than 3.0 mm and tensile strength is lower than 550 MPa, the prediction accuracy of the Keeler model is comparatively reliable. However, when the sheet thickness is greater than 3.0 mm or the tensile strength is greater than 550 MPa, the prediction accuracy of the above five empirical models is significantly poor.

2.3. Critical Mechanical Properties for FLD0 Prediction

In a standard uniaxial tensile test, the digital image correlation (DIC) method is used to measure strain and elongation. Figure 4a shows the engineering stress–strain curve of hot-rolled high-strength steel S550MC with a thickness of 2.5 mm. The engineering stress–strain curve is transformed into the real stress–strain curve by fitting the index; then, the n-value can be obtained. The n-value is often used to describe sheet formability as an important parameter, such as in the Keeler model. However, as shown in Figure 4b, there is a yield plateau on the engineering stress–strain curve of hot-rolled steel. Due to that, it is not sufficient to describe the stress–strain behavior with the power law equation [42], which further means the n-value obtained from the power law equation cannot accurately capture the actual strain hardening behavior.
Figure 4c shows the true plastic strain measurement of the local necking point in the uniaxial tensile test. It can be seen that the plastic strain in the whole narrow region of the test specimen is uniform before the time of the maximum uniaxial tensile force. During this period, the strain ratio of the true plastic width strain to the true plastic longitudinal strain is stable at −1/2. It is worth mentioning that, at the time of maximum uniaxial tensile force, the uniform elongation is approximately equal to the true plastic longitudinal strain. After the time of maximum uniaxial tensile force, the region where the plastic strain increment continues reduces gradually until local instability and fracture occur. Simultaneously, the strain state changes from uniaxial tension to a plane strain condition. It is clear that the distance between the uniform elongation point and local onset necking point is long for hot-rolled steel sheets. The local onset necking point is much closer to the specimen fracture point which corresponds to the total elongation. Therefore, a stronger correlation between At and the FLC0 rather than Ag is verified.
Furthermore, the effect of sheet thickness on the FLC0 has been widely reported in the literature [43,44]. And the influence of thickness on the FLC is significant especially for hot-rolled steel sheets; the explanation for this is that as thickness increases, local necking becomes more diffuse and the time to reach the critical depth of fracture which is defined as failure is increased [45].

3. Establishment of New Prediction Model for FLC0

According to the above analysis, the key parameters of At and thickness t were adopted to establish the prediction model of the FLC0 in this work. The experimental data of FLC0 with different thicknesses of hot-rolled steel sheets were extracted to investigate the mathematical relation between the FLC0 and the total elongation At, along with the thickness.
As illustrated in Figure 5, the correlation of At with the FLC0 was studied from the typical thicknesses 2.0 mm, 2.5 mm and 3.0 mm. It is obvious that the correlation of the FLC0 with At is not linear at different thicknesses. Then, a cubic polynomial equation was used to regress the correlation of FLC0 with At:
F L C 0 = A 0 + A 1 A t + A 2 A t 2 + A 3 A t 3
A 0 A 1 A 2 A 3 are the constant parameters of the cubic polynomial equation. Then the parameters were fitted from the data in Figure 5a–c, and are summarized in Table 1.
Furthermore, the influence of thickness on the FLC0 was studied based on the FLC0 of 2.0 mm thickness. The correlation ratio between strain and thickness (CRST) is defined as the ratio of the FLC0 of other thicknesses to the FLC0 of 2.0 mm. The individual experimental FLC0 of QStE600 TM at different thicknesses, including 2.0 mm, 2.5 mm, 3.5 mm and 5.0 mm, was employed to determine the CRST. The calculating results for the CRST are shown in Table 2.
Then the CRST data were plotted vs. the thickness, as shown in Figure 6. The CRST can be calculated as
C R S T = 1.05 t 1.31 0.142
Combining Equations (6) and (7), the predictive model of the FLC0 with total elongation and thickness is established as
F L C 0 = 0.491 3.88 A t + 16.11 A t 2 17.20 A t 3 1.05 t 1.31 0.142
Figure 7 shows the prediction capability of the proposed model compared with the Keeler model. It is clear that the predicted deviations calculated with the Keeler model mostly exceed 10% and even exceed 30% for high-strength steel sheets. By contrast, the predicted deviations calculated with the proposed model are almost under 10%. The results show that the proposed model can accurately predict the FLC0 for hot-rolled steel sheets and the performance is much better than that of the Keeler model in the area of high-strength steel sheets especially.

4. Determination of Phenomenological Model for Complete FLC

A complete FLC consists of two limit curves located in the tension–tension and tension–compression domains, respectively. The FLC covers almost the entire deformation domain in the sheet metal forming processes. In general, the strain ratio spans between those induced by uniaxial and equi-biaxial loads.
Levy [46] described that the slope of the uniaxial tensile strain path in the forming limit diagram depends upon the r-value. The higher the r-value, the greater the slope of the uniaxial tensile strain path in the forming limit diagram. However, hot-rolled steels have the characteristic of weak anisotropy. The distribution of r-values for hot-rolled high-strength steel is concentrated between 0.7 and 0.9 (Figure 2f). As a result, the shapes of FLC curves for hot-rolled steel sheets are almost similar and the main difference is the height of the curves.
According to [24], the left side of the strain-based FLC can be calculated with an equation with a slope of −1:
ε 1 = F L C 0 ε 2
where ε 1  is major strain and ε 2  is minor strain.
Then the right side of the strain-based FLC can be calculated with [26]:
ε 1 = 1 + F L C 0 1 + ε 2 p 1
where p is a material constant.
To determine the parameter of p, the experimental FLCs of SAPH440, QStE460TM, QStE600TM and S700MC were employed for further analysis.
FLC0 can be calculated with Equation (8), the left side of the FLC can be calculated with Equation (9) and the right side of the FLC can be calculated with Equation (10) with different values of p, as shown in Figure 8. The right side of the predicted FLC agrees well with the experimental FLC. Therefore, p was determined as 0.45 for hot-rolled steels. Furthermore, the left side of the FLC predicted with Equation (9) also agrees well with the experimental FLC.
Therefore, the new phenomenological model for the complete FLC of hot-rolled steel sheets can be determined with the combination of Equations (8)–(10). This phenomenological model is named the IMR-Baosteel model, which is shortened to the IB model. Then several experimental FLCs were collected to validate the reliability of the IB model compared with well-known models including the Keeler model, Arcelor model and Tata Steel model.
On one hand, the experimental FLCs of SAPH440, S550MC, S700MC and FB780 tested in the laboratory were employed to verify the reliability of the IB model. As illustrated in Figure 9a, the experimental FLC of low-strength steel SAPH440 can be predicted well with both the Keeler model and the proposed IB model. The left-hand side of the Tata Steel model and the right-hand side of the Arcelor model agree with the experimental points. In contrast, the slopes of the right-hand side of the Tata Steel model and the left-hand side of the Arcelor model deviate from the experimental points. Then the prediction of FLCs with the proposed IB model can agree with the experimental FLCs for hot-rolled high-strength steels (Figure 9b–d). However, these three classic empirical models cannot accurately predict the FLC, and the main deviation derives from the underestimated prediction of the FLC0, especially for high-strength steels.
On the other hand, the experimental FLCs of the SAPH370, QStE340TM, QStE550TM, 580DP, 700DP and Q-P-T steels were collected from the literature to verify the reliability of the IB model. As shown in Figure 10a–c, for SAPH370, QStE340TM and QStE550TM steel, it is clear that the Keeler model slightly underestimates the height of the FLC, while it overestimates the height of the FLC for 580DP and 700DP steel (Figure 10d,e). The Arcelor model and Tata Steel model both underestimate all of these FLCs except SAPH370. Generally, the curve slope of the left-hand side predicted with the Tata Steel model overestimates the measuring points while the curve slope of the right-hand side predicted with the Tata Steel model underestimates the measuring points. Obviously, the prediction with the IB model agrees better with the experimental FLCs (Figure 10a–e). For quenching–partitioning–tempering (Q-P-T) steel (Figure 10f), the IB model can predict the FLC0 well. However, the IB model overestimates the left side of the FLC and underestimates the right side of the FLC, perhaps due to the low r-value, about 0.27, of Q-P-T steel.

5. Conclusions

In this work, a new phenomenological model named the IB model based on tensile properties is proposed to predict the FLC for hot-rolled steel sheets accurately and efficiently. The main conclusions can be summarized as follows:
(1)
The effect of tensile properties on the plane strain forming limit (FLC0) was studied with experimental results of eighty hot-rolled steel sheets under various thicknesses and strengths. Classic empirical models were employed to verify the prediction reliability. The results show that when the sheet thickness is less than 3.0 mm and the tensile strength is lower than 550 MPa, the Keeler model has the best prediction accuracy, with a deviation of less than 10%, which is better than other empirical models. However, there are distinct deviations in predicting hot-rolled high-strength steel sheets with all of the current empirical models. For high-strength hot-rolled steel sheets, the Keeler model almost underestimates the FLC, due to the fact that hot-rolled steels have the characteristics of a low strain hardening exponent and higher thickness.
(2)
For hot-rolled steels, there is a yield plateau on the engineering stress–strain curve. Due to this, it is not sufficient to describe the stress–strain behavior with the power law equation, which means the n-value obtained from power law equation fitting cannot describe the hardening behavior accurately. Combined with the correlation analysis and DIC measurement during the tensile test, it was found that there is a stronger regression relationship between the total elongation and the FLC0. Then the IB model, combining the cubic polynomial equation and power equation, was proposed to regress the correlation of the FLC0 with total elongation and thickness. The errors calculated for the FLC0 with the proposed model are mainly under 10% compared with the errors calculated with the Keeler model, which exceed 30–50% for hot-rolled high-strength steels. Additionally, the IB model is applicable for thicknesses between 1.5 mm and 6.0 mm, which covers most hot-rolled steels being employed. And its reliability for hot-rolled steels out of this thickness range is not verified with effective experimental data.
(3)
In the IB model, the left side of the FLC can be calculated via a line with a slope of –1 for the majority of hot-rolled steels with r-values between 0.7 and 0.9, while the right side of the FLC can be obtained via a modified Keeler model with the exponent (p) determined as 0.45 for hot-rolled steels. Ten complete experimental FLCs of hot-rolled steels from measurements and the literature were used to validate the prediction reliability. The results show that the prediction of the complete FLC with the IB model matches much better with the experimental FLC than those with the other empirical models. However, for Q-P-T steel, the IB model can predict the FLC0 well but cannot predict the left and right sides of the FLC accurately, due to the low r-value of about 0.27.

Author Contributions

Conceptualization, H.-W.S., X.-H.P. and S.-H.Z.; methodology, W.-J.C., H.-W.S. and S.-Y.D.; investigation, W.-J.C. and Z.C.; writing—original draft preparation, W.-J.C. and S.-Y.D.; writing—review and editing, H.-W.S., S.-F.C., X.-H.P., Z.C. and Y.X.; supervision, S.-H.Z.; funding acquisition, S.-F.C. and Y.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China (grant number 2022YFB3304800), the National Natural Science Foundation of China (grant number 52111530293 and grant number 52105412), the project of the Institute of Metal Research, Chinese Academy of Sciences (grant number E055A501) and the scientific and technological cooperation project between Jilin Province and the Chinese Academy of Sciences (grant number 2022SYHZ0007).

Data Availability Statement

All data are contained within the article and in the Refs. [20,21,47].

Acknowledgments

The authors would like to thank Qin-Bo Shen and Hua Zhang from Baosteel Research Institute for their help with the experimental phase of this work.

Conflicts of Interest

Authors Zheng Cai and Xin-Hua Pei were employed by the company Research Institute, Baoshan Iron and Steel Co., Ltd. The remaining authors declare that the re-search was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Appendix A

Table A1. Uniaxial tensile properties and experimental FLC0 for data analysis.
Table A1. Uniaxial tensile properties and experimental FLC0 for data analysis.
Steelt (mm)Rp (MPa)Rm (MPa)n-Valuer-ValueAgAtExp FLC0Cal FLC0Keeler FLC0
SPHC22283330.19510.240.440.4250.4360.389
SPHC2.52223410.191.120.2420.4420.450.4720.424
SPHC32453470.1931.150.2420.4450.480.4960.471
SAPH37023264130.1581.060.170.350.40.3680.326
SAPH3702.53234100.1570.950.1710.3550.410.4030.366
SAPH37033094030.1610.920.1730.3640.430.4340.406
SAPH40023414430.1530.910.1710.3450.360.3620.317
SAPH40033354400.1540.830.1740.3470.40.4140.392
SAPH4402.33524770.1480.80.1680.3210.340.3530.330
SAPH4402.53664850.1470.790.1690.3220.350.3640.342
SAPH4403.53434650.1490.830.1710.3280.380.4040.415
SPFH54024675650.1230.780.1240.2730.270.2810.263
SPFH54034935880.1210.80.1260.2760.320.3230.320
SPFH59025156150.10.810.090.2450.260.2530.218
SPFH5902.55326220.1020.860.0940.2460.280.2750.249
SPFH5904.55426460.0950.760.0990.2560.310.3280.331
QStE340TM2.53795160.1320.850.1630.3250.350.3670.312
QStE340TM33835010.1300.870.1610.330.370.3920.347
QStE380TM24015050.1260.810.1430.3010.30.3120.268
QStE380TM2.54095100.1260.80.1410.3080.330.3460.300
QStE380TM3.753785020.1270.850.1470.3120.40.3890.378
S355MC2.53804950.1290.770.1560.3220.340.3640.306
S355MC33925150.1320.870.1570.3230.370.3830.344
S355MC43695030.1360.810.1590.3310.390.4210.415
S355MC63784890.1430.860.1680.3350.410.4610.550
S420MC24675780.1220.810.1250.270.270.2780.261
S420MC34695870.1290.820.1270.2790.320.3270.338
S420MC3.54835920.1170.730.1290.2820.340.3430.339
S420MC54776060.1160.760.1270.3010.370.3960.417
QStE460TM24976230.1030.780.1150.2530.260.2610.224
QStE460TM2.54836030.1070.820.1160.2540.280.2830.260
QStE460TM3.64916110.110.740.1120.2630.30.3210.327
QStE500TM1.85446190.0970.820.1010.2270.220.2270.202
QStE500TM2.55656360.0980.830.1090.2370.270.2660.241
QStE500TM35516590.0960.780.1080.2320.290.2750.262
QStE500TM4.55426430.1020.840.1120.2480.310.3180.351
QStE500TM65536320.1030.810.1160.2520.320.3410.424
QStE550TM25746650.090.80.0980.2160.230.2300.199
QStE550TM2.55956870.0910.760.1230.2110.260.2450.225
QStE550TM2.86046900.0850.820.1120.2190.270.2590.226
QStE550TM35876820.0820.860.1040.2210.2750.2650.227
QStE550TM3.55916640.0890.790.1090.2310.290.2840.268
QStE600TM26337320.0790.780.0940.2070.220.2250.176
QStE600TM2.56357380.080.860.0960.2080.240.2430.201
QStE600TM3.56227160.0810.850.0950.2120.260.2680.246
QStE600TM56277270.0720.810.1030.210.280.2870.265
QStE650TM26747900.0680.770.0910.190.220.2170.154
QStE650TM2.56657820.0660.820.0920.1960.240.2370.173
QStE650TM36617760.0700.830.0950.2050.250.2540.197
QStE700TM1.57398040.0620.780.0820.1810.190.1780.123
QStE700TM1.87257900.0600.830.080.1820.1950.2040.130
QStE700TM27378020.0630.760.0840.1710.20.2120.143
QStE700TM2.57478200.0580.80.0810.1780.230.2300.149
QStE700TM37247960.0610.810.0860.1840.250.2440.174
QStE700TM47327840.0620.730.0870.2040.260.2710.211
BR440/580HE35145740.1680.750.1310.2770.350.3240.421
580DP3.53896360.1750.80.1750.310.40.3800.472
700DP2.54257580.1380.760.1250.240.270.2690.325
780DP3.25778480.1250.820.1210.210.280.2610.341
FB5902.25126210.1020.770.1280.2570.30.2750.233
FB78036658130.0710.810.0840.1870.260.2450.200
FB78046728100.0750.750.0790.1940.290.2650.250
B780NP37568010.0820.880.0850.250.30.2930.227
B780NP3.57668120.0820.830.0910.260.330.3160.249
B780SF2.57848560.0800.890.0810.2450.280.2740.201
B510L34355620.1320.90.1520.2770.330.3240.344
B510L44135350.1350.930.1550.2850.380.3570.413
B510L54205440.1360.890.1580.2940.410.3860.474
B510L64465330.1430.690.1510.310.410.4240.550
B530L34265620.1320.860.1480.290.320.3410.344
B550L54805750.1250.810.1450.2780.350.3640.443
B610L35566370.1040.820.1050.2320.280.2750.271
B610L4.55726520.1030.780.1210.240.310.3090.354
B650L35976880.0930.830.1050.220.260.2650.254
B700L36547320.0870.780.1060.2120.250.2590.252
B750L2.57337850.0620.770.0840.1830.230.2310.159
B750L3.57357930.0640.810.0860.1840.240.2530.200
BWP7501.57207950.0630.750.0840.210.220.1890.125
BWP7503.57188040.0680.840.0960.2050.260.2630.211
BWP75046958050.0710.780.0920.210.280.2750.238
B980285910470.060.750.0670.140.230.2160.137

References

  1. Jha, G.; Das, S.; Lodh, A.; Haldar, A. Development of hot rolled steel sheet with 600MPa UTS for automotive wheel application. Mater. Sci. Eng. A 2012, 552, 457–463. [Google Scholar] [CrossRef]
  2. Hu, J.; Du, L.X.; Wang, J.J.; Sun, Q.Y. Cooling process and mechanical properties design of hot-rolled low carbon high strength microalloyed steel for automotive wheel usage. Mater. Des. 2014, 53, 332–337. [Google Scholar] [CrossRef]
  3. Chen, W.J.; Song, H.W.; Lazarescu, L.; Xu, Y.; Zhang, S.H.; Banabic, D. Formability analysis of hot-rolled dual-phase steel during the multistage stamping process of wheel disc. Int. J. Adv. Manuf. Technol. 2020, 110, 1563–1573. [Google Scholar] [CrossRef]
  4. Sansot, P.; Frédéric, B.; Vitoon, U.; Surasak, S.; Suwat, J. Experimental and theoretical formability analysis using strain and stress based forming limit diagram for advanced high strength steels. Mater. Des. 2013, 51, 756–766. [Google Scholar]
  5. Keeler, S.P.; Backhofen, W.A. Plastic instability and fracture in sheet stretched over rigid punches. ASM Trans. 1963, 56, 25–48. [Google Scholar]
  6. Marciniak, Z.; Kuczynski, K. Limit strains in the processes of stretch forming sheet metal. Int. J. Mech. Sci. 1967, 9, 609–612. [Google Scholar] [CrossRef]
  7. Nakazima, K.; Kikuma, T.; Hasaku, K. Study on the formability of steel sheet. Yawata Technol. Rep. 1968, 264, 8517–8530. [Google Scholar]
  8. Xia, L.L.; Xu, Y.; El-Aty, A.A.; Zhang, S.H.; Nielsen, K.B.; Li, J.M. Deformation characteristics in hydro-mechanical forming process of thin-walled hollow component with large deformation: Experimentation and finite element modeling. Int. J. Adv. Manuf. Technol. 2019, 104, 4705–4714. [Google Scholar] [CrossRef]
  9. Banabic, D.; Aretz, H.; Paraianu, L.; Jurco, P. Application of various FLD modelling approaches. Mod. Sim. Mater. Sci. Eng. 2005, 13, 759–769. [Google Scholar] [CrossRef]
  10. Ma, Y.; Chen, S.F.; Chen, D.Y.; Banabic, D.; Song, H.W.; Xu, Y.; Zhang, S.H.; Fan, X.S.; Wang, Q. Determination of the forming limit of impact hydroforming by frictionless full zone hydraulic forming test. Int. J. Mater. Form. 2021, 14, 1221–1232. [Google Scholar] [CrossRef]
  11. Kasaei, M.M.; Oliveira, M.C. Influence of the contact with friction on the deformation behavior of advanced high strength steels in the Nakajima test. J. Strain Anal. Eng. 2022, 57, 193–207. [Google Scholar] [CrossRef]
  12. Banabic, D.; Lazarescu, L.; Paraianu, L.; Ciobanu, I.; Nicodim, I.; Comsa, D.S. Development of a new procedure for the experimental determination of the Forming Limit Curves. Ann. CIRP 2013, 62, 255–258. [Google Scholar] [CrossRef]
  13. Yang, Q.B.; Min, J.Y.; Carsley, J.E.; Wen, Y.Y.; Kuhlenkötter, B.; Stoughton, T.B.; Lin, J.P. Prediction of plane-strain specimen geometry to efficiently obtain a forming limit diagram by Marciniak test. J. Iron Steel Res. Int. 2018, 25, 539–545. [Google Scholar] [CrossRef]
  14. Holmberg, S.; Enquist, B.; Thilderkvist, P. Evaluation of sheet metal formability by tensile tests. J. Mater. Process Technol. 2004, 145, 72–83. [Google Scholar] [CrossRef]
  15. Swift, H.W. Plastic instability under plane stress. J. Mech. Phys. Solids 1952, 1, 1–18. [Google Scholar] [CrossRef]
  16. Hill, R. On discontinuous plastic states with special reference to localized necking in thin sheets. J. Mech. Phys. Solids. 1952, 1, 19–30. [Google Scholar] [CrossRef]
  17. Bleck, W.; Deng, Z.; Papamantellos, K.; Gusek, C.O. A comparative study of the forming limit diagram models for sheet steels. J. Mater. Process Technol. 1998, 83, 223–230. [Google Scholar] [CrossRef]
  18. Banabic, D.; Kami, A.; Comsa, D.S.; Eyckens, P. Developments of the Marciniak-Kuczynski model for sheet metal formability: A review. J. Mater. Process Technol. 2021, 287, 116446. [Google Scholar] [CrossRef]
  19. Hao, Q.G.; Qin, S.W.; Liu, Y.; Zuo, X.W.; Chen, N.L.; Rong, Y.H. Relation between microstructure and formability of quenching-partitioning-tempering martensitic steel. Mater. Sci. Eng. A 2016, 671, 135–146. [Google Scholar] [CrossRef]
  20. Chen, W.J.; Yin, S.; Pei, X.H. Mechanical property and formability of 580DP and 700DP hot-rolled dual phase steel. Mater. Mech. Eng. 2020, 44, 92–97. [Google Scholar]
  21. Ma, B.L.; Wan, M.; Zhang, H.; Gong, X.L.; Wu, X.D. Evaluation of the forming limit curve of medium steel plate based on non-constant through-thickness normal stress. J. Manuf. Process. 2018, 33, 175–183. [Google Scholar] [CrossRef]
  22. Kim, W.; Koh, Y.; Kim, H. Formability evaluation for hot-rolled HB780 steel sheet based on 3-D non-quadratic yield function. Met. Mater. Int. 2017, 23, 519–531. [Google Scholar] [CrossRef]
  23. Sing, W.M.; Rao, K.P. Prediction of sheet-metal formability using tensile-test results. J. Mater. Process Technol. 1993, 37, 37–51. [Google Scholar] [CrossRef]
  24. Keeler, S.P.; Brazier, S.G. Relationship between laboratory material characterization and press-shop formability. Proc. Microalloying 1977, 75, 517–530. [Google Scholar]
  25. Raghavan, K.S.; Van Kuren, R.C.; Darlington, H. Recent progress in the development of forming limit curves for automotive sheet steel. SAE Technol. Pap. 1992, 920437. [Google Scholar] [CrossRef]
  26. Paul, S.K. Prediction of complete forming limit diagram from tensile properties of various steel sheets by a nonlinear regression based approach. J. Manuf. Process. 2016, 23, 192–200. [Google Scholar] [CrossRef]
  27. Paul, S.K.; Manikandan, G.; Verma, R.K. Prediction of entire forming limit diagram from simple tensile material properties. J. Strain Anal. Eng. 2013, 48, 386–394. [Google Scholar] [CrossRef]
  28. Cayssials, F. A new method for predicting FLC. In Proceedings of the 20th IDDRG Congress, Brussels, Belgium, 17–19 June 1998; pp. 443–454. [Google Scholar]
  29. Cayssials, F.; Lemoine, X. Predictive model for FLC (arcelor model) upgraded to UHSS steels. In Proceedings of the 24th IDDRG Conference, Besançon, France, 20–22 June 2005; pp. 17.1–17.8. [Google Scholar]
  30. Abspoel, M.; Scholting, M.E.; Droog, J.M.M. A new method for predicting forming limit curves from mechanical properties. J. Mater. Process Technol. 2013, 213, 759–769. [Google Scholar] [CrossRef]
  31. Abspoel, M.; Scholting, M.E.; Lansbergen, M.; An, Y.; Vegter, H. A new method for predicting advanced yield criteria input parameters from mechanical properties. J. Mater. Process Technol. 2017, 248, 161–177. [Google Scholar] [CrossRef]
  32. Gerlach, J.; Kessler, L.; Kohler, A. The forming limit curve as a measure of formability is an increase of testing necessary for robustness simulations. In Proceedings of the IDDRG 50th Anniversary Conference, Graz, Austria, 31 May–2 June 2010; pp. 479–488. [Google Scholar]
  33. Gerlach, J.; Kessler, L.; Kohler, A.; Paul, U. Method for the approximate calculation of forming limit curves using tensile test results. Stahl. Und Eisen 2010, 130, 55–61. [Google Scholar]
  34. Pimentel, A.M.F.; de Carvalho Martins Alves, J.L.; de Seabra Merendeiro, N.M.; Vieira, D.M.F. Comprehensive benchmark study of commercial sheet metal forming simulation softwares used in the automotive industry. Int. J. Mater. Form. 2018, 11, 879–899. [Google Scholar] [CrossRef]
  35. ISO 6892-1:2019; Metallic Materials—Tensile Testing—Part 1: Method of Test at Room Temperature. ISO: Geneva, Switzerland, 2019.
  36. ISO 10113:2020; Metallic Materials—Sheet and Strip—Determination of Plastic Strain Ratio. ISO: Geneva, Switzerland, 2020.
  37. ISO 10275:2020; Metallic Materials—Sheet and Strip—Determination of Tensile Strain Hardening Exponent. ISO: Geneva, Switzerland, 2020.
  38. ISO 12004-2: 2008; Metallic Materials—Sheet and Strip—Determination of Forming-Limit Curves—Part 2: Determination of Forming-Limit Curves in Laboratory. ISO: Geneva, Switzerland, 2008.
  39. Shi, M.F.; Gelisse, S. Issues on the AHSS forming limit determination. In Proceedings of the 25th IDDRG Conference, Porto, Portugal, 19–21 June 2006; pp. 19–25. [Google Scholar]
  40. Ghazanfari, A.; Assempour, A. Calibration of forming limit diagrams using a modified Marciniak–Kuczynski model and an empirical law. Mater. Des. 2012, 34, 185–191. [Google Scholar] [CrossRef]
  41. Chezan, A.R.; Khandeparkar, T.V.; ten Horn, C.H.L.J.; Sigvant, M. Accurate sheet metal forming modeling for cost effective automotive part production. In Proceedings of the 38th IDDRG Conference, Enschede, The Netherlands, 3–7 June 2019; Volume 651, p. 012007. [Google Scholar]
  42. Yun, X.; Gardner, L. Stress-strain curves for hot-rolled steels. J. Constr. Steel Res. 2017, 133, 36–46. [Google Scholar] [CrossRef]
  43. Hosford, W.F.; Duncan, J.L. Sheet metal forming: A review. JOM 1999, 51, 39–44. [Google Scholar] [CrossRef]
  44. Hashemi, R.; Ghazanfari, A.; Abrinia, K.; Assempour, A. The effect of the imposed boundary rate on the formability of strain rate sensitive sheets using the M-K method. J. Mater. Eng. Perform. 2013, 22, 2522–2527. [Google Scholar] [CrossRef]
  45. Keeler, S.P. Forming limit criteria sheets. In Advances in Deformation Processing; Burke, J.J., Weiss, V., Eds.; Plenum Press: New York, NY, USA, 1989; pp. 127–157. [Google Scholar]
  46. Levy, B.S.; Van Tyne, C.J. An approach to predicting the forming limit stress components from mechanical properties. J. Mater. Process Technol. 2016, 229, 758–768. [Google Scholar] [CrossRef]
  47. Lu, W.Q.; Chen, J.S.; Chen, J.; Shi, L.; Xiao, H.; Tong, G. Application of minimum thickness criterion in forecasting FLC of high-strength hot rolled-sheet metal. Mater. Sci. Technol. 2010, 18, 387–395. [Google Scholar]
Figure 1. Mechanical properties of eighty measured hot-rolled steels: (a) total elongation with ultimate tensile strength, and (b) ultimate tensile strength with sheet thickness.
Figure 1. Mechanical properties of eighty measured hot-rolled steels: (a) total elongation with ultimate tensile strength, and (b) ultimate tensile strength with sheet thickness.
Metals 14 00168 g001
Figure 2. Correlation between FLC0 and tensile properties: (a) yield strength (Rp), (b) ultimate tensile strength (Rm), (c) uniform elongation (Ag), (d) total elongation (At), (e) strain hardening exponent (n-value), (f) plastic strain ratio (r-value) and (g) sheet thickness (t).
Figure 2. Correlation between FLC0 and tensile properties: (a) yield strength (Rp), (b) ultimate tensile strength (Rm), (c) uniform elongation (Ag), (d) total elongation (At), (e) strain hardening exponent (n-value), (f) plastic strain ratio (r-value) and (g) sheet thickness (t).
Metals 14 00168 g002aMetals 14 00168 g002b
Figure 3. Prediction of plane strain forming limit (FLC0) of hot-rolled steel sheets by different empirical models under the following conditions: (a) sheet thickness less than 3 mm, (b) sheet thickness greater than 3 mm, (c) sheet tensile strength lower than 550 MPa and (d) sheet tensile strength higher than 550 MPa.
Figure 3. Prediction of plane strain forming limit (FLC0) of hot-rolled steel sheets by different empirical models under the following conditions: (a) sheet thickness less than 3 mm, (b) sheet thickness greater than 3 mm, (c) sheet tensile strength lower than 550 MPa and (d) sheet tensile strength higher than 550 MPa.
Metals 14 00168 g003
Figure 4. Measurement during uniaxial tensile test for S550MC: (a) engineering stress–strain curve, (b) true stress–strain curve and fitting and (c) strain path tracking at local necking point using DIC.
Figure 4. Measurement during uniaxial tensile test for S550MC: (a) engineering stress–strain curve, (b) true stress–strain curve and fitting and (c) strain path tracking at local necking point using DIC.
Metals 14 00168 g004
Figure 5. Regression of FLC0 with total elongation: (a) sheet thickness 2.0 mm, (b) sheet thickness 2.5 mm and (c) sheet thickness 3.0 mm.
Figure 5. Regression of FLC0 with total elongation: (a) sheet thickness 2.0 mm, (b) sheet thickness 2.5 mm and (c) sheet thickness 3.0 mm.
Metals 14 00168 g005
Figure 6. Correlation ratio with strain and thickness (CRST) vs. thickness.
Figure 6. Correlation ratio with strain and thickness (CRST) vs. thickness.
Metals 14 00168 g006
Figure 7. Prediction of plane strain forming limit (FLC0) of various steel sheets with (a) Keeler model and (b) the proposed model.
Figure 7. Prediction of plane strain forming limit (FLC0) of various steel sheets with (a) Keeler model and (b) the proposed model.
Metals 14 00168 g007
Figure 8. The determination of p from experimental FLCs: (a) SAPH440 2.3 mm, (b) QStE460TM 2.5 mm, (c) QStE600TM 3.5 mm and (d) S700MC 3.0 mm.
Figure 8. The determination of p from experimental FLCs: (a) SAPH440 2.3 mm, (b) QStE460TM 2.5 mm, (c) QStE600TM 3.5 mm and (d) S700MC 3.0 mm.
Metals 14 00168 g008
Figure 9. Prediction of complete forming limit diagram with proposed model for various steel sheets: (a) SAPH440, (b) S550MC, (c) S700MC and (d) FB780. Experimental FLCs are collected from measurement.
Figure 9. Prediction of complete forming limit diagram with proposed model for various steel sheets: (a) SAPH440, (b) S550MC, (c) S700MC and (d) FB780. Experimental FLCs are collected from measurement.
Metals 14 00168 g009
Figure 10. Prediction of complete forming limit diagram by various models for various steel sheets: (a) SAPH370, (b) QStE340TM, (c) QStE550TM, (d) 580DP, (e) 700DP and (f) Q-P-T steel. Experimental FLCs data were adapted from Refs. [20,21,47].
Figure 10. Prediction of complete forming limit diagram by various models for various steel sheets: (a) SAPH370, (b) QStE340TM, (c) QStE550TM, (d) 580DP, (e) 700DP and (f) Q-P-T steel. Experimental FLCs data were adapted from Refs. [20,21,47].
Metals 14 00168 g010
Table 1. The fitting parameters of the cubic polynomial equation.
Table 1. The fitting parameters of the cubic polynomial equation.
t (mm) A 0 A 1 A 2 A 3
2.00.491−3.8816.11−17.20
2.50.521−4.1017.15−18.62
3.00.552−4.3618.18−19.15
Table 2. The determination of CRST for QStE600TM.
Table 2. The determination of CRST for QStE600TM.
t (mm)FLC0CRST
2.00.221
2.50.241.091
3.50.261.182
5.00.281.273
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chen, W.-J.; Song, H.-W.; Chen, S.-F.; Xu, Y.; Deng, S.-Y.; Cai, Z.; Pei, X.-H.; Zhang, S.-H. A New Phenomenological Model to Predict Forming Limit Curves from Tensile Properties for Hot-Rolled Steel Sheets. Metals 2024, 14, 168. https://doi.org/10.3390/met14020168

AMA Style

Chen W-J, Song H-W, Chen S-F, Xu Y, Deng S-Y, Cai Z, Pei X-H, Zhang S-H. A New Phenomenological Model to Predict Forming Limit Curves from Tensile Properties for Hot-Rolled Steel Sheets. Metals. 2024; 14(2):168. https://doi.org/10.3390/met14020168

Chicago/Turabian Style

Chen, Wei-Jin, Hong-Wu Song, Shuai-Feng Chen, Yong Xu, Si-Ying Deng, Zheng Cai, Xin-Hua Pei, and Shi-Hong Zhang. 2024. "A New Phenomenological Model to Predict Forming Limit Curves from Tensile Properties for Hot-Rolled Steel Sheets" Metals 14, no. 2: 168. https://doi.org/10.3390/met14020168

APA Style

Chen, W. -J., Song, H. -W., Chen, S. -F., Xu, Y., Deng, S. -Y., Cai, Z., Pei, X. -H., & Zhang, S. -H. (2024). A New Phenomenological Model to Predict Forming Limit Curves from Tensile Properties for Hot-Rolled Steel Sheets. Metals, 14(2), 168. https://doi.org/10.3390/met14020168

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop