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

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

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

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

ФDn(М) dSM

--------q У

с-------

 

 

Поэтому, для выполнения равенства (1.198), необходимо:

СV гх&уг2.

Отсюда и из выражений (1.198) и (1.200) получаем:

 

Ф(М ) =

У ехеуег Rqm

( 1.202)

 

 

 

 

 

D„ (М) =

q

(rQM’ пм)

(1.203)

 

4 л V Ь хЬ уЬ г

Rqm

 

 

 

 

Попутно доказано следующее соотношение:

 

Ф

(rQM' п М)

cISm — I

(Rqm' пм) d%M —

 

 

 

Rqm

 

 

0.

если Q вне поверхности S;

 

 

2п]^гхгугг,

если Q на поверхности 5;

(1.204)

 

У^ехгугг,

если Q внутри

поверхности S;

 

Из выражения (1.202) согласно принципу суперпозиции

для потенциала объемнораспределенного заряда

и потен­

циала простого слоя зарядов, распределенных по некоторой поверхности S, получаем:

Ф(Q) =

____1_

Р (М) dSм,

(1.205)

 

4 л У ех еУ&2 ^

R qm

 

Ф (Q) =

Ф

а (М) dSм-

(1.206)

4л | ' ехеуег

 

Rqm

 

Объемный потенциал (1.205) удовлетворяет уравнению (1.195) в области V, занятой зарядами, и уравнению (1.196) вне V. Потенциал простого слоя (1.206) везде вне S удовле­ творяет уравнению (1.196). На поверхности5 этот потенциал непрерывен, а нормальные составляющие вектора смещения поля простого слоя зарядов претерпевают скачок. Найдем

нормальные составляющие вектора D с наружной и внутрен­ ней сторон поверхности S. Для этого вокруг некоторой точки Q б S вырежем элемент поверхности А5 (рис. 17)

80

столь малых размеров, чтобы можно было пренебречь его кривизной и считать его плоским. Все поле заряженной

поверхности S разобьем на два слагаемых: поле D(I) от за­

рядов на поверхности S — AS и поле D<2) — элементар­ ного заряда на AS. Согласно соотношениям (1.203) и (1.206)

W (Q) =

Уехеуег sj AS 0 {М)

(rQMп м^

cISm-

(1.207)

 

 

Щм

 

 

 

Картина

поля

D{1) симметрична

 

 

 

 

 

по обе стороны от элемента AS вслед­

 

 

 

 

ствие пренебрежения его кривиз­

 

 

 

 

ной.

Поэтому

из

соотношения

 

 

 

 

 

Dnl D(ni

= о (Q)

получаем

 

 

 

 

 

 

 

 

 

 

 

 

(1-208)

 

 

 

 

В

действительности симметрия

 

 

 

 

 

поля Dl2)

несколько

нарушается

 

 

 

 

 

вследствие

реально существующей

 

 

 

 

 

малой

кривизны элемента

AS, по-

 

 

Рис.

17.

 

этому в соотношении (1.208) по­

 

Из

выражений

ставлен приближенный знак равенства.

(1.207) и (1.208) получаем:

 

 

 

 

 

 

 

 

 

Den(Q) =

V (Q) +

(Q) ~

 

+

 

 

+

4*VW z

l

 

<*№ ■

(rQM'

q)

■dSM;

(1.209)

 

s_ &s

 

 

Щм

 

 

 

 

Dln (Q) =

D'n (Q) +

 

Dg> (0 ):

p(Q)

+

 

 

+

 

l

j

 

 

or(AJ)

(rQM> n Q)

cISm.

( 1. 210)

 

 

У ехеуег S—AS

 

R q m

 

 

 

Формулы (1.208) — (1.210) тем точнее, чем меньше кривизна элемента AS, т. е. чем меньше сам элемент. По­ этому в пределе получаем:

Den (Q)

 

g(Q)

I

1

(rQM< п м)

cISm',

 

2

~r

У exeyez

Щм

 

 

( 1. 211)

Dn (Q)

g(Q)

 

irQM' пм)

(ISm-

2

 

У гхгув2

Rqm

 

 

 

( 1. 212)

6 4-691

81

Пусть среда внутри поверхности S является проводящей.

Тогда Dln (Q) э 0 и из соотношения (1.212) для поверхнос­ тной плотности электрических зарядов HaS получаем следую­ щее однородное интегральное уравнение (уравнение задачи Робэна) для анизотропной среды:

« Р Э - 6П у ЬхЬуЬг V

KqM dSM~ 0- ( |'213)

Уравнение (1.213) имеет ненулевые решения, каждое из которых однозначно определяется, если задать полный (суммарный) заряд q проводника:

 

 

j ) o { Q ) d S Q =

q.

 

(1.214)

 

 

s

 

 

 

 

Приведем уравнение

(1.213),

использовав

равенство

(1.214),

к виду

 

 

 

 

 

 

<*(<2) —

 

(^)а(М)

 

(rQM’ nQ>

 

 

2п V вхвувг

 

Rqm

 

 

 

 

 

 

 

2п У Bx & yBz

 

т т - ,

(1.215)

 

 

 

d S u =

или

 

 

 

 

 

 

 

a(Q)

=

ф a (М)

 

(rQM’ V

 

 

 

V вхвувг

 

 

Rqm

 

_ 1_

(rQP>nQ) dSp

q

^

rQP■(nQ)

d S P.

S

R 3q p

 

 

 

Ri,

(1.216)

 

 

 

 

 

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

82

В случае плоского поля соотношения (1.205), (1.206), (1.211) и (1.212) имеют вид:

ф (Q) =

2 л /

 

|р (М ) In

■dSM;

(1.217)

 

гхгу

 

 

R,QM

 

Ф (Q) =

J-

Ф О (М) In ~

dlM\

(1.218)

 

2 я V

ехеу

J

 

 

K q m

 

Dn (Q) =

 

 

-фа(Л1)-

(rQM' пм ) й1м\

(1.219)

 

2 я /

ехеу

 

 

Щм

ct(Q)

 

= ^ & о ( М )

{rQM' Пм),. й1м,

DUQ)

 

 

2я /

8хЕ(у

 

г?2..

(1.220)

 

 

 

гхги J

 

 

где

 

 

 

 

 

 

 

^?QM=

 

(*м ~ V 2

+

(Ум — Уо?

 

 

 

 

е*

 

 

%

 

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

ности 5 раздела сред, а другое—в области V+, занятой ани­ зотропной средой. Это связано с тем, что потенциал ф

удовлетворяет различным уравнениям в областях У- и V+. В области V~ потенциал ф удовлетворяет обычному уравне­

нию Лапласа или Пуассона, а в области —уравнению (1.195) или (1.196). Поэтому при использовании для ф еди­ ного во всем пространстве интегрального представления

приходится в области V+ вводить в качестве вторичных источников объемные заряды. В самом деле уравнение (1.196) можно преобразовать к виду

+ .

•<Ф

d2(p+

I /с+

с+' д2Ф+

ду г

(е г

г х ) д г 2

Аф+

 

ео

 

( 1.221)

 

 

 

 

из которого видно, что в качестве вторичных источников необходимо в области V+ ввести объемные заряды с

6*

8 3

плотностью

02ф.+

32ф

э+ = •

и+ (8^ - «ф ду> + (8,+-

8t) ~дг*~

а это приводит к интегральному уравнению по области Предложим другой подход к составлению интегральных

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

из областей V~ и V+ выражения для ср. Иначе говоря, откажемся от единых для всего поля вторичных источников, и напротив для поля в каждой из областей V~ и V+ вводим на поверхности S свои вторичные источники, распределение плотностей которых должно быть таковым, чтобы на поверх­ ности S выполнялись краевые условия для касательных

и нормальных составляющих векторов Е и D. Эти краевые условия и приводят к системе двух интегральных уравне­ ний на поверхности 5 относительно плотностей вторичных источников.

Для простоты изложения будем считать, что объемные заряды находятся только в области V~. Потенциал в облас­

ти V~ обозначим через ф~, в области

—через ср+ . Тогда:

32ср

. д2<р

, д2ф _

дхг

ду%

3z2

- pfe в области Vk k = 1, 2, . . . , п

» ( 1.222)

О в области V Vh\

+ д2ф +

+ д \ +

+ д2Ф+

= 0 в области V+ (1.223)

х дх1

у ду2 ^ 2

дг2

 

Условия непрерывности касательных составляющих на­ пряженности и нормальных составляющих смещения будут выполнены, если на поверхности 5 будут удовлетворяться краевые условия:

 

 

Ф + _ Ф-

= С;

 

(1.224)

4

cos (л, *) + е+

cos (я,

у) +

 

+

4

cos (я,

г) = е-

,

(1.226)

где С — произвольная постоянная.

84

Потенциал ф в области V будем искать в виде

Ф (Q) =

V

ГPfe (At)

dVм -(-

 

4пе~ / й

l rQM

 

 

 

 

(r Q M ' n M>

dSM i

(1.226)

 

 

,3

rQM

т. e. в качестве вторичных источников для поля вне S выби­ раем двойной слой электрических зарядов — слой диполей.

Потенциал в области V+ представим в виде

Ф + (Q ) =

а (М) dSм,

(1.227)

4я ]/'е+ е+ е+

Rqm

 

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

При таком выборе выражений для потенциалов ф+

и ф_ уравнения (1.222) и (1.223) будут удовлетворяться. В качестве вторичных источников для поля в изотроп­ ной среде можно было бы выбрать простой слой электри­ ческих зарядов. При этом получили бы систему двух ин­ тегральных уравнений, одно из которых было бы уравне­ нием 1-го рода. Сделанный нами выбор вторичных источни­ ков продиктован стремлением получить систему двух ин­ тегральных уравнений 2-го рода. Для вывода системы воспользуемся краевыми условиями (1.224) и (1.225). Из выражения (1.224) и теоремы о предельных значениях

потенциала двойного слоя находим

Т (Q) + - ^

ф т (М) —

!f -n

dSM+

 

in

I

 

 

ГЬм

 

 

Н------- - J L

= ( f )

q(M) dS,

 

2 n V 6 + 8 + 8 J+

Rqm

м

 

 

 

1

у

f

_£i B

dVм

C.

(1.228)

M

,

rQM

К

85

Из выражений (1.225), (1.212) и формулы для нормаль­ ной производной потенциала двойного слоя [6] выводим

<т(<2)

X

 

 

 

&а(М )

Rqm

dSM +

 

У

е+е+е+

*

 

 

 

s

 

 

 

 

 

*

°г

 

 

 

 

+ ■

ф

[т (М )

Т (Q)] х

 

 

(п М< V

 

3 (rQM’ n Q)

(rQM< п м )

<1Sm

 

гз

 

 

'QM

 

 

rQM

 

 

 

 

 

= - h r

2

 

('q-r

dV*-

(i -229)

2lt

1

 

 

гЬм

 

 

Уравнения (1.228) и (1.229) образуют систему интеграль­ ных уравнений 2-го рода, решение которой необходимо для определения плотностей вторичных источников т (М) и о (М). Когда эти источники найдены, то потенциал,

а вместе с ним Е и D могут быть найдены в изотропной и анизотропной части пространства при помощи соотноше­

ний (1.226) и (1.227)

соответственно.

 

разрешима,

то

Если система уравнений

(1.228) и (1.229)

ее решение неединственно. В самом деле,

если

т (Л1)

и а (М )—какое-либо решение этой системы,

то т (М) +

А

и

а (М)+ а0 (М),

где А — произвольная

константа,

а

а0 (М ) — ненулевое

решение

уравнения (1.213),

также

будут решением системы (1.228) и (1.229). Последнее объяс­ няется тем, что как потенциал двойного слоя постоянной плотности, распределенной по замкнутой поверхности S, не создает поля в окружающем пространстве, так и потен­ циал (1.227) с плотностью а0 (М ) не создает поля напряжен­ ности внутри поверхности 5. Учитывая это, преобразуем уравнения (1.228) и (1.229) методом, изложенным в гл. I, 4, 6 так, чтобы

(j) т (М) dSM= 0;

(1.230)

s

 

$ о (М) dSM= 0.

(1.231)

s

 

86

Преобразованные

уравнения

имеют вид:

 

 

т (Q) +

 

 

ф т (М)

(rQM< Пм)

 

 

 

 

 

 

м

 

 

 

 

 

 

s

 

 

 

 

____ 1_ £

 

(ГР М п М )

dSF dSм ■

 

 

 

S

J

 

r3Di

 

 

 

 

 

 

S

 

р м

 

dSc

 

 

 

 

а (М)

 

dS,

2 я ] / е+е+е+

 

 

 

 

'■PM

м

S

 

 

 

 

 

у ° г

 

 

 

 

 

 

р * ji

 

 

l

 

 

dVM; (1.232)

 

 

'QM -

ж

PM

 

 

 

 

a(Q)

2л ]/

 

 

ф а (УИ)

(rQM> wq)

__

 

 

8+8+8+ s-

 

Rqm

 

2 пУ 8+8+8+

^Sm +

■Ф

 

 

 

 

 

 

 

It (M ) - x ( Q ) ] x

M

(n M> Hq )

 

3 ( r Q M ’ ” q ) (r Q M ’ n M>

dSм

'Q M

 

 

,5

 

 

 

 

 

^M

 

 

 

 

= - 2Г 2

f

Pk(M) ^

'

nQ) - dFM.

(1.233)

 

A=1 i/

 

rQM

 

 

Справедлива следующая теорема:

Дели решение системы

интегральных уравнений (1.232) и (1.233) существует, то это решение единственно.

Для того, чтобы решение этой системы было единствен­ ным, необходимо и достаточно, чтобы соответствующая однородная система:

то ( Q )

+

 

 

irQM’ Пм)

 

 

 

Г3

 

 

 

 

 

 

rQM

 

 

1

{Грм> пм)

dSM

 

s

 

ГРМ

 

 

•S

 

 

 

 

2л ]/ 8+8+8+

о0 (М)

1

 

1 eg

dSM = 0;

£

 

Rqm

s J

Rpm

(1.234)

87

<у0(Q)

 

('■qm- wq)

2" У е+е+8+ ,

а 0 ( М )

 

 

QM

 

2яК в+е+е+

^ К (Л4) т0(Q)] х

 

-j-

 

 

S

 

{пм , 4 q )

3 ( r Q M . n Q)

(rQM, п м )

О (1.235)

X

 

d S M =

rQM

'QM

 

имела только нулевое решение. Пусть т0 (М ) и ст0 — какое-либо решение системы (1.234) — (1.235). Рассмотрим функции:

фо (Q) - - 4 ^

ф to (М) -— Q-M-

- Ci;

(1.236)

 

'qm

 

 

Фо (Q) =

&

^ - d S M,

(1.237)

J

"■ом

 

 

4я]/ е+е+е+ /

R QM

 

где

c . = 4 s h r < h < M )

(rPM■Пм1 dS, d5^ (1.238)

r 3

r PM

— постоянная.

Функция Ф0 (Q) удовлетворяет во всем пространстве, за исключением поверхности S , уравнению Лапласа:

 

 

АФо(0) =

0,

 

 

(1.239)

а Фо (Q) — уравнению

 

 

 

 

 

 

 

с+ б2Фо

I

+

<^2фо

|

+

<Э2фо _

0

(1.240)

дх*

h

^ Г -

+

е*

-

и-

Из уравнений (1.234) и (1.235) вытекают краевые усло­

вия для Ф0 (Q) и ф0 (Q) на поверхности S:

 

 

 

фо (Q )

— Фо (Q ) = С 2;

 

(1.241)

(5Фе.

d(fi

 

 

 

 

е<Р*

 

 

8

 

 

 

x) + Bt ~ d r cos(n>У) +

 

+ j 4

 

 

 

г),

 

(1.242)

 

+ 8?

дг cos (я,

 

где

—1

<£cx0(M )U -(j)

dSt

С2=

dSjM

4Л V 8 + 8 + 8 +

S

' РМ

X °у °2

 

— постоянная.

Воспользуемся теоремой Гаусса—Остроградского:

j div adV =

(j) (a, n) dS.

 

(1.243)

Полагая a ф0е grad cp0,

/е+

0 0

\

где матрица 8 = 1

0

О

К

 

V

0

0

е + /

находим

 

 

 

 

 

 

 

 

 

div а = ф0div (е grad ф0) + (grad ф0) е grad ф0) =

 

„ /г-и дг% | + 52ф

, <Э2Ф0\

 

 

~ Фо 11ч ~ й г +

Ц ~5$г

+ Ч - Щ - ) +

 

 

+ еТ

5фо_

 

+ (

дфо

,+ /

<?Фо

 

 

дх

+ ч

 

 

+ е?

dz

 

Отсюда, учитывая (1.240), согласно соотношению (1.243)

для

области V+ находим

 

 

 

 

 

 

 

ффо

 

cos (n>

*) +

cos (п> ^) +

 

 

 

 

 

 

 

 

 

J

+

е*" ~ W C0S (п> г)j

d S ^

l

[

8+ ( ^ ] + в + № - ] +

 

дх

 

 

 

 

 

,+

 

I \ 2

dV.

 

(1.244)

 

 

 

 

 

 

 

 

 

 

+ е? { 4

 

 

 

 

Аналогично для функции Ф0 (Q) находим

 

 

 

дФ‘

= - в - |

| grad Фц |2dF,

(1.245)

 

^ ^ o - w

d S

где знак «—» поставлен потому, что нормаль направлена внутрь V~.

39

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