Міністерство освіти і науки України
НУ ”Львівська політехніка”
“ІКТА”
Кафедра “Автоматики і телемеханіки”
КУРСОВА РОБОТА
з курсу
“ТЕОРІЯ АВТОМАТИЧНОГО КЕРУВАННЯ”
Тема № 6
Варіант № 9
Львів-2006р.
Тема № 6. Варіант № 9.
Завдання.
1.Описати призначення і принцип роботи схеми.
2.Зобразити структурну схему системи.
3.Записати вирази коефіцієнтів передачі окремих ланок, розімкненої і замкнутої системи.
4.Зобразити статичні характеристики окремих ланок системи.
5.Записати вирази диференціальних рівнянь для окремих ланок системи.
6.Знайти рішення диференціальних рівнянь для окремих ланок системи.
7.За результатами рішень диференціальних рівнянь побудувати перехідні характеристики окремих ланок системи.
8.Записати вирази диференціальних рівнянь розімкненої і замкнутої системи.
9.Привести схеми електронного моделювання окремих ланок і замкненої системи.
10.Користуючись одним з алгебраїчних критеріїв стійкості визначити стійкість системи і знайти граничний коефіцієнт підсилення.
11.Записати вирази передаточних функцій для окремих ланок системи.
12.Записати вирази передаточних функцій розімкненої і замкнутої системи.
13.Записати вирази для комплексних коефіцієнтів передачі окремих ланок системи, розімкненої і замкнутої системи.
14.Розрахувати аналітично і побудувати АФХ, ЛАЧХ і ФЧХ окремих ланок системи і розімкненої системи.
15.По АФХ розімкненої системи визначити стійкість системи. Знайти запаси стійкості по амплітуді і фазі.
16.По ЛАЧХ і ФАХ розімкненої системи визначити стійкість системи і знайти запаси стійкості по амплітуді і фазі.
17.Побудувати графік перехідного процесу при одиничній стрибкоподібній дії вхідної величини.
18.По графіку перехідного процесу визначити якісні показники системи.
Слідкуюча система витратоміра літака:
Рівняння ланок:
А) вимірювальна схема
Б) електронний підсилювач
В) двигун
Г) редуктор
Вихідні дані:
Параметри
Розмірність
9
TM
Сек
0.3
TE
Сек
0.02
C
Рад/сек
4
Kn
10
S
В/рад
10
і
12
1.Опис роботи системи.
Слідкуюча система витратоміра літака здійснює контроль над витратами літака. Під впливом газового потоку заслонка зміщується. Так як заслонка зв‘язана з повзунком реостата R1, то повзунок також переміщається на певну відстань. В результаті цього (вх стає відмінним від (вих, що приводить до утворення похибки розузгодження (U. Значення (U подається на електронний підсилювач, де підсилюється, і підсилена напруга Up поступає на двигун. Вал двигуна, при надходженні цієї напруги, починає обертатись. Оберти з вала двигуна подаються на редуктор. Редуктор, який зв‘язаний з повзунком реостата R2, переміщує його доти, поки (U=0 і (вх=(вих.
2.Структурна схема системи.
(вх Uвх (U Up ( (вих
Uвих Рис.2.1 Структурна схема системи
В структурній схемі використані такі позначення:
ЕП-електронний підсилювач;
ВП1,ВП2-вимірювальні перемворювачі;
ДВ-двигун;
Ред-редуктор;
3.Коефіцієнти передачі ланок, розімкненої і замкненої системи.
Коефіцієнтом передачі називається відношення вихідної величини до вхідної в усталеному режимі, тобто коли всі часові похідні дорівнюють нулю.
Коефіцієнт передачі вимірювальних перетворювачів(ВП1,ВП2):
(В/рад);
(В/рад);
Коефіцієнт передачі електричного підсилювача:
;
Коефіцієнт передачі двигуна:
(Рад/в.сек);
Коефіцієнт передачі редуктора:
;
Для розімкненої системи (розімкнена система – система, в якій відсутній контроль за регульованою (вихідною) величиною, в даному випадку для розімкненої системи ) коефіцієнт передачі є добутком коефіцієнтів передачі ланок ВП1, ЕП, ДВ, Ред., оскільки для розімкненої системи вони включені послідовно. Отже коефіцієнт передачі розімкненої системи:
=33.33;
Для замкненої системи (замкнена система – система, в якій відбувається постійний контроль за регульованою (вихідною) величиною) коефіцієнт передачі є добутком коефіцієнта передачі ланки ВП1 та коефіцієнта передачі ланок ЕП, ДВ, Ред з врахуванням охоплення їх від’ємним оберненим зв’язком:
.
Статичні характеристики окремих ланок
Ланки системи:
вимірювальна схема
b) електронний підсилювач
двигун
редуктор
Вирази диференціальних рівнянь для окремих ланок системи
У даній системі перший і другий вимірювальні перетворювачі, а також електронний підсилювач, двигун і редуктор є без інерційними елементами, тому вони описуються алгебраїчними рівняннями. Двигун же є інерційним елементом, тому він описується диференційним рівнянням:
де:
;
C=4.
підставивши ці значення у диференційне рівняння отримуємо остаточне рівняння:
Розв’язок диференціальних рівнянь ланок системи.
Розв’язок неоднорідного диференціального рівняння є сумою повного розв’язку однорідного диференціального рівняння (вільна складова) і часткового розв’язку неоднорідного диференціального рівняння (усталений режим). При розв’язку однорідного диференціального рівняння шукана функція замінюється виразом , де - змінна інтегрування, p – стала. Після відповідних підстановок рівняння зводиться до алгебраїчного і розв’язується відносно p. Маючи p можна знайти розв’язок розв’язку однорідного рівняння. Знаходження розв’язку для усталеного режиму (статика) досягається прирівнюванням всіх похідних до нуля в неоднорідному диференціальному рівнянні.
Двигун. Однорідне диференціальне рівняння:
,
нехай, тоді , скоротивши на отримаємо:
дискримінант , , тоді , .
Отже , де - сталі.
Частковий розв’язок неоднорідного рівняння .
Отже, .
При подачі на вхід одиничної стрибкоподібної функції розв’язок записується як. Перша похідна .
Тоді при нульових незалежних початкових умовах:
Розв’язуючи систему рівнянь отримаємо: .Таким чином при одиничній стрибкоподібній функції на вході ланка описується виразом:
.
Перехідні характеристики ланок системи
Рис.7.1 Перехідна характеристика двигуна, як розв'язок диференціального рівняння двигуна: =>
t
W(t)
0
0
0,005
0,00860882
0,01
0,030319985
0,015
0,062153637
0,02
0,101751004
0,025
0,147245882
0,03
0,197162698
0,035
0,250335647
0,04
0,30584456
0,045
0,362964023
0,05
0,421123029
0,055
0,479872969
0,06
0,538862245
0,065
0,597816131
0,07
0,656520803
0,075
0,714810674
0,08
0,772558341
0,085
0,829666623
0,09
0,88606224
0,095
0,941690799
0,1
0,996512824
0,105
1,050500602
0,11
1,103635688
0,115
1,155906923
0,12
1,207308864
0,125
1,257840538
Інші ланки системи описуються алгебраїчними рівняннями, тобто є безінерційними, тому їх перехідні характеристики є стрибкоподібною функцію (повторює форму вхідного сигналу) з амплітудою рівною коефіцієнту передачі відповідної ланки.
Вирази диференціальних рівнянь розімкненої і замкненої системи
Диференціальне рівняння замкненої системи.
Підставимо рівняння редуктора у рівняння двигуна отримуємо рівняння . Отримане рівняння підставимо і рівняння електронного підсилювача , отримуємо рівняння . Це отримане рівняння підставляємо у рівняння вимірювальної схеми , отримуємо рівняння . Розкриваємо усі дужки у отриманому рівнянні . Замінюємо на і отримуємо таке рівняння . Замінюємо коефіцієнти на відповідні a і b, отримуємо рівняння .
Числові значення коефіцієнтів a і b будуть такими:
Диференціальне рівняння розімкненої системи.
Для розімкненої системи рівняння вимірювальної схеми набуде вигляду , тому диференційне рівняння майже таким як і для замкненої системи тільки без . Отже диференціальне рівняння набуде такого вигляду . Де коефіцієнти будуть аналогічними до коефіцієнтів замкненої системи.
9.Схеми електронного моделювання ланок і замкненої системи.
Електронне моделювання ланок і систем здійснюється при певній комбінації включення типових ланок САУ. Безінерційні елементи системи (вимірювальний перетворювач, електронний підсилювач, редуктор) моделюються за допомогою безінерційних ланок з відповідними коефіцієнтами передачі.
Рис.9.1 Схема електронного моделювання ВП1 вимірювального перетворювача;
Рис.9.2 Схема електронного моделювання ЕП електронного підсилювача;
Рис.9.3 Схема електронного моделювання Ред. редуктора;
Рис.9.4 Схема електронного моделювання ВП2 вимірювального перетворювача;
Двигун описується рівнянням третього порядку. Його моделювання здійснюється за допомогою коливної і інтегруючої ланок:
Рис.9.5 Схема електронного моделювання ДВ двигуна;
Схему порівняння можна промоделювати за допомогою підсумовуючої безінерційної ланки при умові, що сигнали на нього подаються у протилежних фазах. Оскільки у всіх ланках використовується інверсне включення операційних підсилювачів, то кожна моделююча ланка змінює знак (фазу) сигналу на протилежний.
Схема електронного моделювання замкненої системи буде мати такий вигляд:
Рис.9.6 Схема електронного моделювання замкненої системи;
Визначення стійкості системи і визначення граничного коефіцієнта підсилення
Визначення стійкості системи за допомогою критерію стійкості Гурвіца.
По критерію стійкості системи Гурвіца система автоматичного управління буде стійкою, якщо визначник Гурвіца і усі діагональні мінори додатні.
Приведемо диференційне рівняння замкненої системи до операторного вигляду і прирівняємо до нуля. Кінцеве рівняння набуде вигляду . Де:
Складаємо визначник Гурвіца і його діагональні мінори:
Визначення коефіцієнта підсилення.
Граничний коефіцієнт підсилення визначається з умови , де . Звідси визначаємо . Таким чином коефіцієнт підсилення при якому система залишається на межі стійкості рівний 600.
Вирази передаточних функцій для окремих ланок системи
перший і другий вимірювальний перетворювач
електронний підсилювач
двигун
редуктор
Вирази передаточних функцій для розімкненої і замкненої системи
Вираз передаточної функції для розімкнена система.
Вираз передаточної функції для замкненої система.
.
Вирази передаточних функцій для комплексних коефіцієнтів передачі окремих ланок системи, розімкненої і замкненої системи.
Вирази передаточних функцій для комплексних коефіцієнтів передачі окремих ланок системи.
перший і другий вимірювальний перетворювач
електронний підсилювач
двигун
редуктор
Вирази передаточної функцій для розімкненої системи.
Вираз передаточної функції для замкненої системи.
.
Аналітичний розрахунок і побудова АФХ, ЛАЧХ і ФЧХ окремих ланок системи і розімкненої системи.
Без інерційні ланки.
Передаточні функції для без інерційних ланок, а саме: перший і другий вимірювальний перетворювач, електронний підсилювач і редуктор мають вигляд дійсного числа і не залежать від частоти.
Амплітудно фазова характеристика для без інерційних ланок являє собою одну точку з координатами , де K – це коефіцієнт передачі відповідної ланки.
Амплітудна частотна характеристика для без інерційних ланок – це пряма, яка паралельна осі частот з координатою рівня K (для лінійної АФХ) або , тому що , де K – це коефіцієнт передачі відповідної ланки.
АФХ для без інерційних ланок:
a) перший і другий вимірювальний перетворювач
U(w) = S = 10, V(w)=0 – вимірювальні перетворювачі мають тільки дійсну складову передаточної функції. A(w) = S = 10
L(w) =20lg S = 20lg 10 = 20
b) електронний підсилювач
U(w) = Кп = 10, V(w)=0 – електронний підсилювач має тільки дійсну складову передаточної функції.
A(w) = Кп = 10
L(w) =20lg Кп = 20lg 10 = 20
c) редуктор
U(w) = 1/і = 0,0833; V(w)=0 – аналогічно, редуктор, як безінерційна ланка в даній системі має тільки дійсну складову передаточної функції.
ЛАЧХ для безінерційних ланок:
перший і другий вимірювальний перетворювач
U(w) = S = 10, V(w)=0 – вимірювальні перетворювачі мають тільки дійсну складову передаточної функції.
A(w) = S = 10
L(w) =20lg S = 20lg 10 = 20
електронний підсилювач
U(w) = Кп = 10, V(w)=0 – електронний підсилювач має тільки дійсну складову передаточної функції.
A(w) = Кп = 10
L(w) =20lg Кп = 20lg 10 = 20
редуктор
U(w) = 1/і = 0,0833; V(w)=0 – аналогічно, редуктор, як безінерційна ланка в даній системі має тільки дійсну складову передаточної функції.
ФЧХ для безінерційних ланок: Фазочастотні характеристики , це пряма яка проходить по осі абсцис частот:
Інерційні ланки
У нашій системі тільки одна інерційна ланка – це двигун, і тому тільки для нього я будую АФХ, ЛАЧХ і ФЧХ.
З комплексного коефіцієнта передачі двигуна я знаходжу і :
=U(w)+jV(w)
Тоді
АФХ
w
V(w)
U(w)
0,1
-39,96642604
-1,199064725
0,2
-19,93300787
-1,196267576
0,3
-13,23323218
-1,191634379
0,4
-9,867247656
-1,185207518
0,5
-7,835196048
-1,177044975
0,6
-6,470543613
-1,167219043
0,7
-5,487698446
-1,155814769
0,8
-4,743913868
-1,142928172
0,9
-4,159922298
-1,128664329
1
-3,688188521
-1,113135368
1,1
-3,298479305
-1,096458459
1,2
-2,970648353
-1,07875384
1,3
-2,690751546
-1,060142953
1,4
-2,44882745
-1,04074671
1,5
-2,237566032
-1,020683948
1,6
-2,051477057
-1,000070064
1,7
-1,886352431
-0,979015875
1,8
-1,738908188
-0,957626684
1,9
-1,606539928
-0,936001552
2
-1,487151982
-0,914232776
2,1
-1,379035699
-0,892405541
2,2
-1,280781192
-0,870597745
2,3
-1,191212339
-0,848879964
2,4
-1,109338209
-0,827315536
2,5
-1,034316298
-0,805960751
ЛАЧХ
Якщо для двигуна з рівняння передаточної функції, то
lg(w)
L(w)
0
11,71486719
0,301029996
4,839355323
0,602059991
-3,535731201
0,903089987
-13,90189883
1,204119983
-25,71984383
1,505149978
-38,80340588
1,806179974
-53,74131743
2,10720997
-70,49421609
2,408239965
-88,15750873
2,709269961
-106,1137723
3,010299957
-124,1487872
3,311329952
-142,2038652
3,612359948
-160,2639829
3,913389944
-178,325362
ЛФЧХ
Якщо, то
lg(w)
Fi(w)
0
1,277679017
0,301029996
0,682818311
0,602059991
0,18615286
0,903089987
0,032072332
1,204119983
-0,006979053
1,505149978
-0,016743227
1,806179974
-0,019183844
2,10720997
-0,019793964
2,408239965
-0,019946492
2,709269961
-0,019984623
3,010299957
-0,019994156
3,311329952
-0,01999654
3,612359948
-0,019997135
3,913389944
-0,019997284
Розімкнена система.
З комплексного рівняння розімкненої я знаходжу і :
Тоді
АФХ
w
V(w)
U(w)
0,1
-333,020245
-9,991206823
0,2
-166,0917881
-9,967899579
0,3
-110,2659072
-9,929293464
0,4
-82,21884109
-9,875741643
0,5
-65,28677107
-9,807727251
0,6
-53,91580466
-9,72585268
0,7
-45,7262473
-9,630826564
0,8
-39,5286623
-9,523448997
0,9
-34,66255255
-9,404595522
1
-30,73183085
-9,275200457
1,1
-27,48457881
-9,136240109
1,2
-24,7529274
-8,988716373
1,3
-22,42068726
-8,833641152
1,4
-20,40485472
-8,672021963
1,5
-18,64451896
-8,504848993
1,6
-17,09393258
-8,333083805
1,7
-15,71803163
-8,15764978
1,8
-14,48945248
-7,979424346
1,9
-13,38649395
-7,799232936
2
-12,39169389
-7,617844605
2,1
-11,49081496
-7,435969168
2,2
-10,67210928
-7,254255713
2,3
-9,925776816
-7,073292301
2,4
-9,24356063
-6,893606701
2,5
-8,61844055
-6,715667961
2,6
-8,044398421
-6,539888652
2,7
-7,51623623
-6,366627639
З комплексного рівняння розімкненої я знаходжу і :
Тоді
ЛАЧХ
lg(w)
L(w)
0
30,13037363
0,301029996
23,25486177
0,602059991
14,87977525
0,903089987
4,513607621
1,204119983
-7,304337382
1,505149978
-20,38789944
1,806179974
-35,32581098
2,10720997
-52,07870964
2,408239965
-69,74200229
2,709269961
-87,69826582
3,010299957
-105,7332808
3,311329952
-123,7883588
3,612359948
-141,8484764
3,913389944
-159,9098556
Тоді ЛФЧХ
lg(w)
Fi(w)
0
1,277679017
0,301029996
0,682818311
0,602059991
0,18615286
0,903089987
0,032072332
1,204119983
-0,006979053
1,505149978
-0,016743227
1,806179974
-0,019183844
2,10720997
-0,019793964
2,408239965
-0,019946492
2,709269961
-0,019984623
3,010299957
-0,019994156
3,311329952
-0,01999654
3,612359948
-0,019997135
3,913389944
-0,019997284
Визначення стійкості системи по АФХ розімкненої системи, запаси стійкості по амплітуді і по фазі
Визначення стійкості системи за виглядом частотних характеристик здійснюється за допомогою частотних критеріїв стійкості, зокрема за критерієм Найквіста: САУ буде стійка у замкненому стані, якщо розімкнена системи на комплексній площині не охоплює точку з координатами (-1; j0).
АФХ розімкненої системи (фрагмент):
З комплексного рівняння розімкненої я знаходжу і :
Тоді
w
V(w)
U(w)
12,1
-0,025378698
-0,757978217
12,2
-0,021795356
-0,745802181
12,3
-0,0183493
-0,733892464
12,4
-0,015035052
-0,722241643
12,5
-0,011847376
-0,710842544
12,6
-0,008781272
-0,699688231
12,7
-0,005831964
-0,688771997
12,8
-0,002994886
-0,678087359
12,9
-0,000265671
-0,667628042
13
0,002359854
-0,657387978
13,1
0,004885683
-0,647361298
13,2
0,007315637
-0,637542318
13,3
0,009653372
-0,627925539
13,4
0,011902387
-0,618505638
13,5
0,014066035
-0,60927746
13,6
0,016147526
-0,600236015
13,7
0,018149934
-0,591376468
13,8
0,020076206
-0,582694136
13,9
0,021929166
-0,574184481
14
0,023711521
-0,565843109
14,1
0,025425867
-0,557665758
14,2
0,027074695
-0,549648297
14,3
0,028660391
-0,541786723
14,4
0,030185249
-0,534077153
14,5
0,031651468
-0,526515823
12,1
-0,025378698
-0,757978217
12,2
-0,021795356
-0,745802181
АФХ
Як видно з малюнку АФХ розімкненої системи не охоплює точку з координатами (-1; j0), отже система в замкненому стані буде стійкою.
Перевіримо за таблицею – V(w) переходить межу нуля, коли U(w) лежить в межах -0.755.. -0.742, отже, запас стійкості по амплітуді – це відстань між -1 до точки перетину характеристики з віссю абсцис, A=0.25.
Запас стійкості по фазі дорівнює φ ≈13 градусів – це кут між радіусом одиничного кола, проведеного в точку перетину з графіком і віссю абсцис.
Визначення по ЛАЧХ і ФЧХ розімкненої системи стійкість системи і значення запасів стійкості по амплітуді і по фазі
Для того щоб САУ була стійкою, необхідно і достатньо, щоб при досягненні ЛАЧХ значення , .
ЛАЧХ і ФЧХ (фрагменти), визначення стійкості системи:
lg(w)
L(w)
Fi(w)
0
30,13037363
1,277679017
0,301029996
23,25486177
0,682818311
0,602059991
14,87977525
0,18615286
0,903089987
4,513607621
0,032072332
1,204119983
-7,304337382
-0,006979053
1,505149978
-20,38789944
-0,016743227
1,806179974
-35,32581098
-0,019183844
2,10720997
-52,07870964
-0,019793964
2,408239965
-69,74200229
-0,019946492
2,709269961
-87,69826582
-0,019984623
3,010299957
-105,7332808
-0,019994156
3,311329952
-123,7883588
-0,01999654
3,612359948
-141,8484764
-0,019997135
3,913389944
-159,9098556
-0,019997284
Як видно з малюнка, при досягненні фази значення 0, амплітуда сигналу L(w) також дорівнює нулю, що свідчить про те, що система знаходиться на межі стійкості.
Для характеристик в логарифмічному масштабі вводяться коефіцієнти: запаси стійкості по амплітуді L – як логарифм амплітуди, взятої при частоті, відповідної точки пересічення фазової характеристики лінією ; запас стійкості по фазі - як ордината фазової характеристики відраховуючи від лінії , при частоті зрізу.
Запас стійкості по амплітуді визначається, як , A=0.
Запас стійкості по фазі визначається, як fi= 0.
Графік перехідного процесу при одиничній стрибкоподібній дії вхідної величини
Для отримання функції перехідного процесу при подачі на вхід одиничної стрибкоподібної функції скористаємося математичним пакетом MatLab 6.1. Введемо необхідні дані передаточної функції замкненої системи в оператор STEP за принципом: чисельник – 25, а знаменник – це коефіцієнти 0.006 0.3 1 33,33.
.
Визначення якісних показників системи
За графіком перехідного процесу визначаємо наступні якісні показники системи.
Час перехідного процесу – це час, по закінченні якого регульована величина не виходить за межі проміжку 2%.
Час перехідного процесу – 7.13 с.
Максимальне значення – 1.83.
Статичне значення .
Пере регулювання – відношення максимального відхилення регульованої величини до свого усталеного значення у відсотках.
.
Використана література.
П.И. Чинаєв, Н.М. Чумаков. Теория автоматического управления. Києвскоє высшее инженерно-авиационноэ военное училище ВВС. Киэв-1969.
Теория автоматического управления: Учеб. для вузов по спец. «Автоматика и телемеханика». В 2-х ч. Ч. І. Теория линейных систем автоматического управления/ Н.А. Бабаков, А.А. Воронов, и др.; Под ред. А.А. Воронова. - 2-е изд., перераб. и доп. – М.:Высш. шк.,1986. – 367 с.,ил.