The Use of the Acoustic Emission Method to Identify Crack Growth in 40CrMo Steel

The article presents the application of the acoustic emission (AE) technique for detecting crack initiation and examining the crack growth process in steel used in engineering structures. The tests were carried out on 40CrMo steel specimens with a single edge notch in bending (SENB). In the tests crack opening displacement, force parameter, and potential drop signal were measured. The fracture mechanism under loading was classified as brittle. Accurate AE investigations of the cracking process and SEM observations of the fracture surfaces helped to determine that the cracking process is a more complex phenomenon than the commonly understood brittle fracture. The AE signals showed that the frequency range in the initial stage of crack development and in the further crack growth stages vary. Based on the analysis of parameters and frequencies of AE signals, it was found that the process of apparently brittle fracture begins and ends according to the mechanisms characteristic of ductile crack growth. The work focuses on the comparison of selected parameters of AE signals recorded in the pre-initiation phase and during the growth of brittle fracture cracking.


Introduction
Diagnosing and monitoring the condition of structures is an important and very topical area in construction. Aging infrastructure and increasing service load of engineering structures are the main drivers of fast-progressing research in the new interdisciplinary field of knowledge called structural health monitoring (SHM), which is closely related to the durability and safe operation over the life of the structural elements.
The load is a typical random load that is difficult to model with currently known fatigue calculation procedures, and can be modelled only for steady loads with low-level amplitudes. As a result, the fatigue life calculation models for steel structures fail to provide accurate information about the risk of fatigue failure. Where calculation methods provide only limited information, the existing structures are subjected to tests. Fatigue damage testing is carried out during inspection. Cracks, which occur for a number of reasons, are detected and, in the case of fatigue cracks, the rate of their growth is assessed.
Depending on the operating conditions, the loading and materials used, the failure takes place according to a ductile or brittle fracture mechanism.
Ductile cracking usually occurs in steels and is preceded by significant plastic deformation that manifests itself long before the failure of the element, thus allowing preventive measures to be undertaken. The brittle fracture mechanism is much more dangerous.
Brittle cracking proceeds without visible deformation of the element and occurs almost immediately. The occurrence of brittle fracture can be expected mainly in high-strength structural steels or in low-temperature operation and overloading of the element with simultaneous increase in the load speed. Welded joints are also exposed to brittle fracture due to their inhomogeneous microstructure, welding defects and residual stresses. The presence of these factors together with the impact of fatigue loads and corrosive environment leads to the degradation of the material and development of micro-cracks, resulting in the brittle failure of steel.
Leading research centers worldwide have devoted a great deal of attention to the development of methods for assessing the condition of structural elements and preventing their failure. As a result, many methods for evaluation of component durability [1][2][3][4][5] are currently in wide use. These methods are based on aversion to change assuming prior knowledge of the element microstructure, load history, and in-service conditions. However, more accurate results will be obtained when as much information as possible can be collected on in-service parameters and material.
Hence, the modelling of bearing capacity and durability of a structure under real operating conditions requires that: • real operational loads are defined; • the material model is defined, in particular welded joints material that changes over time under the influence of operating conditions; • the load around the defect after its initiation and during development (redistribution of stresses) is determined; and • the interaction of various damage types is identified.
As determining the factors above is very difficult, if not impossible, most analyses are characterized a high degree of conservatism, which significantly reduces assessment accuracy. Appropriate solution should include developing an NDT-based monitoring system able to signal the structural safety risk.
Finding and determining the "destructive characterization" of all hot spots in large-sized structures using the NDT methods is nearly impossible. In steel structures, the volume of "hot spot" is of the order of cubic centimeters, with dimensions of the structure reaching tens and more meters. In our opinion, the solution to the problems of diagnosing engineering structures is the use of continuous, long-term monitoring using the passive NDT methods. One of them is the acoustic emission (AE) method [4][5][6][7][8][9][10][11][12][13][14].
The paper presents a proposal for the use of acoustic emission for the diagnosis of steel structures vulnerable to brittle fracture.
The choice of acoustic emission as the research method was determined by its advantages in relation to other non-destructive methods. These are: The main task of the measurement system consisting of an AE processor is to detect, record, filter and analyze the signals generated by AE sources.
Each destructive process is a source of acoustic emission. The source is described by the parameters of the recorded AE signal. The values of the parameters are used to classify the signals (and, thus, damage processes). The similarity of signals is used to attribute particular signals to the defects caused by specific damage processes, and then by applying statistical grouping methods to identify the existing defects.
For the statistical methods used in identification, an important issue is the optimal selection of recorded 13 AE parameters (counts, counts to peak, amplitude, RMS-root means square voltage, ASL-average signal level, energy, absolute energy, signal strength, rise time, duration, initiation frequency, average frequency, reverberation frequency). The parameters must be characterized by low mutual correlation. A set of diagnostic variables must describe the most important aspects of the studied phenomenon [15][16][17][18]. Techniques such as hierarchical and non-hierarchical clustering methods and Kohonen's neural networks are used to build the reference signals data base for identification of destructive processes in steel structures (IPDKS) [4,5].
This paper uses unsupervised pattern recognition methods to characterize different AE activities corresponding to different fracture mechanisms. A sequential feature selection method based on a k-means clustering algorithm is used to achieve high classification accuracy. Fatigue damage propagation represents the main failure factor. To study the contributions of different types of damage at progressive fatigue stages, a tool with the ability to detect the damage initiation and to monitor failure progress online is needed. Acoustic emission (AE) testing has become a recognized suitable and effective non-destructive technique to investigate and evaluate failure processes in different structural components. The main advantage of AE over other condition monitoring techniques is that detected AE signals can be used to characterize the different damage mechanisms. The approach of using parameter distribution has the advantage of real-time damage detection, but can also lead to false conclusions due to noise effects in the AE signals [18]. This is especially true for materials that are working under the fatigue conditions which are usually present in damage mechanism interactions. Therefore, AE techniques are needed to account for more intricate wave propagation features caused by the anisotropic nature of materials and to enable the identification of a large variety of failure modes. It is now possible to detect and capture very large numbers of AE signals, driving a trend for seeking computationally complex algorithms, such as pattern recognition, to determine the onset of significant AE. For a real structure it is not possible to provide a set of training patterns belonging to multiple damage mechanisms, which thus makes the use of unsupervised pattern recognition techniques more appropriate for these studies [18].
The first step of developing an unsupervised pattern recognition process is to classify signals into groups based on similarities. This process involves statistical effects, and the key point of successful feature selection to construct fine classification accuracy. A paper by Doan et al. [16] presented a feature selection method that introduced a sequential method based on the Gustafson-Kessel clustering algorithm. In these the method, the subset of features is selected by minimizing the Davies-Bouldin (DB) index, which is a metric for the evaluation of classification algorithms. The signals are classified into four groups by comparing their features and deciding upon their similarity.
The assignment of the clustering results to the fracture mechanisms is achieved by a detailed analysis of the physical meaning of the data [19][20][21][22]. This is the first time that the pattern recognition technique has been applied to a database acquired from such a complex structure in a fatigue testing environment [16]. The applied feature selection algorithm proves to be a powerful tool providing relevant clustering when used together with a k-means algorithm. When a crack occurs in the material, it results in a rapid release of energy, transmitting in the form of an elastic wave, namely acoustic emission (AE). The AE-based detection method has been intensively used also in non-destructive assessments of cracks in papers [23][24][25]. Additionally, Rabiei and Modarres revealed a log-linear relationship between the AE features and crack growth rate, and presented an end-to-end approach for structural health management [10]. Qu et al. presented a comparative study of the damage level diagnostics of gearbox tooth using AE and vibration measurements; the results indicated that vibration signals were easily affected by mechanical resonance, while the AE signals showed a more stable performance [26]. Zhang et al. have studied defect detection of rails using AE and wavelet transform at a high speed [27]. Li and al. created a template library of cracking sounds and designed a detection device using voiceprint recognition with an accuracy of 77% [28]. Hase et al. combined the Hurst exponent and the neural network to develop a crack detection algorithm of carbide anvils [27]. In the above methods, the AE sensors are usually attached to the surface of the monitoring object, what makes measurement difficult.
A more practical crack identification and detection method is still lacking. Aiming to improve recognition accuracy and generalization, in papers [29][30][31] a novel crack identification method based on acoustic emission and pattern recognition is proposed. In these methods, the sound pulses from cracks are firstly separated from the original signal by pre-processing. The high-dimensional features are reduced adaptively by using principal component analysis (PCA). The algorithm combines a k-nearest neighbor (kNN) classifier with a support vector machine (SVM) to refine the classification outcome. While debris monitoring does not require any electronics, it is simple to interpret and has excellent sensitivity to wear-related failure; this method is insufficient to non-benign cracks as no debris is produced. Acoustic emission (AE) and vibration signals have more quantitative results to detect the earliest stage of damage in rotating machinery. AE and vibration methods are based on recording transient signals in two different frequency spectrums. While vibration method is based on features that are extracted from time and frequency domain signals recorded by low frequency accelerometers in order to assess the changes in vibrational properties as related to the damage [29][30][31], the AE method is based on detecting propagating elastic waves released from active flaws. Once transient signals are collected, signal processing methods, such as wavelet decomposition [11], empirical mode decomposition [11], and multivariate pattern recognition [20], are applied. Typical parameters extracted from the transient signals are root mean square value, frequency domain characteristics, energy, spectral kurtosis, and peak-to-peak vibration level. Due to the difference in the frequency bandwidth, the AE method is more sensitive to microcracks as compared to the vibration method. Typical AE data acquisition approach is based on threshold: an AE signal is detected when the signal level is above threshold. As crack growth is a stochastic process, it is considered that while some data will be lost due to the idle time of data acquisition system between waveform recording intervals, crack information will be stochastically detected. However, it is important to identify how to analyze long-duration signals in order to reduce the influence of background noise from the extracted features. In studies, the AE signals, which accompanying of the fatigue crack growth is important obtained from the scaled laboratory experiments. Acoustic emission (AE) is a health monitoring approach which acts as a passive receiver to record internal activities in structures. This method is capable of continuous monitoring, which is not the case with most traditional methods. In addition, AE sensors are very sensitive and can capture signals due to micro-scale defect formations coming from the internal regions of structures rather than only those at the surface [18,23]. An unsupervised and supervised pattern recognition algorithm was employed to classify the AE signals. Different damage mechanisms for specimens during cracking were identified using reference database AE signals created from signals registered in the tests described below.

Materials and Methods
The tests discussed below were carried out using Zwick-100 testing machine (ZwickRoel Ulm, Germany ) on SENB (Single Edge Notch in Bending) specimens ( Figure 1) made of 40CrMo steel. The specimens with dimensions 12.5 mm × 25 mm × 110 mm were subjected to heating at 850 • C (15 min), quenching in oil, tempering at 250 • C (3 h) and cooling in oil. As a result, a material with a tempered martensite microstructure was obtained ( Figure 2a). The signals of loading, specimen deflection and acoustic emission (AE) were recorded during testing. SENB type specimens with a total fracture length, which includes a notch + previously fatigue crack, equal to about 0.5•a/W are recommended by the ASTM and PN-EN standards [32][33][34][35] and commonly used to determine fracture toughness characteristics-critical values of the stress intensity factor (SIF) KIC or the critical J-integral, JIC. If during the entire process of loading the SENB specimen in the net-section area before the crack tip there is a definite predominance of the linear-elastic nature of the stress and strain fields before the crack tip are described by SIF and the fracture toughness is characterized by the critical value of KIC. The condition allowing the use of SIF is to limit the size of the yielding before the crack tip: rp ≤ 0.01•a, where for a plane strain rp = (1/6π)•(KI/σy) 2 . If the plastic zone is larger, it means that the material in front of the fracture tip is elastic-plastic and for the description of mechanical fields an integral J must be used, and fracture toughness is represented by the critical integral J -JIC [36,37].
For determining material strength characteristics, cylindrical specimens 10 mm in diameter were prepared with the measuring section of 50 mm. For determining the critical value of fracture toughness, KIC, the specimens with dimensions 12.5 mm × 25 mm × 110 mm with one-sided notch and crack fatigue (SENB) with a total length of about 12.5 mm were made. Figure 2b shows an example plot of loading a cylindrical specimen. On the basis of the graphs, strength characteristics of 40CrMo steel were determined: Re = σy = 1475 MPa; Rm = σuts = 1800 MPa; E = 205 GPa. During the loading of the SENB specimen, it was determined that the crack growth process takes place according to the brittle fracture mechanism. All standard conditions were met. The obtained critical value of fracture toughness is a material characteristic and is equal to KIC = 39 MPa•m 1/2 . SENB type specimens with a total fracture length, which includes a notch + previously fatigue crack, equal to about 0.5·a/W are recommended by the ASTM and PN-EN standards [32][33][34][35] and commonly used to determine fracture toughness characteristics-critical values of the stress intensity factor (SIF) K IC or the critical J-integral, J IC . If during the entire process of loading the SENB specimen in the net-section area before the crack tip there is a definite predominance of the linear-elastic nature of the stress and strain fields before the crack tip are described by SIF and the fracture toughness is characterized by the critical value of K IC . The condition allowing the use of SIF is to limit the size of the yielding before the crack tip: r p ≤ 0.01·a, where for a plane strain r p = (1/6π)·(K I /σ y ) 2 . If the plastic zone is larger, it means that the material in front of the fracture tip is elastic-plastic and for the description of mechanical fields an integral J must be used, and fracture toughness is represented by the critical integral J -J IC [36,37].
For determining material strength characteristics, cylindrical specimens 10 mm in diameter were prepared with the measuring section of 50 mm. For determining the critical value of fracture toughness, K IC , the specimens with dimensions 12.5 mm × 25 mm × 110 mm with one-sided notch and crack fatigue (SENB) with a total length of about 12.5 mm were made. Figure 2b shows an example plot of loading a cylindrical specimen. On the basis of the graphs, strength characteristics of 40CrMo steel were determined: R e = σ y = 1475 MPa; R m = σ uts = 1800 MPa; E = 205 GPa. During the loading of the SENB specimen, it was determined that the crack growth process takes place according to the brittle fracture mechanism. All standard conditions were met. The obtained critical value of fracture toughness is a material characteristic and is equal to K IC = 39 MPa·m 1/2 . SENB type specimens with a total fracture length, which includes a notch + previously fatigue crack, equal to about 0.5•a/W are recommended by the ASTM and PN-EN standards [32][33][34][35] and commonly used to determine fracture toughness characteristics-critical values of the stress intensity factor (SIF) KIC or the critical J-integral, JIC. If during the entire process of loading the SENB specimen in the net-section area before the crack tip there is a definite predominance of the linear-elastic nature of the stress and strain fields before the crack tip are described by SIF and the fracture toughness is characterized by the critical value of KIC. The condition allowing the use of SIF is to limit the size of the yielding before the crack tip: rp ≤ 0.01•a, where for a plane strain rp = (1/6π)•(KI/σy) 2 . If the plastic zone is larger, it means that the material in front of the fracture tip is elastic-plastic and for the description of mechanical fields an integral J must be used, and fracture toughness is represented by the critical integral J -JIC [36,37].
For determining material strength characteristics, cylindrical specimens 10 mm in diameter were prepared with the measuring section of 50 mm. For determining the critical value of fracture toughness, KIC, the specimens with dimensions 12.5 mm × 25 mm × 110 mm with one-sided notch and crack fatigue (SENB) with a total length of about 12.5 mm were made. Figure 2b shows an example plot of loading a cylindrical specimen. On the basis of the graphs, strength characteristics of 40CrMo steel were determined: Re = σy = 1475 MPa; Rm = σuts = 1800 MPa; E = 205 GPa. During the loading of the SENB specimen, it was determined that the crack growth process takes place according to the brittle fracture mechanism. All standard conditions were met. The obtained critical value of fracture toughness is a material characteristic and is equal to KIC = 39 MPa•m 1/2 . The reference signal database in the IPDKS method for steel structures was developed using the Fuzzy k-means algorithm. The k-means algorithm belongs to a group of non-hierarchical clustering approaches. Its essence lies in the random choice of initial centres. In each iteration further approximations of the patterns are searched for and enumerated using the given methods. Depending on the assumptions made, the reference element may be one of the elements of the X population or belong to a certain universe U ⊇ X.
In metric spaces, a reference element can be calculated as an arithmetic mean to represent the centre of gravity of the cluster.
In general, the k-clustering algorithm input is a set of X objects and the expected number of clusters k, and the output is the division into subsets {C 1 , C 2 , . . . , C k }. Frequently, k-clustering algorithms belong to the category of optimization algorithms. Optimization algorithms assume that there is a loss function k: {x | X ⊆ S} → R + specified for each S subset. The aim is to find a branch due to minimizing the sum of losses described by Equation (1): In order to apply the iterative algorithm, the distance measure used should be determined in the grouping process. Instead of the median, you can use a point that is the resultant point for a given cluster (representation element) calculated, e.g., as a geometric or arithmetic mean. Then we deal with k-means (k-means) algorithms described by Equation (2): where: x i -the arithmetic mean (or geometric) of the cluster.
Doing so means that the centroids search for their correct positions using Equation (3) [4,5]: where P ω i x j is the conditional probability of belonging j-th element to the i-th group, b is the parameter that has to take values other than 1, x _ j is the j-th element. The probability function is normalized according to Equation (4): The probability of the membership of the element in each of the clusters P ω i x j is calculated from Equation (5): where d ij 2 = x j − µ i 2 is the distance of a data point x j from the center of the group µ i . The k-means algorithm consists of the following steps:

2.
Compute the distances between the data points and the cluster centroids.

5.
If: • there are no changes in µ i and P ω i x j -return µ _ 1 . . . µ _ c , • otherwise go back to Step 2.
When this algorithm is used, the number of clusters is pre-determined. However, the speed of computation compensates for this inconvenience. The reference signals database is created with 13 correlated AE parameters (rise time, durability, counts, counts to peak, energy, RMS, ASL, amplitude, average frequency, initiation frequency, reverberation frequency, absolute energy, signal strength). Reference signals obtained in this way make it possible to identify individual destructive processes in all tested structures. The database can be supplemented as per the diagram above whenever new measurement data are obtained.
The use of the grouping method is particularly important for the identification of fracture mechanisms as it allows for timely response to the possible averse effects.
The study used a 24-channel "µSamos" acoustic emission processor, 40 dB preamplifiers and two resonance sensors with flat characteristics in the 30-80 kHz range and two broadband sensor in the 100-1200 kHz range. The use of four sensors with flat characteristics was aimed at determining the lower and upper discrimination thresholds of the AE signal, while the broadband sensor had the task of controlling whether during the measurement there are no signals generated in other frequency bands.

Results and Discussion
During the loading of the SENB specimen, the crack growth process takes place according to the brittle fracture mechanism and critical value of fracture toughness is a material characteristic equal to K IC = 39 MPa·m 1/2 . At the critical moment, the length of the plastic zone according to the Irwin model for plane strain equals: 2r p = 2(1/6π)·(K IC /σ y ) 2 = 0.074 mm. According to Williams' formulas in front of the plastic zone in the crack plane direction, the level of stress tensor components is σ xx = σ yy = 1808 MPa, σ zz = ν·(σ xx +σ yy ) = 1193 MPa. Directly in the plastic zone the stress components are significantly higher and the stress opening component, σ yy , can obtain values higher than 3σ y , i.e., over 4500 MPa [38]. The σ zz and σ xx components are lower, but also exceed the level of 2σ y .
High stress levels in front of the crack tip cause material damage processes, cracking development through various mechanisms: brittle fracture on the cleavage planes, brittle intergranular cracking, ductile fracture through nucleation, growth and coalescence of voids. The predominance of any fracture mechanism depends on the stress-strain state in the tested element, on the microstructure type, loading method, the temperature of the environment, and on other effects. Here, the situation was seemingly simple-SENB specimens crack according to the brittle mechanism. The AE method was used to monitor the processes during loading of steel specimens.
Five specimens were tested. The load plots and the character of AE signals distribution were similar in all tested samples. Figure 3b shows the characteristic loading curve and numerous points of the absolute energy of AE signal and the graph illustrating the cumulative absolute energy vs. load on time from 3.5 to 6.2 s. The graphs show AE parameters without using the clustering method (identification) of destructive processes. Figure 3a shows the change in the absolute energy parameter (aJ), and the specimen loading as a function of time and COD (crack opening displacement). Acoustic emission signals generated during the study showed an increase in their numbers and increase in the value of absolute energy parameter at COD of 0.045 (mm) and 4.2 kN of load. It can be noticed that the increase in the value of the absolute energy parameter occurs when the cracking process begins (initiation, propagation). parameter at COD of 0.045 (mm) and 4.2 kN of load. It can be noticed that the increase in the value of the absolute energy parameter occurs when the cracking process begins (initiation, propagation). During the observations of the fracture surface of the specimens examined using the scanning microscope, three basic mechanisms in the crack development process were identified. The images shown in Figures 4 are presents the same fragment the fracture surface, but at different magnitude. In Figure 4a, at 500×, we can see the general character of brittle fracture surface. In Figure 4b, at 1500×, the details of fracture mechanisms are shown. At the fracture surface, two types of brittle fracture were present: intercrystallite cracking and transcrystallite cleavage. (Figure 4a,b). There were also areas where cracking developed according to the ductile mechanism, that is, where voids nucleated, grew, and coalesced around the small particles of precipitations (Figure 4b).
(a) 500× increasing the image (b) 1500× increasing the image Observations of the process of cracking over time shows that at first cracking proceeds according to brittle mechanisms, while cracking due to sudden immediate coalescence of voids is the last stage of crack growth. These observations were confirmed by the analysis of AE signals generated in the tested SENB specimen. However, in order to be able to detect the differences in the processes accompanying the cracking of steel specimens, the signal grouping method must be used. Figure 4a, b shows the selected AE parameters (absolute energy) not analyzed using the clustering methods. It can be noticed that it is difficult to interpret the changes in the work of the tested element. It is During the observations of the fracture surface of the specimens examined using the scanning microscope, three basic mechanisms in the crack development process were identified. The images shown in Figure 4 are presents the same fragment the fracture surface, but at different magnitude. In Figure 4a, at 500×, we can see the general character of brittle fracture surface. In Figure 4b, at 1500×, the details of fracture mechanisms are shown. At the fracture surface, two types of brittle fracture were present: intercrystallite cracking and transcrystallite cleavage. (Figure 4a,b). There were also areas where cracking developed according to the ductile mechanism, that is, where voids nucleated, grew, and coalesced around the small particles of precipitations (Figure 4b).
Materials 2019, 12, x FOR PEER REVIEW 8 of 14 parameter at COD of 0.045 (mm) and 4.2 kN of load. It can be noticed that the increase in the value of the absolute energy parameter occurs when the cracking process begins (initiation, propagation).
(a) (b) Figure 3. The graph of specimen loading (a) with absolute energy of AE signals (aJ) and (b) graph cumulative absolute energy vs. load on time from 3.5 to 6.2 s. The graphs show AE parameters without using the clustering method (identification) of destructive processes.
During the observations of the fracture surface of the specimens examined using the scanning microscope, three basic mechanisms in the crack development process were identified. The images shown in Figures 4 are presents the same fragment the fracture surface, but at different magnitude. In Figure 4a, at 500×, we can see the general character of brittle fracture surface. In Figure 4b, at 1500×, the details of fracture mechanisms are shown. At the fracture surface, two types of brittle fracture were present: intercrystallite cracking and transcrystallite cleavage. (Figure 4a,b). There were also areas where cracking developed according to the ductile mechanism, that is, where voids nucleated, grew, and coalesced around the small particles of precipitations (Figure 4b).
(a) 500× increasing the image (b) 1500× increasing the image Observations of the process of cracking over time shows that at first cracking proceeds according to brittle mechanisms, while cracking due to sudden immediate coalescence of voids is the last stage of crack growth. These observations were confirmed by the analysis of AE signals generated in the tested SENB specimen. However, in order to be able to detect the differences in the processes accompanying the cracking of steel specimens, the signal grouping method must be used. Figure 4a, b shows the selected AE parameters (absolute energy) not analyzed using the clustering methods. It can be noticed that it is difficult to interpret the changes in the work of the tested element. It is Observations of the process of cracking over time shows that at first cracking proceeds according to brittle mechanisms, while cracking due to sudden immediate coalescence of voids is the last stage of crack growth. These observations were confirmed by the analysis of AE signals generated in the tested SENB specimen. However, in order to be able to detect the differences in the processes accompanying the cracking of steel specimens, the signal grouping method must be used. Figure 4a, b shows the selected AE parameters (absolute energy) not analyzed using the clustering methods. It can be noticed that it is difficult to interpret the changes in the work of the tested element. It is impossible to assess the moment of initiation of the crack or its subsequent propagation. It can be said that individual parameters not analyzed by clustering methods do not help in the interpretation of changes in the tested samples. That is why it is so important to use Big Data analysis to create and identify destructive processes on the basis of the pattern recognition methods.
Using one of the grouping methods, namely the fuzzy k-means algorithm, we developed the reference signal data base describing the processes of brittle fracture during destructive tests of steel structures. The AE signals generated and recorded in the crack initiation and development area were subjected to a grouping analysis using the NOESIS 5.8 program, which allowed to separate four basic classes. Figure 5a-c present the effects of using method of grouping signals-a non-hierarchical k-means method to identify mechanisms accompanying the brittle fracture.
Materials 2019, 12, x FOR PEER REVIEW 9 of 14 impossible to assess the moment of initiation of the crack or its subsequent propagation. It can be said that individual parameters not analyzed by clustering methods do not help in the interpretation of changes in the tested samples. That is why it is so important to use Big Data analysis to create and identify destructive processes on the basis of the pattern recognition methods. Using one of the grouping methods, namely the fuzzy k-means algorithm, we developed the reference signal data base describing the processes of brittle fracture during destructive tests of steel structures. The AE signals generated and recorded in the crack initiation and development area were subjected to a grouping analysis using the NOESIS 5.8 program, which allowed to separate four basic classes. Figure 5a-c present the effects of using method of grouping signals-a non-hierarchical kmeans method to identify mechanisms accompanying the brittle fracture.  Analyzing the division of AE signals into the classes marked in Figure 5 with colors and figures and using photographs (Figure 4) from the scanning microscope, it can be noticed that cracking due to sudden immediate combination of the voids, breaking bridges between them are pink signals marked with the number 3. The same type of signals recorded in the last stage of cracking the sample tensile steel S355JR, when there was no cracking at the brittle mechanisms ( Figure 6). Signals marked in blue and number 1 come mainly from elastic and plastic deformations, as well as from individual noise generated by the strength machine and the loading system. These signals, therefore, probably can be equated with the complex, rapid dislocation motion process in ferrite and the act of fracture.  We can, therefore, assume that the signals of green color characterize the process of dislocation motion, and cracking on the cleavage surfaces in ferrite, i.e., brittle cleavage fracture. The red signals show different frequencies at different stages of the load frequency band. Parameter (IF) signals Analyzing the division of AE signals into the classes marked in Figure 5 with colors and figures and using photographs (Figure 4) from the scanning microscope, it can be noticed that cracking due to sudden immediate combination of the voids, breaking bridges between them are pink signals marked with the number 3. The same type of signals recorded in the last stage of cracking the sample tensile steel S355JR, when there was no cracking at the brittle mechanisms ( Figure 6). Signals marked in blue and number 1 come mainly from elastic and plastic deformations, as well as from individual noise generated by the strength machine and the loading system. These signals, therefore, probably can be equated with the complex, rapid dislocation motion process in ferrite and the act of fracture. Analyzing the division of AE signals into the classes marked in Figure 5 with colors and figures and using photographs (Figure 4) from the scanning microscope, it can be noticed that cracking due to sudden immediate combination of the voids, breaking bridges between them are pink signals marked with the number 3. The same type of signals recorded in the last stage of cracking the sample tensile steel S355JR, when there was no cracking at the brittle mechanisms ( Figure 6). Signals marked in blue and number 1 come mainly from elastic and plastic deformations, as well as from individual noise generated by the strength machine and the loading system. These signals, therefore, probably can be equated with the complex, rapid dislocation motion process in ferrite and the act of fracture.  Figure 7a,b presents the parameters of the acoustic emission signal: initiation frequency (IF) (kHz) and reverberation frequency (RF) (kHz) vs. time and signal classes. Analysing these graphs, it can be seen that green signals have a similar frequency range in the initial stage (IF) parameter as well as in the later stage (RF) parameter during signal duration. This range falls within the frequency range of pink signals.
We can, therefore, assume that the signals of green color characterize the process of dislocation motion, and cracking on the cleavage surfaces in ferrite, i.e., brittle cleavage fracture. The red signals show different frequencies at different stages of the load frequency band. Parameter (IF) signals  Figure 7a,b presents the parameters of the acoustic emission signal: initiation frequency (IF) (kHz) and reverberation frequency (RF) (kHz) vs. time and signal classes. Analysing these graphs, it can be seen that green signals have a similar frequency range in the initial stage (IF) parameter as well as in the later stage (RF) parameter during signal duration. This range falls within the frequency range of pink signals.
We can, therefore, assume that the signals of green color characterize the process of dislocation motion, and cracking on the cleavage surfaces in ferrite, i.e., brittle cleavage fracture. The red signals show different frequencies at different stages of the load frequency band. Parameter (IF) signals occurring in the earlier stage of the load have much higher than others, which allows to assume that these signals are generated in a different environment than ferrite. These may be signals from the cracking of precipitates segregated at the ferrite grain boundaries or connected with the process of de-bonding from the ferrite matrix. These signals can be imitated, in ferrite, and in precipitations. That is why there is a different frequency in parameter (IF).
The signal imitated in the precipitations should be short. However, the signal imitated in ferrite is longer and therefore parameter (RF) is similar. In the next stage of loading, red signals characterize the development of intercrystallite fracture growth and are imitated in ferrite.
Materials 2019, 12, x FOR PEER REVIEW 11 of 14 occurring in the earlier stage of the load have much higher than others, which allows to assume that these signals are generated in a different environment than ferrite. These may be signals from the cracking of precipitates segregated at the ferrite grain boundaries or connected with the process of de-bonding from the ferrite matrix. These signals can be imitated, in ferrite, and in precipitations. That is why there is a different frequency in parameter (IF). The signal imitated in the precipitations should be short. However, the signal imitated in ferrite is longer and therefore parameter (RF) is similar. In the next stage of loading, red signals characterize the development of intercrystallite fracture growth and are imitated in ferrite.

General Conclusions
To raise AE technology to the next level, we must improve it from several fronts. Only then can AE become a serious player in structural integrity assessment. We need to better characterize AE

General Conclusions
To raise AE technology to the next level, we must improve it from several fronts. Only then can AE become a serious player in structural integrity assessment. We need to better characterize AE sources; to extract more from AE signals reaching sensors; to devise effective AE parameters; to integrate new analysis methods, such as wave transformation analysis(WT), moment tensor analysis (MTA), and pattern recognition analysis (PRA); to include simulation tools in analyses; to accumulate basic data on structures with standardized procedures; and to devise combinatorial approach between localized damage evaluation and long-range detection and to develop regional or global database under international cooperation. Using systematic approach along with NDT, fracture mechanics, etc., improved AE can offer more value to industrial users and to contribute substantially to structural integrity assessment. Development of the proposed AE monitoring technique reported in this paper facilitates for the prognostics and life predictions of the structure. The developed methodology can be utilized for continuous in-service monitoring of structures and has proven to be promising for use in practice.

1.
Application of the k-means grouping method based on the analysis in the 13 parameter space of AE signals provides positive results in the classification of destructive processes in fracture toughness K Ic tests.

2.
Application of the k-means grouping method allows differentiating between brittle fractures mechanisms. 3.
The process of brittle fracture was shown to be a sequential process composed of inter-grain cracking, brittle-cleavage and ductile fracture.

4.
Acoustic emission technology was found to be able to successfully monitor and detect real-time crack initiation and also separate and identify AE signals that correspond to different fracture mechanisms 5.
The grouping analysis shows that some AE parameters noticeably change after the initiation of the crack (absolute energy, initial frequency, reverberation frequency), which indicates their high independence and suitability for monitoring the cracking process. 6.
The developed base of reference signals forms the theoretical basis for using pattern-based AE technology to detect and monitor crack propagation in structures that do not have a "load history" to assess the safety of steel structures using the principles of fracture mechanics.