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

Математический анализ II Учебное Пособие

.pdf
Скачиваний:
87
Добавлен:
14.04.2015
Размер:
7.4 Mб
Скачать

2. Вычисление тройного интеграла

Пусть тело T есть простая область (рис. 10). Допустим, что оно ограничено снизу поверхностью z =z1(x ,y ) , сверху поверхностью z =z 2 (x ,y ) , а с

боков - цилиндрической поверхностью, образующие которой параллельны оси O z , а направляющей служит граница простой области D , расположен-

ной в плоскости xO y , причём функции z1(x ,y )

и z 2 (x ,y ) непрерывны в

области D . Пусть, кроме того, функция f (x ,y ,z )

интегрируема в теле T .

Тогда можно сказать, что

 

 

 

 

 

 

 

z 2

(x ,y )

 

 

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

 

f (x ,y ,z )dz dxdy ,

 

 

1

 

 

 

 

T

D z

 

(x ,y )

 

 

причём интеграл, стоящий справа, записывается так:

∫∫∫f (x ,y ,z )dxdy dz = ∫∫dxdy

z 2 (x ,y )

 

f (x ,y ,z )dz .

T

D

 

 

z1 (x ,y )

 

В том случае, если область D

ограничена

снизу непрерывной кривой

y =y1(x ) , сверху - непрерывной кривой y =y 2 (x ) , а с боков прямыми x =a и x =b , то последнюю формулу можно записать так:

 

b

y 2 (x )

z 2

(x ,y )

∫∫∫f (x ,y ,z )dxdy dz = dx

dy

 

f (x ,y ,z )dz .

T

a

y1 (x )

z1 (x ,y )

Интеграл, стоящий справа, называется трехкратным или повторным. Заметим, что выбирая внешнее интегрирование по переменной y или z , можно написать ещё пять различных трехкратных интегралов, через которые выражается данный интеграл I . Порядок выполнения операций интегрирования зависит от вида области, по которой выполняется интегрирование.

3.Геометрический смысл тройного интеграла

Втом случае, если подынтегральная функция f (x ,y ,z ) 1, то очевидно,

что тройной интеграл I = ∫∫∫dxdydz даёт нам объём тела T , по которому

T

ведётся интегрирование.

Пример 1. Вычислить объём тела, ограниченного координатными плоскостями и плоскостью x +y +z 1 = 0 (P ) (рис. 11).

Решение. Искомый объём v = ∫∫∫dxdydz , где тело T есть пирамида,

T

ограниченная плоскостью P и координатными плоскостями.

Очевидно, что этот тройной интеграл можно выразить шестью различными способами через трёхкратный:

64

 

 

z

 

 

1xy

1

1x

1xy

1

1y

1xy

 

 

 

 

v =∫∫dxdy

dz =dx dy

dz =dy dx dz

 

1

 

 

Dxy

0

0

0

0

0

0

0

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1yz

1

1y

1yz

1

1z

1yz

 

D xz

D y z

y

v =∫∫dydz

dx =dy dz

dx =dz dy

dx

 

Dyz

0

0

0

0

0

0

0

 

D xy

1

или

 

 

1x 1xz

 

1z 1xz

x

 

 

 

1xz

1

1

1

 

 

v =∫∫dxdz

dy =dx dz

dy =dz dx dy

 

рис. 11

 

 

 

 

Dxz

0

0

0

0

0

0

0

 

 

 

 

Проведём вычисления по последней формуле (нетрудно убедиться, что остальные интегралы вычисляются аналогично), получим

v = 1 dz 1z dx 1xz dy .

0

0

0

 

 

 

 

 

 

 

 

 

 

 

 

Имеем

 

 

 

 

 

 

 

 

 

 

1z

 

 

 

1x z

 

 

 

 

1z

 

 

2

 

 

 

1 z

2

 

 

 

 

 

 

 

I 1Вн = dy =1 x z ; I 2Вн =

 

 

 

 

=

) .

(1 x z )dx = x x

 

xz

 

 

(

0

 

 

 

 

0

 

2

 

 

0

 

2

 

 

 

 

 

 

 

 

 

 

 

Наконец, v = 1

(1 z )2

dz =

1

куб. ед.

 

 

 

 

 

 

 

 

2

6

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

§3. Приложения двойных и тройных интегралов

1.Вычисление площади плоской области D

Мы установили, что площадь плоской области может быть вычислена по формуле

SD = ∫∫dxdy .

D

2. Вычисление площади кривой поверхности

Ранее мы установили, что площадь кривой поверхности S , заданной уравнением z = f (x ,y ) и расположенной над областью D в плоскости

xO y , вычисляется по формуле

 

S = ∫∫

p 2 (x ,y ) +q 2 (x ,y ) +1dxdy ,

 

z (x ,y )

D

 

z (x ,y )

 

где p (x ,y ) =

, q(x ,y ) =

.

x

 

 

 

 

y

65

3. Вычисление объёма телаT

Пусть тело T ограничено снизу простой областью D , сверху - поверхностью z = f (x ,y ), f (x ,y ) 0, f (x ,y ) непрерывна в области D , а с боков

цилиндрической поверхности, образующие которой параллельны оси O z , а направляющей служит контур области D , то объём тела можно вычислить с помощью двойного интеграла так:

vT = ∫∫f (x ,y )dxdy

D

или так:

vT = ∫∫∫dxdy dz .

T

4. Вычисление массы поверхности

Пусть в каждой точке поверхности S , заданной уравнением z = f (x ,y ) плотность равна ρ(x ,y ,z ) , где ρ(x ,y ,z ) - непрерывная функция в каждой

точке поверхности S , а функция f (x ,y ) - непрерывна в области D

плоско-

сти xO y и имеет в ней непрерывные частные производные

 

 

f (x ,y )

= p (x ,y ) и

f (x ,y )

=q(x ,y ) .

 

 

 

y

 

 

x

 

 

 

Разбивая произвольным образом поверхность S на n частей,

заметим,

что масса k -й ячейки приблизительно равна

m k ρ(xk ,y k ,z k )

Sk , где

Sk - площадь k -й ячейки, тогда масса всей поверхности S :

 

 

n

n

 

 

 

M = m k ρ(xk ,y k ,zk ) Sk .

 

 

k =1

k =1

 

 

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

M = ∫∫ρ(x ,y ,z ) p 2 (x ,y ) +q 2 (x ,y ) +1dxdy .

D

Вчастности, если поверхность S лежит в плоскости xO y , т.е. совпадает

собластью D , то M D = ∫∫ρ(x ,y )dxdy .

D

5. Вычисление массы тела

Если в каждой точке тела T , ограниченного простой поверхностью, задана плотность ρ = ρ(x ,y ,z ) , где ρ - непрерывная функция в каждой точке

телаT , то, проведя аналогичные рассуждения, получим, что масса тела

66

M T = ∫∫∫ρ(x ,y ,z )dxdy dz .

T

6. Моменты плоской фигуры

Из курсов теоретической механики известно, что статическим моментом Sl материальной точки массы m относительно оси l называется произве-

дение массы этой точки на расстояние до оси l , т.е. Sl =m rl .

Моментом инерции I l материальной точки относительно оси l называется произведение массы m этой точки на квадрат её расстояния от оси l ,

т.е. Sl =m rl 2 .

Статическим моментом системы материальных точек относительно оси l называется сумма статических моментов относительно этой оси всех материальных точек, входящих в систему.

Пусть в плоской области D распределена масса с плотностью ρ(x ,y ) . Разобьём область D на n частей, Sk - площадь k -й ячейки (k =1, 2,...,n ). В ячейке Dk возьмём произвольную точку (xk ,y k ) , тогда в силу сделанного выше определения можем считать, что

n

n

Sx ρ(xk ,y k ) y k

Sk ; Sy ρ(xk ,y k ) xk Sk .

k =1

k =1

Измельчая дробление, в пределе получим точные значения для статических моментов области D о осей O x и O y :

Sx = ∫∫ρ(x ,y )y dxdy ;

Sy = ∫∫ρ(x ,y )x dxdy .

D

D

Проводя аналогичные рассуждения для моментов инерции области отно-

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

 

I x = ∫∫ρ(x ,y )y 2 dxdy ;

I y = ∫∫ρ(x ,y )x 2 dxdy .

D

D

7. Координаты центра масс

Пусть D - плоская область, в которой распределена масса с плотностью ρ(x ,y ) . По определению центром масс плоской области называется точка

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

xc = MSy , y c = MSx ,

где M - масса плоской области D , а Sx и Sy - статические моменты.

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

67

r
ρ

x

 

=

∫∫ρ(x ,y )xdxdy

, y

 

=

∫∫ρ(x ,y )y dxdy

c

D

c

D

 

∫∫ρ(x ,y )dxdy

∫∫ρ(x ,y )dxdy

 

 

 

 

 

 

 

D

 

 

 

D

Рассмотрим теперь некоторое тело T , ограниченное простой поверхно-

стью, и пусть в нём распределена масса,

плотность которой ρ(x ,y ,z ) , то

для координат центра масс этого тела получим совершенно аналогичные выражения:

x

 

=

∫∫ρ(x ,y ,z )xdv

, y

 

=

∫∫ρ(x ,y ,z )y dv

, z

 

=

∫∫ρ(x ,y ,z )zdv

.

c

T

c

T

c

T

∫∫ρ(x ,y ,z )dv

∫∫ρ(x ,y ,z )dv

∫∫ρ(x ,y ,z )dv

 

 

 

 

 

 

 

 

 

 

T

 

 

 

T

 

 

 

T

 

§4. Криволинейные координаты и замена переменных

вкратных интегралах

1.Криволинейные координаты

Относительно пространственной декартовой системы координат O xyz положение точки M определяется её тремя декартовыми координатами x , y и z . Если z > 0 , то положение точки M можно определить, задав такие три параметра: ρ - расстояние

z

M

0

y

точки M

от оси O z ; r

-расстояние точки

x

φ

 

M от начала координат;

ϕ - двухгранный

рис. 12

угол между плоскостью xO z и плоскостью,

 

 

 

 

проходящей через ось O z и точку M (рис.

12). Очевидно,

что параметры

ρ , r и ϕ можно выразить через декартовы координаты x , y

и z . Парамет-

ры ρ , r

и ϕ мы можем также назвать координатами точки M , ϕ [0, 2π[ .

И вообще, за координаты точки M мы можем принять любые три функции:

ξ =ξ(x ,y ,z ), η =η(x ,y ,z ), ζ =ζ (x ,y ,z )

 

(1)

лишь бы только соотношениями (1) координаты x ,y и z определялись од-

нозначно:

 

x =x (ξ,η,ζ ), y =z (ξ,η,ζ ), z =z (ξ,η,ζ )

(2)

Т.е. ни одно из соотношений (1) или (2) не должно противоречить другим или быть следствием других. Заметим, что из соотношений (2) в этом случае параметры ξ, η и ζ также будут определяться однозначно. Можно

доказать, что эти условия выполняются, если определить I (ξ,η,ζ ) , назы-

68

ваемый определителем Якоби или якобианом преобразования, отличен от нуля, т.е.

 

xξ

xη

xζ

 

 

 

 

 

 

I (ξ,η,ζ ) =

y ξ

yη

y ζ

 

0

(3)

 

z ξ

zη

z ζ

 

 

 

2. Координатные поверхности

Зафиксируем какую-нибудь координату, определённую соотношениями (1), например ξ , положив ξ =c1 , тогда получим ξ (x ,y ,z )=c1 (рис. 13).

С геометрической точки зрения этому уравнению в пространстве соответствует некоторая поверхность ξ . Аналогично можно определить коорди-

натные поверхности η и ζ соответственно:

η(x ,y ,z ) =c 2 и ζ (x ,y ,z )) =c3 .

Координатные поверхности ξ =c1 , η =c 2 и ζ =c3 при соблюдении усло-

вия (2) пересекаются в некоторой точке M . Таким образом, точка M определяется как точка пересечения координатных поверхностей (рис. 13).

3. Координатные линии

 

z

к.л.ζ

 

ξ

 

 

η =c 2

 

ξ0

ζ

 

 

M

 

 

r

ζ0

ξ =c1

к.л.ξ

η0

y

 

 

ζ =c 3

0

η

к.л.η

 

 

 

 

x

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

 

(

x ,y ,z

)

=c

1

 

ξ

 

 

 

 

η(x ,y ,z )

=c 2

.

 

 

 

 

 

 

 

 

Очевидно, что кривая, по которой пересекаются эти поверхности, обладает таким свойством, что вдоль этой кривой координаты ξ и η постоянны, а меняется

одна только координата ζ , поэто-

му эта кривая называется координатной линией ζ . Аналогично пе-

рис. 13 ресечение поверхностей

 

(

x ,y ,z

)

=c

1

 

 

 

(

x ,y ,z

)

=c

2

 

ξ

 

 

 

 

и

η

 

 

 

 

ζ (x ,y ,z )=c3

 

ζ (x ,y ,z )=c3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

даёт нам соответственно координатные линии η и ξ .

69

Очевидно, что в общем случае координатные линии представляют собой некоторые кривые, поэтому координаты ξ, η и ζ называются криволиней-

ными координатами. Проведем к координатным линиям, пересекающимся в точке M , касательные, направления которых соответствуют направлениям возрастания координат. Орты этих осей называются ортами криволинейных координатных осей и обозначаются соответственно ξ, η и ζ . Систему кри-

волинейных координат называют ортогональной, если ортогональны орты ξ0 , η0 и ζ0 , т.е. если выполняются условия

ξ0 η0 = η0 ξ0 = ζ0 ξ0 = 0 .

Заметим, что в декартовой системе координат O xy z координатными поверхностями будут являться плоскости, параллельные координатным плоскостям xO y , xO z , y O z , а координатными линиями - прямые, параллельные осям O x , O y и O z .

4. Коэффициенты Ламе

Рассмотрим вектор r - радиус-вектор точки M (рис. 13). Очевидно, что производная ξr лежит на касательной к годографу вектора r , построенно-

му в предположении, что меняется только координата ξ , а координаты η и ζ остаются неизменными в сторону возрастания координаты ξ . Очевидно,

что ξ0 =

r

/

 

r

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ξ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ξ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

r

 

r

 

 

 

 

 

 

0

 

 

 

 

 

 

r

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Аналогично η

=

 

 

 

/

 

 

 

 

 

,

 

 

ζ

 

=

 

 

 

 

 

/

 

 

 

 

 

 

.

 

 

 

 

 

 

 

η

 

η

 

 

 

 

 

ζ

 

ζ

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

r

 

 

 

 

Выражения

H ξ =

 

 

 

, H η

=

 

 

 

, H ζ

 

=

 

 

 

называются коэффициентами

ξ

 

 

 

η

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ζ

 

 

 

 

 

Ламе. Так как радиус-вектор точки M .

 

 

 

r =x i +y j +z k , то, выполняя диф-

ференцирование, получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

x i

 

 

y

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

+

j +

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ξ

 

 

ξ

ξ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ξ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

x

 

 

 

y

 

 

 

z

 

 

(4)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

i

+

 

 

 

 

 

 

 

j +

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

η

 

η

η

η

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

x

 

 

 

y

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

i

+

 

 

 

 

 

 

 

j +

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ζ

ζ

ζ

ζ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Принимая во внимание эти соотношения, нетрудно получить такие вы-

ражения для коэффициентов Ламе:

 

H ξ = xξ2 +′y ξ2 +zξ2 , H η = xη2 +′yη2 +zη2 , H ζ = xζ2 +′y ζ2 +z ζ2

(5)

70

 

Примем во внимание соотношения (4) и (5), тогда выражение для ортов криволинейных координат запишутся так:

ξ0 = (xξi +y ξj +z ξk )/ xξ2 +y ξ2 +z ξ2 η0 = (xηi +yηj +zηk )/ xη2 +yη2 +zη2 ζ0 = (xζi +y ζj +z ζk )/ xζ2 +y ζ2 +z ζ2

5. Сферические координаты

Сферическими координатами точки M называются параметры ρ , θ и ϕ

(рис. 14), где

z

ρ

- расстояние от точки M до начала ко-

 

z

θ ρr M (x ,y ,z )

 

 

 

ординат (0 ρ < +∞);

 

 

 

θ - угол, отсчитываемый от оси O z до век-

 

0

y y

 

тора ρ (0 θ π) ;

 

ϕ

- угол между плоскостью xO z и плос-

 

x

M

x

ϕ

 

костью, проходящей через точку M и

 

 

 

ось O z ( 0 ϕ 2π ).

 

 

рис. 14

 

 

 

 

 

Очевидно, что в сферических координатах координатными поверхностями являются также поверхности : - ρ =c1 - сфера радиуса c1 с

центром в начале координат; θ =c 2 - конус с вершиной в начале координат,

образующие которого составляют

угол

θ =c 2 с осью

O z ;

ϕ =c3

-

плоскость, проходящая через ось

O z

и образующая

угол

ϕ =c 2

с

координатной плоскостью xO z .

 

 

 

 

 

Сфера ρ =c1 и конус θ =c 2 пересекаются по окружности, которая представляет собой координатную линию ϕ . Сфера ρ =c1 и плоскость ϕ =x3 пересекаются по окружности (координатная линия θ ). Конус θ =c 2 и плоскость ϕ =c3 пересекаются по прямой (координатная линия ρ ). Орты криволинейной сферической системы координат изображены на рис. 14.

Непосредственно из рис. 14 можно установить связь между декартовы-

ми и сферическими координатами.

x = ρsinθ cosϕy = ρsinθ sinϕ .

z = ρcosθ

71

Вычислим теперь Якобиан преобразования для сферических координат

sinθ cosϕ J (ρ,θ,ϕ) = sinθ sinϕ

cosθ

Найдём теперь r , r , r

ρ θ ϕ

с соотношениями (4):

ρcosθ cosϕ

ρsinθ sinϕ

 

= ρ2 sinθ .

 

ρcosθ sinϕ

ρsinθ cosϕ

 

ρsinθ

0

 

 

для сферических координат в соответствии

r

 

 

= sinθ cosϕi +sinθ sinϕj +cosθk

 

ρ

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

= ρcosθ cosϕi + ρcosθ sinϕj ρsinθk

θ

 

 

 

 

 

r

 

 

= −ρsinθ sinϕi + ρsinθ cosϕj

 

 

 

 

 

ϕ

 

 

 

 

 

 

 

 

 

 

А теперь нетрудно написать выражения для коэффициентов Ламе в сферических координатах: H ρ =1, H θ = ρ, H ϕ = ρsinθ . Очевидно, что орты

ρ0 , θ0 и φ0 запишутся так:

 

0

 

 

r

 

 

ρ

 

 

=

 

 

 

/ H ρ = sinθ cosϕi +sinθ sinϕj +cosθk

 

 

ρ

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

0

 

 

 

 

θ

 

 

=

 

 

 

/ H θ = cosθ cosϕi +cosθ sinϕj sinθk

 

 

 

θ

 

 

 

 

 

 

 

φ

0

=

 

r

/ H ϕ = −sinϕi +cosϕj

 

 

 

 

 

 

 

 

ϕ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Нетрудно убедиться, что криволинейная сферическая система координат ортогональна. Действительно

ρ0θ0 = ρ0φ0 = θ0φ0 = 0

6. Цилиндрические и полярные координаты

 

Цилиндрическими

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

 

 

 

 

z

 

z0

 

 

точки M

(рис.

15) называются

пара-

 

 

 

 

ρ

 

 

 

метры:

 

 

 

 

 

 

 

 

 

 

 

φ

0

 

точки M от

оси O z

 

 

 

 

 

 

 

 

ρ

- расстояние

 

 

 

0

 

 

ρ0

 

y

 

 

 

 

 

 

 

 

( 0 ρ < +∞);

 

 

 

 

 

 

 

 

 

 

ϕ

- угол, отсчитываемый от плоскости

 

 

 

 

 

 

 

 

 

 

 

xO z

до плоскости,

проходящей через

x

 

 

φ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

точку M и осью O z (0 ϕ < 2π) ;

 

 

 

 

рис. 15

 

 

 

 

 

 

 

72

 

 

 

 

 

 

 

 

 

 

z - координата, совпадающая с декартовой координатой z (−∞<z <+∞).

Очевидно, что в цилиндрических координатах координатными поверхностями являются:

ρ =c1 - цилиндрическая поверхность, образующие которой параллельны оси O z , а направляющей служит окружность в плоскости xO y с центром в начале координат радиуса ρ =c1 ;

ϕ =c 2 - плоскость, проходящая через ось O z и образующая с плоскостью xO y угол ϕ ;

z =c3 - плоскость, параллельная координатной плоскости xO y .

Координатные поверхности ρ =c1 и ϕ =c 2 пересекаются по прямой, па-

раллельной оси O z (координатная линия z ). Координатные поверхности ρ =c1 и z =c 2 пересекаются по окружности (координатная линия ϕ ). Ко-

ординатные поверхности ϕ =c 2 и z =c3 пересекаются по прямой (координатная линия ρ ). Орты цилиндрической системы координат и координатные линии ρ , ϕ и z изображены на рис. 15.

Непосредственно из рис. 15 нетрудно установить связь между декарто-

выми и цилиндрическими координатами:

x = ρcosϕ

 

 

y = ρsinϕ .

z =z

 

 

Найдём Якобиан преобразования для цилиндрических координат:

 

cosϕ

ρsinϕ

0

 

 

 

 

 

 

J (ρ,ϕ,z ) =

sinϕ

ρcosϕ

 

0

 

 

= ρ .

 

0

 

0

 

 

1

 

 

 

Нетрудно получить выражения для

r ,

rr

 

и

r

, а также коэффициен-

ϕ

z

 

 

ρ

 

 

ты Ламе H ρ , H ϕ и H z :

 

 

 

 

 

r

= cosϕi +sinϕj,

r

= −ρsinϕi + ρcosϕj,

r

=k

 

ρ

ϕ

z

 

 

 

 

H ρ =1, H ϕ = ρ, H z =1.

С учётом этих соотношений получим:

ρ0 = cosϕi + sinϕj, φ0 = −sinϕi +cosϕj, . z0 = k

73