Numerical Evaluation of Seismic Soil Pressures on Rigid Walls Fixed to the Bedrock

Seismic soil pressures developed on a 7 m rigid retaining wall fixed to the bedrock are investigated using a finite element model that engages nonlinear soil intended materials available in OpenSees. This allows incorporation of the inelastic behavior of the soil and wave propagation effects in the soil-wall system seismic response. The nonlinear response of the soil was validated using the well-stablished, frequency-domain, linear-equivalent approach. An incremental dynamic analysis was implemented to comprehensively examine the effect of soil nonlinearity and input motion on the induced seismic pressures and to evaluate current code equations/methodologies at different levels of earthquake intensity. The results show that soil nonlinearity and seismic wave amplification may play an important role in the response of the soil-wall system. Therefore, methodologies that rely only on peak ground acceleration may introduce large bias on the estimated seismic pressures in scenarios where high nonlinearity and site amplification are expected.     

IDENTIFICATION AND VERIFICATION OF SEISMIC DEMAND FROM DIFFERENT HYSTERETIC MODELS

The goal of this paper is to develop a modified Bouc-Wen hysteretic model from cyclic loading test data for reinforced columns, including the behavior of stiffness degradation, strength deterioration, pinching and softening effects of RC members. Seismic demands on this inelastic single degree of freedom system when subjected to both near-fault ground motion and far-field ground motion excitations were examined.

The cyclic loading test of reinforced concrete columns was experimentally observed and a system identification computer program was developed to solve each control parameter of the hysteretic model. A least-squared method for identifying parameters of the model is proposed in this paper. The hysteretic constitutive law produces a smoothly varying hysteresis such as the control-parameters for strength deterioration, stiffness degradation, pinching and softening effects. Two implementations of (1) flexure damage and (2) shear damage were conducted to provide better understanding of hysteretic behavior of RC structural members. A pseudo-dynamic experiment was also developed to verify the model parameters.

Based on the developed hysteretic model, the seismic demand of this inelastic model was investigated by using both near-fault ground motion data and far-field ground motion data as input motion. An R-μT inelastic response spectrum from different hysteretic models was generated.     

A NON-PARAMETRIC APPROACH TO ATTENUATION RELATIONS

Abstract

A non-parametric multidimensional regression method is proposed for the prediction of seismic ground motion parameters. The main features which distinguish the method from standard regression procedures are: (1) The relationship between the input and output variables is not selected a priori by a prediction law, (2) an arbitrary number of input variables Can be taken into account, provided that an appropriate data base exists, and (3) the computational procedure is very simple. The results can be easily updated when new information becomes available. The method has been applied for the derivation of attenuation relations by using a combination of databases compiled by other researchers. In the majority of the cases discussed in this paper, the method was used for the prediction of horizontal peak ground acceleration as a function of magnitude and distance. In some cases, ground conditions were also taken into account. Some results on the attenuation relations of peak ground velocity and displacement, as well as Arias intensity, are also presented.     

Estimation of Strong Ground Motions in Southwest Western Australia with a Combined Green's Function and Stochastic Approach

As only a very limited number of earthquake strong ground motion records are available in southwest Western Australia (SWWA), it is difficult to derive a reliable and unbiased strong ground motion attenuation model based on these data. To overcome this, in this study a combined approach is used to simulate ground motions. First, the stochastic approach is used to simulate ground motion time histories at various epicentral distances from small earthquake events. Then, the Green's function method, with the stochastically simulated time histories as input, is used to generate large event ground motion time histories. Comparing the Fourier spectra of the simulated motions with the recorded motions of a ML6.2 event in Cadoux in June 1979 and a ML5.5 event in Meckering in January 1990, provides good evidence in support of this method. This approach is then used to simulate a series of ground motion time histories from earthquakes of varying magnitudes and distances. From the regression analyses of these simulated data, the attenuation relations of peak ground acceleration (PGA), peak ground velocity (PGV), and response spectrum of ground motions on rock site in SWWA are derived.     

EVALUATION OF CAPACITY SPECTRUM METHOD FOR ESTIMATING THE PEAK INELASTIC RESPONSES

The accuracy of the capacity spectrum method (CSM) depends on the precise estimation of equivalent period and damping ratio as well as the modification of the demand spectrum. In this paper, the CSM provided in ATC-40 for estimating the peak inelastic responses is evaluated. First, the effect of equivalent period and damping ratio estimation on the accuracy of the CSM is assessed. Analyses results indicate that the difference between estimation methods is large when the structural nonlinearity is large, but becomes negligible as the hardening ratio increases. Next, the reduction factors provided in ATC-40 and Eurocode are evaluated. It is found that the acceleration responses obtained using the factor of Eurocode is closer to the actual ones than those obtained using the factors of ATC-40. Finally, the demand spectrum is constructed using the peak absolute acceleration and pseudo-acceleration. The results obtained using the peak absolute acceleration is found to be generally larger than those obtained using the pseudo-ones. Since the original CSM generally underestimates the response, the use of peak absolute acceleration in the construction of demand spectrum produces the response relatively closer to the exact one. However, the use of peak absolute acceleration overestimates the response more when the original CSM overestimates the response.     

SEISMIC HAZARD ASSESSMENT IN FLORENCE CITY ITALY

A seismic hazard analysis of Florence city was performed in the frame of a project concerning the dynamic behaviour of cable-stayed bridges. Both a probabilistic approach and a methodology based on the use of a local macroseismic catalogue were applied. A local catalogue was expressly compiled for this purpose, to collect the macroseismic intensities actually observed at the site as a result of past earthquakes. This sort of catalogue is an independent tool to verify the assumptions of the probabilistic approach (seismic zoning, earthquake recurrence relation, attenuation model), though it can supply results in terms of macroseismic intensity only and reflects the effective seismic history at the site, without taking into account any variability. The Cornell' methodology was used to assess probabilistic hazard in terms of macroseismic intensity, peak ground acceleration, peak ground velocity, and pseudovelocity uniform response spectra. The local catalogue points out level VII of the Mercalli-Cancani-Sieberg scale (MCS) as the maximum intensity historically observed in Florence. The probabilistic approach leads to the consideration of intensity VIII MCS as the maximum credible for the city. The probabilistic analysis in terms of ground motion was performed using attenuation relations estimated for alluvium sites, since the geology of Florence area is represented by fluvial and lacustrine deposits of various thickness. Peak ground acceleration values with 90% non exceedence probability in 50 and 500 years are respectively 145 and 219 cm/s's for a shallow alluvium site, and 95 and 157 cm/s's for a deep alluvium site; the corresponding peak ground velocity values for sites located on alluvium are 6.41 and 11.76 cm/s. Uniform response spectra are provided for shallow and deep alluvium sites, according to frequency-dependent attenuation relations estimated from strong Italian earthquakes.     

Non-Iterative Equivalent Linearization Based on Secant Period for Estimating Maximum Deformations of Existing Structures

The capacity spectrum method of ATC-40 uses the secant period as the equivalent period of equivalent linear systems. Therefore, it results in a direct graphical comparison. The maximum inelastic displacement and acceleration demands of structures can be simultaneously obtained from the intersection of the demand and capacity diagrams. However, for evaluation of existing structures, the demands need to be determined through iterations since the equivalent period and damping of the equivalent linear systems currently available are both a function of the (displacement) ductility ratio, which is unknown and is the target of evaluation. In addition, the equivalent damping used in the capacity spectrum method is independent of periods of vibration. It may lead to poor estimations of maximum responses especially for short-period systems. This article proposes two equivalent linear systems based on the secant period to estimate the maximum displacement and acceleration responses of existing structures. Both the recommended equivalent period and damping are defined by the strength ratio (elastic lateral strength/yield lateral strength), rather than the ductility ratio. Because the strength ratio of existing structures is a known parameter, the maximum displacement and acceleration responses of these structures can be determined without iterations. Besides, effects of periods of vibration on the equivalent linear systems are also included in this study. The equivalent damping is derived from statistical analyses for bilinear single-degree-of-freedom (SDOF) systems with different periods of vibration, strength ratios and post-yield stiffness based on 72 earthquake ground motions recorded on firm sites. Procedures and examples for applications of the proposed equivalent linear systems on nonlinear static analysis procedures are also provided.     

Floor Spectra of Inelastic RC Frame Buildings Considering Ground Motion Characteristics

A range of reinforced concrete frame buildings with different levels of inelasticity as well as periods of vibration is analyzed to study the floor response. The derived floor acceleration response spectra are normalized by peak ground acceleration, peak floor acceleration, and ground response spectrum. The normalization with respect to ground response spectrum leads to the lowest coefficients of variation. Based on this observation as well as previous studies, an amplification function is proposed that can be used to develop design floor spectra from the ground motion spectrum, considering the building’s dynamic characteristics and level of inelasticity.     

Maximum Possible Ground Motion for Linear Structures

This article develops a method to generate ground motion time histories that maximize the response of a given linearly elastic structure. The root mean square (RMS) level of the input power spectral density (PSD) is used as a strong motion parameter. It is related to seismological data that is readily available. An empirical relation to estimate RMS value of the PSD from peak ground acceleration, magnitude, rupture distance, and shear wave velocity is derived from world-wide strong motion data. The ground motion is obtained by solving the inverse problem such that the structural response is maximized under the constraint of fixed value of RMS level of the input PSD enforced using a Lagrange multiplier. The proposed methodology is illustrated for a single-degree of freedom system, a six storey building and an earthen dam. It is shown that the critical PSD obtained in all the cases is a narrow band process resulting in stochastic resonance and not a Dirac-delta function with the entire energy of the system concentrated at its natural frequency. Moreover, the critical excitation samples generated using this critical PSD resembles actual earthquake acceleration time histories.     

A PROCEDURE FOR COMBINING VERTICAL AND HORIZONTAL SEISMIC ACTION EFFECTS

The vertical component of earthquake ground motion has generally been neglected in the earthquake-resistant design of structures. This is gradually changing due to the increase in near-source records obtained recently, coupled with field observations confirming the possible destructive effect of high vertical vibrations.

In this paper, simple procedures are suggested for assessing the significance of vertical ground motion, indicating when it should be included in the determination of seismic actions on buildings. Proposals are made for the calculation of elastic and inelastic vertical periods of vibration incorporating the effects of vertical and horizontal motion amplitude and the cross-coupling between the two vibration periods. Simplified analysis may then be used to evaluate realistic vertical forces by employing the vertical period of vibration with pertinent spectra without resorting to inelastic dynamic analysis.

Finally, a procedure is suggested for combining vertical and horizontal seismic action effects which accounts for the likelihood of coincidence, or otherwise, of peak response in the two directions.     

SIGNIFICANT GROUND MOTION PARAMETERS FOR EVALUATION OF THE SEISMIC PERFORMANCE OF SLENDER MASONRY TOWERS

The evaluation of seismic risk of masonry monuments requires to study the combination of vulnerability and hazard. In the present work, the global seismic response of slender masonry towers has been studied by means of a specific 3-D fibre model. Accounting for the particular behaviour of such structures, the hazard should also be described by a suitable measure of intensity of the seismic action. A variety of different parameters relating with the ground acceleration recordings have been investigated for what regards their correlation with the damage indicators of the model. The combination of the peak ground velocity of the horizontal component and of the significant duration is an effective measure of intensity. This measure can be improved by considering the accord of the frequency content of the ground motion with the dynamical characteristics of the tower. Since in some cases the effect of the vertical component proved to be important, a further improvement can be obtained by taking into account also the vertical ground motion intensity.     

Shaking Table Testing of a Full-Scale Prefabricated Three-Story Timber-Frame Building

Shake table tests were carried out on a 7 m × 5 m three-story, timber light-frame building (7.5 m height) at the TreesLab laboratory (Eucentre) in Pavia. The aim of the research was to evaluate the seismic behavior of a typical Italian prefabricated timber building and to study the interaction between the individual structural components tested in quasi-static manner in a previous experimental study. The 1979 Montenegro Earthquake ground motion, recorded at the Ulcinj-Hotel Albatros station, was selected as the ground motion for seismic tests. The maximum peak ground acceleration was scaled to 0.07 g, 0.27 g, 0.5 g. 0.7 g, and 1 g in order to evaluate the building’s performance at different levels of seismic input. More than 100 instruments were used to monitor the behavior of the building during seismic tests measuring acceleration, displacement, and forces. The visual inspection shows that the building did not show any damage during all seismic tests. However the data analysis (dynamic identification, capacity spectrum, inter-story drift) confirm that during the 1.00 g test the structure went beyond its linear elastic limit. The results obtained from this experimental study suggest that the design hypotheses commonly adopted in practice for seismic analysis (e.g., in terms of force distributions between the walls, and also the behavior factor q) are not always consistent with the real behavior of timber frame multi-story buildings, and should be backed by more accurate knowledge of the contributions of the individual structural components.     

Inelastic Displacement Demand of Strength-Degraded Structures

This article presents findings from parametric studies involving nonlinear time-history analyses of inelastic systems with and without strength degradation. Results showed that estimates based on the equal-displacement and equal-energy propositions can be exceeded significantly by the inelastic displacement demands in the acceleration and velocity-sensitive regions of the response spectrum. The displacement demand behaviour is sensitive to the strength degradation and the frequency properties of the ground shaking. With a modest strength reduction factor of 2, the inelastic displacement demand would typically be constrained by the Peak Displacement Demand as indicated on the elastic displacement response spectrum for 5% damping.     

MODIFICATION OF SITE RESPONSE IN PRESENCE OF LOCALISED SOFT LAYER

Experimental and numerical simulations are performed to evaluate the modification of ground response resulting from either the presence of soft layers or occurrence of partial liquefaction. Results from two densely instrumented dynamic centrifuge tests are presented to show the ambiguous role played by the presence of a soft layer. It was found that the lateral extent of the soft layer has significant influence on the overall response of the layered strata and any structure founded on it. The experimental observations are supported by simplified numerical analysis. The amplification or deamplification of the input motion is found to be a function of the ratio of the width of soft layer to the wave length. Based on the numerical analysis, a general function describing the site amplification is presented which may be used as a guide in seismic design of foundations in such layered strata.     

AN EVALUATION OF THE STRONG GROUND MOTION RECORDED DURING THE MAY 1, 2003 BINGÖL TURKEY,EARTHQUAKE

An important record of ground motion from a M6.4 earthquake occurring on May 1, 2003, at epicentral and fault distances of about 12 and 9 km, respectively, was obtained at a station near the city of Bingöl, Turkey. The maximum peak ground values of 0.55 g and 36 cm/s are among the largest ground-motion amplitudes recorded in Turkey. From simulations and comparisons with ground motions from other earthquakes of comparable magnitude, we conclude that the ground motion over a range of frequencies is unusually high. Site response may be responsible for the elevated ground motion, as suggested from analysis of numerous aftershock recordings from the same station. The mainshock motions have some interesting seismological features, including ramps between the P-and S-wave that are probably due to near- and intermediate-field elastic motions and strong polarisation oriented at about 39 degrees to the fault (and therefore not in the fault-normal direction). Simulations of motions from an extended rupture explain these features. The N10E component shows a high-amplitude spectral acceleration at a period of 0.15 seconds resulting in a site specific design spectrum that significantly overestimates the actual strength and displacement demands of the record. The pulse signal in the N10E component affects the inelastic spectral displacement and increases the inelastic displacement demand with respect to elastic demand for very long periods.     

INSTANTANEOUS OPTIMAL CONTROL FOR SEISMIC ANALYSIS OF NON-LINEAR STRUCTURES

In this paper, an effective active control algorithm is developed for the vibration control of non-linear structural systems subjected to earthquake excitation. It is an attempt to include the non-linear characteristics of the structural behaviour throughout the entire analysis (design and validation), accounting for the eventual cumulative structural damage and energy dissipation. This is a very important factor since, in current design practice, structures are assumed to behave nonlinearly when subjected to strong ground motion. The proposed algorithm focusses on the instantaneous optimal control approach for the development of the control algorithm where the nonlinearities are brought into the analysis through a non-linear state vector and a non-linear open loop term. A performance index that is quadratic in the control force and in the non-linear states and is subjected to a non-linear constraint equation, is minimised at every time step. The effectiveness of the proposed non-linear instantaneous optimal control (NIOC) strategy is critically evaluated in comparison with currently available active control techniques. Numerical simulation results indicate that the proposed approach provides a significant reduction of the peak response quantities, such as maximum response deformation, maximum response acceleration, ductility of the system, associated with a reduced maximum control force.     

ANALYTICAL MODELS FOR BRICK MASONRY INFILLED R/C FRAMES UNDER LATERAL LOADING

Proposed in this paper are two analytical models for predicting the inelastic response of unreinforced brick masonry infills in reinforced concrete frames subjected to mono-tonic and reversed cyclic loading. The first model is based on the traditional diagonal strut concept, while the second one is a simple isoparametric element with shear deformation only. All the essential characteristics of the hysteretic behaviour of the panel, including strength and stiffness degradation, pinching and slippage, are explicitly taken into account. The models are implemented in a general-purpose program for the inelastic time-history analysis of structures, and are used for studying the seismic behaviour of typical multistorey frames with various arrangements of infill panels, including structures with an open ground storey. The results of the analysis are in agreement with both experimentally observed behaviour and with experience regarding seismically damaged buildings.     

NONLINEAR SHAKE TABLE IDENTIFICATION AND CONTROL FOR NEAR-FIELD EARTHQUAKE TESTING

The primary focus of a structural shake table system is the accurate reproduction of acceleration records for testing. However, many systems deliver variable and less than optimal performance, particularly when reproducing large near-field seismic events that require extreme table performance. Improved identification and control methods are developed for large hydraulic servo-actuated shake table systems that can exhibit unacceptable tracking response for large, near-field seismic testing. The research is presented in the context of a 5-tonne shake table facility at the University of Canterbury that is of typical design. The system is identified using a frequency response approach that accounts for the actual magnitudes and frequencies of motion encountered in seismic testing. The models and methods developed are experimentally verified and the impact of different feedback variables such as acceleration, velocity and displacement are examined.

The methods show that shake table control in testing large near-field seismic events is often a trade off between accurate tracking and nonlinear velocity saturation of the hydraulic valves that can result in severe acceleration spikes. Control methods are developed to improve performance and include both acceleration and displacement feedback to reduce the acceleration spikes, and record modification, where the reference signal is modified to conform to the shake table's operational parameters. Results show record modification gives exact tracking for near-field ground motions, and optimal system response for reference signals with velocity components greater then the system capabilities. Overall, the research presents a methodology for simple effective identification, modelling, diagnosis and control of structural shake table systems that can be readily generalised and applied to any similar facility.     

Design Earthquake Ground Motion Prediction for Perth Metropolitan Area with Microtremor Measurements for Site Characterization

Perth is the largest city in Western Australia and home to three-quarters of the state's residents. In recent decades, there have been a lot of earthquake activities just east of Perth in an area known as the South-West Seismic Zone. Previous numerical results of site response analyses based on limited available geology information for PMA indicated that Perth Basin might amplify the bedrock motion by more than 10 times at some frequencies and at some sites. Hence, more detailed studies on site characterization and amplification are necessary. The microtremor method using spatial autocorrelation (SPAC) processing is a useful tool for gaining thickness and shear wave velocity (SWV) of sediments and has been adopted in many previous studies. In this study, the response spectrum of rock site corresponding to the 475-year return period for PMA is defined according to the probabilistic seismic hazard analysis (PSHA) based on the latest ground motion attenuation model of Southwest Western Australia. Site characterization in PMA is performed using two microtremor measurements, namely SPAC technique and H/V method. The clonal selection algorithm (CSA) is introduced to perform direct inversion of SPAC curves to determine the soil profiles of representative PMA sites investigated in this study. Using the simulated bedrock motion as input, the responses of the soil sites are estimated using numerical method based on the shear-wave velocity vs. depth profiles determined from the SPAC technique. The response spectrum of the earthquake ground motion on surface of each site is derived from the numerical results of the site response analysis, and compared with the respective design spectrum defined in the Australian Earthquake Loading Code. The comparison shows that the code spectra are conservative in the short period range, but may slightly underestimate the response spectrum at some long period range.     

A THEORETICAL ATTENUATION MODEL FOR EARTHQUAKE-INDUCED GROUND MOTION

A theoretical attenuation model of earthquake-induced ground motion is presented and discussed. This model is related directly to physical quantities such as source and wave motion parameters. An attenuation formula for rms acceleration of ground motion is derived and verified using acceleration data from moderate-sized earthquakes recorded in Iceland from 1986 to 1997. The source parameters and the crustal attenuation are computed uniformly for the applied earthquake data. Furthermore, attenuation formulas for peak ground acceleration are put forward.