Розробка лінійних систем автоматичного управління (САУ) з цифровою корегуючою ланкою

Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

Вступ

Мета курсової роботи - набуття студентами навичок у розрахунку лінійних систем автоматичного управління (САУ) з цифровою корегуючою ланкою, роль якої може виконувати мікропроцесор, що управляє обчислювальною машиною, чи будь-який спеціалізований цифровий управляючий пристрій.

Відповідно до завдання необхідно розробити слідкуючу систему, що задовольняє визначеним технічним умовам. Система повинна забезпечувати синхронне і синфазне обертання двох вісей, механічно не зв'язаних між собою. Входом системи є кут повороту сельсина-датчика, а виходом - кут повороту вихідного вала редуктора, механічно зв'язаного з робочим механізмом і з ротором сельсина-приймача.

Слідкуючі системи, розглянутого типу, широко застосовуються для дистанційного управління різними механізмами, а також при побудові САУ в різних галузях промисловості.

1. Завдання на виконання курсової роботи

1.1. Мета курсової роботи - проектування слідкуючої системи, яка задовольняє заданим технічним умовам.

1.2. Вимірювальний пристрій - сельсинні пари.

1.3. Виконавчий двигун - двигун постійного струму серії МІ.

1.4. Підсилювач потужності - електромашинний підсилювач з поперечним полем.

1.5. Вихідні дані для проектування системи.

1.5.1. Статичний момент навантаження об'єкта управління - Мос, Нм

1.5.2. Момент інерції об'єкта управління - Jо, кгм²

1.5.3. Максимальна кутова швидкість об'єкта управління - ?о max.

1.5.4. Максимальне кутове прискорення об'єкта управління - ?о max

1.5.5. Вимоги, які запропоновані до якості процесу управління: максимальне перерегулювання - ? мах,; час регулювання - tр ; максимальна кінетична похибка - xmax, рад.

Варіанти вихідних даних приведені в додатоку 1.

1.6. Курсова робота повинна містити такі розділи.

1.6.1. Розробка функціональної схеми.

1.6.2. Вибір елементів системи: виконавчого двигуна (серії МІ) і електромашинного підсилювача потужності (ЕМП), розрахунок пере1.6.3. Складання передаточних функцій елементів незкорегованої слідкуючої системи.

1.6.4. Побудова логарифмічних частотних характеристик (ЛАЧХ ) незкорегованої системи, бажаної системи і послідовної корегуючої ланки.

1.6.5. Побудова на ЕОМ перехідної функції H(t) і визначення по ній показників якості перехідного процесу для системи з безперервною послідовною корегуючою ланкою.

1.7. Розрахунково-пояснювальна записка повинна містити такі розділи.

1.7.1. Вступ (мета виконання роботи, опис слідкуючої системи, принцип її роботи, опис функціональної схеми системи).

1.7.2. Вихідні дані для проектування системи.

1.7.3. Розрахункова частина.

1.7.4. Висновок (основні характеристики зпроектованої системи).

1.7.5. Список літератури.

1.8. Додатком до розрахунково-пояснювальної записки повинні бути функціональна і структурна електричні схеми слідкуючої системи; ЛАЧХ і ЛФЧХ зкорегованої і незкорегованої системи, графік перехідній функції системи , таблиці і графіки, необхідні для виконання даної роботи.

1.9. Записка повинна бути оформлена відповідно до вимог ЄСКД (написана чорним чорнилом чи пастою, або надрукована на принтері на аркушах формату 11). Орієнтовний обсяг записки - 20 - 25 сторінок.

2. Порядок розрахунку слідкуючої системи

2.1 Розробка функціональної схеми

В слідкуючій системі, яка проектується, в якості виконавчого двигуна (Д) використовується двигун постійного струму серії МІ, в якості підсилювача потужності - електромашинний підсилювач з поперечним полем (ЕМП). Для вимірювального пристрою (ВП) рекомендується використовувати сельсинну пару: сельсин-датчик і сельсин-трансформатор (приймач). Оскільки вимірювальний пристрій працює на змінному струмі, а підсилювач потужності і виконавчий двигун - на постійному струмі, то після вимірювального пристрою повинний бути застосований фазовий детектор (ФД). Крім зазначених елементів у функціональну схему входять: корегувальний пристрій (КП), підсилювач напруги (П), редуктор (Р), за допомогою якого виконавчий двигун з'єднується з об'єктом управління і ротором сельсина-трансформатора, та об'єкт управління (ОУ).

Функціональна схема цифрової слідкуючої системи, наведена на рис.1. Рис.1. Функціональна схема цифрової слідкуючої системи

Рис.1. Функціональна схема слідкуючої системи

слідкуюча цифрова корегуюча ланка

2.2 Вибір виконавчого двигуна

Вибір двигуна починають з розрахунку необхідної потужності, що повинна бути достатньої для забезпечення заданої швидкості і прискорення об'єкта управління при заданому навантаженні.

Необхідна потужність, Вт:

,

де ?_р???ККД редуктора, ?_{р =} 0,72 .

З каталогу (додаток 2) вибираємо найближчий двигун більшої потужності Р_{н >} Р_тр і виписуємо його паспортні дані:

Р_н - номінальна потужність (Вт);

n_н - номінальна швидкість обертання (об/хв);

U_{н -} номінальна напруга (В);

I_н - номінальний струм якоря (А);

R_д- опір ланцюга обмотки якоря (Ом);

J_д - момент інерції якоря (кг.м²);

?_д- ККД двигуна.

Потім послідовно визначаємо такі величини:

номінальна кутова швидкість двигуна ?_н (с^-1) -

?_н = ?n_н/30 ;

номінальний момент двигуна М_н (Н.м) -

М_н = 9,55Р_н/n_{н ;}

оптимальне передаточне число редуктора i_р -

,

де J_р = 1.10^-4 кг.м² - момент інерції редуктора.

Визначаємо необхідний момент на валу двигуна:

.

Для обраного двигуна необхідно перевірити вимогу чи задовольняє він моменту і швидкості у відповідності до таких умов:

М_тр/М_н ???????????????_{o max} i_{р /} ?_н a ;

де ??????коефіцієнт припустимого перевантаження двигуна по моменту (для двигуна постійного струму ?????????????a - коефіцієнт припустимого короткочасного збільшення швидкості двигуна понад номінальний, за звичай a = 1,2 - 1,5.

Приклад. Вибрати виконавчий двигун слідкуючої системи, якщо відповідно до технічного завдання: Jо=100кг.м²; М_ос= 120 Нм; ?_омах= 0,7 з^-1; ?_omax= 0,44 з^-2 ; ?_р = 0,72.

Необхідна потужність двигуна

Р_тр = 2(120 + 100. 0,44). 0,7/0,72 = 319 Вт .

Вибираємо двигун типу МІ-22 (додаток2) з такими параметрами: Р_н = 370Вт; U_н = 110 В; n_н = 3000 об/хв; J_д = 0,004 кг.м².

Задамо момент інерції редуктора, який підведений до вала двигуна J_р=1.10^-4 кг.м^2, і знайдемо передаточне число редуктора

.

Необхідний момент

.

Номінальний момент обраного двигуна дорівнює:

.

Проведемо перевірку двигуна: по моменту - 1,28/1,18< 10 ;

по швидкості - 0,7.355/314 = 0,79 < 1,5 .

У результаті перевірок двигуна по моменту і швидкості видно, що він не перевантажений. Отже, двигун МІ-22 обраний правильно.

2.3 Вибір підсилювача потужності

В якості підсилювача потужності використовуємо ЕМП з поперечним полем. При виборі підсилювача необхідно дотримуватися таких умов.

1. Номінальна потужність підсилювача повинна задовольняти нерівності

Рпн Р_н / ?_д ,

де ??_д - ККД двигуна.

2. Номінальна напруга підсилювача повинна бути не менше номінальної напруги виконавчого двигуна.

3. Номінальний струм підсилювача повинний бути не менше, ніж номінальний струм двигуна.

З огляду на зазначені умови, вибираємо тип ЕМП (див. додаток 3).

2.4 Складання передаточних функцій елементів слідкуючої системи

2.4.1 Виконавчий двигун

Передаточна функція виконавчого двигуна по куту повороту має вид (якщо знехтувати індуктивністю ланцюга якоря)

де К_д - коефіцієнт підсилення двигуна, рад/В.с : К_д = ?_{н /} U_{н ;}

Т_{д -} електромеханічна постійна часу:

В останній формулі ??=1,2 - постійний коефіцієнт;

J_c - сумарний момент інерції, приведений до вала двигуна

;

2.4.2 Електромашинний підсилювач (ЕМП)

Передаточна функція ЕМП

К_ЕМП - коефіцієнт підсилення ЕМП по напрузі

де U_ЕМП - напруга на виході ЕМП ;

- напруга обмотки управління ЕМП ;

Р_у, R_у - відповідно потужність і опір обмотки управління ЕМП,

Ту, Т_кз - постійні часу обмотки управління і короткозамкненої обмотки

якоря ЕМП.

2.4.3 Підсилювач

Підсилювач (на функціональній схемі - П) слугує для узгодження вихідного сигналу ЦАП із вхідним опором обмотки управління ЕМП. Його можна вважати безінерційною ланкою з передаточною функцією W_п(Р) = К_п.

У розрахунках прийняти К_п= 1 _.

2.4.4 Фазовий детектор

Передаточна функція фазового детектора W_фд (P) = К_фд,

де К_фд=1 - коефіцієнт підсилення фазового детектора.

2.4.5 Вимірювальний пристрій

Передаточна функція вимірювального пристрою W_вп(Р)=К_вп ,

де К_вп=1 - коефіцієнт підсилення вимірювального пристрою.

2.4.6 Редуктор

Передаточна функція редуктора W_ред(Р)=К_ред=1/i_р .

Структурна схема незкорегованої слідкуючої системи, наведена на рис.2.

Рис.2. Структурна схема незкорегованої слідкуючої системи

2.5 Розрахунок послідовної безперервної корегуючої ланки методом ЛАЧХ

Побудова логарифмічної амплітудно-частотної характеристики (ЛАЧХ) послідовної корегуючої ланки проводиться в такій послідовності.

1. Будується ЛАЧХ заданої (незкорегованої) системи.

2. Будується бажана ЛАЧХ за заданими показниками якості перехідного процесу.

3. Будується ЛАЧХ послідовної корегуючої ланки шляхом графічного вирахування ЛАЧХ заданої системи з бажаної ЛАЧХ.

4. По виду ЛАЧХ корегуючої ланки визначається його передаточна функція (безперервна).

2.5.1 Побудова ЛАЧХ заданої системи по виду передаточної функції

Передаточну функцію розімкненої системи потрібно представити у виді добутку передаточних функцій типових динамічних ланок (обмежимося випадком, коли в системі відсутні коливальні ланки і ланки із запізненням). Наприклад, нехай передаточна функція розімкненої системи має вид:

.

Побудову зручно проводити в такій послідовності.

1. Визначити частоти, при яких відбувається зміна нахилу ЛАЧХ і які розраховуються як w_i=1/Т_i, та відкласти їх по вісі абсцис у логарифмічному масштабі (Т_i - постійні часу передаточної функції К_з(Р) ).

2. Відкласти точку A1 з координатами w_A1=1c^-1 і L(w_A1)=20lgК_з ( див. рис. 3). Через точку A1 провести пряму з нахилом -20 дБ/дек. Побудована в такий спосіб пряма лінія збігається з ЛАЧХ при частотах, менших першої зупрягаючої частоти (один по одному їхнього розташування на вісі частот ліворуч-праворуч).

3. На частоті сполучення w_i характеристика змінює свій нахил: або на +20 дБ/дек, якщо постійна часу Т_i=1/--w_i знаходиться в чисельнику вихідної передаточної функції, або на -20 дБ/дек, якщо постійна часу Т_i знаходиться в знаменнику передаточної функції.

Для дослідження системи на стійкість за амплітудно-фазовим критерієм стійкості за допомогою логарифмічних частотних характеристик, необхідно крім ЛАЧХ побудувати ще логарифмічну фазочастотну характеристику (ЛФЧХ). По вісі абсцис відкладається частота в логарифмічному масштабі (використовують ту ж вісь частот, що і для побудови ЛАЧХ), а по вісі ординат відкладають аргумент амплітудно-фазової характеристики j(w)--чи у градусах чи у радіанах в лінійному масштабі. Для розглянутого приклада j(w) розраховується по формулі

j(w)=-----9_-----arctgT_у--w-----arctgT_кзw-----arctgT_дw--------[град.].

2.5.2 Побудова бажаної ЛАЧХ

При побудові бажаної ЛАЧХ виділяють три області: область низьких частот, область середніх частот і область високих частот. Вид ЛАЧХ у кожній з областей по різному впливає на якість системи. В області низьких частот вид ЛАЧХ визначає точність роботи системи в сталих режимах. Область середніх частот визначає динамічні властивості системи (швидкодія, коливальність). Вид ЛАЧХ в області високих частот практично мало впливає на якість системи.

Побудову бажаної ЛАЧХ зручно починати з області середніх частот у такій послідовності.

По заданих величинах s_max і t_р визначаємо за допомогою табл. 1 частоту зрізу w_с.

Наприклад, нехай задано, що --s_max=30 %-- ,----t_р=1,5 с. З табл.1 для s_max =30% --визначаємо t_рw_с/p=3,2, звідки випливає:

w_с= 3,2p/1,5=6,7--c^-1 .

Таблиця 1

?max , %	10	15	20	25	30	35	40
tр?с/?	5	4,4	4	3,6	3,2	3	2,8
L1, дБ	18	15	13,5	12	11	10,5	10
?, гр	85	80	65	55	45	40	35



Рис.3. а-а-а - ЛАЧХ заданої (незмінної частини) системи;

б-б-б - ЛАЧХ бажана;

c-c-c - ЛАЧХ послідовної корегуючої ланки.

Наносимо на вісь абсцис частоту зрізу w_с і проводимо через неї пряму лінію з нахилом -20 дБ/дек (рис.3).

Частота w₂, яка обмежує область середніх частот бажаної ЛАЧХ ліворуч, визначається величиною відрізка L₁, що може бути знайдена в залежності від заданої величини s_max з табл.1. Частота w₃, яка обмежує область середніх частот праворуч, визначається величиною відрізка L₂, при цьому .

В області високих частот бажану ЛАЧХ потрібно будувати у виді прямолінійних відрізків з нахилом, кратним -20 дБ/дек. (тобто -40, - 60, -80 і т.д.), таким чином, щоб різниця характеристик бажаної і заданої (у межі при w--®--Ґ)----складала пряму лінію, рівнобіжну вісі частот.

В області низьких частот бажана ЛАЧХ будується в такий спосіб. По заданій величині коефіцієнта підсилення системи К_с=w_оmах/x_max визначаємо величину L_A2=20LgK_c і відзначаємо на кресленні точку A₂ c координатами w_A2=1 c^-1 і L_A2 (див. рис.3). Через точку A₂ проводимо пряму лінію з нахилом -20 дБ/дек.

Від точки М, що обмежує область середніх частот ліворуч, проводимо пряму лінію з нахилом -40дБ/дек до перетинання з низькочастотною частиною бажаної ЛАЧХ.

Якщо в завданні на розробку слідкуючої системи, зазначена максимально припустима похибка спостереження Хмах за умови, що вхідний сигнал може змінюватися з максимальною кутовою і з максимальним кутовим прискоренням eomax , то для виконання цих вимог необхідно, щоб бажана ЛАЧХ не попадала у заборонену область.

Заборонена область будується в такий спосіб. Визначаємо на кресленні точку В з координатами:

.

Від точки B вправо проводимо пряму лінію з нахилом -40 дБ/дек, а вліво - пряму лінію з нахилом -20 дБ/дек.

Якщо ЛАЧХ, побудована по заданому коефіцієнту Кс, попадає в заборонену область, то це означає, що при даному коефіцієнті Кс задана точність спостереження не може бути забезпечена і коефіцієнт Кс потрібно збільшити, (бажану ЛАЧХ зрушувати вправо), щоб бажана ЛАЧХ не попадала в заборонену область.

По виду ЛАЧХ бажаної можна записати передаточну функцію безперервної зкорегованої (бажаної) системи. Для розглянутого приклада (крива б-б-б... рис.3) передаточна функція має вид:

. (1)

Для визначення передаточної функції бажаної системи можна скористатися програмою, що наведена в додатку 4. Програма написана мовою BASIC і дозволяє знайти постійні часу T1, Т2, Т3 бажаної передаточної функції по показниках якості smах і tp. При цьому передаточна функція записується у виді

. (2)

Типова бажана ЛАЧХ, по якій записана передаточна функція (2), показана на рис. 4.

Показник ступеня k визначається нахилом заданої ЛАЧХ в області високих частот.

Рис 4. Типова ЛАЧХ бажаної слідкуючої системи

Для забезпечення заданих показників якості перехідного процесу скорегована система повинна мати визначений запас стійкості по фазі. Необхідна величина запасу сталості по фазі (q)--для заданої величини максимального перерегулювання зазначена в табл.1. Після побудови бажаної ЛАЧХ потрібно розрахувати і побудувати ЛФЧХ зкорегованої системи і визначити q--.

2.5.3 Розрахунок послідовної корегуючої ланки

ЛАЧХ послідовного безперервної корегуючої ланки будується шляхом графічного вирахування з ЛАЧХ бажаної ЛАЧХ заданої частини системи ( на рис.3 - лінія с-с-с).

По виду ЛАЧХ потрібно записати передаточну функцію безперервної послідовної корегуючої ланки. Для ЛАЧХ с-с-с на рис.3 передаточна функція має вид:-

. (2)

Коефіцієнт К_кор визначається зі співвідношення: 20lgК_кор=L₃ .

Цю ж передаточну функцію можна одержати, якщо бажану передаточну функцію К_Б(Р) поділити на передаточну функцію заданої частини системи К_з(Р).

2.6 Моделювання слідкуючої системи з безперервною послідовною корегуючою ланкою

Для того щоб переконатися, що корекція системи проведена правильно і скорегована система має показники якості перехідного процесу не гірше заданих, потрібно провести моделювання. Рекомендується моделювати скореговану систему на ПЕОМ, використовуючи пакет Matlab simulink.

Якщо виявиться, що показники якості скорегованої системи гірше заданих, то потрібно вносити в корекцію відповідні зміни. Так, якщо перерегулювання виявилося більше заданого, то потрібно збільшувати відрізки L1 і L2. Якщо час регулювання виявилося більше заданого, потрібно збільшувати частоту зрізу бажаної ЛАЧХ.

Література

1. Бесекерский В.А., Попов Е.П. Теория систем автоматического регулирования. М.: Наука, 2001. 766 с.

2. Батоврин А.А., Дашевский П.Г. и др. Цифровые следящие системы судовой автоматики. Л.: Судостроение, 1999. 445 с.

Додаток 1

Варіанти вихідних даних для проектування слідкуючої системи

Вариант	Статичний момент навантаження Мос, Нм	Момент інерції об'єкта управління Jо, кгм2	Максим. кутова швидкість ?о max, с-1	Максим. кутове прискорення ?о max, с-2	Максим. кинетич. похибка Xmax, рад	Максим. пере -регулювання ?max, %
1	2	3	4	5	6	7
1	50	20	0,5	0,02	0,01	20
2	100	70	1,0	0,06	0,02	25
3	200	100	1,5	0,06	0,01	25
4	300	90	2,0	0,08	0,02	30
5	350	120	2,5	0,08	0,02	30
6	50	30	3,0	0,05	0,03	35
7	100	60	0,5	0,02	0,01	20
8	150	90	1,0	0,06	0,02	20
9	200	80	1,5	0,05	0,02	25
10	250	120	2,0	0,08	0,03	25
11	300	100	0,5	0,02	0,01	30
12	350	100	1,0	0,04	0,02	30
13	50	50	0,5	0,05	0,01	35
14	100	100	0,2	0,08	0,02	35
15	150	120	0,5	0,02	0,05	35
16	200	150	1,0	0,03	0,02	20
17	250	180	1,0	0,04	0,01	25
18	300	150	1,5	0,02	0,02	25
19	350	150	0,5	0,02	0,05	30
20	50	70	2,5	0,05	0,01	30
21	100	50	3,0	0,05	0,01	20

1	2	3	4	5	6	7
22	150	100	2,5	0,08	0,05	25
23	200	120	1,5	0,05	0,02	30
24	70	25	0,75	0,025	0,01	20
25	120	50	0,8	0,01	0,02	20
26	110	60	0,5	0,04	0,01	25
27	300	120	0,5	0,01	0,02	25
28	150	100	1,0	0,02	0,05	30
29	100	70	1.2	0,05	0,02	20
30	75	50	1,0	0,01	0,05	25

Учбові групи	Час регулювання, tp, с
АТ - хх1	2
АТ - хх2	2,5
АТ - хх3	3

Додаток 2

Технічні дані двигунів серії МІ

Тип двигуна	Потужність	Швидкість обертання	Напруга живлення	Струм якоря	Опір ланцюга якоря	ККД	Момент інерц.
	Рн, кВт	nн, об/хв	Uн, В	Iн, А	Rд, Ом	?д, %	Jд, кгм2
1	2	3	4	5	6	7	8
МІ-11	0,12	3000	60	2,87	0,46	62	0,0015
	0,1	2000	60	2,27	0,94	63	0,0015
	0,12	3000	100	1,53	1,48	62	0,0015
	0,1	2000	110	1,22	3,60	63	0,0015
МІ-12	0,2	3000	60	4,57	0,23	66	0,002
	0,12	2000	60	2,72	0,52	64	0,002
	0,2	3000	110	2,46	0,765	66	0,002
	0,12	2000	110	1,46	1,74	64	0,002
МІ-21	0,25	3000	60	5,6	0,284	67	0,0035
	0,2	2000	60	4,3	0,645	68	0,0035
	0,25	3000	110	3,05	0,945	67	0,0035
	0,2	2000	110	2,33	2,20	68	0,0035
МІ-22	0,37	3000	60	8,2	0,195	71	0,04
	0,25	2000	60	5,5	0,360	75	0,004
	0,12	1000	60	2,6	1,44	64	0,004
	0,37	3000	110	4,4	0,546	72	0,004
	0,25	2000	110	2,9	1,29	70	0,004
	0,12	1000	110	1,4	4,58	64	0,004
МІ-31	0,45	3000	60	10,3	0,204	68	0,009
	0,37	2000	60	8,2	0,405	70	0,009
	0,45	3000	110	5,6	0,585	68	0,009
	0,2	1000	60	4,4	1,32	66	0,009
	0,37	2000	110	4,4	1,16	70	0,009
	0,2	1000	110	2,4	3,9	66	0,009
МІ-32	0,76	2500	110	8,2	0,368	80	0,0132
	0,45	1500	110	5,0	0,975	75	0,0132
	0,37	1000	110	4,2	2,21	73	0,0132
	0,76	2500	220	4,1	1,36	80	0,0132
	0,37	1000	220	2,1	8,37	73	0,0132
	1,6	2500	110	19,5	0,249	73	0,035
	1,1	1500	110	13,0	0,67	74	0,035
	0,76	1000	110	9,0	1,3	72	0,035
	1,6	2500	220	9,5	0,93	73	0,035
	1,1	1500	220	6,4	2,63	75	0,035
	0,76	1000	220	4,5	5,32	72	0,035
МІ-42	3,2	2500	110	36,3	0,1	78	0,065
	1,6	1500	110	18,2	0,32	78	0,065
	1,1	1000	110	12,6	0,75	75	0,065
	3,2	2500	220	18,0	0,376	79	0,065
	1,6	1500	220	9,1	1,28	78	0,065
	1,1	1600	220	6,3	2,95	75	0,065
МІ-51	5,5	2500	220	27,2	0,164	82	0,125
	3,2	1500	220	17,1	0,46	82	0,125
	1,6	1000	220	8,7	1,1	79	0,125
МІ-52	7,0	2500	220	37,0	0,088	84	0,15
	4,5	1500	220	23,3	0,26	85	0,15
	2,5	1000	220	13,1	0,569	82	0,15

Додаток 3

Технічні дані ЕМП

Тип ЕМП	Потужність ЕМП	Потужність управління	Напруга	Струм якоря	Опір обмотки управління	Постійні часу
						Ту	Ткз
	кВт	Вт	В	А	Ом	с	с
ЕМП-3А3	0,2	0,4	115	1,75	1000	0,005	0,018
ЕМП-5А3	0,5	0,4	115	4,35	1000	0,01	0,033
ЕМП-12А3	1,0	0,4	115	8,7	2200	0,015	0,06
ЕМП-25А3	2,0	0,4	230	9,1	1500	0,02	0,1
ЕМП-50А3	4,0	0,5	230	17,4	2200	0,03	0,17
ЕМП-70А3	6,0	0,5	230	26	1500	0,04	0,22
ЕМП-100А3	8,5	0,5	230	37	1000	0,06	0,28

Програма LOGHAR.

Визначення постійних часу передаточної функції бажаної системи

10 PRINT "Визначення пост. часу передат. функції бажаної системи"

20 PRINT " Перед. ф-ція, яка відшукується, має вид:"

30 PRINT " Kc(T(2)*P+1)/((T(1)*P+1)*(T(3)*P+1)^(n- m)*(T(4)P+1)*..(T(m+3)P+1)*P)"

40 PRINT "n-ступінь полінома Q(P) - знаменника перед. ф-ції заданої сист."

50 PRINT "m-кількість пост. часу Q(P), менших, за Т(3)"

60 PRINT "Q(P)=(Tз(1)*P+1)*(Tз(2)*P+1)*...*(Tз(N)*P+1)"

70 PRINT "Введіть порядок полінома знаменника Q(P) заданої перед.ф-ції"

80 PRINT "N="

90 INPUT N

100 DIM T3(5), T(8)

110 FOR I = 1 TO 8

120 T(I) = 0

130 NEXT I

140 PRINT "Введіть пост. часу знаменника заданої перед. ф-ції Q(P)"

150 FOR I = 1 TO N

160 PRINT "Tз("; I; ")="

170 INPUT T3(I)

180 NEXT I

190 PRINT "Введіть величину макс. перерегулювання"

200 PRINT "Сигма макс. ,% ="

210 INPUT SM

220 IF SM = 10 THEN C = 5: L1 = 18

230 IF SM = 15 THEN C = 4.4: L1 = 15

240 IF SM = 20 THEN C = 4: L1 = 13.5

250 IF SM = 25 THEN C = 3.6: L1 = 12

260 IF SM = 30 THEN C = 3.2: L1 = 11

270 IF SM = 35 THEN C = 3: L1 = 10.5

280 IF SM = 40 THEN C = 2.8: L1 = 10

290 PRINT "Порядок астатизму NU="

300 INPUT NU

310 PRINT "Коефіцієнт підсилення бажаної системи Кс="

320 INPUT KC

330 PRINT "Час регулювання TR="

340 INPUT TR

350 M1 = .434

360 OC = C * 3.14 / TR

370 XC = LOG(OC) * M1

380 B = 20 * XC

390 X2 = (B - L1) / 20

400 T(2) = 1 / (10 ^ X2)

410 A = 20 * LOG(KC) * M1

420 X1 = (L1 + 40 * X2 - A) / 20

430 T(1) = 1 / (10 ^ X1)

440 X3 = (L1 + 20 * XC) / 20

450 T(3) = 1 / 10 ^ X3

460 IF NU = 1 THEN GOTO 490

470 X1 = (40 * X2 + L1 - A) / 40

480 T(1) = 1 / (10 ^ X1)

490 M = 0

510 FOR I = 1 TO N

520 IF T3(I) <= T(3) * .75 THEN M = M + 1: T(3 + M) = T3(I)

540 NEXT I

550 S = N - M

560 IF S = 2 THEN T(3) = .8 * T(3)

570 IF S > 2 THEN T(3) = .7 * T(3)

580 X3 = LOG(1 / T(3))

590 PRINT "T(1)="; T(1); "T(2)="; T(2); "T(3)="; T(3)

600 PRINT "Вираз у дужках знаменника (Т(3)*Р+1) повинен зводитися в ступінь (n-m)"

610 PRINT "n-m="; S

620 FOR I = 1 TO M

630 K = 3 + M

640 PRINT "T("; K; ")="; T(K)

650 NEXT I

660 PRINT "Wc="; OC; "LgWc="; XC

670 PRINT "Lg(1/T1)="; X1; "Lg(1/T2)="; X2; "Lg(1/T3)="; X3

680 END

Размещено на Allbest.ru