Учебное пособие 800455
.pdfОкончание табл. 16
8 |
1 |
0,001 |
10 |
1,8 |
0,55 |
9 |
1 |
0,001 |
10 |
1,9 |
0,6 |
10 |
1 |
0,001 |
10 |
2,1 |
0,65 |
11 |
1 |
0,001 |
10 |
2,2 |
0,7 |
12 |
1 |
0,001 |
10 |
2,3 |
0,75 |
13 |
1 |
0,001 |
10 |
2,4 |
0,8 |
14 |
1 |
0,001 |
10 |
2,5 |
0,85 |
15 |
1 |
0,001 |
10 |
2,6 |
0,9 |
КОНТРОЛЬНЫЕ ВОПРОСЫ
1.Составные участки ЖЛАЧХ.
2.Свойства ЖЛАЧХ при одновременном изменении постоянных времени.
3.ЖЛАЧХ статической системы.
4.ЖЛАЧХ астатической системы.
5.Определение параметров качества регулирования.
6.Схема, реализующая ЖЛАЧХ статической системы.
7.Схема, реализующая ЖЛАЧХ астатической системы.
ЛАБОРАТОРНАЯ РАБОТА № 12
ИССЛЕДОВАНИЕ АДЕКВАТНЫХ МОДЕЛЕЙ АВТОМАТИЧЕСКОЙ СИСТЕМЫ
Цель работы: получение и исследование адекватных моделей системы автоматического управления.
ТЕОРЕТИЧЕСКИЕ СВЕДЕНИЯ
Синтез осуществляется в пространстве состояний.
В систему автоматического управления вводится модальный регулятор, на вход которого подаются соответству-
81
ющие переменные (выходные сигналы в структурной схеме системы) с коэффициентом обратной связи.
Для расчета в среде Matlab нужно иметь навыки получения моделей:
–структурная схема в системе Simulink на основе использования блоков Transfer Fnc;
–развернутая структурная схема в системе Simulink на основе использования блоков Integrator;
–lti-объект в рабочей области среды Matlab, представленный системой пространства состояний объекта управления.
ПРЕДВАРИТЕЛЬНОЕ ЗАДАНИЕ
1.Используется дифференциальное уравнение 2-го порядка для системы автоматического управления.
2.Cтатическая система автоматического управления (САУ) в составе предварительного усилителя – (пропорциональное звено с коэффициентом Kpu), усилителя мощности (апериодическое звено с коэффициентом Ку и постоянной времени Ту), двигателя (апериодическое звено с коэффициентом Kd и постоянной времени Td=Tm), датчика обратной связи (пропорциональное звено с коэффициентом Kdos) – рис. 33.
МЕТОДИЧЕСКИЕ УКАЗАНИЯ К ВЫПОЛНЕНИЮ ЛАБОРАТОРНОЙ РАБОТЫ
1.Параметры САУ представлены в табл. 17.
2.Структурная схема статической САУ – рис. 33.
|
|
Ky |
Kd |
1 |
1 |
Kpu |
|
|
|
Ty.s+1 |
Tds+1 |
|
||
|
|
Out1 |
||
In1 |
|
|
|
|
Gain |
Transfer Fcn |
Transfer Fcn1 |
|
|
|
|
Kdos
Gain1
Рис. 33. Структурная схема статической САУ
82
2.1.Для исследований используется в рабочей области Matlab следующая последовательность команд:
»[a,b,c,d] =linmod(‘имя файла’);
»h=ss(a,b,c,d); »step(h).
2.2.Получить, изменяя Kpu, переходную характеристику при перерегулировании 5% (5.01%) – колебательную и при перерегулировании 0.05% – апериодическую на границе с колебательной.
2.3.Определить параметры качества регулирования
по переходной характеристике: время регулирования; перерегулирование и установившееся значение позиционной ошибки.
3. Исследовать структурную схему САУ с доступом к переменным состояния – рис. 34.
1 |
Kpu |
Ky |
1 |
|
s |
||
In1 |
|
Ty |
|
|
|
||
Gain |
|
Integrator |
|
|
Transfer Fcn2 |
||
|
|
|
|
|
|
|
1/Ty |
|
|
|
Gain2 |
Kd |
1 |
1 |
|
|
s |
||
Td |
Out1 |
||
|
|||
|
Integrator1 |
||
Transfer Fcn3 |
|
||
|
|
||
|
1/Ty |
|
|
|
Gain3 |
|
|
Kdos |
|
|
|
Gain1 |
|
|
Рис. 34. Структурная схема САУ с доступом к переменным состояния
Таблица 17
Параметры автоматической системы
Вариант |
Предвар. |
Усилитель |
Двигатель |
Датчик |
||
усилитель |
мощности |
|
|
ОС |
||
|
|
|
|
|
|
|
Кpu |
Кy |
Тy |
Kd |
Td |
Кdos |
|
1 |
2 |
11 |
0.0011 |
2 |
0.011 |
0.12 |
2 |
2 |
11 |
0.0012 |
2 |
0.01 |
0.12 |
3 |
2 |
11 |
0.0013 |
2 |
0.01 |
0.12 |
4 |
2 |
11 |
0.0014 |
2 |
0.01 |
0.12 |
5 |
2 |
11 |
0.0015 |
2 |
0.01 |
0.12 |
6 |
3 |
12 |
0.0021 |
3 |
0.011 |
0.11 |
83
Окончание табл. 17
7 |
3 |
12 |
0.0022 |
3 |
0.011 |
0.11 |
8 |
3 |
12 |
0.0023 |
3 |
0.011 |
0.11 |
9 |
3 |
12 |
0.0024 |
3 |
0.011 |
0.11 |
10 |
3 |
12 |
0.0025 |
4 |
0.011 |
0.11 |
11 |
1 |
10 |
0.0016 |
4 |
0.012 |
0.13 |
12 |
1 |
10 |
0.0017 |
4 |
0.012 |
0.13 |
13 |
1 |
10 |
0.0018 |
4 |
0.012 |
0.13 |
14 |
1 |
10 |
0.0019 |
4 |
0.011 |
0.13 |
15 |
1 |
10 |
0.002 |
4 |
0.012 |
0.13 |
3.1.Ввести поочередно полученные значения Kpu
(см. пп. 2.2) для переходных характеристик (колебательная и апериодическая 2 – го порядка) от предыдущей схемы.
3.2.Определить параметры качества регулирования
по переходной характеристике: время регулирования; перерегулирование и установившееся значение ошибки.
4. Получить векторно-матричное уравнение по структурной схеме с доступом к переменным состояния.
4.1.Ввести команды в рабочую область Matlab:
»a1=[ ]; »b1=[ ]; »c1=[ ]; »d1=[0];
»h1=ss(a,b,c,d); »step(h1).
4.2. Определить параметры качества регулирования
по переходной характеристике: время регулирования; перерегулирование и установившееся значение ошибки.
5. Исследование передаточной функции системы в zpk –
форме.
5.1. Ввести в рабочую область Matlab команду получения zpk – формы передаточной функции структурной схемы
»h2=zpk(h1);
84
»step(h2).
5.2. Определить параметры качества регулирования
по переходной характеристике: время регулирования; перерегулирование; установившее значение позиционной ошибки.
6. Сделать выводы по работе об адекватности схем.
КОНТРОЛЬНЫЕ ВОПРОСЫ
1.Структурная схема САУ в области Simulink.
2.Структурная схема САУ с использованием интегрирующих звеньев.
3.Представление структурной схемы САУ в lti-объекте.
4.Представление простейших САУ с помощью интегрирующих звеньев.
5.Структурная схема с доступом к переменным состояния апериодического типового динамического звена 1-го порядка.
6.Структурная схема с доступом к переменным состояния апериодического типового динамического звена 2-го порядка.
7.Структурная схема с доступом к переменным состояния колебательного типового динамического звена.
ЛАБОРАТОРНАЯ РАБОТА № 13
ИССЛЕДОВАНИЕ АВТОМАТИЧЕСКОЙ СИСТЕМЫ С МОДАЛЬНЫМ РЕГУЛЯТОРОМ
Цель работы: исследование системы автоматического регулирования скорости с модальным регулятором.
ТЕОРЕТИЧЕСКИЕ СВЕДЕНИЯ
При моделировании в среде программы Matlab (переменная s = p), математические модели представлены tf – формой или zpk – формой.
85
Для tf – формы ММ представляется в виде передаточной функции, где числитель и знаменатель это полиномы, например, 2-го порядка
W (s) = (b2 s² + b1 s + bo) / (а2 s² + а1 s + аo),
где bi, аi – коэффициенты полиномов; s – переменная.
В этом случае дифференциальное уравнение 1-го и 2-го порядка стандартных форм представляется в общем виде
а2*d²y(t)/dt² + а1*dy(t)/dt + аo*y(t) = b2*d²x(t)/dt² +
+b1*dx(t)/dt + bо*x(t).
Для zpk – формы ММ уравнения n-го порядка представляется в виде нулей, полюсов и обобщенного коэффициента передачи
H (s) = K*(s-z1)*(s – z2) …(s - zm) / (s – p1)*(s – p2)…(s – pn),
где zi, pi и К – нуль, полюс и обобщенный коэффициент передачи.
Третьей формой МО элементов является ss-форма, то есть векторно-матричное уравнение.
Векторно-матричное уравнение стационарной системы запишется:
dx/dt = A*x + B*u y = C*x + D*u.
А – матрица размера n x n (квадратная матрица n-го порядка), называемая матрицей системы. Она характеризует динамические свойства системы;
В – прямоугольная n x r матрица, называемая матрицей управления. Она характеризует воздействие входных переменных uj на переменные состояния xi;
86
С – m x n матрица измерения. Она характеризует связь выходных координат yk (как правило, это измеряемые переменные) с переменными состояния xi;
D – m x r матрица. Она характеризует непосредственное воздействие входов иj на выходы уk.
На рис. 35 изображена структурная схема, построенная по уравнению состояний.
u |
dx/dt |
x |
y |
B |
|
1/S |
C |
A
D
Рис. 35. Структурная схема, построенная по уравнению состояний
При получении уравнений пространства состояний типового динамического звена 2-го порядка (например, двигателя постоянного тока с возбуждением от постоянных магнитов) воспользуемся следующими преобразованиями:
Wд(s) = Кд/(Тм*Тя*s²+Тм*s+1) = (1/Тм*s)*[1/(Тя*s+1)]/
/[1 + (1/Тм*s)*[1/(Тя*s + 1)]].
Структурная схема, реализующая полученную передаточную функцию типового звена 2-го порядка (на примере передаточной функции двигателя постоянного тока) с учетом преобразования апериодического типового динамического звена 1-го порядка и переноса элемента структурной схемы через сумматор против направления передачи сигнала, изображена на рис. 36.
87
u |
dx2/dt |
x2 |
dx1/dt |
x1 = y |
Kд |
1 |
|
|
1 |
Tя*Tм |
S |
|
|
S |
|
1 |
|
|
|
|
Tя |
|
|
|
|
1 |
|
|
|
|
Tя*Tм |
|
|
|
Рис. 36. Структурная схема типового динамического звена 2-го порядка
Матрицы уравнения состояний типового динамического звена 2-го порядка (например, двигателя постоянного тока с
возбуждением от постоянных магнитов) определяются выражениями
|
х1 |
х2 |
|
u |
dx1/dt |
0 |
1 |
dx1/dt |
0 |
А = |
|
|
В = |
|
dx2/dt |
–1/(ТяТм) –1/Тя , |
dx2/dt |
Кд/(ТяТм) , |
|
х1 |
х2 |
|
|
|
С = у 1 |
0 , |
D = 0 |
|
|
ПРЕДВАРИТЕЛЬНОЕ ЗАДАНИЕ
1. Структурная схема исходной системы автоматического регулирования скорости, состоящая из регулятора - П-звено, двигателя постоянного тока –типовое динамическое звено 2-го порядка и датчика – П-звено звено изображена на рис. 37.
88
|
|
Kd |
1 |
1 |
1 |
|
1 |
Kreg |
|
s |
s |
||
T a*T m |
Out1 |
|||||
|
|
|
|
|||
In1 |
|
|
Integrator |
Integrator1 |
||
Gain |
|
|
||||
|
Transfer Fcn3 |
|
|
|
||
|
|
|
1 |
|
|
|
|
|
|
Ta |
|
|
|
|
|
|
Transfer Fcn1 |
|
|
|
|
|
|
1 |
|
|
|
|
|
|
Ta*Tm |
|
|
|
|
|
|
Transfer Fcn2 |
|
|
|
|
|
Ktg |
|
|
|
|
|
|
Gain1 |
|
|
|
Рис. 37. Структурная схема исходной системы автоматического регулирования скорости
2. Двигатель заменяется структурной схемой изображенной на рис. 36 .
3. Структурная схема с доступом к переменным состояния представлена на рис. 38.
4. Преобразовать структурную схему системы автоматического регулирования скорости (рис. 37) в схему для получения переменных уравнения пространства состояний – рис. 38.
1 |
|
|
|
|
|
|
Kreg |
|
|
|
|
Kd |
|
|
1 |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
Tm*Ta.s2+Tm.s+1 |
|
|
|||||||||||
|
|
|
||||||||||||||||
In1 |
|
|
|
|
|
|
|
|
|
|
|
Out1 |
||||||
|
|
|
|
|
|
|
Gain1 |
Transfer Fcn |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
Ktg |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
Gain |
|
|
|
|||
Рис. 38. |
|
|
|
Структурная |
схема системы автоматического |
регулирования скорости с доступом к переменным состояния
89
Таблица 18
Параметры автоматической системы
Вариант |
Регулятор |
Двигатель |
Датчик обратной связи |
|
|
|
|
|
|
|
|
|
|
Кreg |
Kd |
Tm |
Ta |
Кtg |
|
|
1 |
2 |
2 |
0.012 |
0.003 |
0.12 |
|
|
2 |
2 |
2 |
0.014 |
0.003 |
0.12 |
|
|
3 |
2 |
2 |
0.016 |
0.004 |
0.12 |
|
|
4 |
2 |
2 |
0.018 |
0.004 |
0.12 |
|
|
5 |
2 |
2 |
0.020 |
0.005 |
0.12 |
|
|
6 |
3 |
3 |
0.022 |
0.005 |
0.11 |
|
|
7 |
3 |
3 |
0.024 |
0.006 |
0.11 |
|
|
8 |
3 |
3 |
0.026 |
0.006 |
0.11 |
|
|
9 |
3 |
3 |
0.028 |
0.007 |
0.11 |
|
|
10 |
3 |
4 |
0.03 |
0.007 |
0.11 |
|
|
11 |
1 |
4 |
0.032 |
0.008 |
0.13 |
|
|
12 |
1 |
4 |
0.034 |
0.008 |
0.13 |
|
|
13 |
1 |
4 |
0.036 |
0.009 |
0.13 |
|
|
14 |
1 |
4 |
0.038 |
0.009 |
0.13 |
|
|
15 |
1 |
4 |
0.04 |
0.01 |
0.13 |
|
5. Матрицы уравнения состояний для |
структурной схемы |
САР скорости с доступом к переменным состояний определяются
|
х1 |
х2 |
dx1/dt |
0 |
1 |
А = |
|
|
dx2/dt [–1/(Тm*Тa) –Ktg*Kreg*Kd/(Tm*Ta)] |
–1/Tа , |
|
|
u |
|
dx1/dt |
0 |
|
B= |
|
|
dx2/dt |
Kreg*Kd/(Tm*Ta) , |
|
90