14_08_18_ТАУ_1,2_Лекционный курс
.pdf
|
Q |
|
0.005 |
|
м3 |
|
|
KQ = |
max |
= |
|
= 0.0127, |
|
. |
(10.8) |
|
|
|
|||||
ϕ p |
|
0.392 |
|
c рад |
|
||
|
max |
|
|
|
|
|
|
Математическая зависимость между высотой воды в баке H(T ) и ее расходом также имеет вид интегральной зависимости:
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
||||||
h(t) = Kh ∫ Qdt = KhKQ ∫ ϕ p(t)dt , |
|
|
(10.9) |
|||||||||||||||||||||||||
|
|
0 |
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
||||||||||
1 |
1 |
|
|
= 1 м−2. |
|
|
|
|
|
|
||||||||||||||||||
где Kh = |
|
|
= |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
a b 1×1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
Передаточная функция этого звена: |
|
|
|
|||||||||||||||||||||||||
Wh(s) |
|
|
|
|
h(s) |
|
|
|
|
|
KhKQ |
|
|
|
|
|
|
|||||||||||
= |
|
|
|
|
|
|
|
|
= |
|
|
|
|
|
|
. |
|
|
|
|
|
(10.10) |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
ϕ р(s) |
|
|
|
|
|
|
|
s |
|
|
|
|
|
|
|||||||||
Рассмотрим процесс заполнения бака водой до высоты hз =1 м . В этом |
||||||||||||||||||||||||||||
случае передаточные функции системы имеют вид. |
|
|||||||||||||||||||||||||||
• Задание частоты вращения двигателя: |
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
Uз |
10 |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
Wзад(s) = |
|
|
= |
|
|
|
=10. |
|
|
|
|
|
|
(10.11) |
||||||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
hз |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
• Передаточная функция преобразователя и двигателя: |
|
|||||||||||||||||||||||||||
Wд(s) |
|
|
|
|
ωд(s) |
|
|
= Kд = 6.28. |
|
|
(10.12) |
|||||||||||||||||
= |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
Uз(s) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
• Передаточная функция задвижки: |
|
|
|
|||||||||||||||||||||||||
|
|
|
ϕ р(s) |
|
|
|
|
|
Kз |
0.0006 |
|
|
|
|
|
|||||||||||||
Wз(s) = |
|
|
= |
|
|
|
|
= |
|
|
|
|
. |
|
|
|
(10.13) |
|||||||||||
|
|
|
s |
|
|
|
s |
|
|
|
||||||||||||||||||
|
|
|
|
|
ωд(s) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
• Передаточная функция бака: |
|
|
|
|||||||||||||||||||||||||
Wh(s) |
|
|
|
|
h(s) |
|
|
|
|
|
KhKQ |
1× 0.0127 |
|
0.0127 |
. |
(10.14) |
||||||||||||
= |
|
|
|
|
|
|
|
|
= |
|
|
|
|
= |
|
= |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
ϕ р(s) |
|
s |
|
s |
s |
|
||||||||||||||||
Структурная схема |
|
исполнительного |
механизма |
представлена на |
||||||||||||||||||||||||
рис. 10.1. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
101
|
Wзад(s) |
Wд(s) |
|
Wз (s) |
||
h (s) |
|
|
|
|
|
ϕ р (s) |
|
|
|
|
|
||
з |
|
|
0.0006 |
|
|
|
|
10 |
6.28 |
|
|
||
|
|
|
|
|
||
|
|
S |
|
|||
|
|
|
|
|
||
Задание |
Задание |
|
|
|
|
|
Двигатель |
Задвижка |
|||||
уровня воды |
скорост и |
|
|
|
|
|
Рисунок 10.1 - Структурная схема исполнительного механизма Эквивалентная передаточная функция исполнительного механизма
(рис. 10.1):
Wem ( s) = |
ϕ р (s) |
0.03768 |
|
|
|
|
= Wзад(s)Wд(s)Wз (s) = |
|
. |
(10.15) |
|
|
|
||||
|
hз (s) |
s |
|
На рис. 10.2 изображена структурная схема исполнительного механизма и
бака.
|
|
WEM (S) |
|
WH(S) |
|||||
hз (s) |
|
|
|
|
ϕ р (s) |
|
|
|
H(S) |
|
0.03768 |
|
0.0127 |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
S |
|
|
|
|
|
|
|
|
|
|
|
S |
|
|||
Задание |
|
|
|
|
|
|
|
|
|
Исполнительный |
|
Бак |
|||||||
уровня воды h |
механизм |
|
|
|
|
Рисунок 10.2 - Структурная схема исполнительного механизма и бака Если выполнить систему управления заполнением бака водой без
регуляторов, то она будет неработоспособна, поскольку содержит цепочку последовательно соединенных интегрирующих звеньев (рис.10.3).
hз (s) |
|
|
|
|
|
|
WEM (S) |
|
WH(S) |
|||||||
|
|
|
|
|
|
|
|
|
ϕ р (s) |
|
|
|
|
H(S) |
||
|
|
|
|
|
|
0.03768 |
|
0.0127 |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
S |
|
|
|
|
|
|
|
Задание |
|
|
|
|
|
|
|
|
|
|
S |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Исполнительный |
|
Бак |
|
|||||||
уровня воды h |
|
|
|
|
|
|
механизм |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рисунок 10.3 - Структурная схема системы заполнения бака водой при полностью открытой задвижке без регуляторов
Каждое звено сдвигает фазу на девяносто градусов, что в сумме дает сто восемьдесят градусов. Сдвиг фазы на сто восемьдесят градусов в замкнутой системе равнозначен введению положительной обратной связи. Структурная
102
схема системы управления заполнением бака водой без регулятора и переходной процесс в системе показаны соответственно на рис. 10.3 и 10.4.
h, |
|
|
|
|
|
|
2.5 |
|
|
|
|
|
|
м |
|
|
|
|
|
|
2 |
|
|
|
|
|
|
1.5 |
|
|
|
|
|
|
1 |
|
TN |
|
|
|
|
|
|
|
|
|
|
|
0.5 |
|
|
|
|
|
|
0 |
|
|
|
|
|
|
0 |
100 |
200 |
300 |
400 |
500 |
600 |
|
|
|
|
|
|
t, с |
Рисунок 10.4 - Переходной процесс в системе В системе регулирования возникают колебания с перерегулированием
100%. Период колебаний определим из передаточной функции замкнутой системы (рис. 10.3):
|
|
|
|
|
|
|
|
|
0.03768 |
|
|
0.0127 |
|
|
|
|||||
|
з |
|
|
|
Wem ( s)Wh ( s) |
|
|
|
|
× |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
W |
( s) = |
|
|
= |
|
|
s |
|
|
s |
|
= |
||||||||
|
1+Wem ( s)Wh ( s) |
|
|
0.03768 |
|
|
|
0.0127 |
|
|||||||||||
|
|
|
1+ |
|
× |
|
|
|||||||||||||
|
|
|
s |
s |
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
= |
|
|
|
|
|
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
(10.16) |
||||
|
|
|
S |
2 +1 |
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
0.000479 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Приравняв полученное выражение для передаточной функции замкнутой системы (9.16) к передаточной функции колебательного звена (9.14), получим:
|
1 |
|
|
|
|
|
|
|
1 |
|
|
. |
|
|
|
|
|
= |
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
1 |
|
|
2ξ |
|
|
|||
|
|
2 +1 |
|
2 |
|
|
|
(10.17) |
|||||
|
|
s |
|
+ |
s +1 |
||||||||
|
|
|
|
|
|||||||||
|
0.000479 |
|
2 |
s |
|
|
|||||||
|
|
|
|
|
ω |
|
|
ω |
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
пп
Из уравнения (10.17) следует, что в системе существуют незатухающие колебания (демпфирование колебаний ξ = 0) с угловой частотой
ωn = 0.000479 = 0.0218 с−1.
Период незатухающих колебаний (рис. 10.4):
103
TN = |
2π |
2 ×3.14 |
= 288.1 с. |
|
|
|
= |
|
(10.18) |
||
|
0.0218 |
||||
|
ωN |
|
|
Постановка задачи оптимизации. При заданной схеме исполнительного механизма и бака (рис. 10.2) построить систему автоматического заполнения бака водой со следующими показателями качества:
•перерегулирование отсутствует,
•время полного открытия задвижки не должно превышать заданного значения TO = 10 с,
•время переходного процесса заполнения бака водой должно быть минимальным.
Выполним систему управления заполнением бака водой многоуровневой (рис. 10.5). Внутренний контур будет управлять положением задвижки, а внешний - заполнением бака.
hз (s) |
WEM ( S) |
W (S) |
|
|
|
ϕ p (s) |
H |
H |
S |
) |
|
|
|||||
К рh |
0.0127 |
( |
|
||
0.03768 |
|
|
|
||
|
|
|
|
||
|
К рb |
S |
|
|
|
|
S |
|
|
|
|
|
|
|
|
|
Рисунок 10.5 - Структурная схема системы заполнения бака водой при полностью открытой задвижке
10.2.1. Расчет параметров регулятора и моделирование переходных процессов в контуре управления положением задвижки
В контуре управления положением задвижки применим регулятор в виде пропорционального звена, поскольку задвижка представлена интегрирующим звеном и в замкнутой системе регулирования статическая ошибка будет равна нулю. Структурная схема модели системы изображена на рис. 10.6.
104
|
|
|
|
|
|
|
|
|
|
|
WEM ( S) |
|
ϕ р ( s) |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ϕ рз ( s) |
|
|
|
|
|
|
|
К рb |
|
0.03768 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
S |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
CONSTANT |
|
GAIN |
TRANSFER FCN |
|
|
SCOPE |
|||||||||||
Угол открытия |
Регулятор |
Исполнительный |
|
|
|
|
|||||||||||
задвижки |
|
механизм |
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рисунок 10.6 - Структурная схема модели контура управления положением задвижки
Главная задача при выборе регулятора в данном случае - обеспечить требуемое время регулирования (TO = 10 с) при отсутствии перерегулирования. Контур регулирования положения задвижки (рис. 10.6) содержит интегратор, охваченный единичной обратной связью, и представляет собой апериодическое звено. Как известно, переходные процессы в апериодическом звене протекают без перерегулирования, ЛАЧХ апериодического звена имеет наклон
дек в области частоты среза.
Поскольку в замкнутой системе управления положением задвижки колебания отсутствуют, то задача выбора регулятора сводится к определению
значения К pb , при котором время регулирования составит TO = 10 с.
Для решения данной задачи воспользуемся приложением SISO Tool пакета Matlab.
В командной строке пакета Matlab зададим передаточную функцию
исполнительного механизма |
|
WEM=TF([0.03768],[1 0]), |
(10.19) |
и, предварительно, регулятора
Kpb=tf(1). (10.20) Дальше в командную строку введем команду Sisotool.
Для ввода передаточных функций исполнительного механизма и регулятора выполним пункт меню File→Import. Откроется окно (рис. 10.7) вкладки Import System Data.
105
Рисунок 10.7 - Окно программы Import System Data
Выбираем в окне SISOModels передаточную функцию исполнительного механизма WEM и с помощью кнопки со стрелкой переносим ее в окно с передаточной функцией объекта управления G. То же самое делаем для регулятора задвижки, перенеся передаточную функцию Kpb в окно с передаточной функцией регулятора C .
После задания передаточных функций всех элементов структурной схемы системы нажимаем кнопку OK и переходим к окну программы SISO Design Tool (рис. 10.8).
106
Рисунок 10.8 - Формирование желаемой ЛАЧХ контура управления положением задвижки
Для контроля длительности процесса открытия задвижки откроем окно Step Respons для LTI Viewer for SISO Design Task.
В окне Open-Loop Bode Editor, смещая ЛАЧХ, корректируем передаточную функцию разомкнутого контура так, чтобы длительность переходного процесса (время регулирования) открытия задвижки (рис. 10.9) не превышала TO = 10 с.
107
Рисунок 10.9 - Переходная функция контура управления положением задвижки
При выполнении этого условия получен коэффициент передачи
регулятора K рb = 7.95 (рис. 10.8). Структурная схема модели системы
управления положением задвижки с рассчитанным коэффициентом усиления изображена на рис. 10.10.
|
|
|
|
|
|
|
|
К рb |
|
|
WEM ( S) |
|
ϕ р ( s) |
||||
ϕ рз ( s) |
|
|
|
|
|
|
|
7.95 |
|
0.03768 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
S |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
CONSTANT |
|
GAIN |
TRANSFER FCN |
|
|
SCOPE |
|||||||||||
Угол открытия |
Регулятор |
Исполнительный |
|
|
|
|
|||||||||||
задвижки |
|
механизм |
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рисунок 10.10 - Структурная схема модели оптимизированной системы управления положением задвижки
108
|
Промоделируем систему (рис. 10.10), подав на вход значение угла, |
||||||||
соответствующее |
максимальному |
открытию |
|
задвижки |
|||||
ϕ рз |
= ϕ р |
= 0.392 рад. |
График |
переходного |
процесса |
в |
контуре |
||
|
max |
|
|
|
|
|
|
|
|
управления положением задвижки представлен на рис. 10.11. |
|
|
|||||||
|
|
ϕ |
рз = 0.392 |
5% |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
ϕ р, 0.4 |
|
|
|
|
|
|
|
|
|
−1 |
|
|
|
|
|
|
|
|
|
с |
|
|
|
|
|
|
|
|
|
|
0.3 |
|
|
|
|
|
|
|
|
|
0.2 |
|
|
|
|
|
|
|
|
|
0.1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
t р |
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
0 |
|
5 |
|
10 |
15 |
t, |
с |
Рисунок 10.11 - Переходной процесс управления положением задвижки Как видно из графика, время регулирования составляет 10 с, что
соответствует требуемому времени полного открытия задвижки TO = 10 с. Преобразуем структурную схему рис. 10.10, определив передаточную
функцию замкнутого контура управления положением задвижки:
з |
K pb(s)Wem(s) |
K pbWem(s) |
|
||
Wem(s) = |
|
= |
|
. |
(10.21) |
|
|
||||
|
1+ K pb(s)Wem(s) 1+ K pbWem(s) |
|
|
|
7.95× |
|
0.03768 |
|
|
|
|
|
||
з |
|
s |
|
|
|
1 |
|
|
|||
|
|
|
|
|
|
|
|||||
Wem(s) = |
|
|
|
|
|
|
|
= |
|
. |
(10.22) |
|
7.95 |
× 0.03768 |
|
|
|||||||
1 |
|
|
3.34s +1 |
|
|||||||
+ |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
s
Получим также эквивалентную передаточную функцию объекта управления высотой воды в баке: последовательное соединение звеньев
|
|
|
|
з |
|
|
|
|
исполнительного механизма Wem(s) и бака WH(S) : |
|
|||||||
|
з |
|
|
|
|
|
(10.23) |
|
Wоh (s) = Wem (s)Wh (s) , |
|
|
|
|||||
1 |
|
0.0127 |
|
0.0127 |
|
|||
Wоh (s) = |
|
|
× |
|
= |
|
. |
(10.24) |
3.34s +1 |
s |
2 |
||||||
|
|
|
3.34s + s |
|
||||
|
|
|
|
|
|
|
109
10.2.2.Расчет параметров регулятора и моделирование переходных процессов в контуре управления высотой воды в баке
Вконтуре управления высотой воды в баке применим регулятор в виде пропорционального звена, поскольку передаточная функция бака описывается интегрирующим звеном и в замкнутой системе регулирования статическая ошибка будет равна нулю. Структурная схема модели системы изображена на рис. 10.12.
|
|
|
|
|
|
|
|
|
|
|
WOH (S) |
|
|
|
H(S) |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
hз (s) |
|
|
|
|
|
|
К рh |
|
0.0127 |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
3.34S |
+S |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
CONSTANT |
|
|
|
|
|
GAIN |
TRANSFER FCN |
|
|
|
SCOPE |
||||||
Задание |
|
|
|
|
|
Регулятор |
|
|
Бак |
|
|
|
|
|
|
|
||
уровня воды |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рисунок 10.12 - Структурная схема модели системы управления высотой воды в баке
Задачу формирования желаемой ЛАЧХ формулируем следующим образом: обеспечить максимально широкую полосу пропускания (максимальное быстродействие) при отсутствии колебаний уровня воды.
В пакете Matlab зададим переходную функцию объекта управления высотой воды в баке, введя в командную строку
WOH=TF([0.0127],[3.34 1 0]) , |
(10.25) |
и, предварительно, регулятора |
|
Kph=tf(1). |
(10.26) |
Вкомандную строку введем команду Sisotool. Для контроля перерегулирования и длительности переходного процесса откроем окно Step Response для LTI Viewer for SISO Design Task.
Вокне Open-Loop Bode Editor (рис. 10.13) корректируем ЛАЧХ разомкнутого контура так, чтобы обеспечить максимальное быстродействие при отсутствии перерегулирования (рис. 10.14).
110