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

книги / Основы автоматики

..pdf
Скачиваний:
1
Добавлен:
19.11.2023
Размер:
41.15 Mб
Скачать

новитьояваоокочаототныеавтоколебания, а при больших началь­ ных условиях - низкочасто!ные автоколебания с большей амплидутой.

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

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

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

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

1. Прямой метод А.М.Ляпунова является общим методом иссле­ дования устойчивости нелинейных систем. Применение его часто вызывает трудности и поэтому в дальнейшем мы его рассматривать не будем.

2. Метод румынского ученого В.М.Попова также предназначен для исследования устойчивости. Подробно этот метод рассматри­ вается в § 12.4.

3. Метод фазовой плоскости (ом.§ 12.3) дает наглядное пред­ ставление не только об устойчивости, но и о характере процесоа в сиотеме.

4 . Метод гармонической линеаризации (си.§ 12.5) - наиболее распространеннный приближенный метод исследования нелинейных систем.

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

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

§ 12.3. МЕТОД ФАЗОВОЙ ПЛОСКОСТИ

Графо-аналитический метод изображения переходных процессов

на фаговой плоскости

был введен в

теорию автоматического управ­

 

 

ления академиком А.А.Ан­

 

 

дроновым в 1944 г . И хо­

 

 

тя этот метод обычно при­

 

 

меняется для исследования

 

 

нелинейных систем

первого

 

 

и второго порядка, а при

 

 

помощи специальных прие­

 

 

мов - лишь некоторых си­

 

 

стем третьего порядка,

 

 

он дает наглядное

пред­

 

 

ставление о характере про­

 

 

цесса управления,

его

 

 

возможных формах,

влиянии

 

 

различных факторов.

 

 

Как известно,

переход­

 

 

ный процесс - это

процесс

 

 

изменения во времени ка­

Ф

 

кой-либо переменной, ха­

 

рактеризующей состояние

Рис.1 2 .6 .Изображение

переходного

системы. Обычно в

каче­

процесса на фазовой плоскости

стве такой переменной

 

 

выбираетоя либо управляемая величина, либо ошибка (рассогла­ сование). Графически переходный процесс изображается в виде некоторой кривой на плоскости, по оси абсцисс которой отклады­ вается независииая переменная - время t , а по оси ординат - переменная, характеризующая состояние системы (рис.12.б ,а ).

При изображении переходного процесса

на фазовой плоокости

по оси абсцисс

обычно откладывается сама

переменная

х (управля­

емая величина,

ошибка), а по оси ординат скорость

ее

изменения

у =

(рис.1 2 .6 ,б ). Вообще же по осям координат

фазовой

плоскооти можно откладывать и любые другие переменные, но так поступают значительно реже.

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

щий колебательный переходный процесс. При

t = О процесс начи­

нается

в точке

I

с положительным начальным

отклонением х. и

положительной начальной

скороотью

у0 = г / т

при

t = 0.

о

На фазовой плоскости начальная точка

I

изобразится

в виде

точки

М, с координатами

х л ,

у.

 

. Далее

на участке процесса 1-2

переменная х

 

 

О

“О

 

 

у =

постепенно

увеличивается,

а

скорооть

уменьшается. При этом изображающая точка

М на фазовой плоско­

сти будет описывать плавную кривую

MjMg. В точке

2 перемен­

ная X

достигает

своего

максимального значения, а

скорость у =

= о. Поэтому на фазовой

плоскости точка

2 изобразится в виде

точки М2, расположенной на оси абсцисс. Эти рассуждения можно продолжить и дальше. Таким образом, каждой точке реального переходного процесса соответствует вполне определенное положе­ ние изображающей точки М на фазовой плоскости.

Кривая, которую описывает изображающая точка на фазовой

плоскости, называется ф а з о в о й

траекторией.

Отметим два основных правила, которые следует учитывать

при построении фазовых траекторий.

1. В верхней полуплоскости изображающая точка всегда дви­

жется слева

направо, а в

нижней полуплоскости - справа

налево.

Это связано

с тем, что в

верхней полуплоскости скорость

О

и, следовательно, сама переменная должна возрастать. £ нижней

полуплоскости у = ^ < 0 и

переменная X должна

убывать.

2 . Фазовые траектории пересекают ось абсцисс обязательно

под прямым углом, так как

на этой оои скорость

у = 0, что со-

ответoraует максимуму или иинииуну оаиой переменой X . Это пра­ вило не раопроотраняетоя на фазовые траектории, которые закан­

чиваются на

оои х

,

не пересекая ее

(ом.табл.12.I ) .

 

 

Для более полного представления об особенностях изображе­

ния процессов на фазовой плоокооти в

табл.12.1 приведены фазо­

 

 

 

вые

траектории,

соответствующие наиболее

 

 

 

характерным переходным

процеооан.

На

 

 

 

основании s i ой таблицы

 

можно сделать сле­

 

 

 

дующие выводы.

 

 

 

 

 

 

 

 

 

 

1. Затухающий апериодический

процеос

 

 

 

изображается на фазовой плоокооти в виде

 

 

 

фазовых траекторий, вливающихся в начало

 

 

 

координат (табл.12.I ,а )

 

или в

точку,

со­

 

 

 

ответствующую новому установившемуся

зна­

Рио.12.7. Отображе­

чению (табл.1 2 .1 ,6 ).

 

 

 

 

 

ние переходного про-

 

2. Расходящийоя апериодический про­

пеооа (рис.12.4) на

цесс

изображается на фазовой плоокооти

фаэовую плоскость

 

в виде фазовых траекторий, удаляющихся

 

 

 

 

 

 

от начала координат (таб л .1 2 .1 ,в,г).

 

3. Изображением затухающего колебательного процеооа являют­

ся сходящиеоя спиралевидные фазовые траектории (рис.12.б,

 

табл .12 .1,д).

 

 

 

 

 

 

 

 

 

 

4. Изображением расходящегося колебательного процесоа

 

являются расходящиеся

спиралевидные

фазовые

траектории

 

(табл .12 .1,3).

 

 

 

 

 

 

 

 

 

5. Изображением периодического колебательного процесса

 

являютоя замкнутые

фазовые траектории

(табл.12.1,е,ж ).

Обра­

зованный ими контур обычно называется

ц и к л

о м.

Для сину­

соидальных колебаний цикл имеет форму эллипса, так как,еоли

х = a sin си t ,

то

у =

= а си cos to t

.

Эти два

урав­

нения и есть

уравнения эллипса с полуосями а я ш . Боли же

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

На рис.12.7 и 12.8 изображены фазовые траектории, соответ­ ствующие возможным процессам в нелинейных системах (см.рис.12.4 и 1 2 .5). Здеоь мы встречаемой с двумя видами фазовых траекто­ рий, которых не было в линейных системах. Контур, к которому

в пределе отягиваются все

фазовые траектории,

называется ^устой­

чивым предельным циклом (рио.12.8). При этом

фазовые траекто­

рии, расположенные внутри

предельного цикла,

являются раоходя-

Рис.12.8. Отображение переходного процесоа (рио.12.5) на фазовую плоокость

щимися, а расположенные вне его - сходящийся. Такой предель­ ный цикл соответствует устойчивый периодический колебаниям в нелинейной сиотеме, названный нами автоколебаниями. На фазовой плоскости может быть и неустойчивый предельный цикл (показав пунктиром на рис.12.7 и 12 .8,6). Такой предельный цикл соот­ ветствует неустойчивым периодическим колебаниям в сиотеме.

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

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

Пример 12.I . Исследуем один из возможных вариантов авто­ матической системы стабилизации углового положения искусствен­ ного спутника Земли (ИСЗ). Эта система обычно строится по схе­ ме, изображенной на рис.1 .5 , и включает элементы, измеряющие угловое положение спутника (например, гиростабилизированную платформу), формирующий элемент (в простейшем случае - усили­ тель) и исполнительные элементы (например, газовые реактивные двигатели). Коли пренебречь влиянием внешней среды, уравнение управляемого объекта (ИСЗ) может быть записано в виде

d2u>

(к . п

э - j j r

где 0 - момент инерции спутника относительно оси вращения; Ч>- угол поворота ИСЗ (управляемая величина);

Т а б л и ц а 12.1

а)

Д)

386

б)

г )

з)

387

М- момент (управляющее воздействие) со отороны исполнитель­ ного элемента.

Структурная охема

системы изображена

на рио.12.9, где

к0=

= -j- . В дальнейшем для простоты мы будем полагать, что про­

граммное значение угла поворота ИСЗ ij^ =

0 и тогда

 

 

 

 

 

 

 

 

ЛЧ, = Ч’лр-Ч’ = -Ч’

 

 

 

 

(12.2)

Допуотим, что в управляющем уотройотве реалиэуетоя следую­

щая зависимость:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

И

= Мо Ф { Ч

' Т $

 

 

 

 

(12.3)

 

 

 

 

 

 

 

 

 

 

 

 

 

М0 -

постоянная положительная величина (момент);

 

 

 

 

некоторая функция от угла

поворота

у

и

угло-

 

'

 

вой скорости

 

.

 

 

 

 

 

 

 

 

Решая уравнения

(12 .I)

-

(12.3) оовмеотно,

получим диффе-

ренциальное уравнение

замкнутой системы

 

 

 

 

 

 

 

 

 

 

 

d ?

 

к Ф ft.!?

 

 

 

 

(12.4)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

ма

 

 

 

 

 

 

 

 

 

 

 

 

 

 

к = - J - .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Уравнение (12.4) имеет второй порядок, и будь оно линейным,

решение его

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

труда. Но входящая в

 

 

 

 

 

 

 

 

 

это уравнение

функция

 

 

 

 

м

 

1

 

 

 

Ф(

tj> ,

 

)

в общем слу­

н - ф —“

 

 

 

 

 

чае,

как будет

показано даль­

Ф(ф.РУ')

э ?

 

 

 

 

 

 

 

 

 

 

 

ше, является

нелинейной,

по­

Рис.12.9. Структурная

схема

этому

 

и уравнение (1 2 .4 ),

вообще

говоря,

является

не­

 

к

примеру 12 .I

 

 

 

линейным. Поэтому мы

будем

 

 

 

 

 

 

 

 

 

решать его методом фазовой плоскооти. Для этого введем новую

переменную -

угловую скорооть

Q =

. Но так

как

на фазовой

плоскости нет

переменной

t

,

сделаем следующее

преобразование:

 

 

 

dzy

~

dO.

~

dCi

dip

 

dO.

 

 

 

 

 

 

 

 

d t z

dt

dy

dt

 

dip О.

 

 

(12.5)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда

уравнение (12.4)

принимает вид

 

 

 

 

 

 

 

QdO. =- к

Cl) dip

(12.6)

Уравнение (12.6) и есть дифференциальное уравнение фаговых траекторий. При построении фазовой плоскости переыеннув if мы будем откладывать по оси абсцисс, а переменную Q - по оси ординат.

Рассмотрим влияние различных законов управления на про­ цессы в оиотеме стабилизации ИСЗ.

Рис.12.10. Фазовые траектории к примеру 12.I при использовании линейного исполнительного устройства

I .

Пуоть управляющее

устройство

предотавляет

ообой обыч­

ный линейный усилитель,

т .е .

Ф ( if , Cl

) - линейная

функция

угла

поворота ij>

(рис.12.10,а ):

 

 

 

 

 

 

 

Ф (if,Q) = /<, if

 

(12.7)

Тогда

дифференциальное уравнение

(12.6)

можно записать в виде

 

 

 

 

QdQ. =-кк, y d y

,

(12.8)

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

торий

 

 

 

Q

 

 

 

 

 

 

 

 

 

 

 

 

 

(12.9)

 

 

 

 

 

 

 

 

 

где

С( -

постоянная

интегрирования, определяемая начальными

условиями

if = if0

,

Q

= Q 0

при

( = 0.

 

Задаваясь различными начальными условиями, получим оемей-

ство фазовых траекторий,

или,

по выражению А.А.Андронова, ф а -

з о в ы й

п о р т р е т

системы. Нетрудно убедиться, что

уравнение (12.9) еоть параметрическое уравнение эллипса. По­ этому фазовый портрет системы будет представлять ообой семейотво концентрических эллипсов (рис.12.10,б). Следовательно, в сиотеме существует периодический колебательный процесс. Так

как мы рассматриваем линейную систему, то наличие таких коле­ баний показывает, что система находится на колебательной грани­ це устойчивости. В этом можно убедиться и другим споообом, ис­ пользуя теорию линейных автоматических сиотем. При Ф( у , С2 )=

-передаточная функция разомкнутой оистемы

W (p) = А ^ 1 о _ =

РР

Характеристическое уравнение замкнутой оистемы

д{р) = рг + к,к =О

имеет два

чисто мнимых

корня

р1г = ± j -/к^к

Это и

еоть признак

колебательной границы устойчивости.

2.

Посмотрим теперь,

как будет вести оебя система стаби­

лизации ИСЗ, еоли в качестве уоилителя использовать трехпози­ ционное поляризованное реле. Характеристику реле примем такой, как на рио.12. 1 ,а .

Управление спутником в этом олучае будет осуществляться по

следующему

принципу. При малых углах

отклонение ИСЗ,

лежащих

в пределах

зоны нечувствительности реле - Ь < ^> < + Ь

,

управ­

ляющий момент М равен

нулю, так

как

Ф( 41 , Q ) =

0.

Но как

только угол

отклонения

превышает

пороговое значение |^ |

э» b ,

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

(рис.12 .I I ) .

Область

I располагается правее

вертикальной пунк­

тирной линии, проходящей чере точку . В

этой облаоти

if>+b,

Ф( ч , Q )

= + I (см .рис.12. I ,а )

и поэтому нелинейное

диф­

ференциальное

уравнение (12. 6) превращается в линейное

 

 

 

 

Q d Q = - k d y

(1 2 .Ю )

 

Проинтегрировав

его, получим уравнение

фазовых траекторий

 

 

 

-y - + Avf =

с г

(12.II)

где

С2~ постоянная

интегрирования.

 

 

Нетрудно убедиться,

что это есть уравнение параболы, симмет­

ричной относительно оси абсцисо и

обращенной вершиной вправо.

Поэтому фазовый портрет

системы в

области I при различных

С

представляет собой семейство парабол(рис.12. I I ) .

2