一类广义积分的算法及实现

Laplace变换与B样条在社会生活的各个领域都有着广泛的应用，考虑将二者联系起来，运用B样条方法求解Laplace积分。对Laplace积分具有导数不连续点的像原函数，运用重节点B样条进行拟合，再引入差商，以两个不同出发点导出Laplace积分的算法。并在均匀网格分划下，给出一系列解析公式，然后用算例编程计算，实例表明简单、快速而且便于计算机实现。     

实际井眼条件下过套管电阻率测井响应的传输线方程正演算法

过套管电阻率测井通过测量套管壁电势实现测量地层视电阻率,不同的地层模型对过套管电阻率测井会有不同的测井响应. 在传输线方程系数中本文考虑了径向含多个界面地层的影响,对轴向（井轴方向）呈层状、径向阶跃变化地层模型,给出了电势分布递推公式,计算了过套管电阻率测井响应,绘出了测井响应曲线. 算例表明,该方法有较强的轴向地层界面分辨能力;对径向均匀地层计算结果非常接近地层模型电导率,能较好地反映实际地层特征;对高电导地层,低电导水泥环有明显的测井响应;当径向存在多个地层界面时计算结果则为径向各层电导率的综合反映,不是某一径向层的电导率;算例显示该方法有较快（算例运行时间在1s内）的计算速度.     

动力方程求解的显式积分格式及其稳定性与适用性

文献（１）给出了一种求解有阻尼体系动力方程的显式积分格式,文中以数值表格的形式给出了格式的稳定性条件,本文对该格式的稳定性问题作了进一步的分析,并给出了其计算稳定性条件的表示式。本文还着重讨论了基于这一显式积分格式的推导过程而派生出的另一形式的积分格式的稳定性,并指出了该派生格式的适用性问题。     

计算断层自发破裂过程的离散网格划分研究

离散化网格的空间步长选取在各种数值算法中都是一个很受关注的问题,在全空间均匀介质模型和简单离散化方案下,利用边界积分方程方法研究在自发破裂求解过程中动力学参数组合Dc和Te (Dc为临界滑动弱化位移,Te为有效的断层破裂面上的初始应力)对计算网格划分的影响,初步得到了Dc和Te参数空间中的有效计算网格的选取规律,对合理、有效地运用边界积分方程方法计算地震震源的破裂过程具有重要的指导意义.      

沿圆曲线的Radon变换数值解

研究具有紧支集且在支集内连续的二元函数沿上半圆曲线的Radon变换反演问题。基于对投影函数的Fourier变换，反演问题可以归结为具有弱奇性及震荡核的Abel积分方程的求解。我们证明了当圆曲线中心及半径在一定范围内变化时，在已知沿上半圆曲线的Radon变换情况下，这个积分方程的解具有唯一性，并给出了消除Abel积分方程弱奇性的数值方法。在考虑投影数据噪声的情况下，给出了多次加权改善系数矩阵条件数稳定的数值方法，并通过数值模拟验证所提出方法的有效性。      

时域谱元法的质量特性模型及其构建方法

研究了构建时域谱单元质量特性模型的数学机制,针对时域切比雪夫谱单元和勒让德谱单元建立了一种直接导出谱单元一致质量矩阵和集中质量矩阵的统一数学方法,对比分析两种谱单元质量特性模型的特征,并从物理角度探讨了谱单元质量特性模型的合理性.研究表明,数值积分点与谱单元节点选取是否一致是决定时域谱单元形成一致质量模型或集中质量模型...     

瞬态弹性波散射的变换域积分方程法

针对瞬态弹性波散射的问题,从弹性动力学问题的积分表示定理出发,采用Laplace变换的方法,得到了变换域内均质体位移场的积分方程表示;在此基础上推导了适合瞬态弹性波对异质体散射求解的变换域位移场积分方程。     

蒸发条件下土壤水分响应模拟

基于有限元(FEM)和改进的积分型Richards方程解法(IRE方法)对蒸发条件下5种土体土壤水分响应进行了研究.数值实验结果表明:在土壤表面潜在蒸发量0.50 cm/d的情况下,5种土体土壤含水率变化曲线均呈现单拐点两阶段的特点,拐点出现在地表下20cm左右,拐点上部区域曲线曲率大于下部区域,两阶段的划分以15d左右为界,前阶段比后阶段的土壤水分变化快;蒸发模拟结果很好的证明了蒸发三阶段理论.总蒸发量和下边界排水量与土壤结构密切相关,而总水量变化量和变化率与土壤质地有关.IRE方法与FEM模拟结果基本一致,解法相对简单,模拟结果可靠性高.     

多层弹性平面问题解的精确刚度矩阵法

从平面弹性力学的基本方程出发,利用Fourier积分变换等数学手段,推导出了单层平面问题的刚度矩阵,然后按有限元法组成总体刚度矩阵。通过求解由总体刚度矩阵所构成的代数方程和Fourier积分逆变换,得到在任意静荷载作用下多层弹性平面问题的精确解。由于刚度矩阵不含有正指数项,计算时不会出现溢出现象,从而克服了传递矩阵法的缺点。推导过程中摒弃了应力函数的选择,使得问题的求解更加合理化,同时也为进一步研究此类问题如温度场、动力学等奠定了理论基础。计算实例证明了推导结果的准确性。     

Lippmann-Schwinger积分方程广义超松弛迭代解法及其收敛特性

为克服Born级数只适用于弱散射和短程传播的弱点,从光学散射理论中引入了针对强光学散射问题的Lippmann-Schwinger方程广义超松弛迭代解法,并对这种方法在反射地震条件下的算法表现、收敛特点、计算效率以及计算精度等问题进行了分析和讨论.理论分析和数值计算结果均表明:(1)与散射级数重整化相比,超松弛法具有更坚实的数学基础;(2)超松弛法可以有效地克服Born级数的弱点,得到与有限差分相当的数值模拟结果;(3)在反射地震条件下LS方程超松弛法表现优异,完全可以解决反射地震数值模拟中的强散射问题.