Воробьев Теория электромагн поля и СВЧ (Кривець)
.pdf31
→ |
|
|
→ |
→ |
уравнения Максвелла: div D = ρсвоб . Так как |
D = εa E , а |
|||
→ |
|
|
|
|
E = −gradϕ , то |
|
|
|
|
divgradϕ = − |
ρсвоб |
. |
|
|
|
|
|
||
|
εa |
|
|
|
Оператор 2 = ∆ = divgrad |
называют |
оператором |
||
Лапласа или лапласианом. Поэтому можно записать, что
∆ϕ = − |
ρсвоб |
. |
(1.17) |
|
|||
|
εa |
|
|
Это уравнение называют уравнением Пуассона.
Если свободных зарядов ρсвоб нет, то уравнение
Пуассона переходит в уравнение Лапласа: |
|
∆ϕ = 0 . |
(1.18) |
Вдекартовой системе координат уравнения Пуассона
иЛапласа имеют вид:
∂2ϕ + |
∂2ϕ + |
∂2ϕ = − |
ρсвоб |
; |
|
|
|||||
∂x2 |
∂y2 |
∂z2 |
|
εa |
|
∂2ϕ |
+ ∂2ϕ |
+ ∂2ϕ |
= 0 . |
|
|
∂x2 |
∂y2 |
∂z2 |
|
|
|
Для цилиндрической и сферической систем координат они представлены в [3].
При получении уравнения Пуассона для векторного
→
потенциала A запишем первое уравнение Максвелла для
→ |
|
→ |
|
|
стационарного поля: rot H =δ . |
|
|
||
Умножив обе части на |
µa и используя материальное |
|||
уравнение (1.2), получим: |
|
|
|
|
→ |
|
→ |
|
→ |
µa rot H |
= rot B |
= µa |
δ . |
|
32
Воспользуемся |
выражением |
магнитной |
индукции |
||
|
|
|
|
→ |
→ |
через векторный потенциал (1.5). Тогда rotrot A = µa δ . |
|||||
С учётом того, |
|
→ |
|
→ |
→ |
что rotrot A = graddiv A− 2 |
A , можно |
||||
записать |
|
|
|
|
|
|
→ |
→ |
→ |
|
|
graddiv A− 2 A |
= µa δ . |
|
|||
Так как для стационарного |
поля |
|
→ |
||
линии вектора A |
|||||
замкнуты сами на себя, то |
→ |
В итоге получаем |
|||
div A = 0 . |
|||||
уравнение Пуассона для векторного потенциала: |
|
||||
|
→ |
|
→ |
|
|
|
2 A = −µa |
δ . |
|
(1.19) |
|
1.4 Уравнение энергетического баланса электромагнитного поля (теорема Умова-Пойнтинга)
Кроме уравнений Максвелла и закона сохранения заряда, большое значение в теории электромагнитного поля имеет теорема Умова-Пойнтинга, которая
описывает энергетические соотношения распределения полей в заданном объёме. Её не сложно получить из уравнений Максвелла.
Для вывода теоремы Умова-Пойнтинга воспользуемся первым и вторым уравнениями Максвелла в дифференциальной форме, записанными с учётом материальных уравнений (1.1), (1.2), (1.6):
rot H = γ E+ε |
|
→ |
; rot E = −µ |
|
→ |
|
|
∂ E |
|
∂ H . |
|||
→ |
→ |
|
|
→ |
|
|
|
|
a |
∂t |
|
a |
∂t |
Домножим почленно первое уравнение Максвелла
→ |
→ |
скалярно на E , а второе – на H и вычтем из первого уравнения второе уравнение. Получим
|
|
|
|
|
|
33 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
→ |
|
|
|
|
|
→ |
|
→ |
→ → |
|
→ |
= |
→ → |
|
|
→ |
∂ E |
+ |
µ |
|
|
→ |
∂ H |
. |
||
E rot H − H rot E |
γ E E+ε |
a |
E |
∂t |
a |
|
H |
∂t |
||||||||||
С учётом того, что |
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
→ |
→ → |
→ |
|
|
→ |
|
→ |
|
|
|
|
|
→ |
→ |
|
|||
E rot H − H rot E = div H×E |
= −div E×H |
, |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
можно записать |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂ |
|
→ |
|
|
|
|
|
|
|
→ |
|
|
→ → |
|
|
|
|
→ |
E |
|
|
|
→ |
|
∂ H |
|
||||
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|||||
−div E×H |
|
= γ E |
|
+εa E |
|
|
|
|
+ |
µa |
H |
|
|
|
. |
|||
|
|
∂t |
|
|
∂t |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
→→
Внеся E и H под знаки операторов дифференцирования, получим
→ → |
|
|
2 |
|
1 ∂(εa E2 ) |
|
1 |
∂(µa H |
2 ) |
|
|||||||||||||||
−div E×H |
= γ |
E |
|
+ |
|
|
|
|
|
|
|
|
|
+ |
|
|
|
|
|
|
|
= |
|||
|
2 |
|
|
|
∂t |
|
|
|
2 |
|
|
|
∂t |
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
= γ E2 + |
|
∂ |
|
εa E |
2 |
|
µa H |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
+ |
|
. |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
∂t |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
2 |
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
В окончательном виде теорема Умова-Пойнтинга в |
|||||||||||||||||||||||||
дифференциальной форме записи имеет вид |
|
|
|
||||||||||||||||||||||
|
→ → |
|
|
|
2 |
|
|
|
∂ |
ε |
|
E2 |
|
µ |
H 2 |
|
|
|
|||||||
−div |
E×H |
= γ |
E |
|
+ |
|
|
|
a |
|
|
|
+ |
|
a |
|
|
|
. |
(1.20) |
|||||
|
|
∂t |
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
2 |
|
|
|
|||||||
В (1.20) |
|
|
εa E2 |
|
|
представляет |
|
|
собой |
|
энергию |
||||||||||||||
|
|
|
2 |
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
H 2 |
|
|
|
|
электрического поля в единице объёма, |
µ |
|
– энергию |
||||||||||||||||||||||
|
a |
|
|
||||||||||||||||||||||
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
магнитного поля в единице объёма |
|
|
dV , |
|
γ E2 |
– тепловые |
|||||||||||||||||||
потери электромагнитной энергии. Для определения энергии во всём объёме проинтегрируем выражение (1.20) по объёму V . С учётом теоремы Остроградского-Гаусса
(1.7) теорема Умова-Пойнтинга в интегральной форме
записи примет вид
34
→ → |
→ |
|
d |
|
εa E |
2 |
|
µa H |
2 |
|
|
−∫ E×H dS |
= ∫γ E2dV + |
∫ |
|
+ |
|
dV . (1.21) |
|||||
dt |
|
|
|||||||||
S |
|
V |
V |
2 |
|
|
2 |
|
|
|
|
Векторное |
произведение |
|
→ |
|
→ |
= |
→ |
||||
E×H |
П называют |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
вектором Пойнтинга, он характеризует значение и направление перемещения энергии, проходящей в единицу
→
времени через единицу площади, перпендикулярной П
→
(рис. 1.3). Если вектор П направлен внутрь поверхности, то его поток, проходящий через поверхность, будет
→→
положительным: |
−∫ П dS >0 |
(при положительном |
|
S |
|
→ |
в сторону |
внешней нормали к |
направлении dS |
||
поверхности). |
|
|
|
dS |
|
|
dS |
E |
|
S |
|
|
H |
|
V
П=[E H]
→ → →
Рисунок 1.3 – Ориентация векторов E , H , П относительно поверхности S, ограничивающей объём V
Физический смысл теоремы Умова-Пойнтинга: энергия электромагнитного поля расходуется на тепловые потери
|
|
|
35 |
|
|
|
|
|
|
γ E2 и на |
приращение электрической |
∂ |
|
εa E |
2 |
и |
|||
|
|
||||||||
|
|||||||||
|
|
|
|
∂t |
2 |
|
|
||
магнитной |
∂ |
µa H 2 |
энергий в заданном объёме. |
|
|
||||
|
|
|
|||||||
|
∂t 2 |
|
|
|
|
|
|
||
Теорема Умова-Пойнтинга в комплексной форме записи используется для описания энергетического баланса гармонических электромагнитных колебаний. При выводе её воспользуемся первым и вторым уравнениями Максвелла (1.14):
• |
• |
• |
• |
• |
→ |
→ |
→ |
→ |
→ |
rot H = γ E+ jωεa E ; rot E = − jωµa H .
При нахождении полной мощности необходимо
•
→
комплекс вектора напряжённости электрического поля E умножить на сопряжённый комплекс вектора
•
→
напряжённости магнитного поля H и проделать последовательность операций, изложенных выше при получении уравнения для мгновенных значений. В итоге получим
• |
+ 2 jω µa H |
2 |
εa E |
2 |
|
→ |
(1.22) |
||||
−div П = γ E2 |
− |
. |
|||
|
2 |
|
2 |
|
|
Для определения энергии во всём объёме проинтегрируем (1.22) по V и применим к левой части теорему Остроградского-Гаусса (1.7):
|
• |
|
µa H |
2 |
|
εa E |
2 |
|
|
→ → |
− |
||||||
−∫ П dS = ∫γ E2dV + j2ω∫ |
|
|
dV . (1.23) |
|||||
S |
V |
V |
2 |
|
|
2 |
|
|
Первое слагаемое правой части представляет собой активную мощность, второе – реактивную мощность. Таким образом, теорема Умова-Пойнтинга в комплексной
36
форме может быть записана и сформулирована следующим образом:
•
→→
−∫ П dS = P + jQ ,
S
•
→
поток комплексного вектора Пойнтинга П сквозь замкнутую поверхность равен комплексной мощности, выделяемой внутри объёма, ограниченного этой поверхностью.
1.5 Законы изменения векторов электромагнитного поля на границе раздела двух сред (граничные условия)
Для областей, содержащих границу раздела двух или более сред, непосредственное решение дифференциальных уравнений Максвелла невозможно. Обычно решают уравнения для каждой среды в отдельности, а полученные решения «сшивают» на границе раздела. Для этого используют так называемые граничные условия – соотношения между значениями векторов поля по обе стороны от границы раздела. Поскольку уравнения Максвелла являются векторными, а решения их обычно находятся в проекциях на оси координат, то граничные условия удобно представить в виде нормальной (проекция на оси у ) и тангенциальной (проекция на оси x )
составляющих (рис. 1.4).
Методика вывода граничных условий базируется на использовании уравнений Максвелла (1.9) – (1.12).
37
y
Kn |
K |
Kτ |
|
|
x |
||
|
|
|
→
Рисунок 1.4 – Представление вектора K в проекциях на оси координат в виде тангенциальной ( Kτ ) и
нормальной ( Kn ) составляющих
Граничные условия для нормальных составляющих поля. Пусть достаточно гладкий элемент поверхности ∆S разделяет две среды 1 и 2 с различными диэлектрическими проницаемостями ε1 и ε2 ; в каждой среде параметры ε1 и
ε2 постоянны (рис. 1.5).
|
dS1 |
D , B |
1 |
|
|
|
|
α1 |
|
|
|
|
1 |
|
1 |
∆S1 |
|
|
|
ε1, µ1 |
∆S |
dS |
h |
|
|
∆ |
|
||
ε2, µ2 |
|
|
α2 |
|
|
|
|
|
|
D2, B2 |
|
|
∆S2 |
|
2 |
dS2 |
|
||
|
|
|||
Рисунок 1.5 – Преломле-
→ →
ние векторов D и B на границе раздела сред
Обозначим вектор электрической индукции (электрического смещения)
→
в среде 1 через вектор D1 , в
→
среде 2 – через вектор D2 .
Построим на плоской границе раздела элемент цилиндра с высотой ∆h →0
и основаниями |
цилиндра |
∆S1 = ∆S2 = ∆S . |
Векторы |
→→
dS1 и dS2 будут направлены перпендикулярно
38
→ → |
→ |
основаниям и поверхности раздела ( dS = n dS , где |
n – |
нормаль к поверхности раздела), на которой в общем случае распределён заряд с поверхностной плотностью
|
→ |
→ |
через |
σ = dq . Обозначим угол между векторами D |
и dS |
||
dS |
1 |
1 |
|
|
|
|
|
→ |
→ |
|
|
α1 , между векторами D2 |
и dS2 – через α2 . |
|
|
Воспользуемся третьим уравнением Максвелла в интегральной форме (1.11):
→ →
∫ D dS = ∫ρdV = ∫σ dS .
S V S
Суммарный интеграл по поверхности цилиндра в левой части будет содержать два интеграла по основаниям и один – по боковой поверхности, который можно исключить с учётом того, что при ∆h →0 площадь боковой поверхности ∆Sбок →0. Тогда для выделенного цилиндра
→ → |
∫ |
→ → |
+ ∫ |
→ → |
= ∫σ dS . |
∫ D dS = |
D1 dS1 |
D2 dS2 |
|||
S |
∆S1 |
|
∆S2 |
|
∆S |
Поскольку ∆S1 = ∆S2 = ∆S , а dS1 = dS2 = dS , то с учётом
равенства подынтегральных скалярных произведений при переходе к записи в проекциях на оси координат получим
D1 cosα1 + D2 cos(180°−α2 )=σ .
Для окончательной записи граничных условий разложим,
→
например, вектор D1 на
нормальную и тангенциальную составляющие (рис. 1.6).
D1τ
D1n α1 D1
Рисунок 1.6 – Пример разложения вектора
→
D на составляющие
39
Из рис. 1.6 видно, что
cosα1 = D1n .
D1
Аналогичное разложение можно провести и для
→
вектора D2 . Тогда
D1 cosα1 = D1n , D2 cos(180°−α2 )= −D2n .
Учитывая, что среды 1 и 2 изотропны, а заряд σ распределён по поверхности раздела равномерно,
получаем окончательное выражение для нормальных составляющих вектора электрической индукции на границе раздела сред:
D1n − D2n =σ . |
(1.24) |
Физический смысл: нормальная составляющая вектора электрической индукции Dn при переходе через границу
раздела двух сред претерпевает скачок, численно равный поверхностной плотности электрического заряда σ .
Граничное условие для нормальной составляющей магнитного поля можно получить из четвёртого уравнения
→ →
Максвелла (1.12): ∫ B dS = 0 .
S
Вэтом случае, вывод граничных условий аналогичен выводу, приведённому выше, для электрического поля (студенты вывод выполняют самостоятельно).
Врезультате получаем
B1n − B2n = 0
или
B1n = B2n . |
(1.25) |
Физический смысл: нормальная составляющая вектора магнитной индукции Bn при переходе через границу
раздела двух сред не изменяется.
40
Граничные условия для тангенциальных составляющих поля выводятся из первого и второго уравнений Максвелла (1.9), (1.10).
Пусть достаточно гладкая поверхность S разделяет две среды 1 и 2 с различными магнитными (диэлектрическими) проницаемостями µ1 и µ2 (ε1 и ε2 );
параметры сред постоянны (рис. 1.7). Обозначим вектором
→ |
→ |
H1 |
напряжённость магнитного поля в среде 1, H2 – в |
среде 2. Охватим границу небольшим контуром длиной и высотой ∆l , состоящим из элементарных отрезков ∆l1 ,
∆l2 и 2 ∆l3 , определяемых по направлению единичными
→ |
→ |
→ |
|
|
векторами dl1 , |
dl2 |
и dl3 . Предположим, что |
∆l3 |
= ∆h →0, |
а по поверхности S в бесконечно тонком слое, помещённом на границе раздела, протекает поверхностный ток с плотностью iпов = dI / dl . Обозначим
|
|
|
|
|
|
|
|
|
|
→ |
|
|
|
|
|
|
|
|
|
|
|
|||
угол падения вектора H1 на границу через α1 , |
а угол |
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
→ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
преломления вектора H2 |
|
– через α2 . |
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
H1, E1 |
|
|
|
|||||
|
|
|
|
|
|
∆ l1 |
|
(dl1) |
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
∆ l3 |
|
µ1, ε1 1 |
|
|||||||||||||
∆ h |
|
|
|
|
|
∆ l |
|
α1 |
|
|
|
|||||||||||||
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
α2 |
|
|
|
|
|
|
|
µ2, ε2 |
2 |
|
||||||||
|
|
|
|
|
|
|
|
(dl3) |
|
|||||||||||||||
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
∆ l2 |
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
(dl2) |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
H |
2, |
E |
|
|
|
|
|
|
|
|
|
|
|
|
|
→ |
→ |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Рисунок 1.7 – Преломление векторов H и E на границе раздела сред
