Feynmann - Feynmann 1

Однако в отличие от предыдущей задачи нам потребуются теперь компоненты ускорения в двух направлениях, которые мы назовем х и у. Положение планеты в данный момент будет определяться координатами х и у, поскольку третья коорди­ната z всегда равна нулю.

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

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

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

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

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

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

Перейдем к дальнейшему. Возьмем интервал времени e=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 3=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. График движения планеты вокруг Солнца.

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

Давайте посмотрим теперь, как вычислить движение Непту­на, Юпитера, Урана и остальных планет. Можно ли сделать подробные расчеты со множеством планет, учитывая к тому же и движение Солнца? Разумеется, можно. Найдем сначала силу, действующую на каждую данную планету, например на ту, которую мы обозначим номером i и координаты которой х i, y iи z i(i=1 может означать Солнце, i=2 — Меркурий, i=3 — Венеру и т. д.). Наша задача — найти координаты всех планет. По закону тяготения x-компонента силы, действующая на i-ю планету со стороны планеты номер j' с координатами х jу j, z j, будет равна —Gm im j(x i-x j)/r 3 jj . Если же учесть силы со стороны всех планет, то получим следующую систему уравнений: