The aim of this poster is to examine the seismic response of two structural systems when subjected to observed and simulated ground motions (GMs) for the 22 February 2011 (22Feb2011) Christchurch earthquake (Razafindrakoto et al. (2018)) via an automated workflow. The layout and technical details of the automated workflow are described at Motha et. al. (2019).
Validation is an essential step to assess the applicability of simulated ground motions for utilization in engineering practice, and a comprehensive analysis should include both simple intensity measures (PGA, SA, etc), as well as the seismic response of a range of complex systems obtained by response history analysis. In order to enable a spectrum of complex structural systems to be considered in systematic validation of ground motion simulations in a routine fashion, an automated workflow was developed. Such a workflow enables validation of simulated ground motions in terms of different complex model responses by considering various ground motion sets and different ground motion simulation methods. The automated workflow converts the complex validation process into a routine one by providing a platform to perform the validation process promptly as a built-in process of simulation post-processing. As a case study, validation of simulated ground motions was investigated via the automated workflow by comparing the dynamic responses of three steel special moment frame (SMRF) subjected to the 40 observed and 40 simulated ground motions of 22 February 2011 Christchurch earthquake. The seismic responses of the structures are principally quantified via the peak floor acceleration and maximum inter-storey drift ratio. Overall, the results indicate a general agreement in seismic demands obtained using the recorded and simulated ensembles of ground motions and provide further evidence that simulated ground motions can be used in code-based structural performance assessments in-place of, or in combination with, ensembles of recorded ground motions.
Results from a series of 1D seismic effective stress analyses of natural soil deposits from Christchurch are summarized. The analysed soil columns include sites whose performance during the 2010-2011 Canterbury earthquakes varied significantly, from no liquefaction manifestation at the ground surface to very severe liquefaction, in which case a large area of the site was covered by thick soil ejecta. Key soil profile characteristics and response mechanisms affecting the severity of surface liquefaction manifestation and subsequent damage are explored. The influence of shaking intensity on the triggering and contribution of these mechanisms is also discussed. Careful examination of the results highlights the importance of considering the deposit as a whole, i.e. a system of layers, including interactions between layers in the dynamic response and through pore water pressure redistribution and water flow.
Validating dynamic responses of engineered systems subjected to simulated ground motions is essential in scrutinising the applicability of simulated ground motions for engineering demand analyses. This paper compares the responses of two 3D building models subjected to recorded and simulated ground motions scaled to the NZS1170.5 design response spectrum, in order to evaluate the applicability of simulated ground motions for use in conventional engineering practice in New Zealand. The buildings were designed according to the NZS1170.5 and physically constructed in Christchurch prior to the 2010-2011 Canterbury earthquakes. 40 recorded ground motions from the 22 February 2011 Christchurch earthquake, along with the simulated ground motions for this event from Razafindrakoto et al. (2018) are considered. The seismic responses of the structures are principally quantified via the peak floor acceleration and maximum inter-storey drift ratio. Overall, the results indicate a general agreement in seismic demands obtained using the recorded and simulated ensembles of ground motions and provide further evidence that simulated ground motions using state-of-the-art methods can be used in code-based structural performance assessments inplace of, or in combination with, ensembles of recorded ground motions.
Despite the relatively low seismicity, a large earthquake in the Waikato region is expected to have a high impact, when the fourth-largest regional population and economy and the high density critical infrastructure systems in this region are considered. Furthermore, Waikato has a deep soft sedimentary basin, which increases the regional seismic hazard due to trapping and amplification of seismic waves and generation of localized surface waves within the basin. This phenomenon is known as the “Basin Effect”, and has been attributed to the increased damage in several historic earthquakes, including the 2010-2011 Canterbury earthquakes. In order to quantitatively model the basin response and improve the understanding of regional seismic hazard, geophysical methods will be used to develop shear wave velocity profiles across the Waikato basin. Active surface wave methods involve the deployment of linear arrays of geophones to record the surface waves generated by a sledge hammer. Passive surface wave methods involve the deployment of two-dimensional seismometer arrays to record ambient vibrations. At each site, the planned testing includes one active test and two to four passive arrays. The obtained data are processed to develop dispersion curves, which describe surface wave propagation velocity as a function of frequency (or wavelength). Dispersion curves are then inverted using the Geopsy software package to develop a suite of shear wave velocity profiles. Currently, more than ten sites in Waikato are under consideration for this project. This poster presents the preliminary results from the two sites that have been tested. The shear wave velocity profiles from all sites will be used to produce a 3D velocity model for the Waikato basin, a part of QuakeCoRE flagship programme 1.
Buildings subject to earthquake shaking will tend to move not only horizontally but also rotate in plan. In-plan rotation is known as “building torsion” and it may occur for a variety of reasons, including stiffness and strength eccentricity and/or torsional effects from ground motions. Methods to consider torsion in structural design standards generally involve analysis of the structure in its elastic state. This is despite the fact that the structural elements can yield, thereby significantly altering the building response and the structural element demands. If demands become too large, the structure may collapse. While a number of studies have been conducted into the behavior of structures considering inelastic building torsion, there appears to be no consensus that one method is better than another and as a result, provisions within current design standards have not adopted recent proposals in the literature. However, the Canterbury Earthquakes Royal Commission recently made the recommendation that provisions to account for inelastic torsional response of buildings be introduced within New Zealand building standards. Consequently, this study examines how and to what extent the torsional response due to system eccentricity may affect the seismic performance of a building and considers what a simple design method should account for. It is concluded that new methods should be simple, be applicable to both the elastic and inelastic range of response, consider bidirectional excitation and include guidance for multi-story systems.
The ultimate goal of this study is to develop a model representing the in-plane behaviour of plasterboard ceiling diaphragms, as part of the efforts towards performance-based seismic engineering of low-rise light timber-framed (LTF) residential buildings in New Zealand (NZ). LTF residential buildings in NZ are constructed according to a prescriptive standard – NZS 3604 Timberframed buildings [1]. With regards to seismic resisting systems, LTF buildings constructed to NZS3604 often have irregular bracing arrangements within a floor plane. A damage survey of LTF buildings after the Canterbury earthquake revealed that structural irregularity (irregular bracing arrangement within a plane) significantly exacerbated the earthquake damage to LTF buildings. When a building has irregular bracing arrangements, the building will have not only translational deflections but also a torsional response in earthquakes. How effectively the induced torsion can be resolved depends on the stiffness of the floors/roof diaphragms. Ceiling and floor diaphragms in LTF buildings in NZ have different construction details from the rest of the world and there appears to be no information available on timber diaphragms typical of NZ practice. This paper presents experimental studies undertaken on plasterboard ceiling diaphragms as typical of NZ residential practice. Based on the test results, a mathematical model simulating the in-plane stiffness of plasterboard ceiling diaphragms was developed, and the developed model has a similar format to that of plasterboard bracing wall elements presented in an accompany paper by Liu [2]. With these two models, three-dimensional non-linear push-over studies of LTF buildings can be undertaken to calculate seismic performance of irregular LTF buildings.
This dissertation addresses a diverse range of topics in the physics-based broadband ground motion simulation, with a focus on New Zealand applications. In particular the following topics are addressed: the methodology and computational implementation of a New Zealand Velocity Model for broadband ground motion simulation; generalised parametric functions and spatial correlations for seismic velocities in the Canterbury, New Zealand region from surface-wave-based site characterisation; and ground motion simulations of Hope Fault earthquakes. The paragraphs below outline each contribution in more detail. A necessary component in physics-based ground motion simulation is a 3D model which details the seismic velocities in the region of interest. Here a velocity model construction methodology, its computational implementation, and application in the construction of a New Zealand velocity model for use in physics-based broadband ground motion simulation are presented. The methodology utilises multiple datasets spanning different length scales, which is enabled via the use of modular sub-regions, geologic surfaces, and parametric representations of crustal velocity. A number of efficiency-related workflows to decrease the overall computational construction time are employed, while maintaining the flexibility and extensibility to incorporate additional datasets and re- fined velocity parameterizations as they become available. The model comprises explicit representations of the Canterbury, Wellington, Nelson-Tasman, Kaikoura, Marlborough, Waiau, Hanmer and Cheviot sedimentary basins embedded within a regional travel-time tomography-based velocity model for the shallow crust and provides the means to conduct ground motion simulations throughout New Zealand for the first time. Recently developed deep shear-wave velocity profiles in Canterbury enabled models that better characterise the velocity structure within geologic layers of the Canterbury sedimentary basin to be developed. Here the development of depth- and Vs30-dependent para-metric velocity and spatial correlation models to characterise shear-wave velocities within the geologic layers of the Canterbury sedimentary basin are presented. The models utilise data from 22 shear-wave velocity profiles of up to 2.5km depth (derived from surface wave analysis) juxtaposed with models which detail the three-dimensional structure of the geologic formations in the Canterbury sedimentary basin. Parametric velocity equations are presented for Fine Grained Sediments, Gravels, and Tertiary layer groupings. Spatial correlations were developed and applied to generate three-dimensional stochastic velocity perturbations. Collectively, these models enable seismic velocities to be realistically represented for applications such as 3D ground motion and site response simulations. Lastly the New Zealand velocity model is applied to simulate ground motions for a Mw7.51 rupture of the Hope Fault using a physics-based simulation methodology and a 3D crustal velocity model of New Zealand. The simulation methodology was validated for use in the region through comparison with observations for a suite of historic small magnitude earthquakes located proximal to the Hope Fault. Simulations are compared with conventionally utilised empirical ground motion models, with simulated peak ground velocities being notably higher in regions with modelled sedimentary basins. A sensitivity analysis was undertaken where the source characteristics of magnitude, stress parameter, hypocentre location and kinematic slip distribution were varied and an analysis of their effect on ground motion intensities is presented. It was found that the magnitude and stress parameter strongly influenced long and short period ground motion amplitudes, respectively. Ground motion intensities for the Hope Fault scenario are compared with the 2016 Kaikoura Mw7.8 earthquake, it was found that the Kaikoura earthquake produced stronger motions along the eastern South Island, while the Hope Fault scenario resulted in stronger motions immediately West of the near-fault region. The simulated ground motions for this scenario complement prior empirically-based estimates and are informative for mitigation and emergency planning purposes.
Beach ridge stratigraphy can provide an important record of both sustained coastal progradation and responses to events such as extreme storms, as well as evidence of earthquake induced sediment pulses. This study is a stratigraphic investigation of the late Holocene mixed sand gravel (MSG) beach ridge plain on the Canterbury coast, New Zealand. The subsurface was imaged along a 370 m shore-normal transect using 100 and 200 MHz ground penetrating radar (GPR) antennae, and cored to sample sediment textures. Results show that, seaward of a back-barrier lagoon, the Pegasus Bay beach ridge plain prograded almost uniformly, under conditions of relatively stable sea level. Nearshore sediment supply appears to have created a sustained sediment surplus, perhaps as a result of post-seismic sediment pulses, resulting in a flat, morphologically featureless beach ridge plain. Evidence of a high magnitude storm provides an exception, with an estimated event return period in excess of 100 years. Evidence from the GPR sequence combined with modern process observations from MSG beaches indicates that a paleo storm initially created a washover fan into the back-barrier lagoon, with a large amount of sediment simultaneously moved off the beach face into the nearshore. This erosion event resulted in a topographic depression still evident today. In the subsequent recovery period, sediment was reworked by swash onto the beach as a sequence of berm deposit laminations, creating an elevated beach ridge that also has a modern-day topographic signature. As sediment supply returned to normal, and under conditions of falling sea level, a beach ridge progradation sequence accumulated seaward of the storm feature out to the modern-day beach as a large flat, uniform progradation plain. This study highlights the importance of extreme storm events and earthquake pulses on MSG coastlines in triggering high volume beach ridge formation during the subsequent recovery period.