Overturning of the façade in single-nave churches under seismic loading

The out-of-plane collapse of the façade represents one of the major threats and the most frequent cause of damages of churches due to strong earthquakes. Due to the slenderness of the façade and the lack of adequate connections to the side walls and the wooden roof, the seismic action can trigger the overturning. A detailed assessment is therefore required to judge whether or not to intervene. This paper presents an approach for the seismic assessment of the stability of the façade, through a discrete element model based on a photographic survey, with the aim of representing the actual geometry and arrangement of the stone units and their effects on the kinematics of the overturning. The collapse mechanism is simulated with both, quasi-static pushover and dynamic pulse-based analyses and the results compared to those of conventional rigid-body kinematics. The proposed approach is then applied to seven masonry churches that suffered severe damages during the 2009 L’Aquila (Italy) earthquake and the failure mode provided by the analyses is compared to the damages caused by the earthquake. The method is able to give a reliable estimate of the expected failure mechanism, taking into account the quality of the masonry and the connections to the side walls, while also providing the seismic acceleration required to trigger the motion and the ultimate displacement beyond which collapse occurs.


Introduction
Masonry churches are not designed to withstand horizontal loads and their vulnerability is considerably higher than that of ordinary masonry buildings. This conjecture is confirmed by the damage that has occurred during the seismic events that have struck Italy in the recent decades (Doglioni, F, Moretti, A, and Petrini 1994;Da Porto et al. 2012;Sorrentino et al. 2014;De Matteis et al. 2019;Penna et al. 2019). Because of their structural morphology, characterised by a prevailing longitudinal dimension, slender perimeter walls and the lack of internal diaphragms, churches rarely show a global response to seismic action, but rather a local response in which the most vulnerable structural elements detach from the rest and collapse with a rigid body motion (Lagomarsino and Podestà 2004a, b). Accordingly, the seismic behaviour of these buildings is often analysed by considering each architectural component (façade, apse, bell-tower, dome, vault, triumphal arche, etc.), named macro-element, which is characterized by a substantially autonomous structural response. For instance, the Italian guideline for monumental buildings (MiBACT 2015) suggests considering twenty-eight collapse mechanisms determined a priori. The verification of each collapse mechanism can be carried out on the basis of the kinematic theorem of limit analysis (Heyman 1966), which provides the capacity curve, and through the identification of an equivalent oscillator to estimate the corresponding seismic demand. Alternatively, the seismic demand can be estimated according to the rocking model (Housner 1963), and the different mechanisms analysed through linearization of the equation of motion, considering pulse-type ground motion (Mauro et al. 2015).
However, the construction details and quality of the masonry significantly influence the collapse mechanism and the adoption of a predefined shape of the macro-element is not always straightforward (Milani and Valente 2015a). For example, the detachment of façades is generally studied through simplified models based on the assumption that the façade behaves like a rigid body, detaching with a vertical crack from the transverse wall. However, this approach hardly predicts the actual behaviour that occurs during seismic action with regard to both, the part undergoing rocking motion and the position of the hinges (Milani and Venturini 2011). The shape and position of the cornerstones that provide the interlock between the façade and the transverse walls, in particular, plays a fundamental role in the structural response in the event of seismic action (de Felice and Mauro 2010).
A refined method to predict the out-of-plane behaviour of façades of the churches is particularly relevant, considering that during the earthquake that struck L'Aquila in 2009, the collapse of façades was the prevailing cause of damage of religious architectural heritage. The survey of damages reported that in a stock of 241 churches, 30% of the sample was damaged due to the activation of façade first-way collapse mechanism (Da Porto et al. 2012). A similar behaviour was also noted during the seismic events in central Italy in 2016, where the overturning of the façades was one of the recurring damage caused by the lack of efficient connections to the longitudinal walls. (Penna et al. 2019). Therefore, it becomes necessary to resort to a refined model that allows to study the collapse mechanism of the façade in masonry churches, by modelling the construction details that affect the structural behaviour. Various numerical methods have been proposed in these years (Casolo 2000;Milani et al. 2013;Milani and Valente 2015b;Pantò et al. 2017Pantò et al. , 2016Asteris et al. 2017;Casapulla et al. 2017;de Felice et al. 2017;Iannuzzo 2019;Lasciarrea et al. 2019;Malena et al. 2019). However, the effective masonry pattern is not taken into consideration in any of these methods.
The present work aims at proposing a new approach to assess the out-of-plane collapse mechanism of masonry buildings, based on a refined modelling of masonry walls with distinct elements (Lemos 2007;de Felice 2011) with the aim of representing the actual geometry and arrangement of the stone units and their effects on the kinematics of the local mechanism. Thanks to the careful reproduction of masonry pattern and to the capability of the DEM in simulating the motion in finite kinematics, the collapse mechanism and the crack pattern induced by the seismic action are detected, without the need to define the mechanism a priori. The effectiveness of the Discrete Element Method (DEM) to capture the seismic behaviour of masonry structures has already been proved in Literature (Azevedo et al. 2000;de Felice 2011;Al Shawa et al. 2012;Malomo et al. 2018; Mordanova and de Felice 2018; Meriggi et al. 2019) but its use for seismic assessment still needs to be settled. Accordingly, a procedure for static pushover and dynamic impulsive analyses has been developed in (Gobbin et al. 2021) and applied here to seven masonry churches located in L'Aquila (Italy) region and damaged during the 2009 earthquake. For each church, the connection of the façade to side walls is discretized into 2D distinct elements and the out-of-plane failure of the façade is modelled with either static or dynamic analyses to provide respectively, the capacity curve and the overturning domain. The comparison between the model predictions and either the earthquake damage or the analytical expressions for the rigid block, allows to assess both the reliability of the proposed approach and the limitation of the conventional rigid-body kinematics.

The conventional rigid-body analysis
As a reference for the following numerical analyses, let us consider the façade simply as a parallelepiped rigid body having thickness Fig. 1), and denote by θ the tilting of the façade towards the outside.
The horizontal seismic acceleration â st 0 that trigger the overturning mechanism and the horizontal displacement of the centre of gravity d st 0 beyond which the façade collapses are respectively: Accordingly, the capacity curve for slender walls is simply a bilinear curve with a first rigid branch up to the maximum acceleration â st 0 and a decreasing linear branch up to the maximum displacement d st 0 (Fig. 2). In the case of dynamic loading, as shown in (Housner 1963), a well-known analytical solution is obtained for a half sinusoidal wave acceleration pulse defined as: a(t) = −âsin( t + ), where is the circular frequency of the sinusoidal wave and â the amplitude acceleration, while = sin −1 ( ĝ a sin ) is a phase angle which insures that at time t = 0 the base acceleration have reached the value â st 0 required to initiate rocking. The equation of motion can be written as: (1) The condition for overturning can be obtained by imposing that = at the end of the pulse, and the resulting equation provides the amplitude â 0 of the sinusoidal pulse with circular frequency required to just overturn the façade: The limit condition provided by Eq. (4) can be summarised by the rocking domain ( Fig. 2), which expresses the amplitude of the pulse as a function of its period. Above the curve, the façade fails, while below it recovers. The horizontal asymptote corresponds to the acceleration â st 0 below which the rocking motion is not activated (Mauro et al. 2015), while the period T * such that the tangent to the overturning domain is equal to −1 separates the long periods for which the domain is mainly controlled by the pulse acceleration from the short periods controlled by the pulse duration.

Discrete element modelling
The proposed numerical model is based on the DEM (Cundall 1971) originally conceived for the study of rocks mechanics. The method is able to take into account the arrangement and interlocking of block units, by explicitly modelling the cracks and discontinuities that may occur during the motion and follows the displacement in finite kinematics, where single blocks can detach from the rest of the construction as the Fig. 2 Capacity curve (left) and rocking overturning domain (right) under sinusoidal pulse, provided by conventional rigid body analysis collapse mechanism evolves. In the present case, the blocks are considered as rigid bodies, while the contacts are modelled by a normal spring and a shear spring located at the edges of the blocks. The joint normal stiffnesses, K N , is determined by dividing the overall elastic modulus of masonry E m to the average blocks height h m . Shear stiffness, K S , is simply considered as half of the normal stiffness (Al Shawa 2012; Meriggi 2019; Mordanova 2018). However, stiffness does not significantly affect the results in the absence of second-order effects (de Felice 2011). The height of the blocks differs from case to case and also between the façade and the side wall, from a minimum of 15 cm, in the side wall of the church of S. Maria ad Cryptas, to 40 cm, for the lower part of the side wall of the church of S. Stefano. The strength of the contact is controlled by a Mohr Coulomb condition |τ|≤ σ·tan(φ), assuming a friction angle φ = 30°, while neglecting cohesion and tensile strength (Meriggi 2019, Mordanova 2018. With the assumption of rigid blocks, the equations of motion are explicitly integrated through the central difference scheme only for the centroid of each block at each cycle, resulting in a lower computational burden. Once blocks velocities are derived, displacement and forces are calculated at the interfaces on the basis of the stiffnesses at the contacts. The cycle ends with the determination of blocks forces. Even though the DEM is based on a dynamic relaxation algorithm, static convergence is possible through a fictitious damping, that dissipates residual vibrations (Lemos 2012). The DEM approach is far more time-consuming than rigid block limit analysis (Portioli 2020), but has the following main advantages: (1) the configuration is updated during the analysis in finite deformation, thus allowing the softening branch of the pushover curve to be captured; (2) the brittle behavior can be accounted for in the analysis, thus allowing to include an initial cohesion or a tensile strength that vanishes with the opening of the joint.

Discretization of masonry
Aiming at reproducing the effective connections with the lateral walls, each stone is modelled according to the actual geometry obtained from the survey of the wall face in a polygonal shape which incorporates the stone and the semi-thickness of the mortar joint. Such a refined model is used for the thickness of the façade and the adjacent part of the longitudinal wall, while the remaining part is modelled as an elastic block since it does not affect the collapse mechanism. A vertical fictitious joint with infinite tensile strength and cohesion, is used to join the blocks belonging to both the façade and the side wall, thus providing the interlocking between the two walls. The model is bidimensional and, accordingly, the blocks densities and contact stiffnesses are multiplied by their effective depth: either half of length of the façade l f , or the thickness of the transverse wall b s , as shown in Fig. 3, where the façade and the side wall heights are respectively denoted as h f and h s , their thickness as b f and b s . The length of the façade is denoted by l f , while l s represents the length of the adjacent side wall that is explicitly modelled. If there are openings in the façade, such as doors, windows or rosettes, their extent is deducted from the calculation of the actual thickness to be considered in the definition of densities and stiffnesses. Finally, the structure is laid on a rigid block that forms the base of the model.

Pushover quasi-static analysis
Following the procedure proposed in (Gobbin et al. 2021) the ascending branch of the pushover curve is simply obtained under force control, by applying the horizontal ground acceleration a g to the centroids of each block in small increments and checking for quasi-static equilibrium. The maximum acceleration a 0 st for which convergence is reached corresponds to the peak of the capacity curve. The acceleration steps are progressively decreased when approaching to a 0 st , as shown by the acceleration and displacement time-histories depicted in Figs. 4 and 5, respectively, where t p st denotes the time at which the maximum acceleration is reached, and d p st the corresponding displacement. A ground acceleration greater than a 0 st would activate the collapse mechanism. Once the mechanism is triggered by slightly increasing the ground acceleration a g over the threshold a 0 st , the softening branch is obtained under displacement control. The effective acceleration a eff of the blocks involved in the collapse mechanism is lower than the applied ground acceleration a g , due to the non-equilibrium state, and accordingly the base shear V is lower than the product of the overall mass m and the ground acceleration a g . For a given displacement of the control point, the effective acceleration can be estimated by detecting those blocks which are failing and computing the corresponding mass m * . Then, the effective acceleration of the mechanism can be estimated as the ratio of the base shear depurated by the contribution of the portion of the structure that is not failing, and the failing mass m * as follows: The detailed procedure implemented into the discrete element code is described in (Gobbin et al. 2021). By repeating this procedure in small steps and progressively decreasing the imposed ground acceleration, the whole softening branch of the capacity

Pulse-based dynamic analysis
The pulse-based analysis is carried out by applying a half sinusoidal acceleration lasting for a time period T = ∕ with amplitude a to the rigid basement in the horizontal direction and without damping: The analysis is repeated for varying amplitude and duration to detect the amplitude a 0 and period T 0 of the acceleration pulse after which the façade overturns. The failure mechanism is directly provided by the analysis and the algorithm automatically recognizes if collapse occurs or if the structure recovers: in the first case the displacement of the control point diverges, while in the second case it reaches its peak and then reduces as shown in Figure. Two different strategies are followed based on the length of the pulse. For long pulses, having period longer than T*, given a pulse duration T 0 , the amplitude is progressively increased starting from the value a 0 st up to the value a 0 (T 0 ) after which collapse takes place. For very long pulses, the dynamic analysis provides the same acceleration a 0 st given by pushover, while an acceleration lower than a 0 st would not activate the motion. Conversely, for short pulses, given an acceleration amplitude a 0 , the period is progressively increased up to the value T 0 (a 0 ) after which collapse occurs (Fig 6).
The sets of acceleration amplitude a 0 and pulse period T 0 obtained by this procedure define the overturning domain, as the frontier beyond which the structure would collapse. The numerical domain obtained by this procedure could be compared with that provided by Eq. (4) referred to the rigid block. It should be noted that, in this case, the structure is also subjected to the initial pulse ramp for values less than â st 0 and furthermore the overturning does not only occur at the end of the pulse (Makris & Konstantinidis, 2003), since when the excitation expires (i.e. for t > T), the model is let run with Time-histories of the horizontal acceleration applied in quasi-static pushover and pulsedynamic analyses horizontal acceleration equal to zero and, especially for short pulses, the displacement could increase and lead the structure to collapse well beyond the duration of the pulse. The analyses also provide the failure mechanism and the displacement capacity, i.e. the maximum displacement d 0 (T) that the structure is able to sustain without collapsing (Fig. 5). It appears, from the analyses, that the displacement capacity does not depend on the pulse period and coincide with the value d 0 st provided by pushover analysis as described in previous section. However, since small increments in the amplitude or in the period can lead the model to collapse or recovery, due to the unavoidable approximations, the maximum displacement d 0 (T) is generally a lower-bound estimate of the  effective displacement capacity, which can be more easily estimated through pushover analysis (Fig. 7c).
As will be shown in the next section, the proposed numerical approach can be directly integrated in a procedure for design and assessment of the churches against the out-ofplane overturning of the façade since: (1) The prediction of the collapse mechanism guides the design of any seismic mitigation measures; (2) The capacity curve can be directly compared with the expected seismic demand spectrum to estimate the seismic risk (Lagomarsino 2015); (3) The overturning domain can be directly compared with the rocking demand spectra to estimate the ground motion required for overturning (Makris and Konstantinidis 2003;Mauro et al. 2015).

Application to single-nave churches
Seven single-nave churches belonging to the area surrounding the city of L'Aquila ( Fig. 8), have been considered for the application of the above-described procedure for seismic assessment. These churches, after the earthquake on April 6th 2009, suffered some damages: some experienced the detachment of the façade and the triggering of the overturning mechanism, without collapsing, some others were secured with chains that prevented from failing out-of-plane. In this section the churches are briefly described, and their geometrical and mechanical properties are provided. An overall view of the churches is shown in Fig. 9, while the detailed characteristics of the connections between the façade and the side wall is shown in Fig. 10. In this latter figure the photographic survey after distortion correction is shown with superimposed recognition of masonry blocks. During the photographic survey, the dimensions of the facade and side walls of the churches are surveyed, detecting the size and arrangement of the stone elements. The discrete element model is finally constructed according to the actual masonry pattern recognized by the survey. It could be argued that the model built as explained, only considers the masonry pattern of the facing of masonry, while the internal arrangement of the stones is not accounted for. In other words, since the model is in 2D and the third dimension is simply considered by providing the thickness, the resulting model considers each block as extruded for the entire thickness of the wall. This is clearly a limitation of the proposed approach which depends on its two-dimensional character, so some care is required in the construction of the model. For example, in cases where the façade differs significantly from the wall core, it is possible to attribute to the model a reduced thickness equal to that of the face alone.
In Table 1 the geometrical characteristics of the churches are listed: l f and b f (l s and b s ) are the length and the thickness of the façade (of the lateral wall) h f and h s their heights. In Table 2 the mechanical properties are reported: γ f and γ s are the specific weight of the façade and the side walls, respectively, K Nf and K Ns the normal stiffness of the joints, K Sf and K Ss the corresponding shear stiffness, while the friction angle φ of the contacts is set equal to 30°.

Description of the churches
The Church of Sant'Andrea in Stiffe is in the municipality of San Demetrio ne'Vestini (AQ). The original structure dates to the X Century a.c., but the current structure derives from the renovation works of the XIV century, when the semicircular apse was replaced with a rectilinear one. The building has almost a rectangular plan, with a front width of about 8 m and a length of about 25 m. The façade has a horizontal crown, the masonry texture is made with rough-hewn stones with regular recurrences.
The Church of Santa Maria degli Angeli is located in the city of L'Aquila, the structure was built on the remains of an existing building. The niches along the nave show signs of subsequent infill and repeated rearrangements. The building has almost a rectangular plan, with a front width of about 10 m and a length of about 16 m. The side walls have two regular rectangular windows on each side, the masonry is in small and poorly worked stones.
The Church of Santa Maria ad Cryptas is in the municipality of Fossa (AQ), the structure was probably built between the XII and XIII Century a.c. and it preserves frescoes of great value and beauty. Inside there is the pre-existence of some vaulted parts no longer present, while the roof with wooden trusses, partially renovated in the XX century, is currently visible. The structure opens into a single hall measuring approximately 10 m by 17 m, ending with an apse not aligned with the walls of the church. The masonry is made of rough-hewn stone blocks, typical of the area around the city of L'Aquila. The Church of Santa Maria del Presepe is in the municipality of Paganica (AQ). The original structure dates around 1605 (Doglioni, F, Moretti, A, and Petrini 1994). Similarly to Sant'Andrea in Stiffe, the façade has a rectangular shape, with a front width of 15.5 m and a height of 17.6 m, covered in square ashlars, while the masonry of the side wall consists of irregular stones.
The Church of San Paolo ad Peltuinum is located in the municipality of Prata d'Ansidonia (AQ). From its first construction between the seventh and eighth centuries, the church has undergone many changes under Lombard, Benedictine, Cistercian, Renaissance and eighteenth-century influences. The church presents a cross plan with a single nave and a transept, without apse, with an overall length of about 24 m. Two different types of masonry can be identified: the lower part is made of large parallelepiped blocks and dates back to the medieval period, while the upper part, made of small irregular stones, is the result of a more recent reconstruction.
The Church of San Sisto dates back to the VIII century a.C. and is located in the city of L'Aquila. The façade has a rectangular shape with an horizontal crown. The building has a single nave structure that borders with lower buildings on the left, which were probably previously part of the church itself. The façade has a rectangular shape, with a front width of 10 m and a height of 13 m, covered in square ashlars, while the masonry of the transversal wall is in small and poorly worked stones. The Church of Santo Stefano is located in the municipality of San Pio delle Camere (AQ). There is no precise information regarding the date of construction of the Church, but an analysis of the walls shows a typical building typology in Abruzzo from the late thirteenth century to the early fourteenth century. After the second half of the XIII century transept and apse should have been realized. The building has a Latin cross plan, a single nave with a transept and apse. The hall has a rectangular plan with a width of 7 m, a length of 20 m and a maximum height of 7 m. The façade is made of square well-arranged stone blocks. The side walls, on the other hand, are made of irregularly shaped stones, with two poorly connected faces. Three different types of masonry were considered for this church, one for the façade and two for the side wall consisting of large blocks in the lower part and masonry with irregular minute stones in the upper part.

Numerical results
In this section, the numerical and analytical results are presented for the seven churches described above. As a first experimental validation of the proposed methodology, the collapse mechanisms detected through the DEM are compared to the the crack pattern caused by the 2009 earthquake for the two churches of San Paolo and S. Maria ad Cryptas, which suffered the most severe damages with the earthquake. The comparison confirms the great capability of the model to foresee the weakest joints where the detachment occurs and the resulting collapse mechanism, as shown in Fig. 11, where the mechanism provided by the analysis is superposed to the actual crack pattern.
In order to compare the Distinct Element (DE) numerical results with the conventional rigid-body (RB) analysis, in Figs. 12,13,14,15,16,17 and 18 for each church, the failure mechanism, the corresponding stability domain under sinusoidal pulses and the pushover curve provided by DEM are shown (in green line) superposed on those (in red line) provided by Eqs. (1) and (4) when considering the façade as a rigid body free to rotate outwards.
The numerical analysis of the church of Sant'Andrea in Stiffe provides the crack pattern shown in Figure A, in which the overturning mechanism involves a non-negligible part of the side wall due to the effectiveness of the connection, ensured by the modelling of the blocks arranged as in reality. The acceleration and displacement capacities are respectively: a st 0 =1.36 m/s 2 and d st 0 = 912 mm, with a substantial increase as compared to the rigid block model (Fig. 12).
The analysis of the church of San Paolo exhibits a similar increase (of about 50%) in acceleration capacity (â 0 st = 1.09 m/s 2 a 0 st = 1.71 m/s 2 ) while almost no increment was found in the displacement capacity, which remains about 1000 mm, for both the RB and DE analyses. For brief pulses (lower than 1 s), RB and DE provide almost the same results, while for longer pulses, a substantial increase in the acceleration capacity provided by DE is observed. The resulting collapse mechanism provides an explanation for this observed. Due to the quality of the masonry, which consists of large square blocks in the basement, the rotation hinge is located at approximately half of the height (Fig. 13a). Therefore, the acceleration capacity increases because of the reduced slenderness of the part affected by the overturning, while the displacement capacity, which is governed by the thickness of the façade, remains almost unchanged.
The church of San Sisto (Fig. 14) has the lower acceleration capacity as compared to other case studies presented in this paper, due to the high slenderness of the façade:   Figure C). Even though the failure mechanism shows a non-negligible part of the lateral wall involved in the collapse mechanism ( Figure  A), thanks to a good interlocking of the stones, its stabilizing effect is only marginal as a consequence of the small thickness of the lateral wall which is only about 5% of the length of the façade. Even though the failure mechanism shows a non-negligible part of the lateral wall involved in the collapse mechanism ( Figure), its stabilizing effect is not relevant as a consequence of the small thickness b s of the lateral wall.
The analysis of the church of Santo Stefano exhibits a slight underestimate (Fig. 18)

Discussion and conclusions
The 2009 earthquake caused severe damage to the architectural heritage with several churches that experienced the out-of-plane overturning of the façade. In current practice, the stability of the façade with respect to overturning, is carried out by assuming that the masonry remains monolithic and neglecting the restraint provided by the lateral walls through the corner stones. Aiming at proposing a general framework to assess the seismic capacity of collapse mechanisms of masonry structures and evaluate the efficacy of seismic prevention measures, an approach based on a refined model of the actual pattern of masonry through the Distinct Element Method (DEM) is illustrated. The seismic response of seven masonry churches located at L'Aquila has been studied through dynamic pulse-based and quasi-static pushover analyses, with both analytical and numerical approach. The influence of the quality of masonry and the beneficial effect of the interlocking of the façade with the lateral wall has been investigated and the results compared with the conventional rigid block (RB) analysis.
The proposed approach provides a reliable estimate of the expected collapse mechanism under seismic loading as shown by comparison with the actual crack pattern experienced by the churches during the earthquake. Accurate prediction of the collapse mechanism proves to be a useful tool for both, seismic assessment and eventual design of earthquake mitigation measures. In particular, the seismic assessment could be conducted either with static pushover, by simply comparing the capacity curve provided by the analysis with the expected seismic demand resulting from seismic hazard, or with dynamic pulses, by comparing the overturning domain provided by the analysis with the expected earthquake rocking demand spectra. With regard to the evaluation of seismic mitigation measures, the analysis can be easily repeated by introducing the actual characteristics of the intervention into the model.
From the sample of churches surveyed, the results in terms crack pattern are not affected by the duration of the dynamic pulse and are quite similar to those obtained with static pushover, apart from very short period pulses, which usually provide a more spread crack pattern. A more smeared crack pattern may also be the result of the assumptions of no cohesion and tensile strength in the DEM. Generally, both the displacement and acceleration capacities determined by the DEM, are greater (Sant'Andrea in Stiffe, Santa Maria degli Angeli, Santa Maria del Presepe San Sisto) or equal (San Paolo, Santa Maria ad Cryptas) than those of the conventional RB model, due to the beneficial effect of façade interlocking. The analytical RB model still represents, a reliable and effective instrument for the seismic assessment of masonry that generally provides a lower bound estimate. When the characteristics of the masonry in terms of bock shape and arrangement are such as to prevent the complete detachment of the façade from the side walls, the latter determine a significant increase in the out-of-plane capacity with the upward displacement of the tilting hinge and a corresponding reduction in the slenderness of the part subject to overturning, and an increase in the stabilising moment of the side wall involved in the mechanism. Such a contribution, which is generally neglected in current practice, would result in a greater dissipation capacity under seismic load and a reduced need for seismic mitigation interventions.
In the sample of churches under study about half of the cases show the fair reliability of the conventional RB model with deviations of less than 10% with respect to the DE model. Moreover, the DE analyses display a mechanism in which the façade remains almost monolithic during overturning, according to the assumption of RB analysis. However, in the remaining half of the cases, the RB model estimates just about 60% of the actual capacity in terms of acceleration and about 80% of the actual capacity in terms of displacement (Table 3), because of neglecting the connection between the façade and the side wall.
For brief pulses, it could be noted from the overturning domains that most of the numerical solutions result in a lower capacity with respect to the conventional RB analysis. However, as recalled in Sect. 3.4, the two analysis are not directly comparable due to the differences in the applied pulse since, contrary to what assumed by Housner (1963), for short period excitations, overturning of the façade in numerical analysis usually occurs long after the pulse action has ended. Therefore, the pulse-based numerical analyses are more severe than the analytical estimates.
In the sample of churches under study, the portion of the structure mobilised by failure only represents a fraction ranging from 15 to 35% of the overall mass, including the lateral walls. This clearly shows the local character of the collapse mechanism and the need to determine the failing mass m * as explained in Sect. 3.3, for the evaluation of the capacity curve and the estimate of the ultimate displacement. When referring to the mass of the façade, the failing mass m * usually ranges between 95 and 115%. The higher mass is mobilised in the mechanism, the greater contribution arises from the lateral walls, resulting in the increase in the out-of-plane capacity. On the contrary, if a smaller mass of the façade is mobilised, there is no direct relationship with the seismic capacity, since it could be either the result of a higher position of the rotation hinge with a corresponding increase in capacity, or a lower contribution of the side walls with a resulting lower capacity. No significant differences emerge between dynamic analyses, when varying the pulse period, nor in terms of failure mechanism neither in terms of displacement capacity, proving that the problem is mainly controlled by overall geometry of the church and the interlocking provided by the masonry, rather than on the characteristics of seismic input. The substantial stability of the solution regardless of the characteristics of the input confirms the robustness of the proposed approach and its reliability in the assessment of local collapse mechanism.
The proposed approach could appear to be cumbersome since it requires the appraisal of the actual stone masonry units connecting the façade to the side wall. In some cases, the stone units are covered by the plaster and their recognition requires instrumental inspections. However, as known to ancient builders, the corner stone plays a substantial role in ensuring the stability of the wall against overturning. Therefore, as digital survey techniques become more available, reproducing the actual masonry pattern in key parts of the building become easier and including it in a computational tool is a step towards the reliable assessment of historic masonry structures.

Data Availability
The authors confirm that the data supporting the findings of this study are available within the article.

Conflict of interest
The authors declare that there is no conflict of interest.
Code Availability Analyses are carried out though ad-hoc algorithms presented in (Gobbin et al. 2021) and developed inside the commercial software UDEC 6.0.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.