иНорматика el_polya
.pdfщью закона полного тока, который еще называют третьим уравнением Максвелла. Для точки используют дифференциальную форму закона:
|
|
|
|
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