Согласно формулам (3.314) и (3.316) получаем:
г-0 |
|
[пм,Н(М)], Vq р |
.~i^erQM |
cISm — |
№РеНе (Q), ___ & L & |
rQM |
4л |
J |
|
|
|
i |
г |
_ ~ |
— 1> |
|
|
|
|
|
, - l b / Q M |
|
|
|
|
'qm |
|
— ([ПЛ1, E (M)], Vq) Vq |
e~ik/QM\ \ , o |
|
—------- |
>аг>м — |
|
|
|
|
r QM |
) ) |
|
|
j®VeHAQ)’ еСЛИ |
Q £Ve< |
|
О, |
если Q £ |
|
ftfl |
’ _ |
/ |
p—ikirQM |
[Пм, н (M)], Vq Pi |
r |
4n |
- ф V— - |
v-/i> • V1 Г* |
и L |
\ |
rQM |
|
1 C ( _ j . |
—ikirQM |
|
+ Т 5Г $ {lnM,E(M)]kl |
^ - - + |
|
+ {[пм, E (M )], V q) V q |
- i k i r Q M |
cIS m = |
|
|
r Q M |
|
|
|
|
|
|
|
|
|
|
= |
{/®MMQ), |
если |
|
|
|
(3.319) |
|
|
1 |
0 |
, |
если |
|
|
|
|
Складывая соотношения |
(3.318) и (3.319), |
находим |
/сор,Я |
|
[Пм, Я (М)], Vq р |
rQM |
|
|
|
|
|
|
|
|
|
|
— Pi |
r QM |
cIS m — |
|
|
|
|
|
|
|
|
|
|
|
4~ Ф [«м,Я(Л1)], |
|
. - / f t / Q M |
, 2 e,— i k i r QM |
dSu ■ |
V |
|
------ft? |
r |
|
4Я |
Ч |
|
r CM |
|
r QM |
|
|
1 |
- - |
|
|
/ a~ lkerQM |
p~ i k i r QM |
dSM = |
у^ -ф апм .Я ^)], VQjVQj-ly-- |
|
rQM |
|
|
J/®реЯ (Q), |
|
|
|
|
|
если |
QgK,; |
|
(3.320) |
|
|
l/copi^i(Q), |
если |
Q£ V{. |
|
|
Аналогично для вектора Е (Q) получаем
im eE°e( Q ) - - g - § [пм, £ |
(М)], Vqк - |
S |
rQM |
|
-- 8; ---------- |
) dSM 4" |
'ПМ |
|
QM |
|
(3.32 i)
Формулы (3.320) и (3.321) дают единое для всего простран ства интегральное представление для векторов поля через их касательные составляющие на границе S раздела сред. При Q £ 5 ядра последних интегралов в соотношениях (3.320) и (3.321) являются непрерывными функциями, так как при вычитании особенности сокращаются. Поэтому эти интегралы изменяются непрерывно при переходе точки Q через поверхность 5. Вторые интегралы в этих соотношениях можно интерпретировать как разности потенциалов про стых слоев токов, поэтому эти интегралы также непрерывны при переходе Q через поверхность 5. Первые интегралы представляют собой разность вихрей потенциалов простых слоев токов. Поэтому эти интегралы претерпевают скач'ок при переходе Q через S. Предельные значения этих интегра лов могут быть определены по формулам, полностью ана логичным формулам (3.205) — (3.207). Учитывая это, умно жим формулу (3.320) векторно на ~nQ и устремим Q к S со
стороны У(. Тогда
/®Ре [nQ, Н°е(Q)] + (ц, [riQ, Н (Q)] — ре [nQ, Н (Q)]) —
|
С - - - |
I , |
e~ i k/Q M |
p~ i kirQM \ |
4 я |
■Ф [nQ, [Пм, Е (Л1)]] (k\ *-г----------k\ е— ------- ) dSM- |
с |
\ |
' QM |
r Q M ) |
|
|
|
|
~ i k erQM |
|
П<2, ([Пм, £ (A f)], V q) V q |
|
4 я ■Ф |
|
|
r QM |
|
s L |
|
|
|
|
~ l k i'Q M |
dSM = / cojjl, [nQ, H(Q)]. |
|
r Q M |
|
|
|
|
Откуда следует |
|
|
|
|
[я<г,Я(<2)] = |
— |
r - w° |
m |
|
|
Цг + |
\nQ, m |
|
|
|
|
|
1 |
f |
Г _ |
2n |
;e■— |
- j (p |
tiQ, |
((i. |
s |
|
|
|
|
e ~ i k irQ M
— IV |
r Q M |
dSni -j- 2mo(iie + iic) X |
|
л | |
|
/ —ikerQM |
„ a— i k ir Q M |
\ |
X <J}|[nQ, |
[«M, Я ( М ) ] ] Ш - ^ ----------- k V |
■ + |
|
|
'Q M |
r Q M |
) |
|
|
,— i k erQ M |
- > k ir Q M J |
|
+ «Q, (l«M, £ (Al)], Vq) Vq
Аналогично выводим
N . £ Ш = - e^ p— [«Q, # m -
e - i k erQ M
[tlM, £ (M )], V q 8,
r Q M
— e, |
|
— i k ( Q M |
X |
|
|
/ |
|
|
|
|
|
|
|
|
Л0М ' |
d5,M — 2л(0(8e + 8() |
X |
Г f - |
- |
^ |
, / |
2 |
e - i hS Q M |
—/'V qm \ |
X ф ([nQ, [nM, Я (Л4)Ц ^ |
£2 |
-------- ft? - ~ 7 ^ — j + |
+ U<5, ([Пм, |
^ |
|
|
/ - i k / Q M |
- i h ' Q M |
я |
(Af)], Vq)Vq |
r Q M |
r Q M |
cISm• |
|
|
|
|
|
|
(3.323)
Соотношения (3.322) и (3.323) образуют систему инте гральных уравнений К. Мюллера. После решения этой си стемы легко находим векторы поля в любой точке простран ства при помощи (3.320) и (3.321).
Перейдем теперь к квазистационарному полю, т. е. бу
дем считать среду в области V+ проводящей с проводи мостью у и магнитной проницаемостью р. Токами смещения пренебрежем и положим, что электромагнитное поле воз буждается заданными электрическими токами проводи мости, расположенными в V~. Пренебрежение токами сме щения в проводящей среде равносильно принятию диэлек-
„ |
„ |
е( —— |
V |
, |
трическои проницаемости |
внутренней среды |
|
пренебрежение токами смещения во внешнем пространстве эквивалентно принятию диэлектрической проницаемости внешней среды, равной нулю (ее — 0).
Теперь, казалось бы, интегральные уравнения для рас чета квазистационарного электромагнитного поля можно получить из уравнений К. Мюллера (3.322) и (3.323), если
положить в последних ег = и ее == 0. Однако изберем
другой путь по следующим причинам. Во-первых, для ква зистационарного электромагнитного поля справедливы бо лее простые уравнения, чем те, которые получаются из урав нений К- Мюллера при помощи указанной подстановки. Во-вторых, указанная подстановка является формальной и нуждается в некотором обосновании. Дело в том, что соотношение (3.315), на котором основан вывод уравнений К- Мюллера (3.322) и (3.323), не справедливо для внешней среды при условиях квазистационарности, т. е. при ге = 0, а следовательно, и сами уравнения К. Мюллера нельзя счи тать безоговорочно справедливыми при таком значении ее.
Чтобы отчетливо уяснить изложенное, обратимся к урав
нениям (3.295) — (3.298), которые при условиях |
квазиста |
ционарности (ее = 0) имеют вид: |
|
rot Нв = б; |
(3.324) |
rot Ёе = — /сорчД» |
(3.325) |
B e = Toi Ат — grad<p„,j |
(3.326) |
Е, = rot Д. — /ощ0Лт — grad <р,. |
(3.327) |
Эти соотношения необходимо еще дополнить уравнением
|
div Ее — 0. |
(3.328) |
Из уравнений (3.324) и (3.326) получаем |
|
|
rot Не = rot rot Ат = 6. |
(3.329) |
Откуда, если |
положить |
|
|
|
|
div Ат = 0, |
(3.330) |
следует |
|
|
|
|
|
дд» |
= - f c |
(3.331) |
Из уравнений (3.325) и (3.327) выводим |
|
rot Ее = rot rot Ае — /сор0 rot Ат = — /юр0# е — |
= |
— /сор0 rot Ат + |
/сор0 grad cpm. |
|
Откуда находим |
|
|
|
если положить |
ДАе = |
0, |
(3.332) |
|
|
|
|
|
div |
= /<вр0фт . |
(3.333) |
Таким образом, в квазистационарном поле связь между векторным электрическим и скалярным магнитным потенциа лами сохраняется, в то время как связь типа (3.300) между векторным магнитным и скалярным электрическим потен циалами отсутствует, т. е. эти потенциалы и их источники являются независимыми. Следовательно, квазистационарное поле во внешней (непроводящей) среде описывается двумя векторными (магнитным и электрическим) и одним скалярным (электрическим) потенциалами. При этом ска
лярный потенциал |
сре согласно формулам (3.327), |
(3.328) |
и (3.330) удовлетворяет уравнению Лапласа: |
|
|
Дфе = 0. |
(3.334) |
Соотношение (3.326) с учетом уравнения (3.333) прини |
мает вид |
|
|
|
Н е = |
rot Ат----- |
r-J— grad div Ае. |
(3.335) |
Чтобы найти выражения для Не и Ее, аналогичные фор мулам (3.314) — (3.317), предположим, что поле вне облас ти Ve отсутствует. Тогда, чтобы поле в Vе осталось преж ним, необходимо ввести на поверхности 5 простые слои электрических и магнитных токов и простой слой электри ческих зарядов, плотности которых соответственно опре делим по следующим формулам:
I е — [ ЯНе,]> Iт ~ 1^> Ее]> &е ==® 0 (^ » Ее)■
(3.336)
Из соотношений (3.336), (3.335), (3.332), (3.331) получа
ем выражение для Не (Q), которое совпадает с равенством (3.318), если в последнем положить ге = 0:
/(op,,#° ■ 4/Ця> § \[ п м ,Н е № \, |
dS/4 — |
S |
\ r QM |
([пм, Ее(М)], |
\ r QM I |
|
jti>\i0He(Q), если |
(3.337) |
L0, если Q £ V{. |
|
Из формул (3.336), (3.334), (3.332), (3.331), (3.327) по
лучаем выражение для Ёе (Q), которое существенно отли чается от соотношений (3.315) и (3.317):
|
E ° e ( Q ) |
|
1 |
[«М , Ее (М)1, VГ |
( I S m |
— |
|
4я |
|
|
S |
|
\ |
r QM |
|
|
|
|
|
|
|
|
|
|
|
|
|
/МЦо |
\nWHe т \_ dSM |
|
|
|
|
|
4я |
|
'О М |
|
|
|
|
4я |
f Й*. |
|
|
|
|
|
|
|
|
£ е (Q), |
если Q £ Ve\ |
(3.338) |
|
|
|
|
,0, |
если |
Q £ Уг. |
|
|
Это отличие |
|
выражении |
(3.338), |
|
состоит |
в том, что в |
кроме касательных составляющих Е и Н на S, входит и
нормальная составляющая Е на S. В этом смысле выра
жение для Ее в квазистадионарном поле является более
сложным, чем в быстропеременном. Выражение же для Е(, на оборот, упрощается для квазистационарного поля. В самом
деле для Е{ справедливо выражение типа (3.338) с той лишь
„ |
что вместо |
1 |
разницей, |
ядра----- нужно использовать ядро |
|
|
rQM |
|
r |
™ v w rQM |
|
|
rQM |
Кроме |
того, для |
внутренней (проводящей) среды |
(пм, Е[ (Л4))==0. Учитывая все это, находим для Е{следующее выражение:
|
|
|
—(1+/) V |
r QM |
+ |
4я |
[Пм, н ( ((И)] |
rQM |
dSм |
|
|
|
E( (Q), если Q £ v c;
(3.339)
[О, если Q £ Vg.
Выражение для Ht (Q) будет совпадать с формулой
(3.319), если в последней положить е* = - ? - : /(О
Jffl. |
lnM, H t (M)], |
VQ |
|ге -П+/)] / |
rQM |
dS/л |
4я |
|
|
|
r Q M |
|
). |
|
£ |
. |
. |
-u +л V T ' v* |
|
4я |
9 |
Iпм, в, (M)] И |
------------------ dS„ + |
|
|
|
|
rQM |
|
|
|
{[ПМ, Et (M)], |
VQ) V0 U |
- 0 +/, V aJF rQM |
s |
---------------------rQ M » |
$ |
|
\ |
(ja>[iHt (Q), если Q^V
10, если Q £ V e.
Учитывая, |
что |
|
|
|
|
[Ям, Tit (Л1)] = [Ям, Tie (М)1 = \пм, П (Л4)[ и |
IЯм, % (М)] - [Ям, X |
(М)] = |
[Ям, Ё (М)[, |
и складывая соотношения (3.337) и (3.340), |
находим |
|
4л |
|
|
|
|
|
<оцу |
|
|
|
|
р,е |
' Q M |
|
|
|
|
|
|
cISm — |
|
r Q M |
|
|
|
|
|
|
|
|
|
-б + / > |
j/" <»МуИУ. |
|
4 я I |
[Пм, .Б (М)| /шиуе |
|
2~ 'QM |
'Q M |
|
dSM — |
|
|
|
|
|
О |
|
|
Q M |
|
|
|
|
|
в- - П + / > |/ s p r Q M |
I m J T e (Q ), |
Q € V ,; |
|
• ) d S M = |
{ /(0}x^ |
( Q) j |
|
|
r QM |
|
|
|
|
|
|
|
|
|
(3.341) |
Используя выражения (3.339) и (3.341) и проводя аб солютно те же рассуждения, что и при выводе уравнений К- Мюллера (3.322) и (3.323), приходим к следующей систе ме интегральных уравнений для расчета квазистационарного поля:
[TiQ, ТТ(Q)i = |
[Яс, |
m |
(О)] - |
|
|
[ |
M-о ~г И |
|
|
|
|
|[пм, Я(М)], VQ(—|ia_ |
|
<£[я<?, |
2л (Но + и) •П |
L |
|
\ |
r QM |
|
-(1+/>|Л-Р r QM ^ |
|
|
|
|
[ie |
) ] ] |
is |
u - |
|
|
|
|
|
TQM |
|
|
|
HY |
|
- o + o y ^ 'o * |
ФCtlQ, [Ям, E (Af)|| |
|
'Q M |
■cISm -£■ |
2л(Но+Н) |
s |
|
|
|
+ |
2то (ц( |
|
1 |
|
|
|
|
|
-<1+/) У |
(onvVГ т |
|
|
|
— TQM |
dSM', |
(3.342) |
|
|
jj |
|
r QM |
|
|
|
[ПС}, i |
(Q)] |
|
|
|
|
- d + /) ] /Sft 'QM |
|
2n |
ПМ, E (M)], Vq(—------- ------------- |
(ISm + |
|
\ |
rQM |
|
|
„ |
- 0 + i ) V n r rQM |
|
+ |
" f r 9 I»«* I«", Я (Л4)]] ± --------— --------- dSM. |
|
S |
и (3.343), |
|
(3.343) |
Решив уравнения (3.342) |
найдем касательные |
составляющие векторов Е и Н на S, по которым согласно
формулам (3.337), (3.339) и (3.340) можем найти Не (Q),
Et (Q) и Н( (Q) в любой точке пространства. Однако для
нахождения Ее (Q) согласно формуле (3.338) решения систе мы (3.342) и (3.343) недостаточно. Необходимо найти еще
(Пм, Ее (Л4)), для чего требуется решить еще одно инте гральное уравнение. Однако вряд ли это целесообразно,
поскольку поле Ее (Q) само по себе не представляет осо бого интереса.
В отличие от уравнений (3.210) и (3.211), которые явля ются сингулярными, уравнения (3.342) и (3.343) являются фредгольмовскими. Решив уравнения (3.342) и (3.343), сра
зу находим касательные составляющие Е и Н на S, которые часто представляют основной интерес, так как позволяют определить суммарные активные потери. После же решения интегральных уравнений (3.210) и (3.211) для нахождения
векторов £ и Я необходимы некоторые дополнительные вы числения. Перечисленное выгодно отличает уравнения
(3.342) и (3.343) от уравнений (3.210) и (3.211). Однако си стема уравнений (3.210) и (3.211) имеет м е н ь ш у ю р а з м е р н о с т ь , так как состоит из одного скалярного и од-
ного векторного уравнений на S, в то время как система (3.342) и (3.343) состоит из двух векторных уравнений на S. Кроме того, ядра уравнений (3.342) и (3.343) несколько сложнее, чем ядра уравнений (3.210) и (3.211). Все это свидетельствует о том, что систему уравнений (3.210) и (3.211) удобнее решать на ЭВМ, чем систему (3.342) и (3.343).
8. ИСПОЛЬЗОВАНИЕ ПРИБЛИЖЕННЫХ ГРАНИЧНЫХ УСЛОВИЙ НА ПОВЕРХНОСТЯХ ПРОВОДЯЩИХ ТЕЛ ДЛЯ РАСЧЕТА КВАЗИСТАЦИ0НАРН0Г0 ЭЛЕКТРОМАГНИТНОГО ПОЛЯ
Обычный подход к расчету квазистационарного поля в пространстве, частично заполненном проводящей средой, состоит в отыскании решений уравнений поля в прово дящих телах и в окружающей их воздушной среде и в «склеивании» этих решений с учетом краевых условий на поверхностях раздела сред. Основные математические труд ности возникают при решении уравнений поля в проводя щей среде, тогда как интегрирование уравнений поля в воздушной среде всегда можно свести к решению уравне ния Лапласа. Поэтому заслуживают внимания те задачи,
вкоторых можно обойти интегрирование уравнений поля
впроводящей среде и свести весь расчет к решению уравне ний поля во всем внешнем (воздушном) пространстве. По следнее обычно удается достичь, используя приближенные граничные условия на поверхностях проводящих тел [34]. Приближенные краевые условия на границе проводящей среды можно ввести в тех случаях, когда проводящее тело является массивным и поверхностный эффект резко прояв лен, либо когда проводящее тело представляет собой тон кую проводящую ферромагнитную оболочку.
Рассмотрим те задачи, в которых проводящее тело яв ляется массивным и поверхностный эффект резко проявлен.
При резком поверхностном эффекте электромагнитная вол на, преломляясь на границе проводящего тела, будет про никать в его толщу в направлении, близком к нормали. Эта волна, вообще говоря, неплоская, но если проводник мас сивный и поле на его поверхности медленно меняется от точки к точке в масштабе длины волны внутри тела, то радиус кри визны фронта электромагнитной волны внутри проводника будет велик по сравнению с длиной волны и глубиной про никновения. Поэтому с высокой степенью точности можно