1. Introduction
Nondestructive experimental techniques and monitoring set-up are increasingly being adopted to acquire and assess the progress of dangerous structural phenomena, such as cracking and damage, and to estimate their successive evolution. The adoption of the better controlling and monitoring technique strictly depends upon the typology of the reinforced concrete or masonry structures under consideration, and on the information to be collected [
1]. For historical buildings, non-destructive techniques (NDT) can be exploited for different purposes: (1) revealing hidden structural elements, such as floor beams, or arches and piers which have been incorporated into the walls; (2) assessing the mechanical properties of masonry and mapping the inhomogeneity of the wall components (e.g., adoption of different bricks throughout the life of a building); (3) estimating the extent of cracking in damaged structures; (4) mapping voids and flaws; (5) monitoring moisture content and water rising due to capillary action; (6) surveying surface decay phenomena; and (7) assessing the mechanical and physical properties of brick or stone and mortar.
The acoustic emission (AE) technique has proved particularly effective in the assessment of structural integrity, in that it allows an assessment of the amount of energy emitted due to fracture propagation and helps derive information on the criticality of the undergoing process.
At the present time, the AE technique can be exploited during experimental tests to investigate on the damage advancement in ductile or brittle materials prior to the final collapse [
2,
3]. Moreover, this non-destructive monitoring technique is successful when studying critical phenomena and to predict the remaining lifetime and durability in full-scale structures [
4].
According to this method, it is feasible to acquire the transient elastic waves related to each stress-induced crack advancement event inside a structure or a specimen. These waves can be detected and recorded by transducers applied on the external surface of samples or structural elements. The transducers are piezoelectric sensors that commute the power of the elastic waves into voltaic signals. A proper analysis of the AE waveform parameters (peak amplitude, duration time and frequency) helps provide detailed information about the damage progression, such as the emitted energy, cracking pattern, and fracture mode [
5].
The critical conditions that anticipate the collapse can be monitored analyzing the
b-value calculated from the Gutenberg-Richter (GR) law. The GR law can fit, with basically the same accuracy, data from both earthquake distributions in seismic areas [
6] and AE technique structural monitoring [
7,
8,
9], even though two different dimensional scales are involved.
The connection between recorded waves and fracture mode depends on various factors like geometry, propagation distances, and relative orientations. Both experimental and numerical results show the crucial influence of heterogeneities in the crack propagation path. Therefore, it should be carefully considered for AE characterization of large structures, while it should not be disregarded even in small-scale sample laboratory studies to improve cracking characterization [
10,
11].
The energy emitted by the monitored structure is strictly connected to the energy detected by AE sensors. The energy dissipation in correspondence of crack nucleation and propagation in quasi-brittle materials and structures plays a fundamental role and influences their mechanical behavior throughout their entire life. In recent times, an ad hoc method has been implemented for structural monitoring by means of the AE technique according to fractal concepts. The fractal theory inherently accounts for the multiscale nature of energy dissipation and for the resulting strong size effects. This energetic approach let to introduce a useful parameter for structural damage assessment that is based on the correlation between AE monitored in the structure and the corresponding emissions detected in specimens of different sizes.
This methodology is effective in quantifying the effects of the environmental configuration on the final response, since it inherently accounts for the size of structures and the signal attenuation and distortion recorded during the AE monitoring.
The effectiveness of the assessment can be greatly improved when NDT monitoring is combined with proper numerical analysis. Numerical modeling is useful both for more efficient structural assessment and for more extensive comprehension of meso and microstructural phenomena of materials.
Analytical and numerical modeling of acoustic emissions in metals dates back to the eighties of last century [
12,
13]. The numerical simulation of progressive failure of rock and associated acoustic emissions or seismicity, depending on the scale of interest, was among others addressed by Tang [
14]. Acoustic emission monitoring in concrete structures started in the sixties of last century [
2], while numerical simulation of the acoustic emission process in concrete, according to the authors’ knowledge, are more recent [
15,
16].
The acoustic emission phenomenon can be modeled numerically following different approaches, depending on the scale and on the aspects of the problem under consideration. Dynamical models are preferred if the acoustic emission source mechanism must be considered in detail. In the literature different source hypotheses have been proposed [
17] (point-like, linear or based upon fracture mechanics description of the crack propagation), while the elastic acoustic wave propagation is affected by the attenuation, dispersion and propagation in the guiding media. If a heterogeneous concrete-like material is considered, together with the presence of another source of wave disturbance like tendons, the problem become very challenging [
18], and numerical results must be compared with additional very innovative experimental techniques like X-ray tomography [
19].
When the geometry of the sample or structure become more complex, quasi-static numerical models are usually preferred, which can provide localization in space and time of the acoustic emission events, based on a detailed description of damage propagation in the media. In this case, the signal acquired by piezoelectric sensors cannot be compared directly with numerical results, and it is necessary to calculate the acoustic emission location with some moment-tensor algorithm. Different numerical strategies can be adapted to this purpose. The so-called lattice-model allows for a detailed discretization of the microstructure, and acoustic emission events can be related to the rupture of lattice elements or bundle of elements [
20]. Alternatively, continuous discretization can be coupled with discontinuities described with the cohesive crack model [
21], and crack advancements can be directly related to the occurrence of acoustic emissions.
Numerical simulations of acoustic emissions in masonry and masonry structures are quite recent and relatively poorly investigated. For this reason, in the following, a number of applications of combined AE and numerical analysis performed by the authors on masonry and masonry structures are reviewed, from real-size structures down to model structures and meso or microstructural material modeling.
2. Historical Masonry Structures
In the last fifteen years, the authors have investigated a number of relevant existing buildings belonging to the Italian cultural heritage. In this section these structures are briefly described, in order to illustrate the combination of numerical modeling and acoustic emission monitoring.
2.1. Historical Masonry Towers of Alba
These medieval masonry towers are the tallest and most iconic constructions from the XIII
th century that are conserved in Alba up to the present time (
Figure 1) [
22,
23].
The three towers have been analyzed by three-dimensional finite element models, set up starting from Autocad® drawings (Autodesk®, San Rafael, CA, USA). The discretization mesh was obtained adopting 20-node isoparametric solid brick elements, and the analysis was performed with the finite element commercial code iDIANA (TNO DIANA BV®, Delft, The Netherlands). The discretization mesh was such that the wall thickness of the towers included five element’s nodes minimum.
The models accounted for the dimension and shape of openings and for the variation of the wall thickness at various levels. Conversely, the presence of timber floors has been neglected. Each structure was mainly subjected to its dead load. The three towers present various damage patterns, ranging from smeared cracking to dominant cracks. The acoustic emission (AE) technique was exploited to monitor the crack propagation in the most relevant structural areas of the three towers. Piezoelectric transducers applied to the inner or outer surface of the tower walls were used to acquire the elastic waves emitted as a consequence of crack propagation within the masonry volume.
From the survey it was found that the Sineo tower is not exactly vertical, and therefore the possible evolution of the eccentricity of the tower must be considered. A direct assessment of the tilt evolution based on the mechanical characteristics of the foundation is not very reliable, due to the geotechnical uncertainties of the underlying soil layers. On the other hand, a loading scenario can be analyzed, where an increasing rotation of the base is considered, together with the self-weight of the masonry and wind thrust acting on the higher lateral surface of the tower (
Figure 2). The crack pattern scenario allowed for localization of the main regions subjected to cracking in case of tilt evolution. Therefore, the location of AE sensors was optimized [
22,
23].
The mechanical response of Torre Astesiano was characterized by the presence of a dominant crack, which was uncertainly ascribed to the impact of a cannon ball during a past conflict. The structural analysis (
Figure 3) allowed it to be proven that thermal stresses due to seasonal and daily temperature oscillation could justify the opening of the crack. Two piezoelectric transducers were placed in the inner masonry layer of the tower, close to the fourth floor level, and near the tip of the main sub-vertical crack, where dynamic analysis provided the main stress intensification.
The Bonino tower was recently subjected to some restructuring works, which involved the opening of new apertures in the masonry, especially in the correspondence of the first floor level. Another cause of concern was the damage detected in the correspondence of the ornamental layer of the upper part of the tower. The structural analysis (
Figure 4) allowed for the localization of stress increases in the correspondence of the new openings, as well as for the determination of tensile stresses due to the higher stiffness of the ornamental layer, which was made with stone whose Young modulus is almost twice that of the surrounding masonry.
The AE equipment was placed in correspondence of these critical areas, in order to assess the influence on damage propagation. An evolutionary release of energy is recorded under the effect of self-weight, which can be ascribed to pseudo-creep behavior of the masonry.
The rate of propagation of the micro-cracks in the three masonry towers (
Figure 5) is correlated to the time dependence of the structural damage observed during the monitoring. The ratio of the cumulative number of AE counts recorded during the monitoring process,
N, to the number obtained at the end of the observation period,
Nd, as a function of time,
t, is equal to:
The parameter
td in Equation (1) refers to the whole structure monitoring time, while the
Ed and
Nd parameters are usually lower than the values attained at critical conditions (
Ed ≤
Emax;
Nd ≤
Nmax). It is possible to obtain an assessment about the stability condition of the structure deriving the exponent
as the linear best-fitting coefficient in the bi-logarithmic diagram (
Figure 5) where AE data from monitoring are reported.
The exponent
indicates that the damaging process is slowing down, and that the structure is evolving towards a stable condition. On the contrary, the parameter
indicates that the process is becoming unstable. Finally, the case
refers to a metastable condition,
i.e., although the phenomenon is linearly increasing over time, it could attain either instability or stability conditions. The Sineo tower analysis yielded a slope
, the Astesiano tower provided
, and the Bonino tower
(
Figure 5). During the monitoring period, these results reveal that in the Sineo and Bonino towers the damage process is stable. On the contrary, as far as the Astesiano tower is concerned, the damage is approaching a metastable condition.
2.2. Cathedrals and Vaulted Structures
A monitoring campaign and numerical analysis has been performed on the ancient Cathedral of Syracuse in Sicily [
24] (
Figure 6a). The acoustic emission (AE) technique is adopted to assess the damage pattern evolution. The localization of the propagating cracks is obtained using six synchronized acoustic emission sensors. A clear correlation between the regional seismic activity and the AE acquisition data has been obtained. In fact the AE count rate presents peaks corresponding to the main seismic events. During the observation period (
Figure 6b), the number of AE counts was
.
In order to assess the rate of the damage propagation, as given in Equation (1), the AE technique data were reported in the bilogarithmic plane to obtain the best-fitting slope, providing . The quasi-linear progression of damage over time confirms that the process in the pillar is in metastable conditions. A detailed geometrical survey of the most damaged pillar allowed for the definition of an accurate 3D model.
The geometry of each block, as well as the presence of masonry inserts, have been considered. Two main loads were considered: the dead load (of the pillar and of the surrounding structure), and a horizontal seismic load provided as horizontal ground acceleration. Due to the horizontal acceleration, cracking can take place in the pillar. A detail of crack nucleation is shown in
Figure 6d. The crack occurrence provided by the analysis agrees quite well with the crack localization provided by the AE recording. Cracking corresponds both to diffuse cracking in the continuum elements of the sandstone blocks and to opening or sliding of the discrete interfaces between blocks.
Another case was the Hospital San Giovanni in Turin (Italy) [
25], a masonry building complex initiated in 1680, under the design of the Italian architect Amedeo di Castellamonte (1610–1683) (
Figure 7). The ground floor of the complex was recently chosen to host an important fossil collection from the Regional Museum of Natural Science.
Due to this change of use, an assessment of the structural load capacity of the masonry vault beneath the first floor was necessary because the fossil collection involves a significant increase in the vault load. During the
in situ load test, we recorded the acoustic emissions from the vault, as well as displacements of the vault and strains in the steel rods. We compared the experimental data with the numerical results obtained from finite element modeling of cracking and crushing. After validation, the model allowed us to assess the ultimate load-bearing capacity of the vault. The 3D model provided a slightly better estimate of the displacements close to the abutment (
Figure 8). From the direct comparison of the cumulative AE counting distribution with the deformation of the structure, we determined that the two diagrams (
Figure 9) match each other in shape quite well. In the present case, if the time scale is plotted without rescaling, the two curves do not match identically. Nevertheless, AE are confirmed as a useful tool to measure damage, specifically in the case of very stiff vaulted structures.