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

книги из ГПНТБ / Основы авиационной автоматики учеб. пособие

.pdf
Скачиваний:
50
Добавлен:
22.10.2023
Размер:
13.63 Mб
Скачать

Оптимальные процессы во времени приведены на рис. 8.25,6 и в. Графики •/(т) характеризуют закон изменения вращающего момента (или напряжения) двигателя.

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

Указанный способ обеспечения минимального времени регу­ лирования в системе является следствием более общего прин­ ципа максимума Понтрягина Л. С., из которого вытекает [1], что для получения максимального быстродействия при любом порядке линейного дифференциального уравнения объекта все­ гда необходимо использовать релейное управление.

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

менять

знак

отклонения.

Последнее

положение

доказано

Фельдбаумом А. А. [5].

что

оптимальный процесс при С

Из рис. 8.25,п

видно,

- const

может

быть только

при какой-либо одной

группе на­

чальных

условий

(например,

когда все

начальные

положения

изображающей точки лежат па отрезке траектории хоцщ). При начальном положении на д'озиіз и на х^пц процессы уже не будут оптимальными. В первом случае процессы будут с пере­ регулированием, во втором — будет скользящий режим.

Для того чтобы при любых начальных отклонениях переход­ ный процесс был оптимальным, линией переключения должна быть сама фазовая траектория отп\р (и соответственно оп'т\'р' — в противоположном квадранте).

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

Уравнение системы стабилизации углового положения кос­ мического летательного аппарата с нелинейной корректирую­ щей цепью будет иметь вид:

Зс + F[x + /(х )] = 0.

(8.35)

Функция f нелинейного корректирующего сигнала зависит от скорости X и свойств реле. Следовательно, для получения опти­

399

мального процесса релейная система должна снабжаться спе­ циальным вычислительным устройством, образующим функ­ цию f. Для уравнения (8.35) линией переключения ріщпо будет парабола

X + — з/2 = 0.

2

Следовательно, искомая функция

 

f[x) =

(jr)2signA\

(3.36^

Множитель sign.v указывает, что квадратичная функция про­ изводной должна быть нечетной функцией (для образования линии переключения р'т\'п'о).

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

тимальных, линию переключения можно

взять в виде прямой,

и оптимальное значение порции сигнала

от тахогенератора бу­

дет при С= 0,3 [1].

 

§ 8.4. МЕТОД ГАРМОНИЧЕСКОЙ ЛИНЕАРИЗАЦИИ

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

ской линеаризации, разработанный

академиками Крыловым

Н. М. и Боголюбовым Н. Н.

Л. С., распространившим

Метод развит далее Гольдфарбом

его на системы автоматического регулирования для исследо­ вания устойчивости и автоколебаний нелинейных систем, и По­ повым Е. П., осуществившим его дальнейшее обобщение и раз­ витие (в частности для исследования качества регулирования нелинейных систем). Метод гармонической линеаризации по­ зволяет исследовать устойчивость систем и исследовать уста­ новившиеся процессы, т. е. определить амплитуду и частоту автоколебаний и получить зависимость основных параметров автоколебаний от параметров системы. Метод является прибли­ женным, но довольно широко применяется, так как по своему математическому аппарату является весьма простым и, кроме того, обеспечивает достаточно высокую для практики точность. Чтобы пояснить сущность метода гармонической линеаризации, представим нелинейную систему в виде двух частей: нелиней­ ной и линейной (рис. 8.26,а).

400

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

X = А sin t,

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

Нелинейная У,

Линейная X

часть

часть

о)

 

б)

Рис . 8.26. К пояснению метода гармонической линеаризации:

а — представление нелинейной системы в виде нелинейной и линейной частей; б

структурная схема нелинейной системы с од­ ним нелинейным элементом

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

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

Этот метод позволяет получить тем более точные результа­ ты, чем ближе линейная часть по своим «подавляющим свой­ ствам» к фильтру низких частот.

Структурную схему нелинейной системы с одним нелиней­ ным элементом представим в виде рис. 8.26,6.

Запишем уравнение линейной части системы

 

X — W (D) у = Q ( £ > )

ba + b^DA-------b bmDm

У, (8.37)

 

 

P {D) У =

c0 + ^ ö - j ----- cnDn

где D

d

символ дифференцирования.

 

 

 

dt

26. И зд. № 5312

401

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

 

 

 

y = F [ - x )

=

- F{x).

 

(8.38)

Знак

минус

обусловлен

учетом

уравнения

рассогласования

е = л;вх — X

= — X ,

так как

входной сигнал

 

полагаем рав­

ным нулю (хвх — 0).

 

сигнал

 

л;= уі sin оі ^

с

неизвестными

Пусть гармонический

ш

нам амплитудой А п частотой

автоколебаний

действует

на

вход

нелинейного

элемента. В

этом случае его

выходная

ве­

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

 

 

 

 

 

у = — F(x) = — F(i4 sin <о£).

 

(8.39)

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

 

у =

F [А sin u>f) ==

— i]) (А) sin id*;,

(8.40)

где

4>(Л) = f

F (А sin ш t) s i n

d[y>t)

коэффициент ря-

 

J .

 

 

 

 

да

Фурье или амплитуда первой

гармоники выходной

величи­

ны нелинейного элемента. ф(/1) зависит от вида характеристи­ ки нелинейного элемента и амплитуды сигнала А, действующе­ го на его входе.

Рассмотрим некоторые типы нелинейных характеристик и со­ ответствующие нм зависимости ф(4).

Если на вход идеального нелинейного элемента (рис. 8.27,а) действует гармонический сигнал а' = і4 зіпшЛ то на его выходе

^ y ip(f) ^

F(Asincöt)

I

\

I

L

 

и

ъ

 

І

S)

 

ФЩ

 

 

 

 

 

■4В.

О)

0

 

1К

 

А

 

 

Ь)

Рис . 8.27. Идеальный релейный элемент:

а — характеристика и

зависимость

х = А sin coif;

б — периодический несинусоидальный сигнал ре­

лейного элемента

F{A sin о> t)

и закон

изменения

амплитуды

первой

гармоники

ф(.Л) во

времени;

в

— зависимость ф(4)

 

402

получим 'Последовательность

знакопеременных

импульсов

(рис. 8.27,6). Амплитуда первой гармоники такой

периодиче­

ской функции постоянна и не

зависит от А (рис.

8.27,в). Она

равна:

 

 

 

45

АI )

ФИ) =

(8

ТС

 

Для нелинейного элемента с зоной нечувствительности (рис. 8.28,а) выходной сигнал представляет собой серию прямоуголь­ ных импульсов, следующих с паузой друг за другом (рис. 8.28,6). Амплитуда первой гармоники в этом случае зависит от амплитуды входного сигнала:

45

(8.42)

ф(А) = — V А 2- я2-

itА

 

Характеристика ф(А) приведена на рис. 8.28,в.

/Ф(А

 

/ F(AsLTi(Ot)

1 Л

/

/ '

'

" Т

 

 

 

 

У

5)

 

Р и с. 8.28. Релейный элемент с зоной нечувствитель­ ности:

а — характеристика н зависимость х = А sin ш і\ б — периодический несииусоидальный сигнал релейного элемента F (А sin ші) и закон изменения амплиту­ ды первой гармоники ф (Л> во времени; в — зави­

симость ф

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

Если входной сигнал нелинейного элемента запишем в комплексной форме

x = A e Jmt,

(8.43)

26*

403

а амплитуду первой гармоники нелинейного элемента предста­ вим как

3> = ф ( Л ) е Уо,<,

( 8 . 4 4 'і

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

 

V

Ф(Л)

 

 

 

 

=

.

( 8 . 4 5 )

или

X

А

 

 

 

 

 

 

 

У =

Л

 

( 8 . 4 6 )

 

 

 

 

Понятие

передаточной функции

^ „(4 )

позволяет

заменить

данное нелинейное звено для

периодического режима и по­

стоянной

амплитуды А некото*рым линейным эвеном

с коэффи­

циентом усиления W H(A), так как нелинейное уравнение (8.33) заменяется линейным (8.46). Такой процесс сведения нелиней­ ной зависимости к линейной носит название гармонической ли­ неаризации нелинейных зависимостей. Это название вытекает из разложения нелинейных колебаний на гармонические состав­ ляющие (гармоники) и пренебрежения высшими гармониками на выходе нелинейного элемента и линейной части системы. В результате гармонической линеаризации мы получаем своеоб­ разное линейное звено, коэффициент усиления которого ІУДА) зависит от амплитуды входного сигнала А.

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

ствительные числа, так как они не

создают фазовых сдвигов

при прохождении гармонических сигналов и ф(А)

для них чис­

ло вещественное.

 

характеристиками

(типа

 

Для звеньев с неоднозначными

рис. 8.3,а, б, в) гармонический выходной сигнал

нелинейного

элемента (первая гармоника) имеет

сдвиг

фазы по отношению

к сигналу на входе

 

 

 

 

 

 

у=Ъ(А) еЛш,- 9,Д)1 =

ф(А)еМ/|)

= ф (/А) e/W ,

(8.47)

где

ф(/А) — комплексная

амплитуда первой гармоники;

 

ср (А) — сдвиг фаз между выходным и входным

сигна­

 

лами.

 

 

 

 

 

та

Следовательно, передаточная функция нелинейного элемен­

с неоднозначной характеристикой

есть

число

кташлекеное

 

 

Фи А)

 

 

(8.48)

 

W H( j A ) = Z -

А

 

 

 

 

X

 

 

 

Как и для нелинейных элементов с однозначными

характери­

стиками, выражение (8.48)

можно

рассматривать

как

комп­

404

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

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

Понятие о передаточной функции нелинейного элемента по­ зволяет представить структурную схему в виде последователь­ ного соединения W H(jA) и W (jus) (рис. 8.29).

Р и с. 8.29. К пояснению метода гармо­ нической. линеаризации. Представление нелинейной системы в виде последова­ тельного соединения Wn{jA) и W (/шj

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

Ф ( і А М)

W » U A ) W U * )

W ( j A , * )

(8.49)

 

\ + W B ( j A ) W ( M

. 1 + IF ( / Л , со)

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

 

I + WUA, Ш) = 0,

 

где

W {j А, со) = W H(JA) W{Ja>) — передаточная

функция ра­

 

зомкнутого контура нелиней­

В

ной системы.

амплитуды

общем случае W (j A , со) является функцией

входных колебаний и частоты. При постоянном значении А и изменении со от 0 до оо вектор W (/А , ш) вычертит на комплекс­ ной плоскости годограф, аналогичный амплитудно-фазовой ха­ рактеристике разомкнутой системы W (/со).

Семейство годографов W{jA, со)

некоторой системы для раз­

личных значений А

приведено на

рис. 8.30.

Если

годограф

W{jA,<ü) при .всех А

не охватывает

точку —

1,/0,

то система

устойчива, и автоколебаний в ней не возникает. В этом случае

при любых А модуль вектора

| W{jA, < o j | < 1 .

Здесь

ча­

стота,

при

которой со (со, Л) =

— іг. Физически

это

будет

озна­

чать,

что

если в результате действия каких-либо

возмущений

в системе возникнут колебания, то после снятия этих возму­

щений колебания будут затухать.

| W [JA, <о_)| > 1,

Если при некоторых значениях А на частоте

то система неустойчива, колебания в ней

будут нарастать.

405

II, наконец, если j W (jA, <uj | = I, система находится па гра­ нице устойчивости, в системе будут автоколебания. Это проис­ ходит при определенных значениях А и ш, которые и являют­ ся искомыми амплитудой и частотой автоколебаний Аа и ша.

Следует заметить, что при однозначной характеристике не­ линейного элемента равенство t f ( с о ) г = — л может быть исполь­ зовано для определения возможной частоты автоколебаний, а

второе

соотношение! W ( jA , ш_) | = 1

позволяет определить ам­

плитуду. Вследствие того, что вектор

W (JA, со)

имеет

зависи­

 

 

мость

не только

от частоты,

но

и от

 

 

амплитуды,

может иметься

несколько

 

 

значений А и со,

при

которых

годо­

 

 

граф

W JA , <о)

проходит

через

точ­

 

 

ку — 1,/0.

 

 

 

 

 

 

 

 

Это значит, что в системе возмож­

 

 

но несколько периодических

режимов

 

 

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

 

 

стотами и амплитудами. Часть из этих

 

 

■режимов действительно может наблю­

 

 

даться

в системе, а

другая

часть не

 

 

может. Наблюдаться могут только ус­

 

 

тойчивые периодические режимы

или

 

 

автоколебания. Другая часть периоди­

 

 

ческих

режимов

является

неустойчи­

 

 

вой и наблюдаться в реальной систе­

Рис .

8.30. Годографы

ме не

может. Возникает таким обра­

зом задача

выделить из всех периоди­

W [jA, ш)для трех значений

 

 

ческих

решений

уравнения

движения

системы устойчивые решения, соответствующие устойчивым пе­

риодическим движениям

или автоколебаниям

системы.

При­

ближенное рассмотрение этого вопроса будет

показано

ниже,

на примерах

 

 

 

 

 

 

 

сто

Для определения амплитуды и частоты автоколебаний вме­

построения

годографа

W (jA,

ш)

обычно

используют

дру­

гой

графический

прием,

связанный

с

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

 

 

1 +

W„ (JA) W (у'ш)

= 0.

 

(8.5Ü)

Это уравнение можно представить в виде

 

 

 

 

— Ен JA) = W (jui),

 

(8.51)

где

E„(jA)==---- ------

—обратная

передаточная функция не-

 

WJJA)

 

линейного элемента.

 

 

Если на комплексной плоскости построить амплитудно-фа­ зовую характеристику линейной части системы W (у'ш) и об­ ратную отрицательную передаточную функцию нелинейного элемента — Ен(/А), то пересечение характеристик будет ука­

4 0 6

зывать на наличие в системе автоколебаний. Пересечение ха­ рактеристик соответствует нахождению системы на границе ус­ тойчивости или прохождению годографа W ( j A , < s > ) через точку

- і,уо.

пересечения характеристик дадут искомые значения

Точки

ша и Л„.

Частота автоколебаний при этом отсчитывается по

W (/со ),

а амплитуда — по характеристике — Е н (/ А).

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

ция есть отношение

А

 

1

(8.52)

W j j A )

У (JA )

 

Поскольку W„(/'A) и En(jA) являются функциями комплекс­ ного аргумента, их можно представить в алгебраической и по­ казательной формах:

W„ (JA) = U„ (А) + jV„ (А) =

W„ (А) ем*);

Ен(./ А ) == S„{A) +

jQ H(А) =

Ен (А)

.

Рассмотрим передаточные

функции

некоторых

нелинейных

звеньев.

З в е н о с и д е а л ь н о й р е л е й н о й х а р а к т е р и с т и ­

кой (см. рис. 8.2,а).

 

 

 

 

 

 

ф ( 4 ) = — ;

W„VA)=W„{A) = - ^ ;

? И ) =

0;

ТС

 

 

 

 

А

 

 

 

 

 

 

 

 

(8.53)

 

 

 

 

 

 

 

Амплитудные зависимости

W H(Л)

и

Ен (А) приведены на

рис. 8.31. График

Ен(Л)

построен

на

комплексной плоскости.

 

 

 

 

 

 

/®н

 

 

 

 

 

ш

й=0

 

 

 

 

 

 

flz

/?7

 

 

 

 

 

 

 

V

 

Рис. 8.31.

Идеальный релейный элемент:

и

а — графики амплитудной характеристики WH(A)

обратной

амплитудной

характеристики

Ен (А)\ б

годограф обратной отрицательной передаточной функ­ ции нелинейного элемента £„ (]А)

407

Однако он построен со знаком минус ( — EH(JA)), что’бы удоб­ нее было его использовать при расчетах, связанных с опреде­

лением амплитуды

и частоты

автоколебаний.

Как

видно из

рис. 8.31,6,

годограф

— £ и (/А)

расположен

вдоль

отрица­

тельной вещественной оси, так

как — En {jA)

в данном случае

не содержит мнимой части.

 

 

 

 

З в е н о

с

з о н о й

н е ч у в с т в и т е л ь н о с т и

(см. рис.

8.2,6).

 

 

 

 

 

 

 

 

Ф(А) = ~

V A 2 -

а2; W H(JA)= W H{A) =

У Ж ^ Г ;

 

■кА

 

 

 

 

кА-1

 

Ф(А) = 0;

 

 

 

- л 2

I

(8.54)

ЕкUА) — £■„ (А) = -—------ = = = = ..

п '

 

н

 

4

’ 4В У А2 — а2

 

На рис. 8.32 приведены амплитудные характеристики и годо­ граф — Е„ (JA ). Условно две сливающиеся ветви показаны раздельно. В действительности они обе располагаются на отри­ цательной вещественной оси.

а

 

Рис. 8.32. Релейный элемент

с зоной

нечувствительности:

— амплитудная

характеристика

WH(Л);

б — обратная амплитуд­

ная

характеристика £„ (Л);

в — годограф

обратной

отрицательной

 

 

 

передаточной функции нелинейного элемента ЕИ(/А)

 

Н е л и н е й н ы й э л е м е н т с н е о д н о з н а ч н о й х а р а к ­

т е р и с т и к о й

(см. рис. 8.3,6).

 

 

 

 

 

Если на вход такого элемента подать синусоидальный сиг­

нал

(рис.

8.33), то

выходная величина

будет серией

знакопе­

ременных

прямоугольных

импульсов (рис.

8.33,а). Изменение

знака

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

будет

опаздывать

на

угол

<р(А) =

 

 

 

а

 

 

 

 

 

 

входной вели­

= arc sin — по отношению к изменению знака

чины.

 

А

 

 

 

 

 

 

 

 

С учетом этого сдвига

 

 

 

 

 

 

 

 

 

 

 

Ф(У'Л)= —

 

 

 

 

(8.55)

 

 

 

 

 

 

TZ

 

 

 

 

 

где

 

<р (Al =

arcsin ——

фазовый сдвиг.

 

 

 

 

 

 

 

 

А

 

 

 

 

 

 

408

Соседние файлы в папке книги из ГПНТБ