This dissertation addresses a diverse range of topics in the area of physics-based ground motion simulation with particular focus on the Canterbury, New Zealand region. The objectives achieved provide the means to perform hybrid broadband ground motion simulation and subsequently validates the simulation methodology employed. In particu- lar, the following topics are addressed: the development of a 3D seismic velocity model of the Canterbury region for broadband ground motion simulation; the development of a 3D geologic model of the interbedded Quaternary formations to provide insight on observed ground motions; and the investigation of systematic effects through ground motion sim- ulation of small-to-moderate magnitude earthquakes. The paragraphs below outline each contribution in more detail. As a means to perform hybrid broadband ground motion simulation, a 3D model of the geologic structure and associated seismic velocities in the Canterbury region is devel- oped utilising data from depth-converted seismic reflection lines, petroleum and water well logs, cone penetration tests, and implicitly guided by existing contour maps and geologic cross sections in data sparse subregions. The model explicitly characterises five significant and regionally recognisable geologic surfaces that mark the boundaries between geologic units with distinct lithology and age, including the Banks Peninsula volcanics, which are noted to strongly influence seismic wave propagation. The Basement surface represents the base of the Canterbury sedimentary basin, where a large impedance contrast exists re- sulting in basin-generated waves. Seismic velocities for the lithological units between the geologic surfaces are derived from well logs, seismic reflection surveys, root mean square stacking velocities, empirical correlations, and benchmarked against a regional crustal model, thus providing the necessary information for a Canterbury velocity model for use in broadband seismic wave propagation. A 3D high-resolution model of the Quaternary geologic stratigraphic sequence in the Canterbury region is also developed utilising datasets of 527 high-quality water well logs, and 377 near-surface cone penetration test records. The model, developed using geostatistical Kriging, represents the complex interbedded regional Quaternary geology by characterising the boundaries between significant interbedded geologic formations as 3D surfaces including explicit modelling of the formation unconformities resulting from the Banks Peninsula volcanics. The stratigraphic layering present can result in complex wave propagation. The most prevalent trend observed in the surfaces was the downward dip from inland to the eastern coastline as a result of the dominant fluvial depositional environment of the terrestrial gravel formations. The developed model provides a benefi- cial contribution towards developing a comprehensive understanding of recorded ground motions in the region and also providing the necessary information for future site char- acterisation and site response analyses. To highlight the practicality of the model, an example illustrating the role of the model in constraining surface wave analysis-based shear wave velocity profiling is illustrated along with the calculation of transfer functions to quantify the effect of the interbedded geology on wave propagation. Lastly, an investigation of systematic biases in the (Graves and Pitarka, 2010, 2015) ground motion simulation methodology and the specific inputs used for the Canterbury region is presented considering 144 small-to-moderate magnitude earthquakes. In the simulation of these earthquakes, the 3D Canterbury Velocity Model, developed as a part of this dissertation, is used for the low-frequency simulation, and a regional 1D velocity model for the high-frequency simulation. Representative results for individual earthquake sources are first presented to highlight the characteristics of the small-to-moderate mag- nitude earthquake simulations through waveforms, intensity measure scaling with source- to-site distance, and spectral bias of the individual events. Subsequently, a residual de- composition is performed to examine the between- and within-event residuals between observed data, and simulated and empirical predictions. By decomposing the residuals into between- and within-event residuals, the biases in source, path and site effects, and their causes, can be inferred. The residuals are comprehensively examined considering their aggregated characteristics, dependence on predictor variables, spatial distribution, and site-specific effects. The results of the simulation are also benchmarked against empir- ical ground motion models, where their similarities manifest from common components in their prediction. Ultimately, suggestions to improve the predictive capability of the simulations are presented as a result of the analysis.
A review of the literature showed the lack of a truly effective damage avoidance solution for timber or hybrid timber moment resisting frames (MRFs). Full system damage avoidance selfcentring behaviour is difficult to achieve with existing systems due to damage to the floor slab caused by beam-elongation. A novel gravity rocking, self-centring beam-column joint with inherent and supplemental friction energy dissipation is proposed for low-medium rise buildings in all seismic zones where earthquake actions are greater than wind. Steel columns and timber beams are used in the hybrid MRF such that both the beam and column are continuous thus avoiding beam-elongation altogether. Corbels on the columns support the beams and generate resistance and self-centring through rocking under the influence of gravity. Supplemental friction sliders at the top of the beams resist sliding of the floor whilst dissipating energy as the floor lifts on the corbels and returns. 1:20 scale tests of 3-storey one-by-two bay building based on an earlier iteration of the proposed concept served as proof-of-concept and highlighted areas for improvement. A 1:5 scale 3-storey one-by-one bay building was subsequently designed. Sub-assembly tests of the beam-top asymmetric friction sliders demonstrated repeatable hysteresis. Quasi-static tests of the full building demonstrated a ‘flat bottomed’ flag-shaped hysteresis. Shake table tests to a suite of seven earthquakes scaled for Wellington with site soil type D to the serviceability limit state (SLS), ultimate limit state (ULS) and maximum credible event (MCE) intensity corresponding to an average return period of 25, 500 and 2500 years respectively were conducted. Additional earthquake records from the 22 February 2011 Christchurch earthquakes we included. A peak drift of 0.6%, 2.5% and 3.8% was reached for the worst SLS, ULS and MCE earthquake respectively whereas a peak drift of 4.5% was reached for the worst Christchurch record for tests in the plane of the MRF. Bi-directional tests were also conducted with the building oriented at 45 degrees on the shake table and the excitation factored by 1.41 to maintain the component in the direction of the MRF. Shear walls with friction slider hold-downs which reached similar drifts to the MRF were provided in the orthogonal direction. Similar peak drifts were reached by the MRF in the bi-directional tests, when the excitation was amplified as intended. The building self-centred with a maximum residual drift of 0.06% in the dynamic tests and demonstrated no significant damage. The member actions were magnified by up to 100% due to impact upon return of the floor after uplift when the peak drift reached 4.5%. Nonetheless, all of the members and connections remained essentially linearelastic. The shake table was able to produce a limited peak velocity of 0.275 m/s and this limited the severity of several of the ULS, MCE and Christchurch earthquakes, especially the near-field records with a large velocity pulse. The full earthquakes with uncapped velocity were simulated in a numerical model developed in SAP2000. The corbel supports were modelled with the friction isolator link element and the top sliders were modelled with a multi-linear plastic link element in parallel with a friction spring damper. The friction spring damper simulated the increase in resistance with increasing joint rotation and a near zero return stiffness, as exhibited by the 1:5 scale test building. A good match was achieved between the test quasi-static global force-displacement response and the numerical model, except a less flat unloading curve in the numerical model. The peak drift from the shake table tests also matched well. Simulations were also run for the full velocity earthquakes, including vertical ground acceleration and different floor imposed load scenarios. Excessive drift was predicted by the numerical model for the full velocity near-field earthquakes at the MCE intensity and a rubber stiffener for increasing the post joint-opening stiffness was found to limit the drift to 4.8%. Vertical ground acceleration had little effect on the global response. The system generates most of its lateral resistance from the floor weight, therefore increasing the floor imposed load increased the peak drift, but less than it would if the resistance of the system did not increase due to the additional floor load. A seismic design procedure was discussed under the framework of the existing direct displacement-based design method. An expression for calculating the area-based equivalent viscous damping (EVD) was derived and a conservative correction factor of 0.8 was suggested. A high EVD of up to about 15% can be achieved with the proposed system at high displacement ductility levels if the resistance of the top friction sliders is maximised without compromising reliable return of the floor after uplift. Uniform strength joints with an equal corbel length up the height of the building and similar inter-storey drifts result in minimal relative inter-floor uplift, except between the first floor and ground. Guidelines for detailing the joint for damage avoidance including bi-directional movement were also developed.
The Acheron rock avalanche is located in the Red Hill valley almost 80 km west of Christchurch and is one of 42 greywacke-derived rock avalanches identified in the central Southern Alps. It overlies the Holocene active Porters Pass Fault; a component of the Porters Pass-Amberley Fault Zone which extends from the Rakaia River to beyond the Waimakariri River. The Porters Pass Fault is a dextral strike-slip fault system viewed as a series of discontinuous fault scarps. The location of the fault trace beneath the deposit suggests it may represent a possible source of seismic shaking resulting in the formation of the Acheron rock avalanche. The rock mass composition of the rock avalanche source scar is Torlesse Supergroup greywacke consisting of massive sandstone and thinly bedded mudstone sequences dipping steeply north into the centre of the source basin. A stability analysis identified potential instability along shallow north dipping planar defects, and steep south dipping toppling failure planes. The interaction of the defects with bedding is considered to have formed conditions for potential instability most likely triggered by a seismic event. The dTositional area of the rock avalanche covers 7.2 x 105 m2 with an estimated volume of 9 x 10 m3 The mobilised rock mass volume was calculated at 7.5 x 106 m3• Run out of the debris from the top of the source scar to the distal limit reached 3500m, descending over a vertical fall of almost 700m with an estimated Fahrboschung of 0.2. The run out of the rock avalanche displayed moderate to high mobility, travelling at an estimated maximum velocity of 140-160 km/hour. The rapid emplacement of the deposit is confirmed by highly fragmented internal composition and burial of forest vegetation New radiocarbon ages from buried wood retrieved from the base of Acheron rock avalanche deposit represents an emplacement age closely post-dating (Wk 12094) 1152 ± 51 years B.P. This differs significantly from a previous radiocarbon age of (NZ547) 500 ± 69 years B.P. and modal lichenometry and weathering-rind thickness ages of approximately 460 ± 10 yrs and 490 ± 50 years B.P. The new age shows no resemblance to an earthquake event around 700- 500 years B.P. on the Porters Pass-Amberley Fault Zone. The DAN run out simulation using a friction model rheology successfully replicated the long run out and velocity of the Acheron rock avalanche using a frictron angle of 27° and high earth pressure coefficients of 5.5, 5.2, and 5.9. The elevated earth pressure coefficients represent dispersive pressures derived from dynamic fragmentation of the debris within the mobile rock avalanche, supporting the hypothesis of Davies and McSaveney (2002). The DAN model has potential applications for areas prone to large-scale instability in the elevated slopes and steep waterways of the Southern Alps. A paleoseismic investigation of a newly identified scarp of the Porters Pass Fault partially buried by the rock avalanche was conducted to identify any evidence of a coseismic relationship to the Acheron rock avalanche. This identified three-four fault traces striking at 078°, and a sag pond displaying a sequence of overbank deposits containing two buried soils representing an earthquake event horizon. A 40cm vertical offset of the ponded sediment and lower buried soil horizqn was recorded, which was dated to (Wk 13112 charcoal in palosol) 653 ± 54 years B.P. and (Wk 13034 palosol) 661 ± 34 years B.P. The evidence indicates a fault rupture occurred along the Porters Pass Fault, west of Porters Pass most likely extending to the Red Lakes terraces, post-dating 700 years B.P., resulting in 40cm of vertical displacement and an unknown component of dextral strike slip movement. This event post dates the event one (1000 ± 100 years B.P) at Porters Pass previously considered to represent the most recent rupture along the fault line. This points to a probable source for resetting of the modal weathering-rind thicknesses and lichen size populations in the Red Hill valley and possibly the Red Lakes terraces. These results suggest careful consideration must be given to the geomorphic and paleoseismic history of a specific site when applying surface dating techniques and furthermore the origin of dates used in literature and their useful range should be verified. An event at 700-500 years B.P did not trigger the Acheron rock avalanche as previously assumed supporting Howard's conclusions. The lack of similar aged rupture evidence in either of the Porters Pass and Coleridge trenches supports Howard's hypothesis of segmentation of the Porters Pass Fault; where rupture occurs along one fault segment but not along another. The new rock avalanche age closely post-dating 1200-1100 years B.P. resembles the poorly constrained event one rupture age of 1700-800 years B.P for the Porters Pass Fault and the tighter constrained Round Top event of 1010 ± 50 years B.P. on the Alpine Fault. Eight other rock avalanche deposits spread across the central Southern Alps also resemble the new ages however are unable to be assigned specific earthquake events due to the large associated error bars of± 270 years. This clustering of ages does represent compelling lines of evidence for large magnitude earthquake events occurring over the central Southern Alps. The presence of a rock avalanche deposit does not signify an earthquake based on the historical evidence in the Southern Alps however clustering of ages does suggest that large Mw >7 earthquakes occurred across the Southern Alps between 1200-900 years BP.
In this thesis, focus is given to develop methodologies for rapidly estimating specific components of loss and downtime functions. The thesis proposes methodologies for deriving loss functions by (i) considering individual component performance; (ii) grouping them as per their performance characteristics; and (iii) applying them to similar building usage categories. The degree of variation in building stock and understanding their characteristics are important factors to be considered in the loss estimation methodology and the field surveys carried out to collect data add value to the study. To facilitate developing ‘downtime’ functions, this study investigates two key components of downtime: (i) time delay from post-event damage assessment of properties; and (ii) time delay in settling the insurance claims lodged. In these two areas, this research enables understanding of critical factors that influence certain aspects of downtime and suggests approaches to quantify those factors. By scrutinising the residential damage insurance claims data provided by the Earthquake Commission (EQC) for the 2010- 2011 Canterbury Earthquake Sequence (CES), this work provides insights into various processes of claims settlement, the time taken to complete them and the EQC loss contributions to building stock in Christchurch city and Canterbury region. The study has shown diligence in investigating the EQC insurance claim data obtained from the CES to get new insights and build confidence in the models developed and the results generated. The first stage of this research develops contribution functions (probabilistic relationships between the expected losses for a wide range of building components and the building’s maximum response) for common types of claddings used in New Zealand buildings combining the probabilistic density functions (developed using the quantity of claddings measured from Christchurch buildings), fragility functions (obtained from the published literature) and cost functions (developed based on inputs from builders) through Monte Carlo simulations. From the developed contribution functions, glazing, masonry veneer, monolithic and precast concrete cladding systems are found to incur 50% loss at inter-storey drift levels equal to 0.027, 0.003, 0.005 and 0.011, respectively. Further, the maximum expected cladding loss for glazing, masonry veneer, monolithic, precast concrete cladding systems are found to be 368.2, 331.9, 365.0, and 136.2 NZD per square meter of floor area, respectively. In the second stage of this research, a detailed cost breakdown of typical buildings designed and built for different purposes is conducted. The contributions of structural and non- structural components to the total building cost are compared for buildings of different usages, and based on the similar ratios of non-structural performance group costs to the structural performance group cost, four-building groups are identified; (i) Structural components dominant group: outdoor sports, stadiums, parkings and long-span warehouses, (ii) non- structural drift-sensitive components dominant group: houses, single-storey suburban buildings (all usages), theatres/halls, workshops and clubhouses, (iii) non-structural acceleration- sensitive components dominant group: hospitals, research labs, museums and retail/cold stores, and (iv) apartments, hotels, offices, industrials, indoor sports, classrooms, devotionals and aquariums. By statistically analysing the cost breakdowns, performance group weighting factors are proposed for structural, and acceleration-sensitive and drift-sensitive non-structural components for all four building groups. Thus proposed building usage groupings and corresponding weighting factors facilitate rapid seismic loss estimation of any type of building given the EDPs at storey levels are known. A model for the quantification of post-earthquake inspection duration is developed in the third stage of this research. Herein, phase durations for the three assessment phases (one rapid impact and two rapid building) are computed using the number of buildings needing inspections, the number of engineers involved in inspections and a phase duration coefficient (which considers the median building inspection time, efficiency of engineer and the number of engineers involved in each assessment teams). The proposed model can be used: (i) by national/regional authorities to decide the length of the emergency period following a major earthquake, and estimate the number of engineers required to conduct a post-earthquake inspection within the desired emergency period, and (ii) to quantify the delay due to inspection for the downtime modelling framework. The final stage of this research investigates the repair costs and insurance claim settlement time for damaged residential buildings in the 2010-2011 Canterbury earthquake sequence. Based on the EQC claim settlement process, claims are categorized into three groups; (i) Small Claims: claims less than NZD15,000 which were settled through cash payment, (ii) Medium Claims: claims less than NZD100,000 which were managed through Canterbury Home Repair Programme (CHRP), and (iii) Large Claims: claims above NZD100,000 which were managed by an insurance provider. The regional loss ratio (RLR) for greater Christchurch for three events inducing shakings of approximate seismic intensities 6, 7, and 8 are found to be 0.013, 0.066, and 0.171, respectively. Furthermore, the claim duration (time between an event and the claim lodgement date), assessment duration (time between the claim lodgement day and the most recent assessment day), and repair duration (time between the most recent assessment day and the repair completion day) for the insured residential buildings in the region affected by the Canterbury earthquake sequence is found to be in the range of 0.5-4 weeks, 1.5- 5 months, and 1-3 years, respectively. The results of this phase will provide useful information to earthquake engineering researchers working on seismic risk/loss and insurance modelling.
Geosynthetic reinforced soil (GRS) walls involve the use of geosynthetic reinforcement (polymer material) within the retained backfill, forming a reinforced soil block where transmission of overturning and sliding forces on the wall to the backfill occurs. Key advantages of GRS systems include the reduced need for large foundations, cost reduction (up to 50%), lower environmental costs, faster construction and significantly improved seismic performance as observed in previous earthquakes. Design methods in New Zealand have not been well established and as a result, GRS structures do not have a uniform level of seismic and static resistance; hence involve different risks of failure. Further research is required to better understand the seismic behaviour of GRS structures to advance design practices. The experimental study of this research involved a series of twelve 1-g shake table tests on reduced-scale (1:5) GRS wall models using the University of Canterbury shake-table. The seismic excitation of the models was unidirectional sinusoidal input motion with a predominant frequency of 5Hz and 10s duration. Seismic excitation of the model commenced at an acceleration amplitude level of 0.1g and was incrementally increased by 0.1g in subsequent excitation levels up to failure (excessive displacement of the wall panel). The wall models were 900mm high with a full-height rigid facing panel and five layers of Microgird reinforcement (reinforcement spacing of 150mm). The wall panel toe was founded on a rigid foundation and was free to slide. The backfill deposit was constructed from dry Albany sand to a backfill relative density, Dr = 85% or 50% through model vibration. The influence of GRS wall parameters such as reinforcement length and layout, backfill density and application of a 3kPa surcharge on the backfill surface was investigated in the testing sequence. Through extensive instrumentation of the wall models, the wall facing displacements, backfill accelerations, earth pressures and reinforcement loads were recorded at the varying levels of model excitation. Additionally, backfill deformation was also measured through high-speed imaging and Geotechnical Particle Image Velocimetry (GeoPIV) analysis. The GeoPIV analysis enabled the identification of the evolution of shear strains and volumetric strains within the backfill at low strain levels before failure of the wall thus allowing interpretations to be made regarding the strain development and shear band progression within the retained backfill. Rotation about the wall toe was the predominant failure mechanism in all excitation level with sliding only significant in the last two excitation levels, resulting in a bi-linear displacement acceleration curve. An increase in acceleration amplification with increasing excitation was observed with amplification factors of up to 1.5 recorded. Maximum seismic and static horizontal earth pressures were recorded at failure and were recorded at the wall toe. The highest reinforcement load was recorded at the lowest (deepest in the backfill) reinforcement layer with a decrease in peak load observed at failure, possibly due to pullout failure of the reinforcement layer. Conversely, peak reinforcement load was recorded at failure for the top reinforcement layer. The staggered reinforcement models exhibited greater wall stability than the uniform reinforcement models of L/H=0.75. However, similar critical accelerations were determined for the two wall models due to the coarseness of excitation level increments of 0.1g. The extended top reinforcements were found to restrict the rotational component of displacement and prevented the development of a preliminary shear band at the middle reinforcement layer, contributing positively to wall stability. Lower acceleration amplification factors were determined for the longer uniform reinforcement length models due to reduced model deformation. A greater distribution of reinforcement load towards the top two extended reinforcement layers was also observed in the staggered wall models. An increase in model backfill density was observed to result in greater wall stability than an increase in uniform reinforcement length. Greater acceleration amplification was observed in looser backfill models due to their lower model stiffness. Due to greater confinement of the reinforcement layers, greater reinforcement loads were developed in higher density wall models with less wall movement required to engage the reinforcement layers and mobilise their resistance. The application of surcharge on the backfill was observed to initially increase the wall stability due to greater normal stresses within the backfill but at greater excitation levels, the surcharge contribution to wall destabilising inertial forces outweighs its contribution to wall stability. As a result, no clear influence of surcharge on the critical acceleration of the wall models was observed. Lower acceleration amplification factors were observed for the surcharged models as the surcharge acts as a damper during excitation. The application of the surcharge also increases the magnitude of reinforcement load developed due to greater confinement and increased wall destabilising forces. The rotation of the wall panel resulted in the progressive development of shears surface with depth that extended from the backfill surface to the ends of the reinforcement (edge of the reinforced soil block). The resultant failure plane would have extended from the backfill surface to the lowest reinforcement layer before developing at the toe of the wall, forming a two-wedge failure mechanism. This is confirmed by development of failure planes at the lowest reinforcement layer (deepest with the backfill) and at the wall toe observed at the critical acceleration level. Key observations of the effect of different wall parameters from the GeoPIV results are found to be in good agreement with conclusions developed from the other forms of instrumentation. Further research is required to achieve the goal of developing seismic guidelines for GRS walls in geotechnical structures in New Zealand. This includes developing and testing wall models with a different facing type (segmental or wrap-around facing), load cell instrumentation of all reinforcement layers, dynamic loading on the wall panel and the use of local soils as the backfill material. Lastly, the limitations of the experimental procedure and wall models should be understood.