Accuracy of three-dimensional seismic ground response analysis in time domain using nonlinear numerical simulations

To provide appropriate uses of nonlinear ground response analysis for engineering practice, a three-dimensional soil column with a distributed mass system and a time domain numerical analysis were implemented on the Open Sees simulation platform. The standard mesh of a three-dimensional soil column was suggested to be satisfied with the specified maximum frequency. The layered soil column was divided into multiple sub-soils with a different viscous damping matrix according to the shear velocities as the soil properties were significantly different. It was necessary to use a combination of other one-dimensional or three-dimensional nonlinear seismic ground analysis programs to confirm the applicability of nonlinear seismic ground motion response analysis procedures in soft soil or for strong earthquakes. The accuracy of the three-dimensional soil column finite element method was verified by dynamic centrifuge model testing under different peak accelerations of the earthquake. As a result, nonlinear seismic ground motion response analysis procedures were improved in this study. The accuracy and efficiency of the three-dimensional seismic ground response analysis can be adapted to the requirements of engineering practice.

Amplification of in-plane seismic ground motion by group cavities in layered half-space (Ⅱ): with saturated poroelastic soil layers

As the continuation study on amplification of in-plane seismic ground motion by underground group cavities in layered half-space, this study extends to the case of poroelastic half-space with dry poroelastic and saturated poroelastic soil layers. The influence of poroelastic layers on the amplification of seismic ground motion is studied both in frequency domain and time domain using indirect boundary element method （IBEM）. It is shown that for the example of a saturated poroelastic site in Tianjin under the excitation of Taft wave and E1 Centro wave, the amplification of seismic ground motion in poroelastic case is slightly smaller than that in the elastic case, and the amplification of PGA （peak ground acceleration） and its PRS （peak response spectrum）.. can be increased up to 38.8% and 64.6%; the predominant period of response spectra in poroelastic case becomes shorter to some extent compared with that in the elastic case. It is suggested that the effect of underground group cavities in poroelastic half-space on design seismic ground motion should be considered.

Orthogonal expansion of ground motion and PDEM-based seismic response analysis of nonlinear structures

This paper introduces an orthogonal expansion method for general stochastic processes. In the method, a normalized orthogonal function of time variable t is first introduced to carry out the decomposition of a stochastic process and then a correlated matrix decomposition technique, which transforms a correlated random vector into a vector of standard uncorrelated random variables, is used to complete a double orthogonal decomposition of the stochastic processes. Considering the relationship between the Hartley transform and Fourier transform of a real-valued function, it is suggested that the first orthogonal expansion in the above process is carried out using the Hartley basis function instead of the trigonometric basis function in practical applications. The seismic ground motion is investigated using the above method. In order to capture the main probabilistic characteristics of the seismic ground motion, it is proposed to directly carry out the orthogonal expansion of the seismic displacements. The case study shows that the proposed method is feasible to represent the seismic ground motion with only a few random variables. In the second part of the paper, the probability density evolution method （PDEM） is employed to study the stochastic response of nonlinear structures subjected to earthquake excitations. In the PDEM, a completely uncoupled one-dimensional partial differential equation, the generalized density evolution equation, plays a central role in governing the stochastic seismic responses of the nonlinear structure. The solution to this equation will yield the instantaneous probability density function of the responses. Computational algorithms to solve the probability density evolution equation are described. An example, which deals with a nonlinear frame structure subjected to stochastic ground motions, is illustrated to validate the above approach.

Amplification of in-plane seismic ground motion by group cavities in layered half-space (Ⅰ)

Amplification of in-plane seismic ground motion by underground group cavities in layered half-space is studied both in frequency domain and time domain by using indirect boundary element method （IBEM）, and the effect of cavity interval and spectrum of incident waves on the amplification are studied by numerical examples. It is shown that there may be large interaction between cavities, and group cavities with certain intervals may have significant amplification to seismic ground motion. The amplification of PGA （peak ground acceleration） and its PRS （peak response spectrum） can be increased up to 45.2% and 84.4%, for an example site in Tianjin, under the excitation of Taft wave and E1 Centro wave; and group cavities may also affect the spectra of the seismic ground motion. It is suggested that the effect of underground group cavities on design seismic ground motion should be considered.

Method of forecasting seismic energy induced by longwall exploitation based on changes in ground subsidence

A method of forecasting total seismic energy induced by longwall exploitation, based on changes in ground subsidence, is presented in the form of a linear regression model with one with one independent variable. In the method, ground subsidence is described with a cross-section area of a subsidence trough Pw along a line of observations in the direction of an advancing longwall front, approximately along the axis of the longwall area. Total seismic energy is determined on the basis of seismic energy data of tremors induced by exploitation. The presentation consists of a detailed method and evaluation of its predictive ability for the area of longwall exploitation within the region of one of the coal mines in the Upper Silesian Coal Basin. This method can be used for forecasting the total seismic energy released by tremors within the area directly connected with the exploitation, in which the seismic activity induced by this exploitation occurs. The estimation of the parameters of the determined model should each time be carried out with investigations of the correctness of the model. The method cannot be applied when the number of recorded phenomena is small and when there is insufficient data to make it possible to calculate the index Pw.

Probabilistic seismic hazard assessment of Kazakhstan and Almaty city in peak ground accelerations

As for many post-soviet countries, Kazakhstan’s building code for seismic design was based on a deterministic approach. Recently, Kazakhstan seismologists are engaged to adapt the PSHA(probabilistic hazard assessment) procedure to the large amount of available geological, geophysical and tectonic Kazakh data and to meet standard requirements for the Eurocode 8. The new procedure has been used within National projects to develop the Probabilistic GSZ(General Seismic Zoning) maps of the Kazakhstan territory and the SMZ(Probabilistic Seismic Microzoning) maps of Almaty city. They agree with the seismic design principles of Eurocode 8 and are expressed in terms of not only seismic intensity,but also engineering parameters(peak ground acceleration PGA). The whole packet of maps has been developed by the Institute of Seismology, together with other Kazakhstan Institutions. Our group was responsible for making analysis in PGA. The GSZ maps and hazard assessment maps for SMZ in terms of PGA for return periods 475 and 2475 years are considered in the article.

The effect of local irregular topography on seismic ground motion

Effects of irregUlar topography on ground motion for incident P, SV and the propagation of Rayleigh waves are studied by combining finite element method with modified transmitting boundary. TheoretiCal models include isolated protrUding topography and similar adjacent Protruding topography. The concluaion drawn from thisstudy is that the effects Of isolated protruding topography are remarkably larger for Rayleigh wave propagation than for P and SV they waves; Considering adjacent irregUlar toography ground motion is amplified, the duration of ground motion becomes longer and the speCtral ratios exhibit narrowband peaks Considering adjacent irregular topography and Rayleigh wave Propagation, the theoretical results wb more approach the results obtained in practice.

On the realization of seismic microzonation of Almaty (Kazakhstan) in ground accelerations based on the “continual” approach

Seismic microzonation for Almaty city for the first time use probabilistic approach and hazard is expressed in terms of not only macroseismic intensity,but also Peak Ground Acceleration(PGA).To account for the effects of local soil conditions,the continual approach proposed by A.S.Aleshin[1,2]was used,in which soil coefficients are a function of the continuously changing seismic rigidity.Soil coefficients were calculated using the new data of geological and geophysical surveys and findings of previous geotechnical studies.The used approach made it possible to avoid using soil categories and a jump change in characteristics of soil conditions and seismic impact.The developed seismic microzonation maps are prepared for further introduction into the normative documents of the Republic of Kazakhstan.

Analytical investigations of seismic responses for reinforced concrete bridge columns subjected to strong near-fault ground motion

Strong near-fault ground motion, usually caused by the fault-rupture and characterized by a pulse-like velocity- wave form, often causes dramatic instantaneous seismic energy （Jadhav and Jangid 2006）. Some reinforced concrete （RC） bridge columns, even those built according to ductile design principles, were damaged in the 1999 Chi-Chi earthquake. Thus, it is very important to evaluate the seismic response of a RC bridge column to improve its seismic design and prevent future damage. Nonlinear time history analysis using step-by-step integration is capable of tracing the dynamic response of a structure during the entire vibration period and is able to accommodate the pulsing wave form. However, the accuracy of the numerical results is very sensitive to the modeling of the nonlinear load-deformation relationship of the structural member. FEMA 273 and ATC-40 provide the modeling parameters for structural nonlinear analyses of RC beams and RC columns. They use three parameters to define the plastic rotation angles and a residual strength ratio to describe the nonlinear load- deformation relationship of an RC member. Structural nonlinear analyses are performed based on these parameters. This method provides a convenient way to obtain the nonlinear seismic responses of RC structures. However, the accuracy of the numerical solutions might be further improved. For this purpose, results from a previous study on modeling of the static pushover analyses for RC bridge columns （Sung et al. 2005） is adopted for the nonlinear time history analysis presented herein to evaluate the structural responses excited by a near-fault ground motion. To ensure the reliability of this approach, the numerical results were compared to experimental results. The results confirm that the proposed approach is valid.

Assessment of seismic ground motion amplification and liquefaction at a volcanic area characterized by residual soils

A seismic hazard was assessed related to site effects at Abbadia San Salvatore, central Italy, on the Mt. Amiata slopes, an ancient volcanic area characterized by residual soils(thick layers of loose to dense sands originated from weathering of the trachydacitic lavas). The seismic ground amplification and soil liquefaction related to these layers were recognized as the major seismic hazards for the area.Geological, geophysical, and geotechnical surveys were carried out on the volcanic rocks. The Horizontal-to-Vertical Spectral Ratio(HVSR) analysis of 252 noise measurements and 29 shear-wave velocity models of the subsoil allowed a seismic microzonation of the studied area, distinguished by thick weathered volcanic sands and shear-wave impedance contrast with respect to the seismic bedrock(volcanic bedrock). The differentiation of classified zones allowed recognition of areas characterized by residual(almost undisturbed) soils from those with soils probably affected by flowing water. The analysis of hazards revealed that peak acceleration by seismic amplification of ground motion exceeded the value set by the national rules(0.175 g) in a restricted area of the zone characterized by the most perturbed soils(Zone D);the potential occurrence of soil liquefaction was also greater in this zone. Finally, the study showed potential high hazards due to site effects of the volcanic mountainous area characterized by residual soils as opposed to an alluvial plain formed by volcanic debris where these effects have generally been more recognized.

Simplified method for simulation of ergodic spatially correlated seismic ground motion

A simplified method for the simulation of the ergodic spatially correlated seismic ground motion is proposed based on the commonly used original spectral representation method. To represent the correlation in the ground motion, the phase angles are given by explicit terms with a clear physical meaning. By these explicit terms, the computational efficiency can be improved by converting the decomposition of the complex cross-spectral matrix into the decomposition of the real incoherence coefficient matrix. Double-indexing frequencies are introduced to simulate the ergodic seismic ground motion, and the ergodic feature of the improved method is demonstrated theoretically. Subsequently, an explicit solution of the elements of the lower triangular matrix under the Cholesky decomposition is given. With this explicit solution, the improved method is simplified, and the computational efficiency can be improved greatly by avoiding the repetitive Cholesky decomposition of the cross-spectral matrix in each frequency step. Finally, a numerical example shows the good characteristic of the improved method.

Seismic Responses of Asymmetric Base-Isolated Structures under Near-Fault Ground Motion

An inter-story shear model of asymmetric base-isolated structures incorporating deformation of each isolation bearing was built, and a method to simultaneously simulate bi-directional near-fault and far-field ground motions was proposed. A comparative study on the dynamic responses of asymmetric base-isolated structures under near-fault and far-field ground motions were conducted to investigate the effects of eccentricity in the isolation system and in the superstructures, the ratio of the uncoupled torsional to lateral frequency of the superstructure and the pulse period of near-fault ground motions on the nonlinear seismic response of asymmetric base-isolated structures. Numerical results show that eccentricity in the isolation system makes asymmetric base-isolated structure more sensitive to near-fault ground motions, and the pulse period of near-fault ground motions plays an import role in governing the seismic responses of asymmetric base-isolated structures.

Comparison of ground truth location of earthquake from InSAR and from ambient seismic noise: A case study of the 1998 Zhangbei earthquake

Because ambient seismic noise provides estimated Green’s function （EGF） between two sites with high accuracy, Rayleigh wave propagation along the path connecting the two sites is well resolved. Therefore, earthquakes which are close to one seismic station can be well located with calibration extracting from EGF. We test two algorithms in locating the 1998 Zhangbei earthquake, one algorithm is waveform-based, and the other is traveltime-based. We first compute EGF between station ZHB （a station about 40 km away from the epicenter） and five IC/IRIS stations. With the waveform-based approach, we calculate 1D synthetic single-force Green’s functions between ZHB and other four stations, and obtain traveltime corrections by correlating synthetic Green’s functions with EGFs in period band of 10–30 s. Then we locate the earthquake by minimizing the differential travel times between observed earthquake waveform and the 1D synthetic earthquake waveforms computed with focal mechanism provided by Global CMT after traveltime correction from EGFs. This waveform-based approach yields a location which error is about 13 km away from the location observed with InSAR. With the traveltime-based approach, we begin with measuring group velocity from EGFs as well as group arrival time on observed earthquake waveforms, and then locate the earthquake by minimizing the difference between observed group arrival time and arrival time measured on EGFs. This traveltime-based approach yields accuracy of 3 km, Therefore it is feasible to achieve GT5 （ground truth location with accuracy 5 km） with ambient seismic noises. The less accuracy of the waveform-based approach was mainly caused by uncertainty of focal mechanism.

Seismic Performance Assessment of Reinforced Concrete Box Culverts under Near and Far Fault Seismic Ground Motion Records

Studying the critical response characteristics of reinforced concrete box culverts with diverse geometrical configurations under seismic excitations is a necessary step to develop a reasonable design method. In this work, numerical analysis and assessment of reinforced concrete box culverts for seismic loading in addition to standard static loading from dead and live loads is conducted, aiming to highlight the critical difference in the seismic performances between two and three cell box culverts under near and far-fault ground motion. The results show how and where the seismic loading alters the responses of seismic loading of the models including the effect on safety and failure. The geometrical configurations of the culvert combined with the loading scenarios also significantly influence the magnitude and distribution of the seismic responses. The findings of this work shed light on the critical role of the geometrical configurations and shaking event in the seismic responses of reinforced concrete box culverts and this procedure can be applied as seismic assessment method to any culvert shape, size, and material.

Condition Assessment of August A. Busch Bridge Deck Using Portable Seismic Property Analyzer and Ground Penetrating Radar

Ground penetrating radar (GPR) and the portable seismic property analyzer (PSPA) have been extensively used in the past two decades for monitoring, quantifying, and mapping the deterioration of bridge decks. Using PSPA and GPR ensures regular monitoring of bridge conditions, leads to the early detection of deterioration. This research is to address the condition of August A. Busch bridge deck owned by the Missouri Department of Conservation. Visual inspection, GPR, and PSPA data were acquired on the bridge deck. Over 90% of the bridge deck was in fair to good condition with an average compressive strength of over 2500 psi. GPR data showed no indication of significant deterioration. The overall bridge deck was determined to be in fair to good condition.

Study on the Relationship of Seismic Design Ground Motion at Different Fortification Levels in the Xinjiang Region

The seismic risk analysis results of 79 cities in Xinjiang are presented, and the bedrock peak ground accelerations under three seismic levels and their ratios are discussed. Then, the relationship between earthquake environments and the seismic risk analysis results of different exceeding probabilities are researched. The results show that minor and major earthquake motion parameters calculated from moderate earthquakes do not have a consistent probability and the ratio of bedrock peak accelerations under different exceedance probabilities are closely correlated with earthquake environments.

The Relationship Among Bedrock Seismic Ground Motion Parameters with Different Exceedance Probabilities in the Panxi Area

Based on the calculation of the bedrock effective peak acceleration (EPA) zoning map in the Panxi area, the ratios of EPA with exceedance probabilities of 63%, 5%, 3%, 2% and 1% over 50 years to that of 10% in 50 years are 0.302, 1.30, 1.55, 1.76 and 2.14, respectively. The seismic effect will be conservative and safe if taking this zoning map as the earthquake resistant fortification level and following the relevant rules of the Code for Seismic Design of Buildings (GBJ11 89) to calculate the seismic effect. Furthermore, the main factors that influence the A10/A63 ratios have been found to be the attenuation relationship of seismic ground motion, the division of seismic potential source regions and the seismicity parameters. These achievements are helpful to the spreading and applying of the zoning map.

Engineering characteristics of near-fault vertical ground motions and their effect on the seismic response of bridges

A wide variety of near-fault strong ground motion records were collected from various tectonic environments worldwide and were used to study the peak value ratio and response spectrum ratio of the vertical to horizontal component of ground motion, focusing on the effect of earthquake magnitude, site conditions, pulse duration, and statistical component. The results show that both the peak value ratio and response spectrum ratio are larger than the 2/3 value prescribed in existing seismic codes, and the relationship between the vertical and horizontal ground motions is comparatively intricate. In addition, the effect of the near-fault ground motions on bridge performance is analyzed, considering both the material nonlinear characteristics and the P-A effect.

Seismic Ground Motion Zoning Maps of the Pangxi Region

The seismotectonic environment and seismic activity in Southwest China region were studied based on new data and new results obtained during the Eighth and Ninth Five-Year Plans, the seismic areas and zones and potential seismic source zones were determined, and the relation between seismic activity parameters and ground motion attenuation was determined. Finally the seismic ground motion zoning maps of the Pangxi region was compiled by using the multi-parameter and multi-scheme method.