close

Вход

Забыли?

вход по аккаунту

;doc

код для вставкиСкачать
 «Вестник ИГЭУ»
Вып. 1 2014 г.
УДК 624.365.3:536.3
Математическая модель оптимизации работы термической печи
для нагрева насыпных садок
О.Б. Колибаба, В.В. Бухмиров, М.Г. Сулейманов
ФГБОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина»,
г. Иваново, Российская Федерация
E-mail: [email protected]
Авторское резюме
Состояние вопроса: В термических печах машиностроительных предприятий распространен нагрев насыпных
садок, образованных изделиями большой номенклатуры. В существующих режимах термообработки в настоящее время применяют далеко не оптимальные скорости нагрева металла. Неучет влияния фильтрации продуктов сгорания через садку приводит к завышенному времени пребывания металла в печи. В связи с этим разработка математической модели термической печи, учитывающей фильтрацию продуктов сгорания через садку,
позволит оптимизировать режимы нагрева насыпных садок.
Материалы и методы: Численная реализация модели осуществлена методом конечных разностей и методом
дискретного удовлетворения краевых условий.
Результаты: Разработана математическая модель сопряженного теплообмена в термической печи для нагрева
насыпных садок, учитывающая фильтрацию продуктов сгорания.
Выводы: Модель позволяет оптимизировать работу печи, используя экономические критерии оптимальности.
Ключевые слова: математическая модель, термическая печь, теплообмен, насыпная садка, оптимальные режимы.
Mathematical model of batch furnace optimization
O.B. Kolibaba, V.V. Bukhmirov, M.G. Suleymanov
Ivanovo State Power Engineering University, Ivanovo, Russian Federation
E-mail: [email protected]
Abstract
Background: It is customary to use the technique of heating batches formed by a great number of products in heating
furnaces of machine-building plants. The speed of metal heating in the thermal treatment conditions currently used is not
high enough. If the influence of burned gas filtration is not taken into account, batch heating takes an unreasonably long
time. Consequently, development of a heating furnace mathematical model taking into account the burned gas filtration
through the batch can improve batch heating conditions.
Materials and methods: The model is implemented by the finite-difference method and the method of discrete satisfaction of boundary conditions.
Results: We have developed a coupled heat transfer mathematical model taking into account burned gas filtration in
batch heating furnaces.
Conclusions: The model can optimize the furnace operation by using economic optimality criteria.
Key words: mathematical model, heating furnaces, heat transfer, batch, optimal conditions.
Термообработка – одна из самых длительных и энергоемких стадий в производстве
металлических изделий, открывающая возможности для решения проблем энерго- и ресурсосбережения на промышленных предприятиях.
Тепловые процессы термической обработки
металла состоят из последовательных циклов,
включающих нагрев металла до заданной температуры, выдержку при постоянной температуре печи и охлаждение с различными скоростями до заданной температуры.
Термообрабатываемые изделия, размещают в термических печах в контейнерах, на
поддонах или насыпью на подине печи, образуя, таким образом, насыпную пористую садку.
Основное требование к нагреву насыпных садок в термических печах – отсутствие местных
перегревов и пережогов при заданных пара-
метрах качества нагрева металла (температуре поверхности изделий и перепаде температур по сечению садки в конце нагрева).
Теория и практика нагрева позволяют
утверждать, что далеко не всегда в существующих тепловых режимах термообработки
реализованы оптимальные скорости нагрева
металла [1]. Необоснованными, как правило,
являются: продолжительность выдержки, общий период нагрева металла и температура
печи при посаде в нее заготовок. Неучет влияния фильтрации продуктов сгорания через
садку приводит к завышению времени пребывания металла в печи, а следовательно, и к
перерасходу энергии.
В связи с этим необходимо исследование возможности сокращения времени пребывания металла в печи (увеличения производи-
 ФГБОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина»
1
 «Вестник ИГЭУ»
Вып. 1 2014 г.
тельности установки) при заданных параметрах качества нагрева металла. Для реализации поставленной цели разработана математическая модель термической печи с учетом
фильтрации продуктов сгорания через садку,
позволяющая находить оптимальные режимы
нагрева насыпных садок.
Математическая модель построена на
базе решения сопряженной задачи теплообмена в системе «газ–кладка–металл», предусматривающего одновременное определение
температурных полей газа, кладки и металла с
учетом их влияния друг на друга.
При моделировании приняты следующие
допущения:
1. Кладка рассмотрена как «идеальная»
для лучистого теплообмена обмуровка. В этом
случае для печи, работающей в режиме, близком к стационарному, затраты теплоты на аккумуляцию кладкой незначительны. Тогда задача сопряженного теплообмена в системе
«газ–кладка–металл» упрощается и сводится к
определению температурного поля металла и
температуры газа, необходимой для поддерживания в рабочем пространстве печи нагрева
металла в целях обеспечения требуемых параметров качества.
2. Излучение газа, кладки и металла
принято «серым». На продолжительность выдержки, а следовательно, на время пребывания металла в печи, механизм поступления
теплоты (излучение или конвекция) не влияет.
Решающая роль в этом случае отводится
внутреннему теплообмену в насыпной садке.
3. Расчетная температура на всей обогреваемой поверхности металла в данный момент времени принята одинаковой. Расчетная
температура на всей внутренней поверхности
кладки также принята одинаковой. Наличие
хорошо развитой турбулизации газов в рабочем пространстве печи позволяет допустить и
изотермичность газового объема.
4. При моделировании насыпная садка
рассмотрена как изотропная пористая неограниченная пластина с эффективными теплофизическими свойствами.
Теплоперенос в садке описан дифференциальным уравнением теплопроводности с
переменными теплофизическими свойствами,
в котором второе слагаемое в правой части
уравнения учитывает фильтрацию газов через
садку:
ρ cмэ
дTм
д  э дTм 
дT
=
λм
+ w f cг м ,
дτ дх 
дх 
дх
(1)
где смэ , λмэ – эффективные удельная теплоемкость и коэффициент теплопроводности садки;
х – координата; τ – время; w – скорость фильтрации продуктов сгорания через садку; f – порозность садки; cг – удельная теплоемкость
газа; Тм – температура металла.
Краевые условия для решения уравнения (1) имеют следующий вид:
– начальное условие
Tм ( х,0) = 0;
– граничные условия
(1, τ)
дT
λмэ м
= σг-м (Tг4 − Tм4 (1, τ)) +
дх
+σкл-м (Tкл4 (1, τ) − Tм4 (1, τ)) + αг-м (Tг − Tм (1, τ)),
λмэ
дTм (0, τ)
= 0,
дx
(2)
(3)
(4)
где Тг – температура газа; σг-м, σкл-м – приведенные коэффициенты излучения в системе
«газ–металл» и «кладка–металл»; αг-м – коэффициент теплоотдачи от газа к металлу.
Теплообмен в кладке описан одномерным дифференциальным уравнением теплопроводности:
скл
дTкл
д 
дT 
=
λкл кл  ,

дτ
дx 
дx 
(5)
где скл , λкл – удельная теплоемкость и коэффициент теплопроводности кладки; Ткл – температура кладки.
Краевые условия для решения уравнения (5) имеют следующий вид:
– начальное условие
Tкл ( x,0) = f ( x );
– граничные условия
дTкл (1, τ)
λкл
= σг-кл [Tг4 − Tкл4 (1, τ)] +
дх
+σм-кл [Tм4 (1, τ) − Tкл4 (1, τ)] + αг-кл [Tг − Tкл (1, τ)],
λкл
дTкл (0, τ)
= αн (Tкл (0, τ) − Tc ),
дх
(6)
(7)
(8)
где Тс – температура окружающей среды; σг-кл ,
σм-кл – приведенные коэффициенты излучения
в системе «газ–кладка», «металл–кладка»;
αг-кл – коэффициент теплоотдачи от газа к
кладке; αн – коэффициент теплоотдачи от
кладки в окружающую среду.
При решении системы дифференциальных уравнений (1), (5) с краевыми условиями
(2)–(4) и (6)–(8) время нагрева является варьируемым параметром, расчетный размер кладки
равен ее толщине, а расчетный размер садки
R при известной емкости печи Е и порозности
садки f рассчитывают по формуле
E
,
(9)
R=
ρ (1 − f ) F
где F – площадь пода, занимаемая садкой; ρ –
плотность металла.
С учетом принятого допущения 1 из системы (1)–(8) можно исключить уравнения
(5)–(8), а уравнение (3) можно записать в виде
q = σгкм (Т г4 − Tм4 ) + αг-м (Tг − Tм ),
(10)
 ФГБОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина»
2
 «Вестник ИГЭУ»
Вып. 1 2014 г.
где q – плотность теплового потока, поступающего на металл; σгкм – приведенный коэффициент излучения газа на металл с учетом излучения кладки.
Расчет приведенного коэффициента излучения σгкм с учетом принятого допущения
выполнен по формуле В.Н. Тимофеева [2]:
σгкм = σ0 ⋅ εм
(1 + ω − εг ) εг
,
1 − εг (εм + (1 − εм ) εг ) + ω
(11)
где σ0 – коэффициент излучения абсолютно
черного тела; εг, εм – степени черноты газа и
металла соответственно; ω – степень развития
обмуровки.
В формуле (11) степень черноты газа εг
зависит от состава газа и его температуры,
которую находят путем решения сопряженной
задачи теплообмена. В нулевом приближении
температуру газа принимают по экспериментальным данным, как средневзвешенную величину за время нагрева.
Коэффициент теплоотдачи от газа к
обогреваемой поверхности металла определен по формуле [3]
αг-м = 0,0558
λг
RV
Кр ,
(12)
где λг – коэффициент теплопроводности газа
(продуктов горения); RV – обобщенный размер
рабочего пространства печи; Кр – безразмерный комплекс [4]:
 Т 
ВVг  г  ρг w г RV3
 273 
,
(13)
Кр =
µг νг Vp.п
где B – расход топлива на печь; Vг – удельный
выход продуктов горения; ρг – плотность продуктов горения; µг – динамический коэффициент вязкости продуктов горения на выходе из
туннеля горелки, определяемый при действительной температуре горения; νг – кинематический коэффициент вязкости продуктов горения в пограничном слое у поверхности металла; wг – скорость выхода продуктов горения из
туннеля горелки; Vр.п – объем рабочего пространства печи.
В уравнении (1) эффективные теплофизические свойства садки аппроксимированы
зависимостями следующего вида:
1


λмэ (T ) = λ 0  1 + bT + b2T 2  ;
(14)
2


cмэ (T ) = С0 Н,
(15)
где Н = (1 + hT
1 ) при T ≤ 1023 К; Н = h2 (1 + h3 T )
при T > 1023 К; b, h1, h2, h3 – коэффициенты
аппроксимации.
При решении сопряженной задачи теплообмена тепловой поток на поверхности металла был аппроксимирован функцией вида

τ 
q k = q0 1 + A1
 ,
τ
k 

(16)
где коэффициент A1 = qk / q0 − 1 ; q0, qk – плотность теплового потока на обогреваемой поверхности в начале и в конце нагрева соответственно.
В начальный момент нагрева плотность
теплового потока рассчитывают по формуле
(10):
4
4
q0 = σгкм (Т г,0
− Т м,0
) + αг-м (Т г,0 − Т м,0 ),
(17)
где Тм,0 – заданная начальная температура
металла; Тг,0 – температура газа в начале нагрева, являющаяся варьируемым параметром
в процессе моделирования.
Величина плотности теплового потока в
конце нагрева определена приближенно по
формуле [4] с последующим уточнением в
процессе решения сопряженной задачи:
э
2∆Tд ⋅ λм
,
(18)
qk =
R
где ∆Тд – заданный допустимый перепад температур по сечению садки в конце нагрева.
По результатам расчета температурного
поля садки находят параметры качества нагрева: температуру поверхности металла в
конце нагрева Тм,к и перепад температур по
сечению садки в конце нагрева ∆Тд. Необходимая температура на поверхности садки в
конце нагрева может быть обеспечена достижением величины средней по сечению конечной температуры металла:
2
(19)
Tср.к = Tм,к − ∆Tд ,
3
Полученное в результате расчета температурного поля садки значение температуры
поверхности садки сравнивают с заданным
значением. Совпадение величин достигают
варьированием плотности теплового потока qk
в конце нагрева.
Допустимый перепад температур по сечению садки в конце нагрева определяют итерационным расчетом, варьируя время нагрева.
При найденных значениях плотности теплового потока на поверхности садки в конце
нагрева, температуры поверхности металла в
конце нагрева и времени нагрева решением
уравнения (10) определяют температуру газа,
которую необходимо поддерживать в рабочем
пространстве печи для обеспечения заданных
параметров качества нагрева металла.
Расход топлива на печь рассчитывают из
уравнения теплового баланса:
Q + Qкл + Qак − QE
B = мр
,
(20)
Qн + QВ + QT − Qух
где QВ – физическая теплота воздуха; QТ – физическая теплота топлива; Qух – потери теплоты с уходящими газами; Qм – затраты теплоты
на нагрев металла; Qкл – потери теплоты через
кладку; Qак – теплота, аккумулированная кладкой; QЕ – теплота экзотермических реакций.
 ФГБОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина»
3
 «Вестник ИГЭУ»
Вып. 1 2014 г.
Для численной реализации математической модели использованы два известных метода: метод конечных разностей с применением чисто неявной схемы [5] для расчета температурного
поля
садки
и
численноаналитический метод дискретного удовлетворения краевых условий (ДУКУ) [6] для решения
сопряженной задачи теплообмена и определения расхода топлива на печь.
С помощью разработанной математической модели были выполнены расчетные исследования процессов нагрева насыпных садок из углеродистых и легированных сталей
(Ст.45, 35Х, 25Х1МФ, Х18Н9ТЛ) в термических
печах, их результаты сопоставлены с данными
промышленного эксперимента (см. рисунок).
Процесс замера температур осуществлялся
непрерывно. Полученное расчетом время нагрева практически совпадает с экспериментальным. Отклонение не превышает 2 %.
Математическая модель тепловой работы садочных термических печей позволяет оптимизировать режимы работы существующих
установок в зависимости от конкретных требований производства.
Заключение
Разработанная математическая модель
сопряженного теплообмена термической печи
для нагрева насыпных садок позволяет находить температуру металла, время пребывания
садки в печи, расход топлива на печь и является базой для выбора оптимального режима
работы печи, исходя из минимальной себестоимости нагрева металла или максимальной
производительности установки.
Список литературы
1. Теплообмен и тепловые режимы в промышленных печах / Ю.И. Розенгарт, Б.Б. Потапов, В.М. Ольшанский, А.В. Бородулин. – Киев: Вищ. шк., 1986. – 290 с.
2. Телегин А.С., Швыдкий В.С., Ярошенко Ю.Г.
Тепломассоперенос. – М.: Академкнига, 2002. – 455 с.
3. Бровкин Л.А., Коптев Б.Г. Расчетные формулы
определения усредненного коэффициента теплоотдачи
конвекцией в камерных печах // Изв. вузов СССР. Черная
металлургия. – 1980. – № 7. – С. 106–107.
4. Бровкин Л.А. Температурные поля тел при нагреве и плавлении в промышленных печах. – Иваново,
1973. – 364 с.
5. Самарский А.А., Вабищевич П.Н. Вычислительная теплофизика. – М.: Либроком, 2009. – 784 с.
6. Крылова Л.С. Проектирование и эксплуатация
теплотехнологических установок кузнечно-термического
производства машиностроительных заводов: учеб. пособие. – Иваново, 2001. – 96 с.
References
Процессов нагрева насыпных садок из стали 25Х1МФ:
–– –– –– – расчетные значения;
– экспериментальные данные; Тг.ад(τ) – температура газа, адаптированная
на показания термопары
Достаточно хорошее совпадение имеют
экспериментальные и расчетные значения
температур поверхности и теплового центра
садки. Расчетные значения температуры газа
отличаются от показаний контрольной термопары в среднем на 5–7 %. Так как в практике
режим нагрева ведется по показаниям контрольной термопары, то полученная моделированием температура газа адаптирована на
условия промышленного эксперимента с помощью функции вида
k ( τ) = κ1τ + κ2τ2 + κ3 τ3 ,
(21)
где k1, k2, k3 – коэффициенты аппроксимации.
1. Rozengart, Yu.I., Potapov, B.B., Ol'shanskiy, V.M.,
Borodulin, A.V. Teploobmen i teplovye rezhimy v promyshlennykh pechakh [Heat transfer and thermal conditions in industrial furnaces]. Kiev, Vishcha shkola, 1986. 290 p.
2. Telegin, A.S., Shvydkiy, V.S., Yaroshenko, Yu.G.
Teplomassoperenos [Heat and mass transfer]. Moscow,
Akademkniga, 2002. 455 p.
3. Brovkin, L.A., Koptev, B.G. Raschetnye formuly
opredeleniya usrednennogo koeffitsienta teplootdachi konvektsiey v kamernykh pechakh [Calculating formulae for definition of the aggregate calorific efficiency in batch furnaces].
Izvestiya vuzov SSSR. Chernaya metallurgiya, 1980, no. 7,
pp. 106–107.
4. Brovkin, L.A. Temperaturnye polya tel pri nagreve i
plavlenii v promyshlennykh pechakh [Temperature fields of
objects being heated and melted in industrial furnaces]. Ivanovo, 1973. 364 p.
5. Samarskiy, A.A., Vabishchevich, P.N. Vychislitel'naya teplofizika [Computational thermal physics]. Moscow,
Librokom, 2009. 784 p.
6. Krylova, L.S. Proektirovanie i ekspluatatsiya teplotekhnologicheskikh ustanovok kuznechno-termicheskogo
proizvodstva mashinostroitel'nykh zavodov [Planning and exploiting of the open-flame reheating furnaces in machinery
construction works]. Ivanovo, 2001. 96 p.
 ФГБОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина»
4
 «Вестник ИГЭУ»
Вып. 1 2014 г.
Колибаба Ольга Борисовна,
ФГБОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина»,
кандидат технических наук, доцент кафедры энергетики теплотехнологий и газоснабжения,
e-mail: [email protected]
Бухмиров Вячеслав Викторович,
ФГБОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина»,
доктор технических наук, профессор, зав. кафедрой теоретических основ теплотехники,
e-mail: [email protected]
Сулейманов Муса Гусейнович,
ФГБОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина»,
аспирант, математик кафедры теоретических основ теплотехники,
e-mail: [email protected]
 ФГБОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина»
5
1/--страниц
Пожаловаться на содержимое документа