三维直流电阻率有限元-无限元耦合数值模拟

为解决传统有限元截断边界所引起的问题,本文提出了一种新的三维直流电阻率有限元-无限元耦合数值模拟方法.首先推导了无限元三维单元映射函数,然后提出了一种全新的最优的无限元形函数并与多种其他形函数进行了对比,随后将其与非结构化四面体有限元相结合,取代了传统的混合边界条件,使得电位在无限域内连续并在无限远处衰减为零,最终形成的左端矩阵稀疏对称并与场源位置无关.数值计算表明,该方法可以在近似测区大小的计算范围内得到与混合边界条件相当的计算精度,优于相同计算范围下齐次边界条件的解,有利于减少计算节点数;由于左端矩阵不随场源位置改变,有利于加速反演计算.