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

книги из ГПНТБ / Митрохин В.Т. Выбор параметров и расчет центростремительной турбины на стационарных и переходных режимах

.pdf
Скачиваний:
96
Добавлен:
24.10.2023
Размер:
9.12 Mб
Скачать

При v = 0 из (3.45)

получим

 

 

 

 

 

8 и

V

"1

0

"

 

1 - М 2

 

 

 

 

 

 

 

 

X

_^21

^22.

 

0

sin %

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

h,(l)l,

 

 

— е

 

 

 

 

 

 

 

 

 

 

 

 

 

 

An

An

1- М2

 

 

 

-4ц

А12

 

X

е

*»>,

 

 

(3.62)

A<i\

А<&

 

 

А2\

А22

 

 

 

 

 

 

 

 

т

т .

 

 

 

 

 

 

 

h,(l)l,

 

 

 

 

 

 

 

 

J V U . V l - e 1 - " 2 -

1

 

 

 

Если

источник

возмущения

расположен

перед

решеткой и

£о=1, то из (3.47)

и (3.48)

получим

 

 

 

 

 

 

 

 

 

Р0а 0

811 "I- 012

 

 

(3.63)

 

 

 

 

 

 

 

21 +

 

 

 

 

 

 

 

Pi «1

 

 

 

 

 

 

 

 

/=•„ =

811 + 812.

 

 

(3.64)

Если

источник

возмущения

расположен за

решеткой,

то при

£о= — 1 получим

соответственно

 

 

 

 

 

 

 

 

 

 

Р0«0

812 — »1

 

 

(3.65)

 

 

 

 

 

р 2

а 2

022 —

012

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рр = ЬЦ — 8 12-

 

 

 

(3.66)

Рассмотрим влияние коэффициента радиальности ц на коэф­ фициенты усиления вращающихся решеток центростремительной турбины.

Для расчета были взяты три решетки с коэффициентами ра­ диальности (х = 0,3; 0,5 и 0,8. На входе в решетки в абсолютном движении было принято: М С о = 1 , ао = 20°, р0 = 90°, Гх = 0,1 м и Q = = 2000 1/с. В выходном сечении решеток был принят осевой вы­ ход стационарного потока в абсолютном движении и в соответ­ ствии с изменением коэффициента радиальности р, от 0,3 до 0,8 изменялось число MW2 от 0,299 до 0,845. Проточные части и тре­ угольники скоростей рассмотренных вариантов приведены на

рис. 3.4. При изменении коэффициента радиальности

р, ИЗМеНЯ-

О.2

dr

Поэтому

лись число МЮ г и величина параметра h (/) =

г — .

в приведенных ниже результатах фактически рассмотрено сов­ местное влияние параметра h(l) и М№ 2 на коэффициенты усиле­ ния вращающихся решеток.

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

130

источник возмущения расположен

перед решеткой

(£о=1).

При

ц = 0,8 (случай, близкий к решетке

рабочего колеса

осевой

тур­

бины) выходное сопротивление решетки близко к так называе­

мому «закрытому концу».

По мере уменьшения коэффициента

радиальности и уменьшения

Ми ,2 выходное сопротивление решет­

ки уменьшается. Во всем

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

импеданса изменяет знак. Это значит, что, например, увеличение

пульсаций давления в данном сечении сопровождается

уменьше­

нием пульсаций

скорости.

 

 

 

Зависимость

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

усиления по давлению

F

=—

 

а0=20°, fi0

= 90°

 

Pi

 

 

 

/1=0,3

р.=0,5

А=0,8

Рис. 3.4. Профили проточных частей и тре­ угольники скоростей

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

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

FpZ Fр с а • Fр р ,,.

Если источник возмущения расположен перед турбиной, ве­ личина F p s существенно меньше единицы.

5*

1$1

Рассмотрим теперь коэффициенты усиления вращающихся решеток рабочего колеса центростремительной турбины в случае, когда источник возмущения расположен за турбиной. Зависи­ мость £2 от (х показана на рис. 3.7. Качественно зависимость импеданса от коэффициента радиальности не изменяется по срав-

0,3 Ofi

0,7

а

0,3 0,4-

0,7

А

Рис. 3.5.

Зависимость

импедан­

Рис. 3.6. Зависимость Fv

от \i

при

са

?2 от ц при £ 1= 1

 

 

 

 

0,3 OA

0,5

0J

А

 

0,3

OA

0,7

/ I

Рис. 3.7. Зависимость импеданса £г

Рис.

3.8. Зависимость Fp

от р.

от (х при £1= —1

 

при £1= —1

 

 

нению с рассмотренным выше случаем:

при увеличении

| i импе­

данс решетки

все более приближается к так называемому «закры­

тому концу»

и не меняет своего знака.

Однако

количественно

значения £2

по сравнению со случаем,

когда

£о=1, отличают­

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

Значения коэффициентов усиления

по давлению Fv при

£о=

= — 1 в зависимости от \i приведены

на рис. 3.8. Величина

Fv в

132

рассматриваемом случае существенно зависит от ц. Если при ц.= = 0,3 колебания давления за решеткой почти в два раза больше, чем колебания давления на входе, то при jx = 0,8 колебания дав­ ления за решеткой составляют ~0,8 от колебаний на входе. Со­

поставление данных,

приведенных

на рис. 3.6 и 3.8, обнаружи­

вает существенное различие в зависимостях Fp от jx при £о = 1 и

ci = — 1 . Это различие

объясняется

тем, что когда источник воз­

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

претерпевал

существенные изменения.

 

 

 

 

 

Заканчивая рассмотрение коэффициентов усиления решеток,

определим их для следующих трех

характерных случаев.

 

1. В решетке

подводится

механическая

энергия так, что

Ас

 

 

 

 

 

 

 

 

 

= 90°;

= 20°; р2

—!•'= — 1 . Это условие выполнено

так:

 

 

 

 

 

= 90° и а2 =

20°. При этом

числа

М

задавались

следующим

образом: ЬАС

= 0,2737; М Ю 1 = 0 , 8 ;

М № г =

0,2595.

 

 

2. Неподвижная

решетка

 

=

0

при

М 1 =

0,8,

М 2 =

= 0,2737; а0 = ро = 2О°; р2 =

а2 =

90°.

 

 

 

 

 

 

3. От потока отводилась механическая

энергия

^ ^ - = 1 ^ .

При этом: а0 ==20°;

р 0 = 9 0 ° ;

р 2 = 2 0 °

и а 2 = 9 0 ° .

Числа

М за­

давались так: Mf„ — 0,8; М Ш 1 = 0,2737;

M ^ 2

= 0,85.

 

 

Вращающиеся

решетки

рассматривались

при ц = 1 , поэтому

значения коэффициентов суммарной матрицы 6,j рассчитывались по следующей формуле:

J12

l A i J 22.

1

о

1 0

Ап

Ап

1 0

An

An

(3.67)

|_0

sin р2 .

0 1

A2i

Ачч.

0 1

All

An

 

Поскольку коэффициенты усиления решеток рассматривались при £i = l и £i= — 1 , расчет £г и Fp проводился по формулам (3.64) —(3.66).

Зависимости выходного импеданса решетки £2 и коэффициен­ та усиления по давлению Fp от Ас„/и в случае, когда источник возмущения расположен перед решеткой (£о=1), приведены на

рис. 3.9 и 3.10. Выходное сопротивление решетки при — = 1

и

близко к характеристическому. В неподвижной решетке сопро­

тивление существенно больше, причем с уменьшением Аси

вы­

ходное сопротивление имеет

максимум. Отметим, что как импе­

данс £2, так и величина Fv

в рассматриваемом

случае меняют

знак.

 

 

 

 

Когда источник

возмущения расположен за

решеткой

(£о =

= — 1), выходной

импеданс

£2 и коэффициент усиления по дав­

лению немонотонно зависят

от Acju (рис. 3.11

и 3.12). Обра-

133

тим внимание, что как £2 , так и Fp меняют при этом свои знаки. При рассмотрении систем, у которых отсутствует среднее движе­ ние, поток энергии определяется осредненным за период значе­ нием p'w' (см. ниже) и изменение знака импеданса свидетель­ ствует о том, что рассматриваемый элемент из пассивного прев-

 

АСи_

 

и

Рис. 3.10. Зависимость Fp от

Аси

при £ ;=1

 

 

 

и

и

 

 

 

Рис. 3.11. Зависимость £2

от АсиРис. 3.12. Зависимость Fp от

Аси/и-

при

£ 1= —1

при £ i= —1

 

ращается в

элемент,

генерирующий волновую энергию.

Для

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

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

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

134

3.8. ОПРЕДЕЛЕНИЕ УСЛОВИЙ ГЕНЕРАЦИИ ВОЛНОВОЙ ЭНЕРГИИ

ВРАЩАЮЩЕЙСЯ РЕШЕТКОЙ

Определенные выше динамические коэффициенты Ьц позво­ ляют вычислить основные'динамические характеристики реше­ ток. Помимо этих характеристик, определяющих динамические реакции, с помощью коэффициентов 6,j можно определить еще •одну важную характеристику. Решетки, как и всякий линейный четырехполюсник, могут подразделяться на два важных класса: пассивные и активные. Начнем с некоторых определений. Пусть на входе или выходе из решетки есть источник возмущения. Рас­ смотрим установившийся режим работы при синусоидальном воз­ буждении. Для такого синусоидального установившегося режима колебаний можно определить средний за период поток колеба­ тельной энергии я. Если поток колебательной энергии на входе в четырехполюсник щ, а на выходе из него я 2 , то пассивным на­ зывается четырехполюсник,, у которого соблюдается условие

л^ — я 2

= я > О,

(3.68)

где я — энергия, потребляемая

элементом.

 

Условию (3.68) эквивалентно соотношение, называемое коэф­

фициентом усиления мощности,

 

F„ =

i ^ . < l .

(3.69)

л1

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

Наша задача — выразить величину коэффициента усиления мощности через динамические коэффициенты 6*j. Для этого прежде всего необходимо определить средний за период поток энергии через сечение.

Выведем соотношение, определяющее баланс энергий рас­ сматриваемой линеаризованной системы, описываемой уравне­ ниями (3.5) и (3.6). Для составления баланса энергий необхо­

димо

образовать

по

переменным р'

и

w'

с учетом

уравнений

(3.5)

и (3.6) квадратичную

форму и

при

этом отделить члены,

содержащие p'w',

р12

и w12

под знаком

производных

d/dt и д/д1.

Такой прием широко употребляется в математической физике *. Для акустических систем, в которых среднее движение не­ обходимо принимать во внимание, такой подход впервые осуще­ ствил В. Л. Эпштейн. Прежде всего умножим скалярно уравне­

ние (3.5) на

р, а уравнение

(3.6) — на

1/р. С учетом

(3.7)

будем

иметь

 

 

 

 

 

 

 

 

d(w'~p)

, - д . -

,. .

-

д г » 2

Жг^л

др'

 

 

 

p ¥ ( w w

) + р

р

ё Г [ ~ — ]+ — =

0 '

( 3 J 0 )

* Для желающих подробно ознакомиться с этим вопросом можно рекомен­ довать книгу С. К. Годунова «Уравнения математической физики». М , «Нау­ ка», 1971, 416 с.

135

где учтено, что р не зависит от i.

 

 

 

 

 

 

 

 

 

Для уравнения

(3.6)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W

 

а? )

dw'

 

 

 

 

dw

 

In A .

dt \"р

а 2

 

 

dl

dl

 

 

 

 

 

a/

 

• да

- ) -f-

 

 

 

 

 

 

рД2

 

 

 

 

 

р

 

 

 

 

 

 

 

dl

1

 

 

 

 

 

 

 

 

, (

d In p

,

d l n

/

 

=

0.

 

 

 

(3.71)

 

 

 

 

 

 

 

 

 

a/

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для образования квадратичной формы надо умножить урав­

нение

(3.70) скалярно

на до'. Рассмотрим

 

произведение

,

d

,

ч

 

 

 

 

 

d ,

 

 

 

dw'

 

да'

(да, p).

да—(да

дар),

 

поскольку(дадар)= дар

 

 

 

 

a/ v

 

dl

v

 

 

 

 

 

dl v

 

 

 

 

^

 

 

 

 

Умножим правую

и левую части этого соотношения на ДО

 

 

 

чю'

 

 

,—-

dw'

1

iil2

 

дар

 

 

 

 

 

(да'дар) = да'дар —

 

- ^

 

 

 

 

 

 

 

с)/

 

 

 

 

а/

 

 

 

 

 

 

 

 

Используя очевидное

равенство

 

 

 

 

 

 

 

 

 

 

 

 

 

а

/

 

да12 ч

, —

dw'

да12

д

 

-

 

 

 

 

 

 

 

— д а р ) =

да'дар—--f-^r-—(да,

р),

 

 

 

 

dl

\

 

2

 

 

 

о'.'

 

Т ~ а/

 

 

 

 

 

 

 

 

 

9

 

 

 

г>.<

 

 

 

 

 

 

 

 

 

из двух последних соотношений

получим

 

 

 

 

 

 

 

 

 

да' —

 

(да'дар) :

 

w\2

да р

 

 

и Д 2

а (дар)

(3.72)

 

 

 

 

dl

 

 

~ 2

 

а/

 

 

 

 

 

а/

 

v

 

 

 

 

 

 

 

 

 

С учетом (3.72) результат скалярного умножения уравнения

(3.70) на до' можно представить в следующем

виде:

 

 

 

_а_

щ | 1 2 - \

 

д

/И)12

 

да 12

 

 

 

 

 

 

 

 

 

dt

т р ) + й ( т

 

 

га/

~г

а;

 

 

 

р ш

 

 

 

 

 

 

 

 

 

.

 

,

-

а г ™г

 

 

 

 

 

ар'

 

 

(3.73)

 

 

 

 

 

 

 

 

dl

 

 

 

 

 

 

а/

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

Для преобразования уравнения неразрывности умножим каж­ дый член (3.71) скалярно на р'. Производя операции, аналогич­ ные по форме описанным выше, получим

 

п\2

 

 

 

 

ада ,

In /

 

2 р а 2 +dl \

2 - р й 2

 

 

 

L щ)

i

dt

I

 

 

az

rf/

 

j _ _ a _

 

j _ _ a _

 

 

 

 

 

а 2 а/

 

2 а/ \ р а 2 /

 

а/

 

 

1

/ / /

д Inр

d In

/

= 0 .

(3.74)

 

 

 

136

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

 

а>12

 

д_

р w

-\- pw

W12

,12

dt

~ 2 ~

 

 

2a2p2

 

dl

 

 

 

 

 

d (рда)

 

 

 

 

 

 

 

70 _ap_

• p w

 

Ш)2

Q2/-2

 

Д2

2

 

 

2dl

dl

 

 

 

 

 

 

 

 

 

 

 

' д In p I

d In /

,12

 

 

 

 

 

 

(5/

 

 

 

 

d/

/ « 2

а/ \ -p J

 

pa2

\

Й

T

 

 

 

 

1

й

 

 

 

 

(3.75)

 

 

 

 

 

 

 

 

Т а ? P a 2 / J

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

_а_

-Р I И)12

,12

_а_ р w

-у р W I

 

\- —^

\

X

 

 

dt

 

2р2Д 2

dl

 

V 2

 

2р2д2

/

 

 

Л.

 

 

 

Xdtdtdf

Owм12

й (рда)

 

_ L J L

И)2

Q2/-2N

+

 

 

dl

•W dl — p w

« 2

а/

 

 

 

dа/In р

 

 

 

 

 

 

 

 

d\nf

,12

 

 

— d In /

 

 

 

dl

 

а/

-да -

 

 

ч

pa2

 

 

 

1

д I w

 

 

 

 

dl

 

 

 

 

 

 

 

 

 

 

 

«2

а/

 

 

 

 

 

 

 

При

интегрировании

по объему

необходимо учитывать,

что

в рассматриваемой одномерной по координате постановке пара­ метры потока не зависят от площади сечения f. Интеграл в пра­ вой части (3.76) не приводится к поверхностному. Этот интеграл определяет изменение волновой энергии в рассматриваемом объеме.

Интегрирование в левой части (3.76), поскольку там стоит интеграл от дивергенции, может быть выполнено просто. При этом результат можно представить состоящим из двух частей.

Первая часть

определяет

изменение по времени волновой

энергии

 

 

 

 

12

*2

„12

,12

dfdldt--

 

 

 

 

 

2а2р2

I,

/,

 

 

 

 

 

137

 

 

 

712

-712

dfdl.

(3.77)

 

\ 2

1

2а2р2

( Т + 12д2

р 2

 

 

 

 

 

 

 

 

 

 

Как

известно,

 

_ / да12

,12

 

 

 

величина р

2а2р2

 

носит

название

 

 

 

 

 

 

 

волновой

энергии

единицы объема,

поскольку ^ р

dfdl —

и

кинетическая энергия газа в данный момент времени t, а

и

(•

п12

dfdl — потенциальная энергия сжатия газа.

\

 

•'

2ра2

 

Вторая часть составляет разность потоков волновых энергий через сечения, ограничивающие решетку,

t* и

 

 

,12

dfdldt

=

 

 

2а2р2

 

 

 

 

р w

pw

,12

 

 

2й 2р2

/ 2 -

 

 

 

 

 

р W

/ щД2

-712

fi\dt.

(3.78)

 

2а2р2

 

 

 

 

Каждое из выражений, стоящих в квадратных скобках (3.78), представляет полный поток волновой энергии через сечения со­ ответственно на входе в решетку и на выходе из нее.

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

p'w'-\- pw ( цу\2

,12

(3.79)

 

 

 

2«2р 2

 

Первое слагаемое представляет собой обычный поток акусти­ ческой энергии. Второе слагаемое определяет конвективный пе­ ренос волновой энергии через сечение со скоростью среднего по­ тока w. Впервые формулу для потока волновой энергии среды, скорость среднего движения которой равна w, в виде (3.79) по­ лучил В. Л. Эпштейн.

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

потока энергии, дает

следующие

значения.

Поскольку р'—

= р cos iot

 

 

 

w' — -^-

I Л I I p

] cos (to/

A N ) ,

pa

 

 

 

138

где = — акустическая проводимость, то, обозначая среднюю

за период величину двумя чертами сверху, получим для первого члена

 

 

 

 

 

р w

 

 

 

\Р\

 

 

 

 

(3.80)

где

\р \ —-амплитуда колебания

давления;

 

 

 

 

 

 

б — действительная

часть

акустической

проводимости.

 

Для второго члена будем иметь

 

 

 

 

 

 

 

 

1 £ Е = = = =

,

 

 

1

рда h

I 2 1 Р I2

1

Mw

| ч | 2 | / ? | 2 .

(3.81)

— pww

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

4

 

 

р 2 а 2

4

ра

 

 

 

Учитывая

(3.80) и

(3.81),

окончательно для среднего

за пе­

риод потока

волновой

энергии через сечение будем

иметь

 

 

 

 

 

f\P\2

Н - у

M T O ( l +

h l 2

)

 

 

(3.82)

 

 

 

 

2 pa

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В сечении на входе в решетку

 

поток

энергии

 

определяется

выражением

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/ о ! Л ) 1 2

+

у М с 0 ( 1 +

| Г ] 0 | 2 )

 

 

(3.83)

 

 

 

 

2рс«о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а в сечении на выходе из решетки

 

при среднем числе М

потока

M = M№ 2 sin р2

так:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/ 2 I P 2 I

2

 

^ _ M T O 2 s i n P 2 ( l + | r ) 2 | 2 )

 

(3.84)

 

 

 

 

2 р 2 а 2

82 +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Коэффициент

усиления

по

мощности

решетки

 

определяется

из

(3.69),

(3.83)

и (3.84)

 

 

 

 

 

 

 

 

 

 

 

 

 

2

М,

 

Re (С2) +

 

4 " M w . sin р 2

(1 4-

U 2

I 2 (3.85)

 

Я0

\ Р2«2 /

М Ш 2 sin h

| С 2 1 2

 

 

 

 

 

 

 

 

 

 

 

 

1 + •

 

 

 

 

 

 

 

 

 

 

I Со I 2

2

CO

I 2

 

 

 

 

 

 

 

 

c

 

 

где учтено, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/ 0

P2«2

 

М д а а

sin p2

 

 

 

 

 

Если задан импеданс в каком-либо сечении (£2 или So), то по известным динамическим коэффициентам по (3.46) или (3.47) определяется импеданс в другом сечении (£о или £2) и по (3.48) — модуль передаточной функции по давлению. Таким образом, при заданном импедансе в сечении на входе или на выходе из решет-

139

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