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

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

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

где: - постоянная интегрирования, определяющаяся начальными условиями и возмущением; - корни характеристического уравнения являются комплексными:

 

 

 

 

 

 

 

 

 

, +1

=

±

 

 

 

 

 

(3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Каждая пара корней дает в выражении (2) составляющую переходного

процесса, равную:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Сi e

i j i t

Ci 1e

i j i t

e

it

Ci e

j it

Ci 1e

j it

/

e

it

sin it i

(4)

 

 

 

 

 

 

 

 

 

 

 

C

 

 

 

C

 

 

 

 

 

 

arctg

Ci 1

.

 

 

 

 

 

 

 

 

где:

 

C 2

C 2

,

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

i

 

i 1

 

 

 

 

 

 

Ci

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Как видно соотношение (4) представляет собой синусоиду (частота i) с изменяющейся экспоненциально амплитудой. Если i <0, то составляющая затухает. Если хотя бы одна составляющая из всех n, где n - порядок дифференциальной системы уравнений, будет иметь i >0, система неустойчива. В частном случае, при i =0 , i= i+1=ai переходный процесс носит апериодический характер. В соответствии с этим, для линейных систем будем различать апериодическую и колебательную устойчивость.

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

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

На рис. 5 представлены две одноконтурные структурные схемы регулирования двух разных параметров на некотором, например, тепловом объекте.

а)

б)

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

11

Предположим для схемы на рис.5.а известны передаточная функция (ПФ) разомкнутой системы Wt (s) от регулирующего воздействия D (например, изменение расхода топлива) к регулируемому параметру T (изменение температуры теплоносителя), а также передаточная функция регулятора в обратной связи t (s). Пусть ПФ и закон регулирования будут простейшими:

Wt

 

 

At

;

t k0 ,

т.е. D k0 T .

(5)

 

 

 

sTt

 

1

 

 

 

 

Корень характеристического полинома нерегулируемой системы λ (знаменателя ПФ разомкнутой системы Wt (p)) полностью определен величиной постоянной времени Tt и не может быть изменен:

W

 

 

At

p 1 .

(6)

 

 

t ( раз)

1

sTt

Tt

 

 

 

Очевидно, λ имеет отрицательное действительное значение и определяет вид переходного процесса после возмущения как апериодически затухающий по экспоненциальному закону e t . При этом интенсивность затухания зависит от величины λ. Соответственно при замыкании контура регулирования характеристический корень будет определяться из полинома знаменателя замкнутой передаточной функции:

 

 

 

 

 

 

 

 

 

 

 

At

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Wt (s)

 

 

 

 

 

 

 

 

At

 

 

At

 

s

1 Atk0

.

 

W

 

 

 

 

 

 

 

1 sTt

 

 

 

 

(7)

t ( зам)

 

W (s)

 

 

 

 

1 st Atk

 

(1 Atk

 

) sT

 

 

1

 

(s)

 

 

 

 

At t (s)

 

 

 

 

 

Tt

 

 

 

t

 

 

 

 

0

 

0

 

 

 

 

 

 

t

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 sTt

 

 

 

 

 

 

 

 

 

 

 

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

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

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

12

нескольких параметров (или сознательном не учете его) система может сохранять свойство односвязности при управлении или регулировании.

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

Рассмотрим более реальный пример регулирования рассмотренного теплового объекта, когда учитывается взаимное влияние выходных параметров T (изменение температуры теплоносителя) и P (изменение давления теплоносителя). На рисунке 6 представлен регулируемый объект с параллельным воздействием двух контуров.

Рис.6 Регулируемый объект с параллельным воздействием двух контуров Двухконтурную схему (рис.6) можно представить в виде

одноконтурной (рис.7).

13

Рис.7 – Неявное представление внутренних контуров в системе

При этом

W

 

 

Wt

.

 

Wt t

 

 

1

 

Если в системе управления имеется несколько замкнутых контуров с более сложными законами регулирования, например пропорционально –

дифференцирующий,

D (k0 k1s) T ,

то

многопараметрический

характеристический полином включает в качестве искомых настроечных

параметров регуляторов несколько пар коэффициентов

k0i k1i .

 

 

 

Тогда характеристический полином можно записать:

 

 

 

A s n

A

s n 1 f

n 1

(k

0i

, k

) A

 

s n 2 f

n 2

(k

k ) ...

n

n 1

 

 

1i

n 2

 

 

0 1

 

 

(8)

 

 

 

 

 

 

 

 

 

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

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

Существуют две основные схемы учета многосвязности:

1.Представление структурно объекта в виде системы с вложенными контурами регулирования, как это предложено в рассмотренном примере

14

(рис.6-7). В более общем виде многосвязную систему с множеством вложенных контуров можно представить следующим образом (рис.8):

Рис.8 Многосвязная система с множеством вложенных контуров Передаточная функция замкнутой по всем контурам системы по

любому из параметров (в частности по ) вычисляется согласно соотношению (9):

Wз

 

W

 

.

(9)

 

W W WE E

 

1

WU U

 

Вторая схема представление объекта предполагает структуру многомерного объекта, имеющего несколько входных и несколько выходных

переменных, через частные передаточные функции, отражающие

15

собственное и взаимное влияние параметров и контуров регулирования (рис.8). На этом рисунке, в качестве примера, приведен двухконтурный тепловой объект, где:

T - изменение температуры теплоносителя - выходной параметр первого контура;

-S - изменение расхода топлива - входной параметр первого контура;

-P -изменение давления теплоносителя - выходной параметр второго

контура;

- C - изменение температуры наружного воздуха - входной параметр второго контура;

-1 , 2 - передаточные функции регуляторов первого и второго контуров;

-Wij (p) собственные и взаимные передаточные функции параметров первого

и второго контуров регулирования

16

а)

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

17

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

L11

 

 

 

W11 (1 W22 2 ) W12

2W21

 

.

(10)

 

 

 

 

 

 

 

W11 1

W22 2 W12 2W21 1

W11 1W22

 

1

2

 

 

2 Методы определения свойств объекта.

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

1. Аналитическое описание объекта.

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

Составим дифференциальное уравнение системы, образованной двумя баками (рис.1), приняв в качестве входной величины изменение расхода жидкости на притоке G1, а выходной - изменение уровня h2 во втором по ходу жидкости баке. Площади поперечного сечения баков A1 и A

A1 h1 G1 G2

G3

 

 

G1 G10

G1

 

 

 

 

 

 

 

 

 

 

 

 

 

G2 G20

G2

 

 

(1.1)

G1

v1

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

G2

kv2

 

 

 

G3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2gh1

 

 

 

 

 

 

 

 

 

 

 

Где

18

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

3: Расход воды от смесителя в втором баке.10: Начальная величина первого расхода.20: Начальная величина второго расхода.1: Изменине первого расхода.

2: Изменине втрового расхода.

К: Постоянная величина. h:Уровень жидкости.

h10:Желаемая величина уровня в смесителе.

Для линеаризации дифференциального уравнения воспользуемся рядом Тейлора:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2g

 

 

h1 ,

 

 

 

 

 

 

 

 

 

 

 

 

2gh1

 

2gh10

 

 

 

 

 

 

 

 

 

 

(1.2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

h10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h1

 

 

1

 

 

(G G

 

G

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

A1

1

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

(1.3)

 

 

h1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g

h

 

 

 

 

 

 

(G G G

20

G

2

2gh

 

 

 

 

 

 

 

 

 

 

 

t

 

A1

10

 

 

 

 

1

 

 

 

 

 

 

 

 

 

10

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2h10

 

При ∆t→0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dh1

= R − Ch ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.4)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Где R =

 

v1+kv2

,

С =

 

 

 

 

 

√g

 

.

 

 

 

 

 

 

 

 

 

A1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A1√2h10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

(s) =

R

,

(1.5)

 

1

s+C

 

 

 

 

 

Для второго бака аналогично получим:

19

 

A h

 

 

G

 

G

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

3

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.1.6)

G4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2gh2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для

 

 

 

линеаризации

нелинейного

 

 

 

дифференциального

уравнения

воспользуемся рядом Тейлора:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2g

 

h2 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

2gh2

 

2gh20

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.6)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

h20

 

 

 

 

 

 

 

 

 

 

 

 

h2

 

 

1

 

(G

 

G

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

3

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

(1.7)

 

h2

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g

h

 

 

 

 

g

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( 2gh

 

2gh

20

 

h

2

)

 

 

 

 

 

 

 

 

 

t

 

 

 

A2

 

 

 

 

 

10

 

 

 

1

 

 

 

 

2h20

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2h10

 

 

 

 

 

 

 

 

 

 

При ∆t→0

dhdt2 = Dh1 − Eh2,

(1.1.8)

Где:

 

 

 

 

 

 

 

 

 

 

 

1

 

g

 

√g

D =

 

 

,E =

 

 

,

A2

2h10

 

 

A2√2h20

 

Воспользовавшись дифференциальным уравнением для контрольного (второго) бака, и применяя преобразование Лапласа получим :

H2(s)

E

 

 

 

=

 

 

,

 

(1.9)

H1(s)

s+D

 

Передаточная функция объект регулирования является H2(s).

 

H2(s) =

 

 

E R

,

(1.10)

 

 

 

 

(s+D)(s+C)

 

 

 

 

 

 

Составим дифференциальное уравнение системы управления по температуре (рис 1.):

 

 

 

=

(− +

 

+

),

(1.11)

 

 

 

 

 

1

1

2

2

 

 

 

 

 

 

 

 

 

 

 

= − + + ,

 

 

(1.12)

 

 

 

 

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

 

= ,

 

 

 

(1.13)

 

 

 

 

 

 

 

 

 

20