Damage Diagnosis of Reactive Powder Concrete under Fatigue Loading Using 3D Laser Scanning Technology

Abstract: Damage mechanisms of Reactive Powder Concrete (RPC) under fatigue loading are investigated using the 3D laser scanning technology. An independently configured 3D laser scanning system is used to monitor the damaging procedure. Texture analysis technique is also applied to enhance the understanding of the damage mechanisms of RPC under fatigue loading. In order to obtain the characteristic parameters of point cloud data, a point clouds projection algorithm is proposed. Damage evolution is described by the change of point cloud data of the damage in the 2D plane and 3D space during fatigue loading. The Gray Level Co-occurrence Matrix (GLCM) method is used to extract the characteristic parameters to evaluate the statue of the structural. Angular Second Moment and Cluster Shadow of typical sensitive characteristic indexes is screened by using the Digital Feature Screening. The reliability of the damage index was verified by image texture analysis and data expansion. Indexes extracted in this paper can be used as a new structural health monitoring indicator to assess health condition.


Introduction
In recent years, considerable attention has been paid to the health monitoring of advanced composite structures. Damage identification through various non-destructive testing techniques predict or diagnose the health of composite structures [1][2][3]. Reactive Powder Concrete (RPC), an emerging materials, are widely concerned in railway engineering, roads, buildings, airports and nuclear industrial facilities-ultra high strength, high toughness and high durability [4][5][6][7][8][9]. However, there are few researches on damage identification and detection of composite structures using RPC as raw materials. It is critical to popularize intelligent damage monitoring technology of RPC structure damage.
Certain damage, which have considerable impact on structure reliability and security, are inevitably retained because of the influence of fatigue load [10][11][12]. The process of causing damage caused by cyclic loading includes the combined results of the cumulative process of crack initiation, expansion and final fracture. When the RPC structure is subjected to a continuous fatigue load, initial microcracks appear on the surface of the structure. The stress generated by the matrix crack increases and continues to aggregate, expanding and eventually forming macroscopic cracks. Scholars evaluate the health status of the structure by monitoring the development of damage [13][14][15]. Finding damage in time and predicting injuries is the key to avoiding temporary accidents. Therefore, it is very necessary to fully grasp the development of crack damage as well as find an intelligent method to quickly and conveniently monitor the health of the structure.
plate. The sensor scan plane is parallel to the ground and perpendicular to the drive direction. The parameters of the laser sensor are set on the surface of the sensor control interface. The lifting speed of the scanning optical center is consistent with the motion output of the linear sliding group.

Principle of GLCM
Texture is an important image function, and because of its powerful discrimination ability, it can be used as certainly effective tool for pattern recognition technology. The most widely applied texture analysis technology is the Gray Level Co-occurrence Matrix (GLCM), in which texture features are extracted from the co-occurrence matrices Pd, θ of gray values f(i) and gray values f(j) based on statistical methods (Haralick, Shanmugam, & Dinstein 1973). The principle diagram of gray symbiosis matrix, as shows in Figure 2. The following explains the GLCM schematic in Figure 2. First, input an image. Then a standard window is specified, and the probability of pixel pairs appearing under the same specified conditions is calculated. The probability distribution values are counted to generate a number matrix. Finally, the matrix is converted to a grayscale image. Additionally, the feature parameters characterizing the image properties are extracted. The characteristic parameters extracted in this study include: Angular Second Moment (ASM); Entropy (ENT); Variance (VAR); Correlation (COR); inverse difference moment (IDM); Cluster shadow (CLS), as shown in Table 1.

Principle of GLCM
Texture is an important image function, and because of its powerful discrimination ability, it can be used as certainly effective tool for pattern recognition technology. The most widely applied texture analysis technology is the Gray Level Co-occurrence Matrix (GLCM), in which texture features are extracted from the co-occurrence matrices P d , θ of gray values f(i) and gray values f(j) based on statistical methods (Haralick, Shanmugam, & Dinstein 1973). The principle diagram of gray symbiosis matrix, as shows in Figure 2. plate. The sensor scan plane is parallel to the ground and perpendicular to the drive direction. The parameters of the laser sensor are set on the surface of the sensor control interface. The lifting speed of the scanning optical center is consistent with the motion output of the linear sliding group.

Figure 1.
Design drawing of 3D laser scanning system.

Principle of GLCM
Texture is an important image function, and because of its powerful discrimination ability, it can be used as certainly effective tool for pattern recognition technology. The most widely applied texture analysis technology is the Gray Level Co-occurrence Matrix (GLCM), in which texture features are extracted from the co-occurrence matrices Pd, θ of gray values f(i) and gray values f(j) based on statistical methods (Haralick, Shanmugam, & Dinstein 1973). The principle diagram of gray symbiosis matrix, as shows in Figure 2. The following explains the GLCM schematic in Figure 2. First, input an image. Then a standard window is specified, and the probability of pixel pairs appearing under the same specified conditions is calculated. The probability distribution values are counted to generate a number matrix. Finally, the matrix is converted to a grayscale image. Additionally, the feature parameters characterizing the image properties are extracted. The characteristic parameters extracted in this study include: Angular Second Moment (ASM); Entropy (ENT); Variance (VAR); Correlation (COR); inverse difference moment (IDM); Cluster shadow (CLS), as shown in Table 1. The following explains the GLCM schematic in Figure 2. First, input an image. Then a standard window is specified, and the probability of pixel pairs appearing under the same specified conditions is calculated. The probability distribution values are counted to generate a number matrix. Finally, the matrix is converted to a grayscale image. Additionally, the feature parameters characterizing the image properties are extracted. The characteristic parameters extracted in this study include: Angular Second Moment (ASM); Entropy (ENT); Variance (VAR); Correlation (COR); inverse difference moment (IDM); Cluster shadow (CLS), as shown in Table 1. p 2 (i, j, d, θ) Gray distribution uniformity and texture thickness.
Texture period size. Note: , k = 2, 3 . . . , 2g. Table 1 introduces formulas and damage characteristics of the parameters. The direction, trend and complexity of fracture damage in gray space were quantitatively described by feature parameters

Bending Fatigue Test
Under the repeated load during the actual engineering, the initial microscopic pores and cracks in the structure generate new micro-cracks. At the same time, the continuous development of the original pores and micro-cracks reduces the durability of the structure and causes fatigue damage. Under cyclic loading, the bending fatigue failure mode of reactive powder concrete is manifested as the failure mode of single critical fatigue main crack. The laser scanning system captures the failure process of the fatigue main crack, monitors the damage, and predicts the structural damage pattern. Therefore, a three-point bending fatigue loading test was carried out to verify the reliability of the 3D laser scanning system.

Material
PO42.5 ordinary Portland cement is used for cement, and its performance index is shown in Table 2. The steel fibers used had a diameter of 0.22 mm, a fiber length of 13 mm, and a tensile strength of more than 2850 MPa. A KY-1 type polycarboxylic acid high performance water reducing agent with a water reduction rate of 35% is selected. The admixture is made of grade I fly ash and micro silica.

Production of Test Pieces
According to the "Standards for Active Powder Concrete" GB/T 31387-2015, the mix ratio is shown in Table 3. The specifications and quantities of the test pieces are shown in Table 4. Under the standard curing conditions, the specimens of the specimens were removed after 24 h of curing. Then, after curing at 90 • C for 72 h, it was allowed to stand under standard curing conditions for 28 days.

Experiment
The RPC bending fatigue specimens are three-point bending loads. The fatigue testing machine uses IPC Global UTM-30 dynamic hydraulic servo multi-function material testing machine. The effective span of the test beam is 120 mm, as shown in Figure 3a. The actual loading situation is shown in Figure 3b. The equal-amplitude sine wave load is selected, the load frequency is 15 Hz, and the fatigue load is performed in a stress control manner.

Production of test pieces
According to the "Standards for Active Powder Concrete" GB / T 31387-2015, the mix ratio is shown in Table 3. The specifications and quantities of the test pieces are shown in Table 4. Under the standard curing conditions, the specimens of the specimens were removed after 24 hours of curing. Then, after curing at 90℃ for 72 hours, it was allowed to stand under standard curing conditions for 28 days.

Experiment
The RPC bending fatigue specimens are three-point bending loads. The fatigue testing machine uses IPC Global UTM-30 dynamic hydraulic servo multi-function material testing machine. The effective span of the test beam is 120 mm, as shown in Figure 3(a). The actual loading situation is shown in Figure 3(b). The equal-amplitude sine wave load is selected, the load frequency is 15 Hz, and the fatigue load is performed in a stress control manner.

Result and analysis
The whole process of damage mechanism of RPC under fatigue loading are investigated using the 3D laser scanning technology as show in Figure 4.

Result and Analysis
The whole process of damage mechanism of RPC under fatigue loading are investigated using the 3D laser scanning technology as show in Figure 4.
The overall process consists of six parts in Figure 4. The establishment of a 3D laser scanning system. The execution of RPC fatigue loading test. Further, the GLCM is generated. Including extracted feature parameters, parameter screening and verification processes. Stack each part as the main line of research. And mark the core content of each part.   The overall process consists of six parts in Figure 4. The establishment of a 3D laser scanning system. The execution of RPC fatigue loading test. Further, the GLCM is generated. Including extracted feature parameters, parameter screening and verification processes. Stack each part as the main line of research. And mark the core content of each part.

Bending strength
Four sets of test pieces were selected for the static load bending test, and the average value of the 28-day bending strength was measured, as shown in Table 5. There was no significant change in the surface of the test piece before the test loading. The initial crack occurs in the test piece when the loading force reaches 70% to 80% of the ultimate bending strength. The flexural strength value did not decrease as the specimen appeared across the crack. It was not until the steel fiber in the cracked cross-section was pulled out that the test specimen was completely damaged.

Bending Strength
Four sets of test pieces were selected for the static load bending test, and the average value of the 28-day bending strength was measured, as shown in Table 5. There was no significant change in the surface of the test piece before the test loading. The initial crack occurs in the test piece when the loading force reaches 70% to 80% of the ultimate bending strength. The flexural strength value did not decrease as the specimen appeared across the crack. It was not until the steel fiber in the cracked cross-section was pulled out that the test specimen was completely damaged.

Data Collection
The structure is located in the scanning sector detection area by adjusting the horizontal position of the holder body. In order to ensure complete structural modeling information, the effective detection range in the vertical direction is greater than the height of the detected object. The data acquisition process as shown in Figure 8. Figure 5 shows the process of obtaining a structural point cloud image. Verification of scanning system accuracy based on bending fatigue load test. A total of 100 scans were performed in each group. Figure 6 is a structural point cloud diagram. The structure is located in the scanning sector detection area by adjusting the horizontal position of the holder body. In order to ensure complete structural modeling information, the effective detection range in the vertical direction is greater than the height of the detected object. The data acquisition process as shown in Figure 8.

Damage parameter acquisition
When the damage degree of the RPC structure is slight, the displacement and deflection of the damage source area are small. However, the location and elevation information of the damaged area is constantly changing as the degree of damage continues to develop. This change can be obtained in  The structure is located in the scanning sector detection area by adjusting the horizontal position of the holder body. In order to ensure complete structural modeling information, the effective detection range in the vertical direction is greater than the height of the detected object. The data acquisition process as shown in Figure 8.

Damage parameter acquisition
When the damage degree of the RPC structure is slight, the displacement and deflection of the damage source area are small. However, the location and elevation information of the damaged area is constantly changing as the degree of damage continues to develop. This change can be obtained in

Damage Parameter Acquisition
When the damage degree of the RPC structure is slight, the displacement and deflection of the damage source area are small. However, the location and elevation information of the damaged area is constantly changing as the degree of damage continues to develop. This change can be obtained in two ways. On the one hand, the continuous development of cracks in 2D plane can be achieved. On the other hand, the displacement of the crack and the change in height information in the 3D space.
The 2D image can be captured by the CCD camera built in the 2D laser ranging sensor. According to the principle of the independently developed 3D laser scanning system, the point cloud data for obtaining the scanned test piece is a digital matrix composed of data points having position information.
This value can be calculated by the Formula (2). Therefore, the damage region of the digital matrix composed of the 2D image and the point cloud elevation information is selected as the Region of Interest (ROI), respectively. The GLCM of the ROI is constructed separately, the feature parameters are extracted, and the texture features are calculated therefrom. The angles are determined as 0 • , 45 • , 90 • , and 135 • , the image grayscale is selected as 265, and the distance between the pixel pairs is 2 pixels. Therefore, when the construction factor is d = 2, g = 2 8 and θ selects 0 • , 45 • , 90 • and 135 • , and six characteristic parameters of the damaged area of the test piece are extracted. In order to satisfy the image rotation invariance, the average of the 4 angles is taken as the final value function parameter.

Parameter Extraction in 2D Plane
The images of the fixed damage area are acquired with the same number of frames, as shown in Figure 7.
two ways. On the one hand, the continuous development of cracks in 2D plane can be achieved. On the other hand, the displacement of the crack and the change in height information in the 3D space.
The 2D image can be captured by the CCD camera built in the 2D laser ranging sensor. According to the principle of the independently developed 3D laser scanning system, the point cloud data for obtaining the scanned test piece is a digital matrix composed of data points having position information. This value can be calculated by the Formula (2). Therefore, the damage region of the digital matrix composed of the 2D image and the point cloud elevation information is selected as the Region of Interest (ROI), respectively. The GLCM of the ROI is constructed separately, the feature parameters are extracted, and the texture features are calculated therefrom. The angles are determined as 0°, 45°, 90°, and 135°, the image grayscale is selected as 265, and the distance between the pixel pairs is 2 pixels. Therefore, when the construction factor is d = 2, g = 2 8 and θ selects 0°, 45°, 90° and 135°, and six characteristic parameters of the damaged area of the test piece are extracted. In order to satisfy the image rotation invariance, the average of the 4 angles is taken as the final value function parameter.

Parameter extraction in 2D plane
The images of the fixed damage area are acquired with the same number of frames, as shown in Figure 7. The ROI of damage images were captured in Figure 7. The number of image acquisition interval frames is 10 images. The direction of the arrow indicates the direction of damage development. As fatigue loads continue to accumulate, damage continues to expand and develop. A nonlinear texture change is formed on a 2D plane. The characteristic parameters of each image were acquired and data recorded. The ROI of damage images were captured in Figure 7. The number of image acquisition interval frames is 10 images. The direction of the arrow indicates the direction of damage development. As fatigue loads continue to accumulate, damage continues to expand and develop. A nonlinear texture change is formed on a 2D plane. The characteristic parameters of each image were acquired and data recorded.

Parameter Extraction in 3D Space
A transformation algorithm for transforming 3D scene data into 2D elevation image is proposed based on the principle of orthogonal projection. Figure 8 is a schematic diagram of a conversion algorithm for 3D point cloud data projection elevation data information.

Parameter extraction in 3D space
A transformation algorithm for transforming 3D scene data into 2D elevation image is proposed based on the principle of orthogonal projection. Figure 8 is a schematic diagram of a conversion algorithm for 3D point cloud data projection elevation data information. In Figure 8, the elevation information of the3D image is represented by the Z-axis in the coordinate system. Therefore, the projection point on the projection plane is the elevation data of the image, which can be represented by a matrix Z(x, y) consisting of discrete quantities: where x is the range of values (0, m) and the range of y is (0, n). All points in the Z(x, y) matrix are represented by one byte of 8 bits, converted to any value in the 0-255 interval, based on the one-dimensional sampling theorem and the principle of linear quantization. Therefore, the matrix Z(x, y) is a grayscale image of 256 gradation values. The converted image is named an elevation projection. An elevation digital matrix corresponding to images of a 2D plane is obtained. Then, convert them into projected images in sequentially. Feature parameters of the projected image are extracted separately. Finally, data statistics are performed. In Figure 8, the elevation information of the3D image is represented by the Z-axis in the coordinate system. Therefore, the projection point on the projection plane is the elevation data of the image, which can be represented by a matrix Z(x, y) consisting of discrete quantities:

Damage diagnosis
z(x m , y 1 ) · · · z(x m , y n ) where x is the range of values (0, m) and the range of y is (0, n). All points in the Z(x, y) matrix are represented by one byte of 8 bits, converted to any value in the 0-255 interval, based on the one-dimensional sampling theorem and the principle of linear quantization. Therefore, the matrix Z(x, y) is a grayscale image of 256 gradation values. The converted image is named an elevation projection. An elevation digital matrix corresponding to images of a 2D plane is obtained. Then, convert them into projected images in sequentially. Feature parameters of the projected image are extracted separately. Finally, data statistics are performed.

Damage Diagnosis
For the purpose of obtaining the damage indicators of fatigue damage model, the parameters is screening use the DFS (Digital Feature Screening) method. The DFS method establishes the selection criteria of the damage characteristic index: obtaining the crack image in a fixed area under the step of increasing the loading level, and sequentially performing image numbering. Finally, damage characteristic index of the crack image is obtained. Further, a digital feature screening analysis chart is established, as shown in Figure 9. Figure 9a shows the variation of the characteristic parameters as the damage grows in 3D space. After analysis, it is known that as the fatigue damage increases, the ASM characteristic parameters will continue to increase according to certain rules. The remaining characteristic parameters change irregularly with the development of the fracture, so the remaining characteristic parameters are excluded. In order to monitor the development of fractures in 3D space, ASM characteristic parameters with good positive correlation with fracture development trends were selected. Figure 9(a) shows the variation of the characteristic parameters as the damage grows in 3D space. After analysis, it is known that as the fatigue damage increases, the ASM characteristic parameters will continue to increase according to certain rules. The remaining characteristic parameters change irregularly with the development of the fracture, so the remaining characteristic parameters are excluded. In order to monitor the development of fractures in 3D space, ASM characteristic parameters with good positive correlation with fracture development trends were selected. Figure 9(b) shows a tendency of characteristic parameters to change as the damage advances in the 2D plane. Analysis of the data changes in the graph shows that as the fatigue damage increases, the CLS characteristic parameters continue to decrease. As the crack develops, the other characteristic parameters show different trends. Therefore, it is excluded. In order to monitor the development of two-dimensional spatial crack damage, the CLS characteristic parameters which are negatively correlated with the crack development trend are selected.  Figure 9b shows a tendency of characteristic parameters to change as the damage advances in the 2D plane. Analysis of the data changes in the graph shows that as the fatigue damage increases, the CLS characteristic parameters continue to decrease. As the crack develops, the other characteristic parameters show different trends. Therefore, it is excluded. In order to monitor the development of two-dimensional spatial crack damage, the CLS characteristic parameters which are negatively correlated with the crack development trend are selected. Figure 9a,b respectively count the changes of the crack development of the six characteristic parameters, and obtain the monitoring basis of the 2D plane and 3D space cracks. The characteristic parameters ASM and CLS can quantify the development of cracks in different spaces and fully describe the development of cracks. Therefore, CLS and ASM characteristic parameters can be used as damage indicators for fatigue damage in RPC bending loading test to detect the development of damage.

Experimental Validation
The reliability of the damage indicator can be verified in two ways. Firstly, it is analyzed whether the texture characteristics characterized by the feature parameters are consistent with the development of the damage. Secondly, multiple sets of experiments were performed, and 30 images of continuous lesion areas were randomly acquired for data verification.

Texture Feature Verification
In order to verify the reliability of the damage index, it is analyzed whether the texture characteristics of the texture are consistent with the development trend of the damage. Table 6 shows texture characteristics of damage indicators. Analysis of the texture characteristics of the damage indicators in Table 5 results in the following conclusions. It shows that in the 2D plane, with the continuous development of fatigue damage, the damage image width with damage as the texture element is increasing. In 3D space, the increase in fatigue damage causes the image texture uniformity to continuously decrease. Therefore, the texture characteristics indicated by the damage index are consistent with the continuous accumulation process of fatigue damage. Finally, the reliability of the damage index is verified.

Data Verification
In order to further verify the broad validity of the damage index. The damage index of the randomly acquired damage continuous image is extracted. Additionally, the fitting relationship is established, as show in Figure 10.  Figure 9(a) and Figure 9(b) respectively count the changes of the crack development of the six characteristic parameters, and obtain the monitoring basis of the 2D plane and 3D space cracks. The characteristic parameters ASM and CLS can quantify the development of cracks in different spaces and fully describe the development of cracks. Therefore, CLS and ASM characteristic parameters can be used as damage indicators for fatigue damage in RPC bending loading test to detect the development of damage.

Experimental validation
The reliability of the damage indicator can be verified in two ways. Firstly, it is analyzed whether the texture characteristics characterized by the feature parameters are consistent with the development of the damage. Secondly, multiple sets of experiments were performed, and 30 images of continuous lesion areas were randomly acquired for data verification.

Texture feature verification
In order to verify the reliability of the damage index, it is analyzed whether the texture characteristics of the texture are consistent with the development trend of the damage. Table 6 shows texture characteristics of damage indicators. Analysis of the texture characteristics of the damage indicators in Table 5 results in the following conclusions. It shows that in the 2D plane, with the continuous development of fatigue damage, the damage image width with damage as the texture element is increasing. In 3D space, the increase in fatigue damage causes the image texture uniformity to continuously decrease. Therefore, the texture characteristics indicated by the damage index are consistent with the continuous accumulation process of fatigue damage. Finally, the reliability of the damage index is verified.

Data verification
In order to further verify the broad validity of the damage index. The damage index of the randomly acquired .  There is a specify positive correlation between ASM and the development trend of injury. Additionally, CLS is identify negatively correlated with the development of fatigue damage. CLS and ASM have the same trend and fantastic stability after comparison with multiple experimental images. The verification results show that the CLS and ASM can used as indicators of damage sensitivity characteristics of crack damage monitoring. This indicator can comprehensively monitor the damage status of cracks and accurately predict the health trend of the structure.

Conclusion
The 3D laser scanning technology of RPC under fatigue load conditions can clearly and effectively reveal the point cloud data features during the damage process. This work can provide a basis for the health monitoring of the RPC structure. Some of the conclusions drawn are as follows: 1. The production principle and process of the inventive 3D laser scanning system was described.
Using the system, the damage process of RPC under fatigue load was studied. In the process, the scanned 3D point cloud data simulates the RPC damage process.
2. When the damage degree of the RPC structure is slight, the displacement and deflection of the damage source area are small. However, as the degree of damage continues to develop, the location and elevation information of the damaged area is constantly changing. The nonlinear variation of the damage region is described by the GLCM. And the damage characteristic parameter is extracted. The system can be popularized and applied to the detection and diagnosis of the health status of other composite materials. The damage mechanism is analyzed by the extracted damage index. However, the scanning system is limited to the detection of damaged areas inside the threedimensional space that can be scanned by laser. Therefore, it is necessary to use methods combined There is a specify positive correlation between ASM and the development trend of injury. Additionally, CLS is identify negatively correlated with the development of fatigue damage. CLS and ASM have the same trend and fantastic stability after comparison with multiple experimental images. The verification results show that the CLS and ASM can used as indicators of damage sensitivity characteristics of crack damage monitoring. This indicator can comprehensively monitor the damage status of cracks and accurately predict the health trend of the structure.