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

ovta-posibnyk

.pdf
Скачиваний:
61
Добавлен:
19.03.2015
Размер:
2 Мб
Скачать
☆

Рівняння векторних ліній (див. (6.1)) для цього випадку мають

вигляд

dx

dy

dz ,

 

 

 

 

x

y

z

 

інтегрування їх дає систему двох рівнянь y C1x ,

z C2 y за

умови r

0 , які у просторі визначають два набори площин із

виключеним початком координат. Перетин їх дає систему променів, що виходять із початку координат. Вони і визначають сім'ю векторних ліній заданого поля.

2.2. Похідна від векторної функції за напрямком

Як і у випадку скалярного поля, уведемо поняття похідної за

напрямком. Напрямок, у якому нас цікавить зміна неперервно

диференційовної функції

 

 

 

 

, характеризуватимемо одиничним

 

 

A

вектором

 

. Розглянемо

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

до

 

 

 

s

 

 

точку r , близьку

 

r , таку, що

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

r s . Тоді границю відношення

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A(r s ) A(r )

 

 

A

 

 

 

 

lim

A(r ) A(r )

 

lim

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

називатимемо похідною векторної функції у точці r

за напрям-

ком s

. Вважаючи

 

 

малим і використовуючи розвинення в ряд

Тейлора

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s , y

s , z s )

 

 

 

 

 

 

 

 

A(r

s ) A(x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 A

 

 

 

 

 

 

 

1

 

A s

1

A

 

 

 

 

 

 

 

 

 

 

A( x, y, z)

 

 

s

 

 

s o( ) ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1! x

 

 

 

одержуємо

 

 

 

 

 

 

 

 

 

1

 

1! y

2

 

 

1! z

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

A

 

 

A

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

s1

 

 

s2

 

 

 

s3

 

 

 

 

s1

 

 

 

s2

 

 

 

s3

 

 

A

(s

) A .

 

s

x

y

 

z

 

 

x

 

y

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Зазначимо, що напрямок необов'язково характеризувати оди-

ничним вектором

 

s

. Можна говорити про похідну від функції

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Узагальнюючи останню

A за вектором v довільної довжини.

 

 

формулу, отримаємо

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v1

 

 

 

v2

 

 

 

 

 

v3

 

 

 

A (v

) A .

 

 

 

 

 

 

 

 

 

 

x

y

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

101

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Приклад 8. Знайти (a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

)r .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Розв'язання. (a

)r ax

 

 

ay

 

 

 

 

az

 

 

 

r

 

 

x

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ax

 

ay

 

 

 

 

az

 

xi

yj

zk

axi

ay j

azk

a.

x

 

y

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

характеризує середовище, точки якого

Нехай поле A(x, y, z,t)

рухаються

зі швидкостями

 

x

,

y

,

z

 

 

 

 

v

 

t

t

 

t

. Знайдемо швид-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

кість зміни поля з часом у заданій точці середовища

dA A A x A y A z A (v ) A . dt t x t y t z t t

Доданок At виражає швидкість зміни величини A із часом у

заданій точці простору і називається локальною похідною. Локальна похідна характеризує нестаціонарність поля. Доданок

(v ) A утворюється за рахунок зміни координат точки внаслі-

док її руху (конвекції), тому його називають конвективною похідною. Конвективна похідна характеризує просторову не-

однорідність поля в даний момент часу.

 

 

 

 

 

 

 

Застосовуючи останню формулу

до

векторного

 

поля

 

dv

 

v

 

 

 

 

швидкостей рідини v x, y, z,t , маємо

dt

 

t

(v

)v

. Як

видно, повне прискорення складається із двох доданків: локального прискорення v та переносного (конвективного)

прискорення (v )v . t

102

2.3. Інваріантні характеристики векторних полів

 

Розглянемо диференціал

векторного стаціонарного поля

A x, y, z

 

 

 

 

 

 

 

 

 

 

 

 

 

dA

 

 

 

dA

dr

 

 

або за компонентами

 

dr

 

 

 

 

Ai

 

 

 

 

dA

 

dx

 

,

 

 

x j

 

 

i

 

 

j

 

тобто в результаті згортки деякого геометричного об'єкта, який

задається 9 компонентами

Ai

,

i, j

 

із вектором, координа-

1,3,

 

 

x j

 

 

 

ти якого дорівнюють {dxi} , утворився вектор із координатами {dAi}. Із наслідку з оберненої тензорної ознаки випливає, що

величини

Ai

,

i, j

 

утворюють тензор другого рангу. У

1,3,

 

 

x j

 

 

 

результаті згортки тензора другого рангу за парою індексів утворюється скалярна величина

Ai ( A) div A ,

xi

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

 

Приклад 9. Обчислити

div r ,

де r ( x, y, z) – радіус-вектор

точки.

 

 

 

 

 

 

 

 

 

 

x

y

z

 

 

Розв'язання. div r

( r )

x

y

z 3 .

 

 

Приклад 10. Обчислити div

ra , де a – сталий вектор,

r =| r | ,

( x, y, z) .

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

103

 

 

 

 

 

 

rax

 

 

rax

 

rax

 

Розв'язання. div ra ra

 

 

 

 

 

 

 

 

y

z

 

 

 

 

 

 

x

 

 

 

 

 

 

 

x

 

y

 

 

 

z

 

 

 

 

 

 

 

ax

ay

az

 

a

r .

 

 

 

 

r

 

r

 

 

 

 

r

 

 

 

 

 

r

 

 

 

 

 

Оскільки довільний тензор можна представити у вигляді суми симетричного sij та антисиметричного aij відносно переста-

новки індексів тензорів, то справедливе таке розкладання:

 

 

 

 

 

 

 

 

 

Ai

 

s

a ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x j

ij

 

 

ij

 

 

 

 

 

 

де

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

A

 

 

Aj

 

 

 

1

 

A

 

Aj

s

 

 

 

 

 

 

i

 

 

 

 

, a

 

 

 

 

i

 

 

.

 

2

x

 

 

x

2

x

 

x

ij

 

 

j

 

 

 

ij

 

 

j

 

 

 

Ai

 

 

 

 

 

 

i

 

 

 

 

 

 

 

i

 

Величини

 

 

утворюють симетричний тензор лише тоді,

x j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

потенціальне (див. п. 2.4). Другий

A(r ) –

інваріант (сума головних мінорів відповідної матриці) антисиметричного тензора другого рангу aij , компоненти якого можна

представити у вигляді (див. (5.6)) aij eijk k , дорівнює 2 . У

 

 

 

 

 

 

 

 

 

 

 

 

1

 

A

 

 

Aj

 

 

 

нашому

 

 

випадку

 

a

 

 

 

 

i

 

 

 

 

 

і

відповідно

 

 

 

2

x

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

ij

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

1

 

A

 

Aj

 

1

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

[ A]

 

 

 

 

[rot A]

,

де

i, j,k утво-

2

x

 

 

2

 

2

k

 

 

j

 

x

 

 

 

k

 

 

 

 

k

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

є позначенням векторної ве-

рюють циклічний порядок, а rot A

личини, яку називають ротором (вихором) векторного поля (інваріантне означення див. п. 2.6) і яка змінюється від точки до точ-

ки. Можемо говорити, що rot A є вектором, який характеризує відхилення заданого векторного поля від потенціального (для

потенціального поля rot A 0 ).

Отже, div A та rot A є інваріантними диференціальними характеристиками векторних полів у кожній точці простору.

104

ного поля.

Приклад 11. Обчислити rot r , де r

( x, y, z) .

 

 

 

 

 

 

 

i

 

 

j

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Розв'язання.

rot r

r

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

y

 

 

z

 

 

 

 

 

 

 

 

 

 

 

x

 

 

y

 

 

 

z

 

 

 

 

 

 

 

 

 

z

 

y

z

 

x

 

 

 

 

y

 

x

0 .

i

 

y

 

j

x

z

 

k

 

x

 

 

 

 

z

 

 

 

 

 

 

 

y

 

Приклад 12. Обчислити rot ra , де a

– сталий вектор, r =| r | .

Розв'язання. Знайдемо x -компоненту rot ra :

[r a]

 

 

 

 

 

(ra

 

)

 

(ray )

 

 

y

 

 

z

 

 

rot(ra)

x

 

 

 

z

 

 

 

a

z

 

a

y

 

 

 

x .

 

y

 

z

r

r

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тому rot ra

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.4.Потенціальне векторне поле. Лінійний інтеграл від вектора

Означення. Векторне поле A(r ) 3 Ai ( x, y, z)ei , задане у про-

i 1

сторовій області G , називається потенціальним, якщо існує така скалярна функція ( x, y, z) , що для всіх точок області G вико-

нується рівність:

A( x, y, z) grad ( x, y, z) .

Зауваження. Для силових полів зазвичай

A(x, y, z) grad ( x, y, z) .

Тому далі враховуватимемо знак мінус перед градієнтом.

Функція ( x, y, z) називається скалярним потенціалом вектор-

Щоб з'ясувати, чи є поле A потенціальним, необхідно перевірити виконання таких умов потенціальності векторного поля:

105

 

 

A1

A2 ,

A1

A3 ,

A3

A2 ,

 

 

y

x

z

x

y

z

або

, які є наслідком незалежності другої мішаної похід-

rot A 0

ної скалярної функції ( x, y, z)

від порядку диференціювання.

Векторними лініями потенціального векторного поля є лінії поля градієнта відповідного скалярного потенціалу (r ) , тобто лінії най-

швидшої зміни скалярної функції (r ) , що перпендикулярні до поверхоньрівняскалярногопотенціалу (r ) const (див. § 1).

Важливу характеристику векторного поля дає лінійний інтеграл від вектора – звичайний криволінійний інтеграл другого роду – інтеграл уздовж напрямленого контуру від скалярного

добутку вектора A(r ) на векторний елемент довжини:

A dr A1dx A2dy A3dz ,

де вектор dr

C

C

має складові dx,dy,dz , а диференціали координат не

незалежні, а є приростами відповідних координат уздовж контуру. Такими інтегралами виражається, наприклад, робота сили A(r ) над матеріальною точкою, яка рухається за заданою траєкторією від точки M0 доточки M1 . Якщовектор A(r ) – потенціальний, то

A1dx A2dy A3dz x dx x dy x dz d

є повним диференціалом скалярної функції. Обчислюючи інтеграл, отримаємо

 

 

M1

 

M1

A dr

 

A dr

d (M0 ) (M1) .

C

 

M0

 

M0

У цьому випадку контурний інтеграл визначається положенням початкової та кінцевої точок інтегрування (рис. 6.1) і не залежить

від форми контуру.

Інтегруючи вздовж замкненого контуру

(рис. 6.2) і враховуючи, що

M1

 

 

M1

 

A dr

d (M0 ) (M1)

M0

 

 

M0

 

 

 

106

значення.
Криволінійний інтеграл другого роду по замкненій кусково-гладкій оріє-
нтованій кривій C
A dr
C
називають циркуляцією векторного поля
у заданому напрямку.
(x, y, z) може набувати іншого
неоднозначна, то після обходу контуру і повернення у точку M0 , функ-
ція
нозначна.

та

 

M0

 

 

 

 

 

 

A dr (M1) (M0 ) ,

 

M1

 

 

 

 

матимемо

 

 

 

 

 

M1

 

M0

 

A dr

 

A dr

 

A dr 0 ,

C

 

M0

 

M1

 

якщо скалярна функція ( x, y, z) одЯкщо функція (x, y, z)

M1

M0

Рис. 6.1. Лінійний інтеграл від потенціального вектора не залежить від шляху інтегрування

M1

M0

Рис. 6.2. Циркуляція потенціального вектора дорівнює нулю

A вздовж контуру C

Ми показали, що циркуляція потенціального вектора A вздовж довільного замкненого контуру C дорівнює нулю за умови, що відповідний скалярний потенціал є однозначною функцією. Зауважимо, що такої важливої властивості довільний вектор не має.

 

Приклад 13. Показати, що напруженість електростатичного по-

ля E , створеного системою зарядів, є потенціальним вектором.

 

Розв'язання. Вектор E

можна записати у вигляді E grad ,

де

(r )

1

 

 

ei

 

 

– однозначна функція. Інтеграл E dr

 

 

 

 

4 0

i

 

r ri

 

 

C

 

 

 

 

має прямий фізичний смисл – він виражає роботу, яку виконує поле при переміщенні одиничного позитивного заряду вздовж кривої

107

від проекції вектора

C . Умова E dr 0 означає те, що робота з переміщення заряду

C

вздовж замкненоїкривої C дорівнює нулю.

Загальних правил обчислення криволінійних інтегралів немає. Їх намагаються звести до обчислення визначених інтегралів, належним чином параметризувавши контур (приклади див. у розд. 8). Якщо контур інтегрування замкнений, можна скористатися формулою Стокса (див. п. 2.6) і звести обчислення до знаходження поверхневого інтегралаі подвійного інтеграла Рімана.

2.5. Потік векторного поля через поверхню

Однією із важливих характеристик векторного поля є потік векторного поля через поверхню. Потоком векторного поля

A(r ) через гладку або кусково-гладку орієнтовану поверхню S називається поверхневий інтеграл другого роду по поверхні S A(r ) на нормаль n( x, y, z) до обраної сто-

рони цієї поверхні:

 

 

 

 

AndS A dS

( A n)dS .

S

S

 

S

Інтерпретацію потоку вектора легко дати у випадку, коли вектором A(r ) є вектор швидкості рухомої рідини. Тоді потік век-

тора A(r ) – це кількість (об'єм) рідини, що проходить за одини-

цю часу через орієнтовану поверхню S .

Існує декілька способів обчислення потоку вектора, використовуючи, зокрема:

1)проектування на одну із координатних площин;

2)проектування на всі три координатні площини;

3)перехід до криволінійних координат на поверхні.

Якщо поверхня, по якій ведеться інтегрування, є замкненою, то обчислення поверхневого інтеграла можна звести до обчислення об'ємного інтеграла, використовуючи теорему Остроград- ського–Гаусса, про яку йтиметься пізніше (див. п. 2.6).

108

2.6. Дивергенція та ротор. Інтегральні теореми

Щоб з'ясувати математичний і фізичний зміст ротора (вихору) та дивергенції векторного поля, розглянемо їх інваріантні

визначення – менш формальні, ніж через оператор набла:

 

 

 

 

 

div A A ,

rot A A ,

але більш наочні. Вони, можливо, є складнішими, але не залежать від вибору системи координат. Почнемо з ротора.

Виберемо точку M , у якій ми хочемо визначити ротор векторного поля. Задамо деякий напрямок одиничним вектором n .

Побудуємо маленьку плоску площадку

 

 

 

S , що проходить через точку M пе-

 

M

n

рпендикулярно до n . Означимо на-

 

S

 

прямок обходу контуру C , що обме-

r

 

 

 

жує цю площадку (рис. 6.3), як такий,

 

 

 

що узгоджується із напрямком n пра-

Рис. 6.3. Ілюстрація

вилом правого гвинта: якщо дивитися

з вістря вектора n , то обхід контуру

правила правого гвинта

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

 

 

 

 

 

 

Означення. Проекція вектора ротора (вихору) rot A на напря-

мок n у точці M – це віднесена до одиниці площі циркуляція век-

тора A по замкненому нескінченно малому контуру, що обмежує

площадку знормаллю n , в якій лежитьдана точка:

 

 

 

 

 

 

 

 

 

 

A dr

 

 

(rot A)n n

rot A lim

C

 

.

 

 

S

 

 

S 0

 

 

 

Ця величина показує, наскільки векторне поле

A r відріз-

няється від потенціального.

Можна переконатися у тому, що це визначення збігається з

означенням rot A A у декартових координатах. Для цього обчислимопроекціїротора, користуючисьінваріантнимозначенням і вибравши площадку у вигляді прямокутника зі сторонами, паралельними координатним осям. Наприклад, виберемо прямоку-

109

тну

площадку

S x y

(x,y+ y,z)

(x + x,y + y,z)

(рис. 6.4) та n уздовж осі Oz .

 

 

 

 

 

 

 

 

 

 

 

 

Враховуючи малі розміри пло-

 

 

 

 

 

 

щадки та використовуючи при

 

 

 

 

 

 

 

 

 

 

 

 

обчисленні контурного інтегра-

 

 

 

 

 

 

(x,y,z)

(x + x,y,z)

ла (рис. 6.4) теорему про середнє

Рис. 6.4. Обхід

на кожній із ділянок, отримаємо

 

 

 

прямокутної площадки

 

 

A dr Ax (x, y, z) x Ay (x x, y, z) y

C

 

Ax (x, y y, z) x Ay (x, y, z) y

Ay ( x x, y, z) Ay (x, y, z) y

A ( x, y, z) A (x, y y, z) x

Ay

 

A

 

 

 

 

 

x

 

S ,

 

x

x

x

 

y

 

 

 

 

 

 

 

де x, y – значення координат у деяких точках відповідних відрі-

зків контуру. Враховано також, що обхід протилежних сторін здійснюється у протилежному напрямку, а при стягуванні площадки S у точку M усі координати набувають значень, що відповідають цій точці. Підставивши останній результат в означення і здійснивши граничний перехід, отримаємо вираз для

 

 

 

Ay

 

A

 

A

 

 

 

x

у декартових координатах. Аналогічно

x

 

z

 

 

y

 

знаходять інші проекції вихору.

Таким чином, вихор вектора в заданому напрямку не залежить від вибору системи координат і дорівнює поверхневій густині циркуляції цього вектора по контуру площадки, що перпендикулярна до цього напрямку. Вихор вектора відрізнятиметься

від нуля, якщо векторні лінії поля A закручуються, – мають замкнену або спіралеподібну форму.

Розглянемо тверде тіло, яке обертається навколо осі Oz із сталою кутовою швидкістю. Векторне поле швидкостей точок цього тіла подамо у вигляді

110

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]