Файл: Основы автоматического управления..pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 09.04.2024

Просмотров: 414

Скачиваний: 15

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

§ 10.5. М ЕТОД Т О Ч Е Ч Н Ы Х О Т О БРА Ж ЕН И Й

435

= zjkl оси абсцисс вертикальную прямую до пересечения с кривой F0 (£); из полученной точки пересечения проводим горизонтальную прямую до

пересечения с кривой t \

(Ѳ); абсцисса 04

этой точки пересечения, умноженная

на kl, будет значением щ функции /

(z) при

z = zt.

Построив кривую

V = / (z) (рис.

10.5.4),

видим, что рассматриваемая

система имеет один устойчивый предельный цикл с амплитудой колебаний скорости, равной z* = kit,*.

Заметим, что в случае неявного уравнения кривой точечного отображе­ ния типа Fi (Ѳ) = F0 (Q нет необходимости строить эту кривую, а доста­ точно построить графики левой и правой частей определяющего ее уравне­ ния (рис. 10.5.3). Точки пересечения кривых F0 (Q и Ft (Ѳ) определяют пре­ дельные циклы, поскольку в этих точках ѵ = z. Предельный цикл устойчив,

если в точке пересечения кривых F0 (£) и F, (Ѳ) кривая F0 (Q имеет мень­ ший по абсолютной величине угол наклона к оси абсцисс, и неустойчив,

если кривая F0 (Ö имеет больший по абсолютной величине угол наклона к оси абсцисс.

Определив амплитуду z* колебаний скорости выходной переменной системы, можно найти и все другие параметры автоколебаний — амплитуду колебаний выходной переменной у* и период Т0 (а следовательно, и частоту).

Для этого положим в формуле (10.4.23) t = Т012, z0 =

—z*, z = z*, v = І.

Тогда получим

 

 

Г»

Го

 

z* = z*e 21' + Ы ( 1 - е

2 T ).

(10.5.3)

Отсюда находим период автоколебаний:

 

 

Т0 = 2Т In kl-\-z*

 

(10.5.4)

kl — z*

 

 

Частота автоколебаний равна coo — 2я/Т0. Амплитуду у* колебаний выходной переменной можно определить по формуле (10.4.19), положив в ней ѵ = 1,

У о = —d, z0 = —z*, у

= — у * ,

2 = 0. Тогда получим

 

j* = d-f Г [ Ѵ - Ш п ( і +

- ^ - ) ] .

(10.5.5)

Полагая, например, d =

0,1; Т =

0,1; к = 1;

I = 1, вычисляем

величину Д

по последней формуле

(10.5.1) и

строим графики функций F0

(Ö и Fi (Ѳ)

28*


436 ГЛ . 10. УСТО Й ЧИ ВО СТЬ Н Е Л И Н Е Й Н Ы Х СИСТЕМ

(рис. 10.5.5) в соответствии с формулами (10.5.2). Точка их пересечения дает амплитуду автоколебаний скорости z* = kit* « 0,96. После этого по форму­ лам (10.5.4) определяем период и частоту автоколебаний: Т0 « 0,36 с, со0 «

« 17,5

с-1.

Наконец, по формуле (10.5.5) определяем амплитуду автоколе­

баний:

у*

0,167.

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

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

Для исследования пери­ одических режимов в ре­ лейных системах, содержащих один релейный элемент и произволь­

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

В связи с трудностью точного исследования нелинейных систем большое значение имеют приближенные методы. Пожалуй, наибо­ лее эффективным и вследствие этого наиболее распространенным приближенным методом исследования нелинейных систем является метод гармонической линеаризации. Для систем, близких к линей­ ным, часто применяется так называемый метод малого параметра, основанный на разложении неизвестных элементов движения системы в степенные ряды относительно малого параметра [35]. Этот метод дает хорошие результаты в случае, когда совокупность нелинейных членов в дифференциальных уравнениях системы


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

437

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

Мы ограничимся здесь изложением метода гармонической линеаризации.

§ 10.6. Приближенное исследование нелинейных систем методом гармонической линеаризации

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

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

системы в виде сумм неизвестных

постоянных

составляющих

и синусоидальных составляющих с

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

амплитудами

и начальными фазами и неизвестной общей для всех переменных частотой. При'этом начальную фазу одной из переменных (обычно интересующей нас выходной переменной системы) можно считать равной нулю. Тогда каждая из входящих в уравнения системы

переменных будет зависеть от трех

неизвестных

параметров:

одна — от постоянной составляющей,

амплитуды

и частоты,

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


438

ГЛ . 10. УСТОЙЧИВОСТЬ Н Е Л И Н Е Й Н Ы Х СИСТЕМ

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

Покажем детально, как применяется изложенный общий метод

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

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

Предполагая в нелинейной системе смещенные автоколебания при х = 0 , пред­ ставим все переменные в виде суммы постоянных составляющих и синусоидаль­

ных слагаемых с неизвестными амплиту­ дой, частотой и начальной фазой одного из них:

У =

У* +

ау sin

z =

z* + az sin (ю0^ + Фг)-

(10.6.1)

Осуществляя

гармоническую

линеаризацию нелинейного звена

и пользуясь формулами (8.4.8), (8.4.9), (8.4.12),

представим пере­

менную z

в виде

 

 

 

 

 

Z =

ф* +

I Фн (ау,

У*) I ау sin (со0* +

Фн)>

(Ю.6.2)

где Фн (ßy, у*) — комплексный гармонический коэффициент уси­ ления, а ф„ = arg Фн (а.у, У*)- Выходную переменную у линей­ ной части системы в силу принципа суперпозиции выразим через переменную z в следующем виде:

у = —Ф (0) z* — |Ф(гсо0) I аі sin

+ ф2 +

Ф (too)].

(Ю.6.3)

Приравнивая правые части в формулах (10.6.1)

и (10.6.2) для z

и в формулах (10.6.1) и (10.6.3) для г/, получим

 

 

z* + az sin (cö0* + ф2)=Ф* + І Фн К ,

У*) I а-ѵsin (а>0*+Фн)>

(Ю.6.4)

У * + ау sin ( ä 0 t = —Ф (0) Z* — I Ф (icöo) I a z sin [®о^+Фг+Ф(®о)] • (10.6.5)

Сравнивая в правых и левых частях уравнений (10.6.4) и (10.6.5) постоянные составляющие, а также амплитуды и фазы синусоидаль­ ных составляющих, получим уравнения

Ф*(ay,y*) = z*,

Ф(0) z * = - y \

(10.6.6)

I Фц iflyi У*) I Q'V=

I Ф (г®о) I O-Z—dyi

фн (a,j, у*) = ф г, ф г + ф Ы = Я.


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

439

Исключая из уравнений (10.6.6) переменные az, z*, ф2, получим три

уравнения для определения переменных у*, ау,

со0:

У* =

—Ф (0) ф* (ау,

у*),

 

(10.6.7)

|Ф(іюо)

М Ф н К , У*)

1=

1)

(Ю.6.8)

arg Ф (ісоо) + arg Фн (ау,

у*)

= я.

(10.6.9)

Решить систему уравнений (10.6.7), (10.6.8) и (10.6.9) в общем виде невозможно. Для решения этих уравнений можно применить различные приближенные методы и, в частности, метод последова­ тельных приближений. Удобно также воспользоваться графи­ ческими способами. Рассмотрим воз­ можные способы графического реше­ ния этой системы уравнений.

Первый способ состоит в следую­ щем. Задаваясь значениями г/*, нахо­ дим из уравнения (10.6.7) соответ­ ствующие значения ау. По результа­ там этих вычислений строим на плоскости параметров у*, ау кривую 1 (рис. 10.6.2). Затем при определен­ ном значении со„ задаемся значениями у* и определяем из уравнений (10.6.8)

и (10.6.9) соответствующие значения ау. По результатам этих вычислений строим кривую 2 в координатах у*, ау, соответствую­ щую уравнению (10.6.8), и кривую 3, соответствующую уравне­ нию (10.6.9). Точка пересечения кривых 2 и 3 является решением уравнений (10.6.8) и (10.6.9) для заданного значения со0. Построив кривые 2 и 3 для ряда значений со0, соединим точки их пересечения кривой 4. Точка пересечения кривой 4 с кривой 1 определяет зна­ чения у* и ау, удовлетворяющие всем трем уравнениям (10.6.7), (10.6.8) и (10.6.9). Для определения частоты автоколебаний со0 следует подставить найденные значения у* и ау в уравнение (10.6.8) или в уравнение (10.6.9). В результате получим уравнение для определения ю0. Частоту автоколебаний можно также определить по шкале значений ю0, получаемой на кривой 4 при ее построении по достаточно большому числу точек. Однако этот способ менее точен, чем предыдущий.

Второй способ аналогичен первому, но отличается только при­ емом отыскания кривой 4. Объединим уравнения (10.6.8) и (10.6.9)

в одно уравнение

 

 

Ф (і(о0) Фн («у. У*)

= —1

(10.6.10)

и перепишем его в виде

 

 

1

 

(10.6. 11)

Ф (іюо) =

У*)'

Фн (“у,