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

DOI: 10.12737/article_5a27cb8de48688.43489417

Ханин С.И., д-р техн. наук, проф., Старченко Д.Н., канд. техн. наук, доц., Мордовская О.С., канд. техн. наук, доц.,

Харин Н.П., инж.

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

ПРИМЕНЕНИЕ СИСТЕМ АВТОМАТИЗИРОВАННОГО ПРОЕКТИРОВАНИЯ ДЛЯ РАСЧЕТА ПРОЧНОСТНЫХ ПАРАМЕТРОВ КЛАССИФИЦИРУЮЩИХ УСТРОЙСТВ ШАРОВЫХ МЕЛЬНИЦ*

dh@intbel.ru

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

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

Введение. Шаровые мельницы распространены в различных отраслях промышленности для грубого и тонкого помола различных материалов. Эти мельницы имеют относительную простоту конструкции и эксплуатации. Их универсальность позволяет использование в условиях сухого и мокрого способа измельчения как одно-компонентных, так и многокомпонентных шихт и шламов, обеспечивая достаточно высокую однородность распределения частиц в измельченном продукте по их крупности и составу. Нерациональность процессов движения шароматери-альной загрузки, выделения из ее среды кондиционных частиц материала приводит к невысокой эффективности процесса помола, следствием чего является снижение производительности и повышение удельного расхода электроэнергии. Это ограничивает преимущества шаровых мельниц перед аналогичными помольными агрегатами [1-5]. Эволюционное совершенствование шаровых мельниц осуществляется на протяжении всей истории их существования. Одним из направлений их модернизации является изменение конструкций внутримельничных устройств, позволяющих осуществлять процесс выделения из шароматериальной загрузки частиц материала, достигших определенной крупности, и подавать на последующее доизмельчение или направлять в готовый продукт [6]. Это значительно снижает нерациональное переизмельчение мелких

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

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

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

Так, в работе [7] предлагается методика расчета сил, действующих на элементы цилиндрического классифицирующего устройства, основывающаяся на графоаналитическом методе [8].

Гк = ь У^ _ р3 ) АУк, (1)

где Ь - длина к-го сектора классифицирующего устройства, м; 7 - номер сектора контура загрузки (7 = 1,2,3), ушм - объемный вес шароматериальной загрузки, Н/м3; g - ускорение свободного падения, м/с2; рок и ры - радиусы, ограничивающие контур мелющей загрузки, м; Avк - угловой сектор, рад.

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

Методика. Заслуживает особого интереса дискретно-событийный подход к процессу движения мелющих тел [9-13]. В рамках этого подхода разработана математическая модель, адекватно описывающая протекающие в цилиндрической или конусообразной камере шаровой мельницы процессы пространственного движения мелющих тел, взаимодействия с бронефутеровкой и друг другом. В основу математической модели положено движение отдельно взятого сферического мелющего тела, которое последовательно взаимодействует с другими мелющими телами сферической формы, а также бронефутеровкой камеры и ограничивающих ее перегородок или днищ. Взаимодействия сферических мелющих тел с бронефутеровкой и друг с другом рассматриваются как мгновенные и вызывают изменения скоростей, которые изменяются также мгновенно. Мелющие тела на интервалах времени между соударениями движутся под действием силы тяжести по параболическим траекториям.

Время движения мелющего тела до броне-футерованной конусообразной поверхности камеры определяется из уравнения:

Ег -(К+№+к1 -у.' - *2в, ■ V;2 ] 2

(2)

+2+ у,'-У '+

' 008 О

'2 й ( Дт/

'008 в

= 0.

где хУ ; ' и V/, V у , V -координаты центра масс

мелющего тела в момент начала отсчета времени и проекции скорости на оси координат, соответственно, м; Дт - диаметр конусообразной камеры

у большего основания, м; вт - угол между образующей и продольной т -го конуса, рад; - радиус -го мелющего тела, м.

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

Введем следующие обозначения. /2 - момент инерции корпуса мельницы относительно оси 2 (рисунок 1); х, у - координаты точки взаимодействия; 5, 5 , 5 - проекции импульса взаимодействия на соответствующие оси; пх, п , П , Т, Т , Т, К, К' К - проекции единичных векторов локальной системы координат (рис. 2); ,5,5 -проекции импульса взаимодействия на соответствующие оси локальной системы координат; / -коэффициент трения скольжения; к - коэффициент восстановления при ударном взаимодействии; для . -го сферического тела: щ - масса, Ii - момент инерции, Уи, Уг, и - проекции скорости центра масс и угловой скорости на соответствующие оси. За локальную принята система координат с началом в точке взаимодействия и направлением осей: п - вдоль нормали взаимодействия, т -вдоль вектора относительной скорости, Ь - перпендикулярно осям п и т . При изменении положения места взаимодействия происходит изменение проекций единичных векторов (пх,пу,пг,тх,ту, тг, Ьх,Ьу, Ьг) локальной системы

координат и координат точки взаимодействия (х1 и У,)-

Обозначим индексы параметров следующими цифрами: 0 и 2 - до и после взаимодействия; 1 - в момент, когда сила трения обращается в ноль.

Импульсы взаимодействия мелющего тела с бронефутеровкой камеры мельницы и эффективный коэффициент трения скольжения V определяются следующими выражениями:

5 =-

-(1 + к)(Уоп -Ар (■ пу -■ пх))

1 (х,пу -УП )2 (Х,Ту -у,Тх )

-т--V-т-

(3)

Т] = Ш1П

(Х, • пу - У, ■ пх )

,(4)

1 (х п - у п ) (X п -у п )

-+ у ^^-(х ■_ _у ■т) (5)

7 (х,Ту-у$тХ)

2 щ I,

I х Т - у тх

(х ■ п - у ■ п

где

V, = Уот + г.Щрь - АрОТуХ, - ту);У = (1 + к)(Урп - Ар(х,Иу -уп)) •

т

+

Рис. 1. Бронефутерованная камера в абсолютной системе координат

Рис. 2. Мелющее тело в локальной системе координат

Нагрузка от воздействия мелющих тел с силой F на фрагмент поверхности бронефутеровки камеры площадью Q вычисляется как д' = F / Q . Касательная ^ и нормальная ^ составляющие силы F рассчитываются как:

/Л?; ^ ^/Л?. (6)

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

Бронефутерованные поверхности цилиндрической или конусообразной (рис. 3) камер разделяются на N участков в направлении оси вращения 2 и на N секторов в поперечных сечениях. Определим номер участка на бронефутерованной поверхности конусообразной камеры, взаимодействующего с мелющим телом в точке, имеющей координаты (х8, у, 23). Конусообразная ка-

мера характеризуется: углом вт наклона образующей, длиной , диаметром Дт у большего основания.

Длина Л1 участка (рис. 3), при их количестве N, равна:

А/ = т / N.

Рис. 3. Схема к разделению конусообразной камеры на участки

Его номер (в направлении оси Z ):

(8)

Угол между ограничивающими сектор радиусами:

" ' (9)

Аа = 2^ Номер этого сектора:

N„ •

N „ =

а

'Аа'

(10)

где а8 - угол между радиус-вектором точки взаимодействия и осью ОХ .

as = arctg—.

(11)

Развёртка выделяемого фрагмента поверхности бронефутеровки камеры представляет форму трапеции. Её высота вычисляется как h = А/ / cos вт, а размеры оснований - нижнего И1 и верхнего И2:

И, =Аа(°1у2-Nsl -А/-tg0mj (122)

И2 -N + 1)-А/- tgdm j. (13)

Площадь выделяемого фрагмента поверхности в форме трапеции:

х

s

Q = h

H + и2

AaAl

Dm/ _ /2

(N,t + 0.5)-A/ • tg0.

(14)

2 2 cos 6m

Основная часть. Для компьютерной реализации процессов движения мелющих тел в камере цилиндрического классифицирующего устройства, их взаимодействия с бронефутеров-кой и друг с другом формируются начальные положения мелющей загрузки [14]. Затем для мелющих тел осуществляется расчёт процессов их движения, взаимодействия с бронефутеровкой и друг с другом. На бронефутеровке устройства устанавливается номер фрагмента и координаты точки взаимодействия с мелющим устройством. Импульсы взаимодействия мелющих тел с броне-футеровкой суммируются по соответствующим фрагментам за назначенный промежуток времени и рассчитываются действующие на фрагменты нагрузки.

Для шаровой мельницы D*L=2*10,5 м осуществлен расчёт напряжённо-деформированного состояния цилиндрического классифицирующего устройства. Была принята рабочая относительная частота вращения барабана мельницы у

= 0,76 укр. Устройство длиной L = 0,9 м, с внутренним диаметром Dцку = 1,8м выполнено из стальных стержней диаметром dc = 0,06 м, расположенных с промежутками и образованием щелевых отверстий шириной ащо = 0,003 м и длиной 1о = 0,2 м. Материалом стержней является сталь 09Г2С ГОСТ 2590-206 ([ат ] = 275 МПа). Устройство загружено мелющими телами с коэффициентом ф = 0,3. Так как оно длительное время эксплуатируется в условиях высокого абразивного износа, то рассматриваем состояние при 50 % износе стержней по их диаметру. Воспользовавшись, в соответствии с [14], компьютерной реализацией процессов движения мелющих тел в камере устройства, их взаимодействия с бронефу-теровкой и друг с другом рассчитали нагрузки на фрагменты конструкции устройства. Полученные результаты, приведенные в таблице 1, применили для расчёта напряжённо-деформированного состояния цилиндрического классифицирующего устройства в CAD/CAM/CAE системе UmgrapЫcx NX [14]. Результаты расчета напряжений и деформаций в конструкции устройства приведены на рисунках 4-6.

Таблица 1

Нагрузки на цилиндрическое классифицирующее устройство мельницы D*L = 2*10,5 м

Нагрузка на фрагменты устройства Р, Н

Угловой номер фрагмента

V чс с- сх о с - с с V чс с- сх о с 1 с с 2 с с 2 yf 2 чс 2 с- 2 сх 2 о 2 с 2 с с 3 с 3 с yf 3 чс 3 с- 3 сх 3 с 3 с 3 5

Угол, град.

с с сх 2 с 2 с- 3 с 3 чс 4 4 V 5 с чс V 6 о 6 6 сх 7 с сх сх сх с о t" о с с чс с - yf с е- V с о с с сх с с сх с 4 с- 5 с чс чс чс F V t" с сх yf сх а сх

5 С с 55 С 1 с 0< чс с- 6 с 2 с с- 2 о чс 4 а сх с- 1 с чс 6 1 с сх с чс V 1 2 с V с сх с 1 3 с V 2 3 ЧС 4 3 F 3 3 1 с-с 3 сх сх с 3 V чс 3 сх 6 2 сх а 3 с с 2 3 о сх 3 сх чс 2 с с 6 2 3 32 12 F 4 yf 7 7 о 4 чс о о с-чс t" с с 7 t" с 5 с с с с

2 с V с с- чс с сх 1 с F сх с 17 V чс 6 с с УГ 1 о о 5 2 с- 3 3 V о 4 5 сх 4 4 с с 5 4 сх 5 с- 1 5 с чс 5 4 2 V 5 5 7 УГ 5 сх 5 5 с чс с- 5 сх 7 4 t" сх 6 ^J 7 с 5 4 с с г 4 сх с 4 с 3 сх 3 3 чс с- 2 с yf 72 0< V чс ^J с 4 сх с 21 V yf 7 с с с а е-

U С с с с с с 3 V чс с с чс сх с 1 с чс с чс сх с 6 чс 3 с сх ЧС 1 чс с- 6 2 чс с- 4 3 сх V 3 53 1 с с 4 сх чс с 4 с с- 3 4 с 3 4 с сх 7 с а 1 4 сх сх 4 с 4 4 с с 5 4 с с 6 4 о yf 54 о о 1 4 V чс 5 3 5 с а 2 ^ чс 2 с с 6 2 чс о V с V 5 с 5 сх о 1 с с 7 36 V о 13 с с

5 5 1 с а с с 4 V о 1 t- 1 сх с- с- 1 сх сх 7 1 с 3 2 1 сх 6 2 с с- V 2 с 53 о с 1 4 о чс 7 4 чс V 5 4 с УГ о 4 4 с чс 5 4 чс с 7 4 6 yf 4 4 4 t" с- 4 с-г- 3 3 с с- 4 3 с с 4 3 с- чс 43 7 а 1 3 V с- 6 2 с 42 о сх 1 F чс с 4 t" с 21 с F w чс 7 с с с

с 5 а с Е- 5 1 и с С < % 5 с с V V о сх V 1 ст сх с с с V 4 с V 1 с-сх 5 2 с-сх 6 3 с 7 3 с-с- чс о 4 6 а сх 4 t" чс 7 V 1 4 чс с- 2 4 V чс t" 4 с о 4 4 сх с 64 с о 1 4 с с 6 3 с 7 3 сх 63 с с 3 о чс 4 с чс 2 4 с с- 5 о 7 о с с чс 11 с с о с- 4 V с 2 о с 2 с

6 с с с V с- 71 с с 4 V с УГ о о 5 со 7 1 с 6 2 1 чс с 3 5 V 2 3 5 с 6 4 V 7 4 с с- 2 УГ чс чс 2 5 t" с- 2 5 чс сх 5 V сх 1 5 чс с 5 5 5 5 о с с- 4 с yf 7 4 с 4 о с г 5 с yf 2 4 с сх сх 3 с с-с 3 сх с- 2 3 с е-yf 3 t" сх 5 с сх с V yf 11 о чс 7 V yf 4 сх с 2 с

7 с с о с V сх 1 с- V 4 yf с 4 с- V с о 3 с сх 5 1 сх с с- 2 о 7 2 с 2 3 с ес с чс 1 3 F с с 1 3 о с 3 3 с 3 4 3 3 V с-с 3 5 с с 3 чс с- 2 с 6 с- чс 2 с с- 6 2 сх о 22 с о 1 с сх 4 а чс сх чс 6 ^J с чс 2 сх с с- с 7 V yf 5 с IT- 3 с чс чс с

Максимальные значения напряжений в цилиндрических стержнях - а = 255,3 Мпа (меньше допускаемого [ат ] на 7 %); перемещений -

8,3-10"4 м характеризуют рабочее состояние устройства после длительной эксплуатации.

Ustro«tvo_assyfemi_siml Solution 1 Результат Subcase - Static Loads 1. Статический шаг 1 Напряжение - По элементам/у злам. Без лееереджа. По Миэесу Мин ООО Макг. 398 62 Единицы = Н/ммл2(МПа) Деформация Перемещение ■ По узлам Величина

Ustrolstvo_aesyfem1_siml : Solution 1 Результат Subcase - Static Loads 1. Статический шаг 1 Перемещение - По узлам. Величина Мин О ООО. Макс О 831, Единицы = мм Деформация Перемещение - По узлам Велич*

Единицы = Н/ммл2(МПа)

Maximum

Элемент 1380314 Узел 1383174 255 304 Н/ммЛ2(МПа)

а) б)

Рис. 4. Распределение напряжений и деформаций в конструкции цилиндрического классифицирующего устройства мельницы БхЬ=2х10,5 м: а) напряжения, б) деформации

а) б)

Рис. 5. Распределение напряжений и деформаций в стержне цилиндрического классифицирующего устройства мельницы БхЬ=2х10,5 м: а) напряжения, б) деформации

а)

б)

Рис. 6. Распределение напряжений и деформаций в дистанционной вставке цилиндрического классифицирующего устройства мельницы БхЬ=2х10,5 м: а) напряжения, б) деформации

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

Выводы. Получена математическая модель, описывающая параметры процессов движения мелющего тела в конусообразной и цилиндрической камерах шаровой мельницы, взаимодействия с их бронефутеровкой и друг другом. Использование математической модели позволяет осуществлять компьютерную реализацию рассматриваемых процессов в камерах мельницы для всей мелющей загрузки в совокупности и устанавливать характеризующие их параметры. Разработана методика расчета нагрузок от воздействия мелющих тел на фрагменты поверхности бронефутеровки камер. С использованием разработанного программного обеспечения выполнен расчет нагрузок, действующих на цилиндрическое классифицирующее устройство промышленной шаровой мельницы. В программной среде CAD/CAM/CAE системы Unigraphicx NX рассчитаны напряжения и деформации в конструкции устройства.

*Работа выполнена в рамках Программы развития опорного университета на базе БГТУ им. В.Г. Шухова.

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

1. Романович А.А. Энергосберегающие агрегаты для измельчения материалов цементного производства с анизотропной структурой: Авто-реф. дис. д-ра. техн. наук: Белгород. 2014. 43 с

2. Науменко Ю.В. Основи теори режимiв роботи барабанних млишв. Рiвне.: СПД Зелент O.I., 2009. 282 с.

3. Иоргенсен. С.В. Помол цемента. Сравнительный анализ вертикальной валковой и шаровой мельниц // Цемент и его применение. 2006.

№ 6. С.39-44.

4. Севостьянов В.С. Энергосберегающие помольные агрегаты. Белгород. 2006. 435 с.

5. Хардер Й. Развитие одноэтапных процессов измельчения в цементной промышленности // Цемент. Известь. Гипс. 2006. № 1. С. 24-38.

6. Пат. 2236298 Российская федерация, МПК ВО2С 17/06. Трубная мельница / С.И. Ханин , В.В. Кайдаш, А.В. Чалов, Д.Н. Солодовников, О.С. Ха-нина; заявитель и патентообладатель Белгор. гос. технолог. акад. строит. материалов. -№2003113249/22; заявл. 05.05.03; опубл. 20.09.04, Бюл. №26. 6 с.

7. Солодовников Д.Н. Цилиндрическое внутримельничное классифицирующее устройство трубной мельницы: дис... канд. техн. наук. Белгород, 2009. 162 с.

8. Крюков Д.К. Футеровки шаровых мельниц. М.: Машиностроение, 1965. 184с.

9. Brilliantov N.V., Spahn F., Hertzsch JanMartin, Thorsten Poschel Model for collisions in granular gases // Phisical review E. 1996. Vol. 53. № 5. P.5382-5392.

10.R. Reichardt and W. Wiechert |/ Event driven simulation of a high energy ball mill. In Proceedings ASIM 2003, page 249. ASIM 2003.

11.Rosa Ramirez, Thorsten Poschel, Nikolai V. Brilliantov, Thomas Schwager // Physical review e. 1999. Vol. 60. № 4. P.4464-4472.

12.Schwager T., Poschel T. Coefficient of restitution for viscoelastic spheres: The effect of delayed recovery // Physical review e. 2008. Vol. 78. № 5. P. 1304-1316.

13. Thorsten Poschel Thomas Schwager // Computational granular dunamics. Models and algoritms. Springer - Verlang Berlin Heidelberg 2005. P. 322.

14. Свидетельство о государственной регистрации программы для ЭВМ №2012615855. Система расчёта и визуализации параметров движения мелющих тел в корпусе шаровой барабанной мельницы / Старченко Д.Н., Ханин С. И., Труха-чев С.С.; заявитель и патентообладатель БГТУ им. В.Г. Шухова; №2012613407, заявлено 27.04.2012; зарегистрировано. 27.06.2012.

15.Ельцов М.Ю., Козлов А.А., Седойкин А.В., Широкова Л.Ю. Проектирование в NX под управлением Teamcenter. 2010.783 с.

Информация об авторах

Ханин Сергей Иванович, доктор технических наук, профессор кафедры механического оборудования. E-mail: dh@intbel.ru

Белгородский государственный технологический университет им. В.Г. Шухова, Россия, 308012, Белгород, ул. Костюкова, д. 46.

Старченко Денис Николаевич, кандидат технических наук, начальник отдела подготовки электронных дидактических средств. E-mail: crovin@mail.ru

Белгородский государственный технологический университет им. В.Г. Шухова, Россия, 308012, Белгород, ул. Костюкова, д. 46.

Мордовская Ольга Сергеевна, кандидат технических наук, доцент кафедры механического оборудования. E-mail: unique.ox@gmail.com.

Белгородский государственный технологический университет им. В.Г. Шухова, Россия, 308012, Белгород, ул. Костюкова, д. 46.

Харин Никита Петрович, инженер кафедры механического оборудования. E-mail: bigpandanick@yandex.ru

Белгородский государственный технологический университет им. В.Г. Шухова, Россия, 308012, Белгород, ул. Костюкова, д. 46.

Поступила в октябре 2017 г.

© Ханин С.И., Старченко Д.Н., Мордовская О.С., Харин Н.П., 2017

Hanin S.I., Starchenko D.N., Mordovskaya O.S., Harin N.P. COMPUTER-AIDED DESIGN SYSTEM APPLICATION FOR CALCULATION PARAMETERS OF STRENGTH OF BALL MILL CLASSIFICATION DEVICES

Constructive improvement of the classification devices allows increasing capacity of ball mills, to reduce the specific energy consumption for the grinding process of the material. Developing methods of calculation that allows you to determine the magnitude and direction of applying loads to the design of the device is actual. The article describes the mathematical model of the process of the spatial motion of grinding media, their interaction with shell liners, and each other in the cone-shaped and cylindrical chambers of the mill. The method of calculating loads on classifying devices from the action of grinding media is considered based on a mathematical model. The results of calculation loads, strains and deformations acting on the cylindrical classifier of the industrial using mill in its design obtained using computer-aided design systems are presented.

Keywords: ball mill, classification device, grinding media, mathematical model of the process, methods, load, deformation.

Information about the authors Hanin Sergey Ivanovich, PhD, Professor. E-mail: dh@intbel.ru

Belgorod State Technological University named after V.G. Shukhov. Russia, 308012, Belgorod, st. Kostyukova, 46.

Starchenko Denis Nikolaevich, PhD, Assistant professor. E-mail: crovin@mail.ru

Belgorod State Technological University named after V.G. Shukhov. Russia, 308012, Belgorod, st. Kostyukova, 46.

Mordovskaya Olga Sergeevna, PhD, Assistant professor. E-mail: unique.ox@gmail.com.

Belgorod State Technological University named after V.G. Shukhov. Russia, 308012, Belgorod, st. Kostyukova, 46.

Harin Nikita Petrovich, Engineer. E-mail: bigpandanick@yandex.ru

Belgorod State Technological University named after V.G. Shukhov. Russia, 308012, Belgorod, st. Kostyukova, 46.

Received in October 2017

© Hanin S.I., Starchenko D.N., Mordovskaya O.S., Harin N.P., 2017