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

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

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

X

sin гj (\] — b) +

cos /■i (Г) — b)j — [zn (b) — zr (6)] X

 

X ~ Г +

1 jsinrj. (ri — b)

Применяя в частных решениях (3.61) теорему о среднем значении, выносим F (л) за знак интеграла. На основании (3.56) имеем

с р

( \ + р

ср)

После интегрирования частные решения принимают следующий вид:

 

2п.ч —■

1 + Р -

ср

2а2Р

 

(Л —b f

 

 

 

 

 

 

 

 

 

 

р +

1

 

2/2

 

 

 

 

 

 

 

г2

1 +

 

№■

[l

— ^ (-y -S i +

Cj)

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 + p -

 

 

j t 2

Г

(л — 6)

«А 1

 

 

Тх

 

 

 

 

 

 

P + 1

 

2а2Р [|

 

Ir

 

} (3-63)

 

 

 

 

 

1 + P —Л;р

 

It 2

 

— bf

 

 

 

 

Zpu

 

 

 

+

 

 

 

 

 

P +

1

20^“

 

 

 

 

 

 

 

 

 

 

2i2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

?1 ( - ~ sl +

C1

 

 

 

 

 

1

 

^Г.Ч

 

 

1 + p-

cp

 

я 2

Г (Л — b)

PeA ,

 

где

 

 

 

dr|

 

 

 

P +

1

 

 

L'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r =

i

 

 

и

и

 

 

л2С (р -j- 1)

.

 

 

"I /”,

/х2

,

ril =

 

rhk,

k =

У - 2a2g—

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2a2|3p

 

 

 

 

 

 

 

h = h1l —

п Ч С ( р +

1)

c

=

 

 

z =

Pi11

 

 

 

 

2a2pp

 

 

 

2cnZ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|(3.64)

 

ft2

 

 

 

 

 

 

 

 

 

-A(ri—6)

 

 

 

j t 2Z 2C

( p +

1)

 

ех =

е

 

 

 

 

ra

A2kz

 

 

 

2a 2Pp

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Sx =

Sin

г (»} — 6)

f

Cj — cos

Г(Л — Ь)

 

 

 

 

 

v1l

'

 

 

^

 

 

 

 

!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Разбив интервал интегрирования на п участков длиной т = 1/п, положим b = тх, г) = т т + т. Введем обозначения

Q (т) = zr (тх) + pzn (тх),

1

dzr (Т])

i P

dzu (л)

п

dr\

dr\

2П(тт) — 2Г(mx),

 

1

dzn (тр

 

dzr (л)

k

dr\

 

dr\

Ч = т т

(3.65)

r ) = m x

132

Подставив в (3.65) общие (3.62) и частные (3.61) решения, получим рекуррентные формулы:

Q (т + 1) = Q (т) + R (т) +

+

4 а 2|3л2 ^

^ — Рщ+\)

 

 

R (т + 1) =

R (т) Н--- 2а2fjn2 ^

Р Рm+i).

S + 1) = 5 (т) е1

sx + сх) + Т (пг) е& —

 

 

 

 

( h

(3.66)

 

2

 

 

 

 

 

1 ~ е1 I ~ S1 + С1

 

 

Щ Г ^ + Р P m+\)

 

Ji

 

т (m + 1) =

Т (т) ег (сх — -4-Sj) —

 

-

S (т)

------^ + Р -

 

“7 ^

где Pm+1 — среднее значение силы Р (л)

при т т

< л < т т + т,

ех = е-А/л,

sx = sin r/n, q = cos

гIn.

 

 

Перемещения и скорости подрессоренной и неподрессоренной

масс будут

 

 

 

 

 

гп (тт)

=

Q (m) + S (m)

 

P + 1

 

 

 

 

 

 

 

гГ (тт)

=

Q (m ) p S

(m)

 

 

 

 

P + 1

 

 

1

&п

к

K (m ) + T (m )

(3.67)

 

 

Tj

dr)

r\—mx

p + I

 

 

 

 

~т~ R (m) — p7 (rn)

 

 

 

 

 

k

 

'"i

dr]

r ] = m T

P+1

 

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

ИЗОГНУТОЙ ОСИ баЛКИ Примет ВИД 2Г (т т + т) = Z (Г), ll)|4=mi:+T

или в обозначениях (3.66) и (3.67)

Q ( m + i) - p s ( m + i)

= %

M {i

C T + 1 )sin J i E(£ ?_+1)„,

р + 1

-Т,

v ’

1 '

п

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

133

К о л е б а н и я п л а с т и н

Уравнение колебаний нагруженной прямоугольной плиты, опер той по всему контуру (рис. 37) имеет вид

д%)

.

d*w

рh

d2w

 

q(x,

y, t)

(3.68)

~дх*

+ 2 дх2ду2

 

ТТ ~dF

 

 

D

 

 

 

 

 

 

где w (x, y,

t) — отклонение точки (x,

у) от положения равнове­

сия; рh — масса плиты,

отнесенная

к

единице

 

поверхности;

 

 

D =

Eh3

а

интенсивность

 

 

ттг-г,------я ;

 

 

 

12(1—-v)2 ’

 

 

 

 

 

 

внешней нагрузки.

 

 

 

 

 

 

Решение уравнения (3.68) записы­

 

 

ваем в виде ряда

 

 

 

 

 

 

 

 

СО

с о

 

 

 

тпх

 

 

W

 

2

 

 

 

 

 

 

 

2 ф/nn (t) sin „ X

 

 

 

 

m=1/i=l

 

 

 

 

 

 

 

 

X

sin

плу

 

 

(3.69)

Рис. 37.

 

? - 2

2

т

п

где w0

=

4Р 0р1а 2/я4О , pj =

a/b, P 0

вес груза,

перемещающегося по плас­

тинке.

Разлагаем внешнюю

нагрузку

в ряд:

 

 

 

</) sin

и

sin S Sи L -

(3.70)

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

,,,

4РД

. mnvxt

. nnvut

Н,w = ~ ± ~ sm — г ~ 51П — f -

где vx, vy — составляющие

скорости

перемещения подвижной

массы. Подставив разложения (3.69) и (3.70) в уравнение колеба­

ний и положив г\=

vxt, получим уравнение относительно обобщен­

ных координат ц>тп:

 

 

 

тп (Т|)

 

 

йц2

' “Ь ^тпфт/г (л)

 

Р (Т)) я2

sm тя?|

/гяг)|л2

а2а2 (1 + Pi)2

 

где

 

 

 

 

л т 2 + " 2Р?

Р(У]) = ^ Ж

 

1 +

pf

Л

а =

vxa У р h

И2 =

----- 7=

 

 

лУ О (1 + р2)

 

134

Интеграл уравнения имеет вид

 

 

фтп ('П) “ А-тп ('По) ^OS &тпТ\“Ь &тп (Но)

&тпЦ “f"

 

Г

тлк

птйлхуУх?

sin km„ (n — A,) dk

 

л \ Р (к) sin-------- sin--------------

 

J

a

a

 

 

_ i _

________________________________________________________________________________________________________

 

 

aa (m2 + n2ц^) (1 + Hi)

 

Если груз

движется

параллельно оси

х, то vy — 0 и выражение

для обобщенной координаты принимает вид

 

 

4>тп (л) — Атп (Ло) cos kmnr\+

Втп (л0) sin kmnx\+

 

л sin

 

Ч

 

 

+ ■aa (m2+ «Vi) i1 + Hi)

j" Р (X) sin - тлХ sin kmn (л — A,) dX,

rj„

 

 

где уг — координата

линии

действия

нагрузки.

Вынося среднее

значение коэффициента динамичности давлений Рср за знак инте­ грала, получаем

 

 

/ \

 

 

/

\

 

I

1

 

^<Ртп

Smn +

 

фт п ( л ) — фт п (Ло) dm n

 

 

 

dTWi=r

 

 

плуф!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P c p sin

 

a

 

 

 

 

тлri

. тлrin

 

+ ■(,m2 +

и2(лj)2 — /n2a2 (1 +

H i)2

Sin—^

-----Sin----—

 

 

 

 

 

 

 

 

 

 

 

 

(1 + m) та

^

 

 

тЩо

 

 

 

 

 

 

 

т2 + /г2ц2

Ътп

 

ЕЗо.1,

 

 

 

 

 

 

 

 

а

У

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Sm n

— sin k m n

(т]

Ло)’

 

Сгпп

COS k.n-n

Ло)*

Разобьем промежуток

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

на участки длиной т = ajr

и обозначим Рср = Р к+1

при kx <

л <

kx +

т. Положив

 

 

 

 

М (ГП,

п,

k) = фmn (kx),

 

 

 

 

 

 

N (т,

п,

k)

=

 

 

dVmn (Л)

 

 

 

 

 

 

 

d\=kx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

получим рекуррентные

формулы

 

 

 

 

 

 

 

 

М (т,

п, k +

l ) =

М (т,

п,

k) стп + N (т,

п,

k) smn - f

 

1 k+l

плу^

 

 

 

тл (k

1)

 

 

тлк

 

 

a

 

 

 

 

 

 

 

p* + iSin'

 

 

 

 

 

 

 

+

(т2+ n2Hi)3 — m2a2 (1 +

212

Sin

 

 

 

 

" Cmn Sin ■

 

h,)1

 

 

 

 

 

 

 

 

 

 

 

 

(1 4- Hi) ma

®яи COS

тлк

 

 

 

 

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

m2 4- a2p,

 

 

 

 

 

 

 

 

 

N (m,

n, k 4-

1) =

N (m, n,

k) cmn — M (m,

n,

k) smn - f

135

+

где

pk+\sin

nnyxV-x

 

' Щ (1

-f |Ij)

 

гая (ft + 1)

 

a

 

COS

(ma + п2Ц|)2 — m2a2 (1 +

M-i)2

m2 +

n2

 

 

 

ma (1 +

Ц])

COS

mrtft

+ Smn Sin - mnk

7

9

2

m3 + я2^|

 

 

 

 

 

я (m2 + n2|X[)

.... = cos-

Я (/Л2 -f- Л2|1[)

= sm-

(l +

p,2)ar

 

(l+(x|)ar

Вертикальная составляющая траектории груза определяется из уравнения

Ро

d?z„

=

P 0 — PR(t),

g

dt2

 

 

где Р д — динамическая

реакция

балки, равная силе давления

движущегося груза. Приняв гд.г =

w0zr, найдем

гг (11) =

гг (т)0) + J

Zr

(Л — По) +

" (Т1

_ Т1о)

- ( 1 - ^

с

 

dTlri=rb

Vl Ч0/

' 8a2p(l

+ M.?)2a2

 

где § =

P^JcP hpg.

При

т]0 =

kx,

т] =

kx + т,

Ргр =

Pk+i,

т) — % =

т = у получим

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

соотношения, аналогич­

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

 

Q ( k + l) =

Q(k) +

R(k) +

8a2f5r3 (1 +

м|)3

(1 — Pk+l),

 

 

 

 

 

 

 

 

 

 

 

 

 

R ' k + l ) ^ R ( k )

4 a 2pr2 (1 +

fxf)2

(1 Pk+1),

 

 

 

 

 

 

 

 

 

 

 

 

где

Q (k) =

zr (kx),

R (k) = т

dz,

 

Условие

пересечения

dr\

ti=*t

 

 

 

 

 

траектории

груза

с

поверхностью

пластины

принимает

вид

 

 

 

2Г (л) =

2

2 фт п

(Л) sin -253L sin

b

 

(3.71)

 

 

 

 

 

 

 

 

а

 

 

 

или в приведенных выше обозначениях для г] =

kx +

т

 

 

Q (k +

1) = 2

2

М {tn, п, k +

1) sin

OTTI(fe+

Ч . sin

ПЩхУ-х .

 

 

 

т

п

 

 

 

 

 

г

 

 

а

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.71а)

Из

этого

условия

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

Рт+ ь

обобщенные

координаты

Фотп (л) и прогибы w.

 

 

 

 

 

 

 

 

 

Вертикальную скорость груза при х\ = kx обозначим

vz. Тогда

 

 

 

 

 

a 2 {kx) =

vz (kx) а У ph

 

 

 

 

 

 

 

 

 

я Y d (i +

Hi)2

 

 

 

 

 

 

 

 

 

 

 

 

 

136

Значение а 2 может быть найдено в процессе вычислений:

a z (kx) =

4а3{5/-(1 Н Нч)

R № ,

 

 

 

я2уд

 

 

 

Уд =

v2/ag.

 

 

При косом ударе вертикальная

v2 и горизонтальная v

составляю­

щие скорости соударения заданы. По v2 определяется

R (0). Да­

лее расчет ведется по приведенной выше

схеме.

 

В качестве примера расчета рассмотрим пластину (см. рис. 37) размерами а = 200 см, b = 100 см, h = 2 см. Примем а = 0,75, Р = 0,25. Груз двигается параллельно оси х вдоль оси симметрии. В разложениях (3.69) удерживаем 30 членов ряда по т и 11 членов ряда по п, соответствующих нечетным гармоникам. Интервал ин­ тегрирования разобьем на 200 участков.

При расчетах используем условие (3.71), дополненное учетом,

контактных деформаций и неровностей аналогично (2.94):

 

Q ( k + 1) - I > ' Z M ( m ,

п, k-\- 1)sin

mJt(* + 1 )

s i n ^ ^

-----

m n

 

 

 

 

zH(kx + t ) kxP k+lq + Zt

= 0,

213,

 

причем zx = kx = kP 4T/w0,

w0 = 0,176

см и примем kx =

0,00242.

Рассмотрим пластину с синусоидальной вогнутой неровностью1 глубиной 0,2 см и протяженностью 20 см и сопоставим несколько

случаев: 1) неровность лежит в

интервале 0,2 < ; г\/1 <

0,3; 2) не­

ровность находится в интервале

0,45 <

г\/1 < ; 0,55;

3)

неровность

лежит в интервале 0,7 < r\jl <

0,8; 4)

пластина без

неровностей.

Обозначим шСт = 0,0496 w0 наибольший

статический

прогиб пла­

стины под действием массы, покоящейся посередине. Сопостав­ ление динамических коэффициентов Wi = w/wcr в сечении х = а / 2

для случая 1 (wx= 1,40), случая 2 ( щ = 1,38) и случая 3

(wx—

= 1,31) с динамическим коэффициентом для случая 4 (wx =

1,30)

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

(при а = 0,25 Ш1= 1,55, при а = 0,5

a>i = l,52). При

а = 0,75

в

случае 1 груз совершает три отскока,

в случаях 2 и 3 — два и

в

случае 4 — один отскок. Первое нарушение контакта

между

грузом и пластиной происходит при входе на неровность, а после каждого соударения наблюдается скачок динамического давле­ ния = 4 -у 8). Перед сходом с пластины груз, соударяясь с пластиной (г)/а = 0,9), вызывает резкий всплеск динамического давления. Напряжения в пластине можно определить, если счи­ тать, что контактная сила распределена по малой площадке [50].

Г л а в а ч е т в е р т а я

НЕСТАЦИОНАРНЫЕ КОЛЕБАНИЯ РОТОРОВ ТУРБОАГРЕГАТОВ

ПРИ ВНЕЗАПНОМ КОРОТКОМ ЗАМЫКАНИИ ГЕНЕРАТОРА

Колебания системы вал— диск — лопатки

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

к. з.— один из наиболее

опасных

с точки зрения воздействия

на систему роторов. При

к. з. на

бочку ротора генератора начи­

нает действовать внезапно приложенный переменный во времени электродинамический крутящий момент [4,6], амплитудное значе­ ние которого в несколько раз превышает номинальное. Этот мо­ мент является источником возникновения интенсивных крутиль­ ных колебаний всей системы роторов, которые с течением времени быстро затухают вследствие внешнего и внутреннего сопротивлений в механической системе и вследствие затуханий токов к. з. Несмот­ ря на кратковременность режима в валопроводе и рабочих лопат­ ках роторов могут развиваться значительные внутренние усилия.

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

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

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

138

нужденные колебания определены в виде разложения по собствен­ ным формам системы Бал — лопатки, найденным методом динами­ ческих жесткостей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

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