Damage Mechanism Evaluation of Large-Scale Concrete Structures Affected by Alkali-Silica Reaction Using Acoustic Emission

: Alkali-silica reaction has caused damage to concrete structures, endangering structural serviceability and integrity. This is of concern in sensitive structures such as nuclear power plants. In this study, acoustic emission (AE) was employed as a structural health monitoring strategy in large-scale, reinforced concrete specimens affected by alkali-silica reaction with differing boundary conditions resembling the common conditions found in nuclear containments. An agglomerative hierarchical algorithm was utilized to classify the AE data based on energy-frequency based features. The AE signals were transferred into the frequency domain and the energies in several frequency bands were calculated and normalized to the total energy of signals. Principle component analysis was used to reduce feature redundancy. Then the selected principal components were considered as features in an input of the pattern recognition algorithm. The sensor located in the center of the conﬁned specimen registered the largest portion of AE energy release, while in the unconﬁned specimen the energy is distributed more uniformly. This conﬁrms the results of the volumetric strain, which shows that the expansion in the conﬁned specimen is oriented along the thickness of the specimen.


Introduction
Alkali-silica reaction (ASR) is a chemical processes that has caused damage in concrete structures such as bridges [1][2][3][4], nuclear power plants [5][6][7][8][9], and concrete dams [5,10]. This reaction usually occurs between alkali hydroxides in the pore solution and the reactive silica in some aggregates. The result is an alkali-silica gel, which is hygroscopic and imbibes water and humidity. The gel expands in humidity exceeding 80% [6]. This expansion induces pressure on the concrete matrix and aggregates and causes micro-cracks and cracks when the pressure exceeds the tensile strength of the concrete [11,12]. Several traditional methods such as visual inspection, coring, and petrographic analysis have been utilized for monitoring and identifying the behavior of damage caused by ASR. Visual inspection is very subjective and depends on experience [13]. In addition, damage can only be detected by using this method, once it has reached the surface. Coring is a destructive method, which is discouraged in sensitive structures like nuclear power plants. Moreover, it is difficult to generalize the global condition of structures through evaluation of a few cores. Petrographic analysis is also

Specimen and Chamber Preparation
In this study, three large-scale reinforced concrete blocks were cast. Two specimens were reactive and one was a control specimen ( Figure 1). Reactive specimens are the specimens that experience expansion due to alkali-silica reaction. Two layers of steel reinforcement mesh were located at the top and bottom surfaces of the specimens. The reinforcement mesh includes US #11 Grade 60 with a nominal diameter of 36 mm (1.41 inches) at 25 cm (10 inches) spacing ( Figure 2). The specimens have 7.62 cm (3 inches) cover at top and bottom of the specimens. Moreover, the square steel plates were installed on the ends of the rebar to achieve the full development length in a relatively short distance inside the specimens. For concrete mixture, 350 kg/m 3 cement, 175 kg/m 3 water, 1180 kg/m 3 coarse aggregate, and 728 kg/m 3 fine aggregate were used for both reactive and control specimens. The main difference between the specimens was using of NaOH or LiNO 3 . In the reactive specimens, 9.8 kg/m 3 NaOH solution was added, while in the control specimen 1.9 kg/m 3 LiNO 3 was included in the mixture [22].
The coarse aggregates in both control and reactive specimens were greenschist coarse aggregates from North Carolina, which were highly-reactive. Non-reactive manufactured sand and low-alkali Portland cement Type II were also utilized in the concrete mixture. Water to cement weight ratio for both control and reactive specimens was 50%. To increase the alkali content of reactive specimens, 50% sodium hydroxide solution was added to the mixture for reactive specimens. Additionally, a 30% lithium nitrate solution was added to the control specimens to mitigate the ASR effect. More details regarding concrete mixture design was presented by Hayes et al. [22].
The specimens shown in Figure 1 have a cubic shape with dimensions of 3.50 m × 3.0 m × 1.0 m (136 inch × 116 inch × 40 inch). There are two reactive specimens; confined and unconfined specimens. A rigid steel brace was utilized to restrain the confined specimen from expanding in the lateral plane. The unconfined specimen does not have a steel brace in its plane. The two reactive specimens resemble the different stress boundary conditions found in real structures. A 1.5 mm layer of polyethylene was employed between the concrete and the steel frame in the confined specimen to reduce friction. More detailed information about the specimens and steel frame design were presented in Hayes, et al. [22]. In this experiment, the specimens stand on four short, steel columns (76 cm height) to provide access to the bottom of the specimens.
Each specimen's formwork was removed 12 days after casting which gave them time for final finishing and concrete setting. After casting, the specimen surfaces were sprayed with curing compound and covered with wet burlap and plastic sheets to reduce moisture loss and cracks from shrinkage during drying. The burlap was kept wet until chamber construction was complete [22].
A large chamber with the dimensions of 16.2 m long, 7.3 m wide, and 3.7 m high was constructed around the specimens to control the temperature and relative humidity of the environment. The chamber was equipped with systems that kept the specimens at the environmental temperature of 38 ± 1 • C and humidity of 95% ± 5%. The environmental chamber was initiated 26 days after casting and operated continuously except for shutdown for measurements and inspections [22]. Appl

Measurement Equipment
To measure the strain inside the specimens, two kinds of strain sensors, alkali-tolerant strain transducers (KM-100B from Tokyo Sokki Kenkyujo) and long-gauge fiber optic extensometers (from SMARTECH/Roctest), were employed. In total, sixty-four 100 mm-gauge strain transducers were placed inside the specimens, along with the specimen dimensions. The sensors were fixed on the support fabricated by smooth steel bars with a diameter of 3 mm. Additionally, five high-precision, fiber optic extensometers were installed inside the specimens. The sensors in the X-Y plane had 1.5 m and 1.8 m gauge lengths and were located at the bottom of the specimens. Two fiber optic sensors were installed along the specimen thickness with a gauge length of 0.8 m [22].
Seven acoustic emission sensors were utilized in each of the confined and unconfined specimens (reactive specimens) as shown in Figures 1 and 2. Two sensors were employed for the control specimen. The internal AE sensors were broadband WDIUC-AST (manufactured by MISTRAS Group, Inc., New Jersey, USA), with an operating frequency range of 200-900 kHz and an internal low-noise 40 dB preamplifier. The external sensors were resonant R6I-UC (manufactured by MISTRAS Group, Inc., New Jersey, USA), with an operating frequency range of 35-100 kHz and an internal 40 dB preamplifier. The special polymer coating on the sensors (internal and external) with an internal waterproof cable makes the sensors insulated, non-conductive, and capable of operating under water. Three broadband sensors were installed inside the reinforcement cages prior to casting for reactive specimens and one sensor was installed in the control specimen. Sensors 1, 2, 4, and 5 were attached under the second layer of top reinforcement mesh (under reinforcement in the Y direction) with epoxy and fastened to the rebar by several cable ties and duct tape, with the sensing surface

Measurement Equipment
To measure the strain inside the specimens, two kinds of strain sensors, alkali-tolerant strain transducers (KM-100B from Tokyo Sokki Kenkyujo) and long-gauge fiber optic extensometers (from SMARTECH/Roctest), were employed. In total, sixty-four 100 mm-gauge strain transducers were placed inside the specimens, along with the specimen dimensions. The sensors were fixed on the support fabricated by smooth steel bars with a diameter of 3 mm. Additionally, five high-precision, fiber optic extensometers were installed inside the specimens. The sensors in the X-Y plane had 1.5 m and 1.8 m gauge lengths and were located at the bottom of the specimens. Two fiber optic sensors were installed along the specimen thickness with a gauge length of 0.8 m [22].
Seven acoustic emission sensors were utilized in each of the confined and unconfined specimens (reactive specimens) as shown in Figures 1 and 2. Two sensors were employed for the control specimen. The internal AE sensors were broadband WDIUC-AST (manufactured by MISTRAS Group, Inc., Princeton Junction, NJ, USA), with an operating frequency range of 200-900 kHz and an internal low-noise 40 dB preamplifier. The external sensors were resonant R6I-UC (manufactured by MISTRAS Group, Inc., Princeton Junction, NJ, USA), with an operating frequency range of 35-100 kHz and an internal 40 dB preamplifier. The special polymer coating on the sensors (internal and external) with an internal waterproof cable makes the sensors insulated, non-conductive, and capable of operating under water.
Three broadband sensors were installed inside the reinforcement cages prior to casting for reactive specimens and one sensor was installed in the control specimen. Sensors 1, 2, 4, and 5 were attached under the second layer of top reinforcement mesh (under reinforcement in the Y direction) with epoxy and fastened to the rebar by several cable ties and duct tape, with the sensing surface facing downward along the negative Z axis (Figures 1d and 2d). Sensors 3 and 6 were in the middle of the specimen thickness, with the sensing surface facing toward the negative X axis (Figure 1d). Sensors 3 and 6 were attached on the supports by epoxy and duct tape and cables were fastened on the support by cable ties (Figure 2c). The supports were fastened to the top and bottom reinforcement meshes. The locations of internal sensors were colored on the top reinforcements to avoid over-vibrating during casting and potentially decoupling the sensors.
Four resonant sensors were attached at the bottom of the reactive specimens and one to the control specimen ( Figure 2b). The resonant sensors were attached with epoxy on the surface of the concrete and fixed by holders as shown in Figure 2b. Internal and external sensors are referred to as broadband and resonant sensors, respectively, in this paper. The holders were made of stainless steel and polyurethane pads were attached between the concrete and the holders to avoid potential corrosion in the holders.
External and internal sensors pose differing technical challenges. Internal sensors are of scientific interest and are not well-suited for implementation in existing structures. One potential technical consideration is the formation of voids around the sensors due to concrete shrinkage. This may result in interruption of data or reduction of sensitivity. Another technical challenge is lack of access after casting. The external sensors are accessible, but decoupling may be a consideration.
A 16-channel Sensor Highway II (SHII), manufactured by MISTRAS Group, Inc. (Princeton Junction, NJ, USA), was utilized as a data acquisition system. The sensitivity of external sensors was checked by applying Hsu-Nielsen sources [23].
The sampling rate was set to 1000 kHz. Pre-trigger time is the period required for data acquisition to save a signal prior to threshold intersection, which was set to 256 µs. HDT (hit definition time) is the time for which if a signal crosses the threshold the signal will continue until the end of HDT without stopping. HDT was set to 400 µs. HLT (hit lockout time) is the time defined at the end of a signal for neglecting any reflected signal that exceeds the threshold, which was set to 200 µs. PDT (peak definition time) is the time which ensures correct identification of the signal peak for rise time and peak amplitude measurements [24]. In this study, PDT was set to 200 µs. The initial threshold was 32 dB.

Pattern Recognition Algorithm
To identify the damage mechanisms in the reactive specimens a pattern recognition algorithm was employed for data classification. Pattern recognition is under the machine learning field and has two major types: unsupervised and supervised. When a background history for a data set is available or the data set has labeled classes, supervised pattern recognition is employed. If there are no labeled classes available, unsupervised pattern recognition can be utilized to identify the potential patterns in the data set based on the selected features.
An agglomerative hierarchical clustering algorithm [25] was employed for classifying the AE data into subsets. The clustering procedure is chronologically illustrated in Figure 3. The first step in the pattern recognition was deriving the frequency-based features for the signals. AE signals were transferred into the frequency domain using the Fast Fourier Transform (FFT). The FFT amplitude spectra were determined for each AE signal and the frequency domain was divided into ten frequency ranges. Then, signal energy in each frequency range was calculated and normalized to the total energy of the signal. These ten signal energies are the signal features. For example, in Figure 3a, the area of the hatched region is the energy of the signal in the frequency band from 200 kHz to 250 kHz. This value is normalized to the total energy of the signal, which is the entire area under the FFT spectrum. Principle component analysis (PCA) was used to reduce redundancy in the data. In this analysis, the original data is projected to the new orthogonal coordinates having high variation. An input for PCA is a matrix with the number of columns and rows equal to the feature numbers and the number of hits. Then a variance-covariance matrix for the features was calculated. The coefficients and variance of a specific principle component were calculated by eigenvalue analysis of the variance-covariance matrix. The principal components were selected in a way that represented more than 90% of the entire data. The principal components of the original features were selected as the input features for the pattern recognition algorithm. The algorithm initially calculated the Euclidian distance between the resulting data from PCA analysis. The result was a proximity matrix that contained distances between the original objects (data).
The objects were initially linked together according to calculated distances in the previous step and Ward's method. Ward's method is based on calculating the total within-cluster sum of squares of the data resulting from combining the clusters [25]. In each level, the data was merged into a binary linkage and the clusters were again merged into new clusters according to the Ward's method. This procedure was continued to form a single cluster which includes all data. The number of the cluster was determined according to the developed dendrogram and the height of each link with respect to the adjacent links [26]. Appl

Results and Discussions
The internal strains of the specimens in X, Y, and Z directions were recorded from the casting date. The results are presented up to 200 days in Figure 4. The data is presented from 51 days in the figures since external sensors were attached in this day and all the figures in this paper remain in a same time scale for a convenient comparison. Figure 4a shows the average expansion along the specimen dimensions for the reactive specimens and the average volumetric expansions for both reactive and control specimens are presented in Figure 4b. The results for volumetric strains show that reactive specimens were expanded. Conversely, the control specimens had no expansion and the internal strain was relatively constant, only showing small shrinkage over time. However, the volumetric strains for both confined and unconfined specimens were very close. This observation indicates that an imposed confinement in the specimens did not influence the expansion of the entire specimens (volumetric strain). Furthermore, the confinement changed the direction of expansion (Figure 4a).

Results and Discussions
The internal strains of the specimens in X, Y, and Z directions were recorded from the casting date. The results are presented up to 200 days in Figure 4. The data is presented from 51 days in the figures since external sensors were attached in this day and all the figures in this paper remain in a same time scale for a convenient comparison. Figure 4a shows the average expansion along the specimen dimensions for the reactive specimens and the average volumetric expansions for both reactive and control specimens are presented in Figure 4b. The results for volumetric strains show that reactive specimens were expanded. Conversely, the control specimens had no expansion and the internal strain was relatively constant, only showing small shrinkage over time. However, the volumetric strains for both confined and unconfined specimens were very close. This observation indicates that an imposed confinement in the specimens did not influence the expansion of the entire specimens (volumetric strain). Furthermore, the confinement changed the direction of expansion (Figure 4a).

Results and Discussions
The internal strains of the specimens in X, Y, and Z directions were recorded from the casting date. The results are presented up to 200 days in Figure 4. The data is presented from 51 days in the figures since external sensors were attached in this day and all the figures in this paper remain in a same time scale for a convenient comparison. Figure 4a shows the average expansion along the specimen dimensions for the reactive specimens and the average volumetric expansions for both reactive and control specimens are presented in Figure 4b. The results for volumetric strains show that reactive specimens were expanded. Conversely, the control specimens had no expansion and the internal strain was relatively constant, only showing small shrinkage over time. However, the volumetric strains for both confined and unconfined specimens were very close. This observation indicates that an imposed confinement in the specimens did not influence the expansion of the entire specimens (volumetric strain). Furthermore, the confinement changed the direction of expansion (Figure 4a).   Generally, the expansion through the thickness of the reactive specimens (Z direction) is much larger than the expansion throughout the X-Y plane, due to partial or complete confinement in the plane and lack of steel reinforcement through the thickness. The maximum strains in the Z direction for confined and unconfined specimens are approximately 3.6 and 2.1 times the strain in the Y direction, respectively, at 200 days. The steel frame in the confined specimen causes a reduction of both in-plane strains. In the confined specimen, the strain along the X direction is less than the Y direction, while in the unconfined specimen, the strain along X and Y are almost the same. The strain along the X direction in the confined specimen is 58% of the strain along the X direction in the unconfined specimen at 200 days. However, the strain along the Y direction in the confined specimen is 84% of the strain along the Y direction in the unconfined specimen at 200 days. Moreover, the confinement in X-Y plane for the confined specimen caused an increased strain rate along the Z direction (thickness). The expansion strain rate in terms of time decreases at the point around 150 days, which is named inflection point as seen in Figure 4. The inflection point in volumetric strain is considered as one of the important point for ASR modeling, where the curvature of volumetric strain is changed [27,28]. The latency and characteristic times (two modeling parameters) are experimentally determined by knowing the location of inflection point [6,28]. The inflection point is shown in Figure 4b and marked in other figures in this paper.
Acoustic emission data was recorded through the internal sensors from the casting day and resonant sensors started recording from the concrete age of 51 days. Filtering AE data is an important step for reducing the amount of non-relevant data for post-processing. The possible sources of false AE data can be friction between the structural components, such as the steel frame and specimens, and water dripping from the chamber ceiling due to high humidity. Two different filtering procedures were developed to minimize the non-genuine data for the internal and external sensors, separately. The filtering is different for the internal and external sensors due to their differences in sensitivity and location. For instance, the internal sensors are much less sensitive than the external sensors and were located inside the specimens, thereby receiving less environmental noise. The AE data below 32 dB and 41 dB for the internal and external sensors was filtered from the data set. Then, the signals were further filtered using a Swansong filtering procedure [29,30]. This method is based on the observation that genuine AE signals with high amplitude should have long duration and vice versa [11]. Therefore, false signals in this method are categorized by long duration with low amplitude and short duration with high amplitude. The data is presented in terms of duration versus amplitude distribution. Signals which did not comply with the characteristics of genuine signals were deleted based on visual observation of the waveforms. The rejection limits for the internal and external sensors are presented in Table 1. In addition, suspicious signals were removed by inspection of waveforms and the chamber activity timetable provided by the University of Tennessee, Knoxville. The filtered AE data from 50 to 195 days after casting for all sensors is presented in Figure 5. The figure illustrates the amplitude and cumulative signal strength (CSS) in terms of specimen age. The time window of 50 to 195 days was selected because the highest strain change rate occurred during this period. The gap in the data from 65 to 100 days is missing data due to a difficulty with the data acquisition system caused by an unexpected energy surge. The amount of data in the control specimen is much lower than the data for the reactive specimens. The difference between AE data for the control and reactive specimens is also observable from plots of cumulative signal strength versus time ( Figure 5). Therefore, the relatively high AE activity for the reactive specimens can likely be associated with expansion caused by the alkali-silica reaction.

Acoustic Emission Energy Release
To enable comparisons between the confined and unconfined specimens only AE data recorded by the internal sensors is discussed. These sensors are less sensitive than the resonant sensors, which results in lower volumes of data. The main reason for choosing the broadband sensors is that their broad frequency response makes them suitable for frequency analysis. The resonant sensors attached on the bottom surface of the specimens are more representative of what may be used in practice for optimized detection and source location. Three dimensional source location of large-scale specimens requires a method to accurately calculate time of arrival for very weak signals, which is a future step of this research. This study involves frequency analysis and therefore focuses on broad band sensors. Analysis of the resonant sensor data is a future consideration.
In Figure 1, the red marks and corresponding labels in red font denotes the schematic sensor locations for both confined and unconfined specimens. Sensor coordinates are presented in Table 2. The last column of the table represents the coordinates of normal vectors, which are perpendicular to

Acoustic Emission Energy Release
To enable comparisons between the confined and unconfined specimens only AE data recorded by the internal sensors is discussed. These sensors are less sensitive than the resonant sensors, which results in lower volumes of data. The main reason for choosing the broadband sensors is that their broad frequency response makes them suitable for frequency analysis. The resonant sensors attached on the bottom surface of the specimens are more representative of what may be used in practice for optimized detection and source location. Three dimensional source location of large-scale specimens requires a method to accurately calculate time of arrival for very weak signals, which is a future step of this research. This study involves frequency analysis and therefore focuses on broad band sensors. Analysis of the resonant sensor data is a future consideration.
In Figure 1, the red marks and corresponding labels in red font denotes the schematic sensor locations for both confined and unconfined specimens. Sensor coordinates are presented in Table 2. The last column of the table represents the coordinates of normal vectors, which are perpendicular to the sensing surfaces of sensors and the directions of vectors are toward the outside of the sensors. The orientation of sensors is shown in Figure 1d. The AE cumulative signal strengths of the internal sensors for the reactive specimens are presented in Figure 6a,b. The total CSS versus time for the confined specimen is much higher than the CSS for the unconfined specimen at 195 days (the CSS for the confined specimen is 2.35 times the value for the unconfined specimen). Moreover, in Figure 5a,c, it can be seen that a significant portion of released AE energy is attributed to sensor 3, which is located at mid-height of the confined specimen. The CSS rate for sensor 3 (5805 pVs/Day) is much larger than for sensors 1 and 2 (1545 and 1212 pVs/Day), thereby increasing the difference in the CSS between the sensors. On the other hand, this trend in the confined specimen is not observed in the unconfined specimen where AE energy release is not concentrated in one specific sensor. The distribution of energy was approximately uniform among the sensors. However, the CSS for the sensors at the top reinforcement mesh (sensors 4 and 5) is larger than in sensor 6. the sensing surfaces of sensors and the directions of vectors are toward the outside of the sensors. The orientation of sensors is shown in Figure 1d. The AE cumulative signal strengths of the internal sensors for the reactive specimens are presented in Figure 6a,b. The total CSS versus time for the confined specimen is much higher than the CSS for the unconfined specimen at 195 days (the CSS for the confined specimen is 2.35 times the value for the unconfined specimen). Moreover, in Figure 5a,c, it can be seen that a significant portion of released AE energy is attributed to sensor 3, which is located at mid-height of the confined specimen. The CSS rate for sensor 3 (5805 pVs/Day) is much larger than for sensors 1 and 2 (1545 and 1212 pVs/Day), thereby increasing the difference in the CSS between the sensors. On the other hand, this trend in the confined specimen is not observed in the unconfined specimen where AE energy release is not concentrated in one specific sensor. The distribution of energy was approximately uniform among the sensors. However, the CSS for the sensors at the top reinforcement mesh (sensors 4 and 5) is larger than in sensor 6. These observations illustrate that the confined specimen has larger AE energy release in the middle layers of the specimen, which is increasing with the progression of the ASR reaction, than the unconfined specimen. This may be due to a larger expansion strain through the thickness of the These observations illustrate that the confined specimen has larger AE energy release in the middle layers of the specimen, which is increasing with the progression of the ASR reaction, than the unconfined specimen. This may be due to a larger expansion strain through the thickness of the confined specimen. This large expansion is expected to cause more damage through the thickness of specimen and consequently more AE energy release at middle of the thickness. In the unconfined specimen the crack distribution is expected to be less anisotropic than in the confined specimen. Therefore, it is expected that the AE energy was more uniformly distributed through the thickness of this specimen.

Pattern Recognition of AE Data
The AE signals were classified according to the agglomerative hierarchical algorithm as explained in Section 3.1. Dendrograms resulting from the analysis are presented in Figure 7. Dashed red lines in the subfigures indicate the desired height of links for clustering. The results of the pattern recognition show three clusters for each reactive specimen (both confined and unconfined). The clusters of confined specimens are indicated by Cluster-1, Cluster-2, and Cluster-3. Accordingly, the clusters of the unconfined specimen are labeled Cluster-4, Cluster-5, and Cluster-6. The horizontal axis in Figure 7 gives the data labels, which show either labels of the original data sets (signal number) or the label number of the clusters that resulted from merging the original data. The height of each link shows the distance between the two objects. Each link between two objects is shown by an upside-down U-shaped line in the figure. The data is shown in terms of the first three principal components (PC) to visualize the distribution of clusters with respect to each other in Figure 8. Although some signals in the cluster 5 and 6 were not ideally separated, in general, the clusters indicate a reasonable separation in the PC space ( Figure 8). The reason for an overlap between clusters 5 and 6 is that some signals in cluster 5 have a similar energy contribution in a specific frequency range to some signals in cluster 6. As seen in Figure 9b, clusters 5 and 6 have similar average energy contribution between 250 to 300 kHz. confined specimen. This large expansion is expected to cause more damage through the thickness of specimen and consequently more AE energy release at middle of the thickness. In the unconfined specimen the crack distribution is expected to be less anisotropic than in the confined specimen. Therefore, it is expected that the AE energy was more uniformly distributed through the thickness of this specimen.

Pattern Recognition of AE Data
The AE signals were classified according to the agglomerative hierarchical algorithm as explained in Section 3.1. Dendrograms resulting from the analysis are presented in Figure 7. Dashed red lines in the subfigures indicate the desired height of links for clustering. The results of the pattern recognition show three clusters for each reactive specimen (both confined and unconfined). The clusters of confined specimens are indicated by Cluster-1, Cluster-2, and Cluster-3. Accordingly, the clusters of the unconfined specimen are labeled Cluster-4, Cluster-5, and Cluster-6. The horizontal axis in Figure 7 gives the data labels, which show either labels of the original data sets (signal number) or the label number of the clusters that resulted from merging the original data. The height of each link shows the distance between the two objects. Each link between two objects is shown by an upside-down U-shaped line in the figure. The data is shown in terms of the first three principal components (PC) to visualize the distribution of clusters with respect to each other in Figure 8. Although some signals in the cluster 5 and 6 were not ideally separated, in general, the clusters indicate a reasonable separation in the PC space ( Figure 8). The reason for an overlap between clusters 5 and 6 is that some signals in cluster 5 have a similar energy contribution in a specific frequency range to some signals in cluster 6. As seen in Figure 9b, clusters 5 and 6 have similar average energy contribution between 250 to 300 kHz. The average energy of signals in terms of frequency ranges is shown in Figure 9 for the reactive specimens. These values, as previously mentioned, were calculated using the FFT amplitude spectrum. Afterward, the calculated values were normalized by the total energy of the signal. The average values for each cluster were then calculated. The energy shift to the higher frequency components for the unconfined specimen is apparent when viewed alongside the confined specimen. The clusters in each specimen can be separated based on the frequency content. In the confined specimen, the low-frequency cluster (Cluster-1) has approximately 69% of its energy in a frequency The average energy of signals in terms of frequency ranges is shown in Figure 9 for the reactive specimens. These values, as previously mentioned, were calculated using the FFT amplitude spectrum. Afterward, the calculated values were normalized by the total energy of the signal. The average values for each cluster were then calculated. The energy shift to the higher frequency components for the unconfined specimen is apparent when viewed alongside the confined specimen. The clusters in each specimen can be separated based on the frequency content. In the confined specimen, the low-frequency cluster (Cluster-1) has approximately 69% of its energy in a frequency The average energy of signals in terms of frequency ranges is shown in Figure 9 for the reactive specimens. These values, as previously mentioned, were calculated using the FFT amplitude spectrum. Afterward, the calculated values were normalized by the total energy of the signal. The average values for each cluster were then calculated. The energy shift to the higher frequency components for the unconfined specimen is apparent when viewed alongside the confined specimen. The clusters in each specimen can be separated based on the frequency content. In the confined specimen, the low-frequency cluster (Cluster-1) has approximately 69% of its energy in a frequency range of 0-100 kHz. The medium-frequency cluster (Cluster-2) has 42% of its energy concentrated in a frequency range of 50-150, while the high-frequency cluster (Cluster-3) has 51% of its energy between 150-300 kHz.
In the unconfined specimen, the low-frequency cluster (Cluster-4) has 62% of its energy in a frequency range of 0-150 kHz and the medium-frequency cluster (Cluster-5) has 49% of its signal energy in the frequency range of 150-300 kHz. The high-frequency cluster (Cluster-6) has 54% of its signal energy concentrated between the frequencies of 300-450 kHz. The Cluster-3 and Cluster-5 share similar frequency content. More signal features for the clusters are illustrated in Figure 10. The average feature values for each cluster were normalized by the maximum feature values for each cluster. In the confined specimen, Cluster-1 initiated at a higher concrete age (data was analyzed through a concrete age of 195 days) than clusters with higher frequencies (Cluster-2 and Cluster-3). The average amplitude of the signals in Cluster-3 (the highest frequency) is higher than the other signals. Average signal strength for Cluster-1 is lower than the values for Cluster-2 and Cluster-3, and duration is higher for the cluster with the lowest frequency content, for example, Cluster-1. range of 0-100 kHz. The medium-frequency cluster (Cluster-2) has 42% of its energy concentrated in a frequency range of 50-150, while the high-frequency cluster (Cluster-3) has 51% of its energy between 150-300 kHz.
In the unconfined specimen, the low-frequency cluster (Cluster-4) has 62% of its energy in a frequency range of 0-150 kHz and the medium-frequency cluster (Cluster-5) has 49% of its signal energy in the frequency range of 150-300 kHz. The high-frequency cluster (Cluster-6) has 54% of its signal energy concentrated between the frequencies of 300-450 kHz. The Cluster-3 and Cluster-5 share similar frequency content. More signal features for the clusters are illustrated in Figure 10. The average feature values for each cluster were normalized by the maximum feature values for each cluster. In the confined specimen, Cluster-1 initiated at a higher concrete age (data was analyzed through a concrete age of 195 days) than clusters with higher frequencies (Cluster-2 and Cluster-3). The average amplitude of the signals in Cluster-3 (the highest frequency) is higher than the other signals. Average signal strength for Cluster-1 is lower than the values for Cluster-2 and Cluster-3, and duration is higher for the cluster with the lowest frequency content, for example, Cluster-1. A clear correlation is present between the frequency content of the signal clusters and the rise angle values (rise time over amplitude ratio) as has been observed by other researchers [31,32]. The higher the frequency components are in a signal, the lower rise angle value the signal possesses. In A clear correlation is present between the frequency content of the signal clusters and the rise angle values (rise time over amplitude ratio) as has been observed by other researchers [31,32]. The higher the frequency components are in a signal, the lower rise angle value the signal possesses. In the unconfined specimen, Cluster-5 and Cluster-6 exhibit higher hit rates at the higher concrete age compared to Cluster-4, and the average amplitude of signals for the clusters with higher frequency components is slightly higher than for signals in Cluster-4. However, the average duration for the signals in the cluster with the low-frequency components (Cluster-4) is much longer than the duration for Cluster-5 and Cluster-6.
In Figure 11, the variation of cumulative signal strength in terms of the age of the concrete for each cluster is presented. The cumulative signal strengths were normalized by the maximum value for each specimen. In the confined specimen, the signals with the highest frequency components (Cluster-3) have dominant CSS from the early age. However, the CSS of Cluster-2 is very close to the CSS of Cluster-3 up to the concrete age of 150 days. After 150 days, the CSS rate for Cluster-3 increases, while the CSS rate of Cluster-2 continues with approximately the same rate. The signals in cluster Cluster-1 have negligible signal strength compared to Cluster-2 and Cluster-3 and initiate primarily after 120 days. In the unconfined specimen, the AE energy is primarily attributed to cluster Cluster-4 up to approximately 150 days. After 150 days, the clusters with the higher frequency components (Cluster-5 and Cluster-6) become prominent in terms of AE energy release. the unconfined specimen, Cluster-5 and Cluster-6 exhibit higher hit rates at the higher concrete age compared to Cluster-4, and the average amplitude of signals for the clusters with higher frequency components is slightly higher than for signals in Cluster-4. However, the average duration for the signals in the cluster with the low-frequency components (Cluster-4) is much longer than the duration for Cluster-5 and Cluster-6. In Figure 11, the variation of cumulative signal strength in terms of the age of the concrete for each cluster is presented. The cumulative signal strengths were normalized by the maximum value for each specimen. In the confined specimen, the signals with the highest frequency components (Cluster-3) have dominant CSS from the early age. However, the CSS of Cluster-2 is very close to the CSS of Cluster-3 up to the concrete age of 150 days. After 150 days, the CSS rate for Cluster-3 increases, while the CSS rate of Cluster-2 continues with approximately the same rate. The signals in cluster Cluster-1 have negligible signal strength compared to Cluster-2 and Cluster-3 and initiate primarily after 120 days. In the unconfined specimen, the AE energy is primarily attributed to cluster Cluster-4 up to approximately 150 days. After 150 days, the clusters with the higher frequency components (Cluster-5 and Cluster-6) become prominent in terms of AE energy release. In the confined specimen, most of the energy contribution is related to Cluster-2 (energy concentration in 50-150 kHz) and Cluster-3 (energy concentration in 150-300 kHz). The SSCF for  SSCF of the clusters with high-frequency components (Cluster-5 and Cluster-6) increases with time. There is no obvious energy concentration in the sensor located at mid-thickness of the specimen (sensor 6), which is different from what is observed in the confined specimen. In both specimens, the SSCF declines in low-frequency signals and increases in high-frequency signals with time. This trend in signal frequency from low to high in the confined specimen is not pronounced before 150 days. The SSCF for cluster Cluster-3 is slightly greater than the Cluster-2 at 150 days. The signal frequency trend initiates primarily after 150 days in this specimen. On the other hand, in the unconfined specimen, the frequency content evolution of AE signals is obvious from an earlier stage of ASR reaction (66 days) and is more significant after 150 days. In the confined specimen, most of the energy contribution is related to Cluster-2 (energy concentration in 50-150 kHz) and Cluster-3 (energy concentration in 150-300 kHz). The SSCF for Cluster-3 is increases with time, particularly after 150 days. Most of the AE energy for Cluster-3 is concentrated in sensor 3 (mid-thickness of the specimen) after 100 days. Cluster-2 and Cluster-3 both have prominent AE energy at sensor 2 before the 66th day. Then the highest AE energy portion moves to sensor 3, whereas SSCF of Cluster-2 is much lower than Cluster-3, especially after 150 days. The SSCF for cluster Cluster-1 is negligible compared to other clusters. In the unconfined specimen, the highest SSCF is for cluster Cluster-4 (energy concentration between 0-150 kHz) at sensor 4 (attached to top reinforcement) at 66 days. However, this energy contribution decreases with time, and the SSCF of the clusters with high-frequency components (Cluster-5 and Cluster-6) increases with time. There is no obvious energy concentration in the sensor located at mid-thickness of the specimen (sensor 6), which is different from what is observed in the confined specimen. In both specimens, the SSCF declines in low-frequency signals and increases in high-frequency signals with time. This trend in signal frequency from low to high in the confined specimen is not pronounced before 150 days. The SSCF for cluster Cluster-3 is slightly greater than the Cluster-2 at 150 days. The signal frequency trend initiates primarily after 150 days in this specimen. On the other hand, in the unconfined specimen, the frequency content evolution of AE signals is obvious from an earlier stage of ASR reaction (66 days) and is more significant after 150 days.
The confined specimen has a higher extensional strain along the Z direction than the unconfined specimen (approximately 42% more at 195 days). This expansion leads to tension concentration through the thickness of the confined specimen. Since there is no confinement through the specimen thickness it is susceptible to crack formation. In the unconfined specimen, the expansion strain is more evenly distributed between the X-Y plane and the thickness. Therefore, the tension is more uniformly distributed in the entire specimen in comparison to the confined specimen. This is also observable from the AE data, where sensor 3, located at the mid-height of the confined specimen, has a larger SSCF than the unconfined specimen (e.g., 65% for the confined specimen versus 35% for the unconfined specimen at 195 days).
As mentioned previously, the frequency of AE signals progresses from low to high as the concrete ages. This may be attributable to the formation of cracking through the coarse aggregate due to ASR progression. The crack formation inside the aggregate is expected to have higher frequency components than the cement matrix and interfacial transition zone (ITZ) as mentioned by Farnam et al. [17]. The transition from low-frequency signals to the high-frequency signals in the confined specimen initiated later than for the unconfined specimen (after 150 days). However, there are different contradictory hypotheses relating to formation of cracks in concrete due to ASR [33]. For instance, osmotic pressure theory was proposed by Hanson to describe the mechanism of expansion [34]. In this theory, the cement paste surrounding to reactive aggregates acts as a semi-permeable membrane, which water solution can pass inside the region around the reactive aggregates, but alkali-silica ions are enclosed in the reactive regions. This causes osmotic pressure and alkali-silica gel swells and exerts pressure to the cement paste. This pressure leads to crack formation in the cement paste [34]. McGowan and Vivian also proposed a similar theory as osmotic theory, which transforming a solid alkali-silica layer on a reactive aggregate to a gel by absorbing moisture from the pore solution was explained as a main reason of cracking in the cement paste due to ASR [35]. Bazant and Steffens suggested that the cracking is caused in the cement paste and interfacial transition zone due to accumulation of alkali-silica gel in the interfacial transition zone and resulting gel pressure [36]. On the other hand, Dron and Brivot assumed that crack formation occurred far away from reactive aggregates due to diffusing dissolved silica away from aggregate into the pores in the cement paste [37]. Some researchers observed that ASR gel initially forms inside a reactive aggregate and causes the pressure and crack formation inside the aggregate and surrounding cement paste [12,[38][39][40][41]. Ponce and Batic [42] related the cracking pattern of concrete due to ASR to the types of reactive aggregate. ASR cracks start to form inside aggregate or in the cement matrix depending on the aggregate type [42].
In the confined specimen, the CSS rate of Cluster-3 (energy concentration 150-300 kHz), started to increase at the age of 150 days. In the unconfined specimen, the CSS rate for the cluster with the higher frequency components also increases around that time. 150 days is close to the inflection point of the volumetric strain curve, after which point expansion rates decrease. In addition, the first visible cracks were observed at the age of 150 days on the sides of the unconfined specimen, but no cracks were visible on the top surface of the unconfined specimen. Cracks could not be traced in the confined specimen on the sides due to the steel confinement frame. From the above observations, 150 days is a significant time period for ASR in the specimens, which generally agrees with trends in the AE data.

Conclusions
Acoustic emission was utilized for monitoring the activities caused by ASR in large-scale reinforced concrete specimens. The specimens resemble common nuclear power plant containments with no shear reinforcement. An agglomerative hierarchical algorithm was used to classify the AE data based on the energy-frequency dependent features to study and identify the damage mechanisms in the specimens with different stress boundary conditions. The conclusions of this study are summarized as follows: • A significant portion of the AE data (in terms of cumulative signal strength) in the confined specimen was recorded by sensor 3, which was located at mid-thickness. However, the portion of cumulative signal strength for the corresponding sensor at mid-thickness of the unconfined specimen was less in comparison to the other two sensors in that specimen. This agrees with expectations, as the confined specimen exhibited increased out-of-plane expansion in comparison to the unconfined specimen, meaning that the crack distribution is expected to be more concentrated near mid-thickness of the confined specimen than near the reinforcement layers.

•
The frequency contents of signals in the confined and unconfined specimens evolved from low to high frequency with the age of concrete although this evolution started later in the confined specimen than the unconfined specimen. Since the high-frequency AE signals have been associated to the cracking in aggregates [17], different crack mechanisms in aggregate for the confined and unconfined specimens are expectable. However, there are different contradictory ASR cracking hypotheses that have been proposed by the researchers [33][34][35][36][37][38][41][42][43].

•
There is a coincident point observed in the strain curves and the CSS of Cluster-3 in the confined specimen and Cluster-6 in the unconfined specimen. The CSS rate in terms of concrete age increases at the time (around 150 days), when the strain rate is decreasing. This point is named as inflection point of strain curve, where the curvature of strain curve changes from positive to negative. The inflection point location in terms of concrete age depends on the kinetic of ASR reaction and diffusion process. Determining the inflection point, latency and characteristic time are experimentally estimated, which are the two important modeling parameters. According to the results of AE data and clustering, the inflection point location in terms of concrete age could be estimated from the variation in CSS rate change of clustered AE data.

•
Monitoring of a structural system with acoustic emission can provide useful information regarding condition based maintenance and/or retrofit. For example, one potential time of action for treating affected structures is around the inflection point in the volumetric strain curve which can be approximated through acoustic emission data. This point coincided with observation of first visible surface cracking. After identifying the time of action, treatment methods may be implemented to mitigate the effects of ASR. Injection of lithium solution is a chemical alternative to mitigate ASR provided that enough solution penetration in the structure can be achieved. Another method is to remove moisture through coatings and sealers such as silane sealers and bituminous or elastomeric coatings. After conducting these methods, structures should be monitored for enough time to evaluate the efficiency of the method or methods.