动力方程求解的显式积分格式对数值模拟波动特性的影响分析

基于李小军等提出的显式积分格式应用于波动有限元模拟时的递推形式,重新定义了波动传递函数;以一维均匀离散体系为模型,对给出的传递函数进行了理论推导和求解,并通过数值试验验证了其正确性;通过对传递函数的分析研究,探讨了这种显式积分格式对离散网格中波动传播规律包括截止频率、频散现象和能量耗散等方面的影响．以期为该显式积分格式在波动问题求解中的应用提供更详细的理论参考．      

三维直流电场数值模拟的拟解析近似法:张量近似

拟解析近似方法是一种解决电磁场散射问题的快速求解积分方程的近似方法,它绕开了传统数值方法中的求解大型代数方程组或大型矩阵问题,适用于强散射和大扰动问题.本文应用孙建国提出的求解异常电场积分方程的张量拟解析近似理论公式,研究用其求解直流电场积分方程.利用接近实际的地电模型对异常电场进行模拟,研究了均匀场中异常球体的张量拟解析近似解;并对均匀场中的立方体异常体进行了数值计算.效果良好并具有很高的计算速度.研究结果为三维直流电场快速正反演模拟打下基础.     

基于递归积分的复频移PML弹性波方程交错网格高阶差分法

常规完全匹配吸收边界(PML)对以近掠射角入射到界面上的波以及低频波、损耗波都会产生虚假边界反射.基于递归积分的不分裂复频移PML算法,利用复频移拉伸函数,极大地改善了PML边界条件的性能,我们进一步推导出基于递归积分的不分裂复频移PML弹性波方程交错网格高阶差分法,对长条形介质模型进行数值模拟,与常规PML算法进行比较说明该算法对以掠射角入射到PML界面的波以及PML层内损耗波的吸收效果.     

基于耦合反射/透射系数单程波传播算子的地震波模拟研究

单程波近似实际上是一种多次前向散射和单次后向散射近似.利用单程波近似来描述波传播可以极大地节省地震数值模拟的计算时间和内存,实现地震波长距离传播模拟和三维地震模拟快速计算.本文基于单程波近似和波动积分方程的分离变量逼近,从广义Lippmann-Schwinger波动积分方程推导出耦合反射/透射系数的单程波传播算子.该算子由两部分构成：分离变量Fourier单程波传播算子和薄板间的反射/透射系数表达.前者将常规的Fourier分裂步单程波传播算子(SSF)推广适应横向强速度变化介质和大角度传播波场.后者是利用垂直波数来表示反射/透射系数,自然耦合到波场传播的计算过程中,其为地质界面倾角的隐式表达,精确描述振幅随入射角的变化,能适应任意复杂的模型.通过两个数值算例和一个实际地质模型的计算,本文将该方法和边界元法进行了比较,结果表明：在算例给出的介质横向速度变化情况下,本文提出的方法在相位和振幅方面与全波数值方法基本吻合.     

泥石流的二维数学模型

泥石流是在重力作用下,由砂粒石块和水等组成的固液混合物,是一种发生于山区的复杂的地质灾害现象。泥石流主要是由暴雨诱发引起的,它沿着复杂的三维地形高速流动,具有流体流动的特性。为了模拟泥石流的运动规律,预测降雨诱发的泥石流的到达距离和泛滥范围,减少和避免泥石流引起的灾害,把泥石和雨水组成的固液混合物假定为遵循均匀、连续、不可压缩的、非定常的牛顿流体运动规律。基于质量守恒方程和Naiver-stokes方程,采用深度积分方法,推导出了一个模拟泥石流运动的二维数学模型。所有方程式可用有限差分法来求解。结合GIS,该模型可用于预测泥石流的流动距离和泛滥范围,以及泛滥范围内的危险房屋和路段,也可以用于泥石流灾害的风险性分析。     

边坡稳定性分析的无条分法

基于滑面应力修正技术并取整个滑体为受力体,实现了对滑体的整体极限平衡分析。同时,利用Green公式将有关域积分转化为边界积分,而无需再对滑体进行条分。该方法数值特性好,亦可方便地将二维的成果向三维推广,从而实现满足全部6个平衡条件的三维严格极限平衡法。本文是二维部分的成果。     

竖向荷载下桩基础弹性分析的改进计算方法

将Randolph剪切位移方法中桩身位移与桩端位移的函数关系简化为一多项式,并与Poulos积分方程中土体柔度系数矩阵相结合,提出了一种竖向受荷单桩弹性分析的改进计算方法,从而避免了Poulos积分方程法中的差分运算以及由此带来的其他矩阵运算,同时比Randolph方法更能准确模拟桩身剪切应力的分布,并将单桩的改进计算方法应用于群桩分析。将改进的计算方法与Poulos、Randolph方法以及Chow混合方法的单桩和群桩计算结果进行了比较,结果表明：改进计算方法是可行的,可满足精度要求。     

边界积分方程用于电阻率Zohdy反演的初步研究

研究利用边界积分方程进行电阻率Ｚｏｈｄｙ反演的有关技术,结果表明在探测区域划分为数千个单元的情况下,为节省正、反演过程中所需要的计算机内存和 ＣＰＵ时间,可采取以下措施：（１）假定电位分块线性变化,可以使节点数目显著减少,这里每一块都包含若干个单元;（２）利用边界积分方程的特点,可事先计算并存储所有必要的边界积分值,供每一步迭代中调用;（３）利用边界积分方程所特有的数值延拓功能,可以进一步减少节点数目．为了说明这些手段的有效性,先给出了一个数值模拟成像结果,它清楚地表明基于边界积分方程的电阻率 Ｚｏｈｄｙ反演是快速而有效的,然后给出一个实测剖面的成像实例,成像结果经钻孔检验正确．     

移动荷载下黏弹性半空间体上双层板的动力响应

采用黏弹性半空间体上无限大双层板模型模拟路面结构,通过三重Fourier积分变换得到单位脉冲荷载作用下路面动力响应的Green函数。根据线性系统的叠加原理,利用广义Duhamel积分推导出移动荷载作用下路面动力响应的解析解。采用自适应Simpson法编制了计算奇异、振荡函数的广义积分计算程序,完成了路面动力响应从波数-频率域到时间-空间域的转换。结合算例,对移动荷载作用下路面的振动规律进行研究,讨论上、下层板厚度和板材料的弹性模量对路面动力响应的影响,明确了路面结构的振动特性,研究结果可为路面结构的设计、施工提供理论指导。     

求解反应扩散方程的紧致隐积分因子方法

积分因子(IF)方法是近年来提出的求解刚性常微分方程组的一种有效的数值方法。本文应用改进的紧致隐积分因子(cIIF)方法求解二维反应扩散方程。在空间离散上采用二阶中心差分方法。数值模拟得到各种斑图结构与所引文献结果相当一致。