Upsetting Analysis of High-Strength Tubular Specimens with the Taguchi Method

In order to obtain input data for numerical simulations of tube forming, the material properties of tubes need to be determined. A tube tensile test can only be used to measure yield stress and ultimate tensile stress. For tubes with a large diameter/thickness ratio (D/t), tensile specimens are cut out and processed in a similar way as with sheet metal. However, for thin tubes with a diameter/thickness ratio below 10, the tensile specimens could not be cut out. The flow curve of the analyzed tube with a small diameter and D/t ratio of 7 was determined with a ring-shaped specimen. The experimental force-travel diagram was acquired. A reverse-engineering method was used to determine flow curves by numerical simulations. Using an L25 orthogonal array of the Taguchi method different flow curve parameters and friction coefficient combinations were selected. Tube upsetting with determined parameter combinations was performed with the finite element method. With analysis of variance influential equations among selected input parameters were determined for the force levels at six upsetting states. With the evaluation of known friction coefficients and flow curve parameters, K, n, and ε0 according to the Swift approximation were determined and proved by the final shape of the workpiece.


Introduction
Thin-walled tubular materials are used in several industrial fields, ranging from the automotive industry [1], piping systems, bent tubular holders, and chassis [2,3] to small products in shoe industries where toroidal products are required [4].Several cases of those applications seem to be technologically less demanding; however, their thin-walled structure increases the complexity of stable production.In order to achieve prescribed production tolerances, geometrical, and material parameters of the used tubular feedstock need to be well-defined.Preliminary studies of planned forming operations, as well as influences of their input parameters, are often studied in a digital environment by using the finite element method (FEM) in combination with sensitivity analysis of the studied process [5][6][7].Material data play a crucial role in: analyses of technological windows of produced parts, • stable mass production.
Therefore, reliable material data of the feedstock used are indispensable.
Considering the standardized testing of material properties, uniaxial tension, and/or compression tests are most commonly used [8][9][10].They often describe material behavior (yield stress, ultimate tensile stress, flow curve, material's anisotropy, etc.) for a majority of forming processes.In some cases, such basic testing procedures insufficiently determine the yield loci of the material used.Both tests could thus be additionally supported by other tests, such as shear test, bi-axial tension test [11], bulge test [12], tube torsion test, and determination of a forming limit diagram (FLD) depending on the application in which the analyzed material will be implemented.
A selection of the proper material testing is further defined by the type of the feedstock used.It is known that tensile, bi-axial tensile, as well as combined tensile-compression tests, are used for sheet metal while a pure compression test is mostly used for bulk materials.
As tubular material is often produced from a sheet metal, it needs to be tested for tension and compression stress states.The properties of a tubular material are still often determined by a uniaxial tensile test [13].In longitudinal directions, tensile specimens are comparable with standardized specimens until the tube radius is too small and the tube curvature influences the accuracy of the measured results.Furthermore, in some cases, a tube expresses a different yield curve in longitudinal and hoop directions due to orthotropic material properties.Wang and Lin have, therefore, proposed a modified tensile test for the hoop direction [14,15].To be able to perform this test, the tube diameter needs to be sufficient to obtain the smallest needed gauge section of 25 mm and a length-to-width ratio of at least 4 [14].
In the previous decade, several authors have analyzed the behavior of a tubular material, in particular in combination with a hydroforming process [15][16][17], the wrinkling of such material in a bending process [18][19][20] or in structural and safety applications for the automotive industry [21][22][23].They have analyzed the tubular material with a diameter-to-thickness (D/t) ratio between 10 and 150.The testing of a material with a diameter-to-thickness ratio lying below 10 is rarely found.Due to a low D/t ratio, such material calls for special testing procedures, in particular when the flow curves need to be determined.The mechanical properties of tubes are analyzed according to the EN 10002-1:2002 standard, in which only yield stress R p , ultimate tensile stress R m and total elongation at fracture A t are determined.However, numerical analyses of a bending process by using the finite element method (FEM) also call for the data of a flow curve, which could not be obtained with a tensile test of the tube itself.A tri-axial stress state appears at the tensile loading of the tubular material.The inner diameter of a tube cannot be measured during the tensile test; the cross-sectional area necessary to calculate the true stress is thus not measurable.Therefore, the real cross-section of the tube could not be determined, and the procedure can be used only for the determination of yield stress R p and ultimate tensile stress R m .This paper presents a determination of a flow curve of a tubular material with an increased strength of E420 + CR2 quality for the D/t ratios of 7. Since the D/t ratio is too small for a tensile specimen to be cut out, a reverse engineering approach was selected aimed at determining the flow curve parameters.

Material Properties
The analyzed tubular material of E420 + CR2 quality, being in line with EN 10305-3:2003, is tested according to the EN 10305-3 standard with uniaxial tensile loading for testing sheet metal.Since the material is used for applications in which an increased strength is demanded, a minimal yield strength of 420 MPa is required as evident from the material nomenclature E420.The CR2 sign represents improved feedstock tolerances being ±0.05 mm of the nominal tube outer diameter for the analyzed tube types and sizes.As already mentioned in the introduction, a conventional tensile test cannot be used for a tubular material with such small D/t ratios.For the experiments and determination of flow curves by reverse engineering the tube with a nominal diameter of φ13.93 ± 0.05 mm × 2 ± 0.08 mm was selected.The D/t ratio was 7 for all analyzed specimens.
Before determining a flow curve with reverse engineering, it is necessary to experimentally obtain the mechanical properties of the tube.The experiments were performed on the universal testing machine AMSLER equipped with force and travel transducers to obtain yield stress R p , ultimate tensile stress R m , and total elongation A t .The mechanical properties of the analyzed material together with  To determine the flow curve of precision tubes of analyzed specimens a tube upsetting test combined with reverse-engineering of this test analyzed with the FEM program ABAQUS and statistical evaluation was performed.

Modified Upsetting Test
The widely used upsetting test allows a precise determination of mechanical and forming properties at a compression stress-strain state.The influence of contact friction is minimized by a special specimen shape for the Rastegaew test [25] or by lubrication with a polytetrafluoroethylene (PTFE) foil when a simple cylindrical specimen is used.However, it is not possible to perform a traditional billet upsetting test on a thin-walled tubular material due to the buckling phenomenon.Preparation of cylindrical billets from a tube wall is very difficult, expensive, and, because of their extremely small size, was also strongly influenced by the size effect.
To obtain a flow curve of the analyzed tube, the upsetting test was substituted with the modified ring compression test.Originally, the ring compression test with a standardized ratio of inner diameter, outer diameter, and specimen's height is used to obtain friction among the formed material and tool surface in bulk metal-forming processes.The geometrical ratio is influenced by the expected level of contact friction [26] having an optimal value of d o :d i :h = 6:3:2 [27,28] or 18:9:6 mm.However, the upsetting of the analyzed thin-walled tube does not correspond to the ratio d o :d i = 6:3.In our analyses of the tubular material, this ratio corresponds to the ratio 6:4.28.Since the area in contact with the tool is much smaller regarding the standardized ring compression test, it can be expected that the influence of the friction coefficient does not have a major impact on the forming process.However, this impact needs to be analyzed in greater detail with numerical simulations.

Experimental Work
The upsetting of the tubular specimen was performed on a universal testing machine AMSLER (Zwick Roell AG, Ulm, Germany) with a nominal force of 300 kN.The forming force was measured with a 500 kN load cell with an accuracy of ±450 N while the displacement was acquired with an incremental transducer having an accuracy of ±0.01 mm.Load-stroke diagrams were obtained by using a specialized program for signal processing [24].
The specimens of the tube with a wall thickness of 2 mm and an initial height of 10 mm were upset up to the height of 5.81 ± 0.02 mm.This upsetting height represents a true equivalent strain of ε e = 0.543.To minimize the friction between the tubular specimen and both tooling plates, the PTFE foil having a thickness of 0.2 mm was applied as a lubricant.
The experiments were performed with three samples.Their outer diameters were in the range of ±0.015 mm and inner diameters in the range of ±0.035 mm.The force-travel diagram used in further analyses represents the average value of all experiments.The standard deviation was in the range of ±2% for the whole range of force course.
After upsetting, the specimens were halved by wire electro-discharge machining (EDM) and sealed into a billet-shaped resin for their easier handling.To acquire the cross-sections of the specimens needed for their comparative evaluation with FEM results, they were scanned on the x-y table with Since the measuring equipment does not have a fix focus and magnification, it is necessary to calibrate it prior to measurements.Therefore, the measurement calliper with a thickness of 1.00 mm was scanned together with the cross-sections of the upset specimens.With the use of image processing [29], the acquired pictures were aligned according to the X-Y coordinate system and contour points were acquired for comparison with the FEM results as shown in Figure 1.Since the measuring equipment does not have a fix focus and magnification, it is necessary to calibrate it prior to measurements.Therefore, the measurement calliper with a thickness of 1.00 mm was scanned together with the cross-sections of the upset specimens.With the use of image processing [29], the acquired pictures were aligned according to the X-Y coordinate system and contour points were acquired for comparison with the FEM results as shown in Figure 1.Friction represents an important factor in the shaping of a specimen in the upsetting process.Its influence increases with the size of the surfaces that are in contact during the forming process.This phenomenon is, in particular, necessary for bulk forming processes.However, in the case of axial compressive loading of a tubular specimen presented in the paper, it is assumed the only minor role of the friction coefficient on the tube shaping resulting from a small contact surface needs to be considered by a numerical method and statistical Taguchi approach.
The friction coefficient of the PTFE foil used in experimental work was determined with a standardized ring upsetting test [30,31].According to the well-known diagram with calibration curves for various Coulomb's friction coefficients μ the applied foil expressed a friction value between μ = 0.03 and μ = 0.04 [32].Therefore, the evaluation of friction coefficient for analysis of variance (ANOVA) with the target value of μ = 0.035 was selected.

FEM Model
The upsetting process of the tubular specimens was analyzed with the ABAQUS commercial program [33] using the static implicit solver.Due to the axis-symmetric nature of the problem only a 2D-analysis with rotational symmetry was modelled (Figure 2).The FEM model consists of a deformable elastic-plastic specimen following Swift's law [34] and analytical rigid tool parts.The analyzed cross-section was meshed with mapped FEM elements with a minimal 0.15 mm size.According to the initial size of the tested tubular specimen, the mesh consists of 952 (68 × 14) nodes and 871 (67 × 13) elements of the CAX4R type.The Coulomb friction model was used.Apart from the variation of the parameters of the flow curve yield stress Rp, material constant K, and hardening exponent n other material parameters used in the FEM model are as follows: density ρ = 7850 kg/m 3 , Young's modulus E = 210 GPa, Poisson's ratio ν = 0.3.According to the experiment, the bottom tool plate was fixed and the upper tool plate moved with a speed of 0.84 mm/s.The maximal tool displacement was 4.19 mm.Friction represents an important factor in the shaping of a specimen in the upsetting process.Its influence increases with the size of the surfaces that are in contact during the forming process.This phenomenon is, in particular, necessary for bulk forming processes.However, in the case of axial compressive loading of a tubular specimen presented in the paper, it is assumed the only minor role of the friction coefficient on the tube shaping resulting from a small contact surface needs to be considered by a numerical method and statistical Taguchi approach.
The friction coefficient of the PTFE foil used in experimental work was determined with a standardized ring upsetting test [30,31].According to the well-known diagram with calibration curves for various Coulomb's friction coefficients µ the applied foil expressed a friction value between µ = 0.03 and µ = 0.04 [32].Therefore, the evaluation of friction coefficient for analysis of variance (ANOVA) with the target value of µ = 0.035 was selected.

FEM Model
The upsetting process of the tubular specimens was analyzed with the ABAQUS commercial program [33] using the static implicit solver.Due to the axis-symmetric nature of the problem only a 2D-analysis with rotational symmetry was modelled (Figure 2).The FEM model consists of a deformable elastic-plastic specimen following Swift's law [34] and analytical rigid tool parts.The analyzed cross-section was meshed with mapped FEM elements with a minimal 0.15 mm size.According to the initial size of the tested tubular specimen, the mesh consists of 952 (68 × 14) nodes and 871 (67 × 13) elements of the CAX4R type.The Coulomb friction model was used.Apart from the variation of the parameters of the flow curve yield stress R p , material constant K, and hardening exponent n other material parameters used in the FEM model are as follows: density = 7850 kg/m 3 , Young's modulus E = 210 GPa, Poisson's ratio ν = 0.3.According to the experiment, the bottom tool plate was fixed and the upper tool plate moved with a speed of 0.84 mm/s.The maximal tool displacement was 4.19 mm.

Design of Experiments
The Taguchi method is an engineering method for analyzing process parameters, which was developed by Genichi Taguchi.It can be used to obtain the relationship between the observed manufacturing process and its products [5,35].The major advantage of the Taguchi method is the ability to analyze the observed process with a small number of experiments.It is a simple and systematic design of experiment tools, which is often used in statistically robust design [6,36,37].
In the presented paper, the Taguchi method was used to analyze the influence of the material parameters K, n, and Rp as well as the Coulomb's friction coefficient μ on the forming force course and the final cross-sectional shape of the specimen.The parameters used in FEM simulations were based on the Taguchi method L25 orthogonal array with four influential factors, in which each factor was presented on five levels.

Selection of Flow Curve Parameters
The analyzed tube is directly rolled from a steel strip to a corresponding diameter and welded with electrical resistance welding.After removing the outer welding burr and cooling the hollow sections are calibrated to precise, standardized dimensions [38].Since the material is already pre-deformed, the use of the Swift flow curve approximation [34] was indispensable.The material description follows the equation: where K is a material constant, n is a hardening exponent, and ε0 is the initial shift of the flow curve to obtain the yield stress Rp at a declared value of εe = 0.002.Since the yield stress is already pre-defined by the experiments, the value of ε0 needs to be recalculated for each pair of parameters K and n in a reverse-engineering analysis.Digital analyses of upsetting the tubular specimens were performed with various flow curve parameters.The range of flow curve parameters was selected around the values known for sheet metal E420.Therefore, the hardening exponent n was selected in a range from 0.08 to 0.36 with an increment of 0.07.Similarly, the material constant K was increased from the lower value of 650 MPa to highest value of 950 MPa with an increment of 75 MPa.In order

Design of Experiments
The Taguchi method is an engineering method for analyzing process parameters, which was developed by Genichi Taguchi.It can be used to obtain the relationship between the observed manufacturing process and its products [5,35].The major advantage of the Taguchi method is the ability to analyze the observed process with a small number of experiments.It is a simple and systematic design of experiment tools, which is often used in statistically robust design [6,36,37].
In the presented paper, the Taguchi method was used to analyze the influence of the material parameters K, n, and R p as well as the Coulomb's friction coefficient µ on the forming force course and the final cross-sectional shape of the specimen.The parameters used in FEM simulations were based on the Taguchi method L 25 orthogonal array with four influential factors, in which each factor was presented on five levels.

Selection of Flow Curve Parameters
The analyzed tube is directly rolled from a steel strip to a corresponding diameter and welded with electrical resistance welding.After removing the outer welding burr and cooling the hollow sections are calibrated to precise, standardized dimensions [38].Since the material is already pre-deformed, the use of the Swift flow curve approximation [34] was indispensable.The material description follows the equation: where K is a material constant, n is a hardening exponent, and ε 0 is the initial shift of the flow curve to obtain the yield stress R p at a declared value of ε e = 0.002.Since the yield stress is already pre-defined by the experiments, the value of ε 0 needs to be recalculated for each pair of parameters K and n in a reverse-engineering analysis.Digital analyses of upsetting the tubular specimens were performed with various flow curve parameters.The range of flow curve parameters was selected around the values known for sheet metal E420.Therefore, the hardening exponent n was selected in a range from 0.08 to 0.36 with an increment of 0.07.Similarly, the material constant K was increased from the lower value of 650 MPa to highest value of 950 MPa with an increment of 75 MPa.In order to analyze the influence of friction, the Coulomb's friction coefficient in a wide range from µ = 0.025 to µ = 0.205 with an interval of 0.045 was selected.Simulations were performed based on the combinations of the parameters K, n, R p , and the friction coefficient µ using the Taguchi method designed on five levels as presented in Table 2.

Determination of Impacts on Force Course
The influence of particular parameters K, n, R p , and µ on the course of the forming force was analyzed at six points (Figure 3) defining different upsetting levels.The forces were labelled from F 1 to F 6 according to the tool displacement (Table 3).According to the selected values of input parameters (marked with blue in Table 4) following the Taguchi L25 matrix, the FEM simulations were performed.Forces F1 to F6 were extracted from each analysis (Table 4).According to the selected values of input parameters (marked with blue in Table 4) following the Taguchi L 25 matrix, the FEM simulations were performed.Forces F 1 to F 6 were extracted from each analysis (Table 4).Using the ANOVA concept, the influences of the parameters K, n, R p , and µ, as well as their two-factor products on particular observed force states were studied, and response surface was determined.As a measure of parameter importance the F-test was performed.The significance of a particular parameter was analyzed with the level of "p-value" below 5%, marked with green fields in Table 5.It can be found that parameter significance varies at different levels of forming force F 1 to F 6 .
It can be observed that the material parameter K is influential at the forces of all six upsetting states, whereas the parameter n is significant only in the first three force states.The influence of the material parameter R p is significant in all six states, but its influence decreases at the last two states.The Coulomb's friction coefficient µ does not have a significant influence on the first force state; however, from the second to the third stated its influence increases, and then from the fourth to the last force state it decreases again, although its influence is still significant.This shows an important tribological influence on the force-travel course despite the relatively small contact surface of the compressed ring specimen.
With all six force states, F 1 to F 6 , the influences of the products K•n and n•R p have also been shown as significant (Table 5).
The performed ANOVA delivers the following equations for the upsetting forces F 1 to F 6 : Metals 2016, 6, 257 9 of 14 The obtained influential equations of forming forces F 1 to F 6 represent the basis for the determination of flow curve parameters when they are compared with experimentally obtained values.

Multi-Parameter Analysis for an Inverse-Engineering Approach
The reverse engineering approach of flow curve determination is based on multi-parameter optimization based on equations F 1 to F 6 (Equations ( 2)-( 7)).The range for particular input parameters can be narrowed down to the experimentally expected range.The real experiment of the tube compression was compared with corresponding FEM analysis.For FEM simulation input data based on the multi-functional parameter determination having the material properties described in Section 2 were selected.Through this, the value of yield stress R p is fixed to 572 MPa.In the ANOVA it was observed that the initial assumption of the negligible influence of friction coefficient was wrong.Therefore, the proper selection of friction coefficient needs to be considered in further analysis.Values between µ = 0.03 and µ = 0.04 are selected according to the friction measurement of PTFE foil from the literature [32].Through this, the unknown parameters are limited to flow curve parameters K and n only.The multifunctional analysis can be solved analytically or graphically.However, the analytical determination delivers more accurate results and was therefore selected as a preferred solution.The graphical solution would be used only in the case of the unsolvable analytic solution.The experimentally obtained target values of F 1 to F 6 are presented in Table 6.To perform the multifunctional analysis, the spread range for each experimentally determined force F 1 to F 6 is to examine the procedure used for determination of K, n, and µ was based on the spread of F 1 in the 4% range, and F 2 to F 6 in the range of 2% according to the experimentally obtained values.The analytical, multifunctional optimization of all six target equations (Equations ( 2)-( 7)) finally delivers four different values of K and n as presented in Table 7.Despite that, all four results in Table 7 fit the prescribed criteria of the multifunctional analysis; the solution in the first row is the most optimal result since the used ANOVA optimization delivers the best solution first.The value of constant ε 0 is calculated from The obtained values of flow curve parameters were inserted into the FEM model, and the upsetting process in a digital environment was re-run.The FEM simulation was re-run with the material parameters following the Swift equation according to the case number 1 from the Table 7: Comparison of experimental and FEM predicted force course during the upsetting process together with their differences are presented in Figure 4.The differences are in a close range lying below the 5% limit.The value of constant ε0 is calculated from The obtained values of flow curve parameters were inserted into the FEM model, and the upsetting process in a digital environment was re-run.The FEM simulation was re-run with the material parameters following the Swift equation according to the case number 1 from the Table 7: σf = 713•(0.321+ εe) 0.195 (9) Comparison of experimental and FEM predicted force course during the upsetting process together with their differences are presented in Figure 4.The differences are in a close range lying below the 5% limit.7.

Comparison of Tube Shaping
Determination of flow curves by the reverse engineering method led to the quality determination of the force-travel course of the upset specimens.However, it is necessary to prove whether the specimens with the Taguchi method, ANOVA, multifunctional optimization, and FEM determined material parameters are in fact shaped in the same manner as the experimental ones.
In addition to the evaluation of forming force F versus punch displacement, the specimen shape was observed.This research analyzed the specimen shape at various compression states for the  7.

Comparison of Tube Shaping
Determination of flow curves by the reverse engineering method led to the quality determination of the force-travel course of the upset specimens.However, it is necessary to prove whether the specimens with the Taguchi method, ANOVA, multifunctional optimization, and FEM determined material parameters are in fact shaped in the same manner as the experimental ones.
In addition to the evaluation of forming force F versus punch displacement, the specimen shape was observed.This research analyzed the specimen shape at various compression states for the determination of material parameters K, n, and R p .The simulated variation of material parameters was first selected according to Equation ( 9) and analyzed for three different values of the Coulomb's friction coefficient: µ = 0.035, µ = 0.045, and µ = 0.15.The influence of material parameters on the specimen shape was also analyzed with 50% higher material parameters than those determined in Equation ( 9): Figure 5 shows that the friction coefficient has a measurable influence on part shaping.According to FEM simulation, the changes of the shape (diameter) are indeed minor regarding the material parameters and they are for the specimens having diameter of 13.93 mm in the range below ±0.05 mm.Based on this information, it was concluded that the determination of material parameters is significantly more reliable and accurate according to the analysis of the forming force.
To obtain comparative contour data of the upset tubular specimen acquired by the CCD camera as presented in Figure 1 and matched with final FEM results, a specialized program was written [24].The contour-recognition program orients the specimen according to the x-y axis, compares the evaluated sample with the scanned calibre of 1.00 mm, takes out the contour points of the specimen, and compares them with the coordinates of the FEM mesh extracted from the nodes after the simulation.σf = 1070•(0.468+ εe) 0.29 (10) Figure 5 shows that the friction coefficient has a measurable influence on part shaping.According to FEM simulation, the changes of the shape (diameter) are indeed minor regarding the material parameters and they are for the specimens having diameter of 13.93 mm in the range below ±0.05 mm.Based on this information, it was concluded that the determination of material parameters is significantly more reliable and accurate according to the analysis of the forming force.To obtain comparative contour data of the upset tubular specimen acquired by the CCD camera as presented in Figure 1 and matched with final FEM results, a specialized program was written [24].The contour-recognition program orients the specimen according to the x-y axis, compares the evaluated sample with the scanned calibre of 1.00 mm, takes out the contour points of the specimen, and compares them with the coordinates of the FEM mesh extracted from the nodes after the simulation.
As is evident from Figure 6, the specimen correlates well with the FEM results for the observed specimen height at the end of the forming process.However, minor discrepancies between the upset specimens and FEM results can be observed at the inner contour of the specimen and upper part of the outer bulged specimen's area.These differences are not significant; however, influential parameters on part shaping are to analyze more into the detail in the future.Currently, in FEM simulations only the fixed value of friction coefficient μ was considered.In particular, the influence of the friction coefficient during the process is to study in order to evaluate its influence on the part shaping.As is evident from Figure 6, the specimen correlates well with the FEM results for the observed specimen height at the end of the forming process.However, minor discrepancies between the upset specimens and FEM results can be observed at the inner contour of the specimen and upper part of the outer bulged specimen's area.These differences are not significant; however, influential parameters on part shaping are to analyze more into the detail in the future.Currently, in FEM simulations only the fixed value of friction coefficient µ was considered.In particular, the influence of the friction coefficient during the process is to study in order to evaluate its influence on the part shaping.

Conclusions
In this paper, the determination of a flow curve for thin-walled tubular specimens was implemented.The analyses were based on a modified upsetting test of tubular specimens with a D/t ratio of 7. Using the Taguchi method and analysis of variance (ANOVA), the material's flow curve parameters were determined.The main conclusions are the following: 1) The forming force versus punch displacement was used for a reverse engineering approach to determine the flow curve parameters since it was found that the variation of material parameters has only a minor role on the specimen shaping.
2) The optimization approach based on the comparison of experimentally and statistically determined values (Taguchi + ANOVA) of forming force at six upsetting levels delivered the parameters of the material's flow curve K = 713 MPa, n = 0.195, and ε0 = 0.321 according to Swift's power law.The friction coefficient μ in this optimization was selected according to the

Conclusions
In this paper, the determination of a flow curve for thin-walled tubular specimens was implemented.The analyses were based on a modified upsetting test of tubular specimens with a D/t ratio of 7. Using the Taguchi method and analysis of variance (ANOVA), the material's flow curve parameters were determined.The main conclusions are the following: 1) The forming force versus punch displacement was used for a reverse engineering approach to determine the flow curve parameters since it was found that the variation of material parameters has only a minor role on the specimen shaping.

2)
The optimization approach based on the comparison of experimentally and statistically determined values (Taguchi + ANOVA) of forming force at six upsetting levels delivered the parameters of the material's flow curve K = 713 MPa, n = 0.195, and ε 0 = 0.321 according to Swift's power law.The friction coefficient µ in this optimization was selected according to the one measured for the PTFE foil.The yield stress of the tube R p was also determined experimentally.

3)
The comparison of the FEM-predicted force course based on the determined material flow curve and the experimental force course showed that the differences are below 5%.

4)
The comparison of the cross-sectional contours of the upset specimens with the numerically obtained results for the determined flow curve according to Swift's law showed only minor shape deviations.

5)
Future analyses should focus on the geometrical optimization of the specimens regarding the d 0 :h 0 ratio for various D/t ratios.The presented methodology needs to be further tested regarding its robustness according to all input parameters, such as the level of the flow curve, the value of R p , and the friction coefficient.

a
CCD camera.Concentric lighting was provided by a flexible optical fiber cable and a SCHOTT FOSTEC DCR III halogen light source.

Figure 1 .
Figure 1.Acquisition of the specimen contour by a CCD camera.

Figure 1 .
Figure 1.Acquisition of the specimen contour by a CCD camera.

Figure 2 .
Figure 2. Rotational-symmetric FEM model of the tube (left) in a 3D-view and the results of the equivalent strain distribution obtained by a FEM simulation in a 3D-view (right).

Figure 2 .
Figure 2. Rotational-symmetric FEM model of the tube (left) in a 3D-view and the results of the equivalent strain distribution obtained by a FEM simulation in a 3D-view (right).

Figure 4 .
Figure 4. Comparison of force course for selected material parameters from Table7.

Figure 4 .
Figure 4. Comparison of force course for selected material parameters from Table7.

Figure 5 .
Figure 5. Influence of the material parameters and friction coefficient on the forming force (left) and specimen shape (right).

Figure 5 .
Figure 5. Influence of the material parameters and friction coefficient on the forming force (left) and specimen shape (right).

Figure 6 .
Figure 6.Comparison between the part shape and FEM results obtained by multi-parameter optimization.

Figure 6 .
Figure 6.Comparison between the part shape and FEM results obtained by multi-parameter optimization.

Table 1 .
As evident from Table1, the tested material is in line with the EN 10305-3:2002 standard regarding the prescribed mechanical properties.

Table 2 .
Combinations of the analyzed parameters K, n, R p , and µ.

Table 4 .
FEM -obtained results of forming forces at six observed upsetting levels.

Table 3 .
Forces F 1 to F 6 and corresponding punch displacements.

Table 4 .
FEM-obtained results of forming forces at six observed upsetting levels.

Table 5 .
The influence of the parameters K, n, R p , and µ, and their products on different levels of upsetting force.

Table 6 .
Experimentally obtained values of F 1 to F 6 .

Table 7 .
Determined combination of material parameters and friction coefficient µ.