Research Article  Open Access
Optimization of Bending Process Parameters for Seamless Tubes Using Taguchi Method and Finite Element Method
Abstract
The threedimensional tube (or pipe) is manufactured by CNC tube bending machine. The key techniques are determined by tube diameter, wall thickness, material, and bending radius. The obtained technique through experience and the trial and error method is unreliable. Finite element method (FEM) simulation for the tube bending process before production can avoid wasting manpower and raw materials. The computeraided engineering (CAE) software ABAQUS 6.12 is applied to simulate bending characteristics and to explore the maximum stress and strain conditions. The Taguchi method is used to find the optimal parameters of bending. The confirmation experiment is performed according to optimal parameters. Results indicate that the strain error between CAE simulation and bending experiments is within 6.39%.
1. Introduction
In the mechanical manufacturing industry, tube bending is very commonly used in processing and its products are widely applied in major industries, such as power plants, aviation, petrochemical industry, and the boiler industry.
The bending process for thin wall tubes is often a technical issue. The basic parameters of the copper elbow bending process include the bending moment, the type and degree of deformation of cross section, and change of the wall thickness. The axial upsetting force on the above mentioned parameters was applied to improve problems associated with the changes of wall thickness [1].
Tubes subjected to a point load result in local oval forming which influences tube strength. Oval forming and the flattening effect of tubes experienced during bending was first researched by Brazier and therefore was called the Brazier effect [2]. An increase in internal pressure, together with a point load condition, results in a change to sectional properties of the tube, thus leading to bending of the tube and ultimately structural failure.
Shaw and Kyriakides [3, 4] experimentally tested the response and stability of elastoplastic pipes under combined bending and external pressure. They designed a small lightweight instrument, which can be used to measure the change in the diameter of the tube at a point in the plane of bending. The same instrument is used to experimentally measure the ovalization of inelastic tubes under cyclic bending. They demonstrated that reverse bending and subsequent repeated cyclic bending cause a gradual growth of the ovalization of the tube’s cross section. They further extended their experimental program to investigate the stability of thinwalled tubes under cyclic bending. The results indicated that, under curvaturesymmetric loading, the tube progressively ovalizes to a critical value at which point it buckles.
Vollertsen et al. [5] gave a systematic overview about the existing processes for the bending of extrusions in both small and large scale production. Some of the important features of working accuracy in bending extrusions were reviewed. This included the accuracy of the curvature and deviations in the cross section. Methods of open and closedloop control based on process modeling for increased accuracy were shown as well as practical methods to avoid cross section deformation.
For pipe bends subjected to outofplane bending, the addition of torsional effects makes it more complex to estimate the plastic collapse moment. The collapse moment is reached first at a section around the middle of the bend part. Mourad and Younan [6, 7] performed an analysis to evaluate the behavior of a pipe bend subjected to outofplane bending and internal pressure. They showed that the loaded end of the bend was the most severely strained cross section and found that considerable plastic deformation occurred before instability was reached, especially in the presence of applied pressure.
Cross section ovalization is one of the flaws in the thinwalled tube NC bending process. A unique quality estimator of the cross section geometry of bent tubes has been developed, which is able to synthetically represent the main functional and aesthetical properties of tubes. The deformation occurring in the cross section of the tube has been studied through a set of FEM simulations, run with an explicit code with shell elements. The model has been validated by comparing the experimental and numerical results of specific cases in terms of geometry of the cross section in a 90° bend. A proposed indicator is able to reliably estimate the difficulty of a bending process with respect to a tendency towards ovalization and collapse of the cross section [8].
Much research, using FEM and experiments, has been carried out on the thinwalled tube bending processes. The stress and strain distribution, tube wall thinning, cross section ovalization, and other defects in the NC bending have been studied [9]. The finite element method (FEM) is a useful tool for optimizing the process layout in metal forming. Numerical simulation helps designers to estimate metal formability and improve process layout.
In this study, the computeraided engineering (CAE) software ABAQUS applied finite element method (FEM) to simulate bending characteristics and to explore the maximum stress and strain conditions. The Taguchi method was used to find the optimal parameters of bending forming. The confirmation experiment using strain gauge was performed according to optimal parameters. The results indicate that the CAE precisely simulated the bending experiment.
2. Elastoplastic Analysis
If an elastoplastic deformable body is in dynamic equilibrium state, the following force balance equation (1) and boundary conditions (2) are established for every material point which is covered in Ω of the body: where = material density, = displacement vector, = Cauchy stress tensor = body force vector, = drag force vector on the body surface, = spring constant, = known displacement vector on the displacement boundary , = known load vector on the load boundary , and = mixed boundary, the sum of and .
Principle of virtual work, finite element, and numerical integration are applied to derive the following nonlinear differential equation:where = total mass matrix, = total damping matrix, = total stiffness matrix, and = total load vectors.
3. Design of Experiments
The Taguchi method was created by Genichi Taguchi in 1949. It is an engineering methodology for obtaining the relationship between products and the manufacturing process.
The Taguchi method is an important tool for robust design, combining experimental design theory and the concept of quality loss function. It offers a simple and systematic approach to optimize performance design, quality, and cost [10–12].
Taguchi’s approach is based on statistical design of experiments, which fulfills the requirements to solve engineering problems and enhance process optimization. Taguchi suggests analyzing variation using an appropriately chosen signaltonoise ratio (). These ratios are derived from the quadratic loss function. The Taguchi analysis of the ratio involves three kinds of quality characteristic, including thenominalthebetter, thesmallerthebetter, and thelargerthebetter. To obtain optimal bending performance, thesmallerthebetter quality characteristic for minimizing maximum stress and strain must be taken. It can be expressed aswhere is the observed data. Regardless of category, is the number of experiments in a trial. A greater ratio corresponds to a better performance. The level of a factor with the highest signaltonoise ratio is the optimum level [13].
In this study, the finite element analysis combined with Taguchi methods L_{18} (2^{1} × 3^{7}) orthogonal array was used to simulate the bending process of threedimensional seamless pipe. Table 1 represents experimental factors and their levels. Table 2 is the orthogonal array. CAE software ABAQUS was applied to simulate bending according to orthogonal array. The optimal bending process parameters were obtained via ratio analysis. Finally, the confirmation experiment was performed to verify the consistency between the simulation and the actual bending process.


In practical threedimensional bending process, CNC tube bending process was illustrated in Figures 1 and 2. axis was the feed direction of the tube; axis was the rotation direction of the tube along its axis, and axis was bending direction of the tube which was bent by wiper die. Three form blocks were fixed and axis of seamless tube were fixed. In accordance with the design parameters of the Taguchi method wiper die was rotated along bending axis of the tube while tube rotated and fed itself to complete , , and three strokes.
The material used in this experiment was the seamless SUS304 stainless steel tube with corrosion resistance, high strength, electrical conductivity, and good mechanical properties (Table 3); its maximum stress was 626 MPa and maximum strain was 0.42. This experimental tube had a length of 800 mm, an inner diameter of 34 mm, and a wall thickness of 1.5 mm.

4. Results and Discussions
The maximum stresses and strains were obtained by ABAQUS simulation analysis. Eighteen data sets are listed in Table 4. Set numbers 2 and 12 in Table 4 show that maximum stresses were very close to the tensile strength of material SUS304 (626 MPa). Set numbers 13 and 18 of plastic strain also exceeded the maximum strain (0.42). Thus, these four sets of bending pipe will be broken.

The experimental results were then transformed into a signaltonoise () ratio. Taguchi recommends the use of the ratio to measure quality characteristics. A greater ratio corresponds to improved quality. In this study, thelowerthebetter quality characteristic (5) was used to calculate the ratio for maximum stress and strain and the results are listed in Table 5.

Table 6 represents the ratio at each level of all control factors for strain. The effect of the control factor was determined by its level of difference values. A larger control factor level difference resulted in a greater effect on strain in the bending process. Figure 3 shows the response graph for strain.

Table 7 expresses the quality characteristic (strain) which was obtained from actual experiment. The factor response table (Table 8) and graph (Figure 4) of the quality characteristic can be made according to the experimental data. Tables 6 and 8 indicated that factor D was the second most effective factor for strain and the most effective factor for quality characteristic, because it had the highest level difference among all control factors.


The optimal combination of bending parameters for minimizing strain was A1B3C2D1E2F3G1H3. The optimal combination of bending parameters for quality characteristics was A2B3C2D1E2F3G1H3. The level of factor A (feed stroke) was the only difference, and A2 was selected for the sake of productivity. The optimum parameters were listed in Table 9.

The optimal parameters were reused through CAE analysis to obtain the maximum stress of 503.5 MPa and maximum plastic strain of 0.327. A comparison was made between these two values and the average values of 18 sets of data (Table 4). The results indicated that the maximum stress of 503.5 Mpa was very close to the average stress of 503.1 Mpa obtained from the Taguchi method of 18 sets of data, and maximum strain of 0.327 was lower than the average strain of 0.353 obtained from the Taguchi method of 18 sets of data, and it can be made sure that the Taguchi method selected optimal parameters are correct.
A confirmation test was done after the optimal combination of process parameters was determined. The measure point of strain obtained by finite element analysis is near the bending radius of the tube; the strain value is 0.1342, as shown in Figure 5. At the same position, the strain obtained by actual tube bending experiments using strain gauges is 0.1428. The results indicate that the error between confirmation test and actual experiment is 6.39% (Table 10).

The completed workpiece using optimal parameters is shown in Figure 6. The surface around the curved corner of the optimum workpiece is obviously smoother than the original workpiece.
5. Conclusion
The CNC bending parameters were selected as Taguchi experimental factors rather than bending radius. This is because the actual bending process had to consider operation convenience. The results show that the selected optimum parameters can actually apply to bending operation without further conversion. The experimental results may be summarized as follows.(1)According to response diagram, factor D (rotation axis stroke 2) is the most effective process factor. The level of this factor should be smaller to avoid product failure.(2)The optimal bending parameters for minimizing maximum strain are A2, B3, C2, D1, E2, F3, G1, and H3. They are as follows: feed strokes 1 and 2 = 140 mm; feed stroke 3 = 120 mm; rotation stroke 1 = −250°; rotation stroke 2 = −225°; rotation stroke 3 = −80°; bending stroke 1= 46.5°; bending stroke 2 = 46.5°; and bending stroke 3 = 49°.(3)The optimal parameters were utilized to execute a confirmation experiment, the results of which indicate that the strain error between experiment and FEM simulation is 6.39%.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
References
 F. Stachowicz, “Bending with upsetting of copper tube elbows,” Journal of Materials Processing Technology, vol. 100, no. 1, pp. 236–240, 2000. View at: Publisher Site  Google Scholar
 L. G. Brazier, “On the flexure of thin cylindrical shells and other ‘Thin’ sections,” Proceedings of the Royal Society of London A, vol. 116, pp. 104–114, 1927. View at: Publisher Site  Google Scholar
 P. K. Shaw and S. Kyriakides, “Inelastic analysis of thinwalled tubes under cyclic bending,” International Journal of Solids and Structures, vol. 21, no. 11, pp. 1073–1100, 1985. View at: Publisher Site  Google Scholar
 S. Kyriakides and P. K. Shaw, “Inelastic buckling of tubes under cyclic bending,” Journal of Pressure Vessel Technology, vol. 109, no. 2, pp. 169–178, 1987. View at: Publisher Site  Google Scholar
 F. Vollertsen, A. Sprenger, J. Kraus, and H. Arnet, “Extrusion, channel, and profile bending: a review,” Journal of Materials Processing Technology, vol. 87, no. 1–3, pp. 1–27, 1999. View at: Publisher Site  Google Scholar
 H. M. Mourad and M. Y. A. Younan, “The effect of modeling parameters on the predicted limit loads for pipe bends subjected to outofplane moment loading and internal pressure,” Transactions of the ASME, Journal of Pressure Vessel Technology, vol. 122, no. 4, pp. 450–456, 2000. View at: Publisher Site  Google Scholar
 H. M. Mourad and M. Y. A. Younan, “Nonlinear analysis of pipe bends subjected to outofplane moment loading and internal pressure,” Journal of Pressure Vessel Technology, vol. 123, no. 2, pp. 253–258, 2001. View at: Publisher Site  Google Scholar
 A. Mentella and M. Strano, “Rotary draw bending of small diameter copper tubes: predicting the quality of the crosssection,” Proceedings of the Institution of Mechanical Engineers Part B: Journal of Engineering Manufacture, vol. 226, no. 2, pp. 267–278, 2012. View at: Publisher Site  Google Scholar
 F. Paulsen and T. Welo, “Application of numerical simulation in the bending of aluminiumalloy profiles,” Journal of Materials Processing Technology, vol. 58, no. 23, pp. 274–285, 1996. View at: Publisher Site  Google Scholar
 E. Bagci and Ş. Aykut, “A study of Taguchi optimization method for identifying optimum surface roughness in CNC face milling of cobaltbased alloy (stellite 6),” The International Journal of Advanced Manufacturing Technology, vol. 29, no. 910, pp. 940–947, 2006. View at: Publisher Site  Google Scholar
 N. Alagumurthi, K. Palaniradja, and V. Soundararajan, “Optimization of grinding process through Design of Experiment (DOE)—a comparative study,” Materials and Manufacturing Processes, vol. 21, no. 1, pp. 19–21, 2006. View at: Publisher Site  Google Scholar
 M. S. Phadke, Quality Engineering Using Robust Design, Prentice Hall, Englewood Cliffs, NJ, USA, 1989.
 H. H. Lee, Taguchi Method, Quality Characteristics Concept and Practices, GaoLi Books, Taipei, Taiwan, 2005.
Copyright
Copyright © 2015 JuiChang Lin and Kingsun Lee. 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.