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

Mech-Slobod

.pdf
Скачиваний:
263
Добавлен:
19.03.2015
Размер:
14.21 Mб
Скачать

261

Якщо припустити, що сила, яка прикладена до яблука з боку Землі й спричинює прискорення його вільного падіння на Землю, має таку саму природу, що й сила, яка прикладена до Місяця з боку Землі й спричинює його рух по орбіті, то для яблука можна записати

m

я

a

я

= G

mяM З

,

(7.11)

R2

 

 

 

 

 

 

 

 

 

 

 

 

З

 

 

де

 

mя

маса яблука, а RЗ радіус Землі, якому

з великою точністю дорівнює відстань між

яблуком і центром Землі.

З рівностей (7.10) та (7.11) випливає просте співвідношення між прискореннями яблука та Місяця

a

я

 

R

З

 

2

 

 

=

 

 

,

(7.12)

 

 

 

 

aM

 

 

 

 

 

 

 

RЗМ

 

 

 

яке можна перевірити на підставі даних спостережень за рухом Місяця та за відомою

величиною прискорення вільного падіння, з яким падає яблуко aя

= g = 9,81 м c2 .

 

 

Період обертання Місяця TM = 27,32 доби,

а радіус його орбіти

RЗМ = 384 тис. км. Величину

доцентрової складової прискорення Місяця aM можна обчислити за формулою (7.3),

 

 

a

M

=

4π 2

R

ÇÌ

= 2,72 103 ì / ñ2 .

 

 

 

 

 

 

 

T 2

 

 

 

 

 

 

 

 

 

 

M

 

 

 

 

 

 

 

 

 

 

Тоді відношення прискорень яких набуває яблуко та Місяць внаслідок притягання з боку Землі

є

aÿ / aM

= (9,81 ì

c2 ) /(2,72 103 c2 ) ≈ 3,6 103 = 602 .

В

той

же

час

відношення

R

/ R = (384 103

км) /(6,36 103 км) = 60,38 .

Таким

чином,

виявляється,

що

відношення

ЗМ

 

З

 

 

 

 

 

 

 

 

 

 

 

прискорень яблука та Місяця з великою точністю дорівнює квадрату оберненого відношення відстаней від центра Землі до яблука і до Місяця у повній відповідності до формули (7.12), що є вагомим аргументом на користь зробленого припущення про однакову природу сил тяжіння між небесними тілами і силою тяжіння поблизу поверхні Землі. Саме це дає підстави називати закон яжіння (7.1) всесвітнім або універсальним6, оскільки він справедливий для всіх без виключення тіл у Всесвіті.

6Термін універсальний походить від латинського universalis всеосяжний, загальний, який, у свою чергу, походить від латинського universum загальне, Всесвіт.

262

Необхідно підкреслити, що так звані гравітаційні маси m1 та m2 , що фігурують у формулі

закону всесвітнього тяжіння (7.1) є мірою гравітаційних властивостей тіл, тобто їх здатності притягуватись одне до одного незалежно від їх інших фізичних властивостей. Вперше силу гравітаційної взаємодії виміряв Кавендіш у 1798 р. у зв’язку з вимірюванням гравітаційної сталої G .

Гравітаційна стала G не може бути визначена з астрономічних спостережень, оскільки для її визначення за формулою (7.1) необхідно знати величину сили, з якою притягуються небесні тіла та гравітаційні маси цих тіл. Кавендішу вдалося безпосередньо виміряти надзвичайно малу силу гравітаційного притягання між двома свинцевими кулями за допомогою винайдених ним крутильних терезів. Той факт, що всі тіла незалежно від будь-яких їхніх властивостей, в тому числі незалежно від величини їх інертної маси, падають біля поверхні Землі з однаковим прискоренням дає можливість припустити, що гравітаційна маса пропорційна інертній. Дійсно, з одного боку на підставі другого закону Ньютона для тіла, що вільно падає біля поверхні Землі, можна записати

mінa = G mгрM З ,

RЗ2

а з іншого боку відомо, що для всіх тіл

 

 

(7.13)

 

 

a = g . Тоді виявляється, що

m = G

M З

m

 

або

gR2

 

 

ін

 

гр

 

 

 

З

 

 

 

mін = kmгр , тобто інертна маса пропорційна гравітаційній і навпаки. Цей факт є чисто експериментальний. Він дозволяє обрати коефіцієнт пропорційності k = 1 ї вважати інертну та гравітаційну маси рівними одна одній, для їх вимірювання використовувати один і той же самий еталон, і встановити для них одну й ту саму одиницю вимірювання, а саме кілограм.

Оскільки одиницею гравітаційної маси є кілограм, а одиниці сили було встановлено раніше на підставі формули другого закону Ньютона, то гравітаційна стала G виявляється розмірною величиною: її розмірність прямо випливає із формули закону всесвітнього тяжіння (7.1). Сучасне значення гравітаційної сталої G = (6,6742 ± 0,0010) ×1011 Н м2 кг2 . Треба сказати, що точність вимірювання G не набагато зросла з часів Кавендіша, оскільки при зовнішній простоті постановки експерименту дуже складно усунути всі фактори, які можуть вносити похибку в обчислену величину

G .

263

7.2. Рух частинки в полі центральної сили, обернено пропорційної квадрату відстані до силового центру

7.2.1 Рух частинки в полі центральної сили

Застосуємо закони збереження моменту імпульсу та енергії до розгляду руху частинки в полі центральної сили.

Оскільки момент центральної сили F = f (r) rr відносно центра силового поля завжди дорівнює

нулю:

M = [r × F ] = [r × ( f (r)

r

)] =

f (r)

[r × r ] = 0 ,

(7.14)

 

r

 

r

 

 

то момент імпульсу частинки в полі центральної сили, визначений відносно центра поля,

зберігається,

L = [r × p] = const .

Оскільки центральна сила є консервативна сила, то повна механічна енергія частинки в полі центральної сили також зберігається:

E = mυ 2 +U (r) = const .

2

Оскільки напрям вектора L є незмінним у просторі, а радіус-вектор частинки r та вектор її імпульсу p , а отже і вектор швидкості, перпендикулярні до L , траєкторія частинки повністю лежить в одній площині, перпендикулярній до вектора L . Рівняння цієї площини є (L r ) = 0 , тобто вона

проходить через центр поля (силовий центр).

Оберемо напрям осі OZ декартової системи координат вздовж вектора L . Тоді відмінною від нуля буде лише проекція вектора L на цю вісь, Lz = L . Траєкторія частинки при такому виборі

системи координат повністю лежить в площині XOY . У координатному записі маємо такі вирази для величини моменту імпульсу та для повної механічної енергії частинки

L = m(xy yx)

 

 

 

1

 

 

 

 

(7.15)

 

m(x

2

+ y

2

E =

2

 

 

) + U (r)

 

 

 

 

 

 

264

Оскільки силове поле центральне, тобто і сила і потенціальна енергія залежать лише від відстані до силового центру, зручно в площині траєкторії ввести полярну систему координат з початком у силовому центрі O.

Якщо полярну вісь полярної системи координат спрямувати вздовж осі OX декартової системи координат, то між x, y та r,ϕ існує зв’язок

x = r cosϕ y = r sin ϕ .

Відповідно

x = r cosϕ r sinϕ ϕ y = r sinϕ + r cosϕ ϕ .

Тоді

L = m[r cosϕ (r sinϕ + r cosϕ ϕ ) r sinϕ (r cosϕ r sinϕ

E = 1 m[(r cosϕ r sinϕ ϕ )2 + (r sinϕ + r cosϕ ϕ )2 ] + U (r2

і після нескладних алгебраїчних перетворень остаточно

(7.16)

(7.17)

ϕ )] = mr 2ϕ

 

1

2

2

2

(7.18)

) =

2 m(r

 

+ r ϕ

 

) + U (r)

L = mr 2ϕ

 

 

 

 

(7.19)

 

1

 

E =

m(r 2

+ r 2ϕ 2 ) + U (r)

2

 

 

 

З виразу для моменту імпульсу (7.19а) випливає, що так звана секторіальна швидкість частинки,

тобто відношення

площі

dσ =

1 r 2 dϕ , яку «замітає» радіус-вектор частинки

за елементарний

 

 

 

 

 

 

 

 

 

2

 

 

проміжок часу dt , до величини цього проміжку є величина стала,

 

dσ

 

d

1

2

 

1

2

 

L

 

dt

=

 

( 2 r

 

dϕ ) = 2 r

ϕ =

 

= const.

(7.20)

dt

 

2m

Це твердження становить зміст так званого другого закону Кеплера, який встановив його з аналізу спостережень за рухом планети Марс по еліптичній орбіті навколо Сонця в центральному гравітаційному полі останнього, функція потенціальної енергії для якого, як було показано в Розділі 5

(див. формулу (5.36)) є U (r) = αr .

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

265

ϕ =

 

 

L

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(7.22)

 

mr 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

де всі величини в правій частині є додатно визначеними.

 

 

 

 

 

За допомогою рівності

(7.22) можна виключити ϕ з виразу для повної енергії частинки (7.19) і

подати останню у вигляді

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

 

2 L

 

2

 

1

 

2

 

L2

 

 

1

 

2

 

 

E =

 

 

m r

 

+ r

 

 

 

 

+ U (r) =

 

mr

 

+

 

 

+ U (r) =

 

mr

 

+ Uеф(r) ,

(7.23)

2

 

 

2

2

 

2mr

2

2

 

 

 

 

 

 

mr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

де введено так звану ефективну потенціальну енергію

 

 

 

 

 

Uеф

(r) =

 

L2

+ U (r) .

 

 

 

 

 

 

 

 

 

 

 

 

(7.24)

2mr2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким чином, вираз для повної енергії частинки в центральному полі має такий самий вигляд як і для повної енергії частинки при одновимірному русі (див. формулу (5.56)), де потенціальну

енергію U (r) замінено на потенціальну енергію Uеф(r) . Отже рух частинки в радіальному

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

Uеф(r) . Швидкість такого руху є

dr = ±

2 (E U еф )

(7.25)

dt

m

 

 

 

 

і вона обертається на нуль при умові

 

E = Uеф(r) =

L2

 

+ U (r) .

(7.26)

2mr2

 

 

 

 

На відміну від випадку одновимірного руху розглянутого в Розділі 5 частинка при цьому не

зупиняється,

оскільки

завжди існує азимутальна складова швидкості

rϕ . Тому радіальний рух

частинки в центральному полі називають квазіодновимірним.

Рівняння траєкторії частинки в полярних координатах можна знайти з рівності (7.22),

виключивши з неї час. Для цього перепишемо її у вигляді

 

dϕ =

L

dt ,

(7.27)

mr 2

 

 

 

підставимо в неї dt , знайдене з рівності (7.25),

266

dt = ±

dr

(7.28)

 

2

(E Uеф)

 

 

m

 

 

 

 

і отримаємо зв’язок між диференціалами dϕ та dr

 

 

 

 

 

 

L

 

dr

 

 

dϕ =

 

 

 

 

mr2

 

 

 

2

 

.

(7.29)

 

 

 

 

 

 

 

 

 

[E Uеф(r)]

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

Для знаходження зв’язку між полярним

кутом ϕ та полярним радіусом r треба

взяти

невизначені інтеграли від лівої та правої частин рівняння (7.29)

 

 

 

 

 

 

 

L

 

dr

 

 

 

 

 

 

 

mr 2

 

 

 

ϕ = ±

2

 

 

+ ϕ 0 ,

(7.30)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[E Uåô (r)]

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

де через ϕ0 позначено сталу інтегрування.

 

 

Таким чином,

 

ми отримали ϕ = ϕ (r) ,

тобто рівняння траєкторії у явному вигляді,

який

залежить від конкретного виразу для потенціальної енергії U (r) .

Ясно, що при наявності залежності ϕ = ϕ (r) можна, в принципі, знайти й обернену залежність

r = r(ϕ ) .

Як і у випадку істинного одновимірного руху, квазіодновимірний рух можливий лише при

таких r , при яких E Uеф 0. При цьому він може бути як інфінітним, якщо область можливих значень r обмежена лише знизу, r rmin , або фінітним, якщо rmax r rmin . При цьому характер

руху і розміри області, в якій він можливий, визначаються конкретним виглядом потенціальної енергії U (r) і повною механічною енергією частинки. Приклади руху частинки в полі з

потенціальною енергією U (r) = rαn при різних значеннях α та n подано на Рис.7.5.

Необхідно відзначити, що навіть при фінітному русі, коли траєкторія частинки повністю лежить в кільцевій області rmax r rmin , вона в загальному випадку не є замкненою. Умова замкненості

може бути отримана з таких міркувань. Одновимірний радіальний фінітний рух частинки є

267

Рис. 7.5. Рух частинки в центральному полі

періодичним. Цей період дорівнює проміжку часу, за який частинка, починаючи рухатись від

внутрішньої границі rmin області доступної для руху досягає зовнішньої границі rmax ,

а потім

повертається назад до внутрішньої границіrmin області доступної для руху.

 

268

За цей час радіус-вектор частинки повернеться на кут

 

 

 

 

L

 

 

rmax

 

 

 

dr

 

 

 

 

2mr2

(7.31)

ϕ = 2

2

.

 

rmin

[E Uеф(r)]

 

 

 

m

 

 

 

 

 

 

 

Після n «ходінь» частинки від rmin до rmax і назад радіус-вектор частинки повернеться на кут

nϕ .

Якщо nϕ = m2π , то після n «ходінь» частинки від rmin до rmax і назад та здійснення m

повних

обертів радіус-вектора частинки (тобто зміні азимутального кута ϕ

на 2πm ) частинка

прийде в ту саму точку простору і траєкторія замкнеться. Таким чином, траєкторія частинки буде

2π

замкненою, якщо відношення ϕ є раціональне число. Відомо два типи центральних полів, у яких

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

віддалі від силового центра, U (r) ~ r1 та поля, в яких потенціальна енергія частинки прямо

пропорційна квадрату віддалі від силового центра U (r) ~ r2 . До першого типу належать

гравітаційне поле точкової маси та електростатичне поле точкового заряду, а до другого поле квазіпружної сили.

7.2.2 Рух частинки під дією центральної сили, обернено пропорційної квадрату відстані до силового центру

Дослідимо тепер рух частинки в полі з U (r) = αr , зокрема знайдемо явний вигляд можливих

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

Uеф(r) = α

+

 

L2

.

 

 

 

 

 

 

(7.31)

 

2mr2

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

записане у вигляді

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

dr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕ ϕ0 = ±

 

 

 

 

r2

 

 

 

.

(7.32)

 

 

 

 

α

 

L2

 

 

 

 

 

 

 

 

 

 

 

 

2m E

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

r

 

2mr

 

 

 

 

 

 

 

 

 

 

 

 

269

Обчислимо інтеграл у правій частині (7.32):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dr

 

 

 

 

 

 

 

 

1

 

 

 

 

dr

 

 

L

 

 

dy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕ = ±

 

 

 

α

 

 

L2

 

 

 

= y = r

, dy = − r2

= 2m

 

 

L2

 

 

2

=

 

 

 

 

 

 

 

 

 

 

 

2m E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E αy

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

2mr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

 

 

dy

 

 

 

 

 

 

 

 

L

y +

α

2m

= z

 

 

L

 

 

 

2m

dz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2m

2

L2

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

L

 

 

=

2m

 

 

 

L

 

y + α

 

2m

 

2

α

2

 

 

 

2m

 

 

=

2m

 

 

α

2

m

 

 

 

 

 

 

 

 

 

 

dy

=

dz

 

 

 

 

E +

 

z

2

 

 

 

 

 

E

 

 

 

 

 

+

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2m

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2L

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

2L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dz

 

= arccos

z

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

2

z

2

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(7.33)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При обчисленні інтеграла (7.33) було зроблено заміни змінних

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y =

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(7.34а)

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

та

z =

L

y + α

2m2 ,

 

2m

2

L

а також позначено

a2 = E +

α 2m2 .

 

 

 

2L

 

Отже, рівняння (7.32) можна подати у вигляді

θ = arccos az ,

де позначено θ = ϕ ϕ0 , тобто перенесено початок відліку полярного кута.

Знайдемо косинус від обох частин рівності (7.35)

az = cosθ .

Після підстановки виразів для z , y і a у рівняння (7.36) одержимо зв’

кутом θ та полярним радіусом r

(7.34б)

(7.34)

(7.35)

(7.36)

язок між полярним

L 1

+

α

2m2 =

E +

α 2 m2 cosθ .

(7.37)

2m r

 

2

L

 

2L

 

270

 

Для приведення (7.37) до вигляду r = r(θ )

виконаємо ряд перетворень

1 =

 

2mE

+ α

2

m

2

αm

=

 

1

2

 

=

1

( 1 + e cosθ ),

 

 

cosθ

α m

+ 1+ 2EL

cosθ

 

r

 

L2

 

L4

 

L2

 

L2

 

α 2m

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

де введено позначення

 

 

 

 

 

 

 

 

 

 

 

α

 

m

=

1

 

 

 

 

 

 

 

 

 

 

 

 

 

(7.38)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L2

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

та

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 +

2EL2

= e .

 

 

 

 

 

 

 

 

 

 

 

(7.39)

 

 

 

 

 

 

α 2m

 

 

 

 

 

 

 

 

 

 

 

 

 

Остаточно маємо рівняння траєкторії в полярних координатах

 

 

 

r =

 

 

p

 

.

 

 

 

 

 

 

 

 

 

(7.40)

1+ ecosθ

 

 

 

 

 

 

 

 

 

Рівняння (7.40) є відомим з аналітичної геометрії рівнянням конічних перерізів, яке зв’язує так

званий фокальний радіус-вектор r з полярним кутом θ . Залежно від значень фокального параметра p та ексцентриситету e це рівняння може визначати еліпс ( e < 1), параболу ( e = 1) або гіперболу

( e > 1).

У випадку притягання частинки до силового центру (α < 0 ) у рівнянні (7.40) необхідно брати

знак «+»

r =

p

,

(7.41)

1+ ecosθ

а у випадку відштовхування частинки від силового центру (α > 0 ) необхідно брати знак «»

r =

p

 

1+ ecosθ .

(7.42)

Аналіз різних випадків руху частинки відповідно до (7.41) і (7.42) подано на Рис. 7.6, у верхній

частині якого наведено графіки залежності ефективної потенціальної енергії Uеф від відстані до

силового центру r , а нижче при збереженні одного й того ж масштабу по осі абсцис зображено траєкторії частинок при різних значеннях їх повної енергії E як для випадку притягання(α < 0 ), так і для випадку відштовхування (α > 0 ). Там же показано додаткові побудови і наведено параметри,

які використовуються для опису кривих другого порядку (асимптоти, півосі, директриси тощо).

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