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

книги из ГПНТБ / Сахарников Н.А. Высшая математика учебник

.pdf
Скачиваний:
98
Добавлен:
27.10.2023
Размер:
24.22 Mб
Скачать

Пользуясь известным выражением длины дуги линии, заданной параметрически (см. и. 159):

U

« = j V х'и2 + у'и2 z[2du, и0

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

dsu ~ V х ’и*

zz 2du.

(12)

Таким образом, получим

 

 

dsu = Hudu,

dsv = Hv du,

dsw~-=Hw'dw,

(13)

где

____________

 

Н и =

Y х'и 2 + У и 2 + Z u 2,

 

Н ѵ =

Ѵ х ' ѵ 2 + У ѵі +

z'v 2,

 

Hw= V x w2-\ y'wz + Zw2.

(14)

Величины Hu, Hu и Hw, определяемые формулами (14), назы­ ваются коэффициентами Ламе * системы координат (и, v, w).

Эти величины суть функции и, ѵ и w. Каждый коэффициент Ламе характеризует скорость изменения длины дуги соответствующей

координатной линии. Например,

Н и —

ds

 

 

 

 

 

- ~ - .

 

 

 

 

 

Вычислим коэффициенты Ламе в двух важных случаях. В слу­

чае цилиндрической

системы

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

(см. п. 97) х =

=

р cos ф, у = р sin ф, z = z и,

полагая

и =

р,

ѵ = ф , w — z,

по формулам (14) получим

 

 

 

 

 

 

 

 

 

 

 

 

 

Яр = 1,

Я ф = р,

Я г =

1.

 

 

 

(15)

 

В

случае сферической

системы

координат

имеем

(см. и. 97)

X = г sin Ѳcos ф,

у =

г sin Ѳsin ф,

z

г cos Ѳ

и, полагая: и =

=

г,

v = ф, w = Ѳ,

по формулам

(14)

получим

 

 

 

 

 

Hr= 1,

Яф = г 8іпѲ,

Яѳ = г.

 

 

(16)

 

Э л е м е н т

о б ъ е м а в

к р и в о л и н е й н ы х к о ­

о р д и н а т а х .

Рассмотрим

трехмерную

область

(рис. 123),

ограниченную тремя

парами

координатных

поверхностей: и =

=

и0,

и = и0 +

du,

v — ѵ0,

v =

ѵ0 -f

dv,

w =

w0,

w == w0 +

+

dw. При достаточно малых | du |, | dv |, | dw | эту область в первом

приближении можно принять за параллелепипед с ребрами, соот­ ветственно равными dsu, dsv, dsw и направленными по касатель­ ным к координатным линиям. Под элементом объема в системе

Габриель Ламе (1795—1870) — французский математик.

координат (и, v, w) понимают объем этого

параллелепипеда.

В частном случае, когда система координат (и, v,

w) ортогональна,

упомянутый параллелепипед — прямоугольный и его объем равен произведению его измерений: dx = dsudsvdsw. С помощью формул (13) получаем выражение элемента объема в ортогональной си­ стеме координат

dx = НиНѵНа du du dw.

(17)

В о б щ е м с л у ч а е координатные линии могут

пересе­

каться не под прямым углом. Обозначим через еи, ev, ew

единич­

ные векторы касательных к соответствующим координатным ли­ ниям в рассматриваемой точке М 0, направленные в сторону возра­ стания координат (см. рис. 123). В общем случае элемент объема

есть объем

параллелепипеда, ребрами ко­

w)

торого

служат

 

векторы

dsueu,

dsvev,

 

 

dswew.

Он

равен

абсолютной

величине

 

смешанного

произведения

этих

векторов

 

(см. п.

96):

 

 

 

 

 

 

 

 

 

d x = ± ([d sueu,

dsv ev],

dswew).

(18)

 

Рассмотрим

радиус-вектор г =

х

і +

 

+ У j + zk, проекции которого х,

у,

z свя­

 

заны с и, v,

w равенствами (7). Дифферен­

 

цируя г пои,

v,

w,

получим три

вектора:

 

 

 

 

 

О

I У и 3 I-

к,

 

 

v 1 4 Уѵ 3 + Zv k,

(19)

 

 

 

 

 

:

Xraj 1

Уѵі

 

 

направленные по касательным к соответствующим координатным линиям (см. и. 185). Длины векторов (19) соответственно равны значениям коэффициентов Ламе:

I г'иI = Y х'и 4- Уи24- Zu =

Ни, IК\ =нѵ, Iг ; I = Hw.

(20)

Поэтому согласно (13) имеем

 

 

= Ние„ = ги,

sveu— r*,, swеш™ Гу,.

 

С помощью равенств (19) получим

 

dsueu = x'ud u i- r y'udu] i~zuduk,

 

dsv ea = x ' v d v i - t y ' v d v i - ^ Z v d v k ,

(21)

dsœ ew=•= x'w dw i (-y'wdw j 4~ z'-^dwk.

 

Из (18) и (21), учитывая представление смешанного произведения векторов через проекции сомножителей (см. п. 96), следует выра­ жение элемента объема в системе координат (и, v, w)

где

n = D (*• У’ 2)

Хц

Хѵ

X'w

У и

Уѵ

y'w

D (и, и, w)

 

Zu

Zv

Zw

Геометрически функциональный определитель D представляет отношение, в котором преобразование (7) изменяет объем.

179. Вычисление тройного интеграла в ‘.(криволинейных коор­ динатах. В цилиндрической системе координат координатными поверхностями, проходящими через точку М 0, являются: цилинд­ рическая поверхность р = р0, полуплоскость ф = ф0 и плоскость

z

— z0. Три пары координатных поверхностей р = р0, Р = Ро +

+

dp, ф =

Фо, ф = Фо + ^Ф.

2 =

z0, z =

z0 +

dz ограничивают

 

і, 1

область

(рис. 124),

которую можно принять в пер­

 

вом приближении за прямоугольный

параллеле­

 

 

пипед с ребрами, длины которых суть

М йА — dp,

 

 

MgB

pdф, MgC = dz.

Поэтому

элемент объема

 

 

цилиндрической системе

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

 

 

 

 

dт = pdpdфdг.

 

(23)

 

 

Этот же

результат получается

по формулам (15)

 

 

и (17) или (22), так как D = р.

 

 

 

Для

вычисления

тройного

интеграла (2) в

 

 

цилиндрической системе координат следует раз­

 

 

бить область В на

элементы координатными по­

 

 

верхностями

системы

и

составить

интеграль­

 

ную сумму (1). В результате предельного пере­

Рис. 124.

хода при

Хп -> 0 с учетом формулы (23) получим

Ш

/< « .

у> z )d t= lim 2 / ( æ*. Уkt zfe)Axfe =

 

В

 

 

Я„->0 в

 

=

Hm

V. / (pÄc°s фй, р^ьіпф*, zk)pk txpk ^ k l4zk.

 

Ь п - * 0

в

 

 

Окончательно

имеем

}Jj/(pсоэф, рэіпф, z) p d p d ф d z .

 

Ш / ( * . У'

z)d T =

(24)

в

 

 

в

 

Если область В ограничена координатными поверхностями си­ стемы (w, V, w), то соответствующий повторный интеграл будет иметь постоянные пределы для каждой из переменных интегриро­ вания.

П р и м е р 1. Вычислить объем тела, вырезанного цилиндром я2-(-у2 = 1

из шара

-|-!/2 _ [ _ z

2

^ 4 .

Решение

 

 

2

П

1

Ѵ і - р г

1

^ =

ЙТ = 2 I

 

^

J

dz = 4л j р V 4 р2 ір = ~ - ( 8 — 3,/г).

н

о

 

о

о

о

В сферической системе координат координатными поверхно­

стями системы,

проходящими через М 0, являются сфера г = г0,

полуплоскость ф = ф0 и коническая поверхность Ѳ =

Ѳ0.

Три пары координатных поверхностей г = г0, г ~~ r 0+

dr, ф = ф0,

Ф = ф о + d(p,

Ѳ = Ѳ0, Ѳ = Ѳ0 + dQ ограничивают область

(рис. 125), которую можно принять приближенно за прямоуголь­

ный параллелепипед с ребрами, длины которых суть М 0А

= rdQ,

М 0В =

г sin0 гіф, М 0С — dr. Поэтому элемент объема в сфериче­

ских

координатах

равен

 

 

 

 

 

 

 

dx — г2sin Ѳdr Ар dQ.

(25)

Тот же

результат

получается

по форму­

с

лам

(16)

и

(17), а также

по формуле (22),

 

так

как

D = г2 sin Ѳ.

 

в

тройном

 

Замена

переменных

 

интеграле (2) приводит согласно (25) к

 

равенству^1

 

 

 

 

 

 

 

 

 

Ш /( * , У,

z)dx =

 

 

 

IJ j /

В

 

 

 

 

 

 

sin Ѳcos ф,

г sin Ѳsin ф, r cos Ѳ) r2sin Ѳdr dq> dQ.

(26)

 

в

 

 

 

 

 

 

 

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

П р и м е р 2. Вычислить объем шара радиусом R. Но формуле (26) имеем

 

2 Л

Я

й

 

V = ^ Л j

dx = j’ Лр

^

sin ѲйѲ ^ /-2 dr = 2л • 2 • ~

л№.

В о б щ е м

с л у ч а е

любой криволинейной системы коор­

динат (и, V, w), связанной с декартовой системой (х, у,

z) равенст­

вами (7), элемент объема представлен формулой (22). Замена пере­

менных (7) в тройном интеграле

(2)

приведет к равенству

/ / j /(z,

У,

z)dx =

= I j / / ( : r (н, и, w), у (и, и, w), z(u, v, w))\ D \dudv dw. (27)

в

Здесь правые части (7) предполагаются непрерывно дифферен­ цируемыми в В, якобиан D — сохраняющим знак в В, f (М ) —

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

Если область В ограничена координатными поверхностями системы (и, V, w) , то в соответствующем повторном интеграле пределы получатся постоянными.

П р и м е р

3.

Вычислить

объем

тела, ограниченного эллипсоидом

Р е ш е н и е .

Вводим обобщенные

сферические координаты (г, ф, Ѳ),

связанные с декартовыми координатами равенствами

X — ar sin 0 cos ф,

у = br sin Ѳsin ф, z = wsin0.

В новой системе координат эллипсоид представляет координатную по­

верхность г — 1.

Якобиан преобразования координат в нашем случае равен

Ü = aber2 sin Ѳ. Поэтому искомый

объем равен

 

 

2п

л

г

V = J ^

^ dx = abc J

 

J

sin Ѳ dO ^ r~dr~ —nabe,

 

в

o

o

 

о

180. Приложения тройного интеграла в механике. Рассмотрим трехмерную область В, в которой непрерывно распределено вещество с плотностью р (М ). Разделим В на элементарные части. Элементарная масса равна dm = р dr.

Элементарные статические моменты относительно координат­ ных плоскостей определяются равенствами

dSxy= zdm, dSyz^xdm , dSxz = y dm.

Поэтому для массы и статических моментов всего тела получим

следующие

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

представления:

 

 

 

 

 

т =

Я1рd x ' s

* y =111zpd r ’

 

 

 

 

 

 

В

 

 

В

 

 

(28)

 

 

Syz =

J J J гр dr,

SX2 -- J j J f/p dr.

 

 

 

 

 

в

 

 

в

 

 

 

Координаты

центра

тяжести

С

(х с ,

у с ,

zc)

тела согласно

определению

этого

 

понятия

удовлетворяют

соотношениям

тхс = SyZ,

myc ~ S xz,

mzc ~ S xy, поэтому имеем

 

 

 

■ —

т

У с

$хг

$ху

(29)

 

 

С

т

т

»

 

где величины,

входящие

в правые

части,

определяются форму­

лами (28).

 

 

 

 

 

 

 

 

 

 

Моментом инерции системы точек относительно оси назы­ вается величина, равная сумме произведений масс этих точек на квадраты их расстояний до оси. Элементарные моменты инер­ ции относительно координатных осей в общем случае равны

dJx — {y2-\-z2)dm, dJy — (x2~r z2)dm, dJz = (x2 -{-у2) dm,

где y2-\-z2, X2 + z2 и X2 + У2 суть квадраты расстояний от точки М (х, у , z) соответственно до осей Ох, Оу и Oz. Поэтому моменты

инерции всего

тела

равны

 

 

 

/* = Ш

(j/M-z2)pdr,

Jy = \ ^ { x 2~'r z2)ç)dx,

 

 

В

 

 

в

 

 

 

j j j ( х 2 +

У 2) P d x -

 

 

 

 

в

 

 

 

 

 

П р и м е р .

Вычислить

момент

инерции

шара ж2 + у2 + z2

R

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

 

Р е ш е н и е .

Вычислим интеграл Jz в сферических координатах. Полу­

чим

 

 

й

2 Я

Я

 

 

 

 

 

Jz = J J J г4 sin3 Ѳdr dcp dd— J r* dr Jdcp

Ц sin3 ѲdG = -^r лЯ5.

 

Глава X I

ВЕКТОРНЫЙ АНАЛИЗ И ТЕОРИЯ ПОЛЯ

181. Понятие поля. Понятие поля лежит в основе многих представлений современной физики. Векторный анализ является средством изучения полей.

При рассмотрении физических явлений встречаются величины трех видов: скалярные, векторные и тензорные. Соответственно этЛіу различают три вида полей: скалярные, векторные и тензор­ ные. В общем случае полем величины а называется облаетъ В трех­ мерного пространства, с каждой точкой которой в каждый мо­ мент времени связано определенное значение величины а. Поэтому а есть функция точки М и времени t. Заметим, что областью В может быть и все трехмерное пространство. Если величина а скалярная (температура, давление, потенциал и т. д.), то и поле называется скалярным. Если а есть вектор (сила, скорость и т. д.), то соответствующее поле называется векторным. Если а — тензор (напряженности, проводимости и т. п.), поле называется тензорным.

Поле величины ос называется стационарным, или установи­ вшимся, если ос не зависит от времени t. В противном случае оно называется неустановившимся.

Пример с к а л я р н о г о п о л я . Пусть данная область заполнена газом, температура которого есть функция точки и времени. Однако со вре­ менем может наступить тепловое равновесие, в результате которого темпера­ тура в каждой точке области стабилизируется во времени. Температура станет функцией лишь точки, поле температур станет стационарным.

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

2. Пусть в некоторой области пространства распределено некоторое вещество. Тогда на материальную точку с единичной массой действует сила

тяготения, зависящая от местоположения точки. Эти силы образуют вектор­ ное поле, называемое полем тяготения, или гравитационным полем, соответ­ ствующим данному распределению масс.

3. Если в некоторой области пространства распределены электрические заряды, то на единичный электрический заряд, находящийся в точке М, эти заряды действуют с некоторой силой Е (М). Эти силы образуют вектор­ ное поле, называемое электростатическим, соответствующим данному рас­ пределению зарядов.

Примером т е н з о р н о г о п о л я является поле тензора упругих напряжений тела, с каждой точкой которого связапо определенное значение тензора напряжения. При этом в каждой системе координат тензор в данной точке можно представлять системой девяти чисел (в то время как скаляр можно представить одним числом, а вектор — системой трех чисел).

§31. СКАЛЯРНЫЕ И ВЕКТОРНЫЕ ПОЛЯ

182.Скалярное поле. Поверхности и линии уровня. Пусть дано стационарное скалярное поле <р (М). Если выбрать систему декартовых координат , у, z), то каждая точка М будет иметь координаты X, у, z и функция точки cp (М) станет функцией трех независимых переменных ф (х, у, z). Мы всегда будем предпола­ гать эту функцию непрерывно дифференцируемой в рассматри­ ваемой области.

Основным вопросом исследования скалярного поля является вопрос об изменении функции ф при переходе из одной точки пространства в другую. Для выяснения этого вопроса рассмотрим прежде всего геометрическое место точек, в которых величина ф сохраняет постоянное значение, оно называется поверхностью уровня. Через каждую точку М 0 области В проходит единственная поверхность уровня. Ее уравнение в выбранной системе коор­ динат имеет вид

ф(ж, у, z) -----ф(х0, у0, д„) г. , (1)

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

Указанный способ изображения поля особенно удобен, если речь идет о поле, заданном в плоской области А двух перемен­ ных. Такое поле описывается функцией двух переменных ф (х, у), определенной в области А. Равенство ф (х, у) = с определяет, вообще говоря, некоторую кривую, называемую линией уровня плоского скалярного поля. С помощью линий уровня изобра­ жается, например, рельеф местности на топографических картах (путем задания поля высот точек местности над уровнем моря). В случае поля температур на плоскости линии уровня называ­ ются изотермами, в случае поля давлений — изобарами и т. д.

183.Производная по направлению. Проследим за изменением

данной скалярной величины

(М)

при

перемещении точки М

в каком-либо заданном направлении

m (cos a, cos ß,

cos 7)

(рис. 126).

 

 

 

 

 

О п р е д е л е н и е . Производной

скалярной

функции

cp (М)

по направлению m в точке М

называется

предел

(если он сущест­

вует) отношения приращения Дф функции при смещении точки М в направлении вектора т к величине

іі

ft

^

 

этого

смещения р = d (М , N), когда

 

1---------——*■ "

 

последнее стремится к нулю; она

 

Рис. 126.

 

 

обозначается символом

 

 

 

 

 

 

■!*- =

П т

^

= 1 іт

(2)

 

дт

p -о

Р

p- о

Р

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

влений,

поэтому

 

функция ф (М ) в каждой точке имеет бес­

конечное

множество

производных,

соответству­

 

ющих различным направлениям.

 

по

направ­

 

Понятие

производной

функции

 

лению

является

 

обобщением

понятия

частной

 

производной,

потому

что

частные

 

производные

 

функции

ф (X, у,

 

z)

являются

производными

 

этой функции в направлении соответствующих

 

координатных

осей.

 

 

 

 

 

 

 

 

Для вывода формулы производной по направ­

 

лению,

удобной

в

вычислительном

отношении,

 

обозначим

через

ж,

у и z

координаты

фиксиро­

 

ванной точки М. Тогда точка N будет иметь коор­

Рис. 127.

динаты х~г р cos а,

у + р cos ß,

z +

р cos у. Вели­

 

чины X ,

у, z, cos а,

 

cos ß

и cosy

фиксированы,

р. Обозна­

поэтому

ф (N) является функцией

 

только смещения

чим

эту

функцию

 

через

 

 

 

 

 

 

 

 

 

ф (р) =

 

ф(ж-:-рсо8<х, y-f-pcosß, z-f-pcosy).

(3)

При

р = 0 имеем ф(0) = ф(ж, у,

z). Следовательно, производная ф

по направлению m равна производной от ф по р в точке р —0:

ІФ

= И т

№ г * ( . 0> = ф*(0).

(4)

° т

о+О

Р

 

Дифференцируя ф(р) по р как сложную функцию, получим

ф' (р) = ц)'х (N) cos а + фу (N ) cos ß -f- фг (N) cos у.

Положив р - - 0, получим окончательно

д(р

ôcp

cos а

0ф cos ß

ôcp

cos y.

(5)

dm

дх

 

~w

dz

 

 

Формула (5) представляет выражение производной скалярной функции ф (М) по любому данному направлению m (cos а, cos ß, cos у) через производные этой функции по трем взаимно перпен­ дикулярным направлениям i, j и к.

И р и м е р. Рассмотрим поле температур cp (М), создаваемое беско­ нечно длинным прямолинейным источником тепла, расположенным по оси Oz

(рис. 127).

Пусть, например, ф = (хг + у2)~]. Поверхностями уровня являются цилиндрические поверхности вращения вокруг источника тепла. Производ­ ную функции ф по направлению m (cos а, cos ß, cos у) найдем по формуле (5)

—2 (х2-ту2)~2coset-)- у cos ß).

В

частности,

если

пк к,

то cos а =- cos ß= 0

и - f^ -= 0

(вдоль ис-

 

 

 

 

 

 

r

dm

 

 

точішка тепла температура не меняется).

Если m —грг, у),

 

X

то cos а — — ,

 

У

Зф

2 (ж2 4-у2)

2

0

(чем

дальше

от

источника

cosß = — и

—Т -= ------—

= -------—<

v

г

dm

 

г5

гз

 

 

 

2

тепла,

тем

холоднее).

Если пі = —г(х, у),

то

 

 

получим

= -}■ — (в на­

правлении к источнику тепла температура возрастает).

184. Градиент скалярного поля. Пусть дано скалярное поле Ф (х, у, z). Среди различных возможных направлений в данной фиксированной точке М пространства выделим направление,

перпендикулярное к

поверхности

уровня

скалярной функции

Ф (х, у,

z), проходящей через точку М. Одно из двух таких напра­

влений

представлено

вектором

 

 

 

 

п(ірі, фу,

Ф^)

(6)

{см. п. 133). Вектор п играет важную роль в теории поля и назы­ вается его градиентом.

Ф о р м а л ь н о е о п р е д е л е н и е г р а д и е н т а : гра­ диентом скалярного поля ф (М ) в точке М называется вектор, обо­

значаемый п — grad ф (М ) и

определяемый

равенством

 

grad ф (Ж);

0ф

дф

 

m

дх

дх i + f r

k -

Если функция ф (х,у, z) имеет частные производные первого по­ рядка, то в каждой точке М эта функция имеет свой градиент. Таким образом, скалярное поле ф (М ) порождает векторное поле grad ф (М).

Соседние файлы в папке книги из ГПНТБ