Mechanical Behavior and Damage of Zinc Coating for Hot Dip Galvanized Steel Sheet DP600

Advanced high strength galvanized steel sheet has been one of the dominant materials of modern automotive panels because of its outstanding mechanical properties and corrosion resistance. The zinc coating thickness of hot dip galvanized steel sheet is only about 10–20 μm, which is a discarded object on the macro level. However, it is obvious to damage and impact on stamping performance. Therefore, this paper takes zinc coating as the research object and builds its mechanical constitutive model based on a nano-indentation test and dimensional analysis theory. We separated the zinc coating from the galvanized steel substrate and constructed a sandwich material model by introducing a cohesive layer to connect the zinc coating and the steel substrate. We obtained the interface binding energy between the zinc coating and the steel substrate through the nano-scratch test. The accuracy of the model is verified by the finite element analysis of hemispherical parts. We used the five-layers element model with 0 thickness cohesive layer to simulate the zinc coating damage of galvanized steel sheet. The hemispherical part drawing experiment is used to verify the feasibility of the finite element analysis results. The results demonstrate that it is more accurate to consider the finite element numerical simulation of the zinc coating, introducing the cohesive element to simulate damage between the coating and the substrate. Drawing depth, stamping force, and the strain of the numerical simulation are closer to the experimental results.


Introduction
Galvanized advanced high strength steel (AHSS) has become the dominant material of automotive lightweight because of its superior strength, elongation, and outstanding corrosion resistance, which can follow the requirements of collision safety, energy saving, emission reduction, and environmental protection [1]. The coating is an ultra-thin material galvanized on the upper and lower surfaces of the steel substrate.. There are two galvanizing methods: hot-dip galvanizing and electro galvanizing [2]. The thickness of the zinc coating is about 8 µm [3], which represents a slight difference between galvanizing methods. It neglects this level of magnitude on macroscopic scales. But zinc coating layer makes the surface characteristics of galvanized steel different from that of nongalvanized steel. These zinc coating layers are easy to damage in severe plastic deformation under external pressure and sliding friction [4][5][6]. It will also affect the formation of welding joints of stamping parts [7]. Typical damages of the coating layers are cracking and powdering flaking. The exfoliated coating adheres to the surfaces of the steel and die, resulting in falling, wearing, and affecting the stamping formability of the steel sheet [8,9]. Therefore, it is of great practical significance to study the mechanical properties of the zinc coating and its influence on the deformation of galvanized AHSS during stamping forming.
Zinc plating is an alloying process, and its element content will influence the performance of the galvanized layer. Vourlias et al. [10] studied the effect of alloying elements on the interfacial reaction of the galvanized coating. They found that the concentration and distribution of alloying elements played an important role in the phase growth. To improve the coating quality of high strength DP steel, a nickel layer was pre-electroplated on the steel substrate before galvanizing, and an enhancement in the coating quality was discovered [11]. Parisot et al. [12] confirmed the damage mechanisms at work in three different microstructures of a zinc coating on an interstitial-free steel substrate under tension, plain strain, and expansion loading. Two main fracture mechanisms, intergranular fracture and transgranular boundary fracture, were created in cold-rolled coatings, tempered cold-rolled coatings, and uncrystallized coatings. Further investigation shows that the damage between the galvanized layer and substrate depends on the different mechanical behavior and thickness of intermetallic compounds [13]. Petit et al. [14] investigated the galvanization ability of a high strength chromium-rich TRIP800 (transformation induced plasticity steel). Under different loads and deformation conditions, zinc coating also shows different mechanics. The damage behavior of different galvanized steel sheets under tension-bending conditions was explored by Yu et al. [15]. They discovered that the exfoliating and scratching of coating are two types of surface damage. To quantify the damage of galvanized coating, Xu et al. [16] suggested a parameter to illustrate the powdering of galvanized coating and its validity by double reverse Olsen test. Song and Sloof [17] introduced a "macroscopic atomic model" to assess the adhesion of the zinc coating for hot-dip galvanized dual phase steel. Adhesion force is an important property of coated steel, which concluded the quality of the coating. The damage of the coating in the stamping process results from the relative sliding friction between the steel sheet and the die surfaces. Wang et al. [9] considered the evolutionary process of galling of galvanized AHSS DP590 (dual-phase steel) under the condition of friction coupling plastic bending deformation. First, the coating is cracked because of tensile deformation. Then, further powdering of the coating resulted in a notable increase in the surface roughness of the contact interface. Die wear and surface galling of stamping parts will be caused because of the adhesion effect of coating. Therefore, the coating layer must be considered as an essential research branch of AHSS stamping.
Numerical simulation is an effective tool and method to study the strain distribution stamping forming [18]. As a non-negligible factor, we must consider the coating in the numerical simulation of stamping forming. The upper and lower surface of the galvanized steel sheet is zinc coating, and the intermediate substrate is a steel sheet. We can understand it as three-layers of composite material, like a sandwich equivalent to the stamping of multiple layers of material. It includes multi-scale modeling of materials [19], stamping formability [20][21][22], spring-back [23,24], and evolution of interface structure [25]. A three-dimensional finite element (FE) simulation of a thin zinc coating on a galvanized steel sheet has been performed by taking the multicrystalline structure of the coating into account [26]. In addition, Ucun and Aslantas [27] simulated the mechanical properties and force status of the cutting tools with different coating types by the method of single layer and multi-layer finite element modeling. Lee et al. [28] simulated the deep drawing of galvanized steel with a multilayer element model, and the interface between substrate and coating was set to be adherent. Kim et al. [29] studied the fracture of the coating layer of galvanized steel using V-bending test and numerical simulation of the solid element. The results show that the coating and substrate begin to failure when the local plastic equivalent plastic reaches 0.26-0.28. Bettaieb et al. [30] proposed a solid shell element to describe the coating under T-bending based on the enhanced assumed strain (EAS) technique and the assumed natural strain (ANS) technique, that can avoid locking problems caused by the large aspect ratio of ultra-thin coating elements. It is shown that a higher accuracy simulation of coating can be achieved by using such a solid shell element.
Although the coating is taken into consideration in the above studies, the element model of numerical simulation still cannot describe the deformation model of the coating and substrate of galvanized AHSS. Therefore, in this paper, a coupled element model based on solid shell element and shell element is introduced to simulate the deformation of galvanized AHSS sheet, which recognized the scale effect of the coating. The mechanical properties of the coating and substrate were obtained by nano-indentation and uniaxial tensile test as the input condition of numerical simulation. It employs the dynamic friction coefficient considering the positive pressure between the stamping interfaces. The stamping forming experiment is developed and implemented, and we compare the numerical simulation results are compared with the experimental results. It is used to prove the numerical simulation accuracy of the coupled element model.

Mechanical Behavior of Zinc Coating
The hot dip galvanized cold rolled steel DP600 with a thickness 1.2 mm was chosen as the research object, which is the most applied steel in the galvanized advanced high strength steel sheet. The composition of the steel substrate is C-0.13%, Si-0.36%, Mn-0.168%, Alt-0.048%, and the composition of the zinc coating is mainly expressed by the thickness of the coating, as shown in the following experimental results. The stress-strain curves of the substrate and the coating layer of the steel sheet will be obtained by the experimental method. The detailed process of the experiment is as follows.

Constitutive Model of Zinc Coating
By inspecting the cross-section electron microscopic (SEM) results of many specimens, the thickness of the zinc layer is between 14.5-15.3 µm. Therefore, the thickness of DP600 zinc coating is set to 15 µm, that is to say, the upper and lower surfaces of the substrate have a zinc coating with a thickness of 15 µm, which is like a sandwich material. SEM measured the cross-section morphology of the specimens displayed in Figure 1, and the thickness of the zinc coating is 15.1 µm. The uniaxial tension test is used to get the mechanical properties of the steel substrate. To study the mechanical properties of the steel substrate, it is necessary to remove the zinc coating. We analyzed the surface morphology of a steel sheet before and after acid pickling, the results as shown in Figure 2. The surface material of the coating layer is much looser than that of the substrate. It can be seen that the zinc coating is a thin film covering the surface of the substrate, and its related to the galvanizing process. This is also one of the main reasons that the coating is easily powdered and damaged when it is deformed and then affects the mechanical properties of the steel. Therefore, the mechanical properties of galvanized steel sheets and non-galvanized steel sheets should be different.  Samples are taken respectively on R00 degrees, R45 degrees, and R90 degrees from rolling direction, and then uniaxial tensile tests are executed. The Z050 electronic universal tensile testing machine was used to carry out the uniaxial tensile tests on the samples. The maximum experimental force is 50 kN, and it is set up with vertical and horizontal extensometer to measure the anisotropic parameters. It offers detailed mechanical parameters in Table 1. The results show that zinc coating involves a definite effect on the mechanical properties of steel sheet. Although the results have not been able to make a conclusion of consistency enhancement or reduction. Therefore, the zinc coating should not be neglected in a galvanized steel sheet, and it should establish its independent mechanical model to reveal accurately the mechanical characteristics of galvanized steel sheet.

15.1μm
Nano-indentation experiment can obtain the young's modulus, hardness and load displacement curve of zinc coating. Based on the dimensional analysis method, the material constitutive model of the zinc coating is constructed by processing these data. The detailed measurement principle of nanoindentation is shown in Figure 3. The hmax is the maximum depth of the indentation, the hf is the residual depth, and the Pmax is the maximum load. Contact depth hc refers to the depth of the contact between the pressure head and the pressed object when pressed into the measured material, as shown in Figure 3b. Depending on the power function law proposed by Oliver and Pharr [31], the loading and unloading curves are analyzed and fitted.
The prototype of the in-situ nano-mechanical tester used in this experiment is TI750, which is produced by Hysitron ® company in the Irvine city, California, USA. The maximum load of the instrument is 30mN, and the minimum load is 70 nN. It makes the indenter of Berkovich diamond, whose elastic modulus is Ei = 1141Gpa and Poisson's ratio is vi = 0.07. The test temperature is 25 °C, the loading mode is load control, and the constant strain rate is 0.05 s −1 . The results of loading and unloading curves are shown in Figure 4. It is assumed that the stress-strain curves of zinc coating are power hardening function models, as shown in Equation (1).  where E is the elastic modulus of the zinc coating which has been obtained by the above analysis. Only yield strength σy and hardening index n need to be solved.
In this paper, the two unidentified parameters are solved by the combination of dimensionless analysis and finite element numerical simulation. The related parameter dimensions are shown in Table 2 as follows. According to the above analysis and test results, the elastic modulus of the zinc coating is 135 GPa. The mathematical model of stress-strain is shown in Equation (2) and the curve is shown in in Figure 5a. The finite element simulation of Nano-indentation is carried out by using the model, and it shows the compared results with the test results in Figure 5b. The deviation of the two curves is within 2%, which shows that the model can reflect the mechanical properties of the coating. The finished constitutive model of the zinc coating is given in Equation (3). The yield stress is 647.72 Mpa and the yield strain is 4.79814 × 10 −3 .

Interface Binding Energy between Zinc Coating and Substrate
The interface binding force between the coating and the substrate refers to the binding strength between the coating and the substrate. The external force required by the coating per unit area from the substrate which is an important index for testing binding strength between the coating and substrate. The sliding module of the UMT-TriboLab ® friction test machine is used in this experiment. The spherical point pressure head with a cone angle of 120° and a radius of 200µm is used. After the surface is polished, DP600 galvanized steel sheet is preloaded with the force of 0.6 kg for a continuous 10s, the linear loading force of the maximum positive load is 2.0 kg, and the duration is 120s. It shows its load condition in Figure 6a. The morphology of the scratch is measured, and the results are displayed in Figure 6b. From the analysis of the results, the transverse displacement of the scratch initiation point and the interface failure point are (−11.12501, −12.663326), the length is about 1.538 mm, and the longitudinal depth is 15.1 µm. Figure 7 is the change of tangential load Fx, friction coefficient µ and acoustic signal AE with time t. The analysis shows that when the interface is broken, the tangential force will decrease sharply or to 0, and then gradually increase, which leads to the clear fluctuation of the friction coefficient and sound signal at the fracture point (the breakpoint in the Figure 7). As the polishing surface of the sample cannot be smooth, there may be impurities in the coating layer, and it does not treat the surface of the substrate which represents the coating and the substrate a staggered state of convex and concave. Therefore, the tangential force, friction coefficient and acoustic signal curve are not linear changes, but small fluctuations. However, there will be large changes in the interface failure. Considering this, we can judge the point of the critical load. Combined with the breakpoints of the curves, interface failure time is 92.22 s, and the corresponding critical load is 15.76 N, and the critical friction coefficient is 0.202. The failure of the coating results from the joint action of three stress states, which are the static pressure stress formed by the loading load, the friction stress of the contact area when the pressure head slips on the surface of the specimen, and the residual stress in the coating. It is well known that the pressure head is spherical and its radius Ri=200µm. The contact radius of the indenter and coating is obtained in Equation (4) [32].
where a is the contact radius between the indenter and the coating, L is the critical load, and V and E are Poisson's ratio and elastic modulus (s, i mark substrate and indenter). When the indenter is rigid, the Ei means infinity, and the Equation (4) can be reformulated as Equation (5)   (1 2 ) 2π where σxx, σyy, and σzz are three normal stress components, while σxy, σyz, and σzx are three shear stress components.
According to the above stress analysis, because the coating is very thin, it can be assumed that the strain of the coating will not change with its thickness. The strain of the coating and the substrate is the same at the interface. The coating strain in front of the indenter can be described as shown in Equation (7) [32].
where εxx, εyy, εzz are three normal strain components and µ is the coefficient of friction.
The strain of Z is expected to be very small and can be ignored. Using the above strain field, the elastic deformation energy stored in each element in the coating can be evaluated and the residual stress in the coating is considered. The detailed equations are described as follows [32]: where Uelastic is the elastic deformation energy per unit volume, ε'xx and ε'yy are residual strains, σres is the residual stress of the zinc coating, and Vf and Ef are the Poisson's ratio and elastic modulus of coating.
The thin coating in front of the indenter releases its elastic deformation energy through interfacial delamination and coating cracking, which leads to the failure of the coating. It is assumed that the coating failure area under the critical load is a semi-annular zone separated from the substrate and the radius is c, then the expression of the total energy stored in the coating is shown as Equation (11).
where U is the total energy stored in the coating, t is the thickness of the coating, c is the radius of the failed part of the coating, W is the interface binding strength of the interface, t is the plastic dissipation at the crack tip in the separation process, and γf is the surface energy of the coating. The failure behavior of the thin coating will continue when the energy released by the failure behaviors such as delamination and spalling of the coating is larger than that of the interface and surface energy. Therefore, the following Equation (12) [33] can be consulted from the Griffith energy balance condition. It will generate the interface binding strength of the interface, as shown in Equation (13) [33]. Depending on the experimental results and the above calculation, and neglecting the residual stress of the coating, the interface binding strength between the coating and the substrate is obtained, as shown in Table 3.
where W is the interface binding energy and γf/c can be obtained by estimating the elastic deformation energy stored in each unit volume of the coating.

Hemispherical Part Drawing Experiment
The hemispherical part drawing model is used for finite element analysis, as shown in Figure 8. The diameter and the thickness of the steel sheet is 100 mm and 1.2 mm, the outer diameter of the die and the blank holder are all 55 mm, the diameter of the punch is 20 mm, the inner diameter of the die and the blank holder are 20 and 27 mm respectively, and the profile radius and outer profile radius are all 0.75mm. In this paper, we use the BUT400 ® machine as the drawing experimental equipment. As shown in Figure 9, the maximum clamping force of the sample is 400 kN, the maximum stroke of the piston is 120 mm, and the working speed is 0.750 mm/min. The test sample is a hot galvanized steel sheet DP600 with a diameter of 100mm and a thickness of 1.2 mm. By monitoring the blank holder force, we observed the influence of drawing depth on the formability of galvanized sheet. Figure 10 shows four samples of the hemispherical part drawing test. We intend eight groups of hemispherical part drawing experiments to test the drawing depth and forming force of different blank holder forces. Detailed results of the test are shown in Table 4. The experimental results show that the drawing depth and the forming force are distinct at different blank holder force, and the drawing depth increases first and then declines with the increase of the blank holder force.  Figure 11 is a tendency analysis of the results of the hemispherical part drawing experiment. With the increase of the blank holder force, the drawing depth increased slightly. However, when the blank holder force increases to a certain extent, the drawing depth tends to be stable and remains unchanged. At the beginning, the forming force increases with the increase of the blank holder force, then decreases with the increase of the blank holder force and stabilizes. The experimental results are correct and consistent with expectations. Spray paint for speckle measurement was sprayed on the surface of the sample, and strain was calculated by GOM Argus ® optical strain analysis system. Figure 12 indicates the strain measurement

Damage Model Based on Cohesion Element
To explore and improve the zinc coating damage numerical simulation of galvanized steel sheet stamping, numerical simulation analysis was performed by using the blank holder force condition and the drawing depth as the stamping stroke. The zinc damage simulation of galvanized steel sheet is analyzed by establishing a finite element model considering the zinc coating. The model of finite element analysis is completely consistent with the experimental model. Abaqus ® 6.14 version software is used to carry out the finite element numerical simulation of hemispherical part drawing. To promote the efficiency of computation, the 1/4 model is adopted to analyze, and it adds symmetric constraints to the nodes on the X axis and the Y axis, as shown in Figure 13a. The final finite element analysis model is presented in Figure 13c. Three models are used to describe the material structure of galvanized steel sheet, as illustrated in Figure 13b. The element type of sheet metal is a solid element. The first model (Finite element model I, single-layer model) is to analyze the coating and substrate as a unified material, and the mechanical parameters are the data of uniaxial tension of galvanized steel sheet. The second (Finite element model II, three-layers model) and third (Finite element model III, five-layers model) models consider the substrate and the coating as two materials. Among them, the mechanical parameters of the substrate are the uniaxial tensile data for removing the zinc coating, but the mechanical parameters of the zinc coating are obtained by the Nano-indentation test and calculation. The adhesion between the coatings and the substrate of the second model are constrained by tie binding. Stiffness data transfer is equivalent to the rigid connection of the two binding surfaces and the binding area serves the absence of relative movement and deformation while the stiffness is large. However, in the third model, a further cohesive element is brought in between the coating and the substrate which is a zero-thickness element layer applied to simulate the generation and expansion of the crack in the binding layer. In the cohesive element model, the cracking of the coating is determined based on the maximum principal stress criterion, and the damage evolution process of the crack is conducted based on the interface binding energy.
Equation (12) shows the constitutive equation for the finite analysis of the zinc coating, and it shows the interfacial binding parameters between the coating and the substrate in Table 3. Three kinds of finite element model are used to simulate the hemispherical part drawing process of nine groups of different blank holder force and drawing depth. The detailed setting parameters are shown in Table 4. Stroke in numerical simulation is controlled by punch stroke, that is, the drawing depth. The total number of numerical simulations for the hemispherical part drawing is 27 times.

Simulation Results and Discussions
The numerical simulation results are analyzed and compared with the experimental results. Figure 14 is the numerical simulation result of hemispherical part drawing of the single-layer model and the three-layers model (binding) when blank holder force is 50 kN and punch stroke is 9.76 mm. The results show that the maximum Mises stress of the substrate is 7 MPa larger than that of the single-layer model. Mises stress distribution in the upper and lower surfaces is basically the same, with only a slight difference. However, their maximum Mises stress is nearly 100 MPa smaller than that of the substrate in the three-layers model.  Figure 15 shows results from the numerical simulation of the five-layers model under blank holder force 50 kN and punch stroke 9.76 mm, which has two 0-thickness cohesive elements to describe the linkage between the coating and the substrate. Figure 15a is the stress distribution simulation result of substrate, which is similar to the simulation result of the three-layers model. Comparing Figure 14c and Figure 15b with Figure 14d and Figure 15d, it is shown that although the maximum stress of the coating layers on the upper and lower surfaces is the same as that of the threelayers model, the stress distribution is quite different. Figure 15c,e are the simulated results of the degradation of upper and lower 0-thickness cohesive layers. The damage coefficient scalar stiffness degradation (SDEG) in almost all deformation zones is 1.0, indicating that the cohesive element has completely degenerated under the action of biaxial tensile stress. It can be concluded that the galvanized coating layers and the substrate have been disconnected. Figure 16a shows the simulation effect of the separation of the zinc coating from the steel substrate. Under the action of tension, the interface cracks gradually, the damage coefficient "SDEG" is different, and interface cracks are different, especially galvanized layer cracks most seriously at the corner of the die and the top of the punch. As shown in Figure 16b, when the simulation time is between (0.05, 0.06), the cohesion element begins to damage, and gradually cracks with the increase of stress. When the simulation time is about 0.08, the cohesion element layer is almost completely damaged, and obvious cracking occurs between the zinc coating and the steel substrate. Therefore, the substrate and the coating layers are not completely binding in the drawing process. Once the galvanized layer is detached from the substrate, owing to its tiny thickness scale, it will be powdered under the dual action of friction and stress. The results of pulverization will further affect the formability and surface quality of sheet metal by friction. Therefore, the influence of the galvanizing layer on forming performance is a dynamic process. If it is necessary to simulate the formability more accurately, it is necessary to change the friction coefficient dynamically according to the degeneration of the cohesive element. According to the above analysis results, the stress and strain between different models cannot be further compared because of different layers of the finite element models. Therefore, this paper focuses on comparing the accuracy of numerical simulation of different models by forming force. Figure 17 shows the comparison between the forming force of different numerical simulation models and experimental under four groups of different blank holder force and stroke. In the first half of the drawing, the numerical simulation results are larger than the experimental results, and the deviation is larger. However, with the drawing process of the latter half stroke, the numerical simulation value is smaller than that of experimental value, and the deviation between the two is decreasing. Moreover, 4 sets of results in Figure 16 show the same trend. However, on the whole, the deviation between the numerical simulation of a five-layer model and the experimental results is the smallest. Therefore, the results show that the forming force of the five-layer model is closer to the experimental results, and the numerical simulation accuracy is the highest. Comparatively, the accuracy of the single-layer model is the worst, while that of the three-layers model is in the middle. It can be observed on the above evaluation that although the forming force of the five-layers model is the closest to the experimental results, there is still some deviation from the experimental value. Especially in the beginning of the 2/3 stamping stroke, the difference is more significant. Therefore, the accuracy of numerical simulation needs to be further improved. Depending on the above analysis, friction may be a meaningful factor affecting the forming force. Especially under the condition of high strength and pressure, the powdering effect of the coating will be more obvious. However, the effect of the powdering of the zinc coating on friction has not been systematically studied. The next step of the paper is commensurate to study the dynamic friction in the stamping process of advanced high strength galvanized steel sheet. In order to further analyze the accuracy of numerical simulation of different finite element models, the maximum forming force was extracted and compared with the experimental results, as shown in Figure 18. Table 5 is the accuracy analysis of the numerical simulation results of the maximum forming force. The average deviation of the maximum forming force of single-layer model, three-layers model and five-layer model is 2.979%, 1.173%, and 0.645%, respectively. The results of the maximum forming force are consistent with the previous analysis, and the five-layers model is the closest to the experimental results and the highest accuracy. Moreover, the variation of the maximum forming force under the different blank holder force of the three models is consistent with the experimental results. This demonstrates that the three finite element models established in this paper can simulate the stamping process of hemispherical part drawing successfully. The material parameters obtained in the experiment and the constitutive model constructed in this paper can accurately describe the mechanical deformation behavior of galvanized steel sheet, which proves the feasibility and accuracy of the above research. In conclusion, interfacial binding energy is one of the important micro-mechanical properties of the zinc coating, which determines the forming quality, damage and service life of the coating steel sheet. In this paper, interfacial binding energy between zinc coating and substrate was achieved by the scratch test and related models, and applied to the numerical simulation of a hemispherical part drawing experiment to predict the damage evolution process of the coating. However, the experimental samples are directly from factories, and the surface morphology of the sample substrate changes greatly. The interface between the coating and the substrate is a state of staggering convexity and concavity, and the thickness distribution of the coating is nonuniform, which leads to the inaccuracy of the scratch test results. Therefore, the interface binding energy of galvanized sheet still needs to be further studied in the later stage of work.
In addition, because the fracture energy of the galvanized coating is difficult to obtain, and the sample size is large, the preset crack location is difficult to set, so the XFEM method is not used to investigate the powdering of the galvanized coating. Therefore, in order to further improve the accuracy of the finite element numerical model, it is necessary to conduct in-depth research based on the cohesion element model.

Conclusions
In this paper, the hot dip galvanized high strength steel DP600 is chosen as the research object. We studied mechanical behavior and its damage with the steel substrate. The experimental and numerical results show that the constitutive model can accurately describe the plastic deformation behavior of the coating, and the finite element model of the hemispherical part drawing has a high accuracy of damage numerical simulation. The main conclusions of this study are as follows: • The effect of the zinc coating on the formability of galvanized steel cannot be ignored. In this paper, we studied the zinc coating and its mechanical properties. The constitutive model of zinc coating based on Nano-indentation experiment and numerical model was proposed, and the constitutive model can accurately describe the plastic mechanical behavior of zinc coating.

•
The interface binding energy is an important index to reflect the zinc coating quality, and also a critical criterion to evaluate the coating damage. Based on the nano-scratch test and Griffith's energy theory, the interface binding strength between the coating and the substrate was obtained, which can be invoked as a criterion for the detachment between the galvanized layer and steel substrate. This provides a theoretical basis for the damage of the zinc coating in the stamping process.

•
The damage evolution of the binding between the zinc coating and the substrate was simulated by introducing a cohesive element with 0 thickness. It can simulate the detached and damage evolution of zinc coating from the substrate accurately. Compared with other multi-layer finite element models, it has improved the accuracy of the numerical simulation of stamping through the comparative analysis of the forming force history and the maximum forming force with the experimental results.
The constitutive model, interface binding energy and damage finite element model achieved in this paper can improve the numerical simulation accuracy of galvanized steel sheet stamping to a certain extent. But there is still some deviation from the experimental results. According to the analysis, this may be caused by the great change of interfacial friction after the zinc coating is damaged and detached from the steel substrate. Therefore, further studies will consider the dynamic friction model after coating damaging, and further improve the numerical simulation accuracy of advanced high strength galvanized steel sheet.