Машиностроение

DOI 10.5862/JEST.231.19
УДК 621.515 (075.8)

Ю.Б. Галеркин, А.А. Дроздов

ОПТИМИЗАцИЯ цЕНТРОбЕжНОЙ КОМПРЕССОРНОЙ СТУПЕНИ

С ОСЕРАДИАЛЬНЫМ РАбОЧИМ КОЛЕСОМ

Yu.B. Galerkin, A.A. Drozdov

OPTIMIZATION OF A CENTRIFUGAL COMPRESSOR STAGE

WITH AXIAL-RADIAL IMPELLER

Представлены обновления 6-й версии программ Метода универсального моделирования для

расчета характеристик центробежных компрессорных ступеней с уточненным моделированием
пространственных осерадиальных колес. Разработанная математическая модель
идентифицирована по результатам модельных испытаний и CFD-расчетов. Произведена
оптимизация центробежной компрессорной ступени с осерадиальным рабочим колесом при
помощи разработанной программы. В результате оптимизации осерадиального рабочего колеса
КПД ступени повышен на 1%. Результаты оптимизации проверены CFD-расчетами в программе
NUMECAFINE/TURBO. Разработаны рекомендации по выбору размеров осерадиальных
рабочих колес.
ЦЕНТРОБЕЖНЫЙ КОМПРЕССОР, МЕТОД УНИВЕРСАЛЬНОГО МОДЕЛИРОВАНИЯ, КПД, КОЭФФИЦИ-
ЕНТ НАПОРА, ОСЕРАДИАЛЬНОЕ РАБОЧЕЕ КОЛЕСО.

The 6-th updated version of the universal modeling method for centrifugal compressor stage calculation
is presented, adjusted for modeling 3D impeller performances. The developed mathematical model was
identified by the results of model tests and CFD-calculations. The efficiency definition error at a design
point was 0,39 %. The efficiency definition error at five flow rate points (except a point of the maximum
flow rate) was 0,88 %. The centrifugal compressor stage was optimized by a new mathematical model.
As a result of the optimization, stage efficiency increased by 1%. The optimization results have been
verified by CFD-calculations in the NUMECA FINE/TURBO program. Design recommendations for
3D impeller design are offered.
CENTRIFUGAL COMPRESSOR, UNIVERSAL MODELING METHOD, EFFICIENCY, HEAD СOEFFICIENT,
3D IMPELLER.

Оптимальное газодинамическое проектиро- ров, агрегатов турбонаддува ДВС и т.д. Наиболее

вание центробежных компрессоров – важная и широкой областью применения центробежных
актуальная задача, что связано с широким рас- компрессоров в России является добыча и транс-
пространением этого типа компрессоров и боль- портировка природного газа. Так, газотран-
шим количеством энергии, затрачиваемым на их спортная система ОАО «Газпром» включает не-
привод. сколько сотен компрессорных цехов и
Актуальность работы обусловлена тем, что несколько тысяч газоперекачивающих агрегатов
осерадиальные рабочие колеса (ОРК) широко (ГПА) [1].
востребованы в проектах газотурбинных двига- Значительная часть применяемых компрес-
телей летательных аппаратов и малой энергети- соров ГПА не соответствует современным тре-
ки, промышленных центробежных компрессо- бованиям из-за низкого КПД или изменивших-

179
 Научно-технические ведомости Cанкт-Петербургского государственного политехнического университета. 4(231)’2015

ся параметров газопроводов. Это приводит к и безлопаточными диффузорами (БЛД) [5–7]. В

необходимости замены ГПА целиком или, что 6-й версии программ Метода универсального
более эффективно, только проточной части в моделирования сделаны существенные усовер-
существующем корпусе с использованием смен- шенствования для уточненного расчета ступеней
ных проточных частей (СПЧ). с осерадиальными рабочими колесами.
Эти факторы накладываются на объявленный
Президентом РФ курс на импортозамещение вы- Уточненный расчет осерадиальных рабочих
сокотехнологичной продукции, к которой отно- колес
сится и компрессорная техника. Необходимо В программах 4-й версии МУМ описание
проектировать и освоить производство тех ком- формы осерадиальных рабочих колес сугубо при-
прессоров, которые ранее закупались у иностран- ближенное. Задаются только относительный
ных производителей. Развитие современных вы-
диаметр втулки Dвт , диаметр входа в РК D0 ,
числительных систем, а так же расширение
высота лопаток на выходе b2 . Эмпирические
представлений о процессах, протекающих в про-
формулы устанавливают осевую протяженность
точной части, изменили подходы к проектирова-
РК, радиусы кривизны, внутренние и внешние
нию и расчету характеристик компрессоров [2].
ограничивающие поверхности меридионально-
Современной тенденция – повышение бы-
го контура. Положение входной кромки счита-
строходности центробежных компрессорных сту-
ется совпадающим с контрольным сечением «0»
пеней ГПА, что приводит к применению осера-
на входе в РК.
диальных рабочих колес [3]. При проектировании
В программах 6-й версии форма проточной
таких ступеней следует принимать во внимание,
части ОРК в меридиональной плоскости задает-
что ОРК имеют сложную геометрическую форму
ся более подробно. На рис. 1 представлена схема,
лопаток и меридионального контура. Это услож-
поясняющая задаваемые размеры ОРК в про-
няет их оптимальное проектирование, делает не-
грамме 6-й версии.
обходимым анализ влияния значительно больше-
го количества геометрических параметров. Помимо указанных выше D0 , Dвт , b2 , β л1 и
Цель настоящей работы – создание такой ин- β л2 , задаются также осевая протяженность ко-
женерной методики, которая в состоянии быстро
и надежно выполнить оптимизацию осерадиаль-
2
ного рабочего колеса и идентификацию усовер-
шенствованной математической модели, а также
произвести оптимизационное исследование
центробежной компрессорной ступени с ОРК,
спроектированной на перспективные параметры
ГПА.
Для быстрого и достоверного определения
характеристик компрессора и их конструирова-
ния профессор СПбПУ Ю.Б. Галёркин создал
метод оптимального газодинамического про-
ектирования, который является результатом
обобщения и анализа большого объема экспе-
риментальных данных, накопленных в СПбПУ.
Комплекс соответствующих компьютерных про-
грамм получил название Метода универсально-
го моделирования (МУМ) [4].
Предшествующие версии программ Метода
универсального моделирования эффективно ис-
пользовались при разработке новых компрессо-
ров по заданиям промышленности на базе сту-
пеней с радиальными рабочими колесами (РК) Рис. 1. Схема ОРК по разработанной
математической модели

180
 Машиностроение

леса Lm , радиусы кривизны основного Rh и по- определяются на средней линии межлопаточно-

крывающего Rs дисков. Положение входной (
го канала S 0 s = S 0h = S 0 ср :)
кромки лопатки в меридиональной плоскости So ОРК
определяется отношением меридиональной про- S 0 s = S 0h = S 0 ср = =
тяженности лопатки к общей меридиональной π 2
D2
протяженности средней линии ОРК 4 (2)
lm л = lm л lm ОРК . Входная кромка лопатки пред-  4 δл 
полагается перпендикулярной к средней линии 
( )
= lm  2 1 + Dср − z РК .
π sin β л ср 
ОРК в этой точке. Кроме того, задается угол 2
наклона покрывающего диска ОРК. Эти раз- В математическую модель ОРК было внесе-
меры полностью описывают меридиональную но множество других изменений, представлен-
форму большинства ОРК. ных в работах [9–11].
Для расчета площади лопатки необходимо
знать изменение лопаточного угла по длине ло- Идентификация математической модели
патки. Для этого был проведен анализ несколь- для высокорасходных ступеней
ких вариантов осерадиальных колес, спроекти- с осерадиальными рабочими колесами
рованных в СПбПУ. В результате линейной Ранее была идентифицирована математиче-
аппроксимации получен следующий закон из- ская модель 6-й версии для ступеней с радиаль-
менения углов по длине лопатки: ными рабочими колесами и различными типами
начальный участок с lm тек от 0 до 0,45 lm л , на диффузоров [10]. Полученный в результате иден-
котором принимается линейное увеличение угла тификации единый набор эмпирических коэф-
от β л1 до β лмакс =1,35 л2; фициентов и программа 6-й версии обеспечива-
средний участок с lm тек от 0,45 до 0,75 lm л и ют определение КПД на расчетном режиме со
постоянным лопаточным углом л макс; средней погрешностью 0,86 %. По пяти рабочим
конечный участок с lm тек от 0,75 до 1,0 lm л , точкам (без точки максимального расхода) сред-
на котором принимается линейное уменьшение няя погрешность составляет 1,22 %.
угла от л макс до л2. Математическая модель МУМ включает в
Заданная таким образом форма ОРК позво- себя эмпирические коэффициенты, применяе-
ляет более точно по сравнению с предыдущими мые только для расчета ступеней с ОРК. Для их
версиями программы определить площадь огра- уточнения использованы результаты испытания
ничивающих поверхностей и лопаток ОРК. Для двухзвенной ступени с полуоткрытым ОРК, вы-
этого меридиональная длина лопатки разделя- полненные профессором СПбПУ А.М. Симоно-
ется на участки, для каждого из которых опре- вым. Ступень с осерадиальным рабочим коле-
деляется его длина ∆lmi , высота лопатки bi и сом, безлопаточным диффузором и сборной
камерой спроектирована на расчетные параме-
диаметр Di . Площадь всей лопатки складывает-
тры: Φ м =м ьшь8;м Ψ м =м ьш74; Mu = 0,78. Для
ся из суммы площадей каждого участка с учетом
сопоставления расчетов с экспериментом были
изменения угла лопатки рабочего колеса по его
использованы характеристики двухзвенной сту-
длине. Площадь лопатки определяется формулой
пени ОРК+БЛД, так как данных о размерах
i ∆lmi bi сборной камеры у авторов нет.
Sл = ∑ , (1)
sin β лi cos χ Так как имеющихся в распоряжении экспе-
0
где – угол навала лопатки. риментальных данных (шесть эксперименталь-
Величина угла навала лопатки зависит от ных точек характеристики) недостаточно для
многих факторов и изменятется по длине и вы- проведения идентификации, а других экспери-
соте лопатки. В упрощенной постановке, при- ментальных данных у авторов не было, то были
меняемой в данном случае, угол навала прини- применены также результаты CFD-расчетов.
мается постоянной величиной, равной 20°[8]. Сопоставление экспериментальных данных с
Площади поверхностей основного и покры- CFD-расчетами показали хорошее совпадение
вающего дисков принимаются одинаковыми и на расчетном режиме работы [12, 13].

181
 Научно-технические ведомости Cанкт-Петербургского государственного политехнического университета. 4(231)’2015

Для идентификации были спроектированы шмΨ

две ступени с ОРК на следующие параметры:
Φ м=мьшэь5;мΨ м=мьш56. Варианты различают-
ся числом лопаток ОРК (15 и 22 шт), соответ-
ственно были изменены входные и выходные
углы и форма лопаток в целом. Одно измеренное
и два рассчитанных значения КПД сопоставле-
ны с расчетом по математической модели с не-
откорректированными и откорректированными
коэффициентами (табл. 1).
Таблица 1

Сопоставление КПД измеренных и рассчитанных Рис. 2. Характеристики двухзвенной ступени

по CFD-программам и новой математической модели ОРК+БЛД сплошные линии – расчет; точки –
эксперимент (А.М. Симонов, ЛПИ)
Эксперимен- CFD- CFD-
КПД тальные данные расчет расчет
А.М. Симонова ступени ступени На рис. 2 экспериментальная характеристи-
для ОРК+БЛД с ОРК с ОРК ка двухзвенной ступени ОРК + БЛД сопостав-
zРК = 15 zРК = 22 лена с расчетом по МУМ.
эксп
87,0% – – Средняя погрешность расчета КПД по пяти
– 86,4 % 84,5 % точкам (за исключением максимального рас-
NUMECA
хода) – 0,88 %. В расчетной точке погрешность
до
МУМ
корректи-
равна 0,39 %. Положительный результат позво-
ровки эм- ляет производить расчетный анализ ступеней с
85,3 % 82,3 % 80,1 % ОРК при помощи разработанного пакета про-
пирических
коэффици- грамм.
ентов
после Оптимизация ОРК
МУМ
корректи- На втором этапе работы Методом универ-
ровки эм- сального моделирования была спроектирована
86,8 % 86,4 % 85,0 %
пирических
центробежная компрессорная ступень промежу-
коэффици-
ентов
точного типа с ОРК на расчетные параметры,
соответствующие перспективной первой ступе-
ни ГПА: Φ м=мьшэь5;мΨ м=мьш56; Mu = 0,7; Reu =
Расчеты ступеней при помощи МУМ пока-
= 5600000; k = 1,4.
зали необходимость учета влияния трения на
Новая математическая модель МУМ позво-
лопатках ОРК на потери КПД. Так как число
ляет изучить влияние ряда геометрических па-
лопаток влияет на величину нагрузки на лопат-
ки, то корректируется и эмпирический коэффи- раметров ОРК: меридионального размера Lm
циент, определяющий скорость в начале задней ОРК; диаметра D0 входа в ОРК; положения
поверхности лопатки. входной кромки лопатки lm л ; радиусов кривиз-
По расчету с уточненными значениями эм- ны поверхностей основного и покрывающего
пирических коэффициентов гидравлический дисков ОРК (Rs , Rh ); высоты лопаток на выходе
КПД ступени с ОРК с 15 лопатками на расчетном b2 ; числа лопаток zРК.
режиме работы равен 86,4 %, т.е. в точности со- План оптимизационных расчетов, целью ко-
ответствует значению, рассчитанному по торых является получение наивысшего КПД
NUMECA Fine/Turbo. Для ступени с ОРК с 22 ступени, предлагает следующее:
лопатками различие между КПД, рассчитанным 1) Соотношения размеров ОРК для исходной
по математической модели и по программе
ступени – Lm = 0,3, zРК = 15, D0 = 0,6 , lm л = 0,792,
NUMECA, составляет 0,5 % (85,0 и 84,5% соот-
Rs = 0,215, Rh = 0,257;
ветственно).

182
 Машиностроение

2) оптимизируются четыре группы ОРК с скорости в каналах. Одновременно происходит

Lm , равными 0,26, 0,28, 0,30, 0,35. Это сделано уменьшение в 7,88 раз потерь смешения (с 0,0335
с целью найти оптимальный размер Lm ; до 0,00425). Потери трения на ограничивающих
3) в пределах каждой группы варьируется поверхностях ОРК закономерно остаются неиз-
число лопаток, относительный диаметр входа и менными. Оптимальным является вариант с
положение входной кромки (именно в перечис- zРК=15, он имеет на 0,8 % больший КПД по срав-
ленной последовательности). нению с наименее эффективным вариантом с
Радиусы кривизны основного и покрываю- zРК= 22.
щего дисков Rs , Rh оказывают малое влияние на Далее оптимизировался относительный диа-
КПД ступени и не оптимизируются. метр входа ОРК. Согласно рекомендациям [4]
Для всех рассмотренных вариантов были по- размеры входа в рабочее колесо выбираются с
добраны лопаточные углы входа и выхода с це- ′ на вхо-
учетом условия минимума скорости w1min
лью обеспечения заданных условного коэффи- де в межлопаточные каналы у покрывающего
циента расхода и коэффициента теоретического диска:
напора. 1
2

Первой производилась оптимизация ступени 2  Ф расч  3

D0 w min = Dвт + 23  ε ' τ  . (3)
с ОРК при Lm = 0,3. Были сопоставлены вари- 1 1
анты с 10, 12, 15, 18 и 22 лопатками. Полученные Потери КПД определяются уровнем w1′ и
значения конструктивных параметров, коэффи- коэффициентом потерь ОРК. Канал с диаметром
циентов потерь, скоростей и КПД ступени – в входа по формуле (3) не всегда обеспечивает ми-
табл. 2. нимальный коэффициент потерь. В качестве
Таблица 2 параметра оптимизации используется коэффи-
циент AD [4]:
Параметры и размеры вариантов ступени с различным D0
числом лопаток zРК AD = . (4)
D0 w 0min
zРК , , w1s′ м+ ,
л1 л2 м м
Для анализируемого ОРК D0 w min = 0,613 , т.е.
град град %
10 27,9 67,9 0,717 0,0597 0,0335 84,8
у исходного ОРК коэффициент AD = 0,979.
Относительный диаметр входа изменяется с
12 28,7 57,9 0,717 0,0674 0,0225 85,2 шагом 0,01 в пределах от 0,53 до 0,63. В табл. 3
15 29,7 51,1 0,719 0,0801 0,0127 85,3 приведены значения конструктивных параме-
18 30,5 48 0,721 0,0930 0,0073 85,0 тров, коэффициентов потерь, скоростей и КПД.
22 31,4 46,2 0,724 0,1100 0,0043 84,5 Максимальному КПД соответствуют значе-
ния AD = 0,897–0,913. При этом КПД ступени на
Максимальный КПД ступени достигается 0,8 % больше, чем у исходного варианта. Ско-
при числе лопаток zРК =15. Уменьшение числа рость на входе у оптимальных вариантов не ми-
лопаток рабочего колеса требует увеличения вы- нимальная, однако за счет меньшей площади
ходного угла л2 из-за роста нагрузки [4]. Боль- лопаток потери трения меньше, чем в вариантах
шая нагрузка увеличивает фактор диффузорно- с большими AD. В вариантах с AD > 1 имеет место
сти, что приводит к более раннему отрыву и значительный рост потерь как за счет увеличения
росту потерь смешения. Увеличение числа лопа- уровня скоростей, так и за счет увеличения пло-
ток приводит к росту смоченной поверхности и щади лопаток.
росту потерь трения. Максимальному КПД сту- Из двух равноценных с точки зрения КПД
пени соответствует оптимальное число лопаток вариантов для дальнейших расчетов принято
ОРК. Потери трения на лопатках превалирует значение D0 = 0,56, так как скорость на входе в
над остальными потерями. Увеличение числа лопаточную решетку w1s′ в этом варианте мень-
лопаток с 10 до 22 приводит к более чем двукрат- ше, чем в равноценном ему по величине КПД
ному увеличению потерь трения – в 2,137 раза варианте с D0 = 0,55.
(с 0,0453 до 0,0968). Росту потерь способствует Далее сопоставлялись 10 вариантов с раз-
загромождение потока и вызванный этим рост личным положением входной кромки – lm л в

183
 Научно-технические ведомости Cанкт-Петербургского государственного политехнического университета. 4(231)’2015

Таблица 3

Параметры и размеры вариантов ступени с различными диаметрами входа ОРК D0 .

, ра , ра AD
w1s′ м
м+ м
,%
D0 л1 л2

ьш5я 4ьш5 5э ьш864 ьш744 ьшь55я ьшьэ8э 85ш9

ьш54 я8ш8 5э ьш88э ьш7я6 ьшь578 ьшьэ56 86шь
ьш55 я7шэ 5э ьш897 ьш7яь ьшь6ь7 ьшьэ4ь 86шэ
0,56 35,4 51 0,914 0,724 0,0640 0,0131 86,1
ьш57 яяш9 5э ьш9яь ьш7юю ьшь67я ьшьээ8 86шь
ьш58 яюш4 5э ьш946 ьш7э9 ьшь7эя ьшьээ7 85ш8
ьш59 яэ 5э ьш96я ьш687 ьшь756 ьшьээ9 85ш6
ьш6ь ю9ш7 5э ьш979 ьш7э9 ьшь8ьэ ьшьэю7 85шя
ьш6э ю8ш4 5э ьш995 ьш7э9 ьшь857 ьшьэяэ 85шь
ьш6ю ю7шя 5э эшьээ ьш7юю ьшь9эя ьшьэя7 84ш8
ьш6я ю6шэ 5э эшью8 ьш7ю5 ьшь977 ьшьэ5э 84шю

Таблица 4

lë
Параметры и размеры вариантов ступени с различными положениями входной кромки лопаток ОРК m
шм
э ю
шм w1s′ м
мчм м
шмс
lm л

ьш6э яьш7 59ш5 ьш7яь ьшь5юю ьшьээю 86шь

ьш64 яэш4 57шя ьш7ю6 ьшь54ю ьшьэээ 86шэ
ьш67 яюшю 55ш9 ьш7ю4 ьшь56ь ьшьэээ 86шэ
ьш7 яя 54шя ьш7юя ьшь579 ьшьээь 86шю
ьш7я яяш8 5я ьш7юэ ьшь6ьь ьшьээь 86шю
0,75 34,3 52 0,722 0,0615 0,0109 86,2
ьш76 я4ш6 5ю ьш7юя ьшь6э9 ьшьэюю 86шэ
ьш79ю я5ш4 5э ьш7ю4 ьшь64ь ьшьэяь 86шэ
ьш8э я5ш9 5ьш5 ьш7ю6 ьшь65э ьшьэяэ 86шь
ьш84 я6ш8 49ш5 ьш7яь ьшь67я ьшьэя5 85ш9

пределах от 0,61 до 0,84. В табл. 4 приведены В результате проведенной оптимизации сту-

конструктивные параметры, коэффициенты по- пени промежуточного типа с Lm = 0,3 максимум
терь, безразмерная скорость на периферии ло- КПД на расчетном режиме – 86,2 % достигнут
патки, КПД сопоставленных вариантов. при следующих размерах проточной части:
Максимальный КПД ступени 86,2% дости- D0 = 0,56; lm л = 0,75; zРК=15; Rs = 0,215; Rh = 0,257.
гается при lm л = 0,7–0,75. По мере приближения Для групп вариантов ступени со значениями
положения входной кромки лопатки к началу Lm 0,35, 0,28 и 0,26 проведено аналогичное пред-
ОРК увеличиваются потери трения на лопатках ставленному выше расчетное исследование. Ре-
(в 1,31 раза), потери на ограничивающих поверх- зультаты незначительно отличаются количе-
ностях (в 1,21 раза). ственно, а качественные закономерности

184
 Машиностроение

идентичны. Геометрические размеры и КПД а) б)

оптимальных вариантов ступеней с различной
осевой протяженностью представлены в табл. 5.
Периферия
Максимальным КПД обладает ступень с Втулка
Lm = 0,28, который на 0,5% выше чем у худшего
варианта с Lm = 0,35. В результате оптимизации
КПД ступени удалось повысить по сравнению с Средняя линия
исходным вариантом на 1,0%.
Таблица 5

Геометрические соотношения, лопаточные углы и КПД

оптимизированных вариантов ОРК Рис. 3. Изменение лопаточных углов на трех линиях
z , , шм тока по длине межлопаточного канала (а) и внешний
Lm D0 lm л л1 л2
вид лопаточной решетки (б)
ра ра с
ьшя5 ьш55 эя ьш69 я5шю 6ьш4 85ш8
Сравнение характеристик КПД исходной и
ьшяь ьш56 э5 ьш75 я4шя 5ю 86шю оптимизированной ступени (рис. 4) показывает,
0,28 0,56 15 0,69 33,7 54,7 86,3 что оптимизация позволила повысить КПД сту-
ьшю6 ьш56 э5 ьш69 я5ш4 54ш5 86шю пени на всех режимах работы ступени, кроме
режима максимального расхода. КПД ступени
Окончательные выводы по результатам про- на расчетном режиме повысился на 1,02 %; по
веденного исследования можно сделать после всей характеристике (9 режимов работы) КПД
проверки полученных результатов путем числен- повысился в среднем на 0,88 %.
ного эксперимента (CFD-расчет). Сопоставление расчетных данных по CFD-
программе и по новой математической модели,
Проверка результатов оптимизации а так же результаты оптимизации ступени по-
центробежной компрессорной ступени зволяют сделать следующие выводы:
оценки КПД ступени на расчетном режиме
Полученные инженерным методом проекти-
по 6-й версии математической модели с иденти-
рования результаты оптимизации нуждаются в
фицированным набором эмпирических коэф-
проверке. Для этого исходная ступень и ступень,
фициентов и по CFD-методам совпадают с точ-
полученная в результате оптимизации с осевой
ностью 0,25 %;
длиной Lm = 0,28 и оптимизированными раз- подтверждена возможность эффективно оп-
мерами, указанными в табл. 5, были рассчитаны тимизировать размеры и форму ОРК и ступени
по CFD-программе NUMECA Fine/Turbo. Ло- в целом при использовании новой математиче-
паточные решетки ОРК были спрофилированы ской модели.
согласно принципам, выработанным и под- В целом по результатам исследования можно
твержденным в научно-исследовательской ла- заключить следующее:
боратории «Газовая динамика турбомашин»
оптимальный размер входного диаметра D0
ОНТИ СПбПУ. Внешний вид оптимизирован-
определяется исходя из того, что параметр AD
ного ОРК и распределение углов по длине ло-
должен лежать в пределах 0,897–0,914;
патки представлены на рис. 3.
Сопоставление рассчитанных по новой ма- оптимальная величина размера lm л лежит в
тематической модели и по CFD-программе зна- пределах 0,69–0,75; это соответствует положе-
чений гидравлического КПД на расчетном ре- нию входной кромки между радиальным и осе-
жиме показывает хорошее совпадение: вым положениями на криволинейном участке
= 87,71%; NUMECA = 87,46 %. Это можно ОРК;
МУМ
считать хорошим результатом, подтверждающим излишне длинные или короткие ОРК с осе-
правильность разработанной новой математи- вой протяженностью Lm равными 0,35 и 0,26
ческой модели и корректность проведенной менее эффективны (∆ =0,5 – 0,1 %). Оптималь-
идентификации. ной длиной является Lm = 0,28.

185
 Научно-технические ведомости Cанкт-Петербургского государственного политехнического университета. 4(231)’2015

ьш96 2
ьш95
ьш94
ьш9я
ьш9ю
ьш9э
ьш9 4

ьш89
ьш88
ьш87
ьш86
ьш85
ьш84
ьш8я
ьш8ю
ьш8э
ьш8
ьш79
ьш78
ьшь75 ьшь8 ьшь85 ьшь9 ьшь95 ьшэ ьшэь5 ьшээ ьшээ5 ьшэю ьшэю5 ьшэя Φ

Рис. 4. Характеристики КПД ступени, двухзвенной ступени и ОРК

Выводы Приложение
Разработана новая, 6-я версия инженерного Условные обозначения: b – высота лопатки;
метода оптимального проектирования центро- D – диаметр; k – показатель адиабаты; Lm – осе-
бежных компрессорных ступеней, в том числе с вая протяженность ОРК; l –длина лопатки;
осерадиальными рабочими колесами. Произ- Mu – число Маха; R – радиус кривизны; Reu –
ведена идентификация математической модели число Рейнольдса; S – площадь; w – относитель-
ОРК по результатам модельных и численных ная скорость потока; z – число лопаток; л – угол
экспериментов. между касательной к средней линии лопатки
Выполнена оптимизация рабочего колеса с рабочего колеса и обратным окружным направ-
условным коэффициентом расхода 0,105 и ко- лением; ε – отношение плотностей; – коэф-
эффициентом теоретического напора 0,56. Вы- фициент полезного действия; м– толщина ло-
работаны рекомендации по выбору осевой дли- патки; м –м коэффициент стеснения; м – угол
ны ступени, положению входной кромки наклона образующей покрывающего диска;
лопатки, размеру входного диаметра. Результаты Φ – условный коэффициент расхода; – угол
оптимизации проверены путем проведения чис- навала лопатки; м– коэффициент потерь; ΨТ –
ленного эксперимента в программе NUMECA. коэффициент теоретического напора.
Разработанная новая версия инженерного Подстрочные индексы: 0, 1, 2, 4 – параметры
метода (Метод универсального моделирования) – потока и размеры контрольных сечений на гра-
эффективный способ решения задачи газодина- нице элементов проточной части; вт – втулка;
мического проектирования высокорасходных л – лопатка; макс – максимальное; о ограничи-
центробежных компрессорных ступеней с осе- вающие поверхности; ОРК – осерадиальное ра-
радиальными рабочими колесами. бочее колесо; расч – расчетное значение; РК ра-
бочее колесо; ср – среднее значение;
Работа выполнена при поддержке гранта Пре- см – смешение; тр – трение; эксп – эксперимен-
зидента Российской Федерации для молодых канди- тальный; m – меридиональный; h – втулка;
датов наук МК-7066.2015.8. min – минимальное; s – периферия.

186
 Машиностроение

Надстрочные индексы: c = c u2 , b = b D2 – ный размер отнесен к характерному линейному

надстрочная черта означает, что скорость отне- размеру (наружному диаметру рабочего колеса);
сена к характерной скорости вращения, линей- T* – параметры торможения; w′ – скорость с
учетом стеснения потока.

СПИСОК ЛИТЕРАТУРЫ

1. Галеркин Ю.Б., Солдатова К. В. Моделирова- ческий проект, результат модельных испытаний). //

ние рабочего процесса промышленных центробеж- Компрессорная техника и пневматика. 2014. № 8.
ных компрессоров. Научные основы, этапы разви- С. 19–24.
тия, современное состояние: Монография. СПб.: 8. Селезнев К.П., Галеркин Ю.Б. Центробежные
Изд-во Политехнического университета, 2011. 327 с. компрессоры // Л.: Машиностроение, 1982.
2. Труды научной школы компрессоростроения 9. Galerkin Y., Drozdov A., Soldatova K. Modern
СПбГПУ / Под ред. Ю.Б. Галеркина. М.: Изд-во state of the universal modeling for centrifugal compressors.
СПбГПУ, 2010. [Electronic] / Y. Galerkin, // International Conference
3. Галеркин Ю.Б. Рекстин А.Ф., Дроздов А.А., on Numerical Methods in Industrial Processes. Paris.
Солдатова К.В. Особенности газодинамического 2015.
проектирования центробежных компрессоров для 10. Галеркин Ю.Б., Дроздов А.А. Моделирование
газовой промышленности. // Труды 19 Междуна- газодинамических характеристик центробежных
родного симпозиума «Потребители-производители компрессорных ступеней с осерадиальными рабочи-
компрессоров и компрессорного оборудования». ми колесами // Научно-технические ведомости
СПб., 2015. С. 112-126. СПбГПУ. : Наука и образование 2014. №3(202).
4. Галеркин Ю.Б. Турбокомпрессоры. Рабочий С. 45–53.
процесс, расчет и проектирование проточной части. 11. Галеркин Ю.Б., Дроздов А.А., Солдатова К.В.
М.: ООО «Информационно-издательский центр Новая версия математической модели для расчета
«КХТ», 2010. 596 с. С ил. ISBN 978–5–9902032–1–1. газодинамических характеристик центробежных
5. Галеркин Ю.Б., Солдатова К.В., Дроздов А.А. компрессоров // Труды 16 Международной научно-
Уточнение алгоритма расчета параметров потока в технической конференции по компрессорострое-
центробежной компрессорной ступени // Научно- нию. Санкт-Петербург. 2014. С. 296–305.
технические ведомости СПбГПУ. 2010. № 4(110). 12. Гамбургер Д.М., Епифанов А.А., Гаев А.В.
С. 150–157. Численное исследование трехмерного течения вяз-
6. Галеркин Ю.Б., Солдатова К.В., Дроздов А.А. кого сжимаемого газа в рабочем колесе центробеж-
Развитие компьютерных программ Метода универ- ного компрессора // Научно-технические ведомости
сального моделирования 1-го уровня // Труды 15 СПбГПУ.: Наука и образование. 2009. №2(78)
Международной научно-технической конференции С. 76–82.
по компрессорной технике. Том 1. Казань, 2011. 13. Мельников В., Прокусов А. Использование
С. 276–285. программного комплекса FlowVision при расчете
7. Галеркин Ю.Б., Рекстин А.Ф., Солдатова К.В., элементов проточной части турбокомпрессоров в
Дроздов А.А. Высокоэффективный одноступенча- ОАО «СКБТ» // САПР и графика. 2005. №4.
тый полнонапорный компрессор ГПА (газодинами- С. 92–96.

REFERENCES

1. Galerkin Y.B., Soldatova K.V. Modelirovaniye 3. Galerkin Y.B., Rekstin A.F., Drozdov A.A.,
rabochego protsessa promyshlennykh tsentrobezhnykh Soldatova K.V. Osobennosti gazodinamicheskogo
kompressorov. Nauchnyye osnovy, etapy razvitiya, proyektirovaniya tsentrobezhnykh kompressorov dlya
sovremennoye sostoyaniye. Monografiya. [Gas dynamic gazovoy promyshlennosti. [Features of gas-dynamic
modeling of industrial centrifugal compressors. Scientific design of centrifugal compressors for the gas industry].
bases, steps of development, modern state.] SPb: Izd-vo Trudy 19 Mezhdunarodnogo simpoziuma «Potrebiteli-
Politekhnicheskogo universiteta, 2011. 327p. (rus) proizvoditeli kompressorov i kompressornogo
2. Galerkina Y.B. Trudy nauchnoy shkoly oborudovaniya». Sankt-Peterburg. 2015. P. 112–126.
kompressorostroyeniya SPbGPU. [Proceedings of the (rus)
scientific school of compressor SPbSTU] M.: Izd-vo 4. Galerkin Y.B. Turbokompressory. Rabochiy
SPbGPU, 2010. (rus) protsess, raschet i proyektirovaniye protochnoy chasti.

187
 Научно-технические ведомости Cанкт-Петербургского государственного политехнического университета. 4(231)’2015

[Turbo compressors]. M.: OOO «Informatsionno- 10. Galerkin Y.B., Drozdov A.A. Modelirovaniye
izdatelskiy tsentr «KKhT». 2010. 596 S. ISBN 978–5– gazodinamicheskikh kharakteristik tsentrobezhnykh
9902032–1–1. (rus) kompressornykh stupeney s oseradialnymi rabochimi
5. Galerkin Y.B., Soldatova K.V., Drozdov A.A. kolesami [Modeling of gas dynamic performances of
Utochneniye algoritma rascheta parametrov potoka v centrifugal compressor stages with 3D impellers].
tsentrobezhnoy kompressornoy stupeni. [Clarification Nauchno-tekhnicheskiye vedomosti SPbGPU.: Nauka i
of the algorithm for the flow parameters calculating in a obrazovaniye. SPb., 2014 . № 3(202). S. 45–53. (rus)
centrifugal compressor stage]. Nauchno-tekhnicheskiye 11. Galerkin Y.B., Drozdov A.A., Soldatova K.V.
vedomosti SPbGPU. 2010. № 4(110). S. 150–157. (rus) Novaya versiya matematicheskoy modeli dlya rascheta
6. Galerkin Y.B., Soldatova K.V., Drozdov A.A. gazodinamicheskikh kharakteristik tsentrobezhnykh
Razvitiye kompyuternykh programm Metoda kompressorov. [The new version of the mathematical
universalnogo modelirovaniya 1-go urovnya. [The model for calculating gas dynamic performances of
development of computer software Universal Method centrifugal compressors]. Trudy 16 Mezhdunarodnoy
simulation 1st level]. Trudy 15 Mezhdunarodnoy nauchno-tekhnicheskoy konferentsii po
nauchno-tekhnicheskoy konferentsii po kompressornoy kompressorostroyeniyu. Sankt-Peterburg. 2014. S. 296–
tekhnike. Tom 1. Kazan. 2011. P. 276–285. (rus) 305. (rus)
7. Galerkin Y.B., Rekstin A.F., Soldatova K.V., 12. Gamburger D.M. Yepifanov A.A. Gayev A.V.
Drozdov A.A. Vysokoeffektivnyy odnostupenchatyy Chislennoye issledovaniye trekhmernogo techeniya
polnonapornyy kompressor GPA (gazodinamicheskiy vyazkogo szhimayemogo gaza v rabochem kolese
proyekt, rezultat modelnykh ispytaniy) [Highly efficient tsentrobezhnogo kompressora. [Numerical investigation
full-head, single-stage linear gas compressor (gas- of three-dimensional viscous compressible gas in the
dynamic design, the result of model testing)]. impeller of a centrifugal compressor]. Nauchno-
Compressors and pneumatics. 2014. № 8. P. 19–24. (rus) tekhnicheskiye vedomosti SPbGPU. Nauka i obrazovaniye.
8. Seleznev K.P. Galerkin Y.B. Tsentrobezhnyye 2009. № 2(78): S. 76–82. (rus)
kompressory. [Centrifugal compressors]. L.: 13. Melnikov V., Prokusov A. Ispolzovaniye
Mashinostroyeniye, 1982. (rus) programmnogo kompleksa FlowVision pri raschete
9. Galerkin Y.B., Drozdov A.A., Soldatova K.V. elementov protochnoy chasti turbokompressorov v OAO
Modern state of the universal modeling for centrifugal «SKBT». [Using the software package FlowVision for
compressors. [Electronic] International Conference on flow calculation of the turbochargers elements in
Numerical Methods in Industrial Processes. Paris. 2015. «SKBT»]. SAPR i grafika. 2005. № 4. S. 92–96. (rus)

СВЕДЕНИЯ Об АВТОРАХ/AUTHORS
е И мЮ чм –м м м м м м щ
м м« м м »м щ м
м м м .м
э95ю5эшм шм .м щ шм м .шмю9.
Eщmail:м[email protected]
GALerkINмYurii B.м–мPeterмtheмGreatмSt.мPetersburgмPolytechnicмUniversity.
ю9мPolitechnicheskayaмSt.шмSt.мPetersburgшмэ95ю5эшмRussia.
Eщmail:м[email protected]
оз овм чм–м м щ м м м
м .м
э95ю5эшм шм .м щ шм м .шмю9.
Eщmail:м[email protected]
droZdovмAleksandr A.м–мPeterмtheмGreatмSt.мPetersburgмPolytechnicмUniversity.
ю9мPolitechnicheskayaмSt.шмSt.мPetersburgшмэ95ю5эшмRussia.
Eщmail:м[email protected]

188