![](/user_photo/_userpic.png)
книги из ГПНТБ / Шишков А.А. Газодинамика пороховых ракетных двигателей. Инженерные методы расчета
.pdfСравнение эффективности различных органов регулирования вектора тяги можно провести по их качеству, т. е. по отношению управляющей силы к потерям тяги Ry/AR, и по результатам ве сового анализа всей ракеты. В табл. 31 приведены соответству ющие данные по ряду органов управления для мощной транс портно-космической ракеты при одинаковом уровне потребных управляющих сил [23].
3. 10. ОБТЕКАНИЕ ВЫДВИЖНОГО ЩИТКА
При нагружении выдвижного щитка в сверхзвуковую струю в области среза сопла часть потока отклоняется от стенки соп ла и обтекает щиток.
Возьмем часть стенки плоского сверхзвукового сопла, на ко торой установлен щиток высотой /г. Поворот сверхзвукового по тока осуществляется в косом скачке, отходящем от линии отры ва. Рассмотрим течение в окрестности отрыва (рис. 42). Погра-
в)
Рис. 42. Обтекание сверхзвуковым потоком препятствия на стейке:
о —схема течения при обтекании |
выдвижного щитка; |
б—распределение давления: |
в- -взаимодействие |
потока со струйным |
препятствием |
ничный слой перед линией отрыва предполагается турбулентным, высота ступеньки (глубина погружения щитка) больше толщи ны пограничного слоя б. В косом скачке линии тока отклоняются от первоначального направления, параллельного стенке, на ко нечный угол влубь основного потока. Тангенциальный разрыв между отклоненным потоком и газом, находящимся у стенки, неустойчив и размывается в турбулентную область [31].
Из условия равновесия системы, состоящей из косого скачка и зоны турбулентного смешения, вычисляется перепад давления в скачке, отходящем от точки турбулентного отрыва сверхзвуко вого потока p2.lpi = f{Mi) ä; 1+0,5 Mi (см. § 3.5), где Mi — число Маха перед плоскостью отрыва хт.
Из хода рассуждений следует, что общая картина турбулент ного отрыва не зависит от того, каким способом он создается:
109
ступенькой (различного вида), противодавлением, падающей ударной волной, поперечным вдувом газа или впрыском жидко сти.
Экспериментальные данные показывают, что давление при отклонении сверхзвукового потока перед ступенькой повышает ся не скачком, а достаточно плавно: на участке длиной -—2,5 ö
перед |
плоскостью |
отрыва |
давление |
увеличивается |
до |
» 0 ,6 (р2— Рі) +Рі [36]. |
После |
плоскости |
отрыва, находящейся |
на расстоянии ~4,5 Іг от ступеньки, давление продолжает плав но расти и достигает значения р2 в плоскости, расположенной иа расстоянии — 2/г от ступеньки. Непосредственно перед ступень кой давление снова возрастает и на ее лицевой стороне дости гает величины » 1 ,6 р 2- Это объясняется натеканием на ступень ку высокоскоростной части клиновидной области турбулентного размыва тангенциального разрыва и образованием у верхней кромки отошедшей ударной волны у. Граница области размыва составляет со стенкой угол ~ 16,5°. Над этой линией располагает ся более или менее однородная область А сверхзвукового тече ния. Линия, соединяющая точку отрыва с верхней кромкой сту пеньки, составляет со стенкой ~13°. Эта граничная линия отде ляет поток, проходящий через ступеньку, от того потока, который попадает в зону обратных токов и образует клиновидную об ласть, поворачивающую поток в отходящем от ее вершины ко
сом скачке |
интенсивности |
p2jpі (р11с. 42, е). |
Боковая сила, |
вызванная |
турбулентным |
отрывом плоского |
сверхзвукового |
потока перед ступенькой, может быть записана в виде суммы
интегралов |
прироста давления по области |
перед |
отрывом |
|
(~ 2,5 б ) и по области отрыва ( — 4,5 /г). |
|
|
|
|
Приближенно |
|
|
|
|
Ry_ |
4,5 (рч— рi) hr (.МО |
|
|
|
R |
Р\Н |
п |
|
|
где R = piH/r(lAi) =РіН (1+kMr) — пустотная |
тяга |
двигателя |
||
|
при |
выведенном щитке; |
||
|
Н — высота поперечного сечения |
|||
|
плоского сопла шириной=1; |
|||
|
Кѵ= 2,3 Міг(Мі) — коэффициент относительной |
|||
|
боковой силы; |
|
||
|
Рщ = Іі/Н-— относительная площадь по |
|||
|
гружения щитка. |
|
||
Потери тяги вследствие погружения щитка в поток (или от |
||||
носительное |
сопротивление щитка) без |
учета |
давления на его |
тыльной стороне равны
дR __ |
1 ,%p2hr (М|) __ — |
R |
Р\Н |
где /Сі = 1,6р2/[/Оіг(М1)]= 1,6 г(Мі) (1+ 0,5 M i)— коэффициент относительных потерь.
по
Коэффициенты относительной боковой силы Кѵ и относитель ных потерь тяги К а , а также отношение RV/\AR\ =Ку/Ка д л я плоского сопла со щитком на срезе приведены в табл. 32 в за висимости от числа Мі потока (/г=1,2).
|
|
|
|
Таблица 32 |
Mj |
Т'д/Т'кр |
Ку |
Ка |
/ э д д я і |
2 |
1.9 |
0,80 |
0,55 |
1,5 |
3 |
6,4 |
0,59 |
0,34 |
1,7 |
4 |
28 |
0,46 |
0,24 |
1,9 |
Таким образом, относительная боковая сила и потери тяги при погружении плоского щитка_ пропорциональны относитель ной омываемой площади щитка Fm, а коэффициенты пропорцио нальности Кѵ и К а зависят от числа Мі и /г. Качество щитка, т. е. отношение Ry/\&R\, равно приблизительно 1,7.
Рис. 43. Схема течения перед цилиндром, установленным на пластине и облекаемым сверхзвуковым потоком:
а—вид |
сверху; |
б—пространственная |
картина |
течения; |
/—косой |
скачок |
уплотнения; 2—отошедшая ударная |
волна; |
Л—скачок уплотнения в отклоненном потоке; 4—граница застой ной зоны; 5—линия отрыва потока; 6—хвостовой скачок уплот нения; 7—застойная зона за цилиндром; 8—линия, проходящая через точку максимального давления в застойной зоне
Выдвижной щиток РДТТ имеет не только верхнюю, но и бо ковые кромки [41]. Турбулентный отрыв потока при натекании на такое трехмерное препятствие имеет более сложный харак тер. В области отрыва возникают течения, параллельные линии пересечения плоскостей щитка и сопла. Трехмерные эффекты становятся особенно существенными, если высота щитка превы шает размах в 5 и более раз (рис. 43). Например, в случае обте кания достаточно высокого цилиндра, установленного на пла стине, длина зоны отрыва (расстояние от передней кромки ци
111
линдра до линии отрыва) зависит только от диаметра цилиндра: /T» 2af«4r; максимальное повышение давления в отрывной зоне
Р-2ІРі ~ 1+0,5 Мі [17].
Рост расхода газов через область перед трехмерным щитком приводит к уменьшению длины зоны отрыва. Растекание потока вдоль лицевой стороны трехмерного препятствия приводит так же к изменению распределения давления в отрывной зоне вбли зи препятствия и на его лицевой стороне, где давление возра стает из-за натекания более скоростной части потока по сравне нию с двухмерным случаем. Необходимые уточнения боковой силы и сопротивления реального щитка определяются экспери
ментально. Опытные характеристики интерцепторной |
системы |
||
управления вектором тяги даны в табл. 33 |
[41]. |
|
|
|
|
Таблица 33 |
|
Относительная площадь щитка в выходном се |
0,05 |
0,10 |
0,15 |
чении сопла Fm |
|
|
|
Угол отклонения вектора тяги |
3° |
6° |
9° |
Относительные потери тяги 1R/R |
0,0'2 |
0,05 |
0,08 |
Обтекание дефлектора аналогнпчно рассмотренной картине течения в сопле около выдвижного щитка. Однако поверхность дефлектора составляет отличный от нуля угол с направлением оси у (перпендикулярной к оси сопла). Поэтому равнодейству ющая сил давления на дефлектор является не только сопротив лением, но и имеет составляющую по оси у. Кольцо дефлектора воспринимает ~2/3 полной боковой силы, а ~1/3 боковой силы приходится на стенку сопла в области отрыва потока.
3. 11. СОЗДАНИЕ БОКОВОЙ СИЛЫ ВТОРИЧНОЙ СТРУЕЙ
Если через отверстие в стенке сопла втекает в поток вторич ная струя жидкости или газа, то часть потока отклоняется от стенки, а вверх по течению от отверстия образуется зона повы шенного давления. При этом возникает боковая сила, которая складывается из реактивной силы вторичной струи и равнодейст вующей сил давления в области отрыва потока.
Обтекание струйного препятствия в сопле сопровождается не только возникновением боковой силы, но и приростом тяги, так как сопротивление вторичной струи не передается соплу, а равно действующая сил повышенного давления по зоне отрыва имеет осевую составляющую (боковая стейка сопла наклонена к оси РДТТ). Особое достоинство таких органов управления вектором
тяги заключается в том, что в данном случае исключается пли сводится к минимуму (внутри регулирующего клапана) взаимо действие между движущимися поверхностями и высокотемпера турным двухфазным потоком продуктов сгорания твердого топ лива.
Вторичная струя жидкости или газа воздействует на основ ной поток как источник массы, количества движения и энергии.
Рассмотрим поток в плоском сверх |
|
|
звуковом сопле (с выходным сече |
|
|
нием высотой Я), а также воздей |
|
|
ствие иа этот поток вторичной струм |
|
|
бесконечно малой интенсивности сЮ |
|
|
(рис. 44). Размеры возмущенной об |
|
|
ласти ограничены линией Маха, |
от |
|
ходящей от начала этой области под |
|
|
углом а и заканчивающейся на про |
|
|
тивоположной кромке выходного сре |
|
|
за. Длина области равна |
Рис. 44. Взаимодействие вто |
|
1 = Н ctg а = Н ctg arcsin — = |
ричной струн бесконечно ма |
|
лой |
интенсивности с пото |
|
М |
ком |
в плоском сверхзвуко- |
= н V м -- 1 .
Повышение давления в потоке вследствие ввода массы clG под углом ß к оси сопла оценим по теории одномерных пото ков [80, 93]:
|
dp |
_ |
А'М2 |
f dh |
[ 2 / j |
I |
к — 1 доз\ |
[___ |
|
|||
|
р |
М2— 1 I СрТ |
|
\ |
|
2 |
j a |
|
|
|||
|
|
- [ 1 + ( / г - |
1)Ма] |
^ |
C0S^ G |
N |
I , |
|
||||
|
|
1 |
1 |
|
; |
' |
|
vG |
I |
|
||
где /? |
и р ■— энтальпия и давление потока; |
|
|
|||||||||
ѵ = М.а — скорость потока; |
|
|
|
|
|
|
||||||
|
ÜBT — скорость вторичной струи; |
|
вторичной |
струи; |
||||||||
увт cos ß — осевая |
составляющая |
скорости |
||||||||||
|
N — молекулярный вес. |
|
по возмущенной области |
|||||||||
Равнодействующая |
сил давления |
|||||||||||
длиной |
l — H ] f |
— 1 равна |
dRyv = ldp, |
а |
соответствующий |
|||||||
удельный импульс определяется соотношением |
|
|
||||||||||
|
dRy р |
._ |
kpH М2 |
Gdh |
=*(■ |
4 м Й' - |
\ |
|||||
|
УР' |
dG |
G |
'М2— 1 |
I cpTdG |
|||||||
|
|
_ |
[ l - L |
( / ? _ |
1 ) M * |
Уat COS ß |
GdN I |
|
||||
|
|
|
|
|
|
|
V |
|
|
NdG I ' |
|
|
Результаты расчета Iy для случая |
вдува |
в основной поток (k= |
||||||||||
= 1,24; |
7о=3000 К; М = 2,5) |
различных газов и жидкостей при |
||||||||||
ведены в табл. 34 по данным работы [80]. |
|
|
|
ИЗ
|
|
Таблица 34 |
|
|
Тепло, поглощае |
|
|
|
мое впрыскиваемым |
|
|
|
веществом при на |
|
|
Вдуваемое рабочее тело |
греве от началь |
I,/, м/с |
ЛЯ*. |
ной температуры |
|||
|
до статической |
|
кг/м3■с |
|
температуры в по |
|
|
|
токе (Т= 1700 К), |
|
|
|
кДж/моль |
|
|
Инертные газы |
|
|
|
Н2 |
41,8 |
7300 |
0,06 |
n 2 |
46 |
2320 |
0,27 |
Продукты сгорания твердого топ |
—54 |
4260 |
7,1 |
лива (от= 1 , 6 6 г/см3) |
|
|
|
Нейтральные жидкости |
|
|
|
НоО |
144 |
1240 |
1 , 2 |
HCI2 F2 (фреон-12) |
156 |
1560 |
2 , 0 |
Вго |
83,6 |
1700 |
5,3 |
Реагирующие жидкости и газы |
|
|
|
(Н2 содержится в потоке) |
—426 |
6600 |
0,9 |
0 2=2Н оО—2Н2 |
—860 |
4700 |
6 , 8 |
N2О4 = N2“Ь4-Н2О—4-Н2 |
2 1 0 |
1470 |
2 , 2 |
N2 0 4 = 2 N0 2 ( д и с с о ц .) |
|
|
|
* 6 — плотность рабочего тела. Газы хранятся при давлении 10,3 МПа |
|||
и впрыскиваются со звуковой скоростью под углом ß=90° |
к оси |
сопла. |
Боковая сила при вдуве вторичной струи складывается из равнодействующей сил давления dRv р и собственной реактивной силы dRm = !BTdG.
|
riF |
Здесь |
/»т = г'втsin 3+ (Рв' — Р) du , |
где рвт— давление во вторичной струе; F-ач ■— площадь сопла вторичной струи.
В соответствии с этим удельный боковой импульс струйного органа управления равен (табл.35)
=/вТ'
Из результатов расчета Іѵ по изложенной линейной теории и сравнения с опытными данными (экстраполированными к зна чениям dG?xO) следует:
114
Таблица 35
|
Газовая постоян |
Удельный импульс |
Коэффициент |
|
Вдуваемый газ |
ная, |
усиления |
||
Дж |
вторичной струп |
|
|
|
|
|
/цт> М/С |
к |
І у |
|
кг ■К |
К |
у - j |
|
|
|
|
1ВТ |
|
н2 |
4160 |
2000 |
|
3,7 |
n 2 |
295 |
500 |
|
4,7 |
Продукты сгора |
314 |
1800 |
|
2,4 |
ния твердого топ лива
—линейная теория дает удовлетворительное описание зави симости боковой силы от газотермохимических факторов;
—удельный боковой импульс, возникающий при вдуве рабо чего тела в закритическую часть сопла, в несколько раз выше собственного удельного импульса вторичной струи /вт; коэф
фициент усиления для газовой вторичной струи равен
(dG— М3)
=2,44-4,7;
'ПТ
—впрыск жидкостей, реагирующих с выхлопными газами, и
вдув горячих продуктов сгорания твердого топлива (из каме ры РДТТ) — способы, обладающие наибольшими потенциальны ми возможностямии.
Картина отклонения части сверхзвукового потока вторичной струей аналогична картине обтекания препятствия, установлен ного на стенке (см. рис. 42, в). Перед проходящей по всему пло скому соплу поперечной щелью, через которую вдувается вто ричная газовая струя, образуется отрывная зона с возвратными течениями. Повышение давления в косом скачке, отходящем от вершины клинообразной отрывной зоны, приближенно равно РгІРі = 1+ 0,5Mi, где Mi — число Маха перед линией отрыва [77].
Длина отрывной зоны, как и прежде, пропорциональна высо те проникновения вторичной струи: Ita (4,2-=-4,5)/гСТр- Высоту проникновения /гстр можно определить из равенства количества движения струи в поперечном направлении ее сопротивлению основному потоку, т. е.
°вт'г'иг= [(!+ ?) Р г ~ аРЛ Астр.
где (l+ ß )p 2 — среднее давление на переднюю грань контроль ной поверхности, ограничивающей вторичную струю, причем
і + ( й — і)М?/4
115
api — среднее давление на заднюю |
грань ( а ^ І ) . |
|
||
Таким образом, для суммарной боковой силы при вдуве газа |
||||
из поперечной щели в основной поток получаем |
|
|||
Ry = Ry и -г |
= 4,2 ( /;2 — Л ) Лстр -f А>„, |
2 , 2 ДА х О ң ' У ң г _____ |
+ ^вт> |
|
ß)(l 4- ÜjöMj)—а |
||||
|
|
|
где Rur — тяга вторичной струп.
Так как в пустоте (RnT)cc = k+\/k[GaKPz(K)]ßT, для коэффици ента усиления в пустоте KU= RUI (Rin)** имеем
д- = |
______________ 2,2Л']/сЛпт_______________ I |
/г„т |
|
у |
(А- + 1) [(1 + Р) (1 -I- 0 , 5 Д \ ,) - а] г (Х„т) |
(/?,„)„ |
' |
Расчет коэффициента усиления по этой формуле |
при вдуве |
||
азота в воздушный поток дает значения /\?у= 2,94-3,2 |
при Мі = |
=2-^6, согласующиеся с данными специальных опытов.
Вреальных соплах РДТТ вдув вторичной струи осуществля ется через отверстия в щели ограниченной длины. При этом воз никают течения в направлении, параллельном щели, и вторичная струя обтекается не только сверху, но п с боков. Картина тече ния с учетом трехмерных эффектов существенно усложняется. Ведущая роль в определении газодинамических характеристик струйных органов управления вектором тяги РДТТ принадлежит эксперименту.
Результаты экспериментальных исследований взаимодействия воздушного потока с поперечной газовой струей, вытекающей из отверстия, проанализированы в работе [19] с позиций теории
подобия. При этом предполагалось, что влияние кривизны сте нок мало, а зоны перераспределения давления, обуславливаю щие возникновение дополнительной боковой силы, полностью умещаются на стенке сопла. В широкой области значений па раметра Gmvj(pL2) удалось провести корреляцию опытных зна чений коэффициента усиления:
К у= |
- |
1,87 lg- ^ + 1 , 1 2 , |
|
|
|
p L °- |
|
где V и р — скорость и давление основного потока |
(в месте |
||
вдува); |
от |
критического сечения до |
центра от |
L — расстояние |
верстия вдува.
В исследованном диапазоне параметра GnTv/(pL2) =0,06-1-0,9 коэффициент усиления изменяется от 3 до 1.
Влияние трехмерных и нелинейных эффектов приводит к то му, что удельный боковой импульс / у и коэффициент усиления Ку уменьшаются при увеличении относительного расхода вторич ной струи GBT/G (табл. 37). Относительная боковая сила Ry/R сначала растет с увеличением GnT/G, достигает максимума при некотором его значении, а затем уменьшается при дальнейшем
116
Таблица 36
|
РДТТ; G = |
|
|
|
|
|
|
РДТТ ракеты |
|
|
|||
|
=•■16,7 кг/с; |
ЖРД; |
G =7,25 |
кг/с; |
|
|
|
||||||
|
|
|
„Титан-ЗС“; |
|
|
||||||||
|
вдув из ПГГ |
вдув из ЖГГ |
[40] |
|
|
|
|
||||||
|
|
впрыск N0O4 [64] |
|
||||||||||
|
[40] |
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
||
Gbt/G |
0 , 0 1 2 |
0,03 |
0,03 0,05 |
0,07 |
0,09 0,18 |
0 , 0 2 |
0,05 |
0 , 1 0 |
0,15 0 , 2 0 |
0,30 |
|||
Ry/R |
|
|
|
|
|
|
|
1 |
|
|
|
0,115 0,14 |
|
0,023 |
0,05 |
0,05 0,065 0,077 0,07 0,045|0,02 0,0-15 0,075 0 , 1 0 |
|||||||||||
ly, м/с |
4300 |
3600 |
3000 2400 |
2 1 0 0 |
1500 |
470 |
2600 2400 |
2 0 0 0 |
1700 1600 |
1 2 0 0 |
|||
Ку |
— |
— |
2 1 , 6 |
1,4 |
0,9 |
0,3 |
— |
— |
— |
— |
— |
— |
|
росте |
GBT/G. |
Например, |
для |
рассмотренного |
в табл. |
36 ЖРД |
(Ry/R)max= 0,077 при GBT/G = 0,07. Немонотонный характер зави
симости R,j/R= f(GrsT/G) объясняется тем, |
|
|
|
|
|
||||||
что при слишком большой |
интенсивности |
|
|
|
|
|
|||||
вдува возмущенная |
область захватывает |
|
|
|
|
|
|||||
газовые слои вблизи противоположной по |
|
|
|
|
|
||||||
верхности сопла, и давление повышается |
|
|
|
|
|
||||||
на таких участках сопла, которые дают |
|
|
|
|
|
||||||
боковое усилие в направлении вторичной |
|
|
|
|
|
||||||
струн. Наибольшее значение боковой силы |
|
|
|
|
|
||||||
получается |
тогда, |
когда |
возмущенная |
Рис. |
45. |
Схема |
тече |
||||
область занимает приблизительно |
поло |
мня |
при |
наибольшем |
|||||||
вину проходного сечения сопла |
(рис. 45). |
допустимом |
интен |
||||||||
На рабочем участке для зависимости |
сивности |
вдува: |
|||||||||
Ry/R = f (GBT/G) может быть использована |
/ —у д а р н а я |
в о л н а ; |
2— |
||||||||
линейная либо степенная |
аппроксимация. |
в о з м у щ е н н а я |
об л а с т ь ; |
||||||||
|
5—м есто |
в д у ва |
|||||||||
В частности, при вдуве фреона Rv/R tt 0,35, |
|
|
|
|
|
||||||
GBT/G |
[90]; |
Ry/R=0,3 (GBT/G)0’5' |
для |
ЖРД |
[40]; |
R,JR = |
|||||
= 0,4 |
( G b t / G ) |
0 ' 7 2 для РДТТ ракеты «Титан-ЗС» [64]. |
|
|
|
||||||
Высокая эффективность системы впрыска в ракете «Титан-ЗС» обеспечи |
|||||||||||
вается |
конструкцией форсунок, |
порядком и местом |
расположения |
их |
в сопле |
||||||
( 4 н р = |
0,96 м, |
(/„ = 2,71 |
м), а |
также |
подбором |
впрыскиваемой |
жидкости. |
||||
Впрыск N2O.1 осуществляется в каждом квадранте через 6 |
форсунок под дав |
лением 5,25 МПа, диаметр проходного сечения в месте впрыска 1,8 м. При
полете непрерывно определяется избыток N20.i, который сливается через |
все |
24 форсунки, не создавая боковой силы. Масса NjO/, рассматривается |
как |
прибавка к топливу с собственным импульсом 100=120 [64]. |
|
Г л а в а IV
НЕУСТАНОВИВШЕЕСЯ ТЕЧЕНИЕ ГАЗА В РДТТ
Неустановившиеся течения газа осуществляются в пороховом двигателе в периоды воспламенения заряда, отсечки тяги и по летного регулирования тяги. Задачи об изменении характеристик РДТТ при неустановившнхся режимах работы являются нели нейными даже в квазистационарном приближении.- Для их реше ния применяются численные методы, а также приближенные, обладающие различной степенью сложности и точности [9, 62, 65, 69, 75]. Полученные аналитические решения для переходных про цессов в двигателе при выходе на режим, отсечке тяги и перете кании газов по соединительному газопроводу сообщающихся двигателей просты и достаточно точны. Время неустаповившегося истечения в рассматриваемых случаях порядка ІО-2 времени работы РДТТ. Изменение геометрических характеристик двига теля (свободного объема и площади проходного сечения) за это время того же порядка малости (10~2) и в дальнейшем не учи тывается.
4.1. ПЕРЕХОДНЫЕ ПРОЦЕССЫ В РДТТ ПРИ ВОСПЛАМЕНЕНИИ ЗАРЯДА
Во время горения твердого топлива со скоростью и в прогре том слое устанавливается распределение температуры, прибли женно описываемое экспоненциальной зависимостью
Т ( x ) ^ T 3Jr (Fs~ T 3)е- д'"/а,
где Ts, Т3— температура поверхности горящего топлива и на чальная температура заряда;
ü = X / ( c q ) — коэффициент температуропроводности; X— расстояние от горящей поверхности.
Всего в прогретом слое аккумулировано количество тепла
со
cpQ ( T - T 3) d x ^ ± ( T s- T 3).
о
Основной запас этого тепла заключен в слое толщиной б= а/и, время прогрева которого порядка 4 = б/« = а/и2 (времена тепло
118