Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Гительман А.И. Динамика и управление судовых газотурбинных установок

.pdf
Скачиваний:
13
Добавлен:
25.10.2023
Размер:
15.08 Mб
Скачать

НосТи измерений и, во-вторых, о сохранении этого свойства при за­ носе проточной части.

Строго говоря, этот вывод относится только к положению изо­ дром, так как по указанным выше причинам подобный анализ в от­ ношении к. п. д. и положения границы помпажа представляет боль­ шие трудности. Однако практически одинаковое положение изодром в статике и динамике позволяет предположить, что и остальные параметры компрессора различаются мало.

Г Л А В А III ХАРАКТЕРИСТИКИ КАМЕР ГОРЕНИЯ,

ТЕПЛООБМЕННЫХ АППАРАТОВ, ТРАКТОВ, ГРЕБНЫХ ВИНТОВ

§ 11. Характеристики камер горения

При решении вопросов управления на установившихся режимах обычно достаточно использовать уравнение стационарного теплового баланса камеры горения без учета внешнего теплообмена:

Q;Br]K.r = c rTrGr - c BTBGB.

(143)

Здесь Q”. ккал/кг — количество тепла, выделяемое одним кило­ граммом топлива при полном сгорании (теплотворная способность топлива); В, кг/с— расход топлива в камере; т]к г — к. п. д. камеры (в современных камерах судовых ГТД т]к г = 0,98 н-0,99); Тг, Тв, °К — соответственно температура газа при выходе из камеры горения и температура воздуха при входе в камеру; сг, св, ккал/(кг-°С) — теплоемкости, соответствующие Тг и Тв; Gr, GB— секундный расход соответственно воздуха и газа через камеру го­ рения.

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

“ ■= -& •

044)

Здесь Lо, кг/кг — количество воздуха, минимально необходимое для полного сгорания одного килограмма топлива. Для нефтяных топлив Ь 0 = 14 -н15.

В камерах горения судовых газотурбинных двигателей коэффи­ циент избытка воздуха, как правило, на любых режимах не выхо­ дит за пределы ав = 3 н-30. Малые значения ав соответствуют ре­ жимам с высокой температурой газа, большие — режимам с низкой температурой [см. (146а)].

9 0

Учитывая, что Gr = GB+ В, можно написать, пользуясь выра­ жением для а:

 

Gr =

GB^.

(145)

где р, ~ 1 +

— коэффициент

увеличения расхода

газа в ка­

мере за счет сгорания топлива. При указанных выше значениях L0

и а в этот коэффициент может изменяться в пределах р = 1,002-И ,02. Следует отметить, что к. п. д. камеры горения при неизменных состоянии ее проточной части, качестве распыливания и сорте топ­ лива может быть представлен в виде зависимости от коэффициента избытка воздуха (рис. 46). Имея такую зависимость, нетрудно с по­ мощью формул (143), (144) представить характеристику камеры го­

рения в графической форме (рис. 47).

Рис. 46. Зависимость к. п. д. камеры горения ГТУ-20 от коэффициента избытка воздуха (экспериментальные данные).

 

Топливо — легкое

дизельное.

 

Приняв т)к. г =

0,98; Qp = 9800

ккал/кг; р =

1,0 и заменив

ст и св усредненной

величиной сср =

0,265, получим

из уравнений

(143), (144) формулы для грубой оценки перепада температур в ка­ мере горения:

 

Д < к . г = * г - * в ~ ю | - ;

(146)

 

ив

 

 

М к, г « 2400/ав,

(146а)

где В — в кг/ч,

GB— в кг/с.

основные ве­

Приведенные

соотношения позволяют определять

личины, характеризующие работу камеры горения на установив­ шихся режимах.

Работа камеры ограничена зоной устойчивого горения, за пре­ делами которой наступают низкочастотные колебания или срыв факела. Эту зону удобно выразить [10] в виде зависимости от рас­ хода воздуха, характеризующего скорости потока, и коэффициента избытка воздуха, характеризующего концентрацию топлива.

Однако, как правило, современные камеры судовых ГТУ могут работать без срыва факела в диапазоне, значительно превышающем значения а в, характерные для любых стационарных и переходных режимов ГТУ. Например, камера горения ГТУ-20 по данным стен­ довых испытаний может работать без срыва факела при а в > 50, т. е. при почти десятикратном обеднении смеси по сравнению с но­

91

минальным режимом. За десятилетний период эксплуатации (вклю­ чая заводские и ходовые испытания) срыва факела из-за обеднения или обогащения смеси в камерах ГТУ-20 не было отмечено ни на стационарных, ни на самых форсированных переходных и пусковых режимах, включая помпажи.

Причиной отдельных срывов факела при наладочных испытаниях служило падение давления топлива ниже давления газа в камере горения, т. е. прекращение (хотя и кратковременное) подачи топлива. Ограничение сброса топлива, вызываемое работой регуляторов скорости ТКВД, исключение кавитационного режима главного топливного насоса и другие мероприятия, устраняющие недопустимое

уменьшение расхода, т. е. давления

топлива, обеспечили надежную много­

летнюю эксплуатацию камер без срыва

факела. Здесь следует отметить, что

помимо

срыва

факела

ограничением

в работе камеры горения может явиться

чрезмерное увеличение

конуса факела,

возникающее при значительных коэф­

фициентах избытка воздуха, т. е. обыч­

но на очень малых режимах. Это мо­

жет привести к перегреву и прогарам

входных участков жаровых труб, к их

короблению. Так, на камерах ГТУ-20

в процессе наладочных испытаний было

замечено, что длительная работа (около

6 ч) на режиме ниже холостого хода

примерно на 20% по частоте вращения

Рис. 47. Расчетная характера ТКВД

при

а «а 2,5ав ноМ приводит

стика

камеры горения.

к прогару дефлекторов. Поэтому

при

режимах

следует

 

выборе программ

управления на малых

своевременно

учитывать

возможные

ограниче­

ния, налагаемые

указанным свойством камеры горения.

В

ряде

случаев

исходя

из

программ

управления

(например,

при

за­

пуске или на самых малых режимах) все же оказывается не­ обходимым обеспечить надежную работу камеры при коэффициенте избытка воздуха, превышающем допустимое максимальное значение. Одним из решений может служить уменьшение доли первичного воздуха на режимах повышенного а в с помощью перепуска части воздуха, поступающего в зону горения. Например, на начальном этапе создания двигателей ГТУ-20 в камере была предусмотрена специальная заслонка, которая до некоторого давления в цикле ГТУ оставалась открытой, увеличивая отношение расходов вторичного воздуха к расходу первичного. Впоследствии необходимость в таком перепуске отпала, так как все режимы, предусмотренные программой

управления, были отработаны в диапазоне а в

5-е12, обеспечиваю­

щем надежную работу камеры горения.

 

Рассматривая камеру горения в связи с вопросами управления, необходимо учитывать возможное изменение ее характеристик в про­

9 2

цессе эксплуатации. В частности, может изменяться неравномерность температурного поля за камерой, что следует учитывать при изме­ рении температуры газа в ГТУ (§ 31). Возможны также снижение к. п. д., сужение зоны устойчивой работы и другие изменения, свя­ занные с заносами, деформацией элементов проточной части камеры, ухудшением распыливания и смесеобразования. Например, при длительной работе может существенно увеличиться неравномерность подачи топлива (рис. 48). Особенно заметно эксплуатационные фак­ торы могут проявляться при работе на тяжелых сортах топлива.

о

т

т

т о eo6u;,«rA

Рис. 48. Изменение неравномерности комплекта форсунок.

• — в начале испытаний; X — после на­ работки 245 ч на моторном топливе ДТ-1; В общ — общий расход топлива через шесть

форсунок; # тах. £ т1п — максимальный и минимальный расход через одну форсунку.

Рис. 49. Осциллограмма давления перед форсунками ГТУ-20 при включении подачи топлива в процессе запуска (по данным испытаний на стенде).

Расчетные соотношения, описывающие работу камеры горения на пусковых и переходных режимах при плавном изменении рас­ хода топлива, в большинстве случаев могут быть приняты такими же, как на установившихся режимах, т. е. по формулам (143) и (146). Заметные отличия возникают обычно на переходных и пусковых режимах при резком изменении расхода топлива (это возмущение обычно является более сильным, чем возможные возмущения по расходу воздуха). В частности, для надежного зажигания камеры, как правило, целесообразно организовать начальный бросок дав­ ления перед форсунками (бросок расхода) по типу показанного на рис. 49. Это улучшает качество распыливания в момент подачи; ко­ нус распыливания становится шире и, как следствие, вероятность отказа в воспламенении рабочего факела от пускового блока умень­ шается. Однако для исключения заброса температуры газа сразу после завершения процесса воспламенения расход топлива необ­ ходимо быстро уменьшить. При этом время тВтах, в течение которого

допустимым является повышенный расход топлива, можно оценить соотношением

^ВОСПЛ < Т * ™ * < Т л шах >

(147)

гДе ^воспл = 0,1 -т-0,3 с — время, необходимое для завершения про­ цесса смесеобразования и воспламенения в современных камерах

9.3

горения; тлтах — допустимое время пребывания лопаток турбинной

группы при температуре газа,

соответствующей 5 тах.

Для определения

тлтах можно по соотношению (143) оценить

7 rBmax при Дтах (т. е.

полагая

режим камеры квазистационарным)

и далее по соотношению (306) вычислить тлшах для ТтВтт.

На форсированных режимах могут проявляться аккумуляция газовоздушной смеси в объеме камеры, аккумуляция тепла металлом жаровых труб и остальными омываемыми поверхностями, запазды­ вание, обусловленное временем, потребным для смесеобразования и горения топлива, запаздывание на участке форсунка—выходное сечение камеры горения и т. п. Влияние аккумуляции, оцененной по формулам, аналогичным (87) и (107), для камер судовых ГТУ в большинстве случаев весьма мало. В частности, для большинства переходных процессов ГТУ-20 отставание по температуре газа вслед­ ствие аккумуляции тепла металлом камеры не превышает 0,5%. Примерно такой же величиной характеризуется влияние объема. Таким образом, оба рассмотренных вида аккумуляции можно, как указывалось выше, не принимать во внимание в большинстве пере­ ходных режимов. Однако при исследовании отдельных участков процессов разгона и сброса, имеющих повышенные ускорения, при анализе режимов с малым запасом устойчивости и в других аналогич­ ных случаях влияние газовоздушной емкости и металла камеры, так же как и процессов смесеобразования, горения и др., следует учи­ тывать. Необходимые для этой оценки постоянные времени могут

быть определены по формулам Ту

к

=

ов

и 7 “ — - ! f K , где

 

 

tau

7V , Тх — постоянные времени объема и металла камеры; уср, GB—

средний удельный вес и расход газовоздушной смеси; VK, рм, см, F — объем камеры, вес, теплоемкость и поверхность металла, участвую­ щего в теплообмене; ам— коэффициент теплоотдачи.

Для камер горения главных судовых ГТУ значения постоянных времени на режиме полного хода обычно находятся в пределах Ту =

= 0,05 -г-0,1 и Г" = 5 ч-15 с. Для авиационных камер горения по данным В. А. Боднера и Ю. А. Рязанова TVk обычно не превышает

0,02 с. Анализ осциллограмм, приведенных в работе [6] для ряда камер стационарных ГТУ, показывает, что с учетом времени пере­ мещения топливного клапана (—0,05 с) время от начала заметного изменения температуры газа до выхода на уровень, примерно соот­ ветствующий новому установившемуся режиму, составляет 0,1— 0,2 с. Аппроксимируя рассматриваемый участок переходного про­ цесса экспонентой и предполагая, что на этом участке инерционность в основном обусловлена газовой и тепловой емкостью камер горения, можно оценить постоянную времени, соответствующую этой емкости, диапазоном 0,03—0,07 с.

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

94

Поэтому Для оценки этого времени наиболее надежно пользоваться экспериментальными данными, полученными на сходных по кон­ струкции камерах горения. В частности, интерес представляют дан­ ные, полученные для ряда камер горения, близких по конструкции к камерам судовых ГТУ. Например, по данным В. А. Боднера и Ю. А. Рязанова, передаточное запаздывание (суммарное время сме­ сеобразования и горения) для авиационных камер горения находится в диапазоне 0,1—0,2 с. Анализ осциллограмм, приведенных в ра­ боте [6], и экспериментальные данные В. И. Федорова показывают, что для камер горения ряда стационарных ГТУ чистое запаздывание составляет 0,1—0,3 с.

§12. Характеристики теплообменных аппаратов

Вкачестве гидравлического сопротивления и аккумулирующей воздушной емкости теплообменные аппараты (регенераторы и возду­ хоохладители) обычно рассматриваются в составе газовоздушных трактов (см. § 14). Поэтому здесь остановимся на их теплотехни­ ческих и теплоинерционных харак­ теристиках.

Регенератор. Основной характе­ ристикой регенератора на стацио­ нарных режимах является степень регенерации:

 

 

 

 

■и

(148)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

индексы

«в»,

«г», 1 и 2 — соот­

 

 

 

 

5, «г/с

ветственно означают воздух, газ,

 

 

 

 

 

 

 

 

 

вход и выход.

 

 

Рис.

50.

Степень регенерации г

Как

показывают расчеты и опыт­

в зависимости

от расхода

воздуха

ные

данные,

г с

приемлемой

точ­

GB

во всем

диапазоне

режимов

ГТУ-20

(стендовые испытания).

ностью

является

однопараметриче­

 

Кривая —

по формуле

(150).

ской

зависимостью от расхода

воз­

 

 

 

 

 

 

духа (газа) через регенератор. Например, на рис. 50 приведены такие данные для всего возможного диапазона режимов судовой ГТУ. При перепуске некоторой части газа в обвод регенератора снижение

степени регенерации

приближенно определяется зависимостью

 

 

 

'

тр

1

и

‘-'тр

) .

 

(149)

 

 

 

ТР

 

 

 

 

 

 

 

где

krp — постоянный

коэффициент;

GTp

и GB— количество

газа

(воздуха),

перепускаемое перед газовой

стороной регенератора,

и

расход

воздуха

по

его

воздушной

стороне.

При перепуске

в количестве GTp =

(0 -э-0,3)

GB коэффициент kTp

составляет

око­

ло 0,7.

 

 

 

 

 

 

 

 

 

95

Для приближенной оценки степени регенерации на стационарных режимах хорошее совпадение с опытными данными дает формула: 1

г

 

(150)

где индекс 0 определяет величины

на

расчетном режиме; kr *=» 0,2

для продольного обтекания, kr 0,3 для поперечного.

Таким образом, при известных г,

trl

и tbl из формулы (148) опре­

деляют температуру iB2, необходимую для расчета камеры горения и характеристик тракта.

Температура газа за регенератором tr2 приближенно определяется из условия равенства тепла, полученного воздухом и отданного га­ зом:

^ г С г ( ^ г 1

^гй)

6 в С в

(^в2 ^Bl)-

(151)

В большинстве случаев GBl «=> Gr

и св «

0,95сг. Поэтому

 

^Г2 ^ ^ri

0,95 (tB2

1вА)

(152)

Расчет теплоинерционных характеристик регенератора обычно представляет значительные трудности. Действительно, попытки представить регенератор в качестве одного звена с сосредоточенными параметрами приводят в отдельных случаях к недопустимым по­ грешностям. Аналитическое же решение уравнений в частных произ­ водных, учитывающих распределенность параметров, связано с боль­ шими трудностями вычислений и недостаточно наглядно. При этом обычно требуется значительная идеализация процессов. Несмотря на это, получаемые с помощью двумерного преобразования Лапласа передаточные функции описываются весьма сложными трансцен­ дентными зависимостями [22], [55]. С другой стороны, анализ и рас­ четы показывают, что эту сложную динамическую структуру можно заменить простой эквивалентной моделью, которая с приемлемой для практики точностью воспроизводит количественные соотношения между входными и выходными параметрами регенератора в динамике.

Представим уравнение (148) в виде

Л^в2 — k n ^ r l +

kr0Af,

(153)

где

В 10

1 Формула является упрощенным вариантом вывода [33].

96

— '

20

 

Ar =

Гр —г .

 

Гр

индекс 0 соответствует режиму, принимаемому за исходный. Такая форма позволяет рассматривать общее изменение темпера­

туры воздуха tB2 как сумму частных изменений, каждое из которых

вызывается только одним из входных параметров:

trl, tBl,

GB или

GTp/GB, т. е.

 

 

^ в 2 — | Д ^ в 2 |<r l = v a r " ф j A ^ B2j / B j= v a r ~ Ь | А ^ в2

| r = v a r .

( 1 5 4 )

Здесь и далее степень регенерации г рассматривается в виде ал­ гебраической функции расходов газа и воздуха как на стационарных, так и на переходных режимах.

Рассмотрим возможный характер реакции tB2 на изолированные возмущения по каждому из входных параметров.

При возмущении по trl передача тепла газом воздуху в условиях обычно выполняемого противотока наиболее эффективна в началь­ ный период, причем основную роль играет ограниченный (выходной) участок поверхности регенератора. Поэтому логично предположить, что реакцию tB2 в этом случае как при продольном, так и при по­ перечном обтекании можно приближенно описать уравнением одного инерционного звена. Запишем это условие в операторной форме:

A<n (s)

( 1 5 5 )

= AViL-1 T n s +

1

По тем же причинам уравнением инерционного звена можно огра­ ничиться и при описании реакции tB2 на изменение расхода газа G, т. е

 

1

Т

$

Т

Т

-

(1 5 6 )

Разбивая воздушную сторону регенератора на п последовательно

расположенных

участков, нетрудно убедиться, что реакция

tB2

на изменение tBl

эквивалентна

выходному сигналу, формируемому

на п последовательно включенных инерционных звеньях. Физи­ чески это связано с тем, что в начальный период после изменения tBl наиболее эффективный теплообмен происходит между воздухом и металлом на входных участках регенератора.

Как известно, п последовательно расположенных инерционных звеньев можно приближенно заменить последовательным включением звена чистого запаздывания и двух-трех инерционных звеньев. Учитывая, что реакция tB2 на возмущение по tBl является значи­ тельно более инерционным процессом, чем все остальные процессы, можно без особого ущерба для точности ограничиться последователь­ ным включением звена чистого запаздывания и одного инерционного

7 А. И. Гительман

97

звена, а в ряде случаев вообще описать процесс уравнением одного инерционного звена, т. е.

 

|д7 I.

Ь т -1

A^Bi (s) е

Ls

/11:74

 

I А Г в2 p B1= va r —

KB1L

T Bls +

1

 

ИЛИ

 

 

 

 

 

 

 

I

|<&1= v a r

=

1 7 ’BlSB1/

1

(1 5 8 )

где L — время

запаздывания

от начала возмущения по tBl

до на­

чала изменения

ts2-

 

 

 

 

 

Постоянные времени Тт1, Тв1, Тг и время чистого запаздывания L в уравнениях эквивалентных звеньев (155), (156) и (157) опреде­ ляются в основном такими величинами, как вес, поверхность, теп­ лоемкость теплообменных элементов и коэффициенты теплоотдачи. При этом последняя величина является существенно переменной, а ее изменение практически зависит только от весового расхода воздуха (газа). Поэтому Тг1, Тв1, Тг и L можно приближенно выразить в виде зависимостей от G. Подставляя (155), (156) и (157) в (154), получаем

At,

kriL-1

n (s)

+

kr0L

в2

 

Т’п* + 1

 

 

 

 

A r ( s )

,

T V s + i

" г

-l-M --1 ■-A^

1L •

(159)

Уравнение (159) позволяет при любых изменениях входных пара­ метров регенератора определить температуру воздуха за регенера­ тором как во времени, так и на новом установившемся режиме. Последний характеризуется условием s — 0, при котором уравнение (159) превращается в уравнение статики (154).

Для использования уравнения (159) в расчетах или в математи­ ческой модели ГТУ необходимо определить зависимость Тг1, Тг, Тв1 и L от расхода воздуха (газа). Для этой цели следует провести отдельный расчет регенератора как изолированного динамического объекта. Наиболее точным и физически наглядным следует признать широко распространенный в практике метод расчета, при котором регенератор разбивают на ряд последовательно или параллельно­ последовательно включенных участков, каждый из которых рассма­ тривается как звено с сосредоточенными параметрами. При расчете на ЭВМ число участков может быть принято таким, чтобы результат практически не отличался от схемы бесступенчатого распределения параметров (в большинстве случаев не требуется больше 5— 10 уча­ стков). Произведя во всем возможном диапазоне режимов серии ти­ повых возмущений по каждому из входных параметров, следует аппроксимировать полученные переходные процессы кривыми, яв­ ляющимися решением уравнений (155)—(157) (например, экспонен­ тами при ступенчатом возмущении). В результате обработки данных во всем диапазоне режимов получаются искомые зависимости Тг1, Тв1, Тг и L от расхода воздуха (газа), необходимые для упрощенного

98

о) til.

описания динамики регенератора в составе ГТУ с помощью урав­ нения (159).

Такой подход к динамическим характеристикам регенератора, по­ мимо сравнительной методической простоты, обеспечивает большую наглядность и отсутствие каких-либо трудностей при моделирова­ нии ГТУ на ЭВМ.

В частности, этот метод используется при математическом моде­ лировании характеристик регенераторов отечественных судовых ГТУ. Поучастковый расчет на ЭЦВМ, выполненный, например, для регенератора, имеющего пла­ стинчатые теплообменные элемен­ ты поверхностью около 1400 м3

Рис. 51. Расчетные динамические ха­

Рис. 52. Зависимость постоянных времени

рактеристики

регенератора: а — при

и

времени

запаздывания

регенератора

одновременном ступенчатом изменении

 

 

от расхода воздуха.

параметров;

б — при раздельном сту­

1,

2

— постоянные времени

при изменении

 

пенчатом изменении

температуры газа и соотношения расходов

1

температуры газа; 2 — соотношения

газа и воздуха; 3 — время запаздывания при

 

изменении

температуры воздуха на входе.

расходов газа и воздуха; 3 — температуры

 

 

 

 

 

воздуха на входе. Сплошные линии — ра*

 

 

 

 

 

счет на ЭВМ с учетом распределенности

с суммарным весом около 11 500 кг,

параметров,

штриховые — по формулам

(155) — (159);

А /в2 — отклонение темпе­

показал

следующее.

При ступен­

ратуры воздуха

от исходного значения.

чатом возмущении по

^г1 измене­

экспонентой

 

ние

tB2 хорошо аппроксимируется

(рис. 51, б). Постоянная

времени этого процесса зави­

сит

от расхода газа и воздуха

(рис. 52).

Такой же характер, но

с другой постоянной времени, наблюдается при ступенчатом

возмущении

по

соотношению воздуха

и газа (рис.

51, б). Само

изменение

расходов

при

неизменном

их соотношении весьма

слабо влияет на

iB2,

т. е.

в ряде случаев может не

учитываться.

Наконец, при ступенчатом изменении tBl

температура

tB2 некоторое

время остается почти неизменной (что может рассматриваться как чистое запаздывание) и только затем начинает изменяться по закону, который с достаточной точностью также можно аппроксимировать

7 *

99

Соседние файлы в папке книги из ГПНТБ