第2章 大地电磁法三维矢量有限元模拟(13)

 

其中, 当i=1、 2、 3、 4时, 

不难看出, 上式定义的基函数具有零散度和非零旋度的特性。而且, 构成离散单元各个小平面上的切向场仅由组成小平面的棱边上的切向场决定。因此上式给出的展开式不仅保证了穿越棱边的切向场连续, 而且保证了穿越离散单元表面时切向场的连续性。

2.3.3 三维矢量单元分析

1.区域剖分

用六面体单元对区域进行剖分(见图2-5), 每个单元的八个角点做为节点, 并给每个单元的电阻率ρ=1/σ赋值, 立方体单元的沿x方向宽度为a, 沿y方向宽度为b, 沿z方向宽度为c, 单元的八个节点的编号和棱边编号如图2-6所示。

图2-5 网格剖分的节点和单元排列示意图

 

读书导航