Теория автоматического управления. Линейные системы управления
.pdf7. В результате подстановки (4.22) и (4.29) в векторное уравнение Ла гранжа (4.14) и последующих преобразований математическая модель стенда описывается системой трех обыкновенных дифференциальных уравнений:
2[w2JCJi:+^3(^ + ^cosa)(i:--/sinad)](pH-[J| + т2х2 + w3(x + /cosa)2]cp =
= Щ |
|
- м сХ, |
|
(т2 + |
cos a |
a 2 - т 31sina a ~[{m2 + m3)x + m3l cosa](p2 = |
|
^ F 2 |
- |
Fc2, |
( |
m3l2a - m 3l sina x - m 3/cosa(x- /sina)cp2 =
= M 3 - M cз - m3 g l cos a .
(4.30)
Уравнения в форме (4.30) обычно преобразуют к нормальной форме (4.15) Коши, разрешая их относительно обобщенных ускорений.
Как следует из (4.30), рассматриваемый механизм является нелиней ным объектом третьего порядка с перекрестными связями по обобщенным координатам и скоростям, однако при определенных условиях (раздел 4.5) он может быть аппроксимирован линейной моделью.
В гл. 5 приведены примеры составления динамических моделей элек тродвигателя постоянного тока в форме обыкновенных дифференциальных уравнений, операторных уравнений, структурных схем, сигнальных графов, а в гл. 9 - в векторно-матричной форме.
Для решения уравнений динамики САУ целесообразно применение широко распространенных математических систем {Maple V R3...87, Matlab
5.0...7.0, MathCAD Pro {Premium, 11, 12), Mathematica 2...5 и др.). Прежде всего, пакеты символьной математики этих математических систем позво ляют преобразовать векторно-матричное описание (4.30) к нормальной фор ме Коши, исключив многочисленные рутинные преобразования. Кроме того, эти математические системы позволяют выполнить численное интегрирова ние систем линейных и нелинейных дифференциальных уравнений, а также наглядно представить результаты вычислений в виде таблиц и графиков и решить целый перечень задач математического моделирования динамиче ских систем.
4.5. Линеаризация САУ
При анализе и синтезе САУ применяются математические модели (ММ), представляемые, как правило, обыкновенными дифференциальными уравнениями (ОДУ). Большинство реальных ОУ в широком диапазоне изме нения их переменных являются нелинейными, однако, как показывает прак
тика, в области достаточно малых отклонений координат (переменных) они могут быть аппроксимированы линейными ММ. Свойство линейности ММ ОУ позволяет при исследовании САУ воспользоваться преобразованием Ла пласа к ММ в форме ОДУ и свести интегрирование ОДУ к простым алгеб раическим преобразованиям. Кроме того, линейность преобразований и по лучаемых линейных подпространств координат лежит в основе векторно матричных моделей САУ и их исследования в пространстве состояний, т. е. во временной области. Последнее обстоятельство позволяет применить при синтезе и анализе САУ упоминаемые ранее компьютерные системы Matlab, MathCAD, Maple V, Mathematica и др., базирующиеся на матричных методах исследования линейных систем.
Любая линейная система удовлетворяет свойствам суперпозиции и го могенности. Первое свойство означает, что произвольная сумма аддитивных воздействий X\(t) + х2(1) на входе САУ дает реакцию y\(t) + y2(t) на выходе САУ. Второе свойство предполагает выполнение условия масштабируемо сти, т. е. при изменении входного сигнала х х в к раз выходной сигнал у\ ли нейной САУ изменится соответственно в к раз.
Подавляющее большинство механических и электрических элементов САУ являются линейными в достаточно широком диапазоне изменения их переменных (координат) относительно стационарного режима, чего нельзя сказать о гидравлических, пневматических, термодинамических и иных эле ментах. Вместе с тем, даже такие элементы САУ можно линеаризовать при условии достаточно малых отклонений координат в окрестности точки ста ционарного режима (рабочей точки).
Любую непрерывную функцию у{х) в окрестности рабочей точки х = х0 можно разложить в ряд Тейлора
, ( х ) . у ( х 0) |
Л \ х . х |
1! |
dx2 'х ~х° |
2! |
(4.31) |
|
7 |
^ 0 |
dx'x ~x° |
|
|||
В окрестности рабочей точки |
при малых отклонениях переменной х от |
|||||
х0 выражение (4.31) можно аппроксимировать линейной формой |
|
|||||
Жх) = у(хо) + ^ | *=*„ (х - х0) = у 0 + к(х - |
х0), |
(4.32) |
||||
где к - |
тангенс угла наклона касательной к кривой у(х) в точке х0. |
|
||||
Выражение (4.32) можно преобразовать к виду |
|
|||||
У(х)~Уо = к ( х - х 0) |
|
|
|
(4.33) |
||
или |
|
|
|
|
|
|
Ду = кАх. |
|
|
|
|
(4.34) |
|
Данный метод линеаризации иногда еще называют методом касатель |
||||||
ной линеаризации в рабочей точке х0 или вдоль рабочей траектории |
(*())• |
Рассмотрим простейший пример линеаризации модели идеального ма ятника как статического элемента САУ (см. рис. 4.4).
|
Уравнение движения маятника, отклоненного на угол ф от вертикали, |
|
в осях “угол - вращающий момент” в соответствие с (4.20) имеет вид |
|
|
|
М =-mgl БШф, |
(4.35) |
где |
М - вращающий момент; |
|
|
g - ускорение свободного падения. |
|
|
Линеаризуем (4.35) в окрестности рабочей точки М0(ф0), где ф0 при |
|
мем равным нулю: |
|
|
|
Ш = -mgl соБфо (Ф - ф0) = -mgl($ - Ф0) |
(4.36) |
или |
AM - -w g/Дф. |
(4.37) |
|
Рассмотрим пример линеаризации нелинейного уравнения, описываю |
щего зависимость электромагнитного момента М двигателя постоянного то ка от тока якоря /я и магнитного потока Ф,
М = С мФ /я, |
(4.38) |
где Смконструктивная постоянная двигателя.
Уравнение (4.38) относится к классу нелинейных уравнений, посколь ку содержит произведение координат электродвигателя - магнитного потока и тока якоря. Линеаризуем (4.38) в окрестности рабочей точки М0(Ф0, /яо), соответствующей, например, номинальному режиму работы двигателя, т. е. При Мо —Мн, Фо~‘ ФН1*я0 *ян •
М0 + АМ = см(ф 0 + АФ)(/я0 + ) • |
(4.39) |
Пренебрегая в (4.39) произведением приращений координат, получим |
|
линеаризованное уравнение в приращениях |
|
АМ = См(Ф0Д*;,+ /я0ДФ). |
(4.40) |
В этом уравнении Ф0 и /я0 предполагаются величинами постоянными, а следовательно, уравнение (4.40) относится к классу линейных (линеаризо ванных в рабочей точке) уравнений.
Если управление двигателем осуществляется одновременно по цепям якоря и магнитного потока (цепи возбуждения двигателя), то рабочая точка в процессе управления будет смещаться относительно начального (номиналь ного) режима, образуя семейство рабочих точек или рабочую траекторию. В этом случае при применении уравнения (4.40) говорят о линеаризации ис ходного нелинейного уравнения (4.38) вдоль рабочей траектории
-Мо ~ Сл,Фо i я0*
Помимо касательной линеаризации при исследовании нелинейных САУ в частотной области применяют метод гармонической линеаризации, а при исследовании стохастических САУ - стохастической линеаризации [5].
5. Структурные методы исследования линейных САУ
Исследование линейных динамических систем, математические моде ли которых представлены в форме обыкновенных дифференциальных урав нений, может быть упрощено, если вместо рассмотрения переменных, харак теризующих состояние системы во времени (оригиналов), рассматривать со ответствующие им изображения, получаемые на основе преобразования Ла пласа. Применение такого подхода позволяет заменить операции дифферен цирования и интегрирования более простыми арифметическими операциями умножения и деления изображения на оператор р и, тем самым, значительно облегчить исследование САУ. Преобразования Лапласа лежат в основе по лучения передаточных функций линейных динамических звеньев, объектов и систем управления. Они также являются основой составления и преобразо вания структурных схем - удобной и наглядной формой структурного опи сания и исследования динамических свойств линейных САУ.
5.1. Преобразование Лапласа, передаточные функции и матрицы
Возможность линеаризации технических объектов и систем предостав ляет в распоряжение исследователей аппарат преобразования Лапласа.
Прямое преобразование Лапласа (определение изображения) непре рывной функции времени^) имеет вид [1, 2]
F{p) =*'\f(t)e-i’,dt = L{f{tj\, |
(5.1) |
|
О |
|
|
где р - оператор Лапласа (комплексная переменная). |
|
|
Обратное преобразование Лапласа (определение оригинала) функции |
||
комплексной переменной F(p) имеет вид |
|
|
1 |
а+j(D |
|
Д 0 = г Ц |
1 F{p)ep'dp. |
(5.2) |
271/ |
(J.JO |
|
При решении большинства практических задач применяются таблицы преобразований Лапласа, полученные на основе выражения (5.1). В табл. 5.1 приведены некоторые наиболее важные соотношения непрерывной функции времени t и функции комплексной переменной р.
Как видим, переменную р в преобразовании Лапласа можно рассмат ривать как оператор дифференцирования при описании САУ во временной области, т. е.
cl
(5.3)
Таблица 5.1
Соотношения непрерывной функции времени / и функции комплексной переменной р
Наименование
функции
Импульсная
Импульсная с чистым запаздыванием Ступенчатая
Линейная
Степенная
Асимптотически зату хающая
Асимптотически нарас тающая
Синусоидальная
Косинусоидальная
Синусоидальная асим птотически затухающая
Косинусоидальная асимптотически зату хающая Дифференцирующая
Оригинал функции
л о 5(0
5 ( t - Т )
КО
t
?
е - ° '
1 - е ~ а'
sin(co/)
COS(CD t)
e~at sin(co/)
e~al cos(co/)
'г? /—v |
II |
Изображение функции
Ftp)
1
е-т'
1
P
1
P2 n\
Pn+X
1
p + a
a
P(P + a )
CO
p 2 + 0)2
P
p2 +io2 CO
(p + a ) 2 +ш 2 p + CL
(p + a ) 2 + oo2
p kF ( p ) - p kAm - Pk-2m -
. . . - p Y k- \ o)
Механизм использования преобразования Лапласа для упрощения ре шения уравнений динамики САУ, т. е. определения реакции выходной коор динаты системы на некоторое входное воздействие, предполагает следую щую последовательность действий:
1) составление дифференциальных уравнений системы;
2) прямое преобразование по Лапласу этих дифференциальных уравнений, используя уравнение (5.1) или таблицы преобразований;
3)решение полученных алгебраических (операторных) уравнений от носительно выходной переменной;
4)получение оригинала выходной переменной как функции времени, используя обратное преобразование Лапласа (5.2) или таблицы преобразова ний (отдельно для каждого входного воздействия, если их несколько);
5)суммирование реакций при необходимости оценки влияния на сис тему одновременно нескольких входных воздействий (реализация принципа наложения, или, иначе, принципа суперпозиции линейных систем).
Преобразования Лапласа лежат в основе получения одной из основных характеристик линейных стационарных (с постоянными параметрами) сис
тем - передаточной функции, описывающей динамическую связь перемен ных отдельных звеньев системы и системы в целом в терминах “входвыход”
Передаточная Функция W(p) линейной САУ - отношение преобразова ния Лапласа выходной переменной У(р) к преобразованию Лапласа входной переменной X(р) при условии, что все начальные условия (К(0), Г (0)...) сис
темы при t = 0 равны нулю, т. е. |
|
|
|
П р ) |
= 0 ... |
(5.4) |
|
Щ р) = Х ( РУ |
|||
I |
|
||
Заметим, что передаточная |
функция системы |
несет информацию о |
взаимосвязи только се выходной и входной переменных и не несет никакой
информации о ее внутренних переменных. |
|
|
||||
Пусть дифференциальное уравнение системы имеет вид |
|
|||||
d nY |
d „-\y |
+ + a0Y - b„_ |
d ”- ' x + Ь„ |
d ”~2X + ...b0X , |
(5.5) |
|
dt” |
||||||
dt" - + а,Л - 1 |
dt”'' |
dt”' 2 |
|
|||
где X(t\ Y(t) - |
соответственно входное воздействие и выходная реакция |
|||||
системы, |
|
|
|
|
|
а, , bj - постоянные коэффициенты, i =1,..., п - \ .
Тогда с учетом (5.3), (5.4) получим передаточную функцию системы
|
|
л-1 |
+ bcs |
|
^ ) = 2М |
= . |
Ъ„~\Р |
(5-6) |
|
Х(р) |
р |
+ап_хр п |
+ +а0 |
|
Если полином, стоящий в знаменателе, приравнять нулю, то мы полу чим характеристическое уравнение, названное так потому, что его корни (полюса) определяют характер движения системы. Корни полинома, стояще го в числителе, называют нулями системы. В полюсах передаточная функция (5.6) обращается в бесконечность, в нулях - в нуль. Расположение полюсов и нулей на комплексной /7-плоскости определяет характер собственного (сво бодного) движения системы. Выбор рационального расположения полюсов и нулей замкнутой системы за счет применения корректирующих звеньев - один из широко применяемых подходов к синтезу САУ (см. гл. 10).
причем входные воздействия могут оказывать влияние на несколько выход ных (управляемых) переменных. Такие системы называют многосвязными многомерными. Для их описания и исследования применяют не передаточ ные функции, а передаточные матрицы.
Пусть число управляемых величин равно и, число задающих воздейст вий т, а число возмущающих воздействий к. Совокупности выходных, за дающих и возмущающих переменных можно записать в виде соответствую щих векторов-столбцов X=coI(x,... х„), Z=col(z,... zm), F=coI(/i.../t). Уравне ние такой н-связной системы можно записать в векторно-матричной опера
торной форме |
|
|
|
А(р)Х(р) = В(рУДр) +C(p)F(p), |
|
(5.7) |
|
где А(р), В(р), С(р) - операторные передаточные матрицы, |
|||
Яц(Р) |
а \п(р) |
i'll(Р) |
К Л р Т |
А(р) = |
; в(/>)= |
|
3 |
_а„| (Р) |
ат,(.Р). |
Al ( р ) |
Кт(Р)_ |
си (р) |
С1к(Р) |
|
|
С(р) = |
|
|
|
Snl(p) |
С„к(Р)_ |
|
|
Разрешая уравнение (5.7) относительно вектора выходных перемен |
|||
ных, получим |
|
|
|
Х(р) = A' 1(p)B(p)Z(p) + A -‘(p)C(p)F(p). |
(5.8) |
Первое слагаемое в (5.8) представляет собой реакцию системы на за дающие воздействия, второе - на возмущающие воздействия. Для получения оригинала выходной переменной как функции времени необходимо исполь зовать обратное преобразование Лапласа (5.2) или таблицы преобразований (отдельно для каждого входного воздействия) и просуммировать реакции при необходимости оценки влияния на систему одновременно нескольких входных воздействий.
5.2. Типовые динамические звенья и структурные схемы САУ
Для наглядного представления структуры сложной динамической сис темы управления как совокупности элементов и связей между ними приме няют структурные схемы и графы. Структурные схемы, как и сигнальные графы, представляют собой графическое изображение структуры САУ. Эле менты структурной схемы САУ представляют в виде типовых динамических звеньев, характеризующихся однонаправленностью, одним входом и одним выходом, что позволяет применить для их описания аппарат передаточных функций. Если реальное динамическое звено не обладает однонаправленно стью, т. е. выход оказывает влияние на вход, то такой элемент представляют
в виде направленного звена с обратной связью. Если у элемента несколько входов, в его структуру включают суммирующие звенья (сумматоры) - спе цифические многовходовые безынерционные звенья с единичными коэффи циентами передачи по каждому входу, причем каждый вход сумматора обо значается знаковой функцией (+ или -). Если у элемента несколько выходов, это означает, что его нельзя рассматривать как элементарное звено и к нему необходимо применить декомпозицию, выделив соответствующие числу вы ходов звенья.
На структурных схемах динамические звенья изображают прямоуголь никами, входные и выходные воздействия - подходящими и отходящими от прямоугольников стрелками и текстовыми надписями, обозначающими их формальный вид. Внутри прямоугольников вводят обозначения передаточ ных функций звеньев. Сумматоры сигналов (переменных САУ), как правило, обозначают окружностями, сами сигналы - подходящими к окружностям стрелками с указанием имен переменных и знаков алгебраического сумми рования (+ или -). Следует отметить, что в технической литературе, компью терных системах автоматизированного проектирования и управления, систе мах сопровождения жизненного цикла САУ (САПР, АСУ ТП, CAD/CAM, SCADA, CALS) встречается множество графических обозначений сумматоров сигналов, однако в силу своей простоты все они интуитивно понятны.
В обобщенной форме структурная схема динамического звена приве дена на рис. 5.3.
Х(Р) |
г«*(р) |
W{p) |
----- ► |
Рис. 5.3. Обобщенная структурная схема
динамического звена
Типовые элементарные динамические звенья, их реакции на единич ные ступенчатое и импульсное воздействия приведены в табл. 5.2.
В качестве примера составления структурных схем динамических объ ектов управления рассмотрим электродвигатель постоянного тока, регули руемый по цепям якоря и возбуждения [4]. Функциональная схема и схемы замещения электродвигателя приведены на рис. 5.4.
В структуре электродвигателя можно выделить три основных цепи (см рис. 5.4,6, 5.4,в, 5.4,г):
- цепь якоря, питаемая регулируемым напряжением С/я; R3, L3 - COOT ветственно эквивалентное активное сопротивление и эквивалентная индук тивность якорной обмотки; £д-э.д.с. электродвигателя; /„ —ток якоря;
- цепь возбуждения, питаемая регулируемым напряжением U0, RB,-Lа' соответственно эквивалентное активное сопротивление и эквивалентная ин дуктивность обмотки возбуждения; /в - ток возбуждения;