Разработка и методика определения параметров вибрационного стенда с консольной просеивающей поверхностью

Известия Уральского государственного горного университета. 2019. Вып. 2(54). С. УДК 622.277:621.927 https://doi.org/10.21440/2307-2091-2019-2-93-99

Разработка и методика определения параметров вибрационного стенда с консольной просеивающей поверхностью

Аркадий Васильевич ЮДИН, Магомед Казбекович АбДуЛКАРИМОВ, Анатолий Григорьевич ПОПОВ, Виктор Степанович ШЕСТАКОВ*

Уральский государственный горный университет, Россия, Екатеринбург

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

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

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

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

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

В ведение

На кафедре горных машин и комплексов УГГУ проводятся исследования по совершенствованию конструкций грохотов и определению их параметров для условий эксплуатации с трудногрохотимым материалом. Теоретические исследования опубликованы в [1-5], а для проведения экспериментальных исследований на кафедре разработан вибрационный стенд. Исследования по повышению эффективности грохочения проводятся и другими научными и производственными коллективами [6-14].

Решение задач исследования

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

Стенд предназначен для экспериментальных исследований вибрационного грохота с каскадной просеивающей поверхностью (ПП), колосники которой закреплены консольно. Назначение грохота - разделение трудногрохотимой горной массы (ГМ), например, содержащей глинистые включения, на фракции ±0,04 м с наибольшим размером загружаемого куска 0,1 м. На кафедре разрабатывается усовершенствованная конструкция на базе серийного грохота ГИТ-51, эффективность грохочения и будет исследована на стенде. Общий вид стенда показан на рис. 1.

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

Направленные колебания РО осуществляются двумя самобалансными вибраторами, связанными упругими муфтами общим валом. Вибраторы установлены на общей балке и закреплены с обеих сторон РО и приводятся во вращение электродвигателем П-42 постоянного тока с регулируемой частотой. Двигатель и РО связаны через упругую лепестковую муфту. Вибраторы имеют два дебалансных вала, установленных на сферических подшипниках в сварном корпусе, на концах расположены дебалансы. Валы смещены относительно друг друга при помощи зубчатой передачи с передаточным отношением, равным 1, что обеспечивает вращение их в противоположном направлении. Упругая система включает четыре опоры, в которых находятся восемь винтовых пружин. Стенд включает пульт управления, измерительную аппаратуру для пуска и регулирования частоты колебаний РО. Просеивающая поверхность закреплена внутри бортовин и представляет собой каскадную колосниковую решетку.

Движение РО на первом этапе исследований рассматривается как направленные колебания одномассной системы с зарезонансным режимом работы и с углом вибрации 30о. Параметры РО подобраны так, что коэффициент режима

"Н shestakov.v.s@mail.ru

ИИрэУ/о rcid.org/0000-0001 -9468-6696

Рисунок 1. Общий вид стенда грохота. Figure 1. General view of the screen stand.

Г соответствует движению ГМ по ПП подбрасыванием. Частица ГМ будет иметь контакт с поверхностью до тех пор, пока составляющая силы инерции частицы, нормальная к поверхности, не превысит нормальную составляющую силы тяжести

2

MAMw sin в • sin wt > G4sin а,

где М - масса РО; Оч - сила тяжести частицы ГМ; АМш2 - амплитуда ускорения РО; в - угол вибрации; а - угол наклона РО к горизонту; w - круговая частота вынужденных колебаний.

Поскольку в условии (1) левая часть уравнения имеет максимальное значение при wt = п/2, то отрыв частицы наступит тогда, когда коэффициент режима Г будет больше 1:

Л„ш2 sin 6 , ч

Г = -М----> 1.

g cos а

Из опыта работы грохотов известно, что необходимое по интенсивности встряхивание ГМ происходит, если угол вибрации в принят не менее 30о, а коэффициент режима Г в диапазоне 2,7-3,0. По выбранным параметрам Ам, Г, а, в и равенству (2), решенному относительно ш, определяется требуемая частота круговых колебаний РО:

1

rg cos а

Лм sin в

Движение одномассной колебательной системы с инерционным приводом описывается приведенным дифференциальным уравнением

цСг Сг то 2

у +-г— у +-г— у =-0— Гш ,

М + то М + то М + то

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

Из (3) вычисляется амплитуда вынужденных колебаний грохота в установившемся режиме по формуле:

2

т гш

А„ =- 0

M

T¡Cr - (M + то)ш2 + |i2Cr2w2

Рисунок 2. Динамическая расчетная схема грохота. Figure 2. Dynamic design diagram of the screen.

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

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

Р = .

-С-

M + т

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

Опыт эксплуатации комплекса очистки известняка от глины на Чаньвинском карьере (Пермский край) показал, что колосниковые ПП виброгрохотов при разделении ГМ на фракции ±0,04 м часто забиваются глинистыми фракциями, при этом эффективность грохочения снижается [1]. Колосниковые ПП серийных грохотов осуществляют колебания совместно с РО, расстояние между колосниками всегда постоянно и не способствует интенсивному встряхиванию ГМ.

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

Просеивающая поверхность стенда включает четыре каскада, закрепленных между бортовинами. Каждый каскад включает опорный блок, на котором посредством болтовых соединений закреплены консольные колосники. Предусмотрено использование отличающихся колосников по массе, жесткости, длине, а также установка на них грузов. Колосники набраны из пружинной листовой стали 65Г толщиной 1-3 мм.

Вторая фаза разработки стенда направлена на реализацию сформулированной идеи, расчетная схема стенда приведена на рис. 2.

Без учета неупругих сопротивлений в системе уравнения движения для масс М и тпр принимают вид:

M-M = -С Ум + С(Ут - Ум ) + P cos wt]

К

т y = -С

пр' т к4' т У М ' J

где yM,yM - перемещение и ускорение массы М; ут,ут - перемещение и ускорение массы т; тп - приведенная масса колосников, учитывающая массу дополнительного груза G и массу единицы длины колосника, С - жесткость колосника; P cos wt - возмущающая сила привода.

Вводя обозначение (Сг + Ск)/М = а; CJM = b; Ск/тпр = с; P/M = d, уравнение (4) можно записать в компактной форме:

у М + аУм - ЬУт = d C0S Wt 1

- =-С+ с Г

У т С—т С-М

Частное решение системы (5) имеет вид:

= A.sin wt, у = A sin wt.

УМ M У т т

Таблица 1. Исходные данные параметров стенда для расчета.

Table 1. Initial data of stand parameters for calculation.

Параметр Значение параметра

Размер РО в плане L х B, м 1,2 х 3

Угол наклона РО а, град. 0

Угол вибрации р, град. 30

Частота вынужденных колебаний ш, кол/мин 0-1500

Амплитуда колебаний РО АМ, мм 3,0-4,0

Общий коэффициент жесткости упругих опор РО С, Н/м 533100

Возмущающая сила привода Р, кН 6,5-8,5

Количество каскадов ПП N, шт. 4

Количество колосников (общее) W, шт. 32

Масса рабочего органа M, кг 300-350

Частота собственных колебаний РО шРО, 1/с 35-40

Длина консоли колосников l, м 0,25-0,3

Толщина колосника h, м 0,002; 0,003; 0,004

Вес дополнительного груза на колосниках Gг, Н_5; 10; 15_

Таблица 2. Результаты численного моделирования колебаний двухмассной системы. Table 2. Results of numerical simulation of oscillations of a dual-mass system.

Параметры стенда Расчетные параметры

l h G г m C к ш г ш к АМ А m

0,25 0,002 0,15 1,5 3276 354 476 -3,9 1,8

0,25 0,003 0,15 1,7 11 059 362 793 -15,5 187,9

0,25 0,004 0,15 2,0 26 214 362 1141 -3,3 -7,0

0,3 0,002 0,15 1,6 1896 314 395 -3,9 0,8

0,3 0,003 0,15 1,9 6400 357 586 -4,1 4,0

0,3 0,004 0,15 2,1 15 170 359 833 -0,3 -44,5

Ам /Аст

Am /Аст

16 11 6

0 -4

-9 -14 -19

250

750

16

11

6

1 0

-4 -9 -14 -19

250

Рисунок 3. Амплитудно-частотная характеристика. а - рабочего органа; б - просеивающей поверхности; Аст - статическое перемещение просеивающей поверхности.

Figure 3. Amplitude frequency response. a - work tool; b - deck plate; Ast - static motion of the deck plate.

После подстановки у1 ,у М,ут,у т в уравнение (5) получим систему для определения амплитуд колебаний Ам и Ат

(а - ш)Лм - ЬАм = й 1 сАм + (с - ш2)Лт = 0 Г

где АМ, Ат - амплитуда колебаний просеивающей поверхности грохота и колосника.

Совместное решение уравнений (6) дает следующее значение амплитуд колебаний масс М и тпр:

Ам = РС / т ) - ш2)/(МВ); Ат = РСк/(МтщБ),

а

б

со

где

В = {[

Выражение (8), входящее в знаменатель выражения (7), представляет собой квадратичное уравнение относительно ш и определяет частотное уравнение системы (6). Решая частотное уравнение относительно ш2, определяют две собственных частоты всей системы ш и ш :

г к

ш2г к = 0,5[

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

г к.

Отношение амплитуд АМ и Аш вынужденных колебаний определяется из выражения

А 2

-м = [

ш

В выражениях (7)-(10) значения параметров определены по формулам:

Рг = [(С/(М + шо)] °'5; рк = [(Ск/шпр)]°5 ; шпр = ш + цш!;

С = Щ /13;С = р2М,

к ' X ' Г г Г '

где рг, рк - частота собственных колебаний массы М и колосников; тпр - приведенная масса колосниковой ПП; ц - коэффициент приведения массы колосника, ц = 33/140 [15]; l - длина консоли колосника; Е - модуль упругости; Jx - момент инерции поперечного сечения колосника.

Разработке стенда предшествовали исследования поведения грохота как двухмассной колебательной системы на примере серийного вибрационного грохота типа ГИТ с измененной просеивающей поверхностью. Анализ и численное моделирование показали, что массы M и m выполняют разноамплитудные колебания AM и Am и отношение амплитуд AM/Am зависит в основном от параметров колосников и высшей частоты резонансных колебаний шк. Результаты анализа приведены в работе [1]. На следующем исследовании поставлена задача - проверить результаты экспериментально.

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

Для моделирования параметров колебаний стенда в соответствии с предложенной методикой разработана компьютерная программа на алгоритмическом языке (фонд кафедры ГМК УГГУ). Исходные данные для расчета параметров стенда сведены в табл. 1.

В результате численного моделирования получены значения параметров экспериментального стенда. Фрагмент расчетных значений приведен в табл. 2. На рис. 3 показаны расположения резонансных кривых.

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

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

ЛИТЕРАТУРА

1. Юдин А. В., Шестаков В. С., Абдулкаримов М. К. Численный анализ колебаний двухмассовой консервативной системы вибрационного грохота для трудногрохотимой горной массы // Изв. вузов. Горный журнал. 2018. № 5. С. 93-101.

2. Юдин А. В. Моделирование вынужденных колебаний просеивающей поверхности инерционного грохота с консольно защемленными колосниками // Изв. вузов. Горный журнал. 2016. № 6. С. 63-70.

3. Stadkowski A., Yudin A., Komissarov A., Lagunova Yu., Akhmetova M., Stolpovskikh I. Calculation of parameters and design of the movable transfer station with vibrating screen feeder for the conveyor of deer queries // International Journal of Engineering and Technology (UAE). 2018. Vol. 7. special issue, № 2. P. 148-151. http://dx.doi.org/10.14419/ijet.v7i2.23.11904

4. Волков Е. Б., Ляпцев С. А. Влияние угла наклона рабочей поверхности вибрационного грохота на эффективность грохочения // Современные проблемы науки и образования. 2013. № 4. С. 8. URL: http://www.science-education.ru/ru/article/view?id=9642

5. Волков Е. Б., Ляпцев С. А. Компьютерное моделирование процесса грохочения // Международный журнал экспериментального образования. 2012. № 4. С. 49-50.

6. Вайсберг Л. А. Проектирование и расчет вибрационных грохотов. М.: Недра, 1986. 144 с.

7. Вайсберг Л. А., Коровников А. Н., Трофимов В. А. Модернизация технологических циклов грохочения на основе инновационного оборудования (к 100-летию института «Механобр») // Горный журнал. 2017. № 1. С. 11-17. https://doi.org/10.17580/gzh.2017.01.02

8. Картавый А. Н. Вибрационные агрегаты для переработки минерального и техногенного сырья. Моделирование и элементы расчета по критериям энерго- и ресурсоэффективности. М.: мГгУ, 2013. 328 с.

9. Газалеева Г. И., Цыпин Е. Ф., Червяков С. А. Рудоподготовка, дробление, грохочение, обогащение. Екатеринбург, 2014. 914 с.

10. Назаров К. С., Фет Ш. Анализ современных конструктивных решений, повышающих эффективность виброклассификации трудногро-хотимых материалов // ГИАБ. 2009. Т. 16, № 12. С. 383-393

11. Шишкин Е. А., Лебедев А. И. Исследование параметров вибрационного грохота с применением инструментов имитационного моделирования // Ученые заметки ТОГУ. 2016. Т. 7. С. 281-286. URL: http://ejoumal.pnu.edu.ru/media/ejoumal/articles-2016/TGU_7_48.pdf

12. Ferrara G., Preti U., Schena G. D. Modelling of screening operations // International Journal of Mineral Processing. 1988. Vol. 22, № 1/4. P. 193-222. https://doi.org/10.1016/0301-7516(88)90064-6

13. Rumyantsev S., Tarasov D. Numerical Simulation of Non-linear Dynamics of Vibration Transport Machines in Case of Three Independently Rotating Vibration Exciters // Recent Advances in Applied Mathematics: Proceedings of the American Conference on Applied Mathematics (AMERICAN-MATH'10). Cambridge: Harvard University, USA, January 27-29, 2010. Р. 191-194.

14. Subasinghe G. K. N. S., Schoap W., Kelly E. G. Modelling screening as a conjugate rate process // International Journal of Mineral Processing. 1990. Vol. 28. P. 289-300. https://doi.org/10.1016/0301-7516(90)90047-3

15. Тимошенко С. П. Колебания в инженерном деле. М.: Физматиздат, 1959. 439 с.

Статья поступила в редакцию 26 января 2019 г.

УДК 622.277:621.927 https://doi.org/10.21440/2307-2091-2019-2-93-99

Development and methods for determining the parameters of a vibration stand with an open-end deck plate

Arkadiy Vasil'evich YUDIN, Magomed Kazbekovich ABDULKARIMOV, Anatoliy Grigor'evich POPOV Viktor Stepanovich SHESTAKOV*

Ural State Mining University, Ekaterinburg, Russia

The relevance of the work is due to the need of the mining and construction industry in high-quality preparation of raw materials, increasing the yield of standard fractions, reducing the loss of mineral resources in the development of carbonate open pit. A problem of screening plants is the difficulty of separating rock mass containing clay inclusions. Serial inertial screens with a grizzly deck plate are often clogged with clay and reduce their effectiveness.

The purpose of the work: improvement of screens deck plates, increasing screening efficiency, reducing the degree of clogging by clay inclusions in the separation of a difficult-to-reach rock mass; development of a vibration stand for testing abstract theorems on the improvement of deck plates under experimental conditions.

Methodology of the research: development of methods for calculating the parameters of a stand on the basis of the vibration system, the implementation of the design solution and the preparation of methods for measuring the vibrations of the work tool and the deck plate. Results. A vibration stand with an improved deck plate has been developed. The stand includes: work tool, controlled-angle vibration exciters, flexible support, supporting frame, spring petal coupling, DC electric drive, control board, measuring equipment of the drive, sensor system for measuring oscillations and recording equipment. The peculiarity of the stand is that there are four supporting blocks between the sides of the work tool on which open-end grizzly are installed with varying parameters. Grizzlies form a four-stage deck plate.

Conclusion. A scheme has been developed for an experimental stand of avibrating grizzly with an improved deck plate. It is hypothesized that the screen with a open-end deck plate can be rightfully regarded as a dual-mass oscillatory system. The method of calculating the parameters of such a system is described. To evaluate the methodology, a program was developed and the results were analyzed. As a result of the work, the characteristics of the stand and the main design parameters are given.

Keywords: vibrating stand, screen, side, deck plate, grizzly, beam, hinge, elastic support, vibrator, model, amplitude, oscillation frequency.

REFERENCES

1. Yudin A. V., Shestakov V. S., Abdulkarimov M. K. 2018, Numerical analysis of oscillations of a dual-mass conservative vibrating screen system for a hard rock mass. Izvestiya vysshikh uchebnykh zavedeniy. Gornyi zhurnal [ News of the Higher Institutions. Mining Journal], no. 5, pp. 93-101.

2. Yudin A. V., 2016, Simulation of forced oscillations of a deck plate of an inertial screen with pinched grizzlies. Izvestiya vysshikh uchebnykh zavedeniy. Gornyi zhurnal [News of the Higher Institutions. Mining Journal], no. 6, pp. 63-70.

3. Stadkowski A., Yudin A., Komissarov A., Lagunova Yu., Akhmetova M., Stolpovskikh I. 2018, Calculation of parameters and design of the movable transfer station with vibrating screen feeder for the conveyor of deer queries. International Journal of Engineering and Technology (UAE), vol. 7, no 2, pp. 148-151. http://dx.doi.org/10.14419/ijet.v7i2.23.11904

4. Volkov E. B., Lyaptsev S. A. 2013, The influence of the angle of inclination of the working surface of the vibrating screen on the efficiency of screening. Sovremennyye problemy nauki i obrazovaniya [Modern problems of science and education], no. 4, p. 8 (In Russ.). URL: http://www. science-education.ru/ru/article/view?id=9642

5. Volkov E. B., Lyaptsev S. A. 2012, Computer simulation of screening process. Mezhdunarodnyy zhurnal eksperimental'nogo obrazovaniya [International journal of experimental education], no. 4, pp. 49-50.

6. Weisberg L. A. 1986, Proyektirovaniye i raschet vibratsionnykh grokhotov [Design and calculation of vibration screens]. Moscow, 144 p.

7. Weisberg L. A., Korovnikov A. N., Trofimov V. A. 2017, Modernization of technological cycles of screening on the basis of innovative equipment (to the 100th anniversary of the Mekhanobr Institute. Gornyy zhurnal [Mining journal], no. 1, pp. 11-17. (In Russ.) https://doi.org/10.17580/gzh.2017.01.02

8. Kartaviy A. N. 2013, Vibratsionnyye agregaty dlya pererabotki mineral'nogo i tekhnogennogo syr'ya. Modelirovaniye i elementy rascheta po kriteriyam energo- i resursoeffektivnosti [Vibration units for the processing of mineral and man-made materials. Modeling and calculation elements according to the criteria of energy and resource efficiency]. Moscow, 328 p.

9. Gazaleeva G. I., Tsypin E. F., Chervyakov S. A. 2014, Rudopodgotovka, drobleniye, grokhocheniye, obogashcheniye [Ore preparation, crushing, screening, enrichment]. Ekaterinburg, 914 p.

10. Nazarov K. S., Fet Sh. 2009, Analysis of modern design solutions that increase the efficiency of vibration classification of hard-to-get materials. Gorny Informatsionno-Analiticheskiy Byulleten

11. Shishkin E. V., Lebedev A. I. 2016, Study of parameters of the vibrating screen using tools of simulation modeling. Uchenyye zametki TOGU [Electronic scientific journal "Scientists notes PNU"], vol. 7, pp. 281-286 (In Russ.). URL: http://ejournal.pnu.edu.ru/media/ejournal/articles-2016/ TGU_7_48.pdf

12. Ferrara G., Preti U., Schena G. D. 1988, Modelling of screening operations. International Journal of Mineral Processing, vol. 22, no. 1/4, pp. 193-222. https://doi.org/10.1016/0301-7516(88)90064-6

13. Rumyantsev S., Tarasov D. 2010, Numerical Simulation of Non-linear Dynamics of Vibration Transport Machines in Case of Three Independently Rotating Vibration Exciters. Recent Advances in Applied Mathematics: Proceedings of the American Conference on Applied Mathematics (AMERICAN-MATH'10). Cambridge, Harvard University, USA, January 27-29, pp. 191-194.

14. Subasinghe G. K. N. S., Schoap W., Kelly E. G. 1990, Modelling screening as a conjugate rate process. International Journal of Mineral Processing, vol. 28, pp. 289-300. https://doi.org/10.1016/0301-7516(90)90047-3

15. Timoshenko S. P. 1959, Kolebaniya vinzhenernom dele [Oscillations in engineering]. Moscow, 439 p.

The article was received on January 26, 2019

" Н shestakov.v.s@mail.ru С https://orcid.org/0000-0001 -9468-6696