Расчет настроек автоматического регулятора
Пермский Государственный Технический Университет
Курсовая работа
по предмету Автоматизация технологических процессов и
производств.
Тема Расчет настроек автоматического регулятора.
Выполнил ст-т гр. АТП-93
Никулина Д. В.
Проверил Бильфельд Н. В.
г. Березники, 1998
Содержание.
1. Координаты кривых разгона.
1.1 Схемы для Ремиконта.
1.2 Координаты и график кривой разгона по возмущению.
1.3 Координаты и график кривой разгона по заданию.
1.4 Координаты и график кривой разгона по управлению.
2. Интерполяция по 3 точкам.
2.1 Линейное сглаживание и график кривой разгона по возмущению.
2.2 Линейное сглаживание и график кривой разгона по заданию.
2.3 Линейное сглаживание и график кривой разгона по управлению.
3. Нормирование кривых разгона.
3.1 Нормирование кривой разгона по возмущению.
3.2 Нормирование кривой разгона по заданию.
3.3 Нормирование кривой разгона по управлению.
4. Аппроксимация методом Симою.
4.1 По возмущению.
4.2 По заданию.
4.3 По управлению.
5. Проверка аппроксимации методом Рунге — Кутта.
5.1 По возмущению.
5.2 По заданию.
5.3 По управлению.
5.4 Сравнение передаточных функций.
5.5 Сравнение кривых разгона.
6. Расчет одноконтурной АСР методом Роточа.
1. Задание
Исследовать работу комбинированной автоматической системы управления в целом и ее отдельных контуров. Провести расчет оптимальных настроечных параметров регуляторов АСР.
2. Координаты кривой разгона
С помощью программы связи ЭВМ с контроллером снимаем координаты кривой разгона.
Для этого сначала поочередно программируем Ремиконт
1. по возмущению
2. по заданию
3. по управлению
В программе тренды меняя задание добиваемся устойчивости систем.
После того как системы установились приступаем к проведению эксперимента. Для этого устанавливаем алгоблок 11 с которого будем снимать кривую разгона, алгоблок 11 на который будем подавать скачек, амплитуду скачка 10 и интервал времени 0,5.
После просмотра полученных точек кривых разгона удаляем одинаковые.
И строим соответствующие графики.
2.1 Координаты и график кривой разгона по каналу возмущения
табл.2.1
1
0,0000
39,0000
22
11,0000
47,1500
2
1,0000
39,1500
23
11,5000
47,4500
3
1,5000
39,3500
24
12,0000
47,6000
4
2,0000
40,0000
25
12,5000
47,7000
5
2,5000
40,4000
26
13,0000
47,8500
6
3,0000
40,8000
27
13,5000
48,0500
7
3,5000
41,2000
28
14,0000
48,1000
8
4,0000
42,0500
29
14,5000
48,2000
9
4,5000
42,5000
30
15,0000
48,2500
10
5,0000
42,9000
31
15,5000
48,4000
11
5,5000
43,3000
32
16,0000
48,4500
12
6,0000
44,0500
33
16,5000
48,5000
13
6,5000
44,4000
34
17,0000
48,5500
14
7,0000
44,7500
35
17,5000
48,6500
15
7,5000
45,1000
36
19,0000
48,7000
16
8,0000
45,6500
37
19,5000
48,7540
17
8,5000
45,9000
38
20,0000
48,8000
18
9,0000
46,1500
39
21,5000
48,8500
19
9,5000
46,4000
40
22,0000
48,9000
20
10,0000
46,8000
41
26,5000
48,9500
21
10,5000
47,0000
42
27,0000
49,0000
2.2 Координаты и график кривой разгона для внутреннего канала
табл.2.2
1
0,0000
58,0000
30
15,0000
65,9500
2
1,0000
58,0500
31
15,5000
66,1000
3
1,5000
58,3000
32
16,0000
66,2000
4
2,0000
58,4500
33
16,5000
66,4000
5
2,5000
58,7000
34
17,0000
66,5000
6
3,0000
59,2000
35
18,5000
66,6000
7
3,5000
59,4500
36
18,0000
66,6500
8
4,0000
59,7000
37
18,5000
66,8000
9
4,5000
60,0000
38
19,0000
66,9000
10
5,0000
60,6000
39
19,5000
66,9500
11
5,5000
60,8500
40
20,0000
67,0500
12
6,0000
61,1500
41
20,5000
67,1500
13
6,5000
61,4500
42
21,0000
67,2000
14
7,0000
62,0000
43
21,5000
67,2500
15
7,5000
62,3000
44
22,0000
67,3000
16
8,0000
62,5500
45
22,5000
67,3500
17
8,5000
62,8000
46
23,0000
67,4000
18
9,0000
63,3000
47
23,5000
67,4500
19
9,5000
63,5500
48
24,0000
67,5000
20
10,0000
63,7500
49
24,5000
67,5500
21
10,5000
64,0000
50
25,0000
67,6000
22
11,0000
64,4000
51
26,0000
67,6500
23
11,5000
64,5500
52
26,5000
67,7000
24
12,0000
64,7500
53
28,0000
67,7500
25
12,5000
64,9500
54
29,0000
67,8000
26
13,0000
65,2500
55
31,5000
67,8500
27
13,5000
65,4000
56
32,0000
67,9000
28
14,0000
65,5500
57
39,0000
67,9500
29
14,5000
65,8500
58
39,5000
68,0000
2.3 Координаты и график кривой разгона основного контура
табл 2.3
1
0,0000
50,9500
30
14,5000
58,4000
2
0,5000
50,9500
31
15,0000
58,5500
3
1,0000
51,0000
32
15,5000
58,7500
4
1,5000
51,0500
33
16,0000
59,0000
5
2,0000
51,1000
34
16,5000
59,1500
6
2,5000
51,2000
35
17,0000
59,2500
7
3,0000
51,5000
36
17,5000
59,3500
8
3,5000
51,6500
37
18,0000
59,6000
9
4,0000
51,8500
38
18,5000
59,6500
10
4,5000
52,0500
39
19,0000
59,7500
11
5,0000
52,5000
40
19,5000
59,8500
12
5,5000
52,8000
41
20,0000
60,0000
13
6,0000
53,0500
42
20,5000
60,1000
14
6,5000
53,6000
43
21,0000
60,1500
15
7,0000
53,9000
44
21,5000
60,2000
16
7,5000
54,1500
45
22,0000
60,3000
17
8,0000
54,4500
46
22,5000
60,3500
18
8,5000
55,0000
47
23,0000
60,4000
19
9,0000
55,3000
48
23,5000
60,5000
20
9,5000
55,5500
49
24,0000
60,5500
21
10,0000
55,8500
50
25,0000
60,6000
22
10,5000
56,3500
51
25,5000
60,6500
23
11,0000
56,6000
52
26,0000
60,7000
24
11,5000
56,8500
53
27,0000
60,7500
25
12,0000
57,0500
54
27,5000
60,8000
26
12,5000
57,5000
55
30,0000
60,8500
27
13,0000
57,7000
56
30,5000
60,9000
28
13,5000
57,9000
57
36,0000
60,9500
29
14,0000
58,0500
58
36,5000
61,0000
3. Интерполяция по трем точкам.
В программе ASR, пользуясь пунктом “интерполировать по 3-м” поочередно считаем кривые разгона и строим соответствующий график.
3.1.Линейное сглаживание и график кривой разгона для внешнего контура
табл. 3.1
1
0,0000
38,9914
22
11,0000
47,2000
2
1,0000
39,1667
23
11,5000
47,4000
3
1,5000
39,5000
24
12,0000
47,5833
4
2,0000
39,9167
25
12,5000
47,7167
5
2,5000
40,4000
26
13,0000
47,8667
6
3,0000
40,8000
27
13,5000
48,0000
7
3,5000
41,3500
28
14,0000
48,1167
8
4,0000
41,9167
29
14,5000
48,1833
9
4,5000
42,4833
30
15,0000
48,2833
10
5,0000
42,9000
31
15,5000
48,3667
11
5,5000
43,4167
32
16,0000
48,4500
12
6,0000
43,9167
33
16,5000
48,5000
13
6,5000
44,4000
34
17,0000
48,5667
14
7,0000
44,7500
35
17,5000
48,6333
15
7,5000
45,1667
36
19,0000
48,7000
16
8,0000
45,5500
37
19,5000
48,7500
17
8,5000
45,9000
38
20,0000
48,8000
18
9,0000
46,1500
39
21,5000
48,8500
19
9,5000
46,4500
40
22,0000
48,9000
20
10,0000
46,7333
41
26,5000
48,9500
21
10,5000
46,9833
42
27,0000
49,0000
3.2. Линейное сглаживание и график кривой разгона для внутреннего контура
табл 3.2
1
0,0000
57,9667
30
15,0000
65,9667
2
1,0000
58,1167
31
15,5000
66,0833
3
1,5000
58,2667
32
16,0000
66,2333
4
2,0000
58,4833
33
16,5000
66,3667
5
2,5000
58,7833
34
17,0000
66,5000
6
3,0000
59,1167
35
18,5000
66,5833
7
3,5000
59,4500
36
18,0000
66,6833
8
4,0000
59,7167
37
18,5000
66,7833
9
4,5000
60,1000
38
19,0000
66,8833
10
5,0000
60,4833
39
19,5000
66,9667
11
5,5000
60,8667
40
20,0000
67,0500
12
6,0000
61,1500
41
20,5000
67,1333
13
6,5000
61,5333
42
21,0000
67,2000
14
7,0000
61,9167
43
21,5000
67,2500
15
7,5000
62,2833
44
22,0000
67,3000
16
8,0000
62,5500
45
22,5000
67,3500
17
8,5000
62,8833
46
23,0000
67,4000
18
9,0000
63,2167
47
23,5000
67,4500
19
9,5000
63,5333
48
24,0000
67,5000
20
10,0000
63,7667
49
24,5000
67,5500
21
10,5000
64,0500
50
25,0000
67,6000
22
11,0000
64,3167
51
26,0000
67,6500
23
11,5000
64,5667
52
26,5000
67,7000
24
12,0000
64,7500
53
28,0000
67,7500
25
12,5000
64,9833
54
29,0000
67,8000
26
13,0000
65,2000
55
31,5000
67,8500
27
13,5000
65,4000
56
32,0000
67,9000
28
14,0000
65,6000
57
39,0000
67,9500
29
14,5000
65,7833
58
39,5000
68,0000
3.3 Линейное сглаживание и график кривой разгона по основному каналу
табл. 3.3
1
0,0000
50,9500
30
14,5000
58,3333
2
0,5000
50,9500
31
15,0000
58,5667
3
1,0000
51,0000
32
15,5000
58,7667
4
1,5000
51,0500
33
16,0000
58,9667
5
2,0000
51,1167
34
16,5000
59,1333
6
2,5000
51,2667
35
17,0000
59,2500
7
3,0000
51,4500
36
17,5000
59,4000
8
3,5000
51,6667
37
18,0000
59,5333
9
4,0000
51,8500
38
18,5000
59,6667
10
4,5000
52,1333
39
19,0000
59,7500
11
5,0000
52,4500
40
19,5000
59,8667
12
5,5000
52,7833
41
20,0000
59,9833
13
6,0000
53,1500
42
20,5000
60,0833
14
6,5000
53,5167
43
21,0000
60,1500
15
7,0000
53,8833
44
21,5000
60,2167
16
7,5000
54,1667
45
22,0000
60,2833
17
8,0000
54,5333
46
22,5000
60,3500
18
8,5000
54,9167
47
23,0000
60,4167
19
9,0000
55,2833
48
23,5000
60,4833
20
9,5000
55,5667
49
24,0000
60,5500
21
10,0000
55,9167
50
25,0000
60,6000
22
10,5000
56,2667
51
25,5000
60,6500
23
11,0000
56,6000
52
26,0000
60,7000
24
11,5000
56,8333
53
27,0000
60,7500
25
12,0000
57,1333
54
27,5000
60,8000
26
12,5000
57,4167
55
30,0000
60,8500
27
13,0000
57,7000
56
30,5000
60,9000
28
13,5000
57,8833
57
36,0000
60,9500
29
14,0000
58,1167
58
36,5000
61,0000
4. Нормирование кривых разгона.
С помощью программы ASR в пункте нормировать последовательно производим нормирование сглаженных кривых и упорядочиваем время начиная с 0,0000, с шагом 0,5 для того чтобы привести полученную динамическую характеристику к единичному виду.
4.1 Нормированная кривая разгона для внешнего контура
табл.4.1
1
0,0000
0,0000
22
10,5000
0,8201
2
0,5000
0,0175
23
11,0000
0,8401
3
1,0000
0,0508
24
11,5000
0,8585
4
1,5000
0,0924
25
12,0000
0,8718
5
2,0000
0,1407
26
12,5000
0,8868
6
2,5000
0,1807
27
13,0000
0,9001
7
3,0000
0,2356
28
13,5000
0,9117
8
3,5000
0,2923
29
14,0000
0,9184
9
4,0000
0,3489
30
14,5000
0,9284
10
4,5000
0,3905
31
15,0000
0,9367
11
5,0000
0,4421
32
15,5000
0,9450
12
5,5000
0,4921
33
16,0000
0,9500
13
6,0000
0,5404
34
16,5000
0,9567
14
6,5000
0,5754
35
17,0000
0,9634
15
7,0000
0,6170
36
18,5000
0,9700
16
7,5000
0,6553
37
18,0000
0,9750
17
8,0000
0,6903
38
18,5000
0,9800
18
8,5000
0,7152
39
19,0000
0,9850
19
9,0000
0,7452
40
19,5000
0,9900
20
9,5000
0,7735
41
20,0000
0,9950
21
10,0000
0,7985
42
20,5000
1,0000
4.2 Нормированная кривая разгона для внутреннего контура
табл.4.2
1
0,0000
0,0000
30
14,5000
0,7973
2
0,5000
0,0150
31
15,0000
0,8090
3
1,0000
0,0299
32
15,5000
0,8239
4
1,5000
0,0515
33
16,0000
0,8372
5
2,0000
0,0814
34
16,5000
0,8505
6
2,5000
0,1146
35
17,0000
0,8588
7
3,0000
0,1478
36
17,5000
0,8688
8
3,5000
0,1744
37
18,0000
0,8787
9
4,0000
0,2126
38
18,5000
0,8887
10
4,5000
0,2508
39
19,0000
0,8970
11
5,0000
0,2890
40
19,5000
0,9053
12
5,5000
0,3173
41
20,0000
0,9136
13
6,0000
0,3555
42
20,5000
0,9203
14
6,5000
0,3937
43
21,0000
0,9252
15
7,0000
0,4302
44
21,5000
0,9302
16
7,5000
0,4568
45
22,0000
0,9352
17
8,0000
0,4900
46
22,5000
0,9402
18
8,5000
0,5233
47
23,0000
0,9452
19
9,0000
0,5548
48
23,5000
0,9502
20
9,5000
0,5781
49
24,0000
0,9551
21
10,0000
0,6063
50
24,5000
0,9601
22
10,5000
0,6329
51
25,0000
0,9651
23
11,0000
0,6578
52
25,5000
0,9701
24
11,5000
0,6761
53
26,0000
0,9751
25
12,0000
0,6993
54
26,5000
0,9801
26
12,5000
0,7209
55
27,0000
0,9850
27
13,0000
0,7409
56
27,5000
0,9900
28
13,5000
0,7608
57
28,0000
0,9950
29
14,0000
0,7791
58
28,5000
1,0000
4.3 Нормированная кривая разгона по основному каналу
табл. 4.3
1
0,0000
0,0000
30
14,5000
0,7579
2
0,5000
0,0050
31
15,0000
0,7779
3
1,0000
0,0100
32
15,5000
0,7977
4
1,5000
0,0166
33
16,0000
0,8143
5
2,0000
0,0315
34
16,5000
0,8259
6
2,5000
0,0498
35
17,0000
0,8408
7
3,0000
0,0713
36
17,5000
0,8541
8
3,5000
0,0896
37
18,0000
0,8673
9
4,0000
0,1177
38
18,5000
0,8756
10
4,5000
0,1493
39
19,0000
0,8872
11
5,0000
0,1824
40
19,5000
0,8988
12
5,5000
0,2189
41
20,0000
0,9088
13
6,0000
0,2554
42
20,5000
0,9154
14
6,5000
0,2919
43
21,0000
0,9221
15
7,0000
0,3201
44
21,5000
0,9287
16
7,5000
0,3566
45
22,0000
0,9353
17
8,0000
0,3947
46
22,5000
0,9420
18
8,5000
0,4312
47
23,0000
0,9486
19
9,0000
0,4594
48
23,5000
0,9552
20
9,5000
0,4942
49
24,0000
0,9602
21
10,0000
0,5290
50
24,5000
0,9652
22
10,5000
0,5622
51
25,5000
0,9701
23
11,0000
0,5857
52
25,5000
0,9751
24
11,5000
0,6153
53
26,0000
0,9801
25
12,0000
0,6434
54
26,5000
0,9851
26
12,5000
0,6716
55
27,0000
0,9900
27
13,0000
0,6899
56
27,5000
0,9950
28
13,5000
0,7131
57
28,0000
1,0000
29
14,0000
0,7347
5. Аппроксимация методом Симою.
С помощью программы ASR в пункту аппроксимации последовательно считаем площади каждой из кривой разгона для последующего получения уравнения передаточной функции.
Для кривой разгона по внешнему контуру для объекта второго порядка получаем следующие данные
Значения площадей
F1= 6.5614
F2= 11.4658
F3= -4.5969
F4= -1.1636
F5= 44.0285
F6= -120.0300
Ограничимся второй площадью. F1
,>
a1 = F1 + b1
a2 = F2 + b2 + b1 F2
a3 = F3 + b3 + b2 F1 + b1 F2
a1 = 6.5614 + b1
a2 = 11.4658 + b1 6.5614
0 = — 4.5969 + b1 11.4658
Решив систему получаем b1 = 0.4
a1 = 6.9614
a2 = 14.0904
Тогда передаточная функция объекта второго порядка по внешнему контуру имеет вид
0.4 s
W(s)=——————————
2
14.0904 s + 6.9614 s + 1
Для кривой разгона по внутреннему контуру для объекта второго порядка получаем следующие данные
Значения площадей
F1= 9.5539
F2= 24.2986
F3= -16.7348
F4= -14.7318
F5= 329.7583
F6= -1179.3989
Для определения передаточной функции решаем систему, так как F3<0.
a1 = 9.5539 + b1
a2 = 24.2986 + 9.5539 b2
0 = -16.7348 + b1 24.2986
Решив систему получаем b1 = 0.6887
a1 = 10.2426
a2 = 30.8783
Тогда передаточная функция объекта второго порядка по внутреннему контуру имеет вид
0.6887 s + 1
W(s) = ——————————
2
30.8783s + 10.2426 s + 1
Для кривой разгона по заданию для объкта третьего порядка с запаздыванием получаем следующие данные
Значения площадей
F1= 10.6679
F2= 38.1160
F3= 30.4228
F4= -46.5445
F5= 168.8606
F6= -33.3020
Так как F3
1
W(s) =——————————-
2
38.1160 s + 10.6679 s + 1
6. Проверка аппроксимации методом Рунге — Кутта.
В программе ASR в пункте передаточная функция задаем полученные передаточные функции. И затем строим графики экспериментальной и аналитической кривых разгона (по полученной передаточной функции).
6.1 Для кривой разгона по внешнему контуру
Устанавливаем для проверки методом Рунге-Кутта конечное время 27c, шаг 0,5с.
6.2 Для кривой разгона по внутреннему контуру
Устанавливаем конечное время 39с, шаг 0,5с.
6.3 Для кривой разгона по основному каналу
При задании передаточной функции учитываем чистое запаздывание 0,08с.
Устанавливаем конечное время 32с, шаг изменения 0,5с.
Получили, что кривые разгона практически одинаковы, следовательно аппроксимация методом Симою сделана верно.
6.4 Сравнение экспериментальных и исходных передаточных функции
объект исходная экспериментальная
передаточная передаточная
функция функция
второго порядка 1 0.6887 s
по возмущению W(s)= —————— W(s)= ——————————
2 2
36 s + 12 s + 1 30,8783 s + 10.2426 s + 1
второго порядка 1 0.4 s
по заданию W(s)= —————————— W(s)= ——————————-
2 2
16,1604 s + 8,04 s + 1 14.0904 s + 6.9614 s + 1
третьего порядка 1 1
с запаздыванием W(s)= ————————————- W(s)= ——————————-
по управлению 3 2 2
91.125 s + 60.75 s + 13.5 s + 1 38.1160 s + 10.6679 s + 1
Анализируя таблицу можно сделать вывод о том, что передаточные функции второго порядка практически одинаковы, а третьего порядка значительно отличаются.
6.5 Сравнение экспериментальных и фактических кривых разгона.
Для исходных передаточных функций с помощью программы ASR, пунктов аппроксимация (создать передаточную функцию и изменить время) получим координаты кривых разгона и сравним их с экспериментальной кривой
— по внешнему контуру
— по внутреннему контуру
— по основному каналу
Полученные значению передаточных функций не значительно отличают от фактических, что говорит о достаточно не большой погрешности между фактическими и экспериментальными данными.
Расчет одноконтурной АСР методом Роточа.
В программе Linreg задаем параметры объекта. Выбираем в качестве регулятора ПИ- регулятор. И рассчитываем его настройки
а) для экспериментальной передаточной функции.
В программе Linreg задаем передаточную функцию объекта второго порядка с запаздыванием. Выбираем ПИ-регулятор и определяем его настройки.
Получаем kp = 1.0796
Tu = 8.0434
В программе SIAM пользуясь следующей схемой для одноконтурной системы
Подаем скачек на сумматор, стоящий после запаздывания и получаем график переходного процесса по заданию
Подаем скачек на сумматор, стоящий перед объектом и получаем график переходного процесса по возмущению
б) для фактической передаточной функции
В программе Linreg задаем передаточную функцию объекта третьего порядка с запаздыванием. Выбираем ПИ-регулятор и определяем его настройки.
Получаем kp = 0.8743
Tu = 8.3924
В программе SIAM пользуясь схемой для одноконтурной системы получаем
— переходный процесс по заданию
Расчет каскадной АСР методом Роточа.
а) для экспериментальной передаточной функции.
Первоначально определим настройки внутреннего регулятора для внутреннего контура с передаточной функцией W1(s).
0.4s + 1
W1(s) = —————————
2
14.0904s + 6.9614s +1
С помощью программы ASR получим АФХ по передаточной функции и определим значения u(m,w), v(m,w), a(m,w), w.
v(m,w)
u(m,u)
a(m,w)
w
kp
Tu
1,0000
0,0000
0,0000
0,0000
0,0000
0,0000
1,0211
-0,0678
1,0234
0,0100
15,0783
0,0109
1,0360
-0,1398
1,0454
0,0200
7,4774
0,0211
1,0439
-0,2151
1,0659
0,0300
4,9709
0,0307
1,0442
-0,2931
1,0845
0,0400
3,7336
0,0395
1,0361
-0,3728
1,1012
0,0500
3,0067
0,0475
1,0194
-0,4531
1,1156
0,0600
2,5367
0,0547
0,9936
-0,5329
1,1275
0,0700
2,2147
0,0609
0,9587
-0,6108
1,1368
0,0800
1,9877
0,0660
0,9147
-0,6857
1,1431
0,0900
1,1826
0,0701
0,8619
-0,7559
1,1464
0,1000
1,1713
4,4754
0,8008
-0,8203
1,1464
0,1100
1,6386
4,5739
0,7323
-0,8775
1,1429
0,1200
1,1584
0,0749
0,6576
-0,9263
1,1360
0,1300
1,5905
0,0737
0,5778
-0,9658
1,1254
0,1400
1,6169
0,0711
0,4945
-0,9953
1,1114
0,1500
1,6842
0,0668
0,4095
-1,0143
1,0938
0,1600
1,8064
0,0609
0,3243
-1,0229
1,0731
0,1700
2,0137
0,0533
0,2407
-1,0214
1,0493
0,1800
2,3750
0,0438
0,1601
-1,0103
1,0229
0,1900
3,0885
0,0324
0,0840
-0,9906
0,9942
0,2000
5,0095
0,0000
0,0134
-0,9635
0,9635
0,2100
26,1125
0,0034
Так как настройки регулятора не могут быть отрицательными то ограничимся 3 квадрантом. И с помощью программы на BASIC рассчитаем оптимальные настройки для ПИ — регулятора методом Стефани по следующим формулам
A^2(m,w) m 1
Tu = ———————— , kp = ———- — ———-
w(m^2+1)* v(m,w) v(m,w) u(m,w)
наибольшее отношение kp/Tu и будет оптимальными настройками.
Получили что kp = 1.712763
Tu = 4.47537
В программе SIAM с помощью схемы для одноконтурной системы без запаздывания получаем переходные процессы по заданию и по возмущению
Сравнивая график кривой разгона по основному каналу и переходный процесс внутреннего контура каскадной системы делаем вывод о том, что за время запаздывания основного контура переходный процесс во внутреннем контуре затухнуть не успевает, следовательно передаточная функция эквивалентного объекта имеет вид
Wоб(s) * Wp1(s)
Wоб(s) = ————————— =
1 + Wоб1(s) * Wp1(s)
1 1
——————————— * (1,7128 + ———- )
2 4,4754s
38,1160s + 10,6679s + 1
————————————————————— =
0,4s + 1 1
1 + ————————— * (1,7128 + ———-)
2 4,4754s
14,0904s + 6,9614s + 1
3 2
107.9987s + 67.4444s + 14.6247s + 1
= —————————————————————————
5 4 3 2
4116.4785s + 3186.9547s + 969.316s + 138.1861s + 15.7294s + 1
Определяем настройки ведущего регулятора. Для ПИ-регулятора получаем
kp = 0.1249
Tu = 5.4148
В программе SIAM с помощью схемы каскадной системы получаем переходный процесс по заданию
С помощью схемы каскадной системы получаем переходный процесс по возмущению
б) для реальной передаточной функции.
Определим настройки внутреннего регулятора для объекта второго порядка с передаточной функцией
1
W1(s) =————————-
2
16,1604s + 8.04s + 1
Получаем следующие настройки регулятора kp = 4.3959
Tu = 6.5957
В программе SIAM пользуясь схемой одноконтурной системы без запаздывания получаем графики переходных процессов по заданию и по возмущению
Сравнивая график кривой разгона по основному каналу и переходный процесс внутреннего контура каскадной системы делаем вывод о том, что за время запаздывания основного контура переходный процесс во внутреннем контуре затухнуть не успевает, следовательно передаточная функция эквивалентного объекта имеет вид
Wоб(s) * Wp1(s)
Wоб(s) = ————————— =
1 + Wоб1(s) * Wp1(s)
1 1
——————————— * (4.3959 + ———- )
3 2 6.5957s
91.125s + 60.75s + 13.5s + 1
————————————————————— =
1 1
1 + ———————— * (4.3959 + ———-)
2 6.5957s
16.1604s + 8.04s + 1
3 2
468.5449s + 249.2673s + 37.0334s + 1
= ———————————————————————————————
6 5 4 3 2
42696.154s + 49705.969s + 25770.6474s + 7229.3112s + 1076.6779s+71.4868s+ 1
Определяем настройки ведущего регулятора. Для ПИ-регулятора получаем
kp = 1.2822
Tu = 6.3952
В программе SIAM с помощью схем для каскадной системы получим переходные процессы по заданию и по возмущению
Расчет комбинированной АСР.
а) для эксперементальной передаточной функции
Расчет компенсирующего устройства
В программе SIAM с помощью смоделированной схемы комбинированной системы без компенсатора получим соответствующий переходный процесс
Определим передаточную функцию фильтра для структурной схемы где выход компенсатора поступает на вход регулятора по формуле
Wов(s)
Wф(s) = ——————— ,
Wоб(s) * Wр(s)
где Wов(s) — передаточная функция канала по возмущению,
Wоб(s) — передаточная функция объекта,
Wp(s) — передаточная функция регулятора
0,6887s + 1
——————————
2
30.8783 s + 10.2426 s + 1
Wф(s) = ———————————————————- =
1 1
——————————- * (1.0796 + ———- )
2 8.0434 s
38.8783 s + 10.6679 s + 1
4 3 2
232.5099 s + 40.1406 s + 98.6173 s + 8.6837 s
= ————————————————————
3 2
268.1379 s + 119.8220 s + 18.9263 s + 1
Настроечные параметры компенсирующего устройства будут оптимальными, если АФХ фильтра равны нулю при нулевой и резонансной частоте.
б) для реальной передаточной функции