
Электромагнитные поля и волны.-4
.pdf
164
данной точке. Величина вектора (интенсивность векторного поля) обычно характеризуют густотой силовых линий.
Согласно определению силовой линии, элемент касательной к ней параллелен вектору a в каждой точке, У параллельных векторов проекции пропорциональны. Это приводит к следующим уравнениям для силовых линий
dx |
= |
dy |
= |
dz |
, |
dr |
= |
rdα |
= |
dz |
, |
dr |
= |
rdθ |
= |
r sin(θ )dα |
(П1.13) |
|
|
|
|
|
|
|
|
|
|||||||||
ax a y az |
ar aα az |
ar aθ |
|
aα |
|
Каждое из этих уравнений распадается на систему из двух дифференциальных уравнений, интегрирование которых приводит к уравнению силовой линии.
Пример П1. 1. Найти векторную линию магнитного поля проводника с
r I r
током H = 2πr α0 , здесь r- расстояние от оси провода до точки М, I – ток в проводе. Выразим в цилиндрической системе координат орт α0 через орты по
|
|
|
|
|
|
r |
, r ] |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
r |
r |
r |
|
|
|
[z0 |
|
|
r r |
r |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
осям z и r α0 |
= [z0 |
, r0 ] |
= |
|
|
|
|
. |
Вектор r = x0 x + y0 y , |
х и у- координаты точки |
||||||||||||||||
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
r |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
наблюдения. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
r |
r |
r |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
r |
]= |
|
x 0 |
y0 |
z 0 |
|
|
r |
|
r |
|||
|
|
|
|
|
|
|
|
|
|
|
|
r |
|
0 |
0 |
1 |
|
|
|
|||||||
Раскрывая векторное произведение, получим [z0 |
, r |
|
= − yx 0 |
+ xy0 . |
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
y |
0 |
|
|
|
|
|
|
|
|
|
r |
|
|
|
I |
r |
I |
r |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
Следовательно, H = − |
|
|
|
|
yx 0 + |
|
xy 0 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
2πr |
|
|
2πr |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Дифференциальное уравнение векторных линий запишется как: |
− |
dx |
= |
dy |
= |
dz |
, |
|||||||||||||||||||
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y |
x 0 |
|||||
|
|
|
|
|
|
откуда |
x 2 |
+ y 2 = R 2 , и |
|
z = c . |
|
|
|
|
|
|
|
|
|
|
Т.е. векторные линии в данном случае являются окружностями с центром на оси z. Конечно, эту задачу можно решить и в цилиндрической системе координат. Предлагается это сделать самостоятельно.
∙ Поток векторного поля. Потоком векторного поля a через поверхность S называется поверхностный интеграл по поверхности S от скалярно-


166
|
r r |
|
r |
r |
r |
|
|
|
|
|
|
|
|
φ1 = ∫ rds |
= |
∫(r |
0 R + z |
0 z)r 0 dS |
= R2πRh = 2πR |
2 h , |
|
|
|
||||
|
S11 |
|
S11 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
φ3 = − ∫ |
r |
|
r |
|
r |
r |
|
= 0, |
|
|
|
|
|
r 0rz 0dS = 0 , |
т.к. r |
0 z |
0 |
||||||
|
|
|
|
|
S 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
φ2 = ∫ |
|
r |
r |
r |
|
|
|
|
|
|
|
|
|
|
(r |
0 r + z 0 h)z 0 dS = hπR |
2 |
|||||
|
|
|
|
|
|
S 3 |
|
|
|
|
|
|
|
. |
Искомый поток |
φ = φ + φ |
3 |
= 3π R 2h . |
|
|
|
|
|||||
|
|
|
|
|
1 |
|
|
|
|
|
|
|
∙Дивергенция векторного поля. Под дивергенцией вектора a в
точке М понимают предел отношения потока вектора a через замкнутую поверхность вокруг т. М к объему, находящемуся внутри этой поверхности при стремлении объема к нулю
|
|
|
r |
r |
|
|
r |
(M ) = lim |
|
∫ ads |
|
||
diva |
V →0 |
|
|
. |
(П1.16) |
|
V |
|
|||||
|
|
|
|
|
|
|
Это определение поясняет физический смысл дивергенции, как объем- |
||||||
ной плотности потока вектора |
a в точке M. |
|
|
|
|
Точка М, в которой diva > 0 называется источником. Если diva < 0
точка М является стоком. Формула (П1.16) |
дает инвариантное определение ди- |
||||||||||||||||||||
вергенции и, как из неё следует, |
|
дивергенция есть скалярная функция точек |
|||||||||||||||||||
пространства. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Выражения дивергенции в трех системах координат имеют вид |
|||||||||||||||||||||
r |
= |
∂a |
x + |
∂a y |
+ |
∂a |
z |
|
|
|
|
|
|
|
|||||||
diva |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
∂y |
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
∂x |
|
|
|
|
|
|
∂z |
|
|
|
|
|
|
|
|||||
r |
|
1 ∂(ra |
|
) |
|
|
|
1 ∂a |
|
∂a |
|
|
|
|
|
||||||
diva |
= |
|
|
∂r |
r |
|
+ |
|
|
|
|
α + |
|
|
z |
|
|
|
(П1.17) |
||
r |
|
|
r |
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
∂α |
∂z |
|
|
|
|
|||||||||
r |
|
1 |
∂(r |
2a |
|
) |
|
|
|
|
1 |
|
∂(a sin(θ )) |
|
1 |
∂a |
|||||
diva |
= |
|
|
|
|
|
r |
|
|
+ |
|
|
|
|
|
|
θ |
+ |
|
α |
|
r 2 |
∂r |
|
|
|
r sin(θ ) |
|
∂θ |
r sin(θ ) |
∂α |
||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||
Свойства дивергенции: |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
r |
r |
|
|
|
|
|
|
|
|
r |
|
r |
|
|
= Cdiva |
|
|||
div(a1 |
+ a2 ) = diva1 |
+ diva2 |
divCa |
|

167
div(ϕ a) = ϕ diva + a gradϕ
div[ab] = b × rota - a × rotb |
|
|
(П1.17) |
||||||||
Пример П1.3. Вычислить div r |
в прямоугольной и сферической сис- |
||||||||||
темах координат, где r - радиус вектор, определяемый по формуле |
|
||||||||||
|
r |
r |
r |
|
|
|
r |
0 z . |
|
(П1.18) |
|
|
r = x 0 x + y |
0 y + z |
|
||||||||
|
r |
∂x ∂y ∂z |
|
|
|
||||||
Тогда, согласно (П1.17) |
divr |
= ∂x + |
∂y |
+ ∂z |
= 3 . |
|
|
||||
|
|
r |
r |
0 r . Тогда в формуле П1.17 для сфе- |
|||||||
Представим радиус вектор как r |
= r |
||||||||||
рических координат проекция ar |
= r и |
r |
= |
|
1 ¶(r 2 × r) |
= 3 |
|
||||
divr |
|
|
|
¶ r |
|
||||||
|
r 2 |
|
|||||||||
|
|
|
|
|
|
|
|
|
|||
Из данных примеров следует, что вычисление дивергенции |
не зависит |
||||||||||
от выбора системы координат. |
|
|
|
|
|
|
|
|
|
|
|
Пример П1.4. Вычислить div от произведения 2x y z r , где |
r - радиус |
||||||||||
вектор. |
|
|
|
|
|
|
|
|
|
|
|
Используя третье выражение (П1.17), получим для заданного условия div(2xyzr) = 2xyz × divr + r × grad(2xyz) . И далее, решая каждое слагаемое, πпо-
лучим, что |
div(2xyzr) =12xyz . |
|
|
|
· |
Теорема Остроградского – Гаусса. Поток вектора |
a через замк- |
||
нутую поверхность S равен интегралу от дивергенции этого вектора, взятого по |
||||
объему V, ограниченному этой поверхностью |
|
|||
|
r |
= |
r |
|
|
∫ a d S |
∫ div a dv . |
(П1.19) |
|
|
s |
|
v |
|
Пример П1.5 Решить задачу примера П1.3 с помощью теоремы Остро- |
||||
градского – |
Гаусса. |
|
|
|
Решение. Как следует из предыдущих примеров дивергенция радиуса – вектора равна 3. Поэтому интеграл в правой части (П1.19) будет равен 3V =3·πR2h, где
V – объем цилиндра. Как видно, вычисление потока с помощью теоремы Остроградского – Гаусса значительно упростило решение задачи. Этот вывод спра-

168
ведлив для большинства встречающихся задач вычисления потока, но, конечно,
только для замкнутых поверхностей. |
|
|
||
· |
Циркуляция вектора. |
Циркуляцией вектора a по замкнутому- |
||
контуру L называется криволинейный интеграл вида |
|
|||
|
|
r |
|
|
|
|
Ц = ∫ adl |
|
(П1.23) |
|
|
L |
|
|
Элемент dl направлен по касательной в каждой точке контура L Направ- |
||||
ление элемента dl указывает направление обхода контура |
|
|||
при интегрировании. |
|
|
|
|
Пример П1.6. Вычислить циркуляцию векторного |
|
|||
r |
r |
|
- окружности |
|
поля a = |
R ×α o + 5sinα ×θ o по контуру L |
|
||
с радиусом R=1 (рис. П1.4). Направление обхода контура - в |
|
|||
положительном направлении сферической координаты ά. |
Рис. П1.4 |
|||
|
||||
|
r |
r |
r |
2dα |
Решение. Определим a × dl = |
(Rα 0 + 5sinα ×θ 0 ) × Rα 0 dα = R |
При интегрировании по окружности угол ά изменяется от 0 до 2π. В результате получим Ц=2πR2.
· Ротор векторного поля. Теорема Стокса.
Ротором вектора a в точке М называется вектор равный пределу отношения следующего вида
|
|
r |
× dS ] |
|
r |
|
∫[a |
||
|
S |
|
|
|
rota(M ) = lim |
V →0 |
|
(П1.21) |
|
|
|
|||
|
|
DV |
На практике чаще определяют ротор как предел отношения циркуляции по контуру в окрестности точки М к площади контура, при стремлении площади к нулю. При этом определяется не сам ротор как вектор, а его проекция на нормаль к контуру.
|
|
|
|
r |
× dl |
|
r |
|
|
|
∫ a |
||
|
n = lim |
|
L |
|
(П1.22) |
|
|
|
|
||||
rota(M ) |
|
S →0 |
|
|||
|
|
|
||||
|
|
|
|
DS |

169
Нормаль к контуру связана с направлением обхода контура правилом правого винта, когда нормаль направлена в сторону движения винта, если сам винт вращается в направлении обхода контура.
Выражения для ротора в различных системах координат обычно представляют в виде определителей
|
|
|
|
|
|
r |
|
|
r |
0 |
|
|
|
|
r |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
x 0 |
|
|
y |
|
|
|
|
z 0 |
|
|
|
|
|
|
|
|
|
||||||||
|
r |
= |
|
¶ |
¶ |
|
¶ |
|
|
|
|
|
|
в прямоугольной системе координат, |
(П1.23) |
|||||||||||||||
rot a |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
¶x |
|
|
¶y |
|
|
|
¶z |
|
|
||||||||||||||||||||
|
|
|
|
|
ax |
|
|
a y |
|
|
|
az |
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
r |
0 |
|
|
|
|
|
|
|
|
|
|
r |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
r |
|
|
|
|
r |
|
|
|
|
z 0 |
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
α 0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
r |
|
|
|
|
r |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
¶ |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
r |
= |
¶ |
|
|
|
|
|
|
¶ |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
rot a |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
- в цилиндрической системе координат, |
(П1.24) |
|||||||||||||
|
¶r |
¶α |
|
|
|
|
¶z |
|
|
|||||||||||||||||||||
|
|
|
ar |
r aα |
|
|
|
|
az |
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
r |
|
|
|
|
|
|
|
r |
|
|
r |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
θ 0 |
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
r 0 |
|
|
|
|
|
|
α 0 |
|
|
|
|
|||||||||||||
|
|
|
r |
2 sinθ |
|
|
|
r sinθ |
|
|
r |
|
|
|
|
|||||||||||||||
r |
= |
|
|
¶ |
|
|
|
|
|
¶ |
|
|
|
|
|
|
¶ |
|
|
|
|
|||||||||
rot a |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
в сферической системе координат (П1.25) |
||
|
|
|
|
¶r |
|
|
|
|
|
|
¶θ |
|
|
¶α |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
ar |
|
|
|
|
|
r aθ |
|
r sinθ |
× aα |
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Раскрываются эти определители по первой строке, как показано в формуле П1.6.
Теорема Стокса утверждает, что циркуляция вектора по замкнутому контуру L равна потоку ротора от этого вектора через поверхность S опирающуюся на этот контур.
r |
× dl |
r |
× dS |
(П1.26) |
∫ a |
= ∫ rota |
|||
L |
|
S |
|
|
Как следует из теоремы Стокса модуль вектора rota является плотностью циркуляции.
Свойства ротора: |
|
|
|
r |
r |
± rotb .2. |
rotCa = Crota |
1. rot(a |
± b ) = rota |
(П1.27)

171
С формальной точки зрения это – вектор, компонентами которого являются символы производных. Применяя этот оператор к скалярной функции φ(х.у.z), по правилу умножения скаляра на вектор получим выражение
gradφ(x,y,z). |
r |
0 |
∂ϕ |
r |
0 |
∂ϕ |
r |
0 |
∂ϕ |
= gradϕ |
Ñϕ = x |
¶x |
+ y |
¶y |
+ z |
¶z |
|||||
|
|
|
|
|
|
|
|
Легко также показать, что скалярное произведение наблы на вектор дает дивергенцию вектора, а их векторное произведение – ротор вектора.
|
|
r |
r |
|
r |
|
r |
|
(П1.29) |
|
|
|
|
|
|||||
|
|
(Ña) = diva, |
|
Ña |
|
= rota |
|
||
|
Оператор Ñ, обладая формально свойствами вектора и одновременно |
||||||||
свойствами дифференциального оператора, |
|
позволяет вычислять некоторые- |
|||||||
сложные выражения векторного анализа. Например, |
|
|
|||||||
r |
r |
r r |
|
|
|
|
r |
r |
(П130). |
rotrota |
= [Ñ[Ña |
]]= Ñ × (Ña) - a(ÑÑ) = grad (diva) - Ñ2 a |
Здесь использовано правило раскрытия двойного векторного произведения (П1.7). Появившийся при этом преобразовании оператор Ñ2 называется оператором Лапласа и играет важную роль в электродинамике и вообще физике. Применительно к скалярной функции координат он выражается следующим образом
|
2 |
¶2ϕ |
¶2ϕ |
¶2ϕ |
|
|
Ñ ϕ (x, y, z) = |
¶x2 |
+ ¶y 2 |
+ ¶z 2 |
(П1.31) |
Применительно к вектору, заданному в прямоугольной системе коорди- |
|||||
|
нат он применяется к каждой проекции вектора, как к скаляру |
||||
r |
r |
|
r |
r |
(П1.32) |
Ñ2a(x, y, z) = x 0Ñ2ax |
(x, y, z) + y |
0Ñ2a y (x, y, z) + z 0Ñ2 az (x, y, z) |
|||
При использовании оператора набла необходимо иметь в виду: |
|
1) Оператор аналогичный оператору Ñ не может быть введен в цилиндрической или сферической системах координат. в этом нетрудно убедиться, ес-
r |
r |
ли ввести его так, чтобы Ñϕ = gradϕ . Но тогда (Ña) ¹ diva и [Ña |
] ¹ rota . |
2) Как отмечалось выше такие характеристики скалярных и векторных полей, как градиент, дивергенция и ротор являются инвариантными к выбору системы координат и поэтому соотношения, полученные с помощью оператора
172
Ñ в прямоугольной системе остаются справедливыми и в другой, но записан-
ные через ротор, дивергенцию и градиент, т. е. без оператора Ñ.
3)Оператор Лапласа (иногда он обозначается символом ) в криволинейных координатах применительно к скалярным функциям определяется как
Ñ2ϕ = div(gradϕ ) , (П1.33)
а к векторным - из соотношение (П1.30)
|
|
|
|
|
|
|
|
|
|
|
r |
|
|
|
|
|
|
|
|
|
r |
|
|
|
r |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(П1.34) |
|
|
|
|
|
|
|
|
|
Ñ2a |
= graddiva - rotrota |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
Оператор Лапласа в цилиндрической системе координат представляется |
|||||||||||||||||||||||||||||||||||||||||
в виде |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Dϕ = Ñ |
2 |
×ϕ |
|
|
1 |
|
|
¶ |
|
¶ϕ |
|
1 |
|
|
¶2ϕ |
|
|
¶2ϕ |
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
= |
|
|
× |
|
|
|
r × |
|
|
+ |
|
|
|
× |
|
|
|
+ |
|
|
|
|
|
|
(П1.35) |
|||||||||
|
|
|
|
|
|
|
|
r |
|
|
|
r |
2 |
|
¶α |
2 |
|
¶z |
2 |
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
¶r |
|
¶r |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
в сферической - |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Dϕ = |
1 |
|
¶ |
|
2 |
|
¶ϕ |
|
|
|
|
|
1 |
|
|
|
|
|
¶ |
|
|
|
¶ϕ |
|
|
|
|
|
|
1 |
|
|
|
¶2ϕ |
|||||||
|
|
× |
|
r |
|
× |
|
+ |
|
|
|
|
|
|
|
|
|
× |
|
|
sinθ × |
|
|
|
+ |
|
|
|
|
|
|
|
× |
|
|
||||||
|
2 |
|
|
|
2 |
|
|
|
|
|
θ |
|
|
|
|
|
|
2 |
|
|
2 |
θ |
α 2 |
||||||||||||||||||
|
r |
|
¶r |
|
|
¶r |
|
|
r |
×sin |
|
¶θ |
|
|
|
¶θ |
|
|
|
r |
×sin |
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
¶ |
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(П1.36) |
· Некоторые тождества векторного анализа. Приведем основные тождества векторного анализа. В их справедливости можно убедиться либо непосредственным вычислением, либо с помощью оператора набла.
|
|
) º 0, |
rot(gradϕ) º 0, |
div(gradϕ) º Ñ2ϕ |
|
||||||
div(rota |
|
||||||||||
r |
r |
|
r |
r r |
|
|
|
r |
|
||
rot(ϕa) º [gradϕ × a |
]+ ϕ × rota, div(ϕa) º a |
× gradϕ + ϕdiva |
|
||||||||
|
|
r |
|
r |
r |
r r |
r |
r |
r |
r |
(П1.37) |
rotrota = grad (diva) - Ñ2a, |
div[ab] º b |
× rota |
- a |
× rotb |
|||||||
Векторные поля, |
в которых rota = 0 называются потенциальными. |
Если выполняется условие diva = 0 во всех точек внутри некоторой области, то поле в ней называется соленоидальным. Скалярные поля, удовлетворяющие уравнению Лапласа, называются гармоническими..