Міністерство освіти і науки України
НУ ”Львівська політехніка”
“ІКТА”
Кафедра “Автоматики і телемеханіки”
КУРСОВА РОБОТА
з курсу
“ТЕОРІЯ АВТОМАТИЧНОГО КЕРУВАННЯ”
Тема № 6
Варіант № 3
Львів-2006р.
Тема № 6 Варіант № 3
Завдання.
1.Описати призначення і принцип роботи схеми.
2.Зобразити структурну схему системи.
3.Записати вирази коефіцієнтів передачі окремих ланок, розімкненої і замкнутої системи.
4.Зобразити статичні характеристики окремих ланок системи.
5.Записати вирази диференціальних рівнянь для окремих ланок системи.
6.Знайти рішення диференціальних рівнянь для окремих ланок системи.
7.За результатами рішень диференціальних рівнянь побудувати перехідні характеристики окремих ланок системи.
8.Записати вирази диференціальних рівнянь розімкненої і замкнутої системи.
9.Привести схеми електронного моделювання окремих ланок і замкненої системи.
10.Користуючись одним з алгебраїчних критеріїв стійкості визначити стійкість системи і знайти граничний коефіцієнт підсилення.
11.Записати вирази передаточних функцій для окремих ланок системи.
12.Записати вирази передаточних функцій розімкненої і замкнутої системи.
13.Записати вирази для комплексних коефіцієнтів передачі окремих ланок системи, розімкненої і замкнутої системи.
14.Розрахувати аналітично і побудувати АФХ, ЛАЧХ і ФЧХ окремих ланок системи і розімкненої системи.
15.По АФХ розімкненої системи визначити стійкість системи. Знайти запаси стійкості по амплітуді і фазі.
16.По ЛАЧХ і ФАХ розімкненої системи визначити стійкість системи і знайти запаси стійкості по амплітуді і фазі.
17.Побудувати графік перехідного процесу при одиничній стрибкоподібній дії вхідної величини.
18.По графіку перехідного процесу визначити якісні показники системи.
Слідкуюча система витратоміра літака:
Рівняння ланок:
А) вимірювальна схема
Б) електронний підсилювач
В) двигун
Г) редуктор
Вихідні дані:
Параметри
Розмірність
3
TM
Сек
0.1
TE
Сек
0.05
C
Рад/сек
0.5
Kn
100000
S
В/рад
0,001
і
4
1.Опис роботи системи.
Слідкуюча система витратоміра літака здійснює контроль над витратами літака. Під впливом газового потоку заслонка зміщується. Так як заслонка зв‘язана з повзунком реостата R1, то повзунок також переміщається на певну відстань. В результаті цього (вх стає відмінним від (вих, що приводить до утворення похибки розузгодження (U. Значення (U подається на електронний підсилювач, де підсилюється, і підсилена напруга Up поступає на двигун. Вал двигуна, при надходженні цієї напруги, починає обертатись. Оберти з вала двигуна подаються на редуктор. Редуктор, який зв‘язаний з повзунком реостата R2, переміщує його доти, поки (U=0 і (вх=(вих.
2.Структурна схема системи.
Uвх
Рис.2.1 Структурна схема системи
В структурній схемі використані такі позначення:
ЕП-електронний підсилювач;
ВП1,ВП2-вимірювальні перемворювачі;
ДВ-двигун;
Ред-редуктор;
3.Коефіцієнти передачі ланок, розімкненої і замкненої системи.
Коефіцієнтом передачі називається відношення вихідної величини до вхідної в усталеному режимі, тобто коли всі часові похідні дорівнюють нулю.
Коефіцієнт передачі вимірювальних перетворювачів(ВП1,ВП2):
(В/рад);
(В/рад);
Коефіцієнт передачі електричного підсилювача:
;
Коефіцієнт передачі двигуна:
(Рад/в.сек);
Коефіцієнт передачі редуктора:
;
Для розімкненої системи (розімкнена система – система, в якій відсутній контроль за регульованою (вихідною) величиною, в даному випадку для розімкненої системи ) коефіцієнт передачі є добутком коефіцієнтів передачі ланок ВП1, ЕП, ДВ, Ред., оскільки для розімкненої системи вони включені послідовно. Отже коефіцієнт передачі розімкненої системи:
=12.5;
Для замкненої системи (замкнена система – система, в якій відбувається постійний контроль за регульованою (вихідною) величиною) коефіцієнт передачі є добутком коефіцієнта передачі ланки ВП1 та коефіцієнта передачі ланок ЕП, ДВ, Ред з врахуванням охоплення їх від’ємним оберненим зв’язком:
.
Статичні характеристики окремих ланок
Ланки системи:
вимірювальна схема
b) електронний підсилювач
двигун
=>
редуктор
Вирази диференціальних рівнянь для окремих ланок системи
У даній системі перший і другий вимірювальні перетворювачі, а також електронний підсилювач, двигун і редуктор є без інерційними елементами, тому вони описуються алгебраїчними рівняннями. Двигун же є інерційним елементом, тому він описується диференційним рівнянням:
де:
;
C=0.5.
підставивши ці значення у диференційне рівняння отримуємо остаточне рівняння:
Розв’язок диференціальних рівнянь ланок системи.
Розв’язок неоднорідного диференціального рівняння є сумою повного розв’язку однорідного диференціального рівняння (вільна складова) і часткового розв’язку неоднорідного диференціального рівняння (усталений режим). При розв’язку однорідного диференціального рівняння шукана функція замінюється виразом , де - змінна інтегрування, p – стала. Після відповідних підстановок рівняння зводиться до алгебраїчного і розв’язується відносно p. Маючи p можна знайти розв’язок розв’язку однорідного рівняння. Знаходження розв’язку для усталеного режиму (статика) досягається прирівнюванням всіх похідних до нуля в неоднорідному диференціальному рівнянні.
Двигун. Однорідне диференціальне рівняння:
,
нехай, тоді , скоротивши на отримаємо:
дискримінант , , тоді , .
Отже , де - сталі.
Частковий розв’язок неоднорідного рівняння .
Отже, .
При подачі на вхід одиничної стрибкоподібної функції розв’язок записується як. Перша похідна .
Тоді при нульових незалежних початкових умовах:
Розв’язуючи систему рівнянь отримаємо: .Таким чином при одиничній стрибкоподібній функції на вході ланка описується виразом:
.
Перехідні характеристики ланок системи
Рис.7.1 Перехідна характеристика двигуна, як розв'язок диференціального рівняння двигуна: =>
t
W(t)
0
0
0,005
0,024385288
0,01
0,047581291
0,015
0,069646012
0,02
0,090634623
0,025
0,110599608
0,03
0,12959089
0,035
0,147655955
0,04
0,164839977
0,045
0,181185924
0,05
0,19673467
0,055
0,211525095
0,06
0,225594182
0,065
0,238977112
0,07
0,251707348
0,075
0,263816724
0,08
0,275335518
0,085
0,286292534
0,09
0,29671517
0,095
0,306629488
0,1
0,316060279
0,105
0,325031125
0,11
0,333564458
0,115
0,341681615
0,12
0,349402894
0,125
0,356747602
Інші ланки системи описуються алгебраїчними рівняннями, тобто є безінерційними, тому їх перехідні характеристики є стрибкоподібною функцію (повторює форму вхідного сигналу) з амплітудою рівною коефіцієнту передачі відповідної ланки.
Вирази диференціальних рівнянь розімкненої і замкненої системи
Диференціальне рівняння замкненої системи.
Підставимо рівняння редуктора у рівняння двигуна отримуємо рівняння . Отримане рівняння підставимо і рівняння електронного підсилювача , отримуємо рівняння . Це отримане рівняння підставляємо у рівняння вимірювальної схеми , отримуємо рівняння . Розкриваємо усі дужки у отриманому рівнянні . Замінюємо на і отримуємо таке рівняння . Замінюємо коефіцієнти на відповідні a і b, отримуємо рівняння .
Числові значення коефіцієнтів a і b будуть такими:
Диференціальне рівняння розімкненої системи.
Для розімкненої системи рівняння вимірювальної схеми набуде вигляду , тому диференційне рівняння майже таким як і для замкненої системи тільки без . Отже диференціальне рівняння набуде такого вигляду . Де коефіцієнти будуть аналогічними до коефіцієнтів замкненої системи.
9.Схеми електронного моделювання ланок і замкненої системи.
Електронне моделювання ланок і систем здійснюється при певній комбінації включення типових ланок САУ. Безінерційні елементи системи (вимірювальний перетворювач, електронний підсилювач, редуктор) моделюються за допомогою безінерційних ланок з відповідними коефіцієнтами передачі.
Рис.9.1 Схема електронного моделювання ВП1 вимірювального перетворювача;
Рис.9.2 Схема електронного моделювання ЕП електронного підсилювача;
Рис.9.3 Схема електронного моделювання Ред. редуктора;
Рис.9.4 Схема електронного моделювання ВП2 вимірювального перетворювача;
Двигун описується рівнянням третього порядку. Його моделювання здійснюється за допомогою коливної і інтегруючої ланок:
Рис.9.5 Схема електронного моделювання ДВ двигуна;
Схему порівняння можна промоделювати за допомогою підсумовуючої безінерційної ланки при умові, що сигнали на нього подаються у протилежних фазах. Оскільки у всіх ланках використовується інверсне включення операційних підсилювачів, то кожна моделююча ланка змінює знак (фазу) сигналу на протилежний.
Схема електронного моделювання замкненої системи буде мати такий вигляд:
Рис.9.6 Схема електронного моделювання замкненої системи;
Визначення стійкості системи і визначення граничного коефіцієнта підсилення
Визначення стійкості системи за допомогою критерію стійкості Гурвіца.
По критерію стійкості системи Гурвіца система автоматичного управління буде стійкою, якщо визначник Гурвіца і усі діагональні мінори додатні.
Приведемо диференційне рівняння замкненої системи до операторного вигляду і прирівняємо до нуля. Кінцеве рівняння набуде вигляду . Де:
Складаємо визначник Гурвіца і його діагональні мінори:
Визначення коефіцієнта підсилення.
Граничний коефіцієнт підсилення визначається з умови , де . Звідси визначаємо . Таким чином коефіцієнт підсилення при якому система залишається на межі стійкості рівний 800.
Вирази передаточних функцій для окремих ланок системи
перший і другий вимірювальний перетворювач
електронний підсилювач
двигун
редуктор
Вирази передаточних функцій для розімкненої і замкненої системи
Вираз передаточної функції для розімкнена система.
Вираз передаточної функції для замкненої система.
.
Вирази передаточних функцій для комплексних коефіцієнтів передачі окремих ланок системи, розімкненої і замкненої системи.
Вирази передаточних функцій для комплексних коефіцієнтів передачі окремих ланок системи.
перший і другий вимірювальний перетворювач
електронний підсилювач
двигун
редуктор
Вирази передаточної функцій для розімкненої системи.
Вираз передаточної функції для замкненої системи.
.
Аналітичний розрахунок і побудова АФХ, ЛАЧХ і ФЧХ окремих ланок системи і розімкненої системи.
Без інерційні ланки.
Передаточні функції для без інерційних ланок, а саме: перший і другий вимірювальний перетворювач, електронний підсилювач і редуктор мають вигляд дійсного числа і не залежать від частоти.
Амплітудно фазова характеристика для без інерційних ланок являє собою одну точку з координатами , де K – це коефіцієнт передачі відповідної ланки.
Амплітудна частотна характеристика для без інерційних ланок – це пряма, яка паралельна осі частот з координатою рівня K (для лінійної АФХ) або , тому що , де K – це коефіцієнт передачі відповідної ланки.
АФХ для без інерційних ланок:
a) перший і другий вимірювальний перетворювач
U(w) = S = 0.001, V(w)=0 – вимірювальні перетворювачі мають тільки дійсну складову передаточної функції. A(w) = S = 0,001
L(w) =20lg S = 20lg 0,001 = -60
b) електронний підсилювач
U(w) = Кп = 100000, V(w)=0 – електронний підсилювач має тільки дійсну складову передаточної функції.
A(w) = Кп = 100000
L(w) =20lg Кп = 20lg 100000 = 100
c) редуктор
U(w) = 1/і =0.25; V(w)=0 – аналогічно, редуктор, як безінерційна ланка в даній системі має тільки дійсну складову передаточної функції.
А(w)=0.25;
L(w)=20Lg(A(w))=-12.04;
ЛАЧХ для безінерційних ланок:
перший і другий вимірювальний перетворювач
U(w) = S = 0.001, V(w)=0 – вимірювальні перетворювачі мають тільки дійсну складову передаточної функції.
A(w) = S = 0.001
L(w) =20lg S = 20lg 0.001 = -60
електронний підсилювач
U(w) = Кп = 100000, V(w)=0 – електронний підсилювач має тільки дійсну складову передаточної функції.
A(w) = Кп = 100000
L(w) =20lg Кп = 20lg 100000 = 100
редуктор
U(w) = 1/і = 0.25; V(w)=0 – аналогічно, редуктор, як безінерційна ланка в даній системі має тільки дійсну складову передаточної функції.
А(w)=0.25;
L(w)=20Lg(A(w))=-12.04;
ФЧХ для безінерційних ланок: Фазочастотні характеристики , це пряма яка проходить по осі абсцис частот:
Інерційні ланки
У нашій системі тільки одна інерційна ланка – це двигун, і тому тільки для нього я будую АФХ, ЛАЧХ і ФЧХ.
З комплексного коефіцієнта передачі двигуна я знаходжу і :
=U(w)+jV(w)
Тоді
АФХ
w
V(w)
U(w)
0,1
-5,00200055
-0,050045028
0,2
-2,501499898
-0,050080078
0,3
-1,667665011
-0,050105108
0,4
-1,250496189
-0,050120088
0,5
-0,999993734
-0,050124999
0,6
-0,832824617
-0,050119837
0,7
-0,713274854
-0,050104606
0,8
-0,623487604
-0,050079326
0,9
-0,553542535
-0,050044026
1
-0,497487563
-0,04999875
1,1
-0,451535107
-0,049943551
1,2
-0,413160208
-0,049878496
1,3
-0,380614906
-0,049803662
1,4
-0,352650733
-0,049719137
1,5
-0,328352222
-0,049625021
ЛАЧХ
Якщо для двигуна з рівняння передаточної функції, то
lg(w)
L(w)
0
-6,020708486
0,301029996
-12,1276273
0,602059991
-18,55553074
0,903089987
-26,01809766
1,204119983
-35,42595601
1,505149978
-46,51599521
1,806179974
-58,31655197
2,10720997
-70,31243268
2,408239965
-82,35071804
2,709269961
-94,39542823
3,010299957
-106,4396229
ЛФЧХ
Якщо, то
lg(w)
Fi(w)
0
-0,003999979
0,301029996
-0,007999829
0,602059991
-0,015998635
0,903089987
-0,031989084
1,204119983
-0,063912833
1,505149978
-0,127307742
1,806179974
-0,250617701
2,10720997
-0,473201458
2,408239965
-0,797255315
2,709269961
-1,116567592
3,010299957
-1,331339927
3,311329952
-1,44932698
3,612359948
-1,509836793
3,913389944
-1,540288217
Розімкнена система
З комплексного рівняння розімкненої я знаходжу і :
Тоді
АФХ
w
V(w)
U(w)
0,1
-125,0500137
-1,2511257
0,2
-62,53749746
-1,252001951
0,3
-41,69162527
-1,2526277
0,4
-31,26240472
-1,253002193
0,5
-24,99984336
-1,25312498
0,6
-20,82061542
-1,252995913
0,7
-17,83187135
-1,252615147
0,8
-15,58719011
-1,251983141
0,9
-13,83856337
-1,251100656
1
-12,43718907
-1,249968751
1,1
-11,28837767
-1,248588783
1,2
-10,32900521
-1,2469624
1,3
-9,515372638
-1,245091538
1,4
-8,816268329
-1,242978415
1,5
-8,208805546
-1,240625523
1,6
-7,675820867
-1,238035624
1,7
-7,204190803
-1,235211736
1,8
-6,783709547
-1,232157133
1,9
-6,406321107
-1,228875326
2
-6,065581806
-1,225370062
2,1
-5,756276339
-1,221645307
2,2
-5,474138553
-1,217705239
2,3
-5,215645053
-1,213554236
2,4
-4,977860412
-1,209196861
З комплексного рівняння розімкненої я знаходжу і :
Тоді
ЛАЧХ
lg(w)
L(w)
0
21,93809169
0,301029996
15,83117287
0,602059991
9,403269432
0,903089987
1,940702517
1,204119983
-7,467155834
1,505149978
-18,55719504
1,806179974
-30,3577518
2,10720997
-42,3536325
2,408239965
-54,39191786
2,709269961
-66,43662806
3,010299957
-78,4808227
Тоді ЛФЧХ
lg(w)
Fi(w)
0
-0,003999979
0,301029996
-0,007999829
0,602059991
-0,015998635
0,903089987
-0,031989084
1,204119983
-0,063912833
1,505149978
-0,127307742
1,806179974
-0,250617701
2,10720997
-0,473201458
2,408239965
-0,797255315
2,709269961
-1,116567592
3,010299957
-1,331339927
3,311329952
-1,44932698
3,612359948
-1,509836793
3,913389944
-1,540288217
Визначення стійкості системи по АФХ розімкненої системи, запаси стійкості по амплітуді і по фазі
Визначення стійкості системи за виглядом частотних характеристик здійснюється за допомогою частотних критеріїв стійкості, зокрема за критерієм Найквіста: САУ буде стійка у замкненому стані, якщо розімкнена системи на комплексній площині не охоплює точку з координатами (-1; j0).
АФХ розімкненої системи (фрагмент):
З комплексного рівняння розімкненої я знаходжу і :
Тоді
АФХ
w
V(w)
U(w)
7,8
-1,005315615
-0,81596897
7,9
-0,982619474
-0,808193008
8
-0,960553279
-0,800461066
8,1
-0,939095629
-0,79277484
8,2
-0,918226047
-0,785135932
8,3
-0,89792493
-0,777545844
8,4
-0,878173495
-0,770005988
8,5
-0,858953741
-0,762517681
8,6
-0,840248399
-0,755082156
8,7
-0,822040902
-0,747700559
8,8
-0,804315341
-0,740373954
8,9
-0,787056434
-0,733103325
9
-0,770249496
-0,725889578
9,1
-0,753880406
-0,718733545
9,2
-0,737935579
-0,711635989
9,3
-0,722401946
-0,704597598
9,4
-0,707266921
-0,697618998
9,5
-0,692518384
-0,69070075
9,6
-0,678144658
-0,683843353
9,7
-0,664134489
-0,677047245
9,8
-0,650477024
-0,670312811
9,9
-0,637161797
-0,663640378
10
-0,624178712
-0,657030223
7,8
-1,005315615
-0,81596897
Як видно з малюнку АФХ розімкненої системи не охоплює точку з координатами (-1; j0), отже система в замкненому стані буде стійкою.
A=0.95.
Запас стійкості по фазі дорівнює φ ≈63 градусів – це кут між радіусом одиничного кола, проведеного в точку перетину з графіком і віссю абсцис.
Визначення по ЛАЧХ і ФЧХ розімкненої системи стійкість системи і значення запасів стійкості по амплітуді і по фазі
Для того щоб САУ була стійкою, необхідно і достатньо, щоб при досягненні ЛАЧХ значення , .
ЛАЧХ і ФЧХ (фрагменти), визначення стійкості системи:
lg(w)
L(w)
Fi(w)
0
21,93809169
-0,003999979
0,301029996
15,83117287
-0,007999829
0,602059991
9,403269432
-0,015998635
0,903089987
1,940702517
-0,031989084
1,204119983
-7,467155834
-0,063912833
1,505149978
-18,55719504
-0,127307742
1,806179974
-30,3577518
-0,250617701
2,10720997
-42,3536325
-0,473201458
2,408239965
-54,39191786
-0,797255315
2,709269961
-66,43662806
-1,116567592
3,010299957
-78,4808227
-1,331339927
Як бачимо з рисунка, ЛФЧХ розімкненої системи не досягає значення –п, отже, це випадок, при якому неможливо визначити стійкість системи за ЛАЧХ та ЛФЧХ РС.
Графік перехідного процесу при одиничній стрибкоподібній дії вхідної величини
Для отримання функції перехідного процесу при подачі на вхід одиничної стрибкоподібної функції скористаємося математичним пакетом MatLab 6.1. Введемо необхідні дані передаточної функції замкненої системи в оператор STEP за принципом: чисельник – 20, а знаменник – це коефіцієнти 0.001 0.1 1 20.
.
step([12.5],[0,005 0,1 1 12.5]),grid
Якісні показники системи
За графіком:
- необхідне значення регульованої величини,
- максимальне значення регульованої величини,
- час перехідного процесу при ,
- час встановлення максимального значення,
Висновок
Система автоматичного регулювання при дослідженнях показала себе досить стійкою, з малим часом перехідного процесу () при відносно великих вимогах до неї (). САК володіє нормальними запасами стійкості, що забезпечить стабільну роботу під час випробувань.
Використана література
П.И. Чинаєв, Н.М. Чумаков. Теория автоматического управления. Києвскоє высшее инженерно-авиационноэ военное училище ВВС. Киэв-1969.
Теория автоматического управления: Учеб. для вузов по спец. «Автоматика и телемеханика». В 2-х ч. Ч. І. Теория линейных систем автоматического управления/ Н.А. Бабаков, А.А. Воронов, и др.; Под ред. А.А. Воронова. - 2-е изд., перераб. и доп. – М.:Высш. шк.,1986. – 367 с.,ил.