Материал: Ракетно-космическая техника и технология 2015. тез. докл. Российской науч.- техн. конф., посвященной 50-летию первого запуска РН «Протон» 16 июля 1965 года. Федеральное космическое агентство

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

как давление в камере сгорания (КС) двигателя (тяга), массовые расходы окислителя и горючего на входе в двигатель при колебаниях давлений окислителя и горючего на входе в двигатель [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

Источник: https://studfile.net/preview/16568800/