共查询到20条相似文献,搜索用时 171 毫秒
1.
2.
3.
三维地面断面电阻率测深有限元数值模拟 总被引:6,自引:0,他引:6
用有限单元法进行了电导率分块均匀的三维点源电场电阻率测探的数值模拟。首先给出了三维构造中点源电场的边值问题、变分问题;然后,用有限单元法求解变分问题,将区域剖分成六面体单元,在单元中进行三线性函数插值,将变分方程化为线性代数方程组;最后解方程组,得各节点的电位值,进而计算出地表的视电阻率。对几例较典型的地电模型进行试算,结果表明本方法是行之有效的。 相似文献
4.
三维地形大地电磁场的边界元模拟方法 总被引:1,自引:0,他引:1
提出了一种用边界元法计算大地电磁场三维地形影响的数值模拟方法.首先用矢量积分理论和电磁场边界条件, 将上半空间(空气)和下半空间(地下介质)两个区域电磁场边值问题变为仅对地形界面的两个矢量面积分方程, 其中一个计算磁场, 称磁场方程; 另一个计算电场, 称电场方程.然后将对地形界面的积分剖分为一系列的三角单元积分.在三角单元积分中, 假设单元中电磁场为水平均匀大地空间电磁场与地形影响的迭加, 并假设地形影响为常项, 这样既保证了计算精度又使得计算方法简便.通过分解和计算, 每一个矢量面积分方程分解为对应3个坐标方向的3个常量线性方程, 这些线性方程组成了对角占优的线性方程组, 可用SSOR方法求解.文中给出了2个三维地形上大地电磁视电阻率曲线的计算结果. 相似文献
5.
《物探化探计算技术》2016,(3)
对三维场源二维地电模型的正演计算称为2.5维数值模拟。这里从麦克斯韦方程组出发,1分别求解电磁场的一次场和二次场,将三维场源降为二维;2利用傅氏变换将空间域方程转化为波数域,应用有限元求解波数域电磁场方程,引入无限元解决无穷远边界的收敛问题;3根据电磁场值实虚部的曲线特征,按对数等间隔选取21个波数。编写代码计算均匀半空间与解析解结果对比,电阻率的均方相对误差均小于0.5%,证明其有效性。计算了三种地电模型的电磁场响应,对单一高、低阻体和高低阻组合体的模拟效果真实,异常中心位置基本吻合。结果证明,这里的CSAMT2.5-D正演算法可以模拟较为复杂的地电模型,并取得良好的效果。 相似文献
6.
径向基函数点插值无网格法(radial point interpolation method,RPIM)是一种新型的无网格法,其形函数具有插值特性,且形式简单,易于施加本质边界条件。文中介绍了径向基函数点插值无网格法的基本原理,推导了三维情况下点插值无网格法的基本公式。从变分原理出发,结合比奥固结理论,建立了流-固耦合的三维点插值无网格法基本方程和数值积分方法,并开发了相应计算程序。通过三维悬臂梁和单向固结问题的数值试验,验证了该方法对三维弹性问题和流-固耦合问题的适用性和有效性 相似文献
7.
8.
9.
水平地形三维电场的边界单元解法 总被引:1,自引:0,他引:1
用格林公式将水平地形三维电场的边值问题转变为不均匀体表面的边界积分方程,然后用三角单元对不均匀体的表面进行剖分,并用高斯求积公式将积分方程转变为线性代数方程组。解方程组,得不均匀体表面的电位及其法向导数。根据电位和法向导数,可用数值方法计算地面的电位和电场强度。 相似文献
10.
11.
有耗媒质中线天线激励下三维异常体的散射 总被引:1,自引:0,他引:1
本文研究均匀有耗媒质中,在半波长偶极天线的场激励下,多种形状和参数的异常体的散射特性。用矩量法分别求解良导电球体散射的磁场积分方程和有限导电的长方形异常体散射的电场积分方程。讨论了与数值法有关的近似、误差和检验等问题。通过数值计算得出媒质电参数改变、散射体大小、形状以及空间位置改变时散射场变化的规律。水中金属球和介质块散射的测量结果与数值计算结果吻合较好,验证了数值模型的正确性。 相似文献
12.
贮量集中有限元方法在渗流计算中的优点 总被引:2,自引:0,他引:2
普通有限元方法和有限差分方法是渗流计算中常用的数值方法,而贮量集中有限元方法则很少被实际应用。为了使渗流计算者全面了解贮量集中有限元方法,并能够较多地应用这种方法,本文将渗流计算中的贮量集中有限元方法、普通有限元方法和有限差分方法进行了对比,具体比较和分析了有限元方程组的系数矩阵特征、计算工作量,通过具体的算例分析了它们的计算精度,由此指出:贮量集中有限元方法的计算精度与有限差分法基本相同,而远好于普通有限元法,它的计算工作量小于其他二种方法 相似文献
13.
电导率分块连续变化的二维MT有限元模拟(Ⅱ) 总被引:1,自引:0,他引:1
用有限单元法进行电导率分块连续变化的二维大地电磁场的模拟,本文的第一部分给出有限单元法系数矩阵的计算公式,这部分给出四个模型的计算结果。在有限元中,分别采用双线性播值和双二次插值,计算结果与解析法或有限差分法进行对比,吻合得很好。算例表明,本方法具有计算精度高、网络剖分简单、模拟功能强的优点。 相似文献
14.
基于球坐标系下有限差分的地磁测深三维正演 总被引:2,自引:0,他引:2
为了计算全球尺度电磁感应的响应,本文介绍地磁测深频率域三维正演。正演算法采用球坐标系下的交错网格有限差分方法,从Maxwell方程的积分形式出发,采用PARDISO对离散后的方程组求解,避免了迭代求解的散度校正。为了验证本文结果的正确性和精度,与前人的有限元和有限差分方法进行了对比,一维层状模型的三维交错网格有限差分数值结果和解析解相对误差小于5%,双半球模型的计算结果与前人的计算结果完全吻合。三维"棋盘模型"计算表明磁场分量对异常体的大小和位置具有很好的分辨能力。 相似文献
15.
电导率分块连续变化的二维MT有限元模拟(Ⅰ) 总被引:3,自引:1,他引:3
用有限元法进行电导率分块连续变化的二维大地电磁的模拟。本文第一部分给了有限元法系数矩阵的计算公式,可用于编制有限元法的计算程度。模型上的计算结果和计算精度的分析交城本文的第二部分介绍。 相似文献
16.
为提高电磁波场数值计算精度,对时间域有限差分法与时间域伪谱法进行了对比研究。时间域有限差分法是一种目前流行的电磁场时域数值计算方法,已被广泛应用于求解与时间有关的偏导数方程。对于大规模数值计算,时间域有限差分法需要较多的内存空间。时间域伪谱法基于时间域有限差分法,该方法使用快速傅里叶度换来计算麦克斯韦方程中的空间导数。由于傅里叶变换的准确性,时间域伪谱法使计算精度提高,数值计算时所需的格子数大大减少,这极大地节省了计算机内存空间,适合于大规模正反演问题的数值计算。 相似文献
17.
18.
19.
对三维问题的分析是数值流形方法发展的必然,在数值流形方法覆盖位移函数的基础上构造了一种六面体有限覆盖的三维流形单元,推导了相应的应变矩阵、刚度矩阵及平衡方程等表达式。同时,由于目前数值流形方法的模拟分析主要是采用线弹性模型,而对于非线性模型分析研究很少;根据数值流形方法的特点和岩土体的本构模型,给出了适用于数值流形方法进行非线性分析的算法。该方法利用中点增量法进行求解,以改变 模型和 模型中弹性模量的方式来反映非线性,其实质是用分段线性来取代非线性。通过地基沉降计算算例表明,数值流形方法在三维岩土体中进行非线性分析中是有效的。 相似文献
20.
基于异常复电位2.5维CR有限元数值模拟 总被引:3,自引:0,他引:3
论述了基于异常复电位2.5维复电阻率法的正演计算技术。首先给出了二维地电断面、三维交流双点电源复电阻法的变分问题;然后利用有限单元法求解变分问题,在矩形网格中增加2条对角线形成三角网格对研究区域进行剖分,在单元中对复电位和岩、矿石复电导率进行线性插值,将变分方程化为线性代数方程组;最后,求解线性方程组,得到各节点的异常复电位值,并根据选定的测量装置计算出地表的视复电阻率,进而计算出复电阻率的振幅、相位或复电阻率的实、虚分量。通过理论模型的计算检验算法的有效性后构建模型算例,给出了其不同测量装置下的计算结果。 相似文献