Теплоэнергетика, 2021, № 5, стр. 78-88

ПИД-регуляторы с исполнительным механизмом постоянной скорости: варианты реализации, автоподстройка длительности импульсов

В. Ф. Кузищин a, В. С. Дробилко a*, А. Я. Стерлин b**

a Национальный исследовательский университет “Московский энергетический институт”
111250 Москва, Красноказарменная ул., д. 14, Россия

b Центральный аэрогидродинамический институт им. Н.Е. Жуковского
140180 г. Жуковский, Московская обл., ул. Жуковского, д. 1, Россия

* E-mail: dvs@retech.ru
** E-mail: andrey.sterlin@tsagi.ru

Поступила в редакцию 01.07.2020
После доработки 03.09.2020
Принята к публикации 23.09.2020

Полный текст (PDF)

Аннотация

Рассматриваются три варианта реализации регуляторов с исполнительными механизмами (ИМ) постоянной скорости, которые применяются в системах автоматизации технологических процессов: с релейно-импульсным регулирующим блоком (РБИ) с инерционной обратной связью и предвключенным дифференциатором (РД – реальное дифференциирующее звено), с расчетом воздействия по ПИД-алгоритму и применением модели ИМ для его реализации, с расчетом скорости регулирующего воздействия и применением широтно-импульсного модулятора. Приведены формулы расчета параметров первичных элементов блоков РБИ и РД из условия приближения к эквивалентному ПИД-алгоритму. Вариант системы с ПИД-регулятором с управлением ИМ по его модели рассматривается в двух модификациях: с отдельными интеграторами для И-составляющей регулятора и модели ИМ и с общим интегратором для И-составляющей и модели ИМ. В последней модификации сигнал от модели ИМ вычитается из И-составляющей регулятора, что устраняет возможность достижения интегратором предельных значений И-составляющей. Сравнение вариантов производится на основе анализа работы автоматической системы регулирования в пульсирующем режиме работы исполнительного механизма с использованием имитационного моделирования для участка регулирования температуры перегретого пара парового котла с поверхностным пароохладителем. При анализе учитывались как обычные показатели динамической точности, так и число включений ИМ за время переходного процесса. Показано, что достоинством варианта с ПИД-алгоритмом и общим интегратором для И-составляющей и модели ИМ является то, что для него не требуется вычислять вторую производную регулируемой переменной, кроме того, отсутствует проблема накопления разности сигналов от ИМ и его модели. Этот вариант можно считать более предпочтительным при реализации ПИД-алгоритма в программном обеспечении программируемых контроллеров. Даны рекомендации для автоматической подстройки длительности импульсов и параметров релейного элемента регулирующего блока в ходе процесса регулирования в зависимости от отдельных составляющих алгоритма. Приведенные рекомендации обеспечивают существенное снижение частоты включений ИМ в процессе регулирования.

Ключевые слова: автоматический регулятор, ПИД-алгоритм, исполнительный механизм, постоянство скорости, длительность импульсов, объект управления, релейный элемент, широтно-импульсный модулятор

В автоматических системах регулирования (АСР) технологических процессов широко применяются исполнительные механизмы на основе электропривода переменного тока с постоянной скоростью перемещения. Они рассчитаны на управление трехуровневым реверсивным логическим сигналом z(t) = (–1; 0; +1) (здесь z – сигнал управления, t – время) и в процессе работы могут иметь только три расчетных состояния по отношению к скорости S перемещения регулирующего органа (РО) в зависимости от номинальной скорости механизма Sм, % УП/c: “Открытие” (S = Sм), “Стоп” (S = 0) и “Закрытие” (S = –Sм). Номинальную скорость механизма можно выразить через паспортный параметр Tм – время полного хода механизма от 0 до 100% указателя положения (сигнала от датчика положения), с:

(1)
${{S}_{{\text{м}}}} = {\text{ }}{{100} \mathord{\left/ {\vphantom {{100} {{{T}_{{\text{м}}}}}}} \right. \kern-0em} {{{T}_{{\text{м}}}}}}.$

Существуют три основных способа формирования линейных законов регулирования (в рамках ПИД-алгоритма) с ИМ постоянной скорости:

применение релейно-импульсного регулирующего блока с инерционной обратной связью [1, 2], реализующего ПИ-алгоритм, в сочетании с предвключенным дифференциатором для приближения к ПИД-регулятору [3];

использование регулирующего блока, формирующего ПИД-алгоритм с аналоговым выходом, в комплекте с позиционером, который реализует расчетное воздействие ИМ на объект на основе сигнала жесткой обратной связи по положению исполнительного механизма или его модели [47];

применение регулирующего блока, использующего алгоритм приращения и вычисляющего скорость (или приращение) регулирующего воздействия, в сочетании с широтно-импульсным модулятором (ШИМ), формирующим импульсы управления интегрирующим исполнительным механизмом, скважность которых пропорциональна рассчитанной скорости регулирования [8, 9].

При практической реализации и эксплуатации регуляторов с ИМ постоянной скорости необходимо решить следующие основные задачи:

недопущение высокочастотных знакопеременных автоколебаний (переброски) в замкнутом контуре регулирующего блока или позиционера;

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

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

Вопросам реализации регуляторов с применением ИМ постоянной скорости в технической литературе уделяется достаточно большое внимание. Это связано как с широким применением подобных механизмов при автоматизации технологических процессов, так и с внедрением микропроцессорных программируемых контроллеров в системах управления [36].

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

ОБЩИЕ СВЕДЕНИЯ О ВАРИАНТАХ РЕАЛИЗАЦИИ ПИД-АЛГОРИТМА И СИСТЕМЕ РЕГУЛИРОВАНИЯ

Исследования выполнены на имитационной модели АСР для участка регулирования температуры перегретого пара парового котла типа БМ35-РФ (50 т/ч, 4.0 МПа, 440°C). Котел имеет поверхностный пароохладитель, конструкция которого не позволяет получить после него достаточно надежный сигнал по температуре пара. В связи с этим для регулирования температуры перегретого пара применяется одноконтурная АСР. Регулируемой переменной y(t) является температура пара за конвективным пароперегревателем, а регулирующим воздействием μ(t) – изменение степени открытия регулирующего клапана на линии подачи охлаждающей воды через пароохладитель. На основании экспериментальных данных для указанного объекта регулирования была получена передаточная функция

(2)
${{W}_{{{\text{об}}}}}\left( s \right) = \frac{{{{K}_{{{\text{об}}}}}\exp ( - s{{\tau }})}}{{\left( {1 + s{{T}_{1}}} \right)\left( {1 + s{{T}_{2}}} \right)\left( {1 + s{{T}_{3}}} \right)}},$
где Kоб – коэффициент передачи рабочего объекта; s – оператор Лапласа; τ – время запаздывания; Т1, Т2, Т3 – постоянные времени.

В выражении (2):

(3)
$\left. \begin{gathered} {{K}_{{{\text{об}}}}} = {\text{ }}0.{\text{57}}{{^\circ {\text{C}}} \mathord{\left/ {\vphantom {{^\circ {\text{C}}} \% }} \right. \kern-0em} \% };\,\,\,\,{{T}_{{\text{1}}}} = {\text{ 25 с}}; \hfill \\ {{T}_{{\text{2}}}} = {\text{ 25 с}};\,\,\,\,{{T}_{{\text{3}}}} = {\text{ 7}}0{\text{ с}};\,\,\,\,{{\tau }} = {\text{ 14 c}}. \hfill \\ \end{gathered} \right\}$

Для регулятора используется ПИД-алгоритм с передаточной функцией

(4)
${{W}_{{\text{р}}}}\left( s \right) = {{K}_{{\text{п}}}}\left[ {1 + \frac{1}{{{{T}_{{\text{и}}}}s}} + \frac{{{{T}_{{\text{д}}}}s}}{{{{{\left( {{{T}_{{\text{ф}}}}s + 1} \right)}}^{2}}}}} \right],$
где Kп – коэффициент пропорциональности регулятора; Ти, Тд, Тф – постоянные времени интегрирования, дифференцирования и фильтрации; Tд = αTи; Tф = Tд/Cф; α = Tд/Tи; Cф = Tд/Tф.

Здесь при Д-составляющей используется демпфер (фильтр второго порядка с постоянной времени Tф) для подавления случайных помех на высоких частотах. Коэффициент Cф может настраиваться с учетом конкретных условий. По умолчанию принимается Cф = 8, как в контроллерах серии “Ремиконт” [8].

Для объекта (2) с параметрами (3) по методике, приведенной в [10], рассчитаны оптимальные параметры ПИД-регулятора, полученные из условия минимизации квадратичного интегрального показателя ошибки регулирования (I2 = min) при ограничении частотного показателя колебательности М ≤ 1.4 [11]:

(5)
$\left. \begin{gathered} {{K}_{{{\text{п}}{\text{.опт}}}}} = 4.16{\% \mathord{\left/ {\vphantom {\% {^\circ {\text{C}}}}} \right. \kern-0em} {^\circ {\text{C}}}};\,\,\,\,{{T}_{{{\text{и}}{\text{.опт}}}}} = 81\,{\text{с}}; \hfill \\ {{T}_{{{\text{д}}{\text{.опт}}}}} = 28.4\,{\text{с}}\,\,\,\,{\text{при}}\,\,\,\,{{\alpha }_{{{\text{опт}}}}} = {{{{T}_{{{\text{д}}{\text{.опт }}}}}} \mathord{\left/ {\vphantom {{{{T}_{{{\text{д}}{\text{.опт }}}}}} {{{T}_{{{\text{и}}{\text{.опт }}}}}}}} \right. \kern-0em} {{{T}_{{{\text{и}}{\text{.опт }}}}}}} = 0.35. \hfill \\ \end{gathered} \right\}$

В имитационной модели ИМ учитывали ограничения предельных положений и постоянство скорости изменения регулирующего воздействия μ(t). Время полного хода механизма (от 0 до 100%) составляло Тм = 25 с.

При сравнении эффективности вариантов рассматривали реакции системы регулирования при возмущении по каналу регулирующего органа λ = 20% при оптимальных параметрах ПИД-алгоритма. Для полученных процессов оценивали такие показатели, как максимальное отклонение, степень затухания и число включений за фиксированное время переходного процесса.

ВАРИАНТ 1. ПИД-РЕГУЛЯТОР НА БАЗЕ РЕЛЕЙНО-ИМПУЛЬСНОГО БЛОКА С ИНЕРЦИОННОЙ ОБРАТНОЙ СВЯЗЬЮ И ПРЕДВКЛЮЧЕННЫМ ДИФФЕРЕНЦИАТОРОМ

Структурная схема АСР с ПИД-регулятором по варианту 1 показана на рис. 1. Основой ПИД-регулятора является релейно-импульсный регулирующий блок РБИ, содержащий релейный элемент, охваченный инерционной обратной связью. Релейный элемент (РЭ) имеет трехпозиционную характеристику (рис. 2). Здесь и далее

(6)
$\Delta = 2.0{{\Delta }_{{\text{в}}}}.$
Рис. 1.

Структурная схема АСР с ПИД-регулятором по варианту 1. 1 – ПИД-регулирующий блок; 2 – РБИ; 3 – рабочий ИМ; 4 – объект управления; 5 – релейный элемент; 6 – инерционная обратная связь; 7 – дифференциатор; 8 – сумматор; u – сигнал задания; y – регулируемая переменная; yд – переменная на выходе дифференциатора; ε, ε1 – входной сигнал; xо.с – сигнал обратной связи; σ, z – сигналы на входе и выходе релейного элемента; μ – регулирующее воздействие; λ – воздействие на входе в объект управления

Рис. 2.

Характеристика релейного элемента. Δ – зона нечувствительности; Δв – зона возврата; C – уровень выхода РЭ

Звено инерционной обратной связи имеет передаточную функцию апериодического звена:

${{W}_{{{\text{о}}{\text{.с}}}}}\left( s \right) = \frac{{{{K}_{{{\text{о}}{\text{.с}}}}}}}{{{{T}_{{{\text{о}}{\text{.с}}}}}\,s + 1}},$
где Kо.с, То.с – коэффициент передачи и постоянная времени обратной связи.

В отечественной практике регуляторы подобной структуры были разработаны и внедрены в теплоэнергетике еще в 50-х годах ХХ в. и известны как электронные регуляторы ВТИ. В [1] приведены результаты исследования динамических характеристик этих регуляторов при синусоидальном воздействии в трех возможных режимах их работы: пульсирующем, постоянной скорости и одного кратковременного включения за полупериод колебаний. Отмечено, что основной режим работы регулятора – пульсирующий, при котором максимальная расчетная скорость регулирующего воздействия |Sрег(t)| = [dμ/dt] не превосходит располагаемую номинальную скорость исполнительного механизма Sм. В пульсирующем режиме ИМ включается несколько раз подряд в одном направлении, а сигнал обратной связи xо.с(t) успевает отслеживать изменение входного сигнала ε1(t). В этом режиме регулирующий РБИ в комплекте с ИМ постоянной скорости реализует приближенно ПИ-алгоритм регулирования с передаточной функцией

${{W}_{{{\text{ПИ}}}}}\left( s \right) = {{K}_{{\text{п}}}}\left( {1 + \frac{1}{{s{{T}_{{\text{и}}}}}}} \right).$

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

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

Для блока РБИ при работе в пульсирующем режиме длительность импульса tи, длительность паузы tп, период следования импульсов T0 и коэффициент заполнения (скважность) γ можно представить в виде следующих выражений [6]:

$\begin{gathered} {{t}_{{\text{и}}}} = {{{{\Delta }_{{\text{в}}}}} \mathord{\left/ {\vphantom {{{{\Delta }_{{\text{в}}}}} {\left| {{{V}_{{{\text{св}}}}} - \left( {\varepsilon {\text{'}} + {\varepsilon \mathord{\left/ {\vphantom {\varepsilon {{{T}_{{{\text{о}}{\text{.с}}}}}}}} \right. \kern-0em} {{{T}_{{{\text{о}}{\text{.с}}}}}}}} \right)} \right|}}} \right. \kern-0em} {\left| {{{V}_{{{\text{св}}}}} - \left( {\varepsilon {\text{'}} + {\varepsilon \mathord{\left/ {\vphantom {\varepsilon {{{T}_{{{\text{о}}{\text{.с}}}}}}}} \right. \kern-0em} {{{T}_{{{\text{о}}{\text{.с}}}}}}}} \right)} \right|}};\,\,\,\,{{t}_{{\text{п}}}} = {{{{\Delta }_{{\text{в}}}}} \mathord{\left/ {\vphantom {{{{\Delta }_{{\text{в}}}}} {\left| {\left( {\varepsilon {\kern 1pt} {\text{'}} + {\varepsilon \mathord{\left/ {\vphantom {\varepsilon {{{T}_{{{\text{о}}{\text{.с}}}}}}}} \right. \kern-0em} {{{T}_{{{\text{о}}{\text{.с}}}}}}}} \right)} \right|}}} \right. \kern-0em} {\left| {\left( {\varepsilon {\kern 1pt} {\text{'}} + {\varepsilon \mathord{\left/ {\vphantom {\varepsilon {{{T}_{{{\text{о}}{\text{.с}}}}}}}} \right. \kern-0em} {{{T}_{{{\text{о}}{\text{.с}}}}}}}} \right)} \right|}}; \\ {{T}_{0}} = {{{{t}_{{\text{и}}}}{{V}_{{{\text{св}}}}}} \mathord{\left/ {\vphantom {{{{t}_{{\text{и}}}}{{V}_{{{\text{св}}}}}} {\left| {\left( {\varepsilon {\kern 1pt} {\text{'}} + {\varepsilon \mathord{\left/ {\vphantom {\varepsilon {{{T}_{{{\text{о}}{\text{.с}}}}}}}} \right. \kern-0em} {{{T}_{{{\text{о}}{\text{.с}}}}}}}} \right)} \right|}}} \right. \kern-0em} {\left| {\left( {\varepsilon {\kern 1pt} {\text{'}} + {\varepsilon \mathord{\left/ {\vphantom {\varepsilon {{{T}_{{{\text{о}}{\text{.с}}}}}}}} \right. \kern-0em} {{{T}_{{{\text{о}}{\text{.с}}}}}}}} \right)} \right|}}; \\ {{\gamma }} = {{{{t}_{{\text{и}}}}} \mathord{\left/ {\vphantom {{{{t}_{{\text{и}}}}} {{{T}_{0}} = \left| {{{\left( {\varepsilon {\kern 1pt} {\text{'}} + {\varepsilon \mathord{\left/ {\vphantom {\varepsilon {{{T}_{{{\text{о}}{\text{.с}}}}}}}} \right. \kern-0em} {{{T}_{{{\text{о}}{\text{.с}}}}}}}} \right)} \mathord{\left/ {\vphantom {{\left( {\varepsilon {\kern 1pt} {\text{'}} + {\varepsilon \mathord{\left/ {\vphantom {\varepsilon {{{T}_{{{\text{о}}{\text{.с}}}}}}}} \right. \kern-0em} {{{T}_{{{\text{о}}{\text{.с}}}}}}}} \right)} {{{V}_{{{\text{св}}}}}}}} \right. \kern-0em} {{{V}_{{{\text{св}}}}}}}} \right|}}} \right. \kern-0em} {{{T}_{0}} = \left| {{{\left( {\varepsilon {\kern 1pt} {\text{'}} + {\varepsilon \mathord{\left/ {\vphantom {\varepsilon {{{T}_{{{\text{о}}{\text{.с}}}}}}}} \right. \kern-0em} {{{T}_{{{\text{о}}{\text{.с}}}}}}}} \right)} \mathord{\left/ {\vphantom {{\left( {\varepsilon {\kern 1pt} {\text{'}} + {\varepsilon \mathord{\left/ {\vphantom {\varepsilon {{{T}_{{{\text{о}}{\text{.с}}}}}}}} \right. \kern-0em} {{{T}_{{{\text{о}}{\text{.с}}}}}}}} \right)} {{{V}_{{{\text{св}}}}}}}} \right. \kern-0em} {{{V}_{{{\text{св}}}}}}}} \right|}}, \\ \end{gathered} $
где $\varepsilon {\kern 1pt} {\text{'}}$ – производная (скорость изменения) переменной ε; ${{V}_{{{\text{св}}}}} = {{{{K}_{{{\text{о}}{\text{.с}}}}}C} \mathord{\left/ {\vphantom {{{{K}_{{{\text{о}}{\text{.с}}}}}C} {{{T}_{{{\text{о}}{\text{.с}}}}}}}} \right. \kern-0em} {{{T}_{{{\text{о}}{\text{.с}}}}}}}$ – скорость связи (параметр РБИ).

Одним из специфических параметров настройки релейно-импульсного регулятора является tи min – минимальная длительность импульсов, выражение для которой при $\varepsilon {\kern 1pt} {\text{'}} \to 0$ и $\varepsilon \to 0$ записывается в виде

(7)
${{t}_{{{\text{и}}\,{\text{min}}}}} = {{{{\Delta }_{{\text{в}}}}} \mathord{\left/ {\vphantom {{{{\Delta }_{{\text{в}}}}} {{{V}_{{{\text{св}}}}}}}} \right. \kern-0em} {{{V}_{{{\text{св}}}}}}}.$

Из приведенных формул следует, что, изменяя Δв и Δ с учетом (6), можно влиять на длительность и период следования импульсов включения ИМ.

Между параметрами ПИ-алгоритма и параметрами РБИ в пульсирующем режиме существуют следующие приближенные соотношения [6]:

(8)
$\left. \begin{gathered} {{T}_{{{\text{о}}{\text{.с}}}}} = {{T}_{{\text{и}}}};\,\,\,\,{{K}_{{{\text{о}}{\text{.с}}}}} = {{{{T}_{{\text{и}}}}} \mathord{\left/ {\vphantom {{{{T}_{{\text{и}}}}} {\left( {{{T}_{{\text{м}}}}{{K}_{{\text{п}}}}} \right)}}} \right. \kern-0em} {\left( {{{T}_{{\text{м}}}}{{K}_{{\text{п}}}}} \right)}}; \hfill \\ {{K}_{{\text{п}}}} = {{{{S}_{{\text{м}}}}} \mathord{\left/ {\vphantom {{{{S}_{{\text{м}}}}} {{{V}_{{{\text{св}}}}}}}} \right. \kern-0em} {{{V}_{{{\text{св}}}}}}} = {{{{S}_{{\text{м}}}}{{T}_{{{\text{o}}{\text{.c}}}}}} \mathord{\left/ {\vphantom {{{{S}_{{\text{м}}}}{{T}_{{{\text{o}}{\text{.c}}}}}} {\left( {{{K}_{{{\text{о}}{\text{.с}}}}}С} \right)}}} \right. \kern-0em} {\left( {{{K}_{{{\text{о}}{\text{.с}}}}}С} \right)}}. \hfill \\ \end{gathered} \right\}$

Как было отмечено, для повышения качества регулирования может потребоваться применение ПИД-алгоритма при работе с ИМ постоянной скорости [3, 4]. С этой целью в схеме АСР (см. рис. 1) параллельно цепи регулируемой переменной y(t) на входе РБИ включен дифференциатор с передаточной функцией реального дифференцирующего (РД) звена

${{W}_{{{\text{р}}{\text{.д}}}}}(s) = {{{{K}_{{{\text{р}}{\text{.д}}}}}{{T}_{{{\text{р}}{\text{.д}}}}}s} \mathord{\left/ {\vphantom {{{{K}_{{{\text{р}}{\text{.д}}}}}{{T}_{{{\text{р}}{\text{.д}}}}}s} {\left( {{{T}_{{{\text{р}}{\text{.д}}}}}s + 1} \right)}}} \right. \kern-0em} {\left( {{{T}_{{{\text{р}}{\text{.д}}}}}s + 1} \right)}},$
где Kр.д, Тр.д – коэффициент передачи и постоянная времени РД-звена.

На параметры дифференциатора накладываются определенные требования из условия приближения комплексной частотной характеристики (КЧХ) РД-звена ${{W}_{{{\text{р}}{\text{.д}}}}}\left( {j{{\omega }}} \right)$ к КЧХ идеального дифференциатора в рабочем диапазоне частот. Для приближения к идеальному Д-звену достаточно выполнить условие: ${{T}_{{{\text{р}}{\text{.д}}}}} = {{{{T}_{{{\text{д1}}}}}} \mathord{\left/ {\vphantom {{{{T}_{{{\text{д1}}}}}} {{{K}_{{{\text{р}}{\text{.д}}}}}}}} \right. \kern-0em} {{{K}_{{{\text{р}}{\text{.д}}}}}}}$ при ${{K}_{{{\text{р}}{\text{.д}}}}} = 10$ (здесь Тд1 – постоянная времени идеального РД-звена на базе дифференциатора). Тогда

${{W}_{{{\text{р}}{\text{.д}}}}}\left( {j{{\omega }}} \right) \approx {{W}_{{\text{д}}}}\left( {j{{\omega }}} \right) = \left( {j{{\omega }}} \right){{T}_{{{\text{д1}}}}},$
где ${{T}_{{{\text{д1}}}}} = {{T}_{{{\text{р}}{\text{.д}}}}}{{K}_{{{\text{р}}{\text{.д}}}}}.$

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

$\begin{gathered} {{W}_{{{\text{рег}}}}}(s) = {{K}_{{\text{п}}}}\left( {1 + s{{T}_{{{\text{д1}}}}}} \right)\left( {1 + \frac{1}{{s{{T}_{{\text{и}}}}}}} \right) = \\ = {{K}_{{\text{п}}}}\left( {1 + \frac{{{{T}_{{{\text{д1}}}}}}}{{{{T}_{{\text{и}}}}}} + \frac{1}{{s{{T}_{{\text{и}}}}}} + s{{T}_{{{\text{д1}}}}}} \right). \\ \end{gathered} $

Это выражение можно привести к передаточной функции эквивалентного ПИД-алгоритма:

(9)
${{W}_{{{\text{рег}}}}}(s) = {{K}_{{{\text{п}}{\text{.э}}}}}\left( {1 + \frac{1}{{s{{T}_{{{\text{и}}{\text{.э}}}}}}} + s{{T}_{{{\text{д}}{\text{.э}}}}}} \right),$
где

(10)
${{K}_{{{\text{п}}{\text{.э}}}}} = {{K}_{{\text{п}}}}\left( {1 + с} \right);$
(11)
${{T}_{{{\text{и}}{\text{.э}}}}} = {{T}_{{\text{и}}}}\left( {1 + с} \right);$
(12)
${{T}_{{{\text{д}}{\text{.э}}}}} = {{{{T}_{{{\text{д1}}}}}} \mathord{\left/ {\vphantom {{{{T}_{{{\text{д1}}}}}} {\left( {1 + с} \right)}}} \right. \kern-0em} {\left( {1 + с} \right)}};$
(13)
$с = {{{{T}_{{{\text{д1}}}}}} \mathord{\left/ {\vphantom {{{{T}_{{{\text{д1}}}}}} {{{T}_{{\text{и}}}}}}} \right. \kern-0em} {{{T}_{{\text{и}}}}}}.$

Таким образом, структурная схема (см. рис. 1) позволяет реализовать ПИД-регулятор на базе исходного ПИ-регулятора. Однако она имеет следующие недостатки:

параметры эквивалентного ПИД-алгоритма сложным образом связаны с параметрами первичных элементов (ПИ-регулятора и дифференциатора);

соотношение постоянных времени α = Tд/Tи для эквивалентных параметров имеет ограничение по максимальному значению

(14)
${{\alpha }_{{{\text{max}}}}} = 0.25.$

Действительно, из (10)–(13) следует

(15)
$\alpha = {{{{T}_{{{\text{д}}{\text{.э}}}}}} \mathord{\left/ {\vphantom {{{{T}_{{{\text{д}}{\text{.э}}}}}} {{{T}_{{{\text{и}}{\text{.э}}}}}}}} \right. \kern-0em} {{{T}_{{{\text{и}}{\text{.э}}}}}}} = {с \mathord{\left/ {\vphantom {с {{{{\left( {1 + с} \right)}}^{2}}}}} \right. \kern-0em} {{{{\left( {1 + с} \right)}}^{2}}}} = \alpha (с).$

Функция α(с) достигает максимума αmax при с = 1, поэтому получается выражение (14).

Способ расчета параметров первичных элементов (ПИ-регулятора и дифференциатора), если заданы параметры эквивалентного ПИД-алгоритма, выглядит следующим образом. Из (15) можно составить уравнение для расчета параметра c через α (при α ≤ 0.25):

${{c}^{2}}\alpha + c\left( {2\alpha - 1} \right) + \alpha = 0.$

Выражение для расчета с имеет вид

(16)
$c = \frac{{1 - 2\alpha - \sqrt {1 - 4\alpha } }}{{2\alpha }}.$

Далее, используя выражения (10)–(13), можно получить формулы для расчета параметров ПИ-регулятора и дифференциатора, которые обозначены индексом 1:

(17)
${{T}_{{{\text{и1}}}}} = {{{{T}_{{{\text{и}}{\text{.э}}}}}} \mathord{\left/ {\vphantom {{{{T}_{{{\text{и}}{\text{.э}}}}}} {\left( {1 + с} \right)}}} \right. \kern-0em} {\left( {1 + с} \right)}};$
(18)
${{T}_{{{\text{д1}}}}} = {{T}_{{{\text{и1}}}}}c;$
(19)
${{K}_{{{\text{п1}}}}} = {{{{K}_{{{\text{п}}{\text{.э}}}}}} \mathord{\left/ {\vphantom {{{{K}_{{{\text{п}}{\text{.э}}}}}} {\left( {1 + c} \right)}}} \right. \kern-0em} {\left( {1 + c} \right)}}.$

При c = 1, т.е. при α = 0.25, формулы упрощаются:

(20)
${{T}_{{{\text{и1}}}}} = 0.5{{T}_{{{\text{и}}{\text{.э}}}}};$
(21)
${{T}_{{{\text{д1}}}}} = {{T}_{{{\text{и1}}}}};$
(22)
${{K}_{{{\text{п1}}}}} = 0.5{{K}_{{{\text{п}}{\text{.э}}}}}.$

Таким образом, формулы (16)(19) позволяют рассчитать параметры ПИ-регулятора и дифференциатора для рассматриваемого варианта АСР (см. рис. 1) из условия приближения к эквивалентному ПИД-алгоритму при ограничении (14).

Далее приведен пример использования этого варианта для объекта (2) с параметрами (3) и ПИД-алгоритма (4).

Сначала вычисляют параметр α = Tд/Tи эквивалентного ПИД-алгоритма: ${{\alpha }_{{{\text{оп}}}}} = {{28.4} \mathord{\left/ {\vphantom {{28.4} {81}}} \right. \kern-0em} {81}} = 0.35.$ Так как это значение превышает максимально возможное, то принимают α = 0.25. Для этого случая по методике [9] определяют оптимальные параметры: ${{K}_{{{\text{п}}{\text{.опт}}}}}$ = 4.32%/°C, ${{T}_{{{\text{и}}{\text{.опт}}}}}$ = 100 c, ${{T}_{{{\text{д}}{\text{.опт}}}}}$ = 25 c.

По формулам (20)(22) вычисляют параметры ПИ-регулятора и идеального дифференциатора: ${{T}_{{{\text{и1}}}}}$ = 50 с, ${{T}_{{{\text{д1}}}}}$ = 50 с, ${{K}_{{{\text{п1}}}}}$ = 2.16. Далее находят значения параметров блоков 7, 6 (см. рис. 1): ${{K}_{{{\text{р}}{\text{.д}}}}} = 10,$ ${{T}_{{{\text{р}}{\text{.д}}}}} = {{{{T}_{{{\text{д1}}}}}} \mathord{\left/ {\vphantom {{{{T}_{{{\text{д1}}}}}} {10}}} \right. \kern-0em} {10}}$ = 5 с, ${{T}_{{{\text{о}}{\text{.с}}}}} = {{T}_{{{\text{и1}}}}}$ = 50 с, ${{K}_{{{\text{о}}{\text{.с}}}}} = {{{{T}_{{{\text{и1}}}}}} \mathord{\left/ {\vphantom {{{{T}_{{{\text{и1}}}}}} {\left( {{{T}_{{\text{м}}}}{{K}_{{{\text{п1}}}}}} \right)}}} \right. \kern-0em} {\left( {{{T}_{{\text{м}}}}{{K}_{{{\text{п1}}}}}} \right)}}$ = 50/(25 × 2.16) = 0.926. На рис. 3 приведены зависимости основных переменных от времени в процессах в АСР при возмущении λ = 20% по каналу μ(t) для варианта 1.

Рис. 3.

Зависимость основных переменных от времени при процессах в АСР для варианта 1. 1 – λ; 2, 4y; 3, 5 – μ; 2, 3 – без подстройки tи min; 45 – с подстройкой tи min (c1 = 0.25)

Расчеты показали, что полученные графики процессов для АСР с РБИ и ИМ (линии 2, 3) практически совпали с соответствующими графиками для аналогового ПИД-регулятора с теми же параметрами. Однако при этом наблюдается большое число включений ИМ, а именно N = 36 (при времени расчета Tрасч = 600 с и tи min = 0.2 с).

Для сокращения частоты включения ИМ предложен алгоритм автоматической подстройки параметра tи min (с соответствующим изменением Δв и Δ) в ходе процесса регулирования по текущим данным измерения входного сигнала, который действует на основе расчета пропорциональной П- и дифференциальной Д-составляющих в виде следующих выражений:

(23)
$\left. \begin{gathered} {{t}_{{{\text{и}}\,{\text{min}}}}} = {{\left( {\left| \varepsilon \right| + \left| {{{y}_{{\text{д}}}}} \right|} \right){{с}_{1}}{{K}_{{\text{п}}}}{{T}_{{\text{м}}}}} \mathord{\left/ {\vphantom {{\left( {\left| \varepsilon \right| + \left| {{{y}_{{\text{д}}}}} \right|} \right){{с}_{1}}{{K}_{{\text{п}}}}{{T}_{{\text{м}}}}} {100}}} \right. \kern-0em} {100}}; \hfill \\ {{\Delta }_{{\text{в}}}} = {{{{t}_{{{\text{и}}\,{\text{min}}}}}100} \mathord{\left/ {\vphantom {{{{t}_{{{\text{и}}\,{\text{min}}}}}100} {\left( {{{T}_{{\text{м}}}}{{K}_{{\text{п}}}}} \right)}}} \right. \kern-0em} {\left( {{{T}_{{\text{м}}}}{{K}_{{\text{п}}}}} \right)}};\,\,\,\,\Delta = 2.0{{\Delta }_{{\text{в}}}}, \hfill \\ \end{gathered} \right\}$
где с1 – настраиваемый коэффициент (оптимальное значение с1 = 0.25).

При больших отклонениях и скоростях регулирования компоненты ε и yд также имеют большие значения, что увеличивает порог срабатывания и зону возврата РЭ и, следовательно, длительность импульсов и снижает частоту включения. В то же время при приближении к установившемуся состоянию компоненты ε и yд имеют малые значения, что уменьшает порог срабатывания и зону возврата РЭ и длительность импульсов и позволяет реагировать на малые отклонения регулируемой переменной. Рассчитанное значение tи min ограничивается по нижнему пределу исходным заданным значением.

На рис. 3 линии 4, 5 изображены при тех же условиях, что и для предыдущего случая (λ = 20%), но с применением автоматической подстройки параметра tи min. Частота включений ИМ существенно сократилась при незначительном снижении точности регулирования.

Основные показатели качества процессов имеют следующие значения:

без подстройки параметра tи min: максимальное отклонение регулируемой переменной A1 = 3.89, степень затухания ψ = 0.97, число включений ИМ в интервале времени расчета (Трасч = 600 с) N = 36;

с подстройкой параметра tи min (при c1 = 0.25): A1 = 4.01, ψ = 0.93, N = 18.

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

ВАРИАНТ 2. РЕГУЛИРУЮЩИЙ БЛОК С ПИД-АЛГОРИТМОМ И ФОРМИРОВАНИЕМ ИМПУЛЬСОВ УПРАВЛЕНИЯ ПО МОДЕЛИ ИСПОЛНИТЕЛЬНОГО МЕХАНИЗМА

Структурная схема АСР с регулятором по варианту 2 показана на рис. 4.

Рис. 4.

Структурная схема АСР с регулятором по варианту 2 с отдельными интеграторами (а) и общим интегратором (б) для И-составляющей и модели ИМ. 1 – блок ПИД; 2 – позиционер; 3 – рабочий ИМ; 4 – объект управления; 5, 6, 7 – П-, Д-, И-компоненты; 8 РЭ; 9 – модель ИМ; 10 – интегратор; 11, 12 – узел масштабирования; 13 – комплексный регулирующий блок; yп.д – суммарный сигнал П- и Д- компонент; xинт, yинт – сигналы на входе и выходе интегратора; zб, zм – сигналы “Больше” и “Меньше” на выходе РЭ

В блоке 1 формируется полностью расчетное значение ПИД-алгоритма регулирования yр = yп + + yд + yи (здесь yп, yд, yи – пропорциональная, дифференциальная и интегральная составляющие). В блоке 2 формируются импульсный сигнал z(t) для управления исполнительным механизмом 3 для реализации расчетного воздействия yр. Этот блок имеет схему позиционера и содержит релейный элемент и модель идеального ИМ, которые охвачены жесткой обратной связью. Модель ИМ выполнена в виде интегратора с передаточной функцией ${{W}_{{{\text{инт}}}}}(s) = {1 \mathord{\left/ {\vphantom {1 {(s{{T}_{{{\text{инт}}}}})}}} \right. \kern-0em} {(s{{T}_{{{\text{инт}}}}})}}$ при Tинт = Tм (здесь Тинт – постоянная времени интегратора модели ИМ).

В процессе работы АСР выходной сигнал модели ИМ от позиционера yм следит за сигналом yр на выходе регулирующего блока с точностью до Δв = Δ/2. При каждом срабатывании РЭ включается рабочий ИМ, действующий на объект 4. При Tинт = Tм приращения регулирующего воздействия Δμ приближенно повторяют приращения выхода позиционера Δyм и регулирующего блока Δyр. Формирование сигнала yр на выходе блока 1 в операторной форме описывается выражением, соответствующим передаточной функции (4) ПИД-регулятора:

(24)
$\begin{gathered} {{y}_{{\text{р}}}} = {{y}_{{\text{п}}}} + {{y}_{{\text{и}}}} + {{y}_{{\text{д}}}} = {{K}_{{\text{п}}}}\left( {u - y} \right) + \\ + \,\,\frac{{{{K}_{{\text{п}}}}\left( {u - y} \right)}}{{{{T}_{{\text{и}}}}s}} + \frac{{{{K}_{{\text{п}}}}{{T}_{{\text{д}}}}sy}}{{{{{\left( {{{T}_{{\text{ф}}}}s + 1} \right)}}^{2}}}}. \\ \end{gathered} $

Минимальная длительность импульса включения tи min (при медленных изменениях сигнала yр) для данного варианта регулятора определяется формулой

(25)
${{t}_{{{\text{и}}\,{\text{min}}}}} = {{{{\Delta }_{{\text{в}}}}} \mathord{\left/ {\vphantom {{{{\Delta }_{{\text{в}}}}} {{{S}_{{\text{м}}}}}}} \right. \kern-0em} {{{S}_{{\text{м}}}}}} = {{{{\Delta }_{{\text{в}}}}{{T}_{{\text{м}}}}} \mathord{\left/ {\vphantom {{{{\Delta }_{{\text{в}}}}{{T}_{{\text{м}}}}} {100}}} \right. \kern-0em} {100}}.$

Так как обычно значение tи min является заданным, то для его обеспечения требуется устанавливать зону возврата по формуле

(26)
${{\Delta }_{{\text{в}}}} = {{{{t}_{{{\text{и}}\,{\text{min}}}}}100} \mathord{\left/ {\vphantom {{{{t}_{{{\text{и}}\,{\text{min}}}}}100} {{{T}_{{\text{м}}}}}}} \right. \kern-0em} {{{T}_{{\text{м}}}}}}.$

Достоинство схемы с позиционером (см. рис. 4, а) – возможность применения регулирующего блока с аналоговым выходным сигналом. Следует отметить, что в процессе эксплуатации рассмотренного варианта АСР может постепенно накапливаться разница между выходами интегрирующих устройств 3 и 9 из-за неизбежного различия характеристик рабочего ИМ и его модели, что может нарушить работоспособность системы.

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

Блоки 8, 10 и 12 на рис. 4, б образуют замкнутый контур, выполняющий те же функции, что и позиционер 2 на рис. 4, а. Подобный принцип работы используется в функциональных блоках PID_ES программного обеспечения контроллеров Simatic S7 фирмы Siemens.

Формирование регулирующего воздействия μ(t) в комплекте с ИМ для рассматриваемого варианта в операторной форме описывается следующим выражением, соответствующим ПИД-алгоритму регулирования:

(27)
${{\mu }} = {{K}_{{\text{п}}}}\left( {u - y} \right) + \frac{{{{K}_{{\text{п}}}}\left( {u - y} \right)}}{{{{T}_{{\text{и}}}}s}} + \frac{{{{K}_{{\text{п}}}}{{T}_{{\text{д}}}}sy}}{{{{{\left( {{{T}_{{\text{ф}}}}s + 1} \right)}}^{2}}}} \cdot $

Формула (27) по своему виду повторяет приведенную ранее формулу (24), но она выражена относительно выходной переменной позиционера μ(t), что отвечает назначению позиционера – формировать воздействие на объект в соответствии с расчетным сигналом yр(t) на выходе алгоритма регулирования.

Следует отметить, что здесь Д-составляющая формируется только для сигнала регулируемой переменной y(t) и не применяется для сигнала задания u(t). Это позволяет исключить броски выходного сигнала при изменении задания. Если требуется оптимизировать процесс при изменении задания, то для этого случая рекомендуется использовать специальный командный блок [11].

При автоматической подстройке длительности импульсов по алгоритму блока PID_ES предусмотрена адаптация порогов срабатывания РЭ (Δ и Δв) в ходе процессов регулирования для уменьшения частоты включения ИМ. Это производится в зависимости от сигналов yп. д и yинт (см. рис. 4) по алгоритму

(28)
${{\Delta }_{{\text{в}}}} = \min \left\{ {\left| {{{y}_{{{\text{п}}{\text{.д}}}}}} \right|;\left| {{{y}_{{{\text{инт}}}}}} \right|} \right\}{{c}_{2}},$
где с2 – настраиваемый коэффициент.

В целях уменьшения частоты включения для варианта 2 был применен предложенный выше альтернативный способ автоматической подстройки длительности импульсов путем изменения порога срабатывания РЭ (при Δ = 2Δв) в зависимости от сигналов yп и yд аналогично выражению (23). Алгоритм подстройки имеет вид

(29)
$\left. \begin{gathered} {{t}_{{{\text{и}}\,{\text{min}}}}} = {{\left( {\left| {{{y}_{{\text{п}}}}} \right| + \left| {{{y}_{{\text{д}}}}} \right|} \right){{с}_{1}}{{T}_{{\text{м}}}}} \mathord{\left/ {\vphantom {{\left( {\left| {{{y}_{{\text{п}}}}} \right| + \left| {{{y}_{{\text{д}}}}} \right|} \right){{с}_{1}}{{T}_{{\text{м}}}}} {100}}} \right. \kern-0em} {100}}; \hfill \\ {{\Delta }_{{\text{в}}}} = {{{{t}_{{{\text{и}}\,{\text{min}}}}}100} \mathord{\left/ {\vphantom {{{{t}_{{{\text{и}}\,{\text{min}}}}}100} {{{T}_{{\text{м}}}}}}} \right. \kern-0em} {{{T}_{{\text{м}}}}}};\,\,\,\,\Delta = 2.0{{\Delta }_{{\text{в}}}}. \hfill \\ \end{gathered} \right\}$

На рис. 5 в качестве примера приведены графики процессов регулирования для варианта 2 при тех же условиях, что и для варианта 1.

Рис. 5.

Зависимость основных переменных от времени при процессах в АСР для варианта 2. 1 – λ; 2, 4y; 3, 5 – μ; 2, 3 – без подстройки параметров; 4, 5 – с подстройкой параметров по алгоритму (28) (а), по алгоритму (29) (б)

Основные показатели качества процессов при tи min = 0.2 c, Tрас = 600 с имеют следующие значения:

без подстройки параметров Δв, Δ, tи minA1 = = 3.74°С, ψ = 0.99, N = 44;

с подстройкой параметров Δв, Δ по алгоритму (28) A1 = 4.02°С, ψ = 0.99, N = 20;

с подстройкой параметра tи min по алгоритму (29) A1 = 3.89°С, ψ = 0.96, N = 19.

Анализ данных качества процессов показал, что подстройка Δв и Δ по алгоритму (28) позволяет существенно снизить число включений ИМ. При этом немного возрастает динамическая ошибка регулирования. Автоматическая подстройка параметра tи min по предлагаемому алгоритму (29) позволила снизить число включений ИМ даже в большей степени, чем по алгоритму (28), предусмотренному в функциональном блоке PID_ES. Динамическая ошибка регулирования оказалась здесь также меньше.

Следует отметить, что система по варианту 2 имеет максимальное отклонение чуть меньше, чем по варианту 1. Это связано с тем, что для варианта 1 параметр α = Tд/Tи был ограничен значением 0.25 (при оптимальном значении αоп = 0.31), а для варианта 2 такое ограничение отсутствует.

ВАРИАНТ 3. РЕГУЛИРУЮЩИЙ БЛОК С РАСЧЕТОМ СКОРОСТИ РЕГУЛИРУЮЩЕГО ВОЗДЕЙСТВИЯ И ШИРОТНО-ИМПУЛЬСНЫМ МОДУЛЯТОРОМ

Структурная схема АСР с ПИД-регулятором по варианту 3 приведена на рис. 6. В схеме используется регулирующий блок 1 с импульсным выходом (блок РИМ) с расчетом приращения Δyр и коэффициентом заполнения γ импульсного сигнала в комплекте с широтно-импульсным модулятором 2. Для управления средней скоростью регулирующего воздействия с ИМ постоянной скорости используется трехуровневый импульсный сигнал z(t), который создается с помощью блока ШИМ под действием аналогового сигнала yр(t) с выхода регулирующего блока 1. Импульсный сигнал характеризуется коэффициентом заполнения

(30)
${{\gamma }} = {{{{t}_{{\text{и}}}}} \mathord{\left/ {\vphantom {{{{t}_{{\text{и}}}}} {\left( {{{t}_{{\text{и}}}} + {{t}_{{\text{п}}}}} \right)}}} \right. \kern-0em} {\left( {{{t}_{{\text{и}}}} + {{t}_{{\text{п}}}}} \right)}} = \frac{1}{{1 + {{{{t}_{{\text{п}}}}} \mathord{\left/ {\vphantom {{{{t}_{{\text{п}}}}} {{{t}_{{\text{и}}}}}}} \right. \kern-0em} {{{t}_{{\text{и}}}}}}}}.$
Рис. 6.

Структурная схема АСР с расчетом приращения воздействия и широтно-импульсным модулятором. 1 – регулирующий блок; 2 – ШИМ; 3 –рабочий ИМ; 4 – объект управления; 57 – блоки вычисления приращений П-, Д-, И-составляющих; 8 – блок вычисления γр

При подаче аналогового сигнала x (от регулирующего блока) формирование импульсного сигнала на выходе ШИМ [9] происходит в соответствии с рис. 7.

Рис. 7.

Схема блока широтно-импульсной модуляции (а), графические зависимости, поясняющие работу ШИМ (б). 1 – интегратор; 2 – релейный элемент; e – сигнал на входе релейного элемента блока ШИМ

Длительности паузы и импульса вычисляются по формулам

(31)
${{t}_{{\text{п}}}} = {{{{\Delta }_{{\text{в}}}}{\text{d}}t} \mathord{\left/ {\vphantom {{{{\Delta }_{{\text{в}}}}{\text{d}}t} {\left| x \right|}}} \right. \kern-0em} {\left| x \right|}};$
(32)
${{t}_{{\text{и}}}} = {{{{\Delta }_{{\text{в}}}}{\text{d}}t} \mathord{\left/ {\vphantom {{{{\Delta }_{{\text{в}}}}{\text{d}}t} {\left( {1 - \left| x \right|} \right)}}} \right. \kern-0em} {\left( {1 - \left| x \right|} \right)}},$
где dt – время цикла расчета.

Минимальная длительность импульса ${{t}_{{{\text{и}}\,{\text{min}}}}}$ (при x << 1) определяется выражением

(33)
${{t}_{{{\text{и}}\,{\text{min}}}}} = {{\Delta }_{{\text{в}}}}{\text{d}}t.$

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

(34)
${{\Delta }_{{\text{в}}}} = {{{{t}_{{{\text{и}}\,{\text{min}}}}}} \mathord{\left/ {\vphantom {{{{t}_{{{\text{и}}\,{\text{min}}}}}} {{\text{d}}t}}} \right. \kern-0em} {{\text{d}}t}}.$

Значение коэффициента заполнения импульсного сигнала зависит от отношения ${{{{t}_{{\text{п}}}}} \mathord{\left/ {\vphantom {{{{t}_{{\text{п}}}}} {{{t}_{{\text{и}}}}}}} \right. \kern-0em} {{{t}_{{\text{и}}}}}},$ выражение для которого в соответствии с (31), (32) выглядит так:

${{{{t}_{{\text{п}}}}} \mathord{\left/ {\vphantom {{{{t}_{{\text{п}}}}} {{{t}_{{\text{и}}}}}}} \right. \kern-0em} {{{t}_{{\text{и}}}}}} = {{\left( {1 - x} \right)} \mathord{\left/ {\vphantom {{\left( {1 - x} \right)} x}} \right. \kern-0em} x} = \frac{1}{x} - 1.$

С учетом формулы (30) γ = yр.

Таким образом, приведенная структура ШИМ обеспечивает формирование импульсного сигнала, коэффициент заполнения которого равен входному сигналу от регулирующего блока. С учетом формул (31)–(33) можно получить выражения для относительной длительности импульсов

${{{{t}_{{\text{и}}}}} \mathord{\left/ {\vphantom {{{{t}_{{\text{и}}}}} {{{t}_{{{\text{и}}\,{\text{min}}}}}}}} \right. \kern-0em} {{{t}_{{{\text{и}}\,{\text{min}}}}}}} = {1 \mathord{\left/ {\vphantom {1 {\left( {1 - {{\gamma }}} \right)}}} \right. \kern-0em} {\left( {1 - {{\gamma }}} \right)}}$

и периода их следования

${{{{T}_{{{\text{сл}}}}}} \mathord{\left/ {\vphantom {{{{T}_{{{\text{сл}}}}}} {{{t}_{{{\text{и}}\,{\text{min}}}}}}}} \right. \kern-0em} {{{t}_{{{\text{и}}\,{\text{min}}}}}}} = {1 \mathord{\left/ {\vphantom {1 {\left[ {\left( {1 - {{\gamma }}} \right){{\gamma }}} \right]}}} \right. \kern-0em} {\left[ {\left( {1 - {{\gamma }}} \right){{\gamma }}} \right]}}.$

Эта функция имеет минимум при γ = 0.5. При этом $0.5{{T}_{{{\text{сл}}}}} = 4{{t}_{{{\text{и}}\,{\text{min}}}}},$ что может привести к повышенной частоте включения ИМ.

Усредненное значение скорости регулирующего воздействия Sрег(t) можно выразить зависимостью

где ${{\mu *}}\left( t \right)$ – выходная переменная.

Таким образом, в импульсном режиме управления исполнительный механизм 3 (см. рис. 6) совместно с ШИМ является линейным интегрирующим звеном относительно усредненной выходной переменной ${{\mu *}}\left( t \right)$ и входной переменной γ(t) с передаточной функцией

(35)
${{W}_{{\text{м}}}}(s) = {{{{\mu *}}\left( s \right)} \mathord{\left/ {\vphantom {{{{\mu *}}\left( s \right)} {{{\gamma }}\left( s \right)}}} \right. \kern-0em} {{{\gamma }}\left( s \right)}} = {{100} \mathord{\left/ {\vphantom {{100} {\left( {s{{T}_{{\text{м}}}}} \right)}}} \right. \kern-0em} {\left( {s{{T}_{{\text{м}}}}} \right)}}.$

С учетом интегрирующих свойств ИМ задача управляющего устройства для схемы, изображенной на рис. 6, состоит в формировании импульсного сигнала z(t), коэффициент заполнения которого γ(t) должен соответствовать производной от расчетного изменения выхода регулятора μ(t). В частности, для получения ПИД-алгоритма относительно входной переменной y(t) блок РИМ должен вычислять сигнал yр.б(t), пропорциональный производной от ПИД-закона.

Для рассматриваемого варианта на выходе регулирующего блока 1 вычисляется приращение регулирующего воздействия Δyр, которое в операторной форме определяется выражением

$\begin{gathered} \Delta {{y}_{{\text{р}}}} = \Delta {{y}_{{\text{п}}}} + \Delta {{y}_{{\text{и}}}} + \Delta {{y}_{{\text{д}}}} = {{K}_{{\text{п}}}}\Delta \left( {u - y} \right) + \\ + \,\,\frac{{{{K}_{{\text{п}}}}\Delta \left( {u - y} \right)}}{{{{T}_{{\text{и}}}}s}} - \frac{{{{K}_{{\text{п}}}}{{T}_{{\text{д}}}}s\Delta y}}{{{{{\left( {{{T}_{{\text{ф}}}}s + 1} \right)}}^{2}}}}. \\ \end{gathered} $

На выходе регулирующего блока 8 вычисляется переменная yр.б – производная регулирующего воздействия ${{\Delta {{y}_{{\text{р}}}}} \mathord{\left/ {\vphantom {{\Delta {{y}_{{\text{р}}}}} {\Delta t}}} \right. \kern-0em} {\Delta t}},$ пересчитанная на относительную скорость регулирования в долях располагаемой скорости исполнительного механизма Sм и в операторной форме имеющая вид

(36)
${{y}_{{{\text{р}}{\text{.б}}}}} = \frac{{{{K}_{{\text{п}}}}{{T}_{{\text{м}}}}}}{{100}}\left[ {s\left( {u - y} \right) + \frac{{\left( {u - y} \right)}}{{{{T}_{{\text{и}}}}}} + \frac{{{{T}_{{\text{д}}}}{{s}^{2}}y}}{{{{{\left( {{{T}_{{\text{ф}}}}s + 1} \right)}}^{2}}}}} \right].$

При этом переменная yр.б представляет собой расчетное значение коэффициента заполнения импульсного управляющего сигнала для ИМ:

Далее с учетом передаточной функции (35) исполнительного механизма с ШИМ-управлением этот сигнал интегрируется. В результате для варианта 3 получается ПИД-алгоритм, для которого регулирующее воздействие μ в операторной форме определяется выражением

${{\mu }} = {{K}_{{\text{п}}}}\left[ {\left( {u - y} \right) + \frac{{\left( {u - y} \right)}}{{{{T}_{{\text{и}}}}s}} - \frac{{{{T}_{{\text{д}}}}sy}}{{{{{\left( {{{T}_{{\text{ф}}}}s + 1} \right)}}^{2}}}}} \right].$

Полученная формула практически совпадает с приведенной ранее формулой (27), что свидетельствует о том, что варианты 2 и 3 реализации ПИД-регулятора с теоретической точки зрения являются эквивалентными.

Подобная структура (см. рис. 6) широко применяется в программном обеспечении промышленных регулирующих контроллеров, в частности Протар, Минитерм, Квинт, Текон, Simatic S7 (Siemens) и др.

Для варианта 3 в целях сокращения частоты включения ИМ аналогично предыдущим случаям рекомендуется использовать алгоритм автоматической подстройки параметра tи min на основе расчета П- и Д-составляющих с коррекцией параметров РЭ. В данном случае этот алгоритм имеет следующий вид:

$\left. \begin{gathered} {{t}_{{{\text{и}}\,{\text{min}}}}} = \left( {\left| {{{y}_{{\text{п}}}}} \right| + \left| {{{y}_{{\text{д}}}}} \right|} \right){{{{с}_{1}}{{T}_{{\text{м}}}}} \mathord{\left/ {\vphantom {{{{с}_{1}}{{T}_{{\text{м}}}}} {100}}} \right. \kern-0em} {100}}; \hfill \\ {{\Delta }_{{\text{в}}}} = {{{{t}_{{{\text{и}}\,{\text{min}}}}}} \mathord{\left/ {\vphantom {{{{t}_{{{\text{и}}\,{\text{min}}}}}} {{\text{d}}t}}} \right. \kern-0em} {{\text{d}}t}};\,\,\,\,\Delta = 2.0{{\Delta }_{{\text{в}}}}. \hfill \\ \end{gathered} \right\}$

При проверке качества процессов регулирования для АСР по варианту 3 было получено, что результаты для варианта 3 полностью совпали с данными варианта 2 при использовании предлагаемого алгоритма автоматической подстройки длительности импульсов (29) по всем приведенным показателям качества, включая число импульсов включения ИМ (при c1 = 0.25). Поэтому отдельные графики процессов для варианта 3 здесь не приводятся (совпадают с графиками на рис. 5, б). Следует отметить, что для варианта 3 в соответствии с выражением (35) требуется вычислять вторую производную входного сигнала.

ВЫВОДЫ

1. Автоматическая система регулирования с релейно-импульсным ПИ-регулятором и предвключенным дифференциатором удобна при модернизации системы с релейно-импульсным регулятором. По динамической точности этот вариант может немного уступать другим вариантам при отношении постоянной времени дифференцирования к постоянной времени интегрирования более 0.25.

2. Достоинством системы с аналоговым ПИД-алгоритмом и общим интегратором для И-составляющей и модели ИМ является то, что для него не требуется вычислять вторую производную регулируемой переменной, кроме того, отсутствует проблема накопления разности сигналов от ИМ и его модели. Этот вариант можно считать более предпочтительным при реализации ПИД-алгоритма в программном обеспечении программируемых контроллеров.

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

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

Список литературы

  1. Давыдов Н.И. Динамические характеристики электронных регуляторов ВТИ // Теплоэнергетика. 1954. № 5. С. 34–45.

  2. Александрова Н.Д., Давыдов Н.И. О настройках импульсных релейных регуляторов на малоинерционных объектах // Теплоэнергетика. 2002. № 5. С. 54–57.

  3. Жигунов В.В. Настройка ПИ-регулятора с предвключенным дифференциатором при ограничении на корневой показатель колебательности // Энергетик. 2014. № 9. С. 36–40.

  4. Денисенко В.В. ПИД-регуляторы: принципы построения и модификации. Ч. 1 // Современные технологии автоматизации. 2006. № 4. С. 66–74.

  5. Денисенко В.В. ПИД-регуляторы: принципы построения и модификации. Ч. 2 // Современные технологии автоматизации. 2007. № 1. С. 78–88.

  6. Олссон Г., Пиани Д. Цифровые системы автоматизации и управления. СПб.: Невский Диалект, 2001.

  7. Беляев Г.Б., Кузищин В.Ф., Смирнов Н.И. Технические средства автоматизации в теплоэнергетике: учеб. пособие. М.: Энергоиздат, 1982.

  8. Технические средства автоматизации. Программно-технические комплексы и контроллеры: учеб. пособие / И.А. Елизаров, Ю.Ф. Мартемьянов, А.Г. Схиртладзе, С.В. Фролов. М.: Машиностроение, 2004.

  9. Бочкарева Е.Ю., Кузищин В.Ф. Настройка длительности импульсов регуляторов с исполнительным механизмом постоянной скорости // Новое в российской электроэнергетике. 2009. № 9. С. 35–47.

  10. Ротач В.Я. Теория автоматического управления. М.: Издательский дом МЭИ, 2008.

  11. Кузищин В.Ф., Петров С.В. Настройка автоматических регуляторов с определением модели объекта второго порядка с запаздыванием по двум точкам комплексной частотной характеристики // Теплоэнергетика. 2012. № 10. С. 50–57.

Дополнительные материалы отсутствуют.