Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Клюев, А. С. Автоматическое регулирование

.pdf
Скачиваний:
137
Добавлен:
22.10.2023
Размер:
17.2 Mб
Скачать

5-5. ИНТЕГРАЛЬНЫЕ ОЦЕНКИ КАЧЕСТВА ПРОЦЕССА

Динамической ошибкой переходного процесса, возни­ кающей при изменении воздействия на систему (см. рис. 5-1), называется переменная во времени раз­ ность

x (t)= h(oо)—h(t),

(5-46)

где h(оо) — установившееся значение выходной

величи­

ны после окончания переходного процесса; h(t)

— теку­

щее переменное значение выходной величины в течение переходного процесса.

Рис. 5-14. Графическое представление линейной

интегральной

оценки

качества

процесса регулирования.

 

 

а —при

отсутствии

перерегулирования;

б — при

наличии

перерегулиро­

вания.

 

 

 

 

 

На рис. 5-14 показана кривая x(t) переходного про­

цесса без перерегулирования, а на

рис. 5-14,6 — такая

же кривая, но с перерегулированием

(см. рис. 5-1). Из

этих кривых видно, что чем

меньше заштрихованная

площадь на рис. 5-14,а, тем быстрее ликвидируется дина­ мическая ошибка.

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

Заштрихованная площадь на рис. 5-14,а измеряется определенным интегралом функции x(t) при изменении t от 0 до оо:

(5-47)

260

Интеграл Л называется линейной интегральной оцен­ кой качества.

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

Преобразуем по Лапласу выражение (5-46), учитывая, что согласно формуле (3-22) для замкнутой системы

L[h(t)\ = - ^ - ф(р) и что изображение постоянной вели­ чины h(оо)(при g0 = l) равно:

L[h(oo)\ — Ь[Ф (0) 1 ] = - ^ ^ . '

С учетом этого в преобразованном по Лапласу виде выражение (5-46) запишется как

Х(р) = фМ ^ і фМ .

Вместе с тем, полагая р— >-0 в выражении (2-7) и учитывая равенство (5-47), имеем:

\

Х (0 )= \x{t)dt = I

о

Таким образом, линейная интегральная оценка каче­ ства Л при поступлении на вход системы единичного ступенчатого задающего воздействия определяется через значения передаточной функции замкнутой системы для установившегося и неустановившегося ее состояний:

/, = lira фМ — ф (РУ.

(5.48)

р->0 Р

Так, например, для системы, являющейся в динами­ ческом отношении апериодическим звеном с передаточ­ ной функцией, определяемой уравнением (2-22), линей­ ная интегральная оценка качества равна:

/, = lim ■

Т р + 1

=lira

k T р

= 1ІШ:

kT

~=kT.

 

{ Т р + \ ) р

 

р->0

 

 

^ ( 7 - / » + ! )

 

261

/, = lim

Линейная интегральная оценка качества реального дифференцирующего звена с передаточной функцией, определяемой выражением (2-53), равна:

О

ІіТр

Т р + 1

р->0

кТ.

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

Линейная интегральная оценка качества звеньев обо­ их рассмотренных видов тем больше, чем больше по­ стоянная времени звена п его коэффициент передачи. Если система имеет передаточную функцию иитегро-

дпффереицирующего звена [см. формулу

(2-57)], то

к

Н ТлР+ 1)

 

 

 

ТдР 4- 1

 

к(Тп~ Т л)

 

/, = lim —

= lim

= А(ГЯ- Г Д).

Р - * о

Р

р->о

ТъР+ 1

 

Линейная интегральная оценка качества этой системы зависит от соотношения постоянных времени Г„ н Тя.

Если в переходном процессе имеются перерегулирова­ ния или этот процесс колебательный, то динамическая ошибка регулирования (см. рис. 5-14,6) за время пере­

ходного

процесса неоднократно меняет знак. Если и

в этих

случаях пользоваться линейной интегральной

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

Например, при установившихся незатухающих коле­ баниях выходной величины положительные и отрицатель­ ные площади, ограниченные кривой переходного процес­ са, будут равны и линейная интегральная оценка каче­ ства будет равна нулю. Если исходить из линейной инте­ гральной оценки качества, следует признать качество та­ кой системы очень хорошим; фактически Же эта система по качественным показателям практически непригодна для использования. В связи с этим при наличии перере­ гулирования в переходном процессе и колебательном

262

его характере применяют

квадратичную

интегральную

оценку качества вида

 

 

со

X2 (t) dt.

 

/, = j

(5-49)

о

 

 

Квадратичная интегральная оценка качества равна площади, ограниченной кривой xz(t) (рис. 5-15,а).

Однако квадратичная интегральная оценка качества процесса регулирования в ряде случаев также не дает объективной оценки характера переходного процесса.

Рис. 5-15. Графическое представление квадратичной интеграль­ ном оценки качества процесса регулирования.

В качестве примера рассмотрим (рис. 5-15,6) кривые апериодического процесса (кривая 1) и колебательного процесса (кривая 2) с большим коэффициентом колеба­ тельности.

Как следует из рис. 5-15,6, интегральная квадратич­ ная оценка качества колебательного переходного процес­ са будет меньше, чем апериодического процесса. Факти­ чески же в данном случае более приемлемым является апериодический процесс.

Недостатки рассмотренных выше интегральных оце­ нок качества обусловили появление так называемого обобщенного интегрального критерия вида

х

СО

 

 

\V (t) dt,

(5-50)

 

6

 

2 6 3

где V(t) — квадратичная форма от отклонения регулш руемой величины и ее -производных, например

(ІХ

d nx \

d t + ь

+ - • - -Ьтг d t n )•

Физический смысл критерия Іѵ состоит в том, что минимизируется не только величина и длительность от­ клонения регулируемой величины от заданного значения,

но и производных от

отклонения

(скорость изменения

отклонения, ускорение

и т. д.).

• • •, уп выбираются

Постоянные коэффициенты уі,

с учетом ограничений отдельных динамических показа­ телей переходного процесса.

Так,

чем больше

коэффициент

у,, тем

больше

роль

 

 

00

 

 

 

второго

слагаемого

^ l (dx/dt)a dt,

тем при равных

зна-

 

 

о

 

плавным, но

чениях Іѵ переходный процесс будет более

зато и более длительным.

Для большинства практических систем к переходно­ му процессу предъявляются требования минимизации квадратичного интеграла динамической ошибки (заштри­ хованных площадей на рис. 5-15) и скорости изменения регулируемой величины.

В этом случае обобщенный интегральный критерий качества запишется в виде

•=П-+Ч£)Чdt.

Это выражение

можно преобразовать:

A , =

J[ *

+

Ѵ ъ

( і г

) *dt — 2 y Y ^ x d x =

 

 

 

 

 

 

f

 

 

 

 

 

' d t+

Y > 2

Так как

x(oo) = 0,

то,

обозначив

х(0) = х 0 и Y b —

= Т, получим:

 

 

00

 

 

 

 

 

 

 

 

 

Іѵ=

Т^0-\- j" J^c+ Г

dt.

264

При заданной величине Т этот интеграл имеет мини­ мум, если подынтегральная функция

Т Т Т + * = ° -

Решая это дифференциальное уравнение, находим вы­ ражение оптимального переходного процесса для дан­

ных систем:

_t_ X = х йе т.

Таким образом, оптимальным переходным процессом в этом случае будет изменение регулируемой величины по экспоненте с постоянной времени Т, определяемой за­ данной максимально допустимой скоростью изменения регулируемой величины.

Методы нахождения интегральных оценок качества переходного процесса изложены, например, в [Л. 5 и 30].

5-6. ПРИБЛИЖЕННОЕ ПОСТРОЕНИЕ ПЕРЕХОДНОГО ПРОЦЕССА ПО ВЕЩЕСТВЕННОЙ ЧАСТОТНОЙ ХАРАКТЕРИСТИКЕ СИСТЕМЫ

Показатели качества системы можно определить, имея график переходного процесса (см. рис. 5-1). Но так как для получения графика переходного процесса необ­ ходимо аналитическое решение дифференциального урав­ нения системы, что является во многих случаях сложной задачей, то при расчетах автоматических систем регули­ рования нашли широкое применение приближенные ме­ тоды анализа переходных процессов.

Согласно формулам (3-'22) и (3-76) переходный про­ цесс замкнутой системы полностью определяется ее пере­

даточной функцией или АФХ,

которую' в соответствии

с (3-38)

можно разделить на

вещественную Uф(ю) и

мнимую

Ѵф(со) части:

 

Ф (/(0 ) — П ф ( а ) + ] У ф ( © ) .

Каждой АФХ однозначно соответствует веществен­ ная частотная характеристика {/$(■©), причем обе эти характеристики, как это видно из выражений (3-28) и (3-36), зависят только от всех коэффициентов дифферен­ циального уравнения; поэтому вещественная частотная характеристика Пф(со) замкнутой - системы однозначно определяет характер ее переходного процесса.

?§5

К переходным процессам в линейных АСР применим закон суперпозиции, как указывалось в гл. 1; это значит, что если входное воздействие представить как сумму со­ ставляющих воздействии и найти уравнения или постро­ ить кривые переходных процессов в системе для каждой составляющей порознь, то переходный-процесс, создавае­

 

 

 

 

мый

входным

 

воздейст­

 

 

 

 

вием

в целом,

будет

ра­

 

 

1\

 

вен

сумме

 

переходных

%

^

 

процессов

для

всех

со­

¥

а.

г г

 

ставляющих

воздействий.

 

11

 

 

 

11

 

На

этом

 

II

 

основана

¥

 

11

 

методика определения пе­

 

 

I 1

 

¥

 

11

 

реходного

процесса

по

¥

 

1{

 

вещественной

частотной

О

 

* , \ ?

 

характеристике

замкну­

 

 

ш , ш 2

 

 

 

 

 

той системы

при поступ­

 

 

 

 

лении на ее вход единич­

 

 

 

 

ного

ступенчатого управ­

 

 

 

 

ляющего воздействия.

так­

 

 

 

 

Можно

было

бы

 

 

 

 

же построить

кривую

пе­

 

 

 

 

реходного

процесса и по

 

 

 

 

мнимой

частотной харак­

 

 

 

 

теристике

замкнутой

си­

 

 

 

 

стемы.

 

замкнутая

си­

 

 

 

 

Пусть

 

 

 

 

стема

 

имеет

 

веществен­

 

 

 

 

ную

частотную

характе­

 

Uffj-Upt

 

ристику,

изображенную

 

 

 

 

па рис. 5-16,а.

 

 

 

Рис. 5-16. Представление веще­

Предположим,

что

ственной частотной

характеристи­

входное воздействие,

при­

ки

замкнутой системы суммой

ложенное

к

системе,

со­

трапецеидальных

вещественных

стоит

из

отдельных

со­

частотных характеристик.

ставляющих,

 

веществен­

ные

частотные

 

 

характеристики для

каждой

из кото­

рых имеют вид трапеции; при этом две стороны трапеции совпадают е осями кобрдинат со и Vф (со), третья парал­ лельна оси со, а четвертая наклонна.

Выбор трапеции в качестве типовой формы состав­ ляющих вещественных частотных характеристик обуслов­ лен следующими соображениями:

1) действительные вещественные частотные характе-

266

ристикп реальных систем регулирования легко расчленя­ ются на небольшое число трапецеидальных составляю­ щих;

2)вычисление ординат кривой переходного процесса для трапецеидальной вещественной частотной характе­ ристики достаточно просто и может быть облегчено ис­ пользованием таблиц;

3)точность расчетов достаточно велика.

Так как сумма переходных процессов от отдельных

составляющих образует

переходный процесс системы

в целом, то и сумма состав­

 

ляющих

вещественных

ча­

 

стотных характеристик, име­

 

ющих

форму

трапеций,

 

образует

вещественную

ча­

 

стотную

характеристику

си­

 

стемы (рис. 5-16,6).

приняв

 

Таким

образом,

 

трапецеидальную

форму

ве­

 

щественной

частотной

ха­

 

рактеристики за

типовую и

Рис. 5-17. Трапецеидальная

составив

таблицы

ординат

вещественная частотная харак-

кривых переходного

процес-

теристика.

са /г(т) для единичных трапеций с различными наклона­ ми четвертой стороны, мы можем с помощью таблиц и не­ сложных пересчетов построить переходные процессы для каждой составляющей вещественной частотной характе­ ристики и, просуммировав их ординаты, получить кри­ вую переходного процесса замкнутой системы.

Любая трапецеидальная вещественная частотная ха­ рактеристика (рис. 5-17) характеризуется высотой ІУфо(сй), интервалом пропускания частот юо, интервалом равномерного пропускания частот соц и коэффициентом

наклона х-и^/іоо-

единичной трапеции принимается

Для типовой

£Лі)о(сй) = 1 и шо= 1;

поэтому единичная трапеция харак­

теризуется только коэффициентом наклона x=cöd.

Для единичных трапеций с различными величинами X могут быть вычислены ординаты переходного процесса в виде Іг(т), где т = ‘соо^ — безразмерный параметр време­ ни. Из сказанного ранее вытекает, что для единичной трапеции т = 1• t — t.

Величины ординат переходных процессов, вычислен­ ные для различных х и х, называются Л-функциями.

267

'Значения /i-функциіі приведены в приложении 2.

Для перехода от /г-фуикции к переходной функции x{t), соответствующей данной составляющей трапецеи­ дальной вещественной частотной характеристике с тем лее %, но с Д ф 0 ( с о ) ^ 1 и соо+ М , необходимо значения /г-функции умноліить на Дф0(ю), а для перехода к но­

вому значению времени необходимо учесть,

что т==то0^

и, следовательно, t = т/со0. Поэтому

x(t) = Пф0(о))/г(т/со0).

На рнс. 5-1б,а выполнена разбивка вещественной частотной ха­

рактеристики системы на трапеции, которые отдельно

представлены

на рис.

5-16,6.

Здесь

они

расположены

 

так, что

их

основания

к— г, 0—е и л—б совпадаю т с осью со.

U,|, (со)

 

 

 

В

данном

случае

характеристика

может

быть за ­

менена

тремя

составляющими трапецеидальными

характеристиками

л—в— г— к, л— б — а и кд— е—0, которые

характеризуются следую ­

щими параметрами:

 

 

 

 

 

 

 

 

 

С/фоі= С/фі—С/фг= 1,2+0,3= 1,5;

ffidi= W2=2;

 

 

 

ш01 =

ш3 =

4; к =

— — =

0,5;

 

 

 

 

 

 

 

 

W01

 

 

 

 

 

 

 

U Ф02~ UфоU ф 1= 1,0— 1,2= — 0,2;

 

 

 

 

 

Wd2=0;

(Оо2=Ы| = 1,5;

и=0;

 

 

 

 

 

’С фоз =!£Уф2= —0,3;

Wd3=с0-і=5;

 

 

 

 

 

w03 =

“ б =

1Л

——

Гі с

 

 

 

 

 

10; X =

= 0 ,5 .

 

 

 

 

 

 

 

 

ш03

 

 

 

 

П о таблицам й-функций (см. Приложение 2) находим значения ординат переходных функций /і(т) для единичных трапецеидальных вещественных частотных характеристик при х = 0 ,5 и х = 0 .

По /і-функцин для х = 0 ,5 находим переходные функции, соответ­ ствующие первой трапеции Хі(/) и третьей трапеции x3(t).

Рис. 5-18. Переходная функция системы с вещественной ча­ стотной характеристикой, изо­ браженной на рис. 5-16,а, вы­ численная методом трапецеи­ дальных характеристик.

268

Йо A-функции Для х = 0 находим переходные функции, соответствующие второй трапеции x2(t).

Результаты рабчетов приведены в табл. 5-1, где т, Л(Х_ 0) и

Л ^—0 ,5 ) взяты из таблиц Л-функций:

 

х і ~ ^ Ф о і А ( х _ о ,5)> х г ~

^ФогЛ(-/==о) >х з — ^ Ф о з ^ —0,5)’

 

 

 

 

Т

 

 

Z

 

Z

 

 

 

 

t l ~

«оі

“ «02 И

І з

~ «03 ’

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а 5-1

 

 

ю4

 

*.«)

 

Х,{()

 

 

 

 

 

 

 

 

 

 

 

 

т

Л(Х=0)

О

 

 

 

 

 

 

 

 

II

t,

 

Хі

 

 

-'•a

 

*3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

0

0

0

0

 

0

0

— 0

1

0,31

0,461

0 ,2 5

0,691

0 ,6 6 7

— 0,062

0,1

— 0,138

2

0 ,5 7 2

0,833

0 ,5

 

1,25

1,33

 

— 0,114

0 ,2

— 0 ,2 5

3

0,755

1,061

0 ,7 5

 

1,591

2

 

— 0,151

0 ,3

— 0,318

4

0 ,8 5 7

1,142

1

 

1,713

2 ,6 6

 

— 0,171

0 ,4

— 0,343

5

0 ,8 9 6

1,118

1,25

 

1,677

3 ,3 3

 

— 0,179

0 ,5

— 0,3 3 4

6

0,903

1,051

1,5'

 

1,577

4

 

— 0,181

0 ,6

— 0,315

7

0 ,9 0 4

0 ,9 9 3

1,75

 

1,490

4 ,6 7

 

— 0,181

0 ,7

— 0 ,2 9 6

8

0,91

0 ,9 9 6

2

 

1,449

5 ,3 3

 

— 0,182

0 ,8

— 0 ,2 9

10

0,939

0 ,9 8 2

2 ,5

 

1,473

6 ,6 7

 

— 0,188

1,0

— 0,299

15

0,956

1,005

3 ,7 5

1,508

10

 

— 0,191

1,5

— 0 ,3

17

0,965

1,012

4 ,2 5

 

1,518

11,3

 

— 0,193

1,7

— 0,303

21

0,969

0,995

5 ,2 5

 

1,492

14

 

— 0 ,1 9 4

2,1

— 0 ,2 9 7

25

0,975

1,0

6,25

1 ,5

16,7

 

— 0,195

2 ,5

— 0 ,3

Построив кривые Х\((), x2(t) и x3(t) переходных процессов для составляющих трапецеидальных характеристик и просуммировав их ординаты графически (рис. 5-18), найдем кривую переходного про­ цесса системы:

*(0=*і(0+*2(0+*зС0.

имеющей вещественную частотную характеристику, вид которой представлен на рис. 5-16,а.

5-7. ПРИМЕРЫ ОЦЕНКИ КАЧЕСТВА АСР

Пример 5-1.

О п р е д е л и м п а р а м е т р ы

н а с т р о й к и АСР

т е м п е р а т у р ы

с у ш и л ь н о г о ш к а ф а

(см . р и с .

1-5 и

п р и м е р 2-1),

о б е с п е ч и в а ю щ и е з а д а н н у ю с т е п е н ь

у с т о й ч и в о с т и а и с т е п е н ь к о л е б а т е л ь н о с т и т.

При требуемой степени устойчивости а, заменив в характеристи­

ческом уравнении

(4-45) замкнутой-системы р на — а + /со ,

получим:

ТщТв(—а+/со)3+ (Гщ+Тп) (—'а+;'«)2—

— а + / м + А О б А р = 0 .

269

Соседние файлы в папке книги из ГПНТБ