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

Московский_государственный_технический_университет

.pdf
Скачиваний:
121
Добавлен:
23.03.2016
Размер:
1.25 Mб
Скачать

Пусть в области G определено векторное поле

 

P; Q; R , L – замкнутый

a

контур, расположенный в области G, σ – поверхность, ограниченная гладким

или кусочно-гладким контуром L, n0 M cos ,cos ,cos - единичный вектор

нормали на выбранной стороне поверхности σ. Пусть функции P; Q; R

непрерывны вместе со своими частными производными. Тогда справедлива

формула Стокса:

 

 

 

R

 

Q

P

 

R

 

Q

 

P

 

Pdx Qdy Rdz

y

 

z

cos

z

 

cos

x

 

y

cos d . (15.9)

 

 

 

 

 

 

 

x

 

 

 

 

L

 

 

 

 

 

 

 

 

 

 

Ориентация контура L согласована с ориентацией поверхности σ по правилу

правого винта. (См.рис. 15.4) Или:

 

 

R

 

Q

P

 

R

Pdx Qdy Rdz

y

 

z

dydz

z

 

x

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

 

Q

 

P

(15.10)

dxdz

 

dxdy .

 

 

x

 

y

 

Левая часть формулы Стокса – это циркуляция векторного поля вдоль a M

контура L, а правая представляет собой поток через поверхность σ

векторного поля . Формулу Стокса можно записать в векторной форме:

rot a

a, dr

L

 

 

 

 

 

 

 

rot a, n0

d rot a, d .

 

 

 

 

 

 

 

 

 

 

 

Рисунок 15.4

(15.11)

Физический смысл формулы Стокса:

циркуляция векторного поля a M

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

натянутую на этот контур.

Замечание. Формула Стокса остаѐтся справедливой и в случае,

когда поверхность σ является плоской областью, параллельной

51

какой-нибудь координатной плоскости. Тогда формула Стокса превращается

в формулу Грина:

 

 

Q

 

P

Pdx Qdy

x

 

dxdy .

L

 

 

y

Пример 15.3. Проверим, является ли потенциальным векторное поле

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a 3x2 y3 z3i 3z3 x3 y2

2y j

3x3 y3z2 2z k .

 

 

 

 

 

Решение. По формуле (15.8) имеем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

j

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rot a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

y

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3x2 y3 z3

3x3 y2 z3 2 y

3x3 y3 z2 2z

 

 

 

 

 

 

 

 

 

 

 

 

 

9x2 y2 z3 9x2 y2 z3

 

 

 

 

 

 

9x3 y2 z2 9x3 y2 z2 i

j 9x2 y2 z3

9x2 y2 z

3

k

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ответ: данное поле является потенциальным, т.к. rot a 0. .

 

Пример 15.4. Вычислим циркуляцию векторного поля

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a 3x

y z i

2x 4y z j

5z 3y k

по контуру треугольника ABC , где A(1,1,0),

B(0,0,2), C(3,0,1), двумя способами: 1) с

помощью криволинейного интеграла, 2) по формуле Стокса.

Выясним, как зависит циркуляция от расположения контура в данном поле.

Рисунок 15.5

Решение.

1) Вычислим циркуляцию по формуле (15.3):

 

 

 

 

 

 

 

 

 

 

Ц

 

a, dr

 

3x y z dx 2x 4y z dy 5z 3y dz

 

 

 

(см.

 

ABCA

 

 

 

L

AB

 

BC

CA

рис.15.5)

2) Составим уравнения прямых, содержащих стороны треугольника.

52

 

 

 

x 1

 

 

 

 

 

y 1

 

 

 

 

z

y x

 

 

 

dy dx

 

AB 1,1, 2 AB :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

;

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

1

 

 

 

 

2

z 2x 2

 

dz

2dx

 

 

 

 

 

x

 

 

y

 

 

 

 

z 2

 

 

 

y 0

 

dy 0

 

 

BC 3,0, 1 AB :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

3 0

 

 

 

 

 

 

 

 

 

 

x 3z 6

 

dz 3dz

 

 

 

 

 

x 1

 

 

 

 

y 1

 

 

 

z

x 2t 1

 

dx 2dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

AC 2,

1,1 AC :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t y t 1;

 

dy dt .

 

 

2

 

 

 

 

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z t

 

 

 

dz dt

 

Вычислим криволинейные интегралы, преобразовав их в определѐнные.

 

0

 

 

 

 

 

 

 

 

 

22x 20 dx 9

 

 

 

 

AB

1

 

 

 

 

 

 

 

 

 

1

 

 

29

 

 

 

 

 

29z 54 dz

 

 

 

Ц 9.

 

 

4

Просуммировав, получим ответ:

2

BC

2

 

 

 

 

 

AB BC CA

 

0

19

 

 

 

 

 

19t 1 dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 1

 

 

 

 

 

 

 

CA

1

 

 

 

 

 

 

 

 

2) Проверим полученный ответ, вычислив циркуляцию по формуле Стокса:

 

 

 

 

ds . (см.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ц rot a, n0

 

(15.11)). В

 

качестве

поверхности

σ, на которой

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

располагается заданный контур,

выберем плоскость

треугольника ABC.

Тогда

 

 

 

 

 

 

 

 

 

 

Единичный вектор

нормали,

составляющий

n

BA BC 1, 5, 3 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

1, 5, 3 .

 

 

 

 

 

 

 

 

 

 

 

острые углы с осями координат: n0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

35

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

по формуле (15.8) и скалярное произведение:

 

 

4,1,3 ,

Найдѐм rot a

rot a

 

 

 

1

 

 

18

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rota, n0

 

 

 

 

 

4 5 9

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

35

35

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

18

 

 

 

18

 

 

 

 

 

18

 

 

1

 

 

 

 

9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Получим

Ц rot a, n0

ds

 

 

 

ds

 

 

 

S ABC

 

 

 

 

 

 

BA BC

 

 

 

 

35 9 ,

 

 

 

 

 

 

 

 

 

2

 

 

 

 

35

 

35

 

35

35

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

т.е. тот же ответ.

3) Исследуем поведение циркуляции при перемещении нашего контура в пространстве. Для этого используем определение скалярного произведения.

Ц rot a, n0 ds

rota

n0 cos ds

27 cos ds , где – угол между

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

векторами

 

 

 

 

 

 

 

 

rot a

и n0 .

 

 

 

 

 

53

, n0

Видим, что циркуляция зависит от этого угла . Отсюда следует, что если треугольник ABC перемещается в пространстве параллельно своему исходному положению, то угол не меняется, и циркуляция по контуру остаѐтся равной 9, а при повороте контура меняется направление вектора что влечѐт изменение циркуляции.

Циркуляция достигнет максимального значения, когда cos 1, т.е. когда

 

 

 

Цmax

 

 

 

 

 

3

 

 

 

 

27ds

27 S ABC

 

 

 

 

rot a

n0

. В этом случае

 

105 .

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ответ: Ц=9, Цmax 32 105 .

Занятие 16.

Оператор Гамильтона. Повторные дифференциальные операции.

Ауд.: ОЛ-6 №№ 2382, 2385(б), 2398(б) или ОЛ-5 №№ 10.122, 126, 126, 127, 149 (а, в).

1.Оператор Гамильтона.

Определение. Оператором Гамильтона называют символический вектор

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(набла):

 

i

 

 

j

 

 

k

или

 

 

;

 

;

 

.

x

y

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

y

 

z

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

2.Дифференциальные операции первого порядка.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1)

Умножая вектор

на скалярную функцию U x, y, z , получим grad U :

 

 

 

 

 

 

 

 

 

 

 

 

 

U

U

U

 

 

 

U

 

i

 

 

j

 

 

 

 

k U

 

 

i

j

k

grad U .

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

y

 

 

 

z

 

 

 

 

x

 

y

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P; Q; R даѐт:

2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Скалярное произведение вектора на вектор-функцию a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

div a

 

, a

 

 

P

 

 

Q

 

 

R div a .

 

 

 

 

 

 

 

y

 

z

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

54

 

 

 

 

 

 

 

 

 

 

 

 

 

P; Q; R даѐт

3) Векторное произведение вектора

на вектор-функцию a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

j

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rot a

:

a

 

 

 

 

 

 

 

rot a.

 

 

x

 

y

 

z

 

 

 

 

 

 

 

P

 

Q

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Применять оператор нужно, соблюдая следующие правила:

1) Линейность оператора . Пусть 1, 2 , 3 , . - вещественные или комплексные числа, F,G, H , - скалярные или векторные функции.

Тогда 1F 2G 3H 1 F 2 G 3 H

2) Действие оператора на произведение функций.

 

 

 

 

 

F G H F GH

F

 

 

 

 

 

 

 

 

 

G H

FG H .

 

 

 

 

(Символ «↓» указывает на объект действия оператора).

Докажем, например, что a

div

 

 

 

 

 

 

 

 

 

 

div a

b

a

b

a

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

b

rota

a

rotb .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a b

a ,b b a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

, b

b

, a

b

rota

a

rotb .

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

 

 

 

 

 

 

 

 

 

 

 

 

векторов p a

, b

p,

a

b .

 

 

 

 

 

 

 

2.Дифференциальные операции второго порядка.

 

 

Пусть в области G заданы скалярное поле U x, y, z и

векторное поле

 

 

 

 

 

 

 

имеют в области G непрерывные

a P; Q; R , такие, что функции U , P, Q, R

 

 

 

 

 

 

 

 

 

 

 

 

частные

производные

 

второго порядка.

Тогда

и

gradU -

 

rot a

дифференцируемые векторные поля,

 

-

дифференцируемое

скалярное

div a

 

 

 

 

 

 

 

 

 

 

 

 

поле.

В этих

случаях к полям

и

gradU применимы

операции

rot a

55

вычисления дивергенции и ротора,

 

 

 

 

 

 

 

 

 

 

а к полю div a - операция вычисления

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

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

div gradU , rot grad U

, div rot a

rot rot a

, grad div a .

 

 

Определение.

 

Операция

 

 

 

 

 

 

 

 

 

 

называется

оператором

Лапласа и

 

 

 

div grad U

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U . Оператор Лапласа

обозначается

 

символом

(лапласиан): div grad U

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

можно записать с помощью оператора набла :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

div grad U

, U

2U U .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

2

2

 

 

При этом

2

 

 

 

i

 

 

j

 

 

 

k

 

 

 

 

 

 

 

 

 

.

 

 

 

y

 

 

 

x

2

y

2

z

2

 

 

 

 

 

 

 

x

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2U

 

 

2U

 

 

2U

 

 

 

 

 

 

 

 

Тогда

 

U 2U x2

 

y2

 

z2 .

 

 

 

 

 

 

 

Определение .

 

Уравнение

 

 

в

 

 

частных

 

производных

U 0 или

2U

2U

2U

0 называют уравнением Лапласа.

 

 

x2

y2

z2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Дважды дифференцируемую функцию

U x, y, z ,

удовлетворяющую этому

уравнению, называют

гармонической в области

G. Для плоской области

уравнение Лапласа примет вид: 2U 2U 0 .

x2 y2

Примеры гармонических функций:U ax by c , U 3x 2xy 7zx ,

U r , r x2 y2 z2 , r 0 .

4.Запись дифференциальных операций второго порядка в операторной

 

 

 

 

 

форме.

 

 

 

 

 

 

rot grad U U 0 - потенциальное поле grad U является безвихревым

полем.

 

 

 

 

 

 

 

 

 

 

 

div rot a ,

 

 

, a

0 - векторное поле роторов соленоидально.

 

 

 

 

 

 

rot rot a

grad div a

a

56

 

 

 

 

 

 

grad diva

rot rota

a ,где

 

 

 

 

 

 

 

 

 

a P

x, y, z i Q x, y, z j

R x, y, z k P Q R

- вектор-функция,

полученная в результате применения оператора Лапласа функциям P,Q, R .

Определение. Векторное поле

 

 

q M называется гармоническим в области G,

если в этой области оно одновременно является потенциальным и соленоидальным, т.е. для любой точки гармонического поля выполняются

 

0

 

0 . Для односвязных областей из первого условия

условия rot q

и div q

следует, что

 

 

 

 

 

q

grad ,

а из второго, что div q

div grad 0 , т. е.

функция φ, потенциал поля – решение уравнения Лапласа и гармоническая функция.

Пример 17.1. Проверим, является ли гармонической функция

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U ln

 

, r x2

y2 , r 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение. Вычислим частные производные второго порядка функции

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2U

 

 

2U

 

 

 

 

 

 

U

 

ln x2 y2 и вычислим U x2

 

y2 .

 

 

 

 

 

 

2

 

 

 

 

 

 

U

 

 

 

 

 

 

x

U

 

 

y

 

 

2U

 

x2 y2

2U

 

y2 x2

 

 

 

 

 

,

 

 

 

 

,

x2

 

 

,

y2

 

 

,

 

x

x2 y2

y

x2 y2

x2 y2 2

x2 y2 2

 

U 2U2

2U2

0 функция

U ln

1

- гармоническая.

Тогда поле

r

 

 

 

x

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

y

 

 

 

 

 

 

2U 2U

 

 

 

q

grad U

 

 

i

 

 

j -гармоническое и div q

x2

y2

0, rot q

0

x2

y2

x2 y2

57