книги / Расчет центробежных и осевых насосов лопастных гидромашин
..pdfУгол поворота потока в решетке направляющею аппарата
А d |
оС4 ~ |
— 30 |
t |
(63) |
|
Скорость кевозмущеннсго |
потока |
|
|
|
|
V |
- - |
/ С - |
( К з / г ) г . |
|
(64) |
Угол невозиущенного |
потока |
|
|
|
|
|
|
|
|
|
(65) |
Число лопаток втравляющего аппарата выбираете;; по условию (46). Расчет треугольников скоростей в решетках профилей направляю щего аппарата также ведется в табличной форме для каждого расчет ного сеченияс
4.4. Пыбор густоты решеток профилей
J:;6op густоты периферийной решетки профилей имеет определяй
те: свойств цредвключешюго рабочего колеса и получения высокой экономичности. Особенности к рекомендации, касающиеся выбора гус тоты периферийной решетки профилей, подробно изложены е гл . 3* Дополнительно можно рекомендовать зависимость, связывающую густоту периферийной решетки профилей о требуемым критическим кавитацион ным коэффициентом быстроходности, полученную в работе [13] и пред ставленную на рис. 15.
Следует иметь в виду, что увеличение густоты решетки профилей приводит к возрастанию профильных потерь и снижает экономичность •насоса. Поэтому нужно стремиться к минимальной густоте на перифе рии, достаточной по антикавитационным требованиям.
Чрезмерное уменьшение густоты решеток профилей осевого насоса в условиях диффузорното характера течения жидкости в межлопастном канала молот привести к отрыву пограничного слоя, р е з к о м у возраста нию потерь ::ч вихреобразование и снижению напора но сравнению с расчетным. Поэтому при выборе закона изменения густоты по размаху лопасти обязательно .чолмна быть обеспечена оезотрпвнооть течв1ш я в мамлоиаотном канале на всех расчетных радиуса::.
Лисфузоркость ^ечеипя в можлопастном канале определяется трв-
41
О |
2 0 |
40 |
6 0 |
SO j i z f г р а д |
Р и с . 4b
42
мя основными параметрами: углом выхода потока ft2 , углом поворота
потока Aft и густотой решетки t/t . Следовательно, для конкретных значений ft2 к Aft ш л ю найти такую густоту, при которой в медлопастном канале наступит отрыв потока. Эти зависимости были полу чены экспериментально Хоуэллом [l, 8] с учетом 2С£-ного запаса по критического значению утла поворота Aft и представлены на рис. 16. По графику Хоуэлла определяются значения (£/£ )^оп в каждом рас
четном сечении.
Затем строится график изменения принятой густоты по размаху
лопасти. X первом приближении принимают его линейным так, чтобы густота на периферии соответствовала принятого ранее значению, а
в остальных /сечениях была бы не меньше допустимой по условиям безотршшости обтекания. В ходе пЬследующих приближений характер этой зависимости уточняется, но всегда она долина иметь монотонный ха рактер.
4.5. Последовательность расчета решетки профилей по методу Вознесенского-Пекина
Расчет решеток профилей выполняется о учетом стеснения потока лопастями, поэтому предварительно необходимо задать изменение мак
симальной относительной толщины профиля по радиусу. Максимальной относительной толщиной называют отношение максимальной толщины к длине хорды лопасти в данном сечении. На периферии рекомендуется принимать $тох f t - 0,025 0,035. .Максимальная относительная
толщина во втулочном сечении должна определяться на основе пове рочного расчета на прочность. Один из возможных способов расчета на прочность лопасти рабочего колеса осевого насоса изложен в мо нографии А.Н.Палира [в]. Закон изменения максимальной относитель ной толщины по радиусу в первом приближении принимается линейным,
ав последующих приближениях уточняется.
Улопаток направляющего аппарата максимальная относительная
толщина мало меняется по радиусу и обычно принимается в пределах
S n o J t e ° ’04 - °.°5 -
Целью расчета лопастной системы осевого насоса является по лучение геометрических характеристик лопасти рабочего колеса и ло патки направляющего аппарата. Геометрические характеристики строят ся на основе расчета решеток профилей • в принятых цилиндрических сечениях последовательно для лопасти рабочего колеса и лопатки на
правляющего аппарата» |
43 |
|
В пат-дог,i расчетном сечепки |
необходимо по определенным ранее |
||||
постоянным величинам £ , t ,Jbw |
, V u z , |
, Н г и |
Smax |
онруде |
|
лить кривизну скелетной дуяки д |
и |
угол |
установки |
дуъека |
|
|
в |
|
|
|
|
Расчет ведется последовательно |
по кагхдому расчетному |
сочопню |
с предварительным учетом стеснения потока профилями решетки, |
коэф |
||||
фициент стеснения подсчитывается |
по формуле |
|
|||
|
|
|
|
|
(6G) |
|
|
У _ |
^ . |
З 'т а х , |
|
|
|
|
3 |
ь Si/7A , |
|
Угол установки |
дужки |
с учетом стеснения • |
|
||
А ~ в,е*К"£?,/*) ' |
(67) |
||||
|
|||||
Относительная |
скорость с учетом стеснения |
|
|||
|
|
|
К г Г я |
(68) |
|
|
U Tоо= |
|
|||
|
|
|
|||
|
|
|
Sin Ji0 |
|
|
Кривизна скелетной дужки подсчитывается по формуле |
|
||||
|
а |
, |
|
Н т____ |
(69) |
|
J 9 |
u > z p K t w „ L |
|
||
Функция L я |
л |
|
~ ^ |
(Го >fit) определена в результате |
|
£ Л |
|
||||
|
|
|
|
|
■численного решения системы интегральных уравнении в зависимости от угла уотановки хорды и относительного шага решетки. Значение функ
ции |
L определяется интерполяцией по графику, представленному |
на |
|
рис. |
17 по известному углу установки дужки с учетом стеснения |
j$ |
|
И относительному шагу решетки TQ - i / £ в данном расчетном |
сече |
||
нии. |
|
|
|
|
Для учета влияния присоединенных вихрей на остальных дужках |
||
решетки вводится поправка к углу установки дужки |
• Фактически |
||
это |
означает, что профиль необходимо довернуть га |
некоторый угол |
|
A fic , чтобы обеспечить безударный вход в условиях |
работы профиля |
врешетке.
Впроцессе численного решения системы интегральных уравнений было установлено, что если угол установки дужки fio < 40°, угол
44
дозорота |
л fiQ |
зависит только от относительного |
шага Т0 и кривиз |
|
ны уз -. 3 этом |
случав его величина определяется |
по графику, |
пред |
|
ставленному на |
рис. 18. |
|
|
|
Золи |
угол |
установки душей fi0 > 40°, что соответствует |
решеткам |
профилен направляющего аппарата и иногда корневым решетка,, рабочего
колеса, сказывается его дополнительное влияние на угол .поворота, который теперь зависит от трех параметров Л Д , = /
В работах [I, б] представлены зависимости, позволяюаще интерполи рованием по значению кривизны скелетной дужки определить требуемый угол доворота.
3 каждом приближении угол установки дужки с учетом доворота
подсчитывается по формуле f i ^ ~ fiQ + Л fiQ . После этого выполня
ется следующее приближение: определяется функция L , подсчитывается кривизна скелетной душки fig , находится угол доворота и его значение сравнивается с предыдущим приближением. Процесс последо вательных приближений прекращается, если разница в углах доворота не превышает 1°.
После этого определяется радиус кривизны скелетной дужки
Ra - |
(70) |
#2 $ 1П f i g
Расчет решетки дужек направляющего аппарата выполняется ана логично.
По завершении расчета строятся геометрические характеристики лопасти и проверяется их плавность и монотонность. Корректировку геометрических характеристик производят изменением принятого зако на распределения густоты и максимальной1 относительной толщины про филя по радиусу. Для облегчения расчетов в последующих приближени ях можно воспользоваться разработанной на кафедре программой "Дон:1 для персональной ЭШ , составленной на языке БЕЙСИК в диалоговом режиме.
После получения удовлетворительных геометрических характе ристик выполняют ’'одевание" скелетных дужек телесным профилем. В осевых насосах для этой цели чаще всего используют стандартный симметричный профиль, относительные координаты которого приведены в табл. 2.
|
|
|
|
Т а б л и ц а |
2 |
|
5 |
ъ |
5 |
Z |
3 |
г |
|
0 |
0 |
0,15 |
0,763 |
0.7 |
0,756 |
|
0,0025 |
0,147 |
0,2 |
0,84 |
0,8 |
0,56 |
|
0,005 |
0,196 |
■ 0,3 |
0,949 |
0,9 |
0,342 |
|
0,01 |
0,294 |
0.4 |
0,998 |
0,95 |
0,222 |
|
0,025 |
0,405 |
0,45 |
1.0 |
0,97 |
0,168 |
|
0,05 |
0,516 |
.0,5 |
0,982 |
0,99 |
0,092 |
|
0,1 |
0,662 |
0,6 |
0,895 |
1.0 |
0 |
|
Здесь размеры длины даны в долях от длины дудки 6^ , которая подсчитывается по формуле
‘ г ° ’о т , т Ь ; ‘ .
Толщина задается в долях от |
максимальнол толщины профиля. |
|||
_ |
.Идя получения |
размеров |
проектируемого |
профиля нудно величины |
5 |
умножить на |
и отложить это значение |
вдоль скелетной дужки |
от входной кромки с учетом принятого масштаба. Соответствующее зна чение z умножается на половину максимальной толщины профиля в дан ном сечении clmax /2, и этим радиусом из полученной точки на ске летной дудке проводится окружность. Искомая поверхность профиля является касательной к полученным окружностям. Построение профиля рекомендуется выполнять не менее, чем в пятикратном масштабе.
5. ОСНОВНЫЕ ВОПРОСЫ АВТОМАТИЗАЦИИ ПРОЕКТИРОВАНИЯ РАБОЧИХ КОЛЕС ЦЕНТРОБЕЖНЫХ НАСОСОВ
Проектирование рабочих колес центробежных турбомашин, в том числе насосов, выполняется в настоящее время либо с использовани ем прототипов, либо по струйной теории с введением ряда эмпирических коэффициентов и рекомендаций. Решение обратной задачи для центро бежных турбомашин (компрессоров, насосов, вентиляторов) до сих пор разработано недостаточно и не дает удовлетворительных результатов при проектировании.
Поэтому на кафедре гидромашиностроения был принят подход, по которому при расчете рабочего колеса решается рад прямых задач.
47
в результате определяется оптимальны: (или несколько конкуренто способных) Еариант, подлог'/дщии дальнейшему исследованию. Для реше ния прямой задачи необходимо иметь геометрические и гидродинами ческие параметры рабочего колеса, т.е. форму меридианного сечения, поверхности тока, распределение скоростей к т.д. оти параметры оп ределяются в первом приближении на основании расчетов по струйной теории. Для автоматизации проектирования расчет рабочего колеса
по струыюй теории был запрограммирован согласно fl, 4]. Косле оп ределения геометрических и гидродинамических параметров рабочего колеса (размеров рабочего колеса, значении' скоростей, углов лопас тей и т.д.) строятся обводы рабочего колеса в меридианной плоскости. Затем производится расчет на ЭК;1 меридианного потока и поверхнос тей тока. На полученных поверхностях тока находят параметры лопаст ной системы и выполняют расчет обтекания лопастей и дисков в слое переменной толщины. На основании полученных результатов вгаолкяется расчет трехмерного пограничного слоя на ограшыивающкх поверх ностях межлопастного канала (с учетом влияния вторичных течений)
и определяются гидравлические потери трения и гпдрашшчнекий КПД рабочего колеса. На основании этих расчетов выбирается оптимальный вариант рабочего колеса. Опыт показывает, что при прочих разных условиях чем выше КПД рабочего колеса, тем выше и КПД ступени.
На кафедре гидромашшюстроения ЛПИ разработаны программы и инструкции по пользованию, ими для расчета рабочего колеса, его раз меров, формы меридианного сечения, расчета меридианного потока, обтекания лопастей и дисков идеальной жидкостью, расчета простран ственного пограничного слоя и потерь с оценкой КПД. Часть, из них используется студентами при курсовом проектировании.
Разработка ступени центробежного насоса начинается с выбора основных размеров рабочего колеса. Методика-расчета, заложенная в программе, совпадает с методикой проектирования рабочего колеса, принятой на кафедре гидромашностроения [I, 4]. Расчет основан на струйной теории, согласно схеме бесконечного числа лопастей с по следующей поправкой на их конечное число. Затем на основании рас чета течения идеальной а вязкой жидкости можно оценить гидравли ческий КПД и проверить соответствие техническому заданию теорети ческого и действительного напоров спроектированного рабочего ко леса. '
48
Входными данни.ш дня расчета на ЗЕД являются подача (расход), напор, частота вращения рабочего колеса, число его потоков и сту
пеней, а также плотность рабочей жидкости. Сначала определяют ис ходные данные для расчета рабочего колеса: подачу и напор проекти
руемой |
ступени, |
коэффициент быстроходности, ,КПД ступени, |
мощность |
|||||
на валу |
|
насоса, |
диаметр |
пала [4]. |
|
|
|
|
На |
следующем |
стало |
выбирают основные размеры рабочего колеса: |
|||||
диаметр |
|
входного |
|
отверстия JB0 , диаметр втулки d 8 , ширину канала |
||||
на входе |
и выходе |
колеса |
5 1 и 5 2 , наружный диаметр колеса-В г . |
|||||
На |
начальном |
этапе |
проектирования также определяются некото |
|||||
рые кинематические величины - переносные |
скорости |
и 1 к |
и 2 и мери |
|||||
дианные |
|
составляющие абсолютной скорости |
1ГП1 без |
учета и с учетом |
стеснения потока лопастями на входе и выходе рабочего колеса. Рас чет ведется последовательными приближениями. Следует отметить, что
всамой методике проектирования рабочих колес центробежных насосов заложен ряд коэффициентов для ■определения того или иного параметра, значения которых изменяются в довольно широких пределах. Выбор ко эффициентов производят на основании опытных данных. По алгоритму программы эти коэффициенты выбирают методом линейной интерполяции
взаьисдглости от некоторого параметра, определенного заранее. На-
пример, диаметр d B определяется по формуле
|
|
|
+ |
|
, |
(71) |
где d Q - диаметр вала. |
|
|
|
|
||
В |
программе при |
cL0 ^ |
20 мм |
к = 1,4, а при |
d 0 ^ 200 мм |
К = |
и 1,2. |
В диапазоне же |
изменения |
d a от 20 мм до |
200 мм значение К |
||
получают линейной интерполяцией. Поскольку такое |
задание К не |
|
||||
всегда |
является возможным |
(например, при проектировании ступеней |
с |
повышенными |
антикавитавдоиными качествами), то при необходимости |
в |
программе предусмотрена возможность независимого введения этого |
|
коэффициента |
по своег.дг усмотрению. |
Все полученные размеры колеса округляют до ближайшего целого значения, а диаметр вала выбирают ближайшим к размеру нормального ряда, рекомендованного стандартом.
Для характеристики потока на входе и выходе колеса находят ряд кинематических величин: относительные скорости потока, окруж ные составлялдоо абсолютной скорости vu при конзчном и бесконеч
ном числох лопастей, абсолютную окорость потока на выходе яз колеса.
49
l) |
^ |
l / ’t ’Nj |
S) |
H
Р и с . i9
Автоматизация проектирования меридианного сечения центробежного рабочего колеса была осуществлена при обеспечении плавности изменения ?леридианных скоростей потока без учета стеснения его ло пастями от входного отверстия рабочего колеса до выхода из него. По заданному закону изменения меридианной скорости вдоль радиуса R
из уравнения расхода находится ширина канала в меридианном сече нии (рис. 19) -8 (R). Затем вдоль линии АВ для построения контура
основного диска в первом приближении строятся окружности диаметром, равным ширине.канала 8 (R) (рис. 19,а). Таким образом находится
форма покрывающего диска в первом приближении. Радиус закругления
на повороте R n - Кп И 0 |
подбирается путем варьирования коэффициента |
К п , который вводится в |
ЭВМ (Кп = 0,1 * 0,2). Форма основного дис |
ка во втором приближении определяется как касательная к окружнос тям, которые в свою очередь касаются покрывающего диска. Форма ос
новного диска в |
таком |
виде принимается окончательной. Затем нахо |
|
дится радиус кривизны |
полученной средней линии в точке на входной |
||
кромке лопасти |
R ep (рис. 19,6) и в |
соответствии с рекомендацией |
|
[ 4 ] проверяется выполнение условия |
|
||
|
|
означает, что |
радиус поворота меридианного |
сечения слишком мал и следует задать другой радиус закругления, выбрав большее значение коэффициента К п .
50