Программа расчета параметров движения

ЦЕНТРА МАСС РН НА АТМОСФЕРНОМ

УЧАСТКЕ ПОЛЕТА

Программа обеспечивает моделирование движения центра масс РН в сферически-скоростной системе координат на участке полета первой ступени.

Параметрами управления на атмосферном участке полета (выделенные в массиве НУ курсивом) являются время выдачи предварительной команды на выключение маршевого ЖРД tкст1 и модуль максимального отрицательного угла атаки на участке разворота РН по тангажу αмах.

Идентификатор: RN1S

Используемые подпрограммы: AER1, LINTQ, ATMF

Программа использует файлы исходных данных IDRN1, AER1.DAT и файл результатов счета REZ1S.

Основной файл исходных данных IDRN1 включает следующие массивы:

К – настройки режимов работы программы, где:

К(1) – схема расчета траектории первой ступени:

0- пассивный участок полета отделившейся ступени не рассчитывается;

1-ПУТ ОЧ1 рассчитывается до момента Н<0;

К(2)-способ задания коэффициента силы лобового сопротивления отделившейся части первой ступени:

1-коэффициент Cx ОЧ не зависит от числа М и равен PD(14),

2-коэффициент переменный и равен Cx РН;

К(3) - 0 (в версии 12.00 не используется);

К(4) - 0 (в версии 12.00 не используется);

К(5) - 1 режим работы программы (в версии 12.00 не используется);

К(6) - режим печати результатов счета:

0-только в конечных точках,

1-в процессе полета с заданным шагом;

К(7) - 0 (в версии 12.00 не используется);

К(8) - режим выдачи результатов счета:

1-в массив REZ1S

2-на терминал и в массив REZ1S.

ST - параметры точки старта и модели Земли, где

ST(1) - широта (град) на общем земном эллипсоиде,

ST(2) - долгота (град),

ST(3) - высота над ОЗЭ (м).

ST(4) -2 (в версии 12.00 не используется).

F - начальные условия интегрирования, где

F(1) - время t (с);

F(2) - скорость V (м/с);

F(3) –тангаж (угол наклона вектора скорости к ЛМГ) (град);

F(4) –курсовой угол (от местной параллели). (град);

F(5) - высота Н. (м);

F(6) - широта (град);

F(7) - долгота (град).

R - краевые условия:

R(1) – 1(в версии 12.00 не используется);

R(2) – 1 (в версии 12.00 не используется).

PD - параметры конструкции РН и ДУ, где

PD(1) - стартовая масса РН (кг);

PD(2) - пустотная тяга маршевого ЖРД (кг);

PD(3) - пустотный удельный импульс маршевого ЖРД (с);

PD(4) - площадь среза сопла маршевого ЖРД (м2).

PD(5) - характерная площадь РН для расчета аэродинамики (м2).

PD(6)- коэффициент коррекции значения коэффициента силы лобового сопротивления, рассчитанного в подпрограмме AER1 (используется для оперативного изменения Сх без редактирования всего объема ИД по аэродинамике);

PD(7) - коэффициент коррекции значения коэффициента подъемной силы, рассчитанного в подпрограмме AER1;

PD(8)- предельное значение продольной перегрузки, выдерживаемой на участке между 8 и 9 точками циклограммы маршевого ЖРД;

PD(9) - коэффициент усиления автомата стабилизации при отработке nx1=const;

PD(10)- пустотная тяга рулевого ЖРД (кг);

PD(11)- пустотный удельный импульс рулевого ЖРД (с);

PD(12)- площадь среза сопла рулевого ЖРД (м2);

PD(13)- время выключения рулевого ЖРД и отделения ступени, относительно времени выключения (обнуления тяги) маршевого ЖРД (с);

PD(14)- коэффициент силы лобового сопротивления ОЧ 1 ступени, (используется при расчете баллистического ПУТ и не зависит от числа Маха);

PD(15)- минимально допустимый угол атаки (град);

PD(16)-максимально допустимый угол атаки (используется в АС для формирования программного угла атаки) (град);

PD(17)- 0 (в версии 12.00 не используется);

PD(18)- 0(в версии 12.00 не используется).

S-настройки блока интегрирования и вспомогательные параметры

S(1) - шаг интегрирования на активном участке траектории (с);

S(2) - шаг интегрирования в конце активного участка траектории (с);

S(3) - время перехода на мелкий шаг (отсчитывается от конца АУТ) (с);

S(4) - шаг интегрирования на пассивном участке траектории (с);

S(5) - шаг интегрирования в конце пассивного участка траектории (с);

S(6) - высота полета на нисходящей части ПУТ при переходе на мелкий шаг (м);

S(7) - шаг печати результатов счета на участке выведения (с);

S(8) - шаг печати результатов счета на пассивном участке полета (с);

S(9) - время окончания интегрирования ПУТ (с);

S(10) - время от момента старта, при котором производится сброс массы конструкции РН, например, головного обтекателя (м/с);

S(11) – величина массы отделяемой на АУТ (кг);

S(12) - масса отделяемой части первой ступени (для счета ПУТ) (кг).

D - настройки для решения краевой задачи, где

D(1) - 1 (в версии 12.00 не используется);

D(2) - 1(в версии 12.00 не используется);

D(3) - 1 (в версии 12.00 не используется);

D(4) - 1 (в версии 12.00 не используется);

D(5) - 1 (в версии 12.00 не используется);

D(6) - 0 (конечная высота интегрирования ПУТ ОЧ 1 ступени) (м).

UP – исходный массив параметров программы управления, где

UP(1) - функциональная команда «время начала перевода маршевого ЖРД на конечную ступень тяги (используется перед началом интегрирования для уточнения 9-12 точек циклограммы маршевого ЖРД) (с);

UP(2) - модуль величины минимального угла атаки (град);

UP(3) - время конца вертикального участка полета РН (с);

UP(4) – интервал времени от конца вертикального участка полета до момента достижения минимального угла атаки (с);

UP(5) - скорость изменения угла тангажа на вертикальном участке полета РН (град/с);

UP(6) – начальный угол атаки на вертикальном участке полета РН (град).

Массивы G(12) и G2(12) - (в версии 12.00 не используется).

RMD1, RMD1– исходная циклограмма работы маршевого ЖРД R(t): таблица времени (с) и таблица режима тяги (в долях номинала) Размерность массивов: (13).

RMD3, RMD4 -таблицы аргументов (с) и функции изменения относительной удельной тяги маршевого ЖРД (приведенной к номинальному режиму Rмд=1.0) при дросселировании двигателя Jуд(Rмд). Размерность массивов: (5).

Вспомогательный файл исходных данных AER1.dat используется для задания аэродинамических характеристик РН и включает следующие массивы размерностью (16):

Т11 – таблица аргументов по углу атаки (град) для зависимостей программа расчета параметров движения - student2.ru .

Т12 – таблица для задания зависимости коэффициента силы лобового сопротивления от угла атаки.

Т13 – таблица для задания зависимости коэффициента подъемной силы от угла атаки.

Т15 – таблица аргументов по числу Маха для зависимостей программа расчета параметров движения - student2.ru .

Т16 – таблица для задания зависимости приведенного коэффициента силы лобового сопротивления от числа Маха.

Т17 – таблица для задания зависимости приведенного коэффициента подъемной силы от числа Маха.

В процессе интегрирования уравнений движения РН формируются в массив результатов счета Y, где

Y(1) - время полета (с);

Y(2) – относительная скорость (м/с);

Y(3) - угол наклона вектора скорости к плоскости горизонта (град);

Y(4) - курсовой угол (град);

Y(5) - высота (км);

Y(6) - широта (град);

Y(7) - долгота (град);

Y(8) - масса РН (т);

Y(9) - число Маха;

Y(10) - угол атаки (град);

Y(11) - дальность от старта до проекции РН на поверхность ОЗЭ (продольная дальность) (км);

Y(12) - продольная перегрузка nx1;

Y(13) - скоростной напор (кг/м2);

Y(14) - тяга ДУ (т);

Y(15) - коэффициент силы лобового сопротивления Сх;

Y(16) - коэффициент подъемной силы Су.

Y(17) - режим работы маршевого ЖРД;

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

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

Наши рекомендации