正交各向异性含液饱和多孔介质中应力波的传播

基于Biot的孔隙介质理论，研究了正交各向异性含液饱和多孔介质中应力波的传播特性．本文引入动态渗透率，导出了整个实频域内应力波传播的复特征方程及其解析解，给出了各种应力波成分的波速和衰减的解析表达武，计算了频散曲线和衰减曲线，并讨论了各类应力波之间的耦合关系及介质的各向异性对应力波传播的影响．     

用二次形函数薄层法分析弹性层状地基中的动力问题

薄层法是分析和模拟弹性波在层状介质中传播的一种半解析半数值方法。本文在土层垂直方向离散中利用Galerkin加权残值法推导出二次形函数薄层元的计算公式。采用薄层单元模拟半空间上的层状场地,模型底面用阻尼器边界或傍轴边界代替半空间。利用点源简谐荷载作用下的土层反应与其它数值分析方法的对比讨论薄层模型的设置指标。并用一次和二次两种形函数离散方法计算了层状地基中的面波频散曲线、圆形均布简谐荷载作用下弹性半空间的位移反应和半无限地基中单桩的竖向阻抗函数,分别讨论其计算精度。     

基于积分方程方法的弹性波多重散射计算

运用积分方程方法计算了含多个随机分布椭圆柱型孔洞的随机非均匀介质中相干波的速度和衰减系数,分析了这种介质的频散特性。首先,建立了散射位移场满足的积分方程,推导了单个椭圆柱孔洞的散射截面计算公式。接着分析了在含多个随机分布椭圆柱型孔洞的随机非均匀介质中弹性波的多重散射,给出在统计平均意义下的相干波的波速和衰减系数计算公式。然后用Matlab进行了编程,给出了一个数值算例,并将计算结果与波函数展开法进行了比较,分析了随机空隙介质的频散特征及其孔洞椭圆偏心率和材料空隙率的影响。     

非均匀吸收介质中地震波的广义传播矩阵解法

本文在复频域内,通过应用混合变量粘弹性波方程和线性常微分方程组的指数矩阵解法,给出了一种计算非均匀吸收介质中地震波传播的广义传播矩阵解法。该方法适用于各种粘弹性模型,可模拟任意震源及所产生的各种体波、面波,数值结果表明具有很高的计算精度。     

弹性波散射与多重散射的T矩阵方法

研究弹性波散射与多重散射的T矩阵方法。首先,基于Helmholtz体内和体外公式推导了对应于圆柱型散射体的T矩阵元素的具体表达式;接着分析了在含多个随机分布圆柱型散射体的随机非均匀介质中弹性波的多重散射并给出在统计平均意义下的相干波的定义以及波速和衰减系数计算公式;最后,针对Ge/Al、Sic/Al复合材料用Matlab进行了编程和数值计算;计算单个柱型散射体的散射截面以及随机非均匀介质中相干波的速度和衰减系数,分析了这种介质的频散特性。     

垂直断层破碎带对Rayleigh波传播与场地地震动反应的影响

采用基于时域粘弹性人工边界的Rayleigh波输入方法,对具有不同宽度、不同剪切波速的垂直断层破碎带的场地进行了Rayleigh波作用下的动力时程反应分析,研究了垂直断层的宽度与剪切波速变化对Rayleigh波传播及场地地震动反应的影响。计算结果表明:在Rayleigh波入射一侧的中远场地表面,断层宽度与剪切波速的变化对地面峰值位移反应影响很小;在Rayleigh波入射一侧的断层角点以及距角点较近处,随着断层宽度的增大或剪切波速的减小,峰值位移反应均增大;在另一侧中远场地面的峰值位移反应随着断层宽度的增大或剪切波速的减小而减小,且水平位移峰值减小的程度略大于竖向位移峰值。由此可以看到,软弱的断层破碎带对通过的Rayleigh波有削弱作用,且破碎带越宽或剪切波速越小削弱作用越明显。     

横观各向同性液体饱和多孔介质中平面波的传播

基于孔隙介质的Ｂｉｏｔ理论１，研究了横观各向同性液体饱和多孔介质中平面波的传播特性。首先导出了波传播的特征方程并给出了其解析解，结果显示：有４种不同波速的平面体波传播；第一准纵波，第二准纵波，准横波和反平面横波。文中给出了波速和衰减的解析表达式，数值计算了频散曲线和衰减曲线，并讨论了各类准体波位移之间的耦合关系。     

利用瑞雷波法评价既有铁路路基

根据瑞雷波频散相速度能够反演地层的剪切波速度结构 ,瑞雷波的特性使得瑞雷波法具有良好的场地适应性 ,在望京至和平里的既有铁路路基探测中 ,结合反射波法的CDP技术和瑞雷波谱分析技术 ,获得了 2 -D剪切波速度剖面 ,根据剪切波速度与轻型动力触探击数的一致性 ,从而能够在垂直和水平两个方向上对铁路路基进行评价     

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

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

流固耦合地震波动问题的显式谱元模拟方法

流固耦合地震波动问题主要研究由流体和固体构成的复杂系统中地震波传播特性及其规律. 传统模拟方法中一般以声波方程、弹性波方程的数值解分别描述理想流体和弹性固体中的波动, 并实时地处理两种不同性质介质之间的相互耦合作用, 数值格式复杂且限制数值模拟精度与计算效率. 本文采用谱元法结合多次透射公式人工边界条件实现了一种流固耦合地震波动问题的高阶显式数值计算方法. 该方法利用了流固耦合问题统一计算框架,可将饱和多孔介质的Biot波动方程分别退化为理想流体的声波方程和弹性固体的弹性波方程. 通过P波垂直入射的水平成层理想流体-饱和多孔介质-弹性固体场地模型、P波斜入射的不规则层状界面以及任意形状界面的理想流体-饱和多孔介质-弹性固体场地模型等三个算例, 与传递函数法解析解以及集中质量有限元法计算结果进行对比分析, 证明了本文方法的正确性与有效性. 数值模拟结果表明, 本文方法相较传统有限元法可以少得多的节点数量获得更高的数值精度, 并且在较宽的频率范围内都能可靠地模拟出流固耦合系统的动力响应, 充分体现出本文方法兼顾高精度、计算效率和复杂场地建模灵活的特点.      

Wave propagation in a piezoelectric rod of 6 mm symmetry

The wave propagation in a piezoelectric rod of 6 mm symmetry is investigated by applying a 3-D piezoelectric elastic model. A self-adjoint method is introduced to solve this problem, this method avoids calculating the generalized eigenvalue equation, it completely draws the dispersion curves in the forms of Quasi-P wave, Quasi-SV wave and Quasi-SH wave under the self-adjoint boundary condition, and it can evaluate the dispersion curves of all kinds of boundary conditions. As an example, the dispersion curves of PLT-5H are completely drawn, we also found the Quasi-SV wave has standing wave phenomenon in the PLT-5H rod. In addition the relation of dispersion curves among different boundary conditions is discussed, and an experiment method is introduce to decide the dispersion curves for another boundary conditions.     

GUIDED CIRCUMFERENTIAL WAVES IN DOUBLE-WALLED CARBON NANOTUBES

A model of guided circumferential waves propagating in double-walled carbon nan- otubes is built by the theory of wave propagation in continuum mechanics,while the van der Waals force between the inner and outer nanotube has been taken into account in the model.The dispersion curves of the guided circumferential wave propagation are studied,and some dispersion characteristics are illustrated by comparing with those of single-walled carbon nanotubes.It is found that in double-walled carbon nanotubes,the guided circumferential waves will propagate in more dispersive ways.More interactions between neighboring wave modes may take place.In particular,it has been found that a couple of wave modes may disappear at a certain frequency and that,while a couple of wave modes disappear,another new couple of wave modes are excited at the same wave number.     

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.     

NUMERICAL SIMULATION OF ELASTIC WAVE PROPAGATION IN ANISOTROPIC MEDIA

A new numerical simulation algorithm is presented for the elastic wave propagation in heterogeneous anisotropic media. We make discretization of the computational domain by using triangular and quadrangular grids. The scheme is based on integral equilibrium at each node to simulate the elastic wave propagation in heterogeneous anisotropic media. The scheme is very flexible in dealing with arbitrary surface topography, inner openings, liquid-solid boundaries and irregular interfaces. Moreover, the free-surface condition of complex geometrical boundaries can be satisfied naturally. This work is an extension of the grid method for the elastic wave propagation in heterogeneous isotropic media, and a quadrangular grid with low computational cost is also introduced. Project supported by the National Natural Science Foundation of China(19672016).     

各向异性介质中弹性波的数值模拟

提出了一种非均匀各向异性介质中弹性波传播的数值模拟算法。该方法可以灵活地运用于具有任意地表形状、内部孔洞、固液边界和不规则内部交界面的介质情况,另外,该方法自然满足复杂几何边界的自由表面条件。这种基于三角形和四边形离散网格的算法使用的是围绕每个节点的积分平衡方程,而不是其它有限差分法中使用的各个节点满足的弹性动力学的微分方程。该文工作是非均匀各向同性介质中弹性波传播格子法研究的继续。除了研究各向异性介质中波的传播以外,还给出了一种能够省时的四边形网格的格子法。     

Prediction on dispersion in elastoplastic unsaturated granular media

This study focuses on the propagation of the plane wave in the elastoplastic unsaturated granular media, and the wave equations and dispersion equations are derived for the media under the framework of Cosserat theory. Due to symmetry, five different wave modes are considered and predicted for the elastoplastic unsaturated granular media based on the Cosserat theory, including two longitudinal waves, one rotational longitudinal wave and two coupled transverse–rotational transverse waves. The correspondence is discussed between these Cosserat wave modes and the classical wave modes. Based on the dispersion equations, the dispersion behaviors are obtained for the five Cosserat wave modes. The results indicated that the different stress-strain stages,including the elastic, hardening and softening stages, have obvious effect on the dispersion behaviors of the Cosserat wave modes.     

Surface wave transference in a piezoelectric cylinder coated with reinforced material

The present study deals with the propagation of a polarized shear horizontal(SH)wave in a pre-stressed piezoelectric cylinder circumscribed by a self-reinforced cylinder.The interface of the two media is assumed mechanically imperfect.For obtaining the dispersion relation,the mathematical formulation has been developed and solved by an analytical treatment.The effects of various parameters,i.e.,the thickness ratio,the imperfect interface,the initial stress,the reinforcement,and the piezoelectric and dielectric constants,on the dispersion curve are observed prominently.The dispersion curves for different modes have been also plotted.The consequences of the study may be used for achieving optimum efficiency of acoustic wave devices.     

WAVE PROPAGATION IN A TRANSVERSELY ISOTROPIC THERMOELASTIC SOLID CYLINDER OF ARBITRARY CROSS-SECTION

The wave propagation in an infinite, homogeneous, transversely isotropic solid cylinder of arbitrary cross-section is studied using Fourier expansion collocation method, within the frame work of linearized, three-dimensional theory of thermoelasticity. Three displacement potential functions are introduced, to uncouple the equations of motion and the heat conduction. The frequency equations are obtained for longitudinal and flexural (symmetric and antisymmetric) modes of vibration and are studied numerically for elliptic and parabolic cross-sectional zinc cylinders. The computed non-dimensional wave numbers are presented in the form of dispersion curves.     

一维声子晶体带隙的非局部辛分析

周期性弹性复合结构(声子晶体)中传播的弹性波存在特殊的色散关系:弹性波只能在某段频率范围内无损耗的传播,该频率范围称为通带.一维声子晶体的色散问题可以看作分层介质中弹性波的传播问题,利用二维弹性理论予以分析.为了研究非局部效应对声子晶体带隙特性的影响,将Eringen的二维非局部弹性理论引入到Hamilton体系下,利用精细积分与扩展的Wittrick Williams算法可获取任意频率范围内的本征解.通过对不同算例的数值计算,分析和对比了非局部理论方法与传统局部理论方法的差别.并进一步指出了该套算法的适用性和优势所在.     

SH-type waves dispersion in an isotropic medium sandwiched between an initially stressed orthotropic and heterogeneous semi-infinite media

The present paper studies the propagation of shear waves (SH-type waves) in an homogeneous isotropic medium sandwiched between two semi infinite media. The upper half-space is considered as orthotropic medium under initial stress and lower half-space considered as heterogeneous medium. We have obtained the dispersion equation of phase velocity for SH-type waves. The propagation of SH-type waves are influenced by inhomogeneity parameters and initial stress parameter. The velocity of SH-type wave has been computed for different cases. We have also obtained the dispersion equation of phase velocity in homogeneous media in the absence of initial stress. The velocities of SH-type waves are calculated numerically as a function of kH (non-dimensional wave number) and presented in a number of graphs. To study the effect of inhomogeneity parameters and initial stress parameter we have plotted the velocity of SH-type wave in several figure. We have observed that the velocity of wave increases with the increase inhomogeneity parameters. We found that in both homogeneous and inhomogeneous media the velocity of SH-type wave increases with the increase of initial stress parameter. The results may be useful for the study of seismic waves propagation during any earthquake and artificial explosions.