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

Tom_2

.pdf
Скачиваний:
59
Добавлен:
18.02.2016
Размер:
3.2 Mб
Скачать

2

 

y2

1

1

2

 

 

 

1

 

 

1

 

 

= ò(2 - 3y +

)dy + ò(2 - 2x)dx - ò(2 - 5x + 2x2 )dx = -

+1-

 

=

 

.

 

3

6

6

0

2

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 2. Вычислить поток векторного поля a = x

 

 

- y

 

+ z

 

 

i

 

j

k

через внешнюю сторону части S конуса z2 = x2 + y2 ,

лежащую выше

плоскости z = 0 и внутри цилиндра x2 + y2 = C2 (C > 0) .

 

 

 

 

 

 

 

 

 

 

 

 

Решение.

Часть

 

S

 

конуса

 

 

z2 = x2 + y2 ,

 

 

лежащую

выше

плоскости z = 0 , зададим явно функцией z =

 

 

x

2

+ y

2

 

 

 

 

 

 

 

 

 

 

 

 

%

 

 

 

 

 

 

 

 

, где (x; y)Î D ,

%

 

 

 

 

2

 

 

 

2

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

+ y < C } . Нормаль к внешней стороне части S конуса

D = {(x; y)

 

 

образует тупой угол с ортом

 

 

 

оси Oz.

 

 

Поэтому это отрицательная

k

сторона S. Предварительно вычислив

¢

=

 

 

 

x

 

 

 

,

 

 

 

¢

 

=

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

zx

 

x2 + y2

 

 

zy

 

 

x2 + y2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x; y) D ,

 

 

где

 

 

 

D = {(x; y)

 

0 < x2 + y2 < C2} ,

 

по

 

формуле

(3)

 

 

 

 

 

 

 

 

находим поток

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

æ

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ö

 

 

 

 

 

 

 

 

 

 

 

 

 

П = -òò

ç

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

÷

 

 

 

 

 

 

 

 

 

 

 

 

 

ç -x ×

 

 

 

 

 

 

 

 

 

 

+ y ×

 

 

 

 

 

 

+

 

x

 

+ y

 

 

 

÷ dxdy

=

 

 

 

 

 

 

 

x

2

 

+ y

2

x

2

+ y

2

 

 

 

 

 

 

 

 

 

 

 

 

 

D è

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ø

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2y2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

 

 

 

 

 

 

 

 

 

 

 

C3

 

 

 

 

 

 

= -òò

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dxdy = -2 ò sin2 ϕ dϕ ò ρ2 d ρ = -

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

x2 + y2

 

 

 

 

 

 

 

 

 

 

D

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример

 

 

3.

 

 

 

 

 

Вычислить

поток

 

П

 

векторного

 

поля

a = xy

i

+ yz

j

+ zx

k

 

 

 

 

через левую сторону поверхности

S

 

(сторону

поверхности

 

будем

 

 

считать

левой,

 

если

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

заданной

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos(n, i ) £ 0 ),

 

параметрически:

 

x = u3 - v3,

 

y = u2v, z = uv2 , "(u;v)ÎW

,

где

 

 

 

 

W

замыкание W = {(u;v)

 

 

 

 

0 < u <1, 0 < v <1} .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение. Чтобы использовать формулу (4) , вычислим

 

 

 

 

 

PQ

=x y u u x y v v

R

uz . z v

Так как P = xy = u5v - u2v4 , Q = yz = u3v3, R = zx = u4v2 - uv5,

будем иметь

154

 

 

u5v - u2v4 u3v3 u4v2 - uv5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D =

 

 

 

 

3u2

2uv

 

 

v2

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-3v2

u2

 

 

2uv

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 4u7v3 - 4u4v6 + 3u8v2 - 3u5v5 - 3u3v7 + 6u5v5 - 6u2v8 - u7v3 +

 

 

 

+u4v6 - 6u6v4 = 3u8v2 + 3u7v3 - 6u6v4 + 3u5v5 - 3u4v6 - 3u3v7 - 6u2v8.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

D( y, z)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Поскольку cos(n, i )

=

n

, где

 

n

=

 

A

+ B

 

+ C

 

, A =

D(u,v)

 

,

 

D(z, x)

 

 

 

D(x, y)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A = 3u

2

v

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B =

 

 

 

,

C =

 

 

 

,

имеем cos(n, i )

> 0, так как

 

 

> 0

D(u,v)

D(u,v)

 

 

(u,v) W.

Следовательно, по формуле (4) поток

11

П= -òòDdu dv = òduò(-3u8v2 - 3u7v3 + 6u6v4 - 3u5v5 +3u4v6 + 3u3v7 +

W

 

 

0 0

 

 

 

 

 

 

+6u2v8 )dv =

1 æ

-u8 - 3 u7 + 6 u6 - 1 u5 + 3 u4 + 3 u3 +

2 u2

ödu =

 

ò

ç

4

5

2

7

8

3

÷

 

è

ø

 

0

 

 

 

 

 

 

 

 

= - 19 - 323 + 356 - 121 + 353 + 323 + 92 = 1260359 .

20. Поток вектора через замкнутую поверхность. Формула Остроградского. Если поверхность S замкнута, то в качестве нормали к ней возьмем внешнюю нормаль.

Теорема 1 (Остроградского). Если в некоторой пространственной области V координаты P, Q, R вектора a

непрерывны и имеют непрерывные частные производные

Px , Qy , Rz , то поток вектора a через замкнутую кусочно-гладкую

поверхность S, расположенную в области V, равен тройному интегралу по области V, ограниченной этой поверхностью, то есть

П =

òò

(a, n)ds =

òòò

æ P

+

Q

+

R ödxdydz .

(5)

 

 

è x

 

y

 

z ø

 

 

 

 

 

ç

 

 

 

÷

 

 

S

 

V

 

 

 

 

 

 

Доказательство. Пусть замкнутая поверхность S ограничена поверхностями S1 и S2 , являющимися графиками функций z1 = f1(x, y) , z2 = f2 (x, y) соответственно, и цилиндрической

155

поверхностью S3 с направляющей − границей области D (проекции поверхности S на плоскость Oxy ) и образующей, параллельной оси

Oz (рис. 2). Пользуясь формулой (4.4.1) приведения тройного интеграла к повторному, будем иметь:

 

R(x, y,z)dxdydz = òòdxdy

f (x,y)

R(x, y,z)dz =òòR(x, y, f2(x, y))dxdy -

òòò

2 ò

V

z

D

f (x,y)

z

D

 

 

 

1

 

 

-òòR(x, y, f1(x, y))dxdy = òòR(x, y,z)dxdy + òòR(x, y,z)dxdy.

D

 

S2

 

 

S1

z

S2 n

f2(x,y)

S3

S1 n f1(x,y)

0

y

D

x

Рис. 2

согласно формуле (5.9.9),

(6)

При выполнении преобразований (6) использованы формулы (см. п. 5.9.20), которые связывают поверхностный интеграл второго рода с двойным интегралом и тот факт, что на поверхности S1

вектор нормали образует с осью Oz тупой угол, а на поверхности S2 − острый.

Поскольку на поверхности S3 вектор нормали есть n = (α; β; 0) , то

òòR(x, y, z)dx dy = òò(0 ×α + 0× β + R ×0)ds = 0.

S3

S3

Прибавляя к поверхностным интегралам в формуле (6) равный нулю интеграл òòR(x, y, z) dx dy, будем иметь

 

S3

 

 

 

òòò R(x, y, z)

dx dy dz = òòR(x, y, z)dx dy +òòR(x, y, z)dx dy +

 

V

z

S2

S1

(6΄)

 

+òòR(x, y, z)dx dy =òòR(x, y, z)dx dy.

S3

 

S

 

 

 

 

Аналогично доказывается, что

 

 

 

 

òòò

P

dxdydz = òòPdydz , òòò

Q

dxdydz = òòQdzdx .

(7)

x

y

V

S

V

S

 

 

 

 

Складывая формулы (6΄) и (7), получим формулу (5).

Формула (5) называется формулой Остроградского.

156

Пусть a – поле скоростей движения несжимаемой жидкости. Если П = 0 , то это означает, что из объема V, ограниченного поверхностью S, вытекает столько жидкости, сколько и прибывает в него. Если П > 0 , то жидкости вытекает больше, чем прибывает. Это

говорит о том, что внутри объема V

имеется источник − место, где

жидкость появляется. Если П < 0 ,

то жидкости вытекает из тела

меньше, чем прибывает. В этом случае говорят, что внутри V имеется

сток − место стекания жидкости.

 

 

 

Пример 4.

Вычислить

поток

векторного

поля

a = 2xyi

 

- y2

 

+ z3

 

через внешнюю сторону замкнутой поверхности

 

j

k

S: {x2 + y2 + z2 = 2Rz, x2 + y2 = z3} .

Решение. Подставляя a в формулу (5), получим, что поток поля a равен: П = òòò3z2dxdydz , где область V изображена на рис. 3.

V

Рис. 3

до 2R cosθ , получим

Полученный тройной интеграл удобно вычислять в сферических координатах ρ, ϕ,θ ,

для

которых

элемент

объема

dv = ρ2 sinθ d ρ dϕ dθ ,

а уравнения

сферы

и

конуса

имеют вид

ρ = 2R cosθ

 

и

 

θ = π ,

 

 

 

 

 

3

соответственно.

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

области V

угол ϕ изменяется от 0

до 2π , угол θ − от 0 до π

, а ρ от 0

 

 

 

 

3

 

 

 

π 3

2Rcosθ

 

 

П = 3 ò dϕ ò dθ

ò

ρ2 cos2 θ × ρ2 sinθ d ρ =

 

0

 

0

 

0

 

 

 

 

 

π

 

 

 

 

2Rcosθ

 

 

 

3

 

 

 

 

 

 

= ò cos2 θ sinθ dθ ò ρ4d ρ =

 

 

0

 

 

π

 

0

 

 

 

 

 

 

 

 

 

 

 

 

192π R5 3

 

 

153π R5

 

=

 

 

 

 

ò cos7

θ sinθ dθ =

 

.

 

 

5

 

32

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

157

30. Дивергенция векторного поля. Количественной характеристикой векторного поля a является дивергенция (расходимость).

Рассмотрим векторное поле a = Pi + Qj + Rk и выделим в нем

малое тело V объемом v , ограниченное замкнутой поверхностью S . Предел отношения потока вектора a через поверхность S к объему v при стягивании объема v в некоторую точку M V

называется дивергенцией вектора a в точке M или расходимостью и обозначается diva(M ) .

Таким образом, по определению

 

 

 

 

 

òò(a, n)ds

 

div a(M ) = lim

S

 

 

.

(8)

 

v

 

V M

 

 

 

Данное определение дивергенции

 

не связано

с системой

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

Вычислим дивергенцию в декартовой прямоугольной системе

координат

Oxyz ,

считая,

что

координаты

вектора

a = (P(x, y, z);

Q(x, y, z);

R(x, y, z))

непрерывно дифференцируемые

по x, y, z , соответственно, функции в области V.

Используя формулу Остроградского, определение (8) и теорему о среднем значении для тройного интеграла, получим

òòò

æ P

 

Q

 

R ö

 

 

P

 

Q

 

R

 

 

è

x

 

y

 

z ø

 

 

 

 

 

 

diva(M ) = lim

V

ç

 

+

 

+

÷ dxdydz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

+

 

+

 

 

.

 

 

 

 

v

 

 

x

y

z

V M

 

 

 

 

 

 

 

 

 

M

 

 

 

 

 

 

 

 

 

Итак, в декартовой системе координат имеем

 

 

 

 

diva(M ) =

P

+

Q +

R .

 

 

 

 

 

(9)

 

 

 

 

 

x

 

y

z

 

 

 

 

 

 

 

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

Исходя из физического смысла потока (условно считается, что a(M ) есть поле скоростей фиктивного стационарного потока

несжимаемой жидкости), можно сказать, что: при div a(M ) > 0 точка

M представляет собой

источник, откуда жидкость вытекает; при

div a(M ) < 0 точка M

есть сток, поглощающий жидкость. Как

следует из равенства (8), величина div a(M ) характеризует мощность (интенсивность, плотность) источника или стока в точке M . В этом

158

состоит физический смысл дивергенции. Понятно, что, если в объеме V , ограниченном замкнутой поверхностью S , нет ни источников, ни стоков, то div a = 0 .

 

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

 

 

 

 

 

 

 

1) Если c

 

– постоянный вектор, то div c = 0 .

 

 

 

 

 

 

2) div

(c × a) = c ×diva , где c = const .

 

 

 

 

 

 

 

 

 

3) div

(a +

 

 

) = diva + div

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4)

Если

 

U

 

 

скалярная

 

 

функция,

a

 

– вектор, то

div(U ×a) =U ×diva + agradU .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Все перечисленные свойства следуют непосредственно из

формулы (9).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

Так как U × a =U × P ×

 

+U ×Q ×

 

 

+U × R ×

 

, то

 

 

 

 

 

 

i

 

j

k

 

 

 

 

 

div(U × a) =

 

 

(U × P) +

 

(U ×Q) +

 

(U × R) =

 

 

 

 

 

 

 

 

 

 

 

 

y

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= U ×

P

+ P × U

+U × Q

 

+ Q × U

 

+U × R

+ R × U

=

 

 

 

 

 

 

x

 

 

 

 

 

 

x

 

 

 

 

 

 

 

y

 

 

 

 

 

y

 

 

 

 

 

z

 

 

 

z

 

 

 

 

 

 

æ

P

+

Q

+

R ö

+ P ×

U

+ Q ×

 

U

 

+ R ×

U

=U diva

 

+ a gradU.

= U ç

x

y

 

÷

 

x

 

y

 

z

 

è

 

 

 

 

 

z ø

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Упражнение 1. Доказать свойства дивергенции

 

1) − 3).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

a(M ) = 2xy2

 

 

 

+ 3 z2

 

 

в точке

 

M0 (1;-2;1) .

 

 

 

 

 

 

i

- yzj

k

 

 

 

 

 

 

 

 

Решение. Имеем P = 2xy2 , Q = - yz, R = 3z2 . Тогда

 

 

 

 

 

P

 

 

 

 

 

 

= 2y2

 

 

 

Q

 

 

 

 

= -z

 

 

= -1;

R

 

= 6z

 

 

= 6.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 8;

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

M0

 

 

 

 

 

 

 

M0

 

 

 

 

 

y

 

M0

 

 

 

 

M0

 

 

z

 

M0

 

M

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

По формуле (9) получим diva(M0 ) = 8 -1+ 6 =13 .

 

 

 

 

Формула Остроградского (5) в векторной форме с

использованием понятия дивергенции принимает вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

òò(a,

n)ds = òòòdiva

 

dxdydz ,

 

 

 

 

(10)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

 

 

 

 

 

 

 

V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

159

 

 

 

Пример 6.

Найти

 

поток

вектора a = x3

 

 

+ y3

 

+ R2 z

 

 

через

 

 

 

i

j

k

всю поверхность S тела

 

H

(x2

+ y2 ) £ z £ H в направлении внешней

 

R2

нормали.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение. Имеем diva = 3(x2 + y2 ) + R2 . Согласно формуле (10)

 

 

 

 

 

 

òò(a,n)ds = òòò(3(x2 + y2 ) + R2 )dxdydz .

 

 

 

 

 

 

 

 

S

 

 

 

 

V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для

вычисления

 

 

тройного

интеграла перейдем

к

цилиндрическим координатам. Уравнение поверхности S примет вид

z =

Hr2

, и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R2

 

 

 

 

 

 

R

 

 

 

 

 

H

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

òòò(3(x2 + y2 ) + R2 )dxdydz = ò dϕ ò(3r2 + R2 )rdr ò

dz =

 

 

V

 

 

 

 

 

 

 

 

 

 

0

 

0

 

 

 

 

Hr2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R2

 

 

 

 

 

 

 

=

R

(3r2

æ

 

 

2

ö

 

 

 

R

(R4 + 2R2r2 -3r4 )rdr = π HR4 .

ò

+ R2 )ç H - Hr

 

÷ rdr = H

ò

 

 

 

 

ç

R

2

÷

 

R

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

40.

è

 

 

ø

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Соленоидальное

(трубчатое)

поле.

Векторное

поле

a = (P;Q; R) называется соленоидальным

или трубчатым во

всех

точках M некоторой области V, если дивергенция этого поля в этих точках равна нулю, т.е.

diva(M ) = 0, M V .

(11)

Соленоидальное поле не имеет в области V ни источников, ни стоков и обладает следующими свойствами.

1) В соленоидальном поле поток вектора a (M ) через любую

замкнутую поверхность S равен нулю.

Это свойство следует из формулы (10).

2) В соленоидальном поле потоки вектора через различные сечения векторной трубки равны между собой (принцип сохранения интенсивности векторной трубки).

Доказательство. Рассмотрим два сечения S1 и S2 векторной трубки (рис. 4), которые образуют вместе с боковой поверхностью S3 замкнутую поверхность S. Вследствие соленоидальности поля a , получим

160

П = òò(a, n)ds + òò(a, n)ds + òò(a, n)ds =

S1

S2

S3

(12)

= òò(a,

n)ds = òòòdivadxdydz = 0,

S

 

V

 

где V − объем, ограниченный поверхностью S .

 

n3

 

Но

òò(a,n)ds = 0 , так как

 

S

S3

 

 

 

S3

S1

n1

Рис. 4

òò(a,

S1

нормаль к боковой поверхности S3 перпендикулярна вектору a , направленному по касательной к векторной трубке. Из равенств (12) получаем

n)ds = -òò(a, n)ds ,

S2

значит поток вектора a через любое сечение векторной трубки имеет одно и то же значение.

С точки зрения гидродинамики этот факт означает следующее: если a поле скоростей несжимаемой жидкости, то ее объем, протекающий через любое сечение векторной трубки в единицу времени, сохраняется неизменным.

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

Упражнение 2. Установить справедливость свойства 3) соленоидального поля.

Пример 7. Определить, являются ли соленоидальными следующие поля:

1)a1(M ) = xi - 2yj - zk ;

2)a2 (M ) = x(z2 - y2 )i + y(x2 - z2 ) j + z(y2 - x2 )k .

Решение.

1)

Согласно

формуле

(11),

находим

div a1(M ) =1- 2 -1 = -2 ¹ 0. Поле не является соленоидальным.

2) Имеем:

div a (M ) = z2 - y2

+ x2 - z2 + y2 - x2 º 0 ,

то есть

 

 

2

 

 

 

поле a2 (M ) соленоидально.

161

§3. Линейный интеграл, циркуляция и ротор векторного поля

10. Линейный интеграл в векторном поле. Пусть в области V

заданы непрерывное векторное поле a = Pi + Qj + Rk и

ориентированная гладкая линия L. Обозначим через τ единичный вектор касательной к линии L, направление которого совпадает с выбранным направлением на линии.

Линейным интегралом вектора a вдоль линии L называется криволинейный интеграл первого рода от скалярного произведения

векторов a и τ

:

 

 

 

ò(a, τ

) dl ,

(1)

 

L

 

г д е d l – д и ф ф е р е н ц и а л д л и н ы д у г и л и н и и L

( dl = (dx)2 + (dy)2 + (dz)2 ).

Линейный интеграл (1) может быть записан в различных видах.

Если кривая задана вектор-функцией r (t) ,

то вектор dr

будет

направлен

по касательной

к заданной кривой

и dr = τ

dl ,

тогда

ò(a, τ

) dl = ò(a, dr ) .

 

 

 

 

 

 

 

 

 

 

L

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если a = Pi

+ Qj

+ Rk

и r = xi + yj + zk ,

то

 

 

 

 

 

 

ò(a, dr ) = òPdx + Qdy + Rdz .

 

 

(2)

 

 

 

 

 

L

L

 

 

 

Правая часть (2) есть криволинейный интеграл второго рода.

Линейный интеграл меняет знак с плюса на минус, если

изменить направление ориентированной кривой.

 

 

 

Если

линия

 

задана параметрическими уравнениями x = x(t),

y = y(t), z = z(t),

t0 t t1 , то имеем

 

 

 

 

 

t1

 

 

 

 

 

 

 

ò(a, dr ) = ò

 

 

 

 

 

 

 

(P(x(t), y(t), z(t))x (t) + Q(x(t), y(t), z(t)) y (t) +

(3)

 

L

t0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ R(x(t), y(t), z(t))z′(t))dt.

Если a − силовое поле, то линейный интеграл равен работе, которую поле совершает по перемещению материальной точки вдоль ориентированной кривой.

Вычисление линейного интеграла в случае плоского поля проводится по формулам (5.2.8).

162

Пример 1. Вычислить линейный интеграл плоского векторного поля a = (2x y2 +1) i + (3x + 2y2 −10) j по линии L, если L − дуга

параболы x = 3 − y2 , заключенная между точками A(–1;–2) и B(2;1),

(рис.1).

Решение. Согласно формуле (2), имеем

ò(2x y2 +1)dx + (3x + 2y2 −10)dy ,

L

где линия L = ACB (рис. 1) задана уравнением x = 3 − y2 . Подставляя

x = 3 − y2 в выражение (2) и учитывая, что y изменяется вдоль L от –2 до 1, получаем

1

ò ((2(3 −y2 ) − y2 +1)(−2y) + (3(3 − y2 ) + 2y2 −10))dy =

−2

1

= ò (6y3 y2 −14y −1)dy = −7,5 .

 

 

 

 

−2

 

 

 

 

 

 

 

Пример 2.

Найти

работу силы

 

 

 

 

 

 

 

 

 

 

 

F

= xi

+ yj + zk

при

перемещении

материальной точки вдоль первого

витка винтовой линии

 

Рис. 1 x = aet cost, y = aet sin t, z = aet

из точки A(0; 0; 0) в точку B(a; 0; a).

Решение. I способ.

Так как dx = aet (cost − sin t)dt,dy = aet (sin t + cost)dt, dz = aet dt

иa2e2t ((cost − sin t)cost + (sin t + cost)sin t +1)dt = 2a2e2t

учитывая,

что t = −∞ в точке A и t = 0 в точке

 

 

 

0

ò(

 

, dr ) = 2a2 ò e2t dt = a2 .

F

L

−∞

dt,

то,

B,

имеем

 

 

 

 

 

 

II

способ.

В данном

случае

 

 

,dr ) = (r , dr ) =

1

d(r2 ) , где

(F

r =

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r = a

 

 

 

 

 

r

 

,

причем

r = 0 в точке A и

 

в точке B. Имеем:

 

 

2

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

r2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

, dr ) =

1

ò

2

d(r2 ) =

= a2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ò(F

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

2

 

0

 

 

2

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

163

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