Simulation-based PSHA for the Canterbury region (Cybershake v17.8)
Research papers, University of Canterbury Library
None
None
At 00:02 on 14th November 2016, a Mw 7.8 earthquake occurred in and offshore of the northeast of the South Island of New Zealand. Fault rupture, ground shaking, liquefaction, and co-seismic landslides caused severe damage to distributed infrastructure, and particularly transportation networks; large segments of the country’s main highway, State Highway 1 (SH1), and the Main North Line (MNL) railway line, were damaged between Picton and Christchurch. The damage caused direct local impacts, including isolation of communities, and wider regional impacts, including disruption of supply chains. Adaptive measures have ensured immediate continued regional transport of goods and people. Air and sea transport increased quickly, both for emergency response and to ensure routine transport of goods. Road diversions have also allowed critical connections to remain operable. This effective response to regional transport challenges allowed Civil Defence Emergency Management to quickly prioritise access to isolated settlements, all of which had road access 23 days after the earthquake. However, 100 days after the earthquake, critical segments of SH1 and the MNL remain closed and their ongoing repairs are a serious national strategic, as well as local, concern. This paper presents the impacts on South Island transport infrastructure, and subsequent management through the emergency response and early recovery phases, during the first 100 days following the initial earthquake, and highlights lessons for transportation system resilience.
This paper presents preliminary results of an experimental campaign on three beam-column joint subassemblies extracted from a 22-storey reinforced concrete frame building constructed in late 1980s at the Christchurch’s Central Business District (CBD) area, damaged and demolished after the 2010-2011 Canterbury earthquakes sequence (CES). The building was designed following capacity design principles. Column sway (i.e., soft storey) mechanisms were avoided, and the beams were provided with plastic hinge relocation details at both beam-ends, aiming at developing plastic hinges away from the column faces. The specimens were tested under quasi-static cyclic displacement controlled lateral loading. One of the specimens, showing no visible residual cracks was cyclically tested in its as-is condition. The other two specimens which showed residual cracks varying between hairline and 1.0mm in width, were subjected to cyclic loading to simulate cracking patterns consistent with what can be considered moderate damage. The cracked specimens were then repaired with an epoxy injection technique and subsequently retested until reaching failure. The epoxy injection techniques demonstrated to be quite efficient in partly, although not fully, restoring the energy dissipation capacities of the damaged specimens at all beam rotation levels. The stiffness was partly restored within the elastic range and almost fully restored after the onset of nonlinear behaviour.
This work investigates the possibility of developing a non-contact, non-line of sight sensor to measure interstorey drift through simulation and experimental validation. • The method uses frequency-modulated continuous wave (FMCW) radar to measure displacement. This method is commonly in use in a number of modern applications, including aircraft altimeters and automotive parking sensors. • The technique avoids numerous problems found in contemporary structural health monitoring methods, namely integral drift errors and structural modification requirements. • The smallest achievable detection error in displacement was found to be as low as 0.26%, through simulated against the displacement response of a single degree of freedom structure subject to ground motion excitation. • This was verified during experimentation, when a corner-style reflector was placed on a shake table running ground motion data taken from the 4th September 2010 earthquake in Christchurch. These results confirmed the conclusions drawn from simulation.
This article presents a quantitative case study on the site amplification effect observed at Heathcote Valley, New Zealand, during the 2010-2011 Canterbury earthquake sequence for 10 events that produced notable ground acceleration amplitudes up to 1.4g and 2.2g in the horizontal and vertical directions, respectively. We performed finite element analyses of the dynamic response of the valley, accounting for the realistic basin geometry and the soil non-linear response. The site-specific simulations performed significantly better than both empirical ground motion models and physics based regional-scale ground motion simulations (which empirically accounts for the site effects), reducing the spectral acceleration prediction bias by a factor of two in short vibration periods. However, our validation exercise demonstrated that it was necessary to quantify the level of uncertainty in the estimated bedrock motion using multiple recorded events, to understand how much the simplistic model can over- or under-estimate the ground motion intensities. Inferences from the analyses suggest that the Rayleigh waves generated near the basin edge contributed significantly to the observed high frequency (f>3Hz) amplification, in addition to the amplification caused by the strong soil-rock impedance contrast at the site fundamental frequency. Models with and without considering soil non-linear response illustrate, as expected, that the linear elastic assumption severely overestimates ground motions in high frequencies for strong earthquakes, especially when the contribution of basin edge-generated Rayleigh waves becomes significant. Our analyses also demonstrate that the effect of pressure-dependent soil velocities on the high frequency ground motions is as significant as the amplification caused by the basin edge-generated Rayleigh waves.
Nel presente articolo si illustra una procedura per il processamento automatizzato di prove CPT, il calcolo di vari indici di liquefazione e la rappresentazione dei dati su mappa. La procedura è applicata al caso studio del terremoto di Christchurch, Nuova Zelanda, del 22 febbraio 2011 (magnitudo momento, Mw = 6.2). Dall’analisi spaziale dei risultati emerge una buona correlazione tra le mappe ottenute per l’indicatore degli effetti al suolo e i danni osservati (su terreni e strutture). Tuttavia, per confermare la validità di tale procedura, sarà necessario esaminare ulteriori casi studio nel mondo.
This paper provides a brief discussion of observed strong ground motions from the 14 November 2016 Mw7.8 Kaikoura earthquake. Specific attention is given to examining observations in the near-source region where several ground motions exceeding 1.0g horizontal are recorded, as well as up to 2.7g in the vertical direction at one location. Ground motion response spectra in the near-source, North Canterbury, Marlborough and Wellington regions are also examined and compared with design levels. Observed spectral amplitudes are also compared with predictions from empirical and physics-based ground motion modelling.
None
We present preliminary observations on three waters impacts from the Mw7.8 14th November 2016 Kaikōura Earthquake on wider metropolitan Wellington, urban and rural Marlborough, and in Kaikōura township. Three waters systems in these areas experienced widespread and significant transient ground deformation in response to seismic shaking, with localised permanent ground deformation via liquefaction and lateral spreading. In Wellington, potable water quality was impacted temporarily by increased turbidity, and significant water losses occurred due to damaged pipes at the port. The Seaview and Porirua wastewater treatment plants sustained damage to clarifier tanks from water seiching, and increased water infiltration to the wastewater system occurred. Most failure modes in urban Marlborough were similar to the 2010-2011 Canterbury Earthquake Sequence; however some rural water tanks experienced rotational and translational movements, highlighting importance of flexible pipe connections. In Kaikōura, damage to reservoirs and pipes led to loss of water supply and compromised firefighting capability. Wastewater damage led to environmental contamination, and necessitated restrictions on greywater entry into the system to minimise flows. Damage to these systems necessitated the importation of tankered and bottled water, boil water notices and chlorination of the system, and importation of portaloos and chemical toilets. Stormwater infrastructure such as road drainage channels was also damaged, which could compromise condition of underlying road materials. Good operational asset management practices (current and accurate information, renewals, appreciation of criticality, good system knowledge and practical contingency plans) helped improve system resilience, and having robust emergency management centres and accurate Geographic Information System data allowed effective response coordination. Minimal damage to the wider built environment facilitated system inspections. Note Future research will include detailed geospatial assessments of seismic demand on these systems and attendant modes of failure, levels of service restoration, and collaborative development of resilience measures.
This report presents the simplified seismic assessment of a case study reinforced concrete (RC) building following the newly developed and refined NZSEE/MBIE guidelines on seismic assessment (NZSEE/MBIE, semi-final draft 26 October 2016). After an overview of the step-by-step ‘diagnostic’ process, including an holistic and qualitative description of the expected vulnerabilities and of the assessment strategy/methodology, focus is given, whilst not limited, to the implementation of a Detailed Seismic Assessment (DSA) (NZSEE/MBIE, 2016c). The DSA is intended to provide a more reliable and consistent outcome than what can be provided by an initial seismic assessment (ISA). In fact, while the Initial Seismic Assessment (ISA), of which the Initial Evaluation Procedure is only a part of, is the more natural and still recommended first step in the overall assessment process, it is mostly intended to be a coarse evaluation involving as few resources as reasonably possible. It is thus expected that an ISA will be followed by a Detailed Seismic Assessment (DSA) not only where the threshold of 33%NBS is not achieved but also where important decisions are intended that are reliant on the seismic status of the building. The use of %NBS (% New Building Standard) as a capacity/demand ratio to describe the result of the seismic assessment at all levels of assessment procedure (ISA through to DSA) is deliberate by the NZSEE/MBIE guidelines (Part A) (NZSEE/MBIE 2016a). The rating for the building needs only be based on the lowest level of assessment that is warranted for the particular circumstances. Discussion on how the %NBS rating is to be determined can be found in Section A3.3 (NZSEE/MBIE 2016a), and, more specifically, in Part B for the ISA (NZSEE/MBIE 2016b) and Part C for the DSA (NZSEE/MBIE 2016c). As per other international approaches, the DSA can be based on several analysis procedures to assess the structural behaviour (linear, nonlinear, static or dynamic, force or displacement-based). The significantly revamped NZSEE 2016 Seismic Assessment Guidelines strongly recommend the use of an analytical (basically ‘by hand’) method, referred to the Simple Lateral Mechanism Analysis (SLaMA) as a first phase of any other numerically-based analysis method. Significant effort has thus been dedicated to provide within the NZSEE 2016 guidelines (NZSEE/MBIE 2016c) a step-by-step description of the procedure, either in general terms (Chapter 2) or with specific reference to Reinforced Concrete Buildings (Chapter 5). More specifically, extract from the guidelines, NZSEE “recommend using the Simple Lateral Mechanism Analysis (SLaMA) procedure as a first step in any assessment. While SLaMA is essentially an analysis technique, it enables assessors to investigate (and present in a simple form) the potential contribution and interaction of a number of structural elements and their likely effect on the building’s global capacity. In some cases, the results of a SLaMA will only be indicative. However, it is expected that its use should help assessors achieve a more reliable outcome than if they only carried out a detailed analysis, especially if that analysis is limited to the elastic range For complex structural systems, a 3D dynamic analysis may be necessary to supplement the simplified nonlinear Simple Lateral Mechanism Analysis (SLaMA).” This report presents the development of a full design example for the the implementation of the SLaMA method on a case study buildings and a validation/comparison with a non-linear static (pushover) analysis. The step-by-step-procedure, summarized in Figure 1, will be herein demonstrated from a component level (beams, columns, wall elements) to a subassembly level (hierarchy of strength in a beam-column joint) and to a system level (frame, C-Wall) assuming initially a 2D behaviour of the key structural system, and then incorporating a by-hand 3D behaviour (torsional effects).
Decision making on the reinstatement of the Christchurch sewer system after the Canterbury (New Zealand) earthquake sequence in 2010–2011 relied strongly on damage data, in particular closed circuit television (CCTV). This paper documents that process and considers how data can influence decision making. Data are analyzed on 33,000 pipes and 13,000 repairs and renewals. The primary findings are that (1) there should be a threshold of damage per pipe set to make efficient use of CCTV; (2) for those who are estimating potential damage, care must be taken in direct use of repair data without an understanding of the actual damage modes; and (3) a strong correlation was found between the ratio of faults to repairs per pipe and the estimated peak ground velocity. Taken together, the results provide evidence of the extra benefit that damage data can provide over repair data for wastewater networks and may help guide others in the development of appropriate strategies for data collection and wastewater pipe decisions after disasters.
The focus of the study presented herein is an assessment of the relative efficacy of recent Cone Penetration Test (CPT) and small strain shear wave velocity (Vs) based variants of the simplified procedure. Towards this end Receiver Operating Characteristic (ROC) analyses were performed on the CPT- and Vs-based procedures using the field case history databases from which the respective procedures were developed. The ROC analyses show that Factors of Safety (FS) against liquefaction computed using the most recent Vs-based simplified procedure is better able to separate the “liquefaction” from the “no liquefaction” case histories in the Vs liquefaction database than the CPT-based procedure is able to separate the “liquefaction” from the “no liquefaction” case histories in the CPT liquefaction database. However, this finding somewhat contradicts the assessed predictive capabilities of the CPT- and Vs-based procedures as quantified using select, high quality liquefaction case histories from the 20102011 Canterbury, New Zealand, Earthquake Sequence (CES), wherein the CPT-based procedure was found to yield more accurate predictions. The dichotomy of these findings may result from the fact that different liquefaction field case history databases were used in the respective ROC analyses for Vs and CPT, while the same case histories were used to evaluate both the CPT- and Vs-based procedures.
The magnitude Mw7.8 ‘Kaikōura’ earthquake occurred shortly after midnight on 14 November 2016. This paper presents an overview of the geotechnical impacts on the South Island of New Zealand recorded during the postevent reconnaissance. Despite the large moment magnitude of this earthquake, relatively little liquefaction was observed across the South Island, with the only severe manifestation occurring in the young, loose alluvial deposits in the floodplains of the Wairau and Opaoa Rivers near Blenheim. The spatial extent and volume of liquefaction ejecta across South Island is significantly less than that observed in Christchurch during the 2010-2011 Canterbury Earthquake Sequence, and the impact of its occurrence to the built environment was largely negligible on account of the severe manifestations occurring away from the areas of major development. Large localised lateral displacements occurred in Kaikōura around Lyell Creek. The soft fine-grained material in the upper portions of the soil profile and the free face at the creek channel were responsible for the accumulation of displacement during the ground shaking. These movements had severely impacted the houses which were built close (within the zone of large displacement) to Lyell Creek. The wastewater treatment facility located just north of Kaikōura also suffered tears in the liners of the oxidation ponds and distortions in the aeration system due to ground movements. Ground failures on the Amuri and Emu Plains (within the Waiau Valley) were small considering the large peak accelerations (in excess of 1g) experienced in the area. Minor to moderate lateral spreading and ejecta was observed at some bridge crossings in the area. However, most of the structural damage sustained by the bridges was a result of the inertial loading, and the damage resulting from geotechnical issues were secondary.
At 00:02 on 14th November 2016, a Mw 7.8 earthquake occurred in and offshore of the northeast of the South Island of New Zealand. Fault rupture, ground shaking, liquefaction, and co-seismic landslides caused severe damage to distributed infrastructure, and particularly transportation networks; large segments of the country’s main highway, State Highway 1 (SH1), and the Main North Line (MNL) railway line, were damaged between Picton and Christchurch. The damage caused direct local impacts, including isolation of communities, and wider regional impacts, including disruption of supply chains. Adaptive measures have ensured immediate continued regional transport of goods and people. Air and sea transport increased quickly, both for emergency response and to ensure routine transport of goods. Road diversions have also allowed critical connections to remain operable. This effective response to regional transport challenges allowed Civil Defence Emergency Management to quickly prioritise access to isolated settlements, all of which had road access 23 days after the earthquake. However, 100 days after the earthquake, critical segments of SH1 and the MNL remain closed and their ongoing repairs are a serious national strategic, as well as local, concern. This paper presents the impacts on South Island transport infrastructure, and subsequent management through the emergency response and early recovery phases, during the first 100 days following the initial earthquake, and highlights lessons for transportation system resilience.
A 3D high-resolution model of the geologic structure and associated seismic velocities in the Canterbury, New Zealand region is developed utilising data from depthconverted seismic reflection lines, petroleum and water well logs, cone penetration tests, and implicitly guided by existing contour maps and geologic cross sections in data sparse subregions. The model, developed using geostatistical Kriging, explicitly represents the significant and regionally recognisable geologic surfaces that mark the boundaries between geologic units with distinct lithology and age. The model is examined in the form of both geologic surface elevation contour maps as well as vertical cross sections of shear wave velocity, with the most prominent features being the Banks Peninsula Miocene-Pliocene volcanic edifice, and the Pegasus and Rakaia late Mesozoic-Neogene sedimentary basins. The adequacy of the modelled geologic surfaces is assessed through a residual analysis of point constraints used in the Kriging and qualitative comparisons with previous geologic models of subsets of the region. Seismic velocities for the lithological units between the geologic surfaces have also been derived, thus providing the necessary information for a Canterbury velocity model (CantVM) for use in physics-based seismic wave propagation. The developed model also has application for the determination of depths to specified shear wave velocities for use in empirical ground motion modelling, which is explicitly discussed via an example.
This study uses 44 high quality liquefaction case histories taken from 22 locations affected by the 2010-2011 Canterbury earthquake sequence to evaluate four commonly used CPT-VS correlations (i.e., Robertson, 2009; Hegazy and Mayne, 2006; Andrus et al., 2007; McGann et al., 2015b). Co-located CPT soundings and VS profiles, developed from surface wave testing, were obtained at 22 locations and case histories were developed for the Mw 7.1, 4 September 2010 Darfield and Mw 6.2, 22 February 2011 Christchurch earthquakes. The CPT soundings are used to generate VS profiles using each of four CPT-VS correlations. These correlated VS profiles are used to estimate the factor of safety against liquefaction using the Kayen et al. (2013) VS-based simplified liquefaction evaluation procedure. An error index is used to quantify the predictive capabilities of these correlations in relation to the observations of liquefaction (or the lack thereof). Additionally, the error indices from the CPT-correlated VS profiles are compared to those obtained using: (1) the Kayen et al. (2013) procedure with surface wave-derived VS profiles, and (2) the Idriss and Boulanger (2008) CPT-based liquefaction evaluation procedure. Based on the error indices, the evaluation procedures based on direct measurements of either CPT or VS provided more accurate liquefaction triggering estimates than those obtained from any of the CPT-VS correlations. However, the performance of the CPT-VS correlations varied, with the Robertson (2009) and Hegazy and Mayne (2006) correlations performing relatively poorly for the Christchurch soils and the Andrus et al. (2007) and McGann et al. (2015b) correlations performing better. The McGann et al. (2015b) correlation had the lowest error indices of the CPT-VS correlations tested, however, none of the CPT-VS correlations provided accurate enough VS predictions to be used for the evaluation of liquefaction triggering using the VS-based liquefaction evaluation procedures.
This poster aims to present fragility functions for pipelines buried in liquefaction-prone soils. Existing fragility models used to quantify losses can be based on old data or use complex metrics. Addressing these issues, the proposed functions are based on the Christchurch network and soil and utilizes the Canterbury earthquake sequence (CES) data, partially represented in Figure 1. Figure 1 (a) presents the pipe failure dataset, which describes the date, location and pipe on which failures occurred. Figure 1 (b) shows the simulated ground motion intensity median of the 22nd February 2011 earthquake. To develop the model, the network and soil characteristics have also been utilized
The M7.8 Kaikoura Earthquake in 2016 presented a number of challenges to science agencies and institutions throughout New Zealand. The earthquake was complex, with 21 faults rupturing throughout the North Canterbury and Marlborough landscape, generating a localised seven metre tsunami and triggering thousands of landslides. With many areas isolated as a result, it presented science teams with logistical challenges as well as the need to coordinate efforts across institutional and disciplinary boundaries. Many research disciplines, from engineering and geophysics to social science, were heavily involved in the response. Coordinating these disciplines and institutions required significant effort to assist New Zealand during its most complex earthquake yet recorded. This paper explores that effort and acknowledges the successes and lessons learned by the teams involved.
SeisFinder is an open-source web service developed by QuakeCoRE and the University of Canterbury, focused on enabling the extraction of output data from computationally intensive earthquake resilience calculations. Currently, SeisFinder allows users to select historical or future events and retrieve ground motion simulation outputs for requested geographical locations. This data can be used as input for other resilience calculations, such as dynamic response history analysis. SeisFinder was developed using Django, a high-level python web framework, and uses a postgreSQL database. Because our large-scale computationally-intensive numerical ground motion simulations produce big data, the actual data is stored in file systems, while the metadata is stored in the database.
The Canterbury earthquakes of 2010 and 2011 caused significant damage and disruption to the city of Christchurch, New Zealand. A Royal Commission was established to report on the causes of building failure as a result of the earthquakes as well as look at the legal and best-practice requirements for buildings in New Zealand Central Business Districts. The Royal Commission made 189 recommendations on a variety of matters including managing damaged buildings after an earthquake, the adequacy of building codes and standards, and the processes of seismic assessments of existing buildings to determine their earthquake vulnerability. In response the Ministry of Business, Innovation and Employment, the agency responsible for administering building regulation in New Zealand, established a work programme to assist with the Canterbury rebuild and to implement the lessons learned throughout New Zealand. The five primary work streams in the programme are: • Facilitating the Canterbury Rebuild • Structural Performance and Design Standards • Geotechnical and structural guidance • Existing Building Resilience • Post Disaster Building Management This paper provides more detail on each of the work streams. There has been significant collaboration between the New Zealand Government and the research community, technical societies, and engineering consultants, both within New Zealand and internationally, to deliver the programme and improve the resilience of the New Zealand built environment. This has presented major challenges for an extremely busy industry in the aftermath of the Canterbury earthquakes. The paper identifies the items of work that have been completed and the work that is still in progress at the time of writing.
Christchurch City Council (Council) is undertaking the Land Drainage Recovery Programme in order to assess the effects of the earthquakes on flood risk to Christchurch. In the course of these investigations it has become better understood that floodplain management should be considered in a multi natural hazards context. Council have therefore engaged the Jacobs, Beca, University of Canterbury, and HR Wallingford project team to investigate the multihazards in eastern areas of Christchurch and develop flood management options which also consider other natural hazards in that context (i.e. how other hazards contribute to flooding both through temporal and spatial coincidence). The study has three stages: Stage 1 Gap Analysis – assessment of information known, identification of gaps and studies required to fill the gaps. Stage 2 Hazard Studies – a gap filling stage with the studies identified in Stage 1. Stage 3 Collating, Optioneering and Reporting – development of options to manage flood risk. This present report is to document findings of Stage 1 and recommends the studies that should be completed for Stage 2. It has also been important to consider how Stage 3 would be delivered and the gaps are prioritised to provide for this. The level of information available and hazards to consider is extensive; requiring this report to be made up of five parts each identifying individual gaps. A process of identifying information for individual hazards in Christchurch has been undertaken and documented (Part 1) followed by assessing the spatial co-location (Part 2) and probabilistic presence of multi hazards using available information. Part 3 considers multi hazard presence both as a temporal coincidence (e.g. an earthquake and flood occurring at one time) and as a cascade sequence (e.g. earthquake followed by a flood at some point in the future). Council have already undertaken a number of options studies for managing flood risk and these are documented in Part 4. Finally Part 5 provides the Gap Analysis Summary and Recommendations to Council. The key findings of Stage 1 gap analysis are: - The spatial analysis showed eastern Christchurch has a large number of hazards present with only 20% of the study area not being affected by any of the hazards mapped. Over 20% of the study area is exposed to four or more hazards at the frequencies and data available. - The majority of the Residential Red Zone is strongly exposed to multiple hazards, with 86% of the area being exposed to 4 or more hazards, and 24% being exposed to 6 or more hazards. - A wide number of gaps are present; however, prioritisation needs to consider the level of benefit and risks associated with not undertaking the studies. In light of this 10 studies ranging in scale are recommended to be done for the project team to complete the present scope of Stage 3. - Stage 3 will need to consider a number of engineering options to address hazards and compare with policy options; however, Council have not established a consistent policy on managed retreat that can be applied for equal comparison; without which substantial assumptions are required. We recommend Council undertake a study to define a managed retreat framework as an option for the city. - In undertaking Stage 1 with floodplain management as the focal point in a multi hazards context we have identified that Stage 3 requires consideration of options in the context of economics, implementation and residual risk. Presently the scope of work will provide a level of definition for floodplain options; however, this will not be at equal levels of detail for other hazard management options. Therefore, we recommend Council considers undertaking other studies with those key hazards (e.g. Coastal Hazards) as a focal point and identifies the engineering options to address such hazards. Doing so will provide equal levels of information for Council to make an informed and defendable decision on which options are progressed following Stage 3.
Recently developed performance-based earthquake engineering framework, such as one provided by PEER (Deierlein et al. 2003), assist in the quantification in terms of performance such as casualty, monetary losses and downtime. This opens up the opportunity to identify cost-effective retrofit/rehabilitation strategies by comparing upfront costs associated with retrofit with the repair costs that can be expected over time. This loss assessment can be strengthened by learning from recent earthquakes, such as the 2010 Canterbury and 2016 Kaikoura earthquakes. In order to investigate which types of retrofit/rehabilitation strategies may be most cost-effective, a case study building was chosen for this research. The Pacific Tower, a 22-storey EBF apartment located within the Christchurch central business district (CBD), was damaged and repaired during the 2010 Canterbury earthquake series. As such, by taking hazard levels accordingly (i.e. to correspond to the Christchurch CBD), modelling and analysing the structure, and considering the vulnerability and repair costs of its different components, it is possible to predict the expected losses of the aforementioned building. Using this information, cost-effective retrofit/rehabilitation strategy can be determined. This research found that more often than not, it would be beneficial to improve the performance of valuable non-structural components, such as partitions. Although it is true that improving such elements will increase the initial costs, over time, the benefits gained from reduced losses should be expected to overcome the initial costs. Aftershocks do increase the predicted losses of a building even in lower intensities due to the fact that non-structural components can get damaged at such low intensities. By comparing losses computed with and without consideration of aftershocks for a range of historical earthquakes, it was found that the ratio between losses due to main shock with aftershocks to the losses due to the main shock only tended to increase with increasing main shock magnitude. This may be due to the fact that larger magnitude earthquakes tend to generate larger magnitude aftershocks and as those aftershocks happen within a region around the main shock, they are more likely to cause intense shaking and additional damage. In addition to this observation, it was observed that the most significant component of loss of the case study building was the non-structural partition walls.
Non-structural elements (NSEs) have frequently proven to contribute to significant losses sustained from earthquakes in the form of damage, downtime, injury and death. In New Zealand (NZ), the 2010 and 2011 Canterbury Earthquake Sequence (CES), the 2013 Seddon and Cook Strait earthquake sequence and the 2016 Kaikoura earthquake were major milestones in this regard as significant damage to building NSEs both highlighted and further reinforced the importance of NSE seismic performance to the resilience of urban centres. Extensive damage in suspended ceilings, partition walls, façades and building services following the CES was reported to be partly due to erroneous seismic design or installation or caused by intervening elements. Moreover, the low-damage solutions developed for structural systems sometimes allow for relatively large inter-story drifts -compared to conventional designs- which may not have been considered in the seismic design of NSEs. Having observed these shortcomings, this study on suspended ceilings was carried out with five main goals: i) Understanding the seismic performance of the system commonly used in NZ; ii) Understanding the transfer of seismic design actions through different suspended ceiling components, iii) Investigating potential low-damage solutions; iii) Evaluating the compatibility of the current ceiling system with other low-damage NSEs; and iv) Investigating the application of numerical analysis to simulate the response of ceiling systems. The first phase of the study followed a joint research work between the University of Canterbury (UC) in NZ, and the Politecnico Di Milano, in Italy. The experimental ceiling component fragility curves obtained in this existing study were employed to produce analytical fragility curves for a perimeter-fixed ceiling of a given size and weight, with grid acceleration as the intensity measure. The validity of the method was proven through comparisons between this proposed analytical approach with the recommended procedures in proprietary products design guidelines, as well as experimental fragility curves from other studies. For application to engineering design practice, and using fragility curves for a range of ceiling lengths and weights, design curves were produced for estimating the allowable grid lengths for a given demand level. In the second phase of this study, three specimens of perimeter-fixed ceilings were tested on a shake table under both sinusoidal and random floor motion input. The experiments considered the relationship between the floor acceleration, acceleration of the ceiling grid, the axial force induced in the grid members, and the effect of boundary conditions on the transfer of these axial forces. A direct correlation was observed between the axial force (recorded via load cells) and the horizontal acceleration measured on the ceiling grid. Moreover, the amplification of floor acceleration, as transferred through ceiling components, was examined and found (in several tests) to be greater than the recommended factor for the design of ceilings provided in the NZ earthquake loadings standard NZS1170.5. However, this amplification was found to be influenced by the pounding interactions between the ceiling grid members and the tiles, and this amplification diminished considerably when the high frequency content was filtered out from the output time histories. The experiments ended with damage in the ceiling grid connection at an axial force similar to the capacity of these joints previously measured through static tests in phase one. The observation of common forms of damage in ceilings in earthquakes triggered the monotonic experiments carried out in the third phase of this research with the objective of investigating a simple and easily applicable mitigation strategy for existing or new suspended ceilings. The tests focused on the possibility of using proprietary cross-shaped clip elements ordinarily used to provide seismic gap as a strengthening solution for the weak components of a ceiling. The results showed that the solution was effective under both tension and compression loads through increasing load bearing capacity and ductility in grid connections. The feasibility of a novel type of suspended ceiling called fully-floating ceiling system was investigated through shaking table tests in the next phase of this study with the main goal of isolating the ceiling from the surrounding structure; thereby arresting the transfer of associated seismic forces from the structure to the ceiling. The fully-floating ceiling specimen was freely hung from the floor above lacking any lateral bracing and connections with the perimeter. Throughout different tests, a satisfactory agreement between the fully-floating ceiling response and simple pendulum theory was demonstrated. The addition of isolation material in perimeter gaps was found effective in inducing extra damping and protecting the ceiling from pounding impact; resulting in much reduced ceiling displacements and accelerations. The only form of damage observed throughout the random floor motion tests and the sinusoidal tests was a panel dislodgement observed in a test due to successive poundings between the ceiling specimen and the surrounding beams at resonant frequencies. Partition walls as the first effective NSE in direct interaction with ceilings were the topic of the final experimental phase. Low-damage drywall partitions proposed in a previous study in the UC were tested with two common forms of suspended ceiling: braced and perimeter-fixed. The experiments investigated the in-plane and out-of-plane performance of the low-damage drywall partitions, as well as displacement compatibility between these walls and the suspended ceilings. In the braced ceiling experiment, where no connection was made between ceiling grids and surrounding walls no damage in the grid system or partitions was observed. However, at high drift values panel dislodgement was observed on corners of the ceiling where the free ends of grids were not restrained against spreading. This could be prevented by framing the grid ends using a perimeter angle that is riveted only to the grid members while keeping sufficient clearance from the perimeter walls. In the next set of tests with the perimeter-fixed ceiling, no damage was observed in the ceiling system or the drywalls. Based on the results of the experiments it was concluded that the tested ceiling had enough flexibility to accommodate the relative displacement between two perpendicular walls up to the inter-storey drifts achieved. The experiments on perimeter-fixed ceilings were followed by numerical simulations of the performance of these ceilings in a finite element model developed in the structural analysis software, SAP2000. This model was relatively simple and easy to develop and was able to replicate the experimental results to a reasonable degree. Filtering was applied to the experimental output to exclude the effect of high frequency noise and tile-grid impact. The developed model generally simulated the acceleration responses well but underestimated the peak ceiling grid accelerations. This was possibly because the peak values in time histories were affected by impact occurring at very short periods. The model overestimated the axial forces in ceiling grids which was assumed to be caused by the initial assumptions made about the tributary area or constant acceleration associated with each grid line in the direction of excitation. Otherwise, the overall success of the numerical modelling in replicating the experimental results implies that numerical modelling using conventional structural analysis software could be used in engineering practice to analyse alternative ceiling geometries proposed for application to varying structural systems. This however, needs to be confirmed through similar analyses on other ceiling examples from existing instrumented buildings during real earthquakes. As the concluding part of this research the final phase addressed the issues raised following the review of existing ceiling standards and guidelines. The applicability of the research findings to current practice and their implications were discussed. Finally, an example was provided for the design of a suspended ceiling utilising the new knowledge acquired in this research.
This dissertation addresses several fundamental and applied aspects of ground motion selection for seismic response analyses. In particular, the following topics are addressed: the theory and application of ground motion selection for scenario earthquake ruptures; the consideration of causal parameter bounds in ground motion selection; ground motion selection in the near-fault region where directivity effect is significant; and methodologies for epistemic uncertainty consideration and propagation in the context of ground motion selection and seismic performance assessment. The paragraphs below outline each contribution in more detail. A scenario-based ground motion selection method is presented which considers the joint distribution of multiple intensity measure (IM) types based on the generalised conditional intensity measure (GCIM) methodology (Bradley, 2010b, 2012c). The ground motion selection algorithm is based on generating realisations of the considered IM distributions for a specific rupture scenario and then finding the prospective ground motions which best fit the realisations using an optimal amplitude scaling factor. In addition, using different rupture scenarios and site conditions, two important aspects of the GCIM methodology are scrutinised: (i) different weight vectors for the various IMs considered; and (ii) quantifying the importance of replicate selections for ensembles with different numbers of desired ground motions. As an application of the developed scenario-based ground motion selection method, ground motion ensembles are selected to represent several major earthquake scenarios in New Zealand that pose a significant seismic hazard, namely, Alpine, Hope and Porters Pass ruptures for Christchurch city; and Wellington, Ohariu, and Wairarapa ruptures for Wellington city. A rigorous basis is developed, and sensitivity analyses performed, for the consideration of bounds on causal parameters (e.g., magnitude, source-to-site distance, and site condition) for ground motion selection. The effect of causal parameter bound selection on both the number of available prospective ground motions from an initial empirical as-recorded database, and the statistical properties of IMs of selected ground motions are examined. It is also demonstrated that using causal parameter bounds is not a reliable approach to implicitly account for ground motion duration and cumulative effects when selection is based on only spectral acceleration (SA) ordinates. Specific causal parameter bounding criteria are recommended for general use as a ‘default’ bounding criterion with possible adjustments from the analyst based on problem-specific preferences. An approach is presented to consider the forward directivity effects in seismic hazard analysis, which does not separate the hazard calculations for pulse-like and non-pulse-like ground motions. Also, the ability of ground motion selection methods to appropriately select records containing forward directivity pulse motions in the near-fault region is examined. Particular attention is given to ground motion selection which is explicitly based on ground motion IMs, including SA, duration, and cumulative measures; rather than a focus on implicit parameters (i.e., distance, and pulse or non-pulse classifications) that are conventionally used to heuristically distinguish between the near-fault and far-field records. No ad hoc criteria, in terms of the number of directivity ground motions and their pulse periods, are enforced for selecting pulse-like records. Example applications are presented with different rupture characteristics, source-to-site geometry, and site conditions. It is advocated that the selection of ground motions in the near-fault region based on IM properties alone is preferred to that in which the proportion of pulse-like motions and their pulse periods are specified a priori as strict criteria for ground motion selection. Three methods are presented to propagate the effect of seismic hazard and ground motion selection epistemic uncertainties to seismic performance metrics. These methods differ in their level of rigor considered to propagate the epistemic uncertainty in the conditional distribution of IMs utilised in ground motion selection, selected ground motion ensembles, and the number of nonlinear response history analyses performed to obtain the distribution of engineering demand parameters. These methods are compared for an example site where it is observed that, for seismic demand levels below the collapse limit, epistemic uncertainty in ground motion selection is a smaller uncertainty contributor relative to the uncertainty in the seismic hazard itself. In contrast, uncertainty in ground motion selection process increases the uncertainty in the seismic demand hazard for near-collapse demand levels.