Search

found 193 results

Images for Loads; more images...

Research papers, University of Canterbury Library

Post-tensioned timber technology was originally developed and researched at the University of Canterbury (UC) in New Zealand in 2005. It can provide a low-damage seismic design solution for multi-storey mass timber buildings. Since mass timber products, such as cross-laminated timber (CLT), have high in-plane stiffness, a post-tensioned timber shear wall will deform mainly in a rocking mechanism. The moment capacity of the wall at the base is commonly determined using the elastic form of the Modified Monolithic Beam Analogy (MMBA). In the calculation of the moment capacity at the wall base, it is critical to accurately predict the location of the neutral axis and the timber compressive stress distribution. Three 2/3 scale 8.6m tall post-tensioned CLT walls were experimentally tested under quasi-static cyclic loading – both uni-directional and bi-directional- in this study. These specimens included a single wall, a coupled wall, and a C-shaped core-wall. The main objective was to develop post-tensioned C-shaped timber core-walls for tall timber buildings with enhanced lateral strength and stiffness. To better understand the timber compressive stress distributions at the wall base, particle tracking technology (PTT) technology was applied for the first time to investigate the behaviour of the compression toe. Previous post-tensioned timber testing primarily used the displacement measurements to determine the timber compressive behavior at the wall base or rocking interfaces. However, by using PTT technology, the timber strain measurements in the compression zone can be much more accurate as PTT is able to track the movement of many particles on the timber surface. This paper presents experimental testing results of post-tensioned CLT walls with a focus on capturing timber compressive behavior using PTT. The PTT measurements were able to better capture small base rotations which occurred at the onset of gap opening and capture unexpected phenomena in core-wall tests. The single wall test result herein presented indicates that while the MMBA could predict the moment rotation behavior with reasonable accuracy, the peak strain response was under predicted in the compression toe. Further detailed study is required to better understand the complex strain fields generated reflective of the inherent cross-thickness inhomogeneity and material variability of CLT.

Research papers, University of Canterbury Library

To reduce seismic vulnerability and the economic impact of seismic structural damage, it is important to protect structures using supplemental energy dissipation devices. Several types of supplemental damping systems can limit loads transferred to structures and absorb significant response energy without sacrificial structural damage. Lead extrusion dampers are one type of supplemental energy dissipation devices. A smaller volumetric size with high force capacities, called high force to volume (HF2V) devices, have been employed in a large series of scaled and full-scaled experiments, as well as in three new structures in Christchurch and San Francisco. HF2V devices have previously been designed using very simple models with limited precision. They are then manufactured, and tested to ensure force capacities match design goals, potentially necessitating reassembly or redesign if there is large error. In particular, devices with a force capacity well above or below a design range can require more testing and redesign, leading to increased economic and time cost. Thus, there is a major need for a modelling methodology to accurately estimate the range of possible device force capacity values in the design phase – upper and lower bounds. Upper and lower bound force capacity estimates are developed from equations in the metal extrusion literature. These equations consider both friction and extrusion forces between the lead and the bulged shaft in HF2V devices. The equations for the lower and upper bounds are strictly functions of device design parameters ensuring easy use in the design phase. Two different sets of estimates are created, leading to estimates for the lower and upper bounds denoted FLB,1, FUB,1, FUB,2, respectively. The models are validated by comparing the bounds with experimental force capacity data from 15 experimental HF2V device tests. All lower bound estimates are below or almost equal to the experimental device forces, and all upper bound estimates are above. Per the derivation, the (FLB,1, FUB,1) pair provide narrower bounds. The (FLB,1, FUB,1) pair also had a mean lower bound gap of -34%, meaning the lower bound was 74% of device force on average, while the mean upper bound gap for FUB,1 was +23%. These are relatively tight bounds, within ~±2 SE of device manufacture, and can be used as a guide to ensure device forces are in range for the actual design use when manufactured. Therefore, they provide a useful design tool.

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

According to TS 1170.5, designing a building to satisfy code-prescribed criteria (e.g., drift limit, member safety, P-Δ stability) at the ultimate limit state and relying on the inherent margins within the design code would lead to an acceptable mean annual frequency of collapse (λ꜀) in the range of 10−⁴ to 10−⁵. Modern performance objectives, such as λ꜀ and expected annual loss (EAL), are not explicitly considered. Although buckling-restrained braced frame (BRBF) buildings were widely adopted as lateral load-resisting systems for office and car park buildings in the Christchurch rebuild following the Canterbury earthquakes in New Zealand, there are currently no official guidelines for their design. The primary focus of this study is to develop a risk-targeted design framework for BRBF buildings that can achieve the performance objectives desired by stakeholders. To this extent, key factors influencing λ꜀ and EAL of BRBF buildings are identified. These factors include gusset plate design, number of storeys, design drift limit, BRBF beam-column connection, brace configuration, brace angle, brace material grade, and analysis method (equivalent lateral force vs. modal response spectrum). A novel 3D BRBF modelling approach capable of simulating out-of-plane buckling failure of buckling-restrained brace (BRB) gusset plates is developed. Prior experimental studies on sub-assemblies conducted elsewhere have demonstrated that gusset plates and end zones may buckle out of plane prematurely, before BRBs reach their maximum axial compression load carrying capacity. Current 2D BRBF macro models, typically used in research, cannot simulate this failure mode. A conventional 2D BRBF model underestimates the λ꜀ of a case-study 4-storey super-X configured steel BRBF building (designed according to NZS-3404) by a factor of two compared to the estimate from the proposed 3D model. These findings suggest that the current NZS-3404 gusset plate design method may undersize gusset plates and that using a 2D BRBF model in this case can significantly underestimate λ꜀. Three improved alternative gusset plate design methods that are easy to implement in practice are identified from the literature. Gusset plates in two case-study 4-storey steel BRBF buildings with super-X and diagonal configurations are designed using both the NZS-3404 method and alternative methods. All three alternative design methods are found to be conservative, resulting in an almost three-fold lower λ꜀ for both case-study BRBF buildings compared to those designed using the NZS-3404 method. Analysis results indicate that (i) bidirectional interaction has no significant effect on gusset plate buckling and (ii) mid-span gusset plates are more susceptible to buckling than corner gusset plates. A framework for seismic loss assessment using incremental dynamic analysis (IDA), called loss-oriented hazard-consistent incremental dynamic analysis (LOHC-IDA), is developed. IDA can be conducted with a generic record set, eliminating the arduous site-specific record selection required to conduct multiple stripe analysis (MSA). Traditional IDA, however, is limited in producing hazard-consistent estimates of engineering demand parameters (EDPs), which LOHC-IDA overcomes. LOHC-IDA improves upon existing methods by: (i) incorporating correlations among engineering demand parameters across intensity levels and (ii) using peak ground acceleration (PGA) to predict peak floor acceleration (PFA). For two case-study steel BRBF buildings, LOHC-IDA estimates the EAL and loss distributions conditioned on the intensity level that closely match the MSA results, with an average absolute error of 5%. The influence of factors beyond gusset plate design on the λ꜀ and EAL of 26 case-study steel BRBF buildings (designed in accordance with TS 1170.5) is examined. Hazard-consistent λ꜀ and EAL for these buildings are estimated using the FEMA P-58 loss and risk assessment framework. Among the 26 case-study buildings, 23 satisfy the maximum code-specified λ꜀ limit of 10−⁴. The EAL, normalised by the total building replacement cost, is highest for 2-storey BRBFs (0.22% on average), followed by 4-storey BRBFs (0.16% on average) and 8-storey BRBFs (0.11% on average). Reducing the design drift limit has the most significant effect on lowering λ꜀ (all BRBF designs were drift governed), followed by transitioning from pinned to moment-resisting beam-column connections, reducing the brace angle, and increasing brace strength. BRBF buildings designed using the equivalent lateral force method, on average, have a lower λ꜀ compared to those designed using the modal response spectrum method. Diagonally configured BRBFs exhibit the lowest λ꜀, followed by super- X and chevron configured BRBFs. Most design variables, apart from drift limit and beam-column connection, have limited influence on EAL. A simple method for EDP-targeted design of steel BRBF buildings is proposed. For this purpose, linear regression and CatBoost machine learning models are developed to predict steel BRBF building EDPs using peak storey drift ratio (PSDR) and PFA estimates from the 26 case-study buildings at intensity levels ranging from 80% to 0.5% probability of exceedance in 50 years. The R²ₐₔⱼ of these models is around 0.98, while the average prediction error is less than 10%. Fundamental period (T₁), total building height (Hₜ), and pseudospectral acceleration at T₁, denoted as Sₐ(T₁), are selected as the features to predict PSDR, while T₁, Hₜ, and PGA are the features selected to predict PFA. The EDP-targeted design has three steps: (i) for a given Hₜ value, the PSDR prediction model is used to identify a suitable T₁ that can achieve a desired PSDR target at the design intensity, (ii) a force-based design is then conducted iteratively to achieve the target T₁ by using an appropriate ductility factor and design drift limit, and (iii) based on the T₁ in the final design iteration, the PFA demand estimated by the PFA prediction models is used as a conservative input for the design of acceleration-sensitive non-structural elements. An equation to predict λ꜀ at the design stage is proposed for collapse risk-targeted seismic design of buildings. This equation comprises three principal components: reserve building strength, a proxy for effective structural stiffness, and reserve building deformation capacity. This equation is calibrated for the collapse risk-targeted design of BRBF buildings in New Zealand using results from 26 case-study BRBF buildings. The validity of this equation is demonstrated with three design verification examples designed to specific λ꜀ targets. Considering λ꜀ from hazard-consistent incremental dynamic analysis as the benchmark, the mean absolute percentage error in the design-stage prediction of λ꜀ of the verification buildings is approximately 10%.

Research papers, University of Canterbury Library

Field surveys and experimental studies have shown that light steel or timber framed plasterboard partition walls are particularly vulnerable to earthquake damage prompting the overarching objective of this research, which is to further the development of low damage seismic systems for non-structural partition walls in order to facilitate their adoption by industry to assist with reducing the losses associated with the maintenance and repair cost of buildings across their design life. In particular, this study focused on the behaviour of steel-framed partition walls systems with novel detailing that aim to be “low-damage” designed according to common practice for walls used in commercial and institutional buildings in New Zealand. This objective was investigated by (1) investigating the performance of a flexible track system proposed by researchers and industry by experimental testing of full-scale specimens; (2) investigating the performance of the seismic gap partition wall systems proposed in a number of studies, further developed in this study with input from industry, by experimental testing of full-scale specimens; and (3) investigating the potential implications of using these systems compared with traditionally detailed partition wall systems within multi-storey buildings using the Performance Based Earthquake Engineering loss assessment methodology. Three full-scale testing frames were designed in order to replicate, under controlled laboratory conditions, the effects of seismic shaking on partition walls within multi-storey buildings by the application of quasi-static uni-directional cyclic loading imposing an inter-storey drift. The typical configuration for test specimens was selected to be a unique “y-shape”, including one angled return wall, with typical dimensions of approximately 2400 mm along the main wall and 600 mm along (approximately) the returns walls with a height of 2405 mm from floor to ceiling. The specimens were aligned within test frames at an oblique angle to the direction of loading in order to investigate bi- directional effects. Three wall specimens with flexible track detailing, two identical plane specimens and the third including a doorway, were tested. The detailing involved removing top track anchors within the proximity of wall intersections, thus allowing the tracks to ‘bow’ out at these locations. Although the top track anchors were specified to be removed the proximity of wall intersections, a construction error was made whereby a single top track slab to concrete anchor was left in at the three-way wall junction. Despite this error, the experimental testing was deemed worthwhile since such errors will also occur in practice and because the behaviour of the wall can be examined with this fixing in mind. The specimens also included an acoustic/fire sealant at the top lining to floor boundary. In addition to providing drift capacities, the force-displacement behaviour is also reported, the dissipated energy was computed, and the parameters of the Wayne-Stewart hysteretic model were fitted to the results. The specimen with the door opening behaved significantly different to the plane specimens: damage to the doorway specimen began as cracking of the wallboard propagating from the corners of the doorway following which the L- and Y- shaped junctions behaved independently, whereas damage to the plane specimens began as cracking of the wallboard at the top of the L-junction and wall system deformed as a single unit. The results suggest that bi-directional behaviour is important even if its impact cannot be directly quantified by the experiments conducted. Damage to sealant implies that the bond between plasterboard and sealant is important for its seismic performance. Careful quality control is advised as defects in the bond may significantly impact its ability to withstand seismic movement. Two specimens with seismic gap detailing were tested: a steel stud specimen and a timber stud specimen. Observed drift capacities were significantly greater than traditional plasterboard partition systems. Equations were used to predict the drift at which damage state 1 (DS1) and damage state 2 (DS2) would initiate. The equation used to estimate the drift at the onset of DS1 accurately predicted the onset of plaster cracking but overestimated the drift at which the gap filling material was damaged. The equation used to predict the onset of DS2 provided a lower bound for both specimens and also when used to predict results of previous experimental tests on seismic gap systems. The gap-filling material reduced the drift at the onset of DS1, however, it had a beneficial effect on the re-centring behaviour of the linings. Out-of-plane displacements and return wall configuration did not appear to significantly impact the onset of plaster cracking in the specimens. A loss assessment according to the PBEE methodology was conducted on four steel MRF case study buildings: (1) a 4-storey building designed for the Christchurch region, (2) a 4-storey building designed for the Wellington region, (3) a 12-storey building designed for the Christchurch region, and (4) a 12- storey building designed for the Wellington region. The fragility parameters for a traditional partition system, the flexible track partition system, and the seismic gap steel stud and timber stud partition systems were included within the loss assessment. The order (lowest to highest) of each system in terms of the expected annual losses of each building when incorporating the system was, (1) the seismic gap timber stud system, (2) the seismic gap steel stud system, (3) the traditional/baseline system, and (4) the flexible track system. For the seismic gap timber stud system, which incurred the greatest reduction in expected annual losses for each case study building, the reduction in expected annual losses in comparison to the losses found when using the traditional system ranged from a 5% to a 30% reduction. This reinforces the fact that while there is a benefit to the using low damage partition systems in each building the extent of reduction in expected annual losses is significantly dependent on the particular building design and its location. The flexible track specimens had larger repair costs at small hazard levels compared to the traditional system but smaller repair costs at larger hazard levels. However, the resulting expected annual losses for the flexible track system was higher than the traditional system which reinforces findings from past studies which observed that the greatest contribution to expected annual losses arises from low to moderate intensity shaking seismic events (low hazard levels).

Research papers, University of Canterbury Library

Between 2010 and 2011, Canterbury experienced a series of four large earthquake events with associated aftershocks which caused widespread damage to residential and commercial infrastructure. Fine grained and uncompacted alluvial soils, typical to the Canterbury outwash plains, were exposed to high peak ground acceleration (PGA) during these events. This rapid increase in PGA induced cyclic strain softening and liquefaction in the saturated, near surface alluvial soils. Extensive research into understanding the response of soils in Canterbury to dynamic loading has since occurred. The Earthquake Commission (EQC), the Ministry of Business and Employment (MBIE), and the Christchurch City Council (CCC) have quantified the potential hazards associated with future seismic events. Theses bodies have tested numerous ground improvement design methods, and subsequently are at the forefront of the Canterbury recovery and rebuild process. Deep Soil Mixing (DSM) has been proven as a viable ground improvement foundation method used to enhance in situ soils by increasing stiffness and positively altering in situ soil characteristics. However, current industry practice for confirming the effectiveness of the DSM method involves specific laboratory and absolute soil test methods associated with the mixed column element itself. Currently, the response of the soil around the columns to DSM installation is poorly understood. This research aims to understand and quantify the effects of DSM columns on near surface alluvial soils between the DSM columns though the implementation of standardised empirical soil test methods. These soil strength properties and ground improvement changes have been investigated using shear wave velocity (Vs), soil behaviour and density response methods. The results of the three different empirical tests indicated a consistent improvement within the ground around the DSM columns in sandier soils. By contrast, cohesive silty soils portrayed less of a consistent response to DSM, although still recorded increases. Generally, within the tests completed 50 mm from the column edge, the soil response indicated a deterioration to DSM. This is likely to be a result of the destruction of the soil fabric as the stress and strain of DSM is applied to the un‐mixed in situ soils. The results suggest that during the installation of DSM columns, a positive ground effect occurs in a similar way to other methods of ground improvement. However, further research, including additional testing following this empirical method, laboratory testing and finite 2D and 3D modelling, would be useful to quantify, in detail, how in situ soils respond and how practitioners should consider these test results in their designs. This thesis begins to evaluate how alluvial soils tend to respond to DSM. Conducting more testing on the research site, on other sites in Christchurch, and around the world, would provide a more complete data set to confirm the results of this research and enable further evaluation. Completing this additional research could help geotechnical DSM practitioners to use standardised empirical test methods to measure and confirm ground improvement rather than using existing test methods in future DSM projects. Further, demonstrating the effectiveness of empirical test methods in a DSM context is likely to enable more cost effective and efficient testing of DSM columns in future geotechnical projects.

Research papers, University of Canterbury Library

Probabilistic Structural Fire Engineering (PSFE) has been introduced to overcome the limitations of current conventional approaches used for the design of fire-exposed structures. Current structural fire design investigates worst-case fire scenarios and include multiple thermal and structural analyses. PSFE permits buildings to be designed to a level of life safety or economic loss that may occur in future fire events with the help of a probabilistic approach. This thesis presents modifications to the adoption of a Performance-Based Earthquake Engineering (PBEE) framework in Probabilistic Structural Fire Engineering (PSFE). The probabilistic approach runs through a series of interrelationships between different variables, and successive convolution integrals of these interrelationships result in probabilities of different measures. The process starts with the definition of a fire severity measure (FSM), which best relates fire hazard intensity with structural response. It is identified by satisfying efficiency and sufficiency criteria as described by the PBEE framework. The relationship between a fire hazard and corresponding structural response is established by analysis methods. One method that has been used to quantify this relationship in PSFE is Incremental Fire Analysis (IFA). The existing IFA approach produces unrealistic fire scenarios, as fire profiles may be scaled to wide ranges of fire severity levels, which may not physically represent any real fires. Two new techniques are introduced in this thesis to limit extensive scaling. In order to obtain an annual rate of exceedance of fire hazard and structural response for an office building, an occurrence model and an attenuation model for office fires are generated for both Christchurch city and New Zealand. The results show that Christchurch city is 15% less likely to experience fires that have the potential to cause structural failures in comparison to all of New Zealand. In establishing better predictive relationships between fires and structural response, cumulative incident radiation (a fire hazard property) is found to be the most appropriate fire severity measure. This research brings together existing research on various sources of uncertainty in probabilistic structural fire engineering, such as elements affecting post-flashover fire development factors (fuel load, ventilation, surface lining and compartment geometry), fire models, analysis methods and structural reliability. Epistemic uncertainty and aleatory uncertainty are investigated in the thesis by examining the uncertainty associated with modelling and the factors that influence post-flashover development of fires. A survey of 12 buildings in Christchurch in combination with recent surveys in New Zealand produced new statistical data on post-flashover development factors in office buildings in New Zealand. The effects of these parameters on temperature-time profiles are evaluated. The effects of epistemic uncertainty due to fire models in the estimation of structural response is also calculated. Parametric fires are found to have large uncertainty in the prediction of post-flashover fires, while the BFD curves have large uncertainties in prediction of structural response. These uncertainties need to be incorporated into failure probability calculations. Uncertainty in structural modelling shows that the choices that are made during modelling have a large influence on realistic predictions of structural response.

Research papers, University of Canterbury Library

The lateral capacity of a conventional CLT shear wall is often governed by the strength and stiffness of its connections, which do not significantly utilize the in-plane strength of the CLT. Therefore, CLT shear walls are not yet being used efficiently in the construction of mass timber buildings due to a lack of research on high-capacity connections and alternative wall configurations. In this study, cyclic experiments were completed on six full-scale, 5-ply cantilever CLT shear walls with high-capacity hold-downs using mixed angle screws and bolts. All specimens exhibited significantly higher strength and stiffness than previously tested conventional CLT shear walls in the literature. The base connections demonstrated ductile failure modes through yielding of the hold-down connections. Based on the experimental results, numerical models were calibrated to investigate the seismic behaviour of CLT shear walls for prototype buildings of 3 and 6-storeys in Christchurch, NZ. As an alternative to cantilever (single) shear walls, a type of coupled wall with steel link beams between adjacent CLT wall piers was investigated. Effective coupling requires the link beam-to-wall connections to have adequate strength to ensure ductile link beam responses and adequate stiffness to yield the link beams at a relatively low inter-storey drift level. To this end, three beam-to-wall connection types were developed and cyclically tested to investigate their behaviour and feasibility. Based on the test results of the critical connection, a 3-storey, 2/3-scale coupled CLT wall specimen with three steel link beams and mixed angle screwed hold-downs was cyclically tested to evaluate its performance and experimentally validate the system concept. The test results showed a relatively high lateral strength compared to conventional CLT shear walls, as well as a high system ductility ratio of 7.6. Failure of the system was characterised by combined bending and withdrawal of the screws in the mixed angle screw hold-downs, yielding and eventual inelastic buckling of the steel link beams, CLT toe crushing, and local CLT delamination. Following the initial test, the steel link beams, mixed angle screw hold-downs, and damaged CLT regions were repaired, then the wall specimen was re-tested. The repaired wall behaved similarly to the original test and exhibited slightly higher energy dissipation and peak strength, but marginally more rapid strength deterioration under cyclic loading. Several hybrid coupled CLT shear walls were numerically modelled and calibrated based on the results of the coupled wall experiments. Pushover analyses were conducted on a series of configurations to validate a capacity design method for the system and to investigate reasonable parameter values for use in the preliminary design of the system. Additionally, an iterative seismic design method was proposed and used to design sample buildings of 6, 8, and 10-storeys using both nonlinear pushover and nonlinear time history analyses to verify the prototype designs. Results of the sample building analyses demonstrated adequate seismic behaviour and the proposed design parameters were found to be appropriate. In summary, high-capacity CLT shear walls can be used for the resistance of earthquakes by using stronger base connections and coupled wall configurations. The large-scale experimental testing in this study has demonstrated that both cantilever and coupled CLT shear walls are feasible LLRSs which can provide significantly greater lateral strength, stiffness, and energy dissipation than conventional CLT shear wall configurations.

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, Lincoln University

Mitigating the cascade of environmental damage caused by the movement of excess reactive nitrogen (N) from land to sea is currently limited by difficulties in precisely and accurately measuring N fluxes due to variable rates of attenuation (denitrification) during transport. This thesis develops the use of the natural abundance isotopic composition of nitrate (δ15N and δ18O of NO₃-) to integrate the spatialtemporal variability inherent to denitrification, creating an empirical framework for evaluating attenuation during land to water NO₃- transfers. This technique is based on the knowledge that denitrifiers kinetically discriminate against 'heavy' forms of both N and oxygen (O), creating a parallel enrichment in isotopes of both species as the reaction progresses. This discrimination can be quantitatively related to NO₃- attenuation by isotopic enrichment factors (εdenit). However, while these principles are understood, use of NO₃- isotopes to quantify denitrification fluxes in non-marine environments has been limited by, 1) poor understanding of εdenit variability, and, 2) difficulty in distinguishing the extent of mixing of isotopically distinct sources from the imprint of denitrification. Through a combination of critical literature analysis, mathematical modelling, mesocosm to field scale experiments, and empirical studies on two river systems over distance and time, these short comings are parametrised and a template for future NO₃- isotope based attenuation measurements outlined. Published εdenit values (n = 169) are collated in the literature analysis presented in Chapter 2. By evaluating these values in the context of known controllers on the denitrification process, it is found that the magnitude of εdenit, for both δ15N and δ18O, is controlled by, 1) biology, 2) mode of transport through the denitrifying zone (diffusion v. advection), and, 3) nitrification (spatial-temporal distance between nitrification and denitrification). Based on the outcomes of this synthesis, the impact of the three factors identified as controlling εdenit are quantified in the context of freshwater systems by combining simple mathematical modelling and lab incubation studies (comparison of natural variation in biological versus physical expression). Biologically-defined εdenit, measured in sediments collected from four sites along a temperate stream and from three tropical submerged paddy fields, varied from -3‰ to -28‰ depending on the site’s antecedent carbon content. Following diffusive transport to aerobic surface water, εdenit was found to become more homogeneous, but also lower, with the strength of the effect controlled primarily by diffusive distance and the rate of denitrification in the sediments. I conclude that, given the variability in fractionation dynamics at all levels, applying a range of εdenit from -2‰ to -10‰ provides more accurate measurements of attenuation than attempting to establish a site-specific value. Applying this understanding of denitrification's fractionation dynamics, four field studies were conducted to measure denitrification/ NO₃- attenuation across diverse terrestrial → freshwater systems. The development of NO₃- isotopic signatures (i.e., the impact of nitrification, biological N fixation, and ammonia volatilisation on the isotopic 'imprint' of denitrification) were evaluated within two key agricultural regions: New Zealand grazed pastures (Chapter 4) and Philippine lowland submerged rice production (Chapter 5). By measuring the isotopic composition of soil ammonium, NO₃- and volatilised ammonia following the bovine urine deposition, it was determined that the isotopic composition of NO₃ - leached from grazed pastures is defined by the balance between nitrification and denitrification, not ammonia volatilisation. Consequently, NO₃- created within pasture systems was predicted to range from +10‰ (δ15N)and -0.9‰ (δ18O) for non-fertilised fields (N limited) to -3‰ (δ15N) and +2‰ (δ18O) for grazed fertilised fields (N saturated). Denitrification was also the dominant determinant of NO₃- signatures in the Philippine rice paddy. Using a site-specific εdenit for the paddy, N inputs versus attenuation were able to be calculated, revealing that >50% of available N in the top 10 cm of soil was denitrified during land preparation, and >80% of available N by two weeks post-transplanting. Intriguingly, this denitrification was driven by rapid NO₃- production via nitrification of newly mineralised N during land preparation activities. Building on the relevant range of εdenit established in Chapters 2 and 3, as well as the soil-zone confirmation that denitrification was the primary determinant of NO₃- isotopic composition, two long-term longitudinal river studies were conducted to assess attenuation during transport. In Chapter 6, impact and recovery dynamics in an urban stream were assessed over six months along a longitudinal impact gradient using measurements of NO₃- dual isotopes, biological populations, and stream chemistry. Within 10 days of the catastrophic Christchurch earthquake, dissolved oxygen in the lowest reaches was <1 mg l⁻¹, in-stream denitrification accelerated (attenuating 40-80% of sewage N), microbial biofilm communities changed, and several benthic invertebrate taxa disappeared. To test the strength of this method for tackling the diffuse, chronic N loading of streams in agricultural regions, two years of longitudinal measurements of NO₃- isotopes were collected. Attenuation was negatively correlated with NO₃- concentration, and was highly dependent on rainfall: 93% of calculated attenuation (20 kg NO₃--N ha⁻¹ y⁻¹) occurred within 48 h of rainfall. The results of these studies demonstrate the power of intense measurements of NO₃- stable isotope for distinguishing temporal and spatial trends in NO₃ - loss pathways, and potentially allow for improved catchment-scale management of agricultural intensification. Overall this work now provides a more cohesive understanding for expanding the use of NO₃- isotopes measurements to generate accurate understandings of the controls on N losses. This information is becoming increasingly important to predict ecosystem response to future changes, such the increasing agricultural intensity needed to meet global food demand, which is occurring synergistically with unpredictable global climate change.

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 nonlinear dynamic soil-foundation-structure interaction (SFSI) can signifi cantly affect the seismic response of buildings, causing additional deformation modes, damage and repair costs. Because of nonlinear foundation behaviour and interactions, the seismic demand on the superstructure may considerably change, and also permanent deformations at the foundation level may occur. Although SFSI effects may be benefi cial to the superstructure performance, any advantage would be of little structural value unless the phenomenon can be reliably controlled and exploited. Detrimental SFSI effects may also occur, including acceleration and displacement response ampli cation and differential settlements, which would be unconservative to neglect. The lack of proper understanding of the phenomenon and the limited available simpli ed tools accounting for SFSI have been major obstacles to the implementation of integrated design and assessment procedures into the everyday practice. In this study concepts, ideas and practical tools (inelastic spectra) for the seismic design and assessment of integrated foundation-superstructure systems are presented, with the aim to explicitly consider the impact of nonlinearities occurring at the soil-foundation interface on the building response within an integrated approach, where the foundation soil and superstructure are considered as part of an integrated system when evaluating the seismic response, working synergically for the achievement of a target global performance. A conceptual performance-based framework for the seismic design and assessment of integrated foundation-superstructure systems is developed. The framework is based on the use of peak and residual response parameters for both the superstructure and the foundation, which are then combined to produce the system performance matrix. Each performance matrix allows for worsening of the performance when different contributions are combined. An attempt is made to test the framework by using case histories from the 2011 Christchurch earthquake, which are previously shown to have been severely affected by nonlinear SFSI. The application highlights the framework sensitivity to the adopted performance limit states, which must be realistic for a reliable evaluation of the system performance. Constant ductility and constant strength inelastic spectra are generated for nonlinear SFSI systems (SDOF nonlinear superstructure and 3DOF foundation allowing for uplift and soil yielding), representing multistorey RC buildings with shallow rigid foundations supported by cohesive soils. Different ductilities/strengths, hysteretic rules (Bi-linear, Takeda and Flag-Shape), soil stiffness and strength and bearing capacity factors are considered. Footings and raft foundations are investigated, characterized respectively by constant (3 and 8) and typically large bearing capacity factors. It is confi rmed that when SFSI is considered, the superstructure yielding force needed to satisfy a target ductility for a new building changes, and that similarly, for an existing building, the ductility demand on a building of a given strength varies. The extent of change of seismic response with respect to xed-base (FB) conditions depends on the class of soils considered, and on the bearing capacity factor (SF). For SF equal to 3, the stiffer soils enhance the nonlinear rotational foundation behaviour and are associated with reduced settlement, while the softer ones are associated with increased settlement response but not signi ficant rotational behaviour. On average terms, for the simplifi ed models considered, SFSI is found to be bene cial to the superstructure performance in terms of acceleration and superstructure displacement demand, although exceptions are recorded due to ground motion variability. Conversely, in terms of total displacement, a signi cant response increase is observed. The larger the bearing capacity factor, the more the SFSI response approaches the FB system. For raft foundation buildings, characterized by large bearing capacity factors, the impact of foundation response is mostly elastic, and the system on average approaches FB conditions. Well de fined displacement participation factors to the peak total lateral displacement are observed for the different contributions (i.e. peak foundation rotation and translation and superstructure displacement). While the superstructure and foundation rotation show compensating trends, the foundation translation contribution varies as a function of the moment-to-shear ratio, becoming negligible in the medium-to-long periods. The longer the superstructure FB period, the less the foundation response is signifi cant. The larger the excitation level and the less ductile the superstructure, the larger the foundation contribution to the total lateral displacement, and the less the superstructure contribution. In terms of hysteretic behaviour, its impact is larger when the superstructure response is more signifi cant, i.e. for the softer/weaker soils and larger ductilities. Particularly, for the Flag Shape rule, larger superstructure displacement participation factors and smaller foundation contributions are recorded. In terms of residual displacements, the total residual-to-maximum ratios are similar in amplitudes and trends to the corresponding FB system responses, with the foundation and superstructure contributions showing complementary trends. The impact of nonlinear SFSI is especially important for the Flag Shape hysteresis rule, which would not otherwise suffer of any permanent deformations. By using the generated peak and residual inelastic spectra (i.e. inelastic acceleration/ displacement modifi cation factor spectra, and/or participation factor and residual spectra), conceptual simplifi ed procedures for the seismic design and assessment of integrated foundation-superstructure systems are presented. The residual displacements at both the superstructure and foundation levels are explicitly considered. Both the force- and displacement-based approaches are explored. The procedures are de fined to be complementary to the previously proposed integrated performance-based framework. The use of participation factor spectra allows the designer to easily visualize the response of the system components, and could assist the decision making process of both the design and assessment of SFSI systems. The presented numerical results have been obtained using simpli ed models, assuming rigid foundation behaviour and neglecting P-Delta effects. The consideration of more complex systems including asymmetry in stiffness, mass, axial load and ground conditions with a exible foundation layout would highlight detrimental SFSI effects as related to induced differential settlements, while accounting for PDelta effects would further amplify the displacement response. Also, the adopted acceleration records were selected and scaled to match conventional design spectra, thus not representing any response ampli cation in the medium-to-long period range which could as well cause detrimental SFSI effects. While these limitations should be the subject of further research, this study makes a step forward to the understanding of SFSI phenomenon and its incorporation into performance-based design/assessment considerations.

Research papers, University of Canterbury Library

Environmental stress and disturbance can affect the structure and functioning of marine ecosystems by altering their physical, chemical and biological features. In estuaries, benthic invertebrate communities play important roles in structuring sediments, influencing primary production and biogeochemical flux, and occupying key food web positions. Stress and disturbance can reduce species diversity, richness and abundance, with ecological theory predicting that biodiversity will be at its lowest soon after a disturbance with assemblages dominated by opportunistic species. The Avon-Heathcote Estuary in Christchurch New Zealand has provided a novel opportunity to examine the effects of stress, in the form of eutrophication, and disturbance, in the form of cataclysmic earthquake events, on the structure and functioning of an estuarine ecosystem. For more than 50 years, large quantities (up to 500,000m3/day) of treated wastewater were released into this estuary but in March 2010 this was diverted to an ocean outfall, thereby reducing the nutrient loading by around 90% to the estuary. This study was therefore initially focussed on the reversal of eutrophication and consequent effects on food web structure in the estuary as it responded to lower nutrients. In 2011, however, Christchurch was struck with a series of large earthquakes that greatly changed the estuary. Massive amounts of liquefied sediments, covering up to 65% of the estuary floor, were forced up from deep below the estuary, the estuary was tilted by up to a 50cm rise on one side and a corresponding drop on the other, and large quantities of raw sewage from broken wastewater infrastructure entered the estuary for up to nine months. This study was therefore a test of the potentially synergistic effects of nutrient reduction and earthquake disturbance on invertebrate communities, associated habitats and food web dynamics. Because there was considerable site-to-site heterogeneity in the estuary, the sites in this study were selected to represent a eutrophication gradient from relatively “clean” (where the influence of tidal flows was high) to highly impacted (near the historical discharge site). The study was structured around these sites, with components before the wastewater diversion, after the diversion but before the earthquakes, and after the earthquakes. The eutrophication gradient was reflected in the composition and isotopic chemistry of primary producer and invertebrate communities and the characteristics of sediments across the sample sites. Sites closest to the former wastewater discharge pipe were the most eutrophic and had cohesive organic -rich, fine sediments and relatively depauperate communities dominated by the opportunistic taxa Capitellidae. The less-impacted sites had coarser, sandier sediments with fewer pollutants and far less organic matter than at the eutrophic sites, relatively high diversity and lower abundances of micro- and macro-algae. Sewage-derived nitrogen had became incorporated into the estuarine food web at the eutrophic sites, starting at the base of the food chain with benthic microalgae (BMA), which were found to use mostly sediment-derived nitrogen. Stable isotopic analysis showed that δ13C and δ15N values of most food sources and consumers varied spatially, temporally and in relation to the diversion of wastewater, whereas the earthquakes did not appear to affect the overall estuarine food web structure. This was seen particularly at the most eutrophic site, where isotopic signatures became more similar to the cleaner sites over two-and-a-half years after the diversion. New sediments (liquefaction) produced by the earthquakes were found to be coarser, have lower concentrations of heavy metals and less organic matter than old (existing) sediments. They also had fewer macroinvertebrate inhabitants initially after the earthquakes but most areas recovered to pre-earthquake abundance and diversity within two years. Field experiments showed that there were higher amounts of primary production and lower amounts of nutrient efflux from new sediments at the eutrophic sites after the earthquakes. Primary production was highest in new sediments due to the increased photosynthetic efficiency of BMA resulting from the increased permeability of new sediments allowing increased light penetration, enhanced vertical migration of BMA and the enhanced transport of oxygen and nutrients. The reduced efflux of NH4-N in new sediments indicated that the capping of a large portion of eutrophic old sediments with new sediments had reduced the release of legacy nutrients (originating from the historical discharge) from the sediments to the overlying water. Laboratory experiments using an array of species and old and new sediments showed that invertebrates altered levels of primary production and nutrient flux but effects varied among species. The mud snail Amphibola crenata and mud crab Austrohelice crassa were found to reduce primary production and BMA biomass through the consumption of BMA (both species) and its burial from bioturbation and the construction of burrows (Austrohelice). In contrast, the cockle Austrovenus stutchburyi did not significantly affect primary production and BMA biomass. These results show that changes in the structure of invertebrate communities resulting from disturbances can also have consequences for the functioning of the system. The major conclusions of this study were that the wastewater diversion had a major effect on food web dynamics and that the large quantities of clean and unpolluted new sediments introduced to the estuary during the earthquakes altered the recovery trajectory of the estuary, accelerating it at least throughout the duration of this study. This was largely through the ‘capping’ effect of the new liquefied, coarser-grained sediments as they dissipated across the estuary and covered much of the old organic-rich eutrophic sediments. For all aspects of this study, the largest changes occurred at the most eutrophic sites; however, the surrounding habitats were important as they provided the context for recovery of the estuary, particularly because of the very strong influence of sediments, their biogeochemistry, microalgal and macroalgal dynamics. There have been few studies documenting system level responses to eutrophication amelioration and to the best on my knowledge there are no other published studies examining the impacts of large earthquakes on benthic communities in an estuarine ecosystem. This research gives valuable insight and advancements in the scientific understanding of the effects that eutrophication recovery and large-scale disturbances can have on the ecology of a soft-sediment ecosystem.