共查询到19条相似文献，搜索用时 265 毫秒
1.
首先从麦克斯韦方程出发,用伽里金方法推导了三维三分量CSAMT法的有限元方程.在研究过程中,认识到加入散度条件的必要性,在公式中强加了散度条件,提高了解的完备性.其次将成功应用于二维线源频率域电磁法有限元模拟中的两种技术推广到三维中,一是边界条件统一采用一阶吸收边界,使线源产生的电磁波在边界上按波的传播规律被吸收,以降低平面波假设造成的影响;二是总体系数矩阵的存储,用两个二维数组分别记录总体系数矩阵的非零元素及其在总体结点编号中所处的位置,使总体系数矩阵的存储量达到最小的同时,物理意义明确,迭代求解时迅速简便.最后用均匀半空间模型进行了验证. 相似文献
2.
3.
4.
波动数值模拟的稳定性是获得可靠结果的前提.透射边界是一类具有高阶、高效等特点的人工边界,其引发的高频失稳是由内域格式和透射边界的不当耦合所致.本文针对P-SV波动有限元模拟中透射边界引发的失稳问题,基于GKS定理的群速度解释,通过对有限元和透射边界的频散分析揭示了数值失稳机理为透射边界和相邻内域格式支持了群速度指向内域的高频P波或SV波,波动能量将从边界进入内域引发数值失稳.同时,对比连续模型频散指出引发失稳的谐波是由有限元离散引入.本文采用修改的数值积分方法调整有限元刚度,以消除有限元中引发边界失稳的高频波动成分,从而稳定实现透射边界.理论分析和数值实验均表明本文稳定措施的有效性. 相似文献
5.
频率域可控源电磁是在大地电磁测深的基础上发展起来的一种人工源电磁测深法,其二维电磁响应的计算须采用数值模拟方法.本文以Matlab为程序编译工具,采用双二次插值的有限单元法,推导出相应的计算公式.为了模拟无穷远边界及满足计算机的内存需求,在保证计算精度的情况下设计了非均匀网格剖分.在程序编制中,只存储有限元系数矩阵的非零元素,大大减少了正演计算的时间.针对频率域可控源电磁法中卡尼亚电阻率在过渡区和近区畸变的问题,给出了全区视电阻率的迭代公式,并对典型的一维层状模型以及简单二维模型进行了计算.过渡区和近区数据经过校正后,可以正确反映出模型的地电特征,证明了线源下近区勘探的可能性. 相似文献
6.
复频移完全匹配层(Complex Frequency-Shifted PML,CFS-PML)在长时间时域计算中对凋落波、倏失波具有好的吸收效果,并被广泛应用于时域有限差分模拟中.而本文采用卷积方法将CFS-PML应用于时域有限元求解GPR波动方程的数值模拟中.论文以TM波为例,推导了基于CPML(Convolutional PML)边界的时域有限元GPR波动方程求解公式,采用Newmark-β方法对时间导数进行离散,有效改善了时域有限元GPR数值计算程序的稳定性.并以狭长模型为例,开展了CPML边界中关键参数m、R和κ的选取实验,通过对比反射误差大小确定了综合最优参数组合.相同时刻UPML与CPML波场快照、3个检测点的反射误差比较,说明CPML较UPML具有更好的吸收效果.最后,采用非规则四边形网格对1个复杂GPR模型进行剖分,应用加载CPML边界条件的FETD程序对该模型进行了正演,得到了二维剖面法、宽角法正演GPR剖面图,说明非规则四边形对复杂模型的良好适应性,基于CPML边界条件的FETD可有效减少边界反射误差,能实现对任意复杂不规则模型的正演模拟. 相似文献
7.
基于双二次插值的探地雷达有限元数值模拟 总被引：3，自引：0，他引：3
从探地雷达(GPR)满足的波动方程出发,详细介绍了二维GPR模型单元剖分、二次插值、数值积分和有限元刚度矩阵总体合成的GPR有限元求解过程.为解决数值模拟时截断边界处的超强反射,采用Clay Bout透射边界条件对雷达波进行衰减,进而压制了来自截断边界处的反射波.在满足时间步长与空间网格差分稳定性前提下,采用中心差分法对GPR有限元方程进行离散,并用不完全LU分解预处理的BICGSTAB算法求解系数方程组,然后编制了基于双二次插值的GPR有限元正演模拟matlab程序.运用该程序分别对矩形和"V"字形两个典型地电模型进行正演计算,得到了正演剖面图,将该正演剖面图与基于线性插值的FEM算法的正演剖面图做了对比分析.结果表明基于双二次插值FEM算法相比基于双线性插值FEM算法异常响应更明显,具有更高的模拟精度,更有利于指导雷达剖面的数据解译. 相似文献
8.
基于非结构网格有限元方法开展了三维复杂地电模型的线源井地电法的高效正演模拟研究,探讨了通过求取电场响应导数来刻画目标体边界范围、采用差异场地形校正技术来消除地形影响等措施对井地电法成像的效果和精度的影响。并通过对比与解析解,验证了本文数值解算法的有效性。模型计算结果表明：积水巷道的空间位置和走向均会引起视电阻率的显著变化,视电阻率变化率的极值准确且清晰地指示了巷道边界的位置;电位的归一化总水平导数极大地提高了井地电法对目标体复杂边界位置的识别能力;地形对井地电场分布的影响也很大,其视电阻率响应与地形形状近似呈对称关系,利用差异场技术能有效地削弱地形对井地电法高精度成像的影响。 相似文献
9.
10.
利用有限元软件Ansys二次开发语言apdl进行二次开发,将多次透射边界(MTF边界)和弹阻边界(V-S边界)添加到软件中,从而实现对土-结构动力相互作用问题中无限域场地的模拟。以野外大比例(1/2)土-箱基-框架结构试验模型的牵引释放自由振动试验(波源问题)和爆破震动试验(散射问题)数据为基础,通过有限元数值计算结果和试验结果的对比可知,施加了人工边界的土体模型可以较好地模拟波在无限域土体中的传播,有效地减少有限元计算模拟时的计算规模。同时多次透射边界在计算精度以及散射问题中地震波输入的便捷性方面要优于弹阻边界。 相似文献
11.
在多源直流电阻率法有限元三维数值模拟中,传统混合边界条件由于刚度矩阵与源点位置相关,无法实现线性方程组的快速回代求解,目前常用齐次边界条件或无限元边界进行替代,虽然实现了快速回代求解,但同时也降低了数值模拟的精度.为了实现快速回代求解,并确保数值模拟的计算精度,本文提出了一种近似边界条件方法.其核心思想是将与场源位置相关的边界系数矩阵从刚度矩阵中分离出来,使得分离后的刚度矩阵与场源位置无关.并将边界系数矩阵与边界处一次电场向量的乘积移到线性方程组右端源项中,当场源位置发生改变时,只需要重新计算右端源项就可以实现快速回代求解.理论模型数值计算表明,在水平地形条件下,本文边界条件数值精度优于混合边界条件;在起伏地形条件下,与齐次边界条件相比,本文边界条件数值结果与混合边界条件吻合度更高. 相似文献
12.
提出了一种用边界元法计算频率域人工源三维地形电磁场的数值模拟方法.首先用矢量积分理论和电磁场边界条件，将上半空间（空气）和下半空间（地下介质）两个区域电磁场边值问题变为仅对地形界面的两个矢量面积分方程.然后将对地形界面的积分剖分为一系列的三角单元积分.在三角单元积分中，假设单元中电磁场为无限大气空间电磁场与地形影响的叠加,并假设地形影响为常项,这样既保证了计算精度又使得计算方法简便.通过分解和计算，每一个矢量面积分方程分解为对应三个坐标方向的三个常量线性方程，这些线性方程组成了对角占优的线性方程组，可用SSOR方法求解.文中给出了垂直磁偶源的垂直磁场地形影响的例子. 相似文献
13.
14.
Haitjema HM 《Ground water》2006,44(1):102-105
The analytic element method, like the boundary integral equation method, gives rise to a system of equations with a fully populated coefficient matrix. For simple problems, these systems of equations are linear, and a direct solution method, such as Gauss elimination, offers the most efficient solution strategy. However, more realistic models of regional ground water flow involve nonlinear equations, particularly when including surface water and ground water interactions. The problem may still be solved by use of Gauss elimination, but it requires an iterative procedure with a reconstruction and decomposition of the coefficient matrix at every iteration step. The nonlinearities manifest themselves as changes in individual matrix coefficients and the elimination (or reintroduction) of several equations between one iteration and the other. The repeated matrix reconstruction and decomposition is computationally intense and may be avoided by use of the Sherman-Morrison formula, which can be used to modify the original solution in accordance with (small) changes in the coefficient matrix. The computational efficiency of the Sherman-Morrison formula decreases with increasing numbers of equations to be modified. In view of this, the Sherman-Morrison formula is only used to remove equations from the original set of equations, while treating all other nonlinearities by use of an iterative refinement procedure. 相似文献
15.
2.5D direct‐current resistivity forward modelling and inversion by finite‐element–infinite‐element coupled method 
下载免费PDF全文
下载免费PDF全文 Yuan Yuan Jianke Qiang Jingtian Tang Zhengyong Ren Xiao Xiao 《Geophysical Prospecting》2016,64(3):767-779
To reduce the numerical errors arising from the improper enforcement of the artificial boundary conditions on the distant surface that encloses the underground part of the subsurface, we present a finite‐element–infinite‐element coupled method to significantly reduce the computation time and memory cost in the 2.5D direct‐current resistivity inversion. We first present the boundary value problem of the secondary potential. Then, a new type of infinite element is analysed and applied to replace the conventionally used mixed boundary condition on the distant boundary. In the internal domain, a standard finite‐element method is used to derive the final system of linear equations. With a novel shape function for infinite elements at the subsurface boundary, the final system matrix is sparse, symmetric, and independent of source electrodes. Through lower upper decomposition, the multi‐pole potentials can be swiftly obtained by simple back‐substitutions. We embed the newly developed forward solution to the inversion procedure. To compute the sensitivity matrix, we adopt the efficient adjoint equation approach to further reduce the computation cost. Finally, several synthetic examples are tested to show the efficiency of inversion. 相似文献
16.
系统地论述了大地电磁三维交错采样有限差分数值模拟算法实现过程中交错网格剖分、积分公式离散化、边界条件、方程组求解、三维张量阻抗的计算等内容. 由于提出了简洁的边界条件，采用了解大型系数矩阵方程组的双共轭梯度稳定解法，所实现的三维交错采样有限差分数值模拟算法具有迭代收敛稳定、计算精度高、速度快等特点. 通过两个理论模型的计算结果检验了算法的正确性和计算精度. 所实现的三维交错采样有限差分数值模拟算法为研究三维反演问题奠定了基础. 相似文献
17.
根据弹性薄板理论和振动理论,运用边界元法(BEM)分析了薄板弯曲自由振动的频率特性,从薄板自由振动的微分方程式出发,采用动态基本解推导出其相应的边界积分方程式,计算中对边界积分方程进行缩减,应用频率扫描的方法研究薄板的动态特性,准确、有效地得到了其相应的固有频率、给出的数值算例证明了本文的方法简单、可行,且具有良好的解析计算精度。 相似文献
18.
基于连续体系模型与传递矩阵法的结合,提出了一种简化方法来计算多肢剪力墙的动力特性。整个结构体系等效模拟为一个夹层结构梁;建立了相应的基本微分方程;通过求解微分方程得到了楼层间变形函数;利用边界条件及层间变形函数导出了层间传递矩阵;据此计算出了体系的振型与自振频率。算例分析表明本文方法计算结果与文献给出的解是比较吻合的。 相似文献
19.
The scaled boundary finite‐element method is extended to simulate time‐harmonic responses of non‐homogeneous unbounded domains with the elasticity modulus and mass density varying as power functions of spatial coordinates. The unbounded domains and the elasticity matrices are transformed to the scaled boundary coordinates. The scaled boundary finite‐element equation in displacement amplitudes are derived directly from the governing equations of elastodynamics. To enforce the radiation condition at infinity, an asymptotic expansion of the dynamic‐stiffness matrix for high frequency is developed. The dynamic‐stiffness matrix at lower frequency is obtained by numerical integration of ordinary differential equations. Only the boundary is discretized yielding a reduction of the spatial dimension by one. No fundamental solution is required. Material anisotropy is modelled without additional efforts. Examples of two‐ and three‐dimensional non‐homogeneous isotropic and transversely isotropic unbounded domains are presented. The results demonstrate the accuracy and simplicity of the scaled boundary finite‐element method. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献