S или точно в противоположную сторону. Это свойство сохраняется во все время движения, что лишний раз подчеркивает условность термина «начальная». Остальные три типа орбит — плоские кривые, не содержащие прямолинейных участков.
2. Эллипс (рис.2). Центральное тело S, как ни странно это звучит, находится не в центре эллипса, а в одном из двух его фокусов. Отличие эллипса от окружности измеряется эксцентриситетом е — отношением расстояния между фокусами к длине большой оси. Эксцентриситет окружности равен нулю. Эллипс тем более вытянут, чем ближе е к единице.
Рис.2
3. Парабола (рис.3). По параболе частица уходит в бесконечность. Скорость частицы уменьшается, неограниченно приближаясь к нулю. Фигурально выражаясь, частица уходит в бесконечность и останавливается там.
Рис.3.
4. Гипербола (рис.4). По гиперболе частица уходит в бесконечность, приближаясь к некоторой прямой, асимптоте. Скорость частицы приближается к некоторой положительной величине υ∞— скорости на бесконечности, оставаясь все время больше нее.
Рис.4.
По какой из трех кривых будет двигаться частица зависит от полной механической энергии Е единицы массы, включающей в себя кинетическую Ек и гравитационную потенциальную Ер. Поскольку трения нет, то Е сохраняется во все время движения. Оказывается, частица движется по эллипсу, если Е<0; по параболе, если Е=0; по гиперболе, если Е>0. Напомню, что потенциальная энергия имеет смысл с точностью до постоянного слагаемого. В физике и астрономии это слагаемое принято фиксировать условием Ер=0, когда частица находится бесконечно далеко от S.
При таком соглашении
Ек=υ2/2 и Ер=—К2/r, (4)
где К =√GM, а М — масса S. Если расстояния измерять в километрах, время — в секундах, то К=364305, если S — Солнце; К = 631,35, если S — Земля. На практике часто вместо Е используют более наглядную величину — скорость υ=√2Ек. Критическому значению Е=0 отвечает вторая космическая скорость υII (называемая также скоростью убегания или параболической скоростью). Понятно, что υII — не число, а зависящая от расстояния до S величина: скажем, для спутника Земли υII=11 км/с вблизи поверхности планеты, но υII=1,5 км/с у орбиты Луны. Полезно знать, что первая космическая (круговая) скорость υI и параболическая скорость υII различаются только множителем √2: υII=υI√2≈1,41υI
Между круговой и параболической скоростями есть принципиальная разница. Чтобы двигаться по окружности, круговую скорость следует направить перпендикулярно радиусу-вектору, соединяющему центральное тело и частицу. Чтобы уйти на бесконечность, достаточно развить параболическую скорость; при этом ее направление безразлично, лишь бы избежать столкновения с S.
За исключением специального случая (когда скорость направлена точно к S или точно в противоположную сторону) орбиты оказались кривыми линиями. К тому же, движение по орбитам неравномерно. Самая большая скорость — в перицентре (ближайшей к S точке орбиты), и чем дальше от перицентра, тем она меньше. Наименьшая скорость в случае эллипса — в апоцентре (наиболее удаленной от S точке орбиты).
Дадим количественные соотношения. Расстояние rр от S до перицентра выражается через большую полуось а (среднее расстояние от движущегося тела до S) и эксцентриситет е по формуле rp=а(1—е). Расстояние rа от S до апоцентра rа=а(1+е). Скорости в экстремальных точках (апсидах) эллипса составляют:
υp=υI(a)√(1+e)/√(1—e) и υa=υI(a)√(1—e)/√(1+e)
Здесь υI(a) — круговая скорость на расстоянии от a до S. В свою очередь υI убывает обратно пропорционально квадратному корню из расстояния до S: υI=K/√r.
Между большой полуосью и периодом обращения существует связь, открытая еще И. Кеплером в начале XVII в.:
Р = 2π(а3/2/K) (5)
Разумеется, выражение постоянной К через G и М — заслуга Ньютона.
Если эллипс близок к окружности, различие скоростей в разных точках орбиты невелико. У Земли в ее движении вокруг Солнца е=0,016, υp=31км/с, υa=29км/с. У кометы Галлея эллипс очень вытянут: е=0,96; так что υp=51км/с, υa=1км/с. Такой характер ускорений и замедлений на орбите понять легко, если воспользоваться аналогией с вращением грузика на стержне вокруг горизонтальной оси. Внизу скорость наибольшая, наверху — наименьшая. В нашей задаче «вниз» — это направление к притягивающему центру, «вверх» — прочь от него. Причина изменений скорости и для планеты, и для маятника одна: закон сохранения энергии. «Наверху» потенциальная энергия гравитации максимальна, «внизу» — минимальна. Для кинетической энергии соотношение противоположно.
Набор орбит оказался небольшим. В век космонавтики мы можем выбирать высоту или период обращения искусственных небесных тел в широких пределах, но в силу (5) по отдельности, а не вместе. Наименьший период обращения ИСЗ — полтора часа — соответствует круговой орбите минимальной высоты. Максимального периода теоретически нет, но подавляющее большинство ИСЗ имеют период не более 24 час.
Многие искусственные спутники Земли (ИСЗ) летают низко, почти царапая Землю: в масштабе школьного глобуса (1:50000000) не далее сантиметра от него. Тут уж даже Землю шаром считать нельзя, хоть на глазок это и незаметно. А вот Юпитер и особенно Сатурн обладают отчетливо видимым сжатием. Одним словом, чтобы идти дальше, надо разобраться с формой небесных тел и их притяжением.
Начнем с последнего. Пусть нам известна форма и строение протяженного небесного тела Т. Как определить силу тяготения, с которой Т притягивает какую-либо частицу Q? Перейдем к ускорению — оно не зависит от массы пробной частицы (уникальное свойство гравитационного поля, открытое Г. Галилеем). Поэтому можно считать, что Т создает вокруг себя (и в себе самом тоже) поле ускорений, математически точное описание гравитационного поля. Как найти его? Разобьем мысленно Т на столь малые кубики, чтобы их размерами можно было бы пренебречь по сравнению с расстоянием до Q (рис.5).
Рис.5
Вектор ускорения ws, сообщаемого Q со стороны s-гo кубика, равен согласно (1)
ws=—(Gms/rs3)rs (6)
Поясним, откуда взялся минус и куб в знаменателе. Модуль ускорения равен Gms/rs2, и он умножен на единичный вектор —rs/rs направления от массы ms к точке Q (рис.5). Полное ускорение равно векторной сумме (6) по всем кубикам. Разумеется, так получается приближенная величина. Чтобы вычислить точную, нужно перейти к пределу, устремляя ребро кубика к нулю. В пределе получим тройной интеграл по телу Т. С помощью хорошего компьютера интеграл взять нетрудно. Но ведь даже для данного тела его нужно считать в огромном количестве точек пространства. Чаще всего идут другим путем. Как уже говорилось, Ньютон сумел вычислить интеграл для шара со сферическим распределением плотности и убедился, что внешние частицы шара притягивают в точности как материальная точка той же массы, помещенная в его центре. А дальше П.-С. Лаплас предложил следующую схему определения гравитационного поля Т. Во-первых, проще вместо векторного поля ускорений иметь дело со скалярным полем гравитационной потенциальной энергии Ер единицы массы Q. Оба поля однозначно определяют друг друга. Во-вторых, представим поле в виде ряда, т.е. суммы бесконечного числа слагаемых:
Ер=V0+V1+V2+… (7)
Здесь начальное слагаемое описывает притяжение шара с центром в центре масс Т и нам уже известно из формулы (4): V0=—К2/r. В отличие от силы, потенциал шара убывает обратно пропорционально первой степени расстояния от центра масс Т. Следующие слагаемые Vs убывают обратно пропорционально rs+1, причем V1=0. Если Q далеко, то достаточно взять несколько первых членов (7) или даже только начальный член, чтобы получить удовлетворительную точность. Иными словами, гравитационное поле любого тела с удалением от него все больше напоминает поле шара, в полном соответствии с наблюдением древних софистов, что издали и квадратная башня кажется круглой. Для близких Q (например, если Т — Земля, Q — ИСЗ) для высокоточного определения гравитации надо брать десятки и сотни слагаемых. Каждое из них представляет не очень сложную функцию координат точки Q. Например,
V2=(A1x2+A2y2—(A1+A2)z2+A3xy+A4yz+A5zx)/r5
Важно, что Vs содержит числовые коэффициенты. Например, в V2 их пять: A1÷А5. Эти коэффициенты можно определить, измеряя гравитационный потенциал, или ускорение на поверхности тела или вблизи нее. А можно следить за движением его искусственных спутников. В любом случае мы получаем систему многих алгебраических уравнений со многими неизвестными (коэффициентами типа