
- •Раздел 2. Нелинейные системы.
- •Глава 10.
- •Основные понятия.
- •§1. Ведение в динамику нелинейных систем автоматического управления.
- •Следующие две задачи имеют особое значение для анализа динамики нелинейных систем.
- •§ 2. Статические характеристики типовых нелинейных элементов.
- •Однозначные нелинейности.
- •Неоднозначные нелинейности
- •Нелинейные характеристики, содержащие люфт.
- •Начальные значения при математическом описании неоднозначных нелинейностей.
- •Регуляторы «включить-выключить».
- •Динамика нелинейных систем автоматического управления.
- •Глава 11. Методы фазового пространства.
- •§1. Основные понятия теории устойчивости.
- •Свойства фазовых траекторий.
- •Все корни характеристического уравнения находятся в левой полуплоскости (имеют отрицательную действительную часть). В этом случае положение равновесия асимптотически устойчиво.
- •Хотя бы один корень характеристического уравнения находится в правой полуплоскости (имеет положительную действительную часть). В этом случае положение равновесия неустойчиво.
- •§2. Фазовая плоскость.
- •Выясним, как отображаются на фазовой плоскости различные функции времени.
- •2. Вычисление интервала времени при движении по фазовой траектории.
- •3. Рассмотрим типы фазовых траекторий для линейного уравнения второго порядка.
- •4. Нелинейное уравнение второго порядка.
- •§3. Исследование динамики релейных систем на многолистной фазовой плоскости.
- •1. Постановка задачи.
- •2. Общие вопросы изучения динамики релейной системы на фазовой плоскости. Многолистная фазовая плоскость.
- •2.3. Запаздывание.
- •3. Приведение дифференциального уравнения системы к простейшей форме.
- •3. Коэффициент воздействия по производной f заменяется величиной: (11)
- •4. Исследование динамики релейной системы управления объектом, имеющим передаточную функцию
- •§4. Система управления с переменной структурой объектом второго порядка.
- •§5. Исследование динамики релейной системы управления объектом, имеющим передаточную функцию .
- •1. Упрощенные принципиальные схемы двух систем данного класса.
- •2.1. Форма фазовых траекторий релейных систем данного класса.
- •2.2. Запаздывание.
- •3. Исследование динамики системы.
- •§6. Точечные преобразования на фазовой плоскости и их применение для расчета релейных систем
- •Рассмотрим вычисление параметров автоколебаний с помощью точечного преобразования.
- •Глава 12.
- •§1. Структурное преобразование нелинейных систем. Приведение структурной схемы к канонической форме.
- •§2. Частотные критерии абсолютной устойчивости положения равновесия и отрезка покоя.
- •§3. Автоколебания в системах управления.
- •1.Гармонический коэффициент линеаризации нелинейного элемента.
- •2.Уравнение гармонического баланса.
- •§4. Численный эксперимент. Опасные явления в нелинейных системах.
§3. Автоколебания в системах управления.
Постановка задачи.
Предположим, что исходная система, содержащая один нелинейный элемент, приведена к виду, показанному на рисунке 1. Входное воздействие будем предполагать импульсным. То есть система выводится из положения равновесия (в некоторую точку фазового пространства), после чего входное воздействие равно нулю.
Рис. 1.
Таким образом, мы будем рассматривать собственные движения системы. Нас будет интересовать, имеет ли место в системе автоколебания, и каковы их параметры. С математической точки зрения задача состоит в том, чтобы определить имеются ли предельные циклы Пуанкаре в фазовом пространстве, каковы их параметры. Решить вопрос об устойчивости предельных циклов (после того, как они найдены). Нужно отметить, что поставленная задача значительно труднее, чем задача об исследовании положения равновесия, и её решение в общем виде на сегодняшний день отсутствует.
Академиком А.Андроновым и его учениками (Горьковская школа по теории колебаний) был разработан в середине сороковых годов метод точечных преобразований и с его помощью были найдены точные решения классических задач авторегулирования второго порядка и в некоторых случаях третьего порядка. Отметим, что нелинейные элементы в этих задачах были кусочно-линейными. Дальнейшее продвижение в этом направлении, т.е использование точечного метода для преобразования систем более высокого порядка, связано со значительными трудностями.
Поэтому ниже мы приведем только один результат точного решения задачи изучения автоколебаний для систем сколь угодно высокого порядка, полученный В.А.Якубовичем.
Частотный критерий автоколебательности В.А.Якубовича.
Предварительно нам потребуется формулировка автоколебательной системы. Подчеркнём, что колебания могут быть сколь угодно сложной формы!
Система называется автоколебательной, если удовлетворяются следующие условия:
1.Сигнал на выходе y(t) системы ограничен.
2.При t→∞ имеется бесконечно много перемен знака сигнала y(t).
3.Пусть tk’ и tk’’ – бесконечные последовательности моментов времени.
Имеют место соотношения:
(1)
4.Время пребывания сигнала y(t) в промежутках (-∞;α) , (α;β), (β;∞) ограничено величиной Т.
Обратим внимание, что числа α и β являются нижними границами автоколебаний.
Система будет автоколебательной, если ее нелинейный элемент имеет график, показанный на рисунке 2.
-
Рис. 2
На рисунке 2 сектор II - это сектор неустойчивости линейной системы сравнения (ЛСС), а сектор I – сектор абсолютной устойчивости исходной системы.
Метод гармонического баланса для исследования автоколебаний, имеющих форму, близкую к синусоиде.
В сороковых годах Л.С.Гольдфарбом был предложен удивительно простой метод [x] определения автоколебаний в системах автоматического управления. Метод не является математически доказанным утверждением и базируется на правдоподобных рассуждениях. Этот метод позднее был с исключительной полнотой развит в работах Е.П.Попова [x]. Сразу же отметим очень важное обстоятельство. Метод не дает информации о жестком или мягком режиме возникновения автоколебаний.
Перейдём к его описанию. Л.С.Гольдфарб обратил внимание, что очень часто автоколебания в технических системах авторегулирования на выходе объекта имеют почти гармонический характер. Запись колебаний выходной координаты y(t) может быть аппроксимировано функцией Asin(Ωt), если соответственно подобрать числа А и Ω. На основе этого наблюдения Л.С.Гольдфарб построил логику своих рассуждений на следующем исходном предположении.
Сигнал на входе нелинейного элемента (x(t)=-y(t)) следует принять равным
x(t)=Asin(Ωt) (2)
где числа А и Ω необходимо вычислить.
Из исходного предположения вытекает, что автоколебания в системе имеют период:
(3)
Сигнал на выходе нелинейного элемента u(t)=N(Asin(Ωt)) является периодической функцией с периодом Т, но форма его колебаний может быть очень далека от гармонической, например, прямоугольная волна и т.д. сигнал u(t) представляет собой ряд Фурье:
(4)
(Для простоты предположим, что постоянная составляющая отсутствует, т.е. D0=0.)
Далее по логике рассуждений необходимо, чтобы линейный блок пропускал только первую гармонику сигнала u(t), т.е. АЧХ W(p) имело форму фильтра высоких частот.
-
Рис. 3
Физически такая АЧХ невозможна у линейного блока, однако допустима ситуация, когда АЧХ имеет форму, показанную на рисунке 4. В этом случае говорят, что линейный блок удовлетворяет «гипотезе фильтра».
Возможна ситуация, когда АЧХ линейного блока удовлетворяет гипотезе резонанса и имеет вид, показанный на рисунке 5, но этот случай часто встречается в теории колебаний и нехарактерен для задач авторегулирования. Далее будем предполагать, что АЧХ удовлетворяет гипотезе фильтра.
Рис. 4а |
Рис. 4b |
На рисунке 4а приведена амплитудно-частотная характеристика в случае, когда W(p) не содержит интегрирующее звено, на рисунке 4b – когда содержит.
-
Рис. 5
Автоколебания – это стационарный процесс, поэтому при его рассмотрении мы предполагаем, что он длится уже достаточно долго, и переходный процесс линейного блока давно затух. Следовательно, на выходе линейного блока мы можем рассматривать только установившийся процесс, который будет равен
(5)
Предполагаем, что для k > 1 величина |W(Ωk)| пренебрежимо мала, что соответствует «гипотезе фильтра».
Если удовлетворяется «гипотеза фильтра», то приближённо можно записать
(6)
Отметим, что идеальный случай, когда |W(Ωk)|=0 при k>1, принципиально невозможен.
Подчеркнем, что в формуле (6) мы отбросили высшие гармоники не потому, что они малы, а потому, что их влияние на выходной сигнал системы пренебрежимо мало благодаря «гипотезе фильтра».
Таким образом, задача состоит в том, чтобы зная N(x) и W(p), вычислить числа А и Ω. Для решения этой задачи предварительно введем два важных понятия.