автореферат диссертации по энергетике, 05.14.04, диссертация на тему:Повышение эффективности энергоиспользования на промышленном предприятии на основе оптимизации теплотехнологических установок

Введение

1. Анализ ТЭБ промышленного предприятия и постановка задачи исследований

1.1 ТЭБ АО "Автокран"

1.2 Анализ ТЭБ АО "Автокран" с позиций интенсивного энергосбережения

1.3 Анализ путей повышения эффективности энергоиспользования на промышленных предприятиях

Выводы по главе

2. Математическая модель режима работы термической печи с выдвижным подом ; ч

2.1 Анализ существующих математических моделей

2.2 Постановка задачи

2.3 Математическое описание сопряженного теплообмена в рабочем пространстве печи

2.4 Алгоритм решения сопряженной задачи теплообмена при нагреве углеродистых сталей

1) период нагрева

2) период выдержки

3) период выгрузки-загрузки

2.5 Алгоритм решения сопряженной задачи теплообмена при нагреве высоколегированных сталей

1) период нагрева

2.6 Сопоставление результатов математического моделирования с данными промышленного эксперимента

Выводы по главе

3.Разработка и математическое моделирование установки по использованию ВЭР за ВТУ периодического действия

3.1 Разработка схемы утилизации теплоты продуктов сгорания для регенеративного и внешнего теплоиспользования

3.2 Математическое описание теплообмена в рекуператоре с переходным режимом работы

3.3 Решение задачи радиационно-конвективного теплообмена в системе газ-кладка-ряд труб

3.4 Алгоритм расчета рекуператора Выводы по главе

4. Оптимизация и исследование работы ТТУ периодического действия

4.1 Исследование работы ТТУ периодического действия с регенеративным теплоиспользованием

4.1.1 Постановка задачи

4.1.2 Алгоритм решения задачи

4.2 Постановка задачи и критерии оптимизации ТТУ

4.3 Алгоритм выбора оптимального варианта теплообменника

4.4 Оптимизация режимов работы ТТУ периодического действия Выводы по главе

5. Практическое применение результатов исследований

5.1 Выбор оптимальной конструкции теплообменника для регенеративного и внешнего теплоиспользования Выводы по главе Выводы по диссертации Список литературы Приложение

АКТУАЛЬНОСТЬ РАБОТЫ. В настоящее время энергопотребление обеспечивается в основном за счет сжигания органического топлива в установках различных мощностей и конструкций, оказывающих негативное влияние на экологическую обстановку.

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

За счет утилизации ВЭР при современном уровне техники можно ежегодно получать до 15,6 ГДж тепловой энергии [1]. Большинство ВЭР отрасли относятся к тепловым с температурой от 400 до 1400 °С, при этом выход продуктов сгорания от отдельной технологической установки составляет 300-10000 м /час. У значительной части производственных агрегатов кузнечных, термических, литейных цехов, потребляющих высококачественное топливо, в частности природный газ, коэффициент использования этого топлива равен 10-30 % [1]. Это обусловливает разработку принципиально новых систем утилизации теплоты на промышленных предприятиях и проведение мероприятий по повышению эффективности энергоиспользования на основе оптимизации теплотехнологических установок (ТТУ).

ЦЕЛЬЮ РАБОТЫ является проведения комплексного энергетического анализа промышленного предприятия с оценкой возможных масштабов энергосбережения и исследования путей повышения эффективности энергоиспользования за счет утилизации 5 отходящих газов, разработка и оптимизация ТТУ с внешним или внутренним теплоиспользованием.

НАУЧНАЯ НОВИЗНА работы:

1. Проведен комплексный энергетический анализ машиностроительного предприятия на основе критериев оценки качества использования топливно-энергетических ресурсов (ТЭР) в теплотехнологических системах (TTC).

2. Известный ранее метод решения сопряженных задач теплообмена - метод дискретного удовлетворения краевых условий (ДУКУ) -распространен на случай решения трехмерной нестационарной сопряженной задачи теплообмена в системе газ- кладка- металл. На основе этого разработана математическая модель оптимизации режима работы термической печи с выдвижным подом при ее эксплуатации в условиях конкретного производства.

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

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

5. Разработаны алгоритмы:

- выбора оптимального варианта ТТУ с внешним или внутренним теплоиспользованием;

- оптимизации рекуператора для действующей ТТУ с внешним или внутренним теплоиспользованием из условия максимального годового экономического эффекта;

- оптимизации ТТУ с внутренним теплоиспользованием отходящих газов с учетом динамической взаимосвязи печи и рекуператора;

- исследования совместной работы печи периодического действия и конвективного петлевого рекуператора с целью обеспечения надежности работы ТТУ.

МЕТОДЫ ИССЛЕДОВАНИЙ. Прикладные результаты работы получены с использованием методов планирования и анализа эксперимента и методов технико-экономической оптимизации. При математическом моделировании сопряженного теплообмена в термической печи периодического действия используется метод ДУКУ, модель радиационно-конвективного теплообмена в системе газ- кладка- ряд труб рекуператора построена на решении системы балансовых уравнений, записанных для каждого расчетного сечения.

ПРАКТИЧЕСКАЯ ЦЕННОСТЬ работы заключается в следующем:

- большие резервы экономии ТЭР, выявленные на примере машиностроительного предприятия, стимулируют на проведение глубокого энергетического анализа топливно-энергетических балансов (ТЭБ) предприятий аналогичного профиля;

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

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

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

- разработаны алгоритмы сопоставления и оптимизации схем ТТУ с внешним и внутренним теплоиспользованием, которые могут быть использованы при проектировании ТТУ для предприятий аналогичного профиля;

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

РЕАЛИЗАЦИЯ РЕЗУЛЬТАТОВ РАБОТЫ. Разработаны конструкция эффективного теплообменника и схема утилизации теплоты высокотемпературных газов. Выполнен рабочий проект установки внешнего теплоиспользования и передан для внедрения на АО "Автокран". Проведена оптимизация режима работы термической печи с выкатным подом.стальцеха АО "Автокран".

ОСНОВНЫЕ ПОЛОЖЕНИЯ ВЫНОСИМЫЕ НА ЗАЩИТУ:

- схема утилизации теплоты высокотемпературных газов для внешнего или внутреннего теплоиспользования;

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

- математическая модель конвективного петлевого рекуператора;

- результаты сопоставления данных расчета нагрева садки металла в термической печи периодического действия с промышленным экспериментом;

- алгоритм сопоставления схем с внешним и внутренним теплоиспользованием;

- результаты исследования совместной работы термической печи периодического действия и конвективного петлевого рекуператора.

ДОСТОВЕРНОСТЬ разработанной математической модели термической печи периодического действия проверена сопоставлением результатов расчета на модели с данными промышленного эксперимента.

АПРОБАЦИЯ РЕЗУЛЬТАТОВ РАБОТЫ. Основные результаты работы докладывались и обсуждались на международной конференции "Экология и теплотехника - 1996" (Днепропетровск, 1996), на международной научно-технической конференции "Состояние и перспективы развития электротехнологии" (IX Бенардосовские чтения, Иваново, 1999г.), на международной конференции "Проблемы печной теплотехники" (Днепропетровск, 1999).

ПУБЛИКАЦИИ. По материалам диссертации опубликовано 8 печатных работ.

ОБЪЕМ И СТРУКТУРА РАБОТЫ. Диссертационная работа состоит из введения, пяти глав, заключения, списка литературы из 152 наименований и 4 стр. приложений. Изложена на 142 стр., содержит 6 таблиц, 13 рисунков.

Выводы по диссертации

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

Наиболее существенными результатами выполненной работы являются:

1. Составлен ТЭБ машиностроительного предприятия "АО Автокран" и выявлены основные виды и количество тепловых отходов. В результате анализа ТЭБ предприятия с позиций интенсивного энергосбережения установлено, что технологический КПД совокупности ТТУ находится на уровне 19,7 %, а интегральный коэффициент полезного использования первичного топлива в TTC производства условного технологического продукта равен 14,7 %. Проведен анализ существующих схем утилизации теплоты высокотемпературных газов и путей повышения эффективности энергоиспользования на промышленных предприятиях.

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

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

Представлена математическая модель оптимизации режима работы термической печи периодического действия с выкатным подом для нагрева насыпных садок из углеродистых и легированных сталей при ее эксплуатации в условиях конкретного производства. Алгоритм модели построен на базе решения трехмерной нестационарной сопряженной задачи в системе газ- кладка- металла. На математической модели выполнен расчет для условий промышленных испытаний печи с выдвижным подом стальцеха АО "Автокран". Сопоставление всех имеющихся экспериментальных данных с расчетными показало, что величина среднего квадратичного отклонения значений температуры от данных промышленного эксперимента не превышает 10 %. Представлена математическая модель конвективного петлевого рекуператора. Модель радиационно-конвективного теплообмена в системе газ- кладка- ряд труб рекуператора построена на решении системы балансовых уравнений, записанных для каждого расчетного сечения.

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

Дан алгоритм выбора оптимального варианта теплообменника. Поиск оптимального варианта осуществляется по величине годового экономического эффекта, который определяется отдельно для внешнего и внутреннего теплоиспользования.

Проведено исследование совместной работы термической печи периодического действия и конвективного петлевого рекуператора. Предложено ввести в стоимость готового продукта неявные затраты на развитие технологического процесса в данном случае - на отопление цеха. Применение предложенного критерия позволит выбирать оптимальное направление теплоиспользования. Рассмотрена задача оптимизации режима работы ТТУ периодического действия. В качестве целевой функции принята стоимость операции нагрева.

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

1. Новгородский Е.Е., Пермяков Б.А. Энергосберегающие установки на машиностроительных заводах // Промышленная энергетика. -1995.- № 7. с.4-6.

2. Вторичные энергоресурсы и энерготехнологическое комбинирование / Семененко H.A., и др. Киев.:Выща школа, - 1979. - 295 с.

3. Григорьев В.Н. Повышение эффективности использования топлива в промышленных печах. М.Металлургия, - 1977. - 287 с.

4. Коротин А.Н., Бровкин Л.А., Лукьянов В.П. К выбору напрвления использования ВЭР // Промышленная энергетика. 1988. - № 10. - с.

5. Бровкин Л.А., Крылова Л.С. К решению уравнения теплопроводности при нелинейных граничных условиях // Изв. вузов. Энергетика. 1971. - № 9. - с. 122-127.

6. Андреев Е.И. Приближенный метод расчета тепломассообмена между газом и пленкой жидкости // Инж. физ. журн. 1987. Т. 53. -№2. - с. 191.

7. Гузов Л.А. Нестационарная теплопроводность в многослойном шаре // Математическое моделирование и оптимизация процессов тепломассообмена в установках промышленной энергетики: Межвуз. сб. науч. тр. / Иван. гос. ун-т. Иваново, 1983. - с. 20-24.

8. Применение асимптотических методов к решению задач нестационарной теплопроводности / Ю.И. Дударев, А.П. Кашин, В.И. Лозбин и др. // Инж. физ. журн. 1982. Т. 42. - № 3. - с. 25. -Деп. в ВИНИТИ 01.10.81, №4665-81.

9. Бровкин JI.A., Глазов B.C. // Изв. вузов. Энергетика. 1987. - № 7. -с. 104-108.

10. Копрински С.Ц. Об одной задаче теплопроводности // Инж. физ. журн. 1980. Т.39.- № 1. - с. 161 - Деп. в ВИНИТИ 28.11.79, № 30180 Деп.

11. Копрински С.Ц. О некоторых задачах теории тепло и массопереноса разрешимых с помощью преобразований Лапласа /У Инж. физ. журн. 1980. Т.

12. Видин Ю.В. Расчет нагрева излучением термически тонкого тела при переменной температуре среды // Изв. вузов. Черная металлургия. -1986. -№3.- с. 117-118.

13. Шабловский О.Н. // Инж. физ. журн. 1981-1982. Т. 40. - № 3. - с. 510-517.

14. Коляно Ю.М., Грицько Е.Г. Пространственная нестационарная задача теплопроводности для призмы при зависящем от координат коэффициенте теплоотдачи // Инж. физ. журн. 1982. Т. 42. - № 1. -с. 133-137.

15. Гайнутдинов Р.Ш., Доронина Л.Н. Температурное поле полуограниченного тела в условиях зависимости коэффициента теплообмена и температуры среды от времени // Инж. физ. журн. 1982 Деп. в ВИНИТИ 11.11.81, № 5174-81.

16. Губа В.М., Постольник Ю.С. Расчет температуры и напряжений в пластине при одностороннем нагреве постоянным тепловым122потоком // Изв. вузов. Черная металлургия. 1985. - № 10. - с. 119122.

17. Видин Ю.В. Расчет начальной стадии несимметричного нагрева пластины // Изв. вузов. Черная металлургия. 1985. - № 3. - с. 112115.

18. Скориков К.И., Сабельников А.Г. Инженерный метод расчета температурного поля сплошных и полых цилиндров // Изв. вузов. Черная металлургия. 1983.- № 7. - с. 130-133.

19. Сурков Г.А. К вопросу о нестационарном лучистом взаимодействии твердых тел // Инж. физ. журн. 1980. Т.38. - № 2,- с. 286-289.

20. Точное аналитическое решение одной задачи теплопроводности / В.И. Антипов, В.В. Лебедев, Б.П. Николаев и др. // Инж. физ. журн. -1980. Т.38 № 2. - с. 349. Деп. в ВИНИТИ 02.04.79, № 3926-79 Деп.

21. Аналитическое решение задач теплопроводности для многослойных конструкций при переменных во времени коэффициентах теплообмена / В.А. Кудинов, Н.В. Дилигенский, Н.И. Лаптев и др.// Изв. АН СССР Энергетика и транспорт. 1996. - № 2. - с. 64-68.

22. Решение нестационарных трехмерных задач теплопроводности для многослойных тел / В.А. Кудинов, Б.В. Воробьев, А.Д. Росляков и др. // Изв. АН СССР Энергетика и транспорт. 1991. - № 3. - с. 151 -157.

23. Видин Ю.В. Расчет нагрева излучением термически тонкого тела при переменной температуре среды // Изв. вузов Черная металлургия. -1986. № 3. - с. 117-118.

24. Метод расчета нагрева массивного металла в камерной печи / В.И. Тимошпольский, И.А. Трусова, Ю.А. Малевич и др. // Научн. и прикл. пробл. энерг. 1985. - № 12. - с. 113-116.

25. Полубесов Ю.С., Талдыкин И.А. К расчету температурного поля полого цилиндра // Изв. вузов. Черная металлургия. 1986. - № 1. -с. 138-141.

26. Малевич Ю.А., Тимошпольский В.И., Лигун А.А. К расчету нагрева массивных тел //Изв вузов. Энергетика. 1985. - № 10. - с. 104-107.

27. Сидоров B.C. Решение задачи теплопроводности с граничными условиями третьего рода при произвольной зависимости критерия Bi от координат и времени // Инж. физ. журн. 1981. Т. 41. - № 3. - с. 536-540.

28. Видин Ю.В. Расчет начальной стадии нагрева плоского тела с переменными свойствами // Инж. физ. журн. 1981. Т. 41. - № 3. - с. 480-484.

29. Использование метода конечной глубины проникновения для расчета нагрева плоской плиты под воздействием лучистого теплового потока / В.М. Боришанский, М.А. Готовский, Н.В. Мизонов и др. //Инж. физ. журн. 1980. Т.39. - № 1,- с. 138-142.

30. Панферов В.И. Об оптимальном управлении нагревом окисляющихся массивных тел при теплообмене со средой через поверхностный слой окалины // Изв. вузов. Черная металлургия. -1984.-№2.-с. 87-91.

31. Панферов В.И. Об оптимальном управлении нагревом металла // Изв. вузов. Черная металлургия. 1981. -№ 8. - с. 117-119.

32. Решение краевых задач методом Монте-Карло / Б.С. Еленов, А.А. Кронберг, Г.А. Михайлов, К.К. Сабельфельд. Новосибирск: Наука. Сибирское отделение, 1980. - 176 с.

33. Маликов Ю.К., Лисиенко В.Г. Численный метод решения задач теплопроводности для двумерных тел сложной формы // Инж. физ. журн, 1981. Т. 40. -№ З.-с. 503-509.124

34. Губа В.М. К расчету температуры и напрялсений в слитке прямоугольного сечения при нагреве излучением // Изв. вузов. Черная металлургия. 1998. - № 2. - с. 67-70.

35. Применение схемы Гаусса при численном решении задачи теплопроводности // Изв. вузов. Энергетика. 1988. - № 6. - с. 75-78.

36. Девочкина С.И. Температурное поле шара с переменными теплофизическими свойствами при малых числах Бо. -1985. № 1. -с. 125-128.

37. К расчету температурных полей при нагреве под термообработку слитков в рециркуляционных печах / Ю.А. Герасимов, В.И. Иванов, П.Г. Краснокутский и др. // Изв. вузов. Черная металлургия. 1983. -№ 8,- с.

38. Шайкина А.Е., Клейман Е.И. Решение смешанной краевой линейной задачи теплопроводности в коротком полом цилиндре локально-одномерным методом // Изв. вузов. Черная металлургия. 1989. - № 8.-с. 118-122.

39. Крылова Л. С. К численному методу расчета нестационарного температурного поля полуограниченного тела с переменными коэффициентами // Изв. вузов. Черная металлургия. -1981. № 5. - с.

40. Бардыбахин А.И. Моделирование нагрева окисляющегося металла / Инж. физ. журн. 1985. - Т. 48. - № 4. - 688-689 с. - Деп. в ВИНИТИ 31.10.84, № 7009-84.

41. Панферов В.И., Паразнкин Б.Н. Моделирование нагрева окисляющихся массивных тел методом сеток с подвижными узлами // Изв. вузов. Черная металлургия. 1982. - № 4. - с. 105-109.

42. Юрчук В.Л., Сариогло В.Г. К определению температурных полей в элементах энергетического оборудования // Теплофизические исследования элементов энергетических установок / Сб. н. тр. Киев.: Наукова думка. 1986. - с. 54-58.125

43. Скибин А.П., Червяков В.В., Югов В.П. Метод конечных элементов, основанный на интегрировании по контрольному объему, для двумерных нестационарных эллиптических задач // Изв. вузов. Энергетика. 1995. -№ 1. - с. 142-151.

44. Котенев В.И. Приближенный метод решения задач нестационарной теплопроводности // Изв. АН СССР Энергетика и транспорт. 1989. - № 3. - с. 111-116.

45. Похорилер B.JI. Решение нестационарных задач теплопроводности методом дискетного Z преобразования // Инж. физ. журн. - 1980. Т.39. - № 5.-с. 909-915.

46. Похорилер B.JI. Решение нелинейных задач теплопроводности методом Z преобразования. // Инж. физ. журн. - 1986. - Т. 51. - № 1.- с. 143-150.

47. Похорилер B.JI. Решение двумерных задач теплопроводности методом Z преобразования // Инж. физ. журн. - 1983. - Т. 44. - № 5. -с. 818-824.

48. Коленда З.С., Гнездов E.H. О зональном методе расчета лучистого теплообмена с введением условных поверхностей // Изв. вузов. Черная металлургия. 1982. - № 1.-е. 138-142.

49. Крупенников С.А. Применение модифицированного зонального метода для расчета сложного теплообмена // Изв. вузов. Черная металлургия. 1995. - № 5. - с. 46-49.

50. Суринов Ю.А., Рубцов В.В. Применение обобщенного зонального метода к решению нестационарных задач теории радиационно-кондуктивного теплообмена // Изв. вузов. Энергетика. 1992. - № 7-8.-е. 50-57.

51. Цирельман Н.М. Вариационное решение многомерных задач нестационарной теплопроводности // Теплопроводность и диффузия: Сб. науч. тр. /Риж. политехи, ин-т. Рига, 1987. с. 28-33.

52. Постольник Ю.С. Повышение точности инженерных расчетов радиационного нагрева массивных тел // Изв. вузов. Черная металлургия. 1988. - № 10. - с. 126-129.

53. Слесаренко А.П., Шеметов Н.Ф. Температурное поле пластин и бесконечных призматических тел сложного сечения при переменном во времени коэффициенте теплообмена // Инж. физ. журн. 1983. - Т. 44. - № 4 . - с. 673-676.

54. Бровкин Л.А., Глазов B.C. Температурное поле полупрозрачной пластины при интенсивном нагреве излучением // Изв. вузов. Энергетика. 1987. - № 7. - с. 104-108.

55. Бублевский М.Л., Бублевский Л.И. Модель расчета температурных полей термически массивных тел // Изв. вузов. Черная металлургия. -1991.-№6.-с. 84-86.

56. Борбоц Ю.С. Численно-анолитический метод решения задач теплопроводности с переменными граничными условиями // Изв. вузов. Черная металлургия. 1981. - № 3. - с. 153-156.

57. Губинский В.И., Радченко Ю.Н. Решение нелинейной задачи теплопроводности цилиндра // Изв. вузов. Черная металлургия. -1989. -№ 1.-е. 124-129.

58. Бровкин JI.A., Гузов JI.A. Инженерный расчет нагрева многослойной пластины при граничных условиях первого рода // Изв. вузов. Энергетика. 1985. - № 9. - с. 94-97.

59. Бровкин Л.А. Температурное поле пластины с переменными коэффициентами в установившемся режиме нагрева // Изв. вузов. Энергетика. 1986. - № 5. - с. 72-76.

60. Нагрев полого цилиндра в регулярном режиме / Ю.К. Маликов, A.A. Скворцов, В.Г. Емельянов и др. // Изв. вузов. Черная металлургия. -1986.-№5.-с. 132-134.

61. Ольшанский В.М. Об одном приближенном решении задачи нагрева массивных тел излучением // Изв. вузов. Черная металлургия. 1984. -№ 11.-с. 116-120.

62. Анализ и возможности применения приближенного аналитического решения в расчетах технологии нагрева стали // Изв. вузвов. Энергетика. 1988. - № 4 . - с. 77-81.

63. Тимошпольский В.И. Инженерный способ расчета массивных тел в условиях лучистого теплообмена // Изв. вузов. Черная металлургия. 1986. -№ 7.-с. 126-128.

64. Маликов Ю.К., Скворцов A.A. Нагрев полого цилиндра в регулярном режиме // Изв. вузов. Черная металлургия. 1986. - № 5.

65. Ольшанский В.М. Об одном приближенном решении задачи нагрева массивных тел одновременно излучением и конвекцией // Изв. вузов. Черная металлургия. 1987. - № 1.-е. 121-124.

66. Бровкин Л.А., Гнездов E.H. К обобщению упорядоченного (регулярного) режима нагрева тел простой формы при постоянной температуре среды // Изв. вузов. Черная металлургия. 1981. - № 7. -с. 130-135.

67. Гаврильев Р.И. Температурное поле ограниченного цилиндра при установившемся квазистационарном тепловом режиме / Инж. физ. журн. 1985. Т. 49. 1. - 150-151 с. - Деп. в ВИНИТИ 20.02.85, № 1323-85.

68. Соколов А.К. Численно-аналитический метод расчета несимметричного нагрева пластины с учетом окалинообразования // Изв. вузов. Энергетика. 1994. - № 5-6. - с. 75-80.79. // Изв. вузов. Энергетика. 1986. - № 5. - с. 72-76.

69. Бровкин Л.А., Глазов B.C. Температурное поле полуограниченной пластины при интенсивном нагреве излучением // Изв. вузов. Энергетика. 1987. - № 7. - с. 104-108.

70. Ю.А. Герасимов, В.И. Иванов, П.Г. Краснокутский и др. / к расчету температурных полей при нагреве под термообработку слитков в рециркуляционных печах // Изв. вузов. Черная металлургия. 1983. -№8.-с. 92-95.

71. Денисов М.А., Шкляр Ф.Р. Михалев Г.А. Математическая модель расчета нагрева металла в печах с шагающим подом. Сообщение 1 // Изв. вузов. Черная металлургия. 1980. - № 10.-е. 96-99.

72. Математическая модель расчета нагрева металла в печах с шагающим подом. Сообщение 2 / М.А. Денисов, Г.А. Михалев, Ф.Р. Шкляр, А,А. Кузовников // Изв. вузов. Черная металлургия. 1980. -№ 12.-е. 97-101.129

73. К расчету температурных полей при нагреве под термообработку слитков в рециркуляционных печах / Ю.А. Герасимов, В.И. Иванов, П.Г. Краснокутский и др. // Изв. вузов. Черная металлургия. 1983. -№ 8. - с. 92-95.

74. Решетняк И.С., Словиковский П.А., Судос П.Т. Математическое моделирование нагрева слитков в регенеративном колодце // Интенсификация и повышение эффективности металлургического производства: Межвуз. сб. науч. тр. / Киев, 1980. с. 87-93.

75. Хилков Б.М., Быков В.В. Разработка инженерной методики расчета режимов нагрева металла в методических печах // Промышленная энергетика. -1981.-№4.-с. 40-43.

76. Лаусмаа Т.М., Минк И.Р. Влияние изотермичности кладки и газового объема термической печи на температурное поле садки // Науч. тр. / Таллин, политехи, ин-т. 1981. - № 501. - с. 47-60.

77. Математическая модель нагрева металла в нагревательных колодцах с одной верхней горелкой / Пятков, Венко, Ахманов и др. /7 Металлургия. 1982. Т. 37. - № 4. с. 10-14.

78. Губинский В.И., Серебрянников А.Л. Математическая модель тепловой работы рециркуляционной термической печи / Днепр, металлург, ин-т. Днепропетровск, 1983. - 18 с. - Деп. в Черметинформация 23.05.83, № 2012чм-д83.

79. Антонов В.И. Математическая модель процесса нагрева слитков под ковку в пламенных камерных печах // Изв. вузов. Черная металлургия. 1979. - № 1. - с. 131-137.

80. Антонов В.И. Номограммы для расчета охлаждения садок на выдвинутой подине при предварительной термической обработке // Изв. вузов. Черная металлургия. 1982. - № 9. - с. 135-137.

81. Аронов М.А., Антонов В.Н. Исследование процесса охлаждения сложных садок на выдвинутой подине при предварительной130термической обработке // Конструкции и строительство тепловых агрегатов: Межвуз. сб. науч. тр. / Москва, 1981.-е. 26-32.

82. Андрийко А.Н., Рядно A.A., Каримов И.К. Сопряженная задача радиационно-конвективного теплообмена в камерных печах косвенного нагрева / Днепр, ун-т. Днепропетровск, 1979. - 12 с.-Деп. в ин-те Черметинформация 10.01.80, № 849.

83. Беляев Н.М., Барабаш Н.М., Рядно A.A. Решение сопряженных задач нагрева металла в нагревательных устройствах / Днепр, ун-т. -Днепропетровск, 1979. 8 с. - Деп. в Черметинформация 13.12.79, № 800.

84. Бровкин Л.А., Коленда З.С., Гнездов E.H. К решению сопряженной задачи теплообмена в проходных печах // Изв. вузов. Черная металлургия. 1981. -№ 11.-е. 125-128.

85. Аронов М.А. Исследование режима принудительного охлаждения поковок в термических печах // Кузнечно-штамповочное производжство. 1985. - № 3. - с. 38-40.

86. Губинский В.И., Ефименко Н.И. Математическая модель тепловой работы печи для отжига прутков / Днепр, металлург, ин-т. -Днепропетровск. 1983. - 13 с. - Деп в Укр НИИНТИ 23.06.83, № 586 Ук-Д83.

87. Бровкин Л.А., Крылова Л.С. К определению температурных полей тел с переменными коэффициентами // Изв. вузов. Черная металлургия. 1984. - № 1.-е. 149-153.

88. Маликов Ю.К., Лисиенко В.Г., Волков В.В. Численный метод решения сопряженной задачи радиационно-конвективного и кондуктивного теплообмена // Инж. физ. журн. 1982. Т. 43. - № 3. -с. 467- 474.

89. Крупенников С.А., Решение сопряженной задачи теплообмена в нагревательной печи // Изв. вузов Черная металлургия. 1991. - № 9. -с. 91-93.

90. Математическое моделирование нагрева заготовок в печи с шагающими балками (сопряженная постановка) / В.И. Завелион, Р.Б. Вайс, И.А. Трусова, И.Л. Духвалов // Изв. вузов. Энергетика. 1993. -№ 1-2.-с. 113-117.

91. Крупенников С.А. Решение сопряженной задачи теплообмена в нагревательной печи при наличии нескольких поверхностей сопряжения // Изв. вузов. Черная металлургия. 1994. - № 9. - с. 6165.

92. Зобнин Б.Ф. Нагревательные печи. М:. Машиностроение, 1964. -с. 137.

93. Бровкин Л.А., Коптев Б.Г. Расчетные формулы определения усредненного коэффициента теплоотдачи конвекцией в камерных печах // Изв. Вузов СССР. Черная металлургия. 1980. - № 7. - с. 106107.

94. Тайц Н.Ю. Технология нагрева стали. М:. Металлургиздат, 1962. -с. 567.

95. Губинский В.И., Пашин И.К. Разработка и испытание конструкций высокотемпературных металлических рекуператоров // Сталь.1993,-№6.-с. 91-93.

96. Новгородский Е.Е., Пермяков Б.А. Энергосберегающие установки на машиностроительных заводах // Промышленная энергетика. 1995. -№ 7. - с. 4-6.

97. Выбор оптимального направления использования вторичных энергоресурсов на промышленном предприятии / А.Н. Коротин, В.Ф. Никишов, A.M. Смирнов и др. // Изв. вузов. Черная металлургия.1994. -№3.~ с. 68-71.

98. В.Г. Лисиенко, Ю.К. Шарнин, Ю.К. Маликов и др. Использование струйных рекуператоров для подогрева воздуха в высокотемпературных участках дымоотводящих трактов металлургических агрегатов // Изв. вузов. Черная металлургия. -1992.-№ 2.-с. 66-69.

99. Использование металлических рекуператоров длявысокотемпературного нагрева воздуха в металлургических агрегатах / В.Г. Лисиенко, Ю.К. Маликов, А.В. Саплин, А.Е.133

100. Востротин // Изв. вузов. Черная металлургия. 1992. - № 5. - с. 8991.

101. Казанцев Е.И. Промышленные печи справочник. М.: Металлургия. -1975.-368 с.

102. Михеев М.А., Михеева И.М. Основы теплопередачи. М.: Энергия. -1973.-320 с.

103. Оптимальное по расходу топлива и приведенным затратам напряжение пода печей с шагающими балками / A.A. Кузовников, Е.В. Прополов, Г.П. Куликова и др. // Сталь. 1983. - № 12. - с. 7983.

104. Оптимизация нагрева при обогащении вентиляторного воздуха кислородом и постоянном расходе топлива / М.П. Ревун, А.Н. Минаев, В.Н. Погорелов и др. // Изв. вузов. Черная металлургия. -1985.-№3.-с. 115-118.

105. Оптимизация нагрева при обогащении вентиляторного воздуха кислородом и переменном расходе топлива / М.П. Ревун, А.Н. Минаев, В.Н. Погорелов и др. // Изв. вузов. Черная металлургия. -1985.- №5.- с. 132-136.

106. Девочкина С.И., Захарова Е.В., Кузнецов Г.Г. Управление нагревом металла в камерных печах периодического действия // Изв. вузов. Черная металлургия. 1989. - № 6. - с. 124-126.

107. Ковалевский В.Б., Панасюк В.И., Седяко О.Ю. Алгоритм решения задачи нагрева тел с минимальным расходом топлива // Изв. вузов. Энергетика. 1990. - № 10. - с. 86-90.

108. Программирование нагрева металла по расходу топлива с учетом технологических ограничений / Е.А. Капустин, Л.Э. Гольдфарб, Г.Н. Сидорин и др. // Изв. вузов. Черная металлургия. 1980. - № 7. - с. 114-118.

109. Трубицын Г.В. // Изв. вузов. Черная металлургия. 1982. - № 3. -с. 136-138.

110. Трубицын Г.В., Сотников Г.В. Эффективность оптимальных по расходу топлива режимов нагрева слитков в нагревательных колодцах // Изв. вузов. Черная металлургия. 1983. - № 10. - с. 148149.

111. Трубицын Г.В., Сотников Г.В. Численное решение задачи оптимального нагрева металла в камерных печах с минимальным расходом топлива // Изв. вузов. Черная металлургия. 1983. - № 12. - с. 103-107.

112. Многокритериальная оптимизация режимов нагрева листовых слитков в проходных печах / В.И. Тимошпольский, П.В. Севастьянов, Н.Л. Мандель и др. // Изв. вузов. Энергетика. 1988. -№ 7. - с. 80-84.

113. Повышение качества нагрева стали в пламенных печах / П.Г. Краснокутский, О.И. Тищенко, E.H. Трикашная и др. // Изв. вузов. Черная металлургия. 1991. - № 5. - с. 86-88.

114. Бардыбахин А.И. Наискорейший и экономичный нагрев металла с учетом влияния окисления на теплообмен // Автоматика и телемеханика. 1986. - № 6. - с. 32-41.

115. Применение метода магистральной оптимизации при нагреве термомассивной пластины / Ф.М. Бабушкин, В.Б. Ковалевский, В.И. Тимошпольский и др. // Изв. вузов. Энергетика. 1992. - № 1.-е. 105-107.

116. Бардыбахин А.И. Условия оптимальности для нагрева металла с минимальным окислением // Изв. вузов. Черная металлургия. 1995. - № 3.-е. 65-69.

117. Проректор ИГЭУ по научной работе1. Щелыкалов Ю.Я.199 г.ожидаемого экономического эффекта от внедрения конвективного рекуператора

118. Стоимость операции нагрева при существующем режиме нагрева в термической печи без теплообменника 73,26 руб/т.

119. Стоимость операции нагрева при установке конвективного^4 теплообменника в летний период (по результатам оптимизации на ПЭВМ)-58,33 руб/т.

120. Годовой экономический эффект от установки теплообменника в летний период1. Эм =(£> вгде: пл количество часов работы печей приходящихся на летний период, ч/год;й производительность печи, т/час.

121. Эл = (73,26 58,33) • 1694 • 0,782 = 19778 руб/год

122. Э3 =111309-23641 = 87667 руб/год1. Оглавление

123. Расчет тешюфизических свойств садки, SADKAN 5

124. Задание массивов теплофизических свойствкомпонентов горения, BLOCK DATA BDI 8

125. Задание массивов теплофизических свойств нагреваемых сталей, METAL 10

126. Выбор марки нагреваемой стали, NAMEMET 14

127. Квадратичная интерполяция по формуле Эйткена, RTFX 15

128. Определение коэффициентов полинома видаt(x,Fo)=t(0,0 )+( 1 +b ix+b2x2+ b3x3+b45^KQFPOL 16

129. Решение СЛАУ методомГаусса, FjäÄÜSS 18

130. Расчет режима работы термической печи периодического действия, TERMP3DL 19

131. Л Задание массивов теплофизических свойств,1. BLOCK DATA BDI 90

132. Определение коэффициентов аппроксимации Ь1^ Ьпш RESCUE 93

133. Определение среднемассовой температуры кладки в начальный момент времени, SREDMT 95

134. Расчет вспомогательных функций метода ДУКУ, FUNK 96

135. Интерполяция табличных значений вспомогательных функций метода ДУКУ, ALI 100

136. Интерполяция табличных значений вспомогательных функций метода ДУКУ, ATSE 10137 Решение СЛАУ, SIMQ 102

137. Расчет вспомогательных функций метода ДУКУ, FÜNF 104

138. Расчет вспомогательных функций метода ДУКУ, FUNG 109

139. Расчет вспомогательных функций метода ДУКУ, FUNX 112

140. Расчет конвективного петлевого рекуператора с подвижными перегородками, PETLR

141. Задание массивов теплофизических свойств, BLOCK DATA BD2

142. Задание массивов теплофизических свойств, BLOCK DATA ВОЗ

143. Расчет конвективного петлевого рекуператора в случае, когда все подвижные перегородки подняты, TRUBKA

144. Расчет конвективного петлевого рекуператора в случае, когда часть подвижных перегородок опущена, PEREGR

145. Расчет коэффициента теплоотдачи конвекцией при движении воздуха в 1-ой трубе, TIPR

146. Расчет коэффициента теплоотдачи конвекцией при обтекании продуктами сгорания 1-ой трубы, ALFA

147. Расчет усредненной температуры на внутренней поверхности кладки борова печи, IZLUCH

148. Расчет значений корней характеристического уравнения, XARKOR

149. Расчет стоимости операции нагрева в термической печи, SEBEST

150. Проверка принятых физических ограничений, LIMITS

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

152. Графическое отображение результатов расчета нагрева углеродистой и низколегированной стали, GRAFYGL

153. Подготовка массивов расчета нагрева углеродистой и низколегированной стали, GRAFNVY

154. Графическое отображение результатов расчета нагрева высоколегированной стали, GRAFLEG

155. Подготовка массивов расчета нагрева высоколегированной стали, GRAFNVL5с--------------------------------------------------------------------------------------------------------------------------

156. С РАСЧЕТ ТЕПЛОФИЗИЧЕСКИХ СВОЙСТВ САДКИ. ФАЙЛ ДАННЫХ ЯВЛЯЕТСЯ С ИСХОДНЫМ ДЛЯ РАСЧЕТА ТЕРМИЧЕСКОЙ ПЕЧИ ПЕРИОДИЧЕСКОГО ДЕЙСТВИЯ С ВЫХОДНЫМИ ДАННЫМИ ПРОГРАММЫ ЯВЛЯЮТСЯ: С 1,- ПЛОТНОСТЬ САДКИ (PLEF)

157. С 2. ТЕПЛОЕМКОСТЬ САДКИ (CEF )

158. С 3. ТЕПЛОПРОВОДНОСТЬ САДКИ (LEF )

159. С 4. ТЕМПЕРАТУРОПРОВОДНОСТЬ (AEF )

160. С ВСЕ ПЕРЕЧИСЛЕННЫЕ ЗНАЧЕНИЯ НАХОДЯТСЯ ДЛЯ ИНТЕРВАЛА ТЕМПЕРАТУР С ОТ Т=0 ДО 1000 С ШАГОМ 50.с--------------------------------------------------------------------------------------------------------------------------

161. DIMENSION А( 10),LAMN2( 11 ),LAM02( 11 ),ТАМН20( 11 ),LAMC02( 11),

162. TPLN2(16),TPL02(16),TPLH20(16),TPLC02(16),

163. ТХ2(10),ТХ1(11),ТХ(16),TXC(24),TXL(24),

164. TXRO(24),RON2( 11 ),R002(11 ),R0H20( 10),ROCO2( 11),

165. CMET(24),ROMET(24),LMET(24),PLEF(20),CEF(20),7 ' TR20X( 15),TL20X(5),TC20X( 13),R20X( 15),L20X(5),

166. C20X(13),LEF(20),AEF(20),TG(20)

167. С ИМПОРТ ДАННЫХ И ОБЩЕГО БЛОКА.

168. COMMON /SS/ NAME,K1 ,FI,SIG,S,DT,N2,02,H20,C02

169. COMMON /DDPR/ LAMN2,TAM02,LAMH20,LAMC02,TPLN2,TPL02,TPLH20,

170. TPLC02,TX2,TX1,TX,TXC,TXL,TXR0,R0N2,R002,

171. ROH20,ROC02,TR20X,TL20X,TC20X,R20X,C20X,L20X

172. С ПРИВЕДЕНИЕ К ВЕЩЕСТВЕННОМУ ВИДУ.

173. REAL LAMN2,LAM02,LAMH20,LAMC02,LMET,K1 ,K2,N2,LEF,L20X

174. CHARACTER text*65 EQUIVALENCE (A(1),NAME)

175. С ОТКРЫТИЕ ФАЙЛА ИСХОДНЫХ ДАННЫХ И РЕЗУЛЬТАТОВ.

176. OPEN (1 ,FILE='sadkaef.dat',STATUS='OLD') OPEN (2,FILE='dan.txt',FORM='FORMATTED', 1 ACCESS='DIRECT,,STATUS='OLD,,RECL=75,ERR=999)

177. С ЧТЕНИЕ ДАННЫХ И ФАЙЛА DAN.TXT

178. DO 10 1=1,10 READ (2,5) NP,TEXT,A(I) С WRITE(*,5) NP,TEXT,A(I) 5 FORMAT (I4,A65,F9.4) 10 CONTINUE1. CLOSE (2)1. (NAME.EQ.9) GOTO 30 IF (NAME.NE.0) GOTO 20

179. С ОБРАЩЕНИЕ К ПОДПРОГРАММЕ ДЛЯ ВЫБОРА НОМЕРА СТАЛИ.1. CALL NAMEMET (NAME)1. (NAME.EQ.9) GOTO 30

180. С ОБРАЩЕНИЕ К ПОДПРОГРАММЕ ДЛЯ ПРИСВОЕНИЯ ТЕПЛОФИЗИЧЕСКИХ СВ.

181. CALL METAL (NAME,СМЕТ,LMET,ROMET)

182. ОРГАНИЗАЦИЯ ЦИКЛА ПО ЗНАЧЕНИЯМ ТЕМПЕРАТУРЫ.1. M=01. DO 100 T=50,1000,50 M=M+1

183. CALL RTFX (11 ,TX 1 ,RON2,T,PLN2) CALL RTFX (11 ,TX1 ,R002,T,PL02) CALL RTFX (10,TX2,ROH2O,T,PLH2O) CALL RTFX (11 ,TX 1 ,R0C02,T,PLC02)1. (NAME.EQ.9) THEN

184. CALL RTFX (15,TR20X,R20X,T,PLMET)1. GOTO 401. ELSE1. END IF

185. CALL RTFX (24,TXRO,ROMET,T,PLMET) TG(M)=T

186. PLGAZ=0.01*(N2*PLN2+02*PL02+H20*PLH20+C02*PLC02) PLEF(M)=PLMET*(1 ,-FI)+PLGAZ*FI

187. CALL RTFX (16,TX,TPLN2,T,TCN2) CALL RTFX (16,TX,TPL02,T,TC02) CALL RTFX (16,TX,TPLH20,T,TCH20) CALL RTFX (16,TX,TPLC02,T,TCC02)1. (NAME.EQ.9) THEN

188. CALL RTFX (13 ,TC20X,C20X,T,TPLMET)1. GOTO 501. ELSE1. ENDIF

189. CALL RTFX (24,TXC,CMET,T,TPLMET)

190. TPLGAZ=0.01*(N2*TCN2+02*TC02+H20*TCH20+C02*TCC02) CEF(M)=TPLMET*PLMET*(1.-FI)+TPLGAZ*FI

191. CALL RTFX (11 ,TX1 ,LAMN2,T,RLMN2) CALL RTFX (11 ,TX1 ,LAM02,T,RLM02) CALL RTFX (11,TX1,LAMH20,T,RLMH20) CALL RTFX (11,TX1,LAMC02,T,RLMC02)1. (NAME.EQ.9) THEN

192. CALL RTFX (5,TL20X,L20X,T,RLMET)1. GOTO 601. ELSE1. ENDIF

193. CALL RTFX (24,TXL,LMET,T,RLMET)

194. RLMGAZ=0.01*(N2*RLMN2+02*RLM02+H20*RLMH20+C02*RLMC02)1. (SIG/S.LT.0.001) THEN K2=1.0

195. ALFA=0.16*((T+273.)/100)**3 TERA=(RLMGAZ+ALFA*S)/RLMET GOTO 3001. ELSE ENDIF

196. ALF A=(4.E-08 *((T+273)**4-((T+273 )-DT)* *4))/DT TERAB=(RLMGAZ+ALFA* SIG)/RLMET K2=(S+SIG)*TERAB/(S*TERAB+SIG)81. BLOCK DATA BDI

197. С МАССИВЫ ТЕПЛОФИЗИЧЕСКИХ СВОЙСТВ ПРОДУКТОВ СГОРАНИЯ.

198. DIMENSION LAMN2( 11 ),LAM02( 11),LAMH20(11),LAMC02(11),

199. TPLN2( 16),TPL02(16),TPLH20( 16),TPLC02( 16),

200. TX2( 10),TX 1(11 ),TX(16),TXC(24),TXL(24),

201. TXRO(24),RON2(l 1),R002(11),ROH20(10),

202. R0C02(11),TR20X(15),TL20X(5),TC20X( 13),5 R20X(15),C20X(13),L20X(5)

203. С ПРИВЕДЕНИЕ К ВЕЩЕСТВЕННОМУ ВИДУ.

204. REAL LAMN2,LAM02,LAMH20,LAMC02,L20X С ЗАПИСЬ ДАННЫХ В ОБЩИЙ БЛОК.

205. COMMON /DDPR/ LAMN2,LAM02,LAMH20,LAMC02,TPLN2,

206. TPL02,TPLH20,TPLC02,TX2,TX 1 ,TX,TXC,TXL,

207. TXRO,RON2,R002,ROH20,ROC02,TR20X,TL20X,3 TC20X,R20X,C20X,L20X

208. С ЗАДАНИЕ ЗНАЧЕНИЙ ТЕПЛОФИЗИЧЕСКИХ СВОЙСТВ ПРОДУКТОВ СГОРАНИЯ.

209. DATA (LAMN2(I),I=1,11) /0.0243,0.0315,0.0385,0.0449,1 0.0507,0.0558,0.0604,0.0642,2 0.0672,0.0701,0.0723/

210. DATA (LAM02(I),1=1,11) /0.0247,0.0329,0.0407,0.0480,1 0.0550,0.0615,0.0675,0.0728,2 0.0777,0.0820,0.0858/

211. DATA (LAMH20(I),1=1,11) /0.0148,0.0237,0.0335,0.0442,2 0.0559,0.0684,0.0818,0.0956,3 0.1103,0.1244,0.1407/

212. DATA (LAMC02(I),1=1,11) /0.0147,0.0228,0.0309,0.0391,1 0.0472,0.0549,0.0621,0.0687,2 0.0751,0.0809,0.0863/

213. DATA (TPLN2(I),I=1,16) /1294.0,1298.0,1298.0,1306.0,1 1315.0,1327.0,1340.0,1352.0,2 1365.0,1377.0,1390.0,1403.0,3 1415.0,1424.0,1436.0,1444.0/

214. DATA (TPL02(I),I=1,16) /1306.0,1319.0,1336.0,1357.0,1 1377.0,1398.0,1415.0,1436.0,2 1449.0,1465.0,1478.0,1491.0,3 1503.0,1511.0,1520.0,1528.0/

215. DATA (TPLH20(I),I=1,16) /1495.0,1507.0,1524.0,1541.0,1 1566.0,1591.0,1616.0,1641.0,2 1666.0,1696.0,1725.0,1750.0,3 1775.0,1805.0,1830.0,1855.0/

216. DATA (TPLC02(I),I=1,16) /1599.0,1700.0,1788.0,1863.0,1 1930.0,1989.0,2043.0,2089.0,2 2098.0,2169.0,2202.0,2236.0,3 2265.0,2290.0,2315.0,2336.0/

217. DATA (TX2(I),I=1,10) /100.0,200.0,300.0,400.0,500.0, 2 600.0,700.0,800.0,900.0,1000.0/

218. DATA (TX1(I),1=1,11) /0.0,100.0,200.0,300.0,400.0,2 500.0,600.0,700.0,800.0,900.0,3 1000.0/

219. DATA (TX(I),I=1,16) /0.0,100.0,200.0,300.0,400.0,2 500.0,600.0,700.0,800.0,900.0,3 1000.0,1100.0,1200.0,1300.0,4 1400.0,1500.0/

220. DATA (TXC(I),I=1,24) /100.0,150.0,200.0,250.0,300.0,1 350.0,400.0,450.0,500.0,550.0,2 600.0,650.0,700.0,750.0,800.0,3 850.0,900.0,950.0,1000.0,1050.0,4 1100.0,1150.0,1200.0,1250.0/

221. DATA (TXL(I),I=1,24) /0.0,50.0,100.0,150.0,200.0,250.0,1 300.0,350.0,400.0,450.0,500.0,2 550.0,600.0,650.0,700.0,750.0,3 800.0,850.0,900.0,950.0,1000.0,4 1050.0,1100.0,1150.0/

222. DATA (TXRO(I),I=l ,24) /0.0,15.0,50.0,100.0,150.0,200.0,1 250.0,300.0,350.0,400.0,450.0,2 500.0,550.0,600.0,650.0,700.0,3 750.0,800.0,850.0,900.0,950.0,4 1000.0,1050.0,1100.0/

223. DATA (RON2(I),1=1,11) /1.25,0.916,0.723,0.597,0.508,1 0.442,0.392,0.352,0.318,0.291,2 0.268/

224. DATA (R002(I),1=1,11) /1.429,1.05,0.826,0.682,0.580,1 0.504,0.447,0.402,0.363,0.333,2 0.306/

225. DATA (R0H20(I),I= 1,10) /0.588,0.464,0.384,0.326,0.284, 1 0.252,0.226,0.204,0.187,0.172/

226. DATA (R0C02(I),1=1,11 ) /1.9767,1.447,1.143,0.944,0.802,1 0.698,0.618,0.555,0.502,0.460,2 0.423/

227. DATA (TR20X(I),I=1,15) /20.0,100.0,200.0,300.0,400.0,1 500.0,600.0,650.0,700.0,750.0,2 800.0,850.0,900.0,950.0,1000.0/

228. DATA (TL20X(I),I=1,5) /100.0,200.0,300.0,400.0,500.0/

229. DATA (TC20X(I),I=1,13) /100.0,200.0,300.0,400.0,500.0,1 600.0,650.0,700.0,750.0,800.0,2 850.0,900.0,950.0/

230. DATA (R20X(I),I=1,15) /7830.0,7810.0,7780.0,7780.0,7710.0,1 7710.0,7640.0,7640.0,7640.0,7640.0,2 7640.0,7640.0,7640.0,7640.0,7640.0/

231. DATA (C20X(I),I=1,13) /486.0,561.0,582.0,615.0,653.0,737.0,1 800.0,871.0,1047.0,930.0,1068.0,2 862.0,712.0/

232. DATA (L20X(I),I=1,5) /50.2,46.0,41.9,40.2,38.1/

233. SUBROUTINE METAL (NAME,СМЕТ,LMET,ROMET)

234. С МАССИВЫ ТЕПЛОФИЗИЧЕСКИХ СВОЙСТВ НАГРЕВАЕМЫХ СТАЛЕЙ.

235. DIMENSION C08KP(24),C08(24),C20(24),C40(24),CY8(24),CY8I(24),

236. CY12(24),L08KP(24),L08(24),L20(24),L40(24),LY8(24),

237. LY8I(24),LY12(24),R008KP(24),R008(24),R020(24),

238. RO40(24),ROY8(24),ROY8I(24),ROY 12(24),

239. CSX15(24),LSX15(24),ROSX15(24),

240. CMET(24),ROMET(24),LMET(24)

241. С ПРИВЕДЕНИЕ К ВЕЩЕСТВЕННОМУ ВИДУ.

242. REAL L08KP,L08,L20,L40,LY8,LY8I,LY12,LSX15,LMET С ЗАДАНИЕ ЗНАЧЕНИЙ: ТЕПЛОЕМКОСТЬ,ТЕПЛОПРОВОДНОСТЬ,ПЛОТНОСТЬ.

243. DATA (C08KP(I),I=1,24) /486.0,594.0,502.0,511.0,519.0,1 528.0,536.0,548.0,561.0,574.0,2 590.0,611.0,628.0,645.0,678.0,3 682.0,695.0,691.0,691.0,691.0,4 691.0,687.0,687.0,687.0/

244. DATA (C08(I),I=1,24) /486.0,494.0,502.0,511.0,519.0,1 528.0,540.0,548.0,561.0,578.0,2 590.0,607.0,628.0,662.0,682.0,3 695.0,703.0,609.0,695.0,695.0,4 695.0,691.0,691.0,691.0/

245. DATA (C20(I),I=1,24) /486.0,494.0,502.0,511.0,519.0,1 528.0,540.0,548.0,561.0,574.0,2 590.0,607.0,628.0,669.0,703.0,3 703.0,703.0,699.0,695.0,691.0,4 691.0,691.0,687.0,687.0/

246. DATA (C40(I),I=1,24) /486.0,494.0,498.0,507.0,515.0,1 523.0,532.0,544.0,557.0,569.0,2 582.0,595.0,607.0,678.0,674.0,3 662.0,657.0,653.0,653.0,653.0,4 649.0,649.0,649.0,653.0/

247. DATA (CY8(I),I=1,24) /486.0,502.0,515.0,523.0,532.0,1 540.0,548.0,557.0,589.0,582.0,2 595.0,603.0,619.0,724.0,716.0,4712.0,708.0,703.0,699.0,695.0, 695.0,691.0,691.0,691.0/

248. DATA (CY8I(I),I=1,24) /502.0,511.0,523.0,523.0,519.0,2