VS
Программы
Пособия
Сейсмонитор
en
 

2. Автоматизированный синтез систем

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

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

Рисунок 1. Структурная схема объекта регулирования

При заданных структурах корректирующих связей и неизвестных в них параметрах осуществляется описание всей системы в виде дифференциальных уравнений, характеризующих свободное поведение системы (1):
 (1)
Тогда характеристическое уравнение для (1) примет вид (2):
 (2)
В выражении (2) коэффициенты a являются функциями синтезируемых параметров. Желаемое характеристическое уравнение для рассматриваемой системы имеет вид (3):
 (3)
Все корни λi характеристического уравнения (3) считаются известными и задаются пользователем или автоматически рассчитываются программой. Следовательно, коэффициенты bi уравнения также известны, так как рассчитываются по известным корням уравнения. Тогда уравнения (2) и (3) представляются в виде системы уравнений (4):
 (4)
Приравнивая полученные уравнения, будут определены значения неизвестных параметров искомой системы (5):
 (5)
Таким образом, все неизвестные коэффициенты корректирующих звеньев могут быть получены путём решения системы уравнений (5). После этого полученные значения подставляются в систему, и выполняется её моделирование. Если результаты моделирования удовлетворяют требованиям, то синтез выполнен успешно. В ином случае возможно внесение изменений в требуемые характеристики системы, состоящие, например, в увеличении времени переходного процесса.

Выбор регулятора.
Основным правилом при выборе порядка передаточной функции регулятора является равенство порядка всей системы и количества неизвестных параметров вводимых регуляторов. Кроме того, вводимый регулятор должен иметь как можно меньший порядок. Передаточные функции регуляторов для 0, 1 и 2 порядков приведены ниже (6):
 (6)
 
Так, если неизменяемая часть системы имеет 3 порядок и необходимо ввести один регулятор, то его передаточная функция может быть только W22(p), т.к. тогда порядок системы будет равен 3+2 и число неизвестных параметров (B0, B1, B2, A0, A1) также равно 5.

При выборе корней желаемого характеристического уравнения действует следующее правило. Для того, чтобы система стала устойчивой, время переходного процесса которой будет приблизительно равно tпп при колебательности μ, необходимо выбрать доминантный корень в левой полуплоскости (7):
 (7)
Значение μ выбирается исходя из технических требований. При отсутствии дополнительной информации значение показателя принимается равным 1,6. Начальное значение действительной части α0 рассчитывается следующим образом (8):
 (8)
Все последующие значения действительной части увеличиваются в 1,5–3,0 раза относительно предыдущего значения, а значение мнимой части пересчитывается согласно значению действительной части.

Пример синтеза простого объекта.
На рисунке 2 представлена структурная схема объекта с последовательным корректирующим устройством, требуется обеспечить установление всех процессов за 1 секунду. Порядок системы без введения корректора равен 1,. Так как выбирается простейший корректор, порядок которого не увеличит порядок системы, а число неизвестных параметров будет равно общему порядку системы с корректором, следовательно, корректирующее устройство для рассматриваемой системы будет равно W0(p)=K.
Рисунок 2. Структурная схема исследуемого объекта

По структурной схеме составляется матричная математическая модель системы, которая имеет вид (9):
 (9)
По матричной математической модели записываются системы дифференциальных и алгебраических уравнений (10):
 (10)
Выражая выходы звеньев и подставляя их в дифференциальное уравнение, будет получено дифференциальное уравнение и матрица A (11):
 (11)
После этого в программе синтеза SintACS вводится матрица A (см. рисунок 3), задаётся время переходного процесса и запускаются алгоритмы расчёта.

Рисунок 3. Диалоговое окно SintACS
 
В результате выполнения синтеза будет получено значение параметра K, которое равно 0,3.

Пример синтеза сложного объекта.
Система с тремя неизвестными регуляторами представлена на рисунке 4. Требуется определить структуру регуляторов и найти такие их параметры, чтобы переходные процессы завершились за 0,9 секунд
.
Рисунок 4. Структурная схема системы в общем виде
 
Тогда выбранные согласно правилам регуляторы представлены в структурной схеме на рисунке 5.
Рисунок 5. Структурная схема системы с выбранными звеньями
 
Путем исключения из последующего рассмотрения внешнее воздействие и сумматор, на которое оно действует, будет получена обобщенная структура системы. Полученная система имеет пятый порядок, а число синтезируемых параметров равно 5. Так как порядок системы и число неизвестных параметров равны, то возможно применение синтеза корректирующих устройств. 

Матричная математическая модель исследуемой системы имеет вид (12):
 (12)
По матричной математической модели составляются дифференциальные уравнения (13):
 (13)
По матричной математической модели составляются алгебраические уравнения (14):
 (14)
После этого в (14) выражаются все сигналы (15):
 (15)
Тогда система дифференциальных уравнений (13) примет вид (16):
 (16)
Полученная система дифференциальных уравнений (16) представляется в матричном виде (17):
 (17)
На основе матрицы (17) производится синтез в программе SintACS, диалоговое окно которой представлено на рисунке 6.

Рисунок 6. Диалоговое окно SintACS для сложной системы

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

Рисунок 7. Модель исследуемой системы в SimACS

Выполняя моделирование системы, будет получен график выходного сигнала системы (см. рисунок 8).

Рисунок 8. График переходного процесса системы

Как следует из приведённого графика, время переходного процесса составляет 0,7 секунды и не превышает заданного в техническом задании времени 0,9 секунды. Следовательно, синтез системы проведён успешно.

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

Тогда дифференциальные уравнения в матричном виде для исследуемой системы примут вид (18):
 (18)
В программе SintACS записывается матрица A, выполняются вычисления, и в строках будут получены искомые коэффициенты характеристического полинома, которые представлены на рисунке 10.
Рисунок 10. Диалоговое окно SintACS при нахождении коэффициентов характеристического полинома

Таким образом, в программе SintACS получено 2 коэффициента характеристического полинома: zT2+zT1 и zT1*zT2+zT1*K1*zT2*K2. Тогда, выполняя обратную подстановку для zT1 и zT2, будет получен характеристический полином в символьном виде (19):
 (19)
Кроме того, характеристический полином может быть получен классическим способом, приравнивая знаменатель передаточной функции замкнутой системы 1+K1K2/(1+pT1)(1+pT2) к нулю. После преобразований рассчитанное решение будет равно решению (19), полученному в программе SintACS.

Получение передаточных функций по структурной схеме.
Для исследуемой системы, представленной на рисунке 11, получить передаточную функцию замкнутой системы по ошибке в символьном виде при помощи программы SintACS можно следующим образом.
Рисунок 11. Структурная схема исследуемого объекта
 
Передаточная функция по ошибке (Y3/U) в программе SintACS может быть вычислена путём расчёта двух матриц. Для их получения следует записать математическую модель системы без переменных состояния в подматрице алгебраических уравнений. Для динамических звеньев по главной диагонали для выходов записываются знаменатели с обратным знаком, а при входах – числители со своим знаком. Для суммирующих и усилительных звеньев по главной диагонали заносятся –1 (20):
 (20)
Тогда по условию задачи необходимо определить передаточную функцию следующего вида (21):
 (21)
Верхний определитель получен путём замены третьего столбца четвёртым, но с обратным знаком. Он отражает передаточную функцию на сумматоре Y3 при подаче воздействия U.  Нижний определитель соответствует замкнутой системе. Последовательно решая определители в SintACS, будет получена передаточная функция замкнутой системы по ошибке в символьном виде (22):
 (22)
Кроме того, передаточная функция замкнутой системы может быть получена классическим способом и равна 1/(1+K1K2/(1+pT1)(1+pT2)). После преобразований рассчитанное решение будет равно решению (22), полученному в программе SintACS.

Рейтинг@Mail.ru
Все материалы сайта принадлежат лично Василию Щербакову,
а также соответствующим авторам при указании ссылки. Администрация сайта не несет никакой ответственности за стороннюю информацию. Вы можете пожаловаться на контент, если он нарушил Ваши права. Для этого свяжитесь с администрацией.
Допускается использование материалов в некоммерческих или учебных целях с указанием ссылки на этот сайт.