共查询到20条相似文献，搜索用时 46 毫秒
1.
本文给出了用边界元法求积分方程的算法。通过算例验算表明，这种算法综合了边界元法精度高和用矩量法解积分方程速度快的优点，是一种求解积分方程十分有效的算法。 相似文献
2.
本文提出用二层介质格林函数作为余量加权函数,归化边界积分方程的边界单元法。通过对有水平覆盖层的点源二维、点源三维正演计算表明,本算法明显地优化了以往的边界元法,计算精度高,方程组规模小,节省内存空间,仅在一般微机上即可实现,扩大了边界单元法在电法勘探中的应用范围。 相似文献
3.
本文以二层介质中存在球体的模型为例，讨论了用边界元方法解决层状介质中存在三维不均匀体情况下点源场视电阻率的计算问题。文中通过对基本解的选取，避免了对水平界面的积分计算，并且在积分方程中消去了位的法向导数，使得积分方程的形式更为简单，也大大地减少了所需的边界剖分节点数和未知量的个数，节省了计算机内存，提高了计算精度。计算的结果表明，边界元法对本问题的研究是成功的。 相似文献
4.
本文讨论了在均匀谐变场激发下导电导磁体电磁响应的边界单元法计算问题。文中首先建立矢量位的边值问题,然后导出其边界积分方程,最后用边界元法进行求解。本方法的计算结果具有精度高、速度快的优点。 相似文献
5.
结合边界元法和离散元法,提出一种可以进行计算颗粒内部应力和破碎路径的方法。该方法利用离散元法求解颗粒的相互作用和每个颗粒上的荷载。然后利用边界元法计算颗粒的应力分布,为了实现动态平衡,将颗粒的加速度视为恒定大小的体力。但体力导致边界积分方程中出现域积分,故采用直线积分法将域积分转化为边界积分,以保证边界元法降维的优势。为了提高边界元的计算效率,对于几何形状相似的颗粒,以其中一个颗粒作为模板颗粒,只需要计算模板颗粒在局部坐标系中的系数矩阵,其他相似颗粒可以通过局部和全局坐标系之间的映射获得。在得到应力后,基于Hoek-Brown准则来判断颗粒是否破碎。此外,将破坏路径简化为直线,并采用最小二乘法拟合得到破坏路径。 相似文献
6.
7.
在边界积分方程离散的基础上,根据远场位移的特性,用衰减函数描述远场位移,阐述了自适应高斯积分算法,并结合前人成果—用解析法求出在远场上的强奇异积分,编制了适合于半无限域线弹性问题的计算机程序,并进行了事例验证,验证表明,该法比传统边界元具有更高的精度。 相似文献
8.
9.
频率测深的地形影响及其校正方法 总被引：3，自引：1，他引：3
从含源麦克斯韦方程出发，按照二维地形条件下的线源非齐次赫姆霍兹方程确定其边界条件，用边界元法计算了起伏地形条件下线源的频测影响曲线，并将其转换为二维地形点源赤道偶极装置频测曲线。总结了地形对频测曲线影响的规律，提出了地形校正的方法。 相似文献
10.
本文给出了瞬变弹性动力问题的一个边界元法。该法是利用威尔逊——θ法的差分公式,把运动方程化为椭圆型微分方程。根据贝蒂定理和动力点荷载的特解,可获得动力问题的边界积分方程。这个解法是在真实时间域内逐步求解的,不需要使用拉氏变换。一个应力波传播的数值算例证实了该方法使用方便,且解答精度较高。 相似文献
11.
本文用边界元法解决非稳定系统中水质预测的计算问题,建立了三维水质预测问题的边界元法计算模型,为简化计算手续、提高计算精度,讨论了Laplace变换边界元法,推广了Laplace变换边界元法中的Schapery方法。 相似文献
12.
一种新的二维电阻率模拟方法 总被引：2，自引：0，他引：2
目前，对地下介质作电阻率成像需要将研究区域网格化。本文以格林函数为数学工具导出了网格化地电模型二维问题的边界积分方程，该方程适合于电阻率成像。用边界单元法求解该方程可以得出空间电位分布，从而进行视电阻率模拟。计算与实验结果的对比表明，这种方法是正确、可行的。 相似文献
13.
边界元法在数值解法中利用拉普拉斯方程的基本解和格林公式,简单易行,可以解决物性参数分区均匀的一类正演问题(电法和重、磁法勘探问题),也可计算向上延拓.特别是将它与有限元法结合,解决复杂边值问题,可取得较理想的效果. 边界元法的基本计算方法有直接法和间接法两种. 边界元法的直接法现以磁法勘探正问题为例来说明. 解磁法正问题可归结为在求解区域Ω上解拉普拉斯方程的边值问题,并可用边界元法求解, 相似文献
14.
三维地形大地电磁场的边界元模拟方法 总被引：1，自引：0，他引：1
提出了一种用边界元法计算大地电磁场三维地形影响的数值模拟方法.首先用矢量积分理论和电磁场边界条件, 将上半空间(空气)和下半空间(地下介质)两个区域电磁场边值问题变为仅对地形界面的两个矢量面积分方程, 其中一个计算磁场, 称磁场方程; 另一个计算电场, 称电场方程.然后将对地形界面的积分剖分为一系列的三角单元积分.在三角单元积分中, 假设单元中电磁场为水平均匀大地空间电磁场与地形影响的迭加, 并假设地形影响为常项, 这样既保证了计算精度又使得计算方法简便.通过分解和计算, 每一个矢量面积分方程分解为对应3个坐标方向的3个常量线性方程, 这些线性方程组成了对角占优的线性方程组, 可用SSOR方法求解.文中给出了2个三维地形上大地电磁视电阻率曲线的计算结果. 相似文献
15.
点源二维起伏界面二层介质边界元解 总被引：3，自引：0，他引：3
本文以二层介质下界面上隆模型为例,讨论了用边界元法在解决层状介质中存在二维异常体情况下点源场视电阻率计算问题。文中对基本解的选取避免了对水平层状界面的积分计算,从而大大减小了所需的边界剖分点和未知量的数目,大量地节省了计算机内存,改善了计算精度,使其从解决二层介质中横向不均匀性问题推广到解决三层以上介质横向不均匀性的问题。结果表明,边界元法对本问题的研究是成功的。 相似文献
16.
介绍了联合运用积分方程法与有限元法(简称混合法)来计算三维非均匀介质中电磁场分布的理论方法,并进行了数值模拟实验。混合法的原理是引入一个包围非均匀目标体的虚构边界,在边界内部的场用有限元法模拟,在边界外部(包括边界)的场用积分方程表达,二者在边界上通过场的连续性耦合起来。数值实验结果表明,混合法既能显著地减小网格规模,又能灵活地模拟复杂的介质情况,且计算精度较高。 相似文献
17.
广义逆函数值Pad逼近(GIPA)被引入来求解积分方程,而广义逆函数值ε 算法是计算GIPA的一种有效递推算法。本文用函数值ε 算法求出GIPA的分母多项式的零点,然后取它的实部作为积分方程特征值的近似值,与其它算法进行了实例比较,显示出估计积分方程特征值的ε 算法的优越性。 相似文献
18.
基于边界元法的边坡矢量和稳定分析 总被引：4，自引：0，他引：4
矢量和法物理力学意义明确,计算简单,且能根据边坡当前的应力分布状态合理地评价其整体稳定性状态。其中边坡的应力状态通常是采用有限元法来求解。由于边界元法具有研究问题降阶、离散化带来的误差值仅产生在边界以及计算量小等优点,在工程中得到了广泛应用;对于平面问题,以源点作为原点,以所积分单元的切向和法向为坐标轴建立局部坐标系,对于线性单元可以得到所有积分的解析解。因此,可以得到计算区域内部任意点的场变量的解析解,这就保证了位于边界附近区域场变量的精度。利用边界元法得到二维边坡体内连续的应力分布状态,使用矢量和法对该边坡进行稳定性分析,并且与基于有限元的矢量和法、极限平衡法进行对比分析。边坡圆弧滑面和折线滑面的计算结果表明,基于边界元法得到的矢量和安全系数和基于有限元的矢量和法、极限平衡法基本一致;边界元法对应的矢量和安全系数对边界单元尺寸不敏感。 相似文献
19.
扶永铭 《成都理工学院学报》1995,22(3):109-114
该文研究了低频人工时变场（似稳场）中任意导电导磁矿体的边界元素法。通过对矢量位的边界积分方程，再以样条边界元素法进行离散，得到实分量的常系数方程组和虚分量的二阶微分方程组。此方程可利用时间积分的样条函数方法求解，解出适量位后，利用磁场强度与矢量位Ａ的相应关系而求得二次磁场的水平分量和垂直分量。 相似文献
20.
利用格林公式及点源场的边界条件,将地下半空间位函数的体积积分变为地面位函教的积分方程以后,用边界元方法解此积分方程便可得到地面各节点的电位值。其计算精度和速度取决于剖分方式和节点密度。本文采用变网格的剖分方式,用较少的点实现了计算,并兼顾了精度和速度。 相似文献