Equivalent force control method for generalized real‐time substructure testing with implicit integration

This paper presents a new method, called the equivalent force control method, for solving the nonlinear equations of motion in a real‐time substructure test using an implicit time integration algorithm. The method replaces the numerical iteration in implicit integration with a force‐feedback control loop, while displacement control is retained to control the motion of an actuator. The method is formulated in such a way that it represents a unified approach that also encompasses the effective force test method. The accuracy and effectiveness of the method have been demonstrated with numerical simulations of real‐time substructure tests with physical substructures represented by spring and damper elements, respectively. The method has also been validated with actual tests in which a Magnetorheological damper was used as the physical substructure. Copyright © 2007 John Wiley & Sons, Ltd.     

Expected loss‐based alarm threshold set for earthquake early warning systems

Earthquake early warning systems (EEWS) seem to have potential as tools for real‐time seismic risk management and mitigation. In fact, although the evacuation of buildings requires warning time not available in many urbanized areas threatened by seismic hazard, they may still be used for the real‐time protection of critical facilities using automatic systems in order to reduce the losses subsequent to a catastrophic event. This is possible due to the real‐time seismology, which consists of methods and procedures for the rapid estimation of earthquake features, as magnitude and location, based on measurements made on the first seconds of the P‐waves. An earthquake engineering application of earthquake early warning (EEW) may be intended as a system able to issue the alarm, if some recorded parameter exceeds a given threshold, to activate risk mitigation actions before the quake strikes at a site of interest. Feasibility analysis and design of such EEWS require the assessment of the expected loss reduction due to the security action and set of the alarm threshold. In this paper a procedure to carry out these tasks in the performance‐based earthquake engineering probabilistic framework is proposed. A merely illustrative example refers to a simple structure assumed to be a classroom. Structural damage and non‐structural collapses are considered; the security action is to shelter occupants below the desks. The cost due to a false alarm is assumed to be related to the interruption of didactic activities. Results show how the comparison of the expected losses, for the alarm‐issuance and non‐issuance cases, allows setting the alarm threshold on a quantitative and consistent basis, and how it may be a tool for the design of engineering applications of EEW. Copyright © 2007 John Wiley & Sons, Ltd.     

Analysis of the weightiness of site effects on reinforced concrete (RC) building seismic behaviour: The Adra town example (SE Spain)

The damage distribution in Adra town (south‐eastern Spain) during the 1993 and 1994 Adra earthquakes (5.0 magnitude), that reached a maximum intensity degree of VII (European Macroseismic Scale (EMS scale)), was concentrated mainly in the south‐east zone of the town and the most relevant damage occurred in reinforced concrete (RC) buildings with four or five storeys. In order to evaluate the influence of ground condition on RC building behaviour, geological, geomorphological and geophysical surveys were carried out, and a detailed map of ground surface structure was obtained. Short‐period microtremor observations were performed in 160 sites on a 100m × 100m dimension grid and Nakamura's method was applied in order to determine a distribution map of soil predominant periods. Shorter predominant periods (0.1–0.3 s) were found in mountainous and neighbouring zones and larger periods (greater than 0.5 s) in thicker Holocene alluvial fans. A relationship T = (0.049 ± 0.001)N, where T is the natural period of swaying motion and N is the number of storeys, has been empirically obtained by using microtremor measurements at the top of 38 RC buildings (ranging from 2 to 9 storeys). 1‐D simulation of strong motion on different soil conditions and for several typical RC buildings were computed, using the acceleration record in Adra town of the 1993 earthquake. It is noteworthy that all the aforementioned results show the influence of site effects in the degree and distribution of observed building damage. Copyright © 2007 John Wiley & Sons, Ltd.     

Stability and accuracy analysis of outer loop dynamics in real‐time pseudodynamic testing of SDOF systems

Real‐time pseudodynamic (PSD) testing is an experimental technique for evaluating the dynamic behaviour of a complex structure. During the test, when the targeted command displacements are not achieved by the test structure, or a delay in the measured restoring forces from the test structure exists, the reliability of the testing method is impaired. The stability and accuracy of real‐time PSD testing in the presence of amplitude error and a time delay in the restoring force is presented. Systems consisting of an elastic single degree of freedom (SDOF) structure with load‐rate independent and dependent restoring forces are considered. Bode plots are used to assess the effects of amplitude error and a time delay on the steady‐state accuracy of the system. A method called the pseudodelay technique is used to derive the exact solution to the delay differential equation for the critical time delay that causes instability of the system. The solution is expressed in terms of the test structure parameters (mass, damping, stiffness). An error in the restoring force amplitude is shown to degrade the accuracy of a real‐time PSD test but not destabilize the system, while a time delay can lead to instability. Example calculations are performed for determining the critical time delay, and numerical simulations with both a constant delay and variable delay in the restoring force are shown to agree well with the stability limit for the system based on the critical time delay solution. The simulation models are also used to investigate the effects of a time delay in the PSD test of an inelastic SDOF system. The effect of energy dissipation in an inelastic structure increases the limit for the critical time delay, due to the energy removed from the system by the energy dissipation. Copyright © 2007 John Wiley & Sons, Ltd.     

Effectiveness of simple approaches in mitigating residual deformations in buildings

Developments in performance‐based seismic design and assessment approaches have emphasized the importance of considering residual deformations. Recent investigations have also led to a proposed direct displacement‐based design (DDBD) approach which includes an explicit consideration of the expected residual deformations as an integral part of the design process. Having estimated the expected residual deformations in a structure, engineers are faced with the problem of reducing them to meet the targeted performance levels under pre‐defined seismic hazard levels. Previous studies have identified the post‐yield stiffness as a primary factor influencing the magnitude of residual deformations in single degree of freedom and multiple degree of freedom structures. In this paper, a series of simple approaches to increase the post‐yield stiffness of traditional framed and braced systems for the purpose of reducing residual deformations are investigated. These methods do not utilize recentring post‐tensioned technology. This contribution addresses the feasibility of altering the lateral post‐yield stiffness of structural systems by: (i) using different reinforcement materials with beneficial features in their stress–strain behaviour; (ii) re‐designing the section geometry and properties of primary seismic‐resisting elements; and (iii) introducing a secondary elastic frame to act in parallel with the primary system. The efficiency of each of these techniques is investigated through monotonic and cyclic moment‐curvature and non‐linear time‐history analyses. Of these approaches the design and introduction of an elastic secondary system was found to be most effective and consistent in reducing residual deformations. A simplified design approach for achieving the desired increase of a system's post‐yield stiffness is also presented. Copyright © 2007 John Wiley & Sons, Ltd.     

Earthquake‐resistant structural design through energy demand and capacity

An energy‐based earthquake‐resistant structural design method is proposed. The proposed method uses specific input energy spectra, modal or time‐history analyses, and energy distribution among structural members. For a given member strength and stiffness, a relationship between the energy attributable to damage absorbed by a member and its cumulative ductility demand can be determined. Member strength, stiffness and energy capacity are design parameters which are simultaneously used in the design. The method can avoid soft‐storey design. The damage is measured based on a cumulative basis considering earthquake magnitude, frequency, and duration. Tests have been carried out to determine energy absorbing capacities of various structural components. More efforts are needed to make the energy‐based earthquake‐resistant structural design practical, but ssimple formulations for this method are possible. Copyright © 2007 John Wiley & Sons, Ltd.     

Acceleration response modification factors for nonstructural components attached to inelastic moment‐resisting frame structures

A statistical analysis of the peak acceleration demands for nonstructural components (NSCs) supported on a variety of stiff and flexible inelastic regular moment‐resisting frame structures with periods from 0.3 to 3.0 s exposed to 40 far‐field ground motions is presented. Peak component acceleration (PCA) demands were quantified based on the floor response spectrum (FRS) method without considering dynamic interaction effects. This study evaluated the main factors that influence the amplification or decrease of FRS values caused by inelasticity in the primary structure in three distinct spectral regions namely long‐period, fundamental‐period, and short‐period region. The amplification or decrease of peak elastic acceleration demands depends on the location of the NSC in the supporting structure, periods of the component and building, damping ratio of the component, and level of inelasticity of the supporting structure. While FRS values at the initial modal periods of the supporting structure are reduced due to inelastic action in the primary structure, the region between the modal periods experiences an increase in PCA demands. A parameter denoted as acceleration response modification factor (R_acc) was proposed to quantify this reduction/increase in PCA demands. Copyright © 2007 John Wiley & Sons, Ltd.     

Improved seismic hazard model with application to probabilistic seismic demand analysis

An improved seismic hazard model for use in performance‐based earthquake engineering is presented. The model is an improved approximation from the so‐called ‘power law’ model, which is linear in log–log space. The mathematics of the model and uncertainty incorporation is briefly discussed. Various means of fitting the approximation to hazard data derived from probabilistic seismic hazard analysis are discussed, including the limitations of the model. Based on these ‘exact’ hazard data for major centres in New Zealand, the parameters for the proposed model are calibrated. To illustrate the significance of the proposed model, a performance‐based assessment is conducted on a typical bridge, via probabilistic seismic demand analysis. The new hazard model is compared to the current power law relationship to illustrate its effects on the risk assessment. The propagation of epistemic uncertainty in the seismic hazard is also considered. To allow further use of the model in conceptual calculations, a semi‐analytical method is proposed to calculate the demand hazard in closed form. For the case study shown, the resulting semi‐analytical closed form solution is shown to be significantly more accurate than the analytical closed‐form solution using the power law hazard model, capturing the ‘exact’ numerical integration solution to within 7% accuracy over the entire range of exceedance rate. Copyright © 2007 John Wiley & Sons, Ltd.     

Prediction of the median IDA curve by employing a limited number of ground motion records

A methodology has been proposed which can be used to reduce the number of ground motion records needed for the reliable prediction of the median seismic response of structures by means of incremental dynamic analysis (IDA). This methodology is presently limited to predictions of the median IDA curve only. The reduction in the number of ground motion records needed to predict the median IDA curve is achieved by introducing a precedence list of ground motion records. The determination of such a list is an optimization problem, which is solved in the paper by means of (1) a genetic algorithm and (2) a proposed simple procedure. The seismic response of a simple, computationally non‐demanding structural model has been used as input data for the optimization problem. The presented example is a three‐storey‐reinforced concrete building, subjected to two sets of ground motion records, one a free‐field set and the other a near‐field set. It is shown that the median IDA curves can be predicted with acceptable accuracy by employing only four ground motion records instead of the 24 or 30 records, which are the total number of ground motion records for the free‐field and near‐field sets, respectively. Copyright © 2007 John Wiley & Sons, Ltd.     

Effects of surface characteristics on infiltration patterns in an arid shrub desert

Precipitation is often the sole source of water replenishment in arid and semi‐arid areas and, thus, plays a pertinent role in sustaining desert ecosystems. Revegetation over 40 years using mainly Artemisia ordosica and Caragana korshinskii at Shapotou Desert Experimental Research Station near Lanzhou, China, has established a dwarf‐shrub and microbiotic soil crust cover on the stabilized sand dunes. The redistribution of infiltrated moisture through percolation, root extraction, and evapotranspiration pathways was investigated. Three sets of time‐domain reflectometry (TDR) probes were inserted horizontally at 5, 10, 15, 20, 30 and 40 cm depths below the ground surface in a soil pit. The three sets of TDR probes were installed in dwarf‐shrub sites of A. ordosica and C. korshinskii community with and without a microbiotic soil crust cover, and an additional set was placed in a bare sand dune area that had neither vegetation nor a microbiotic soil crust present. Volumetric soil moisture content was recorded at hourly intervals and used in the assessment of infiltration for the different surface covers. Infiltration varied greatly, from 7·5 cm to more than 45 cm, depending upon rainfall quantity and soil surface conditions. In the shrub community area without microbiotic soil crust cover, infiltration increased due to preferential flow associated with root tunnels. The microbiotic soil crust cover had a significant negative influence on the infiltration for small rainfall events (～10 mm), restricting the infiltration depth to less than 20 cm and increasing soil moisture content just beneath the soil profile of 10 cm, whereas it was not as strong or clear for larger rainfall events (～60 mm). For small rainfall events, the wetting front depth for the three kinds of surface cover was as follows: shrub community without microbiotic soil crust > bare area > shrub community with microbiotic soil crust. In contrast, for large rainfall events, infiltration was similar in shrub communities with and without microbiotic soil crust cover, but significantly higher than measured in the bare area. Soil water extraction by roots associated with evapotranspiration restricted the wetting front penetration after 1 to 3 h of rainfall. Copyright © 2006 John Wiley & Sons, Ltd.