Владимир Сурдин - Солнечная система (Астрономия и астрофизика)

Рис.7. Формы вращающихся тел. Указаны последовательности фигур равновесия несжимаемых, «жидких» тел (сплошные линии) и сжимаемых, газовых тел (пунктир). Оси вращения у всех фигур на рисунке расположены вертикально.

Небесные тела лунных и более размеров резко неоднородны: плотность в центре существенно превышает плотность у поверхности. Для Земли — на порядок, для Юпитера — на 4-5 порядков, для Солнца — на 7 порядков. Так что однородные фигуры равновесия служат лишь крайне упрощенными моделями. Но в случае медленного вращения форму поверхности можно представить аналогичным (7) рядом Ляпунова:

ƒ(φ)= R[ƒ 0(φ)+ƒ 1(φ)+ƒ 2(φ) +…] (8)

Тут требуются пояснения. Форму поверхности вращения естественно задавать уравнением r=ƒ(φ), связывающим широту φ с расстоянием от поверхности до центра масс r функциональной зависимостью ƒ . Таков смысл левой части (8). В правой части R — характерный размер тела, например, радиус равновеликого шара. Тогда ƒ 0 тождественно равна единице, так что в нулевом приближении тело является шаром r=R — const . Остальные члены ряда дают малые поправки, причем ƒ s пропорциональна q s . Здесь q=ω 2R 3/(GM) представляет собой безразмерный малый параметр, равный отношению центробежной силы к силе тяготения на экваторе шара массы М и радиуса R . Для Земли, Юпитера, Солнца q равно соответственно 0,0034; 0,083; 0,00002. Наибольшим значением q =0,139 в Солнечной системе обладает Сатурн.

Функция ƒ 1 имеет вид ƒ 1(φ)= Aq(1—3sin 2φ) , где число А определяется распределением масс внутри тела Т . Для однородного тела А=5/12 . Для противоположного крайнего случая сосредоточенной в центре массы, окруженной невесомой атмосферой, А=1/6 . Остальные ƒ s можно найти последовательно методом Ляпунова.

Функция ƒ , представляющая поверхность сжатого эллипсоида вращения Е , также может быть разложена в ряд (8), причем ƒ 0=1 . ƒ 1=е 2(1—3sin 2(φ))/6 , где е — эксцентриситет меридионального сечения. Подбирая его так, чтобы Aq=е 2/6 , добьемся совпадения Rƒ 0 и Rƒ 1 у Т и Е . Таким образом, любая фигура равновесия в нулевом приближении — шар, в первом — сжатый эллипсоид вращения.

Как рассчитывают трассы небесных тел в сложных гравитационных полях? Простых формул, подобных выведенным Кеплером и Ньютоном для описания движения частицы вокруг шара, для сложных полей не существует. Более того, за редчайшими исключениями вообще не существует абсолютно точных формул. Это следствие реальной сложности движений. Какими же средствами располагает современная наука? В самых общих чертах их можно разделить на две группы.

1. Аналитические методы.С их помощью сложное движение можно представить как наложение бесконечного числа простых движений. До предела упрощенный пример — знакомая по школьным учебникам формула суммы бесконечного числа членов геометрической прогрессии

1/x=1+(1-x)+(1-x) 2+(1-x) 3+… (9)

Предположим, что марсиане умеют складывать, вычитать и умножать числа, представленные десятичными дробями, и знают, что есть и обратное умножению действие — деление, но делить не научились. Так вот, левую часть (9) марсиане смогут вычислить, складывая большое количество чисел из правой части, а каждое из них получается умножением (1—х) самого на себя. Уже на этом простейшем примере видны две особенности аналитического подхода.

Во-первых, для получения точного ответа нужно проделать бесконечно много операций, что невозможно. Но для достижения заданной точности нужно произвести уже конечное число операций. Последнее тем больше, чем выше требования к точности. Пусть, например, х=1,1 . Чтобы ошибка (9) составила не больше 0,01 , следует взять два слагаемых справа; четыре слагаемых гарантируют погрешность менее 0,0001 .

Во-вторых, формулы аналитического метода работают не для одного какого-то набора значений переменных величин, а для любых их значений из некоторой области. Так, равенство (9) можно применять не только при х=1,1 , а для всех значений х от нуля до двух. Но для х=—1,5 формула (9) не годится и приходится прибегать к другим соотношениям. Напимер, для —2<���х<0 можно применить формулу 1/х=—[1+(1+х)+(1+х) 2+(1+х) 3+…] . В небесной механике орбиты разного типа также описываются разным набором аналитических формул.