CC BY
27
УДК 631.316.02-001.5
Центробежное фракционирование длинностебельных материалов
И. В. Игнатенко, А. В. Погорелов
(Донской государственный технический университет)
Для фракционирования длинностебельных сельскохозяйственных материалов по длине предложен принцип центробежного фракционирования без применения сит. Описано устройство бесситового центробежного фракционера и его расчётная схема. Определены основные процессы: разгон частиц центробежным диском, разбрасывание в кольцевом аэродинамическом канале и разделение разделительными цилиндрами. Произведено математическое моделирование процессов разгона частиц центробежным диском и движения сброшенных частиц в восходящем воздушном потоке. Проведён компьютерный эксперимент. Получены траектории движения частиц в кольцевом аэродинамическом канале. Отличие траекторий так значительно, что можно говорить о сепарационных свойствах сочетания центробежного поля сил и воздушного потока. Приводятся рекомендации по компоновке разделительных цилиндров. Результаты моделирования свидетельствуют о принципиальной возможности реализации процесса фракционирования дпинносгебельных материалов центробежным рабочим органом за счёт различий траекторий разброса частиц в воздушном потоке. Ключевые слова: переработка сельскохозяйственного сырья, центробежное фракционирование, воздушный поток, математическое моделирование, компьютерный эксперимент.
Введение. В сельскохозяйственном производстве возникает необходимость разделения длинностебельных материалов типа соломы по длине на фракции [1]. Широкое распространение получило использование соломистой подстилки при напольном выращивании бройлеров. Однако солома, прошедшая простое измельчение, не соответствует зоотехническим требованиям по длине резки (рекомендованная длина — от 100 до 180 мм). Слишком короткие фракции представляют опасность для птицы при заглатывании. Длинная фракция не образует слоя нужной плотности.
Состояние проблемы. Возникает задача механизации фракционирования соломистой массы, получаемой при уборке зерновых культур в качестве отходов. Фракционирование длинностебельной массы не имеет должного обеспечения техническими средствами. Типовые технологии фракционирования продуктов основаны на применении различных сит. Для длинностебельных материалов ситовые технологии неработоспособны. Требуется разработка бесситовых технологий.
Цель статьи — обоснование использования для фракционирования соломистой массы центробежного бесситового рабочего органа.
Принцип центробежного бесситового фракционирования. Принцип центробежного бесситового фракционирования основан на аэродинамической сепарации потока частиц, разбрасываемого быстровращающимся диском. В нём совмещены процессы разгона частиц на центробежном диске и воздушная сепарация сбрасываемых с диска частиц в аспирационном кольцевом канале. Функциональная схема центробежного фракционирования соломистых материалов представлена на рис. 1.
Основными рабочими органами машины являются быстровращающийся диск 2 и разделительные цилиндры 5. Длинностебельная масса после измельчения через питатель 1 подаётся в центр диска 2, приводимый во вращение от электродвигателя 3, и под действием центробежных сил движется по нему, набирая скорость. Частицы, сбрасываемые с диска, попадают в кольцевой зазор, продуваемый потоком воздуха от вентилятора 4. Траектории частиц в зазоре зависят от их свойств — в частности от парусности. Лёгкие частицы будут уноситься с воздухом в отводящие патрубки, остальные будут попадать в разные цилиндры 5, распределяясь по фракциям.
Ворох
Рис. 1. Функциональная схема центробежного фракционирования:
1 — питатель; 2 — диск; 3 — электродвигатели; 4 — вентилятор; 5 — разделительные цилиндры
Расчётная схема процесса центробежного фракционирования. Расчётная схема (рис. 2) представляет собой плоский диск, вращающийся на вертикальном валу с угловой скоростью со. Поток стебельной массы считается состоящим из отдельных длинных частиц, не взаимодействующих друг с другом, имеющих разную массу т и размеры: длину /, ширину Ь, толщину б. Частица считается длинной, если её толщина и ширина намного меньше длины I Длинная частица поступает на диск в точку с начальным радиусом г0 (положение 1). Ориентация её характеризуется углом ф! случайной величины. Силами трения частица увлекается вращающимся диском в сложное движение: переносное вращательное вместе с диском и относительное поступательное по диску.
Помимо сил трения на частицу действуют центробежные силы, распределённые по её длине неравномерно, так что образуется крутящий момент Мк. Под действием центробежных сил частица движется к краю диска (положение 2), одновременно разворачиваясь, уменьшая угол ориентации (ф2 < Ф1). Её устойчивое состояние наступает, когда частица расположится по радиусу (положение 3). Тогда угол ориентации ф3 = 0, закручивающий момент Мк = 0.
По достижении края диска частица сбрасывается с него с абсолютной скоростью 1/а и движется в воздушном канале между диском и разделительными конусами, подвергаясь воздействию воздушного потока со скоростью и. Под действием аэродинамических сил частица (точнее, её центр тяжести) описывает траекторию. Так как частицы имеют разную длину, массу и парусность, то и траектории их в воздушном потоке будут разные. Лёгкие частицы будут выдуваться наверх, тяжёлые, длинные — падать вниз.
Разделение на фракции по длине осуществляется кольцевыми каналами, расположенными на пути потока сброшенных с диска частиц.
Рис. 2. Расчётная схема центробежного фракционирования
В расчётной схеме присутствуют три разных процесса:
1) разгон — движение длинной шероховатой частицы по вращающемуся диску;
2) разбрасывание — движение частицы в воздушном потоке с начальной скоростью по разным траекториям;
3) разделение потока частиц на фракции по кольцевым каналам.
Математическая модель фракционирования включает модели этих трёх процессов.
Математическое моделирование движения частицы по диску. Движение частицы по вращающемуся диску — хорошо проработанный в литературе вопрос [2]. Его теория, разработанная академиком П. М. Василенко, представлена в учебниках. Однако теория предполагает, что частица круглая, и требует доработки, если речь идёт о длинной частице.
Представим относительное движение частицы по диску как смещение центра тяжести длинной частицы и её вращение в плоскости диска вокруг центра тяжести.
Смещение центра тяжести ^ рассматриваем аналогично задаче П. М. Василенко как движение точечной массы т по горизонтальному диску с прямыми лопастями, установленными под углом ф к радиусу.
Уравнение движения центра масс по Василенко П. М.
5" + - о)2$ = [л0со2 cos (ф0 + ф)/ cos ф ] - fg
Здесь f— коэффициент трения частицы о диск; ср — угол трения; tgcp = f.
Уравнение (1) является неоднородным дифференциальным уравнением второго порядка. Его характеристическое уравнение:
Л2 + 2/соЛ - со2 = О
имеет корни
\ = со • ^(l + f2) - Ао, Л2 = -со • Л/(ГТ77) - Ао.
Полное решение уравнения (1) с начальными условиями: при t = 0 Ъ, = % =0 имеет вид:
1
(1)
(2)
(3)
fg__r cos(ф0 ±ф) со2 0 COS ф
Л2 Aj
(A2eA‘f
(4)
(5)
Переносная скорость представлена окружной тангенциальной скоростью 14 = до. Абсолютная скорость рассева частиц
Скорость относительного движения частицы вдоль направления лопасти:
5 = и =
fg ^ COS (ф0 +ф)" ^1^2 (a\t aK,t \
со2 Г0 coscp 1 J J *< 1 *< 1
V, =V, +К .
(6)
Вращение длинной частицы вокруг центра тяжести описывается дифференциальным уравнением вида:
Лр = Мк-Мв, (7)
где 3 — центральный момент инерции длинной частицы; МК = %яхо2Дл0 + ^¡пф — момент распределённой центробежной нагрузки; Мв= Ка1\соэф — момент аэродинамических сил сопротивления; Ка— коэффициент аэродинамического сопротивления.
Уравнения (1) и (7) составляют полную математическую модель движения длинной частицы по диску.
Анализ уравнений показывает, что к моменту сброса с диска длинная частица стремится сориентироваться по направлению действия центробежных сил по радиусу диска. Этот центробежный эффект упорядочивает последующее движение сброшенных частиц в воздушном потоке. Анализ не обнаруживает дестабилизирующих факторов, способных вызвать кувыркание частицы. Поворот длинных сориентированных частиц просматривается только в горизонтальной плоскости. Моделирование процессов разбрасывания. После сброса с диска частица движется в равномерном воздушном потоке со скоростью и/ под действием силы тяжести тд и силы лобового сопротивления Я, направление которой противоположно направлению скорости V относительно воздушного потока (рис. 3). Для плоской задачи начальные условия: при Ь = 0 щ = 1/Л= 0. Уравнения движения будут:
тх = -Л соэ 0 , (8)
тг = -ЛэтЭ + тд , (9)
где 0 — угол между осью Х\л направлением относительной скорости V
Сила сопротивления Л противоположна направлению относительной скорости частицы относительно потока V. Её величина пропорциональна квадрату относительной скорости, т. е.
Л = ккУ2. Вместе с силой веса С сила сопротивления Я образует равнодействующую силу Р,
направленную против скорости движения частицы и вызывающую торможение частицы. На рис 3. показано равновесие сил в момент сброса частицы с диска и в текущем положении. Сравнение
показывает, что по мере движения в воздушном потоке сила ^ уменьшается, но эффект её остаётся тормозящим.
Проекции вектора относительной скорости на оси координат будут:
Vx = х + u cos 3; Vz =z + í/sinp . (10)
Выражение для модуля вектора скорости V
И = № + Vz = л!ц2 + *'2 + Z2 + • (х ■ COS 3 + Z • sin 3) (11)
Коэффициент сопротивления kR будем считать пропорциональным площади лобового сечения 5 частицы. kR = kS. При скорости воздушного потока, близкой к скорости витания частицы, можно принять k= const.
Полное выражение для силы сопротивления:
R=kS V2, (12)
где Vнаходится по выражению (11).
Выражения для sin0 и cos0 получим через проекции V:
sine = 14/1/, cose = Vz/V. (13)
С учётом соотношений (10) получим:
Sin 0 = (z + Í7 Sin 3) /1/, COS0 = (л- + Í/COS3) /V . (14)
Полученные соотношения (12) и (14) приводят дифференциальные уравнения задачи (8) и (9) к
виду:
тх = —к • S - V (х + и cos 3) (15)
mz = mg -к S -V {z + u - sin 3) (16)
Если учесть, что объём частицы О = a'b'i, то масса частицы
т= Оу = a'b'íy (17)
Тогда площадь лобового сечения 5 можно выразить через массу частицы и её толщину, являющуюся квалитетом сепарации:
S = ашЬ = т/(у{) (18)
С учётом соотношений (11) и (18) и уравнений (15) и (16) получим математическую модель процесса разброса в виде системы нелинейных дифференциальных уравнений второго порядка:
х = ~—(х + и соэ^и2 + х2 + г2 + 2у-(л--С05р + г-5тр) (19)
у/7
г = д -—{г+ и%\пр)^и2 + х2 + ¿2 + 2и ■ (х ■ соэ^ + г ■ Бтр) (20)
р(
Отличительной особенностью уравнений является наличие в них членов, зависящих от длины частицы I
Уравнения трансцендентны и алетического решения не имеют. Для построения траекторий относительного движения частиц полученные уравнения решали методом компьютерного моделирования в среде 51ти1тк Ма^аЬ.
Рассчитаны траектории частиц соломистой массы. На рис. 4 представлены траектории абсолютного движения ДА) частиц в равномерном воздушном потоке при режиме сепарации: и = 8 м/с, V, = 20 м/с, 3 = 70°.
-@-Ряд1 -0—Ряд2 -6-РядЗ -е-Ряд4 -х-Ряд5
Рис. 4. Траектории движения соломистых частиц в воздушном потоке при центробежной сепарации. Длины фракций: ряд 1 — 250 мм; ряд 2 — 180 мм; ряд 3 — 100 мм; ряд 4—15 мм; ряд 5 — 5 мм
Результаты моделирования говорят о большой разнице в траекториях частиц лёгкой и тяжёлой фракций. Траектории отличаются так сильно, что можно говорить о сепарационных свойствах сочетания центробежного поля сил и воздушного потока.
Следует отметить быстротечность процесса сепарации. Время сепарации при скорости воздушного потока, принятой в пневмосепарирующих системах для многокомпонентного вороха, составляет по расчётам 0,1—0,3 с, что намного меньше, чем при сепарации на решётах. Следовательно, можно ожидать высокой производительности центробежного фракционирования.
Моделирование процесса разделения на фракции. Различия траекторий частиц различной массы и размеров в центробежном поле позволяют осуществлять разделение частиц длинностебельного материала. Разницу в дальности выброса самой лёгкой и тяжёлой частиц следует считать мерой фракционных свойств центробежной установки.
Процесс разделения разлетающихся частиц по фракциям осуществляется разделительными цилиндрами. Входные параметры — гранулометрический состав, радиусы сепарирующих цилиндров, высота кольцевых зазоров. Выходные — гранулометрический состав в каждом кольцевом зазоре. Зависимость между ними описывается методами статистики при помощи гистограмм и является наиболее информативным показателем фракционных свойств установки. Однако получение таких показателей достаточно трудоёмко.
Более наглядным для практики показателем является толщина потока. Величина толщины потока равна расстоянию между точками пересечения крайних траекторий с наклонной образующей конической поверхности разделителей. Угол наклона рекомендуется принимать 30 градусов.
Толщина потока увеличивается с увеличением диаметров конуса разделителя. Однако минимальный радиус, определяющий начало образующей конуса, не должен превышать дальность полёта самой тяжёлой фракции на уровне условного пола.
Проведённые расчёты по структурным моделям в среде Simulink показали, что при изменении длины соломистых частиц от 250 до 100 мм дальность полёта фракции составляет примерно 1,5 высоты диска над уровнем условного пола. Этого позволяет делать окна разделительных конусов достаточно большими, чтобы избежать забивания. Устанавливая разделительные перегородки через 300 мм, можно легко получать необходимое разделение на 3 фракции, что в большинстве случаев вполне приемлемо для практики.
Выводы. Результаты моделирования разброса частиц свидетельствуют о принципиальной возможности реализации процесса фракционирования длинностебельных материалов центробежным бесситовым разбрасывающим рабочим органом.
Библиографический список
1. Московский, М. Н. Вопросы экономического применения соломистой подстилки в ПТФ в условиях ЮФО / М. Н. Московский, А. Д. Гришков, А. В. Погорелов: мат-лы междунар. научн.-практ. конф. 4—5 марта 2010 г., Ростов-на-Дону, 2010. — С. 346—350.
2. Теория, конструкция и расчёт сельскохозяйственных машин / Е. С. Босой и др. — Москва: Машиностроение, 1977. — 568 с.
Материал поступил в редакцию 08.04.2012.
References
1. Moskovskij, М. N. Voprosy' e'konomicheskogo primeneniya solomistoj podstilki v PTF v usloviyax YuFO / M. N. Moskovskij, A. D. Grishkov, A. V. Pogorelov: mat-ly' mezhdunar. nauchn.-prakt. konf. 4—5 marta 2010 g., Rostov-na-Donu, 2010. — S. 346—350. — In Russian.
2. Teoriya, konstrukciya i raschyot sel'skoxozyajstvenny'x mashin / E. S. Bosoj i dr. — Moskva: Mashinostroenie, 1977. — 568 s. — In Russian.
CENTRIFUGAL FRACTIONATION OF LONG-STALKED MATERIALS I. V. Ignatenko, A. V. Pogorelov
(Don State Technical University)
The screenless centrifugal fractionation approach is offered for the fractionation of the long-stalked agricultural materials lengthways. The screenless centrifugal fractionation device and its analytical model are described. The basic processes are specified - accelerating particles with the centrifugal disc, spreading in the annular aerodynamic channel, and assorting with the partitioning cylinders. The mathematical simulation of the particle acceleration with the centrifugal disc and the dropped particles motion in the un ward airflow processes is carried out. The computer-based experiment is conducted. Paths of particles in the annular aerodynamic channel are obtained. The paths differ so greatly that it is possible to state the separating properties of the centrifugal force field and airflow combination. The recommendations on the arrangement of the partitioning cylinders are given. The simulation results imply the feasibility of assorting the long-stalked materials with the centrifugal tool at the cost of the differences in particle spreading paths in the airflow.
Keywords: agricultural raw materials processing, centrifugal fractionation, air flow, mathematical simulation, com-puter-based experiment.
