Файл: Филиппов, А. П. Воздействие динамических нагрузок на элементы конструкций.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' + |
2л л г ) ) |
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 ) 0ч |
—=z- |
||
|
sf + (и/ + Р)2 |
|
Здесь 2«/ = /г + xsyJ S/ |
= si ~ |
nh si ~- собственная круговая час- |
тота системы вал — диск — лопатки; щ — приведенный коэффициент демпфирования;
' П ( |
/„А2/+ 2я - М |
Г |
Rj = 2 |
й£ (р) р3 (ф;/)2 ф + |
|
Ь=1 |
ё |
г„ |
143