Міністерство освіти і науки України
Національний університет “Львівська політехніка”
Кафедра АТХП
Лабораторна робота №6
ДОСЛІДЖЕННЯ ЧАСТОТНИХ ХАРАКТЕРИСТИК
ТИПОВИХ ЛАНОК НА ЕЛЕКТРИЧНИХ СХЕМАХ
Виконала: ст. гр. АВ-32
Прийняв: Матіко Ф.Д.
Львів – 2006 р.
Мета роботи: дослідити стійкість систем автоматичного регулювання (САР) за допомогою критеріїв Гурвіца та Найквіста, проаналізувати фактори, що можуть викликати втрату стійкості САР, дослідити вплив зміни параметрів елементів САР на її стійкість.
Необхідна теоретична підготовка: поняття стійкості динамічних систем, критерії стійкості систем.
ПОСЛІДОВНІСТЬ ВИКОНАННЯ РОБОТИ
Змістом виконання роботи є дослідження стійкості заданої системи автоматичного регулювання. За допомогою критеріїв Гурвіца (частина 1) і Найквіста (частина 2) необхідно визначити критичні значення коефіцієнтів настроювання автоматичного регулятора, а також значення параметрів настроювання, що забезпечують заданий запас стійкості за амплітудою або за фазою. Далі в середовищі Matlab-Simulink складають структурну схему системи і в результаті моделювання отримують перехідні функції замкнутої системи. Вихідними даними для виконання лабораторної роботи є функція передачі об’єкту регулювання з конкретними числовими значеннями параметрів, та функція передачі автоматичного регулятора з невідомими параметрами настроювання.
В частині 1 лабораторної роботи студентові необхідно:
1. Скласти характеристичне рівняння заданої САР, розрахувати числові значення його коефіцієнтів.
2. Скласти визначник Гурвіца та розрахувати критичні значення параметрів настроювання регулятор.
3. Скласти структурну схему замкнутої системи в Simulink, встановити задані та розраховані значення коефіцієнтів моделей.
4. Здійснити моделювання САР при одиничній стрибкоподібній зміні вхідної величини та отримати перехідну функцію САР.
5.Дослідити вплив зміни параметрів настроювання регулятора на перехідний процес в замкнутій системі. Порівняти перехідні функції, отримані принаймні при трьох значеннях параметрів настроювання регулятора (рівних, менших та більших за критичне значення), визначаючи кожний раз характер та вигляд перехідної функції, а також коефіцієнт заникання перехідного процесу, де А1, А3 – відповідно перша та третя амплітуди коливного процесу.
В частині 2 лабораторної роботи студентові необхідно:
1. На основі критерію Найквіста, розрахувати критичні значення параметрів настроювання заданого регулятора.
2. Розрахувати значення параметрів настроювання регулятора на заданий викладачем запас стійкості по фазі або по амплітуді.
3. Побудувати в середовищі Matlab амплітудно-фазові характеристики розімкнутої САР при розрахованих значеннях параметрів настроювання регулятора та оцінити графічно запас стійкості системи за фазою та амплітудою.
4. Скласти структурну схему замкнутої САР в Simulink, встановити задані значення параметрів об’єкта регулювання та розраховані значення параметрів регулятора.
5. Здійснити моделювання САР при одиничній стрибкоподібній зміні вхідної величини та отримати перехідні функції САР для критичного значення параметру регулятора та для параметра регулятора, що забезпечує заданий запас стійкості.
6. Дослідити вплив зміни параметрів об’єкта регулювання (наприклад, збільшення запізнення, збільшення-зменшення сталих часу, збільшення-зменшення коефіцієнта передачі) на перехідний процес в САР. Для цього порівняти перехідні функції при різних значеннях параметрів системи, кожний раз визначаючи коефіцієнт заникання та запас стійкості САР за амплітудою і фазою.
ВИКОНАННЯ ЗАВДАННЯ
а) Дослідження стійкості САР за допомогою критерію Гурвіца
Завдання: Дослідити стійкість системи автоматичного регулювання з І-регулятором, структурна схема якої показана на рис.1. Функція передачі об’єкту регулювання
Функція передачі автоматичного І- регулятора
,
де - коефіцієнт передачі регулятора.
Для дослідження стійкості системи автоматичного регулювання за критерієм Гурвіца необхідно отримати характе-ристичне рівняння замкнутої САР. Функція передачі заданої САР має вигляд
Характеристичне рівняння одержимо, прирівнявши знаменник отриманої функції передачі (характеристичний поліном системи) до нуля
.
Загальний вигляд характеристичного рівняння третього порядку є таким
.
Тобто для досліджуваної системи ; ; ; . Як бачимо, всі коефіцієнти є одного знаку (додатними) при будь-яких значеннях коефіцієнта передачі І-регулятора, отже необхідна умова стійкості виконується.
Для визначення умов стійкості сформуємо головний визначник Гурвіца
та визначники діагональних мінорів нижчих порядків
,
.
Досліджувана система буде стійкою, якщо >0 і >0 і >0. Визначник є додатний, знак визначається тільки знаком , тому для того, щоб САР була стійкою достатньо щоб >0. Якщо , то САР буде на границі стійкості, при - нестійка.
Знайдемо критичне значення коефіцієнта передачі І-регулятора, при якому система буде знаходитися на границі стійкості з рівняння
.
Для забезпечення стійкості САР, значення коефіцієнта передачі І-регулятора повинно бути меншим від критичного . Якщо коефіцієнт передачі регулятора виявиться більшим за критичне значення , то система буде нестійкою.
Змоделюємо досліджувану САР в середовищі Simulink:
Рис 2. Структурна модель САР з І-регулятором.
Перехідна функція, отримана при критичному значенні =0,134, показана на рис.3.
Рис. 3. Перехідна функція САР, отримана при =0,134.
Перехідна функція є коливною із незаникаючими коливаннями. Усталені коливання здійснюються відносно значення вихідної величини 0. Амплітуда усталених коливань А= 0.9 , період коливань 51 с.
Змоделюємо перехідну функцію САР із значенням коефіцієнта передачі I-регулятора, більшим за критичне значення =0,134 , наприклад при =0,3 (рис. 4).
Отриманий перехідний процес в САР є коливним, розбіжним. Коливання здійснюються відносно значення регульованої величини, що дорівнює 0. Період коливань 45 с.
Рис. 4. Перехідна функція САР, отримана при =0,3.
Далі змоделюємо перехідну функцію САР із значенням коефіцієнта передачі І-регулятора, меншим за критичне значення, наприклад при =0,05 (рис. 5).
Усталене значення вихідної величини після завершення перехідного процесу (статична похибка регулювання) дорівнює 0. Період коливань становить 79 с, коефіцієнт заникання коливань 0,81.
Рис. 5. Перехідна функція САР, отримана при =0,05.
Висновок: Таким чином, розраховане критичне значення параметра настроювання І-регулятора =0,134 дійсно виводить систему на границю стійкості. Якщо =0,3, то САР є нестійкою, при =0,05<0,134 САР стає стійкою. Отримані перехідні функції САР підтверджують правильність розрахунку критичного значення коефіцієнта передачі І-регулятора за критерієм Гурвіца. Перехідні процеси в САР є коливними, період коливань зменшується із збільшенням коефіцієнта передачі І-регулятора. Дослідження заданої САР показує, що із зменшенням коефіцієнта передачі І-регулятора статична похибка регулювання залишається сталою.
б) Дослідження стійкості САР за допомогою критерію Найквіста
Завдання: Об’єктом дослідження в другій частині лабораторної роботи є та сама САР, що і в першій частині.
Задане значення запасу стійкості САР: запас стійкості по фазі γ=π/6.
Задані зміни параметрів: (збільшення коефіцієнта передачі на 20%).
За допомогою критерію Найквіста розрахувати критичне значення параметра настроювання І-регулятора, а також таке його значення, що забезпечує заданий запас стійкості САР. Дослідити перехідну функцію САР із заданим запасом стійкості, а також дослідити вплив збільшення коефіцієнта передачі .
Для того щоб спростити знаходження виразів частотних характеристик, функцію передачі об’єкта регулювання доцільно представити як добуток двох аперіодичних ланок, тобто представити послідовним з’єднанням цих ланок, як показано на рис.6.
Рис. 6. Структурна схема досліджуваної САР.
Згідно з критерієм Найквіста амплітудно-фазова-характеристика (АФХ) САР, що знаходиться на границі стійкості проходить через точку (-1,і0), тобто виконується умова (25)
або
,
де - амплітудно-фазова характеристика розімкнутої системи; амплітудно-частотні характеристики об’єкта регулювання і автоматичного регулятора відповідно; фазочастотні характеристики відповідно об’єкта регулювання і автоматичного регулятора.
Визначимо необхідні амплітудно- і фазочастотні характеристики. Відомо, що АЧХ послідовного з’єднання елементів дорівнює добутку амплітудно-частотних характеристик елементів, що складають з’єднання, а ФЧХ – сумі їх фазочастотних характеристик. Відповідно АЧХ і ФЧХ об’єкту регулювання матимуть вигляд
,
АЧХ і ФЧХ І-регулятора -
.
Знайдені вирази характеристик підставимо в систему
З рівнянь системи визначимо частоту та критичне значення коефіцієнта передачі І-регулятора k за допомогою числового методом із застосуванням функції fsolve в Matlab.
x0=[0.1;0.1];y=fsolve('fff',x0,3e-6);
де
function y=fff(x);
k=x(1); w=x(2);
y1=atan(x(2)*7)+atan(8*x(2))-pi/2;
y2=2*x(1)/x(2)/sqrt((49*x(2)^2+1)*(64*x(2)^2+1))-1;
y=[y1;y2]
x0=[0.1;0.1]
В результаті виконання програми отримаємо розв’язок системи: ,
Для побудови АФХ виконаємо наступну програму в Matlab:
%побудова АФХ розімкнутої САР на границі стійкості
w=[0:0.01:5];
k=0.05143435954093;
A1=2./sqrt(49.*w.^2+1);
A2=1./sqrt(64.*w.^2+1);
A3=k./w;
fi1=-atan(7*w);
fi2=-atan(8*w);
fi3=-pi/2;
A=A1.*A2.*A3;
fi=fi1+fi2+fi3;
polar(fi,A),grid;
На рис. 7 показана амплітудно-фазова характеристики розімкнутої системи, яка підтверджує, що розрахунок критичного значення I-регулятора виконаний правильно, оскільки АФХ проходить через точку з радіусом 1 і кутом 1800, що відповідає точці (-1, і0) в декартових координатах.
Рис. 7. АФХ розімкнутої САР (замкнута САР на границі стійкості).
Для розрахунку значення коефіцієнта та частоти зрізу , при якому забезпечується заданий запас стійкості системи по фазі γ=π/6 складають систему
.
Для розрахунків скористаємося наведеною вище програмою, яка відрізнятиметься від неї лише виглядом допоміжної функції y2, а саме
function y=fff1(x);
k=x(1); w=x(2);
ga=pi/6;
y1=atan(x(2)*7)+atan(8*x(2))-pi/2+ga;
y2=2*x(1)/x(2)/sqrt((49*x(2)^2+1)*(64*x(2)^2+1))-1;
y=[y1;y2]
x0=[0.1;0.1]
В результаті виконання програми отримуємо значення частоти зрізу і коефіцієнта передачі І-регулятора k= =0.05143435954093.
Рис. 8. АФХ розімкнутої САР із заданим запасом стійкості по фазі (крива 1) та АФХ системи із коефіцієнтом передачі (крива 2).
На рис. 8 (крива 1) показана АФХ розімкнутої САР при розрахованому значенні параметра регулятора, з якої видно, що запас стійкості САР становить π/6.
Побудуємо також перехідну функцію розрахованої САР.
Рис. 9. Перехідна функція САР із заданим запасом стійкості по фазі.
Статична похибка регулювання дорівнює 0, А1 = 1,19, А3 = 0,23 і коефіцієнт заникання коливань .
Далі оцінимо вплив збільшення коефіцієнта передачі (0.06172123144912) на перехідну функцію. Як видно із рис. 10, для перехідної функції системи регулювання із статична похибка регулювання дорівнює 0, перша амплітуда , третя амплітуда , коефіцієнт заникання коливань . Для порівняння на рис. 10 пунктирною лінією показана також перехідна функція з рис. 9.
Рис. 10. Перехідні функції САР з коефіцієнтом передачі k=0.05143435954093 (штрихова лінія) та коефіцієнтом передачі k=0.06172123144912 (суцільна лінія).
Порівнюючи ці два графіка, бачимо, що при збільшенні коефіцієнта передачі статична похибка регулювання не змінилася, а запас стійкості САР зменшився, оскільки коефіцієнт заникання зменшився з 0,80 до 0,73. Зменшення запасу стійкості САР із збільшеним коефіцієнтом передачі на 20% проілюстровано також за допомогою АФХ розімкнутої системи на рис.8 (крива 2), з якої видно, що запас стійкості по фазі зменшився до 250.
Висновок: Знайдене критичне значення параметру І-регулятора за допомогою критерію Найквіста збігається з розрахованим за критерієм Гурвіца, що підтверджує правильність розрахунків. Застосування критерію Найквіста дозволяє розраховувати САР на заданий запас стійкості. Дослідження заданої САР показують, що із збільшенням коефіцієнта передачі запас стійкості САР зменшується і по фазі, і по амплітуді.