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

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

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

динамические коэффициенты

которой

определяются через А ц

следующим образом:

_

1

^22

Mi

Z12 —

М\

 

 

 

(3.25)

Z1\- — "

 

M\

i l2

 

где AA=AnA22— Al2A2\— определитель передаточной матрицы. Если поставить в выражения для определителя значения Ац,

то легко получить следующее соотношение:

д Д = е т , ' е 7 * ' .

(3.26)

Рассмотрим более подробно динамические коэффициенты Ац. Они не зависят от граничных условий и определяются свойства­ ми рассматриваемой системы. Отметим прежде всего предель­ ные выражения коэффициентов при со->-0:

 

Km

 

 

 

2

 

 

 

• е

 

 

 

I - M 1

 

 

А,.

0;

 

(3.27)

 

h{l)l

 

 

 

А2

• М

1 - М 2

»

 

 

А 22

• 1.

 

 

Для вращающейся решетки коэффициенты Ац зависят от числа М и фактора вращения h(l). Если положить р = а (h(l)=0), то коэффициенты Ац и матрица [А] будут описывать нестационарное течение в элементе неподвижного канала и в канале вращающейся решетки осевой турбомашины. В этом слу­ чае, как нетрудно видеть, будем иметь

 

ik

 

ik

 

У 1 = - -

— ^

У2=-

1 +

м

1

— м

 

Ац — А22

1

т'Че"');

 

1

(3.28)

А;

т'гы)

 

А 21-

1 1

т'1-

 

?1а1

Выражения (3.28) определяют также нестационарное течение в прямой трубе с постоянной скоростью среднего движения (М = =.const).

Наконец, полагая в (3.28) М-УО, получим известные акусти­ ческие зависимости

ПО

An = A22 = cos kl; \

Al2

= ip:a1

si n

An

=

sin

kl;

(3.29)

 

kl.

 

Необходимо отметить одно важное обстоятельство. Если рас­ смотреть решение для вращающейся решетки в несжимаемой жидкости М->0, то из (3.21) и (3.24) видно, что в этом случае \Ч,2 = ± i k и коэффициенты Ац совпадают с (3.29).

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

Выражения (3.29) с соответствующей заменой параметра ра на индуктивность и емкость описывают процесс в электротехнике так называемой длинной линии. Динамические коэффициенты в

общем случае являются комплексными функциями

переменно­

го s. В дальнейшем будем вместо круговой частоты

употреблять

 

<1>

 

величину v =

.

 

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

3.4.2. Стыковка решений на участках

Итак, необходимо связать параметры на выходе из i-ro участ­ ка (обозначим эти параметры индексом 1) с параметрами на входе в ( J +1) участок (индекс 2). Для связи этих параметров необходимо воспользоваться общими уравнениями сохранения в интегральной форме (см., например, [58]):

— неразрывности

 

 

P2w2F2=

 

plw1F1

JL\pdv

.

(3.30)

 

 

 

 

 

 

 

 

ot

J

 

— и энергии

 

 

 

 

 

 

 

v

 

 

 

 

 

 

 

 

 

 

 

 

Р2Щ?2

 

1

 

 

=

Р Л ^

 

 

 

'

2

У.

1

P2

 

 

 

 

 

 

 

д

r

J

«

+

J _

A

]

^

.

( 3.31

 

dt

J

\

2

x — 1

p

 

 

 

V

В связи с тем, что между рассматриваемыми сечениями 1 и 2 объем 1>->0, можно пренебречь интегральными членами в (3.30) и (3.31), которые превращаются в обычные стационарные одно­ мерные уравнения. Отметим также, что в этих интегральных со-

* Точнее это влияние имеет порядок числа М.

Ш

отношениях естественно положить U[=u2. Тогда уравнения не­ разрывности и энергии примут вид

 

 

 

Рч

(3.32)

х —

1 Р!

 

 

 

1 Р2

 

Линеаризуем систему (3.32), после несложных преобразова­

ний получим

 

 

 

 

1

Pi

 

Pi

 

 

Pi

 

P2

Ш)2

 

 

 

 

(3.33)

Pi

' —

P2

' —

 

Pi

-WiW1

— —

\-W2W2.

)

 

P2

 

Система (3.35)

связывает пульсации скоростей и давлений

на выходе из t-го участка с параметрами

на входе в £+1 участок.

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

p2 = A u p 1 J r A n w 1 ;

(3.34)

Если разрешить систему (3.33) относительно р% и w2 и срав­ нить полученные коэффициенты при р\ и W\ с системой (3.34), то получим

 

 

1 —-

Ai

Р2

 

 

1 — м :

 

Pi

Л,2

=

 

Р2

 

 

P ^ i z r

1 - М * .

 

 

 

(3.35)

 

 

 

Pl

 

 

 

 

 

 

 

 

1 —•

 

А2

 

1

fl2 М0

 

 

 

 

P l ^ l

й 1

 

 

 

 

 

1 —

м :

 

Л

=

"2 М 2

 

 

А-22

 

 

 

 

1 — м :

П2

Коэффициенты A-,j пригодны для стыковки решений на от­ дельных участках как во вращающихся, так и в неподвижных каналах. Из системы уравнений (3.35) можно видеть еще одно подтверждение ранее сделанному выводу о том, что решение можно построить только в области дозвуковых скоростей средне­ го течения МШ<С1-

Отметим особо важный частный случай несжимаемой жидко­

сти.

Положим в (3.35)

Mt», —>0, М^, —*0.

Возникающая

при

этом

неопределенность

легко

устраняется,

если заменить М =

— — и положить ai = a2. В этом предельном

случае будем

иметь

а

 

 

 

 

 

 

 

А'и=\;

л ; 2 = Л 2 1 =

0;

j

 

 

 

Wo

}

 

(3.36)

 

 

w l

)

 

 

 

Из (3.36) видно, что отличие нестационарного течения через вращающуюся решетку в несжимаемой жидкости от акустиче­ ских зависимостей проявляется лишь в факторе изменения про­ ходного сечения

=^=-1А.

 

\

Wi

J 2 I

 

Вычисляем, наконец,

определитель матрицы [А]. После

неслож­

ных преобразований

получим

 

 

A ^ ' =

j g . J ^ -

( 1 - M » . ) .

(3.37)

 

P l f l l

MWi

( l - M ^ )

 

3.4.3. Общее решение

Итак, истинная зависимость параметров среднего стационар­ ного течения для получения решения заменена ступенчатой зави­ симостью (рис. 3.1). Необходимо связать между собой парамет­ ры во входном сечении канала (1) с параметрами в выходном сечении (2), при этом в разд. 3.4.1 мы получим решение на участ­

ках

k ' i - \ — k h а в разд. 3.4.2 на

участках

k^ — k i .

Эти решения

выражены

в виде матричных

соотношений. Решению

на

участ­

ках

k'i-ik[ соответствует

матрица [А\.

В

связи

с

тем, что

параметры

на выходе

из участков

k \ - i ~ k i

равны

параметрам

на

входе

в участки

для

связи

входных

параметров

участков

— ki с

выходными

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

ki — k \

надо пергмножить матрицы [А] и [А'].

Такое соединение четырехполюсников между собой называет­ ся в теории линейных цепей каскадным [29]. Если мы длину

113

/ от входа в канал (1) до выхода из него (2) разбили на п уча­ стков, то для связи pi и w{ с р2 и ш)2 необходимо выполнить п перемножений

Рг

^ [ л л ] . . . [ л ; ] [ л 1

]

Р\

w4

 

 

 

 

Обозначая результат перемножения матриц [б], окончатель­ но получим

Pi

"8п

§12

Pi

= [8'] Pi

(3.38)

= _ ^21

 

jw2

§22 _ _W1

-wi

_

Рис. 3.1. Зависимость числа М стационарного течения от коор­ динаты /:

— истинная зависимость:

разбиение на 7 участков; разбиение на 4 участка; разбиение на 2 участка;

— один участок с постоянным числом М

Выражение (3.38) дает решение задачи о нестационарном одномерном течении идеальной сжимаемой жидкости через ка­ налы решеток турбомашины (вращающиеся и неподвижные).

Рассмотрим некоторые свойства полученного решения. Преж­ де всего представляет интерес предельный случай <в-»-0. Как вид­

но

из (3.27), в этом предельном случае решение на

участках

ki-i

kt для неподвижных каналов и каналов осевых

турбома-

шин определяется единичными матрицами, поэтому связь между параметрами на входе и на выходе определяется квазистацио­ нарным решением перемножением матриц [А'], не зависящих от о. Матрицы [А'] определяются по скачкам параметров на от-

114

дельных участках. При го-»-0 произведение отдельных малых скачков соответствует одному общему скачку, определяемому по

параметрам на входе и на

 

 

 

выходе

из

решетки.

Пока­

\А'„

А'П

 

жем

это

на

примере

непод­

 

 

 

вижной

решетки.

 

 

 

 

 

 

 

Рассмотрим

решетку

при

 

 

 

значении

Л1] на

входе

и

М2

м,

 

на

выходе. Разобьем

интер­

 

 

 

вал

 

на

два

участка

(рис.

 

 

 

3.2,

а) и на три

участка

(см.

 

 

 

рис. 3.2, б). Необходимо по­

 

 

 

казать, что при со—>-0 (матри­

 

 

 

ца

[Л] — единичные

матри­

 

 

 

цы)

 

результат

двойного

пе­

 

 

 

ремножения

(б)

совпадает с

 

 

 

матрицей

(а).

 

 

 

 

 

 

 

В

случае

(б),

перемно­ Рис.

3.2.

Схема

перемножения матриц

жая

 

матрицы, будем

иметь

(разбиение

на три участка)

'п

Am

~Ап

А'п

 

AuAu

+ ArtA'iu

АцА~2-{-А12А 22

1\

Ло2

2\

л2 "2

 

 

 

 

Д л я подтверждения сказанного достаточно показать, что

А'и^АиА'п + АпАп;

А22 — А21А12

Л 2 2 Л2 2-

Рассмотрим

первое равенство.

А "

 

А

" ' , будем иметь

Подставляя

явные выражения для

и

 

 

 

 

 

 

 

12

 

 

 

Р2

а

1

? '

\

Д 1

+

А п А п

+

л 1 2 л 2 1 = ^ -

 

 

 

 

^

 

 

 

1 - М 2

 

Р 1

 

( 1 - М 1 2 )

 

I / / Р2 -f- pa —

р'

М ' 2 - М 2 ^ Л

\

а

/

а-2

- М

\

а

 

 

 

а\

 

 

М ' (1 — М | ) piflj

 

 

(1 — м 1 2 )

 

1 — М2, А.

 

 

 

 

2

 

2

 

 

 

 

Р2

\

й

\

 

- А п .

 

1 -

м2

 

 

Pi

2

 

 

 

 

115

Точно так же можно показать и равенство других коэффи­ циентов. Таким образом, в квазистационарном случае ы-*-0 мож­ но не перемножать матрицы [А'] между собой, а прямо опреде­ лять результат, подставляя значения параметров на входе в ре­ шетку и на выходе из нее в систему (3.35).

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

зовании

ЭЦВМ

такое

перемножение

не представляет

осо­

бых трудностей,

желательно ограничить

число

участков

раз­

биения.

 

 

 

 

 

 

 

 

Как

известно,

последовательные

приближения

интегрального

уравнения Вольтерра второго рода

достаточно быстро

сходятся.

В связи

с тем, что наша

задача сведена

к уравнению

Вольтерра

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

Число

М с среднего стационарного течения

изменяется

по ко­

ординате /, как показано на рис. 3.1.

 

 

Этот

сопловой аппарат

рассчитывался в

следующих

вари­

антах.

 

 

М с у + М

 

 

 

1. Число М с постоянно

и равно среднему значению —*^

2. Зависимость МС = М(/) заменена двумя участками равной длины, числа М на которых равны соответственно входному М Г |

ивыходному числам МС з (две ступени).

3.Зависимость МС = М(/) заменена ступеньками с четырьмя участками.

4.Зависимость МС = М(/) заменена ступеньками с семью уча­ стками.

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

Результаты расчета: зависимости действительной Re (йп) и мнимой Im(sji) частей коэффициента суммарной переда­ точной матрицы Sii от частоты v приведены на рис. 3.3. Для

остальных

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

b\j

зависимости

протекают

анало­

гично.

 

 

 

 

при п>4

 

 

 

 

Прежде

всего

отметим,

что

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

переда­

точной матрицы

8/; практически

совпадают.

 

 

 

Важно

также

отметить,

что

при л = 2

качественные

законо­

мерности,

свойственные точному

решению

( п > 4 ) , полностью

повторяются. Количественно значения Re (6')

и Im(6')

при

п = 2

достаточно

близки к

точному

решению.

Поскольку при

числе

участков п — 2 нужны

всего два перемножения, можно

получить

аналитические зависимости

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

передаточной

мат-

116

1,0

Jm(d„)

--

0,1

 

 

0

0,5

-

I

-0,5 I

—•

0 WO 200 300

\

ч ^

--

 

-0,1

V

 

 

 

 

 

 

 

 

 

 

\ \

 

 

-0,2

\

 

 

 

 

 

 

 

 

 

 

 

\

-0,3

 

 

 

\

 

 

 

 

 

 

 

 

 

 

4

 

-0,4

 

 

 

 

\ \ I

 

 

-

 

 

Л

 

 

 

 

 

\

\ \

-0,5 — —

 

 

 

 

 

 

 

 

 

 

 

 

 

 

\

-0,6

 

 

 

 

 

\

 

 

 

 

 

 

 

 

 

 

 

\ к

V

 

 

 

 

 

\ 4

 

 

 

 

 

 

\

 

 

 

400 500 600 700.

\

0,8'

 

 

900 v.ljc

0 100 200 300

~—

400 500

I

1

 

/

 

~г—

 

—i /

 

1

 

I

 

/

/

/

/

 

/

 

/

У

/ '

 

"""

у. У

 

У

600 700 800 900 v1/с

a)

5)

Рис. 3.3. Зависимости действительной Re(6n ) и мнимой

1т(6ц) частей коэффициента суммарной передаточной мат­

рицы бц от частоты v:

система с сосредоточенными

параметрами;

разбиение на 7 участков;

 

разбиение на 4 участка;

 

разбиение на 2 участка; ••— один участок с постоянным числом М

рицы

8,-у. При необходимости этим

обстоятельством

мы

будем

пользоваться. При выборе длины участка разбиения

необходимо

учитывать, что длина участка должна

быть существенно

меньше

длины

волны.

 

 

 

 

 

Наконец, из приведенного сопоставления можно сделать еще

один важный

вывод. Основные трудности решения

заключались

в том, что исходные дифференциальные уравнения были

уравне­

ниями с переменными

коэффициентами, при постоянных по дли­

не средних

параметрах потока

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

случае

M = M l

М 2 = const j

дифференциальные уравнения

превраща­

ются в уравнения с постоянными коэффициентами и их решения достаточно просты. Приведенное сопоставление показывает, что

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

/ не­

обходим: решение при M = const и качественно (при со-»-0)

и ко­

личественно не совпадает с решением, полученным при

учете

изменения средних параметров потока от координат.

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

Рассмотрим два случая.

1. На входе в канал подается возмущение и изучается реак­ ция на это возмущение на выходе.

2. Возмущение подается со стороны выхода из канала и изу­ чается реакция на входе.

Если при одинаковых возмущениях в рассмотренных случаях реакции также одинаковы, то такой четырехполюсник назы­ вается обратимым, если реакции различны — необратимым четы­ рехполюсником. Как указано в [29], критерием обратимости яв­ ляется величина определителя передаточной матрицы

Д § ' = 1 .

Вычислим определитель передаточной матрицы, описываю­ щий нестационарное течение в канале.

Как известно, определитель суммарной матрицы при каскад­ ном соединении равен произведению определителей матриц, вхо­

дящих

в это соединение: Д6' = ЛЛ-АЛ',

с учетом (3.25) и (3.36)

будем

иметь

 

 

 

=

О ^ * 4 » . )

е - '«(т . +т,) .

Очевидно, что Д6== 1 и канал турбомашины является необра­ тимым четырехполюсником. Поэтому при изучении динамических реакций важно различать с какой стороны подается возмущение: со стороны входа или со стороны выхода.

118

3.5. ОДНОМЕРНОЕ НЕСТАЦИОНАРНОЕ ТЕЧЕНИЕ В РЕШЕТКАХ

Необходимо отыскать связь между нестационарными значе­ ниями параметров в абсолютном и относительном движениях и тем самым перейти от течения в канале к течению в решетке. Для неподвижных решеток, обтекаемых стационарным потоком без угла атаки, полученное в предыдущем параграфе решение формальной заменой М Ю | на Мс , « M.W2 на МС 2 при h(l) = 0 описы­ вает нестационарное течение в решетке. Перейдем к построению решения для вращающихся решеток.

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

жесткими, линеаризуя соотношение C2a = ^ 2 s i n

62, получим

c2a = w2 sin fS2-

(3.39)

При этом можно принять, что пульсации статического давле­ ния в абсолютном и относительном движениях совпадают. Таким образом, связь между нестационарными параметрами в абсолют- v ном и относительном движениях в выходном сечении вращаю­

щейся решетки записывается так:

 

"1

0

Pi

О

sin 82

Wo

или

 

 

(3.40)

 

 

 

 

Pi

Pi

 

 

С

1™1

 

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

Свяжем теперь параметры на входе в решетку, т. е. в зазоре между сопловым аппаратом и рабочим колесом (обозначим это сечение индексом «0») с параметрами во входном сечении вра­ щающегося канала (индекс «1»). В абсолютном движении пуль­ сации параметров в зазоре определяются значениями Ро и с0 - Скорость с0 в этом сечении связана со скоростью WQ очевидным соотношением

2

2

(3.41)

WU =

CQ-{- и — 2c0u cos а0 .

 

Поскольку лопатки соплового аппарата считаются жесткими,, при нестационарном течении принимается ao=const.

Так как расстояние между сечениями «0» и «1» невелико, для связи параметров в этих сечениях можно использовать урав­ нения сохранения массы и энергии в виде (3.32):

р0 с0 sin a0 r= plwl sin

(3.42)

119

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