Метод гармонической линеаризации: Методические указания к лабораторной работе. Гармоническая линеаризация Коэффициенты гармонической линеаризации

14.08.2021

Министерство образования и науки Российской Федерации

Саратовский государственный технический университет

Балаковский институт техники, технологии и управления

Метод гармонической линеаризации

Методические указания к лабораторной работе по курсу «Теория автоматического управления» для студентов специальности 210100

Одобрено

редакционно –издательским советом

Балаковского интститута техники,

технологии и управления

Балаково 2004

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

ОСНОВНЫЕ СВЕДЕНИЯ

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

Предполагается, что исследуемая нелинейная САУ может быть представлена в следующем виде

причем нелинейная часть должна иметь одну нелинейность

. (1)

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

Любую функцию или сигнал можно разложить в ряд по системе линейно-независимых, в частном случае ортонормированных функций. В качестве такого ортогонального ряда может быть использован ряд Фурье.

Разложим в ряд Фурье выходной сигнал нелинейной части системы

, (2)

здесь - коэффициенты Фурье,

. (3)

Таким образом, сигнал согласно (2) может быть представлен в виде бесконечной суммы гармоник с возрастающими частотами и т. д. Этот сигнал поступает на вход линейной части нелинейной системы.

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

, (4)

причем степень полинома числителя должна быть меньше степени полинома знаменателя. В этом случае АЧХ линейной части имеет вид

где 1 - не имеет полюсов, 2 - имеет полюс или полюса.

Для АЧХ справедливо записать

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

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

Этот сигнал проходит по всему замкнутому контуру системы Рис.1 и на выходе нелинейного элемента без учета более высоких гармоник, согласно (2) имеем

. (7)

При исследовании нелинейных систем с помощью метода гармонической линеаризации возможны случаи симметричных и несимметричных колебаний. Рассмотрим случай симметричных колебаний. Здесь и.

Введем следующие обозначения

Подставив их в (7), получим . (8)

С учетом того, что

, где ,

. (9)

Согласно (3) и (8) при

. (10)

Выражение (9) является гармонической линеаризацией нелинейности устанавливает линейную связь входной переменной и выходной при . Величины и называются коэффициентами гармонической линеаризации.

Необходимо отметить, что уравнение (9) является линейным для конкретных величин и (амплитуды и частоты гармонических колебаний в системе). Но в целом оно сохраняет нелинейные свойства, так как коэффициенты различны для различных и . Эта особенность и позволяет исследовать с помощью метода гармонической линеаризации свойства нелинейных систем [ Попов Е.П.].

В случае несимметричных колебаний гармоническая линеаризация нелинейности приводит к линейному уравнению

. (12)

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

Уравнения (9) и (11) позволяют получить передаточные функции гармонически линеаризованных нелинейных элементов. Так для симметричных колебаний

В этой главе будет изложен метод гармонической линеаризации для приближенного определения периодических решений (автоколебаний) и устойчивости нелинейных систем любого порядка, который по идее близок к методу эквивалентной линеаризации или методу гармонического баланса Н. М. Крылова и Н. Н. Боголюбова, а по результатам - также и к методу малого параметра Б. В. Булгакова.

Рассматриваемый приближенный метод является мощным средством исследования нелинейных автоматических систем в смысле простоты и довольно большой универсальности его аппарата в применении к самым разнообразным нелинейностям. Однако надо иметь в виду, что он решает задачу приближенно. Имеются определенные ограничения его применимости, о которых будет сказано ниже. Эти ограничения обычно хорошо соблюдаются в задачах теории автоматического регулирования. Практические расчеты и эксперимент показывают приемлемость этого метода для многих видов нелинейных систем.

Пусть дано какое-нибудь нелинейное выражение вида

Разложив функцию в правой части выражения (18.1) в ряд Фурье, получим

что означает отсутствие постоянной составляющей в данном разложении. В настоящей главе будет везде предполагаться выполнение условия отсутствия постоянной составляющей (18.5). Впоследствии (глава 19) будет дан метод исследования автоколебаний при наличии постоянной составляющей, т. е. в случае невыполнения условия (18.5).

Если принять во внимание, что из (18.2) и (18.3)

то формулу (18.4) при условии (18.5) можно будет записать в виде

где q - коэффициенты гармонической линеаризации, определяемые формулами:

Итак, нелинейное выражение (18.1) при заменяется выражением (18.6), которое с точностью до высших гармоник аналогично линейному. Эта операция и называется гармонической линеаризацией. Коэффициенты постоянны при постоянных значениях , т. е. в случае периодического процесса. В переходном колебательном процессе с изменением а и со коэффициенты q и изменяются (см.гл. 20). Для разных амплитуд и частот периодических процессов коэффициенты выражения (18.6) будут различны по величине. Это очень важное для дальнейшего обстоятельство является существенным отличием гармонической линеаризации, по сравнению с обычным способом линеаризации (§ 3.1), приводящим к чисто линейным выражениям, которые применялись в предыдущих разделах книги. Указанное обстоятельство позволит путем применения к выражению (18.6) линейных методов исследования проанализировать основные свойства нелинейных систем, которые не могут быть обнаружены при обычной линеаризации.

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

Здесь возможны два варианта: 1) кривая имеет гистерезисную петлю (например, рис. 16.18, в, рис. 16.22, г, д), и 2) кривая не имеет гистерезисной петли (рис. 16.8, б, рис. 16.22, а и др.).

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

при условии отсутствия постоянной составляющей:

Если же кривая не имеет гистерезисной петли, то так как при будет

(при гистерезисной петле этот интеграл не был нулем вследствие различия в очертании кривой при возрастании и убывании

Следовательно, при отсутствии гистерезисной петли нелинейное выражение (18.8) заменяется более простым:

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

Гистерезисная же петля вводит согласно (18.9), кроме того, еще производную, дающую отставание по фазе, так как Таким образом, нелинейное отставание по координате в виде гистерезисной петли превращается при гармонической линеаризации в эквивалентное линейное отставание по фазе.

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

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

Например, при гармонической линеаризации второго из уравнений (16.3) придется иметь дело с функцией при . В этом случае получаем

при условии

Если же функция или функция будет единственной нелинейной функцией в уравнении нелинейного звена, то при гармонической

линеаризации можно положить и

аналогично прежним формулам (18.6) и (18.7). Но при этом величина а во всех выкладках будет амплитудой колебаний скорости а не самой координаты х. Последняя же будет иметь тогда амплитуду

При вычислении коэффициентов гармонической линеаризации по формулам (18.10) надо иметь в виду, что при симметричных нелинейных характеристиках интеграл можно получить удвоением интеграла , т. е.

а для симметричных относительно начала координат безгистерезисных характеристик при вычислении можно писать

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

Коэффициенты гармонической линеаризации релейных звеньев. Найдем коэффициенты и уравнений наиболее типичных релейных звеньев по формулам (18.10). Возьмем общий вид характеристики релейного звена изображаемой графиком рис. 18.1, а, где есть любое дробное число в интервале

Как частные случаи будут получены уравнения других типов релейных звеньев.

Если колебания входной величины имеют амплитуду то согласно рис. 18.1, а движения в системе не будет. Если амплитуда то переключения реле происходят в точках А, В, С, D (рис. 18.1, б), в которых имеем

Следовательно, после использования свойств каждый из интегралов (18.10) разбивается на три слагаемых:

причем первое и третье из них согласно рис. 18.1, а и будут нулями. Поэтому выражения (18.10) принимают вид

а уравнение релейного звена с характеристикой вида рис. 18.1, а будет иметь вид (18.9) с полученными здесь значениями и .

Рассмотрим частные случаи.

Для релейного звена с характеристикой без гистерезисной петли, но с зоной нечувствительности (рис. 18.1, а), полагая из вышенаписанных формул получаем

Для релейной характеристики с гистерезисной петлей типа рис. полагая имеем

Наконец, для идеального релейного звена (рис. 18.1, е), полагая находим

На последнем примере легко видеть смысл гармонической линеаризации релейной характеристики. Написанное выражение для q означает замену ломаной характеристики прямолинейной (рис. 18.1, е) с таким наклоном, чтобы эта прямая приблизительно заменяла собой тот участок ломаной который охватывается заданной амплитудой а. Отсюда становится вполне понятной обратно пропорциональная зависимость от а, даваемая формулой (18.18), так как чем больше амплитуда а колебаний входной величины тем более пологой должна быть прямая приблизительно заменяющая ломаную

Аналогично обстоит дело и с релейной характеристикой на рис. 18.1, г для которой наклон заменяющей ее прямой дается формулой (18.16). Следовательно, всякое безгистерезисное релейное звено в колебательном процессе эквивалентно такому «линейному» звену, передаточное число (коэффициент усиления) которого уменьшается с увеличением амплитуды колебаний входной величины, начиная с

Что касается релейного звена с гистерезисной петлей, то согласно (18.9) и (18.17) оно заменяется линейным звеном с аналогичным прежнему коэффициентом усиления , но, кроме того, еще с введением отрицательной производной в правой части уравнения. Введение отрицательной производной в противовес положительной (см. § 10.2) вносит отставание по фазе в реакции звена на входное воздействие. Это служит «линейным эквивалентом», заменяющим эффект действия нелинейности в виде гистерезисной петли. При этом

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

Коэффициенты гармонической линеаризации других простейших нелинейных звеньев. Рассмотрим нелинейное звено с зоной нечувствительности и с насыщением (рис. 18.2, а). Согласно рис. 18.2, б, где

интеграл (18.10) на участке разбивается на пять слагаемых, причем два из них равны нулю. Поэтому

откуда с заменой получаем

где определяются формулами (18.19). Ввиду отсутствия гистерезисной петли здесь

Итак, уравнение нелинейного звена с характеристикой вида рис. 18.2, а будет где определяется выражением (18.20).

Как частный случай отсюда получается значение для звена с зоной нечувствительности без насыщения (рис. 18.2, в). Для этого в предыдущем решении нужно положить и, следовательно, Тогда

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

Введение

Релейные системы получили широкое распространение в практике автоматического регулирования. Достоинством релейных систем является простота конструкции, надежность, простота обслуживания и настройки. Релейные системы представляют собой особый класс нелинейных АСР.

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

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

Структура исследуемой системы

Исследуемая система (рис.) 1, включает в себя элемент сравнения ЭС, релейный элемент РЭ, исполнительный двигатель (идеальный интегратор с коэффициентом усиления =1), объект регулирования (апериодическое звено с тремя постоянными времени , , и коэффициентом усиления ). Значения параметров системы приведены в табл. 1 приложения А.

Статические характеристики (характеристики вход-выход) исследуемых релейных элементов приведены на рис. 2.

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

Исследуемая АСР может быть смоделирована с помощью известных пакетов моделирования, например, SIAM или VisSim.

Замечание. В некоторых пакетах моделирования значение выходного

сигнала реле может принимать лишь значения ±1 вместо ±В, где В произвольное число. В таких случаях необходимо коэффициент усиления интегратора принять равным .

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

Для выполнения работы каждый студент получает от преподавателя вариант исходных данных (см. раздел 2).

Работа выполняется в два этапа.

Первый этап – расчетно-исследовательский (может быть выполнен вне лаборатории).

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

Необходимый теоретический материал изложен в разделе 4; в разделе 5 приведены контрольные вопросы.

3.1. Расчетно – исследовательская часть

1. Получить выражения для амплитудно-частотной и фазо-частотной, вещественной и мнимой характеристик линейной части системы.

2. Рассчитать и построить амплитудно-фазовую характеристику линейной части системы. Для расчета использовать программы из пакета ТАУ. Обязательно распечатать значения вещественной и мнимой частотной характеристик (10 – 15 точек, соответствующих третьему и второму квадрантам).

4. Используя графоаналитический метод Гольдфарба, определить амплитуду и частоту автоколебаний и их устойчивость для всех четырех реле. Расчет параметров автоколебаний можно осуществить и аналитически. Качественно изобразить фазовый портрет системы для каждого из случаев.

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

Экспериментальная часть

1. Используя один из доступных пакетов моделирования, собрать схему моделирования исследуемой АСР. По разрешению преподавателя можно воспользоваться готовой схемой. Настроить параметры схемы в соответствии с заданием.

2. Исследовать переходный процесс в системе с идеальным реле (вывести на печать), подавая на вход скачкообразное воздействие x(t)=40*1(t). Измерить амплитуду и частоту автоколебаний, сравнив их с расчетными значениями. Повторить эксперимент, установив не нулевые начальные условия (например, у(0)=10, у (1) (0)=-5).

3. Исследовать переходный процесс в системе с трехпозиционным реле для двух различных значений амплитуды входного сигнала x(t)= 40*1(t) и x(t)=15*1(t). Вывести на печать переходные процессы, измерить амплитуду и частоту автоколебаний (если они существуют), сравнить их с расчетными значениями, сделать выводы.

4. Исследовать переходные процессы в системе с трехпозиционным реле для других значений коэффициента усиления линейной части (см. п.5, раздел 3.1).

5. Исследовать переходные процессы в системе с двухпозиционными реле с гистерезисом при нулевых и не нулевых начальных условиях и x(t)=40*1(t). Вывести на печать переходные процессы, измерить амплитуду и частоту автоколебаний (если они существуют), сравнить их с расчетными значениями, сделать выводы.

Теоретическая часть

Широко распространенным методом расчета нелинейных систем является метод гармонической линеаризации (описывающих функций) .

Метод позволяет определять параметры автоколебаний (амплитуду и частоту), устойчивость автоколебаний, устойчивость положения равновесия нелинейной АСР. На базе метода гармонической линеаризации разработаны методы построения переходных процессов, анализа и синтеза нелинейных АСР .

Метод гармонической линеаризации

Как уже отмечалось, в нелинейных и в особенности релейных АСР часто наблюдаются устойчивые периодические колебания постоянной амплитуды и частоты, так называемые автоколебания . Причем автоколебания могут сохраняться даже при значительных изменениях параметров системы. Практика показала, что во многих случаях колебания регулируемой величины (рис. 3) близки к гармоническим.

Близость автоколебаний к гармоническим позволяет использовать для определения их параметров – амплитуды A и частоты w 0 – метод гармонической линеаризации. В основе метода лежит предположение, что линейная часть системы является фильтром низких частот (гипотеза фильтра). Определим условия, при которых автоколебания в системе могут быть близки к гармоническим. Ограничимся системами, которые как на рис. 3 могут быть приведены к последовательному соединению нелинейного элемента и линейной части. Предположим, что сигнал задания величина постоянная, для простоты примем его равным нулю. А сигнал ошибки (рис 3) является гармоническим:

(1)

Выходной сигнал нелинейного элемента как всякий периодический сигнал – на рисунке 3 это прямоугольные колебания – может быть представлен в виде суммы гармоник ряда Фурье.

Допустим, что линейная часть системы является фильтром низких частот (рис. 4) и пропускает только первую гармонику с частотой w 0 . Вторая с частотой 2w 0 и более высокие гармоники отфильтровываются линейной частью. В этом случае на выходе линейной части будет существовать практически только первая гармоника , а влиянием высших гармоник можно пренебречь

Таким образом, если линейная часть системы является фильтром низких частот, а частота автоколебаний w 0 удовлетворяет условиям

, (4)

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

(5)

не меньше двух

Условие (6) выполняется для многих реальных систем. Примером могут служить апериодическое звено второго порядка и реальное интегрирующее

. (7)

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

где и – действительная и мнимая части ,

– аргумент ,

– модуль .

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

(9)

численно равен отношению амплитуды первой гармоники на выходе нелинейного элемента к амплитуде входного гармонического сигнала.

Аргумент

(10)

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

Для неоднозначных нелинейностей, рис. 2,в, 2,г, определяется по формуле

где S – площадь петли гистерезиса. Площадь S берется со знаком плюс, если петля гистерезиса обходится в положительном направлении (рис. 2,в) и со знаком минус в противном случае (рис. 2,г).

В общем случае и вычисляются по формулам

, (12)

где , – нелинейная функция (характеристика нелинейного элемента).

С учетом вышеизложенного, при исследовании автоколебаний, близких к гармоническим, нелинейная АСР (рис. 3) заменяется эквивалентной с коэффициентом гармонической линеаризации вместо нелинейного элемента (рис. 5). Выходной сигнал нелинейного элемента на рис. 5 обозначен как , это

подчеркивает, что нелинейный элемент генерирует только

первую гармонику колебаний. Формулы для коэффициентов гармонической линеаризации для типовых нелинейностей можно найти в литературе, например, в . В таблице приложения В приведены характеристики исследуемых релейных элементов, формулы для и их годографы. Там же приведены формулы и годографы для обратного коэффициента гармонической линеаризации , определяемого выражением

, (13)

где и действительная и мнимая часть . Годографы и строятся в координатах , и , соответственно.

Запишем теперь условия существования автоколебаний. Система на рис. 5 эквивалентна линейной. В линейной системе существуют незатухающие колебания, если она находится на границе устойчивости. Воспользуемся условием границы устойчивости по критерию Найквиста:

. (14)

Уравнение (14) естьусловие существования автоколебаний, близких к гармоническим. Если существуют действительные положительные решения А и w 0 уравнения (14), то в нелинейной АСР существуют автоколебания близкие к гармоническим. В противном случае автоколебания отсутствуют или не являются гармоническими. Уравнение (14) распадается на два – относительно действительной и мнимой части:

;

Поделив обе части уравнения (14) на и принимая во внимание формулу (13), получим условие существования автоколебаний в форме Гольдфарба Л.С.:

. (17)

Уравнение (17) также распадается на два:

(18)

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

Гольдфарб предложил графоаналитический метод решения системы (17) и определения устойчивости автоколебаний.

В координатах , и , строятся годографы и (рис. 6,а). Если годографы пересекаются, то автоколебания существуют. Параметры автоколебаний – А и w 0 определяются в точках пересечения – частота w 0 по годографу , амплитуда по годографу . На рис. 6,а – две точки пересечения, что указывает на наличие двух предельных циклов.

б)

Для определения устойчивости автоколебаний согласно Гольдфарбу штрихуется левая сторона АФХ линейной части при движении вдоль АФХ в направлении возрастания частоты (рис. 6).

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

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

На рис. 6,б качественно изображен фазовый портрет соответствующий двум предельным циклам на рис. 6,а. Точке пересечения с параметрами и на рис. 6,а соответствует не устойчивый предельный цикл на рис. 6,б, точке с параметрами и и добиться срыва автоколебаний , в этом случае годографы и не пересекаются. Этого же эффекта можно добиться, увеличив зону нечувствительности d или уменьшив амплитуду выходного сигнала реле В. Существует некоторое предельное значение К л, при котором АФХ линейной части касается Ошибка! Ошибка связи. при этом , а значение амплитуды равно . Естественно, это приводит к качественному изменению фазового портрета системы.

, (4)

не меньше двух

где и – действительная и мнимая части ,

– аргумент ,

– модуль .

Аргумент

Для неоднозначных нелинейностей, рис. 2,в, 2,г, определяется по формуле

В общем случае и вычисляются по формулам

где , – нелинейная функция (характеристика нелинейного элемента).

Подчеркивает, что нелинейный элемент генерирует только

где и действительная и мнимая часть . Годографы и строятся в координатах , и , соответственно.

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

Нелинейная система представляется в виде соединения линейной и нелинейной части (рис. 5).

Рис. 5 Схема нелинейной системы

Выходной сигнал нелинейной части системы в общем случае определяется выражением

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

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

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

В этом выражении a i и b i - коэффициенты Фурье. Для симметричных нелинейностей F 0 =0.

Основным условием, которое накладывает метод на линейную часть системы, является условие фильтра нижних частот. Считается, что линейная часть пропускает только первую гармонику колебаний. Данное допущение позволяет считать высшие гармоники в (7.19) несущественными и ограничиться рассмотрением только первой гармоники сигнала y(t).

то выражение (7.20) можно переписать в виде

Первое уравнение системы (7.17) примет вид

В этом выражении

Результат замены нелинейности F(x,sx) выражением

и называется гармонической линеаризацией. Величины q и q 1 называются коэффициентами гармонической линеаризации или просто гармоническими коэффициентами. Для однозначных нелинейностей обычно q 1 =0 . Формулы для гармонических коэффициентов, соответствующих типовым нелинейностям, приводятся в приложениях.

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

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

1). В характеристическом уравнении системы, полученном из (7.22) делаем замену s=j и получим

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

3).Решение этой системы дает частоту и значения гармонических коэффициентов. Если эти значения вещественны и положительны, то в системе существует предельный цикл. По значениям гармонических коэффициентов можно определить амплитуду предельного цикла.
4). Общим признаком устойчивости предельного цикла, т.е. существования автоколебаний, является равенство нулю предпоследнего определителя Гурвица при полученных значениях амплитуды и частоты предельного цикла. Часто более удобно использовать условие устойчивости предельного цикла, в основе которого лежит критерий устойчивости Михайлова.

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

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

Рис.6 Схема нелинейного рулевого привода

Зададим следующие параметры нелинейности скоростной характеристикм рулевого привода: b = 0.12, k 1 = tg =c/b = 6.7. Коэффициенты гармонической линеаризации этой нелинейности определяются выражениями

Заменив в схеме нелинейную характеристику гармоническим коэффициентом, получим передаточную функцию рулевого привода

Подставим эту передаточную функцию в структурную схему системы стабилизации угла тангажа и определим передаточную функцию замкнутой системы

В характеристическом уравнении замкнутой системы сделаем замену s = j и выделим вещественную и мнимую части.

Из второго уравнения системы получим выражение для частоты: , и подставив его в первое уравнение, после преобразований получим

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

По этим значениям можно вычислить для двух случаев все коэффициенты характеристического уравнения и определить частоты, соответствующие каждому значению q(А). Получим:

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

Теперь оценим устойчивость предельных циклов. Используем неравенство, полученное из критерия Михайлова, для чего определим

Производная от коэффициента гармонической линеаризации, входящая в полученные выражения, вычисляется по формуле

Расчеты по выше приведенным формулам показывают, что первый предельный цикл не устойчив и возникает он при (0) 0.1166(6.7 0 ). Если начальное отклонение меньше указанного, то процесс на входе нелинейного элемента затухает (рис.7. 7) и система устойчива.

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

Рис. 8

Путем моделирования определено, что область притяжения устойчивого предельного цикла лежит приблизительно в пределах (0) 0.1167 - 1.4 (6.71 0 - 80.2 0 ).

samarapedsovet.ru