Next Article in Journal
Ultra-Precision Manufacturing Technology for Difficult-to-Machine Materials
Next Article in Special Issue
Influence of Manganese Content on Martensitic Transformation of Cu-Al-Mn-Ag Alloy
Previous Article in Journal
Study on Degradation Law and the Equivalent Thickness Model of Steel Subjected to Sulfate Corrosion
Previous Article in Special Issue
Hot Deformation Behaviour of Additively Manufactured 18Ni-300 Maraging Steel
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of Lankford Coefficients on Springback Behavior during Deep Drawing of Stainless Steel Cylinders

1
School of Material and Energy, Guangdong University of Technology, Guangzhou 510006, China
2
Guangdong Provincial Key Laboratory of Metal Forming and Forging Equipment Technology, Foshan 528300, China
*
Author to whom correspondence should be addressed.
Materials 2023, 16(12), 4321; https://doi.org/10.3390/ma16124321
Submission received: 9 May 2023 / Revised: 5 June 2023 / Accepted: 6 June 2023 / Published: 11 June 2023
(This article belongs to the Special Issue Progress in Plastic Deformation of Metals and Alloys (Second Volume))

Abstract

:
Accurate prediction of springback is increasingly required during deep-drawing formation of anisotropic stainless steel sheets. The anisotropy of sheet thickness direction is very important for predicting the springback and final shape of a workpiece. The effect of Lankford coefficients (r00, r45, r90) with different angles on springback was investigated using numerical simulation and experiments. The results show that the Lankford coefficients with different angles each have a different influence on springback. The diameter of the straight wall of the cylinder along the 45-degree direction decreased after springback, and showed a concave valley shape. The Lankford coefficient r90 had the greatest effect on the bottom ground springback, followed by r45 and then r00. A correlation was established between the springback of workpiece and Lankford coefficients. The experimental springback values were obtained by using a coordinate-measuring machine and showed good agreement with the numerical simulation results.

1. Introduction

As an important metal-forming process, sheet metal stamping is widely applied in the modern industry [1,2]. Springback is an inevitable physical phenomenon during the metal sheet-forming process [3,4,5]. The influence of springback on the accuracy and tolerance of a dimension is remarkable. The traditional trial-and-error and empirical methods for weakening springback and obtaining height-precision parts are time-consuming and expensive. The occurrence of defects, such as wrinkling, cracking, and springback, during sheet formation can be predicted with numerical simulations [6,7,8,9]. However, the predictions of springback and the final shape of the workpieces have a low accuracy rate because of the strong plastic anisotropy in the thickness direction.
A lot of research has been carried out in order to understand the influence of material properties and process parameters on springback behavior. Huang [10] analyzed the effects of different process parameters on springback during the stamping process using finite element numerical simulations. Minh [11] also used finite element simulation to analyze the effects of various factors—such as the blank holder force, friction coefficient, and blank thickness—on the springback of high-strength steel. Based on the numerical simulation results, it was evident that the blank holder force and blank thickness were the main factors affecting springback. Hashem and Roohi [12] utilized a numerical simulation to determine the effect of die and punch profile radii, as well as blank holder force on the springback and thinning percentage in the deep-drawing process of the cylindrical parts. The results show that an increased springback is observed due to an increased punch radius, and punch corner radius has been identified as the most significant effect on springback. Lajarin [13] found the blank holder force to be the most influential parameter for the springback of high-strength steel, followed by the die radius and friction conditions. Starman [14] proposed a numerical method to optimize the blank shape and tool geometry in a 3D sheet-metal-forming operation, with the effects of sheet-metal edge geometry and springback after forming and trimming being considered throughout the optimization process. Huang et al. [15] studied the defect behavior during the stamping of thin-walled semicircular shells with a bending angle via an analytical model, experiments, and a finite element (FE) simulation. The springback decreased with the increasing of blank holder force and decreasing of stamping speed. Aydın et al. [16] investigated the formability and springback behavior of dualphase (DP600) and high-strength low-alloy (HSLA) sheets bonded with laser-beam welding. The results show that springback behavior changes depending on the die angle and holding time. Where the die angle is up to 45°, the springback angle increased by 10.4%. When holding time was increased by 10-s, the springback angle decreased by 21.2%, on average. Saito et al. [17] carried out a springback experiment of 980 MPa on high-strength steel sheets with V-shaped and U-shaped bending at temperatures ranging from room temperature to 973 K. The amount of springback decreased with the temperature rise, especially at temperatures above 573 K. Springback was much reduced at lower forming speeds. The influence of stress relaxation on springback was investigated using the V-shaped bending springback test and viscoplastic stress analysis. Chang et al. [18] studied the bending springback of medium-Mn steel, a third-generation automobile steel, under different working conditions through experiment and simulation, and they also analyzed the influences of rolling direction, bending angle, and punch fillet radius on springback. The effect of the rolling direction on the springback angle was negligible. The bending angle had a positive effect on the springback angle, while the punch fillet radius had a negative effect.
During recent decades, many researchers have explored the influence of anisotropy on springback and its prediction. Ragai et al. [19] provided an experimental and computational study of springback during draw-bending of stainless steel 410. The effect of several parameters such as blank holding force, lubrication, and anisotropy on springback were discussed. Parsa et al. [20] studied the springback of hyperboloid sheet metal formation theoretically, numerically, and experimentally. Emphatically, they analyzed the influence of thickness and curvature radius on springback. The experimental results showed that the influence of material anisotropy on the forming springback of hyperboloid sheet metal is related to material parameters. Gomes et al. [21] investigated the variation in springback in high-strength steels due to material anisotropy. They analyzed and compared material models based on different yield criteria using the geometry of a standard U-shape. The results showed a discrepancy between springback predicted by the various material models and the variations in springback from 0- and 90-degree material orientation. Leu and Zhuang [22] developed a simplified approach by considering the thickness ratio, normal anisotropy, and the strain-hardening exponent to estimate the springback angle in the vee bending process for high-strength steel sheets. The numerical simulation showed that the springback ratio increased as normal anisotropy increased or as the thickness ratio and the strain-hardening exponent decreased. Verma and Haldar [23] investigated the effect of anisotropy on springback for the benchmark problem in Numisheet-2005. An analytical model was developed to cross-check the prediction from the finite element analysis. Both of the models predicted that higher anisotropy leads to more springback. Lee et al. [24] investigated the influences of the anisotropy and friction models for high-strength steel sheets during U-draw/bending and suggested the optimum selection of the models for springback simulations. Jung et al. [25] developed an elastoplastic material constitutive model with anisotropic evolution, which they then applied to a U-bending process. By comparing the measured springback angle with the predicted springback angle, they showed that the model could accurately predict the springback angle in different directions.
Stainless steel sheets are widely used in modern industry, and springback during the stamping process is a common problem. Various process parameters can affect springback during the deep-drawing process of anisotropic stainless steel sheets. Although many scholars have conducted investigations on the springback problem during sheet metal formation using the finite element numerical simulation method and experimental methods, the research on springback of anisotropic sheet metal is still mostly confined to the forming of V/U-shaped parts, and there is little research on other common forming parts, for example, cylindrical cups. The influence of Lankford coefficients with different angles on springback during the cylinder deep-drawing process has not been clearly researched.
In this paper, a cylinder deep-drawing process with anisotropic stainless steel sheets was simulated based on the Barlat–Lian 1989 anisotropy yield criterion [26] by using Dynaform 5.9 software, and was used to predict springback. The Taguchi and ANOVA techniques were utilized to establish the correlation between springback at different angles from the rolling direction and Lankford coefficients (r00, r45, r90) of 304 stainless steel. The ANOVA showed that the Lankford coefficient had a significant effect on springback. This research shows that each Lankford coefficient has an obvious influence on springback in diffident angles from the rolling direction by using the experimental and numerical simulation.

2. Finite Element Simulation (FEM) Analysis

2.1. FEM Simulation Procedure

In this paper, the finite element numerical simulation was carried out on Dynaform. Dynaform software is a special piece of software jointly developed by ETA and LSTC for numerical simulation of sheet metal formation. It is a combination of LS-DYNA solver and ETA/FEMB front and back processor, and it is one of the most popular CAE tools for sheet metal formation and die design. Figure 1 shows the cylinder deep-drawing die; the actual object of the model is thecooking pot. The dimensions of the blank, die, punch, and blank holder are given in Table 1. One quarter of the 3D numerical model can be applied to the FEM model. The simulations require a large amount of computational time if they are not simplified, but they can provide a greater degree of precision. In this paper, a complete 3D numerical model was used. The 3D numerical model is shown in Figure 2. Table 2 shows the Lankford coefficients of metal sheets in different rolling directions with two horizontal factors set. The other mechanical properties of the materials were imported from the materials library in Dynaform. The punch and die were set as rigid, and the velocity of the punch was set at 2000 mm/s. The friction coefficient between the tools and the blank were set to 0.125. The contact-one-way surface-to-surface mode was employed to determine the friction type, and the adaptive meshing method was adopted to mesh the geometry model [27]. The full integrated planar shell was used, and the element type was defined as the time-efficient full-order integral Belytschko–Tsay shell element. This allowed for the adoption of four-point integration in order to avoid the appearance of “hourglassing” mode. A dynamic explicit algorithm was used to calculate the forming process. The implicit algorithm was applied to calculate the springback process.
The material was modeled as an elastic–plastic material. The anisotropic characteristic was described by the Barlat–Lian 1989 anisotropic yield criterion [28]. The Barlat–Lian 1989 anisotropic yield criterion and the Hosford series’ yield criterion were used to analyze the plastic flow law of the drawing process [29,30,31]. Three stress–strain curves were obtained from the tensile test for the model material, as shown in Figure 3. The different curves were determined according to the ratio of the Lankford coefficients in each direction of the actual material.

2.2. Taguchi Technique

The Taguchi technique was applied to the design scheme of the numerical simulation [32]. The two levels of the three-parameter orthogonal design, considering interactions (27), are presented in Table 3. The springback of different angles from the rolling direction was the process response. In order to understand the influence of Lankford coefficients, the ANOVA technique was applied to illustrate the degree of significance of each Lankford coefficient, including interactions.

2.3. Measurement Set-Up

Figure 4 shows the typical shape characteristics and measurement locations of the cylindrical cup. After formation, the workpiece was measured using CMM. Angles (α) were measured every 45 degrees from the rolling direction, and diameters were measured every 15 mm in the five sections along the height. A diagrammatic sketch of angles from the rolling direction is shown in Figure 4.

2.4. Formation Analysis

The forming limit diagram (FLD) and thickness change diagram can intuitively show the dynamic drawing process of the sheet metal and predict the formation of defects, such as cracking and wrinkling, and the thickness distribution of the sheet metal [33]. Figure 5a shows the forming limit diagram of the cylindrical cup after deep-drawing formation, Figure 5b shows the forming limit diagram of the cylindrical cup after springback, and Figure 5c is the cloud diagram of springback change in the cylindrical cup after springback calculation. It can be seen that the cylindrical cup fluctuates after springback with different degrees in the flange. The springback is apparent at 0°, 45°, and 90° positions, which shows a cyclical trend of first decreasing and then increasing along the rolling direction. The straight wall of the cylinder also showed uneven springback.
Due to the uneven springback deformation of the straight wall of the cylindrical cup, the sections with heights of 45 mm and 60 mm were selected for measurement, and 120 coordinate points were measured for each section. The difference between coordinate values of data points before and after springback was calculated. The cross-section difference point cloud diagrams are shown in Figure 6. The co-ordinates only represent the position of data points on the section of the cylindrical drawing section. The distance between each point and the origin represents the springback value. It can be seen that the springback difference between the two heights is similar, and is in the range of 0.150–25 mm. Within the angle of 0–45° from the rolling direction, the springback difference firstly decreases, and then it increases. At the position of the maximum plastic strain value of r45, that is, at the positions 45°, 135°, 225°, and 315° from the rolling direction, the springback difference of the cylinder drawing part reaches its maximum value.

2.5. Stress–Strain Analysis

The straight wall of the cylindrical cup is an area of force transmission during deep drawing, and no more plastic deformation occurs. The straight wall experiences a single axial tensile stress. There is a small amount of axial elongation and deformation. The state of stress and strain during deep drawing is shown in Figure 7. The first principal stress and strain of the model of the straight wall model’s middle layer was extracted to analyze the reasons for uneven springback.
The stress–strain analysis diagrams of cylindrical deep drawing at the heights of 45 mm and 60 mm are shown in Figure 8. The stress–strain data of 60 points on the circumference of the straight wall were extracted, and the red circle represents the average stress–strain value of all points. It can be seen that, at the height of 45 and 60 mm, the first principal stress was greater than the other two directions at the position of 45° from the rolling direction, while the first principal strain was smaller than the other two directions.
The stress–strain values for the three rolling directions of 0°, 45°, and 90° were compared and analyzed, and the results are shown in Table 4. The stress at 45° at the height of 45 mm is 23% higher, and the stress at 45° at the height of 60 mm is 19.37% higher than that of the other rolling directions. This is because the hardening curves are for different rolling directions. The value of the hardening curve at 45° from the rolling direction was larger, and the stress value required during deep drawing was larger. The strain in the 45° direction was smaller and contained more elastic stress in the deformation process, resulting in greater springback deformation after unloading.

2.6. Boundary Inflow Analysis

The diagram of inflow of the cylindrical cup boundary material is shown in Figure 9. It shows a cyclical trend of first increasing and then decreasing between 0° and 90° from the rolling direction. At the positions 45°, 135°, 215°, and 315°, there was a larger inflow, and the maximum value was 40.76 mm. The Lankford coefficients in the 45° direction were greater than those for the 0° and 90° directions. When the Lankford coefficients were large, the deformation resistance of the flange of the metal sheet was reduced, and the material flowed more easily. The flow stress value in the 45° direction was large, so the inflow of material was larger. The flow stress in the 0° and 90° directions was smaller, elongation deformation was easier, and the inflow was smaller. This may be one of the reasons for the greater springback difference in this direction.
The sheet firstly underwent elastic deformation, and then plastic deformation occurred after the stress exceeded the flow stress during the deformation process. After unloading, the internal stress was redistributed, and then springback occurred. The deformation and plastic deformation of the cylindrical cup drawing in the 45° direction was less severe than the other two directions. The circumferential stress in the 45° direction was relatively large, resulting in larger springback deformation at this position. The diameter of the cylinder after springback in the 45° direction was smaller, showing a greater springback difference.

3. Experimental Procedures

3.1. Experimental Set-Up

Two different stainless steel sheets with the same thickness were selected for the experimental test. Based on the previous experimental tests [34,35], strong anisotropic properties were present in the two materials. Lankford coefficients of r00, r45, and r90 are listed in Table 2. For cylinder deep drawing, circular blanks with a diameter of 315 mm and a thickness of 0.6 mm were prepared. Figure 10 shows the drawing die designed and fabricated based on the simulation model.

3.2. Experimental Results

After calculating the weighted average of the co-ordinates of two types of stainless steel workpieces at different heights, the radius values at different angles were obtained, as shown in Figure 11. The average values of the cross-section point cloud can be compared at the height of 30 mm, 45 mm, and 60 mm. The valley shape of the depression was clearer and more obvious in the material with the larger Lankford coefficients at the positions of 45°, 135°, 225°, and 315°. At the height of 15 mm near the bottom of the cylinder, the low-Lankford-coefficient material showed a more rounded cross-section. The high-Lankford-coefficient material showed a less rounded cross-section after drawing, and it was close to an ellipse along the long axis of the rolling direction and the short axis perpendicular to the rolling direction. At a height of 75 mm near the flange, the sections of both kinds of stainless steel showed an elliptical shape after deep drawing and springback. The experimental results are in good agreement with the FEM simulation results.
Table 5, Table 6, Table 7 and Table 8 compare the radius values along the three rolling directions at five section heights obtained from numerical simulation and experiments. It can be seen that the difference in radius between different rolling directions is more obvious in the simulation, and it was the largest at the heights of 45 mm and 60 mm. In the material with high-level Lankford coefficients, the differences reached 0.433 mm and 0.318 mm, respectively. In the material with low-level Lankford coefficients, the differences reached 0.387 mm and 0.32 mm, respectively. The experimental difference between radii in different rolling directions was close to the simulation result in the material with high-level Lankford coefficients. The difference at the heights of 30 mm and 45 mm reached 0.158 mm and 0.204 mm, respectively. The experimental radius difference between different rolling directions in material with low-level Lankford coefficients was small, reaching 0.127 mm and 0.08 mm at the heights of 45 mm and 60 mm, respectively.
Figure 12 shows experimental measurements of diameter at five sections along the height. They have a similar trend. The section at the height of 15 mm was close to the radius of the punch nose. The section at the height of 75 mm was similar to the radius of the die shoulder. The fillet radius had a great influence on the workpiece diameter. The cross-section of the cylinder after deep drawing showed an oval shape after springback. The diameters at the height of 30, 45, and 60 mm showed a similar trend. The results showed that the r45 Lankford coefficient is the maximum value. In addition, as the Lankford coefficient increased, the diameter decreased.

4. Results and Discussions

4.1. Application of ANOVA

Table 9 shows the results of springback prediction by FEM simulation, which shows that the springback of every angle from the rolling direction is without symmetrical characteristics.
To investigate the degree of significance of the Lankford coefficients, the ANOVA technique was used to analyze the springback. The mean overall value S / N S / N ¯ is expressed as Equation (1), where k is the number of simulations. The range of two levels (SRj) is shown in Equation (2). The sum of squares owing to the variations of the overall mean (SS) and the mean of the Lankford coefficients with interactions (SSj) are expressed as Equations (3) and (4), respectively. The percentage values (%p-Valuej) were calculated using Equation (5), which is generally applied when measuring the degree of significance of each Lankford coefficient [36].
S / N ¯ = 1 8 k = 1 8 S / N k
S R j = j = 1 7 S / N 1 j S / N 2 j
S S = i = 1 8 S / N i j S / N ¯ 2
S S j = j = 1 7 S / N i j S / N ¯ 2
% p V a l u e j = S S j S S × 100

4.2. Effects of Process Parameters on Springback

The results of the range analysis and variance analysis are shown in Table 10. It is revealed that the influence of Lankford coefficients on springback is different at different angles. The source of r00 had a critical effect on the springback at θ000 from the rolling direction, the r45 is the key factor causing springback at θ045 and θ315, and r90 is key at θ135, θ180, and θ225. Furthermore, r45 × r00 has significant values for springback at θ090 because of interactions with Lankford coefficients. Meanwhile, r90 × r45 is the key factor of springback at θ270. The measurement error also easily affected the range analyses, and so the ANOVA technique was used to analyze the springback. The ANOVA results shown in Table 9 correspond well with the range analysis results. Based on these analysis results, it has been found that the interactions of Lankford coefficients at different angles from the rolling direction have a clear effect on the springback.
Figure 13a–h shows the sensitivity analysis of the effect of Lankford coefficients on springback. When the interactions of the Lankford coefficients were not considered, the amount of springback decreased with r90 and r45 at all angles, and the amount of springback decreased with increasing r00, except at the angle of θ000. When considering the interactions, the amount of springback increased with increasing r90 × r45. When the r45 and r90 increased simultaneously, the interaction of r90 × r45 hindered springback and caused it to decrease. The amount of springback increased with increasing r45 × r00, except at the angle of θ180. The amount of springback decreased with increasing r90 × r00, except at the angles of θ045 and θ180. The results show that the influences of the Lankford coefficient on springback at different angles are interrelated and interact with each other.

4.3. Comparing the FEM Simulation and Experimental Results

The experimental results of two kinds of 304 stainless steel were measured using CMM, as shown in Table 11. The No. 1 and No. 8 FEM simulation results are shown in Table 12. Figure 14 shows the comparison of an average bottom fillet of the FEM simulation and experimental results at different angles from the rolling direction. The material with high-level Lankford coefficients had a larger amount of springback at 0, 90, and 270 degrees from the rolling direction. The experimental results have good agreement with the FEM simulation results, with the bottom fillet showing a similar trend. The springback of the cylinder bottom fillet occurred along the rolling direction, and there was an increasing trend with every 45° decrease.
The comparison between the bottom fillets of materials with high- and low-level Lankford coefficients is shown in Figure 15. The experimental results have good agreement with the FEM simulation results, showing a similar trend in springback. The amount of the bottom fillet decreased with an increase in the Lankford coefficient at all locations, except for 0, 90, and 270 degrees. The trend is more apparent in the FEM simulation results.

5. Conclusions

The influence of Lankford coefficients on stainless steel cylindrical cups was investigated using both experiments and numerical simulation. The conclusions are as follows:
The simulation and experimental results show that the Lankford coefficients had an obvious effect on the cross-section diameter. The flow velocity of the blank was different because of the anisotropy of the metal sheet, which makes the stress–strain values accumulate in different directions during the deep-drawing process, and finally causes a clear difference in springback. The simulated springback value for the straight wall was between 0.15 mm and 0.25 mm. The maximum springback value was at the position of 45° from the rolling direction. Specifically, the diameters at different height sections were related to the Lankford coefficients at different angles from the rolling direction, which were characterized by a concave valley in the 45 degree direction of the straight wall. The radius difference between the 45 degree rolling direction and the other two directions at each section height was between 0.1 mm and 0.3 mm.
The ANOVA results illustrated the influence of Lankford coefficients on the springback of the bottom fillet. The Lankford coefficient has different levels of effects on springback depending on the angle from the rolling direction. The 90-degree angle had the greatest influence, followed by the 45-degree, with 0 degrees having the least influence. The experimental results showed a similar trend to the simulation results. In addition, the springback of the bottom fillet decreased with the increasing overall Lankford coefficients.
The combination of the FEM simulation, the ANOVA technique, and the experimental study of the cylinder deep-drawing process is an effective method for studying the influence of Lankford coefficients on springback and predicting the final shape with high precision.
In this paper, the study of the effect of Lankford coefficients on springback of a cylindrical cup during the deep-drawing process remained at the macroscopic stage, and further analysis on the microscopic aspects was not carried out. The mechanism by which metal anisotropy influences the springback of a cylinder needs to be explored further. Although some characteristic rules regarding springback and cylindrical cup properties during deep drawing were obtained in this study, the analysis of the specific degree of influence of Lankford coefficients on springback properties is still in the preliminary stage. Therefore, the quantitative analysis of the influence of Lankford coefficients on the springback of the cylindrical cup during deep drawing will remain the focus of future research. Based on the previous research on the cylindrical cup, the springback of large complex thin-walled parts in deep drawing will be further explored.

Author Contributions

Conceptualization, F.W.; Methodology, Software, Visualization, and Writing, Y.H. and C.H.; Validation, Z.H.; Supervision, Project Administration, and Funding Acquisition, Z.Z. and F.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No. 52175294 and No. 51705085), the State Engineering Research Center for Metallic Materials Net-shape Processing, South China University of Technology (Grant No. 2020011), and the Core Technology Research Foundation of Foshan City (Grant No. 1920001001369).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy.

Acknowledgments

The authors gratefully acknowledge the financial support from the National Natural Science Foundation of China (Grant No. 52175294 and No. 51705085), the State Engineering Research Center for Metallic Materials Net-shape Processing, South China University of Technology (Grant No. 2020011), and the Core Technology Research Foundation of Foshan City (Grant No. 1920001001369).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Trzepieciński, T. Recent Developments and Trends in Sheet Metal Forming. Metals 2020, 10, 779. [Google Scholar] [CrossRef]
  2. Centeno, G.; Silva, M.B. Tube and Sheet Metal Forming Processes and Applications. Metals 2022, 12, 553. [Google Scholar] [CrossRef]
  3. Guo, S.; Tian, C.; Pan, H.; Tang, X.; Han, L.; Wang, J. Research on the Springback Behavior of 316Ln Stainless Steel in Micro-Scale Bending Processes. Materials 2022, 15, 6373. [Google Scholar] [CrossRef] [PubMed]
  4. Gia Hai, V.; Thi Hong Minh, N.; Nguyen, D.T. A Study on Experiment and Simulation to Predict the Spring-Back of Ss400 Steel Sheet in Large Radius of V-Bending Process. Mater. Res. Express 2020, 7, 16562. [Google Scholar] [CrossRef]
  5. Kut, S.; Pasowicz, G.; Stachowicz, F. On the Springback and Load in Three-Point Air Bending of the Aw-2024 Aluminium Alloy Sheet with Aw-1050a Aluminium Cladding. Materials 2023, 16, 2945. [Google Scholar] [CrossRef]
  6. Solfronk, P.; Sobotka, J.; Kolnerova, M.; Zuzanek, L. Spring-Back Prediction for Stampings from the Thin Stainless Sheets. MM Sci. J. 2016, 2016, 1090–1094. [Google Scholar] [CrossRef]
  7. Ha, G.X.; Oliveira, M.G.; Andrade-Campos, A.; Manach, P.Y.; Thuillier, S. Prediction of Coupled 2D and 3D Effects in Springback of Copper Alloys After Deep Drawing. Int. J. Mater. Form. 2021, 14, 1171–1187. [Google Scholar] [CrossRef]
  8. Ul Hassan, H.; Maqbool, F.; Güner, A.; Hartmaier, A.; Ben Khalifa, N.; Tekkaya, A.E. Springback Prediction and Reduction in Deep Drawing Under Influence of Unloading Modulus Degradation. Int. J. Mater. Form. 2016, 9, 619–633. [Google Scholar] [CrossRef]
  9. Trzepiecinski, T.; Lemu, H.G. Effect of Computational Parameters on Springback Prediction by Numerical Simulation. Metals 2017, 7, 380. [Google Scholar] [CrossRef] [Green Version]
  10. Huang, Y.; Chou, I.; Jiang, C.; Wu, Y.; Lee, S. Finite Element Analysis of Dental Implant Neck Effects on Primary Stability and Osseointegration in a Type Iv Bone Mandible. Bio-Med. Mater. Eng. 2014, 24, 1407–1415. [Google Scholar] [CrossRef]
  11. Minh, N. Effect of Forming Parameters on Springback of Advanced High Strength Steel DP800. Appl. Mech. Mater. 2014, 703, 182–186. [Google Scholar] [CrossRef]
  12. Hashemi, S.J.; Roohi, A.H. Minimizing Spring-Back and Thinning in Deep Drawing Process of St14 Steel Sheets. Int. J. Interact. Des. Manuf. 2022, 16, 381–388. [Google Scholar] [CrossRef]
  13. Lajarin, S.F.; Marcondes, P.V. Influence of Process and Tool Parameters on Springback of High-Strength Steels. Proceedings of the Institution of Mechanical Engineers, Part B. J. Eng. Manuf. 2014, 229, 295–305. [Google Scholar] [CrossRef]
  14. Starman, B.; Cafuta, G.; Mole, N. A Method for Simultaneous Optimization of Blank Shape and Forming Tool Geometry in Sheet Metal Forming Simulations. Metals 2021, 11, 544. [Google Scholar] [CrossRef]
  15. Huang, X.; Guan, B.; Zang, Y.; Wang, B. Investigation of Defect Behavior During the Stamping of a Thin-Walled Semicircular Shell with Bending Angle. J. Manuf. Process. 2023, 87, 231–244. [Google Scholar] [CrossRef]
  16. Aydın, K.; Karaağaç, İ.; Uluer, O. The Formability and Springback Characterization of Laser-Welded Dp–Hsla Sheets. Appl. Phys. A 2019, 125, 525. [Google Scholar] [CrossRef]
  17. Saito, N.; Fukahori, M.; Hisano, D.; Hamasaki, H.; Yoshida, F. Effects of Temperature, Forming Speed and Stress Relaxation on Springback in Warm Forming of High Strength Steel Sheet. Procedia Eng. 2017, 207, 2394–2398. [Google Scholar] [CrossRef]
  18. Chang, Y.; Wang, N.; Wang, B.T.; Li, X.D.; Wang, C.Y.; Zhao, K.M.; Dong, H. Prediction of Bending Springback of the Medium-Mn Steel Considering Elastic Modulus Attenuation. J. Manuf. Process. 2021, 67, 345–355. [Google Scholar] [CrossRef]
  19. Ragai, I.; Lazim, D.; Nemes, J.A. Anisotropy and Springback in Draw-Bending of Stainless Steel 410: Experimental and Numerical Study. J. Mater. Process. Technol. 2005, 166, 116–127. [Google Scholar] [CrossRef]
  20. Parsa, M.H.; Nasher Al Ahkami, S.; Pishbin, H.; Kazemi, M. Investigating Spring Back Phenomena in Double Curved Sheet Metals Forming. Mater. Des. 2012, 41, 326–337. [Google Scholar] [CrossRef]
  21. Gomes, C.; Onipede, O.; Lovell, M. Investigation of Springback in High Strength Anisotropic Steels. J. Mater. Process. Technol. 2005, 159, 91–98. [Google Scholar] [CrossRef]
  22. Leu, D.; Zhuang, Z. Springback Prediction of the Vee Bending Process for High-Strength Steel Sheets. J. Mech. Sci. Technol. 2016, 30, 1077–1084. [Google Scholar] [CrossRef]
  23. Verma, R.K.; Haldar, A. Effect of Normal Anisotropy on Springback. J. Mater. Process. Technol. 2007, 190, 300–304. [Google Scholar] [CrossRef]
  24. Lee, J.; Barlat, F.; Lee, M. Constitutive and Friction Modeling for Accurate Springback Analysis of Advanced High Strength Steel Sheets. Int. J. Plast. 2015, 71, 113–135. [Google Scholar] [CrossRef]
  25. Jung, J.; Jun, S.; Lee, H.; Kim, B.; Lee, M.; Kim, J.H. Anisotropic Hardening Behaviour and Springback of Advanced High-Strength Steels. Metals 2017, 7, 480. [Google Scholar] [CrossRef] [Green Version]
  26. Barlat, F.; Lian, J. Plastic behavior and stretch-ability of sheet metals. Part I A yield function for orthotropic sheet under plane stress conditions. Int. J. Plast. 1989, 5, 51–66. [Google Scholar] [CrossRef]
  27. Zhang, D.; Bai, D.; Liu, J.; Guo, Z.; Guo, C. Formability Behaviors of 2a12 Thin-Wall Part Based on Dynaform and Stamping Experiment. Compos. Part B Eng. 2013, 55, 591–598. [Google Scholar] [CrossRef]
  28. Panich, S.; Uthaisangsuk, V.; Suranuntchai, S.; Jirathearanat, S. Investigation of Anisotropic Plastic Deformation of Advanced High Strength Steel. Mater. Sci. Eng. A 2014, 592, 207–220. [Google Scholar] [CrossRef]
  29. Zhang, F.F.; Chen, J.S.; Chen, J.; Huang, X.Z.; Lu, J. Review on development and experimental validation for anisotropic yield criterions. Adv. Mech. 2012, 42, 68–80. [Google Scholar]
  30. Hosford, W.F. A Generalized Isotropic Yield Criterion. J. Appl. Mech. 1972, 39, 607–609. [Google Scholar] [CrossRef]
  31. Huang, M.; Man, C. A Generalized Hosford Yield Function for Weakly-Textured Sheets of Cubic Metals. Int. J. Plast. 2013, 41, 97–123. [Google Scholar] [CrossRef]
  32. Eklarkar, S.V.; Nandedkar, V.M. Fe Simulation and Optimization of Drawbead Parameters Using Taguchi Method for Hemispherical Cup. Mater. Today Proc. 2021, 44, 4709–4716. [Google Scholar] [CrossRef]
  33. Shinge, V.R.; Dabade, U.A. Experimental Investigation on Forming Limit Diagram of Mild Carbon Steel Sheet. Procedia Manuf. 2018, 20, 141–146. [Google Scholar] [CrossRef]
  34. Ailinei, I.; Galatanu, S.; Marsavina, L. Influence of Anisotropy on the Cold Bending of S600Mc Sheet Metal. Eng. Fail. Anal. 2022, 137, 106206. [Google Scholar] [CrossRef]
  35. Džoja, M.; Cvitanić, V.; Safaei, M.; Krstulović-Opara, L. Modelling the Plastic Anisotropy Evolution of Aa5754-H22 Sheet and Implementation in Predicting Cylindrical Cup Drawing Process. Eur. J. Mech. A Solids 2019, 77, 103806. [Google Scholar] [CrossRef]
  36. Thipprakmas, S. Application of Taguchi Technique to Investigation of Geometry and Position of V-Ring Indenter in Fine-Blanking Process. Mater. Des. 2010, 31, 2496–2500. [Google Scholar] [CrossRef]
Figure 1. Geometry of drawing die.
Figure 1. Geometry of drawing die.
Materials 16 04321 g001
Figure 2. Model in FEM.
Figure 2. Model in FEM.
Materials 16 04321 g002
Figure 3. The stress–strain curves in different angles to the rolling direction: (a) high-level Lankford coefficient material hardening curve and (b) low-level Lankford coefficient material hardening curve.
Figure 3. The stress–strain curves in different angles to the rolling direction: (a) high-level Lankford coefficient material hardening curve and (b) low-level Lankford coefficient material hardening curve.
Materials 16 04321 g003
Figure 4. The Schematic diagram of the measurement position: (a) deformation and measurement locations and (b) diagrammatic sketch of angles from the rolling direction.
Figure 4. The Schematic diagram of the measurement position: (a) deformation and measurement locations and (b) diagrammatic sketch of angles from the rolling direction.
Materials 16 04321 g004
Figure 5. Forming limit diagram of the cylindrical cup and the cloud diagram of springback: (a) limit diagram of the cylindrical cup, (b) Limit diagram of the cylindrical cup after springback, and (c) cloud map of the cylindrical cup after springback.
Figure 5. Forming limit diagram of the cylindrical cup and the cloud diagram of springback: (a) limit diagram of the cylindrical cup, (b) Limit diagram of the cylindrical cup after springback, and (c) cloud map of the cylindrical cup after springback.
Materials 16 04321 g005
Figure 6. The springback difference point cloud of the cross-section of the cylindrical cup: (a) high–level Lankford coefficient material and (b) low–level Lankford coefficient material.
Figure 6. The springback difference point cloud of the cross-section of the cylindrical cup: (a) high–level Lankford coefficient material and (b) low–level Lankford coefficient material.
Materials 16 04321 g006
Figure 7. The stress and strain state of straight wall area.
Figure 7. The stress and strain state of straight wall area.
Materials 16 04321 g007
Figure 8. Stress and strain states on sections with different heights: (a) stress values at the 45 mm height position, (b) strain values at the 45 mm height position, (c) stress values at the 60 mm height position, and (d) strain values at the 60 mm height position.
Figure 8. Stress and strain states on sections with different heights: (a) stress values at the 45 mm height position, (b) strain values at the 45 mm height position, (c) stress values at the 60 mm height position, and (d) strain values at the 60 mm height position.
Materials 16 04321 g008
Figure 9. The boundary inflow of material and the trend of material inflow: (a) schematic diagram of the material boundary inflow and (b) numerical change diagram.
Figure 9. The boundary inflow of material and the trend of material inflow: (a) schematic diagram of the material boundary inflow and (b) numerical change diagram.
Materials 16 04321 g009
Figure 10. Die assembly for experimental drawing.
Figure 10. Die assembly for experimental drawing.
Materials 16 04321 g010
Figure 11. Average vertical wall shell’s dimension of the experiments and simulations: (a) high-level Lankford coefficient material simulation results, (b) low-level Lankford coefficient material simulation results, (c) high-level Lankford coefficient material experimental results, and (d) low-level Lankford coefficient material experimental results. Note: Black: H = 15 mm Red: H = 30 mm Green: H = 45 mm Blue: H = 60 mm Light Blue: H = 75 mm.
Figure 11. Average vertical wall shell’s dimension of the experiments and simulations: (a) high-level Lankford coefficient material simulation results, (b) low-level Lankford coefficient material simulation results, (c) high-level Lankford coefficient material experimental results, and (d) low-level Lankford coefficient material experimental results. Note: Black: H = 15 mm Red: H = 30 mm Green: H = 45 mm Blue: H = 60 mm Light Blue: H = 75 mm.
Materials 16 04321 g011aMaterials 16 04321 g011b
Figure 12. Diameter of different sections of experimental measurements: (a) the material with high-level Lankford coefficients and (b) the material with low-level Lankford coefficients.
Figure 12. Diameter of different sections of experimental measurements: (a) the material with high-level Lankford coefficients and (b) the material with low-level Lankford coefficients.
Materials 16 04321 g012
Figure 13. The sensitivity analysis of the effect of Lankford coefficients on springback.
Figure 13. The sensitivity analysis of the effect of Lankford coefficients on springback.
Materials 16 04321 g013aMaterials 16 04321 g013b
Figure 14. Comparison of the average of experimental results and the FEM simulation results.
Figure 14. Comparison of the average of experimental results and the FEM simulation results.
Materials 16 04321 g014
Figure 15. Comparison of the variation of angles between the FEM simulation results and the experimental results.
Figure 15. Comparison of the variation of angles between the FEM simulation results and the experimental results.
Materials 16 04321 g015
Table 1. Basic geometrical parameters.
Table 1. Basic geometrical parameters.
ParameterDimension in mm
blank size diameter (BD)315
blank thickness (t)0.6
punch diameter (PD)180
punch nose radius (rp)8
die shoulder radius (rd)4
die diameter (DD)181.32
radial clearance between punch and die (wc)0.66
height of drawing (h)80
Table 2. Test factors and their levels.
Table 2. Test factors and their levels.
Factorsr
Level r00r45r90
Low-level0.991.260.92
High-level1.071.361.03
Table 3. Experimental design of orthogonal considering interactions (27) for FEM simulation.
Table 3. Experimental design of orthogonal considering interactions (27) for FEM simulation.
Simulation
No.
Factors
r90r45r90 × r45r00 r90 ×
r00
r45 × r00Error
1LowLowLowLowLowLowLow
2LowLowLowHighHighHighHigh
3LowHighHighLowLowHighHigh
4LowHighHighHighHighLowLow
5HighLowHighLowHighLowHigh
6HighLowHighHighLowHighLow
7HighHighLowLowHighHighLow
8HighHighLowHighLowLowHigh
Table 4. The stress and strain in three directions at the heights of 45 mm and 60 mm of high-level Lankford coefficient material.
Table 4. The stress and strain in three directions at the heights of 45 mm and 60 mm of high-level Lankford coefficient material.
Direction45°90°Difference
(Max–Min)
Height
H45
Stress/Mpa770.475947.854782.773177.379 (23%)
Strain0.4070.3490.3960.058
H60
Stress/Mpa973.371162997.662188.63 (19.37%)
Strain0.4680.3680.4630.1
Table 5. The radius values of the simulations with high-level Lankford coefficients on the three directions.
Table 5. The radius values of the simulations with high-level Lankford coefficients on the three directions.
Direction0°/mm45°/mm90°/mmDifference
/mm
Height/mm
H1589.56189.49289.5370.068
H3089.98789.74089.9830.246
H4590.50690.07290.4860.434
H6090.61790.29990.5920.318
H7590.59990.57590.5170.024
Table 6. The radius values of the simulations with low-level Lankford coefficients on the three directions.
Table 6. The radius values of the simulations with low-level Lankford coefficients on the three directions.
Direction0°/mm45°/mm90°/mmDifference
/mm
Height/mm
H1589.54189.49089.4950.051
H3089.88489.71089.8830.174
H4590.42890.04190.4040.387
H6090.60390.28390.5570.320
H7590.65290.58590.5150.067
Table 7. The radius values of the experiments with high-level Lankford coefficients on the three directions.
Table 7. The radius values of the experiments with high-level Lankford coefficients on the three directions.
Direction0°/mm45°/mm90°/mmDifference
/mm
Height/mm
H1590.54690.50690.452−0.007
H3090.78190.56390.6600.158
H4590.84190.59790.7600.204
H6090.81190.64290.7080.118
H7590.90590.83790.7750.003
Table 8. The radius values of the experiments with low-level Lankford coefficients on the three directions.
Table 8. The radius values of the experiments with low-level Lankford coefficients on the three directions.
Direction0°/mm45°/mm90°/mmDifference
/mm
Height/mm
H1590.51790.49890.5160.019
H3090.69490.61090.6450.060
H4590.78190.62890.7290.127
H6090.78090.65490.6970.085
H7590.89590.82690.7670.005
Table 9. The results of springback prediction by FEM simulation.
Table 9. The results of springback prediction by FEM simulation.
The Amount of Springback/(°)
θ000θ045θ090θ135θ180θ225θ270θ315Difference (Max–Min)
190.20590.05290.31390.62390.62990.64589.99990.3030.44
290.29490.74890.64590.83890.61790.72990.36390.3150.544
390.25590.56690.77991.02190.44590.93190.60190.4070.766
490.21790.24290.02790.34290.65490.41290.26890.1120.627
590.20590.51890.37990.52990.62690.51490.33890.5770.421
690.29590.50790.65290.72390.17590.48190.41890.4470.428
790.02790.1490.24690.00690.06989.96890.14889.8840.362
890.22589.96890.30390.18890.24990.16190.27990.0580.257
average90.21590.37590.38590.53490.43390.48090.30290.2630.319
Table 10. Results of range analysis and variance analysis.
Table 10. Results of range analysis and variance analysis.
Range AnalysisVariance Analysis (The Key Factor and Contribution Value)
θ000r00r00(28.79%), r45(18.94%)
θ045r45r45(37.53%), r45 × r00(23.21%)
θ090r45 × r00r45 × r00(54.46%), r45(17.11%)
θ135r90r90(29.35%), r45(20.66%)
θ180r90r90(48.00%), r45 × r00(23.13%)
θ225r90r90(47.89%), r45(15.19%)
θ270r90 × r45r90 × r45(38.81%), r45 × r00(23.17%)
θ315r45r45(47.64%), r90 × r45(33.00%)
Table 11. Average experimental results of two kinds of stainless steel.
Table 11. Average experimental results of two kinds of stainless steel.
MaterialThe Amount of Springback/(°)
θ000θ045θ090θ135θ180θ225θ270θ315Difference (Max–Min)
High–level Lankford coefficient90.82390.40190.96290.50490.94290.62090.68390.4600.561
Low–level Lankford coefficient90.75690.52990.84590.70091.00990.79290.67190.5650.48
Difference (Max–Min)0.067−0.1280.117−0.196−0.389−0.1720.012−0.105
Table 12. FEM simulation results of springback.
Table 12. FEM simulation results of springback.
Simulation
No.
The Amount of Springback/(°)
θ000θ045θ090θ135θ180θ225θ270θ315Difference (Max–Min)
190.20590.05290.31390.62390.62990.64589.99990.3030.44
890.22589.96890.30390.18890.24990.16190.27990.0580.257
Difference (Max–Min)0.02−0.084−0.01−0.435−0.38−0.4840.28−0.245
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

Wu, F.; Hong, Y.; Zhang, Z.; Huang, C.; Huang, Z. Effect of Lankford Coefficients on Springback Behavior during Deep Drawing of Stainless Steel Cylinders. Materials 2023, 16, 4321. https://doi.org/10.3390/ma16124321

AMA Style

Wu F, Hong Y, Zhang Z, Huang C, Huang Z. Effect of Lankford Coefficients on Springback Behavior during Deep Drawing of Stainless Steel Cylinders. Materials. 2023; 16(12):4321. https://doi.org/10.3390/ma16124321

Chicago/Turabian Style

Wu, Fei, Yihao Hong, Zhengrong Zhang, Chun Huang, and Zhenrong Huang. 2023. "Effect of Lankford Coefficients on Springback Behavior during Deep Drawing of Stainless Steel Cylinders" Materials 16, no. 12: 4321. https://doi.org/10.3390/ma16124321

APA Style

Wu, F., Hong, Y., Zhang, Z., Huang, C., & Huang, Z. (2023). Effect of Lankford Coefficients on Springback Behavior during Deep Drawing of Stainless Steel Cylinders. Materials, 16(12), 4321. https://doi.org/10.3390/ma16124321

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