УДК 621.65
В. Г. МОХНАТКИН, П. Н. СОЛОНЩИКОВ, А. С. ФИЛИНКОВ
ТЕОРЕТИЧЕСКОЕ ОПРЕДЕЛЕНИЕ ГИДРАВЛИЧЕСКИХ ХАРАКТЕРИСТИК ЛОПАСТНОГО КОЛЕСА УСТАНОВКИ ДЛЯ ПРИГОТОВЛЕНИЯ ЖИДКИХ КОРМОВЫХ СМЕСЕЙ
Ключевые слова: напор, подача, скорость, сопротивление, установка, характеристика.
Аннотация. В данной статье представлено теоретическое определение напорно-расходной характеристики в установке для приготовления жидких кормовых смесей. На основании проведенных исследований получены уравнения, позволяющие определить напор и расход при различной частоте вращения.
Линии тока относительного движения и представляют собой произвольные прямые 1/-2/ и 1//-2// (рисунок 1). В нашем случае при данной расчётной схеме удельная энергия И" > И[ на величину
удельной энергии И', которую среда приобретает в рабочем колесе за счет воздействия на нее лопастей и которая определяется по уравнению:
Н' = ^2^pL+с^ + z; —^ (1)
Pg 2g
где ———— удельная энергия давления в точках 2/ и 1/;
Рё
- удельная энергия положения в точках 2/ и 1/;
_2 2'
с — с ъ/ 1 /
— ;-удельная кинетическая энергия в точках 2 и 1.
2 я
Для второй ступени И" > И" поэтому удельная энергия, приобретаемая в ней средой, будет равна
© Мохнаткин В. Г., Солонщиков П. Н., Филинков А. С., 2014
79
Н'' = р2 р 1 + с2 с12 + ^ — 2", (2)
Рё
где ———— удельная энергия давления в точках 2// и 1//;
Рё
г.1 — - удельная энергия положения в точках 2// и 1//;
2 2
со — ^ ^1/„л1/
2Я
----удельная кинетическая энергия в точках 2 и 1 .
Рисунок 1 - Параллелограмм скоростей лопастного колеса установки
Определим составляющие уравнения (2.5 и 2.6). Используем для этого центр прямоугольной системы координат, который жестко свяжем с центром рабочего колеса (рис. 2.7).
Пользуясь системой дифференциальных уравнений движения жидкости (уравнения Эйлера), получим следующую систему:
^ ; др
Ф----^ = ;
р дх dt
ф — ; др = ^; (3)
у р дх dt
^ ; др duv
Ф----— = —- ,
2 р дх dt
где Фх, Фу, Фх - проекции единичных векторов действующих массовых сил на оси ох, оу,07;
их, иу, их - проекции относительных скоростей на оси ох, оу,
07;
р - гидростатическое давление в точке ; р - плотность жидкости.
Дифференциальные уравнения (3) преобразуется к виду интегрирования, умножением на производные йх; йу; йх
; (др др др — \ — dx +----dy +---<
Р дУ дг У (4)
du duz
= —- dx +---- dy +--- dz .
dt dt dt
Заменим в правой части полученного уравнения (4) величины
йх; йу; йх на выражения йх=ий; йу=иуйК й2=ихА:
; du duv du
ФСх + Ф С + ФСх-------dp =—-dx +------— dy +---сЬ . (5)
у р А А А
Сокращая уравнение (5) на А и подставляя в правую часть
уравнения сумму дифференциалов суммы ^ = ы1 , получим
ФхСх + ФуСу + ФС — ~ Ср = С \ ^. (6)
Проекции единичных векторов массовых сил в выбранной системе координат равны:
Фх = а)2х; Фг = (оО2 — ё); Ф2 = 0. (7)
фСх + фСу + ФС —| ~Сх + — Су + — Сх | =
Подставляя значения проекций единичных векторов массовых сил в уравнение (7), получим:
а хСх + а уСг — gdz------Ср = С \~^\. (8)
Разделив переменные и проинтегрировав, получим
р а2 х2 т2у2 и2
- = — + Ч~ — 82 — ~ + к, (9)
Р 2 2 2
где к - постоянная интегрирования.
Заменив х2 +у2=К2 и поделив на ускорение свободного падения g получим уравнение, которое позволяет определить величину гидростатического давления в любой точке канала:
р_ и2 .
— = —— 2—^~+к1, (10)
Р 2 2ё
где к]=^.
Используя полученные уравнения для определения гидростатического давления в любых точках, можно определить его величину в точках 1 и 2 , а также 1" и 2//. Таким образом, значение гидростатического давления в любых точках будет равно:
р1 = ^ — 21 —^ + К , (11)
Р 2 2ё
р2 = — 2- — «;[_ + , (12)
Р 2 2ё 1
р1 _ 0 ^ _ — 2" — ^ , (13)
Р2 = — 2» — ^ +^ . (14)
Р 2 2ё
Подставим вычисленные значения гидростатических давлений в уравнения (1) и (2), определяющие величину напора в каждой ступени:
Н = + к — ^ —
2 2 2я 1 2 1 2я
г 2' „2'
—к, + -
H’’ = - z - ^ + k - - + z - U-----
2 2g 1 2 1 2g
S S (16)
-k + —--— + z - z'.
2 g 2 b
После несложных преобразований получим значения удельной энергии для каждой ступени:
я, ш2 R - R )Ц - и\ | с2 - с2 2 g 2g 2g ’
;г ш (R - R’’ ),u22’’ - uf | —2’’ - —i2’’
2g 2g 2g '
Применяя данную теорему косинусов для определения u2 = о/ + с2 - 2 • т • с • cos а выражения (17) и (18), приводим к виду, известному как формула лопастных машин:
я, °2c2cos а2 - °1—1cos «1 g ’
н„ = m2—2cos «2 - т1 с1cos а1 (20)
g
Общий напор для двух ступеней будет равен: т’с2 cos а2 - т. —’ cos а[ + т2с2 cos а2 - т1 CY cos а’
H =---------------------------------------------------• (21)
g
За счет использования двухступенчатой схемы рабочего колеса поток приобретает повышенную удельную энергию.
Подача установки определяется из объема жидкости, поступающей из рабочего колеса за единицу времени. При этом ее значение будет определяься следующей зависимостью:
02 = 2•*• R • Ь • с;2, (22)
Q = 2 •*• R • К • с'г2, (23)
где Q'2 и Q2 - соответственно подача устройства на 1 и 2 ступени;
R и Я'2 - соответственно наружный радиус рабочего колеса 1 и 2 ступени;
К' и К'' - соответственно ширина колеса в выходном сечении 1 и 2 ступени;
с г2 и с” 2- соответственно радиальные (нормальная
выходному сечению) составляющие абсолютной скорости 1 и 2 ступени.
Объем жидкости, поступающей в единицу времени в рабочее колесо устройства ввода и смешивания, без учета стеснения потока лопастями на входе в колесо, можно определить по формуле:
е; = 2-я-я; • ь; • с;,, (24)
б; = 2-я-Я, • % • сп , (25)
где б и б - соответственно подача устройства на 1 и 2 ступени;
Я и Я - соответственно наружный радиус рабочего колеса 1 и 2 ступени;
Ь[ и Ъ\ - соответственно ширина колеса в выходном сечении 1 и 2 ступени;
с л и с” 1- соответственно радиальные (нормальная входному сечению) составляющие абсолютной скорости 1 и 2 ступени.
Вследствие неразрывности (сплошности) потока подача устройства Q на каждой ступени будет равна:
е = б, = 62 = е;= 62. (26)
Конструктивное исполнение рабочего колеса принимает следующий вид Ь = Ь = Ь2 = Ь = Ь2. При этом радиальные составляющие абсолютной скорости могут быть определены по формулам:
е ' е
2яЯ1Ь ; 2яЯ2Ъ
е ■ е
2яЯ Ь ’ ; 2яЯ
(27)
(28)
'2
Из параллелограмма скоростей (рисунок 1) тангенциальные составляющие абсолютной скорости могут быть определены как
сг; = ®; - С;Рё0; сг2 = Ю2 - С;2^8 02
(29)
<; = - СгР80; , сг2 = - С;2^802 . (30)
с;1 =
с;1 =
Подставим значения тангенциальных составляющих абсолютной скорости (29) и (30) в основное уравнение лопастного насоса (21) [1 с.5]:
H = ■
®2-I ®2 -q| +®" -U" -Q 1 +
(31)
ё
Приводим подобные значения с учетом того, что
V = ю[Я, У2 = ю2Я и V = ОЯ, V = ю2Я . При этом а\= т2 = со , величина напора определится выражением:
12)
н «2 -[(r22 -r"2)-(R22 -r;2)] +
® - [(ctgP; - CtgP2 ) + (ctgP; - CtgP2 )]
(32)
2nbg '
Или
H=A +BQ, (33)
где A W - [R - R;2) - R - R;2)^ V - [(Ctg0; - CtgP2) + (ctgPl - CtgP'j] g ’ 2nbg
Получены уравнения для определения напора установки при постоянной частоте вращения рабочего колеса п. Получены коэффициенты которые постоянны A = const и B = const, а напорная линия представляет собой прямую. Таким образом, при, Q=0 напор Н равен
я to1 -[(r22 -r;2)-(r22 -r;2)] (34)
g ’ а при напоре H=0
w2 -[(r22 - r;2)+(r22 - r;2)] - inb (ctg0 - ctg02 ) - (ctg Д" - Ctg02 ) '
Как видно из полученной формулы, вид зависимости И=/ (Q) будет главным образом зависеть от величины углов в2 и Pi. Поэтому линия, характеризующая напор, будет горизонтальной, если будет сохраняться условие р2=р1, то есть радиальное расположение лопаток.
Построим теоретические напорные характеристики для суммы напоров двух ступеней рабочего колеса (рис 2).
Как видно из рисунка 2, сумма напоров удваивает общий напор. Если имеем большое число лопаток, идеальная характеристика представляет прямую линию, причем характеристики, полученные экспериментальным путем, имеют криволинейный вид, близкий к параболе, это объясняется тем, что в основном уравнении Бернулли не учитываются сопротивления межлопаточных каналов, обусловленные трением о стенки и завихренностью жидкости [2 с.5].
Местные потери напора связаны с отставанием угловой частоты вращения потока от угловой частоты вращения лопастного колеса, вызванного конечным числом лопастей и неподвижных направляющих лопаток.
Рисунок 2 - Теоретическая напорная характеристика для первой сту-
і И = ДС н=Ж>)
пени и! = ЛО), второй ступени И = _ДО рабочего колеса и их сумме
Для учета сопротивлений в межлопастных каналах примем зависимость от расхода в виде, аналогичном для трубопровода:
К = К + +^202- (36)
Тогда зависимость Н=ДО) с учетом сопротивления запишется
как
Н = (Л - К) - - к&2, (37)
где к0 - удельные потери в каналах, м,
кі и к2 - экспериментальные коэффициенты, полученные в результате исследования характеристики установки (кі=В).
Ввиду того, что угол рі определяется геометрией профиля лопасти, он будет больше выходного угла р2, который можно определить
по значению величины B, полученной при экспериментальных исследованиях, решив уравнение (37).
Коэффициенты A, kj и k2 в уравнении (37) определяем, используя электронные таблицы «Microsoft Office Excel 2007» и метод наименьших квадратов, при этом получим:
Я750 = 2,31-0,012• Q-0,103-Q2, Я1000 = 3,11 -0,327• Q-0,009 • Q2, Я1500 = 6,97-0,292• Q-0,044• Q2.
(38)
(39)
(40)
По полученным уравнениям (38-40) строим напорную характеристику установки при частотах вращения рабочего колеса 750, 1000 и 1500 мин"1 (рис. 3).
Напорно-расходные характеристики (рис.3) получены по теоретическим расчетам по уравнениям (38...40). Видно, что для соответствующей частоты вращения напор будет меньше, так как они учитывали потери, принятые при экспериментальных исследованиях.
- теоретическая
- экспериментальная
Рисунок 3 - Напорные характеристики установки для приготовления жидких кормовых смесей
По составляющим уравнения (37) можно спрогнозировать напорную характеристику устройства при различной частоте вращения вала рабочего колеса с учетом соответствующих сопротивлений.
ЛИТЕРАТУРА
1. Русских В. М. Совершенствование конструктивно-технологической схемы молочного насоса с целью расширения его функциональных возможностей: автореф. дис. канд. техн. наук. Киров: 2002. 20 с.
2. Солонщиков П. Н., Совершенствование конструкции и оптимизация параметров установки для приготовления жидких кормовых смесей на базе лопастного насоса: автореф. дис. канд. техн. наук. Киров: 2013. 20 с.
DETERMINATION OF HYDRAULIC PERFORMANCE TEOERETICHESKOE IMPELLER USTANVOKI COOKING LIQUID FEED MIXTURES
Keywords: characteristic impedance, installation, pressure, speed, supply.
Annotation. This paper presents a theoretical definition of pressure-flow characteristics in a plant for making liquid feed mixtures. Based on the studies obtained equations to determine the pressure and flow rate at different speeds.
МОХНАТКИН ВИКТОР ГЕРМАНОВИЧ - д. т. н., профессор, ФГБОУ ВПО ВятГСХА, Россия, Киров, ([email protected]).
MOKHNATKIN VICTOR HERMANOVITCH - doctor of technical sciences, professor, Vyatka State Agricultural Academy, Russia, Kirov, ([email protected]).
СОЛОНЩИКОВ ПАВЕЛ НИКОЛАЕВИЧ - к. т. н., ФГБОУ ВПО ВятГСХА, Россия, Киров, ([email protected]).
SOLONSCHIKOV PAVEL NIKOLAEVICH - Ph.D., Vyatka State Agricultural Academy, Russia, Kirov, ([email protected]).
ФИЛИНКОВ АНДРЕЙ СЕРГЕЕВИЧ - к. т. н., ФГБОУ ВПО ВятГСХА, Россия, Киров, ([email protected]).
FILINKOV ANDREW SERGEEVICH - Ph.D., Vyatka State Agricultural Academy, Russia, Kirov, ([email protected]).