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

Синтез систем автоматического управления в биотехнологии методом стандартных коэффициентов

.pdf
Скачиваний:
0
Добавлен:
09.08.2026
Размер:
454 Кб
Скачать

Решетняк Е.П., Рыбин М.Б., Лямина Т.В.

СИНТЕЗ СИСТЕМ АВТОМАТИЧЕСКОГО УПРАВЛЕНИЯ В БИОТЕХНОЛОГИИ МЕТОДОМ СТАНДАРТНЫХ КОЭФФИЦИЕНТОВ

Саратов 2006

1

УДК 62

ББК 32.96

Решетняк Е.П., Рыбин М.Б., Лямина Т.В.

Синтез систем автоматического управления в биотехнологии методом стандартных коэффициентов / Решетняк Е.П., Рыбин М.Б., Лямина Т.В. – Саратов: Саратовский государственный аграрный университет имени Н.И. Вавилова, 2006. – 17 с. — опубл. 15.06.2012. — Электрон. версия печ. публ. — Доступ с сайта ЭБС IPRbooks.

В задачах синтеза САУ в биотехнологии в настоящее время используются два метода: графоаналитический и интегральных оценок. В работе предложена достаточно простая и эффективная методика синтеза САУ биотехнологическими объектами, основанная на использовании метода стандартных коэффициентов. Ее высокая эффективность показана на примерах синтеза двух САУ. Она проста для понимания и может быть рекомендована для использования студентами при курсовом и дипломном проектировании.

2

В в е д е н и е

С и н т е з с и с т е м ы а в т о м а т и ч е с к о г о у п р а в л е н и я ( С А У ) —

это определение структуры, состава элементов, и значений параметров САУ, при которых она удовлетворяет предъявляемым к ней требованиям.

Данное определение характеризует общую задачу синтеза.

Задача определения типа, схемы и параметров закона управления (ЗУ)

или корректирующих устройств (КУ) называется ч а с т н о й з а д а ч е й с и н т е з а С А У .

Разработано большое количество методов решения частной задачи синтеза САУ, в частности, в биотехнологии в настоящее время использу-

ются два метода: графоаналитический [1,2] и интегральных оценок [3].

Ниже для решения задачи синтеза САУ биотехнологическими объектами управления (ОУ) предлагается использовать метод стандартных коэффи-

циентов [4].

3

1. Постановка задачи синтеза САУ методом

стандартных коэффициентов

Большинство (ОУ) в биотехнологии может быть представлено следую-

щей структурной схемой [5]:

Регулирующий

 

 

 

 

 

 

орган

 

 

F(p)

 

G(p)

 

 

 

 

 

Y(p)

 

 

 

 

k

k=1

 

 

 

 

 

 

 

 

 

 

 

 

 

( T1 p 1 )( T2 p 1 )

 

 

 

 

 

 

 

 

 

 

 

 

Параметры передаточных функций ОУ, как правило, получаются в ре-

зультате идентификации на основании экспериментально полученной кри-

вой разгона. Задающее воздействие и возмущение принимаются постоян-

ными, регулирующий орган – безынерционным (его передаточная функция

равна 1).

 

 

Т р е б у е т с я с и н т е з и р о в а т ь

САУ, удовлетворяющую следующим

требованиям: время регулирования tр

tр зад , перерегулирование

зад,

установившаяся ошибка САР при воздействии задающего сигнала

и воз-

мущения должна бать равна нулю; должны быть ограничены так же мак-

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

на при воздействии возмущения.

В методе стандартных коэффициентов структура и параметры регулято-

ра: k0 , k1, k2, … должны быть выбраны так, чтобы передаточная функция замкнутой САУ равнялась желаемой (стандартной передаточной функции

или передаточной функции эталонной модели - ЭМ):

 

WСАУ ( p, k0 , k1 ,...) WЖ ( p ) WЭМ ( p )

(1)

или

BСАУ ( p, k0

, k1 ,...)

 

BЖ ( p )

 

BЭМ ( p )

.

 

 

 

 

 

 

AСАУ ( p, k0

, k1 ,...)

AЖ ( p )

AЭМ ( p )

 

 

 

 

 

 

 

Для решения полиномиального уравнения (1) обычно выписывают сис-

тему уравнений для коэффициентов при одинаковых степенях p. Эта зада-

4

ча будет иметь решение, если структура замкнутой САУ и ЭМ выбраны так, чтобы выполнялись условия «соответствия передаточных функций»

(СПФ):

1) совпадали размерности полиномов числителей и знаменателей пере-

даточных функций замкнутой САУ и ЭМ:

dim BСАУ ( p, k0 , k1 ,...) dim BЭМ ( p ) ,

(2)

dim АСАУ ( p, k0 , k1 ,...) dim АЭМ ( p ) ;

 

2) число n искомых параметров регулятора равнялось числу уравнений,

получаемых из уравнения (1):

dim AЭМ ( p ) dim BЭМ ( p ) n .

(3)

Выбор структуры замкнутой САУ должен обеспечивать также выпол-

нение требований качества и устойчивости. В рассматриваемой задаче синтеза замкнутая САУ должна быть астатической первого порядка по от-

ношению к задающему воздействию и возмущению. Обеспечение астатиз-

ма достигается использованием И, ПИ, ПИД – регуляторов. Для обеспе-

чения заданного времени регулирования необходимо использовать после-

довательно включенный П – регулятор.

В работе выбраны две пары структур замкнутой САУ и ЭМ, отвечаю-

щих перечисленным выше требованиям. Рассмотрим особенности синтеза этих САУ методом стандартных коэффициентов и их свойства.

5

2. Синтез САУ №1

Структурная схема САУ №1 представлена на рис.1.

 

 

 

 

 

 

 

 

 

 

F(p)

G(p)

 

(p)

 

 

 

U(p)

 

 

 

 

 

 

Y(p)

 

 

 

 

 

k

 

 

 

k1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( T1 p 1 )( T2 p 1 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k0

p

1 k2 p

Рис. 1. Структурная схема САУ №1

Р А С Ч Е Т П Е Р Е Д А Т О Ч Н О Й Ф У Н К Ц И И С А У № 1

Передаточная функция замкнутой САУ

 

 

 

( k0 k1 p )k

 

 

 

WСАУ ( p )

 

p( T1 p 1 )( T2 p 1 )

 

 

1

( k0 k1 p )( 1 k2 p )k

 

 

 

 

 

p( T p

1 )( T p 1 )

 

 

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

 

( k0

k1 p )k

 

 

 

 

 

 

 

 

T T p3

( T

T ) p2

p k k k

2

p2 ( k k

2

k )k p k k

 

 

 

 

1

2

 

 

1

2

 

 

 

 

1

 

 

 

0

1

0

 

 

 

 

 

 

 

 

 

 

 

(

k1

 

 

p 1 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k0

 

 

 

 

 

 

 

.

(4)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T T

 

 

p3

 

( T

T k k

k )

p2

[( k k

2

k

 

)k 1]

 

 

 

1 2

 

 

1

2

1 2

 

 

 

 

 

0

1

 

 

p 1

 

 

 

k0k

 

 

 

k0k

 

 

 

 

 

 

 

 

k0k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Как видно, в знаменателе передаточной функции САУ (4) необходимо задавать три коэффициента, для этого введены три передаточных числа. Однако, что делать с параметрами числителя ПФ (2)?

6

В Ы Б О Р Э Т А Л О Н Н О Й М О Д Е Л И ( Э М )

Выбор эталонной модели (ЭМ) осуществляется исходя из требований к качеству замкнутой системы. Так как разомкнутая система имеет порядок равный двум, ЭМ выберем в виде инерционного звена 2-го порядка

 

Wэм

( p )

 

 

 

 

k ЭМ

 

 

 

 

 

 

 

 

 

 

 

k ЭМ

 

 

,

 

 

T

2

p 2 2T

 

d

 

 

 

 

 

 

 

 

p

2 a

 

 

 

 

 

 

 

ЭМ

ЭМ

p 1 a

M 2

M 1

p 1

 

 

 

 

 

 

ЭМ

 

 

 

 

 

 

 

 

 

 

 

 

(5)

 

 

 

T 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

M 2

 

;

 

a

M 1

2T

ЭМ

d

ЭМ

;

 

 

k

ЭМ

1.

 

 

 

 

 

 

ЭМ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

с

оптимальным

 

 

значением

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

 

коэффициента

затухания

dэм

= 0.7.

В этом

случае

время

регулирования tР

3 Tэм ,

а пере-

регулирование < 5%. Постоянную времени ЭМ выбираем из условия:

TЭМ tР / 3 .

(6)

С целью обеспечения требования астатизма ЭМ по отношению к задаю-

щему воздействию выбираем

kЭМ = kСАУ = 1.

Приведем структуру передаточной функции ЭМ (5) к виду (4), помножив ее числитель и знаменатель на множитель ( k1 / k0 p 1) :

 

 

 

 

 

(

k1

p 1 )

 

 

 

 

 

 

 

 

(

k1

p 1 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Wэм( p )

 

 

 

 

 

 

 

k0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k0

 

 

 

 

 

 

 

 

 

 

 

p2 a

 

 

 

 

k1

 

 

 

 

 

 

p2 a

 

 

 

 

 

 

k1

 

 

( a

 

 

 

 

 

p 1 )(

 

p 1 ) ( a

 

 

 

 

 

p 1 )(

p 1 )

M 2

M 1

k0

M 2

M 1

k0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

k1

p 1 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k0

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

aM 2k1

p3 (

aM 1k1

a

 

) p2

( a

 

 

k1

 

) p 1

 

 

 

 

 

 

M 2

M 1

 

 

 

 

 

 

 

 

 

k0

 

 

 

k0

 

 

 

 

 

 

 

k0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Сравнивая полученную передаточную функцию с передаточной функ-

цией (4), приравниваем коэффициенты знаменателей при одинаковых сте-

пенях р:

 

T1T2

 

aM 2k1

;

 

 

 

 

 

 

 

 

 

k0k

 

 

k0

 

( T1 T2 k1k2k )

 

aM 1k1

a

;

 

 

k0k

 

 

 

 

 

M 2

 

 

 

 

k0

 

7

[( k0k2 k1 )k 1]

a

 

 

k1

 

 

M 1

k .

k k

 

 

 

0

 

 

0

 

Отсюда последовательно определяем:

k1 T1T2 ; aM 2k

T1 T2 k1k2 k aM 1k1k aM 2 k0 k; k0 k2 k 1 aM 1k0 k;

k

 

 

T1 T2 aM 1k1k

 

aM 2

k

 

.

2

 

 

0

 

 

k1k

 

k1

 

 

 

 

 

 

 

Подставляя k2 в предпоследнее выражение, получим:

 

aM 2

k k

 

2 ( k

 

T1 T2 aM 1k1k

a

 

 

k )k

 

1

0;

 

 

 

 

 

M 1

0

 

 

k1

 

 

0

 

 

 

 

 

k1k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

M 2

k k 2

( T

T

2

)k

0

k

1

0.

 

 

 

 

 

 

 

 

 

 

0

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда получаем два решения:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( T

T )

 

( T

T )2

4a

 

k k

 

 

 

 

 

 

 

 

 

k0

 

 

1

 

 

2

 

 

 

 

1

 

 

2

 

M 2

 

1

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2aM 2k

 

 

 

 

 

 

 

 

 

 

 

 

12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M A P L E - П Р О Г Р А М М А Р А С Ч Е Т А П Е Р Е Д А Т О Ч Н Ы Х Ч И С Е Л Р Е Г У Л Я Т О Р А

Ввод параметров передаточной функции ОУ:

W(p)=k/[(T1p+)(T2p+1)]

> T1:=20.0; T2:=40.0; k:=10.0;

T1 := 20.0

T2 := 40.0

k := 10.0

Ввод заданного времени регулирования

> tр:=60.0;

 

:= 60.0

 

Определение параметров передаточной функции ЭМ:

 

Wэм(p)=1/(Am2*p*p+Am1*p+1),

где Tm = tР /3;

Am2=Tm*Tm;

Am1=2*Tm*d;

(7)

(8)

(9)

8

> Tm:=tр/3; d:=0.7;

Tm := 20.00000000

d := .7

> Am2:=Tm*Tm;

Am1:=2*Tm*d;

Am2 := 400.0000000

Am1 := 28.00000000

Определение передаточных чисел регулятора:

П е р в о е р е ш е н и е :

> k1:=T1*T2/Am2/k;

k1 := .2000000000

> k0:=((T1+T2)+sqrt((T1+T2)^2-4*Am2*k1*k))/(2*Am2*k);

k0 := .01000000000

> k2:=(-T1-T2+Am1*k1*k)/k1/k+Am2*k0/k1;

k2 := 18.00000000

Определение передаточных чисел регулятора:

Вт о р о е р е ш е н и е :

>k0:=((T1+T2)-sqrt((T1+T2)^2-4*Am2*k1*k))/(2*Am2*k);

k0 := .005000000000

> k2:=(-T1-T2+Am1*k1*k)/k1/k+Am2*k0/k1;

k2 := 8.000000000

9

h(t)

t, c

а) 1 и 2-й вар. регулятора: g(t) = 1(t), f(t) = 0, h( ) = 1,

tр = 56 с, = 4,3%.

h(t)

t, c

б) 1-й вар. регулятора: g(t) = 0, f(t) = 0,5, h( ) = 0, hmax = 0,85

10

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]