CC BY
37
Научный журнал КубГАУ, №112(08), 2015 года
1
УДК 681.513.66
05.00.00 Технические науки
СОЗДАНИЕ АДАПТИВНОЙ СИСТЕМЫ УПРАВЛЕНИЯ С ЭТАЛОННОЙ МОДЕЛЬЮ ДЛЯ СТАБИЛИЗАЦИИ ТЕМПЕРАТУРЫ АГЕНТА СУШКИ ДЛЯ ЗЕРНОСУШИЛКИ ТИПА ДСП16
Берестин Николай Константинович аспирант
Пугачев Василий Иванович к.т.н., доцент
Пиотровский Дмитрий Леонидович д.т.н., профессор, заведующий кафедрой автоматизации производственных процессов ФГБОУ ВПО «Кубанский государственный технологический университет» Краснодар, Россия
В статье описывается создание адаптивной системы управления с эталонной моделью для стабилизации температуры агента сушки. Показано, что работа замкнутых систем с управлением по модели вполне удовлетворительна и может быть внедрена в производство. Увеличение коэффициента усиления регулятора позволяет сократить на треть длительность переходного процесса, практически не изменив динамический заброс
Ключевые слова: АДАПТИВНАЯ СИСТЕМА, СТАБИЛИЗАЦИЯ, ДИНАМИЧЕСКИЕ ХАРАКТЕРИСТИКИ, УСТОЙЧИВОСТЬ, ЦИФРОВАЯ СИСТЕМА, ПЕРЕХОДНЫЙ ПРОЦЕСС
UDC 681.513.66 Technical sciences
CREATING AN ADAPTIVE CONTROL SYSTEM WITH A REFERENCE MODEL TO STABILIZE THE TEMPERATURE OF DRYING AGENT FOR GRAIN DRYER TYPE OF DSP16
Berestin Nikolay Konstantinovich Postgraduate student
Pugachev Vasiliy Ivanovich Cand.Tech.Sci., assistant professor.
Piotrovskiy Dmitriy Leonidovich Dr.Sci.Tech., professor
Kuban State Technological University, Krasnodar, Russia
The article describes the process of creation of an adaptive control system with a reference model to stabilize the temperature of drying agent. It is shown that the closed systems work with the operated models is quite satisfactory and therefore can be implemented in production. The increase of the gain of the controller allows reducing the duration of the transition process by one third, almost without changing the dynamic cast
Keywords: ADAPTIVE SYSTEM, STABILIZATION, DYNAMIC PROPERTIES, STABILITY, DIGITAL SYSTEM, TRANSITION PROCESS
Для управления объектами с переменными параметрами в [1-3] предлагалось использовать адаптивную систему управления с эталонной моделью. Это достигается введением в контур управления модели-эталона, которая соединена с основной системой управления.
На рисунке приведена схема включения модели-эталона параллельно объекту.
Выходная координата объекта управления x измеряется и сравнивается с выходным сигналом модели x м. Разность между ними вводится в цепь обратной связи, выходной сигнал u^ которой вычитается из управляющего сигнала u(t), действующего на объект.
http://ej.kubagro.ru/2015/08/pdf/28.pdf
Научный журнал КубГАУ, №112(08), 2015 года
2
Рисунок 1 - схема включения модели-эталона
Для рассматриваемой схемы справедливы соотношения: x (p ) = W 0 (p )[u(p) - (x(p) - Wm (p )u(p) )W oc (p )J = W 0 (p )[u(p) - W oc (p )x(p) + W oc (p )W m (p )- u(p) ];
x (p )|1 + W 0 (p )W oc (p )J= W 0 (p )|l + W
_ ( ) x(p) W0(p)-[1 + Woc(p)Wm(p)J
We (p)
u(p)
1 + W0 (p)Woc (p)
oc (p )W m (p )1 u (p )
We(p)» Wm (p).
Такая система обеспечивает ослабление влияния изменения динамических характеристик объекта на динамические характеристики замкнутой системы в целом.
Реализация управления по модели сталкивается с проблемой обеспечения устойчивости эквивалентной модели и всей замкнутой системы управления.
Как подобрать параметры обратной связи, которые обеспечат устойчивость эквивалентной модели при заданном диапазоне изменения параметров объекта.
Ранее [2] было показано, что интегро-дифференцирующее звено наилучшим образом подходит для этой цели. Но как выбрать значение Кос
http://ej.kubagro.ru/2015/08/pdf/28.pdf
Научный журнал КубГАУ, №112(08), 2015 года
3
и Тос, чтобы модель была устойчивой? Для этого используем расширенную амплитудно-фазовую характеристику.
Параметры системы управления:
3
1 1
Wc(p) = ------ Wi(p) =
Wm(p) :=
Womin(p) =
ААААААЛАЛЛААААЛ ^ '
Wo1(p) =
ЛЛЛЛЛЛЛЛЛЛ^ Г s
Tc • p ' Ti • p + 1
48 • p + 4.8
2
200 • p + 30 • p + 1
2
2 • p + 11.2 • p + 1 Wo1(p) = Wc(p) • Wi(p) • Womin(p)
4.800 • p + .4800
4 3 2
40. • p4 + 226. • p3 + 31.20 • p2 + 1. • p
Womax(p) =
Wo2(p) =
12 • p + 1.2
2 • p + 11.2 • p + 1 Wo2(p) = Wc(p) • Wi(p) • Womax(p)
1.200 • p + .1200
4 3 2
40. • p4 + 226. • p3 + 31.20 • p2 + 1. • p
3
Wm(p) = ------------
aaaaaaaaa4^ '
200 • p2 + 30 • p + 1,
Hm(t) = 3. • e(- 1000)4 - 6. • e(- 5000e-1)4 + 3.
ЛЛЛЛЛЛЛА4 /
Woc(p) =
Toc • p + Koc
0.01 • p + 1
Для получения некоторого запаса устойчивости для области параметров примем:
ln( 1 -y)
(__ m —
i = V-1, У = 01,
m=
ЛЛЛЛЛ
2 • p m = 0.017
Wob(p) =
1
Wo1(p)
Wob(i>w)
A1(i, w) + A2(i, w) + A3(i, w) A4(i, w)
где:
http://ej.kubagro.ru/2015/08/pdf/28.pdf
Научный журнал КубГАУ, №112(08), 2015 года
4
A1(i, w)
A2(i ,w) = A3(i, w)
3 2 4 4
= L(—261.) • i3 + 6.57 • i2 - .735e-1 • i + .308e-3 + .389e4 • i4J • w4
(.220e5 i3 — .111e4 i2 + 18.5- i —. ш4 • w3 + (.303e4 i2 — 102.- i + .854 • w2
[(—1.63) + 97.2 • i] • w A4(i,w) := (467. • i — 7.84) • w + 46.7
Koc(w) = —Re(Wob(i,w)) — m • Im(Wob(i,w))
—Im(Wob(i, w))
Toc(w) = ---------
w
w := 0,0 + 0.01.. 5.
50 40 30
Koc(w)
----- 20
10
0
0 2 4 6 8 10
Toc(w)
Рисунок 2 - область параметров для минимальной нагрузки
50 40 30
Koc(w)
----- 20
10
0
Рисунок 3 - Область параметров для максимальной нагрузки
Wo2(p)
0 2 4 6 8 10
Toc(w)
http://ej.kubagro.ru/2015/08/pdf/28.pdf
Научный журнал КубГАУ, №112(08), 2015 года
5
Использование расширенной амплитудно-фазовой характеристики облегчает и упрощает выбор параметров обратной связи систем с эквивалентной моделью. Примем:
/T°c 7,
Koc = 10
ЛЛЛЛЛЛЛЛЛ ?
Woc(p)
Toc • p + Koc
0.01 • p + 1
Woc(p) ®
7 • p + 10
.1e-1 • p + 1
Передаточные функции эквивалентных моделей:
We1(p)
Wo1(p)(Wm(p) • Woc(p) + 1) Wo1(p) • Woc(p) + 1
We2(p)
Wo2(p)(Wm(p) • Woc(p) + 1) Wo2(p) • Woc(p) + 1
WJp)
We2p)
3 2
9.600 p3 + 961.4- p2 + 244.8 p + 292.8 . 1523e5 p3 + .2072e5 p2 + 2924, p + 80. • p6 + 8456, p5 + .4569e5 p4 + 96.,
32
2.4000 p3 + 240.36 p2 + 61.212 p + 73.200 10185: p3 + 5564.1- p2 + 738.40 p + 80. • p6 + 8456, p5 + 45685: p4 + 24.,
Hel(t) = 3.05 - . 145e — 5 ■ f/400 '1*- 745e-3 ■ t[r5A3yi+ m-t(r"h~V)
—.148e — 1 ■e("-552^1)'t-cos(.663-t)-.130e- 1 ■ ■ an(.663 ■ t)
-5.99 • e(- .500e-1)•1
и , n, 1П (- 1DD )t ... - (-5.49)1 (-.9Pde-l)
He2(t) = 3.05 - ,330e- ]0 ■ & ' - .18le - 3 ■ e1 ^ +- 2.S7 - *
- 6.00- S~ ^^Vsile- 1 ■f£--293*-1>t-cos{33Q-t)-
-.390e - 1 • e(- .293e-1)•t • sin(.330 • t).
http://ej.kubagro.ru/2015/08/pdf/28.pdf
Научный журнал КубГАУ, №112(08), 2015 года
6
Рисунок 4 - Графики переходных функций эквивалентных объектов
Проведем оптимизацию параметров управляющего устройства по расширенным АФХ.
1 -ln( 1 - у)
Wob(p) = ------ m = -----------------------------
We1(p), y= 0.99, 2 -p , m = 0.733
Wob(i, w) =
C1(i, w) + C2(i, w) + C3(i, w) + C4(i, w + C5(i, w))
C6(i, w) + C7(i, w) a = 0. a(w) = w - (m2 + 1)
C0(w) = 0.5 -(a(w) - Im(Wob(i,w)) +yja(w)2 - Im(Wob(i,w))2 + 4 - a(w) - w - a)
2 - m - w - a
C1j» = m - Im(Wob(i,w)) - Re(Wob(i,w)) +
C0(w)
w = 0.0 ,0.0 + 0.001.. 0.2
http://ej.kubagro.ru/2015/08/pdf/28.pdf
Научный журнал КубГАУ, №112(08), 2015 года
7
Рисунок 5 - Линия требуемой относительной степени затухания у = 0.99
Оптимальные параметры управляющего устройства:
Kp = 0.64
Ti =
Kp
0.0124
Ti := 51.6
Проведя расчет переходных функций при различных нагрузках объекта получим графики (рисунок 6)
1.2 1
Hipil(t) 0 8
----- 0.6
Hzpi2 (t)
♦♦♦♦♦ о.4
0.2 0
0 10 20 30 40 50 60 70
t
Рисунок 6 - Переходные функции замкнутой системы при различных нагрузках объекта.
Найдем период квантования, обеспечивающий измерение регулируемой величины в дискретные моменты времени без потери информации.
.uLU
+ + > + Ж t 4- j/r
Ж4
Ж 4-
http://ej.kubagro.ru/2015/08/pdf/28.pdf
Научный журнал КубГАУ, №112(08), 2015 года
8
Амплитудно-фазовая характеристика замкнутой системы при минимальной нагрузке:
где:
44 33 55 22
F1(i,w) = 765.• w • i + 262.• w • i + 7.63-w • i + 251.-w • i + 20.9-w• i + .460
4 .4 3 .3 2 .2 7 .7
F2(i,w) = .110e5 • w • i + .583e4 • w • i + 990. • w • i + 80. • w • i
F3(i,w) = .846e4 • w6 • i6 + .457e5 • w5 • i5 + 44.9 • w • i + .460
ЛЧХ^)
2 2 Re(Wzpi1(i, w)) + Im(Wzpi1(i, w))
Рисунок 7 - амплитудно-фазовая характеристика замкнутой системы при минимальной нагрузке:
Результаты расчета:
p • 0.5
T = ------
Wc = °.86, /ww Wc , X = 182. Рекомендуем Т= 1.
Найдем цифровые модели отдельных элементов системы управления. Охватим сервомотор единичной обратной связью. Тогда
http://ej.kubagro.ru/2015/08/pdf/28.pdf
Научный журнал КубГАУ, №112(08), 2015 года
9
1
Wc(p) = -
10 • p + 1.
Будем искать цифровую модель одновременно для трех элементов:
объекта управления, измерителя и сервомотора.
48 • p + 4.8 12 • p + 1.2
Womin(p) = ----------------- Womax(p) =
Wi(p) =
2
2 • p2 + 11.2 • p + 1
1
2
2 • p2 + 11.2 • p + 1
20 • p + 1
Wo1(p) =
/WWWWWr '
Wo1(p) = Wc(p) • Wi(p) • Womin(p)
1.
3 2
8.33333 • p3 + 47.0833 • p2 + 6.50000 • p + .208333.
Wo2(p) = Wc(p) • Wi(p) • Womax(p)
1.
W2(p) = --------------------------------------------
33.3333 • p3 + 188.333 • p2 + 26. • p + .833333.
Примем период квантования Т = 1, полученный в непрерывной системе управления по модели. Передаточная функция модели:
3
Wm(p) = ---------------------
/WWWW\vr' ry
200 • p2 + 30 • p + 1, её переходная функция:
Hm(t) = 3. • e(- 1000» - 6. • e(- 5000e-1>t + 3.
лллллллА Z
Решетчатая функция:
тт , ч 0 (- .5000e-1>n , (- .25000e-1fn _
Hm(n) = 3. • e - 6. • e + 3.
-1'
Wm(z) = Hm(z) • \1 - z
ЛЛЛЛЛЛЛЛЛ^ ’ \ / '
) • (1 - z 1)
Найдя Hm(z) как прямое Z -преобразование и умножив на (1 - z 1) получаем:
Wm(z) =
ЛЛАЛАЛАЛА4 Z
.182880e-2 • z + .178365e-2
2
z2 - 1.926539 • z + .9277436
http://ej.kubagro.ru/2015/08/pdf/28.pdf
Научный журнал КубГАУ, №112(08), 2015 года
10
Правильность полученного результата подтверждает предел:
lim Wm(z) float,3 ® 3.00
z ® 1
Поскольку в цифровых системах легко вычислить производную, то для звена обратной связи:
Toc = 7
Koc = 10
Woc(p) =
Toc • p + Koc
0.01 • p + 1
Toc •
uJ
1(n) - u1(n - f'
T
+ 10 • u1(n) = 17. • u1(n) - 7. • u1(n - 1.)
У
Woc(z)
u1(z)
e2(z)
Woc(z) = 17 - 7 • z
-1
lim Woc(z) float,3 ® 10.
z ® 1
Аналогично получены дискретные передаточные функции для объекта при различных нагрузках:
Wo1(z)
/wwwwvy '
2
.734400e-2 • z2 + .123040e-1 • z + .582000e-3
3 2
z3 - 1.86851 • z2 + .876243 • z - .351800e-2
Wo2(z) =
/wwwwvy '
2
.183600e-2 • z2 + .307600e-2 • z + .145000e-3
3 2
z3 - 1.86850 • z2 + .876231 • z - .351900e-2
Дискретные передаточные функции эквивалентных моделей:
We1(z)
Wo1(z)(Wm(z) • Woc(z) + 1) Wo1(z) • Woc(z) + 1
„ s (-.183e-|) ■ z +. 173e-3 z + .93 le-4+ .622e-l- z - .729 e-1- z4+ ,290e-l- z Wel(z) =-----------------------------------------------------------
9.57z3-3.94z2+9.71z5- 12.4- z+ z -4.55- z + .570- z+ ,484e-l
http://ej.kubagro.ru/2015/08/pdf/28.pdf
Научный журнал КубГАУ, №112(08), 2015 года
11
lim We1(z) float,3 ® 3.04
z ® 1
We2(z)
Wo2(z)(Wm(z) • Woc(z) + 1) Wo2(z) • Woc(z) + 1
(-.lS2e-b-z4 + .155e-l-s3-.457e-2-z2 + .724e-2-z5 + .434e-4z + .233e-4 We3^z) =-----------------------------------------------------------
8.62-z - 111-z4 + 4.93-z3- 1.44-z2 + z7-4.55-z+. 142 z + . 12le-1
lim We2(z) float,3 ® 2.86
z ® 1
Рисунок 8 - графики переходных функций эквивалентных моделей.
Проведем расчет переходных функций замкнутой системы, работающей по модели при различных нагрузках.
Для оптимальных параметров регулятора, полученных для непрерывной системы управления по модели:
Kp = 0.64 Ti = 51.6
переходные фу4нкции замкнутой цифровой системы управления по модели имеют вид:
http://ej.kubagro.ru/2015/08/pdf/28.pdf
Научный журнал КубГАУ, №112(08), 2015 года
12
Рисунок 9 - переходные функции замкнутой цифровой системы управления по модели при минимальной и максимальной нагрузках
Попробуем сократить длительность переходного процесса увеличив коэффициент усиления регулятора. Получим новые графики переходных функций цифровой системы с эталонной моделью при минимальной Hz1(n) и максимальной Hz2(n) нагрузках (рисунок 10)
Рисунок 10 - Графики переходных функций цифровой системы с эталонной моделью при минимальной Hz1(n) и максимальной Hz2(n) нагрузках
Работа замкнутых систем с управлением по модели вполне удовлетворительна и может быть внедрена в производство. Увеличив коэффици-
http://ej.kubagro.ru/2015/08/pdf/28.pdf
Научный журнал КубГАУ, №112(08), 2015 года
13
ент усиления регулятора, удалось сократили на треть длительность переходного процесса практически не изменив динамический заброс.
Выводы:
1. При создании системы управления по модели необходимо в первую очередь обеспечить устойчивость эквивалентной модели во всем диапазоне изменения параметров объекта.
2. Весьма эффективным является введение производной в звено обратной связи эквивалентного объекта.
3. Легкость реализации производной в цифровых системах позволяет заменить интегро - дифференцирующее звено непрерывной системы на пропорционально-дифференцирующее звено в цифровой.
4. Оптимальные параметры непрерывной системы могут быть использованы при реализации цифровой системы.
5. Применение системы управления по модели обеспечивает практически одинаковые качества переходного процесса при различных параметрах объекта.
6. Применение системы управления по модели позволяет исключить из адаптивной системы необходимость устанавливать дополнительное оборудование для измерения нагрузки печи, проводить пересчет параметров настройки управляющего устройства, а также дополнительные предварительные исследования.
Литература
1. Пугачев В.И. Исследование возможностей адаптивной системы управления с эталонной моделью // Изв. Вузов. Пищевая технология. -2007. -№ 2. -С 80-83.
2. Пугачев В.И. Устойчивость эквивалентной модели системы управления с эталонной моделью и объектом третьего порядка // Изв. Вузов. Пищевая технология. -2007. -№ 4. -С 88-91.
3. Пугачев В. И. Влияние чистого запаздывания на устойчивость эквивалентной модели адаптивной системы управления с эталонной моделью // Изв. Вузов. Пищевая технология. -2007. -№ 5-6. -С 76-78.
http://ej.kubagro.ru/2015/08/pdf/28.pdf
Научный журнал КубГАУ, №112(08), 2015 года
14
4. Кудрявцев Е.М. Mathcad 11. Полное руководство по русской версии. - М.: МДК Пресс, 2005. -592 с.
References
1. Pugachev V.I. Issledovanie vozmozhnostej adaptivnoj sistemy upravlenija s jetalon-noj model'ju // Izv. Vuzov. Pishhevaja tehnologija. -2007. -№ 2. -S 80-83.
2. Pugachev V.I. Ustojchivost' jekvivalentnoj modeli sistemy upravlenija s jeta-lonnoj model'ju i ob#ektom tret'ego porjadka // Izv. Vuzov. Pishhevaja tehnologija. -2007. -№ 4. -S 88-91.
3. Pugachev V.I. Vlijanie chistogo zapazdyvanija na ustojchivost' jekvivalentnoj modeli adaptivnoj sistemy upravlenija s jetalonnoj model'ju // Izv. Vuzov. Pishhevaja tehnologija. -2007. -№ 5-6. -S 76-78.
4. Kudrjavcev E.M. Mathcad 11. Polnoe rukovodstvo po russkoj versii. - M.: MDK Press, 2005. -592 s.
http://ej.kubagro.ru/2015/08/pdf/28.pdf
