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

Теория систем автоматического управления. В.А. Бесекерский, Е.П. Попов, 1975

.pdf
Скачиваний:
6761
Добавлен:
21.03.2016
Размер:
27.19 Mб
Скачать

Термин «коэффициент передачи» можно пояснить следующим образом. Если подать на вход звена только постоянное значение x10 (рис. 3.3, б) и найти установившееся значение выходной величины x30 (рис. 3.3, в), то из (3.9) получимx30 = k1x10 . Таким

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

Следовательно, коэффициент передачи определяет собой наклон (с учетом масштабов по осям) линейной статической характеристики звена (рис. 3.3, а). Заметим, что нелинейную характеристику звена часто называют характеристикой с переменным по входной величине коэффициентом передачи. Из (3.9) очевидно, что

размерностьk1 = размерность выходной ввеличин х3 . размерность входной величиных1

В размерность коэффициента передачи может входить также время t. Так, из уравнения (3.9) следует, что

размерностьk3 = размерность х3 × размерность t размерность х2

а из уравнения (3.10) следует, что для такого звена

 

 

 

 

размерностьk3 =

 

 

 

 

 

размерность х3

 

 

 

 

 

 

 

 

размерность х × размерность t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

Постоянные времени T1 T2 и Т3, как следует из уравнений (3.9) и (3.10),имеют

размерность времени.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вторая форма записи.

Считая условно оператор дифференцирования

p =

d

 

алгебраической

величиной,

решим

уравнение

(3.11)

относительно

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

выходной величины:

 

 

 

(k2 + k3 p)x2 (t)

 

 

 

 

 

x3 =

 

k1x1 (t)

 

+

 

+

 

k4 f1 (t)

 

. (3.13)

1 +T p +T 2 p 2 +T 3 p3

 

1 +T p +T 2 p 2 +T 3 p3

1+T p +T 2 p 2

+T 3 p3

 

 

 

1

2

 

 

 

3

 

 

 

 

1

2

 

 

 

3

1

2

3

 

Выражения

 

 

 

 

 

 

k1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W1 ( p) =

 

 

 

 

 

 

 

,

 

 

 

3.14

 

 

 

 

 

1 +T p +T 2 p2

+T 3 p3

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W2 ( p) =

 

k2

+ k3 p

 

 

 

,

 

 

 

3.15

 

 

 

 

 

1 +T p +T 2 p2

+T

3 p3

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

W3 ( p) =

 

 

 

 

k4

 

 

 

.

 

 

 

3.16

 

 

 

 

 

1 +T p +T 2 p2

+T

3 p3

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

 

3

 

 

 

 

 

 

 

 

 

называются в теории регулирования передаточными функциями. Уравнение (3.13) можно представить в виде

x3 =W1 ( p)x1 (t) +W2 ( p)x2 (t) +W f ( p)f1 (t). (3.17)

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

Передаточные функции, формулы для которых устанавливаются выражениями (3.14)

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

Более строго передаточная функция определяется через изображения Лапласа или Карсона — Хевисайда (см. главу 7). Если ввести изображения, например по Лапласу, входных и выходных величин звена:

X1 (s) = L[x1 (t)], X 2 (s) = L[x2 (t)],

X 3 (s) = L[x3 (t)], F1 (s) = L[f1 (t)],

где s = с + jw — комплексная величина, то передаточную функцию (3.14) можно строго определить как отношение изображений выходной и входной величин звена:

W

(s) =

X 3

(s)

=

 

 

 

k1

 

, (3.18)

X

 

(s)

1

+T s +T 2 s2

+T 3 s3

1

 

1

 

 

 

 

 

 

 

 

 

1

2

3

 

при нулевых начальных условиях и равных нулю остальных воздействиях на звено;

X 2 (s) = 0 , F1 (s) = 0 . Аналогичным образом можно определить передаточные функции (3.15) и (3.16). Поэтому вместо дифференциального уравнения (3.17), куда входят

функции, времени

x1 (t) ,

x2 (t) , x3 (t) ,

можно написать

при

нулевых

начальных

условиях уравнение для изображений в виде, совпадающем по форме с (3.17):

 

 

 

X 3 (s) =W1 (s)X1 (s) +W2 (s)X 2

(s) +W f (s)F1 (s).

(3.19)

 

 

или в развернутом виде:

 

 

(k2 + k3 s)X 2 (t)

 

 

 

 

 

X 3 (s) =

 

k1X1 (s)

 

k4 F1 (s)

 

 

 

 

+

 

 

 

+

 

 

 

. (3.20)

1 +T s +T 2 s2

+T 3 s3

1 +T s

+T 2 s2

+T 3 s3

1 +T s +T 2 s2

+T 3 s3

 

1

2

3

1

2

3

 

 

1

2

3

 

 

Вдвух последних выражениях фигурируют не функции времени, а их изображения: X1 (s), X 2 (s), X 3 (s), и F1 (s), где s = с + jw — комплексная величина.

Визображениях Лапласа и Карсона — Хевисайда комплексная величина часто обозначается той же буквой р, что и оператор дифференцирования, причем р = с + jw. В этом случае уравнение (3.19) будет иметь вид

X 3 ( p) =W1 ( p)X1 ( p) +W2 ( p)X 2 ( p) +W f ( p)F1 ( p). (3.21)

Здесь, как и в уравнении (3.19), фигурируют изображения функций X1 ( p), X 2 ( p),

X 3 ( p), и F1 ( p) .

В дальнейшем

будет употребляться

символ

дифференцирования p =

 

d

 

для

 

dt

 

 

 

 

 

 

 

символической записи дифференциальных

уравнений, куда входят функции

времени

x1 (t) , x2 (t) и т. д., и

комплексная

величина р = с + jw для записи уравнений

с

изображениями функций времени по

Лапласу или

Карсону — Хевисайду

X1 ( p),

X 2 ( p) и т. д. Запись передаточных функций звена, в том и в другом случае сливается в одну: W1 ( p), W2 ( p) и т. д. Однако в передаточных функциях буква р будет означать

символ дифференцирования p = dtd или комплексную величину р = с + jw в зависимости от того, рассматриваются ли функции времени или их изображения.

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

схемы, показанной на рис. 3.4. Передаточные функции звеньев или отдельных участков схемы позволяют легко получить общее уравнение всей системы в виде (3.13) или (3.20), а в дальнейшем в случае необходимости перейти к исходному дифференциальному уравнению вида (3.9). Подобным же образом могут быть получены передаточные функции и структурные схемы и для других дифференциальных уравнений звеньев, например для рассмотренного выше уравнения (3.10). Подробнее этот вопрос изложен в §

5.4.

ГЛАВА 4 ДИНАМИЧЕСКИЕ ЗВЕНЬЯ И ИХ ХАРАКТЕРИСТИКИ

§ 4.1. Общие понятия

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

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

.

Обозначим входную величину звена через х1 а выходную через х2 (рис. 4.1). Возмущение, действующее на звено, в соответствии с изложенным выше обозначим f(t).

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

Взвеньях позиционного, или статического, типа линейной зависимостью х2 = kx1 связаны выходная и входная величины в установившемся режиме (рис. 4.2, а). Коэффициент пропорциональности k между выходной и входной величинами представляет собой коэффициент передачи звена.

Взвеньях интегрирующего типа линейной зависимостью - dxdt2 = kx1 связаны

производная выходной величины и входная величина в установившемся режиме (рис. 4.2,6). В этом случае для установившегося режима будет справедливым равенство

x2 = k x1dt , откуда и произошло название этого типа звеньев. Коэффициент

пропорциональности k в этом случае также является коэффициентом передачи звена. Если входная и выходная величины звена имеют одинаковую размерность, то коэффициенту передачи соответствует размерность [сек-1].

В звеньях дифференцирующего типа линейной зависимостью x2 = k dxdt1 связаны в

установившемся режиме выходная величина и производная входной (рис. 4.2, в), откуда и произошло название этого типа звеньев. Коэффициент пропорциональности k; является коэффициентом передачи звена. Если входная и выходная величины имеют одинаковую размерность, то коэффициенту передачи соответствует размерность [сек].

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

дифференциальным уравнением, представленным в стандартной форме:

T 2

d 2 x

2

+T

dx

2

+ x

 

= k x + k

 

dx

+ k

 

f (t)

 

 

 

 

1

 

 

 

 

 

 

2 dt

 

2 dt 2

 

1 dt

 

2

1 1

 

3

 

.

При нулевых начальных условиях, т. е. в том случае, если для t <0 входная и выходная величины, а также их производные тождественно равны нулю, и при отсутствии внешнего возмущения (f(t)=0) может быть найдена передаточная функция звена как отношение изображений по Лапласу (или Карсону) выходной и входной величин:

W ( p) =

X 2

( p)

=

 

k1

+ k2 p

=

 

 

k1 (1 +T3 p)

(4.1)

X

1

( p)

1 +T p +T 2 p2

1

+T p +T 2 p2

 

 

 

 

 

 

 

 

1

2

 

 

 

1

2

 

где k1 — коэффициент передачи звена, T3 =k2/k1 - постоянная времени.

При известной передаточной функции выходная величина (точнее ее изображение по Лапласу или по Карсону) может находиться из выражения

X 2 ( p) =W ( p)X1 ( p).

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

W

F

( p) =

X 2 ( p)

=

 

k3

 

(4.2)

F( p)

1 +T p +T 2 p2

 

 

 

 

 

 

 

 

 

1

2

 

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

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

Втабл. 4.1 не приводятся сведения о большой группе так называемых корректирующих звеньев, используемых для улучшения динамических качеств автоматических систем. Эти звенья будут рассмотрены в главе 10.

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

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

Переходная функция, или переходная характеристика, h (t) представляет собой переходный процесс на выходе звена, возникающий при подаче на его вход скачкообразного воздействия при величине скачка, равной единице (рис. 4.3). Такое входное воздействие называется ступенчатой единичной функцией и обозначается x1(t)=1(t), что соответствует х1 = 0 при t<0 и x1=1 при t >0. Предполагается, что единица имеет ту же размерность, что и физическая величина на входе звена.

Если входное воздействие представляет собой неединичную ступенчатую функцию х1=N1(t), выходная величина будет равна х2 = N1 (t).

Более строго переходную функцию можно определить как отношение выходной величины звена х2 (t) к высоте ступенчатого скачка х1 (t) = N1 (t) на его входе, т. е. h (t)= N-1x2(t) . При этом размерность h(t) соответствует размерности передаточной функции звена.

Ступенчатая функция представляет собой распространенный вид входного воздействия в автоматических системах. К такому виду сводятся мгновенное изменение нагрузки электрического генератора, мгновенное возрастание нагрузки на валу двигателя, мгновенный поворот входного валика следящей системы и т. п.

Умножение какой-либо функции времени х (t) на ступенчатую единичную функцию 1 (t) означает, что функция времени х (t) будет существовать только при t>= 0, при t < 0 она обращается в нуль. Это иллюстрируется рис. 4.4.

.

Функции веса w(t) представляет собой реакцию звена на единичную импульсную функцию, поданную на его вход рис. 4.5). Единичная импульсная функция, или дельтафункция, представляет собой производную от ступенчатой единичной функции: δ (t) =1'(t) . Дельта-функция тождественно равна нулю повсюду, кроме точки t= 0, где она

стремится к бесконечности.

Основное свойство дельта-функции заключается в том, что

+∞

δ(t)dt =1 (4.3)

−∞

т. е. она имеет единичную площадь.

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

Дельта-функция может быть представлена как предел некоторого выражения, например:

δ(t) = limσ eσt 1(t)

Нетрудно установить связь между переходной функцией и функцией веса. Рассмотрим входное воздействие звена в виде конечного по высоте и ширине импульса с площадью Nε = 1, прикладываемого при t = 0 (рис. 4.6). Такой импульс может быть заменен двумя ступенчатыми функциями N1 (t) и — N1 (t -ε ), прикладываемыми ко входу звена со сдвигом во времени ε . Тогда выходная величина звена будет равна

x2 (t) = N[h(t) h(t ε)] (4.4)

Будем теперь увеличивать высоту импульса N, одновременно уменьшая его ширину ε , но так, чтобы все время площадь импульса равнялась единице,

т. е. Nε = 1. Помножив и поделив правую часть равенства (4.4) на ε и перейдя к пределу, получим. функцию веса

w(t) = lim

Nε[h(t) h(t ε)]

=

dh(t)

(4.5)

ε

 

dt

 

 

 

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

.

В случае, если на вход звена поступает неединичная импульсная функция х1 = Gδ (t), на выходе звена получится х2 = Gω (t).

Более строго функцию веса можно определить как отношение выходной величины звена х2(t) к площади поданного на его вход импульса х1 = Gδ (t), т. е. w(t) =G-1х2 (t). При этом размерность w(t) соответствует размерности передаточной функции звена, деленной на время.

Импульсная функция также представляет собой распространенный вид входного воздействия в автоматических системах. К такому виду можно свести например, кратковременный удар нагрузки на валу двигателя, кратковременный ток короткого замыкания генератора, отключаемый плавкими предохранителями, и т. п. В действительности реальные импульсные воздействия на автоматическую систему всегда будут конечными по величине и продолжительности. Однако в случае, если их продолжительность весьма мала по сравнению с временем переходного процесса звена или автоматической системы, то с большой степенью точности реальный импульс может быть заменен дельта-функцией с некоторым масштабирующим коэффициентом, что позволяет оценить переходный процесс по виду функции веса.

Функция веса звена связана с его передаточной функцией преобразованием Лапласа, а именно: передаточная функция есть изображение функции веса и связана с ней интегральным преобразованием

W ( p) = w(t)ept dt (4.6)

0

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

W ( p) = ph(t)ept dt (4.7)

0

Для входного воздействия произвольного типа, прикладываемого в момент t= 0, переходный процесс на выходе звена при нулевых начальных условиях может быть определен на основании интеграла Дюамеля — Карсона по переходной функции:

x2 (t) = x1 (0)h(t) + x1 (τ)h(t τ)dτ, (4.8)

0

или по функции веса:

x2 (t) = x1 (τ)w(t τ)dτ,

0

еде τ — вспомогательное время интегрирования, изменяющееся в пределах от нуля до рассматриваемого текущего момента времени t.

Более подробно методика нахождения переходного процесса при произвольном входном воздействии будет рассмотрена в главе 7.

§ 4.3. Частотная передаточная функция и частотные характеристики

Важнейшей характеристикой динамического звена является его частотная передаточная функция. Для получения ее рассмотрим динамическое звено (рис. 4.1) в случае, когда возмущение f(t) = 0, а на входе имеется гармоническое воздействие х1 =X1Mсоs wt, где Х1M — амплитуда, а w — угловая частота этого воздействия.

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

x2 = X 2M cos(wt +ψ ).

Воспользуемся формулой Эйлера и представим входную и выходную величины в виде суммы экспоненциальных функций:

x =

 

X1M

[e jwt + e

jwt ] = x '+x ",

 

 

 

 

 

 

1

 

2

 

 

 

 

1 1

 

 

 

 

 

 

 

 

 

 

(4.10)

 

 

 

X 2M

 

 

 

 

 

 

 

x2

=

 

 

[e

j(wt+ψ )

+ e

j(wt +ψ )

] = x2

'+x2

 

 

 

 

 

 

".

2

 

 

 

 

 

 

 

 

 

 

 

 

 

В линейной системе на основании принципа суперпозиции можно рассмотреть

отдельно прохождение составляющих x1 '

и

x1" . Кроме того, можно легко показать, что

достаточно рассмотреть прохождение только составляющейx1 ' , которая в выходной

величине дает составляющуюx2 ' .

Соотношение

между

составляющими x1" и x2 "

получается таким же, как между

x1 ' и x2 ' .

Поэтому в

дальнейшем

рассмотрении

воспользуемся символической записью соs wt = e jwt . Тогда

 

 

 

x1

= X1M e

jwt

,

 

 

 

 

 

 

 

 

(4.11)

 

 

 

 

 

 

 

 

 

 

 

 

x2

= X 2M e

j(wt+ψ )

 

 

 

 

 

 

.

 

 

 

 

Символичность этой сокращенной

записи

заключается

в

отбрасывании

составляющих с множителем e jwt .

 

 

 

 

 

 

 

 

 

Для нахождения соотношения между входной и выходной гармоническими величинами звена воспользуемся его дифференциальным уравнением в виде

T 2

d 2 x

2

+T

dx

2

 

+ x

 

= k x

+ k

 

 

dx

 

 

 

 

 

 

 

1

 

 

 

 

dt

 

 

 

 

dt

2 dt 2

 

 

 

 

1

 

 

2

1 1

 

2

 

Из выражений 4.11 определим производные:

 

 

 

 

 

 

 

 

 

dx1

 

 

= jwX

1M e jwt ,

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx2

 

= jwX

2M e j(wt +ψ ) ,

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

2 x

2

= ( jw)2 X 2M e j(wt +ψ ) .

 

 

 

 

 

 

dt 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставив значения входной и

выходной величин

и

 

их производных в исходное

дифференциальное уравнение, получим:

T22 (jw)2 X 2M e j(wt +ψ ) +T1 jwX 2M e j(wt+ψ ) + X 2M e j(wt+ψ ) = k1 X1M e jwt + k2 jwX1M e jwt ,

откуда после сокращения на общий множитель e jwt найдем:

X 2M

e jψ

=

 

k1 + k2 jw

=W ( jw) (4.12)

X

1M

1+T jw +T 2

( jw)2

 

 

 

1

2

 

 

Это выражение называется частотной передаточной функцией звена» Таким образом, частотная передаточная функция W(jw) представляет собой комплексное число, модуль которого равен отношению амплитуды выходной величины к амплитуде входной, а аргумент — сдвигу фаз выходной величины по отношению к входной:

modW ( jw) = W ( jw) = X 2M ,

X1M (4.13)

argW ( jw) =ϕ.

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

W ( jw) =

X 2 (iw)

=W ( p)

 

p= jw (4.14)

 

X1 ( jw)

 

 

 

 

что непосредственно вытекает из формулы (4.1) при переходе от изображения Лапласа к изображению Фурье; следовательно, частотная передаточная функция легко получается из обычной передаточной функции подстановкой р = jw.

Частотная передаточная функция звена есть изображение Фурье его функции веса, т. е. имеет место интегральное преобразование

W ( jw) = w(t)ejwt dt. (4.15)

0

Частотная передаточная функция может быть представлена в следующем виде:

W ( jw) = A(w)e jψ =U (w) + jV (w), (4.16)

где А (w) — модуль частотной передаточной функции, ψ (w) — аргумент или фаза, U

(w) и V (w) — вещественная и мнимая составляющие частотной; передаточной функции. Модуль частотной передаточной функции находится как отношение модулей числителя и знаменателя. Для рассмотренного выше примера (4.12), аргумент или фаза частотной передаточной функции находится как разность аргументов числителя и

знаменателя. Для (4.12) имеем:

ψ (w) = arctg

k2 w

arctg

 

 

T1w

.

 

1

T 2 w2

 

k

1

 

 

 

 

 

 

 

2

 

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

U (w) =

 

k

(1 T 2 w2 ) + k T w2

,

 

1

2

2 1

 

(1T 2 w

2 ) +T 2 w2

 

 

 

 

 

 

2

1

 

V (w) =

k

w(1 T 2 w2 ) k T w

.

2

2

1 1

 

(1 T 2 w

2 ) +T 2 w2

 

 

 

 

 

 

2

1

 

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

Амплитудно-фазовая частотная характеристика (а. ф. х.) строится на комплексной плоскости. Она представляет собой геометрическое место концов векторов (годограф), соответствующих частотной передаточной функции W ( jw) =U ( jw) + jV (w) при

изменении частоты от нуля до бесконечности (рис. 4.7). По оси абсцисс откладывается вещественная часть U (w) = ReW ( jw) и по оси ординат — мнимая часть V (w) = ImW ( jw) .

Для каждой частоты на комплексной плоскости наносится точка. Полученные точки