CC BY
73
Р. Р. Якупов, В. Н. Налимов, И. Г. Хисамеев
К ВОПРОСУ ПРОФИЛИРОВАНИЯ РОТОРОВ ВИНТОВОГО КОМПРЕССОРА
Ключевые слова: численные методы, профиль ротора, винтовой компрессор.
Описаны численные методы применяемые при профилировании роторов винтового компрессора и критериальном анализе этих профилей, такие как интерполяция кубическим сплайном и методы решения системы нелинейных уравнений. Показаны достоинства и недостатки этих методов, а также варианты уменьшения влияния данных недостатков.
Key wordbi: numerical methods, rotor shape, screw compressor.
Numerical methods have been described, used for shaping rotors of the screw compressor and for criterial analysis of the shapes obtained such as cubic spline interpolation and non-lineas equation system solution methods. Advantages and drawbacks of these methods have been shown as well as the way how to minimize the effect of the above shortcomings.
Наиболее сложным элементом конструкции винтовых компрессоров являются профильные поверхности роторов. Они в наибольшей мере определяют энергетические характеристики компрессоров, такие как КПД и коэффициент подачи. Многие фирмы, в том числе и ЗАО «НИИтурбокомпрессор им В. Б. Шнеппа», уделяют большое внимание совершенствованию профилей роторов [1]. Открытие новых возможностей при разработке профилей так же объяснятся развитием вычислительной техники, которая позволяет использовать численные методы для решения ряда задач, возникающих при профилировании роторов.
Ротора в винтовом компрессоре исполняют роль пары вращающихся поршней и зубчатой пары, поэтому профили роторов винтового компрессора должен отвечать определённым условиям, а именно:
1. Линия контакта должна изолировать рабочие полости друг от друга и от полости всасывания и нагнетания.
2. При вращении роторов линия контакта должна перемещаться в пространстве, обеспечивая уменьшение парной полости и сжатие газа.
3. Геометрия роторов должна удовлетворять основной теореме зацепления.
4. Профиль ротора должен обеспечивать технологичность изготовления (позволять сформировать профиль фрезы), при этом необходимо плавное сопряжение кривых, из которых состоит профиль ротора. Наиболее широко для описания данных кривых используется параметрическая форма записи функции.
X = ^ (1)
У = ^2 ()
ли
Для плавного сопряжения кривых необходимо решить следующую систему уравнений.
X = )
у = Г2 (ti.Pi Рп) Х = ^3 (2’^1’...’^п )
У = /ф^,..^ ) (2)
С^2 (11)/ / (12 )/
/а^ /а12
/ (М/ = С/ (12 )/ .
. А /^2
где X, у - координаты точки сопряжения; 11, 12 - параметры кривых в точке сопряжения; Р, Б - дополнительные параметры, определяющие форму используемых кривых. Часть дополнительных параметров может быть рассчитана при решение системы уравнений введением дополнительных условий при сопряжении, либо могут варьироваться для анализа и подбора оптимальных значений.
При использовании сложных кривых решение данной системы аналитически вызывало большие трудности, что приводило к тому, что для профилирования роторов применялись «простые» кривые, такие как циклоида, окружность, отрезок прямой. С развитием компьютерной техники, появилась возможность решать данную систему численными методами, что позволяет использовать в профилировании «сложные» кривые. Для решения системы нелинейных уравнений используются следующие методы: метод Ньютона, метод простых итераций, метод спуска. При этом начальное приближение легко выбирается из построения данных кривых на плоскости.
Вторым ограничением на использование сложных кривых при профилировании роторов винтового компрессора является необходимость поиска сопряженной кривой, отвечающей теореме зацепления. Использование численных методов позволяет получить сопряженный профиль в виде координат точек и значений касательных в этих точках.
Так через произвольную точку А (рис. 1) проводим нормаль к профилю п-п, и находим току пересечения нормали с начальной окружностью К. Определяем угол а, который определяет на какой угол повернется ротор, когда данная точка войдет в зацепление. Зная угол а и переведя координаты данной точки из подвижной системы координат ро-
тора в неподвижные координаты, мы получим координаты точки линии зацепления, а переведя координаты в подвижную систему координат ответного ротора, получим координату ответного профиля.
Рис. 1
филя
Определение координат ответного про-
Для определения угла а необходимо решить следующую систему.
(х - хА )
У к + Уа
х2 + У2 = г2н (3)
а = - агСап^—^
Современные теоретические профиля роторов винтовых компрессоров состоят из множества кривых, изменяя форму которых можно влиять на эффективность профиля. Для получения оптимальных характеристик необходимо анализировать профиль ещё на этапе проектирования. Можно выделить несколько критериев по которым производят анализ и оптимизацию профилей.
• Относительная длина линии контакта;
• Относительная площадь треугольной щели;
• Относительная величина защемленного объема;
• Относительная площадь впадины между зубьями;
• Угол давления.
При анализе данных критериев часто встает задача определения площади сложных фигур. При использовании в качестве исходных данных координат профилей, данная задача решается интерполяцией профиля и численным интегрированием [2]. Для уменьшения трудоемкости можно воспользоваться формулой расчета площади произвольного п-угольника:
(4)
Рп =(хіУ2 + х2Уз +•
- (уіх2 + у2х3 + ••• +
(5)
•• + хпУ1)-Упх1)
где Х|, у, - координаты вершин п-угольника.
Как показывает проверка, при достаточно полном задании профиля в виде точек, погрешность данного метода составляет сотые доли процента (рис. 2). При этом данный метод лишен ограничений, накладываемых интерполяцией сплайном.
5 = 220 мм2 , Аб = 0,12 мм2
Рис. 2 - Расчет площади многоугольника
Хоть замена профильных кривых ломаной и позволяет облегчить процесс анализа профилей, но в некоторых случаях это приводит в увеличению погрешности до недопустимых значений. Примером может служить расчет зазора в зацеплении (рис.к 3). Погрешность составляет до 20 %.
Рис. 3 - Определение профильного зазора в зацеплении
В этом случае целесообразнее применять интерполяцию сплайном, которая позволяет описывать плавные кривые с большой точностью, обеспечивая непрерывность самой кривой и её первой и второй производной. Непрерывность второй производной позволяет получать очень гладкие кривые и при этом является дополнительным условием для определения коэффициентов кубических парабол, из которых составляется сплайн. Но если при стыковки двух кривых не обеспечивается условие непрерывности второй производной, то при интерполяция сплайном в данном месте будет вносить дополнительную погрешность (рис. 4).
Рис. 4 - Погрешность при интерполяции сплайном
Чтобы исключить данную погрешность, необходимо при интерполяции сплайном в качестве дополнительного условия использовать, вместо непрерывности второй производной, значение первой производной в узловых точках. Это позволит избежать дополнительных погрешностей в местах сты-
ковки и увеличит точность расчета первой производной, которая используется при профилировании фрез.
Численный методы позволяют разработать универсальные методики расчета и анализа профилей, которые ускорят процесс разработки и внедрения новых профилей для роторов винтовых компрессоров.
Литература
1. Хисамеев, И.Г. Анализ направлений в разработке профилей винтовых компрессоров / И. Г. Хисамеев, Р. Р. Якупов, В. Н. Налимов // Вестник Казан.технолог.ун-та. - 2011. - № 16. - С. 245-249.
2. Чекушкин, Г.Н. Численный метод профилирования и анализа геометрии рабочих камер роторных компрессоров / Г. Н. Чекушкин, М. С. Хамидуллин, В. Н. Налимов // Казань. - 1988. -
С. 1-25.
3. Хисамеев, И.Г. Создание винтовых Компрессорных установок нового поколения/ И.Г. Хисамеев, М.Г. Абд-реев, Ю.А. Паранин, М.Т. Садыков, Р.Р. Яку-пов//Вестник Казан. технол. ун-та. - 2011. - №17. - С. 100-203.
© Р. Р. Якупов - асс. каф. холодильной техники и технологий КНИТУ, [email protected]; В. Н. Налимов - нач. отд.расчетов и испытаний отделения винтовых компрессоров ЗАО «НИИтурбокомпрессор»; И. Г. Хисамеев - д-р техн. наук, проф., зав. каф. холодильной техники и технологий КНИТУ, [email protected].
