скачать З ![]() Исследовать устойчивость состояния равновесия нелинейной следящей системы.
Исходные данные. ![]() ![]() ![]() ![]() ![]() Введение. Линеаризацию нелинейностей будем проводить, используя метод гармонической линеаризации. Важное преимущество метода гармонической линеаризации заключается в возможности его применения к системам высокого порядка с любой сложностью линейных частей системы и с самыми разнообразными комбинациями мест включения нелинейных звеньев. В более или менее простых нелинейных системах метод гармонической линеаризации позволяет решить задачу полностью аналитически или с применением графо - аналитических построений. В основу метода гармонической линеаризации исследования периодических процессов в нелинейных автоматических системах положено гармоническое представление сигналов системы. Далее будем рассматривать системы автоматического управления с одним нелинейным звеном. Кроме этого полагаем, что линейная часть системы имеет передаточную функцию ![]() Тогда структурную схему системы можно представить в виде, представленном на рисунке 1. ![]() И ![]() Звено с передаточной функцией ![]() ![]() ![]() Вычислим передаточную функцию линейной части. ![]()
Передаточная функция линейной части системы будет выглядеть следующим образом: ![]() По определению передаточная функция ![]() ![]() ![]() ![]() На основании свойств преобразования Лапласа из последнего равенства получаем дифференциальное уравнение относительно регулируемой величины x(t): ![]() Подставляя численные значения, получаем: ![]() 2. Вычисление коэффициентов гармонической линеаризации: Нелинейное звено типа «Ограничение». Его характеристика имеет следующий вид: ![]() где ![]() ![]() ![]() Вычислим коэффициент гармонической линеаризации ![]() ![]() ![]() Учитывая график сигнала на выходе нелинейного звена, получаем: ![]() А так как ![]() ![]() то преобразуя выражение q(a) получим: ![]() А так как характеристика нелинейного звена не имеет гистерезиса, то q(α)=0. В результате имеем: ![]() Подставляя численные значения, получаем: ![]() ![]() Построим график функции q(a): ![]() ^ Пусть теперь структурная схема гармонически линеаризованной системы управления имеет вид, показанный на рисунке 2, и пусть линейная часть системы с передаточной функцией ![]() обладает свойством фильтра. Тогда уравнение гармонически линеаризованной системы запишется в виде ![]() где амплитуда ![]() ![]() ![]() ![]() ![]() ![]() Линеаризованное уравнение будет выглядеть следующим образом ![]() ![]() Подставляя численные значения, получаем: ![]() ^ Уравнение ![]() можно рассматривать как линейное уравнение с постоянными коэффициентами, для которого характеристическое уравнение будет выглядеть следующим образом: ![]() в нашем случае ![]() Обозначим ![]() В результате характеристическое уравнение имеет вид: ![]() Подставляя численные значения, получаем: ![]() ![]() ^ Если в системе возможен периодический процесс ![]() ![]() ![]() ![]() тогда характеристическое уравнение для нашего случая будет выглядеть следующим образом ![]() ![]() Выделим действительную и мнимую части: ![]() ![]() где ![]() ![]() Приравниваем действительную и мнимую части к нулю: ![]() Мы получили систему уравнений, решением которой и будут частота автоколебаний ![]() ![]() Решая второе уравнение системы относительно ![]() ![]() Решая первое уравнение системы относительно переменной ![]() ![]() ![]() ![]() Подставляя численные значения, получаем: ![]() ![]() Построим график функции ![]() с помощью программы MathCad: ![]() ![]() ![]() ![]() В результате получаем, что ![]() Условие устойчивости После того, как были определены амплитуда ![]() ![]() В основе анализа устойчивости периодических процессов в системе управления лежит исследование расположения корней характеристического полинома гармонически линеаризованной системы ![]() Периодический процесс в системе имеет место в том случае, если амплитудно-фазовая характеристика разомкнутой гармонически линеаризованной системы ![]() проходит через точку ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() Итак, требуется, чтобы на данной частоте ![]() ![]() ![]() или ![]() Отсюда следует, что на рис. ![]() положительный отсчет амплитуды ![]() ![]() ![]() Аналитические зависимости для анализа устойчивости периодических процессов в нелинейных системах управления можно получить с помощью критерия Михайлова для гармонически линеаризованной системы. Для этой цели в характеристический полином (1) гармонически линеаризованной замкнутой системы подставим ![]() ![]() В приведённом выше выражении выделим действительную и мнимую части ![]() Так как мы исследуем периодический процесс в системе управления, то годограф Михайлова должен пройти через начало координат плоскости ![]() ![]() ![]() ![]() ![]() Пусть теперь в (2) ![]() системе будут иметь место расходящиеся колебательные процесс- ![]() сы. Таким образом, можно сделать следующий вывод. Для устойчивого периодического процесса в нелинейной системе управления требуется, чтобы при ![]() ![]() Условие устойчивости периодических процессов в нелинейной системе управления формулируется следующим образом. В системе возникают устойчивые, периодические процессы (автоколебания), если: 1. Выполняется условие ![]() 2. Все остальные корни характеристического полинома (1) гармонически линеаризованной системы (кроме использованной нами пары чисто мнимых) имели отрицательные вещественные части, т. е. полином ![]() удовлетворял критерию Михайлова (Гурвица). Т. к. в нашем случае действительная и мнимая части будут равны: ![]() Тогда частные производные будут равны: ![]() ![]() ![]() ![]() То условие устойчивости примет следующий вид: ![]() ![]() То есть: ![]() Вычислим ![]() ![]() ![]() При подстановке чисел получим следующее неравенство: ![]() При подстановке чисел, получаем следующие решения данного неравенства: a > 1 Cледовательно, при значении амплитуды ![]() ^ Коэффициент передачи линейной части в нашем случае будет ![]() В нашем случае характеристический полином линеаризованной системы выглядит следующим образом: ![]() Выделим действительную и мнимую части: ![]() ![]() где ![]() ![]() Исходя из выше изложенного, критический коэффициент передачи будем определять, решая следующую систему: ![]() В результате имеем следующее выражение для ![]() ![]() Из этого выражения определим критический коэффициент передачи линейной части при значении амплитуды ![]() ![]() ^ . На выходе системы получаем: ![]() Теперь покажем, как при изменении коэффициента передачи линейной части относительно его критического значения будет меняться график на выходе. Возьмём коэффициент передачи линейной части больше его критического значения, например ![]() ![]() Т.е. при ![]() ![]() Теперь возьмём коэффициент передачи линейной части меньше его критического значения, например ![]() ![]() Т.е. при ![]() ![]() При к=11: ![]() А схема моей системы в MatLab выглядит следующим образом: ![]() С помощью пекета MatCad составим программу для построения фазового портрета системы: ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() С помощью пекета MatCad составим программу для построения фазового портрета системы при критическом значении коэффициента передачи линейной части: ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() Вывод: В этом курсовом проекте я провела исследование устойчивости состояния равновесия нелинейной следящей системы с одним нелинейным звеном типа «ограничение». Для этого я выполнила линеаризацию нелинейности, используя метод гармонической линеаризации. Затем провела анализ устойчивости системы и определила значение критического коэффициента передачи линейной части. И в заключении, я смоделировала систему в пакете программ MatLab 6.5. Моделирование показало, что полученное мной аналитически значение коэффициента передачи линейной части, верно. Список используемой литературы: 1. Курс лекций по ТАУ Пузанова В.П за 5-6 семестры. 2. Методическое пособие Пузанов В.П, Челышев В.А «Расчет параметров автоколебательных процессов в нелинейных системах автоматического управления и регулирования» 3. Ануфриев И.К «Самоучитель MatLab 6.5» 4. Гультяев А «Визуальное моделирование в среде MatLab»
|