Evaluation of Peak Ground Velocity as a “Good” Intensity Measure for Near-Source Ground Motions

This paper investigates the “goodness” of peak ground velocity as a dependable intensity measure for the earthquake shaking of civil structures. The paper stresses the importance of distinguishing between acceleration pulses and velocity pulses, and identifies two classes of near-source ground motions: those where the peak ground velocity is the integral of a distinguishable acceleration pulse and those where the peak ground velocity is the result of a succession of high-frequency, one-sided acceleration spikes. It is shown that the shaking induced by the former class is in general much more violent than the shaking induced by the latter class of motions even when motions that belong to the former class may be generated by significantly smaller-magnitude earthquakes. Building on the dimensional analysis introduced in the companion papers this paper shows that both linear and nonlinear structural responses from a variety of records which exhibit distinguishable pulses scale better with the peak pulse acceleration than with the peak pulse velocity, indicating that the peak pulse acceleration is a more representative intensity measure of the earthquake shaking. This conclusion is further supported from the response analysis of linear and bilinear single-degree-of-freedom oscillators subjected to selected records from the 1999 Chi-Chi Taiwan earthquake that exhibit unusually high and long period velocity pulses. The paper shows that these high velocity pulses alone do not impose unusual demands on most civil structures. What is more detrimental are local, distinguishable acceleration pulses that override the long period velocity pulses.     

Pseudostatic Approach for Seismic Analysis of Single Piles

This paper evaluates a simple approximate methodology for estimating the maximum internal forces of piles subjected to lateral seismic excitation. The method involves two main steps: computation of the free-field soil movements caused by the earthquake and the analysis of the response of the pile to the maximum free-field soil movements (considered as static movements) plus a static loading at the pile head, which depends on the computed spectral acceleration of the structure being supported. The applicability of this approach has been verified by an independent benchmark analysis developed by the writers. It is demonstrated that the proposed method yields reasonable estimates of the pile maximum moment and shear. The methodology is then used to obtain the response of the Ohba-Ohashi bridge in Japan to one of the earthquakes that occurred in the 1980s. Good agreement is found between the computed and measured pile moments.     

Empirical Predictive Models for Earthquake-Induced Sliding Displacements of Slopes

Earthquake-induced sliding displacement is the parameter most often used to assess the seismic stability of slopes. The expected displacement can be predicted as a function of the characteristics of the slope (yield acceleration) and the ground motion (e.g., peak ground acceleration), yet there is significant aleatory variability associated with the displacement prediction. Using multiple ground motion parameters to characterize the earthquake shaking can significantly reduce the variability in the prediction. Empirical predictive models for rigid block sliding displacements are developed using displacements calculated from over 2,000 acceleration–time histories and four values of yield acceleration. These empirical models consider various single ground motion parameters and vectors of ground motion parameters to predict the sliding displacement, with the goal of minimizing the standard deviation of the displacement prediction. The combination of peak ground acceleration and peak ground velocity is the two parameter vector that results in the smallest standard deviation in the displacement prediction, whereas the three parameter combination of peak ground acceleration, peak ground velocity, and Arias intensity further reduces the standard deviation. The developed displacement predictive models can be used in probabilistic seismic hazard analysis for sliding displacement or used as predictive tools for deterministic earthquake scenarios.     

Earthquake Simulator Testing and Seismic Evaluation of Suspended Ceilings

Past seismic events have demonstrated the vulnerability of suspended ceilings, classified as nonstructural components, to earthquake damage. Engineers, architects, and manufacturers all participate to help ensure that the units perform satisfactorily during earthquakes. To address the seismic susceptibility, the U.S. national codes and federal and regulatory guidelines recommend two distinct approaches. The first mandates capacity and installation requirements and the second addresses damage states, introducing performance based. Application of these methods entails either accurate structural analysis or seismic qualification by experimentation. Until recently, scant data were available to evaluate the adequacy of the required installations or the damage states. To address this issue, full-scale earthquake laboratory tests of suspended ceiling systems have been undertaken by researchers and manufacturers. Experimentation showed that ceilings meeting the code requirements performed well. The only mode of failure observed during the tests was the loss of panels. Laboratory data were also used to construct fragility curves for the specimen. Large vertical accelerations, typically not observed in the field, dislodged panels close to the center of the test frame. This type of response differs from that presented in past earthquake reconnaissance reports and thus necessitates further examination.     

Investigation of effect of near-fault motions on substandard bridge structures

The objective of this study was to investigate the effects of near-fault ground motions on substandard bridge columns and piers. To accomplish these goals, several large scale reinforced concrete models were constructed and tested on a shake table using near- and far-field ground motion records. Because the input earthquakes for the test models had different characteristics, three different measures were used to evaluate the effect of the input earthquake. These measures are peak shake table acceleration, spectral acceleration at the fundamental period of the test specimens, and the specimen drift ratios.For each measure, force-displacement relationships, strains, curvatures, drift ratios, and visual damage were evaluated.Results showed that regardless of the measure of input or response, the near-fault record generally led to larger strains,curvatures, and drift ratios. Furthermore, residual displacements were small compared to those for columns meeting current seismic code requirements.     

Seasonally Frozen Soil Effects on the Seismic Site Response

Several large-magnitude earthquakes, including the Prince William Sound earthquake of March 1964 and the Denali earthquake of November 2002, occurred in the state of Alaska and caused considerable damages to its transportation system, including damage to several highway bridges and related infrastructure. Some of these damages are related to frozen soil effects. However, only limited research has been carried out to investigate the effects of frozen soils on seismic site responses. A systematic investigation of seasonally frozen soil effects on the seismic site response has been conducted and is presented in this paper. One bridge site in Anchorage, Alaska, was selected to represent typical sites with seasonally frozen soils. A set of input ground motions was selected from available strong-motion databases and scaled to generate an ensemble of hazard-consistent input motions. One-dimensional equivalent linear analysis was adopted to analyze the seismic site response for three seismic hazard levels, i.e., maximum considered earthquake (MCE), AASHTO design, and service design level hazards. Parametric studies were conducted to assess the sensitivity of the results to uncertainties associated with the thickness and shear-wave velocity of seasonally frozen soils. The results show that the spectral response of ground motions decreases as the thickness of seasonally frozen soil increases, and the results are insensitive to the shear-wave velocity of seasonally frozen soils. In conclusion, it is generally conservative to ignore the effects of seasonally frozen soils on seismic site response in the design of highway bridges.     

Dynamic Testing of a Masonry Structure on a Passive Isolation System

Lightly reinforced and unreinforced masonry buildings have not performed well in earthquakes. Evaluation of past performance of masonry structures has led to more stringent design and construction requirements in the current building codes, and has raised concerns about the performance of existing lightly reinforced and unreinforced masonry buildings in future earthquakes. Base isolation has been shown to be effective in reducing damage to large building structures, and appears to be particularly effective in protecting stiff masonry structures. Using the base isolation principle, Kansas State University’s stiffness decoupler for the base isolation of structures (SDBIS) was designed to effectively reduce the acceleration and force transferred into a building superstructure during a seismic event. The sliding system uses a passive method to provide damping and to dissipate some of the kinetic energy to reduce relative displacements. In addition, the SDBIS system includes a self-centering element that will recover the majority of the induced displacement and provide resistance to overturning. In order to apply the SDBIS system to the masonry building industry, dynamic testes were performed to evaluate the structural response of a full-size one-story masonry model that was supported by the SDBIS system. Acceleration time-history results are presented for dynamic tests using the July 21, 1952 Kern County earthquake, Station 1095 Taft Lincoln School record, the May 19, 1940 Imperial Valley earthquake, Station 117 El Centro Array #9 record, the February 9, 1971 San Fernando earthquake, Station 279 Pacoima Dam record, and the January 17, 1994 Northridge earthquake, Station 24436 Tarzana Cedar Hill record ground motions. Test results show the system is effective when used with a masonry structure.     

Design of Architectural Glazing to Resist Earthquakes

A state of the practice summary is presented with regard to the design of architectural glazing to resist earthquakes. Highlighted are new model building code provisions that have been adopted recently by a number of code bodies including the 2003 International Building Code and the 2002 NFPA 5000 Building Code. These model building codes reference seismic design provisions for architectural glazing in ASCE 7-02 and industry-accepted seismic test methods for architectural glazing developed and published by the American Architectural Manufacturers Association. Designers, specifiers, building owners, and glazing system manufacturers need to be aware of these significant changes in the way architectural glazing components are to be designed and tested for earthquake resistance.     

ONDA: Computer Code for Nonlinear Seismic Response Analyses of Soil Deposits

This paper describes a newly developed computer code for performing one-dimensional nonlinear dynamic analysis (ONDA) of soil deposits. The code has been developed by revisiting the 1982 work by Ohsaki with the purpose of simulating the ground response to an earthquake of moderate intensity (i.e., values of peak ground acceleration on stiff soil on the order of 0.15 to 0.25g, which are typical of many sites in Italy). In the Ohsaki model a horizontally stratified soil deposit is idealized as a discrete mechanical system composed of a finite number of lumped masses connected with a series of springs and dashpots. Nonlinearity is modeled by assuming (1) a “backbone” curve that describes the initial monotonic loading of the stress-strain curve, and (2) a “rule” that simulates the unloading-reloading paths and stiffness degradation undergone by soil as seismic excitation progresses. Typically, the backbone curve is obtained from conventional cyclic undrained loading laboratory tests. The rule generally used is the so-called Masing criterion, which assumes that the unload-reload branches of the stress-strain curve have the same shape as the initial loading curve but are affected by a scale factor (n) equal to 2. In this work, the Masing criterion has been modified by assuming a scale factor (n) not necessarily equal to 2. It turns out that a factor n greater than 2 allows the simulation of cyclic hardening, while cyclic softening can be modeled by assuming decreasing values of n even smaller than 2. Pyke proposed in 1979 to use a scale factor (n) lower than 2 to simulate cyclic degradation. According to Pyke, the n parameter is a function of the mobilization factor. The generalization of the Masing criterion allows ONDA to properly simulate the phenomena of soil hardening and soil degradation, giving it the capability to compute the permanent strains developed during a seismic event. The procedure required to evaluate the model parameters is also described in the paper. Note that the laboratory tests examined gave values of n between 2 and 6 for a strain level not greater than 0.3%. In ONDA the numerical solution of the nonlinear equations of motion is obtained using the unconditionally stable Wilson θ algorithm (with θ ≥ 1.37). The new method has been used to predict the seismic response at two sites in Italy. For these case studies, the maximum input acceleration was not greater than 0.3g and the computed shear strains were less than 0.2%. The ONDA results have been compared with those computed with SHAKE, EERA (equivalent-linear analysis), and NERA (nonlinear analysis).     

Small-Strain Behavior of Granular Soils. II: Seismic Response Analyses and Model Evaluation

Using the recorded response at two vertical array sites, the SimSoil model presented in the companion paper is evaluated. The SimSoil model, which describes the small strain nonlinear behavior of granular materials, is implemented as a material model in AMPLE2000, a nonlinear, one-dimensional site response analysis code. Shear wave velocity profiles and laboratory test data available for both the La Cienega site, which was instrumented over 250?m, and the Lotung site, which was instrumented over 47?m, were used to determine SimSoil model parameters. Predictions from AMPLE2000 are compared with the measured response at several elevations for earthquakes that resulted in both nonlinear and nearly linear soil behavior. Using the available laboratory data and known input motions, the predictions of the response at these sites matched the recorded response well for varied magnitudes of shaking with a single set of parameters for each site.     

A comparison of IBC with 1997 UBC for modal response spectrum analysis in standard-occupancy buildings

This paper presents a comparison of the seismic forces generated from a Modal Response Spectrum Analysis (MRSA) by applying the provisions of two building codes, the 1997 Uniform Building Code (UBC) and the 2000-2009 International Building Code (IBC), to the most common ordinary residential buildings of standard occupancy. Considering IBC as the state of the art benchmark code, the primary concern is the safety of buildings designed using the UBC as compared to those designed using the IBC. A sample of four buildings with different layouts and heights was used for this comparison. Each of these buildings was assumed to be located at four different geographical sample locations arbitrarily selected to represent various earthquake zones on a seismic map of the USA, and was subjected to code-compliant response spectrum analyses for all sample locations and for five different soil types at each location. Response spectrum analysis was performed using the ETABS software package. For all the cases investigated, the UBC was found to be significantly more conservative than the IBC. The UBC design response spectra have higher spectral accelerations, and as a result, the response spectrum analysis provided a much higher base shear and moment in the structural members as compared to the IBC. The conclusion is that ordinary office and residential buildings designed using UBC 1997 are considered to be overdesigned, and therefore they are quite safe even according to the IBC provisions.     

Influence of Input Motion and Site Property Variabilities on Seismic Site Response Analysis

Seismic site response analysis evaluates the influence of local soil conditions on earthquake ground shaking. There are multiple sources of potential uncertainty in this analysis; the most significant pertaining to the specification of the input motions and to the characterization of the soil properties. The influence of the selection of input ground motions on equivalent-linear site response analysis is evaluated through analyses performed with multiple suites of input motions selected to fit the same target acceleration response spectrum. The results indicate that a stable median surface response spectrum (i.e., within ±20% of any other suite) can be obtained with as few as five motions, if the motions fit the input target spectrum well. The stability of the median is improved to ±5 to 10% when 10 or 20 input motions are used. If the standard deviation of the surface response spectra is required, at least 10 motions (and preferably 20) are required to adequately model the standard deviation. The influence of soil characterization uncertainty is assessed through Monte Carlo simulations, where variations in the shear-wave velocity profile and nonlinear soil properties are considered. Modeling shear-wave velocity variability generally reduces the predicted median surface motions and amplification factors, most significantly at periods less than the site period. Modeling the variability in nonlinear properties has a similar, although slightly smaller, effect. Finally, including the variability in soil properties significantly increases the standard deviation of the amplification factors but has a lesser effect on the standard deviation of the surface motions.     

Simulating Seismic Response of Cantilever Retaining Walls

Many failures of retaining walls during earthquakes occurred near waterfront. A reasonably accurate evaluation of earthquake effects under such circumstance requires proven analytical models for dynamic earth pressure, hydrodynamic pressure, and excess pore pressure. However, such analytical procedures, especially for excess pore pressure, are not available and hence comprehensive numerical procedures are needed. This paper presents the results of a finite-element simulation of a flexible, cantilever retaining wall with dry and saturated backfill under earthquake loading, and the results are compared with that of a centrifuge test. It is found that bending moments in the wall increased significantly during earthquakes both when backfill is dry or saturated. After base shaking, the residual moment on the wall was also significantly higher than the moment under static loading. Liquefaction of backfill soil contributed to the failure of the wall, which had large outward movement and uneven subsidence in the backfill. The numerical simulation was able to model quite well the main characteristics of acceleration, bending moment, displacement, and excess pore pressure recorded in the centrifuge test in most cases with the simulation for dry backfill slightly better than that for saturated backfill.     

Liquefaction Resistance of Soils: Summary Report from the 1996 NCEER and 1998 NCEER∕NSF Workshops on Evaluation of Liquefaction Resistance of Soils

Following disastrous earthquakes in Alaska and in Niigata, Japan in 1964, Professors H. B. Seed and I. M. Idriss developed and published a methodology termed the “simplified procedure” for evaluating liquefaction resistance of soils. This procedure has become a standard of practice throughout North America and much of the world. The methodology which is largely empirical, has evolved over years, primarily through summary papers by H. B. Seed and his colleagues. No general review or update of the procedure has occurred, however, since 1985, the time of the last major paper by Professor Seed and a report from a National Research Council workshop on liquefaction of soils. In 1996 a workshop sponsored by the National Center for Earthquake Engineering Research (NCEER) was convened by Professors T. L. Youd and I. M. Idriss with 20 experts to review developments over the previous 10 years. The purpose was to gain consensus on updates and augmentations to the simplified procedure. The following topics were reviewed and recommendations developed: (1) criteria based on standard penetration tests; (2) criteria based on cone penetration tests; (3) criteria based on shear-wave velocity measurements; (4) use of the Becker penetration test for gravelly soil; (4) magnitude scaling factors; (5) correction factors for overburden pressures and sloping ground; and (6) input values for earthquake magnitude and peak acceleration. Probabilistic and seismic energy analyses were reviewed but no recommendations were formulated.     

Liquefaction Resistance of Soils: Summary Report from the 1996 NCEER and 1998 NCEER∕NSF Workshops on Evaluation of Liquefaction Resistance of Soils

Following disastrous earthquakes in Alaska and in Niigata, Japan in 1964, Professors H. B. Seed and I. M. Idriss developed and published a methodology termed the “simplified procedure” for evaluating liquefaction resistance of soils. This procedure has become a standard of practice throughout North America and much of the world. The methodology which is largely empirical, has evolved over years, primarily through summary papers by H. B. Seed and his colleagues. No general review or update of the procedure has occurred, however, since 1985, the time of the last major paper by Professor Seed and a report from a National Research Council workshop on liquefaction of soils. In 1996 a workshop sponsored by the National Center for Earthquake Engineering Research (NCEER) was convened by Professors T. L. Youd and I. M. Idriss with 20 experts to review developments over the previous 10 years. The purpose was to gain consensus on updates and augmentations to the simplified procedure. The following topics were reviewed and recommendations developed: (1) criteria based on standard penetration tests; (2) criteria based on cone penetration tests; (3) criteria based on shear-wave velocity measurements; (4) use of the Becker penetration test for gravelly soil; (4) magnitude scaling factors; (5) correction factors for overburden pressures and sloping ground; and (6) input values for earthquake magnitude and peak acceleration. Probabilistic and seismic energy analyses were reviewed but no recommendations were formulated.     

Effect of Asynchronous Earthquake Motion on Complex Bridges. I: Methodology and Input Motion

Based on observed damage patterns from previous earthquakes and a rich history of analytical studies, asynchronous input motion has been identified as a major source of unfavorable response for long-span structures, such as bridges. This study is aimed at quantifying the effect of geometric incoherence and wave arrival delay on complex straight and curved bridges using state-of-the-art methodologies and tools. Using fully parametrized computer codes combining expert geotechnical and earthquake structural engineering knowledge, suites of asynchronous accelerograms are produced for use in inelastic dynamic analysis of the bridge model. Two multi-degree-of-freedom analytical models are analyzed using 2,000 unique synthetic accelerograms with results showing significant response amplification due to asynchronous input motion, demonstrating the importance of considering asynchronous seismic input in complex, irregular bridge design. The paper, Part 1 of a two-paper investigation, presents the development of the input motion sets and the modeling and analysis approach employed, concluding with sample results. Detailed results and implications on seismic assessment are presented in the companion paper: Effect of Asynchronous Motion on Complex Bridges. Part II: Results and Implications on Assessment.     

In-Plane Seismic Response of URM Walls Upgraded with FRP

Recent earthquakes have shown the vulnerability of unreinforced masonry (URM) buildings, which have led to an increasing demand for techniques to upgrade URM buildings. Fiber reinforced polymer (FRP) can provide an upgrading alternative for URM buildings. This paper presents results of dynamic tests investigating the in-plane behavior of URM walls upgraded with FRP (URM-FRP). These tests represent pioneer work in this area (dynamic and in-plane). Five half-scale walls were built, using half-scale brick clay units, and upgraded on one face only. Two moment/shear ratios (1.4 and 0.7), two mortar types (M2.5 and M9), three composite materials (carbon, aramid, and glass), three fiber structures (plates, loose fabric, and grids), and two upgrading configurations (diagonal “X” and full surface shapes) were investigated. The test specimens were subjected to a series of synthetic earthquake motions with increasing intensities on a uniaxial earthquake simulator. The tests validate the effectiveness of the one side upgrading: the upgrading technique improved the lateral resistance of the URM walls by a factor ranging from 1.3 to 2.9; however, the improvement in the lateral drift was less significant. Moreover, no uneven response was observed during the test due to the single side upgrading. Regarding the upgrading configurations, the bidirectional surface type materials (fabrics and grids) applied on the entire surface of the wall (and correctly anchored) can help postpone the three classic failure modes of masonry walls: rocking (“flexural failure”), step cracking, and sliding (“shear failures”). Additionally, in some situations, they will postpone collapse by “keeping the bricks together” under large seismic deformations. On the other hand, the diagonal “X” shape was less successful and premature failure was developed during the test.     

Hydrodynamic Pressures on Arch Dam during Earthquakes

A complete three-dimensional finite difference scheme has been developed and used to analyze the nonlinear hydrodynamic pressures on arch dam during earthquakes. Both free-surface waves and nonlinear convective acceleration are included in the analysis. Various dam shapes and reservoir banks were studied and the characteristic of nonlinear contribution due to curved geometry and free-surface waves are discussed. Numerical experiments have been made to determine the desirable mesh size arrangements and time increments. The effects of surface wave and convective acceleration on hydrodynamic pressure could augment the hydrodynamic pressure to 10% larger than that of a linear analysis. For an earthquake with large ground displacement, a large rise of the water surface could probably happen and be a severe threat to human safety in recreational areas and even cause overflow to affect the safe operation of the arch dam. An empirical formula is given to predict approximate hydrodynamic force.     

Dynamic Experiments and Analyses of a Pile-Group-Supported Structure

Experimental data on the seismic response of a pile-group-supported structure was obtained through dynamic centrifuge model tests, and then used to evaluate a dynamic beam on a nonlinear Winkler foundation (BNWF) analysis method. The centrifuge tests included a structure supported on a group of nine piles founded in soft clay overlying dense sand. This structure was subjected to nine earthquake events with peak accelerations ranging from 0.02 to 0.7g. The centrifuge tests and dynamic analysis methods are described. Good agreement was obtained between calculated and recorded structural responses, including superstructure acceleration and displacement, pile cap acceleration and displacement, pile bending moment and axial load, and pile cap rotation. Representative examples of recorded and calculated behavior for the structure and soil profile are presented. Sensitivity of the dynamic BNWF analyses to the numerical model parameters and site response calculations are evaluated. These results provide experimental support for the use of dynamic BNWF analysis methods in seismic soil-pile-structure interaction problems involving pile-group systems.     

Use of Exact Solutions of Wave Propagation Problems to Guide Implementation of Nonlinear Seismic Ground Response Analysis Procedures

One-dimensional nonlinear ground response analyses provide a more accurate characterization of the true nonlinear soil behavior than equivalent-linear procedures, but the application of nonlinear codes in practice has been limited, which results in part from poorly documented and unclear parameter selection and code usage protocols. In this article, exact (linear frequency-domain) solutions for body wave propagation through an elastic medium are used to establish guidelines for two issues that have long been a source of confusion for users of nonlinear codes. The first issue concerns the specification of input motion as “outcropping” (i.e., equivalent free-surface motions) versus “within” (i.e., motions occurring at depth within a site profile). When the input motion is recorded at the ground surface (e.g., at a rock site), the full outcropping (rock) motion should be used along with an elastic base having a stiffness appropriate for the underlying rock. The second issue concerns the specification of viscous damping (used in most nonlinear codes) or small-strain hysteretic damping (used by one code considered herein), either of which is needed for a stable solution at small strains. For a viscous damping formulation, critical issues include the target value of the viscous damping ratio and the frequencies for which the viscous damping produced by the model matches the target. For codes that allow the use of “full” Rayleigh damping (which has two target frequencies), the target damping ratio should be the small-strain material damping, and the target frequencies should be established through a process by which linear time domain and frequency domain solutions are matched. As a first approximation, the first-mode site frequency and five times that frequency can be used. For codes with different damping models, alternative recommendations are developed.