In this thesis, focus is given to develop methodologies for rapidly estimating specific components of loss and downtime functions. The thesis proposes methodologies for deriving loss functions by (i) considering individual component performance; (ii) grouping them as per their performance characteristics; and (iii) applying them to similar building usage categories. The degree of variation in building stock and understanding their characteristics are important factors to be considered in the loss estimation methodology and the field surveys carried out to collect data add value to the study. To facilitate developing ‘downtime’ functions, this study investigates two key components of downtime: (i) time delay from post-event damage assessment of properties; and (ii) time delay in settling the insurance claims lodged. In these two areas, this research enables understanding of critical factors that influence certain aspects of downtime and suggests approaches to quantify those factors. By scrutinising the residential damage insurance claims data provided by the Earthquake Commission (EQC) for the 2010- 2011 Canterbury Earthquake Sequence (CES), this work provides insights into various processes of claims settlement, the time taken to complete them and the EQC loss contributions to building stock in Christchurch city and Canterbury region. The study has shown diligence in investigating the EQC insurance claim data obtained from the CES to get new insights and build confidence in the models developed and the results generated. The first stage of this research develops contribution functions (probabilistic relationships between the expected losses for a wide range of building components and the building’s maximum response) for common types of claddings used in New Zealand buildings combining the probabilistic density functions (developed using the quantity of claddings measured from Christchurch buildings), fragility functions (obtained from the published literature) and cost functions (developed based on inputs from builders) through Monte Carlo simulations. From the developed contribution functions, glazing, masonry veneer, monolithic and precast concrete cladding systems are found to incur 50% loss at inter-storey drift levels equal to 0.027, 0.003, 0.005 and 0.011, respectively. Further, the maximum expected cladding loss for glazing, masonry veneer, monolithic, precast concrete cladding systems are found to be 368.2, 331.9, 365.0, and 136.2 NZD per square meter of floor area, respectively. In the second stage of this research, a detailed cost breakdown of typical buildings designed and built for different purposes is conducted. The contributions of structural and non- structural components to the total building cost are compared for buildings of different usages, and based on the similar ratios of non-structural performance group costs to the structural performance group cost, four-building groups are identified; (i) Structural components dominant group: outdoor sports, stadiums, parkings and long-span warehouses, (ii) non- structural drift-sensitive components dominant group: houses, single-storey suburban buildings (all usages), theatres/halls, workshops and clubhouses, (iii) non-structural acceleration- sensitive components dominant group: hospitals, research labs, museums and retail/cold stores, and (iv) apartments, hotels, offices, industrials, indoor sports, classrooms, devotionals and aquariums. By statistically analysing the cost breakdowns, performance group weighting factors are proposed for structural, and acceleration-sensitive and drift-sensitive non-structural components for all four building groups. Thus proposed building usage groupings and corresponding weighting factors facilitate rapid seismic loss estimation of any type of building given the EDPs at storey levels are known. A model for the quantification of post-earthquake inspection duration is developed in the third stage of this research. Herein, phase durations for the three assessment phases (one rapid impact and two rapid building) are computed using the number of buildings needing inspections, the number of engineers involved in inspections and a phase duration coefficient (which considers the median building inspection time, efficiency of engineer and the number of engineers involved in each assessment teams). The proposed model can be used: (i) by national/regional authorities to decide the length of the emergency period following a major earthquake, and estimate the number of engineers required to conduct a post-earthquake inspection within the desired emergency period, and (ii) to quantify the delay due to inspection for the downtime modelling framework. The final stage of this research investigates the repair costs and insurance claim settlement time for damaged residential buildings in the 2010-2011 Canterbury earthquake sequence. Based on the EQC claim settlement process, claims are categorized into three groups; (i) Small Claims: claims less than NZD15,000 which were settled through cash payment, (ii) Medium Claims: claims less than NZD100,000 which were managed through Canterbury Home Repair Programme (CHRP), and (iii) Large Claims: claims above NZD100,000 which were managed by an insurance provider. The regional loss ratio (RLR) for greater Christchurch for three events inducing shakings of approximate seismic intensities 6, 7, and 8 are found to be 0.013, 0.066, and 0.171, respectively. Furthermore, the claim duration (time between an event and the claim lodgement date), assessment duration (time between the claim lodgement day and the most recent assessment day), and repair duration (time between the most recent assessment day and the repair completion day) for the insured residential buildings in the region affected by the Canterbury earthquake sequence is found to be in the range of 0.5-4 weeks, 1.5- 5 months, and 1-3 years, respectively. The results of this phase will provide useful information to earthquake engineering researchers working on seismic risk/loss and insurance modelling.
Low Damage Seismic Design (LDSD) guidance material being developed by Engineering NZ is considering a design drift limit for multi-storey buildings of 0.5% at a new damage control limit state (DCLS). The impact of this new design requirement on the expected annual loss due to repair costs is investigated for a four-storey office building with reinforced concrete walls located in Christchurch. The LDSD guidance material aims to reduce the expected annual loss of complying buildings to below 0.1% of building replacement cost. The research tested this expectation. Losses were estimated in accordance with FEMA P58, using building responses from non-linear time history analyses (performed with OpenSees using lumped plasticity models). The equivalent static method, in line with NZS 1170.5 and NZS 3101, was used to design the building to LDSD specifications, representing a future state-of-practice design. The building designed to low-damage specification returned an expected annual loss of 0.10%, and the building designed conventionally returned an expected annual loss of 0.13%. Limitations with the NZS 3101 method for determining wall stiffness were identified, and a different method acknowledging the relationship between strength and stiffness was used to redesign the building. Along with improving this design assumption, the study finds that LDSD design criteria could be an effective way of limiting damage and losses.
"You displayed great fortitude in the face of sudden, overwhelming loss."
A pdf transcript of Ian's second earthquake story, captured by the UC QuakeBox Take 2 project. Interviewer: Samuel Hope. Transcriber: Josie Hepburn.
This thesis presents the application of data science techniques, especially machine learning, for the development of seismic damage and loss prediction models for residential buildings. Current post-earthquake building damage evaluation forms are developed for a particular country in mind. The lack of consistency hinders the comparison of building damage between different regions. A new paper form has been developed to address the need for a global universal methodology for post-earthquake building damage assessment. The form was successfully trialled in the street ‘La Morena’ in Mexico City following the 2017 Puebla earthquake. Aside from developing a framework for better input data for performance based earthquake engineering, this project also extended current techniques to derive insights from post-earthquake observations. Machine learning (ML) was applied to seismic damage data of residential buildings in Mexico City following the 2017 Puebla earthquake and in Christchurch following the 2010-2011 Canterbury earthquake sequence (CES). The experience showcased that it is readily possible to develop empirical data only driven models that can successfully identify key damage drivers and hidden underlying correlations without prior engineering knowledge. With adequate maintenance, such models have the potential to be rapidly and easily updated to allow improved damage and loss prediction accuracy and greater ability for models to be generalised. For ML models developed for the key events of the CES, the model trained using data from the 22 February 2011 event generalised the best for loss prediction. This is thought to be because of the large number of instances available for this event and the relatively limited class imbalance between the categories of the target attribute. For the CES, ML highlighted the importance of peak ground acceleration (PGA), building age, building size, liquefaction occurrence, and soil conditions as main factors which affected the losses in residential buildings in Christchurch. ML also highlighted the influence of liquefaction on the buildings losses related to the 22 February 2011 event. Further to the ML model development, the application of post-hoc methodologies was shown to be an effective way to derive insights for ML algorithms that are not intrinsically interpretable. Overall, these provide a basis for the development of ‘greybox’ ML models.
This study analyses the Earthquake Commission’s (EQC) insurance claims database to investigate the influence of seismic intensity and property damage resulting from the Canterbury Earthquake Sequence (CES) on the repair costs and claim settlement duration for residential buildings. Firstly, the ratio of building repair cost to its replacement cost was expressed as a Building Loss Ratio (BLR), which was further extended to Regional Loss Ratio (RLR) for greater Christchurch by multiplying the average of all building loss ratios with the proportion of building stock that lodged an insurance claim. Secondly, the total time required to settle the claim and the time taken to complete each phase of the claim settlement process were obtained. Based on the database, the regional loss ratio for greater Christchurch for three events producing shakings of intensities 6, 7, and 8 on the modified Mercalli intensity scale were 0.013, 0.066, and 0.171, respectively. Furthermore, small (less than NZD15,000), medium (between NZD15,000 and NZD100,000), and large (more than NZD100,000) claims took 0.35-0.55, 1.95-2.45, and 3.35-3.85 years to settle regardless of the building’s construction period and earthquake intensities. The number of claims was also disaggregated by various building characteristics to evaluate their relative contribution to the damage and repair costs.
Glazing systems are non-structural elements in a building that, more often than not, appear to be given little consideration in seismic design. Recent experimental work into glazing systems at the University of Canterbury, however, has shown that glazing systems can be very susceptible to serviceability damage, defined as loss of water-tightness. The focus of this paper is to highlight the difference in vulnerability of standard and seismic glazing systems and consider the implications of this for future repair costs and losses. The paper first describes the damage states chosen for glazing units according to the repair strategies required and expected repair costs. This includes three damage states: DS1: Water Leakage, DS2: Gasket Failure and DS3: Frame/Glass Failure. Implementing modern performance-based earthquake engineering, the paper proceeds to highlight a case study comparing costs and expected losses of a standard glazing unit and a seismic glazing unit installed on a case study building. It is shown that the use of seismic glazing units is generally beneficial over time, due to the early onset of serviceability damage in standard glazing units. Finally, the paper provides suggestions for designers aimed at reducing costs related to earthquake induced repairs of glazing.
A pdf transcript of Sara Green's second earthquake story, captured by the UC QuakeBox Take 2 project. Interviewer: Laura Moir. Transcriber: Sarah Woodfield.
A pdf transcript of John's second earthquake story, captured by the UC QuakeBox Take 2 project. Interviewer: Samuel Hope. Transcriber: Natalie Looyer.
A pdf transcript of Alvin Wade's second earthquake story, captured by the UC QuakeBox Take 2 project. Interviewer: Joshua Black. Transcriber: Josie Hepburn.
A pdf transcript of Betty and Michael's second earthquake story, captured by the UC QuakeBox Take 2 project. Interviewer: Samuel Hope. Transcriber: Sarah Woodfield.
A pdf transcript of Part 1 of Tracey Waiariki's second earthquake story, captured by the UC QuakeBox Take 2 project. Interviewer: Lucy Denham. Transcriber: Lucy Denham.
A pdf transcript of Vic Bartley's second earthquake story, captured by the UC QuakeBox Take 2 project. Interviewer: Samuel Hope. Transcriber: Sarah Woodfield.
A pdf transcript of Bev McCashin's second earthquake story, captured by the UC QuakeBox Take 2 project. The interview was conducted via Zoom. Interviewer: Laura Moir. Transcriber: Lauren Millar.
As damage and loss caused by natural hazards have increased worldwide over the past several decades, it is important for governments and aid agencies to have tools that enable effective post-disaster livelihood recovery to create self-sufficiency for the affected population. This study introduces a framework of critical components that constitute livelihood recovery and the critical factors that lead to people’s livelihood recovery. A comparative case study is employed in this research, combined with questionnaire surveys and interviews with those communities affected by large earthquakes in Lushan, China and in Christchurch and Kaikōura, New Zealand. In Lushan, China, a framework with four livelihood components was established, namely, housing, employment, wellbeing and external assistance. Respondents considered recovery of their housing to be the most essential element for livelihood diversification. External assistance was also rated highly in assisting with their livelihood recovery. Family ties and social connections seemed to have played a larger role than that of government agencies and NGOs. However, the recovery of livelihood cannot be fully achieved without wellbeing aspects being taken into account, and people believed that quality of life and their physical and mental health were essential for livelihood restoration. In Christchurch, New Zealand, the identified livelihood components were validated through in-depth interviews. The results showed that the above framework presenting what constitutes successful livelihood recovery could also be applied in Christchurch. This study also identified the critical factors to affect livelihood recovery following the Lushan and Kaikōura earthquakes, and these include community safety, availability of family support, level of community cohesion, long-term livelihood support, external housing recovery support, level of housing recovery and availability of health and wellbeing support. The framework developed will provide guidance for policy makers and aid agencies to prioritise their strategies and initiatives in assisting people to reinstate their livelihood in a timely manner post-disaster. It will also assist the policy makers and practitioners in China and New Zealand by setting an agenda for preparing for livelihood recovery in non-urgent times so the economic impact and livelihood disruption of those affected can be effectively mitigated.
Floor systems with precast concrete hollow-core units have been largely used in concrete buildings built in New Zealand during the 1980’s. Recent earthquakes, such as the Canterbury sequence in 2010-2011 and the Kaikoura earthquake in 2016, highlighted that this floor system can be highly vulnerable and potentially lead to the floor collapse. A series of research activities are in progress to better understand the seismic performance of floor diaphragms, and this research focuses on examining the performance of hollow core units running parallel to the walls of wall-resisting concrete structures. This study first focused on the development of fragility functions, which can be quickly used to assess likelihood of the hollow-core being able to survive given the buildings design drift, and secondly to determine the expected performance of hollow-core units that run parallel to walls, focusing on the alpha unit running by the wall. Fragility functions are created for a range of different parameters for both vertical dislocation and crack width that can be used as the basis of a quick analysis or loss estimation for the likely impact of hollow-core floors on building vulnerability and risk. This was done using past experimental tests, and the recorded damage. Using these results and the method developed by Baker fragility curves were able to be created for varying crack widths and vertical dislocations. Current guidelines for analysis of hollow-core unit incompatible displacements are based on experimental vertical displacement results from concrete moment resisting frame systems to determine the capacity of hollow-core elements. To investigate the demands on hollow-core units in a wall-based structure, a fibre-element model in the software Seismostruct is created and subject to quasi-static cyclic loading, using elements which are verified from previous experimental tests. It is shown that for hollow-core units running by walls that the 10 mm displacement capacity used for hollow-core units running by a beam is insufficient for members running by walls and that shear analysis should be used. The fibre-element model is used to simulate the seismic demand induced on the floor system and has shown that the shear demand is a function of drift, wall length, hollow-core span, linking slab length and, to a minor extent, wall elongation.
In major seismic events, a number of plan-asymmetric buildings which experienced element failure or structural collapse had twisted significantly about their vertical axis during the earthquake shaking. This twist, known as “building torsion”, results in greater demands on one side of a structure than on the other side. The Canterbury Earthquakes Royal Commission’s reports describe the response of a number of buildings in the February 2011 Christchurch earthquakes. As a result of the catastrophic collapse of one multi-storey building with significant torsional irregularity, and significant torsional effects also in other buildings, the Royal Commission recommended that further studies be undertaken to develop improved simple and effective guides to consider torsional effects in buildings which respond inelastically during earthquake shaking. Separately from this, as building owners, the government, and other stakeholders, are planning for possible earthquake scenarios, they need good estimates of the likely performance of both new and existing buildings. These estimates, often made using performance based earthquake engineering considerations and loss estimation techniques, inform decision making. Since all buildings may experience torsion to some extent, and torsional effects can influence demands on building structural and non-structural elements, it is crucial that demand estimates consider torsion. Building seismic response considering torsion can be evaluated with nonlinear time history analysis. However, such analysis involves significant computational effort, expertise and cost. Therefore, from an engineers’ point of view, simpler analysis methods, with reasonable accuracy, are beneficial. The consideration of torsion in simple analysis methods has been investigated by many researchers. However, many studies are theoretical without direct relevance to structural design/assessment. Some existing methods also have limited applicability, or they are difficult to use in routine design office practice. In addition, there has been no consensus about which method is best. As a result, there is a notable lack of recommendations in current building design codes for torsion of buildings that respond inelastically. There is a need for building torsion to be considered in yielding structures, and for simple guidance to be developed and adopted into building design standards. This study aims to undertaken to address this need for plan-asymmetric structures which are regular over their height. Time history analyses are first conducted to quantify the effects of building plan irregularity, that lead to torsional response, on the seismic response of building structures. Effects of some key structural and ground motion characteristics (e.g. hysteretic model, ground motion duration, etc.) are considered. Mass eccentricity is found to result in rather smaller torsional response compared to stiffness/strength eccentricity. Mass rotational inertia generally decreases the torsional response; however, the trend is not clearly defined for torsionally restrained systems (i.e. large λty). Systems with EPP and bilinear models have close displacements and systems with Takeda, SINA, and flag-shaped models yield almost the same displacements. Damping has no specific effect on the torsional response for the single-storey systems with the unidirectional eccentricity and excitation. Displacements of the single-storey systems subject to long duration ground motion records are smaller than those for short duration records. A method to consider torsional response of ductile building structures under earthquake shaking is then developed based on structural dynamics for a wide range of structural systems and configurations, including those with low and high torsional restraint. The method is then simplified for use in engineering practice. A novel method is also proposed to simply account for the effects of strength eccentricity on response of highly inelastic systems. A comparison of the accuracy of some existing methods (including code-base equivalent static method and model response spectrum analysis method), and the proposed method, is conducted for single-storey structures. It is shown that the proposed method generally provides better accuracy over a wide range of parameters. In general, the equivalent static method is not adequate in capturing the torsional effects and the elastic modal response spectrum analysis method is generally adequate for some common parameters. Record-to-record variation in maximum displacement demand on the structures with different degrees of torsional response is considered in a simple way. Bidirectional torsional response is then considered. Bidirectional eccentricity and excitation has varying effects on the torsional response; however, it generally increases the weak and strong edges displacements. The proposed method is then generalized to consider the bidirectional torsion due to bidirectional stiffness/strength eccentricity and bidirectional seismic excitation. The method is shown to predict displacements conservatively; however, the conservatism decreases slightly for cases with bidirectional excitation compared to those subject to unidirectional excitation. In is shown that the roof displacement of multi-storey structures with torsional response can be predicted by considering the first mode of vibration. The method is then further generalized to estimate torsional effects on multi-storey structure displacement demands. The proposed procedure is tested multi-storey structures and shown to predict the displacements with a good accuracy and conservatively. For buildings which twist in plan during earthquake shaking, the effect of P-Δλ action is evaluated and recommendations for design are made. P-Δλ has more significant effects on systems with small post- yield stiffness. Therefore, system stability coefficient is shown not to be the best indicator of the importance of P-Δλ and it is recommended to use post-yield stiffness of system computed with allowance for P-Δλ effects. For systems with torsional response, the global system stability coefficient and post- yield stiffness ration do not reflect the significance of P-Δλ effects properly. Therefore, for torsional systems individual seismic force resisting systems should be considered. Accuracy of MRSA is investigated and it is found that the MRSA is not always conservative for estimating the centre of mass and strong edge displacements as well as displacements of ductile systems with strength eccentricity larger than stiffness eccentricity. Some modifications are proposed to get the MRSA yields a conservative estimation of displacement demands for all cases.
