CC BY
77
Электронное научно-техническое издание
НАУКА и ОБРАЗОВАНИЕ
Эя Н?ФС 77 - 30569. Государственная регистрация №0421100025. ISSN 1994-0408_
Математическое моделирование рабочих процессов винтового маслозаполненного компрессора
77-30569/260341
# 01, январь 2012 Филиппов И. В.
УДК 621.514
Калужский филиал МГТУ им. Н.Э. Баумана
fiv61@yandex.ru
Значительную часть из общего количества винтовых компрессоров (ВК), используемых в различных отраслях промышленности, составляют винтовые маслозаполненные компрессора (ВКМ). В 2005 году на долю ВКМ приходилось около 85% из всего объёма мирового производства ВК. В России ВКМ изготавливаются на базе унифицированного типоразмерного ряда, который охватывает диапазон по производительности от 4 до 50 м /мин при давлении внутреннего сжатия 0,5...0,8 МПа. Зарубежные фирмы, например, «Atlas Copco» производит ВКМ серии GA производительностью до 90 м /мин и давлением внутреннего сжатия до 1,3 МПа
Отсутствие надежных методик по определению параметров рабочего процесса ВК, а также проблемы, возникающие при эксплуатации ВК в анормальных режимах, послужили причиной проведения работы, результаты которой изложены в данной статье. Целью работы является формирование математической модели ВКМ, которая, в дальнейшем, используется в качестве основы предлагаемых методик для определения термодинамических параметров рабочих процессов компрессора.
ВКМ относятся к классу объёмных компрессорных машин, рабочим веществом которых является масловоздушная смесь. Впрыскивание масла в полости сжатия происходит после отсоединения их от камеры всасывания, т.е. в начале фазы внутреннего сжатия. Основная функция впрыскиваемого масла - это отвод теплоты сжатия и улучшение термодинамических показателей работы компрессора. Оптимальное массовое соотношение «масло-воздух» находится в пределах от 6 до 8. Изменение этого соотношения в ту или иную стороны оказывает существенное влияние
на технико-экономические показатели ВКМ, а как показали проведенные исследования [3], может стать одной из причин аварии.
При впрыскивании масла в полости сжатия часть подводимой к компрессору энергии тратится на разбрызгивание масла и на сообщение ему скорости направленного движения. Поэтому удельная мощность сжатия в ВКМ может быть как меньше, так и больше удельной мощности сухого сжатия в зависимости от условий работы и количества впрыснутого масла [1 ].
Как известно из теории ВК [1], рабочий процесс ВК состоит из трёх фаз: всасывания (линия 0-1), внутреннего сжатия (линия 1-А) и нагнетания (линии А-НА, А-В-НВ и А-С-НС) (рис.1).
Давлениер НВ
НА НС
Л в
А
\с
Объём полостей сжатия № Рис.1. рЖ-диаграмма рабочего процесса ВК.
Впрыскивание масла в полости сжатия происходит после их отсоединения от камеры всасывания, т.е. масло поступает в закрытую термодинамическую систему. Для таких систем справедлив первый закон термодинамики [2, 4]:
аи = dQv + dQd - сЛЬ, (1)
где аи - изменение внутренней энергии закрытой термодинамической системы; dQv -тепло, подводимое (+) или отводимое (-) от системы во внешнюю среду; dQd - тепло, выделяющееся в результате преобразования диссипативных потерь; dL - работа, совершаемая системой при изменении её объёма.
В отличие от винтовых компрессоров сухого сжатия (в дальнейшем ВКС), для которых можно пренебречь dQv, рассмотрение процесса внутреннего сжатия в ВКМ, без учета теплообмена с окружающей средой, может привести к ошибочным результатам [4]. Наличие в сжимаемом рабочем теле значительного количества мелкодисперсного вещества с большой теплоёмкостью (каковым и является масло) несомненно, сказывается на результатах процесса сжатия, поэтому следует принимать dQv ^ 0. Одним из принимаемых допущений при исследовании фазы внутреннего сжатия является то, что впрыснутое масло (в общем случае - жидкость) рассматривается как «внешняя среда», в которую отводится тепло внутреннего сжатия dQv. Влияние теплообмена рабочего тела со стенками камеры сжатия по сравнению с влиянием теплообмена между рабочим телом и жидкостью на рабочий процесс компрессора пренебрежимо мало и, в дальнейшем, оно не рассматривается и не учитывается. Поэтому можно записать следующее выражение для отводимой теплоты:
^ = Мтст ЪТт (2)
где Мт - масса впрыснутого в полости сжатия масла, кг; с т - удельная теплоёмкость масла, Дж/кг*К; ЪТт - повышение температуры масла в результате теплообмена со сжимаемым воздухом, ЪТт = 15...25 К [1].
В качестве идеального рабочего процесса внутреннего сжатия в ВКМ принимается изоэнтропное сжатие, а отклонение от фактического процесса от идеального обуславливается наличием диссипативных потерь и теплообмена. По аналогии с ВКС [4], диссипативные потери, имеющие место в ВКМ, можно разделить на функционально зависимые и независимые от параметров воздуха в камере всасывания. Потери первой группы включают газодинамические потери на трение винтов о сжимаемый газ, газа о корпус и т.п. К независимым потерям следует отнести потери на разбрызгивание и на сообщение жидкости скорости направленного движения, а также потери на трение по плоскости «гребень винта - масляная плёнка -корпус» [4]. Тепло, получающееся в результате преобразования функционально зависящих потерь обозначим dQdу, а от преобразования потерь второй группы - dQd с.
Таким образом, уравнение (1) можно представить в виде:
Мvdu = МтОтЪТт + + ^.с + ^ , (3)
или в конечной форме:
Мг(сг.АТА - с, Т) _ МтатЪТт + а, + & с + (4)
где Му - масса воздуха, заключенного в полостях сжатия, кг; с у л и с у. А - удельная теплоемкость воздуха в начале и после внутреннего сжатия, Дж/кг*К; Т1 и ТА -температура воздуха в начале и после внутреннего сжатия, К; - работа изоэнтропного сжатия, Дж. Работа изоэнтропного сжатия, совершаемая системой при изменении её объёма, определяется из выражения [4]:
Г к-1 Л
Ь _ —^ рЩ г к -1 1 1
П 2г - 1
V У
(5)
где к - показатель адиабаты сжимаемого газа; р 1 - давление воздуха в начале фазы внутреннего сжатия, Па, Пд.г - изоэнтропная степень повышения давления (определяется из паспортных данных ВКМ); Ж1 - объём полостей сжатия в начале фазы внутреннего сжатия, м2.
Выражая тепло dQd.у и dQd.с. как некоторую часть этой работы с использованием коэффициентов диссипативных потерь kd.у и kd.с, соответственно, и выполнив преобразования, получим выражение для вычисления температуры воздуха после внутреннего сжатия:
Р Ж
(1 + ка.у)КР1 + ка. сКР1. пот + ^ МтТт. 1 + с„. 1 Т __'
А~ РЩ л* '
су. + стМ т
(6)
где Тт. 1 - начальная температура впрыскиваемого масла, К; р 1пот - давление воздуха в полостях сжатия после окончания фазы всасывания при работе в номинальном режиме; К - коэффициент, зависящий от паспортных данных компрессора:
Г к-1 Л
К
к -1 1
п А; - 1
V У
(7)
Остальные параметры фазы внутреннего сжатия определяются из выражений: средний показатель политропы т:
1п
Г Т Л
Т А
Т
V У
т _—+ 1
1п е * , (8)
давление воздуха р А:
т
Рл=р [ % ]т-1. (9)
степень повышения давления пА:
п - Ра
п А --
р1 (10)
где 8г - геометрическая степень сжатия ВКМ.
Воспользовавшись перечисленными выражениями, можно получить достаточно точные численные значения всех основных параметров фазы внутреннего сжатия в ВКМ, что на практике означает возможность устанавливать, а при необходимости, и прогнозировать режим работы компрессора при изменении условий эксплуатации. Уравнения (6) - (10) представляет собой методику определения параметров фазы внутреннего сжатия, а комплексное решение названных уравнений позволяет получить математическую модель термодинамических процессов, происходящих во время рассматриваемой фазы.
Для проверки значимости расчётных данных, получаемых с помощью предлагаемой методики, были проведены испытания компрессоров типа 6ВКМ-25/8. Некоторые результаты экспериментальных исследований представлены в таблице 1. Для сравнения, в этой же таблице, приведены расчётные значения параметров внутреннего сжатия, полученные по предлагаемой методике. При проведении испытаний, давление воздуха в камере всасывания р 1 изменялось дросселированием во всасывающем трубопроводе, температура воздуха Т1 - определялась температурой атмосферного воздуха, а начальная температура впрыскиваемого масла Тт .1 -изменялась в зависимости от условий охлаждения масла. Таким образом, была получена достаточно большая область варьирования значений исходных величин, что позволило повысить объективность результатов как экспериментальных, так и теоретических исследований.
В качестве дополнительного замечания к таблице 1, следует прояснить вопрос о том, как же устанавливался основной режим работы компрессора, т.е. какое давление необходимо было поддерживать в нагнетательном трубопроводе, для того чтобы при изменениир 1, Т1 или Тт1 сохранялся режим работы «А». Как показали все дальнейшие исследования, ответ на этот вопрос и является, по сути, решением проблемы повышения эффективности работы винтовых компрессоров.
Таблица 1.
Параметры фазы внутреннего сжатия ВКМ типа 6ВКМ-25/8
Исходные данные Расчётные данные Экспериментальные данные
р и Па Т1, К Ти.1, К Та, К р а, МПа Па Та, К р а, МПа Па
92500 281 331 348,8 0,5741 6,21 347 0,5700 6,20
83800 281 331 347,2 0,5178 6,18 345 0,5150 6,15
75200 281 331 345,7 0,4626 6,15 345 0,4590 6,10
71000 281 331 344,9 0,4358 6,14 344 0,4400 6,19
65500 281 331 343,9 0,4008 6,12 342 0,4000 6,11
93800 292 353 370,6 0,5952 6,35 372 0,6000 6,39
90500 292 353 370,0 0,5734 6,34 370 0,5700 6,30
83200 292 353 368,7 0,5253 6,31 370 0,5150 6,19
77700 292 353 367,8 0,4893 6,29 368 0,4900 6,31
70600 292 353 366,5 0,4431 6,27 365 0,4400 6,23
Фаза нагнетания в ВКМ следует непосредственно за внутренним сжатием и в зависимости от соотношения давлений: в полостях сжатия -рА и в камере нагнетания -р н , компрессор будет работать в одном из режимов: при р А = р НА - в основном режиме «А», при ра < рНВ или ра > рНС - в анормальных режимах, соответственно, дожатия «В» или расширения «С» (рис. 1) [1, 4]. Процесс выравнивания давлений рА и рН происходит практически мгновенно, вследствие чего, изменением геометрического объёма полостей сжатия во время этого процесса, а также теплообменом между рабочим телом, корпусом и винтами можно пренебречь. Так и было сделано при исследовании ВКС: процесс выравнивания давлений во время фазы нагнетания рассматривался как адиабатное дожатие (режим «В») и адиабатное расширение (режим «С») [4]. Однако, при изучении фазы нагнетания в ВКМ, нельзя не учитывать теплообмен между воздухом и впрыснутым маслом и, по всей видимости, процессы дожатия или расширения следует рассматривать как политропные.
При работе компрессора в режиме «С» температура расширяющегося воздуха уменьшается: температурное равновесие между воздухом и впрыснутым маслом нарушается. Следовательно, можно предположить, что в данном случае, теоретически
возможен некоторый возврат тепла от масла к воздуху, несмотря на значительную разницу их теплоёмкостей. Последнее замечание позволяет пренебречь теплообменом в рассмотренном случае и не учитывать влияние масла на процесс выравнивания давлений. Дальнейшие теоретические и экспериментальные исследования подтвердили обоснованность принятого допущения. Поэтому, как и для ВКС, процесс выравнивания давлений во время фазы нагнетания в ВКМ для режима «С» рассматривается как необратимое адиабатное расширение из полости постоянного объёма [4]. Дифференциальное уравнение первого закона термодинамики для такого процесса запишется в виде:
Mdu + иМА-Нdx - hA.mMA_Нdx = 0, (П)
где Ми u - текущие значения, соответственно, массы и удельной энергии воздуха в полостях сжатия, кг и Дж/кг; МА-Н - массовый расход воздуха, поступающего из полостей сжатия в камеру нагнетания, кг/с; h Ат - удельная энтальпия торможения воздуха поступающего из полостей сжатия, Дж/кг^К; т - время, с.
Учитывая, что МА-Н dт = - dM, а также h - u = pW = RT, можно уравнение (11) переписать в виде:
cvMdT - RTdM = 0, (12)
Разделяя переменные в уравнении (12) и интегрируя полученное дифференциальное уравнение, имеем:
М ^
— = Const,
Т (13)
а с учетом М = рW/RT и пренебрегая незначительным изменением объема полостей сжатия за время выравнивания давлений
= Const
T (14)
Состояние равновесия, наступившее после выравнивания давлений, характеризуется температурой нагнетания ТНС и давлением р НС. Следовательно, можно записать
к-1 к-1 РАк _ РнС
T T
А ±НС (15)
или
к-1
Т = Т
1 нс 1 а
ГРнс_ V к Ра
Как видно из последнего выражений, температура нагнетаемого воздуха при работе ВКМ в режиме «С» зависит от параметров воздуха после внутреннего сжатия и в камере нагнетания.
При работе компрессора в режиме «В» температура нагнетаемого воздуха вследствие внегеометрического дожатия повышается. Исследование выравнивания давлений как адиабатного процесса для ВКМ может привести к ошибочным результатам. Наличие в системе мелкодисперсного вещества с большой теплоёмкостью может оказать существенное влияние на конечные результаты процесса выравнивания давлений. Площадь контактного теплообмена воздуха с маслом для ВКМ значительна: так, например, если принять, что всё впрыснутое масло распыляется на капли шарообразной формы, то суммарная площадь капель при диаметре в 1 мм для компрессоров типа 6ВКМ-25/8 составит 0,1217 м , а при диаметре капель 0,5 мм -0,2435 м . Поэтому для ВКМ, работающего в режиме «В», необходимо учитывать теплообмен между дожимаемым воздухом и маслом, а процессы, происходящие во время фазы нагнетания, следует рассматривать как политропное сжатие.
Рассмотрим открытую термодинамическую систему, рабочим веществом которой является дисперсная система: масло - дисперсная фаза, воздух - дисперсная среда. Первый закон термодинамики для открытой термодинамической системы с учётом «окружающей» среды (в виде мелкодисперсного масла) можно записать в виде:
Mdu + иМН-Аdт - кНтМН-Аdт - -dQv.
(17)
где dQv - тепло, переданное маслу от дожимаемого воздуха, Дж; МН- А - массовый расход воздуха, поступающего из камеры нагнетания в полости сжатия, кг/с; к н.т -удельная энтальпия торможения воздуха поступающего из камеры нагнетания в полости сжатия, Дж/кг^К.
Теплота dQv может быть определена из выражения (2), следовательно:
сМТ + е^М - СрТн^М + МтСтЪТт - 0. (18)
где ТНл - температура адиабатного торможения поступающего из камеры нагнетания воздуха, К.
После преобразования данного выражения, разделяя переменные, получим линейное дифференциальное уравнение первого порядка:
сСТ сСМ
cpTH ,t - cvT cvM + cmMm (19)
или, разделив знаменатель и левой, и правой частей на сv и принимая сm/cv = Кmv, получаем:
dT _ dM kTH ,t - T _ M + KmvMm ' (20)
Принимая температуру торможения потока ТН. t постоянной и выполнив интегрирование, получим:
(kTH, t - Т )(KmMm + М ) _ Const, (21)
Записав выражение (21) для начального и конечного состояния воздуха в течение фазы нагнетания («начало» - р А и ТА; «конец» - р НВ и Тнв), принимая М = p/RT, получим следующее уравнение:
f \ / Л
(kTH . t - Та ) Kmv М m _ (kTH . t - Тнв )
v rta j
KM- + Рнв
mv m
V RThb J
(22)
В результате преобразования данного выражения, получим квадратичное уравнение типа:
АТНв + ВТНВ + С = (23)
где А, В и С - коэффициенты квадратичного уравнения:
А = К^ МmR,
B = ^Н.РАРа - АТА + РНВ - РА, С = . г рнв . (24)
Решая уравнение (23) (с учётом (24)) относительно ТНВ и принимая во внимание только положительные корни этого уравнения, можно получить искомые значения температуры нагнетаемого воздуха при работе ВКМ в режиме «В».
Последовательное решение уравнений (6) и (16), или (6) и (23), позволяет получить значение температуры нагнетаемого воздуха при работе ВКМ в любом из режимов, а значит, получить её зависимость от условий работы. В качестве примера на рисунке 2 приведены расчётная и экспериментальная зависимости температуры нагнетаемого воздуха для различных режимов работы компрессора 6ВКМ-25/8.
Рис. 2. Расчетные и экспериментальные зависимости температуры Тн от давления
воздуха ркн:
♦ и А - при р 1 = 93592 Па; Т1 = 288 К; ТтЛ = 328 К р = 563300 Па; Та = 346 К); • и о - при р = 93791 Па; Т1 = 267 К; ТшЛ = 307 К (р А = 571600 Па; ТА = 325 К)
Используя полученную математическую модель можно строить 3Б зависимости температуры нагнетаемого воздуха от других параметров процесса. Например, на рисунке 3 представлена 3Б зависимость температуры нагнетаемого воздуха в режиме «В» от давления в камере нагнетания и температуры воздуха после внутреннего сжатия.
Рис. 3. 3Б зависимость Тнв от рнв и ТА.
Подводя итог, можно сделать следующие выводы:
предлагаемая математическая модель ВКМ представляет собой многофакторную функцию переменных величин, являющихся термодинамическими параметрами рабочего процесса;
математическая модель ВКМ позволила разработать методику определения параметров воздуха после внутреннего сжатия и в конце фазы нагнетания для ВКМ;
сравнивая представленные на рисунке 2 зависимости, можно отметить хорошую согласованность расчетных и экспериментальных данных, а это означает, что предлагаемая методика позволяет достаточно точно прогнозировать режим работы компрессора и рассчитывать соответствующие ему основные параметры воздуха.
Список литературы
1. Сакун И.А. Винтовые компрессоры. Л.: Машиностроение, 1970. 400 с.
2. Цейтлин Ю.А., Мурзин В.А. Пневматические установки шахт. М.: Недра, 1985. 352 с.
3. Филиппов И.В. Ударные волны в винтовых маслозаполненных компрессорах.// Изв. Вузов. Машиностроение. 1995, №1-3, С. 26-32.
4. Филиппов И.В. Работа винтовых компрессоров при переменных параметрах газа во всасывающем и нагнетательном трубопроводах.// Изв. Вузов. Горный журнал. 1989, №2, С.109-112.
electronic scientific and technical periodical
SCIENCE and EDUCATION
_EL № KS 77 - 3Ü56'». .V;II421100025, ISSN 1994-jMOg_
Mathematical modeling of screw oil-filled compressor's working processes
77-30569/260341
# 01, January 2012 Filippov I.V.
Bauman Moscow Technical University, Kaluga Branch
fiv61@yandex.ru
Nowadays there is no method reliable enough to calculate the set of important screw oil-filled compressor parameters analitically, in particular, parameters of inner compression and injection. The results of the survey, included in the article, have allowed creating a mathematical model of different processes, taking place in screw oil-filled compressor, and developing the methods for the calculating of inner compression and injection parameters. The comparative analysis of theoretical and experimental results showed good correlation between each other, that means this method can be recommended for the forecasting of the compressor working mode and calculating the basic parameters of this mode.
Publications with keywords: mathematical model, oil-injected rotary screw compressor, operating conditions Publications with words: mathematical model, oil-injected rotary screw compressor, operating conditions
Reference
l.Sakun I.A., Screw compressors, Lrningrad (L.), Mashinostroenie, 1970, 400 p.
2. Tseitlin Iu.A., Murzin V.A., Pneumatic installation of mines, Moscow, Nedra, 1985, 352 p.
3. Filippov I.V., Shock waves in the oil-injected screw compressors, Izv. Vuzov. Ser. Mashinostroenie 1-3 (1995) 26-32.
4. Filippov I.V., Work screw compressors at variable parameters of the gas in the suction and discharge piping, Izv. Vuzov. Ser. Gornyi zhurnal 2 (1989) 109-112.
