大地电磁二维正演中的有限元-径向基点插值法

作者：李俊杰; 严家斌

来源：中国有色金属学报, 2015, 25(05): 1314-1324.

DOI：10.19476/j.ysxb.1004.0609.2015.05.026

• 径向基点插值法
• 有限元
• 无网格方法
• 大地电磁

径向基点插值法(RPIM)作为一种插值型无网格方法,为改善无网格点插值法(PIM)在形函数构造过程中可能出现的矩阵奇异性问题而提出的一种方法,该算法支持域无量纲尺寸的选择区间大,能更好地处理各类工程与科学计算问题。介绍了RPIM的近似原理,给出了径向基函数形状参数的推荐值;从大地电磁二维变分问题出发利用Galerkin法结合高斯积分公式推导出相应的系统矩阵离散表达式;为提高RPIM的计算效率,将RPIM与有限元法(FEM)耦合,提出了有限元-径向基点插值法(FE-RPIM),多个模型的数值计算验证了RPIM精度高、处理复杂模型便利及耦合法计算复杂模型高效的特点。

• 单位
  浙江省水利水电勘测设计院; 有色资源与地质灾害探查湖南省重点实验室; 物理学院; 中南大学