# Optimization of the single point incremental forming process for titanium sheets by using response surface

Badreddine Saidi, Laurence Moreau-Giraud, Abel Cherouat

► To cite this version:

Badreddine Saidi, Laurence Moreau-Giraud, Abel Cherouat. Optimization of the single point incremental forming process for titanium sheets by using response surface. NUMIFORM 2016 - 12th International Conference on Numerical Methods in Industrial Forming Processes, Jul 2016, Troyes, France. pp.10011, 10.1051/matecconf/20168010011. hal-03623276

HAL Id: hal-03623276

https://hal.science/hal-03623276

Submitted on 24 Oct 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Optimization of the single point incremental forming process for titanium sheets by using response surface

Badreddine Saidi<fnref n="1" />,<fnref n="2" />*, Laurence Giraud-Moreau<fnref n="1" />, Abel Cherouat<fnref n="1" />

**Abstract.** The single point incremental forming process is well-known to be perfectly suited for prototyping and small series. One of its fields of applicability is the medicine area for the forming of titanium prostheses or titanium medical implants. However this process is not yet very industrialized, mainly due its geometrical inaccuracy, its not homogeneous thickness distribution... Moreover considerable forces can occur. They must be controlled in order to preserve the tooling. In this paper, a numerical approach is proposed in order to minimize the maximum force achieved during the incremental forming of titanium sheets and to maximize the minimal thickness. A surface response methodology is used to find the optimal values of two input parameters of the process, the punch diameter and the vertical step size of the tool path.

## 1 Introduction

The single point incremental forming process (ISF) is an emerging process with a high industrial interest. This process is based on the use of a spherical tool, which is moved along CNC controlled tool path (figure 1). During the incremental forming process, the sheet blank is fixed in sheet holder. The tool follows a certain tool path and progressively deforms the sheet. Because incremental forming is a dieless process, it is perfectly suited for prototyping and small volume production [1, 2].

Medical implants are one of the potential fields of application of the incremental forming process, due to the need of product customization to each patient. One of the major advantages of the ISF process is the ability to program several punch trajectories on the same machine in order to obtain different shapes perfectly suited to each patient. Titanium alloys are generally used by the biomedical industries because of their corrosive resistance and great biocompatibility with the human body. Titanium is recognized by its high strength-to-weight ratio. But its high cost and difficult formability compared with aluminum make that the incremental forming of titanium sheet is little studied. Before the incremental forming process of titanium sheets is really industrialized, studies are necessary. The ISF process indeed suffers from a big slowness, geometrical inaccuracy, a not homogeneous thickness distribution which reduces its industrial suitability [1, 3, 4, 5].

During this process, considerable forces can occur [5]. The control of the forming force is particularly important

when titanium sheets are used to ensure the safe use of a CNC milling machine and preserve the tooling and machinery [5, 6].

The aim of this paper is to determine input parameters of the process in order to minimize the maximum force achieved during the incremental forming of titanium sheets and to optimize the thickness reduction. Several studies were developed to optimize the incremental forming of aluminum sheets [2, 8]. A little less was done for titanium sheets.

In this paper, a numerical procedure is developed in order to study the influence of two input parameters, the punch diameter and the vertical step on the forming force and the minimal thickness. To conduct this study, an experimental design combined to response surface methodology is used [2]. The experimental design is based on finite element analysis realized with ABAQUS/Explicit. Mathematical relations are defined from the response surfaces to predict the forming force and the thinning rate as function of the two input parameters and to define the input parameters allowing to maximize the forming force and maximize the minimal thickness.

## 2 Incremental forming of titanium sheets

The single point incremental sheet forming (SPIF) process is illustrated on the figure 1. Different trajectories along which the tool moves have been studied in the literature. The tool path most frequently used and that the

<footnote n="1">¹Université de Technologie de Troyes France-GAMMA3, France</footnote>

<footnote n="2">²Ecole Nationale d'Ingénieurs de Tunis-MAI Tunisia</footnote>

<footnote>ª Corresponding author: laurence.moreau@utt.fr</footnote>
authors of the paper have already studied during experimental works [3] is composed of rotational movements in the horizontal plane combined with small steps in the vertical direction (see figure 2). This tool path is used in this study.

<figure><img src="image_2.png" /><figcaption>**Figure 1.** Incremental forming process.</figcaption></figure>

This paper concerns the incremental forming of titanium ASTM grade 2 sheets with a thickness of 0.5 mm. A truncated cone with a wall angle of 60°, a maximum diameter of 100 mm and a depth of 20 mm was chosen for this preliminary study about the optimization of the incremental forming process with titanium sheets. The sheet has a size of 300 mm x 300 mm.

<figure><img src="image_3.png" /><figcaption>**Figure 2.** Tool path.</figcaption></figure>

## 3 Numerical simulation

During the incremental forming process, the sheet is held firmly by the clamping system so that there is no slipping. A full model of blank holder modeled by rigid surfaces has been chosen to represent the experimental clamping system. A vertical pressure of 4.3MPA [9] has been applied. Figure 3 presents the numerical model at the beginning of the process. The initial and final mesh of the sheet for a punch diameter of 10mm and a step size of 1 mm is shown on figure 4.

<figure><img src="image_4.png" /><figcaption>**Figure 3.** Numerical model.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 4.** Initial and final mesh of the sheet.</figcaption></figure>

In what concerns material modeling, an elasto-plastic model was chosen for the simulation. The titanium behavior was modeled by means of a Swift type hardening law :

$$ \sigma = k (\varepsilon_0 + \varepsilon_p)^n \quad (1) $$

A tensile test allowed generating the strain-stress curve for the titanium sheet. Figure 5 presents this stress-strain curve.

<figure><img src="image_6.png" /><figcaption>**Figure 5.** Stress-strain curve.</figcaption></figure>

An inverse analysis using the data from this tensile test has been used to identify the mechanical parameters. Table 1 summarizes the mechanical properties of the material obtained thanks to the inverse method.

<table><caption>**Table 1.** Mechanical properties of the Titanium grade 2.</caption><tr><td>Modulus of elasticity : <i>E</i></td><td>100 Gpa</td></tr><tr><td>Yield strength: <i>σ</i></td><td>210 Mpa</td></tr><tr><td>Hardening coefficient: <i>n</i></td><td>0.084</td></tr><tr><td>Coefficient of proportionality: <i>K</i></td><td>697,7 MPa</td></tr><tr><td>Strain : ε<sub>0</sub></td><td>0.00157</td></tr></table>

## 4 Experimental design

The goal of this present study is to determine the values of the punch diameter (dp), and the vertical step size (Δz) which allow to minimize the forming force and maximize the minimal thickness. As the number of input variables is low, a full factorial design at three levels has been chosen to realize the experimental design. Table 2 presents the forming factors and their respected levels.

The numerical simulation of the incremental forming process has been performed with the ABAQUS explicit solver. The finite elements chosen to mesh the sheet are shell elements S4R with 4 nodes and reduced integration. The spherical tool is modeled by a rigid surface.
<table><caption>**Table 2.** Forming factors and levels.</caption><thead><tr><th>Factors</th><th>Factor notation</th><th>Level 1</th><th>Level 2</th><th>Level 3</th></tr></thead><tbody><tr><td>dp : Tool diameter (mm)</td><td>A</td><td>10</td><td>12.5</td><td>15</td></tr><tr><td>Δz : Step size (mm)</td><td>B</td><td>0.5</td><td>0.75</td><td>1</td></tr></tbody></table>

Nine numerical simulations of the forming process were carried out. The responses values are the minimal thinning and maximal forming force. Table 3 presents the numerical results performed according to the full factorial design.

<table><caption>**Table 3.** Numerical results.</caption><thead><tr><th>n°</th><th>Dp (mm)</th><th>Δz (mm)</th><th>Minimal thickness (mm)</th><th>Maximal forming force (N)</th></tr></thead><tbody><tr><td>1</td><td>10</td><td>0.5</td><td>0,255909</td><td>870</td></tr><tr><td>2</td><td>10</td><td>0.75</td><td>0,222006</td><td>950</td></tr><tr><td>3</td><td>10</td><td>1</td><td>0,172874</td><td>1030</td></tr><tr><td>4</td><td>12.5</td><td>0.5</td><td>0,258584</td><td>920</td></tr><tr><td>5</td><td>12.5</td><td>0.75</td><td>0,106811</td><td>980</td></tr><tr><td>6</td><td>12.5</td><td>1</td><td>0,197069</td><td>1110</td></tr><tr><td>7</td><td>15</td><td>0.5</td><td>0,205474</td><td>1060</td></tr><tr><td>8</td><td>15</td><td>0.75</td><td>0,205505</td><td>1120</td></tr><tr><td>9</td><td>15</td><td>1</td><td>0,274236</td><td>1180</td></tr></tbody></table>

## 5 Optimization results

A quadratic polynomial equation has been considered in order to fit the numerical results and get a relationship between the response values and the input parameters.

The two polynomial models (2) and (3) of second degree considered for the forming force and minimal thickness are the following [10]:

$$F_{\max} = a_0 + a_1 d_p + a_2 \Delta z + a_3 d_p^2 + a_4 d_p \Delta z + a_5 \Delta z^2 \quad (2)$$

$$STH = b_0 + b_1 d_p + b_2 \Delta z + b_3 d_p^2 + b_4 d_p \Delta z + b_5 \Delta z^2 \quad (3)$$

Table 4 presents the values of the coefficients, which allow to predict the forming force and the thinning rate according to the two input parameters. These numerical results have been obtained with ANOVA.

<table><caption>**Table 4.** Coefficients of the fit functions.</caption><thead><tr><th>a<sub>0</sub></th><th>a<sub>1</sub></th><th>a<sub>2</sub></th><th>a<sub>3</sub></th><th>a<sub>4</sub></th><th>a<sub>5</sub></th></tr></thead><tbody><tr><td>1082,2</td><td>- 80,7</td><td>233,3</td><td>5,1</td><td>- 16,0</td><td>186,7</td></tr><tr><th>b<sub>0</sub></th><th>b<sub>1</sub></th><th>b<sub>2</sub></th><th>b<sub>3</sub></th><th>b<sub>4</sub></th><th>b<sub>5</sub></th></tr><tr><td>1,48</td><td>- 0,12</td><td>- 1,40</td><td>0,0037</td><td>0,03</td><td>0,6</td></tr></tbody></table>

The 3D response surfaces showing the variations of the forming force and the minimal thickness with respect to the punch diameter and the step size are given in figure 6 and 7. These surfaces are useful to understand the

interaction properties between both the punch diameter and step size and the output parameters. We can see on these figures that the increase of the tool diameter and the step size induces an increase of force forming. As regards the minimal thickness, a nonlinear response is observed. The maximal thickness for the considered input parameters is obtained for the minimal tool diameter and minimal step size.

These response surfaces have been used to determine the optimal values of the tool diameter and the vertical step size which allow to minimize the forming force and maximize the minimal thickness. The results are presented in the table 5.

<figure><img src="image_7.png" /><figcaption>**Figure 6.** Response surface associated to the forming force.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 7.** Response surface associated to the minimal thickness.</figcaption></figure>

<table><caption>**Table 5.** Coefficients of the fit functions.</caption><thead><tr><th>Factor</th><th>Low</th><th>High</th><th>Optimum</th></tr></thead><tbody><tr><td>dp (mm)</td><td>10,0</td><td>15,0</td><td>10,0</td></tr><tr><td>Δz (mm)</td><td>0,5</td><td>1,0</td><td>0,5</td></tr><tr><td></td><td colspan="2">STH max (mm)</td><td>0,273923</td></tr><tr><td></td><td colspan="2">Fres min (N)</td><td>865,556</td></tr></tbody></table>

An analysis of variance was conducted with ANOVA in order to analyze the effects of the input parameters on the value of the forming force and minimal thickness. They are presented in figures 8 and 9.

A quasi linear effect of the punch diameter and the step size on the forming force is observed. Figure 9 shows a quadratic effect of the input parameters on the thickness response. Figures 10 and 11 presenting the standardized PARETO chart show that the punch diameter and the vertical step size have a similar effect on the forming force but the vertical step size has a more important effect than the punch diameter on the thickness.
<figure><img src="image_9.png" /><figcaption>**Figure 8.** Effects of the input parameters on the value of the forming force.</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Figure 9.** Effects of the input parameters on the value of minimal thickness.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Figure 10.** Standardized PARETO chart for force forming.</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Figure 11.** Standardized PARETO chart for the thickness.</figcaption></figure>

Figures 12 and 13 present the normal probability plot. The process parameters are considered significant to the experiment as there is a large angled slope of the straight line. We can see that all the points are relatively

close to the straight line and thus that the data follow a normal distribution [8].

<figure><img src="image_13.png" /><figcaption>**Figure 12.** Normal probability distribution for force forming.</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>**Figure 13.** Normal probability distribution for thickness.</figcaption></figure>

## 6 Conclusion

In this paper, a numerical approach based on finite element analysis and response surface methodology has been applied to optimize the incremental forming process of titanium sheets. A second-order quadratic function has been defined in order to predict the maximal forming force and minimal thickness as function of the punch diameter and the vertical step size.

The impact of the punch diameter and the vertical step size on the maximal forming force and minimal thickness has been studied.

The optimal values of the input parameters allowing minimizing the forming force and maximizing the minimal thickness have been defined from the response surfaces.

This study about the optimization of the incremental forming process with titanium sheets is a preliminary study. Only two input parameters were considered and the material low used in the finite element model doesn't take into account the general anisotropy. This study will be pursued in future works.

## References

1. G.Ambrogio, I. Corstantino, L. De Napoli, L. Filice, L. Fratini, M. Muzzupappa, J. of Materials Proc. Technol., **1153–154**, p. 501–507, (2004)
2. V. Gulati, A. Aryal, P. Katyal, A. Goswami, J.I. of The Institution of Engineers (India): Series C, (2015)

3. L. Giraud-Moreau, A. Cherouat, A. Merat *IDDRG 2014*, (2014)

4. J. Duflou, Y. Tunckol, A. Szekeres, P. Vanherck, J. of Mat. Proc. Technol., **1189**, p. 65–72, (2007)

5. R. Aerens, P. Eyckens, A. Van Bael, J. R. Duflou, Int. J. Adv. Manuf. Technol., (2009)

6. L. Giraud-Moreau, A. Cherouat, J. Zhang, H. Borouchaki, Key Engin. Mater., **554-557**, pp 1375-1381, (2013)

7. B. Saidi, A. Boulila, A. Mahfoubh, N. Rachid, Mechanics & Industry, **16**, n° 1410, (2015)

8. S. Kurra, N. H. Rahman, S. P. Regalla, A. K. Gupta, Comput. Science and Software Engin., (2015)

9. J. Belchior, Thèse de l'INSA Rennes, France, (2013)

10. C. Raju, S. Narayanan, Measurement, **78**, 296–308, (2016)