Regional seismic risk assessment is paramount in earthquake-prone areas, for instance, to define and implement prioritisation schemes for earthquake risk reduction. As part of the Indonesia School Programme to Increase Resilience (INSPIRE), this paper proposes an ad hoc rapid-visual-survey form, allowing one to (1) calculate the newly proposed INSPIRE seismic risk prioritisation index, which is an empirical proxy for the relative seismic risk of reinforced concrete (RC) buildings within a given building portfolio; (2) calculate the Papathoma Tsunami Vulnerability Assessment (PTVA) index, in any of its variations; (3) define one or more archetype buildings representative of the analysed portfolio; (4) derive detailed numerical models of the archetype buildings, provided that the simulated design is used to cross-check the model assumptions. The proposed INSPIRE index combines a baseline score, calibrated based on fragility curves, and a performance modifier, calibrated through the analytic hierarchy process (AHP) to minimise subjectivity. An attempt to define a multi-hazard prioritisation scheme is proposed, combining the INSPIRE and PTVA indices. Such a multi-level framework is implemented for 85 RC school buildings in Banda Aceh, Indonesia, the most affected city by the 2004 Indian Ocean earthquake–tsunami sequence. As part of the proposed framework, two archetype buildings representative of the entire portfolio are defined based on the collected data. Their seismic performance is analysed by means of non-linear static analyses, using both the analytical simple lateral mechanism analysis (SLaMA) method and numerical finite-element pushover analyses to investigate the expected plastic mechanisms and derive displacement/drift thresholds to define appropriate damage states. Finally, non-linear dynamic analyses are performed to derive fragility curves for the archetype buildings. This paper demonstrates the effectiveness of the INSPIRE data collection form and proposed index in providing a rational method to derive seismic risk prioritisation schemes and in allowing the definition of archetype buildings for more detailed evaluations/analyses.
Regional seismic risk assessment is paramount in highly earthquake-prone areas. In fact, in several countries around the world, a large portion of the building stock has been designed and/or constructed according to obsolete structural codes, which include little-to-no provisions for earthquake resistance and detailing. Several reinforced concrete (RC) buildings fall in this category, and they often represent the highest share for both residential and commercial occupancy in many countries (e.g. in Italy 48 % of the buildings constructed after 1971 are made of RC; ISTAT 2011). RC structural systems are also widely used in the design of critical infrastructure, such as hospitals and school facilities. Those are the focus of this paper. Clearly, it is desirable that any risk-mitigation strategy designed by governmental agencies should be based on a rational understanding of the risk of large building groups – or portfolios – at a country level (or in a smaller region). However, it is cost-ineffective to perform detailed structural simulations for a large amount of structures, given the shortage of both financial and technical and/or computational resources. Therefore, a multi-level approach is usually preferred, starting with a screening based on simplified and rapid methods and performing more detailed structural analyses only for selected groups of structures at higher relative risk and for which an archetype (or index) building can be identified (e.g. FEMA P-154, 2015; Benedetti and Petrini, 1984; Grant et al., 2007).
Common approaches for regional seismic risk assessment of RC buildings (see Sect. 2 for more details) refer to typological approaches based on predetermined building categories (e.g. Giovinazzi and Lagomarsino, 2004) or the use of rapid-visual-survey (RVS) forms and calibrated empirical seismic vulnerability/risk indices (e.g. Uva et al., 2016). Although these approaches rely on various assumptions, such simplified methods provide valuable proxies to develop prioritisation schemes (i.e. performing a ranking of the buildings in a given portfolio based on their relative vulnerability or risk-related metrics). As discussed, such simplified methods include some degree of subjectivity by the analyst, mainly reflected in the choice and the assigned relative importance of the parameters involved in the analysis. Moreover, given the low amount of information required, such methods do not allow one to further refine the analysis, providing a more detailed, second-level seismic risk assessment. Finally, those methods mostly refer to seismic hazard, which in some countries might not be enough for the development of a rational multi-hazard prioritisation scheme.
Numerous evidences of previous natural-hazard events have highlighted the vulnerability of school infrastructure to natural hazards and particularly to earthquake-induced ground shaking. From the structural and architectural points of view, school buildings are especially vulnerable given structural characteristics that typically include large rooms, large windows (particularly in tropical climates), and corridors, all of which may represent seismic vulnerability factors. At the same time, schools play a critical role in the education of a community's next generation, with school children being one of the most vulnerable components of the society due to their age and their developmental stage. A safer and resilient school can save valuable lives of children and help to bring normality back to society in times of disaster. These considerations set school buildings apart from their peers in terms of priority for assessment and resource allocation for structural retrofitting. In fact, some of the world disaster reduction campaigns by the United Nations International Strategy for Disaster Reduction (UNISDR) were carried out together with various partner organisations under the theme of “Disaster Risk Reduction Begins at School” (UNCRD, 2009). Recently, the Comprehensive School Safety Framework (CSSF; GADRRRES, 2017) has proposed an integrated approach to reduce disaster risk and promote resilience in the education sector. The CSSF is funded on three pillars: “Safe Learning Facilities” (including “implementing assessment and prioritisation plans for retrofitting or replacing unsafe schools, including relocation”), “School Disaster Management” and “Risk Reduction and Resilience Education”.
Based on the above discussion, a new RVS form and a seismic risk prioritisation index for RC buildings are proposed in this study to address the above-mentioned gaps. Both the RVS form and the seismic index are the first outcomes of the INSPIRE project (Indonesia School Programme to Increase Resilience). INSPIRE looks to develop an advanced, harmonised and science-based risk assessment framework for school infrastructure in Indonesia subjected to cascading earthquake–tsunami hazards. It also assesses the effectiveness of different soft (e.g. risk-reduction education) and hard (e.g. structural retrofitting) mitigation measures in reducing casualties, economic loss and disruption to school infrastructure, thus increasing community resilience.
INSPIRE rapid-visual-survey form completed for a building with ID 14A (Sect. 4.1).
The INSPIRE RVS form (Fig. 1) is designed to be completed by trained engineers in approximately 20–30 min – depending on the size of the building – by means of a sidewalk survey. This is a one-page form including various sections related to the general identification and geolocation of the building, its geometric properties (including space for sketching the building's shape and footprint), and its structural characteristics and deficiencies, including the structural typology and the details of the main structural members. It is also possible to assign a confidence level for each parameter, allowing for a better classification and weighting of the data after a campaign of RVSs. The back of the form is used to provide definition of both the parameters and the confidence levels and provides a blank space that can be used to register extra information. The collected data are fully compatible with both the Global Earthquake Model (GEM) Building Taxonomy (Brzev et al., 2013) and the Hazard United States (HAZUS) model (Kircher et al., 2006). The collected information allows one to (1) calculate the proposed INSPIRE seismic risk prioritisation index, introduced in Sect. 3; (2) calculate the Papathoma Tsunami Vulnerability Assessment index (PTVA4, Dall'Osso et al., 2016, in any of its variations, described in Sect. 4.2); (3) define one or more representative archetype buildings consistent with the local building codes and practice; and (4) derive detailed numerical models of the archetype buildings, provided that the simulated design is used to cross-check the model assumptions. The unique features of the INSPIRE form include (1) the possibility of considering both seismic and tsunami hazards while requiring a reasonable amount of time to complete the survey; (2) the consideration of confidence levels for the parameters, which is particularly useful in deriving statistics; and (3) the possibility to be expanded to include other hazards with simple modifications/customisations.
In particular, the INSPIRE seismic risk prioritisation index (Sect. 3) aims at providing a simple method to derive a prioritisation scheme, minimising the subjectivity involved in the calculation. In fact, mechanics-based fragility functions are used to define a baseline score. A performance modifier is defined based on parameters that can jeopardise the seismic performance of a building (e.g. presence of short columns, pounding potential). The weight assigned to each parameter is defined through the analytic hierarchy process (AHP; Saaty, 1980), providing a mathematically consistent and rational solution to the weighting process.
In this study, the INSPIRE RVS form and proposed multi-hazard risk prioritisation index are applied to a portfolio of 85 RC school buildings in Banda Aceh, Indonesia, highlighting the simplicity and rapidity of the whole process. Moreover, consistently with the proposed multi-level framework, a detailed analytical and numerical seismic fragility assessment is provided for the identified archetype school building, demonstrating the effectiveness of the INSPIRE RVS form in providing inputs and allowing more detailed analyses.
Various methodologies to derive prioritisation schemes for buildings based on their relative seismic vulnerability/risk are available in the scientific literature and or international standards/guidelines. Each of these is characterised by a different underlying approach, basic assumptions and/or applicability conditions. Most of these are based on the calculation of a seismic vulnerability/risk index for a portfolio of buildings. The results can be used to rank those buildings, defining a priority list of structures requiring further investigation. A comprehensive overview of this previous research is not within the scope of this paper; however, a briefing of some relevant past initiatives is presented here.
The procedure proposed in the guidelines by the Applied Technology Council (ATC 3-06, 1978) uses a strength-based approach to define an earthquake capacity ratio, comparing the actual strength of the building to the code requirement for new ones. Adjustments are also adopted to consider in situ material properties and insufficient detailing (compared to modern design). Such a capacity-to-demand ratio is defined as the earthquake capacity ratio, and it is calculated as the minimum of the component-by-component strength ratios.
The Alaska's Department of Education (1997), among others, has produced surveying forms to assess the structural conditions of buildings and the associated seismic vulnerability, with focus on school buildings. Such forms mainly consist of checklists investigating areas of potential concerns for seismic vulnerability. On the other hand, the procedure introduced by the Federal Emergency Management Agency (FEMA P-154, 2015) is based on a rapid visual screening of buildings and a two-level approach for a fast assignment of a seismic vulnerability index (which requires no mechanical calculation from the user). The companion FEMA P-155 (2015) describes the rationale behind the scoring system, which is directly connected to the probability of collapse of archetype building categories. Such a method is based on the HAZUS framework (and typological force-displacement curves) to define the building categories and to derive a seismic-only assessment.
The New Zealand Society for Earthquake Engineering (NZSEE) defines an initial evaluation procedure (IEP) in the NZSEE guidelines of 2017, providing a broad indication of the seismic rating of a building based on a sidewalk survey. The score is expressed in terms of the percentage of the new building standard (%NBS), which is the ratio of the displacement capacity of the building for the life safety limit state over the minimum capacity required for a new building for the same limit state. Firstly, a baseline %NBS is calculated using specifically tabulated coefficients relating to year of design, strengthening interventions, importance of the structure, assumed ductility capacity, site hazard, presence of near-fault effects, soil type, etc. It is assumed that the capacity of the building cannot be lower than the minimum specified by the code valid for the year of design. This is only true if building codes are legally enforced. If this is not the case, the seismic vulnerability assessment may be characterised by a higher level of epistemic uncertainty. The baseline value is reduced if structural weaknesses are present (e.g. pounding potential, soil characteristics, presence of vertical and plan irregularity, presence of short columns).
The Italian National Group for Earthquake Defence (Gruppo Nazionale per la Difesa dai Terremoti, GNDT) has also provided
a seismic vulnerability index (Benedetti and Petrini, 1984; Angeletti et
al., 1988) based on simple assessment forms including, among other
parameters, the structural material, the typology of the lateral-load-resisting system (LLRS), the quality of the building materials and the
overall construction, and the existing damage level (if any). The
vulnerability index is based on a weighted sum of such parameters and is
defined in the range [0, 100] for masonry and [
Other rapid surveying forms and rapid procedures have been proposed by different authorities and organisations, such as the World Health Organization (WHO) and the United Nations (UN), with special focus on developing countries. For instance, Dhungel et al. (2012) collected and assessed the physical condition of 1381 school building units in Nepal. The data were collected by mobilising the school teachers; school vulnerability, calculated on the basis on the empirical weightage on different factors (e.g. structural material, number of storeys, shape of the roof), was used to estimate the possible damage, casualties and injuries for earthquakes of different seismic intensities.
Finally, a broader perspective is provided in the work by Grant et al. (2007), which proposes a four-level prioritisation framework focused on school buildings, filtering the buildings with increased level of detail, allowing policy-makers to choose the filtering thresholds based on the available resources. The idea is to firstly check the code-based demand deficit of the buildings, comparing the new structural code demand to the code appropriate for the year of construction, in terms of PGA. For the buildings with a deficit above a given threshold, the GNDT index is calculated. The buildings with the highest rating are assessed with a simplified mechanics-based procedure, and finally the most critical ones are assessed with structural models providing full details.
The INSPIRE index (
The HAZUS MH4 framework provides, among other information, an extensive set
of fragility curves representing the seismic performance of archetype
buildings which are classified based on four parameters: material (Mat), BSS,
building height and seismic code level (defined according to the Uniform Building Code, UBC, 1994; ICBO, 1994). Adopting the HAZUS fragility database as a reference is
further motivated by the fact that several countries around the world have
adopted seismic provisions which are, to different extents, consistent with
the recommendations of the UBC 1994. For each archetype building category,
four fragility functions are provided in HAZUS, respectively for the slight,
moderate, extensive and complete damage states (DS), or limit states. The
fragility functions are defined as log-normal cumulative distribution
functions, or CDFs (Eq. 2), in terms of different
intensity measures, among which is the peak ground acceleration (PGA). The
curves are defined by a median PGA (
For the scope of this paper, the HAZUS fragility database has been filtered to consider only the curves related to RC buildings (namely, categories C1, C2, C3). Moreover, only the fragility curves corresponding to the extensive damage limit state (DS3) have been considered, which are mainly related to the life safety performance objective according to modern seismic codes (e.g. ASCE 7–10 – ASCE, 2010; Eurocode 8, or EC8 – CEN, 2004; NZSEE, 2017). The selected fragility curves are reported in Table 1. Details of the involved parameters are provided in Table 2.
Selected fragility curves from the HAZUS MH4 framework (Kircher et al., 2006).
n.a.: not assigned.
Description of the HAZUS MH4 categories involved in the INSPIRE index (modified after Kircher et al., 2006).
To define the baseline score of the INSPIRE index, for each considered
archetype building category, and its corresponding fragility curve, the
probability to exceed DS3 (extensive damage) is calculated
(Fig. 2a) for three levels of PGA, 0.1, 0.25 and
0.4 g, respectively corresponding to low, moderate and high seismicity
levels. The analyst will select the seismicity level appropriate for the
considered building portfolio/geographic area. It is worth noting that, in
modern seismic codes, DS3 is related to the life safety performance
objective. The process above is repeated for each archetype building
category in the HAZUS database, such that it is possible to map the building
basic parameters to the exceeding probability of DS3, conditional to the
considered PGA value (
Eight secondary parameters are used to define the performance modifier. These parameters are not explicitly considered in the HAZUS framework (and
therefore in the baseline score) but, if present, can
jeopardise the seismic performance of a given building. Firstly, based on
Table 3, a score is assigned for each secondary
parameter (
Secondary parameters: definition, available alternatives, related scores and weights.
The secondary parameters defining the performance modifiers are deemed to complement the information in the HAZUS fragility curves. The idea is that the baseline score (HAZUS fragility database) provides the (conditional) seismic risk of a given building class, while the secondary parameters are related to building-specific vulnerability factors. The secondary parameters have been selected based on the fundamental rules of seismic design (e.g. Paulay and Priestley, 1992) and the commonly observed post-earthquake damage on RC structures (e.g. Elnashai and Di Sarno, 2008; Palermo et al., 2017; De Luca et al., 2018). For each of them, Table 3 provides guidance on the selection of the alternatives. The score for each alternative has been defined based on a uniform partitioning of the range [0 %, 100 %].
Evidently, the secondary parameters defining the performance modifier do not have the same influence on the overall vulnerability and risk. For example, the absence of infill walls at the ground storey can lead to a soft-storey mechanism, which in turn can results in local or global collapse. On the other hand, the addition of storeys to a building can somehow increase its vulnerability, but it is generally less likely this would lead to collapse. Therefore, the weight of the former parameter should be higher than the latter, to reflect such a different effect on the overall seismic risk.
To this extent, the AHP, originally proposed by Saaty (1980), is used to calibrate the weights in the proposed procedure. This allows one to have a rational and mathematically consistent assignment of the weights, which is based on pairwise comparisons between the secondary parameters and uses eigenvalues theory. Hence, the selected weights can reflect the relative importance of each parameter with respect to the others in the determination of the performance modifier. It is worth mentioning that such an approach has been successfully adopted in other earthquake engineering applications, such as the selection of the optimal seismic retrofitting strategy for case-study buildings (Caterino et al., 2008). Also, it is worth mentioning that, in its current version, the expert judgement defining the weights adopted in the procedure is provided by the authors. However, such weights can be further updated considering the opinion of a group of experts in the field of structural and earthquake engineering.
After the definition of the parameters involved in the analysis, the first
step of the process consists of expressing an expert opinion about every
possible pairwise comparison of those parameters. Given two parameters
Scale of relative importance of the secondary parameters (Saaty, 1980).
Once the pairwise comparisons have been performed, the vector
Judgement matrix adopted for the calibration of the weights.
Therefore, the consistency of the pairwise comparison is measured by
calculating the consistency index (CI) using Eq. (5),
where
In many common situations, it is very likely that seismic considerations alone are not enough to define a robust prioritisation scheme for the decision-making process by a governmental agency or the owner of a large building portfolio. To deal with such a scenario, this paper presents a simple methodology to include other hazards in the prioritisation scheme.
Simplified analytical indices are available for the estimation of the
vulnerability of buildings to natural hazards different than seismic hazard,
such as tsunami (e.g. Dall'Osso et al., 2016), flood (e.g. Stephenson and
D'Ayala, 2014; Pazzi et al., 2016; Nassipour et al., 2018) and wind (e.g.
Womble et al., 2016). Once the desired single-hazard vulnerability
indices (
Definition of a two-dimensional (
According to the definition of the seismic risk prioritisation index given
in Sect. 3, it appears evident that a degree of subjectivity is always involved in the calculation.
The baseline score is based on the fragility functions available in the
HAZUS MH4 guidelines. Although such a model is largely based on numerical
analyses of building models, the results are provided for a limited number
of structural categories. The user should carefully select the category that
best matches with the characteristics of the considered buildings, with
special reference to the code level parameter. The weights needed for the calculation of the performance modifier are based
on a subjective set of pairwise comparisons between the secondary
parameters, although this is derived in a mathematically consistent and
rational way that allows us to minimise the chance of having randomly assigned
weights.
The ratio between the maximum possible baseline score and the maximum
possible performance modifier is arbitrarily set to unity, i.e.
It is worth noting that the subjectivity is an intrinsic component of any
simplified vulnerability index, including the INSPIRE index. Therefore,
rather than assessing a single building, these methods should only be used
to derive a proxy for the relative vulnerability (or risk) of buildings in a
given portfolio and to define a prioritisation scheme for possible risk-mitigation actions. However, in this section some measures are proposed to
control and minimise the subjectivity involved with the definition of the
index. An example of such measures is given in Sect. 4, in relation to a real building portfolio application.
It is virtually impossible to perfectly match the considered buildings with the archetype buildings in HAZUS. However, a careful examination of the characteristics of the considered buildings should be carried out to better select the appropriate HAZUS categories. Characteristics such as the presence of strong beams vs. weak columns should be considered, which can lead to a pre-code classification or documented structural retrofit measures, which can lead to higher code level classification. A review of the (history of the) structural and/or seismic codes appropriate for the considered buildings can considerably reduce the level of subjectivity. These can be compared to the different provisions in UBC1994, defining equivalence relationships. Any information related to the year of construction (or design) of the considered buildings is fundamental in such a process. Overall, it is deemed that the assignment of the HAZUS categories to the considered buildings should reflect their expected differences in their seismic performance rather than perfectly match the properties of the archetype in the HAZUS category definition.
As an alternative to the HAZUS definition, the baseline score can be
redefined adopting, if available, a portfolio-specific set of fragility
curves. To this aim, the OpenQuake platform (open-source,
Map of the school buildings in Banda Aceh, including non-surveyed ones. Background map by © OpenStreetMap contributors 2019. Distributed under a Creative Commons BY-SA License.
An illustrative set of weights needed to calculate the performance modifier is given in this paper. However, a case-specific AHP, for instance involving groups of local experts, can be performed to derive new weights that match the characteristics of the considered building portfolio. Such a procedure cannot remove the subjectivity related to the performance modifier but provides the user with a tool to reduce it and have a close match between the analysed building portfolio and the adopted weights.
Finally, the subjectivity related to the ratio between the maximum possible baseline score and performance modifier can be tested through a sensitivity analysis. For a given building portfolio, the idea is to repeat the calculation of the seismic index with different values of the maximum possible performance modifier. This allows one to check the reliability of the priority list to this assumption. If slight modifications in the maximum possible performance modifier lead to high differences in the resulting priority list, engineering judgement should be adopted to justify the final choice.
The case-study portfolio selected for this study consists of 85 RC buildings
belonging to 44 school compounds located in Banda Aceh, the capital and
largest city in the province of Aceh, Indonesia
(Fig. 4). Banda Aceh is located in the island of
Sumatra and, according to the 2000 census, has a population of 219 070 people (Seta, 2000). The city was severely affected by the 26 December 2004
Indian Ocean earthquake (moment magnitude,
Sample photos of the building portfolio (taken on 16–19 October 2018).
A RVS campaign through the INSPIRE form (Fig. 1) was carried out to collect administrative, geometric and mechanical data related to the investigated buildings. The RVS campaign (Fig. 5) was conducted by three teams composed of one experienced engineer and two final-year (one undergraduate and one postgraduate) local engineering students working for 4 d. The surveys were conducted in all the suburbs of Banda Aceh (Fig. 6a) to obtain a clear overview of the construction practice in the city. For all the surveyed buildings, the BSS is a reinforced concrete frame with infills. The majority (81 %) of the buildings in the portfolio have a rectangular plan, with the remaining 19 % composed of L-, C- or T-shaped plans (Fig. 6b). A total of 68 % of the surveyed buildings are two storeys high, while one- and three-storey buildings represent 15 % and 16 % of the portfolio, respectively (Fig. 6c). The year of construction for each building was retrieved from the school signboard, registry or by interview of the school principal. Figure 6d shows that few buildings (18 %) survived the 2004 earthquake–tsunami sequence; hence, the majority of the portfolio (57 %) was constructed between 2005 and 2011, while the remaining 25 % was built from 2013 onwards.
Statistics for the 85 surveyed school buildings.
The accurate knowledge of the year of construction is essential for the structural characterisation of the portfolio. In fact, such information can be coupled with the history/evolution of the available structural and/or seismic codes to derive minimum-by-law dimensions of the structural members, reinforcement detailing, level of considered vertical and lateral forces in the design, etc. In this specific case, the appropriate structural code for the whole portfolio is the SKBI 1.3.53.1987 (Badan Standardasi Nasional, BSN, 1987). The first seismic code appropriate for the region is the SNI 1726:2002 (BSN, 2002), which is inspired by the American Uniform Building Code 1997 (or UBC, 1997; ICBO, 1997), which also facilitates the compatibility with the HAZUS framework. Stricter seismic provisions were adopted when the seismic code was updated in the SNI 1726:2012 (BSN, 2012), which is fully consistent with the American ASCE 7–10 (ASCE, 2010). Therefore, apart from a minority (12) of buildings constructed before 2002, approximately half of the buildings are constructed according to the SNI 1726:2002 (lower) standards while the other half refers to the updated SNI 1726:2012, with nominally better nominal seismic performance.
The analysis above, based only on the information of the INSPIRE forms,
allows one to identify an archetype building which represents the construction
practice for school building in Banda Aceh. The archetype building is a
two-storey, rectangular RC building. Figure 7 shows
the geometrical characteristics of the archetype building geometry, which
are based on the modal (most frequent) values of the empirical distributions
(histograms) considering only the rectangular, two-storey buildings in the
portfolio (69). The archetype building has 10 bays in the longitudinal direction (one
for the staircase, three for each classroom). In the transverse direction
there are two 3.5 bays and a 2 m corridor bay. The storey height is equal to
3.5 m. Columns of
Two-storey, rectangular buildings: geometry trends and adopted
values for the archetype building.
Considering the same geometry, two different sets of structural details are provided, to reflect the SNI 1726:2002 and SNI 1726:2012 seismic codes. Table 6 provides examples of structural reinforcement for typical members. It is worth mentioning that, for the post-2012 archetype structure, the cross-section height of beams and columns are respectively 5 cm and 10 cm larger than the corresponding members in the pre-2012 archetype structure. The reinforcement of the structural members is selected by cross-checking visual information (see Fig. 5) with the outcome of both gravity-based and seismic-based simulated design according to the SNI codes. To this aim, acting loads are calculated considering permanent dead loads (according to the suggested material properties in SKBI 1.3.53.1987) and live load equal to 5 kPa (1 kPa for the roof). Gravity axial load on columns is calculated based on a tributary area approach. This cross-checking exercise has shown that the real observed amount of longitudinal reinforcement is greater than the minimum by code. On the other hand, based on the limited visual information available for the transverse reinforcement, no joint stirrups are conservatively considered for both the pre-2012 and post-2012 classes, regardless of the requirement in both codes.
Typical structural details for the archetype building(s).
Note:
Based on the data collected with the forms, the INSPIRE seismic risk prioritisation index is calculated for the whole portfolio. Moreover, the tsunami PTVA4 index (Dall'Osso et al., 2016) is calculated, finally combining these results to derive a multi-hazard index. It is worth mentioning that the resulting index values are arbitrarily categorised in three groups, respectively green, yellow and red tags by defining two threshold values for the various indices. The definition of such thresholds is essentially a subjective (often political) choice that shapes the prioritisation scheme, based for instance on resources availability. For a governmental agency, those can be calibrated estimating the average structural retrofit (or relocation) cost per building and defining the amount of available public funding in two or more time windows (e.g. 1 and 5 years) to obtain specified risk-reduction objectives. As a proof of concept, in this paper the thresholds are selected to be equal to 33 % and 66 % for the calculated seismic, tsunami or multi-hazard indices.
Application to the seismic and tsunami indices to 85 RC school
buildings in Banda Aceh, Indonesia.
The PTVA index, similarly to the proposed INSPIRE index, allows one to derive
the relative tsunami vulnerability of a building. It is calculated as a
weighted sum of two parameters: the structural vulnerability and the
so-called water vulnerability. The first parameter depends on three
factors: the type of the lateral resisting system, the depth of the flood
water at the building location and the degree of external protection (e.g.
presence of seawalls). The water vulnerability depends on the ratio of
the inundated-to-total number of storeys. This parameter is calculated using the inundation data from the 2004 Indian Ocean
tsunami. In particular, Iemura et al. (2012) provide 85 field-measurements of
the tsunami height (from the ground) for the city of Banda Aceh. For the
purpose of this paper, a linear regression (
Sensitivity analysis for the baseline-to-performance modifier
ratio. Maximum performance modifier equal to
Figure 8a and b show the application of the INSPIRE seismic risk prioritisation index to the considered portfolio. To calculate the baseline scores, the HAZUS fragility curves related to the C1 category “concrete moment frame” are used. Since the infills of the investigated frames are made of a single layer of poor-quality clay bricks, their presence is neglected. According to the analysis of the year of construction and the history of the structural and seismic codes in Indonesia, the categories pre-code and low-code are adopted for the pre-2012 and post-2012 buildings, respectively. Given the particularly small differences in the characteristics of the considered buildings, the INSPIRE index is particularly similar for the whole portfolio [32 %, 64 %]. The slight differences in the value of the index are due to the performance modifiers, mainly governed by short columns and/or pounding for the majority of the schools. This is a further confirmation of the uniformity of the construction practice for school buildings in Banda Aceh, also observed in other countries (e.g. Nassirpour et al., 2018; Zhou et al., 2018). It is worth mentioning that, for this particular portfolio, soil condition is not influencing the performance modifier, since the unfavourable soil parameter is set to “yes” for all the buildings. On the other hand, the results for the PTVA tsunami index (Fig. 8c, d) show a larger variability [7 %, 75 %] and a clear correlation between the distance from the coast and the relative tsunami risk. This result, although preliminary, might suggest that the distance from the coast can be used as a very simple proxy for the water vulnerability parameter in the tsunami PTVA index.
Given the above-mentioned trends, the combination of the seismic and tsunami indices (Fig. 8e, f) clearly indicates that the tsunami hazards play a substantial role in determining the prioritisation scheme for the school buildings in the city of Banda Aceh. Indeed, the overall trend of the multi-hazard results is practically equal to the trend of the tsunami index results. For instance, the developed maps could be used to identify safer areas where strategic buildings (e.g. schools or hospitals) should be located, increasing the awareness of vulnerabilities that could be integrated in emergency planning for critical infrastructure disruption (e.g. Pescaroli and Alexander, 2016).
Numerical modelling strategy.
To control the role of the subjectivity in calculating the INSPIRE seismic index, a sensitivity study is conducted. The seismic index is applied to the entire portfolio four times (Fig. 9), by considering that the maximum value of the performance modifier is equal to 20 %, 30 %, 40 % and 50 %. The maximum value of the baseline score is respectively equal to 80 %, 70 %, 60 % and 50 %.
The results of the seismic index are ranked in descending order of risk. The results clearly show that, for this portfolio, the overall priority list is rather insensitive to the maximum baseline-to-performance modifier ratio. Indeed, a small number of swaps in the priority list is sought and, with the same definition of the thresholds for the tags, the number of buildings in each category has a negligible dependency on the ratio above.
As discussed above, the INSPIRE form provides enough data to build refined numerical models for one or more selected buildings, in a multi-level approach. While the previous sections illustrate the portfolio-level approach (relative-risk prioritisation), this section illustrates a possible structure-specific application of the proposed framework. To this aim, the archetype building(s) defined in Sect. 4.1 is (are) further analysed by means of non-linear static and non-linear dynamic analyses to derive structure-specific fragility curves. Those are derived considering a two-dimensional representation of the longitudinal and transverse frames that compose the archetype building. As discussed above, due to their small thickness and poor quality, infills are not considered in the models. Both the pre-2012 and the post-2012 archetype buildings are analysed, leading to four different computational models.
Numerical modelling strategy (from Gentile et al., 2019a, Taylor & Francis, all rights reserved).
Non-linear static analyses results for the archetype building.
The simple lateral mechanism analysis (SLaMA) method is firstly used to obtain a first estimation of the non-linear force-displacement curve of the building. SLaMA (NZSEE, 2017; Gentile et al., 2019a, b, c, d) is a tool to derive both the expected plastic mechanism and the capacity curve of RC frame, wall and dual-system buildings by using a “by-hand” procedure (i.e. using an electronic spreadsheet). This allows one to identify potential structural weaknesses in the lateral resisting mechanism and to test the reliability of numerical computer models in capturing the most probable behaviour of a structure. It is worth mentioning that each beam and column in the system is characterised considering many possible failure mechanisms (i.e. flexure, bar buckling, lap-splice failure, shear), considering that the weakest link will govern the overall structural behaviour.
The results of SLaMA are compared to refined numerical pushover analyses carried out with the finite-element modelling software Ruaumoko (Carr, 2016). The adopted modelling strategy, previously validated on experimental results (Magenes and Pampanin, 2004), is based on a lumped plasticity approach and is summarised in Table 7 and Fig. 10. The characterisation of the structural members is consistent with the approach adopted for SLaMA. P-delta effects are not modelled, given that the building is just two storeys high and made of RC. Fully fixed boundary conditions are considered at the base, and floor diaphragms are modelled as rigid in plane. A uniform force profile is adopted.
Archetype building: fragility curves based on cloud analysis,
calculated for
Figure 11 shows the results of the non-linear static analyses. Firstly, the fundamental periods of such frames in the transverse and longitudinal directions are equal to 0.55 and 0.47 s for the pre-2012 archetype and to 0.5 and 0.44 s for the post-2012 one. The plastic mechanism (numerically based) of the pre-2012 archetype frames is very similar to the one related to the post-2012 ones, in both the transverse and longitudinal directions. In the transverse direction, the plastic mechanism, calculated at the onset of the extensive damage state (DS3), is characterised by the development of plastic hinges for the roof beams and base columns, as well as joint shear hinges at the first storey (mixed-sway mechanism). No brittle shear failure is registered for beams and columns (both considering SLaMA and the pushover analysis). The beam causing DS2 (yielding) is highlighted in a blue circle in Fig. 11a, while the joint panel causing DS3 is highlighted with a red circle. In the longitudinal direction, the DS3 plastic mechanism is characterised by a soft storey at the first storey, with one internal column causing DS2, but the DS3 displacement of the structure is limited by the ultimate drift in the beam–column joint, highlighted with a red circle in Fig. 11b and d. According to the similarities in the plastic mechanism, the capacity curves for the pre-2012 and post-2012 archetype frames are particularly consistent. It is worth mentioning that the DS3 base shear for the post-2012 archetype building is 16 % and 19 % greater than the pre-2012 archetype, respectively in the transverse and longitudinal directions.
On the other hand, negligible differences are recorded for the displacement at DS3, for both the pre- and post-2012 archetypes and both transverse and longitudinal directions. This is because, in all four cases, the failure of one joint panel causes the attainment of such a damage state. The displacement in the force-displacement curves is calculated at the effective height (Priestley et al., 2007), which is approximately equal to 5.50 m. It is evident that the structure has a particularly low displacement capacity, since this is limited by the low drift capacity of the joint panels. Finally, there is a satisfactory agreement between SLaMA and the numerically based pushover, especially considering that the ultimate displacement in the SLaMA curves represents DS3.
Fragility curves of the archetype frames based on cloud analysis
(
The same two-dimensional models are analysed by means of a cloud analysis (Jalayer and Cornell, 2009), which consists of a series of non-linear time history analyses considering a large database of unscaled ground-motion records. For this application, the SIMBAD database (Selected Input Motions for displacement-Based Assessment and Design, Smerzini et al., 2014) is adopted. SIMBAD includes 467 tri-axial accelerograms, generated by 130 worldwide seismic events (shallow crustal earthquakes with moment magnitudes ranging from 5 to 7.3 and epicentral distances ranging up to 35 km). A subset of 150 records is considered here to provide a statistically significant number of strong-motion records of engineering relevance for the applications presented in this paper. As in Rossetto et al. (2016), these records are selected by first ranking the 467 records in terms of their PGA values (by using the geometric mean of the two horizontal components) and then (arbitrarily) keeping the component with the largest PGA value (for the 150 stations with highest mean PGA). This record selection strategy is compatible with the adopted cloud analysis procedure and with the lack of specific, freely available, ground-motion databases for the considered case-study region. A rigorous site-specific, hazard-consistent record selection could be used for alternative non-linear demand estimation methods for probability-based seismic risk assessment (Jalayer and Cornell, 2009), such as multiple-stripe analysis. The models used for the numerical pushover analysis are also adopted for the cloud analysis. The hysteretic behaviour of beams and columns is characterised be the revised Takeda model (Saiidi and Sozen, 1979), with the columns having a thinner loop. On the other hand, the hysteretic behaviour of the beam–column joints is modelled using the modified Sina model (Saiidi and Sozen, 1979), which is able to capture their pinching behaviour.
The results of the dynamic analyses (150 non-linear time history simulations
for each of the four models) are used to plot a cloud of points in the
plane inter-storey drift (chosen as the engineering demand parameter, EDP) vs.
pseudo-spectral acceleration at the first fundamental period (
The results show that the seismic fragility of the post-2012 structures is reduced with respect to the pre-2012 ones (Table 8). On average, the fragility median of the post-2012 frames are 5.6 % and 4.3 % higher than the pre-2012 ones, respectively in the transverse and longitudinal directions. The related dispersion is reduced, on average, by 16.1 % and 14.3 % for transverse and longitudinal directions.
It is worth noting that the structure-specific fragilities obtained with the analyses presented herein could be used to redefine the INSPIRE index. However, this would require the derivation of numerical fragility curves for a much larger set of archetype buildings, consistently with the building typologies identified in Table 1. In the authors' opinion, adopting the refined fragilities as an input for the INSPIRE index definition would be inappropriate. Indeed, the proposed index is defined as a quick and practical tool for seismic vulnerability prioritisation of large building portfolios (level 1 analysis). A more time-consuming fragility analysis (level 2) should instead be used to derive quantitative seismic risk estimates for one or more selected building in the database and to design structure-specific risk-mitigation strategies (e.g. structural retrofitting).
This paper introduces the INSPIRE index, which is an empirical proxy for the relative seismic risk of buildings within a given building portfolio and allows us to define prioritisation schemes for risk-mitigation measures. The definition of such an index represents the first step of the wider framework of the Indonesia School Programme to Increase Resilience (INSPIRE), which aims to develop an advanced, harmonised and science-based risk assessment framework for school infrastructure in Indonesia subjected to cascading earthquake–tsunami hazards.
To this aim, a data collection form used for the rapid visual inspection of RC buildings is first developed and presented. Such a form allows one to calculate the INSPIRE seismic risk prioritisation index and the tsunami PTVA index (in any of its variations), to obtain a level of geometrical and mechanical information sufficient to define one or more archetype buildings (representative of the portfolio) and/or to build refined numerical models, provided that the simulated design is adopted to cross-check the available information.
The INSPIRE index is specifically calibrated for RC buildings, and consists of two parts: a baseline score and a performance modifier. The baseline score is based on the HAZUS MH4 fragility curves, while the performance modifier is based on the score of the building with regard to eight secondary parameters, which, if present, are deemed to jeopardise the building performance. To minimise subjectivity, the relative weight of the secondary parameters is defined according to the analytic hierarchy process. This allows us to have a rational and mathematically consistent assignment of the weights, which is based on pairwise comparisons between the secondary parameters and eigenvalues theory.
The INSPIRE form and seismic risk prioritisation index are adopted for the analysis of 85 RC school buildings in the city of Banda Aceh, Indonesia, which is located on the island of Sumatra, the area most affected by the 2004 Indian Ocean earthquake–tsunami sequence. The joint application of the INSPIRE seismic risk prioritisation index and the PTVA tsunami index allow one to define a clear and transparent rationale for any prioritisation schemes for school buildings in Banda Aceh. In fact, the relative seismic risk of the considered buildings is particularly similar, while the relative tsunami risk shows a strong dependence with the distance from the coast. Indeed, the results show that a multi-hazard-based priority list is mostly governed by the tsunami risk for the case-study portfolio.
The advantages of using the INSPIRE form are further demonstrated by defining two archetype buildings, representative of the portfolio, based on the RVS results. The seismic performance of the archetype buildings is firstly analysed by means of non-linear static analyses, both analytically using the simple lateral mechanism analysis (SLaMA) method and numerically using refined finite-element models. Finally, the archetypes are analysed by means of cloud analysis, performing non-linear dynamic analyses using 150 unscaled natural ground motions and deriving fragility curves.
The results in this paper demonstrate the effectiveness of both the INSPIRE RVS form and INSPIRE seismic risk prioritisation index in providing a rational method to derive a prioritisation scheme, which can be extended including multi-hazard considerations, and in allowing the definition of an archetype building for more detailed evaluations/analyses.
This study represents a first step toward a comprehensive framework for earthquake and tsunami vulnerability and risk assessment and the selection of optimal retrofitting strategies for school facilities in Indonesia, through a to a multi-criteria decision-making analysis. Future research will investigate the numerically based tsunami fragility of the archetype buildings adopting different approaches (e.g. Petrone et al., 2017) and a full seismic loss analysis considering non-structural components, which often represent the highest share of the seismic losses.
Data will be made available on request.
RG and CG conceived and designed the study, developed the proposed framework and wrote the manuscript; RG performed the analyses; RG, YI and IR collected the data with EM supervising the data collection; all the authors discussed the results.
The authors declare that they have no conflict of interest.
This study was performed in the framework of the INSPIRE (Indonesia School Programme to Increase Resilience) and i-RESIST (Increasing REsilience of Schools in Indonesia to earthquake Shaking and Tsunami) projects, funded by the British Council through the Newton Institutional Links scheme and Research England through the University College London (UCL) Global Challenges Research Fund (GCRF) Small Research Grants scheme.
The Tsunami and Disaster Mitigation Research Centre (TDMRC) at Syiah Kuala University is acknowledged for the technical support during the field work. Peter Sammonds, Karim Aljawhari and Nurmalahayati Nurdin from UCL are gratefully acknowledged for helping in planning and executing the data collection. Rifqi Irvansyah, Putri Istiqamah, Sausan Zahrah, Natasya Amalia, Fajarul Aulia and Ilham M. Siddiq are gratefully acknowledged for their help in the data collection.
This research has been supported by the British Council (Newton Institutional Links) and Research England (UCL Global Challenges Research Fund Small Research Grants).
This paper was edited by Maria Ana Baptista and reviewed by two anonymous referees.