Моделирование тепломассообменных процессов и разработка жомосушильной установки с компрессией вторичного пара тема диссертации и автореферата по ВАК РФ 05.14.04, кандидат технических наук Рудаков, Юрий Игнатьевич

• Рудаков, Юрий Игнатьевич
• кандидат технических наук
• 2002

Оглавление диссертации кандидат технических наук Рудаков, Юрий Игнатьевич

Введение.

Обозначения.

Глава 1. Литературный обзор.

1.1. Общие сведения о свекловичном жоме.

1.2. Барабанная жомосушилка.

1.3. Паровая жомосушилка фирмы ВВ8.

1.4. Чанный испаритель (тостер) фирмы СКЕТ.

1.5. Постановка задач исследования.

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

2.1. Физические характеристики жома.

2.2. Исходные предпосылки и балансовые уравнения для секции испарителя.

2.3. Флегмовый режим.

2.4. Переходный режим.

2.5. Режим внутреннего испарения.

2.6. Определение входных и выходных параметров процесса в испарителе.

Глава 3. Устройство и работа паровых жомосушильных установок с компрессией вторичного пара.

3.1. Исходные предпосылки.

3.2. Установка с пароструйным компрессором.,.

3.3. Установка с турбокомпрессорным агрегатом.

3.4. Установки с рециркуляцией и перегревом пара в выносном калорифере.

3.5. Комбинированный способ сушки.

Глава 4. Расчет процесса в установке и анализ результатов.

4.1. Геометрические характеристики частиц жома.

4.2. Коэффициенты теплообмена в испарителе.

4.3. Оптимальный режим работы испарителя.

4.4. Тепловой расчет системы очистки вторичного пара.

4.5. Расчет коэффициента инжекции пароструйного компрессора.

4.6. Тепловой расчет системы охлаждения сжатого пара.

4.7. Пример расчета совместной работы испарителя и системы кондиционирования вторичного пара.

4.8. Расчетный анализ работы установки при комбинированном способе сушки.

4.9 Расчет экономических показателей.

Выводы.

Введение диссертации (часть автореферата) на тему «Моделирование тепломассообменных процессов и разработка жомосушильной установки с компрессией вторичного пара»

Актуальность темы. В настоящее время на российских сахарных заводах жом сушат в барабанных сушилках. При этом для получения 1 т сухого жома сжигают примерно 0.6 т условного топлива [1] и в атмосферу в смеси с отработавшими топочными газами выбрасывают от 4 до 5 т сокового пара при температуре 150° С. В связи с большой энергоемкостью такого способа сушки и постоянно растущими ценами на энергоносители многие сахарные заводы прекратили выпуск сухого жома. Между тем, сухой жом является хорошим кормом для сельскохозяйственных животных, удобным для транспортирования и длительного хранения. Известно, что свекловичный жом, высушенный в среде перегретого водяного пара, является экологически чистым продуктом, пригодным для приготовления различных пищевых добавок, используемых в кондитерской, мясной и других отраслях пищевой промышленности [2]. На рынках западных стран качественный жом продается по цене 100 $ за тонну. Возросшие цены на продукцию животноводства, а также расширившиеся возможности экспорта сухого жома на выгодных условиях возобновляют интерес производителей сахара к сушке жома. Однако тенденция повышения цен на энергоносители сохраняется [3]. В проекте государственной энергетической стратегии России на период до 2020 года предусматривается быстрый рост цен на газ: в 2.5 раза (при льготном налогообложении) к 2003 году и еще в 1.4 раза к 2005 году (без учета инфляции), а к 2007 году предполагается вывести их на уровень равновесия с ценами газа на европейском рынке, что означает рост примерно в 3.7-4 раза по сравнению с ценами 2001 года. Из сказанного выше очевидна актуальность разработки энергосберегающей и экологически чистой технологии сушки жома.

Наряду с этим важной задачей является совершенствование всего теплосилового хозяйства сахарного завода. Типовой сахарный завод располагает собственной ТЭЦ, оснащенной парогенераторами высокого давления (3.5 МП а) и турбогенераторами, снабжающими производство электроэнергией и ретурным (отработавшим в турбинах) паром низкого давления (0.2-0.3 МПа). Дефицит ретурного пара (достигающий 40%) устраняют с помощью редук-ционно-охладительной установки (РОУ), в которой пар высокого давления дросселируют и охлаждают водой [4]. Использование турбин с противодавлением и РОУ также приводит к излишнему расходу топлива.

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

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

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

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

Апробация. Основные положения диссертации обсуждались на отчетных научных конференциях ВГТА (2001, 2002 г.), региональных межвузовских научно-тематических семинарах "Моделирование процессов тепло- и массообмена" и "Физико-технические проблемы энергетики и экологии" при ВГТУ (2001, 2002 г.), а также на техническом совете ЗАО "Кристалл Групп".

Публикации. По материалам диссертации опубликовано 7 работ, в том числе 3 патента на изобретения.

Структура и объем работы. Диссертация содержит введение, список обозначений, четыре главы, выводы, библиографический список и приложения. Она изложена на 133 страницах, включает 23 рисунка и 7 таблиц. Библиографический список содержит 60 наименований.

Заключение диссертации по теме «Промышленная теплоэнергетика», Рудаков, Юрий Игнатьевич

ВЫВОДЫ

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

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

3. Производительность испарителя по высушенному жому зависит от параметров и величины подачи греющего пара. Потребление глухого греющего пара ограничено сверху допустимой величиной давления в рубашках испарителя (1.0 МПа) и коэффициентом теплопередачи (кА~ 125 Вт/{мАК)), поэтому основным средством повышения производительности является увеличение подачи в испаритель перегретого острого пара. Величина подачи острого пара ограничена явлением уноса мелких частиц жома в мокрый циклон. Оптимальный режим работы испарителя определяется с учетом указанных факторов.

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

5. На сахарных заводах, располагающих барабанными жомосушилка-ми, экономически целесообразно применение комбинированного способа сушки, согласно которому в барабанных сушилках жом высушивают топочными газами до промежуточной влажности 40-60%, после чего жом перегружают в испаритель с паровым обогревом и досушивают до конечной влажности 12%. Этот способ позволит при сохранении той же производительности по сухому жому существенно сократить расход топлива на сушку или при том же расходе топлива значительно повысить производительность. При этом следует ожидать повышения качества и снижения потерь жома, так как выгоранию подвержены достаточно сухие частицы, которые в комбинированном способе сушки не вступают в контакт с горячими топочными газами, а находятся в среде перегретого водяного пара.

6. Реализация данного способа сушки на Чернянском сахарном заводе даст прибыль ~ 0.25 млн. $ в год, что соизмеримо с затратами на проектирование, приобретение необходимого оборудования и монтаж установки, следовательно, срок окупаемости затрат равен одному году. При использовании современного жомового пресса типа 81ог(1 2500 размер прибыли увеличится до 1.2 млн. $ в год при том же сроке окупаемости.

Список литературы диссертационного исследования кандидат технических наук Рудаков, Юрий Игнатьевич, 2002 год

1. Сшт n.M. Технология сахара. М.: Пищевая пром-сть, 1967. - 625 с.

2. Щеренко А.П. Сушка жома перегретым водяным паром: экономично, перспективно // Сахар. 2001. - №6. - С. 18-20.

3. Белобров В. Реформирование отечественной электроэнергетики. Проблемы энергоснабжения промышленных объектов // Аналитический журнал РЦБ. Товарный рынок.-2001. №6 (07). - С.43-47.

4. Колесников В.А., Нечаев Ю.Г. Теплосиловое хозяйство сахарных заводов. -М.: Пищевая пром-сть, 1980. 392 с.

5. Технология свеклосахарного производства / Под ред. Р.А. Мак-Джинниса. М.: Пищепромиздат, 1958. - 488 с.

6. Знаменский Г.М. Технологическое оборудование сахарных заводов. М.: Пищепромиздат, 1952. - 424 с.

7. Гришин М.А. Установки для сушки пищевых продуктов / М.А. Гришин, В.И. Анатазевич, Ю.Г. Семенов. -М.: Агропромиздат, 1989. 215 с.

8. Стандарт предприятия, пуск, наладка и испытание жомосушильных установок (СТП-0561 -37-83). Киев: Укрсахпром, 1983.

9. Технология производства растительных масел / Под ред. В.М. Копейков-ского и СИ. Данилъчука. М,: Легкая и пищевая пром-сть, 1982, - 416 с.

10. Масликов В.А. Технологическое оборудование производства растительных масел. М.: Пищевая пром-сть, 1974. - 440 с.

11. Парфенопуло М.Г. Исследование процесса сушки свекловичного жома: Дисс. канд. техн. наук / Воронежский технол. ин-т. Воронеж, 1967.

12. Парфенопуло М.Г Исследование процесса сушки свекловичного жома: Автореф. дисс. . канд. техн. наук / Воронежский технол. ин-т. Воронеж, 1967.

13. ЛыковA.B. Теория сушки. М.: Энергия, 1968. - 472 с.

14. Гинзбург A.C. Технология сушки пищевых продуктов. М.: Пищевая пром-сть, 1976. - 248 с.

15. Харин В.М., Рудаков Ю.И. Объемная усадка и изменение плотности влажных материалов при сушке. // Хран. и перераб. сельхозсырья. -2002. № . -С.

16. Антипов С. Т. Интенсификация процесса сушки масличных семян в аппарате с вращающимся барабаном: Дисс. . канд. техн. наук. Воронежский технол. ин-т, 1983.

17. Шишацкий Ю.И. Сушка хлебопекарных дрожжей. Воронеж: Изд-во Воронежского гос. ун-та, 1992. - 180 с.

18. Чубик И.А., Маслов A.M. Справочник по теплофизическим константам пищевых продуктов и полуфабрикатов.-М.: Пищевая пром-сть, 1965.-156 с.

19. Аравии В.И., Нумеров С.Н. Теория движения жидкостей и газов в не деформируемой пористой среде. М.: Гостехтеориздат, 1953 .-618с.

20. Харин В.М., Агафонов ГВ. Теоретические основы тепло- и влагообмен-ных процессов пищевой технологии. М: Пищевая пром-сть, 2001. - 344 с.

21. Харин В.М., Агафонов Г.В. Теория гигро- и гидротермической обработки капиллярнопористых тел. Воронеж: изд-во ВГТА, 2000. - 184 с.

22. Харин В.М., Агафонов Г.В. Внешний влаго- и теплообмен капиллярнопо-ристого тела с газо-паро-жидкостной средой // Теор. основы хим. технол. -1999. Т.ЗЗ. - №3. - С.252-258.

23. Харин В.М., Агафонов ГВ., Рудаков Ю.И. Кинетика влаго- и теплообмена между капиллярнопористым телом и аэрозолем // Теор. основы хим. тех-нол. 2002. - Т.36. - №4. - С. 1-9.

24. Кафаров В.В. Методы кибернетики в химии и химической технологии.-М.: Химия, 1976.-464 с.

25. Закгейм А. Ю. Введение в моделирование химикотехнологических про-цессов.-М.: Химия, 1973,-223 с.

26. Федоткин И.М. Методы расчета реакторов пищевой технологии.- Киев: Вища школа, 1978.-200 с.

27. Безденежных A.A. Математические модели химических реакторов.- Киев: Техшка, 1970.-322С.

28. Харин В.М., Агафонов Г.В. Внешний влаго- и теплообмен капиллярнопо-ристого тела с газо-паровой средой // Теор. основы хим. технол. 1999. -Т.ЗЗ.-№2.-С. 144-152.

29. Харин В.М., Агафонов Г.В., Горяинов A.A. Внутренний влаго- и теплопе-ренос в капиллярнопористых телах // Теор. основы хим. технол. 2000. -Т.34.-№5.-С.520-525.

30. Харин В.М., Агафонов Г.В., Горяинов A.A. Кинетика гигротермической обработки капиллярнопористых материалов // Теор. основы хим. технол. -2001.-Т35.-№1.-С.12-20.

31. ЛыковA.B., МихайловЮ.А. Теория тепло- и массопереноса. М. - Л.: Гос-энергоиздат, 1963. - 536 с.

32. Ъ2. Харин СЕ. Физическая химия.-Киев: Изд-во Киевск. ун-та, 1961. 408 с.

33. Курс физической химии, т. 1 / Я.И. Герасимов, В.П. Древинг, Е.Н.Еремин и др. М.: Госхимиздат, 1963. - 624 с.

34. Вукалович МП. Таблицы теплофизических свойств воды и водяного пара / М.П. Вукалович, СЛ. Ривкин, A.A. Александров. М.: Изд-во стандартов, 1969. 408 с.

35. Михайлов Ю.А. Сушка перегретым паром. М.: Энергия, 1967. - 200 с.

36. Харин В.М., Рудаков Ю.И. Повышение энергетической эффективности жомосушильных установок // Вестник ВГТУ. Сер. Энергетика. 2001. ВЫП.7.1.-С.51-56.

37. Харин В.М., Рудаков Ю.И. Тепловой расчет жомосушильной установки с чанным испарителем и пароструйным компрессором // Там же. С.61-63.

38. Установка для сушки жома: Пат.2178866 РФ / В.М. Харин, Ю.И. Рудаков, М.И. Кобрисев, М.В. Харин №2001109627/06; Заявл. 10.04.2001. Опубл. 27.01.2002.-Бюл.№3.

39. Харин В.М., Рудаков Ю.И. Тепловой расчет жомосушильной установки с чанным испарителем и паровым турбокомпрессором // Вестник ВГТУ. Сер. Энергетика. 2001. Вьш.7.1. С.57-60.

40. Установка для сушки жома: Пат.2178867 РФ / В.М. Харин, Ю.И. Рудаков, М.И. Кобрисев, М.В. Харин №2001109701/06; Заявл. 10.04.2001. Опубл. 27.01.2002.-Бюл.№3.

41. Паровая сушильная установка / В.М: Харин, Ю.И. Рудаков, М.И. Кобрисев, М.В. Харин. Заявка №2002106376/13(006731).

42. Паровая сушильная установка / В.М. Харин, Ю.И. Рудаков, М.И. Кобрисев, М.В. Харин. Заявка №2002106377/13(006732).

43. Способ сушки свекловичного жома / В.М. Харин, Ю.И. Рудаков, М.И. Кобрисев, М.В. Харин №2001112208/13(012627); Заявл. 03.05.2001. - Решение ФИПС о выдаче патента 13.03.2002.

44. Аэров М.Э., Тодес О.М. Гидравлические и тепловые основы работы аппаратов со стационарным и кипящим зернистым слоем.- Л.: Химия, 1968.512 с.

45. Кондратьев Г.М. Регулярный тепловой режим. М.: Гостехиздат, 1954. -408 с.

46. Кондратьев Г.М. Тепловые измерения. М. - Л.: Машгиз, 1964. - 244 с.

47. Лыков A.B. Теория теплопроводности. М.: Высшая школа, 1967. - 599 с.

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

49. ПерриДж. Справочник инженера-химика. Т.2.- Л.: Химия, 1969. 504 с.114

50. Исаченко В.П., Осипова В.А., Сукомел А.С. Теплопередача.- М.: Энергия,1969.-440 с.

51. Теплотехника /У4.Я. Баскаков, Е.В. Берг, O.K. Витт и др. М.: Энергоиз-дат, 1982. - 264 с.

52. Теплотехника / A.M. Архаров, СИ. Исаев, И.А. Кожинов и др. М.: Машиностроение, 1986. - 432 с.5 3. Коган В.Б. Теоретические основы типовых процессов химической технологии. Л.: Химия, 1977. - 592 с.

53. Касаткин А.Г. Основные процессы и аппараты химической технологии. -М.: Химия, 197 3.-752 с.

54. Стабников В.Н., Попов В.Д., Лысянский В.М., Редъко Ф.А. Процессы и аппараты пищевых производств. М.: Пищевая пром-сть, 1976. - 663 с.

55. Бакластов A.M. Проектирование, монтаж и эксплуатация теплоиспользующих установок. М.: Энергия, 1970. - 568 с.

56. Волощук В.М., Седунов ЮС. Процессы коагуляции в дисперсных системах. Л.: Гидрометеоиздат, 1974. - 3 20 с,5 8. Седунов Ю.С Физика образования жидко-капельной фазы в атмосфере. -Л.: Гидрометеоиздат, 1972. 208 с.

57. ХргианА.Х. Физика атмосферы. Т.2.- Л.: Гидрометеоиздат, 1978. 248 с.

58. Соколов Е.Я., Зингер КМ. Струйные аппараты. М.: Энергия, 1970.-288 с.

59. Программа расчета параметров процесса в испарителе

60. GLS :'ПрограммаPRIL1 (QuickBASIK4.5)

61. Постоянные READ pi, G, Ry, Fr DATA 3.14,9.81,8314,1

62. Физические свойства воды и водяного пара1. READk,m,pt,Tt,rt,Ts

63. DATA 1.3,18,610.8,273.15,2501000,373.151. READ Aw, Bw, Cw, Pw, Qw

64. DATA9.248,27.098,2.005,-.2224,-.1201

65. DEFfiir(T) = rt*(H-Pw*(T/Tt- l) + Qw*(T/Tt- 1)Л2)

66. DEF fops (T) = pt * (Tt / T) л Aw * EXP(Bw * (1 Tt / T) - Cw * (1 - T / Tt))ps = fnps(Ts): r = fiir(Ts): Dp = m / Ry / Ts2 * ps * r

67. DEF fopv (T) = ps + Dp * (T Ts)

68. DEF folc (T) = .0595 * SQR(T / Tt) / (1 + 2.46 * Tt / T)

69. DEF fonc (T) = .0000291 * (T / Tt) л (3 / 2) / (T / Tt + 2.38)

70. DEF foroc (p, T) = m / Ry * p / T

71. DEF fovc (p, T) = fonc(T) / foroc(p, T)

72. DEF fovl (T) = 2.42 * 10 л -(8 * (T 171) / (T - 140))

73. Физические свойства материала READ omz, bv, cO, cpr, usl, us2, bt DATA 1360,1.34,2116,le-14,3, .15,.035

74. DEF foro (uv) = omz * (uv + 1) / (bv * uv + 1)

75. DEF foom (uv) = omz / (1 + bv * uv) DEF foD (uv) = 4.63E-08 + 2.37E-06 * uv DEF fodelta (uv) = .072 * (1 (1 - uv / 2.15) 2) DEF foups2 = us2 - bt * (T - Ts)

76. Геометрические характеристики частиц READ11,12,13, Wstr DATA .00075,.004,.04,75 ustr = Wstr/(100-Wstr)

77. Характеристики чанного испарителя (тостера)1. READZ,D,Ft,kt1. DATA 10,2.1,80,125

78. S = pi * D л 2 / 4: kgm = kt: kgc = kt: pc = ps

79. Характеристики слоя материала READ Hsl, eps, kp DATA .4,.4,.8

80. Vsl = S * Hsl: Vm = (1 eps) * Vsl: Fmz = Vm * Fz / Vzdez = 4 * eps / (1 eps) * Vz / Fz

81. Ftl = Ft / Z: Fgm = (1 eps) * Ftl: Fgc = eps * Ftl

82. DEF fiiFm (uv) = Fmz / (1 + bv * uv) л (1 / 3)

83. DEF fnFcm (uv) = kp * fnFm(uv)

84. DEF fiide (uv) = dez * (1 + bv * uv) (1 / 3)

85. Параметры рабочего, глухого и острого греющего пара READpl,Tl,hp, pg, Tg,Tos DATA 3500000,708.15,3303400,1000000,691,685 hg = hp: ho = hp

86. Параметры конденсата греющего пара READ Tgs,hkg, his DATA 453,762600,417510

87. Параметры насыщенного пара и конденсата атмосферного давления1. READ hs, cs, hk, cl

88. DATA 2676500,2040,417510,4190roc = fnroc(ps, Ts): vc = fnvc(ps, Ts): Ic = fiilc(Ts)1. Pre = cs * roc * vc / Ic

89. Удельная теплоемкость острого пара с2 = (ho hs) / (Tos - Ts)

90. CLS : PRINT "Введите (через запятую) входные параметры материала-" PRINT "температуру Тп в градусах Цельсия, влажность Wn в %" PRINT "и нажмите клавишу Enter"

91. PRINT "Tn, Wn": INPUT Tn, Wn•РАСЧЕТ ПРОЦЕССА В ИСПАРИТЕЛЕ Tn = Tn + Tt:un = Wn/(100-Wn) Wk=12:uk = Wk/(100-Wk)

92. Допусимый размер уносимых фракций 1у = .15 * СШТ(1000 * Iz) /1000 'Расчет параметра optJcгок = fiiro(uk): Ary = G * ly л 3 * (rok roc) / vc 2 / roc

93. Rey = Ary / (18 + .61 * SQR(Aiy)): wy = Rey * vc / lyoptJc = S * wy * roc: ws = optJc / Z / roc / S1. Qgt = Ft * kt * (Tgs Ts)

94. Расчет параметров optJo и optGswczn = cO + cl * un: rz = r + czn * (Ts Tn) / (un - uk)optJo = (rz * optic Qgt) / (rz + ho - hs)

95. OptGsw = (optJc optlo) / (un - uk)

96. Итерационный расчет процесса в испарителе

97. Ga = .7 * OptGsw: Gb = 1.3 * optGswm4: Gsw = (Ga + Gb) / 2: GOTO ml4ml 5: IF uv(Z) > uk THEN Gb = Gsw ELSE Ga = Gsw1. (Gb Ga) / Gsw > .00001 THEN GOTO m41. GOTO m60ml4: 'Ядро расчета

98. N = l:uO = un:TO = Tn:ifO = 0m5: io = optJo / Z / Gsw: V = fiiV(uO): F = niF(uO)

99. Fm = niFm(uO): Fcm = fnFcm(uO)de = fnde(uO): dz = fndz(uO): om = fnom(uO)vl = fnvl(TO): v2 = fnvc(ps, TO)

100. Al = cz + xo + r * ks * (del2 + H2 * delp * Dp)

101. D1 = cz * TO r * ks * (uO - us2 - (bt + H2 * delp * Dp) * Ts)1. Zn = Al *B2-A2*B 1

102. Tv = (B2*Dl-Bl *D2)/Zn:Tf=(Al * D2 A2 * Dl)/Znuf = us2 bt * (Tf - Ts): pv = ps + H2 * Dp * (Tv - Ts)uv = (uO ks2 * (del2 * (Tv - Tf) + H2 * delp * (pv - ps) - uf)) / (1 + ks2)

103. N < Z THEN N = N + 1: uO = uv: TO = Tv: ifO = ifw ELSE GOTO ml 51. RETURNm60: 'Суммирование массовых и тепловых потоков Jc = 0: Jo = 0: Jg = 0:Qc = 0 FORN = lTOZ

104. Jc = Jc + Gsw * ic(N): Jo = Jo + Gsw * io(N): Jg = Jg + Gsw * ig(N) Qc = Qc + Gsw * ic(N) * hc(N)1. NEXTNhe = Qc / Jc: Tc = Ts + (he hs) / cs Gk = Gsw * (1 + uv(Z)): Jk = Jg

105. Проверка теплового баланса испарителя Qg = Jg * (hg hkg): Qo = Jo * (hg - he) DQn = Jg * (hg - hkg) + Jo * hg - Jc * he DQs = Gsw * (cO * (Tn - Tv(Z)))

106. DQw = Gsw * (cl * (un * (Tn Tt) - uv(Z) * (Tv(Z) - Tt))) eQ = (DQn + DQs + DQw) / (Qo + Qg)

107. Печать таблицы данных расчета CLS : LOCATE 1, 55: PRINT "Таблица"

108. PRINT" Основные параметры процесса в испарителе при:"

109. PRINT " рс="; рсо; "Па,"; " Тн="; Тп Tt; "С, WH="; Wn; "% (ин="; uno; "),"; PRINT MWK=M; Wk;"%, ly="; 1000 * ly; "мм" FOR J = 1 TO 73: PRINT "";: NEXT J: PIONT "" PRESrr " N"; " P "; " H";

110. PRINT" uv Tv,C Tc,C dpv,na 12 ic if ig кем"

111. FOR J = 1 TO 73: PRINT "-"; : NEXT J: PRINT "-"1. FORN = lTOZ1. PRINT USING "##";N;

112. PRINT " "; r$(N);""; H(N);"";

113. PRINT USING "#.###"; uvo(N);

114. PRLNTTUSING" ###.#"; tvo(N);

115. PRINT USING" ###.#"; tco(N);

116. РРШТ USING " #####"; dpvo(N);

117. PRINT USING" ####"; i2o(N);

118. PRINT USING "#.###"; ico(N);

119. PRINT USING " ####"; ifo(N);

120. PRINT USING "#.###"; igo(N);

121. PRINT USING "###.#"; kcmo(N)1. NEXTN

122. FOR J = 1 TO 73: PIONT "-"; : NEXT J: PRINT "-"121

123. PRINT "N-номера секций (начиная с верхней); Р-название режима: Ф-флегмовый,"

124. PRINT "П-переходный / номер версии, В-режим внутреннего испарения; Н=Н1,Н2."

125. PRINT "Потребление греющего пара:";

126. PRINT "глухого Jr="; CINT(3600 * Jg) / 1000; ",";

127. PRINT " острого Jo="; CINT(3600 * Jo) / 1000; "т/ч."

128. PRINT "Выход вторичного пара Jc="; CINT(3600 * Jc) /1000; "т/ч,";

129. PRINT" его температура Тс="; CINT(100 * (Тс Tt)) / 100; "С."

130. PRINT "Производительность испарителя GCB="; CINT(3600 * Gsw) / 1000; " и1. GK = ";

131. PRINT CINT(3600 * Gk) /1 ООО; "т/ч."

132. PRINT "Теплоподвод от глухого и острого пара:";

133. PRINT " Qr="; CINT(Qg / 1000); "и Qo="; CINT(Qo / 1000); "кВт."

134. PRINT "(Размерность кем Вт/(мЛ2*К)."1. END

135. Программа расчета производительности испарителя и расхода острого пара в зависимости от начальной влажности и температуры жома

136. CLS : 'Программа PRIL2 (QuickBASIC 4.5) READ pi, G, Ry DATA 3.14,9.81,8314

137. Физические свойства воды и водяного пара1. READm,pt,Tt,rt,Ts,cl

138. DATA 18,610.8,273.15,2501000,373.15,41901. READ Aw, Bw, Cw, Pw, Qw

139. DATA9.248,27.098,2.005,-.2224,-.1201

140. DEF for (T) = rt * (1 + Pw * (T/Tt- l) + Q w * (T/Tt-1 )Л2)

141. DEF fops (T) = pt * (Tt / T) л Aw * EXP(Bw * (1 Tt / T) - Cw * (1 - T / Tt))ps = f^ps(Ts): r = for(Ts)

142. DEF fonc (T) = .0000291 * (T / Tt) л (3 / 2) / (T / Tt + 2.38)

143. DEF foroc (p, T) = m / Ry * p / T

144. DEF fovc (p, T) = fonc(T) / foToc(p, T)

145. Физические свойства материала1. READ omz, bv, cO1. DATA 1360,1.34,2116

146. DEF foro (uv) = omz * (uv + 1) / (bv * uv + 1)

147. Характеристики чанного испарителя1. READ D, Ft,kt1. DATA 2.1,80,125

148. S = pi * D л 2 / 4: pc = ps1. Параметры греющего пара1. READ pi, Tl, hp

149. DATA 3500000,708.15,3303400hg = hp: ho = hp

150. Параметры конденсата греющего пара READ Tgs, hkg DATA 453,762600

151. Параметры насыщенного пара атмосферного давления1. READhs1. DATA 2676500

152. РАСЧЕТ ПРОЦЕССА В ИСПАРИТЕЛЕ 'Задаваемые постоянные параметры READ Тп, Wk, Iz DATA323.15,12,.001973

153. Допусимый размер уносимых фракций 1у = .15 * CINT(1000 * Iz) / 1000 'Расчет параметров Jc,Qg,Jg DEF feuv (W) = W / (100 W)uk = fhuv(Wk): rok = fhro(uk): roc = fiiroc(ps, Ts): vc = fnvc(ps, Ts)

154. Ary = G * ly а 3 * (rok roc) / vc а 2 / roc

155. Rey = Ary / (18 + .61 * SQR(Ary)): wy = Rey * vc / ly

156. Jc = S * wy * roc: Qg = Ft * kt * (Tgs Ts): Jg = Qg / (hg - hkg)

157. DEF fncz (W) = cO + cl * fnuv(W)

158. DEF fnrz (W) = r + &cz( W) * (Ts Tn) / (fnuv(W) - uk)

159. DEF fnJo (W) = (fnrz(W) * Jc Qg) / (fhrz(W) + ho - hs)

160. DEF fiiGsw (W) = (Jc * (ho hs) + Qg) / (fiiuv(W) - uk) / (fhrz(W) + ho - hs)

161. DEF fhGk (W) = (1 + uk) * fhGsw(W)

162. Wl = 40: W2 = 90: Tn = 323.15: maxJo = fnJo(Wl): maxQo = maxJo * (ho hs)minJo = (r * Jc Qg) / (r + ho - hs): minQo = minJo * (ho - hs)1. Графика

163. KEY OFF: LOCATE 25, 1: CLS : SCREEN 12 mGk=12

164. VIEW (48, 20)-(200,270): WINDOW (Wl, 1.2 * mGk)-(W2, 0)

165. NE (Wl, 1.2 * mGk)-(Wl, 0): LINE -(W2, 0)

166. CATE 1,6: PRINT "Gk,Jo, т/ч"

167. CATE 17, 5: PRINT "O": LOCATE 18, 5: PRINT Wl

168. NE (Wl, mGk)-(1.03 * Wl, mGk): LINE (Wl, mGk / 2)-(1.03 * Wl, mGk / 2) LOCATE 11,4: PRINT mGk / 2: LOCATE 4, 3: PRINT mGk LINE (W2, 0)-(W2, .02 * mGk): LOCATE 18, 14: PRINT " Wn, % "; W2 DW = (W2-Wl)/20

169. GOSUB ml: Tn = Tn + 50: GOSUB ml: GOTO m2ml: LINE (Wl, 3.6 * fhJo(Wl))-(Wl, 3.6 * fnJo(Wl))

170. FOR W = Wl TO W2 STEP DW: LINE -(W, 3.6 * fiiJo(W)): NEXT W

171. AE (Wl, 3.6 * fnGk(Wl))-(Wl, 3.6 * fnGk(Wl))

172. PmST "Qo="; CINT(minQo / 1000);"-"; CINT(maxQo / 1000); "кВт." END

173. Программа расчета основных параметров пароструйного компрессора в зависимости от давления сжатого пара

174. CLS : Программа PRIL3 (QuickBASIC 4.5)

175. Физические свойства воды и водяного пара1. READ Ry, к, m1. DATA83 14,1.3481. Параметры рабочего пара1. READpl,Tl,hl

176. DATA 3500000,708.15,3303400

177. Параметры инжектируемого пара1. READps, Ts, hs

178. DATA 101291,373.15,26765001. Газодинамические функции1. Ш Р Ж ( 1 ) = (1 -(1/lm^

179. DEF fiiq (1) = bz * 1 * (fim(l)) л (1 / к)

180. DEF fiiFz (о) = (1 о л 2) * о ( к -1)1. Постоянные параметры1. READfil,fi2,fi3,fi41. DATA .95,.975,.9,.925

181. Массив параметров сжатого пара

182. READ р(0), р(1), р(2), р(3), р(4), р(5)

183. DATA.20,.21,.22,.23,.24,.25

184. READ р(6), р(7), р(8), р(9), р(10)1. DATA .26,.27,.28,.29,.30

185. READ Ts(0), Ts(l), Ts(2), Ts(3), Ts(4), Ts(5)

186. DATA 120.23,121.78,123.27,124.71,126.09,127.43

187. READ Ts(6), Ts(7), Ts(8), Ts(9), Ts( 10)

188. DATA 128.73,129.98,131.20,132.39,133.54

189. READ hs(0), hs(l), hs(2), hs(3), hs(4), hs(5)

190. DATA 2706900,2709200,2711300,2713300,2715300,2717200

191. READ hs(6), hs(7), hs(8), hs(9), hs(lO)

192. DATA 2719000,2720700,2722300,2723900,2725500

193. FOR i = 0 TO 10: p(i) = p(i) * 1000000: Ts(i) = Ts(i) + Tt: NEXT i

194. Тасчет коэффициента инжекции1. FORi = 0Td0p2 = p(i): Ts2 = Ts(i): hs2 = hs(i)nCH = ps / p2: nPC = p2/ pi: nHC = 1 / nCH: mqC3 = pi / p2 * qPS 'Итерации l a=l: lb = .4ml: 1C3 = (la + lb) / 2: GOSUB m91. и > Unp2 THEN la = 1C3 ELSE lb = 1C3

195. ABS(U Unp2) > .0001 THEN GOTO ml ELSE GOTO ml9m9: qC3 = fnq(lC3): nC3 = fim(lC3)

196. КЗ = I + fi3 * nPC * (nC3 nCH) / к / nz / 1C3 / qPH

197. K4 = 1 + fi3 * nHC * (nC3 nC2) / к / nz / IC3 / qH2и = (Kl * IPH КЗ * 1C3) / (K4 * 1C3 - K2 * 1H2) / kot1. RETURNml8: 0 = (dm + db) / 2: IF fiiFz(o) > A THEN db = 0 ELSE dm = о

198. ABS(db dm) > .000001 THEN GOTO ml 8:1. RETURNml9: hen = (hi + U * hs) / (1+ U)

199. Табуляция результатов расчета

200. U(i) = U: hcn(i) = hen: dh(i) = hcn(i) hs(i)1. NEXT i1. Графика

201. KEY OFF: LOCATE 25, 1: CLS : SCREEN 12 mh = 3.1

202. VIEW (48, 30)-(200, 190): WINDOW (p(0), mh)-(p(10), 2.5)

203. NE (p(0), mh)-(p(0), 2.5): LINE -(p(10), 2.5)

204. CATE 1, 15: PRINT "(a)": LOCATE 2, 5: PRINT 'Ъ,МДж/кг"

205. NE (p(0), hs(0) / 1000000)-(p(0), hs(0) / 1000000)

206. FOR i = 0 TO 10: LINE -(p(i), hs(i)/ 1000000): NEXT i

207. NE (p(0), hcn(O) / 1000000)-(p(0), hcn(O) / 1000000)

208. FOR i = 0 TO 10: LINE -(p(i), hcn(i) / 1000000): NEXT i

209. NE (p(0), 3)-(1.01 * p(0), 3): LINE (p(10), 0)-(p(10), 2.515)

210. NE (p(10), 0)-(p(10), .03): LINE (p(8), 0)-(p(8), .03)

211. NE (p(6), 0)-(p(6), .03): LINE (p(4), 0)-(p(4), .03)

212. NE (p(2), 0)-(p(2), .03): LINE (p(0), U(0))-(p(0), U(0))

213. FOR i = 0 TO 10: LINE -(p(i), U(i)): NEXT i

214. CATE 13,32: PMNT "2 p2, бар 3"

215. Печать таблицы данных расчета1. Округление данных

216. Т1о = СШТ(Т1 Tt): Tso = CINT(Ts - Tt)pto = CINT(pl / 10000) / 100: pso = CINT(ps / 10000) / 100

217. FOR i = 0 TO 10: Uo(i) = CINT(100 * U(i)) / 100hcno(i) = CINT(hcn(i) / 10000) /100: hso(i) = CINT(hs(i) / 10000) / 1001. NEXTi

218. CATE 16, 57: PIONT "Таблица"

219. PRINT" Основные параметры пароструйного компрессора в зависимости"; PRINT" от давления сжатия при:";

220. PRINT " р1 = "; pio;", ps="; pso; "МПа, TI = "; Tio;", Ts = "; Tso; "C" FOR J = 1 TO 76: PRINT"";: NEXT J: PRINT "" PRINT " Значения | ",, "Давление сжатого пара p2, МПа" PIONT "параметров!";

221. FOR i = О TO 10: PRINT USING " #.## "; CINT(p(i) /10000) / 100;: NEXT i PIUNT : FOR J = 1 TO 76: PRINT"-";: NEXT J: PRINT"-" PRESrT" и ";

222. FOR i = 0 TO 10: PRINT USING "?¥.##"; Uo(i);: NEXT i PRINT : PRINT " Ьсп,МДж/кг ";

223. FOR i = 0 TO 10: PRINT USING " #.## "; hcno(i);: NEXT i PRINT '^,МДж/кг";

224. FOR i = 0 TO 10: PRINT USING " #.## "; hso(i); : NEXT i: PMNT FOR J = 1 TO 76: PRESIT "";: NEXT J: PRINT ""1. END

225. Программа расчета основных параметров работы установки

226. CLS : Программа PRIL4 (QuickBASIC 4.5)1. Постоянные1. READ pi, g, Ry1. DATA 3.14,9.81,8314

227. Физические свойства воды и водяного пара1. READ к, т, pt, Tt, rt, Ts

228. DATA 1.3,18,610.8,273.15,2501000,373.151. READ Aw, Bw, Cw, Pw, Qw

229. DATA 9.248,27.098,2.005,-.2224,-. 1201

230. DEFfnr(T) = rt*(l + Pw*(T/Tt-l) + C ! w*(T/Tt-1 )Л2)

231. DEFihps (T) = pt * (Tt/T)лAw* EXP(Bw* (1 -Tt/T)-Cw* (1 -T/Tt))ps = fnps(Ts): r = &ir(Ts)

232. DEF fimc (T) = .0000291 * (T / Tt) л (3 / 2) / (T / Tt + 2.38)

233. DEF fiiroc (p, T) = m / Ry * p / T

234. DEF fevc (p, T) = fnnc(T) / fiiroc(p, T)

235. Физические свойства материала1. READ omz, bv, cO1. DATA 1360,1.34,2116

236. DEF fhro (Uv) = omz * (Uv + 1) / (bv * Uv + 1) 'Характеристики чанного испарителя READ D, Ft, kt DATA 2.1,80,125 S = pi * D 2 / 4: pc = ps

237. Параметры рабочего, глухого и острого греющего пара1. READpl,Tl,hl

238. DATA 3500000,708.15,3303400hp = hl: hg = hi: ho = hl

239. Параметры конденсата греющего пара1. READ Tgs, hkg1. DATA 453,762600

240. Параметры насыщенного пара и конденсата атмосферного давления1. READ hs, cs, hk, cl

241. DATA 2676500,2040,417510,4190

242. Задаваемые постоянные параметры

243. READ Tn, Wn, Wk, ly DATA 373.15,50,12,.00031. Tc = Ts

244. РАСЧЕТ ОСНОВНЫХ ПАРАМЕТРОВ РАБОТЫ ИСПАРИТЕЛЯun = Wn / (100 Wn): uk = Wk / (100 - Wk)rok = friro(uk): roc = fhroc(ps, Ts): vc = fnvc(ps, Ts)

245. Ary = g * ly 3 * (rok roc) / vc A 2 / roc

246. Rey = Ary / (18 + .61 * SQR(Ary)): wy = Rey * vc / ly

247. Jc = S * wy * roc: Qg = Ft * kt * (Tgs Ts): Qo = Jc * (ho - hs): Q = Qg + Qo

248. Jg = Qg / (hg hkg): rh = r + ho - hs: Jo = (r * Jc - Qg) / rh

249. Gsw = (Jc Jo) / (un - uk): Gk = Gsw * (1 + uk)

250. РАСЧЕТ СИСТЕМЫ ОЧИСТКИ ПАРА1. READ Jpw, Tpw, Txw1. DATA .278,363.15,293.15hxw = cl * (Txw Tt): Jkl = Jpw * (Tpw - Txw) / (Tgs - Txw)

251. Jwl = Jpw * (Tgs Tpw) / (Tgs - Txw): U l = Jwl / Jkl

252. Js = (Jc * (r + cs * (Tc Ts)) + Jpw * c l * (Tpw - Ts)) / r

253. Jsl = (Jc * cs * (Ts Tc) + Jpw * (r - c l * (Tpw - Ts))) / r

254. РАСЧЕТ ОСНОВНЫХ ПАРАМЕТРОВ ПАРОСТРУЙНОГО1. КОМПРЕССОРА1. Газодинамические функции

255. DEFfhn(l) = (l-(l/lm)A2)A(k/(k-l))

256. DBF fnq (1) = bz * 1 * (fim(l)) a (1 / к)

257. DEF fiiFz (о) = (1 о 2 ) * о ( к - 1)1. Постоянные параметры1. READfil,fi2,fi3,fi41. DATA .95,.975,.9,.925

258. H = Im * SQR(1 nPH a ((k - 1) / k)): qPH = bz * IPH * nPH a (1 / k)nPS = nPH*nz:qPS = bz*lm* SQR(1 -nPSA((k-1)/^) * nPS a (1/к)qH2m = bz * (1 oz 2) (1 / (k - 1))

259. Массив параметров сжатого пара

260. READ р(0), р(1), р(2), р(3), р(4), р(5)

261. DATA.20,.21,.22,.23,.24,.25

262. READp(6),p(7),p(8),p(9),p(10)1. DATA .26,.27,.28,.29,.30

263. READ hs(0), hs(l), hs(2), hs(3), hs(4), hs(5)

264. DATA 2706900,2709200,2711300,2713300,2715300,2717200

265. READ hs(6), hs(7), hs(8), hs(9), hs(lO)

266. DATA 2719000,2720700,2722300,2723900,2725500

267. FOR i = О ТО 10: p(i) = p(i) * 1000000: NEXT i 'Расчет

268. FOR i = 0 TO 10: p2 = p(i): hs2 = hs(i)nCH = ps / p2: nPC = p2 / pi: nHC = 1/nCH: mqC3 = p 1 / p2 * qPS1. Итерацииla = l: lb = .4m7:1C3 = (la + lb)7 2: GOSUB m91. и > Unp2 THEN la = 1C3 ELSE lb = 1C3

269. ABS(U Unp2) > .0001 THEN GOTO m? ELSE GOTO ml9m9: qC3 = fhq(lC3): nC3 = fnn(lC3)

270. КЗ = 1 + fi3 * nPC * (nC3 nCH) / к / nz / 1C3 / qPH

271. K4 = 1 + fi3 * nHC * (nC3 nC2) / к / nz / 1C3 / qH2

272. U = (K1 *1PH-K3 *lC3)/(K4*lC3-K2*IH2)/kot1. RETURNml8: 0 = (dm + db) / 2: IF fnFz(o) > A THEN db = о ELSE dm = о

273. ABS(db dm) > .000001 THEN GOTO ml8:1. RETURNml9: Jp = Js/U: Jcn = (l + 1/U) * Js: hen = (hp + U * hs)/ (1 +U) 'РАСЧЕТ СИСТЕМЫ ОХЛАЖДЕНИЯ СЖАТОГО ПАРА Jk2 = Jg Jkl: Jks = Jen * (hen - hs2) / (hs2 - hkg)

274. U(i) = U: Jp(i) = Jp: Jcn(i) = Jen: hcn(i) = hen: Jks(i) = Jks: Jw2(i) = Jw2

275. J2(i) = J2: U2(i) = U2: Jcm(i) = Jem: Jl(i) = Jl: Ke(i) = Ke1. NEXTi1. Графика

276. KEY OFF: LOCATE 25, 1: CLS : SCREEN 12 mJ = CINT(3.6* J2(10))

277. VIEW (48, 150)-(200, 395): WINDOW (p(0), mJ)-(p(10), 0) LINE (p(0), mJ)-(p(0), 0): LESOE-(p(10), 0)

278. CATE 9, 4: PRINT "JB2,Jp,Jl,Jcn,J2, т/ч": LOCATE 25, 5: PRINT "O"

279. NE (p(0), mJ)-(1.01 * p(0), mJ): LOCATE 10, 3: PRESTTmJ

280. NE (p(0), mJ / 2)-(1.01 * p(0), mJ / 2): LOCATE 18, 3: PRINT mJ / 2

281. NE (p(10), 0)-(p(10), .015 * mJ): LINE (p(8), 0)-(p(8), .015 * mJ)

282. NE (p(6), 0)-(p(6), .015 * mJ): LINE (p(4), 0)-(p(4), .015 * mJ)

283. NE (p(2), 0)-(p(2), .015 * mJ): LOCATE 26, 6: PRINT "0.2 p2, МПа 0.3

284. NE (p(0), 3.6* Jw2(0))-(p(0), 3.6* Jw2(0))

285. FOR i = 0 TO 10: UNE -(p(i), 3.6 * Jw2(i)): NEXT i1.CATE 23, 21: PRBvJTT

286. NE (p(0), 3.6* Jp(0))-(p(0), 3.6 * Jp(0))

287. FOR i = 0 TO 10: LINE -(p(i), 3.6 * Jp(i)): NEXT i1.CATE 17,21: PREЛJT "2"

288. NE (p(0), 3.6 * Jl(0))-(p(0), 3.6 * J1(0))

289. FOR i = 0 TO 10: UNE -(p(i), 3.6 * Jl(i)): NEXT i1.CATE 13, 27: Р1ШТ "3"

290. SIE (p(0), Jcn(0))-(p(0), Jcn(O))

291. FOR i = 0 TO 10: LINE -(p(i), 3.6 * Jcn(i)): NEXT i

292. CATE 12, 27: PRINT "4": LINE (p(0), 3.6 * J2(0))-(p(0), 3.6 * J2(0))

293. FOR i = 0 TO 10: LINE -(p(i), 3.6 * J2(i)): NEXT i1.CATE 11, 21: PRINT "5"

294. CATE 20, 35: PREvfT "1 JB2": LOCATE 21, 35: PRINT "2 - Jp" LOCATE 22, 35: PRINT "3 - Л ": LOCATE 23,35: PIONT "4 - Jen" LOCATE 24, 35: PMNT "5 - J2" 'Округления

295. CATE 28, 1: PRINT "Заданные параметры: ";

296. PRINT " TH="; Tn Tt; "C, WH="; Wn;"%, WK="; Wk; "%,"

297. PRINT," JnB="; CINT(360 * Jpw)/100; "т/ч, Тпв="; Tpw Tt; "C,";

298. PRINT "TxB="; Txw Tt; "C."

299. PRINT "Расчетные параметры:";

300. PRINT "Jr="; CINT(360 * Jg) /100; "т/ч, Jo="; CINT(360 * Jo) / 100; "т/ч, ";

301. PRINT "Jc="; CINT(360 * Jc) /100; "т/ч,";

302. PRINT "Js="; CINT(360 * Js) / 100; "т/ч, "

303. PRINT ," JB1="; CINT(360 * Jwl) /100; "т/ч,";

304. PRINT "JK1="; CINT(360 * Jkl) /100; "т/ч,";

305. PRINT "1шл="; CINT(360 * Jsl) / 100; "т/ч,";

306. PRINT "Jk2="; CINT(360 * Jk2) /100; "т/ч"131

307. PRINT, " GcB="; CINT(3 60 * Gsw)/100; "т/ч,"; PRINT "GK="; CINT(360 * Gk)/ 100; "т/ч."

308. PRINT : PRINT " При изменении давления вторичного пара от 0.2 до 0.3 МПа:"

309. PRINT "Jp="; Jpo(O);"-"; Jpo(lO); "т/ч,"; PRINT "Jcn="; Jcno(O);"-"; Jcno(lO); "т/ч,"; PRINT "JB2="; Jw2o(0);"-"; Jw2o(10); "т/ч," PRINT "Jl="; Jlo(O); "-"; Jlo(lO); "т/ч,"; PRINT "J2="; J2o(0); "-"; J2o(10); "т/ч,";

310. PRINT "Ke="; CINT(1 000 * Ke(0)) / 1000; "-"; CINT(1000 * Ke(lO)) /1 000 END

311. Программа расчета основных параметров работы комбинированной установки

312. CLS : 'Программа PRIL5 (QuickBASIC 4.5)

313. Параметры рабочего, глухого и острого греющего пара1. READpl,Tl,hp

314. DATA 3500000,708.15,3303400hg = hp:ho = hp

315. Параметры конденсата греющего пара1. READTgs,hkg1. DATA 453,762600

316. Параметры насыщенного пара атмосферного давления READTs,hs,r

317. DATA 373.15,2676500,2257109

318. Характеристики барабанной сушилки и испарителя1. READ Jwb, Ft, kt1. DATA 2.222,80,1251. Постоянные параметры1. READZb,Zi,Wk1. DATA 2,1,12uk = Wk/(100-Wk)

319. Qg = kt * Ft * (Tgs Ts): Jg = Qg / (hg - hkg)1. Определение функций

320. DEF fiiug (Wg) = Wg / (100 Wg)

321. DEF J5iun (Wg) = (Zb * Jwb * uk + Zi * Jwi * ftiug(Wg)) / Jw

322. DEF &Wn (Wg) = 100 * &un(Wg) / (1 + fiiun(Wg))

323. DEF fiiGsw (Wg) = Jw / (fiiug(Wg) uk)

324. DEF fiiGk (Wg) = fnGsw(Wg) * (1 + uk)

325. Построение графиков ¥н(¥ж)

326. Wgl = 70: Wg2 = 85: Wnl = 30: Wn2 = 60

327. KEY OFF: LOCATE 25, 1: CLS : SCREEN 12

328. VIEW (48, 30)-(200, 220): WEAJDOW (Wgl, Wn2 + l)-(Wg2 + 1, Wnl)

329. NE (Wgl, Wn2 + l)-(Wgl, Wnl): LINE -(Wg2 + 1, Wnl)

330. CATE 3, 3: PRESTT Wn2: LOCATE 14, 3: PRINT Wnl

331. CATE 15, 5: PRINT Wgl: LOCATE 15, 23: PRINT Wg21.CATE 2, 4: PRINT "wh,%"1.CATE 15, 15: PMNT '^ж,%"1. FORj = lT 0 3 : d W = 5*j1.sfE (Wgl + dW, Wnl)-(Wgl + dW, Wnl + .5): NEXT j1. FORj = lTO3 : dW= 10 *j

332. NE (Wgl, Wnl + dW)-(Wgl + .5, Wnl + dW): NEXT j LOCATE 3, 23: PIONT "1": LOCATE 7, 23: PRINT "4" Ъариируемая вличина подачи острого пара в испаритель (т/ч) READ Jot(l), Jot(2), Jot(3), Jot(4) DATA 8,10,12,14 'Пересчет Jo в кг/с

333. FOR i = 1 TO 4: Jo(i) = Jot(i) /3.6: NEXT i

334. FOR i = 1 TO 4: Jo = Jo(i): GOSUB m5: GOTO mlOm5:

335. Qo = Jo * (ho hs): Q = Qg + Qo

336. Jc = Jo + Q / г: Jwi = Q / r: Jw = Zb * Jwb + Zi * Jwi1. RETURNmlO: LINE (Wgl, fhWn(Wgl))-(Wgl, &Wn(Wgl)) FOR Wg = Wgl TO Wg2 STEP 1 LINE -(Wg, fnWn(Wg)): NEXT Wg NEXT!

337. Построение графиков gk(wk)" mG= 10

338. VIEW (265, 30)-(417, 220): WINDOW (Wgl, mG + 1.2)-(Wg2 + .5, 4)

339. NE (Wgl, mG + 1.2)-(Wgl, 4): LINE -(Wg2 + .5,4)

340. CATE 4, 30: PRINT mG: LOCATE 14, 31: РИШТ 4

341. CATE 15,32: PRE4T Wgl: LOCATE 15,51: PRESTT Wg2 .

342. CATE 2, 31: PRINT "Ок,т/ч"1.CATE 15, 42: PIONT "¥ж,%"1. FORj = lT0 3 :dW = 5*j

343. ЕЕЖ (Wgl + dW, 4)-(Wgl + dW, 4.15): NEXT j1. FORj = l T04:dG = 2*j

344. ЕЕЖ (Wgl, 4 + dG)-(Wgl + .5,4 + dG): NEXT j

345. FOR i = 1 TO 4: Jo = Jo(i): GOSUB m5

346. NE (Wgl, 3.6 * fnGk(Wgl))-(Wgl, 3.6 * fhGk(Wgl))

347. FOR Wg = Wgl TO Wg2 STEP 1: LINE -(Wg, 3.6 * foGk(Wg)): NEXT Wg1. NEXT i

348. CATE 9, 40: PRINT " 1": LOCATE 6, 40: PRINT "4"

349. CATE 18, 2: PRINT "Основные характеристики работы комбинированной установки"

350. PRESTT " Z6="; Zb;", ZH="; Zi;", Jwb="; CINT(3.6 * Jwb); "т/ч";

351. PIONT ", Jr="; CINT(360 * Jg) /100; "т/ч"

352. PRESIT " Jo= 1 -"; Jot(l);", 2 -"; Jot(2);

353. PRESIT 3 -"; Jot(3);", 4 -"; Jot(4); "т/ч"1. END