Sensitivity Analysis of the Square Cup Forming Process Using PAWN and Sobol Indices

: This study investigates the sensitivity of the square cup forming process. It analyses how the uncertainties in the material properties, friction and process conditions affect the results of the square cup, such as equivalent plastic strain, geometry change, thickness reduction, punch force and springback. The cup flange and the die curvature region are identified as highly affected areas, while the cup bottom is least affected by the uncertainties. Two sensitivity analysis techniques, PAWN and Sobol indices, are compared. In particular, the study shows that PAWN indices require a significantly smaller number of simulations than Sobol indices, making them a more efficient choice for sensitivity analysis. While both PAWN and Sobol indices generally give comparable results, discrepancies arise in the analysis of springback, where PAWN indices show superior accuracy, particularly when dealing with multimodal distributions. This observation highlights the importance of selecting the appropriate sensitivity analysis method based on the nature of the data being analysed. These results provide insights for optimizing stamping processes to reduce production time and costs.


Introduction
Metal forming processes are widely used in automotive, aerospace, and metalworking industries due to their ability to produce complex shapes with high precision and efficiency [1].These processes cover a wide range of techniques, including stamping, forging and extrusion, each tailored to specific applications and material properties.With the increasing demand for lightweight components with improved mechanical properties, the optimization of metal forming processes has become critical to maintaining competitiveness in the global marketplace.
Finite Element Analysis (FEA) plays a key role in predicting deformation processes and identifying factors that limit component formability.It enables the simulation of complex forming processes and the performance evaluation of different tool configurations and process parameters.However, conventional FEA approaches often overlook the inevitable uncertainties present in real-world industrial environments, a concern highlighted by several researchers [2,3].These uncertainties arise from a variety of sources, including material variability, geometric imperfections, and process variations, and can have a significant impact on the reliability and robustness of forming simulations [4][5][6].In addition, the authors in [2] pointed out that these uncertainties can change over time (e.g., due to tool wear).Ignoring these uncertainties leads to sub-optimal or unreliable forming process designs [7], Metals 2024, 14, 432 2 of 18 which in turn has a significant impact on component quality, scrap rates and manufacturing costs, for example due to downtime for troubleshooting and maintenance [8].
With the increasing availability of large amounts of data and improved computational capabilities, interest in stochastic modelling and quantification of uncertainties in sheet metal forming processes have been growing [9][10][11][12].Several researchers have modelled and quantified the influence of different uncertainty sources on the final product, by resorting to distinct methods, such as Monte Carlo simulation [13,14], the design of experimental techniques [15,16] and metamodels [17,18].In [12], an uncertainty analysis combining FEA and Monte Carlo simulation was designed to evaluate the influence of uncertainty sources on the end product.Also, in [16], the influence of the material and friction variability on the geometry of a U-channel was studied using an analysis of variance (ANOVA), with them reaching the conclusion that the parameters C and n of the Swift hardening law were the most influential factors in the springback and maximum thinning results.In another piece of research, [18], a metamodel was used based on machine learning techniques to predict the occurrence of defects, considering the variability in the material properties and process parameters.In [17], the authors compared the performance of various metamodel techniques by considering the uncertainties in the material behaviour and sheet thickness.The authors concluded that some metamodel techniques are able to accurately predict the forming process results.The application of any of the presented methods requires the execution of several numerical simulations of the forming process.To reduce the computational cost of these strategies, screening techniques can be used to select the main sources of uncertainty to be considered in the stochastic model [16].
A sensitivity analysis is traditionally performed using analysis of variance (factorial ANOVA).However, this approach assumes a unimodal distribution of results, which is not necessarily true.Factorial ANOVA also requires a model to be predefined to establish a correlation between the factors and the results [16].Sobol indices (also referred to as functional ANOVA [19]) can address this problem by not requiring a hypothesis in the form of a model, but have the limitation of only being reliable for unimodal data.The PAWN indices appear as a technique that can overcome both of these vulnerabilities [20,21], as they do not rely on a statistical second-order moment (i.e., variance) [22].Instead, PAWN indices use the Cumulative Distribution Functions to evaluate the influence of the different factors on the results.
Although PAWN sensitivity analysis can present a promising alternative to the abovementioned approaches, its application in the sheet metal forming process is still unexplored.The goal of this work is to explore the use of PAWN sensitivity analysis, by investigating the sensitivity of the square cup forming process to material properties, friction, and process conditions.The results of the PAWN indices will be compared to the Sobol indices, obtained in a previous sensitivity analysis [23].By employing sensitivity analysis techniques such as PAWN and Sobol indices, the aim is to identify key parameters influencing the forming process of a square cup and provide insights for optimizing process parameters to mitigate the effects of uncertainties.

Numerical Model
A schematic representation of the square cup forming process is shown in Figure 1, with tool geometries based on the benchmark proposed by the NUMISHEET' 93 conference [24].Throughout the process, a constant blank holder force (BHF) is used to control the material flow.The punch is then advanced to a total displacement of 40 mm.At the end, the square cup is extracted from the tools, promoting the springback of the square cup.To optimize the computational efficiency, only a quarter of the model is simulated, taking advantage of material and geometric symmetries.The square blank, with an initial thickness of t 0 and a side length of 75 mm, is discretized into 1800 elements using 8-node hexahedral solids, with 2 elements in thickness and 30 elements per side.This mesh guarantees an error inferior to 1% in the maximum yield stress and a simulation time inferior to 5 min (as can be concluded from Figure 1c).The interaction between the blank and the tools is governed by Coulomb's law, using a constant coefficient of friction µ 0 .The contact with friction problem is solved using the augmented Lagrangian method.The tools are rigid, represented using Nagata patches [25], and are only allowed to displace on the vertical direction, being the displacement of the blank holder controlled by the imposed force (BHF).Computational simulations were performed using the in-house code DD3IMP (Deep Drawing 3D Implicit Code) [26], which uses an updated Lagrangian scheme to integrate the constitutive law in an implicit way.All numerical simulations were performed on computers equipped with an Intel ® Core™ i7-8700K Hexa-Core processor 4.7 GHz.On average, each simulation took approximately 4 min and 34 s to complete.
Metals 2024, 14, x FOR PEER REVIEW 3 of 21 hexahedral solids, with 2 elements in thickness and 30 elements per side.This mesh guarantees an error inferior to 1% in the maximum yield stress and a simulation time inferior to 5 min (as can be concluded from Figure 1c).The interaction between the blank and the tools is governed by Coulomb's law, using a constant coefficient of friction  .The contact with friction problem is solved using the augmented Lagrangian method.The tools are rigid, represented using Nagata patches [25], and are only allowed to displace on the vertical direction, being the displacement of the blank holder controlled by the imposed force ( ).Computational simulations were performed using the in-house code DD3IMP (Deep Drawing 3D Implicit Code) [26], which uses an updated Lagrangian scheme to integrate the constitutive law in an implicit way.All numerical simulations were performed on computers equipped with an Intel ® Core™ i7-8700K Hexa-Core processor 4.7 GHz.On average, each simulation took approximately 4 min and 34 s to complete.A low-carbon steel (DC06) was considered as the material of the sheet.The plastic behaviour is described by the Hill'48 orthotropic yield criterion [27] and the Swift work hardening law [28].The Hill'48 orthotropic yield criterion is defined as follows: where  is the yield stress; , , , ,  and  are the parameters that define the shape of the yield surface; and  ,  ,  ,  ,  and  are components of the Cauchy stress tensor, wri en in the orthotropic coordinate system 0xyz.In this work, it is assumed that  =  = 1.5 (identical to von Mises) and the condition  +  = 1, meaning that the yield stress, , is comparable to the uniaxial tensile stress aligned with the rolling direction.
The anisotropy coefficients for 0°, 45° and 90° (relative to the rolling direction), respectively denoted by  ,  and  , can be determined by following the equations: A low-carbon steel (DC06) was considered as the material of the sheet.The plastic behaviour is described by the Hill'48 orthotropic yield criterion [27] and the Swift work hardening law [28].The Hill'48 orthotropic yield criterion is defined as follows: where Y is the yield stress; F, G, H, L, M and N are the parameters that define the shape of the yield surface; and σ xx , σ yy , σ zz , τ xy , τ xz and τ yz are components of the Cauchy stress tensor, written in the orthotropic coordinate system 0 xyz .In this work, it is assumed that L = M = 1.5 (identical to von Mises) and the condition G + H = 1, meaning that the yield stress, Y, is comparable to the uniaxial tensile stress aligned with the rolling direction.The anisotropy coefficients for 0 • , 45 • and 90 • (relative to the rolling direction), respectively denoted by r 0 , r 45 and r 90 , can be determined by following the equations: The yield stress (Y ) evolution with plastic deformation is described by the Swift hardening law: where C, ε 0 and n are material constants, and ε p represents the equivalent plastic strain derived from the Hill'48 orthotropic yield criterion by assuming an associated flow rule.The initial yield stress, Y 0 , is given by: This study will focus on investigating the influence of the uncertainty in the material properties, friction and the forming process conditions on the stamping of the square cup.Therefore, the input parameters of the numerical model whose influence will be studied are Young's modulus (E), Poisson's ratio (ν), anisotropy coefficients (r 0 , r 45 and r 90 ,), the initial yield strength (Y 0 ), the hardening coefficient (n), the parameter C of the Swift law, the initial sheet thickness (t 0 ), the coefficient of friction (µ 0 ) and the blank holder force (BHF).The uncertainty in the input parameters, associated with the material properties, friction, and process conditions, is assumed to follow a normal distribution with a given mean (µ) and standard deviation (σ), as shown in Table 1.The mean values of the constitutive parameters, blank thickness and friction coefficient were obtained from [16] based on the NUMISHEET' 93 Benchmark.The mean value of the blank holder force was optimized in a previous work [17].The standard deviations were obtained from [16] based on a literature review, except the Poisson ratio, friction coefficient and blank holder force that were based on empirical assumptions made in [17].The influence of the uncertainty in these input parameters was studied on the numerical results of the square cup, i.e., the output parameters of the model: punch force (PF), equivalent plastic strain (EPS), thickness reduction (TR), geometry changes (GC) and springback (SB).The PF and the EPS values are directly obtained from the numerical simulation, while the TR, the GC and SB are defined by: where t 0 and t f are the initial and final sheet thickness, respectively, evaluated for each node of the square cup numerical model; (x f , y f , z f ) and (x f , y f , z f ) are, respectively, the final spatial position of a given node for the numerical simulation with and without uncertainty (i.e., using the mean values of Table 1); and (x bs , y bs , z bs ) is the spatial position of a given node before the springback, i.e., before removing the tools.The GC quantifies the positional difference of a given node between the deterministic and the stochastic simulation, i.e., it is a measure of the shape accuracy between the desired cup shape and the final shape that is affected by the uncertainty.The SB quantifies the positional difference of a given node before and after removing the tools, i.e., it is a measure of the springback.The variability in the forming results was evaluated using a quasi-Monte Carlo method coupled with Sobol sequence to generate the sample of input parameters according to the normal distributions given in Table 1.A sample size of 3000 simulations of the square cup was chosen to ensure the convergence of the statistical measures, mean and standard deviation.Figure 2 shows the mean and standard deviation for the results of EPS, TR, Metals 2024, 14, 432 5 of 18 GC and SB along the square cup [23].From this figure, it can be seen that the equivalent plastic strain reached its maximum at the wall and side edges of the flange, with the most pronounced dispersion of its values occurring in the latter region.In the case of thickness reduction, the highest mean values occurred at the cup wall, with significant standard deviation values only at the flange edge.The mean and standard deviation value for the geometry changes were higher at the flange edge and wall.The highest mean and standard deviation values for the springback occurred at the flange edge.The bottom of the cup is the only area where the uncertainty of the input parameters did not affect the cup results, making a sensitivity analysis in this area unnecessary.Figure 3 shows the evolution of the mean and standard deviation of the force applied by the punch as a function of its displacement.The force reaches maximum values between 17 and 30 mm of displacement and the standard deviation was higher within this range of values, indicating that it is during this phase of displacement that the uncertainty of the input parameters has the greatest effect on the force variability of the punch.Figure 4 shows the distributions of the maximum values of each output parameter (i.e., forming results).All outputs follow central tendency distributions, except for the springback (Figure 4e).It is important to consider the distributions of each output variable because, as noted above, the distributions of the results limit the sensitivity analysis that can be applied.In particular, the distribution associated with springback (Figure 4e) is an obstacle to the use of variance-based approaches (e.g., Sobol indices and factorial ANOVA), as these are not appropriate for data that do not follow a distribution with central tendency.Figure 4 shows the distributions of the maximum values of each output parameter (i.e., forming results).All outputs follow central tendency distributions, except for the springback (Figure 4e).It is important to consider the distributions of each output variable because, as noted above, the distributions of the results limit the sensitivity analysis that can be applied.In particular, the distribution associated with springback (Figure 4e) is an obstacle to the use of variance-based approaches (e.g., Sobol indices and factorial ANOVA), as these are not appropriate for data that do not follow a distribution with central tendency.

Sensitivity Analysis
The PAWN indices are a sensitivity analysis technique based on the cumulative distribution function of a given data sample, as opposed to other techniques (Sobol and factorial ANOVA) based on the variance of the probability density functions [20].This allows PAWN indices to be applied to samples that follow more complex distributions, i.e., multimodal distributions.The basic concept behind these indices is that the influence of an

Sensitivity Analysis
The PAWN indices are a sensitivity analysis technique based on the cumulative distribution function of a given data sample, as opposed to other techniques (Sobol and factorial ANOVA) based on the variance of the probability density functions [20].This allows PAWN indices to be applied to samples that follow more complex distributions, i.e., multimodal distributions.The basic concept behind these indices is that the influence of an input is related to the change produced in the cumulative distribution function of an output.
The first approach proposed for deriving PAWN indices, for a given output y, is to compare cumulative distribution functions [20].Specifically, the cumulative distribution function obtained by assuming that all input parameters vary, F y (y), is compared with other n functions, F y|x i (y), obtained by assuming a fixed value for a particular input parameter, x i , while varying the remaining input parameters.The n functions F y|x i (y) are derived for different fixed values of x i .The PAWN index, P i , for the parameter x i was obtained through a statistical measure, such as the average of the Kolmogorov-Smirnov statistics between the function F y (y) and the n functions F y|x i (y) with the input parameter x i fixed.The PAWN index, P i , is defined by the following equation: where KS(x i ) is the Kolmogorov-Smirnov statistic that quantifies the difference between the cumulative distribution functions F y (y) and one of the n cumulative distribution functions F y|x i (y) when the input parameter x i is fixed at a given value x i (n) .The statistic KS(x i ) is obtained by the following expression: Subsequently, a more efficient approach [21] was proposed that allows the use of generic sampling (i.e., it is not mandatory to assume fixed values for the input parameter x i ).The new approach allows PAWN indices to be calculated, starting by dividing the input parameter domain x i into n equal intervals and obtaining the cumulative distribution functions for each of the n subdomains, F y|x i (y).In addition, the cumulative distribution, F y (y), was calculated for the entire data set, i.e., the entire domain of x i .Similarly to the first approach, the maximum vertical difference was calculated (using the Kolmogorov-Smirnov statistic) between the cumulative distribution curve of each data interval, F y|x i (y), and the cumulative distribution curve of the entire data set, F y (y).The value of the PAWN index, P i , for the parameter x i was calculated using Equation (9).
The Sobol indices are a sensitivity analysis technique based on the calculation of the variance (i.e., based on the probability distribution function), unlike the PAWN indices which are based on the calculation of the cumulative distribution function.The variance is a parameter of central tendency that represents the average of the squares of the difference between the value of the mean of a given distribution and the value of each of its constituent points.The Sobol Indices (total sensitivity indices) were calculated using the following equation [29]: where V(y) is the unconditional variance of the output parameter y and V[E(y|x ∼i )] is the conditional variance of the expected value of the output parameter y when varying all the input parameters except x ∼i .The fact that Sobol indices are based on variance means that they only have accurate results for unimodal distributions.Bearing this in mind, it is expected that the results obtained by the PAWN indices for the springback output parameter will be different from the results obtained by the Sobol indices for the same parameter.On the other hand, PAWN indices have the disadvantage of requiring the calculation of cumulative distribution functions, whereas Sobol indices only require the calculation of variance.Another disadvantage of PAWN indices is that it is necessary to define the number of intervals into which the data are divided, which is another factor that can influence the results.Both indices quantify the influence of the parameter uncertainty on the forming results in the same way; i.e., the higher their values are, the higher is the influence of a given input uncertainty on the result variability.For indices close to zero there is no influence on the result.

Stabilization Analysis
The indices are considered to be stabilized when their variation is equal to or less than 5%.This stabilization criterion is the same for all input parameters and is validated for successive intervals of 500 simulations, i.e., for 500, 1000, 1500, 2000, 2500 and 3000 simulations.The stabilization analysis of the indices was calculated considering the maximum values of the output parameters.From Figure 5, in the case of punch force, equivalent plastic strain and thickness reduction, the value of the PAWN and Sobol indices stabilized with the use of 2500 and 1500 base simulations, respectively.For the geometry change, the PAWN and Sobol indices stabilized their value when more than 2000 simulations were used.For the springback, the PAWN and Sobol indices stabilized when more than 1500 simulations were used.It is worth noting that in the case of the Sobol indices, a base sample of m simulations requires a total of m × (11+ 2) simulations to evaluate the sensitivity indices for the 11 input parameters, according to the traditional procedure [30].It is therefore evident that the computation of the PAWN indices is clearly more efficient.
Figure 6 shows the results of the PAWN and Sobol indices for the maximum values of the output parameters.The figure shows that the uncertainty in the C parameter of Swift's law has the greatest influence on the variability of the punch force (Figure 6a).The variability in the equivalent plastic strain is most influenced by the uncertainty in the hardening coefficient n, and in the anisotropy coefficient, r 90 (Figure 6b).The variability in the thickness reduction is most affected by the hardening coefficient n, and the anisotropy coefficient, r 90 (Figure 6c).As for the variability in the geometry changes (Figure 6d), the anisotropy coefficients r 90 and r 0 are the ones that have greater importance for the PAWN indices.For the Sobol indices, r 90 and r 45 are the most important inputs.In the springback (Figure 6e), the results are similar in terms of which parameters are more influential, although there is a significant discrepancy in the values of the indices.The most influential parameters in this case are the hardening coefficient n, the parameter C and the blank holder force.The uncertainties that most influence the variability in the punch force, equivalent plastic strain and thickness reduction, are similar for the PAWN and Sobol indices.However, for the geometry change and springback, the most influential uncertainties are different for the PAWN and Sobol indices.

Pareto Analysis
In this section, a Pareto analysis is performed on the results of the indices to separate the more important from the less important input parameters.Figures 7 and 8 show the graphs corresponding to the Pareto analysis for both sensitivity analysis techniques, where the vertical bars represent the absolute value of the index for the parameter in question and the black curve represents the accumulated percentage.According to Pareto, the most important parameters are those that contribute to 80% of the variability of the output parameters.From Figures 7 and 8, it can be observed that the results of the Sobol indices consider that a small number of uncertainty factors (at most six input parameters) are responsible for 80% of the output variability.On the other hand, the PAWN indices consider that a large number of uncertainty factors (at least six input parameters) are responsible for 80% of the output variability.the springback (Figure 6e), the results are similar in terms of which parameters are more influential, although there is a significant discrepancy in the values of the indices.The most influential parameters in this case are the hardening coefficient , the parameter  and the blank holder force.The uncertainties that most influence the variability in the punch force, equivalent plastic strain and thickness reduction, are similar for the PAWN and Sobol indices.However, for the geometry change and springback, the most influential uncertainties are different for the PAWN and Sobol indices.

Pareto Analysis
In this section, a Pareto analysis is performed on the results of the indices to separate the more important from the less important input parameters.Figures 7 and 8 show the graphs corresponding to the Pareto analysis for both sensitivity analysis techniques,    From Figures 7 and 8, it can be concluded that for the punch force, according to the Sobol indices, 80% of the variability is associated with the uncertainties in the  parameter of Swift's law and the initial thickness (Figure 8a).In addition to these parameters, the PAWN indices also identify the uncertainties in the hardening coefficient, the anisotropy coefficients at 45° and 90°, and the Poisson ratio as being responsible.As for the equivalent plastic strain, the Sobol indices (Figure 8b), 80% of the variability comes from the uncertainties in the hardening coefficient and the anisotropy coefficient at 90°.In addition to these parameters, the PAWN indices (Figure 7b) show that 80% of the variability in the equivalent plastic strain is caused by the uncertainty in the anisotropy coefficients at 0° and 45°, the initial yield stress and the initial blank thickness.
The variability in the thickness reduction is significantly affected by the uncertainty in the hardening coefficient, the anisotropy coefficients at 0°, 45° and 90°, the  and the initial thickness, according to the PAWN indices (Figure 7c).For the Sobol indices (Figure 8c), only the uncertainties of the hardening coefficient and the anisotropy coefficients at 45° and 90° account for 80% of the thickness reduction variability.For the PAWN indices, the variability in the geometry changes is due to the uncertainties in the anisotropy coefficients at 0°, 45° and 90° and the hardening coefficient, the Poisson ratio and the  (Figure 7d).In the case of the Sobol indices (Figure 8d), the uncertainty in the anisotropy coefficients and the hardening coefficient are the most important.For the springback, both indices (Figures 7e and 8e) assume that 80% of the variability is due to the uncertainty in the parameter  of Swift's law, the  , the hardening coefficient, the anisotropy From Figures 7 and 8, it can be concluded that for the punch force, according to the Sobol indices, 80% of the variability is associated with the uncertainties in the C parameter of Swift's law and the initial thickness (Figure 8a).In addition to these parameters, the PAWN indices also identify the uncertainties in the hardening coefficient, the anisotropy coefficients at 45 • and 90 • , and the Poisson ratio as being responsible.As for the equivalent plastic strain, the Sobol indices (Figure 8b), 80% of the variability comes from the uncertainties in the hardening coefficient and the anisotropy coefficient at 90 • .In addition to these parameters, the PAWN indices (Figure 7b) show that 80% of the variability in the equivalent plastic strain is caused by the uncertainty in the anisotropy coefficients at 0 • and 45 • , the initial yield stress and the initial blank thickness.
The variability in the thickness reduction is significantly affected by the uncertainty in the hardening coefficient, the anisotropy coefficients at 0 • , 45 • and 90 • , the BHF and the initial thickness, according to the PAWN indices (Figure 7c).For the Sobol indices (Figure 8c), only the uncertainties of the hardening coefficient and the anisotropy coefficients at 45 • and 90 • account for 80% of the thickness reduction variability.For the PAWN indices, the variability in the geometry changes is due to the uncertainties in the anisotropy coefficients at 0 • , 45 • and 90 • and the hardening coefficient, the Poisson ratio and the BHF (Figure 7d).In the case of the Sobol indices (Figure 8d), the uncertainty in the anisotropy coefficients and the hardening coefficient are the most important.For the springback, both indices (Figures 7e and 8e) assume that 80% of the variability is due to the uncertainty in the parameter C of Swift's law, the BHF, the hardening coefficient, the anisotropy coefficients at 90 • , and the Young's modulus (for PAWN indices) or the initial yield stress (for Sobol indices).In summary, the main uncertainty factors are similar for both sensitivity analyses, although the PAWN sensitivity shows that the results are influenced by more factors.

Sensitivity Indices per Region of the Square Cup
Sensitivity indices are calculated for the different regions of the square cup.The aim of this analysis is to understand, for the different output parameters, which uncertainties in the input parameters are most influential in different regions of the cup.For the equivalent plastic strain, the uncertainty in the hardening coefficient (Figure 9a) has a significant influence in the central region of the wall and at the radius of curvature of the punch.Since the maximum variability of the equivalent plastic strain is observed in the central region of the wall (see Figure 2a), it is therefore essential to control the uncertainty in the hardening coefficient.The uncertainty in the anisotropy coefficient at 90 • significantly influences the equivalent plastic strain variability in the radius of curvature of the die and in the upper region of the cup wall, which is also a critical region with a high variability in the equivalent plastic strain values (see Figure 2a).The results obtained are qualitatively similar for both indices.
Metals 2024, 14, x FOR PEER REVIEW 14 of 21 coefficients at 90°, and the Young's modulus (for PAWN indices) or the initial yield stress (for Sobol indices).In summary, the main uncertainty factors are similar for both sensitivity analyses, although the PAWN sensitivity shows that the results are influenced by more factors.

Sensitivity Indices per Region of the Square Cup
Sensitivity indices are calculated for the different regions of the square cup.The aim of this analysis is to understand, for the different output parameters, which uncertainties in the input parameters are most influential in different regions of the cup.For the equivalent plastic strain, the uncertainty in the hardening coefficient (Figure 9a) has a significant influence in the central region of the wall and at the radius of curvature of the punch.Since the maximum variability of the equivalent plastic strain is observed in the central region of the wall (see Figure 2a), it is therefore essential to control the uncertainty in the hardening coefficient.The uncertainty in the anisotropy coefficient at 90° significantly influences the equivalent plastic strain variability in the radius of curvature of the die and in the upper region of the cup wall, which is also a critical region with a high variability in the equivalent plastic strain values (see Figure 2a).The results obtained are qualitatively similar for both indices.From Figure 10, it can be concluded that the variability in the thickness reduction, at the bo om of the cup, along the radius of curvature of the punch and in the central region of the wall (Figure 10a) is significantly influenced by the uncertainty in the strain hardening coefficient.Although the variability at the bo om of the cup is not significant (see Figure 2b), along the radius of curvature of the punch and in the central region of the wall there is a significant variability in the thickness reduction (see Figure 2b), so the uncertainty in the hardening coefficient is essential to control.The uncertainty in the anisotropy coefficient at 90° (Figure 10b) significantly affects the variability of the thickness reduction in the cup wall, the radius of curvature of the die and the flange.However, only the wall zone has a significant variability in thickness reduction (see Figure 2b), which means that the uncertainty of the anisotropy coefficient at 90° is particularly critical in this region.The results obtained are qualitatively similar for both indices.From Figure 10, it can be concluded that the variability in the thickness reduction, at the bottom of the cup, along the radius of curvature of the punch and in the central region of the wall (Figure 10a) is significantly influenced by the uncertainty in the strain hardening coefficient.Although the variability at the bottom of the cup is not significant (see Figure 2b), along the radius of curvature of the punch and in the central region of the wall there is a significant variability in the thickness reduction (see Figure 2b), so the uncertainty in the hardening coefficient is essential to control.The uncertainty in the anisotropy coefficient at 90 • (Figure 10b) significantly affects the variability of the thickness reduction in the cup wall, the radius of curvature of the die and the flange.However, only the wall zone has a significant variability in thickness reduction (see Figure 2b), which means that the uncertainty of the anisotropy coefficient at 90 • is particularly critical in this region.The results obtained are qualitatively similar for both indices.Figure 11a shows that the variability in the geometry changes is influenced by the hardening coefficient at the bo om of the wall and at the bo om of the cup.The uncertainty in the anisotropy coefficients at 0° and 90° (Figure 11b,c) largely influence the variability of the geometry changes in the upper part of the square cup.These areas coincide with the regions of the cup where the variability in the geometry change is significant (see Figure 2c).Therefore, to control the variability of the geometry change, it is more important to control the uncertainty in the anisotropy coefficients at 0° and 90° than in the hardening coefficient.The results obtained were qualitatively similar for both indices.Figure 11a shows that the variability in the geometry changes is influenced by the hardening coefficient at the bottom of the wall and at the bottom of the cup.The uncertainty in the anisotropy coefficients at 0 • and 90 • (Figure 11b,c) largely influence the variability of the geometry changes in the upper part of the square cup.These areas coincide with the regions of the cup where the variability in the geometry change is significant (see Figure 2c).Therefore, to control the variability of the geometry change, it is more important to control the uncertainty in the anisotropy coefficients at 0 • and 90 • than in the hardening coefficient.The results obtained were qualitatively similar for both indices.The springback variability throughout the cup is mainly influenced by the uncertainty in the blank holder force, , the anisotropy coefficient at 90° and the  parameter of Swift's law (Figure 12a).However, it can be noted that among these parameters, the parameter  of Swift's law is the most critical (Figure 12b), since it has a greater effect in areas where the springback variability reaches higher values.For springback, it is clear The springback variability throughout the cup is mainly influenced by the uncertainty in the blank holder force, BHF, the anisotropy coefficient at 90 • and the C parameter of Swift's law (Figure 12a).However, it can be noted that among these parameters, the parameter C of Swift's law is the most critical (Figure 12b), since it has a greater effect in areas where the springback variability reaches higher values.For springback, it is clear that the distributions of the Sobol and PAWN indices are different.This is to be expected as the springback distribution (see Figure 4) is bimodal and therefore a sensitivity analysis using Sobol indices is not appropriate in this case.

Maximum Sensitivity Indices per Region of the Square Cup
Figure 13 shows the input parameters for which the PAWN or Sobol indices are maximum in each zone of the square cup.For the equivalent plastic strain (Figure 13a), the geometry changes (Figure 13b), the thickness reduction (Figure 13c), the uncertainty in the anisotropy coefficient at 90° and the hardening coefficient were the most important in a large part of the cup; the uncertainty in the anisotropy coefficient at 0° also presented some areas of significant importance.The springback variability throughout the cup was significantly affected by the uncertainty in the blank holder force and the parameter , with the uncertainty in the anisotropy coefficients at 45° and 90° also influential in certain locations of the cup (Figure 13d).As expected, the results for the PAWN and Sobol indices were identical for all output parameters except springback.The main difference was at the bo om and wall of the cup, where the PAWN indices indicate that the uncertainty in the anisotropy coefficient at 45° is the most important parameter, while the Sobol indices consider the  to be the most important parameter.

Maximum Sensitivity Indices per Region of the Square Cup
Figure 13 shows the input parameters for which the PAWN or Sobol indices are maximum in each zone of the square cup.For the equivalent plastic strain (Figure 13a), the geometry changes (Figure 13b), the thickness reduction (Figure 13c), the uncertainty in the anisotropy coefficient at 90 • and the hardening coefficient were the most important in a large part of the cup; the uncertainty in the anisotropy coefficient at 0 • also presented some areas of significant importance.The springback variability throughout the cup was significantly affected by the uncertainty in the blank holder force and the parameter C, with the uncertainty in the anisotropy coefficients at 45 • and 90 • also influential in certain locations of the cup (Figure 13d).As expected, the results for the PAWN and Sobol indices were identical for all output parameters except springback.The main difference was at the bottom and wall of the cup, where the PAWN indices indicate that the uncertainty in the anisotropy coefficient at 45 • is the most important parameter, while the Sobol indices consider the BHF to be the most important parameter.

Evolution of Sensitivity Indices for the Punch Force
Figure 14 shows the evolution of the PAWN and Sobol indices for the punch force as a function of its displacement.It can be seen that the influence of the uncertainty in the parameter  on the punch force variability increases with the punch displacement.The influence of the uncertainties in the hardening coefficient, , and the yield stress,  , on the punch force variability increased initially but lost influence with the punch

Evolution of Sensitivity Indices for the Punch Force
Figure 14 shows the evolution of the PAWN and Sobol indices for the punch force as a function of its displacement.It can be seen that the influence of the uncertainty in the parameter C on the punch force variability increases with the punch displacement.The influence of the uncertainties in the hardening coefficient, n, and the yield stress, Y 0 , on the punch force variability increased initially but lost influence with the punch displacement.The uncertainty in the initial thickness has a constant influence on the punch force throughout its displacement.The primary distinction between the results obtained from the two-sensitivity analysis is the influence of the uncertainty in the anisotropy coefficient, r 0 .However, the r 0 significant variations between 0 and 12 mm (for the Sobol indices) is artificial, i.e., it comes from errors in the calculation of Equation ( 12), when both variances in the dividend and divisor were approximately 0. displacement.The uncertainty in the initial thickness has a constant influence on the punch force throughout its displacement.The primary distinction between the results obtained from the two-sensitivity analysis is the influence of the uncertainty in the anisotropy coefficient,  .However, the  significant variations between 0 and 12 mm (for the Sobol indices) is artificial, i.e., it comes from errors in the calculation of Equation ( 12), when both variances in the dividend and divisor were approximately 0.

Conclusions
The main objective of this study was to quantify the influence of the uncertainty in material, friction, and process parameters on the forming results of a square cup.For this, a sensitivity analysis was carried out using two different techniques.It was concluded that the bo om of the cup is the area least affected by the uncertainty in the material, friction and process parameters.On the other hand, the flange, the radius of curvature of the die and the wall are the areas that show the greatest variability.For the equivalent plastic strain, the variability is more pronounced at the flange, as for geometry changes and springback.The uncertainty in the parameter  of Swift's law, the hardening coefficient and the initial sheet thickness significantly influence the variability of the punch force.The variability of equivalent plastic strain, thickness reduction and geometry changes are significantly influenced by the uncertainty in the hardening coefficient and the anisotropy coefficients.The variability of springback is significantly influenced by the uncertainty of parameter  of Swift's law, the hardening coefficient, the anisotropy coefficient at 90° and the blank holder force.It is worth noting that there is a difference between the results obtained using PAWN and Sobol indices, particularly in the case of the springback response.The study emphasizes that the type of distribution of the results plays a crucial

Conclusions
The main objective of this study was to quantify the influence of the uncertainty in material, friction, and process parameters on the forming results of a square cup.For this, a sensitivity analysis was carried out using two different techniques.It was concluded that the bottom of the cup is the area least affected by the uncertainty in the material, friction and process parameters.On the other hand, the flange, the radius of curvature of the die and the wall are the areas that show the greatest variability.For the equivalent plastic strain, the variability is more pronounced at the flange, as for geometry changes and springback.The uncertainty in the parameter C of Swift's law, the hardening coefficient and the initial sheet thickness significantly influence the variability of the punch force.The variability of equivalent plastic strain, thickness reduction and geometry changes are significantly influenced by the uncertainty in the hardening coefficient and the anisotropy coefficients.The variability of springback is significantly influenced by the uncertainty of parameter C of Swift's law, the hardening coefficient, the anisotropy coefficient at 90 • and the blank holder force.It is worth noting that there is a difference between the results obtained using PAWN and Sobol indices, particularly in the case of the springback response.The study emphasizes that the type of distribution of the results plays a crucial role in the quality of the results, and it is therefore crucial to select the most appropriate sensitivity analysis technique.For example, springback has a bimodal distribution, making PAWN indices more suitable for sensitivity analysis.The PAWN sensitivity analysis proved to be more computationally efficient than the Sobol sensitivity analysis, with improved accuracy for results following multimodal distributions.Although the PAWN indices are more efficient, evaluating them still requires a significant number of numerical simulations.To reduce the computational cost without compromising the quality of the results, we propose exploring the use of metamodelling techniques to evaluate both indices.This will allow to perform sensitivity analysis on more complex numerical models, that for instance consider deformable tools or the use of more advanced constitutive models that can describe texture evolution and the Bauschinger effect.

Figure 1 .
Figure 1.Square cup forming process: (a) dimensions of the tools in mm; (b) numerical model.adapted with permission from [18]; (c) mesh sensitivity.

Figure 1 .
Figure 1.Square cup forming process: (a) dimensions of the tools in mm; (b) numerical model.adapted with permission from [18]; (c) mesh sensitivity.

Figure 3 .
Figure 3. Evolution of the mean (a) and standard deviation (b) of the force applied by the punch.

21 Figure 3 .
Figure 3. Evolution of the mean (a) and standard deviation (b) of the force applied by the punch.

Figure 13 .
Figure 13.Parameters corresponding to the maximum PAWN index (left) and Sobol index (right) in each zone of for (a) EPS; (b) GC; (c) TR; and (d) SB.

Figure 13 .
Figure 13.Parameters corresponding to the maximum PAWN index (left) and Sobol index (right) in each zone of for (a) EPS; (b) GC; (c) TR; and (d) SB.

Figure 14 .
Figure 14.Indices of punch force as a function of displacement for (a) PAWN indices and (b) Sobol indices.

Figure 14 .
Figure 14.Indices of punch force as a function of displacement for (a) PAWN indices and (b) Sobol indices.