как давление в камере сгорания (КС) двигателя (тяга), массовые расходы окислителя и горючего на входе в двигатель при колебаниях давлений окислителя и горючего на входе в двигатель [1]. Обычно рассматриваемый частотный диапазон не превышает 50 Гц (низкие частоты). Из полученных АФЧХ ЖРД необходимо определить, в каких пределах может изменяться коэффициент усиления амплитудных характеристик.
АФЧХ двигателя получаются двумя способами:
1)экспериментальным, т.е. проведением специального стендового огневого испытания двигателя с пульсатором, который создаёт гармонические колебания давления на входе в двигатель по магистралям окислителя и горючего;
2)расчётным, т.е. разработкой математической модели ЖРД и созданием соответствующего программного обеспечения.
Рассмотрим однокамерный ЖРД с дожиганием генераторного газа при работе на основном режиме тяги (рис. 1).
На рис. 1 показана упрощённая схема взаимодействия агрегатов, гидравлических и газовых трактов, которые имеют следующие обозначения:
KR – критическая часть КС;
KC_вх– зоны смесеобразования и воспламенения продуктов сгорания в КС;
FF - газовые форсунки КС;
БНГ, БНО – бустерные насосы горючего и окислителя;
НО, НГ1, НГ2 – насосы окислителя, горючего I-ой и II-ой ступеней ТНА; РР – регулятор расхода горючего по линии газогенератора
(регулятор тяги двигателя); ГГ_вх – часть газогенератора после форсунок, где происходит основной процесс горения;
ГГ_вр – выравнивающая решётка газогенератора (ГГ); hm_g1 - hm_g12 – гидравлические магистрали горючего; hm_ok1 - hm_ok6 – гидравлические магистрали окислителя; gm_1-gm_5 – газовые магистрали.
Для проведения расчёта АФЧХ математическая модель двигателя, описывающая работу и физические процессы во всех агрегатах, гидравлических и газовых трактах, линеаризуется в окрестности рассматриваемого режима работы. В этой системе
85
уравнений коэффициенты при неизвестных зависят от конструктивных и режимных параметров. Полученную систему уравнений решают последовательно для каждого значения частоты колебаний из рассматриваемого частотного диапазона. Для наглядности результаты расчёта АФЧХ представляют в графическом виде (строят графики АЧХ и ФЧХ) [2, 3].
Рис. 1. Расчетная схема двигателя
Рассмотрим процесс получения уравнений для расчёта АФЧХ на примере вывода уравнения расхода через гидротурбину (1)
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Pвых) , |
|
|
(1) |
|||||
|
|
|
|
|
|
|
m F 2 (Pвх |
|
|
||||||||
где Pвх |
– давление на входе в турбину, |
|
|
|
|
|
|
|
|||||||||
Pвых |
– давление на выходе турбины, |
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
m – массовый расход жидкости через гидротурбину, |
|
|
|
||||||||||||||
F – условная площадь соплового аппарата гидротурбины, |
|
||||||||||||||||
– плотность рабочей жидкости. |
|
|
|
|
|
|
|
||||||||||
Продифференцируем по времени t уравнение (1), получим |
|
|
|||||||||||||||
|
|
|
1 |
1 |
|
|
|
dP |
dP |
|
|
||||||
|
|
dm |
|
|
|
|
|||||||||||
|
|
|
F |
|
|
|
|
|
|
2 |
вх |
|
вых |
|
. (2) |
||
|
|
|
|
|
|
|
|
|
|||||||||
|
|
dt |
2 |
|
|
2 (Pвх Pвых) |
|
dt |
dt |
|
|
||||||
86
Умножим обе части уравнения (2) на t и перейдём к записи уравнения в отклонениях, разделим обе части уравнения (2) наm , получим
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
m |
|
|
|
|
|
|
m |
2 P |
|
P |
|
Pвх Pвых . |
(3) |
||||
Введём коэффициент |
|
|
|
вх |
вых |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
t1 |
|
|
m |
|
. |
(4) |
||||
|
2 |
P |
P |
|
|||||||
|
|
|
|
|
вх |
|
вых |
|
|
|
|
Запишем переменные уравнения (3) в гармоническом виде через |
|||||||||||
cos t и sin t , получим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
m1 |
cos t m2 sin t t1 Pвх,1 cos t t1 Pвх,2 sin t |
|
|||||||||
|
t1 Pвых,1 cos t t1 Pвых,2 sin t . |
(5) |
|||||||||
Выделим переменные, стоящие в уравнении (5) передcos t |
и sin t |
||||||||||
. Получим два уравнения, которые войдут в систему уравнений |
|||||||||||
ЖРД для расчёта АФЧХ |
t1 |
Pвх,1 t1 Pвых,1 0, |
(6) |
||||||||
|
|
||||||||||
|
m1 |
||||||||||
|
|
t1 |
Pвх,2 t1 Pвых,2 0. |
(7) |
|||||||
|
m2 |
||||||||||
Для рассматриваемой схемы и режима работы ЖРД разработана математическая модель, включающая 396 линеаризованных алгебраических уравнений, и алгоритм программы на языке программирования MicrosoftVisualBasic для
MicrosoftExcel [5, 6].
Процесс расчёта АФЧХ в программе сводится к заданию исходных данных и параметров возмущений. Расчёт выполняется в программной среде MicrosoftExcel 2010. Полученные в программе результаты расчёта могут представляться как в табличном, так и в графическом видах.
Зададим возмущения в виде гармонических колебаний давления окислителя (магистраль hm_ok1) на входе в двигатель со значением амплитуды равным 1 кгс/см2. Получим следующие результаты расчёта АФЧХ в графическом виде (рис. 2 – 4). На рис. 5 представлены результаты расчёта АФЧХ при подаче на вход двигателя возмущения в виде гармонических колебаний давления горючего (магистраль hm_g1) со значением амплитуды равным 1 кгс/см2.
87
Графики на рис. 2 – 5 позволяют определить уровни частот и амплитуд, на которых повышается вероятность нарушения устойчивости и появления резонанса в магистралях и агрегатах двигателя.
По результатам расчётов проводится анализ, делаются соответствующие выводы. По предварительной оценке значение
коэффициента усиления амплитудных характеристик может находиться в пределах 1 – 3,89.
|
3.5 |
|
|
|
|
|
|
|
|
АЧХ |
0 |
|
|
|
|
|
|
|
|
|
|
|
-100 |
|
|
|
3 |
|
|
|
|
|
|
|
|
ФЧХ |
|
|
|
|
|
|
|
|
|
|
|
|
-200 |
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Амплитуда, кгс/см |
2.5 |
|
|
|
|
|
|
|
|
|
-300 |
Фаза, градусы |
2 |
|
|
|
|
|
|
|
|
|
-400 |
||
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
-500 |
||
1.5 |
|
|
|
|
|
|
|
|
|
-600 |
||
|
|
|
|
|
|
|
|
|
|
|||
1 |
|
|
|
|
|
|
|
|
|
-700 |
||
|
|
|
|
|
|
|
|
|
|
-800 |
||
|
0.5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
-900 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
-1000 |
|
|
0 |
5 |
10 |
15 |
20 |
25 |
30 |
35 |
40 |
45 |
50 |
|
|
|
|
|
|
Частота, Гц |
|
|
|
|
|
||
Рис.2. АФЧХ для давления в КС при колебаниях давления |
||||||||||||
окислителя
Амплитуда, кг/с
1.2 |
300 |
АЧХ
ФЧХ 200
1 



100
0.8 


0 0.6 




-100
0.4 |
|
|
|
|
|
|
|
|
|
-200 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
-300 |
0.2 |
|
|
|
|
|
|
|
|
|
-400 |
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
-500 |
0 |
5 |
10 |
15 |
20 |
25 |
30 |
35 |
40 |
45 |
50 |
Частота, Гц
Фаза, градусы
Рис. 3. АФЧХ для массового расхода горючего на входе в двигатель при колебаниях давления окислителя
88
|
4.5 |
АЧХ |
|
|
|
|
|
|
|
|
100 |
|
|
4 |
|
|
|
|
|
|
|
|
80 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
ФЧХ |
|
|
|
|
|
|
|
|
|
|
|
|
3.5 |
|
|
|
|
|
|
|
|
60 |
|
|
|
|
|
|
|
|
|
|
|
|
40 |
|
|
|
3 |
|
|
|
|
|
|
|
|
|
градусы |
|
|
|
|
|
|
|
|
|
|
|
20 |
||
кг/с |
|
|
|
|
|
|
|
|
|
|
||
2.5 |
|
|
|
|
|
|
|
|
|
0 |
||
2 |
|
|
|
|
|
|
|
|
|
-20 |
||
Амплитуда, |
|
|
|
|
|
|
|
|
|
Фаза, |
||
1.5 |
|
|
|
|
|
|
|
|
|
-40 |
||
|
|
|
|
|
|
|
|
|
-60 |
|||
1 |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
-80 |
|||
0.5 |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
-100 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
-120 |
|
|
0 |
5 |
10 |
15 |
20 |
25 |
30 |
35 |
40 |
45 |
50 |
|
|
|
|
|
|
Частота, Гц |
|
|
|
|
|
|
|
Рис. 4. АФЧХ для массового расхода окислителя на входе в |
||||||||||||
|
двигатель при колебаниях давления окислителя |
|
|||||||||||
|
0.35 |
|
|
|
|
|
|
|
|
АЧХ |
100 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0.3 |
|
|
|
|
|
|
|
|
ФЧХ |
50 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
2 |
0.25 |
|
|
|
|
|
|
|
|
|
0 |
градусы |
|
, кгc/см |
|
|
|
|
|
|
|
|
|
|
|||
0.2 |
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
-50 |
|||
|
0.15 |
|
|
|
|
|
|
|
|
|
|
|
|
Амплитуда |
0.1 |
|
|
|
|
|
|
|
|
|
-100 |
Фаза, |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
-150 |
|||
|
0.05 |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
0 |
|
|
|
|
|
|
|
|
|
-200 |
|
|
|
0 |
5 |
10 |
15 |
20 |
25 |
30 |
35 |
40 |
45 |
50 |
|
|
|
|
|
|
|
Частота, Гц |
|
|
|
|
|
|
|
Рис. 5. АФЧХ для давления в КС при колебаниях давления |
|||||||||||||
|
|
|
|
|
|
горючего |
|
|
|
|
|
||
Разработанная математическая модель, программное обеспечение и полученные результаты расчёта АФЧХ ЖРД используются для анализа продольной устойчивости ракеты, а также в САПР предметного назначения для повышения качества проектирования и доводки узлов, агрегатов и магистралей ЖРД и ракеты в целом.
89