- •Введение
- •Использование преобразование Лапласа для исследования динамических систем
- •Раздел 1. Элементарные звенья систем управления
- •Понятие звена
- •Усилительное звено
- •1.3. Интегрирующее звено
- •1.4. Апериодическое звено
- •1.5. Колебательное звено
- •1.6. Дифференцирующее звено
- •1.7. Звено запаздывания
- •1.8. Частотные характеристики типовых звеньев
- •1.8.1. Частотные характеристики усилительного звена
- •1.8.2. Частотные характеристики интегрирующего звена
- •1.8.3. Частотные характеристики апериодического звена
- •1.8.4. Частотные характеристики колебательного звена
- •1.8.5. Частотные характеристики дифференцирующего звена
- •1.8.6. Частотные характеристики звена запаздывания
- •1.9. Характеристики элементарных звеньев в табличной форме
- •Раздел 2. Линейный анализ и синтез систем управления с использованием математической лаборатории matlab
- •2.1 Анализ и синтез динамической системы в пространстве состояний
- •2.2 Устойчивость систем 1го, 2го порядка и более высокий порядок
- •2.3 Инерционная динамическая система 2-го порядка
- •2.4 Синтез инерционных систем первого и второго порядка
- •2.4.1 Выбор коэффициента усиления в инерционной системе I порядка.
- •2.4.2 Выбор коэффициентов усиления и гибкой обратной связи для инерционной системы второго порядка.
- •2.4.3 Применение методов аналитического синтеза для структуризации систем управления.
- •2.5 Исследование устойчивости для систем порядка n
- •2.5.1 Корневой анализ системы порядка n
- •2.6 Метод синтеза системы управления по желаемому значению частот и коэффициентов относительного демпфирования форм движения изолированной системы
- •2.6.1 Требование к распределению частот и относительных коэффициентов демпфирования изолированных форм движения
- •2.6.2 Выбор коэффициентов усиления по заданным требованиям к собственным значениям, частотам и коэффициентам относительного демпфирования в замкнутой системе
- •2.6.3 Математическая модель в пространстве состояний
- •2.6.4 Анализ характеристик замкнутой системы управления
- •2.7 Аналитическое конструирование оптимального регулятора (акор)
- •3. Алгебраические критерии устойчивости Рауса и Гурвица.
- •4. Частотные критерии анализа устойчивости и синтеза систем уравнения
- •4.1 Метод частотного анализа устойчивости Михайлова
- •4.2 Критерий устойчивости Найквиста-Михайлова и запасы устойчивости по амплитуде и по фазе для систем с единичной отрицательной обратной связью
- •4.3 Выбор коэффициентов усиления к в прямой и обратной цепи по заданному запасу устойчивости.
- •4.4 Метод наращивания устойчивых контуров. Выбор коэффициентов усиления позиционно-дифференцирующего и позиционно-интегрально-дифференцирующего контроллера
- •4.5 Расчет позиционного интегрально-дифференцирующего контроллера(pid)
- •5. Синтез корректирующих звеньев в прямой цепи
- •5.1 Расчет параметров передаточной функции фильтра средних частот (фсч)
- •5.2 Расчет параметров фильтра низких частот (фнч)
- •5.3 Расчет параметров фильтра высоких частот (фвч)
- •6. Синтез корректирующего звена в цепи обратной связи.
- •Процедура улучшения решения за счёт выбора параметров во втором приближении.
- •Раздел 7. Аналитический обзор методов исследования устойчивости систем автоматического управления
- •7.1. Критерии устойчивости [1, 7, 11, 12].
- •7.2. Частотные критерии устойчивости [7, 8, 11].
- •7.3. Запас устойчивости систем [7].
- •7.4. Точность систем [8].
- •7.5. Качество систем [1, 2, 8, 12].
- •7.6. Случайные процессы в системах [8].
- •Литература
1.5. Колебательное звено
Колебательное звено имеет дифференциальное уравнение
.
(13)
Передаточная функция звена
.
(14)
Характер переходного процесса звена или соединения, определяемого дифференциальным уравнением (13), зависит от расположения корней его характеристического уравнения
(15)
на комплексной плоскости.
Корни характеристического уравнения (15)
(16)
С учетом (14), (11) и таблицы преобразований Лапласа, находим изображение выходной величины:
.
(17)
В зависимости от знака подкоренного выражения (16) при нахождении оригинала по его изображению (17) могут возникнуть три случая:
1. При T1/T2 > 2 оба корня характеристического уравнения вещественные отрицательные: р1 = - α1, р2 = - α2. С учетом этого запишем выражение (17) в виде
.
По этому изображению согласно таблицы преобразований Лапласа находим оригинал:
.
(18)
Таким образом, при T1/T2 > 2 переходный процесс определяется двумя экспонентами и в этом случае дифференциальное уравнение (13) характеризует переходные процессы соединения, состоящего из двух соединенных последовательно апериодических звеньев. Это видно также непосредственно из передаточной функции соединения, если ее записать в виде
или
,
где T3=1/α1 и T4=1/ α2.
Следовательно, при T1/T2 > 2 нет необходимости вводить понятия нового типового звена, хотя на практике часто такое соединение называют инерционным звеном второго порядка.
2. При T1/T2 = 2 характеристическое уравнение имеет два одинаковых вещественных отрицательных корня
p1 = p2 = - α = - 1/T2.
С учетом этого запишем выражение (17) в виде
.
По таблице преобразования Лапласа находим:
.
(19)
Переходный процесс периодический. Так как при этом передаточная функция (14) может быть представлена в виде
,
где T = 1/α, то при T1/T2 = 2, так же как и при T1/T2 > 2, нет необходимости вводить понятия нового типового звена.
3. При T1/T2 < 2 характеристическое уравнение имеет два сопряженных комплексных корня
,
где
;
.
(20)
С учетом этого запишем выражение (17) в виде
.
(21)
Обозначив в (21)
и
,
найдем оригиналы:
и
.
Находим характер изменения выходной величины звена:
,
или
.
(22)
Таким
образом, переходный процесс звена при
T1/T2
< 2,
характеризуемый уравнением (22), периодичен
и представляет собой затухающую
синусоиду, амплитуда которой убывает
от полупериода к полупериоду по
экспоненциальному закону
.
В этом случае звено нельзя представить
в виде соединения из других звеньев. В
связи с этим элементарное звено,
динамические качества которого
определяются дифференциальным
уравнением (13), при T1/T2
< 2
относится к типовым звеньям и называется
колебательным звеном. Переходные
процессы колебательного звена в
зависимости от отношения T1/T2
представлены на рис. 7.
Рис. 7. Передаточная функция и переходные процессы колеба-тельного звена при различных значениях отношения постоянных времени.
Как следует из выражения (22), мнимая составляющая ω корней характеристического уравнения является круговой частотой колебательного звена. Период колебаний Т = 2π/ω. Оценкой переходного процесса колебательного звена служит степень затухания колебаний. Степенью затухания ψ называется отношение разности двух соседних амплитуд одного знака (взятых относительно среднего положения kx0вх) к первой из них (рис. 8,а):
.
(23)
Как следует из рис. 8,а,
,
.
Так как t2 – t1 = T, то подставив значения A1 и A2 в (23), получим:
.
(24)
Чем ближе к единице величина ψ, тем быстрее затухают колебания переходного процесса.
Рис. 8. Переходные процессы колебательного звена при 1>ψ>0 и ψ = 0.
Степень затухания зависит от отношения вещественной составляющей комплексных корней характеристического уравнения α к их мнимой составляющей ω. В свою очередь это отношение определяется отношением постоянных времени T1/T2:
.
4. При T1 = 0 T1/T2 = 0 вещественная и мнимая составляющие корней характеристического уравнения будут равны:
;
.
Подставив эти значения в выражение (22) для переходного процесса колебательного звена, получим
.
(25)
Такое колебательное звено называется консервативным.
Переходный процесс будет в этом случае незатухающим колебательным (так как ψ = 0) с частотой ω0 = 1/T2, периодом T = 2πT2 и амплитудой A=kx0вх (рис.8,б).
Чем больше T1 и меньше Т2, тем больше степень затухания колебательного звена.
Следовательно, для уменьшения колебательности систем регулирования в колебательных звеньях необходимо увеличивать постоянную времени T1 и уменьшать T2. Однако это целесообразно делать лишь и определенных пределах, так как при чрезмерном увеличении отношения T1/T2 переходный процесс затягивается (см. рис. 7) и время регулирования увеличивается.
На рис. 9 даны примеры колебательных звеньев.
Входной величиной мембранного пневматического клапана (рис. 9,а) является давление ΔРвх, а выходной – перемещение ΔSвых штока клапана (отсчет ведется в малых приращениях от равновесного состояния).
Рис. 9. Примеры колебательных звеньев.
Если нельзя пренебречь инерцией подвижной системы клапана и силами трения, то условие равновесия сил, действующих на клапан, запишется как
.
Входное усилие при площади F мембраны равно:
.
Сила инерции fи равна произведению массы m подвижной системы на ускорение a = d2(ΔSвых)/dt2:
.
Учитывая только силу вязкого трения, которая пропорциональна скорости перемещения подвижной системы, получим:
.
Сила противодействия пружины пропорциональна ее сжатию
,
где с - жесткость пружины.
Подставив значения сил в уравнение равновесия, получим:
.
В настоящее время принято составлять дифференциальные уравнения звеньев в безразмерных (относительных) единицах.
Безразмерной единицей давления будем считать отношение ΔРвх к максимальной величине давления Рмакс на мембрану, при котором клапан полностью закрывается; безразмерной единицей перемещения штока клапана примем отношение ΔSвых к полному ходу Sмакс
;
,
откуда
;
;
.
Подставив эти значения в дифференциальное уравнение, получим выражение его в безразмерных единицах:
.
С учетом того, что сSмакс = РмаксF, можно записать:
.
Таким
образом, при учете инерции подвижной
системы и вязкого трения мембранный
пневматический клапан при b/
<2
является колебательным звеном.
Постоянные времени и коэффициент передачи его равны:
;
;
.
Из этого примера следует, что в элементах систем регулирования вязкое трение не всегда является нежелательным. В данном случае достаточно высокое вязкое трение обеспечивает устойчивую работу клапана, так как постоянная времени T1 пропорциональна коэффициенту вязкого сопротивления b.
Практически, когда силы вязкого трения в механических элементах недостаточны, применяют дополнительное демпфирование подвижной системы, т. е. вводят дополнительную силу, противодействующую перемещению подвижной системы и пропорциональную скорости этого перемещения.
Если пневматический клапан применяется в системе с инерционным объектом, в котором переходные процессы протекают медленно, т. е. скорости изменения рвх и sвых небольшие, то величина ускорения d2Sвых/dt2 с точностью, достаточной для практических расчетов, может быть принята равной нулю. Тогда дифференциальное уравнение клапана примет вид:
.
Следовательно, в этом случае можно пренебречь инерционностью подвижных частей пневматического клапана и представлять его в динамическом отношении как апериодическое звено с передаточной функцией; определяемой формулой (10).
Па рис. 9,б приведена электрическая схема, переходный процесс которой также описывается дифференциальным уравнением второго порядка.
Постоянные времени и коэффициент передачи в этом случае равны:
;
;
.
При T1/T2 < 2 схема представляется колебательным звеном. Все три параметра схемы выражаются через одни и те же величины четырех сопротивлений и двух емкостей. Это является ее недостатком, так как параметры настройки, определяющие динамические свойства звена, взаимозависимы. Поэтому установка оптимальной величины одного из трех параметров настройки в большинстве случаев не дает возможности получить оптимальные значения также для двух остальных параметров. Кроме этого, такая настройка трудоемка и требует высокой квалификации наладчика.
