Tool Wear Prediction in the Forming of Automotive DP980 Steel Sheet Using Statistical Sensitivity Analysis and Accelerated U-Bending Based Wear Test

The forming process of ultra-high-strength steel (UHSS) may cause premature damage to the tool surface due to the high forming pressure. The damage to and wear of the tool surface increase maintenance costs and deteriorate the surface quality of the sheet products. Hence, a reliable prediction model for tool wear can help in the efficient management of the quality and productivity of formed sheet products of UHSS. In this study, a methodology is proposed for predicting the wear behavior of stamping tools that are used in the forming process of DP980 steel sheet. Pin-on-disk tests were conducted based on the Taguchi method to develop the tool wear prediction model. Using statistical analysis based on the signal-to-noise (S/N) ratio and ANOVA result, the contact pressure and the sliding distance were selected as the major contact parameters for tool wear. The Archard wear model has a limitation in predicting the nonlinear behavior of tool wear. Therefore, the power-law nonlinear regression model as a function of the contact pressure and the sliding distance was constructed. To verify the reliability of the constructed tool wear prediction model, the U-draw bending tests were designed. The modified U-draw bending test, which accelerates tool wear, is newly designed to evaluate the tool wear for different contact pressures and sliding distances. The modified die generated a contact pressure four times higher than that of the conventional die from the finite element (FE) simulation results. The tool wear prediction model was validated by comparing the predicted results with the experimental results of DP980 sheets formed using the physical vapor deposition (PVD) CrN-coated STD11 tool steel.


Introduction
The demand for weight reduction in an auto-body has been continuously increasing in the automotive industry over recent years. For this purpose, non-ferrous metal sheets, such as carbon-fiber-reinforced plastics (CFRPs), and magnesium and aluminum alloys have been regarded as alternatives to conventional steel sheets owing to their low density and high specific stiffness compared to steel. Nevertheless, steel sheets are still the major material candidates for car body design due to their low cost and eco-friendly manufacturing process, in addition to efficient control over their strength and ductility using properly designed heat treatments. These materials are effective in improving fuel efficiency and crashworthiness. For example, Li et al. [1] proposed a typical method for lightweighting by applying thinner and higher strength steel sheets to automotive parts based on the simulation of the body part impact while satisfying crashworthiness. Steel makers have been developing a variety of ultra-high-strength steel (UHSS) sheets that are To improve productivity and reduce maintenance costs, process management based on the accurate prediction of tool wear is required. However, the sheet metal forming process involves various parameters such as lubrication, binder force, tool geometry and contact characteristics. All of which affect tool wear. Therefore, establishing reliable tool management that relies solely on wear experiments based on trial-and-error and the engineer's experience is very challenging and increases both time and cost. To overcome these challenges, numerous research studies have been performed to understand the contact and wear mechanisms by using numerical modeling and simulations. Podra and Andersson [4] analyzed a conical spinning contact considering surface topography both experimentally and with simulation and proposed a modeling and simulation procedure. Kim et al. [5] verified the simulation accuracy of the wear depth in the block-on-ring experiment considering the continuous wear propagation. Lengiewicz and Stupkiewicz [6] proposed a computationally efficient model for wear evolution based on the Archard wear model and verified with the pin-on-disk test problem modeled by three-dimensional finite elements.
Some researchers have tried to develop methodologies for predicting the tool wear behavior in the actual manufacturing processes. For the metal cutting operation, Yen et al. [7] proposed a methodology to predict the tool wear evolution over long cutting periods in the cutting process based on the finite element (FE) analysis. Binder et al. [8] applied the FE-based tool wear simulation method combined with the modified Usui tool wear model in order to predict the wear behavior of the coated cutting tool with a complex shape. Binder et al. [9] extended the constructed simulation method to the tool and process design for accurate tool wear prediction in metal cutting. Pimenov et al. [10] developed a simulation model of stress conditions for the flank wear on the tool's flank surface using the stress intensity equation. They used FE modeling to simulate the evolutionary behav- Figure 1. Example of tool wear occurrence in the sheet metal forming of the B-pillar part using DP980 steel sheets [3].
To improve productivity and reduce maintenance costs, process management based on the accurate prediction of tool wear is required. However, the sheet metal forming process involves various parameters such as lubrication, binder force, tool geometry and contact characteristics. All of which affect tool wear. Therefore, establishing reliable tool management that relies solely on wear experiments based on trial-and-error and the engineer's experience is very challenging and increases both time and cost. To overcome these challenges, numerous research studies have been performed to understand the contact and wear mechanisms by using numerical modeling and simulations. Podra and Andersson [4] analyzed a conical spinning contact considering surface topography both experimentally and with simulation and proposed a modeling and simulation procedure. Kim et al. [5] verified the simulation accuracy of the wear depth in the block-on-ring experiment considering the continuous wear propagation. Lengiewicz and Stupkiewicz [6] proposed a computationally efficient model for wear evolution based on the Archard wear model and verified with the pin-on-disk test problem modeled by three-dimensional finite elements.
Some researchers have tried to develop methodologies for predicting the tool wear behavior in the actual manufacturing processes. For the metal cutting operation, Yen et al. [7] proposed a methodology to predict the tool wear evolution over long cutting periods in the cutting process based on the finite element (FE) analysis. Binder et al. [8] applied the FE-based tool wear simulation method combined with the modified Usui tool wear model in order to predict the wear behavior of the coated cutting tool with a complex shape. Binder et al. [9] extended the constructed simulation method to the tool and process design for accurate tool wear prediction in metal cutting. Pimenov et al. [10] developed a simulation model of stress conditions for the flank wear on the tool's flank surface using the stress intensity equation. They used FE modeling to simulate the evolutionary behavior of tool wear over a long cutting period and adapted the wear rate equation for estimating tool wear for a three-dimensional simulation that considered changes in local thermal, frictional, and wear properties. For the tool wear evaluation of the forging process, the Archard wear model [11] is usually used to the FE-based simulation for predicting the wear behavior of the forging tool. Lee and Jou [12] applied the Archard wear model considering the temperature-dependent hardness of warm forging tools for predicting the wear amount of the forging die. Behrens and Schaefer [13] also used a similar approach for the wear prediction of hot forging tools. They concluded that the Archard wear model shows good performance in predicting the wear phenomenon of warm and hot forging dies. The Archard wear model has been also used to predict the wear behavior of the stamping die in the sheet metal forming process. Hoffmann et al. [14] established an advanced simulation method using the Archard wear model and the geometry update scheme for the sheet metal forming process of the B-pillar reinforcement component. Hoffmann and Nurnberg [15] suggested a new approach to investigate the quantitative tool wear depth in the cylindrical cup drawing process and evaluated the tool wear for various grades of automotive steel sheets. As mentioned above, most wear simulations have adopted the Archard model, considering its efficiency and simplicity. Because the Archard model has reported that the wear amount is linearly proportional to the contact pressure and the sliding distance, there is a limitation in the prediction accuracy of nonlinear evolutionary behaviors of wear. Some researchers have presented the power-law-type prediction models expressed as a function of normal load and sliding distance to improve the prediction accuracy of the wear behavior. Rhee [16] proposed a power law empirical wear equation for predicting the weight loss of the asbestos-reinforced polymers sliding against metal surfaces. Bayer [17] also used the power law equation for developing a general model of surface wear mechanisms for both abrasive and adhesive wear. Hsu et al. [18] summarized various wear models for metals, including the power law equation and complicated equations based on the wear mechanism map. These research results support the idea that a power-law-based wear model can provide improved results for the nonlinear evolutionary wear behavior in the metal forming process. Although powerlaw-based wear models have improved prediction accuracy compared to the conventional Archard-based models, there are insufficient studies on modeling the wear of tool dies used in the stamping of automotive steel sheets, particularly UHSS sheets. Eriksen [19] reported that tool wear is generally due to the abrasive wear mechanism that causes geometrical changes over the die edge for the deep drawing process.
The aim of this study is to propose a new methodology to predict the abrasive wear behavior of stamping tools for dual-phase (DP) steel with the strength of 980 MPa. The power law wear model is implemented in this study to improve the predictive capability for the estimation of abrasive tool wear. First, a pin-on-disk test was performed to analyze the wear behavior and obtain an experimental database for developing the tool wear prediction model. Then, the Taguchi method was employed for an efficient and systematic experimental design. The sensitivity analysis was performed using the signal-to-noise (S/N) ratio and the analysis of variance (ANOVA) result to select the main control parameters affecting the wear behavior. Based on these results, the parameters of the power-law-based tool wear model were identified with selected major control parameters. Finally, the reliability of the tool wear prediction model was verified by comparing the predicted results of the conventional and modified dies for the drawing test with experimental results.

Materials and Tool Wear Characteristics
The contact characteristics acting on the tool surface in the sheet metal forming process were identified through numerical simulations. A B-pillar (or central pillar) of an autobody was used to investigate the range of contact pressure and position where the tool wear is mostly expected. The material used in this study is DP980 steel sheet with the thickness of 1.4 mm, which is produced by POSCO Gwangyang steelworks in South Korea. The mechanical properties are listed in Table 1, which are measured by MTS 810 (MTS Systems Corporation, Eden Prairie, MN, USA). The FE simulation was performed using the dynamic explicit software PAM-STAMP v2017.0 (ESI-Group, Rungis, France). To simplify the simulation procedure, Hill's48 yield function with isotropic hardening law was assumed. As the amount of the tool wear increases, the friction coefficient at the contact interfaces between the sheet and tool surfaces does not change noticeably [20,21]. Therefore, the Coulomb friction coefficient was selected as 0.15, which showed good correlation between the experiment and the numerical simulation [22]. The tool was discretized by rigid shell elements because the elastic deformation of the tool has no significant effect on the contact pressure, as reported by Erosoy-Nurnberg et al. [23]. The blank was modeled with deformable shell elements. The element size at the contact interface between the blank and the tool has to be sufficiently fine to predict the accurate contact pressure. Mesh refinement was used because the smaller the mesh size, the higher the computation time. The initial mesh size was 10 mm, and the mesh refinement method was set to 5 levels so that the minimum mesh size was 0.625 mm. Stamping simulation was performed by Inter Core i7-6700 (8.0 GB of RAM). The computation took about 6 h. Figure 2 shows the contact pressure distribution after the forming process of the B-pillar. The simulation results indicate that the maximum contact pressure of 1.85 GPa was predicted at the sharp curvature of the tool. The potential wear locations due to high contact pressure could be confirmed from the observation of a real stamping tool for the same B-pillar. The preliminary simulation and experimental analysis indicated that the contact pressure was the most significant factor affecting the tool wear. Based on the simulation results, the experimental conditions of the pin-on-disk tests were established for developing a tool wear prediction model of DP980. thickness of 1.4 mm, which is produced by POSCO Gwangyang steelworks in South Korea. The mechanical properties are listed in Table 1, which are measured by MTS 810 (MTS Systems Corporation, Eden Prairie, MN, USA). The FE simulation was performed using the dynamic explicit software PAM-STAMP v2017.0 (ESI-Group, Rungis, France). To simplify the simulation procedure, Hill's48 yield function with isotropic hardening law was assumed. As the amount of the tool wear increases, the friction coefficient at the contact interfaces between the sheet and tool surfaces does not change noticeably [20,21]. Therefore, the Coulomb friction coefficient was selected as 0.15, which showed good correlation between the experiment and the numerical simulation [22]. The tool was discretized by rigid shell elements because the elastic deformation of the tool has no significant effect on the contact pressure, as reported by Erosoy-Nurnberg et al. [23]. The blank was modeled with deformable shell elements. The element size at the contact interface between the blank and the tool has to be sufficiently fine to predict the accurate contact pressure. Mesh refinement was used because the smaller the mesh size, the higher the computation time.
The initial mesh size was 10 mm, and the mesh refinement method was set to 5 levels so that the minimum mesh size was 0.625 mm. Stamping simulation was performed by Inter Core i7-6700 (8.0 GB of RAM). The computation took about 6 h. Figure 2 shows the contact pressure distribution after the forming process of the B-pillar. The simulation results indicate that the maximum contact pressure of 1.85 GPa was predicted at the sharp curvature of the tool. The potential wear locations due to high contact pressure could be confirmed from the observation of a real stamping tool for the same B-pillar. The preliminary simulation and experimental analysis indicated that the contact pressure was the most significant factor affecting the tool wear. Based on the simulation results, the experimental conditions of the pin-on-disk tests were established for developing a tool wear prediction model of DP980.

Experimental Setup for a Pin-on-Disk Test
A pin-on-disk test is the most widely used experimental method to analyze wear characteristics quantitatively according to the process parameters and environmental conditions. In this study, this test method was used to estimate the quantitative amount of wear to develop a wear prediction model for DP980 steel sheets. Figure 3 shows the pin-

Experimental Setup for a Pin-on-Disk Test
A pin-on-disk test is the most widely used experimental method to analyze wear characteristics quantitatively according to the process parameters and environmental conditions. In this study, this test method was used to estimate the quantitative amount of wear to develop a wear prediction model for DP980 steel sheets. Figure 3 shows the pin-on-disk tester comprising of the pin, the disk, and a weight. The TRB 3 Tribometer (CSM Instruments SA, Peseux, Switzerland) was used as the testing machine. The pin comes in contact with the metallic disc with a designated normal force applied by a weight, and wear occurs at the contact interface between the static pin and the rotating disk during the test. The pin was prepared using the same processing conditions as that for the B-pillar dies. The material was STD (Steel Tool Dies)-11 tool steel with proper heat treatment, which is produced by POSCO Pohang steelworks in South Korea. As shown in Figure 4, quenching heat treatment was used to harden the STD11 tool steel by introducing the martensite microstructure. The heat treatment includes multiple heating processes at temperatures of 600, 850, and 1030 • C, for 90 min at each temperature, which was followed by a two-step tempering treatment for 2 h to reduce the brittleness arising from the quenching process standardized in KS D 3753 [24]. The Vickers hardness HV 0.1 of the STD11 tool steel (pin) was 787.9 ± 13.5, measured at a normal force of 0.1 kgf. The Vickers hardness of DP980 steel is about 320 Hv, as reported by Farabi et al. [25]. The pin surfaces were coated with CrN using the physical vapor deposition (PVD) method. The wear test was performed according to various contact parameters using this experimental setup. The chemical composition of the DP980 steel sheet and STD11 tool steel is presented in Table 2. on-disk tester comprising of the pin, the disk, and a weight. The TRB 3 Tribometer (CSM Instruments SA, Peseux, Switzerland) was used as the testing machine. The pin comes in contact with the metallic disc with a designated normal force applied by a weight, and wear occurs at the contact interface between the static pin and the rotating disk during the test. The pin was prepared using the same processing conditions as that for the B-pillar dies. The material was STD (Steel Tool Dies)-11 tool steel with proper heat treatment, which is produced by POSCO Pohang steelworks in South Korea. As shown in Figure 4, quenching heat treatment was used to harden the STD11 tool steel by introducing the martensite microstructure. The heat treatment includes multiple heating processes at temperatures of 600, 850, and 1030 °C, for 90 min at each temperature, which was followed by a two-step tempering treatment for 2 h to reduce the brittleness arising from the quenching process standardized in KS D 3753 [24]. The Vickers hardness HV0.1 of the STD11 tool steel (pin) was 787.9 ± 13.5, measured at a normal force of 0.1 kgf. The Vickers hardness of DP980 steel is about 320 Hv, as reported by Farabi et al. [25]. The pin surfaces were coated with CrN using the physical vapor deposition (PVD) method. The wear test was performed according to various contact parameters using this experimental setup. The chemical composition of the DP980 steel sheet and STD11 tool steel is presented in Table 2.
(a) (b)    on-disk tester comprising of the pin, the disk, and a weight. The TRB 3 Tribometer (CSM Instruments SA, Peseux, Switzerland) was used as the testing machine. The pin comes in contact with the metallic disc with a designated normal force applied by a weight, and wear occurs at the contact interface between the static pin and the rotating disk during the test. The pin was prepared using the same processing conditions as that for the B-pillar dies. The material was STD (Steel Tool Dies)-11 tool steel with proper heat treatment, which is produced by POSCO Pohang steelworks in South Korea. As shown in Figure 4, quenching heat treatment was used to harden the STD11 tool steel by introducing the martensite microstructure. The heat treatment includes multiple heating processes at temperatures of 600, 850, and 1030 °C, for 90 min at each temperature, which was followed by a two-step tempering treatment for 2 h to reduce the brittleness arising from the quenching process standardized in KS D 3753 [24]. The Vickers hardness HV0.1 of the STD11 tool steel (pin) was 787.9 ± 13.5, measured at a normal force of 0.1 kgf. The Vickers hardness of DP980 steel is about 320 Hv, as reported by Farabi et al. [25]. The pin surfaces were coated with CrN using the physical vapor deposition (PVD) method. The wear test was performed according to various contact parameters using this experimental setup. The chemical composition of the DP980 steel sheet and STD11 tool steel is presented in Table 2.

Experimental Design Parameters
An efficient experimental design method is required to construct a reliable wear prediction equation. In this study, the Taguchi method presented by Taguchi and Konishi [26] was used for the experimental design of the pin-on-disk test with various parameters. This method is a powerful and widely used design scheme in engineering analysis because a small number of experimental cases can be used to evaluate the response of control parameters and experimental results in the design space.
To formulate a tool wear prediction model, an orthogonal array should be constructed by investigating the actual sheet metal forming process. Pereira et al. [27] reported that the contact pressure and the sliding distance are essential for evaluating the wear behavior and lifetime of the tool. Rhee [16] confirmed that the sliding velocity can be an important parameter in constructing an empirical wear equation. Based on references, the contact pressure, the sliding distance, and the velocity were selected as the major factors affecting the wear behavior. For three contact parameters, an orthogonal array L 25 (5 3 ) can be used for the experimental design of the pin-on-disk test. The range of the contact parameters should be confirmed for constructing the orthogonal array.
In the previous section, the maximum contact pressure of the B-pillar predicted by the FE simulation was 1.85 GPa, as shown in Figure 2. For the pin-on-disk test, the contact pressure should be transformed into a normal force. Therefore, the normal force corresponding to this contact pressure was calculated and set as the force range of the pin-on-disk test. Because the tool structure is in the range of the elastic deformation during the sheet metal forming process [28,29], the normal force acting on the contact surface was derived using the Hertz contact theory [30]. The Hertz contact theory states that when two elastic bodies are loaded by contact, elastic bodies create a contact area in which the contact pressure is distributed. As shown in Figure 5, the elastic sphere of radius R is indented to a depth d, which induces a contact surface of radius a.

Experimental Design Parameters
An efficient experimental design method is required to construct a reliable wear prediction equation. In this study, the Taguchi method presented by Taguchi and Konishi [26] was used for the experimental design of the pin-on-disk test with various parameters. This method is a powerful and widely used design scheme in engineering analysis because a small number of experimental cases can be used to evaluate the response of control parameters and experimental results in the design space.
To formulate a tool wear prediction model, an orthogonal array should be constructed by investigating the actual sheet metal forming process. Pereira et al. [27] reported that the contact pressure and the sliding distance are essential for evaluating the wear behavior and lifetime of the tool. Rhee [16] confirmed that the sliding velocity can be an important parameter in constructing an empirical wear equation. Based on references, the contact pressure, the sliding distance, and the velocity were selected as the major factors affecting the wear behavior. For three contact parameters, an orthogonal array L25 (5 3 ) can be used for the experimental design of the pin-on-disk test. The range of the contact parameters should be confirmed for constructing the orthogonal array.
In the previous section, the maximum contact pressure of the B-pillar predicted by the FE simulation was 1.85 GPa, as shown in Figure 2. For the pin-on-disk test, the contact pressure should be transformed into a normal force. Therefore, the normal force corresponding to this contact pressure was calculated and set as the force range of the pin-ondisk test. Because the tool structure is in the range of the elastic deformation during the sheet metal forming process [28,29], the normal force acting on the contact surface was derived using the Hertz contact theory [30]. The Hertz contact theory states that when two elastic bodies are loaded by contact, elastic bodies create a contact area in which the contact pressure is distributed. As shown in Figure 5, the elastic sphere of radius R is indented to a depth d, which induces a contact surface of radius a.

= √
(1) The depth d is generated by the applied normal force.
where Er is the reduced elastic modulus, which can be calculated as follows: where ν1 and ν2 are Poisson's ratios and E1 and E2 are the elastic moduli for the elastic bodies denoted by subscripts 1 and 2, respectively, in Figure 5. The distribution of the The depth d is generated by the applied normal force.
where E r is the reduced elastic modulus, which can be calculated as follows: where ν 1 and ν 2 are Poisson's ratios and E 1 and E 2 are the elastic moduli for the elastic bodies denoted by subscripts 1 and 2, respectively, in Figure 5. The distribution of the contact pressure acting on the elastic body p is expressed as a function of the distance r from the center of the circle at the contact area.
where p 0 is the maximum contact pressure, which is applied at the center of the contact surface, as reported by Johnson [31].
The radius of the contact surface is given by In this study, elastic bodies 1 and 2 are the CrN coating of the pin and the DP980 steel sheet, respectively. Ichimura and Ando [32] measured the elastic modulus of CrN coating (E 1 ), which can be determined to be 375 GPa. Poisson's ratio of the CrN coating (ν 1 ) is selected as 0.28, obtained from Fabis et al. [33]. The elastic modulus of the DP980 steel sheet (E 2 ) and Poisson's ratio of the DP980 steel sheet (ν 2 ) are 210 GPa and 0.3, respectively, which are the general elastic properties of the steel sheet. As shown in Figure 6, an analytical solution for the normal force is calculated as a function of the contact pressure by using the Hertz contact theory. From this result, the normal force equivalent to the contact pressure of 1.85 GPa was calculated as 9.50 N for the given geometry of the pin. Therefore, a maximum normal force of 10.0 N was selected for the pin-on-disk test.
where p0 is the maximum contact pressure, which is applied at the center of the co surface, as reported by Johnson [31].

= 3 2
The radius of the contact surface is given by In this study, elastic bodies 1 and 2 are the CrN coating of the pin and the DP980 sheet, respectively. Ichimura and Ando [32] measured the elastic modulus of CrN co (E1), which can be determined to be 375 GPa. Poisson's ratio of the CrN coating ( selected as 0.28, obtained from Fabis et al. [33]. The elastic modulus of the DP980 sheet (E2) and Poisson's ratio of the DP980 steel sheet (ν2) are 210 GPa and 0.3, respecti which are the general elastic properties of the steel sheet. As shown in Figure 6, an lytical solution for the normal force is calculated as a function of the contact pressu using the Hertz contact theory. From this result, the normal force equivalent to the co pressure of 1.85 GPa was calculated as 9.50 N for the given geometry of the pin. There a maximum normal force of 10.0 N was selected for the pin-on-disk test.
The mechanical press is generally used for the sheet metal forming process of t pillar. The conventional stamping speed is approximately 10 spm in the B-pillar mass duction process, and the maximum forming velocity is approximately 180 mm/s, co ering the sliding motion of the sheet metal. Therefore, the maximum sliding velocity set to 188.4 mm/s (120 rpm), with a rotation radius of 15 mm. The sliding distance selected as 1000 revolution (94.20 m) for sufficient tool wear to occur. Table 3 show levels of the three contact parameters based on the above maximum test values. F these values, an orthogonal array L25 (5 3 ) can be constructed as shown in Table 4.  The mechanical press is generally used for the sheet metal forming process of the B-pillar. The conventional stamping speed is approximately 10 spm in the B-pillar mass production process, and the maximum forming velocity is approximately 180 mm/s, considering the sliding motion of the sheet metal. Therefore, the maximum sliding velocity was set to 188.4 mm/s (120 rpm), with a rotation radius of 15 mm. The sliding distance was selected as 1000 revolution (94.20 m) for sufficient tool wear to occur. Table 3 shows the levels of the three contact parameters based on the above maximum test values. From these values, an orthogonal array L 25 (5 3 ) can be constructed as shown in Table 4.

Measurement of Pin Wear Depth
The pin-on-disk test was performed in various experimental conditions based on the orthogonal array constructed in the previous section. Three repetitive tests were performed for all design cases. After the pin-on-disk tests, the tool wear areas were measured with an optical electron microscope (OEM) ECLIPSE MA200 (Nikon, Tokyo, Japan), as shown in Figure 7. From the tool wear area, the wear scar diameter d w can be obtained by assuming the measured wear area of the pin as a circle standardized in ASTM (American Society for Testing and Materials) G99-17 [34].
Then, the tool wear depth h is related to the wear scar diameter d w as follows: where r p is the pin end radius. By using the following equation, the pin volume losses V loss can be calculated from the wear scar diameter d w and the tool wear depth h derived from the previous equations.
In Table 4, the measured wear depth h and pin volume loss V loss are listed for all experimental design cases. In general, the coating thickness of CrN for the stamping die is approximately 10 µm. The test range of the pin-on-disk is well designed to evaluate the wear amount considering the CrN coating thickness.
Metals 2021, 11, x FOR PEER REVIEW 9 of 20 By using the following equation, the pin volume losses Vloss can be calculated from the wear scar diameter dw and the tool wear depth h derived from the previous equations.
In Table 4, the measured wear depth h and pin volume loss Vloss are listed for all experimental design cases. In general, the coating thickness of CrN for the stamping die is approximately 10 μm. The test range of the pin-on-disk is well designed to evaluate the wear amount considering the CrN coating thickness.

Sensitivity Analysis of Wear Parameters
The S/N ratio and the ANOVA result were used to check the significant control factors affecting the wear characteristics of the tool. The S/N ratio was used to determine the robustness of the control factors of the product or process. It measures how the response changes relative to the target values under various noise conditions. The performance characteristics can be quantified in terms of the respective response to noise and signal factors. During the response analysis, the deviation between the experimental value and the desired value can be computed using the loss function and converted into the S/N ratio. There are three kinds of S/N ratios for analyzing and identifying performance characteristics, depending on the objective of the experiments: higher-the-better, lower-thebetter, and nominal-the-best [35]. Among them, the higher-the-better ratio is suitable for evaluating the characteristics of the tool wear problem because tool wear shows an increasing tendency according to the contact parameters. The S/N ratio (η) is calculated as follows: where yi is the data measured in the i-th experiment and n is the total number of experiments. As you can see from Figure 8, both pin wear depth and pin volume loss show variation in the higher-the-better ratio with respect to the contact parameters. From this figure, it was confirmed that the sliding velocity might not be a significant factor in the current tool wear during the investigated forming process. The ANOVA is another sensitivity analysis method to evaluate the statistical significance of independent variables that generate variance in dependent variables. It is used to identify the relationship between independent and dependent variables. Table 5 shows the ANOVA results of the pin wear depth and the pin volume loss according to the three contact parameters. The F-value represents the significance level of the control factors. The p-value can also be used to check the significance of control factors with the criterion of 0.05 for the 95% confidence level. This table also indicates that the sliding velocity is a statistically insignificant factor for both pin wear depth and pin volume loss, as the F-

Sensitivity Analysis of Wear Parameters
The S/N ratio and the ANOVA result were used to check the significant control factors affecting the wear characteristics of the tool. The S/N ratio was used to determine the robustness of the control factors of the product or process. It measures how the response changes relative to the target values under various noise conditions. The performance characteristics can be quantified in terms of the respective response to noise and signal factors. During the response analysis, the deviation between the experimental value and the desired value can be computed using the loss function and converted into the S/N ratio. There are three kinds of S/N ratios for analyzing and identifying performance characteristics, depending on the objective of the experiments: higher-the-better, lowerthe-better, and nominal-the-best [35]. Among them, the higher-the-better ratio is suitable for evaluating the characteristics of the tool wear problem because tool wear shows an increasing tendency according to the contact parameters. The S/N ratio (η) is calculated as follows: where y i is the data measured in the i-th experiment and n is the total number of experiments. As you can see from Figure 8, both pin wear depth and pin volume loss show variation in the higher-the-better ratio with respect to the contact parameters. From this figure, it was confirmed that the sliding velocity might not be a significant factor in the current tool wear during the investigated forming process. The ANOVA is another sensitivity analysis method to evaluate the statistical significance of independent variables that generate variance in dependent variables. It is used to identify the relationship between independent and dependent variables. Table 5 shows the ANOVA results of the pin wear depth and the pin volume loss according to the three contact parameters. The F-value represents the significance level of the control factors. The p-value can also be used to check the significance of control factors with the criterion of 0.05 for the 95% confidence level. This table also indicates that the sliding velocity is a statistically insignificant factor for both pin wear depth and pin volume loss, as the F-value of the sliding velocity is relatively much smaller than that of other contact parameters.
Moreover, the p-value is also larger than 0.05. Therefore, it was concluded that the sliding velocity can be excluded for developing the prediction model because of its negligible effect on wear behavior. value of the sliding velocity is relatively much smaller than that of other contact parameters. Moreover, the p-value is also larger than 0.05. Therefore, it was concluded that the sliding velocity can be excluded for developing the prediction model because of its negligible effect on wear behavior.

Tool Wear Prediction Model
The tool wear prediction model was constructed by approximating the pin-on-disk test results for various design cases generated by the Taguchi method. The normal force was converted again into contact pressure because the measured value for tool wear prediction in the numerical simulation is the local contact pressure. After excluding the sliding velocity based on the statistical sensitivity analysis described in the Section 3.4, the tool wear depth was formulated by the power law equation as follows: where P and L are the contact pressure and the sliding distance, respectively, and K, α, and β are model constants. From the fitting of the proposed power law model to the experimental data, the exponents of the contact pressure and the sliding distance are α = 1.7570 and β = 0.5073, respectively, and the parameter K = 0.01661 for the investigated DP980 steel and tool material. The identified model parameters show that the contact pressure is a more sensitive parameter to the tool wear than the sliding distance. To verify the reliability of the tool wear prediction model, the comparison between the measured and predicted wear depths was examined for each run of the orthogonal array L25 (5 3 ) as shown in Figure 9. The least-squares R 2 value and R 2 Adj. value in the estimation of the tool wear depth were found to be 87.01% and 86.15%, respectively, which indicates a high level

Tool Wear Prediction Model
The tool wear prediction model was constructed by approximating the pin-on-disk test results for various design cases generated by the Taguchi method. The normal force was converted again into contact pressure because the measured value for tool wear prediction in the numerical simulation is the local contact pressure. After excluding the sliding velocity based on the statistical sensitivity analysis described in the Section 3.4, the tool wear depth was formulated by the power law equation as follows: where P and L are the contact pressure and the sliding distance, respectively, and K, α, and β are model constants. From the fitting of the proposed power law model to the experimental data, the exponents of the contact pressure and the sliding distance are α = 1.7570 and β = 0.5073, respectively, and the parameter K = 0.01661 for the investigated DP980 steel and tool material. The identified model parameters show that the contact pressure is a more sensitive parameter to the tool wear than the sliding distance. To verify the reliability of the tool wear prediction model, the comparison between the measured and predicted wear depths was examined for each run of the orthogonal array L 25 (5 3 ) as shown in Figure 9. The least-squares R 2 value and R 2 Adj. value in the estimation of the tool wear depth were found to be 87.01% and 86.15%, respectively, which indicates a high level of linear relationship between the measured and predicted results. Thus, the power law equation is a suitable formulation to develop the tool wear prediction model. of linear relationship between the measured and predicted results. Thus, the power law equation is a suitable formulation to develop the tool wear prediction model.

Verification of Tool Wear Using the U-Draw Forming Process
The U-draw bending test was adopted to evaluate the wear amount during the drawing process. A modified punch shape was designed by numerical simulation to cause contact pressure localization and accelerate tool wear. Based on the conventional and modified tool shapes, the drawing wear tests were performed to investigate the actual wear behavior according to the die contact pressure and the sliding distance. Finally, the tool wear prediction model was verified by comparing the tool wear depth obtained from the experiments in terms of the angle of the die radius.

Drawing Tool Design by Numerical Simulation
A numerical simulation was performed for the conventional U-draw bending test to investigate the contact conditions during the test. The simulation was performed using PAM-STAMP v2017.0 (ESI-Group, Rungis, France). For simulation efficiency, a half model was constructed by applying the symmetric conditions of the tool and material properties. The tool was modeled with rigid elements because the elastic deformation of the die has little effect on the contact pressure, as reported by Ersoy-Nurnberg et al. [23]. The blank was discretized by deformable solid elements with three layers. As shown in Figure 10, fine meshes were applied in the contact interface between the tool and the blank to predict the contact pressure more accurately at the interface by excluding the mesh size effect. However, finer meshes increase the computation time. After optimizing the mesh sizes of the simulations, the determined element sizes of the tool and the blank were 750,000 and 53,988, respectively. The same blank material from the previous section was used again, and its mechanical properties are summarized in Table 1.
The friction coefficient between the tool and the blank was selected as 0.15, which achieved good correlation between the experiment and the numerical simulation, reported by Pereira et al. [22]. In the numerical simulation, the contact pressure distribution was observed during a punch stroke of 40 mm at the die shoulder where a severe contact pressure occurred. Figure 11a shows the contact pressure distribution according to the punch stroke along the die shoulder. Figure 11b,c shows the maximum contact pressure and its location with respect to the punch stroke. For the conventional die, the maximum contact pressure stabilizes after the 12 mm punch stroke. The peak and stabilized values were 0.940 and 0.502 GPa, respectively. From this result, it can be noted that the conventional drawing test has a limitation for describing the severe contact condition of the

Verification of Tool Wear Using the U-Draw Forming Process
The U-draw bending test was adopted to evaluate the wear amount during the drawing process. A modified punch shape was designed by numerical simulation to cause contact pressure localization and accelerate tool wear. Based on the conventional and modified tool shapes, the drawing wear tests were performed to investigate the actual wear behavior according to the die contact pressure and the sliding distance. Finally, the tool wear prediction model was verified by comparing the tool wear depth obtained from the experiments in terms of the angle of the die radius.

Drawing Tool Design by Numerical Simulation
A numerical simulation was performed for the conventional U-draw bending test to investigate the contact conditions during the test. The simulation was performed using PAM-STAMP v2017.0 (ESI-Group, Rungis, France). For simulation efficiency, a half model was constructed by applying the symmetric conditions of the tool and material properties. The tool was modeled with rigid elements because the elastic deformation of the die has little effect on the contact pressure, as reported by Ersoy-Nurnberg et al. [23]. The blank was discretized by deformable solid elements with three layers. As shown in Figure 10, fine meshes were applied in the contact interface between the tool and the blank to predict the contact pressure more accurately at the interface by excluding the mesh size effect. However, finer meshes increase the computation time. After optimizing the mesh sizes of the simulations, the determined element sizes of the tool and the blank were 750,000 and 53,988, respectively. The same blank material from the previous section was used again, and its mechanical properties are summarized in Table 1.
The friction coefficient between the tool and the blank was selected as 0.15, which achieved good correlation between the experiment and the numerical simulation, reported by Pereira et al. [22]. In the numerical simulation, the contact pressure distribution was observed during a punch stroke of 40 mm at the die shoulder where a severe contact pressure occurred. Figure 11a shows the contact pressure distribution according to the punch stroke along the die shoulder. Figure 11b,c shows the maximum contact pressure and its location with respect to the punch stroke. For the conventional die, the maximum contact pressure stabilizes after the 12 mm punch stroke. The peak and stabilized values were 0.940 and 0.502 GPa, respectively. From this result, it can be noted that the conventional drawing test has a limitation for describing the severe contact condition of the DP980 for a B-pillar. Therefore, a modified die design is required to generate a higher contact pressure. DP980 for a B-pillar. Therefore, a modified die design is required to generate a higher contact pressure. A new die design was proposed for generating a higher contact pressure than the conventional U-draw bending test. To introduce a more severe contact condition, a small bead with a height of 0.1 mm was added at the die angle of 59°, where the peak contact pressure of the conventional die occurs. A similar level of contact pressure was generated at the start regions of the die radius. Then, the modified die resulted in a higher contact pressure at the bead position, unlike the conventional die. The maximum contact pressure of the modified die was 1.866 GPa, as presented in Figure 11b. This value is similar to the contact pressure investigated in the forming simulation of the B-pillar. Figure 11c also confirms the position of the maximum contact pressure after stabilization. The simulation results indicate that the modified die design is effective for describing the severe contact behavior of UHSS.

Drawing Wear Test
The drawing wear test was performed using the conventional and modified U-bending dies proposed by the numerical simulation. The geometric and process parameters DP980 for a B-pillar. Therefore, a modified die design is required to generate a higher contact pressure. A new die design was proposed for generating a higher contact pressure than the conventional U-draw bending test. To introduce a more severe contact condition, a small bead with a height of 0.1 mm was added at the die angle of 59°, where the peak contact pressure of the conventional die occurs. A similar level of contact pressure was generated at the start regions of the die radius. Then, the modified die resulted in a higher contact pressure at the bead position, unlike the conventional die. The maximum contact pressure of the modified die was 1.866 GPa, as presented in Figure 11b. This value is similar to the contact pressure investigated in the forming simulation of the B-pillar. Figure 11c also confirms the position of the maximum contact pressure after stabilization. The simulation results indicate that the modified die design is effective for describing the severe contact behavior of UHSS.

Drawing Wear Test
The drawing wear test was performed using the conventional and modified U-bending dies proposed by the numerical simulation. The geometric and process parameters represent a typical sheet metal stamping operation that is commonly prone to wear. The A new die design was proposed for generating a higher contact pressure than the conventional U-draw bending test. To introduce a more severe contact condition, a small bead with a height of 0.1 mm was added at the die angle of 59 • , where the peak contact pressure of the conventional die occurs. A similar level of contact pressure was generated at the start regions of the die radius. Then, the modified die resulted in a higher contact pressure at the bead position, unlike the conventional die. The maximum contact pressure of the modified die was 1.866 GPa, as presented in Figure 11b. This value is similar to the contact pressure investigated in the forming simulation of the B-pillar. Figure 11c also confirms the position of the maximum contact pressure after stabilization. The simulation results indicate that the modified die design is effective for describing the severe contact behavior of UHSS.

Drawing Wear Test
The drawing wear test was performed using the conventional and modified U-bending dies proposed by the numerical simulation. The geometric and process parameters represent a typical sheet metal stamping operation that is commonly prone to wear. The test configuration is summarized in Table 6, which was reported from Weiss et al. [36]. Figure 12 shows the experimental setup and schematic diagrams of the drawing wear test. It was designed as an insert-type tool for efficient measurement of the wear profile on the tool surface. The tool was made using STD11 tool steel and heat-treated under the same conditions as the pin of the pin-on-disk tester, as shown in Figure 3. The tool surface was coated by the PVD method with CrN to improve the wear resistance. The coating thickness was measured using the OEM ECLIPSE MA200 (Nikon, Tokyo, Japan), and it was approximately 9.9 µm, as shown in Figure 13. The Vickers hardness HV 0.1 of the PVD CrN-coated layer was 2137.6 ± 35.7, measured at a normal force of 0.1 kgf. test configuration is summarized in Table 6, which was reported from Wei Figure 12 shows the experimental setup and schematic diagrams of the drawi It was designed as an insert-type tool for efficient measurement of the wear p tool surface. The tool was made using STD11 tool steel and heat-treated und conditions as the pin of the pin-on-disk tester, as shown in Figure 3. The tool coated by the PVD method with CrN to improve the wear resistance. The co ness was measured using the OEM ECLIPSE MA200 (Nikon, Tokyo, Japan) approximately 9.9 μm, as shown in Figure 13. The Vickers hardness HV0.1 of th coated layer was 2137.6 ± 35.7, measured at a normal force of 0.1 kgf.   The experiments with the conventional die performed 5000 CrN-coated layer. The modified die only required 2500 hits beca pressure. The tool surface at the centerline of the die shoulder was ECLIPSE MA200 (Nikon, Tokyo, Japan) to confirm the wear state observed segmentally at an interval of 10° in the area of interest. A to magnify the tool surface. The individual photographs were atta complete view of the die radius surface. Figure 14 shows the scratches captured on the die radius surface after the wear test. For scratches (or defects on the surface) mostly occurred at the start o 0.50 GPa contact pressure was applied on the die radius of 7.5° for in Figure 11. In the case of the modified die, it was found that th pletely worn out, and the boundary between the coating layer an be observed at the die angle of 59.0°, where a pressure of 1.866 G time, as shown in Figure 11. The worn surface of the modified die was further observed by OEM. Microscratches were observed at were parallel to the sliding direction, as shown in Figure 15. It wear in the sheet metal forming process was generally due to abr results confirm again that the contact pressure and the sliding parameters that have a significant influence on the tool wear beh The experiments with the conventional die performed 5000 hits to wear out the PVD CrN-coated layer. The modified die only required 2500 hits because of the higher contact pressure. The tool surface at the centerline of the die shoulder was observed using the OEM ECLIPSE MA200 (Nikon, Tokyo, Japan) to confirm the wear state of the tool surface. It was observed segmentally at an interval of 10 • in the area of interest. A 10× optical lens was used to magnify the tool surface. The individual photographs were attached together to create a complete view of the die radius surface. Figure 14 shows the amount and locations of scratches captured on the die radius surface after the wear test. For the conventional die, the scratches (or defects on the surface) mostly occurred at the start of the die radius, where a 0.50 GPa contact pressure was applied on the die radius of 7.5 • for a long distance, as shown in Figure 11. In the case of the modified die, it was found that the coating layer was completely worn out, and the boundary between the coating layer and the tool material could be observed at the die angle of 59.0 • , where a pressure of 1.866 GPa was applied for a long time, as shown in Figure 11. The worn surface of the modified die at the die angle of 59.0 • was further observed by OEM. Microscratches were observed at the worn surface, which were parallel to the sliding direction, as shown in Figure 15. It was found that the tool wear in the sheet metal forming process was generally due to abrasive wear. The ovserved results confirm again that the contact pressure and the sliding distance are the contact parameters that have a significant influence on the tool wear behavior.
The wear depths for the conventional die and the modified die were quantitatively measured to estimate the wear behaviors of the respective tool surfaces. The tool wear behavior was evaluated by analyzing the relationship between the measured wear depth and the contact condition of the respective tool surface. The tool wear depth was measured using Leitz Infinity 12.10.6 (Hexagon AB, Stockholm, Sweden), which is a contact-type three-dimensional measuring instrument capable of recognizing relative coordinates with a precision of 0.3 + L/1000 [µm] (L = mm). First, the relative coordinates were recognized in the area where the tool wear did not occur. Next, the wear depth was measured by the absolute difference of the surface profile at the wear occurrence position before and after the wear test. As shown in Figure 16a, the relative coordinate was recognized at the blue dashed line where there was no contact with the blank. Subsequently, the origin of the coordinate was generated. Based on the origin of the coordinate, the tool surface profile was recognized by measuring the solid red line, which is the centerline of the tool, in contact with the blank during the drawing process. As shown in Figure 16b, the point data were measured along with the profile of the tool surface and each point was connected using the spline option in AUTOCAD 2020 (Autodesk, Inc., San-Rafael, CA, USA). The maximum wear depth of the conventional die was 2.25 µm at a die angle of 17 • . For the modified die, the wear depth drastically increased up to 10.42 µm because of the high contact pressure. The wear location also moved to the die angle at approximately 59 • . From this result, the drawing wear test condition was established successfully to induce the high contact pressure predicted by the numerical simulation.
scratches (or defects on the surface) mostly occurred at the start of the die radius, wher 0.50 GPa contact pressure was applied on the die radius of 7.5° for a long distance, as sho in Figure 11. In the case of the modified die, it was found that the coating layer was co pletely worn out, and the boundary between the coating layer and the tool material co be observed at the die angle of 59.0°, where a pressure of 1.866 GPa was applied for a lo time, as shown in Figure 11. The worn surface of the modified die at the die angle of 59 was further observed by OEM. Microscratches were observed at the worn surface, wh were parallel to the sliding direction, as shown in Figure 15. It was found that the t wear in the sheet metal forming process was generally due to abrasive wear. The ovserv results confirm again that the contact pressure and the sliding distance are the cont parameters that have a significant influence on the tool wear behavior.  The wear depths for the conventional die and the modified measured to estimate the wear behaviors of the respective too behavior was evaluated by analyzing the relationship between t were measured along with the profile of the tool surface and each point was connected using the spline option in AUTOCAD 2020 (Autodesk, Inc., San-Rafael, CA, USA). The maximum wear depth of the conventional die was 2.25 μm at a die angle of 17°. For the modified die, the wear depth drastically increased up to 10.42 μm because of the high contact pressure. The wear location also moved to the die angle at approximately 59°. From this result, the drawing wear test condition was established successfully to induce the high contact pressure predicted by the numerical simulation.

Validation of the Tool Wear Prediction Procedure
The validation process was performed to verify the reliability of the tool wear prediction model and its identified parameters. The experimental and predicted wear depth profiles were compared to evaluate the accuracy of the tool wear prediction procedure. The sliding distance and the contact pressure were investigated to predict the wear depth using the prediction model. The numerical simulation method can be used for investigating the sliding distance and the contact pressure along the die radius. Figure 17a shows the sliding distance with respect to the die radius angle. The decrease in the sliding distance was caused by the sequential onset of the contact between the die and the blank as the punch stroke went upward. Figure 17b shows the contact pressure-sliding distance curves for the two die designs and die angles. By applying the contact pressure-sliding distance histories to the tool wear prediction model, a continuous wear depth profile can be obtained, as shown in Figure 18. For the conventional and modified dies, the maximum wear depths by the prediction model were 2.56 and 10.88 μm, respectively. For the conventional die, the predicted value showed a difference of 0.31 μm from the experimental value, which is an error of 13.78%. For the modified die, the difference between the predicted value and the experimental value was 0.46 μm, i.e., an error of 4.41%. The locations of the maximum wear depth were similar between the experimental and predicted profiles. For comparison, the commonly used Archard wear model was used to predict the wear behavior of the drawing wear tests. For the conventional and modified tools, the maximum predicted wear depths were 13.01 and 8.48 μm, respectively, as shown in Figure 18. Then, the prediction errors of the maximum wear depth were 478.22 and 18.62%, respectively. In the case of the conventional die, the predicted wear depth at a die angle of 7.5° by the Archard model was significantly over-predicted compared to that of the power law prediction model. This is because the predicted die angle of 7.5° has a relatively long sliding distance (Figure 17b) and the Archard wear model is directly proportional to the sliding distance, while the exponent term for the sliding distance is 0.0573 (<1.0) in the

Validation of the Tool Wear Prediction Procedure
The validation process was performed to verify the reliability of the tool wear prediction model and its identified parameters. The experimental and predicted wear depth profiles were compared to evaluate the accuracy of the tool wear prediction procedure. The sliding distance and the contact pressure were investigated to predict the wear depth using the prediction model. The numerical simulation method can be used for investigating the sliding distance and the contact pressure along the die radius. Figure 17a shows the sliding distance with respect to the die radius angle. The decrease in the sliding distance was caused by the sequential onset of the contact between the die and the blank as the punch stroke went upward. Figure 17b shows the contact pressure-sliding distance curves for the two die designs and die angles. By applying the contact pressure-sliding distance histories to the tool wear prediction model, a continuous wear depth profile can be obtained, as shown in Figure 18. For the conventional and modified dies, the maximum wear depths by the prediction model were 2.56 and 10.88 µm, respectively. For the conventional die, the predicted value showed a difference of 0.31 µm from the experimental value, which is an error of 13.78%. For the modified die, the difference between the predicted value and the experimental value was 0.46 µm, i.e., an error of 4.41%. The locations of the maximum wear depth were similar between the experimental and predicted profiles. For comparison, the commonly used Archard wear model was used to predict the wear behavior of the drawing wear tests. For the conventional and modified tools, the maximum predicted wear depths were 13.01 and 8.48 µm, respectively, as shown in Figure 18. Then, the prediction errors of the maximum wear depth were 478.22 and 18.62%, respectively. In the case of the conventional die, the predicted wear depth at a die angle of 7.5 • by the Archard model was significantly over-predicted compared to that of the power law prediction model. This is because the predicted die angle of 7.5 • has a relatively long sliding distance (Figure 17b) and the Archard wear model is directly proportional to the sliding distance, while the exponent term for the sliding distance is 0.0573 (<1.0) in the proposed power law model. For the modified tool die, the predicted wear depth at the die angle of 59.0 • by the Archard model was quite under-predicted compared to that of the power law prediction model. This is attributed to the high contact pressure at the die angle of 59.0 • (Figure 17b). The Archard wear model has a linear relationship with the contact pressure, whereas the proposed model is a nonlinear function with an exponent of 1.7570 (>1.0) for the contact pressure. Therefore, the Archard wear model is proven to have a limitation in predicting the nonlinear behavior of wear in the sheet forming process of DP980 steel. On the other hand, the power-law-type wear model provides more reliable results than the linear model for predicting the wear behavior of DP980 sheets under realistic sheet forming processes.

Conclusions
This study proposes a prediction method for abrasive tool wear behavior in the drawing process using a DP980 steel sheet by identifying the relationship between the wear depth and the contact conditions on the tool surface.
(1) The contact conditions applied to the die surface were evaluated by the sheet metal forming of the B-pillar based on FE simulation. The contact pressure 1.852 GPa was concentrated at the die corner radius where the actual wear occurred. It was conformed that the contact pressure has a significant influence on the tool wear, by comparing the wear status of the B-pillar die and simulation results. (2) The pin-on-disk test was conducted for the quantitative evaluation of the wear behavior. Based on the contact pressure of 1.852 GPa of the B-pillar simulation result,

Conclusions
This study proposes a prediction method for abrasive tool wear behavior in the drawing process using a DP980 steel sheet by identifying the relationship between the wear depth and the contact conditions on the tool surface.
(1) The contact conditions applied to the die surface were evaluated by the sheet metal forming of the B-pillar based on FE simulation. The contact pressure 1.852 GPa was concentrated at the die corner radius where the actual wear occurred. It was conformed that the contact pressure has a significant influence on the tool wear, by comparing the wear status of the B-pillar die and simulation results. (2) The pin-on-disk test was conducted for the quantitative evaluation of the wear behavior. Based on the contact pressure of 1.852 GPa of the B-pillar simulation result,

Conclusions
This study proposes a prediction method for abrasive tool wear behavior in the drawing process using a DP980 steel sheet by identifying the relationship between the wear depth and the contact conditions on the tool surface.
(1) The contact conditions applied to the die surface were evaluated by the sheet metal forming of the B-pillar based on FE simulation. The contact pressure 1.852 GPa was concentrated at the die corner radius where the actual wear occurred. It was conformed that the contact pressure has a significant influence on the tool wear, by comparing the wear status of the B-pillar die and simulation results. (2) The pin-on-disk test was conducted for the quantitative evaluation of the wear behavior. Based on the contact pressure of 1.852 GPa of the B-pillar simulation result, the range of test conditions was established. According to statistical analyses using the S/N ratio and the ANOVA result, it was found that contact pressure and sliding distance had a significant influence on the tool wear of the drawing process but the effect of the sliding velocity was minor. (3) The tool wear prediction model was constructed by the power-law-type regression model to predict the nonlinear evolution of the tool wear depth as a function of the contact pressure and the sliding distance based on the pin-on-disk test results. The U-draw bending test was employed to verify the reliability of the wear prediction model. The modified die was newly designed to accelerate the wear behaviors of the tool surface, which generated a contact pressure four times higher than that of the conventional die. The drawing wear test was conducted using conventional and accelerated drawing tools to measure the wear depth under different contact conditions. For the conventional die and the modified die, the predicted values showed an error of 13.78 and 4.41%, respectively, compared with the experimental values. From the validation process, it can be concluded that the power-law-type wear prediction model could provide reliable prediction for the continuous wear depth profile of the DP980 material and the PVD CrN-coated tool. (4) The suggested prediction procedure can be further generalized to obtain tool wear predictions by constructing an experimental database of the pin-on-disk tests of combinations of various sheet materials and tool coatings. It has the advantage of saving time and cost because it can predict the tool wear in the drawing process by using pin-on-disk test results. This prediction procedure will help reduce maintenance costs and productivity. Funding: This study has been conducted with the support of the Korea Institute of Industrial Technology as "Development of intelligent root technology with add-on modules (KITECH EO-21-0009)" and "WC300 R&D Project (S2482274)".

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.

Acknowledgments:
The authors appreciate the support of POSCO Global R&D Center for providing the test material and the research equipment. η Signal-to-noise (S/N) ratio K, α, β Coefficient of the wear depth prediction model