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

книги из ГПНТБ / Толшин В.И. Основы автоматики и автоматизации энергетических установок учебник

.pdf
Скачиваний:
21
Добавлен:
25.10.2023
Размер:
14.19 Mб
Скачать

так как все члены, входящие в Q(со), содержат со; следовательно, Q(0) =0. Поэтому

y(t) = P ( 0)

Ф ( » - Я ( 0 ) е м йш

>

 

 

— во

Произведем преобразования подынтегрального выражения,

учитывая, что ejmt = cos at+j sin at:

£

1

= | p c ) + ; q w -

p (0)] (cosa< + • S I n

_

 

J®

 

>

 

 

 

_

[Р(ш) — P(0)]COSCD#

Q (id) COS <o£

[P(m)—P(0)] Sin at

|

 

j (O

 

0)

Ш

 

 

 

. Q(co)sinco£_[P(co) — P(0)] sinm^-f- Q (со) cos at

 

 

 

CO

 

 

0)

 

 

 

 

Q (со) sin co£ — [P (u>) — P (0)] cos at

 

 

 

'

 

(0

^ "

 

 

 

Учитывая проведенные преобразования, получаем

 

 

 

 

Q (со) cos at -j- \P (со) — P (0)] sin at ^

 

 

,

W _ p ( 0 , +

^ J

 

 

 

a

 

 

 

 

 

 

 

 

 

j

F \P (c u )

P ( 0)] C O S at — Q ( с о ) sin at ^

(6.13)

 

2~

|

со

 

 

 

 

 

Подынтегральная функция 1-го интеграла четная, поэтому его нижний предел можно взять равным нулю и удвоить значение интеграла. Подынтегральная функция 2-го интеграла нечетная, поэтому второй интеграл обращается в нуль. Таким образом,

cos » tdw + ±

f P (« )sin * f_ d<D_

Ш . J JSH ^ d a .

(6.14)

Известно, что последний интеграл выражения (6.14) равен я/2, поэтому

» <Я - 4 Р (0)

+

Л + У j Р ( ^ ± Л. (6.15)

170

Если принять нулевые начальные условия, то до приложения внешнего воздействия (при t<0) y (t ) —0. Заменив в (6.14) t на —t, получим

 

‘Q (со) cos at

 

оо

 

 

o-4-P(0) + i- |

d<o

Р(со) Sin o)^

da. (6.16)

О)

те

СО

Вычтем (6.16) из (6.15), в результате получим

 

 

y(t) =

h (t) =

Р (®) sin at ^

 

(6.17)

(1)

 

 

 

 

 

 

 

О

Следовательно, если входное воздействие представляет собой единичную ступеньку, то переходная функция системы может быть получена по известному выражению ВЧХ Р(со) с помощью фор­ мулы (6.17).

Однако вычисление h{i) при реальном криволинейном харак­ тере зависимости Я (со) сопряжено со значительными трудностями. Поэтому обычно используют приближенное решение задачи. Это

решение 'основано на том, что переходная

функция (6.17) может

быть относительно просто рассчита­

 

 

 

 

на,

если ВЧХ

имеет

вид

прямо­

 

 

 

 

угольной трапеции.

 

 

 

 

 

 

 

 

Прямоугольная трапеция (рис.

 

 

 

 

6.9)

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

следующими

 

 

 

 

величинами:

значением Р(0)

(при

 

 

 

 

ю=0), частотой

 

среза

соП)

коэффи­

 

 

 

 

циентом наклона

x=a>J mп- Найдем

 

 

 

 

значение h(t)

 

для

прямоугольной

 

 

 

 

трапеции как функцию от величин

 

 

 

 

Я(0), %и соп.

 

 

 

 

 

 

 

 

 

 

 

Очевидно, что на участке 0—шя

 

 

 

 

 

Р И

=

Р(0),

 

 

Рис.

6.9.

Прямоугольная трапеция

на участке соя — шп

 

 

 

(D

-

а

\

 

 

 

 

 

 

 

 

Тогда

 

 

Я (со) = Р (0 ) 1 шп —

 

J

-

 

 

 

 

 

 

 

"п

 

 

 

 

h (t) = — 1 — — sin

 

р (0)

 

 

2

г

 

 

. . . .

 

2 Г*

Р (

0) .

 

 

 

ш

 

 

f

 

j

 

 

 

7ZJ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

Я(0)(со — сда) sin atdu).

 

 

 

 

 

 

те

ш

«>К

— <°а)

 

 

 

 

 

 

 

 

 

 

а

 

 

 

 

 

 

171

Рис. 6.10. Замена истинной ВЧХ прибли­ женной

Введем понятие интегрального синуса si («о = Р ^

инового времени т=£шп. Получим

Л[Я(0), X. ■=] =

# - Я (0 )(81(Хт)Ч 1

C O S t —

C O S y T

. (6.18)

 

 

i — X si СО—si(xT) +

 

 

 

Так как значения интегрального синуса определяются числен­ ным методом, то выражение (6.18) при Я(0) = 1 вычислено заранее для различных значений х и т (приложение 2).

Для получения переходного процесса в реальном времени и для случая, когда Я (0)^1, необходимо перестроить переходную функ-

цию /г(т, х) подстановкой

х

табличного зна­

t = — и умножением

чения Л(т,

X) на Р (0).

Шп

 

Однако

реальная ВЧХ отличается от прямоугольной трапеции.

Поэтому ее представляют

не одной трапецией,

а в виде суммы

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

больше число отрезков у ломаной, т. е. чем точнее ломаная линия аппроксимирует реаль­

ную зависимость Я (со). Затем Я (со) в виде ло­

маной линии заменяют трапецеидальными харак­ теристиками по следую­ щим правилам:

а) вертикальная боко­ вая сторона каждой тра­ пеции должна совпадать с осью ординат ВЧХ;

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

На рис. 6.10 показана замена реальной ВЧХ суммой трапеце­ идальных характеристик. Первоначально кривая заменена ломаной линией P(0)bcdef. Затем строятся трапеции в соответствии с изло­

женными правилами.

 

Участок ломаной Р(0)Ьст заменяется двумя

трапециями:

Og\h = aP(0)6 (в данном случае — треугольник) и

Опгса. Сумма

172

ординат этих трапеций равна ординатам ломаной Р(0)Ьст. Участок ломаной mdef заменяется трапециями Ogef и Ogdm.

Так как трапеция Опгса, а также Ogdm в сумме составляют трапецию aiOd\cx, то окончательно ВЧХ заменяется тремя

трапециями:

1— Og\h; 2 — Ogef; 3 — aiOd\C\.

Сумма ординат

трапеций 1, 2

и 3 при любом значении

со равна

ординате лома­

ной P(0)bcdef при том же значении со.

вычисленных для всех

Суммирование переходных функций

трапецеидальных характеристик, можно производить графически. В этом случае при построении hL(t) следует учитывать масштабы времени.

Пример 6.2. В качестве примера рассмотрим определение показателей пере­ ходных процессов в системе регулирования частоты вращения дизель-генератора с регулятором непрямого действия и упруго присоединенным катарактом, схема которой была представлена на рис.. 6.4. Принцип действия регулятора был изло­

жен в примере 6.1. Постоянные уравнения составляют: Та = \,Ъ2 сек,

7^=0,025 сек,

У; =0,037 сек; 6=0,03; Зс =0,04. Требуется

определить показатели

переходных

процессов при сбросе 100% нагрузки.

 

 

Для этой цели составляем передаточную

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

Ф ( р )

Из у

-угсб(Р)

Й з Т

1 + W06(p) w p (p) ■

 

Раскрывая W0c,(p) и

(р),

получаем

 

Ф ( р )

 

Айр* -\г

р \

В0р 3 + BlPt + BiP + 1/8 ’

 

где A0= T iT sb^lo; A^Tfi^jb;

B0= T aTi TsbY:lo;

5 i= TiTab^lb+TsTa; B ^ T a + T J o ;

Подставляя в эту формулу /со и учитывая, что исследование ведется для сброса нагрузки (Х=—1), получаем выражение ВЧХ

 

р . . _

С0и)4+С]Ц)-+1/5

 

 

К ’

Соив+СзоН +С^-И /й2

где Cq= A qBiAjД,у С\—А^В^—Ао/о—Дц

С^—В^,

С^=В^

27ДAq; С^=В^—2Дуо.

ВЧХ САР

представлена на рис. 6.11, а, а отдельные прямоугольные трапе­

ции, на которые она разбита, — на рис. 6.11,6.

 

 

Трапеция 1 соответствует Р(0)аЬВ, yi=0,27 шп,=15

Яг (0 )= —0,026.

Трапеция 2

соответствует

P(0)dcB,

у2=0,26

шп,= 2 0

Л (0)=0,026.

Трапеция 3

соответствует

P(0)Eed,

у3=0,77

а>п3= 26

Я3 (0)=0,022.

Трапеция 4

соответствует

EOfe,

=0,68

мП4=38

Я, (0)=0,008.

По таблицам приложения

2 определяются /г-функции, соответствующие зна­

чениям единичных трапецеидальных характеристик Л/ (xi, т) с величинами xi, Х2>

Хз, Xj при разных т. Значения

Лг- {yj, z)

вписываются

в

столбцы

2, 5, 8 и II

табл. 6.2.

В столбцы 3, 6,

9

и

12

этой

таблицы

вписываются

произведения

Р/(0)/ц(х,-,

т),

а в столбцы

4,

7,

10

и 13

— значения истинного

времени t, пере-

считанные

по

 

Т

 

 

 

 

 

 

 

производится

формуле t ——

. Суммирование функций Л/Я,-(0)

графическим путем. На рис. 6.12

представлены функции Л,-Я/ (0)

и функция h(t).

По

графику определяем показатели

переходных

процессов:

нерегулирование

°п ~

0,015, время переходного

процесса до вхождения

в

зону

нестабильности

±0,0027п s: 0,54

сек.

 

 

 

 

 

 

 

 

 

 

173

а)

Рис. 6.11. Разбиение ВЧХ на отдельные прямоугольные трапеции

Рис. 6.12. Построение приближенной переходной характеристики

174

Т а б л и ц а 6.2

 

 

 

 

 

Расчет Л-функций

 

 

 

 

 

Трапеция

1

Трапеция 2

Трапеция 3

Трапеция 4

 

h\

Ро(0) hi

t

h.

Р2(0) hn

t

 

Рз (0) As

t

К Л (0 )Л 4

t

1

2

3

4

5

6

7

8

9

10

11

12

13

1,0

-0 ,3 9

-0,010

0,067

0,56

0,014

0,65

0,53

0,011

0,038

0,52

0,0042

0,026

2,0

-0,71

-0,018

0,133

0,97

0,025

0,10 0,94

0,021

0,077

0,92

0,0074

0,053

4,0

—1,04

—0,027

0,270

1,15

0,030

0,20 1,16

0,025

0,154

1,16

0,0043

0,105

6,0

—1,05

-0,027

0,400

0,93

0,024

0,30

0,96

0,021

0,23

0,98

0,0078

0,158

8,0

—1,02 -0,027

0,530

0,96

0,025

0,40

0,94

0,021

0,31

0,94

0,0075

0,216

10,0

—1,03

—0,027

0,670

1,06

0,027

0,50

1,04

0,023

0,38

1,02

0,0082

0,26

12,0

—1,02 —0,027

0,800

1,00

0,026

0,60

1,03

0,023

0,46

1,03

0,0083

0,310

§ 6.4. Косвенные методы оценки качества переходных процессов

Приблиоюенная оценка вида переходного процесса по вещественной частотной характеристике

Можно выделить следующие основные свойства ВЧХ, по ко­ торым оценивается вид переходного процесса:

а) Приблизительно одинаковым вещественным частотным ха­ рактеристикам соответствуют приблизительно одинаковые по виду переходные процессы.

б) Если имеются две сходные вещественные характеристики, различающиеся в основном значениями и>п, то переходный процесс, определяемый ВЧХ, у которой ве-

личл-ша шп

‘больше, протекает быст­

рее, чем переходный процесс, опре­

деляемый

ВЧХ,

у которой

величи­

на соп

меньше.

Это положение

сле­

дует из выражения

(6.18)..

Таким

образом,

длительность

переходного

процесса

будет

тем

меньше,

чем

более полого идет ВЧХ.

 

функ­

в)

Значение

переходной

ции при 1 = 0 и

t =

оо

может быть

установлено по соотношениям

Рис. 6.13. ВЧХ, имеющая

y t-o = p (m)

и У/=«, = -р (и>)-

разрыв

 

 

1 0 = со

 

 

СО-я 0

Это вытекает из теоремы о начальном и конечном значении ориги­ нала.

г) Если при каком-либо <b = coi, отличном от нуля, ВЧХ стре­ мится к бесконечности и в ней имеется разрыв непрерывности (рис. 6.13, кривая 1), то это означает, что система регулирования

175

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

Если ВЧХ стремится к бесконечности при со = 0 (рис. 6.13, кри­ вая 2), то это означает, что характеристическое уравнение имеет нулевой корень и система находится на границе устойчивости.

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

d P ™ <0; Р {ш)> 0 при всех ш, то величина перерегулирования не

превышает 18% статического отклонения.

е) Часть вещественной частотной характеристики, в которо ее ординаты становятся меньше 0,1—0,2 начального значения Р(0), можно не принимать во внимание, так как это приводит к несу­ щественным погрешностям при анализе переходной функции в ее начальной стадии.

Пример 6.3. Требуется оценить качество регулирования САР частоты враще­ ния, показанной на рис. 6.4, при дополнительном оборудовании регулятора уст­ ройством измерения активной мощности. Сигнал измерителя активной мощности

Рис. 6.14. Структурная схема САР с комбиниро­ ванным принципом регулирования

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

Принято, что изменению нагрузки на 100% (А,= I) соответствует перемеще­ ние хода рейки на 50% (г2=0,5). Постоянные времени обоих сервомоторов оди­ наковы: Т. — Т' .

176

С о ст а в и м п е р ед а т о ч н у ю ф у н к ц и ю за м к н у т о й си стем ы Ф к (р):

где Ф(р) — передаточная функция замкнутой САР, реализующей регулирование только по отклонению регулируемой величины.

Ошибка САР на установившихся режимах при & =0,5 составит wKmax= =0,5Ф (0), где Ф '(0)=итах=5=0,03.

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

Р

0,03

О

 

Ю

20

О) 30

 

Рис. 6.15. ВЧХ замкнутых систем:

-----------------

 

ВЧХ САР по отклонению регулируемой

величины;

 

ВЧХ САР

с дополнительным импуль­

 

 

сом по возмущению

 

Выражение для ВЧХ получим в виде

Р к (ш) = Р(ш) — kPap И + ^рJ s ^ Q W

1 +

где Р(со) — ВЧХ САР по отклонению регулируемой величины.

Так как величина Ts^ мала по сравнению с величинами. Та и ■7’,- (см. при­

мер 6.2), то, пренебрегая ею,

получаем Р к (а>) = 0ДР(ш).

можно судить о том,

На рис. 6.15 даны ВЧХ

обеих САР. По виду ВЧХ

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

Частотные критерии

Частотные критерии качества определяются по виду ампли­ тудно-фазовой характеристики разомкнутой системы, и поэтому их удобно использовать в случае применения частотных методов исследования.-Для этой цели вводятся понятия: запас устойчивости по модулю (амплитуде) и запас устойчивости по фазе.

12 В. И. Толшнн

177

Запасом устойчивости по модулю называется величина, опре­ деляемая в общем случае для условно устойчивой САР по выра­ жениям

L j = 20 lg — дб и Z-2= 20 lg и2 дб,

 

 

 

 

и1

 

 

 

 

 

 

 

где «1 и «2 — расстояния от начала

координат до точек пересече­

ния АФХ с отрицательной полуосью абсцисс

(рис. 6.16).

 

 

 

Запасом устойчивости по

 

фазе

'называется

величина

 

I-L=180 -Ьяр, где

 

ф — аргу­

 

мент

частотной

передаточ­

 

ной функции при модуле,

 

равном

единице.

ц

может

 

Геометрически

 

быть

получена

следующим

 

путем.

 

Строят

окружность

 

радиусом, равным единице,

 

и €

 

центром

в

начале

 

координат. Угол между от­

 

рицательной полуосью

абс­

 

цисс и радиус-вектором,

 

направленным

в

точку

пе­

Рис. 6.16. К определению запаса устойчи­

ресечения окружности

еди­

вости по АФХ

ничного

радиуса

с

АФХ,

 

есть ц.

 

 

 

 

 

Запасы устойчивости по модулю и фазе

могут быть косвенно

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

В хорошо демпфированных системах запас устойчивости по мо­ дулю составляет 6—20 дб, а запас устойчивости по фазе 30—60°.

Интегральные оценки

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

Наиболее простая интегральная оценка качества

оо

(6.19)

h = \ y d t ,

о

 

где y=y(t) — отклонение регулируемой величины от нового уста­ новившегося значения (у/=„ = 0).

178

интерпретация интегра­ льной оценки качества/]

Геометрически величина / j соответствует площади под кривой переходного процесса (рис. 6.17). Чем меньше площадь, тем быст­ рее затухает переходный процесс, тем меньше величина перерегу­ лирования и выше точность регулирования в переходном режиме. Однако для колебательного процесса (рис. 6.18) суммирование положительных и отрицательных площадей может при­ вести к условию /1= 0, в то время как качествр регулирования будет низким.

Поэтому более предпочтительна другая, так называемая «квадратичная», оценка качества

оо

(6.20)

/« = J y 2dt.

о

■ Рис. 6.17. Графическая

л а оценка не зависит от знаков отклонений и, следовательно, от формы

процесса. Однако величина /2 также не зависит от частоты коле­ баний.

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

ния может быть использована улучшенная квадратичная инте­ гральная оценка качества / 3:

 

 

h

= \

(У2+ T * y ' 2)dt, (6.21)

 

 

4

О

 

 

 

 

 

 

где у ' — .производная

регулируе­

 

 

 

мой

величины,

которая

 

 

 

тем больше, чем больше

 

 

 

частота колебаний;

 

Рис. 6.18. График колебатель­

Т — постоянная

времени.

 

ного переходного

процесса

Существуют методы

расчета

ин­

зволяют избежать

 

тегральных оценок, которые по­

непосредственного

интегрирования.

Так,

на­

пример, величина /] может быть определена следующим образом:

I\ — \i> (t)dt =

\[m Г у {t)e

pidt = lim Y(p),

О

^-° 0

p-*-0

где У(p) — изображение по Лапласу величины у.

Подобным образом могут быть получены выражения для /2 и /3.

12*

179

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