Расчет настроек автоматического регулятора (работа 3)

Министерство общего и профессионального образования РФ

Пермский государственный технический университет

Березниковский филиал

Курсовая работа

по предмету: Автоматизация технологических процессов и

производств.

Тема: Расчет настроек автоматического регулятора.

Выполнил: ст-т гр. АТП-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. по заданию

3. по управлению

объект 2-го порядка

01 1.1 21 1.2 21 1.3

313 311 312

1 11 1 11 1 11

Н1=-100 Н1=-100 Н1=-100

Н2=100 Н2=100 Н2=100

к6= к6=1 к6=1

Т1= Т1= Т1= по заданию

Т1= Т1= Т1= по возмущению

Объект 3-го порядка с запаздыванием по управлению

01 1.1 21 1.2 21 1.3 21 1.4 26 1.5

315 311 312 313 314

1 11 1 11 1 11 1 11 1 11

Т1= Т1= Т1= Т1= Т1=

В программе тренды меняя задание добиваемся устойчивости систем.

После того как системы установились приступаем к проведению эксперемента. Для этого устанавливаем алгоблок 11 с которого будем снимать кривую разгона, алгоблок 11 на который будем подавать скачек, амплитуду скачка 10 и интервал времени 0,5.

После просмотра полученных точек кривых разгона удаляем одинаковые.

И строим соответствующие графики.

Координаты и график кривой разгона по возмущению.

1

0,0000

30,0000

2

1,0000

30,2000

3

2,0000

31,5000

4

3,0000

32,6000

5

4,0000

33,7000

6

5,0000

35,2500

7

6,0000

36,1000

8

7,0000

36,8500

9

8,0000

37,4500

10

9,0000

38,2000

11

10,0000

38,5500

12

11,0000

38,8500

13

12,0000

39,2000

14

13,0000

39,4000

15

14,0000

39,5500

16

15,0000

39,6500

17

16,0000

39,7500

18

17,0000

39,8000

19

18,0000

39,8500

20

19,0000

39,9000

21

20,0000

39,9500

22

21,0000

39,9500

23

22,0000

39,9500

24

23,0000

40,0000


Координаты и график кривой разгона по заданию.

1

0,0000

50,0000

2

1,0000

50,1500

3

2,0000

51,0000

4

3,0000

51,8000

5

4,0000

52,6500

6

5,0000

53,5000

7

6,0000

54,7000

8

7,0000

55,4000

9

8,0000

56,1000

10

9,0000

56,9000

11

10,0000

57,4000

12

11,0000

57,8000

13

12,0000

58,1500

14

13,0000

58,6000

15

14,0000

58,8500

16

15,0000

59,0500

17

16,0000

59,2500

18

17,0000

59,4000

19

18,0000

59,5000

20

19,0000

59,6500

21

20,0000

59,7000

22

21,0000

59,7500

23

22,0000

59,8000

24

23,0000

59,8500

25

24,0000

59,9000

26

25,0000

59,9000

27

26,0000

59,9500

28

27,0000

59,9500

29

28,0000

59,9500

30

29,0000

59,9500

31

30,0000

60,0000


Координаты и график кривой разгона по управлению.

1

0,0000

40,0000

21

20,0000

47,4000

2

1,0000

40,0000

22

21,0000

48,0000

3

2,0000

40,0000

23

22,0000

48,4500

4

3,0000

40,0000

24

23,0000

48,8000

5

4,0000

40,0000

25

24,0000

48,8000

6

5,0000

40,0000

26

25,0000

49,1000

7

6,0000

40,0000

27

26,0000

49,3500

8

7,0000

40,0000

28

27,0000

49,5000

9

8,0000

40,0000

29

28,0000

49,5000

10

9,0000

40,2500

30

29,0000

49,5000

11

10,0000

40,9000

31

30,0000

49,7000

12

11,0000

40,9000

32

31,0000

49,8000

13

12,0000

41,8000

33

32,0000

49,8000

14

13,0000

42,8500

34

33,0000

49,8500

15

14,0000

43,9500

35

34,0000

49,9000

16

15,0000

44,9500

36

35,0000

49,9500

17

16,0000

44,9500

37

36,0000

49,9500

18

17,0000

45,9000

38

37,0000

49,9500

19

18,0000

46,7500

39

38,0000

49,9500

20

19,0000

47,4000

40

39,0000

50,0000


Интерполяция по трем точкам.

В программе ASR, пользуясь пунктом “интерполировать по 3-м” поочередно считаем кривые разгона и строим соответствующий график.

Линейное сглаживание и график кривой разгона по возмущению.

1

0,0000

29,8167

2

1,0000

30,5667

3

2,0000

31,4333

4

3,0000

32,6000

5

4,0000

33,8500

6

5,0000

35,0167

7

6,0000

36,0667

8

7,0000

36,8000

9

8,0000

36,5000

10

9,0000

38,0667

11

10,0000

38,5333

12

11,0000

38,8667

13

12,0000

39,1500

14

13,0000

39,3833

15

14,0000

39,5333

16

15,0000

39,6500

17

16,0000

39,7333

18

17,0000

39,8000

19

18,0000

39,8500

20

19,0000

39,9000

21

20,0000

39,9333

22

21,0000

39,9500

23

22,0000

39,9667

24

23,0000

39,9917

Линейное сглаживание и график кривой разгона по заданию.

1

0,0000

49,8833

2

1,0000

50,3833

3

2,0000

50,9833

4

3,0000

51,8167

5

4,0000

52,6500

6

5,0000

53,6167

7

6,0000

54,5333

8

7,0000

55,4000

9

8,0000

56,1333

10

9,0000

56,8000

11

10,0000

57,3667

12

11,0000

57,7833

13

12,0000

58,1833

14

13,0000

58,5333

15

14,0000

58,8333

16

15,0000

59,0500

17

16,0000

59,2333

18

17,0000

59,3833

19

18,0000

59,5167

20

19,0000

59,6167

21

20,0000

59,7000

22

21,0000

59,7500

23

22,0000

59,8000

24

23,0000

59,8500

25

24,0000

59,8833

26

25,0000

59,9167

27

26,0000

59,9333

28

27,0000

59,9500

29

28,0000

59,9667

30

29,0000

59,9833

31

30,0000

59,9833


Линейное сглаживание и график кривой разгона по управлению.

1

0,0000

40,0000

21

20,0000

48,4167

2

1,0000

40,0000

22

21,0000

48,6833

3

2,0000

40,0000

23

22,0000

48,9000

4

3,0000

40,0000

24

23,0000

49,0833

5

4,0000

40,0000

25

24,0000

49,3167

6

5,0000

40,0000

26

25,0000

49,4500

7

6,0000

40,0000

27

26,0000

49,5333

8

7,0000

40,0000

28

27,0000

49,6000

9

8,0000

40,0833

29

28,0000

49,7000

10

9,0000

40,6833

30

29,0000

49,7667

11

10,0000

41,2000

31

30,0000

49,8167

12

11,0000

41,8500

32

31,0000

49,8500

13

12,0000

42,8667

33

32,0000

49,9000

14

13,0000

43,9167

34

33,0000

49,9333

15

14,0000

44,6167

35

34,0000

49,9500

16

15,0000

45,2667

36

35,0000

49,9500

17

16,0000

45,8667

37

36,0000

49,9667

18

17,0000

47,1833

38

37,0000

49,9917

19

18,0000

47,6000

39

38,0000

49,9917

20

19,0000

47,9500

40

39,0000

50,0000


Нормирование кривых разгона.

С помощью программы ASR в пункте нормировать последовательно производим нормирование сглаженных кривых и упорядочиваем время начиная с 0,0000, с шагом 1.0

Нормированная кривая разгона по возмущению.

1

0,0000

0,0000

2

1,0000

0,0737

3

2,0000

0,1589

4

3,0000

0,2735

5

4,0000

0,3964

6

5,0000

0,5111

7

6,0000

0,6143

8

7,0000

0,6863

9

8,0000

0,7551

10

9,0000

0,8108

11

10,0000

0,8567

12

11,0000

0,8894

13

12,0000

0,9173

14

13,0000

0,9402

15

14,0000

0,9550

16

15,0000

0,9664

17

16,0000

0,9746

18

17,0000

0,9812

19

18,0000

0,9861

20

19,0000

0,9910

21

20,0000

0,9943

22

21,0000

0,9959

23

22,0000

0,9975

24

23,0000

1,0000

Нормированная кривая разгона по заданию.

1

0,0000

0,0000

2

1,0000

0,0494

3

2,0000

0,1086

4

3,0000

0,1909

5

4,0000

0,2733

6

5,0000

0,3687

7

6,0000

0,4593

8

7,0000

0,5449

9

8,0000

0,6173

10

9,0000

0,6831

11

10,0000

0,7391

12

11,0000

0,7802

13

12,0000

0,8198

14

13,0000

0,8543

15

14,0000

0,8840

16

15,0000

0,9053

17

16,0000

0,9235

18

17,0000

0,9383

19

18,0000

0,9514

20

19,0000

0,9613

21

20,0000

0,9745

22

21,0000

0,9794

23

22,0000

0,9909

24

23,0000

0,9926

25

24,0000

0,9942

26

25,0000

0,9942

27

26,0000

0,9975

28

27,0000

1,0000

Нормированная кривая разгона по управлению.

1

0,0000

0,0000

21

20,0000

0,7606

2

1,0000

0,0000

22

21,0000

0,7957

3

2,0000

0,0000

23

22,0000

0,8424

4

3,0000

0,0000

24

23,0000

0,8691

5

4,0000

0,0000

25

24,0000

0,8907

6

5,0000

0,0000

26

25,0000

0,9091

7

6,0000

0,0000

27

26,0000

0,9324

8

7,0000

0,0000

28

27,0000

0,9458

9

8,0000

0,0083

29

28,0000

0,9541

10

9,0000

0,0384

30

29,0000

0,9608

11

10,0000

0,0684

31

30,0000

0,9708

12

11,0000

0,1201

32

31,0000

0,9775

13

12,0000

0,1852

33

32,0000

0,9825

14

13,0000

0,2869

34

33,0000

0,9858

15

14,0000

0,3920

35

34,0000

0,9908

16

15,0000

0,4621

36

35,0000

0,9942

17

16,0000

0,5271

37

36,0000

0,9958

18

17,0000

0,5872

38

37,0000

0,9958

19

18,0000

0,6689

39

38,0000

0,9975

20

19,0000

0,7189

40

39,0000

1,0000


Аппроксимация методом Симою.

С помощью программы ASR в пункту аппроксимации последовательно считаем площади каждой из кривой разгона для последующего получения уравнения передаточной функции.

Для кривой разгона по возмущению для объекта второго порядка получаем следующие данные:

Значения коэффициентов:

F1= 6.5614

F2= 11.4658

F3= -4.5969

F4= -1.1636

F5= 44.0285

F6= -120.0300

Ограничимся второй площадью. F1>F2, тогда передаточная функция по возмущению для объекта второго порядка имеет вид:

1

W(s)=---------------------------

2

11,4658s + 6.5614s + 1

Для кривой разгона по заданию для объекта второго порядка получаем следующие данные:

Значения коэффициентов:

F1= 9.5539

F2= 24.2986

F3= -16.7348

F4= -14.7318

F5= 329.7583

F6= -1179.3989

Ограничимся второй площадью , с учетом того что F1>F2. Тогда передаточная функция по управлению для объекта второго порядка имеет вид:

1

W(s)=----------------------------

2

24,2986s + 9.5539s +1

Для кривой разгона по заданию для объекта третьего порядка с запаздыванием получаем следующие данные:

Значения коэффициентов:

F1= 10.6679

F2= 38.1160

F3= 30.4228

F4= -46.5445

F5= 168.8606

F6= -33.3020

Ограничимся третьей площадью и учтем что каждая последующая площадь больше предыдущей. Тогда передаточная функция по заданию для объекта третьего порядка с запаздыванием имеет вид:

1

W(s)=----------------------------------------

3 2

30,4228s + 38.1160s + 10.7769 + 1

Проверка аппроксимации методом Рунге - Кутта.

В программе ASR в пункте передаточная функция задаем полученные передаточные функции. И затем строим графики экспериментальной и аналитической кривых разгона (по полученной передаточной функции).

Для кривой разгона по возмущению.

Устанавливаем для проверки методом Рунге-Кутта конечное время 35,5с, шаг 0,5с.

Для кривой разгона по заданию.

Устанавливаем конечное время 55с, шаг 0,5с.

Для кривой разгона по управлению.

При задании передаточной функции учитываем чистое запаздывание 0,08с.

Устанавливаем конечное время 39с, шаг изменения 0,5с.

Получили, что кривые разгона практически одинаковы, следовательно аппроксимация методом Симою сделана верно.

Для объекта второго порядка по возмущению имеем погрешность метода около 25%, по заданию - около 15%, а для объекта третьего порядка с запаздыванием по управлению - около 5%.

Сравним экспериментальные и исходные передаточные функции:

объект исходная экспериментальная

передаточная передаточная

функция функция

второго порядка 1 1

по возмущению W(s)= -------------------- W(s)= -----------------------------

2 2

0.01s + 0.2s + 1 11.465s + 6.5614s +1


второго порядка 1 1

по заданию W(s)= ----------------------- W(s)= -----------------------------

2 2

0.4489s + 1.34s +1 24.2986s + 9.5539s +1


третьего порядка 1000 1

с запаздыванием W(s)= ------------------------------------- W(s)= -------------------------------------

по управлению 3 2 3 2

4.2188s + 168.75s + 2250s + 1 30.4228s + 38.116s + 10.7769s + 1

Полученные значению передаточных функций отличают на 1000 - 7500, что говорит о достаточно большой погрешности между фактическими и экспериментальными данными.

Для исходных передаточных функций с помощью программы ASR, пунктов аппроксимация (создать передаточную функцию и изменить время) получим координаты кривых разгона и сравним их с экспериментальной кривой:

1. по возмущению:

0,0000

0,0000

0,1000

0,2917

0,2000

0,6094

0,3000

0,8066

0,4000

0,9099

0,5000

0,9596

0,6000

0,9824

0,7000

0,9925

0,8000

0,9968

0,9000

0,9987

1,0000

0,9995

1,1000

0,9998

1,2000

0,9999

1,3000

1,0000

2. по заданию:

1

0,0000

0,0000

17

4,0000

0,9822

2

0,2500

0,0547

18

4,2500

0,9871

3

0,5000

0,1723

19

4,5000

0,9907

4

0,7500

0,3083

20

4,7500

0,9933

5

1,0000

0,4399

21

5,0000

0,9951

6

1,2500

0,5565

22

5,2500

0,9965

7

1,5000

0,6549

23

5,5000

0,9975

8

1,7500

0,7350

24

5,7500

0,9982

9

2,0000

0,7987

25

6,0000

0,9987

10

2,2500

0,8484

26

6,2500

0,9991

11

2,5000

0,8867

27

6,5000

0,9993

12

2,7500

0,9158

28

6,7500

0,9995

13

3,0000

0,9378

29

7,0000

0,9997

14

3,2500

0,9542

30

7,2500

0,9998

15

3,5000

0,9665

31

7,5000

0,9999

16

3,7500

0,9755

32

7,7500

1,0000

3. по управлению:

Сравнивая экспериментальные и фактические кривые разгона видим, что они отличаются очень сильно. Фактическая кривая разгона приходит к 1 на много быстрее, чем экспериментальная.

Расчет АФХ передаточных функций.

1. Объект второго порядка по возмущению:

а) экспериментальная:

1

1,0000

0,0000

34

-0,3751

-0,5372

68

-0,1944

-0,0241

2

1,0211

-0,0678

35

-0,3828

-0,5004

69

-0,1891

-0,0210

3

1,0360

-0,1397

36

-0,3877

-0,4653

70

-0,1839

-0,0182

4

1,0438

-0,2150

37

-0,3903

-0,4320

71

-0,1790

-0,0155

5

1,0440

-0,2927

38

-0,3909

-0,4006

72

-0,1741

-0,0131

6

1,0359

-0,3720

39

-0,3897

-0,3709

73

-0,1695

-0,0108

7

1,0191

-0,4516

40

-0,3871

-0,3431

74

-0,1650

-0,0087

8

0,9935

-0,5304

41

-0,3832

-0,3170

75

-0,1606

-0,0067

9

0,9591

-0,6072

42

-0,3783

-0,2927

76

-0,1564

-0,0049

10

0,9161

-0,6805

43

-0,3725

-0,2699

77

-0,1524

-0,0032

11

0,8649

-0,7492

44

-0,3661

-0,2488

78

-0,1484

-0,0017

12

0,8062

-0,8121

45

-0,3592

-0,2291

79

-0,1446

-0,0003

13

0,7408

-0,8681

46

-0,3518

-0,2108

80

-0,1410

0,0011

14

0,6700

-0,9163

47

-0,3442

-0,1939

81

-0,1374

0,0023

15

0,5948

-0,9560

48

-0,3363

-0,1781

82

-0,1340

0,0034

16

0,5166

-0,9868

49

-0,3283

-0,1636

83

-0,1306

0,0045

17

0,4367

-1,0085

50

-0,3202

-0,1501

84

-0,1274

0,0055

18

0,3565

-1,0211

51

-0,3121

-0,1376

85

-0,1243

0,0064

19

0,2774

-1,0249

52

-0,3040

-0,1260

86

-0,1213

0,0072

20

0,2003

-1,0203

53

-0,2960

-0,1153

87

-0,1184

0,0079

21

0,1265

-1,0081

54

-0,2880

-0,1054

88

-0,1156

0,0086

22

0,0567

-0,9890

55

-0,2802

-0,0962

89

-0,1128

0,0093

23

-0,0083

-0,9640

56

-0,2726

-0,0877

90

-0,1102

0,0099

24

-0,0680

-0,9339

57

-0,2651

-0,0799

91

-0,1076

0,0104

25

-0,1222

-0,8997

58

-0,2577

-0,0726

92

-0,1052

0,0109

26

-0,1708

-0,8924

59

-0,2505

-0,0659

93

-0,1028

0,0114

27

-0,2136

-0,8228

60

-0,2435

-0,0597

94

-0,1004

0,0118

28

-0,2509

-0,7817

61

-0,0237

-0,0540

95

-0,0982

0,0122

29

-0,2829

-0,7399

62

-0,2301

-0,0487

96

-0,0960

0,0126

30

-0,3098

-0,6978

63

-0,2237

-0,0437

97

-0,0939

0,0129

31

-0,3322

-0,6562

64

-0,2175

-0,0392

98

-0,0918

0,0132

32

-0,3502

-0,6153

65

-0,2114

-0,0350

99

-0,0898

0,0134

33

-0,3644

-0,5756

66

-0,2056

-0,0311

100

-0,0879

0,0134

67

-0,1999

-0,0275

б) фактическая

2. Объект второго порядка по заданию.

а) экспериментальная

б) фактическая

3. Объект третьего порядка с запаздыванием по управлению

а) экспериментальная

б) фактическая

Расчет одноконтурной АСР методом Роточа.

В программе Linreg задаем параметры объекта с учетом оператора Лапласса.