Enhancing Multistage Deep-Drawing and Ironing Manufacturing Processes of Axisymmetric Components : Analysis and Experimentation

An optimization technique for combined processes of deep-drawing and ironing has been created in order to improve the total process time and cost in manufacturing procedures of axisymmetric components. The initial solution is optimized by means of an algorithm that minimizes the total time of the global process, based on relationship between lengths, diameters, and velocities of each stage of a multistage process and subject to constraints related to the drawing ratio. The enhanced solution offers a significant reduction in time and cost of the global process. The final results, applied to three cases, are compared with experimental results, showing the accuracy of the complete solution.


Introduction
The industry of metallic components manufacturing requires developments to be more efficient, in particular in the deepdrawing procedures, where it is important to decrease the process times and costs as in mass production.Thus, it is necessary to devise specific algorithms that will satisfy these demands.These algorithms should be based on technological and scientific basis that will provide solutions that are ready for transfering to the industries.The deep-drawing process has been analysed with this objective in mind due to its convenient nature as a global model which includes all stages of the process, namely, drawing, redrawing, and ironing.
The majority of literature contributions are focused on the study of properties of process, in particular, the prediction of the limiting drawing ratio (LDR) [1][2][3], the blank design using different methodologies, such as parametric NURBS surfaces [4], upper bound method [5], or artificial neural network [6,7], the effect of die radius on the blank holder force and drawing ratio [8], the predicted thickness distribution of the deep drawn circular cup of stainless steel [9], the improvement of drawability by means of technological parameters [10,11] or the formability with different thickness [12].However, some efforts have been realised about the parts design [13] or generation of algorithms, mainly related to the process planning; Ramana and Rao [14] developed a framework based on knowledge related to design-process planning integration for sheet metal components, although there is no evidence of its application.Also, Vosniakos et al. [15] devised an intelligent system to process design of sheet parts.As can be seen, the researches of deepdrawing processes are not focused on the reduction of time, despite frequently being used on mass production due to the characteristics of the parts.
This paper presents a model that provides a comprehensive analysis of those phenomena occurring in the multistage processes of axisymmetric geometry and applied to the manufacturing of this type of components.The scientific development stems from the work done by Leu [1] and Sonis et al. [2] which provides LDR solutions based on normal International Journal of Manufacturing Engineering anisotropy value, strain hardening exponent, and others, applied to the drawing and redrawing stages.The model [16,17] has a scientific foundation based on the literature covering the plastic deformation processes and, in particular, drawing [18], redrawing [19], and ironing.This study focuses on analyzing geometries of axisymmetric components, manufactured by a multistage deep-drawing process.
The drawing processes from a blank [6,8] have been researched and their governing equations are well known [3,5].However, the redrawing and ironing phenomena have received relatively less attention [2,20].This work contributes to the definition of a global process of drawing, redrawing, and ironing combined process, together with a common focus on the drawing process in order to obtain a global solution.
The model is used to carry out a quantitative and integrated analysis in a multistage deep-drawing process based on a scientific solution and real measurements.Also, the model permits the modification for some process variables to predict their influence on the process [9].In addition, it is based on the definition of limiting conditions to guarantee the stability of the process.The simultaneous accomplishment of limiting situations of each process (drawing, redrawing, and ironing) allows for fixing a boundary for the values of each stage which gives the initial solution.The model permits the optimisation of the initial solution from several points of view: total process time, manufacturing cost, among others, such as punch and die wear.This optimisation is based on the resolution of an algorithm by means of recursive functions, which explores all the possibilities of the process and selects the most adequate one.The algorithm is supported by a software tool that provides the optimized solution [21].The optimized solution is compared with the experimental results.
The remainder of the paper is organised as follows.In Section 2, the methodology used is formulated and applied to some industrial cases.In Section 3, the initial solution found from the scientific model is presented.Section 4 defines the optimization process.In Section 5, the optimization algorithm for solving the multistage process is proposed.In Sections 6 and 7, the simulated and experimental results are compared in order to show the reliability of the complete solution found by the proposed algorithm.

Methodology
The methodology used in the definition and resolution of the algorithm is based on a model that departs from an initial solution, defined by the technological constraints that characterize the multistage deep-drawing processes (see Section 3).The use of these processes in mass production will require the minimization of the total process time and the quantity of raw material, thus acting as constraints.
The algorithm resolution has been carried out by the software called deep-drawing tool (DDT) developed by the author.This tool allows the user to perform the total multistage process; it allows the evaluation of the most important variables of the process and the implementation of the optimization algorithm which is the objective of this paper.The algorithm results are contrasted experimentally by means of the resolution of several industrial cases.The applied product type is an axisymmetric component produced in brass by multistage deep-drawing and ironing processes.The use of this material is justified due to its good properties suitable for these multistage processes.Table 1 shows the characteristics of brass UNS C26000 used on the resolution of the industrial cases presented in Section 5.
The solutions provided by this algorithm have been tested against three industrial cases and geometries as shown in Table 2.
The resolution of these types of cases requires the combination of different processes of drawing, redrawing, and ironing.Figure 1 shows the drawing and redrawing/ironing press used on the manufacturing processes of the industrial cases.Figure 2 shows some of the manufacturing stages for Case A.

Model Approach
Figure 3 shows the flow diagram of the model providing the initial solution: model phases, input variables required to define the "new project, " and the flow information between different stages.The model develops the required solution from the input data.These input data correspond to the dimensions and material of the final piece.The required input data are the external diameter (  ), the length (  ), bottom thickness (  ), wall thickness (  ), and the material type.The model permits the selection of the material type and its mechanical characteristics needed to form the part: density,  (kg/m 3 ); ultimate tensile strength,   (MPa); yield tensile strength,   (MPa); material rigid-plastic constant,  (MPa); strain hardening exponent, ; and normal anisotropy value, .From the input data, the dimensions of the initial part or blank are calculated based on the incompressibility condition of the plastic deformation.This condition considers that the volume of the piece is unchanged throughout the deformation process [18].A detailed formulation about the previous model that determines the initial solution can be found in  Ramirez et al. [22].The dimensions of this blank are the starting point to develop the drawing, redrawing, and ironing processes until the required final dimensions are achieved.In this way, data from the blank becomes the input for the next step: drawing.

Drawing.
The model presents a calculation procedure for each process: drawing, redrawing, and ironing.The determination of the initial solution is treated independently for each of the threads that may occur.It is important to note that, depending on the type of piece, in some cases, its geometry is not necessary to consider redrawing or ironing stages.For determining the initial solution of the drawing, a hypothesis that would later be amended in the adjustment of the technological process is used, and it is the size of the diameter of the punch corresponding to that stage.The LDR is a measure of deep draw ability which defines the largest blank that can be drawn without tearing.It is the ratio between the maximum blank diameter and punch diameter.The design of this punch will depend on the LDR and the final dimensions of the piece.It is also possible that this diameter is changed later in the adjustment of the technological process in order to allow a clearance between the punch and the workpiece for subsequent stages of redrawing and ironing [20].The solution for the initial stage of drawing is calculated by considering two limiting conditions as described below.The model selects the drawing diameter or the die, , by using the largest diameter ( 1,1 and  1,2 ) obtained from the two drawing conditions.Once the diameter is known, the model determines all the dimensions needed to define the drawing stage.Verification of the drawing stage dimension is conducted using two experimental conditions based on empirical studies and collected field data.The corresponding flow diagram is shown in Figure 4.
The first limiting drawing condition is related to the maximum force ( ,max ) executed by the punch on the workpiece during the drawing process, which must be less than the cracking load of the material ( cr ), according to (1).These forces establish the range in which the force can be obtained using from the efficiency coefficient ( def ), ultimate tensile strength, blank thickness ( 0 ), and the initial ( 0 ) and final ( 1 ) diameters (note that the mean wall diameter is  ,1 =  1 +  0 ), and they are based on the experimental expression from Siebel and Beisswänger [23].This condition requires the following: cr =  ⋅  ,1 ⋅  0 ⋅   . ( The second drawing limiting condition is calculated by using the expression developed by Leu [1], which showed good agreement between theoretical and experimental results [2,3].Considering the condition of constant volume throughout the process of plastic deformation, LDR can be defined as follows: ( This expression estimates the LDR as a function of normal anisotropy, the strain hardening exponent, and efficiency ().
In this manner, upper and lower bounds are established, considering the materials proprieties and the drawing capacity.The two limiting conditions use new data from the blank geometry and process efficiency ( def and ).

Redrawing.
Once the dimensions of the drawing stage are obtained, the model provides these data as input values for the next phase of the process: redrawing.The goal of this step is to obtain the final dimensions of the piece needed to perform the next step: ironing.The solution for the initial stage of redrawing is calculated from the consideration of three limiting conditions as described below.The model selects the diameter of the die from the largest of the solutions obtained in the three redrawing constraints ( 1⋅,1 ,  1⋅,2 , and  1⋅,3 ).Once the diameter is known, the model determines all dimensions needed to define this stage.Figure 5     , intemediate stages The first redrawing limiting condition is related to the redrawing load of the punch during the process of redrawing ( re,max ), which must be less than the cracking load of the material ( cr ), according to (3) and (4), respectively.Consider The first condition is based on Siebel and Beisswänger [23], where  is the friction coefficient,  is the die half-angle,   is the punch diameter, and   is the die radius.
The second redrawing limiting condition assumes the rigid-plastic condition of the material.Assuming that the material behaves according to a parabolic law that approximates the potential behaviour of metallic materials coldmanufactured, then it is possible to determine an expression relating the external diameter of the part of a generic step of the process of redrawing,  1⋅ , to the diameter of the previous stage.Assuming the material is annealed (  = 0) and considering the total deformation between the initial and final states, the condition takes the following form: (5) The third redrawing limiting condition is referred to as the limiting drawing ratio.In applying the third limiting condition of drawing, the model applies the formulation defined by Sonis et al. [2] about the LDR study in the operations of redrawing.The model considers the effects of the normal anisotropy of the material, friction coefficient, coefficient of strain hardening, and the radius of the input die (  ).The LDR is used in this model as a variable to determine the required number of redrawing steps and size of the stages.It is assumed that the material is rigid-plastic [2].Moreover, assuming that the material is rotationally symmetric, the same properties are based on the existence of normal anisotropy and planar isotropy.The Sonis model [2] is based on the tension that is created in the area of the radius of the die redrawing causing instability in the plastic wall of the cup, which is equal to the radial tension in the drawing area of the flange, due to the continuity of tension throughout the piece.Based on the Sonis model, the expression for LDR  can be written as +   (−1) /LDR (−1) + √ (  +   (−1) ] ] ] −  5 ln ((   +  (−1) International Journal of Manufacturing Engineering where  1 ,  3 ,  4 , and  5 are constants,   is the corner radius, and   is the die opening radius [2].
Based on the previous expressions in the function under the LDR, an expression depending on the LDR is obtained, and the nonlinear equations are solved by the Newton-Raphson method.In this way, it is possible to determine the values of LDR for each stage of redrawing, starting from an initial die radius and by the consideration that the die radius of each redrawing step will be reduced until 80% of the corresponding value of the previous phase.Once the LDR for each redrawing phase is found, the model determines the diameter of each stage.

Ironing.
The number of stages of ironing depends on the size of the final part.The model is based on the performance of three limiting conditions to be drawn from stage 2 to stage .The calculation of the process variable is the external diameter of the piece, according to the block diagram, illustrated in Figure 6.
By considering separately the merits of the three limiting conditions for ironing, the solution can be determined for phases 2 to .
The First Ironing Limiting Condition.The maximum ironing load ( ,max ) must be less than the cracking load (  ) of the material.This condition is expressed as the following: Second Ironing Limiting Condition.This condition assumes that the maximum load in the ironing process (  ) must be less than the yield tensile strength ((  )).This condition provides good results in parts of drawing [24] and it requires Third Ironing Limiting Condition.This refers to the value of the thickness reduction ratio limit ().This ratio is widely used in the calculation and design of the processes of drawing [5].The model performs the calculation for each  stage, which mainly depends on the stage of drawing and the type of material used.

Optimization Process
The model is based on the minimization of the total process time, which is the time to approach the punch from the initial stage, plus the time of the operation, plus the time to recover the punch to its initial position, according to the following objective function: where   is the objective function,   is the total process time,   is the total process time at stage ,  , is the approach time at stage ,  , is the operation time at stage , and  , is the recovery time at stage .
The process time at stage  is given by the expression where  −1 is the part length in stage  − 1,   is the punch length in stage ,   is part length in stage , V , is the approach velocity of the press in the stage , V , is the ironing velocity in the stage , and V , is the recovery velocity of the press in the stage .
As can be seen, the time is defined by means of the velocities involved in each operation and the punch or part length of each stage.Moreover, the part length in stage  is related to the bottom thickness, the blank diameter, the punch diameter, and the optimized diameter.Therefore, the model presents an optimization problem with an objective function which is subject to some constraints as follows: subject to where and   is the optimized diameter in stage ,   is the bottom thickness in stage ,  , is the punch diameter in stage , and Δ  is the drawing surplus ratio between the stages  and , defined by the expression  The optimization requires a distribution between the stages 1 to  of the drawing surplus ratio, Δ  .This distribution is conditioned to the minimization of the total process time as the sum of each stage time.Thus, an improvement is achieved in each drawing stage.In this manner, the optimized drawing conditions are separate from the maximum limits fixed by the initial drawing conditions.

International Journal of Manufacturing Engineering
The resolution of the recursive function is carried out according to the values of all the possible diameters between   and   (1 + Δ  ).The algorithm identifies which combination is the most suitable, and it selects the optimized diameter   in each stage, such that the condition of the objective function is satisfied; that is, the combination of diameters must give the minimum total process time.
The algorithm also permits the modification of the velocity parameters in each press required in the multistage process (approach velocity, operation velocity, and recovery velocity of the punch).Accordingly, it is possible to realise a more adequate distribution of available presses in the facility.

Algorithm Resolution Process
The algorithm proposed in Figure 7 performs the resolution in the following steps.
Step 1. Definition of incremental factor , given by the following expression: where Δ  is the surplus ratio of drawing and   is the number of times that the drawing surplus ratios are fractioned.
Step 2. Progressive increment of diameter of each process stage, from   to   (1 + ), where  is a parameter that varies from 0 to factor   .
Step 3. Resolution of the recursive function is based on obtaining all possible substages, by means of the modification of each stage diameters, from   to   (1 + Δ  ).
Step 4. Once all the possible diameters of each stage are defined, the algorithm searches the arrangement that allows for the minimizing of the total process time.
Step 5.The search stops and the best solution up to the current iteration is given as the output.

Experimental Results
The deep-drawing tool has allowed for the verification of this algorithm's integrity.Computational and experimental tests have been carried out in brass, in particular in UNS C26000 alloy, applied to three parts.The dimensions of the parts are shown in Table 2.
For the experimental results of the three cases presented in Table 2, its analytical resolution has been carried out by the software DDT according to the flowchart presented in Figure 3.The final result of using the software tool is an improved solution that optimizes the total process time and reduces the manufacturing cost of the presented industrial cases.These costs use the information from the results provided by the previous steps and are based on the total time of the deformation process and the labour cost per hour.Figure 8 shows the results of the tool for the initial solution corresponding to Case A.
Figure 9 presents the calculation of the optimized process by means of the DDT using the proposed algorithm.Values for   and   in operations of redrawing and ironing are given.
The resolution of the optimization process using the proposed algorithm requires the introduction of the machinery operation parameters involved in the process.
The computer tool requires the input of these data (Figure 10), according to the flowchart of the algorithm presented in Figure 3.
One of the most important variables in the analysis of such processes is the variable "ratio of wall thickness reduction." The optimization process allows the user to obtain a greater stability of the process maintaining the balance of the different phases.Table 3 shows the evolution of the drawing load for the industrial Cases A, B, and C under study.Table 4 shows the evolution of the wall thickness.The results show that the optimized process achieves a better balance so that it is more stable.This stability is transmitted in a better distribution of the capabilities of drawing, redrawing, and ironing between the different stages, subsequently resulting in improved processing time and reduced manufacturing cost.
Once the evolution of different variables that influence the process has been analyzed, the results for "process time" and "manufacturing cost" are presented in order to study the improvements that could be achieved by the implementation of the proposed algorithm.
Tables 5 and 6 show the evolution of these improvements on the variables "total process time" (  ) and "manufacturing cost" (  ) in terms of number of iterations performed by the algorithm.
Similarly, the evolution of "total process time" and the "manufacturing cost" for Cases B and C has been analyzed.Figures 11 and 12 show the improvements evolution obtained for the three Cases A, B, and C. It can be observed that the greater improvements are achieved in the first iterations (first five according to Case A), but the results improve further as the iteration increases.

Analysis and Discussion of the Results
The resolution of the industrial Cases A, B, and C by means of the DDT Algorithm presents important advantages.As shown in Tables 7, 8, and 9, the DDT algorithm produces improvements between 5.17% and 8.18% compared with the initial solution for the process time variable and between 4.40% and 7.78% for the manufacturing cost.In the same manner, the experimental results obtained an improvement between 6.55% and 9.34% for the process time and from 6.60% to 11.55% for the manufacturing cost.These improvements are more relevant in mass production.The most
presents the flow diagram of this redrawing algorithm.International Journal of Manufacturing Engineering

Figure 3 :
Figure 3: Flow diagram of the model.

Figure 4 :
Figure 4: Flow diagram of the drawing algorithm.

Figure 5 :
Figure 5: Flow diagram of the redrawing algorithm.

1 Figure 6 :
Figure 6: Flow diagram of the ironing algorithm.

Table 2 :
Dimensions of parts (in mm).
-Stage n (d n , l n , e n , DR n , K n , F e ) -Stage m (d m , l m , e m , DR m , K m , F e ) xe -Stages 1 to n (D n , L n , E n , DR n , K n and F e ) -Total optimized time, t e -Total optimized cost, C f Input L ui , ai , ei and s