Earing Reduction by Varying Blank Holding Force in Deep Drawing with Deep Neural Network

In the present study, we propose a novel method of varying blank holding force (BHF) with the segmental blank holder and investigated its influence on the earing reduction in the circular deep drawing process of an aluminum alloy sheet. Based on the analysis of cup height profile, the principle of varying BHF using segmental blank holder was presented and analyzed by analytical theory and numerical simulation. The optimal varying BHF was reasonably determined and compared by using the analytical model and deep neural network (DNN) model integrated with genetic algorithm (GA). The integrated DNN-GA model revealed an accurate prediction and optimization of varying BHF for the minimum earing height variation, which showed a superior result to the analytical model. The optimal varying BHF exhibited a significant influence on the earing formation, resulting in the noticeable decrease of earing height variation. For volume consistency, it was found that an increase in thickness at the cup wall region predicted with the optimal varying BHF was achieved in the transverse direction, which implies an improvement of deep-drawability. Such results indicate that the varying BHF is more reasonable and effective than the uniform BHF. Furthermore, the material properties of the blank sheet also affected the reduction of earing in the deep drawing with varying BHF. The present study revealed that the lower the material strength, the more significant the earing reduction in the deep drawing with varying BHF will be.


Introduction
Deep drawing is a widely used sheet metal forming process in the automotive, packaging, and household appliances industries. During the deep drawing process, the punch slowly pushes the blank sheet down into the die, and the material of the blank sheet flows into the die cavity. This forming process step is characterized by the material plastic-flow phase, which may be affected by blank holding force, drawbeads, type, and amount of lubricant as well as shape and size of the initial blank. Hence, it is essential to control the plastic-flow phase of the sheet material into the die cavity, in order to obtain the final component shape without any defects.
One of the major defects that commonly occurs in the circular cup deep drawing process is called earing. Earing is defined as the formation of wavy edges at the open end of the cup during deep drawing due to the anisotropy of the sheet metal. Comprehensive studies have been conducted to investigate earing defects in deep drawing [1][2][3][4]. In the industry, earing is undesirable since an additional processing step is required to trim some metal from the top of the deep-drawn cup. This wastes material, is time-consuming, and has a high manufacturing cost. In the recent past, various different approaches have been proposed to eliminate or minimize earing defects in the deep drawing process [5][6][7][8][9][10][11]. The optimization of the initial blank shape is the most widely used method that can effectively minimize the formation of earing. Kishor and Ravi Kumar [5] reported that earing in deep-drawn cups was significantly reduced by using a modified initial blank shape with consideration of the material flow in different directions during deep drawing and planar anisotropy of the sheet metal. Recently, Singh et al. [6] proposed a two-stage blank modification technique combining directionalities of both the yield strength and plastic strain ratio to minimize earing defects. They presented a noticeable reduction in earing height by 83% with the modified blank. Aretz and Engler [7] studied the analytical earing compensation approach to produce new blank geometries, resulting in a deep-drawn cup without ears. In other approaches, investigations of various texture components [9], different processing routes [10], and the use of drawbeads [11] have also been proposed to reduce the formation of earing in the deep drawing process.
In the traditional deep drawing process, a uniform blank holding force (BHF), which that is kept constant throughout the process, is generally used for shallow drawing of parts with relatively simple shapes. For complex parts, uniform BHF could lead to poor control of the material plastic-flow phase and low accuracy of deep drawing [12]. Recently, variable BHF has received more attention as a sheet metal forming technique used in advanced manufacturing technologies. A number of comprehensive studies have numerically and experimentally investigated variable BHF to improve the formability as well as the wrinkling and tearing defects in the deep drawing process [12][13][14]. To date, however, there have been few studies on the use of variable BHF to address earing defects in deep drawing. Previous studies [15,16] have reported that the earing height increased with the increase of BHF. In order to attenuate such a plastic anisotropy, in the present study, a non-uniform BHF distribution is proposed so that the locally varying BHF is reduced and increased in the hill and valley regions, respectively, by using the blank holder segments for earing reduction. Zhang et al. [17] presented a novel radial segmental blank holder technique based on electro-permanent magnet technology to improve the formability of the deep drawing process. Hassan et al. [18] developed and fabricated a new blank holder divided into eight segments consisting of four drawing segments and four small wedges to investigate the effect of aided friction on the deep drawing of the circular cup. They found that the earing height obtained by the aided friction using the eight-segment blank holder was in good agreement with the theoretical curve. Han et al. [19] reported that an earless circular cup could be achieved by applying a distributed BHF. However, they did not present the method for the determination of distributed BHF and there was no discussion on the thickness variation in a deep-drawn cup when the distributed BHF was used. In the present study, therefore, an analytical model based on the analysis of cup height was proposed to determine the optimization of varying BHF, leading to the greatest reduction in earing.
Recently, finite element (FE) simulation has been demonstrated as a robust and effective tool in the analysis of the metal forming process. Kim et al. [20] revealed the influence of the individual orientation crystals and the interaction with neighboring crystals on the rotations of crystals in the polycrystalline interstitial free (IF) steel during channel die compression, based on the fully three-dimensional crystal plasticity finite element method (CPFEM). Neto et al. [21] investigated the effect of boundary conditions on the prediction of springback and wrinkling in sheet metal forming. They reported that the symmetry conditions and the blank material significantly influenced the springback and wrinkling behaviors, and the FE analysis results had good consistency with the experimental data. In recent years, the FE simulation has also been widely used in studies of deep drawing, for example, for evaluation of a novel flexible technique for a micro deep drawing process utilizing a floating ring [22], investigation of a warm hydromechanical deep drawing process to improve the formability of lightweight materials [23], and the prediction of wrinkling [24,25], especially the prediction of earing in the deep drawing process [26][27][28]. Shi et al. [26] developed a crystal plasticity-based finite element model to simulate the earing in deep drawing. They found that the earing formed early at a very low strain during deep drawing and intensified with the increased drawing, but retained the overall earing profiles. Tikhovskiy et al. [27] showed that the prediction of earing in the deep drawing of Nowadays, artificial neural network (ANN) has been employed as a flexible and efficient model to express the non-linear and complex relationships between input and output parameters [29][30][31]. Wang et al. [30] successfully developed the ANN model to predict the mechanical properties of austenitic stainless steel. They found that the yield strength and ultimate tensile strength were accurately predicted as a function of chemical composition, heat treatment conditions, and test temperature. Furthermore, genetic algorithm (GA) has been widely used and combined with the ANN model in solving optimization problems in various research fields [32][33][34][35]. This hybrid approach is considered as a powerful and empirical technique to predict and optimize the input variables with respect to the maximum or minimum objectives.
The present study investigated the effect of varying BHF by using a segmental blank holder on the earing formation in the circular deep drawing process of an aluminum alloy sheet. An analytical model was proposed to determine the optimized varying BHF. In another approach, the deep neural network (DNN) model coupled with GA is well developed to predict and optimize the varying BHF for minimization of earing height variation. Then, the prediction capability of the integrated DNN-GA model was compared with the proposed analytical model in this work. The earing height profiles and thickness variations along the rolling direction (RD), diagonal direction (DD), and transverse direction (TD) of the deep-drawn cup by using uniform and varying BHFs are discussed in detail. Furthermore, due to the predominant dependence of the earing on the material properties, there are limitations on further reductions of earing formation in deep drawing. Therefore, the effect of the material properties on the earing reduction in deep drawing with varying BHF using a segmental blank holder was also investigated in this study.

Design of Segmental Blank Holder
A schematic of the axis-symmetric deep drawing process is shown in Figure 1. Note that R p and r p are the radius and corner radius of the punch, respectively. R d and r d are the radius and corner radius of the die, respectively. R b and t are the radius and thickness of the blank sheet, respectively. the overall earing profiles. Tikhovskiy et al. [27] showed that the prediction of earing in the deep drawing of a ferritic stainless steel sheet was in good agreement with the experimental measurements for the position and the height of the earing, based on the texture component CPFEM. Nowadays, artificial neural network (ANN) has been employed as a flexible and efficient model to express the non-linear and complex relationships between input and output parameters [29][30][31]. Wang et al. [30] successfully developed the ANN model to predict the mechanical properties of austenitic stainless steel. They found that the yield strength and ultimate tensile strength were accurately predicted as a function of chemical composition, heat treatment conditions, and test temperature. Furthermore, genetic algorithm (GA) has been widely used and combined with the ANN model in solving optimization problems in various research fields [32][33][34][35]. This hybrid approach is considered as a powerful and empirical technique to predict and optimize the input variables with respect to the maximum or minimum objectives.
The present study investigated the effect of varying BHF by using a segmental blank holder on the earing formation in the circular deep drawing process of an aluminum alloy sheet. An analytical model was proposed to determine the optimized varying BHF. In another approach, the deep neural network (DNN) model coupled with GA is well developed to predict and optimize the varying BHF for minimization of earing height variation. Then, the prediction capability of the integrated DNN-GA model was compared with the proposed analytical model in this work. The earing height profiles and thickness variations along the rolling direction (RD), diagonal direction (DD), and transverse direction (TD) of the deep-drawn cup by using uniform and varying BHFs are discussed in detail. Furthermore, due to the predominant dependence of the earing on the material properties, there are limitations on further reductions of earing formation in deep drawing. Therefore, the effect of the material properties on the earing reduction in deep drawing with varying BHF using a segmental blank holder was also investigated in this study.

Design of Segmental Blank Holder
A schematic of the axis-symmetric deep drawing process is shown in Figure 1. Note that Rp and rp are the radius and corner radius of the punch, respectively. Rd and rd are the radius and corner radius of the die, respectively. Rb and t are the radius and thickness of the blank sheet, respectively.  In fact, the BHF provided by the conventional blank holder is primarily applied in a narrow region of the flange edge due to the increase in edge thickness during the deep drawing process [17]. In addition, because of the unsteady plastic-flow phase of the material, the BHF may not be effectively applied to the different flange regions of the sheet metal. This means that the sheet material in distinct flange regions would flow differently into the die cavity caused by the anisotropy of the sheet metal, resulting in the occurrence of earing. In order to overcome these drawbacks, a segmental blank holder was proposed to reasonably apply the varied BHF region by region on the flange of the sheet metal to make the material plastic-flow phase uniform in the die cavity, leading to the reduction of earing formation. For the sake of earing reduction, the greater the number of blank holder segments used, the more control over the earing reduction can be obtained [19]. In the present study, due to the orthogonal symmetry, only a quarter blank holder was generated to divide into six segments for introducing the varying BHF. The blank holder segments were determined and identified by the sensitivity analysis, in which the effect of the number of segments on the reduction of earing or improvement of formability was investigated. For the complex parts, the sensitivity analysis should be carried out to determine the appropriate number of segments. The schematic of the blank holder segment is shown in Figure 2, in which the varying BHF can be applied effectively on different flange regions. For this objective, the quarter blank holder consists of six segments with an accurate angle of 15 • . Corresponding to each segment, the applied BHF is needed to properly determine the configuration that results in the greatest reduction of earing formation. In fact, the BHF provided by the conventional blank holder is primarily applied in a narrow region of the flange edge due to the increase in edge thickness during the deep drawing process [17]. In addition, because of the unsteady plastic-flow phase of the material, the BHF may not be effectively applied to the different flange regions of the sheet metal. This means that the sheet material in distinct flange regions would flow differently into the die cavity caused by the anisotropy of the sheet metal, resulting in the occurrence of earing. In order to overcome these drawbacks, a segmental blank holder was proposed to reasonably apply the varied BHF region by region on the flange of the sheet metal to make the material plastic-flow phase uniform in the die cavity, leading to the reduction of earing formation. For the sake of earing reduction, the greater the number of blank holder segments used, the more control over the earing reduction can be obtained [19]. In the present study, due to the orthogonal symmetry, only a quarter blank holder was generated to divide into six segments for introducing the varying BHF. The blank holder segments were determined and identified by the sensitivity analysis, in which the effect of the number of segments on the reduction of earing or improvement of formability was investigated. For the complex parts, the sensitivity analysis should be carried out to determine the appropriate number of segments. The schematic of the blank holder segment is shown in Figure 2, in which the varying BHF can be applied effectively on different flange regions. For this objective, the quarter blank holder consists of six segments with an accurate angle of 15°. Corresponding to each segment, the applied BHF is needed to properly determine the configuration that results in the greatest reduction of earing formation.  Figure 3 shows the variation of the BHF corresponding to the variation of earing height. The principle of varying BHF is to reduce the BHF in regions where earing occurs and increase the BHF in the valley regions. For a better understanding, as shown in Figure  3, it can be observed that the earing was formed at 45° and 135° to the RD. Hence, the applied BHF in these regions should be relatively smaller than those at 0° and 90° to the RD. According to the study by Yoon et al. [36], the analytical formula for the prediction of earing can be successfully derived by considering both the r-value and yield stress directionalities simultaneously, as follows:   Figure 3 shows the variation of the BHF corresponding to the variation of earing height. The principle of varying BHF is to reduce the BHF in regions where earing occurs and increase the BHF in the valley regions. For a better understanding, as shown in Figure 3, it can be observed that the earing was formed at 45 • and 135 • to the RD. Hence, the applied BHF in these regions should be relatively smaller than those at 0 • and 90 • to the RD. According to the study by Yoon et al. [36], the analytical formula for the prediction of earing can be successfully derived by considering both the r-value and yield stress directionalities simultaneously, as follows:

Analytical Model for Varying Blank Holding Force (BHF) with Segmental Blank Holder
where H cup (θ) is the cup height and θ is the circumferential angle to the RD. The definition of other parameters in Equation (1) can be found in [36]. In this study, based on the variation of BHF as indicated in Figure 3, the analytical model for the determination of varying BHF applied on the blank holder segment, F(θ), is proposed as where α > 0, β < 0 and γ > 0 are the fitting parameters, while [f ] is expressed as Metals 2021, 11,395 where H cup (θ) is the cup height and θ is the circumferential angle to the RD. The d of other parameters in Equation (1) can be found in [36]. In this study, based on ation of BHF as indicated in Figure 3, the analytical model for the determination of BHF applied on the blank holder segment, F(θ), is proposed as where α > 0, β < 0 and γ > 0 are the fitting parameters, while [f] is expressed as Here, H cup is the cup height as a function of θ + T/2, where T is the period tha obtained from the earing height profile, as shown in Figure 3. The subscripts "avg" and "min" correspond to the average, maximum, and minimum cup heights, resp By substituting Equation (4) Here, H cup is the cup height as a function of θ + T/2, where T is the period that can be obtained from the earing height profile, as shown in Figure 3. The subscripts "avg", "max", and "min" correspond to the average, maximum, and minimum cup heights, respectively. By substituting Equation (4) into Equation (2), the magnitude of BHF for each segment can be determined analytically. According to Equation (4), f > 0 when the H cup is larger than H cup avg , and Equation (2) becomes F (θ) = α[f ] + γ to illustrate the gradual increase of BHF applied above H cup avg . Additionally, Equation (2) becomes F (θ) = β[−f ] + γ in the case of f ≤ 0 when the H cup is smaller than H cup avg , indicating the gradual decrease of BHF provided below H cup avg . Thus, a non-uniform BHF distribution is proposed so that the locally varying BHF is reduced and increased in the hill and valley regions, respectively, in order to induce uniform plastic-flow phase of the material.

Finite Element Simulation of Deep Drawing
Using Hill's 1948 anisotropic material model [37], available in the library of commercial software ABAQUS, a FE simulation of the deep drawing of an aluminum alloy sheet was carried out. Due to the orthogonal symmetry of the tool and blank geometries, only a quarter of the FE model is generated for the deep drawing simulation in order to save the computational time, as shown in Figure 4a. The blank sheet was initially meshed using a shell element (S4R) with the reduced integration. Figure 4b shows the FE mesh of the blank sheet with the total number of elements of 1600. The dimensions of the tool and the blank sheet for the FE simulation of deep drawing are given in Table 1 [36]. The punch, die, and blank holder are considered as rigid bodies whose material properties do not need to be defined. The die is fixed in all directions, while the punch can move in the vertical direction. The blank holder is used to apply force on the flange region of the blank sheet. A penalty contact is used between the tool and blank sheet by the surface-to-surface model. The Coulomb friction coefficient between the tool and the blank sheet is assumed to be 0.1 [38]. To investigate the effect of varying BHF on earing reduction in deep drawing, in the present study, the friction coefficient was assumed to be the same value in different flange regions of the sheet metal.
Using Hill's 1948 anisotropic material model [37], available in the library of c cial software ABAQUS, a FE simulation of the deep drawing of an aluminum all was carried out. Due to the orthogonal symmetry of the tool and blank geometr a quarter of the FE model is generated for the deep drawing simulation in orde the computational time, as shown in Figure 4a. The blank sheet was initially mesh a shell element (S4R) with the reduced integration. Figure 4b shows the FE mes blank sheet with the total number of elements of 1600. The dimensions of the too blank sheet for the FE simulation of deep drawing are given in Table 1 [36]. Th die, and blank holder are considered as rigid bodies whose material propertie need to be defined. The die is fixed in all directions, while the punch can mov vertical direction. The blank holder is used to apply force on the flange region of t sheet. A penalty contact is used between the tool and blank sheet by the surface-to model. The Coulomb friction coefficient between the tool and the blank sheet is a to be 0.1 [38]. To investigate the effect of varying BHF on earing reduction in dee ing, in the present study, the friction coefficient was assumed to be the same valu ferent flange regions of the sheet metal.  To investigate the effectiveness of varying BHF on earing reduction and t variation in deep drawing, an AA1050-O aluminum alloy sheet was used as th sheet material. The material properties of the AA1050-O blank sheet were chos identical to those in Yoon et al. [39], as given in Table 2. Then, several FE simul deep drawing were carried out to determine the conditions under which a prop form BHF with less earing can be formed. An additional FE simulation was a formed using varying BHF to observe the reduction of earing formation in deep d   To investigate the effectiveness of varying BHF on earing reduction and thickness variation in deep drawing, an AA1050-O aluminum alloy sheet was used as the blank sheet material. The material properties of the AA1050-O blank sheet were chosen to be identical to those in Yoon et al. [39], as given in Table 2. Then, several FE simulations of deep drawing were carried out to determine the conditions under which a properly uniform BHF with less earing can be formed. An additional FE simulation was also performed using varying BHF to observe the reduction of earing formation in deep drawing. In this work, the earing height variation, ∆H, was used to assess the earing reduction in the deep-drawn cup. Figure 5 shows a schematic of calculating ∆H, estimated as where h i is the height of the ith node at the edge of the deep-drawn cup, while h avg is the average height of these nodes, and n is the number of nodes at the edge of the deepdrawn cup.
properties of the blank sheet varied. Accordingly, several material properties wit ent strength coefficients, K, and planar anisotropy, Δr, were chosen to perform the ulations of deep drawing by using the uniform and varying BHFs. For more comp three different strength coefficients (132 MPa, 250 Mpa, and 500 MPa) and three planar anisotropies (0.53, 0.77, and 1.04) were generated to investigate the effect rial properties on the earing reduction in deep drawing with varying BHF.
In order to achieve the optimal varying BHF by using the analytical mod factorial design is applied to determine the fitting parameter of α, β, and γ, as ind Equation (2). In this work, three levels of each factor were assumed as α: 1500, 29 β: -650, -250, -450; and γ: 720, 820, 920. The levels of the three factors above were a so that 100 ≤ BHF ≤ 5000. The upper and lower limits were chosen to avoid the oc of wrinkling or tearing during the deep drawing process. Figure 6 shows the al for the determination of the varying BHF based on the analytical model in con with FE simulation. In the first step, the FE simulation was carried out by using a BHF to obtain the earing formation in the deep-drawn cup. Then, the cup hei found to determine the parameter f using Equation (4). Accordingly, by substitu [f] into Equation (2), the BHF for each segment was derived in terms of α, β, and given values of α, β, and γ, the full factorial design was used to generate matrix c tions for calculating the varying BHF. Next, the FE simulations were conducted w ferent varying BHFs. The criterion used to assess the effect of the varying BHF reduction of earing height variation. Finally, the optimal varying BHF, F(θ), was fully determined in terms of the proper combination of α, β, and γ. Furthermore, for the sake of saving costs, FE simulations of deep drawing were carried out to generate the dataset for developing DNN model, which will be described in Section 3, and also study the effects of the material properties on the earing reduction in deep drawing by using the varying BHF. The same FE model was used while the material properties of the blank sheet varied. Accordingly, several material properties with different strength coefficients, K, and planar anisotropy, ∆r, were chosen to perform the FE simulations of deep drawing by using the uniform and varying BHFs. For more comparisons, three different strength coefficients (132 MPa, 250 Mpa, and 500 MPa) and three different planar anisotropies (0.53, 0.77, and 1.04) were generated to investigate the effect of material properties on the earing reduction in deep drawing with varying BHF.
In order to achieve the optimal varying BHF by using the analytical model, a full factorial design is applied to determine the fitting parameter of α, β, and γ, as indicated in Equation (2). In this work, three levels of each factor were assumed as α: 1500, 2900, 4300; β: -650, -250, -450; and γ: 720, 820, 920. The levels of the three factors above were assumed so that 100 ≤ BHF ≤ 5000. The upper and lower limits were chosen to avoid the occurrence of wrinkling or tearing during the deep drawing process. Figure 6 shows the algorithm for the determination of the varying BHF based on the analytical model in conjunction with FE simulation. In the first step, the FE simulation was carried out by using a uniform BHF to obtain the earing formation in the deep-drawn cup. Then, the cup height was found to determine the parameter f using Equation (4). Accordingly, by substituting the [f ] into Equation (2), the BHF for each segment was derived in terms of α, β, and γ. With given values of α, β, and γ, the full factorial design was used to generate matrix combinations for calculating the varying BHF. Next, the FE simulations were conducted with different varying BHFs. The criterion used to assess the effect of the varying BHF was the reduction of earing height variation. Finally, the optimal varying BHF, F(θ), was successfully determined in terms of the proper combination of α, β, and γ.

Deep Neural Network (DNN) Model
In fact, it is very difficult to establish the functional relationship between the earing height variation output and varying BHF input parameters by conventional mathematical methods. Figure 7 shows the DNN architecture, which is an efficiently predictive model for the prediction of earing height variation in this study. The three-layer DNN model including input, hidden, and output layers was performed by using a feedforward backpropagation neural network in the MATLAB environment. The transfer functions in the hidden and output layers were in the form of hyperbolic tangent sigmoid (tansig) and linear (purelin), respectively. The varying BHF of six segments and the earing height variation were defined as the input and output parameters, respectively. Thus, the number of neurons in the input and output layers were six and one, respectively. Figure 7 shows the DNN model consisting of three hidden layers with seven neurons for each layer. These were determined through the iteration method indicated in Figure 8. The performance of the DNN model was evaluated by mean square error (MSE) and coefficient of determination (R 2 ) as follows: Figure 6. Algorithm for determination of optimal varying BHF with full factorial design of α, β, γ based on the analytical model by using the FE simulation.

Deep Neural Network (DNN) Model
In fact, it is very difficult to establish the functional relationship between the earing height variation output and varying BHF input parameters by conventional mathematical methods. Figure 7 shows the DNN architecture, which is an efficiently predictive model for the prediction of earing height variation in this study. The three-layer DNN model including input, hidden, and output layers was performed by using a feedforward backpropagation neural network in the MATLAB environment. The transfer functions in the hidden and output layers were in the form of hyperbolic tangent sigmoid (tansig) and linear (purelin), respectively. The varying BHF of six segments and the earing height variation were defined as the input and output parameters, respectively. Thus, the number of neurons in the input and output layers were six and one, respectively. Figure 7 shows the DNN model consisting of three hidden layers with seven neurons for each layer. These were determined through the iteration method indicated in Figure 8. The performance of the Metals 2021, 11, 395 9 of 22 DNN model was evaluated by mean square error (MSE) and coefficient of determination (R 2 ) as follows: where y,ŷ, and y are the target, predicted, and average target values of earing height variation respectively, and n is the number of data samples used in the DNN model.
where y, ŷ, and ӯ are the target, predicted, and average target values of earing height variation respectively, and n is the number of data samples used in the DNN model.  The full factorial design of varying BHF input parameters was used to generate the number of data samples for developing the DNN model. Three levels with low-middlehigh values of each varying BHF (N) were chosen as 100, 2500, and 4900, respectively. Accordingly, the total of 729 models were simulated by the FE analysis for obtaining the where y, ŷ, and ӯ are the target, predicted, and average target values of earing height variation respectively, and n is the number of data samples used in the DNN model.  The full factorial design of varying BHF input parameters was used to generate the number of data samples for developing the DNN model. Three levels with low-middlehigh values of each varying BHF (N) were chosen as 100, 2500, and 4900, respectively. Accordingly, the total of 729 models were simulated by the FE analysis for obtaining the The full factorial design of varying BHF input parameters was used to generate the number of data samples for developing the DNN model. Three levels with low-middlehigh values of each varying BHF (N) were chosen as 100, 2500, and 4900, respectively. Accordingly, the total of 729 models were simulated by the FE analysis for obtaining the earing height variation output. Prior to training, the input and output data were normalized in the range of 0 and 1 to increase accuracy and speed of the network. Then, the DNN was trained and validated by 442 and 189 data samples of input and output, respectively. The remaining data samples were not used in the DNN model because the earing height variation output could not be obtained from the simulation due to the occurrence of tearing during the deep drawing process.

DNN-Genetic Algorithm (GA) Optimization Model
GA has been widely used to solve the optimization problems based on the mechanics of natural genetics and natural selection. In the present study, the developed DNN model in conjunction with GA was proposed and carried out for the prediction and optimization of varying BHF parameters in the objective of minimum earing height variation. The uniqueness of the integrated DNN-GA model is capable of finding the optimized variables without giving the explicitly mathematical function of the objective, since the trained DNN model is used as an objective function in GA to calculate the fitness value. The procedure for implementing the DNN-GA optimization algorithm is presented in Figure 8. First, the algorithm is started with the development of the DNN model. By changing the number of hidden neurons and layers, a number of network architectures are investigated. The architecture that shows the highest performance in terms of MSE and R 2 for prediction is considered as the best DNN model. Next, the GA is performed by randomly generating the initial population. Based on the best trained DNN model, the fitness is calculated for all initial individuals. As the criteria optimization is not satisfied, the new generation is reproduced by using three operators of selection, crossover, and mutation. The process is repeated until the maximum population size or generation size is reached. In this work, the population and generation sizes were given as 100 and 200, respectively. The probability of crossover rate and mutation were set as 0.85 and 0.01, respectively, to establish the DNN-GA model and then optimize the varying BHFs for minimum earing height variation.

Verification of the Present Finite Element Analysis (FEA)
In order to verify the accuracy of the present finite element analysis (FEA), the analytical and experimental results presented by Yoon et al. [36] for aluminum AA2090-T3 sheet material were used. The FE model was generated with the same geometry and material properties as given in Tables 1 and 2. Additionally, the same blank holding force of 22.2 kN was used in this verification. As shown in Figure 9, the prediction of the earing height profile using the present FEA is represented in comparison with the previous report [36]. Figure 9 shows that the number of predicted ears correlated well with the analytical and experimental findings [36]. The result based on the present FEA showed that the cup height at 0 • was higher than the height at 90 • , which was in good agreement with the previous experiment. This demonstrates that the present FEA predicted the earing trend correctly and showed good consistency with the earlier works [36]. In addition, the predicted earing could also be confirmed by the theoretical analysis using planar anisotropy, ∆r, which is defined as ∆r = r 0 − 2r 45 + r 90 2 where r 0 , r 45 , and r 90 are the anisotropy coefficients of the sheet metal at 0 • , 45 • , and 90 • to the RD, respectively. According to the theoretical analysis, ears are formed at 45 • to the RD as the negative planar anisotropy; in contrast, ears are observed at 0 • and 90 • to the RD if the planar anisotropy is positive [40]. Here, the anisotropy coefficients of AA2090-T3 are given in Table 2 [38]. As the planar anisotropy was negative (∆r = −1.125), ears were found at 45 • to the RD, as indicated in Figure 9. Therefore, the present FEA revealed that the prediction of earing was in good agreement with the previous analytical and experimental results [36] as well as the theoretical analysis. This means that the present FEA was successfully verified and could be used in this study. the prediction of earing was in good agreement with the previous analytical and experimental results [36] as well as the theoretical analysis. This means that the present FEA was successfully verified and could be used in this study. Figure 9. Prediction of earing height profile of AA2093-T3 using the present FEA (Hill's1948) compared to the analytical formula and experimental data in the previous study [36] for verification.

Prediction of Earing by Using the Uniform BHF
First, the FE simulations of deep drawing were conducted using BHFs of 0.05 kN and 15 kN, respectively. It was observed that wrinkling occurred when the BHF was not enough to prevent the buckling of the blank sheet, and tearing appeared due to the insufficient material flow when the BHF was excessive [19]. Next, the predictions of earing in deep drawing using various uniform BHFs of 1 kN, 3 kN, and 5 kN were carried out, as shown in Figure 10. The present FEA predicted that the height of the deep-drawn cup monotonically decreases from 0° to 45° to the RD and monotonically increases from 45° to 90° to the RD. This is due to the contribution of the yield stress directionality [36]. As the planar anisotropy was positive (Δr = 0.530), the formation of earing was observed at 0° and 90° to the RD, which well correlated with the theoretical analysis [40]. It can be seen that four main ears formed at 0°, 90°, 180°, and 270° in the cup deep drawing process. This is partially attributed to the slower flow of sheet metal into the die cavity in these regions than in the other ones. The results show that the uniform BHFs of 1 kN and 3 kN exhibited two bigger ears at 0° and 180° as well as two smaller ears at 90° and 270°. Meanwhile, the uniform BHF of 5 kN predicted almost the same four ears at 0°, 90°, 180°, and 270°. The results revealed that the height of earing, the peaks of the waves shown in Figure 10, increased when the BHF increased [15]. A significant increase of earing height was observed at 90° and 270° with the larger BHF. By using Equation (5), the earing height variation, ΔH, can then be determined. The numerical values of ΔH for uniform BHFs of 1 kN, 3 kN, and 5 kN were 12.38 mm, 12.94 mm, and 13.76 mm, respectively. This result illustrates that the ΔH also increased with the increase in BHF. Figure 9. Prediction of earing height profile of AA2093-T3 using the present FEA (Hill's1948) compared to the analytical formula and experimental data in the previous study [36] for verification.

Prediction of Earing by Using the Uniform BHF
First, the FE simulations of deep drawing were conducted using BHFs of 0.05 kN and 15 kN, respectively. It was observed that wrinkling occurred when the BHF was not enough to prevent the buckling of the blank sheet, and tearing appeared due to the insufficient material flow when the BHF was excessive [19]. Next, the predictions of earing in deep drawing using various uniform BHFs of 1 kN, 3 kN, and 5 kN were carried out, as shown in Figure 10. The present FEA predicted that the height of the deep-drawn cup monotonically decreases from 0 • to 45 • to the RD and monotonically increases from 45 • to 90 • to the RD. This is due to the contribution of the yield stress directionality [36]. As the planar anisotropy was positive (∆r = 0.530), the formation of earing was observed at 0 • and 90 • to the RD, which well correlated with the theoretical analysis [40]. It can be seen that four main ears formed at 0 • , 90 • , 180 • , and 270 • in the cup deep drawing process. This is partially attributed to the slower flow of sheet metal into the die cavity in these regions than in the other ones. The results show that the uniform BHFs of 1 kN and 3 kN exhibited two bigger ears at 0 • and 180 • as well as two smaller ears at 90 • and 270 • . Meanwhile, the uniform BHF of 5 kN predicted almost the same four ears at 0 • , 90 • , 180 • , and 270 • . The results revealed that the height of earing, the peaks of the waves shown in Figure 10, increased when the BHF increased [15]. A significant increase of earing height was observed at 90 • and 270 • with the larger BHF. By using Equation (5), the earing height variation, ∆H, can then be determined. The numerical values of ∆H for uniform BHFs of 1 kN, 3 kN, and 5 kN were 12.38 mm, 12.94 mm, and 13.76 mm, respectively. This result illustrates that the ∆H also increased with the increase in BHF.
Since the r-values were used in Hill's 1948 yield function for FE simulation of deep drawing, the strain state of the deep-drawn cup was observed. Figure 11 shows the contour plots of the equivalent plastic strain of the deep-drawn cup using the uniform BHFs of 1 kN, 3 kN, and 5 kN. Note that the labels of SNEG and SPOS in the legend represent the interior and exterior surfaces of the deep-drawn cup, respectively. Due to the quadratic Hill yield function, the equivalent plastic strain distributions are predominantly concentrated in four regions, leading to four major ears [37]. As shown in Figure 11, the maximum (Max) equivalent plastic strains were obtained to be 0.78, 0.81, and 1.03 by using the uniform BHFs of 1 kN, 3 kN, and 5 kN, respectively. The result revealed that the equivalent plastic strain increased when the BHF increased. Moreover, comparing the various uniform BHFs of deep drawing illustrated that the equivalent plastic strain in the case of BHF using 5 kN was more critical than that of using the uniform BHFs of 1 kN and 3 kN. This high level of plastic strain can induce tearing defects in deep drawing when a larger BHF is introduced [14]. With these drawbacks, it can be seen that there is a limit to reducing earing formation in deep drawing with uniform BHF. Therefore, varying the BHF by using the segmental blank holder is proposed to reduce earing in deep drawing. Since the r-values were used in Hill's 1948 yield function for FE simulation of deep drawing, the strain state of the deep-drawn cup was observed. Figure 11 shows the contour plots of the equivalent plastic strain of the deep-drawn cup using the uniform BHFs of 1 kN, 3 kN, and 5 kN. Note that the labels of SNEG and SPOS in the legend represent the interior and exterior surfaces of the deep-drawn cup, respectively. Due to the quadratic Hill yield function, the equivalent plastic strain distributions are predominantly concentrated in four regions, leading to four major ears [37]. As shown in Figure 11, the maximum (Max) equivalent plastic strains were obtained to be 0.78, 0.81, and 1.03 by using the uniform BHFs of 1 kN, 3 kN, and 5 kN, respectively. The result revealed that the equivalent plastic strain increased when the BHF increased. Moreover, comparing the various uniform BHFs of deep drawing illustrated that the equivalent plastic strain in the case of BHF using 5 kN was more critical than that of using the uniform BHFs of 1 kN and 3 kN. This high level of plastic strain can induce tearing defects in deep drawing when a larger BHF is introduced [14]. With these drawbacks, it can be seen that there is a limit to reducing earing formation in deep drawing with uniform BHF. Therefore, varying the BHF by using the segmental blank holder is proposed to reduce earing in deep drawing.  Since the r-values were used in Hill's 1948 yield function for FE simulation of deep drawing, the strain state of the deep-drawn cup was observed. Figure 11 shows the contour plots of the equivalent plastic strain of the deep-drawn cup using the uniform BHFs of 1 kN, 3 kN, and 5 kN. Note that the labels of SNEG and SPOS in the legend represent the interior and exterior surfaces of the deep-drawn cup, respectively. Due to the quadratic Hill yield function, the equivalent plastic strain distributions are predominantly concentrated in four regions, leading to four major ears [37]. As shown in Figure 11, the maximum (Max) equivalent plastic strains were obtained to be 0.78, 0.81, and 1.03 by using the uniform BHFs of 1 kN, 3 kN, and 5 kN, respectively. The result revealed that the equivalent plastic strain increased when the BHF increased. Moreover, comparing the various uniform BHFs of deep drawing illustrated that the equivalent plastic strain in the case of BHF using 5 kN was more critical than that of using the uniform BHFs of 1 kN and 3 kN. This high level of plastic strain can induce tearing defects in deep drawing when a larger BHF is introduced [14]. With these drawbacks, it can be seen that there is a limit to reducing earing formation in deep drawing with uniform BHF. Therefore, varying the BHF by using the segmental blank holder is proposed to reduce earing in deep drawing.   Figure 12 shows the variation of BHF and its effect on the earing height profiles with three cases of varying BHF obtained by using the analytical model. According to Equation (2), three combinations of α, β, and γ were selected from the full factorial design to calculate three cases of varying BHF, as indicated in Figure 12a. It was observed that the BHF continuously varied region by region in the flange of the sheet metal with a greater amount of BHF at the valley region corresponding to the earing profile to control the plastic-flow phase of material into the die cavity. This confirms that the proposed analytical model was properly established to determine the varying BHF for investigation of earing reduction in deep drawing. Figure 12b shows a comparison of earing height profiles by using the uniform and three-case varying BHFs. For the sake of comparison, the predicted earing using the uniform BHF of 3 kN was used as a reference. It was also noted that all the FEA conditions including material properties, type and size of mesh, and friction coefficient were the same in the case of using uniform BHF, except for applying different varying BHFs by using a segmental blank holder. The result reveals that the varying BHF had a significant effect on the cup height as well as the reduction of earing. As shown in Figure 12b, the cup height increased by using the varying BHF, while the earing height variation noticeably decreased. In comparison with the uniform BHF, it was observed that the cup height considerably increased at 0 • and 45 • to the RD while little change of the cup height was observed at 90 • to the RD with the varying BHF, leading to the decreased earing height variation. The result shows that with the different varying BHFs, the reduction of earing also varied. The numerical values of earing height variation, ∆H, were calculated to be 10.21 mm, 9.57 mm, and 8.82 mm by using the different varying BHFs with case 1, case 2, and case 3, respectively. These values were less than that of 12.94 mm in the case of using the uniform BHF. This result demonstrates that the earing in deep drawing can be reduced with varying BHF. In order to obtain the greatest reduction in earing, the optimal varying BHF was determined based on the analytical model in combination with the FE simulation. According to the procedure as indicated in Figure 6, the optimal combination of α, β, and γ can be found as α = 4300, β = −650, and γ = 720. By substituting these values into Equation (2), the final equation for the calculation of optimal varying BHF is obtained as Figure 12 shows the variation of BHF and its effect on the earing height profiles with three cases of varying BHF obtained by using the analytical model. According to Equation (2), three combinations of α, β, and γ were selected from the full factorial design to calculate three cases of varying BHF, as indicated in Figure 12a. It was observed that the BHF continuously varied region by region in the flange of the sheet metal with a greater amount of BHF at the valley region corresponding to the earing profile to control the plastic-flow phase of material into the die cavity. This confirms that the proposed analytical model was properly established to determine the varying BHF for investigation of earing reduction in deep drawing. Figure 12b shows a comparison of earing height profiles by using the uniform and three-case varying BHFs. For the sake of comparison, the predicted earing using the uniform BHF of 3 kN was used as a reference. It was also noted that all the FEA conditions including material properties, type and size of mesh, and friction coefficient were the same in the case of using uniform BHF, except for applying different varying BHFs by using a segmental blank holder. The result reveals that the varying BHF had a significant effect on the cup height as well as the reduction of earing. As shown in Figure 12b, the cup height increased by using the varying BHF, while the earing height variation noticeably decreased. In comparison with the uniform BHF, it was observed that the cup height considerably increased at 0° and 45° to the RD while little change of the cup height was observed at 90° to the RD with the varying BHF, leading to the decreased earing height variation. The result shows that with the different varying BHFs, the reduction of earing also varied. The numerical values of earing height variation, ΔH, were calculated to be 10.21 mm, 9.57 mm, and 8.82 mm by using the different varying BHFs with case 1, case 2, and case 3, respectively. These values were less than that of 12.94 mm in the case of using the uniform BHF. This result demonstrates that the earing in deep drawing can be reduced with varying BHF. In order to obtain the greatest reduction in earing, the optimal varying BHF was determined based on the analytical model in combination with the FE simulation. According to the procedure as indicated in Figure 6, the optimal combination of α, β, and γ can be found as α = 4300, β = -650, and γ = 720. By substituting these values into Equation (2), the final equation for the calculation of optimal varying BHF is obtained as

Effect of Varying BHF on Earing Reduction
650   Figure 13 shows the DNN model prediction for earing height variation with the coefficient of determinations and compared the results between the target and predicted values. The results show that the DNN model with three hidden layers of seven neurons in each layer accurately reproduced the earing height variation with the best value of overall R 2 = 0.99. The comparison revealed that the DNN prediction was in good agreement with the FE results. This indicates that the developed DNN model was reasonable and able to be used as the fitness function in GA for optimization of the varying BHF parameters. efficient of determinations and compared the results between the target and predicted values. The results show that the DNN model with three hidden layers of seven neurons in each layer accurately reproduced the earing height variation with the best value of overall R 2 = 0.99. The comparison revealed that the DNN prediction was in good agreement with the FE results. This indicates that the developed DNN model was reasonable and able to be used as the fitness function in GA for optimization of the varying BHF parameters.

Optimization of Varying BHF with Integrated DNN-GA Model
The optimized varying BHF for minimization of earing height variation was successfully obtained by using the integrated DNN-GA model, as given in Table 3.
The results before and after optimizations were well compared. Note that the varying BHF indicated by before optimization was selected as the best result among all of the full factorial designs. Figure 14 shows the earing height profile by using uniform BHF in a comparison of varying BHFs between before and after optimizations. It can be seen that there was a negligible difference between the results before and after optimization. This means that the integrated DNN-GA can accurately predict the optimization of varying BHF. In order to obtain the result before optimization, it is time-consuming and has a high manufacturing cost for implementing a great number of experiments by using full factorial design. Thus, the integrated DNN-GA model is a robust and efficient alternative for overcoming these disadvantages. In addition, compared to the uniform BHF, the results

Optimization of Varying BHF with Integrated DNN-GA Model
The optimized varying BHF for minimization of earing height variation was successfully obtained by using the integrated DNN-GA model, as given in Table 3.
The results before and after optimizations were well compared. Note that the varying BHF indicated by before optimization was selected as the best result among all of the full factorial designs. Figure 14 shows the earing height profile by using uniform BHF in a comparison of varying BHFs between before and after optimizations. It can be seen that there was a negligible difference between the results before and after optimization. This means that the integrated DNN-GA can accurately predict the optimization of varying BHF. In order to obtain the result before optimization, it is time-consuming and has a high manufacturing cost for implementing a great number of experiments by using full factorial design. Thus, the integrated DNN-GA model is a robust and efficient alternative for overcoming these disadvantages. In addition, compared to the uniform BHF, the results after optimization by the DNN-GA model were noticeably decreased in earing height variation. The significant increase of cup height at 45 • and nearly retained cup height at 90 • to the RD resulted in the reduction of the earing height variation. In comparison with the ∆H of 12.94 mm from the uniform BHF, the use of the varying BHF predicted by the DNN-GA model showed a considerable decrease of ∆H by~52% to 6.24 mm, as shown in Table 3. Furthermore, Figure 14b compares the relative earing height between the uniform and varying BHFs. The relative earing height is determined as where H cup and H cup min are the cup height and minimum cup height of the deep-drawn cup, respectively. In the comparison between before and after optimization, the results also showed the insignificant difference in the percentage of earing height above the minimum cup height. However, in the case of using the uniform BHF, it revealed the result of 12-13%, which decreased significantly to 4-7% with the optimized varying BHF by the integrated DNN-GA model. The result reveals the remarkable reduction of earing height variation by using the varying BHF. This means that the DNN-GA model performed in this study is capable of the prediction and optimization of input parameters with minimization of the objective. after optimization by the DNN-GA model were noticeably decreased in earing height variation. The significant increase of cup height at 45° and nearly retained cup height at 90° to the RD resulted in the reduction of the earing height variation. In comparison with the ΔH of 12.94 mm from the uniform BHF, the use of the varying BHF predicted by the DNN-GA model showed a considerable decrease of ΔH by ~52% to 6.24 mm, as shown in Table  3. Furthermore, Figure 14b compares the relative earing height between the uniform and varying BHFs. The relative earing height is determined as where H cup and H cup min are the cup height and minimum cup height of the deep-drawn cup, respectively. In the comparison between before and after optimization, the results also showed the insignificant difference in the percentage of earing height above the minimum cup height. However, in the case of using the uniform BHF, it revealed the result of 12-13%, which decreased significantly to 4-7% with the optimized varying BHF by the integrated DNN-GA model. The result reveals the remarkable reduction of earing height variation by using the varying BHF. This means that the DNN-GA model performed in this study is capable of the prediction and optimization of input parameters with minimization of the objective.

Comparison of Analytical and DNN-GA Models
In the present study, a comparison of analytical and DNN-GA models used to obtain the optimal varying BHF was well discussed. With the analytical model, the optimal varying BHF was calculated by using Equation (9) and was also determined from the DNN model using GA. The magnitudes of two optimal varying BHFs are given in Table 3 for the purpose of comparison. This indicates that the two approaches can provide a similar result of varying BHF. The earing height variation obtained by two optimal results also exhibited a negligible difference. Figure 15 shows the results of earing height and relative earing height profiles of the deep-drawn cup by using the optimal varying BHFs in a comparison between the analytical and DNN-GA models. It can be obviously seen that the optimal varying BHFs obtained by the analytical and DNN-GA models resulted in the same trend of earing height profiles with a slight discrepancy. This implies that the developed DNN-GA model in this study is comparable and efficient to achieve the optimization of varying BHF, which leads to the reduction of earing in deep drawing. Figure 15a shows that both optimal varying BHFs had a substantial influence on the earing formation. In fact, compared to the uniform BHF, the optimal varying BHFs revealed a considerable decrease in earing height variation by~45% and~52% with the analytical and DNN-GA models, respectively, as indicated in Table 3.
oped DNN-GA model in this study is comparable and efficient to achieve the optimization of varying BHF, which leads to the reduction of earing in deep drawing. Figure 15a shows that both optimal varying BHFs had a substantial influence on the earing formation. In fact, compared to the uniform BHF, the optimal varying BHFs revealed a considerable decrease in earing height variation by ~45% and ~52% with the analytical and DNN-GA models, respectively, as indicated in Table 3.
In addition, Figure 15b shows the relative earing height in comparison between the analytical and DNN-GA models. It can be observed that the percentage of earing height above the minimum cup height by using two optimal results was not much different apparently. However, it revealed a remarkable decrease compared to the uniform BHF. The result implies that the varying BHF, which can be optimized by both the analytical and integrated DNN-GA models, using the segmental blank holder is more reasonable for the reduction of earing formation than the uniform BHF using the conventional blank holder in the deep drawing process. Figure 15. Comparison of (a) earing height and (b) % earing height profiles of AA1050-O predicted by the optimal varying BHFs using the analytical and DNN-GA models.

Thickness Variation in Different Directions
The current results have clearly shown the significant effect of the varying BHF on earing reduction in deep drawing. The ability and reliability of the proposed method in the thickness variations of deep-drawn cups are discussed. Figure 16 shows the sheet thickness distribution of deep-drawn cups by the use of uniform BHF in comparison to the optimal varying BHFs obtained by the analytical and DNN-GA models. The results Figure 15. Comparison of (a) earing height and (b) % earing height profiles of AA1050-O predicted by the optimal varying BHFs using the analytical and DNN-GA models.
In addition, Figure 15b shows the relative earing height in comparison between the analytical and DNN-GA models. It can be observed that the percentage of earing height above the minimum cup height by using two optimal results was not much different apparently. However, it revealed a remarkable decrease compared to the uniform BHF. The result implies that the varying BHF, which can be optimized by both the analytical and integrated DNN-GA models, using the segmental blank holder is more reasonable for the reduction of earing formation than the uniform BHF using the conventional blank holder in the deep drawing process.

Thickness Variation in Different Directions
The current results have clearly shown the significant effect of the varying BHF on earing reduction in deep drawing. The ability and reliability of the proposed method in the thickness variations of deep-drawn cups are discussed. Figure 16 shows the sheet thickness distribution of deep-drawn cups by the use of uniform BHF in comparison to the optimal varying BHFs obtained by the analytical and DNN-GA models. The results reveals a slight decrease in thickness by using both optimal varying BHFs compared to the uniform BHF. It can be observed that the maximum thinning was found in the corner region of the deep-drawn cup. As shown in Figure 16, the numerical values of maximum thinning using the optimal varying BHFs obtained by the analytical and DNN-GA models were 1.41 and 1.44, respectively, which were only a little lower than that of 1.48 with the result of using uniform BHF. The result showed an increase of thinning in the deep-drawn cup by using the optimal varying BHF, but was very negligible.
In this study, the maximum thinning rate, θ max was used to assess the occurrence of tearing in deep drawing, expressed as [41] θ max = 1 − t min t initial (11) where t initial is the initial thickness of the blank sheet (t initial = 1.6 mm), and t min is the minimum thickness of the deep-drawn cup (t min = 1.48 mm, 1.41 mm, and 1.44 mm for using the uniform and two optimal varying BHFs, respectively). By substituting these values into Equation (11), the maximum thinning rates were obtained as 0.08, 0.12, and 0.10, respectively. Generally, the maximum thinning rate cannot exceed the critical thinning rate of the sheet material, which was taken as 0.25, according to prior experience [12]. The result revealed that the maximum thinning rates by using the uniform and optimal varying BHFs were much lower than the critical thinning rate. This implies that there would be no occurrence of tearing in the deep-drawn cup by using the varying BHF.
Metals 2021, 11, 395 17 of 23 reveals a slight decrease in thickness by using both optimal varying BHFs compared to the uniform BHF. It can be observed that the maximum thinning was found in the corner region of the deep-drawn cup. As shown in Figure 16, the numerical values of maximum thinning using the optimal varying BHFs obtained by the analytical and DNN-GA models were 1.41 and 1.44, respectively, which were only a little lower than that of 1.48 with the result of using uniform BHF. The result showed an increase of thinning in the deep-drawn cup by using the optimal varying BHF, but was very negligible. In this study, the maximum thinning rate, θmax was used to assess the occurrence of tearing in deep drawing, expressed as [41] min max 1 initial t t θ = − (11) where tinitial is the initial thickness of the blank sheet (tinitial = 1.6 mm), and tmin is the minimum thickness of the deep-drawn cup (tmin = 1.48 mm, 1.41 mm, and 1.44 mm for using the uniform and two optimal varying BHFs, respectively). By substituting these values into Equation (11), the maximum thinning rates were obtained as 0.08, 0.12, and 0.10, respectively. Generally, the maximum thinning rate cannot exceed the critical thinning rate of the sheet material, which was taken as 0.25, according to prior experience [12]. The result revealed that the maximum thinning rates by using the uniform and optimal varying BHFs were much lower than the critical thinning rate. This implies that there would be no occurrence of tearing in the deep-drawn cup by using the varying BHF.
In order to further observe the thickness variation, the thickness of the deep-drawn cup was estimated along the RD, DD, and TD at three different regions of the cup: (I) cup bottom, (II) cup corner, and (III) cup wall, as shown in Figure 17. Note that the reference point is taken from the center of the cup bottom. Figure 18 shows the thickness variations of the deep-drawn cup along the RD, DD, and TD using the uniform and two optimal varying BHFs. The results show a non-uniform thickness distribution at different regions of the deep-drawn cup due to the different strain states [42]. At the top of the cup wall region, the thickness distribution in the RD, DD, and TD by using the uniform and optimal varying BHFs was higher than the initial thickness (tinitial = 1.6 mm). This is due to the radial tensile and circumferential compressive stress states in the flange region. According to the previous reports [42,43], the compressive stress along the circumferential direction can contribute to the increase in the thickness of the blank sheet. A decreased thickness was observed at the cup bottom due to the experiences of the biaxial tensile stress state and the friction between the punch and the blank sheet. Meanwhile, the maximum thinning In order to further observe the thickness variation, the thickness of the deep-drawn cup was estimated along the RD, DD, and TD at three different regions of the cup: (I) cup bottom, (II) cup corner, and (III) cup wall, as shown in Figure 17. Note that the reference point is taken from the center of the cup bottom. Figure 18 shows the thickness variations of the deep-drawn cup along the RD, DD, and TD using the uniform and two optimal varying BHFs. The results show a non-uniform thickness distribution at different regions of the deep-drawn cup due to the different strain states [42]. At the top of the cup wall region, the thickness distribution in the RD, DD, and TD by using the uniform and optimal varying BHFs was higher than the initial thickness (t initial = 1.6 mm). This is due to the radial tensile and circumferential compressive stress states in the flange region. According to the previous reports [42,43], the compressive stress along the circumferential direction can contribute to the increase in the thickness of the blank sheet. A decreased thickness was observed at the cup bottom due to the experiences of the biaxial tensile stress state and the friction between the punch and the blank sheet. Meanwhile, the maximum thinning of the deep-drawn cup was found at the cup corner region due to the plane strain condition; however, it was very negligible. In the comparison between the results of using the analytical and DNN-GA models, the thickness in different directions revealed a similar variation. Additionally, it showed a slight decrease in thickness compared to the uniform BHF due to the increase in the deep-drawn cup height [44]. The interesting point is the increase in thickness in the case of using the optimal varying BHF observed at the cup wall region in the TD. This is due to the decrease in varying BHF toward the end of the deep drawing process, leading to the lower tensile stress in the wall region of the blank sheet. Overall, the difference in thickness of the deep-drawn cup between the uniform and varying BHFs was insignificant. It can be concluded that the varying BHF is an effective method to minimize earing formation in deep drawing with little influence on the thickness variation in the deep-drawn cup. Although the blank holding force of a segmental holder was held constant in the present study, the time dependence of the BHF seems to be important for further study.
deep drawing process, leading to the lower tensile stress in the wall region of the blank sheet. Overall, the difference in thickness of the deep-drawn cup between the uniform and varying BHFs was insignificant. It can be concluded that the varying BHF is an effective method to minimize earing formation in deep drawing with little influence on the thickness variation in the deep-drawn cup. Although the blank holding force of a segmental holder was held constant in the present study, the time dependence of the BHF seems to be important for further study.

Effect of the Material Properties on Earing Reduction
Even though the varying BHF has been utilized, there is also a limit on the further reduction of the earing since the formation of earing predominantly depends on the material properties. In fact, the material properties, especially the plastic anisotropy, significantly affect the quality of the deep-drawn cup, which results in the formation of earing. For example, the prediction of more than four ears was observed by involving more anisotropy coefficients of yield function in the FE model [45]. It has been also reported that the anisotropy parameters are correlated with the enlarging, number, and direction of ears, while the planar anisotropy is associated with the enlarging and size of ears [40]. Moreover, the contribution of r-value and yield stress directionalities plays an important role in predicting the earing of the deep-drawn cup [36].
Furthermore, FE simulations of deep drawing were also carried out to study the effect of the material properties on the earing reduction in deep drawing by using the varying BHF. Figure 19 shows the prediction of earing height profiles by using the uniform and varying BHFs for AA2090-T3, which is a highly anisotropic material. The results showed that significant earing was predicted due to the high anisotropy. Since the planar anisotropy was negative (∆r = −1.125), the earing formation was observed at 45 • to the RD, which is in good agreement with the theoretical analysis [40]. By using the uniform BHF, it can be seen that the AA2090-T3 sheet material showed a higher discrepancy between the maximum and minimum cup height than that of the AA1050-O sheet material, which is a low-anisotropy material. It has been also observed that the highly anisotropic material of AA2090-T3 exhibited a higher cup height variation than that of the mildly anisotropic material of AA6022-T4 [46]. By introducing the varying BHF, which is determined by using the analytical model, the earing height profile was predicted. As shown in Figure 19a, the result of using the varying BHF revealed an insignificant influence on the reduction of earing formation for the highly anisotropic material of AA2090-T3. By varying the BHF, it was observed that the cup height had an increase at 0 • to the RD; however, very little change on the cup height was observed at 45 • and 90 • to the RD. This may lead to a negligible effect on the reduction of earing height variation by using the varying BHF. Indeed, the earing height variations, ∆H, were determined to be 26.45 mm and 24.04 mm for the uniform and varying BHFs, respectively. This result reveals that the reduction of earing height variation was less than 10%. This can be better observed by the result of relative earing height, as shown in Figure 19b. By using the uniform BHF, the percentage of earing height above the minimum cup height was obtained as~27%, while a negligible decrease to~26% was observed by using the varying BHF. The result showed that the reduction of earing in deep drawing by varying the BHF for the highly anisotropic material of AA2090-T3 was inconsiderable.

Effect of the Material Properties on Earing Reduction
Even though the varying BHF has been utilized, there is also a limit o reduction of the earing since the formation of earing predominantly depend terial properties. In fact, the material properties, especially the plastic anisot cantly affect the quality of the deep-drawn cup, which results in the formati For example, the prediction of more than four ears was observed by involvi sotropy coefficients of yield function in the FE model [45]. It has been also r the anisotropy parameters are correlated with the enlarging, number, and ears, while the planar anisotropy is associated with the enlarging and size Moreover, the contribution of r-value and yield stress directionalities plays role in predicting the earing of the deep-drawn cup [36].
Furthermore, FE simulations of deep drawing were also carried out to fect of the material properties on the earing reduction in deep drawing by u ying BHF. Figure 19 shows the prediction of earing height profiles by using and varying BHFs for AA2090-T3, which is a highly anisotropic material showed that significant earing was predicted due to the high anisotropy. Sin anisotropy was negative (Δr = -1.125), the earing formation was observed at 4 which is in good agreement with the theoretical analysis [40]. By using the u it can be seen that the AA2090-T3 sheet material showed a higher discrepa the maximum and minimum cup height than that of the AA1050-O sheet ma is a low-anisotropy material. It has been also observed that the highly anisotr of AA2090-T3 exhibited a higher cup height variation than that of the mildl material of AA6022-T4 [46]. By introducing the varying BHF, which is deter ing the analytical model, the earing height profile was predicted. As shown i the result of using the varying BHF revealed an insignificant influence on t of earing formation for the highly anisotropic material of AA2090-T3. By vary Moreover, several different material properties with strength coefficient, K, and planar anisotropy, ∆r, were chosen to further investigate the effect of the material properties on the earing reduction in deep drawing by using the varying BHF. Table 4 shows the reduction of earing height variation corresponding to the K and ∆r. The reduction of earing height variation is estimated as where ∆H and ∆H are the earing height variations calculated using the varying and uniform BHFs, respectively. It can be observed that a material with a greater combination of K and ∆r showed less influence on the reduction of earing height variation. As shown in Figure 20, the reduction of earing height variation monotonically increased with the decrease of the K and ∆r simultaneously. This result demonstrates that the lower the material strength, the greater the effect on the reduction of earing formation in deep drawing by using the varying BHF that can be obtained. 21, 11, 395 20 of 23 earing in deep drawing by varying the BHF for the highly anisotropic material of AA2090-T3 was inconsiderable. Figure 19. Results of (a) earing height and (b) % earing height profiles in a comparison between uniform and varying BHFs for a highly anisotropic sheet material of AA2090-T3.
Moreover, several different material properties with strength coefficient, K, and planar anisotropy, Δr, were chosen to further investigate the effect of the material properties on the earing reduction in deep drawing by using the varying BHF. Table 4 shows the reduction of earing height variation corresponding to the K and Δr. The reduction of earing height variation is estimated as ' 100 H reduction H Δ = × Δ (12) where ΔH' and ΔH are the earing height variations calculated using the varying and uniform BHFs, respectively. It can be observed that a material with a greater combination of K and Δr showed less influence on the reduction of earing height variation. As shown in Figure 20, the reduction of earing height variation monotonically increased with the decrease of the K and Δr simultaneously. This result demonstrates that the lower the material strength, the greater the effect on the reduction of earing formation in deep drawing by using the varying BHF that can be obtained.  Figure 19. Results of (a) earing height and (b) % earing height profiles in a comparison between uniform and varying BHFs for a highly anisotropic sheet material of AA2090-T3.

Conclusions
In the present study, the varying blank holding force (BHF) by using the segmental blank holder was proposed to investigate its effect on the earing reduction in the circular deep drawing process of aluminum alloy sheet. The following conclusions of this study can be summarized as follows: (i) With an interval holder-segment angle of 15°, the different varying BHFs obtained by the analytical model resulted in the different influences on earing reduction. The noticeable decrease of earing height variation by ~45% was revealed by the optimal varying BHF, indicating a more reasonable technique with varying BHF than the uniform BHF for earing reduction in deep drawing. (ii) The developed deep neural network (DNN) model with the integration of genetic algorithm (GA) revealed the optimization of varying BHF with a considerable reduction in earing height by ~52%, which showed a superior result to the analytical model. This confirmed that the proposed DNN-GA model effectively predicted and optimized the process parameters in this study.

Conclusions
In the present study, the varying blank holding force (BHF) by using the segmental blank holder was proposed to investigate its effect on the earing reduction in the circular deep drawing process of aluminum alloy sheet. The following conclusions of this study can be summarized as follows: (i) With an interval holder-segment angle of 15 • , the different varying BHFs obtained by the analytical model resulted in the different influences on earing reduction. The noticeable decrease of earing height variation by~45% was revealed by the optimal varying BHF, indicating a more reasonable technique with varying BHF than the uniform BHF for earing reduction in deep drawing.
(ii) The developed deep neural network (DNN) model with the integration of genetic algorithm (GA) revealed the optimization of varying BHF with a considerable reduction in earing height by~52%, which showed a superior result to the analytical model. This confirmed that the proposed DNN-GA model effectively predicted and optimized the process parameters in this study. (iii) For a volume consistency, the non-uniform thickness distribution was observed with an increase of thickness in the transverse direction at the cup wall region by varying BHF, which indicated an improvement in deep-drawability. Although BHF of a segmental holder was held constant in the present study, the time dependence of the BHF seems to be important for further study. (iv) Furthermore, the present study revealed that, due to the effect of material properties on earing formation, the lower the material strength, the greater the effect on the earing reduction in deep drawing with the varying BHF will be. Funding: This research received no external funding.