CC BY
68
УДК 681.52.01
К ОБОСНОВАНИЮ МЕТОДА СТРУКТУРНО-ПАРАМЕТРИЧЕСКОГО СИНТЕЗА АВТОМАТИЧЕСКИХ РЕГУЛЯТОРОВ
С.В. Панферов, А.И. Телегин, В.И. Панферов
ON JUSTIFICATION OF THE METHOD OF STRUCTURE PARAMETRIC SYNTHESIS OF AUTOMATIC REGULATORS
S.V. Panferov, A.I. Telegin, V.l. Panferov
Обосновывается вид эталонной передаточной функции замкнутой системы автоматического регулирования. Решается задача структурно-параметрического синтеза автоматических регуляторов для типовых промышленных объектов управления. Анализируется качество выбора структур и настроек регуляторов.
Ключевые слова: синтез регулятора, структурно-параметрический синтез, выбор настроек, выбор структуры.
A type of reference closed-loop transfer function of automatic regulation is justified. The problem of structure parametric synthesis of automatic regulators for standard industrial control objects is solved. The quality of structure selection and regulator adjustments is analyzed.
Keywords: regulator synthesis, structure parametric synthesis, superstructure selection, structure selection.
В настоящее время в промышленности для автоматического регулирования различных переменных технологических процессов широко используются ПИД-регуляторы и их частные варианты. Накоплен огромнейший опыт применения таких регуляторов, в частности, разработаны и апробированы различные способы их настройки. Вместе с тем, как неоднократно отмечает В.Я. Ротач в своем сравнительно недавно изданном учебнике [1], П-, ПИ- и ПИД-«...алгоритмы были получены чисто эвристическим путем» [1, с. 82] и что «... достаточно убедительное формальное доказательство целесообразности их применения ... до сих пор получить не удалось» [1, с. 24].
В настоящей работе предлагается метод структурно-параметрического синтеза регуляторов по желаемой (эталонной) передаточной функции замкнутой системы. Работа является развитием результатов, полученных в работах [2, 3] и, как нам представляется, может служить некоторым формальным обоснованием целесообразности применения ПИД-регуляторов.
Панферов Сергей Владимирович - аспирант кафедры систем управления и математического моделирования Миасского филиала ЮУрГУ; [email protected].
Телегин Александр Иванович - д.ф.-м.н., профессор, заведующий кафедрой систем управления и математического моделирования Миасского филиала ЮУрГУ; [email protected].
Панферов Владимир Иванович - д.т.н., профессор, заведующий кафедрой теплогазоснабжения и вентиля-ции ЮУрГУ; [email protected].
1. Решение задачи синтеза регуляторов
Задача решалась методом выбора передаточной функции замкнутой системы в виде
fV3C(p)=—-—ехр(-тоб/?), где 0 - достаточно ма-вр+1
лая постоянная времени, р - комплексная переменная, а параметр приравнивался ко времени запаздывания объекта управления, который при этом представлялся одной из следующих типовых передаточных функций:
1 -ехр(-тоб/г);
То5Р
об
TrfP+l
к,
ехр(-тоб/>);
Об
-exp(-To6jp),
С12Р +ЩР+1
где ко6, Гоб, - соответственно коэффициент передачи, постоянная времени и время запаздывания объекта, ах, а2 - коэффициенты дифференциального уравнения объекта. В результате приме-
Panferov Sergey Vladimirovich - post-graduate student of control system and mathematical modeling department of Miass branch of SUSU; [email protected].
Telegin Alexander Ivanovich - PhD, professor, head of control system and mathematical modeling department of Miass branch of SUSU; [email protected].
Panferov Vladimir Ivanovich - PhD, professor, head of heat and gas supply and ventilation department of SUSU; [email protected].
нения данного подхода получили, что для объектов первого типа передаточная функция квазиоп-тимального регулятора будет равна
т
Ж,(р) =—2^—, т.е. близким к идеальному являет-тоб
ся П-регулятор с коэффициентом передачи
к„ = ■
об
Для объекта второго типа квазиопти-
тоб
мальным будет ПИ-регулятор, его передаточная функция будет иметь вид
Щр) =
об
^(Тоб+0)
1+-
1
ТобР
Передаточная
функция регулятора для объекта третьего типа
Оі
будет такой И? (р)=-
•об(хоб +е)
1 а2 1+--+—р
«і Р
•Ч
т.е. целесообразно применение ПИД-регулятора.
2. Обоснование вида эталонной передаточной функции
При выборе эталонной передаточной функции замкнутой системы руководствовались следующими соображениями. Очевидно, что идеальной передаточной функцией замкнутой системы по задающему воздействию является передаточная функция вида Игзс(р) = 1 [4]. В этом случае САР абсолютно точно отрабатывает задание, а также полностью исключает влияние возмущений на процесс управления [4]. Однако, как это достаточно хорошо известно научной общественности, добиться такой передаточной функции совершенно невозможно. Поэтому есть смысл попытаться за счет выбора регулятора получить такую передаточную функцию замкнутой системы, которая в определенной мере будет близка к идеальной. Нетрудно видеть, что при малом значении параметра 0 следующие передаточные функции близки к 1:
в * ^, е~вр, причем при 0 0 предел этих пере-
даточных функций будет точно равен 1. Кроме 1
того
вр+1
*е вр, т.е. данные передаточные функ-
ции приближенно равноценны. Отметим также, что по данным работы [5] система с передаточной
функцией —-— является оптимальной по робаст-вр+1
ности и точности. В этой работе, в частности, приводятся следующие интересные для рассматриваемой проблемы данные: «...чем ближе к отрицательной вещественной полуоси располагаются корни характеристического уравнения системы», «...тем большую робастность имеет система» и что «...если все полюсы системы находятся на отрицательной вещественной полуоси, то она» «...обладает» «...наиболее высоким потенциалом по робастности», «...из множества чисто инерци-
онных систем наиболее структурно-робастным является простейшее инерционное звено первого порядка».
Кроме того, дополнительным обоснованием для выбора данной передаточной функции в качестве эталона являются следующие соображения.
Известно [6], что достаточно рациональным (вполне предпочтительным) является следующий критерий качества переходных процессов в САР:
1= |[є2 +а2(</є/<*)2]<*,
(1)
где е(0=х (0~х(0 - ошибка регулирования (рассогласование), дг3(0 и *(/) - соответственно заданное и действительное значение регулируемой величины, * - время, а2 - некоторый весовой коэффициент.
Известно также [6], что оптимальным по минимуму этого критерия переходным процессом является экспоненциальный процесс, т.е. процесс вида
е(0=е(0)ехр(-г/а), (2)
где е(0) - значение ошибки регулирования при *=0.
Если при этом считать, что такой переходный процесс должен иметь место при отработке САР
единичного ступенчатого задания х3(0= 1(0, то в
этом случае е(0)=1 и выходной сигнал САР будет
иметь вид
х(/)=1-ехр(-^/а). (3)
В связи с этим Лапласово изображение выходной величины запишется так 1
Цх(0}=Х(р)=
р(ар+1)
(4)
(5)
далее, учитывая, что
цх\т=х3(Р)=-р
найдем передаточную функцию образцовой (эталонной) САР
Х(р)_ 1
Щр)=
х3(р) ар+1
(6)
что соответствует вышеотмеченному.
Для выбора величины а имеются следующие рекомендации [6]: так как длительность переходного процесса в САР с передаточной функцией (6) составляет примерно (3-5-4)а, поэтому если задано время регулирования /Р, то, допуская известный запас, а следует вычислять по соотношению
а=—— ■ (7)
(5+6)
Понятно, что параметр а в этом случае играет роль постоянной времени 0 эталонной передаточной функции САР.
Если известна передаточная функция замкнутой САР по заданию РУЗС (р), то передаточная функция регулятора ^(р) находится по формуле:
*Х(Р)
ад=-
(8)
КвШ-ЯкСрЯ
где ^ (р) - передаточная функция объекта.
Если динамические свойства объектов управления описывать вышеприведенными передаточными функциями и при этом Щх;(р) выбрать в
1
виде
0/7 + 1
, то в соответствии с формулой (8) бу-
дут получаться физически нереализуемые структуры регуляторов из-за наличия в числителе РУр(р) сомножителя вида ехр(тобр). Поэтому,
учитывая все вышеизложенное, при решении задачи полагали, что Щс(р)=—-—ехр(-тоб/>)- Таким
0/?+1
образом, было и установлено, что для каждого конкретного объекта управления, принадлежащего множеству типовых динамических объектов целесообразно применение конкретного регулятора из ПИД-семейства.
3. Обсуждение подхода к решению задачи
Следует отметить, что данный подход к решению задачи синтеза (посредством предварительного выбора передаточной функции замкнутой системы) вообще-то известен уже давно. По-видимому, впервые в отечественной литературе идея этого метода достаточно ясно изложена в работе [7, с. 7-8]. Многие методы синтеза САР, в частности, и методы [4] по существу основаны на предварительном выборе желаемой передаточной функции замкнутой системы. В известной мере этот подход используется даже при выборе желаемой логарифмической амплитудно-частотной характеристики разомкнутой системы [8-10]. Обсуждаются и используются такие подходы и в настоящее время [5, 11, 12]. Однако непосредственно этот метод в литературе практически не рассматривается, во всяком случае, авторам неизвестны работы конкретно с вышеизложенной постановкой и решением задачи синтеза. Вместе с тем, метод достаточно прост и отличается предельной ясностью. Очень важно, что метод приводит к явному формульному решению задачи, дает однозначный ответ на вопрос: какой регулятор следует применить на данном объекте и каковы должны быть его параметры настройки. Последнее замечание достаточно значимо, по крайней мере, в учебных целях, так как в литературе по существу нет формализованных процедур, позволяющих однозначно решать задачу синтеза регуляторов. Так, например, в [7, с. 85], в частности, отмечается, что каждый из типовых переходных процессов «...может быть обеспечен любым регулятором
непрерывного действия». При этом в работах [7, 13] для любого типового переходного процесса предлагаются настройки любого из типовых регуляторов. Если учесть, что требование простоты закона регулирования в настоящее время уже не является существенным, то следует признать неоднозначность рекомендаций работ [7, 13] по выбору структуры регуляторов в современных условиях. Определенная неоднозначность прослеживается также и в известных процедурах параметрической настройки уже выбранного каким-то образом закона регулирования. Здесь достаточно упомянуть о показателе колебательности М, ведь в работе [9] требуется, чтобы он не превышал значения, равного двум, а в работе [14, с. 59-60] вообще утверждается, что «...обоснованные рекомендации по выбору М устанавливают на основании опыта эксплуатации каждого класса систем», т.е. уже после проектирования и внедрения САР. Таким образом, рассмотренный метод, на наш взгляд, незаслуженно обделен вниманием.
Примечательно то, что в ряде случаев полученные настройки достаточно близки к рекомендациям работ [7, 13]. Для объекта второго типа полученные настройки ПИ-регулятора при 0 = 0 полностью совпадают с настройками работы [13] для переходного процесса с минимальной квадратичной площадью отклонения от идеального результата. Для объекта первого типа настройка П-регулятора для апериодического переходного процесса и переходного процесса с 20 % перерегулированием в работе [13] отличается от полученной только на постоянный множитель. Для переходного процесса с минимальной квадратичной площадью отклонения настройка П-регулятора в работе [13] отсутствует, причем никаких разъяснений на этот счет не приводится. Для объекта третьего типа в работах [7, 13] вообще нет никаких рекомендаций. Подчеркнем, однако, что в работах [7, 13] не указывается, как была найдена структура приводимых формул для настройки.
4. Свойства сконструированных систем
Нетрудно видеть, что во всех рассмотренных случаях передаточная функция разомкнутой системы получается одинаковой и равной
ехрО-Тоб/?)
(б + 'Гоб)^
функция замкнутой САР по задающему воздействию будет равна 1¥зс(р) =—------С*р( -----------
(0+*об)Р+езФ(-тобР) Отсюда следует, что дифференциальное уравнение замкнутой САР будет иметь вид:
^рс(Р) =-
поэтому передаточная
(6 + тоб)“^Г^ + л:(/“тоб)=:х3(?-тоб) • ш
(9)
Решая уравнение (9) для единичного ступенчатого воздействия, найдем переходную функцию замкнутой системы. На рисунке приведены переход-
ные функции САР для = 5 с и различных значений параметра 0 : для кривой 1 0 = 0 с; для кривой 2 0 = 3 с; для кривой 3 0 = 5 с; для кривой 4 0 = 8 с.
Как видно из рисунке во всех случаях САР устойчивая, причем с увеличением параметра 0 колебательность переходной функции уменьшается, переходный процесс все в большей степени приобретает апериодический характер (демпфиро-ванность САР возрастает), причем при f —> оо, х —»1, т.е. САР астатическая.
Переходные функции САР
Проанализировали влияние параметра на показатели качества переходных процессов в рассматриваемых САР, определялись перерегулирование ст, время регулирования fP, а также значения следующих критериев:
'к
/,= J|e(f)|<*; о 'к
/2= j в2 (*)<#,
(10)
(П)
где fK - конечное время оценки переходного процесса. Здесь время регулирования /Р определялось как время, по истечение которого отклонение регулируемой величины от задания не будет превышать 5 %. Результаты приведены в табл. 1.
Как видно из табл. 1, перерегулирование ст с увеличением параметра растет, начиная с величины 45 %, однако рост достаточно медленный, в наших расчетах эта величина не превышала 50 %, fP также растет, причем практически прямо пропорционально увеличению , растут и величины /) и /2, которые оценивают близость переходных процессов к идеальному результату - единичной функции 1(f). Таким образом, можно сделать вывод, что чем больше запаздывание объекта управления Tgg, тем хуже переходный процесс в системе, сконструированной по рассматриваемому методу. Вместе с тем удалось установить достаточно интересное влияние параметра 0, небольшое положительное значение этого параметра приближает САР к идеалу, уменьшаются все перечисленные показатели (табл. 2). Отметим также, что при некоторых сочетаниях параметров и 0 показатели качества переходного процесса, особенно перерегулирование ст резко меняются, так, например, если т = 1 с, то при 0 = 0с ст = 45%, а при 0 = 1 с уже ст = 2,5%, при 0 = 2с, ст = 0%. При этом очевидно, что во многих случаях за счет выбора параметра 0 всегда можно добиться приемлемого качества процесса регулирования.
Таблица 1
Показатели качества переходных процессов в САР при 0 = 0 с и возмущении по заданию
тоб>с 1,0 2,0 3,0 5,0 7,0 10,0 15,0 20,0
СТ, % 45,0 47,5 48,3 48,9 49,2 49,4 49,6 49,7
fp, С 8,2 16,8 25,4 42,5 70,4 101,6 153,2 204,8
11, ед. регулируемой величиныхс 2,57 5,45 8,35 14,12 19,75 28,63 43,12 57,62
/2, ед. регулируемой величины2 ХС 15,62 32,58 49,60 83,66 117,72 168,84 254,04 339,25
Таблица 2
Показатели качества переходных процессов в САР при 0 = 1 с и возмущении по заданию
тоб> с 1,0 3,0 5,0 10,0
СТ, % 2,5 24,1 32,6 40,4
fp,C 3,4 19,4 33,5 84,6
1г, ед. регулируемой величиныхс 2,10 6,51 11,67 25,46
/2, ед. регулируемой величины2 ХС 15,39 45,19 76,69 158,95
Если проанализировать передаточные функции САР по каналу «возмущение со стороны регулирующего органа - рассогласование», то несложно установить, что квазиоптимальные системы с объектами второго и третьего типов также будут астатическими и в этом случае, однако система с объектом первого типа будет только статической, причем установившаяся ошибка регулирования будет равна 1/кр = (хоб+в)/Тоб, т.е. она тем
меньше, чем больше коэффициент передачи 11-регулятора.
Используя критерий устойчивости Найквиста для систем с запаздыванием [8], нашли, что критическая частота будет равна = *
, по-
б+^об
этому условие устойчивости замкнутой САР в общем случае запишется так: тоб/(в+хо6)<п/2 . В нашем же случае всегда тоб/(0+тоб)<1, поэтому рассматриваемая САР устойчива.
Отметим также, что в данной работе структуру и настройки автоматических регуляторов выбирали, полагая, что передаточная функция замкнутой системы должна быть следующей
Мзс(р) = ехР( товР) Однако на самом деле полу-6р+1
чилось, что fV3C(p) = -
exp(~W)
. Не-
(9 + *обХР + в>ф(-‘гоб/0 трудно видеть, что поставленная цель все таки достигнута, но с той степенью точности, с какой выполняется соотношение ехр(-тобр)»1-х^р.
5. Некоторые замечания
В.Я. Ротачем в [1, с. 188] предлагается так называемое «...формальное обоснование применимости ПИД-регулятора»: оптимальная передаточная функция регулятора выбирается из условия минимума дисперсии выходной величины САР при отработке эквивалентного возмущения, «...приведенного непосредственно к выходу объекта» [1, с. 180]. Здесь получаемый ответ зависит от вида основного возмущающего воздействия, для возмущения, являющегося, например, белым шумом, оптимальная передаточная функция регулятора будет одной, а для возмущения другого типа совсем другой. Заметим, что, применяя этот метод, иногда можно получить результаты, достаточно близкие к тем, что следуют из вышеизложенной методики. Так, например, для объекта управления с передаточной функцией третьего типа ВЛ. Рогач, полагая, что САР обладает высоким показателем технологической работоспособности, т.е. при низкочастотном характере основных возмущений получает совершенно такой же результат, что получен нами для И'зс(р)=ехр(-хо6р). Вместе с тем, как нам представляется, вышеизложенный способ синтеза, во-первых, заметно проще по содержанию, чем способ [1], во-вторых, не требует опреде-
ления характера возмущающих воздействий и, в-третьих, предпочтительнее потому, что задача у САР вообще-то двоякая: как можно точнее отрабатывать задание, влияние же возмущений на регулируемую величину при этом должно быть как можно менее заметным (САР должна отрабатывать два вида воздействий - это задающие и возмущающие воздействия). При этом хорошо известно, что «...реакция регулируемой величины на скачок задающей позволяет судить о том, насколько хорошо оптимизирован контур и по отношению к возмущениям» [4, с. 6], т.е. безусловно предпочтительнее использовать способ синтеза, ориентированный на оптимальную отработку задающего воздействия. Кроме того, аргументы, приведенные ВЛ Ротачем [1, с. 190, с. 191]), являются, по существу, одним из вариантов обоснования структуры регулятора Ресвика.
В работах А.З. Сметаны [15-18] утверждается, что для объекта с передаточной функцией в виде звена чистого запаздывания fFo6(p)=exp(-xo6p) найдена (на оптимизаторе!?) передаточная функция оптимального регулятора, эта передаточная функ-
. 0,75(0,5 тбр+1) ция имеет вид WJp)=-------------—------, очевидно,
что это ПИ-регулятор. По-видимому, никакой формализованной процедуры выбора структуры регулятора у автора нет, речь может идти лишь об определении оптимальных параметров настройки регулятора, при этом критерием оптимальности здесь
00
является функционал / = jj е| dt. Вместе с тем заме-
о
там, что в работе [19, с. 242] существенно аргументировано (достаточно обоснованно) показывается, что для звена чистого запаздывания идеальным регулятором является И-регулятор, такой же результат получается и при применении рассматриваемого способа синтеза САР: действительно, подставляя ^об (Р)= ехР(~тобP) в (8)и учитывая используемые
допущения, получим, что Wn (р)=--------^---.
р (0+ТобХр
Далее, сообщается, что если объект имеет передаточную функцию вида
Wo6(p)=exp(-zo6p)W1(p), то передаточная функция оптимального регулятора будет такой:
0,75(0,5 тоб/7+1)
. При этом автор предла-
*ОбРЩР)
гает описывать широкий класс объектов управления следующей, якобы «нестандартной», передаточной функцией: крб
^об=-
-ехр(-тоб/>),
Т0бр+ехр(-^р)
где А, - некоторый параметр. Выбор такой передаточной функции автор в работе [15] объясняет тем, что «...большинство объектов регулирования в теплоэнергетике описываются передаточными
функциями, характеристические уравнения которых имеют комплексные корни», «...известные модели типа ехр(-хо6р)/(То6р +1) или более
сложные ехрС-Тодр)/^^/? +1)(7’1р + 1)"] имеют характеристические уравнения только с действительными корнями». Нетрудно видеть, что данный предлог является достаточно надуманным, никакой необходимости во включении в этом случае в знаменатель передаточной функции трансцендентного члена е~Хр нет, достаточно записать знаменатель передаточной функции объекта в виде многочлена общего вида относительно р, как известно, всякий такой многочлен с действительными коэффициентами и со степенью, начиная со второй и выше, может иметь как действительные, так и комплексные корни. Затем же А.З. Сметана указывает, что данная передаточная функция «...может быть преобразована в передаточную функцию третьего или второго порядка с запаздыванием», т.е. происходит возврат к так называемому стандартному представлению динамических свойств объектов управления. Поэтому естественно возникает вопрос: так в чем же заключается оригинальность и новизна работы?
Если использовать метод А.З. Сметаны, то передаточная функция разомкнутой САР будет иметь вид:
цг _ °>75(0,5ТрбР+ЦехрН:,^)
150 *об Р
поэтому передаточная функция замкнутой САР по задающему воздействию будет равна:
цг ( 0>75 (0,5 Х<*Р +1) ехР(-тобР)
^обР+'О,75 (0,5 хобр+1) ех^-т^р) ’ Отсюда следует, что дифференциальное уравнение замкнутой САР будет следующим:
(,Л/0,75)^+0>5 +
ш ш
х(/-тоб)=0,5тоб щ-тоб), (12)
ш
т.е. получается довольно близкая к (9) структура. Однако следует заметить, что все результаты получены автором работ [15-18] достаточно сложным, запутанным и не совсем «прозрачным» путем. В работе [16] приводится условие устойчивости систем с такими передаточными функциями, это условие применительно к нашему случаю полностью совпадает с вышеприведенным неравенством. В этой же работе [16] указывается, что если при этом 0<тоб/(6+тоб)<0,37, то переходный процесс в САР монотонный, если
О,37<тоб/(0 + тоб)<я/2 - колебательный затухающий, при / (0+Тоб) > я/2 переходный процесс колебательный расходящийся. Численное моделирование показало также, что область моно-
тонных переходных процессов, выделяемая автором работы [16] неравенством
О<тоб/(0+тоб)<О,37 подтверждается, однако, только приближенно.
Как показал анализ литературных источников, некоторые из вышеприведенных результатов при определенных допущениях можно получить и методами работ [1, 19]. Однако, в частности, подход работы [19] излишне усложнен: задача решается методом эвристического конструирования алгоритмической структуры идеальной САР, при этом в процессе решения (задачи структурнопараметрического синтеза) используются структуры практически нереализуемых регулятора Ресви-ка и предиктора Смита. Найденная таким образом структура идеального регулятора, как отмечает и сам автор [19, с. 240, с. 241], имеет ряд недостатков. Так в работе [19, с. 241], в частности, утверждается, что «...существенным недостатком системы с регулятором Ресвика является ее критичность или сильная чувствительность к малым вариациям запаздывания объекта: система устойчива только при точном равенстве запаздывания объекта и запаздывания» модели. Однако при этом никаких конкретных данных не приводится, по-видимому, этот вопрос - вопрос исследования грубости структур с регуляторами типа Ресвика все-таки еще требует детального исследования. Предложенная процедура решения задачи, на наш взгляд, заметно проще, причем на всех этапах решения получаются структурно устойчивые регуляторы и системы, кроме того, конкретно робастность предлагаемых нами структур и настроек регуляторов, насколько это нам известно, еще не исследована.
В работе [20] задача структурнопараметрического синтеза автоматических регуляторов решается методом В.Я. Ротача [1]. Если найденная таким образом передаточная функция регулятора не выражается через передаточные функции П-, И-, ПИ-, ПИД-регуляторов, то методом разложения в цепную дробь, а в сложных случаях, дополнительно, проводя предварительную декомпозицию задачи, эту передаточную функцию представляют в виде соединения «...двух-трех типовых законов» [20]. Данная процедура применяется авторами также и к так называемому восста-новительно-прогнозирующему регулятору (ВП-регулятору). Таким образом, устанавливается приближенная связь между квазиоптимальными и типовыми законами регулирования.
Понятно, что недостатки метода [1] целиком и полностью присутствуют и в работе [20], кроме того, метод работы [20] нельзя назвать полностью формализуемым.
Следует заметить, что представление динамических свойств реальных объектов управления приведенными передаточными функциями, как известно, является достаточно приближенным.
Неоптимальными могут быть как структура математической модели объекта, так и численные значения ее параметров, т.е. возможно неоднозначное решение этих вопросов. Так, например, в работах [21, 22] утверждается, что «...во многих случаях ... динамика объектов регулирования с достаточной точностью может быть представлена передаточными функциями без учета транспортного (чистого) запаздывания» [22, с. 56] и что запаздывание авторами называется «кажущимся» [22, с. 56]. Вместе с тем следует отметить, что подавляющее большинство исследователей все-таки отдают предпочтение моделям с запаздыванием, здесь следует отметить работы [5-7, 9, 13, 14, 19]. Звено чистого запаздывания, если оно присутствует в модели объекта управления, переносится в характеристическое уравнение замкнутой системы, что делает невозможным применение некоторых аналитических методов синтеза автоматических регуляторов. Этим, как нам представляется, в частности, и объясняется стремление некоторых авторов исключить это затруднение. Кстати говоря, структура прогнозатора Смита такова, что позволяет уничтожить именно запаздывание в характеристическом уравнении замкнутой системы.
Так как качество выбора и настройки автоматических регуляторов напрямую зависит от решения задачи построения математической модели объекта управления, то качество переходных процессов в реальных системах может отличаться от ожидаемого и весь вопрос состоит в том, насколько в каждом конкретном случае будет допустимым это отличие. Сконструированная САР, в частности, может и не удовлетворять такому требованию как робастность [25], этот вопрос, как и в целом вопрос о приемлемости для практики рассмотренного способа синтеза САР требуют конкретных исследований и достаточно тщательной проверки. Вместе с тем нельзя не упомянуть, что уже частичное одобрение и применение способа можно усмотреть в работах [7, 13], а также в работах [9,
14, 23, 24] и других. В данном же случае ограничимся исследованием того, как точность определения времени запаздывания объекта сказывается на свойствах системы. Выполненные расчеты показали, что если параметры объекта определены точно, то при 0 = 0 с перерегулирование ст = 49 %. Если же параметры объекта определены неточно, причем если время запаздывания модели тмод меньше времени запаздывания, имеющего
место на реальном объекте, то6 то перерегулирование возрастает. Так, например, если тмод =4 с, а = 5 с , то а = 72 % . САР в такой ситуации, когда из-за погрешности параметрической идентификации получилось, что тмод < то6, может и
потерять устойчивость. В самом деле, если структура математической модели объекта и все ее па-
раметры, кроме времени запаздывания, определены достаточно точно, то во всех рассмотренных случаях передаточная функция разомкнутой системы получается одинаковой и равной
\¥ж(р) _ехР( тобР) ^ поэхому сдр будет устой-
(0 + тмод)Р
чива лишь при удовлетворении следующего неравенства: Тоб /(0+тмод ) < 71 / 2 .
Влияние тмод > тоб аналогично увеличению параметра 0, т.е. демпфированность САР возрастает.
В заключение заметим, что в работах [23, с. 132, 26, с. 182] указывается, что передаточная функция замкнутой системы может быть выбрана в виде инерционного звена с чистым запаздыванием, однако при этом порядок инерционного звена не оговаривается и никаких аргументов, обоснований для такого выбора не приводится, что позволяет заключить, что у авторов по данному вопросу нет никаких особых предпочтений. По нашему мнению, данные предложения - это ничто иное, как один из вариантов задания передаточной функции замкнутой системы.
Заключение
Приводится некоторое формальное обоснование применимости семейства ПИД-регуляторов на промышленных объектах. Анализируются свойства сконструированных систем автоматического регулирования.
Литература
1.Ротач,В.Я. Теория автоматического управления: учебник для вузов / В. Я. Ротач. - М.: Изд-во МЭИ, 2004. - 400 с.
2. Панферов, В. И. Об одном подходе к решению задачи выбора и настройки автоматических регуляторов / В. И. Панферов // Известия Челябинского научного центра. - 2004. - Вып. 4(26). -
С. 139-144.
3. Панферов, В. И. Выбор и настройка автоматических регуляторов в системах теплоснабжения /
В. И. Панферов // Вестник ЮУрГУ. Серия «Строительство и архитектура». — Челябинск: Изд-во ЮУрГУ, 2005. -Вып. 3, № 13(53). - С. 81-84.
4. Фрер, Ф. Введение в электронную технику регулирования / Ф. Фрер, Ф. Орттенбургер. - М.: Энергия, 1973. - 190 с.
5. Мань, Н. В. Оптимальный синтез робастной каскадной автоматической системы управления / Н. В. Мань // Теплоэнергетика. - 2000. - № 9. -С. 22-28.
6. Воронов, А. А. Основы теории автоматического управления. Автоматическое регулирование непрерывных линейных систем / А. А. Воронов. -М.: Энергия, 1980. — 312 с.
7. Копелович, А. П. Инженерные методы расчета при выборе автоматических регуляторов /
А. П. Копелович. - М.: Металлургиздат, 1960. -190 с.
8. Теория автоматического управления: ч. 1: Теория линейных систем автоматического управления/ Под ред. А. А. Воронова. — М. : Высш. школа, 1977.-304 с.
9. Теория автоматического управления: учебник для вузов / Под ред. А. В. Нетушша. - М.: Высшая школа, 1976. - 400 с.
10. Цыпкин, Я. 3. Основы теории автоматических систем /Я. 3. Цыпкин. - М.: Наука, 1977. -560 с.
11.Гончаров, В.И. Синтез робастных регуляторов низкого порядка / В. И. Гончаров, А. В. Лиепинъш, В. А. Рудницкий // Изв. РАН. Теория и системы управления. - 2001. -№ 4. - С. 36-43.
12. Лозгачев, Г. И. Построение модальных робастных регуляторов по передаточной функции замкнутой системы / Г. И. Лозгачев, Л. А. Тютюнникова // Изв. РАН. Теория и системы управления. -2006. -№ 4.- С. 5-8.
13. Проектирование систем контроля и автоматического регулирования металлургических процессов / Г. М. Глинков, В. А. Маковский,
С. Л. Лотман, М. Р. Шапировский; под ред. Г. М. Глинкова. - М. : Металлургия, 1986. — 352 с.
14. Кулаков, Г. Т. Инженерные экспресс-мето-ды расчета промышленных систем регулирования: справ, пособие / Г. Т. Кулаков. - Минск: Вышэйшая школа, 1984. -192 с.
15. Сметана, А. 3. Методика расчета параметров настройки систем автоматического регулирования теплоэнергетических процессов / А. 3. Сметана // Теплоэнергетика. — 2002. —№10. -С. 40-45.
16. Сметана, А. 3. Методика определения передаточной функции линейного динамического объекта по его переходной характеристике / А. 3. Сметана // Изв. РАН. Энергетика. — 1998. — №2. -С. 142-155.
17. Сметана, А. 3. Методика определения параметров настройки регуляторов теплоэнергетических процессов / А. 3. Сметана // Изв. РАН. Энергетика. — 2001. -№ 2. —С. 80-87.
18. Сметана, А. 3. Автоматическая и автоматизированная настройка регуляторов теплоэнергетических процессов / А. 3. Сметана // Теплоэнергетика. - 2004. - №11.- С. 47-52.
19. Лукас, В. А. Теория автоматического управления/В. А. Лукас. -М.: Недра, 1990. - 416с.
20. Кулаков, С. М. О приближенном соответствии между квазиоптималъными и типовыми законами управления / С. М. Кулаков, В. В. Штефан, С. П. Огнев, И. А. Штефан // Изв. вузов. Черная металлургия. -1999. -№ 4.- С. 33-40.
21. Чертков, Н. К. Аналитические формулы оптимальной настройки авторегуляторов / Н. К. Чертков, С. В. Корябина // Теплоэнергетика. -1969.-№9. -С. 28-30.
22. Чертков, Н. К. Пакет компьютерных программ для настройки систем автоматического регулирования / Н. К. Чертков, В. Н. Чертков // Теплоэнергетика. -2007. -№ 9.- С. 56-60.
23. Скаржепа, В. А. Тиристорные цифровые регуляторы температуры / В. А. Скаржепа, К. В. Шелехов, А. С. Герасимов. - Киев: Техника, 1979. -144 с.
24. Петелин, Д. П. Автоматизация производственных процессов текстильной промышленности: учебник для вузов, в 5 кн.; кн. 1: Основы автоматики и технические средства автоматизации в текстильной промышленности/Д. П. Петелин, Э. М. Ромаш, А. Б. Козлов и др. - М.: Лег-промбытиздат, 1992. — 240 с.
25.Дорф, Р. Современные системы управления / Р. Дорф, Р. Бишоп - М.: Лаборатория Базовых Знаний, 2002. - 832 с.
26. Медведев, Р. Б. АСУ ТП в металлургии / Р. Б. Медведев, Ю. Б. Бондарь, В. Д. Романенко. -М.: Металлургия, 1987. - 256 с.
Поступила в редакцию 2 октября 2008 г.
