We present initial results from a set of three-dimensional (3D) deterministic earthquake ground motion simulations for the northern Canterbury plains, Christchurch and the Banks Peninsula region, which explicitly incorporate the effects of the surface topography. The simu-lations are done using Hercules, an octree-based finite-element parallel software for solving 3D seismic wave propagation problems in heterogeneous media under kinematic faulting. We describe the efforts undertaken to couple Hercules with the South Island Velocity Model (SIVM), which included changes to the SIVM code in order to allow for single repetitive que-ries and thus achieve a seamless finite-element meshing process within the end-to-end ap-proach adopted in Hercules. We present our selection of the region of interest, which corre-sponds to an area of about 120 km × 120 km, with the 3D model reaching a depth of 60 km. Initial simulation parameters are set for relatively high minimum shear wave velocity and a low maximum frequency, which we are progressively scaling up as computing resources permit. While the effects of topography are typically more important at higher frequencies and low seismic velocities, even at this initial stage of our efforts (with a maximum of 2 Hz and a mini-mum of 500 m/s), it is possible to observe the importance of the topography in the response of some key locations within our model. To highlight these effects we compare the results of the 3D topographic model with respect to those of a flat (squashed) 3D model. We draw rele-vant conclusions from the study of topographic effects during earthquakes for this region and describe our plans for future work.
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 2010-2011 Christchurch earthquakes generated damage in several Reinforced Concrete (RC) buildings, which had RC walls as the principal resistant element against earthquake demand. Despite the agreement between structural engineers and researchers in an overall successfully performance there was a lack of knowledge about the behaviour of the damaged structures, and even deeper about a repaired structure, which triggers arguments between different parties that remains up to these days. Then, it is necessary to understand the capacity of the buildings after the earthquake and see how simple repairs techniques improve the building performance. This study will assess the residual capacity of ductile slender RC walls according to current standards in New Zealand, NZS 3101.1 2006 A3. First, a Repaired RC walls Database is created trying to gather previous studies and to evaluate them with existing international guidelines. Then, an archetype building is designed, and the wall is extracted and scaled. Four half-scale walls were designed and will be constructed and tested at the Structures Testing Laboratory at The University of Auckland. The overall dimensions are 3 [m] height, 2 [m] length and 0.175 [m] thick. All four walls will be identical, with differences in the loading protocol and the presence or absence of a repair technique. Results are going to be useful to assess the residual capacity of a damaged wall compare to the original behaviour and also the repaired capacity of walls with simpler repair techniques. The expected behaviour is focussed on big changes in stiffness, more evident than in previously tested RC beams found in the literature.
Eccentrically Braced Frames (EBFs) are a widely used seismic resisting structural steel system. Since their inception in the late 1970s, they have been a viable option with an available stiffness that is between simple braced systems and moment resisting systems. A similar concept, the linked column frame (LCF), uses shear links between two closely spaced columns. In both cases, the key component is the active link or the shear link, and this component is the objective of this study. The performance of high rise EBF buildings in the 2010 and 2011 Christchurch earthquakes was beyond that which was expected, especially considering the very high accelerations recorded. As the concrete high-rises were torn down, two EBF buildings remained standing and only required some structural repair. These events prompted a renewed interest in bolted shear links, as well as their performance. While some research into replaceable shear links had already been done (Mansour, 2011), the objectives of this study were to improve on the shear link itself, with the consideration that links built in the future are likely to be bolted. The main components of this study were to: 1. Reduce or eliminate the requirements for intermediate web stiffeners, as they were suspected of being detrimental to performance. Furthermore, any reduction in stiffening requirements is a direct fabrication cost saving. Links with low web aspect ratios were found to achieve exceptional ductilities when no stiffeners were included, prompting new design equations. 2. Ensure that the stresses in the ends of links are adequately transferred into the endplates without causing fractures. Although most of the experimental links had web doubler plates included, four had varied lengths of such doubler plates from 0.0 in. to 8.0 in. The link without any doubler plates performed to a similar level to its peers, and thus it is likely that links with quality end details may not need web doubler plates at all. 3. Evaluate the performance of a link with double sided stiffeners without the use of web welds, as opposed to conventional single sided, welded stiffeners. This link performed well, and web-weld-less double sided stiffeners may be an economical alternative to conventional stiffeners for deeper sections of links. 4. Evaluate the performance of a link with thin endplates that are made efficient with the use of gusset plates. This link performed to an acceptable level and provides evidence for a cost effective alternative to thick endplates, especially considering the high overstrength end moments in links, typically requiring 16-bolt connections. 5. Examine the potential use of an alternative EBF arrangement where the collector beam is over sized, and the link section is formed by cutting out parts of the beam's web. After running a series of finite element models each with a unique variation, a number of approximate design rules were derived such that future research could develop this idea further experimentally. 6. Ensure that during testing, the secondary elements (members that are not the shear link), do not yield and are not close to yielding. None of the instrumented elements experienced any unexpected yielding, however the concerns for high stresses in the collector beam panel zone during design were warranted. The use of an existing New Zealand design equation is recommended as an extra check for design codes worldwide. The above objectives were mainly conducted experimentally, except: the data set for item 1 was greatly expanded through the use of a calibrated numerical model which was then used in an extensive parametric study; item 5 was purely finite element based; and, a small parametric study was included for item 3 in an attempt to expand on the trends found there.
Tsunami events including the 2004 Indian Ocean Tsunami and the 2011 Tohoku Earthquake and Tsunami confirmed the need for Pacific-wide comprehensive risk mitigation and effective tsunami evacuation planning. New Zealand is highly exposed to tsunamis and continues to invest in tsunami risk awareness, readiness and response across the emergency management and science sectors. Evacuation is a vital risk reduction strategy for preventing tsunami casualties. Understanding how people respond to warnings and natural cues is an important element to improving evacuation modelling techniques. The relative rarity of tsunami events locally in Canterbury and also globally, means there is limited knowledge on tsunami evacuation behaviour, and tsunami evacuation planning has been largely informed by hurricane evacuations. This research aims to address this gap by analysing evacuation behaviour and movements of Kaikōura and Southshore/New Brighton (coastal suburb of Christchurch) residents following the 2016 Kaikōura earthquake. Stage 1 of the research is engaging with both these communities and relevant hazard management agencies, using a survey and community workshops to understand real-event evacuation behaviour during the 2016 Kaikōura earthquake and subsequent tsunami evacuations. The second stage is using the findings from stage 1 to inform an agent-based tsunami evacuation model, which is an approach that simulates of the movement of people during an evacuation response. This method improves on other evacuation modelling approaches to estimate evacuation times due to better representation of local population characteristics. The information provided by the communities will inform rules and interactions such as traffic congestion, evacuation delay times and routes taken to develop realistic tsunami evacuation models. This will allow emergency managers to more effectively prepare communities for future tsunami events, and will highlight recommended actions to increase the safety and efficiency of future tsunami evacuations.
This poster presents work to date on ground motion simulation validation and inversion for the Canterbury, New Zealand region. Recent developments have focused on the collection of different earthquake sources and the verification of the SPECFEM3D software package in forward and inverse simulations. SPECFEM3D is an open source software package which simulates seismic wave propagation and performs adjoint tomography based upon the spectral-element method. Figure 2: Fence diagrams of shear wave velocities highlighting the salient features of the (a) 1D Canterbury velocity model, and (b) 3D Canterbury velocity model. Figure 5: Seismic sources and strong motion stations in the South Island of New Zealand, and corresponding ray paths of observed ground motions. Figure 3: Domain used for the 19th October 2010 Mw 4.8 case study event including the location of the seismic source and strong motion stations. By understanding the predictive and inversion capabilities of SPECFEM3D, the current 3D Canterbury Velocity Model can be iteratively improved to better predict the observed ground motions. This is achieved by minimizing the misfit between observed and simulated ground motions using the built-in optimization algorithm. Figure 1 shows the Canterbury Velocity Model domain considered including the locations of small-to-moderate Mw events [3-4.5], strong motion stations, and ray paths of observed ground motions. The area covered by the ray paths essentially indicates the area of the model which will be most affected by the waveform inversion. The seismic sources used in the ground motion simulations are centroid moment tensor solutions obtained from GeoNet. All earthquake ruptures are modelled as point sources with a Gaussian source time function. The minimum Mw limit is enforced to ensure good signal-to-noise ratio and well constrained source parameters. The maximum Mw limit is enforced to ensure the point source approximation is valid and to minimize off-fault nonlinear effects.
As damage and loss caused by natural hazards have increased worldwide over the past several decades, it is important for governments and aid agencies to have tools that enable effective post-disaster livelihood recovery to create self-sufficiency for the affected population. This study introduces a framework of critical components that constitute livelihood recovery and the critical factors that lead to people’s livelihood recovery. A comparative case study is employed in this research, combined with questionnaire surveys and interviews with those communities affected by large earthquakes in Lushan, China and in Christchurch and Kaikōura, New Zealand. In Lushan, China, a framework with four livelihood components was established, namely, housing, employment, wellbeing and external assistance. Respondents considered recovery of their housing to be the most essential element for livelihood diversification. External assistance was also rated highly in assisting with their livelihood recovery. Family ties and social connections seemed to have played a larger role than that of government agencies and NGOs. However, the recovery of livelihood cannot be fully achieved without wellbeing aspects being taken into account, and people believed that quality of life and their physical and mental health were essential for livelihood restoration. In Christchurch, New Zealand, the identified livelihood components were validated through in-depth interviews. The results showed that the above framework presenting what constitutes successful livelihood recovery could also be applied in Christchurch. This study also identified the critical factors to affect livelihood recovery following the Lushan and Kaikōura earthquakes, and these include community safety, availability of family support, level of community cohesion, long-term livelihood support, external housing recovery support, level of housing recovery and availability of health and wellbeing support. The framework developed will provide guidance for policy makers and aid agencies to prioritise their strategies and initiatives in assisting people to reinstate their livelihood in a timely manner post-disaster. It will also assist the policy makers and practitioners in China and New Zealand by setting an agenda for preparing for livelihood recovery in non-urgent times so the economic impact and livelihood disruption of those affected can be effectively mitigated.
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 full scale, in-situ investigations of instrumented buildings present an excellent opportunity to observe their dynamic response in as-built environment, which includes all the real physical properties of a structure under study and its surroundings. The recorded responses can be used for better understanding of behavior of structures by extracting their dynamic characteristics. It is significantly valuable to examine the behavior of buildings under different excitation scenarios. The trends in dynamic characteristics, such as modal frequencies and damping ratios, thus developed can provide quantitative data for the variations in the behavior of buildings. Moreover, such studies provide invaluable information for the development and calibration of realistic models for the prediction of seismic response of structures in model updating and structural health monitoring studies. This thesis comprises two parts. The first part presents an evaluation of seismic responses of two instrumented three storey RC buildings under a selection of 50 earthquakes and behavioral changes after Ms=7.1 Darfield (2010) and Ms=6.3 Christchurch (2011) earthquakes for an instrumented eight story RC building. The dynamic characteristics of the instrumented buildings were identified using state-of-the-art N4SID system identification technique. Seismic response trends were developed for the three storey instrumented buildings in light of the identified frequencies and the peak response accelerations (PRA). Frequencies were observed to decrease with excitation level while no trends are discernible for the damping ratios. Soil-structure interaction (SSI) effects were also determined to ascertain their contribution in the seismic response. For the eight storey building, it was found through system identification that strong nonlinearities in the structural response occurred and manifested themselves in all identified natural frequencies of the building that exhibited a marked decrease during the strong motion duration compared to the pre-Darfield earthquakes. Evidence of foundation rocking was also found that led to a slight decrease in the identified modal frequencies. Permanent stiffness loss was also observed after the strong motion events. The second part constitutes developing and calibrating finite element model (FEM) of the instrumented three storey RC building with a shear core. A three dimensional FEM of the building is developed in stages to analyze the effect of structural, non-structural components (NSCs) and SSI on the building dynamics. Further to accurately replicate the response of the building following the response trends developed in the first part of the thesis, sensitivity based model updating technique was applied. The FEMs were calibrated by tuning the updating parameters which are stiffnesses of concrete, NSCs and soil. The updating parameters were found to generally follow decreasing trends with the excitation level. Finally, the updated FEM was used in time history analyses to assess the building seismic performance at the serviceability limit state shaking. Overall, this research will contribute towards better understanding and prediction of the behavior of structures subjected to ground motion.
