Тяга Falcon 9 Merlin 1d рассчитана на каждый момент полета

В настоящее время я пытаюсь воссоздать запуск falcon 9 в единстве. Для расчетов мне нужна тяга первой ступени (делаю пока только для одной ступени).

Как рассчитать его в зависимости от высоты? Из значений у меня только от SpaceX (на уровне моря и в вакууме), но мне нужны значения в каждый момент полета. Какие-либо предложения?

Я делаю простую 2D-симуляцию, где ракета просто запускается прямо вверх.

космические значения

Смотрите здесь space.stackexchange.com/a/43845/6944 и ссылки в нем.
Но я имею в виду, как получить значения в любой момент полета? Для расчетов нужна колонка изменения тяги на протяжении всего полета. Любые предложения по этому поводу?
Обратите внимание, что двигатели не работают на полных оборотах в течение всего полета. Они регулируются в районе max-q, а также могут снижаться ближе к концу записи, чтобы уменьшить перегрузку.

Ответы (2)

Я собираюсь выразить это в английских единицах, потому что это то, что я чувствую, и поэтому у меня меньше шансов сделать неловкую ошибку.

Давайте воспользуемся уравнением тяги для невоздушных двигателей, скопированным отсюда Из Общего уравнения тяги в сторону Циолковского, как объяснить отбрасывание этих членов по ходу дела?

      Ф "=" м ˙ е В е + ( п е п 0 ) А е

Первое слагаемое справа от = — импульсная тяга . То, что мы здесь делаем, является постоянным для данной настройки дроссельной заслонки.

Вы даете:

  • Ф в = 1 849 500 фунтов силы
  • Ф с = 1 710 000 фунтов силы

Мы можем найти площадь выходной плоскости, подставив то, что мы знаем, в два уравнения и вычтя одно из другого. Это дает нам площадь выходной плоскости (для всех девяти двигателей) 66 футов 2 .

Ответ (без источника) на этот вопрос. Температура и давление выхлопа ракеты дают давление в выходной плоскости Мерлина, равное 0,7 атм.

Теперь мы можем рассчитать термин импульсной тяги м ˙ е В е быть 1 751 703 фунтов силы.

Быстрая проверка работоспособности: рассчитываемая вами импульсная тяга должна быть больше, чем тяга на уровне моря, и меньше, чем тяга в вакууме.

С этим и удобным столиком с атмосферой денди у нас есть все, что нам нужно.

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

Высота Давление внешней среды Толкать
0 футов 2117 фунт силы футов 2 1 709 793 фунта-силы
50 тыс. футов 241 фунт силы футов 2 1 833 609 фунтов силы
100 тыс. футов 23 фунт силы футов 2 1 847 997 фунтов силы
бесконечность и за ее пределами 0 фунт силы футов 2 1 849 515 фунтов силы

График зависимости тяги (фунт-сила) от высоты (фут)

введите описание изображения здесь

И, кстати, у меня есть еще один вопрос. Как вы получили площадь выходной плоскости для всех двигателей, равную 66 ft^2? Используя уравнение F (вакуум) - F (уровень моря) = p0 * Ae, я получил Ae, равный 6,119 м ^ 2, что намного меньше вашего Ae. Также, если мы посмотрим на цифры, мы увидим, что falcon 9 площадь поперечного сечения +- 42 м^2 - намного БОЛЬШЕ, чем у моего Ae, но меньше, чем у вас. Я что-то упустил в уравнении?
@mad.redhead помнишь, я использовал английские единицы измерения. 66 футов ^ 2 довольно близко к 6 м ^ 2.
но почему такая большая разница между Ae и площадью поперечного сечения??
@mad.redhead В этой статье в Википедии диаметр выхода Мерлина указан как 3 фута. Это чуть более 7 футов ^ 2, поэтому 9 из них составляют около 63 фута ^ 2. Чертовски близко к моим 66! en.wikipedia.org/wiki/SpaceX_Merlin
Оооо, это были не метры. Хорошо... я действительно запутал других людей и себя тоже. надо быть внимательнее(

Что касается зависимости тяги от высоты, то при решении аналогичной задачи я использую барометрическую формулу для получения давления окружающей среды на данной высоте. Вы можете получить разницу между известной тягой (или удельным импульсом, если хотите) в вакууме ( Th1 ) и на уровне моря ( Th0 ). Эта разница dTh = Th1 - Th0 будет умножена на зависящий от высоты коэффициент K , так что на уровне моря этот коэффициент будет равен 1, а в вакууме равен нулю. Затем из тяги (или удельного импульса) в вакууме вычитается произведение коэффициента на разность тяги: Тяга_на_заданной_высоте = Th1 - dTh * K.

Вам нужно получить этот коэффициент K . Воспользуемся для этого барометрической формулой: K = давление_рассчитанное_по_барометрической_формуле / давление_на_уровне_моря .

Я использую обе верхние формулы для давления в статье Википедии , в зависимости от базовой высоты, и таблицу для базовой высоты, поэтому зависимость тяги от высоты имеет ступенчатый характер.

Я не уверен, что это оптимальный путь, в моем случае он дает правдоподобные результаты (я использую программу GMAT для моделирования подъема ракеты). Но я все еще думаю о методической правильности этого пути.

Другие полезные ссылки: Плотность воздуха , Международная стандартная атмосфера .

Это немного отличается от того, что я ожидал рассчитать, но я тоже попробую! Спасибо за ответ, вы мне очень помогаете!