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

книги из ГПНТБ / Богомолов С.И. Взаимосвязанные колебания в турбомашинах и газотурбинных двигателях

.pdf
Скачиваний:
23
Добавлен:
23.10.2023
Размер:
6.31 Mб
Скачать

I I

 

в а р и а н т .

 

На торцах

участка

оболочки

должны

быть

выполнены

условия

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W

 

0.

 

 

 

 

 

 

 

 

у _

 

 

у

 

_

 

 

 

 

 

 

 

 

 

 

 

 

 

 

После преобразований, аналогичных выполненным ра­

нее,

получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-1-1

 

 

 

 

 

 

(114)

При

этом

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а-> о 2

 

 

 

 

 

 

 

 

 

 

 

г

 

ей у

 

 

 

I I I

а2

=

ац

 

 

 

 

 

 

 

 

 

 

в а—•

С1о[ •а2 .ійз.,1 а:зі,

 

d2

 

 

 

 

 

 

 

С2

 

 

 

 

 

 

 

 

 

 

р и а н т .

На

торцах

участка

оболочки

должны

быть

выполнены

условия

 

 

 

 

 

 

 

 

 

 

 

 

 

W

 

Pu

 

 

0.

 

 

 

 

 

 

 

 

у _ і

Pz_

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Соответствующее

матричное

уравнение четвертого по­

рядка

может

быть

записано

в

 

форме

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(115)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Rl

 

a3

Ô3

 

 

 

 

 

 

 

 

 

 

Сз

ö?3.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а з

=

« u

 

 

 

 

Ьз

= Й12

• a i . a 4 2 ;

 

 

 

 

 

 

 

 

О 4 1 !

 

 

 

 

 

-i

 

 

 

cl

=

a2\

Й24Й44

 

^3

= Û22

Û24Û44 Û42-

 

Матричное

уравнение

колебаний

системы облопаченных

дисков,

связанных

в

пакет

оболочкой

вращения.

Запишем

уравнения колебании конструктивных элементов ротора.

Матричное уравнение колебаний /-го диска

пакета

имеет вид

 

X.*, , =П.Х'..

(116)

Матричное уравнение колебаний участка оболочки между двумя дисками или оболочек — крайних участков ротора

НО

при жестком защемлении или шарнирном оппрашш (двух типов) на концах пакета запишется как

х ^ - е д .

(П7)

Матричное уравнение колебаний левого крайнего

участка

оболочки при свободном

левом конце пакета

 

х ; +

| = в д .

( l i s )

Матричное уравнение колебаний правого крайнего участ­

ка оболочки при

свободном правом конце

пакета

 

х ; + , =

ад.

( и з )

Рекуррентные

матричные

соотношения

(116) — (119)

позволяют записать матричное уравнение колебаний па­

кета дисков как единого целого.

Это

уравнение связы­

вает параметры

напряженного и

деформированного со­

стояния в крайних сечениях оболочек

пакета.

Если пакет

из /( дисков начинается

и заканчивается

участками оболочки со свободными торцами (рис. 37),

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

целого запи­

сывается следующим образом:

 

или

 

x ; . H = s x ;

(120)

где

 

При других рассмотренных выше граничных условиях на торцах ротора, начинающегося и заканчивающегося участком оболочки, матрица 5 будет иметь вид

S = К(к+і)ПкР

• • • n ^ u

• • • ^ * 2 П і ВН­

ЕСЛИ пакет

дисков

начинается

п заканчивается дис­

ком, матрица

S имеет

структуру

 

 

 

5

=

 

вд^...п,я;,...ед.

 

Если

пакет

дисков

начинается

участком

оболочки и

заканчивается

диском,

формирование матрицы 5 выпол­

няется аналогично.

 

 

 

 

Таким

образом,

разработанный

алгоритм

позволяет

получить

матричное

уравнение колебаний пакета дисков,

m

который может как начинаться, так и заканчиваться участком оболочки либо диском.

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

летворяем граничным условиям на торцах

крайних обо­

лочек. Если

пакет начинается

или заканчивается

диском,

необходимо

ввести фиктивный

радиус (гх—г^

или

г 2 = г 2 ф )

и удовлетворить граничным условиям на окружности последнего. При этом

''о < Пф < Гд, Го < Оф < ГА.

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

Запишем уравнение (120) следующим образом:

Рх

 

^12

s13

S,4

Px

ми

S2 1

S,,

Sr.

5 2 4

My

и

 

 

s-M

 

и

d\V

 

s,.

si3

 

dW

dx

' - H Su

Su

dx

Подчиняя это уравнение граничным условиям на ле­ вом и правом концах ротора

Рх

I Рх

= о,

My

M

У і

получаем систему линейных алгебраических уравнений

 

dW,

 

 

dx '

 

 

dW,

(121)

0 =

ЗД+S24^.

Условия нетривиальности решения системы (121) дает частотное уравнение пакета

 

Д ( Р )

•43 °14

_ Q

 

 

Sv-J S24

 

 

 

 

 

 

Вид

частотного уравнения в зависимости

от граничных

условий

на концах ротора показан

в табл.

31.

112

 

 

Т а б л и ц а 31

 

Правый край

 

Свобод­

Защемлен

Оперт со

Шарннрно

скользящим

ный

оперт

 

шарниром

Левый край

 

 

Оболочка

=0

'' + 1

о

s

Свободный

Рх

= 0 Рх

МУ.

Защемлен

иU 1

d\V = 0 dU = 0

dx J1

dr . i

Оперт со скользя­ щим шарниром

=0

Шарнирно оперт

M,

= 0 M y

U

и It

Г Р..

 

 

 

=

0

 

_ '+1

•Si3

5 1

4

S23

S24

S3

S 1 4

S23

S24

Sl3

5j4

 

•S23

S24

S21

S24

 

-

и

-

 

 

 

dW

=

Ü

= 0

.

dx

-

- М У .

'•+1

 

 

 

Диск

 

-

 

 

 

dU

 

 

 

_ dr

_

 

 

 

S33

S 3 4

^13

^14

 

 

 

S23

S 24

•Ss3 ^34

•^13

^14

Sj 3

S44

S23

S24

S33

S34

•S]3

Sj4

S 4 3

S 1 4

S 2 3

S24

 

 

Su

Su

S 4 1

SA

*^2l

"^24

~Мѵ~у = 0

и' + 1

=0

(7

5гд S24

•^ЗЗ S34

S23 S24

•^33 S34

•^23 ^24

•^33 S34

*S°1 5g4

•Ssi S34

113

§ 2. СПЕКТР ЧАСТОТ И ФОРМ КОЛЕБАНИЙ ПАКЕТА ДИСКОВ

Спектр

частот

и

форм

колебаний

пакета

 

одинаковых

дисков.

Пакет

из одинаковых дисков,

хотя

и

не

встреча­

ется в практике, позволяет

качественно

проанализировать

особенности спектра частот пакета.

 

 

 

 

 

 

 

 

 

Одним из существенных факторов, определяющих

появление

взаимосвязанных

колебаний

дисков

в

пакете,

 

 

 

 

пn- ~ — —n

 

является

геометрия

уча­

 

 

 

 

 

стка

оболочки

между

 

 

 

 

 

дисками.

 

Немаловаж­

 

 

 

 

 

 

Y

 

ную

роль

играет

при

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

этом т а к ж е

соотношение

 

 

 

 

 

 

 

 

между упругими

и инер­

 

 

 

 

 

 

 

 

ционными

 

свойствами

 

 

 

 

 

It

 

участков

оболочки

и

 

 

 

 

 

 

 

 

дисков.

 

Поэтому

гео­

 

 

 

 

 

 

 

 

метрические

 

 

размеры

 

 

 

 

 

 

 

 

этого

пакета

соответст­

 

 

 

 

 

 

 

 

вуют

размерам

рабоче­

 

 

 

 

 

 

 

 

го

колеса

компрессора

 

 

 

 

 

 

 

 

одного

из

ГТД (рис.

 

 

 

 

 

 

 

 

41, я): l t

= / я =0,0835 м;

 

 

 

 

 

 

 

 

S = 0,45

1 0 - 2 л ( ,

г 1

=

 

 

 

 

 

 

 

 

=

Го =

 

17,6

10--' м,

 

 

 

 

 

 

 

 

гА

="22,95 • Ю - 2

л/,

г0

=

Рис. 41. Рабочее колесо

компрессора

=

1,25-

Ю - 2 м.

Измене­

Г Т Д

и пакеты

для двух

и

пяти

оди­

ние

толщины

диска

 

наковых

рабочих

колес.

 

вдоль

 

радиуса

приве­

 

 

 

 

 

 

 

 

дено в табл . 32. Модуль

упругости

материала

колеса /5 =

1,96-

101 1

н/м2,

плот­

ность

р =

7,84 • 103

н •

сек2/м*.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а

32

 

г • ІО2, м

0

2,50

3,10

17,00 17,38

17,53 18,20

21,10

21,85

21,85

h • 102,

м

1,15

1,15

0,42

0,42

1,2

1,2

0,50

0,50

4,20

4,20

Поскольку основное внимание уделялось выявлению особенностей, которые вносят в колебания системы обо­ лочкой вращения, ступень рассматривалась без лопаток. Формирование пакета соответственно из двух и пяти дисков показано на рис. 41, б, в.

114

Д л я анализа спектра

частот

и форы колебаний па­

кета был определен спектр

частот

отдельного колеса при

граничных условиях свободных торцов оболочки, не

сопрягающихся с

диском,

 

 

 

 

 

 

 

 

 

 

 

Рх

= M,, = Ру

= Рг

= О

 

 

 

и условиях скользящего

шарнира

 

 

 

 

 

 

 

 

Рх

 

= ми = W = V = 0.

 

 

 

Спектр

низших

 

частот

колебаний

ступени

при

двух

узловых

диаметрах

 

на диске

приведен в табл.

33.

 

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а

33

Номер

частоты

Свободные

кон­

 

Скользящий

 

цы

оболочки

 

шарнир

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

593

 

 

 

 

 

 

 

 

 

2

 

 

 

 

1261

 

1164

 

 

 

 

 

3

 

 

 

1876

 

1383

 

 

Табл.

 

33

показывает,

что спектр

частот

колеса

(час­

тота дана

в

герцах) изменяется с изменением

граничных

условий

на торцах

 

оболочки. Это означает,

что участки

оболочки

 

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

системы

 

являются

короткими

и при включении данного рабочего

колеса

в

пакет его

колебания

будут

 

зависеть

от упругих и

инерционных

свойств

дисков,

присоединенных

к

торцам

оболочки.

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

Д л я

сравнения спектра частот колеса с парциальными

частотами

его

отдельных

элементов

были

определены

частоты

внешней

 

и внутренней

частей

диска,

жестко

закрепленных

по

 

окружности

 

радиуса

барабана .

Дл я

наружной

части

диска

(гА = 0,2295

м,

г0 = 0,1760 м)

при

п=2

 

низшие

частоты

^ = 1980 гц;

f%

= 9600

гц.

Д л я

внутренней

части,

жестко закрепленной

по

наруж ­

ному

контуру

(гА = 0,1760 м;

г 0

= 0,0125

м),

 

 

 

 

 

 

 

! х

=

1340 гц,

/ 2

=

3220

гц.

 

 

 

 

115

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

Данные расчетов колеса

позволяют

проанализировать

спектр

частот

пакета

из

двух одинаковых дисков. На

рис. 42, б представлены

результаты

расчета низших

частот

пакета

из двух

дисков с условиями скользящего

/,-№и>

/,,-тзгц

Л'/Wen

 

 

/2-Жга

^/гг/Шгц

і

/ а - l W i a

Рис.

42. Низшие

частоты н формы колебаний одиночного

диска,

подкрепленного оболочками вращения

н пакета

из двух

 

 

дисков.

 

 

шарнира на концах оболочки. Сравнение спектра частот пакета со спектром отдельного колеса показывает, что в пакете происходит расслоение частот колеса. Частоте колеса 1164 гц в пакете соответствуют две близкие час­ тоты 1163 и 1170 гц с различной ориентацией форм коле­ баний. Соответственно частота колеса 1383 гц расслаива­ ется здесь на частоты 1370 и 1410 гц. На рис. 43 приведен спектр частот пакета нз пяти одинаковых дисков, соот­ ветствующий частоте колеса 1383 гц.

 

В реальных

пакетах,

облопаченные

диски и

участки

подкрепляющей оболочки которых имеют

разные

упругие

и

инерционные

свойства,

такой строгой

закономерности

не

наблюдается,

однако

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

116

место.

Определить

число частот

в спектре

реального

пакета

для

рассматриваемого диапазона, а

также оце­

нить

их

качественно

 

 

 

можно

лишь

непосред­

 

і

 

ственно путем

расчета

 

 

 

 

 

взаимосвязанных

коле­

)

 

 

баний

дисков

в

па­

1

 

кете.

 

 

 

 

 

лебания пакета дисков сложной конструкции.

Как уже отмечалось, роторы некоторых ГТД представляют собой вал с укрепленным на нем

диском, который на подкреплягощей оболочке вращения несет систему облопаченных дисков (рис. 1, в). Рассмотрим результаты расчета вза­ имосвязанных колеба­ ний дисков ротора ГТД

такой конструкции. Расчетная схема рас­

сматриваемого пакета представлена на рис. 44. Пакет состоит из двух облопаченных и одного необлопаченного дисков. Необлопаченный диск непосред­ ственно связан с валом. Условиями закрепле­ ния его внутреннего контура являются усло­ вия закрепления си­ стемы на валу. В качестве последних при-

няты

условия жестко­

го

закрепления

на

линейное

перемещение

упругого

на угловое.

 

I

 

f!2 'fflßjW

1

 

 

 

 

t)

Ja-ПЮгь

 

 

14

/ 1

1 к _

J Р и с_ 43.. Спектр частот и форм коле- б а н и й пакета из пяти дисков, соот­

ветствующих колебаниям рабочего колеса с частотой 1383 гц.

внутреннего контура диска и

117

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

 

 

 

 

 

 

 

 

 

X i - f i =

SXL,

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

5

= П э / ? п і П 3 / ? ц П 1 .

 

 

 

 

 

 

 

П р и

этом

r i j ,

П у

— матрицы

перехода

через

облопачеи-

 

 

 

 

 

 

 

 

ные диски;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

П а

— матрица

перехода

через

пеоблопачен-

 

 

Ru.

 

 

 

 

ный диск;

 

 

 

 

 

 

 

 

 

 

 

 

 

Л I I L

 

м а т р и ц ы

участков

 

цилиндрической

 

 

 

 

 

 

 

оболочки, расположенных"между дис­

 

 

 

 

 

 

 

 

ками.

 

 

 

 

 

 

 

 

 

 

 

Исходными

для

расчета

с л у ж и л и

реальные

размеры

 

 

 

 

 

 

 

 

 

одной

из

конструкций

 

по­

 

 

 

 

 

 

 

 

 

добного типа. Радиусы внут­

 

 

 

 

 

 

 

 

 

р е н н и х

контуров

облопачен-

 

 

 

 

 

 

 

 

 

ных

 

дисков

 

составляют

 

 

 

 

 

 

 

 

 

0,07 м,

 

необлопаченного —

 

 

 

 

 

 

 

 

 

0,06

м.

Радиус

наружного

 

 

 

 

 

 

 

 

 

контура

первого

диска —

 

 

 

 

 

 

 

 

 

0,38

м.

На

диске

закрепле­

 

 

 

 

 

1!

 

 

ны

лопатки

длиной

0,25

м.

 

 

 

 

 

 

 

 

 

Р а д и у с наружного

контура

 

 

 

 

 

 

 

 

 

 

третьего

диска

— 0,37 м. Н а

 

 

 

 

 

 

 

 

 

 

диске

установлены

лопатки

 

 

 

 

 

 

 

 

 

 

длиной

0,295 м. Радиус

на­

 

 

 

 

 

 

 

 

 

 

ружного

контура

необлопа­

Рис.

 

44.

 

Расчетная

 

схема

ченного

диска

совпадает

с

пакета

дисков

ротора

ГТД .

радиусом

 

подкрепляющей

 

м.

 

 

 

 

 

 

 

оболочки

 

и

 

составляет

0,175

Все диски

являются

круглыми

пластинами

пере­

менной вдоль радиуса толщины. Лопатки

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

собой

закрученные

стержни переменного

поперечного

сечения.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

На рис. 45 приведены частоты (при я =

0)

и

соответ­

ствующие

формы

колебаний

облопаченных

дисков,

жест­

ко закрепленных по окружности

радиуса

подкрепляющей

оболочки

 

(частоты

и

формы

парциальных

систем),

а

на

рис.

46 — частоты

и

формы

взаимосвязанных

колебаний

системы,

когда

диски совершают

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

коле­

бания

(?г =

0).

Сравнительный

анализ

частот

и

форм

колебаний показывает, что спектр частот

взаимосвязан­

ных

колебаний

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

узла

реального

ротора

118

качественно н количественно отличается от спектра

парциальных

частот

системы.

 

 

(п >

 

 

При веерных

колебаниях

дисков

2)

спектр час­

тот взаимосвязанных

колебаний

состоит

в

основном из

парциальных

частот.

Вместе

с

тем, кроме

этих частот,

 

a

at/at

 

 

,f

foc g

 

 

Рис. 45. Спектр частот и форм колебаний отдельных облопаченных дисков Г Т Д при /і = 0.

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

§ 3. ИСПОЛЬЗОВАНИЕ АЛГОРИТМА РАСЧЕТА ПАКЕТА ДИСКОВ ДЛЯ РАСЧЕТА СОБСТВЕННЫХ КОЛЕБАНИЙ ОДИНОЧНЫХ ДИСКОВ

Диск, подкрепленный

оболочками с

двух сторон. Алго­

ритм расчета пакета

облопаченных

дисков позволяет

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

119

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