Abstract
In the forming process of galvanized sheet, the friction between the die and the blank often causes the zinc coating of galvanized sheet to peel off, scratch, and crack. The aim of this study is to evaluate and calculate the fractal characteristics of the surface morphology of galvanized sheet and the effect of pressure on the interfacial friction behavior. Two steel plates, GA and GI, produced by Shanghai Baosteel Company, were used as materials to conduct tribological experiments, measure the surface profile and threedimensional shape of the galvanized sheet, and calculate the fractal dimension and fractal roughness parameters. According to the analysis results of friction surface damage of galvanized sheet, the damage failure parameters of galvanized sheet are calculated. On this basis, according to the adhesive friction theory, the total surface friction value of galvanized sheet is obtained, and the fractal calculation model of galvanized sheet friction is established. The simulation results show that the galvanized sheet has fractal characteristics. The average values of fractal dimension and scale factor of SP781BQ alloy hotdip galvanized sheet are 1.52 and 0.23 µm, respectively. The average fractal dimension and scale coefficient of HC420/780DPD + Z hotdip galvanized sheet are 1.60 and 0.11 µm, respectively. The friction coefficient calculated by the proposed method is consistent with the theoretical value, and the error is less than 10%, which proves the accuracy and feasibility of the friction fractal calculation method.
1 Introduction
The surface of any friction pair in construction machinery is rough and uneven on the micro scale. The basic reason of friction and wear is that the two rough surfaces contact each other under the action of load, resulting in heat, deformation, adhesion, groove, and other effects. Therefore, the study of the contact behavior of micro rough surface has been one of the important topics in tribology [1 2 3].
In the physical field, in the literature [4], two different methods (GA and GI) were selected to evaluate the influence of surface morphology, pressure, and lubrication conditions on the friction behavior of galvanized steel plate interface. The surface morphology of the sample is measured by profilometer, fractal parameters are calculated by fractal theory, and the relationship between the actual contact area and the load is obtained. The contact model of galvanized plate and die is established; According to the contact model, the friction model of the relationship between the surface morphology and the load of the galvanized plate is established. Literature [5] introduced the deep learning representation and modeling method, characterized and modeled the galvanized sheet surface and fractal surface, and studied the modeling process of friction surface morphology. The results showed that the surface tribological parameters could not represent the surface treatment method, and the surface tribological parameters could not be reconstructed by the specified processing methods or parameters. Moreover, fractal surface modeling cannot achieve both fractal parameters and surface tribological parameters. Therefore, the subsequent simulation process cannot meet the requirements of scale and parameters. The depth learning representation and modeling method are introduced to characterize and model the surface of galvanized sheet and fractal surface. In the field of engineering, the main research contents of literature [6] include optimizing fractal surface modeling method and improving SSK distribution of fractal surface morphology. The multiscale modeling and simulation of the surface topography of turning milling are carried out using the multidirectional continuous wavelet, and the compression sensing simulation is carried out to reduce the sampling rate. The superresolution reconstruction process of rough surface topography scanning is complex. A method of fractal dimension calculation of 3D rough surface based on wavelet transform is proposed in document [7]. The method decomposes the 3D surface topography data by wavelet transform. The threedimensional fractal dimension of rough surface is obtained by fitting the mean value of the square of wavelet coefficients under different decomposition scales. To verify the correctness and accuracy of the method, the 3D fractal surface is constructed under the condition of known parameters, and the results of fractal dimension calculation under different methods are compared and analyzed. The results provide the parameter basis for the establishment of fractal contact model of rough surface. However, in practical application, the friction coefficient calculated by the above traditional method is different from the theoretical value, and the fractal dimension and proportion coefficient of galvanized sheet are unreasonable.
To solve the defects of traditional methods, this paper studies the genetic algorithm and genetic algorithm of galvanized sheet. Based on the fractal contact theory, the surface morphology of the galvanized plate is measured by optical profilometer, the damage and failure parameters of the galvanized plate are calculated, and the fractal calculation model of friction of the galvanized plate is established. The results are compared with the experimental data. The innovation is that according to the damage analysis results of the friction surface of the galvanized plate, the damage failure parameters of the galvanized plate are calculated. According to the adhesion friction theory, the total surface friction value of the galvanized plate is obtained, and the fractal calculation model of the friction of the galvanized plate is established.
2 Surface morphology and fractal characteristics of galvanized sheet
Two kinds of steel plates, GA and GI, produced by Shanghai Baosteel Co., Ltd., were used as experimental materials to study the tribological properties of the two kinds of steel plates, SP781BQ and HC420/780DPD + Z, respectively [8,9,10]. The two kinds of galvanized steel sheets were cut into rectangular samples with the size of 30 mm × 100 mm. The surface of all samples shall be cleaned with absolute ethanol, kept clean, and then dried with cold air. The surface morphology of galvanized sheet was measured by Veeco Wyko nt1100 noncontact optical profilometer [11]. Three micro regions on the surface of the sample were randomly selected as the measurement surface to obtain the surface profile curve, and determine the surface roughness value and threedimensional surface morphology. The sliding friction experiment was carried out. The experiment only studies the friction of the galvanized sheet under different mold roughness, different pressure, and no lubrication conditions. The surface roughness of the die is 0.4, 0.8, and 1.2 µm, respectively. Under the positive pressure of 1,000 N, the contact surface between the upper and lower die and the sample is about 2 mm × 30 mm rectangular area [12,13,14]. The friction sliding speed is set at 150 mm min^{−1}. To reduce the experimental error, the average value of all experimental data was taken after three measurements.
Three samples of GA and GI galvanized sheet were prepared. Randomly select two microzones on the sample surface as the measurement area (Figure 1), and use Myko NT1000 profiler to measure the surface profile and roughness of the sample (Figure 2); a better area was selected on the galvanized steel sample, and its threedimensional morphology was measured by laser confocal microscope.
Figure 1
Figure 2
Figures 3 and 4 show the measured surface profiles of GA and GI plates (position 1 on the left and position 2 on the right).
Figure 3
Figure 4
It can be seen from the comparison of Figures 3 and 4 that the surface of the GA plate of the alloyed hotdip galvanized steel sheet is relatively rough, and the degree of unevenness is relatively large. The surface of the pure zinc galvanized steel GI plate is relatively flat, and the degree of unevenness is small. The average surface roughness measured is 0.889 and 0.525 µm, respectively. This is because the coating of GA plate is composed of hard brittle phase Zn Fe alloy, while the coating of GI plate is composed of soft Zn. The roughness of hard GA plate is larger than that of soft GI plate.
The fractal dimension and fractal roughness parameters calculated based on Figures 3 and 4 are shown in Table 1.
Brand  Fractal dimension

Scale coefficient

Average fractal dimension  Average scale coefficient (


GA (ST781BQ)  1.522  0.227  1.52  0.23 
1.528  0.229  
1.497  0.232  
GA (ST781BQ)  1.597  0.114  1.6  0.11 
1.604  0.111  
1.601  0.112 
It can be seen from Table 1 that the fractal dimensions of the two kinds of galvanized sheets meet the requirements of 1 < d < 2, so the two kinds of sheets have fractal characteristics [15]. The average fractal dimension of GA plate is smaller than that of GI plate, but the average value of scale coefficient is larger than that of GI plate. This is because the smooth surface has a deeper selfsimilar fine structure, and thus has a higher fractal dimension.
3 Relationship between surface morphology and load of galvanized sheet
3.1 Damage analysis of load friction surface of galvanized sheet under different load directions
The friction interface element of the galvanized sheet is only subjected to normal and two tangential load stresses, and its constitutive correlation [16] expression is recorded as:
Among them,
In a single form, in the bilinear constitutive model of the interface element on the surface of galvanized sheet, the analytical expression of the corresponding interface displacement
The interface relative displacement of the global failure point control [17] is as follows:
In general, the damage initiation and broadening degree of interface elements are generated under mixed loading. It is necessary to consider the effect of coupling interface damage between the stress components of three different cracking states. If the shear strength of the two interfaces is equal, the energy release rate of crack strain and the interface stiffness are the same.
The analytical expression of the control displacement in the mixed mode is recorded as follows:
Assuming the parameter
The final failure uses the principle of secondary strain energy release rate to obtain the global damage displacement in the mixed state:
For a representative interface, its area A can be regarded as the sum of undamaged area n and damaged area m, which is specifically expressed as [18]:
The stress in the representative interface is:
Among them,
In the case of neglecting the constitutive model of friction interface element, it is considered that the interface element will not generate normal interface damage under normal pressure. Therefore, when the interface is damaged, only normal compressive stress can exist in the interface area where debonding damage has occurred. However, when the interface unit is debonded, the debonded surface becomes rougher. At this time, when the interface is subjected to normal compression load, the corresponding slip of the interface with respect to the other phases will generate strong friction force in the debonding interface [19,20,21].
Considering the external friction of galvanized sheet at debonding interface, the representative interface internal stress in equation (9) can be transformed into:
Among them,
3.2 Fractal calculation of friction damage failure parameters
According to the fractal calculation results of friction damage on the surface of galvanized sheet, the improved Hashin’s friction damage failure principle is used as the judgment standard of fiber bundle failure, because the material strength evaluation value obtained using Hashin friction damage failure principle is far less than the experimental value. By adding the weight element of shear stress component [22,23], the bad effect of shear stress component on the failure of fiber element is reduced. The judgment basis of the improved threedimensional Hashin friction damage failure principle is as follows.
The expression of fiber cracking failure in galvanized sheet fiber bundle is as follows:
The analytical formula of resin fracture failure in galvanized sheet fiber bundle is as follows:
Among them,
The mechanical properties of different volumes of galvanized sheet were tested. The residual stiffness and strength models of longitudinal and surface shear friction damage of unidirectional galvanized sheet with different fiber volume components were obtained by fitting.
The galvanized sheet unidirectional area longitudinal tensile friction damage test can obtain the residual stiffness and strength models related to the number of friction damage
The friction damage times
The surface of galvanized sheet is pushed by friction damage load. With the increase in the number of friction damage loading cycles
With the increase in the number of friction damage, the properties of unit materials in the cell gradually decline. After a certain degree of friction damage, the element will generate damage failure, and the material properties of the element will degrade suddenly after failure. To accurately judge the friction damage of materials, it is necessary not only to evaluate the failure mode principle of single cell, but also to define the basis of final damage on the surface of galvanized sheet. Under the action of friction damage load, when the damage is widened to a fixed level, the surface of galvanized sheet will not have the bearing capacity, and the global structure will eventually be damaged.
When simulating progressive damage under static load, when the stress–strain curve has an inflection point, the average unit cell strain will reach the cracking strain. If the damaged unit has covered the surface of all fiber bundles, the single cell structure cannot be sustained. The static cracking strain is chosen as the basis to judge the failure of single cell structure. That is to say, under the action of friction damage load, when the average cell strain reaches the static cracking strain, the global failure of the surface structure of galvanized sheet is judged. In the friction damage calculation program, only the influence of the maximum friction damage stress load on the single friction damage performance is considered. The residual strength model and stiffness model also complete degradation under the highest friction damage stress load.
3.3 Fractal calculation model for friction of galvanized sheet
According to the adhesive friction theory, the total friction can be expressed as follows:
When the shear strength of the zinc layer of the alloyed hotdip galvanized sheet GA is 113 MPa, the friction coefficient is
The relationship between the friction coefficient and the positive pressure of two different galvanized sheets is shown in Figure 5. It can be seen from the figure that the friction coefficient of the two galvanized sheets decreases with the increase in load. The friction coefficient of GA plate decreases greatly in the initial stage, but decreases slowly when the load reaches 3,000 N. In contrast, the friction coefficient of GI plate is relatively gentle in the whole process of decline.
Figure 5
4 Simulation experiment design and result analysis
4.1 The accuracy of friction damage judgment of galvanized sheet surface under the proposed method
Simulation experiments are designed to verify the accuracy of the fractal calculation method of friction damage. According to the static test of the specimen, the static strength used in this paper is 480 MPa. According to the load movement curve of the static test, the cracking strain can be calculated to be 1.68%.
Figure 6 shows the S–N (friction stress) curve of friction damage calculation.
Figure 6
It can be seen from Figure 6 that the friction damage and stress level curves of the calculation method in this paper have good linearity, and the S–N curve close to the surface of the galvanized sheet under various stress tracking results is linear. In the simulation experiment, it can be seen that the shear load has more influence on the friction damage. In this paper, the improved threedimensional Hashin principle is used to reduce the influence of shear model on the evaluation of friction damage, so as to enhance the accuracy of friction damage calculation.
Figure 7 shows the transformation of the ratio of friction damage elements of galvanized sheet with normalized damage under different stress levels.
Figure 7
It can be seen from Figure 7 that under the friction damage stress level, the element expansion rate of galvanized sheet is faster, and the element damage ratio increases exponentially. The reason is that when a certain element has damage failure, after the stiffness matrix is reduced by the principle of sudden drop of stiffness, because of the phenomenon of stress accumulation. As a result, the adjacent elements will grow rapidly under the next cycle load, and then the expansion rate will increase significantly. Under different stress levels, it is reliable to judge the structural failure using the failure damage ratio of single cell structure, and it also shows that the crack strain used in this method as the evaluation standard of structural failure is more reasonable. At the same time, it can be seen that the higher the friction damage stress level, the faster the fiber bundle from the initial damage to the final failure of the structure, and the faster the failure rate.
4.2 The accuracy test of the friction fractal calculation result of the proposed method
To carry out the sliding friction experiment of galvanized sheet, the friction coefficient of galvanized sheet under different loads was measured. SP781BQ galvanized sheet and HC420/780DPD + Z galvanized sheet were used in this experiment. The diversity of galvanized sheet materials can enhance the reliability of the experimental results. The friction coefficient is calculated using the friction model established above. The comparison between the theoretical value and the experimental value is shown in Tables 2 and 3.
Category  Load (N)  

1,000  1,500  2,000  
Theoretical value  0.1853  0.1625  0.1451 
Experimental value  0.2058  0.1751  0.1554 
Error (%)  7.5  9.1  7.5 
Category  Load (N)  

1,000  1,500  2,000  
Theoretical value  0.1350  0.1388  0.1276 
Experimental value  0.1501  0.1459  0.1365 
Error (%)  10  7.4  5.7 
It can be seen from Tables 2 and 3 that the theoretical value of friction coefficient of the two kinds of galvanized sheet is close to the experimental value, and the error between them is less than 10%. The theoretical value of the friction coefficient is slightly less than the experimental value because the influence of ploughing force is ignored.
5 Conclusion
 (1)
The galvanized sheet has fractal characteristics. The average fractal dimension and scaling factor of SP781BQ hotdip galvanized sheet were 1.52 and 0.23 µm respectively, and the average fractal dimension and scaling coefficient of HC420/780DPD Z hotdip galvanized sheet were 1.60 and 0.11 µm respectively.
 (2)
In this paper, a friction model is established and the stressfriction curve and the loaddamage unit curve obtained are consistent with the general law of damage expansion. Therefore, the frictional damage fractal calculation method proposed in this paper can accurately and effectively judge the friction damage of galvanized sheet surface, and has high applicability and robustness.
 (3)
The test results are verified by comparing the sliding friction of the plate. The results show that the theoretical values of friction coefficients are in good agreement with the experimental values. The average error is less than 10%, which proves the accuracy and feasibility of the frictional fractal model.
Acknowledgment
The study was supported by “Youth Teacher Growth Program of Hunan University (2018).”
Conflict of interest: Authors state no conflict of interest.
Data availability statement: All data generated or analysed during this study are included in this article.
References
[1] Zhao Y, Luo Y, Zhang Z, Zhang H, Guo X, Wang S, et al. Fractal dimension characterization of joint surface morphology on dissimilar friction stir lap welding of Al/Mg. Materials. 2019;12(23):3941. Search in Google Scholar
[2] Goritskii VM, Kulemin AM, Silina NG. Structure and mechanical properties of colddeformed galvanized rolled sheets for light thinwalled constructions. Russian Metall (Metally). 2019;2019(4):431–5. Search in Google Scholar
[3] Papanikolaou M, Salonitis K. Fractal roughness effects on nanoscale grinding. Appl Surf Sci. 2019;467–8(15):309–19. Search in Google Scholar
[4] Ivanovici M. Fractal dimension of color fractal images with correlated color components. IEEE Trans Image Process. 2020;29(12):8069–82. Search in Google Scholar
[5] Sun G, Zhu H, Ding C, Zhou Y. Multifractal detrended fluctuation analysis on friction coefficient during the friction process. J Tribol. 2018;140(1):011601. Search in Google Scholar
[6] Zhao C, Li J. Equilibrium selection under the Bayesbased strategy updating rules. Symmetry. 2020;12(5):13–6. Search in Google Scholar
[7] Cui H, Wang Q, Lian Z, Li L. Theoretical model and experimental research on friction and torque characteristics of hydroviscous drive in mixed friction stage. Chin J Mech Eng. 2019;32(5):69–79. Search in Google Scholar
[8] Cao W, Wang WD, Xu HS, Sergeyev IV, Struppe J, Wang X, et al. Dynamic states recognition of friction noise in the wear process based on moving cut dataapproximate entropy. J Tribol. 2018;140(5):36–44. Search in Google Scholar
[9] Li G, Zhang K, Gong J, Jin X. Calculation method for fractal characteristics of machining topography surface based on wavelet transform. Procedia CIRP. 2019;79(7):500–4. Search in Google Scholar
[10] Jinasena A, Kaasa GO, Sharma R. Improved realtime estimation of return flow rate of drilling fluids by model adaptation for friction parameter. IEEE Sens J. 2019;19(20):9314–23. Search in Google Scholar
[11] Stütz M, Buzolin R, Pixner F, Poletti C, Enzinger N. Microstructure development of molybdenum during rotary friction welding. Mater Charact. 2019;151(9):506–18. Search in Google Scholar
[12] Simoni L, Beschi M, Legnani G, Visioli A. Modelling the temperature in joint friction of industrial manipulators. Robotica. 2019;37(5):906–27. Search in Google Scholar
[13] Wei HY, Hussain G, Iqbal A. Surface roughness as the function of friction indicator and an important parameterscombination having controlling influence on the roughness: recent results in incremental forming. Int J Adv Manuf Technol. 2019;101(9–12):2533–45. Search in Google Scholar
[14] Lin HL, Chen HY, Wang YJ. Effects of helium addition to Arbased shielding gas on the lapjoint performance of 5052 aluminum alloy to galvanized steel sheet during the MIG weldbrazing process. Solid State Phenom. 2020;305(1):129–35. Search in Google Scholar
[15] Liu Y, Hu B, Wu S, Wang M, Zhang Z, Cui B, et al. Hierarchical nanocomposite electrocatalyst of bimetallic zeolitic imidazolate framework and MoS2 sheets for nonPt methanol oxidation and water splitting. Appl Catal B Environ. 2019;258(4):117970. Search in Google Scholar
[16] AlSmadi M, Freihat A, Khalil H, Momani S, Ali Khan R. Numerical multistep approach for solving fractional partial differential equations. Int J Comput Methods. 2016;14(2017):1750029. Search in Google Scholar
[17] AlSmadi M, Arqub OA, Momani S. Numerical computations of coupled fractional resonant Schrödinger equations arising in quantum mechanics under conformable fractional derivative sense. Phys Scr. 2020;95(7):075218 (21pp). Search in Google Scholar
[18] AlSmadi M, Arqub OA. Computational algorithm for solving Fredholm timefractional partial integrodifferential equations of Dirichlet functions type with error estimates. Appl Math Comput. 2019;342(5):280–94. Search in Google Scholar
[19] Asa B, Oaac D, Dz E. Fuzzy fractional differential equations under the MittagLeffler kernel differential operator of the ABC approach: theorems and applications. Chaos Solitons Fractals. 2021;146:110891. Search in Google Scholar
[20] Cheng HH, Kung PT, Wang BR, Chiu LT, Tsai WC. A new regularization method for dynamic load identification. Sci Prog (1916). 2020;103(3):3685042093128–3. Search in Google Scholar
[21] Dagdeviren OE. Exploring load, velocity, and surface disorder dependence of friction with onedimensional and twodimensional models. Nanotechnology. 2018;29(31):315704–15. Search in Google Scholar
[22] Wu LJ, Li Z, Wang D, Yin H, Huang X, Zhu ZQ. Onload field prediction of surfacemounted PM machines considering nonlinearity based on hybrid field model. IEEE Trans Magn. 2019;55(3):1–11. Search in Google Scholar
[23] Wang L, Peng Y, Xie Y, Chen B, Du Y. A new iteration regularization method for dynamic load identification of stochastic structures. Mech Syst Signal Process. 2021;156(4):107586. Search in Google Scholar
© 2021 Xiaoyong Yang and Tan Jin, published by De Gruyter
This work is licensed under the Creative Commons Attribution 4.0 International License.