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

Основы высшей математики для инженеров 2009

.pdf
Скачиваний:
226
Добавлен:
12.03.2016
Размер:
10.04 Mб
Скачать

6.2. Векторное поле

291

 

 

 

DD(F,

n

)dS

 

 

 

 

S

 

div F

lim

.

V

 

 

V 0

 

Формулу Остроградского – Гаусса можно записать в виде

DD(F, n)dS DDDdiv F dV.

SV

Пусть векторное поле — поле скоростей v движущейся жидкости с плотностью + 1. W — область, ограниченная замкнутой поверхно стью S. Тогда

DD(v, n)dS DDDdiv v dW.

SW

DD(v, n)dS — количество жидкости, протекающей через S в едини

S

цу времени, точнее, алгебраическая сумма количеств втекающей и вы текающей через S жидкостиM62(поток берется по внешней стороне по верхности, поэтому если вектор скорости v направлен наружу, то (v, n) 0, если — внутрь, то (v, n) 0).

Если

DD(v, n)dS 0,

S

то из объема W вытекает жидкости больше, чем втекает. Если же

DD(v, n)dS 0,

S

то больше жидкости втекает.

Предположим, что в некоторой точке M div v 0. В силу непре рывности частных производных div v 0 и в достаточно малой окрестности W, ограниченной сферой S с центром в точке М. Но тогда

DDDdiv v dW 0

W

и, следовательно,

DD(v, n)dS 0.

S

Поэтому из W через поверхность S жидкости вытекает больше, чем втекает. В этом случае точку М называют источником. Если div v 0,

292

Глава 6. Скалярное поле

то точка М — сток. Если div v 0, то поток равен 0, следовательно, ко личества втекающей и вытекающей жидкости равны.

П р и м е р. Для векторного поля F r xi yj zk вычислить его поток через произвольную замкнутую поверхность S.

Решение.

 

 

 

 

 

 

 

HP

 

HQ

 

HR

 

 

 

 

 

 

 

 

DD(F,

n

)dS DDDdiv F dV

DDD

 

 

 

 

 

dV 3V.

 

Hy

 

S

 

 

 

V

Hx

 

 

Hz

 

 

 

V

 

 

 

 

 

6.3. Формула Стокса. Ротор и циркуляция векторного поля. Оператор Гамильтона

Формула Стокса. Условия независимости криволинейного интеграла от пути интегрирования. Ротор и циркуляция векторного поля. Вектор ная запись формулы Стокса. Оператор Гамильтона и его применение.

6.3.1. Формула СтоксаM62. Условия независимости криволинейного

интеграла от пути интегрирования

Формула Стокса является обобщением формулы Остроградского – Грина и позволяет свести вычисление криволинейного интеграла по замкнутому контуру L к вычислению потока через поверхность, огра ниченную этим контуром.

Теорема. Пусть функции P(x, y, z), Q(x, y, z), R(x, y, z) непрерывны вме2 сте со своими частными производными первого порядка на поверхно2 сти S, ограниченной замкнутым контуром L.

Тогда

DPdx Qdy Rdz

L

 

7HR HQ

HP HR

HQ HP

:

(1)

DD9

 

 

 

cos)

 

 

 

cos*

 

 

 

cos L< dS.

 

 

 

 

 

 

S

8 Hy

 

Hz

Hz

 

Hx

Hx

 

Hy

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если при этом сторона поверхности выбрана (т. е. выбрано направление единичного вектора нормали), то направление обхода контура L выбира2 ется так, чтобы поверхность оставалась слева.

Доказательство этой теоремы мы опускаем.

6.3. Формула Стокса. Ротор и циркуляция векторного поля…

293

З а м е ч а н и е. Если поверхность S — часть плоскости, параллельной плоскости XOY, то dz 0, и мы получим формулу Грина:

 

 

HQ

 

HP

DPdx Qdy DD

 

 

 

dxdy.

 

 

L

S

Hx

 

Hy

Из формулы (1) следует, что если

HR

 

HQ

,

HP

 

HR

,

HQ

 

HP

,

(2)

 

 

 

 

 

 

Hy Hz

Hz Hx

Hx Hy

 

то криволинейный интеграл по любой пространственной замкнутой кривой равен 0, т. е.

DPdx Qdy Rdz 0.

(3)

L

 

Таким образом, условия (2) являются достаточными для того, что бы криволинейный интеграл (3) не зависел от пути интегрирования.

Как и в случае плоской кривой, можно показать, что условия (2) являются не только достаточными,M62но и необходимыми условиями не зависимости криволинейного интеграла (3) от пути интегрирования. При выполнении этих условий подынтегральное выражение есть пол ный дифференциал некоторой функции U(x, y, z):

dU(x, y, z) Pdx Qdy Rdz

и, следовательно, для любых двух точек М и N имеем

(N )

DdU U(N) U(M).

(M )

6.3.2. Ротор и циркуляция векторного поля. Векторная запись формулы Стокса

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

F P(x, y, z)i Q(x, y, z)j R(x, y, z) k.

Вектор

HR

 

HQ

HP HR

HQ

 

HP

 

(4)

 

 

 

 

 

 

 

i

 

 

 

j

 

 

 

k

 

 

 

 

 

 

Hy

 

Hz

Hz

 

Hx

Hx

Hy

 

294

Глава 6. Скалярное поле

называется ротором (или вихрем) векторного поля и обозначается rot F . Для того чтобы запомнить выражение (4), введем определитель

i

j

k

H H H .

Hx Hy Hz

P Q R

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

j

 

k

 

 

 

 

H

 

H

 

H

 

.

rot F

Hx

 

Hy

Hz

 

 

 

 

 

 

 

 

 

P

 

Q

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

Можно показать, что rotM62F не зависит от системы координат, а оп ределяется исходным векторным полем. Полю вектора F соответству ет новое векторное поле — поле его вихрей.

П р и м е р. Пусть дано поле скоростей абсолютно твердого тела, вращающегося с постоянной угловой скоростью I. Если тело вращает ся вокруг оси OZ, то v Iyi Ixj.

Найдем rot v. В данном случае P Iy, Q Ix, R 0, поэтому

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

j

k

 

 

 

 

 

rot

 

 

 

H

 

H

 

H

 

 

 

 

 

 

 

 

 

 

 

k (I I) 2kI.

v

Hx

Hy

Hz

 

 

 

 

 

 

 

 

 

 

 

 

Iy

Ix

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Итак, rot v 2Ik . Модуль вектора 2I, следовательно, равен удво енной угловой скорости, с которой тело вращается вокруг оси OZ. От сюда происхождение названия «ротор», т. е. «вращение».

Возвратимся к формуле Стокса. Заметим, что

 

7 HR HQ

HP HR

HQ HP

:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

DD9

 

 

 

cos)

 

 

 

cos*

 

 

 

cos L< dS DD(rot F, r)dS.

 

 

 

 

 

 

S

8 Hy

 

Hz

Hz

 

Hx

Hx

 

Hy

;

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6.3. Формула Стокса. Ротор и циркуляция векторного поля…

295

 

 

Интеграл

DPdx Qdy Rdz называется

циркуляцией

векторного

 

 

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

поля F по

замкнутому

 

контуру L и

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

D(F, dr

), где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

dx i dy j dz k .

 

 

 

 

 

 

 

 

 

dr

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5)

 

 

 

 

 

 

 

D(F

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, dr

) DD(rot F, n)dS.

 

 

 

 

 

 

 

L

 

 

 

S

 

 

 

 

 

 

Соотношение (5) означает, что циркуляция вектора F по замкну тому контуру L равна потоку ротора вектора F через поверхность S, ограниченную контуром L.

П р и м е р. Вычислить циркуляцию вектора F x2 yi j z k вдоль контура x2 y2 R2 , z 0, приняв в качестве поверхности полусферу z R2 x2 y2 .

Решение.

M62i j k

rot F H H H x2 k . Hx Hy Hz

x2 y 1

z

D(F, dr) DD(rot F, n)dS DD x2 cos L dS DD x2 dxdy

L

S

S

 

 

 

 

Nxy

 

 

2'

R

R

4 2'

1 cos2

(

 

'R

4

 

Dd(Dr2 cos2 ( rdr

 

D

d(

 

.

4

2

 

 

 

0

0

0

4

 

 

 

 

 

 

 

 

 

 

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

Мы рассмотрели основные дифференциальные величины вектор ного анализа: градиент, дивергенцию и ротор. Эти величины могут быть записаны короче, если ввести символ, называемый оператором Гамильтона, или оператором S («набла»).

S H i H j H k.

Hx Hy Hz

296

Глава 6. Скалярное поле

Будем рассматривать этот оператор условно как вектор. Установим правила действия с ним.

1. Пусть U(x, y, z) — скалярная функция. Под произведением век тора S на U условимся понимать следующий вектор:

H

 

H

 

H

 

 

 

HU

 

HU

 

HU

 

 

 

 

 

 

 

k.

SU

 

i

 

j

 

k U

 

i

 

j

 

 

Hy

 

Hx

Hy

Hz

Hx

 

 

Hz

 

 

 

 

 

Но HU i HU j HU k gradU , следовательно, SU gradU .

Hx Hy Hz

2. Под скалярным произведением вектора S на вектор функцию F Pi Qj Rk условимся понимать следующую величину:

 

 

 

 

 

 

 

 

H

 

 

 

H

 

 

 

 

H

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

HP

 

 

 

HQ

 

HR

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(S,

F)

 

 

 

i

 

 

 

j

 

 

 

k,

Pi

Qj

Rk

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Hx Hy Hz

 

 

 

 

 

 

 

 

 

Hx Hy Hz

 

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

 

 

 

 

M62Qj Rk :

 

 

 

 

 

 

 

 

 

 

 

 

F Pi

 

 

 

 

 

 

 

 

 

Но

 

HP

 

HQ

 

 

HR

 

div F , следовательно,

(S, F) div F .

 

 

 

 

 

Hx

 

Hy

 

 

 

Hz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3. Рассмотрим,

наконец,

 

векторное

 

 

произведение

вектора S на

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

j

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H

 

H

 

H

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S ! F

 

 

rot F

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Hx

Hy

 

Hz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P

Q

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Итак, операции с вектором S производятся по правилам вектор

ной алгебры. При этом следует учесть, что произведения H , H , H на

Hx Hy Hz

скалярные функции заменяются частными производными этих функ ций по x, y, z соответственно.

6.4. Потенциальное и соленоидальное поля. Уравнение Лапласа

297

6.4. Потенциальное и соленоидальное поля. Уравнение Лапласа

Потенциальное поле. Отыскание потенциала. Соленоидальное поле. Уравнение Лапласа.

6.4.1. Потенциальное поле. Отыскание потенциала

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

F P(x, y, z)i Q(x, y, z) j R(x, y, z) k называется

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.

 

 

 

 

безвихревым, если во всех его точках rotF

 

 

 

 

 

 

 

 

 

0 получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из условия rotF

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

Q

,

 

 

P

 

R

,

 

 

 

Q

 

P

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

z

 

 

z

x

 

 

x

y

Но в этом случае, как было показано ранее,

 

 

 

 

 

 

 

P dx

 

 

 

P dx Q dy Rdz 0

 

 

 

 

 

 

 

 

Q dyM62( ) ( ).

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и поэтому

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(N )

 

 

 

 

 

 

 

 

(N )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Rdz

dU U N U M

 

(M )

 

 

 

 

 

 

 

 

(M )

 

 

 

 

 

 

 

 

 

 

С другой стороны

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dU

U

dx

U

dy

U

dz

U

P,

U

Q,

U

R.

 

 

 

 

 

 

 

x

y

 

z

 

x

 

 

 

y

 

 

z

Учитывая, что

gradU U i U j U k ,

x y z

получим

gradU Pi Qj Rk F .

Таким образом, безвихревое поле является градиентом некоторой функции U. Функция U, градиент которой равен F , называется потен# циальной функцией (потенциалом) поля F . Поле, имеющее потенци# ал, называется потенциальным.

298 Глава 6. Скалярное поле

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

Верно и обратное: всякое потенциальное поле является безвихре#

вым. Действительно, пусть

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U

 

 

U

 

 

 

U

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

 

 

i

j

 

k .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

y

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

j

 

 

 

k

 

 

 

 

 

 

 

 

 

rotF rot grad U

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

y

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U

 

 

U

 

 

U

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

y

 

 

 

z

 

 

 

 

 

2U

 

 

2U

 

2U

 

 

2U

 

 

 

 

 

 

2U

 

 

2U

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

j

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z y y z

 

 

M62

y x x y

 

 

z x

 

 

x

z

 

 

 

 

 

 

Пусть

F

Pi Qj Rk

 

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

поле. Работа

сил поля по перемещению единичной массы есть

 

A P dx Q dy Rdz.

L

Поскольку поле потенциально, т. е. P dx Q dy Rdz dU , то этот инте# грал не зависит от пути интегрирования, при этом потенциал (при вы# боре пути интегрирования специальным образом), может быть найден

по формуле

(x, y, z)

U (x, y, z)

 

P dx Q dy Rdz

 

(x0, y0, z0)

x

y

z

P(x, y0, z0)dx Q(x, y, z0)dy 2(x, y, z)dz.

x0

y0

z0

Тогда работа по перемещению единичной массы из точки M1(x1, y1, z1) в точку M2 (x2 , y2 , z2 ) есть

(x2, y2, z2)

A

P dx Q dy Rdz U (x2 , y2 , z2 ) U (x1, y1, z1).

(x1, y1, z1)

6.4. Потенциальное и соленоидальное поля. Уравнение Лапласа

299

Для плоского поля

(x2, y2)

A P dx Q dy U (x2 , y2 ) U (x1, y1).

(x1, y1)

П р и м е р. Рассмотрим материальную точку массы m, на которую действует сила тяжести, равная mq. Пусть точка перемещается в верти# кальной плоскости: F mg j.

 

 

P 0,

Q mg,

R 0.

i

j

k

 

 

rotF

 

i 0 j 0 k 0 0.

x

y

z

 

 

0

mg

0

 

Рис. 6.2

 

 

 

 

Таким образом, поле сил тяжестиM62потенциально. Найдем его потенци# ал U (x, y):

dU 0 dx mg dy 0 dz

(x, y)

U (x, y) mg dy mgy.

(0, 0)

6.4.2.Соленоидальное поле

Пусть задано векторное поле F Pi Qj Rk . Линию, в каждой точке которой вектор поля является касательным, называют вектор#

ной линией.

 

 

Пусть в данном векторном поле дана замк#

 

нутая кривая L, через каждую точку которой

 

проходит векторная линия. Множество всех

 

векторных линий, проходящих через кривую L,

 

образует

поверхность, называемую векторной

 

трубкой.

 

 

 

Пусть

в любой

точке поля divF 0. Рас#

 

смотрим

область V,

ограниченную частью S

Рис. 6.3

300

Глава 6. Скалярное поле

векторной трубки и сечениями S1 и S2 . По формуле Остроградского – Гаусса

divF dV (F , n)dS (F , n)dS (F , n)dS,

V

S

S 1

S 2

где вектор нормали n во всех случаях направлен наружу. Поскольку divF 0, то

divF dV 0 (F , n)dS (F , n)dS (F , n)dS 0.

V

S

S 1

S 2

Заметим, что поскольку векторная трубка состоит из векторных линий, то вектор F лежит в плоскости, касательной к S и, следова# тельно (F , n) 0. Тогда (F , n)dS 0. Поэтому

S

 

 

 

 

M6

 

2

 

 

 

 

 

 

(F

,

 

n

)dS (F ,

n

)dS 0,

 

S 1

 

 

 

 

 

 

 

 

 

 

 

 

S

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

,

 

 

 

)dS,

 

(F

)dS (F

 

n

n

 

 

S 1

 

 

 

 

 

 

 

 

 

 

S 2

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

(F

 

 

n1)dS (F , n)dS,

 

S 1

 

 

 

 

 

 

 

 

 

 

S 2

где n1 n.

Таким образом, если векторное поле таково, что во всех его точках divF 0, то поток вектора через любое сечение векторной трубки име# ет одно и то же значение.

Если векторное поле — поле скоростей жидкости, и div v 0 (т. е. отсутствуют источники и стоки), то количество жидкости, протекаю# щее через любое поперечное сечение векторной трубки, имеет одно и то же значение.

Векторное поле, в любой точке которого divF 0, называется соле ноидальным или трубчатым.