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

иНорматика el_polya

.pdf
Скачиваний:
20
Добавлен:
01.05.2015
Размер:
4.32 Mб
Скачать

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

 

 

 

 

E

.

 

rot H = jЭ + а

 

(1.19)

 

 

 

 

 

 

t

 

а для поля от тока, распределенного по некоторой площади S, ограниченной кривой L используют интегральную форму закона.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

 

 

 

 

 

rot H dS

H d l IЭ а

 

dS ,

(1.20)

 

 

 

S

 

 

 

 

L

 

 

 

 

 

S

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Эти соотношения выглядят несколько сложнее выражений для расчета электрического поля (1.7) и (1.8). Это связано с двумя причинами.

Во-первых, в левой части уравнения (1.19) стоит операция нахождения ротора, а эта операция сложнее операции расчета дивергенции (см.1.7). Магнитное поле направлено перпендикулярно к току и направлению от точки, в которой проводится расчет. Плотность тока не является источником расходящегося из этой точки магнитного поля, а определяет вращательную его составляющую. Эта операция ближе к векторному произведению, а не скалярному. Действительно, если в

прямоугольной системе координат ввести фиктивный вектор с про-

екциями

 

,

 

,

 

, то левую часть уравнения (1.7) можно предста-

x

y

z

 

 

 

 

 

 

 

 

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

а

 

 

 

E , а левую часть уравнения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.19) в виде векторного произведения

 

 

H .

Во-вторых, в уравнении (1.19) в правой части кроме плотности то-

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

jЭ

 

=

а

E

.

(1.21)

см

 

 

 

 

 

t

 

Из структуры выражения (1.20) следует считать эту величину еще одним видом тока. Ее называют током смещения. Таким, образом, магнитное поле создается ток проводимости (1.12) или током смеще-

ния (1.21).

Процесс возникновения магнитного поля от второго слагаемого в правой части выражения (1.19) можно рассмотреть и по-другому. Магнитное поле возникает тогда, когда есть ток или изменяющееся электрическое поле. Изменяющееся электрическое поле создает перемен-

11

ное магнитное поле. Такой процесс можно назвать индукцией. Впрочем, электромагнитной индукцией принято называть обратный процесс – возникновение электрического поля при изменении магнитного.

Выражения для расчета магнитного поля (1.19) и (1.20) можно записать для переменных токов, изменяющихся с определенной круговой частотой , воспользовавшись методом комплексных амплитуд.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.22)

 

rotH

 

jЭ

i

 

n E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H d l

IЭ

i

 

n E dS

(1.23)

L

 

 

 

 

 

 

 

 

 

 

 

 

S

Принцип перестановочной двойственности

Мы получили два уравнения Максвелла, позволяющие рассчитать электрическое поле, создаваемое электрическим зарядом и магнитное поле, создаваемое электрическими токами проводимости и смещения. Существует еще два уравнения Максвелла, второе и четвертое. Чтобы записать эти уравнения предполагают, что электрические и магнитные явления однотипны. Наряду с электрическими зарядами и токами вводятся магнитные заряды qm и плотность магнитных токов проводимо-

 

 

 

 

 

сти jМ и смещения a

H

. Магнитные заряды создают расходящееся

 

t

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

E

H ;

D

– B ;

jэ – jм ;

 

n

 

n ;

э

м; σэ– σм . (1.24)

 

 

 

 

 

 

 

 

 

 

 

 

 

В этом и состоит принцип перестановочной двойственности, анало-

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

Закон непрерывности магнитных силовых линий (второе уравнение Максвелла)

Применим принцип перестановочной двойственности к первому уравнению Максвелла и запишем закон непрерывности магнитных силовых линий последовательно в дифференциальной (1.25), инте-

12

гральной (1.26) форме и для комплексной амплитуды напряженности полей (1.27), (1.28).

 

 

 

 

 

 

 

 

М .

(1.25)

div

а H

 

 

 

 

 

 

 

 

 

qм

(1.26)

S

а H dS

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

м .

 

div

n H

(1.27)

 

 

 

 

 

 

 

q м .

(1.28)

 

 

 

 

 

 

S

n H dS

 

 

 

 

 

 

 

 

 

 

Здесь введена комплексная магнитная проницаемость, которая получается из (1.11)

 

=

a +

м

=

a - i

м

(1.29)

n

i

 

 

 

 

 

 

 

Закон электромагнитной индукции ( четвертое уравнение Максвелла)

Применим принцип перестановочной двойственности к третьему уравнению Максвелла и запишем закон электромагнитной индукции последовательно в дифференциальной (1.30), интегральной (1.31) форме и для комплексной амплитуды напряженности полей (1.32), (1.33).

 

 

 

 

 

H

.

(1.30)

rotE j

0

 

 

 

 

 

m

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

H

 

 

 

 

 

Edl Iм

0

 

 

dS .

(1.31)

L

S

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rotE

jм

i

n H .

(1.32)

 

 

 

Iм

i

 

 

 

 

 

 

 

(1.33)

 

 

 

 

 

 

 

 

L

Edl

 

 

n H dS .

 

 

 

 

 

 

 

 

 

 

 

 

S

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

13

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица1.1

 

 

 

 

 

 

 

Уравнения Максвелла

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для функций времени

 

 

 

 

 

 

Для комплексных амплитуд

 

 

Вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Дифференциальная

 

Интегральная

 

 

 

Дифференциаль-

 

Интегральная

 

Название

форма

 

 

 

 

форма

 

 

 

 

ная форма

 

 

 

 

форма

 

 

(номер)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Закон Гаусса

 

а E

 

 

 

 

а E dS

q э

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q Э

 

(1)

div

Э

 

 

S

 

 

 

div

n

E

 

Э

 

 

S

n E dS

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Закон непре-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

рывности маг.

 

а H

 

 

 

 

а H dS

qм

 

 

 

 

 

 

 

 

 

 

 

 

 

q м

 

силовых. линий

div

м

 

 

S

 

 

div

n

H

 

м

 

 

S

n H dS

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Закон полно-

 

 

 

E

 

 

 

 

E

 

 

 

 

 

 

 

 

 

 

 

IЭ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L H dl

i

n E dS

го тока (3)

rot H =

jЭ +

а

 

 

H dl

IЭ

а

 

 

dS

rotH

jЭ

i

n E

 

 

 

 

 

 

S

 

 

t

L

S

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Закон элек-

 

 

 

 

 

 

 

 

 

H

 

 

 

 

 

 

 

 

 

I

 

 

 

 

тромагнит-

 

 

 

H

 

L Edl

 

 

 

 

 

 

 

Edl

 

 

i

 

H dS

rotE

j

0

 

Iм

 

а S

 

t

dS

rotE

jм

i

n H

L

 

 

 

M

 

n

 

ной индук-

 

m

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ции (4)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для электрической проницаемости:

0

= 10-9/36

(Ф/м) ;

 

 

 

а =

0 ;

 

n =

a +

Э

=

a - i

 

Э .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

Для магнитной проницаемости:

 

 

0 = 4 10-7 Гн/м;

 

 

 

а =

0 ;

 

n

=

a +

м

=

 

a - i

м .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

14

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Лабораторная работа 2. Расчет электрического поля от дискретного набора электрических зарядов

Цель работы: Научиться рассчитывать напряженность электрического поля от произвольного набора электрических зарядов в произвольной точке трехмерного пространства.

Используемое программное обеспечение: MathCAD, версия не ниже MathCAD13.

Лабораторное задание

1.Освоить методику составления программы для расчета электрического поля от набора точечных источников.

2.Составить программу расчета электрического поля в наборе точек трехмерного пространства, для которого = 1. Поле созда-

ется набором точечных зарядов (до 10 штук, в задании даются значения в 4 точках, см. табл.1.2). Заряд, величина которого задана в кулонах*10-9, располагается в точке с координатами (х1, y1, z1).

Рассчитать электрическое поле в 10 точках, равномерно распределенных на заданном отрезке прямой, начало и конец которой заданы в табл.1.3. Рассчитывается модуль вектора электрического поля, его проекции на координатные оси. В отчете привести график изменения электрического поля и его проекций на координатные оси.

Методические указания к выполнению лабораторного задания

Лабораторное задание выполняется по вариантам. В таблицах 1.2 и 1.3 приведены исходные данные. В каждом варианте задается 4 заряда. Для каждого заряда указано его положение – координаты в прямоугольной системе координат (x1, y1, z1) – и величина точечного заряда (q) в

кулонах *10-9. (За величину точечного заряда принимается действительная часть приведенного в таблице значения) . Столбцы x2 и y2 в

этой лабораторной работе не используются.

В таблице 1.3 заданы начало и конец отрезка прямой, на котором лежат точки, в которых рассчитывается электрическое поле. Расчет

15

проводится в 11 точках. Эти точки должны равномерно распределяться по прямой.

Таблица 1.2

.

Варианты задания, заряды, токи .

16

Продолжение таблицы 1.2 .

17

 

 

 

 

 

 

Таблица 1.3.

 

 

Координаты отрезка прямой .

 

 

 

Вариант

 

Начало

 

 

Конец

 

 

 

 

 

 

 

 

 

x1

y1

z1

x2

y2

z2

 

 

 

1

12

14

30

16

18

40

 

2

8

14

30

12

18

40

 

3

4

14

30

8

18

40

 

4

0

14

30

4

18

40

 

5

-4

14

30

0

18

40

 

6

-8

14

30

-4

18

40

 

7

-12

14

30

-8

18

40

 

8

12

14

30

12

18

40

 

9

8

14

30

12

18

40

 

 

 

 

 

 

 

 

 

10

4

14

30

8

18

40

 

11

0

14

30

4

18

40

 

 

 

 

 

 

 

 

 

12

-4

14

30

0

18

40

 

 

 

 

 

 

 

 

 

13

-8

14

30

-4

18

40

 

14

-12

14

30

-8

18

40

 

 

 

 

 

 

 

 

 

15

12

14

30

16

18

40

 

16

8

14

30

12

18

40

 

17

4

14

30

8

18

40

 

18

0

14

30

4

18

40

 

19

-4

14

30

0

18

40

 

20

-8

14

30

-4

18

40

 

21

-12

14

30

-8

18

40

 

22

12

14

30

16

18

40

 

23

8

14

30

12

18

40

 

24

4

14

30

8

18

40

 

 

 

 

 

 

 

 

 

25

0

14

30

4

18

40

 

26

-4

14

30

0

18

40

 

 

 

 

 

 

 

 

 

27

-8

14

30

-4

18

40

 

 

 

 

 

 

 

 

 

28

-12

14

30

-8

18

40

 

 

 

 

 

 

 

 

 

 

 

 

18

 

 

 

 

Выполняя лабораторную работу, следует сначала в программе MathCAD составить файл, с помощью которого можно рассчитать поле в произвольной точке, если число зарядов будет достигать 10. Отладить файл на тестовом примере (вариант 28), а затем, убедившись, что расчет проводится без ошибок, ввести исходные данные по варианту и получить результат.

Обозначение переменных:

xp, yp, zp – координаты точки, в которой рассчитывается поле;

r – двумерный массив расстояний между точками, в которые помещен заряд и точками, в которых рассчитывается поле;

fix, fiy,fiz – углы, которые составляет радиус вектор r c осями координат x, y, z;

i – номер точки, в которую помещается заряд;

xn, yn, zn, q – координаты и величина i –го заряда;

dx, dy, dz – проекция расстояния между точкой, в которой находится заряд и точкой, для которой рассчитывается электрическое поле на оси x, y, z соответственно,

E, EX, EY,EZ – электрическое поле от одного заряда и его проекции на координатные оси,

EP, EPX, EPY, EPZ – модуль суммарного вектора электрического поля от всех зарядов и его проекций на координатные оси.

sp – расстояние между точками, в которых рассчитывается электрическое поле. Это расстояние рассчитывается в плоскости хoy

10 9

0

 

Ф/м – диэлектрическая проницаемость вакуума;

36

 

 

 

– относительная диэлектрическая проницаемость;

а = * 0 – абсолютная диэлектрическая проницаемость.

Для расчета воспользуемся выражением (1.2) которое позволяет найти поле точечного заряда. В каждой точке рассчитаем поле от каждого заряда и сложим их.

Алгоритм расчета поля.

1.Введите исходные данные:

a. диэлектрическую проницаемость

19

 

0 :

10

9

:= 1;

 

 

 

 

 

36

 

 

 

 

 

b.

номера точек, в которых располагаются заряды, i и в которых

 

надо рассчитать электрическое поле, k

c.

i := 0..9; k := 0..9;

 

 

 

 

d.координаты и величины зарядов. Переменные xn, yn, zn, q вводятся, как одностолбцовые матрицы. Столбец содержит 10 элементов. Первые 4 элемента всех столбцов вводятся в соответствии с вариантом. Остальным элементам присваивается значение 0. Можно после отладки программы некоторые или все нули заменить значениями координат и величинами зарядов, если вы решили рассчитать электрическое поле для большего их числа.

e.координаты начальной и конечной точки прямой, на которой находятся точки для расчета электрического поля x1, y1, z1, x2, y2, z2.

2.Введите формулы для расчета

a.Координат точек, в которых рассчитывается электрическое поле и диэлектрической проницаемости:

b.

c. проекции на координатные оси расстояния между точками, в которых находятся заряды и точка, для которой рассчитывается электрическое поле, само расстояние и углы с координатными осями. Будет рассчитываться двумерные массивы значений, поэтому формулы должны выглядеть так:

Индексы k, i означают, что рассчитывается величина (например, расстояние rk,i) в k – ой точке от i-ой точки;

20