Search

found 475 results

Images for model; more images...

Research papers, University of Canterbury Library

Liquefaction-induced lateral spreading in Christchurch and surrounding suburbs during the recent Canterbury Earthquake Sequence (2010-2011) caused significant damage to structures and lifelines located in close proximity to streams and rivers. Simplified methods used in current engineering practice for predicting lateral ground displacements exhibit a high degree of epistemic uncertainty, but provide ‘order of magnitude’ estimates to appraise the hazard. We wish to compare model predictions to field measurements in order to assess the model’s capabilities and limitations with respect to Christchurch conditions. The analysis presented focuses on the widely-used empirical model of Youd et al. (2002), developed based on multi-linear regression (MLR) of case history data from lateral spreading occurrence in Japan and the US. Two issues arising from the application of this model to Christchurch were considered: • Small data set of Standard Penetration Test (SPT) and soil gradation indices (fines content FC, and mean grain size, D50) required for input. We attempt to use widely available CPT data with site specific correlations to FC and D50. • Uncertainty associated with the model input parameters and their influence on predicted displacements. This has been investigated for a specific location through a sensitivity analysis.

Research papers, University of Canterbury Library

This study explores the nature of smaller businesses’ resilience following two major earthquakes that severely disrupted their place of doing business. Data from the owners of ten smaller businesses are qualitative and longitudinal, spanning the period 2011 through 2018, providing first-hand narrative accounts of their responses in the earthquakes’ aftermath. All ten owners showed some individual resilience; six businesses survived through to 2018, of which three have recovered strongly. All three owned their premises; operated business-tobusiness models; and were able to adapt and continue to follow path-extension strategies. All the other businesses had direct business-to-customer models operating from leased premises, typically in major retail malls. Four eventually recognised path-exhaustion at different times and so did not survive through to 2018. We conclude however that post-disaster recovery is best explained in terms of business model resilience. Even the most resilient of individual owners will struggle to survive if their business model is either not resilient or cannot be made so. Individual resilience is necessary but not sufficient.

Images, Alexander Turnbull Library

Scientists stand before a model that will predict earthquakes. The model is a large arm attached to wires and switches with a thumb that flicks coins '"Heads" we have big quake at five-o-clock... "tails" we don't'. Context: short term predictions (hours to days) are in general unlikely to be possible, at present. Relates to the Christchurch earthquakes which experts have said could go on for years. Quantity: 1 digital cartoon(s).

Research papers, University of Canterbury Library

This paper summarizes the development of a high-resolution surficial shear wave velocity model based on the combination of the large high-spatial-density database of cone penetration test (CPT) logs in and around Christchurch, New Zealand and a recently-developed Christchurch-specific empirical correlation between soil shear wave velocity and CPT. This near-surface shear wave velocity model has applications for site characterization efforts via the development of maps of time-averaged shear wave velocities over specific depths, as well as use in site response analysis and ground motion simulation.

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

Supplemental energy dissipation devices are increasingly used to protect structures, limit loads transferred to structural elements and absorbing significant response energy without sacrificial structural damage. Lead extrusion dampers are supplemental energy dissipation devices, where recent development of smaller volumetric size with high force capacities, called high force to volume (HF2V) devices, has seen deployment in a large series of scaled and full-scaled experiments, as well as in three new structures in Christchurch, NZ and San Francisco, USA. HF2V devices have previously been designed using limited precision models, so there is variation in force prediction capability. Further, while the overall resistive force is predicted, the knowledge of the relative contributions of the different internal reaction mechanisms to these overall resistive forces is lacking, limiting insight and predictive accuracy in device design. There is thus a major need for detailed design models to better understand force generation, and to aid precision device design. These outcomes would speed the overall design and implementation process for uptake and use, reducing the need for iterative experimental testing. Design parameters from 17 experimental HF2V device tests are used to create finite element models using ABAQUS. The analysis is run using ABAQUS Explicit, in multiple step times of 1 second with automatic increments, to balance higher accuracy and computational time. The output is obtained from the time- history output of the contact pressure forces including the normal and friction forces on the lead along the shaft. These values are used to calculate the resistive force on the shaft as it moves through the lead, and thus the device force. Results of these highly nonlinear, high strain analyses are compared to experimental device force results. Model errors compared to experimental results for all 17 devices ranged from 0% to 20% with a mean absolute error of 6.4%, indicating most errors were small. In particular, the standard error in manufacturing is SE = ±14%. In this case, 15 of 17 devices (88%) are within ±1SE (±14%) and 2 of 17 devices (12%) are within ±2SE (±28). These results show low errors and a distribution of errors compared to experimental results that are within experimental device construction variability. The overall modelling methodology is objective and repeatable, and thus generalizable. The exact same modelling approach is applied to all devices with only the device geometry changing. The results validate the overall approach with relatively low error, providing a general modelling methodology for accurate design of HF2V devices.

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

An extensive research program is on-going at the University of Canterbury, New Zealand to develop new technologies to permit the construction of multi-storey timber buildings in earthquake prone areas. The system combines engineered timber beams, columns and walls with ductile moment resisting connections using post-tensioned tendons and eventually energy dissipaters. The extensive experimental testing on post-tensioned timber building systems has proved a remarkable lateral response of the proposed solutions. A wide number of post-tensioned timber subassemblies, including beam-column connections, single or coupled walls and column-foundation connections, have been analysed in static or quasi-static tests. This contribution presents the results of the first dynamic tests carried out with a shake-table. Model frame buildings (3-storey and 5-storey) on one-quarter scale were tested on the shake-table to quantify the response of post-tensioned timber frames during real-time earthquake loading. Equivalent viscous damping values were computed for post-tensioned timber frames in order to properly predict their response using numerical models. The dynamic tests were then complemented with quasi-static push and pull tests performed to a 3-storey post-tensioned timber frame. Numerical models were included to compare empirical estimations versus dynamic and quasi-static experimental results. Different techniques to model the dynamic behaviour of post-tensioned timber frames were explored. A sensitivity analysis of alternative damping models and an examination of the influence of designer choices for the post-tensioning force and utilization of column armouring were made. The design procedure for post-tensioned timber frames was summarized and it was applied to two examples. Inter-storey drift, base shear and overturning moments were compared between numerical modelling and predicted/targeted design values.

Research Papers, Lincoln University

Globally, the maximum elevations at which treelines are observed to occur coincide with a 6.4 °C soil isotherm. However, when observed at finer scales, treelines display a considerable degree of spatial complexity in their patterns across the landscape and are often found occurring at lower elevations than expected relative to the global-scale pattern. There is still a lack of understanding of how the abiotic environment imposes constraints on treeline patterns, the scales at which different effects are acting, and how these effects vary over large spatial extents. In this thesis, I examined abrupt Nothofagus treelines across seven degrees of latitude in New Zealand in order to investigate two broad questions: (1) What is the nature and extent of spatial variability in Nothofagus treelines across the country? (2) How is this variation associated with abiotic variation at different spatial scales? A range of GIS, statistical, and atmospheric modelling methods were applied to address these two questions. First, I characterised Nothofagus treeline patterns at a 15x15km scale across New Zealand using a set of seven, GIS-derived, quantitative metrics that describe different aspects of treeline position, shape, spatial configuration, and relationships with adjacent vegetation. Multivariate clustering of these metrics revealed distinct treeline types that showed strong spatial aggregation across the country. This suggests a strong spatial structuring of the abiotic environment which, in turn, drives treeline patterns. About half of the multivariate treeline metric variation was explained by patterns of climate, substrate, topographic and disturbance variability; on the whole, climatic and disturbance factors were most influential. Second, I developed a conceptual model that describes how treeline elevation may vary at different scales according to three categories of effects: thermal modifying effects, physiological stressors, and disturbance effects. I tested the relevance of this model for Nothofagus treelines by investigating treeline elevation variation at five nested scales (regional to local) using a hierarchical design based on nested river catchments. Hierarchical linear modelling revealed that the majority of the variation in treeline elevation resided at the broadest, regional scale, which was best explained by the thermal modifying effects of solar radiation, mountain mass, and differences in the potential for cold air ponding. Nonetheless, at finer scales, physiological and disturbance effects were important and acted to modify the regional trend at these scales. These results suggest that variation in abrupt treeline elevations are due to both broad-scale temperature-based growth limitation processes and finer-scale stress- and disturbance-related effects on seedling establishment. Third, I explored the applicability of a meso-scale atmospheric model, The Air Pollution Model (TAPM), for generating 200 m resolution, hourly topoclimatic data for temperature, incoming and outgoing radiation, relative humidity, and wind speeds. Initial assessments of TAPM outputs against data from two climate station locations over seven years showed that the model could generate predictions with a consistent level of accuracy for both sites, and which agreed with other evaluations in the literature. TAPM was then used to generate data at 28, 7x7 km Nothofagus treeline zones across New Zealand for January (summer) and July (winter) 2002. Using mixed-effects linear models, I determined that both site-level factors (mean growing season temperature, mountain mass, precipitation, earthquake intensity) and local-level landform (slope and convexity) and topoclimatic factors (solar radiation, photoinhibition index, frost index, desiccation index) were influential in explaining variation in treeline elevation within and among these sites. Treelines were generally closer to their site-level maxima in regions with higher mean growing season temperatures, larger mountains, and lower levels of precipitation. Within sites, higher treelines were associated with higher solar radiation, and lower photoinhibition and desiccation index values, in January, and lower desiccation index values in July. Higher treelines were also significantly associated with steeper, more convex landforms. Overall, this thesis shows that investigating treelines across extensive areas at multiple study scales enables the development of a more comprehensive understanding of treeline variability and underlying environmental constraints. These results can be used to formulate new hypotheses regarding the mechanisms driving treeline formation and to guide the optimal choice of field sites at which to test these hypotheses.

Research papers, University of Canterbury Library

This paper summarizes the development of a region-wide surficial shear wave velocity model based on the combination of the large high-spatial-density database of cone penetration test (CPT) logs in and around Christchurch, New Zealand and a recently-developed Christchurch-specific empirical correlation between soil shear wave velocity and CPT. The ongoing development of this near-surface shear wave velocity model has applications for site characterization efforts via the development of maps of time-averaged shear wave velocities over specific depths, and the identification of regional similarities and differences in soil shear stiffness.

Research papers, University of Canterbury Library

This poster aims to present fragility functions for pipelines buried in liquefaction-prone soils. Existing fragility models used to quantify losses can be based on old data or use complex metrics. Addressing these issues, the proposed functions are based on the Christchurch network and soil and utilizes the Canterbury earthquake sequence (CES) data, partially represented in Figure 1. Figure 1 (a) presents the pipe failure dataset, which describes the date, location and pipe on which failures occurred. Figure 1 (b) shows the simulated ground motion intensity median of the 22nd February 2011 earthquake. To develop the model, the network and soil characteristics have also been utilized

Research papers, University of Canterbury Library

This document reviews research-based understandings of the concept of resilience. A conceptual model is developed which identifies a number of the factors that influence individual and household resilience. Guided by the model, a series of recommendations are developed for practices that will support individual and household resilience in Canterbury in the aftermath of the 2010-2011 earthquakes.

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

The purpose of this thesis is to conduct a detailed examination of the forward-directivity characteristics of near-fault ground motions produced in the 2010-11 Canterbury earthquakes, including evaluating the efficacy of several existing empirical models which form the basis of frameworks for considering directivity in seismic hazard assessment. A wavelet-based pulse classification algorithm developed by Baker (2007) is firstly used to identify and characterise ground motions which demonstrate evidence of forward-directivity effects from significant events in the Canterbury earthquake sequence. The algorithm fails to classify a large number of ground motions which clearly exhibit an early-arriving directivity pulse due to: (i) incorrect pulse extraction resulting from the presence of pulse-like features caused by other physical phenomena; and (ii) inadequacy of the pulse indicator score used to carry out binary pulse-like/non-pulse-like classification. An alternative ‘manual’ approach is proposed to ensure 'correct' pulse extraction and the classification process is also guided by examination of the horizontal velocity trajectory plots and source-to-site geometry. Based on the above analysis, 59 pulse-like ground motions are identified from the Canterbury earthquakes , which in the author's opinion, are caused by forward-directivity effects. The pulses are also characterised in terms of their period and amplitude. A revised version of the B07 algorithm developed by Shahi (2013) is also subsequently utilised but without observing any notable improvement in the pulse classification results. A series of three chapters are dedicated to assess the predictive capabilities of empirical models to predict the: (i) probability of pulse occurrence; (ii) response spectrum amplification caused by the directivity pulse; (iii) period and amplitude (peak ground velocity, PGV) of the directivity pulse using observations from four significant events in the Canterbury earthquakes. Based on the results of logistic regression analysis, it is found that the pulse probability model of Shahi (2013) provides the most improved predictions in comparison to its predecessors. Pulse probability contour maps are developed to scrutinise observations of pulses/non-pulses with predicted probabilities. A direct comparison of the observed and predicted directivity amplification of acceleration response spectra reveals the inadequacy of broadband directivity models, which form the basis of the near-fault factor in the New Zealand loadings standard, NZS1170.5:2004. In contrast, a recently developed narrowband model by Shahi & Baker (2011) provides significantly improved predictions by amplifying the response spectra within a small range of periods. The significant positive bias demonstrated by the residuals associated with all models at longer vibration periods (in the Mw7.1 Darfield and Mw6.2 Christchurch earthquakes) is likely due to the influence of basin-induced surface waves and non-linear soil response. Empirical models for the pulse period notably under-predict observations from the Darfield and Christchurch earthquakes, inferred as being a result of both the effect of nonlinear site response and influence of the Canterbury basin. In contrast, observed pulse periods from the smaller magnitude June (Mw6.0) and December (Mw5.9) 2011 earthquakes are in good agreement with predictions. Models for the pulse amplitude generally provide accurate estimates of the observations at source-to-site distances between 1 km and 10 km. At longer distances, observed PGVs are significantly under-predicted due to their slower apparent attenuation. Mixed-effects regression is employed to develop revised models for both parameters using the latest NGA-West2 pulse-like ground motion database. A pulse period relationship which accounts for the effect of faulting mechanism using rake angle as a continuous predictor variable is developed. The use of a larger database in model development, however does not result in improved predictions of pulse period for the Darfield and Christchurch earthquakes. In contrast, the revised model for PGV provides a more appropriate attenuation of the pulse amplitude with distance, and does not exhibit the bias associated with previous models. Finally, the effects of near-fault directivity are explicitly included in NZ-specific probabilistic seismic hazard analysis (PSHA) using the narrowband directivity model of Shahi & Baker (2011). Seismic hazard analyses are conducted with and without considering directivity for typical sites in Christchurch and Otira. The inadequacy of the near-fault factor in the NZS1170.5: 2004 is apparent based on a comparison with the directivity amplification obtained from PSHA.

Research papers, University of Canterbury Library

Well-validated liquefaction constitutive models are increasingly important as non-linear time history analyses become relatively more common in industry for key projects. Previous validation efforts of PM4Sand, a plasticity model specifically for liquefaction, have generally focused on centrifuge tests; however, pore pressure transducers installed at several free-field sites during the Canterbury Earthquake Sequence (CES) in Christchurch, New Zealand provide a relatively unique dataset to validate against. This study presents effective stress site response analyses performed in the finite difference software FLAC to examine the capability of PM4Sand to capture the generation of excess pore pressures during earthquakes. The characterization of the subsurface is primarily based on extensive cone penetration tests (CPT) carried out in Christchurch. Correlations based on penetration resistances are used to estimate soil parameters, such as relative density and shear wave velocity, which affect liquefaction behaviour. The resulting free-field FLAC model is used to estimate time histories of excess pore pressure, which are compared with records during several earthquakes in the CES to assess the suitability of PM4Sand.

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

Seismic behaviour of typical unreinforced masonry (URM) brick houses, that were common in early last century in New Zealand and still common in many developing countries, is experimentally investigated at University of Canterbury, New Zealand in this research. A one halfscale model URM house is constructed and tested under earthquake ground motions on a shaking table. The model structure with aspect ratio of 1.5:1 in plan was initially tested in the longitudinal direction for several earthquakes with peak ground acceleration (PGA) up to 0.5g. Toppling of end gables (above the eaves line) and minor to moderate cracking around window and door piers was observed in this phase. The structure was then rotated 90º and tested in the transverse (short) direction for ground motions with PGA up to 0.8g. Partial out-of-plane failure of the face loaded walls in the second storey and global rocking of the model was observed in this phase. A finite element analysis and a mechanism analysis are conducted to assess the dynamic properties and lateral strength of the model house. Seismic fragility function of URM houses is developed based on the experimental results. Damping at different phases of the response is estimated using an amplitude dependent equivalent viscous damping model. Financial risk of similar URM houses is then estimated in term of expected annual loss (EAL) following a probabilistic financial risk assessment framework. Risks posed by different levels of damage and by earthquakes of different frequencies are then examined.

Research papers, University of Canterbury Library

© 2017 The Royal Society of New Zealand. This paper discusses simulated ground motion intensity, and its underlying modelling assumptions, for great earthquakes on the Alpine Fault. The simulations utilise the latest understanding of wave propagation physics, kinematic earthquake rupture descriptions and the three-dimensional nature of the Earth's crust in the South Island of New Zealand. The effect of hypocentre location is explicitly examined, which is found to lead to significant differences in ground motion intensities (quantified in the form of peak ground velocity, PGV) over the northern half and southwest of the South Island. Comparison with previously adopted empirical ground motion models also illustrates that the simulations, which explicitly model rupture directivity and basin-generated surface waves, lead to notably larger PGV amplitudes than the empirical predictions in the northern half of the South Island and Canterbury. The simulations performed in this paper have been adopted, as one possible ground motion prediction, in the ‘Project AF8’ Civil Defence Emergency Management exercise scenario. The similarity of the modelled ground motion features with those observed in recent worldwide earthquakes as well as similar simulations in other regions, and the notably higher simulated amplitudes than those from empirical predictions, may warrant a re-examination of regional impact assessments for major Alpine Fault earthquakes.

Research papers, University of Canterbury Library

In 2010 and 2011 Christchurch, New Zealand experienced a series of earthquakes that caused extensive damage across the city, but primarily to the Central Business District (CBD) and eastern suburbs. A major feature of the observed damage was extensive and severe soil liquefaction and associated ground damage, affecting buildings and infrastructure. The behaviour of soil during earthquake loading is a complex phenomena that can be most comprehensively analysed through advanced numerical simulations to aid engineers in the design of important buildings and critical facilities. These numerical simulations are highly dependent on the capabilities of the constitutive soil model to replicate the salient features of sand behaviour during cyclic loading, including liquefaction and cyclic mobility, such as the Stress-Density model. For robust analyses advanced soil models require extensive testing to derive engineering parameters under varying loading conditions for calibration. Prior to this research project little testing on Christchurch sands had been completed, and none from natural samples containing important features such as fabric and structure of the sand that may be influenced by the unique stress-history of the deposit. This research programme is focussed on the characterisation of Christchurch sands, as typically found in the CBD, to facilitate advanced soil modelling in both res earch and engineering practice - to simulate earthquake loading on proposed foundation design solutions including expensive ground improvement treatments. This has involved the use of a new Gel Push (GP) sampler to obtain undisturbed samples from below the ground-water table. Due to the variable nature of fluvial deposition, samples with a wide range of soil gradations, and accordingly soil index properties, were obtained from the sampling sites. The quality of the samples is comprehensively examined using available data from the ground investigation and laboratory testing. A meta-quality assessment was considered whereby a each method of evaluation contributed to the final quality index assigned to the specimen. The sampling sites were characterised with available geotechnical field-based test data, primarily the Cone Penetrometer Test (CPT), supported by borehole sampling and shear-wave velocity testing. This characterisation provides a geo- logical context to the sampling sites and samples obtained for element testing. It also facilitated the evaluation of sample quality. The sampling sites were evaluated for liquefaction hazard using the industry standard empirical procedures, and showed good correlation to observations made following the 22 February 2011 earthquake. However, the empirical method over-predicted liquefaction occurrence during the preceding 4 September 2010 event, and under-predicted for the subsequent 13 June 2011 event. The reasons for these discrepancies are discussed. The response of the GP samples to monotonic and cyclic loading was measured in the laboratory through triaxial testing at the University of Canterbury geomechanics laboratory. The undisturbed samples were compared to reconstituted specimens formed in the lab in an attempt to quantify the effect of fabric and structure in the Christchurch sands. Further testing of moist tamped re- constituted specimens (MT) was conducted to define important state parameters and state-dependent properties including the Critical State Line (CSL), and the stress-strain curve for varying state index. To account for the wide-ranging soil gradations, selected representative specimens were used to define four distinct CSL. The input parameters for the Stress-Density Model (S-D) were derived from a suite of tests performed on each representative soil, and with reference to available GP sample data. The results of testing were scrutinised by comparing the data against expected trends. The influence of fabric and structure of the GP samples was observed to result in similar cyclic strength curves at 5 % Double Amplitude (DA) strain criteria, however on close inspection of the test data, clear differences emerged. The natural samples exhibited higher compressibility during initial loading cycles, but thereafter typically exhibited steady growth of plastic strain and excess pore water pressure towards and beyond the strain criteria and initial liquefaction, and no flow was observed. By contrast the reconstituted specimens exhibited a stiffer response during initial loading cycles, but exponential growth in strains and associated excess pore water pressure beyond phase-transformation, and particularly after initial liquefaction where large strains were mobilised in subsequent cycles. These behavioural differences were not well characterised by the cyclic strength curve at 5 % DA strain level, which showed a similar strength for both GP samples and MT specimens. A preliminary calibration of the S-D model for a range of soil gradations is derived from the suite of laboratory test data. Issues encountered include the influence of natural structure on the peak-strength–state index relationship, resulting in much higher peak strengths than typically observed for sands in the literature. For the S-D model this resulted in excessive stiffness to be modelled during cyclic mobility, when the state index becomes large momentarily, causing strain development to halt. This behaviour prevented modelling the observed re- sponse of silty sands to large strains, synonymous with “liquefaction”. Efforts to reduce this effect within the current formulation are proposed as well as future research to address this issue.

Research papers, University of Canterbury Library

Natural catastrophes are increasing worldwide. They are becoming more frequent but also more severe and impactful on our built environment leading to extensive damage and losses. Earthquake events account for the smallest part of natural events; nevertheless seismic damage led to the most fatalities and significant losses over the period 1981-2016 (Munich Re). Damage prediction is helpful for emergency management and the development of earthquake risk mitigation projects. Recent design efforts focused on the application of performance-based design engineering where damage estimation methodologies use fragility and vulnerability functions. However, the approach does not explicitly specify the essential criteria leading to economic losses. There is thus a need for an improved methodology that finds the critical building elements related to significant losses. The here presented methodology uses data science techniques to identify key building features that contribute to the bulk of losses. It uses empirical data collected on site during earthquake reconnaissance mission to train a machine learning model that can further be used for the estimation of building damage post-earthquake. The first model is developed for Christchurch. Empirical building damage data from the 2010-2011 earthquake events is analysed to find the building features that contributed the most to damage. Once processed, the data is used to train a machine-learning model that can be applied to estimate losses in future earthquake events.

Research papers, University of Canterbury Library

This poster presents work to date on ground motion simulation validation and inversion for the Canterbury, New Zealand region. Recent developments have focused on the collection of different earthquake sources and the verification of the SPECFEM3D software package in forward and inverse simulations. SPECFEM3D is an open source software package which simulates seismic wave propagation and performs adjoint tomography based upon the spectral-element method. Figure 2: Fence diagrams of shear wave velocities highlighting the salient features of the (a) 1D Canterbury velocity model, and (b) 3D Canterbury velocity model. Figure 5: Seismic sources and strong motion stations in the South Island of New Zealand, and corresponding ray paths of observed ground motions. Figure 3: Domain used for the 19th October 2010 Mw 4.8 case study event including the location of the seismic source and strong motion stations. By understanding the predictive and inversion capabilities of SPECFEM3D, the current 3D Canterbury Velocity Model can be iteratively improved to better predict the observed ground motions. This is achieved by minimizing the misfit between observed and simulated ground motions using the built-in optimization algorithm. Figure 1 shows the Canterbury Velocity Model domain considered including the locations of small-to-moderate Mw events [3-4.5], strong motion stations, and ray paths of observed ground motions. The area covered by the ray paths essentially indicates the area of the model which will be most affected by the waveform inversion. The seismic sources used in the ground motion simulations are centroid moment tensor solutions obtained from GeoNet. All earthquake ruptures are modelled as point sources with a Gaussian source time function. The minimum Mw limit is enforced to ensure good signal-to-noise ratio and well constrained source parameters. The maximum Mw limit is enforced to ensure the point source approximation is valid and to minimize off-fault nonlinear effects.

Research papers, University of Canterbury Library

The aim of this thesis was to examine the spatial and the temporal patterns of anxiety and chest pain resulting from the Canterbury, New Zealand earthquaeks. Three research objectives were identified: examine any spatial or termporal clusters of anxiety and chest pain; examine the associations between anxiety, chest pain and damage to neighbourhood; and determine any statistically significant difference in counts of anxiety and chest pain after each earthquake or aftershock which resulted in severe damage. Measures of the extent of liquefaction the location of CERA red-zones were used as proxy measures for earthquake damage. Cases of those who presented to Christchurch Public Hospital Emergency Department with either anxiety or chest pain between May 2010 and April 2012 were aggregated to census area unit (CAU) level for analysis. This thesis has taken a unique approach to examining the spatial and spatio-temporal variations of anxiety and chest pain after an earthquake and offers unique results. This is the first study of its kind to use a GIS approach when examining Canterbury specific earthquake damage and health variables at a CAU level after the earthquakes. Through the use of spatio-termporal scan modelling, negative and linear regression modelling and temporal linear modelling with dummy variables this research was able to conclude there are significant spatial and temporal variations in anxiety and chest pain resulting from the earthquakes. The spatio-termporal scan modelling identified a hot cluster of both anxiety and chest pain within Christchurch at the same time the earthquakes occurred. The negative binomial model found liquefaction to be a stronger predictor of anxiety than the Canterbury Earthquake Recovery Authority's (CERA) land zones. The linear regression model foun chest pain to be positively associated with all measures of earthquake damage with the exception of being in the red-zone. The temporal modelling identified a significant increase in anxiety cases one month after a major earthquake, and chest pain cases spiked two weeks after an earthquake and gradually decreased over the following five weeks. This research was limited by lack of control period data, limited measures of earthquake damage, ethical restrictions, and the need for population tracking data. The findings of this research will be useful in the planning and allocation of mental wellbeing resources should another similar event like the Canterbury Earthquakes occur in New Zealand.

Research papers, University of Canterbury Library

Recent global tsunami events have highlighted the importance of effective tsunami risk management strategies (including land-use planning, structural and natural defences, warning systems, education and evacuation measures). However, the rarity of tsunami means that empirical data concerning reactions to tsunami warnings and tsunami evacuation behaviour is rare when compared to findings about evacuations to avoid other sources of hazard. To date empirical research into tsunami evacuations has focused on evacuation rates, rather than other aspects of the evacuation process. More knowledge is required about responses to warnings, pre-evacuation actions, evacuation dynamics and the return home after evacuations. Tsunami evacuation modelling has the potential to inform evidence-based tsunami risk planning and response. However to date tsunami evacuation models have largely focused on timings of evacuations, rather than evacuation behaviours. This Masters research uses a New Zealand case study to reduce both of these knowledge gaps. Qualitative survey data was gathered from populations across coastal communities in Banks Peninsula and Christchurch, New Zealand, required to evacuate due to the tsunami generated by the November 14th 2016 Kaikōura Earthquake. Survey questions asked about reactions to tsunami warnings, actions taken prior to evacuating and movements during the 2016 tsunami evacuation. This data was analysed to characterise trends and identify factors that influenced evacuation actions and behaviour. Finally, it was used to develop an evacuation model for Banks Peninsula. Where appropriate, the modelling inputs were informed by the survey data. Three key findings were identified from the results of the evacuation behaviour survey. Although 38% of the total survey respondents identified the earthquake shaking as a natural cue for the tsunami, most relied on receiving official warnings, including sirens, to prompt evacuations. Respondents sought further official information to inform their evacuation decisions, with 39% of respondents delaying their evacuation in order to do so. Finally, 96% of total respondents evacuated by car. This led to congestion, particularly in more densely populated Christchurch city suburbs. Prior to this research, evacuation modelling had not been completed for Banks Peninsula. The results of the modelling showed that if evacuees know how to respond to tsunami warnings and where and how to evacuate, there are no issues. However, if there are poor conditions, including if people do not evacuate immediately, if there are issues with the roading network, or if people do not know where or how to evacuate, evacuation times increase with there being more bottlenecks leading out of the evacuation zones. The results of this thesis highlight the importance of effective tsunami education and evacuation planning. Reducing exposure to tsunami risk through prompt evacuation relies on knowledge of how to interpret tsunami warnings, and when, where and how to evacuate. Recommendations from this research outline the need for public education and engagement, and the incorporation of evacuation signage, information boards and evacuation drills. Overall these findings provide more comprehensive picture of tsunami evacuation behaviour and decision making based on empirical data from a recent evacuation, which can be used to improve tsunami risk management strategies. This empirical data can also be used to inform evacuation modelling to improve the accuracy and realism of the evacuation models.

Research papers, The University of Auckland Library

The influence of nonlinear soil-foundation-structure interaction (SFSI) on the performance of multi-storey buildings during earthquake events has become increasingly important in earthquake resistant design. For buildings on shallow foundations, SFSI refers to nonlinear geometric effects associated with uplift of the foundation from the supporting soil as well as nonlinear soil deformation effects. These effects can potentially be beneficial for structural performance, reducing forces transmitted from ground shaking to the structure. However, there is also the potential consequence of residual settlement and rotation of the foundation. This Thesis investigates the influence of SFSI in the performance of multi-storey buildings on shallow foundations through earthquake observations, experimental testing, and development of spring-bed numerical models that can be incorporated into integrated earthquake resistant design procedures. Observations were made following the 22 February 2011 Christchurch Earthquake in New Zealand of a number of multi-storey buildings on shallow foundations that performed satisfactorily. This was predominantly the case in areas where shallow foundations, typically large raft foundations, were founded on competent gravel and where there was no significant manifestation of liquefaction at the ground surface. The properties of these buildings and the soils they are founded on directed experimental work that was conducted to investigate the mechanisms by which SFSI may have influenced the behaviour of these types of structure-foundation systems. Centrifuge experiments were undertaken at the University of Dundee, Scotland using a range of structure-foundation models and a layer of dense cohesionless soil to simulate the situation in Christchurch where multi-storey buildings on shallow foundations performed well. Three equivalent single degree of freedom (SDOF) models representing 3, 5, and 7 storey buildings with identical large raft foundations were subjected to a range of dynamic Ricker wavelet excitations and Christchurch Earthquake records to investigate the influence of SFSI on the response of the equivalent buildings. The experimental results show that nonlinear SFSI has a significant influence on structural response and overall foundation deformations, even though the large raft foundations on competent soil meant that there was a significant reserve of bearing capacity available and nonlinear deformations may have been considered to have had minimal effect. Uplift of the foundation from the supporting soil was observed across a wide range of input motion amplitudes and was particularly significant as the amplitude of motion increased. Permanent soil deformation represented by foundation settlement and residual rotation was also observed but mainly for the larger input motions. However, the absolute extent of uplift and permanent soil deformation was very small compared to the size of the foundation meaning the serviceability of the building would still likely be maintained during large earthquake events. Even so, the small extent of SFSI resulted in attenuation of the response of the structure as the equivalent period of vibration was lengthened and the equivalent damping in the system increased. The experimental work undertaken was used to validate and enhance numerical modelling techniques that are simple yet sophisticated and promote interaction between geotechnical and structural specialists involved in the design of multi-storey buildings. Spring-bed modelling techniques were utilised as they provide a balance between ease of use, and thus ease of interaction with structural specialists who have these techniques readily available in practice, and theoretically rigorous solutions. Fixed base and elastic spring-bed models showed they were unable to capture the behaviour of the structure-foundation models tested in the centrifuge experiments. SFSI spring-bed models were able to more accurately capture the behaviour but recommendations were proposed for the parameters used to define the springs so that the numerical models closely matched experimental results. From the spring-bed modelling and results of centrifuge experiments, an equivalent linear design procedure was proposed along with a procedure and recommendations for the implementation of nonlinear SFSI spring-bed models in practice. The combination of earthquake observations, experimental testing, and simplified numerical analysis has shown how SFSI is influential in the earthquake performance of multi-storey buildings on shallow foundations and should be incorporated into earthquake resistant design of these structures.

Research papers, University of Canterbury Library

Recent advances in timber design at the University of Canterbury have led to new structural systems that are appropriate for a wide range of building types, including multi-storey commercial office structures. These buildings are competitive with more traditional construction materials in terms of cost, sustainability and structural performance. This paper provides seismic design recommendations and analytical modelling approaches, appropriate for the seismic design of post-tensioned coupled timber wall systems. The models are based on existing seismic design theory for precast post-tensioned concrete, modified to more accurately account for elastic deformation of the timber wall systems and the influence of the floor system. Experimental test data from a two storey post-tensioned timber building, designed, constructed and tested at the University of Canterbury is used to validate the analytical models.