Time-harmonic boussinesq problem for a continuously non-homogeneous soil

An analytical solution is presented for the response of a non-homogeneous, compressible, elastic half-space to a time-harmonic vertical point load on its surface. The shear modulus is assumed to increase continuously with depth. The model is chosen so as to describe uniformly deposited cohesionless soils. Expressions for displacements and stresses in the interior of the half-space medium are derived by means of Hankel transforms and contour integration. Selected numerical results are presented to demonstrate the influence of non-homogeneity. Finally, some effects are pointed out to be used in connection with vibration tests for subsoil investigation.     

Dynamic analysis of a layered half-space subjected to moving line loads

Using a thin-layer method enhanced by continued-fraction absorbing boundary conditions, dynamic responses of a layered half-space subjected to a series of constant and time-harmonic line loads moving at a constant speed are studied. The thin-layer method for moving line loads is formulated for plane-strain as well as antiplane-shear conditions and is verified by comparison of computed responses of a homogeneous half-space subjected to a single constant load on its surface against available analytical solutions. Next, time-harmonic loads on a homogeneous half-space are examined. The study continues with both constant and time-harmonic loads on a layered half-space. Finally, multiple constant and time-harmonic loads are considered. The formulation and results demonstrate the effectiveness and versatility of the method in problems of dynamic response of layered media to moving loads.     

Time harmonic point load and dynamic contact problem of contacting fluid and poroelastic half-spaces

The dynamic response of contacting fluid and fluid-saturated poroelastic half- spaces to a time-harmonic vertical point force or a point pore pressure is investigated. The solutions are formulated using the boundary conditions at the fluid-porous medium interface. The point load solutions are then used to solve the dynamic problem of the vertical vibration of a rigid disc (both permeable and impermeable discs are included) on the surface of the poroelastic half-space. The contact problems are solved by integrating the point force and point pore pressure solutions over the contact area with unknown discontinuous force and pore pressure distributions, which are determined from the boundary conditions. The solutions are expressed in terms of dual integral equations, which are converted to Fredholm integral equations of the second kind and solved numerically. Selected numerical results for the vertical dynamic compliance coefficient for the cases with or without fluid overlying the poroelastic half-space are presented to show the effects of the fluid. The influence of the permeability condition of the disc on the compliance of the poroelastic half-space is investigated. The displacement, vertical stress, pore pressure in the poroelastic half-space and water pressure in the fluid half-space are also examined for different poroelastic materials and frequencies of excitation. The present results are helpful in the study of the dynamic response of foundations on the seabed under seawater.     

Parameters determination for a linearly inhomogeneous half-space using characteristics of the time-harmonic surface waves

The paper presents a study of time-harmonic surface waves in a linearly inhomogeneous half-space. The study is based on the solution of that problem for an arbitrary (from 0 to 1/2) value of Poisson's ratio. Vertical vibrations due to a vertical harmonic force, which at large distances from the force represent Rayleigh-type waves, and transverse horizontal vibrations due to a horizontal force, which at large distances form waves of Love's type, are considered in detail. Material damping is taken into consideration. Inhomogeneity significantly affects relationships connecting wave characteristics and the frequency of vibration, and it is shown in the paper how this fact can be used for determining material properties (surface shear modulus, degree of inhomogeneity, damping ratio) with the help of experimental results concerning wave propagation over the surface of the half-space. It is shown that for forced waves the relationship between the wave phase angle and distance can significantly differ from a straight line, i.e. the wave number varies with distance. Therefore, it is desirable to relate experimental and theoretical results to such parts of wave propagation line, which correspond to same phase angle intervals. Copyright © 2000 John Wiley & Sons, Ltd.     

层状半空间中洞室对平面SH波的放大作用

利用间接边界元法,求解了弹性层状半空间中洞室对入射平面SH波的放大作用问题,并以基岩上单一土层中洞室对入射平面SH波的放大作用为例进行了数值计算分析。本文模型的特点之一是考虑了层状场地的动力特性,因而更接近于实际工程;特点之二是计算精度非常高。研究表明,层状半空间中洞室对波的放大作用与均匀半空间中情况有着本质的差别;层状半空间中洞室附近地表动力响应由土层动力特性和洞室对波的散射二者共同决定。土层动力特性不仅影响洞室附近地表位移的幅值,还会影响地表位移的频谱。在土层的前几阶共振频率附近,随着基岩与土层剪切波速比的增大,土层的影响随之增大,而随着土层厚度的增加,土层的影响随之减小,并逐渐趋于均匀半空间情况。     

Estimation of Elastic Moduli in a Compressible Gibson Half-space by Inverting Rayleigh-wave Phase Velocity

A Gibson half-space model (a non-layered Earth model) has the shear modulus varying linearly with depth in an inhomogeneous elastic half-space. In a half-space of sedimentary granular soil under a geostatic state of initial stress, the density and the Poisson’s ratio do not vary considerably with depth. In such an Earth body, the dynamic shear modulus is the parameter that mainly affects the dispersion of propagating waves. We have estimated shear-wave velocities in the compressible Gibson half-space by inverting Rayleigh-wave phase velocities. An analytical dispersion law of Rayleigh-type waves in a compressible Gibson half-space is given in an algebraic form, which makes our inversion process extremely simple and fast. The convergence of the weighted damping solution is guaranteed through selection of the damping factor using the Levenberg-Marquardt method. Calculation efficiency is achieved by reconstructing a weighted damping solution using singular value decomposition techniques. The main advantage of this algorithm is that only three parameters define the compressible Gibson half-space model. Theoretically, to determine the model by the inversion, only three Rayleigh-wave phase velocities at different frequencies are required. This is useful in practice where Rayleigh-wave energy is only developed in a limited frequency range or at certain frequencies as data acquired at manmade structures such as dams and levees. Two real examples are presented and verified by borehole S-wave velocity measurements. The results of these real examples are also compared with the results of the layered-Earth model.     

垂直地震力作用下弹性半空间上扁壳基础的解析解

用布希涅斯克定义的弹性半空间内的垂直位移包括两项积分,除了积分号前面系数的差别之外,第一项积分是单层位势而第二项积分为双层位势。若扁壳基础是正高斯曲率的几何曲面,则壳底与半空间表面间的挤压强度就是半空间表面作用的分布垂直荷载。当越过边界时,双层势位的函数值和单层势位的法向导数值发生跳跃。利用这些性质,本文得出布希涅斯克积分的反演公式,从而避开要求解偏微分—积分方程组的巨大数学困难而易于得出解析解。以椭园抛物面扁壳为例说明本文方法的应用。     

THE INSENSITIVITY OF REFLECTED SH WAVES TO ANISOTROPY IN AN UNDERLYING LAYERED MEDIUM1

Propagation in the plane of mirror symmetry of a monoclinic medium, with displacement normal to the plane, is the most general circumstance in anisotropic media for which pure shear-wave propagation can occur at all angles. Because the pure shear mode is uncoupled from the other two modes, its slowness surface in the plane is an ellipse. When the mirror symmetry plane is vertical the pure shear waves in this plane are SH waves and the elliptical SH sheet of the slowness surface is, in general, tilted with respect to the vertical axis. Consider a half-space of such a monoclinic medium, called medium M, overlain by a half-space of isotropic medium I with plane SH waves incident on medium M propagating in the vertical symmetry plane of M. Contrary to the appearance of a lack of symmetry about the vertical axis due to the tilt of the SH-wave slowness ellipse, the reflection and transmission coefficients are symmetrical functions of the angle of incidence, and further, there exists an isotropic medium E with uniquely determined density and shear speed which gives exactly the same reflection and transmission coefficients underlying medium J as does monoclinic medium M. This means that the underlying monoclinic medium M can be replaced by isotropic medium E without changing the reflection and transmission coefficients for all values of the angle of incidence. Thus no set of SH seismic experiments performed in the isotropic medium in the symmetry plane of the underlying half-space can reveal anything about the monoclinic anisotropy of that underlying half-space. Moreover, even when the underlying monoclinic half-space is stratified, there exists a stratified isotropic half-space that gives the identical reflection coefficient as the stratified monoclinic half-space for all angles of incidence and all frequencies.     

Quasi-static deformation of a viscoelastic half-space by a displacement dislocation

The correspondence principle is used to get the Laplace transformed solution of the problem of the quasi-static deformation of a viscoelastic half-space by a shear displacement dislocation from the corresponding static elastic results. The transformed solution is inverted for the Voigt and the Maxwell viscoelastic models. It is shown that, for a vertical dip-slip fault, the surface displacements for the viscoelastic case are identical with the elastic displacements. In the case of a vertical strike-slip fault, detailed numerical results are obtained for both a point source and a finite rectangular source. It is found that the results for the viscoelastic models differ significantly from the corresponding elastic results.     

Solution of dynamic Green׳s function for unsaturated soil under internal excitation

The closed form three-dimensional Green׳s function of a semi-infinite unsaturated poroelastic medium subjected to an arbitrary internal harmonic loading is derived, with consideration of capillary pressure and dynamic shear modulus varying with saturation. By applying the Fourier expansion techniques and Hankel integral transforms to the circumferential and radial coordinates, respectively, the general solution for the governing partial differential equations is obtained in the transformed domain. A corresponding boundary value problem is formulated. The integral solutions for the induced displacements, pore pressure and net stress are then determined considering the continuity conditions. The formulas are compared with the degenerated solution of saturated soils and confirmed. Numerical results reveal that the response of the unsaturated half-space depends significantly on the saturation by altering dynamic shear modulus to account for the effects of matric suction on soil stiffness. Slight differences between the results occur if only the saturation is taken into account. Moreover, a large source-depth results in a pronounced contribution to the reduction of surface displacement amplitudes. The analytical solutions concluded in the study offer a broader application to dynamic response associated with axi-symmetric and asymmetric conditions.     

竖向非均匀介质中的Love面波

本利用ＫＷＢＪ２（即几何近似）理论研究介质参数随深度作连续变化的竖向非均匀弹性半空间上覆盖一层厚度为Ｈ的元首中向同性的弹性介质时Ｌｏｖｅ面波的频散问题。给出了频散方程。中以剪切弹性模量和质量密度随深度呈抛物线变化的非均匀介质为例,给出其最低阶振型的频散曲线     

Lateral dynamic compliance of pile embedded in poroelastic half space

The time-harmonic response of a pile in a poroelastic half space and under lateral loadings is studied. By treating the pile as a one-dimensional structure and the half-space as a three-dimensional poroelastic continuum, the dynamic interaction between a pile and a poroelastic medium is formulated as a Fredholm integral equation of the second kind. Green's functions for a distributed lateral force field acting inside a poroelastic half space is an important ingredient of this paper. Numerical results for lateral dynamic compliance functions are presented to illustrate the dynamic characteristics of a pile in a poroelastic half space.     

Vertical vibration of an embedded rigid foundation in a poroelastic soil

This paper considers time-harmonic vertical vibration of an axisymmetric rigid foundation embedded in a homogeneous poroelastic soil. The soil domain is represented by a homogeneous poroelastic half space that is governed by Biot's theory of poroelastodynamics. The foundation is subjected to a time-harmonic vertical load and is perfectly bonded to the surrounding half space. The contact surface can be either fully permeable or impermeable. The dynamic interaction problem is solved by employing an indirect boundary integral equation method. The kernel functions of the integral equation are the influence functions corresponding to vertical and radial ring loads, and a ring fluid source applied in the interior of a homogeneous poroelastic half space. Analytical techniques are used to derive the solution for influence functions. The indirect boundary integral equation is solved by using numerical quadrature. Selected numerical results for vertical impedance of rigid foundations are presented to demonstrate the influence of poroelastic effect, foundation geometry, hydraulic boundary condition along the contact surface and frequency of excitation.     

Experimental study on the dynamic modulus of compacted loess under bidirectional dynamic load

The dynamic characteristics of compacted loess are of great significance to the seismic construction of the Loess Plateau area in Northwest China, where earthquakes frequently occur. To study the change in the dynamic modulus of the foundation soil under the combined action of vertical and horizontal earthquakes, a hollow cylindrical torsion shear instrument capable of vibrating in four directions was used to perform two-way coupling of compression and torsion of Xi'an compacted loess under different dry density and deviator stress ratios. The results show that increasing the dry density can improve the initial dynamic compression modulus and initial dynamic shear modulus of compacted loess. With an increase in the deviator stress ratio, the initial dynamic compression modulus increases, to a certain extent, but the initial dynamic shear modulus decreases slightly. The dynamic modulus gradually decreases with the development of dynamic strain and tends to be stable, and the dynamic modulus that reaches the same strain increases with an increasing dry density. At the initial stage of dynamic loading, the attenuation of the dynamic shear modulus with the strain development is faster than that of the dynamic compression modulus. Compared with previous research results, it is determined that the dynamic modulus of loess under bidirectional dynamic loading is lower and the attenuation rate is faster than that under single-direction dynamic loading. The deviator stress ratio has a more obvious effect on the dynamic compression modulus. The increase in the deviator stress ratio can increase the dynamic compression modulus, to a certain extent. However, the deviator stress ratio has almost no effect on the dynamic shear modulus, and can therefore be ignored.     

Response of a randomly inhomogeneous layered media to harmonic excitations

We present a unified formulation of an analytical method for evaluating the response of a random soil medium to surface or earthquake excitations. Specifically, we are interested in this study to the case of a horizontally stratified layered soil profile. Soil properties, mass density and shear modulus of each layer are modeled as spatial random fields. The soil profile is laying on a homogeneous half-space. Integro-differential equations are formulated and solved using the Laplace transform method. Numerical results are, firstly, obtained in terms of effects of soil properties on distributions of mean dimensionless displacement and stress with depth. Then, the amplification function of the surface layer of a soil profile is formulated and a parametric study is conducted to examine the effects of stochastic variations of mass density and shear modulus on the amplification function. We, lastly, computed the seismic response of a site located in the basin of Metidja in Algeria to base rock accelerations which were recorded at Keddara's station during the May 21, 2003 Boumerdes earthquake, in terms of free surface Fourier spectrum amplitudes and accelerations. Comparison between computed accelerations and recorded ones at Dar Elbeida site has proved the effectiveness of the approach.     

Stochastic dynamic analysis of a layered half-space

A stochastic thin-layer method is developed for the analysis of wave propagation in a layered half-space. A random field of shear moduli in the layered system is considered in terms of multiple correlated random variables. Expanding the random moduli and uncertain responses by means of Hermite polynomial chaos expansions and applying the Galerkin method in the spatial as well as stochastic domains, stochastic versions of thin-layer methods for a layered half-space in plane strain and antiplane shear are obtained. In order to represent the infinite half-space, continued-fraction absorbing boundary conditions are included in the thin-layer models of the half-space. Using these stochastic methods, dynamic responses of a layered half-space subjected to line loads are examined. Means, coefficients of variance, and probability density functions of the half-space responses with a varying correlation coefficient of the shear moduli are computed and verified by comparison with Monte Carlo simulations. It is demonstrated that accurate probabilistic dynamic analysis is possible using the developed stochastic thin-layer methods for a layered half-space.     

Out-of-plane (SH) soil-structure interaction: Semi-circular rigid foundation revisited

The model studied in this paper presents an extension of previous work for a shear wall on a semi-circular rigid foundation in an isotropic homogeneous and elastic half-space. The objective is to develop a soil-structure interaction model that can later be applied to the case of a flexible foundation. As shown in the Introduction below, Luco considered the case of a rigid foundation subjected to vertical incident plane SH waves, and Trifunac extended the solution for the same rigid foundation subjected to SH waves but for arbitrary angles of the incidence. In this paper, a new approach and model are presented for the same semi-circular rigid foundation with a tapered-shape (instead of rectangular) superstructure. The analytical expression for the deformation of the semi-circular rigid foundation below this tapered shear wall with soil-structure interaction in an isotropic homogeneous and elastic half-space is thus derived. Results are then compared with those of Trifunac discussed in the section below. This problem formulation can and will later be extended in the case of a flexible foundation that is semi-circular or arbitrarily shaped.     

Dynamic analysis of a rigid circular foundation on a transversely isotropic half-space under a buried inclined time-harmonic load

The dynamic analysis of a surface rigid foundation in smooth contact with a transversely isotropic half-space under a buried inclined time-harmonic load is addressed. By virtue of the superposition technique, appropriate Green׳s functions, and employing further mathematical techniques, solution of the mixed-boundary-value problem is expressed in terms of two well-known Fredholm integral equations. Two limiting cases of the problem corresponding to the static loading and isotropic medium are considered and the available results in the literature are fully recovered. For the static case, the results pertinent to both frictionless and bonded contacts are obtained and compared. With the aid of the residue theorem and asymptotic decomposition method, an effective and robust approach is proposed for the numerical evaluation of the obtained semi-infinite integrals. For a wide range of the excitation frequency, both normal and rotational compliances are depicted in dimensionless plots for different transversely isotropic materials. Based on the obtained results, the effects of anisotropy are highlighted and discussed.     

Cone models for nearly incompressible soil

Cones can be used to model soil in a unified strength-of-materials approach. For the vertical and rocking motions involving predominantly compressional-extensional deformation, the corresponding dilatational wave velocity tends to infinity for Poisson's ratio approaching 1/2. Based on the rigorous solution for the dynamic stiffness of a rigid disk for all frequencies, whereby the partition of the power among P-, S- and Rayleigh waves is also discussed, two special features are necessary for the vertical and rocking motions for nearly incompressible soil with Poisson's ratio between 1/3 and 1/2: (1) The appropriate wave velocity is selected as twice the shear wave velocity and not as the dilatational wave velocity; (2) A trapped mass which increases linearly with Poisson's ratio is introduced. The trapped mass can be assigned to the base mat, allowing the cone model to be constructed in the same way for all Poisson's ratios. The realization of cone models for surface foundations on a homogeneous half-space and on a layer on a flexible half-space and for embedded and pile foundations is addressed.     

Time-harmonic response of transversely isotropic multilayered half-space in a cylindrical coordinate system

With the aid of the analytical layer-element method, a comprehensive analytical derivation of the response of transversely isotropic multilayered half-space subjected to time-harmonic excitations is presented in a cylindrical coordinate system. Starting with the governing equations of motion and the constitutive equations of transversely isotropic elastic body, and based on the Fourier expansion, Hankel and Laplace integral transform, analytical layer-elements for a finite layer and a half-space are derived. Considering the continuity conditions on adjacent layers׳ interfaces and the boundary conditions, the global stiffness matrix equations for multilayered half-space are assembled and solved. Finally, some numerical examples are given to make a comparison with the existing solution and to demonstrate the influence of parameters on the dynamic response of the medium.