Передаточные функции выводятся путем линеаризации переходного энергетического баланса нагревателя с мешалкой вокруг установившейся рабочей точки. Для идеально перемешиваемого резервуара с постоянной массой жидкости, температурой на входе (T_i), температурой пара (T_{st}), расходом (F), площадью теплообмена (A) и общим коэффициентом теплопередачи (U) полученная модель является моделью первого порядка. Она позволяет получить отдельные передаточные функции от температуры входящей жидкости и температуры пара к температуре в резервуаре.

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

Начните с нестационарного энергетического баланса

Определите допущения для резервуара и потока

Предположим, что резервуар идеально перемешивается, поэтому температура на выходе равна температуре в резервуаре (T). Также предположим, что объем жидкости, плотность, теплоемкость, расход и свойства теплопередачи постоянны.

Резервуар содержит массу жидкости (M), а технологическая жидкость имеет теплоемкость (C_p). Тепло поступает через паровой змеевик в соответствии с уравнением

[ Q = UA(T_{st}-T) ]

где (U) — общий коэффициент теплопередачи, а (A) — эффективная площадь теплообмена.

Запишите член накопления

Тепловая энергия, накопленная в резервуаре, изменяется со скоростью

[ MC_p\frac{dT}{dt} ]

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

Запишите члены потока и нагрева

Чистая энергия, поступающая через технологический поток, составляет

[ FC_p(T_i-T) ]

Первая часть, (FC_pT_i), представляет энергию, поступающую с питанием. Вторая часть, (FC_pT), представляет энергию, покидающую резервуар с идеально перемешанным потоком на выходе.

Добавление тепловой мощности парового змеевика дает полный баланс:

[ MC_p\frac{dT}{dt}

FC_p(T_i-T)+UA(T_{st}-T) ]

Преобразуйте динамическую модель

Раскройте члены теплопередачи

Раскрытие правой части дает

[ MC_p\frac{dT}{dt}

FC_pT_i-FC_pT+UAT_{st}-UAT ]

Сбор членов с температурой резервуара дает

[ MC_p\frac{dT}{dt}+(FC_p+UA)T

FC_pT_i+UAT_{st} ]

Нормализуйте уравнение

Деление на (MC_p) дает

[ \frac{dT}{dt} + \left( \frac{F}{M}+\frac{UA}{MC_p} \right)T

\frac{F}{M}T_i + \frac{UA}{MC_p}T_{st} ]

Определим

[ a=\frac{F}{M}+\frac{UA}{MC_p} ]

[ b_i=\frac{F}{M} ]

[ b_{st}=\frac{UA}{MC_p} ]

Модель принимает вид

[ \frac{dT}{dt}+aT=b_iT_i+b_{st}T_{st} ]

Коэффициент (a) определяет скорость, с которой резервуар приближается к новому равновесию. Его обратная величина,

[ \tau=\frac{1}{a}

\frac{MC_p}{FC_p+UA} ]

является эффективной тепловой постоянной времени.

Перейдите к переменным отклонения

Определите установившуюся рабочую точку

Пусть начальные значения установившегося состояния равны (T_s), (T_{i,s}) и (T_{st,s}). В установившемся состоянии

[ \frac{dT}{dt}=0 ]

поэтому баланс установившегося состояния имеет вид

[ aT_s=b_iT_{i,s}+b_{st}T_{st,s} ]

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

Определите отклонения температуры

Определим переменные отклонения как

[ T'=T-T_s ]

[ T_i'=T_i-T_{i,s} ]

[ T_{st}'=T_{st}-T_{st,s} ]

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

[ \frac{dT'}{dt}+aT'

b_iT_i'+b_{st}T_{st}' ]

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

Примените преобразование Лапласа

Преобразуйте дифференциальное уравнение

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

[ sT'(s)+aT'(s)

b_iT_i'(s)+b_{st}T_{st}'(s) ]

Вынесение преобразования температуры резервуара за скобки дает

[ (s+a)T'(s)

b_iT_i'(s)+b_{st}T_{st}'(s) ]

Следовательно,

[ T'(s)

\frac{b_i}{s+a}T_i'(s) + \frac{b_{st}}{s+a}T_{st}'(s) ]

Определите две передаточные функции

Передаточная функция от температуры входящей жидкости к температуре в резервуаре равна

[ G_i(s)

\frac{T'(s)}{T_i'(s)}

\frac{b_i}{s+a}

\frac{F/M}{s+F/M+UA/(MC_p)} ]

при условии, что отклонение температуры пара равно нулю.

Передаточная функция от температуры пара к температуре в резервуаре равна

[ G_{st}(s)

\frac{T'(s)}{T_{st}'(s)}

\frac{b_{st}}{s+a}

\frac{UA/(MC_p)}{s+F/M+UA/(MC_p)} ]

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

Таким образом, полная зависимость входа от выхода имеет вид

[ T'(s)=G_i(s)T_i'(s)+G_{st}(s)T_{st}'(s) ]

Физическая интерпретация модели

Реакция на температуру на входе

Путь изменения температуры на входе имеет числитель (F/M). Более высокий расход заставляет возмущения на входе сильнее влиять на резервуар и сокращает время пребывания жидкости.

В установившемся состоянии коэффициент усиления по температуре на входе равен

[ K_i=G_i(0)=\frac{b_i}{a} ]

или

[ K_i

\frac{F/M}{F/M+UA/(MC_p)} ]

Реакция на температуру пара

Путь изменения температуры пара имеет числитель (UA/(MC_p)). Увеличение площади теплообмена или общего коэффициента теплопередачи усиливает влияние изменений температуры пара на резервуар.

Его коэффициент усиления в установившемся состоянии равен

[ K_{st}=G_{st}(0)=\frac{b_{st}}{a} ]

или

[ K_{st}

\frac{UA/(MC_p)}{F/M+UA/(MC_p)} ]

Для этой упрощенной модели

[ K_i+K_{st}=1 ]

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

Постоянная времени и полюс

Обе передаточные функции имеют одинаковый знаменатель и, следовательно, одинаковый полюс:

[ s=-a ]

Общая постоянная времени равна

[ \tau=\frac{1}{a} ]

Используя (\tau), передаточные функции также можно записать как

[ G_i(s)=\frac{K_i}{\tau s+1} ]

[ G_{st}(s)=\frac{K_{st}}{\tau s+1} ]

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

Прогнозирование ступенчатых реакций

Ступенчатое изменение температуры на входе

Для ступенчатого изменения величиной (\Delta T_i) температуры на входе, при неизменной температуре пара,

[ T'(t)

K_i\Delta T_i \left(1-e^{-t/\tau}\right) ]

Температура в резервуаре приближается к своему новому установившемуся значению экспоненциально.

Ступенчатое изменение температуры пара

Для ступеньки температуры пара величиной (\Delta T_{st}), при неизменной температуре на входе,

[ T'(t)

K_{st}\Delta T_{st} \left(1-e^{-t/\tau}\right) ]

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

Одновременные изменения

Если обе температуры изменяются одновременно,

[ T'(t)

K_i\Delta T_i \left(1-e^{-t/\tau}\right) + K_{st}\Delta T_{st} \left(1-e^{-t/\tau}\right) ]

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

Связь параметров с измерениями пилотной установки

Оценка теплопередающей способности

Величина (UA) контролирует, насколько эффективно изменения температуры пара передаются в резервуар. Ее можно оценить из

[ Q=UA\Delta T_m ]

где (\Delta T_m) — соответствующая средняя разность температур, часто логарифмическая средняя разность температур для расчетов теплообменников.

Измеренное значение (U) включает сопротивления конвекции, теплопроводности стенки и загрязнений.

Оценка влияния перемешивания

Перемешивание улучшает смешивание и может увеличить коэффициент теплопередачи со стороны жидкости. Оно также подтверждает допущение об идеальном перемешивании, используемое в модели.

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

Осторожное использование безразмерных параметров

Отношение

[ B=\frac{UA}{FC_p} ]

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

Большое значение (B) указывает на то, что паровой нагрев силен по сравнению с эффектами пропускной способности жидкости. Малое значение (B) указывает на то, что входной поток оказывает большее влияние на температуру на выходе.

Понимание компромиссов

Модель первого порядка — это приближение

Модель предполагает постоянные (U), (A), (F), (M) и (C_p). На реальной пилотной установке эти параметры могут варьироваться в зависимости от расхода, температуры, вязкости, загрязнения и скорости мешалки.

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

Температура пара может не равняться тепловой нагрузке

Модель рассматривает температуру пара как входной параметр и использует

[ Q=UA(T_{st}-T) ]

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

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

Динамика транспортировки и датчиков может иметь значение

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

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

Линеаризация ограничивает диапазон применимости

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

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

Сделайте правильный выбор для своей цели

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

• Если ваш основной фокус — возмущения температуры на входе: Используйте (G_i(s)), делая упор на отношение расхода к объему (F/M) и коэффициент усиления ступенчатой реакции на входе (K_i).
• Если ваш основной фокус — управление температурой пара: Используйте (G_{st}(s)), делая упор на (UA/(MC_p)), площадь теплообмена, перемешивание и общий коэффициент (U).
• Если ваш основной фокус — валидация пилотной установки: Выполните ступенчатые испытания, оцените коэффициенты усиления и постоянную времени, а затем сравните измеренную реакцию с прогнозами первого порядка.
• Если ваш основной фокус — высокоточное моделирование: Сохраните нелинейный энергетический баланс и добавьте динамику транспортировки, змеевика, клапана и датчиков там, где данные показывают, что модели первого порядка недостаточно.

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

Сводная таблица:

Моделируйте и проверяйте свою пилотную установку с LABPARK

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

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

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

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

Связанные товары

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

Люди также спрашивают

• Как пилотные установки реакторов определяют пространство проектирования QbD? Мастерство масштабирования химических процессов.
• При настройке пилотной установки с реакторным узлом для высокоэкзотермических реакций, как предотвратить тепловой разгон?
• Как механизмы реакций, катализируемых ионами металлов, влияют на проектирование пилотных установок в химической технологии? От кинетики до масштабирования
• Какую роль играют учебные пилотные установки в масштабировании «зеленой» химии? Преодоление разрыва между лабораторией и промышленностью
• Какие существуют барьеры для интенсификации процессов и как пилотные установки помогают их преодолеть?