Операционная форма записи дифференциальных уравнений

Операционная форма записи дифференцированных уравнений, примеры.

Операционное управление – есть следствие прямого преобразования Лопласа диф Ур-й САР в основе этого лежит интеграл Лопласа.

Операционная форма записи

Где все а и b – постоянные коэф для САР (n≥m)

Y,x – функции соответствующие входной и выходной величине.

Передаточная ф-ия последовательно и параллельно соединенных звеньев. Примеры

Передаточная ф-ия представляет собой отклонение изоброжения выходной и входной величины.

1) последовательное соединение передаточных чисел.

х1 – промежуточное значение

Y=W2X1 Y=W1W2X=WэквX => Wэкв = W1W2

X1=W1X Wэкв= интеграл от i=1 до 1 по Wi

2) параллельное соединение

Эквивалентная передаточная ф-ия звена охваченного обратной связью пример.

Передаточная ф-я – отношение изображения входной и выходной величины.

Дата добавления: 2015-07-26 ; просмотров: 6 ; Нарушение авторских прав

Формы записи дифференциальных уравнений

Стационарные линейные непрерывные САУ наиболее часто описываются линейными дифференциальными уравнениями с постоянными коэффициентами:

. (2.9)

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

Применяется также операторная форма записи уравнения (2.9):

. (2.10)

В этом уравнении через « » обозначен оператор дифференцирования .

Заметим, что по сложившейся традиции символ « » используется также в преобразованиях Лапласа и Карсона-Хевисайда, но является комплексным числом .

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

(2.11)

или в операторной форме

. (2.12)

Первая стандартная символическая форма записи уравнения (2.11) имеет следующий вид:

, (2.13)

где ; ; ; .

Форма (2.13) представляет собой операторно-структурное описание системы, т.е. в виде операторов звеньев, составляющих структурную схему системы (далее эти понятия разъясняются), и связей между ними. В этой форме — постоянные времени звена, измеряемые в секундах; — передаточный коэффициент звена.

Из изложенного выше следует, что уравнение (2.9) в этой форме перепишется в следующем виде:

, (2.14)

где ; ; .

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

.

Передаточная функция САУ, поведение которой во времени описывается уравнением (2.9), имеет следующий вид :

.

В формуле (2.15) через и обозначены изображения (по Лапласу) выходной и входной переменных САУ при нулевых начальных условиях и равенстве нулю внешних возмущений, а через и — полиномы относительно комплексной переменной .

Вторая стандартная форма записи дифференциального уравнения имеет следующий вид:

или . (2.16)

В (2.16) и являются полиномами (символическими) относительно оператора .

Из сравнения первой и второй стандартных форм записи дифференциальных уравнений следует, что с математической точки зрения различие между этими формами весьма несущественно и состоит лишь в различном представлении коэффициентов уравнений. В ТАУ принято называть уравнения вида (2.9) — (2.14), (2.16) уравнениями типа «вход-выход».

Третья стандартная форма записи дифференциального уравнения принципиально отличается от форм записи, описанных выше. В этой форме записи используются переменные состояния. Отметим, что понятие «состояние» является базовым в современной ТАУ (СТАУ). Переменные состояния — это промежуточные переменные системы (рис.2.2), число которых равно ее порядку . В общем случае входные и выходные переменные могут быть векторными величинами размерности и соответственно.

Координаты состояния х1, х2 , . ,хn u y Рис.2.2 — Состояние системы

Переменные состояния называют также координатами состояния, так как их совокупность задает вектор состояния .

Множество возможных положений этого вектора образует векторное пространство , называемое пространством состояний системы. В переменных состояния САУ описывается векторно-матричным уравнением

, (2.17)

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

— вектор переменных состояния — внутренних координат системы;

— вектор входных переменных (управляющих и возмущающих);

— вектор наблюдаемых или выходных переменных; размерности матриц , , , соответственно, ( ), ( ), ( ).

Процессы в САУ в свободном движении (без внешних воздействий) согласно уравнению (2.17) описываются векторно-матричным уравнением с характеристическим уравнением , где — единичная матрица, или в развернутом виде системой дифференциальных уравнений

с характеристическим уравнением

. (2.18)

Эти уравнения при определенных начальных условиях дают возможность изучить процессы в системе путем их решения численными методами с использованием ЭВМ.

Разработаны различные способы перехода от уравнений типа «вход-выход» к уравнениям состояния вида (2.17) и наоборот. Один из наиболее распространенных способов состоит в следующем. Пусть САУ описывается уравнением (2.9). Введем обозначения

, , . , ,

.

С помощью этих обозначений преобразуем уравнение (3.1) к следующему виду:

, (2.19)

где ; ;

; .

В нашем примере и являются скалярными величинами. В общем случае (2.17) — это, соответственно, вектор наблюдаемых или выходных переменных и вектор входных переменных (управляющих и возмущающих), поэтому в (2.19) матрицы и выродились в вектор-столбец и вектор-строку соответственно.

Система уравнений (2.19) представляет собой описание линейной непрерывной системы в пространстве состояний . Уравнения (2.19) с матрицей называют уравнениями в форме Фробениуса.

Если , то

; .

Форма уравнений (2.19) с подобными матрицами и называется в ТАУ канонической формой фазовой переменной.

Задание 1

1.1. По дифференциальному уравнению системы:

Для каждого типового звена 1 – 12 (таблицы 2.1) в соответствии с его параметрами вывести дифференциальное уравнение, операторное уравнение, и выражение передаточной функции.

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

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

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

Третья стандартная форма записи дифференциального уравнения — переменные состояния.

Таблица 2.1 – Исходные коэффициенты

Задание 2

2.1Для каждого звена (таблицы 2.2) по его передаточной функции записать дифференциальное уравнение.

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

Вар	Передаточная функция	Значения параметров передаточной функции
	а0=1; а1=5; а2 =1,2; а3 =0,9; а4=0,5;в0=1;в1=3; в2=0,8;в3=0,3
	а0=1; а1=5; а2 =1,2; а3 =0,9; в0=1;в1=3; в2=0,8
	а0=1; а1=5; а2 =1,2; в0=1; в1=3;
	а0=1; а1=5; а2 =1,2; а3 =0,9; а4=0,5; в0=10
	а0=1; а1=5; а2 =1,2; а3 =0,9; в0=10
	а0=1; а1=5; а2 =1,2; а3=0,9;а4=0,5;в0=1;в1=3; в2=0,8;в3=0,3
	Т0=2; Т1=4; Т2=1,1;Т3=0,9
	Т0=2; Т1=4; Т2=1,1;Т3=1,1;Т4=,9
	К= 10;Т1=4; Т2=1,1;Т3=0,9
	К= 10; Т2=1,1;Т3=0,9 Т4=0,9
	Т0=0,7; Т1=3;Т2=1,2;Т3=0,9;Т4=0,8;Т5=0,5
	К=10 Т0=0,7; Т1=3;Т2=1,2;Т3=0,9;Т4=0,8;Т5=0,5;

Задание №3

3.1 Для заданной схемы необходимо составить операторное уравнение для каждого элемента схемы САУ.

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

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

3.4 Сформировать математическое описание систем автоматического регулирования в виде третьей стандартной формы записи дифференциального уравнения — В переменных состояния САУ описываемых векторно-матричным уравнением.

Схема, показанная на рисунке 2.2, представляет собой САР температуры в помещении. Объектом регулирования (ОР) в дан­ной системе является помещение, для которого регулируемая ве­личина — температура внутри помещения Ө, регулирующее (уп­равляющее) воздействие — температура воздуха Ө_К, поступающего из калорифера, возмущающее воздействие — изменения внешних факторов f(в общем случае изменение температуры атмосферного воздуха, его влажности, скорости ветра). При исследовании сис­темы в качестве основного возмущения следует рассматривать из­менение температуры окружающего воздуха.

Воспринимающим органом — ВО (датчиком, чувствительным элементом) в данной САР является терморезистор R_Д, включен­ный в мостовую схему, обеспечивающую с помощью резистора R_Озадание необходимого значения температуры в помещении и выполняющую также функции сравнивающего органа — СО (эле­мента сравнения). Усиление сигнала разбалансаΔU(сигнала рас­согласования) измерительной мостовой схемы обеспечивается посредством усилителя. Усиленный сигнал Uобеспечивает вра­щение двухфазного исполнительного двигателя, который изменя­ет перемещение клапана (заслонки) на трубопроводе подачи парав калорифер, чем достигается изменение температуры воздуха на входе калорифера — регулирующего воздействия на объектерегулирования.

1 — помещение; 2 — теплообменник (калорифер), 3 — измерительная мостовая схема; 4 — двухфазный ис­полнительный двигатель, 5 — дифференциальный магнитный усилитель; 6 — клапан (заслонка)

Рис. 2.2. Схема САР температуры

Динамические свойства объекта регулирования и элементов системы описываются следующими уравнениями:

где То, Т2, Т3, Т4 — постоянные времени, с; Ө — значение температуры воздуха в помещении, °С, Ө к — значение температуры воздуха на выходе калорифера, °С; к, к₁, к₂, к₃, к₄— коэффициенты передачи; f— возмущающее воздействие на объекте регулирования; Uд —падение напряжения на термодатчике, В; ΔU— напряжение на выходе мостовой схемы (сигнал рассогласования), В; μ. — линейное перемещение клапана, см; U₀ — задающий сигнал, В.

Значения параметров элементов САР по вариантам даны в таб­лице 2.3.

Заданное значение температуры в помещении Ө = 20 °С.

Значения параметров элементов САР

Примечание. Для всех вариантов постоянные времени Т3 = 20 с, Т4=0,5 с.

Схема САР, приведенная на рисунке 2.3, обеспечивает стаби­лизацию угловой скорости электродвигателя постоянного тока который совместно с рабочим механизмом является объектом ре­гулирования. Регулируемая величина объекта — угловая скорость двигателя ω, регулирующее воздействие — напряжение Uг,пода­ваемое от генератора на якорь двигателя. Возмущающее воздейст­вие на объекте регулирования — момент сопротивления М_с, соз­даваемый рабочим механизмом. Угловая скорость двигателя ωконтролируется тахогенератором, сигнал которого Uтг, пропор­циональный скорости, сравнивается с задающим сигналом U₃. Сигнал рассогласования ΔU = U₃— U_Tг усиливается магнитным усилителем и воздействует на обмотку возбуждения генератора, выполняющего функции исполнительного органа (элемента).

Динамические свойства объекта регулирования и элементов САР описываются следующими уравнениями:

гдеТ_д, Ту, T_v — постоянные времени, с; К_д, К_м, К_тг, К_у, К_г — коэффициенты передачи соответствующих элементов систем

1 — задающий потенциометр; 2 — магнитный усилитель; 3 — генератор; 4 — двигатель; 5 — тахогенератор; 6 — рабочий механизм

Рис. 2.3. Схема САР угловой скорости электродвигателя

Значения параметров элементов САР

Значения параметров объекта регулирования и элементов сис­темы для различных вариантов указаны в таблице 2.4. Заданное значение угловой скорости ω = 40 рад/с.

На рисунке 2.4 изображена схема САР давления Р в ресивере (воз­духосборнике) 1, который является в данной системе объектом регу­лирования. Давление в ресивере регулируется посредством изменения количества воздуха Q, зависящего от положения заслонки 2, т.е. от ее линейного перемещения Х₃, которое можно рассматривать как регу­лирующее воздействие на входе объекта регулирования. Внешним возмущением, вызывающим отклонение регулируемой величины — давления Р, является изменение расхода сжатого воздуха Q_c.

Рис 2.4 Схема САР давления Р в ресивере

Давление в данной системе контролируется с помощью сильфонного датчика 3, выход­ная величина которого — пере­мещение Х_с сильфона 5 одно­значно зависит от разности сил ΔF= F₀— F_p, где F_p— сила, соз­даваемая давлением Р, F₀— си­ла натяжения пружины 6, кото­рое можно изменять винтом 7.

Перемещение сильфона Х_сс помощью потенциометрического преобразователя 4 преобразуется в электрический сигнал — напряжение U, которое усиливается электронным усилителем 8. Выходной сигнал усилителя U_yуправляет электромагнитным при­водом 9, связанным с заслонкой 2,

В данной САР сильфонный датчик выполняет функции вос­принимающего, задающего и сравнивающего органов. Как вос­принимающий орган он контролирует давление Р, преобразуя его в силу Fp. Задание требуемого давления в ресивере обеспечивается посредством силы F0. Как сравнивающий орган сильфон обеспе­чивает сравнение величин F0 и Fp, в результате чего, как отмеча­лось ранее, получается ΔF= F0 — Fp — сигнал рассогласования.

Динамические свойства объекта регулирования и элементов САР описываются следующей системой уравнений:

Физическая сущность переменных, входящих в уравнения, от­ражена выше в описании схемы САР. Параметры T0, T1, T2, T3 и К0, Кq, Кв, Кc, Кп, Ку, К₃ — соответственно постоянные времени и ко­эффициенты передачи. Их размерности и значения по вариантам даны в таблице 2.5. Требуемое значение давления Р = 500 кПа.

Значения параметров элементов САР

На электрических станциях при производстве электроэнергии предъявляют определенные требования к стабильности частоты f генерируемой ЭДС. Частота f однозначно определяется угловой скоростью ω рабочего колеса гидротурбины. В связи с этим гид­ротурбины на электростанциях оснащают САР угловой скорости. На рисунке 2.5 показана схема одного из вариантов такой САР.

В данной системе объектом регулирования является гидротур­бина 1, регулируемой величиной — угловая скорость ω .Она при постоянном расходе воды изменяется в зависимости от нагрузки на валу турбины, т. е. от мощности Р, которая потребляется от ге­нератора 2 (с увеличением мощности угловая скорость снижается, с уменьшением — возрастает). Таким образом, мощность Р явля­ется внешним возмущающим воздействием на объекте регулиро­вания. Для регулирования угловой скорости предусмотрена за­слонка 3, с помощью которой изменяется расход воды через тур­бину. Он однозначно зависит от вертикального перемещения X заслонки. Следовательно, перемещение заслонки X можно рас­сматривать как регулирующее воздействие объекта регулирова­ния. Угловая скорость ω контролируется посредством тахогенератора 4, ЭДС Е которого сравнивается с задающим напряжением U0. Сигнал рассогласования Δ U через усилитель 5 управляет по­средством электродвигателя 6 и редуктора 7 заслонкой 3.

Рис. 2.5 Схема САР угловой скорости рабочего колеса гидротурбины

Динамические свойства элементов САР описываются следую­щей системой уравнений:

Дифференциальные уравнения САУ. Форма вход-выход, операторная форма вход-выход и форма Коши описания САУ с сосредоточенными параметрами

Вернёмся к уравнению прямолинейного движения точечного объекта с переменной массой (4).

,

Где — сила, действующая на объект (вход);

— координата положения объекта – (выход);

— масса объекта (параметр).

Рассматриваемое уравнение называется уравнением в форме вход-выход.

В ТАУ есть неписанные установившиеся традиционные правила, касающиеся символьного описания объекта безотносительно к его природе.

Коэффициенты уравнения принято обозначать греческими буквами

; .

Сила, если она выступает в роли внешнего воздействия, то её обозначают F, а если в роли управляющего, то U.

Таким образом, уравнение (4) примет вид.

. (5)

Любое дифференциальное уравнение n-го порядка можно представить виде системы n уравнений первого порядка или в форме Коши (иначе говорят уравнениями в пространстве состояний). Продемонстрируем это на примере (5). Введём обозначения , , где х_i=(i=1,2) – компоненты вектора состояний. С учётом этих обозначений уравнение (5) можно переписать в виде систем уравнений в форме Коши.

.

Достаточно часто встречаются звенья, имеющие нелинейную зависимость между входной и выходной координатами. Если для малых отклонений от установившегося режима нелинейность несущественна, то в этом случае до составления исходных дифференциальных уравнений САУ выполняют процедуру линеаризации.
Линеаризацией называется замена реальных нелинейных уравнений статических характеристик элементов близкими к ним линейными уравнениями. Линеаризация возможна, если нелинейная характеристика непрерывна и имеет непрерывные частные производные. На рис.2.1. приведена геометрическая интерпретация линеаризации по методу малых отклонений.

Рис.2.1. Геометрическая интерпретация линеаризации

Разложив функцию y=f(x) в ряд Тейлора, получим

где y₀— значение выхода, соответствующее входу x₀; d k y/dx k — значения производных, взятых в точке А(x₀;y₀). Тогда для малых отклонений x:

или

где при x=x₀.

Если выходная величина элемента зависит от нескольких входных воздействий, то при линеаризации по методу малых приращений следует определять частные производные по всем воздействиям, а приращение выхода является суммой частных приращений, т.е.

где x₁, x₂, …, x_n — приращения входных воздействий; — частные производные.

4. Обобщение уравнений динамики САУ

К линеаризованной САУ можно применить принцип суперпозиции: реакция системы на несколько одновременно действующих входных воздействий равна сумме реакций на каждое воздействие в отдельности. Это позволяет звено с двумя входами u и f разложить на два звена, каждое из которых имеет один вход и один выход (рис.27). Поэтому в дальнейшем мы ограничимся изучением поведения систем и звеньев с одним входом, уравнение динамики которых имеет вид:

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

Обычно n m, так как при n n = d n /dt n . Это лишь другое обозначение операции дифференцирования. Обратная дифференцированию операция интегрирования записывается как 1/p. В операторной форме исходное дифференциальное уравнение записывается как алгебраическое:

Не надо путать эту форму записи с операционным исчислением хотя бы потому, что здесь используются непосредственно функции времени y(t), u(t) (оригиналы), а не их изображения Y(p), U(p), получаемые из оригиналов по формуле преобразования Лапласа. Вместе с тем при нулевых начальных условиях с точностью до обозначений записи действительно очень похожи. Это сходство лежит в природе дифференциальных уравнений. Поэтому некоторые правила операционного исчисления применимы к операторной форме записи уравнения динамики. Так оператор p можно рассматривать в качестве сомножителя без права перестановки, то есть py yp. Его можно выносить за скобки и т.п.

Поэтому уравнение динамики можно записать также в виде:

Дифференциальный оператор W(p) называют передаточной функцией. Она определяет отношение выходной величины звена к входной в каждый момент времени: W(p) = y(t)/u(t), поэтому ее еще называют динамическим коэффициентом усиления. В установившемся режиме d/dt = 0, то есть p = 0, поэтому передаточная функция превращается в коэффициент передачи звена K = b_m/a_n.

Знаменатель передаточной функции D(p) = a_op n + a₁p n — 1 + a₂p n — 2 + . + a_n называют характеристическим полиномом. Его корни, то есть значения p, при которых знаменатель D(p) обращается в ноль, а W(p) стремится к бесконечности, называются полюсами передаточной функции.

Числитель K(p) = b_op m + b₁p m — 1 + . + b_m называют операторным коэффициентом передачи. Его корни, при которых K(p) = 0 и W(p) = 0, называются нулями передаточной функции.

Звено САУ с известной передаточной функцией называется динамическим звеном. Оно изображается прямоугольником, внутри которого записывается выражение передаточной функции. То есть это обычное функциональное звено, функция которого задана математической зависимостью выходной величины от входной в динамическом режиме. Для звена с двумя входами и одним выходом должны быть записаны две передаточные функции по каждому из входов. Передаточная функция является основной характеристикой звена в динамическом режиме, из которой можно получить все остальные характеристики. Она определяется только параметрами системы и не зависит от входных и выходных величин. Например, одним из динамических звеньев является интегратор. Его передаточная функция W_и(p) = 1/p. Схема САУ, составленная из динамических звеньев, называется структурной.

6. Элементарные динамические звенья

Динамика большинства функциональных элементов САУ независимо от исполнения может быть описана одинаковыми по форме дифференциальными уравнениями не более второго порядка. Такие элементы называют элементарными динамическими звеньями. Передаточная функция элементарного звена в общем виде задается отношением двух полиномов не более чем второй степени:

W_э(p) = .

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

где p₁, p2, . p_n — корни полинома D(p). Аналогично