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

Системы управления (курсовой проект) / Задание на курсовое проектирование СУ_2015

.pdf
Скачиваний:
14
Добавлен:
25.09.2019
Размер:
2.49 Mб
Скачать

E: Энергия жидкости.

m:Полная масса жидкости в смесителе.

Т: Температура жидкости в смесителе .

ср: Теплоемкость жидкости.

Т1 : Температура первого расхода жидкости.

Т2 : Температура второго расхода жидкости.

1: Расход холодной воды.

2: Расход горячей воды.

3: Расход воды от смесителя в вторую баку

Дифференцируя уравнение (1.1.13)

и получим:

 

= с (

 

+

 

),

(1.14)

 

 

 

 

 

 

 

При условии, что управление температурой происходит при постоянном (зафиксированном) уровне получим, что уравнение (1.12) равно нулю и m (масса жидкости в баке) постоянная величина, следует:

 

(− +

 

+

) = с

 

 

,

(1.15)

 

 

1

1

2

2

 

 

 

 

 

 

 

 

 

 

В данном случае температура на входе бака 1 = 0, и при постянном уровне,и постянные величины. Применяя преобразование Лапласа, получим:

T(s)

 

=

2

.

(1.16)

( )

 

 

+

 

2

 

 

 

 

 

Составим дифференциальное уравнение перекрестных каналов 12 и 21 (рис.1). Чтобы найти зависимости температуры от расхода 1, используем

уравнение (1.1.15) и поставим вместо 2 = − 1

и получим

 

 

+ = ( − ).

(1.17)

 

 

 

 

 

2

1

 

 

 

 

 

 

 

 

 

 

 

Продифференцировав это уравнение, получим:

 

2

 

 

 

 

 

 

 

 

+

 

=

1

.

(1.18)

2

 

 

 

 

2

 

 

Предположим, что все начальные условия равны нулю применяя преобразование Лапласа, получим:

21

2 ( ) + ( ) =

( ),

(1.19)

 

 

 

2

1

 

( )

 

=

2

.

 

(1.20)

( )

 

 

 

+

 

 

Чтобы найти зависимости уровня в баке от расхода 2, используем уравнение (1.1.3). В этом случае мы хотим найти отношение, которое описывает зависимость уровня в баке от второго расхода, поэтому можно считать, что первый расход постоянен, и тогда получим 1 = 0:

 

 

h1

 

1

 

(G G

 

 

G

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

t

A1

1

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

(1.21)

 

 

h1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g

h

 

 

 

 

 

(G G

20

G

2

2gh

 

 

 

 

 

 

 

 

 

 

t

A1

10

 

 

 

 

 

 

 

 

10

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2h10

 

 

 

 

 

 

Устремив ∆t→0, получим:

 

 

 

 

 

1

= − С ,

 

 

 

 

 

 

 

 

 

 

 

(1.22)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где =

√g

 

 

, =

 

1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A1√2h10

 

 

 

 

1

 

 

 

 

 

 

 

Применив преобразование Лапласа к выражению (1.22) получим:

1( ) = .2( ) +

Подставив выражение (1.23) в (1.9) получим:

H2(s)

=

 

E

=

 

 

2( ),

H1(s)

s+D

( + )( + )

H2(s)

=

 

 

 

.

 

( ),

( + )( + )

 

2

 

 

 

 

 

 

 

Снятие кривой разгона объекта управления.

(1.23)

(1.24)

(1.25)

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

Снимают кривую разгона следующим образом: объект приводят в равновесное состояние, при котором все входные и выходные величины постоянны. Затем быстрым перемещением регулирующего органа (клапана, заслонки, регулятора мощности и т.д.) вносят скачкообразное возмущение.

22

После этого записывают периодически результаты измерения выходной величины до тех пор, пока выходная величина не примет нового установившегося значения или пока не установится постоянная скорость ее изменения. По точкам отсчета строят кривую в координатах: выходная величина - время, которая и является кривой разгона. Для снятия кривой разгона удобно использовать самопишущие регистрирующие приборы с ленточной картограммой.

Различные возможные виды кривых разгона изображены на рис.10 .

Рис.10 Кривые разгона.

а) - без точки перегиба и с ненулевой начальной скоростью;

б) - с точкой перегиба и нулевой начальной скоростью; в) - с чистым запаздыванием; г) - для объектов без самовыравнивания.

Аналитическое выражение для кривой разгона соответствует переходной функции, которая определяется известным преобразованием Лапласа:

Y s

1

c j

W s

e

st

 

 

 

 

ds,

(2.1)

 

 

 

 

2 c j

 

s

 

 

где W(p) - передаточная функция объекта.

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

23

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

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

-метод площадей;

-метод дополнительных членов;

-метод последовательного логарифмирования;

-аппроксимация дифференциальными уравнениями;

-аппроксимация суммой элементарных звеньев;

-графический метод.

1.Определение передаточной функции одноемкостных объектов.

Объекты одноемкостные описываются дифференциальным уравнением первого порядка:

T

dY

Y KXвх .

(2.2)

 

 

dt

 

 

 

Передаточная функция такого объекта:

 

W s

K

.

(2.3)

 

Ts 1

 

 

 

 

 

В этом случае передаточная функция полностью определена, если известны коэффициент усиления К и постоянная времени Т. Решение уравнения (26) представляет собой экспоненту и имеет вид:

 

1 e

t

.

 

Y KX

T

(2.4)

 

вх

 

 

 

 

 

Таким образом, если экспериментально полученная кривая разгона хорошо аппроксимируется экспонентой, то непосредственно по этой кривой можно получить параметры k и Т (рис.11).

24

Рис.11 Кривая разгона системы первого.

Для описания астатического объекта достаточно знать время запаздывания и установившуюся скорость изменения выходной величины К (рис.1.1г). Передаточная функция объекта без самовыравнивания.

W s

Ke s

.

(2.5)

 

 

s

 

3. Определение параметров передаточной функции системы второго порядка.

Для кривой разгона (рис. 10 - б) объект относится к системе управления второго порядка и описывается дифференциальным уравнениям соответствующего порядка. В первом приближении можно попытаться описать данную кривую дифференциальным уравнением второго порядка,

T T

d2Y

T T

 

dY

Y KX

 

.

(2.6)

 

 

 

1 2 dt2

1 2

 

dt

 

вх

 

 

Здесь задача сводится к определению постоянных времени T1 и Т Решение дифференциального уравнения будет иметь вид:

 

 

T1

 

t

 

 

Y KX 1

 

e

T1

 

 

 

 

T1 T2

 

 

 

 

T2

T1 T2

 

t

 

 

 

 

e

T2

.

(2.7)

 

 

 

 

 

 

 

Для определения значений T1 и Т2 проводят касательную к кривой разгона в точке перегиба А (рис.12).

25

Рис.12 Кривая разгона для системы второго порядка.

Величина под касательной ВС равна сумме постоянных времени

BC T1 T2 .

 

 

(2.8)

Величину Т2

можно определить из выражения:

 

 

KX

 

 

 

ta

 

T

вх

T

 

 

 

e

2 .

(2.9)

 

 

2

tg

 

 

 

 

 

 

 

 

 

 

Это выражение трансцендентно, поэтому прямо определить Т2 нельзя. Оно может быть найдено как абсцисса точки пересечения двух функций:

t1 T2 ln T2 C

 

 

 

T ta

 

t

 

(2.10)

 

2

2

T2

 

 

 

 

 

 

26

C ln

tg

(2.11)

KXвх

 

 

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

27

3 Методы настройки двухсвязных систем регулирования

Из общего числа систем регулирования около 15% составляют двухсвязные системы регулирования (рис.13). В таких системах даже при наличии устойчивой автономной работы двух регуляторов вся система может стать неустойчивой за счет действия перекрестной связи в объекте управления.

рис.13 Схема двухсвязной системы управления

Объект управления в двухсвязной системе представлен в Р канонической форме. Удобство такого представления заключается в том, что путем активного эксперимента можно определить все передаточные функции по соответствующим каналам. Промежуточные сигналы x 1 , x 2 , x 3 , x 4 обычно недоступны для измерения, поэтому управление ведется по вектору

выхода Y:

 

 

1

],

(3.1)

= [

2

 

 

На практике довольно большое число систем являются двухсвязными. Для объективной настройки регуляторов двухсвязных систем формируется критерий качества вида:

0 = 1 1 + 2 2,

(3.2)

28

где 1 и 2 — коэффициенты веса (штрафа), J1 и J 2 — критерии качества первого и второго контуров.

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

Задача настройки регулятора состоит в том, чтобы при заданных 1 и2 обеспечить минимальное значение J 0 системы, где

 

= ∫2

( ), = 1, .,

(3.3)

 

0

 

 

 

 

 

 

 

Рассмотрим различные методы настройки регуляторов в двухсвязных системах.

3.1Метод автономной настройки регуляторов

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

регулятор Р 2 переводится в ручной режим работы;

настраивается регулятор Р 1 так, чтобы критерий J 1 был минимален;

отключается настроенный регулятор Р 1 и включается регулятор Р 2 ;

настраивается Р 2 , обеспечивая минимум J 2 ;

оба регулятора включаются в работу.

Такой подход рекомендуется использовать если:

наблюдается малое взаимное влияние контуров;

быстродействие одного контура значительно выше другого (контуры разнесены по частотам);

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

29

3.2 Метод итеративной настройки регуляторов

Этот метода аналогичен предыдущему, но здесь осуществляется многократная настройка регуляторов Р 1 и Р 2 (последовательная подстройка) с целью обеспечения минимального значения критерия качества J 0 всей системы.

Следует учитывать, что только метод итеративной настройки регуляторов обеспечивает качественную работу двухсвязной системы даже при наличии сильных перекрестных связей. Это объясняется тем, что оптимизация критерия качества J 0 системы происходит при включенных Р 1 и Р 2 .

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

3.2 Метод аналитического конструирования регуляторов

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

Структурная схема оптимального регулятора состояния, содержащего наблюдающее устройство, приведена на рис.14. Схема содержит следующие элементы: Н — наблюдатель, ОУ — объект управления, МОУ — модуль объекта управления, ОРС — оптимальный регулятор состояния, Е Н — ошибка наблюдения, X М — вектор состояния модели, X зад.— вектор задания, U — вектор входа ОУ, Y — вектор выхода ОУ, Y М — вектор выхода модели.

30