# Multifeature Extraction of Three-Dimensional Topography of Carbon Steel Specimen during Fatigue Process

Authors: Tao Liu; Jingxiong Wu; Jingfa Lei; Xue Wang; Bingqi Zhang; Shu Zhang
Journal: Advances in Materials Science and Engineering
 (2021)

---

## Abstract

In order to investigate the variation of three-dimensional metal surface topography during fatigue process, a three-dimensional (3D) topography acquisition platform was built with an in situ tensile tester and a three-dimensional profilometer. Q235 steel specimens were chosen as research objects, and the three-dimensional surface topography information at various stages of fatigue damage was obtained. Through the characterization of three-dimensional roughness, combined with surface height distribution and multifractal analysis, the variations of metal surface topography in the fatigue process were described. Results show that the arithmetic mean deviation of the surface (Sa), the width of the multifractal spectrum (Δα), and the mean value of surface height distribution (μ) and its standard deviation (δ) increase nonlinearly with the increase of fatigue cycles. The rate of fatigue damage is slow in the early stage and high in the middle and late stages. The surface height distribution amplitude (A) decreases with the increase of fatigue cycles, which indicates that the height data concentration decreases, and the metal surface becomes uneven. The Bayesian data fusion method was applied to establish a nonlinear mapping between the topography features and the damage, with the above five characteristic parameters (Sa, Δα, A, μ, and δ) as the data layer. Finally, a surface topography feature fusion method is proposed, and a case study is conducted to verify its applicability. The research results can provide reference for fatigue damage assessment.


---

## Body

## 1. Introduction

Metal surface topography is closely related to its fatigue strength and wear resistance, which affects the service life and reliability of components. When metal components, especially polycrystalline material components, are subjected to cyclic loading, the slip generated by the persistent slip band in the grains of the surface can be largely transferred, which makes the surface rough and causes corresponding changes in topographical characteristics [1]. The variation of surface topography is applied to describe fatigue damage, which is expected to bring new direction and ideas for fatigue damage assessment.At present, several methods have been proposed to describe the change of metal surface topography, such as two-dimensional (2D) gray-scale image processing of metal surfaces and feature extraction of energy and entropy parameters from the gray level co-occurrence matrix [2–4]. In addition, the laser speckle pattern is also applied to characterize the surface roughness [5]. In the research of surface topography evolution during the fatigue damage process, 2D and 3D roughness parameters are adopted to explore the inherent relationship between damage degree and topography [6–9]. It is found that, with the increase of fatigue damage degree, the roughness, frequency distribution, and fractal dimension of the surface profile are also increased [10]. The linear correlations between fatigue damage accumulation and roughness and between roughness and stiffness degradation are elucidated [11]. A composite damage index, which can reflect the accumulation of surface topography and topological changes in each fatigue damage state [12], is proposed, and a crack path prediction algorithm is established on this basis [9]. However, in the process of 3D surface topography extraction, even if a high-precision image acquisition instrument is utilized, the gray-scale image obtained often suffers from poor accuracy after 2D gray-scale processing, and the characteristic data is susceptible to environmental interference. The roughness parameters that characterize the topography, such as Ra (the arithmetic mean deviation of the profile), are scale-dependent quantities, and it is arduous to describe the self-similarity, hierarchical characteristics, and detailed changes of surface topography.Surface height information and its distribution provide the basis for extracting surface topography features. Meanwhile, fractal description has the advantage of scale invariance, which does not depend on the observation scale and resolution of the instrument, and can reflect the self-similarity of the rough surface structure at different scales [13, 14]. In this paper, Q235 steel specimens are selected to perform the axial-tensile fatigue test. The surface height information under various fatigue cycles is collected to obtain 3D roughness parameters such as the surface arithmetic mean deviation. The amplitude, mean value, and standard deviation characteristics of surface height distribution in the fatigue process are analyzed, and the multifractal spectrum characteristics from the early stage of fatigue to fracture are obtained. Based on the Bayesian data fusion method, a multifeature data fusion evaluation model of surface topography in the fatigue damage process is established.


## 2. Multifeature Analysis of 3D Surface Topography during the Fatigue Damage Process

### 2.1. Basic Characteristic of 3D Surface Topography

In the field of mechanical engineering, 3D surface topography refers to the microscopic characteristics of the surface of parts, such as peak-valley characteristics, texture direction, and shape deviation. These characteristics are often directly related to the surface properties of parts (fatigue strength, wear resistance, corrosion resistance, etc.). When the metal material/component is subjected to cyclic loading, the cyclic slip (extrusion and intrusion) leads to crack nucleation, and the irreversibility of shear displacement along the slip band makes the material surface rough. The rough surface is shown in the form of micropeak (extrusion) and valley (intrusion).The traditional measurement and evaluation of surface topography is mainly based on 2D and 3D surface roughness. The 3D roughness can reflect the microtopography and performance characteristics of the target surface more comprehensively and intuitively. ISO 25178 divides 3D surface roughness parameters into six types: height parameters, spacing parameters, function and related parameters, mixing parameters, characteristic parameters, and other parameters [15, 16]. The most commonly used is the arithmetic mean deviation Sa of height parameters. The schematic diagram of the surface arithmetic mean deviation Sa (unit: μm) is shown in Figure 1, which refers to the arithmetic average height in the sampling area.Figure 1
Schematic diagram of the surface arithmetic mean deviationSa.The mathematical expression is as follows:(1)Sa=1A∬Azx,ydxdy,where A is the sampling area and zx,y is the height of the point at (x, y).


### 2.2. Multifractal Description of Metal Surface Topography during the Fatigue Process

Multifractal, also known as fractal measure, is a singular combination of multiple subsets, and each subset has different fractal behavior. The fractal dimension distribution of multifractal is uneven, which is based on the extension of the original fractal concept to the nonuniform fractal, and its description of the fine structure of the fractal is hierarchical. The surface topography of the metal fatigue process is highly nonlinear and complex, and the use of multifractal can describe its local characteristics more effectively. The calculation process of multifractal parameters is as follows.The surface height matrix of the metal fatigue process is divided inton regions Si (i = 1, 2, 3, ..., n). Let εi be the scale size of the ith region and Pi be the generating probability of the region Si. For each region with scale size ε, the generating probability is as follows:(2)Piε=Siε∑Siε,where Siε represents the sum of the height values at all positions of the ith area and ∑Siε indicates the sum of heights at the whole area.In a scale-free self-similar area, the following relationship exists for probabilityPi(ε):(3)Pi∝εiαi,where α is the singularity scale index, which is a quantity reflecting the degree of singularity of each small area in fractal.The relationship between the number of probability subsetsNα(ε) (corresponding to α) and ε under the scale ε is defined as the following:(4)Nαε∝ε−fα,where Nαε is the number of probability subsets corresponding to α under the scale ε and fα is the fractal dimension with the same α subset, and the multifractal spectrum is composed of α and fα under different scales.Let the partition functionχq(ε) be the weighted sum of the probability Pi(ε) to the power of q:(5)χqε=∑Pqε,where q is the weighting factor.Using the unified parameter measurementu(q, ε) for calculation,(6)uq,ε=Pq∑Pqε.Ifq > 1, the high probability set plays the main role, and if q < 1, the low probability set plays the main role. The weighting function is to have a more comprehensive understanding of the fine structure on the target surface.The calculation formula of the multifractal spectrum is as follows:(7)αq=limε⟶0∑uq,εlnPεlnε,fq=limε⟶0∑uq,εlnuq,εlnε.The multifractal spectrum image can be obtained by using all theα and f(α) data. The common characteristic parameter of the multifractal spectrum α-f(α) is multifractal spectral width Δα = αmax-αmin. Δα characterizes the uneven distribution of objects with fractal characteristics, and it can reflect the fluctuation of the elevation signal. The larger the Δα value, the more severe the elevation fluctuation, and the smaller the Δα value, the smaller the elevation fluctuation.


### 2.3. Metal Surface Height Distribution and Numerical Characteristics during Fatigue

The random surface topography is determined by the height distribution with position coordinates as variables. It is difficult to describe surface topography using strict mathematical expressions. However, the statistical parameters of surface topography can be used for the evaluation of surface topography. Statistical frequency distribution refers to the interval division of each observation value according to a certain interval length so that each observation value forms a different interval distribution. Due to the difference in the nature of objective phenomena, their frequency distribution curves will show different shapes. The height distribution is a reflection of the surface profile shape, which plays an important role in the study of the size, shape, and distribution of the surface profile. Therefore, the surface height distribution can be adopted to statistically describe the changes in metal surface topography during fatigue.Normal distribution, also known as Gaussian distribution, is the most common probability distribution of continuous random variables. The common physical quantity that obeys normal distribution includes the measurement error of object size, the measurement of biological shape, and the quality index of products under normal production conditions. On the basis of normal distribution, combined with the specific parameter of surface height data, the statistical analysis of surface height distribution in the metal fatigue process is carried out. The details are as follows.If the distribution of surface height (z) of the fatigue process satisfies(8)fz=12πσe−z−μ2/2σ2,−∞<z<+∞,then the metal surface height in the fatigue process obeys the normal distribution, which is recorded as follows:(9)z∼Nμ,σ2.In equation (9), μ is the mathematical expectation (mean value) and the position parameter, which determines the location of the distribution, that is, the dominant numerical value of the surface height in the process of metal fatigue. σ (the standard deviation) is the scale parameter. When μ is a certain value, the shape of the curve is determined by σ, and the smaller the σ value is, the higher and thinner the curve is, which indicates that the metal surface height distribution is more concentrated, the uniformity is better, and the surface is smoother. The larger the σ value is, the more squat the curve is, which indicates that the metal surface height distribution is more dispersed, the uniformity is worse, and the surface is more rough and uneven. In addition, the amplitude A of surface height distribution is the frequency corresponding to the highest point of the normal distribution curve of surface height.


### 2.4. Surface Topography Feature Fusion Method

Naive Bayesian classifier is a simple and effective linear classifier. When the features are independent of each other, especially in the case of small-scale samples, the Bayesian model can show a good classification effect because of the small amount of sample data, and the feature variables are independent of each other, so the Bayesian method is applied for surface feature fusion.There are three kinds of Bayesian models: the polynomial model, Gaussian model, and Bernoulli model. The polynomial model and Bernoulli model are suitable for discrete data variables, while the Gaussian model is suitable for continuous data variables. In this paper, the Gaussian model is used as the fusion classification model of sample surface features.The Naive Bayesian method is a learning method using the Bayesian formula. The Bayesian method forms prior distribution by quantifying prior information and applies prior distribution to statistical inference to improve the accuracy of statistical inference. Naive means that it assumes the components of input variables (vectors) are independent of each other:(10)PXY=PX/Y∗PYPX,where P(Y) is the prior probability, PX/Y is conditional probability, obtained through historical data statistics, and PY/X is the posterior probability, which is the target.The Gaussian model is a kind of the Naive Bayesian model. It assumes that each featurexi obeys normal distribution under each category y. The probability density function of normal distribution is adopted to calculate the probability in the algorithm as follows:(11)Pxici=12πσ2exp−xi−μ22σ2,where μ represents the mean value of feature xi in the sample of category ci, i = 1, 2, 3, …:(12)μ=∑inxiN.σ2 stands for the variance of feature xi in the sample of category y:(13)σ2=∑inxi−μ2N,whereN is the number of training samples.It is known from the Naive Bayesian formula that, for the posterior probability values of different categoriesci, the denominator of the Bayesian formula is always the same, so only the numerator is required to be solved. Taking category c1 as an example, the posterior probability formula is as follows:(14)Pc1|xi=Pc1⋅∏i=1nPxi/c1i=1,2,3....In the same way, the posterior probabilities of other categories can be obtained.As described in Sections2.1–2.3, the five features describing the surface topography are arithmetic mean deviation Sa, surface height distribution amplitude A, mean value μ, standard deviation σ, and multifractal spectral width Δα, respectively. The above features correspond to xi in equation (11). The features are independent of each other and meet the requirements of the Naive Bayesian method. The feature data are all continuous variables, so the Gaussian model is selected as the fusion model. In the Gaussian model, the mean and variance of fatigue damage states (i.e., cyclic loading cycles, corresponding to ci of equation (11)) are calculated, and then, the probability density function of normal distribution is obtained. According to the probability density function, the probability density value of continuous test data is calculated to obtain the prior probability. The posterior probability value corresponding to each fatigue damage state is obtained by substituting the probability density value into the posterior probability formula equation (14). The fatigue damage state corresponding to the maximum posterior probability value is the prediction result.


## 2.1. Basic Characteristic of 3D Surface Topography

In the field of mechanical engineering, 3D surface topography refers to the microscopic characteristics of the surface of parts, such as peak-valley characteristics, texture direction, and shape deviation. These characteristics are often directly related to the surface properties of parts (fatigue strength, wear resistance, corrosion resistance, etc.). When the metal material/component is subjected to cyclic loading, the cyclic slip (extrusion and intrusion) leads to crack nucleation, and the irreversibility of shear displacement along the slip band makes the material surface rough. The rough surface is shown in the form of micropeak (extrusion) and valley (intrusion).The traditional measurement and evaluation of surface topography is mainly based on 2D and 3D surface roughness. The 3D roughness can reflect the microtopography and performance characteristics of the target surface more comprehensively and intuitively. ISO 25178 divides 3D surface roughness parameters into six types: height parameters, spacing parameters, function and related parameters, mixing parameters, characteristic parameters, and other parameters [15, 16]. The most commonly used is the arithmetic mean deviation Sa of height parameters. The schematic diagram of the surface arithmetic mean deviation Sa (unit: μm) is shown in Figure 1, which refers to the arithmetic average height in the sampling area.Figure 1
Schematic diagram of the surface arithmetic mean deviationSa.The mathematical expression is as follows:(1)Sa=1A∬Azx,ydxdy,where A is the sampling area and zx,y is the height of the point at (x, y).


## 2.2. Multifractal Description of Metal Surface Topography during the Fatigue Process

Multifractal, also known as fractal measure, is a singular combination of multiple subsets, and each subset has different fractal behavior. The fractal dimension distribution of multifractal is uneven, which is based on the extension of the original fractal concept to the nonuniform fractal, and its description of the fine structure of the fractal is hierarchical. The surface topography of the metal fatigue process is highly nonlinear and complex, and the use of multifractal can describe its local characteristics more effectively. The calculation process of multifractal parameters is as follows.The surface height matrix of the metal fatigue process is divided inton regions Si (i = 1, 2, 3, ..., n). Let εi be the scale size of the ith region and Pi be the generating probability of the region Si. For each region with scale size ε, the generating probability is as follows:(2)Piε=Siε∑Siε,where Siε represents the sum of the height values at all positions of the ith area and ∑Siε indicates the sum of heights at the whole area.In a scale-free self-similar area, the following relationship exists for probabilityPi(ε):(3)Pi∝εiαi,where α is the singularity scale index, which is a quantity reflecting the degree of singularity of each small area in fractal.The relationship between the number of probability subsetsNα(ε) (corresponding to α) and ε under the scale ε is defined as the following:(4)Nαε∝ε−fα,where Nαε is the number of probability subsets corresponding to α under the scale ε and fα is the fractal dimension with the same α subset, and the multifractal spectrum is composed of α and fα under different scales.Let the partition functionχq(ε) be the weighted sum of the probability Pi(ε) to the power of q:(5)χqε=∑Pqε,where q is the weighting factor.Using the unified parameter measurementu(q, ε) for calculation,(6)uq,ε=Pq∑Pqε.Ifq > 1, the high probability set plays the main role, and if q < 1, the low probability set plays the main role. The weighting function is to have a more comprehensive understanding of the fine structure on the target surface.The calculation formula of the multifractal spectrum is as follows:(7)αq=limε⟶0∑uq,εlnPεlnε,fq=limε⟶0∑uq,εlnuq,εlnε.The multifractal spectrum image can be obtained by using all theα and f(α) data. The common characteristic parameter of the multifractal spectrum α-f(α) is multifractal spectral width Δα = αmax-αmin. Δα characterizes the uneven distribution of objects with fractal characteristics, and it can reflect the fluctuation of the elevation signal. The larger the Δα value, the more severe the elevation fluctuation, and the smaller the Δα value, the smaller the elevation fluctuation.


## 2.3. Metal Surface Height Distribution and Numerical Characteristics during Fatigue

The random surface topography is determined by the height distribution with position coordinates as variables. It is difficult to describe surface topography using strict mathematical expressions. However, the statistical parameters of surface topography can be used for the evaluation of surface topography. Statistical frequency distribution refers to the interval division of each observation value according to a certain interval length so that each observation value forms a different interval distribution. Due to the difference in the nature of objective phenomena, their frequency distribution curves will show different shapes. The height distribution is a reflection of the surface profile shape, which plays an important role in the study of the size, shape, and distribution of the surface profile. Therefore, the surface height distribution can be adopted to statistically describe the changes in metal surface topography during fatigue.Normal distribution, also known as Gaussian distribution, is the most common probability distribution of continuous random variables. The common physical quantity that obeys normal distribution includes the measurement error of object size, the measurement of biological shape, and the quality index of products under normal production conditions. On the basis of normal distribution, combined with the specific parameter of surface height data, the statistical analysis of surface height distribution in the metal fatigue process is carried out. The details are as follows.If the distribution of surface height (z) of the fatigue process satisfies(8)fz=12πσe−z−μ2/2σ2,−∞<z<+∞,then the metal surface height in the fatigue process obeys the normal distribution, which is recorded as follows:(9)z∼Nμ,σ2.In equation (9), μ is the mathematical expectation (mean value) and the position parameter, which determines the location of the distribution, that is, the dominant numerical value of the surface height in the process of metal fatigue. σ (the standard deviation) is the scale parameter. When μ is a certain value, the shape of the curve is determined by σ, and the smaller the σ value is, the higher and thinner the curve is, which indicates that the metal surface height distribution is more concentrated, the uniformity is better, and the surface is smoother. The larger the σ value is, the more squat the curve is, which indicates that the metal surface height distribution is more dispersed, the uniformity is worse, and the surface is more rough and uneven. In addition, the amplitude A of surface height distribution is the frequency corresponding to the highest point of the normal distribution curve of surface height.


## 2.4. Surface Topography Feature Fusion Method

Naive Bayesian classifier is a simple and effective linear classifier. When the features are independent of each other, especially in the case of small-scale samples, the Bayesian model can show a good classification effect because of the small amount of sample data, and the feature variables are independent of each other, so the Bayesian method is applied for surface feature fusion.There are three kinds of Bayesian models: the polynomial model, Gaussian model, and Bernoulli model. The polynomial model and Bernoulli model are suitable for discrete data variables, while the Gaussian model is suitable for continuous data variables. In this paper, the Gaussian model is used as the fusion classification model of sample surface features.The Naive Bayesian method is a learning method using the Bayesian formula. The Bayesian method forms prior distribution by quantifying prior information and applies prior distribution to statistical inference to improve the accuracy of statistical inference. Naive means that it assumes the components of input variables (vectors) are independent of each other:(10)PXY=PX/Y∗PYPX,where P(Y) is the prior probability, PX/Y is conditional probability, obtained through historical data statistics, and PY/X is the posterior probability, which is the target.The Gaussian model is a kind of the Naive Bayesian model. It assumes that each featurexi obeys normal distribution under each category y. The probability density function of normal distribution is adopted to calculate the probability in the algorithm as follows:(11)Pxici=12πσ2exp−xi−μ22σ2,where μ represents the mean value of feature xi in the sample of category ci, i = 1, 2, 3, …:(12)μ=∑inxiN.σ2 stands for the variance of feature xi in the sample of category y:(13)σ2=∑inxi−μ2N,whereN is the number of training samples.It is known from the Naive Bayesian formula that, for the posterior probability values of different categoriesci, the denominator of the Bayesian formula is always the same, so only the numerator is required to be solved. Taking category c1 as an example, the posterior probability formula is as follows:(14)Pc1|xi=Pc1⋅∏i=1nPxi/c1i=1,2,3....In the same way, the posterior probabilities of other categories can be obtained.As described in Sections2.1–2.3, the five features describing the surface topography are arithmetic mean deviation Sa, surface height distribution amplitude A, mean value μ, standard deviation σ, and multifractal spectral width Δα, respectively. The above features correspond to xi in equation (11). The features are independent of each other and meet the requirements of the Naive Bayesian method. The feature data are all continuous variables, so the Gaussian model is selected as the fusion model. In the Gaussian model, the mean and variance of fatigue damage states (i.e., cyclic loading cycles, corresponding to ci of equation (11)) are calculated, and then, the probability density function of normal distribution is obtained. According to the probability density function, the probability density value of continuous test data is calculated to obtain the prior probability. The posterior probability value corresponding to each fatigue damage state is obtained by substituting the probability density value into the posterior probability formula equation (14). The fatigue damage state corresponding to the maximum posterior probability value is the prediction result.


## 3. An Application Case Study

### 3.1. Experimental Materials

Q235 steel specimens (GB/T 3075–2008) are selected as the experimental objects. The specimens are obtained by laser cutting of the same batch of materials. The geometric size of the specimen (unit: mm) is shown in Figure2.Figure 2
Specimen size (unit: mm).The material composition is as follows: C:0.22%, Mn:0.162%, Co:0.10%, Cr:0.01%, Zn:0.015%, Pb:0.062%, Mo:0.003%, P:0.112%, S:0.101%, Si:0.411%, K:0.518%, and Fe: Bal. The tensile strength is 383 MPa.


### 3.2. Experimental Process

The experimental platform (Figure3) consists of an in situ tensile tester (Mtest5000-F-K, Sinotest Equipment Co., Ltd.) and a three-dimensional profilometer (VR3200, Keyence). The in situ tensile tester is loaded using a piezoelectric stacking device with a deformation measurement accuracy of 1 μm, displacement resolution of 0.01 μm, and test force error of ±1% of the indicated value. The 3D profilometer is based on the principle of grating projection for three-dimensional measurement, with a measurement accuracy of 3 μm in height, 2 μm in width, and a repeatability of 0.5 μm in both height and width.Figure 3
The experimental platform.Four specimens were selected for the experiment, numbered 1, 2, 3, and 4. The fiber laser marking machine (JD-GQ20) was used for treating the surface of specimens in the raster scanned format to keep the surface state consistent. The processing area is shown in Figure4.Figure 4
Dimensions of the processing area of the specimen (unit: mm).Low-cycle fatigue tests were carried out on four specimens by the in situ tensile tester. The loading conditions are as follows: constant frequency 13 Hz, loading average value 740 N, and amplitude 60 N. In the process of the fatigue test, under the same conditions, the surface topographic characteristics of four specimens were obtained after 0, 2000, 4000, 6000, and 8000 cycles and fracture, respectively. It is found that all the specimens are fractured within the range of 9700 ± 200 cycles.In order to distinguish from surface waviness and surface shape error when collecting 3D surface topography, the standard Gaussian filtering method is used to remove components with the wave distance greater than 1 mm. When extracting surface topography data, data points are taken from the center of the specimen outward to form a 128∗128 data region. The distance between any two data points is 7.39μm, and the region size is 0.939 mm × 0.939 mm.


## 3.1. Experimental Materials

Q235 steel specimens (GB/T 3075–2008) are selected as the experimental objects. The specimens are obtained by laser cutting of the same batch of materials. The geometric size of the specimen (unit: mm) is shown in Figure2.Figure 2
Specimen size (unit: mm).The material composition is as follows: C:0.22%, Mn:0.162%, Co:0.10%, Cr:0.01%, Zn:0.015%, Pb:0.062%, Mo:0.003%, P:0.112%, S:0.101%, Si:0.411%, K:0.518%, and Fe: Bal. The tensile strength is 383 MPa.


## 3.2. Experimental Process

The experimental platform (Figure3) consists of an in situ tensile tester (Mtest5000-F-K, Sinotest Equipment Co., Ltd.) and a three-dimensional profilometer (VR3200, Keyence). The in situ tensile tester is loaded using a piezoelectric stacking device with a deformation measurement accuracy of 1 μm, displacement resolution of 0.01 μm, and test force error of ±1% of the indicated value. The 3D profilometer is based on the principle of grating projection for three-dimensional measurement, with a measurement accuracy of 3 μm in height, 2 μm in width, and a repeatability of 0.5 μm in both height and width.Figure 3
The experimental platform.Four specimens were selected for the experiment, numbered 1, 2, 3, and 4. The fiber laser marking machine (JD-GQ20) was used for treating the surface of specimens in the raster scanned format to keep the surface state consistent. The processing area is shown in Figure4.Figure 4
Dimensions of the processing area of the specimen (unit: mm).Low-cycle fatigue tests were carried out on four specimens by the in situ tensile tester. The loading conditions are as follows: constant frequency 13 Hz, loading average value 740 N, and amplitude 60 N. In the process of the fatigue test, under the same conditions, the surface topographic characteristics of four specimens were obtained after 0, 2000, 4000, 6000, and 8000 cycles and fracture, respectively. It is found that all the specimens are fractured within the range of 9700 ± 200 cycles.In order to distinguish from surface waviness and surface shape error when collecting 3D surface topography, the standard Gaussian filtering method is used to remove components with the wave distance greater than 1 mm. When extracting surface topography data, data points are taken from the center of the specimen outward to form a 128∗128 data region. The distance between any two data points is 7.39μm, and the region size is 0.939 mm × 0.939 mm.


## 4. Result Analysis

### 4.1. 3D Surface Topography Evolution and Its Analysis

In order to visually observe the changes of three-dimensional surface topography of specimens subjected to cyclic loading, the three-dimensional surface topography maps of specimens under different fatigue cycles are drawn (taking no. 1 specimen as an example, the number of cycles is 0, 2000, 4000, 6000, 8000, and fracture), as shown in Figure5. It can be seen that the surface of the specimen without cyclic loading is relatively uniform, with fewer prominent peaks and valleys. With the increase of the number of cycles, the density of the surface peaks and valleys gradually increase, and the peak and valley values also increase correspondingly, showing obvious intrusion and extrusion effects.Figure 5
3D surface topography of the specimen during fatigue.After the fatigue fracture of the specimen, macro- and microanalyses of the fracture were performed. Figure6(a) shows the macroscopic morphology of the fatigue fracture, which is composed of the fatigue crack initiation region, stable fatigue crack propagation region, and final fracture region. The crack originates from the position of stress concentration on the surface of the specimen and propagates uniformly to the inside. Due to repeated extrusion and friction, the fatigue crack initiation region is brighter, and the fracture surface of the final fracture region is dark gray and rough fibrous.Figure 6
Macro- and micromorphology of fatigue fracture. (a) The macroscopic morphology of the fatigue fracture. (b) The fatigue crack propagation region (1000X). (c) The final fracture region (1000X). (d) The final fracture region (2000X).
(a)(b)(c)(d)Scanning electron microscopy (Zeiss EVO18) was used to observe the fatigue fracture of the specimen, and SEM images of the fracture were obtained, as shown in Figures6(b)–6(d).Figure6(b) shows the micromorphology of the fatigue crack propagation region, in which there are obvious wavy fatigue striations extending from the initiation region to the interior of the specimen. Figure 6(c) shows the micromorphology of the final fracture region. With the application of cyclic loading, the microcracks continue to expand and connect to form macrocracks, and the fracture morphology is similar to that of static load fracture. There are a large number of dimples in the fracture of the final fracture region, which is ductile fracture. When the fracture is magnified 2000 times, as shown in Figure 6(d), there are second phase particles at the bottom of dimples. The toughness of the second phase particles is poor, which will cause higher stress concentration during crack propagation and may lead to the appearance of secondary cracks.


### 4.2. Change in the Basic Characteristics of 3D Surface Topography during Fatigue

At the initial stage of the metal subjected to cyclic loading, dislocations are arranged in the material matrix as dislocation walls. With the increase of fatigue cycles, the dislocation arrangement becomes orderly; meanwhile, the thin straight and evenly distributed slip lines appear on the crystal grain surface. When the cyclic loading continues to act, the slip line widens and forms slip bands. With the widening of the slip band, the metal surface roughness changes due to the intrusion and extrusion phenomenon.Figure7 shows the variation curve of the surface arithmetic mean deviation Sa during the fatigue process. Sa presents a nonlinear growth trend with the increase of fatigue cycles.Figure 7
Changes ofSa in the fatigue process.


### 4.3. Multifractal Spectrum Analysis

Figure8 presents the multifractal spectrum of four specimens from the initial state to fatigue fracture. The value range of the weight factor q is [−55, 55], and the interval is 1. It can be seen from the figure that all the spectra are convex unimodal functions, indicating that the surface topography of the specimens exhibits various fractal characteristics under different cycles, that is, the metal surface topography has the multifractal characteristic during the fatigue process. Therefore, the multifractal spectrum can be applied to describe the change of metal surface topography.Figure 8
The surface multifractal spectrum of the metal fatigue process. (a) Specimen no. 1. (b) Specimen no. 2. (c) Specimen no. 3. (d) Specimen no. 4.
(a)(b)(c)(d)Multifractal spectrum width Δα is a characteristic parameter of the multifractal spectrum, αmax is the maximum value of the singularity index, and αmin is the minimum value of the singularity index, which corresponds to the minimum and maximum height probabilities, respectively. Δα reflects the surface height fluctuation of the metal fatigue process. The larger the Δα value, the wider the range of height fluctuation. At this time, the surface of the specimen fluctuates violently, and the multifractal spectrum becomes wider. Conversely, the smaller the Δα value, the smaller the height fluctuation range. The specimen surface is flat, and the multifractal spectrum becomes narrower.Figure9 shows the multifractal spectrum width variation curve of four specimens in the fatigue process. It is clear that Δα increases with the increase of fatigue cycles, indicating that the height fluctuation range of surface topography increases with the increment of fatigue cycles, and the surface changes from smooth to rough.Figure 9
Variation curves of multifractal spectrum width.


### 4.4. Height Frequency Distribution and Digital Characteristic Analysis

The metal surface height data collected for each damage state are divided into intervals of 5μm. For example, when the fatigue cycle of no. 1 specimen is 6000, the maximum and minimum surface height data are 15.09 μm and -27.2 μm, respectively. Therefore, the surface height is divided into the following 10 intervals: (−30, −25), (−25, −20), (−20, −15), (−15, −10), (−10, −5), (−5, 0), (0, 5), (5, 10), (10, 15), and (15, 20), and calculate the number of height data in accordance with the abovementioned intervals. The frequency distribution curve of metal surface height is shown in Figure 10. With the increase of fatigue cycles, the concentration of metal surface height data decreases, and the span becomes larger.Figure 10
The surface height frequency distribution curve of the metal fatigue process. (a) Specimen no. 1. (b) Specimen no. 2. (c) Specimen no. 3. (d) Specimen no. 4.
(a)(b)(c)(d)The normality of the surface height distribution of the four specimens during the fatigue process is tested, and it is found that the surface height distribution of each damage state obeys the normal distribution. Figure11 shows the normality test results of the surface height distribution of no. 1 specimen in each damage state.Figure 11
The normality test of surface height distribution during fatigue. (a) 0 cycle. (b) 2000 cycles. (c) 4000 cycles. (d) 6000 cycles. (e) 8000 cycles. (f) Fracture.
(a)(b)(c)(d)(e)(f)Figures12–14 are the variation curves of the numerical characteristics (amplitude A, mean μ, and standard deviation δ) of the surface height distribution of the four specimens during fatigue. The amplitude A decreases with the increase of fatigue cycles, while the mean μ and standard deviation δ increase with the increase of fatigue cycles, indicating that the concentration of metal surface height data decreases during the fatigue process, and the proportion of the main height value is increasing, the surface becomes uneven, and the roughness increases.Figure 12
The variation curves of amplitude (A).Figure 13
The variation curves of mean valueμ.Figure 14
The variation curves of standard deviationδ.


## 4.1. 3D Surface Topography Evolution and Its Analysis

In order to visually observe the changes of three-dimensional surface topography of specimens subjected to cyclic loading, the three-dimensional surface topography maps of specimens under different fatigue cycles are drawn (taking no. 1 specimen as an example, the number of cycles is 0, 2000, 4000, 6000, 8000, and fracture), as shown in Figure5. It can be seen that the surface of the specimen without cyclic loading is relatively uniform, with fewer prominent peaks and valleys. With the increase of the number of cycles, the density of the surface peaks and valleys gradually increase, and the peak and valley values also increase correspondingly, showing obvious intrusion and extrusion effects.Figure 5
3D surface topography of the specimen during fatigue.After the fatigue fracture of the specimen, macro- and microanalyses of the fracture were performed. Figure6(a) shows the macroscopic morphology of the fatigue fracture, which is composed of the fatigue crack initiation region, stable fatigue crack propagation region, and final fracture region. The crack originates from the position of stress concentration on the surface of the specimen and propagates uniformly to the inside. Due to repeated extrusion and friction, the fatigue crack initiation region is brighter, and the fracture surface of the final fracture region is dark gray and rough fibrous.Figure 6
Macro- and micromorphology of fatigue fracture. (a) The macroscopic morphology of the fatigue fracture. (b) The fatigue crack propagation region (1000X). (c) The final fracture region (1000X). (d) The final fracture region (2000X).
(a)(b)(c)(d)Scanning electron microscopy (Zeiss EVO18) was used to observe the fatigue fracture of the specimen, and SEM images of the fracture were obtained, as shown in Figures6(b)–6(d).Figure6(b) shows the micromorphology of the fatigue crack propagation region, in which there are obvious wavy fatigue striations extending from the initiation region to the interior of the specimen. Figure 6(c) shows the micromorphology of the final fracture region. With the application of cyclic loading, the microcracks continue to expand and connect to form macrocracks, and the fracture morphology is similar to that of static load fracture. There are a large number of dimples in the fracture of the final fracture region, which is ductile fracture. When the fracture is magnified 2000 times, as shown in Figure 6(d), there are second phase particles at the bottom of dimples. The toughness of the second phase particles is poor, which will cause higher stress concentration during crack propagation and may lead to the appearance of secondary cracks.


## 4.2. Change in the Basic Characteristics of 3D Surface Topography during Fatigue

At the initial stage of the metal subjected to cyclic loading, dislocations are arranged in the material matrix as dislocation walls. With the increase of fatigue cycles, the dislocation arrangement becomes orderly; meanwhile, the thin straight and evenly distributed slip lines appear on the crystal grain surface. When the cyclic loading continues to act, the slip line widens and forms slip bands. With the widening of the slip band, the metal surface roughness changes due to the intrusion and extrusion phenomenon.Figure7 shows the variation curve of the surface arithmetic mean deviation Sa during the fatigue process. Sa presents a nonlinear growth trend with the increase of fatigue cycles.Figure 7
Changes ofSa in the fatigue process.


## 4.3. Multifractal Spectrum Analysis

Figure8 presents the multifractal spectrum of four specimens from the initial state to fatigue fracture. The value range of the weight factor q is [−55, 55], and the interval is 1. It can be seen from the figure that all the spectra are convex unimodal functions, indicating that the surface topography of the specimens exhibits various fractal characteristics under different cycles, that is, the metal surface topography has the multifractal characteristic during the fatigue process. Therefore, the multifractal spectrum can be applied to describe the change of metal surface topography.Figure 8
The surface multifractal spectrum of the metal fatigue process. (a) Specimen no. 1. (b) Specimen no. 2. (c) Specimen no. 3. (d) Specimen no. 4.
(a)(b)(c)(d)Multifractal spectrum width Δα is a characteristic parameter of the multifractal spectrum, αmax is the maximum value of the singularity index, and αmin is the minimum value of the singularity index, which corresponds to the minimum and maximum height probabilities, respectively. Δα reflects the surface height fluctuation of the metal fatigue process. The larger the Δα value, the wider the range of height fluctuation. At this time, the surface of the specimen fluctuates violently, and the multifractal spectrum becomes wider. Conversely, the smaller the Δα value, the smaller the height fluctuation range. The specimen surface is flat, and the multifractal spectrum becomes narrower.Figure9 shows the multifractal spectrum width variation curve of four specimens in the fatigue process. It is clear that Δα increases with the increase of fatigue cycles, indicating that the height fluctuation range of surface topography increases with the increment of fatigue cycles, and the surface changes from smooth to rough.Figure 9
Variation curves of multifractal spectrum width.


## 4.4. Height Frequency Distribution and Digital Characteristic Analysis

The metal surface height data collected for each damage state are divided into intervals of 5μm. For example, when the fatigue cycle of no. 1 specimen is 6000, the maximum and minimum surface height data are 15.09 μm and -27.2 μm, respectively. Therefore, the surface height is divided into the following 10 intervals: (−30, −25), (−25, −20), (−20, −15), (−15, −10), (−10, −5), (−5, 0), (0, 5), (5, 10), (10, 15), and (15, 20), and calculate the number of height data in accordance with the abovementioned intervals. The frequency distribution curve of metal surface height is shown in Figure 10. With the increase of fatigue cycles, the concentration of metal surface height data decreases, and the span becomes larger.Figure 10
The surface height frequency distribution curve of the metal fatigue process. (a) Specimen no. 1. (b) Specimen no. 2. (c) Specimen no. 3. (d) Specimen no. 4.
(a)(b)(c)(d)The normality of the surface height distribution of the four specimens during the fatigue process is tested, and it is found that the surface height distribution of each damage state obeys the normal distribution. Figure11 shows the normality test results of the surface height distribution of no. 1 specimen in each damage state.Figure 11
The normality test of surface height distribution during fatigue. (a) 0 cycle. (b) 2000 cycles. (c) 4000 cycles. (d) 6000 cycles. (e) 8000 cycles. (f) Fracture.
(a)(b)(c)(d)(e)(f)Figures12–14 are the variation curves of the numerical characteristics (amplitude A, mean μ, and standard deviation δ) of the surface height distribution of the four specimens during fatigue. The amplitude A decreases with the increase of fatigue cycles, while the mean μ and standard deviation δ increase with the increase of fatigue cycles, indicating that the concentration of metal surface height data decreases during the fatigue process, and the proportion of the main height value is increasing, the surface becomes uneven, and the roughness increases.Figure 12
The variation curves of amplitude (A).Figure 13
The variation curves of mean valueμ.Figure 14
The variation curves of standard deviationδ.


## 5. Topography Feature Data Fusion and Fatigue Damage Assessment Based on Bayesian Estimation

### 5.1. Prior Probability Analysis

Through the low-cycle fatigue test, 24 sets of data with 6 damage states (0, 2000, 4000, 6000, 8000 cycles, and fracture) were obtained, of which 18 sets of data were as training data and 6 sets of data were as test data. Taking 0 cycle as an example, the mean and variance (i.e., prior data) of its characteristic parameters were calculated, as shown in Table1. The prior data of each damage state can also be obtained. Figure 15 shows the probability density diagram of each characteristic parameter under different fatigue damage states.Table 1
Mean and variance values of each characteristic parameter (0 cycle).
Characteristic parameterMeanVarianceμ0.22941.57 × 10−4δ3.59361.67 × 10−2A130132.22 × 105Sa1.67891.02 × 10−2Δα0.98396.32 × 10−6Figure 15
The priori probability distribution of each characteristic parameter under different cycles. (a)μ. (b) δ. (c) A. (d) Sa. (e) Δα.
(a)(b)(c)(d)(e)


### 5.2. Posterior Probability Analysis

Table2 shows the 6 sets of test data for each fatigue damage state, and the posterior probability values of the test data were obtained based on the prior probability from the training data (Table 3).Table 2
Six sets of test data values for each damage state.
Fatigue cycleμδASaΔα00.2243.487134001.61520.981320000.90224.309103702.07770.9940001.5315.83575983.73571.036560002.728.37551045.94851.117580006.28510.1538989.50931.22Fracture12.4219.1937231414.59851.4112Table 3
Posterior probability values of test data.
Test dataPredicted loading cycles0 cycle2000 cycles4000 cycles6000 cycles8000 cyclesFacture02.146 × 10−28.75 × 10−1612.624 × 10−7498.717 × 10−13959.498 × 10−4302≈020001.779 × 10−126.915 × 10−51.425 × 10−151.0 × 10−911.96 × 10−3919.313 × 10−149040002.06 × 10−3626.656 × 10−2124.947 × 10−34.33 × 10−4314.436 × 10−2808≈060002.13 × 10−7134.3 × 10−4559.0 × 10−1424.776 × 10−38.35 × 10−3084.799 × 10−190180003.9 × 10−7272.312 × 10−4001.558 × 10−1782.15 × 10−492.59 × 10−71.196 × 10−81Fracture≈03.57 × 10−34591.19 × 10−32312.83 × 10−17432.877 × 10−6555.78 × 10−4From the posterior probability of the test data in Table3, a series of posterior probability values were obtained for each group of test data, and the prediction loading cycles corresponding to the maximum value were taken as the prediction result. It is found that the predicted results matched the actual categories, and the Naive Bayesian method could be used for the data fusion of surface topography features in the fatigue damage process.


## 5.1. Prior Probability Analysis

Through the low-cycle fatigue test, 24 sets of data with 6 damage states (0, 2000, 4000, 6000, 8000 cycles, and fracture) were obtained, of which 18 sets of data were as training data and 6 sets of data were as test data. Taking 0 cycle as an example, the mean and variance (i.e., prior data) of its characteristic parameters were calculated, as shown in Table1. The prior data of each damage state can also be obtained. Figure 15 shows the probability density diagram of each characteristic parameter under different fatigue damage states.Table 1
Mean and variance values of each characteristic parameter (0 cycle).
Characteristic parameterMeanVarianceμ0.22941.57 × 10−4δ3.59361.67 × 10−2A130132.22 × 105Sa1.67891.02 × 10−2Δα0.98396.32 × 10−6Figure 15
The priori probability distribution of each characteristic parameter under different cycles. (a)μ. (b) δ. (c) A. (d) Sa. (e) Δα.
(a)(b)(c)(d)(e)


## 5.2. Posterior Probability Analysis

Table2 shows the 6 sets of test data for each fatigue damage state, and the posterior probability values of the test data were obtained based on the prior probability from the training data (Table 3).Table 2
Six sets of test data values for each damage state.
Fatigue cycleμδASaΔα00.2243.487134001.61520.981320000.90224.309103702.07770.9940001.5315.83575983.73571.036560002.728.37551045.94851.117580006.28510.1538989.50931.22Fracture12.4219.1937231414.59851.4112Table 3
Posterior probability values of test data.
Test dataPredicted loading cycles0 cycle2000 cycles4000 cycles6000 cycles8000 cyclesFacture02.146 × 10−28.75 × 10−1612.624 × 10−7498.717 × 10−13959.498 × 10−4302≈020001.779 × 10−126.915 × 10−51.425 × 10−151.0 × 10−911.96 × 10−3919.313 × 10−149040002.06 × 10−3626.656 × 10−2124.947 × 10−34.33 × 10−4314.436 × 10−2808≈060002.13 × 10−7134.3 × 10−4559.0 × 10−1424.776 × 10−38.35 × 10−3084.799 × 10−190180003.9 × 10−7272.312 × 10−4001.558 × 10−1782.15 × 10−492.59 × 10−71.196 × 10−81Fracture≈03.57 × 10−34591.19 × 10−32312.83 × 10−17432.877 × 10−6555.78 × 10−4From the posterior probability of the test data in Table3, a series of posterior probability values were obtained for each group of test data, and the prediction loading cycles corresponding to the maximum value were taken as the prediction result. It is found that the predicted results matched the actual categories, and the Naive Bayesian method could be used for the data fusion of surface topography features in the fatigue damage process.


## 6. Discussion

By analyzing the data of 3D surface topography of the Q235 steel specimen during the fatigue damage process, it is found that the statistical characteristics of the above 3D topography, such as arithmetic mean deviationSa, surface height distribution mean μ, and standard deviation δ. show obvious regular changes with the increase of fatigue cycles. Although the experiment was carried out on a Q235 material, the feature extraction method and corresponding rules can also be applied to other carbon steel materials, such as 45 steel and Q345.In the process of low-cycle fatigue damage, significant plastic deformation will be generated during each cycle of loading, and its fatigue behavior can be related to plastic deformation and limited fatigue cycles. The fatigue damage process consists of three stages: microstructure evolution, crack initiation, and propagation, and the crack initiation life accounts for the majority of the total life. The low-cycle fatigue cracks preferentially nucleate in the persistent slip band (PSB) with the highest slip inhomogeneity and maximum strain localization. The slip produced by the PSB in the surface grains can be transferred in a large amount, and the phenomenon of material intrusion and extrusion caused by this slip band can be observed on the free surface of the specimen. So far, the mechanism of crack initiation caused by this intrusion and extrusion is still unclear, which may be related to the dynamic equilibrium between dislocation proliferation and annihilation in the PSB. The whole process of crack initiation, propagation, and fracture was described by fatigue fracture analysis, such as the nucleation of microcracks at the specimen surface, fatigue striations caused by crack propagation, dimple morphology, and second phase particles in the transient fracture region. The nature of fracture is verified, and the necessary means are provided to reveal the fracture mechanism.This paper presents the statistical results of surface topography features and describes the quantitative relationship between them and fatigue cycles, aiming to establish a data fusion assessment model for fatigue damage by constructing a nonlinear mapping between surface topography features and damage. In this way, a feasible operational method for fatigue damage assessment in engineering is sought to avoid the current situation where the micromechanism (such as crack nucleation, vacancy dipole, and interface dislocation.) of the fatigue damage process is unclear and difficult to quantify the damage. It is becoming a hot issue to correlate the microscale deformation mechanism of metal materials/components with fatigue failure at the macroscale.


## 7. Conclusion

Q235 steel specimens were selected for tensile fatigue tests, and the three-dimensional surface topography information was extracted. The characteristics of metal surface topography during the fatigue process, such as arithmetic mean deviationSa, surface height distribution amplitude A, mean μ, standard deviation δ, and multifractal spectrum width Δα, were obtained.It is found thatSa shows a nonlinear growth trend with the increase of fatigue cycles. Due to the initiation and expansion of cracks in the middle and late stages of fatigue, the prominent peaks and valleys are higher and deeper, and the growth rate of Sa becomes faster. In the digital features of surface height distribution, the amplitude A decreases with the increase of fatigue cycles, and the mean value μ and standard deviation δ increase with the increase of fatigue cycles. The multifractal spectrum width Δα increases with the increment of fatigue cycles, indicating that the fluctuation range of metal surface height data increases with the increase of fatigue damage, and the surface becomes uneven.The arithmetic mean deviationSa, digital features of surface height distribution (amplitude A, mean μ, and standard deviation δ), and multifractal spectral width Δα were taken as the multifeature quantities of surface topography. The data fusion of the above multifeature quantities was performed based on the Naive Bayesian algorithm, and its validity was preliminarily verified.This study can provide methodological support for fatigue damage assessment of mechanical components. However, the amount of experimental data in the modeling process is still small, which needs to be increased. In addition, in order to improve the prediction accuracy and reduce the uncertainty, the subsequent research will conduct fusion assessment based on topographic features and combined with multiple physical parameters, which is also an important direction for future nondestructive testing.



---

*Source: 6680855-2021-01-27.xml*