CC BY
51
Взаимодействие частиц минеральных удобрений с криволинейными лопатками вертикального ротора
И.А. Казачков, А.Ю. Попов
Дробление частиц минеральных удобрений, как правило, приводит к снижению их эффективности [1]. Оно сопровождается увеличением доли мелких частиц в общей массе, что ведет к снижению равномерности рассева
[2, 3].
Для устранения неблагоприятных условий удара частиц удобрений о рабочие органы ротора рассмотрим возможность применения лопаток в форме окружности в плане. Этот вопрос неоднократно рассматривался в работах целого ряда авторов, применительно к центробежно-дисковым аппаратам [4 - 6].
Представим рабочую поверхность лопатки в виде части окружности радиуса гл, и отклоним назад на угол у о. Для улучшения процесса движения частиц и увеличения скорости метания, расположим рассматриваемую лопатку вогнутой поверхностью по направлению вращения ротора.
Схема сил, действующих на движущуюся частицу удобрения по вогнутой поверхности лопатки, показана на рис. 1.
Дифференциальное уравнение движения частицы по поверхности лопатки имеет вид [7]:
= м>2Яп соБу + gсоб(у±^) - /2Яп Бту + gБт(у±^) - 2]. (1)
Преобразуем уравнение (1), для чего рассмотрим совокупность сил, действующих на частицу при движении по вогнутой поверхности лопатки (рисунок 1), из которого следует, что
Яп собу = СбшЗ; Яп Бту = СсобЗ- гл,
7 2 2
гл + Го + 2гл Го бш уо - радиус кривизны поверхности
Е
лопатки, м; 3 = —+ 30 - центральный угол, пройденного пути частицы, град;
Г
л
а0 = агеэт 0 ^—- - начальное значение центрального угла &, град; угловая скорость.
Рис. 1. - Схема сил, действующих на частицу при движении по вогнутой поверхности лопатки
Подставив полученные значения в уравнение (1), получим нелинейное дифференциальное уравнение второго порядка следующего вида:
4 = g[с[ ± А)) -1 ± Ао)] +
с sinа+/
с cos а - гл -
24
/Л
(2)
Каши. Для этого введем обозначения
-/
в результате получим систему
Для получения численного решения данного дифференциального уравнения (2) необходимо его преобразовать в систему двух уравнений первого порядка [8], то есть преобразовать нелинейное уравнение к форме
У1 = 4 ~
[
У2 = 4
уравнений первого порядка относительно функций пути 4 и относительной скорости 4 = ¥г:
' у! = V;
У 2 = g + Мо) -1 + Ао)]
+ м
с sinа+/
С cosа- г„ -
24
/
(3)
Решение полученной системы уравнений (3) базируется на методе численного интегрирования Рунге-Кутты второго порядка с коррекцией в средней точке. Основная идея данного метода заключается в том, что производные аппроксимируются через значение самой искомой функции, в интервале выбранного шага, для перехода из точки в другую требуется дважды вычислить правую часть дифференциального уравнения, чтобы обеспечить наибольшую близость к ряду Тейлора.
Решение системы уравнений при начальных условиях t = 0, % = ¡1,
Е = иау дает возможность вычислить время движения частицы по
поверхности лопатки на пути, равном % = гл3л и скорость Уг до момента
схода частицы с конца лопатки. Для определения абсолютной скорости и ее направления в момент схода с лопатки, введем дополнительное условие
>1 ^ I.
Для определении угла и направления схода необходимо учесть тот факт, что время t1 соответствует продолжительности движения частицы
удобрения по лопатке рассеивателя с радиальными лопатками. У рассеивателей с наклонными лопатками, тем более с лопатками в виде части окружности в плане, значение времени t1 необходимо понимать как условное время. С действительным временем t движения частицы удобрения угол и направление схода связаны соотношением:
wt = wtl ±(у-3л + ак), (4)
С вт(3л +3о)
где аК = агссоБ--угол наклона лопатки к радиусу ротора,
К
с2 + г2 - К2
град; 3Л = агссоБ-л--30 - центральный угол раствора лопатки,
2С • г„
л
град.
Решение уравнений (3), (4) и анализ полученных зависимостей показывает, что на абсолютную скорость и величину угла схода частиц удобрений с конца лопатки оказывают влияние кинематические параметры,
начальная скорость поступления частиц на лопатку Уп, скорость вращения ротора /г0, а так же геометрические параметры, угол наклона лопаток относительно начального радиуса ц0, радиус кривизны поверхности гл и центральный угол раствора Зп.
Уг.
м / с 1 0
9
8
7
6
5
4
^ , град
70 65 60 55 50
^Уг
У?
1 - ' V = 0
Ф 45 п
40 2 - = 0,16
'Х 35 3 - у?/ 30
/Я
= 0,34
ц 0, град
Рис. 2. - Зависимость абсолютной скорости и угла схода от геометрических параметров лопаток при различных значениях кинематического показателя В результате анализа можно видеть, что увеличение кинематического
показателя Уу/я вызывает уменьшение угла схода и рост абсолютной скорости (рис 2). Изменение угла схода wt от кинематического показателя Уу/я близка к линейной зависимости. При отклонении лопатки назад, то
/ н
есть увеличении угла наклона лопатки к начальному радиусу, эта зависимость постепенно смещается в сторону уменьшения углов схода. Примерно тоже, происходит с зависимостью абсолютной скорости, которая с ростом угла отклонения ц0 вызывает увеличения значений абсолютной
скорости. При этом вид зависимости абсолютной скорости, позволяет сделать вывод о наличии экстремума, то есть дальнейшее увеличение угла наклона вызовет снижение абсолютной скорости.
Анализ влияния геометрических параметров кривизны лопаток (рис. 3) на величину угла схода и абсолютной скорости показал, что отклонение лопаток назад по ходу вращения ротора, вызывает уменьшение угла схода и увеличение абсолютной скорости. Необходимо отметить, что зависимость угла отклонения лопатки у0 от абсолютной скорости Уг имеет экстремум, положение которого определяется кинематикой работы роторного рассеивателя [9].
м / с 10
град 60
55
50
45
40
35,
30
-х-Ч £
1 -3л = 38,2 '2 -3 = 20,83
X 3 -&л = 18,09 ■Р 25
у 0, град
Рис. 3. - Зависимость абсолютной скорости и угла схода от геометрических параметров криволинейной поверхности лопатки Полученные зависимости позволяют выявить влияние на величину абсолютной скорости и угла схода при работе рассеивателя с криволинейными лопатками, таких факторов как радиус кривизны поверхности, центральный угол раствора, угол отклонения лопатки, а так же параметров поступления частиц на лопатки. Для того, чтобы учесть влияние последних, рассмотрим процесс поступления частиц на поверхность лопатки Процесс поступления частицы удобрений с туконаправителя на лопатки ротора происходит с некоторой начальной скоростью Уп на радиусе г0 - Яп с учетом свободного падения до момента встречи с лопаткой на расстоянии £ [10, 11].
Для определения максимального расстояния ? (рис. 4), на котором произойдет встреча частицы при свободном падении, воспользуемся формулой равноускоренного движения.
? = У^ вш х +
2
(5)
где t - время падения частицы до встречи с лопаткой, с; % - угол наклона туконаправителя к горизонту, град.
Рис. 4. - Схема взаимодействия частиц удобрений с криволинейной поверхностью с учетом наклона лопатки и ширины потока
wt = — ± Л + т, 2
(6)
где Л = ц0 - агевт
Я
центральный угол, учитывающий
г0в1п ^0
г
отклонение лопатки к начальному радиусу 0, град; т = ¡и0 - агевт
Я
центральный угол, учитывающий ширину потока, град; 2 - количество лопаток.
С учетом кривизны поверхности, решение поставленной задачи выглядит в виде системы уравнений:
--
2п
г
2 • 3
2 Г 81П—1 л 2
3
У0 - у + ^0 - агсзт-£
Г0 81И
У0
33 2
Яп
- агсБт- 0
Яп
сов
31
У0 - — + М0
где Яп = д/С2 + гл2 - 2Сгл соб(30 + 31) - радиус точки встречи частицы с лопаткой, м.
Решение полученной системы уравнений дает время свободного падения частиц до момента встречи с лопаткой и центральный угол 31 этой встречи.
Таким образом, применение лопаток с формой поверхности в виде части окружности в плане позволяет уменьшить дробление частиц минеральных удобрений за счет уменьшения удара.
Литература:
1. Забродин, В.П. Анализ взаимодействия частиц с лопаткой наружного диска распределительного рабочего органа [Текст] // Совершенствование технологических процессов и конструкций сельскохозяйственных машин. - Краснодар, 1989. - Вып. 294(322). - С. 8491.
2. Петренко С. С. Определение конструктивных параметров шнекового смесителя сыпучих материалов [Электронный ресурс] // «Инженерный вестник Дона», 2013, №1 - Режим доступа: http://www.1vdon.ru/magaz1ne/arch1ve/n1y2013/1536 (доступ свободный) -Загл. с экрана. - Яз. рус.
3. Воронин В.В., Адигамов К.А., Петренко С.С., Сизякин Р.А. Критерии и способы оценки качества смешивания сыпучих материалов [Электронный ресурс] // «Инженерный вестник Дона», 2012, №4 (часть 2) -
Режим доступа: http://www.ivdon.ru/magazine/archive/n4p2y2012/1400 (доступ свободный) - Загл. с экрана. - Яз. рус.
4. Забродин, В.П. Распределительный рабочий орган разбрасывателей удобрений [Текст] // Механизация и электрификация сел. хоз-ва. - 1985. - №7. - С. 25-27.
5. Olieslagers R., Ramon H., Baerdemaeker J. De. Calculation of Fertilizer Distribution Patterns from a Spinning Disc Spreader by means of a Simulation Model [Text] // Journal of Agricultural Engineering Research. - 1996. - №63 (2), pp. 137-152.
6. Van Liedekerke P., Tijskens E., Dintwa E. DEM simulations of the particle flow on a centrifugal fertilizer spreader [Text] // Powder Technology. -2009. - №190 (3), pp. 348-360.
7. Хаджиев А.Х., Дадаходжаев А. Некоторые вопросы теории криволинейного движения частиц минеральных удобрений по неподвижной горизонтальной плоскости [Текст] // Вопросы механизации и электрификации сел. хоз-ва. - Ташкент, 1982. - Вып. 23. - С. 3-13.
8. Мудров А.Е. Численные методы для ПЭВМ на языке Бейсик, Фортран и Паскаль [Текст] / А.Е. Мудров. - Томск: МП «Раско», 1991. - 272 с.
9. Черноволов, В. А. К определению угла сектора рассева центробежного дискового аппарата с радиальными лопатками [Текст]. // Тр. АЧИМСХа. Сер. Механизация сельского хозяйства. - 1971. - Вып. 20.- С. 4648.
10. Жилин, А.П. Исследование движения минеральных удобрений в роторе с горизонтальной осью вращения [Текст] / Исследование, проектирование и производство рабочих органов сельскохозяйственных машин. - Ростов на Дону, 1980. - С. 35-38.
11. Жилин А.П. К определению некоторых параметров роторной разгрузочной тележки для склада минеральных удобрений [Текст] /
Комплексная механизация и автоматизация сельскохозяйственного производства. - Ростов на Дону, 1981. - С. 130-134.
