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

14_08_18_ТАУ_1,2_Лекционный курс

.pdf
Скачиваний:
38
Добавлен:
12.02.2016
Размер:
5.42 Mб
Скачать

K1 < K2. Пусть первая САУ устойчива в замкнутом состоянии, вторая - неустойчива (рис. 8.6).

Im

RE

ωπ

ω

 

π

ωπ

Рисунок 8.6 - Устойчивость двух систем

Если W1( P) - передаточная функция первой САУ, то передаточная функция второй САУ W2 ( P) = KW1( P) . Вторую САУ можно представить последовательной цепочкой из двух звеньев с передаточными функциями K

(усилительное звено) и W1( P) , поэтому результирующие ЛАЧХ строятся как сумма ЛАЧХ каждого из звеньев.

Поэтому ЛАЧХ второй САУ: L2 (ω) = 20lg K + L1(ω) , а ЛФЧХ:

ϕ2 (ω ) = ϕ1(ω ) .

Пересечениям АФХ вещественной оси соответствует значение фазы

ϕ = −π . Это соответствует точке пересечения ЛФЧХ линии −π координатной

сетки. При

этом,

как видно на АФХ (рис. 8.6, а),

амплитуды

A1(ω) < 1, и

A2 (ω) > 1,

что

соответствует на ЛАЧХ

(рис. 8.6, б)

значениям

L1(ω) = 20lg A1(ω) < 0 и L2 (ω) = 20lg A2 (ω) > 0 .

Сравнивая ЛАЧХ и ЛФЧХ можно заключить, что система в замкнутом

состоянии будет устойчива, если значению ЛФЧХ ϕ = −π будут

соответствовать отрицательные значения ЛАЧХ и наоборот. Запасам устойчивости по модулю H1 и H2 , определенным по АФХ соответствуют расстояния от оси абсцисс до ЛАЧХ в точках, где ϕ = −π , но в логарифмическом масштабе.

81

Особыми точками являются точки пересечения АФХ с единичной

окружностью и осью вещественных чисел. Частоты ωC1 и ωC2, при которых АФХ пересекает единичную окружность, называют частотами среза. Частоту

ωπ при которой АФХ пересекает ось вещественных чисел, называют

частотой фазового сдвига.

(

 

)

(

 

)

 

ЛАЧХ пересекает горизонтальную ось в точках A

ω

 

= 1 L

ω

 

= 0 .

Если при частоте среза ωC1 фаза АФХ ϕC1 > −π (рис. 8.6, а кривая 1), то замкнутая САУ устойчива. На рис. 8.6, б это выглядит так, что пересечению ЛАЧХ горизонтальной оси соответствует точка ЛФЧХ, расположенная выше

линии ϕ = −π . И, наоборот, для неустойчивой замкнутой САУ (рис. 8.6 а кривая 2) ϕC2 < −π , поэтому при ω = ωC2 ЛФЧХ проходит ниже линии

ϕ = −π . Угол ϕ1 = ϕC1 ( −π ) является запасом устойчивости по фазе. Этот

угол соответствует расстоянию от линии ϕ = −π до ЛФЧХ.

Исходя из сказанного, критерий устойчивости Найквиста по логарифмическим частотным характеристикам, в случаях, когда АФХ только один раз пересекает отрезок вещественной оси [−∞; −1] , можно сформулировать следующим образом. Чтобы замкнутая САУ была устойчива необходимо и достаточно, чтобы частота фазового сдвига, при которой ЛФЧХ пересекает линию −π , была больше частоты среза ωC .

Рисунок 8.7 - Устойчивость системы со сложным видом АФХ Если АФХ разомкнутой САУ имеет сложный вид (рис. 8.7), то ЛФЧХ

может несколько раз пересекать линию ϕ = −π . В этом случае применение критерия Найквиста несколько усложняется. Однако во многих случаях данной формулировки критерия Найквиста оказывается достаточно.

82

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

Рассмотрим в качестве примера замкнутую систему регулирования, изображенную на рис. 8.8.

X(S)

W1(S)

 

W2(S)

Y(S)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

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

W1( P) =

2

и W2

( P) =

 

1

.

 

 

 

 

4S2

 

 

0.1S +1

 

+ 0.5S +1

Для определения запаса устойчивости по амплитуде и фазе воспользуемся пакетом Matlab. В командной строке пакета зададим две передаточные функции:

H1=TF([2],[0.1 1]) ,

H2=TF([1],[4 0.5 1]) .

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

h=h1×h2.

Для построения ЛАЧХ и ЛФЧХ воспользуемся функцией margin (h ) . Полученные характеристики изображены на рис. 8.9.

83

L(ω )

ωC H

ϕ (ω )

ϕ ωπ

Рисунок 8.9 - ЛАЧХ и ЛФЧК системы, изображенной на рис. 8.9

Вданном примере запас по модулю составляет 8.09 дБ и запас по фазе –

7.53градуса.

9.Расчет регуляторов в системах подчиненного регулирования

9.1.Общие сведения

Рассмотрим расчет регуляторов в системе подчиненного регулирования на примере повода постоянного тока независимого возбуждения.

Для поводов постоянного тока, как правило, применяется многоконтурная система регулирования последовательного действия (рис. 9.1).

X3

X2

X1

Y1

Y2

Y3

Wр3

 

Wр2

Wр1 Wо1

Wо2

Wо3

Рисунок 9.1 - Обобщенная структурная схема системы подчиненного

регулирования

 

 

 

 

 

84

WPI (S) , который

Каждый контур системы имеет отдельный регулятор

настраивается соответственно передаточной функции объекта регулирования

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

В приводах с двигателями постоянного тока независимого возбуждения первый (внутренний) контур предназначен для регулирования тока якоря; второй (внешний по отношению к контуру тока) - для регулирования частоты обращения; третий - для регулирования положения и т.д.

Регулятор внешнего контура W(s) управляет основным параметром (положением). Работа регулятора WP2 (S) , также как и WP1(S) , подчинена той

же главной цели - регулированию основного параметра.

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

звена с передаточной функцией WOI _ K (S) , которая содержит большую

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

W

(s) =W

( s) ×

1

,

m

 

оi

oi_ k

 

(9.1)

 

 

 

(τnis +1)

 

 

 

 

 

 

 

n=1

 

где WOI _ K (S) - передаточная функция звеньев, действие которых устраняется

регулятором.

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

9.2. Настройка контура регулирования на модульный оптимум

Рассмотрим более подробно вопрос компенсации постоянных времени.

85

Представим объект регулирования в виде апериодического звена,

переходная функция которого изображена на рис. 9.2:

 

 

 

Wо

( р) =

 

 

 

1

 

,

 

 

 

 

Тоs +1

 

 

 

(9.2)

 

 

 

 

 

 

 

 

 

 

 

 

 

где TO - постоянная времени объекта регулирования.

 

 

 

Y

 

 

 

 

 

 

 

 

 

x

 

1

Y

 

 

 

 

 

 

 

 

 

WP ( S)

х

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1.0

 

 

 

 

 

 

 

 

 

 

TOS +1

A)

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

TO S + 1

x

 

1

Y

 

 

 

 

 

 

 

WP ( S)

0.5

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

T S +1

 

 

 

 

 

 

 

 

 

 

 

 

O

 

б)

 

 

Tи s + 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рисунок 9.3 -

Структурная

схема

0.0

 

 

 

 

 

 

 

 

T

апериодического

звена

охваченного

 

 

 

 

 

 

 

 

 

Рисунок 9.2 - Переходные функции

отрицательной обратной связью

апериодических звеньев

 

 

 

 

Включим последовательно с объектом регулирования WO(S) регулятор

Wp (s)

- рис. 9.3, а. Очевидно, что для получения идеального переходного

процесса передаточная функция контура регулирования должна иметь вид:

W (S) = W

P

(S)W (S) = 1.

 

 

 

(9.3)

 

 

 

 

 

 

O

 

 

 

 

Для обеспечения условия (9.3) регулятор должен быть выполнен в виде

форсирующего звена с передаточной функцией:

 

 

 

W

р

(s) =

Т

s +1.

 

 

 

 

(9.4)

 

 

о

 

 

 

 

 

 

 

 

Звенья такого типа физически нереализуемы из-за ограничений на

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

амплитудой, которая сопоставима с амплитудой полезного сигнала.

 

Поэтому в реальных системах для объектов, которые имеют

передаточную функцию в виде апериодического звена, применяют ПИ-

регулятор:

 

 

 

 

 

 

 

 

 

 

 

 

 

( s) =

Т

 

s

+1

 

 

 

 

 

 

W

р

 

о

 

 

,

 

 

 

 

(9.5)

 

 

 

Ти s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где Tи - постоянная времени интегрирования регулятора.

86

i − го

Передаточная функция разомкнутого контура (рис. 9.3, б), в котором есть объект с передаточной функцией (9.2) и регулятор с передаточной функцией (9.5) имеет вид:

p

( s) = Wр ( s)W о( s) =

Tos +1

 

Tos +1

 

1

.

(9.6)

W

×

=

 

Tиs

 

 

Tиs

Tos +1

 

Передаточная функция замкнутого контура (рис. 9.3, б):

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

W р(s)

 

 

 

 

 

 

 

 

 

 

 

W з(s) =

 

 

Т

s

 

 

1

 

 

 

 

=

 

 

и

 

 

 

=

 

.

(9.7)

 

р(s)

 

 

 

1

 

Ти s+1

1 +W

1 +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ти s

Таким образом, включение ПИ-регулятора с передаточной функцией (9.5) последовательно с объектом управления вида (9.2) позволяет получить передаточную функцию замкнутого контура, который описывается уравнением (9.7). Переходная функция контура изображена на рис. 9.2. Как следствие, разомкнутый контур с постоянной времени TO заменен замкнутым контуром с

меньшей постоянной времени Tи .

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

р

(

s

)

=

 

1

 

.

 

 

 

 

 

 

 

 

 

Wi

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т

s ∏ (τ

 

s +1)

 

 

(9.8)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ui

n=1

ni

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При τ

NI

<<1 можно приближенно умножение заменить суммой:

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

m

 

 

 

 

 

 

(τnis+1) ≈ s∑ τni

+1 = Т

µi

s +1,

(9.9)

n=1

 

 

 

 

 

 

n=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

где Тµ i = ∑ τni

- сумма малых постоянных времени (некомпенсированная

 

 

n=1

 

 

 

 

 

 

 

постоянная времени).

С учетом (9.9) передаточные функции разомкнутого и замкнутого контуров запишутся в виде:

87

р

 

 

 

 

 

 

 

1

 

 

,

 

Wi (

s

)

=

 

 

 

 

 

 

 

 

 

Т

 

s

(Т

s +1)

(9.10)

 

 

 

 

 

 

 

 

 

ui

 

µi

 

 

 

 

з

 

 

 

 

 

 

 

1

 

 

.

Wi ( s) =

 

 

 

 

 

 

 

 

Т

s

(Т

1

+

1

(9.11)

 

 

 

 

s + )

 

 

 

 

 

 

 

ui

 

 

µi

 

 

 

 

Уравнение (9.11) описывает колебательное звено (рис. 9.4).

Y

3

 

δ

1.0

 

 

 

 

 

2

 

 

0.5

1

 

 

 

 

 

0.0

tн tм1

tр

T

 

 

 

Рисунок 9.4 - Переходные функции колебательного звена Характер переходных процессов, которые протекают в колебательном

звене, определяется значениями TµI и Tиi . Обозначим TUI =AITµI . Тогда,

изменяя аi можно получить необходимое качество регулирования (рис. 9.4). При аi>4- процесс, близкий к апериодическому (рис. 9.4, кривая 1); при

аi = 2- колебательный процесс с перерегулированием 4.3% (рис. 9.4, кривая 2); при аi = 1- колебательный процесс с перерегулированием 15% (рис. 9.4, кривая 3).

В электроприводе оптимальным считается процесс, изображенный на рис. 9.4, кривая 2 (настройка на модульный или технический оптимум):

перерегулирование δ = 4.3%,

время первого достижения установившегося значения tн = 4.7tµ ,

• время регулирования t

р

= 8.4T .

 

µ

При настройке контуров регулирования на модульный оптимум передаточная функция 1-го замкнутого контура примет вид ( A1 = 2):

W

з

(

s

)

 

1

 

.

 

=

 

 

 

 

2Tµs(Tµs +1) +1

(9.12)

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

88

С определенным допущением передаточная функция 2-го замкнутого

контура (а1= 2,

а2 = 2 ):

 

 

W2з( s) =

 

1

 

.

 

 

 

 

(9.13)

 

4Tµs(2Tµs +1) +1

 

 

 

 

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

Определим соотношение между а и ξ для передаточной функции колебательного звена (9.12):

W кол( p) =

 

 

 

1

 

=

 

 

 

 

1

 

 

 

,

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

aT

 

p

+ aT

p +1

 

1

 

2

+

s

+

1

(9.14)

 

 

 

 

2 s

 

 

 

µ

 

µ

 

 

 

 

 

ωп

 

 

 

 

 

 

 

 

 

 

ωп

 

 

 

 

 

 

где ξ - коэффициент демпфирования колебаний; ωN - угловая частота собственных колебаний, имеющих место при ξ = 0.

С возрастанием коэффициента демпфирования частота, при которой

наступает резонанс, уменьшается ω р =ωп 1−ξ 2 . При ξ = 1 колебаний нет –

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

1

 

 

 

1

 

 

,

ω п =

 

 

 

 

=

 

 

 

 

 

 

 

 

(9-15)

 

2

 

 

 

аТ

 

 

Тµ а

 

 

 

 

 

 

 

 

 

 

µ

 

 

 

 

 

 

ξ =

 

aTµωn

 

=

 

а Тµ

 

=

 

а

 

.

(9-16)

 

 

2

 

 

2 Т

 

 

2

 

 

 

 

 

а

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

µ

 

 

 

 

 

 

 

Период собственных колебаний:

 

Т

 

=

=

 

 

 

 

.

 

 

 

 

 

 

 

п

 

 

 

 

 

 

 

 

 

 

 

 

(9-17)

 

 

 

 

 

 

 

 

 

 

 

 

 

ωп

Т

 

 

а

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

µ

 

 

 

 

 

 

 

 

 

 

 

 

89

9.3.Особенности настройки контуров регулирования

9.3.1.Объект регулирования в виде колебательного звена

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

W о( s) =

 

1

 

 

.

 

 

 

 

 

(9.18)

а Т

s(Т s+

) +

1

1 µ

µ

1

 

 

 

 

 

 

Тогда структурная схема второго контура регулирования примет следующий вид:

х

з2

 

 

х

з1

1

х

2

 

WP2

(S)

 

 

 

 

 

 

 

A1TµS (TµS +1) +1

 

 

 

 

 

 

 

 

 

 

Рисунок 9.5 - Контур регулирования с объектом регулирования в виде колебательного звена

Рассмотрим правомерность замены (для расчета регулятора) данного

1

колебательного звена апериодическим звеном

 

, т.е. отбрасыванием

 

а1Тµ s+1

а1Тµ2s2 .

Условие настройки второго контура на модульный оптимум:

 

 

 

1

 

 

 

1

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

Wр2

(s)

 

=

а а Т

s (а Т

s+1)

(9.19)

 

 

 

 

а Т

s+1

 

 

 

1

µ

 

 

2 1 µ

1

µ

 

 

Передаточная функция регулятора второго контура:

1

Wр2(s) = . (9.20)

а1а2Тµ s

При передаточной функции регулятора вида (9.20) второй контур будет настроен на модульный оптимум при а2 = 2.

Если 1≤ξ1

≤∞

(а1≥4) настройки второго контура мало отличаются от

оптимальных. При

ξ1

 

(

1

)

колебания растут. Для

ξ1

 

(

1

)

 

<1

 

а <4

 

 

= 0.707

а

=2

(т.е. в случае, когда объектом регулирования является контур, оптимизированный по модульному оптимуму) δ = 8%. При наличии в объекте

90