复杂条件下多面体单元直流电三维有限元数值模拟.pdfVIP

10血Chi眦International 2011 The Geo-ElectromagneticWorkshop 复杂条件下多面体单元直流电三维有限元数值模拟 张钱江1，戴世坤2 中南大学 l引言 ， 对于复杂条件下地球物理问题的有限元数值模拟，单元的选取需要满足两个基本条件：1)能够模拟出任 意物性复杂分布的地质模型；2)能够计算出单元积分的解析解。根据这两个要求我们构思出了多面体单元有 限单元法，它在不增加节点个数和计算量的前提下，能有效的提升有限单元法模拟复杂条件下三维地质模型 的能力和精度。 在有限单元法算法中我们采用二维数组对多面体单元节点线性相关的非零元素进行了压缩存储，有效的 减少了计算机内存和计算量，提高了计算效率和并行能力。 在解线性方程组的过程中，我们分别采用了迭代法和直接法进行求解。其中迭代法主要为带预置条件的 共轭梯度法，其中预优矩阵主要有：对角线阵，三对角阵，快速松弛法(SSoR)，带修正的不完全CHOLESKY 分解，不完全LDLT分解；直接法主要是基于SUPERLU分解原理的快速求解方法。通过对多种方法计算时问和 并行能力的分析对比，总结出每种方法的适用范围以及求解三维直流电有限元法的最佳方案。 最后通过具体算例验证多面体单元有限元法算法的准确性并讨论其实用性。 2三维直流电场问题 2．1边值问题 将双点电源置于水平地面，电位的边值问题可归纳如下： EQ V·(oVu)=-2z,￡(a)+2zS(B) ％=o ∈L (1) ∈L 哆乞+夕乙-rD[％∞双，=．，一，一么∞双白，一，=。 式子中A，B为地面上双电源位置；，=．和％分别是测点到电源点A，B的距离；L为区域Q的地面边 界；L为区域Q的地下边界；厅为边界的外法向方向；盯为介质的电导率；u为电位。 2．2变分问题 与边值问题(1)等价的变分问题为： F(材)=』[吉印”)2—2叫一-+2叫口-]m+％f．％一‘) 阮咄¨)_r，An：c。s(rn冲订 (2) 留0)=o 3多面体单元有限单元法 多面体单元有限单元法是指用多面体单元对整个区域进行离散划分。多面体单元如图lb所示。由十二个 四面体单元合并成，共有九个节点，其中第九个节点为单元内部公共节点。它具有如下四个优点：1)单元外 部八个节点随着节点坐标的改变，多面体单元形状也相应发生变化，可以是六面体单元到十二面体单元之间 的任意多面体，从而能模拟出任意复杂的地质模型。2)单元内部公共节点的存在，赋予了多面体单元内存在 物性界面的能力，提升了对于复杂物性界面的模拟能力。3)多面体单元内部公共节点的位置可以任意变化， 适用于具有多层物性界面的复杂地质模型。4)多面体单元由四面体单元按照一定的规则组装而成，我们可以 计算出四面体单元积分的解析解，因此也能得到多面体单元的解析解。 第十届中国国际地球电磁学术讨论会2011 7 l 3 8 2 (a》 (b)