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

Андрущак С.В., аспирант, Беседин П.В., д-р техн. наук, проф., Алексеев С.В., канд. техн. наук, доц.

Белгородский государственный технологический университет им. В.Г. Шухова

РАЗРАБОТКА ЭКСТРЕМАЛЬНОЙ СИСТЕМЫ УПРАВЛЕНИЯ ПОДАЧЕЙ ШЛАМА

В ЦЕМЕНТНУЮ ПЕЧЬ

sg-bel@mail.ru

Рассмотренная методика экстремального управления подвижностью шлама при механическом воздействии позволяет эффективно снизить энергозатраты на получение клинкера за счет уменьшения количества воды, используемой для разжижения шлама на 5...9 % при воздействии вибрирующего органа в диапазоне частот от 10 до 35 Гц. При этом снижается расход тепла на испарение влаги в печном агрегате, а также уменьшается количество электроэнергии, потребляемой насосами при перекачивании шлама. Особенно важно то, что значительно увеличивается подвижность шламов, приготовленных на основе сырьевых материалов с высокой гигроскопичностью. Рассмотренную методику управления подвижностью можно применить не только для шлама, но и для других неньютоновских бингамовских жидкостей, таких как очень густые суспензии и пасты.

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

Введение. В цементной технологии мокрого способа производства цемента около 40 % топлива, расходуемого на обжиг клинкера, затрачивается на испарение воды, которая обеспечивает подвижность шлама [1]. Известно также, что снижение влажности шлама на 1% приводит к экономии топлива, расходуемого на обжиг клинкера, на 5 % [2, 3, 4], при этом снижение влажности шлама приводит к увеличению производительности печей.

На подвижность шламов оказывает заметное влияние воздействие вибрирующих органов [5, 6]. Особенно важно то, что вибрирующее воздействие значительно увеличивает подвижность шламов с высокой гигроскопичностью. Например, шлам с растекаемостью 52 мм в результате вибровоздействия увеличил свою подвижность более чем в 11 раз.

В результате экспериментальных исследований было установлено, что виброобработка позволяет снизить влажность транспортируемых шламов на 5.. .9 % при воздействии в диапазоне частот от 10 до 35 Гц, т.е. данная зависимость вязкости шлама имеет экстремальный характер [7]. При этом заранее неизвестно точное оптимальное значение частоты, так как цементные заводы имеют сырьевые материалы различной минералогической природы, а также это значение может «смещаться» на каждом конкретном заводе из-за изменений реологических свойств шламов. Для решения поставленной задачи был разработан способ регулирования подачи шлама в цементную печь [8], а также методология нечеткого вывода в задаче управления дозированием шлама [9, 10].

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

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

Структура объекта управления и его математическая модель Задача оптимального управления подвижностью шлама состоит в определении закона изменения во времени частоты вибрации, при котором перепад давления в вибрирующем патрубке, через который протекает шлам, имеет максимальное значение (рис. 1). При максимальном перепаде давления достигается наилучшая подвижность шлама для его перекачивания с помощью шламового насоса, а также существуют значения частоты вибрации, при которых достигается максимальное значение перепада давления в условиях дрейфа, вызванного изменением влажности шлама м!^). В общем случае уравнение безынерционного экстремального объекта регулирования имеет вид:

Ц = f , wl,...,wn ).

(1)

f(®*,...-Pm * Wl,...,Wn ) ^ f(P\,...-Pm, W1,...W„ ). (2)

Функция / обладает тем свойством, что для каждого фиксированного набора чисел С

Далее стоит отметить, что для любого набора

w¡,...,wn набор Ю]

единственный (функ-

Qw существует набор

,Ю„

при котором f

достигает минимума или максимума. Это означает в случае максимума, что

ция f имеет только одну точку экстремума-минимума).

Рис. 1 Экстремальная зависимость перепада давления от частоты вибрации и влажности сырьевого шлама

*

*

Безынерционность объекта позволяет упростить структуру адаптивного регулятора, что математически означает: rnk(t)=pk (k=1,m). Другими словами, управляющее воздействие формируется как настраиваемые параметры (из условия max f), т.е. управляющими параметрами.

Таким образом, уравнение экстремального объекта принимает вид

л = f (Pi,...,Pm, wi,...wn X

(3)

где (1=1,т) - управляющие настраиваемые параметры.

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

Уравнение однопараметрического объекта имеет вид:

v = f (A w).

(4)

В структуру математической модели объекта управления (рис. 2), созданной в среде Simulink, для исследования динамики экстремального управления подвижностью шлама входят следующие компоненты: электродвигатель (Engine), вибратор (Vibrator), датчик (Sensor) и экстремальная характеристика объекта (Extreme), входными сигналами для которой являются частота вибрации (frequency) и влажность шлама (humidity), а выходным - перепад давления на участке трубопровода с вибратором.

Vibfatcr

Рис. 2 Структурная модель объекта управления виброобработкой шлама в среде Matlab Simulink

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

характеристики объекта, поэтому рассматриваемую систему можно отнести к классу статических самонастраивающихся систем.

Существуют следующие методы поиска экстремума в однопараметрических экстремальных системах:

• измерение производной;

• запоминание экстремума;

• периодический поисковый сигнал.

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

Запишем приближенное выражение производной:

д/ ^ / (р + Ар, у) - / (1-АД у) др ~ 2Ар !

(5)

где Ав - достаточно малое число пробного шага.

При Ав —>0 выражение в правой части совпадает с производной /

Значения /(в+Ав, м) и /(в-Ав, м) доступны непосредственному измерению, поэтому приближенное значение производной / может быть получено при вычитании измерений для двух значений управляющего параметра (различающихся на величину 2Ав) и делении разности на число 2Ав. Так как воздействия на объект управления пробными шагами требует времени Т, то алгоритм поиска экстремума можно описать дискретно в виде:

р[(к + 1)Т] = р(кТ) - ак / 2Ар[/(р(кТ) + Ар, у) - /(р(кТ) - Ар, у)] (к = 0,1,2,...).

(6)

Представим параметр Ьк = ак/2Ав (к=1, 2, 3, ...), который назовем параметром рабочего шага.

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

1. в момент времени t0=0 к объекту прикладывается пробное воздействие в0+Ав, в0г Ав и определяются значения выхода объекта;

2. на основе выражения 6 формируется первый рабочий шаг в(Т)= в0гЬ1[[(в0+Ар, w) -/(в-Ав, w)];

3. к объекту прикладывается управляющее воздействие в(Т);

4. на второй итерации к объекту вновь прикладываются пробные воздействия в(Т)+Ав и в(Т)-Ав и измеряется выход объекта;

5. вычисляется значение управляющего параметра на втором шаге в(2Т)=в(Т)-Ь2[[(в(Т)+Ав, -/(в(Т)-Ав, w)];

6. к объекту прикладывается управляющее воздействие в(2Т).

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

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

В качестве алгоритма адаптации величины рабочего шага можно принять алгоритм

к+1

где АРв(кТ), АРв[(к-1)Т] - значения к-го и (к-1)-го рабочих шагов алгоритма (6) экстремального управления; V - некоторое положительное число, выбирая из условия устойчивости процесса

(7).

При сходимости решения к числу в* предположим, что пробный шаг Ав достаточно мал:

К+1 = К (1+V ■ А рр(кТ) а р(к - 1)Т]),

(7)

Наложим на функцию г (в, м) ограничения, чтобы эта функция удовлетворяла неравенству:

|г(Р, У) < Р-Р*

(10)

др

= г(Р(кТ), у) = гк

(8)

Тогда величина гк вычисляется достаточно точно. Следовательно, (7) примет вид:

р[(к + 1)Т] = р(кТ)-аЛ (к = 0,1,2,...). (9)

где 5 - неизвестное положительное число.

Определим эффективность к-го шага следующим числом:

р[(к +Щ]-р* (к = 0,1,2,...), (11) к р(кТ)-р*

модуль которого характеризует скорость приближения к экстремуму при к-м рабочем шаге. Если \ук\>1, то происходит удаление от точки в*, а если \ук\<1, то приближение к экстремуму. Подставляя данное условие в (9) получим:

к

Р(кТ) -р* —аЛ

"к ~-_ 1 — ак

к Р(кТ) —р* к

Из полученного выражения можно получить условие сходимости:

1 - a,.

PikT) -р-

< 1

или

0 < a.

PikT) -р*

< 2.

(13)

Так как дробь в выражении (13) всегда больше 0, то можно не рассматривать левую часть меньше 0 в двойном неравенстве.

На основе (10) получим условие сходимости поиска экстремума:

aks < 2 (k = 0,1,2,...).

(14)

Задавая величину 5 согласно условию, можно обеспечить устойчивость процесса поиска экстремума:

ак < -(к = 0,1,2,...). (15)

P(kT) -р*

(k = 0,1,2,...).

(12)

Из уравнений (49) и (50) получим выражения для расчета параметра 5:

riß, w) f (Р + ДР, w) - f (Р-ДР, w)

P-P*

2Дрр - р *

(16)

Основная часть

Разработка экстремального регулятора в

среде МайаЬ БтиНпк Рассмотрим реализацию описанного алгоритма экстремального управления в современном средстве имитационного моделирования Matlab Simulink, которое позволяет значительно упростить анализ динамики сложных экстремальных систем при различных алгоритмах работы экстремальных регуляторов [12]. В структуре модели системы экстремального управления (рис. 3) лежит модель объекта управления, созданная в среде Simulink, для исследования динамики экстремального управления подвижностью шлама [13].

Step ccnection

Рис. 3 Структурная модель системы экстремального управления в среде Matlab Simulink

Алгоритм управления заложен в блоке (Extreme Controller), который имеет один выходной параметр - частота вибрации, и семь входных параметров: Input - перепад давления на участке вибрации, является обратной связью; tau - постоянная времени объекта управления, включая привод вибратора и датчики; deltaB - значение пробного шага (5); a - величина рабочего шага (6); v - коэффициент устойчивости адаптации рабочего шага (7); Transient correction - включение отслеживания переходного процесса; Step correction - включение адаптации рабочего шага.

Экстремальный регулятор (Extreme Controller) воздействует через исполнительный механизм на объект управления. Объект

управления представлен в виде нелинейной экстремальной статической характеристики (Extreme), представленной на рисунке 1, которая при дрейфе является также функцией времени. Инерционные свойства объекта и измерительных устройств отражает динамическая часть в виде вибратора (Vibrator) с двигателем (Engine) и дифференциального манометра (Sensor). Дрейф влажности шлама задается блоком (Humidity) и монотонно изменяется на протяжении всего времени моделирования.

Далее рассчитаем входные параметры экстремального регулятора. Согласно выражению (5) число пробного шага Ар (deltaB) должно быть минимально, однако, реальный электропривод на базе частотного преобразо-

r

k

r

k

s

"k =

r

k

вателя сможет достаточно точно изменять обороты двигателя с шагом 2Ав = 120 об/мин, т.е. Ав = 1 Гц.

W (5) = Wд (я) ■ Wв (я) ■ Wм (я) =

где Т2=0,0165 с, 1=0,385 с.

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

Постоянную времени объекта управления (tau) для рассматриваемой модели можно определить аналитически из общей передаточной функции динамических звеньев: 54

(17)

0,01675в2 + 0,3855 +1

Далее рассчитаем величину рабочего шага согласно условию устойчивости процесса поиска экстремума (15). Для этого примем в=2 Гц, Ав=1 Гц, влажность сырьевой смеси ю=40%, а также прогнозируемое значение рабочего шага в точке экстремума в*=23.

Определим значения выхода объектов, про-

ведя моделирование: f (ß - Aß, w) = 1,669 Па; f (ß + Aß, w) = 2,298 Па.

По формуле (16) рассчитаем значение 5:

2,298 -1,669 ппле 1 ггпп

s = —-,—-—т = 0,015; — = 66,77.

2 • 1 • 2 - 23 s

В итоге, согласно выражению (15) установим величину рабочего шага ак равной 65.

По условию \ук\<1 примем коэффициент устойчивости адаптации рабочего шага V равным 0,1.

Представим алгоритм экстремального управления в виде структурной модели регулятора (рис. 4), разработанного в среде Simulink, и построенного в соответствии в рассмотренными алгоритмами управления [14].

Step ccnection

Рис. 4 Структурная схема экстремального регулятора в среде Matlab Simulink

Входной сигнал перепада давления, а также значение tau, поступают в блок анализа переходного процесса (Transient analysis). Данный блок выдает импульсные сигналы в моменты времени, когда переходный процесс пробного воздействия завершился, или через равные промежутки времени 4 х tau, если отключено отслеживание переходного процесса (Transient correction).

Согласно алгоритму поиска экстремума (6) к объекту управления должны прикладываться пробные воздействия [вк+А[в, [вк-А[в. Синхронизацию между всеми блоками регулятора для переключения между пробными воздействиями осуществляет блок (Control generator), ключевым компонентом которого является RS-триггер. Сигналом к переключению триггера является импульс, поступающий от блока (Transient analysis).

Следующий блок (Test memory) хранит в своей памяти значения выхода объекта для пробных воздействий P+Afi, рк-Ар. Сигналом к сохранению результата положительного fik+Afi воздействия является спад импульсного входного сигнала, а к сохранению отрицательного рк-Afi фронт импульсного сигнала. Выходным значением данного блока является разность между значениями выхода объекта для пробных воздействий P+Afi и pk-Ap одной итерации алгоритма.

В блоке (Initial step parameter) осуществляется расчет начального параметра рабочего шага b0= a0/2A^, который поступает на вход блока (Work step correction). Алгоритм адаптации параметра рабочего шага (7) заложен в блоке (Work step correction). Блок (Work step) согласно выражению (6) выполняет произведение параметра рабочего на разность значений выхода для

пробных воздействий ftk+Aft и fa-Aft. Вычисленный сигнал поступает на выход блока. Для хранения в памяти и подачи на выход экстремального регулятора пробных воздействий fik+tfi и pk-tp используется блок (B memory).

Переключение между пробными воздействиями осуществляется подачей на вход блока импульсов от блока (Control generator). При поступлении фронта импульса, происходит переключение от отрицательного пробного воздействия pk-tp к положительному в+kft. Когда на вход поступает спад импульсного сигнала, выполняется расчет k+1 рабочего шага согласно выражению (6), и сохраняется в блоке (Memory B). Для этого вычисляется разность текущего рабочего шага fik и рассчитанного в блоке (Work step) значения корректировки. Также при спаде импульсного сигнала осуществляется переключение к отрицательному воздействию pk-tp.

Анализ графиков переходных процессов при движении системы к экстремуму

Результаты исследования динамики функционирования экстремального регулятора представлены в виде графиков переходных процессов в двух режимах поиска: с отслеживанием завершения переходного процесса и адаптацией рабочего шага (рис. 6), и без них (рис. 5).

По полученным кривым графически определим параметры качества работы экстремальной системы, а именно частоту и амплитуду автоколебаний, а также величину потерь на поиск экстремума [15].

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

После выхода системы на экстремальное значение 24 Гц ее движение характеризуется наличием поисковых колебаний ("рысканий") около экстремума с амплитудой автоколебаний 1 Гц и частотой около 3 с"1. При отключенной коррекции выход системы на экстремальное значение в 7 Па произошел приблизительно за 7 секунд модельного времени, а с коррекцией за 4 секунды, следовательно потери на поиск составляют 7 и 4 секунды соответственно.

На графике смоделированного переходного процесса можно проследить заданный дрейф влажности от 41% до 40,3% в результате чего перепад давления на патрубке с вибратором упал от 7 до 6 Па. При этом также происходит смещение искомой экстремальной частоты в меньшую сторону.

¿30

°0 10 20 30 40 50 60 70

Время, с

Рис. 5 График переходного процесса движения системы к экстремуму без отслеживания завершения переходного процесса и адаптации рабочего шага в среде МайаЬ 81ши1тк

-

Г

F -

i i i i i i

О 10~ 20- ЗС 40 50 60 70

"о 10 20 30 40 50 60 70

Вреня. с

Рис. 6 График переходного процесса движения системы к экстремуму с отслеживанием завершения переходного процесса и адаптацией рабочего шага в среде МайаЬ 8ти1тк

Выводы

1. Механическое воздействие вибрирующего органа с частотой вибрации 10 -35 Гц позволяет увеличить подвижность цементного шлама, а за счет увеличения его подвижности - снизить влажность шлама на 5 -9 %. При этом, предложенная в работе математическая модель экстремальной системы управления подвижностью шлама, позволяет обеспечить оптимальную частоту вибрации вибрирующего органа, которая в свою очередь обеспечивает максимальную подвижность шлама.

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

например, неньютоновских и бингамовских жидкостей.

3. Разработанный алгоритм реализован в среде имитационного моделирования МаАаЬ Simulink, логика программирования в которой схожа с написанием программ для программируемых логических контроллеров.

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

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

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Беседин П.В., Трубаев П.А., Панова О.А., Гришко Б.М. Некоторые направления энергосбережения в технологии цемента // Цемент и его применение. 2011. № 2. С. 130134.

2. Беседин П.В., Андрущак С.В., Козлов

B.К. Методы физико-механического регулирования реологических свойств сырьевых шламов // Цемент и его применение. 2015. № 3.

C.70-73.

3. Пономарев И.Ф. Эффективные способы снижения влажности сырьевого шлама // Цемент. 1974. № 2. С. 9-11.

4. Богомолов Б.Н. Интенсификация процесса обжига путем предварительной виброобработки сырьевых шламов // Цемент. 1973. №7.С. 9-11.

5. Беседин П.В., Андрущак С.В., Козлов В.К. Исследование и моделирование процесса движения цементного шлама в лабораторных

условиях // Вестник БГТУ им. В.Г. Шухова. 2015. №4. С. 113-119.

6. Ткачев, В. В. Снижение влажности шлама путем разрушения его коагуляционной структуры. // Цемент. 1974. № 5. С. 12-13.

7. Касаткин А.Г. Основные процессы и аппараты химической технологии // 8-ое изд. -М.: Химиздат, 1971. 783 с.

8. Пат. 2014112882. Российская Федерация, МПК F 27 B 7/32. Способ регулирования подачи шлама в цементную печь / Беседин П.В., Ан-друщак С.В., Шарапов Р.Р., Козлов В.К., Новичков С.Г.; заявитель и патентообладатель Белгородский государственный университет; заявл. 02.04.2014; опубл. 2015.06.30.

9. Беседин П.В., Андрущак С.В., Козлов В.К. Методы нечеткой логики в задачах управления дозированием шлама // Цемент и его применение. 2014. № 6. С. 74-77.

iü.Besedin P.V., Andrushhak S.V., Kozlov V.K. Fuzzy Inference Technique in the Task of Sludge Batching Management // International Journal of Soft Computing. 2015. Vol.10-13.

11. Александров А.Г. Оптимальные и адаптивные системы: Учеб. пособие для вузов по спец. «Автоматика и упр. в техн. системах». М.: Высш. шк, 1989. 263 с.

12. Рубанов В.Г., Бушуев Д.А. Моделирование экстремальных систем управления в среде matlab и simulink как средство анализа динамики. // Научные ведомости Белгородского государственного университета. Серия: История. Политология. Экономика. Информатика. 2012. № 19-1 (138). . 169-175.

13.Кривилев А. Основы компьютерной математики с использованием системы MATLAB. М: Лекс-Книга, 2005. 492 с.

14.Черных, И.В. SIMULINK: среда создания инженерных приложений. М.: ДИАЛОГ-МИФИ, 2003. 496 с.

15.Казакевич В.В., Родов А.Б. Системы автоматической оптимизации. М.: Энергия, 1977. 288 с.

Andruschak S.V., Besedin P.V., Alekseev S.V.

DEVELOPMENT MANAGEMENT SYSTEM IN EXTREME FEED SLUDGE CEMENT KILN

Our procedure extreme mobility management cuttings by mechanical action to effectively reduce power consumption for clinker production by reducing the amount of water used for diluting the slurry to 5 ... 9% when subjected to a vibrating body in the frequency range from 10 to 35 Hz. This reduces the heat consumption for the evaporation of moisture in the furnace unit, as well as the amount of electric power is reduced, consumption of pumps for pumping sludge. It is particularly important that the mobility is greatly increased slurries prepared on the basis of raw materials with a high hygroscopic property. Consideration of mobility management technique can be applied not only to the sludge, but also for other non-Newtonian Bingham fluids such as very dense suspensions and pastes.

Key words: extreme control, the cement slurry, rheological properties of the slurry, mathematical model, the mechanical impact.

Андрущак Сергей Викторович, аспирант кафедры технологии стекла и керамики. Белгородский государственный технологический университет им. В.Г. Шухова. Адрес: Россия, 308012, Белгород, ул. Костюкова, д. 46. E-mail: sg-bel@mail.ru

Беседин Павел Васильевич, доктор технических наук, профессор кафедры технологии стекла и керамики. Белгородский государственный технологический университет им. В.Г. Шухова. Адрес: Россия, 308012, Белгород, ул. Костюкова, д. 46. E-mail: pvbesedin@yandex.ru

Алексеев Сергей Вячеславович, кандидат технических наук, доцент кафедры технологии стекла и керамики. Белгородский государственный технологический университет им. В.Г. Шухова. Адрес: Россия, 308012, Белгород, ул. Костюкова, д. 46. E-mail: Aleks_SB@List.ru