论文部分内容阅读
在当前的有限元三维直流电阻率法中数值模拟中,存在两个明显的问题。第一个问题是,由于采用正方体或六面体等规则单元所引起的几何模型离散化误差不能够被有效地消除,其会引起较大的误差存在于最后的视电阻率数值解中,尽管有一些学者曾尝试着用一些四面体网格来消除这些较大的误差,但是其采用的方法仅仅限于通过划分的每个六面体网格成为几个四面体单元来形成简单的四面体网格,就其本质而言,这些网格还是基于规则的六面体块结构的,其不可避免的会产生较大的误差在曲面边界等不规则边界处;第二个问题是,标准的有限单元法有非常可靠的精度收敛理论,但是其仅仅是通过要求有限单元无限小且有限元单元形状函数足够高来得到保证.显而易见地,这种无限的要求必要造成巨大的计算消耗,这种过大的消耗使得三维的直流电阻率模拟变成不实现的事情。特别地,对于复杂的无解析解的三维模型来说,可靠的数值解就不可能由这些无限要求来保证。在本文中,我们呈述了一种全新的基于非结构化网格的自适应有限元数值模型方法。非结构化网格的引入可以有效地消除由结构化网格产生的离散化误差;而迭代的自适应局部网格加密技术,使得我们的迭代数值解渐近地逼近未知的真实值。每一步网格的近似最优化,加快了迭代过程的收敛性,具有单元形状质量保证的同时,网格内节点与单元呈梯度分布,这两种性质即保证了迭代过程的快速收敛,降低了计算时间又保证了计算工作量的降低。非结构化网格是由Delaunay加密算法自动生成,每一个单元的形状的质量有所保证,使由单元形状引起的误差降至最低。自适应迭代算法为h型的自适应算法,因为其执行简单并且具有高收敛性。我们采用Zienkiewicz andZhu提出的后验误差估计技术与局部单元自适应加密或加粗技术来进行迭代过程,其核心技术为梯度恢复技术。实验数据表明,经过4-5步迭代后,我们可以得到相对误差小于0.5%-1.0%的高精度结果,并且电源点的奇异性可以完全有效地消除。由此可见,本文提出的自适应有限元算法不仅能够大幅度提高数值度的精度,而且还可以消除直流电场内的奇异性现象,并且大幅度节省计算消耗,从而使本文的技术可以PC个人计算机上进行运行,有效地缓解计算精度与计算时间的矛盾,从而更好的服务于直流电阻率勘探。