Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Хасанов лекции / Эл-статика практика.pdf
Скачиваний:
256
Добавлен:
10.02.2015
Размер:
647.7 Кб
Скачать

25

Соответственно выражение для напряженности поля есть

E~ =

grad ϕ =

 

~ar

+ ~aϕ

 

1

 

+ ~aθ

1 ∂

(

 

 

Er cos θ) 1

 

R33 =

 

"

 

 

 

 

 

 

·

 

 

 

 

 

 

 

#

 

 

 

 

 

 

∂r

 

 

r sin ϕ

∂ϕ

 

 

r ∂θ

 

 

r

 

 

 

 

1 + 2

 

3

cos θ

 

 

1

 

 

 

3

 

sin θ.

 

 

 

 

= ~arE0

R

~aθE0

 

 

R

 

 

 

 

 

3

 

3

 

 

 

 

 

 

 

 

 

r

 

 

 

r

 

 

 

 

Поверхностная плотность заряда равна: σ = −41π ∂ϕ∂r |r=R = 43π E0 cos θ.

Комментарий: так как интеграл от поверхностной плотности заряда равен нулю, безразлично, заряжена сфера или изолирована.

Задание на дом: решить задачи 2.16–2.18 на стр. 33.

1.7Дипольный момент. Тензор квадрупольного момента. Поле системы зарядов на больших расстояниях

Скалярный потенциал произвольной системы зарядов на больших расстояниях определяется следующим выражением:

 

q

~

 

1

3

 

 

 

 

 

 

(dr~ )

X

 

 

 

 

(84)

ϕ =

 

+

 

 

 

+

 

Qαβxαxβ + . . .

 

 

r

 

r3

 

2r5

α,β=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

Здесь q — полный заряд системы, d — дипольный момент системы:

 

~d = Z ~rρ(~r )dV

N

 

 

 

(85)

=1 qi~rl,

 

 

 

 

 

 

 

 

 

 

 

 

iX

 

 

 

 

Qαβ — тензор квадрупольного момента.

 

 

 

 

 

Qαβ = Z (3xαxβ − r2δαβ)ρ(r)dV

N

− r(i)2

δαβ

.

(86)

=1 qi 3xα(i)xβ(i)

 

 

 

 

 

 

 

 

 

iX

 

 

 

 

Пример 1.7.1 Плотность заряда в одном из возбужденных состояний атома водорода в сферических координатах имеет вид:

ρ = ρ0r4e−λr sin4 θ,

где , . Вычислить ~ и .

ρ0 λ = const d Qα,β

По определению дипольного момента (85):

Z Z 3

~ ~ ~ ~ ~ X ~

d = ρ(r)r dV = ρ[ix + jy + kz] dV = eαdα.

α=1

Вычисление компонент dα дает значение, равное нулю, т.к. вычисляются

интегралы от нечетных функций в симметричных пределах, т.е. ~ . d = 0

26

Для вычисления компонент тензора квадрупольного момента воспользуемся формулой (86). При α = β = 1 получим:

Q11 = Qxx = Z0∞ Z0π Z0ρ(r)[3x2 − r2]r2dr sin θ dθ dϕ =

(87)

= ρ0

Z0r8e−λrdr Z0π Z0sin4 θ(3 sin2 θ cos2 ϕ − 1) sin θ dθ dϕ.

 

Для интегрирования по углу ϕ воспользуемся выражением:

 

 

1

 

 

 

Z0

cos2 ϕ dϕ =

 

Z0

(1 + cos 2ϕ) = π.

 

 

2

 

В результате для (87) после интегрирования по ϕ и замены переменных интегрирования u = λr, x = cos θ находим:

 

ρ

+1

ρ0

 

32

 

Qxx =

 

Z0u8e−u du Z−1 (1 − x2)2[3π(1 − x2) − 2π]dx =

 

8! 2π

 

.

λ9

λ9

15 · 7

Аналогичные вычисления для Qyy дают: Qyy = Qxx, а для Qzz = −2Qxx при α 6= β Qαβ равен нулю.

Пример 1.7.2. Полусфера радиуса R заряжена с постоянной объемной плотностью ρ. Определить компоненты дипольного и тензора квадрупольного моментов.

Выбирая систему координат в центре полусферы, получим dx = dy = 0. Соответственно для dz:

 

R

π/2

 

3

 

dz = Z

ρzr2dr sin θ dθ dϕ = ρ Z0

r3 dr Z0

Z0

cos θ sin θ dθ =

 

QR,

8

где Q — полный заряд полусферы: Q = ρ23 πR3, Qxx — компонента тензора квадрупольного момента:

Z R Z π/2 Z

Qxx = ρ(3x2 − r2)r2dr sin θ dθ dϕ =

0 0 0

= ρ Z R r4dr Z π/2 Z 2π(3 sin2 θ cos2 ϕ − 1) sin θ dθ dϕ = 0.

00 0

Аналогично находим, что равны нулю и другие компоненты тензора Qαβ.

Пример 1.7.3 Цилиндр высоты H и радиуса R заряжен с объемной плотностью ρ = ρ0 = const. Определить Qxx элемент тензора квадрупольного момента такой системы зарядов.

По определению для Qxx в системе координат, связанной с осью цилиндра, находим:

Qxx = ρ0

 

R

 

H

 

[3ρ2 cos2 α

 

2+z2)]ρ dρ dz dα =

ρ0πHR4

1

 

4

H

2 .

 

0

 

0

 

 

 

 

 

 

 

 

 

Z

 

Z

 

Z

 

4

 

3

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

27

Пример 1.7.4. Определить поле системы точечных зарядов (см. рис.12) на больших расстояниях.

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

 

 

~

~

~

~

~

 

 

d = q(−ia) − q(ja) + q(ia) − q(ja) = 0.

 

 

Определим компоненты квадрупольного момента,

Рис. 12:

 

для примера:

 

 

 

 

 

 

 

 

 

Q11 = Qxx =

4

3xi2 − ri2

= q2a2

+ qa2 + q2a2 + qa2 = 6a2q.

=1 qu

iX

 

 

 

 

 

Аналогично вычисляются другие компоненты. В результате Q22 = −Q11, а остальные компоненты равны нулю. Таким образом, на больших расстояниях поле имеет вид:

ϕ =

1

(x2Qxx + y2Qyy) = −

3a2q

sin2 θ(1 − 2 cos2 ϕ),

2r5

r3

где r, θ, ϕ — переменные сферической системы координат.

Пример 1.7.5. Начало декартовой системы координат совпадает с центром тонкого кольца радиуса R. Найти поле на больших расстояниях от кольца с точностью до квадрупольных слагаемых, если кольцо заряжено с линейной плотностью τ.

Полный заряд кольца равен: q = 2πτR. Дипольный момент

~ Z ~ Z ~ ~

d = rρdV = (ix + jy)τR dα = 0.

0

Компоненты тензора квадрупольного момента:

Z Z

Qxx = [3x2 − R2]τR dα = [3R2 cos α − R2]τR dα = τR3π.

Аналогично Qyy = τR3π, Qzz = −τR32π и Qαβ = 0 при α 6= β. В результате

для скалярного потенциала находим:

 

 

 

 

ϕ(~r ) =

Q

+

1 3

 

Q

+

1

[x2Qxx + y2Qyy] =

r

2r5

Qαβxαxβ ' r

2r5

 

 

 

 

 

 

X

 

 

 

 

 

 

 

αβ=1

28

=Q + τR3π h1 − cos2 θi r 2r3

Пример 1.7.6. Найти поле

 

точечного

диполя

~

 

= const,

 

d

 

расположенного в начале координат.

 

 

 

 

 

 

~

~

 

 

 

 

 

 

 

 

 

 

В этом случае P =

dδ(~r ), и из решения уравнения (10) получим для

вектора Герца электрического типа:

~

~

 

 

 

 

 

Π = d/r. Отсюда:

 

 

 

 

~

 

 

 

 

~

 

~

 

E~ = 4π~dδ(~r ) + rot rot

d

=

~dδ(~r ) +

3~n(~nd) − d

.

 

 

 

r

 

 

 

 

r3

 

 

 

Задание на дом: решить задачи 2.19–2.24 на стр. 34.

1.8 Силы и энергия в электростатическом поле

 

 

 

 

 

 

 

 

Величина

 

 

 

Формула

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Энергия

единицы

объема

 

~ ~

; ε = Z

 

 

 

w и полная энергия ε

w =

ED

w dV.

(88)

 

 

 

электростатического поля

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Сила,

действующая

на

 

 

 

~

~

 

(89)

 

 

~

 

 

 

 

 

F = qE

 

 

 

 

 

 

 

 

 

 

 

 

 

 

точечный заряд E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Сила,

действующая

на

 

 

F~ = Z

ρE~ dV

(90)

 

распределенный в пространстве

 

 

 

 

 

 

 

 

 

заряд

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Сила, действующая на единицу

 

 

 

 

E2

 

(91)

 

поверхности

проводника

 

 

 

~p =

 

~n

 

 

 

 

 

(давление)

 

 

 

~n - нормаль к поверхности

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

29

 

Величина

 

 

 

Формула

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Выражение

для

силы

с

 

F= I

 

 

 

 

k

MklnkdS,

(92)

 

 

использованием

тензора

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

 

 

натяжений Максвелла

 

Mkl — тензор натяжений;

 

 

 

 

 

 

 

 

Mkl =

1

 

"EkEi

1

E2 δki#

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

~n — нормаль к dS.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Сила,

действующая

на

 

 

~

 

 

 

 

 

 

 

~

·

 

 

 

 

~

 

(93)

 

 

дипольный момент в поле

 

 

 

F = (d

grad )E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Знергия

заряда во

внешнем

 

 

 

Z

 

 

 

ρ(~r )ϕ(~r 0)dV

(94)

 

 

поле

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Энергия системы зарядов

 

 

 

 

1

 

 

 

 

 

 

 

qiqj

 

 

 

 

 

(95)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

=j ~ri

~rj

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

iX

 

 

 

|

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

|

 

 

 

 

 

 

 

Энергия системы проводников

 

 

 

 

 

 

 

1

Xi qiϕi

 

 

 

 

 

(96)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Энергия

диполя во

внешнем

 

 

 

 

 

 

 

 

 

 

~

~

 

 

 

 

 

 

(97)

 

 

поле

 

 

 

 

 

 

 

 

 

 

 

− (d

· E)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Энергия

 

квадрупольного

 

1

3

 

 

 

 

 

 

 

 

 

2

ϕ(0)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Qα,β

(98)

 

 

момента

во

внешнем

поле

со

 

6

α,β=1

∂x

α

∂xβ

 

 

скалярным потенциалом ϕ

 

 

 

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 1.8.1. Заряд Q равномерно распределен по поверхности

сферы радиуса R. Найти абсолютную величину силы, разрывающей

сферу на две равные половины. Решить задачу двумя методами с использованием формул (91) и (92).

Данная система зарядов создает поле, равное нулю для r < R, и поле, совпадающее с полем точечного заряда Q/r2 для r > R.

1 способ решения. Выделим произвольную половину сферы. На элемент поверхности dS действует сила E2~n dS/8π. Результирующая сила равна:

 

E2

π/2

E2

F~ = Z

 

~n dS =

Z0

Z0

 

~nR2 sin θ dθ dϕ.

Вычислим проекции этой силы на декартовы оси с учетом: nx = sin θ cos ϕ,

30

ny = sin θ sin ϕ. Например, для Fz имеем:

π/2

E2

Q2

fz = Z0

Z0

 

cos θR2 sin θ dθ dϕ =

 

.

8R2

Проекции Fx = Fy = 0.

2 способ решения. В соответствии с (92) имеем:

I

Fz = [Mxznx + Myzny + Mzznz] dS.

Поверхность интегрирования включает половину сферы и основание. Однако интеграл по основанию равен нулю в силу Mik=0 для точек r < R. Таким образом,

Fz =

π

 

ExEz

sin θ cos ϕ +

EyEz

sin θ sin ϕ +

Ez2 21 E2

Z0

Z0

 

 

 

 

 

 

 

 

cos θ sin θ dθ dϕ.

 

 

(99)

По определению на сфере Ex =

Q sin θ cos ϕ/R2, Ey

= Q sin θ sin ϕ/R2,

Ez = Q cos θ/R2. Подставляя

Ek в (99), находим

после выполнения

интегрирования: Fz = Q2/8R2, как и при первом способе решения.

Пример 1.8.2. Шар радиуса R заряжен с объемной плотностью ρ = const. Найти силу, разрывающую шар на две равные половины. Решить задачу двумя методами.

На основании решения примера 1.2.2 известно, что поле такой системы

зарядов есть

 

 

 

 

 

 

 

 

~

4

πρ~r,

~

Q~r

 

(100)

E =

3

r < R; E =

r2

 

r

, r > R

1 способ решения. В соответствии с (90) сила, действующая на единицу

объема равна

~

ρE. Следовательно, для определения суммарной силы,

действующей на полусферу, надо проинтегрировать данное выражение по объему полусферы. При интегрировании получим, что Fx и Fy равны нулю. Соответственно для Fz находим:

R

π/2

 

3 Q2

 

Fz = Z0

Z0

Z0

ρE(r) cos θr2dr sin θ dθ dϕ =

 

 

 

,

16

R2

где Q — полный заряд шара.

2 способ решения. На основании (92) необходимо интегрировать по замкнутой поверхности, которая в данном случае состоит из основания и поверхности полусферы. То есть интеграл разбивается на сумму двух

31

интегралов. При интегрировании по основанию для нормали к поверхности

имеем nx

= 0, ny = 0, nz

= −1. Интегрирование удобно выполнять в

цилиндрической системе координат, в которой dS = ρ dρ dα,

 

1

 

R

1

Fz(1)

= Z MzznzdS = −

 

Z0

Z0

"Ez2

 

(Ex2 + Ey2 + Ez2)# ρ dρ dα.

2

С учетом (100) и соотношений Ex = E0 cos α, Ey = E0 sin α, Ez = 0, E0 = 43 πρ2, θ = π/2 находим после интегрирования: Ez(1) = Q2/16R2.

При интегрировании выражения (92) по поверхности сферы необходимо подставить nx = sin θ cos ϕ, ny = sin θ sin ϕ, nz = cos θ. В результате

 

π/2

Q

 

Fz(2) =

Z

Z [Mxz sin θ cos ϕ + Myz sin θ sin ϕ + Mzz cos θ]2R sin θ dθ dϕ =

,

 

8R2

 

0

0

 

 

где использованы выражения для тензора натяжений Максвелла:

Mxz =

E12

sin θ cos ϕ cos θ; Myz =

E12

sin θ sin ϕ cos θ; Mzz = "cos2

θ

1

#

E1

 

 

 

 

 

2

 

 

 

 

 

 

Суммарная сила: Fz = F (1)

+ F (2)

= 3Q2/16R2.

 

 

 

 

 

 

 

 

z

z

 

 

 

 

 

 

 

 

Комментарий: может показаться, что метод, основанный на использовании тензора натяжений Максвелла, весьма громоздок и неудобен. Однако преимущество этого метода проявляется там, где нет простой симметрии.

Пример 1.7.3 Шар радиуса R однородно заряжен с объемной плотностью ρ = const. Определить энергию поля, создаваемого зарядами шара.

Напряженность поля определяется выражением (100). В соответствии с

(88) полная энергия поля равна:

 

 

 

 

 

 

 

 

ε =

1

R π 2π

 

Qr

 

2 r2dr sin θdθdϕ +

∞ π 2π

 

Q

 

2 r2dr sin θdθdϕ

=

3Q

 

 

 

 

 

 

 

 

 

 

" r3

#

Z Z Z

"r2

#

 

5R

 

Z Z Z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 0 0

 

 

 

 

R 0 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 1.8.4. Потенциал внешнего электростатического поля имеет вид ϕ(x, y) = α exp[−β(x2 − y2)] и мало меняется по размеру квадруполя (см. рис.12 на стр 27). Найти энергию квадруполя во внешнем поле.

В соответствии с (98) и решением, приведенном в примере (1.7.4) , находим:

ε = 61

 

 

2ϕ

=

1

Qxx

2ϕ(0)

2ϕ(0)

 

= 4αβa2q.

X

Qαβ ∂x ∂x(0)

6

∂x2

+ Qyy ∂y2

 

 

 

 

α β

 

 

 

 

 

 

 

 

 

 

 

 

 

αβ