Appropriate viscous damping for nonlinear time‐history analysis of base‐isolated reinforced concrete buildings

There is no consensus at the present time regarding an appropriate approach to model viscous damping in nonlinear time‐history analysis of base‐isolated buildings because of uncertainties associated with quantification of energy dissipation. Therefore, in this study, the effects of modeling viscous damping on the response of base‐isolated reinforced concrete buildings subjected to earthquake ground motions are investigated. The test results of a reduced‐scale three‐story building previously tested on a shaking table are compared with three‐dimensional finite element simulation results. The study is primarily focused on nonlinear direct‐integration time‐history analysis, where many different approaches of modeling viscous damping, developed within the framework of Rayleigh damping are considered. Nonlinear direct‐integration time‐history analysis results reveal that the damping ratio as well as the approach used to model damping has significant effects on the response, and quite importantly, a damping ratio of 1% is more appropriate in simulating the response than a damping ratio of 5%. It is shown that stiffness‐proportional damping, where the coefficient multiplying the stiffness matrix is calculated from the frequency of the base‐isolated building with the post‐elastic stiffness of the isolation system, provides reasonable estimates of the peak response indicators, in addition to being able to capture the frequency content of the response very well. Furthermore, nonlinear modal time‐history analyses using constant as well as frequency‐dependent modal damping are also performed for comparison purposes. It was found that for nonlinear modal time‐history analysis, frequency‐dependent damping, where zero damping is assigned to the frequencies below the fundamental frequency of the superstructure for a fixed‐base condition and 5% damping is assigned to all other frequencies, is more appropriate, than 5% constant damping. Copyright © 2013 John Wiley & Sons, Ltd.     

Constitutive model for soil amplification of ground shaking: Parameter calibration,comparisons, validation

A one-dimensional constitutive model, developed for the nonlinear ground response analysis of layered soil deposits, is calibrated and validated experimentally in this paper. The small number of parameters renders the model easily implementable, yet quite flexible in effectively reproducing almost any type of experimentally observed hysteretic soil behavior. In particular, the model generates realistic shear modulus and damping curves as functions of shear strain, as well as stress–strain hysteresis loops. The model is calibrated against three sets of widely-used published shear modulus and damping (G : γ and ξ : γ) curves and a library of parameter values is assembled to facilitate its use. The model, along with a developed explicit finite-difference code, NL-DYAS, for analyzing the wave propagation in layered hysteretic soil deposits, is tested against established constitutive models and numerical tools such as Cyclic1D [12] and SHAKE [42], and validated against experimental data from two centrifuge tests. Emphasis is given on the proper assessment of the V_s profile in the centrifuge tests, on the role of soil nonlinearity, and on comparisons of two inelastic codes (NL-DYAS and Cyclic1D) with equivalent linear (SHAKE) analysis.     

Simulating the Taipei basin response by numerical modeling of wave propagation

Taipei, the capital of Taiwan, suffered from destructive earthquakes four times during the 20th century (M _L = 7.3 on April 15, 1909; M _L = 6.8 on November 15, 1986; the Chi–Chi M _L = 7.3 earthquake on September 21, 1999; and M _L = 6.8 on March 31, 2002). Analysis of recorded data shows a strong dependence of spectral amplification in the Taipei Basin on earthquake depth and azimuth. At low frequencies (f < 3 Hz) significant larger amplifications are observed for shallow earthquakes as compared to intermediate depth events. The former ones also display strong azimuthal dependence. As structures with large response periods such as bridges and tall buildings are sensitive to these low frequencies the understanding of the associated wave effects within the basin and their role for site effect amplification is critical. The tool we employ is 3D finite-difference modeling of wave propagation of incident wave fronts. The available detailed model of the basin allows studying the wave effects. Modeling clearly reveals that basin edge effects as observed in data are related to surface wave generation at the basin edges with a high degree of azimuthal dependency. The reproduced site amplification effects are in qualitative agreement with the observations from strong motion data.     

瑞利阻尼矩阵对深厚场地地震反应的影响分析

介绍场地时域分析基本理论和几种常见的阻尼矩阵确定方法,并基于有限差分方法,采用实际地震记录,考虑输入地震动强度和频谱特性,分别选择5种不同的瑞利阻尼矩阵形式,比较分析其对深厚场地时域动力反应的影响。结果表明:不同阻尼矩阵对加速度和反应谱均有影响,影响程度和输入地震动的频谱特性有关;当土层基频接近输入地震动卓越频率时,只根据土层基频确定阻尼矩阵是可行的,但当输入地震动卓越频率高于或远高于土层基频时,该方法会严重低估地表加速度峰值和反应谱谱值,时域分析中阻尼的确定应综合考虑土层特性和地震动频谱特性。     

Earthquake ground motion modeling of induced seismicity in the Groningen gas field

A physics‐based numerical approach is used to characterize earthquake ground motion due to induced seismicity in the Groningen gas field and to improve empirical ground motion models for seismic hazard and risk assessment. To this end, a large‐scale (20 km × 20 km) heterogeneous 3D seismic wave propagation model for the Groningen area is constructed, based on the significant bulk of available geological, geophysical, geotechnical, and seismological data. Results of physics‐based numerical simulations are validated against the ground motion recordings of the January 8, 2018, M_L 3.4 Zeerijp earthquake. Taking advantage of suitable models of slip time functions at the seismic source and of the detailed geophysical model, the numerical simulations are found to reproduce accurately the observed features of ground motions at epicentral distances less than 10 km, in a broad frequency range, up to about 8 Hz. A sensitivity analysis is also addressed to discuss the impact of 3D underground geological features, the stochastic variability of seismic velocities and the frequency dependence of the quality factor. Amongst others, results point out some key features related to 3D seismic wave propagation, such as the magnitude and distance dependence of site amplification functions, that may be relevant to the improvement of the empirical models for earthquake ground motion prediction.     

Centrifuge modeling of seismic response of layered soft clay

Centrifuge modeling is a valuable tool used to study the response of geotechnical structures to infrequent or extreme events such as earthquakes. A series of centrifuge model tests was conducted at 80g using an electro-hydraulic earthquake simulator mounted on the C-CORE geotechnical centrifuge to study the dynamic response of soft soils and seismic soil–structure interaction (SSI). The acceleration records at different locations within the soil bed and at its surface along with the settlement records at the surface were used to analyze the soft soil seismic response. In addition, the records of acceleration at the surface of a foundation model partially embedded in the soil were used to investigate the seismic SSI. Centrifuge data was used to evaluate the variation of shear modulus and damping ratio with shear strain amplitude and confining pressure, and to assess their effects on site response. Site response analysis using the measured shear wave velocity, estimated modulus reduction and damping ratio as input parameters produced good agreement with the measured site response. A spectral analysis of the results showed that the stiffness of the soil deposits had a significant effect on the characteristics of the input motions and the overall behavior of the structure. The peak surface acceleration measured in the centrifuge was significantly amplified, especially for low amplitude base acceleration. The amplification of the earthquake shaking as well as the frequency of the response spectra decreased with increasing earthquake intensity. The results clearly demonstrate that the layering system has to be considered, and not just the average shear wave velocity, when evaluating the local site effects.     

Propagation ofSH-waves in a heterogeneous visco-elastic layer overlying a homogeneous voigt medium

Summary The propagation ofSH-waves in a vosco-elastic layer having variable material properties has been studied to explain the damping phenomenon of the wave motion. The damping coefficient may be calculated and the effect of heterogeneity on such type of wave motion may be observed from the obtained frequency equation.     

Evaluation of site effects on ground motions based on equivalent linear site response analysis and liquefaction potential in Chennai,south India

We study local site effects with detailed geotechnical and geophysical site characterization to evaluate the site-specific seismic hazard for the seismic microzonation of the Chennai city in South India. A Maximum Credible Earthquake (MCE) of magnitude 6.0 is considered based on the available seismotectonic and geological information of the study area. We synthesized strong ground motion records for this target event using stochastic finite-fault technique, based on a dynamic corner frequency approach, at different sites in the city, with the model parameters for the source, site, and path (attenuation) most appropriately selected for this region. We tested the influence of several model parameters on the characteristics of ground motion through simulations and found that stress drop largely influences both the amplitude and frequency of ground motion. To minimize its influence, we estimated stress drop after finite bandwidth correction, as expected from an M6 earthquake in Indian peninsula shield for accurately predicting the level of ground motion. Estimates of shear wave velocity averaged over the top 30 m of soil (V_S30) are obtained from multichannel analysis of surface wave (MASW) at 210 sites at depths of 30 to 60 m below the ground surface. Using these V_S30 values, along with the available geotechnical information and synthetic ground motion database obtained, equivalent linear one-dimensional site response analysis that approximates the nonlinear soil behavior within the linear analysis framework was performed using the computer program SHAKE2000. Fundamental natural frequency, Peak Ground Acceleration (PGA) at surface and rock levels, response spectrum at surface level for different damping coefficients, and amplification factors are presented at different sites of the city. Liquefaction study was done based on the V_S30 and PGA values obtained. The major findings suggest show that the northeast part of the city is characterized by (i) low V_S30 values (<?200 m/s) associated with alluvial deposits, (ii) relatively high PGA value, at the surface, of about 0.24 g, and (iii) factor of safety and liquefaction below unity at three sites (no. 12, no. 37, and no. 70). Thus, this part of the city is expected to experience damage for the expected M6 target event.     

Viscous damping formulation and high frequency motion propagation in non-linear site response analysis

Non-linear time domain site response analysis is widely used in evaluating local soil effects on propagated ground motion. This approach has generally provided good estimates of field behavior at longer periods but has shortcomings at relatively shorter periods. Viscous damping is commonly employed in the equation of motion to capture damping at very small strains and employs an approximation of Rayleigh damping using the first natural mode only. This paper introduces a new formulation for the viscous damping using the full Rayleigh damping. The new formulation represents more accurately wave propagation for soil columns greater than 50 m thick and improves non-linear site response analysis at shorter periods. The proposed formulation allows the use of frequency dependent viscous damping. Several examples, including a field case history at Treasure Island, California, demonstrate the significant improvement in computed surface response using the new formulation.     

Time-domain forward and inverse modeling of lossy soils with frequency-independent Q for near-surface applications

In this paper, we are concerned with a full-waveform-based methodology that allows the simultaneous imaging of the soil's stiffness and attenuating properties, using solely the soil's surficial response to probing waves.To date, field observations of small-strain wave attenuation in geomaterials at moderate spatial scales suggest that a commonly used metric of intrinsic and apparent attenuation, the seismic quality factor Q, is frequency-independent for a wide part of the frequency spectrum, including the frequency range of interest to seismic applications. We discuss first the forward simulation of waves in near-surface soil deposits directly in the time-domain using simplified models that adequately approximate nearly frequency-independent Q. To this end, we first review various attenuation models that aim at reproducing the frequency-independent Q behavior, and conclude, supported by site analyses, that, even though a generalized Maxwell body with eight Maxwell elements in parallel (GMB8) provides the best fit to frequency-independent Q, we favor a version of it with fewer parameters (GMB2), in order to reduce modeling complexity, while still retaining good agreement with the GMB8 model.We report on forward site analyses that lend credence to the choice of the GMB2 simplified model. We, then, use the GMB2 constitutive relation in the context of full-waveform inversion, and report on numerical experiments that lead to the imaging of the soil's properties in heterogeneous semi-infinite domains.     

Numerical analysis of seismic wave amplification in Nice (France) and comparisons with experiments

The analysis of site effects is very important since the amplification of seismic motion in some specific areas can be very strong. In this paper, the site considered is located in the centre of Nice on the French Riviera. Site effects are investigated considering a numerical approach (Boundary Element Method) and are compared with experimental results. The experimental results are obtained thanks to real earthquakes (weak motion) and microtremor measurements. The investigation of seismic site effects through numerical approaches is interesting because it shows the dependency of the amplification level on such parameters as wave velocity in surface soil layers, velocity contrast with deep layers, seismic wave type, incidence and damping.In this specific area of Nice, experimental measurements obtained for weak motion lead to strong site effects. A one-dimensional (1D) analytical analysis of amplification does not give a satisfactory estimation of the maximum reached levels. A boundary element model is then proposed considering different wave types (SH, P, SV) as the seismic loading. The alluvial basin is successively assumed as an isotropic linear elastic medium and an isotropic linear viscoelastic solid with Zener type behaviour (standard solid). The influence of frequency and incidence is analysed. The thickness of the surface layer, its mechanical properties, its general shape as well as the seismic wave type involved have a great influence on the maximum amplification and the frequency for which it occurs. For real earthquakes, the numerical results are in very good agreement with experimental measurements for each motion component. The boundary element method leads to amplification values very close to the actual ones and much larger than those obtained in the 1D case. Two-dimensional basin effects are then very strong and are well reproduced numerically.     

Implementation of configuration dependent stiffness proportional damping for the dynamics of rigid multi-block systems

The distinct element method (DEM) has been used successfully for the dynamic analysis of rigid block systems. One of many difficulties associated with DEM is modeling of damping. In this paper, new procedures are proposed for the damping modeling and its numerical implementation in distinct element analysis of rigid multi-block systems. The stiffness proportional damping is constructed for the prescribed damping ratio, based on the non-zero fundamental frequency effective during the time interval while the boundary conditions remain essentially constant. At this time interval, the fundamental frequency can be estimated without complete eigenvalue analysis. The damping coefficients will vary while the damping ratio remains the same throughout the entire analysis. A new numerical procedure is developed to prevent unnecessary energy loss that can occur during the separation phases. These procedures were implemented in the development of the distinct element method for the dynamic analyses of piled multi-block systems. The analysis results for the single-block and two-block systems were in a good agreement with the analytic predictions. Applications to the seismic analyses of piled fourblock systems revealed that the new procedures can make a significant difference and may lead to much-improved results.     

有限元算法在声波方程数值模拟中的频散分析

针对有限元算法在地震波数值模拟中的数值频散问题,利用集中质量矩阵双线性插值有限元算法,推导了二维声波方程的频散函数.在此基础上采用定量分析方法,对比分析了网格纵横长度比变化时的入射方向、空间采样间隔、地震波频率以及地层速度对数值频散的影响.数值算例和模型正演结果表明:当采用集中质量矩阵双线性插值有限元算法时,为了有效地压制数值频散,在所使用震源子波的峰值频率对应的波长内,采样点数目应不少于20个;减小网格长度的纵横比可以有效地抑制入射角(波传播方向与z轴的夹角)较小的地震波的数值频散;地震波频率越高,传播速度越慢,频散越严重,尤其是当相速度与其所对应的频率比值小于2倍空间采样间隔时,不仅会出现严重的数值频散,还会出现假频现象.      

The effects of torsion and motion coupling in site response estimation

Soil amplification characteristics are investigated using data from the Chibaken‐Toho‐Oki earthquake and its aftershocks recorded at Chiba dense array in Japan. The frequency‐dependent amplification function of soil is calculated using uphole‐to‐downhole spectral ratio analysis, considering the horizontal components of shear wave. The identified spectral ratios consistently demonstrate the splitting of peaks in their resonance frequencies and low amplification values in comparison with a 1D model. The torsional behaviour and horizontal ground motion coupling are clarified as the reasons for these phenomena at the site. To prove the hypothesis, the torsional motion is directly evaluated using the data of the horizontal dense array in different depths at the site. The comparison between Fourier spectra of torsional motion and identified transfer functions reveals the peaks at the same frequencies. The wave equation including torsion and horizontal motion coupling is introduced and solved for the layered media by applying wave propagation theory. Using the developed model, the effects of torsional motion with horizontal motion coupling on soil transfer function are numerically examined. Splitting and low amplification at resonance frequencies are confirmed by the results of numerical analysis. Furthermore, the ground motion in two horizontal directions at the site is simulated using site geotechnical specification and optimizing the model parameters. The simulated and recorded motions demonstrate good agreement that is used to validate the hypothesis. In addition, the spectral density of torsional ground motions are compared with the calculated one and found to be well predicted by the model. Finally, the results are used to explain the overestimation of damping in back‐calculation of dynamic soil properties using vertical array data in small strain level. Copyright © 2003 John Wiley & Sons, Ltd.     

Caustics in a periodically layered transversely isotropic medium with vertical symmetry axis

Wave propagation in a finely layered medium is a very important topic in seismic modelling and inversion. Here we analyse non‐vertical wave propagation in a periodically layered transversely isotropic (VTI) medium and show that the evanescent (attenuation) zones in the frequency‐horizontal slowness domain result in caustics in the group velocity domain. These caustics, which may appear for both the quasi‐compressional (qP) and quasi‐shear (qSV) wave surfaces are frequency dependent but display weak dependence at low frequencies. The caustics computed for a specific frequency differ from those observed at the low‐ and high‐frequency limits. We illustrate these caustics with a few numerical examples and snapshots computed for both qP‐ and qSV‐wave types.     

Finite element simulations of wave propagation in soils using a Viscoelastic model

Wave propagation in soil is dependent on both the stiffness and the material damping of the soil. While some researchers have performed finite element modelling of resonant column tests and wave propagation in soil, most do not describe the methodology in detail and there is little or no verification of the correctness of the model. Viscoelastic model has been used to model wave propagation in soil. However, the determination of the parameters in the viscoelastic model is complicated and the parameters may not be related to the soil properties. This paper presents a simplified viscoelastic model with soil parameters obtainable from advanced geotechnical testing to simulate wave propagation in soil medium taking into account of material damping. The viscoelastic material model was first calibrated by replicating torsional, longitudinal and flexural modes resonant column tests. The relationships between the parameters of the simplified viscoelastic model and their corresponding stiffness and damping properties were investigated. An equation was proposed to correlate the decay constant used in the simplified viscoelastic model and the material damping ratio obtained through the application of the logarithmic decrement method on the modelled resonant column test results. The simplified viscoelastic model was then evaluated by modelling wave propagation in a semi-infinite medium. Results indicated that the viscoelastic model with parameters as proposed in this paper is able to model wave propagation in soils.     

Dynamic stiffness of foundation embedded in layered halfspace based on wave propagation in cones

An Erratum has been published for this article in Earthquake Engineering & Structural Dynamics 33(6) 2004, 793. The dynamic stiffness of a foundation embedded in a multiple‐layered halfspace is calculated postulating one‐dimensional wave propagation in cone segments. In this strength‐of‐materials approach the sectional property of the cone segment increases in the direction of wave propagation. Reflections and refractions with waves propagating in corresponding cone segments occur at layer interfaces. Compared to rigorous procedures the novel method based on cone segments is easy to apply, provides conceptual clarity and physical insight in the wave propagation mechanisms. This method postulating one‐dimensional wave propagation in cone segments with reflections and refractions at layer interfaces is evaluated, calculating the dynamic stiffness of a foundation embedded in a multiple‐layered halfspace. For sites resting on a flexible halfspace and fixed at the base, engineering accuracy (deviation of ±20%) is achieved for all degrees of freedom with a vast parameter variation. The behaviour below the cut‐off frequency in an undamped site fixed at its base is also reliably predicted. The accuracy is, in general, better than for the method based on cone frustums, which can lead to negative damping. Copyright © 2003 John Wiley & Sons, Ltd.     

Numerical estimation of site effects in the city of Arica due to natural soil variability using the spectral elements method

In recent years, because of a seismic gap of 136 years in the north of Chile, several studies have been carried out with the purpose of characterizing the dynamic properties of the soils of the most populated cities in this area. The purpose of this work is to estimate, by numerical modeling, the seismic amplification of a densely populated area of the city of Arica. Spatial distribution of main soil dynamic properties has been obtained from a detailed geophysical survey, including surface-wave based methods and gravimetry. To estimate the site effects, we solve the wave propagation equation in full heterogeneous media by the spectral element method. This method allows including the topography, the irregular contact between soils and bedrock and heterogeneities of main materials’ properties along the computational domain. The inelastic behavior of the soil has been considered, using stiffness degradation and damping curves. The computations were done using the high performance open-source numerical code SPEED. Results from the distribution of the peak ground acceleration are compared against standard 1D horizontally layered modeling using equivalent linear and nonlinear viscoelastic approximations of nonlinear soil behavior. Linear viscoelastic results indicate that the north area of the city of Arica is not especially prone to tridimensional site effects induced by soil’s natural variability and the apparent bedrock geometry. Nevertheless, non-conservative differences of about 20% between 1D v/s 3D analyses were found under inelastic material behavior. In general, the two 3D approaches followed to include soils’ nonlinearity provide similar site amplification estimates.     

基于等效线性化方法的一维土层地震反应通用计算程序对比

基于等效线性化的一维土层地震反应计算是目前国内外普遍采用的方法,国外的SHAKE91、DEEPSOIL和我国的LSSRLI-1即是根据这一方法编制的通用计算程序。本文采用这3个程序进行了不同地震波、不同输入地震动幅值下不同场地类型的土层地震反应计算,并对三者的结果进行了全面的比较分析。结果表明:①SHAKE91和DEEPSOIL程序的计算结果完全相同;②当土层最大剪应变均采用时域计算时,LSSRLI-1程序的计算结果与SHAKE91和DEEPSOIL程序基本相同,但有微小差别,其原因是:在基于等效剪应变通过离散形式的剪切模量和阻尼比随等效剪应变变化的关系曲线确定等效剪切模量和阻尼比时,DEEPSOIL和SHAKE91采用的插值方法与LSSRLI-1不同;③当LSSRLI-1程序采用频域经验关系计算土层最大剪应变时,特别是在强地震动输入下得到的土层地表加速度峰值和加速度反应谱与另外两个程序的计算结果有差别,且土层最大剪应变随着输入加速度的增大出现较大的差别。因此,本文建议:当采用LSSRLI-1程序计算土层地震响应时,应使用程序中的时域解方法代替以往默认的频域经验关系方法。     

钻井中首波的共振与衰减

本文从P波首波的严格割线积分出发,首次给出了揭示其共振与衰减特性的渐近展开式,通过数值计算,对这些公式进行了检验,证明它们与严格解符合的相当好.解析的和数值结果一致表明:首波具有共振特性,频谱是连续的;振幅随源距z按1/z规律衰减;各共振峰幅度随岩石泊松比v的增大按[v/（1-v）]²规律增大.时域计算表明,所谓“滤波效应”正是首波共振激发特性的反映,传播过程中不频散.