CC BY
15
УДК 621.514.54
Математическое описание рабочих процессов винтового компрессора на фазе регулирования производительности
А.В. АНУФРИЕВ, В.И. ПЕКАРЕВ
СПбГУНиПТ
The change in the section area of a relief port of a screw compressor, the volumetric capacity of which is regulated by a slide valve, is considered. The standard slide valve and the slide valve with changed rims located under the angles, equal to the turning angles of the screw lines of the corresponding rotors are compared. The analysis of change of the relief port area by the turning angles of the driving rotor shows that if the valve with changed rims is used, then this area in the phase of control is always larger, i.e., the gas dynamic resistance is less. Besides, the mathematical model of the compressor BX-130 (R22) shows that compressor efficiency with the changed slide valve during capacity control is 2 - 8 % higher, than if the standard slide valve is used.
Большой комплекс работ, проведенных на кафедре холодильных машин и низкопотенциальной энергетики СПбГУНиПТ, показал значительное снижение эффективности работы винтового компрессора при регулировании холодопроизводительности золотниковым способом. Эффективный КПД при уменьшении производительности до 50 % снижается на 30 — 40 % |3, 4]. Это происходит главным образом вследствие уменьшения геометрической степени сжатия 8Г компрессора при снижении холодопроизводительности, а также из-за появления дополнительных газодинамических потерь при выталкивании “лишнего” количества пара из парной полости в камеру всасывания.
Процессы в винтовом компрессоре при регулировании холодопроизводительности можно условно разбить на четыре этапа. Первый — процесс всасывания, при котором весь объем парной полости заполняется паром рабочего вещества. Второй — процесс выталкивания части пара рабочего вещества, не участвующего в процессе сжатия, из парной полости обратно в камеру всасывания. Третий этап — непосредственно процесс сжатия и, наконец, четвертый — нагнетание (выталкивание пара из парной полости в камеру нагнетания).
Для создания адекватной математической модели процессов, происходящих в винтовом компрессоре при регулировании его объемной производительности, необходимо знать площадь сечения перепускного окна, через которое “лишний” пар перепускается из парной полости в камеру всасывания при любом угле поворота ведущего ротора. Однако размер и конфигурация этого окна зависят не только от угла поворота ротора, но и непосредственно от положения золотника регулятора. В связи с этим необходимо получить зависимости плошади перепускного окна /рсг от угла поворота ведущего ротора <р, при фиксированных значениях положения золотника.
Рассмотрим более подробно изменение площади перепускного окна экспериментального компрессора ВХ-130 при положении золотника, соответствующем 40 % объемной производительности, т.е. положение регулятора дол -
(2)
жно быть таким, чтобы сжатие началось в момент, когда объем парной полости соответствует заданному значению величины объемной производительности.
При производительности 40 % объем парной полости в начале сжатия должен составлять 40 % от ее полного объема Щ,:
40Г
\¥. =-----2-. (И
о,4 юо и;
В общем случае
ХІЇ
IV =
’ 100
где х — заданное значение объемной производительности компрессора, %.
Для того чтобы правильно определить положение золотника, необходимо определить угол поворота ротора в момент, когда объем парной полости составит х процентов от полного объема.
Для экспериментального компрессора ВХ-130 полный объем парной полости равен И/0= 429 см\ значит, при производительности 40 % сжатие должно начаться в момент, когда объем полости составит IV,, 4= 171,6 см3. Такому значению соответствует угол поворота ведущего ротора ф] = 224,6°.
Далее, воспользовавшись зависимостью между углом поворота ротора и продольным перемещением винтовой линии 2[, определяем положение золотника регулятора производительности:
2 _ У Л /гі
р 360 ’
(3)
где Zp — продольное положение золотника от торца всасывания, мм;
/г, — осевой шаг ведущего винта, мм.
Изменение площади перепускного окна при регулировании холодопроизводительности компрессора с помощью золотникового регулятора с кромками, расположенными под углами, равными углам наклона винтовых линий соответствующих винтов (измененный золотник) I I |,
Рис. 1. Зависимость площади сечения перепускного окна /рсг от угла поворота ведущего ротора ф( при регулировании измененным золотником
дано на рис. 1. Изменение конфигурации перепускного окна происходит в три этапа. На первом этапе имеет место линейное увеличение площади перепускного окна (участок 1 — 2 на рис. 1). В конце этого этапа линейная зависимость сменяется выпуклой (участок 2 — 3 на рис. 1), т.е. производная с!/р^/с1 ^р, уменьшается, что объясняется тем обстоятельством, что кромки торца всасывания расположены под прямым углом к продольной оси компрессора, а винтовые линии роторов — нет, а значит, одинаковым приращениям г/ф! будут соответствовать разные приращения с!/рег
На втором этапе площадь перепускного окна остается постоянной, о чем свидетельствует горизонтальный участок 3—4зависимости/рег =Л ср,) на рис. 3. На этом этапе парная полость уже отделилась от торца всасывания, но еще не достигла кромок золотника. При этом площадь перепускного окна равна полной проекции парной полости на горизонтальную плоскость, а значит, площадь перепускного окна имеет максимально возможное значение (как будет показано ниже, зависимость^,,^/(ср,) может и не иметь горизонтального участка).
Третий этап — уменьшение площади сечения перепускного окна. На рис. 1 этому этапу соответствует участок 4—7. На большей его части (5 — 6) площадь сечения перепускного окна уменьшается линейно, однако на участках 4—5и 6 — 7имеет место нелинейность, обусловленная наличием поперечного участка кромки золотника регулятора со стороны ведомого винта. Несколько иначе выглядит зависимость/рег-/( ср^ при регулировании холодопроизводительности стандартным золотником, кромки которого со стороны всасывания расположены под прямым углом к продольной оси компрессора (рис. 2). На начальном этапе, как и в первом случае, имеет место линейное увеличение площади сечения перепускного окна (участок / — 2 на рис. 2). Далее зависимость У|Ч., =_Д (р|) становится нелинейной (участок 2—4). Нелинейность, как и ранее, объясняется различными углами кромок торца всасывания и наклона винтовых линий винтов. Естественно, это обстоятельство зависит только от геометрии винтов и является общим для золотниковых регуляторов обоих видов. Как видно из рис. 2, функцияУрС,-—У( ф)) имеет экстремум в точке 3 на участке 2 — 4 в отличие от случая с измененным золотником (рис. 1), когда максимальному значению функции соот-
,
2/
л
у
1
ч. 6
0 50 100 150 200 Ф,,*
Рис. 2. Зависимость площади сечения перепускного окна / от угла поворота ведущего ротора ср, при регулировании стандартным золотником
ветствовал целый участок 3 — 4. Это объясняется тем, что при таком положении золотника, которое показано на рис. 1, в момент, когда задние по направлению движения парной полости зубья отделились от торца всасывания, передние уже достигли кромки регулятора. Вследствие этого максимальное значение площади сечения перепускного окна меньше, чем полная площадь проекции парной полости на горизонтальную плоскость. Следует отметить, что при объемной производительности 30 % от полной и ниже горизонтальный участок все же будет иметь место, однако его протяженность по углу поворота ведущего ротора будет значительно меньшей, чем при регулировании измененным золотником. Участок 3—4 (рис. 2) характеризуется нелинейным уменьшением площади перепускного окна. Этот участок более протяженный, чем для измененного золотника (см. 4 - 5 на рис. 1), так как кромки стандартного золотника располагаются перпендикулярно оси компрессора, а у измененного золотника присутствует лишь небольшой поперечный участок.
Далее следует линейное уменьшение площади перепускного окна на участке 4 — 5. При достижении задними по направлению движения зубьями винта кромок золотника зависимость вновь становится нелинейной (участок 5 — 6 на рис. 2).
Таким образом, проведенный анализ показывает, что в случае использования золотника с измененными кромками имеет место увеличение проходного сечения на всем протяжении фазы регулирования, кроме начального участка, когда передние по ходу движения винтов зубья парной полости еще не достигли прямой кромки стандартного золотника. Таким образом, в случае использования измененного золотника при перепуске пара имеет место уменьшение величины газодинамических сопротивлений.
Термодинамическое состояние хладагента в рабочей полости компрессора при регулировании производительности может быть описано с помощью следующих дифференциальных уравнений [2|:
_Г1+0 ЛЦ +
с!(р с IVсо (И V ) с/ф
+ Х'р-х ,<Л.х , - Х/(Лых _, - 'Чег
гН
сЮ с/И'
—--сор--------+
Л с/ф
сГГ_
с/ф сутО)
+Х('»* ■ “'К'.* , + /П'((Лх , -вщш., - Ср« )]•
(5)
гдер — давление в рабочей полости;
Ф - угол поворота ротора;
I - фактор сжимаемости газа;
К — газовая постоянная;
с„— удельная теплоемкость при постоянном объеме; ]У — объем парной области;
Т— температура;
(2 — теплота;
/ — энтальпия;
(7— массовый расход;
V — объем в рабочей полости; ш — масса.
Для математического описания процесса перепуска рабочего вещества из парной полости в камеру всасывания вводятся следующие допущения:
Стечение газа одномерное изоэнтропное;
/для описания течения в каждый момент времени могут быть использованы зависимости для установившегося потока газа.
Тогда массовый расход перепускаемого хладагента можно записать в следующем виде:
<Л«г = аер/р*>/2РО>-Ас) ’ <6>
где (X - коэффициент расхода;
е - коэффициент расширения;
/рс[ - площадь проходного сечения перепускного окна, м2;
р — плотность рабочего вещества в парной полости, кг/м1;
Р и Рве — соответственно давление в парной полости и в камере всасывания, Па.
Вследствие отсутствия данных по стационарным продувкам перепускных окон винтовых компрессоров значения коэффициента расхода определяли путем анализа индикаторных диаграмм экспериментального компрессора ВХ-130с помощью разработанной математической модели 121.
Для этого уравнение (4) решалось относительно величины расхода перепускаемого пара (7рег:
Ар с\}У о) (10
(, _ Дф г/? ск
(|+Й
/
сПУ V-,
а
(7)
При этом величина Дд/Дф, записанная в конечных разностях, определялась непосредственно из опытной индикаторной диаграммы 13). Д;иее вычислялось значение ко-эффициента расхода перепускного окна:
а =-----, 6рег ■ (8)
Зависимость коэффициента расхода от числа Рейнольдса, полученная в результате вычислений, изображена на рис. 3.
За характерный параметр при определении числа Рейнольдса принимали эквивалентный диаметр перепускного окна
Рис. 3. Зависимость коэффициента расхода перепускного окна а от числа Рейнольдса
IX =
4/Р'
/7„„
(9)
где /7рег — периметр перепускного окна.
Коэффициент расширения может быть найден из уравнения
е = 1
С р~рь
(10)
к р
где С — экспериментальный коэффициент.
Так как относительная разность давлений в парной полости и в камере всасывания на фазе регулирования невелика, принимаем кр равным 1.
На рис. 4 и 5 представлены зависимости давления пара хладагента в парной полости от ее объема в процессе регулирования при объемной производительности компрессора 50 %.
Сжатие после фазы регулирования начинается при давлении, превышающем давление всасывания. Эта разность давлений обусловлена наличием газодинамических сопротивлений при прохождении пара из парной полости в камеру всасывания. Увеличение индикаторной работы компрессора определяется площадью индикаторной диаграммы под кривой процесса на участке регулирования производительности. На рис. 5 давление в начале сжатия ниже, чем на рис. 4, а значит, меньше и потери работы на выталкивание пара в камеру всасывания. Таким образом, рисунки наглядно демонстрируют энергетический эффект от использования золотника с измененными кромками.
Как показали расчеты, проведенные с помощью математической модели для компрессора ВХ-130, работающе-
Р, бар \
ь\
300
400
\Л/, см
Рис. 4. Зависимость давления в парной полости р от объема (V в ходе регулирования производительности ст андартным золотником: сплошная линия — расчет; пунктир — жеперимент
Рис. 5. Зависимость давления я парной полости р от объема IV в ходе регулирования прои тодительности ишененным плотником: сплошная линия — расчет; пунктир — эксперимент
го на 422 в режимах парокомпрессорных холодильных машин, использование золотника с кромками, расположен ными под углами, равными углам наклона винтовых ли-
ний соответствующих винтов, позволяет повысить КПД компрессора при регулировании производительности на
2 — 8 % в зависимости от режима. При этом наибольший эффект наблюдается при работе на режимах с невысокими значениями внешней степени повышения давления. Расчетные зависимости внутреннего изоэнтропного КПД от внешней степени повышения давления качественно и количественно согласуются с результатами экспериментальных исследований |3|, что говорит об адекватности разработанной математической модели.
Список литературы
1. А.с. 1691558СССР, МКИ Г 04 С 18/16. Винтовой компрессор / АН. Носков, А.П. Алексеев, В.И. Пекарев, В.И. Ведайко. 1991. Бюл. К» 42.
2. Пекарев В.И.. Ануфриев А.В. Исследование процесса сжатия холодильного винтового компрессора с регулируемой производительностью //Турбины и компрессоры. 2004. № 3, 4 (28, 29).
3. Повышение эффективности работы холодильного пинтового компрессора при регулировании производительности. - Л.: ЛТИХП. Отчет о ПИР по теме 948, № гос. регистрации 01900049636.
4. Холодильные машины: Учебник / Под общ. ред. Л.С. Тимофее веко га. — СПб.: Политехника, 1997.
