Many buildings with relatively low damage from the 2010-2011 Canterbury were deemed uneconomic to repair and were replaced [1,2]. Factors that affected commercial building owners’ decisions to replace rather than repair, included capital availability, uncertainty with regards to regional recovery, local market conditions and ability to generate cash flow, and repair delays due to limited property access (cordon). This poster provides a framework for modeling decision-making in a case where repair is feasible but replacement might offer greater economic value – a situation not currently modeled in engineering risk analysis.
The 2010-2011 Canterbury earthquake sequence, and the resulting extensive data sets on damaged buildings that have been collected, provide a unique opportunity to exercise and evaluate previously published seismic performance assessment procedures. This poster provides an overview of the authors’ methodology to perform evaluations with two such assessment procedures, namely the P-58 guidelines and the REDi Rating System. P-58, produced by the Federal Emergency Management Agency (FEMA) in the United States, aims to facilitate risk assessment and decision-making by quantifying earthquake ground shaking, structural demands, component damage and resulting consequences in a logical framework. The REDi framework, developed by the engineering firm ARUP, aids stakeholders in implementing resilience-based earthquake design. Preliminary results from the evaluations are presented. These have the potential to provide insights on the ability of the assessment procedures to predict impacts using “real-world” data. However, further work remains to critically analyse these results and to broaden the scope of buildings studied and of impacts predicted.
The Canterbury Earthquake Sequence (CES) of 2010-2011 produced large seismic moments up to Mw 7.1. These large, near-to-surface (<15 km) ruptures triggered >6,000 rockfall boulders on the Port Hills of Christchurch, many of which impacted houses and affected the livelihoods of people within the impacted area. From these disastrous and unpredicted natural events a need arose to be able to assess the areas affected by rockfall events in the future, where it is known that a rockfall is possible from a specific source outcrop but the potential boulder runout and dynamics are not understood. The distribution of rockfall deposits is largely constrained by the physical properties and processes of the boulder and its motion such as block density, shape and size, block velocity, bounce height, impact and rebound angle, as well as the properties of the substrate. Numerical rockfall models go some way to accounting for all the complex factors in an algorithm, commonly parameterised in a user interface where site-specific effects can be calibrated. Calibration of these algorithms requires thorough field checks and often experimental practises. The purpose of this project, which began immediately following the most destructive rupture of the CES (February 22, 2011), is to collate data to characterise boulder falls, and to use this information, supplemented by a set of anthropogenic boulder fall data, to perform an in-depth calibration of the three-dimensional numerical rockfall model RAMMS::Rockfall. The thesis covers the following topics: • Use of field data to calibrate RAMMS. Boulder impact trails in the loess-colluvium soils at Rapaki Bay have been used to estimate ranges of boulder velocities and bounce heights. RAMMS results replicate field data closely; it is concluded that the model is appropriate for analysing the earthquake-triggered boulder trails at Rapaki Bay, and that it can be usefully applied to rockfall trajectory and hazard assessment at this and similar sites elsewhere. • Detailed analysis of dynamic rockfall processes, interpreted from recorded boulder rolling experiments, and compared to RAMMS simulated results at the same site. Recorded rotational and translational velocities of a particular boulder show that the boulder behaves logically and dynamically on impact with different substrate types. Simulations show that seasonal changes in soil moisture alter rockfall dynamics and runout predictions within RAMMS, and adjustments are made to the calibration to reflect this; suggesting that in hazard analysis a rockfall model should be calibrated to dry rather than wet soil conditions to anticipate the most serious outcome. • Verifying the model calibration for a separate site on the Port Hills. The results of the RAMMS simulations show the effectiveness of calibration against a real data set, as well as the effectiveness of vegetation as a rockfall barrier/retardant. The results of simulations are compared using hazard maps, where the maximum runouts match well the mapped CES fallen boulder maximum runouts. The results of the simulations in terms of frequency distribution of deposit locations on the slope are also compared with those of the CES data, using the shadow angle tool to apportion slope zones. These results also replicate real field data well. Results show that a maximum runout envelope can be mapped, as well as frequency distribution of deposited boulders for hazard (and thus risk) analysis purposes. The accuracy of the rockfall runout envelope and frequency distribution can be improved by comprehensive vegetation and substrate mapping. The topics above define the scope of the project, limiting the focus to rockfall processes on the Port Hills, and implications for model calibration for the wider scientific community. The results provide a useful rockfall analysis methodology with a defensible and replicable calibration process, that has the potential to be applied to other lithologies and substrates. Its applications include a method of analysis for the selection and positioning of rockfall countermeasure design; site safety assessment for scaling and demolition works; and risk analysis and land planning for future construction in Christchurch.
Predicting building collapse due to seismic motion is critical in design and more so after a major event. Damaged structures can appear sound, but collapse under following major events. There can thus be significant risk in decision making after a major seismic event concerning the safe occupation of a building or surrounding areas, versus the unknown impact of unknown major aftershocks. Model-based pushover analyses are effective if the structural properties are well understood, which is not valid post-event when this risk information is most useful. This research combines Hysteresis Loop Analysis (HLA) structural health monitoring (SHM) and Incremental Dynamic Analysis (IDA) methods to determine collapse capacity and probability of collapse for a specific structure, at any time, a range of earthquake excitations to ensure robustness. The nonlinear dynamic analysis method presented enables constant updating of building performance predictions using post-event SHM results. The resulting combined methods provide near real-time updating of collapse fragility curves as events progress, quantifying the change of collapse probability or seismic induced losses for decision-making - a novel, higher resolution risk analysis than previously available. The methods are not computationally expensive and there is no requirement for a validated numerical model. Results show significant potential benefits and a clear evolution of risk. They also show clear need for extending SHM toward creating improved predictive models for analysis of subsequent events, where the Christchurch series of 2010-2011 had significant post-event aftershocks after each main event. Finally, the overall method is generalisable to any typical engineering demand parameter.
On 22 February 2011，a magnitude Mw 6.3 earthquake occurred with an epicenter located near Lyttelton at about 10km from Christchurch in Canterbury region on the South Island of New Zealand (Figure 1). Since this earthquake occurred in the midst of the aftershock activity which had continued since the 4 September 2010 Darfield Earthquake occurrence, it was considered to be an aftershock of the initial earthquake. Because of the short distance to the city and the shallower depth of the epicenter, this earthquake caused more significant damage to pipelines, traffic facilities, residential houses/properties and multi-story buildings in the central business district than the September 2010 Darfield Earthquake in spite of its smaller earthquake magnitude. Unfortunately, this earthquake resulted in significant number of casualties due to the collapse of multi-story buildings and unreinforced masonry structures in the city center of Christchurch. As of 4 April, 172 casualties were reported and the final death toll is expected to be 181. While it is extremely regrettable that Christchurch suffered a terrible number of victims, civil and geotechnical engineers have this hard-to-find opportunity to learn the response of real ground from two gigantic earthquakes which occurred in less than six months from each other. From geotechnical engineering point of view, it is interesting to discuss the widespread liquefaction in natural sediments, repeated liquefaction within short period and further damage to earth structures which have been damaged in the previous earthquake. Following the earthquake, an intensive geotechnical reconnaissance was conducted to capture evidence and perishable data from this event. The team included the following members: Misko Cubrinovski (University of Canterbury, NZ, Team Leader), Susumu Yasuda (Tokyo Denki University, Japan, JGS Team Leader), Rolando Orense (University of Auckland, NZ), Kohji Tokimatsu (Tokyo Institute of Technology, Japan), Ryosuke Uzuoka (Tokushima University, Japan), Takashi Kiyota (University of Tokyo, Japan), Yasuyo Hosono (Toyohashi University of Technology, Japan) and Suguru Yamada (University of Tokyo, Japan).
© 2019, Springer-Verlag GmbH Germany, part of Springer Nature. Prediction of building collapse due to significant seismic motion is a principle objective of earthquake engineers, particularly after a major seismic event when the structure is damaged and decisions may need to be made rapidly concerning the safe occupation of a building or surrounding areas. Traditional model-based pushover analyses are effective, but only if the structural properties are well understood, which is not the case after an event when that information is most useful. This paper combines hysteresis loop analysis (HLA) structural health monitoring (SHM) and incremental dynamic analysis (IDA) methods to identify and then analyse collapse capacity and the probability of collapse for a specific structure, at any time, a range of earthquake excitations to ensure robustness. This nonlinear dynamic analysis enables constant updating of building performance predictions following a given and subsequent earthquake events, which can result in difficult to identify deterioration of structural components and their resulting capacity, all of which is far more difficult using static pushover analysis. The combined methods and analysis provide near real-time updating of the collapse fragility curves as events progress, thus quantifying the change of collapse probability or seismic induced losses very soon after an earthquake for decision-making. Thus, this combination of methods enables a novel, higher-resolution analysis of risk that was not previously available. The methods are not computationally expensive and there is no requirement for a validated numerical model, thus providing a relatively simpler means of assessing collapse probability immediately post-event when such speed can provide better information for critical decision-making. Finally, the results also show a clear need to extend the area of SHM toward creating improved predictive models for analysis of subsequent events, where the Christchurch series of 2010–2011 had significant post-event aftershocks.
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.
This dissertation addresses a diverse range of applied aspects in ground motion simulation validation via the response of complex structures. In particular, the following topics are addressed: (i) the investigation of similarity between recorded and simulated ground motions using code-based 3D irregular structural response analysis, (ii) the development of a framework for ground motion simulations validation to identify the cause of differences between paired observed and simulated dataset, and (iii) the illustration of the process of using simulations for seismic performance-based assessment. The application of simulated ground motions is evaluated for utilisation in engineering practice by considering responses of 3D irregular structures. Validation is performed in a code-based context when the NZS1170.5 (NZS1170.5:2004, 2004) provisions are followed for response history analysis. Two real buildings designed by engineers and physically constructed in Christchurch before the 2010-2011 Canterbury earthquake sequence are considered. The responses are compared when the buildings are subjected to 40 scaled recorded and their subsequent simulated ground motions selected from 22 February 2011 Christchurch. The similarity of recorded and simulated responses is examined using statistical methods such as bootstrapping and hypothesis testing to determine whether the differences are statistically significant. The findings demonstrate the applicability of simulated ground motion when the code-based approach is followed in response history analysis. A conceptual framework is developed to link the differences between the structural response subjected to simulated and recorded ground motions to the differences in their corresponding intensity measures. This framework allows the variability to be partitioned into the proportion that can be “explained” by the differences in ground motion intensity measures and the remaining “unexplained” variability that can be attributed to different complexities such as dynamic phasing of multi-mode response, nonlinearity, and torsion. The application of this framework is examined through a hierarchy of structures reflecting a range of complexity from single-degree-of-freedom to 3D multi-degree-of-freedom systems with different materials, dynamic properties, and structural systems. The study results suggest the areas that ground motion simulation should focus on to improve simulations by prioritising the ground motion intensity measures that most clearly account for the discrepancies in simple to complex structural responses. Three approaches are presented to consider recorded or simulated ground motions within the seismic performance-based assessment framework. Considering the applications of ground motions in hazard and response history analyses, different pathways in utilising ground motions in both areas are explored. Recorded ground motions are drawn from a global database (i.e., NGA-West2 Ancheta et al., 2014). The NZ CyberShake dataset is used to obtain simulations. Advanced ground motion selection techniques (i.e., generalized conditional intensity measure, GCIM) are used for ground motion selection at a few intensity levels. The comparison is performed by investigating the response of an example structure (i.e., 12-storey reinforced concrete special moment frame) located in South Island, NZ. Results are compared and contrasted in terms of hazard, groundmotion selection, structural responses, demand hazard, and collapse risk, then, the probable reasons for differences are discussed. The findings from this study highlight the present opportunities and shortcomings in using simulations in risk assessment. i
