Переходный процесс в системе описывается уравнением:
Если все корни простые действительные ( ), то каждая из составляющих переходного процесса изменяется по закону экспоненты.
При этом переходная составляющая процесса стремится к нулю при (рис.1.2 а) или уходит в бесконечность от установившегося значения при (рис.1.2 б).
y y
0 t 0 t
а) б)
Рис. 1.2
Если в уравнении имеется комплексный сопряженный корень ,то соответствующая ему составляющая переходного процесса изменяется по синусоидальному закону (рис.1.3):
,
где
амплитуда колебаний;
частота колебаний;
сдвиг фаз (начальная фаза).
y y
i i
0 t 0 t
а) б)
Рис.1.3.
Эти величины играют роль постоянных интегрирования. При этом при (рис.1.3 а) и уходит в бесконечность при (рис.1.3. б).
АСР, переходные процессы в которых затухают с течением времени, называют устойчивыми.
АСР, у которых переходные процессы расходятся с течением времени, называют неустойчивыми.
Работоспособные системы должны быть устойчивыми.
Из изложенного вытекает необходимое условие устойчивости системы -отрицательность действительной части корней характеристического уравнения, соответствующего дифференциальному уравнению, которым описывается динамика системы.
1.3. Принцип суперпозиции. Типовые возмущения.
Для линейной системы справедлив принцип суперпозиции, заключающийся в том, что сумме любых возмущений соответствует сумма выходных реакций, каждая из которых определяется соответствующим воздействием; при любом изменении входного возмущения без изменения его формы выходная величина претерпевает такое же изменение, также не изменяя формы.
Принцип суперпозиции применим не только к суммам, но и к интегралам. Если входное возмущение в системе представляет собой бесконечно большое число бесконечно малых элементарных возмущений, то выходная величина линейной системы представляет собой сумму бесконечно малых реакций на эти бесконечно малые возмущения.
Принцип суперпозиции дает возможность выразить реакцию системы на любое возмущение через ее реакцию на определенный вид элементарных возмущений. Для этого достаточно представить произвольное возмущение элементарными воздействиями выбранного типа.
В качестве типовых возмущений чаще всего применяют единичную скачкообразную функцию, единичную импульсную функцию, единичную линейную функцию, единичное гармоническое колебание.
1. Единичная скачкообразная функция описывает мгновенное изменение какого-то воздействия от 0 до 1 (рис.1.4).
f
t
0 t
Рис.1.4.
Аналитически скачкообразную функцию записывают как :
0 при и
1 при и
2. Единичная импульсная функция описывает кратковременное возмущение, имеющее характер кратковременного импульсного толчка (рис. 1.5 ).
f
0 t
Рис.1.5.
Единичная импульсная функция, называемая - функцией Дирака, представляет собой первую производную от единичной скачкообразной функции:
и равна нулю везде, кроме , где она принимает бесконечное значение, причем при условии, что интеграл от нее по любому интервалу, содержащему , равен единице.
Функцию, обладающую такими свойствами, можно получить как предел положительного прямоугольного импульса, имеющего единичную площадь, когда длительность этого импульса стремится к нулю (рис. 1.6).
f
S(=1)
h
t
0
Рис.1.6.
3. Единичную линейную функцию при называют еще рамповым возмущением (Рис. 1.7).
f
0 t
Рис.1.7.
Такое возмущение является типичным для следящих систем регулирования.
4. Единичное гармоническое возмущение чаще всего записывают как функцию, изменяющуюся по синусоидальному закону ( Рис. 1.8):
f
1
0 t
Рис. 1.8
Такой тип возмущений применяют при частотных методах анализа АСР.
2. ПРИМЕНЕНИЕ ПРЕОБРАЗОВАНИЙ ЛАПЛАСА В ТАУ
2.1. Прямое и обратное преобразования Лапласа
Уравнения динамики АСР в символической Форме имеют вид:
,
или
.
Для нахождения интеграла функции нужно решить характеристическое уравнение:
Уравнения пятой и выше степеней в радикалах не решаются, а решения уравнений третьей и четвертой степеней громоздки. Поэтому при анализах АСР переходят от классических методов решения дифференциальных уравнений к их решениям с помощью преобразований Лапласа.
Преобразование Лапласа – функциональное преобразование, при котором функция вещественного переменного преобразуется в функцию комплексного
переменного .
Функция преобразуема по Лапласу, если она определена и однозначна для всей области и если
(конечен),
т.е. если она удовлетворяет условиям Дирихле.
Минимум значения , при котором указанный интеграл конечен, называют абсциссой абсолютной сходимости. Для большинства функций =0. Функцию называют оригиналом. Функцию - изображением (по Лапласу).
Преобразование оригинала в изображение называют прямым преобразованием Лапласа и осуществляют с помощью интеграла:
Преобразование изображения в оригинал называют обратным преобразованием Лапласа и осуществляют с помощью интеграла:
,
Операции прямого и обратного преобразований Лапласа будем обозначать через L и L-1 соответственно.
Связь между оригиналом и изображением записывается в виде:
.
Пример 1. Найти изображение единичной скачкообразной функции .
L L .
Пример 2. Найти изображение показательной функции .
L L ;
2.2. Основные теоремы преобразования Лапласа
2.2.1. Теоремы линейности
L L , ,
L L L .
2.2.2. Теоремы о конечном и начальном значениях оригинала
,
.
2.2.3. Теорема запаздывания в области вещественного переменного.
Если функция смещена на величину от начала координат (Рис. 2.1.), то такую функцию называют функцией с запаздывающим аргументом и записывают как .
f
0 t
Рис.2.1.
Изображение такой функции имеет вид:
L .
2.2.4. Теорема смещения в области комплексного переменного.
L ,
где - величина смещения в комплексной плоскости.
2.2.5. Теоремы масштабов
L ,
L ,
2.2.6. Теорема дифференцирования при нулевых начальных условиях
L L ,
где n-порядок производной.
2.2.7. Теорема интегрирования при нулевых начальных условиях
L ,
Здесь n-кратность интеграла.
2.2.8. Теорема свертки оригиналов
Если
и
то
,
где -переменная интегрирования.
2.3. Передаточная функция
Динамику АСР относительно регулируемой величины по каналу возмущающего воздействия можно записать в виде дифференциального уравнения в символической форме:
,
или
Если это уравнение преобразовать по Лапласу, используя теоремы линейности и дифференцирования при нулевых начальных условиях, то преобразованное по Лапласу уравнение по форме записи будет совпадать с символичной формой. Отличие заключается в том, что вместо оригиналов и нужно записывать их изображения и и под символом дифференцирования понимать комплексную переменную. С учетом сказанного, преобразованное по Лапласу при нулевых начальных условиях исходное уравнение динамики примет вид:
,
где
и -операторные многочлены левой и правой частей уравнения,
и -изображения выходной и входной величин.
Последнее выражение можно записать в виде
Отношение изображения выходной величины системы к изображению входного воздействия при нулевых начальных условиях называют передаточной функцией системы.
Обычно, эту функцию записывают как
2.4. Переходный процесс в АСР
Изображение регулируемой величины , т.е. изображение переходного процесса можно выразить через передаточную функцию системы относительно регулируемой величины по каналу возмущающего воздействия как:
,
где -изображение возмущающего воздействия.
Переходя к оригиналу, получим переходный процесс в виде:
L-1 L-1 .
Однако, производить непосредственное интегрирование затруднительно, поэтому для нахождения переходного процесса по его изображению для случаев, когда изображение переходного процесса является правильной рациональной дробью
,
следует применить теоремы Хевисайда (теоремы разложения).
Например, если корни характеристического уравнения действительные простые, отличные от нуля, то
,
где
- -ый корень характеристического уравнения ,
,
3. ДИНАМИЧЕСКИЕ ХАРАКТЕРИСТИКИ АСР
3.1. Комплексная частотная функция
Предположим, что АСР описывается уравнением:
,
где
- регулируемая величина;
- приложенное к системе возмущающее воздействие.
Допустим, что входное воздействие – синусоидальное колебательное, описываемое уравнением:
,
где
-амплитуда колебаний;
-частота колебаний, .
В этом случае, в соответствии с принципом суперпозиции, спустя некоторое время после приложения входного воздействия на выходе системы установятся синусоидальные колебания с такой же частотой, но отличные по модулю и фазе.
где
амплитуда;
сдвиг фаз.
Продифференцируем входную и выходную функции:
Подставив далее значения функций и их производных в исходное уравнение и разделив выходную величину на входное воздействие, получим:
Величину называют комплексной частотной функцией.
Комплексная частотная функция является отношением вынужденного движения на выходе системы к входному гармоническому воздействию.
Сравнивая выражение комплексной частотной функции с передаточной функцией, можно увидеть, что комплексная частотная функция является частным случаем передаточной функции, когда комплексная переменная является чисто мнимой величиной .
Обозначив комплексную частотную функцию через , можно записать:
Комплексная частотная функция представляет собой вектор в комплексной плоскости (Рис.3.1)
0
Рис. 3.1
Вектор характеризуется амплитудой (А) и фазой ( ). Амплитудой (модулем) вектора является его длина. Фаза вектора – угол между положительной вещественной полуосью и направлением вектора.
Если изменять частоту входных колебаний от 0 до , то получим ряд векторов в комплексной плоскости, концы которых опишут некоторую кривую, называемую годографом. Годограф вектора комплексной частотной функции при изменяющейся частоте в пределах , называют амплитудно-фазовой частотной характеристикой (АФХ).
В соответствии с рисунком 3.1 различают частотные характеристики:
1. Вещественная частотная характеристика - проекция вектора на вещественную полуось.
2. Мнимая частотная характеристика - проекция вектора на мнимую полуось.
3. Амплитудно-частотная характеристика .
4. Фазо-частотная характеристика .
Между приведенными частотными характеристиками имеются однозначные зависимости:
3.2 Переходная функция
Переходной функцией называют изменение во времени выходной величины системы в результате приложения к ней единичного скачкообразного входного воздействия.
Обозначив переходную функцию через и соответствующее ей изображение через , можно записать:
;
.
3.3 Весовая функция
Весовой функцией называют реакцию системы на импульсную возмущающую функцию .
Так как , то , т.е. изображение импульсной функции равно единице.
Тогда , т.е. изображение весовой функции равно передаточной функции. Весовая функция записывается как:
.
4. ТИПОВЫЕ ЗВЕНЬЯ АСР
4.1. Введение
Любая АСР состоит из реальных физических элементов, отличающихся друг от друга своими физическими свойствами. Однако, многие из них, будучи различными по своей физической сущности, могут быть описаны одинаковыми математическими уравнениями. Это дает возможность классифицировать объекты и отдельные элементы системы в целом по их математическому описанию.
4.2. Безынерционное звено
К таким звеньям относят физические элементы, описываемые уравнением:
.
Коэффициент пропорциональности называют коэффициентом передачи.
Его размерность определяется размерностями выходной и входной величин:
.
Преобразовав уравнение по Лапласу:
,
получим
,
т.е. передаточная функция звена равна его коэффициенту передачи.
Амплитудно-фазовая характеристика звена:
представляет собой точку на вещественной полуоси комплексной плоскости на расстоянии от начала координат (Рис. 4.1).
K
0
Рис.4.1.
Переходная функция:
представляет собой скачкообразную функцию высотою (Рис.4.2)
h
k
0 t
Рис. 4.2.
Весовая функция представляет собой импульсную функцию (Рис.4.3).
V
0 t
Рис.4.3.
4.3. Идеальное дифференцирующее звено
К ним относят элементы АСР, у которых выходная величина пропорциональна производной входного воздействия:
,
где Т - коэффициент пропорциональности, имеющий размерность времени и
называемый постоянной дифференцирования.
Преобразовав уравнение по Лапласу, получим:
и .
АФХ звена
является прямой, совпадающей с верхней мнимой полуосью комплексной плоскости (Рис.4.4).
)
900
0
Рис.4.4.
Переходная функция звена:
представляет собой импульсную функцию (Рис.4.5).
h
0 t
Рис.4.5.
Весовая функция звена является импульсной функцией второго рода (Рис.4.6).
V
0 t
Рис.4.6.
4.4. Идеальное интегрирующее звено
Такими звеньями называют элементы АСР, у которых выходная величина пропорциональна интегралу от входного воздействия:
.
Это уравнение можно записать и в виде:
.
Отсюда вытекает физический смысл коэффициента пропорциональности:
.
Коэффициент пропорциональности является скоростью изменения выходной величины, приведенной к единице возмущающего воздействия. В связи с этим величину называют приведенной скоростью.
Преобразовав исходное уравнение по Лапласу, получим:
и .
АФХ звена
является прямой, совпадающей с нижней мнимой полуосью комплексной плоскости (Рис. 4.7).
0 900
ω
0
Рис.4.7.
Переходная функция звена
является прямой линией с угловым коэффициентом (Рис. 4.8)
h
t
0
Рис.4.8.
Весовая функция звена представляет собой скачкообразную функцию высотой (Рис.4.9).
V
0 t
Рис.4.9.
4.5. Апериодическое звено первого порядка
К таким звеньям относят физические элементы, описываемые дифференциальным уравнением первого порядка:
,
где
-коэффициент передачи;
-постоянная времени.
Преобразовав уравнение по Лапласу:
,
получим передаточную функцию звена:
.
АФХ звена
,
где
, .
Амплитудно-фазовая характеристика представляет собой полуокружность , расположенную в четвертом квадранте комплексной плоскости (Рис.4.10)
K
0 900
Рис.4.10.
Переходная функция звена :
,
Так как
= ,
то
.
График переходной функции приведен на рисунке 4.11 а. Она представляет собой возрастающую экспоненту.
h v
K
0 t 0 t
а) б)
Рис.4.11.
Весовая функция звена является положительной убывающей по модулю экспонентой (Рис.4.11б).
4.6. Реальное дифференцирующее звено
К таким звеньям относят элементы АСР, описываемые дифференциальным уравнением:
.
Преобразовав его по Лапласу
,
получим передаточную функцию
.
Отсюда
,
где
,
.
При , . При , .
АФХ является полуокружностью, расположенной в первом квадранте комплексной плоскости (Рис.4.12).
900
0
K
Рис.4.12.
Переходная функция звена
является положительной убывающей экспонентой (Рис. 4.13а).
h V
0 t
0 t
а) б)
Рис4.13.
Весовая функция звена представляет собой отрицательную убывающую по модулю экспоненту (Рис.4.13б).
4.7. Инерционные звенья второго порядка
К таким звеньям относят элементы АСР, описываемые дифференциальным уравнением второго порядка:
,
где
-коэффициент передачи;
, -постоянные времени.
Преобразуем уравнение по Лапласу:
Отсюда
.
Корни характеристического уравнения
,
т.е. корни многочлена в знаменателе передаточной функции, выражаются как:
.
Обозначим , .
Тогда
и
Рассмотрим случаи:
1. Если , то корни характеристического уравнения действительные простые и переходная функция
Так как
то
и
,
где
; .
Производные функции
;
.
При , , .
Переходная функция при асимптотически стремится к значению , т.е. к установившемуся значению.
При
.
.
.
В соответствии со значениями при и графики функции и ее
производных принимают вид, приведенный на рисунке 4.14.
h(t) Точку на кривой ,
соответствующую
максимальной первой
K производной и переходу через
нуль второй производной,
0 t называют точкой перегиба
(V(t)) кривой.
0 t
K
T2
0 t
Рис.4.14.
2. Если , то корни характеристического уравнения действительные равные .
Следовательно,
.
И в этом случае , кривая функции имеет вид, изображенный на рисунке 4.14.Следовательно, при процесс апериодический.
3. Если , то корни комплексные сопряженные.
Тогда
и
Переходный процесс колебательный, затухающий по экспоненциальному закону.
В связи с этим, элементы АСР, описываемые уравнением второго порядка, называют колебательными звеньями в случае и апериодическими в случаях .
Аналитическое выражение АФХ звеньев второго порядка имеет вид:
,
где
;
.
Графики АФХ апериодического и колебательного звеньев приведены на
рисунке 4.15 а,б.
К К
1800 1800
рез
а) апериодическое б) колебательное
Рис.4.15.
АФХ колебательного звена (рис. 4.15б) имеет промежуточный экстремум. Частоту, соответствующую этому экстремуму (максимуму), называют резонансной.
4.8. Общее свойство типовых звеньев
В общем случае передаточная функция звена может быть записана в виде правильной рациональной дроби.
Раскладывая многочлены и на множители, можно записать:
где
- общий коэффициент передачи звеньев;
- -ые корни многочленов и .
Корни многочлена числителя передаточной функции называют ее нулями, корни многочлена знаменателя передаточной функции называют полюсами.
Так как система устойчива, то полюсы ее передаточной функции должны располагаться в левой полуплоскости, нули же могут располагаться и в левой и в правой полуплоскостях комплексной плоскости, например рисунок 4.16.
0
Рис.4.16.
Из рисунка 4.16 видно, что фазы векторов АФХ меньше для случаев, когда нули располагаются в левой полуплоскости, по сравнению со случаем, когда нуль располагается в правой полуплоскости.
Звенья, у которых все полюсы и нули передаточной функции располагаются в левой полуплоскости (или на мнимой оси), называют минимально-фазовыми.
Рассмотренные типовые звенья являются минимально-фазовыми.
4.9. Звено транспортного (чистого) запаздывания
В ряде технологических процессов имеет место отставание во времени появления выходного сигнала по сравнению с моментом приложения входного воздействия.
Классическим примером такого явления служат транспортные средства в виде различных транспортеров.
Уравнение, описывающее это явление, имеет вид:
,
где -величина запаздывания выходного сигнала по сравнению с входным воздействием.
Преобразовав уравнение по Лапласу, получим:
АФХ является окружностью единичного радиуса с центром в начале координат комплексной плоскости (рис. 4.17).
1
0
Рис.4.17.
Переходная функция звена запаздывания
является ступенчатой функцией, смещенной на величину запаздывания от начала координат (рис. 4.18).
h
1
0 t
Рис.4.18.
Весовая функция представляет собой импульсную функцию, смещенную на величину запаздывания от начала координат (рис. 4.19).
V
0 t
Рис.4.19.
Звено транспортного запаздывания является неминимально-фазовым, его передаточная функция имеет бесконечное число нулей и полюсов, расположенных как в левой, так и в правой полуплоскостях комплексной плоскости.
Однако, принципиально возможно осуществить приближение передаточной функции передаточными функциями минимально-фазового вида. Необходимость этого возникает при моделировании динамических свойств системы. Такое приближение можно осуществить, используя метод разложения функции в быстро сходящиеся ряды при ограничении сравнительно небольшим числом членов ряда.
Наиболее универсальным методом является приближение с помощью дробно-рациональной функции.
Разложение функции в дробный ряд Пада является примером такого разложения:
,
где и - функции комплекса и ν , причем
Например, при
.
-
Типовые соединения звеньев
4.10.1. Введение
В структурных схемах АСР отдельные ее элементы (звенья) определенным образом взаимодействуют между собой, получая и выдавая соответствующие сигналы.
Возможны три типа взаимодействия, иначе говоря, три типа соединения звеньев.
4.10.2. Последовательное соединение звеньев
Последовательным соединением называют тот случай взаимодействия, когда выходной сигнал предыдущего звена является входной величиной последующего
(рис. 4.20).
f(t) y1(t) y2(t) yn-2(t) yn-1(t) yn(t)
1 2 … n-1 n
Рис.4.20.
Обозначив передаточные функции звеньев через в соответствии с рисунком 4.20 находим общую передаточную функцию соединения звеньев.
Передаточная функция последовательного соединения звеньев равна произведению передаточных функций отдельных звеньев, входящих в соединение.
4.10.3. Параллельное согласное соединение звеньев
Таким соединением называют случай взаимодействия элементов АСР, когда имеется один вход, а выходы алгебраически суммируются (рис. 4.21).
y1(t)
1
y2(t)
2
f(t) y(t)=∑yi(t)
yn-1(t)
n-1
yn(t)
n
Рис.4.21.
В соответствии с рисунком 4.21 имеем:
Передаточная функция параллельного согласного соединения звеньев равна сумме передаточных функций звеньев, входящих в соединение.
4.10.4. Параллельное встречное соединение звеньев.
Такое соединение звеньев чаще называют соединением по типу обратной связи. Это такое соединение двух групп звеньев, когда сигнал первой группы звеньев передается на вход этой же группы, соответствующим образом измененный другой группой звеньев (рис. 4.22).
f(t) + z(t) y(t)
Wпр(p)
± yос(t)
Wос(p)
Рис.4.22.
В соответствии с рисунком 4.22 имеем:
Y(p)=Wпр(p)Z(p),
Z(p)=F(p)±Yос(p),
Yoc(p)=Woc(p)Y(p).
Подставив второе и третье уравнения в первое, получим:
Y(p)=Wпр(p)[F(p)± Woc(p)Y(p)]= Wпр(p)F(p)± Wпр(p )Woc(p)Y(p);
Y(p) Wпр(p )Woc(p)Y(p)=Wпр(p)F(p).
Отсюда
Передаточная функция соединения звеньев по типу обратной связи определяется дробью, в числителе которой записывается передаточная функция прямой цепи ( ) передачи сигнала, а в знаменателе – единица, алгебраически просуммированная с произведением передаточных функций прямой цепи и цепи обратной связи ( ).
Знак ″-″ в знаменателе относится к случаю положительной обратной связи, т.е. когда сигнал от обратной связи совпадает по знаку с основным входным сигналом.
Знак ″+″ относится к отрицательной обратной связи, т.е. к случаю, когда сигнал от обратной связи и основной входной сигнал противоположны по знаку.
В системах регулирования для обеспечения устойчивости их работы применяют отрицательную обратную связь.
В таких случаях:
.
Положительная обратная связь применяется в качестве местных обратных связей в различного рода усилителях с целью повышения коэффициента усиления последних.
Передаточную функцию прямой цепи передачи сигналов можно записать в виде:
где
Кпр - коэффициент передачи прямой цепи,
- передаточная функция элементов прямой цепи при единичном
коэффициенте передачи.
Тогда
.
При достаточно большом значении коэффициента передачи прямой цепи (теоретически при Кпр ) получим:
.
Контур с обратной связью при достаточно большом коэффициенте передачи прямой цепи называют предельной системой.
Динамические свойства предельной системы определяются динамическими свойствами обратной связи.
5. Основные свойства типовых объектов регулирования
5.1 Статические и астатические объекты регулирования
С учетом явления запаздывания передачи входного сигнала уравнение, описывающее динамику одномерного объекта, можно записать в виде:
,
где
-регулируемая величина;
-основное возмущающее воздействие (нагрузка объекта);
-величина запаздывания.
Преобразовав уравнение по Лапласу, получим:
.
Запишем уравнение в виде:
или
где
-коэффициент передачи объекта;
-постоянные времени.
Левый многочлен уравнения можно записать в виде сомножителей первой степени.
Тогда
и
.
Следовательно, в этом случае объект структурно может быть представлен в виде последовательного соединения звена транспортного запаздывания и апериодических звеньев первого порядка (рис. 5.1).
q(t) y(t)
W(p) W1(p) W2(p) … Wn-1(p) Wn(p)
Рис.5.1.
Объект статичен, так как в установившемся режиме имеется однозначная зависимость между выходным сигналом и входным воздействием
Если в левой части исходного уравнения отсутствует свободный член , то в установившемся режиме не будет однозначной зависимости между выходом и входом, т.е. объект астатический. В этом случае
.
Это уравнение можно записать в виде:
или
где - постоянные времени;
Отсюда
.
Следовательно, астатический объект регулирования может быть структурно представлен в виде последовательного соединения звена транспортного запаздывания, интегрирующего звена и апериодических звеньев первого порядка
(рис. 5.2).
q(t) y(t)
W(p) Wu(p) W1(p) … Wn(p)
Рис.5.2.
Подставив в передаточные функции , получим аналитические выражения АФХ объектов.
Графики АФХ объектов регулирования приведены на рисунке 5.3 а,б.
Коб
0
0 0
а) статический ОР б) астатический ОР
Рис.5.3.
Переходные характеристики статического и астатического объектов приведены на рисунке 5.4 а,б.
h h
0 t 0 t
а) статический ОР б) астатический ОР
Рис.5.4.
-
Экспериментальное определение временных динамических характеристик объектов регулирования
Для получения временной динамической характеристики объекта регулирования организуют эксперимент.
В какой-то момент времени к объекту регулирования, находящемуся в установившемся режиме работы, прикладывают некоторое возмущающее воздействие и затем регистрируют отклонения регулируемой величины во времени до ее восстановившегося значения с помощью соответствующего измерительного прибора.
Перед началом эксперимента нужно убедиться в том, что объект регулирования находится в установившемся режиме работы. Для достоверности измеренных отклонений регулируемой величины необходимо один и тот же опыт повторять 2-3 раза.
В найденную путем эксперимента динамическую характеристику объекта включены характеристика собственно объекта и характеристика регистрирующего прибора.
В связи с этим желательно в качестве датчиков отклонений физических величин от их номинальных значений применять те, которые входят в измерительный блок регулирующего устройства, поддерживающего на каком-то значении данную физическую величину.
В качестве возмущающих воздействий могут быть использованы типовые функции, некоторые из которых приведены в таблице 5.1.
Таблица.5.1.
№п/п Возмущение
f(t) График
f(t) Изображение
F(p) 1.
2.
3.
4.
Ступенчатое при
высоте ступени h
П-образное при высоте h и при длительности
Ступенчатое при высоте h и времени нанесе-ния возмущения
Треугольное при высоте h и времени нанесе-ния возмущения
2 f
h
0 t
f
h
0 t
f h
0
t
f
h
0 t
Наиболее часто при определении динамических свойств объекта регулирования применяют ступенчатое возмущающее воздействие (строка 1 табл.5.1), т.е. экспериментально находят переходную функцию объекта.
5.3. Определение передаточной функции объекта регулирования по его экспериментальной переходной функции
-
Статический объект регулирования
5.3.1.1. Структурное представление статического объекта последовательным соединением апериодического звена первого порядка и звена запаздывания
В общем виде переходный процесс в объекте, как реакция на ступенчатое возмущение, приведен на рисунке 5.5 а.
Для удобства использования этого процесса при нахождении передаточной функции объекта, его (процесс) представляют приведенным к единичному ступенчатому возмущению, т.е. переходной характеристикой, как реакцией на единичное ступенчатое возмущение (рис. 5.5 б).
f
h
0 t
y h(t)
y(t)
= A
0 t 0 t
а) б)
Рис.5.5.
При представлении объекта простейшей моделью в виде последовательного соединения звена запаздывания и апериодического первого порядка передаточная функция запишется как:
.
Параметры передаточной функции и об=0+е находятся из графика переходной функции, как показано на рисунке 5.5 б.
5.3.1.2. Структурное представление статического объекта последовательным соединением звена запаздывания и апериодического звена второго порядка
При такой модели объекта
.
Значения Коб и 0 (чистое запаздывание) определяются сразу же из графика переходной функции (рис. 5.6).
Для определения коэффициентов передаточной функции Т1 и Т22 на графике переходной функции находят точку перегиба (т. А) и определяют ее координаты (0+n),
;проводят касательную к кривой в точке перегиба; определяют численные значения Т1, Тоб и рассчитывают F1, F2, F=F1+F2, причем площадями являются:
F1 –площадь, ограниченная кривой , ее асимптотой h(∞) и прямой, проведенной
через точку перегиба параллельно оси ординат;
F2 –площадь, ограниченная кривой , ее асимптотой h(∞) и прямыми,
проведенными параллельно оси ординат на расстояниях t=0 и t=0+n от начала
координат;
h
Tоб
F2 T1 F1 F1+F2=F
B C
h(t)
A
hп
0 t
Рис.5.6.
После произведенных построений находят
или .
Коэффициент Т1 может быть определен, как длина (в размерности времени) проекции касательной к переходной характеристике в точке перегиба на асимптоту кривой .
5.3.1.3. Определение передаточной функции статического объекта по его экспериментальной переходной характеристике методом Симою (методом площадей)
Статический объект регулирования в общем случае может быть описан линейным дифференциальным уравнением с постоянными коэффициентами или передаточной функцией вида
которая с желаемой степенью точности может быть аппроксимирована дробно-рациональной функцией
Из этих уравнений имеем
.
Отсюда
или
где
Здесь - нормированная переходная функция.
Введя безразмерное время , и заменив переменную, получим:
.
Это выражение позволяет вычислить любой коэффициент .
Например, при к=1
при к=2
при к=3
Из этих уравнений вытекает геометрический смысл коэффициентов F1, F2, …,Fк.
Коэффициент F1 является площадью, ограниченной нормированной кривой переходной функции, осью ординат и линией установившегося значения функции
(рис. 5.7).
F1
F1 F2
1 F1(∞)
0 t 0 t
Рис.5.7. Рис.5.8.
Коэффициент F2 является площадью, ограниченной функцией F1(t), осью ординат и линией установившегося значения функции F1(∞) (рис. 5.8).
Коэффициент F3 – площадь, ограниченная функцией F2(t),осью ординат и линией установившегося значения функции F2(∞) (рис. 5.9.).
F2
F3
F2(∞)
0 t
Рис.5.9.
Коэффициент Fк – площадь, ограниченная функцией Fк-1(t), осью ординат и линией установившегося значения функции Fк-1(∞).
Иначе говоря, коэффициенты Fк являются интегралами к-ой кратности функции . Из геометрических представлений коэффициентов Fк вытекает простой метод нахождения безразмерной передаточной функции объекта, заключающийся в последовательности выполнения операций:
-
Ось абсцисс переходной функции разбивают на отрезки t, исходя из условия, что на протяжении всего графика функция мало отличается от прямой в пределах двух t (рис. 5.10).
h
0 t
Рис.5.10.
2. Функцию н приводят к нормированному (тарированному) виду путем деления (t ) на установившееся значение .
.
-
Заполняют столбцы 1,2,3 таблицы 5.2.
Таблица 5.2.
1234Xxxxxxxxxxxx
Определяют площадь F1
где - сумма третьего столбца таблицы 5.2.
-
Заполняют столбец 4 таблицы 5.2 и строят функцию в масштабе безразмерного времени (рис. 5.11).
1
0
Рис.5.11.
-
Предполагая, что можно ограничиться тремя коэффициентами знаменателя передаточной функции, заполняют таблицу 5.3 и определяют площади F2 и F3.
Таблица.5.3.
123456XxxxxxXxxxxx
где
, - суммы элементов четвертого
и шестого столбцов таблицы 5.3.
В этом случае
.
6. Находят размерную передаточную функцию объекта в виде:
где - коэффициент передачи объекта.
7. Если в объекте имеется чистое запаздывание 0, то передаточная функция объекта записывается как:
-
Астатический объект регулирования
5.3.2.1. Структурное представление объекта регулирования, обладающего интегрирующими свойствами, последовательным соединением звена запаздывания и интегрирующего звена.
В наипростейшем виде структурная модель объекта представляется последовательным соединением звена запаздывания и интегрирующего звена.
Такому представлению структурной модели объекта соответствуют последовательные участки переходной характеристики на рисунке 5.12: отрезок прямой АВС и асимптота СД.
h
Д
C
A B
0 t
Рис.5.12.
В таком случае передаточная функция объекта запишется как:
,
где
об=0+е – общее запаздывание, равное сумме транспортного 0 и емкостного е
запаздываний;
Коб – коэффициент передачи объекта.
Коэффициентом передачи объекта регулирования называют отношение установившегося значения выходной величины к установившемуся значению входного воздействия.
В данном случае под коэффициентом передачи объекта понимают отношение установившегося значения скорости изменения выходной величины к установившемуся значению входного воздействия .
В соответствии с этим:
.
Для единичного скачкообразного входного воздействия =1, а установившееся значение скорости изменения выходной величины, в соответствии с рис. 5.12, есть отношение приращения h ординаты асимптоты к соответствующему приращению t оси абсцисс.
Таким образом, численное значение Коб для астатического объекта регулирования находится как угловой коэффициент асимптоты кривой h(t), т.е. как тангенс угла наклона асимптоты к оси времени, например
.
-
Структурное представление объекта регулирования, обладающего интегрирующими свойствами, последовательным соединением звена транспортного запаздывания, апериодического звена первого порядка и интегрирующего звена.
Более точный способ аппроксимации динамических свойств астатического объекта вытекает из представления структурной модели объекта в виде последовательного соединения трех типовых звеньев: звена транспортного запаздывания, апериодического первого порядка и интегрирующего.
Такому представлению структурной модели объекта соответствуют последовательные участки переходной характеристики h(t) на рисунке 5.13: отрезок прямой АВ, кривая ВЕ, отрезок прямой ЕД, являющийся частью асимптоты ДС.
h
F Д
h(∞)
h(t)
h1(t) E h1(∞)
α
A B K
0 t
C
Рис.5.13.
В этом случае
где Коб – коэффициент передачи, определяемый, как и ранее, через коэффициент
наклона асимптоты
0 – транспортное запаздывание;
Тоб – постоянная времени апериодического звена, определяемая из рисунка 5.13
как .
-
Структурное представление объекта регулирования, обладающего интегрирующими свойствами, параллельным соединением интегрирующего и статического звеньев, последовательно соединенным со звеном транспортного запаздывания
В частных случаях технологических процессов кривые переходных функций могут быть знакопеременными от момента приложения воздействия до установившегося значения.
Например, переходные характеристики относительно уровня воды в барабане парового котла, отражающие явления уплотнения или набухания объема котловой воды в барабане котла при изменениях его нагрузки по пару (рис. 5.14).
h h
0 t
0 t
а) уменьшение б)увеличение
расхода пара расхода пара
Рис.5.14.
Такие кривые переходных функций могут быть представлены двумя функциями:
прямой hи (t) путем переноса асимптоты кривой в точку t=0 и кривой hст(t), получаемой вычитанием исходной характеристики h(t) из прямой hи (t), например рисунок 5.15.
h
hu(t)
h(t)
hcт(t)
α
0 t
Рис.5.15.
В соответствии с этим объект структурно может быть представлен в виде схемы, приведенной на рисунке 5.16.
+
f(t) h(t)
-
Рис.5.16.
Измерив на исходном графике переходной функции величину запаздывания 0, можно записать передаточную функцию, соответствующую явлению чистого запаздывания.
.
Передаточную функцию интегрирующего звена находят, как и в предыдущих случаях, по тангенсу угла наклона асимптоты кривой h(t)
.
Передаточную функцию статического звена находят из кривой hст(t) одним из способов нахождения передаточной функции статических объектов.
При этом
Коб.ст (при .
Тогда
.
6. Основные линейные законы регулирования и их реализация в промышленных регуляторах
|