伽遼金法加權余量積分式的推導, 比用廣義變分原理推倒變分方程簡單得多, 所以本書中三維情況下的大地電磁變分方程采用了伽遼金方法。但是, 伽遼金方法有一個缺點, 微分算符￡出現(xiàn)在式(2.30)中的積分號內, 與變分法的泛函相比, 具有更高階的微商, 形函數(shù)應該具有較高的階次的連續(xù)性, 這往往給形函數(shù)的構造帶來極大的困難, 為了避免上述困難, 可以對式(2.30)進行分部積分, 或應用格林公式, 來降低微分項的階次。

2.2.2 三維MT邊值問題的Galerkin有限元方程

式(2.15)和(2.16)定義了MT問題的邊值問題。借助于矢量恒等式B·(×A)=A·(×B)+·(A×B)和Gauss分部積分公式, 其對應的變分為

b(E, V)=f(V), V∈H(curl)    (2.31)

上式中, H(curl)={V‖V∈L2(Ω), V×n=n×V0}, b為雙線性表達式, f 為單線性表達式: 

b=∫Ω(×E·×V-iωμσE·V)dΩ    (2.32)

如果在此就考慮方程(2.16)所示的Dirichlet邊界條件作為本書的源項, 右端項f可表示為: 

f=∮ΩV·E0dΓ    (2.33)

式(2.33)表示在邊界上對已知的電場E0作邊界積分。

為了求解方程(2.32)所表示的電場分布, 本書在此采用矢量有限元的方法。在矢量有限元法中, 需要首先把整個求解區(qū)域離散化成一系列的六面體單元。與節(jié)點型有限元不同之處在于, 在網(wǎng)格六面體單元中的每條邊的切向上定義單標量的切向電場值E。考慮任意的六面體單元e, 在e中電場的近似表達式為: