Например, Бобцов

К расчету протечек в рабочей части винтового компрессора

Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование» УДК 621.514

№ 3, 2014

К расчету протечек в рабочей части винтового компрессора
Докукин Владимир Николаевич, dockvn@mail.ru Д-р техн. наук, проф. Пронин Владимир Александрович
Университет ИТМО 191002, Санкт-Петербург, ул. Ломоносова, 9

В статье рассмотрены математическая модель и методика расчета течения сжимаемой среды через щели в рабочей части винтового маслозаполненного компрессора, с учетом подвижности их стенок. Анализ различных типов щелей, проведенный авторами работы показал, что в большей части щелей винтового компрессора имеет место встречное движение стенок и рабочего вещества, а стенки на линии контакта и со стороны нагнетания имеют попутное движение. Вследствие того, что при совпадении векторов скоростей протечки увеличиваются, возникает необходимость их учета, что и было сделано в предложенной математической модели. Благодаря предложенному методу расчета протечек повышается точность расчета коэффициента подачи винтового маслозаполненного компрессора. Адекватность предложенной методики подтверждена экспериментальными данными, путем сравнения коэффициента подачи, рассчитанного с учетом предлагаемой методики с коэффициентом подачи полученным в ходе испытания экспериментального компрессора.
Ключевые слова: коэффициент подачи ВКМ, протечка рабочего вещества, классификация щелей.

To leakage calculation in working part of a screw compressor
Dokukin V.N., D.Sc., prof. Pronin V.A.
University ITMO 191002, Russia, St. Petersburg, Lomonosov str., 9

In the article it is considered a mathematical model and calculation method of refrigerant-oil mixture flow through gaps in the working part of the oil screw compressor with movability of their walls. Analysis of different types of gaps, which was made by the authors, has shown that reciprocal movement of the walls and working fluid takes place in most of the gaps of the screw compressor, but walls on the contact line and from the side of discharge has the same direction movement. Due to the fact of velocity vectors matching leakage increases, therefore it becomes necessary to take them into account, which was in the proposed mathematical model. Due to the proposed method of calculation of leakage, the accuracy of volume coefficient calculation of oil screw compressor improves. The adequacy of the proposed method has been confirmed by experimental data, by comparing the calculated volume coefficient taking into account the proposed methodology with volume coefficient obtained during the test of experimental compressor.
Keywords: volume coefficient, leakage of working fluid, gaps classification.
В винтовых маслозаполненных компрессорах (ВКМ), также, как и в других компрессорах объемного принципа действия различают протечки внешние и внутренние.

17

Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование»

№ 3, 2014

Протечки через концевые уплотнения (сальники) вала ротора, в виде утечек рабочего вещества или подсоса атмосферного воздуха называются внешними. Протечки рабочего вещества через зазоры между деталями в рабочем пространстве ВКМ со стороны нагнетания на сторону всасывания называются внутренними протечками и оказывают существенное влияние на рабочий процесс.
Причем, протечки рабочего вещества в момент всасывания, называются утечками и влияют в основном на коэффициент подачи, а в меньшей степени на эффективный КПД и мощность двигателя. Утечки происходят из полостей с повышенным давлением в полости всасывания, а также в полости, соединенные в этот момент с камерой всасывания. Заполнение полости всасывания рабочим веществом утечек ведет к уменьшению объема, вновь поступающего (всасываемого) рабочего вещества. Температура пара рабочего вещества утечек более высокая чем температура вновь всасываемого пара, что приводит к повышению температуры и снижению плотности газовой смеси находящейся в полости всасывания. Соответственно масса пара вновь всасываемого рабочего вещества уменьшается.
Влияние различных объемных потерь на производительность ВКМ оценивается коэффициентом подачи:
=1- у- б- w- м- вс (1),

где у – объемные потери при утечках рабочего вещества через зазоры в парные полости винтов;
б – объемные потери при заполнении полости балластным рабочим веществом; w – объемные потери вследствие подогрева рабочего вещества называются коэффициентом подогрева;
м – объемные потери от масла, которое поступает в полости винтов и занимает часть объема; вс – объемные потери вследствие сопротивления во всасывающем тракте и окне всасывания. Известно, что основными видами объемных потерь являются потери, которые связаны с утечками и балластным рабочим веществом [1-3]. Причем потери через утечки значительно превосходят балластные потери. Поверхности, образующие зазоры между роторами винтового маслозаполненного компрессора, движутся. Величина протечек зависит от длины зазоров, от угла поворота винтов, от перепада давления в щели, скорости точек поверхности винта, параметров потока перед щелью [4]. Следовательно, течение в щелях нестационарное. В классической газовой динамике рассматривается течение совершенного газа, для которого в качестве термического уравнения состояния применяется уравнение Клайперона P=ρRT (2), а в качестве калорического уравнения состояния – зависимость i=cpT (3), причем удельная теплоемкость при постоянном давлении предполагается неизменной. Среда, сжимаемая ВКМ, отличается по своим термодинамическим свойствам от совершенных газов. Поэтому представляет интерес рассмотрение течений реальных газов с конкретными уравнениями состояния через щели в ВКМ. Анализ различных типов щелей показал, что в большей части щелей ВКМ имеет место встречное движение стенок и рабочего вещества, а стенки на линии контакта и с торца нагнетания имеют попутное движение. Вследствие того, что при совпадении векторов скоростей протечки увеличиваются, возникает необходимость их учета, что и было сделано в предложенной математической модели. При вращении роторов винтового маслозаполненного компрессора изменяются параметры щелей и векторы скоростей взаимного движения поверхностей их образующих. Таким образом движение сжимаемого рабочего вещества через такие зазоры следует считать нестационарным.

18

Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование»

№ 3, 2014

Течение рабочего вещества во всех рассматриваемых щелях будем считать ламинарным. Характерную для каждого зазора высоту δi будем считать малой по сравнению с его глубиной l. Также
считаем малым комплекс i / l Re .
Теплообменом потока рабочего вещества с ограничивающими его твердыми поверхностями пренебрегаем.
Для упрощения рассмотрим движение сжимаемой среды по отношению к зубьям ведущего (ВЩ) винта, что позволит упростить граничные условия на поверхностях зубьев и выявить влияние вращения зубьев ВЩ винта на течение в щелях. Затем от рассмотрения движения в системе координат, вращающейся вместе с ВЩ винтом, перейдем к новым координатным системам, связанным непосредственно с зубом ВЩ винта или с его поверхностями.
Рассмотрим исходную систему основных уравнений относительного движения. К основным уравнениям, описывающим движение сжимаемой среды через щели, будем относить уравнение неразрывности, три уравнения движения и уравнение сохранения энергии.
Цилиндрическая система координат r, θ, z, которая вращается с постоянной угловой скоростью w1 вместе с ведущим (ВЩ) винтом ось z направлена вдоль оси винта.
Уравнение неразрывности для нестационарного течения сжимаемой сплошной среды имеет вид:

t

( r

Wr )

r

( W)

( Wz )
z

Wr 0 (4) r

где ρ – плотность;
t – время;
Wr ,W ,Wz – радиальная, окружная и осевая составляющие относительно скорости потока по
отношению к ведомому (ВМ) винту. Уравнения нестационарного движения сплошной среды во вращающейся с постоянной угловой
скоростью w2 системе цилиндрических координат имеет следующий вид:

dWr W dt r

1

P r

w22r

2w2W

Fr (5)

dW WrW dt r

1P r 2w2Wr F (6)

dWz dt

1P z Fz (7)

где P – давление;

Fr , F , Fz – величины, характеризующие влияние вязкой среды.
Уравнение сохранения энергии может быть записано в форме справедливой для реального газа

в виде:

di dP dt dt Nдис di ( gradT ) (8),

19

Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование»

№ 3, 2014

где i – удельная энтальпия; T – температура; Nдис – мощность сил трения, диссипируемая в теплоту.
Если воспользоваться уравнением неразрывности (4) и уравнениями движения (5) и (7), то уравнение (8) можно представить в виде:

( i* P) 1 tr

( r

rWr i* )

r

( W i*)

Wr Fr W F

Wz Fz Wzw22r Nдис

( z

Wzi* )

(9),

di ( gradT )

где i* – полная удельная энтальпия,

i*

i

1 2

(Wr2

W2

Wz2 ) (10)

На базе полученной математической модели была разработана методика расчета протечек.

Последовательность расчета следующая:

1) Устанавливается дискретность точек по углу поворота ВЩ винта и угловой координа-

те на кромке зуба.

2) Цикл интегрирования осуществляется по , в ходе которого производятся все вычисле-

ния, изложенные ранее. 3) Определяется объем полости канавки за зубом ВЩ винта, соответствующий объему рабо-
чей полости сжатия при данном положении зуба:

где S

V

max
z S2

R

d (11)

z1

- площадь сечения канавки винта, перекрытая зубом,

R - расстояние от оси винта до центра тяжести сечения зуба.
Для зуба окружного профиля:

S

r32 2

sin 2 2

(12)

R

R

r3

4 3

2

sin3 sin2

cos (13)

где - половина центрального угла при зацеплении ВЩ винта с ВМ винтом:

arccos r1 1 cos r3

(14)

20

Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование»

4) Вычисляются:
перепад давления - P между полостями перед и за зубом ВЩ винта;
плотность паромасляной смеси в полости сжатия ;

температура смеси T ,

n

P

Pвс

V V

max

1 (15)

вс

V - max мм V

(16)

n1

T

V Tвс V

max

(17)

№ 3, 2014

где Pвс , вс ,Tвс - параметры рабочей среды при всасывании;
V max - полный объем полости сжатия;
n - показатель политропы процесса сжатия в рабочей полости;
м - плотность масла;
м - концентрация масла в паромасляной смеси. 5) Производится цикл интегрирования по , в ходе которого определяются геометрические и кинематические характеристики щели. Расчет ведется в соответствии с системой координат. При расчетах определяются проекции характерных скоростей относительно зуба ВЩ винта:

Ux

w
1

z2 z1

r3

(18)

U z w1 R1 1 cos r3 cos (19)

Относительная скорость точки на поверхности винта:

U

U

2 x

U

2 z

(20)

Глубина щели L определяется конфигурацией вершины зуба винта.

Минимальная высота щели равна:

h0 r0 1 esin

(21)

гдеr0 - расстояние между центрами сечений зуба, e - эксцентриситет. При расчетах задается численное значение эксцентриситета e в пределах от 0,1 до 0,9 и расчеты
выполняются для ряда значений в этой области.
Значения W0 - средний расход при фиксированном угле и различных значениях определялись методом хорд, после чего подсчитывался удельный расход по дуге d :

21

Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование»

№ 3, 2014

q r3h0 W0 (22)
Далее удельный расход суммировался по углу и по углу . Составляющая коэффициента подачи ВКМ, соответствующая протечкам через рассмотренную щель, определялась выражением:

11

q, всV

dd
max

(23)

После классификации всех характерных щелей в рабочей части винтового компрессора и определения протечек рабочей среды через них на основании изложенной методики, определяется коэффи-
циент подачи винтового маслозаполненного компрессора. Предложенная авторами математическая модель апробирована путем сравнения данных числен-
ного эксперимента с данными полученными экспериментально [5-8].

Список литературы

1. Канышев Г.А., Чистяков Ф.М. Коэффициент подачи винтового фреонового маслозаполненного компрессора // Холодильная техника: М., 1979. вып. 12. с. 12-16. 2. Пронин В.А., Верболоз А.П. Оценка влияния подвижности стенок щелей на протечки компримируемой среды в винтовом однороторном компрессоре (ВКО) // Вестник международной академии холода. 2012. № 1. с. 31-33. 3. Сакун И.А. Винтовые компрессоры: Основы теории, методы расчета, конструкции. – Л.: Машиностроение,1970. – 400 с. 4. Пронин В.А. Тенденции развития компрессорной техники // Научный журнал НИУ ИТМО. Серия: Холодильная техника и кондиционирование. 2014. № 1. С. 9. 5. Докукин В.Н., Емельянов А.Л., Носков А.Н.. Результаты испытаний маслозаполненного винтового компрессора малой производительности в высокотемпературных режимах // Вестник международной академии холода. 2009. № 1. с. 6-8. 6. Докукин В.Н., Носков А.Н., Муштуков Д.Ю. Результаты испытаний холодильного масозаполненного винтового компрессора малой производительности с новым профилем зубьев в высокотемпературных режимах // Вестник международной академии холода. 2011. № 1. с. 17-21. 7. В.Н. Докукин. Испытание винтового маслозаполненного компрессора в режиме теплового насоса // Тр. XV Международной научно-технической конференции по компрессорной технике / ЗАО «НИИтурбокомпрессор им. В. Б. Шнеппа». – Казань, 2011. – Т. 1. – с. 74-81. 8. Емельянов А.Л., Горбатов К.М., Гаранов С.А. Гибридная испарительно-компрессионная установка кондиционирования воздуха // Вестник международной академии холода. 2013. № 4. с. 34-37. 9. Шамеко С.Л., Докукин В.Н. Современные методы энергоаудита турбоагрегатов // Холодильная техника и кондиционирование. 2012. № 1.

22

Научный журнал НИУ ИТМО. Серия «Холодильная техника и кондиционирование»

№ 3, 2014

References
1. Kanyshev G.A., Chistyakov F.M. Koeffitsiyent of giving of the screw freon mas-lozapolnenny compressor // Kholodil'naya tekhnika: M., 1979. vyp. 12. P. 12-16. 2. Pronin V.A., Verboloz A.P. Otsenka of influence of mobility of walls of cracks on leakages of the compressed environment in the screw one-rotor compressor (VKO) // Vestnik mezhdunarodnoi akademii kholoda. 2012. № 1. P. 31-33. 3. Sakun I.A. Screw compressors: Theory bases, methods of calculation, design.. – L.: Mashinostroenie,1970. – 400 s. 4. Pronin V.A. Tendencies of development of compressor equipment // Nauchnyi zhurnal NIU ITMO. Seriya: Kholodil'naya tekhnika i konditsionirovanie. 2014. № 1. P. 9. 5. Dokukin V.N., Emel'yanov A.L., Noskov A.N.. Results of tests of the maslozapolnenny screw compressor of small productivity in the high-temperature modes // Vestnik mezhdunarodnoi akademii kholoda. 2009. № 1. P. 6-8. 6. Dokukin V.N., Noskov A.N., Mushtukov D.Yu. Results of tests of the refrigerating masozapolnen-ny screw compressor of small productivity with a new profile of teeths in the high-temperature modes // Vestnik mezhdunarodnoi akademii kholoda. 2011. № 1. P. 17-21. 7. V.N. Dokukin. Test of the screw maslozapolnenny compressor in the mode of the thermal pump// Tr. XV Mezhdunarodnoi nauchno-tekhnicheskoi konferentsii po kompressornoi tekhnike / ZAO «NIIturbokompressor im. V. B. Shneppa». – Kazan', 2011. – T. 1. – P. 74-81. 8. Emel'yanov A.L., Gorbatov K.M., Garanov S.A. Gibridnaya vaporizing and compression установ air conditionings // Vestnik mezhdunarodnoi akademii kholoda. 2013. № 4. P. 34-37. 9. Shameko S.L., Dokukin V.N. Sovremennye metody energoaudita turboagregatov // Kholodil'naya tekhnika i konditsionirovanie. 2012. № 1.

23