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

第25卷第 5期 中国有色金属学报 2015年5月 Volume25Number5 TheChineseJournalofNonferronsMetals M av 2015 文章编号：1004．0609(2015)一05—1314—11 大地电磁二维正演中的 有限元一径向基点插值法 李俊杰 ，严家斌 (1．浙江省水利水电勘测设计院，杭州 310002； 2．中南大学 地球科学与信息物理学院 有色资源与地质灾害探查湖南省重点实验室，长沙 4100831 摘 要：径向基点插值法(RPIM)作为一种插值型无网格方法，为改善无网格点插值法(PIM)在形函数构造过程中 可能出现的矩阵奇异性问题而提出的一种方法，该算法支持域无量纲尺寸的选择区间大，能更好地处理各类工程 与科学计算问题。介绍了RPIM 的近似原理，给出了径向基函数形状参数的推荐值；从大地电磁二维变分问题出 发利用Galerkin法结合高斯积分公式推导出相应的系统矩阵离散表达式；为提高RPIM 的计算效率，将RPIM与 有限元法(FEM)耦合，提出了有限元一径向基点插值法(FE．RPIM)，多个模型的数值计算验证了RPIM精度高、处 理复杂模型便利及耦合法计算复杂模型高效的特点。 关键词：径向基点插值法；有限元；无网格方法；大地电磁 中图分类号：P631 文献标志码：A M agnetotellurictwo-dimensionalforwardby finiteelement—radialpointinterpolationmethod LIJun-jie，YANJia-bin2 (1．ZhejiangDesignInstituteofWaterConservancyandHydroelectricPower,Hnagzhou310002，China； 2．KeyLaboratoryofNonferrousResourcesandGeologicalHazardDetection， SchoolofGeosciencesnadInfo·Physics，CentralSouthUniversity,Changsha410083，China) Abstract：Polynomialbasisinterpolationmethod(RPIM)，asakindoftypicalinterpolationmeshfreemethod，was proposedtoovercomethedefectsofpointinterpolationmehtod(PIM)thattheconstructionprocessoftheshapefunction involvesthematrixinverseoperation．Thismethodovercomeshtematrixinverseproblem，nadsupportsthewiderdomain dimensionless size intervalto better dealwith allkindsofengineering and scientific computing problems．The approximateprincipleofRPIM wasinrtoducedindetail，andhtediscretesystem matrixexpressioncorrespondingtohte magnetotellurictwo-dimensionalvariationalproblem bycombiningtheGalerkinmehtodandthega