книги из ГПНТБ / Тозони О.В. Расчет трехмерных электромагнитных полей
.pdfСледовательно,
ФDn(М) dSM |
--------q У |
с------- |
|
|
Поэтому, для выполнения равенства (1.198), необходимо:
СV гх&уг2.
Отсюда и из выражений (1.198) и (1.200) получаем:
|
Ф(М ) = |
4я У ехеуег 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) |
||
|
4л У^ехгугг, |
если 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) = |
4лУехеуег 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) + |
D® (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) |
||
|
|
4л У ехеуег S—AS |
|
R q m |
|
|
|
Формулы (1.208) — (1.210) тем точнее, чем меньше кривизна элемента AS, т. е. чем меньше сам элемент. По этому в пределе получаем:
Den (Q) |
|
g(Q) |
I |
1 |
(rQM< п м) |
cISm', |
|
2 |
~r |
4я У exeyez |
Щм |
||
|
|
( 1. 211) |
||||
Dn (Q) |
— |
g(Q) |
|
irQM' пм) |
(ISm- |
|
2 |
|
4л У гхгув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 |
|
|
|
2я 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) |
2я |
M |
, |
rQM |
К
85
Из выражений (1.225), (1.212) и формулы для нормаль ной производной потенциала двойного слоя [6] выводим
<т(<2)
X
|
|
|
&а(М ) |
Rqm |
dSM + |
|
|
2л У |
е+е+е+ |
* |
|
|
|
||
s |
|
|
|
|
|||
|
* |
°г |
|
|
|
|
|
+ ■2я |
ф |
[т (М ) — |
Т (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 |
|
|
|
|
||||
|
у ° г |
|
|
|
|
|
|
|
2л р * 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