Оптимизация переходных процессов изменения уровня воды в барабане паровых котлов Текст научной статьи по специальности «Строительство и архитектура»

УДК 621.181.6

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

Докт. техн. наук, проф. КУЛАКОВ Г. Т.1), асп. КУХОРЕНКО А. Н.2)

'■'Белорусский национальный технический университет,

22Командно-инженерный институт МЧС Республики Беларусь

E-mail: [email protected]

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

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

Ключевые слова: переходный процесс, уровень воды, барабан паровых котлов.

Ил. 7. Библиогр.: 10 назв.

OPTIMIZATION OF TRANSIEN PROCESSES OF WATER LEVEL VARIATION IN DRUM OF STEAM BOILERS KULAKOV G. T., KUHARENKO A. N.

Belorussian National Technical University, "Institute of Command Engineers" of the Ministry for Emergency Situations of the Republic of Belarus (ICE)

The work of regulator in general three-impulse automatic control system of water level in drum of boiler doesn't supply quality of internal and external distur-

bance attack (presentation of regulation mistakes). That is why it is needed to improve. Different methods of proportional plus reset controller regulation of three-phase automatic feed control system are considered. There were suggested new methods to improve the quality of regulation of water level in boilers. Here the step system of automatic regulation was determined, on the base of transfer function.

It is noticed that optimal transient processes supply calculation of numerical value of transmission factor of regulator at у =2,618, it is more then was recommended, but statistic mistakes remain. The transient simulation method in fast-time scale is recommended, this allow to determine early the value of statistic mistake of regulation by disturbances of reheated steam consumption and properly change the task to compensating device of step automatic control system. And numerical value of time constant criteria should be calculated on the base of numbers of golden section(Phi), taking into account the definite time constant of lead section and time-lag, time-lag on controlled influence channel, and also taking into account maximum value of controlled influence. This method allow to reduce in two times the total time of regulation, to decrease absolute mistake of regulation in three times, and maximum value of regulation influence by feedwater in 1,7 times.

Keywords: transient process, water level, drum of steam boilers.

Fig. 7. Ref.: 10 titles.

Максимально допустимые отклонения уровня воды в барабане котла конкретного типа устанавливаются заводом-изготовителем. Поэтому поддержание заданного уровня воды в барабане котла является важным технологическим параметром, от которого также зависят надежность, безопасность и экономичность работы котельного оборудования электростанции. Это требование - важная и актуальная задача при их эксплуатации. Паровые барабанные котельные агрегаты, оставшиеся в наследство Республике Беларусь от распада СССР, получили широкое распространение в ее энергетике, как и трехимпульсные системы автоматического регулирования (САР), а также их модификации [1].

Регулятор трехимпульсной САР обеспечивает хорошую отработку скачка задания, но не обеспечивает качественную отработку внутреннего и крайнего внешнего возмущений (наличие статических ошибок регулирования). Такая работа регулятора снижет экономичность, надежность и долговечность работы котла, а также делает актуальными вопросы ее модернизации. Динамика уровня воды в барабане котла при возмущении расходом питательной воды характеризуется отсутствием самовыравнивания, а также явлением «набухания уровня» при возмущении расходом перегретого пара [2].

Простейший статический регулятор (П-регулятор) имеет угол опережения, равный ф = п, что позволяет ему обеспечивать минимальную динамическую погрешность регулируемого параметра в переходных режимах для объектов без самовыравнивания. Простейший идеальный астатический регулятор (И-регулятор) имеет малый угол опережения, равный ф = п/2, что не пригодно для работы на объектах без самовыравнивания, так как замкнутая система становится нейтральной или неустойчивой. Угол опережения у реального интегрирующего регулятора занимает промежуточное положение между этими двумя регуляторами, но динамическая погрешность регулирования у него больше, чем у статического. Поэтому в трехимпульсной САР уровня воды в барабане котла с целью повышения устойчивости регулирования применяют ПИ-регулятор с жест-

кой обратной связью. При этом САР состоит из двух контуров: внутреннего, включающего регулятор и устройство жесткой обратной связи, и внешнего, включающего объект регулирования и эквивалентный регулятор, состоящий из устройства жесткой обратной связи и регулятора [3].

Схема моделирования переходных процессов трехимпульсной САР питания барабанных паровых котлов [2] приведена на рис. 1.

Рис. 1. Схема моделирования САР питания барабанного парового котла: у(г) - основная регулируемая величина; хзд - заданное значение основной регулируемой

величины; хр(Г) - регулирующее воздействие; Л - внутреннее возмущение; /2 - возмущение нагрузкой (расходом перегретого пара); %р(р) - передаточная функция

ПИ-регулятора; %1(р) - передаточная функция объекта регулирования по каналу регулирующего воздействия; %в*(р) = W2(р) - %3(р) - передаточная функция внешнего возмущения расходом перегретого пара; ГОС - главная обратная связь; ВОС - внутренняя обратная связь

Передаточная функция объекта регулирования представлена в виде реального интегрирующего звена

%*( р) =-1-,

Т1 Р(т1 Р +1) (!)

где Т1 - постоянная времени идеального интегрирующего звена; т1 - постоянная времени инерционного звена первого порядка, численно равная величине запаздывания по каналу регулирующего воздействия; р - оператор Лапласа.

Передаточная функция ПИ-регулятора имеет вид

К (ТиР+1) Ти Р

% (Р) = , (2)

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

Численное значение коэффициента передачи ПИ-регулятора трехим-пульсной САР обычно рассчитывают с учетом передаточной функции объекта (1) по формуле

Т

К = у ^, т1

где у - доля отношения постоянных времени Т1/т1, численное значение которой у различных авторов переменное.

Например, в [2] рекомендуется принимать у = 1,0, в [4] - у = 0,8 при степени затухания переходных процессов у = 0,75, в [5] - у = 0,55 при степени затухания переходных процессов у = 0,9.

Передаточная функция внешнего возмущения расходом перегретого пара с явлением «набухания уровня» воды в барабане котла имеет следующий вид [6]:

Ш* (р) = W2 (р) - Ш (р) = - , (4)

Т3 Р + 1 Т4 Р

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

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

(р) Т п +1

ш (р) =_р =____(5)

Пжв1\Р) , т„ , > ГлтЛ ' ^ '

1 + шр (р) Т I 1+к л

К

р +1

т. е. получаем интегро-дифференцирующее звено.

Так как Т„

^1 + к Л кр

> Ти, передаточная функция эквивалентного регуля-

тора (5) в динамике представляет собой звено медленного реагирования.

На примере парогенератора БКЗ-210-140, работающего на общую паровую магистраль (максимальная производительность 210 т/ч; минимальная - 140 т/ч [2]), осуществим моделирование переходных процессов трех-имульсной САР питания при основных возмущениях с использованием пакета Simulink программного обеспечения МаЛаЬ. Динамика объекта регулирования по уровню при возмущении расходом питательной воды и перегретого пара парогенератора БКЗ-210-140 имеет следующие параметры: т1 = 21 с; Т1 = 37 с; к3 = 1,5; Т3 = 20 с; Т4 = 67 с.

При моделировании переходных процессов численное значение у в формуле (3) изменяли в пределах от 0,146 до 2,382 в соответствии с рядом чисел правила золотого сечения, приняв за целое отношение Т1/т1. Полученные при этом графики переходных процессов трехимпульсной САР питания при основных возмущениях приведены на рис. 2.

а

б

в г

Рис. 2. Переходные процессы трехимпульсной САР питания барабанного парового котла

БКЗ-210-140 при отработке основных возмущений с различными значениями у: а - отработка скачка задающего воздействия хзд; б - отработка внутреннего возмущения /1; в - отработка внешнего возмущения расходом пара с явлением «набухания уровня» /2; г - изменение регулирующего воздействия хр (расхода питательной воды) при отработке внешнего возмущения расходом пара; Г - время

Из приведенных на рис. 2 графиков переходных процессов видно, что при отработке внутреннего возмущения появляется статическая ошибка регулирования, которая устанавливается на одном уровне при любых значениях у. Вместе с тем, лучшие прямые показатели качества регулирования обеспечиваются при численном значении коэффициента у > 1: сокращается время регулирования, уменьшаются максимальная величина перерегулирования и отклонения расхода питательной воды. Статическая ошибка регулирования также возникает при отработке внешнего возмущения расходом пара, появление которого характеризуется явлением «набухания уровня» в барабане котла и стабилизируется на одном уровне. Следует отметить, что при топочном возмущении статическая ошибка регулирования отсутствует.

Аналогичное моделирование проводили с использованием ПИ-регулятора, настроенного по методу симметричного оптимума [7]. Однако при этом скачок задания отрабатывали с величиной перерегулирования 43,2 %, а максимальная динамическая ошибка регулирования при отработ-

ке внутреннего возмущения превысила допустимые значения, но в пределах допустимой степени затухания переходного процесса [2].

Так как динамика переходных процессов типовой каскадной и САР с дифференциатором при основных воздействиях адекватна [8], то дальнейшее исследование проводили с использованием каскадной двухконтурной САР (КСАР), стабилизирующий регулятор которой оптимизирует внутреннее возмущение, но при этом отрабатывает задающее воздействие хзд1 без перерегулирования, а корректирующий - внешнее возмущение (рис. 3).

.Л

.11

Л'зд]

■Тр(0

-у

"р| (Р)

"'„„ (р) Н'пн(р)

-1 ).зСС

-1

гос

к(р)

МО

Рис. 3. Схема моделирования каскадной двухконтурной САР уровня воды в барабане котла: у(Г) - основная регулируемая величина (уровень воды в барабане котла); хзд1, хзд2 - заданное значение промежуточной и основной регулируемой величины соответственно; хр(Г) - регулирующее воздействие (расход питательной воды); /ь/ - внутреннее и внешнее возмущения (расходом перегретого пара); ^р1(р), ^р2(р) - передаточная функция стабилизирующего и корректирующего регулятора; Шъ* - передаточная функция внешнего возмущения расходом перегретого пара; ГОС - главная обратная связь; ВОС - внутренняя обратная связь;

^оп(р), Wшl(р) - передаточная функция опережающего и инерционного участков объекта регулирования соответственно

Передаточная функция опережающего участка имеет вид

<п (р) = к

С р+1

(6)

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

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

^ (р) = -

1

^ (Р )

(7)

(р)1 - Wздl (р)

При этом оптимальная передаточная функция замкнутой САР внутреннего контура по задающему воздействию (критерий оптимальности)

Ьд1

= 1( ^) соответствует инерционному звену первого порядка

^ (Р) =

1

Тзд! Р + 1

(8)

где Тзд1 - заданная постоянная времени критерия оптимальности по промежуточной регулируемой величине.

Подставив (6) и (8) в (7), получим

Т* _1_ 1

Ур1 (Р) = . (9)

копТ зд1Р

Видно, что в качестве стабилизирующего регулятора получился ПИ-регулятор с одним параметром оптимальной динамической настройки Тзд1.

Выбор значений Тзд1 осуществляют в соответствии с правилом золотого

сечения, используемым для оптимизации КСАР в [9, 10]. При этом за целое принимали постоянную времени разгона по расходу питательной воды

(С = 5 с)

Тзд1 е [0,146; 0,236; 0,382; 0,44; 0,56; 0,618; 1,0]Го*п. (10)

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

У>2(Р) = ^--^д^, (11)

р У,„2 (р)1 - ^ (Р)

где передаточная функция эквивалентного объекта У,кв2 (р) с учетом заданной передаточной функции внутреннего контура (8) и передаточной функции инерционного участка (1) примет следующий вид:

1

Уэкв (Р) = ^ (Р)Щ (Р) = + -—. (12)

Т1 Р(Т1 Р + 1)(Тзд1 р + 1)

При этом заданная передаточная функция системы при отработке скачка задания хзд2 = 1(^ (критерий оптимальности) согласно структуре эквивалентного объекта будет равна

Узд2 (Р) =7-Цт (13)

(Тзд2 Р + 1)

где Тзд2 - заданная постоянная времени критерия оптимальности по основной регулируемой величине.

Подставив (12) и (13) в (11), получим значение оптимальной передаточной функции корректирующего устройства

ТЛ т, Р +1) (Т зд1 р +1)

Ур2(Р) = зд у (14)

3Тзд2

^ Р 2 + Тзд2 Р + 1

V У

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

времени: Тзд1 = 0,146Го*п; Тзд2 = 0,618т1; Тзд2 = 1т1. Графики переходных процессов КСАР при основных воздействиях приведены на рис. 4.

1,22

У, хзд 1,0

0,55

0,(6 0,44 0,22

1 г

I 1

/ /

1

-Г„,=0,14бС 1^ = 0,618т, — Я-=0,14И£; Г«, = 1т,

0,55

у,!г 0,41 0,22 0 -0,2 -0,11

-0,(6 -0,8 -1,0

100 200 300 1400 500 6, с 700 в

--■Т^О.ИбС; ¡(„пМ^йг, — ■-?'о1 = 0,14СТ;; Т^ - 1т,

Д\

|

|

|

к

\

V

1(6 -у, / -12 -10 -

. х10-3

0,(6

х» /2 0,14

0,2

0 -0,2

-0,14

---Г,.,=0,14(0;;,; У'73 , ----' = 0,146/^,; Т^ - 1т,

А

\ \ \\

[

\

\

100 200 300 400 500 г, с 700 г

-0,(6 -0,8

/

у

, 1 /

1 1 /

/ \ /

1 / V

„,= 0,146/;; 7 ^ = 61 Кг, ГИ1 = 0,146У^; Т^ - 1т,

V

0 100 200 300 400 500 г, с 700

0 100 200 300 400 500 г, с 700

б

а

Рис. 4. Переходные процессы каскадной двухконтурной САР уровня воды в барабане котла

при отработке основных возмущений с различными значениями времени Тзд1 и Тзд2: а - отработка скачка задающего воздействия хзд; б - отработка внутреннего возмущения /1; в - отработка внешнего возмущения расходом пара /2; г - изменения регулирующего воздействия хр при отработке внешнего возмущения расходом пара

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

Для устранения этого недостатка дополнительно измеряли расход перегретого пара (при моделировании крайнего внешнего возмущения /), устанавливали время стабилизации уровня г1, определяли величину статической ошибки регулирования, как разность заданного и установившегося значения ошибки регулируемой величины (уровень воды в барабане), и в момент времени стабилизации уровня изменяли задание корректирующему регулятору на величину Ау = хзд2 — у) с противоположным знаком.

В результате чего устраняли статическую ошибку регулирования при отработке внешнего возмущения с явлением «набухания уровня» (рис. 5), однако при этом увеличивалось время регулирования.

б

0,8 У, /2 0,44 0,22 0 -0,2 -0,44 -0,(3 -008 -1,0

а

—-т„,=0.1467;,,; з^о^кг, -Т'„, = 0,1467;; Т^ - 1т,

1,5

^ /2 1,0

0,5

-0,5

-1,0

V V

/

\г

V --Г„, = 0,146У;;; Г^ = Я5618Т1 — - 7„, = 0,14бС Т^ = 1Т|

0 100 200 300 400 500 г, с 700

100 200 300 400 500 г, с 700

Рис. 5. Переходные процессы каскадной двухконтурной САР уровня воды в барабане котла с устройством корректировки ошибки регулирования при возмущении расходом пара: а - отработка внешнего возмущения расходом пара /2; б - величина регулирующего воздействия хр при отработке внешнего возмущения расходом пара

Причем увеличение численного значения постоянной времени Тзд2 с 0,618т1 до 1,000т1 приводило к увеличению отклонений уровня воды в барабане котла от заданного значения (рис. 5а), однако несколько уменьшало отклонения расхода питательной воды (рис. 5б). Вместе с тем из-за того, что коррекцию задания проводили в моменты времени стабилизации уровня, т. е. после окончания переходных процессов (рис. 4), общая площадь ошибки регулирования по уровню была достаточно велика (рис. 5). Для устранения этого недостатка одновременно с изменением нагрузки (расхода пара) запускали программы моделирования САР уровня воды в барабане котла при возмущении расходом пара в ускоренном масштабе времени, например с уменьшением реального времени в 20 и 40 раз (рис. 6).

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

0

0,8г

У, Я -0,44 -0,2 -

0 --0,2 --0,44 --006 -

-0,8 --1,0 _

\ \ —т .1467 Г

\\

1 \ \

\\

\

-----

008г ^ /2 -

0,44 -0,22 -

0 -

-0Д --0Д --0^ --0^ --1,0 _

— Г„, =0.14йГ^,; Г^ = 1г,

Л

/ \\

1 \\

\\

--...... -----

0 1 2 3 4 5 6 7 8 9 г, с 11

0

1

2

3

4

5 г, с 7

Рис. 6. Отработка внешнего возмущения расходом пара / каскадной двухконтурной САР уровня воды в барабане котла в ускоренном масштабе времени: а - в 20 раз; б - в 40 раз

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

рис. 7 (для постоянных времени: Тзд1 = 0,146Го*п; Т 2 = 0,618т1).

008 У, / 0,44 0,22 0 -0,22 -0,44

-0,(5

1,5

^ /2 1,0

0,5

-0,5

0 50 100 150 200 250 г, с 350

0 50 100 150 200 250 г, с 350

Рис. 7. Отработка внешнего возмущения расходом пара / каскадной двухконтурной САР уровня воды в барабане котла: а - с корректировкой величины статической ошибки регулирования при г е 0; 5; 10; 15; 20 с и при г = 200 с; б - изменения регулирующих воздействий хр при отработке внешнего возмущения расходом пара

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

В Ы В О Д Ы

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

б

а

б

а

1

те численного значения коэффициента у обеспечивает отработку скачка задания без статической ошибки регулирования. Причем оптимальные переходные процессы обеспечивает расчет численных значений коэффициента передачи регулятора при у = 2,618, что больше чем рекомендуемые в [2, 4, 5]. Однако при этом сохраняются статические ошибки регулирования при отработке внутреннего и внешнего возмущений с явлением «набухания уровня».

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

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

Л И Т Е Р А Т У Р А

1. Г у р е в и ч, Э. З. Состояние, перспективы развития и технико-экономическая оценка автоматизации питания барабанных парогенераторов / Э. З. Гуревич. - Минск, 1974. - 60 с.

2. К л ю е в, А. С. Наладка систем автоматического регулирования барабанных паровых котлов / А. С. Клюев, А. Т. Лебедев, С. И. Новиков. - М.: Энергоатомиздат, 1985. - 280 с.

3. Н о в и к о в, С. И. Рекомендации по обучению приемам настройки регуляторов на специальном стенде / С. И. Новиков, Б. Я. Жихарев, М. Н. Г. Мац. - М.: СЩТИ ОРГРЭС, 1973. - 43 с.

4. С т е ф а н и, Е. П. Основы расчета настройки регуляторов теплоэнергетических процессов / Е. П. Стефани. - М.: Энергия, 1972. - 376 с.

5. С п р а в о ч н и к по наладке автоматических устройств контроля и регулирования / В. А. Дубровный [и др.]. - Киев: Наукова думка, 1981. - Ч. 2. - 940 с.

6. Ш и р о к и й, Д. К. Расчет параметров промышленных систем регулирования: справ. пособие / Д. К. Широкий, О. Д. Куриленко. - Киев: Техшка, 1972. - 232 с.

7. Ф р е р, Ф. Введение в электронную технику регулирования: пер. с нем. / Ф. Фрер, Ф. Орттенбургер. - М.: Энергия, 1973. - 193 с.

8. К у з ь м и ц к и й, И. Ф. Теория автоматического управления: учеб. / И. Ф. Кузьмиц-кий, Г. Т. Кулаков. - Минск: БГТУ, 2010. - 574 с.

9. К у л а к о в, Г. Т. Методика структурно-параметрической оптимизации каскадных систем автоматического регулирования на основе модифицированного упредителя Смита / Г. Т. Кулаков, А. Т. Кулаков, В. В. Кравченко // Энергетика... (Изв. высш. учеб. заведений и энерг. объединений СНГ). - 2012. - № 2. - С. 40-45.

10. С о р о к о, Э. М. Золотые сечения, процессы самоорганизации и эволюции систем: введение в общую теорию гармонизации систем / Э. М. Сороко. - М.: КомКнига, 2006. -264 с.

R E F E R E N C E S

1. G u r e v i c h, E. Z. (1974) State, Long-Term Plan of Development and Technical-Economical Analysis of Feed Automation of Drum Heat Generators. Minsk. 60 p. (in Russian).

2. K l u e v, A. S., L e b e d e v, A. G., & N o v i k o v, S. I. (1985) Ajustment of Automatic Regulation System for Drum Steam Boilers. Moscow, Energouzdat. 280 p. (in Russian).

3. N o v i k o v, S. I., Z h i h a r e v, B. Ya., & M a z, M. N. G. (1973) Training Recommendations on Regulators Adjustment on Special Bench. Moscow, SShchTI ORGRES. 43 p. (in Russian).

4. S t e f a n i, E. P. (1972) Basis of Adjustment Calculation for Regulators of Heat Power Processes. Moscow, Energy. 376 p. (in Russian).

5. D u b r o v n y i, V. A., Z a b o k r i t s k i i, E. I., T r e g u b, V. G., & K h o l o d o vs k i i, B. A. (1981) Ref. Book on Adjustment of Automatic Control and Regulation Devices. Kiev, Navukova Dumka, Part 1. 940 p. (in Russian).

6. S h i r o k i, D. K., & K u r i l e n k o, O. D. (1972) Parameters Calculation of Industrial Regulation Systems. Kiev, Texshka. 232 p. (in Russian)

7. F r e r, F., O r t t e n b u r g e r, F. (1973) Introduction in Electronic Devices of Regulation. Moscow, Energy. 193 p. (in Russian).

8. K u z m i t s k y i, I. F., & K u l a k o v, G. T. (2010) Theory of Automatic Control. Minsk: Belorussian State Technical University. 574 p. (in Russian).

9. K u l a k o v, G. T., K u l a k o v, A. T., K r a v c h e n k o, V. V. (2012) Methods of Structural Parameters Optimization of Step Automatic Regulation on the Base of Modified Predictor Smit. Izvestiia Vysshikh Uchebnykh Zavedenii i Energeticheskikh Ob 'edinenii - Energetika. [Proceedings of the Higher Education Institutions and Power Engineering Associations - Power Engineering], 2, 40-45 (in Russian).

10. S o r o k o, E. M. (2006) Gold Cross-Sections, Processes of Self-Organization, Evolution of Systems: Introduction Into General Theory of Harmonization System. Moscow, KomKniga. 264 p. (in Russian).

Представлена кафедрой ТЭС Поступила 07.04.2014