Tom_2
.pdf2 |
|
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 |
|
|
|
|
|
|
|
|
|
|
2π |
|
|
|
|
|
|
|
|
C |
|
|
|
|
|
|
|
|
|
|
|
2π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 |
|
|
2π |
|
π 3 |
2Rcosθ |
|
|
|||
П = 3 ò dϕ ò dθ |
ò |
ρ2 cos2 θ × ρ2 sinθ d ρ = |
|||||||
|
0 |
|
0 |
|
0 |
|
|
|
|
|
|
π |
|
|
|
|
2Rcosθ |
|
|
|
3 |
|
|
|
|
|
|
||
= 6π ò 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 |
|
|
|
|
|
2π |
|
R |
|
|
|
|
|
H |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
òòò(3(x2 + y2 ) + R2 )dxdydz = ò dϕ ò(3r2 + R2 )rdr ò |
dz = |
|
|
|||||||||||||||||||||||
V |
|
|
|
|
|
|
|
|
|
|
0 |
|
0 |
|
|
|
|
Hr2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
R2 |
|
|
|
|
|
|
|
|
= 2π |
R |
(3r2 |
æ |
|
|
2 |
ö |
|
|
|
R |
(R4 + 2R2r2 -3r4 )rdr = π HR4 . |
□ |
|||||||||||||
ò |
+ R2 )ç H - Hr |
|
÷ rdr = 2π 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