覆水饱和双相介质圆弧场地P波激励下散射问题的理论解

采用波函数展开法对平面P波入射复杂水域地形的空间变异性地震动场进行研究,该水域地形具有覆水层、饱和双相介质、场地非平坦以及第二类分层(场地跨越分层界面)等属性.首先,依据地震波反射和透射特性推导直角坐标系下的自由波场分布;然后,根据场地属性并引入大圆弧法分析极坐标系下的含有待定系数的散射波场;进而,结合土-水分界面和饱和土层分界面边界条件,求解散射波场中的待定系数;最后,通过自由波场和散射波场得到覆水饱和双相介质圆弧场地波函数理论解.基于理论解,通过算例验证了理论推导的合理性及可靠性,分析了地表位移在不同入射条件下的差异性.结果 表明,相对于均匀介质,饱和双相介质会显著影响地表位移分布.此外,入射波频率和角度对地震地面运动特性也有较大的影响.     

圆弧形层状沉积谷地对入射平面SV波散射解析解

利用Fourier—Bessel级数展开法给出了圆弧形层状沉积谷地对入射平面SV波散射问题的一个解析解，并分析了谷中层状沉积排列顺序及相对刚度等因素对散射的影响．     

大延伸非均匀介质中地震波全弹性散射理论Ⅰ--弹性波单次散射理论

所描述的工作聚焦于大延伸非均匀介质中非均匀弹性地震波散射问题的研究.应用Born近似及等效源原理,推出了来自连续横向无界非均匀层的弹性散射波的通解.这一工作是解决大延伸非均匀介质的弹性地震波多次散射问题的基础.在上述通解的基础上,建立了适用于大延伸非均匀介质的全弹性散射理论.该理论可包容小尺度非均匀体、大延伸非均匀介质全弹性波单次弱散射理论及标量波单次弱散射理论,因此可视其为一个更为广义和统一的弱散射理论.     

大延伸非均匀介质中地震波全弹性散射理论I-弹性波单次散射理论

所描述的工作聚焦于大延伸非均匀介质中非均匀弹性地震波散射问题的研究。应用Born近似及等效源原理,推出了来自连续横向无界非均匀层的弹性散射波的通解。这一工作是解决大延伸非均匀介质的弹性地震波多次散射问题的基础。在上述通解的基础上,建立了适用于大延伸非均匀介质的全弹性散射理论。该理论可包容小尺度非均匀体,大延伸非均匀介质全弹性波单次弱散射理论及标量波单次弱散射理论,因此可视其为一个更为广义和统一的弱散射理论。     

曲面曲率对Rayleigh波传播特性的影响

对任意形状的均匀各向同性线弹性曲面物体,用 WKB～（1）方法求解了沿曲面传播的Rayleigh表面波的运动微分方程,同时考虑了波传播方向及其垂直方向曲面曲率对波的穿透性的影, 所获波动方程的势函数解答表明,在一般情况下垂直波传播方向的曲面曲率对波的穿透深度的影响是不容忽视的.进而以同种介质平面表面情况下的Rayleigh面波的传播特性为基准,给出了曲面曲率引起波数或波速变化的解析表达式.通过理论分析和数值算例,描述了曲面上Rayleigh面波传播行为的一些基本特征.     

海洋地震工程流固耦合问题统一计算框架

海底地震动的模拟以及海洋工程结构的地震反应分析中,涉及到海水、饱和海床、弹性基岩、结构之间的相互耦合.传统的方法分别采用声波方程描述理想流体、Biot方程描述饱和海床、弹性波方程描述基岩和结构,分别考虑相互之间的耦合,十分不便.本文基于理想流体、固体分别为饱和多孔介质的特殊情形(孔隙率分别为1和0),由饱和多孔介质的Biot方程可退化得到理想流体的声波方程和固体的弹性波方程.然后,以饱和多孔介质方程为基础,经集中质量有限元离散,考虑不同孔隙率的饱和多孔介质之间耦合的一般情形,建立了该耦合情形的求解方法.进一步论证了该一般情形的耦合计算方法可分别退化到流体与固体、流体与饱和多孔介质、固体与饱和多孔介质之间的耦合计算,从而将流体、固体、饱和多孔介质间的耦合问题纳入到统一计算框架,并编制了相应的三维并行分析程序.以P-SV波垂直入射时,半无限层状海水-饱和海床、海水-弹性基岩、海水-饱和海床-弹性基岩三种情形的动力分析为例,采用统一计算框架结合透射边界条件进行求解,并与传递矩阵方法得到的解进行对比,验证了该统一计算框架的有效性以及并行计算的可行性.      

气饱和与水饱和多孔介质的弹性波场数值模拟

基于Biot介质理论,对气饱和多孔介质与水饱和多孔介质中弹性波的传播进行了数值模拟.通过饱和多孔介质的一阶双曲型速度-应力弹性波波场分离方程,采用交错网格高阶有限差分法实现了气饱和与水饱和介质中的高精度数值模拟,并利用完全匹配层(PML)吸收边界来处理边界反射问题,取得了较好的效果.模拟实例表明,水饱和介质中快纵波速度要远高于气饱和介质中的速度,并且频散较小,衰减较弱;含水饱和夹层的气饱和介质对弹性波振幅、能量有减弱作用;且气固耦合作用比水固耦合作用要小,气饱和情况下流相与固相基本达到完全解耦.     

波在双孔隙与单孔隙介质界面上的反射与透射

基于Biot理论和双重孔隙介质理论研究了弹性波在双重孔隙介质与流体饱和单一孔隙介质 界面的反射和透射问题,在界面上假定裂缝孔隙流体相对于固体骨架的位移为零,推导了反 射系数和透射系数的计算公式,数值讨论了反射系数和透射系数随入射角和频率的变化关 系. 同时,讨论了双重孔隙介质中3种压缩波(P-1, P-2和P-3波)和一种剪切波(S波) 的频散和衰减特性.     

用动光弹方法研究隧道岩爆的爆破扰动机理

利用动态光弹性方法,以Lamb模型为依据,研究了不同荷载形式下,模型内应力波的传播、相互作用规律以及刻槽的拦截效应。实验结果与弹性动力学解答基本一致:(1)倾斜冲击荷载作用下,模型内传播的有P波、S波和Rayleigh波,其中P波和Rayleigh波对岩爆有突出贡献;(2)等差条纹图的分布规律,可以合理解释岩爆高发区段与开挖面的位置关系;(3)模型表面的刻槽,可以实现对Rayleigh波的大量拦截,说明隧道表面切槽的岩爆控制方法理论上是可行的。     

波由单相弹性介质向横观各向同性液体饱和多孔介质传播时的反射与透射

研究了弹性波从单相各向同性弹性介质向双相各向异性液体饱和介质的传播问题。基于Biot的孔隙介质理论，并考虑动态渗透率，对以任意角度，任意频率入射的弹性波，解析地导出了在界面上的透射角，透射波波速和衰减系数，给出了各波在界面上的能量分配，数值计算了反射，透射系数与入射角，孔隙度等的关系曲线。     

AN ANALYTICAL SOLUTION FOR THREE-DIMENSIONAL DIFFRACTION OF PLANE P-WAVES BY A HEMISPHERICAL ALLUVIAL VALLEY WITH SATURATED SOIL DEPOSITS

An analytical solution for the three-dimensional scattering and diffraction of plane P-waves by a hemispherical alluvial valley with saturated soil deposits is developed by employing Fourier-Bessel series expansion technique. Unlike previous studies, in which the saturated soil deposits were simulated with the single-phase elastic theory, in this paper, they are simulated with Biot's dynamic theory for saturated porous media, and the half space is assumed as a single-phase elastic medium. The effects of the dimensionless frequency, the incidence angle of P-wave and the porosity of soil deposits on the surface displacement magnifications of the hemispherical alluvial valley are investigated. Numerical results show that the existence of a saturated hemispherical alluvial valley has much influence on the surface displacement magnifications. It is more reasonable to simulate soil deposits with Biot's dynamic theory when evaluating the displacement responses of a hemispherical alluvial valley with an incidence of P-waves.     

Seismic responses of a hemispherical alluvial valley to SV waves： a three-dimensional analytical approximation

An analytical solution to the three-dimensional scattering and diffraction of plane SV-waves by a saturated hemispherical alluvial valley in elastic half-space is obtained by using Fourier–Bessel series expansion technique. The hemispherical alluvial valley with saturated soil deposits is simulated with Biot’s dynamic theory for saturated porous media. The following conclusions based on numerical results can be drawn: (1) there are a significant differences in the seismic response simulation between the previous single-phase models and the present two-phase model; (2) the normalized displacements on the free surface of the alluvial valley depend mainly on the incident wave angles, the dimensionless frequency of the incident SV waves and the porosity of sediments; (3) with the increase of the incident angle, the displacement distributions become more complicated; and the displacements on the free surface of the alluvial valley increase as the porosity of sediments increases.The project was supported by the National Natural Science Foundation of China (50478062 and 10532070) and Open Fund at the Key Laboratory of Urban Security and Disaster Engineering (Beijing University of Technology), Chinese Ministry of Education. The English text was polished by Keren Wang.     

SERIES SOLUTION FOR SCATTERING OF PLANE SH－WAVES BY MULTIPLE SHALLOW CIRCULAR－ARC CANYONS

ＳＥＲＩＥＳＳＯＬＵＴＩＯＮＦＯＲＳＣＡＴＴＥＲＩＮＧＯＦＰＬＡＮＥＳＨ－ＷＡＶＥＳＢＹＭＵＬＴＩＰＬＥＳＨＡＬＬＯＷＣＩＲＣＵＬＡＲ－ＡＲＣＣＡＮＹＯＮＳＦａｎｇＹｉｎｇｇｕａｎｇ（房营光）（Ｇｕａｎｇｄｏｎｇ｀ＩｎｓｔｉｔｕｔｅｏｆＴｅｃｈｎｏｌ...     

Theoretical analysis of generalized loadings and image forces in a planar magnetoelectroelastic layered half-plane

The problem of a planar transversely isotropic magnetoelectroelastic layered half-plane subjected to generalized line forces and edge dislocations is analyzed. The complete solutions consist only of the simplest solutions for an infinite magnetoelectroelastic medium with applied loadings. The physical meaning of this solution is the image method. It is shown that the explicit solutions include Green's function for originally applied singularities in an infinite medium and the other image singularities are induced to satisfy free surface and interface continuity conditions. The mathematical method used in this study provides an automatic determination for the locations and magnitudes of all image singularities. The locations and magnitudes of image singularities are dependent on the roots of the characteristic equation which is related to the material constants of the layered half-plane. With the aid of the generalized Peach-Koehler formula, the explicit expressions of image forces acting on dislocations are easily derived from the full-field solutions of the generalized stresses. Numerical results for the full-field distributions of stresses, electric fields, and magnetic fields in the layered half-plane medium are presented based on the analytical solutions. The image forces and equilibrium positions of one dislocation, two dislocations, and an array of dislocations are presented by numerical calculations and are discussed in detail.     

Overall constitutive description of symmetric layered media based on scattering of oblique SH waves

This papers investigates the scattering of oblique shear horizontal (SH) waves off finite periodic media made of elastic and viscoelastic layers. It further considers whether a Willis-type constitutive matrix (in temporal and spatial Fourier domain) may reproduce the scattering matrix (SM) of such a system. In answering this question the procedure to determine the relevant overall constitutive parameters for such a medium is presented. To do this, first the general form of the dispersion relation and impedances for oblique SH propagation in such coupled Willis-type media are developed. The band structure and scattering of layered media are calculated using the transfer matrix (TM) method. The dispersion relation may be derived based on the eigen-solutions of an infinite periodic domain. The wave impedances associated with the exterior surfaces of a finite thickness slab are extracted from the scattering of such a system. Based on reciprocity and available symmetries of the structure and each constituent layer, the general form of the dispersion and impedances may be simplified. The overall quantities may be extracted by equating the scattering data from TM with those expected from a Willis-type medium. It becomes evident that a Willis-type coupled constitutive tensor with components that are assumed independent of wave vector is unable to reproduce all oblique scattering data. Therefore, non-unique wave vector dependent formulations are introduced, whose SM matches that of the layered media exactly. It is further shown that the dependence of the overall constitutive tensors of such systems on the wave vector is not removable even at very small frequencies and incidence angles and that analytical considerations significantly limit the potential forms of the spatially dispersive constitutive tensors.     

多水平成层的非均匀介质圆弧形峡谷地下多点地震动模拟

开发了考虑圆弧场地水平多分层效应的地下多点地震动模拟可视化程序,验证其可靠性。首先,描述了给出的圆弧场地多分层效应的SH波动散射频域解,其本质是突破了传统均匀介质和分层不穿越峡谷前提的局限性;其次,利用从水平地表到峡谷表面再到地下土层的传递函数,通过两次传递得到了圆弧峡谷地下的目标功率谱和反应谱;进而,具体描述所导出的圆弧场地平-凹相关性的相干函数及其推导过程;最后,基于上述理论,开发了多层状的非均匀圆弧场地的多点地震动模拟程序,给出了模拟结果与目标功率谱、反应谱和相干函数的拟合对比,验证了该理论和程序的合理性和可靠性,便于所发展理论和方法的实用化和推广。     

圆弧形凹陷地形表面覆盖层对入射平面P波的影响

利用Fourier-Bessel级数展开法给出了表面具有覆盖层的圆弧形凹陷地形对入射平面P波散射问题的一个解析解，并利用该解分析了不同形状凹陷地形表面覆盖层刚度和厚度对入射P波的影响。结果表明，凹陷地形表面覆盖层的存在，即使厚度很薄，对入射P波的散射也具有很大覆盖，覆盖层刚度和厚度的变化可显著改善凹陷地形场地的动力特性。     

A modified saddle point method for predicting sound fields above a non-locally reacting porous medium

This paper examines an asymptotic analysis for predicting sound fields above a rigid-frame porous medium, the so-called non-locally reacting porous medium. Their solutions can be represented by a direct wave term, a reflected wave term and a diffraction wave term. Exact and analytical solutions are available for the direct wave and the reflected wave from a perfectly hard ground. In the contrary, the diffraction wave term can only be cast in an integral form that is amenable to approximate analysis. A modified saddle-point method is explored to evaluate the diffraction integral asymptotically. Three different types of non-locally reacting surfaces, which are an extended reaction, a hard-backed layer, and an impedance-backed layer, were considered. The sound fields above these porous surfaces have the same form but they are different by an augmented diffraction term in the solutions. The analytical formula for the total sound fields, which can be stated in a closed form, offer a physically interpretable solution comprising of a direct wave and ground reflected wave terms. This latter term can further be split into a specularly reflected plane wave and ground wave components. A series of numerical comparisons have been conducted to validate the asymptotic analyses presented in this study. It has been demonstrated that the overall sound fields can be predicted well by the formula for all incidence angles and for a wide range of non-locally reacting porous surfaces.     

Wave dispersion in randomly layered materials

Wave scattering in materials composed of two kinds of alternating layers with different elastic properties and randomly distributed thicknesses has been modeled. The general form of the dispersion equation is derived for the unbounded layered medium. It defines two basic macroscopic characteristics of the scattered wave: phase velocity and attenuation, which are explicit functions of wave frequency and microscopic parameters of the system: acoustic properties of the layers and stochastic characteristics of their thickness distributions. The analytical expressions are derived for three special cases: for long waves; for a periodic medium composed of layers with constant thicknesses and for random medium with uniform distribution of layer thicknesses. Special attention is paid to the analysis of the frequency dependence of the wave parameters. It was shown that the predictions of the model for long waves and for periodic medium are compatible with the results obtained in the literature.Moreover, comparison of theoretical results for frequency dependent wave parameters with numerical simulations of pulse transmission through the slab of the randomly layered medium shows good qualitative and quantitative agreement in wide frequency range.