CC BY
71
С. В. Визгалов, А. М. Ибраев ПЕРИОДИЗАЦИЯ РАБОЧЕГО ПРОЦЕССА КОМПРЕССОРА ТИПА РУТС И ОПРЕДЕЛЕНИЕ ЗАВИСИМОСТИ ОБЪЕМА ПОЛОСТИ ОТ УГЛА ПОВОРОТА РОТОРА
Ключевые слова: компрессор типа Рутс, объем рабочей полости.
Предложена схема периодизации рабочего процесса в компрессоре внешнего сжатия типа Рутс, учитывающая 5 стадий рабочего процесса, различающихся как характером изменения объема полостей, так и схемой протечек через щелевые зазоры. Численным интегрированием определены зависимости объемов полостей на всасывании для двух типов профилей роторов - эвольвентного с подрезкой и окружного без подрезки, которые представлены в графическом и аналитическом виде. Получено, что скорость роста объема полостей на всасывании практически не зависит от типа профиля.
Keywords: Roots type compressor, the volume of working chamber.
A scheme of periodization of the workflow in the compressor external compression type Roots, taking into account the 5 stages of the workflow, as distinct nature of the change of volume the chamber and scheme of leaks through clearance. Numerical integration are defined depending on the volume of the chamber absorption for the two types of rotors - involute with trimming and county without trimming, which are presented in graphical and analytical form. It was found that the rate of growth of cavities on the suction is practically independent of the type ofprofile.
Введение
Компрессор типа Рутс показанный на рис. 1, представляет собой двухроторную машину с прямозубыми роторами, имеющими одинаковый наружный диаметр, одинаковое число зубьев (лопастей) на ведущем и ведомом роторе, которое как правило равно 2 или 3, профили зубьев на ведущем и ведомом роторе также идентичны, что существенно упрощает технологию изготовления роторов [1]. Передача крутящего момента от ведущего к 5 3
Рис. 1 - Схема компрессора типа Рутс: ведущий ротор, 2- ведомый ротор, 3- корпус, 4- патрубок всасывания, 5- патрубок нагнетания, 6- шестерни синхронизации роторов
ведомому ротору производится с помощью синхронизирующих шестерен, насаженных на консольные участки роторов вне рабочего объёма компрессора. Всасывание и нагнетание газа осуществляется через патрубки и окна, расположенные на цилиндрической части корпуса компрессора в месте пересечения его расточек. Роторы компрессора вращаются в противоположные стороны с равными скоростями, не касаясь друг друга и стенок корпуса, что объясняется наличием зазоров между ними.
Следует отметить, что машина типа Рутс относится к компрессорам внешнего сжатия газа, в которых газ переносится со стороны всасывания на сторону нагнетания в замкнутых рабочих полостях без изменения их объема. Повышение давления газа происходит при раскрытии такой рабочей полости на стороне нагнетания и связано главным образом с обратным натеканием газа через патрубок нагнетания в рабочую полость - процессом выхлопа. Эта особенность рабочего процесса компрессора типа Рутс приводит к довольно существенному снижению его производительности и КПД, росту пульсаций давления газа и уровня шума с ростом частоты вращения роторов п0 и отношения давлений в компрессоре П=рк/рн. Следствием этого является ограничение его диапазона работы по указанным параметрам [2].
Наиболее адекватным методом расчета и исследования характеристик компрессора является составление математической модели рабочего процесса, учитывающей основные значимые факторы. Основные положения модели и расчетные уравнения даны в работе [3] в виде дифференциальных зависимостей изменения давления р и температуры Т в полости от угла поворота ротора ф - dp/dф и ^^ф.
Периодизация рабочего процесса
Система уравнений математической модели dp/dф и dT/dф дополняется зависимостями объема рабочей полости от угла ф - V=f(ф) и dV/dф, которые для компрессора типа Рутс определяются численным методом для каждой из полостей: всасывающей, изолированной рабочей полости и полости нагнетания.
Так как рабочий процесс в роторном компрессоре разбит на ряд стадий, отличающихся как характером изменения объема рабочей полоти, так и схемой протечек между по-
лостями, важно четко представлять схему периодизации этих стадий.
Выделим следующие пять стадий рабочего процесса компрессора типа Рутс:
- всасывание в парную полость увеличивающегося объема;
- выделение рабочей полости на всасывании;
- перенос «изолированной» рабочей полости постоянного объема на сторону нагнетания;
- раскрытие рабочей полости на нагнетании;
- парная полость уменьшающегося объема на нагнетании.
Одновременно в машине существуют несколько сообщающихся между собой через щелевые зазоры полостей, в которых протекают различные стадии рабочего процесса, смещенные друг относительно друга по углу поворота ротора на углы, кратные 90° (для двухлопастных машин). Принято моделирование проводить для одной рабочей полости, объемом V(ф) т.к. остальные представляют эту же полость на углах поворота ф±90°, ф±180° и т.д.
Принятая при моделировании компрессора Рутс схема периодизации рабочих процессов приведена на рис. 2, там же условно стрелками нанесена схема протечек между рассматриваемой и смежными полостями.
е г
Рис. 2 - Схема периодизации рабочих процессов и протечек
Расчет рабочего процесса начинается с периода всасывания (рис. 2а). Начальный угол поворота роторов составляет ф0=- (у<5ис+уп). В этот период рабочая полость представляет собой парную полость на всасывании (ППВ) увеличивающегося объема V (заштрихована на рис.2) которая сообщается со всасывающим патрубком (1) через окно всасывания, с полостью на угле Ф+90° объемом V+90o через радиальную (5) и торцевую (6) щели и, полостью нагнетания на угле ф+180° объемом V+180° через радиальную (2), торцевую (3) и про-
фильную (4) щели. Углы уэис, уп, ус, уаи - показаны на рис. 3.
Начиная с угла фф1= ус-уп расчет ведется для двух полостей, объемом V-! и V (включающем V-!), рис. 2Ь до момента когда давления в этих полостях не будут отличаться более чем на 3%. С этого угла и до угла фгр2=90-уэис-уп расчет ведется для рабочей полости на всасывании (РПВ) объемом VI, сообщающейся со всасывающимся патрубком (1) через окно, образованное нормалью АА1, проведенной из точки А (патрубок) на профиль ротора; с полостью
V+180° через радиальную (2) и торцевую (3) щели. Параллельно также продолжается расчет перетечек через щели 2, 3, 4, 5, 6 в парную полость V для определения коэффициента подачи компрессора.
При повороте ротора на угол фгр2=90-уэис-уп рассматриваемая полость отсоединяется от всасывающего патрубка и до угла фгрз=90+у(1и+уп существует в виде изолированной рабочей полости (ИРП) объем которой не изменяется (рис. 2с). Эта полость сообщается через радиальный (2) и торцевой (3) зазоры с полостью на нагнетании V+90 и через радиальный (1) и торцевой (3) зазоры с полостью на всасывании V-90. Характеристика СМСф для нее есть 0.
Раскрытие рабочей полости на нагнетании - РПН (рис. 2ф происходит с момента фгр3=90+у^+уп до момента ф=фгр4, при котором наблюдается выравнивание давлений в полостях объемом V и V+90o. Рассматриваемая полость, сообщается с нагнетательным патрубком через окно, образованное нормалью АА1 , подобно окну на всасывании, с полостью на всасывании V-90° через радиальную (2) и торцевую (3) щели и полостью на нагнетании V+90° через торцевую щель (3). С момента выравнивания давлений в полостях V и V+90° эти полости объединяются и рассматриваются как одна единая парная полость (рис. 2е) объем которой уменьшается (ППН). Эта полость сообщается с патрубком нагнетания (1), полостью всасывания V-180° через радиальную (2), торцевые (3) и профильные (4) щели и полостью V-90° через радиальную (5) и торцевые (6) щели. На угле фф5=180+у^<5+уп начинается раскрытие второй рабочей ячейки, (рис. 2І} объемом У-90°, на нагнетании и появляется утечка (5) в эту полость. Расчет ведется до угла фк=фгр2+90°.
Расчет объема рабочей полости
Для построения универсальной методики определения объема рабочей полости на каждой стадии процесса, независимо от кривых, очерчивающих профиль ротора удобно аппроксимировать его профиль, кроме участков подрезки кубическим сплайном и представить уравнением вида в системе координат ротора
Ур = Бр!(Хр) = 8,Хр + Ь,х2 + С;Хр + сі,
где ai, bi, ci - массивы коэффициентов сплайна,; j = 1, n,
n- число точек на профиле.
Точность аппроксимации позволяет заменить реальный профиль аппроксимированным профилем-эквивалентом, имеющем единое уравнение для всех участков.
Также необходимо иметь координаты линии зацепления в неподвижной системе X1zO1Y1z, связанной с расточкой ведущего ротора в виде
fX1z = Xlz(ф); [ Y1z = ^(ф);
Расчетная схема представлена на рис. З. Используя уравнения связи неподвижной системы X1zO1Y1z и подвижных систем координат, связанных с вращающимися роторами переводим координаты точки контакта К в систему координат первого и второго ротора XP1O1YP1 и Xp2O2Yp2.
|XP1K = -X1z) ■ sin ф-Y1z ■ cos ф;
I Yp1K = X1z ■ cosф- Y1z ■ sinф;
Г XP2K = (A - X1z)■ cosф- Y1z ■ sinф; |YP2K =-(A - X1z) ■ sin ф- Y1z ■cos ф;
Торцовое сечение парной полости на вса-
Рис. 3 - Расчетная схема
сывании (ППВ) представляет собой фигуру А2КБАіСА2. Интегрируя численным методом область, ограниченную полученными участками профиля А2К на ведомом роторе и КБЛі на ведущем и расточкой статора определяем торцевую площадь полости на угле поворота ф как сумму площадей:
^ППВ = 2СОКД2 + ^А1ОЕД1 + ^ОйЕС .
В данной статье ввиду ограниченности объема уравнения для определения ее составляющих не приводятся, полностью уравнения приведены в работе [4].
Рабочая полость существует в изолированном от патрубков всасывания и нагнетания виде (ИРП) в диапазоне углов ротора ф=л/2-у8ис-ур ...л/2+уаи+ур, ее объем неизменный, равный
где Рр, Ьр- радиус, длина ротора, кэ- коэффициент использования объема корпуса.
Результаты расчета
Результаты расчета изменения объемов полостей на всасывании для двух типов профилей представлены на рис. 4. Рассмотрены роторы с эвольвентным подрезанным профилем (кривая 1), такой профиль применяется на компрессорах типа 1А11, 2А22 и других, выпускаемых Мелитопольским компрессорным заводом (Украина) и роторы с окружным профилем без подрезки (кривая 2), характерном для вакуум-насосов типа ДВН.
V
ИРП
^ИРП ■ Lp 0,5^ ■ Rp ■ ks ■ Lp
Рис. 4 - Изменение объема рабочей полости: 1 - эвольвентный профиль (у8ис=12°, Ур=10°); 2 -окружной профиль ротора(у8ис=4,78°, ур=0°);
Как видно из графиков объем парной полости УППВ увеличивается линейно по углу поворота ротора ф, точки разрыва на графике возникают в момент отсоединения изолированной полости при фгр2=90-узис-уп. Графики рис. 4б построены для периода, изображенного на рис. 2б, когда рабочая полость начинает условно отсоединяться от парной полости, ее объем вначале растет достаточно быстро, затем рост замедляется и кривая выходит на постоянное значение объема, равное УИРП. Снижение объема происходит с момента раскрытия полости на нагнетании, т.е. с угла фгр3=90+у^3+уп рис. 2d (правые спадающие ветки графиков на рис. 4б).
В целом законы изменения объема полостей
подобны, с той разницей, что для окружного профиля объем изменяется более плавно и период существования полости в изолированном состоянии существенно меньший, чем для эвольвентного подрезанного профиля.
Закон изменения объема полости на всасывании определяет ее наполнение газом и таким образом влияет на рабочий процесс, т.е на давление, температуру, массу газа в полости. Это влияние прослеживается через скорость изменения объема, выражаемую производной dV/dф. Зависимости У(ф) рис. 4а показывают, что скорости изменения объема парной полости на всасывании практически одинаковы для двух рассмотренных типов профилей. С момента выделения рабочей полости на всасывании скорости изменения V(ф) начинают различаться. В связи с более коротким периодом существования рабочей полости в изолированном состоянии у машины с окружным профилем без подрезки будет наблюдаться более резкий рост давления в момент раскрытия полости на нагнетании, т.к. в процессе переноса изолированной полости давление в ней начинает расти за счет притечек газа через щели.
Для использования в расчетах, полученные зависимости изменения объема полости обобщены в виде следующих уравнений:
- для эвольвентного профиля:
= 1,091 + 0,0111 • ф - для ф= -22,425.. .67,575°;
Чирп
= 0,713 + 7,478 • 10-3 • ф . для
^ИРП
ф= 36,28.67,575°;
- для окружного профиля:
= 1,01 + 0,0108 • ф - для ф= -4,78 .85,22°;
^ИРП
= 0,79 + 5,186• 10-3 •ф - для ф= 51.85,22°.
^ИРП
Эвольвентный профиль рассматривался
применительно к геометрии машины 1А11, окружной - ДВН50.
Литература
1. А.Г. Головинцев, В.А. Румянцев, В.И. Ардашев, Ротационные компрессоры. М.: Машиностроение, 1964. - 315с.
2. С.В. Визгалов, А.М. Ибраев, И.И. Шарапов, Коэффициент подачи роторного компрессора внешнего сжатия. Вестник Казан. технол. ун-та. - 2010, №10, с. 535-541.
3. А.М. Ибраев, С.В. Визгалов, И.И. Шарапов, Выбор геометрических параметров ротора при проектировании шестеренчатого компрессора внешнего сжатия. Вестник Казан. технол. ун-та. - 2011, №6, с. 187-193.
4. С.В. Визгалов, Дисс. канд. технич. наук. Казан. гос. технол. ун-т. Казань, 2004. - 212 с.
© С. В. Визгалов - доц. каф. холодильной техники и технологии КНИТУ, [email protected], А. М. Ибраев - доц. каф. холодильной техники и технологии КНИТУ, [email protected].
