Разработку математической модели переходных процессов в наиболее простом варианте начнем с описания процессов, происходящих в приточной вентиляционной системе, обслуживающей некоторое помещение. Теплообменник регулируется по сигналу от датчика температуры воздуха, находящегося в этом помещении и реагирующего на ее отклонение от заданной уставки. Следовательно, контур регулирования является замкнутым.При этом предусматривается качественно-качественный способ регулирования, т.е. колебания теплопоступлений и теплопотерь в помещении устраняются за счет изменения температуры приточного воздуха при постоянном его расходе. В свою очередь, температура притока изменяется вследствие подмешивания того или иного количества охлажденной воды из обратного трубопровода через трехходовой клапан к горячей воде, поступающей в теплообменник, также при постоянном общем ее расходе. На рис. 1 представлена функциональная схема автоматизации приточной системы. Структурная схема САУ для системы обеспечения микроклимата в помещении представлена на рис. 2. При ее построении учтено, что в силу принятого способа регулирования входным параметром для помещения, выбранного в качестве объекта регулирования, являются переменные теплопоступления или теплопотери Q, Вт, а выходным— температура воздуха в помещении tв, °С. Собственно САУ вместе с системой вентиляции в этом случае играют роль отрицательной обратной связи для помещения по каналу «Q–tв». Такая схема была предложена авторами ранее в работе [2]. Тогда математическая модель переходного процесса в САУ будет иметь вид [1]: (1) В соответствии с ранее данным определением и схемой САУ, приведенной на рис. 2, по физическому смыслу Wсист здесь представляетизменение температуры tв с течением времени при единичном тепловом воздействии, т.е. при Q = 1 Вт. Следовательно, размерность Wсист должна быть [К/Вт]. Тогда передаточная функция САУ при использовании П-регулятора в линейном варианте будет выглядеть следующим образом: (2) где a1, …, a7 — коэффициенты, получающиеся при подстановке в (1) передаточных функций элементов САУ с учетом их коэффициентов передачи и постоянных времени, р — некоторый комплексный параметр, имеющий размерность с–1. Выражение (2) представляет переходный процесс в виде изображения, получаемого из переходной функции-оригинала с помощью интегрального преобразования Лапласа. Для того, чтобы по изображению восстановить исходную переходную функцию САУ, необходимо воспользоваться методом обратного преобразования Лапласа. Иначе этот метод называется интегральным преобразованием Карсона [3].Однако в большинстве случаев, кроме самых элементарных, для этого необходимы численные методы, требующие применения ЭВМ и достаточно серьезного программного обеспечения. Поэтому можно попытаться получить значение оригинала переходной функции, используя приближенное аналитическое моделирование. Его сущность заключается в формальной замене оператора p на 1/ι, где ι— время с момента воздействия на систему, с. Такая замена базируется на соображениях, вытекающих из анализа размерностей. Как показывают расчеты на простейших примерах, данный прием позволяет достаточно точно определить характер поведения переходной функции, применяя несложный математический аппарат. Погрешность вычисления максимального отклонения при этом не превышает 15–20%, что вполне достаточно для инженерных расчетов. Используя рассмотренный метод, можно легко решить основную задачу проектирования САУ— подобрать коэффициент передачи регулятора. Для этого нужно вычислить допустимую относительную динамическую ошибку регулирования δ, К/Вт, равную допустимому отклонению температуры ∆tв, деленному на величину теплового воздействия на помещение Q. Сравнивая полученный результат с максимальным значением переходной функции Wmax, найденным с использованием того или иного коэффициента передачи Крег, делаем вывод о пригодности принятого регулятора для рассматриваемой системы, если Wmax ≤ δ . Поскольку в развернутом выражении передаточной функции типа (2) единственным неизвестным является именно Крег, задача имеет однозначное решение. Его легко найти построением приближенного графика оригинала переходной функции, получаемого заменой оператора p на1ι, с помощью, например, электронных таблиц Excel, и добиваясь за счет подбора величины Крег выполнения неравенства Wmax ≤ δ. Однако представляет интерес исследование вопроса о влиянии динамических свойств обслуживаемого помещения на характер переходного процесса и, соответственно, на требуемое значение Крег. На рис. 3 и рис. 4 показаны графики переходных функций для однократного теплового воздействия, полученные расчетом по упрощенному методу для разных Крег в трех вариантах. Черными линиями изображены зависимости для простейшего приближения, когда помещение рассматривается как линейное инерционное звено 1-го порядка. Красные линии характеризуют развитие переходного процесса с учетом так называемой объемной нелинейности, или, точнее говоря, того обстоятельства, что в начальный период времени помещение можно описать как интегрирующее звено с передаточной функцией (формула) здесь V — объем помещения, м3; c, Дж/(кг •К), и ρ, кг/м3 — соответственно удельная теплоемкость и плотность внутреннего воздуха. Синие линии демонстрируют поведение переходной функции для другого приближения, учитывающего поверхностную нелинейность, связанную с процессом распространения температурной волны в толще материальных слоев ограждающих конструкций, обращенных в помещение. В этом случае помещение можно аппроксимировать «полуинтегрирующим» звеном с передаточной функцией вида (формула) потому что скорость начального прогрева ограждений пропорциональна ι0,5 [3]. Здесь 0,886 — значение специальной гамма-функции для аргумента, равного 3/2, а коэффициент Когр определяется через площади поверхностей и их показатели теплоусвоения Y, Вт/(м2 •К), вычисляемые в соответствии с [4]. Сравнивая поведение кривых на рис. 3 и рис. 4, можно обнаружить, что учет нелинейности, связанной с аккумуляцией теплоты в объеме воздуха, заключенного в помещении, не приводит к существенному изменению характера переходной функции. Однако скорость повышения или понижения температуры и ее конечное отклонение при этом увеличивается, причем вне зависимости от коэффициента передачи регулятора. Формально здесь дело обстоит так, как будто бы действующий или эффективный коэффициент передачи регулятора оказывается меньше номинального. Тем не менее, конечное отклонение при расчете с использованием данного варианта нас мало интересует, поскольку при больших интервалах времени с момента теплового воздействия такое приближение будет уже несправедливым. Зато время регулирования в силу тех же обстоятельств при учете объемной теплоаккумуляции сокращается, и мы быстрее достигаем конечных параметров процесса при том же значении Крег. В то же время учет нелинейности, связанной с аккумуляцией теплоты в поверхностных слоях ограждений помещения, приводит к заметному искажению формы получаемого графика, а именно к ускоренному приращению в первые моменты времени, пока температурная волна еще не распространилась глубоко в толщу конструкций, и резкому замедлению в последующий период. Причем при более высоких коэффициентах передачи регулятора (рис. 4) возможно даже некоторое снижение отклонения температуры за счет того, что суммарное поглощение теплоты ограждениями и системой вентиляции становится больше, чем тепловое воздействие на помещение. Конечное значение отклонения при этом также уменьшается по сравнению с линейным вариантом, поэтому учет поверхностной нелинейности дает возможность уменьшить требуемое значение Крег. Теперь нужно еще отметить, что объемная аккумуляция теплоты основную роль играет в начальный момент времени, а поверхностная — в более отдаленные. И тогда можно прийти к выводу, что действительный переходный процесс является промежуточным между вариантами с объемной и поверхностной нелинейностью. Поэтому скорость нарастания отклонения температуры и соответственно требуемое и фактическое время регулирования должны определяться, исходя из предположения, что помещение является интегрирующим звеном, а конечное отклонение и коэффициент передачи регулятора — с учетом распространения температурной волны. Таким образом, мы получили приближенную методику построения переходной функции САУ с учетом нелинейных свойств обслуживаемого помещения как звена контура регулирования. Данная методика имеет достаточно простой вид и при дальнейшем развитии позволяет разработать инженерные рекомендации по подбору регуляторов для систем обеспечения микроклимата, включающие в себя все основные факторы, влияющие на переходные процессы, и доступные для использования в практике массового проектирования.
1.Автоматика и автоматизация систем теплогазоснабжения и вентиляции. / Под ред. В.Н. Богословского.— М.: Стройиздат, 1986. 2.О совершенствовании автоматического регулирования систем обеспечения микроклимата // Журнал «С.О.К.», №5/2005. 3.Беляев Н.М., Рядно А.А. Методы теории теплопроводности. Ч. 2.— М.: Высшая школа, 1982. 4.СП 23-101–2004 «Проектирование тепловой защиты зданий». — М.:ГУП ЦПП, 2004.