裸眼井中弹性波传播的非对称模式的数值研究

本文从数值上研究了裸眼井中弹性波传播的非对称模式,给出了合成微地震图和导波（弯曲波）的频散曲线,发现在“硬”地层和“软”地层的井中,导波都是高度频散的,其最大相速度等于地层的横波速度,其截止频率低于对称模式的伪瑞利波的截止频率;在低频（2-3kHz）和长源距（3-4m）的条件下,由非对称的声源（如声偶极子）所产生的微地震图中,初至信号是以横波速度传播的,而以纵波速度传播的信号被抑制。本文的结果对研制横波速度测井仪是有意义的。     

Seismic earth pressures on rigid and flexible retaining walls

While limiting-equilibrium Mononobe–Okabe type solutions are still widely used in designing rigid gravity and flexible cantilever retaining walls against earthquakes, elasticity-based solutions have been given a new impetus following the analytical work of Veletsos and Younan [23]. The present paper develops a more general finite-element method of solution, the results of which are shown to be in agreement with the available analytical results for the distribution of dynamic earth pressures on rigid and flexible walls. The method is then employed to further investigate parametrically the effects of flexural wall rigidity and the rocking base compliance. Both homogeneous and inhomogeneous retained soil is considered, while a second soil layer is introduced as the foundation of the retaining system. The results confirm the approximate convergence between Mononobe–Okabe and elasticity-based solutions for structurally or rotationally flexible walls. At the same time they show the beneficial effect of soil inhomogeneity and that wave propagation in the underlying foundation layer may have an effect that cannot be simply accounted for with an appropriate rocking spring at the base.     

Plain strain soil–structure interaction model for a building supported by a circular foundation embedded in a poroelastic half-space

A simple theoretical model for soil–structure interaction in water saturated poroelastic soils is presented, developed to explore if the apparent building–foundation–soil system frequency changes due to water saturation. The model consists of a shear wall supported by a rigid circular foundation embedded in a homogenous, isotropic poroelastic half-space, fully saturated by a compressible and inviscid fluid, and excited by in-plane wave motion. The motion in the soil is governed by Biot's theory of wave propagation in fluid saturated porous media. Helmholtz decomposition and wave function expansion of the two P-wave and the S-wave potentials is used to represent the motion in the soil. The boundary conditions along the contact surface between the soil and the foundation are perfect bond (i.e. welded contact) for the skeleton, and either drained or undrained hydraulic condition for the fluid (i.e. pervious or impervious foundation). For the purpose of this exploratory analysis, the zero stress condition at the free surface is relaxed in the derivation of the foundation stiffness matrix, which enables a closed form solution. The implications of this assumption are discussed, based on published comparisons for the elastic case. Also, a closed form representation is derived for the foundation driving forces for incident plane (fast) P-wave or SV wave. Numerical results and comparison with the full-scale measurements are presented in the companion paper, published in this issue.     

Embedded foundation in layered soil under dynamic excitations

The critical step in the substructure approach for the soil–structure interaction (SSI) problem is to determine the impedance functions (dynamic-stiffness coefficients) of the foundations. In the present study, a computational tool is developed to determine the impedance functions of foundation in layered soil medium. Cone frustums are used to model the foundation soil system. Cone frustums are developed based on wave propagation principles and force-equilibrium approach. The model is validated for its ability to represent the embedded foundation in layered medium by comparing the results with the rigorous analysis results. Various degrees of freedom, such as, horizontal, vertical and rocking are considered for this study.     

孔隙地层轻质水泥套管井中偶极弯曲波的激发与传播特性

针对软地层套管井中弯曲波频散移向高频,以致超出了现行偶极声波测井仪器的激发频带的问题,本文在井外为孔隙地层时,采用Biot模型对套管井多极源激发的声波场进行了理论推导,对套管井偶极弯曲波的频散特性进行了数值模拟,考察了高密度水泥(或快速水泥)和低密度水泥(或轻质水泥)情况下的弯曲波频散,分析了不同水泥环对弯曲波主频散区和激发谱偏移的影响,重点对轻质水泥套管井中偶极源激发的模式波频散与激发谱及临界折射P、S波激发特性进行了研究,并考察了渗透率、孔隙度对弯曲波频散及衰减的影响,利用实轴积分法计算了偶极声源激发的时域全波波形.分析对比结果表明,在快速水泥情况下弯曲波频散曲线随着地层特征横波速度的减小会迅速向高频移动,随着快速水泥向轻质水泥变化,弯曲波频散曲线向高频移动将会减缓,对于特征横波速度低于1400m·s~(-1)的软地层,偶极弯曲波基础模式主频散区(或截止频率)可以由原来13kHz向低频移动至4kHz;在轻质水泥套管井中,无论是硬地层还是软地层,在目前偶极声波测井仪器声源主频激发下,接收波形中弯曲波均占主导地位.     

Antiplane SH-deformations near a surface rigid foundation above a subsurface rigid circular tunnel

The problem on the dynamic response of a rigid embedded foundation in the presence of an underground rigid tunnel and subjected to excitation of incident anti-plane SH waves is analyzed. By using the exact analytical solution for the two-dimensional SH-wave propagation in and around both the surface rigid foundation and subsurface rigid tunnel, those aspects of the resulting ground motions that are of special interest and importance for seismic resistant design in earthquake analyses have been examined. The computed amplitudes of the resulting periodic ground motions display a very complicated wave-interference between the surface foundation and underground tunnel that lead to observed standing wave patterns, together with abrupt changes in the wave amplitudes and large amplification of the incident motions. Supported by: National Science Foundation grant CMS 97-14859     

横向各向同性地层(VTI)井孔声弹效应对弯曲波的影响

本文发展了建立在地层参考状态为各向同性介质假定下的现行井孔声弹性理论,就井外为横向各向同性面与井轴垂直的、具有9个独立三阶弹性模量的横向各向同性介质(VTI井况),水平面内受双轴应力作用下给出了一个简洁的与井内压力、应力差、应力和以及多极源偏振方位角有关的井孔弯曲波声弹公式,并且导出了平面纵、横波速度的声弹公式.数值考察了弯曲波速度之改变量的灵敏系数随频率的变化、受井外水平双轴应力作用时两种偏振的偶极弯曲波频散曲线以及对应不同方位径向偏振的横波速度.研究结果表明弯曲波声弹公式与5个二阶弹性系数以及7个独立的三阶弹性模量有关;而且由内压引起的井孔弯曲波声弹性公式中的三阶弹性模量仅与6个独立的弹性模量有关.横向各向同性介质井孔弯曲波速度的交叉现象仍是判断地应力存在的标志;一个重要的认识是受双轴应力作用的弯曲波速度变化在低频区主要与c₁₄₄和c₁₅₅两个三阶弹性模量有关,而且此认识与径向偏振的平面横波一致.在缺乏足够实验条件情况下,对VTI情况,以c₁₄₄,c₁₅₅和c₁₂₃三个独立的量进行测量,然后可暂不考虑三阶弹性模量的各向异性,建立简化的应力反演公式.反之,如果已知地层的地应力信息,由简化的声弹公式可以反演三阶弹性模量c₁₄₄,c₁₅₅和c₁₂₃.     

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.     

Fluid–foundation interaction in the seismic response of gravity dams

The seismic response of a dam is strongly influenced by its interaction with the water reservoir and the foundation. The hydrodynamic forces in the reservoir are in turn affected by radiation of waves towards infinity, wave absorption at the reservoir bottom, and cross-coupling between the foundation below the dam and the reservoir bottom. The fluid–foundation interaction effect, i.e. the wave absorption along the reservoir bottom, can be accounted for by using either an approximate one-dimensional (1D) wave propagation model or a rigorous analysis of interaction between the flexible soil along the base and the water. The rigorous approach requires enormous computational effort because of (a) cross-coupling between the foundation of the dam and the soil below the reservoir and (b) frequency dependence of the boundary condition along the fluid-foundation interface. The analysis can be simplified by ignoring the cross-coupling and by using the approximate 1D wave propagation model. The effects of each of these two simplifications on the accuracy and computational efficiency of the procedure used for the seismic response analysis of a dam are examined. Analytical results are presented for the complex frequency-response functions as well as the time histories of the response of Pine Flat dam to Taft and E1 Centro ground motions.     

Effects of rainfall on soil–structure system frequency: Examples based on poroelasticity and a comparison with full-scale measurements

In this paper, a simple two-dimensional soil–structure interaction model, based on Biot's theory of wave propagation in fluid saturated porous media, is used to explain the observed increase of the apparent frequencies of Millikan library in Pasadena, California, during heavy rainfall and recovery within days after the rain. These variations have been measured for small amplitude response (to microtremors and wind excitation), for which Biot's linear theory is valid. The postulated hypothesis is that the observed increases in frequency are due to the water saturation of the soil. The theoretical model used to explore this hypothesis consists of a shear wall supported by a circular foundation embedded in a poroelastic half-space. This rigid foundation model may be appropriate only for the NS response of Millikan library. This paper presents results for the foundation stiffness, and for the system response for model parameters similar to those for Millikan library (located on alluvium with shear wave velocity of about 300 m/s). The foundation impedance matrix, foundation input motion and system response are compared for dry and fully saturated half-space, with permeable and impermeable foundation. The results show that for embedded foundations, the effects of saturation on the horizontal foundation stiffness are as significant as for the vertical stiffness, contrary to what has been known for surface foundations investigated by other authors. Further, the results suggest a 1–2% increase in system frequency of the first two modes of vibration, depending on the drainage condition along the foundation–soil interface. Such increases agree qualitatively with the observations.     

孔、裂隙并存地层中的声波测井理论及多极子声场特征

在油、气储层的勘探和开发中观察到的一个现象是储层岩石中普遍存在孔隙和裂隙.随着近年来孔、裂隙介质弹性波动理论的进展,我们可以将此理论应用于测井技术,以此来指导从声波测井中测量孔、裂隙地层的声学参数.本文计算了孔、裂隙地层里充流体井眼中的多极子声场,分析了声场随裂隙介质的两个主要参数(即裂隙密度和裂隙纵横比)的变化特征.井孔声场的数值计算表明裂隙密度可以大幅度地降低井中声波纵、横波的波速和振幅.随着裂隙密度的增加,在测井频段内也可以看到纵、横波速的频散现象(这种频散在孔隙地层中一般是观察不到的).本文还研究了多极子模式波 (即单极的Stoneley波、伪瑞利波以及偶极的弯曲波)随裂隙参数的变化特征.结果表明,这些模式波的振幅激发和速度频散都受裂隙密度的影响.裂隙密度越高影响越大.此外,裂隙还对模式波的传播造成较大的衰减.相对裂隙密度而言,裂隙纵横比是一个频率控制参数,它控制裂隙对声场影响的频率区间.本文的分析结果对裂缝、孔隙型地层的声波测井具有指导意义.     

Dynamic stiffness and kinematic response of single piles in inhomogeneous soil

The effect of soil inhomogeneity on dynamic stiffness and kinematic response of single flexural elastic piles to vertically-propagating seismic SH waves is explored. A generalized parabolic function is employed to describe the variable shear wave propagation velocity in the inhomogeneous stratum. A layered soil with piece-wise homogeneous properties is introduced to approximate the continuous inhomogeneity in the realm of a Beam-on-Dynamic-Winkler-Foundation model. The problem is treated numerically by means of a layer transfer-matrix (Haskell–Thompson) formulation, and validated using available theoretical solutions and finite-element analyses. The role of salient model parameters such as pile-head fixity conditions, pile-to-soil stiffness ratio, surface-to-base shear wave velocity ratio and rate of inhomogeneity is elucidated. A new normalization scheme for inertial and kinematic response of such systems is presented based on an average Winkler wavenumber. With reference to long piles in moderately inhomogeneous soils, results indicate that: (a) kinematic pile response is essentially governed by a single dimensionless frequency parameter accounting for pile-to-soil stiffness ratio, pile slenderness and soil inhomogeneity and (b) definition of a characteristic pile wavelength allows an approximate estimation of pile elastodynamic response for preliminary design or analysis. Issues related to domain discretization and Winkler moduli are discussed.     

地震波传播的哈密顿表述及辛几何算法

地震波传播过程本质上是能量在传播过程中逐步损耗直至殆尽的过程，而在实际应用中，常在无能量损耗假设下，用弹性波动方程或标量波动方程描述它.在哈密顿(Hamilton)体系表述下，地震波传播过程即为一个无限维的哈密顿系统随时间的演化过程.若不计能量损耗，波场演化过程实质上为一个单参数连续的辛变换，因而对应的数值算法应为辛几何算法.本文首先从地震波标量方程出发，给出哈密顿体系下地震波传播的表述，即任意两个时刻的波场是通过辛变换联系起来的.随后，把波场在时间和相空间离散化后，给出了用于波场计算的一些辛格式，如显式辛格式、隐式辛格式和蛙跳辛格式.并进一步讨论了有限差分格式和辛格式的异同.然后，应用显式辛格式和同阶的有限差分方法给出了同一理论速度模型下的波场和Marmousi速度模型下的单炮记录.数值结果表明，辛算法是一类可行的波场模拟的数值算法.在时间步长较小时，有限差分方法是辛算法的一个很好近似.文中的理论和方法，为地震波传播理论及实际应用研究提供了新的途径.     

Can ground densification improve seismic performance of the soil‐foundation‐structure system on liquefiable soils?

Over the past few decades, soil densification has been widely employed to reduce the liquefaction hazard or consequences on structures. The decision to mitigate and the design of densification specifications are typically based on procedures that assume free‐field conditions or experience. As a result, the influence of ground densification on the performance of structures and the key mechanisms of soil‐structure interaction remains poorly understood. This paper presents results of four centrifuge tests to evaluate the performance of 3‐ and 9‐story, potentially inelastic structures on liquefiable ground with and without densification. Densification was shown to generally reduce the net excess pore pressures and foundation permanent settlements (although not necessarily to acceptable levels), while amplifying the accelerations on the foundation. The influence of these demands on the performance of the foundation and superstructure depended on the structure's strength and dynamic properties, as well as ground motion characteristics. In addition, densification tended to amplify the moment demand at the beam and column connections, which increased permanent flexural deformations and P‐Δ effects (particularly on the heavier and weaker structure) that could have an adverse effect on foundation rotation. The experimental results presented aim to provide insight into the potential tradeoffs of ground densification, which may reduce foundation permanent settlement, but amplify shaking intensity that can result in larger foundation rotation, flexural drifts, and damage to the superstructure, if not considered in design. These considerations are important for developing performance‐based strategies to design mitigation techniques that improve performance of the soil‐foundation‐structure system in a holistic manner.     

A micromechanics-based analytical method for wave propagation through a granular material

This paper describes the creation and use of the stiffness and mass matrices of a granular sample to calculate its modes of vibration and natural frequencies, and obtain its transfer function. This is an important result as it is a necessary prerequisite for the application of linear systems theory to granular materials. The analytical solution for wave propagation through the material is then presented and simulations of bender element tests on an idealised two-dimensional granular sample are used as an example of the applicability of the method. Analysing wave propagation through sand and sandstone in this novel way that is explicitly particle-based can offer significant insights into phenomena such as wave dispersion. It can also offer insight into the interpretation of laboratory bender element tests, which is often challenging. Bender element tests are used to obtain the small-strain stiffness of soil for engineering design and the results generated by applying the analytical method proposed here give insight into the effects of sample aspect ratio and boundary conditions on the observed response. The model presented here allows for comparisons between the static stiffness of the sample and the stiffness value calculated through wave velocities, as obtained by a variety of methods used by experimentalists. This leads to the identification of the best procedure for extracting the small-strain stiffness value for the dry granular sample considered here.     

SH waves in a layer with temperature dependent properties

The paper deals with the problem of SH harmonic wave propagation in an elastic layer with temperature dependent properties. The shear modulus and mass density are linearly dependent on temperature. The layer is rested on a rigid foundation and the upper boundary is free of loadings. The boundary planes are kept at different constant temperatures. The wave velocity and amplitude of stresses are analysed.     

Dynamic response of a group of flexible foundations to incident seismic waves

The dynamic response of a finite number of flexible surface foundations subjected to harmonic incident Rayleigh or SH waves is presented. The foundations are assumed to be resting on an elastic half-space. The results show that the foundation stiffness has a marked effect on the vertical response, while there is only a minor effect on the horizontal displacements. In general, the dynamic response decreases with increasing foundation stiffness. In cases of Rayleigh wave incidence, the existence of an adjacent foundation generates a certain amount of horizontal response in the direction perpendicular to the incident wave and subsequently causes the system to undergo a torsional motion; while in cases of horizontally incident SH waves, a vertical response has been observed and its magnitude is comparable to the response in the direction of the incident wave.     

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.     

横向非均匀介质中弹性地震波传播的数值模拟

应用混合变量弹性动力学方程和线性常微分方程组的矩阵指数解法，将层状介质中广泛应用的弹性波传播矩阵解法推广至横向非均匀介质，给出了一种可计算复杂地质体中弹性波传播的广义传播矩阵数值解法。该方法可模拟任意震源及所产生的各种体波、面波，数值结果表明具有很高的计算精度。     

Seismic analysis of dam-foundation-reservoir system including the effects of foundation mass and radiation damping

One of the main concerns in using commercial software for finite element analyses of dam-foundation-reservoir systems is that the simplifying assumptions of the massless foundation are unreliable. In this study, an appropriate direct finite element method is introduced for simulating the mass, radiation damping and wave propagation effect in foundations of damfoundation-reservoir systems using commercial software ABAQUS. The free-field boundary condition is used for modeling the semi-infinite foundation and radiation damping, which is not a built-in boundary condition in most of the available commercial software for finite element analysis of structures such as ANSYS or ABAQUS and thus needs to be implemented differently. The different mechanism for modeling of the foundation, earthquake input and far-field boundary condition is described. Implementation of the free-field boundary condition in finite element software is verified by comparing it with analytical results. To investigation the feasibility of the proposed method in dam-foundation-reservoir system analysis, a series of analyses is accomplished in a variety of cases and the obtained results are compared with the substructure method by using the EAGD-84 program. Finally, the massed and massless foundation results are compared and it is concluded that the massless foundation approach leads to the overestimation of the displacements and stresses within the dam body.