Damage Pattern Recognition and Crack Propagation Prediction for Crumb Rubber Concrete Based on Acoustic Emission Techniques

: In this study, the clustering method of the concrete matrix rupture and rubber fracture damages as well as the prediction of the ultimate load of crumb rubber concrete using the acoustic emission (AE) technique were investigated. The loading environment of the specimens was a four-point bending load. Six clustering methods including k-means, fuzzy c-means (FCM), self-organizing mapping (SOM), Gaussian mixture model (GMM), hierarchical model, and density peak clustering method were analyzed; the results illustrated that the density peak clustering has the best performance. Next, the optimal clustering algorithm was used to cluster AE signals so as to study the evolution behavior of different damage modes, and the ultimate load of crumb rubber concrete was predicted by an artiﬁcial neural network. The results indicated that the combination of AE techniques and appropriate clustering methods such as the density peak clustering method and the artiﬁcial neural network could be used as a practical tool for structural health monitoring of crumb rubber concrete.


Introduction
In recent years, the production of automobiles has led to the accumulation of numerous waste tires that are difficult to dispose of, resulting in a large amount of space occupied and a high fire hazard.The impact of waste tires on the environment and how to effectively deal with them has become an important issue.Therefore, scholars in different fields have proposed several methods to recycle waste tires.In addition to the method of retreading and reusing [1], waste tires can also be converted into other degradable or usable polymeric materials through a series of chemical methods [2].Gnanaraj et al. [3] pyrolyzed waste tires and applied the products to energy storage devices such as batteries.Derakhshan et al. [4] found that in wastewater treatment systems, the use of waste tires as carriers of organisms is a promising application.Recently, in the direction of civil engineering, the application of grinding waste rubber tires into particles as aggregates has become increasingly widespread [5][6][7][8][9][10][11][12].Many scholars have investigated the mechanical properties of rubber concrete.According to the experimental results, they have found that the addition of rubber particles has significantly improved the ductility and durability of the concrete, but it results in a significant decrease in strength [13][14][15][16][17]. Based on the current urgency of finding ways to improve strength defects and increasing the recycling rate of waste tires, scholars in related fields have stimulated in-depth research.Follow-up scholars have also studied the fatigue and impermeability behaviors of rubber concrete [18][19][20][21][22]. Currently, crumb rubber concrete has been widely used in tennis courts, park areas, and bridge road constructions.Therefore, it is important to monitor the structure of these materials to avoid large losses due to large damage inside or outside the structure.For the study of material microcracking and fracture monitoring methods, Joosten et al. [23] proposed a new computational method for the simulation of mesh division of composite materials, which allows accurate prediction of crack propagation and fracture load when damage occurs in the material.Wu et al. [24] proposed a damage model for describing material fracture, which achieves the adaptive and capture of arbitrary crack propagation in solid materials.Lin et al. [25] analyzed the connection between mechanical properties and waveforms and frequencies of four rocks during fracture to improve the accuracy of crack prediction in rocks.Recently, in order to understand the development of microcracks and the monitoring of fracture-causing behaviors, Xu et al. [26] have conducted experimental studies based on acoustic emission (AE) technology.
For crumb rubber concrete, there have been some studies on its microcrack and fracture characteristics [27][28][29][30][31], but there are few studies on the monitoring of microcrack propagation and fracture behavior under different damage modes [32][33][34].Hence, if the acoustic emission signals of the microcracks at various stages in different modes can be monitored and a reasonable analysis made, it will bring benefits to practical engineering applications.For this reason, this paper used AE technology to conduct damage monitoring research on crumb rubber concrete.To investigate the different damage patterns of crumb rubber concrete under a four-point bending test, it is critical to categorize different AE signals.In recent years, a variety of different clustering methods have been adopted for the classification of AE signals obtained in the damage and failure process of concrete materials [35][36][37][38][39].Although the types of these methods have increased, their inherent limitations have not been studied in depth.At the same time, no one used the clustering method to classify the AE signals for crumb rubber concrete.Although some researchers used the k-means clustering method to make a good distinction between the tension mode and the shear mode, the acoustic emission data of crumb rubber concrete are relatively discrete, and the k-means clustering method has limitations [40].Therefore, the acoustic emission signal classification performance of different clustering methods needs to be analyzed and selected so as to combine acoustic emission technology to get a dependable and accurate structural health monitoring system.At the same time, a simple and efficient method is needed to predict the failure loads of the crumb rubber concrete by monitoring the crack development in concrete engineering using acoustic emission techniques.While moving ahead with computer technology, the artificial neural network algorithm has emerged.By training a large number of input data, an output result can be obtained, and by comparing it with the real value, it can be found that the calculation result is relatively reliable [41].
The research of this paper consists of three main parts: firstly, six clustering methods for analyzing AE signals are introduced in Section 2, including k-means, FCM, GMM, SOM, hierarchical model, and density peak clustering; secondly, the optimal clustering method (density peak clustering) is derived by comparative analysis in Section 3.1 and used to identify damage patterns of crumb rubber concrete; finally, in Section 3.2, the artificial neural network algorithm with input data as acoustic emission parameters is used in this paper to predict the ultimate load of crumb rubber concrete for the first time.

Materials and Apparatus
In the four-point bending test, the raw materials of the specimens included ordinary Portland cement, sand, coarse aggregate, and rubber crumb, among which the maximum particle size of coarse aggregate was 20 mm, and the particle size of rubber crumb was between 1 and 2 mm.The mixing ratio of the respective components can be found from our previous work [42].The content of crumb rubber was divided into four groups of 0%, 5%, 10%, and 15% (according to the principle of volume percentage method).The size of the specimen was 400 × 100 × 100 mm.The four sets of specimens were all designed to C30 compressive strength level, which is between 35 and 45 MPa. Figure 1 and Table 1 show the specific information on the location of the crumb rubber concrete specimen and acoustic emission sensors.
our previous work [42].The content of crumb rubber was divided into four groups of 0 5%, 10%, and 15% (according to the principle of volume percentage method).The size the specimen was 400 × 100 × 100 mm.The four sets of specimens were all designed to C compressive strength level, which is between 35 and 45 MPa. Figure 1 and Table 1 sho the specific information on the location of the crumb rubber concrete specimen and acou tic emission sensors.For the apparatus in this test (Table 2), the 100 kN range of electronic universal testi machine and the Micro-II Digital AE System were adopted (Figure 2a).R6α resonant-ty sensors were used.In order to ensure that the sensors were absolutely adhered closely the surface of the specimen, a couplant was used and fixed with a specific fixture (Figu 2b).The threshold for acoustic emission detection was set to 40 dB and the loading rate the specimen was 0.167 kN/s.After setting-up the test apparatus, the processes of d acquisition are: first, during the four-point bending load, the crumb rubber concrete sp imen generates internal cracks, thus forming stress waves (AE source); second, the A sensors collect the propagated AE stress waves and convert them; third, the preamplif amplifies the AE signals, thus ensuring the appropriate signal strength; fourth, the A signals are stored and analyzed by the data acquisition system; and, finally, the AE ch acteristic parameters obtained are analyzed and studied, such as hit, amplitude, cou duration, etc.In addition, the multichannel monitoring system AE win was used as the so ware for AE signal acquisition and analysis, and the Matlab 10 Neural Network Toolb was used for subsequent processing of the data.For the apparatus in this test (Table 2), the 100 kN range of electronic universal testing machine and the Micro-II Digital AE System were adopted (Figure 2a).R6α resonanttype sensors were used.In order to ensure that the sensors were absolutely adhered closely to the surface of the specimen, a couplant was used and fixed with a specific fixture (Figure 2b).The threshold for acoustic emission detection was set to 40 dB and the loading rate of the specimen was 0.167 kN/s.After setting-up the test apparatus, the processes of data acquisition are: first, during the four-point bending load, the crumb rubber concrete specimen generates internal cracks, thus forming stress waves (AE source); second, the AE sensors collect the propagated AE stress waves and convert them; third, the preamplifier amplifies the AE signals, thus ensuring the appropriate signal strength; fourth, the AE signals are stored and analyzed by the data acquisition system; and, finally, the AE characteristic parameters obtained are analyzed and studied, such as hit, amplitude, count, duration, etc.In addition, the multichannel monitoring system AE win was used as the software for AE signal acquisition and analysis, and the Matlab 10 Neural Network Toolbox was used for subsequent processing of the data.(1) K-means clustering: The main steps of the k-means clustering are: first, the number of clustering categories k (k cluster centers) is artificially determined, and then the similarity between each object and the cluster center are analyzed according to the distance function, and the sample points are classified into the most similar clusters.Next, the centroid of each cluster (the cluster center) is recalculated and this process is repeated until the centroid no longer changes.There are two phases of the k-means clustering algorithm [40]: Distribution phase: Use the nearest cluster center to assign each datum to the same cluster.
among them C i (t) is the cluster i, m i (t) is the center of cluster i, and A n represents the data points.(b) Update phase: Recalculate the cluster center and treat it as the average of the whole cluster.
Next, iterate through the above phases until the cluster no longer changes.Due to the initial randomly selected clustering center playing a vital role on the results of the entire calculation, k-means clustering cannot converge to the optimal solution to a large extent [43].
(2) FCM clustering The main idea of the fuzzy c-means (FCM) clustering follows.Assume a sample set and divide it into c fuzzy groups, and then find the cluster center of each group.Comparing with k-means clustering, this clustering method integrates the essence of fuzzy theory, thereby providing more flexible clustering results.Each object is not directly assigned to a specific cluster, but a weight is assigned to each object and each cluster.If the weight is larger, the object is closer to this cluster.This method attempts to find the minimum of the following objective functions [44]: Through comparison, the difference can be found between these functions and the function of k-means clustering: the fuzzy degree of clustering is determined by the fuzzy coefficient (∂) and the member value (γ ij ) together.The minimum value of the function can be obtained if [44]: There are three stages in the FCM clustering algorithm: (a) Determine cluster centers according to Equations ( 5) and ( 6).(b) Calculate the distance by Equation (4).(c) Update the objective function according to Equation (3).
The FCM clustering has better performance than the k-means clustering, but some limitations remain.For example, the determination of the initial cluster center and member values have a great impact on the clustering results [45], which may lead to the limitation of FCM clustering at local minimum values.
(3) GMM clustering GMM is an abbreviation for Gaussian mixture model (or mixed Gaussian model).According to many phenomena in nature, Gaussian (normal) distribution is followed; nevertheless, the reasons for affecting a distribution are multiple, and each influence factor may determine one.A single effect is a Gaussian distribution; a combination of multiple effects is a multiple Gaussian distribution.
Therefore, the Gaussian mixture model can also be understood as a large model with K single models.The Gaussian distribution is generally combined with the EM algorithm as its likelihood estimation algorithm.The probability model GMM can be expressed as [46]: where w k is the mixed weight function, A are the data, and the mean vector µ k and the covariance matrix ∑ k are represented by a Gaussian density function.There are three steps in the GMM clustering algorithm [46]: First, the algorithm estimates two quantities: the mean vector and the covariance matrix.It then calculates all the data and all the blending weight functions.(b) New parameters are calculated using membership and data.If N k = ∑ wik , the new blending weight can be expressed as: The new mean vector and covariance matrix can be expressed as: (c) Repeat the above two steps to meet the stopping criteria.
(4) SOM clustering The self-organizing mapping (SOM) neural network is a self-organizing and competitive network, and it is also a type of unsupervised learning model.SOM is unique in that it adopts a competitive learning approach to optimize the network.The weight vector of the "winning neuron" of the output layer and its neighboring neurons are calculated according to the following Equation [47]: where u represents the evaluation parameter of the neurons of A i , t represents the time, θ(u, v, t) represents the neighborhood function, r u and r v represent locations of neurons, u and v, and σ(t) are used to describe the width of the function θ(u, v, t).After this calculation, the configuration of the data can be obtained through the SOM clustering lattice, so as to identify and classify the internal data.
(5) Hierarchical model clustering The hierarchical clustering method refers to the hierarchical decomposition of the data set until a certain condition is met.The traditional hierarchical clustering algorithm mainly has two parts: agglomerative and divisive hierarchical clustering [48].

(a)
Agglomerative hierarchical clustering A bottom-up method which starts by determining each object as a single cluster, and these clusters are gradually merged according to certain criteria (the similarity measurement between the two clusters), and the clustering process repeats until all objects meet the number of clusters.(b) Divisive hierarchical clustering This part uses a strategy that is contrary to agglomerative hierarchical clustering: top-down, which first puts all objects into a cluster and then subdivides them into smaller clusters according to a certain metric until the distance between the two nearest clusters exceeds a certain threshold or reaches the number of clusters.
Although hierarchical clustering is relatively simple, it often encounters the problem of selecting merge points or split points.This choice is important, otherwise the process of merging or splitting will be irreversible, which means the processing completed cannot be undone, and the objects cannot be exchanged between clusters.If the decision to merge or split is not appropriate, a low quality clustering result may result.Moreover, hierarchical clustering algorithms do not have good scalability, and a large number of objects and clusters need to be examined and evaluated before deciding to merge or split.

Introduction to Density Peak Clustering Method
In this section, a new clustering method is applied to analyze AE data.The density peak clustering method is an improved algorithm.For the original algorithm, the determination of cluster centers can only be manually selected [49].The improved density peak clustering algorithm can achieve automatic clustering.The specific flow chart is shown in Figure 3.The following describes the specific steps of the density peak clustering method:

Introduction to Artificial Neural Network Method
An artificial neural network is a program that can effectively extract nonlinear data characteristics.It can provide computerized positioning and identify differential faults without requiring a large amount of understanding of sample behavior.It mainly has three parts: input, output, and hidden layer (Figure 4); the three layers are connected to form a network.First, the density of data points is calculated.The density of data point i is calculated as: where d c is the cutoff distance, and then calculate the distance δ from each data point to the data point with greater local density and the closest distance, as follows: Second, for purpose of automatically identifying the cluster center, the density peak S is used to measure the weight of each point as a candidate for the cluster center, as follows: Finally, there are several large values of S i , and data points are divided into several clusters based on the density method.

Introduction to Artificial Neural Network Method
An artificial neural network is a program that can effectively extract nonlinear data characteristics.It can provide computerized positioning and identify differential faults without requiring a large amount of understanding of sample behavior.It mainly has three parts: input, output, and hidden layer (Figure 4); the three layers are connected to form a network.Artificial neural networks typically use the gradient descent method [50][51][52] to perform error minimization efforts.The random weight is calculated by the given network connection method and the output initial value.The network performance can be improved through continuous loop iterations, and finally determine the most reliable system structure and the appropriate number of cycles [53].

Results and Discussion
In this section, the characteristic parameters of the AE signals collected from the crumb rubber concrete specimens in the four-point bending test were analyzed, mainly involving the following 10 parameters: rise time, count, energy, duration, amplitude, RMS, ASL, signal strength, center frequency, and peak frequency.The rise time and duration refer to the time interval between when the AE signal crosses the threshold and the maximum amplitude, and when the AE signal crosses the threshold and drops to the threshold, respectively.Count refers to the number of oscillations of the AE signal crossing the threshold, and energy refers to the area enclosed by the coordinates below the waveform.Additionally, center frequency and peak frequency refer to the information in the spectrum obtained from the FFT of the AE signal.Artificial neural networks typically use the gradient descent method [50][51][52] to perform error minimization efforts.The random weight is calculated by the given network connection method and the output initial value.The network performance can be improved through continuous loop iterations, and finally determine the most reliable system structure and the appropriate number of cycles [53].

Results and Discussion
In this section, the characteristic parameters of the AE signals collected from the crumb rubber concrete specimens in the four-point bending test were analyzed, mainly involving the following 10 parameters: rise time, count, energy, duration, amplitude, RMS, ASL, signal strength, center frequency, and peak frequency.The rise time and duration refer to the time interval between when the AE signal crosses the threshold and the maximum amplitude, and when the AE signal crosses the threshold and drops to the threshold, respectively.Count refers to the number of oscillations of the AE signal crossing the threshold, and energy refers to the area enclosed by the coordinates below the waveform.Additionally, center frequency and peak frequency refer to the information in the spectrum obtained from the FFT of the AE signal.

Clustering Results of Common Clustering Methods
The five common clustering methods described above have one aspect in common: they all need to be classified first and then the optimal number of clusters determined.Therefore, it is necessary to find the cluster number before using the five common clustering methods to cluster the crumb rubber concrete acoustic emission signals.To solve this problem, four evaluation indexes, Davies-Bouldin, Tou, Silhouette [54], and Dunn, were used in this paper.
Appl.Sci.2021, 11, 11476 9 of 24 When using the five common clustering methods and four clustering evaluation indexes in the clustering analysis, a total of six acoustic emission parameters, including rise time, count, energy, duration, amplitude, and peak frequency, were used for analysis [55,56].The clustering evaluation index results of these five commonly used clustering methods can be seen in Figure 5.According to the description of the four clustering evaluation indicators in this section, the smaller the Davies-Bouldin (DB) index, the better the classification effect.The Dunn (D) index, the Silhouette index, and the Tou index are numerical values indicating classification; the higher their value, the better the classification effect.According to Figure 5, the optimal clustering number corresponding to these five clustering methods is three.

Clustering Results of Common Clustering Methods
The five common clustering methods described above have one aspect in common: they all need to be classified first and then the optimal number of clusters determined.Therefore, it is necessary to find the cluster number before using the five common clustering methods to cluster the crumb rubber concrete acoustic emission signals.To solve this problem, four evaluation indexes, Davies-Bouldin, Tou, Silhouette [54], and Dunn, were used in this paper.
When using the five common clustering methods and four clustering evaluation indexes in the clustering analysis, a total of six acoustic emission parameters, including rise time, count, energy, duration, amplitude, and peak frequency, were used for analysis [55,56].The clustering evaluation index results of these five commonly used clustering methods can be seen in Figure 5.According to the description of the four clustering evaluation indicators in this section, the smaller the Davies-Bouldin (DB) index, the better the classification effect.The Dunn (D) index, the Silhouette index, and the Tou index are numerical values indicating classification; the higher their value, the better the classification effect.According to Figure 5, the optimal clustering number corresponding to these five clustering methods is three.In addition, according to the clustering results of the five commonly used clustering methods shown in Figure 6, it is found that some data points cannot be clearly distinguished, and the division of the three clusters is not particularly obvious.Therefore, it is crucial to find a better clustering method to analyze the acoustic emission data of crumb rubber concrete.In addition, according to the clustering results of the five commonly used clustering methods shown in Figure 6, it is found that some data points cannot be clearly distinguished, and the division of the three clusters is not particularly obvious.Therefore, it is crucial to find a better clustering method to analyze the acoustic emission data of crumb rubber concrete.

Clustering Results of the Density Peak Clustering Method
Based on the introduction of the density peak clustering method provided previously, this method was adopted to process the experimental data of crumb rubber concrete.Similarly, six acoustic emission parameters-rise time, count, energy, duration, amplitude, and peak frequency-were used as input data, and the peak density of each data point is shown in Figure 7.The 2210th data point, the 5176th data point and the 14,373rd data point have significantly higher density peaks than the other data points.Based on the principle of density peak clustering, these three points were used as the density center.

Clustering Results of the Density Peak Clustering Method
Based on the introduction of the density peak clustering method provided previously, this method was adopted to process the experimental data of crumb rubber concrete.Similarly, six acoustic emission parameters-rise time, count, energy, duration, amplitude, and peak frequency-were used as input data, and the peak density of each data point is shown in Figure 7.The 2210th data point, the 5176th data point and the 14,373rd data point have significantly higher density peaks than the other data points.Based on the principle of density peak clustering, these three points were used as the density center.The cluster center of the peak cluster, and all the remaining data points were classified separately.The cluster center of the peak cluster, and all the remaining data points were classified separately.After clustering by using the density peak method, the method described in this section also uses the average frequency and peak frequency to evaluate the clustering results.Based on Figure 8, the data points are well divided into three clusters, which can be clearly distinguished, and the average frequency between different clusters can be clearly discerned.The clustering performance is obviously better than the clustering performance of the five commonly used clustering methods noted above.Therefore, the density peak clustering method is the best clustering method for analyzing the acoustic emission data of crumb rubber concrete.This method can be used to perform deeper damage pattern recognition of crumb rubber concrete.After clustering by using the density peak method, the method described in this section also uses the average frequency and peak frequency to evaluate the clustering results.Based on Figure 8, the data points are well divided into three clusters, which can be clearly distinguished, and the average frequency between different clusters can be clearly discerned.The clustering performance is obviously better than the clustering performance of the five commonly used clustering methods noted above.Therefore, the density peak clustering method is the best clustering method for analyzing the acoustic emission data of crumb rubber concrete.This method can be used to perform deeper damage pattern recognition of crumb rubber concrete.The cluster center of the peak cluster, and all the remaining data points were classified separately.After clustering by using the density peak method, the method described in this section also uses the average frequency and peak frequency to evaluate the clustering results.Based on Figure 8, the data points are well divided into three clusters, which can be clearly distinguished, and the average frequency between different clusters can be clearly discerned.The clustering performance is obviously better than the clustering performance of the five commonly used clustering methods noted above.Therefore, the density peak clustering method is the best clustering method for analyzing the acoustic emission data of crumb rubber concrete.This method can be used to perform deeper damage pattern recognition of crumb rubber concrete.According to previous studies, it can be shown that the peak frequency of AE signals generated by the structural failure of concrete changes with the change of the concrete damage mode.Thus, various damage patterns of concrete can be identified and characterized through the peak frequency.The statistical analysis and comparison of the peak frequencies of different events in two kinds of concrete (ordinary concrete and crumb rubber concrete) were carried out (Figures 9 and 10).For ordinary concrete, the peak frequency distribution is mainly concentrated in one band, about 35~50 kHz.Hence, Figure 9 shows that the peak frequency of the main frequency band corresponding to the destruction of the concrete matrix, which is about 35~50 kHz.For crumb rubber concrete, due to the incorporation of rubber, another frequency band with a peak frequency of 85~125 kHz is added to Figure 10, which also represents another new failure mode in addition to the failure of the concrete matrix.
Appl.Sci.2021, 11, x FOR PEER REVIEW 13 of 25 ber concrete) were carried out (Figures 9 and 10).For ordinary concrete, the peak frequency distribution is mainly concentrated in one band, about 35~50 kHz.Hence, Figure 9 shows that the peak frequency of the main frequency band corresponding to the destruction of the concrete matrix, which is about 35~50 kHz.For crumb rubber concrete, due to the incorporation of rubber, another frequency band with a peak frequency of 85~125 kHz is added to Figure 10, which also represents another new failure mode in addition to the failure of the concrete matrix.Based on the results of density peak clustering, the optimal clustering number of crumb rubber concrete is three.Therefore, the acoustic emission signals of crumb rubber concrete were classified into three clusters.According to Figure 11, for the peak frequency of the crumb rubber concrete AE signals, the distribution range of cluster 1 is 10 to 66 kHz, and its main frequency band is 35 to 50 kHz.According to the previous analysis, the damage pattern represented by the cluster 1 signal is the cracking of the concrete matrix.The peak frequency of cluster 2 ranges from 5 to 170 kHz, of which the two main frequency bands are 35 to 50 kHz and 85 to 125 kHz, respectively, which shows that the signal of cluster 2 represents two failure modes: concrete matrix cracking accompanied by rubber particle fracture.The peak frequency distribution of cluster 3 ranges from 20 to 150 kHz, and it has only a main frequency band of 85 to 125 kHz, indicating that the signal of cluster 3 represents a new damage mode after rubber incorporation.ber concrete) were carried out (Figures 9 and 10).For ordinary concrete, the peak frequency distribution is mainly concentrated in one band, about 35~50 kHz.Hence, Figure 9 shows that the peak frequency of the main frequency band corresponding to the destruction of the concrete matrix, which is about 35~50 kHz.For crumb rubber concrete, due to the incorporation of rubber, another frequency band with a peak frequency of 85~125 kHz is added to Figure 10, which also represents another new failure mode in addition to the failure of the concrete matrix.Based on the results of density peak clustering, the optimal clustering number of crumb rubber concrete is three.Therefore, the acoustic emission signals of crumb rubber concrete were classified into three clusters.According to Figure 11, for the peak frequency of the crumb rubber concrete AE signals, the distribution range of cluster 1 is 10 to 66 kHz, and its main frequency band is 35 to 50 kHz.According to the previous analysis, the damage pattern represented by the cluster 1 signal is the cracking of the concrete matrix.The peak frequency of cluster 2 ranges from 5 to 170 kHz, of which the two main frequency bands are 35 to 50 kHz and 85 to 125 kHz, respectively, which shows that the signal of cluster 2 represents two failure modes: concrete matrix cracking accompanied by rubber particle fracture.The peak frequency distribution of cluster 3 ranges from 20 to 150 kHz, and it has only a main frequency band of 85 to 125 kHz, indicating that the signal of cluster 3 represents a new damage mode after rubber incorporation.Based on the results of density peak clustering, the optimal clustering number of crumb rubber concrete is three.Therefore, the acoustic emission signals of crumb rubber concrete were classified into three clusters.According to Figure 11, for the peak frequency of the crumb rubber concrete AE signals, the distribution range of cluster 1 is 10 to 66 kHz, and its main frequency band is 35 to 50 kHz.According to the previous analysis, the damage pattern represented by the cluster 1 signal is the cracking of the concrete matrix.The peak frequency of cluster 2 ranges from 5 to 170 kHz, of which the two main frequency bands are 35 to 50 kHz and 85 to 125 kHz, respectively, which shows that the signal of cluster 2 represents two failure modes: concrete matrix cracking accompanied by rubber particle fracture.The peak frequency distribution of cluster 3 ranges from 20 to 150 kHz, and it has only a main frequency band of 85 to 125 kHz, indicating that the signal of cluster 3 represents a new damage mode after rubber incorporation.In addition, the cumulative energy of AE parameters and the length of the propagating cracks generated when concrete is destroyed have also been studied, and scholars have proved that there is a certain correlation between the two.According to Figure 12, the cumulative energy of the signal of cluster 2 over time is much higher than that of the signals of cluster 1 and cluster 3, which indicates that the signal of cluster 2 corresponds to the larger damage inside the concrete matrix, also proving that the signal of cluster 2 includes not only the cracking of the concrete matrix, but also the new failure mode caused by rubber.At the same time, it should be noted that for concrete structures, the acoustic waves generated when the internal fractures are continuously attenuated during the propagation process, and considering that the distance between each acoustic emission sensor and the damage location is different, there is a certain error between the peak frequency received by sensors and the peak frequency at the original location of the concrete damage, which cannot be ignored.Thus, for different failure modes, it is crucial to find another method to more finely determine their corresponding peak frequencies.In addition, the cumulative energy of AE parameters and the length of the propagating cracks generated when concrete is destroyed have also been studied, and scholars have proved that there is a certain correlation between the two.According to Figure 12, the cumulative energy of the signal of cluster 2 over time is much higher than that of the signals of cluster 1 and cluster 3, which indicates that the signal of cluster 2 corresponds to the larger damage inside the concrete matrix, also proving that the signal of cluster 2 includes not only the cracking of the concrete matrix, but also the new failure mode caused by rubber.At the same time, it should be noted that for concrete structures, the acoustic waves generated when the internal fractures are continuously attenuated during the propagation process, and considering that the distance between each acoustic emission sensor and the damage location is different, there is a certain error between the peak frequency received by sensors and the peak frequency at the original location of the concrete damage, which cannot be ignored.Thus, for different failure modes, it is crucial to find another method to more finely determine their corresponding peak frequencies.
. Fast Fourier transform (FFT) is a calculation method that transforms the time domain into the frequency domain.Due to its high efficiency and speed, this method is commonly used to process acoustic emission waveform signals.For crumb rubber concrete, the peak frequencies of the cluster 1 and the cluster 3 signals have two identical main peaks: 40 and 100 kHz.Among them, the peak frequency 40 kHz of cluster 1 corresponds to a larger amplitude, and the peak frequency 100 kHz of cluster 3 corresponds to a larger amplitude (Figure 13).This indicates that the failure mode reflected by the cluster 1 signal and the cluster 3 signal corresponds to these two main peaks, which are mentioned above: the damage of the concrete matrix and the new failure mode caused by the incorporation of rubber.Meanwhile, the cluster 1 signal and cluster 3 signal are based on the failure mode of the concrete matrix (40 kHz) and the new failure mode (100 kHz) caused by rubber incorporation, respectively.Therefore, for different damage modes, the FFT method can be used to make their corresponding peak frequencies more accurate.Through FFT calculation, it is found that for rubber concrete, the cluster 1 signal is mainly the destruction of the concrete matrix, which also includes the new failure mode caused by rubber fracture.Meanwhile, the cluster 3 signal is mainly caused by the rubber fracture but also includes the damage of the concrete matrix.According to Table 3, the average value of the rise time, count, energy, duration, and amplitude of the cluster 3 signal (mainly caused by a new failure mode caused by rubber) is higher than the cluster 1 signal (mainly caused by the destruction of the concrete matrix), respectively.Thus, the intensity of the cluster 3 signal is slightly greater than the cluster 1 signal.Fast Fourier transform (FFT) is a calculation method that transforms the time domain into the frequency domain.Due to its high efficiency and speed, this method is commonly used to process acoustic emission waveform signals.For crumb rubber concrete, the peak frequencies of the cluster 1 and the cluster 3 signals have two identical main peaks: 40 and 100 kHz.Among them, the peak frequency 40 kHz of cluster 1 corresponds to a larger amplitude, and the peak frequency 100 kHz of cluster 3 corresponds to a larger amplitude (Figure 13).This indicates that the failure mode reflected by the cluster 1 signal and the cluster 3 signal corresponds to these two main peaks, which are mentioned above: the damage of the concrete matrix and the new failure mode caused by the incorporation of rubber.Meanwhile, the cluster 1 signal and cluster 3 signal are based on the failure mode of the concrete matrix (40 kHz) and the new failure mode (100 kHz) caused by rubber incorporation, respectively.Therefore, for different damage modes, the FFT method can be used to make their corresponding peak frequencies more accurate.Through FFT calculation, it is found that for rubber concrete, the cluster 1 signal is mainly the destruction of the concrete matrix, which also includes the new failure mode caused by rubber fracture.Meanwhile, the cluster 3 signal is mainly caused by the rubber fracture but also includes the damage of the concrete matrix.Fast Fourier transform (FFT) is a calculation method that transforms the time do into the frequency domain.Due to its high efficiency and speed, this method is comm used to process acoustic emission waveform signals.For crumb rubber concrete, the frequencies of the cluster 1 and the cluster 3 signals have two identical main peaks: 4 100 kHz.Among them, the peak frequency 40 kHz of cluster 1 corresponds to a l amplitude, and the peak frequency 100 kHz of cluster 3 corresponds to a larger ampl (Figure 13).This indicates that the failure mode reflected by the cluster 1 signal an cluster 3 signal corresponds to these two main peaks, which are mentioned above damage of the concrete matrix and the new failure mode caused by the incorporati rubber.Meanwhile, the cluster 1 signal and cluster 3 signal are based on the failure m of the concrete matrix (40 kHz) and the new failure mode (100 kHz) caused by ru incorporation, respectively.Therefore, for different damage modes, the FFT method be used to make their corresponding peak frequencies more accurate.Through FFT c lation, it is found that for rubber concrete, the cluster 1 signal is mainly the destructi the concrete matrix, which also includes the new failure mode caused by rubber frac Meanwhile, the cluster 3 signal is mainly caused by the rubber fracture but also inc the damage of the concrete matrix.According to Table 3, the average value of the rise time, count, energy, duration amplitude of the cluster 3 signal (mainly caused by a new failure mode caused by ru is higher than the cluster 1 signal (mainly caused by the destruction of the concrete trix), respectively.Thus, the intensity of the cluster 3 signal is slightly greater tha cluster 1 signal.According to Table 3, the average value of the rise time, count, energy, duration, and amplitude of the cluster 3 signal (mainly caused by a new failure mode caused by rubber) is higher than the cluster 1 signal (mainly caused by the destruction of the concrete matrix), respectively.Thus, the intensity of the cluster 3 signal is slightly greater than the cluster 1 signal.The average frequency (AF) refers to the ratio of the ringing counts to the duration time.Many scholars have found that the average frequency of acoustic emission signals corresponding to different damage modes is different.For example, the average frequency of acoustic emission signals produced by tensile damage is higher, while the average frequency of acoustic emission signals produced by shear damage is lower [52].The calculation Equation of the AF follows: AF = AE ringdown counts/the duration time (18) According to Table 3, the average frequency of cluster 1 and 3 signals are quite different.By comparing the average frequencies of these two clusters of signals, the different damage modes in crumb rubber concrete are further characterized.Figure 14 shows that the average frequency of cluster 1 (mainly caused by the destruction of the concrete matrix) is 1000 kHz, and the cluster 3 signal caused by the new failure mode has an average frequency band of 0 to 100 kHz.Hence, the main failure mode of the cluster 1 signal is the tensile failure mode, and the main failure mode of the cluster 3 signal is the shear failure mode.The average frequency (AF) refers to the ratio of the ringing counts to the duration time.Many scholars have found that the average frequency of acoustic emission signals corresponding to different damage modes is different.For example, the average frequency of acoustic emission signals produced by tensile damage is higher, while the average frequency of acoustic emission signals produced by shear damage is lower [52].The calculation Equation of the AF follows: AF AE ringdown counts the duration time  According to Table 3, the average frequency of cluster 1 and 3 signals are quite different.By comparing the average frequencies of these two clusters of signals, the different damage modes in crumb rubber concrete are further characterized.Figure 14 shows that the average frequency of cluster 1 (mainly caused by the destruction of the concrete matrix) is 1000 kHz, and the cluster 3 signal caused by the new failure mode has an average frequency band of 0 to 100 kHz.Hence, the main failure mode of the cluster 1 signal is the tensile failure mode, and the main failure mode of the cluster 3 signal is the shear failure mode.

Prediction of Ultimate Damage Load
In the previous section, the damage patterns of crumb rubber concrete were successfully identified using the AE technique combined with the selected optimal clustering method.At the same time, it is important to find a simple and effective method for predicting the ultimate damage load for crumb rubber concrete in order to ensure the safety of this material for a wide range of applications in practical engineering.On the other hand, the combination of damage pattern recognition and ultimate load prediction is more beneficial to characterize the damage mechanism of crumb rubber concrete.Therefore, in this section, a prediction model of ultimate load was proposed from the function of the AE technique with prediction of crack propagation combined with the artificial neural network.

Prediction of Ultimate Damage Load
In the previous section, the damage patterns of crumb rubber concrete were successfully identified using the AE technique combined with the selected optimal clustering method.At the same time, it is important to find a simple and effective method for predicting the ultimate damage load for crumb rubber concrete in order to ensure the safety of this material for a wide range of applications in practical engineering.On the other hand, the combination of damage pattern recognition and ultimate load prediction is more beneficial to characterize the damage mechanism of crumb rubber concrete.Therefore, in this section, a prediction model of ultimate load was proposed from the function of the AE technique with prediction of crack propagation combined with the artificial neural network.

Prediction of Artificial Neural Network of Acoustic Emission Single-Parameter Data
In this four-point bending test, nine crumb rubber concrete specimens were used: three specimens in which the rubber content was 5%, 10%, and 15% for each.The specific dimensions of the experimental specimens can be seen in the experimental part.In predicting the damage ultimate load of crumb rubber concrete using the artificial neural network, the AE signals and load data of one specimen from each of the three sets of specimens with different rubber contents were extracted as the prediction set, and the AE signals and load data of the other two specimens from the same set were used as the training set to predict the damage ultimate load of the specimens.This operation was repeated three times with different specimens, and the predicted and true values of the three sets of results were averaged and analyzed for errors.
In this section, 10 acoustic emission parameters-rise time (1), count (2), energy (3), duration (4), amplitude (5), RMS (6), ASL (7), signal strength (8), center frequency (9), and peak frequency (10)-predict the damage ultimate load of three different crumb rubber concretes, and compare the prediction results with their true values to calculate the prediction error.The Matlab 10 Neural Network Toolbox was used to train the networks.The Levenberge-Marquardt algorithm was used for back propagation during training.At the same time, the learning rate was set to 0.9 in the calculation process.
In the process of obtaining the prediction results shown in Tables 4-13, the input data are a matrix composed of the single acoustic emission parameter column vectors of the six crumb rubber concrete specimens, and the hidden layer is composed of nine neurons, and the output end is the predicted value of the ultimate load of the crumb rubber concrete.The prediction result of the artificial neural network with a single parameter of rise time follows: when the number of iterations reaches 80, the average error (MSE) reaches the minimum value (Figure 15a); that is, the optimal number of training iterations is 80 times.At the same time, the best iteration numbers of the other nine acoustic emission parameters can also be determined from Figure 15b-j.Table 4 shows that the maximum prediction error of the specimens with three different rubber contents is -2.39%.It is concluded that the single-parameter neural network with rise time predicts that the damage ultimate load error of crumb rubber concrete is small, and the results using different acoustic emission parameters is different from what is shown in Tables 4-13.   Figure 16 shows the absolute values of the prediction errors of the 10 AE singleparameter neural networks.The prediction errors of the two single-parameter neural networks with rise time (1) and amplitude ( 5) are relatively small.The predictions with center frequency (9) and peak frequency (10) are moderate, all within 8%.The remaining parameters of the neural network prediction results are poor.Therefore, rise time is the most suitable acoustic emission parameter for single-parameter neural network prediction of crumb rubber concrete failure load.Figure 16 shows the absolute values of the prediction errors of the 10 AE single-parameter neural networks.The prediction errors of the two single-parameter neural networks with rise time (1) and amplitude ( 5) are relatively small.The predictions with center frequency (9) and peak frequency (10) are moderate, all within 8%.The remaining parameters of the neural network prediction results are poor.Therefore, rise time is the most suitable acoustic emission parameter for single-parameter neural network prediction of crumb rubber concrete failure load.

Prediction of Artificial Neural Network of Acoustic Emission Two-Parameter Data
In this section, using the rise time (1) and amplitude (5) acoustic emission parameters, the damage ultimate loads of three different crumb rubber concretes were predicted by the acoustic emission two-parameter data artificial neural network.The prediction result was calculated by comparing the predicted result with its true value.The training tools and algorithms are the same as the prediction using AE single parameter.
In the process of deriving the prediction results shown in Table 14, the input data consist of a two-parameter matrix composed of the rise time and amplitude column vectors of six specimens, and the hidden layer consisting of nine nerves.The composition of the element is the predicted value of the ultimate load of the crumb rubber concrete.The artificial neural network prediction result with rise time and amplitude double parameters provides an optimal number of 39 training iterations (Figure 17).At the same time, the maximum prediction error of the specimens with three different rubber contents is 2.11%, shown in Table 14.It can be concluded that the rise-time-amplitude two-parameter neural network is useful to predict the damage ultimate load of crumb rubber concrete.

Prediction of Artificial Neural Network of Acoustic Emission Two-Parameter Data
In this section, using the rise time (1) and amplitude (5) acoustic emission parameters, the damage ultimate loads of three different crumb rubber concretes were predicted by the acoustic emission two-parameter data artificial neural network.The prediction result was calculated by comparing the predicted result with its true value.The training tools and algorithms are the same as the prediction using AE single parameter.
In the process of deriving the prediction results shown in Table 14, the input data consist of a two-parameter matrix composed of the rise time and amplitude column vectors of six specimens, and the hidden layer consisting of nine nerves.The composition of the element is the predicted value of the ultimate load of the crumb rubber concrete.The artificial neural network prediction result with rise time and amplitude double parameters provides an optimal number of 39 training iterations (Figure 17).At the same time, the maximum prediction error of the specimens with three different rubber contents is 2.11%, shown in Table 14.It can be concluded that the rise-time-amplitude two-parameter neural network is useful to predict the damage ultimate load of crumb rubber concrete.After using these two acoustic emission parameters to combine the two-parameter neural network prediction of the ultimate load of crumb rubber concrete, the error of prediction with this set of two parameters is the smallest and less than the respective prediction errors of these two parameters.Thus, the rise-time-amplitude two-parameter combination is the most suitable for the two-parameter neural network to determine the acoustic emission parameters of the crumb rubber concrete failure load.Furthermore, the prediction error of all acoustic emission two-parameter neural network models is significantly smaller than that of the corresponding acoustic emission single-parameter neural network parameter model.Figure 17 shows the two-parameter optimal combination rise-time-amplitude (RA) neural network training performance.After using these two acoustic emission parameters to combine the two-parameter neural network prediction of the ultimate load of crumb rubber concrete, the error of prediction with this set of two parameters is the smallest and less than the respective prediction errors of these two parameters.Thus, the rise-time-amplitude two-parameter combination is the most suitable for the two-parameter neural network to determine the acoustic emission parameters of the crumb rubber concrete failure load.Furthermore, the prediction error of all acoustic emission two-parameter neural network models is significantly smaller than that of the corresponding acoustic emission single-parameter neural network parameter model.Figure 17 shows the two-parameter optimal combination rise-time-amplitude (RA) neural network training performance.

Prediction of Artificial Neural Network of Acoustic Emission Multiparameter Data
In this section, five acoustic emission parameters were used together, including rise time (1), duration (4), amplitude (5), center frequency (9), and peak frequency (10) as the input data for three different rubber contents.The ultimate load of the crumb rubber concrete was predicted by the acoustic wave multiparameter data artificial neural network, and the prediction result was compared with its real value to calculate the prediction error.The training tools and algorithms are the same as the prediction using AE single parameter.
In the process of reporting the prediction results in Table 15, the input data are the rise time (1), duration (4), amplitude (5), center frequency (9), and peak frequency (10) of the six crumb rubber concrete specimens.Therefore, it is a multiparameter matrix composed of five column vectors.The hidden layer is composed of nine neurons, and the output end is the predicted value of the ultimate load.As for the prediction results of the artificial neural network using AE multiparameters, the optimal number of training iterations is 52 (Figure 18).At the same time, Table 15 shows that the maximum prediction error of the specimens with three different rubber contents is only 0.89%.Therefore, the multiparameter neural network is useful to predict the ultimate load of crumb rubber concrete.

Prediction of Artificial Neural Network of Acoustic Emission Multiparameter Data
In this section, five acoustic emission parameters were used together, including rise time (1), duration (4), amplitude (5), center frequency (9), and peak frequency (10) as the input data for three different rubber contents.The ultimate load of the crumb rubber concrete was predicted by the acoustic wave multiparameter data artificial neural network, and the prediction result was compared with its real value to calculate the prediction error.
The training tools and algorithms are the same as the prediction using AE single parameter.
In the process of reporting the prediction results in Table 15, the input data are the rise time (1), duration (4), amplitude (5), center frequency (9), and peak frequency (10) of the six crumb rubber concrete specimens.Therefore, it is a multiparameter matrix composed of five column vectors.The hidden layer is composed of nine neurons, and the output end is the predicted value of the ultimate load.As for the prediction results of the artificial neural network using AE multiparameters, the optimal number of training iterations is 52 (Figure 18).At the same time, Table 15 shows that the maximum prediction error of the specimens with three different rubber contents is only 0.89%.Therefore, the multiparameter neural network is useful to predict the ultimate load of crumb rubber concrete.Ten acoustic emission single parameters, rise-time-amplitude (RA) two-parameter and multiparameter neural network prediction error absolute values are also shown in Figure 19, which can be summarized by comparison: the acoustic emission multiparameter neural network prediction results are the most accurate and the prediction error is no more than 1%.The rise-time-amplitude (RA) two-parameter neural network prediction results are also small.Ten acoustic emission single parameters, rise-time-amplitude (RA) two-parameter and multiparameter neural network prediction error absolute values are also shown in Figure 19, which can be summarized by comparison: the acoustic emission multiparameter neural network prediction results are the most accurate and the prediction error is no more than 1%.The rise-time-amplitude (RA) two-parameter neural network prediction results are also small.

Conclusions
(1) In this work, six clustering algorithms were compared to find the best one to analyze the acoustic emission data of crumb rubber concrete.Based on the results and analysis provided, the density peak clustering algorithm is considered to be an optimal algorithm based on the reliability of its computational results.(2) Therefore, the density peak clustering method was adopted to cluster the AE signals of the specimen and to investigate the variation of different damage modes during   Ten acoustic emission single parameters, rise-time-amplitude (RA) two-parameter and multiparameter neural network prediction error absolute values are also shown in Figure 19, which can be summarized by comparison: the acoustic emission multiparameter neural network prediction results are the most accurate and the prediction error is no more than 1%.The rise-time-amplitude (RA) two-parameter neural network prediction results are also small.

Conclusions
(1) In this work, six clustering algorithms were compared to find the best one to analyze the acoustic emission data of crumb rubber concrete.Based on the results and analysis provided, the density peak clustering algorithm is considered to be an optimal algorithm based on the reliability of its computational results.(2) Therefore, the density peak clustering method was adopted to cluster the AE signals of the specimen and to investigate the variation of different damage modes during

Conclusions
(1) In this work, six clustering algorithms were compared to find the best one to analyze the acoustic emission data of crumb rubber concrete.Based on the results and analysis provided, the density peak clustering algorithm is considered to be an optimal algorithm based on the reliability of its computational results.(2) Therefore, the density peak clustering method was adopted to cluster the AE signals of the specimen and to investigate the variation of different damage modes during loading.The damage modes of concrete matrix fracture and rubber fracture correspond to peak frequencies of 40 and 100 kHz, respectively.Based on the results provided, combining AE technology with density peak clustering, the damage pattern identification of crumb rubber concrete could be carried out at a deeper level.(3) For the prediction of the ultimate load of crumb rubber concrete, ten different acoustic emission parameter data were selected to establish an acoustic emission singleparameter artificial neural network prediction model.The prediction results of these ten single-parameter models were compared with the real values, and the prediction errors were calculated separately.Furthermore, five optimal acoustic emission parameter data were selected as the artificial neural network data input end to establish the acoustic emission multiparameter artificial neural network prediction model so as to predict the ultimate load.According to the results, the prediction error is only 0.89%.(4) These methods and conclusions adopted in this study are general and practical for the application and monitoring of crumb rubber concrete in practical engineering.However, some imperfections remain, and in the future, it will be necessary to estab-

Figure 1 .
Figure 1.Configuration of the four-point bending specimen and sensors.

Figure 1 .
Figure 1.Configuration of the four-point bending specimen and sensors.

Figure 2 .
Figure 2. (a) AE experimental equipment and (b) labeled image of the four-point bending test.

Figure 2 .
Figure 2. (a) AE experimental equipment and (b) labeled image of the four-point bending test.

2. 2 .
Introduction to Common Clustering Methods This section introduces five common clustering methods including k-means, fuzzy c-means (FCM), Gaussian mixture model (GMM), self-organizing mapping (SOM), and hierarchical model clustering to analyze acoustic emission signals.The following five clustering methods are introduced separately: Appl.Sci.2021, 11, x FOR PEER REVIEW 8 of 25

Figure 3 .
Figure 3. Flow chart of density peak clustering.

Figure 3 .
Figure 3. Flow chart of density peak clustering.

2. 4 .
Introduction to Artificial Neural Network MethodAn artificial neural network is a program that can effectively extract nonlinear data characteristics.It can provide computerized positioning and identify differential faults without requiring a large amount of understanding of sample behavior.It mainly has three parts: input, output, and hidden layer (Figure4); the three layers are connected to form a network.

Figure 3 .
Figure 3. Flow chart of density peak clustering.

Figure 7 .
Figure 7. Peak density of sample points.

Figure 8 .
Figure 8. Density peak clustering results.3.1.3.Damage Pattern Recognition of Crumb Rubber Concrete by the Density Peak Clustering Method According to previous studies, it can be shown that the peak frequency of AE signals generated by the structural failure of concrete changes with the change of the concrete damage mode.Thus, various damage patterns of concrete can be identified and characterized through the peak frequency.The statistical analysis and comparison of the peak frequencies of different events in two kinds of concrete (ordinary concrete and crumb rub-

Figure 7 .
Figure 7. Peak density of sample points.

Figure 7 .
Figure 7. Peak density of sample points.

Figure 8 .
Figure 8. Density peak clustering results.3.1.3.Damage Pattern Recognition of Crumb Rubber Concrete by the Density Peak Clustering Method According to previous studies, it can be shown that the peak frequency of AE signals generated by the structural failure of concrete changes with the change of the concrete damage mode.Thus, various damage patterns of concrete can be identified and characterized through the peak frequency.The statistical analysis and comparison of the peak frequencies of different events in two kinds of concrete (ordinary concrete and crumb rub-

Figure 8 .
Figure 8. Density peak clustering results.3.1.3.Damage Pattern Recognition of Crumb Rubber Concrete by the Density Peak Clustering Method

Figure 9 .
Figure 9. Peak frequency distribution of ordinary concrete.

Figure 10 .
Figure 10.Peak frequency distribution of crumb rubber concrete.

Figure 9 .
Figure 9. Peak frequency distribution of ordinary concrete.

Figure 9 .
Figure 9. Peak frequency distribution of ordinary concrete.

Figure 10 .
Figure 10.Peak frequency distribution of crumb rubber concrete.

Figure 10 .
Figure 10.Peak frequency distribution of crumb rubber concrete.

Figure 12 .
Figure 12.Cumulative energy of three clusters.

Figure 16 .
Figure 16.Prediction error of each single-parameter neural network.

Figure 16 .
Figure 16.Prediction error of each single-parameter neural network.

Table 1 .
Specific dimensions of the specimen and sensors.

Table 1 .
Specific dimensions of the specimen and sensors.

Table 2 .
Information on test apparatus.

Table 2 .
Information on test apparatus.

Table 3 .
Average AE parameters of cluster 1 and cluster 3.

Table 3 .
Average AE parameters of cluster 1 and cluster 3.
Rubber Content Actual Value (kN) Predictive Value (kN)