
ТОЭ, ТЭМП, Шмелёв, Сбитнев, 2003 г. / Глава5
.pdfГлава 5. ПЕРЕМЕННОЕ ГАРМОНИЧЕСКОЕ ЭЛЕКТРОМАГНИТНОЕ ПОЛЕ
______________________________________________________________________
______________________________________________
§ 5.1. Основные уравнения электромагнитного поля в комплексной форме. Уравнения Максвелла в комплексной форме
Электромагнитное поле, в котором токи, заряды, потенциалы и составляющие векторов меняются по гармоническому закону с одной и той же заданной часто-
той, называется гармоническим электромагнитным полем. Для анализа таких по-
лей целесообразно использовать метод комплексных амплитуд, который применяется в теории цепей синусоидального тока. Пусть в декартовой системе координат задан некоторый вектор N(Q,t), составляющие которого меняются по гармоническому закону с одинаковой частотой
N(Q,t) = 1xNxm(Q) sin(ωt+α) + 1yNym(Q) sin(ωt+β) +1zNzm(Q) sin(ωt+γ) (1)
Амплитуды и начальные фазы составляющих могут быть функциями пространственных координат, но не зависят от времени. При совпадении фаз всех трех составляющих вектор будет меняться по закону синуса, не меняя направления в пространстве. В этом случае говорят, что вектор N(Q,t) линейно поляризован. В общем случае, когда α≠β≠γ, вектор будет вращаться в пространстве, опи-
76

сывая при этом эллипс. В этом случае говорят, что вектор N(Q,t) эллиптически поляризован.
Комплексной амплитудой вектора N(Q,t) будем называть выражение
N&m =1x N xm e jα+1y N ym e jβ+1z Nzm e jγ =1x N&xm +1y N& ym +1z N&zm
С учетом введенного обозначения выражение (1) можно записать в виде
& |
jωt |
|
- jωt |
|
|
|
− N me |
|
) |
(2) |
|
N(Q,t) = 1/(2j)( Nm e |
|
|
Комплексным действующим значением вектора N будем называть выражение
N& = N&2m
С учетом этого обозначения выражение (2) для мгновенного значения вектора можно записать в виде
1 |
& |
jωt |
|
- jωt |
|
N(Q,t) = j 2 |
( N e |
|
− Ne |
|
) |
Векторы N&m (Q) и N& (Q) являются комплексными представлениями гармони-
чески изменяющегося вектора N(Q,t). Используя свойства комплексных представлений, можно из уравнений Максвелла в пространственно-временной форме получить уравнение Максвелла в комплексной (пространственно-частотной) форме. Запишем эти уравнения для действующих значений векторов электромагнитного
поля |
|
|
|
|
|
rot H& |
= δ&п, |
где δ&п= δ& + δ&пр + δ&см , |
δ&см = jωD& |
||
|
|
|
|
|
(3) |
& |
& |
& |
& |
& |
& |
rot( E − Ec ) = – jωB , |
div B = 0, |
div D =ρ |
|||
Система уравнений (3) дополняется уравнениями материальной связи, кото- |
|||||
рые в линеаризованном виде записываются следующим образом |
|||||
|
B& = µa H& , D& = εa E& , δ&пр = γE& |
(4) |
77

Теорема Умова-Пойнтинга в комплексной форме
Объемная плотность комплексной мощности, потребляемой материальной
точкой в гармоническом ЭМП, равна |
|
|
|
|
|
|||
~ |
& & |
|
& |
|
& & |
|
|
& & |
s |
= ( E − Ec ) δп+ jωB · H = ( E − Ec )rot H |
– H rot( E − Ec ) = |
= – div(( E& − E&c )×H )
Комплексная мощность, потребляемая внутри объема V, равна
|
– ∫ |
|
|
|
|
|
div(( E& − E&c )×H )dV = – ∫ (( E& − E&c )×H )dS |
||||
|
V |
|
|
|
S |
|
Используя соотношения (3) и (4), можно доказать, что |
||||
– ∫ (( E& − E&c )×H )dS = ∫ γE2dV – ∫ ( E&с δп – E& δ) dV + ∫ ( E& (jωD& )+ H jωB& )dV= |
|||||
S |
|
V |
V |
V |
|
= ∫ |
~ |
|
2 |
2 |
~ |
pтdV – ∫ sист dV + jω ∫ |
(µaHm |
/2 + εaEm |
/2)dV = Pт - Sист + 2jω(Wм ср – Wэ ср) = |
||
V |
V |
V |
|
|
|
= (Pт – Pист) + j(Qэм – Qист) |
|
|
(5) |
||
|
(5) – это уравнение баланса активных и реактивных мощностей для объема V. |
Уравнение (5) иначе называют теоремой Умова-Пойнтинга в комплексной форме. Левая часть уравнения (5), содержащая поверхностный интеграл, равна комплексной электромагнитной мощности, поглощаемой объемом V из окружающего пространства.
Иначе уравнение (5) можно записать так:
~ |
~ |
, |
Sист |
=Pист + jQист = Pт + jQэм + Sизл |
~
Sизл – комплексная мощность, излучаемая объемом V в окружающее пространство.
Комплексная форма теоремы Умова-Пойнтинга имеет важное практическое значение. Ее используют при расчете электромагнитных излучателей и направ-
78

ляющих систем в радиоэлектронной аппаратуре. Пользуясь уравнением (5), можно определить внутреннее активное и реактивное сопротивление проводника, если в результате анализа ЭМП известно поверхностное распределение комплексных действующих значений напряженностей электрического и магнитного поля.
– ∫ |
( E& |
|
×H )dS = (r + jx)I2 |
||
S |
|
|
Теорема о единственности
Пусть в некотором объеме V, ограниченном замкнутой поверхностью S, известно распределение параметров электрофизических свойств среды и распределение сторонних источников ЭМП. Пусть в этом объеме установился синусоидальный режим ЭМП, в создании которого могли участвовать источники поля, расположенные вне объема V. Пусть, кроме того, известны комплексные значения тангенциальных составляющих вектора E& или H& на граничной поверхности S.
Тогда объемное распределение векторов E& и H& , удовлетворяющее уравнениям Максвелла и граничным условиям, является единственным решением задачи анализа ЭМП.
Доказательство:
Допустим, что существует два различных решения ( E&′, H& ′) и ( E&′′, H& ′′), удовлетворяющих уравнениям Максвелла и граничным условиям. Ввиду линейности
уравнений поля разность этих решений E& = E&′– E&′′и |
H& = H& ′– H& ′′также удовле- |
|||
творяет уравнениям Максвелла при следующих дополнительных условиях: |
||||
а) E&c = 0 и δ& = 0; |
|
|
|
|
б) во всех точках поверхности S или E& |
= 0 или |
H& |
t |
= 0 т.к. по предположению |
t |
|
|
|
E&t′= E&t′′. К разностному полю ( E& , H& ) применим теорему Умова-Пойнтинга
79

– ∫ |
( E& |
|
γ| E& |2dV + jω ∫ (µa| H& |2 – εa| E& |2)dV |
|
×H )dS = ∫ |
(6) |
|||
S |
|
V |
V |
|
Поверхностный интеграл в левой части равен нулю в соответствии с услови-
ем (б), значит соотношение (6) может выполняться только при E& = 0 и H& = 0 во всех точках объема V, из этого следует, что оба решения ( E&′, H& ′) и ( E&′′, H& ′′) тождественны.
Контрольные вопросы
1.Что называется гармоническим электромагнитным полем?
2.Что называют линейно и эллиптически поляризованным вектором?
3.Что называют комплексным амплитудным и комплексным действующим значением гармонически изменяющегося вектора?
4.Как связано мгновенное значение вектора с комплексным амплитудным и комплексным действующим?
5.Какой вид имеют уравнения Максвелла в комплексной форме?
6.Как записываются уравнения материальной связи в комплексной форме?
7.Чему равна объемная плотность комплексной мощности, потребляемой материальной точкой в гармоническом ЭМП?
8.Чему равна комплексная мощность, поглощаемая некоторым объёмом из окружающего пространства?
9.Как записывается уравнение баланса комплексных электромагнитных мощностей для некоторого объёма?
10.Как, пользуясь теоремой Умова–Пойнтинга, можно определить внутреннее активное и реактивное сопротивление проводника?
11.Как формулируется теорема о единственности гармонического электромагнитного поля?
80
12. Как доказывается теорема о единственности гармонического электромагнитного поля с помощью теоремы Умова–Пойнтинга?
§ 5.2. Уравнения математической физики относительно потенциалов гармонического электромагнитного поля
Комплексные параметры электрофизических свойств среды
Изменение электрической поляризованности или намагниченности вещества по гармоническому (в общем случае по периодическому) закону обычно сопровождается тепловыми потерями энергии. В этом случае составляющие вектора электрического смещения отстают по фазе от соответствующих составляющих вектора напряженности электрического поля. При периодическом перемагничивании ферромагнетиков составляющие вектора магнитной индукции отстают по фазе от соответствующих составляющих вектора напряженности магнитного поля. Указан-
ные сдвиги фаз между парами векторов E& и D& , H& и B& можно учесть путем введения комплексных параметров электрофизических свойств среды в уравнения материальной связи
B& = µ&a H& , D& = ε&a E&
где µ&a =µ′a – jµ′a′ = µa exp(–j∆м), ε&a = ε′a – jε′a′ = εa exp(–j∆e).
Аргументы комплексной магнитной и диэлектрической проницаемости ∆м и ∆e называются углами магнитных и диэлектрических потерь. В справочной литературе для различных электротехнических материалов даются значения тангенсов угла этих потерь
tg(∆м) = µ′a′ /µ′a ,
81
tg(∆e) = ε′a′ / ε′a .
При постановке и решении задач анализа гармонических ЭМП возможно объединение токов проводимости и токов смещения
δ&пр + δ&см= δ&и,
где δ&и – индуцированная плотность тока.
С учетом введенного обозначения закон полного тока можно записать в виде rot H& = δ&п= δ& + δ&и
δ&и= jωε&a E& +γE& = (γ+jωε&a ) E& = γ&E& ,
где γ& – комплексная удельная проводимость вещества на данной частоте ω.
Во многих случаях вместо индуцированной плотности тока в расчеты вводят эффективный вектор электрического смещения
D&эф= δ&и/(jω) = ( ε&a +γ/(jω)) E& = ε&a эф E& .
Последнее соотношение для краткости записывают как D& = ε&a E& , т.е. немно-
го изменяют систему обозначений векторов гармонического ЭМП. В дальнейшем уравнения математической физики для гармонического ЭМП будем записывать с учетом последнего обозначения.
Системы электродинамических потенциалов и уравнения математической физики для гармонического электромагнитного поля
Для вектора магнитной индукции всегда выполняется условие div B& = 0, поэтому
B& = rot A& ,
где A& – комплексный векторный магнитный потенциал.
rot( E& − E&c ) = – jωB& = – jω rot A& = – rot(jωA& ), поэтому E& − E&c = – jωA& – grad ϕ& ,
82

где ϕ& – комплексный скалярный электрический потенциал,
A& и ϕ& – система электродинамических потенциалов.
Во многих случаях корректное задание поля вектора E&c в качестве объемно распределенного источника ЭМП вызывает значительные затруднения при поста-
новке задачи анализа поля. В этих случаях вместо E&c задают векторное поле
– rot E&c = δ&м , которое называют полем сторонней плотности магнитного тока. С
учетом этого обозначения закон электромагнитной индукции можно записать в виде:
rot E& = −δ&м − jωB& div B& = − j1ωdivδ&м
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
& |
|
|
& |
|
|
|
& |
|
|
1 |
& |
|
|
& |
|
|
||
|
|
E = − jωA |
−grad |
ϕ B = − |
|
|
|
|
δм + rot A |
|
|||||||||||
|
|
|
|
|
|
|
|
|
& |
|
|
|
jω |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Теперь можно получить систему уравнений математической физики относи- |
|||||||||||||||||||||
тельно потенциалов. За основу можно взять закон полного тока |
|
||||||||||||||||||||
|
& |
|
& |
|
|
& |
|
|
& |
& |
& |
& |
|
|
|
|
& |
& |
& |
|
|
|
rot H |
= δ + jωD ; |
|
H = νа |
B = |
νаrot A – |
νа δм/(jω) |
|
|||||||||||||
|
|
|
|
& |
& |
& |
|
|
& |
& |
& |
|
|
|
|
& |
|
|
|
||
|
|
|
|
D = |
εa E = – jωεa |
A – |
|
εa gradϕ |
|
|
|
||||||||||
|
|
& |
|
2 |
|
& |
|
|
|
|
& |
|
|
1 |
|
|
& |
|
|
||
rot(νa rot |
A) |
− ω |
|
εa |
A |
+ jωεagrad ϕ=δ |
+ |
|
|
|
|
rot(νaδм) |
|
||||||||
|
& |
|
|
|
& |
|
|
& |
|
& |
|
|
|
jω |
& |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(1) |
||||||
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
& |
|
|
& |
|
|
|
|
|
|
|
|
|
|
|
|
div(εagrad ϕ+ jωεa A) = |
|
divδ |
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
& |
& |
|
|
& |
|
|
jω |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(1) – комплексная форма системы уравнений математической физики относительно векторного магнитного и скалярного электрического потенциалов. Не трудно заметить, что уравнения в системе (1) линейно зависимые (второе можно получить из первого взятием дивергенции от обеих частей и делением их на jω). Поэтому для обеспечения единственности решения системы уравнений (1), кроме гранич-
83

ных условий, нужно вводить условие калибровки электродинамических потенциалов.
Рассмотрим систему (1) для случая однородной по электрофизическим свойствам среды внутри расчетной области. Тогда скалярные поля параметров ε&a и ν&а
можно вынести за знак дифференциальных операторов и умножить обе части первого уравнения на µ&а.
|
& |
− ω |
2 & & |
& |
& |
& |
& |
& & |
1 |
|
& |
|
|
|
|
|
|
|
|||||||||
rot(rot A) |
µa εa A + jωεa µa gradϕ=µa δ+ |
jω |
rot δм |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
(2) |
||
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
& |
& |
|
|
|
|
|
|
|
|
div(gradϕ) + jωdiv A = |
& |
divδ |
|
|
|
|
|
|
|
||||
|
& |
|
|
|
jωεa |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Если к системе (2) применить условие калибровки Лоренца div A& = – jωε&aµ&a ϕ& ,
то из (2) можно получить два независимых уравнения:
2 & |
ω2 |
& |
|
& |
1 |
& |
|
||
A+ |
|
|
|
& |
|
|
|
rot δм, |
(3) |
& |
2 |
A = – µaδ – |
jω |
||||||
|
|
|
|
|
|
||||
|
vф |
|
|
|
|
|
|
|
|
|
|
2 |
|
ω2 |
|
1 |
|
& |
|
|
|
|
|
ϕ – |
|
2 |
ϕ = |
|
|
divδ, |
(4) |
|
|
|
|
|
& |
|
|
& |
jωε |
|
|
|
|
|
|
|
|
vф |
|
& |
a |
|
||
|
|
|
|
|
& |
|
|
|
|
|
|
& |
= |
1 |
– фазовая скорость электромагнитной волны. |
|
|||||||
где vф |
ε&aµ&a |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
Уравнение (3) называется векторным уравнением Даламбера, уравнение (4) - скалярным уравнением Даламбера. Если источники ЭМП отсутствуют в расчетной области, то правая часть этих уравнений равна нулю:
|
2 |
& |
ω2 |
|
|
|
|
|
& |
(5) |
|
|
A+ |
|
2 A = 0 |
||
|
|
|
vф |
|
|
|
|
|
& |
|
|
84
|
|
2 |
& |
ω2 |
& |
(6) |
|
|
|
|
ϕ – |
|
2 |
ϕ = 0, |
|
|
|
|
|
vф |
|
|
|
|
|
|
|
& |
|
|
|
где |
ω |
& |
|
|
|
|
|
|
= K – пространственная частота ЭМП. |
|
v&ф
Уравнения (5) и (6) называют векторным и скалярным волновыми уравнениями. Они широко применяются на практике для расчета разнообразных электротехнических и радиотехнических устройств, входящих в состав различного радиоэлектронного оборудования и приборов.
Излучатель Герца
Условно можно считать, что этот излучатель представляет собой малый отрезок провода, по которому течет гармонически изменяющийся ток. При расчетах такой излучатель можно считать материальной точкой с гармонически изменяющимся электрическим дипольным моментом. Аналитические выражения для распределения электродинамических потенциалов вокруг этого излучателя являются фундаментальными решениями уравнений (5) и (6)
|
|
|
|
& |
|
µ&a |
& |
& |
|
|
|
µ&a |
& |
|
|
|
& |
||||
|
|
|
|
A = − |
4πR |
jωP exp(− jKR) = |
4πR |
Il exp(− jKR), |
|||||||||||||
|
|
|
|
|
|
|
|
& |
I&l |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
P = |
jω |
, |
|
|
|
|
|
|
|
|
|
|
|
|
1 R |
|
∂ |
|
-1 & |
& |
|
|
1 |
|
|
P&R |
|
& |
P&R |
& |
|||||
ϕ= − |
& |
|
|
|
|
(R |
P exp(− jKR)) = |
|
& |
( |
|
3 |
|
+ jK |
|
2 |
)exp(− jKR) . |
||||
& |
4πεa R |
|
∂R |
|
|
|
|
4πεa |
|
|
R |
|
|
|
|
R |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Элементарный магнитный излучатель
Условно можно считать, что этот излучатель представляет собой контур с гармонически изменяющимся током. При расчетах такой излучатель можно считать материальной точкой с гармонически изменяющимся магнитным диполь-
85