Файл: Филиппов, А. П. Воздействие динамических нагрузок на элементы конструкций.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 23.10.2024

Просмотров: 72

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

нужденные колебания определены в виде разложения по собствен­ ным формам системы Бал — лопатки, найденным методом динами­ ческих жесткостей1 .

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

Рис. 38.

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

Указанный метод содержит следующие принципиальные по­ ложения:

1.Реальный валопровод заменяется дискретной, по отношению

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

гофа — Клебша.

2. Уравнения движения расчетной модели получены на осно­ вании обобщенного принципа Остроградского— Гамильтона. Внут­ реннее и внешнее демпфирование учитывается по Фойгту и по данным работы [41 ].

3. Решение уравнений вынужденных колебаний при нулевых начальных условиях получено в виде разложения по собственным формам колебаний системы вал — диск — лопатки с помощью преобразования Лапласа.

4. Собственные частоты и формы колебаний системы вал — диск — лопатки определяются методом, основанным на разложе­ нии форм собственных колебаний по формам элементов системы (лопаток, дисков, валов); частотное уравнение при этом не содер­ жит «антирезонансов», характерных для разветвленных систем.

Приведенные ниже результаты относятся к расчетным исследо­ ваниям турбоагрегатов большой мощности. Программа, реализую­ щая метод, составлена на языке программирования «АЛГОЛ — 60» Представляя бочку ротора генератора одной массой, на которую действует сосредоточенный крутящий момент к. з., приведем рас­

четную схему к виду, условно показанному на рис. 38.

1 Расчеты выполнены на ЭЦВМ «БЭСМ-6».

139

Здесь /,• и С[ — суммарный момент инерции участка и соответст­ вующая жесткость; /В1 — момент инерции участка вала под дис­ ком; М (t) — крутящий электродинамический момент к. з.

Выражение для крутящего электродинамического момента в наиболее общем случае двухфазного к. з. может быть представлено в следующем виде [6];

М (t) = ^ 2

А{е~к1^ sinotf -+-

 

/5

/

6

 

+ 2

В 1е~¥

sin 2соt + 2 D fi~kf .

(4.1)

\i=0

'

1=0

 

где со — круговая частота тока в сети. В случае трехфазного ко­ роткого замыкания в выражении (4.1) отсутствует вторая гармо­ ника, а значения амплитуд At и Dh а также коэффициенты kit k-L вычисляются по другим формулам.

Равновесным состоянием системы является вращение с постоян­ ной угловой скоростью со. Предполагается, что при совместных колебаниях диски совершают крутильные движения в своей плос­ кости, а лопатки — пространственные изгибные колебания в поле центробежных сил. .Влиянием центробежных и кориолисовых сил на колебания диска будем пренебрегать, что допустимо для дисков паровых турбин [50]. Лопатки трактуются как закрученные стерж­ ни. Предполагается совпадение центра тяжести и центра изгиба сечения лопатки. Сдвиг, инерция поворота и депланация попереч­ ного сечения, а также влияние кориолисовых сил не учитывается.

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

для общего случая, когда все участки системы имеют облопаченные диски и ко всем массам приложены внешние моменты М ( (t), получены на основании вариационного принципа Остроградского — Гамильтона и могут быть приведены к следующему виду

д

 

h + ~ sr)

(/в‘'9' +

л г ) )

at И р3хМ р

d

t

mt \fit- (р) ру* dp

+

U + х

d t

X

X

{ -

C ;_i ( 0 ,_ ! -

0 ()

+

С, (0 , -

0 ц .,)} = M t (t),

~W

{h +

~ sr) (2lt ~T a‘ (p) P3^* — m№i (p) py*)j +

+

( 1 + X4

- ) ( ^ ^

- «

l-(p)P3G ^ f ') = 0 ,

 

 

 

 

 

 

 

 

d p j

 

 

h +

- J - )

l*( (р) Уi

 

 

d

d

t

 

 

 

140


 

d2zt

d

к и

dyt

 

Q 2^,- ( p ) у i

=

o,

 

+ E l f

dp

„p

 

 

~W

o'p2

^

 

E l yz

 

(*+■4 fti (p) 2i | +

^ +

dp2

dp'1

d

 

 

 

 

 

 

 

(t)

&yi

 

 

 

 

dzi

=

0

 

 

 

+ £/ ® - 0

dp

N i (p )

5p

 

 

 

(i =

1, 2, 3,

, П ) .

 

 

 

 

 

 

 

Здесь

 

 

 

 

 

 

 

 

(4.2)

 

 

 

 

 

 

 

 

 

 

— Gi +

tyi, У* = — 9;P + ФгР +

Do

2* = zt (4.3)

— перемещения элементов системы в абсолютном движении (пере­ мещения диска и лопаток показаны в двух проекциях на рис. 39);

0; = 0 (I, 0 .

Фг = Ф; (Р> *),

У1 =УА Р. 0 .

= zt (Р» 0 — пе­

ремещения элементов системы

в относительном движении;

р— текущий радиус (г0 < р < < R)\ t — время; at (р) —

толщина диска;

p,f

(р)

и i f ,

 

 

 

 

 

 

 

 

lyl, i f

— погонная

масса

и

 

 

 

 

 

 

 

 

моменты инерции сечений ло­

 

 

 

 

 

 

 

 

паток; N( (р) — сила инерции

 

 

 

 

 

 

 

 

части

лопатки

над сечением

 

 

 

 

 

 

 

 

р; Q — угловая

скорость вра­

на диске; у и

 

 

 

 

 

щения;

т (- — число

лопаток

G — удельный

вес

и модуль

сдвига

 

материала

диска;

g — ускорение

свободного

падения;

Е — модуль

упругости

материала лопаток;

X и

h

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

демпфирования

при

рассеянии

энергии

колеба­

ний в

материале элементов и

во

внешнюю среду.

Заметим,

что

Рг(Р) =

i f

= lyl

=

i f

= 0

при р < г и а-г (р) =

0

при р >

г0.

На

основании

принципа

Остроградского — Гамильтона полу­

чены также и естественные граничные условия для функций ф4,

yit zu которые используются при решении уравнений

(4.2):

 

 

_а_

 

 

 

JN H |й =

0

,

 

8ф<

1 + х

dt

\ 2 n G d i ( р ) р 3 -

/ к

 

 

 

d_

 

W

d2yi

 

 

 

fyi

 

 

 

EI

+

V)

 

 

 

1 +

х -

 

. др

 

Ф 2

EI yz

 

dp2

 

 

 

 

 

 

 

 

 

 

 

■Nt (Р)

dyt

 

= 0

,

 

 

 

 

 

 

 

dp

 

 

 

 

 

 

 

 

bzt

1 - f

X-

 

a

E l f

&Vl

+

W ~ w

(4.4)

 

ap

ap2

Ni (P)

QZj

 

= 0,

 

 

 

 

 

 

 

dp

 

 

 

 

 

 

 

 

141



dyt

1 +

X

d

(O

dp2

+ e i ;

а*г« '

= 0,

dp

dt

EI\

aP2 i

dzi

1 +

X

d

EI%

 

 

6 2 Z ;

= 0.

dp

dt

-i f

- + £ / :

ap2

Уравнения (4.2) описывают колебания диссипативной системы, которые, как известно, в общем случае не допускают решения в виде разложения по собственным формам. Для неконсервативных систем метод разложения по собственным формам применим к си­ стемам, допускающим разделение переменных. Существуют при­ ближенные методы разделения переменных для диссипативных си­ стем, а также вообще для неконсервативных систем [10].

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

Уравнения (4.2), описывающие малые колебания дискретно­ континуальной системы, учитывает рассеяние энергии в материале в соответствии с моделью Рэлея — Фойгта. Здесь принято, что внешнее демпфирование пропорционально коэффициенту h, а внут­ реннее — коэффициенту X в любом элементе системы.

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

[41].

Решение уравнений вынужденных колебаний

Применив метод Фурье разделения переменных, представим ре­ шение уравнения (4.2)' в виде

 

6» — 2

Тft ij,

yt =

2 Tjijij,

^

 

to = 2

 

zt

=2

TjZij.

 

Здесь

Tj = Tj (t) — неизвестные

функции

времени; 0;;- =

0,- (г),

%■ =

Фц (р), Уу Уч (р),

Zij =

Z4 (р) — формы колебаний,

соот­

ветствующие абсолютным перемещениям вала и относительным перемещениям диска и лопаток.

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

вим решение (4.5) в (4.2) и

введем изображения функций Tj (t)

и Mt (0

 

«t (р) ■» Г/ (0;

u i (р ) - ^ м с(0 -


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

0£е, tyie, Hie. z£e,

интегрируя

три

последние

уравнения

по р

от

г0 до R и суммируя все уравнения по г, приводим их к следующему

виду:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

П

 

Щ (Р)

,

 

hp) h .fiifiie +

2 п -L-

с

 

 

,

,

 

2

2

(Р2 +

\ а,- (р) р3ф1/ф(еФ +

 

г=!

/

 

1

 

 

 

 

 

g

i

 

 

 

 

 

+ trii | Pi (p) (УцУи + ZijZie) dp

+

(1 + ф ) |c£ (0£/— 0;+l./) (0<e —

 

Q i + i . e )

+ 2n,G J

at (p) P3

dp

dp

 

dp + mt J

F ,

r

d2yg-

diyie

,

 

*

J „ 9 .

dp2

 

 

 

 

 

 

 

 

 

 

 

 

dp2

 

 

 

I

Cl

I

d%i

d%ei't/ie

,

d2yq

 

d%zie

 

1 i

 

 

 

 

”4 .

 

"ls

I

•/‘/ .

dp2

'

 

 

 

 

 

 

,-/г \

dp2

dp2

^

dp2

 

^

 

 

+ E Iy

dza

+dp

d2zn

—~ ' dp2

dzie

d p

dp

d%„ 1 .

R

'c . . . .

-------- ^dp -f- mi j N£(p)

dp2 i

1

Jro‘м

R

m ( Q 2 f p £ (p ) y i j y i e d

i

( e = l . 2, . . . , л).

I

dyu

dy,

+

 

?t;

“Vie

V dp

dp

 

p

= H ( /

i (p) 0(1

 

г=1

 

 

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

а,- (р) =

1

R i^ V iiP )

—=z-

 

sf + (и/ + Р)2

Здесь 2«/ = /г + xsyJ S/

= si ~

nh si ~- собственная круговая час-

тота системы вал — диск — лопатки; щ — приведенный коэффициент демпфирования;

' П (

/„А2/+ 2я - М

Г

Rj = 2

й£ (р) р3 (ф;/)2 ф +

Ь=1

ё

г„

143