Mechanical-Analytical Soil-Dependent Fragility Curves of Existing RC Frames with Column-Driven Failures

: In seismic risk estimation, among the different types of fragility curves used (judgement-based, mechanical, empirical/observational, hybrid), the mechanical ones have the twofold advantage of allowing a better control over the basic parameters and of representing a validation test of the consistency of empirical/observational ones. In this study, fragility curves of RC frames with column-driven failures are obtained from a simpliﬁed analytical pushover method implemented in a simple spreadsheet, thus allowing the user to perform a large number of analyses. More importantly, the proposed method introduces the concept that Limit States at the structural level are obtained consequent to the attainment of the same Limit States at the local level, in the columns’ sections. This avoids using additional criteria, such as interstorey drift thresholds. This simple analytical model allows for rapid development of fragility curves, for any Limit State, of different building typologies identiﬁed by a set of global quantities (number of storeys, story heights, number of spans and span lengths) and by a set of local quantities (element sizes, reinforcement, and material properties). It also allows for a straightforward treatment of the inﬂuence of the soil class on the fragility curves parameters, which is another critical issue addressed in this work that helps when interpreting some literature results using empirical/observational methods.


Introduction
In the recent past, the topic of seismic risk analysis and mitigation has attracted the research community worldwide. Large-scale territorial studies have emphasized the need for developing risk analysis tools to be utilized by both civil protection agencies and insurance companies [1,2]. One of the main components of seismic risk is the vulnerability of structures. Dunand and Guneguen [3] showed that, even in regions with moderate seismic hazard, a high seismic risk can be found if the seismic vulnerability of the construction stock is high.
A convenient and widely adopted method for defining seismic vulnerability is the use of Fragility Curves (FCs). Such curves provide the probability of exceeding a Limit State (LS) or a Damage Level (DL) as function of a seismic intensity measure.
Different classifications have been proposed of the methods adopted to derive FCs [4,5]. A meaningful classification identifies four different groups of procedures, based on the main source of information: (a) judgement-based, (b) mechanical, (c) empirical/observational, and (d) hybrid.
Judgement-based FCs are derived from statistical treatment of estimates provided by experts about the average DL that various types of structures may undergo when subjected

Simplified Analytical Model of A 2D Frame
To describe the global geometry of a 2D frame, representative of a building typology, only a small yet significant set of parameters is needed, such as: number of storeys, storey heights, number of spans, and span lengths. Permanent and variable loads are generally known from the building category. As far as materials are concerned, they can be deduced from the codes enforced at the time of construction of the building. The same codes can be used to deduce the element sizes and the amount of longitudinal and transverse reinforcement, either using the minimum amounts thereby defined or through a simulated design according to the then-enforced rules. The simplified analytical model of a 2D frame rests on its interpretation as a series system of storeys, which in turn are considered as parallel systems of columns, which in turn are considered as series systems of flexural and shear-resisting mechanisms. Notwithstanding the simplified system logic, such models represent a reasonable compromise between typological representativeness and computational efficiency. The method is conceived to be applied to existing buildings designed for gravity loads, where failures are generally column-driven, since columns are often weaker than the framing beams. Additionally, due consideration is given to failure of beam-column joints by limiting the capacity of the framing columns accordingly.
One main issue arising when developing mechanical FCs is the global Damage Levels (DL) definition. Here, some correspondence is established with the analytically defined local Limit States (LS) at both the sectional and the element level. This step is essential in that DLs, which have been defined mainly to facilitate and homogenize observational data, are usually given in descriptive terms and, as such, suffer from being qualitative and subjective. On the other hand, LSs are quantitively and objectively stated and are therefore analytically usable.
In this study, the 2D frame response at the global level is described in terms of base-shear vs. top-displacement by a piece-wise trilinear curve, where three objectively measured LSs are identified: Damage limitation Limit State (DLS), Life safety Limit State (LLS) and Collapse prevention Limit State (CLS).
These three global LSs are directly obtained from the corresponding LSs of each single story, which in turn are derived from the LSs of each single column, which, again, are derived from the local LSs of the sections at the column's end plastic hinges, with criteria that will be explained in the following. Following this approach, any global LS is directly identified based on the attainment of the corresponding LS at the story level, at the column level, and ultimately at the section level. Once LS criteria are established at the section level, no additional LS criteria are needed at the column level, at the story level, and at the global level. This ensures consistency between local LSs and global LSs and avoids introducing a different metric for the global LSs, such as, for example, interstorey drifts.

Frame Section Capacity at x = D, L, C
The sectional moment-curvature diagram can be effectively described by a trilinear curve passing through three Limit States: DLS, corresponding to yielding of the tensile steel bars, LLS, corresponding to crushing of the concrete cover, and CLS, corresponding to crushing of the concrete core ( Figure 1). The coordinates of these three points can be described by closed-form equations for curvature capacity φ xLS and moment capacity M xLS , where x = D, L, C.
For rectangular sections with symmetrical reinforcement, these equations are given in Table 1 (further elaborated from [17]), with the notation in Figure 2, where d is the section effective depth, b is the section width, d c is the section effective depth after concrete cover spalling, A s is the tension area of the longitudinal reinforcement, A sw is the area of the transverse reinforcement, and s is the transverse reinforcement spacing.
Where ε y is the steel yield strain, ε cu is the concrete strain, f c is the concrete strength, µ s is the longitudinal reinforcement mechanical ratio, and n S is the normalized axial load, the latter two given as, respectively: Moreover, σ 2 is the confining stress normalized with respect to the characteristic concrete strength f ck , and η f is the concrete strength increase due to confinement, given as: For rectangular sections with symmetrical reinforcement, these equations are given in Table 1 (further elaborated from [17]), with the notation in Figure 1, where is the section effective depth, is the section width, is the section effective depth after concrete cover spalling, is the tension area of the longitudinal reinforcement, is the area of the transverse reinforcement, and is the transverse reinforcement spacing.

Damage Limitation Limit State (DLS)
Moment capacity For rectangular sections with symmetrical reinforcement, these equations are given in Table 1 (further elaborated from [17]), with the notation in Figure 1, where is the section effective depth, is the section width, is the section effective depth after concrete cover spalling, is the tension area of the longitudinal reinforcement, is the area of the transverse reinforcement, and is the transverse reinforcement spacing.   If needed, the sectional curvature ductility at CLS can be calculated as:

Moment capacity
2.3. Column Capacity at x = D, L, C A column is considered as a series system of flexural plastic hinges, column shear, and joints shear.
The flexural force-displacement diagram of a column is obtained by considering, both, the rotations of the two end plastic hinges, which follow the above-described sectional moment-curvature diagram, and the elastic deformation of the column between them. The so-obtained purely flexural response of the column is then compared to its shear capacity and the joints', to obtain the actual capacity of the column.
The column shear and displacement corresponding to the three Limit States x = D, L, C are obtained from equilibrium and compatibility as: where the elastic stiffness is: where M xLS is the flexural capacity at xLS, H is the element height, φ xLS is the curvature capacity at xLS, E cm is the concrete mean elastic modulus, I is the section moment of inertia, and L p is the plastic hinge length, given, for example, as: The column shear capacity V R is computed as: where ω w is the transverse reinforcement mechanical ratio: with A w being the transverse reinforcement area, while ν = 0.5 is the concrete compressive strength reduction factor, and α c is a coefficient depending on the normalized axial load n S , as follows: The joint shear capacity V jR is computed as: where d c is in Figure 2, b j is the joint effective width, and η = α j 1 − f c 250 with f c in MPa and α j = 0.6 for internal joints and α j = 0.48 for external joints.

Story Capacity at x = D, L, C
A story is considered as a parallel system of j = 1 . . . m columns. The i-th interstorey drift capacity d xLS,i and the corresponding resultant force V xLS,i at x = D, L, C is determined when the first column attains the corresponding Limit State: where u xLS,j is the j-th column displacement at xLS and v j (d xLS,i ) is the j-th column force at d xLS,i . The i-th interstorey secant stiffness is: 2.5. Frame Capacity at x = D, L, C The frame can be represented as a series system (a so-called "stick model") of i = 1, . . . , n storeys, where each storey mass m i is located at height z i ( Figure 3) and where the connecting elements follow the respective trilinear story behavior, obtained from the trilinear columns behaviors, in turn obtained from the trilinear sectional behaviors, as explained in the previous sections.
where is in Figure 1, is the joint effective width, and = α 1 − with in MPa and α = 0.6 for internal joints and α = 0.48 for external joints.

Story Capacity at = , ,
A story is considered as a parallel system of = 1 … columns. The -th interstorey drift capacity , and the corresponding resultant force , at = , , is determined when the first column attains the corresponding Limit State: where , is the -th column displacement at and , is the -th column force at , . The -th interstorey secant stiffness is:

Frame Capacity at = , ,
The frame can be represented as a series system (a so-called "stick model") of = 1 … storeys, where each storey mass is located at height ( Figure 2) and where the connecting elements follow the respective trilinear story behavior, obtained from the trilinear columns behaviors, in turn obtained from the trilinear sectional behaviors, as explained in the previous sections. In such system, the lateral forces are applied along the frame height according to the shape ψ i , which can be taken as either one of the following (height-proportional or mass-proportional): where it should be noticed that ∑ ψ i = 1 for both. From the applied force shape ψ i , the interstorey shear shape τ i is found as: The elastic interstorey displacement shape can be found as: Buildings 2021, 11, 278 7 of 20 with k DLS,i being the i-th story elastic stiffness at DLS.
The frame displacement capacity d xLS and resulting base shear V xLS at x = D, L, C are determined when the first storey attains the corresponding Limit State: where d xLS,i is the i-th interstorey drift in Equation (12) at xLS and V xLS,i is the interstorey shear in Equation (13), taken at the interstorey satisfying Equation (18), denoted with imin.

Simplified Modal Analysis
The first modal shape φ can be taken equal to the displacement shape. In "stick models" with evenly distributed masses along the height, this is a reasonable approximation that yields negligible errors (lower than 3% for both the eigenvalue ratio and the eigenvector distance norm). Therefore: where it should be noticed that, by replacing Equation (17), the maximum value at the top (i = n) is naturally equal to 1. The participation factor Γ and the fundamental period T are then found as usual as: It should be noticed that M and K are diagonal and banded, respectively, as follows (posing for conciseness k i = k DLS,i ): This allows rewriting Equations (21) and (22) using only summations, thus simplifying their use, for example, in spreadsheets or hand calculations: where the effective mass m * and the effective stiffness k * are:

Equivalent SDOF System
The "stick model" of the frame is then transformed into an equivalent single degree of freedom (SDOF) system. The top displacement (d * ) and the base shear (V * ) of the equivalent SDOF system are obtained from Equations (18) and (19), as follows: Having elastic stiffness:

Bilinearization
The SDOF system capacity curve, which has a trilinear shape through DLS, LLS, and CLS, is bilinearized to determine its equivalent yield point, which is found as: The elastic period of such a bilinearized system is:

Displacement Demands at x = D, L, C
The displacement demands on the frame system pertaining to x = D, L, C are obtained from the elastic displacement spectral ordinates S De,xLS (T * ) computed at T * on the displacement spectra pertaining to x = D, L, C, as: where T C,xLS is the period at the onset of the descending branch of the acceleration spectrum, and: where S e,xLS (T * ) are the elastic spectral ordinates computed at T * on the acceleration spectra pertaining to x = D, L, C. It is important to notice that the spectral shape changes based on the soil class.

Definition of Fragility Curve
A fragility curve (FC) is defined as the exceedance probability of a threshold that identifies the attainment of a Limit State for a given intensity measure (im): where the demand is obtained from Equation (35) and the capacity pertaining to different Limit States is in Equation (18). This is one of the merits of this study: instead of using a purposely defined metric, such as, for example, a limit interstorey drift, the global displacement capacity is consistently obtained from mechanics-based considerations, as explained in the previous sections, and is therefore related to Limit States. To compare the so-derived fragility curves with those developed in the literature, which are based on qualitative descriptions of damage levels, some correspondence should be found with the Limit States here adopted.
In empirical procedures, descriptive damage scales are preferred in reconnaissance activities to produce post-earthquake damage statistics. Some of the most frequently used damage scales are: HCR [5], HAZUS99 [18], Vision2000 [19], EMS98 [20], and ATC-13 [21]. Table 2, adapted from Ahmad et al. [22], shows a comparison between these damage scales and the Limit States here adopted. Analytical fragility assessment methodologies are commonly based on two main components: (a) the intensity-measure-to-structural-response function, and (b) the structuralresponse-to-damage-state functions, which result from structural and damage analysis procedures, respectively. Each procedure is affected by uncertainty. When developing FCs from non-linear static analyses, as done in this work, the uncertainty relevant to (b) is explicitly accounted for in all structural parameters: material properties, element sizes, and reinforcement details. The uncertainty related to (a), which is mainly due to record-torecord variability, can be added a posteriori using, for example, the statistically assessed values suggested in FEMA P695 [23]. On the other hand, when developing FCs from non-linear dynamic analyses, the record-to-record uncertainty is explicitly accounted for, whereas including the uncertainty relevant to (b) requires significant computational effort and is usually avoided, under the common assumption that uncertainties of (a) largely exceed uncertainties of (b).

Selection of Frame Typologies
The frame typologies considered in this study are classified based on the number of storeys, ranging from 1 to 5, and then sub-classified based on the number of bays and bay lengths as given in Table 3 and Figure 4. For each frame typology, two sub-typologies have been considered, representing two different construction periods, 1991-2000 denoted as "new", and 1961-1970 denoted as "old". Each sub-typology is characterized by a different range of material properties, concrete strength f cm , steel strength f ym , and of shear reinforcement, stirrup spacing s st , stirrup diameter φ and flexural reinforcement A s (Table 4). Table 3. Global geometry parameters of the frame typologies in Figure 4.    Table 3. Table 3. Global geometry parameters of the frame typologies in Figure 3. Type 1  Type 2  Type 3  Type 4  Type Table 3.

Effects of Soil Class and Location
Fragility Curves (FC) of a given building typology are affected by the soil class. For the purpose of this study, three soil classes have been considered, following the definition in the Italian Building Code (NTC 2018) [24], depending on the equivalent shear wave velocity V S,eq (in m/s): A (V S,eq > 800), B (360 < V S,eq ≤ 800), and C (180 < V S,eq ≤ 360). To each soil category corresponds a set of coefficients, namely, F o and T * c , that define the elastic spectral shape. The values of these coefficients are provided in the Italian Building Code (NTC 2018) [24] for nine different return periods, with linear interpolation to obtain F o and T * c for intermediate periods. This implies that the spectral shape change with the soil class. Therefore, fragility curves developed on a certain soil class need to be adjusted if used in a different one. Such a spectral-shape-changing effect is accounted for in the approach here presented.
As an example, to show such an effect, the FCs pertaining to a given building and location are shown in Figure 5 as obtained on three different soil classes.
in the Italian Building Code (NTC 2018) [24], depending on the equivalent shear wave velocity , (in m/s): A ( , > 800), B (360 < , ≤ 800), and C (180 < , ≤ 360). To each soil category corresponds a set of coefficients, namely, and * , that define the elastic spectral shape. The values of these coefficients are provided in the Italian Building Code (NTC 2018) [24] for nine different return periods, with linear interpolation to obtain and * for intermediate periods. This implies that the spectral shape change with the soil class. Therefore, fragility curves developed on a certain soil class need to be adjusted if used in a different one. Such a spectral-shape-changing effect is accounted for in the approach here presented.
As an example, to show such an effect, the FCs pertaining to a given building and location are shown in Figure 4 as obtained on three different soil classes.
It is worth noticing that the soil class not only affects the median of the FCs, but also the dispersion, so that, for the same value of the PGA, largely different exceedance probabilities are found. Moreover, it should be noticed that the relationships between and * and the return period change from location to location, based on the local hazard. This implies that the spectral shape pertaining to a given soil class for a certain return period varies throughout the territory. Therefore, fragility curves developed at a certain location must be adjusted if used at a different one, even though they refer to the same soil class. It is worth noticing that the soil class not only affects the median of the FCs, but also the dispersion, so that, for the same value of the PGA, largely different exceedance probabilities are found.
Moreover, it should be noticed that the relationships between F o and T * c and the return period change from location to location, based on the local hazard. This implies that the spectral shape pertaining to a given soil class for a certain return period varies throughout the territory. Therefore, fragility curves developed at a certain location must be adjusted if used at a different one, even though they refer to the same soil class.
As an example, to show such effect, three different locations in Italy have been chosen: Ferrara, Reggio Calabria, and L'Aquila. For each location, FCs are developed for soil A, B and C. The results in Figure 6 show that, for each soil class considered, the FCs pertaining to a given building change significantly among the three chosen locations.
Buildings 2021, 11, x FOR PEER REVIEW 13 of 21 As an example, to show such effect, three different locations in Italy have been chosen: Ferrara, Reggio Calabria, and L'Aquila. For each location, FCs are developed for soil A, B and C. The results in Figure 5 show that, for each soil class considered, the FCs pertaining to a given building change significantly among the three chosen locations.

Monte Carlo Analyses
The fragility curves are built for each selected typology (see Section 4.2) and for each selected soil class (see Section 4.3) through Monte Carlo analysis [25,26]. It is worth mentioning that the location effect has not been accounted for and that a specific location has

Monte Carlo Analyses
The fragility curves are built for each selected typology (see Section 4.2) and for each selected soil class (see Section 4.3) through Monte Carlo analysis [25,26]. It is worth mentioning that the location effect has not been accounted for and that a specific location has been selected, i.e., Rome, for the purpose of focusing on the building typologies and the soil effects. The cumulative probability of exceedance of a Limit State xLS (with x = D, L, C) is expressed as function of the horizontal peak ground acceleration, a g , for reasons of compatibility with the commonly used hazard maps in Italy. For each LS considered, for a given intensity measure a g,k , N it Monte Carlo analyses are performed and a set of displacement-based capacity/demand ratios ρ = d xLS /d i a g,k is obtained and subsequently statistically elaborated as: where P xLS is the exceedance probability of the xLS (with x = D, L, C) for the kth seismic intensity a g,k , with N it = 1000 and k = 1, . . . , 52, corresponding to a g = 0.01g, . . . , 1.02g. Once the probability for each a g,k is obtained, LSE non-linear regression [27][28][29] is carried out to determine the parameters µ and V T of the lognormal cumulative distribution function expressed in terms of a g : where Φ(·) is the standard normal cumulative distribution function, µ is the median of the capacity, and V T is the total coefficient of variation, including both capacity and demand variability.

Resulting Fragility Curves
The resulting soil-dependent FCs for soil classes A, B and C are shown in Figures 7-9, respectively, where it should be noticed that, for each frame typology, a fragility fuse is obtained. The fuse is the result of the variation of the global geometry (Table 3), while the dispersion of each FC is the result of the uncertainties in material properties, beam/column sectional sizes, shear and flexural reinforcement. The effect of the soil class on the FCs at various Limit States (LS) is shown in Figure 10.
Building height/number of stories is another parameter that influences the shape of the fragility curves, whereby the fragility of the frames increases with the increasing number of stories. Table 5 reports the FC parameters (mean µ G and coefficient of variation V G of the PGA-based capacity) for three Limit States, given for different soil classes (A, B, C), different construction types ("Old" and "New") and with number of storeys ranging from 1 to 5. The mean is given as a range around its mean value ± the distance from the fuse boundaries.
As a further step, the FCs are aggregated across the soil classes, so as to obtain wider fuses as shown in Figure 10. On passing, these fuses are akin to the observational FCs, which are usually constructed without disaggregating the soil class. Table 6 summarizes these FC parameters (mean µ SG and coefficient of variation V SG of the PGA-based capacity) for three Limit States, including soil type variation, for different construction types ("Old" and "New") and with number of storeys ranging from 1 to 5. The mean is given as a range around its mean value ± the distance from the fuse extremes and V SG is determined as the average of the V G of soil A, B and C. mentioned in Section 4.1, the demand variability should be added as follows: where is the coefficient of variation due to record-to-record variability, which, according to FEMA P695 [23], can be estimated as 0.2 for DLS and 0.4 for LLS and CLS. Table 7 reports the FC parameters and .
(a) (b)      Table 5. Fragility curve parameters (mean , fuse width around the mean, and coefficient of variation of the PGAbased capacity) for three Damage/Limit States, for three soil types (A, B, C), two construction types ("Old", Gravity Load Design, and "New", code-based design) and with number of storeys ranging from 1 to 5.   Table 5. Fragility curve parameters (mean µ G , fuse width around the mean, and coefficient of variation V G of the PGA-based capacity) for three Damage/Limit States, for three soil types (A, B, C), two construction types ("Old", Gravity Load Design, and "New", code-based design) and with number of storeys ranging from 1 to 5. To this point, the capacity-related variability V SG stems directly from the analysis. As mentioned in Section 4.1, the demand variability V D should be added as follows:

Soil Type Design Type
where V D is the coefficient of variation due to record-to-record variability, which, according to FEMA P695 [23], can be estimated as 0.2 for DLS and 0.4 for LLS and CLS. Table 7 reports the FC parameters µ SG and V T . Table 6. Fragility curve parameters including soil type (mean µ SG , fuse width around the mean, and coefficient of variation V SG of the PGA-based capacity), for three Damage/Limit States, two construction types ("Old", Gravity Load Design, and "New", code-based design) and with number of storeys ranging from 1 to 5.  Table 7. Fragility curve parameters including soil type and demand variability (mean µ SG , fuse width around the mean, and total coefficient of variation V T ), for three Damage/Limit States, two construction types ("Old", Gravity Load Design, and "New", code-based design) and with number of storeys ranging from 1 to 5.

Comparison with Literature Fragility Curves
The validity of the developed analytical FCs is appraised by comparing them with some observational and analytical FCs available in the literature (Figure 11).
Del Gaudio et al. [30] derived observational-based FCs from the damage distribution data of 7857 buildings in L'Aquila region collected by the Italian Civil Protection Department after the 2009 earthquake. In particular, FCs for low-and medium-rise structures of different types are described. The damage states considered are DS1, DS2, and DS3, which in this work are made to correspond to DLS, LLS and CLS.
De Luca et al. [31] proposed FCs based on the damage data of 131 buildings in the neighborhood of Pettino after L'Aquila earthquake in 2009. sions in the fragility fuses than those obtained in this study by disaggregating the soil classes (see Figure 6 -9). A more consistent comparison is therefore done with respect to the fragility fuses obtained in this study by aggregating the soil classes (see Figure 9).
Such comparison is done in Figure 11, where the blue, orange and red solid lines represent the soil-aggregated fuses for the three LSs. With reference to Figure 9, the lower boundary is determined as the rightmost FC pertaining to soil class A and the upper boundary is determined as the leftmost FC pertaining to soil class C. In Figure 11, having in mind the correspondence between DS1, DS2 and DS3, to DLS, LLS and CLS, respectively, shown in Table 2, it can be observed that: (a) for DLS, two out of four literature FCs fully fall within the fragility fuse, with the notable exception of that by Del Gaudio et al. [30] which significantly underestimates the DLS-exceedance Rossetto and Elnashai [5] developed FCs for European RC structures based on the damage data collected after 19 seismic events, following the HRC damage scale.
Ahmad et al. [22] developed analytical FCs using displacement-based pushover methodology for earthquake loss assessment.
All observational FCs are developed based on statistical analysis of the damage data collected from a set of buildings after seismic events, without distinguishing among different soil classes. This heterogenous collection of data leads to significantly wider dispersions in the fragility fuses than those obtained in this study by disaggregating the soil classes (see . A more consistent comparison is therefore done with respect to the fragility fuses obtained in this study by aggregating the soil classes (see Figure 10).
Such comparison is done in Figure 12, where the blue, orange and red solid lines represent the soil-aggregated fuses for the three LSs. With reference to Figure 10, the lower boundary is determined as the rightmost FC pertaining to soil class A and the upper boundary is determined as the leftmost FC pertaining to soil class C. All observational FCs are developed based on statistical analysis of the damage data collected from a set of buildings after seismic events, without distinguishing among different soil classes. This heterogenous collection of data leads to significantly wider dispersions in the fragility fuses than those obtained in this study by disaggregating the soil classes (see Figure 6-9). A more consistent comparison is therefore done with respect to the fragility fuses obtained in this study by aggregating the soil classes (see Figure 9).
Such comparison is done in Figure 11, where the blue, orange and red solid lines represent the soil-aggregated fuses for the three LSs. With reference to Figure 9, the lower boundary is determined as the rightmost FC pertaining to soil class A and the upper boundary is determined as the leftmost FC pertaining to soil class C. In Figure 11, having in mind the correspondence between DS1, DS2 and DS3, to DLS, LLS and CLS, respectively, shown in Table 2, it can be observed that: (a) for DLS, two out of four literature FCs fully fall within the fragility fuse, with the notable exception of that by Del Gaudio et al. [30] which significantly underestimates the DLS-exceedance In Figure 12, having in mind the correspondence between DS1, DS2 and DS3, to DLS, LLS and CLS, respectively, shown in Table 2, it can be observed that: (a) for DLS, two out of four literature FCs fully fall within the fragility fuse, with the notable exception of that by Del Gaudio et al. [30] which significantly underestimates the DLS-exceedance probability under high PGA values (as a matter of fact, their "light" FC crosses the "extensive" FC, as shown in Figure 11, which is clearly inconsistent), (b) for LLS, basically three out of four literature FCs fall within the fragility fuse, again with the notable exception of that by Del Gaudio et al. [30], which also in this case significantly underestimates the LLS-exceedance probability under high PGA values, (c) for CLS, two out of four literature FCs fall within the fragility fuse, where Ahmad et al. [22] significantly overestimate and Del Gaudio et al. [30] significantly underestimate the CLS-exceedance probability.

Conclusions
Fragility Curves are a useful tool for conducting risk analyses at an urban, provincial, or regional scale of territorially diffused building stocks. One of the methods to develop fragility curves relies on statistical processing of the data numerically produced by means of models representing the construction typologies in the region of interest. The so-obtained analytical fragility curves have several advantages over other fragility curves obtained from the statistical treatment of observational damage data, such as, for example: (a) ability to quantitatively define the attainment of certain damage state, as opposed to the qualitative description of observational methods, (b) objectivity of the database, as opposed to the surveyor-dependent observational database, (c) completeness of the database, which is built spanning throughout the range of all basic variables, as opposed to the incomplete database of observational data, (d) flexibility of the database, since analyses can be conducted on any region, even for those where no observational data are available. So far, in the literature, these fragility curves have been obtained: (i) by non-linear dynamic analysis conducted on detailed numerical finite element models, or (ii) with simplified models.
In this study, the second approach was undertaken, by developing a simplified mechanical model that ensures a good compromise between accuracy, representativeness, and computational effort. The simplified model assumes shear-type behavior, which is nonetheless acceptable for most of the existing reinforced concrete buildings. The global displacement-based capacity is analytically derived directly from the local capacity of structural elements and resisting mechanisms, through a simplified closed-form push-over analysis that identifies the attainment of the three global Limit States of interest: Damage, Life Safety, and Collapse. The advantage of the developed method over other simplified pushover analyses in the literature [11,22] is both in its non-iterative approach and in its closed-form nature, which renders the method computationally more efficient, making it suitable for conducting risk studies on a large scale. This allows implementation of the developed assessment methodology, either by hand calculations or in a spreadsheet. The potential of this simplified mechanical model has been shown by applying it on selected typologies of reinforced concrete frames, defined based on number of storeys and construction periods. The latter have been chosen to represent different seismic design compliance levels, from none to full. Within each typology, frames have been randomly generated by varying global and local geometry parameters, reinforcement quantities and material properties, within ranges that are representative of the Italian construction practice in different periods.
Moreover, in this study, a fundamental issue that has been dealt with is the following: fragility curves pertaining to a given typology/building differ depending on the location and the soil class. This is because different locations are characterized by different hazard curves and because different soil classes are characterized by different spectral shapes. These combined effects produce significant changes in the fragility curves, which cannot be accounted for in a simplistic way.
Bearing this in mind, analytical fragility curves for different building typologies have been derived for different soil classes and the parameters of the fitted lognormal curves are summarized in Table 5, where only the within-typology variability is considered, Table 6, where the soil influence is also included, and Table 7, where the demand variability is also included. Such curves have been compared with some observational and analytical curves from the literature, showing acceptable agreement.
The main conclusions drawn from this study can be summarized as follows: Location and soil class influence: when developing analytical fragility curves, the influence of the local hazard curve and of the local soil class must be considered. FCs pertaining to the same typology/building change when used at different locations and/or on different soil classes. This induces significant errors on risk and scenario studies at the territorial level. To carry out this study in a more effective manner, a strategy is under development aiming at transforming, through analytical closed-form functions, FCs developed on a certain location and soil class to another one. This is beyond the scope of the present study and will be discussed and presented in a future article. Construction age: RC frames fragility is significantly dependent on the construction age when it spans from pre-seismic-code to seismic-code periods. The fragility curves of the two epochs provide insightful information about the vulnerability features of the structures with respect to their construction age. Building height: the building height/number of stories is a crucial parameter to the evaluation of the fragility curves, since the LS-exceedance probability increases with the height/number of stories (here, studied only up to five). Thus, it is effectively used as a key parameter to define different building typologies.