запишем уравнение (8.11) в относительной форме
d2x |
/ |
к dx \ |
= |
(8.12> |
d\ 2 |
+ F , z[x |
---- |
1 |
dx) |
|
|
где |
/гш/г, . |
|
|
|
t = |
|
|
|
|
Jz ’ |
|
Jz |
|
Оба рассмотренных примера привели нас к нелинейным диф ференциальным уравнениям, которые не могут быть исследо ваны излагавшимися выше методами линейной теории. Для ис следования процессов регулирования в релейных системах ши рокое применение получил метод фазовой плоскости, отличаю щийся достаточно хорошей наглядностью. Однако естественно, - что этот метод может быть использован для исследования си стем, описываемых дифференциальными уравнениями второго порядка. Исследование систем более высоких порядков анало гичными приемами вызывает необходимость отображения про текающих в них процессов в фазовом пространстве, и нагляд ность метода резко уменьшается. Рассмотрим основные поло жения метода.
§ 8.2. МЕТОД ФАЗОВОЙ ПЛОСКОСТИ
Методом фазовой плоскости исследуется поведение откло нений координат или переменных системы от их невозмущен ных значений. В данном случае в качестве невозмущенного движения принимается состояние покоя или равновесия систе мы, когда все ее координаты равны нулю (или постоянны).
Фазовой плоскостью называется плоскость, на которой по двум осям координат (х, у) откладываются переменные, харак теризующие состояние системы. В качестве таких переменных
обычно принимаются |
регулируемая |
величина |
(или отклонение) |
X и ее производная ѵ = |
d x |
|
|
|
---- . |
|
|
|
dt |
|
|
|
Пусть поведение некоторой системы описывается нелиней |
ным уравнением второго порядка |
|
|
|
d2x |
( |
dx |
(8.13) |
|
dt2 |
I |
dt |
|
|
Введя обозначение |
— |
dx |
|
уравнение в виде |
— .представим это |
|
|
dt |
|
|
двух уравнении первого порядка: |
|
|
dx |
~ = F { x , y ) . |
(8.14) |
|
|
dt
Уравнения (8.14) дают закон изменения во времени регулируе мой величины (отклонения) х и скорости у. Оба уравнения можно рассматривать как параметрические уравнения некото рой кривой в плоскости X, у, где параметром является время t. Деля уравнения (8.14) почленно друг на друга (исключая параметр /), получим дифференциальное уравнение этой кри вой на фазовой плоскости
Осуществив интегрирование (8.15), получим уравнение интег ральной кривой у = F (х, с) на фазовой плоскости. Кривая у = F(x, с) называется фазовой кривой, где с — постоянная интегрирования, определяемая из начальных условий.
Бесчисленное множество возможных начальных условий оп ределяет бесчисленное количество фазовых кривых. Фазовые кривые называются также фазовыми траекториями так называе мой изображающей точки т. Изображающая точка — это точ ка с координатами х т= х т{() и y m— y m{t). Положение изо бражающей точки на фазовой траектории указывает на значе ние отклонения х и скорости у в данный момент времени. По скольку при определенных начальных, условиях дифференци альные уравнения имеют единственное решение, то изображаю щая точка будет двигаться по вполне определенной траектории и при отсутствии внешних сил не может перейти на другую траекторию. Примерный вид фазовых траекторий для разных случаев приведен на рис. 8.9. На рис. 8.9,а и б показаны фазовые траектории устойчивой системы. В устойчивой систе ме, независимо от того, описывается ли ее поведение линейны ми или нелинейными уравнениями, переходный процесс зату хает, следовательно, при t ос, х -> 0 и у -- 0.
Это означает, что у устойчивой системы изображающая точ ка по фазовой траектории всегда будет двигаться к началу ко ординат. Если в системе возможны незатухающие колебания, т. е. уравнение (8.13) системы имеет периодическое решение, то изображающая точка будет непрерывно двигаться по неко торой замкнутой кривой (рис. 8.9,б) (замкнутому циклу), так как значение отклонений и скорости повторяется через период. И наоборот, замкнутым фазовым траекториям будут соответ ствовать периодические колебания в нелинейной системе, кото рые мы в дальнейшем будем называть автоколебаниями. Время одного оборота изображающей точки, очевидно, будет равно периоду автоколебаний. В неустойчивой системе изображающая точка уходит от начала координат, т. е. х и у нарастают бес предельно (рис. 8.9,г, д). Теперь можем провести исследова ние приведенных выше двух релейных систем автоматического регулирования.
Исследование проведем для различных видов релейных ха рактеристик. Отдельно оценим влияние запаздывания в релей ном усилителе и влияние корректирующих цепей (тахогенератора в следящей системе и скоростного гироскопа в системе ста билизации) .
Рис. 8.9. Примеры фазовых траекторий:
аи б — фазовые траектории устойчивых систем; в — фазовые траектории при наличии незатухаю щих колебаний в системе; г и б — фазовые тра
ектории неустойчивых систем
1. Релейно-контактная следящая система с идеальной релейной характеристикой
Характеристику, приведенную на рис. 8.10, будет иметь ре лейное устройство с очень высокой чувствительностью. Пола гаем, что реле срабатывает мгновенно и, следовательно, ts = 0, а тахогенератор отключен.
Структурная |
схема |
упрощенной |
системы |
приведена |
на |
рис. 8.11. |
|
|
(см. |
рис. |
8.10) |
обычно |
за |
Уравнение нелинейного элемента |
|
писывают в форме: |
|
|
|
|
|
|
/ - ' ( в ) |
= S i g n s . |
( 8 . 1 6 ) |
1 |
|
|
|
|
|
|
|
|
|
|
л . |
(Tptpp |
|
|
|
|
|
|
|
|
Р и с. 8.10. Характери |
Рис . 8 11. Структурная |
схема |
упрощен |
|
стика идеального |
ре |
ной следящей системы с идеальным релей |
|
лейного элемента |
|
ным элементом |
без запаздывания |
|
F ( e ) = sign г |
|
|
|
|
|
|
|
Данное выражение является символическим обозначением скачкообразной функции, изображенной на рис. 8.10, и читает ся «знак е ». Выражение (8.16) можно заменить эквивалент ным выражением вида:
С учетом |
введенного |
обозначения |
уравнение |
системы (см. |
рис. 8.11) |
в относительных единицах при f = 0 будет иметь вид: |
|
|
|
d 2x |
dx |
, |
|
|
|
|
|
|
d i 2 |
------Ь sign л:=0. |
|
( 8 . 1 7 ) |
|
|
|
dx |
|
|
|
|
Введем |
обозначение sign .ѵ = |
— |
следовательно, |
х = 1 |
при X |
0 |
и |
/. = — 1 при л'^>0. Это позволяет записать урав |
нение (8.17) |
в виде системы двух уравнений: |
|
|
|
|
|
|
dy 1 ■^ — х ; |
|
’ |
( 8 . 1 8 ) |
|
|
|
|
dx |
|
|
|
|
Исключив из (8.18) независимую переменную т, получим уравнение, связывающее х и у.
d y _•/. — у |
dx |
у |
(8.19) |
|
После интегрирования (8.19) |
при |
начальных условиях т = 0, |
X = х0 и у = уо получим уравнения фазовых траекторий |
* = * 0+ Уо- У + |
V __ X |
* In —---- • |
|
|
.1' — х |
При |
* = 1 уравнение фазовой траектории будет |
|
|
X = Jc0 + ^o— У + In У- ^ - . |
( 8.20) |
|
У - 1 |
|
При |
X — — 1 соответственно будем иметь |
|
|
X = л'0Н- Уо — у — Іи 1 1 Уо • |
(8.21) |
|
1+ У |
|
На рис. 8.12,0 по уравнениям (8.20) и (8.21) построены фа зовые траектории при различных начальных условиях для х = 1
и х = — 1.
Рис. 8.12. К пояснению построения фа зовых траекторий:
а — фазовые траектории при различных начальных условиях для х= 1 и х = —1; б — фазовая траектория упрощенной системы без запаздывания с идеальным релейным элементом
Фазовая траектория, отображающая конкретный процесс в
системе |
со |
структурной схемой (см. рис. 8.11) |
при |
каких-либо |
начальных |
условиях, составляется из отрезков |
траекторий для |
х=г- 1 |
и |
X — — 1, изображенных на рис. 8.12,а. |
Переход с |
траектории |
для |
х — 1 па траекторию |
для |
х = |
— 1 будет |
происходить на |
оси у, когда х изменяет |
знак |
и соответственно |
X изменяется от значения + 1 до значения — 1 или наоборот. При достижении фазовой траектории оси у в системе происхо дит переключение реле. В связи с этим ось у является линией
переключения. На рис. 8.12,6 |
построена фазовая |
траектория |
для начальных условий у -= 0, |
н х |
= |
оаа. При этих начальных |
условиях по |
уравнению |
(8.20) |
построен |
участок |
фазовой |
траектории айЬ\, далее участок |
b\axb-2 построен по |
уравнению |
(8.21), участок Ь»афз — вновь |
по |
уравнению |
(8.20). |
у = |
ob\, |
Начальные условия для |
второго |
участка |
.ѵ = |
0, |
для третьего |
участка х = |
0, |
у |
= |
ob2 и т. д. |
Продолжая |
пост |
роение фазовой траектории по участкам, можно заметить, что получается кривая, сходящаяся к началу координат. Это озна чает, что система устойчива и в ней имеет место колебательнозатухающий переходной процесс.
Поясним физические процессы, протекающие в следящей си стеме (см. рис. 8.4) с учетом упрощающих предположений (от сутствует тахогенератор, и реле срабатывает мгновенно). Пусть в мо.мент времени / — 0 мы скачком повернули входную ось
следящей |
системы |
на |
некоторый угол |
f0- |
Это |
будет |
означать, |
что при t |
= 0 |
система имеет |
отклонение |
от |
заданного поло |
жения на угол |
г = |
уо |
(отрезок |
оап на |
оси |
абсцисс). |
В резуль |
тате отклонения срабатывает поляризованное реле и один из контакторов (k\ или k2). Сервомотор запускается, и скорость его нарастает. Отрезок аф\ и характеризует нарастание окорости у и уменьшение первоначального отклонения Хо. При х=0 скорость у равна отрезку оЬІШИмея запас кинетической энер гии, якорь сервомотора проскочит нейтральное положение, и в системе возникнет рассогласование противоположного знака. В это'М случае изменится полярность напряжения, подводимого к обмотке поляризованного реле. В результате срабатывания по ляризованного реле и соответствующего контактора (/г, или /г2) изменится полярность напряжения на якоре сервомотора. Послед ний войдет,в режим торможения противотоком (отрезком йіПі). За время торможения скорость уменьшится до нуля, однако на капливается рассогласование противоположного знака (отрезок оаі). Сервомотор останавливается и затем начинает вращать ся в противоположную сторону, выбирая накопленное рассог ласование оа.\. Этот режим вновь является режимом пуска (от резок а\Ь2) . Скорость отработки рассогласования возрастает, а само рассогласование уменьшается. Далее процесс протека ет аналогично рассмотренному выше. В результате выходной вал колеблется около своего нейтрального положения. Однако
с каждым колебанием скорость, благодаря наличию потерь энергии в системе, уменьшается. Уменьшается амплитуда ко лебаний выходной оси (рассогласование), система стремится к состоянию покоя. Изображающая точка по спирали прибли жается к началу координат. Фазовая траектория имеет вид спирали с изломами на оси ординат, являющейся линией пере ключения реле. Изломы соответствуют моменту переключения реле и изменению знака момента сервомотора. Найдем зави симость выходной величины х и скорости ее изменения у от времени.
Воспользуемся первым уравнением системы (8.18). Так как х постоянна для правой и левой полуплоскости, то мы можем решить данное уравнение, считая х = const, и затем учесть изменение знака х при изменении знака х:. Решение уравне
ния |
у = х при начальных |
условиях т= О, у = у0 будет |
|
d‘с |
|
|
|
иметь следующий вид: |
|
|
(8.22) |
|
у = j '0e - '! -f х(1 — е - т ). |
Проинтегрировав полученное выражение, найдем х как |
функ |
цию |
т при начальных |
значениях |
т = 0, х = х0: |
|
|
х = х0 + уо (1 |
- е- " X |
[т — (1 — е~' )].* |
(8.23) |
На основании (8.22) и (8.23) методом припасовывания началь ных условий можно построить графики изменения х и у в функ ции безразмерного времени т.
Рис. 8.13. Временные характеристики системы:
а — переходные функции х (т) и у (т); б — изменение * (график переключе ния реле)
На рис. 8.13 построены графики изменения х и у в функ ции времени. Эти графики соответствуют фазовому портрету
* Если из (8.22) и (8.23) исключить параметр т, то получим уравнения фазовых траектории (8.20) и (8.21).
системы (см. рис. 8.12). На этом же рисунке 8.13 показан про цесс изменения х, т. е. процесс переключения реле. Графики показывают, что частота переключения непрерывно растет, стремясь в пределе к бесконечности. Таким образом, в состоя нии покоя рассматриваемая система будет характеризоваться колебаниями бесконечно большой частоты и бесконечно малой амплитуды. Этот результат объясняется допущением, сделан ным выше, о возможности мгновенного срабатывания релей ного элемента системы.
2. Система стабилизации углового положения космического летательного аппарата с идеализированными релейными характеристиками
Рассмотрим процессы в упрощенной системе стабилизации углового положения летательного аппарата с релейным управ лением, структурная схема которой приведена на рис. 8.14.
|
J 1 |
г |
|
К |
iF |
|
|
|
|
|
|
|
|
|
J p 2 |
|
|
|
Рис. 8.14. Структурная |
схема |
упро |
|
|
щенной |
системы |
стабилизации |
углового |
|
|
положения космического |
|
летательного |
|
|
|
аппарата |
|
|
|
|
|
В относительных |
переменных х ■ |
|
|
dx _ |
J г |
db |
|
— 9; у |
k |
dt |
уравнение системы будет иметь вид: |
|
к |
dt |
|
|
|
|
|
|
d2x |
_ |
|
0, |
|
(8.24) |
|
------ \-F (х) = |
|
|
dt2 |
' |
|
|
|
|
|
где F(x) — характеристика |
с зоной |
нечувствительности |
(см. |
рис. 8.2,6).
Величину зоны нечувствительности для релейной характе
ристики F(x) обозначим о = — а, где а — угол рассогласо- k
вания, при котором происходит включение реактивного сопла. F(x) принимает следующие значения:
F(x) = |
— 1 |
при |
X |
< |
— |
а ; |
|
0 при |
—о < |
х < э; |
(8.25) |
|
1 |
П р и |
X |
> |
о . |
|
|
Обозначим F(x) = |
— у, |
считая, что у = |
1 |
при X < |
— а, у = |
= — 1 при А ')>а |
И У. |
= 0 при — |
О |
X < о. |
Уравнение |
(8.24) представим в виде двух уравнений:
dy
(8.26)
öx
(8.27)
Деля первое уравнение на второе, найдем дифференциальное уравнение фазовых траекторий:
dy _ ■/.
(8.28)
dx у
Интегрируя (8.28), получим уравнение фазовых траекторий
У2Уо2= 2х { к -г х0). |
(8.29) |
При у ф 0 — это уравнение парабол, а при у = |
0 — уравне |
ние прямых. |
|
Рис. 8.15. Семейства фазовых |
траекторий: |
а — о фО; б— з = 0 |
|
На рис. 8.15,а и б построены семейства фазовых траекторий |
при о Ф 0 и а = 0. Фазовые траектории |
при любых началь |
ных условиях х0 и уо оказываются замкнутыми, так как пара
болы симметричны относительно оси х, а фазовые |
траектории |
в зоне нечувствительности — прямые, параллельные |
оси абсцисс. |
Интегрирование уравнений (8.26) и (8.27) дает |
законы изме |
нения отклонения и скорости в |
предел ах‘каждого |
участка: |
У = |
+ у 0; |
|
(8.30) |
х==~2 |
“Ь -^о- |
|
(8.31) |
При *■Ф 0 скорость изменяется по линейному закону, откло нение— по параболическому. При х = О скорость постоянная, отклонение изменяется по линейному закону. Графики времен ных зависимостей ij(f) и x(f)
приведены на рис. 8.16,а и б. Данной системе свойственно бесконечное множество перио дических колебаний, соответ ствующих бесконечному мно жеству начальных условий. Период колебаний определяет ся начальными условиями. Си стема с такими динамическими свойствами к эксплуатации не пригодна как неустойчивая.
|
|
|
|
|
3. |
Процессы в релейных |
|
Р ы с. |
S.16. |
Законы изм енения |
ско |
системах при |
наличии |
|
запаздывания |
|
рости |
y{t) |
и |
отклонения |
x(t): |
|
|
а — а = 0 ; |
б — 5 = у ( ) |
В |
практике |
всегда имеется |
|
|
|
|
|
|
|
|
|
|
'запаздывание |
при |
прохожде |
нии сигнала через релейный усилитель. Изучим влияние запаз дывания в срабатывании реле и переключении регулирующего органа в системе стабилизации углового положения летательного аппарата. Для отдельных участков уравнения фазовых траек торий для данного случая и законы изменения координат х и у
во времени будут прежними [см. |
(8.29), (8.30) и (8.31)]. |
Влияние запаздывания будет сказываться в том, что изме |
нение X |
с — 1 на 0 и с 0 иа 1 |
будет происходить не на пря |
мых X = |
± з |
(например, точка |
1), а правее, например, в не |
которой |
точке 2 |
(рис. 8.17), так |
как благодаря запаздыванию |
в отключениях и включениях соплового аппарата скачкообраз ное изменение вращающего момента будет происходить после того, как ]х | достигает значения а.
Зная величину запаздывания, всегда можно определить точ
ки переключения. |
Для небольших |
Дт (Дт < 1) |
(см. |
рис. 8.17,а) |
уравнения л и н и й |
переключения (точнее, отключения |
и включе |
ния) |
могут быть приближенно взяты в виде прямых: |
|
|
|
|
у — і з-Г — л-, |
|
|
(8.32) |
|
|
|
Д-с |
|
|
|
где |
Д т характеризует наклон линии переключения. |
|
|
Фазовые траектории, построенные с учетом |
запаздывания |
для каких-либо конкретных начальных условий |
(см. |
рис. |
8.17,а |
и б), будут всегда удаляться от |
начала координат. |
Это |
озна |