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.
A building boom in the 1980s allowed pre-stressed hollow-core floor construction to be widely adopted in New Zealand, even though the behaviour of these prefabricated elements within buildings was still uncertain. Inspections following the Canterbury and Kaikōura earthquakes has provided evidence of web-splitting, transverse cracking and longitudinal splitting on hollow-core units, confirming the susceptibility of these floors to undesirable failure modes. Hollow-core slabs are mainly designed to resist bending and shear. However, there are many applications in which they are also subjected to torsion. In New Zealand, hollow-core units contain no transverse reinforcement in the soffit concrete below the cells and no web reinforcement. Consequently, their dependable performance in torsion is limited to actions that they can resist before torsional cracking occurs. In previous work by the present authors, a three-dimensional FE modelling approach to study the shear flexural behaviour of precast pre-stressed hollow core units was developed and validated by full-scale experiments. This paper shows how the FE analyses have been extended to investigate the response of HC units subjected to torsional actions. Constitutive models, based on nonlinear fracture mechanics, have been used to numerically predict the torsional capacity of HC units and have been compared with experimental results. The results indicate that the numerical approach is promising and should be developed further as part of future research.
Observations of out-of-plane (OOP) instability in the 2010 Chile earthquake and in the 2011 Christchurch earthquake resulted in concerns about the current design provisions of structural walls. This mode of failure was previously observed in the experimental response of some wall specimens subjected to in-plane loading. Therefore, the postulations proposed for prediction of the limit states corresponding to OOP instability of rectangular walls are generally based on stability analysis under in-plane loading only. These approaches address stability of a cracked wall section when subjected to compression, thereby considering the level of residual strain developed in the reinforcement as the parameter that prevents timely crack closure of the wall section and induces stability failure. The New Zealand code requirements addressing the OOP instability of structural walls are based on the assumptions used in the literature and the analytical methods proposed for mathematical determination of the critical strain values. In this study, a parametric study is conducted using a numerical model capable of simulating OOP instability of rectangular walls to evaluate sensitivity of the OOP response of rectangular walls to variation of different parameters identified to be governing this failure mechanism. The effects of wall slenderness (unsupported height-to-thickness) ratio, longitudinal reinforcement ratio of the boundary regions and length on the OOP response of walls are evaluated. A clear trend was observed regarding the influence of these parameters on the initiation of OOP displacement, based on which simple equations are proposed for prediction of OOP instability in rectangular walls.
Major earthquakes, such as the Canterbury and Kaikoura events recorded in New Zealand in 2010 and 2016 respectively, highlighted that floor systems can be heavily damaged. At a reduced or full scale, quasi-static experimental tests on structural sub-assemblies can help to establish the seismic performance of structural systems. However, the experimental performance obtained with such tests is likely to be dependent on the drift protocol adopted. This paper provides an overview of the drift protocols which have been assumed in previous relevant experimental activities, with emphasis on those adopted for testing floor systems. The paper also describes the procedure used to define the loading protocol applied in the testing of a large precast concrete floor diaphragm as part of the Recast floor project at the University of Canterbury. Finally, major limits of current loading protocols, and areas of future research, are identified.
To reduce seismic vulnerability and the economic impact of seismic structural damage, it is important to protect structures using supplemental energy dissipation devices. Several types of supplemental damping systems can limit loads transferred to structures and absorb significant response energy without sacrificial structural damage. Lead extrusion dampers are one type of supplemental energy dissipation devices. A smaller volumetric size with high force capacities, called high force to volume (HF2V) devices, have been employed in a large series of scaled and full-scaled experiments, as well as in three new structures in Christchurch and San Francisco. HF2V devices have previously been designed using very simple models with limited precision. They are then manufactured, and tested to ensure force capacities match design goals, potentially necessitating reassembly or redesign if there is large error. In particular, devices with a force capacity well above or below a design range can require more testing and redesign, leading to increased economic and time cost. Thus, there is a major need for a modelling methodology to accurately estimate the range of possible device force capacity values in the design phase – upper and lower bounds. Upper and lower bound force capacity estimates are developed from equations in the metal extrusion literature. These equations consider both friction and extrusion forces between the lead and the bulged shaft in HF2V devices. The equations for the lower and upper bounds are strictly functions of device design parameters ensuring easy use in the design phase. Two different sets of estimates are created, leading to estimates for the lower and upper bounds denoted FLB,1, FUB,1, FUB,2, respectively. The models are validated by comparing the bounds with experimental force capacity data from 15 experimental HF2V device tests. All lower bound estimates are below or almost equal to the experimental device forces, and all upper bound estimates are above. Per the derivation, the (FLB,1, FUB,1) pair provide narrower bounds. The (FLB,1, FUB,1) pair also had a mean lower bound gap of -34%, meaning the lower bound was 74% of device force on average, while the mean upper bound gap for FUB,1 was +23%. These are relatively tight bounds, within ~±2 SE of device manufacture, and can be used as a guide to ensure device forces are in range for the actual design use when manufactured. Therefore, they provide a useful design tool.
Supplemental energy dissipation devices are increasingly used to protect structures, limit loads transferred to structural elements and absorbing significant response energy without sacrificial structural damage. Lead extrusion dampers are supplemental energy dissipation devices, where recent development of smaller volumetric size with high force capacities, called high force to volume (HF2V) devices, has seen deployment in a large series of scaled and full-scaled experiments, as well as in three new structures in Christchurch, NZ and San Francisco, USA. HF2V devices have previously been designed using limited precision models, so there is variation in force prediction capability. Further, while the overall resistive force is predicted, the knowledge of the relative contributions of the different internal reaction mechanisms to these overall resistive forces is lacking, limiting insight and predictive accuracy in device design. There is thus a major need for detailed design models to better understand force generation, and to aid precision device design. These outcomes would speed the overall design and implementation process for uptake and use, reducing the need for iterative experimental testing. Design parameters from 17 experimental HF2V device tests are used to create finite element models using ABAQUS. The analysis is run using ABAQUS Explicit, in multiple step times of 1 second with automatic increments, to balance higher accuracy and computational time. The output is obtained from the time- history output of the contact pressure forces including the normal and friction forces on the lead along the shaft. These values are used to calculate the resistive force on the shaft as it moves through the lead, and thus the device force. Results of these highly nonlinear, high strain analyses are compared to experimental device force results. Model errors compared to experimental results for all 17 devices ranged from 0% to 20% with a mean absolute error of 6.4%, indicating most errors were small. In particular, the standard error in manufacturing is SE = ±14%. In this case, 15 of 17 devices (88%) are within ±1SE (±14%) and 2 of 17 devices (12%) are within ±2SE (±28). These results show low errors and a distribution of errors compared to experimental results that are within experimental device construction variability. The overall modelling methodology is objective and repeatable, and thus generalizable. The exact same modelling approach is applied to all devices with only the device geometry changing. The results validate the overall approach with relatively low error, providing a general modelling methodology for accurate design of HF2V devices.
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.
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.
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.
The recent instances of seismic activity in Canterbury (2010/11) and Kaikōura (2016) in New Zealand have exposed an unexpected level of damage to non-structural components, such as buried pipelines and building envelope systems. The cost of broken buried infrastructure, such as pipeline systems, to the Christchurch Council was excessive, as was the cost of repairing building envelopes to building owners in both Christchurch and Wellington (due to the Kaikōura earthquake), which indicates there are problems with compliance pathways for both of these systems. Councils rely on product testing and robust engineering design practices to provide compliance certification on the suitability of product systems, while asset and building owners rely on the compliance as proof of an acceptable design. In addition, forensic engineers and lifeline analysts rely on the same product testing and design techniques to analyse earthquake-related failures or predict future outcomes pre-earthquake, respectively. The aim of this research was to record the actual field-observed damage from the Canterbury and Kaikōura earthquakes of seismic damage to buried pipeline and building envelope systems, develop suitable testing protocols to be able to test the systems’ seismic resilience, and produce prediction design tools that deliver results that reflect the collected field observations with better accuracy than the present tools used by forensic engineers and lifeline analysts. The main research chapters of this thesis comprise of four publications that describe the gathering of seismic damage to pipes (Publication 1 of 4) and building envelopes (Publication 2 of 4). Experimental testing and the development of prediction design tools for both systems are described in Publications 3 and 4. The field observation (discussed in Publication 1 of 4) revealed that segmented pipe joints, such as those used in thick-walled PVC pipes, were particularly unsatisfactory with respect to the joint’s seismic resilience capabilities. Once the joint was damaged, silt and other deleterious material were able to penetrate the pipeline, causing blockages and the shutdown of key infrastructure services. At present, the governing Standards for PVC pipes are AS/NZS 1477 (pressure systems) and AS/NZS 1260 (gravity systems), which do not include a protocol for evaluating the PVC pipes for joint seismic resilience. Testing methodologies were designed to test a PVC pipe joint under various different simultaneously applied axial and transverse loads (discussed in Publication 3 of 4). The goal of the laboratory experiment was to establish an easy to apply testing protocol that could fill the void in the mentioned standards and produce boundary data that could be used to develop a design tool that could predict the observed failures given site-specific conditions surrounding the pipe. A tremendous amount of building envelope glazing system damage was recorded in the CBDs of both Christchurch and Wellington, which included gasket dislodgement, cracked glazing, and dislodged glazing. The observational research (Publication 2 of 4) concluded that the glazing systems were a good indication of building envelope damage as the glazing had consistent breaking characteristics, like a ballistic fuse used in forensic blast analysis. The compliance testing protocol recognised in the New Zealand Building Code, Verification Method E2/VM1, relies on the testing method from the Standard AS/NZS 4284 and stipulates the inclusion of typical penetrations, such as glazing systems, to be included in the test specimen. Some of the building envelope systems that failed in the recent New Zealand earthquakes were assessed with glazing systems using either the AS/NZS 4284 or E2/VM1 methods and still failed unexpectedly, which suggests that improvements to the testing protocols are required. An experiment was designed to mimic the observed earthquake damage using bi-directional loading (discussed in Publication 4 of 4) and to identify improvements to the current testing protocol. In a similar way to pipes, the observational and test data was then used to develop a design prediction tool. For both pipes (Publication 3 of 4) and glazing systems (Publication 4 of 4), experimentation suggests that modifying the existing testing Standards would yield more realistic earthquake damage results. The research indicates that including a specific joint testing regime for pipes and positioning the glazing system in a specific location in the specimen would improve the relevant Standards with respect to seismic resilience of these systems. Improving seismic resilience in pipe joints and glazing systems would improve existing Council compliance pathways, which would potentially reduce the liability of damage claims against the government after an earthquake event. The developed design prediction tool, for both pipe and glazing systems, uses local data specific to the system being scrutinised, such as local geology, dimensional characteristics of the system, actual or predicted peak ground accelerations (both vertically and horizontally) and results of product-specific bi-directional testing. The design prediction tools would improve the accuracy of existing techniques used by forensic engineers examining the cause of failure after an earthquake and for lifeline analysts examining predictive earthquake damage scenarios.
