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.
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.
A video of a tractor on a farm near River Road in Lincoln. The tractor is making a several passes over a sand volcano with a rotary hoe. This was one of several soil-remediation techniques tested on farms affected by the 4 September 2010 earthquake.
A video of a tractor on a farm near River Road in Lincoln. The tractor is passing over a large sand volcano with a power harrow. This was one of several soil-remediation techniques tested on farms affected by the 4 September 2010 earthquake.
A video of a dug-out soil pit on a farm near River Road in Lincoln. There is a thick layer of saturated sand between the topsoil and the subsoil where the grass roots end. The saturated soil is collapsing as a result of contractor Tony Fisher, of Fisher Agriculture Ltd, jumping lightly up and down on the edge of the pit.
A video of a dug-out soil pit on a farm near River Road in Lincoln. There is a thick layer of saturated sand between the topsoil and the subsoil where the grass roots end. The saturated soil is collapsing as a result of contractor Tony Fisher, of Fisher Agriculture Ltd, having jumped lightly up and down on the edge of the pit.
A video of contractor Tony Fisher, of Fisher Agricultural Ltd, jumping lightly up and down on the edge of a dug-out soil pit on a farm near River Road in Lincoln. There is a thick layer of saturated sand between the topsoil and the subsoil where the grass roots end. The saturated soil starts to collapse as a result of the movement.
A non-destructive hardness testing method has been developed to investigate the amount of plastic strain demand in steel elements subjected to cyclic loading. The focus of this research is on application to the active links of eccentrically braced frames (EBFs), which are a commonly used seismic-resisting system in modern steel framed buildings. The 2010/2011 Christchurch earthquake series, especially the very intense February 22 shaking, which was the first earthquake worldwide to push complete EBF systems fully into their inelastic state, generating a moderate to high level of plastic strain in EBF active links, for a range of buildings from 3 to 23 storeys in height. This raised two important questions: 1) what was the extent of plastic deformation in active links; and 2) what effect does that have to post-earthquake steel properties? This project comprised determining a robust relationship between hardness and plastic strain in order to be able to answer the first question and provide the necessary input into answering the second question. A non-destructive Leeb (portable) hardness tester (model TH170) has been used to measure the hardness, in order to determine the plastic strain, in hot rolled steel universal sections and steel plates. A bench top Rockwell B was used to compare and validated the hardness measured by the portable hardness tester. Hardness was measured from monotonically strained tensile test specimens to identify the relationship between hardness and plastic strain demand. Test results confirmed a good relationship between hardness and the amount of monotonically induced plastic strain. Surface roughness was identified as an important parameter in obtaining reliable hardness readings from a portable hardness reader. A proper surface preparation method was established by using three different cleaning methods, finished with hand sanding to achieve surface roughness coefficients sufficiently low not to distort the results. This work showed that a test surface roughness (Ra) is not more than 1.6 micron meter (μm) is required for accurate readings from the TH170 tester. A case study on an earthquake affected building was carried out to identify the relationship between hardness and amount of plastic strain demand in cyclically deformed active links. Hardness was carried out from active links shown visually to have been the most affected during one of the major earthquake events. Onsite hardness test results were then compared with laboratory hardness test results. A good relationship between hardness from onsite and laboratory was observed between the test methods; Rockwell B bench top and portable Leeb tester TH170. Manufacturing induced plastic strain in the top and bottom of the webs of hot rolled sections were discovered from this research, an important result which explains why visual effects of earthquake induced active link yielding (eg cracked or flaking paint) was typically more prevalent over the middle half depth of the active link. The extent of this was quantified. It was also evident that the hardness readings from the portable hardness tester are influenced by geometry, mass effects and rigidity of the links. The final experimental stage was application of the method to full scale cyclic inelastic tested nominally identical active links subjected to loading regimes comprising constant and variable plastic strain demands. The links were cyclically loaded to achieve different plastic strain level. A novel Digital Image Correlation (DIC) technique was incorporated during the tests of this scale, to confirm the level of plastic strain achieved. Tensile test specimens were water jet cut from cyclically deformed webs to analyse the level of plastic strain. Test results show clear evidence that cyclically deformed structural steel elements show good correlation between hardness and the amount of plastic strain demand. DIC method was found to be reliable and accurate to check the level of plastic strain within cyclically deformed structural steel elements.
The Avon-Heathcote Estuary, located in Christchurch, New Zealand, experienced coseismic deformation as a result of the February 22nd 2011 Christchurch Earthquake. The deformation is reflected as subsidence in the northern area and uplift in the southern area of the Estuary, in addition to sand volcanoes which forced up sediment throughout the floor of the Estuary altering estuary bed height and tidal flow. The first part of the research involved quantifying the change in the modern benthic foraminifera distribution as a result of the coseismic deformation caused by the February 22nd 2011 earthquake. By analysing the taxa present immediately post deformation and then the taxa present 2 years post deformation a comparison of the benthic foraminifera distribution can be made of the pre and post deformation. Both the northern and the southern areas of the Estuary were sampled to establish whether foraminifera faunas migrated landward or seaward as a result of subsidence and uplift experienced in different areas. There was no statistical change in overall species distribution in the two year time period since the coseismic deformation occurred, however, there were some noticeable changes in foraminifera distribution at BSNS-Z3 showing a landward migration of taxa. The changes that were predicted to occur as a result of the deformation of the Estuary are taking longer than expected to show up in the foraminiferal record and a longer time period is needed to establish these changes. The second stage involved establishing the modern distribution of foraminifera at Settlers Reserve in the southern area of the Avon-Heathcote Estuary by detailed sampling along a 160 m transect. Foraminifera are sensitive to environmental parameters, tidal height, grainsize, pH and salinity were recorded to evaluate the effect these parameters have on distribution. Bray-Curtis two-way cluster analysis was primarily used to assess the distribution pattern of foraminifera. The modern foraminifera distribution is comparable to that of the modern day New Zealand brackish-water benthic foraminifera distribution and includes species not yet found in other studies of the Avon-Heathcote Estuary. Differences in sampling techniques and the restricted intertidal marshland area where the transect samples were collected account for some of the differences seen between this model and past foraminifera studies. xiii The final stage involved sampling a 2.20 m core collected from Settlers Reserve and using the modern foraminiferal distribution to establish a foraminiferal history of Settlers Reserve. As foraminifera are sensitive to tidal height they may record past coseismic deformation events and the core was used to ascertain whether record of past coseismic deformation is preserved in Settlers Reserve sediments. Sampling the core for foraminifera, grainsize, trace metals and carbon material helped to build a story of estuary development. Using the modern foraminiferal distribution and the tidal height information collected, a down core model of past tidal heights was established to determine past rates of change. Foraminifera are not well preserved throughout the core, however, a sudden relative rise in sea level is recorded between 0.25 m and 0.85 m. Using trace metal and isotope analysis to develop an age profile, this sea level rise is interpreted to record coseismic subsidence associated with a palaeoseismic event in the early 1900’s. Overall, although the Avon-Heathcote Estuary experienced clear coseismic deformation as a result of the 22nd of February 2011 earthquake, modern changes in foraminiferal distribution cannot yet be tracked, however, past seismic deformation is identified in a core. The modern transect describes the foraminifera distribution which identifies species that have not been identified in the Avon-Heathcote Estuary before. This thesis enhances the current knowledge of the Avon-Heathcote Estuary and is a baseline for future studies.
