Next Article in Journal
Effect of Needle Type, Number of Layers on FPAFC Composite against Low-Velocity Projectile Impact
Next Article in Special Issue
The Journal Buildings: A Bibliometric Analysis (2011–2021)
Previous Article in Journal / Special Issue
Numerical Simulation and Design Recommendations for Web Crippling Strength of Cold-Formed Steel Channels with Web Holes under Interior-One-Flange Loading at Elevated Temperatures
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Physics-Based Ground Motion Simulations for the Prediction of the Seismic Vulnerability of Masonry Building Compounds in Mirandola (Italy)

1
Faculty of Architecture and Urbanism, Politehnica University of Timişoara, 300223 Timişoara, Romania
2
Department of Engineering and Architecture, University of Trieste, 34127 Trieste, Italy
3
Department of Mathematics and Geosciences, University of Trieste, 34128 Trieste, Italy
4
Department of Structures for Engineering and Architecture, University of Naples “Federico II”, 80125 Napoli, Italy
5
Department of Civil and Environmental Engineering, University of Perugia, 06123 Perugia, Italy
*
Author to whom correspondence should be addressed.
Buildings 2021, 11(12), 667; https://doi.org/10.3390/buildings11120667
Submission received: 22 November 2021 / Revised: 13 December 2021 / Accepted: 14 December 2021 / Published: 20 December 2021
(This article belongs to the Special Issue Buildings: 10th Anniversary)

Abstract

:
The current paper aims at investigating the seismic capacity of a masonry building aggregate in the historical centre of Mirandola based on a reliable ground motion simulation procedure. The examined clustered building is composed of eleven structural units (SUs) mutually interconnected to each other, which are made of brick walls and are characterized by wooden floors poorly connected to the vertical structures. Non-linear static analyses are performed by adopting the 3Muri software to characterize the seismic capacity of both the entire aggregate and the individual SUs. In this framework, a multi-scenario physics-based approach is considered for the definition of the seismic input in terms of broadband seismic signals inclusive of source and site effects. Finally, the incidence of the seismic input variability is discussed for the prediction of the global capacity response of the case study building.

1. Introduction

The repeated occurrence of devastating seismic events has raised the technical-scientific sensitivity towards the issue of the protection of human lives and exposed buildings [1,2]. Generally, existing buildings represent the weakest vulnerable part of urban systems since they are susceptible to damage under seismic phenomena. Throughout the various construction eras, they have preserved techniques and structural details typical of a systematic design procedure conceived to resist, predominantly, gravitational actions only [3,4]. Therefore, existing unreinforced masonry (URM) buildings intrinsically present potential vulnerability factors due to multiple aspects related to the nature of the basic material, geometric complexity, construction heterogeneity and absence of appropriate seismic protection devices, that alter the entire architectural, functional and structural system [5,6].
In general, it was observed that several historical centres are characterized by an ancient building heritage having an inadequate safety level against earthquakes. In the urban centres, the existing construction typologies are frequently grouped in aggregate configurations, and the constructive matrix of structural units (SUs) is the living representation of a well-consolidated design practice that is very vulnerable to seismic actions [7,8,9,10].
Generally, the seismic capacity offered by single SUs is considerably different from that of the entire building aggregate. So, accounting for the highly non-linear response of SUs only, without considering the interactions with adjacent constructions, can inevitably lead to a wrong seismic assessment [11,12,13,14,15,16,17]. For this reason, it is mandatory to consider, preferably with simplified structural models, the intrinsic peculiarities of the SUs and their mutual interactions due to the aggregate effect to correctly foresee their vulnerability level. Different engineering computational models are developed to predict the mechanical behaviour of these historic masonry structures subjected to seismic loads. These evaluation techniques are based on sophisticated approaches, such as the discrete elements method (DEM), which is a mechanical method providing a realistic prediction of the structural seismic response. This method estimates well load capacity, stiffness, and failure hierarchy of masonry structures through highly effective mathematical computational solutions [18,19,20].
However, the behaviour of masonry structures of historical centres also depends on the seismic hazard level of the site. Principally, the type of earthquake (near-field or far-field), the type of seismic signal and related accelerations (vibratory or impulsive), and the type of fault mechanism directly influence the structural performance of this construction technology in each geographical area. Generally, the seismic hazard analyses, whether based on either a probabilistic or a deterministic approach, are based on the use of simplified correlations (i.e., ground motion prediction equations (GMPE)) for the evaluation and prediction of the ground shaking [21,22,23,24]. These formulations, obtained by data-regression analysis deriving from historical events, provide a statistical distribution of the ground motion intensity measure (e.g., spectral acceleration) as a function of multiple parameters, such as the magnitude of the earthquake, Mw, and the site-source distance, R. Furthermore, the effects of the ground motion variability (e.g., inter- or intra-event) plays an essential role in the seismic hazard analysis [25,26]. Therefore, in the first instance (inter-event phenomenon), having identified the peculiarities of the scenario (e.g., magnitude, focal mechanism and the geometry of the fault plane), the variability of the ground motion is correlated to the different possible sources rupture processes. Afterwards, having fixed the reference scenario and a specific rupture process, the spatial variability (site-to-site) of the intensity measure depends mainly on the propagation of the seismic wave through different soil layers. Overall, it is clear that the spatial variability of seismic motion is a fundamental component for the seismic risk analysis of both urban and infrastructural systems for the quantification of expected losses [27,28,29].
With such an attempt, the present work aims to provide an exhaustive evaluation of the seismic vulnerability of masonry building aggregates through a case study placed in Mirandola. To this purpose, a refined modelling procedure for the propagation of seismic waves is used. It starts from the knowledge of the seismogenic sources and the geomorphological properties of the investigated site, taking into account the rupture process of the seismic source, as well as the site effects associated with the simulated scenarios, using a web application developed by the Department of Mathematics and Geosciences of the University of Trieste. The performed ground motion simulations, being independent of the definition of a ground motion intensity measure (e.g., macroseismic intensity, PGA, Housner intensity, etc.), result to be independent of the standard approach based on the use of GMPE equations. Different kinematic rupture processes on the fault plane are parameterized by varying the slip field, nucleation point and rupture velocity to take also into consideration the source stochastic component in the distribution of the ground motion at the site(s), and the implicit consequences on the global seismic response of the SUs placed in aggregated condition.

2. Activity Framework

2.1. The Historical Centre of Mirandola

The historic compartment of Mirandola is one of the oldest urban centres on the Italian peninsula. It developed over the centuries based on a spontaneous urban planning process that gave it the first “Forma Urbis”, apparently casual, but responding to accurate settlement details prepared for the creation of roads and public spaces necessary for civil and public life, where important bureaucratic activities were carried out for the entire community. The historic centre of Mirandola is a tangible representation of the spatial configuration of a well-defined urban project starting from the early Middle Ages. The city was built as a fortified city with a well-structured and organized conformation of the historical settlements. The urban conformation of the municipality started from the construction of buildings of public interest (e.g., churches and palaces around the Castle of Pico), conceived with new architectural and functional elements to create a global planimetric structure more meticulous and responsive to the concept of the city. The city was connected to the neighbouring villages while maintaining the conformation of the main roads. In general, the buildings were arranged in a grouped configuration and bordered on the main façade, thus presenting a free front for external view and access, and an opposite front overlooking the internal pertinence area (Figure 1) [30].
Nowadays, the union of the villages and the organization of new districts have kept the original planimetric configuration of the unchanged octagonal city (Figure 2).
The historic centre is mainly characterized by masonry buildings, erected in a clustered conformation, in which the coexistence of construction heterogeneities, sometimes incongruous, has led to the manifestation of material degradation.

2.2. Main Features of the Case Study Building Compound

The historic centre of Mirandola is made up of 43 urban areas composed of buildings erected and conceived in the aggregate configuration.
The configuration of the current building stock highlights how the size of the built-up area is correlated to the evolution of the city foreseen in the Recovery Plan approved in July 2001. However, although Mirandola did not have a radial expansion, it is noted that the dimensional variation of the building aggregates is linked to the transition from the first expansion phase of the city with a quadrangular plan to the second one with an octagonal urbanized distribution [30,31]. In particular, the buildings located in the historic centre are mainly characterized by a construction uniformity of vertical structures, made of solid brick walls, with the presence of semi-rigid and deformable floors, which represents the most widespread construction typology of the historic centres of the Emilia-Romagna Region [10].
Concerning the case study aggregate, it is an existing complex masonry building with a non-regular shape composed of 11 SUs (Figure 3a). The clustered units forming the aggregate occupy different structural positions. In particular, as reported in Figure 3b, two SUs (called SU1 and SU11) occupy corner positions, nine cells (from SU2 to SU10) are placed in intermediate positions, and only SU7 occupies the heading position.
On average, the number of floors of the various SUs varies from 2 to 4 above ground, with an average inter-floor height of 2.80 m [30].
The brick masonry walls have an average thickness ranging from 0.24 m to 0.50 m. The horizontal structures are mainly made up of wooden elements with upper double planking. For these elements, a design vertical load, Fd, of 5.5 KNm−2 has been estimated, considering the fundamental load combination at the ultimate limit state (ULS) [30].
Regarding the physical conditions of the building aggregate, it is noted that SUs have common walls simply juxtaposed to each other. In addition, the pushing effects of the roofing systems could trigger the overturning of the main façades. This phenomenon is very common in historical buildings since, in the case of earthquakes, the masonry walls, poorly connected to the structure top level, are subjected to horizontal pushing actions deriving from roof, which usually activate out-of-plane mechanisms.
Concerning the material degradation, along the facades facing the main streets, even if partial detachments of the plaster are evident, they do not compromise the functional statics of the building compound.
The mechanical properties of the structural elements characterizing the case study building are deduced from the indications prescribed in Table C8.5.I of the Italian NTC18 standard [32]. Due to the absence of accurate on-site test procedures, the mechanical characteristics of the masonry are reduced by assuming a confidence factor, FC, equal to 1.35, corresponding to a LC1 (limited) knowledge level.
Specifically, for the mechanical parameters (fm, τ0, fv0), the minimum values of the intervals are reported in Table C8.5.I of the Italian standard are considered [32], while for the elastic moduli (E and G), the mean values are used. The mechanical parameters of solid brick masonry walls are summarized in Table 1.

3. Ground Motion Selection

3.1. Seismological Structural Model

To take into account the entire variability of structural behaviour from elastic to inelastic, up to global collapse [33,34], as well as the record-to-record variability (variability related to the mechanism of the seismic source, path attenuation effects and local site), ground motion records must be properly selected and, if necessary, scaled to a certain seismic intensity level to provide relatively stable estimates of the median values of the damage thresholds for the investigated structure.
For the seismic input definition, with the physics-based approach followed here, the modelling of the propagation of seismic waves for the estimation of the seismic hazard requires a prior characterization of the mechanical properties of the soil layers crossed by the waves generated by a seismogenic source. In this regard, the XeRIS web application, as proposed in [35,36], is used to compute the seismic wavefield responsible for structural damage. The modelling of the ground motion scenarios is based on the modal summation (MS) technique, as described in [37], which considers the vibrational modes of the layered medium excited by a seismic source, providing a rapid computational procedure for an adequate simulation of the ground motion in the case of a far-field event. Therefore, in the case of near-field events, the discrete wave number (DWN) technique is used according to the theory proposed in [38]. Nevertheless, these techniques allow for a very efficient calculation of synthetic seismograms in laterally homogeneous layered models, taking into account anelastic attenuation. Thus, the proposed site-specific geological model has been defined starting from the soil structure provided by the Italian Accelerometric Archive, ITACA, [39] considering a soil class C, then adopting the mechanical layers parameters (thickness, density, P and S waves velocity, and attenuation) discussed in [40], which represent the geological setting of the investigated area.

3.2. Fault Scenarios

The calculation of a fault-based ground-shaking scenario allows modelling, at a specific site, the ground motion caused by an earthquake with an extended fault. The computational technique used is based on the theory proposed by [41], where the seismogenic sources are modelled as a set of sub-sources, with a specific scaling source spectra law, also to predict the effects related to the kinematic rupture process (i.e., directivity).
Moreover, some of the kinematic parameters are made by varying the nucleation point, rupture velocity and slip field, and also take into consideration the source stochastic component [42,43].
Four different sources (Mirandola (MIR), Finale Emilia (FE), Veronese (VR) and Ferrara (FE)) are selected from the Web-GIS application DISS (Database of Individual Seismogenic Sources) [44] to compute the seismic input at the selected site. The sources have rectangular geometry of length, L, and width, W, with a variable distance between 1.9 km up to 44 km about the reference site, identified as MR (Lon (°): 11.06 and Lat (°): 44.88—geographical position of the study aggregate). In particular, Mirandola’s source descends solely on the evidence of the recent tectonic activity of the Arc of Ferrara [45]. The 29 May 2012 earthquake in the Emilia-Romagna region of Italy, which was the second main shock of the Emilian sequence, activated this seismogenic source, producing a detectable uplift of the buried Mirandola anticline. The final geometry of the fault adopted is the one validated in the study proposed by [40]. The Finale Emilia seismic source was formed after the occurrence of the 20th May 2012 earthquake, which was the first mainshock of the Emilian sequence giving rise to a detectable uplift of the associated buried anticline. The Veronese source was generated following the 3rd January 1117 Verona earthquake, which was perhaps the strongest event that occurred in the Po Valley. The data collected in both historical and instrumental catalogues have shown that the seismicity of the Po Valley is concentrated mainly along the foothills and the buried thrust fronts of the northern Apennines and the southern Alps [46]. Finally, concerning the existence of the Ferrara’s source, it dates back to the earthquake that occurred on 17th November 1570 based on data relating to both the recent tectonic activity of the Ferrara Arch [45], and from the subsoil geological structure [47]. Thus, based on the above considerations, some of the characteristics of the four considered seismic sources are reported in Table 2.
Furthermore, a map of the location of the case study building (MR) to source configuration has been presented in Figure 4, where the dark blue line identifies the top of the fault projection. The light blue circle (positioned in the centre of the fault projection) is the fault reference point, and the orange area is associated with the area where the position of the nucleation point of the rupture (red dot) is allowed to vary.
Subsequently, starting from the definition of the above-introduced seismic sources, the modelled seismic scenarios consider the moment magnitude, the position of the nucleation point, and the mean rapture velocity, Vs. The simulation proposed is based on the Monte-Carlo technique for the quantification of the distribution of the final slip, considering 100 statistical realizations and assuming the random variability of the nucleation point of coordinates (X, Y) (the maximum allowable variation is ±1/3 of the fault length (X direction) and ±1/3 of the fault width (Y direction)).
The simulation of the source rupture processes has been performed using the PULSYN06 algorithm based on the theory presented in [41,48].
For a sake of example, the fault rupture scenarios generated by four realizations (i.e., #1, #50, #90, #100) out of 100 for the Mirandola (MIR) source, are shown in Figure 5. For each simulated rupture scenario, a 2D final slip map function is derived. The sub-source grid is made up of 17 × 9 crosses in which the red dot represents the nucleation point, the dark blue indicates the high slip of the fault, and the white contours (isochrones) represent the time-domain slip rupture propagation of a non-stationary random process.

3.3. Maximum Credible Seismic Input (MCSI)

For the assessment of the seismic vulnerability of the considered case study, an appropriate calculation procedure is presented to predict the expected hazard for the site of interest.
For earthquake engineering purposes, the maximum credible seismic input (MCSI) represents a reliable estimation of the expected ground shaking level for a specific site, independently of the occurrence rate of earthquakes that could affect the investigated area. In particular, the proposed methodology [37,43] does not use empirical equations, such as GMPE, to derive the intensity measure (e.g., PGA or SA, PGV, PGD engineering parameters), but takes into consideration the seismic scenarios generated by independent seismic sources providing realistic ground motions in the time domain (either as a response spectrum or as a set of seismograms). Nevertheless, the use of source spectra [41] introduces a stochastic element in the proposed methodology; therefore, to define the MCSI spectrum, the uncertainties correlated to the hazard estimation must be considered.
Physics-based synthetic engineering parameters can be calculated through knowledge of the earthquake generation process and the propagation of seismic waves in an inelastic soil profile.
In this regard, the structural and topographical heterogeneities, as well as the influence of the source rupture process (slip distribution and rupturing velocity of the fault), were evaluated using the XeRIS web application, as proposed in [35].
Regarding the physical definition of the MCSI response spectrum, it is described in [29,42] where, for each seismogenic source, n-scenarios (in terms of magnitude, epicentre distance and focal mechanism of the earthquake) have been considered, the obtained n-spectral acceleration values (SA) have been compared, and the maximum median has been selected, as depicted in Figure 6.

4. A Mechanical Vulnerability Assessment Model

4.1. Global Behaviour of the Case Study Aggregate

The case study compound has been numerically analysed using the 3Muri software founded on the Frame Macro-Elements (FME) theory (Figure 7) [49]. This methodological approach assumes that masonry walls are considered as a set of one-dimensional macro-elements (columns, beams, and nodes) mutually interconnected. The resistance criteria of the deformable elements have been assumed based on the requirements of EN 1998-3 [50], which establish maximum threshold values of 0.4% and 0.6% for shear failure and flexural collapse, respectively. The analyses were carried out along the two main global directions of the building, X and Y.
The above-introduced simulations were interrupted when the shear strength decay was equal to 20% of the maximum shear capacity, as reported in the NTC18 standard [32].
To this purpose, two different distributions of seismic forces and four combinations (±X, ±Y) have been considered, neglecting the influence of accidental eccentricities in the calculation procedure. The results in terms of capacity curves are presented in Figure 8.
From the gotten results, it was observed how the aggregate provides a better structural capacity in the transverse direction Y, rather than in the orthogonal one X. In particular, the ultimate displacements achieved in the positive Y direction (Y+) are approximately two times higher than the corresponding ones detected in the perpendicular X-direction referred to the uniform distribution of the forces.
Conversely, in the negative X direction (X), an increase of the displacements was observed, compared to the Y direction.
In particular, the static condition (X 1st Mode) has shown an ultimate displacement, du, equal to −5.18 cm, which is approximately 61% greater than the analogous displacement attained in the direction Y (du = −3.18 cm). Regarding the maximum base shear threshold, it was noted that in the Y direction, the maximum strength is approximately two times greater than the corresponding value in direction X. This discrepancy in terms of global behaviour (e.g., displacements and maximum shear thresholds) in the two analyses directions is due to the irregular planimetric conformation of the entire building conjunctly with a higher presence of masonry resistant area.
Finally, for a global view of the vulnerability threshold and expected damage concerning the examined building, two indices were analysed: the vulnerability index (VI) and the damage index (DI). With this intention, the vulnerability index, VI, is intended as the ratio between the acceleration capacity, PGACi, intended as the peak acceleration that determines the achievement of the ultimate limit state (ULS) condition based on the return periods, Tr, and the corresponding acceleration demand, PGADi, related to the spectral coordinates of the site where the construction is located. In the specific case, the 3 Muri software [49] automatically estimates the maximum spectral acceleration, starting from the definition of the elastic response spectrum (Tr = 475 year) for the site of interest (Mirandola) having assumed a soil class C, as reported in the Italian Accelerometric Archive, ITACA [39]. Consequently, the damage index has been evaluated, taking into account the mathematical formulation proposed in [51]:
DI = μ max - 1 μ d - 1 = d max - d y d u - d y
where µmax and µd are the required and available ductility associated with the structural system, respectively, which can be expressed in terms of displacements; dmax represents the required displacement, evaluated for Ti < TC, based on the area equivalence criterion; and dy and du are the yielding and ultimate displacements associated with the SDoF system, respectively.
The numerator of the above-mentioned formulation indicates how the maximum displacement, dmax, deviates from the elastic one, while the denominator provides indications about the type of failure mechanisms (e.g., when the ultimate displacement, du, is comparable to dy, the fragile mechanism is activated, vice versa, ductile ones). Thus, congruently to what has been introduced, the synthesis of the results has been depicted in Figure 9.
First, the achieved results presented in Figure 9a have shown that the worst-case scenario is attained for the direction X+ (uniform), which corresponds to a C/D ratio equal to 0.67. Similarly, in the direction Y (Y+ 1st Mode), the expected vulnerability level is 0.40 (e.g., medium-low level).
Concerning the damage index DI, as reported in Figure 9b, it has been observed that the value of 0.57 is reached for a uniform distribution of forces in the X+ direction, since the difference between the maximum displacement demand, dmax, and the corresponding yielding displacement, dy, is influenced by the achievement of a premature plastic displacement threshold.
Subsequently, the vulnerability indexes have been associated with the estimated damage index, DI, and then correlated to the damage thresholds provided by the European Macroseismic Scale, EMS-98 [52,53], as reported in Table 3.
In Figure 10, the validation between the real damages detected after the Emilia-Romagna earthquake and the corresponding ones derived from the mechanical analysis procedure evaluated at the ultimate displacement, du, referring to the ULS (ultimate limit state) combination, has been provided.
The damages detected [30] are characterized by the widespread distribution of cracks in the main façade that affect the vertical structures.
More specifically, there is a propensity to damage identified mainly in the spandrel elements, since they have no tensile strength (instead guaranteed by a steel tie-rod or reinforced concrete ring beam) towards the combined effects induced by both bending and shear.
Congruently to what was discussed in the above-introduced comparison, it is possible to appreciate how the simulated structural model provides a global damage level not dissimilar to that really detected in the case study aggregate.

4.2. Capacity Response of Building Structural Units

This section analyses the seismic behaviour of the individual structural units in the aggregate configuration using a numerical procedure based on the capacity curve reconstruction process derived from that of the whole aggregate. This methodological approach allows estimation of the propensity to damage of the SUs, considering their interaction in the aggregate configuration and, therefore, accounting for their contribution in terms of both stiffness and mass. Thus, after assessing the seismic response of the whole aggregate along with the main analysis directions, the behaviours of the reference cells included in the aggregate have been evaluated. In particular, the clustered SUs response has been achieved by adopting a step-by-step simulation process by considering as base shear, Vbi, the analytical sum of the base reactions of the piers of the reference case study unit and, as displacement, the average value of the displacements of the selected monitoring nodes, Dmi, at the top floor. The analytical formulations, adopted for the definition of the MDoF capacity curves, are presented according to the study proposed in [5,54,55]:
{ V b i = j = 1 N | V j | D m i = | j = 1 N D j | | N |
having assumed |Vj| as the absolute value of the base shear force associated with the j-th node at the i-th step, Dj as the absolute value of the top displacement associated with the j-th node at the i-th step, and N as the total number of the nodes at the base of the masonry piers.
The resulting curves in both analysis directions, X and Y, after a bi-linearization process, have been plotted in Figure 11.
From the acquired results, it appears that the behaviour of the SUs in the two analysis directions is quite heterogeneous. In particular, the head units have been denoted with black-grey, the corner units in red, and a black scale colour has been adopted for the other intermediate units. The capacity parameters expressed in terms of forces and displacements for each of the analysed SUs are presented in Table 4.
From the obtained results, it has been observed that, in the X direction, the head unit (SU7) has higher stiffness and strength than the other building units which occupy intermediate and corner positions, respectively.
This circumstance has been revealed since the corner and intermediate units are affected by the beneficial effect of mutual structural confinement, while the head unit is constrained only on one side, thus resulting in greater susceptibility to displacements. In terms of displacements, it was noted that the SU7 (head) provides displacements comparable to the SU5 (intermediate), which has a higher mass due to the presence of contiguous SUs.
Moreover, in the Y direction, it has been observed that the SU7 provides a higher contribution in terms of displacement capacity than the other cells due to the interconnected confinement effect provided by the adjacent SUs. Furthermore, by comparing the maximum shear threshold of the SU7 in the two directions, it can be seen that, in the Y direction, the structural cell has an increased shear strength of about 44% compared to the corresponding threshold achieved in the orthogonal direction.

4.3. Analysis of the Capacity Parameters

To interpret the different global responses of the examined structural units, the main capacity parameters (stiffness, ductility, shear strength, mass and vibration period) have been plotted and correlated to each other, as reported in Figure 12.
First of all, the vibration period of the single structural units analysed, taking into account their contribution in terms of mass and stiffness, has been appropriately compared with the simplified formulation derived from the NTC18 standard [32] (see Figure 12a). From the comparison, it has been noted that in the two analysis directions, the periods associated with the single SUs are heterogeneous with a higher incidence in the X direction, where they exceed the corresponding period provided by the NTC18 standard [32].
However, it is worth pointing out that the NTC18 standard [32] provides a simplified formulation for estimating the vibration period of buildings considered as isolated structures, without considering the aggregate effect. Thus, it is noted that the SUs 1, 2, 3, 6b, 8, 9 and 10 have a low seismic mass-to-stiffness ratio, presenting a reduced vibration period compared to the other structural cells examined in the same analysis direction. Conversely, in Y direction, the vibration period is, on average, 24% lower compared to the analogue values examined in the orthogonal direction.
This circumstance denotes that, in the transverse direction, the stiffness associated with the single SUs is on average 79% higher than that obtained in the X direction, as depicted in Figure 12b. As far as the displacement ductility is concerned, it has been correlated to the wall resistant area (see Figure 12c) to consider the incidence of the structural resistant area (WR.A) on the expected ductility (µ). From the results achieved, it was observed that in the X direction, the point cloud is quite homogeneous for 2.0 < µ < 7.0, which corresponds to a wall resistant area WR.A enclosed in the range between 2.60 m2 and 21.0 m2. On the contrary, due to the higher number of wall panels, in the Y direction, the wall resistant area is within 7.40 m2 and 32.0 m2 and reduces the corresponding displacement.
Finally, in Figure 12d, the correlation between the maximum shear force and the total volume of SUs is shown. The analysis of the results in both analysis directions has revealed how the linear regression function presents a direct proportionality between the monitored parameters. The regression coefficient (R2) is robust, assuming values greater than 0.5 in both examined cases, taking into account all the building-to-building variability in terms of mass, stiffness, and in-elevation interaction.

4.4. Vulnerability Assessment for the Bilinear Elasto-Plastic SUs Models

Following the above-introduced sections, the generation of MCSI spectra was implemented according to the soil structural model and fault scenario selected [56]. In particular, the MCSI has been defined in terms of the response spectrum to accomplish the engineering seismic verification of the case study masonry compound.
With this intent, according to [29,42], the MCSI response spectrum, taking into account the stochastic nature of the algorithm implemented by the XeRIS tool [35], is evaluated at the 95th percentile, based on the number of realizations of the fault rupture process considered. Therefore, by defining the maximum credible seismic input, the maximum contribution in terms of seismic hazard is offered by the Mirandola source (MIR, see Table 4) due to a reduced site-source distance. In these circumstances, the source spectrum derived from the rapture process of the Mirandola source represents the maximum expected MCSI for the generated scenario, enveloping, for periods of interest, the spectra generated by the other sources considered. The derived spectra for both the 50th and 95th percentile (median and maximum), are presented in Figure 13a, with their regularized shape according to the Italian seismic code as illustrated in Figure 13b. Finally, in Figure 13c, the design response spectra for the Mirandola site [32] have been presented, assuming soil class C and Tr = 475 years.
Consequently, the N2 method [57,58,59] has been adopted for the evaluation of the seismic performance of the single SUs in aggregate conditions by adopting the above-introduced MCSI spectrum. The method combines non-linear static analysis and the response spectrum approach to derive the ADRS (acceleration displacement response spectrum) domain, as reported in Figure 14. Moreover, the seismic verification is carried out by estimating the vulnerability index, appropriately evaluated as the ratio between the seismic demand, Di, and the relative capacity, Ci, provided by the studied structures.
The results presented have shown how the expected vulnerability is heterogeneous in both analysis directions, depending on the structural position of SUs, as discussed in Section 3.1 and Section 3.2. For the sake of example, a detailed application of the N2 method to the head unit is depicted in Figure 15, where also the ADRS obtained using the source spectra is compared with the corresponding spectrum derived from the NTC18 standard [32].
From the obtained results, it is noted that, in the transverse direction, there is an increase of both strength (55%) and stiffness (40%) concerning the longitudinal direction response. Subsequently, the results obtained by applying the ADRS spectrum, derived from the NTC18 standard [32], have revealed how the verification in the X direction provides a D/C ratio equal to 0.83 while, considering MCSI at the 50th percentile, the mechanical vulnerability, D/C, is equal to 0.96, with an increase of 14%. In the Y direction, it is noted as D/C|NTC18 = 0.7, while the corresponding value obtained by MCSI is equal to 0.9, with an increase of 22%. Moreover, the ADRS spectrum achieved from the MCSI evaluated at the 95th percentile provides unsafe conditions, since the seismic checks are not satisfied (D/C > 1) in both analysis directions.
In conclusion, the N2 method verifications implemented through the proposed approach to derive seismic spectra have provided results on the safe side, in comparison to those provided by the standard code.

4.5. Fragility Assessment

In this section, the mean fragility curves of the different SUs have been derived. The fragility curves define the probability that a given damage state, dSi, is equalled or exceeded by a specific intensity measure. The adopted methodology considers as intensity measure the spectral displacement, Sd, according to the following Equation (3) [10,55,60]:
P ( d d Si | S d ) = Φ [ 1 β ln ( S d S d ¯ ) ]
where Φ [•] stands for the cumulative normal distribution, while the denominator and β are the median and the standard deviation of the corresponding normal distribution, respectively. The damage thresholds have been assumed according to the theory proposed by [3,60] as a function of the structure yielding and ultimate displacements.
Moreover, to reduce the number of uncertainties deriving from the deterministic approach, a single value has been assumed for the standard deviation.
Specifically, the standard deviation has been set as equal to the displacement ductility, µ, multiplied for a fixed correlation coefficient equal to 0.45. Thus, the damage states and the standard deviations have been identified, as reported in Table 5.
The expected damage levels are in line with the assumption proposed in the framework of the Risk-UE project [61], which is generally used in many seismic risk studies on European building stocks in prone cities. Therefore, the fragility functions have been derived, taking into consideration the average displacements of the above-mentioned case study SUs, as depicted in Figure 16.
The obtained results clearly show how the corner structural unit provides damage probabilities higher than those of the other SUs occupying intermediate and head positions. In the specific case, the corner unit, despite having a number of floors less than the contiguous structural cells, has a reduced lateral constraint guaranteed by the adjacent structural units, so that the torsional phenomena limit the displacement capacity and, consequently, the expected damage is much more pronounced. Regarding the torsional effects, they are strictly correlated by the addition of structural mass rates offered by the contiguous SUs, resulting in a non-uniform distribution of the resistant elements in the plane.
Therefore, the more pronounced eccentricity between the centre of gravity and centre of stiffness triggered the premature failure of the masonry wall panels, confirming the damage level depicted in Section 4.1.

5. Conclusions

The research work analysed the seismic vulnerability of an unreinforced masonry building aggregate located in the historic centre of Mirandola using a physics-based simulation approach.
The first step aims to characterize the expected hazard for the site of interest, based on the physics-based ground motion simulations. To this end, the web application XeRIS has been adopted. The starting point was the definition of the soil structure profile to model the propagation of the seismic waves for the estimation of the seismic hazard. In this regard, parametric studies were conducted both on the geological structural model and on the configurations of the selected seismic sources.
The proposed structural model contains the specific geological properties of the investigated area. Subsequently, four different sources (Mirandola (MIR), Finale Emilia (FE), Veronese (VR) and, Ferrara (FE)), were selected to derive the MCSI source spectrum by adopting a stochastic seismic scenario that considers all the possible variability related to the moment magnitude, the position of the nucleation point, and the mean slip rapture velocity.
The second step concerned the evaluation of the global seismic behaviour of the examined masonry compound through non-linear static analyses. The different global behaviour (e.g., displacements and maximum shear thresholds) in the two analysis directions is mainly due to the irregular planimetric configuration of the entire building conjunctly to the higher presence of resistant wall areas. Subsequently, concerning the global response, two synthetic parameters, namely the vulnerability index (VI) and the damage index (DI), were analysed. To this end, the vulnerability index, VI, was mechanically estimated directly by extrapolating the results from the 3Muri software, which once defined the elastic response spectrum for the site of interest. From the results achieved, it was shown that the X+ (uniform) direction provides a PGAC/PGAD ratio equal to 0.70 while, in the transverse direction (Y+ 1st Mode), the vulnerability stands at 0.40 (e.g., medium-low level). As far as the damage index (DI) is concerned, it has been used for establishing a failure hierarchy based on the analyses performed. In particular, it has been found that the selected parameter varied in the range between 0.1 and 0.6, which mainly corresponds to a ductile (bending) crisis of the structural elements and the corresponding damage level according to the EMS-98 scale.
Consequently, the seismic behaviour of the individual structural units in the aggregate configuration was analysed using a numerical procedure based on the reconstruction process of the capacity curves derived from those of the entire aggregate. This methodological approach allows estimation of the damage propensity of the single SUs, considering their interaction in the aggregate, so taking into account the contribution given by SUs adjacent to the investigated one in terms of both stiffness and mass.
The performed analyses provided important insight into the correlation between shear force and ductility. It has been observed that, by correlating the base shear force and the total volume of the single SUs, there is a direct proportionality among the monitoring parameters. However, the ductility tends to decrease as the total resistant area of walls in Y direction increases up to 35 m2.
Consecutively, the seismic verification was implemented using the bilinear elasto-plastic curves of SUs. To this purpose, the ADRS domain was adopted for the MCSI spectrum (derived at 50th and 9th percentile). In the first instance, the acquired results showed that the expected vulnerability is heterogeneous in the analysis directions, providing unsatisfactory results using the ADRS at the 95th percentile (worst-case scenario). Furthermore, to propose a suitable example regarding the performance of the SUs towards the seismic action, the investigated head unit was considered as a reference case for conducting the seismic verification in terms of demand and capacity displacements. In this regard, the ADRS obtained using the source spectra were compared with the corresponding spectrum derived according to the Italian design code. From the obtained results, it was noticed that, in the transverse direction, an increase of both strength (55%) and stiffness (40%) concerning the longitudinal direction was found.
The seismic verification in the longitudinal direction, obtained by applying the ADRS domain derived from the Italian NTC18 standard, provided a D/C ratio equal to 0.83, while, considering MCSI at the 50th percentile, the mechanical vulnerability, D/C, was equal to 0.96, with an increase of 14%. Moreover, in the Y direction, it was seen as D/C|NTC18 = 0.7, while the corresponding value obtained by MCSI was equal to 0.9, with an increase of 22%.
This aspect demonstrated how the displacement demands achieved by adopting the simulated scenario provide conforming results concerning the prediction of the code spectrum vulnerability procedure.
Finally, the mean fragility curves were derived for the different SUs investigated. The results showed how the corner structural unit provides higher damage probabilities than the other SUs occupying intermediate and head positions. Firstly, this outcome was due to the reduced lateral constraint guaranteed by the adjacent structural units, and, secondly, by the induced effect of the torsional phenomena strictly correlated to the addition of structural mass rates offered by the contiguous SUs, resulting in a non-uniform distribution of the resistant elements in the plane.
Overall, based on the results presented, this research topic provided complete information concerning the adoption of a specific procedure based on ground motion simulations for the seismic assessment of URM buildings.
For this reason, this study represents an initial first step that can provide clear insights into the causes of collapse observed in masonry structures after destructive earthquakes. In this framework, future developments should be oriented to the simulation of other seismic scenarios characterized by different fault mechanisms, epicentre distances, hypocentre depths, and ground amplification, to have a greater consistency on the effects caused by a possible earthquake on a larger building dataset.

Author Contributions

Conceptualization, N.C., F.R., M.F. and A.F.; methodology, N.C., F.R., M.F. and A.F.; software, N.C. F.R. and M.F.; validation, N.C., F.R., M.F. and A.F.; formal analysis, N.C.; investigation, N.C., F.R., M.F. and A.F.; resources, A.F.; data curation, A.F. and G.M.; writing—original draft preparation, N.C., F.R. and M.F.; writing—review and editing, N.C., F.R., A.F. and G.M; visualization, N.C.; supervision, F.R. and A.F.; project administration, A.F. and G.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data deriving from the current study can be provided to the readers based upon their explicit request.

Acknowledgments

The first author thanks Franco Vaccari for his help in the seismic computations, which have been partially made using the XeRiS facilities (http://www.xeris.it/WebApp/, accessed on 10 October 2021). This work was developed under the financial support of the Italian Civil Protection Department, within the ReLUIS-DPC 2019–2021 research project, which is gratefully acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dolce, M.; Prota, A.; Borzi, B.; da Porto, F.; Lagomarsino, S.; Magenes, G.; Moroni, C.; Penna, A.; Polese, M.; Speranza, E.; et al. Seismic risk assessment of residential buildings in Italy. Bull. Earthq. Eng. 2021, 19, 1–34. [Google Scholar] [CrossRef]
  2. da Porto, F.; Donà, M.; Rosti, A.; Rota, M.; Lagomarsino, S.; Cattari, S.; Borzi, B.; Onida, M.; De Gregorio, D.; Perelli, F.L.; et al. Comparative analysis of the fragility curves for Italian residential masonry and RC buildings. Bull. Earthq. Eng. 2021, 19, 1–44. [Google Scholar] [CrossRef]
  3. Lagomarsino, S.; Cattari, S.; Ottonelli, D. The heuristic vulnerability model: Fragility curves for masonry buildings. Bull. Earthq. Eng. 2021, 19, 3129–3163. [Google Scholar] [CrossRef]
  4. Masi, A.; Lagomarsino, S.; Dolce, M.; Manfredi, V.; Ottonelli, D. Towards the updated italian seismic risk assessment: Ex-posure and vulnerability modelling. Bull. Earthq. Eng. 2021, 19, 3253–3286. [Google Scholar] [CrossRef]
  5. Formisano, A.; Chieffo, N.; Vaiano, G. Seismic Vulnerability assessment and strengthening interventions of structural units of a typical clustered masonry building in the campania region of Italy. GeoHazards 2021, 2, 101–119. [Google Scholar] [CrossRef]
  6. Clementi, G.-E.F.; Formisano, A.; Milani, G.; Ubertini, F. Structural health monitoring of architectural heritage: From the past to the future advances. Int. J. Arch. Herit. 2021, 15, 1–4. [Google Scholar] [CrossRef]
  7. Grillanda, N.; Valente, M.; Milani, G. Adaptive NURBS based local failure analyses of retrofitted masonry aggregates. In AIP Conference Proceedings; AIP Publishing LLC: Melville, NY, USA, 2021; Volume 2343, p. 080003. [Google Scholar]
  8. Bernardini, C. Identification of Minimum Unit of Analsis for Seismic Performance Assessment of Masonry Buildings in Ag-gregate. Ph.D. Thesis, University of Florence, Florence, Italy, 2019. [Google Scholar]
  9. Battaglia, L.; Ferreira, T.M.; Lourenço, P.B. Seismic fragility assessment of masonry building aggregates: A case study in the old city Centre of Seixal, Portugal. Earthq. Eng. Struct. Dyn. 2021, 50, 1358–1377. [Google Scholar] [CrossRef]
  10. Chieffo, N.; Formisano, A.; Mochi, G.; Mosoarca, M. Seismic vulnerability assessment and simplified empirical formulation for predicting the vibration periods of structural units in aggregate configuration. Geosciences 2021, 11, 287. [Google Scholar] [CrossRef]
  11. Cima, V.; Tomei, V.; Grande, E.; Imbimbo, M. Seismic vulnerability of italian historical towns: The case of borgo san rocco. In Advances in Science and Technology; Trans Tech Publications Ltd.: Freienbach, Switzerland, 2021; Volume 109, pp. 1–7. [Google Scholar]
  12. Jiménez, B.; Saloustros, S.; Pelà, L. Seismic vulnerability index method for hybrid timber–masonry structures. Numerical calibration and application to the city of Valparaíso, Chile. J. Build. Eng. 2021, 44, 103185. [Google Scholar] [CrossRef]
  13. Angiolilli, M.; Lagomarsino, S.; Cattari, S.; Degli Abbati, S. Seismic fragility assessment of existing masonry buildings in ag-gregate. Eng. Struct. 2021, 247, 113218. [Google Scholar] [CrossRef]
  14. Leggieri, V.; Ruggieri, S.; Zagari, G.; Uva, G. Appraising seismic vulnerability of masonry aggregates through an automated mechanical-typological approach. Autom. Constr. 2021, 132, 103972. [Google Scholar] [CrossRef]
  15. Ruggieri, S.; Cardellicchio, A.; Leggieri, V.; Uva, G. Machine-learning based vulnerability analysis of existing buildings. Autom. Constr. 2021, 132, 103936. [Google Scholar] [CrossRef]
  16. Rapone, D.; Brando, G.; Spacone, E.; De Matteis, G. Seismic vulnerability assessment of historic centers: Description of a pre-dictive method and application to the case study of scanno (Abruzzi, Italy). Int. J. Archit. Herit. 2018, 12, 7–8. [Google Scholar] [CrossRef]
  17. Roca, P.; Cervera, M.; Gariup, G.; Pelà, L. Structural analysis of masonry historical constructions. Classical and advanced approaches. Arch. Comput. Methods Eng. 2010, 17, 299–325. [Google Scholar] [CrossRef] [Green Version]
  18. Ferrante, A.; Giordano, E.; Clementi, F.; Milani, G.; Formisano, A. FE vs. DE modeling for the nonlinear dynamics of a historic church in central Italy. Geosciences 2021, 11, 189. [Google Scholar] [CrossRef]
  19. Ferrante, A.; Schiavoni, M.; Bianconi, F.; Milani, G.; Clementi, F. Influence of stereotomy on discrete approaches applied to an ancient church in Muccia, Italy. J. Eng. Mech. 2021, 147, 04021103. [Google Scholar] [CrossRef]
  20. Salachoris, G.P.; Magagnini, E.; Clementi, F. Mechanical characterization of “Scaglia Rossa” stone masonry through exper-imental and numerical analyses. Constr. Build. Mater. 2021, 303, 124572. [Google Scholar] [CrossRef]
  21. Anderson, J.G.; Cotton, F.; Bindi, D. A ground motion based procedure to identify the earthquakes that are the most relevant for probabilistic seismic hazard analysis. Earthq. Spectra 2021, 37, 762–788. [Google Scholar] [CrossRef]
  22. Rodríguez-Castellanos, A.; Ruiz, S.E.; Bojórquez, E.; Reyes-Salazar, A. Influence of spectral acceleration correlation models on conditional mean spectra and probabilistic seismic hazard analysis. Earthq. Eng. Struct. Dyn. 2021, 50, 309–328. [Google Scholar] [CrossRef]
  23. Ordaz, M.; Salgado-Gálvez, M.A.; Giraldo, S. R-CRISIS: 35 years of continuous developments and improvements for proba-bilistic seismic hazard analysis. Bull. Earthq. Eng. 2021, 19, 2797–2816. [Google Scholar] [CrossRef]
  24. Dávalos, H.; Miranda, E. A ground motion prediction model for average spectral acceleration. J. Earthq. Eng. 2021, 25, 319–342. [Google Scholar] [CrossRef]
  25. Romanelli, F.; Altin, G.; Indirli, M. Spreading NDSHA application from Italy to other areas. In Earthquakes and Sustainable Infrastructure; Elsevier BV: Amsterdam, The Netherlands, 2022; pp. 175–194. [Google Scholar]
  26. Brandmayr, E.; Vaccari, F.; Panza, G.F. Neo-deterministic seismic hazard assessment of Corsica-Sardinia block. Rend. Lincei. Sci. Fis. E Nat. 2021, 1–10. [Google Scholar] [CrossRef]
  27. D’Ayala, D.; Meslem, A.; Vamvatsikos, D.; Porter, K.; Rossetto, T.; Silva, V. Guidelines for Analytical Vulnerability Assessment-Low/Mid-Rise. GEM Tech. Rep. 2014. Available online: https://storage.globalquakemodel.org/resources/publications/ (accessed on 20 October 2021).
  28. Barnaba, M. Physics-Based Seismic Input for Engineering Design and Risk Assessment. Ph.D. Thesis, University of Trieste, Trieste, Italy, 22 October 2021. [Google Scholar]
  29. Panza, G.F.; Vaccari, F.; Romanelli, F.; Amadio, C.; Fasan, M.; Magrin, A. A seismological and engineering perspective on the 2016 Central Italy earthquakes. Int. J. Earthq. Impact Eng. 2016, 1, 395. [Google Scholar] [CrossRef] [Green Version]
  30. Ricchiuto, F. Strategie Valutative Della Vulnerabilità Sismica Degli Aggregati Edilizi in Relazione Alla Caratterizzazione Costruttiva: Risultanze Delle Analisi Sugli Aggregati Dell’Intero Centro Storico Di Mirandola; University of Bologna: Bologna, Italy, 2013. (In Italian) [Google Scholar]
  31. Formisano, A.; Mochi, G.; Chieffo, N. Empirical and mechanical analysis methods for seismic vulnerability assessment of clustered buildings of historical centres: A case study. In Proceedings of the 8th International Conference on Computational Methods in Structural Dynamic Analysis and Its Application to Performance-Based Earthquake Engineering, Athens, Greece, 28–30 June 2021. [Google Scholar] [CrossRef]
  32. Ministry of Infrastructure and Transport. Updating Technical Standards for Construction; Official Gazette: Rome, Italy, 2018. (In Italian) [Google Scholar]
  33. Vamvatsikos, D.; Cornell, C.A. The Incremental Dynamic Analysis and Its Application to Performance-Based Earthquake Engineering. In Proceedings of the 12th European Conference on Earthquake Engineering, London, UK, 9–13 September 2002; p. 479. [Google Scholar]
  34. Lantada, N.; Pujades, L.G.; Barbat, A.H. Vulnerability index and capacity spectrum based methods for urban seismic risk evaluation. A comparison. Nat. Hazards 2008, 51, 501–524. [Google Scholar] [CrossRef]
  35. Vaccari, F.; Magrin, A. A user-friendly approach to NDSHA computations. In Earthquakes and Sustainable Infrastructure; Elsevier BV: Amsterdam, The Netherlands, 2022; pp. 215–237. ISBN 978-0-12-823503-4. [Google Scholar]
  36. Magrin, A.; Gusev, A.A.; Romanelli, F.; Vaccari, F.; Panza, G.F. Broadband NDSHA computations and earthquake ground motion observations for the Italian territory. Int. J. Earthq. Impact Eng. 2016, 1, 131. [Google Scholar] [CrossRef]
  37. Panza, G.F.; Romanelli, F.; Vaccari, F. Seismic wave propagation in laterally heterogeneous anelastic media: Theory and ap-plications to seismic zonation. In Advances in Geophysics; Dmowska, R., Saltzman, B., Eds.; Elsevier: Amsterdam, The Netherlands, 2001; Volume 43, pp. 1–95. [Google Scholar]
  38. Pavlov, V.M. Matrix impedance in the problem of the calculation of synthetic seismograms for a layered-homogeneous isotropic elastic medium. Izv. Phys. Solid Earth 2009, 45, 850–860. [Google Scholar] [CrossRef]
  39. Italian Accelerometric Archive v3.1. Available online: http://itaca.mi.ingv.it/ItacaNet_31/#/home (accessed on 6 October 2021).
  40. Paolucci, R.; Mazzieri, I.; Smerzini, C. Anatomy of strong ground motion: Near-source records and three-dimensional physics-based numerical simulations of theMw6.0 2012 May 29 Po Plain earthquake, Italy. Geophys. J. Int. 2015, 203, 2001–2020. [Google Scholar] [CrossRef]
  41. Gusev, A.A. Broadband kinematic stochastic simulation of an earthquake source: A refined procedure for application in seismic hazard studies. Pure Appl. Geophys. 2010, 168, 155–200. [Google Scholar] [CrossRef]
  42. Fasan, M.; Magrin, A.; Amadio, C.; Panza, G.F.; Romanelli, F.; Vaccari, F.; Noè, S. A possible revision of the current seismic design process to overcome the limitations of standard probabilistic seismic input definition. In Proceedings of the 16th World Conference on Earthquake, 16WCEE 2017, Santiago, Chile, 9–13 January 2017; p. 12. [Google Scholar]
  43. Rugarli, P.; Amadio, C.; Peresan, A.; Fasan, M.; Vaccari, F.; Magrin, A.; Romanelli, F.; Panza, G.F. Neo-Deterministic Scenario-Earthquake Accelerograms and Spectra: A NDSHA approach to seismic analysis. In Engineering Dynamics and Vibrations: Recent Developments; Jia, J., Paik, J.K., Eds.; CRC Press: Boca Raton, FL, USA, 2018; pp. 187–241. [Google Scholar]
  44. Database of Individual Seismogenic Sources, V. 3.0. Available online: http://diss.rm.ingv.it/diss/ (accessed on 6 October 2021).
  45. Vannoli, P.; Burrato, P.; Valensise, G. The seismotectonics of the Po Plain (Northern Italy): Tectonic diversity in a blind faulting domain. Pure Appl. Geophys. 2015, 172, 1105–1142. [Google Scholar] [CrossRef] [Green Version]
  46. Guidoboni, E.; Ferrari, G.; Mariotti, D.; Comastri, A.; Tarabusi, G.; Valensise, G. Catalogue of Strong Earthquakes in Italy (461 B.C.-1997) and Mediterranean Area (760 B.C.-1500). CFTI4Med. 2007. Available online: http://storing.ingv.it/cfti/cfti5/html/info_EN.html (accessed on 20 October 2021).
  47. Mantovani, E.; Viti, M.; Babbucci, D.; Tamburelli, C.; Vannucchi, A.; Baglione, M.; D’Intinosante, V.; Cenni, N. Stima Aggiornata della Pericolosità Sismica in Toscana e aree Circostanti; Universitas Studiorum: Milan, Italy, 2018; ISBN 978-88-336-9014-8. [Google Scholar]
  48. De Natale, G.; Petrazzuoli, S.; Romanelli, F.; Troise, C.; Vaccari, F.; Somma, R.; Peresan, A.; Panza, G.F. Seismic risk mitigation at Ischia island (Naples, Southern Italy): An innovative approach to mitigate catastrophic scenarios. Eng. Geol. 2019, 261, 105285. [Google Scholar] [CrossRef]
  49. STA DATA. 3Muri-Seismic Calculation of Masonry Structures; STA DATA Srl: Turin, Italy. Available online: https://www.3muri.com/documenti/brochure/en/3Muri10.9.0_ENG.pdf (accessed on 20 October 2021).
  50. EN 1998-3. Eurocode 8: Design of structures for earthquake resistance. In Part 3: Assessment and Retrofitting of Buildings; CEN: Bruxelles, Belgium, 2004; pp. 1–97. [Google Scholar]
  51. Chieffo, N. Seismic Vulnerability Assessment of Historical Masonry Building under Near-Field Earthquake. Ph.D. Thesis, Politehnica University of Timisoara, Romania, Balkans, 27 August 2021. [Google Scholar]
  52. Grünthal, G. The European macroseismic scale 1992: (up-dated MSK-scale). In Cahiers du Centre Européen de Géodyna-Mique et de Séismologie 7; Centre Européen de Géodynamique et de Séismologie: Luxembourg, 1998; p. 99. [Google Scholar]
  53. Hadzima-Nyarko, M.; Morić, D.; Pavić, G.; Mišetić, V. Spectral functions of damage index (DI) for masonry buildings with flexible floors. Teh. Vjesn. 2018, 25, 181–187. [Google Scholar] [CrossRef] [Green Version]
  54. Chieffo, N.; Formisano, A. Comparative seismic assessment methods for masonry building aggregates: A case study. Front. Built Environ. 2019, 5, 123. [Google Scholar] [CrossRef]
  55. Chieffo, N.; Clementi, F.; Formisano, A.; Lenci, S. Comparative fragility methods for seismic assessment of masonry buildings located in Muccia (Italy). J. Build. Eng. 2019, 25, 100813. [Google Scholar] [CrossRef]
  56. Hassan, H.M.; Fasan, M.; Sayed, M.A.; Romanelli, F.; ElGabry, M.N.; Vaccari, F.; Hamed, A. Site-specific ground motion modeling for a historical Cairo site as a step towards computation of seismic input at cultural heritage sites. Eng. Geol. 2020, 268, 105524. [Google Scholar] [CrossRef]
  57. Fajfar, P.; Fischinger, M. N2-A Method for Non-Linear Seismic Analysis of Regular Buildings. In Proceedings of the Ninth World Conference in Earthquake Engineering, Tokyo, Japan, 2–9 August 1988; Volume 5, pp. 111–116. [Google Scholar]
  58. Mosoarca, M.; Onescu, I.; Onescu, E.; Azap, B.; Chieffo, N.; Szitar-Sirbu, M. Seismic vulnerability assessment for the historical areas of the Timisoara city, Romania. Eng. Fail. Anal. 2019, 101, 86–112. [Google Scholar] [CrossRef]
  59. Mosoarca, M.; Onescu, I.; Onescu, E.; Anastasiadis, A. Seismic vulnerability assessment methodology for historic masonry buildings in the near-field areas. Eng. Fail. Anal. 2020, 115, 104662. [Google Scholar] [CrossRef]
  60. Cattari, S.; Lagomarsino, S.; Ottonelli, D. Fragility curves for masonry buildings from empirical and analytical models. In Proceedings of the 2nd Conference on Earthquake Engineering and Seismology, Istanbul, Turkey, 25–29 August 2014; pp. 25–29. [Google Scholar]
  61. Mouroux, P.; Le Brun, B. RISK-UE project: An advanced approach to earthquake risk scenarios with application to different european towns. In Assessing and Managing Earthquake Risk; Springer: Berlin/Heidelberg, Germany, 2008; pp. 479–508. [Google Scholar]
Figure 1. Historical urban conformation of the city of Mirandola [30]. On the left, a bird-eye view of the city, in which the highlighted areas indicate the urban expansion process, is illustrated. On the right, the urban development started from the “ancient village” (XIII Century) to the complete configuration of the citadel (XV Century).
Figure 1. Historical urban conformation of the city of Mirandola [30]. On the left, a bird-eye view of the city, in which the highlighted areas indicate the urban expansion process, is illustrated. On the right, the urban development started from the “ancient village” (XIII Century) to the complete configuration of the citadel (XV Century).
Buildings 11 00667 g001
Figure 2. Octagonal compartment and urban blocks plan of the city (a,b) [30]; (c,d) some building typologies in the city centre.
Figure 2. Octagonal compartment and urban blocks plan of the city (a,b) [30]; (c,d) some building typologies in the city centre.
Buildings 11 00667 g002
Figure 3. Localization of the case study aggregate (a) and identification of clustered structural units (SUs) (b).
Figure 3. Localization of the case study aggregate (a) and identification of clustered structural units (SUs) (b).
Buildings 11 00667 g003
Figure 4. The investigated source-site configurations; (a) MIR: Mw= 6.0, edi = 1.9 km; (b) FE: Mw= 6.1, edi = 20.7 km; (d) VR: Mw= 6.7, edi = 43.8 km; (c) FR: Mw= 5.5, edi = 34.6 km.
Figure 4. The investigated source-site configurations; (a) MIR: Mw= 6.0, edi = 1.9 km; (b) FE: Mw= 6.1, edi = 20.7 km; (d) VR: Mw= 6.7, edi = 43.8 km; (c) FR: Mw= 5.5, edi = 34.6 km.
Buildings 11 00667 g004
Figure 5. Four examples of fault-slip rupture processes were generated for the Mirandola source, varying the position of the nucleation point and the number of realization processes; (a) MIR: Mw = 6.0, #1; (b) MIR: Mw = 6.0, #50; (c) MIR: Mw = 6.0, #90; (d) MIR: Mw = 6.0, #100.
Figure 5. Four examples of fault-slip rupture processes were generated for the Mirandola source, varying the position of the nucleation point and the number of realization processes; (a) MIR: Mw = 6.0, #1; (b) MIR: Mw = 6.0, #50; (c) MIR: Mw = 6.0, #90; (d) MIR: Mw = 6.0, #100.
Buildings 11 00667 g005
Figure 6. Simulated maximum credible seismic input (MCSI) for the selected sources. The grey band identifies the MCSI, as defined in [43], which is controlled by three sources: curves 1, 2 and 3 are the median spectra obtained from one hundred realizations of the rupture process for the sources MIR, FE and VR, respectively.
Figure 6. Simulated maximum credible seismic input (MCSI) for the selected sources. The grey band identifies the MCSI, as defined in [43], which is controlled by three sources: curves 1, 2 and 3 are the median spectra obtained from one hundred realizations of the rupture process for the sources MIR, FE and VR, respectively.
Buildings 11 00667 g006
Figure 7. Representation of the case study aggregate: solid (a) and macro-element (b) models.
Figure 7. Representation of the case study aggregate: solid (a) and macro-element (b) models.
Buildings 11 00667 g007
Figure 8. Capacity curves in (a) ±X and (b) ±Y directions.
Figure 8. Capacity curves in (a) ±X and (b) ±Y directions.
Buildings 11 00667 g008
Figure 9. Estimated mechanical vulnerability index (VI) (a) and damage index (DI) (b).
Figure 9. Estimated mechanical vulnerability index (VI) (a) and damage index (DI) (b).
Buildings 11 00667 g009
Figure 10. Validation of the mechanical model by comparison between real damages and numerical ones (a,b).
Figure 10. Validation of the mechanical model by comparison between real damages and numerical ones (a,b).
Buildings 11 00667 g010
Figure 11. Bilinear curves for the selected SUs in X (a) and Y (b) analysis directions.
Figure 11. Bilinear curves for the selected SUs in X (a) and Y (b) analysis directions.
Buildings 11 00667 g011
Figure 12. Capacity of the SUs in terms of (a) vibration period, (b) stiffness and mass, (c) ductility displacements vs. wall resistant area, and (d) base shear vs. total volume.
Figure 12. Capacity of the SUs in terms of (a) vibration period, (b) stiffness and mass, (c) ductility displacements vs. wall resistant area, and (d) base shear vs. total volume.
Buildings 11 00667 g012aBuildings 11 00667 g012b
Figure 13. MCSI spectra at (a) the 50th and 95th percentile, (b) regularization response spectra, and (c) design response spectrum for the Mirandola site.
Figure 13. MCSI spectra at (a) the 50th and 95th percentile, (b) regularization response spectra, and (c) design response spectrum for the Mirandola site.
Buildings 11 00667 g013
Figure 14. Capacity curves of SUs in the ADRS plane; (a) Corner SUs, X dir.; (b) Corner SUs, Y dir.; (c) Intermediate SUs, X dir.; (d) Intermediate SUs, Y dir.; (e) Head SU, X dir.; (f) Head SU, Y dir.
Figure 14. Capacity curves of SUs in the ADRS plane; (a) Corner SUs, X dir.; (b) Corner SUs, Y dir.; (c) Intermediate SUs, X dir.; (d) Intermediate SUs, Y dir.; (e) Head SU, X dir.; (f) Head SU, Y dir.
Buildings 11 00667 g014
Figure 15. Seismic checks of the head SU in the ADRS plane in (a) X and (b) Y directions, considering the worst mechanical analysis case; (a) Head SU, X dir.; (b) Head SU, Y dir.
Figure 15. Seismic checks of the head SU in the ADRS plane in (a) X and (b) Y directions, considering the worst mechanical analysis case; (a) Head SU, X dir.; (b) Head SU, Y dir.
Buildings 11 00667 g015
Figure 16. Fragility curves for the examined structural units in X (ae) and Y (bf) directions; (a) Head SUs; (b) Head SUs; (c) Intermediate SUs; (d) Intermediate SUs; (e) Corner SUs; (f) Corner SUs.
Figure 16. Fragility curves for the examined structural units in X (ae) and Y (bf) directions; (a) Head SUs; (b) Head SUs; (c) Intermediate SUs; (d) Intermediate SUs; (e) Corner SUs; (f) Corner SUs.
Buildings 11 00667 g016aBuildings 11 00667 g016b
Table 1. Mechanical characteristics of masonry material for inspected SUs [32].
Table 1. Mechanical characteristics of masonry material for inspected SUs [32].
Mechanical Properties
Average compressive strengthfm2.60 Nmm−2
Average shear strengthfv00.13 Nmm−2
Young modulusE1500 Nmm−2
Shear elasticity modulusG600 Nmm−2
Dead weightW18 KNm−3
Table 2. Some characteristics of the four selected seismogenic sources [40,44]. The acronyms are related to: edi = epicentral distance; depth = focal depth of the earthquake occurred; strike = defines the orientation of the fault counter-clockwise till the North; dip = represents the inclination of the fault plane; rake = the direction the hanging wall moves during rupture measured relative to the fault strike; sre = represents the angle formed between the focal point and the site.
Table 2. Some characteristics of the four selected seismogenic sources [40,44]. The acronyms are related to: edi = epicentral distance; depth = focal depth of the earthquake occurred; strike = defines the orientation of the fault counter-clockwise till the North; dip = represents the inclination of the fault plane; rake = the direction the hanging wall moves during rupture measured relative to the fault strike; sre = represents the angle formed between the focal point and the site.
IDSourceMwL (km)W (km)edi (km)Depth (km)Strike (°)Dip (°)Rake (°)Sre (°)
1MIR6.022121.908.995609069
2FE6.1106.420.76.21154390195
3VR6.7301443.87.5248409074
4FR5.55.14.034.62.9885090179
Table 3. Correlation between DI and European Macroseismic Scale (EMS-98) damage levels.
Table 3. Correlation between DI and European Macroseismic Scale (EMS-98) damage levels.
DIDamage DescriptionDamage Grade EMS-98
0 < DI < 0.3NullD0
0.3 < DI < 0.6SlightD1
0.6 < DI < 0.8ModerateD2
0.8 < DI < 1.0Near-CollapseD3
1.0 < DICollapseD4-D5
Force DistributionPGAC/PGADDIEMS-98
[-][-]
X+ (Uniform)0.70.6D2
X- (1st Mode)0.50.4D1
X- (Uniform)0.40.2D0
X- (1st Mode)0.40.3D1
Y+ (Uniform)0.20.1D0
Y+ (1st Mode)0.40.2D0
Y- (Uniform)0.30.1D0
Y- (1st Mode)0.20.0D0
Table 4. Capacity parameters for the investigated SUs.
Table 4. Capacity parameters for the investigated SUs.
PositionSUs IDDir. XDir. Y
dy [cm]du [cm]Fmax [daN]dy [cm]du [cm]Fmax [daN]
Corner10.20.949,1200.10.280,529
Intermediate20.31.739,1750.20.260,121
Intermediate30.21.327,6250.20.251,357
Intermediate40.31.838,5780.40.588,310
Intermediate50.21.756,0300.30.4190,361
Intermediate60.21.348,1220.30.4219,884
Intermediate6b0.11.316,1330.10.314,598
Head70.31.6112,8310.41.1251,410
Intermediate80.10.276940.10.222,418
Intermediate90.10.266150.20.232,505
Intermediate100.10.232,7230.30.360,920
Corner110.10.284,4760.20.441,891
Table 5. Definition of damage levels and the corresponding standard deviation.
Table 5. Definition of damage levels and the corresponding standard deviation.
Damage Levels DSi
D10.7·dySlight
D21.5·dyModerate
D30.5·(dy + du)Near Collapse
D4-D5duCollapse
Standard deviation β
0.45·ln(µ)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chieffo, N.; Fasan, M.; Romanelli, F.; Formisano, A.; Mochi, G. Physics-Based Ground Motion Simulations for the Prediction of the Seismic Vulnerability of Masonry Building Compounds in Mirandola (Italy). Buildings 2021, 11, 667. https://doi.org/10.3390/buildings11120667

AMA Style

Chieffo N, Fasan M, Romanelli F, Formisano A, Mochi G. Physics-Based Ground Motion Simulations for the Prediction of the Seismic Vulnerability of Masonry Building Compounds in Mirandola (Italy). Buildings. 2021; 11(12):667. https://doi.org/10.3390/buildings11120667

Chicago/Turabian Style

Chieffo, Nicola, Marco Fasan, Fabio Romanelli, Antonio Formisano, and Giovanni Mochi. 2021. "Physics-Based Ground Motion Simulations for the Prediction of the Seismic Vulnerability of Masonry Building Compounds in Mirandola (Italy)" Buildings 11, no. 12: 667. https://doi.org/10.3390/buildings11120667

APA Style

Chieffo, N., Fasan, M., Romanelli, F., Formisano, A., & Mochi, G. (2021). Physics-Based Ground Motion Simulations for the Prediction of the Seismic Vulnerability of Masonry Building Compounds in Mirandola (Italy). Buildings, 11(12), 667. https://doi.org/10.3390/buildings11120667

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop