Васильев В.Г. Типовые звенья систем автоматического управления (колебательное звено). Методические указания
Подождите немного. Документ загружается.
Федеральное агентство по образованию
ТВЕРСКОЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ
Кафедра автоматизации технологических процессов
Типовые звенья систем автоматического управления
(колебательное звено)
Для студентов 3 курса специальностей 1905500
"Биотехнические и медицинские аппараты и системы"
и 1906600 "Инженерное дело в медико-биологической практике "
по курсу “Управление в медико-биологических системах ”
Tверь 2005
УДК [681.326+681322](075.8)
ББК 32.965 Я7
Методические указания предназначены для студентов 3 курса специальностей
1905500 "Биотехнические и медицинские аппараты и системы" и 1906600
"Инженерное дело в медико-биологической практике" по курсу “Управление в
медико-биологических системах ”. Содержит основные сведения о колебательном
звене теории систем автоматического управления, его переходных и частотных
характеристиках.
Обсуждено на заседании кафедры и рекомендовано к печати (протокол N 14
от 21 июня 2005 г.).
Составитель В.Г. Васильев.
©Тверской государственный
технический университет, 2005
3
Содержание
1. Временные характеристики колебательного звена…………………………3
2. Частотные характеристики колебательного звена………………………….8
3. Построение переходной характеристики колебательного звена
с помощью цифровых ЭВМ………………………….…………………………11
Использованная литература………………………………………..……...11
1. Временные характеристики колебательного звена
Звено называется колебательным, если при подаче на его вход единичного
ступенчатого воздействия процесс изменения его выходной величины будет
иметь форму затухающих амплитудных колебаний. Типичная для колебательного
звена переходная характеристика приведена на рис. 1.
Рис.1. Переходная характеристика колебательного звена
Связь между величиной на выходе и величиной на входе такого звена
определяется следующим дифференциальным уравнением второго порядка:
,
2
2
2
2
1
Kxy
dt
dy
T
dt
yd
T
=++
(1.1)
где К- коэффициент усиления, а Т
1
и Т
2
- постоянные времени звена. При этом
должно быть выполнено условие
12
2TT
<
.
Так что корни характеристического уравнения – комплексные.
Из (1.1) следует, что передаточная функция звена имеет вид
1
)(
2
22
1
++
=
pTpT
k
pW
.
(1.2)
Часто передаточную функцию записывают в виде
1
22
)(
0
2
0
22
00
2
2
0
++
=
++
=
p
p
k
pp
k
pW
ω
ξ
ω
ωξ ω
ω
,
(1.3)
где ω
0
- собственная частота незатухающих колебаний звена, а ξ - степень
затухания (успокоения).
Возможна и такая запись передаточной функции:
4
12
)(
22
++
=
TppT
k
pW
ξ
.
(1.4)
Сравнивая (1.2) , (1.3) и (1.4) можно записать:
• постоянная времени звена
TT
==
0
1
1
ω
;
• постоянная времени звена
0
2
2
ω
ξ
=
T
;
• степень затухания
1
2
2T
T
=
ξ
.
Запишем передаточную функцию (1.1) в виде
)1)(1(
1
)(
43
43
2
22
1
++
=
++
=
pTpT
kk
pTpT
k
pW
,
(1.5)
где
43
2
1
TTT
=
и
432
TTT
+=
, и найдем корни полиномов знаменателей этого
выражения:
3
2
1
2
1
2
22
1
1
2
4
T
T
TTT
p
−=
−+−
=
;
(1.6)
4
2
1
2
1
2
22
2
1
2
4
T
T
TTT
p
−=
−−−
=
,
(1.7)
откуда находим
2
1
2
22
2
1
3
4
2
TTT
T
T
−−
=
;
2
1
2
22
2
1
4
4
2
TTT
T
T
−+
=
.
(1.8)
Далее запишем (1.5) как
+
+
=
++
=
43
43
43
43
43
11
)1)(1(
)(
T
p
T
p
TT
kk
pTpT
kk
pW
.
Обозначим
43
43
TT
kk
Z
=
3
1
T
=
α
4
1
T
=
β
.
Тогда
))((
)(
βα
++
=
pp
Z
pW
.
Данному изображению соответствует оригинал
αββα
βα
−
−
↔
++
−−
tt
ee
pp ))((
1
.
5
Найдем весовую функцию звена
−
−
=
−
−
=
−
−
=≡
−
−
−
−
−−
43
4
3
43
43
34
43
43
11
)()(
t
t
t
t
t
t
t
t
tt
ee
TT
kk
TT
ee
TT
kk
ee
Ztytg
αβ
βα
.
Обозначим
43
43
TT
kk
L
−
=
и найдем переходную функцию звена:
=−++−=+−=−=
−−−−−−
∫∫
44330403
00
4
3
4
3
4
3
)( LTeLTLTeLTeLTeLTdeLdeLth
T
t
T
t
t
T
t
T
t
T
t
T
τ
τ
τ
τ
ττ
−
+
−
−=
+−−
−
=
−−−−
4
3
4
3
43
3
43
3
434343
43
43
1)(
T
t
T
t
T
t
T
t
e
TT
T
e
TT
T
kkeTeTTT
TT
kk
.
(1.9)
Далее, используя соотношения (1.8), найдем
2
1
4
2
1
2
1
2
2
2
43
3
+
−
=
−
TT
T
TT
T
;
2
1
4
2
1
2
1
2
2
2
43
4
−
−
=
−
TT
T
TT
T
.
(1.10)
Обозначим
2
2
2
1
4 TTd
−=
и, учитывая (1.10), можно записать
2
1
2
2
43
3
+=
−
jd
T
TT
T
;
2
1
2
2
43
4
−=
−
jd
T
TT
T
.
(1.11)
Далее с учетом (1.8) находим
2
1
2
1
2
2
1
2
3
22
2
T
dt
j
T
tT
T
jdT
t
T
t
eeee
−
−
−
−
==
;
2
1
2
1
2
2
1
2
4
22
2
T
dt
j
T
tT
T
jdT
t
T
t
eeee
−−
+
−
−
==
.
(1.12)
Подставляя (1.11) и (1.12) в (1.10), получим
.
2
1
2
2
2
1
2
cos
2
1
2
sin
2
1
43
2
1
2
2
2
1
2
2
1
2
1
2
2
1
2
1
2
2
2
1
2
2
2
2
1
2
2
2
1
43
2
1
2
2
1
2
2
2
1
2
2
2
1
2
2
1
2
2
2
1
2
2
1
43
)(
−
+−=
=
−−
+−
−−
−−−=
=
−−
−+
−
+−=
T
tT
e
T
dt
T
dt
d
T
kk
T
tT
e
T
dt
j
e
T
dt
j
e
T
tT
e
T
dt
j
e
jd
T
T
dt
j
e
jd
T
kk
T
dt
j
e
T
tT
e
jd
T
T
dt
j
e
T
tT
e
jd
T
kkty
6
Или в окончательном виде
+−=
−
2
1
2
1
2
2
2
cos
2
sin1)(
2
1
2
T
dt
T
dt
d
T
ekty
T
tT
,
(1.13)
где к = к
1
к
2.
Обозначим:
2
1
2
2
T
T
T
=
- постоянная времени затухания амплитуды колебаний (время затухания
колебаний будет t
c
≈ 3T).
В выражении (1.13) под знаком sin (.) и cos (.) множитель
2
1
2T
d
имеет смысл
частоты ω
0
незатухающих колебаний звена
.
2
4
2
2
1
2
2
2
1
2
1
0
T
TT
T
d
−
==
ω
Умножим левую и правую части этого выражения на постоянную времени
затухания колебаний
2
1
2
2
T
T
T
=
. В результате получим
2
2
2
2
1
0
4
T
TT
T
−
=
ω
.
Обратная величина
d
T
TT
T
T
2
2
2
2
1
2
0
4
1
=
−
=
ω
, что соответствует множителю в
выражении (1.13) перед знаком sin(.).
Тогда выражение (1.13) принимает вид
+−=
−
t
T
teKty
T
t
0
0
0
sin
1
cos1)(
ω
ω
ω
.
Это выражение можно записать и в другом виде. Для этого сделаем следующие
преобразования
1
2
1
2
1
2
2
1
1
2
2
2
1
1
2
1
2
2
2
1
2
1
0
1
2
1
2
2
1
4
2
1
2
4
2
TT
T
T
T
T
TT
T
T
TT
T
d
ξ
ω
−
=
−
=
−
=
−
==
;
(1.14)
2
2
2
2
1
1
1
2
2
2
2
1
22
1
4
2
1
2
4
ξ
ξ
−
=
−
=
−
=
TT
T
T
T
TT
T
d
T
,
где
1
2
2T
T
=
ξ
- параметр затухания. С учетом полученных соотношений будем иметь
7
−
−
+
−
−=
−
t
T
t
T
eKty
T
t
1
2
2
1
2
1
sin
1
1
cos1)(
1
ξ
ξ
ξ
ξ
ξ
.
Анализ показывает, что при ξ=0 (T
2
= 0) колебания являются незатухающими
[ ]
)cos(1)(
0
tkty
ω
−=
и частота, с которой будут происходить колебания,
1
0
1
T
=
ω
. Отсюда становится
ясным справедливость записи
2
00
2
2
0
2
1
2
4
2
2
2
1
2
22
1
2
1
1
)(
ωξ ω
ω
++
=
++
=
++
=
pp
k
T
p
T
T
p
T
k
pTpT
k
pW
и понятна роль постоянной времени Т
2
. Постоянная Т
2
характеризует
демпфирование (ослабление) колебаний (чем она больше, тем больше степень
затухания амплитуд переходного процесса). Постоянная Т
1
”раскачивает”
колебания (чем больше Т
1,
тем меньше степень затухания
1
2
2T
T
=
ξ
).
Говорят, что при
ξ
=0 система не демпфирована. При
ξ
>0 – система не
дедемпфирована; при
ξ
=1 – система обладает критическим демпфированием, при
ξ
>1 – система передемпфирована.
При ξ = 1 2Т
1
= T
2
и колебания вырождаются в апериодический процесс
(ω
0
= 0).
При Т
2
= 0 характеристическое уравнение звена принимает вид
01
22
1
=+
pT
и будет иметь чисто мнимые корни
0
1
2,1
1
ω
j
T
jp
±=±=
,
что и объясняет режим незатухающих колебаний в звене.
При Т
1
= 0 характеристическое уравнение звена принимает вид как и для
апериодического звена с одним действительным корнем. Переходный процесс
будет иметь апериодический характер.
Характеристическое уравнение
01
2
22
1
=++
pTpT
имеет комплексно-сопряженные корни
02,1
1
ω
j
T
p
±−=
с отрицательной
вещественной частью, чем и объясняется режим затухающих колебаний.
8
2. Частотные характеристики колебательного звена
Из (1.2) следует, что амплитудная фазовая частотная характеристика звена
( )
[ ]
.
)()()(
)(
2
2
22
2
2
1
2
2
2
22
2
2
1
2
2
1
2
2
22
2
2
1
2
2
2
1
2
2
2
12
2
2
1
11
1
1
1
11
)()(
ωω
ω
ωω
ω
ωω
ωω
ωω
ωω
ω
ωω
T
J
TT
T
kT
T
T
k
T
jkT
T
k
TT
k
j
T
j
T
k
j
W
+
−
+
=
+
−
=
=
++
=
++
=
−−
−
−
−
−
(2.1)
Годограф амлитудно-фазочастотной характеристики приведен на рис. 2.
Рис. 2. Годограф амлитудно-фазочастотной характеристики
колебательного звена
Находим амплитудно-частотную характеристику
( )
( )
ωω
ω
ω
2
2
2
2
2
2
1
1
T
T
k
j
W
A
+
==
−
(2.2)
и фазочастотную характеристику звена
( )
( )
−
=
−
ω
ω
ω
ϕ
2
2
1
2
2
2
1
T
T
arctg
.
(2.3)
9
Вид характеристик приведен на рис. 3.
Рис. 3.a. Амплитудно-частотная
характеристика
Рис. 3.б. Фазо-частотная характеристика
Существование максимума на амлитудно-частотной характеристике означает, что
система обладает резонансными свойствами, а в переходном процессе будет
наблюдаться перерегулирование. Для систем более высокого порядка чем больше
максимум АЧХ, тем сильнее выражен резонанс и тем больше перерегулирование
в переходном процессе. Поэтому желательно, чтобы САР не обладала ярко
выраженными резонансными свойствами.
Так как ω
0
=
T
1
1
, то амплитудно-частотную характеристику можно записать в
виде
( )
.
22
0
0
2
1
+
−
=
ω
ω
ξ
ω
ω
ω
k
A
(2.4)
Найдем при какой частоте ω характеристика A(ω) имеет максимум. Обозначим
эту резонансную частоту -
ω
2
.
Положим x =
ω
ω
0
. Тогда выражение (2.4) можно представить в виде
)(xkyx
A
=
,
где
( )
( )
.
2
2
2
2
1
1
)(
x
x
xy
ξ
+
=
−
Найдем производную y(x), а затем для того чтобы найти экстремум функции,
приравняем ее нулю.
10
( ) ( )( )
=
+−⋅
+−−=
−
−
ξξξ
22
2
2
12
2
1
2
1
)(
2
2
3
22
2
x
xx
x
x
x
y
I
( )
( )
( )
.
322
2
2
1
1
2
2
2
+
−
−−
=
x
x
xx
ξ
ξ
Приравняем числитель этого выражению к нулю:
2x (1 - 2ξ
2
- x
2
) = 0,
откуда находим: первый максимум при x = 0. Второй максимум будет иметь
место при
x =
ξ
2
21
−
,
или с учетом обозначения
ξ
ωω
2
0
21
−=
r
.
Подставим это выражение в (2.4). Тогда
( )
=
+
=
=
+
−
=
−
−
−
2
0
2
0
2
2
2
0
2
2
0
21
2
21
1
2
1
0
2
0
ω
ξ
ω
ξ
ω
ξ
ω
ω
ω
ξ
ω
ω
ω
k
k
A
rr
ξξξξξ
2424
1
2844
−
=
−+
=
kk
.
Из последнего выражения следует, что при
ξ
, стремящемся к нулю, амплитуда
колебаний будет стремиться к бесконечности, что с физической точки зрения
означает явление резонанса. Вид кривых АЧХ в зависимости ζ приведен на рис. 4.