Модель «вход-выход» для нагревателя с мешалкой непрерывного действия (CSTH) получается путем применения нестационарного баланса энергии, где температура на входе и тепловая нагрузка определяются как входные параметры, а температура в резервуаре — как выходной. Для резервуара постоянного объема с постоянными физическими свойствами модель имеет вид:
[ \frac{dT}{dt}
\frac{F}{V}(T_i-T) + \frac{Q}{\rho c_p V}, ]
где (T) — температура в хорошо перемешанном резервуаре, (T_i) — температура на входе, (F) — объемный расход, (V) — объем жидкости, а (Q) — чистая скорость теплопередачи. Линеаризация вокруг установившегося состояния дает модель «вход-выход» первого порядка, подходящую для анализа передаточных функций и проектирования обратной связи.
Основная формулировка представляет собой динамический энергетический баланс, выраженный в отклонениях от переменных. При (T_i) и (Q) в качестве входов и (T) в качестве выхода, CSTH постоянного объема ведет себя приблизительно как процесс первого порядка с учетом вкладов температуры на входе и тепловой нагрузки.
Определение задачи управления
Выбор регулируемого выхода
Обычно регулируемой переменной является температура жидкости в резервуаре, (T(t)), поскольку она определяет термическое состояние технологического потока, покидающего нагреватель.
Предполагается, что резервуар идеально перемешивается, поэтому температура на выходе равна температуре в резервуаре:
[ T_{\text{out}}=T. ]
Определение входных параметров процесса
Основным управляющим воздействием является тепловая нагрузка, (Q(t)), подаваемая через пар, электрический нагреватель или нагревательную рубашку.
Температура на входе (T_i(t)) обычно рассматривается как возмущающее воздействие, хотя в некоторых конфигурациях пилотных установок ею также можно управлять. Расход (F(t)) также может выступать в качестве возмущения или управляемой переменной, поскольку он изменяет как время пребывания, так и энергетический поток.
Установление допущений моделирования
Полезная первичная модель обычно предполагает:
- Идеальное перемешивание жидкости
- Постоянный объем жидкости
- Постоянную плотность (\rho)
- Постоянную теплоемкость (c_p)
- Постоянный расход на входе и выходе, (F_i=F)
- Пренебрежимо малые потери тепла, если они не моделируются явно
- Пренебрежимо малые эффекты кинетической, потенциальной энергии, работы вала и работы расширения
Эти допущения позволяют создать компактную модель, сохраняя при этом доминирующую тепловую динамику.
Формулировка динамических балансов
Запись баланса общей массы
Общий баланс объема жидкости имеет вид:
[ \frac{dV}{dt}=F_i-F. ]
Для CSTH с постоянным объемом:
[ F_i=F, \qquad \frac{dV}{dt}=0. ]
Если динамика уровня важна, объем резервуара должен оставаться переменной состояния, а не считаться постоянным.
Запись баланса энергии
Для нереакционного нагревателя с мешалкой нестационарный энергетический баланс выглядит так:
[ \rho c_p V\frac{dT}{dt}
\rho c_p F(T_i-T)+Q. ]
Деление на (\rho c_p V) дает уравнение состояния температуры:
[ \frac{dT}{dt}
\frac{F}{V}(T_i-T) + \frac{Q}{\rho c_p V}. ]
Эквивалентно:
[ \frac{dT}{dt} + \frac{F}{V}T
\frac{F}{V}T_i + \frac{1}{\rho c_p V}Q. ]
Первый член справа представляет конвективную энергию, поступающую с потоком, а (Q) — чистую теплоту, добавленную к жидкости. Охлаждение может быть представлено как (Q<0).
Включение эффектов реакции при необходимости
Если CSTH также функционирует как реактор CSTR, энергетический баланс может включать теплоту реакции:
[ \frac{dT}{dt}
\frac{F_i}{V}(T_i-T) + \frac{(-\Delta H_r)r_A}{\rho c_p} + \frac{Q}{\rho c_p V}, ]
где (r_A) — скорость реакции.
В этом случае температура связана с концентрацией через кинетику реакции. Для реакции первого порядка:
[ r_A=k_0\exp\left(-\frac{E}{RT}\right)c_A, ]
и полная нелинейная модель должна также включать баланс компонентов для (c_A).
Преобразование модели в переменные отклонения
Определение рабочей точки установившегося режима
Пусть установившиеся значения равны (T_s), (T_{i,s}) и (Q_s). В установившемся режиме:
[ \frac{dT}{dt}=0, ]
поэтому рабочая точка удовлетворяет условию:
[ 0= \frac{F}{V}(T_{i,s}-T_s) + \frac{Q_s}{\rho c_p V}. ]
Эта рабочая точка является эталоном, относительно которого анализируются динамические изменения.
Определение отклонений
Определим:
[ T'=T-T_s, ]
[ T_i'=T_i-T_{i,s}, ]
[ Q'=Q-Q_s. ]
Для постоянных (F), (V), (\rho) и (c_p) вычитание уравнения установившегося состояния из динамического уравнения дает:
[ \frac{dT'}{dt} + \frac{F}{V}T'
\frac{F}{V}T_i' + \frac{1}{\rho c_p V}Q'. ]
Используя определение времени пребывания
[ \tau=\frac{V}{F}, ]
модель принимает вид:
[ \boxed{ \frac{dT'}{dt} + \frac{1}{\tau}T'
\frac{1}{\tau}T_i' + \frac{1}{\rho c_p V}Q' } ]
Это стандартная линейная модель «вход-выход» первого порядка для нагревателя с мешалкой непрерывного действия постоянного объема.
Связь модели с общими параметрами
Модель также может быть записана в общем виде:
[ \frac{dT'}{dt}+aT'
\frac{1}{\tau}T_i' + K_QQ', ]
где
[ a=\frac{F}{V}=\frac{1}{\tau}, ]
и
[ K_Q=\frac{1}{\rho c_p V}. ]
Если (Q) обозначает сигнал расхода пара, положение клапана или команду нагревателя, а не фактическую скорость теплопередачи, (K_Q) должен включать соответствующий коэффициент усиления исполнительного механизма и теплопередачи.
Вывод передаточных функций
Применение преобразования Лапласа
Применение преобразования Лапласа при условии нулевого начального отклонения температуры дает:
[ sT'(s)+\frac{1}{\tau}T'(s)
\frac{1}{\tau}T_i'(s) + K_QQ'(s). ]
Следовательно,
[ T'(s)
\frac{1/\tau}{s+1/\tau}T_i'(s) + \frac{K_Q}{s+1/\tau}Q'(s). ]
Получение передаточной функции по температуре на входе
Передаточная функция от возмущения температуры на входе к температуре в резервуаре:
[ \boxed{ \frac{T'(s)}{T_i'(s)}
\frac{1}{\tau s+1} } ]
Таким образом, температура на входе влияет на температуру в резервуаре через динамический отклик первого порядка с постоянной времени (\tau=V/F).
Получение передаточной функции по тепловой нагрузке
Передаточная функция от тепловой нагрузки к температуре в резервуаре:
[ \boxed{ \frac{T'(s)}{Q'(s)}
\frac{K_Q}{s+1/\tau}
\frac{K_T}{\tau s+1} } ]
где
[ K_T=\frac{\tau}{\rho c_p V} =\frac{1}{\rho c_p F}. ]
Таким образом, для прямого ввода тепловой мощности:
[ \boxed{ \frac{T'(s)}{Q'(s)}
\frac{1}{\rho c_p F}\frac{1}{\tau s+1} } ]
Общий температурный отклик — это сумма откликов на возмущение и управляющее воздействие:
[ T'(s)
\frac{1}{\tau s+1}T_i'(s) + \frac{K_T}{\tau s+1}Q'(s). ]
Расширение модели для динамики уровня
Добавление уровня жидкости как состояния
Если уровень в резервуаре не является постоянным, определим (h(t)) как высоту жидкости, а (A) как площадь поперечного сечения резервуара. Баланс массы принимает вид:
[ A\frac{dh}{dt}=F_i-F. ]
Энергетический баланс может быть записан как:
[ Ah\frac{dT}{dt}
F_i(T_i-T) + \frac{Q}{\rho c_p}. ]
Множитель (Ah) — это объем жидкости. Следовательно, изменения уровня изменяют тепловую емкость и, следовательно, динамику температуры.
Решение о необходимости связи с уровнем
Модель только по температуре подходит, когда регулирование уровня поддерживает (V) почти постоянным, а изменения уровня малы.
Связанная модель уровня-температуры требуется, когда изменения расхода приводят к значительным изменениям объема, когда регулирование уровня является частью цели управления или когда время пребывания существенно меняется во время работы.
Понимание смысла управления моделью
Интерпретация постоянной времени
Постоянная времени
[ \tau=\frac{V}{F} ]
— это номинальное время пребывания жидкости. Больший объем резервуара приводит к более медленным изменениям температуры, в то время как больший расход приводит к более быстрому обновлению содержимого резервуара.
После ступенчатого изменения процесс первого порядка достигает примерно 63% своего конечного изменения температуры по истечении одной постоянной времени и асимптотически приближается к конечному значению.
Интерпретация коэффициента усиления процесса
Для изменений тепловой нагрузки коэффициент усиления в установившемся режиме равен:
[ K_T=\frac{1}{\rho c_p F}. ]
Повышение температуры, вызванное заданным подводом тепла, уменьшается при увеличении расхода, плотности или теплоемкости.
Для изменений температуры на входе коэффициент усиления в установившемся режиме равен единице в идеальной модели с постоянными свойствами:
[ \lim_{s\to 0}\frac{T'(s)}{T_i'(s)}=1. ]
В установившемся режиме температура в резервуаре в конечном итоге следует за устойчивым изменением температуры на входе, с учетом влияния регулятора тепловой нагрузки.
Использование модели для проектирования обратной связи
Модель определяет ключевые элементы, необходимые для контура регулирования температуры:
- Регулируемая переменная: температура в резервуаре (T)
- Управляемая переменная: тепловая нагрузка (Q)
- Основное возмущение: температура на входе (T_i)
- Постоянная времени процесса: (\tau=V/F)
- Коэффициент усиления процесса: (K_T=1/(\rho c_p F))
- Измерение: датчик температуры
- Исполнительный элемент: паровой клапан, нагреватель или привод тепловой нагрузки рубашки
Эти параметры поддерживают настройку ПИД-регулятора, анализ подавления возмущений, моделирование и сравнение прогнозируемых и измеренных откликов пилотной установки.
Понимание компромиссов
Не рассматривайте нелинейный реактор как систему первого порядка без проверки
Модель первого порядка справедлива для нереакционной или слабо связанной тепловой системы вблизи фиксированной рабочей точки.
Для экзотермического реактора скорость реакции Аррениуса вносит сильную температурную зависимость. Полученная модель может иметь несколько установившихся состояний, нелинейные коэффициенты усиления или тепловую нестабильность, которые простая передаточная функция первого порядка не может описать.
Не пренебрегайте потерями тепла, когда они важны
Если тепло теряется через стенки резервуара, изоляцию или паровое пространство, может потребоваться дополнительный член:
[ Q_{\text{loss}}=UA(T-T_{\text{amb}}). ]
Тогда энергетический баланс содержит чистый подвод тепла:
[ Q_{\text{net}}=Q-UA(T-T_{\text{amb}}). ]
Игнорирование существенных потерь тепла приводит к ошибкам как в коэффициенте усиления установившегося режима, так и в параметрах настройки регулятора.
Различайте скорость теплопередачи и команду нагревателя
Переменная (Q) должна представлять фактическую скорость теплопередачи, если коэффициент усиления (1/(\rho c_p V)) используется напрямую.
Если измеряемым входом является положение клапана, расход пара, команда электрической мощности или температура рубашки, связь исполнительного механизма и теплопередачи должна быть смоделирована отдельно или включена в экспериментально определенный коэффициент усиления процесса.
Учитывайте динамику измерений и исполнительных механизмов
Реальные пилотные установки часто включают задержку датчика, динамику клапанов, транспортные задержки и динамику конденсации пара.
Если эти эффекты значительны, модель процесса может принимать вид:
[ G(s)= \frac{K_T e^{-\theta s}} {(\tau s+1)(\tau_m s+1)}, ]
где (\theta) — задержка, а (\tau_m) — постоянная времени измерения или исполнительного механизма.
Правильный выбор для вашей цели
Используйте сложность модели, соответствующую вопросу управления и условиям эксплуатации.
- Если ваш основной фокус — базовое регулирование температуры: Используйте модель первого порядка с постоянным объемом, где (T_i) — возмущение, а (Q) — управляемый вход.
- Если ваш основной фокус — регулирование уровня и температуры: Сохраняйте баланс массы и энергии, с (h) или (V) и (T) в качестве переменных состояния.
- Если ваш основной фокус — поведение реактора: Добавьте баланс компонентов и зависящую от температуры кинетику реакции, затем линеаризуйте только вокруг четко определенной рабочей точки.
- Если ваш основной фокус — настройка регулятора пилотной установки: Включите измеренные эффекты исполнительного механизма, датчика, потерь тепла и задержки, если они достаточно велики, чтобы влиять на наблюдаемый отклик.
- Если ваш основной фокус — моделирование или обнаружение неисправностей: Проверьте модель на основе данных ступенчатых тестов перед использованием для обучения операторов, оптимизации или диагностики.
Хорошо сформулированная модель «вход-выход» CSTH связывает законы сохранения с измеримым поведением процесса, предоставляя инженерам обоснованную базу для проектирования управления и безопасной эксплуатации пилотной установки.
Сводная таблица:
| Элемент модели | Определение | Значение для управления |
|---|---|---|
| Регулируемый выход | Температура в резервуаре, T | Определяет термическое состояние выхода |
| Управляемый вход | Тепловая нагрузка, Q | Регулирует температуру процесса |
| Основное возмущение | Температура на входе, Ti | Изменяет тепловую нагрузку |
| Динамический баланс | dT/dt = F/V(Ti - T) + Q/(rho cp V) | Описывает поведение температуры |
| Постоянная времени | tau = V/F | Представляет время пребывания жидкости |
| Коэффициент тепловой нагрузки | KT = 1/(rho cp F) | Связывает подвод тепла с изменением температуры |
| Передаточная функция | T'(s)/Q'(s) = KT/(tau s + 1) | Поддерживает настройку ПИД и проектирование обратной связи |
| Расширения модели | Уровень, реакция, потери тепла, задержка, динамика датчиков | Повышает точность пилотной установки при необходимости |
Планируйте более надежный эксперимент по регулированию температуры с помощью LABPARK. Наши учебно-производственные пилотные установки помогают университетам, исследовательским институтам и предприятиям обучать, тестировать и проверять процессы химической технологии. Выбирайте практические системы для химической инженерии, биопроцессов и биотехнологий или приложений по очистке воды, с реалистичной динамикой процессов для энергетических балансов, анализа передаточных функций, настройки ПИД-регуляторов и обучения операторов. Свяжитесь с LABPARK, чтобы обсудить требуемый масштаб, конфигурацию, приборы и цели управления.
Связанные товары
- Учебный пилотный завод для определения распределения времени пребывания и эффективности перемешивания в каскаде многокаскадных перемешивающих резервуаров
- Многореакторная учебная пилотная установка для изучения технологических процессов реакторной техники
- Многофункциональная учебная пилотная установка для изучения каталитических реакций и оценки реакторов в технологических процессах
- Учебная пилотная установка для определения коэффициента теплопередачи кожухотрубного теплообменника
Люди также спрашивают
- Множественность стационарных состояний CSTR: какие физические переменные соответствуют универсальному развертыванию?
- Как конфигурация мешалки влияет на РВП и работу реактора? Оптимизация перемешивания и конверсии в опытных установках.
- Как последовательное соединение нескольких реакторов CSTR влияет на стабильность? Освойте сложную динамику пилотных установок
- Как модель MRF описывает вращение мешалки в CFD? Основные граничные условия
- Как опытная установка РИС-ПР может продемонстрировать ступенчатую схему реактора? Визуализация экономии объема.