论文部分内容阅读
推导了同时考虑电阻率与磁导率变化的大地电磁二维方程,并应用有限单元法进行数值模拟.为了提高计算精度与效率、简化计算节点生成,采用格林定理处理二次场方程源项,并设计实现了一种基于二叉树结构的收缩网格剖分算法;采用基于最少填入元思想的稀疏矩阵符号分析方法,实现了稀疏线性方程组的LDLT求解.利用二次场算法进行模型试算,结果表明所采用的新计算方法大大减少了计算单元数量,提高了计算精度与效率.