Построение годографов Михайлова при помощи пакета MATHCAD (стр. 1 из 2)
построение годографов михайлова при помощи пакета «mathcad»
Цель работы заключается в необходимости получения простого и наглядного инструмента для решения задач расчёта устойчивости систем автоматического управления, что является обязательным условием работоспособности любого промышленного робота и манипулятора.
1 Понятие об устойчивости системы
Как видно из цели исследования, необходимым условием работоспособности системы автоматического управления (САУ), является её устойчивость. Под устойчивостью принято понимать свойство системы восстанавливать состояние равновесия, из которого она была выведена под влиянием возмущающих факторов после прекращения их воздействия [1]. Если система не способна возвращаться в состояние равновесия, которое было нарушено в процессе работы, то для практического использования она непригодна.
На практике для определения устойчивости САУ используют критерии устойчивости, то есть правила, с помощью которых можно определить устойчива ли система, не прибегая к решению дифференциальных уравнений. Одним из таких критериев, есть критерий устойчивости Михайлова.
2 Критерий устойчивости Михайлова
Данный критерий основан на связи характера переходного процесса системы с амплитудой и фазой вынужденных колебаний, устанавливающихся в системе при синусоидальном воздействии. Анализ устойчивости системы этим методом сводится к построению по характеристическому многочлену замкнутой системы (знаменатель передаточной функции), комплексной частотной функции (характеристического вектора):
(1)где
и 
– соответственно вещественная и мнимая части знаменателя передаточной функции, по виду которой можно судить об устойчивости системы.Если задаваться различными значениями частоты
и откладывать по горизонтальной, а по вертикальной осям декартовой системы координат, то будет получена кривая, называемая годографом характеристического вектора или годографом Михайлова.В таком случае, критерий устойчивости Михайлова может быть сформулирован следующим образом: замкнутая САУ устойчива, если комплексная частотная функция
, начинаясь на действительной положительной оси, при изменении частоты от 0 до ∞ огибает против часовой стрелки начало координат, проходя последовательно n квадрантов, где n – порядок характеристического уравнения системы, т. е. 
(2) 
Рисунок 1 – Амплитудно-фазовые характеристики (годографы) критерия Михайлова: а) – устойчивой системы; б) – неустойчивой системы (1, 2) и системы на границе устойчивости (3)
На рис. 1 показаны примеры перемещения годографов Михайлова для различных систем с изменяющимся порядком n характеристического уравнения.
3 Алгоритм построения годографа Михайлова
Рассмотрим последовательность расчёта критерия устойчивости Михайлова и сформируем алгоритм построения годографа, используя математический пакет «MathCad», на приведенных ниже примерах.
Пример 1. Используя критерий Михайлова, определим устойчивость системы автоматического управления электроприводом манипулятора промышленного робота (МПР). Структурная схема САУ электроприводом МПР изображена на рис. 2.
Рисунок 2 – Структурная схема САУ электроприводом МПР
Передаточная функция данной САУ имеет следующее выражение [2]:
(3)или
(4)где kу – коэффициент усиления усилителя, kм – коэффициент пропорциональности частоты вращения двигателя величине напряжения на якоре, Tу – электромагнитная постоянная времени усилителя, Tм – электромеханическая постоянная времени двигателя с учётом инерции нагрузки (по своим динамическим характеристикам двигатель представляет собой передаточную функцию последовательно соединённых инерционного и интегрирующего звеньев), kдс – коэффициент пропорциональности между входной и выходной величинами датчика скорости, K – коэффициент усиления главной цепи:
.Подставим численные значения в выражение передаточной функции:
K = 100 град / (В∙с); kдс = 0,01 В / (град∙с); Tу = 0,01 с; Tм = 0,1с.
Получим:
(5)Далее запишем характеристический многочлен замкнутой системы
заменив s на 
: 
(6)С помощью (1) выделим вещественную и мнимую части и подставим численные значения в полученную комплексную частотную функцию:
(7)Имея данные в виде (7), перейдём непосредственно к использованию математического пакета «MathCad».
Для этого в верхнем меню выберем «Новый…» – «Пустой документ», в котором будем формировать программу построения годографа Михайлова, используя нижеприведенный алгоритм.
Шаг 1. Задать разрешение годографа диапазоном значений индекса i. Например:
(8)Шаг 2. Определить исследуемый диапазон и шаг частоты
, используя значения индекса i (обычно, для практических расчётов, максимальная величина частоты 
не превышает значения 1000, в нашем же примере – достаточно принять 
с частотным шагом 0,1):
(9)Шаг 3. Полученные вещественную
и мнимую 
части характеристического уравнения, зададим численными значениями (в данном случае используя (7)) в виде: 
(10) 
иРисунок 3 – Массивы значений , и , рассчитанные в «MathCad»
(11)Шаг 4. В результате вычислений (9), (10) и (11), получаются массивы значений частоты
, а также вещественной 
и мнимой 
частей (рис. 3).Шаг 5. Далее имея рассчитанные массивы значений
и 
, переходим к построению годографа Михайлова, используя встроенную функцию «MathCad» – «Инструменты графиков», выбрать «Декартов график». Здесь необходимо определить идентификаторы осей (в данном случае ось абсцисс соответствует вещественной части , а ординат – мнимой части 
) и параметры графика в подменю «Формат…». В результате получим график комплексной частотной функции, приведенный на рис. 4. 
Рисунок 4 – Годограф Михайлова для САУ электроприводом МПР
Используя функцию «Трассировка…» (пунктирные линии на рис. 4), можно определить, в соответствующем трассировке окне, точные значения годографа в любой точке рассчитанных массивов.
Таким образом, по рассчитанным данным, построенный годограф Михайлова, начинаясь на действительной положительной полуоси, огибает против часовой стрелки начало координат, проходя последовательно три квадранта, что соответствует порядку характеристического уравнения. Следовательно, данная САУ электроприводом МПР – устойчива.
В соответствии с изложенным алгоритмом, рассмотрим ещё один пример расчёта критерия устойчивости Михайлова и построения комплексной частотной функции.
Пример 2. На современных автомобильных заводах широко применяются большие сварочные роботы (рис. 5). Наконечник сварочного узла (НСУ) подводится к различным местам кузова автомобиля, быстро и точно совершает необходимые действия. Требуется определить устойчивость по критерию Михайлова САУ позиционированием НСУ, структурная схема которой изображена на рис. 6.