Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Lab_rab_3

.pdf
Скачиваний:
12
Добавлен:
17.03.2015
Размер:
677.05 Кб
Скачать

Министерство образования Российской Федерации Уфимский государственный нефтяной технический университет Кафедра автоматизации химико-технологических процессов

Методическое руководство

к лабораторной работе

«Исследование динамических характеристик типовых звеньев АСР»

по курсу

«Теория автоматического управления»

Уфа 2002

В работе описана методика исследования временных характеристик типовых звеньев. Описан порядок проведения работы, требования к оформлению отчета. Приведены необходимые сведения по динамическим характеристикам звеньев.

Составители: Аязян Г.К., канд. техн. наук., профессор; Хоробров В.Р., доцент, канд. техн. наук

Рецензент: Веревкин А.П., д-р техн. наук, профессор

Уфимский государственный нефтяной технический университет, 1987 Электронная редакция: О.В. Кирюшин, 2002

1. Цель работы.

Целью работы является исследование динамических свойств типовых звеньев автоматических систем регулирования (АСР) и определение параметров звеньев по экспериментальным характеристикам. В процессе работы студент должен закрепить теоретический материал по разделам «Динамические характеристики линейных систем» и «Типовые динамические звенья», уметь определять по виду переходной характеристики уравнения звеньев и их параметры.

2. Описание установки.

Экспериментально-исследовательская часть работы проводится на компьютерах с использованием пакетов «ТАУ», «SIAM», «VisSim» или др.

3. Порядок выполнения работы.

Каждый студент проводит исследование трех типовых звеньев. Типы звеньев, численные значения коэффициентов их уравнений и начальные условия задаются преподавателем из табл. 1.

3.1.Получить индивидуальное задание от преподавателя (см. табл. 1, согласно варианту).

3.2.Подготовить ответы на контрольные вопросы. В процессе подготовки студент должен изучить необходимый для проведения работы теоретический материал. Необходимо знать способы описания динамических свойств линейных звеньев. Уметь по передаточной функции или дифференциальному уравнению определять аналитические выражения временных характеристик, освоить методику качественного анализа переходных характеристик по передаточным функциям. Для проверки готовности в разделе 6 приведены контрольные вопросы.

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

3.4.В соответствующем пакете имитационного моделирования (по выбору студента) составить расчетные схемы моделей, произвести расчет при заданных начальных условиях, построить переходные кривые.

3.5.Получить аналитические выражения переходных процессов указанных звеньев при нулевых начальных условиях.

3.6.Оформить отчет и защитить работу.

4. Теоретическая часть.

4.1. Способы описания динамических свойств АСР.

Динамические свойства линейного звена или системы могут быть описаны при помощи дифференциальных уравнений (ДУ), передаточных функций, частотных, временных характеристик и т.д.

 

 

 

 

 

 

 

 

Таблица 1

 

Параметры звеньев АСР

 

 

 

 

Номер звена

Параметры звена

 

Начальные условия

вар.

по приложению 3

K

T

T1

 

T2

y(0)

 

y’(0)

1

2

3

4

5

 

6

7

 

8

1

379

1

1

5

 

2

±2

 

±1

2

469

2

0,9

2

 

3

±1,5

 

±0,5

3

379

3

0,8

6

 

2

±1,5

 

±0,8

4

469

4

3

3

 

4

±1

 

±1

5

379

2,5

1

7

 

3

±2

 

±0,5

6

469

3

4

2

 

2

±1,5

 

±0,8

7

379

3,5

3

8

 

3

±1

 

±1

8

469

2

1,6

3

 

3

±1,4

 

±0,3

9

379

3

1

9

 

4

±1,2

 

±0,5

10

469

2,5

2

4

 

5

3±1,5

 

±0,4

11

379

3

3

10

 

4

2±1

 

±0,35

12

469

2,5

4

5

 

6

±1,2

 

±0,2

13

379

3,5

5

11

 

5

±0,8

 

±0,5

14

469

3

2

4

 

4

±1,3

 

±0,25

15

379

2,5

3

12

 

5

±1

 

±0,6

16

469

4

4

5

 

5

±0,5

 

±0,3

17

379

3

2

13

 

6

±1

 

±1

18

469

2

1

6

 

7

±1,5

 

±0,2

19

379

1

0,8

14

 

6

±2

 

±0,25

20

469

2

0,8

7

 

8

±1,5

 

±0,25

21

379

3

1,5

15

 

7

±1,3

 

±0,4

22

469

4

6

6

 

6

±0,8

 

±0,3

23

379

2,5

5

16

 

7

±1

 

±0,4

24

469

2,6

1,3

7

 

7

±1

 

±0,2

25

379

3

3

17

 

8

±0,9

 

±0,3

26

469

3,5

1,5

10

 

12

±0,6

 

±0,15

27

379

3

1

18

 

8

±1

 

±0,2

28

469

2,5

3

12

 

14

±1,5

 

±0,2

29

379

3,5

4

19

 

9

±0,8

 

±0,3

30

469

3

5

10

 

10

±1

 

±0,2

31

379

3,5

3

20

 

9

±1

 

±0,3

32

469

3

2

12

 

12

±1,5

 

±0,15

33

379

2,5

8

18

 

8

±1

 

±0,2

34

469

4

3

8

 

8

±0,5

 

±0,15

Продолжение табл. 1

1

2

3

4

5

6

7

8

 

 

 

 

 

 

 

 

35

379

3

1,5

16

4

±1

±0,2

36

469

2

0,8

9

9

±2

±0,1

37

379

3,5

0,9

14

6

±2,5

±0,4

38

469

3

2

8

6

±1,5

±0,3

39

379

3,5

4

12

3

±1,5

±0,5

40

469

2,5

2

6

4

±0,8

±0,6

41

379

3

4

10

4

±1

±0,4

42

469

3,5

3

4

6

±1,5

±0,15

43

379

3

2

8

2

±1,5

±0,5

44

469

2,5

1,5

3

5

±1,2

±0,3

45

379

3

4

6

2

±2

±1

46

469

3,5

3

12

10

±1

±0,25

47

379

3

4

4

1

±1,5

±1

48

469

2,5

2,5

12

8

±1,2

±0,35

49

379

3,5

3

3

1

±1,5

±0,7

50

469

3

2

10

8

±1

±0,3

51

379

2,5

2

11

3

±0,7

±0,8

52

469

2

2

10

6

±2

±0,5

53

379

3

1

13

4

±2

±0,5

54

469

2

2

14

10

±1

±0,25

55

379

1

2,5

15

5

±2

±0,3

56

469

2,5

3

14

12

±1,5

±0,1

57

379

3

3

17

7

±1

±0,25

58

469

4

5

16

10

±0,75

±0,2

59

379

3,5

1,5

20

8

±1,2

±0,3

60

469

3

2

15

10

±1,5

±0,2

4.1.1 Дифференциальные уравнения.

Как правило, АСР представляет собой соединение нескольких звеньев, каждое из которых может быть описано дифференциальным уравнением. Таким образом, АСР описывается системой линейных дифференциальных уравнений. Исключая промежуточные переменные, можно получить дифференциальные уравнения n-го порядка:

a n

d n y

+ a n −1

d n −1 y

+ ... + a1

 

d y

+ a 0 y =

dt n

dt n −1

 

dt

 

 

 

 

 

 

 

 

(1)

 

 

d m x

 

 

d m−1 x

 

 

 

 

 

 

= b m

+ b m−1

+ ...

+ b1

 

dx

+ b

0 x

dt m

dt m

 

dt

 

 

 

 

 

 

 

 

 

 

где y(t), x(t) – выходная и входная величины соответственно;

y0 = y(0), y (0i) = y(i) (0) – начальные условия выходной величины и ее про-

изводных.

4.1.2 Передаточная функция.

Динамические свойства линейных звеньев (систем) удобнее описы-

вать передаточной функцией W(s), которая представляет собой отношение изображения по Лапласу выходной величины Y(s) к изображению вход- ной величины X(s) при нулевых начальных условиях, т.е.

W(p) =

L{y(t)}

=

Y(p)

(2)

 

 

 

L{x(t)}

X(p)

где L{f(t)} = f (t)e −pt dt - оператор преобразования Лапласа.

0

Преобразуем по Лапласу обе части уравнения (1) при нулевых начальных условиях, используя теорему о дифференцировании оригинала

 

d

i

 

 

 

 

 

d

j

 

 

 

 

 

y

 

 

 

 

x

 

L a i

 

 

 

= a i p i

Y(p)

и

L b j

 

 

 

= b j p j X(p) ,

 

 

 

 

 

 

 

 

 

dt

i

 

 

 

dt

j

 

 

 

 

 

 

 

 

 

 

вместо уравнения (1) получим:

 

 

 

 

 

 

 

(an pn + … + a1 p + a0) Y(p) = (bm pm + … + b1 p + b0) X(p),

(3)

откуда в соответствии с формулой (2) получим выражение для передаточной функции:

 

 

 

 

 

 

 

 

 

W(p) =

b m p m + ... + b1p + b 0

=

B(p)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

(4)

 

 

 

 

 

 

 

a n p n + ... + a1p + b 0

 

 

 

 

 

 

 

 

 

 

 

 

 

A(p)

 

Если начальные условия ненулевые, то изображения производных

(1) определяются по формулам:

 

 

 

 

 

 

 

 

L{a 0 y} =a0 Y(p),

 

 

 

 

 

 

 

 

 

 

dy

 

 

 

 

 

 

 

 

 

 

 

 

L a

1

 

 

 

=a1 p Y(p) – a1 y0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

L a

2

 

 

 

 

 

 

=a2 p2 Y(p) – a2

p y0

– a2 y’0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

L a n

 

 

 

 

 

 

 

=an pn Y(p) – an

pn-1

y0 – an pn-2 y’0 - … - an y(n-1)0,

 

 

 

 

 

 

 

 

 

 

 

dt

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

L b k

 

 

 

 

 

 

 

=bk pk X(p),

 

k = 1, m .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Преобразуя по Лапаласу обе части уравнения (1) с учетом последних

формул, получено:

(an pn + … a1 p + a0) Y(p) – Pн(p) = (bm pm + … + b1 p + b0) X(p), (5) где Pн(p) = an y0pn-1 (an-1 y0 + an y’0) pn-2 + … + (a2 y0 + a3 y’0 + … + an yo(n-2))p +

+ (a1 yo + a2 y’o +… + an yo(n-1) )

 

 

 

 

 

 

 

 

 

 

 

(6)

 

 

n −1

n

 

 

 

 

 

 

 

 

n −1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(i− k −1)

k

 

= C k p

k

 

 

 

 

 

 

 

 

 

или Рн(р) =

a i y o

 

p

 

 

 

.

 

 

 

 

 

 

 

 

 

k =0 i=k +1

 

 

 

 

 

 

 

 

k =0

 

 

 

 

 

 

 

 

 

 

Для системы 1-го порядка (n = 1):

 

Рн(р) = а1 уо.

 

 

 

 

 

Для системы 2-го порядка (n = 2):

 

Рн(р) = а2 уо р + а1 уо + а2 у’о

Для системы 3-го порядка (n = 3):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рн(р) = а3 уо р2 + (а2 уо + а3 у’о)р + а1 уо + а2 у’о + а3 у”о.

Преобразуя выражение (5), получим операторное уравнение звена

при ненулевых начальных условиях

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y(p) =

b m p m + ... + b1p + b 0

X(p) +

 

 

 

Рн (р)

 

 

(7)

 

 

 

 

 

 

p n

 

 

 

 

 

 

 

 

 

a

n

p n + ... + a

1

p + a

0

 

 

 

a

n

+ ... + a

1

p + a

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

откуда, согласно (4), можно записать операторное уравнение звена

 

 

 

 

 

 

 

 

 

 

Y(p) = W(p) X(p) + Fн(p),

 

 

 

(8)

где Fн(р) =

 

 

 

Рн (р)

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

n

p n

+ ... + a

1

p + a

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Слагаемое Fн(р) в уравнении (8) можно рассматривать как изображение некоторого воздействия fн(t), эквивалентного начальным условиям системы и приведенного к выходу (рис. 1).

 

 

 

 

 

 

 

Fн(р)

 

 

 

Х(р)

 

 

 

 

 

 

 

Y(p)

 

 

 

 

W(p)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если же начальные условия нулевые, то из (8) следует известное

уравнение

 

 

 

 

 

 

 

 

 

 

 

 

 

Y(p) = W(p) X(p).

 

 

(9)

Аналитическое выражение для y(t) можно получить из 7), используя

обратное преобразование Лапласа

 

 

 

 

 

 

 

−1

B(p)

 

 

−1

Pн (p)

 

 

 

y(t) = L

 

 

X(p)

+ L

 

 

= y1 + f

н

(t) ,

(10)

 

 

 

A(p)

 

 

 

 

A(p)

 

 

 

где первое слагаемое определяет выходной сигнал при нулевых начальных условиях, второе – реакцию на воздействие, эквивалентное начальным условиям. Формулы обратного преобразования Лапласа приведены в Приложении 1. Пример определение выходного сигнала звена, описываемого дифференциальным уравнением с начальными условиями, приведен в Приложении 2.

4.1.3 Временные характеристики.

Динамические свойства систем автоматического управления и их звеньев могут быть однозначно определены временными характеристиками: переходной и импульсной (весовой). Для получения указанных характеристик на вход системы (звена) подают определенного вида воздействие x(t) и исследуют реакцию системы (звена) y(t) на это воздействие (рис. 2).

x(t)

САУ,

y(t)

 

 

 

звено

 

 

 

 

Рис. 2

Переходной характеристикой системы (звена) h(t) называется ре-

акция звена на единичное ступенчатое воздействие x(t) = 1 (t) при нулевых начальных условиях.

Изображение единичного ступенчатого воздействия по Лапласу рав-

но L{1(t)} = 1 . Изображение переходной функции находится из (9) при p

Х(р) =

1

:

 

 

 

p

 

 

 

 

 

 

 

 

1

 

(11)

 

 

Н(р) = W(р)

 

.

 

 

 

 

 

 

 

p

 

Замечание. При экспериментальном исследовании динамических характеристик на вход объекта подают скачкообразное воздействие x(t) = A.1(t), где А – амплитуда воздействия. Реакцию объекта на это воздействие иногда называют кривой разгона. Кривая разгона отличается от переходной кривой лишь масштабом.

При x(t) = 1(t) и ненулевых начальных условиях выходной сигнал

y(t), согласно (10) может быть записан следующим образом:

 

y(t) = h(t) + fн(t).

(12)

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

Изображение по Лапласу δ-функции:

 

L{δ(t)} = 1.

(13)

Подставляя в уравнение (9) Х(р) = 1, находим:

 

ω(t) = L-1{W(p)},

(14)

т.е. импульсная характеристика является оригиналом передаточной функции. Известно, что импульсная характеристика равна производной пере-

ходной функции, т.е.

 

ω(t) = h’(t).

(15)

 

4.2. Типовые звенья АСР.

Реальные элементы АСР можно аппроксимировать с достаточной степенью точности отдельными типовыми элементарными звеньями или их соединениями.

Типовыми элементарными звеньями принято называть линейные звенья, динамические свойства которых описываются дифференциальными уравнениями не выше второго порядка с постоянными действительными коэффициентами.

В зависимости от свойств все звенья можно разбить на три группы: статические (пропорциональные), дифференцирующие и интегрирующие.

Передаточная функция статических звеньев имеет вид

W(p) =

b

m

p m + ... + b

1

p + b

0

,

(16)

a

n

p n + ... + a

1

p + b

0

 

 

 

 

 

 

 

 

 

 

 

 

при этом ао 0 и bо 0. Тогда W(0) = K

b 0

 

 

и переходная функция h(t) про-

a 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

порционального звена стремится к некоторому пропорциональному значению (рис. 3)

hуст = lim h(t) = W(0) = K .

t →∞

h(t)

t

Рис. 3 Примерами пропорциональных звеньев являются звенья с переда-

точными функциями

 

 

K

K(T1p +1)

W(p) = K,

W(p) =

 

, W(p) =

 

и т.д.

 

 

 

 

Tp +1

T2 p +1

Передаточная функция дифференцирующего звена имеет вид: W(p) = pr W1(p), (17)

где r 1, а W1(p) соответствует пропорциональному звену. Таким образом, W(p) имеет сомножителем pr. Дифференцирующее звено можно получить из передаточной функции (16), если в последней b0 = 1 (r = 1), b0 = b1 = 0 (r = 2) и т.д.

Из (17) следует, что hуст = W(0) = 0.

Примерный вид переходных кривых дифференцирующих звеньев показан на рис. 4.

h(t)

t

Рис. 4 Примеры дифференцирующих звеньев:

W(p) = p, W(p) =

Kp

,

W(p) =

Kp

 

 

 

Tp + 1

(T1p + 1)(T2 p + 1)

 

 

 

Передаточная функция интегрирующего звена имеет вид:

W(p) = W1 (p) , p r

ит.д.

(18)

где r ≥ 1, а W1(p) соответствует пропорциональному звену. Дифференцирующее звено можно получить из передаточной функции (16), если в послед-

ней а0 = 1 (r = 1), а0 = а1 = 0 (r = 2) и т.д.

Из (17) следует, что hуст = W(0) → ∞.

h(t)

t

Рис. 5 Примеры интегрирующих звеньев:

 

K

 

K

K(T1p + 1)

W(p) =

 

,

W(p) =

 

, W(p) =

 

и т.д.

 

 

 

 

p

 

p(Tp + 1)

p(T2 p + 1)

Передаточные функции типовых звеньев, их переходные и импульсные характеристики приведены в Приложении 3. Располагая переходной характеристикой звена, можно определить его параметры.

5. Требования к оформлению отчета.

Отчет должен быть выполнен аккуратно на листах формата А4 и содержать следующие пункты:

1)задание,

2)вывод аналитических выражений переходных процессов указанных звеньев при нулевых начальных условиях,

3)необходимые расчеты и схемы моделирования заданных звеньев,

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]