二维介质参数的大扰动反演方法

对非均匀介质参数反演问题进行了研究，并提出了用于反演二维介质参数的广义射线近似方法．利用参考场量和扰动变量对声波方程中的介质参数进行处理，并利用Ｇｒｅｅｎ函数理论得到扰动参数比的积分方程．基于非均匀介质中波函数的局部理论和射线理论，引入了全波场的广义射线近似形式，通过定义介质参数函数，把反演目标归结为其第一类Ｆｒｅｄｈｏｌｍ积分方程．利用积分变换方法得到二维介质的介质参数函数，从而得到介质参数，在Ｂｏｒｎ近似方法中，反演的介质参数扰动不能超过２０％，但是在本文中介绍的方法能够有效地反演其扰动比不超过５０％的变化情况     

介质参数反演的广义射线近似方法

在对无粘性介质参数反演问题进行的研究中，引入一种全波场广义射线近似形式，提出一种新的反演参数的方法，文中，首先对由弹性波动方程演变成的声波方程进行分析，引入背景场量和扰动量，并结合Ｇｒｅｅｎ函数理论，得到了介质参数的积分方程；然后结合前人对非均匀介质中波函数局部理论的定性分析，引入一种全波场广度射线近似形式，把问题归结为一个第一类Ｆｒｅｄｈｏｌｍ积分方程；最后对半空间问题层状介质模型进行了反演，算     

三维非均匀介质中弹性波传播的数值模拟

提出了一种三维非均匀介质中弹性波传播数值模拟的方法,文中称为三维格子法。该算法是二维格子法（一种二维非均匀介质中P－SV波传播的数值模拟算法）向三维非均匀介质情况的推广。在空间离散上该文方法与有限元方法类似,容许根据连续体的形状和介质分界面任意剖面网格,且自然满足自由表面边界条件。不同于常规有限差分法在各个节点上满足动力学微分方程,该算法通过满足各节点周围格子的整体平衡（积分平衡方程）来对问题进行求解,三维格子法所需的计算机内存及计算耗时与同阶精度的规则网格有限差分法相当。算例表明,该文提出的三维格子法具有较高的精度且可很好地模拟三维复杂形状地表对弹性波的反射和绕射。     

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

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

被介质包围的非线性弹性杆波动方程的精确解

应用推广的tanh函数方法对被介质包围的非线性弹性杆波动方程进行了求解．得到了该方程的精确孤波解,对部分结果通过数值模拟得到了解的图像．本文获得的结果完善了相关文献已有的研究成果．从求解过程可以看出,推广的tanh函数方法是一种简便有效的方法,可用于求解其它非线性方程或方程组．     

非均匀损伤介质中波传播的数值解

对弹性波在非均匀损伤介质中的传播理论进行了研究。通过将非均匀损伤区域离散成分层均匀的区域，结合相邻区域交界面处的连续条件，推导出了以右行波、左行波为状态向量的波动方程和传递矩阵。对几种非均匀损伤介质中波的传播进行了实例数值计算，并和其解析解的结果进行了比较，讨论了弹性波在非均匀损伤介质中传播的一般性质。     

有限散射信号下二维缺陷形状识别的罚函数方法

研究在有限照射角度和频带宽度下二维缺陷的形状识别问题。首先,通过引进介质参数扰动函数,建立介质参数扰动函数和弹性波散射场之间的非线性关系,并将所关心的缺陷的形状识别问题转化为关于扰动函数的反演;然后,利用变分技术和优化方法求解,为了弥补散射数据的不足,在总的目标函数中,采用附加度量函数作为罚函数;最后,对后场散射远场测量时有限照射角度和频带宽度下几种典型缺陷进行了模拟识别,表明了;表明了罚函数法的有效性。     

周期性非均匀介质中气相爆轰波演变模式研究

气相爆轰波在周期性非均匀介质中的起爆, 稳态传播和失效机制都极为复杂, 很多物理机制尚不明确, 是当前爆轰物理领域研究的热点和难点. 本文使用反应欧拉方程和两步化学反应模型对爆轰波在非均匀介质中的传播机理进行了数值模拟研究, 非均匀性由横向周期性分布的温度扰动体现, 重点分析不同波长、不同幅度的温度扰动对波阵面波系结构的影响. 计算结果表明, ZND爆轰波在温度扰动下向胞格爆轰波的转变主要受制于两种竞争性因素: 一是爆轰波内在的不稳定性; 二是温度扰动的波长和幅度, 前者是内因, 后者是外因. 温度扰动的存在抑制横波的发展, 延迟了ZND爆轰波向胞格爆轰波的演化, 并且内在不稳定性的增加可以减慢这种延迟现象. 这说明, 温度扰动可以在一定的范围内抑制胞格不稳定性的发展, 但是不能够终止这一过程. 温度的不连续性使得爆轰波阵面更为扭曲, 并在横波附近存在较弱的三波点结构, 即温度扰动可增加爆轰波固有的不稳定性, 改变爆轰波阵面的传播机理. 幅值较大的人工温度扰动可抑制爆轰波的传播和爆轰波自身的不稳定性. 爆轰波阵面胞格结构的形成取决于温度扰动与其自身的不稳定性.      

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

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

各向异性展状介质中弹性波的传播矩阵解法

基于广义胡克定律及混和变量弹性波方程，解析求得各层介质内位移、应力传递矩阵，给出了直角坐标系下各向异性层状介质中弹性波的传播矩阵解法．该方法适用于非轴对称各向异性和点源作用，较好地解决了数值计算中有效数字精度损失问题．数值结果表明，计算效率、准确性及稳定性均较好．     

The parameter perturbation method on elastic wave equation in inhomogeneous medium

I.IntroductionTheelasticwaveininhomogeneousmediumiscomplicatedbecauseofthediffracting,scattering,andtransmutingofthewavetapes.Exceptforsomesimpleandregularmediummode1s,thesolutionofelasticwavehasnotbeengotyet.Nowadays,theresearchoftheelasticwavescattering…     

Elastodynamic analysis of inhomogeneous anisotropic bodies

The integral equation method is presented for elastodynamic problems of inhomogeneous anisotropic bodies. Since fundamental solutions are not available for general inhomogeneous anisotropic media, we employ the fundamental solution for homogeneous elastostatics. The terms induced by material inhomogeneity and inertia force are regarded as body forces in elastostatics, and evaluated in the form of volume integrals. The scattering problems of elastic waves by inhomogeneous anisotropic inclusions are investigated for some test cases. Numerical results show the significant effects of inhomogeneity and anisotropy of materials on wave propagations.     

An analytical study on the Rayleigh wave generation in a stratified structure

In this paper, an analytical method is used to investigate the Rayleigh wave generation in a stratified structure and the wave generation in a dry sandy layer constrained between the couple stress and inhomogeneous orthotropic half-spaces. This study is devoted to analyzing the impact of various effective parameters associated with the media on the phase velocities of the wave. The displacement components for each medium are derived by implementing the separable variable method. The frequency equation is secured by using the displacement components in the boundary conditions, imposed at the interfaces between the layer and half-spaces. Moreover, the secured equation is the relation between the phase velocity and the wave number. Numerical computations are performed, and graphical representations are demonstrated between the phase velocity and the wave number for both phase velocities with different values of the parameters.The comparison between the phase velocities is observed for the same value of each parameter.     

Love wave propagation in heterogeneous micropolar media

The present paper framed to study the impact of heterogeneity on propagation of Love wave in a heterogeneous micropolar layer over an elastic inhomogeneous stratum, when both rigidity and density are assumed to vary linearly with depth. The equations of motion have been formulated separately for layer and half-space under suitable boundary conditions. Analytical solution for the dispersion equation has been obtained using method of separation of variables by means of the Airy function and Whittaker function. Some particular cases have also been investigated. Further, as a special case the velocity equation for isotropic layer over a homogeneous half-space coincides with the standard result of Love wave. Numerical calculations of frequency relation have been performed and depicted by means of graphs to exhibit the substantial impact of heterogeneity, micropolar parameters and wave number on the phase velocity of Love wave. The wave velocity is strongly influenced by these parameters.     

Finite element-boundary integral formulation for electromagnetic scattering

A new formulation is described which combines the most robust attributes of the volume finite element and surface integral equation approaches to electromagnetic boundary value solutions. The result is a numerical technique which may be applied to scattering problems involving configurations having metallic surfaces and inhomogeneous penetrable material situated in open spatial regions. This is accomplished by way of coupling internal region finite element modal field solutions to equivalent currents on the surrounding boundary surface through an appropriate surface integral equation. The method is demonstrated for the special case of scattering by axisymmetric inhomogeneous penetrable objects. Example numerical calculations are presented for validation of the procedure and potential problem areas are discussed.     

GENERALIZED RAY APPROXIMATE METHOD TO INVERSE MEDIUM PARAMETERS

In this paper,the inverse problem of the medium parameters in an inhomogeneousmedium is studied and a generalized ray approximate form of the total wave field is described.First,the acoustic wave equation derived from the elastic wave equation is studied,the referential variablesand perturbational variables are introduced,and the integral equation of the medium perturbational pa-rameters is obtained.Then from the point of view of the local principles of the wave function in an in-homogeneous medium,a generalized ray approximate form of the total wave field in an inhomoge-neous medium is described,and attention is focused on the Fredholm integral equation of the firstkind.Finally,the medium parameters in half-plane are inversed.Numerical examples show when theperturbations of the medium parameters are about 0.5,this method can effectively inverse its varia-tion.Apparently,this method is better than the conventional Born weak scattering approximation.     

各向异性体对瞬态SH波散射问题的边界元方法

本文利用Fourier变换和加权残数法建立了正交各向异性体反平面瞬态波散射问题的边界积分方程,构造了在边界上能满足波动方程的边界元函数,它比常用的二次元和高次元具有较好的适应性。数值结果表明利用该形函数计算高频情形下的散射问题具有很高的精度。