CC BY
23
УДК 631.3
DOI 10.18286/1816-4501-2019-1-6-11
ТЕОРЕТИЧЕСКИЕ ИССЛЕДОВАНИЯ ПРОЦЕССА ДОЗИРОВАНИЯ СЫПУЧЕГО
МАТЕРИАЛА
Семашкин Николай Михайлович1, кандидат технических наук, доцент кафедры «Математика и физика»
Исаев Юрий Михайлович1, доктор технических наук, профессор, заведующий кафедрой «Математика и физика»
Крючин Николай Павлович2, доктор технических наук, профессор, заведующий кафедрой «Механика и инженерная графика»
Крючин Александр Николаевич2, кандидат технических наук, доцент кафедры «Механика и инженерная графика»
ФГБОУ ВО Ульяновский ГАУ
432017, г. Ульяновск, бульвар Новый Венец, 1; тел.: 8 (84231) 55-95-49, e-mail: emotion.snm@
mail.ru
2ФГБОУ ВО Самарская ГСХА
446442, Самарская область, п.г.т. Усть-Кинельский, ул. Учебная, 2, e-mail:kryuchin@inbox.ru
Ключевые слова: частица, сыпучий материал, уравнение Навье-Стокса, дозирование, перемещение, пропускная способность.
В статье рассмотрено движение сыпучего материала на основании гидродинамической модели с использованием уравнения Навье-Стокса для вязкой жидкости в векторной форме. Рассмотрен слой перемещаемого материала, толщина которого значительно меньше по сравнению с его длиной и шириной. В теоретической части статьи описываются исследования движения сыпучего материала в дозирующем устройстве, где учитываются площадь и масса материала. А средняя скорость частиц материала ограничивается высотой поднятия штифтов, которая отличается от скорости вращения диска, так как одной из основных характеристик сыпучего материала является коэффициент внутреннего трения частиц материала. При этом для определения пропускной способности дозатора учитывалось выражение для определения средней скорости перемещения частиц материала. Получено уравнение, представляющее собой математическую модель неустановившегося режима движения, и проанализировано с учетом граничных и начальных условий. Найдено отношение скорости перемещения сыпучего материала к линейной скорости транспортера в различных точках по толщине слоя перемещаемого материала в различные моменты времени. В итоге выявлено, что использование штифтов для истечения сыпучего материала из загрузочного окна в дисково-штифтовом дозирующем устройстве позволяет обеспечить значительное увеличение его пропускной способности. Полученные теоретические исследования для определения скорости перемещения сыпучего материала вращающимся диском с подпружиненными штифтами позволяют получить оптимальную пропускную способность устройства при дозировании частиц материала в зависимости от конструктивных параметров этого устройства и физико-механических свойств материала. При максимальной высоте штифтов в средней части зоны загрузки материала пропускная способность устройства не превышала 0,55 кг/с.
Работа выполнена в рамках гранта Президента Российской Федерации для государственной поддержки молодых российских ученых - кандидатов наук МК-6675.2018.8
Введение
В сельскохозяйственном производстве большое значение имеет повышение качества дозирования сыпучих материалов. Одними из
основных направлений улучшения функционирования дозирующих устройств является теоретическое исследование и моделирование процессов переноса материала. Главным и суще-
ственным недостатком существующих средств дозирования следует считать высокие энергетические и материальные затраты на выполнение технологического процесса. Это обуславливает необходимость совершенствования средств дозирования сыпучих материалов, оптимизации их параметров и режимов работы, в том числе и на основе анализа математических выражений, описывающих процесс их работы.
Работа рассматриваемого устройства для дозирования сыпучих материалов заключается в следующем. Под действием вращающегося диска штифты движутся по поверхности направляющей шайбы, причем траектория их движения повторяет форму шайбы за счёт подпружиненных опорных концов (рис. 1). Во время движения штифты приподнимаются выше уровня козырька и далее поочерёдно входят в перемещаемый в устройстве материал. Вследствие этого своды, образовавшиеся возле загрузочного окна, разрушаются, и сыпучий материал равномерно распределяется по вращающемуся диску. Затем материал ровным слоем перемещается через выходное окно и скребком сбрасывается в эжекторное устройство [1].
объекты и методы исследований
Выполним теоретические исследования движения сыпучего материала в дозирующем устройстве. Учтем, что площадь и масса материала известны. Средняя скорость частиц материала, которая ограничена высотой поднятия штифтов, отлична от скорости вращения диска, так как одной из основных характеристик сыпучего материала является коэффициент внутреннего трения частиц материала. Для определения пропускной способности дозатора необходимо учитывать выражение для определения средней скорости перемещения частиц материала.
Движение сыпучего материала совершается за счет вращения диска (рис. 1). Ось X совпадает с направлением перемещения материала диском, а ось Z направлена вертикально вверх.
Для исследования движения исследуемого сыпучего материала рассмотрим его как вязкую жидкость с насыпной плотностью у, а кинематический коэффициент внутреннего трения частиц материала обозначим через п. Динамику движения материала выразим через уравнение Навье-Стокса для вязкой жидкости в векторной форме [2, 3, 4]:
^ 1 cUtS^ F--gradP + vV2s =—,
Y dt
(1)
Рис. 1 - Схема распределения скорости материала, перемещаемого по диску:
1 - ступеньки; 2 - подпружиненные штифты; 3 - вращающийся диск; 4 - козырёк
где F - сила тяжести, действующая на частицу, Н; у - плотность материала, кг/м3; P -нормальное давление на частицу, Па; v = п/y -коэффициент вязкости, м2/с; v2 - дифференциальный оператор Лапласа; и - скорость частицы, м/с; t - время перемещения, с.
Дифференциальный оператор Лапласа имеет вид:
д2 д2 + д2 дх2 ^
V2 =■
ду дх
2 '
где х, у, z - текущие координаты частицы. Спроецировав силы, действующие на частицу сыпучего материала по осям системы координат, перепишем уравнение в следующем виде:
' ёи
_х
ж
= X
TJ2 1 дР
v-V и---
р дх
du
^ = Y + v- Vu - 1дР dt у р ду
dvz ^ 1 дР
—^ = Z + v- V2uz---
dt z р дх
(2)
где и , и , и
1—1 x y z
проекции скоростей частицы сыпучего материала на соответствующие координатные оси, м/с; X, У, Z - проекция массовых сил, Н.
Поскольку движение материала совершается при вращении диска, то уравнение движения частицы материала рассматривается относительно подвижной системы координат. Тогда силы инерции, действующие на частицу, учитываются как массовые силы тяжести.
Так как ось Х направлена по направлению движения, а ось Z перпендикулярна к оси Х, то для режимов движения дозирующего устройства уравнения (2) запишем следующим обра-
зом:
и
Л
1 дР
а и
2
Я = 8
{р & . (3)
Рассмотрим особенности процесса дозирования в рассматриваемом устройстве:
а) относительная скорость перемещения частиц сыпучего материала во время движения зависит от сил внутреннего трения;
б) относительное перемещение материала зависит от параметров движения диска.
Движение потока частиц описывается первым уравнением, входящим в систему (3), если при этом не учитывать силы инерции перенос-
Рис. 2 - Зависимость отношений скоростей £ от координаты г вдоль оси по толщине перемещаемого слоя в различные моменты времени: 1 -1 = 1 с; 2 -1 = 3 с; 3 -1 = 10 с
Рис. 3 - Изменение отношения скоростей £ от координаты г вдоль оси по толщине перемещаемого слоя и времени (.
ного движения. Преобразуем пространственную задачу в одномерную, тогда дифференциальное уравнение примет вид:
а и а2 и
Ж
V
2
(4)
Зададим граничные условия движения диска, учитывающие факторы переносного движения. Обозначим толщину слоя перемещаемого материала через Л [5, 6, 7].
В качестве граничного условия движения материала примем, что при 1 = 0 их = и0, где и0 - линейная скорость на максимальном радиусе вращающегося диска, м/с.
Запишем граничные условия для перемещаемого слоя материала: при 1 = Л; их = 0. Для начальных условий (при t = 0) примем их = 0 (0 < 1 < Л).
Рассматриваемое выражение с определёнными начальными и граничными условиями является уравнением параболического типа. Тогда уравнение (4) с учётом граничных и начальных условий является математической моделью неустановившегося процесса дозирования.
Далее заменим искомую функцию суммой двух функций:
и(ъ,т) = w( х) + и (ъ,т) (5) w,, (ъ) = 0 . w(z)
где у 7 ; функция у 7 удовлетворяет граничным условиям:
w (0) = и0; w (И) = 0
а функция
и (х, т)
(6)
удовлетворяет уравне-
нию
д и ~дт
V-
д и д?
(7)
при следующих однородных граничных условиях:
и (0,т) = 0; и (Ь,т) = 0
(8)
и при начальном условии
и (ъ,0) = и(ъ,0) + w(z) = 0 и(ъ,0) = - w(z)
Окончательно решение
откуда (9)
уравнения
(4)
примет вид:
(
и( 2, * ) = и0
1 2 2 ^ 1 -1----У— ■ е
И п к
ук п
И1
. кп
81П-2
Л
(10)
С учётом того, что можно ограничиться двумя членами данного ряда, получим:
0.4
1-д
V *' У\3
\
\ » 1 1 • »
у ■2 V % * 4 •
\ » 1 • • » » *
\ \ \ \ 1 \ » - ч. _ V
г, мм
мм
Рис. 4 - Зависимость отношения скоростей £ от координаты г в различные моменты времени: 1 -1 = 1 с; 2 -1 = 3 с; 3 -1 = 10 с.
и (2,1) = и
, 2 2 . п -
1----81П— 2 • е
И п И
П 2 (
1 + е
п
• 008 — 2 И
(11)
Для анализа зависимости (11) определим отношение скоростей перемещения сыпучего материала к линейной скорости диска
е( 2, *) = ох (г, *)/ и0
по толщине слоя семян
в различные моменты времени:
~ п2 ( -■
, 2 2 . П -
1----81П—2 • е
И п И
Е
. (2, *) .
и
1 + е
п
• 008 — 2 И
(12)
Из решения выражения (12) и условия, что t находится пределах 0...10 с, следует, что режим движения материала является установившимся, а скорость по оси 1
\ - ^ И
( 2 , *) = и0
V
У
(13)
Полученные выражения объясняют тот факт, что в дозирующее окно под козырёк вначале попадают частицы, располагающиеся рядом с окном, а затем - следующие за ними [8, 9, 10].
Результаты исследований
Результаты расчетов для различных случаев расположения штифтов дозирующего устройства приведены на рисунках 2.7. На рисунке 2 представлена зависимость отношения скоро- Е = их / и0
стей х 0 от координаты 1 вдоль оси по
толщине перемещаемого слоя в различные моменты времени.
На рис. 3 показана зависимость изменения отношения скоростей Е(2, *) = их (2, *)/и от координаты 1 вдоль оси по толщине перемещаемого слоя и времени t.
Рис. 5 - Зависимость отношения скоростей £ от координаты г при ( = 1 с и различных высотах выхода штифтов Л: 1 - h = 2 мм;
2 - h = 4 мм; 3 - h = 6 мм; 4 - h = 8 мм
1
4
**
О 1.6 4.3 6.4 ¡} ш
Ь
Рис. 6 - Изменение пропускной способности 6 от высоты выхода штифтов Л: где 1, 4
- теоретические зависимости; 2, 3, 5 - экспериментальные зависимости с различными вариантами расположения штифтов: 2 - в начале окна; 3 - в середине окна; 5 - в конце окна
Зависимость
отношения
скоростей
Е = V / и от координаты 1 в различные моменты времени при выдвижении штифтов на высоту 6 мм показано на рис. 4.
Зависимость отношения скоростей
Е = их / и0 от координаты 1 при t = 1 с и выхода штифтов h на высоту от 2 мм до 8 мм представлено на рис. 5.
На рис. 6 показано изменение пропускной способности дозирующего устройства
О = Q / Q0 от высоты выхода штифтов Л.
В результате анализа представленных зависимостей выявлено, что большая пропускная способность дозирующего устройства обеспечивается при высоте штифтов более 6 мм и их расположении в начале и в конце зоны подачи материала.
Выводы
Использование штифтов для истечения сыпучего материала из загрузочного окна в дисково-штифтовом дозирующем устройстве позволяет обеспечить значительное увеличение его пропускной способности. Проведенные теоретические исследования по определению скорости перемещения сыпучего материала вращающимся диском с подпружиненными штифтами позволяют получить оптимальную пропускную способность устройства при дозировании частиц материала в зависимости от конструктивных параметров этого устройства и физико-механических свойств материала. При максимальной высоте штифтов в средней части зоны загрузки материала пропускная способность устройства не превышает 0,55 кг/с.
Библиографический список
1. Крючин, Александр Николаевич. Повышение качества посева семян трав самоходной пневматической мини-сеялкой применением дисково-штифтового высевающего аппарата: ав-тореф. дис. ... канд. технических наук: 05.20.01 / А.Н. Крючин. - Пенза, 2016. - 18 с.
2. Груздев, И.Э. Теории шнековых устройств / И.Э. Груздев, Р.Г. Мирзоев, В.И. Яиков. - Л.: Изд-во Ленинградского университета, 1978. - 144 с.
3. Гутьяр, Е.Я. Элементарная теория вертикального винтового транспортёра / Е.Я. Гутьяр // Научные труды МИМЭСХ им. В.М. Молотова. -М.: Машгиз, 1956. - Том 2. - С. 8-12.
4. Желтов, В.П. Расчёт спиральных винтовых конвейеров / В.П. Желтов // Вестник машиностроения. - 1975. - № 5. - С. 18-21.
5. SolidWorks 2007/2008. Компьютерное моделирование в инженерной практике / А.А. Алямовский, А.А. Собачкин, Е.В. Одинцов, А.И. Харитонович, Н.Б. Пономарев. - Спб.: БХВ-Петербург, 2008. - 1040 с.
6. Математическое моделирование при оптимизации параметров сушильной установки контактного типа / В.И. Курдюмов, А.А. Павлу-шин, С.Г. Мударисов, В.И. Долгов // Пермский аграрный вестник. - 2016. - № 2(14). - С. 107- 112.
7. Начальные скорости движения частицы материала при перемещении спиральным винтом / Ю.М. Исаев, Х.Х. Губейдуллин, Н.М. Семаш-кин, И.И. Шигапов // Аграрная наука. - 2014. - № 10. - С. 28 - 30.
8. Semashkin, N.M. Revision of the cosmic distance duality tests / N.M. Semashkin, A.V. Nikolaev // Проблемы теоретической и наблюдательной космологии. 5-ая Ульяновская международная школа-семинар. - Ульяновск: ГПУ им. И.Н. Ульянова, 2016. - 69 с.
9. Optimization of research and methodology work at university in terms of the process approach / I.D. Ibragimov, R.R. Iskhakova, M.A. Galeeva, M.M. Kalashnikova, Yu.V. Ryseva, I.I. Galimzyanova, I.A. Sharonov // Journal of Sustainable Development. -2015. - Том 8, № 3. - С. 234-241.
10. Васильев, В.Ф. Водометные движители / В.Ф. Васильев. - М.:МАДИ ГТУ, 2006. - 45 с.
THEORETICAL RESEARCH OF THE PROCESS OF BULK MATERIAL DOSING
Semashkin N.M.1, Isaev Yu.M.1, Kryuchin N.P.2, Kryuchin A.N.2 1FSBEI HE Ulyanovsk State Agrarian University 2 FSBEI HE Samara State Agricultural Academy 1432017, Ulyanovsk, Novyi Venets Boulevard, 1; tel.: 8 (84231) 55-95-49, e-mail: emotion.snm@mail.ru 2446442, Samara region, Ust-Kinelskyv., Uchebnaya st., 2, e-mail: kryuchin@inbox.ru
Key words: particle, bulk material, Navier-Stokes equation, dosing, movement, throughput.
The article describes the movement of bulk material based on the hydrodynamic model using the Navier-Stokes equation for viscous fluid in vector form. A layer of transported material is studied, its thickness is significantly smaller compared to its length and width. The theoretical part of the article describes studies of the movement of bulk material in a dosing device, which takes into account the area and mass of the material. The average speed of the particles of the material is limited by the height of the raising of pins, which differs from the speed of disk rotation, since one of the main characteristics of the bulk material is the coefficient of internal friction of the particles of the material. At the same time, to determine the throughput of the measuring hopper, an equasion was taken into account to determine the average speed of movement of particles of a material. An equation was derived that represents a mathematical model of an unsteady motion regime and it was analyzed taking into account boundary and initial conditions. The ratio of the velocity of the bulk material to the linear velocity of the conveyor at various points along the thickness of the layer of the material being moved at different time is found. As a result, it was revealed that the use of pins for the flow of bulk material from the loading chamber of the disk-pin dosing device allows to provide a significant increase in its throughput. The obtained theoretical studies of the movement speed specification of the bulk material by a rotating disk with spring pins allow to obtain a suitable throughput of the device when dosing particles of the material depending on the design parameters of this device and the physical and mechanical properties of the material. With the maximum height of the pins in the middle part of the material loading zone, the throughput of the device did not exceed 0.55 kg /s.
Bibliography
1. Kryuchin, Alexander Nikolaevich. Improving the sowing quality of grass seeds with self-propelled pneumatic mini-seeder using disc-pin seeding device: author's abstract of dissertation of Candidate of Technical Sciences: 05.20.01 /A.N. Kryuchin. - Penza, 2016. -18 p.
2. Gruzdev, I.E. Theories of screw devices / I.E. Gruzdev, R.G. Mirzoev, V.I. Yaikov. - L.: Publishing House of the University of Leningrad, 1978. -144 p.
3. Gutyar, E.Ya. Elementary theory of vertical screw conveyor / EYa. Gutyar // Scientific works of Moscow Institute of Mechanization and Electrification of Agriculture named after V.M. Molotov. - M.: Mashgiz, 1956. - Volume 2. - P. 8-12.
4.10. Zheltov, V.P. Calculation of spiral screw conveyors / V.P. Zheltov // Vestnik of mechanical engineering. -1975. - № 5. - P. 18-21.
5. SolidWorks 2007/2008. Computer modeling in engineering practice / A.A. Alyamovsky, A.A. Sobachkin, E.V. Odintsov, A.I. Kharitonovich, N.B. Ponomarev. - SPb.: BKhV-Petersburg, 2008. -1040 p.
6. Mathematical modeling when improving the parameters of a contact-type drying installation / V.I. Kurdyumov, A.A. Pavlushin, S.G. Mudarisov, V.I. Dolgov//Perm Agrarian Journal. - 2016. - № 2 (14). - P. 107-112.
7. The initial velocity of the particle of the material when moved by a spiral screw / Yu.M. Isaev, Kh.Kh. Gubeidullin, N.M. Semashkin, I.I. Shigapov// Agrarian science. - 2014. - № 10. - P. 28 - 30.
8. Semashkin, N.M. Revision of the cosmic distance duality tests / N.M. Semashkin, A.V. Nikolaev // Problems of theoretical and observational cosmology. 5th Ulyanovsk International training Seminar. - Ulyanovsk: Teacher Training University named after I.N. Ulyanov, 2016. - 69 p.
9. Optimization of research and methodology work at university in terms of the process approach / I.D. Ibragimov, R.R. Iskhakova, M.A. Galeeva, M.M. Kalashnikova, Yu.V. Ryseva, I.I. Galimzyanova, I.A. Sharonov // Journal of Sustainable Development. - 2015. - Volume 8, No. 3. - P. 234-241.
10. Vasiliev, V.F. Water jets / V.F. Vasiliev. - M.: Moscow Automobile and Highway State Technical University, 2006. - 45 p.
