基于非结构矢量有限元法的地下UXO探测模拟_毕业论文.pdf

I

基于非结构矢量有限元法的地下 UXO 探测模拟

摘 要

未爆弹（UXO）是一种已武装、发射且因故障而未爆炸的弹药，主要来源于

历史战争或者军事活动，包括子弹、手榴弹、炮弹、炸药、地雷和导弹等。这些

未爆弹往往埋藏在地下较浅位置甚至出露于地表之上，给百姓的人身财产安全带

来巨大威胁，阻碍了社会和经济的发展。据联合国报道，全球未爆弹埋藏量超过

本文基于非结构矢量有限元法进行地下 UXO 探测模拟，分析多种模拟参数

对 UXO 响应特征曲线的影响，从中发现一定规律，为反演计算和野外探测提供

参考依据。非结构四面体网格有效地刻画了复杂 UXO 形态和不同介质分界面，

可灵活地对局部网格加密。矢量插值基函数的使用能够避免伪解的出现，二次场

分离法能够解决源的奇异性问题。本文从电磁场的基本方程麦克斯韦方程组出发，

推导出有限元方法的弱形式方程，再利用 Galerkin 法对控制方程进行有限元离

散，采用四面体矢量单元的插值基函数，获得每个单元的局部方程，接着将局部

方程组装成大型稀疏全局矩阵，利用直接求解法计算目标的正演响应。本文基于

约束 Delaunay 三角剖分算法构建约束 Delaunay 四面体，能够保证网格质量，对

网格局部加密处理，然后设计了空气、地下空间、天线系统以及八种不同形态的

II

控件回调函数进行数据输入输出，最后将运行无误的代码进行打包，形成可下载

的 EXE 软件包。这种方法提供了开发图形用户界面的一种思路，可为其他基于

MATLAB 的软件开发工作带来一定的参考意义。

本文将非结构网格和矢量有限元法引入 UXO 探测之中，成功实现了正演响

应计算，并对各种探测参数进行了模拟研究，其结果为野外探测提供了理论依据，

同时也为反演工作奠定了基础。

关键词：

频率域电磁法，非结构网格，矢量有限元，UXO 正演模拟，MATLAB 可视

化界面软件开发

III

including bullets, grenades, artillery shells, explosives, mines and missiles. These

unexploded bombs are often buried in shallow locations underground or even exposed

above the surface, posing a great threat to the safety of people and property and

hindering social and economic development. According to the United Nations, there are

more than 100 million unexploded bombs buried worldwide, killing tens of thousands

of people every year and contaminating a large amount of land with unexploded bombs.

Therefore, the detection and precise exclusion of underground UXO is of great

significance. UXO detection technology based on electromagnetic induction has

existed for many years and can serve all phases. However, the current electromagnetic

induction simulation generally only considers the response of the metal target body,

ignoring the influence of the surrounding media, and lacks research and discussion of

parameters on the characteristic curves of UXO response, and find certain laws from

them to provide a reference basis for inversion calculation and field detection. The

unstructured tetrahedral mesh effectively portrays the complex UXO morphology and

different media partition interfaces, and can be flexibly encrypted for the local mesh.

The use of vector interpolation basis functions can avoid the appearance of pseudo-

solutions, and the quadratic field separation method can solve the source singularity

IV

problem. In this study, the weak form equations of the finite element method are

developed based on Maxwell equations of electromagnetic field. Then Galerkin's

method is utilized to discretize the control equations for finite element analysis and the

interpolation basis functions of tetrahedral vector cells are used to obtain the local

The results demonstrate that the algorithm exhibits excellent horizontal positioning

resolution. Meanwhile, this paper also studies the effects of UXO attitude, transmit

signal frequency, soil conductivity, magnetic permeability and coil arrangement,

compares the differences of UXO model response under different parameters, and

analyzes various influencing factors and their characteristic laws. This paper also

carries out the technical flow, structural design and visual interface design of software

development based on the App Designer function of MATLAB computing software,

using a series of visual interface design tools and components to complete the interface

construction of the application, using control callback functions for data input and

output, and finally packaging the running error-free code to form a downloadable EXE

theoretical basis for field exploration, and also lay the foundation for inversion work.

V

Key words:

Frequency domain electromagnetic method, Unstructured grid, Vector finite

element, UXO forward simulation, MATLAB visual interface software development.