层状半空间中周期分布凸起地形对平面SH波的散射

提出了一种新的以层状半空间中周期分布斜线荷载动力格林函数为基本解的间接边界元方法,研究了周期分布凸起地形对平面SH波的散射问题.方法将散射波场分解为凸起内部散射波场和凸起外部散射波场.凸起内部散射波场通过在凸起闭合边界上施加虚拟斜线荷载产生的动力响应来模拟,而凸起外部散射波场则通过在凸起与半空间交界面上施加虚拟周期分布斜线荷载产生的动力响应来模拟.周期分布斜线荷载动力格林函数的引入,使得本文方法仅需针对一个凸起进行边界单元的离散和求解,便可完成问题的求解,避免了通过截断无限边界求解而引入的误差,方法具有较高精度的同时显著降低了求解自由度.文中通过与已有结果的比较,验证了方法的正确性,并以均匀半空间和基岩上单一土层中周期分布凸起为例进行了数值计算分析.研究表明,凸起间距对凸起地形间的动力相互作用有着显著的影响,同时层状半空间中周期分布凸起地形对SH波的散射与均匀半空间情况也有着显著的差别.     

曲线坐标在弹性波散射研究中的应用──SH波对不等深度凹陷地形的散射

本文利用曲线坐标的分析方法研究ＳＨ波对凹陷地形散射的边值问题。首先给出了在曲线坐标系中波动方程解的逼近函数序列，并使此解满足由曲线坐标所描述的具有不同深度的凹陷地形上的边界条件，最后将待解的问题归结为对一组无穷代数方程组的求解．本文提供的方法有如下的特点：（１）可将通常求解凹陷地形对ＳＨ波散射的波动方程的混合边值问题（指在不同边界区域上给出不同的应力分量的边界条件）化为一应力边值问题；（２）本文所使用的曲线坐标，可描述具有不同深度与宽度比的拟半圆形凹陷地形散射问题，并给出一般解；（３）本文方法可推广用于求解Ｐ（－ＳＶ）波对凹陷地形的散射问题。作为算例，本文给出了ＳＨ波对具有不同深度与宽度比的凹陷地形散射的数值结果，并进行了分析比较。     

Application of the weighted residual method to diffraction by 2-D canyons of arbitrary shape: I. Incident SH waves

A solution for the two-dimensional scattering and diffraction of plane SH waves by canyons of arbitrary shape in an elastic half space is presented. The wave field for arbitrary geometry in this paper is computed numerically by the method of weighted residues (moment method). The wave displacement field computed by the present residual method for the case of a semi-circular canyon was shown to agree analytically and numerically with that computed by the exact closed form series solution. The same observations about ground amplifications, their dependence on frequencies and orientations of the incident waves, can be stated here for canyons of arbitrary shape as previously made for circular canyons.     

Scattering of SV waves by a canyon in a fluid-saturated, poroelastic layered half-space, modeled using the indirect boundary element method

The scattering of SV waves by a canyon in a fluid-saturated, poroelastic layered half-space is modeled using the indirect boundary element method in the frequency domain. The free-field responses are calculated to determine the displacements and stresses at the surface of the canyon, and fictitious distributed loads are then applied at the surface of the canyon in the free field to calculate the Green's functions for displacements and stresses. The amplitudes of the fictitious distributed loads are determined from the boundary conditions, and the displacements arising from the waves in the free field and from the fictitious distributed loads are summed to obtain the solution. The effects of fluid saturation, boundary conditions, porosity, and soil layers on the surface displacement amplitudes and phase shifts are discussed, and some useful conclusions are obtained. It is shown that the surface displacement amplitudes due to saturation and boundary conditions, different porosities, or the presence of a soil layer can be very dissimilar, and large phase shifts can be observed. The resulting wavelengths for an undrained saturated poroelastic medium are slightly longer than those for a drained saturated poroelastic medium; and are longer for a drained saturated poroelastic medium than those for a dry poroelastic medium. As porosity increases, the wavelengths become longer; and a layered half-space produces longer wavelengths than a homogeneous half-space.     

Surface motion of multiple alluvial valleys for incident plane SH-waves by using a semi-analytical approach

In this paper, the degenerate kernels and Fourier series expansions are adopted in the null-field integral equation to solve the exterior Helmholtz problems with alluvial valleys. The main gain of using degenerate kernels in integral equations is free of calculating the principal values for singular integrals by locating the null-field point exactly on the real boundary. An adaptive observer system is addressed to fully employ the property of degenerate kernels for circular boundaries in the polar coordinate. Image concept and technique of decomposition are utilized for half-plane problems. After moving the null-field point to the real boundary and matching the boundary conditions, a linear algebraic system is obtained without boundary discretization. The unknown coefficients in the algebraic system can be easily determined. The present method is treated as a “semi-analytical” solution since error only attributes to the truncation of Fourier series. Earthquake analysis for the site response of alluvial valley or canyon subject to the incident SH-wave is the main concern. Numerical examples including single and successive alluvial valleys are given to test our program. Limiting cases of a single canyon and two successive canyons are also addressed. Amplification of soft basin is also observed in this study. The validity of the semi-analytical method is verified. Our advantages, well-posed model, principal value free, elimination of boundary-layer effect and exponential convergence and mesh-free, by using the present method are achieved.     

Three-dimensional analysis of spatially varying ground motions around a uniform canyon in a homogeneous half-space

Department of Civil Engineering, University of California, Berkeley, CA 94720, U.S.A. A direct boundary element method to determine the three-dimensional seismic response of an infinitely-long canyon of arbitrary but uniform cross-section cut in a homogeneous viscoelastic half-space is presented. The seismic excitation is represented by P, SV, SH or Rayleigh waves at arbitrary angles with respect to the axis of the canyon. The accuracy of the procedure and implementing computer program is demonstrated by comparison with previous solutions for the limiting case of two-dimensional response, recently obtained three-dimensional response results for infinitely-long canyons, and three-dimensional boundary method solutions presented in this paper for finite canyons.     

A new boundary element solution to evaluate the geometric effects of the canyon site on the displacement response spectrum

This paper presents a step-by-step procedure using the three-dimensional boundary element approach to study the behavior of semi-circular canyons under seismic shear waves. The boundary element code TDASC allows utilization for various canyon geometries, evaluation of concurrent seismic waves and calculation of the ground motions on canyons due to an excitation at any arbitrary point of the incident field. Considering the widening ratio of the canyon(including prismatic, semi-prismatic and non-prismatic canyons), wave characteristics(wavelength, dimensionless period, direction) and maximum amplification pattern, the solution was applied to carry out a series of parametric studies. It was shown that canyon form can significantly affect the displacement amplification, especially at the points located on its edges. By increasing the wave dimensionless frequency(η 1), the amplification pattern becomes more complex. On the basis of the results from a variety of considered cases, a new expression has been presented for the limiting wavelength beyond which the widening of the canyon will not have a major effect on the displacement amplification. To verify the reliability of the proposed approach, the obtained results, expressed in terms of displacement amplitude, were compared with those from the available published literature and a reasonably good agreement was observed.     

Scattering of plane harmonic SH,SV, P and rayleigh waves by non-axisymmetric three-dimensional canyons: A wave function expansion approach

Scattering of elastic waves by three-dimensional canyons embedded within an elastic half-space is investigated by using a wave function expansion technique. The geometry of the canyon is assumed to be non-axisymmetric. The canyon is subjected to incident plane Rayleigh waves and oblique incident SH, SV and P waves. The unknown scattered wavefield is expressed in terms of spherical wave functions which satisfy the equations of motion and radiation conditions at infinity, but they do not satisfy stress-free boundary conditions at the half-space surface. The boundary conditions are imposed locally in the least-squares sense at several points on the surface of the canyon and the half-space. Through a comparative study the validity and limitations of two-dimensional approximations (antiplane strain and plane strain models) have been examined. It is shown that scattering of waves by three-dimensional canyons may cause substantial change in the surface displacement patterns in comparison to the two-dimensional models. These results emphasize the need for three-dimensional modelling of realistic problems of interest in strong ground motion seismology and earthquake engineering.     

饱和层状场地中任意凹陷地形对入射SV波的散射

利用饱和土层的精确动力刚度矩阵和动力格林(Green)函数,采用间接边界元法,在频域内求解了层状饱和场地中任意凹陷地形对入射SV波的散射问题.通过自由场反应分析,求得凹陷地形表面各点的位移和各单元的应力响应;同样计算了虚拟分布荷载的格林影响函数,求得相应的位移和应力响应;根据边界条件确定虚拟分布荷载,将自由场位移响应和...     

Ground motion at canyons of arbitrary shape under incident sh waves

A method for calculating the two-dimensional scattering of incident SH waves by canyons of arbitrary shape is presented. The problem is formulated in terms of a Fredholm integral equation of the first kind with the integration path outside the boundary. Point-source discretization and a least-squares scheme are used. Numerical results are compared with the known analytic solution for a semi-cylindrical canyon. Spatial variations of surface amplitudes are computed for triangular and half-cycle sinusoidal canyons as well.     

Impedance functions for three-dimensional foundations supported on an infinitely-long canyon of uniform cross-section in a homogeneous half-space

A direct boundary element procedure is presented to determine the impedance matrix for a three-dimensional foundation supported on an infinitely-long canyon of uniform cross-section cut in a homogeneous half-space. The uniform cross-section of the canyon permits analytical integration along the canyon axis leading to a series of two-dimensional boundary problems involving Fourier transforms of the full-space Green's functions. Solution of these two-dimensional boundary problems leads to a dynamic flexibility influence matrix which is inverted to determine the impedance matrix. The accuracy of the procedure is demonstrated by comparison with previous solutions for a surface-supported, square foundation and results obtained by a three-dimensional boundary element method (BEM) for a foundation of finite-width supported on an infinitely-long canyon. Compared with the three-dimensional BEM, the present method requires less computer storage and is more accurate and efficient. The foundation impedance matrix determined by this procedure can be incorporated in a substructure method for earthquake analysis of arch dams.     

风化半圆形河谷对柱面SH波的散射解析解

风化河谷地震效应对建在此类场地上工程结构(如大坝、桥梁)的动力稳定性具有潜在的威胁,为揭示河谷风化层对地面运动的影响,利用波函数展开法推导风化半圆形河谷对线源柱面SH波散射问题的解析解。计算不同震源位置条件下风化半圆形河谷的地表位移反应和地面运动放大因子,分析地面运动的幅值和形态,发现随着震源距离的增大,在河谷附近地震动的放大和衰减的交替更频繁。     

三维复杂山谷地形SV波垂直输入地震反应分析

本文基于显式有限元法研究了地震波垂直入射时三维复杂山谷地形对地震地面运动的影响,在数值分析中应用了三维化二维的解法和黏弹性人工边界的处理方法,实现了地震波垂直输入下三维复杂场地地震动数值模拟,并验证了该方法的合理性.以四川桃坪地区一山谷地形作为研究对象,基于地表高程数据分别建立了二维和三维场地模型,对比研究表明:在复杂地形情况下考虑二、三维模型时具有明显差异,三维模型能更真实地反映地形变化对地震动的影响,复杂地形条件下有必要考虑三维实际场地模型.本文对边界自由场的处理方法也可用于处理三维复杂场地地震动斜入射问题,为三维复杂地形场地地震效应研究提供参考.     

2.5D scattering of incident plane SV waves by a canyon in layered half-space

This paper presents a 2.5D scattering of incident plane SV waves by a canyon in a layered half-space by using the indirect boundary element method (IBEM). A free field response analysis is performed to provide the displacements and stresses on the boundary of the canyon where fictitious uniform moving loads are applied to calculate the Green’s functions for the displacements and stresses. The amplitudes of the loads are determined by the boundary conditions. The free field displacements are added to the fic...     

含峭壁V形峡谷对地震SH波散射的解析解

地表地形常引起地震动的局部放大,这是由于地震波传播至局部地形时产生了散射现象.本文利用波函数展开方法和区域匹配技术,提出了含峭壁V形峡谷对平面SH波散射问题的解析解,并进行了退化验证.通过频域内的参数分析,揭示了峭壁深度、入射波频率和角度等因素对峡谷场地地面运动的影响规律,发现上部峭壁会增强峡谷对地震动的地形放大效应.研究结果不仅为数值方法提供了验证基准,还可为含峭壁峡谷周边建筑物的抗震设计提供顺河向地震动输入.     

Anti-plane deformations around arbitrary-shaped canyons on a wedge-shape half-space: moment method solutions

The wave propagation behavior in an elastic wedge-shaped medium with an arbitrary shaped cylindrical canyon at its vertex has been studied. Nunerical computation of the wave displacement field is carried out on and near the canyon surfaces using weighted-residuals (moment method). The wave displacement fields are computed by the residual mcthod for the cases of elliptic, circular, rounded-rectangular and flat-elliptic canyons. The analysis demonstrates that thc resulting surface displacemcnt depends, as in similar previous analyses, on several factors including, but not limited, to the angle of thc wedge, thc geometry of thc vertex, the frcquencies of thc incident waves, the angles of incidence, and the material properties of the media. The analysis provides intriguing results that help to explain geophysical observations regarding the amplification of seismic energy as a function of site conditions.     

Lateral response of dams in semi-elliptical rigid canyons

An analytical closed-form solution is developed for the lateral response of earth and rockfill dams built in semi-elliptical canyons. The dam is idealized as a linearly hysteretic elastic body deforming only in shear, whereas the canyon is assumed to be rigid. The solution for the dam response is given in terms of prolate spheroidal radial and angular functions of the first kind and zero order. Results are presented for natural frequencies, modal displacement shapes, participation factors, and response to transient and steady-state harmonic base excitation for various dam length-to-height ratios. Comparisons are made of the effects of the length-to-height ratio and the canyon shape on the response of dams built in semi-elliptical and rectangular canyons. A subsequent study (Dakoulas, P. & Hsu, C.H., Response of earth dams in semi-elliptical flexible canyons to oblique SH waves, Report, Rice University, Houston, Texas, 1993) extends this model to a semi-elliptical canyon consisting of flexible elastic rock, subjected to obliquely incident harmonic SH waves.     

Scattering of plane sh waves by a semi-elliptical canyon

In this paper we analyse the two-dimensional scattering and diffraction of plane SH waves by a semi-elliptical canyon. The exact series solution of the problem, for general angle of incidence of the plane SH waves, has been used to examine the dependence of surface amplifications inside and near the canyon. The nature of ground motion has been found to depend on two key parameters: (a) The angle of incidence. (b) The ratio of the canyon width to the wave length of incident SH waves. For short incident waves surface displacement amplitudes change rapidly from one point to another, while for the long waves and shallow canyons displacement amplitudes display only minor departure from the uniform half-space amplification of 2. For shallow canyons and long incident waves, the angle of incidence introduces only minor changes into the overall behaviour of surface amplitudes. For deep canyons and nearly grazing incidences, a prominent shadow zone is realized behind the canyon.     

圆弧形凹陷饱和土场地对平面P波散射问题的解析解

以Biot 饱和多孔介质动力学理论为基础，利用Fourier Bessel级数展开法, 得到饱和多孔介质半空间中圆弧形凹陷地形对平面P波的散射问题的解析解. 数值计算给出地表位移幅值，分析了入射波波长、入射角、圆弧高宽比对地表位移幅值的影响，并与现存的在单相介质情况下得到的结论进行对比.     

Sediment delivery,resuspension, and transport in two contrasting canyon environments in the southwest Gulf of Lions

Multiple canyons incise the continental slope at the seaward edge of the continental shelf in the Gulf of Lions and are actively involved in the transfer of sediment from shelf to deep sea. Two canyons in the southwest region of the Gulf of Lions, Lacaze-Duthiers Canyon and Cap de Creus Canyon, were instrumented with bottom-boundary-layer tripods in their heads to evaluate the processes involved in sediment delivery, resuspension and transport. In both canyons, intense cold, dense-water flows carry sediment across the slope. In the Lacaze-Duthiers canyon head (located ∼35 km from the shoreline), dense-water cascading into the canyon was episodic. Currents were highly variable in the canyon head, and responded to interactions between the along-slope Northern Current and the sharp walls of the canyon. Inertial and other high-frequency fluctuations were associated with suspended-sediment concentrations of ∼5 mg/l. In Cap de Creus canyon head (located ∼14 km from the shoreline), downslope currents were higher in magnitude and more persistent than in Lacaze-Duthiers canyon head. Greater suspended-sediment concentrations (peaks up to 20 mg/l) were observed in Cap de Creus Canyon due to resuspension of the canyon seabed during dense-water cascading events. The similarities and contrasts between processes in these two canyon heads emphasize the importance of the interaction of currents with sharp canyon bathymetry. The data also suggest that cold, dense-water flows have more potential to carry sediment to the slope on narrow shelves, and may more efficiently transfer that sediment to the deep sea where a smooth transition between shelf and slope exists.