УДК 534.832.524:621.671
ТЕОРЕТИЧЕСКОЕ ИССЛЕДОВАНИЕ СИЛ ГИДРОДИНАМИЧЕСКОГО ПРОИСХОЖДЕНИЯ, ВЫЗЫВАЮЩИХ ВИБРАЦИЮ ЦЕНТРОБЕЖНЫХ НАСОСОВ
канд. техн. наук В.Е. ПИТОЛИН (Полоцкий государственный университет)
Изложены результаты теоретических исследований сил гидродинамического происхождения, действующих на рабочее колесо центробежного насоса, вызывающих вибрацию его корпуса и излучение звука на лопастной частоте и гармониках. Разработана методика расчета этих сил на различных режимах работы насоса, позволяющая производить оценку амплитуд и фаз вибрации и звука на лопастной частоте по его конструктивным и режимным параметрам уже на стадии проектирования. Результаты расчетов по методике, выполненные для восьми типов консольных центробежных насосов, показали хорошее совпадение с результатами натурных испытаний.
При выполнении расчетов гидродинамических сил, вызывающих вибрацию центробежного насоса (ЦН), а также являющихся основным источником звука на лопастной частоте, необходимо рассмотреть и рассчитать поля давлений и скоростей перекачиваемой среды в спиральном отводе в их взаимодействии с ядрами потоков в выходных сечениях межлопастных каналах рабочего колеса ЦН. До настоящего времени вопрос теоретического анализа этих сил остается открытым. Результаты исследований, выполненных в последнее время [1], несмотря на их основательность и полноту, не позволяют разработать методики, обеспечивающие расчет ожидаемых амплитуд и фаз вибрации ЦН на стадии проектирования.
Постановка задачи исследования. В предыдущих статьях автора была представлена кинематическая схема расчета амплитуд и фаз вибрации корпуса ЦН под действием динамических составляющих поперечных сил Ек, приложенных к рабочему колесу, на различных режимах его эксплуатации [2], а также изложен порядок расчета этих сил и определен метод расчета амплитуд и фаз звукового давления, излучаемого в гидравлический тракт на лопастной частоте [3]. Но для расчета динамических составляющих поперечных сил Ек, как это следует из системы уравнений (36), требуются аналитические выражения для расчета статического давления по периметру рабочего колеса Р2(ф) и радиальной составляющей скорости ядра потока на выходе из рабочего колеса при его движении относительно спирального отвода Ую(ф, /).
Аналитические методы расчета характеристик полей статического давления и радиальных составляющих скоростей потока на выходе из рабочего колеса можно получить решая совместно уравнения, составленные для описания структуры течений в спиральном отводе ЦН при различных подачах и потока рабочей среды в выходном сечении межлопастного канала рабочего колеса.
Исследование структуры течений в спиральном отводе ЦН на различных режимах его работы. Спиральный отвод является важнейшим конструктивным элементом ЦН. Его геометрические параметры существенно влияют на режим работы насоса.
В настоящее время наиболее широко распространены спиральные отводы, построенные с соблюдением закона постоянства момента скорости (рис. 1):
(1)
где Уи - окружная составляющая скорости потока в отводе ЦН на произвольном радиусе г. Поэтому площади поперечных сечений спирального отвода выбираются таким образом, чтобы обеспечить осевую симметрию потока за рабочим колесом на расчетном режиме работы ЦН [4]:
Ор
Г/2(Ф)= „р =соп«, (2)
2-П-Щ -£>2
где индекс «р» обозначает расчетный режим или расчетный геометрический параметр; Я2 - наружный радиус рабочего колеса; О -расчетный расход (подача) рабочей среды через насос; Ь2 - ширина межлопастных каналов в районе выхода из рабочего колеса.
На нерасчетных режимах отдельный межлопастной канал рабочего колеса находится под действием постоянного статического давления на входе, а статическое давление на выходе из канала периодически изменяется в зависимости от углового расположения канала относительно спирального отвода. Такое изменение давления, наложенное спиральным отводом, работающим в режиме конфузора (О > О) или диффузора (О < ОР), а также гидродинамическим следом обтекания языка спирального отвода (О ~ ОР) приводит к периодическому изменению относительной скорости ядра потока в каналах рабочего колеса.
Это и является причиной периодического изменения величины и направления суммарной радиальной силы, действующей на рабочее колесо со стороны потока перекачиваемой среды [3].
Рис. 1. Расчетная схема спирального отвода ЦН: Кс - радиус стенки спирального отвода в рассчитываемом сечении; К - радиус средней линии спирального отвода в рассчитываемом сечении; К3 - радиус кольцевого зазора (расстоянием от оси вращения до языка спирального отвода);
К2 - наружный радиус рабочего колеса насоса; К - произвольное значение радиуса расчетных параметров потока в спиральном отводе; Ь3 - ширина (или эквивалентная ширина) спирального отвода; Ф - координата расчетного сечения спирального отвода 3 - 4; Б - исследуемые сечения выделенного элементарного объема 1 - 2 - 3 - 4 спирального отвода
Из уравнения Эйлера для лопастных машин при условии осевой симметрии потока на входе можно определить окружную составляющую скорости потока в спиральном отводе в районе выхода из рабочего колеса [4]:
Г/2(Ф) = ii■ (Щ ■ сор -З&М) = const, (3)
tg$2
где д = 0,74 - 0,84 - коэффициент, учитывающий влияние конечного числа лопастей (порядок расчета этого коэффициента достаточно подробно изложен в работах [4 и 9]); ю - угловая скорость вращения рабочего колеса ЦН; VR2 - радиальная составляющая скорости потока на выходе из рабочего колеса; Рг - угол установки лопастей рабочего колеса в выходном сечении.
Местное значение окружной составляющей средней скорости потока в отводе на расчетном режиме может быть найдено из уравнения постоянства момента скорости (1):
Rp
• (4)
г
Расход через произвольное сечение спирального отвода с учетом выражения (4) может быть определен следующим образом:
Rc Op) Rc OP)
Qp(ф) = \ Vur(<p)-b3(r)-dr = ^(ф)-R2?- J-3^-dr. (5)
Rp Rp
Для спирального отвода прямоугольного сечения при b3(r) = const после интегрирования выражения (5) получим:
Qp (<p) = VU2 (Ф) * R2PA-in Здесь А - постоянная интегрирования.
A RM . Rp .
Согласно выражению (2)
Qp (Ф) =
271
следовательно,
Rp
Rc (ф) = ^р-ехр
Qp-<P
г■ ь
где А - постоянная интегрирования; Г - постоянная спирального отвода (Г = 2 • к ■ ¥/2 ■ Щ). Постоянная интегрирования А определяется из граничных условий при ф = 0 - Яс = Я3:
, R2p
A = —2-
1
(7)
(8)
(9)
где § = - коэффициент кольцевого зазора.
Тогда для прямоугольного в поперечном сечении спирального отвода (см. рис. 1) текущая площадь сечения 5(ф) будет определяться по формуле
s (ф)=v rc (ф)-r2i v r2
8exp
-1
а расстояние до средней линии спирального отвода
(10)
(11)
Если форма поперечного сечения спирального отвода существенно отличается от прямоугольной, то с достаточной для прикладных вычислений точностью можно воспользоваться эквивалентным значением Яс(ф), полученным из сравнения площадей:
rc (ф)= r2
1+
S (Ф) R2 b
(12)
С учетом выражений (8) и (9) можно получить выражение для радиуса стенки спирального отвода Лс(ф), входящее в формулу (11)
Rc (Ф) = 5 • R2 ' exp
кг-b
= 5-r2- b<p ,
(13)
(
где В = ехр
Qp
уг-ь3)
= const- расчетный параметр спирального отвода прямоугольного сечения.
Тогда выражение для Я(ф) (11) можно привести к виду
д(ф) = 1.д2.(5.5Ф+1Х
а выражение (10) для вычисления площади поперечного сечения спирального отвода к виду
(14)
(15)
Как уже отмечалось, течение в спиральном отводе имеет сложную пространственную структуру. На его характеристики оказывает определенное влияние пристеночный пограничный слой, а также величина расхода рабочей среды через насос, отличающаяся от расчетной.
Для расчета структуры потока в спиральном отводе ЦН на нерасчетных режимах его работы воспользуемся методикой, изложенной в работах Б.И. Боровского [5; 6].
Методика построена на результатах исследования течений в спиральном отводе на основе соотношений спутных турбулентных потоков. Поток на входе в спиральный отвод и поток на выходе из рабочего колеса имеют наибольшие углы относительно нормали к радиусу, схематично эти потоки рассматриваются в каждом последующем контрольном сечении как две спутные плоскопараллельные турбулентные струи (рис. 2).
Рис. 2. Расчетная схема течений в сечениях спирального отвода ЦН: Vuc - окружная составляющая скорости потока в районе наружной стенки спирального отвода; VuR - окружная составляющая скорости потока на средней линии спирального отвода; Vur - окружная составляющая скорости потока в произвольном месте сечения; Vu2 - окружная составляющая скорости потока на выходе из рабочего колеса
Полагая, что один из потоков имеет скорость Км2(ф) = const, а другой - Ус(ф), при которой в силу кинематических соотношений, определенных для сечения турбулентного пограничного слоя, в середине рассматриваемого сечения отвода реализуется скорость ¥иг(ф). Эта скорость определяется из уравнений сохранения энергии, момента количества движения и уравнения неразрывности.
Автором этой схемы расчета Боровским Б.И. [6] получено и исследовано для отводов различной формы, в том числе и для лопаточного отводящего аппарата, соотношение для окружной скорости потока в отводе в рассматриваемом контрольном сечении 3 - 4 (см. рис. 1) в функции от радиуса отвода:
г;дф)
к
u 2
к
u2
1-0,385-
1 —
r-R,
2-(ЛОрЬЛ,)
(16)
Усредненное по периметру рабочего колеса значение окружной составляющей скорости потока на выходе из колеса Уи2 рассчитывается по формуле, аналогичной формуле (3), но с рабочими значениями расхода Q и внешнего радиуса колеса Я2, которые могут отличаться от расчетных:
V„2=V
Rf • со -
Q
ln-b2-R2-tg$2jj
(17)
Выражение для расчета статического давления в контрольном сечении 3 - 4 на входе в спиральный отвод со стороны рабочего колеса получено из условия радиального равновесия
(
Pr (ф) = Pr (ф)+р -к!
0,25 + 0,75-
vur (ф)
\r(<p)J'
(18)
а выражение для расчета статического давления на средней линии отвода, входящее в это уравнение, полученное из уравнения Бернулли, имеет вид
рй(ф) = р1+р-\ïïk-hom-kàm-ka
r( ф)
V Р2
\
-1
0,25—^г
У
u2
Л
UK V т / _
У
u2
-0,03
(19)
Значение скорости потока на средней линии спирального отвода Кмй(ф), входящее в уравнения (16), (18) и (19), находится из краевых условий при ф = 0 и ф = 2п
[;л(Ф) = р;л(2тг)
R( Ф)
R
-1
R(2n)
R
-1
пКл(2 л) Vu2
R(0) ,
(20)
Эти выражения не содержат зависимости R(ф) в явной форме, так как предназначены для численных расчетов параметров течений в отводах с поперечными сечениями произвольной формы.
2
R
2
Замена формы сечения спирального отвода на эквивалентную прямоугольную вносит определенную погрешность в расчеты структуры потока (по некоторым оценкам, до 10 %). Однако погрешность измерительных трактов, используемых для оценки абсолютных уровней виброускорений корпуса насоса, обычно составляют 25 - 30 % [7], поэтому принятые упрощения можно считать допустимыми.
После подстановки выражения (14) в уравнение (19) получим формулу для расчета окружной составляющей скорости движения жидкости в спиральном отводе на режимах, отличающихся от расчетного Уиг(ф), уже в явной форме от координаты ф
ГиК(ц>) = ГиК(2л)
' 5 . _ 1 ,
У,2
\ £ о2тг
Ь-В2п-1
5-1
(20)
Окружную составляющую скорости потока в «горле» спирального отвода Уик(2%), входящую в уравнение (20), можно оценить из уравнения расхода
тогда
Кя(2л) =
0(271) = 0 + 0(0), б(2п) _б + УиК(0)^(0)
5(271)
5(271)
(21)
(22)
После подстановки в эту формулу выражений для площадей поперечного сечения спирального отвода (15) можно получить
Ь3-Я2-(8-В2л-1)
(23)
Начальное значение скорости VuR(0), входящее в формулу (23), определяется краевыми условиями и может быть оценено, как
К*(0)*>К2-
(24)
Формула для расчета статического давления на средней линии спирального отвода может быть получена аналогичным образом из выражения (19):
_ 1 1 _
Рк (Ф) = Р„ + р • • а> • Ки2 - Иот - - ■ Ги2к (Ф) - - • Ки22 ■ (8 • В* -1)
025Лд(ф) 1 ^я(ф)
V -
Л
V.
-1
-0,03
(25)
Потери в спиральном отводе йОЯ! можно ориентировочно оценить по формуле:
^от ~ 4' 2 '
(26)
где 4 - коэффициент гидравлических потерь.
Согласно известным опытным данным [4] коэффициент потерь в спиральном отводе и коническом диффузоре 4 зависит от режима работы отвода: на расходе меньше расчетного он определяется потерями в отводе, а при расходе больше расчетного - потерями в коническом диффузоре.
В диапазоне расходов (0,5 - 1,6)^р коэффициент потерь может быть оценен, как 4 ~ 0,2.
Усредненная во времени величина радиальной составляющей скорости на выходе из рабочего колеса может быть найдена из уравнения неразрывности
Г*2(Ф) =
1
А.
Клад-ад".
(27)
После подстановки в это уравнение выражений для скорости ^(ф) (20), площади поперечного сечения отвода 5(ф) (10) и последующего дифференцирования можно получить формулу для вычисления усредненных во времени радиальных составляющих скорости среды в спиральном отводе на выходе из рабочего колеса:
1п
V,,
1п
5-1
- + 1
2%
Уик( 2л)
[1-
8-1
(28)
и 2
Статическое давление в районе выхода из колеса Р2(ф) можно рассчитать по формуле (18) при г = К2 и подстановке в эту формулу выражения для Я(ф) (14):
(
р2(ф) = рк (ф) + р "к!"
КК(Ф)
0,25 + 0,75—¡К
К
и 2
■1п
2
+1
(29)
Как показывает практика, расчет параметров течения в спиральном отводе в районе выхода из рабочего колеса, выполняемый по формулам (20), (25), (28) и (29), позволяет получить хорошее совпадение с экспериментальными данными, однако нуждается в дальнейшем совершенствовании для анализа структуры течений в районе «языка» спирального отвода и влиянии кольцевого зазора на режимах, близких к расчетным.
Рассмотрение реальной структуры течений на выходе из рабочего колеса ЦН. Известно, что поток рабочей среды в выходном сечении межлопастного канала рабочего колеса не распределяется равномерно по периметру, а обладает хорошо выраженным ядром, расположенным в районе рабочей (напорной) стороны лопасти [8], так называемой области вытеснения.
В работах академика Г.Ф. Проскуры [9] достаточно наглядно показано влияние обратного вихря, развивающегося в межлопастном канале рабочего колеса ЦН, приводящего к перераспределению среды в канале вплоть до отрыва потока от тыльной стороны лопасти на малых расходах и образованию ядра в области вытеснения. Согласно утверждению ученого, реальная картина течений в межлопастных каналах рабочего колеса «... может сопровождаться разрывом сплошности и не отвечает предположениям, положенным в основание гидродинамической теории» [9, с. 318]. По этой причине для расчета реальной структуры течений на выходе из рабочего колеса используются различные эмпирические зависимости. Например, можно воспользоваться моделью структуры потока в выходном сечении канала от координаты активной зоны, полученной А.И. Тимшиным [8], которая имеет следующий вид
прие<\|/<у - кй2(у) = 0,
г
(30)
где ф - координата точки активной зоны потока; е - длина по шагу активной зоны потока; дг= - расход рабочей среды через межлопастной канал (р - число лопастей); у - угловое значение шага установки лопастей. Длина активной зоны потока е существенно зависит от расхода через насос
е = 2 ■ Я, ■ а
8-и2 -Я2 у
(31)
где а - коэффициент пропорциональности, определяемый опытным путем для конкретной конструкции рабочего колеса (обычно находится в пределах 1,8 - 2,6); и2 = шК2 - окружная составляющая скорости наружной кромки рабочего колеса.
Отсчет координаты ф осуществляется от точки пересечения активной (лицевой) стороны лопасти с окружностью к2 по направлению вращения рабочего колеса. Поэтому возможны некоторые трудности при переходе от системы отсчета угла ф, жестко связанного с рабочим колесом (рис. 3), к системе отсчета угла ф, жестко связанной с корпусом насоса, относительно которой система отсчета ф вращается вокруг оси 02 с постоянной угловой скоростью со.
\ ■
Рис. 3. Расчетная схема эпюры радиальных составляющих скоростей потока на выходе из межлопастных каналов рабочего колеса ЦН
2
0.44
Для упрощения процесса перехода от одной системы отсчета к другой относительно некоторой точки на корпусе насоса с координатами (Я2, ф) можно записать
ф=у-ф = у-ю-?, (32)
полагая, что начало отсчета времени t совпадает с моментом прохождения относительно этой точки выходной кромки рабочей стороны лопасти рабочего колеса.
Исследование взаимодействия течений на выходе из рабочего колеса и в спиральном отводе ЦН.
Согласно формуле (31) величина усреднена по каналам рабочего колеса, а величина ^(ф) имеет осевую симметрию, что не соответствует реальной картине течений.
Для обеспечения закона неразрывности потока среды через насос можно записать для канала
Чг (ф)= V К2
КЛ2(ю t )■ Л .
(33)
Полагая значение постоянной интегрирования равным нулю, после интегрирования выражения (33)
(34)
Из выражений (28), (30) и (34) можно получить недостающую формулу для расчета мгновенных значений радиальных составляющих скорости потока на выходе из рабочего колеса
2-<?г(ф)- 1(ф)~ф_
е(ф) -Ь2
(35)
где у = ф-со-/ —
п-2-п
е(ф) = 2-К2 -а-
8-ю-Л
0,44
Отсчет угла юt осуществляется относительно начала отсчета координаты ф (см. рис. 3).
Граничные условия для формулы (35):
- в диапазоне 0 < ф < е(ф) значения УЯ2(ф, Г) вычисляются по формуле (35);
- в диапазоне е(ф) < ф < у значение УЕ2(ф, 0 = 0;
- для определения значений ф в интервале 0 < ф < у = 2п/г множитель п последовательно принимает значения ряда п = -1, 0, 1, ..., г - 1;
- координата ф изменяется в диапазоне от 0 до 2п.
Процесс расчета параметров Ую(ф, () по формуле (35) осуществляется следующим образом:
1) при заданном значении определяющем мгновенное положение рабочего колеса ЦН относительно языка спирального отвода, рассчитываются значения ^2(ф, t0) с учетом закона неразрывности среды (33) и граничных условий, действующих по периметру рабочего колеса;
2) затем колесо поворачивается на угол dф = юdt, и расчет значений КД2(ф, t0 + Л) по формуле (35) повторяется;
3) при достижении колесом положения, соответствующего углу поворота ф = ю^ + у, процесс вычислений прекращается, так как результаты последующих вычислений будут повторять результаты предыдущих с периодом Т = у/ю.
Вычисления значений статического давления Р2(ф) по формуле (29) достаточно выполнить один раз вне зависимости от положения рабочего колеса.
Результаты вычислений вводятся в систему уравнений для расчета значений поперечных сил, действующих на рабочее колесо ЦН со стороны жидкости, вывод которых выполнен в работе [3]:
2тг ^
Рцх (0 = ~Ь2К2 1 Р2 (ф) с0!3 ф^ф--г"1
о 2
Ркг (/) = -Ь2К21 Р2 (ф) втфс/ф -
1 ^(ф^совфЛр- | КЛ2(ф,/')Ки2(ф)8тфй?ф
1 ^(ф^^Пф^ф- i ^2(ф'Ж,2(ф)С08фЛр
(36)
Затем выполняется численное преобразование Фурье последовательности полученных значений, в результате которого вычисляются амплитуды и фазы динамических составляющих поперечных сил, возбуждающих вибрацию корпуса и генерирующих звуковое давление в гидравлический тракт ЦН на лопастной частоте.
0
3
г
о
о
о
о
о
Формулы для расчета уровней вибрации и звука приведены в работах автора [2; 3].
Результаты вычислений вибрации корпуса, выполненных для 8-ми типов ЦН консольной конструкции различного назначения (ЦНР-25/35, ЦН-100/30, НЦКВ-100/35, ЦН-104, ЦН-160/80, БЭН-38, 10ВФ-18, 10ВФ-24), достаточно хорошо совпадают с результатами натурных испытаний. Более того, расчетами подтверждаются и результаты изменения конструктивных и режимных факторов: уменьшение или увеличение частоты вращения рабочего колеса ЦН, уменьшение наружного диаметра рабочего колеса, изменение формы языка спирального отвода и т.д.
Заключение. Разработаны аналитические методы расчета гидродинамических параметров структуры течений в центробежном насосе на различных подачах и скоростях вращения рабочего колеса.
Создана методика расчета амплитуд и фаз динамических составляющих поперечных сил FR, действующих на рабочее колесо ЦН, на различных режимах его эксплуатации, в которой изложен порядок расчета этих сил и определен метод расчета амплитуд и фаз вибрации корпуса ЦН на лопастной частоте. Разработанная методика может использоваться для оценки конструктивных параметров ЦН, направленных на снижение вибрации на этапе проектирования насосов. Результаты расчетов амплитуд и фаз вибрации корпуса ЦН на лопастной частоте подтверждаются результатами натурных испытаний 8-ми типов насосов.
ЛИТЕРАТУРА
1. Перевощиков, С.И. Разработка научных основ управления вибрацией гидродинамического происхождения в центробежных насосах магистральных нефтепроводов: автореф. ... д-ра техн. наук / С.И. Перевощиков; ТПИ. - Тюмень, 2004. - 47 с.
2. Питолин, В.Е. Построение компьютерной модели центробежного насоса для расчета вибрации его корпуса под действием динамических составляющих поперечной силы / В.Е. Питолин // Вестн. Полоц. гос. ун-та. Серия F. Строительство. Прикладные науки. - 2009. - № 12. - С. 95 - 101.
3. Питолин, В.Е. Построение математической модели процесса возбуждения динамических составляющих поперечных сил в рабочем колесе центробежного насоса силы / В.Е. Питолин // Вестн. Полоц. гос. ун-та. Серия C. Фундаментальные науки. - 2011. - № 4. - С. 52 - 58.
4. Ломакин, А.А. Центробежные и осевые насосы / А.А. Ломакин. - М., Л.: Машиностроение, 1966. -264 с.
5. Боровский, Б.И. Особенности течения в спиральном отводе центробежного насоса на нерасчетных режимах / Б.И. Боровский // Изв. ВПИ. Гидродинамика лопаточных машин и общая механика. - Воронеж, 1977. - С. 42 - 49.
6. Боровский, Б.И. Особенности расчета поперечных (радиальных) сил в центробежных насосах с направляющим аппаратом и спиральным отводом / Б.И. Боровский // Изв. ВПИ. Гидродинамика лопаточных машин и общая механика. - Воронеж, 1977. - С. 11 - 17.
7. Клюкин, И.И. Акустические измерения в судостроении / И.И. Клюкин, А.Е. Колесников. - Л.: Судостроение, 1982 - 256 с.
8. Тимшин, А.И. Структура потока на выходе из колеса и ее влияние на характеристики центробежного насоса: автореф. ... канд. техн. наук / А.И. Тимшин; ХГУ. - Харьков, 1972. - 32 с.
9. Проскура, Г.Ф. Гидродинамика турбомашин / Г.Ф. Проскура. - Киев: Машгиз, 1954. - 423 с.
Поступила 22.11.2012
THEORETICAL EXAMINATION OF FORCES OF THE HYDRODYNAMIC ORIGIN, SETTING UP VIBRATIONS OF CENTRIFUGAL PUMPS
V. PITOLIN
Results of theoretical examinations offorces of the hydrodynamic origin acting on the vane wheel rotor of the centrifugal pump, setting up vibrations of its housing and radiation of sound at vane frequencies and harmonics. Method of calculation of these forces at various modes of pump operation is developed, allowing to estimate amplitudes and phases of vibration and sound at vane frequencies on its constructive and mode parameters on the design stage already. Results of calculations made according to the method for 8 types of cantile-vered centrifugal pumps have displayed a good match with the results of full-scale tests.