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

Применение комфортных отображений в решении некоторых задач электро- и магнитостатики

.pdf
Скачиваний:
42
Добавлен:
01.05.2014
Размер:
270.41 Кб
Скачать

G2, такие, что пара векторов

 

~

 

 

правая.

Тогда граничные

 

 

t,~n

 

условия для векторов напряженности магнитного поля и магнит

-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ной индукции записываются в виде (см. [1], § 60):

 

 

 

 

 

 

 

 

 

B1n = B~1,~n γ

 

= B2n = B~2,~n γ ,

 

 

(28)

 

 

 

 

H1t

=

 

~

 

,t~

 

 

 

= H2t =

 

~

 

 

 

 

 

 

(29)

 

 

 

 

H1

 

 

 

 

H2,t~

.

 

 

 

 

 

 

 

 

 

 

 

 

 

γ

 

 

 

 

 

 

 

 

 

γ

 

 

 

 

 

 

Заметим

,

что граница

 

 

 

 

γ0

 

γ+

 

 

γ

,

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γ =

 

S

 

 

 

 

 

 

 

 

γ =

 

x

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

0

: y

 

[0,1] .

 

 

 

:

 

x > 0 , γ =

 

 

 

±

 

 

 

 

0

 

n 0

 

| |

 

 

o

 

 

 

 

±

y

 

 

 

 

 

 

 

 

 

 

 

 

 

0

,

~t =

1

. Следовательно, соотношения

На γ0 имеем ~n = 1

0

(28), (29) на γ0 принимают вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂A

γ0

= B2n = 0,

 

 

 

 

 

(30)

 

 

 

 

 

 

B1n =

∂y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H1t =

 

 

 

 

 

 

= H2t

= 0.

 

 

 

 

 

(31)

 

 

 

 

 

 

 

∂x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γ0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

На γ

 

имеем ~n =

 

±

1

 

 

~

 

 

 

0

 

 

. Следовательно, соотно-

±

 

 

 

, t =

1

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

шения (28), (29) на γ± принимают вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B1n = ±

∂A

γ

 

 

= B2n = 0,

 

 

 

 

 

(32)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂x

±

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H

 

=

 

 

 

 

 

 

 

 

 

= H

 

 

=

 

4πJ.

 

 

 

 

(33)

 

 

 

 

 

∂y

 

 

 

 

 

±

 

 

 

 

 

 

 

 

 

 

1t

 

γ

±

 

 

2t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Введем в рассмотрение скалярную

функцию

F : G1

→ R,

гар

-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

моническую в G1 и сопряженную гармонической функции A (см. [2],

§ 5, определение 5).

Тогда комплекснозначная функция f(z) =

= A(x,y)+iF (x,y) регулярна в области G1, причем f0(z) = A0

(x,y)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

− iAy0 (x,y). В силу соотношений (30), (31), (32), (33) получаем

 

 

 

 

 

 

 

 

 

 

 

f

 

 

 

 

 

= 0,

 

 

 

 

 

 

 

 

 

 

 

(34)

 

 

 

 

 

 

 

 

 

 

f00 γ0

 

= 4πiJ.

 

 

 

 

 

 

 

 

(35)

 

 

 

 

 

 

 

 

 

 

 

 

γ

±

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

11

Рассмотрим конформное отображение области G1 на верхнюю полуплоскость, которое дается функцией h(z) = s(z2 + 1), где s(z) z – регулярная ветвь многозначной функции z в области C\[0, + ∞), такая, что s(−1) = i. Приведем явный вид этой функции:

 

 

s(z) =

|z| exp

 

2 arg z .

 

 

 

 

 

 

 

 

i

Здесь arg z

 

[0,2π) –

 

p

 

 

 

 

 

 

главное значение аргумента комплексного

числа z. В силу теоремы 1 из § 29 [2] получаем, что функция f представима в виде f(z) = g(h(z)) для всех z G1, где функция g регулярна в верхней полуплоскости. Определим граничные условия для функции g на границе верхней полуплоскости = h(γ). Обозначим

0 = h(γ0) = {u + i0 : |u| > 1} ,

S

± = h(γ+ γ) = {u + i0 : |u| < 1} .

Так как f0(z) = g0(w)h0(z) = g0(w)h(zz), где w = u + iv = h(z), что равносильно z = s(w2 − 1), то в силу (34), (35) получаем

 

 

s(w2 g0(1)

0

= 0,

 

 

 

 

 

 

 

 

 

w)

 

 

 

 

 

 

 

 

 

g0(w)

 

 

 

 

= 4πiJ.

 

 

 

 

 

 

w

 

±

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как для w = u + i0 ± справедливо

равенство

s(w2 − 1) ± = p1 − u2 exp

 

2 arg(u2

− 1) =

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= p1 − u

2

exp

 

 

 

 

 

 

2

то для |u| < 1 получаем

(36)

(37)

p

= i 1 − u2,

g0(u + i0) = (Re g)0u(u,0) − i(Re g)0v(u,0) = √4πJu 2 1 − u

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

12

(Re g)u0 (u,0) =

 

0,4πJu

|u| > 1,

 

 

 

 

 

 

 

 

 

, |u| < 1.

 

 

 

 

1

u2

 

 

Отсюда получаем

 

 

 

 

 

 

 

 

 

 

 

 

const ,

 

 

 

 

 

 

 

 

 

u

> 1,

 

Re g(u,0) = ( const

 

 

 

 

 

 

 

|u|

< 1.

(38)

4πJ 1

u2,

 

 

 

 

 

 

 

 

| |

 

 

p

Поскольку для определения векторного поля магнитной индукции в области G1 функцию f достаточно знать с точностью до некоторой константы, то в соотношениях (38) без ограничения общности считаем const = 0.

Итак, для определения гармонической в верхней полуплоско-

сти функции Re g получаем задачу Дирихле:

 

 

 

 

 

 

 

 

 

(Re g(u,v)) = 0,

0,

v > 0,

u > 1,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Re g(u,0) =

 

 

 

 

 

 

 

 

| |

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(39)

 

 

 

(

−4πJ 1 − u2, |u| < 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

см

.

теорему

5

Решение

этой задачи дается формулой Пуассона

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

§ 29 [2]). Получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Re g(u,v) = −4Jv Z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 t2

 

 

 

 

 

 

 

 

 

 

 

dt,

 

v > 0.

(40)

 

(t − u)2 + v2

 

 

 

−1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пользуясь вычисленным при решении задачи 1 интегралом (18),

значение которого дается формулой (24), получаем

 

 

p|

 

 

|

 

2

Re g(u,v) = 4πJ v

 

 

 

 

 

 

cos

ψ(u,v) − φ(u,v)

 

 

 

1

 

w2

 

 

,

где функции φ(u,v) и ψ(u,v) определены в формулах (22) и (23). Таким образом, получаем вид магнитного потенциала в обла-

сти G1:

− | |

 

2

 

 

A(z) = 4πJ Im h(z)

z

cos

ψ(z) − φ(z)

,

 

где

 

 

 

 

13

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h(z) = p|1 + z2| exp

 

arg 1 + z2

,

 

 

 

2

!

 

 

 

 

 

(1

 

Re h(z))2 + (Im h(z))2

 

φ(z) = arccos

p

 

 

 

1 − Re h(z)

 

 

,

 

 

 

 

 

 

 

ψ(z) = arccos

 

 

1 + Re h(z)

 

 

! .

 

 

 

 

 

 

 

 

(1 + Re h(z))2 + (Im h(z))2

-

 

 

силовые линии магнитного поля в обла

На рис. 5 изображеныp

 

 

 

 

 

 

 

сти G1 вида A(z) = const .

Рис. 5. Силовые линии магнитного поля в области G1

14

Задача 3. На плоскости R2 идеальный проводник заполняет замкнутую область:

G2

= n y

: x2 + y2 ≤ 1 o ,

 

x

 

как показано на рис Вектор намагниченности проводника ~

. 5. J

сонаправлен с осью ординат. В области G1 = R2\G2 проводники отсутствуют. Требуется найти магнитное поле в области G1.

G1

y

 

i

 

 

 

 

G2

~

 

 

J

−1

0

1 x

 

−i

 

Рис. 6.

В области 1 векторы напряженности магнитного поля ~ 1 и

G H

~

совпадают и удовлетворяют уравнениям

магнитной индукции B1

 

~

~

= grad A,

Максвелла (25), (26). Векторное поле B1

ищем в виде B1

где скалярная функция A: G1 → R имеет смысл магнитного потенциала. При этом уравнение (26) автоматически выполняется, а уравнение (25) превращается в уравнение Лапласа для A в обла-

сти G1 (27).

Осталось определить, каким граничным условиям на границе γ = ∂G1 области G1 должна удовлетворять функция A. Пусть

~

~

 

 

 

 

H2 и B2 – векторы напряженности магнитного поля и магнитной

индукции в области G2 соответственно.

~

~

~

Тогда B2

= H2

+ 4πJ =

~

~

~

 

 

 

= 0.

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

= −4πJ.

 

 

 

Пусть ~n – единичный вектор внешней нормали к границе обла-

сти 2 а ~ единичный касательный вектор к границе области

G , t –

2 такие что пара векторов ~ правая Граничные условия

G , , t,~n .

для векторов напряженности магнитного поля и магнитной индукции записываются в виде (28), (29).

15

Граница области G1 представляет собой окружность с цен-

тром в нуле радиуса 1 вида γ =

e: ϕ [0,2π)

. На γ имеем

cos ϕ

~

 

 

 

sin ϕ

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

 

соотношения

 

 

~n = sin ϕ , t =

− cos ϕ .

,

(28),

 

 

 

 

 

 

 

(29) принимают вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂A

 

∂A

 

 

 

 

 

 

 

 

 

 

B1n =

 

 

cos ϕ +

 

 

sin ϕ |γ

= B2n = 0,

 

 

(41)

∂x

 

∂y

 

 

H1t =

 

∂A

 

 

∂A

 

 

= H2t = 4πJ cos ϕ.

(42)

 

∂x sin ϕ −

 

∂y cos ϕ γ

Введем в

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

рассмотрение скалярную

функцию

F : G1 → R,

гар

-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

моническую в G1 и сопряженную гармонической функции A (см. [2], § 5, определение 5). Тогда комплекснозначная функция f(z) =

= A(x,y)+iF (x,y) регулярна в области G1, причем f0(z) = A0

(x,y)

− iAy0 (x,y). В силу соотношений (41), (42) получаем

x

 

 

f0

γ e= 4πiJ cos ϕ.

 

 

(43)

Рассмотрим конформное

отображение области G

1

на верхнюю

 

 

 

 

 

полуплоскость, которое дается дробно-линейной функцией:

 

h(z) = zz + 11 i .

В силу теоремы 1 из § 29 [2] получаем, что функция f представима в виде f(z) = g(h(z)) для всех z G1, где функция g регулярна в верхней полуплоскости. Определим граничные условия для функции g на границе верхней полуплоскости = h(γ). Так

как

 

 

 

 

 

 

2i

 

f0(z) = g0(w)h0(z) = g0(w)

 

,

 

 

(z + 1)2

где w = u + iv = h(z), то из соотношения (43) получаем

g0(w) h0

 

 

ϕ

(44)

e

 

e= 4πiJ cos ϕ.

Так как w| = u + i0 = h(e ) = − tg

 

, а

 

2

 

h0 ee

=

 

i

 

 

,

 

1 + cos ϕ

 

то соотношение (44) перепишется в виде

16

g

(u + i0)

 

ϕ = 4πJ cos ϕ(1 + cos ϕ) = 8πJ

1 − u2

.

0

 

u=−tg 2

 

 

 

 

1 + u2

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

для любого u R получаем

 

 

 

 

 

 

(Re g)0 (u,0) = 8πJ

1 − u2

,

 

 

 

 

 

u

1 + u2

 

 

 

откуда

 

 

 

 

 

 

 

 

 

u

 

 

 

 

 

 

 

 

 

 

 

 

 

Re g(u,0) = 8πJ

 

+ const .

(45)

 

 

 

1 + u2

Поскольку для определения векторного поля магнитной индукции в области G1 функцию f достаточно знать с точностью до некоторой константы, то в соотношении (45) без ограничения общности считаем const = 0.

Итак, для определения гармонической в верхней полуплоско-

сти функции Re g получаем задачу Дирихле:

 

 

 

 

 

 

(Re g(u,v)) = 0,

 

v > 0,

 

 

(46)

 

 

 

 

 

 

u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Re g(u,0) = 8πJ

 

 

 

.

 

 

 

 

 

1 + u2

 

см теорему

 

Решение этой

задачи дается формулой Пуассона

 

 

 

 

 

 

 

 

 

 

 

(

.

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

§ 29 [2]). Получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

+∞

 

 

 

 

 

 

 

 

 

 

Re g(u,v) = 8Jv Z

 

t

dt,

 

v > 0.

(47)

(1 + t2)((t − u)2 + v2)

 

 

−∞

 

 

 

 

 

 

 

 

 

 

Нашей целью теперь является вычисление интеграла

 

 

 

 

+∞

 

 

 

 

 

 

 

 

 

 

 

I = Z

t

dt,

 

 

(48)

 

(1 + t2)((t − u)2 + v2)

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

для произвольных u R и v > 0.

Для вычисления этого ин-

теграла воспользуемся теорией вычетов.

Введем комплексную

функцию:

 

z

 

 

 

L(z) =

 

 

 

.

 

 

 

(1 + z2)((z − u)2 + v2)

 

Тогда в силу формулы (21)

из § 13 [2] получаем

 

 

 

 

 

 

 

I = 2πi

res L + res L

.

 

(49)

 

 

u+iv

i

 

 

 

 

 

17

 

 

 

 

В силу леммы 2 § 13 [2] находим

 

 

 

 

 

= (i − u)2

+ v2 2i .

 

u+iv

1 + (u + iv)2 2iv ,

 

 

i

 

 

 

res L =

 

u + iv

 

 

 

 

 

 

1

 

 

 

 

res L

 

 

 

 

 

 

 

i

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя полученные значения вычетов в (49), получаем

 

I = π (1 + (u + iv)2)v

+ (i − u)2 + v2 =

 

 

 

 

 

 

 

 

 

u + iv

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= u + iv − i

(i + u + iv)v

i − u + iv

=

 

 

 

 

 

π

 

 

 

 

 

 

 

u + iv

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

=

 

 

 

 

π(iu − u2 − iuv)

=

 

 

 

πu

.

 

 

 

 

 

 

 

 

 

 

 

 

 

(u2 + (v + 1)2)v

 

 

 

(u + iv

i)((v + 1)2 + u2)v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Следовательно, из (47) получаем для любых u R и v > 0:

 

 

 

 

 

Re g(u,v) =

 

 

 

8πJu

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

u2 + (v + 1)2

 

 

 

 

 

 

 

 

Так как u + iv = h(z) =

z−1

i для z

 

G

, то

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z+1

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u =

Im

z − 1

, v = Re

z − 1

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z + 1

 

 

 

 

 

 

z + 1

 

 

 

 

 

 

 

то получаем вид магнитного потенциала в области G1:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Im z−1

 

 

 

 

 

 

 

 

 

 

A(z) = −8πJ

 

 

 

 

 

 

 

 

z+1

 

 

 

.

 

 

 

 

 

 

z 1

 

2

 

 

 

 

 

 

 

z 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z+1

 

+ 1 + 2 Re z+1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

На рис. 7 изображены силовые линии магнитного поля в обла-

сти G1 вида A(z) = const .

18

Рис. 7. Силовые линии магнитного поля в области G1

Задачи для самостоятельного решения

Задача 4. На плоскости R2 идеальный проводник заполняет замкнутую область:

G2

= n y

: y ≤ 0 o S n y

: x2 + y2 ≤ 1, y ≥ 0 o ,

 

x

x

 

как показано на Вектор намагниченности проводника ~ сона

8. J -

правлен с осью абсцисс. В области G1 = R2\G2 проводники отсутствуют. Требуется найти магнитное поле в области G1.

Задача 5. На плоскости R2 идеальный проводник заполняет замкнутую область:

19

y

 

i

 

G1

 

−1

1

0

x

G2

~

 

J

G2

= n y

Рис. 8.

: x2 + y2 ≤ 1, y ≥ 0 o ,

: y ≤ 0 o S n y

 

x

x

 

как показано на Вектор намагниченности проводника ~ сона

9. J -

правлен с осью ординат. В области G1 = R2\G2 проводники отсутствуют. Требуется найти магнитное поле в области G1.

y

 

i

 

G1

 

−1

1

0

x

G2

~

 

J

Рис. 9.

Задача 6. На плоскости R2 поляризованный диэлектрик заполняет замкнутую область:

G2

= n y

: y ≤ 0 o S n y

: x2 + y2 ≥ 1, y ≥ 0 o ,

 

x

x

 

как показано на Вектор поляризации диэлектрика ~ сонапра

10. P -

влен с осью ординат. В открытой области G1 = R2\G2 вещество отсутствует. Требуется найти электрическое поле в области G1.

Задача 7. На плоскости R2 идеальный проводник заполняет замкнутую область:

G2

= n y

: x2 + y2 ≤ 1, y ≥ 0 o ,

 

x

 

20