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

книги из ГПНТБ / Тозони О.В. Расчет трехмерных электромагнитных полей

.pdf
Скачиваний:
66
Добавлен:
24.10.2023
Размер:
11.32 Mб
Скачать

Согласно формулам (3.314) и (3.316) получаем:

г-0

 

[пм,Н(М)], Vq р

.~i^erQM

cISm

№РеНе (Q), ___ & L &

rQM

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

^ - - +

(3.318)

cISm -\-

 

+ {[пм, 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

 

Ч

 

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),

если

V{.

 

 

2 5 0

Аналогично для вектора Е (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)]) —

251

 

С - - -

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

Г -

/ e_/

[Пм, н (Af)], Vq^ —

Q 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

r Q M

r Q M

dSM-

(3,322)

Аналогично выводим

N . £ Ш = - e^ p— [«Q, # m -

l

2 я (e£ + e,) s

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)

252

Соотношения (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)

253

Эти соотношения необходимо еще дополнить уравнением

 

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)

264

Чтобы найти выражения для Не и Ее, аналогичные фор­ мулам (3.314) — (3.317), предположим, что поле вне облас­ ти Ve отсутствует. Тогда, чтобы поле в осталось преж­ ним, необходимо ввести на поверхности 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

 

S

 

\

r QM

 

 

 

 

 

 

 

 

 

 

 

/МЦо

\nWHe т \_ dSM

 

 

 

 

 

'О М

 

 

 

f Й*.

 

 

 

 

 

 

 

£ е (Q),

если Q £ Ve\

(3.338)

 

 

 

,0,

если

Q £ Уг.

 

Это отличие

 

выражении

(3.338),

состоит

в том, что в

кроме касательных составляющих Е и Н на S, входит и

нормальная составляющая Е на S. В этом смысле выра­

жение для Ее в квазистадионарном поле является более

2 5 5

сложным, чем в быстропеременном. Выражение же для Е(, на­ оборот, упрощается для квазистационарного поля. В самом

деле для Е{ справедливо выражение типа (3.338) с той лишь

что вместо

1

разницей,

ядра----- нужно использовать ядро

 

 

rQM

 

r

™ v w rQM

 

 

rQM

Кроме

того, для

внутренней (проводящей) среды

(пм, Е[ (Л4))==0. Учитывая все это, находим для Е{следующее выражение:

 

 

 

—(1+/) V

r QM

+

[Пм, н ( ((И)]

rQM

dSм

 

 

 

E( (Q), если Q £ v c;

(3.339)

[О, если Q £ Vg.

Выражение для Ht (Q) будет совпадать с формулой

(3.319), если в последней положить е* = - ? - : /(О

Jffl.

lnM, H t (M)],

VQ

|ге -П+/)] /

rQM

dS/л

 

 

 

r Q M

 

).

 

£

.

.

-uV T ' v*

 

9

Iпм, в, (M)] И

------------------ dS„ +

 

 

 

 

rQM

 

 

1

 

{[ПМ, Et (M)],

VQ) V0 U

- 0 +/, V aJF rQM

s

---------------------rQ M »

$

 

\

(ja>[iHt (Q), если Q^V

10, если Q £ V e.

dSM —

(3.340)

256

Учитывая,

что

 

 

 

 

[Ям, Tit (Л1)] = [Ям, Tie (М)1 = \пм, П (Л4)[ и

IЯм, % (М)] - [Ям, X

(М)] =

[Ям, Ё (М)[,

и складывая соотношения (3.337) и (3.340),

находим

 

 

 

 

 

 

<оцу

 

 

 

 

р,е

' 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

 

 

 

17 4-691

257

+

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) имеет м е н ь ш у ю р а з ­ м е р н о с т ь , так как состоит из одного скалярного и од-

258

ного векторного уравнений на 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]. Приближенные краевые условия на границе проводящей среды можно ввести в тех случаях, когда проводящее тело является массивным и поверхностный эффект резко прояв­ лен, либо когда проводящее тело представляет собой тон­ кую проводящую ферромагнитную оболочку.

Рассмотрим те задачи, в которых проводящее тело яв­ ляется массивным и поверхностный эффект резко проявлен.

При резком поверхностном эффекте электромагнитная вол­ на, преломляясь на границе проводящего тела, будет про­ никать в его толщу в направлении, близком к нормали. Эта волна, вообще говоря, неплоская, но если проводник мас­ сивный и поле на его поверхности медленно меняется от точки к точке в масштабе длины волны внутри тела, то радиус кри­ визны фронта электромагнитной волны внутри проводника будет велик по сравнению с длиной волны и глубиной про­ никновения. Поэтому с высокой степенью точности можно

17*

259

Соседние файлы в папке книги из ГПНТБ