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

книги из ГПНТБ / Юсупбеков Н.Р. Автоматизация технологических процессов производства растительных масел

.pdf
Скачиваний:
41
Добавлен:
23.10.2023
Размер:
7.52 Mб
Скачать

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

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

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

Квозмущающим воздействиям относят масличность

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

сового цеха.

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

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

При отсутствии прибора, который бы непрерывно и автоматически измерял масличность шрота, концентра­ цию мисцеллы можно измерять специально созданными для этой цели автоматическими плотномерами (18; 19; 20]. В связи с этим целесообразно рассматривать кон­ центрацию мисцеллы как параметр автоматического контроля и регулирования масличности шрота с по­ мощью соответствующего изменения расхода раствори­ теля.

Успешную автоматизацию данного объекта ослож­ няет значительное транспортное запаздывание экстрак­ тора.

При регулировании таких объектов изодромными регуляторами с двумя параметрами настройки пере­ ходные процессы из-за значительного запаздывания в регулируемом объекте получают с большой динами­

120

ческой погрешностью [ 16J. Использование астатических регуляторов с тремя параметрами настройки повышает качество регулирования.

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

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

В описываемой системе автоматического

регулиро­

вания [2 1 ] один из датчиков концентрации

мисцеллы

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

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

12I

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

Для синтеза САР необходимо располагать динами­ ческими характеристиками объекта регулирования. Эксперименты проводили на Янгиюльском масложиро­ вом комбинате. Методика экспериментального опреде­ ления переходных функций объекта сводилась к сле­ дующему. На выходе экстрактора и в нижней части его загрузочной колонны установили весовые пневма­ тические плотномеры [22], разработанные в Ташкент­ ском политехническом институте. Вторичные приборы датчиков плотности регистрировали кривые разгона объекта. Одновременно с равными интервалами време­ ни в колбы отбирали пробы мисцеллы, концентрации которых при строго постоянной температуре определя­ ли стандартным пикнометрическим методом в централь­ ной заводской лаборатории.

Устанавливали постоянный расход жмыха, стабили­ зировали расход и температуру экстракционного бен­ зина. Затем наносили скачкообразное возмущение по расходу органического растворителя от его номиналь­ ного (10 ж3/час) до минимального значения \Ъм3)час). Таким же образом снимали динамические характерис­ тики экстрактора по интересующим нас каналам, на­ нося возмущения путем скачкообразного увеличения расхода от 5 м3/час до 10 м3/час.

На рис. 26 представлена кривая разгона объекта регулирования по каналу расход бензина — концентра­ ция мисцеллы на выходе экстрактора и расход бензи­ на — концентрация мисцеллы в нижней части загрузоч­ ной колонны. Сплошной линией показана запись вторичного прибора 1РЛ-29А, а точками проставлены результаты измерения концентрации мисцеллы стан­ дартным весовым методом.

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

Нормирование осуществляли путем отнесения теку­ щих значений концентрации мисцеллы к ее номиналь­ ной величине. В качестве номинальной принята

122

концентрация мисцеллы СНОм. = 15%. Величина номи­ нального расхода бензина GH0U. — 10 м3/час.

В первом и втором столбцах табл. 4 даны резуль­ таты эксперимента по определению кривой разгона объекта регулирования по каналу расход бензина — концентрация мисцеллы на выходе экстрактора. Для

с%

6 м3/час.

а мЧчас-

Рис. 26. Кривая разгона объекта по основному и промежуточному каналам.

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

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

5 яг31/час

-*вх.

=

10 м?!час .

 

ао=

*----

= 1.25

 

 

 

 

15%

 

О - " { § [1 -

 

(“

01 -

0,5 [1 -

(0)]j =

123

= 15 [(0,987 + 0,934 + 0,84010,734 + 0,620 +

+ 0,600)

-0,5]

= 63,225

Новый масштаб времени равен:

л =

Д =

- i i -

— 0,25.

v

t \

63,225

 

Для дальнейших вычислений составим вспомогатель­ ную табл. 5.

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а 5

6

1 — С „ „ „

 

1 - 0

d - W

+ x

1 _ 2 0 + Д

0 - " + т . ) х

 

 

 

 

 

 

 

X (1 -

0)

 

^ 2!

 

- С „ о Р М- )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

По основному каналу

 

 

 

0,00

 

1,00

 

 

1,00

1,00

1,000

 

1,000

0,25

 

0,987

 

 

0,75

0,74

0,531

 

0,520

0,50

 

0,934

 

 

0,50

0,46

0,120

 

0,110

0,75

 

0,840

 

 

0,25

0,21

-0,220

 

—0,180

1,00

 

0,734

 

 

0,00

0,00

—0,500

 

-0,360

1,25

 

0,620

 

—0,25

-0,155

-0,720

 

-0,440

1,05

 

0,600

 

-0,50

-0,300

-0,880

 

-0,520

 

 

 

По промежуточному каналу

 

 

 

0,00

 

1,000

 

 

1,00

1,000

1,00

 

1,000

0,33

 

0,973

 

 

0,67

0,650

0,29

 

0,280

0,66

 

0,900

 

 

0,34

0,300

-0,11

 

-0,099

0,99

 

0,830

 

 

0,01

0,083

—0,49

 

-0,400

1,32

 

0,800

 

—0,32

- 0,250

-0 ,8 0

 

-0,640

 

F.z = F* &Q j

V

[1 -

* вых

<гД0)] (1 -

Ш )

-

 

 

 

 

[

i= 0

 

 

 

 

 

 

 

— 0,5 [1 - *

ВЬ1Х- (0)]

J = 63,2252 • 0,25 (0,74 +

0,46 +

 

 

+ 0,21 -

0,155 -

0,3) — 0,5 -1 =

918,6

 

 

7-3 =

7-J АО

J

V

f1 _ *вых. (7А0)] [1 -

2Ш +

 

 

+

* ¥ " 1 -

 

° ’5 I1 -

- W (0)]]} - 63,225s-

 

•0,25 (0,52 + 0,11 -

0,18 -

0,36 -

0,44 —

 

 

 

-0 ,5 2 -0 ,5 ) = -77803

 

 

 

124

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

Таким образом, для исследуемого канала регулиро­ вания справедливо следующее дифференциальное урав­ нение:

918>6 <ых. + 63’225 *вых. +

= *вх .

которое в операторной форме выглядит так:

918,6 Do

,

63,225 г, ,

1

_

1

G(P)

1,25

+

1,25 И +

1

~

1,25 ‘

С (Р)'

На основании этого уравнения можно следующим образом записать передаточную функцию объекта ре­ гулирования по основному каналу:

^ofl, ^ ^ 734,9р2 + 50,58/>+1 в

'

( 173^

где -соб = 18 — время чистого запаздывания объекта регулирования.

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

 

 

 

 

Т а б л и ц а 6

/, мин

дс, %

снорм*

1 —Снорм-

0 А'.

0

0,00

0,000

1,000

0,00

15

0,40

0,027

0,973

0,33

30

1,50

0,100

0,900

0,66

45

2,60

0,170

0,839

0,99

60

3,00

0,200

0,800

1,32

На основании данных таблицы рассчитали коэффи циенты:

п _ 5 м(час (10 м3) час_п с

0

3% (15%)

~ 2,9

Fx = 15 (0,973 +

0,9 + 0,83 +

0,8 - 0,5) == 45.

125

Новый масштаб времени равен:

15

= 0,33

Р2 =

452-0,33 (0,65 +0,3 + 0,0083

-

0,25 - 0,5) = 133,65

F3 =

4 5 ' • 0,33 (0,28 - 0,09 -0 ,4 0 -

 

0,64 — 0,5) =* — 40095.

Итак, искомое дифференциальное уравнение объекта регулирования по промежуточному каналу имеет вид:

1 3 3 ’6 3 * ™ . Пр. + 4 5 < u x . n p . + 2 -5 -* вЫх.пр. = * в х . ( ^ - 8 ) ,

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

^об. г ( Р ) = 53,46/>2+ 18р+ 1 в ’ ( 174^

где хоб 2 = 8 мин. — время чистого запаздывания объек­

та по промежуточному каналу. На рис. 27 изображена структурная схема двухкон­

турной системы регулирования с вводом производной от промежуточной регулируемой величины. Для опре­

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

настройки

основного

регулятора

и дифференциатора

[2 2 ] систему на рис.

26 заменим

эквивалентной системой, показанной на рис. 28.

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

ного регулятора Рг обозначена:

 

WpAP) = Wp (P )W g (P).

(175)

В то же время передаточная

функция основного ре­

гулятора Р2 обозначена:

 

 

 

 

 

076)

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

через М70б экв (Pi.

Найти

ее можно из

следующей

системы уравнении:

Wo6, ( p ) x p(p)-,

 

•МР) =

 

* (Р) =

Wo6 г (р ) х р (р);

( 177)

■ *1 (Р) =

(Р).

 

 

126

г

+ Ч 1

 

I

X

Щ( Р}

>

- 1

-o—

 

 

Wg(P)

 

L

 

-J /

x,

х,

 

W0s.JP}

Рас. 27. Структурная схема двухконтурной системы регулирования с корректирующим и стабилизирующим регулятором.

Рис. 28. Двухконтурная схема регулирования

сдополнительным вводом воздействия от производной промежуточной величины.

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

медленно

и

практически

регулятор

Р,

успевает под­

держивать

вспомогательную величину х 1 почти точно

на заданном

значении,

иными словами,

х, ^ х .

Если исключить х,

)

и х (р) из

 

У,

системы уравне­

ний (177),

то можно записать:

 

 

об. экв

(р)

х(р)

 

Wq6J P ) xp (p )

 

Ху, (Р)

Wo6.AP)xp ip)

W0^ { P )

 

 

 

В рассматриваемой

системе регулирования исполь­

зуют реальное (инерционное) дифференцирующее зве­

но. Уравнение

движения этого звена

можно

записать

в виде:

 

 

dxn

 

7'

dxa

+ Хв

(178)

dt

= К'Т* dt

или в операторной

форме

 

 

 

Tg x B^ . ( Р) р +

^ых. ( Р)

= № , « .

Р•

(179)

где Kg — коэффициент усиления дифференцирующего звена.

Отсюда передаточная функция звена имеет вид:

(180)

Заменив в (180) р на (у'ш), получим выражение для амплитудно-фазовой характеристики звена:

W A ju) _ KgPgJm

_

V (<0^iy)a~bg> 1

oj ( | ------arctgoTg )

_

TgJw+

1

 

 

 

 

kg Tga

j arctg

шT„

(181)

\ П

' T

\

~ " 1

"

~‘s

V

(“ ^ )

2+ l

 

 

 

Приняв во внимание выражение (176), запишем:

тg р + i

(182)

 

128

В то же самое время передаточная функция ПИрегулятора определяется по формуле:

Wp (Р) = kp +

= T'V'i— p ~ kP

(183)

ИЗ.

1 и з . 1

 

Сопоставив уравнения (182) и (183), можем заклю­ чить:

ks

Т = Т

(184)

 

kp

 

Методика определения параметров настройки регу­ ляторов, входящих в состав рассматриваемой системы регулирования, таким образом, сводится к построению амплитудно-фазовых характеристик объекта регулиро­ вания по обоим каналам Wo6 (у'ш) и Wo6 (/ш). Выше

было показано, что

 

Wо б . ЭКВ.. (Р) = ^об„ (Р)

 

^об., (Р)

Следовательно,

чтобы построить амплитудно-фазовую

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

эквивалентного объекта W o6 экв з (уш),

необходимо при одинаковых частотах разделить моду­ ли характеристик, а фазы вычесть. Располагая АФХ эквивалентного объекта, определяют параметры наст­ ройки ПИ-регулятора Р 2 — коэффициент передачи и время изодрома. Затем по формулам (184) вычисляют параметры настройки ПД-регулятора (Kg и Tg). Для определения параметров настройки регулятора с пере­ даточной функцией Wp (P) пользуются тем же прие­ мом, что и для определения настроек вспомогательного регулятора Рг. Из рис. 27 следует, что передаточную функцию эквивалентного регулируемого объекта для регулятора с передаточной функцией Wp (Р) можно определить как передаточную функцию двух парал­ лельно включенных звеньев:

W*. э к в . (Р) = ^ 06., (Р)

(Р) +

(Р).

(185)

Таким образом, векторы амплитудно-фазовых ха­ рактеристик Wo6а (До), Wg (J(p) и Wo6 i (уш) складыва­

ются по правилу параллелограмма и получается АФХ W б экв (у'ш). Затем по последней характеристике обыч­

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

9—341

129

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