Электронная автоматика и техника. Ч. 2. Работа в среде Electronics Workbench
.pdf
Получаем следующую расчетную структурную схему:
W245(p)
2
Рассчитываем общую передаточную функцию W0(p) для соединения с отрицательной обратной связью по выражению (7.4):
W0 |
( p) |
|
W245 ( p) |
|
|
p T25 |
. |
|
W245 ( p) W2 |
|
1 ( p T24 1) ( p T5 1) k2 |
||||
|
1 |
( p) |
|
||||
Для получения окончательного вида находим корни квадратного уравнения в знаменателе:
k2 p2 T24 T5 k2 p (T24 T5 ) (k2 1) 0;
2 p2 1 0,01 2 p (1 0,01) (2 1) 0;
0,02 p2 2,002 p 3 0; |
p |
|
|
2,002 |
(2,002)2 4 0,02 3 |
; |
|
|
|
|
|
|
|||||
|
|
1,2 |
|
|
2 0,02 |
|
|
|
|
|
|
|
|
|
|
|
|
p 2,002 |
4,008004 0,24 |
; |
|
p 2,002 1,9411347196 |
; |
|||
1,2 |
0,04 |
|
|
|
1,2 |
0,04 |
|
|
|
|
|
|
|
|
|
||
|
p1 0,152162; |
|
|
p2 98,5784. |
|
|
||
Окончательный вид общей передаточной функции будет таким:
WO ( p) |
|
p T25 |
|
|
|
p T25 |
|
|
|
|
|
2 ( p 0,152162) ( p 98,5784) |
|||||
|
1 ( p T24 1) ( p T5 1) k2 |
|
||||||
|
|
|
p T25 |
|
|
|
|
|
|
2 0,152162 98,5784 ( p / 0,152162 1) ( p / 98,5784 1) |
|
||||||
|
|
|
p T1 |
|
. |
|
||
|
|
( p T 1) ( p T 1) |
|
|||||
|
|
|
2 |
|
3 |
|
|
|
где T1 T25 /(2 0,152162 98,5784) 0,0066667;
70
T2 1/ 0,152162 6,571943;
T3 1/ 98,5784 0,00101442.
Полученное выражение для передаточной функции W0(p) характеризует собой инерционное дифференцирующее звено.
Рассчитываем комплексную передаточную функцию W (i ) и находим действительную P(w) и мнимую Q(w) части выражения
К комплексной передаточной функции переходим формальной заменой p i w :
W (i ) (T2 i 1)(T3 i 1) .
Для устранения комплексного значения i в знаменателе умножаем числитель и знаменатель на сопряженные комплексные выражения и делим действительную и мнимые части полученного выражения:
W (i ) |
|
T1 i |
|
|
T1 i (T2 i 1)(T3 i 1) |
|
|||||||
(T |
i 1)(T |
i 1) |
2 |
2 |
|
2 |
2 |
|
1) |
||||
|
|
|
|||||||||||
2 |
3 |
|
|
( T |
|
|
1)( T |
|
|
|
|||
|
|
|
|
|
2 |
|
|
3 |
|
|
|
|
|
(T1 i2 2 T1 i )(T3 i 1) (T22 2 1)(T32 2 1)
T1 i3 3T2T3 T1 i2 2T2 T1 i2 2T3 T1 i
(T22 2 1)(T32 2 1)
T1 i 3T2T3 T1 2T2 T1 2T3 T1 i P( ) iQ( ). (T22 2 1)(T32 2 1)
После преобразований получим
P( ) |
|
2T (T |
T ) |
|
; |
||
1 |
3 |
|
2 |
||||
|
(T 2 |
2 1)(T 2 |
2 |
1) |
|||
2 |
|
|
3 |
|
|
|
|
Q( ) |
T (1 2T T ) |
|
|
||||
|
1 |
|
2 |
3 |
|
. |
|
(T 2 |
|
1)(T 2 |
|
|
|
||
|
2 |
2 |
1) |
||||
|
2 |
|
3 |
|
|
|
|
71
Используя полученные выражения для P( ), Q( ) , рассчитыва-
ем остальные частотные характеристики.
АФХ (годограф) строится на комплексной плоскости в координатах P( ) и Q( ) для различных значений циклической частоты [0, ) .
Логарифмическая частотная характеристика (ЛЧХ) рассчитывается по выражению и строится в логарифмических координатах по оси частот по декадам, т.е. при десятикратном изменении частоты
La ( ) 20 lg A( ) , дБ.
Логарифмическая частотная характеристика (ЛАЧХ) La ( ) стро-
ится на одном графике (совместно) с логарифмической частотной характеристикой L( ) . Для ее построения понадобятся частоты со-
пряжения ci , которые определяются из скобок ( pT 1) выражения для общей передаточной функции W(p):
|
|
1 |
0,152162 и |
|
|
1 |
98,5784. |
c1 |
|
T2 |
c2 |
|
T3 |
||
|
|
|
|
||||
Переходная характеристика h(t) реакция системы на единичное ступенчатое воздействие и ее изображение по Лапласу H ( p) рассчитываются следующим образом:
H ( p) W ( p) |
1 |
|
|
T1 p |
|
|
|
|
|
|
|
T1 |
|
|
|
. |
p |
p ( pT |
1)( pT |
1) |
T T |
|
( p |
|
1 |
)( p |
1 |
) |
|||||
|
|
2 |
3 |
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
2 |
3 |
|
|
|
T2 |
|
T3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Корни характеристического выражения знаменателя:
p |
1 |
; |
p |
|
|
1 |
. |
1 |
T2 |
|
|
2 |
|
T3 |
|
|
|
|
|
|
|||
Используем для перехода от изображения к оригиналу для расчета h(t) обратное преобразование Лапласа – разложение Хевисайда:
|
|
|
|
|
t |
|
|
|
|
|
t |
|
h (t) |
|
T1 |
|
e |
T2 |
|
T1 |
|
e |
T3 |
||
|
T |
|
|
|
|
|||||||
|
T |
|
|
|
|
T T |
|
|
||||
2 |
3 |
|
|
|
3 |
2 |
|
|
|
|||
0,001014468(e |
t / 6,571943 |
e |
t / 0,00101442 |
|||||||||
|
|
|
|
|
|
|
). |
|||||
72
Результаты расчетов в Excel
73
74
75
76
Из графика h(t) определяем на уровне 5% от максимального значения h = 0,00075 время регулирования Tр, которое примерно равно 17 с.
Построение логарифмической асимптотической частотной характеристики (ЛАЧХ) делается по виду передаточной функции и изменению текущего угла наклона ЛАЧХ в точках сопряжения. Передаточную функцию нужно перед этим привести к одному из трех возможных видов:
W ( p) K |
( p T1 |
1) ( p T2 |
1) ... ( p Tl |
1) |
|
|
|||||||||
|
|
|
|
|
|
|
; |
|
(7.11) |
||||||
( p T |
1) ( p T |
1) ... ( p T |
1) |
|
|||||||||||
|
|
|
|
|
m |
|
n |
|
z |
|
|
|
|
|
|
W ( p) p T |
( p T1 |
1) ( p T2 |
1) ... ( p Tl |
1) |
; |
(7.12) |
|||||||||
( p T |
1) ( p T |
1) ... ( p T |
1) |
||||||||||||
|
|
|
|
|
|
m |
|
n |
|
z |
|
|
|
|
|
1 |
|
|
( p p1) ( p |
p2 ) ... ( p pl ) |
|
|
|||||||||
W ( p) |
|
|
|
|
. |
|
(7.13) |
||||||||
p T |
( p pm ) p |
pn ) ... ( p pz ) |
|
||||||||||||
Количество скобок вида ( p pi ) может быть любым (от 0).
При этом начальный угловой коэффициент наклона определяется так:
– для передаточной функции W(p) вида (7.11) угловой коэффициент наклона равен 0 дБ/дек (децибел на декаду), т.е. ЛАЧХ проходит горизонтально при очень низких частотах ( 0). Началь-
ное положение ЛАЧХ равно 20 lg(K) , K – из выражения (7.11);
– для вида (7.12) угловой коэффициент наклона равен 20 дБ/дек, т.е. ЛАЧХ вначале при очень низких частотах ( 0) возрастает.
Начальное положение ЛАЧХ неизвестно, т.е. будет построен общий характер изменения ЛАЧХ. В реальном случае кривая будет смещена по вертикали вниз или вверх. Для определения точного начального положения нужен или расчет, или ЛАЧХ проводят, приближая
краcсчитанной ЛАХ;
–для вида (7.13) угловой коэффициент наклона равен – 20 дБ/дек, т.е. ЛАЧХ вначале при очень низких частотах ( 0) уменьшает-
ся. Начальное положение ЛАЧХ как и в случае (7.12) неизвестно,
77
т.е. будет построен лишь общий характер изменения ЛАЧХ. В реальном случае кривая будет смещена по вертикали вниз или вверх. Для определения точного начального положения нужен или расчет, или ЛАЧХ проводят, приближая к рассчитанной ЛАХ.
Во всех случаях дальнейшее изменение угла наклона происходит в точках сопряжениях. Они определяются для каждой из скобок ви-
да ( p Ti 1) :
|
|
1 |
. |
(7.14) |
|
||||
C i |
|
Ti |
|
|
|
|
|
||
В каждой из таких частот сопряжения C i |
происходит измене- |
|||
ние текущего угла наклона на 20 дБ/дек, причем если скобка из которой определена частота сопряжения C i находится в числителе
передаточной функции W(p) (7.11)–(7.13), то изменение текущего угла наклона ЛАЧХ происходит на +20 дб/дек, а если в знаменателе – на –20 дБ/дек.
В нашем примере находим все частоты сопряжения и расставляем их на оси частот (рис. 7.5):
|
|
|
1 |
|
|
1 |
|
0,152162 – из знаменателя, |
||||
|
|
|
|
|
||||||||
C 2 |
|
|
T2 |
|
|
6,571943 |
|
|
||||
|
|
|
|
|
|
|
|
|||||
|
|
|
1 |
|
|
|
1 |
|
|
98,5784 – из знаменателя. |
||
|
|
|
|
|
|
|
||||||
C 3 |
|
|
T3 |
0,001014420 |
||||||||
|
|
|
||||||||||
Так как передаточная функция относится к виду (7.12), то ее начальный угловой коэффициент наклона равен +20 дБ/дек. При частоте C 2 он изменится на –20 дБ/дек, так как найден для скобки
из знаменателя (показано стрелкой вниз справа от C 2 ), и станет равным +20 – 20 = 0 (дБ/дек). В точке C 3 он также изменится на
–20 дБ/дек и станет равным 020 = –20 (дБ/дек). Для построение углов наклона строим вспомогательный прямоугольник со сторонами 1 декада и 20 дБ и проводим диагонали с угловыми коэффициента-
78
ми наклона +20 и –20 дБ/дек. Параллельно этим диагоналям и проводим прямые на ЛАЧХ. Вначале с угловым коэффициентом наклона +20 дБ/дек до пересечения со значением частоты C 2 , да-
лее горизонтально с угловым коэффициентом 0 дБ/дек до пересечения с частотой C 3 , после чего с наклоном –20 дБ/дек. Строим
ЛАЧХ на графике совместно с ЛАХ (рис. 7.5).
При построении приближаем, перемещая по вертикали, ЛАЧХ к ЛАХ (на полученном рисунке кривая Lа (ЛАЧХ) незначительно смещена вверх от кривой L (ЛАХ) для удобства иллюстрации построения).
Проверяем аналитические расчеты, моделируя САР в Workbench. Для этого составляем общую электрическую принципиальную схему САР и набираем ее в Workbench. Подключаем body plotter для частотного анализа и генератор прямоугольных колебания для получения реакции на прямоугольное ступенчатое воздействие, т.е. для получения формы переходной характеристики. Для примера участок схемы одного из вариантов выглядит так (рис. 7.6).
Рис. 7.6. Исследование САР в Workbench
79
