Синтез систем автоматического управления в биотехнологии методом стандартных коэффициентов
.pdfРешетняк Е.П., Рыбин М.Б., Лямина Т.В.
СИНТЕЗ СИСТЕМ АВТОМАТИЧЕСКОГО УПРАВЛЕНИЯ В БИОТЕХНОЛОГИИ МЕТОДОМ СТАНДАРТНЫХ КОЭФФИЦИЕНТОВ
Саратов 2006
1
УДК 62
ББК 32.96
Решетняк Е.П., Рыбин М.Б., Лямина Т.В.
Синтез систем автоматического управления в биотехнологии методом стандартных коэффициентов / Решетняк Е.П., Рыбин М.Б., Лямина Т.В. – Саратов: Саратовский государственный аграрный университет имени Н.И. Вавилова, 2006. – 17 с. — опубл. 15.06.2012. — Электрон. версия печ. публ. — Доступ с сайта ЭБС IPRbooks.
В задачах синтеза САУ в биотехнологии в настоящее время используются два метода: графоаналитический и интегральных оценок. В работе предложена достаточно простая и эффективная методика синтеза САУ биотехнологическими объектами, основанная на использовании метода стандартных коэффициентов. Ее высокая эффективность показана на примерах синтеза двух САУ. Она проста для понимания и может быть рекомендована для использования студентами при курсовом и дипломном проектировании.
2
В в е д е н и е
С и н т е з с и с т е м ы а в т о м а т и ч е с к о г о у п р а в л е н и я ( С А У ) —
это определение структуры, состава элементов, и значений параметров САУ, при которых она удовлетворяет предъявляемым к ней требованиям.
Данное определение характеризует общую задачу синтеза.
Задача определения типа, схемы и параметров закона управления (ЗУ)
или корректирующих устройств (КУ) называется ч а с т н о й з а д а ч е й с и н т е з а С А У .
Разработано большое количество методов решения частной задачи синтеза САУ, в частности, в биотехнологии в настоящее время использу-
ются два метода: графоаналитический [1,2] и интегральных оценок [3].
Ниже для решения задачи синтеза САУ биотехнологическими объектами управления (ОУ) предлагается использовать метод стандартных коэффи-
циентов [4].
3
1. Постановка задачи синтеза САУ методом
стандартных коэффициентов
Большинство (ОУ) в биотехнологии может быть представлено следую-
щей структурной схемой [5]:
Регулирующий |
|
|
|
|
|
|
орган |
|
|
F(p) |
|
||
G(p) |
|
|
|
|
|
Y(p) |
|
|
|
|
k |
||
k=1 |
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
( T1 p 1 )( T2 p 1 ) |
|
|
|
|
|
|||
|
|
|
|
|
|
|
Параметры передаточных функций ОУ, как правило, получаются в ре-
зультате идентификации на основании экспериментально полученной кри-
вой разгона. Задающее воздействие и возмущение принимаются постоян-
ными, регулирующий орган – безынерционным (его передаточная функция
равна 1). |
|
|
Т р е б у е т с я с и н т е з и р о в а т ь |
САУ, удовлетворяющую следующим |
|
требованиям: время регулирования tр |
tр зад , перерегулирование |
зад, |
установившаяся ошибка САР при воздействии задающего сигнала |
и воз- |
|
мущения должна бать равна нулю; должны быть ограничены так же мак-
симальные значения выходного сигнала и отклонения регулирующего орга-
на при воздействии возмущения.
В методе стандартных коэффициентов структура и параметры регулято-
ра: k0 , k1, k2, … должны быть выбраны так, чтобы передаточная функция замкнутой САУ равнялась желаемой (стандартной передаточной функции
или передаточной функции эталонной модели - ЭМ):
|
WСАУ ( p, k0 , k1 ,...) WЖ ( p ) WЭМ ( p ) |
(1) |
||||||||
или |
BСАУ ( p, k0 |
, k1 ,...) |
|
BЖ ( p ) |
|
BЭМ ( p ) |
. |
|||
|
|
|
|
|
|
|||||
AСАУ ( p, k0 |
, k1 ,...) |
AЖ ( p ) |
AЭМ ( p ) |
|||||||
|
|
|
||||||||
|
|
|
|
|||||||
Для решения полиномиального уравнения (1) обычно выписывают сис-
тему уравнений для коэффициентов при одинаковых степенях p. Эта зада-
4
ча будет иметь решение, если структура замкнутой САУ и ЭМ выбраны так, чтобы выполнялись условия «соответствия передаточных функций»
(СПФ):
1) совпадали размерности полиномов числителей и знаменателей пере-
даточных функций замкнутой САУ и ЭМ:
dim BСАУ ( p, k0 , k1 ,...) dim BЭМ ( p ) , |
(2) |
dim АСАУ ( p, k0 , k1 ,...) dim АЭМ ( p ) ; |
|
2) число n искомых параметров регулятора равнялось числу уравнений,
получаемых из уравнения (1):
dim AЭМ ( p ) dim BЭМ ( p ) n . |
(3) |
Выбор структуры замкнутой САУ должен обеспечивать также выпол-
нение требований качества и устойчивости. В рассматриваемой задаче синтеза замкнутая САУ должна быть астатической первого порядка по от-
ношению к задающему воздействию и возмущению. Обеспечение астатиз-
ма достигается использованием И, ПИ, ПИД – регуляторов. Для обеспе-
чения заданного времени регулирования необходимо использовать после-
довательно включенный П – регулятор.
В работе выбраны две пары структур замкнутой САУ и ЭМ, отвечаю-
щих перечисленным выше требованиям. Рассмотрим особенности синтеза этих САУ методом стандартных коэффициентов и их свойства.
5
2. Синтез САУ №1
Структурная схема САУ №1 представлена на рис.1.
|
|
|
|
|
|
|
|
|
|
F(p) |
||||
G(p) |
|
(p) |
|
|
|
U(p) |
|
|
|
|
|
|
Y(p) |
|
|
|
|
|
|
k |
|||||||||
|
|
|
k1 |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
( T1 p 1 )( T2 p 1 ) |
|
||||||
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k0
p
1 k2 p
Рис. 1. Структурная схема САУ №1
Р А С Ч Е Т П Е Р Е Д А Т О Ч Н О Й Ф У Н К Ц И И С А У № 1
Передаточная функция замкнутой САУ
|
|
|
( k0 k1 p )k |
|
|
|
|
WСАУ ( p ) |
|
p( T1 p 1 )( T2 p 1 ) |
|
|
|||
1 |
( k0 k1 p )( 1 k2 p )k |
|
|||||
|
|
|
|||||
|
p( T p |
1 )( T p 1 ) |
|
|
|||
|
|
|
1 |
2 |
|
|
|
|
|
|
|
|
|
|
|
( k0 |
k1 p )k |
|
|
|
|
|
|
|
|
|||||||||
T T p3 |
( T |
T ) p2 |
p k k k |
2 |
p2 ( k k |
2 |
k )k p k k |
|||||||||||||||||||
|
|
|
|
|||||||||||||||||||||||
1 |
2 |
|
|
1 |
2 |
|
|
|
|
1 |
|
|
|
0 |
1 |
0 |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
( |
k1 |
|
|
p 1 ) |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
k0 |
|
|
|
|
|
|
|
. |
(4) |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
T T |
|
|
p3 |
|
( T |
T k k |
k ) |
p2 |
[( k k |
2 |
k |
|
)k 1] |
|
||||||||||||
|
|
1 2 |
|
|
1 |
2 |
1 2 |
|
|
|
|
|
0 |
1 |
|
|
p 1 |
|
||||||||
|
|
k0k |
|
|
|
k0k |
|
|
|
|
|
|
|
|
k0k |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Как видно, в знаменателе передаточной функции САУ (4) необходимо задавать три коэффициента, для этого введены три передаточных числа. Однако, что делать с параметрами числителя ПФ (2)?
6
В Ы Б О Р Э Т А Л О Н Н О Й М О Д Е Л И ( Э М )
Выбор эталонной модели (ЭМ) осуществляется исходя из требований к качеству замкнутой системы. Так как разомкнутая система имеет порядок равный двум, ЭМ выберем в виде инерционного звена 2-го порядка
|
Wэм |
( p ) |
|
|
|
|
k ЭМ |
|
|
|
|
|
|
|
|
|
|
|
k ЭМ |
|
|
, |
|
||
|
T |
2 |
p 2 2T |
|
d |
|
|
|
|
|
|
|
|
p |
2 a |
|
|
|
|||||||
|
|
|
|
ЭМ |
ЭМ |
p 1 a |
M 2 |
M 1 |
p 1 |
|
|||||||||||||||
|
|
|
|
|
ЭМ |
|
|
|
|
|
|
|
|
|
|
|
|
(5) |
|||||||
|
|
|
T 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a |
M 2 |
|
; |
|
a |
M 1 |
2T |
ЭМ |
d |
ЭМ |
; |
|
|
k |
ЭМ |
1. |
|
|
|
|
||||
|
|
ЭМ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
с |
оптимальным |
|
|
значением |
относительного |
|
коэффициента |
затухания |
|||||||||||||||||
dэм |
= 0.7. |
В этом |
случае |
время |
регулирования tР |
3 Tэм , |
а пере- |
||||||||||||||||||
регулирование < 5%. Постоянную времени ЭМ выбираем из условия:
TЭМ tР / 3 . |
(6) |
С целью обеспечения требования астатизма ЭМ по отношению к задаю-
щему воздействию выбираем
kЭМ = kСАУ = 1.
Приведем структуру передаточной функции ЭМ (5) к виду (4), помножив ее числитель и знаменатель на множитель ( k1 / k0 p 1) :
|
|
|
|
|
( |
k1 |
p 1 ) |
|
|
|
|
|
|
|
|
( |
k1 |
p 1 ) |
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
Wэм( p ) |
|
|
|
|
|
|
|
k0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k0 |
|
|
|
|
|
|
|
|||
|
|
|
|
p2 a |
|
|
|
|
k1 |
|
|
|
|
|
|
p2 a |
|
|
|
|
|
|
k1 |
|
||||||||||
|
( a |
|
|
|
|
|
p 1 )( |
|
p 1 ) ( a |
|
|
|
|
|
p 1 )( |
p 1 ) |
||||||||||||||||||
M 2 |
M 1 |
k0 |
M 2 |
M 1 |
k0 |
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
( |
k1 |
p 1 ) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
k0 |
|
|
|
|
|
|
|
|
|
|
|
. |
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
aM 2k1 |
p3 ( |
aM 1k1 |
a |
|
) p2 |
( a |
|
|
k1 |
|
) p 1 |
|
|
|
||||||||||||||||||
|
|
|
M 2 |
M 1 |
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
k0 |
|
|
|
k0 |
|
|
|
|
|
|
|
k0 |
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
Сравнивая полученную передаточную функцию с передаточной функ-
цией (4), приравниваем коэффициенты знаменателей при одинаковых сте-
пенях р:
|
T1T2 |
|
aM 2k1 |
; |
|
|
||
|
|
|
|
|
|
|||
|
k0k |
|
|
k0 |
|
|||
( T1 T2 k1k2k ) |
|
aM 1k1 |
a |
; |
||||
|
|
|||||||
k0k |
|
|
|
|
|
M 2 |
|
|
|
|
|
k0 |
|
||||
7
[( k0k2 k1 )k 1] |
a |
|
|
k1 |
|
|
M 1 |
k . |
|||
k k |
|
||||
|
|
||||
0 |
|
|
0 |
|
|
Отсюда последовательно определяем:
k1 T1T2 ; aM 2k
T1 T2 k1k2 k aM 1k1k aM 2 k0 k; k0 k2 k 1 aM 1k0 k;
k |
|
|
T1 T2 aM 1k1k |
|
aM 2 |
k |
|
. |
2 |
|
|
0 |
|||||
|
|
k1k |
|
k1 |
|
|||
|
|
|
|
|
|
|||
Подставляя k2 в предпоследнее выражение, получим:
|
aM 2 |
k k |
|
2 ( k |
|
T1 T2 aM 1k1k |
a |
|
|
k )k |
|
1 |
0; |
|||||||||||||||
|
|
|
|
|
M 1 |
0 |
||||||||||||||||||||||
|
|
k1 |
|
|
0 |
|
|
|
|
|
k1k |
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
a |
M 2 |
k k 2 |
( T |
T |
2 |
)k |
0 |
k |
1 |
0. |
|
|
|
|
|
|
|
|
||||||||||
|
|
0 |
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
Отсюда получаем два решения: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
( T |
T ) |
|
( T |
T )2 |
4a |
|
k k |
|
|
|
|
|
|||||||||
|
|
|
|
k0 |
|
|
1 |
|
|
2 |
|
|
|
|
1 |
|
|
2 |
|
M 2 |
|
1 |
|
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2aM 2k |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
12 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
M A P L E - П Р О Г Р А М М А Р А С Ч Е Т А П Е Р Е Д А Т О Ч Н Ы Х Ч И С Е Л Р Е Г У Л Я Т О Р А
Ввод параметров передаточной функции ОУ:
W(p)=k/[(T1p+)(T2p+1)]
> T1:=20.0; T2:=40.0; k:=10.0;
T1 := 20.0
T2 := 40.0
k := 10.0
Ввод заданного времени регулирования
> tр:=60.0;
|
tð := 60.0 |
|
Определение параметров передаточной функции ЭМ: |
||
|
Wэм(p)=1/(Am2*p*p+Am1*p+1), |
|
где Tm = tР /3; |
Am2=Tm*Tm; |
Am1=2*Tm*d; |
(7)
(8)
(9)
8
> Tm:=tр/3; d:=0.7;
Tm := 20.00000000
d := .7
> Am2:=Tm*Tm; |
Am1:=2*Tm*d; |
Am2 := 400.0000000
Am1 := 28.00000000
Определение передаточных чисел регулятора:
П е р в о е р е ш е н и е :
> k1:=T1*T2/Am2/k;
k1 := .2000000000
> k0:=((T1+T2)+sqrt((T1+T2)^2-4*Am2*k1*k))/(2*Am2*k);
k0 := .01000000000
> k2:=(-T1-T2+Am1*k1*k)/k1/k+Am2*k0/k1;
k2 := 18.00000000
Определение передаточных чисел регулятора:
Вт о р о е р е ш е н и е :
>k0:=((T1+T2)-sqrt((T1+T2)^2-4*Am2*k1*k))/(2*Am2*k);
k0 := .005000000000
> k2:=(-T1-T2+Am1*k1*k)/k1/k+Am2*k0/k1;
k2 := 8.000000000
9
h(t)
t, c
а) 1 и 2-й вар. регулятора: g(t) = 1(t), f(t) = 0, h( ) = 1, |
tр = 56 с, = 4,3%. |
h(t)
t, c
б) 1-й вар. регулятора: g(t) = 0, f(t) = 0,5, h( ) = 0, hmax = 0,85
10
