Search

found 110 results

Images for motion; more images...

Research papers, University of Canterbury Library

This paper provides a summary of the ground motions observed in the recent Canterbury, New Zealand earthquake sequence. The sequence occurred in a region of relatively moderate seismicity, 130km to the east of the Alpine Fault, the major plate-boundary in the region. From an engineering perspective, the sequence has been primarily comprised of the initial 04/09/2010 Darfield earthquake (Mw7.1) followed by the 22/02/2011 Christchurch earthquake (Mw6.3), and two aftershocks on 13/06/ 2011 (Mw5.3 and 6.0, respectively). The dense spacing of strong motions in the region, and their close proximity to the respective causative faults, has resulted in strong ground motions far exceeding the previous catalogue of strong motion observed in New Zealand. The observed ground motions have exhibited clear evidence of: (i) near-source directivity; (ii) sedimentary basin focusing, amplification and basin effect refraction; (iii) non-linear site response; (iv) cyclic mobility postliquefaction; and (v) extreme vertical ground motions exceeding 2g, among others.

Research papers, University of Canterbury Library

This paper concerns the explicit consideration of near-fault directivity in conventional ground motion prediction models, and its implication for probabilistic seismic hazard analysis (PSHA) in New Zealand. The proposed approach utilises recently developed models by Shahi & Baker (2011), which account for both the 'narrowband' nature of the directivity pulse on spectral ordinates, and the probability of pulse occurrence at the site of interest. Furthermore, in order to correctly consider directivity, distributed seismicity sources are considered as finite-faults, as opposed to their (incorrect) conventional treatment as point-sources. The significance of directivity on hazard analysis results is illustrated for various vibration periods at generic sites located in Christchurch and Otira, two locations whose seismic hazard is comprised of notably different seismic sources. When compared to the PSHA results considering directivity and distributed seismicity as finite faults, it is shown that the NZS1170.5:2004 directivity factor is notably unconservative for all vibration periods in Otira (i.e. high seismic hazard region); and unconservative for Christchurch at short-to-moderate vibration periods ( < 3s); but conservative at long periods ( > 4s).

Research papers, University of Canterbury Library

This presentation summarizes the development of high-resolution surficial soil velocity models in the Canterbury, New Zealand basin. Shallow (<30m) shear wave velocities were primarily computed based on a combination of a large database of over 15,000 cone penetration test (CPT) logs in and around Christchurch, and a recently-developed Christchurch-specific empirical correlation between soil shear wave velocity and CPT. Large active-source testing at 22 locations and ambient-wavefield surface wave and H/V testing at over 80 locations were utilized in combination with 1700 water well logs to constrain the inter-bedded stratigraphy and velocity of Quaternary sediments up to depths of several hundred meters. Finally, seismic reflection profiles and the ambient-wavefield surface wave data provide constraint on velocities from several hundred meters to several kilometres. At all depths, the high resolution data illustrates the complexity of the soil conditions in the region, and the developed 3D models are presently being used in broadband ground motion simulations to further interpret the observed strong ground motions in the 2010-2011 Canterbury earthquake sequence.

Research papers, University of Canterbury Library

The Canterbury Earthquake Sequence (CES) of 2010-2011 caused widespread liquefaction in many parts of Christchurch. Observations from the CES highlight some sites were liquefaction was predicted by the simplified method but did not manifest. There are a number of reasons why the simplified method may over-predict liquefaction, one of these is the dynamic interaction between soil layers within a stratified deposit. Soil layer interaction occurs through two key mechanisms; modification of the ground motion due to seismic waves passing through deep liquefied layers, and the effect of pore water seepage from an area of high excess pore water pressure to the surrounding soil. In this way, soil layer interaction can significantly alter the liquefaction behaviour and surface manifestation of soils subject to seismic loading. This research aimed to develop an understanding of how soil layer interaction, in particular ground motion modification, affects the development of excess pore water pressures and liquefaction manifestation in a soil deposit subject to seismic loading. A 1-D soil column time history Effective Stress Analysis (ESA) was conducted to give an in depth assessment of the development of pore pressures in a number of soil deposits. For this analysis, ground motions, soil profiles and model parameters were required for the ESA. Deconvolution of ground motions recorded at the surface during the CES was used to develop some acceleration time histories to input at the base of the soil-column model. An analysis of 55 sites around Christchurch, where detailed site investigations have been carried out, was then conducted to identify some simplified soil profiles and soil characteristics. From this analysis, four soil profiles representative of different levels of liquefaction manifestation were developed. These were; two thick uniform and vertically continuous sandy deposits that were representative of sites were liquefaction manifested in both the Mw 7.1 September 2010 and the Mw 6.3 February 2011 earthquakes, and two vertically discontinuous profiles with interlayered liquefiable and non-liquefiable layers representative of sites that did not manifest liquefaction in either the September 2010 or the February 2011 events. Model parameters were then developed for these four representative soil profiles through calibration of the constitutive model in element test simulations. Simulations were run for each of the four profiles subject to three levels of loading intensity. The results were analysed for the effect of soil layer interaction. These were then compared to a simplified triggering analysis for the same four profiles to determine where the simplified method was accurate in predicting soil liquefaction (for the continuous sandy deposits) and were it was less accurate (the vertically discontinuous deposits where soil layer interaction was a factor).

Images, UC QuakeStudies

A motion-blurred photograph of houses, with the Port Hills in the background. The photographer comments, "This I hope gives you a feel of what it feels like in an earthquake. When you spend your whole life thinking that you and your home are built on solid ground, it can be quite a shock when you find it is not. You can feel the house shaking like a dog with a toy, rising up violently underneath you or the most gentle form which is when the ground moves gently like a wave moving under a rowing boat. It is not just the movement, you often get a rumbling sound which can precede a violent shake or can result in no movement at all. This means that some vehicles can sound like the rumbling initially and in the early days would get your heart racing. Another form of stress is when big excavators as heavy as a tank move as you can feel the ground shake from streets away, but you do not always hear the engine. For most of us the problem when the shaking starts, is wondering if this is the start of an extremely violent earthquake or will it peter out".

Research papers, University of Canterbury Library

This poster provides a summary of the development of a 3D shallow (z<40m) shear wave velocity (Vs) model for the urban Christchurch, New Zealand region. The model is based on a recently developed Christchurch-specific empirical correlation between Vs and cone penetration test (CPT) data (McGann et al. 2014a,b) and the large high-density database of CPT logs in the greater Christchurch urban area (> 15,000 logs as of 01/01/2014). In particular, the 3D model provides shear wave velocities for the surficial Springston Formation, Christchurch Formation, and Riccarton gravel layers which generally comprise the upper 40m in the Christchurch urban area. Point-estimates are provided on a 200m-by- 200m grid from which interpolation to other locations can be performed. This model has applications for future site characterization and numerical modeling efforts via maps of timeaveraged Vs over specific depths (e.g. Vs30, Vs10) and via the identification of typical Vs profiles for different regions and soil behaviour types within Christchurch. In addition, the Vs model can be used to constrain the near-surface velocities for the 3D seismic velocity model of the Canterbury basin (Lee et al. 2014) currently being developed for the purpose of broadband ground motion simulation.

Research papers, University of Canterbury Library

A major hazard accompanying earthquake shaking in areas of steep topography is the detachment of rocks from bedrock outcrops that subsequently slide, roll, or bounce downslope (i.e. rockfalls). The 2010-2011 Canterbury earthquake sequence caused recurrent and severe rockfall in parts of southern Christchurch. Coseismic rockfall caused five fatalities and significant infrastructural damage during the 2011 Mw 6.2 Christchurch earthquake. Here we examine a rockfall site in southern Christchurch in detail using geomorphic mapping, lidar analysis, geochronology (cosmogenic 3He dating, radiocarbon dating, optically stimulated luminescence (OSL) from quartz, infrared stimulated luminescence from K-feldspar), numerical modeling of rockfall boulder trajectories, and ground motion prediction equations (GMPEs). Rocks fell from the source cliff only in earthquakes with interpolated peak ground velocities exceeding ~10 cm/s; hundreds of smaller earthquakes did not produce rockfall. On the basis of empirical observations, GMPEs and age chronologies we attribute paleo-rockfalls to strong shaking in prehistoric earthquakes. We conclude that earthquake shaking of comparable intensity to the strongest contemporary earthquakes in Christchurch last occurred at this site approximately 5000 to 7000 years ago, and that in some settings, rockfall deposits provide useful proxies for past strong ground motions.

Research papers, University of Canterbury Library

Timber-based hybrid structures provide a prospective solution for utilizing environmentally friendly timber material in the construction of mid-rise or high-rise structures. This study mainly focuses on structural damage evaluation for a type of timber-steel hybrid structures, which incorporate prefabricated light wood frame shear walls into steel moment-resisting frames (SMRFs). The structural damage of such a hybrid structure was evaluated through shake table tests on a four-story large-scale timber-steel hybrid structure. Four ground motion records (i.e., Wenchuan earthquake, Canterbury earthquake, El-Centro earthquake, and Kobe earthquake) were chosen for the tests, with the consideration of three different probability levels (i.e., minor, moderate and major earthquakes) for each record. During the shake table tests, the hybrid structure performed quite well with visual damage only to wood shear walls. No visual damage in SMRF and the frame-to-wall connections was observed. The correlation of visual damage to seismic intensity, modal-based damage index and inter-story drift was discussed. The reported work provided a basis of knowledge for performance-based seismic design (PBSD) for such timber-based hybrid structures.

Research papers, University of Canterbury Library

In order to provide information related to seismic vulnerability of non-ductile reinforced concrete (RC) frame buildings, and as a complementary investigation on innovative feasible retrofit solutions developed in the past six years at the University of Canterbury on pre-19170 reinforced concrete buildings, a frame building representative of older construction practice was tested on the shake table. The specimen, 1/2.5 scale, consists of two 3-storey 2-bay asymmetric frames in parallel, one interior and one exterior, jointed together by transverse beams and floor slabs. The as-built (benchmark) specimen was first tested under increasing ground motion amplitudes using records from Loma Prieta Earthquake (California, 1989) and suffered significant damage at the upper floor, most of it due to lap splices failure. As a consequence, in a second stage, the specimen was repaired and modified by removing the concrete in the lap splice region, welding the column longitudinal bars, replacing the removed concrete with structural mortar, and injecting cracks with epoxy resin. The modified as-built specimen was then tested using data recorded during Darfield (New Zealand, 2010) and Maule (Chile, 2010) Earthquakes, with whom the specimen showed remarkably different responses attributed to the main variation in frequency content and duration. In this contribution, the seismic performance of the three series of experiments are presented and compared.

Research papers, University of Canterbury Library

This thesis documents the development and demonstration of an assessment method for analysing earthquake-related damage to concrete waste water gravity pipes in Christchurch, New Zealand, following the 2010-2011 Canterbury Earthquake Sequence (CES). The method is intended to be internationally adaptable to assist territorial local authorities with improving lifelines infrastructure disaster impact assessment and improvements in resilience. This is achieved through the provision of high-resolution, localised damage data, which demonstrate earthquake impacts along the pipe length. The insights gained will assist decision making and the prioritisation of resources following earthquake events to quickly and efficiently restore network function and reduce community impacts. The method involved obtaining a selection of 55 reinforced concrete gravity waste water pipes with available Closed-Circuit Television (CCTV) inspection footage filmed before and after the CES. The pipes were assessed by reviewing the recordings, and damage was mapped to the nearest metre along the pipe length using Geographic Information Systems. An established, systematic coding process was used for reporting the nature and severity of the observed damage, and to differentiate between pre-existing and new damage resulting from the CES. The damage items were overlaid with geospatial data such as Light Detection and Ranging (LiDAR)-derived ground deformation data, Liquefaction Resistance Index data and seismic ground motion data (Peak Ground acceleration and Peak Ground Velocity) to identify potential relationships between these parameters and pipe performance. Initial assessment outcomes for the pipe selection revealed that main pipe joints and lateral connections were more vulnerable than the pipe body during a seismic event. Smaller diameter pipes may also be more vulnerable than larger pipes during a seismic event. Obvious differential ground movement resulted in increased local damage observations in many cases, however this was not obvious for all pipes. Pipes with older installation ages exhibited more overall damage prior to a seismic event, which is likely attributable to increased chemical and biological deterioration. However, no evidence was found relating pipe age to performance during a seismic event. No evidence was found linking levels of pre-CES damage in a pipe with subsequent seismic performance, and seismic performance with liquefaction resistance or magnitude of seismic ground motion. The results reported are of limited application due to the small demonstration sample size, but reveal the additional level of detail and insight possible using the method presented in this thesis over existing assessment methods, especially in relation to high resolution variations along the length of the pipe such as localised ground deformations evidenced by LiDAR. The results may be improved by studying a larger and more diverse sample pool, automating data collection and input processes in order to improve efficiency and consider additional input such as pipe dip and cumulative damage over a large distance. The method is dependent on comprehensive and accurate pre-event CCTV assessments and LIDAR data so that post-event data could be compared. It is proposed that local territorial authorities should prioritise acquiring this information as a first important step towards improving the seismic resilience of a gravity waste water pipe network.

Research papers, University of Canterbury Library

This thesis documents the development and demonstration of an assessment method for analysing earthquake-related damage to concrete waste water gravity pipes in Christchurch, New Zealand, following the 2010-2011 Canterbury Earthquake Sequence (CES). The method is intended to be internationally adaptable to assist territorial local authorities with improving lifelines infrastructure disaster impact assessment and improvements in resilience. This is achieved through the provision of high-resolution, localised damage data, which demonstrate earthquake impacts along the pipe length. The insights gained will assist decision making and the prioritisation of resources following earthquake events to quickly and efficiently restore network function and reduce community impacts. The method involved obtaining a selection of 55 reinforced concrete gravity waste water pipes with available Closed-Circuit Television (CCTV) inspection footage filmed before and after the CES. The pipes were assessed by reviewing the recordings, and damage was mapped to the nearest metre along the pipe length using Geographic Information Systems. An established, systematic coding process was used for reporting the nature and severity of the observed damage, and to differentiate between pre-existing and new damage resulting from the CES. The damage items were overlaid with geospatial data such as Light Detection and Ranging (LiDAR)-derived ground deformation data, Liquefaction Resistance Index data and seismic ground motion data (Peak Ground acceleration and Peak Ground Velocity) to identify potential relationships between these parameters and pipe performance. Initial assessment outcomes for the pipe selection revealed that main pipe joints and lateral connections were more vulnerable than the pipe body during a seismic event. Smaller diameter pipes may also be more vulnerable than larger pipes during a seismic event. Obvious differential ground movement resulted in increased local damage observations in many cases, however this was not obvious for all pipes. Pipes with older installation ages exhibited more overall damage prior to a seismic event, which is likely attributable to increased chemical and biological deterioration. However, no evidence was found relating pipe age to performance during a seismic event. No evidence was found linking levels of pre-CES damage in a pipe with subsequent seismic performance, and seismic performance with liquefaction resistance or magnitude of seismic ground motion. The results reported are of limited application due to the small demonstration sample size, but reveal the additional level of detail and insight possible using the method presented in this thesis over existing assessment methods, especially in relation to high resolution variations along the length of the pipe such as localised ground deformations evidenced by LiDAR. The results may be improved by studying a larger and more diverse sample pool, automating data collection and input processes in order to improve efficiency and consider additional input such as pipe dip and cumulative damage over a large distance. The method is dependent on comprehensive and accurate pre-event CCTV assessments and LIDAR data so that post-event data could be compared. It is proposed that local territorial authorities should prioritise acquiring this information as a first important step towards improving the seismic resilience of a gravity waste water pipe network.

Research papers, University of Canterbury Library

The Mw 7.8 Kaikōura earthquake ruptured ~200 km at the ground surface across the New Zealand plate boundary zone in the northern South Island. This study was conducted in an area of ~600 km2 in the epicentral region where the faults comprise two main non-coplanar sets that strike E-NE and NNE-NW with mainly steep dips (60о-80°). Analysis of the surface rupture using field and LiDAR data provides new information on the dimensions, geometries and kinematics of these faults which was not previously available from pre-earthquake active faults or bedrock structure. The more northerly striking fault set are sub-parallel to basement bedding and accommodated predominantly left-lateral reverse slip with net slips of ~1 and ~5 m for the Stone Jug and Leader faults, respectively. The E-NE striking Conway-Charwell and The Humps faults accrued right-lateral to oblique reverse with net slips of ~2 and ~3 m, respectively. The faults form a hard-linked system dominated by kinematics consistent with the ~260° trend of the relative plate motion vector and the transpressional structures recorded across the plate boundary in the NE South Island. Interaction and intersection of the main fault sets facilitated propagation of the earthquake and transfer of slip northwards across the plate boundary zone.

Research papers, University of Canterbury Library

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.

Research papers, University of Canterbury Library

This paper presents a seismic velocity model of Canterbury, New Zealand based on 3D geologic surfaces and velocities from a range of data sources. The model provides the 3D crustal structure in the region at multiple length scales for seismic wave propagation simulations, such as broadband ground motion and shallow site response analyses related to understanding the ground motions and site responses during the 2010- 2011 Canterbury earthquakes. Pre-Quaternary geologic horizons are calculated based on the reinterpretation of a comprehensive network of seismic reflection surveys from seven different campaigns over the past 50 years, as well as point constraints across an array of petroleum industry drill holes. Particular attention is given to a detailed representation of Quaternary stratigraphy, representing shallow (z<250m) near-surface layers in the model. Seismic velocities are obtained from seismic reflection processing (for Vp) and also recently performed active and passive surface wave analyses (for Vs). Over 1,700 water wells in the region are used to constrain the complex inter-bedded Quaternary stratigraphy (gravels, sands, silts, organics etc.) near the coastline, including beneath urban Christchurch, which has resulted from fluvial deposition and marine regression and transgression. For the near-surface Springston and Christchurch Formations in the Christchurch urban area (z<50m), high-spatial resolution seismic velocities (including Vs30 ) were obtained from over 13,000 cone penetration tests combined with a recently developed CPT-Vs correlation.

Research papers, University of Canterbury Library

The 4th of September 2010 Mw 7.1 Darfield (Canterbury) earthquake had generated significant ground shaking within the Christchurch Central Business District (CBD). Despite the apparently significant shaking, the observed structural damage for pre-1970s reinforced concrete (RC) buildings was indeed limited and lower than what was expected for such typology of buildings. This paper explores analytically and qualitatively the different aspects of the "apparent‟ good seismic performance of the pre-1970s RC buildings in the Christchurch CBD, following the earthquake reconnaissance survey by the authors. Damage and building parameters survey result, based on a previously established inventory of building stock of these non-ductile RC buildings, is briefly reported. From an inventory of 75 buildings, one building was selected as a numerical case-study to correlate the observed damage with the non-linear analyses. The result shows that the pre-1970s RC frame buildings performed as expected given the intensity of the ground motion shaking during the Canterbury earthquake. Given the brittle nature of this type of structure, it was demonstrated that more significant structural damage and higher probability of collapse could occur when the buildings were subjected to alternative input signals with different frequency content and duration characteristics and still compatible to the seismicity hazard for Christchurch CBD.

Research papers, University of Canterbury Library

Predicting building collapse due to seismic motion is critical in design and more so after a major event. Damaged structures can appear sound, but collapse under following major events. There can thus be significant risk in decision making after a major seismic event concerning the safe occupation of a building or surrounding areas, versus the unknown impact of unknown major aftershocks. Model-based pushover analyses are effective if the structural properties are well understood, which is not valid post-event when this risk information is most useful. This research combines Hysteresis Loop Analysis (HLA) structural health monitoring (SHM) and Incremental Dynamic Analysis (IDA) methods to determine collapse capacity and probability of collapse for a specific structure, at any time, a range of earthquake excitations to ensure robustness. The nonlinear dynamic analysis method presented enables constant updating of building performance predictions using post-event SHM results. The resulting combined methods provide near real-time updating of collapse fragility curves as events progress, quantifying the change of collapse probability or seismic induced losses for decision-making - a novel, higher resolution risk analysis than previously available. The methods are not computationally expensive and there is no requirement for a validated numerical model. Results show significant potential benefits and a clear evolution of risk. They also show clear need for extending SHM toward creating improved predictive models for analysis of subsequent events, where the Christchurch series of 2010-2011 had significant post-event aftershocks after each main event. Finally, the overall method is generalisable to any typical engineering demand parameter.

Research papers, University of Canterbury Library

Earthquakes cause significant damage to buildings due to strong vibration of the ground. Levitating houses using magnets and electromagnets would provide a complete isolation of ground motion for protecting buildings from seismic damage. Two types of initial configuration for the electromagnet system were proposed with the same air gap (10mm) between the electromagnet and reluctance plate. Both active and passive controller are modelled to investigate the feasibility of using a vibration control system for stabilizing the magnetic system within the designed air gap (10mm) in the vertical direction. A nonlinear model for the magnetic system is derived to implement numerical simulation of structural response under the earthquake record in Christchurch Botanic Gardens on 21 February 2011. The performance of the uncontrolled and the controlled systems are compared and the optimal combination of control gains are determined for the PID active controller. Simulation results show both active PID controller with constant and nonlinear attracting force are able to provide an effective displacement control within the required air gap (+/-5mm). The maximum control force demand for the PID controller in the presence of nonlinear attracting force is 4.1kN, while the attracting force in equilibrium position is 10kN provided by the electromagnet. These results show the feasibility of levitating a house using the current electromagnet and PID controller. Finally, initial results of passive control using two permanent magnets or dampers show the structural responses can be effectively reduced and centralized to +/-1mm using a nonlinear centring barrier function.

Research papers, University of Canterbury Library

The objective of this study is to examine the influence of near-fault motions on liquefaction triggering in Christchurch and neighboring towns during the 2010-2011 Canterbury earthquake sequence (CES). The CES began with the 4 September 2010, Mw7.1 Darfield earthquake and included up to ten events that triggered liquefaction. However, most notably, widespread liquefaction was induced by the Darfield earthquake and the Mw6.2, 22 February 2011 Christchurch earthquake. Of particular relevance to this study is the forward directivity effects that were prevalent in the motions recorded during the Darfield earthquake, and to a much lesser extent, during the Christchurch earthquake. A 2D variant of the Richart-Newmark fatigue theory was used to compute the equivalent number of cycles (neq) for the ground motions, where volumetric strain was used as the damage metric. This study is unique because it considers the contribution and phasing of both the fault-normal and fault-parallel components of motion on neq and the magnitude scaling factor (MSF). It was found that when the fault-normal and fault-parallel motions were treated individually, the former yielded a lower neq than the latter. Additionally, when the combined effects of fault-normal and fault-parallel components were considered, it was found that the MSF were higher than those commonly used. This implies that motions containing near-fault effects are less demanding on the soil than motions that do not. This may be one of several factors that resulted in less severe liquefaction occurring during the Darfield earthquake than the Christchurch earthquake.

Research papers, University of Canterbury Library

Live monitoring data and simple dynamic reduced-order models of the Christchurch Women’s Hospital (CWH) help explain the performance of the base isolation (BI) system of the hospital during the series of Canterbury earthquakes in 2011-2012. A Park-Wen-Ang hysteresis model is employed to simulate the performance of the BI system and results are compared to measured data recorded above the isolation layer and on the 6th story. Simplified single, two and three degree of freedom models (SDOF, 2DOF and 3DOF) show that the CWH structure did not behave as an isolated but as a fixed-base structure. Comparisons of accelerations and deflections between simulated and monitored data show a good match for isolation stiffness values of approximately two times of the value documented in the design specification and test protocol. Furthermore, an analysis of purely measured data revealed very little to no relative motion across the isolators for large events of moment magnitude scale (Mw) 5.8 and 6.0 that occurred within 3 hours of each other on December 23, 2011. One of the major findings is that the BI system during the seismic events on December 23, 2011 did not yield and that the superstructure performed as a fixed-base building, indicating a need to reevaluate the analysis, design and implementation of these structures.

Research papers, University of Canterbury Library

Rapid, accurate structural health monitoring (SHM) assesses damage to optimise decision-making. Many SHM methods are designed to track nonlinear stiffness changes as damage. However, highly nonlinear pinched hysteretic systems are problematic in SHM. Model-based SHM often fails as any mismatch between model and measured response dynamics leads to significant error. Thus, modelfree methods of hysteresis loop tracking methods have emerged. This study compares the robustness and accuracy in the presence of significant measurement noise of the proven hysteresis loop analysis (HLA) SHM method with 3 emerging model-free methods and 2 further novel adaptations of these methods using a highly nonlinear, 6-story numerical structure to provide a known ground-truth. Mean absolute errors in identifying a known nonlinear stiffness trajectory assessed at four points over two successive ground motion inputs from September 2010 and February 2011 in Christchurch range from 1.71-10.52%. However, the variability is far wider with maximum errors ranging from 3.90-49.72%, where the second largest maximum absolute error was still 19.74%. The lowest mean and maximum absolute errors were for the HLA method. The next best method had mean absolute error of 2.92% and a maximum of 10.51%. These results show the clear superiority of the HLA method over all current emerging model-free methods designed to manage the highly nonlinear pinching responses common in reinforced concrete structures. These results, combined with high robustness and accuracy in scaled and fullscale experimental studies, provide further validation for using HLA for practical implementation.

Research papers, University of Canterbury Library

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.

Research papers, University of Canterbury Library

The Canterbury earthquake and aftershock sequence in New Zealand during 2010-2011 subjected the city’s structures to a significant accumulated cyclic demand and raised significant questions regarding the low-cycle fatigue demands imposed upon the structures. There is a significant challenge to quantify the level of cumulative demand imposed on structures and to assess the percentage of a structure's fatigue life that has been consumed as a result of this earthquake sequence. It is important to be able to quantify the cumulative demand to determine how a building will perform in a subsequent large earthquake and inform repair and re-occupancy decisions. This paper investigates the cumulative fatigue demand for a structure located within the Christchurch Central Business District (CBD). Time history analysis and equivalent cycle counting methods are applied across the Canterbury earthquake sequence, using key events from September 4th 2010 and February 22nd , 2011 main shocks. The estimate of the cumulative fatigue demand is then compared to the expected capacity of a case study reinforced concrete bridge pier, to undertake a structure-specific fatigue assessment. The analysis is undertaken to approximate the portion of the structural fatigue capacity that has been consumed, and how much residual capacity remains. Results are assessed for recordings at the four Christchurch central city strong motion recording sites installed by the GeoNet programme, to provide an estimate of variation in results. The computed cyclic demand results are compared to code-based design methods and as assessment of the inelastic displacement demand of the reinforcing steel. Results are also presented in a fragility context where a de minimis (inconsequential), irreparable damage and full fatigue fracture are defined to provide a probabilistic assessment of the fatigue damage incurred. This methodology can provide input into the overall assessment of fatigue demands and residual capacity.

Research papers, University of Canterbury Library

© 2019, Springer-Verlag GmbH Germany, part of Springer Nature. Prediction of building collapse due to significant seismic motion is a principle objective of earthquake engineers, particularly after a major seismic event when the structure is damaged and decisions may need to be made rapidly concerning the safe occupation of a building or surrounding areas. Traditional model-based pushover analyses are effective, but only if the structural properties are well understood, which is not the case after an event when that information is most useful. This paper combines hysteresis loop analysis (HLA) structural health monitoring (SHM) and incremental dynamic analysis (IDA) methods to identify and then analyse collapse capacity and the probability of collapse for a specific structure, at any time, a range of earthquake excitations to ensure robustness. This nonlinear dynamic analysis enables constant updating of building performance predictions following a given and subsequent earthquake events, which can result in difficult to identify deterioration of structural components and their resulting capacity, all of which is far more difficult using static pushover analysis. The combined methods and analysis provide near real-time updating of the collapse fragility curves as events progress, thus quantifying the change of collapse probability or seismic induced losses very soon after an earthquake for decision-making. Thus, this combination of methods enables a novel, higher-resolution analysis of risk that was not previously available. The methods are not computationally expensive and there is no requirement for a validated numerical model, thus providing a relatively simpler means of assessing collapse probability immediately post-event when such speed can provide better information for critical decision-making. Finally, the results also show a clear need to extend the area of SHM toward creating improved predictive models for analysis of subsequent events, where the Christchurch series of 2010–2011 had significant post-event aftershocks.

Research papers, University of Canterbury Library

During the past two decades, the focus has been on the need to provide communities with structures that undergo minimal damage after an earthquake event while still being cost competitive. This has led to the development of high performance seismic resisting systems, and advances in design methodologies, in order respect this demand efficiently. This paper presents the experimental response of four pre-cast, post-tensioned rocking wall systems tested on the shake-table at the University of Canterbury. The wall systems were designed as a retrofit solution for an existing frame building, but are equally applicable for use in new design. Design of the wall followed a performance-based retrofit strategy in which structural limit states appropriate to both the post-tensioned wall and the existing building were considered. Dissipation for each of the four post-tensioned walls was provided via externally mounted devices, located in parallel to post-tensioned tendons for re-centring. This allowed the dissipation devices to be easily replaced or inspected following a major earthquake. Each wall was installed with viscous fluid dampers, tension-compression yielding steel dampers, a combination of both or no devices at all – thus relying on contact damping alone. The effectiveness of both velocity and displacement dependant dissipation are investigated for protection against far-field and velocity-pulse ground motion characteristics. The experimental results validate the behaviour of ‘Advanced Flag-Shape’ rocking, dissipating solutions which have been recently proposed and numerically tested. Maximum displacements and material strains were well controlled and within acceptable bounds, and residual deformations were minimal due to the re-centring contribution from the post-tensioned tendons. Damage was confined to inelastic yielding (or fluid damping) of the external dampers.

Research papers, University of Canterbury Library

The purpose of this thesis is to evaluate the seismic response of the UC Physics Building based on recorded ground motions during the Canterbury earthquakes, and to use the recorded response to evaluate the efficacy of various conventional structural analysis modelling assumptions. The recorded instrument data is examined and analysed to determine how the UC Physics Building performed during the earthquake-induced ground motions. Ten of the largest earthquake events from the 2010-11 Canterbury earthquake sequence are selected in order to understand the seismic response under various levels of demand. Peak response amplitude values are found which characterise the demand from each event. Spectral analysis techniques are utilised to find the natural periods of the structure in each orthogonal direction. Significant torsional and rocking responses are also identified from the recorded ground motions. In addition, the observed building response is used to scrutinise the adequacy of NZ design code prescriptions for fundamental period, response spectra, floor acceleration and effective member stiffness. The efficacy of conventional numerical modelling assumptions for representing the UC Physics Building are examined using the observed building response. The numerical models comprise of the following: a one dimensional multi degree of freedom model, a two dimensional model along each axis of the building and a three dimensional model. Both moderate and strong ground motion records are used to examine the response and subsequently clarify the importance of linear and non-linear responses and the inclusion of base flexibility. The effects of soil-structure interaction are found to be significant in the transverse direction but not the longitudinal direction. Non-linear models predict minor in-elastic behaviour in both directions during the 4 September 2010 Mw 7.1 Darfield earthquake. The observed torsional response is found to be accurately captured by the three dimensional model by considering the interaction between the UC Physics Building and the adjacent structure. With the inclusion of adequate numerical modelling assumptions, the structural response is able to be predicted to within 10% for the majority of the earthquake events considered.

Research papers, University of Canterbury Library

The 22 February 2011, Mw6.2-6.3 Christchurch earthquake is the most costly earthquake to affect New Zealand, causing 181 fatalities and severely damaging thousands of residential and commercial buildings, and most of the city lifelines and infrastructure. This manuscript presents an overview of observed geotechnical aspects of this earthquake as well as some of the completed and on-going research investigations. A unique aspect, which is particularly emphasized, is the severity and spatial extent of liquefaction occurring in native soils. Overall, both the spatial extent and severity of liquefaction in the city was greater than in the preceding 4th September 2010 Darfield earthquake, including numerous areas that liquefied in both events. Liquefaction and lateral spreading, variable over both large and short spatial scales, affected commercial structures in the Central Business District (CBD) in a variety of ways including: total and differential settlements and tilting; punching settlements of structures with shallow foundations; differential movements of components of complex structures; and interaction of adjacent structures via common foundation soils. Liquefaction was most severe in residential areas located to the east of the CBD as a result of stronger ground shaking due to the proximity to the causative fault, a high water table approximately 1m from the surface, and soils with composition and states of high susceptibility and potential for liquefaction. Total and differential settlements, and lateral movements, due to liquefaction and lateral spreading is estimated to have severely compromised 15,000 residential structures, the majority of which otherwise sustained only minor to moderate damage directly due to inertial loading from ground shaking. Liquefaction also had a profound effect on lifelines and other infrastructure, particularly bridge structures, and underground services. Minor damage was also observed at flood stop banks to the north of the city, which were more severely impacted in the 4th September 2010 Darfield earthquake. Due to the large high-frequency ground motion in the Port hills numerous rock falls and landslides also occurred, resulting in several fatalities and rendering some residential areas uninhabitable.

Research papers, University of Canterbury Library

Rock mass defect controlled deep-seated landslides are widespread within the deeply incised landscapes formed in Tertiary soft rock terrain in New Zealand. The basal failure surfaces of deep-seated slope failures are defined by thin, comparatively weak and laterally continuous bedding parallel layers termed critical stratigraphic horizons. These horizons have a sedimentary origin and have typically experienced some prior tectonically induced shear displacement at the time of slope failure. The key controls on the occurrence and form of deep-seated landslides are considered in terms of rock mass defect properties and tectonic and climatic forcing. The selection of two representative catchments (in southern Hawke's Bay and North Canterbury) affected by tectonic and climatic forcing has shown that the spatial and temporal initiation of deep-seated bedrock landslides in New Zealand Tertiary soft rock terrain is a predictable rather than a stochastic process; and that deep-seated landslides as a mass wasting process have a controlling role in landscape evolution in many catchments formed in Tertiary soft rock terrain. The Ella Landslide in North Canterbury is a deep-seated (~85 m) translational block slide that has failed on a 5 - 10 mm thick, kaolinite-rich, pre-sheared critical stratigraphic horizon. The residual strength of this sedimentary horizon, (C'R 2.6 - 2.7 kPa, and Ѳ'R = 16 - 21°), compared to the peak strength of the dominant lithology (C' = 176 kPa, and Ѳ' = 37°) defines a high strength contrast in the succession, and therefore a critical location for the basal failure surface of deep-seated slope failures. The (early to mid Holocene) Ella Landslide debris formed a large landslide dam in the Kate Stream catchment and this has significantly retarded rates of mass wasting in the middle catchment. Numerical stability analysis shows that this slope failure would have most likely required the influence of earthquake induced strong ground motion and the event is tentatively correlated to a Holocene event on the Omihi Fault. The influence of this slope failure is likely to affect the geomorphic development of the catchment on a scale of 10⁴ - 10⁵ years. In deeply incised catchments at the southeastern margin of the Maraetotara Plateau, southern Hawke's Bay, numerous widespread deep-seated landslides have basal failure surfaces defined by critical stratigraphic horizons in the form of thin « 20 mm) tuffaceous beds in the Makara Formation flysch (alternating sandstone and mudstone units). The geometry of deep-seated slope failures is controlled by these regularly spaced (~70 m), very weak critical stratigraphic horizons (C'R 3.8 - 14.2 kPa, and Ѳ'R = 2 - 5°), and regularly spaced (~45 m) and steeply dipping (-50°) critical conjugate joint/fault sets, which act as slide block release surfaces. Numerical stability analysis and historical precedent show that the temporal initiation of deep-seated landslides is directly controlled by short term tectonic forcing in the form of periodic large magnitude earthquakes. Published seismic hazard data shows the recurrence interval of earthquakes producing strong ground motions of 0.35g at the study site is every 150 yrs, however, if subduction thrust events are considered the level of strong ground motion may be much higher. Multiple occurrences of deep-seated slope failure are correlated to failure on the same critical stratigraphic horizon, in some cases in three adjacent catchments. Failure on multiple critical stratigraphic horizons leads to the development of a "stepped" landscape morphology. This slope form will be maintained during successive accelerated stream incision events (controlled by long term tectonic and climatic forcing) for as long as catchments are developing in this specific succession. Rock mass defect controlled deep seated landslides are controlling catchment head progression, landscape evolution and hillslope morphology in the Hawke's Bay study area and this has significant implications for the development of numerical landscape evolution models of landscapes formed in similar strata. Whereas the only known numerical model to consider deep seated landslides as an erosion process (ZSCAPE) considers them as stochastic in time and space, this study shows that this could not be applied to a landscape where the widespread spatial occurrence of deep-seated landslides is controlled by rock mass defects. In both of the study areas for this project, and by implication in many catchments in Tertiary soft rock terrain, deep-seated landslides controlled by rock mass defect strength, spacing and orientation, and tectonic and climatic forcing have an underlying control on landscape evolution. This study quantifies parameters for the development of numerical landscape evolution models that would assess the role of specific parameters, such as uplift rates, incision rates and earthquake recurrence in catchment evolution in Tertiary soft rock terrain.

Research papers, University of Canterbury Library

Advanced seismic effective-stress analysis is used to scrutinize the liquefaction performance of 55 well-documented case-history sites from Christchurch. The performance of these sites during the 2010-2011 Canterbury earthquake sequence varied significantly, from no liquefaction manifestation at the ground surface (in any of the major events) to severe liquefaction manifestation in multiple events. For the majority of the 55 sites, the simplified liquefaction evaluation procedures, which are conventionally used in engineering practice, could not explain these dramatic differences in the manifestation. Detailed geotechnical characterization and subsequent examination of the soil profile characteristics of the 55 sites identified some similarities but also important differences between sites that manifested liquefaction in the two major events of the sequence (YY-sites) and sites that did not manifest liquefaction in either event (NN-sites). In particular, while the YY-sites and NN-sites are shown to have practically identical critical layer characteristics, they have significant differences with regard to their deposit characteristics including the thickness and vertical continuity of their critical zones and liquefiable materials. A CPT-based effective stress analysis procedure is developed and implemented for the analyses of the 55 case history sites. Key features of this procedure are that, on the one hand, it can be fully automated in a programming environment and, on the other hand, it is directly equivalent (in the definition of cyclic resistance and required input data) to the CPT-based simplified liquefaction evaluation procedures. These features facilitate significantly the application of effective-stress analysis for simple 1D free-field soil-column problems and also provide a basis for rigorous comparisons of the outcomes of effective-stress analyses and simplified procedures. Input motions for the analyses are derived using selected (reference) recordings from the two major events of the 2010-2011 Canterbury earthquake sequence. A step-by-step procedure for the selection of representative reference motions for each site and their subsequent treatment (i.e. deconvolution and scaling) is presented. The focus of the proposed procedure is to address key aspects of spatial variability of ground motion in the near-source region of an earthquake including extended-source effects, path effects, and variation in the deeper regional geology.

Research papers, University of Canterbury Library

The Amuri Earthquake of September 1, 1888 (magnitude M = 6.5 to 6.8) occurred on the Hope River Segment of the Hope Fault west of Hanmer Plains. The earthquake was felt strongly in North Canterbury and North Westland and caused considerable property damage and landsliding in the Lower Hope Valley. However, damage reports and the spatial distribution of felt intensities emphasize extreme variations in seismic effects over short distances, probably due to topographic focusing and local ground conditions. Significant variations in lateral fault displacement occurred at secondary fault segment boundaries (side-steps and bends in the fault trace) during the 1888 earthquake. This historical spatial variation in lateral slip is matched by the Late Quaternary geomorphic distribution of slip on the Hope River Segment of the Hope Fault. Trenching studies at two sites on the Hope Fault have also identified evidence for five pre-historic earthquakes of similar magnitude to the 1888 earthquake and an average recurrence interval of 134 ± 27 years between events. Magnitude estimates for the 1888 earthquake are combined with a. strong ground motion attenuation expression to provide an estimate of potential ground accelerations in Amuri District during-future earthquakes on the Hope River Segment of the Hope Fault. The predicted acceleration response on bedrock sites within 20 km of the epicentral region is between 0.23 g and 0.34 g. The close match between the historic, inferred pre-historic and geomorphic distribution of lateral slip indicates that secondary fault segmentation exerts a strong structural control on rupture propagation and the expression of fault displacement at the surface. In basement rocks at depth the spatial variations in slip are inferred to be distributed within zones of pervasive cataclastic shear, on either side of the fault segment boundaries. The large variations in surface displacement across fault segment boundaries means that one must know the geometry of the fault in order to evaluate slip-rates calculated from individual locations. The average Late Quaternary slip-rate on the Hope Fault at Glynn Wye Station is between 15.5 mm/yr and 18.25 mm/yr and the rate on the subsidiary Kakapo Fault is between 5.0 mm/yr and 7.5 mm/yr. These rates have been determined from sites which are relatively free of structural complication.

Research papers, University of Canterbury Library

Extended Direct Analysis (EDA), developed at the University of Canterbury, is an advance on the AISC Direct Analysis method for the analysis of frames subjected to static forces. EDA provides a faster, simple and more rational way to properly consider the second-order effects, initial residual stresses (IRS) and the initial imperfections or steel structures under one directional loading than conventional analysis methods. This research applied the EDA method to quantify the effect of member overstrength on frame behaviour for a single storey frame. Also, the effects of IRS, which were included in the EDA static analysis, but which are not considered explicitly in non-linear seismic analysis, were evaluated in two ways. Firstly, they were considered for simple structures subject to increasing cyclic displacement in different directions. Secondly, incremental dynamic analysis with realistic ground motion was used to quantify the likely effect of IRS in earthquakes. It was found that, contrary to traditional wisdom and practice, greater member strengths can result in lower frame strengths for frames under monotonic lateral loading. The structural lateral capacity of the overstrength case was reduced by 6% compared to the case using the dependable member strengths. Also, it resulted significantly different in member demands. Therefore, it is recommended that when either plastic analysis or EDA is used, that both upper and lower bounds on the likely member strength should be considered to determine the total frame strength and the member demands. Results of push-pull analysis under displacement control showed that for IRS ratio, gamma < 0.5 and axial compressive force ratio, N*/Ns, up to 0.5, IRS did affect the structural behaviour in the first half cycle. However, the behavior in the later cycles was not significantly affected. It also showed that the effect of initial residual stresses in the frame was less significant than for the column alone when the column was subjected to similar axial compressive force. The incremental dynamic analysis results from both cantilever column and the three-storey steel frame showed that by increasing gamma = 0 to 0.5, the effect of IRS on seismic responses, based on the 50% confidence level, was less than 3% for N*/Ns, up to 0.5.