有限单元法非均匀剖分对数值模拟影响研究  

The influence of FEM non-uniform division on numerical modelling

作　　者：程勃[1,2] 黄俊革[1] 阮百尧[2] 

机构地区：[1]上海应用技术学院土木建筑与安全工程学院,上海200234 [2]桂林工学院资源与环境工程系,桂林541004

出　　处：《地球物理学进展》2010年第5期1856-1862,共7页Progress in Geophysics

基　　金：国家自然科学基金项目(40774057);上海市教委创新基金资助(08YZ168);上海应用技术学院科研启动基金资助

摘　　要：本文研究了有限单元法非均匀剖分对直流电阻率法数值模拟结果的影响.通过非均匀剖分对半空间地质体的视电阻率异常与均匀剖分的视电阻率对比,分析这种模拟方法的可行性,研究计算结果差异产生的原因.通过相邻单元的不同边长比例剖分的对比,研究在误差允许范围内非均匀剖分的单元边长比例极限.计算结果表明,测线方向非均匀剖分是误差产生的主要原因,XYZ三方向非均匀剖分所带来的误差总和是各方向非均匀剖分产生误差之和;当测线方向剖分单元不均匀剖分比例小于0.49∶0.01时,由于非均匀剖分产生的误差在容许范围内.A Study on influence of numerical modeling with non-uniform division of FEM on DC resistivity is done in the paper. Through the compare between multi-dimensional non-uniform division on the half-space subdivision of the geological apparent resistivity anomalies and evenly subdivision way, analysis the feasibility of numerical stimulation and the causes for the differences of simulation results. Through the contrast of different element deviation inferred the limitation of non-uniform division. The results show that non-uniform division upon the profile line is the main reason of error. The error of three dimensional non-uniform division is generally the sum of each single direction. When the ratio of non-uniform division upon profile line is less than 0.49∶0.01, the error caused by non-uniform division is within the range of tolerance.

关 键 词：有限单元法 数值模拟 非均匀剖分 视电阻率 

分 类 号：P631[天文地球—地质矿产勘探]