
книги из ГПНТБ / Переходные процессы в газотурбинных установках
..pdfосновными его параметрами. Однако полученную детальным расчетом или экспериментально и представленную графически характеристику можно аппроксимировать простыми аналитическими зависимостями с некоторыми коэффициентами. Значения коэффициентов берутся в каждом случае по конкретной характеристике, которая пред полагается заданной. Уравнения, полученные при такой аппрокси мации, по форме универсальны, т. е. одинаковы для любого компрес
сора, и |
отличаются |
лишь |
постоянными. |
Ниже |
предлагается |
простой |
метод представления характеристик |
компрессора аналитическими зависимостями, удобными для расчетов
на ЦВМ. Отметим, |
что |
для |
|
|
|
|||
использования на АВМ полу |
|
|
|
|||||
ченные ниже зависимости от |
|
|
|
|||||
носительно сложны, так как |
|
|
|
|||||
их реализация требует |
боль |
|
|
|
||||
шого |
числа нелинейных |
ре |
|
|
|
|||
шающих |
элементов. |
|
|
|
|
|
||
В |
большинстве |
случаев |
|
|
|
|||
достаточную точность обеспе |
|
|
|
|||||
чивает аппроксимация |
|
изо- |
|
|
|
|||
дром на характеристике |
ком |
|
|
|
||||
прессора |
уравнением |
|
|
|
|
|
||
8—1 |
с [ ' - ( ^ - ) Ъ |
Рис. II.5. |
Схема аппроксимации |
изодром |
||||
компрессора |
|
|||||||
|
|
|
|
(11.33) |
|
G—относительный |
|
|
где |
є — степень повышения |
давления; |
(при |
|||||
веденный) |
расход; |
А, |
В |
я С—постоянные |
для каждой изодромы |
|||
величины |
(рис. II.5). |
|
|
|
|
|
Показатель степени е принимается равным двум или трем, в ред ких случаях четырем (большие значения для более изогнутых изо дром).
Для двух крайних значений относительного |
(приведенного) |
||||||
числа оборотов пг и п г |
коэффициенты Аг, Вг, Сх, |
А2, |
В2, С 2 |
берутся |
|||
по характеристике (рис. II.5) |
соответственно в |
масштабе |
G и є. |
||||
Для любого другого числа оборотов п коэффициенты |
определяются |
||||||
по интерполяционным |
зависимостям: |
|
|
|
|
||
Л = |
Л1 |
+ ( = ^ Г ( Л , - Л 1 ) ; |
|
|
(11.34) |
||
в = |
вг |
П — И х |
(В,-В,); |
|
|
(11.35) |
|
|
|
|
|
||||
|
|
«2 — " l |
|
|
|
|
|
с = |
с х |
+ п - |
•пх |
(С2-Сх). |
|
|
(11.36) |
|
|
« а |
|
|
|
|
|
Здесь показатель степени |
т также легко определить по заданной |
характеристике, записав, например, уравнение (11.35) для оборотов пх
61
и nt. Обычно можно принять т. = 1,25. Отметим, что величину т иногда целесообразно округлять до значения, более удобного для
расчетов. Так, х1'25 можно представить как х х. Если не требуется высокая точность или характеристика аппроксимируется для узкой
области |
изменения п, можно принять |
|
т — 1. |
|
|
|
|
|
|
|||||||
|
Для обеспечения лучшего приближения обычно аппроксимиро |
|||||||||||||||
ванные изодромы |
пг |
и п 2 целесообразно |
провести не через крайние |
|||||||||||||
щ |
|
|
|
|
|
|
точки /, 2 исходных |
изодром, |
||||||||
|
|
|
|
|
|
а через какие-то промежу |
||||||||||
1 |
|
|
|
|
|
|
точные точки |
3, |
4 |
(штрихо |
||||||
|
|
|
|
|
|
вая изодрома |
п2 |
на рис. П.5). |
||||||||
|
|
|
|
|
|
|
||||||||||
и |
|
о, |
|
|
|
|
В |
этом |
случае |
Аъ |
Bt, |
Си |
||||
|
|
|
|
|
|
Л 2 , В2 , |
|
С 2 следует |
брать |
по |
||||||
|
|
|
|
|
п=0,8 |
|
|
|||||||||
Ш |
|
|
|
|
|
концам |
аппроксимированных |
|||||||||
|
|
|
|
|
|
|
изодром. |
|
|
|
|
|
|
|||
|
|
|
|
|
/0,9 |
|
|
Границу помпажа |
можно |
|||||||
|
|
|
|
|
\у10 |
|
представить |
зависимостью |
||||||||
|
|
|
0, |
|
|
|
є — 1 = CG + DG2 |
+ EG3, |
||||||||
0,90 |
|
|
|
|
|
|||||||||||
|
0, S |
\ |
|
|
|
|
|
|
|
|
|
(11.37) |
||||
|
|
|
|
|
\\ |
где постоянные |
С, |
D, |
Е |
не |
||||||
|
|
|
|
|
сложно |
подобрать |
|
по |
задан |
|||||||
|
|
|
|
|
ной характеристике. |
Обычно |
||||||||||
|
|
|
|
|
С |
и Z) > О, |
а |
£ |
< |
0. |
Для |
|||||
|
0JS |
1,00 |
105 |
1J0 |
|
области |
режимов |
при оборо |
||||||||
|
|
|
|
|
Л лІНопт |
тах ниже номинальных чаще |
||||||||||
|
|
|
|
|
Щпт |
Н |
всего |
|
можно |
ограничиться |
||||||
Рис. II . 6 . Зависимость |
т) = |
гт,/ї)ш а х от г для |
лишь |
первыми |
двумя |
чле |
||||||||||
различных относительных оборотов n(попт |
= |
нами |
уравнения |
(11.37). |
|
|||||||||||
|
|
= |
0,8) |
|
|
|
Более сложным |
представ |
||||||||
|
|
|
|
|
|
|
ляется |
|
аппроксимировать с |
высокой точностью зависимость к. п. д. компрессора от его параметров. Для этого нужно представить основные зависимости в таких пара метрах, которые позволили бы получить серию достаточно близких друг другу кривых, удобных для аппроксимации. Эту задачу можно
решить, представив для различных |
значений |
п = const отношение |
к. п. д. к его максимальному значению г\ |
в зависимости |
|
от параметра |
|
|
П„пт Г |
" |
|
где Нопт1Н — отношение адиабатного напора в компрессоре на оптимальном и данном режимах; в качестве оптимального прини мается режим, соответствующий точке максимального к. п. д. на характеристике компрессора; п0пт — обороты, соответствующие ма
ксимальному к. п. д. т)т а х (на рис. II.6 и II.7 п о т = 0,8).
62
На рис. II . 6, на котором представлены такие зависимости для компрессора ГТУ-20, видно, что их можно аппроксимировать пара болой, смещая ее вдоль осей координат для различных значений п.
Выражение для т] имеет вид
r) = l2(z-Q)-(z-Q)d]P, (11.38)
где Q, Р и d постоянные.
Константы Q и Р учитывают смещение кривой ц соответственно
вдоль осей абсцисс и ординат. При аппроксимации |
характеристики |
||||||||
в зоне между номинальными и какими-то минимальными оборотами пх |
|||||||||
для любого промежуточного п значения Q и Р легко получить по |
|||||||||
интерполяционным |
формулам: |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
(11.39) |
|
|
|
|
|
|
2 |
|
|
(11.40) |
|
|
|
|
|
|
|
|
|
||
где |
попт |
— относительные обороты, |
при которых |
достигается |
ма |
||||
ксимальный |
К. П. Д. Tlm a x . |
|
|
|
например, |
||||
|
Построив |
по заданной характеристике для оборотов, |
|||||||
попт |
= |
0,8 и |
tii ~ |
0,5 кривые т) (рис. II.6), |
легко |
найти |
Qx и |
Мх |
|
(берутся абсолютные величины), а |
также |
подобрать значение |
d. |
Для определения d достаточно подставить в уравнение (11.38) зна
чение Г[ для любой отдаленной |
от максимума точки на кривой п = |
= п0пт (рис. II . 6, Q = 0, Р = |
1). Обычно d = 2,2-^2,8 и для дан |
ного компрессора может быть принято одинаковым при любом числе оборотов. В случае необходимости можно было бы дополнительно
учесть зависимость d от п.
Таким образом, уравнения (11.33), (11.37) и (11.38) позволяют аппроксимировать характеристику любого компрессора.
Уместно отметить, что в первом приближении, а также при
исследовании достаточно |
узкой |
зоны |
режимов |
ГТУ |
изменением |
|
к. п. д. компрессора ті обычно можно |
пренебречь. |
Аппроксимация |
||||
кривых т] по уравнению |
(11.38) |
при постоянном |
d |
(а иногда и при |
||
Q = 0, Р = 1,0), как правило, обеспечивает достаточную точность |
||||||
расчета. Еще более высокую точность |
можно получить, |
применив |
такой прием. Вначале нужные режимы просчитываются при постоян ном к. п. д. компрессора. Далее аппроксимируется зависимость г\ от оборотов лишь вдоль полученной линии режимов, затем расчет повторяется уже с учетом изменения к. п. д. Однако такой расчет более трудоемок.
Аппроксимация характеристик различных компрессоров под тверждает достаточную точность описанного метода. На рис. II.7 для иллюстрации показана экспериментальная характеристика
компрессора и рассчитанная в результате аппроксимации по уравне ниям (11.33) и (11.38). Совпадение расчетных и экспериментальных кривых хорошее.
Аппроксимация универсальных характеристик турбины к виду, удобному для расчета на ЭВМ
Для точных расчетов и исследования показателей газотурбинных и паротурбинных установок на частичных и переходных режимах обычно используют полученные экспериментально или детальным поступенчатым расчетом характеристики турбины. Такие харак теристики графически изображают зависимость между расходом, числом оборотов, параметрами перед и за турбиной, к. п. д. турбины.
Анализ показывает, что трудно получить простые универсальные зависимости, одинаковые для любой турбины и позволяющие опре делять основные параметры на переменных режимах с высокой точностью. Однако полученные детальным поступенчатым газодина мическим расчетом или экспериментально и заданные в виде графи ков характеристики можно аппроксимировать простыми уравнениями с некоторыми коэффициентами. Форма этих уравнений универсальна, т. е. одинакова для любой турбины. Изменяются лишь коэффициенты, которые легко находятся по заданным характеристикам данной конкретной турбины. Полученные аппроксимацией зависимости можно использовать для анализа и расчета на ЦВМ характеристик данной установки. Ниже предложен способ такой аппроксимации.
Для осуществления аппроксимации необходимо найти такую форму представления характеристик, выбрать такую систему пара метров и координат, которые позволили бы выразить характеристики простыми аналитическими зависимостями. Расходные характери-
стики турбины удобно представить, как на рис. П.8, в виде зависи мостей между относительным приведенным расходом
и"Р~ |
G 0 Pl У Т10 ' |
отношением давлений в турбине р = р2/рх веденными оборотами
Здесь и ниже дополнительным индексом на расчетном режиме, черточкой сверху —
Расходную характеристику при постоянных приведенных оборотах можно с достаточной точностью аппроксимировать дугою эллипса
и относительными при
0 обозначены параметры относительные величины.
|
пр |
і |
(Р — PJU) |
|
|
|
|
|
|
+ |
б2 |
|
|
|
|
откуда, |
учитывая b |
|
|
|
|
||
|
|
|
(р - |
рЛ )2 |
|
|
|
При отношении |
давлений |
в турбине |
|
|
|||
Рм и постоянном ппр |
приведенный |
Рис. II . 8 . |
Схема аппроксимации |
||||
расход |
Gnp |
достигает |
максимума и |
расходных |
характеристик турбины |
||
остается |
неизменным |
для любого |
|
|
р < pV Назовем этот режим критическим для турбины. Он соответствует возникновению критической скорости в одной (любой) из ступеней турбины. Отметим, что иногда при построении расходных характеристик представляется удобным приводить рас ход не к расчетному, а к максимальному (критическому) при рас четных оборотах.
Для измененных оборотов ппр характеристика выразится дугою эллипса с измененными полуосями. Учитывая сравнительно малое изменение полуосей, такая аппроксимация обеспечивает высокую точность.
Приняв
# = |
(1 — kjfinp) |
а0; РЛ = Р л 0 |
k%nnp, |
получим |
|
|
|
Gnp = (1 |
—кіПпр) |
- Р , |
Кппр)' |
at |
|
Здесь коэффициенты k x и k 2 учитывают изменение полуосей эллипса. Их нетрудно определить по заданным характеристикам турбины (рис. II.8):
К = - А |
- |
; К = - |
В - |
> |
(Н.42) |
1гпр" |
ппр |
12Пр |
- |
ппр |
|
И. В . Котляр |
65 |
где п2пр |
и п1пр — максимальные и минимальные |
приведенные обо |
роты на аппроксимируемых характеристиках; А |
и В — диапазон |
|
изменения |
величин а и р л для этих оборотов (рис. |
II.8). |
Необходимо указать, что при аппроксимации следует подобрать
такие значения $м |
для п1пр и п2пр, |
при которых получается наиболь |
шее приближение |
дуги .эллипса |
к действительной характеристике. |
Эти значения не обязательно должны точно совпадать с точками действительных критических режимов, когда в одной из ступеней возникает критическая скорость (часто удобно принять несколько заниженные значения $м ).
Интересно отметить, что в многоступенчатой турбине при неоди наковых на критическом режиме отношениях давлений в ступенях критический расход при различных ппр оказывается неодинаковым. Это объясняется изменением пропускной способности ступеней с докритическим истечением в зависимости от оборотов и соответ ствующим перераспределением перепадов между ступенями. Для одноступенчатой (активного типа) турбины, а также для турбины, в которой критическая скорость возникает одновременно во всех ступенях, А = 0.
Кроме параметров Gnp, р\ ппр на характеристики турбины обычно наносятся еще кривые для определения ее к. п. д. т]. Для аппрокси мации эти кривые удобно представить как зависимость относитель ного к. п. д. ц — г\/г\тах от среднего относительного характеристи ческого числа
|
|
лср |
— |
JL |
і |
/ я о |
|
|
|
|
|
|
V Н ' |
|
|
||||
|
|
|
|
хсро |
Лп |
|
|
||
|
|
ппр |
|
|
|
|
|||
для |
ряда значений |
= |
const |
(рис. |
|
II.9). Здесь |
|
||
|
|
|
|
|
|
|
823 8 0г,2"Я2' ' |
|
|
и — окружная скорость ї-й ступени; Н — изоэнтропийный |
пере |
||||||||
пад |
в турбине. |
|
|
|
|
|
|
|
|
|
Кривые к. п. д. |
хорошо аппроксимируются |
уравнением |
пара |
|||||
болы |
|
|
|
|
|
|
|
|
|
|
|
Ч = |
[2Сх — (Cx)*]D; |
(11.43) |
|||||
|
С = \ + |
ka{\—ппр); |
D = |
1 — M l — |
пПр). |
|
Здесь величины G и D учитывают смещение кривых к. п. д. вдоль осей абсцисс и ординат. Коэффициенты k3 и kit постоянные для данной турбины, легко находятся по характеристикам (рис. II.9)
|
k3 |
=; - |
Е |
- |
; |
К = - |
1 |
+ |
• |
(П.44) |
|
|
п2ПР |
— |
П1ПР |
|
n2ltp |
— |
Пщр |
|
|
Возможность |
линейной |
интерполяции |
в интервале |
оборотов |
||||||
п2пр—п1пр |
определяется |
|
здесь |
малыми |
смещениями Е и |
Fx + F%. |
66
• — по |
р е з у л ь т а т а м д е т а л ь н о г о п о с т у п е н ч а т о г о |
расчета; |
|
— по |
а п п р о к с и м а ц и о н |
ным у р а в н е н и я м (11.41), |
(11.43), |
|
Р Р — точка |
расчетного р е ж и м а |
|
Отметим, что в.качестве п1пр можно принять любое малое ппр, |
кроме |
|||||
ппр |
= 0, |
когда |
хср |
обращается |
в неопределенность. |
|
При |
умеренном |
рабочем диапазоне изменения ппр обычно |
можно |
|||
принять |
44 = |
0 и О = 1 , а в |
некоторых случаях (для приближен |
|||
ных |
расчетов) |
и |
С = 1. |
|
|
Таким образом, по уравнениям (П.41) и (11.43) можно довольно просто аппроксимировать для аналитических исследований и рас четов на ЭВМ характеристики любой турбины, которые должны быть заданы.
На рис. 11.10 сравниваются характеристики двухступенчатой турбины активного типа, полученные в результате детального поступенчатого расчета с характеристиками, рассчитанными по урав
нениям |
(11.41) и |
(11.43). |
Совпадение указанных |
характеристик |
хорошее, |
что подтверждает |
высокую точность аппроксимации. |
||
В заключение |
укажем, |
что расходные характеристики турбины |
||
вместо уравнения |
(11.41) можно аппроксимировать |
еще более про |
стым и удобным для расчетов на ЭВМ уравнением (с теми же значе ниями коэффициентов)
\3п
(11.45)
\1 — Рло + М я р / J
Детальный анализ показывает, что это уравнение обеспечивает практически такую же точность аппроксимации, как и уравнение (11.41), кроме области, где |3 близко к единице, которая обычно прак тического интереса не представляет.
Г Л А В А I I I
ИССЛЕДОВАНИЕ ПЕРЕХОДНЫХ ПРОЦЕССОВ С ПОМОЩЬЮ ЛИНЕАРИЗОВАННЫХ УРАВНЕНИЙ.
АНАЛИЗ СОБСТВЕННОЙ УСТОЙЧИВОСТИ ГТУ
11. ЛИНЕАРИЗОВАННЫЕ УРАВНЕНИЯ ДВИЖЕНИЯ ГТУ
Изложенный в п. 7 метод расчета переходных процессов сводится к численному решению системы нелинейных дифференциальных уравнений, характеризующих работу всех элементов ГТУ. Этот метод сравнительно трудоемкий, но обеспечивает достаточно высо кую точность расчета.
Линеаризация уравнений движения ГТУ позволяет получить систему, которая довольно просто решается. Линеаризованные уравнения в дифференциальной форме или в конечных разностях удобно использовать для исследования переходных процессов в ГТУ и решения некоторых задач динамики, например расчета собственной устойчивости установки.
Будем |
предполагать, что характеристики |
компрессоров |
и тур |
||||||||||||
бин заданы в |
виде |
зависимостей: |
|
|
|
|
|
|
|
||||||
|
|
|
Gnp. к — Gfip. к (Є , Л-Пр. к)і |
I |
|
^ |
|||||||||
|
|
|
|
Лк = |
-Лк^,П'пр.к)\ |
|
I |
|
|
|
|||||
|
|
|
Gnp. т= |
Gnp. T (e> nnp. |
T)\ |
I |
|
|
2^ |
||||||
где n n p .K = пк/Утг ^ |
|
= |
|
|
'W. r)> |
|
J |
|
|
|
|||||
|
пк |
(в течение |
переходного |
процесса можно |
|||||||||||
полагать |
7 \ |
const), |
|
|
|
|
|
|
|
|
|
|
|
|
|
Здесь |
и Т 3 — температуры |
перед компрессором.и турбиной. |
|||||||||||||
Для практических расчетов удобно дополнительно построить |
|||||||||||||||
характеристики |
приведенных |
крутящих |
|
моментов: |
|
|
|||||||||
|
Мпр. к |
= -=— = = |
р |
|
.0 |
2 8 |
ч - |
= |
М.пРшк (є, nn ; , к ) ; |
(III.3) |
|||||
|
|
Pi |
«ЯР. к(ео' |
— 1 ) % |
|
|
|
|
|
||||||
|
^пр. г= - |
Г = |
=~^~ Т. |
—О 25 \ Т |
~ М-пр. т (е> ппр. т)> |
(HI -4) |
|||||||||
|
|
04 |
|
Ппр. т \ 1 —е0 ' |
j |
|
|
|
|
|
|||||
где р х и р 4 — давление |
перед компрессором и за турбиной. |
|
|||||||||||||
В линеаризованном |
виде |
р а с х о д н ы е |
и м о м е н т и ы е |
||||||||||||
х а р а к т е р и с т и к и |
|
|
к о м п р е с с о р а |
и |
т у р б и н ы |
||||||||||
(III.1)—(III.4) запишутся: |
|
|
|
|
|
|
|
|
|
|
|
||||
|
AGnp,K |
= eAs + PAnnp,K, |
{Аппр.к^Апк)\ |
|
(III.5) |
||||||||||
|
|
ДОл р .T |
= dA~e + S {Апт — 0,5 ДТ3 ); |
|
(Ш.6) |
||||||||||
|
|
|
АМ„Р .к |
= 9К Де + Рк |
Апк; |
|
|
(ІІІ.7) |
|||||||
|
|
АМПр. |
T = dTAe |
+ ST |
{Апт — 0,5 АТа). |
|
(Ш.8) |
||||||||
Здесь |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dGnp. |
к \ |
|
. |
|
p_f_^Gnp1K |
|
|
|||||
|
|
|
|
дъ |
|
Jnnp. |
к |
|
|
|
\дпПр.к/е |
|
|
||
|
|
|
\ |
де |
|
/ п |
п р |
. к |
|
|
|
\дпПр.к/е |
|
|
|
|
|
|
dGjip. г |
\ |
|
. |
S = |
t дйпр. |
т |
|
|
||||
|
|
|
|
де |
|
1ппр.т |
|
|
|
\дппр.т/е |
|
|
|||
|
|
|
\ |
де |
|
/"пр. т |
|
|
\ |
дппр. |
т /е |
|
|
||
Эти частные производные |
определяются |
по графикам |
зависимостей |
||||||||||||
(ПІЛ)—(ІІІ.4). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Линеаризованное |
у р а в н е н и е |
р о т о р а |
|
имеет |
вид |
||||||
*> Чіт |
+ ^ 7 |
А ^ |
+ |
^ |
^ |
+ |
- |
д |
^ |
= |
о. < ш - 9 ) |
где Rp = |
—— |
время |
ротора; |
Мк, |
Мд, МТр, Мт |
— |
отнесенный |
к значению на начальном установившемся режиме момент компрес сора, внешнего потребителя (если он размещен на данном валу), трения, турбины.
С индексом 0 здесь и ниже — параметры на начальном установив
шемся |
режиме. |
|
|
|
|
|
|
Линеаризованное |
у р а в н е н и е |
г а з о в о з д у ш н о й |
|||||
е м к о с т и |
между компрессором и турбиной в ГТУ с регенерацией |
||||||
имеет |
вид |
|
|
|
|
|
|
|
|
|
R |
v ^ - |
+ |
AGK-AGT==0, |
|
где Rv |
Vv |
— время |
емкости; |
V — ее |
объем. |
||
= ^ |
|||||||
|
"го |
|
|
|
|
|
|
В соответствии с уравнением состояния, учитывая, что темпе |
|||||||
ратура |
воздуха |
за регенератором |
изменяется мало, запишем |
||||
|
|
|
Ду = Ар — AT |
Ар |
А~ёобщ, |
||
где р—давление |
в емкости; |
гобщ— |
общая |
степень сжатия компрес |
соров; учитывая малое изменение температуры воздуха за регене
ратором, принято АГ |
0. Подставив Ау в предыдущее |
уравнение, |
|
получим |
|
|
|
Rvd { А |
^ щ ) -\- AGK — AGT = 0. |
|
(ШЛО) |
С в я з ь м е ж д у т е м п е р а т у р а м и |
п е р е д |
д в у м я |
|
п о с л е д о в а т е л ь н о |
в к л ю ч е н н ы м и |
т у р б и н а м и |
(или группами турбин) при отсутствии промежуточных камер
сгорания |
имеет вид |
[40] |
|
|
|
|
|
||
|
|
|
A n |
= |
ДТз - |
0,2 ДІ! А Та, |
(III. 11) |
||
где |
Тз — температура |
перед |
второй |
турбиной; ех — |
отношение |
||||
давлений |
в первой |
турбине. |
|
|
|
|
|||
У р а в н е н и е |
б а л а н с а |
о б щ е й с т е п е н и |
с ж а |
||||||
т и я |
и |
р а с ш и р е н и я |
запишется |
|
|||||
|
|
|
|
|
|
п |
|
ш |
|
|
|
|
^общ |
= 2 Д і , = |
2 Аё,. |
(III. 12) |
Здесь т и п — число турбин и компрессоров, включенных после довательно.
Если расходная характеристика турбины (III.1) и моментные характеристики турбины и компрессора (II 1.3), (II 1.4) отсутствуют, можно воспользоваться для расчетов следующими линеаризован ными уравнениями [40];