Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Моделирование систем управления. Учебно-методическое пособие

.pdf
Скачиваний:
1
Добавлен:
07.09.2026
Размер:
1 Мб
Скачать
☆
ДалеенеобходимовкоманднойстрокеMatlabпереопределить
настроечныекоэффициентырегулятора,приравнявихкединице:
>> kp=1;ki=1;
В элементе PID Controller также следует задать символьное обозначение в строках «Proportional» – kp и «Integral» – ki соответственно. Открытие двойным щелчком блока «Signal Constraint»приведеткпоявлениюокнанастройки(рисунок1.3.8).
Рисунок1.3.8–Окнонастройкирегулятораблока
«SignalConstraint»
Черные ограничительные линии показывают границы, за которые не должна выходить переходная характеристика после расчетапараметровkp,ki.
Если нажать двойным щелчком на пустом поле, то появится окно для задания данных пределов (рисунок 1.3.9). В поле «Start» (начало процесса) в столбце «Time» по умолчанию определены следующиечисловыезначения:
1строка(0)–времяначалапереходногопроцесса;
2строка(0.3000)–времянарастания;
3строка(0.9000)–времярегулирования.
21
Столбец «Amplitude» отображает параметры характеристики,
соответствующиеданномузначению«Time».
Аналогично задаются переменные в поле «End» (конец
процесса):
1строка(0.3000)–времянарастания; 2строка(0.9000)–времярегулирования; 3строка(3)–времямоделирования.
С помощью варьирования числовых значений в столбце «Slope» можно менять наклон ограничительных линий. В колонку «Weight»вносятсяпараметры0,еслинужноотключитьотображение границпроцессаили1,когдаонидолжныприсутствоватьнаграфике.
Рисунок1.3.9–Окнонастройкиграницпереходного
процессаблока«SignalConstraint»
По умолчанию заданы следующие ограничения переходного процесса: время регулирования – 0,9 с, перерегулирование 20 %. Установившеесязначениевблоке«SignalConstraint» по умолчанию приняторавным1.
Для запуска процесса оптимизации следует открыть вкладку «Optimization»,выбрать«TunedParameters»и«Add».Впоявившемся окне«AddParameters»выделитькоэффициенты,которыеотносятсяк регулированиюпотокосцепленияротораkp,ki,затемнажать«OK» (рисунок1.3.10).
В окне «Tuned Parameters» слева можно будет увидеть символьные обозначения введенных коэффициентов. Если их выделять, то справа напротив надписи «Value» отображаются их значения (рисунок 1.3.11). Выбор данных параметров регулятора подтверждаетсянажатиемкнопки«OK».
22
Рисунок1.3.10–Окновыборапараметроврегулятора
блока«SignalConstraint»
Рисунок1.3.11Окноспараметрамирегулятораблока
«SignalConstraint»
Запуск процесса оптимизации осуществляется при помощи
пиктограммы
.
Появится окно с пошаговойоптимизацией(рисунок 1.3.12) и переходными характеристиками, каждая из которых соответствуюет настройкерегуляторанакаждомэтапе(рисунок1.3.13).Врезультате полученыследующиекоэффициентыkp= 168,0098,ki=282,0504. Ихнеобходимоввестивблок«PIDController».
23
Рисунок1.3.12–Окнопошаговойоптимизацииблока
«SignalConstraint»принастройкеконтурарегулирования
потокосцепленияротора
Рисунок1.3.13–Переходныехарактеристики
потокосцепленияроторапринастройкерегулятора
накаждомэтапеблока«SignalConstraint»
24
Окончательный вид переходного процесса потокосцепления
'
k S Sx R m Rx
Li kω − ρω ψ
'
() .
=+⇒
( 1)
1
() .
pp
р пi иi
p
ii
i
i
k Tp
Wp k k
Tp p
+
= = +
iS
TT
ротора АКЗ после настройки коэффициентов регулятора блоком «Signal Constraint» показан на рисунке 1.3.14. Величина перерегулирования составляет 4,29 %, а время регулирования равно0,1с.Результатыпрактическинеотличаютсяотнастройкипо «модульномуоптимуму».
Рисунок1.3.14–Переходнойпроцесспотокосцепления
ротораАКЗпосленастройкикоэффициентоврегулятора
блоком«SignalConstraint»
Контур регулирования тока статора
Для синтеза модели контура регулирования тока статора
iSy воспользуемся вторым уравнением системы (1.2.26) без учета
слагаемых будетопределятьсяаналогично(1.3.1):
.Тогдапередаточнаяфункцияобъекта
u p r T pi p
Sy S Sy
⇒= =
ВыберемПИ-регулятордлякомпенсациипостояннойвремени
T’
спередаточнойфункцией:
S
Приэтомпараметррегулятора:
() (1 ) ()
Wp
i
ip
Sy
()
up
Sy
= (1.3.8)
p
25
( 1)
Tp
'
.
1/
'
S
()
(1.3.6)
r
+
(1.3.7)
Тогда,подставив(1.3.8) в(1.3.7)и,соединивпоследовательно
'' '
() () .
pp
SS S
W pWp
'
'
1
() .
1
p
p
p
i
i
i
S
S
i
k
r
Hp
k
Tr
Tp
p
r
k
= =
+
+
'
'
SS
k
kr
Wpi(p)иWi(p),получим
⋅= =
рi i
Замкнутаясистема(1.3.9)приметвид
'
+
( 1) /
k Tp k r
iS i
Tp Tp Tp
1/
r
( 1)
+
 (1.3.9)
(1.3.10)
Из(1.3.10)следует, чтонастроечныйкоэффициентрегулятора
стоитвыбиратьизусловия k
>>T’S·r.Вэтомслучаепередаточная
pi
функцияприобретаетсвойстваидеальноговоспроизведениявходного сигналасмалойинерционностьюистановитсяравной1.Например, пусть kpi=100T’S·r.
Коэффициент kpi определяется как отношение максимальной
постоянной времени T’S к времени интегрирования Tи [9, 10].
Следовательно,
100
pi
= = =
и
i
''
TT
Tr
Tk=
и
S
100
. Другой коэффициент будет равен
pi
.
T
S
Для расчета параметров контура регулированиятока статора,
необходимовкоманднойстрокеMatlabввестиследующийкод:
>> kpi=100*Ts_*r;kii=100*r;
В приложении Simulink соберем замкнутую систему объекта
управления Transfer fcn с регулятором «PID Controller» (рисунок
1.3.15)согласно(1.3.6)и(1.3.7).
Рисунок1.3.15–Модельконтурастабилизациитока
статора
26
Из полученной временной диаграммы тока статора
()
Mp
()
1
() .
p
ω
(рисунок 1.3.16) нетрудно заметить, что регулятор настроен таким образом, что переходной процесс не имеет перерегулирования, а время регулирования принимает значение примерно t
 = 0,05 с. В
рег
связисэтим,небудемдобиватьсяегоулучшениянастройкойблоками «PIDtuner»и«SignalConstraint».
Рисунок1.3.16–ПереходнойпроцесстокастатораАКЗ
Контур регулирования угловой скорости вала
двигателя
Модель стабилизации угловой скорости вращения АКЗ являетсясистемой подчиненногорегулирования,внутреннийконтур которойподдерживаетназаданномуровнетокстатораiSy.
Передаточная функция по каналу «iSy – электромагнитный момент М» определяется из пятого уравнения системы (1.2.26) при допущении того, что потокосцепление стабилизируется на уровне ψRx=1Вб:
( ) 1, 5 .
Wp k
РазностьмоментовMимоментанагрузкиM
= = ρ (1.3.11)
1
()
ip
Sy
R
связанасугловой
н
скоростьювращенияваладвигателяωmоператором,которыйможно вычислитьнаоснованиишестогоуравнения(1.2.26):
Wp
2
= =− (1.3.12)
m
M M p Jp
( )
27
н
()
Внешний контур образуется последовательным соединением
1, 5 1, 5
11
() .
RR
i
kk
ρρ
( 1)
1
() .
pp
р пи
p
k Tp
Wp k k
Tp p
ωω
ω ωω
ω
+
= = +
1
() ,
1
o
î
k
Wp
Tp Tp
ω
Σ
=
+
оo
T TJ k k
Tk
замкнутойсистемы(1.3.10),W1(p)иW2(p):
Wp
= =
ω
'
Tr
S
k
п
p
Jp p Jp
+
1
0,01 1
+
(1.3.13)
Данный объект содержит интегрирующеезвено, поэтому для достижениязаданноговоспроизведениявходногосигнала настройку регулятора следует осуществлять по критерию «симметричного оптимума». При этом логарифмическая амплитудная частотная характеристика системы приобретает вид, симметричный относительно частоты среза, и максимум запаса устойчивости приходитсянаееокрестность[9,10].
НастроимПИ-регулятордлязаданногоконтура:
(1.3.14)
Объект(1.3.13)относитсяквиду
'
Tr
S
где
= = = = = ρ  Если TΣ < To, то
Σ
0,01; 0, 4; 1,5 .
k
p
i
настроечныепараметрыk
иkиωопределяютследующимобразом:
пω
(1.3.15)
R
TT
4 4 0,04;
= = =
p
ωΣ
T
k
p
ω
о
= =
kT k
2 2 1, 5 0, 01
o
k
= =
и
ω
Схемамоделированияконтурастабилизацииугловойскорости вращенияваладвигателяпредставленанарисунке1.3.17.Изграфика
Tr
100
⋅⋅ ρ
Σ
k
p
ω
p
ω
28
'
S
Tr
S
J
⋅
12
'
J
100
ρ
(1.3.16)
;
R
2
.
R
полученного переходного процесса данной величины определены показатели качества: перерегулирование σ = 43,46 %, время регулированияt
=0,2с(рисунок1.3.18).
рег
Рисунок1.3.17–Модельконтурастабилизацииугловой
скоростивращенияваладвигателя
Рисунок1.3.18–Переходнойпроцессугловойскорости
вращенияваладвигателя
Настройка регулятора с помощью приложения
PID tuner в Matlab
Настройка данного контура регулирования с помощью приложения PID tuner при варьировании параметра «Response time» до значения 0,101 дает переходной процесс со следующими результатами: kpw = 1,77569776876195, kpi = 11,9613114038094. Перерегулирование составляет 20,65 %, что меньше полученного значения в выше приведенном методе, а время регулирования увеличилосьдо0,5с(рисунок1.3.19).
29
Рисунок1.3.19–Переходнойпроцессугловойскорости
вращенияваладвигателяпринастройкеприложением
PIDtuner
Настройка регулятора с помощью приложения
Signal Constraint в Matlab
При настройке блоком «Signal Constraint» в элементе «PID Controller»задаютсясимвольныепеременныеkpw,kiw.Вкомандном окнеимприсваиваетсязначение1.
Рисунок1.3.20–Окнопошаговойоптимизацииблока
«SignalConstraint»принастройкеконтурарегулирования
угловойскорости
30
Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]