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

图2-6 单元中的棱边号次序

2.单元插值

用六面体单元对整个区域进行剖分。对于图2-6的立方体单元, 矢量插值基函数即为式(2.48)。

3.单元分析

将式(2.30)中的区域积分分解为各单元积分之和

其中EFGH表示图2-4中的三维空间中的底面。考虑插值函数(2-47)式。将式(2.51)代入式(2.24), 有

E=NT(Exex+Eyey+Ezez)

δE=NT(δExex+δEyey+δEzez)    (2.52) 

其中δEx, δEy, δEz分别是Ex, Ey, Ez的变分。

将式(2.52)代入式(2.51)的第一项中积分

∫e×δE·×EdΩ=δETeK1eEe    (2.53) 

其中

K1e=∫e(×Ni)×(×Nj)dΩ    (2.54) 

或者

 

读书导航