Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Проектирование спецмашин Часть 4. Стабилизация машин.pdf
Скачиваний:
90
Добавлен:
15.11.2022
Размер:
12.95 Mб
Скачать

для системы пятого порядка - аАаъ> а5а2;

(а4а3- а 5а2)(а2аЛ- а 3а0) - (а лаха5а0)2 >0 .

Следует заметить, что с ростом порядка системы трудоемкость использования критерия Гурвица возрастает, поэтому практически при п > 5 прибегают к более простым и менее трудоемким крите­ риям, например, алгебраическому критерию Рауса или частотным критериям Найквиста и Михайлова.

1.6.3. Критерий Рауса

При исследовании устойчивости по критерию Рауса также должно выполняться необходимое (но недостаточное) условие ус­ тойчивости: все коэффициенты а.х характеристического уравнения должны быть положительны. Для окончательной проверки устой­ чивости составляется таблица-схема (табл. 1.3) из коэффициентов характеристического уравнения.

Т абли ца 1.3. Схема Рауса

 

а п

 

 

Я„-2

 

а п- 1

 

 

я» -з

I, _ а п-\а п-2 ~ а па п-3

U а п-\а п-\ ~ а па п-5

°\ ~

°п-\

 

1/-)

--

 

 

 

°п-\

сГ- 1

' S 1

сГ-

С

 

 

ь,

 

 

1

 

 

 

^

с хЬг - Ь

хс 2

 

j

“ 1 “

с ,

 

 

 

 

 

 

C I

а п-А

°п-Ъ

К _ а п-\а п-Ь ~ а па п-1

Оу -

 

 

°п-\

<3= С>-

<3= 1

 

К

С,

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

меняются нулями. Всего таблица содержит п +1 строку. Согласно критерию Рауса, для устойчивости системы необходимо и доста­ точно, чтобы при а( >0 все коэффициенты правого столбца табли­ цы были положительны:

ап > 0; ап_х> 0; Ъх>0; сх>0; dx> 0...

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

” \ I* )

9

 

р(Тхр + 1)(Т2р + 1)

где К - коэффициент передачи системы, К > 0; Т\ и Т2- постоянные времени, причем Тх= 0,05 с ; Т2 = 0,5 с .

Р е ш е н и е . Характеристическое уравнение системы по из­ вестной передаточной функции разомкнутой системы определится

 

ТхТ2р ъ+ (Тх+ 7,) р “+ р + К —0

или

я3/?3 + а~,р2+ ахр + и0 = 0 ,

где аъ - ТХ2; а2 = Тх+ Т2; ах= 1; а0 = К

Видно, что все коэффициенты характеристического уравнения ах> 0, то есть необходимое условие выполняется. Достаточное ус­ ловие по Гурвицу определяется выполнением неравенства для сис­ темы третьего порядка

а2ах>а3а0, то есть (Г, Ч-7^)• 1 > Т^З^ЛГ или АГ< — + — .

Т

т

11

12

С помощью последнего неравенства можно построить границу устойчивости системы по двум параметрам К и Т\ при Т2 = 0,5 с (рис. 1.49, а) и К и Т2 при Т\ = 0,5 с (рис. 1.49, б). На этих рисунках области устойчивой работы системы заштрихованы.

Таблица 1.4. Расчетная схема-таблица

аъ -

ТХТ2 = 0,05 • 0,5 = 0,025

 

°\ = 1

 

 

а2 = ТХ+Т 2 =0,55

 

о

*

О

 

 

 

 

 

 

О II

II

 

 

 

0,55 1-0,025

10

0,3

L

0,55 0-0,025 0

Л

1~

0,55

 

~ 0,55

u-i

0,55

 

—U

 

 

 

 

 

0 3

 

 

 

0 3

 

 

 

-•10 -0,55 -0

_

3

 

— •0-0,55-0

 

_ 0,55

с - 0’5-5

 

-

-0

С|

0,55

 

0,55 0,55

0,55

 

с, °-Ь,

° = 0

 

A = < v ° - V °=0

 

0,55

 

 

 

 

0,55

 

 

0

0

0

0

В табл. 1.4 представлена расчетная таблица-схема для провер­ ки на устойчивость (по критерию Рауса) рассматриваемой системы при следующих параметрах: Тх=0,05 с ;Т2=0,5 с ;К = 10. Из нее видно, что все коэффициенты первого столбца таблицы положи­ тельны: о3 > 0 ; а2 >0; Ьх>0; с, > 0. При взятых параметрах рас­ сматриваемая система устойчива.

Таблица 1.5. Расчетная схема-таблица

 

а3 =0,025

 

 

 

а2 =0,55

 

 

 

0,55 1-0,025 100

 

1,95

1_

0,55

~

0,55

 

195

 

 

с _

— •100-0,55-0

_

195

0,55

С'

0,55

 

0,55

 

О

 

 

 

II •хГ

 

 

 

а, =1

Оо

 

и

II

 

fj3- II о

 

с2=0

 

d 2 = 0

0

0

В табл. 1.5 представлена расчетная таблица-схема при Г, = 0,05 с ; Т2 = 0,5 с ; К = 100. В этом случае система неустойчи­

ва, так как о3 > 0 ; а2 > 0, а Ьх и с, отрицательны.

Пример 2. Построить область устойчивости следящей системы наведения установки по азимуту в плоскости двух параметров: ко­

эффициента усиления первого каскада усиления К у и коэффици­ ента гибкой обратной связи Кос.

Функциональная схема следящей системы представлена на рис. 1.50. Сельсинная пара (СП) формирует сигнал главной жест­ кой отрицательной обратной связи (ГЖООС) U ^ it) путем элек­

трического соединения двух сельсинов, работающих в трансфор­ маторном режиме, из которых: сельсин-датчик размещен в пульте управления, а его ротор оператором разворачивается на требуемый азимут Увых(0 ; сельсин-приемник размещен на установке (объек­

те), а его ротор поворачивается только вместе с объектом на теку­

щий азимут \|/

(г).

 

 

 

</„(')

Ис„(0

">(') и-(0

«г(0

С0д(г)

о)д(0 = о уя(0

-----СП |—

Усилитель

.. /Т\Г

Генератор|----►

|Элсктродвшитс.1ь|— |—»-|

Редуктор |-|---------- 1

г Т

------

 

-------------

---------

Тахогснсритор

ГЖООС

Рис. 1.50. Функциональная схема системы

Если YBX(0 = v|' вых(0 , то 8(') = Ч'ВХ( 0 - '|'ВЬ1Х(0 = 0, а, следова­ тельно, и электрический сигнал сельсинной пары Ucn(t) = 0, то

есть система наведения отработала требуемый азимут и находится в состоянии покоя; если хр вх (?)*У вых (0 . то 6 (0 * 0 , а, следова-

тельно, и Ucu(t)&0, и система наведения разворачивает установку в сторону уменьшения угла рассогласования, то есть отрабатывает требуемый азимут. Сигнал Ucu(t) усиливается на первом каскаде (усилитель) t/ (г), далее алгебраически суммируется с сигналом

Uoc(t) гибкой (демпфирующей) отрицательной обратной связи и поступает на управляемый генератор электрической энергии U0(t)

(второй каскад усиления). После генератора Ur(t) поступает на якорную обмотку исполнительного электродвигателя, скорость ко­ торого С0д(г) пропорциональна Ur(t) и углу рассогласования 8(0-

Редуктор понижает скорость двигателя, а с его выходным валом жестко связана установка. Гибкая отрицательная связь формирует­ ся с помощью тахогенератора.

Рис. 1.51. Структурная схема системы

На рис. 1.51 представлена структурная схема следящей систе­ мы наведения, составленная по передаточным функциям отдель­ ных звеньев.

Р е ш е н и е . Прежде всего, необходимо определить передаточ­ ную функцию разомкнутой или замкнутой системы, откуда затем получить характеристическое уравнение системы. Для этого необ­ ходимо заменить участок прямой цепи структурной схемы, охва­ ченный обратной связью (/Сос), одним эквивалентным звеном:

Кг

Тгр + 1 7 > + 1

К ( р ) =-

1 + - ^

к,ос

Trp + 1

Тдр + 1

После преобразования:

 

к г /

W3(p) =—

(Тгр + Ш аР + 1) +^

ЗЕ-

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

 

K cn tfy f1

 

 

^ р ( Р ) - KcnKyW3(P)' .

к к

 

Р [(Тгр + 1)(7> +1) +

г

0С1; "

 

 

С.

% Р

 

 

 

 

D

 

 

И Л И

w ( p ) =

 

 

 

[(ГгГдр 3 + (Гг + Гд)р2+ (1 + Г^

ос)р ]<р

где О —КспКу ^ .

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

W ( p )

 

D

У/ ( р ) = ---------- р

_

-------------------------------------------------------------------------------------------------------------------------------------------------.

l +Wp(p) [(Г[ТдР3 + (Гг + Гд )р2+ (1+ * г-*ос )p]ip + D

Таким образом, характеристическое уравнение системы при­ мет вид:

 

[ТгТдр 3 + (Гг + Гд )р2+ (1+^ К

0С )PVp + D = 0,

 

 

 

 

с .

 

 

 

 

И Л И

 

 

аър 3 + а2р 2+ ахр + а0 = 0,

 

 

 

где

а3 = TrTaip;а2 = (Тг +Тд)/р;а, = (1 +^ - К 0С)'Р^ 0=D = КсиКу .

 

Пусть для примера:

 

 

 

 

Тг =0,05с; Гд =0,1с \К Г= 5; гр =1000 ; Сс =0,01

 

 

 

 

Следовательно:

 

 

 

 

 

а3 = 0,05 • 0,МООО = 5 ; а2 = (0,05 + 0,1)1000 = 150 ;

 

 

 

а. = (1+ —

К ж ) • 1000 = 1000 +5 105Кос \ а0 = 0,5—

Kv = 250Kv.

i

0>01

ос,

ос

о

0>01

у

у

Видно, что все коэффициенты характеристического уравнения положительны (так как К у >0; К ^ > 0), необходимое условие ус­

тойчивости выполняется. Достаточным условием по Гурвицу для

системы третьего

порядка является выполнение неравенства

^2^1 ^ a3a0*

 

150(1000 + 5

>5-250Ку или 1,5-102 + 7,5104/ГОС > 1,25/<Гу.

Если неравенство обращается в равенство, то система находит­ ся на границе устойчивости:

1,5 • 102 + 7,5 • 104 /TQP = 1,25Ку (граница устойчивости).

Последнее равенство является уравнением прямой, разделяю­ щей область работы системы на две: устойчивую и неустойчивую (рис. 1.52). Прямая проходит через точки с координатами (К у = 0;

Кос = -0,002) и ( АГу =120; /Гос=0).

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

статическую ошибку системы, обу­

К о

словленную влиянием сухого тре­

^Область

ния, зоны неустойчивости и других

/ устойчивости'

факторов, не учитываемых при рас­

 

смотрении

линеаризованной

систе­

 

мы. Статическая ошибка следящей

Область

системы уменьшается (п. 1.7) с уве­

неустойчивости

личением

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

коэффи­

 

циента усиления (Ку). Однако, как

Рис. 1.52. Определение области

следует и

из графика (рис.

1.52),

устойчивости

увеличение

Ку приводит

к

сниже­

 

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

При исследовании САР на устойчивость часто приходится оп­ ределять влияние на устойчивость не одного (двух) параметров, а целого их комплекса. В последнем случае при достаточно высо­ кой степени характеристического уравнения (/? > 4) может быть несколько границ устойчивости.

Пример 3. Определить с помощью критерия Гурвица области устойчивости следящей системы в функции параметра G, если изве­ стны характеристическое уравнение системы и ряд ее параметров:

а5р 5+ а4р 4 + агр ъ+ а2р 2+ ахр + а0 = 0, где а5 = ТГТДТУ;

а4 =ТГТ^ +Гу(Гг + ГД) ; я3 = ТГ+ГД +ГУ+ ТУТЭ(КГ+ 1)G;

а2=1 + (Тэ + Ту)G ; a,=(l +Ty)G\ a0= G , где G = Кси* ^ , Ч'р

Данная следящая система является системой пятого порядка. Ее можно получить, если в предыдущем примере вместо генерато­ ра взять электромашинный усилитель ЭМУ (двухкаскадный гене­ ратор электрической энергии) с постоянными времени его обмоток управления Ту и якорной Тг , а исполнительный электродвигатель считать инерционным звеном второго порядка (п. 1.2.3) с постоян­ ными времени Тэ и Гд. Электромашинный усилитель можно рас­ сматривать как два последовательно включенных генератора со своими постоянными времени Ту и Тг .

Р е ш е н и е . Пусть Тг = 0,12с; Гд = 0,25с; Ту = 0,01с; Кг = 20;

Т3 = 0,01с, тогда а5=0,0003; а4 =0,0337; а3=0,38 + 0,0021G;

а2 = 1 + 0,02G; а, = 1,010; a0 =G . При этом G > 0.

Для системы пятого порядка условие устойчивости по Гурвицу сводится к требованию положительности всех коэффициентов ха­ рактеристического уравнения и выполнения неравенств:

 

а4аг > а5а2 и {а4аъ - a 5a2)(a2at - а 3а0) - ( а 4а, - а 5а0)2 > 0 .

 

Необходимое условие устойчивости выполняется.

 

Проверка первого неравенства а4аъ > а5а2:

 

0,0337(0,38 + 0,0021G) > 0,0003(1 + 0,02G)

или

0,0125 + 6,48-1(T5G > 0 .

 

Неравенство выполняется.

 

Проверка второго неравенства:

[0,0337(0,38 + 0,0021G) - 0,0003(1 + 0,02G)] • [(1+ 0,02G) • 1,01G -

-(0,38 + 0,0021G0G] - (0,0377 • 1,01G - 0,0003G)2 > 0

При равенстве нулю есть границы.

После преобразований получается уравнение третьей степени:

G(l,175-10"6 -G2 -8,73-10‘4G + 78,8-10“4) = 0,

которое и определяет границы устойчивости по назначению пара­ метра G.

В решении последнего уравнения есть три корня: G, = 0; G, =6,93; G3 = 7,35.

Полученные значения G определяют границы устойчивости следящей системы. Чтобы выяснить, какая область значений G яв­ ляется областью устойчивости, необходимо проверить устойчи­ вость системы при любом значении G. Например, при G = 100 ле­ вая часть уравнения отрицательна, следовательно, система неус­ тойчива. Значит, области неустойчивой работы находятся в пределах 6,93<G<7,35 и при G < 0, а области устойчивости соот­

ветствуют значениям 0<G <6,93

и G>7,35.

 

 

 

/

/

/

/

 

На рис. 1.53 выделены

1

/ / J

Область

 

участки числовой оси значений

устойчивости

соответствующие устойчи-*

н/уст.

уст. х

неустой -

 

 

чивости

*

'

'

'

G

вой и неустойчивой работе сле-

о

6.93

7,55

 

I

 

 

дящей системы. Ввиду малой

 

 

 

 

 

С

 

величины интервала 0<G<6,93

Рис. 1.53. Определение участков

 

практическое значение имеет

 

 

устойчивости

 

 

 

 

лишь область G > 7,35.

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

Рис. 1.54. Структурная схема

На рис. 1.54 представлена структурная схема системы. Пара­ метры системы имеют следующие значения:

* сп = 0 Л 5 ^ а д Гг =0,08с ;ГД=0,2с ;С, =0,0101 —

•;/р=148.

К г = 100; К0 = 0,2; Т0 = 0,05с; Ки = 1 с ' ; Ад = 0,1с

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

О(1 + Т0р )(1 + Ка р + ^ )

Wp(p) = _________________________ р______

(1+ 7 » (1 + Тгр)(1+ та р)р + KrK J Qp 2

~ АГПАГ где D = сп г-.

Се г'р

Передаточная функция замкнутой системы после преобразова­ ния и обозначений имеет вид:

ц, (р) <?вых^_

WP

_

blP3 +b2p 2+b]P+b0

.

qBX(p)

l +Wp(p)

 

аьр ь +a4p* +a}p' +a2p 2 +а,р + а0 '

где b, = DKaT0 ■b2 =D(Ka +T0) -, bt = D{1 + K J 0); b0 = DKU.

a5= ТГТДТ0;a4 = TrTa +T0(Tr +Ta ); a3 = Tr +Ta + T0 +T0(KvKn + KaD);

a2= 1 + D{K^ +TQ)\ ax—D(1 + KnTQ) ; a0 = DKn.

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

 

а5р 5 + а4р 4 +а3р 3

+ a 2p 2atp + a0 = 0.

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

 

0,75-100 = 50с"';

 

0,0101

148

а5 = 0,08

• 0,2 • 0,05 = 0,0008; а4 = 0,08 • 0,2 + 0,05(0,08 + 0,2) = 0,03;

аг = 0,08

+ 0,2 + 0,05 + 0,05(100 • 0,2 + 0,1 • 50) = 1,58;

а2 = 1 + 50(0,1 + 0,05) = 8,5 ; ах= 50(1 +1 • 0,05) = 52,5 ; а0 = 1 • 50 = 50.

С учетом полученных значений характеристическое уравнение имеет вид:

0,0008р5+ 0,03р4 +1,58р3 +8,5р2 + 52,5р + 50 = 0.

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

значениях параметров устойчива.

 

 

 

Т абли ца

1.6.

 

 

 

 

 

 

 

as = 0,0008

 

 

«з=1,58

 

« ,= 52,5

0

 

 

«4=0,03

 

 

« ,= 8 ,5

 

а„ = 5 0

0

L

_ 1.58-0.03-8.5-0.0008 _ , ^

,

52,5 0 ,0 3 -5 0 0,0008

О II сГ-

 

L?i

““

0,03

“ 1)«30

1/л -

0,03

- J l.Zr

 

 

'

 

 

 

 

 

 

 

8,5 1,36-51,2 0,03

„ „

 

1.36-50

 

с3 = 0

 

с, = ---------------------------- = 7,36

 

с, = -----------= 50

 

 

 

 

1

1,36

 

 

1,36

 

 

 

 

51,2

7,36 -50 -1,36

^

 

о II

 

 

 

 

«. = --------------------------- = 42

 

 

 

 

 

1

7,36

 

 

 

 

 

 

 

/, = 5 0 -4 2 =50

 

 

/,= о

 

 

 

 

'

42

 

 

 

 

 

 

Пример 5. Проверить по критерию Рауса устойчивость предыду­ щей следящей системы при следующих значениях ее параметров:

* сп = 5 0 * /

.\ТГ=0,12с; Гд =0,25с; Сс = 0,7— ^ ; i„=143;.

сп

/гр ад

град/

Кг = 200; /Со =0,1; Го =0,01с; # ц=1с’1; # д =0,01с;

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