Propagation of transient SH-waves in a dipping structure by finite- and boundary element methods

Summary The boundary and the finite element formulations for the equations of elasticity are presented and applied to the problem of propagation of transient SH-waves in dipping layers overlying a half-space. When the finite element formulation is used, appropriate boundary conditions are imposed on the additional boundary dividing the half-space into a finite and an infinite region. These conditions ensure the transmission of waves across this boundary. When the boundary element method is applied, it is necessary to satisfy the radiation conditions. Theoretical seismograms for the displacement on the surface of the half-space are presented. They show that, for a specific case, the agreement between the two methods is satisfactory. The results can be compared with those found by the exact method of generalized rays in order to check the validity of the finite and the boundary element methods for the specific problem studied in this paper.     

Seismic response analysis of alluvial valley with two dipping layers by a boundary element-finite element coupling method

A coupling process of time-domain boundary element method and finite element method is constituted for analysis of the dynamic response of an alluvial valley. Since the characteristics of these two methods are different from each other, the time increment of the boundary element method is set to be an integer multiple of the finite element method, and a combination of a linear interpolation of boundary element equations at two successive time steps with a finite element equation is proposed. Using this procedure, responses of alluvial valleys with two dipping layers are analysed for a vertically incident SV wave, and influence of the lower layer on these responses is investigated. It is found that the presence of the additional lower layer affects not only the reflected waves but also the surface waves. Furthermore, when an outcropping lower layer exists at the basin edges, the valley suffers an amplification of motion and an extreme increase in strains. These results show that a small difference in the layer at the basin edge is significant in analysing dynamic responses of an alluvial valley.     

Acoustic emission from a surface-breaking crack under cyclic loading

Summary An isotropic, homogeneous, elastic half-space is subjected to a uniform stress system composed of a constant tensile stress with a superimposed cyclic tensile stress, both parallel to the free surface. The cyclic stresses are assumed to generate a surface-breaking crack of length l(t) which propagates normal to the surface. The unloading of the crack faces generates acoustic emission, which is primarily composed of surface waves. The elastodynamic reciprocity theorem for time-harmonic waves is used to determine the radiated system of transient and steady state surface waves. Dedicated to Professor Franz Ziegler on the occasion of his 70th birthday     

层状弹性场地基岩斜入射地震动二维反演

提出了层状弹性场地基岩斜入射地震动二维反演的一个方法,方法利用土层和半空间的精确动力刚度矩阵,根据地表某一点的水平和竖向两个地震加速度记录,进行基岩地震加速度时程和入射角的反演,并以SV波二维斜入射为例对方法进行了验证。研究表明,根据水平和竖向两个地震加速度记录可以进行基岩地震加速度时程和入射角的准确反演。同时,对文献中仅根据地表水平加速度记录进行基岩面水平地震动的一维反演方法进行了误差分析,研究表明,一维反演因忽略了波型转换,在模型上存在一定的误差,且误差随着实际地震动斜入射角的增大而显著增大。     

局部褶皱对层状地基马蹄形孔洞散射的影响分析

水平岩层在构造作用下会产生局部褶皱,研究褶皱对层状地基马蹄形孔洞散射的影响,对地表结构地震安全性评价具有重要意义。基于子结构法建立了复杂场地散射问题的控制方程,将地震波散射问题的求解转化具有规则边界条件的层状地基(自由场)的动力刚度和波动响应的求解。通过Fourier变换和引入对偶变量,将波动方程转化为一阶常微分方程,采用精细积分算法对土层可实现高效合并,施加边界条件可得到内部节点的格林函数,进一步得到动力刚度。同时,采用精细积分算法代替原传递矩阵法的层间合并,可得到层状自由场的波动响应。这种改进传递矩阵法对土层厚度和层数没有任何限制。通过与文献中的结果对比,验证了方法的正确性,并分析了局部褶皱对层状地基中马蹄形孔洞散射场的影响。结果表明:局部褶皱对地表位移幅值的影响与入射波类型、入射波频率以及局部褶皱几何构造等因素均有关系;地表位移峰值受马蹄形孔洞和局部褶皱共同作用的影响,其影响特性与入射波类型无明显关系。     

Reflection of plane waves at the boundary of a pre-stressed elastic half-space subject to biaxial stretch and simple shear

This paper is concerned with the effect of finite pure homogeneous biaxial stretch together with simple shear deformation on the reflection from a plane boundary of elastic waves propagating in a half-space of incompressible isotropic elastic material. This generalizes the previous work in which, separately, pure homogeneous strain and simple shear were considered. For a special class of constitutive law, it is shown that an incident plane harmonic wave propagating in the considered plane gives rise to a surface wave in addition to a reflected wave for every angle of incidence. The amplitude of the surface wave may vanish at certain discrete angles depending on the state of stress, biaxial stretch and simple shear deformation and then specialized to recover results obtained previously. The amplitude of the reflected wave is independent of the pre-stress but does depend upon the magnitude of deformation under consideration. The dependence of the reflected and surface wave behavior on the angle of incidence, amount of shear strain, biaxial stretch and the state of stress is illustrated graphically.     

Transient elastic waves in a wedge-shaped layer

Summary The theory of generalized ray is applied to analyzing transient elastic waves in a layered half-space with non-parallel interface. The propagation, reflection and refraction of longitudinal (P-) and transverse (SV-) waves which are generated by a line source in the surface layer of a two layer model are considered, each of the two homogeneous and isotropic layers having different density and inverse of wave speeds. Generalized ray integrals for multi-reflected rays in the top layer are formulated by using two rotated coordinate systems, one for each interface, and are expressed in terms of local wave slowness along each interface. Through a series of transformations of the local slowness, all ray integrals are expressible in a common slowness variable. Special attention is given to wave mode changes during reflection. The arrival time of each ray is then determined from the stationary value of the phase function with common slowness of the ray integral. Arrivals of head waves corresponding to rays refracted at a fast bottom are calculated from proper branch points of the Cagniard-mapping.With 3 FiguresPresented by F. Ziegler at the 16 th IUTAM Congress ICDAM Lyngby, Denmark, August 19–25, 1984.     

Propagation of Love waves in an elastic layer with void pores

The paper presents a study of propagation of Love waves in a poroelastic layer resting over a poro-elastic half-space. Pores contain nothing of mechanical or energetic significance. The study reveals that such a medium transmits two types of love waves. The first front depends upon the modulus of rigidity of the elastic matrix of the medium and is the same as the love wave in an elastic layer over an elastic half-space. The second front depends upon the change in volume fraction of the pores. As the first front is well-known, the second front has been investigated numerically for different values of void parameters. It is observed that the second front is many times faster than the shear wave in the void medium due to change in volume fraction of the pores and is significant     

层状饱和半空间中沉积谷地对斜入射平面P1波的三维散射

在作者给出层状饱和场地三维精确动力刚度矩阵和层状饱和半空间中移动荷载动力格林函数基础上,采用间接边界元方法在频域内求解了层状流体饱和场地中沉积谷地对斜入射平面P1波的三维散射问题。该方法的特点在于虚拟移动均布荷载和斜线孔隙水压可以直接施加在沉积与层状饱和半空间交界面而不存在奇异性。该文通过与已有结果的比较验证了方法的正确性,并以均匀饱和半空间和弹性基岩上单一饱和土层中沉积谷地为例进行了数值计算分析。研究表明,沉积谷地对平面P1波的三维散射与二维散射之间存在本质差别,入射角度、孔隙率、饱和土层刚度和饱和土层厚度等参数对沉积谷地附近动力响应有着显著影响。     

An approximate secular equation of Rayleigh waves in an isotropic elastic half-space coated with a thin isotropic elastic layer

In this paper, we are interested in the propagation of Rayleigh waves in an isotropic elastic half-space coated with a thin isotropic elastic layer. The contact between the layer and the half-space is assumed to be welded. The main purpose of the paper is to establish an approximate secular equation of the wave. By using the effective boundary condition method, an approximate secular equation of fourth order in terms of the dimensionless thickness of the layer is derived. It is shown that this approximate secular equation has high accuracy. From the secular equation obtained, an approximate formula of third order for the velocity of Rayleigh waves is established.     

Interface waves along an anisotropic imperfect interface between anisotropic solids

First and second order asymptotic boundary conditions are introduced to model a thin anisotropic layer between two generally anisotropic solids. Such boundary conditions can be used to describe wave interaction with a solid-solid imperfect anisotropic interface. The wave solutions for the second order boundary conditions satisfy energy balance and give zero scattering from a homogeneous substrate/layer/substrate system. They couple the in-plane and out-of-plane stresses and displacements on the interface even for isotropic substrates. Interface imperfections are modeled by an interfacial multiphase orthotropic layer with effective elastic properties. This model determines the transfer matrix which includes interfacial stiffness and inertial and coupling terms. The present results are a generalization of previous work valid for either an isotropic viscoelastic layer or an orthotropic layer with a plane of symmetry coinciding with the wave incident plane. The problem of localization of interface waves is considered. It is shown that the conditions for the existence of such interface waves are less restrictive than those for Stoneley waves. The results are illustrated by calculation of the interface wave velocity as a function of normalized layer thickness and angle of propagation. The applicability of the asymptotic boundary conditions is analyzed by comparison with an exact solution for an interfacial anisotropic layer. It is shown that the asymptotic boundary conditions are applicable not only for small thickness-to-wavelength ratios, but for much broader frequency ranges than one might expect. The existence of symmetric and SH-type interface waves is also discussed.     

Reflection and refraction of plane waves at interface between two piezoelectric media

This paper analyses reflection and refraction of plane waves at a perfect interface between two anisotropic piezoelectric media. The equations of elastic waves, quasi-static electric field, and constitutive relationships for the piezoelectric media are derived. A solution based on the inhomogeneous wave theory is developed to address the inconsistency between the numbers of independent wave modes in the media and the numbers of interfacial boundary conditions to obtain accurate reflection and refraction coefficients in case of strong piezoelectric media, where all the elastic and electric continuity conditions across the interface are satisfied simultaneously. The study shows that there exist independent and zero energy wave modes satisfying the general Snell’s law and propagating along the interface for any incident wave angle. These waves can be treated as pseudo surface waves. It is further found that all the reflection/refraction waves including the pseudo surface waves obey the energy conservation law at the interface boundary. In addition, the analysis also reveals that the reflection and refraction elastic waves can turn into pseudo surface waves at some critical incident angles.     

Three-dimensional fundamental solutions for transient heat transfer by conduction in an unbounded medium,half-space,slab and layered media

This paper presents analytical Green's functions for the transient heat transfer phenomena by conduction, for an unbounded medium, half-space, slab and layered formation when subjected to a point heat source. The transient heat responses generated by a spherical heat source are computed as Bessel integrals, following the transformations proposed by Sommerfeld [Sommerfeld A. Mechanics of deformable bodies. New York: Academic Press; 1950; Ewing WM, Jardetzky WS, Press F. Elastic waves in layered media. New York: McGraw-Hill; 1957]. The integrals can be modelled as discrete summations, assuming a set of sources equally spaced along the vertical direction. The expressions presented here allow the heat field inside a layered formation to be computed without fully discretizing the interior domain or boundary interfaces.The final Green's functions describe the conduction phenomenon throughout the domain, for a half-space and a slab. They can be expressed as the sum of the heat source and the surface terms. The surface terms need to satisfy the boundary conditions at the surfaces, which can be of two types: null normal fluxes or null temperatures. The Green's functions for a layered formation are obtained by adding the heat source terms and a set of surface terms, generated within each solid layer and at each interface. These surface terms are defined so as to guarantee the required boundary conditions, which are: continuity of temperatures and normal heat fluxes between layers.This formulation is verified by comparing the frequency responses obtained from the proposed approach with those where a double-space Fourier transformation along the horizontal directions [Tadeu A, António J, Simões N. 2.5D Green's functions in the frequency domain for heat conduction problems in unbounded, half-space, slab and layered media. CMES: Computer Model Eng Sci 2004;6(1):43–58] is used. In addition, time domain solutions were compared with the analytical solutions that are known for the case of an unbounded medium, a half-space and a slab.     

Love waves on a half-space with a gradient piezoelectric layer by the geometric integration method

The propagation of Love waves on an elastic homogeneous half-space with a piezoelectric gradient covering layer is studied by the geometric integration method in this article. First, the state transfer equation of a Love wave is derived from the governing equations and constitutive relations. Then, the transfer matrix of the state vector is obtained by solving the state transfer equation of a Love wave and then the stiffness matrix is obtained. By combining transfer matrices and the stiffness matrices of the gradient covering layer and the homogeneous half-space, the total surface stiffness matrix of a Love wave is obtained. Lastly, the application of the electrically open circuit and short circuit conditions and mechanically traction-free conditions gives the frequency dispersive relation of a Love wave. For the gradient covering layer, the material constants at the bottom of the covering layer may be greater or smaller than that at the top of the covering layer. The two situations and three kinds of gradient profiles for each of these two situations are investigated. The numerical results show that the Love wave speed is sensitive to not only the material constants at the bottom and the top of the covering layer, but also the gradient profiles of the covering layer.     

Green's function for two-and-a-half dimensional elastodynamic problems in a half-space

 This paper presents an analytical solution, together with explicit expressions, for the steady state response of a homogeneous three-dimensional half-space subjected to a spatially sinusoidal, harmonic line load. These equations are of great importance in the formulation of three-dimensional elastodynamic problems in a half-space by means of integral transform methods and/or boundary elements. The final expressions are validated here by comparing the results with those obtained with the boundary element method (BEM) solution, for which the free surface of the ground is discretized with boundary elements.     

Time‐domain BEM for three‐dimensional site response analysis of topographic structures

This paper presents a formulation of a time‐domain three‐dimensional boundary element method for site response analysis of topographic structures. The boundary element algorithm that uses the presented time‐convoluted traction kernels is applied to site response analyses of topographic structures. The seismic responses of canyon and ridge subjected to incident P and S waves are analyzed to demonstrate the accuracy of the kernels and the applicability of the presented boundary element algorithm for site response analysis of topographic structures. Seismic response analyses of three‐dimensional Gaussian‐shaped ridges show that the three‐dimensional axisymmetric ridge has a more amplification potential compared with three‐dimensional non‐axisymmetric elongated and two‐dimensional ridges, if the ridge is impinged by incident waves with wavelength of about the ridge's width. Copyright© 2009 John Wiley & Sons, Ltd.     

Surface waves in a periodic two-layered elastic half-space

Summary The paper deals with surface waves propagating through a periodic two-layered elastic half-space. The analysis is performed on the basis of a homogenized model with microlocal parameters. The velocity of the surface wave is obtained as a function of geometric and dynamic properties of the subsequent layers. The numerical examples illustrating the variations of the surface wave velocity are presented.     

Propagation of axisymmetric waves on the surface of a cylindrical cavity in elastic medium

A nondamped axisymmetric mode that propagates in an elastic cylindrical waveguide representing an extended cavity with a circular cross section in an infinite homogeneous medium is described. The wave dispersion in this system is analyzed and the similarity with and differences from other elastic media with one boundary are considered, including an infinite round rod and the surface of a half-space (Rayleigh wave). It is shown that, for axisymmetric waves in the cavity, a boundary frequency dependent on the curvature radius always exists, below which the waves are evanescent. A physical interpretation of results is given.     

Dynamic response of rigid strip-foundations by a time-domain boundary element method

The dynamic response of rigid strip-foundations placed on or embedded in a homogeneous, isotropic, linear elastic half-space under conditions of plane strain to either external forces or obliquely incident seismic waves of arbitrary time variation is numerically obtained. The above mixed boundary-value problems are treated by the time domain boundary element method which is used in a step-by-step timewise fashion to provide the foundation response to a rectangular impulse. Numerical examples are presented in detail to demonstrate the use and importance of the proposed method. The method appears to be more advantageous than frequency domain techniques, because it provides the transient foundation response in a natural and direct way and can form the basis for extension to the non-linear case.     

位于弹性半空间上的理想流体层动力反应--平面SV波入射

根据弹性固体与理想流体动力学方程,导出了固-液介质交界面上波的透射与反射系数的理论公式,分析了平面SV波从弹性半空间入射到与理想流体层的交界面时,弹性半空间的弹性模量和密度,平面SV波的入射角和频率,以及流体的体积模量和密度对理想流体层动压力的影响。