2539
.pdf
|
|
|
|
∂U G |
+ |
∂U |
G |
+ |
∂U |
G |
= |
||
div gradU = div |
∂x |
i |
∂y |
j |
∂z |
k |
|||||||
|
|
|
|
|
|
|
|
|
|
||||
= ∂2U + ∂2U |
+ |
∂2U |
= U. |
|
|
|
|
|
|||||
∂x2 |
∂y2 |
|
∂z2 |
div r = ( , rG) и |
rot aG =[ , aG], |
||||||||
в) Используя |
выражение |
|
получим
div rot aG = ( , rot aG)= ( ,[ , aG])=
= |
∂ |
|
∂a |
z − |
∂ay |
+ |
|
|
|
|
|
|
|||
|
|
∂z |
|||||
|
∂x |
∂y |
|
|
∂ |
∂a |
x − |
∂a |
|
+ |
∂ |
∂ay |
− |
∂a |
|
= 0. |
||
|
|
|
|
z |
|
|
|
|
x |
||||
|
|
|
|
∂x |
|
||||||||
∂y |
∂z |
∂x |
|
∂z |
|
∂y |
|
Данная операция приводит к нулю, поскольку мы имеем смешанное произведение компланарных векторов.
г) Действительно
|
|
|
|
|
∂U G |
∂U G |
+ |
∂U G |
|
= rot |
( U )=[ , ]U = |
|||||||||||
rot (grad U )= rot |
∂x |
i + |
∂y |
j |
|
∂z |
k |
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
∂2U |
|
∂2U G |
|
∂2U |
|
∂2U |
|
G |
|
|
∂2U |
|
∂2U G |
|
|||||||
= |
|
− |
|
i + |
|
|
− |
|
|
|
|
|
j |
+ |
|
|
− |
|
k |
= 0. |
||
∂y∂z |
|
∂z∂x |
∂x∂z |
∂x∂y |
|
|||||||||||||||||
|
|
∂z∂y |
|
|
|
|
|
|
|
∂y∂x |
|
Данная операция приводит к нулю, поскольку мы имеем векторное произведение двух коллинеарных векторов.
д) Поскольку rot aG =[ , aG], то rot rot aG = ,[ , aG] .
е) Действительно |
|
∂( |
ax ) |
|
∂( |
ay ) |
|
∂( |
az ) |
|
|
G |
|
|
|
||||||
div ( |
a)= |
|
|
+ |
|
|
+ |
|
|
. |
|
∂x |
|
∂y |
|
∂z |
|||||
|
|
|
|
|
|
|
|
|||
Используя оператор Лапласа |
|
, приведем это выражение к |
виду
131
|
|
|
|
|
G |
|
|
∂ |
|
|
|
|
∂ax |
|
|
|
∂ay |
|
∂az |
|
|
|
∂ |
|
|
|
∂ax |
|
|
∂ay |
|
|
∂az |
|
|
|
|||||||||||
div ( |
|
a)= |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
+ |
|
+ |
|
|
|
+ |
|
|
|
|
|
|
|
+ |
|
|
|
|
+ |
|
|
+ |
|
|||||||
|
|
∂x |
2 |
|
|
∂x |
|
|
∂y |
|
∂z |
|
∂y |
2 |
|
∂x |
∂y |
|
|
∂z |
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
∂ |
|
|
∂ax |
|
|
|
∂ay |
|
|
|
∂az |
|
|
|
|
|
G |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
+ |
|
|
|
|
|
+ |
|
|
|
|
|
|
|
+ |
|
|
|
|
|
|
= div a. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
∂z |
2 |
∂x |
|
|
∂y |
|
|
|
∂z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
grad div aG |
|
||||||||||||||
4.2. Найти: a) |
|
div grad U , если U = x3 y3 z ; б) |
|
, |
|||||||||||||||||||||||||||||||||||||||||||
|
G |
|
|
|
4 G |
|
|
4 |
G |
|
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
G |
|
2 |
yi +z |
2 |
G |
2 |
G |
|
||||||||||
если a |
= x i + y |
|
|
j |
|
+ z |
k ; в) rot rot a , если a = x |
zj +z xk . |
|
||||||||||||||||||||||||||||||||||||||
Решение. а) Найдем сначала градиент функции U = x3 y2 z |
|
||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
grad U = 3x2 y2 zi + 2x3 yzjG+ x3 y3k . |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
Отсюда div grad U = 6xy2 z + 2x3 z = 2xz (3y2 + x2 ). |
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||
б) Найдем сначала дивергенцию вектора a |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
div aG = 4x3 + 4 y3 +4z3 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
Отсюда градиент |
|
|
|
|
|
|
|
|
|
G |
|
|
|
|
|
|
|
|
2 G |
|
|
|
G |
|
|
|
G |
|
|
|
|
|
|
||||||||||||||
|
|
|
|
G |
=12x |
2 G |
+ |
12 y |
2 |
+12z |
2 |
k =12 |
|
+ y |
2 |
+ z |
2 |
|
|
|
|
|
|
||||||||||||||||||||||||
grad div a |
|
|
i |
|
|
j |
|
(x i |
|
j |
k ). |
|
|
|
|||||||||||||||||||||||||||||||||
в) Полагая a =x2 y,a =y2z,a =z2x найдем ротор вектора aG |
|
||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
y |
|
G |
z |
|
|
|
2 G |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
i − z |
|
|
2 |
k . |
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Отсюда rot rot aG = 0 , |
|
rot a = −y |
|
|
j − x |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||
т. е. векторное поле вектора rot aG потен- |
|||||||||||||||||||||||||||||||||||||||||||||||
циально. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
div grad (uv); |
|
|
|
|
|
|
|
|
|
|
G |
|
|
|
|
|
|
|
G |
|
|||||||||||
4.3. |
|
|
Найти: |
|
|
a) |
|
|
|
|
б) |
|
|
|
|
|
|
|
|
|
|
|
|
если |
- |
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
rot rot (cu), |
|
c |
постоянный вектор Решение. а) Найдем сначала градиент произведения двух
скалярных функций
|
∂v |
+v |
∂u G |
|
∂v |
+v |
∂u |
G |
|
∂v |
+v |
∂u |
G |
grad (uv)= u |
∂x |
i |
+ u |
∂y |
|
j |
+ u |
∂z |
k . |
||||
|
|
∂x |
|
|
∂y |
|
|
|
∂z |
|
Дивергенция этого вектора равна
132
div grad (uv)= u |
∂2v + |
|
∂u ∂v |
|
+v ∂2u |
+ |
∂u ∂v |
+u |
∂2v + ∂u ∂v |
+ |
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂x2 |
|
|
|
|
|
∂x ∂x |
|
|
|
|
∂x2 |
|
|
∂x ∂x |
|
|
|
|
|
∂y2 |
∂y ∂y |
|
|
||||||||||||||||||||||||||||
+v |
∂ |
2u |
+ |
|
∂u ∂v |
|
+u |
|
∂ |
2v |
|
+ |
∂u ∂v |
+v |
∂2u |
+ |
∂u |
|
∂v |
= |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||
∂y2 |
|
|
∂y |
∂y |
|
|
∂z2 |
|
∂z ∂z |
∂z2 |
∂z ∂z |
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
2 |
u2 |
|
|
|
|
|
2 |
u2 |
|
|
|
|
|
|
|
|
2 |
u2 |
|
+2 |
|
∂u ∂v |
|
|
|
|
|
∂u ∂v |
|
|
|
∂u ∂v |
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||
= u |
|
∂ |
|
|
+ |
∂ |
|
+ |
|
∂ |
|
|
|
|
|
+ |
|
+ |
|
|
+ |
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
∂x |
|
|
|
|
∂y |
|
|
|
|
|
|
|
|
∂z |
|
|
|
|
|
∂x ∂x |
|
|
|
|
|
∂y ∂y |
|
|
∂z ∂z |
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||
+v |
|
|
|
|
2 |
u2 |
|
|
|
|
2 |
u2 + |
|
|
|
|
|
2 |
u2 |
|
= u div grad v +2(grad u,gradv)+v div gradu. |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
∂ |
|
+ |
∂ |
|
∂ |
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
∂x |
|
|
|
|
|
∂y |
|
|
|
|
|
|
|
|
∂z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
б) Найдем сначала rot (cu) , учитывая, что c - постоянный |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
вектор cG = c (iG |
+ |
|
Gj + k ) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
G |
|
|
= |
|
|
∂u |
|
− |
∂u |
|
G |
+ |
∂u |
− |
∂u |
|
|
G |
+ |
|
∂u |
− |
|
∂u |
G |
|
|
|
|||||||||||||||||||||||||||||||||||||||
|
rot (cu) |
|
|
∂y |
|
|
∂z |
ci |
|
∂z |
|
∂x |
cj |
|
∂x |
|
∂y |
ck . |
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||
Отсюда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
G |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
G |
|
|
|
|
|
|
|
|
|
∂2u |
|
|
|
|
∂2u |
|
∂2u |
|
|
|
|
|
∂2u |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
rot rot (cu) |
= |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
− |
|
|
2 − |
|
|
|
2 |
+ |
|
|
|
|
|
|
|
|
|
ci |
|
+ |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
∂x∂y |
|
|
∂y |
∂z |
|
∂x∂z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
∂2u |
|
|
|
|
|
∂2u |
|
|
|
|
|
|
∂2u |
|
|
|
|
|
∂2u |
|
|
G |
|
|
|
|
|
∂2u |
|
|
|
|
∂2u |
|
|
|
∂2u |
∂2u |
|
G |
|
||||||||||||||||||||||||||||||||||
+ |
|
|
|
|
|
|
|
|
− |
|
|
|
2 |
|
|
− |
|
|
|
|
|
|
|
2 |
|
+ |
|
|
|
|
cj |
|
+ |
|
|
|
|
|
|
|
|
|
− |
|
|
|
|
2 |
|
− |
|
|
|
2 + |
|
|
ck |
= |
||||||||||||||||||
∂y∂z |
∂z |
|
|
|
∂x |
|
∂x∂y |
|
∂x∂z |
∂x |
|
∂y |
∂y∂z |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||
|
∂2u |
|
|
|
|
|
∂2u |
|
|
|
|
|
|
|
|
∂2u |
|
G |
|
|
∂2u |
|
|
|
|
|
∂2u |
|
|
|
∂2u |
|
|
G |
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||
= |
|
|
|
|
2 |
+ |
|
|
|
|
|
|
|
|
|
+ |
|
|
|
|
|
|
|
|
|
|
|
|
ci |
|
+ |
|
|
|
|
|
2 |
+ |
|
|
|
|
|
|
|
|
|
+ |
|
|
|
|
|
|
|
|
cj |
+ |
|
|
|
|
|
|||||||||||||
∂x |
|
∂x∂y |
|
∂x∂z |
|
∂y |
|
|
∂x∂y |
|
∂y∂z |
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||
|
∂ |
2 |
u2 + |
|
|
|
2 |
u |
|
|
|
|
|
|
|
|
|
|
|
2 |
u |
|
|
|
|
G |
|
|
∂ |
2 |
u2 + |
|
|
|
|
2 |
u2 + |
∂ |
2 |
u2 |
|
|
|
|
G |
|
|
G |
G |
|
|
|
||||||||||||||||||||||||
+ |
|
|
|
|
∂ |
|
|
|
+ |
|
|
|
∂ |
|
|
|
|
ck |
|
− |
|
|
|
|
∂ |
|
|
(ci |
+cj |
+ck )= |
|
|||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
∂x∂z |
|
|
|
∂y∂z |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
∂z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂x |
|
|
|
|
|
∂y |
|
|
|
|
∂z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
|
|
|
∂ |
|
|
|
|
|
|
∂ |
|
|
|
|
|
|
|
|
|
|
|
∂ |
|
∂u |
G |
|
|
∂u |
G |
|
|
|
|
∂u |
|
G |
|
|
|
|
|
|
|
|
|
G |
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
= c |
|
|
|
|
|
+c |
|
|
|
|
|
|
+c |
|
|
|
|
|
|
|
|
|
|
i |
+ |
|
|
|
|
|
j |
+ |
|
|
|
|
k |
|
− uc |
= |
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
∂x |
∂y |
|
|
|
∂z |
∂x |
∂y |
|
∂z |
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
= (cG, ) grad u − |
|
|
|
|
ucG. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
133
2.5. Интегралы теории поля и теории потенциала
1°. Поток вектора. Потоком векторного поля a(M ) через
поверхность S называется поверхностный интеграл 2-го рода, который может быть сведен к вычислению поверхностного ин-
теграла 1-го рода
Q = ∫∫(aG, nG)dS = ∫∫andS = ∫∫axdydz + ay dxdz + az dxdy =
|
S |
|
|
|
S |
|
S |
|
(1) |
|
= |
∫∫( |
a |
x |
cosα +a |
y |
cos β + a |
z |
cosγ |
) |
|
|
|
|
|
dS, |
|
S
где nG{cosα, cos β, cosγ} - единичный вектор нормали к повер-
хности S.
Вычисление интеграла (1) может быть представлено в виде
суммы трех двойных интегралов |
|
|
|
|
|
|
|
||||||||||||
∫∫(aG, nG)dS =∫∫ax (x (y, z), y, z)dydz + |
|
|
|
|
|||||||||||||||
S |
|
|
y ( |
|
|
S1 |
|
|
|
) |
|
|
|
z ( |
|
|
|
)) |
(2) |
+ |
∫∫ |
a |
x, y |
( |
x, z |
) |
, z |
dxdz + |
∫∫ |
a |
x, y, z |
( |
x, y |
||||||
|
|
|
|
|
|
|
|
dxdy. |
|
||||||||||
|
S2 |
|
|
|
|
|
|
|
|
|
|
S3 |
|
|
|
|
|
|
|
где S1, S2 , S3 |
- проекция поверхности S на плоскости |
||||||||||||||||||
Oyz,Oxz,Oxy ; |
переменные |
x = x (y, z), y = y (x, z), z = z (x, y) |
|||||||||||||||||
находятся из уравнения поверхности S0 разрешая его |
|||||||||||||||||||
относительно соответствующих координат. |
|
|
|||||||||||||||||
Если вектор aG |
определяет поле скоростей текущей жидко- |
сти или газа, то поверхностный интеграл определяет количество жидкости или газа, протекающей за единицу времени через поверхность S в заданном направлении. Переход к другой стороне поверхности S меняет направление нормали на противоположное, а следовательно, и знак поверхностного интеграла.
При явном задании поверхности z = f (x, y) |
элемент пло- |
щади имеет вид |
|
dS = 1+(z′x )2 +(z′y )2 dxdy , |
(3) |
134
а направляющие косинусы нормали к поверхности находятся по формулам
cosα = |
|
−z′x |
, cos β = |
|
−z′y |
|
, |
||
± |
1+(z′x )2 +(z′y )2 |
± |
1+(z′x )2 |
+(z′y )2 |
|||||
cosγ = |
|
|
1x |
|
|
|
|
|
|
|
|
. |
|
|
(4) |
||||
± |
1+(z′x )2 +(z′y )2 |
|
|
Подставляя в формулу (1) выражения (3), (4), находим еще одну формулу для вычисления потока вектора a через поверхность
Q = ∫∫ ax (−z′x )+ ay (−z′y )+ az dxdy |
(5) |
S |
|
Здесь предполагается, что нормаль образует острый угол с |
|
осью Oz , следовательно, cosγ > 0 и в формулах (4) |
берется |
знак «+». В случае cosγ < 0 в формулах (4) следует брать знак
«-».
2°. Формула Остроградского-Гаусса. Если V - некоторая замкнутая область пространства, ограниченная гладкой поверхностью S , и функции P(x,y,z), Q(x,y,z), R(x,y,z) непрерывны со своими частными производными первого порядка в области V, включая границу S , то справедлива формула Остроградского - Гаусса
|
|
|
∂P |
|
∂Q |
|
∂R |
|
∫∫Pdydz +Qdzdx + Rdxdy = ∫∫∫ |
+ |
+ |
dxdydz . (6) |
|||||
S |
V |
|
∂x |
|
∂y |
|
∂z |
|
Если поверхностный интеграл второго типа заменить поверхностным интегралом первого типа, то формула (2) примет вид
|
|
|
∂P |
|
∂Q |
|
∂R |
|
∫∫(P cosα +Q cos β + R cosγ )dS = ∫∫∫ |
+ |
+ |
dxdydz . (7) |
|||||
S |
V |
|
∂x |
|
∂y |
|
∂z |
|
где α, β,γ - углы между внешней нормалью к поверхности S и координатными осями.
135
В векторной форме формула Остроградского - Гаусса име-
ет вид |
∫∫(aG, nG)dS = ∫∫∫div aG dxdydz , |
|
|
|
|
|||
|
|
|
|
|
(8) |
|||
G |
G |
Gs |
G |
V |
G |
|
|
G |
где a |
= axi |
+ ay j + az k |
или a |
=P(x, y, z)i +Q(x, y, z) j |
+R(x, y, z)k . |
|||
Формула (8) определяет поток векторного поля |
G |
через |
||||||
a |
замкнутую поверхность S в направлении ее внешней нормали
п.
3°. Формула Грина. Пусть L - граница некоторой области S в плоскости Оху и функции P (x, y) и Q (x, y) непрерывны,
вместе со своими частными производными ∂∂Qx и ∂∂Py в
области S, включая границу L. В этом случае справедлива формула Грина.
v∫ |
|
∂Q |
|
∂P |
|
(9) |
Pdx +Qdy = ∫∫ |
− |
dxdy , |
||||
L |
S |
∂x |
|
∂y |
|
|
которая преобразует криволинейный интеграл, взятый по замкнутому контуру, в двойной интеграл по области S, ограниченный этим контуром. Обход контура L выбирается так, чтобы область S оставалась слева, т. е. против часовой стрелки.
4°. Формула Стокса. Если L — некоторый замкнутый контур поверхности S и функции P(x, y, z) , Q(x, y, z) и R(x, y, z)
прерывны вместе со своими частными производными первого порядка на поверхности S, включая границу L, то справедлива формула Стокса
∫Pdx +Qdy + Rdz = |
|
|
|
|
|
|
||||
L |
|
|
|
|
|
|
|
|
|
(10) |
|
|
∂Q |
|
∂P dxdy + |
∂R |
|
∂R dydz + |
|
|
|
= |
∫∫S |
− |
− |
∂P |
− |
∂R dxdz. |
||||
|
∂x |
|
∂y |
∂y |
|
∂x |
∂z |
|
∂x |
|
|
|
|
|
|
|
|
|
|
|
|
Следует заметить, что сторона поверхности и направление обхода контура определяются так же, как и в случае формулы Грина. Если в качестве поверхности S взять плоскую область
136
на плоскости Оху, так что z = 0, то формула (10) преобразуется в формулу Грина.
Если поверхностный интеграл второго типа заменить поверхностным интегралом первого типа, то формула Стокса примет вид
|
|
|
∫Pdx +Qdy + Rdz = |
|
|
|
(11) |
|||
|
|
|
L |
|
|
|
|
|
|
|
= |
∂Q |
− |
∂P cosα + |
∂R |
− |
∂Q cos β + |
∂P |
− |
∂R cosγ ds, |
|
|
∫∫S ∂x |
|
∂y |
∂y |
|
∂z |
∂z |
|
∂x |
|
|
|
|
|
|
|
|
|
|
|
|
|
G |
, |
cos β |
|
G |
, |
где cosα = cos xn |
= cos yn |
|||||
|
|
|
|
|
|
|
|
G |
— |
cosγ = cos zn |
||
|
|
|
направляющие косинусы внешней нормали n к поверхности S. |
||
5°. Работа поля. Циркуляция вектора. Работа поля aG |
вдоль |
|
кривой L определяется интегралом. |
|
|
∫(aG, drG)= ∫axdx + ay dy + az dz . |
(12) |
|
L |
L |
|
Если кривая L - замкнутая, то интеграл (12) называется
циркуляцией векторного поля a вдоль кривой L, обозначается
C = v∫(aG, drG)
L
и характеризует вращательную способность поля на контуре L. Если замкнутая кривая L ограничивает поверхность S, то имеет место формула Стокса, которая в векторной форме при-
мет вид |
|
C = v∫(aG, drG)= ∫∫(rot aG, nG)dS |
(13) |
L
и означает, что циркуляция векторного поля a по замкнутому контуру L равна потоку вектора rot а через поверхность S, ограниченную этим контуром L.
6°. Потенциал. Скалярная величина U называется потенциалом векторного поля a , если grad U = a . Само же
векторное поле aG в этом случае называется потенциальным. Необходимым и достаточным условием потенциальности поля
137
aG является равенство нулю вихря этого поля rot aG = 0 . Поскольку
dU = |
∂U dx + |
∂U dy + |
∂U dz = axdx +ay dy + az dz, |
|
∂x |
∂y |
∂z |
то потенциальная функция U с точностью до произвольного постоянного слагаемого определяется линейным интегралом
U = ∫ (aG, drG)= ∫ axdx + ay dy + az dz = ∫B dU =U (B) −U ( A) . |
(14) |
||
AB |
AB |
A |
AB и |
Линейный интеграл не зависит от формы кривой |
определяется разностью значений потенциального поля в конце и начале пути интегрирования. Поэтому за путь проще всего взять ломаную со звеньями, параллельными координатным
осям, а за начальную точку A - начало координат. |
|
||||||
|
G |
2 |
i − y |
2 G |
2 |
k |
через часть |
5.1. Найти поток вектора a |
= x |
j |
+ z |
||||
сферы x2 + y2 + z2 = R2 , |
x ≥ 0, y ≥ 0, z ≥ 0 , |
в |
|
направлении |
|||
внешней нормали. |
|
|
|
|
|
|
|
Решение. Используя формулу (1), имеем |
|
|
|
||||
∫∫(aG, nG)dS = ∫∫x2dydz − y2dxdz + z2dxdy . |
|||||||
S |
S |
|
|
|
|
|
|
Таким образом, задача сводится к вычислению трех интегралов. Преобразуя эти поверхностные интегралы в двойные по формуле (2), будем иметь
∫∫(aG, nG)dS = ∫∫x2dydz −∫∫y2dxdz + ∫∫z2dxdy =
S |
S1 |
S2 |
S3 |
= ∫∫(R2 − y2 − z2 )dydz −∫∫(R2 − x2 − z2 )dxdz +
S1 |
S2 |
+∫∫(R2 − x2 − y2 )dxdy.
S2
Переходя к полярным координатам, получим
138
|
π |
π |
|
|
∫∫(aG, nG)dS = ∫2 dϕ∫R (R2 − ρ2 )ρd ρ − ∫2 dϕ∫R (R2 − ρ2 )ρd ρ + |
||||
S |
0 0 |
0 |
0 |
|
π |
|
|
|
|
+∫2 dϕ∫R (R2 − ρ2 )ρd ρ = |
πR4 . |
|
|
|
0 |
0 |
8 |
G |
3 G |
|
|
G |
||
5.2. Найти поток векторного поля a |
= xyi + zj |
+ 2 y k через |
||
часть поверхности 4x2 + y2 |
= 4 − z , лежащую в первом октанте |
x ≥ 0, y ≥ 0, z ≥ 0 , в направлении нормали, образующей острый
угол с осью Oz .
Решение. Для нахождения потока воспользуемся формулой (1). Разрешим уравнение поверхности относительно z и найдем производные
z = 4 −4x2 − y2 , z′x = −8x , z′y = −2 y .
Элемент площади поверхности и направляющие косинусы нормали к поверхности S находятся по формулам (3), (4). Поверхность S представляет эллиптический параболоид (рис. 2.5).
Рис. 2.5
Спроектируем его на плоскость Оху и рассмотрим четвертую часть D. По условию нормаль образует острый угол с осью Oz , следовательно, cosγ > 0 и в формулах (4) перед
корнями следует брать знак (+). Расчетная формула примет вид (5). Координаты вектора a по условию равны ax = xy ,
ay = 4 −4x2 − y2 , az = 2 y3 , отсюда
139
Q = ∫∫(8x2 y +(4 −4x2 − y2 )2 y + 2 y3 )dxdy =
D
=8∫1 |
dx |
2 |
1∫−x2 |
ydy = 4∫1 |
y2 |
|
02 |
1−x2 |
dx =16∫1 (1− x2 )dx = 32 . |
|||
|
|
|||||||||||
0 |
|
|
0 |
0 |
|
|
|
|
|
0 |
|
3 |
5.3. Найти поток вектора |
G |
G |
1 G |
2 |
G |
|||||||
a = zi |
+ 6 xj |
+ 9 yk через часть |
||||||||||
поверхности x +2 y +3z =1, |
лежащую в первом октанте x ≥ 0 , |
y ≥ 0 , z ≥ 0 , в направлении нормали, образующей острый угол
с осью Oz .
Решение. Воспользуемся формулой (5). Для этого найдем
проекции |
вектора |
на оси координат: ax |
= z = 13 (1− x −2 y), |
|||||||||||
ay = 1 x , |
az = 2 y и производные z′x = −1 , z′y |
= − 2 |
: |
|||||||||||
6 |
9 |
|
|
|
|
|
|
3 |
|
|
3 |
|
||
|
Q = |
∫∫ |
1 |
(1− x −2 y) |
1 |
+ |
1 |
x |
2 |
+ |
2 |
|
|
|
|
3 |
3 |
6 |
3 |
9 |
|
|
|
||||||
|
|
|
|
|
|
|
y dxdy . |
|||||||
|
|
Sxy |
|
|
|
|
|
|
|
|
|
|
|
|
Проекция поверхности на плоскость Ox y |
представляет об- |
|||||||||||||
ласть Sx y , ограниченную прямыми |
x = 0, y = 0 |
и x +2 y =1 |
||||||||||||
(рис. 2.6). Отсюда |
|
|
|
|
|
|
|
|
|
|
|
|
|
∫1 |
y= |
1 |
(1−x) |
Q = |
1 |
dx |
2∫ |
||
|
9 |
0 |
|
|
0 |
|
|
|
1 |
(1−x) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
1 |
2 |
|
1 1 |
|
|
|
|
|
x2 |
|
1 |
1 |
|
|||
|
|
|
(1 |
− x)dx = |
1 |
|
|
|||||||||||
dy = |
|
y |
|
dx = |
|
|
|
|
|
x − |
|
|
|
= |
|
. |
||
|
|
|
|
|
|
|
|
|
||||||||||
|
9 |
∫0 |
|
|
18 |
∫0 |
|
|
18 |
|
2 |
|
|
0 |
36 |
|
||
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 2.6
140