Скачиваний:
44
Добавлен:
10.05.2014
Размер:
1.75 Mб
Скачать

Пример 5 (4365).

 

 

 

 

 

 

 

 

2

 

 

2

 

 

2

 

Вычислить

 

dydz

+

dzdx

+

dxdy

, Ф-внешняя сторона эллипсоида

x

 

+

y

 

+

z

 

= 1 .

 

 

 

 

2

 

2

 

2

 

∫∫

 

 

 

 

 

 

 

 

 

 

 

 

Φ x

 

y

 

z

 

a

 

 

b

 

 

c

 

 

z = c 1 −

x2

y2

. p = −

 

c

 

 

 

x

 

 

 

 

 

= −

c2 x

 

, q = −

 

c

 

 

 

y

 

 

 

= −

c2 y

.

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

a

2

 

b

2

 

 

 

a

 

 

 

 

x

2

 

 

 

y

2

 

 

 

a

2

z

 

 

 

 

b

 

 

 

x

2

 

y

2

 

b

2

z

 

 

 

 

 

 

 

 

 

 

 

 

 

1 −

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 −

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a2

b2

 

 

 

 

 

 

 

 

 

 

 

 

 

a2

b2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c2

 

c2

 

 

1

 

 

 

c2

1

 

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N = (− p,−q,1), (V , N ) =

 

 

 

 

+

 

 

+

 

 

=

 

 

 

 

 

+

 

 

+

 

 

 

. Обозначим через - верхний

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a2 z b2 z z

 

 

 

z

a2

 

 

 

b2

 

c2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

полуэллипсоид, а через D его проекцию на плоскость xOy. Учитывая ранее сделанное замечание и симметрию относительно координатных осей, получим

 

 

 

 

 

 

 

dzdx

 

 

 

dxdy

 

 

 

 

 

 

 

dzdx

 

 

 

dxdy

 

 

 

 

 

 

 

 

 

 

 

 

dydz

 

+

+

=2

 

 

dydz

+

+

= 2

∫∫

(V , N )dxdy =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∫∫

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∫∫

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Φ x

 

 

 

 

y

 

 

 

 

 

 

z

 

 

 

 

 

x

 

y

 

 

 

 

 

 

 

 

z

 

 

 

D

 

 

 

 

 

1

 

 

 

1

 

 

 

1

 

 

 

 

 

 

dxdy

 

 

 

 

 

 

1

 

 

 

 

 

1

 

 

 

 

 

1 2π

1

 

rdr

2c

 

 

+

 

 

 

 

+

 

 

 

∫∫

 

 

 

 

 

 

 

 

 

 

 

 

 

= 2abc

 

 

 

 

+

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

=

 

2

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

2

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

b

 

 

 

 

c

 

 

 

D

 

1 −

y2

 

 

a

 

 

 

 

 

b

 

 

 

 

 

 

c

 

 

 

0

0

1 − r 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a2

b2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

 

1 1

1

 

 

du

 

 

 

 

 

 

1

 

 

 

 

1

 

 

 

 

 

1

 

 

 

 

 

 

4πabc

 

 

+

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

= 4πabc

 

 

+

 

 

 

 

 

 

+

 

 

 

 

.

 

 

 

 

 

 

 

 

2

b

2

c

2

 

 

 

 

 

 

 

 

 

 

2

b

2

 

 

c

2

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

2

0

 

u

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

§3. Формула Стокса

1. Поверхность, заданная уравнением z = z(x, y)

Рассмотрим ориентированную непрерывно дифференцированную поверхность Ф, однозначно проектируемую не все координатные плоскости. Пусть эта поверхность имеет задание z = z(x, y), (x, y) Dz относительно переменных x,y. Через Г обозначим край этой поверхности с согласованной ориентацией. Согласованной ориентацией называется такая, что при обходе края поверхности в этом направлении с выбранной нормалью, поверхность остается слева.

Пусть P(x,y,z) задана и непрерывна на Ф и имеет там непрерывные частные производные

P

 

P

y

,

z .

Тогда имеет место равенство

 

 

 

 

 

 

 

 

 

 

P dzdx

 

 

=

 

P( x, y, z)dx .

 

 

 

 

P dxdy

 

 

 

∫∫

z

y

 

 

 

 

 

 

Φ

 

 

 

 

 

 

 

Доказательство проведем для положительно ориентированной поверхности

61

Обозначим P*(x,y)=P(x,y,z(x,y)). Отметим, что P( x, y, z) dx = P* ( x, y) dx . Это следует из формулы

γ

вычисления криволинейного интеграла. Действительно, рассмотрим параметризацию кривой γ.

 

x(t)

 

x(t)

 

γ :

:

y(t)

 

t [α, β] тогда

t [α, β].

 

y(t)

 

 

 

 

 

 

z( x(t), y(t))

β

Тогда P* ( x, y) dx = P( x, y, z( x, y)) dx = P( x(t), y(t), z( x(t), y(t))) x' (t)dt ,

 

 

γ

γ

 

 

α

 

 

β

 

 

 

 

P( x, y, z) dx = P( x(t), y(t), z( x(t), y(t))) x' (t)dt .

 

 

α

 

 

 

 

По формуле Грина

 

 

 

 

P( x, y, z) dx = P* ( x, y) dx = ∫∫

P* ( x, y)dxdy =

 

 

 

γ

 

D

y

 

 

 

 

 

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

 

 

 

 

P( x, y, z( x, y)) ∂z dxdy =

 

∫∫

y

 

z

 

 

 

D

 

 

y

∫∫P(x, y, z( x, y))dxdy ∫∫

P( x, y, z( x, y))

z

dxdy = ∫∫

P( x, y, z)dxdy

y

D

y

D

z

Φ

y

 

 

 

 

= ∫∫

P( x, y, z)dxdy ∫∫P( x, y, z) q cosγ dS =

 

 

Φ

y

Φ

z

 

 

 

 

 

 

 

 

∫∫

P( x, y, z)

z

dxdy

 

Φ

z

y

 

 

 

∫∫

P( x, y, z)dxdy + ∫∫

P( x, y, z) cos β

Φ

y

Φ

z

 

 

dS =

 

P( x, y, z) dzdx

 

 

 

P(x, y, z) dxdy . Здесь

 

∫∫

z

y

 

 

Φ

 

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

cos β =

 

 

q

 

, cos γ =

1

, откуда q cos γ = - cos β .

 

 

 

 

 

 

 

 

 

 

 

 

1

+ p2

+ q2

1 + p2 + q2

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

2. Формула Стокса для векторного поля.

62

Пусть Ф допустимая, ориентированная поверхность и V=(P,Q,R) непрерывное на Ф поле, Г-край этой поверхности с согласованной ориентацией.

 

 

P dzdx

 

 

 

 

Из доказанной формулы

∫∫

 

 

=

P( x, y, z)dx формальной заменой z на y,

 

P dxdy

 

 

Φ

z

y

 

 

 

 

 

 

R dydz

 

 

 

 

x на z, P на R (см. рисунок, заменяем 1 на 2) получим

∫∫

 

 

=

R( x, y, z)dz

 

R dzdx

 

 

Φ

y

x

 

 

 

 

Точно так же заменой z на x, y на z, x на y, P на Q (см. рисунок, заменяем 1 на 3) получим

y на x,

.

 

Q

dxdy

Q

 

= Q( x, y, z)dy .

Складывая полученные выражения, получим

∫∫

x

z

dydz

Φ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R Q

 

P R

Q P

 

 

 

Pdx

+ Qdy + Rdz =

∫∫

 

 

 

 

 

 

 

 

 

 

 

 

dydz +

dzdx +

dxdy .

 

 

 

 

 

 

 

 

Φ

y

 

z

 

z

x

x

 

y

Векторная запись формулы Стокса. Ротор определяется по формуле

 

 

r

 

r

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

j

 

 

k

 

 

 

 

 

 

 

 

 

 

rot V =

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

y

 

 

z

 

 

 

 

 

 

 

P

 

Q

 

 

R

 

 

 

 

 

 

 

 

 

Тогда формула Стокса запишется в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cosα cos β

cosγ

 

 

 

 

 

 

 

 

(V, ds) = ∫∫ (rot V, dS)= ∫∫

 

 

 

 

 

 

 

 

 

 

 

dS .

 

x

 

 

 

 

 

y

z

 

Φ

 

Φ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P

 

 

 

 

 

Q

 

R

 

 

Циркуляция векторного поля по краю поверхности равна потоку ротра через эту поверхность. Подробнее о смысле этой терминологии будет сказано позже.

Пример 1.(4367) Вычислить ydx + zdy + xdz , С- окружность x2+y2+z2=a2, x+y+z=0 , проходимая

C

против часовой стрелки, если смотреть с положительной стороны оси Ox.

 

r

 

r

 

r

 

 

 

 

 

 

i

 

j

 

k

 

rot V =

 

 

=(-1,-1,-1). В качестве поверхности с краем C выберем круг сечения плоскости

x

 

y

 

z

 

y

 

z

 

x

 

 

 

 

 

 

 

 

x+y+z=0 шара x2+y2+z2a2 , ориентированный нормалью (1,1,1). Тогда

63

ydx + zdy + xdz = ∫∫

(rot V, dS)= ∫∫ (rot V, n) dS=

3

∫∫

dS=

3

µФ=

3

π a2.

C

 

 

 

Φ

 

 

 

 

 

 

 

 

Φ

 

 

 

 

 

 

 

Φ

 

 

 

 

 

Пример 2.(4368) Вычислить ( x2 yz)dx + ( y2 xz)dy + ( z2 xy)dz , взятый по отрезку винтовой

 

 

 

 

 

AmB

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

линии x=a cos ϕ , y=a sin ϕ , z=

 

h

 

ϕ от А(а,0,0) до B(a,0,h).

 

 

 

 

 

 

2π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

r

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

j

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rot V =

 

 

 

 

 

=(0,0,0),

поэтому интеграл не зависит от пути интегрирования и

 

 

 

 

 

 

 

 

 

 

 

 

x

 

y

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2 yz y2 xz z 2 xy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x = a

вместо винтовой линии выберем отрезок, соединяющий точки A, B. y = 0 t [0, h] .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z = t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h

h

3

 

 

 

 

 

 

 

 

 

 

( x2 yz)dx + ( y2 xz)dy + ( z2 xy)dz = t 2 dt =

 

.

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

AmB

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 3.(4369) Доказать формулу

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

dx

dy

 

 

dz

 

1

([n, r], ds),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

µΦ =

 

cosα

cos β

cos γ

=

 

 

 

 

 

 

 

 

 

 

 

 

 

2

∂Φ

 

x

y

 

 

2

∂Φ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

где Ф- область, лежащая в плоскости с единичной нормалью n , ограниченная кривой ∂Φ ,

согласованно ориентированной с нормалью n .

V = [n, r] = ( z cos β y cosγ , x cosγ z cosα , y cosα x cos β ) ,

 

 

 

r

 

 

r

 

r

 

 

 

 

 

 

 

 

 

 

 

i

 

 

j

 

k

 

rot V =

 

 

 

 

 

 

=(2cosα,2cosβ,2cosγ)=2 n .

 

x

 

 

y

 

z

 

 

 

 

 

 

 

 

 

z cos β y cosγ x cosγ z cosα

y cosα x cos β

 

 

dx

dy

dz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cosα

cos β

cosγ

= (V , ds)= ∫∫(rot V ,dS )= ∫∫(2n, n)dS = 2∫∫dS = 2µΦ .

∂Φ

x

y

z

∂Φ

 

Φ

Φ

Φ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 4. Вычислить ( y z)dx + ( z x)dy + ( x y)dz . С- контур x=a sin2t, y=a sin t cos t, z=a cos2t

C

, t [0,π].

 

x

 

y

 

a

2

a2

Контур лежит в плоскости x+z=a , далее

 

= tg t =

 

, y2=x z , y2=x (a – x) , или x

 

 

+ y2 =

 

.

 

 

 

 

 

y

 

z

 

2

 

4

 

Таким образам, этот контур является эллипсом с полуосями a , a .

2 2

64

 

 

r

 

r

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

j

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= (−2,−2,−2) , n =

1

 

(1,0,1)

, (rot V , n)= −2

 

 

 

 

 

 

rot V =

 

 

 

 

 

 

2 ,

 

 

 

x

 

y

 

z

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y z

z x

x y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( y z)dx + ( z x)dy + ( x y)dz = ∫∫(rot V ,dS )= −2

 

∫∫dS = −2

 

π

a

 

 

a

 

= −πa2

2

2

 

 

 

 

 

 

 

 

C

 

 

 

 

 

 

 

 

 

 

Φ

 

 

 

 

 

Φ

2

 

2

 

 

3. Условия независимости криволинейного интеграла от пути интегрирования в пространстве

Лемма. Для того, чтобы интеграл

(V, ds)

(1)

 

 

не зависел от пути интегрирования Г (соединяющего две фиксированные точки А , В ) необходимо и достаточно, чтобы интеграл (1) был равен нулю для любого замкнутого контура, лежащего в области (циркуляция поля равна нулю по любому кортуру из области).

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

Определение. Область D называется поверхностно односвязной, если для любой кусочно гладкой замкнутой кривой Г (контура Г), лежащей в D можно указать ориентированную допустимую поверхность Ф, расположенную в D, краем которой является Г.

Примеры. Шар является поверхностно односвязной областью. Тор не является поверхностно односвязной областью.

Теорема 1. Для того, чтобы интеграл (1) не зависел от пути в поверхностно односвязной области D, необходимо и достаточно, чтобы rot V =0 в области D.

Доказательство. Достаточность (дано: rot V =0 ) следует из формулы Стокса и сформулированной леммы.

Необходимость. Предположим противное, существует точка М0 (её можно считать внутренней точкой области D), где rot V 0. Следовательно, одна из координат этого вектора в этой точке будет

отлична от нуля. Пусть, например, R(M 0 ) Q(M 0 ) = h > 0. Найдется окрестность этой точки, в

y z

65

которой будет выполняться условие

R(M )

Q(M ) >

h

и которая будет лежать в D.

 

 

y

z

2

 

Сечение этой окрестности плоскостью перпендикулярной оси x и проходящей через точку М0 обозначим K (круг радиуса ε, ориентированный ортом оси x) , а его границу – через C (окружность с согласованной ориентацией). D-проекция K на плоскость yOz

Используя формулу Стокса, получим противоречие:

0 =

(V , ds) =

 

R

Q

 

 

R(

 

dydz =

 

 

 

∫∫

y

 

 

∫∫

 

 

C

 

K

z

D

 

x0 , y, z)

Q( x0 , y, z)

 

h

πε 2 .

dydz >

 

y

z

 

2

 

 

 

Теорема 2. Для того, чтобы интеграл (1) не зависел от пути интегрирования в поверхностно односвязной области D , необходимо и достаточно, чтобы подинтегральное выражение было полным дифференциалом

Pdx+Qdy+Rdz = du.

 

 

 

 

Доказательство. Достаточность. Если Pdx+Qdy+Rdz = du , то P =

u , Q =

u ,

R =

u . Откуда

 

x

y

 

z

следует, что rot V =0 .

 

 

 

 

Необходимость. Определим функцию u по формуле

 

 

 

 

M ( x, y , z )

 

 

 

 

u(x,y,z) = (V, ds) ,

 

 

 

 

M 0

 

 

 

 

где М0(x0, y0, z0) – фиксированная точка в области D , а интегрирование ведется по некоторому пути, соединяющему точки М0 и М. Так как интеграл не зависит от пути интегрирования, то определение

корректно. Докажем, что таким образом определенная функция является искомой, т. е. P = u ,

x

Q =

u , R =

u . Вычислим производную

u

непосредственно по определению.

 

y

z

x

 

Для отрезка ММиспользуем параметризацию

66

x +

xt

y

,t [0,1] . Тогда

u( x + x, y, z) − u( x, y, z)

 

1

M '

1

1

=

P( x, y, z)dx =

P( x + xt, y, z) xdt =P(x+θΔx,y,z), откуда и

 

x

x

x

M

0

 

 

 

 

u

следует требуемое соотношение для частной производной x . Аналогично проводится доказательства для других производных.

§4. Формула Остроградского Гаусса

Определение. Объемно односвязной областью называется область D, удовлетворяющая следующему свойству. Любая замкнутая, кусочно-гладкая, не самопересекающаяся поверхность, расположенная в D , является границей области целиком лежащей в D. Можно сказать, что внутри области нет полостей.

Рассмотрим объемно односвязную область W и функцию R , определенную в этой

R

области и имеющую там непрерывную производную z . Границу этой области, ориентированную положительно, обозначим ∂ W .

Тогда справедлива формула Остроградского Гаусса

∫∫∫

R

dxdydz =

∫∫Rdxdy .

 

W

z

W

При доказательстве этого равенства будем предполагать, что область W выпукла по z ( любая вертикаль пересекает W по отрезку или по пустому множеству). В этом случае область W можно описать, как геометрическое место точек следующего вида:

W = {(x,y,z):z [z1(x,y),z2(x,y)] для любых (x,y) D},

где z1(x,y),z2(x,y) – две непрерывные функции, определенные на D. В этом случае

 

R

 

 

 

z2

( x, y )

R

 

∫∫

 

 

 

dxdydz =

∫∫

dxdy

 

dz =

{R[ x, y, z ( x, y)] − R[ x, y, z ( x, y)]}dxdy =

 

 

 

 

∫∫∫ z

 

 

z

2

1

 

 

 

 

 

 

 

 

 

 

W

D

 

z

( x, y )

 

 

D

 

 

 

 

 

 

 

1

 

 

 

 

 

 

∫∫R( x, y, z)dxdy + ∫∫R( x, y, z)dxdy + ∫∫R( x, y, z)dxdy = ∫∫R( x, y, z)dxdy .

Φ1

Φ 2

 

 

 

 

Φ3

 

W

67

Делая циклические перестановки переменных x→ z→ y, y→ x→ z, z→ y→ x

можно получить еще две формулы для поверхностей выпуклых по другим осям.

∫∫∫

Q

dxdydz =

∫∫Qdzdx , ∫∫∫

P

dxdydz =

∫∫Pdydz .

 

 

W

y

W

W

x

W

 

 

 

 

Если область W удовлетворяет всем трем условиям одновременно и в области задано поле V=(P,Q,R) c непрерывными частными производными па соответствующим переменным, то эти три формулы можно собрать в одну

 

∫∫

 

 

 

 

∫∫∫

P +

Q +

R

 

Pdydz + Qdzdx + Rdxdy .

 

dxdydz =

 

W

x

y

z

W

 

 

 

 

Дивергенция векторного поля определяется по формуле div V =

P +

Q +

R .

 

 

 

 

 

 

 

x

y

z

Тогда, используя векторные обозначения, формулу Остроградского-Гаусса можно записать в виде

∫∫∫

div V dW = ∫∫ (V,dS).

W

W

Формула Остроградского Гаусса будет верна и для областей допускающих разбиение на конечное число областей указанного типа.

Пример 1 (4389). Вычислить I = ∫∫ (x-y+z)dydz+(y-z+x)dzdx+(z-x+y)dxdy, Ф :

Φ

x y + z + y z + x + z x + y = 1 .

u = x y + z v = y z + x

w = z x + y

в системе координат u,v,w

По формуле Остроградского Гаусса I =3 ∫∫∫dxdydz . Сделаем замену переменных

B

u = x y + z

= − +

v y z x , в этом случае в новых координатах граница области будет определяться уравнением w = z x + y

: |u|+|v|+|w|=1. Якобиан отображения равен

68

D(u, v, w)

 

− 1

1

 

D( x, y, z)

 

1

 

3

 

3

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∫∫∫dudvdw =

 

 

2 2

 

=

1

1

−1

= 2 + 2 = =4,

=

, поэтому I =

2

 

=1.

D( x, y, z)

 

 

 

 

 

 

 

− 1

1

1

 

D(u, v, w) 4

4

4

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 2 (4381) Доказать, что если Ф – замкнутая простая положительно ориентированная поверхность, то ∫∫(a, dS )=0, где a - постоянное векторное поле. Утверждение непосредственно

Φ

следует из формулы Остроградского-Гаусса.

Пример 3 (4382). Объем тела равен

W = 1 ∫∫(r, dS ), где W – положительно ориентированная граница области W.

3 W

Утверждение непосредственно следует из формулы Остроградского-Гаусса.

Пример 4 (4383). Объем конуса, ограниченного гладкой конической поверхностью F(x,y,z)=0 и

плоскостью равен W = 1 h Φ , где Φ - площадь основания, h – высота.

3

Поместим начало координат в вершину конуса. Боковую поверхность конуса, ориентированную

внешней нормалью обозначим Ф1 , а основание, ориентированное нормалью m обозначим Ф2 . Тогда

3 W = ∫∫∫3dxdydz =∫∫∫div r dxdydz = ∫∫(r, dS )= ∫∫(r, n)dS )+ ∫∫(r, n)dS Для боковой

W W W Φ1 Φ 2

поверхности конуса скалярное произведение (r, n)= 0 и ∫∫(r, n)dS = 0 . Для поверхности основания

Φ1

конуса (r, n)= (r,−m)= h , поэтому ∫∫(r, n)dS = h∫∫dS = h Φ .

Φ 2 Φ 2

Пример 5 (4390). Вычислить ∫∫(V , dS ), где Ф – часть конической поверхности x2+y2=z2 , 0z a ,

Φ

ориентированной внешней нормалью, а поле V = ( x2 , y2 , z 2 ) . Дополним поверхность до замкнутой.

Основание, ориентированное нужным образом обозначим Ф0 .

69

 

 

 

 

 

 

 

2π

a

h

∫∫(V , dS )= ∫∫∫div V dxdydz = 2∫∫∫(x + y + z) dxdydz = 2 dhr(r cosϕ + r sin ϕ + h)dr =

Φ + Φ 0

 

 

W

 

W

 

0

0

0

2π

a

h

a

h

a

π a4 . ∫∫(V , dS )= ∫∫a2dxdy = a2 µD = πa2 . Таким

2 dhrhdr = 4π hdhrdr = 2π h3dh =

0

0

0

0

0

0

2

Φ 0

 

D

 

 

образом, ∫∫(V , dS )= π a4 πa4

= − π a4 .

 

 

 

 

 

 

Φ

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

§5. Элементы теории поля

1. Введение

Для вектора в будет использоваться обозначение a или a . Функция u(x,y,z) , заданная в области D будет называться скалярным полем. В случае задания трех функций P,Q,R можно говорить о векторном поле V=(P,Q,R). Градиент скалярного поля u , определяется как векторное поле V = grad u

 

u ,

u

,

u

 

=

. Функция u называется потенциалом векторного поля, а само поле называется

x

y

 

z

 

потенциальным. Связь между потенциалом и координатами векторного поля задается соотношением du=Pdx+Qdy+Rdz . Интеграл (V, ds) для замкнутой кривой С называется циркуляцией векторного

C

поля по C. Замкнутая кривая называется контуром, а интеграл по контуру обозначается (V, ds) и

C

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

Определение. Поле V называется соленоидальным, если для него существует векторное поле W такое, что V = rot W. Такое векторное поле W называется векторным потенциалом поля V .

Ранее доказанные утверждения можно сформулировать в виде теоремы

Теорема (Условия потенциальности поля) . Пусть в поверхностно односвязной области D задано непрерывно дифференцируемое поле V=(P,Q,R). Тогда эквивалентны следующие три условия

1. Циркуляция векторного поля (V, ds) равна нулю вдоль любого контура, лежащего в D.

C

2. Поле V потенциальное, т. е. существует дважды непрерывно дифференцируемая функция,

градиентом которой и является данное поле. При этом (V, ds) = u(B) - u(A).

AB

3. Поле V безвихревое.

2.Поток векторного поля

Будем считать, что V=(P,Q,R) – это поле скоростей (рассматривается стационарный поток жидкости). Векторной линией поля V называется линия, касательная к которой в любой точке совпадает по направлению с вектором V .

70