Search

found 29 results

Research papers, University of Canterbury Library

Using case studies from the 2010-2011 Canterbury, New Zealand earthquake sequence, this study assesses the accuracies of paleoliquefaction back-analysis methods and explores the challenges, techniques, and uncertainties associated with their application. While liquefaction-based back-analyses have been widely used to estimate the magnitudes of paleoearthquakes, their uncertain efficacies continue to significantly affect the computed seismic hazard in regions where they are relied upon. Accordingly, their performance is evaluated herein using liquefaction data from modern earthquakes with known magnitudes. It is shown that when the earthquake source location and mechanism are known, back-analysis methods are capable of accurately deriving seismic parameters from liquefaction evidence. However, because the source location and mechanism are often unknown in paleoseismic studies, and because accurate interpretation is shown to be more difficult in such cases, new analysis techniques are proposed herein. An objective parameter is proposed to geospatially assess the likelihood of any provisional source location, enabling an analyst to more accurately estimate the magnitude of a liquefaction-inducing paleoearthquake. This study demonstrates the application of back-analysis methods, provides insight into their potential accuracies, and provides a framework for performing paleoliquefaction analyses worldwide.

Research papers, University of Canterbury Library

This paper presents the ongoing development of a new 3D seismic velocity model of Canterbury, New Zealand. The model explicitly represents the Canterbury sedimentary basin, and other significant geologic horizons, which are expected to have important implications on observed ground motions. The model utilizes numerous sources of data, including 3D regional tomography with a variable-depth inferred Moho, seismic reflection survey lines, geotechnical boreholes and well logs, spectral analysis of surface waves, and CPT logs which provide velocity constraints over their respective ranges of application. The model provides P- and S-wave velocity and density (i.e. Vp, Vs and p) over a grid of input points, and is presently being utilized in broadband ground motion simulations of the 2010-2011 Canterbury earthquakes. Comparison of simulated ground motions with those observed in the 2010-2011 Canterbury earthquakes will help provide a better understanding of the salient physical processes which characterized the unique set of strong ground motions recorded in this sequence of earthquake events.

Research papers, University of Canterbury Library

This paper develops representative ground motion ensembles for several major earthquake scenarios in New Zealand. Cases considered include representative ground motions for the occurrence of Alpine, Hope, and Porters Pass earthquakes in Christchurch, and the occurrence of Wellington, Wairarapa, and Ohariu, fault ruptures in Wellington. Challenges in the development of ground motion ensembles for subduction zone earthquakes are also highlighted. The ground motions are selected based on the generalized conditional intensity measure (GCIM) approach, ensuring that the ground motion ensembles represent both the mean, and distribution of ground motion intensity which such scenarios could impose. These scenario-based ground motion sets can be used to complement ground motions which are often selected in conjunction with probabilistic seismic hazard analysis, in order to understand the performance of structures for the question “what if this fault ruptures?”

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

The abundance of cone penetration test (CPT) data from subsurface explorations in Christchurch and the surrounding areas provides a useful source of information for a characterization of the near surface shear wave velocity ( ) profile for the region. A portion of the investigations were conducted using seismic CPT, enabling the comparison of measured shear wave velocity with CPT data, and subsequently the evaluation of existing CPT- correlations for applicability to Canterbury-specific soils. The existing correlations are shown to be biased, generally over-predicting the observed with depth, thus demonstrating the need for a Canterbury-specific CPT- correlation.

Research papers, University of Canterbury Library

The 2010–2011 Canterbury earthquake sequence began with the 4 September 2010, Mw7.1 Darfield earthquake and includes up to ten events that induced liquefaction. Most notably, widespread liquefaction was induced by the Darfield and Mw6.2 Christchurch earthquakes. The combination of well-documented liquefaction response during multiple events, densely recorded ground motions for the events, and detailed subsurface characterization provides an unprecedented opportunity to add well-documented case histories to the liquefaction database. This paper presents and applies 50 high-quality cone penetration test (CPT) liquefaction case histories to evaluate three commonly used, deterministic, CPT-based simplified liquefaction evaluation procedures. While all the procedures predicted the majority of the cases correctly, the procedure proposed by Idriss and Boulanger (2008) results in the lowest error index for the case histories analyzed, thus indicating better predictions of the observed liquefaction response.

Research papers, University of Canterbury Library

The performance of conventionally designed reinforced concrete (RC) structures during the 2011 Christchurch earthquake has demonstrated that there is greater uncertainty in the seismic performance of RC components than previously understood. RC frame and wall structures in the Christchurch central business district were observed to form undesirable cracks patterns in the plastic hinge region while yield penetration either side of cracks, and into development zones, were less than theoretical predictions. The implications of this unexpected behaviour: (i) significantly less available ductility; (ii) less hysteretic energy dissipation; and (iii) the localization of peak reinforcement strains, results in considerable doubt for the residual capacity of RC structures. The significance of these consequences has prompted a review of potential sources of uncertainty in seismic experimentation with the intention to improve the current confidence level for newly designed conventional RC structures. This paper attempts to revisit the principles of RC mechanics, in particular, to consider the influence of loading history, concrete tensile strength, and reinforcement ratio on the performance of ‘real’ RC structures compared to experimental test specimens.

Research papers, University of Canterbury Library

The paper presents preliminary findings from comprehensive research studies on the liquefaction-induced damage to buildings and infrastructure in Christchurch during the 2010-2011 Canterbury earthquakes. It identifies key factors and mechanisms of damage to road bridges, shallow foundations of CBD buildings and buried pipelines, and highlights the implications of the findings for the seismic analysis and design of these structures.

Research papers, University of Canterbury Library

This is an interim report from the research study performed within the NHRP Research Project “Impacts of soil liquefaction on land, buildings and buried pipe networks: geotechnical evaluation and design, Project 3: Seismic assessment and design of pipe networks in liquefiable soils”. The work presented herein is a continuation of the comprehensive study on the impacts of Christchurch earthquakes on the buried pipe networks presented in Cubrinovski et al. (2011). This report summarises the performance of Christchurch City’s potable water, waste water and road networks through the 2010-2011 Canterbury Earthquake Sequence (CES), and particularly focuses on the potable water network. It combines evidence based on comprehensive and well-documented data on the damage to the water network, detailed observations and interpretation of liquefaction-induced land damage, records and interpretations of ground motion characteristics induced by the Canterbury earthquakes, for a network analysis and pipeline performance evaluation using a GIS platform. The study addresses a range of issues relevant in the assessment of buried networks in areas affected by strong earthquakes and soil liquefaction. It discusses performance of different pipe materials (modern flexible pipelines and older brittle pipelines) including effects of pipe diameters, fittings and pipeline components/details, trench backfill characteristics, and severity of liquefaction. Detailed breakdown of key factors contributing to the damage to buried pipes is given with reference to the above and other relevant parameters. Particular attention is given to the interpretation, analysis and modelling of liquefaction effects on the damage and performance of the buried pipe networks. Clear link between liquefaction severity and damage rate for the pipeline has been observed with an increasing damage rate seen with increasing liquefaction severity. The approach taken here was to correlate the pipeline damage to LRI (Liquefaction Resistance Index, newly developed parameter in Cubrinovski et al., 2011) which represents a direct measure for the soil resistance to liquefaction while accounting for the seismic demand through PGA. Key quality of the adopted approach is that it provides a general methodology that in conjunction with conventional methods for liquefaction evaluation can be applied elsewhere in New Zealand and internationally. Preliminary correlations between pipeline damage (breaks km-1), liquefaction resistance (LRI) and seismic demand (PGA) have been developed for AC pipes, as an example. Such correlations can be directly used in the design and assessment of pipes in seismic areas both in liquefiable and non-liquefiable areas. Preliminary findings on the key factors for the damage to the potable water pipe network and established empirical correlations are presented including an overview of the damage to the waste water and road networks but with substantially less detail. A comprehensive summary of the damage data on the buried pipelines is given in a series of appendices.

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

Generalized conditional intensity measure (GCIM) method is extended to ground motion selection for scenario ruptures. Using different rupture scenarios and site conditions, various aspects of the GCIM methodology are scrutinized, including: (i) implementation of different weight vectors and the composition of the IM vector; (ii) quantifying the importance of replicate selections for different number of desired ground motions; and (iii) the effect of considering bounds on the implicit causal parameters of the prospective ground motions. Using the extended methodology, representative ground motion ensembles for several major earthquake scenarios in New Zealand are developed. Cases considered include representative ground motions for the occurrence of Alpine, Hope, and Porters Pass earthquakes in Christchurch city, and the occurrence of Wellington, Wairarapa, and Ohariu fault ruptures in Wellington city. Challenges in the development of ground motion ensembles for subduction zone earthquakes are also highlighted. The selected scenario-based ground motion sets can be used to complement ground motions which are often selected in conjunction with probabilistic seismic hazard analysis, in order to understand the performance of structures for the question “what if this fault ruptures?”

Research papers, University of Canterbury Library

Generalized conditional intensity measure (GCIM) method is extended to ground motion selection for scenario ruptures. Using different rupture scenarios and site conditions, various aspects of the GCIM methodology are scrutinized, including: (i) implementation of different weight vectors and the composition of the IM vector; (ii) quantifying the importance of replicate selections for different number of desired ground motions; and (iii) the effect of considering bounds on the implicit causal parameters of the prospective ground motions. Using the extended methodology, representative ground motion ensembles for several major earthquake scenarios in New Zealand are developed. Cases considered include representative ground motions for the occurrence of Alpine, Hope, and Porters Pass earthquakes in Christchurch city, and the occurrence of Wellington, Wairarapa, and Ohariu fault ruptures in Wellington city. Challenges in the development of ground motion ensembles for subduction zone earthquakes are also highlighted. The selected scenario-based ground motion sets can be used to complement ground motions which are often selected in conjunction with probabilistic seismic hazard analysis, in order to understand the performance of structures for the question “what if this fault ruptures?”

Research papers, University of Canterbury Library

In recent Canterbury earthquakes, structures have performed well in terms of life safety but the estimated total cost of the rebuild was as high as $40 billion. The major contributors to this cost are repair/demolition/rebuild cost, the resulting downtime and business interruption. For this reason, the authors are exploring alternate building systems that can minimize the downtime and business interruption due to building damage in an earthquake; thereby greatly reducing the financial implications of seismic events. In this paper, a sustainable and demountable precast reinforced concrete (RC) frame system in which the precast members are connected via steel tubes/plates or steel angles/plates and high strength friction grip (HSFG) bolts is introduced. In the proposed system, damaged structural elements in seismic frames can be easily replaced with new ones; thereby making it an easily and quickly repairable and a low-loss system. The column to foundation connection in the proposed system can be designed either as fixed or pinned depending on the requirement of strength and stiffness. In a fixed base frame system, ground storey columns will also be damaged along with beams in seismic events, which are to be replaced after seismic events; whereas in a pin base frame only beams (which are easy to replace) will be damaged. Low to medium rise (3-6 storey) precast RC frame buildings with fixed and pin bases are analyzed in this paper; and their lateral capacity, lateral stiffness and natural period are scrutinized to better understand the pros and cons of the demountable precast frame system with fixed and pin base connections.

Research papers, University of Canterbury Library

In the last century, seismic design has undergone significant advancements. Starting from the initial concept of designing structures to perform elastically during an earthquake, the modern seismic design philosophy allows structures to respond to ground excitations in an inelastic manner, thereby allowing damage in earthquakes that are significantly less intense than the largest possible ground motion at the site of the structure. Current performance-based multi-objective seismic design methods aim to ensure life-safety in large and rare earthquakes, and to limit structural damage in frequent and moderate earthquakes. As a result, not many recently built buildings have collapsed and very few people have been killed in 21st century buildings even in large earthquakes. Nevertheless, the financial losses to the community arising from damage and downtime in these earthquakes have been unacceptably high (for example; reported to be in excess of 40 billion dollars in the recent Canterbury earthquakes). In the aftermath of the huge financial losses incurred in recent earthquakes, public has unabashedly shown their dissatisfaction over the seismic performance of the built infrastructure. As the current capacity design based seismic design approach relies on inelastic response (i.e. ductility) in pre-identified plastic hinges, it encourages structures to damage (and inadvertently to incur loss in the form of repair and downtime). It has now been widely accepted that while designing ductile structural systems according to the modern seismic design concept can largely ensure life-safety during earthquakes, this also causes buildings to undergo substantial damage (and significant financial loss) in moderate earthquakes. In a quest to match the seismic design objectives with public expectations, researchers are exploring how financial loss can be brought into the decision making process of seismic design. This has facilitated conceptual development of loss optimisation seismic design (LOSD), which involves estimating likely financial losses in design level earthquakes and comparing against acceptable levels of loss to make design decisions (Dhakal 2010a). Adoption of loss based approach in seismic design standards will be a big paradigm shift in earthquake engineering, but it is still a long term dream as the quantification of the interrelationships between earthquake intensity, engineering demand parameters, damage measures, and different forms of losses for different types of buildings (and more importantly the simplification of the interrelationship into design friendly forms) will require a long time. Dissecting the cost of modern buildings suggests that the structural components constitute only a minor portion of the total building cost (Taghavi and Miranda 2003). Moreover, recent research on seismic loss assessment has shown that the damage to non-structural elements and building contents contribute dominantly to the total building loss (Bradley et. al. 2009). In an earthquake, buildings can incur losses of three different forms (damage, downtime, and death/injury commonly referred as 3Ds); but all three forms of seismic loss can be expressed in terms of dollars. It is also obvious that the latter two loss forms (i.e. downtime and death/injury) are related to the extent of damage; which, in a building, will not just be constrained to the load bearing (i.e. structural) elements. As observed in recent earthquakes, even the secondary building components (such as ceilings, partitions, facades, windows parapets, chimneys, canopies) and contents can undergo substantial damage, which can lead to all three forms of loss (Dhakal 2010b). Hence, if financial losses are to be minimised during earthquakes, not only the structural systems, but also the non-structural elements (such as partitions, ceilings, glazing, windows etc.) should be designed for earthquake resistance, and valuable contents should be protected against damage during earthquakes. Several innovative building technologies have been (and are being) developed to reduce building damage during earthquakes (Buchanan et. al. 2011). Most of these developments are aimed at reducing damage to the buildings’ structural systems without due attention to their effects on non-structural systems and building contents. For example, the PRESSS system or Damage Avoidance Design concept aims to enable a building’s structural system to meet the required displacement demand by rocking without the structural elements having to deform inelastically; thereby avoiding damage to these elements. However, as this concept does not necessarily reduce the interstory drift or floor acceleration demands, the damage to non-structural elements and contents can still be high. Similarly, the concept of externally bracing/damping building frames reduces the drift demand (and consequently reduces the structural damage and drift sensitive non-structural damage). Nevertheless, the acceleration sensitive non-structural elements and contents will still be very vulnerable to damage as the floor accelerations are not reduced (arguably increased). Therefore, these concepts may not be able to substantially reduce the total financial losses in all types of buildings. Among the emerging building technologies, base isolation looks very promising as it seems to reduce both inter-storey drifts and floor accelerations, thereby reducing the damage to the structural/non-structural components of a building and its contents. Undoubtedly, a base isolated building will incur substantially reduced loss of all three forms (dollars, downtime, death/injury), even during severe earthquakes. However, base isolating a building or applying any other beneficial technology may incur additional initial costs. In order to provide incentives for builders/owners to adopt these loss-minimising technologies, real-estate and insurance industries will have to acknowledge the reduced risk posed by (and enhanced resilience of) such buildings in setting their rental/sale prices and insurance premiums.

Research papers, University of Canterbury Library

The Canterbury earthquakes, which involved widespread damage in the February 2011 event and ongoing aftershocks near the Christchurch central business district (CBD), presented decision-makers with many recovery challenges. This paper identifies major government decisions, challenges, and lessons in the early recovery of Christchurch based on 23 key-informant interviews conducted 15 months after the February 2011 earthquake. It then focuses on one of the most important decisions – maintaining the cordon around the heavily damaged CBD – and investigates its impacts. The cordon displaced 50,000 central city jobs, raised questions about (and provided new opportunities for) the long-term viability of downtown, influenced the number and practice of building demolitions, and affected debris management; despite being associated with substantial losses, the cordon was commonly viewed as necessary, and provided some benefits in facilitating recovery. Management of the cordon poses important lessons for planning for catastrophic urban earthquakes around the world.

Research papers, University of Canterbury Library

The majority of Christchurch’s stormwater has historically been discharged untreated directly into urban surface waterways. These receiving waterways have become adversely affected by the contaminants carried in the stormwater, particularly sediment and heavy metals. An event-based contaminant load model was developed to identify the distribution and magnitude of contaminant loads entering the waterway, as well as to assess the reduction in TSS and heavy metal loads that can be achieved by various stormwater management options. The GIS-Excel based model estimates contaminant loads from an individual storm event based on different contributing impervious surfaces and key rainfall characteristics (rainfall intensity, duration, pH and antecedent dry days). It then calculates contaminant reduction loads that could be achieved through source reduction (e.g. green roofs, repainting) as well as from treatment (e.g. raingardens, wet ponds) applied to different surfaces within the catchment. This model differs from other annual load models as it is event-based and accounts for storm characteristics in its calculation of contaminant loads. Christchurch is a valuable case setting due the unique opportunity for retrofitting improved stormwater management in the post-earthquake rebuild. It is anticipated that this modelling approach could later be adapted for use in other urban settings outside of Christchurch.

Research papers, University of Canterbury Library

This presentation discusses recent empirical ground motion modelling efforts in New Zealand. Firstly, the active shallow crustal and subduction interface and slab ground motion prediction equations (GMPEs) which are employed in the 2010 update of the national seismic hazard model (NSHM) are discussed. Other NZ-specific GMPEs developed, but not incorporated in the 2010 update are then discussed, in particular, the active shallow crustal model of Bradley (2010). A brief comparison of the NZ-specific GMPEs with the near-source ground motions recorded in the Canterbury earthquakes is then presented, given that these recordings collectively provide a significant increase in observed strong motions in the NZ catalogue. The ground motion prediction expert elicitation process that was undertaken following the Canterbury earthquakes for active shallow crustal earthquakes is then discussed. Finally, ongoing GMPE-related activities are discussed including: ground motion and metadata database refinement, improved site characterization of strong motion station, and predictions for subduction zone earthquakes.

Research papers, University of Canterbury Library

Deformational properties of soil, in terms of modulus and damping, exert a great influence on seismic response of soil sites. However, these properties for sands containing some portion of fines particles have not been systematically addressed. In addition, simultaneous modelling of the modulus and damping behaviour of soils during cyclic loading is desirable. This study presents an experimental and computational investigation into the deformational properties of sands containing fines content in the context of site response analysis. The experimental investigation is carried on sandy soils sourced from Christchurch, New Zealand using a dynamic triaxial apparatus while the computational aspect is based on the framework of total-stress one-dimensional (1D) cyclic behaviour of soil. The experimental investigation focused on a systematic study on the deformational behaviour of sand with different amounts of fines content (particle diameter ≤ 75µm) under drained conditions. The silty sands were prepared by mixing clean sand with three different percentages of fines content. A series of bender element tests at small-strain range and stress-controlled dynamic triaxial tests at medium to high-strain ranges were conducted on samples of clean sand and silty sand. This allowed measurements of linear and nonlinear deformational properties of the same specimen for a wide strain range. The testing program was designed to quantify the effects of void ratio and fines content on the low-strain stiffness of the silty sand as well as on the nonlinear stress-strain relationship and corresponding shear modulus and damping properties as a function of cyclic shear strains. Shear wave velocity, Vs, and maximum shear modulus, Gmax, of silty sand was shown to be significantly smaller than the respective values for clean sands measured at the same void ratio, e, or same relative density, Dr. However, the test results showed that the difference in the level of nonlinearity between clean sand and silty sands was small. For loose samples prepared at an identical relative density, the behaviour of clean sand was slightly less nonlinear as compared to sandy soils with higher fines content. This difference in the nonlinear behaviour of clean sand and sandy soils was negligible for dense soils. Furthermore, no systematic influence of fines content on the material damping curve was observed for sands with fines content FC = 0 to 30%. In order to normalize the effects of fines on moduli of sands, equivalent granular void ratio, e*, was employed. This was done through quantifying the participation of fines content in the force transfer chain of the sand matrix. As such, a unified framework for modelling of the variability of shear wave velocity, Vs, (or shear modulus, Gmax) with void ratio was achieved for clean sands and sands with fines, irrespective of their fines content. Furthermore, modelling of the cyclic stress-strain behaviour based on this experimental program was investigated. The modelling effort focused on developing a simple constitutive model which simultaneously models the soil modulus and damping relationships with shear strains observed in laboratory tests. The backbone curve of the cyclic model was adopted based on a modified version of Kondner and Zelasko (MKZ) hyperbolic function, with a curvature coefficient, a. In order to simulate the hysteretic cycles, the conventional Masing rules (Pyke 1979) were revised. The parameter n, in the Masing’s criteria was assumed to be a function of material damping, h, measured in the laboratory. As such the modulus and damping produced by the numerical model could match the stress-strain behaviour observed in the laboratory over the course of this study. It was shown that the Masing parameter n, is strain-dependent and generally takes values of n ≤ 2. The model was then verified through element test simulations under different cyclic loadings. It was shown that the model could accurately simulate the modulus and the damping simultaneously. The model was then incorporated within the OpenSees computational platform and was used to scrutinize the effects of damping on one-dimensional seismic site response analysis. For this purpose, several strong motion stations which recorded the Canterbury earthquake sequence were selected. The soil profiles were modelled as semi-infinite horizontally layered deposits overlying a uniform half-space subjected to vertically propagating shear waves. The advantages and limitations of the nonlinear model in terms of simulating soil nonlinearity and associated material damping were further scrutinized. It was shown that generally, the conventional Masing criteria unconservatively may underestimate some response parameters such as spectral accelerations. This was shown to be due to larger hysteretic damping modelled by using conventional Masing criteria. In addition, maximum shear strains within the soil profiles were also computed smaller in comparison to the values calculated by the proposed model. Further analyses were performed to study the simulation of backbone curve beyond the strain ranges addressed in the experimental phase of this study. A key issue that was identified was that relying only on the modulus reduction curves to simulate the stress-strain behaviour of soil may not capture the actual soil strength at larger strains. Hence, strength properties of the soil layer should also be incorporated to accurately simulate the backbone curve.

Research papers, University of Canterbury Library

Seismically vulnerable buildings constitute a major problem for the safety of human beings. In many parts of the world, reinforced concrete (RC) frame buildings designed and constructed with substandard detailing, no consideration of capacity design principles, and improper or no inclusion of the seismic actions, have been identified. Amongst those vulnerable building, one particular typology representative of the construction practice of the years previous to the 1970’s, that most likely represents the worst case scenario, has been widely investigated in the past. The deficiencies of that building typology are related to non-ductile detailing in beam column joints such as the use of plain round bars, the lack of stirrups inside the joint around the longitudinal reinforcement of the column, the use of 180° end hooks in the beams, the use of lap splices in potential ‘plastic hinge’ regions, and substandard quality of the materials. That type of detailing and the lack of a capacity design philosophy create a very fragile fuse in the structure where brittle inelastic behaviour is expected to occur, which is the panel zone region of exterior beam column joints. The non-ductile typology described above was extensively investigated at the University of Canterbury in the context of the project ‘Retrofit Solutions for New Zealand Multi-Storey Buildings’ (2004-2011), founded by the ‘Foundation for Research, Science and Technology’ Tūāpapa Rangahau Pūtaiao. The experimental campaign prior to the research carried out by the author consisted of quasi-static tests of beam column joint subassemblies subjected to lateral loading regime, with constant and varying axial load in the column. Most of those specimens were representative of a plane 2D frame (knee joint), while others represented a portion of a space 3D frame (corner joints), and only few of them had a floor slab, transverse beams, and lap splices. Using those experiments, several feasible, cost-effective, and non-invasive retrofit techniques were developed, improved, and refined. Nevertheless, the slow motion nature of those experiments did not take into account the dynamical component inherent to earthquake related problems. Amongst the set of techniques investigated, the use of FRP layers for strengthening beam column joints is of particular interest due to its versatility and the momentum that its use has gained in the current state of the practice. That particular retrofit technique was previously used to develop a strengthening scheme suitable for plane 2D and space 3D corner beam column joints, but lacking of floor slabs. In addition, a similar scheme was not developed for exterior joints of internal frames, referred here as ‘cruciform’. In this research a 2/5 scale RC frame model building comprising of two frames in parallel (external and internal) joined together by means of floor slabs and transverse beams, with non-ductile characteristics identical to those of the specimens investigated previously by others, and also including lap splices, was developed. In order to investigate the dynamic response of that building, a series of shake table tests with different ground motions were performed. After the first series of tests, the specimen was modified by connecting the spliced reinforcement in the columns in order to capture a different failure mode. Ground motions recorded during seismic events that occurred during the initial period of the experimental campaign (2010) were used in the subsequent experiments. The hierarchy of strengths and sequence of events in the panel zone region were evaluated in an extended version of the bending moment-axial load (M-N) performance domain developed by others. That extension was required due to the asymmetry in the beam cross section introduced by the floor slab. In addition, the effect of the torsion resistance provided by the spandrel (transverse beam) was included. In order to upgrade the brittle and unstable performance of the as-built/repaired specimen, a practical and suitable ad-hoc FRP retrofit intervention was developed, following a partial retrofit strategy that aimed to strengthen exterior beam column joints only (corner and cruciform). The ability of the new FRP scheme to revert the sequence of events in the panel zone region was evaluated using the extended version of the M-N performance domain as well as the guidelines for strengthening plane joints developed by others. Weakening of the floor slab in a novel configuration was also incorporated with the purpose of reducing the flexural capacity of the beam under negative bending moment (slab in tension), enabling the damage relocation from the joint into the beam. The efficacy of the developed retrofit intervention in upgrading the seismic performance of the as-built specimen was investigated using shake table tests with the input motions used in the experiments of the as-built/repaired specimen. Numerical work aimed to predict the response of the model building during the most relevant shake table tests was carried out. By using a simple numerical model with concentrated plasticity elements constructed in Ruaumoko2D, the results of blind and post-experimental predictions of the response of the specimen were addressed. Differences in the predicted response of the building using the nominal and the actual recorded motions of the shake table were investigated. The dependence of the accuracy of the numerical predictions on the assumed values of the parameters that control the hysteresis rules of key structural members was reviewed. During the execution of the experimental campaign part of this thesis, two major earthquakes affected the central part of Chile (27 of February 2010 Maule earthquake) and the Canterbury region in New Zealand (22 February 2011 Canterbury earthquake), respectively. As the author had the opportunity to experience those events and investigate their consequences in structures, the observations related to non-ductile detailing and drawbacks in the state of the practice related to reinforced concrete walls was also addressed in this research, resulting in preliminary recommendations for the refinement of current seismic code provisions and assessment guidelines. The investigations of the ground motions recorded during those and other earthquakes were used to review the procedures related to the input motions used for nonlinear dynamic analysis of buildings as required by most of the current code provisions. Inelastic displacement spectra were constructed using ground motions recorded during the earthquakes mentioned above, in order to investigate the adequacy of modification factors used to obtain reduced design spectra from elastic counterparts. Finally a simplified assessment procedure for RC walls that incorporates capacity compatible spectral demands is proposed.

Research papers, University of Canterbury Library

Since the early 1980s seismic hazard assessment in New Zealand has been based on Probabilistic Seismic Hazard Analysis (PSHA). The most recent version of the New Zealand National Seismic Hazard Model, a PSHA model, was published by Stirling et al, in 2012. This model follows standard PSHA principals and combines a nation-wide model of active faults with a gridded point-source model based on the earthquake catalogue since 1840. These models are coupled with the ground-motion prediction equation of McVerry et al (2006). Additionally, we have developed a time-dependent clustering-based PSHA model for the Canterbury region (Gerstenberger et al, 2014) in response to the Canterbury earthquake sequence. We are now in the process of revising that national model. In this process we are investigating several of the fundamental assumptions in traditional PSHA and in how we modelled hazard in the past. For this project, we have three main focuses: 1) how do we design an optimal combination of multiple sources of information to produce the best forecast of earthquake rates in the next 50 years: can we improve upon a simple hybrid of fault sources and background sources, and can we better handle the uncertainties in the data and models (e.g., fault segmentation, frequency-magnitude distributions, time-dependence & clustering, low strain-rate areas, and subduction zone modelling)? 2) developing revised and new ground-motion predictions models including better capturing of epistemic uncertainty – a key focus in this work is developing a new strong ground motion catalogue for model development; and 3) how can we best quantify if changes we have made in our modelling are truly improvements? Throughout this process we are working toward incorporating numerical modelling results from physics based synthetic seismicity and ground-motion models.

Research papers, University of Canterbury Library

Recent experiences from the Darfield and Canterbury, New Zealand earthquakes have shown that the soft soil condition of saturated liquefiable sand has a profound effect on seismic response of buildings, bridges and other lifeline infrastructure. For detailed evaluation of seismic response three dimensional integrated analysis comprising structure, foundation and soil is required; such an integrated analysis is referred to as Soil Foundation Structure Interaction (SFSI) in literatures. SFSI is a three-dimensional problem because of three primary reasons: first, foundation systems are three-dimensional in form and geometry; second, ground motions are three-dimensional, producing complex multiaxial stresses in soils, foundations and structure; and third, soils in particular are sensitive to complex stress because of heterogeneity of soils leading to a highly anisotropic constitutive behaviour. In literatures the majority of seismic response analyses are limited to plane strain configuration because of lack of adequate constitutive models both for soils and structures, and computational limitation. Such two-dimensional analyses do not represent a complete view of the problem for the three reasons noted above. In this context, the present research aims to develop a three-dimensional mathematical formulation of an existing plane-strain elasto-plastic constitutive model of sand developed by Cubrinovski and Ishihara (1998b). This model has been specially formulated to simulate liquefaction behaviour of sand under ground motion induced earthquake loading, and has been well-validated and widely implemented in verifcation of shake table and centrifuge tests, as well as conventional ground response analysis and evaluation of case histories. The approach adopted herein is based entirely on the mathematical theory of plasticity and utilises some unique features of the bounding surface plasticity formalised by Dafalias (1986). The principal constitutive parameters, equations, assumptions and empiricism of the existing plane-strain model are adopted in their exact form in the three-dimensional version. Therefore, the original two-dimensional model can be considered as a true subset of the three-dimensional form; the original model can be retrieved when the tensorial quantities of the three dimensional version are reduced to that of the plane-strain configuration. Anisotropic Drucker-Prager type failure surface has been adopted for the three-dimensional version to accommodate triaxial stress path. Accordingly, a new mixed hardening rule based on Mroz’s approach of homogeneous surfaces (Mroz, 1967) has been introduced for the virgin loading surface. The three-dimensional version is validated against experimental data for cyclic torsional and triaxial stress paths.

Research papers, University of Canterbury Library

Structural engineering is facing an extraordinarily challenging era. These challenges are driven by the increasing expectations of modern society to provide low-cost, architecturally appealing structures which can withstand large earthquakes. However, being able to avoid collapse in a large earthquake is no longer enough. A building must now be able to withstand a major seismic event with negligible damage so that it is immediately occupiable following such an event. As recent earthquakes have shown, the economic consequences of not achieving this level of performance are not acceptable. Technological solutions for low-damage structural systems are emerging. However, the goal of developing a low-damage building requires improving the performance of both the structural skeleton and the non-structural components. These non-structural components include items such as the claddings, partitions, ceilings and contents. Previous research has shown that damage to such items contributes a disproportionate amount to the overall economic losses in an earthquake. One such non-structural element that has a history of poor performance is the external cladding system, and this forms the focus of this research. Cladding systems are invariably complicated and provide a number of architectural functions. Therefore, it is important than when seeking to improve their seismic performance that these functions are not neglected. The seismic vulnerability of cladding systems are determined in this research through a desktop background study, literature review, and postearthquake reconnaissance survey of their performance in the 2010 – 2011 Canterbury earthquake sequence. This study identified that precast concrete claddings present a significant life-safety risk to pedestrians, and that the effect they have upon the primary structure is not well understood. The main objective of this research is consequently to better understand the performance of precast concrete cladding systems in earthquakes. This is achieved through an experimental campaign and numerical modelling of a range of precast concrete cladding systems. The experimental campaign consists of uni-directional, quasi static cyclic earthquake simulation on a test frame which represents a single-storey, single-bay portion of a reinforced concrete building. The test frame is clad with various precast concrete cladding panel configurations. A major focus is placed upon the influence the connection between the cladding panel and structural frame has upon seismic performance. A combination of experimental component testing, finite element modelling and analytical derivation is used to develop cladding models of the cladding systems investigated. The cyclic responses of the models are compared with the experimental data to evaluate their accuracy and validity. The comparison shows that the cladding models developed provide an excellent representation of real-world cladding behaviour. The cladding models are subsequently applied to a ten-storey case-study building. The expected seismic performance is examined with and without the cladding taken into consideration. The numerical analyses of the case-study building include modal analyses, nonlinear adaptive pushover analyses, and non-linear dynamic seismic response (time history) analyses to different levels of seismic hazard. The clad frame models are compared to the bare frame model to investigate the effect the cladding has upon the structural behaviour. Both the structural performance and cladding performance are also assessed using qualitative damage states. The results show a poor performance of precast concrete cladding systems is expected when traditional connection typologies are used. This result confirms the misalignment of structural and cladding damage observed in recent earthquake events. Consequently, this research explores the potential of an innovative cladding connection. The outcomes from this research shows that the innovative cladding connection proposed here is able to achieve low-damage performance whilst also being cost comparable to a traditional cladding connection. It is also theoretically possible that the connection can provide a positive value to the seismic performance of the structure by adding addition strength, stiffness and damping. Finally, the losses associated with both the traditional and innovative cladding systems are compared in terms of tangible outcomes, namely: repair costs, repair time and casualties. The results confirm that the use of innovative cladding technology can substantially reduce the overall losses that result from cladding damage.

Research papers, University of Canterbury Library

In most design codes, infill walls are considered as non-structural elements and thus are typically neglected in the design process. The observations made after major earthquakes (Duzce 1999, L’Aquila 2009, Christchurch 2011) have shown that even though infill walls are considered to be non-structural elements, they interact with the structural system during seismic actions. In the case of heavy infill walls (i.e. clay brick infill walls), the whole behaviour of the structure may be affected by this interaction (i.e. local or global structural failures such as soft storey mechanism). In the case of light infill walls (i.e. non-structural drywalls), this may cause significant economical losses. To consider the interaction of the structural system with the ‘non-structural ’infill walls at design stage may not be a practical approach due to the complexity of the infill wall behaviour. Therefore, the purpose of the reported research is to develop innovative technological solutions and design recommendations for low damage non-structural wall systems for seismic actions by making use of alternative approaches. Light (steel/timber framed drywalls) and heavy (unreinforced clay brick) non-structural infill wall systems were studied by following an experimental/numerical research programme. Quasi-static reverse cyclic tests were carried out by utilizing a specially designed full scale reinforced concrete frame, which can be used as a re-usable bare frame. In this frame, two RC beams and two RC columns were connected by two un-bonded post tensioning bars, emulating a jointed ductile frame system (PRESSS technology). Due to the rocking behaviour at the beam-column joint interfaces, this frame was typically a low damage structural solution, with the post-tensioning guaranteeing a linear elastic behaviour. Therefore, this frame could be repeatedly used in all of the tests carried out by changing only the infill walls within this frame. Due to the linear elastic behaviour of this structural bare frame, it was possible to extract the exact behaviour of the infill walls from the global results. In other words, the only parameter that affected the global results was given by the infill walls. For the test specimens, the existing practice of construction (as built) for both light and heavy non-structural walls was implemented. In the light of the observations taken during these tests, modified low damage construction practices were proposed and tested. In total, seven tests were carried out: 1) Bare frame , in order to confirm its linear elastic behaviour. 2) As built steel framed drywall specimen FIF1-STFD (Light) 3) As built timber framed drywall specimen FIF2-TBFD (Light) 4) As built unreinforced clay brick infill wall specimen FIF3-UCBI (Heavy) 5) Low damage steel framed drywall specimen MIF1-STFD (Light) 6) Low damage timber framed drywall specimen MIF2-TBFD (Light) 7) Low damage unreinforced clay brick infill wall specimen MIF5-UCBI (Heavy) The tests of the as built practices showed that both drywalls and unreinforced clay brick infill walls have a low serviceability inter-storey drift limit (0.2-0.3%). Based on the observations, simple modifications and details were proposed for the low damage specimens. The details proved to be working effectively in lowering the damage and increasing the serviceability drift limits. For drywalls, the proposed low damage solutions do not introduce additional cost, material or labour and they are easily applicable in real buildings. For unreinforced clay brick infill walls, a light steel sub-frame system was suggested that divides the infill panel zone into smaller individual panels, which requires additional labour and some cost. However, both systems can be engineered for seismic actions and their behaviour can be controlled by implementing the proposed details. The performance of the developed details were also confirmed by the numerical case study analyses carried out using Ruaumoko 2D on a reinforced concrete building model designed according to the NZ codes/standards. The results have confirmed that the implementation of the proposed low damage solutions is expected to significantly reduce the non-structural infill wall damage throughout a building.

Research papers, University of Canterbury Library

Active faults capable of generating highly damaging earthquakes may not cause surface rupture (i.e., blind faults) or cause surface ruptures that evade detection due to subsequent burial or erosion by surface processes. Fault populations and earthquake frequency-­‐magnitude distributions adhere to power laws, implying that faults too small to cause surface rupture but large enough to cause localized strong ground shaking densely populate continental crust. The rupture of blind, previously undetected faults beneath Christchurch, New Zealand in a suite of earthquakes in 2010 and 2011, including the fatal 22 February 2011 moment magnitude (Mw) 6.2 Christchurch earthquake and other large aftershocks, caused a variety of environmental impacts, including major rockfall, severe liquefaction, and differential surface uplift and subsidence. All of these effects occurred where geologic evidence for penultimate effects of the same nature existed. To what extent could the geologic record have been used to infer the presence of proximal, blind and / or unidentified faults near Christchurch? In this instance, we argue that phenomena induced by high intensity shaking, such as rock fragmentation and rockfall, revealed the presence of proximal active faults in the Christchurch area prior to the recent earthquake sequence. Development of robust earthquake shaking proxy datasets should become a higher scientific priority, particularly in populated regions.

Research papers, University of Canterbury Library

This paper presents the preliminary findings of a study on the resilience and recovery of organisations following the Darfield earthquake in New Zealand on 4 September 2010. Sampling included organisations proximal and distal to the fault trace, organisations located within central business districts, and organisations from seven diverse industry sectors. The research captured information on the challenges to, the impacts on, and the reflections of the organisations in the first months of recovery. Organisations in central business districts and in the hospitality sector were most likely to close while organisations that had perishable stock and livestock were more heavily reliant on critical services. Staff well-being, cash flow, and customer loss were major concerns for organisations across all sectors. For all organisations, the most helpful factors in mitigating the effects of the earthquake to be their relationship with staff, the design and type of buildings, and critical service continuity or swift reinstatement of services.