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

telnov-machanika-and-TO

.pdf
Скачиваний:
57
Добавлен:
30.05.2015
Размер:
2.14 Mб
Скачать

Границы области движения по радиусу соответствуют r = 0 в (69.14), т.е. r находится из уравнения

E =U(r) +

L2

.

(69.19)

2mr2

 

 

 

Если имеется две границы rmin и rmax , то движение называется финитным (ограниченным). Угол поворота при движении от rmin до rmax и обратно

rmax

 

L

dr

 

 

 

 

 

mr2

 

 

 

 

Dj = 2m ò

 

 

 

 

.

(69.20)

 

 

 

2

rmin

E -U(r) -

L

 

 

2mr2

 

 

 

 

 

 

Траектория замкнута, если Dj = 2p nk , где k,n - целые числа (траек-

тория замыкается после n колебаний по радиусу и k полных оборотов по j). В общем случае траектория незамкнута. Оказывается, она замк-

нута только для

потенциалов

 

U µ -

1

(кулоновское

поле) и

 

 

 

 

 

 

 

 

 

r

 

 

U µ r2 (пространственный осциллятор).

 

 

 

Падение на центр возможно, если r2(r = 0) > 0 , т.е.

 

E -U(r) -

L2

 

> 0

 

при r 0 ,

(69.21)

2mr2

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

L2

 

2

= 0 (при r = 0).

(69.22)

r U(r) +

 

< Er

 

2m

 

Поскольку E – конечная величина, но последний член в выражении равен нулю. Отсюда следуют условия падения на центр

1)

L = 0 ;

 

 

(следует из (69.21))

 

(69.23)

2)

U = -

a

 

 

и a >

L2

; ( a =

L2

- круговая орбита)

(69.24)

r2

 

 

 

2m

 

 

 

 

 

2m

 

 

3)

U = -

a

 

,

n > 2, a > 0.

 

 

(69.25)

rn

 

 

 

 

 

 

 

 

 

 

 

 

 

.

171

Uэф(r)

В слабом гравитационном поле U µ -r1 , поэтому частицы не па-

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

U µ -ar - rb2 ..., что приводит к падению материи на центр тяготения.

Примеры эффективного потенциала

1) Притягивающий кулоновский (гравитационный) потенциал

U = -ar ,

Uэф(r)

E2

r

E1

rmin rmax

Рис. 49

U

 

= -a

+

 

L2

.

(69.26)

эф

 

2mr2

 

r

 

 

 

 

На

малых

расстояниях

доминирует

второй член, на больших – первый член. При E > 0 (линия E2 ) частица,

прилетевшая из бесконечности, отразится от потенциального барьера и улетит обратно (инфинитное движение). Для финитного движения необходимо

E < 0 (линия E1 ), при этом движение по радиусу ограничено rmin и rmax .

2) Пространственный осциллятор

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

E

U = kr2

, U

 

= kr2

+

L2

. (69.27)

 

эф

2mr2

 

2

 

2

 

 

Движение финитно при любой физически возможной энергии.

rmin

rmin rmax

Рис. 50

172

Движение по круговым орбитам

Частица движется по окружности, если

rmin = rmax . При этом

¢

. Пусть

 

 

 

 

 

 

Uэф = 0

 

 

 

 

 

 

 

U = -

a

,

F = -

 

na

.

(69.28)

 

 

 

 

 

 

 

rn

 

 

rn+1

 

Частица будет двигаться по окружности при

 

 

 

mv2

=

na

.

(69.29)

 

 

 

r

 

 

 

 

 

rn+1

 

Это получается как из приравнивания ускорения вызванного силой ус-

корению, необходимому для движения по окружности (r = -v2

), а

 

r

 

¢

, где Uэф дается (69.26).

 

так же изUэф = 0

 

Казалось бы, все просто, при каждом радиусе найдется скорость, необходимая для кругового движения. Но будет ли это движение устойчиво? Движение устойчиво, когда на круговой орбите эффективный потенциал имеет минимум, и неустойчиво (падение на центр) когда в

этой точке Uэф имеет максимум.

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

a > 0,

0 < n < 2 ;

(69.30)

a < 0,

n < 0.

(69.31)

При устойчивом движении небольшой радиальный толчок тела приведет к колебаниям вблизи исходной круговой орбиты, например, при n = 1 окружность превращается в эллипс.

Гравитационный потенциал удовлетворяет условию (69.30), а пространственный осциллятор условию (69.31). Однако, если потенциал

U(r) = -ra10 , то тело будет двигаться по окружности, но неустойчиво.

Малейшее начальное радиальное движение будет нарастать, и тело упадет на центр или улетит на бесконечность.

173

§ 70. Кеплерова задача

Немецкий математик и астроном Кеплер (1571-1630) из анализа данных датского астронома Тихо Браге (1546-1601) установил три эмпирических закона:

1)Каждая планета Солнечной системы обращается по эллипсу, в одном из фокусов которого находится Солнце.

2)Каждая планета движется в плоскости, проходящей через центр Солнца, причём за равные времена радиус-вектор, соединяющий Солнце и планету, описывает равные площади.

3)Квадраты периодов обращения планет вокруг Солнца относятся, как

кубы больших полуосей орбит планет:

T 2

 

a3

 

1

=

 

1

.

T 2

 

 

 

a3

 

2

 

2

 

Задача о круговом движении в гравитационном поле довольно простая и еще до Ньютона было известно, что, если F µ -r12 , то получа-

ются законы Кеплера. Исаак Ньютон (1642-1727) математически показал, что, действительно, при F µ -r12 планеты движутся по эллипсам

и при этом выполняется третий закон Кеплера. Рассмотрим эту задачу.

Траектория движения

Итак, имеем

U = -a

,

U

 

= -a

+

L2

, a > 0,

эф

2mr2

 

 

r

 

 

r

 

 

(70.1)

U

 

= -

ma2

при r =

L2

 

,

эф, min

 

 

 

 

 

 

2L

 

*

 

ma

 

 

 

 

 

 

 

 

 

где последнее соотношение соответствует круговому движению. Из (69.18) получаем

174

 

 

 

L

dr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j = ò

 

 

r2

 

 

 

 

 

 

 

 

 

+const =

 

æ

 

2ma

 

 

 

 

 

2

ö

 

 

 

ç

 

 

 

 

L

÷

 

 

 

 

 

 

ç2mE +

 

 

 

 

-

 

 

 

2

÷

 

 

 

 

 

 

 

 

 

 

 

 

 

÷

 

 

 

 

 

 

ç

 

 

 

r

 

 

 

r

 

÷

 

 

 

 

 

 

è

 

 

 

 

 

 

 

ø

 

 

 

 

 

 

 

 

 

 

 

æ

 

 

 

am

ö

 

 

 

 

 

 

 

 

 

çL

-

÷

 

 

 

 

 

 

 

 

d ç

 

 

 

 

 

 

÷

 

 

 

 

 

 

 

 

 

ç

 

 

 

 

 

L

÷

 

 

 

= -ò

 

 

 

 

èr

 

 

 

 

 

ø

 

 

+const.

æ

 

 

 

a2m2

÷ö

 

 

æL

 

am ö2

 

ç

 

+

 

2

÷

 

 

ç

 

-

÷

 

 

 

ç2mE

 

÷

 

÷

 

 

ç

 

 

 

L

 

÷

 

 

ç

 

 

÷

 

 

è

 

 

 

 

ø

 

 

èr

 

L ø

 

Учитывая, что ò

 

dx

 

= -arccos x

, получаем

 

2

2

 

a

-x

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

L

-

am

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j = arccos

r

L

 

 

 

 

+const.

 

 

 

 

 

 

 

2mE +

a2m

2

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

(70.2)

(70.3)

Выбором начала отсчета делаем const = 0 . Отсюда траектория движения

 

p

= 1 +e cosj ,

(70.4)

 

r

 

 

 

 

где введены обозначения

 

 

 

 

 

p =

L2

, e =

1 +

2EL2 .

(70.5)

ma

 

 

 

ma2

 

Это кривые, соответствующие коническим сечениям (кривые, по которым поверхность круглого конуса пересекается плоскостью, см. курс аналитической геометрии). Величины p и e называются

параметром и эксцентриситетом орбиты.

В аналитической геометрии показывается, что при e = 0 – это круг;

e < 1 – эллипс (E < 0 -финитная траектория); e = 1– парабола (E = 0 );

175

e > 1 – гипербола (E > 0 ).

Рассмотрим случай эллиптической траектории. Как известно, у эллипса имеется два фокуса, и сумма расстояний от фокусов до любой

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

точки постоянна. Уравне-

 

 

 

 

p

 

 

 

 

 

 

 

 

 

ние эллипса в декартовых

b

 

a

 

 

 

 

 

 

 

 

 

координатах

 

 

 

 

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2b

 

 

x2

y2

= 1 .

(70.6)

O¢

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ b2

C c

O

 

 

 

 

 

 

 

 

 

 

 

a2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Центр тяготения находит-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ся в одном из фокусов,

 

2a

 

 

rmin

 

 

 

 

 

 

пусть в правом, рис. 51.

r

 

 

 

 

 

 

 

 

 

 

 

Из

 

(70.4)

 

минимальное

max

 

 

 

 

 

 

 

 

 

 

 

 

 

 

расстояние от фокуса (пе-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 51

 

 

 

 

 

 

 

 

 

 

 

ригелий) и

 

максимальное

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

расстояние (афелий)

 

 

 

 

r

 

=

 

 

p

 

, r

=

 

p

 

;

 

 

 

(70.7)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

min

 

1 +e

 

 

max

 

1 -e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

параметр p равен радиусу при j

= p

. Большая полуось

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

a =

1

(r

+r

 

 

)

=

 

p

 

 

= -

a

 

=

a

 

,

(70.8)

2

 

 

1 -e2

2E

 

2 | E |

 

min

 

max

 

 

 

 

 

 

 

 

(при финитном движении энергия отрицательна, поэтому здесь использован модуль E , чтобы знак минус не вводил в заблуждение),

отсюда

rmin = a(1 -e),

rmax

= a(1 +e).

(70.9)

Расстояние между фокусом и центром эллипса

 

 

 

c = 0.5(rmax -rmin ) = ae;

 

 

(70.10)

малая полуось

 

 

 

 

 

 

 

b = a2 -c2 = a 1 -e2

=

 

p

=

L

.

(70.11)

 

 

 

 

 

 

1 -e2

2m | E |

 

176

Как мы видим, большая полуось a зависит только от полной энергии, а малая полуось зависит также от момента импульса.

Заметим, что выражение для rmax,rmin и a легко получить из законов

сохранения энергии и момента импульса. Действительно, в перигелии и афелии нет радиального движения, тогда

L = mvr,

 

E =

mv2

 

a

,

 

 

(70.12)

 

2

-

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

отсюда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

æ

 

 

 

 

 

 

2

÷ö

 

 

 

 

 

a ç

 

 

 

2EL

÷

 

 

r

= -

 

 

ç1

 

1

+

 

 

 

2

÷

,

(70.13)

 

 

 

 

 

max,min

 

 

 

ç

 

 

 

 

 

 

÷

 

 

 

 

 

ç

 

 

 

ma

 

÷

 

 

 

 

 

2E ç

 

 

 

 

 

 

 

 

 

 

è

 

 

 

 

 

 

 

ø

 

 

a =

1

(r

 

+r

 

) = -

a

,

 

 

(70.14)

2

 

 

2E

 

 

 

min

 

max

 

 

 

 

 

 

 

что совпадает с (70.8) и (70.9).

Период обращения

Из сохранения момента импульса

2

dS

 

 

L

(70.15)

L = mr j

= 2m dt

,

S =

 

T,

2m

где T – период обращения. Учитывая, что для эллипса S = pab получаем

T =

2pmab

= pa

 

m

 

= 2pa3/2

m .

(70.16)

L

 

2 | E |3

 

 

 

 

a

 

Отсюда следует третий закон Кеплера

 

 

 

 

 

T 2 =

4p2m a3 .

 

 

(70.17)

 

 

 

 

a

 

 

 

Этот закон легко получить в случае кругового движения, но оказывается формула остается справедливой и для эллиптического движения, если вместо радиуса в формуле использовать большую полуось эллипса. Для гравитационного взаимодействия тел с массами M и m подставляем

177

 

a =GMm,

 

m m =

 

Mm

,

(70.18)

 

M +m

 

 

 

 

 

 

 

 

 

 

 

 

и получаем

T 2 =

 

 

4p2a3

 

.

 

 

(70.19)

G(M +m)

 

 

 

 

 

 

 

 

 

Случай e = 1,

E = 0 -парабола

 

 

 

 

 

 

 

 

 

Dj

Из (70.4) следует, что в этом случае

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

r

= p

=

 

L2

 

 

(70.20)

 

 

2ma

 

 

min

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и угол поворота при облете центра тяготения

 

Dj = 2p, т.е. траектория

поворачивает

в об-

 

ратную сторону (cos 1

при r ). Выра-

 

жение (70.20) также легко получить из законов

Рис. 52

сохранения энергии и момента импульса, пола-

гая в них E = 0 , тогда rmin

находится из

 

 

 

 

 

 

L = mvr,

E

=

mv2

-

a

= 0.

(70.21)

 

 

2

r

 

 

 

 

 

 

 

 

 

 

 

 

 

Случай e > 1, E > 0 , — гипербола

Гипербола – геометрическое место точек, для которых абсолютная величина разности расстояний до двух заданных точек (фокусов) (-c, 0),

(c, 0)постоянна и равна 2a < 2c . Уравнение гиперболы в декартовых координатах

x2

- y2

= 1.

(70.22)

a2

 

b2

 

 

 

В полярных координатах (70.4)

 

p/r = 1 +e cosj,

e > 1 . Минималь-

ное расстояние равно

 

 

 

 

 

 

 

r

 

=

 

 

p

.

(70.23)

 

 

 

 

min

 

e +1

 

 

 

 

 

Полагая r = ¥, находим cosj = -

1 . Половина угла поворота траек-

 

 

 

 

 

e

 

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

178

 

 

 

 

p

 

æ

 

 

ö

 

p

 

 

Dj = j

-

 

ç

 

1÷

-

 

(70.24)

 

= arccosç-

 

÷

 

 

 

 

 

 

2

 

ç

 

e

÷

 

2

 

 

 

 

 

 

 

è

 

ø

 

 

 

отсюда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

æ

 

ö

 

 

 

1

 

 

 

 

 

1

 

ç

p÷

 

 

 

 

 

 

 

 

(70.25)

-cosçDj +

 

÷

= sin Dj =

 

=

 

 

 

 

 

 

 

 

 

 

 

 

ç

2

÷

 

 

 

e

 

 

 

 

 

 

2

 

è

ø

 

 

 

 

 

 

 

1 +

2EL

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ma2

 

Если частица имела на

бесконечности

скорость

v¥ и прицельный

параметр (расстояние между продолжением исходной траектории и центром тяготения) r, то подставляя в (70.25)

L = mrv , E =

mv2

,

¥

2

¥

 

 

 

нетрудно получить

tg Dj = 2raE = GMrv2 .

¥

(70.26)

(70.27)

 

Проследим эволюцию орбит при измерении

 

параметров. Пусть тела пролетают на одном

 

и том же расстоянии от центра тяготения, но

 

с разными скоростями, рис. 53. При малой

 

скорости центр тяготения находится в пра-

 

вом фокусе эллиптической орбиты. При уве-

 

личении скорости орбита становится круго-

 

вой. При еще большей скорости орбита ста-

 

новится снова эллиптической, но центр тяго-

 

тения находится уже в левом фокусе эллип-

 

са. При еще большей скорости траектория

Рис. 53

становится гиперболической.

Случай отталкивающего потенциала

Пусть

U = a

, U

 

= a

+

L2

,

a>0 .

(70.28)

эф

2mr2

r

 

r

 

 

 

 

После вычислений аналогичных (70.2) получается уравнение траектории, отличающееся от (70.4) только знаком минус перед единицей

179

где

O

p

= -1 +e cosj,

(70.29)

r

 

 

p и e даются прежними формулами (70.5). Расстояние в перигелии

 

 

 

r =

p

.

 

 

(70.30)

 

 

 

 

 

 

 

 

 

r

 

min

e -1

 

 

 

 

Dj

 

 

 

 

 

 

Половина полного угла поворота при пролете

 

j

 

 

 

 

 

 

 

 

Dj = p

 

 

 

мимо отталкивающего центра

-j.

 

 

 

 

 

 

 

 

 

2

 

 

 

При r = ¥ из (70.29) имеем

cosj = 1

, то-

Рис. 54

гда

 

 

 

 

 

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

æ

ö

 

 

 

1

 

 

 

 

 

çp

÷

= cosj =

.

(70.31)

 

 

sin Dj = sinç

-j÷

 

 

 

ç

÷

 

 

 

e

 

 

 

 

 

è2

ø

 

 

 

 

 

 

Из сравнения (70.31) и (70.25) следует, что углы отклонения для притягивающего и отталкивающего потенциалов одинаковы! Это весьма неочевидный результат. Половинный угол поворота для частицы, летящей с бесконечности, будет даваться той же формулой (70.27). Разница только в том, что в случае притяжения частица отклоняется в сторону центра притяжения, а в случае отталкивания в обратную сторону.

Конические сечения.

Выше было декларировано, что формула (70.4) описывает эллипсы, параболы, гиперболы. Покажем это. Учитывая, что cos x / r , запи-

шем (70.4) в виде

 

r = p -xe .

(70.32)

Возводя в квадрат, получаем

 

 

x2 y2

p2 2 pxe x2e2 .

(70.33)

Далее, после перегруппировки членов получается

1) e 1 (эллипс)

180

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]