Research Article  Open Access
Optimization of Process Parameters of Stamping Forming of the Automotive Lower Floor Board
Abstract
There are many process parameters which have great effect on the forming quality of parts during automobile panel stamping forming process. This paper took automotive lower floor board as the research object; the forming process was analyzed by finite element simulation using Dynaform. The influences of four main process parameters including BHF (blank holder force), die corner radius, friction coefficient, and die clearance on the maximum thinning rate and the maximum thickening rate were researched based on orthogonal experiment. The results show that the influences of each value of various factors on the target are not identical. On this basis, the optimization of the four parameters was carried out, and the high quality product was obtained and the maximum thinning rate and maximum thickening rate were effectively controlled. The results also show that the simulation analysis provides the basis for the optimization of the forming process parameters, and it can greatly shorten the die manufacturing cycles, reduce the production costs, and improve the production efficiency.
1. Introduction
Automotive lower floor board has an important effect on the transportation process of vehicle; the performance of lower floor determines the stability and safety of the whole vehicle [1, 2]. Automobile lower floor, engine cover, driving cab, and various specialshaped surface and the internal sheet parts of the body all belong to the automotive panels. Automobile panel is generally composed of the space free curved surface, and the production of it is a process of multistep and multiprocess due to its complex shape, small relative thickness, large size, and high surface quality requirements. The forming processes of automobile covering parts generally include drawing, trimming, stamping, shaping, flanging, and other processes, of which the stamping process is the most important step of all, and the results of the stamping are directly related to the product quality and subsequent process. In order to get qualified stamping automotive panel to meet customer requirements, multiple modifying and repeated debugging are required for stamping die, stamping process, technological parameters, and so forth, while “trial and error” method is used in the traditional way for repeated testing, which not only require a longer production cycle, but also require high production costs.
With the continuous development of numerical analysis technology, plastic forming theory, stamping technology, and computer technology, CAE (computer aided engineering) technology is gradually applied in stamping process and die design process for automotive panels [3, 4]; CAE related software are more and more widely used in panel die debugging process. The formability and possible defects of parts can be analyzed through the simulation, which provide the reliable basis for die debugging, so that it can reduce the die manufacturing problems, shorten the production cycle, save production costs, and offer more economic benefits for enterprises [5, 6].
In the stamping forming process of automotive panels, many process parameters have important influence on the quality of forming parts, and many scholars have done many studies in the optimization of process parameters [7–10]. In order to study the effects of several process parameters which have great influence on the forming quality of parts, which can be optimized, this paper adopts automotive lower floor board as the research object. Under the draw bead distribution, the shape and parameters remain unchanged; use the orthogonal experiment and ETA (engineering technology associates)/Dynaform simulation software to analyze BHF (blank holder force), die corner radius, die clearance, and friction coefficient to study impact level of these four parameters with the results of maximum thinning rate and maximum thickening rate and optimize these parameters.
2. Establishment of Simulation Analysis Model
The simulation process of stamping forming mainly includes the establishment of 3D model and finite element models, defining the forming tools, setting the process parameters, and analyzing the results. The flowchart of the process simulation is shown in Figure 1.
2.1. Establishment of 3D Model of Part Automotive Lower Floor Board
According to the requirements of the analysis, draw threedimensional model of automotive lower floor board parts by using Pro/E firstly. Analysis of the parts structure shows that the parts have large transverse size, small longitudinal size, large drawing depth, and complex shapes, so final shape should be finished; go through drawing → trimming → punching → shaping → flanging → lateral punching → trimming and other processes. Depending on the actual production experience, the cracking usually occurs at the lower portion of both ends of the parts; the flat ends and the mid plane may produce insufficient forming. On the basis of the defects the parts may arise; necessary process supplement and binder surface shape design are required to reduce defects in the process of parts stamping. Drawing die model is shown in Figure 2; the part is defined as die and used as the forming tool during the simulation, and the model is saved as .igs format for subsequent finite element modeling.
2.2. Establishment of Finite Element Models
Dynaform software has a total of 16 element types; currently the most commonly used are HuhgesLiu (HL) shell element and BelytschkoTsay (BT) shell element, and these two elements are constructed based on the Mindlin plate theory with the advantages of fast calculation speed and high simulation accuracy [11]. The default element type is BelytschkoTsay shell element in Dynaform software, which is the most efficient element in explicit algorithm. In this paper, the die and sheet metal are meshed by using BT shell element. Comprehensively considering simulation accuracy and computation time, sheet mesh size 12 mm and adaptive grade 4 are accepted. And the finite element model of automotive lower floor board is shown in Figure 3. In order to ensure the overall forming quality of parts, draw bead distribution, shape, and parameters remain unchanged. The sheet metal use stamping steel WSSM1A365A12; the mechanical properties of the materials are given in Table 1, and the plate thickness is 0.7 mm. The plate is located on binder; the constitutive model of sheet material applied in this paper is the three parameters BarlatLian thick anisotropic yield criterion [12]; the criterion is where , , is equivalent to yield stress, is sheet stress along the rolling direction, is sheet stress vertical to the rolling direction, , , , and are the coefficient of material plastic anisotropy, and is BarlatLian index related to the crystal structure.

3. Orthogonal Experiment Design
Orthogonal experiment is a design method to study multiple factors level which is often used for the performance prediction of stamping forming of automotive panels [13–16]. During the process of automotive lower floor board stamping, the main process parameters affecting the forming performance include BHF, die corner radius, die clearance, and friction coefficient. For the above process parameters, orthogonal experiment factors were designed for the BHF (factor ), die corner radius (factor ), die clearance (factor ), and friction coefficient (factor ). The maximum thinning rate and maximum thickening rate are evaluation index, and the maximum thinning rate is as follows: where is the thickness of the sheet before drawing and is the thickness of unit of sheet metal after drawing.
Change () in formula (2) into () for maximum thickness ratio, and maximum thinning rate and maximum thickening rate of products should be as small as possible during the optimization in order to ensure product strength.
BHF is the force on the blank holder in drawing process; it can be used to solve the wrinkling problem in the process of drawing, and it can be calculated by the following formula: where is the area of the pressure side (mm^{2}); is for unit BHF (MPa), and its value is determined by experiment; unit BHF of different materials is selected according to Table 5.3 in [17]. From formula (3), the value of BHF was calculated to be about 500 kN, and the three levels of experiments were chosen for 400 kN, 500 kN, and 600 kN.
There is virtually no effect of the punch speed on the forming performance [18]; we took the punch speed as 3000 mm/s during all the simulation. The die corner radius can affect the stamping process and we took 3 mm, 6 mm, and 9 mm for the experiment. Die clearances were , , and with as sheet metal thickness. We set friction coefficient as 0.10, 0.12, and 0.14 for the experiment. The above factors were listed in Table 2; combination schemes of orthogonal experiment process parameters were designed by using orthogonal array L_{9} (3^{4}) as shown in Table 3.


4. Results and Analysis of Orthogonal Experiment
Due to the large number of experiment data, not all simulation results can be displayed; only the simulation results of process parameters scheme through the orthogonal experimental design corresponding to maximum thinning rate and maximum thickening rate were counted in Table 4.

From the simulation results in Table 4, it can be seen that the maximum thinning rate of most experiment parts is below 30%. In general, the reqkuirement of maximum thinning rate is not more than 30%, while some parts of the automotive lower floor board thinning rate has exceeded the maximum value, so that the part has a very large risk of cracking; maximum thickening rate for all parts of the experiment program is between 4% and 7%; there are small gaps among the results of each experiment scheme; the maximum thickening region is mainly in binder surface, which has a large risk of wrinkling. There are two common methods for solving wrinkling of the binder surface; one is to increase the BHF, and one is to reduce the intensity of drawbeads, and the degree of solving wrinkling should be in accordance with customer requirements. The analysis results of the maximum thinning rate and maximum thickening rate are shown in Table 5, with , , and , respectively, corresponding to the sum of 1, 2, and 3 level of () factors, , , and , respectively, corresponding to the average value of 1, 2, and 3 level of () factors, and stands for range corresponding to each factor. The range in Table 5 represents the influence degree of various factors on the results. It can be seen that the range values of each factor are not equal, which indicates that the factors have different impacts on forming performance of automotive lower floor board. The relationship between the maximum thinning rate and experimental factors is shown in Figure 4, and the relationship between the maximum thickening rate and experimental factors is shown in Figure 5.

We can see clearly from the relationship shown in Figures 4 and 5 that the change trend of the maximum thinning rate and the maximum thickening is almost in the opposite condition with the experimental factors. The maximum thinning rate of automotive lower floor board increases along with the increases of BHF; the maximum thickening rate decreases along with the increases of BHF. So the large BHF can avoid the wrinkle and also may lead to the risk of cracking. The influence of friction coefficient is the same as BHF. The maximum thinning rate decreases along with the increases of die corner radius, and the maximum thickening rate increases first and decreases then along with the increases of die corner radius. The influence of die clearance is minor on these two evaluation indexes.
5. Optimization Schemes and Simulation Results
Due to some parts of automotive lower floor board higher thinning exists during stamping process; these areas would have a greater risk of cracking, and some parts of larger thickening would have a greater risk of wrinkling, double indexes optimization for process parameters were conducted by integrating these two reasons respectively. Single index of maximum thinning rate and maximum thickening rate corresponding to process parameters optimization results are given in Table 5. For the two factors, the maximum thinning rate and maximum thickening rate are the smaller the better, and the smaller the value is, the smaller the risk of cracking and wrinkling on the corresponding parts becomes; the defects of parts also reduce.
For the maximum thinning rate, die corner radius is the biggest influencing factor which determines the structure of die, friction coefficient is the secondary influencing factor which can be changed by adding the lubricant and changing the surface roughness, BHF can be modified by adjusting the parameters of mechanical presses, and die clearance can be adjusted by die repair. According to optimization scheme on the maximum thinning rate, the best optimization scheme is A1B3C3D1; with these parameters we take the simulation with ETA/Dynaform and the results of thickness distribution and forming limit diagram of A1B3C3D1 optimization scheme are, respectively, shown in Figures 6 and 7. After the optimization of process parameters, the maximum thinning rate of automotive lower floor board is 28.53%, the maximum thickening rate is 6.34%, the two rates are controlled in the required range, and the overall forming quality of parts is better without obvious defects.
For the maximum thickening rate, the biggest influencing factor is BHF and is followed by friction coefficient; die corner radius and die clearance have little effect on the maximum rate of thickening. The best optimization scheme for maximum thickening rate is A3B1C1D3, and the simulation results of thickness distribution and forming limit diagram are, respectively, shown in Figures 8 and 9. In this case, the maximum thinning rate of automotive lower floor board is 40.10%; the maximum thickening rate is 4.71%. Obviously, the maximum thinning rate exceeds the forming limit and the automotive lower floor board will be cracked.
According to the orthogonal experiment and simulation by Dynaform, we chose the friction coefficient with 0.10, BHF with 400 kN, die clearance with 0.84 mm, and die corner radius with 9 mm to take the actual production of automotive lower floor board for experiment. The stamping forming part is shown in Figure 10; we can find the lubricant on the surface of automotive lower floor board, which is used to reduce the friction coefficient of the parts to avoid the risk of cracking. The wrinkle only occurs in the outer plate where it will be cut off in the next process. The final product part is shown in Figure 11. We can find that the automotive lower floor board is in high quality without any defects.
6. Conclusions
In order to study the influence of several process parameters which have great influence on the forming quality on the quality of automotive panels, this paper adopted automotive lower floor board as the research object, selected the BHF, die corner radius, die clearance, and friction coefficient with maximum thinning rate and maximum thickening rate as evaluation indexes, used orthogonal experiment method for simulation analysis on the effects of these four parameters, and carried on optimization. It can be obtained from the experiment that the impact of the four process parameters on the maximum thinning rate from strong to weak is die clearance, friction coefficient, the BHF, and die clearance, while on the maximum thickening, the greatest impact is BHF, followed by the friction coefficient and die clearance and die clearance has little effect. The maximum thinning rate and maximum thickening rate could be effectively controlled through orthogonal experiment optimization, and the high quality forming parts can be obtained without obvious defects. This paper can also verify the feasibility of ETA/Dynaform simulation software in the mold manufacturing and debugging stage for actual production guidance.
Conflict of Interests
The authors declare that they do not have any commercial or associative interest that represents a conflict of interests in connection with the work submitted.
Acknowledgments
This work was supported by the Innovation Foundation for Graduate Students of Shandong University of Science & Technology (Grant nos. YC140105 and YC140205). The authors express sincere thanks to the reviewers for their helpful comments and suggestions for improving this paper.
References
 M. Cavazzuti, A. Baldini, E. Bertocchi, D. Costi, E. Torricelli, and P. Moruzzi, “High performance automotive chassis design: a topology optimization based approach,” Structural and Multidisciplinary Optimization, vol. 44, no. 1, pp. 45–56, 2011. View at: Publisher Site  Google Scholar
 F. Bi, S. J. Hu, and B. Song, “Application of nonGaussian closure method in automotive chassis nonlinear vibration analysis,” Forschung im Ingenieurwesen, vol. 75, no. 3, pp. 153–163, 2011. View at: Publisher Site  Google Scholar
 N.M. Wang and B. Budiansky, “Analysis of sheet metal stamping by a finite element method,” Journal of Applied Mechanics, vol. 45, no. 1, pp. 73–82, 1978. View at: Publisher Site  Google Scholar
 N. Wang and S. Tang, “Analysis of bending effects in sheet forming operation,” International Journal for Numerical Methods in Engineering, vol. 25, no. 1, pp. 253–267, 1988. View at: Google Scholar
 M. Kawka, T. Kakita, and A. Makinouchi, “Simulation of multistep sheet metal forming processes by a static explicit FEM code,” Journal of Materials Processing Technology, vol. 8081, pp. 54–59, 1998. View at: Google Scholar
 J. Mackerle, “Finite element analyses and simulations of sheet metal forming processes,” Engineering Computations, vol. 21, no. 8, pp. 891–940, 2004. View at: Publisher Site  Google Scholar
 Y. Wen, W. Zhong, and Y. Liu, “Optimization design of process parameters in sheet metal drawing,” Journal of Plasticity Engineering, vol. 20, no. 3, pp. 31–36, 2013. View at: Publisher Site  Google Scholar
 J. Pan, Y. Zhong, and C. Yuan, “Process parameters optimization for sheet metal forming during drawing with a multiobjective genetic algorithm,” Journal of Tsinghua University, vol. 47, no. 8, pp. 1267–1269, 2007. View at: Google Scholar
 E. Gao, H. Li, H. Kou, H. Chang, J. Li, and L. Zhou, “Influences of material parameters on deep drawing of thinwalled hemispheric surface part,” Transactions of Nonferrous Metals Society of China, vol. 19, no. 2, pp. 433–437, 2009. View at: Publisher Site  Google Scholar
 W. Liu and Y. Yang, “Study on process optimization and objective functions of sheet metal forming based on FEA,” Material Science and Technology, vol. 14, no. 2, pp. 159–161, 2006 (Chinese). View at: Google Scholar
 W. Chen, Sheet Metal Forming CAE Analysis Tutorial, China Machine Press, Beijing, China, 2005, (Chinese).
 F. Barlat and J. Lian, “Plastic behavior and stretchability of sheet metals. Part I: a yield function for orthotropic sheets under plane stress conditions,” International Journal of Plasticity, vol. 5, no. 1, pp. 51–66, 1989. View at: Publisher Site  Google Scholar
 D. Ko, S. Cha, S. Lee, C. Lee, and B. Kim, “Application of a feasible formability diagram for the effective design in stamping processes of automotive panels,” Materials and Design, vol. 31, no. 3, pp. 1262–1275, 2010. View at: Publisher Site  Google Scholar
 M. Firat, O. H. Mete, U. Kocabicak, and M. Ozsoy, “Stamping process design using FEA in conjunction with orthogonal regression,” Finite Elements in Analysis and Design, vol. 46, no. 11, pp. 992–1000, 2010. View at: Publisher Site  Google Scholar
 K. Park and Y. Kim, “The effect of material and process variables on the stamping formability of sheet materials,” Journal of Materials Processing Technology, vol. 51, no. 1–4, pp. 64–78, 1995. View at: Publisher Site  Google Scholar
 L. Chen, W. Chen, J. D. Li, S. N. Heng, and J. Wu, “Optimal design about parameters of cooling pipes in hot stamping die,” Advanced Materials Research, vol. 988, pp. 263–267, 2014. View at: Google Scholar
 S. Wu, Sheet Metal Forming Technology and Mold Design, Northwestern Polytechnical University Press, Xi’an, China, 2002 (Chinese).
 K. Chung, K. Hyunki, and C. Lee, “Forming limit criterion for ductile anisotropic sheets as a material property and its deformation path insensitivity, part i: deformation path insensitive formula based on theoretical models,” International Journal of Plasticity, vol. 58, no. 7, pp. 3–34, 2014. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Guoying Ma and Binbing Huang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.