原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
中国科技论文在线
基于局部加密等级网格的2.5 维直流电法有
限元模拟#
*
潘克家，汤井田，肖晓
5 （中南大学有色金属成矿预测教育部重点实验室, 地球科学与信息物理学院, 长沙 410083 ）
摘要：合理截取半圆形计算区域，采取局部加密的 等级网格，结合对称行索引存贮格式
(CSR)及并行稀疏求解器PARDISO ，提出一种高效、高精度的2.5 维直流电阻率法有限元正
演方案，并编制了相应的Fortran 程序，对具有解析解的三个典型地电模型进行了计算与详
10 细的分析。研究结果表明：圆形截断边界不仅便于在径向方向上采取 等级剖分，而且能
大大简化有限元模拟中单元刚度矩阵的计算。结构化的等级网格避开了通常非结构化网格有
限元计算时繁琐的网格剖分及总体刚度阵的集成过程，且能在不增大问题规模的前提下，显
著提高2.5 维直流电法正演问题源点附近的模拟精度。Intel MKL 的PARDISO 求解器能在
普通PC 机上5 秒内求解电法正演有限元离散得到的100 万阶稀疏线性方程组，可广泛用于
15 各种地球物理正、反演问题。
关键词：直流电阻率; 等级网格; 有限元; 局部加密; PARDISO; 行索引存贮格式
中图分类号：P631
Finite element modeling of 2.5-D DC resistivity based on
20 locally refined graded mesh
PAN Kejia, TANG Jingtian, XIAO Xiao
(Key Laboratory of Metallogenic Prediction of Nonferrous Metals, Ministry of Education/School
of Geosciences and Info-Physics, Central South University, Changsha 410083)
Abstract: Truncating semi-circular computational domain reasonably, and combining the compressed
25 sparse row (CSR) format and the parallel direct sparse solver known as the PARDISO solver, an
efficient, high-precision finite element method for 2.5-D DC resistivity modeling is proposed based on
locally refined graded mesh, and corresponding FORTRAN program are developed. The numerical
results and detailed analysis are given for three typical geoelectric models which have analytical
文档评论（0）