Міністерство освіти і науки України
НУ ”Львівська політехніка”
“ІКТА”
Кафедра “Автоматики і телемеханіки”
КУРСОВА РОБОТА
з курсу
“ТЕОРІЯ АВТОМАТИЧНОГО КЕРУВАННЯ”
Тема № 6
Варіант № 8
Львів-2006р.
Тема № 6. Варіант № 8.
Завдання.
1.Описати призначення і принцип роботи схеми.
2.Зобразити структурну схему системи.
3.Записати вирази коефіцієнтів передачі окремих ланок, розімкненої і замкнутої системи.
4.Зобразити статичні характеристики окремих ланок системи.
5.Записати вирази диференціальних рівнянь для окремих ланок системи.
6.Знайти рішення диференціальних рівнянь для окремих ланок системи.
7.За результатами рішень диференціальних рівнянь побудувати перехідні характеристики окремих ланок системи.
8.Записати вирази диференціальних рівнянь розімкненої і замкнутої системи.
9.Привести схеми електронного моделювання окремих ланок і замкненої системи.
10.Користуючись одним з алгебраїчних критеріїв стійкості визначити стійкість системи і знайти граничний коефіцієнт підсилення.
11.Записати вирази передаточних функцій для окремих ланок системи.
12.Записати вирази передаточних функцій розімкненої і замкнутої системи.
13.Записати вирази для комплексних коефіцієнтів передачі окремих ланок системи, розімкненої і замкнутої системи.
14.Розрахувати аналітично і побудувати АФХ, ЛАЧХ і ФЧХ окремих ланок системи і розімкненої системи.
15.По АФХ розімкненої системи визначити стійкість системи. Знайти запаси стійкості по амплітуді і фазі.
16.По ЛАЧХ і ФАХ розімкненої системи визначити стійкість системи і знайти запаси стійкості по амплітуді і фазі.
17.Побудувати графік перехідного процесу при одиничній стрибкоподібній дії вхідної величини.
18.По графіку перехідного процесу визначити якісні показники системи.
Рівняння ланок:
А) вимірювальна схема
Б) електронний підсилювач
В) двигун
Г) редуктор
Вихідні дані:
Параметри
Розмірність
8
TM
Сек
0.3
TE
Сек
0.03
C
Рад/сек
3
Kn
100
S
В/рад
1
і
12
1.Опис роботи системи.
Слідкуюча система витратоміра літака здійснює контроль над витратами літака. Під впливом газового потоку заслонка зміщується. Так як заслонка зв‘язана з повзунком реостата R1, то повзунок також переміщається на певну відстань. В результаті цього (вх стає відмінним від (вих, що приводить до утворення похибки розузгодження (U. Значення (U подається на електронний підсилювач, де підсилюється, і підсилена напруга Up поступає на двигун. Вал двигуна, при надходженні цієї напруги, починає обертатись. Оберти з вала двигуна подаються на редуктор. Редуктор, який зв‘язаний з повзунком реостата R2, переміщує його доти, поки (U=0 і (вх=(вих.
2.Структурна схема системи.
(вх Uвх (U Up ( (вих
Uвих Рис.2.1 Структурна схема системи
В структурній схемі використані такі позначення:
ЕП-електронний підсилювач;
ВП1,ВП2-вимірювальні перемворювачі;
ДВ-двигун;
Ред-редуктор;
3.Коефіцієнти передачі ланок, розімкненої і замкненої системи.
Коефіцієнтом передачі називається відношення вихідної величини до вхідної в усталеному режимі, тобто коли всі часові похідні дорівнюють нулю.
Коефіцієнт передачі вимірювальних перетворювачів(ВП1,ВП2):
(В/рад);
(В/рад);
Коефіцієнт передачі електричного підсилювача:
;
Коефіцієнт передачі двигуна:
(Рад/в.сек);
Коефіцієнт передачі редуктора:
;
Для розімкненої системи (розімкнена система – система, в якій відсутній контроль за регульованою (вихідною) величиною, в даному випадку для розімкненої системи ) коефіцієнт передачі є добутком коефіцієнтів передачі ланок ВП1, ЕП, ДВ, Ред., оскільки для розімкненої системи вони включені послідовно. Отже коефіцієнт передачі розімкненої системи:
=25;
Для замкненої системи (замкнена система – система, в якій відбувається постійний контроль за регульованою (вихідною) величиною) коефіцієнт передачі є добутком коефіцієнта передачі ланки ВП1 та коефіцієнта передачі ланок ЕП, ДВ, Ред з врахуванням охоплення їх від’ємним оберненим зв’язком:
.
Статичні характеристики окремих ланок
Ланки системи:
вимірювальна схема
b) електронний підсилювач
двигун
редуктор
Вирази диференціальних рівнянь для окремих ланок системи
У даній системі перший і другий вимірювальні перетворювачі, а також електронний підсилювач, двигун і редуктор є без інерційними елементами, тому вони описуються алгебраїчними рівняннями. Двигун же є інерційним елементом, тому він описується диференційним рівнянням:
де:
;
підставивши ці значення у диференційне рівняння отримуємо остаточне рівняння:
6.Розв’язок диференціальних рівнянь ланок системи.
Розв’язок неоднорідного диференціального рівняння є сумою повного розв’язку однорідного диференціального рівняння (вільна складова) і часткового розв’язку неоднорідного диференціального рівняння (усталений режим). При розв’язку однорідного диференціального рівняння шукана функція замінюється виразом , де - змінна інтегрування, p – стала. Після відповідних підстановок рівняння зводиться до алгебраїчного і розв’язується відносно p. Маючи p можна знайти розв’язок розв’язку однорідного рівняння. Знаходження розв’язку для усталеного режиму (статика) досягається прирівнюванням всіх похідних до нуля в неоднорідному диференціальному рівнянні.
Двигун. Однорідне диференціальне рівняння:
,
нехай, тоді , скоротивши на отримаємо:
дискримінант , , тоді , .
Отже , де - сталі.
Частковий розв’язок неоднорідного рівняння .
Отже, .
При подачі на вхід одиничної стрибкоподібної функції розв’язок записується як . Перша похідна .
Тоді при нульових незалежних початкових умовах:
Розв’язуючи систему рівнянь отримаємо: .Таким чином при одиничній стрибкоподібній функції на вході ланка описується виразом:
.
7.Перехідні характеристики ланок системи.
Рис.7.1 Перехідна характеристика двигуна.
t
w(t)
0
0
0,001
3,41179E-06
0,002
1,21713E-05
0,003
2,61032E-05
0,004
4,50379E-05
0,005
6,88113E-05
0,006
9,72648E-05
0,007
0,000130245
0,008
0,000167603
0,009
0,000209195
0,01
0,000254883
0,011
0,000304533
0,012
0,000358013
0,013
0,000415199
0,014
0,000475969
0,015
0,000540206
0,016
0,000607794
0,017
0,000678626
0,018
0,000752593
0,019
0,000829594
0,02
0,000909528
0,021
0,000992299
0,022
0,001077815
0,023
0,001165984
0,024
0,00125672
0,025
0,001349939
Інші ланки системи описуються алгебраїчними рівняннями, тобто є безінерційними, тому їх перехідні характеристики є стрибкоподібною функцію (повторює форму вхідного сигналу) з амплітудою рівною коефіцієнту передачі відповідної ланки.
Вирази диференціальних рівнянь розімкненої і замкненої системи
Диференціальне рівняння замкненої системи.
Підставимо рівняння редуктора у рівняння двигуна отримуємо рівняння . Отримане рівняння підставимо і рівняння електронного підсилювача , отримуємо рівняння . Це отримане рівняння підставляємо у рівняння вимірювальної схеми , отримуємо рівняння . Розкриваємо усі дужки у отриманому рівнянні . Замінюємо на і отримуємо таке рівняння . Замінюємо коефіцієнти на відповідні a і b, отримуємо рівняння .
Числові значення коефіцієнтів a і b будуть такими:
Диференціальне рівняння розімкненої системи.
Для розімкненої системи рівняння вимірювальної схеми набуде вигляду , тому диференційне рівняння майже таким як і для замкненої системи тільки без . Отже диференціальне рівняння набуде такого вигляду . Де коефіцієнти будуть аналогічними до коефіцієнтів замкненої системи.
9.Схеми електронного моделювання ланок і замкненої системи.
Електронне моделювання ланок і систем здійснюється при певній комбінації включення типових ланок САУ. Безінерційні елементи системи (вимірювальний перетворювач, електронний підсилювач, редуктор) моделюються за допомогою безінерційних ланок з відповідними коефіцієнтами передачі.
Рис.9.1 Схема електронного моделювання ВП1 вимірювального перетворювача;
Рис.9.2 Схема електронного моделювання ЕП електронного підсилювача;
Рис.9.3 Схема електронного моделювання Ред. редуктора;
Рис.9.4 Схема електронного моделювання ВП2 вимірювального перетворювача;
Двигун описується рівнянням третього порядку. Його моделювання здійснюється за допомогою коливної і інтегруючої ланок:
Рис.9.5 Схема електронного моделювання ДВ двигуна;
Схему порівняння можна промоделювати за допомогою підсумовуючої безінерційної ланки при умові, що сигнали на нього подаються у протилежних фазах. Оскільки у всіх ланках використовується інверсне включення операційних підсилювачів, то кожна моделююча ланка змінює знак (фазу) сигналу на протилежний.
Схема електронного моделювання замкненої системи буде мати такий вигляд:
Рис.9.6 Схема електронного моделювання замкненої системи;
Визначення стійкості системи і визначення граничного коефіцієнта підсилення
Визначення стійкості системи за допомогою критерію стійкості Гурвіца.
По критерію стійкості системи Гурвіца система автоматичного управління буде стійкою, якщо визначник Гурвіца і усі діагональні мінори додатні.
Приведемо диференційне рівняння замкненої системи до операторного вигляду і прирівняємо до нуля. Кінцеве рівняння набуде вигляду . Де:
Складаємо визначник Гурвіца і його діагональні мінори:
Визначення коефіцієнта підсилення.
Граничний коефіцієнт підсилення визначається з умови , де . Звідси визначаємо . Таким чином коефіцієнт підсилення при якому система залишається на межі стійкості рівний 400.
Вирази передаточних функцій для окремих ланок системи
перший і другий вимірювальний перетворювач
електронний підсилювач
двигун
редуктор
Вирази передаточних функцій для розімкненої і замкненої системи
Вираз передаточної функції для розімкнена система.
Вираз передаточної функції для замкненої система.
.
Вирази передаточних функцій для комплексних коефіцієнтів передачі окремих ланок системи, розімкненої і замкненої системи.
Вирази передаточних функцій для комплексних коефіцієнтів передачі окремих ланок системи.
перший і другий вимірювальний перетворювач
електронний підсилювач
двигун
редуктор
Вирази передаточної функцій для розімкненої системи.
Вираз передаточної функції для замкненої системи.
.
Аналітичний розрахунок і побудова АФХ, ЛАЧХ і ФЧХ окремих ланок системи і розімкненої системи.
Без інерційні ланки.
Передаточні функції для без інерційних ланок, а саме: перший і другий вимірювальний перетворювач, електронний підсилювач і редуктор мають вигляд дійсного числа і не залежать від частоти.
Амплітудно фазова характеристика для без інерційних ланок являє собою одну точку з координатами , де K – це коефіцієнт передачі відповідної ланки.
Амплітудна частотна характеристика для без інерційних ланок – це пряма, яка паралельна осі частот з координатою рівня K (для лінійної АФХ) або , тому що , де K – це коефіцієнт передачі відповідної ланки.
перший і другий вимірювальний перетворювач
електронний підсилювач
редуктор
Фазочастотні характеристики , це пряма яка проходить по осі частот.
АФХ для без інерційних ланок:
ЛАЧХ для без інерційних ланок:
ФЧХ для без інерційних ланок:
Інерційні ланки.
У нашій системі тільки одна інерційна ланка – це двигун, і тому тільки для нього я будую АФХ, ЛАЧХ і ФЧХ.
З комплексного коефіцієнта передачі двигуна я знаходжу і :
Тоді
АФХ
w
V(w)
U(w)
0
#ДЕЛ/0!
#ДЕЛ/0!
1
-8,31933408
-2,518466422
2
-3,364626897
-2,094166119
3
-1,666303025
-1,6318528
4
-0,886439954
-1,242672833
5
-0,489366367
-0,94716071
6
-0,274278221
-0,730326624
7
-0,152189979
-0,571733375
8
-0,080306059
-0,454562596
9
-0,036803446
-0,366676395
10
-0,009988901
-0,299667037
11
0,006681842
-0,247753677
12
0,017014603
-0,206934355
13
0,023298368
-0,17440237
14
0,026950782
-0,148158749
15
0,028872392
-0,126756844
16
0,029647844
-0,109133169
17
0,029663748
-0,094494138
18
0,029179689
-0,082239207
19
0,028372104
-0,071907956
20
0,027362021
-0,063143124
21
0,026232957
-0,055664409
22
0,025042696
-0,04924964
23
0,023831152
-0,043721071
24
0,022625721
-0,038935275
ЛАЧХ
w
lg(w)
L(w)
1
0
9,240149102
2
0,301029996
2,418298828
4
0,602059991
-5,868844334
8
0,903089987
-16,25707532
16
1,204119983
-28,47404454
32
1,505149978
-42,59260797
64
1,806179974
-58,76880402
128
2,10720997
-76,20466313
256
2,408239965
-94,09534439
512
2,709269961
-112,1133249
1024
3,010299957
-130,1641025
2048
3,311329952
-148,2231425
4096
3,612359948
-166,2842521
8192
3,913389944
-184,3458792
ФЧХ
w
lg(w)
Fi(w)
1
0
-0,321447799
2
0,301029996
-0,600347655
4
0,602059991
-0,995486977
8
0,903089987
-1,411550188
16
1,204119983
-1,812920913
32
1,505149978
-2,231996819
64
1,806179974
-2,609381309
128
2,10720997
-2,860798569
256
2,408239965
-2,999092686
512
2,709269961
-3,070069911
1024
3,010299957
-3,105796864
2048
3,311329952
-3,123690446
4096
3,612359948
-3,13264101
8192
3,913389944
-3,137116765
Розімкнена система.
З комплексного рівняння розімкненої я знаходжу і :
Тоді
АФХ
w
V(w)
U(w)
0
#ДЕЛ/0!
#ДЕЛ/0!
1
-23,10926133
-6,995740061
2
-9,346185825
-5,817128107
3
-4,628619515
-4,532924443
4
-2,462333206
-3,45186898
5
-1,35935102
-2,631001973
6
-0,761883947
-2,028685066
7
-0,422749943
-1,588148264
8
-0,223072385
-1,262673879
9
-0,102231793
-1,018545541
10
-0,027746948
-0,832408435
11
0,018560671
-0,688204658
12
0,047262785
-0,574817653
13
0,064717689
-0,484451028
14
0,074863283
-0,41155208
15
0,08020109
-0,352102344
16
0,082355123
-0,303147692
17
0,0823993
-0,262483717
18
0,081054691
-0,228442241
19
0,0788114
-0,199744322
20
0,076005613
-0,175397568
21
0,072869326
-0,154623359
22
0,069563043
-0,136804555
23
0,066197643
-0,121447418
24
0,062849225
-0,108153543
25
0,05956933
-0,096598913
26
0,056392059
-0,086518107
ЛАЧХ
w
lg(w)
L(w)
1
0
27,65652418
2
0,301029996
20,83467391
4
0,602059991
12,54753074
8
0,903089987
2,159299761
16
1,204119983
-10,05766946
32
1,505149978
-24,17623289
64
1,806179974
-40,35242894
128
2,10720997
-57,78828805
256
2,408239965
-75,67896931
512
2,709269961
-93,69694978
1024
3,010299957
-111,7477275
2048
3,311329952
-129,8067674
4096
3,612359948
-147,867877
8192
3,913389944
-165,9295042
ФЧХ
w
lg(w)
Fi(w)
1
0
-0,321447799
2
0,301029996
-0,600347655
4
0,602059991
-0,995486977
8
0,903089987
-1,411550188
16
1,204119983
-1,812920913
32
1,505149978
-2,231996819
64
1,806179974
-2,609381309
128
2,10720997
-2,860798569
256
2,408239965
-2,999092686
512
2,709269961
-3,070069911
1024
3,010299957
-3,105796864
2048
3,311329952
-3,123690446
4096
3,612359948
-3,13264101
8192
3,913389944
-3,137116765
Визначення стійкості системи по АФХ розімкненої системи, запаси стійкості по амплітуді і по фазі
Визначення стійкості системи за виглядом частотних характеристик здійснюється за допомогою частотних критеріїв стійкості, зокрема за критерієм Найквіста: САУ буде стійка у замкненому стані, якщо розімкнена системи на комплексній площині не охоплює точку з координатами (-1; j0).
АФХ розімкненої системи (фрагмент):
Як видно з малюнку АФХ розімкненої системи не охоплює точку з координатами (-1; j0), отже система в замкненому стані буде стійкою.
Запас стійкості по амплітуді – це відстань між -1 до точки перетину характеристики з віссю абсцис. A=0.21
Запас стійкості по фазі дорівнює φ ≈13 градусів.
Визначення по ЛАЧХ і ФЧХ розімкненої системи стійкість системи і значення запасів стійкості по амплітуді і по фазі
Для того щоб САУ була стійкою, необхідно і достатньо, щоб при досягненні ЛАЧХ значення , .
ЛАЧХ і ФЧХ (фрагменти), визначення стійкості системи:
Як видно з малюнка, при досягненні фази значення , амплітуда сигналу , що свідчить про стійкість системи.
Для характеристик в логарифмічному масштабі вводяться коефіцієнти: запаси стійкості по амплітуді L – як логарифм амплітуди, взятої при частоті, відповідної точки пересічення фазової характеристики лінією ; запас стійкості по фазі - як ордината фазової характеристики відраховуючи від лінії , при частоті зрізу.
Запас стійкості по амплітуді визначається, як , A=90.
Запас стійкості по фазі визначається, як fi= 1.5.
Графік перехідного процесу при одиничній стрибкоподібній дії вхідної величини
Для отримання функції перехідного процесу при подачі на вхід одиничної стрибкоподібної функції скористаємося математичним пакетом MatLab 6.1. Введемо необхідні дані передаточної функції замкненої системи в оператор STEP за принципом: чисельник – 25, а знаменник – це коефіцієнти 0.009 0.3 1 25.
.
Визначення якісних показників системи
За графіком перехідного процесу визначаємо наступні якісні показники системи.
Час перехідного процесу – це час, по закінченні якого регульована величина не виходить за межі проміжку .
Час перехідного процесу – 9.61 с.
Максимальне значення – 1.84.
Статичне значення .
Пере регулювання – відношення максимального відхилення регульованої величини до свого усталеного значення у відсотках.
.
Використана література.
П.И. Чинаєв, Н.М. Чумаков. Теория автоматического управления. Києвскоє высшее инженерно-авиационноэ военное училище ВВС. Киэв-1969.
Теория автоматического управления: Учеб. для вузов по спец. «Автоматика и телемеханика». В 2-х ч. Ч. І. Теория линейных систем автоматического управления/ Н.А. Бабаков, А.А. Воронов, и др.; Под ред. А.А. Воронова. - 2-е изд., перераб. и доп. – М.:Высш. шк.,1986. – 367 с.,ил.