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

Ми10 крив и пов

.pdf
Скачиваний:
8
Добавлен:
09.05.2015
Размер:
2.67 Mб
Скачать

Если функция f (x.y, z) непрерывна вдоль кусочно-гладкой поверхности Q , то поверхностный интеграл первого рода существует.

(Доказательство теоремы смотрите, например,[3]).

Поверхностный интеграл первого рода обладает теми же свойствами,

что и двойной интеграл. В частности, выполняется свойство аддитивности: если поверхность Q разбита на части Q1 и Q2 , то

∫∫ f (x, y, z)dq = ∫∫ f (x, y, z)dq + ∫∫ f (x, y, z)dq.

Q

Q1

Q2

Читателю предлагаем самостоятельно сформулировать остальные свойства.

2.2.1. Вычисление поверхностного интеграла I-го рода

Пусть задана гладкая поверхность Q уравнением z = z(x, y) . Так как поверхность гладкая, то z(x, y) – непрерывная функция вместе со своими частными производными. И пусть на поверхности Q определена непрерывная функция f (x, y, z) . Требуется вычислить поверхностный интеграл ∫∫ f (x, y, z)dq (Нельзя смешивать интегрируемую функцию f (x,

Q

y, z) с функцией, входящей в уравнение поверхности Q ).

Предварительно займемся выводом формулы для вычисления площади qi элементарной части qi поверхности Q . Обратимся к рис. 2.2.2 и 2.2.3, где выделен участок qi разбиения области Q на элементарные части с выбранной на нем точкой Mi (xi , yi , z(xi , yi )) . Через точку Mi qi проведем касательную плоскость к поверхности Q . Уравнение касательной плоскости, как известно (формула (2.7)) , имеет вид

zx' (x, y) |Mi (x xi ) + z'y (x, y) |Mi ( y yi ) (z zi ) = 0.

Рис.2.2.2.

71

На этой плоскости выделим элемент qi с площадьюqi (рис. 2.2.3), который проектируется на плоскость xOy в ту же элементарную область σi , что и элемент qi (рис. 2.13). Заменим криволинейный элемент qi

плоским элементом qi , тогда

qi ≈ ∆

qi

.

(2.11)

Рис.2.2.3

Обозначим через γi двугранный угол между касательной плоскостью и плоскостью xOy. Следовательно, σi = ∆qi | cosγi | ( как известно из

аналитической геометрии, площадь проекции плоской фигуры равна площади самой этой фигуры, умноженной на абсолютную величину косинуса двугранного угла между плоскостями). Откуда в силу формулы (2.11) получаем

q

σi

.

(2.12)

 

i

| cosγi |

 

 

 

Линейный угол двугранного угла γi есть в то же время угол между осью

Oz и единичным вектором нормали n к касательной плоскости. И поэтому, см. формулу (2.8),

| cosγi

|=

 

 

 

1

.

(2.13)

 

 

1 +

(zx' (xi , yi ))2 +(z'y (xi , yi ))2

 

 

 

 

 

 

 

 

Таким образом, с учетом (2.12) и (2.13) получаем

 

 

qi

σi

 

 

=

1

+(zx' (xi , yi ))2 +(z'y (xi , yi ))2 σi .

(2.14)

| cosγi |

 

 

 

 

 

 

 

А теперь, для решения поставленной задачи о вычислении поверхностного интеграла, вернемся к интегральной сумме (2.9), соответствующей данному разбиению поверхности Q на части qi и выбору точек

72

Mi

( i =1, 2,..., n).

Принимая во внимание полученное выражение для qi

(2.14), запишем формулу (2.10) в виде:

 

n

n

∫∫f(x, y,z)dq=limd0 f(xi, yi,zi )qi =limd0 f(xi, yi,z(xi, yi )) 1+(zx' (xi, yi ))2 +(z'y(xi, yi ))2σi.

Q

i=1

i=1

Очевидно, если d 0 , ( d — максимальный из диаметров di частичных

областей qi ,i =1,2,..., n,) то и δ 0 , где δ – максимальный из диаметров

проекций σi , i =1, 2,..., n, частичных областей qi ,

на плоскость xOy (рис.

2.13). Отсюда

n

 

 

 

∫∫ f (x, y, z)dq = limδ0

f (xi , yi , z(xi , yi )) 1+(zx' (xi , yi ))2

+(z'y (xi , yi ))2 σi . (2.15)

Q

i=1

 

Сумма, стоящая в правой части равенства (2.15), является интегральной

суммой

для

непрерывной

функции

f (x, y, z(x, y))

1 +(zx' (x, y))2

+(z'y (x, y))2

по плоской

области Sxy

проекции поверхности Q

на плоскость

xOy . Поэтому

ее предел есть

двойной интеграл от указанной функции двух переменных по области Sxy . Следовательно,

∫∫ f (x, y, z)dq = ∫∫ f (x, y, z(x, y)) 1 +(zx' (x, y))2 +(z'y (x, y))2 dσ,

Q

Sxy

иначе

 

∫∫ f (x, y, z)dq = ∫∫ f (x, y, z(x, y)) 1 +(zx' (x, y))2 +(z'y (x, y))2 dxdy. (2.16)

Q

Sxy

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

Если

гладкая поверхность

Q задана

уравнением y = y(x, z) или

x = x( y, z) , то проектируя поверхность Q ,

соответственно, на плоскости

xOz и yOz , из равенства (2.16), можно получить следующие формулы:

∫∫ f (x, y, z)dq = ∫∫ f (x, y(x, z), z))

1 +(yx' (x, z))2

+(yz' (x, z))2 dxdz,

(2.17)

Q

Sxz

 

 

 

 

∫∫ f (x, y, z)dq = ∫∫ f (x( y, z), y, z))

1 +(x'y ( y, z))2

+(xz' ( y, z))2 dydz.

(2.18)

Q

Syz

 

 

 

 

где Sxz , Syz – проекции поверхности Q на координатные плоскости xOz и yOz .

Замечание. В более сложных случаях, когда поверхность кусочно-гладкая или неоднозначно проектируется на координатные плоскости, ее

73

разбивают на гладкие части, однозначно проектирующиеся на одну из координатных плоскостей, тогда интеграл разобьется на сумму интегралов по этим частям, к каждому из которых применимы формулы (2.16) – (2.18). Пример 2.3. Вычислить площадь Q части гиперболического параболоида

z = xy , вырезаемой прямым круговым цилиндром x2 + y2 = 8 (рис.2.2.4).

Рис.2.2.4

 

 

 

 

 

 

Решение. В данном случае поверхность

является

графиком

функции

z(x, y) = xy с областью определения

S

xy

=

{

 

}

В этой

 

 

(x, y) :x2 + y2 8 .

области функция z(x, y) = xy непрерывна и имеет непрерывные частные производные zx' = y и z'y = x . Следовательно имеем (формула (2.10*)):

SQ = ∫∫dq = ∫∫

1 +(zx' (x, y))2 +(z'y (x, y))2 dxdy =∫∫ 1 + y2 + x2 dxdy.

Q

Sxy

Sxy

Для вычисления этого двойного интеграла удобно перейти к полярным координатам x = ρ cosϕ, y = ρsinϕ , в которых

dxdy = ρdρdϕ, Sxy

D ={(ρ,ϕ) : ρ [0, 8],ϕ [0,2π]}.

 

 

Тогда получим

 

 

 

 

 

 

 

 

 

SQ = ∫∫ 1 + y2 + x2 dxdy =

2π dϕ 8

ρ 1+ ρ2 dρ = 2π

(1 + ρ2 )3/ 2

 

 

8

52

π.

 

 

 

=

3

3

Sxy

0

0

 

 

 

0

 

 

 

 

 

 

 

2.2.2. Механические приложения поверхностного интеграла первого рода

Пусть по поверхности Q распределена масса с поверхностной плотностью ρ(x, y, z) , представляющей собой непрерывную функцию на Q . Такую поверхность называют материальной поверхностью. Разобьем поверхность на n частичных областей qi ,i =1, 2,..., n, выбрав на каждой из них произвольную точку Mi (xi , yi , zi ) qi . При этом полагаем, что при малых диаметрах di частичных областей qi поверхностная плотность

74

распределенной по поверхности массы в пределах каждой частичной области постоянна и равна значению ρ(Mi ) = ρ(xi , yi , zi ) .

При этих предположениях, как и в случае материальной кривой (см. 1.2.3.), приходим к следующим выводам:

а) масса всей поверхности Q равна

 

 

mQ = ∫∫ρ(x, y, z)dq;

 

 

 

(2.19)

 

 

 

Q

 

 

 

 

 

 

 

б) статические моменты этой поверхности относительно

 

координатных плоскостей yOz, xOz , и xOy равны

 

 

 

M yOz = ∫∫xρ(x, y, z)dq, MxOz = ∫∫yρ(x, y, z)dq, MxOy = ∫∫zρ(x, y, z)dq,

(2.20)

Q

 

Q

 

 

 

 

Q

 

где x, y, z – координаты точки

M Q;

 

 

 

 

 

в) координаты центра масс поверхности Q имеют вид

 

x

=

M yOz

, y =

M

xOz

, z =

MxOy

 

;

(2.21)

 

 

 

 

 

C

 

mQ

C

mQ

C

mQ

 

 

 

 

 

 

 

 

 

 

г) моменты инерции этой поверхности относительно, например, плоскости yOz , оси Ox и начала координат равны соответственно

JyOz = ∫∫x2 ρ(x, y, z)dq, JOx = ∫∫( y2 + z2 )ρ(x, y, z)dq,

Q Q

JO = ∫∫(x2 + y2 + z2 )ρ(x, y, z)dq.

Q

Аналогично выражаются моменты инерции относительно других осей. Пример 2.4. Вычислить статический момент относительно плоскости xOy и координаты центра масс однородной полусферы радиуса R .

Решение. Выберем систему координат так, чтобы полусфера стояла на плоскости xOy , а начало координат находилось в ее центре. Тогда

уравнение полусферы имеет вид x2 + y2 + z2 = R2 (z 0) .

На основании формулы (2.20), где полагаем ρ(x, y, z) =1, так как

полусфера однородная, имеем MxOy

= ∫∫zdq.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Разрешим уравнение поверхности относительно z : z =

 

R2 x2 y2 Тогда

поверхность

является

графиком

функции

 

z(x, y) =

R2 x2 y2

с

областью определения

Sxy ={(x, y) :x2 + y2

R2 }.

 

На

основании

с

формулы (2.16) получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

MxOy = ∫∫zdq = ∫∫ R2 x2 y2

1+

 

 

 

x

2

 

 

 

+

 

 

y

2

 

 

dxdy =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

2

x

2

y

2

R

2

x

2

y

2

 

Q

Sxy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

75

= ∫∫ R2 x2 y2

 

 

R

 

 

 

dxdy = R∫∫dxdy =R πR2 =πR3 .

R

2

x

2

y

2

Sxy

 

 

 

Sxy

Так как материальная сфера однородная, то из соображений симметрии ее центр тяжести должен находиться на оси Oz , поэтому xC = 0, yC = 0 , а

zC = MxOy . По формуле (2.19) найдем массу полусферы mQ

mQ = ∫∫dq = ∫∫

1 +

 

 

x2

 

 

+

 

 

y2

 

 

dxdy =∫∫

 

 

R

 

 

 

dxdy.

R

2

x

2

y

2

R

2

x

2

y

2

R

2

x

2

y

2

Q

Sxy

 

 

 

 

 

 

 

 

Sxy

 

 

 

 

Для вычисления двойного интеграла перейдем к полярным координатам: x = ρ cosϕ, y = ρsinϕ , в которых

dxdy = ρdρdϕ, Sxy D ={(ρ,ϕ) : ρ [0, R],ϕ [0, 2π]}.

Тогда

 

 

 

R

 

 

 

dxdy = R2π dϕR

 

 

ρdρ

 

= −R2π

 

R

mQ = ∫∫

 

 

 

 

 

 

 

 

R2 ρ2

dϕ = 2πR2 .

 

2

2

 

2

 

 

 

2

2

Sxy

R

 

x

y

 

 

MxOy

 

0 0

 

R ρ

 

0

 

0

 

 

 

 

 

 

Таким образом,

z =

 

=

πR3

=

 

R

.

 

 

 

 

 

m

2πR2

2

 

 

 

 

 

 

 

 

C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q

 

 

 

 

 

 

 

 

 

 

 

 

 

2.3. Поверхностный интеграл второго рода

 

 

 

 

 

 

 

 

2.3.1. Векторные поверхности

 

G

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пусть в пространственной области Q задано векторное поле a(M ) .

Возьмем в

Q гладкую кривую

L, для которой в любой ее точке

M

касательныйG вектор кривой не является коллинеарным вектору

G

a(M )

частности,

a(M ) 0G). Такую кривую будем называть трансверсалъной

векторному

полю a(M ) . Через

каждую точку M L можно

провести

векторную линию рассматриваемого поля. При этом множество всех таких линий, проходящих через точки M L , образует поверхность S G(Рис.

2.3.1), называемую векторной поверхностью векторного поля a(M ) .

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

к векторной линии в точке M S , а значит, и вектор a(M ) будут

JJJJG

перпендикулярны единичному вектору n(M ) нормали к S в этой точке

76

Рис. 2.3.1

Выбрав в качестве кривой L простой замкнутый контур, получим векторнуюG поверхность, называемую векторной трубкой векторного

поля a(M ) (Рис. 2.3.2). Любая векторная линия, не проходящая через точки контура L, целиком лежит либо внутри, либо вне векторной трубки.

Рис. 2.3.2

2.3.2. Задача о потоке векторного поля. Понятие поверхностного интеграла II-го рода

Пусть в пространственной области задано векторное поле

JJJJJJG

G

G

G

F(M ) = X (x, y, z) i

+Y(x, y, z) j +Z(x, y, z) k , где X = X (x, y, z);

Y =Y(x, y, z);

Z = Z(x, y, z) – непрерывные функции.

Рис. 2.3.3

77

На гладкой двусторонней поверхности Q , ограниченной гладким замкнутым контуром L, выберем одну из сторон поверхности Q с

помощью единичного

вектора

нормали

n(M ) к этой поверхности.

Координатами вектора

G

являются

его направляющие косинусы

n(M )

cosα,cos β,cosγ , определяемые формулой (2.8) (Рис.2.3.3) .

JJJJJJG

Вектор F(M ) может определять направление скоростей: потока жидкости,

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

JJJJJJG

определенности, вектор F(M ) определяет скорость потока жидкости, протекающей через поверхность Q . Разобьем поверхность Q

произвольным образом на n элементарных частей, площадь которых обозначим через qi . На каждой из частей qi выберем произвольную

точку

 

 

M

(x , y , z ) ,

и вычислим

значение

вектора

скорости

жидкости

JJJJJJJG

 

i

i i

i

Если выбранное разбиение поверхности Q достаточно

F(M

)

 

(Рис.2.3.4).

i

 

 

 

 

 

 

 

 

 

 

что значение вектора скорости

JJJJJJG

на

мелкое, то можно допустить,

F(M )

площадке

qi

 

постоянное

 

и

 

 

JJJJJJG

 

Рассмотрим

скалярное

 

 

равно F(Mi ) .

 

 

 

 

 

 

 

JG

 

G

 

))

векторов

 

JJJJJJG

и

n(M

) .

Так

как

произведение (F(M

), n(M

 

F(M

)

JG

 

 

G

 

)) | q

i

 

i

 

 

 

 

i

 

 

i

 

 

 

| (F(M

 

), n(M

можно рассматривать как объем призмы с площадью

 

 

i

 

i

 

i

 

 

 

G

G

)) | (Рис.2.3.4), то произведение

основания

q

и высотой

| (F(M

), n(M

JG

 

 

G

 

i

 

 

 

 

 

i

i

 

 

 

 

 

 

 

 

(F(Mi ), n(Mi ))qi

приближенно

определяет

 

количество

жидкости,

протекающей через площадку qi за единицу времени в направлении вектора нормали n(Mi )

Рис.2.3.4

78

Тогда общее количество жидкости, протекающей в единицу времени через

поверхность Q в направлении вектора нормали n(M ) к поверхности Q ,

JJJJJJG

что естественно назвать потоком PQ

векторного поля F(M )

через

n

JG

G

 

поверхность Q , приближенно равно (F(Mi ), n(Mi ))qi , то есть

 

i=1

G

 

 

n JG

 

 

PQ (F(Mi ), n(Mi ))qi .

(2.20)

i=1

Сумма (2.20) называется интегральной суммой. При различных разбиениях поверхности Q на элементарные части qi ,i =1, 2,..., n, и различном выборе

точек Mi

qi

получаем различные значения интегральной суммы (2.20).

Пусть d

максимальный

из

диаметров

di частичных

областей

qi ,i =1, 2,..., n.

Естественно

в

качестве

точного значения

потока

векторного поля принять предел интегральных сумм в правой части (2.20)

при d 0,

т.е.

 

JG

G

JG

G

 

 

 

n

 

PQ

= limd 0

(F(Mi ) n(Mi ))qi

= ∫∫(F(M ) n(M ))dq,

(2.21)

 

 

i=1

 

 

Q

 

 

который не зависит ни от способа разбиения, ни от выбора точек Mi qi .

Интеграл, стоящий в правой части формулы (2.21), как известно, представляет собою поверхностный интеграл первого рода.

Замечание 2.3.1. Если поверхность Q является поверхностью векторной

JJJJJG

трубки

векторного поля F(M ) , то

поток векторного поля через

поверхность Q

 

равен нулю, так как, в этом случае,

JJJJJJG

 

векторы F(M ) и

G

 

 

 

 

 

 

 

 

 

n(M ) ортогональны.

 

 

 

 

Таким образом, нахождение потока векторного поля привело нас к

нахождению предела интегральных сумм вида (2.20) .

 

 

 

 

G

 

 

 

 

 

 

 

Так

какG n(Mi ) = (cosαi ,cos βi ,cosγi ) ,

где αi , βi и

γi – углы между

вектором

n(M

i

)

и направлениями осей

Ox,Oy и Oz соответственно,

JJJJJJJG

 

 

 

 

 

 

 

F(Mi )

= (X (Mi ),Y(Mi ), Z(Mi )) , то

 

 

 

n

JG

G

 

 

 

n

 

 

lim

(F(Mi ) n(Mi ))qi =lim(X(Mi )cosαi +Y(Mi )cos βi

+ Z(Mi )cosγi )qi =

d0

i=1

 

 

 

 

d0

i=1

G

G

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=∫∫(X(x, y, z)cosα +Y(x, y, z)cosβ +Z(x, y, z)cosγ)dq =∫∫(F(M) n(M))dq. (2.22)

Q

 

 

 

 

 

 

 

Q

 

79

n

Рассмотрим предел интегральной суммы lim Z(Mi )cosγiqi . Как

d 0 i=1

известно, угол γi между вектором n(Mi ) и осью Oz равен углу между площадкой разбиения qi и плоскостью xOy (рис. 2.2.3).

 

 

Рис.2.3.5

Тогда cosγiqi

= ∆si , где | si | – площадь проекции si площадки разбиения

q в плоскость

xOy (рис. 2.3.5). Здесь s берется со знаком плюс, если

i

G

i

угол между вектором n(Mi )

и осью Oz острый, и со знаком минус, если

угол тупой (рис. 2.3.5). Так как предел интегральной суммы (2.21), если он существует, не зависит от способа разбиения, то мы можем разбить поверхность Q плоскостями, параллельными координатным плоскостям.

Тогда cosγiqi = ∆si = ∆xiyi и

n

lim

d 0 i=1

n

Z(Mi )cosγiqi = lim Z(xi , yi , zi )xiyi . (2.22*)

d 0 i=1

Если

предел

интегральных

сумм

 

n

 

,

где

d

lim Z(xi , yi , zi )xiyi

 

 

 

 

d 0

i=1

 

 

 

 

 

 

 

di

 

 

qi ,i =1, 2,..., n,

максимальный

из диаметров

частичных

областей

существует, не зависит от способа разбиения поверхности Q

и выбора

точек

Mi (xi , yi , zi ) , то такой предел называется поверхностным

интегралом второго рода ∫∫Z(x, y, z)dxdy. То есть

 

 

 

 

 

Q

 

 

 

 

 

 

 

 

 

n

 

 

 

n

 

 

 

∫∫Z(x, y, z)dxdy = limd 0 Z(xi , yi , zi )xiyi = limd 0 Z(Mi )cosγiqi

=

Q

 

i=1

 

 

 

i=1

 

 

 

 

 

= ∫∫Z(x, y, z)cosγdq.

 

 

 

(2.23*)

 

 

Q

 

 

 

 

 

 

 

80