Response of a randomly inhomogeneous layer overlying a homogeneous half-space to surface harmonic excitation

The paper is concerned with the study of the effect of the randomness of material parameters on mean wave propagation in a semi-infinite viscoelastic medium. The medium considered in this paper is a configuration of a randomly non-homogeneous layer overlying a homogeneous half-space and is loaded harmonically on the top surface. The method used is that of Karal and Keller and is based on the idea of fundamental matrix and Bourret approximation. The integro-differential equation obtained in this paper is solved by the Laplace transform method. By using boundary and continuity conditions the mean wave solution is obtained. Numerical results show that the correlation functions, which introduce long-range interactions, can be the source of the wave amplification.     

Permanent deformations and strains in a shear building excited by a strong motion pulse

Examples of non-linear wave propagation in an elasto-plastic building are presented for excitation by pulses of strong ground motion characteristic of the near-field shaking near earthquake faults. Conditions that lead to the occurrence of permanent deformations in the building are investigated, and the amplitudes and wavelengths of incident pulses that lead to non-linear response are shown. Because the building can fail during the first passage of the incident wave pulse up and down the building (during a period that is shorter than the first natural period of the building), it is concluded that for the analysis and the design of structures in the near-field of earthquake shaking the wave propagation method of analysis must be used in place of the response spectrum method, which is based on the vibrational solution of the same governing equations.     

Harmonic wave propagation through viscoelastic heterogeneous media exhibiting mild stochasticity - I. Fundamental solutions

This work examines the propagation of time harmonic, horizontally polarized shear waves through a naturally occurring heterogeneous medium that exhibits viscous behaviour as well as random fluctuations of its elastic modulus about a mean value. As a first step, the governing equation, which is a heterogeneous Helmholtz equation, is solved using algebraic transformations and the relevant Green's function is obtained for two sets of boundary conditions, one corresponding to a finite depth layer and the other to an infinite layer. Viscous material behaviour is introduced by considering the depth-dependent elastic modulus to be a complex quantity. Subsequently, material stochasticity in the medium is handled through the perturbation approach by assuming that the elastic modulus has a small random fluctuation about its mean value. The final results are closed-form expressions for the mean value and covariance matrix of both the wave speed profile in the medium and the corresponding Green's function. In Part II, (Soil Dynam. Earth. Engng, 1996,15, 129-39), two examples concerning seismic wave propagation in soft topsoil and in sandstone serve to illustrate the methodology and comparisons are made with Monte Carlo simulations.     

The discrete wave number formulation of boundary integral equations and boundary element methods: A review with applications to the simulation of seismic wave propagation in complex geological structures

We review the application of the discrete wave number method to problems of scattering of seismic waves formulated in terms of boundary integral equation and boundary element methods. The approach is based on the representation of the diffracting surfaces and interfaces of the medium by surface distributions of sources or by boundary source elements, the radiation from which is equivalent to the scattered wave field produced by the diffracting boundaries. The Green's functions are evaluated by the discrete wave number method, and the boundary conditions yield a linear system of equations. The inversion of this system allows the calculation of the full wave field in the medium. We investigate the accuracy of the method and we present applications to the simulation of surface seismic surveys, to the diffraction of elastic waves by fractures, to regional crustal wave propagation and to topographic scattering.     

Elastic waves in nonhomogeneous media under 2D conditions: II. Numerical implementation

The purpose of this work is to present three methods of analysis for elastic waves propagating in two dimensional, elastic nonhomogeneous media. The first step, common to all methods, is a transformation of the governing equations of motion so that derivatives with respect to the material parameters no longer appear in the differential operator. This procedure, however, restricts analysis to a very specific class of nonhomogeneous media, namely those for which Poisson's ratio is equal to 0.25 and the elastic parameters are quadratic functions of position. Subsequently, fundamental solutions are evaluated by: (i) conformal mapping in conjunction with wave decomposition, which in principle allows for both vertical and lateral heterogeneities; (ii) wave decomposition into pseudo-dilatational and pseudo-rotational components, which results in an Euler-type equation for the transformed solution if medium heterogeneity is a function of one coordinate only; and (iii) Fourier transformation followed by a first order differential equation system solution, where the final step involving inverse transformation from the wavenumber domain is accomplished numerically. Finally, in the companion paper numerical examples serve to illustrate the above methodologies and to delineate their range of applicability.     

A METHOD OF TACKLING SOME SEISMIC PROBLEMS AND ITS APPLICATIONS *

The theory is explained and practical applications are shown for a numerical procedure in seismology. Particularly the problems concerning the generation of waves under the action of external pressures, and their propagation, in non-homogeneous, both elastic and absorbing media, have been carried out. These problems have been assumed mono-dimensional and refer to plane and spherical waves. The procedure is based on the solution, by means of series, of the wave differential equation, non-homogeneous, and with non-constant coefficients. It is a direct numerical method whose advantage is, mainly, the possibility of tackling, without great difficulties, problems regarding non-homogeneous elastic and absorbing media. On the contrary the methods which require the theoretical expression of the solutions by means of formulae, generally, present conceptual and numerical difficulties. As examples of application of this procedure, the following cases have been carried out by means of numerical calculations. a) Propagation of a wave, initially of symmetrical shape, in a viscoelastic medium; from the results it appears that the wave propagates without losing its symmetry, i.e. without sensible dispersion. A theoretical analysis has been carried out to justify this result, showing that the dispersion in viscoelastic media is noticeable only for relatively high frequencies. It seems that the practical absence of dispersion in field experiments do not exclude necessarily the viscoelastic character of absorption. b) Generation of plane waves under the effect of a uniform pressure distributed on the plane surface of a medium. The way the length and the shape of the generated wave depends not only on the type of pressure acting on the surface but also on the near surface impedance variations has been studied. c) Generation of a spherical wave under the action of a pressure in a spherical hole. The examples treated show how the length and shape of the wave depends on the radius of the hole. Particularly the frequencies of the wave spectrum are proportional to this radius, for a given type of pressure acting in the hole. The characteristics of this procedure would also permit the study of media for which the stress-strain relations are not univocal and linear (non linear absorption). This study, interesting for the wave propagation in rocks, is worth while to be carried out in a special paper.     

Wave propagation in semi-infinite media with topographical irregularities using Decoupled Equations Method

In this paper, a novel semi-analytical method, called Decoupled Equations Method (DEM), is presented for modeling of elastic wave propagation in the semi-infinite two-dimensional (2D) media which are involved surface topography. In the DEM, only the boundaries of the problem are discretized by specific subparametric elements, in which special shape functions as well as higher-order Chebyshev mapping functions are implemented. For the shape functions, Kronecker Delta property is satisfied for displacement function. Moreover, the first derivatives of displacement function with respect to the tangential coordinates on the boundaries are assigned to zero at any given node. Employing the weighted residual method and using Clenshaw–Curtis numerical integration, coefficient matrices of the system of equations are transformed into diagonal ones, which leads to a set of decoupled partial differential equations. To evaluate the accuracy of the DEM in the solution of scattering problem of plane waves, cylindrical topographical features of arbitrary shapes are solved. The obtained results present excellent agreement with the analytical solutions and the results from other numerical methods.     

Three-dimensional analysis of rows of piles as passive barriers for ground vibration isolation

This paper presents a theoretical study on the ground vibration isolation efficiency by multi-rows of piles as passive barrier in a three-dimensional context. Integral equations governing Rayleigh wave scattering are derived according to the Green's solution of Lamb problem. The integral equations are solved accurately and efficiently with an iteration technique. They are used to predict the complicated Rayleigh wave field that is generated by a number of irregular scatterers embedded in an elastic half-space solid. The method is verified with a numerical solution available in the literature for a simplified Rayleigh wave scattering problem. Passive isolation effectiveness of ground vibration with two or three rows of small piles is further studied in detail. Effects of relevant parameters on the effectiveness of vibration isolation are analyzed and presented.     

Representation of microscale ocean wavenumber spectrum correct to the second order

Up to now the finest synthetic aperture radar(SAR) sea surface image is still the state-of-art oceanimage provided by the satellite microwave remotesensing due to its high spatial resolution and dynamicremote sensing mechanism. Early in 1978, Valenzueladerived the analytical representation of SAR image(backscaterring cross section distribution) on the tem-porally stationary and locally homogeneous assump-tion[1] σ0 (θ)IJ = 4πkm cos4θ gIJ (θ) E0(k,0), (1) (1) …     

Fundamental solutions to Helmholtz's equation for inhomogeneous media by a first-order differential equation system

Helmholtz's equation with a variable wavenumber is solved for a point force through use of a first-order differential equation system approach. Since the system matrix in this formulation is non-constant, an eigensolution is no longer valid and recourse has to be made to approximate techniques such as series expansions and Picard iterations. These techniques can accommodate in principle any variation of the wavenumber with position and are applicable to scalar wave propagation in one, two and three dimensions, with the latter two cases requiring radial symmetry. As shown in the examples, good solution accuracy can be achieved in the near field region, irrespective of frequency, for the particular case examined, namely a wavenumber which increases (or decreases) as the square root of the radial distance from source to receiver. Finally, the resulting Green's functions can be used as kernels within the context of boundary element type solutions to study scalar wave scattering in inhomogeneous media.     

云南地区的背景弹性波场分析

利用云南省地震台网44个台站记录的2008年1-9月连续波形数据进行互相关计算,得到了台站间的格林函数,并获取对应的频散曲线,据此分析了该地区的区域背景弹性波场的来源及分布.研究发现,该地区的区域背景弹性波场有着明显的方向性,15s信号的总体传播方向是从东南向西北,也就是说主要来自云南省的东南方向,据此推测区域背景弹性...     

Anti-plane foundationless soil–structure interaction

This paper presents a closed-form wave function analytic solution of two-dimensional scattering and diffraction of anti-plane SH-waves by a two-dimensional foundationless structure that corresponds to a shear wall on an elastic half-space. A wave-function expansion method is used to solve this model by first prescribing a set of wave functions with undetermined coefficients and then assembling them together based on the stress and displacement boundary conditions on the surface between the structure and half space. This results in a set of infinite equations to be solved by truncating to a finite set. The amplitudes and residuals of the displacement and stress distributions around the structure and nearby ground surface will be discussed carefully. While the solution is analytical, the computation of the numerical results involves the evaluation of complicated integrals. This analytic solution will be helpful to the understanding of propagation of seismic or other stress waves within the superstructure(s) undergoing earthquakes or other blast loads.     

基于Kelvin-Voigt黏弹性骨架的含非饱和流体孔隙介质BISQ模型(英文)

本文综合考虑了在波传播过程中孔隙介质的三种重要力学机制——"Biot流动机制一squirt流动机制-固体骨架黏弹性机制",借鉴等效介质思想,将含水饱和度引入波动力学控制方程,并考虑了不同波频率下孔隙流体分布模式对其等效体积模量的影响,给出了能处理含粘滞性非饱和流体孔隙介质中波传播问题的黏弹性Biot/squirt(BISQ)模型。推导了时间-空间域的波动力学方程组,由一组平面谐波解假设,给出频率-波数域黏弹性BISQ模型的相速度和衰减系数表达式。基于数值算例分析了含水饱和度、渗透率与频率对纵波速度和衰减的影响,并结合致密砂岩和碳酸盐岩的实测数据,对非饱和情况下的储层纵波速度进行了外推,碳酸盐岩储层中纵波速度对含气饱和度的敏感性明显低于砂岩储层。     

Viscoelastic representation of surface waves in patchy saturated poroelastic media

Wave-induced flow is observed as the dominated factor for P wave propagation at seismic frequencies. This mechanism has a mesoscopic scale nature. The inhomogeneous unsaturated patches are regarded larger than the pore size, but smaller than the wavelength. Surface wave, e.g., Rayleigh wave, which propagates along the free surface, generated by the interfering of body waves is also affected by the mesoscopic loss mechanisms. Recent studies have reported that the effect of the wave-induced flow in wave propagation shows a relaxation behavior. Viscoelastic equivalent relaxation function associated with the wave mode can describe the kinetic nature of the attenuation. In this paper, the equivalent viscoelastic relaxation functions are extended to take into account the free surface for the Rayleigh surface wave propagation in patchy saturated poroelastic media. Numerical results for the frequency-dependent velocity and attenuation and the time-dependent dynamical responses for the equivalent Rayleigh surface wave propagation along an interface between vacuum and patchy saturated porous media are reported in the low-frequency range (0.1–1,000 Hz). The results show that the dispersion and attenuation and kinetic characteristics of the mesoscopic loss effect for the surface wave can be effectively represented in the equivalent viscoelastic media. The simulation of surface wave propagation within mesoscopic patches requires solving Biot’s differential equations in very small grid spaces, involving the conversion of the fast P wave energy diffusion into the Biot slow wave. This procedure requires a very large amount of computer consumption. An efficient equivalent approach for this patchy saturated poroelastic media shows a more convenient way to solve the single phase viscoelastic differential equations.     

Finite element equations and numerical simulation of elastic wave propagation in two—phase anisotropic media

Based on Biot theory of two-phase anisotropic media and Hamilton theory about dynamic problem,finite element equations of elastic wave propagation in two-phase anisotopic media are derived in this paper.Numerical solution of finite element equations is given.Finally,properties of elastic wave propagation are observed and analyzed through FEM modeling.     

双相各向异性介质中弹性波传播有限元方程及数值模拟

基于Biot双相各向异性介质理论和动态问题的哈密顿原理,推导出任意双相各向异性介质中弹性波传播的有限元方程,并给出双相各向异性介质中弹性波有限元方程的数值解法．最后进行有限元法的数值模拟,对双相各向异性介质中弹性波传播特征进行了模拟与分析．      

Green''''s functions for uniformly distributed loads acting on an inclined line in a poroelastic layered site

Based on one type of practical Biot's equation and the dynamic-stiffness matrices of a poroelastic soil layer and half-space, Green's functions were derived for uniformly distributed loads acting on an inclined line in a poroelastic layered site. This analysis overcomes significant problems in wave scattering due to local soil conditions and dynamic soil-structure interaction. The Green's functions can be reduced to the case of an elastic layered site developed by Wolf in 1985. Parametric studies are then carried out through two example problems.     

Analytical representation of Green’s functions of the Biot equations for elastic waves in a homogeneous two-phase medium

The system of Biot vector equations in the frequency space includes two elliptic-type vector partial differential equations with unknown displacement vectors in the solid and liquid phases. Considering the Biot equations, alongside with Pride’s equations, the key approaches to the theoretical study of the elastic waves in the two-phase fluid-saturated media, the author suggests an analytical solution for the inhomogeneous Biot equations in the frequency space, which is reduced to finding its fundamental solution (Green’s function). The solution of this problem consists of solutions for two systems of Biot equations. In the first system, only the first equation is inhomogeneous, while in the second system, only the second equation is inhomogeneous and, as it is shown, its right-hand side is exclusively a potential function. The fundamental solution of the full system of inhomogeneous Biot equations (in which both equations are inhomogeneous) is represented in the form of Green’s matrix-tensor, for the scalar elements of which the analytical relations are presented. The obtained formulas describing the elastic displacements of both the solid and liquid phases reflect three wave types, namely, compressional waves of the first and the second kind (the fast and the slow waves, respectively) and shear waves. Similar terms (those describing the same type of the elastic waves in the solid and liquid phases) in the expressions for Green’s functions are linked with each other through the coefficient that links the components of the displacement vectors of the solid and liquid phases corresponding to the given wave type.

     

Vibration isolation by trenches in continuously nonhomogeneous soil by the BEM

Vibration isolation of structures from ground-transmitted waves by open trenches in isotropic, linearly elastic or viscoelastic soil with a shear modulus varying continuously with depth is numerically studied. Both an exponential and a linear shear modulus variation with depth are used in this work. Waves produced by the harmonic motion of a rigid surface machine foundation are considered. The problem is solved by the frequency domain boundary element method employing the Green's function of Kausel-Peek-Hull for a thin layered half-space. Thus only the trench perimeter and the soil-foundation interface need essentially to be discretized. The proposed methodology is first tested for accuracy by solving two Rayleigh wave propagation problems in nonhomogeneous soil with known analytical solutions and/or for which experimental results are available. Then the method is applied to vibration isolation problems and the effect of the inhomogeneity on the wave screening effectiveness of trenches is studied.