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

IvanovMaksyuta

.pdf
Скачиваний:
37
Добавлен:
14.02.2016
Размер:
2.63 Mб
Скачать

 

 

L

dρ

 

 

 

Δϕ = 2ρρminmax

 

ρ2

 

 

.

 

 

 

 

 

 

2

 

2m E U (ρ)

L

 

ρ2

 

 

 

 

Умова замкненості траєкторії вимагає, щоб цей кут дорівнював раціональній частині числа 2π , тобто має бути Δϕ = 2πmn , де

m і n – цілі числа. Тоді через n періодів зміни ρ (від ρmin до ρmax і потім до ρmin ) радіус-вектор ρG, зробивши m повних

обертів, збігається зі своїм початковим значенням. Наведемо приклади (рис. 2, а–г).

ρGmax

ρGmax

ρGmin

ρGmin

a) m =1, n =1, Δϕ = 2π

б) m =1, n = 2, Δϕ = π

 

ρGmax

ρGmax ρGmin

G

 

ρmin

в) m = 2, n =1, Δϕ = 4π

г) m = 2, n = 3, Δϕ = 4π 3

Рис. 2. Приклади замкнених траєкторій для різних центральних полів

101

У загальному випадку, коли Δϕ ≠ 2πmn , траєкторії фінітно-

го руху за нескінченний час заповнюють все кільце між двома граничними колами. На рис. 3 це зображено схематично у вигляді розетки.

Зауважимо, що існують лише два типи центральних полів, у яких всі траєкторії фінітних рухів замкнені. Це поля, у яких по-

тенціальна енергія частинки пропорційна 1r або r2 . Причину цього ми з'ясуємо в наступному розділі.

ρGmax

ρGmin

Рис. 3. Приклад незамкненої траєкторії

Із наведених міркувань можна зробити висновок, що траєкторія руху в центральному полі симетрична щодо апсидальних векторів. У точці повороту квадратний корінь

2

E U (ρ)

L2

 

m2ρ2

m

 

змінює знак. Якщо відраховувати кут ϕ від напрямку радіус-

вектора, проведеного в точку повороту, то відрізки траєкторії з двох сторін будуть відрізнятися лише знаком ϕ за однакових

значень ρ, тобто траєкторія симетрична щодо цього напрямку. Це стосується й інфінітних траєкторій, які утворюються двома симетричними гілками, що йдуть від точки повороту ρmin до

нескінченності. Це можна довести і більш строго. Випишемо рівняння Лагранжа

102

mρ−

L2

= −

U

= f (ρ) .

mρ3

∂ρ

 

 

 

З урахуванням d dt = (L mρ2 )ddϕ це рівняння щодо змінної u =1ρ переписується у такому вигляді:

L2u2 d 2u

 

= − f (u) .

 

 

 

 

+u

 

 

2

m

 

 

 

dϕ

 

 

 

Обираючи систему координат таким чином, щоб полярна вісь проходила через один із апсидальних напрямків (кут ϕ дорів-

нює нулю в цьому напрямку), виконуємо відбиття від полярної осі, тобто замінюємо ϕ на −ϕ . Одержане диференціальне рів-

няння інваріантне щодо такого перетворення, що й доводить симетричність траєкторії щодо апсидальних векторів.

"Падіння" частинки в центр можливе лише тоді, коли потенціальнаенергіядостатньошвидкопрямуєдо −∞ за ρ → 0 . Ізнерівності

mρ2

= E U (ρ)

L2

> 0

2

2mρ2

 

 

або

ρ2U (ρ) + L2 < Eρ2

2m

випливає, що ρ може набувати малих значень лише за умови

lim ρ2U (ρ) < −L2 2m , тобто U (ρ) має прямувати до −∞ або як

ρ→0

−αρ2 з α > L22m , або як 1ρn з n > 2 .

Зауважимо, що випадок центрально-симетричного поля є частинним випадком важливої задачі з двома ступенями вільності, в якій у зв'язку з наявністю циклічної змінної існує, як мінімум, два незалежних інтеграли руху. Нехай функція Лагранжа записується таким чином:

L(q1, q1, q2 ) = 12 a11 (q1 )q12 + 2a12 (q1 )q1q2 + a22 (q1 )q22 U (q1 ) ,

тобто q2 – циклічна координата. Тоді відповідний узагальнений імпульс p2 = ∂Lq2 = a22 (q1 )q2 + a12 (q1 )q1 буде інтегралом

103

руху. Виражаючи з нього швидкість q2 і підставляючи її в інший інтеграл руху – повну енергію, одержуємо

E = aеф2(q1 ) q12 +Uеф (q1 ) ,

де

aеф (q1 ) = (a11a22 a122 )a22 , Uеф (q1 ) =U (q1 ) + p22 2a12 (q1 ) .

При цьому слід зазначити, що aеф (q1 ) має бути додатним. Ця

умова, очевидно, виконується, оскільки кінетична енергія є додатно визначеною квадратичною формою узагальнених швидкостей, що за критерієм Сильвестра можливо лише за умов

a11 (q1 ) > 0, a22 (q1 ) > a122 (q1 )a11 (q1 ) . Для центрального поля

U (r) маємо q1 , q2 = ϕ , p2 = L = mρ2ϕ, a12 = 0 , aеф = a11 = m , a22 = mρ2 , і ми повертаємось до вже одержаних вище формул.

§2. Задача Кеплера

Убагатьох випадках рух відбувається під дією сили, що має вигляд FG = krGr3 . Таким є рух у ньютоновій теорії гравітації,

коли точка з масою m рухається в полі тяжіння іншої, нерухомої точки з масою M , розташованої в центрі системи координат. У цьому випадку k = GMm . Таким є і рух частинки з електричним зарядом q у електричному полі іншого, нерухомого

заряду Q , розташованого в центрі системи координат. Сили та-

кого типу називаються ньютонівськими (або кулонівськими).

Задача про рух матеріальної точки в полі ньютонівських сил називається задачею Кеплера. Уперше її розв'язав у 1740 р. швейцарець Й. Бернуллі.

Розглянемоспочаткуполепритягуваннязпотенціаломвзаємодії

U (r) = − αr ,

де α > 0 . Тоді у площині ( x, y) "ефективна" потенціальна енергія записується у вигляді

104

Uеф (ρ) = −

α

+

L2

.

ρ

2mρ2

 

 

 

Вона зображена на рис. 4, з якого видно, що за ρ → 0 функція Uеф (ρ) → +∞ , а за ρ → +∞ функція Uеф (ρ) 0 . За ρ0 = L2αm вона набуває мінімального значення E0 =Uеф (ρ0 ) = −α2m2L2 .

Очевидно, що за E 0 рух буде інфінітним, а за E0 E < 0 – фінітним. Це можна простежити і за поведінкою фазових траєкторій у площині (ρ,ρ) при зміні енергії від E0 до +. Вони також зображені на рис. 4. За енергії E0 величина ρ = 0 , а ρ = const . Зро-

зуміло, що фазовій траєкторії (точка 0) за цієї енергії відповідає рівномірний рух по колу. Покажемо тепер, що фазовій траєкторії 1 за енергії E1 відповідає рух по еліпсу, фазовій траєкторії 2 за

E2 = 0 відповідає рух по параболі, а траєкторії 3 за E3 відповідає рух по гіперболі. Ця відповідність відображена також на рис. 4.

 

Uеф (ρ)

 

 

 

E2

= 0

E3

 

ρ

 

E1

 

3

 

 

 

E0

 

 

 

 

 

2

 

 

 

ρ

 

1

0

 

 

 

 

 

 

0

1 2 3

ρ

 

 

 

 

 

 

Рис. 4. Схематичне зображення ефективного потенціалу Uеф (ρ)

фазових траєкторій і відповідних їм траєкторій руху у випадку ньютонівського центрального поля

105

Формутраєкторіїможназнайтизадопомогоюзагальноїформули

ϕ−ϕ0 =

Ldρ

 

 

 

.

 

 

 

 

ρ2 2m E U

(ρ)

L2

 

 

ρ2

 

 

 

Підставляючи в неї U (ρ) = −αρ , виконуємо таке інтегрування:

 

 

 

 

 

 

L

dρ

 

 

 

 

 

 

 

 

 

 

 

 

ϕ−ϕ0 =

 

 

 

 

ρ2

 

 

 

 

 

 

 

 

 

L2

 

2mα

 

 

 

 

2mE

 

 

 

 

 

 

 

 

ρ2

ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

L

dρ

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

ρ2

 

 

 

 

 

 

 

m2α2

 

L

 

mα

2

 

2mE +

 

 

 

2

 

 

 

 

 

 

ρ

 

 

 

 

L

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

dρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

ρ2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L2

 

 

2mα

 

m2α2

 

 

m2α2

 

 

 

2mE

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

ρ2

 

 

ρ

L2

L2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

mα

 

 

m2α2

 

 

 

 

 

d

 

 

 

 

 

L

 

 

2mE +

 

 

L2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

= −

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

mα

 

 

 

m2α2

2

 

 

 

 

 

 

2mE +

 

 

 

 

1

 

ρ

L

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

mα

 

 

 

 

 

 

 

 

 

= −

 

dx

 

 

,

 

x =

 

 

 

ρ

 

 

L

 

 

.

 

 

 

 

 

 

 

1 x2

 

 

2mE +

m2α2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

Обираючи початок відліку кута ϕ так, щоб ϕ0 = 0 , одержуємо

 

 

 

x = cosϕ , або

 

L

 

mα

 

 

2mE +

m2α2

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

cos ϕ.

 

 

 

 

ρ

 

L

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

Розв'язуючи далі щодо ρ, дістаємо

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L2

 

 

 

 

 

p

 

 

ρ =

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

mα

 

 

 

 

 

=

,

mα + 2mE + m

α2

 

 

 

 

 

 

 

 

1 + 2EL2

 

 

1 + ecos ϕ

 

cosϕ 1 +

cosϕ

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

L

 

L

 

 

 

 

 

 

 

 

 

 

 

 

mα

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

106

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

де p = L2mα і e = 1 + 2EL2mα2 – т. зв. параметр і ексцентри-

ситет траєкторії. Із останньої формули видно, що ρ = p = L2mα за e = 0 . Це відповідає енергії E = E0 = −α2m2L2 , тобто маємо фінітний рух по колу радіуса ρ = p = L2mα . За E0 < E < 0 ексцентриситет 0 < e <1 , тобто також маємо фінітний рух, який відбувається по еліпсу (рис. 5). Із формули ρ = p(1 + ecosϕ) і

рис. 5 видно, що найменша та найбільша відстані (відповідно перигелій та афелій) дорівнюють

ρmin =

 

 

p

= a

(1 e),

ρmax =

 

 

p

 

= a (1 + e) ,

1

+ e

1

e

де a = p (1 e2 ) = α 2 E

– велика піввісь еліпса. Мала ж піввісь

визначається із трикутника, виділеного на

рис. 5, і дорівнює

b = p 1 e2 = L 2m E .

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

x

2b

 

 

 

 

 

ϕ

 

 

2ae

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρmax

ρmin

 

Рис. 5. Схематичне зображення еліптичної орбіти

За E 0 маємо інфінітний рух. Якщо E > 0 , то e >1 і траєкторією руху є гіпербола, яка огинає центр поля, розташований

у фокусі

O

(рис. 6). Використовуючи формули c2 = a2 +b2 ,

 

1

одержуємо перигелій ρmin = p (e +1) = a (e 1) ,

e = c a >1 ,

 

 

107

де a = p (e2 1) = α 2E – "піввісь" гіперболи. У випадку

E = 0 ексцентриситет e =1 , тобто частинка рухається по параболі, перигелій якої дорівнює ρmin = p2 . Цей випадок ре-

алізується, якщо частинка починає свій рух зі стану спокою на нескінченності.

y

O1

b c

x

 

a

 

Рис. 6. Схематичне зображення гіперболічної траєкторії

Звертаємось тепер до руху матеріальної частинки в полі відштовхування U (ρ) = αρ, де α > 0 . У цьому випадку ефективна потенціальна енергія

Uеф (ρ) =

α

+

L2

ρ

2mρ2

 

 

монотонно зменшується від +∞ до нуля при зміні ρ від нуля до +∞ . Енергія частинки завжди більша нуля, а рух завжди інфінітний. Траєкторією є гіпербола ρ = p(1 + ecosϕ) . Вона проходить

повз центр поля навідстані ρmin = p(e 1) = a (e +1) віднього. Знайдемо тепер залежність координат частинки від часу:

 

m

 

 

ρdρ

 

 

 

 

 

m

ρdρ

 

t =

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

,

2

 

E

−ρ

2

+

αρ

L2

2

 

E

a2e2 (ρ −a)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

 

2m

 

E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

108

оскільки a = α 2 E , b = L 2m E , a2 b2 = a2e2 . За допомогою підстановки ρ− a = −aecos ξ цей інтегралприводиться до вигляду t = maα3 (1 ecos ξ)dξ = maα3 (ξ− esin ξ) +C .

Обираючи початок відліку часу так, щоб C = 0 , одержуємо параметричне представлення залежності ρ від t :

ρ = a(1 ecos ξ),

t =

ma

3

(ξ −esin ξ),

0 ≤ ξ ≤ 2π .

α

 

Через цей параметр

ξ можна виразити і декартові координати

частинки x = ρcosϕ ,

y = ρsin ϕ :

 

 

 

 

ρ(1 + ecos ϕ) = p,

ρ + ex = p,

ex = p −ρ =

= a (1 e2 ) a (1 ecos ξ) = ae(cos ξ− e),

x = a (cos ξ− e),

y = ρ2 x2

= a 1

e2 sin ξ.

Особливістю задачі Кеплера є те, що в полі U (ρ) = αr (незалежно від α ) існує ще один інтеграл руху, специфічний для

цього поля. Це т. зв. вектор Лапласа–Рунге–Ленца:

AG = vG×LG+ αrG . r

Щоб продемонструватиG сказане, знайдемо повну похідну від вектора A за часом

G

G

G

αvG

αrG(vG rG)

G

G G

 

αvG

 

αrG(vG rG)

 

A

= v

×L +

r

 

= m(v ×(r ×v )) +

r

 

 

=

r3

r3

 

 

 

G G G

G G G

αvG

αrG(vG rG)

 

 

 

 

 

= mr (v v ) mv (v r ) +

r

 

 

.

 

 

 

 

r3

 

Ураховуючи рівняння руху mvG = αrGr3 , одержуємо AG = 0 , тобто цей вектор зберігається. Знайдемо величину та напрямок цього вектора, обчислюючи скалярний добуток вектора εG = Aα з вектором rGmin = p(1 + e) при α < 0 :

109

G G

1 G

G

G

rG 2

 

L2

ε rmin =

 

rmin ×(rmin ×L) +

min

=

 

+ rmin =

α

r

mα

 

 

 

 

 

min

 

 

 

= −p + rmin = rmin rmin (1 e) = rmine,

G

тобто модуль вектора A дорівнює A = eα, а за напрямком він

паралельний rGmin . Наявність цього додаткового інтеграла (неза-

лежного від двох стандартних інтегралів – енергії E і моменту імпульсу L) визначає той факт, що всі траєкторії фінітного руху для задачі Кеплера є замкненими. Механічні системи, які мають додаткові інтеграли руху і замкнені траєкторії, називаються

системами з випадковим виродженням.

Ще один приклад такої задачі у випадку центрально-симетрич- ного поля відповідає потенціалу U (r) = kr2 2 . Походження до-

даткового інтеграла в цьому випадку легко зрозуміти, розглянувши цю задачу в декартових координатах, у яких вона зводиться до випадку ізотропного двовимірного осцилятора (див. лекцію9). Тут можна сказати, що два інтеграли руху – це енергії E1 і E2 , які збе-

рігаються незалежно, а додатковим інтегралом руху є проекція моменту імпульсу навісь z, що пропорційна xy yx .

Поняття про випадкове виродження відіграє велику роль у фізиці й за межами класичної механіки. Справа в тому, що відповідні інтеграли руху існують і у квантовій механіці. У результаті обидві вищевказані задачі у квантовій механіці (особливо важлива кван- тово-механічна задача Кеплера, рішення якої визначає енергетичний спектр атомаводню) мають важливі особливості.

Контрольні запитання та завдання

1.Інтеграли руху у випадку центральних полів.

2.Типи центральних полів, в яких всі траєкторії фінітних рухів є замкненими. Випадкове виродження.

3.Умови "падіння" частинки в центр.

4.Види траєкторій руху у випадку задачі Кеплера.

5.Вектор Лапласа–Рунге–Ленца.

110

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