第二类非线性Fredholm型积分方程数值解

配置法研究了地球物理中常见的第二类非线性 Fredholm 型积分方程的数值解法,将第二类非线性 Fredholm 型积分方程转化为非线性代数方程组进行求解,采用高斯数值积分公式,给出了数值计算的具体实例.利用Matlab软件的符号运算功能编程计算,克服了非线性方程难于变成求解的困难,数值例子表明该方法编程简便有效.对非线性积分方程和非线性代数方程组的求解都有重要价值.     

越流含水层井群系统地下水流的边界元数值模拟

本文提出一种求解越流含水层井群系统的边界元方法。该法把抽水井井壁作为内边界,通过积分处理,不需要对井壁内边界进行剖分,通过推导求得水头H及其法向导数的边界积分方程,离散化后得到求解该积分方程的一组线性代数方程组、该法用于实例计算表明,计算精度较高。     

同伦摄动法在一类线性积微分方程初值问题中的应用

积微分方程定解问题在数学与其他科学领域里有着重要的应用,利用积分,将一类积微分方程定解问题转为与之等价的第二类Fredholm-Volterra积分方程,然后利用同伦摄动方法求解第二类Fredholm-Volterra积分方程,可得积微分方程定解问题的解,最后利用matlab符号计算功能对实例进行计算,验证了同伦摄动法在求解积微分方程定解问题中是有效的.     

用运界元法求积分方程数值

本给出了用边界元法积分方程的算法。通过处例验算表明，这种算法综合了边界元法精度高和用矩最法解积分方程速度快的优点，是一种求解积分方程十分有效的算法。     

用边界元法求积分方程数值

本文给出了用边界元法求积分方程的算法。通过算例验算表明，这种算法综合了边界元法精度高和用矩量法解积分方程速度快的优点，是一种求解积分方程十分有效的算法。     

复合单元法在求解有自由面渗流问题中的应用

论文指出了在有限元分析过程中地下水自由面上、下的区别,提出可采用复合单元法处理常应变交截单元,并以四面体单元为例,推导出了自由面与单元交截的截面方程,详细说明了复合单元法的原理、系数矩阵的求解方法、推导出体积分和曲面积分的求解公式。     

挠曲线复位的微分方程解法求梁的位移

梁位移求解的方法主要有两种：积分法和叠加法。积分法的困难是当多个荷载同时作用时，要按控制面分段列弯矩方程，这样，两次积分时带来的积分常数较多，导致依边界条件、连续条件定这些常数时要解多个联立方程，计算繁杂。梁位移求解的叠加法是借用典型荷载下梁已知积分结果来进行对应位置位移的叠加，具体是查转角方程、挠曲线方程表格。从表中知，梁任一位置处的位移方程大多是多次多项式，很难记忆，故离开表格便无法采用。本文找到了一种新的方法求梁的位移，通过设置比拟梁，使比拟梁的挠曲线复位于原梁，从而建立了相应的微分方程，求解此方程后便可得梁的位移，避开了分段多次积分或查表的弱点。     

二维位场向上延拓和向下延拓的样条函数法

本文用样条函数法实现了二维位场的向上延拓和向下延拓,向上延拓用样条函数法计算泊松积分,向下延拓用样条函数法求解第一类Ｆｒｅｄｈｏｌｍ积分方程。模型和实例计算表明,本文方法正确可行且精度高。     

起伏地表条件下的声波散射数值模拟的积分方程法

从散射理论的角度来看,起伏地表可以看作是一种特殊的扰动介质,因此应用散射积分方程求解起伏地表条件下的散射场在理论上是可行的。从三维频率域声波方程出发,由格林函数定理,得出起伏地表条件下的散射积分方程。散射积分方程为关于起伏地表的面积分和与速度扰动体有关的体积分之和,同时给出了格林函数在奇异点的积分方法。由于数值离散求解积分方程存在着计算时间太长和存储内存不足的问题,采用电磁散射积分方程的拟解析近似的方法。在假设反射函数为缓变函数的基础上,最终得到其近似表达式,因此散射场的数值求解不必再借助于代数方程组,只要进行数值积分即可。这种方法避开了传统数值计算方法存在的问题,为地震散射波场快速正演模拟打下了基础。理论分析表明,这种方法适用于小扰动的问题。当扰动较大时,拟解析近似会产生较大的误差。     

二维位场向上延拓与向下延拓的样条函数法

本文用样条函数法实现了二维位场的向上延拓和向下延拓，向上延拓用样条函数法计算泊松积分，向下延拓用样条函数法求解第一类Ｆｒｅｄｈｏｌｍｔｍ积分方程。模型和实例计算表明，本文方法正确可行且精度高。     

三维直流电场标量拟解析近似法数值模拟精度分析与评价

拟解析近似方法是一种可以处理强散射或大扰动问题的求解积分方程的近似方法,该方法可以避免用传统数值方法解决问题时,所遇到的大型矩阵或大型代数方程组问题。根据孙建国在文章[14]中的求解异常电场积分方程的标量拟解析近似理论公式,对标量拟解析近似进行研究。对均匀场中的异常球体标量拟解析近似解进行计算,并与理论解进行了对比,验证了其精度,证明拟解析近似方法求解直流电场积分方程是可行的。通过对标量拟解析近似方法在直流电场数值模拟的研究,为三维直流电场快速正反演模拟打下了基础。     

国外瞬变电磁法中正反演问题研究现状

国外在瞬变电磁法正反演问题研究中,运用有限差分法、有限单元法及积分方程法直接在时间域中求解以及运用富氏变换法、衰减谱法和数字滤波法从频率域变换等数值计算方法,并对其在瞬变电磁资料解释中的效果作了评述.     

二重网格方法用于曲面位场转换

本文介绍了二重网格方法及在求解曲面位场转换的积分方程中的应用。给出该积分方程的二重网格方法的具体算法和二重网格迭代算子的结构图。并提出了三点注意事项。     

积分方程解广义逆函数值ε-算法的一种新用法

广义逆函数值Pad逼近(GIPA)被引入来求解积分方程,而广义逆函数值ε 算法是计算GIPA的一种有效递推算法。本文用函数值ε 算法求出GIPA的分母多项式的零点,然后取它的实部作为积分方程特征值的近似值,与其它算法进行了实例比较,显示出估计积分方程特征值的ε 算法的优越性。     

复合地基相互作用系数解法及参数分析

采用积分方程法推导出求解桩-桩、桩-土、土-土之间的相互作用系数所需要的第2类Fredholm积分方程,利用叠加原理对复合地基进行计算分析。通过与试验结果的比较,验证了计算方法的正确性,并对复合地基的沉降、各桩以及土之间的荷载分担情况进行了参数分析。该方法能够用来分析大规模复合地基,具有一定的工程应用前景。     

大规模群桩基础相互作用近似解耦方法

针对传统弹性理论法过高地估计了桩-桩相互作用效应,基于杆件有限单元法建立了半无限土体中群桩基础的桩侧剪应力求解方程,通过简化桩-桩的相互作用效应,将群桩基础桩侧摩阻力的求解方程近似解耦,实现了群桩基础桩身剪应力和位移的快速求解目标。通过两桩相互作用系数以及柔性承台下群桩基础差异沉降的对比分析,验证了该简化方法的合理性。参数分析结果表明,该简化方法计算得到的相互作用系数与严格的边界积分方程法解答较为接近,稍小于Poulos弹性理论法的计算结果;柔性承台下群桩的差异沉降在桩间距较小时与经典解答较为接近,而在桩间距较大时则存在一定的差别。该简化方法大幅减少了群桩计算工作量,适用于大规模群桩基础的快速计算要求。     

时变场中导电导磁体的边界元素法

该文研究了低频人工时变场（似稳场）中任意导电导磁矿体的边界元素法。通过对矢量位的边界积分方程，再以样条边界元素法进行离散，得到实分量的常系数方程组和虚分量的二阶微分方程组。此方程可利用时间积分的样条函数方法求解，解出适量位后，利用磁场强度与矢量位Ａ的相应关系而求得二次磁场的水平分量和垂直分量。     

瞬变弹性动力问题的一个边界元法

本文给出了瞬变弹性动力问题的一个边界元法。该法是利用威尔逊——θ法的差分公式,把运动方程化为椭圆型微分方程。根据贝蒂定理和动力点荷载的特解,可获得动力问题的边界积分方程。这个解法是在真实时间域内逐步求解的,不需要使用拉氏变换。一个应力波传播的数值算例证实了该方法使用方便,且解答精度较高。     

边界元法计算均匀谐变场中导电导磁体的电磁响应

本文讨论了在均匀谐变场激发下导电导磁体电磁响应的边界单元法计算问题。文中首先建立矢量位的边值问题,然后导出其边界积分方程,最后用边界元法进行求解。本方法的计算结果具有精度高、速度快的优点。     

多重网格准线性近似技术在三维航空电磁正演模拟中的应用

系数矩阵存储和线性方程组求解是限制三维电磁积分方程方法发展的主要因素。Zhdanov提出准线性（QL）近似技术,建立了复杂散射场与背景场的线性关系,有效地避免了积分方程中大型线性方程组的求解,但是该算法用于多源问题航空电磁正演模拟时精度不高。因此,本文提出一种基于多重网格准线性（MGQL）近似的算法,并利用系数矩阵的Toeplitz性质存储矩阵和快速傅里叶变换,实现了矩阵与向量的快速乘积、降低了计算复杂度,采用多重网格结合了积分方程方法和准线性近似解法的优点,在保证精度的条件下提高计算速度、减少存储量。针对不同类型网格的模拟实验表明,相比于传统积分方程方法,本文算法在保证计算精度的同时,可以将计算速度极大地提高（>10倍）。