Я моделирую систему из двух двойных звезд, и мне интересно, так ли это.
То, как я это делаю сейчас, заключается в том, что я сначала вычисляю центр масс, а затем радиус от каждой планеты до центра масс.
Затем я вычисляю ускорение для первой планеты с помощью:
где и , вместе составляют полное расстояние между планетами.
Затем я вычисляю скорость
То же самое можно сказать и о другой планете:
где и , вместе составляют полное расстояние между планетами.
Затем я вычисляю скорость
Расстояние по всему кругу (путь, который они пройдут, если они будут двигаться по кругу) равно и
Затем я вычисляю время, необходимое для поворота одного периода, путем деления расстояния на скорость, а затем вращаю планеты вокруг радиуса (от планеты до центра масс) таким образом.
Теперь я прочитал, что, поскольку их массы неравны, они не должны вращаться по кругу. Я сделал это неправильно?
Вот как это выглядит на данный момент:
Пересмотренный ответ в ответ на ваши комментарии:
Ваши расчеты верны.
Обычно орбиты будут эллиптическими, но каждая из двух планет может вращаться вокруг ЦМ по кругу с постоянной скоростью по кругу. (Ускорение — это скорее изменение направления, чем изменение скорости.) Если массы равны, они будут двигаться по одной и той же окружности; если они разные, окружности будут концентрическими. В любом случае они всегда находятся на противоположных сторонах ЦМ и поэтому имеют одинаковый период.
Ваша анимация создает впечатление, что планеты имеют одинаковый размер (и, следовательно, массу). Если да, то их орбиты неверны, потому что КМ вращается.
Сила, в вашем случае вы записали уравнение силы для ньютоновской гравитации , не то же самое, что ускорение. Так что либо приравнивайте силы, и в этом случае вам придется использовать центростремительную силу , а не ускорение, либо приравнивайте ускорения, и в этом случае вы можете использовать второй закон Ньютона, чтобы найти ускорение из силы, полученной с помощью ньютоновской гравитации. Оба метода приведут к одному и тому же ответу.
И вообще двойные звезды будут вращаться друг вокруг друга по эллиптическим орбитам, но круговая орбита является частным случаем эллиптической орбиты, а именно с нулевым эксцентриситетом . Поскольку для эллиптических орбит в целом нет решения в замкнутой форме для положения как функции времени (это потребует от вас итеративного решения уравнения Кеплера ), поэтому в вашем случае может быть проще просто придерживаться круговых орбит.
Qмеханик
AccidentalTaylorРасширение
Дэвид Лунд
Сэмми Песчанка
Дэвид Лунд
ПрофРоб