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

книги / Расчет центробежных и осевых насосов лопастных гидромашин

..pdf
Скачиваний:
30
Добавлен:
12.11.2023
Размер:
6.44 Mб
Скачать

Угол поворота потока в решетке направляющею аппарата

А 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

,fig,J50) .

дозорота

л 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

Соседние файлы в папке книги