Газодинамический расчёт проточной части турбинной ступени активного типа

 

В результате предварительного теплового расчёта паротурбинной установки мы находим следующие параметры рабочего тела перед рассчитываемой ступенью:

·        расход рабочего тела G0 ,

·        начальные температура T0 и давление p0 рабочего тела,

·        давление потока на выходе из ступени p2, или теплоперепад ступени H0,

В процессе газодинамического расчёта турбинной ступени в ряде случаев приходится прибегать к выбору таких параметров, как:

·        форма профилей сопловых и рабочих решёток;

·        высоты решёток ;

·        углы установок лопаток в решётках;

·        конструкция бандажей рабочих лопаток и лабиринтных уплотнений.

·        направление потока при входе в ступень α0,

·        выбор степени реактивности ступени ρ,

·        выбор отношение скоростей u/сф.

 

Выбор степени реактивности

 

Ступени могут быть либо активные (ρ = 0.02 – 0.25), либо реактивные (ρ > 0.4).

Активные ступени выбираются в том случае, когда  d/l > 8, или когда подвод рабочего тела в ступени выполняют парциальным (e < 1.0)

НЕ нашли? Не то? Что вы ищете?

В ступени активного типа основная часть теплоперепада  перерабатывается в сопловой решётке. Профили сопловых и рабочих лопаток такой ступени существенно отличаются друг от друга. В рабочей решётке происходит незначительное ускорение потока при большом угле его поворота.

В ступени реактивного типа характер обтекания сопловых и рабочих решёток практически одинаков и сами профили по своей форме геометрически подобны.

С увеличением ρ улучшается обтекание рабочей решётки, где поток становиться более конфузорным. В связи с этим относительный лопаточный КПД возрастает.

Особенностью ступени с повышенной реактивностью является увеличение усилия, действующего на лопатки и диск в осевом направлении.

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

Таким образом, в конечном счёте выбор степени реактивности ρ является технико-экономической задачей.

 

 Выбор отношения скоростей u/cф

 

В зависимости от степени реактивности определяется оптимальное отношения скоростей u/cф , обеспечивающее максимальную величину относительного лопаточного КПД ступени.

Следует также учитывать, что уменьшение отношения скоростей u/cф<(u/cф)опт , позволяющее при той же окружной скорости переработать больший теплоперепад в ступени, с одной стороны, снижает КПД, а с другой – уменьшает число ступеней или диаметр ступени и тем самым удешевляет изготовление турбины.

В общем случае:

 . (1)

 

Определение основных размеров ступени

 

.

Рис 1

 

При заданном теплоперепаде ступени и выбранном значении отношения скоростей u/cф  диаметр ступени равен:

 

, [м] . (2)

 

Иногда по заданному значению диаметра можно найти располагаемый теплоперепад ступени:

 

, [Дж/кг]. (3)

 

Выходная площадь сопловой решётки для дозвукового режима (M1t < 1) или можно найти из уравнения неразрывности:

 

, [м2], (4)

 

где скорость , а удельный объём  v1t определяется по h, s- диаграмме в конце изоэнтропного расширения рабочего тела в сопловой решётке (рис. 1).

Коэффициенты расхода для сопловой решётки в первом приближении можно определить по формуле:

 

, (5)

 

или принять в первом приближении μ1 = 0.97.

 

При сверхзвуковых скоростях потока (M>1) или ε1 < ε* обычно также применяются суживающиеся решётки, но выходная площадь находится из уравнения:

 

; [м2], (6)

здесь v1* и с* соответствуют критическому отношению давлений ε* (рис. 2) или критическому теплоперепаду , где:

 

. [м/с], (7)

 

Выходная высота сопловой решётки l1 находится из выражения :

 

; [м], (8)

 

здесь e – степень парциальности – длина дуги, занятой сопловой решёткой, отнесённой ко всей окружности:

. (9)

 

Рис. 2

 

Эффективным углом выхода α1э следует задаться, учитывая, что, с одной стороны, желательно уменьшить  α1э для того, чтобы увеличить высоту лопаток и повысить КПД ступени; а, с другой стороны, уменьшение (α1э < 11o) ведёт к росту профильных потерь в решётках.

По величине α1э , заданному углу входа  α0 и числу М1t выбирается профиль сопловой решётки, а по аэродинамическим характеристикам выбранной решётки определяются угол её установки αу и относительный шаг t.

Хорда профиля решётки b1 (рис.1) выбирают с таким расчётом, чтобы обеспечить достаточную прочность лопатки и жёсткость диафрагмы.

Обычно для активных ступеней b1 = 40 – 80 мм ; для реактивных ступеней - составляет b1 = b2 = 20 – 60 мм.

После выбора b1 должна быть подсчитана относительная высота сопловой решётки  и уточнены относительный шаг решётки  и число лопаток  zc .

Для вычисления действительной скорости  с1 необходимо определить коэффициент потери энергии  ζс по приближённой формуле (10).

 

, (10)

Для уточнения значений ζс и μ1 необходимо учитывать поправки на число Re1 = c1tb1/v1 .

Следующим этапом расчёта ступени является построение входного треугольника скоростей, определение относительной скорости входа рабочего тела в рабочую решётку w1 и угла её направления β1 (рис. 3).

 

Рис. 3

 

Скорость потока может быть определена по формуле с1 = φc1t . Действительный угол выхода потока из сопловой решётки можно определить по формуле (11)

. (11)

 

Для расчёта рабочей решётки необходимо знать состояние рабочего тела перед ней, для чего следует подсчитать потери энергии в сопловой решётке:

 

, [Дж/кг]. (12)

 

Высоту рабочей лопатки можно определить по формуле:

 

, [м], (13)

 

где Δlп = 1.5 – 2.5 мм и Δlк = 0 – 1.5 мм соответственно являются перекрышами по периферии и корневому сечений лопатки.

В последних ступенях конденсационных турбин допускается увеличение перекрыши до 20 мм.

Выходную площадь рабочей решётки для докритического режима , т.е. при или  , где  - давление торможения в относительном движении (рис. 2), находим из уравнения неразрывности:

 

, [м2]. (14)

 

Скорость w2t находим по формуле:

 

, [м/с] , (15)

 

а удельный объём v2t находим по h, s- диаграмме в конце изоэнтропного расширения в решётке ( рис. 2).

Коэффициент расхода μ2 можно определить по приближённой формуле (16). В первом приближении μ2 = 0.93.

 

. (16)

 

При сверхзвуковой скорости потока выходную площадь рабочей решётки находим по формуле:

,[м2] ; (17)

здесь  v2* и w* соответствуют критическому отношению давлений или критическому теплоперепаду  H* = 0.5w2* , где:

 

. (18)

 

В большинстве ступеней l2 = l`2 ; в последних ступенях конденсационных турбин принимают обычно l2 > l`2 .

При заданном значении l2 можно определить эффективный угол выхода для рабочей решётки:

. (19)

 

По величине β2э , примерному значению угла входа  β1 , которое может немного отличаться от β1опт, и числу М2t выбирается профиль рабочей решётки, а по аэродинамическим характеристикам выбранной решётки определяются угол её установки βу и относительный шаг t2.

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

 

. (20)

 

Усилие R, действующее на лопатки, можно подсчитать по формуле:

 

, (21)

 

где Ru и Ra определяем по (22) и (23); Wмин – момент сопротивления профиля лопатки.

Окружное усилие может быть найдено на основании уравнения количества движения, записанного для оси u при массовом расходе рабочего тела, равном G, кг/с:

. (22)

 

Осевое усилие Rа может быть найдено из уравнения количества движения в направлении оси а, учитывая при этом разность давлений р1 - р2 , действующих на кольцевую площадь рабочих лопаток Ω=πd2l2:

 

(23)

 

В выражениях (22) и (23) α1 и  - углы направления скоростей с1 и с2 (рис. 4).

Рис. 4

 

Для нержавеющих сталей обычно применяют .

Если выбранный размер профиля не удовлетворяет требованиям прочности, то при сохранении подобия всех размеров решётки профилей хорду следует увеличить в соответствии с выражением:

 

. (24)

 

Определение КПД ступени

Потери энергии в рабочей решётке определяют по формуле:

 

. (25)

 

Коэффициент потерь энергии ζр и коэффициент скорости ψ можно определить по формулам (26) и (27).

 

. (26)

. (27)

 

Для уточнения значений μ2 и ψ необходимо учитывать поправки на число Re2 , где Re2 = w2t b2 / v2 .

Построение выходного треугольника (рис. 3) производится по скорости w2= ψw2t и углу  β2 . Для дозвуковых скоростей и перегретого пара приближённо можно принять β2 = β2э . Для сверхзвуковых скоростей (М2t >1) с учётом отклонения потока в косом срезе рабочей решётки угол

β2 = β2э + δ определяется по формуле (28).

 

. (28)

 

Из треугольника скоростей определяют выходную скорость с2 и потери с выходной скоростью ΔHв.с. .

Далее строим процесс в h, s- диаграмме (рис. 5).

Рис. 5

По формуле (29) определяем окружную работу Hu и по формуле (30) - мощность на лопатках Nu .

 

, [Дж/кг]. (29)

 

Nu = GHu, [Вт]. (30)

 

По формулам (31) и (32) определяем относительный лопаточный КПД ступени ηо.л. ;

 

. (31)

 

. (32)

 

Пример расчёта активной ступени паровой турбины

 

Рассчитать ступень активного типа паровой турбины при следующих условиях:

·        расход пара G = 60 кг/с;

·        давление пара перед ступенью p0 = 5,00 МПа;

·        температура пара перед ступенью t0 = 5000 C ;

·        средний диаметр ступени d1 = d2 = 1.1 м;

·        частота вращения n = 50 c-1 ;

·        скорость пара на входе в ступень с0 = 50 м/с;

 

Чтобы получить максимальный теплоперепад,

вырабатываемый на одновенечной турбинной ступени с максимальным КПД, зададимся степенью реактивности ρ = 0.1; коэффициентом скорости сопловай решётки φ = 0.96; углом входа потока в сопловую решётку α0= 900; углом выхода потока из сопловой решётки α1= 120 и коэффициентом расхода μ1 = 0.97.

Тогда оптимальное отношение скоростей будет равно:

 

= 0.495.

 

Определяем окружную скорость:

u = u1 = u2 = πdn = 172.8 м/с.

 

Находим фиктивную скорость: = 349.1 м/с.

Определяем располагаемый теплоперепад ступени:

 = 60932 Дж/кг.

Зная начальные параметры пара p0 = 5,00 МПа и  t0 = 5000 C, по h,s-диаграмме определяем энтальпию h0 = 3434.655 кДж/кг и υ0 = 0.0686 м3/кг.

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

 

= 3434655 + 1250 = 3435905 Дж/кг.

 

Теоретическая энтальпия пара на выходе из ступени равна:

 

= 3374973 Дж/кг

 

Располагаемый теплоперепад в сопловой решётке равен:

 

= 54839 Дж/кг

 

Определяем теоретическую скорость пара на выходе из сопловой решётки:

 

= 331.18 м/с.

 

Теоретическая энтальпия пара на выходе из сопловой решётки равна:

 

 = 3381066 Дж/кг.

 

По h,s-диаграмме по энтальпии h1t = 3381066 Дж/кг нахадим p1 = 4.26501 МПа и υ1t = 0.07767 м3/кг.

Скорость звука на выходе из сопловой решётки равна:

 

=656.2 м/с.

 

Т.к. M1t = c1t /a1 = 0.50 < 1, то площадь сопловой решётки определяется по формуле:

=0.0145 м2.

Приняв полный подвод пара (e=1), найдём высоту сопловой лопатки:

 

=0.0202 м.

Согласно полученным данным, выбираем из атласа профилей турбинных решёток профиль С-90-12А

По конструктивным соображениям выбираем хорду профиля b1 =0.05 м, толщину выходной кромки Δкр =0.8 мм и относительный шаг =0.8.

Тогда число сопловых лопаток будет равно:

 

= 86.

 

Определяем число = 7.65х106, где ν1 = 2.1375 х 10-6 м2/с – кинематическая вязкость пара в точке 1 h,s-диаграммы.

Уточняем коэффициент расхода по формуле (5):

= 0.9696,

что приблизительно равно ранее заданному значению.

Коэффициент скорости находим по формуле (10):

 

 = 0.960

Находим действительную скорость пара на выходе из сопловой решётки:

c1 = φ c1t = 318.0 м/с

Зная c1 , u и угол α1 ,строим треугольник скоростей и находим w1 и угол β1 .

Однако w1 и угол β1 можно найти по формулам:

 

 = 153.25 м/с;

 

 = 25.560 .

 

Находим теоретическую относительную скорость по формуле (15):

 

 = 188.83 м/с.

 

По формуле (12) находим потери в сопловой решётке:

 

 = 4299 Дж/кг.

 

Определяем энтальпию пара на выходе из сопловой решётки:

 

 = 3385365 Дж/кг.

По h,s-диаграмме, зная h1 и р1, находим энтрапию s1 = 6.9838369 кДж/кг .К .

Далее по h,s-диаграмме из точки 1 производим изоэнтропное расширение пара в рабочей решётке ступени до энтальпии пара h2t = 3379272 Дж/кг и определяем параметры пара в точке 2t :

р2 = 418735 Па; ν = 2.17926 м2/с ; υ2t = 0.0796 м3/кг.

 

Определяем  - давление торможения в относительном движении (рис. 2).

Для этого находим величину w12/2 = 11743 Дж/кг, тогда =3397108 Дж/кг и  = 4417870 Па.

Определяем р1* = ε* = 241083 Па , что на много меньше p2 = 4155250 Па.

Определяем М2t :

 = 0.28.

Задаём μ2 =0.95 и определяем площадь рабочей решётки по формуле (14):

 

 = 0.02663 м2.

 

По формуле (19) находим эффективный угол выхода потока из рабочей решётки:

 

 =0.3588; β2э =21.030 , где l2 = l1+∆l =0.022 м.

Т.к. поток дозвуковой , то угол выхода потока из рабочей решётки можно определить по формуле:

= 0.35018; β2 =20.50 , где задаёмся : ψ = 0.94 .

Находим действительную относительную скорость потока на входе в рабочую решётку:

= 176.9 м/с .

 

Зная w2 , u и угол β2 ,строим треугольник скоростей и находим c2 и угол α2 .

Однако c2 и угол α2 можно найти по формулам:

 

 = 62.35 м/с;

 

 = 96.540 .

Согласно полученным данным, выбираем из атласа профилей турбинных решёток профиль Р-30-21А, который имеет следующие параметры:

 = 0.63 ; Wmin =2.34х10-5 м3 ; b2 =0.0256 м ;

 

Уточняем по формуле (16) :  = 0.953 и

 

По формуле (27) =0.944

По формуле (22) определяем окружное усилие, действующее на рабочую решётку:

 = 18237 Н;

 

по формуле (23) находим осевое усилие, действующее на рабочую решётку:

 

 = - 5900 Н.

 

Определяем суммарное усилие по формуле (21):

 

 =19167.6 Н.

 

По формуле (20) определяем изгибающие напряжения на рабочих лопатках:

 

 = 42105 Па,

где =214 шт.; что ниже допустимого  МПа

Уточняем коэффициент скорости для рабочей решётки:

 

. =0.944

По формуле (25) определяем потери энергии в рабочей решётке:

 

 = 1941 Дж / кг.

Определяем потери энергии с выходной скорорстью:

 

 =1944 Дж / кг.

Определяем относительный лопаточный КПД ступени:

 

= 0.865

 

Проверяем относительный лопаточный КПД ступени по формуле (32):

 

 = 0.866