Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Кураев А.А., Попкова Т.Л., Синицын А.К. Электродинамика и распространение радиоволн

.pdf
Скачиваний:
340
Добавлен:
02.05.2014
Размер:
5 Мб
Скачать
☆

r

r

 

r

 

r

− E rot H

 

 

 

 

div [E, H ]= H rot E

 

 

 

 

r r

 

 

r r

 

 

↓r r

r ↓r

=

 

div[E, H

]= [E, H ]= E, H

+ E,H

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

r

r

 

 

r

r

r

r

r

= H

,

E

− E

,

H = Hrot

E −

Erot

H .

Здесь - вектор-оператор типа

= Xr0 ∂∂x +Yr0 ∂∂y + Zr0 ∂∂z ,

стрелки сверху указывают вектор, на который действует .

Составимr баланс энергии. Для этого первое УМ умножим на ( − E ),

второе на H и сложим, используя затем доказанное векторное тождество. При этом получим

r r r r

r r

r

∂D

r

∂B

 

r r

 

Hrot E − Erot H = div[E, H ]= − E

+ H

 

−δ E .

(3.4)

∂t

∂t

 

 

 

 

 

 

 

Перейдем к интегральной формулировке полученного уравнения, интегрируя левую и правую части по объему V с граничной поверхностью S и используя в левой части теорему Остроградского-Гаусса:

r r

r

r

r

 

r

∂B

 

 

r r

=

∫[E, H ]dS

= −∫ E

∂D

+ H

dV − ∫δ EdV

s

 

V

 

∂t

 

 

∂t

 

V

 

 

 

 

r

 

r

 

r

r

 

 

− PП6 .

(3.5)

 

 

= −∫ E

∂D

+ H ∂B

dV + Pcт

 

 

V

 

∂t

 

 

∂t

 

 

 

 

В записанном уравнении нам известен смысл только последних двух членов справа – это мощность сторонних источников и омических потерь в объеме V. Для того чтобы установить физический смысл остальных членов, рассмотрим два следующих предельных случая.

Первый случай.

Пусть имеется объем V, в котором

1) Рсm= 0,

21

2)V изолирован от rвнешнегоr пространства непрозрачной S, так что на внешней стороне S ∫[E, H ]nrdS = 0 .

s

Тогда имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

PП6 = −

 

 

 

 

r

r

 

r

∂B

 

 

 

 

 

 

 

 

 

 

 

 

 

E

∂D

+ H

dV .

 

 

 

 

 

 

 

 

 

 

 

 

 

V∫

 

 

 

∂t

 

 

 

∂t

 

 

 

 

 

 

 

 

 

 

Но мощность потерь в изолированном объеме V должна быть численно

равна скорости убывания запасенной энергии − dW

. Следовательно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

r

r

 

 

 

r

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

∂D + H

∂B dV = ∂W = ∂wdV ,

 

 

 

 

 

 

 

 

 

∂t

 

 

 

 

 

 

∂t

 

 

 

 

 

 

dt

 

V∫ ∂t

 

 

 

 

 

V∫

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где w – удельная энергия ЭМП.

 

 

 

 

 

Таким образом,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂w

 

 

 

 

∂D

 

 

 

r

 

∂w

 

∂w

 

 

 

 

 

 

 

 

 

 

r

 

 

r

∂B

 

 

м

 

 

 

 

 

 

 

 

 

= E

 

 

 

 

 

+ H

 

 

 

=

 

 

Э

+

 

.

 

 

 

 

(3.6)

 

∂t

 

∂t

 

∂t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂t

 

∂t

 

 

 

 

 

Положим, что εa ≠ f (t) и µa ≠ f (t) (диэлектрических и магнитных

потерь нет). Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂wЭ

 

 

r

 

∂

 

 

r

 

 

 

r

 

 

r

εa ∂

 

r r

∂

 

 

2

 

 

= E

 

εa E

= ε

a

E

 

∂E

=

(EE)=

 

εa E

 

.

 

 

 

 

 

 

 

 

 

 

∂t

 

 

 

 

 

 

 

∂t

 

 

 

∂t

2∂t

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂t

 

 

Таким образом, удельная энергия электрического поля из УМ определяется как

wЭ =

ε

а

Е2

.

(3.7)

 

 

2

 

 

 

 

 

Аналогичным путем находим формулу удельной энергии магнитного

поля:

wМ =

µ

а

Н2

.

(3.8)

 

 

2

 

 

 

 

 

Второй случай:

1) Pст ≠ 0 ;

22

2)процесс установился, т.е. dWd t = 0 ;

3)PП6 = 0;

4)S прозрачна для ЭМП.

 

В этом случае имеем

 

 

 

s∫

 

r

r

r

 

 

 

 

[E, H ]dS = Pcm .

 

(3.9)

 

 

 

 

С другой стороны, поскольку в рассматриваемом случае в объеме V нет

потерь и нет накопления электромагнитной энергии, Pcm = PΣ , где PΣ -

мощность излучения через граничную поверхность S. Следовательно,

 

r

 

v

 

r

 

r r

S0 - плотность потока энергии через

 

∫ E,

H

dS

=

P∑ = ∫S0dS где

S. Та-

s

 

 

 

 

r

r S r

 

 

ким образом,

S0

= [E, H ]. Эта величина называется вектором Умова-

 

Пойнтинга. Этот вектор указывает величину плотности потока энергии и направление распространения энергии ЭМП.

Возвращаясь к исходному уравнению баланса энергии ЭМП, имеем

s∫

r

r

 

 

dW

 

 

S

0 dS

+ PП6

+

dt

= Pcт.

(3.10)

Этот закон сохранения энергии для ЭМП носит название теоремы Умо- ва-Пойнтинга.

В дифференциальной форме он выражается в виде

div Sr0 + pП6 +

∂w

 

= pст .

(3.11)

∂t

 

 

 

 

 

 

 

 

 

 

Получим теперь выражение для диэлектрических и магнитных потерь,

что соответствует случаю εa (t) и

µa (t) . Будем иметь в виду периодические

электромагнитные процессы с периодом Т. Рассмотрим случай, когда

s∫Sr0dSr = 0, PП6 = 0 , среднее за период значение

 

T постоянно, поэтому сред-

W

нее значение

 

T

 

 

 

 

 

 

dW

= 0.

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

Таким образом имеем

 

 

 

r

 

r

 

r

r

 

 

 

 

pcт = E

∂D

+ H

∂B

.

 

 

 

 

∂t

 

 

 

∂t

 

 

 

 

 

 

 

 

 

 

 

 

 

23

Найдем среднее за период Т значение величин:

 

1

Т

 

 

 

 

 

1

 

T

r

 

 

 

 

 

 

r

 

∂B

 

 

 

 

 

 

∫pст dt = pcmT =

 

∫

E

∂D

dt + H

 

 

dt =

 

 

 

T

 

T

 

 

 

 

 

 

∂ t

 

 

 

 

0

 

 

 

 

 

0

 

 

 

∂

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

1

∫

 

r

r

 

1

 

 

∫

 

 

 

r r

 

 

T

 

 

 

Т

 

 

 

 

=

 

 

EdD +

 

 

 

 

 

 

 

 

HdB = pПε

+ pПµ ,

(3.12)

 

 

 

T

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lэ

 

 

 

 

 

 

 

 

 

 

 

 

lM

 

 

 

 

 

 

 

 

 

 

 

 

где lЭ - контур диэлектрического гистерезиса, lМ

- магнитного.

Таким образом,

 

pТ

 

=

1

 

S

 

,

 

рТ

=

1

 

S

 

, т.е. среднее за период Т удель-

 

 

 

 

 

 

Т

 

 

 

 

 

 

 

 

 

Пε

 

Т

 

 

Э

 

 

 

Пµ

 

 

 

 

М

 

 

 

ные диэлектрические потери

 

ТПε

пропорциональны «площади» петли

р

диэлектрического гистерезиса Sэ = ∫EdD , магнитные -

 

ТПµ соответственно

p

 

 

 

r

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lэ

 

 

 

 

 

 

 

 

 

 

площади SМ = ∫HdB петли магнитного гистерезиса.

 

 

lМ

24

ГЛАВА IV

КОМПЛЕКСНЫЕ АМПЛИТУДЫ И ТЕОРЕМА

ОКОМПЛЕКСНОЙ МОЩНОСТИ

4.1.УМ в комплексной форме

Вбольшинстве случаев электродинамические задачи решаются, как и в радиотехнике, в спектральной области и решение ищется для гармонических процессов. Предполагается, что ЭМП имеет монохроматический характер, т.е. частота колебаний векторов ЭМП постоянна. В этом случае используется метод комплексных амплитуд (МКА). При этом имеется в виду

ψ =ψm cos(ω t +ϕ)= Reψme

j(ω t+ϕ )

&

jω t

&

(4.1)

 

= Reψme

 

= Reψ .

Величина ψ&m =ψme jϕ называется комплексной амплитудой функции ψ . Символ Re обозначает выделение действительной части комплексного числа ψ& .

Аналогично можно записать в комплексной форме и векторную вели-

чину

Vr =Vr

cos(ωt +ϕ)= Re Vr&

e jω t , Vr&

=Vr

e jϕ .

(4.2)

m

m

m

m

 

 

Для дальнейшего анализа важны следующие свойства комплексных представлений гармонических процессов.

1.Дифференцирование Vr&(t) по времени t равносильно умножению Vr& на jω :

r

∂t = jω Vr& .

2.Интегрирование Vr&(t) по t равносильно делению Vr на jω :&

∫Vr&dt = Vr& .

jω

3. Справедливо следующее соотношение:

(4.3)

(4.4)

25

Vr1(ω t)Vr2 (ω t)T = 12 Re(Vr&1Vr&2* )= 12 Re(Vr&1*Vr&2 ),

(4.5)

 

 

 

T =

2π

, знак

 

T означает усреднение по t в интервале t [0,T ].

 

ω

 

 

 

 

В дальнейшем будем пользоваться комплексными векторами

Er&, Dr&, Hr& , Br& .

Обратимся к вычислению диэлектрических и магнитных потерь в комплексной форме.

Диэлектрические потери

В этом случае вектор D запаздывает относительно E на угол диэлектрических потерь ∆Е:

 

 

r&

 

r

 

− j∆E

 

 

 

ε

=

D

=

Dm

e

=ε'

− jε" .

(4.6)

r

 

 

r

 

&

 

E&

 

 

 

 

 

 

1a

 

 

Em

 

 

1a

1a

 

Индекс «1» здесь поставлен затем, чтобы выделить диэлектрические потери из общих электрических потерь, куда входят и джоулевы.

Магнитные потери

r

В этом случае вектор B запаздывает относительно H на угол магнитных потерь ∆Н:

 

r

r

 

 

 

 

 

B&

Bm

e

− j∆H

= µ'a − jµ"a .

 

µa =

r

=

 

 

(4.7)

r

 

&

H&

 

 

 

 

 

Hm

 

 

 

 

Таким образом, в комплексной форме УМ сводятся к первым двум и записываются в следующей форме:

 

r&

 

 

 

r&

r&

 

 

 

 

rot H = jωε

 

E +

δ ,

 

 

 

 

 

r

 

&

r

 

 

 

 

(4.8)

 

 

1a

 

 

 

 

 

&

 

 

&

&

 

 

 

 

 

rot E

= − jω µa H.

 

 

 

 

 

 

 

 

 

 

 

 

 

r

r

r

Преобразуем первое УМ с учетом того, что δ&

=σ E&

+δ&ст:

 

r&

 

&

 

σ

r&

r&

 

 

rot H = jω

ε1a − j

E +δст.

 

 

 

 

 

 

 

 

ω

 

 

 

 

26

Обозначим

 

&

 

− j

σ

 

=

&

 

= ε '

 

− jε

"

 

= ε'

 

ε"

+

σ

 

ε

 

 

ω

 

ε

 

 

 

− j

,

 

 

 

1a

 

 

 

 

 

a

 

a

 

σ .

 

a

1a

 

1a

 

ω

ε'

 

 

= ε'

 

 

, ε"

 

= ε"

+

 

 

 

 

 

 

 

 

 

a

 

 

1a

 

 

 

a

 

1a

 

 

ω

 

 

 

 

 

 

 

С учетом этого обозначения УМ в комплексной форме приобретают

вид

 

r&

 

r&

r&

,

 

rot H = jωε

E +δ

 

 

r

&

a

r cm

 

(4.9)

 

&

 

&

&

 

 

rot E

= − jω µa H.

 

 

4.2. Теорема о комплексной мощности

Составим баланс энергии в комплексной форме. Для этого умножим

первое комплексно сопряженное УМ на (− Er&), второе – на Hr& * и сложим. В результате получим

div[Er&,Hr&*]=−jω{µ&a Hm2 −ε&a* Em2 }− Er&δr&cm* =

= − jω{µ'a Hm2 −ε'a Em2 }−ω{µ"a Hm2 + ε"a Em2 }−δr&cm* Er& . (4.10)

Определим среднее за период Т величины:

 

r

T

 

 

 

r

r

 

 

r

r

 

 

r

r

 

 

 

 

S0

 

 

= 12 Re[E&, H& * ]= Re S&0

, S&0

=

12 [E&

, H& * ],

 

 

 

 

 

 

 

 

 

 

 

r&

r&*

′

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&

 

 

 

 

 

 

 

 

 

 

 

 

 

w

Т

 

=

 

1

Re

εa E E

=

 

εa Em

 

,

 

 

 

 

 

 

э

 

2

 

4

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4.11)

 

 

 

 

 

 

 

 

 

r&

r&

 

 

 

 

 

 

 

 

 

 

 

 

wT = 1

 

 

 

 

 

2

 

 

 

 

 

Re

µa H H

 

=

µ'a Hm ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

M

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

r

 

 

 

 

 

 

 

 

 

 

 

r& r&

*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

pcmT =

1

Re (−δ&*E&)= Re

P&cm

, pcm

= −

Eδ

 

= pa + jpr ,

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

pПТ =

1

ω {µ"a Hm2 +ε"a Em2 } .

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

27

Учитывая перечисленное, баланс энергии в дифференциальной форме

можно записать в виде

div Sr&0 + рПТ + 2 jω{wMT − wэТ}= pcm .

В интегральной форме баланс энергии имеет вид s∫Sr&0dSr + PT + 2 jω {W TM −W Tэ}= Pcm .

Разделяя действительную и мнимую части уравнения, получаем теорему об активной мощности :

Re s∫Sr&0dSr + PTП = Pa ,

и теорему о реактивной мощности :

Jms∫Sr&0dSr + 2ω {W TM −W Tэ}= Pr .

Теорему об активной мощности можно более подробно записать так :

Re

Sr&dSr +

ω

{µ"

 

H 2 +ε"

E 2 }dV = P .

(4.12)

s∫

 

2 V∫

 

a

m

a m

a

 

 

 

 

 

 

Смысл теоремы об активной мощности : активная мощность Ра , развиваемая в объеме V, расходуется на потери в V и активную мощность излучения через граничную поверхность S.

Для выяснения смысла теоремы о реактивной мощности рассмотрим случай, когда поток реактивной мощности через граничную поверхность S отсутствует. Тогда

2ω {W TM −W Tэ}= Pr

Если частота источника ω соответствует одной из резонансных частот объема V −ωoi , W TM =W Тэ и реактивная мощность Рr=0.

При ω < ωoi W TM > W Tэ и нужен сторонний накопитель электрической

энергии, которым становится источник питания данного объема (причем обмен энергии происходит дважды за период). Таким образом, реактивная мощность Рr количественно описывает колебания энергии между источником и нагрузкой при ω ≠ωoi .

28

ГЛАВА V

ВОЛНОВЫЕ УРАВНЕНИЯ И ЭЛЕКТРОДИНАМИЧЕСКИЕ ПОТЕНЦИАЛЫ

5.1. Волновые уравнения

Запишем четыре УМ в действительной форме, полагая εa , µa , σ =Const . Последнее означает, что среда однородна, а диэлектриче-

ских и магнитных потерь нет, есть только джоулевы потери. Систему отсчета будем считать неподвижной.

rot H = ∂∂Dt +δr,

rot Er = −∂∂Bt ,

r

div D = ρ ,

r div B = 0.

Применим к обеим частям уравнения (5.2) операцию rot:

rot rot Er = −∂∂t rot Br = −µa ∂∂t rot Hr .

Используя (5.1), тождество rot rot E = grad div E − 2 E , а также чаем

r

∂

2

r

 

1

 

∂δ

 

2 E − µa εa

 

E

=

grad ρ + µa

.

∂t 2

 

 

 

 

εa

∂t

Учтем проводимость среды : δ =σ E +δcm . Тогда имеем

r

∂

2

r

 

∂E

 

1

 

∂δcm .

2 E − µaεa

 

E

− µaσ

=

grad ρ + µa

∂t 2

∂t

 

 

 

 

εa

∂t

(5.1)

(5.2)

(5.3)

(5.4)

(5.3), полу-

(5.5)

Векторноеr уравнение Д’Аламбера (5.5) является волновым уравнением для вектора E ЭМП. Заметим, что в правой части уравнения (5.5) стоят толь-

29

ко функции источников поля, поскольку токи проводимости сами по себе не создают ρпров . Докажем это.

Всоответствии с уравнением непрерывности

∂ρ∂провt = −divδrпров = −σ divEr = −εσa ρпров .

Впоследнем приравнивании учтено (5.3) при условии ρстор = 0 . Таким образом, для ρпров имеем следующее уравнение:

 

∂ρпров

+

σ

ρпров = 0 ,

(5.6)

 

∂ t

εa

 

 

 

 

 

 

 

 

 

−

σ

t

 

 

 

 

 

которое имеет такое решение: ρпров = ρ0e

εa , где ρ0 - начальное возмуще-

ние ρпров . Таким образом, в установившемся режиме (t → ∞) ρпров = 0. Поэтому ρпров ≠ 0 только в присутствии сторонних источников, тогда в уравне-

нии (5.6) появляется ненулевая величина − σ ρ стор.

εа

Получим волновое уравнение для вектора Н . С этой целью применим к обеим частям уравнения (5.1) операцию rot:

r

∂

r

r

rot rot H =

 

rot D + rotδ .

∂t

 

 

 

Далее воспользуемся уравнениями (5.2), (5.4), а также тождествомr

rot rot Hr = grad div Hr − 2 H . При этом получим, учитывая, что δ =σ E +δcm

r

∂

2

r

 

∂H

r

 

2 H −εa µa

 

H

−σ µa

= −rotδcm .

(5.7)

∂t 2

∂t

 

 

 

 

Левые части волновых уравнений для вектора E (5.5) и вектора Н (5.7) идентичны. Это свидетельствует о том, что в пространстве, свободном от ис-

точников (т.е. когда правые части (5.5) и (5.6) равны нулю) вектора E и Н имеют одинаковый волновой характер. Как следует из теории волновых

уравнений, коэффициент при вторых производных по времени равен Vф−2 (Vф

- фазовая скорость волнового процесса в среде). Таким образом, фазовая скорость волнового электромагнитного процесса в свободном пространстве

30