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

pismennyy_d_t_konspekt_lekciy_po_vysshey_matematike_polnyy_k

.pdf
Скачиваний:
86
Добавлен:
12.05.2015
Размер:
11.99 Mб
Скачать

Расчленим этот поверхностный интеграл на три слагаемых, затем

сведем их вычисление к вычислению двойных интегралов. Нормаль к верхней стороне треугольника образует с осью Ох тупой угол, с осью

Оу - тупой, а с осью Oz - острый угол. (Единичный вектор данной

плоскости есть n = ±(~z + ~J - ~k); на верхней стороне cos1 > О,

поэтому надо выбрать знак «минус»; получим: cosa = -~, cos/3 = -~,

COS/ = ~-)

Итак, К= К1 + К2 + Кз. Находим К1, К2, К:3:

К1 = jj z dy dz = - jj

 

 

1

 

о

 

3

z dy dz = - j

dy j

 

z dz = ·" = 2,

 

s

 

 

 

 

вое

 

 

о

зu-з

 

К2

 

jJх

 

 

jj

 

 

2

 

 

о

 

= -

dx dz =

xdxdz

= J

 

J

dz = ···= 2,

 

 

 

 

xdx

 

 

 

 

S

 

 

 

АОС

 

 

О

 

Зх-6

 

 

 

 

 

 

 

 

 

 

 

 

- 2 -

 

 

 

 

 

 

 

 

 

 

 

6-Зх

 

 

Кз = jj

 

 

= jj

 

 

2

- 6 -

 

 

 

 

 

 

!

 

 

ydxdy

ydxdy

= J

ydy = ... = з1·

 

 

 

 

 

 

dx

 

 

 

 

S

 

 

АОВ

 

 

О

о

 

 

В результате имеем: К=~+ 2 +-! =

3~.

 

 

 

Пример 71.З.

Найти поток радиус-вектора f

z

 

 

 

 

 

 

 

 

 

 

 

 

 

через внешнюю сторону поверхности прямого кону­

са, вершина которого совпадает с точкой 0(0; О; О),

если известны радиус основания R и высота конуса

Н (см. рис. 275).

Q Решение:

jj

rпds + jj rпds = К1 + Kz.

К= fj rnds =

S

бок пов

осн

 

Очевидно, что К1 = О, т. к. пpn;f = О;

 

К2= jj rпds=H· jj ds=H·1ГR2 ,

осн

 

осн

т. к. пpn;f = Н. Итак, К= 1ГНR2 .

 

71.З. Дивергенция поля. Формула Остроградского-Гаусса

Важной характеристикой вРкторного поля (71.1)

у

Рис. 275

является так на­

зываемая дивергенция, характеризующая распределение и интенсив­

ность источников и стоков поля.

510

~Дuвергенцuеil (или расходимостью) векторного пол.я

а(М) = Р(х; у; z)l + Q(x; у; z)] + R(x; у; z)k

в точке М называется скаляр вида

дР + f!O + дR и обозначается

символом diva(M), т. е.

 

дх

ду

дz

 

 

 

 

. _

дР

дQ

дR

 

d1va(M)

= дх + ду

+ дz.

(71.6)

Отметим некоторые своi1ства дивергенции.

1.Если а - постоянный вектор, то div а= О.

2.div(c. а) =с. div а, где с= const.

3.div(a +Ь) = div а+ div Ь, т. е. дивергенция суммы двух векторных

функций равна сумме дивергенции слагаемых.

4. Если И - скалярная функция, а -- вектор, то

div(U ·а)= И· diva + agradU.

Эти свойства легко проверить, используя формулу (71.6). Дока­

жем, например, справедливость свойства 4.

Q Так как И · а = И · Р · l + И · Q · J+ И · R · k, то

div(U ·а)= ~(И· Р} +~(И· Q} +~(И·R) =

 

дх

 

ду

дz

 

 

= идР +Раи +идQ +QдИ +идR +RдИ =

дх

дх

ду

ду

дz

дz

= и(аР + дQ + дR) +Раи +Qau +Rau =

дх

ду

дz

дх

ду

дz

= и .div а + а .grad и.

Используя понятия потока и дивергенции векторного поля, запи­

шем известную в анализе (см. (58.9)} формулу Остроградского-Гаусса

lf Pdydz+Qdxdz+Rdxdy= JJJ(~= + ~~ + ~~)dv (71.7) s v

в так называемой векторной форме.

Рассматривая область V, ограниченную замкнутой поверхностью

S, в векторном поле (71.1), можно утверждать, что левая часть форму­

лы (71.7) есть поток вектора а через поверхность S; подынтегральная

функция правой части формулы есть дивергенция вектора а. Следова­

тельно, формулу (71.7) можно записать в виде

 

#an ds = jjj div а· dv

(71.8)

s v

(в котором она чаще всего и встречается).

511

где М0 -

liJ Формула Остроrрадского-Гаусса означает, что поток векторно­ го nоля 'Через замкнутую nоверхность S (в направлении внешнеiJ

нормали, т. е. изнутри) равен тройному интегралу от дuвергенции это­

го поля no обr.ему V, ограни'Ченному данноiJ поверхностью.

Используя формулу (71.8), можно дать другое определение дивер­

генции векторного поля ii(M) в точке М (не связанное с выбором ко­ ординатных осей).

По теореме о среднем для тройного интеграла (см. п. 54.1) имеем:

JJJ divii(M) dv = V · divii(M0 ), v

некоторая (средняя) точка области V. Тогда формулу (71.8)

можно переписать в виде /f andS = V · diva(M0 ). Отсюда s

div а(Мо) = ~ нan ds.

s

Пусть поверхность S стягивается в точку. Тогда V -t О, М0 -t М,

и мы получаем выражение для div ii(M) в точке М:

 

div а(М) = lim Vl fj an ds.

(71.9)

V-tO

s

~Дивергенцией векторного пол.я в точке М называется предел отношения потока поля через (замкнутую) поверхность S, окружа­

ющую точку М, к объему тела, ограниченного этой поверхностью, при

условии, что вся поверхность стягивается в точку М (V -t О).

Определение (71.9) дивергенции эквивалентно (можно показать) определению (71.6).

Как видно из определения, дивергенция векторного поля в точке является скалярной величиной. Она образует скалярное поле в данном

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

Исходя из физического смысла потока (~ычно условно счита­

ют, что ii(M) есть поле скоростей фиктивного стационарного потока несжимаемой жидкости), можно сказать, что: при divii(M) > О точ­

ка М представляет собой источник, откуда жидкость вытекает; при div ii(M) < О точка М есть сток, поглощающий жидкость. Как следует из равенства (71.9), величина div ii(M) характеризует мощность (интен­

сивность, плотность) источника или стока в точке М. В этом состоит физический смысл дивергенции.

Понятно, что если в объеме V, ограниченном замкнутой поверхно­

стью S, нет ни источников, ни стоков, то div ii = О.

Векторное поле, р каждой точке которого дивергенция поля равна

нулю, т. е. divii(M) =О, называется соленоидальным (или трубчатым).

512

Пример 11.4. Найти дивергенцию поля линейных скоросте&!: V

жидкости, вращающейся как твердое тело вокруг неподвижной оси с постоянной угловой скоростью iiJ.

О Решение: Примем ось вращения жидкости за ось Oz. Тогда, как по­

казано ранее (см. пример 69.2), V = -wyl + wx] +О· k. Имеем:

 

-

д

д

д

 

div V(M) = дх(-wy) + ду (wx) + дz (О) =О.

Поле V -

соленоидальное.

 

 

 

 

 

71.4. Циркуляция поля

 

 

Пусть

векторное

поле

образовано

z

 

вf'ктором (71.1). Возьмем н этом поле НЕ'­

которую замкнутую кривую L и выберем

на ней опр<'деленнос направлени<'.

Пусть f

= xz + yJ + zk -

радиус­

 

вектор точки М на кою уре L. Известно,

 

что вектор dт = dx·z+dy· J+dz·k напра­

 

влРн по касательной к кривой в на11ра­

 

влении ее обхода (см. рис. 276)

и ldrl =

Рис. 276

= dl, где dl -

дифференциал дуги кри­

 

вой (dl = ..j(dx)2 + (dy) 2 + (dz)2).

~Криволинейный интеграл по замкну~му контуру L от скалярного произведения вектора li на вектор dr, касательный к контуру L,

называется цuркул.яцuеil вектора ii вдоль L, т. е.

(71.10)

Рассмотрим различные формы записи циркуляции. Так как

а ·dт = jdтl ·пpdra = ат · dl = Рdx + Q dy + R dz,

где ат - проекция вектора li на касательную т, проведенную в напра­

влении обхода кривой L, то равенство (71.10)

можно записать в виде

С= f ат ·dl,

(71.11)

L

или

(71.12)

17 Конспектлекций по высшеи математике ПОJJныи курс

513

~Циркуляция С, записанная в виде (71.12) имеет простой физиче­

ский смысл: если кривая L расположена в силовом поле, то цирку­

ляция - это работа силы а(А1) поля при перемещении материальной

точки вдоль L (п.56.5).

Отметим, что вдоль замкнутых векторных линий циркуляция от­

лична от нуля, потому что в каждой точке векторной линии скалярное

произведение adr сохраняет знак: положительный, если направление

вектора а совпадает с направлением обхода векторной линии; отрица­

тельный - в противном случае.

При.мер 71. 5. Найти циркуляцию вектора ноля линейных ско­

ростей вращающегося тела (см. пример 69.2) V = -r.,;yz + r.,;xJ вдоль

замкнутой кривой L, лежащей в плоскости а, перпендикулярной оси

вращения.

О Решение: Будем считать, что направление нормали к плоскости а

совпадает с направлением оси Oz. Согласно формуле (71.12), имеем:

С = f -wy dx + wx dy = w f dx + хdy =

L

L

= 21.V ( ~ f dx + хdy) = 2w ·S,

 

 

L

где S - площадь поверхности, ограниченной кривой L (см. 56.17).

Заметим, чго если нормаль к поверхности S образует угол 'У с осью Oz, то циркуляция будет равна С= 2w · S · cos1; с изменением угла 'У

величина С изменяется.

 

 

8

 

 

 

z

 

Пример 71. 6. Вычислить циркуляцию век-

 

1

С

 

 

 

 

торного поля

 

 

 

 

а= (х -

2z)z + (х +Зу+ z)J + (5х + y)k

 

 

 

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

 

 

 

A(l;O;O), В(О; 1;0), С(О;О; 1) (см. рис. 277).

 

 

 

О Решение: Согласно формуле (71.12), имеем:

 

Рис. 277

С = f (х -

2z) dx + (х + Зу + z) dy + (5х + у) dz = j

+ j

+ J

L

АВ

ВС

СА

На отрезке АВ: х +у= 1, z =О, следовательно,

 

 

 

J

о

 

з

 

= J- О) tl,x + + З - Зх +О)· (-dx)

+О=

 

АВ

1

 

 

 

514

На отрезке ВС: у+ z = 1, х =О, следовательно,

 

f

о

2 + 2у) ·О+ (О+ Зу+ 1 - у)dy + (О+ у) · (-dy) =

3

= f -

--z:

вс

1

 

 

На отрезке СА: х + z = 1, у= О, следовательно,

 

 

 

1

 

 

J = j(x-2+2x)dx+0-1(5x+O)·(-dx)=-3.

 

 

СА

О

 

Следовательно,

 

 

С= f = ! + ! + ! = ~+ (-~) + (-3) = -3.

 

АВСА АВ ВС СА

71.5. Ротор поля. Формула Стокса

~Ротором (или вихрем) векторного пол.я

а= Р(х; у; z)z + Q(x; у; z)J + R(x; у; z)k

называется вектор, обозначаемый rota(M) и определяемый формулой

_

(дR

-

дQ),i

(дР

-

8R)..,.

+

(дQ

8P)-

(71.13)

rot а(М) =

-ду

-дz

+ -дz

-дх J

-дх

- -ду

k.

 

Формулу (71.13) можно записать с помощью символического опре­

делителя в виде, удобном для запоминания:

i

j

k

rota(M) = lx

а

а

ду

дz

р

Q

R

 

Отметим некоторые своil.ства ротора.

1.Если а - постоянный вектор, то rota =О.

2.rot(с · а) =с · rot а, где с =const.

3.rot(ii + Б) = rot а+ rot Б, т. е. ротор суммы двух векторов равен

сумме роторов слагаемых.

4. Если И - скалярная функция, а а(М) - векторная, то

rot(U ·а)= И rot а+ grad Их а.

Эти свойства легко проверить, используя формулу (71.13). Пока­

жем, например, справедливость свойства 3:

515

Используя понятия ротора и циркуляции, векторного поля, за­

пишем известную в математическом анализе (см. п. 58.4) формулу

Стокса:

f Pdx+Qdy+Rdz= JJ(~R _ ~~)dydz+

L S у

дР

дR)

dx dz +

(дQ

-

дР)

 

+ (-дz

- -д:r

-ОХ

-ду

cl.r cly. (71.14)

Левая часть формулы (71.14) представляет собой циркуляцию вектора

а по контуру L, т. i:>. f

Pdx+Qdy+Rdz =

f aтdl (гм. (71.11)). Интеграл

L

 

L

в правой части формулы (71.14) предгтанляст собой поток вектора rot а через поверхность S, ограниченную контуром L (см. (71.3)), т. е.

дR дQ)

(дР дR)

(дQ дР)

dx dy =

!.!(-ду - -дz

dy dz + -дz

- -дх

dx dz + -дх - -ду

s

 

 

= JJrotniids.

z

 

 

 

 

 

s

 

 

Следовательно, формулу Стокса можно запи-

 

сать в виде

 

 

 

 

f aтdl = jj rot11 ii ds.

(71.15)

0'~-----

L

S

 

 

 

 

уТакое представление формулы Стокса на­

х

Рис. 278

зывают ее векторноiJ, формоiJ.. В этой формуле

 

 

 

положительное направление на контуре L и вы­

бор стороны у поверхности S согласованы между собой так же, как в

теореме Стокса.

 

liJ

Формула (71.15) показывает, что циркуJLЯ'ЦUЯ вектора ii ваоль ,Jа­

мкнутого контура L равна потоку ротора этого вектора ii -ч,ерез поверхность S, лежащущ в поле вектора ii и ограниченную контуром

L (натянутую на контур) (см. рис. 278).

516

Используя формулу (71.14), можно дать другое определение рото­

ра поля, эквивалентное первому и не зависящее от выбора координат­

ной системы.

Для этого применим формулу Стокса (71.15) для достаточно малой

плоской площадки S с контуром L, содержащей точку М.

По теореме о среднем для поверхностного интеграла (п. 57.1, свой­

ство 7) имеем:

jj rotn ads = rotn а(Мо) ·S, s

где М0 - некоторая (средняя) точка площа.л;ки S (см. рис. 279).

Тогца формулу (71.15) можно записать в виде

f a,,.dl = rotn а(Р) ·S.

L

Отсюда:

rot" а(Р) = ~ f а,,.dl.

[,

Пусть контур L стягиваRтся в точку М. Тогда

Рис 279

М0 --+ М, а S --+ О. Перейдя к прРделу, получаем:

 

rotn а(М) = lim - f а,,.dl.

 

1

 

8-+О 8

 

L

 

~Ротором вектора а в точке М называется вектор, проекция

которого на каждое направление равна пределу отношения цирку­

ляции вектора а по контуру L плоской площадки S, перпендикулярной

этому направлению, к площади этой площадки.

Как видно из определения, ротор вектора а(М) есть векторная

величина, образующая собственное векторное поле.

Дадим физическое истолкование понятия ротора векторного поля. Найдем ротор поля линейных скоростей твердого тела, вращающегося

вокруг оси Oz с постоянной угловой скоростью (пример 69.2) w, т. е.

ротор вектора V = -w ·у · i + w ·х · J.

По определению ротора

 

 

i

j

k

 

 

rota(M) =

lx

ly

lz

 

 

 

 

-yUJ

XUJ

0

 

 

(

o_

д(wx))l-(o-8(-yw))J+ (8(xw) _ 8(-yw))k=

 

8z

 

8z

дх

 

 

 

 

 

= О -

О + 2w · k = 2CiJ.

Ротор этого поля направлен параллельно оси вращения, его мо­

дуль равен удвоенной угловой скорости вращения.

517

С точностью до числового множителя ротор поля скоростей V

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

связано само название «ротор» (лат. «вращатель»).

Заме'Чание. Из определения (71.13) ротора вытекает, что напра­

вление ротораэто направление, вокруг :к;9торого циркуляция имеет

наибольшее значение (плотность) по сравнению с циркуляцией вокруг любого направления, не совпадающего с нормалью к площадке S.

Так что связь между ротором и циркуляцией аналогична связи

между градиентом и производной по направлению (см. п. 70.3).

§ 72. ОПЕРАТОР ГАМИЛЬТОНА 72.1. Векторные дифференциальные операции первого

порядка

Основными дифференциальными операциями (действиями) над скалярным полем И и векторным полем а являю1ся grad И, div а, rot а.

Действия взятия градиента, дивергенции и ротора называются вектор­

ными операциями первого порядка (в них участвуют только первые производные).

Эти операции удобно записывать с помощью так называемого опе­

ратора Гамильтона

 

 

д,

д-

д­

\7 = -'/,

+ - J

+ - k .

дх

ду

дz

Этот символический вектор называют также оператором \7 (чита­ ется «набла»); он приобретает определенный смысл лишь в комбинации

со ска,.,'lярными или векторными функциями. Символическое «умноже­

ние» вектора У' на скаляр И или вектор а производится по обычным

правилам векторной алгебры, а «умножение» символов lx, gy, lz на

величины И, Р, Q, R понимают как взятие соответствующей частной

производной от этих величин.

Применяя оператор Гамильтона, получим дифференциальные опе­

рации первого порядка:

1. У'И = (д_l + _Q_J + д_JС) ·И= дИ;, + дИJ + дИJе = gradU.

 

дх

ду

дz

дх ду

 

дz

 

2. У'а =

(д_l+_g_J+_o_k)·(P·l+Q·J+R-k) =

дР+00..+дR = diva.

 

&

~

&

 

&

~

&

 

i

J

k

 

 

 

 

t"'7

д

д

д

=rota.

 

 

 

3.vxa= дж

ду

дz

 

 

 

 

р

Q

R

 

 

 

 

Оператор ГамильтонА применяется для записи и других операций и для вывода различных формул в теории поля. При действиях с ним

518

ведение двух одинаковых векторов равно нулю (нуль-вектор). Это
означает, что поле градиента есть поле безвихревое.
519
2. rot grad И = \1 х (\!И)
= (\1 х \l)U = О, так как векторное произ­
(который тоже называют оператором Лапласа).
liJ

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

ренцирования.

В частности, производная по направлению (70.2) может быть за-

писана в виде

аи

= \!И . е: = (е:. \1) . и,

 

 

ал

где е = (cosa:;cos,В;cos1).

72.2. Векторные дифференциальные операции второго

порядка

После применения оператора Гамильтона к скалярному или век­

торному полю получае1ся новое поле, к которому можно снова приме­

нить этот оператор. в результате получаются оифференчиалъные опе­ рации второго порядка. Нетрудно убедиться, что ИМ<'<'ТСЯ лишь пять дифференциальных опера1~ий второго порядка: div grad И, rot gra<l И,

grad div а, div rot а, rot rot а.

(Понятно, что операция div div а, например, не имеет смысла: div а - скаляр, говорить о дивергенции скаляра, т. е. о div div а, бес­

смысленно.)

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

ственно за оператором.

 

 

liJ

2

2

2

1. div grad И = \l(\lU) = (\1 ·\l)U = (tx2

+ :у2

+ tz2) ·И =

 

-~:Ч + ~:Ч + ~:Ч. Правая часть этого равенства называется

оператором Лапласа скалярной функции И и обозначается дИ. Таким

образом,

а2 и

а2 и

а2 и

 

 

(72.1)

 

divgradU = дИ = --2

+ --2

+ -- 2 ·

 

8z

 

Дифференциальное уравнение Лапласа дИ = О играет важную роль в различных разделах математической физики. Решениями уравнения Лапласа являются так называемые гар.мони'Ческие функции.

За.ме'Чание. К равенству (72.1) можно прийти, введя в рассмотрение

скалярный оператор дельта:

а2 а2 а2

д = \1 . \1 = 8х2 + 2 + 8z2