close

Вход

Забыли?

вход по аккаунту

?

Стендовая система имитации подачи гелия для наддува баков жидкостных ракетных двигателей..pdf

код для вставкиСкачать
УДК 681.518
СТЕНДОВАЯ СИСТЕМА ИМИТАЦИИ ПОДАЧИ ГЕЛИЯ
ДЛЯ НАДДУВА БАКОВ ЖИДКОСТНЫХ РАКЕТНЫХ ДВИГАТЕЛЕЙ
О.П. Свиридов, А.Ю. Матюхин, С.А. Курьянов, Г.И. Скоморохов
Представлены опыт отработки системы подачи в теплообменники охлажденного гелия для наддува баков на этапе доводочных испытаний двигателей на стендах Испытательного комплекса ОАО КБХА, методики настройки системы подачи на заданный режим и расчета ожидаемых изменений параметров гелия в процессе испытания
Ключевые слова: стенд, имитация, система подачи, гелий
Одной из важнейших задач при наземной
отработке реактивных двигательных установок
космических летательных аппаратов является
создание стендовых систем, имитирующих работу ракетных систем, таких как система подачи компонентов топлива, система наддува топливных баков. В предлагаемой статье рассмотрена стендовая система, имитирующая систему
наддува ракетных топливных баков.
Процесс отработки системы подачи гелия в наземных условиях
На этапе доводочных испытаний двигателей РД-0124 на испытательных стендах КБХА
была впервые проведена отработка системы
подачи в теплообменники (ТО) предварительно
охлажденного гелия для наддува топливных
баков горючего и окислителя.
Функционирование штатной системы подачи гелия заключается в следующем. Резервуары с гелием погружены в бак с жидким кислородом. Поэтому гелий до начала работы
системы находится в захоложенном состоянии
с начальной температурой, равной температуре
жидкого кислорода (~ 90К). Гелий подается на
вход в ТО по расходному трубопроводу, оснащенному двумя дроссельными устройствами,
расположенными параллельно и обеспечивающими массовый расход в заданном диапазоне
значений в течение всего времени работы двигателя.
В начальный момент времени подача гелия
осуществляется через одно дроссельное устройство, которое обеспечивает величину расхода по верхнему пределу заданного диапазона.
По мере опорожнения резервуаров и снижения
давления, происходит снижение расхода гелия.
Свиридов Олег Петрович – КБХА, канд. техн. наук,
нач. отдела, тел. (4732) 77-14-40
Матюхин Александр Юрьевич – КБХА, инженер,
тел. (4732) 49-04-00
Курьянов Сергей Александрович – КБХА, инженер,
аспирант, тел. (4732) 70-39-61
Скоморохов Геннадий Иванович – ВГТУ, д-р техн. наук,
профессор, тел. (4732) 43-76-38
При снижении расхода до нижнего уровня
заданного диапазона подключается второе
дроссельное устройство, что приводит к ступенчатому увеличению расхода.
Для проведения отработки теплообменников, разработанных и изготовленных в КБХА,
потребовалось создать стендовую систему подачи гелия наиболее полно имитирующую работу объектовой системы. Стендовая система
подачи гелия должна была удовлетворять следующим требованиям:
1.
Подача гелия к ТО во время огневого испытания двигательной установки на стенде производится из замкнутого объема.
2.
Начальная температура гелия на
входе в ТО в момент подачи и на режиме во
время испытания не должна превышать 93 К.
3.
Время работы системы подачи гелия
не менее 250 с.
4.
В процессе работы система подачи
гелия должна обеспечивать изменение массового расхода в заданных пределах (4÷80 г/с).
5.
При длительной работе системы подачи гелия (более100 с) необходимо подключение дополнительных дроссельных устройств
для поддержания расхода в заданных пределах.
В ходе разработки и проектирования была
принята и внедрена следующая система, схема
которой представлена на рис. 1.
Термостат предназначен для охлаждения
гелия до низких температур порядка 70-80 К,
набора начального давления и выдачи компонента к ТО. Конструкция термостата состоит из
баллона высокого давления, рубашки охлаждения и заборного устройства.
В процессе подготовки к испытанию газообразный гелий с температурой окружающей
среды закачивается в термостат с одновременным охлаждением жидким азотом, подаваемым
из ТРЖА в рубашку охлаждения. Из-за особенностей наземных испытаний, планировки стенда и невозможности размещения термостата в
огневом отсеке, он размещается на открытой
площадке рядом со стендом. При этом длина
расходного трубопровода подачи гелия составила порядка ~30 м. Для сравнения суммарная
длина штатных трубопроводов на РН составляет ~1,5 м. Из-за большой длины расходного
трубопровода к гелию поступало значительное
количество тепла из окружающей среды, что
создавало проблемы по обеспечению и поддержанию заданной температуры гелия на входе в ТО как в момент запуска, так и при работе
на режиме.
Рис. 1. Схема системы подачи гелия
1 – термостат; 2 - транспортный резервуар жидкого азота (ТРЖА);
3 - расходный трубопровод; 4 - система дроссельных устройств; 5 - трубопроводы
подключения к ТО; 6 - блок ТО; 7 - расходомерные участки
Поэтому было принято решение установить на трубопроводы подачи гелия рубашку
охлаждения и подавать в нее жидкий азот. Однако это привело к тому, что в определенный
момент жидкий азот из охладителя превращался в источник подогрева гелия, находящегося в
трубопроводе. В процессе адиабатического истечения гелия из термостата происходило ожидаемое снижение температуры, определяемое
по известным расчетным зависимостям. Но,
проходя по длинному трубопроводу, на входе в
дроссельные устройства гелий уже имел постоянную температуру, приблизительно равную
температуре жидкого азота, находящегося в
рубашке охлаждения. Из-за конструктивных
особенностей обвязки термостата и конфигура-
ции трубопроводов имелись дополнительные
участки с внешними теплопритоками как от
действия окружающей среды, так и от теплового воздействия работающего двигателя.
Все эти особенности сказывались на обеспечении заданного массового расхода гелия,
подаваемого в ТО. Настройка системы обеспечивалась комплектом дроссельных шайб, с
предварительно определенными характеристиками и набором в термостате расчетного начального давления в зависимости от задания на
испытание.
При проведении испытаний измерение
расхода гелия производилось критическими
расходомерами, в которых использовались
сверхзвуковые сопла с цилиндрической горло-
виной с погрешностью замера менее 1%. Расходомеры были установлены на выходе из ТО,
где гелий имел давление до 2,0 МПа (20
кгс/см2), а температуру от 400 до 550 К. Анализ
результатов отработки выявил расхождение
более 20% между ожидаемым расчетным расходом гелия и определенным с помощью расходомеров, что не позволяло обеспечить требования технического задания.
С учетом вышеизложенного, необходимо
было скорректировать и с достаточной точностью производить следующие расчеты:
- дроссельных устройств;
- начального давления гелия в термостате;
- текущих параметров гелия на всем времени
работы системы.
Для решения этих задач был разработан
математический аппарат, позволяющий проводить расчеты указанных параметров как функций от времени с учетом внешних теплопритоков и отклонений параметров реального газа от
идеального.
Математическая модель системы
Для вывода были использованы известные
формулы, описывающие изменение давления и
плотности газа в баллоне от времени:
dP k − 1 ⎡ dQ
dm
k
dV ⎤
, (1)
P
=
−h
−
⎢
dτ
V ⎣ dτ
dτ k − 1 dτ ⎥⎦
dρ
1 ⎛ dm
dV ⎞
=− ⎜
−ρ
⎟,
dτ
V ⎝ dτ
dτ ⎠
(2)
где - объем баллона, м3;
h - энтальпия газа, Дж/кг;
Q - количество тепла, подводящегося к
баллону, Дж;
m - масса газа в баллоне, кг;
Р - давление газа в баллоне, Па;
ρ - плотность газа в баллоне, кг/м3;
K - показатель адиабаты гелия.
Поскольку объем баллона постоянный,
следовательно, dV/dτ=0.
Поскольку гелий в баллоне перед испытанием охлажден до температуры жидкого азота,
то есть до предельной температуры, следовательно, dQ/dτ=0.
С учетом выше сказанного уравнения (1) и
(2) принимают вид:
dP
k − 1 dm
;
= −h
dτ
V dτ
dρ
1 dm
.
=−
dτ
V dτ
(3)
(4)
Уравнение газового состояния реального
газа имеет вид:
ρ=
Р
,
zRT
(5)
где z- коэффициент сжимаемости газа;
R- газовая постоянная, Дж/кг К.
Как видно из приведенных уравнений у реального газа четыре параметра (ρ, Р, z, T) для
определения которых имеются три уравнения
(3), (4), (5).
Чтобы система уравнений была замкнутая
необходимо записать четвертое уравнение.
Представим уравнение связывающее параметры Р, Т и z в виде:
z = a+b
P
.
T
(6)
Используя справочные данные [1], в диапазоне температур от 60 до 130 К и давлений от
4 до 40 МПа, при помощи метода наименьших
квадратов, получаем а=1,0; b = 1,423 ⋅10 −6
К/Па.
Тогда
bP = ( z − 1)T .
(7)
Погрешность определения по формуле (6)
одного из параметров по известным двум в
указанном диапазоне температур и давлений не
превышает 4%.
Подставив значение Т из формулы (7) в
формулу (5) получаем:
ρ=
z −1
zbR
ρ z 0 (z − 1)
.
=
ρ 0 z ( z0 − 1)
(8)
Из (8) видно, что плотность реального газа
однозначно определяется значениями коэффициента сжимаемости и газовой постоянной.
Определим остальные параметры газа в зависимости от коэффициента сжимаемости. Для
этого разделим (3) на (4):
dP
= (k − 1)h .
dρ
(9)
Для реального газа значение энтальпии, с
незначительной погрешностью, можно определить по формуле:
h=
kz
RT .
k −1
(10)
Из (9) и (10) получаем:
dP
= kzRT .
dρ
(11)
Продифференцировав формулы (7) и (8) по
z и поделив полученные выражения, друг на
друга получаем:
dP
dT ⎤
⎡
= z 2 R ⎢T + ( z − 1) ⎥ .
dρ
dz ⎦
⎣
(12)
Из формул (11) и (12) получаем дифференциальное уравнение, позволяющее определить
температуру как функцию от коэффициента
сжимаемости:
(k − z )T = z (z − 1) dT .
⎛ 2
С = k⎜
⎟
⎝ k +1⎠
где
m& =
Интегрируя формулу (13), получаем:
T =C
(z − 1)k −1 .
(14)
zk
Постоянную интегрирования определяем
из начальных условий, после чего окончательно получаем:
⎛z ⎞
T = T0 ⎜ 0 ⎟
⎝ z ⎠
k
⎛ z −1 ⎞
⎟⎟
⎜⎜
⎝ z0 − 1 ⎠
k −1
.
(15)
k
⎡ z (z − 1)⎤
P = P0 ⎢ 0
⎥ .
⎣ z ( z0 − 1)⎦
(16)
На основании формул (8),(15) и (16) получаем связь между параметрами Р, Т, ρ и z при
адиабатическом истечении гелия:
k
z ⎛ ρ ⎞
P ⎛ ρ ⎞
T
⎟
⎟⎟ ;
= 0 ⎜⎜
= ⎜⎜
P0 ⎝ ρ 0 ⎠ T0
z ⎝ ρ 0 ⎟⎠
k −1
⎞ k
z ⎛ P
T
= 0 ⎜⎜ ⎟⎟
T0
z ⎝ P0 ⎠
d
dτ
dm
= m& = µ ⋅ f ⋅ P k ⎜
⎟
dτ
⎝ k +1⎠
где:
А=
(18)
где µ - коэффициент расхода;
f площадь проходного сечения дроссельного устройства, м2;
zш - коэффициент сжимаемости газа на
входе в дроссельное устройство;
Тш - температура газа на входе в дроссельное
устройство, К.
Обозначим
f = µ ⋅ f ⋅C ,
(21)
(k − 1) f ⋅ z0 RT0 ;
V z ш RTш
bP0
bP
; zш = 1 + 0 .
T0
Tш
(22)
(23)
Из формул (17) и (21) получаем:
Р
= (1 + Аτ )
Р0
−
k
k −1
.
(24)
Из формул (19) и (24) получаем:
k
m&
= (1 + Aτ )− k −1 .
m& 0
(25)
Из формул (8) и (21) получаем:
z=
1
,
z ш RTш
k
⎛ ρ ⎞
⎟⎟ . (20)
⎜⎜
ρ
⎝ 0⎠
1
ρ
= (1 + Аτ )− к −1 .
ρ0
z0 = 1 +
Как видно из формулы (17) для реального
гелия, при адиабатическом истечении, зависимости температуры от давления и плотности
отличаются от зависимостей для идеального
гелия на коэффициент z0/z, который имеет переменную величину.
Расход реального газа определяется по
формуле:
k +1
⎛ 2 ⎞ k −1
⎛ ρ ⎞
f ⋅ z 0 RT0
⎟⎟ = −
⎜⎜
V z ш RTш
⎝ ρ0 ⎠
;
(17)
(19)
Интегрирование полученного уравнения
дает:
k −1
.
f ⋅P
.
zш RTш
Как отмечалось выше, температура газа на
входе в дроссельное устройство величина постоянная (Тш=const).
Зависимость плотности гелия от времени
определяется с помощью уравнений (4), (17),
(19). В результате несложных преобразований
получаем:
Воспользовавшись уравнением (7) получаем:
.
Окончательно получаем:
(13)
dz
k +1
⎞ k −1
1
z −1
1 − 0 (1 + Aτ )
z0
−
1
k −1
.
(26)
Полученные формулы позволяют производить расчет всех параметров газа при критическом истечении через дроссельное устройство.
Экспериментальные исследования
При подготовке системы к испытанию [3]
требуется рассчитать коэффициент расхода
дроссельного устройства f и величину начального давления на выходе из гелиевого баллона (Ро), которые необходимы для получения
заданного расхода газа. При этом, как правило,
( )
& max , m& min и время
задан диапазон расходов m
работы системы τ .
Расчет производится по следующему алгоритму:
1. Используя формулу (25) определяется
величина А:
А=
k −1
⎛ m& max ⎞ k
⎜⎜
⎟⎟
−1
⎝ m& min ⎠
τ
.
(27)
AV zш RTш
m&
; P0 = max
(k − 1) ⋅ z0 RT0
f
zш RTш . (28)
4. По формулам (23) определяются zo2 и
zш2 второго приближения.
5. Производится расчет, f , P0 и полученные величины сравниваются с величинами, полученными в первом приближении.
6. Расчет повторяется до тех пор, пока погрешность между двумя приближениями будет
меньше одного процента.
7. После этого из формулы f = µfC определяется диаметр дроссельного устройства:
d =2
f
µ ⋅π ⋅ C
.
(29)
При определении диаметра по этой формуле принимаются следующие значения входящих величин: µ=0,95; С=0,7252; π=3,14.С учетом этого:
d = 1,36 f .
(30)
8. После изготовления, дроссельное устройство продувается на аттестованном стенде и
определяется истинный коэффициент расхода
fи .
9. Производится расчет по алгоритму п. 26, в котором вместо f и и Ро определяются А
и Ро.
10. Производится проверочный расчет параметров гелия. Исходные данные для расчета:
f , Po , V , Tш , Тбо .
Алгоритм расчета следующий:
( )
( )
10.1. По формулам z0 = 1 +
zш = 1+
bP0
;
T0
bP0
. (23)
Tш
А=
10.2. По формуле
(k − 1) f ⋅ z0 RT0
V z ш RTш
(22) определяется параметр А;
Р
= (1 + Аτ )
Р0
По формуле
2. В первом приближении задаются значения: z01 = zш1 = 1 .
3. Используя формулы (19) и (22), рассчитываются величины:
f =
определяются коэффициенты сжимаемости;
−
k
k −1
(24)
определяется давление в конце испытания (τк);
По формуле
zш = 1+
bP0
Tш
(23)
определяется коэффициент сжимаемости zш . в
конце испытания;
Определяется средний коэффициент сжимаемости по формуле:
z
+z
zш.ср = ш.н ш.к .
2
По формуле
А=
(k − 1) ⋅ f ⋅ z0 RT0
V z ш RTш
(31)
(22)
определятся параметр А;
По формуле
Р
= (1 + Аτ )
Р0
−
k
k −1
(24)
определяется давление как функция времени;
По формуле
k
m&
= (1 + Aτ )− k −1
m& 0
(25)
определяется расход как функция времени;
По формуле
1
z=
z −1
1 − 0 (1 + Aτ )
z0
−
1
л −1
(26)
определяется коэффициент сжимаемости;
По формуле
⎛z ⎞
T = T0 ⎜ 0 ⎟
⎝ z ⎠
k
⎛ z −1 ⎞
⎜⎜
⎟⎟
−
z
1
⎝ 0
⎠
k −1
(15)
определяется температура гелия в баллоне;
По результатам расчета строятся графики
изменения этих параметров в зависимости от
времени. Проверка полученных результатов
производилась при огневых испытаниях двигателей на стенде ИК КБХА. По заданным величинам расходов и времени работы был проведен расчет по пунктам 1÷9.
Температура гелия в гелиевом баллоне
Тбо= 78 К.
Объем гелиевого баллона (термостата)
V=0,3 м3.
По этим исходным данным и по приведенным выше формулам был произведен расчет
параметров гелия. По результатам расчета построены расчетные графики (рис. 2), здесь же
приведены экспериментальные данные. Для
сравнения, произведен расчет параметров гелия
по формулам для идеального газа (при z0=zш=1)
(рис. 3).
Сравнительный анализ результатов расчета
и экспериментальных данных показал:
1. Погрешность расчета по приведенным
формулам менее 10%;
2. Погрешность расчета по формулам для
идеального газа порядка 15%.
В результате расчета были получены величины:
Для первого случая (τ = 65 с):
d1=2,2мм;
f1 = 2,456 ⋅ 10 −6 м2;
Ро=15 МПа.
Для второго случая (τ = 250 с с подключением второго дроссельного устройства)
d1=2,2мм; f1 = 2,456 ⋅ 10 −6 м2;
d2=1,4мм; f 2 = 1 ⋅10 −6 м2; Ро=17,85 МПа. При
подготовке к огневому испытанию в гелиевом
баллоне (термостате) были получены следующие температуры:
Температура гелия перед дроссельным
устройством Тш=83 К;
Температура гелия на выходе из гелиевого
баллона для первого случая (τ = 65 с) Т0=102 К,
для второго случая (τ = 250 с) Т0=95 К.
Давление, кгс/см2 ; Температура, К; Расход, г/с
180
Дав ление расчетное (Ррасч)
Дав ление экспериментальное (Рэ)
Р расч
160
Температура расчетная (Трасч)
Температура экспериментальная (Тэ)
140
Расход расчетный (mрасч)
120
Расход экспериментальный (mэ)
Рэ
100
Трас
Тэ
80
60
40
mрас
ч
mэ
20
0
0
20
40
60
80
100
120
140
160
180
200
220
240 260
Время, с
Рис. 2. Сравнение экспериментальных и расчетных параметров гелия (реальный газ)
с подключением второй дроссельной шайбы (dш1=2,2 мм, dш2=1,4 мм)
Отсюда можно сделать выводы
1. Полученные формулы для идеального и
реального гелия достаточно точно описывают
качественный характер изменения параметров
гелия во времени.
2. Формулы для реального гелия более
точно описывают количественный характер изменения параметров гелия во времени.
3. Разработанная методика позволяет производить определение параметров гелия, при
проведении огневых испытаний, с точностью
достаточной для инженерных расчетов.
Заключение
На этапе доводочных испытаний двигателей в КБХА проведена отработка систем подачи в теплообменники предварительно охлажденного гелия для наддува баков окислителя и
горючего в стендовых условиях. В процессе
отработки были решены следующие задачи:
- разработан математический аппарат, позволяющий производить расчет параметров
(давление, температура, массовый расход) гелия в процессе огневого испытания в стендовых условиях с точностью до 10%;
Давление, кгс/см2; Температура, К; Расход, г/с
180
Дав ление расчетное (Ррасч)
Дав ление экспериментальное (Рэ)
160
Температура расчетная (Трасч)
Ррасч
140
Температура экспериментальная (Тэ)
Расход расчетный (mрасч)
120
Рэ
Расход экспериментальный (mэ)
mрас
100
80
Тэ
60
mэ
40
Трас
20
0
0
20
40
60
80
100
120
140
160
180
200
220
240 260
Время, с
Рис. 3. Сравнение экспериментальных и расчетных параметров гелия (идеальный газ)
с подключением второй дроссельной шайбы (dш1=2,2 мм, dш2=1,4 мм)
- создана технология получения заданных
температур гелия в диапазоне 80÷90К путем
его предварительного захолаживания жидким
азотом в специально разработанном термостате;
- создано программное обеспечение, позволяющее проводить расчеты ожидаемых изменений параметров гелия при различных условиях работы системы подачи;
- обеспечены заданные температуры гелия
на входе в ТО в момент запуска и на режиме в
условиях большой удаленности (до 30 м) термостата при огневых испытаниях;
- обеспечены заданные массовые расходы
гелия в различных заданных диапазонах от 4 г/с
до 80 г/с;
- получена удовлетворительная сходимость
характера изменения входных параметров ге-
лия в условиях стенда и в штатных условиях
при работе системы подачи до 250 с.
В результате проведенной отработки и решения вышеизложенных задач была достигнута
достаточно полная имитация работы системы
подачи гелия в стендовых условиях штатным,
что было подтверждено прошедшими летными
испытаниями двигателя РД-0124.
Литература
1. Ривкин С.Л. «Термодинамические свойства газов»,
М., 1973. – 387 с.
2. Варгафтик Н.Б. «Справочник по теплофизическим
свойствам газов и жидкостей», М.,1972. – 496 с.
3. Клюев В.В. Технические средства диагностирования: Справочник / Под общей ред. В.В. Клюева. М.: Машиностроение, 1989. – 672 с.
Конструкторское бюро химавтоматики (г. Воронеж)
Воронежский государственный технический университет
BENCH SYSTEM OF IMITATION OF SUBMISSION
HELIUM FOR PRESSURIZATION OF TANKS OF LIQUID ROCKET ENGINES
O.P. Sviridov, A.J. Matjuhin, S.A. Kurjanov, G.I. Skomorohov
Are submitted experience of working off of system of submission in heatexchangeres the cooled helium for pressurization
of tanks at a stage of honing tests of engines at stands of the Test complex of KBkHA, a technique of adjustment of system of submission on the set mode and calculation of expected changes of parameters of helium during test
Key words: the stand, imitation, system of submission, helium
Документ
Категория
Без категории
Просмотров
8
Размер файла
378 Кб
Теги
баков, гелий, подачи, система, двигателей, pdf, стендовая, наддува, жидкостные, имитация, ракетный
1/--страниц
Пожаловаться на содержимое документа