Читать «K читателям русского издания» онлайн - страница 108

Действительно, координатная плоскость ху выбрана нами таким образом, что z-компоненты как силы, так и начальной скорости равны нулю, а поэтому нет никаких причин, которые бы заставили планету выйти из этой плоскости. Сила при этом будет направлена по линии, соединяющей планету с Солнцем, как это показано на фиг. 9.5.

Фиг. 9.5. Сила притяжения, действующая на планету.

Из этого рисунка видно, что горизонтальная компонента силы так относится к полной ее величине, как координата х относится к расстоянию r. Это сразу следует из подобия тре­угольников. Кроме того, если х положительна, то F_x отрица­тельна, и наоборот.

Таким образом, F_xF=-x/r, или F_я=-Fxlr=-GM mx/r³ и соответственно F_y=-GMmy/r³. Теперь можно воспользо­ваться динамическими законами (9.7) и написать, что х- или y-компонента ускорения, умноженная на массу планеты, равна соответственно х- или y-компоненте силы:

Это именно та система уравнений, которую мы должны решить. Для того чтобы упростить вычисления, предположим, что либо единицы измерения времени или массы выбраны соответствую­щим образом, либо нам просто повезло, словом, получилось так, что GM=1. Для нашего случая предположим, что в на­чальный момент t=0 планета находилась в точке с координа­тами х=0,500 и у=0,000, а скорость ее в этот момент направ­лена параллельно оси у и равна 1,6300. Как же в этом случае делаются расчеты? Снова составляется таблица со столбцами для времени t, координаты х, x-компонент скорости v_x и уско­рения а_х. Затем идут отделенные чертой три колонки: для координаты y, у-компонент скорости и ускорения. Однако, для того чтобы подсчитать ускорения, мы должны воспользо­ваться уравнением (9.17), согласно которому его компоненты равны –х/r³ и –у/r³, а r=(x²+y²). Так что, получив х и у, мы должны где-то в сторонке провести небольшие вы­числения – извлечь квадратный корень из суммы квадра­тов и получить расстояние. Удобно также отдельно вычис­лить и 1/r³.

После этого все готово, чтобы определить компоненты ус­корения. Всю эту работу можно сильно облегчить, если поль­зоваться таблицами квадратов, кубов и обратных величин. На нашу долю останется тогда только умножение х на 1/r³, которое легко выполняется на логарифмической линейке.

Перейдем к дальнейшему. Возьмем интервал времени =0,100. В начальный момент t=0

x(0)=0,500,. у(0)=0,000,

v_x(0) = 0,000, v_y(0)=+1,630.

Отсюда находим

r(0)=0,500, 1/r³=8,000,

a_x=-4,000, а_у=0,000.

После этого можно вычислять компоненты v_x (0,05) и v_y (0,05):

v_я (0,05)=0,000-4,000•0,050 = -0,200,

v_y(0,05)=1,630+0,000-0,100=1,630.

А теперь начнем наш основной расчет:

и т. д.

В результате мы получим числа, приведенные в табл. 9.2, где приблизительно за 20 шагов прослежена половина пути нашей планеты вокруг Солнца. На фиг. 9.6 отложены коорди­наты планеты х и y, приведенные в табл. 9.2.

Фиг. 9.6. График движения планеты вокруг Солнца.

Точки представ­ляют собой последовательные положения планеты через каж­дую десятую долю выбранной нами единицы времени. Видно, что сначала она двигалась быстро, а затем – все медленней и медленней. Видна также и форма кривой движения планеты. Итак, вы теперь знаете, как реально можно вычислять движе­ние планет!