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

05 Методы решения задач электростатики

.pdf
Скачиваний:
112
Добавлен:
07.02.2015
Размер:
287.25 Кб
Скачать

Лекция № 5

Методы решения задач электростатики.

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

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

1.Решение уравнения Пуассона для бесконечного равномерно заряженного круглого цилиндра

Вкачестве простого примера решения уравнения Пуассона найдем потенциал бесконечного цилиндра радиусом a, равномерно заряженного с объемной плотностью

=const

Направим ось Z вдоль оси цилиндра. Из-за аксиальной симметрии распределения заряда потенциал также аксиально симметричен, т.е. = r . Поэтому удобно использовать цилиндрическую систему координат, аксиальный угол в которой обозначим . Оператор Лапласа имеет вид

2

 

2

 

1

1 2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

(1)

r2

 

 

r2

 

2

z2

 

 

 

r r

 

 

 

 

Так как потенциал зависит только от r, то в данном случае выражение существенно упрощается:

 

 

 

 

 

 

 

 

2

 

 

 

 

2

1

 

1 d

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

r2

 

r r

 

 

 

и уравнение Пуассона приобретает вид

 

 

 

 

 

 

 

 

r dr

dr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 d

 

d

1

 

 

 

 

 

 

 

при 0<r<a;

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dr

 

0

 

 

 

r dr

 

 

 

 

 

 

 

 

 

 

 

 

 

1 d

 

d

2

 

 

 

 

при r>a.

 

 

 

 

 

 

 

 

 

r

 

 

 

0

 

 

 

 

 

 

 

 

 

dr

 

 

 

 

r dr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Общие решения этих уравнений находят непосредственным интегрированием:

1 = r2 + A1lnr + B1, 4 0

(2)

(3)

2 = A2 lnr + B2.

(4)

Поскольку потенциал во всех точках должен быть конечным, а lnr при r то в решении (4) необходимо положить А1=0. Удобно нормировать потенциал условием 1(0) = 0. При этом

В1=0 (5)

Методы решения задач электростатики

2

________________________________________________________________________________

Условия непрерывности потенциала и его производной ( '1 =

 

 

r

,

 

'2

= A2/r) при r = a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 0

 

 

 

 

дают два алгебраических уравнения относительно неизвестных А2

и В2:

 

 

 

 

 

A2lna + B2 =

 

a2 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A2

 

 

 

 

 

 

 

 

 

a,

 

 

 

 

 

 

(6)

 

 

 

 

 

 

 

 

 

a

2 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A2

 

a2 ,

 

B2 =

 

 

a2 -A2lna =

 

 

a2

+

 

 

a2 lna

 

 

4 0

4 0

2 0

Тогда

2 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1(r)

=

 

 

 

 

 

r2 при 0<r<a,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2(r) =

 

 

 

a2

 

 

ln

a

 

 

 

 

 

a2 при r>a.

 

 

 

(7)

 

 

 

2 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

4 0

 

 

 

 

 

 

 

Напряженность электрического поля выражается формулами:

 

 

 

 

 

 

 

 

 

Er

 

1

 

 

 

 

 

r

 

 

 

при 0<r<a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

r

 

 

 

2 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Er2

 

2

 

 

 

 

 

 

 

a2

 

 

 

при

r>a

 

 

 

(8)

 

 

 

 

 

 

2 0

r

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Напряженность поля внутри цилиндра растет прямо пропорционально радиусу, а вне цилиндра - убывает обратно пропорционально расстоянию от оси. Этот же результат может быть просто получен с помощью теоремы Гаусса.

2. Формулы Грина

Рассмотрим эффективный метод решения уравнения Пуассона при заданных граничных условиях, сводящийся к замене уравнения Пуассона интегральным уравнением с помощью формул Грина. Исходной при получении формул Грина является формула Остроградского-Гаусса

divAdV AdS.

(9)

V

S

 

Пусть имеются две скалярные функции и . Возьмем в качестве вектора А в формуле (9) вектор grad . Тогда

div( grad ) = divgrad +grad grad = 2 + (10)

С учетом (10) выражение (9) принимает вид

( 2 )dV

 

 

dS,

(11)

 

V

S

 

n

 

где ngrad = , если n - единичный вектор, направленный по нормали к поверхности.

n

Положительное направление выбрано по внешней нормали к поверхности, ограничивающей объем V. Формула (11) называется первой формулой Грина.

Если из (11) вычесть почленно аналогичную формулу, в которой поменять местами функции и , то получится вторая формула Грина

( 2 2 )dV

(

 

 

 

)dS.

(12)

n

 

V

S

 

 

n

 

Интегрирование уравнения для потенциала. Для того, чтобы формуле (12) придать вид, удобный для использования, возьмем в качестве функцию

Методы решения задач электростатики

3

________________________________________________________________________________

r) = 1/ r-r' = 1/R,

где r - радиус-вектор точки наблюдения, r' - переменная интегрирования. Примем во внимание формулу

 

2

1

 

 

4 (r-r')

 

r-r'

 

 

 

 

 

 

 

 

 

где (r-r') - функция Дирака. Правильность этой формулы строго доказывается в математической физике. Однако проверка ее справедливости осуществляется очень просто. Функция в формуле (13) является решением уравнения Пуассона, если заряды расположены в конечной области пространства

(r) =

1

V

(r')dV'

.

(13)

4 0

 

r-r'

 

 

 

 

 

 

 

 

Применим к обеим частям формулы (13) оператор . Тогда

2 (r)

(r)

 

1

(r') 2

 

1

dV'

(14)

 

 

 

r r'

 

0

4 0

 

 

Это равенство тождественно справедливо для всех (r), что возможно лишь при

2

 

 

1

 

 

4 (r-r'),

 

 

r-r'

 

 

 

 

 

 

 

 

 

 

когда равенство (14) сводится к определению функции:

(r) (r') (r-r')dV'.

В качестве в формуле (12) выберем потенциал, удовлетворяющий уравнению Пуассона.1 Тогда (12) принимает вид

[ 4 (r') (r-r')

1

 

(r')

]dV' [ (r')

 

(

1

)

1

 

 

]dS'

(15)

0 R

 

 

 

 

V

S

n' R

R n'

 

где R= r-r' - расстояние между точкой наблюдения и текущей точкой интегрирования.

Выполняя в (15) интегрирование члена с

- функцией, окончательно2 получим

 

(r)

1

 

(r')

 

1

 

1

 

 

 

1

 

 

 

 

 

dV'

 

 

 

 

 

 

 

(

 

)

dS'.

(16)

 

R

4

 

 

 

 

 

4 0 V

 

S

R n'

 

n' R

 

 

 

Эта формула выведена для точек, расположенных внутри объема V. Для точек вне объема правая часть (16) равна нулю.

Если ограничивающая объем поверхность S удаляется на бесконечность и электрическое поле при этом убывает быстрее чем 1/R, то поверхностный интеграл в выражении (16) стремится к нулю, и оно переходит в формулу (13).

На первый взгляд может показаться, что формула (15) дает решение уравнения Пуассона с заданными граничными условиями, поскольку значения функций и / n входят в уравнение только в точках границы объема. Однако такое заключение является неправильным. Дело в том, что на границе области для уравнения Пуассона нельзя задать одновременно и значения функции и ее производной по нормали / n, так как в этом случае задача оказывается переопределенной и, вообще говоря, не имеет решения. Поэтому формула (15) является интегральным уравнением.

1

Вместо = 1/R, = -

4 (r r'), = - / 0

 

 

 

 

2

4 (r)

1

 

(r')

dV' [ (r')

 

(

1

)

1

 

 

]dS'

 

0

R

n'

 

 

 

 

 

 

S

 

R

R n'

Методы решения задач электростатики

4

________________________________________________________________________________

3. Метод изображений

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

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

Поле заряженной оси, расположенной вблизи проводящей плоскости. Заряженная ось

( - заряд на единицу длины) расположена в диэлектрике параллельно поверхности проводящей среды.

h

Проводящая среда

h

m -

Проводящей средой может быть какая-либо металлическая стенка или, например, земля. Требуется определить характер поля в верхней полуплоскости (диэлектрике).

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

Поместим в точку m фиктивный заряд обратного знака (- ). Расстояние h от точки m до плоскости раздела сред такое же, как и расстояние от действительного заряда до плоскости раздела. В этом смысле осуществлено зеркальное изображение.

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

Найдем поле заряженной оси. Применяя теорему Гаусса, найдем напряженность электрического поля

 

 

 

EdS

q

E2 rl

q

 

 

 

 

 

;

 

 

 

 

 

 

 

0

0

 

или окончательно E

q

 

 

, где -

линейная плотность заряда. Для нахождения

2 r 0l

2 r 0

 

 

 

 

 

 

 

 

 

 

 

 

 

потенциала необходимо проинтегрировать значение напряженности

Edl Edr

 

 

 

dr

 

 

lnr B,

2 0

 

2 0

 

 

 

 

 

 

 

 

r

 

где В – постоянная интегрирования.

Легко убедиться в том, что потенциал удовлетворяет уравнению Лапласа, для этого подставим полученное выражение в выражение для оператора Лапласа в цилиндрических координатах (зависимость от полярного угла и от координаты z отсутствует)

Методы решения задач электростатики

5

________________________________________________________________________________

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

=0.

 

 

 

r

 

 

 

 

 

 

 

(

 

lnr)

 

 

 

 

 

 

 

r

r

2 0

 

r r

r

 

r r

 

 

 

 

Так как потенциал каждой из осей удовлетворяет уравнению Лапласа и одновременно удовлетворено граничное условие, то на основании теоремы единственности полученное решение является единственным.

4. Теорема взаимности

При расчетах электростатических полей часто используется следующее утверждение, носящее название теоремы взаимности. Потенциал поля, создаваемый в точке r2 зарядом q, находящимся в точке r1, равен потенциалу, создаваемому в точке r1 тем же зарядом, помещенным в точке r2.

В случае одного точечного заряда эта теорема имеет тривиальный характер

(r1) =

1

 

 

 

 

q

 

=

1

 

 

 

 

q

 

= (r2).

4

 

 

 

r

 

 

4

 

 

 

r

r

 

 

0

 

 

r

2

 

 

0

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

2

1

 

 

Теорема взаимности допускает обобщение. Пусть в точке r1 находится заряд q1, а в точке r2 заряд q2. Тогда имеем, очевидно:

q1 2 (r1) = q2 1 (r2).

Действительно, потенциал, создаваемый в точке r1 зарядом q1, равен

(r2) =

1

 

 

 

 

q1

 

 

,

4

0

 

 

r

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

 

 

тогда как

 

 

 

 

q2

 

 

 

(r1) =

1

 

 

 

 

 

 

.

4

0

 

 

r

r

 

 

 

 

 

 

 

 

 

 

 

 

 

2

1

 

 

 

Пусть, наконец, имеется произвольное множество проводников, имеющих заряды qi при потенциалах i и заряды q'i при потенциалах 'i. Тогда имеет место равенство

qi i' = q'i i . i i

Доказательство проводится путем суммирования потенциалов точечных зарядов. Теорема взаимности весьма полезна для вычисления зарядов, индуцированных на

проводниках.

Пусть, например, к заземленной проводящей сфере радиуса a приближается заряд q и помещается в точке, отстоящей на расстоянии r от начала координат, находящегося в центре сферы. На ее поверхности индуцируется заряд (-q'), причем на основании теоремы взаимности

(a)q' = - q (r),

или

q' q (r) qa .(a) r

5. Метод конформных преобразований

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

Электростатическое поле в пустоте удовлетворяет двум уравнениям: rotE = 0 и divЕ =0,

Методы решения задач электростатики

6

________________________________________________________________________________

такое поле мы называли лапласовым. Первое из них позволяет ввести потенциал поля согласно уравнению Е = - grad . Второе же уравнение показывает, что наряду со скалярным потенциалом можно ввести и векторный потенциал поля А согласно E = rotA.

В плоском случае вектор Е лежит в плоскости ху и зависит только от этих двух координат. Соответственно, вектор А можно выбрать так, чтобы он был везде направлен

перпендикулярно к плоскости ху и имел только составляющую АZ.

Тогда компоненты

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

можно выразить в виде производных от или

А следующим образом EX

 

 

 

rotXA

AZ

 

AY

 

A

 

 

 

 

 

 

 

y

 

 

 

x

 

 

 

 

y

 

z

 

 

 

EX

 

 

A

,

EY

 

 

 

A

.

(17)

 

x

 

y

 

 

 

 

 

 

y

 

 

 

 

 

 

x

 

Но такие соотношения между производными функций и А с математической точки зрения совпадают с условиями Коши – Римана, выражающими тот факт, что комплексное выражение

w iA

(18)

 

 

является аналитической функцией комплексного аргумента

z x iy. Это значит, что

функция w(z) имеет в каждой точке определенную производную, не зависящую от

направления, в котором оно берется. Так, дифференцируя в направлении оси X, найдем, что

 

dw

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

,

 

 

dz

 

 

 

 

 

 

 

 

x

x

 

или

 

 

 

 

 

 

 

 

 

 

dw

 

EX iEY .

(19)

 

 

 

 

 

 

 

 

 

dz

 

 

 

 

 

 

 

 

 

 

Функция называется wкомплексным потенциалом.

Вектор напряженности поля Е коллинеарен вектору dr, который направлен по касательной к силовой линии. Векторы коллинеарны, если их проекции пропорциональны.

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

dx dy . EX EY

Выражая ЕХ и ЕУ через производные от А (EX A иEY A ), перепишем это уравнение

y x

в виде dx A dy A = 0, откуда А(х,у) = const. Таким образом, линии постоянных значений

x y

мнимой части функции w(z) представляют собой силовые линии поля. Линии же

постоянных значений ее вещественной части являются эквипотенциальными линиями. Взаимная ортогональность этих двух семейств линий обеспечивается соотношениями (17), согласно которым (ЕХЕуХЕУ=0)

 

 

 

A

 

 

 

A

0.

 

 

 

 

 

 

 

 

 

 

x x

y y

 

 

Как вещественная, так и мнимая части аналитической функции w(z) удовлетворяют

 

 

 

 

 

 

 

 

 

 

 

уравнению Лапласа. Поэтому с тем же

успехом можно принять Imw(z) в качестве

 

 

 

 

 

 

 

 

 

 

 

Методы решения задач электростатики

7

________________________________________________________________________________

потенциала поля.

Соответственно силовые линии будут

тогда даваться уравнениями

Rew(z)=const. Вместо (18) при этом будем иметь w A i .

 

 

 

 

Простейшим

примером комплексного потенциала

является потенциал поля

заряженной прямой нити (совпадающей с осью z). Напряженность этого поля дается формулами

Er

 

,

E

0,

2 0r

где r, полярные координаты в плоскости ху, а - заряд единицы длины нити. Соответствующий комплексный потенциал

w

 

lnz

 

lnrei

 

(lnr i ) .

 

 

 

 

2 0

 

2 0

 

2 0

 

 

 

 

 

Напряженность поля находится как модуль производной комплексного потенциала

E

dw

 

d

(

 

 

lnz)

 

 

 

 

1

 

 

 

 

1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dz

 

dz

2

0

 

 

 

2 0 z

 

2 0 r

 

 

 

 

 

 

Если же заряженная нить проходит не через начало координат, а через точку (хо, уо ), то комплексный потенциал

w ln(z - z0 ), 2 0

где z0 x0 iy0 .

С математической точки зрения функциональное соотношение w=w(z)

осуществляет конформное отображение плоскости комплексного переменного z на

плоскость комплексного переменного w. Пусть С есть контур сечения проводника в

плоскости ху, а 0 – потенциал этого проводника. Тогда задача об определении поля,

создаваемого этим проводником, сводится к нахождению такой функции w(z), которая

отображала бы контур С в плоскости z на линию

w= 0, параллельную оси ординат в

 

 

плоскости w; тогда вещественная часть Rew дает потенциал рассматриваемого поля. Если

 

 

же функция w(z) отображает контур С на линию параллельную оси абсцисс, то потенциал

дается формулой Imw.

Иногда поступают по другому – находят, какому полю соответствует заданная аналитическая функция (комплексный потенциал). Ценность метода невелика, так как лишь в ряде случаев можно угадать комплексный потенциал. Необходимо отметить, что хотя техника комплексных переменных часто оказывается очень мощной, она ограничена только двумерными задачами; к тому же это все-таки косвенный метод.