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

книги / Практические занятия по высшей математике. Численное решение алгебраических и трансцендентных уравнений, матричное исчисление, векторный анализ и интегрирование линейных дифференци

.pdf
Скачиваний:
20
Добавлен:
12.11.2023
Размер:
15.42 Mб
Скачать

Координатная поверхность. Координатной называется поверх­ ность, на которой одна из координат точки остается постоянной. Координатными поверхностями в прямоугольной системе координат являются плоскости, параллельные координатным плоскостям хОу,

хОг и уОг (фиг. 17,1).

На каждой

из этих плоскостей одна ко*

ордината сохраняет

постоянное значение.

Координатными

поверхностями

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

ординат являются:

I)

плоскости,

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

т. е. поверхности,

на которых координата г остается постоянной:

2) поверхности прямых круговых цилиндров, общей осью которых является ось Ог; на этих поверхностях постоянное значение имеет

координата г; 3)полуплоскости, проходящие через ось Ог н огра­ ниченные ею; на них координата <р сохраняет постоянное значение (фиг. 17,2).

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

являются: 1) сферы с центром в начале координат; на ннх коорди­ ната р сохраняет постоянное значение; 2) полуплоскости, прохо­ дящие через ось Ог и ограниченные ею; на этих кооодинатных по­ верхностях постоянное значение сохраняет координата <р и 3) кру­

говые конусы, общей

осью которых является ось Ог; на каждой

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

поверхностей координата 0 сохраняет посто­

янное значение (фиг. 17,3).

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

В прямоугольной системе координат такими линиями являются прямые, параллельные координатным осям. Напрнмер, координат­ ная х-линия параллельна оси Ох (фиг. 17,4).

Координатными

линиями в

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

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

являются: г-линня

— прямая,

параллельная оси

Ог, <р-линия — ок-

У

г

Фиг. 17,3

ружность, лежащая в горизонтальной плоскости с центром на оси Ог, г-линия — луч, выходящий из произвольной точки оси Ог, парал­ лельный плоскости хОу (фиг. 17,5).

Координатными линиями в сферической системе координат яв­ ляются: р-линия — луч, выходящий изначала координат, 0-линия — полуокружность с центром в начале координат, соединяющая две

точки на оси Ог, и |р-линия — окружность с центром на оси Ог, лежащая в плоскости, параллельной плоскости хОу (фиг. 17,6).

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

которые

являются

криволинейными ортогональными,

существуют

и другие.

 

 

 

 

 

 

 

 

 

Рассмотрим в общем

виде систему ортогональных

криволиней­

ных координат (фиг. 17,7).

 

 

 

Из формул (17,1) и (17

 

 

 

прямоугольные координаты х

 

 

 

являются

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

 

 

 

динат функциями г,

9

и г ,

 

 

 

 

 

 

Х= х(г,

9);

 

 

 

 

 

 

у =

у</,

9);

 

 

 

 

 

 

г =

г,

 

 

 

 

 

 

а в

сферической

системе к

 

 

 

и г

является

функциями

р,

 

 

 

 

 

 

* =

* ( р

, 0 ,

9 );

 

 

 

 

 

 

у * = у ( р.

в,

9 );

 

 

 

 

 

 

г = г(р. 6).

 

 

 

 

Пусть

в

общем

случае

в трехмерном *

 

 

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

Ф“г- 17,6

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

также

криволинейная

 

 

система координат (и, V, ы>) и между ними и прямоугольными коорди­

натами х, у и г установлено взаимно однозначное

соответствие, ко­

торое описывается формулами, аналогичными формулам (А) и (6 ):

х** х(и,

V,

ш);

 

у = у(и,

V.

ш);

(17,3)

г = г {и, о, ш).

Каждая пара координатных поверхностей, проходящих через фик­ сированную точку М, образует в пересечении координатную линию. В фиксированной точке М проведем касательные к координатным линиям Ми, Мо и Мьи. Будем рассматривать только ортогональ­ ные системы криволинейных координат. Это означает, что векторы

21 г. Ц. лежащие на этих касательных, будут попарно пер­ пендикулярны.

Выражения для

элементов длины, площ ади объема

в ортогональных

криволинейных координатах

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

щими

значениям

координат

и,

V,

ю и и

+

Аи,

V + Ао. и» + А ».

Точка

М имеет

координаты

и,

V,

и>: М (и,

V,

ш).

«Гранями» этого

параллелепипеда являются координатные поверхности: ММ^зМ* ЛШ*М»М,; ЛШ,Л14Л1,;

М7М4УИвУИ4; МзМзМцМз, М2М3УИвМ,.

Его «ребрами» являются координатные линии, которые получаются в пересечении указанных «граней».

Квадрат элемента длины в прямоугольной системе координат

<№= № + йуг + йг\

117,4)

Из формул (17,3) следует, что

*9 =

^

( В Д

Возведем каждое из этих равенств в квадрат, почленно их сложим и подставим в (17,4). Тогда _

-[(к) +(!) +(а!) ] а“*+Н и т

д а

т(я)’]**+К)* +(йУ+(в)’]

+ 2 [к •%+

+ ! •,!+ а! •

^

+ 2 [ !

• ! +

Й ■ Й + !

• ! ] ““ ^ +

+ 2 [й ! +

! •

! + я

• » ]* * ’ *•■

 

Рассмотрим три

вектора Ъх,

Ъг, и

1 „

лежащие на касатель­

ных, проведенных из точки М к координатным линиям ММи ММ% и ЛШ3 соответственно:

 

 

 

 

 

 

 

и

 

д х?

.

ду 7

,

дг-г

 

 

 

 

 

 

 

 

 

 

 

Ъ 'и 1

+

й

 

/

ди

 

 

 

 

 

 

 

 

 

 

 

и

дхт

,

ду- .

дг т

 

 

 

(17,6)

 

 

 

 

 

 

 

Я5‘ + 3 5 / + Ж Л:

 

 

 

 

 

 

 

 

т

 

дх т

,

ду т

,

дг т

 

 

 

 

 

 

 

 

 

 

*

 

001 ^

* дш^

'

За*

*

 

 

 

где

^ ,

=~,

^

— проекции вектора

2-ж, на

оси Ох,

Оу.Ог;

 

дх

ду

дг

 

 

 

 

вектора

^

2 на

те

же

оси;

 

дй'

де »

айпроекции

ь

 

Й ’

Й

 

— лроеишш

вектора

1 3

на

те

же

оси;

а I,

|,

и

4 — орты

этих

осей.

 

 

 

 

 

 

 

 

 

 

Так

как векторы

1 ь

1 2,

и

 

попарно перпендикулярны, то их

скалярные произведения Ъг • 1 2,

7.2 7-3

и

7^ • 7-3 равны нулю. Но

на основании

(17,6)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дх

дх

. ди

ду

 

дг

дг

 

 

 

 

 

 

 

 

 

* ди'

Тк>'

д и ' д и '

ды

Ъй~~

 

*

 

 

 

 

и

и =

дх

 

дх

х ду

ду

 

дг

 

дг __

 

 

 

 

 

 

(&' дм

'

ди ' ди)'

Ъи ' Тки ~

 

 

 

 

 

 

7*2 • 7* =

дх

 

дх

,

ду

ду

 

,

дг

 

дг_

 

,

 

 

 

 

ди ' Ш )'

ди ' д и)'

ди ' ди) = 0

поэтому

в

формуле для

Й52 три

последние

слагаемые равны нулю

и тогда

оказывается,

что в

случае

ортогональных

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

координат квадрат элемента длины содержит только квадраты диф­ ференциалов с1иг, йЬ%н <1ьи, а

(17,7)

Выражения, стоящие здесь в квадратных скобках, есть квадраты

длин векторов Ъи Ъ3 и 13, определенных формулами (17,6), т. е. эти выражения равны

(17.8)

поэтому

еЬ* = Ь\йи* +

На основании формулы (17,8) длины «ребер»

рассматриваемого

бесконечно малого параллелепипеда следует принять , равными:

 

 

(<&)„ =» ММХ=*

 

 

 

 

{<Щ = Ш г ~1.гйи-,

 

(17,9]

 

 

(<к)ш— ММ7 = Ьа йи).

 

 

Этн формулы и выражают элементы длины

в ортогональных

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

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

 

 

Множители

Ц,

Ьг, Ц,_ которые входят в эти

формулы, есть

длины векторов

Ъи

Ь3 и Г3, введенных формулами

(17,6), равные

1' - \ /

вЭ’ +вЭ’ +Й)'*

 

(17,10)

с * =

( Й ) + ( Й ) + ( Й ) *

Величины 7-ь Ц, Ц называются коэффициентами Ляме, отвечаю­ щими точке с координатами и, V и да. От точки к точке коэф­ фициенты Ляме будут изменяться, поэтому они являются функция­ ми координат точки. Вследствие этого коэффициенты Ляме назы­ вают также единицами локальной (местной) длины.

Элемент площади в ортогональных криволинейных координатах. Так как рассматриваемая система криволинейных координат орто­ гональна, то площадь «грани» ММ\М3Мг, которую с точностью до бесконечно малых высшего порядка будем считать прямоуголь-

НИКОМ

Щио - Щи • №)«-

(М«) • (М«);

Щ ио = / * М “ Лг}-

(17,11)

Аналогично площади двух других граней ММгМ4М, и ММХМ4М, равны соответственно

Щиш = Щ о Щ и, = ^«М°

 

(17,12)

Щиш = Щ и *Щ и, = ^

<*Ш.

(17.13)

Формулы (17,11), (17,12) и (17,13) и определяют элемент площади в ортогональных криволинейных координатах. .

Элемент объема в ортогональных криволинейных координатах. Чтобы определить элемент объема в ортогональных криволинейных координатах, найдем объем рассматриваемого бесконечно малого криволинейного параллелепипеда. Считая с точностью до бесконеч­ но малых высшего порядка этот параллелепипед прямоугольным, определим, что его объем

= (<&)„ • (4$)„

Л» (Мы) • (Мы) *(М«>)-

Окончательно элемент объема в ортогональных криволинейных ко­ ординатах

IIV= М г М “ Л».

(17,14)

Выражение проекций вектора А на касательные к координат­ ным линиям Ми, М„ Мш, проведенным в точке М (х, у, г) через его проекции на оси прямоугольной системы координат.

Проекции вектора А на оси прямоугольной системы координат

обозначим

через Ах, Ау, Аг, а на касательные к линиям Ми, Мо и

Мш — через Аи, Аи и Аш.

 

 

 

 

Тогда по известным формулам векторной алгебры имеем

 

Аи = Ахсо5 (х,

и) +

А„соз (у, и) +

Агсоз (г,

ы);

 

А0=

Ахсоз (х,

о) +

А,соз (у,о) +

Агсоз (г,

V) ;

(17,15)

А„ =

Ажсоз (х,

ш) +

А^соз {у, ш)

+ А, соз (г, ш).

 

Входящие сюда девять косинусов легко определятся из формул

(17,6). Косинус угла между вектором а и осью I равен проекции вектора на эту ось, разделенной на модуль вектора» т. е.

соз (а, Ц — — .

Так как косинус угла между

касательной к линии Ми и осью Ох

равен косинусу того угла,

который вектор

составляет с осью

Ох, то, учитывая, что

~ , имеем

 

дх

СОЗ (х>

дх

' Зи*

 

Аналогично получаем и значения остальных косинусов. Таким об­ разом,

С05 (х ,

«) =

Ц

дх

со${у,

=

соз (г,

4 =

ц

Я*.

Т и 9

ди *

С05 (X,

V) —

Т~

дх #

005 (у,

=

соз (г,

о) =

у-

дг

со$ (*,

 

ь2

до 9

 

 

соз (г,

о>) «

г

• Ё

и>) = г

в*. с о в ( у , ® ) = Л . ^ ;

 

 

 

дш9

 

 

 

 

 

дг

 

 

 

 

 

 

 

 

д» *

В дальнейшем рассматриваются только Ортогональные криво­ линейные координаты, коюрые для сокращения записей мы будем называть просто криволинейными, опуская слово «ортогональные».

Вывод формул для вычисления градиента, дивергенции, ротора и оператора Лапласа в криволинейных координатах

Задача

17,1. Найти коэффициенты Ляме Ц,

и Ц для цилин­

дрических

координат.

 

 

 

Р е ш е н и е .

В цилиндрических

координатах

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

ординаты

и, V ,

о) имеют значения

и = г,

о = <р, ш = г, а по фор­

мулам (17,1)

х = гсоз95 {/ =

 

 

 

 

 

Г 51П <р;

= г.

ддг

дх

35"

дх

дш

дх

дг =

дх 0

дх дг =

соз <р;

1 (/5

4

5

0 ;

 

ду

ду

 

дг

дг

ди " д г 3* $1п<р;

 

д г

ду

д у _

гсоз <р;

дг

д г

То

Зу

3? *

5?

ду

ду

0,

дг

дг

дш ~

д г ~

дш

дг

По формулам (17,10) получим

 

 

 

=

Т/со5г <р +

51пг 9 =

1 ;

 

/ -а=

|/ ( г $1Пф)®-|- (г соз ф)4 = т\

(17,16)

* * - 1.

 

 

 

Итак, в цилиндрических

координатах

 

 

~ 1 ,

л,

- 1.

(17,17)

378

 

 

 

 

Задача

17,2.

Найти коэффициенты Ляме для сферических ко­

ординат.

 

 

 

Р е ш е н и е . В

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

динаты и,

V и а) имеют значения

и = р; ц == 0; ш = <р.

По формулам (17,2)

 

 

 

х = р з т

0соз <р;

«/= р з т 6 $!п<р;

г= р С 0 5 6;

дх

д»

 

. .

 

ду __ ду

 

Л п ^ зт # ;

| _

|

соз0;

й =

3?==СО5 <Р5Ш0‘

 

ди

др

 

дх

дх

 

 

.

 

ду __ ду

 

р з т ? соз»;

 

 

— рзш 0;

ЭБ =

=

Р со& Т с о 5 1

 

до

дЬ

 

 

 

 

 

 

 

 

 

 

дх

дх

 

.

.

-

 

 

 

 

дг = — = О

 

«КЛ =

=

— Р51пср31Пв;

2 , = |

=

р С05<р 2*п6;

 

Зй)

ду

 

 

 

7,, =

У (соз у зш О)2 +

(51п з1п 0)* 4- (соз 0)2 =

1;

 

 

 

=

У (р соз <р С05 0)а 4- (р 8,п <Рс08 ®)4 4- (—Р 51п О)2 =

р;

 

 

Ь3шжУ (—р $ т

51П 0)* 4- (р СОЗ 9 ЗШ 0)2 = рЗШ 0.

 

Итак, в сферических

координатах

 

 

 

 

 

 

 

 

 

7-1 =

1;

7,* =

р;

/.э = р з т 0 .

 

(17,18)

Задача 17,3. Определить градиент функции V

V (х, у,

г) в

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

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

причем

 

 

 

 

 

х = х ( и ,

V,

а>)\

у = у ( и ,

V, и>); г =

г {и,

о,

а»).

 

Р е ш е н и е .

Известно, что

 

 

 

 

 

 

По формулам

(17,15) проекция §;га<1иУ на направление касательной

в точке М к

дуге ММ3 (фиг. 17,7)

равна

 

прыгай V = прыгай V •соз (х, и) 4- прыгай V • соз (у,

и) 4 -

 

4- прыгай У • соз (г, и).

■(А)

Учнтывая, что

 

 

 

пР*6гас1 V =

прыгай

прыгай У = ^ ,

 

а на основании формул

(17,16)

 

 

соз (х, и) - ^

соз 0,, и) = ^

соз (г, и) = ^

дг

ди 9

 

 

 

 

переписываем (А) в виде

„„

V=

дУ

1

дх ,

дУ

1

ду

, дУ

1

дг

пр„§га<1

Тх • ^

 

 

 

 

Гц

Тг

^

Ти =

 

 

± ( Ж дх

,<№

ду

,

дУ

дг\

 

 

 

 

^

1.1\дх

 

дц'~ ду ' ди

 

дг

ди)'

 

 

 

Выражение в скобках есть нечто иное, как производная по и от сложной функции V (д:, у, г), где

х = х(и, V, а»); у - у { и , V, ш); г — г {и, V, ш),

дУ Т- е‘Ти-

Поэтому

п р ^ г а с Л '- !

(17,19)

Точно так же найдем проекции градиента функции V на касатель­ ные, проведенные в точке М к дугам ММ2 и ЛШ3:

 

 

 

пр„бга<1 =

ц

^ ;

 

(17,20)

 

 

 

пр„бга(1 =

ц

^ .

 

(17,21)

Поэтому йгаб V в криволинейных

координатах

 

 

 

 

«"“■У = 1% ** + Гл%ёи + ^

е т

(17,22)

где еи, еи, еш— орты

осей Ои, Оо и Ош соответственно.

 

Задача

17.

Найти

выражение

дивергенции

вектора А в

криво­

линейных

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

 

 

 

 

Р е ш е н и е .

Известно, что дивергенцией векторного поля век­

тора А в точке

М называется предел,

к которому стремится

отно­

шение потока этого вектора через замкнутую поверхность 5 , окру­ жающую точку М, к объему V области, ограниченной этой поверх­ ностью, при условии, что объем V стягивается в точку М, т. е.

И^ 9

(11у Л = — ♦ (17,23)

Поэтому мы можем дивергенцию вектора А в точке М вычислить, применяя формулу (17,23) к элементарному объему в криволиней­ ных координатах, т. е. применяя ее к нашему бесконечно малому

криволинейному параллелепипеду. Вектор А в криволинейных коор­

Соседние файлы в папке книги