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.
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.
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.
In the last two decades, New Zealand (NZ) has experienced significant earthquakes, including the 2010 M 7.2 Darfield, 2011 M 6.2 Christchurch, and 2016 M 7.8 Kaikōura events. Amongst these large events, tens of thousands of smaller earthquakes have occurred. While previous event and ground-motion databases have analyzed these events, many events below M 4 have gone undetected. The goal of this study is to expand on previous databases, particularly for small magnitude (M<4) and low-amplitude ground motions. This new database enables a greater understanding of regional variations within NZ and contributes to the validity of internationally developed ground-motion models. The database includes event locations and magnitude estimates with uncertainty considerations, and tectonic type assessed in a hierarchical manner. Ground motions are extracted from the GeoNet FDSN server and assessed for quality using a neural network classification approach. A deep neural network approach is also utilized for picking P and S phases for determination of event hypocentres. Relative hypocentres are further improved by double-difference relocation and will contribute toward developing shallow (< 50 km) seismic tomography models. Analysis of the resulting database is compared with previous studies for discussion of implications toward national hazard prediction models.
The overarching goal of this dissertation is to improve predictive capabilities of geotechnical seismic site response analyses by incorporating additional salient physical phenomena that influence site effects. Specifically, multidimensional wave-propagation effects that are neglected in conventional 1D site response analyses are incorporated by: (1) combining results of 3D regional-scale simulations with 1D nonlinear wave-propagation site response analysis, and (2) modelling soil heterogeneity in 2D site response analyses using spatially-correlated random fields to perturb soil properties. A method to combine results from 3D hybrid physics-based ground motion simulations with site-specific nonlinear site response analyses was developed. The 3D simulations capture 3D ground motion phenomena on a regional scale, while the 1D nonlinear site response, which is informed by detailed site-specific soil characterization data, can capture site effects more rigorously. Simulations of 11 moderate-to-large earthquakes from the 2010-2011 Canterbury Earthquake Sequence (CES) at 20 strong motion stations (SMS) were used to validate simulations with observed ground motions. The predictions were compared to those from an empirically-based ground motion model (GMM), and from 3D simulations with simplified VS30- based site effects modelling. By comparing all predictions to observations at seismic recording stations, it was found that the 3D physics-based simulations can predict ground motions with comparable bias and uncertainty as the GMM, albeit, with significantly lower bias at long periods. Additionally, the explicit modelling of nonlinear site-response improves predictions significantly compared to the simplified VS30-based approach for soft-soil or atypical sites that exhibit exceptionally strong site effects. A method to account for the spatial variability of soils and wave scattering in 2D site response analyses was developed and validated against a database of vertical array sites in California. The inputs required to run the 2D analyses are nominally the same as those required for 1D analyses (except for spatial correlation parameters), enabling easier adoption in practice. The first step was to create the platform and workflow, and to perform a sensitivity study involving 5,400 2D model realizations to investigate the influence of random field input parameters on wave scattering and site response. Boundary conditions were carefully assessed to understand their effect on the modelled response and select appropriate assumptions for use on a 2D model with lateral heterogeneities. Multiple ground-motion intensity measures (IMs) were analyzed to quantify the influence from random field input parameters and boundary conditions. It was found that this method is capable of scattering seismic waves and creating spatially-varying ground motions at the ground surface. The redistribution of ground-motion energy across wider frequency bands, and the scattering attenuation of high-frequency waves in 2D analyses, resemble features observed in empirical transfer functions (ETFs) computed in other studies. The developed 2D method was subsequently extended to more complicated multi-layer soil profiles and applied to a database of 21 vertical array sites in California to test its appropriate- ness for future predictions. Again, different boundary condition and input motion assumptions were explored to extend the method to the in-situ conditions of a vertical array (with a sensor embedded in the soil). ETFs were compared to theoretical transfer functions (TTFs) from conventional 1D analyses and 2D analyses with heterogeneity. Residuals of transfer-function- based IMs, and IMs of surface ground motions, were also used as validation metrics. The spatial variability of transfer-function-based IMs was estimated from 2D models and compared to the event-to-event variability from ETFs. This method was found capable of significantly improving predictions of median ETF amplification factors, especially for sites that display higher event-to-event variability. For sites that are well represented by 1D methods, the 2D approach can underpredict amplification factors at higher modes, suggesting that the level of heterogeneity may be over-represented by the 2D random field models used in this study.
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.
A buckling-restrained braced frame (BRBF) is a structural bracing system that provides lateral strength and stiffness to buildings and bridges. They were first developed in Japan in the 1970s (Watanabe et al. 1973, Kimura et al. 1976) and gained rapid acceptance in the United States after the Northridge earthquake in 1994 (Bruneau et al. 2011). However, it was not until the Canterbury earthquakes of 2010/2011, that the New Zealand construction market saw a significant uptake in the use of buckling-restrained braces (BRBs) in commercial buildings (MacRae et al. 2015). In New Zealand there is not yet any documented guidance or specific instructions in regulatory standards for the design of BRBFs. This makes it difficult for engineers to anticipate all the possible stability and strength issues within a BRBF system and actively mitigate them in each design. To help ensure BRBF designs perform as intended, a peer review with physical testing are needed to gain building compliance in New Zealand. Physical testing should check the manufacturing and design of each BRB (prequalification testing), and the global strength and stability of each BRB its frame (subassemblage testing). However, the financial pressures inherent in commercial projects has led to prequalification testing (BRB only testing) being favoured without adequate design specific subassemblage testing. This means peer reviewers have to rely on BRB suppliers for assurances. This low regulation environment allows for a variety of BRBF designs to be constructed without being tested or well understood. The concern is that there may be designs that pose risk and that issues are being overlooked in design and review. To improve the safety and design of BRBFs in New Zealand, this dissertation studies the behaviour of BRBs and how they interact with other frame components. Presented is the experimental test process and results of five commercially available BRB designs (Chapter 2). It discusses the manufacturing process, testing conditions and limitations of observable information. It also emphasises that even though subassemblage testing is impractical, uniaxial testing of the BRB only is not enough, as this does not check global strength or stability. As an alternative to physical testing, this research uses computer simulation to model BRB behaviour. To overcome the traditional challenges of detailed BRB modelling, a strategy to simulate the performance of generic BRB designs was developed (Chapter 3). The development of nonlinear material and contact models are important aspects of this strategy. The Chaboche method is employed using a minimum of six backstress curves to characterize the combined isotropic and kinematic hardening exhibited by the steel core. A simplified approach, adequate for modelling the contact interaction between the restrainer and the core was found. Models also capture important frictional dissipation as well as lateral motion and bending associated with high order constrained buckling of the core. The experimental data from Chapter 2 was used to validate this strategy. As BRBs resist high compressive loading, global stability of the BRB and gusseted connection zone need to be considered. A separate study was conducted that investigated the yielding and buckling strength of gusset plates (Chapter 4). The stress distribution through a gusset plate is complex and difficult to predict because the cross-sectional area of gusset plate is not uniform, and each gusset plate design is unique in shape and size. This has motivated design methods that approximate yielding of gusset plates. Finite element modelling was used to study the development of yielding, buckling and plastic collapse behaviour of a brace end bolted to a series of corner gusset plates. In total 184 variations of gusset plate geometries were modelled in Abaqus®. The FEA modelling applied monotonic uniaxial load with an imperfection. Upon comparing results to current gusset plate design methods, it was found that the Whitmore width method for calculating the yield load of a gusset is generally un-conservative. To improve accuracy and safety in the design of gusset plates, modifications to current design methods for calculating the yield area and compressive strength for gusset plates is proposed. Bolted connections are a popular and common connection type used in BRBF design. Global out-of-plane stability tends to govern the design for this connection type with numerous studies highlighting the risk of instability initiated by inelasticity in the gussets, neck of the BRB end and/or restrainer ends. Subassemblage testing is the traditional method for evaluating global stability. However, physical testing of every BRBF variation is cost prohibitive. As such, Japan has developed an analytical approach to evaluate out-of-plane stability of BRBFs and incorporated this in their design codes. This analytical approach evaluates the different BRB components under possible collapse mechanisms by focusing on moment transfer between the restrainer and end of the BRB. The approach have led to strict criteria for BRBF design in Japan. Structural building design codes in New Zealand, Europe and the United States do not yet provide analytical methods to assess BRB and connection stability, with prototype/subassemblage testing still required as the primary means of accreditation. Therefore it is of interest to investigate the capability of this method to evaluate stability of BRBs designs and gusset plate designs used in New Zealand (including unstiffened gusset connection zones). Chapter 5 demonstrates the capability of FEA to study to the performance of a subassemblage test under cyclic loading – resembling that of a diagonal ground storey BRBF with bolted connections. A series of detailed models were developed using the strategy presented in Chapter 3. The geometric features of BRB 6.5a (Chapter 2) were used as a basis for the BRBs modelled. To capture the different failure mechanisms identified in Takeuchi et al. (2017), models varied the length that the cruciform (non-yielding) section inserts into the restrainer. Results indicate that gusset plates designed according to New Zealand’s Steel Structures Standard (NZS 3404) limit BRBF performance. Increasing the thickness of the gusset plates according to modifications discussed in Chapter 4, improved the overall performance for all variants (except when Lin/ Bcruc = 0.5). The effect of bi-directional loading was not found to notably affect out-of-plane stability. Results were compared against predictions made by the analytical method used in Japan (Takeuchi method). This method was found to be generally conservative is predicting out-of-plane stability of each BRBF model. Recommendations to improve the accuracy of Takeuchi’s method are also provided. The outcomes from this thesis should be helpful for BRB manufacturers, researchers, and in the development of further design guidance of BRBFs.
