IvanovMaksyuta
.pdf
|
|
L |
dρ |
|
|
|
|
Δϕ = 2∫ρρminmax |
|
ρ2 |
|
|
. |
||
|
|
|
|||||
|
|
|
2 |
||||
|
2m E −U (ρ) − |
L |
|
||||
ρ2 |
|||||||
|
|
|
|
||||
Умова замкненості траєкторії вимагає, щоб цей кут дорівнював раціональній частині числа 2π , тобто має бути Δϕ = 2πm
n , де
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πm
n , траєкторії фінітно-
го руху за нескінченний час заповнюють все кільце між двома граничними колами. На рис. 3 це зображено схематично у вигляді розетки.
Зауважимо, що існують лише два типи центральних полів, у яких всі траєкторії фінітних рухів замкнені. Це поля, у яких по-
тенціальна енергія частинки пропорційна 1
r або r2 . Причину цього ми з'ясуємо в наступному розділі.
ρGmax
ρGmin
Рис. 3. Приклад незамкненої траєкторії
Із наведених міркувань можна зробити висновок, що траєкторія руху в центральному полі симетрична щодо апсидальних векторів. У точці повороту квадратний корінь
2 |
E −U (ρ) − |
L2 |
|
|
m2ρ2 |
||
m |
|
||
змінює знак. Якщо відраховувати кут ϕ від напрямку радіус-
вектора, проведеного в точку повороту, то відрізки траєкторії з двох сторін будуть відрізнятися лише знаком ϕ за однакових
значень ρ, тобто траєкторія симетрична щодо цього напрямку. Це стосується й інфінітних траєкторій, які утворюються двома симетричними гілками, що йдуть від точки повороту ρmin до
нескінченності. Це можна довести і більш строго. Випишемо рівняння Лагранжа
102
mρ− |
L2 |
= − |
∂U |
= f (ρ) . |
|
mρ3 |
∂ρ |
||||
|
|
|
З урахуванням d dt = (L mρ2 )d
dϕ це рівняння щодо змінної 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 з α > L2
2m , або як −1
ρn з n > 2 .
Зауважимо, що випадок центрально-симетричного поля є частинним випадком важливої задачі з двома ступенями вільності, в якій у зв'язку з наявністю циклічної змінної існує, як мінімум, два незалежних інтеграли руху. Нехай функція Лагранжа записується таким чином:
L(q1, q1, q2 ) = 12 a11 (q1 )q12 + 2a12 (q1 )q1q2 + a22 (q1 )q22 −U (q1 ) ,
тобто q2 – циклічна координата. Тоді відповідний узагальнений імпульс p2 = ∂L
∂q2 = 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 = krG
r3 . Таким є рух у ньютоновій теорії гравітації,
коли точка з масою 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 ) = −α2m
2L2 .
Очевидно, що за 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 = L2
mα і e = 1 + 2EL2
mα2 – т. зв. параметр і ексцентри-
ситет траєкторії. Із останньої формули видно, що ρ = p = L2
mα за e = 0 . Це відповідає енергії E = E0 = −α2m
2L2 , тобто маємо фінітний рух по колу радіуса ρ = p = L2
mα . За 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 = p
2 . Цей випадок ре-
алізується, якщо частинка починає свій рух зі стану спокою на нескінченності.
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 = αrG
r3 , одержуємо 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
