автореферат диссертации по энергетическому, металлургическому и химическому машиностроению, 05.04.03, диссертация на тему:Особенности процесса сжатия холодильного винтового маслозаполненного компрессора с впрыском жидкого рабочего вещества

ДИССЕРТАЦИЯ.

Введение:.

1 Моделирование винтовых маслозаполненных компрессоров.

1.1 Расчет винтовых компрессоров. Теория стационарного теплообмена.

1.2 Расчет винтовых компрессоров. Нестационарный теплообмен.

1.3 Процесс сжатия поршневого компрессора с впрыском воды в цилиндр.

1.4 Впрыск воды в винтовой компрессор сухого сжатия.

1.5 Инженерный метод расчета процесса сжатия с впрыском жидкости.

1.6 Зарубежные исследования процесса сжатия с впрыском жидкости.

2 Математическая модель винтового маслозаполненного компрессора.

2.1 Математическое моделирование.

2.2 Математическое описание процесса сжатия винтового компрессора.

2.3 Алгоритм расчета процесса сжатия винтового маслозаполненного компрессора.

2.4 Компьютерная программа расчета процесса сжатия.

3 Математическая модель винтового маслозаполненного компрессора с впрыском жидкого рабочего вещества в рабочую полость.

3.1 Допущения и характер процесса.

3.2 Формулы и компьютерная программа.

3.3 Алгоритм расчета процесса сжатия винтового маслозаполненного компрессора с впрыском жидкого хладагента.

3.4 Теоретический эксперимент.

4 Теоретические эксперименты с математическими моделями.;.

4.1 Работа сжатия.

4.2 Оптимальная точка впрыска.

Оглавление

4.3 Рабочие вещества.

4.4 Количество и размер капель жидкости.

Выводы:.

Актуальность работы. Винтовые компрессоры часто используются в холодильных системах большой мощности (30.450 кВт). Они отличаются низкой пульсацией скручивающего момента и отсутствием необходимости во всасывающем и нагнетательном клапанах, что делает их легкими, прочными, к тому же у них низкая вибрация.

Экономичность работы винтового компрессора, как компрессора объемного действия, во многом определяется термодинамической эффективностью процесса сжатия, во время которого расходуется основная часть внешней работы, подводимой к компрессору. Из этого следует, что одним из основных путей повышения термодинамической эффективности процесса сжатия будет улучшение охлаждения сжимаемого газа за счет отвода тепла к дополнительным поверхностям теплообмена (при достаточно высоком коэффициенте теплообмена и значительной разнице температур между охлаждаемым газом и охлаждающим телом) или за счет испарительного охлаждения. Это достигается путем впрыска жидкости в сжимаемый газ. Впрыск жидкости, кроме уменьшения работы сжатия, также снижает тепловую нагрузку на детали и уплотняет щелевые зазоры, уменьшая перетечки.

В настоящее время известен ряд работ, посвященных термодинамике процесса сжатия объемного компрессора с двухфазным рабочим телом. К недостаткам этих работ следует отнести либо невысокую точность, полученных результатов, либо значительную сложность проводимых расчетов (дающих, кстати, не всегда верный результат). Так же необходимо отметить, что все относительно современные и интересные (с точки зрения понимания процесса) исследования, проводившиеся с винтовыми компрессорами, проводились в основном с компрессорами сухого сжатия, сжимающими воздух, где в качестве впрыскиваемой жидкости использовалась вода. Исследования винтовых маслозаполненных компрессоров основывались в большинстве на несовершенной теории стационарного процесса теплопередачи. Исследований возможности охлаждения винтовых компрессоров впрыском жидкого хладагента в рабочую полость в полной мере в нашей стране не проводилось вовсе. Такой способ охлаждения сжимаемого газа дает возможность получить схему холодильной машины, работающей без маслоотделителя. От впрыска жидкого хладагента в процессе сжатия, температура нагнетания понижается настолько, что становится возможным масло, после маслоотделителя, подавать непосредственно в компрессор, без дополнительного охлаждения.

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

Цель работы; накопление и обобщение научных материалов, посвященных впрыску жидкости в полость сжатия винтового компрессора, составление математической модели термогазодинамических процессов в полости винтового компрессора и разработка рекомендаций по расчету и проектированию винтовых компрессоров с впрыском жидкости.

В соответствии с этой целью необходимо было:

• провести обработку и обобщение ранее известных экспериментальных исследований.

• разработать математическую модель термогазодинамических процессов в полости маслозаполненного винтового компрессора.

• провести серию теоретических экспериментов с данной моделью для ее проверки по результатам опубликованных экспериментов.

• усовершенствовать математическую модель ВМК для расчета процесса сжатия ВМК с впрыском жидкого хладагента.

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

• разработать методику по расчету и проектированию ВМК с впрыском жидкости, составить рекомендации по использованию разработанных компьютерных программ для расчета внутренних процессов винтовых компрессоров.

Введение

Научная новизна работы: заключается в следующем - собрана большая база данных по вопросу впрыска жидкости в компрессор, которая позволяет всесторонне рассмотреть процесс впрыска жидкости в полость сжатия.

• разработана математическая модель и компьютерная программа впрыска жидкого хладагента в полость сжатия ВМК.

• проведен расчет холодильной машины работающей без маслоохладителя с впрыском жидкого хладагента.

• созданы компьютерные программы, позволяющие использовать их для расчетов внутренних процессов при разработке и моделировании ВМК.

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

• методику моделирования, позволяющую моделировать внутренние процессы винтового компрессора, принимаемые допущения течения физического процесса.

• математическую модель винтового компрессора с впрыском масла и жидкого хладагента.

• математическую модель винтового маслозаполненного компрессора работающего без маслоохладителя.

• методику моделирования, позволяющую уверенно моделировать все внутренние процессы винтового компрессора.

Апробация работы: основные положения диссертации докладывались и обсуждались на Международной научно-технической конференции "Холод и пищевые производства", Санкт-Петербург, 1997 г.

Публикации: по материалам работы опубликовано три работы

Структура и объем работы: Диссертация состоит из введения, четырех глав, выводов, списка литературы, включающего 47 наименований, и трех приложений. Работа изложена на 131 машинописных страницах, содержит 19 рисунков, 2 таблицы, 15 страниц приложений.

выводы:

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

При впрыске масла в рабочую полость винтового маслозаполненного компрессора эффект от охлаждения сжимаемого газа проявляется в увеличении энергетического КПД на 3.5%.

Установлено, что при впрыске жидкого хладагента в полость сжатия в обычных режимах эксплуатации среднетемпературной холодильной машины улучшения энергетических характеристик не наблюдается, т.к. выигрыш от испарительного охлаждения сжимаемого газа полностью перекрывается затратами мощности на сжатие образовавшегося пара. Впрыск испаряемой жидкости становится энергетически выгоден, только при условии образования капель хладагента диаметром меньше 25.30 мкм. Проведенный теоретический эксперимент показал возможность создания схемы холодильной машины работающей без маслоохладителя. Увеличение потребляемой мощности составляет при этом не более 5.6%. Эксплуатация такой машины может оказаться эффективной при определенных условиях, например при дефиците охлаждающей воды.

При использовании впрыска жидкого хладагента в рабочую полость винтового компрессора из-за образования дополнительного пара давление в полости возрастает быстрее и поэтому необходимо учесть это при расчете площади и момента открытия окна нагнетания. Созданные модели позволяют достаточно точно провести этот расчет.

Впрыск жидкого хладагента позволит снизить температуру рабочего вещества в конце процесса сжатия на 12. 16%.

Вследствие уменьшения конечной температуры сжатия теплонапряженность узлов и деталей заметно снижается.

Место впрыска жидкого хладагента в полость сжатия влияет на энергетические показатели винтового компрессора, расчет показал, что отверстие для впрыска целесообразно разместить ближе к окну нагнетания.

Хладагенты с высокой теплотой парообразования более эффективны для винтового компрессора работающего без маслоохладителя, т.к. эффективное испарительное охлаждение позволяет сократить потери от сжатия образовавшегося пара.

1. Алешин В. И. "Исследование винтового маслозаполненного вакуум-компрессора" Автореф. дис. на степень канд. техн. наук. - М., 1977. - 24 с.

2. Амелин А. Г. "Теоретические основы образования тумана при конденсации пара" М. Химия - 1972г.

3. Берман Я. А., Булыгин В. Г., Рафалович А. П., Жилкин А. Н. "Охлаждение компрессора ТВ-80-1,6 впрыском жидкости в проточную часть" Хим. и нефт. машиностроение - 1984 - №{12} - с.27-28

4. Богданов С. Н., Иванов О. П., Куприянова А. В. "Холодильная техника. Свойства веществ: Справочник" М. - Агропромиздат - 1985г.

5. Верный А. Л., "Исследование и метод расчета винтовых маслозаполненных компрессоров" В сб. Тр. ВНИИкриогенмаш: Процессы, технология и контроль в криогенном машиностроении. - Балашиха - 1978г. - с.72-82

6. Горбис В. Р. "Теплообмен дисперсных сквозных потоков" М. - JI.Энергия -1964г. - 296с.

7. J. Грин X., Лейн В. "Аэрозоли, пыли, дымы и туманы" Л.Химия - 1969г. - 428с.

8. S. Кабаков А. Н., Щерба В. Е. "Математическое моделирование рабочего процесса с впрыскиванием воды в поток сжимаемого воздуха" Изв. вузов - Горный журнал - 1981г. -№{1} -с.71-73

9. Калугин Г. Н. "Винтовые компрессоры с подачей жидкости в рабочую полость" Краснодар - 1994г.

10. Калугин Г. Н., Амелин В. И., Мишин Ю. П. "Влияние впрыскиваемойжидкости на уплотнение щелей в винтовых машинах" Холодильные машины- 1978г. -№{Ю}

11. Канышев Г. А., Чистяков Ф. М. "Влияние свойств масел на энергетические характеристики фреоновых маслозаполненных винтовых компрессоров" -Холодильная техника 1980г. - №{7}

12. Канышев Г. А., Чистяков Ф. М. "Коэффициент подачи винтового фреонового маслозаполненного компрессора" Холодильные машины - 1979г. - №{12}

13. Карслоу Г. Е., Еэгер Д. "Теплопроводность твердых тел" М. Наука - 1964г.

14. Кафаров В. В. "Основы массопередачи" М. - Высш. Ш - 1972г.

15. Коренысов В. И., Немировский С. К. "Термодинамический расчет винтового маслозаполненного компрессора" сб. Тр. СОАН - 1976г.

16. Кошкин Н. Н., Сакун И. А., Бамбушек Е. М. "Холодильные машины" Л. Машиностроение - 1985г.

17. Кутателадзе С. С., Стрикович М. А., "Газодинамика гидрожидкостных систем"- М. Энергия 1976г. - 296с.

18. Ландау Л. Д., Лившиц Е. М. "Механика сплошных сред" М. Гостехиздат -1954г.

19. Левенцов А. А. "Методика расчета процесса сжатия винтового маслозаполненного компрессора с впрыском жидкости в полость сжатия" Сб. тезисов к научной конференции: Проблемы тепломассопереноса - СПбГАХПТ -1995г.

20. Левич В. Г. "Физикохимическая гидродинамика" М. Физматиз - 1959г.

21. Лытевский А. С. "Движение жидких капель в газовом потоке" Изв. вузов СССР - Энерг. - 1963г. - №{7} - с.75-81

22. Назаров О. И., Поваров О. И., Ятцени И. А. "Удар капли о плоскую движущуюся пластину" Изв. вузов СССР - Энерг. - 1984г. - №{8} - с.53-56

23. Пекарев В. И. "Исследование работы винтового компрессора в режимах паровых холодильных машин" Диссертация на соиск. ученой степени канд. техн. наук - 1969г.

24. Перелыптейн И. И., Парушин Е. Б. "Термодинамические и теплофизические свойства рабочих веществ холодильных машин и тепловых насосов" М. -Легкая и пищевая промышленность - 1984г.

25. Плотников В. Т., Филаткин В. И. "Разделительные вымораживающие установки" Агропромиздат - 1987г.

26. Рамм В. М. "Абсорбция газов" М. - Химия - 1976г.

27. Ю.Ребриков В. Д., Фотин Б. С., Хрусталев Б. С., Сидора Н. Н. и д.р. "Влияние впрыска жидкости на рабочий процесс объемного компрессора" В сб. Тр. ЦКТИ: Винтовые компрессоры в энергомашиностроении. - Л. - 1975г. - Вып. 127 - с.82-88

28. Я.Сакун И. А. "Винтовые компрессоры" Л. Машиностроение - 1970г.

29. Слободянюк Л. И., Гогин Ю. Н. "Охлаждение компрессора впрыском воды в цилиндр" Изв. вузов СССР - Энерг. - 1961г. - №{9} - с.62-663."'Теоретические основы тепло- и хладотехники" под ред. проф. Э. И. Гуйго -Изд-во Ленингр. ун-та 1976г.

30. Фукс Н. А. "Испарение и рост капель в газообразной среде" М. Издательство АН СССР - 1958г. - 90с.

31. Шервуд Т. К., Рид Р. "Свойства газов и жидкостей" Л. - 1982г.

32. Шнейдер П. Д. "Теплопроводность. Форсуночные аппараты в химической технологии" М. Ин. лит. - Инж. пробл. технологии - 1960г.

33. Щерба В. Е., Кабаков А. Н., Юша В. Л. "Эффективность рабочего процесса сжатия компрессора объемного действия с впрыском испаряющейся жидкости"

34. Юша В. Л., Щерба В. Е., Кабаков А. Н. "Анализ рабочего процесса всасывания винтового компрессора с впрыском жидкости" Изв. вузов. Энерг. - 1985г. -№{9} - с.81-85

35. Cramer E. U. "Die Einspritztachnik bei Schranubenverdichtern" Chem Prod. -1978 -7 -№{9} -34-36-38

36. La Monica G., Enrile A. "Theoretical and experemintal results on an oil and liquid refrigerant injected screw compressor" Mfilanges fluid, frigerigenes: propr. et apple. -Paris- 1980 - 179-185

37. Persson J. G. "Heat-exchnge in liquid-injected screw-compressors" VDI-Ber. -1987-№{640} - 121-135

38. Shimoji Mihoko, Koda Toshihide "Simulation Technology for a liquid-injection Single-Screw Compressor" Mitsubushi Elec. Adv. - 1995 - 70, March - 25-27

39. Si dem an S., Isenberg J. Direct Contact Heat Transfer with change of Phase: Evaporation of Drops Immischible Liquid Medium. Int. Journal Heat and Mass Transfer 1964 - №{7} - p. 1273-1289

40. Wasser ersetrt ul olfreie Schraubenkompressoren mit Wassereinspritzung // Production 1989-№{23}1. Триложения:

41. Программа расчета процесса сжатия винтового маслозаполненногокомпрессора на языке Фортран-4

42. ПРОГРАММА РАСЧЕТА ПРОЦЕССА СЖАТИЯ МАСЛ03АП0ЛНЕН0Г0 ВК

43. РАБ. ВЕЩЕСТВО АММИАК ,МАСЛО ХС 4 0 ,Т КИП=-10 ,Т КОНД=35,Т МАСЛА=35.

44. DOUBLE PRECISION RO,TG,CV,TGI,TG2,Р,PI INTEGER I,К

45. DIMENSION ZTG(290),TM(290),VRAM(290),DKM(290) ,ALF1(290) , *PG(2 9 0),DTG(290),DTM(290),CVG(290),QM(290),NKM(290),ZR(290) EXTERNAL SALF,ALI,ATSG1. WRITE(*,*) 'НАЧАЛО'1. RG=4 8 8.16 RAB=0.

46. ZTG(1)=2 63.15 TM(1)=308.15 VKAM(l)=4.68E-4 ZMM=(8.78E-4) VM=ZMM/827. DKM=1.E-4 VKM=5.24E-13 FKM=3.14E-8 CPM=2172. NKM(1)=VM/VKM

47. WRITE(*,*) ' ЧИСЛО КАПЕЛЬ МАСЛА=',NKM(1) DNKM=01. DTAU=0.011/2901. Z DV=(3.85E-6)/2.91. WRITE(*,*) ' Z DV=',Z DV1. CALL R7171. TG1=258.15 CCCC1. TG1=2 63.151. TG2=ZTG(1)1. CALL PHAC (TGI,P)1. CALL PL(TG2,P,RO)

48. WRITE(*,*) ' ПЛОТНОСТЬ ГАЗА НА ВСАСЫВАНИИ=',RO*1000 WRITE(*,*) ' ДАВЛЕНИЕ НА ВСАСЫВАНИИ=',Р*100000 ZMG=VKAM(1)*RO*l.E3

49. WRITE(*,*) ' МАССА НА ВСАСЫВАНИИ=',ZMG DO 1 1=1,290 WRITE(6,11) I 11 FORMAT( 2X,'ИДЕТ РАСЧЕТ ТОЧКИ НОМЕР ',13) RO=(ZMG/VKAM(I))*1.Е-3 CALL SALF(ZTG(I),RO,DKM,ALFl(I)) WRITE(*,*) ' ALFA=',ALF1(I) TG=ZTG(I)

50. CALL CVROT(RO, TG,CV) CVG(I)=CV*10 0 0

51. WRITE(*,*) ' ТЕПЛОЕМКОСТЬ=',CVG(I)

52. CALL PF(RO,TG,Z,PI) PG(I)=P1*100000

53. WRITE(*,*) ' ДАВЛЕНИЕ=' , PG(I)

54. QM(I)=NKM(I)*ALF1(I)*FKM*(TM(I)-ZTG(I))*DTAU

55. DTG(I) = (QM(I)+PG(I)* ZDV)/ (ZMG* CVG(I) ) WRITE(*,*) ' ИЗМЕНЕНИЕ ТЕМПЕРАТУРЫ ГАЗА=',DTG(I) ZTG(I+1)=ZTG(I)+DTG(I)

56. WRITE(*,*) ' ТЕМПЕРАТУРА ГАЗА=',ZTG(I+1) DTM(I)=QM(I)/(ZMM*CPM)

57. WRITE(*,*) ' ИЗМЕН.ТЕМПЕРАТУРЫ МАСЛА=',DTM(I) TM(1+1)=TM(I)-DTM(I)

58. WRITE(*,*) ' ТЕМПЕРАТУРА МАСЛА=',TM(I+1)1. VKAM(1+1)=VKAM(I)-ZDV1. NKM(1+1)=NKM(I)-DNKM1. RAB=RAB+PG(I)*ZDV1. ZR(I)=RAB

59. WRITE(*,*)' РАБОТА=',RAB IF (PG(I).GT.1167500) GO TO 77 CCCC IF (PG(I).GT.1348500) GO TO 77 1 CONTINUE7 7 WRITE(*,*) ' СОЗДАНИЕ ФАЙЛА С РЕЗУЛЬТАТАМИ'

60. OPEN(9, FILE="C:\NAUKA\OTVSTD.DAT",STATUS='OLD') DO 9 K=1,2 90

61. WRITE(9,12) K,ZTG(K),PG(K),ZR(K) 12 FORMAT (2X,13,2X, 'T=',F5.1,2X, 'P=' , E9.3, 2X, 'L=',F10.4) 9 CONTINUE CLOSE (9)1. WRITE(*,*) ' КОНЕЦ!'1. STOP END

62. SUBROUTINE SALF(TG,ZN,D,ALF) DOUBLE PRECISION T,RO,ETA,ZN, PR T=TG RO=ZN1. CALL LMDPAP(T,RO,ETA)1. BETA=0.51. VI=0.8E-61. UO=14.241. RE=(UO*D)/VI1. WRITE(*,*) ' RE=',RE1. CALL PRDTLPAP(T,RO,PR)1. WRITE(*,*) ' PR=',PR

63. ZNU=2 *(1+(BETA*(SQRT(RE) ) *(PR**(1/3) )) )1. WRITE(*,*) ' NU=',ZNU1. ALF=ZNU*ETA/D1. RETURN1. END

64. Программа расчета процесса сжатия винтового маслозаполненного омпрессора, с впрыском жидкого рабочего вещества в полость сжатия, наязыке Фортран-4.

65. ПРОГРАММА РАСЧЕТА ПРОЦЕССА СЖАТИЯ ВК С ВПРЫСКОМ ЖИДКОГО РАБОЧЕГО ВЕЩЕСТВА АММИАКА ,МАСЛО ХС 4 0 ,Т КИП=-15 ,Т КОНД=ЭО

66. DOUBLE PRECISION RO,TG,CV,TGI,TG2,Р,PI,ROJ,Z,TV,Р2,ROP, *EG,EJ,ZES,TS,ROS,EP,PS1,TKP,RP INTEGER I,К

67. DIMENSION ZTG(290),TM(2 90),TX(290),VKAM(290),ZMX(290), *VX(290),DKM(290),VKX(290),ALF1(290),ALF2(290),ZMG(290), *FKX(290),PG(290),DTG(290),DTM(290),CPX(290) ,DKX(290) , * CVG(290),ZR(290),DTX(290),QX(290),QM(290),NKM(290), *NKX(290),ZMISP(290)

68. EXTERNAL ALI,ATSG,SALF,TROPE

69. WRITE(*,*) 'НАЧАЛО' RAB=0.

70. ТЕМПЕРАТУРА ГАЗА НА ВСАСЫВАНИИ

71. ZTG(l)=263.15 ТЕМПЕРАТУРА МАСЛА ПОДАВАЕМОГО В КОМПРЕССОР ТМ(1)=328.15

72. ПОЛНЫЙ ОБЪЕМ КАМЕРЫ ЗАПОЛНЯЕМЫЙ ГАЗОМ (БЕЗ УЧЕТА КОЭФ. ЗАПОЛНЕНИЯ) VKAM(1)=4.68Е-4

73. МАССА ВСЕГО МАСЛА ВПРЫСКИВАЕМОГО В ОДНУ РАБОЧУЮ ПОЛОСТЬ В ПРОЦЕССЕ СЖАТИЯ

74. ZMM=(8.78Е-4) ОБЪЕМ МАСЛА ПРИ НАЧАЛЬНОЙ ТЕМПЕРАТУРЕ

75. VM=(ZMM)/827. ДИАМЕТР КАПЛИ МАСЛА, ПРИНИМАЕТСЯ ЭМПЕРИЧЕСКИ DKM=1.Е-4

76. ОБЪЕМ ОДНОЙ СРЕДНЕСТАТИСТИЧЕСКОЙ КАПЛИ МАСЛА

77. VKM=5.24Е-13 ПЛОЩАДЬ ПОВЕРХНОСТИ ОДНОЙ КАПЛИ МАСЛА FKM=3.14Е-8

78. ТЕМПЕРАТУРА МАСЛА В ПРОЦЕССЕ СЖАТИЯ ИЗМЕНЯЕТСЯ МАЛО (2.3 К) ПОЭТОМУ СЧИТАЕМ ТЕПЛОЕМКОСТЬ МАСЛА ПОСТОЯННОЙ ВЕЛИЧИНОЙ

79. СРМ=2172. ЧИСЛО КАПЕЛЬ МАСЛА NKM(1)=VM/VKM

80. WRITE(*,*) ' ЧИСЛО КАПЕЛЬ МАСЛА=',NKM(1) ПРОЦЕСС СЖАТИЯ В ПОЛОСТИ УСЛОВНО РАЗДЕЛЕН НА 2 90 ШАГОВ (ЭТО ОБЕСПЕЧИВАЕТ НЕОБХОДИМУЮ ТОЧНОСТЬ ВЫЧИСЛЕНИЙ И УДОБНО ТЕМ, ЧТО ОДИН ШАГ СОООТВЕТСТВУЕТ ОДНОМУ ГРАДУСУ ПОВОРОТА ВЕДУЩЕГО РОТОРА). DNKM=0

81. ВРЕМЯ ЗА КОТОРОЕ ВАЛ ВЕДУЩЕГО РОТОРА ПОВОРАЧИВАЕТСЯ НА ОДИН ГРАДУС

82. DTAU=0.011/290 ИЗМЕНЕНИЕ ОБЪЕМА ПОЛОСТИ ЗА ОДИН ШАГ ZDV=(3.85Е-6)/2.9 WRITE(*,*) ' Z DV=',Z DV ВЫЗЫВАЕТСЯ ПОДПРОГРАММА ОПИСЫВАЮЩАЯ ТЕРМОДИНАМИЧЕСКИЕ СВОЙСТВА ХЛАДАГЕНТА CALL R717 TG1=258.15

83. ОТСУТСТВУЕТ ПЕРЕГРЕВ НА ВСАСЫВАНИИ ТЕМПЕРАТУРА ГАЗА НА ВСАСЫВАНИИ РАВНА ТЕМПЕРАТУРЕ КИПЕНИЯ TG1=2 63.15 TG2=ZTG(1) ДАВЛЕНИЕ ВСАСЫВАЕМОГО ГАЗА

84. CALL РНАС(TGI,Р) ПЛОТНОСТЬЬ ВСАСЫВАЕМОГО ГАЗА CALL PL(TG2,Р,RO)

85. WRITE(*, *) ' ПЛОТНОСТЬ ГАЗА НА ВСАСЫВАНИИ=',RO*1ООО WRITE(*, *) ' ДАВЛЕНИЕ НА ВСАСЫВАНИИ=',Р*100000 ZMG(1)=VKAM(1)*RO*1.ЕЗ

86. WRITE(*,*) ' МАССА НА ВСАСЫВАНИИ=',ZMG(1)

87. DO 1 1=1,290 WRITE(6,11) I 11 FORMAT( 2X,'ИДЕТ РАСЧЕТ ТОЧКИ НОМЕР ',13) RO=(ZMG(I)/VKAM(I))*1.E-3 CALL SALF(ZTG(I),RO,DKM,ALF1(I)) WRITE(*,*) ' ALFA=',ALF1(I) TG=ZTG(I)

88. CALL CVROT(RO,TG,CV) CVG(I)=CV*10 0 0

89. WRITE(*,*) ' ТЕПЛОЕМКОСТЬ=',CVG(I)

90. CALL PF(RO,TG,Z,PI) PG(I)=P1*100000

91. WRITE(*,*) ' ДАВЛЕНИЕ=',PG(I)1. (I.LT.90) GOTO 3551. (I.EQ.90) GOTO 5551. (I.GT.90) GOTO 5 65355 QX(I)=0.1. ZMG(1+1)=ZMG(1) GOTO 24

92. HA 90 ШАГЕ ВПРЫСКИВАЕТСЯ ЖИДКОСТЬ КОЛИЧЕСТВО ВПРЫСКИАЕМОГО ХЛАДАГЕНТА РАВНО VPR 555 DKX(90)=8.Е-5

93. CALL ТНАС(PI,ТКР) ТХ(I)=ТКР

94. WRITE(*,*) ' TEMPER XLAD=',ТХ(I) МАССА ВПРЫСКИВАЕМОЙ ЖИДКОСТИ VPR=(8.78Е-4)*1.1 Х=0.13

95. ТОЛЬКО 30 ПРОЦЕНТОВ ЖИДКОСТИ ОСТАЕТСЯ В ПОТОКЕ ГАЗА, ЭТА ЖИДКОСТЬ

96. УЧАСТВУЕТ В ПРОЦЕССЕ ТЕПЛООБМЕНА.1. ZMX(I)=(VPR*(1-Х)) *0 . 31. ZMDG=VPR*X1. TV=TX(I)

97. CALL PLJIDK(TV,ROJ) CALL PHAC(TV,P2)

98. WRITE(*,*) ' ДАВЛЕНИЕ ВПР. ПАРА=',P21. CALL PL(TV,P2,ROP)

99. WRITE(*,*) ' ROJ=',ROJ*1000

100. WRITE(*,*) ' ROP=',ROP*10001. CALL IROT(ROP,TV,EP)1. ZEP=EP*10 0 01. CALL IROT(ROJ,TV,EJ)1. ZEJ=EJ*10001. CALL IROT(RO,TG,EG)1. ZEG=EG*10 00

101. ZES=( (ZEG*ZMG(I)+ZEP*ZMDG)/(ZMG(I)+ZMDG) )/1000 WRITE(*,*) ' ЭНТАЛЬПИЯ CMECH=',ZES

102. ROS=((ZMG(I)+ZMDG)/(VKAM(I)-(ZMX(I)/(ROJ*100 0))))/1000 WRITE(*,*) ' ПЛОТНОСТЬ ГАЗА ПОСЛЕ ВПРЫСКА=',ROS*1000

103. WRITE(*,*) ' ИДУ ВЫПОЛНЯТЬ "TROI"'

104. CALL TROI1(ROS,ZES, TS) WRITE(*,*) ' ВЫПОЛНИЛ '

105. WRITE (*,*) ' ТЕМПЕРАТУРА ПОСЛЕ BnPbICKA=',TS1. CALL PF(ROS,TS,Z,PS1)

106. WRITE(*,*) ' ДАВЛЕНИЕ СМЕСИ=',PS1*1000001. ZTG(I)=TS1. PG(I)=PS11. VKX(9 0)=2.6 8 Е-131. VX(I)=ZMX(I)/(ROJ*1000)1. VKAM(I)=VKAM(I)-VX(I)

107. WRITE(*,*) ' ОБЪЕМ КАМЕРЫ СЖАТИЯ ПОСЛЕ ВПРЫСКА=',VKAM(I)1. ZMS=ZMG(I)+ZMDG

108. WRITE(*,*) ' МАССА СМЕСИ=',ZMS1. ZMG(I)=ZMS1. NKX(I)=VX(I)/VKX(I)

109. WRITE(*,*) ' ЧИСЛО КАПЕЛЬ ХЛАДАГЕНТА=',NKX(I)

110. FKX(I)=4.*3.14159*(DKX(I)/2)*(DKX(I)/2) CALL SALF(ZTG(I),RO,DKX(I),ALF2(I)) WRITE(*,*) ' ALFA=',ALF2(I)

111. QX(I)=NKX(I)*ALF2(I)*FKX(I)*(TKP-ZTG(I))*DTAU WRITE(*,*) ' QX=1,QX(I) CALL RPAP(TV,RP)

112. ZMISP(I)=(-1*QX(I))/(RP*1000) WRITE)*,*) ' ZMISP=',ZMISP(I)

113. TX(1+1)=TX(I) ZMG(1+1)=ZMG(I)+ZMISP(I)

114. WRITE(*,*) ' ZMG ПОСЛЕ ИСПАРЕНИЯ ЧАСТИ ХЛАДАГЕНТА=',ZMG(I+1)

115. ZMX(1+1)=ZMX(I)-ZMISP(I) VX(1+1)=ZMX(1+1)/(ROJ*10 0 0) VKX(1+1)=VX(1+1)/NKX(I)

116. DKX(1 + 1)=2*EXP(LOG( (3*VKX(I+1) ) /(4*3.141492) )/3)1. DNKX=NKX(I)/8 0 01. NKX(1+1)=NKX(I)-DNKX1. GOTO 24

117. НАГРЕВ И ИСПАРЕНИЕ ЖИДКОСТИ 365 CALL ТНАС(PI,ТКР)1. (TX(I).GE.TKP) GOTO 556 WRITE(*,*) ' ЖИДКОСТЬ НАГРЕВАЕТСЯ' WRITE(*,*) ' TX=',TX(I) TV=TX(I)1. CALL PLJIDK(TV,ROJ)

118. FKX(I)=4,*3.1415 9*(DKX(I)/2)*(DKX(I)/2) WRITE(*,*) ' FKX=1,FKX(I) CALL SALF(ZTG(I),RO,DKX(I),ALF2(I)) WRITE(*,*) ' ALFA=',ALF2(I)

119. QX (I) =NKX(I) *ALF2 (I) * FKX ( I ) * (TX (I) -ZTG (I) ) *DTAU1. WRITE(*,*) ' QX=',QX(I)1. CALL SCPX(PG(I),CPX(I))1. WRITE(*,*) ' MX=',ZMX(I)

120. WRITE(*,*) ' CP X=',CPX(I)

121. DTX(I)=QX(I)/(CPX(I)* ZMX(I) )

122. WRITE(*,*) ' DTX=',-DTX(I)1. TX(1+1)=TX(I)-DTX(I)1. (TX(I+1).GT.TKP) GOTO 5561. ZMG(1+1)=ZMG(I)1. VKX(I+1)=VKX(I)1. DKX(1+1)=DKX(I)1. NKX(1+1)=NKX(I)-DNKX

123. ZMX(1+1)=NKX(1+1)*VKX(1+1)*ROJ*10001. GOTO 2 424 CONTINUE

124. QM(I)=NKM(I)*ALF1(I)*FKM*(TM(I)-ZTG(I))*DTAU WRITE(*,*) ' ТЕПЛОТА К МАСЛУ=',QM(I) DTG(I)=(QX(I)+QM(I)+PG(I)*ZDV)/(ZMG(I)*CVG(I)) WRITE(*,*) ' ИЗМЕНЕНИЕ ТЕМПЕРАТУРЫ ГАЗА=',DTG(I) ZTG(1+1)=ZTG(I)+DTG(I)

125. WRITE(*,*) ' ТЕМПЕРАТУРА ГАЗА=',ZTG(1+1) DTM(I)=QM(I)/(ZMM*CPM)

126. WRITE(*,*) ' ИЗМЕН.ТЕМПЕРАТУРЫ МАСЛА=',DTM(I) TM(I+1)=TM(I)-DTM(I)

127. WRITE(*,*) ' ТЕМПЕРАТУРА МАСЛА=',TM(I+1) VKAM(I+1)=VKAM(I)-ZDV NKM(1+1)=NKM(I)-DNKM РАБОТА СЖАТИЯ

128. WRITE(*,*) ' P*dV=',PG(I)* ZDV RAB=RAB+PG(I)*ZDV ZR(I)=RAB

129. WRITE(*,*) ' РАБОТА В ТОЧКЕ=',RAB WRITE(*,*) ' РАБОТА СЖАТИЯ=',ZR(I) ,1 IF (PG(I).GT.1167500) GOTO 771. (PG(I).GT.1351500) GOTO 77 1 CONTINUE

130. GOTO 7 7 356 WRITE(*,*) ' БРЕД ! '

131. WRITE(*,*) ' СОЗДАНИЕ ФАЙЛА С РЕЗУЛЬТАТАМИ'

132. OPEN(9, FILE="C:\NAUKA\OTVVPRl.DAT",STATUS='OLD') DO 91 K=1,2 9 0

133. WRITE(9,121) K,ZTG(K),PG(K),ZR(K) L21 FORMAT (2X,I3,2X,'T=',F5.1,2X,'P=',E9.3,2X,'L=',F10.4) 91 CONTINUE CLOSE (9)1. WRITE{*,*) ' КОНЕЦ!'1. STOP1. END

134. SUBROUTINE SALF(TG,ZN,D,ALF) DOUBLE PRECISION T,RO,ETA,ZN,PR T=TG RO=ZN1. CALL LMDPAP(T,RO,ETA)1. BETA=0.51. VI=0.8E-61. U0=14.241. RE=(UO*D)/VI1. WRITE(*,*) ' RE=1,RE1. CALL PRDTLPAP(T,RO,PR)1. WRITE(*,*) ' PR=',PR

135. ZNU=2 * (1+(BETA*(SQRT(RE) )*(PR**(1/3) ) ) )1. WRITE(*,*) ' NU=',ZNU1. ALF=ZNU*ETA/D1. RETURN1. END1. SUBROUTINE SCPX(T,C)

136. DIMENSION Z(5),F(5),W(5),A(5),V(5)1. EXTERNAL ALI, ATSG

137. OPEN(6,FILE="C:\NAUKA\CPX.DAT",STATUS='OLD')1. READ(6,*) Z,F,EP

138. CALL ATSG(T,Z,F,W,5,1,A,V, 5)

139. CALL ALI(T,A, V, Y, 5, EP, IE)1. C=Y*1.E31. CLOSE (6)1. RETURN1. END

140. SUBROUTINE ATSG(X,Z,F,WORK,IROW,ICOL,ARG,VAL,NDIM) DIMENSION Z(1),F(1),W0RK(1),ARG(1),VAL(1) IF(IROW)11,11,11 N=NDIM1.(N-IROW)3,3,22 N=IROW3 B=0 .

141. DO 5 1=1,IROW DELTA=ABS(Z(I)-X) IF(DELTA-B)5,5,44 B=DELTA5 WORK(I)=DELTA B=B+1.1. DO 10 J=1,N1. DELTA=B1. DO 7 1=1,IROW1.(WORK(I)-DELTA)6,7,76 11 = 1 DELTA=WORK(I)7 CONTINUE ARG(J)=Z(II)1.(ICOL-1)8,9,8

142. VAL(2 * J-1)=ff(11) III=II+IROW1. VAL(2 * J)=F(III) GOTO 109 VAL(J)=F(II)10 WORK(II)=B11 RETURN END

143. SUBROUTINE ALI(X,ARG,VAL,Y,NDIM, EPS,IER)1. DIMENSION ARG(1),VAL(1)1.R=21. DELT2=0.1.(NDIM-l)9,7,1

144. DO 6 J=2,NDIM DELT1=DELT2 IEND=J-1

145. DO 2 1=1,IEND H=ARG(I) -ARG(J) IF(H)2,13,2

146. VAL(J)=(VAL(I)*(X-ARG(J))-VAL(J)*(X-ARG(I)))/H DELT2=ABS(VAL(J)-VAL(IEND))1.(J-2)6,6,33 IF(DELT2-EPS)10,10,44 IF(J-5)6,5,55 IF(DELT2-DELT1)6,11,116 CONTINUE7 J=NDIM8 Y=VAL(J)9 RETURN10 IER=0 GOTO 811 IER=112 J=IEND GOTO 813 IER=3 GOTO 12 END

147. Программа расчета отношения работы, затрачиваемой на сжатие в тнтовом маслозаполненном компрессоре, при впрыске в полость жидкого тбочего вещества к работе компрессора без впрыска жидкого хладагента.

148. ПРОГРАММА РАСЧЕТА ЗАВИСИМОСТИ ИЗМЕНЕНИЯ РАБОТЫ СЖАТИЯ ОТ ТЕМПЕРАТУРЫ КИПЕНИЯ И КОНДЕСАЦИИ.

149. РАБОЧЕЕ ВЕЩЕСТВО АММИМК, МАСЛО ХС-4 0, ВПРЫСКИВАЕТСЯ ЖИДКОЕ РАБОЧЕЕ ВЕЩЕСТВО С ТЕМПЕРАТУРОЙ +5 С1.TEGER I,J,К,L

150. DIMENSION РР(3,5),ТО(3,5),ТК(3,5),WCV(3,5),WCS(3,5) EXTERNAL SVPR,SSTD,ALI,ATSG, SALF, TROPE DO 3 K=1,5 TK(1,K)=25+273.15 TK(2,K)=35+273.15 TK(3,K)=40+273 . 15 3 CONTINUE DO 4 L=1,3

151. ТО(L,l)=-10+273.15 ТО(L,2)=-15+273.15 ТО(L,3)=-20+2 73.15 TO(L,4)=-25+273.15 ТО(L,5)=-30+273.15 4 CONTINUE1. DO 11 1=1,3 DO 11 J=l, 5

152. CALL SVPR(TO(I, J) ,TK(I,J),WCV(I, J) ) CALL SSTD(TO(I,J),TK(I,J),WCS(I,J)) PP (I, J) = ((WCV(I,J)/WCS(I,J))-1)*100 WRITE(*,*) ' PP=',PP(I,J),1,J11 CONTINUE

153. WRITE(*,*) ' СОЗДАНИЕ ФАЙЛА С РЕЗУЛЬТАТАМИ' OPEN(3,FILE="C:\NAUKA\REZGRF.DAT",STATUS=1 OLD') DO 91 K=l,5

154. WRITE(3,121) K,PP(1,K),PP(2,K),PP(3,K) L21 FORMAT (2X,I3,2X,F7.4,2X,F7.4,2X,F7.4) 91 CONTINUE CLOSE (3)1. WRITE(*,*) ' КОНЕЦ!'1. STOP1. END

155. SUBROUTINE SVPR(ZZ,ZTK,ZS)

156. ПРОГРАММА РАСЧЕТА ПРОЦЕССА СЖАТИЯ ВК С ВПРЫСКОМ ЖИДКОГО РАБОЧЕГО ВЕЩЕСТВА АММИАКА ,МАСЛО ХС 4 0 ,Т КИП=-15 ,Т КОНД=ЭО

157. DOUBLE PRECISION RO,TG,CV,TGI,TG2,Р,PI,ROJ,Z,TV,Р2,ROP, *EG,EJ,ZES,TS,ROS,EP,PS1,TKP,RP,ZKK,PK INTEGER I

158. WRITE(*,*) ' ЧИСЛО КАПЕЛЬ МАСЛА=',NKM(1) DNKM=01. DTAU=0.011/2901. ZDV=(3.8 5E-6)/2.91. WRITE(*,*) ' ZDV=',ZDV1. CALL R7171. TG1=ZTG(1)1. TG2=ZTG(1)1. ZKK=ZTK

159. CALL PHAC(ZKK, PK) CALL PHAC (TGI,P) CALL PL(TG2,P,RO)

160. WRITE(*,*) ' ПЛОТНОСТЬ ГАЗА НА ВСАСЫВАНИИ=',RO*1000 WRITE(*,*) ' ДАВЛЕНИЕ НА ВСАСЫВАНИИ=',Р*1000001. ZMG(1)=VKAM(1)*RO*l.E3

161. WRITE(*,*) ' МАССА НА ВСАСЫВАНИИ=ZMG(1)

162. DO 1 1=1,290 WRITE( 6,11) I 11 FORMAT( 2X,'ИДЕТ РАСЧЕТ ТОЧКИ НОМЕР ',13) RO=(ZMG(I)/VKAM(I))*l.E-3 CALL SALF(ZTG(I),RO,DKM,ALF1(I)) WRITE(*,*) ' ALFA=',ALF1(I) TG=ZTG(I)

163. CALL CVROT(RO,TG,CV) CVG(I)=CV*10 00

164. WRITE(*,*) ' ТЕПЛОЕМКОСТЬ=',CVG(I)

165. CALL PF(RO,TG,Z,PI) PG(I)=P1*100000

166. WRITE(*,*) ' ДАВЛЕНИЕ=',PG(I) IF (I.LT.90) GOTO 355 IF (I.EQ.90) GOTO 555 IF (I.GT.90) GOTO 5 6555 QX(I)=0.1. ZMG(1+1)=ZMG(1) GOTO 2 4

167. HA 9 0 ШАГЕ ВПРЫСКИВАЕТСЯ ЖИДКОСТЬ КОЛИЧЕСТВО ВПРЫСКИАЕМОГО ХЛАДАГЕНТА РАВНО VPR >55 DKX(90)=8.Е-5

168. CALL ТНАС(PI,ТКР) ТХ(I)=ТКР

169. WRITE(*,*) ' TEMPER XLAD=',ТХ(I) МАССА ВПРЫСКИВАЕМОЙ ЖИДКОСТИ VPR=(8.7 8Е-4)*1.1 Х=0.13

170. ТОЛЬКО 30 ПРОЦЕНТОВ ЖИДКОСТИ ОСТАЕТСЯ В ПОТОКЕ ГАЗА, ЭТА ЖИДКОСТЬ

171. УЧАСТВУЕТ В ПРОЦЕССЕ ТЕПЛООБМЕНА.1. ZMX(I)=( VPR*(1-Х))*0.31. ZMDG=VPR*X1. TV=TX(I)

172. CALL PLJIDK(TV,ROJ) CALL PHAC(TV,P2)

173. WRITE(*,*) ' ДАВЛЕНИЕ ВПР. ПАРА=',P21. CALL PL(TV,P2,ROP)

174. WRITE(*,*) ' ROJ=',ROJ*10 00

175. WRITE(*,*) ' ROP=',ROP*10 001. CALL IROT(ROP,TV,EP)1. ZEP=EP*10 0 01. CALL IROT(ROJ,TV,EJ)1. ZEJ=EJ*10001. CALL IROT(RO,TG,EG)1. ZEG=EG*1000

176. ZES=( (ZEG* ZMG(I)+ZEP*ZMDG)/(ZMG(I)+ZMDG) )/1000 WRITE(*,*) ' ЭНТАЛЬПИЯ СМЕСИ=',ZES

177. ROS=((ZMG(I)+ZMDG)/(VKAM(I)-(ZMX(I)/(ROJ*100 0))))/1000 WRITE(*,*) ' ПЛОТНОСТЬ ГАЗА ПОСЛЕ ВПРЫСКА=',ROS*1000

178. WRITE(*,*) ' ИДУ ВЫПОЛНЯТЬ "TROI"'

179. CALL TROU (ROS, ZES,TS) WRITE(*,*) ' ВЫПОЛНИЛ '

180. WRITE (*,*) ' ТЕМПЕРАТУРА ПОСЛЕ BnPbICKA=',TS1. CALL PF(ROS,TS, Z, PS1)

181. WRITE(*,*) ' ДАВЛЕНИЕ СМЕСИ=',PS1*100000

182. ZTG(I)=TS PG(I)=PS1 VKX(90)=2.68E-13

183. VX(I)=ZMX(I)/(ROJ*10 0 0) VRAM(I)=VKAM(I)-VX(I)

184. WRITE(*,*) ' ОБЪЕМ КАМЕРЫ СЖАТИЯ ПОСЛЕ ВПРЫСКА=',VKAM(I)1. ZMS=ZMG(I)+ZMDG

185. WRITE(*,*) ' МАССА СМЕСИ=',ZMS1. ZMG(I)=ZMS1. NKX(I)=VX(I)/VKX(I)

186. WRITE(*,*) ' ЧИСЛО КАПЕЛЬ ХЛАДАГЕНТА=',NKX(I)

187. FKX(I)=4. *3.14159*(DKX(I)/2)* (DKX (I)/2) CALL SALF(ZTG(I),RO,DKX(I),ALF2(I)) WRITE(*,*) ' ALFA=',ALF2(I)

188. QX(I)=NKX(I)*ALF2(I)* FKX(I)*(TKP-ZTG(I) )*DTAU1. WRITE(*,*) ' QX=',QX(I)1. CALL RPAP(TV,RP)

189. ZMISP(I) = (-1*QX (I) )/(RP*10 00)

190. WRITE(*,*) ' ZMISP=',ZMISP(I)1. TX(I+1)=TX(I)1. ZMG(1+1)=ZMG(I)+ZMISP(I)

191. WRITE(*,*) ' ZMG ПОСЛЕ ИСПАРЕНИЯ ЧАСТИ ХЛАДАГЕНТА=',ZMG(1+1)

192. ZMX(1+1)=ZMX(I)-ZMISP(I) VX(1+1)=ZMX(1+1)/(R0j*1000) VKX(1+1)=VX(1+1)/NKX(I)

193. DKX(1 + 1)=2*EXP(LOG( (3*VKX(I + 1) )/ (4*3.141492) )/3)1. DNKX=NKX(I)/8 001. NKX(1+1)=NKX(I)-DNKX1. GOTO 2 4

194. НАГРЕВ И ИСПАРЕНИЕ ЖИДКОСТИ >65 CALL ТНАС(PI,ТКР)1. (TX(I).GE.TKP) GOTO 556 WRITE(*,*) ' ЖИДКОСТЬ НАГРЕВАЕТСЯ' WRITE(*,*) ' TX=',TX(I) TV=TX(I)1. CALL PLJIDK(TV,ROJ)

195. FKX (I) =4 . *3. 1415 9* (DKX (I) /2) * (DKX (I) /2) WRITE(*,*) ' FKX=',FKX{I) CALL SALF(ZTG(I) ,RO,DKX(I) ,ALF2(I) ) WRITE(*,*) ' ALFA=',ALF2(I)

196. QX(I)=NKX(I)*ALF2(I)*FKX(I)*(TX(I)-ZTG(I))*DTAU1. WRITE(*,*) ' QX=',QX(I)1. CALL SCPX(PG(I),CPX(I))1. WRITE(*,*) ' MX=',ZMX(I)

197. WRITE(*,*) ' CP X=',CPX(I)

198. DTX(I)=QX(I)/(CPX ( I)* ZMX ( I) )

199. WRITE(*,*) ' DTX=',-DTX(I)1. TX(I+1)=TX(I)-DTX(I)1. (TX(I+1).GT.TKP) GOTO 5561. ZMG(1+1)=ZMG(I)1. VKX(I+1)=VKX(I)1. DKX(1+1)=DKX(I)1. NKX(1+1)=NKX(I)-DNKX

200. ZMX(1+1)=NKX(1+1)*VKX(1+1)*ROJ*10001. GOTO 2 424 CONTINUE

201. QM(I)=NKM(I)*ALF1(I)*FKM*(TM(I)-ZTG(I))*DTAU

202. WRITE(*,*) ' ТЕПЛОТА К МАСЛУ=',QM(I)

203. DTG(I) = (QX(I)+QM(I)+PG(I)* ZDV)/(ZMG(I)*CVG(I))

204. WRITE(*,*) ' ИЗМЕНЕНИЕ ТЕМПЕРАТУРЫ ГАЗА=',DTG(I) ZTG(1 + 1 )=ZTG(I)+DTG(I)

205. WRITE(*,*) ' ТЕМПЕРАТУРА ГАЗА=',ZTG(1+1) DTM(I)=QM(I)/(ZMM*CPM)

206. WRITE(*,*) ' ИЗМЕН.ТЕМПЕРАТУРЫ МАСЛА=',DTM(I) TM(I+1)=TM(I)-DTM(I)

207. WRITE( * , * ) ' ТЕМПЕРАТУРА МАСЛА=',ТМ(1+1) VKAM(1+1)=VKAM(I)-ZDV NKM(1+1)=NKM(I)-DNKM РАБОТА СЖАТИЯ

208. WRITE(*,*)' P*dV=',PG(I)*ZDV RAB=RAB+PG(I)*ZDV ZR(I)=RAB

209. WRITE(*,*) ' РАБОТА В ТОЧКЕ=',RAB WRITE(*,*) ' РАБОТА СЖАТИЯ=',ZR(I),I IF (PG(I).GT.1167500) GOTO 77 IF (PG(I).GT.PK*1.E5) GOTO 77 1 CONTINUE

210. GOTO 7 7 >56 WRITE(*,*)' БРЕД !'

211. WRITE(*,*) ' ВЫПОЛНЕНА ПОДПРОГРАММА SVPR' Z S=RAB RETURN END

212. SUBROUTINE SALF(TG,ZN,D, ALF) DOUBLE PRECISION T,RO,ETA,ZN,PR T=TG RO=ZN1. CALL LMDPAP(T,RO, ETA)1. BETA=0.51. VI=0.8E-61. UO=14.241. RE=(UO*D)/VI1. WRITE(*,*) ' RE=',RE1. CALL PRDTLPAP(T,RO,PR)1. WRITE(*,*) ' PR=',PR

213. ZNU=2 *(1+(BETA*(SQRT(RE) )*(PR**(1/3) ) ) )1. WRITE(*,*) ' NU=1,ZNU1. ALF=ZNU*ETA/D1. RETURN1. END1. SUBROUTINE SCPX(T,C)

214. DIMENSION Z(5),F(5),W(5),A(5),V(5)1. EXTERNAL ALI, ATSG

215. OPEN(6,FILE="C:\NAUKA\CPX.DAT",STATUS=1 OLD')1. READ(6,*) Z,F,EP

216. CALL ATSG(T,Z,F,W,5,1,A,V,5)1. CALL ALI(T,A,V,Y,5,EP,IE)1. C=Y*1.E31. CLOSE (6)1. RETURN1. END

217. SUBROUTINE ATSG(X,Z,F,WORK,IROW,ICOL,ARG,VAL,NDIM) DIMENSION Z(1) ,F(1) ,WORK(l),ARG(1) ,VAL(1) IF(IROW)11,11,11 N=NDIM1.(N-IROW)3,3,22 N=IROW3 B=0.

218. DO 5 1=1,IROW DELTA=ABS(Z(I)-X)1.(DELTA-B)5,5,44 B=DELTA5 WORK(I)=DELTA B=B+1.1. DO 10 J=1,N1. DELTA=B1. DO 7 1=1,IROW1.(WORK(I)-DELTA)6,7,76 11 = 1 DELTA=WORK(I)

219. CONTINUE ARG(J)=Z(II) IF(ICOL-l)8,9,8

220. VAL(2 * J-1)=F(II) II1=1I+1ROW

221. VAL (2 * J)=F(111) GOTO 109 VAL(J)=F(II)10 WORK(II)=B11 RETURN END

222. SUBROUTINE ALI(X,ARG,VAL,Y,NDIM,EPS,IER)1. DIMENSION ARG(1),VAL(1)1.R=21. DELT2=0.1.(NDIM-1)9,7,1

223. DO 6 J=2,NDIM DELT1=DELT2 IEND=J-1

224. DO 2 1=1,IEND H=ARG(I)-ARG(J) IF(H)2,13,2

225. VAL(J)=(VAL(I)*(X-ARG(J))-VAL(J)*(X-ARG(I)))/Н DELT2=ABS(VAL(J)-VAL(IEND))1.(J-2)6,6,33 IF(DELT2-EPS)10,10,44 IF(J-5)6,5,55 IF(DELT2-DELT1)6, 11, 116 CONTINUE7 J=NDIM8 Y=VAL(J)9 RETURN10 IER=0 GOTO 811 IER=112 J=IEND GOTO 813 IER=3 GOTO 12 END

226. SUBROUTINE SSTD(ZZ,ZTK,ZS)

227. ПРОГРАММА РАСЧЕТА ПРОЦЕССА СЖАТИЯ МАСЛ03АП0ЛНЕН0Г0 ВК

228. РАБ. ВЕЩЕСТВО АММИАК ,МАСЛО ХС 4 0 ,Т КИП=-10 ,Т КОНД=35,Т МАСЛА=35.

229. DOUBLE PRECISION RO,TG,CV,TGI,TG2,Р,PI, ZKK,РК, Z INTEGER I

230. DIMENSION ZTG(290),TM(290),VKAM(290),DKM(290),ALF1(290) , *PG(290),DTG(290),DTM(290),CVG(290),QM(290),NKM(290) ,ZR(290) EXTERNAL SALF,ALI,ATSG

231. WRITE(*,*) 'НАЧАЛ О' RG=48 8.16

232. RAB=0. ZTG(1)=ZZ TM(1)=30 8.15 VKAM(l)=4.68E-4 ZMM=(8.78E-4) VM=ZMM/82 7. DKM=1.E-4 VKM=5.24E-13 FKM=3.14E-8 CPM=2172. NKM(1)=VM/VKM

233. WRITE(*,*) ' ЧИСЛО КАПЕЛЬ МАСЛА=',NKM(1) DNKM=01. DTAU=0.011/2901. Z DV=(3.85E-6)/2.91. WRITE(*,*) ' Z DV=1,Z DV1. CALL R7171. TG1=ZTG(1)1. TG2=ZTG(1)1. ZKK=ZTK

234. CALL PHAC(ZKK,PK) CALL PHAC(TGI,P) CALL PL(TG2,P,RO)

235. WRITE(*,*) ' ПЛОТНОСТЬ ГАЗА НА ВСАСЫВАНИИ=',RO*1000 WRITE(*,*) ' ДАВЛЕНИЕ НА ВСАСЫВАНИИ=',P*100000 ZMG=VKAM(1)*RO*l.ЕЗ

236. WRITE(*,*) 1 МАССА НА ВСАСЫВАНИИ=',ZMG DO 1 1=1,290 WRITE(6,11) I 11 FORMAT( 2X,'ИДЕТ РАСЧЕТ ТОЧКИ НОМЕР ',13) RO=(ZMG/VKAM(I))*1.E-3 CALL SALF(ZTG(I),RO,DKM,ALF1(I)) WRITE(*,*) ' ALFA=',ALF1(I) TG=ZTG(I)

237. CALL CVROT(RO,TG,CV) CVG(I)=CV*1000

238. WRITE(*,*) ' ТЕПЛОЕМКОСТЬ=',CVG(I)

239. CALL PF(RO,TG,Z,PI) PG(I)=P1*100000

240. WRITE(*,*) ' ДАВЛЕНИЕ=1,PG(I)

241. QM(I)=NKM(I)*ALF1(I)*FKM*(TM(I)-ZTG(I))*DTAU

242. DTG(I)=(QM(I)+PG(I)*ZDV) * /(ZMG*CVG(I)) WRITE(*,*) ' ИЗМЕНЕНИЕ ТЕМПЕРАТУРЫ ГАЗА=',DTG(I) ZTG(I+1)=ZTG(I)+DTG(I)

243. WRITE(*,*) ' ТЕМПЕРАТУРА ГАЗА=',ZTG(1+1) DTM(I)=QM(I)/(ZMM*CPM)

244. WRITE(*,*) ' ИЗМЕН.ТЕМПЕРАТУРЫ МАСЛА=1,DTM(I) TM(I+1)=TM(I)-DTM(I)

245. WRITE(*,*) ' ТЕМПЕРАТУРА МАСЛА=',TM(I+1)1. VKAM(1+1)=VKAM(I)-ZDV1. NKM(1+1)=NKM(I)-DNKM1. RAB=RAB+PG(I)*ZDV1. ZR(I)=RAB

246. WRITE(*,*)' РАБОТА=',RAB IF (PG(I).GT.PK*100000) GO TO 77 1 CONTINUE

247. WRITE(*,*) ' ВЫПОЛНЕНА ПОДПРОГРАММА SSTD' ZS=RAB RETURN END