Сборник статей
1 1 1 1 1 1 1 1 1 1 Рейтинг 4.00 (1 Голос)

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

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

Цель настоящей статьи – выполнить расчетные оценки показателя качества сепарации комплексной аэромеханической безрешетной сепарирующей машины и рассмотреть способы повышения эффективности процесса сепарации в потоке воздуха. Прототипом проточной части сепарирующей машины выбрана проточная часть машины «АЛМАЗ» [6]. Именно эта машина по заявленным характеристикам осуществляет одновременное удаление примесей и разделение зерна на фракции, благодаря различному горизонтальному смещению частиц свободно падающего вороха зерна в горизонтальном сносящем потоке воздуха (см. рис.1).

Схема работы сепарирующей машины

Рис.1. Схема работы сепарирующей машины «АЛМАЗ» [5].

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

image002_43_0f3eba13c1c094c01192061966801cea Расчетный анализ внутренней аэродинамики пневмосепарирующей машины.

; (1)

image003_31_83cd2bdf1bb788cb285cf5a4c874a1f5 Расчетный анализ внутренней аэродинамики пневмосепарирующей машины.

, (2)

где h – высота падения, м; l – величина горизонтального смещения, м; t – время падения, с; Wx, Wy – горизонтальная и вертикальная составляющие скорости движения частицы, м/с; Wxо, Wyо – горизонтальная, направленная против движения сносящего потока, и вертикальная составляющие скорости движения частицы в начальный момент времени, м/с; Wcx, Wcy - горизонтальная и вертикальная составляющие скорости сносящего потока, м/с; k=(Cx· ρ ·S)/(2·m), ρ – плотность воздуха, кг/м3; Cx – коэффициент аэродинамического сопротивления частицы; m – масса частицы, кг; S – площадь поперечного сечения сферической частицы, м2.

При постоянстве во времени скорости сносящего потока и величины коэффициента аэродинамического сопротивления решением уравнений (1) и (2) будут выражения:

image004_29_fbd4bf608d50a29de011def3070a3a1e Расчетный анализ внутренней аэродинамики пневмосепарирующей машины.

; (3)

image005_26_11e54a4ae3e7fc91bc100e7926bb6417 Расчетный анализ внутренней аэродинамики пневмосепарирующей машины.

. (4)

Вычисленные по (3) и (4) высота h и длина l позволяют строить траектории движения сферических частиц во времени. Решение уравнения (3) относительно t определяет время падения частицы для заданной высоты h:

image006_20_c79e68e0216a0af1bffc3895cb7868d6 Расчетный анализ внутренней аэродинамики пневмосепарирующей машины.

,где image007_18_1da9972174bc4704d1b89ba925ece059 Расчетный анализ внутренней аэродинамики пневмосепарирующей машины.. (5)

Подстановка (5) в (4) дает возможность вычислять величину горизонтального смещения l для заданного значения высоты h падения частиц различной плотности, размера и значений скорости движения.

Для определения численных значений величины l зададимся диапазонами физических величин, входящих в выражения (4), (5) и характерных для устройства-прототипа:

- высота падения h: 800 – 1000мм;

- скорость сносящего потока Wc: 5 – 15м/с;

- диаметр сферической частицы : 2 – 5мм;

- плотность частицы ρч: 600 – 1400кг/м3.

Давление в потоке воздуха примем равным атмосферному. Для свободно падающих частиц примем, что их движение можно описывать как движение частицы, обтекаемой бесконечным потоком воздуха. Значение коэффициента аэродинамического сопротивления для обтекания бесконечным потоком, выбранных диапазонов скорости сносящего потока и диаметра частиц в соответствии с данными [6] примем равным Cx = 0,4.

На рис.2 приведена зависимость горизонтального смещения свободно падающих с высоты 800мм частиц в горизонтальном сносящем потоке от плотности и диаметра частицы (коэффициент k). Для выбранного диапазона параметров частиц k изменяется в пределах 0,053 – 0,306.

Маркерами на рис.2 выделены значения для частиц радиусом 2мм и плотностью 1120кг/м3, 1220кг/м3 и 1320кг/м3. Такие параметры частиц характерны, например, для зерна пшеницы, сепарируемой по плотности. Как следует из графика, разница горизонтального смещения указанных фракций не превышает 100мм при скорости сносящего потока Wcx=15м/с. Учитывая, что в процессе работы сепарирующей машины скорость сносящего потока может незначительно изменяться, практически осуществить отбор каждой отдельной фракции на длине 100мм представляется затруднительным.

Увеличить разницу горизонтального смещения различных фракций сепарируемого потока частиц позволяет увеличение времени пребывания частицы в сносящем потоке (см. выражение (4)). К увеличению времени t (см. выражение (5)), приводит наличие вертикальной составляющей Wcy сносящего потока и начальной скорости Wyo движения частицы в вертикальном направлении, направленными против действия силы тяжести. На рис.3 приведена зависимость времени падения частиц (данные и обозначения соответствуют рис.2). Выражение (5) и график на рис.3 демонстрируют нелинейную зависимость времени от Wcy сносящего потока. К существенному (более 10%) увеличению времени падения частицы приводит вертикальная составляющая Wcy>4м/с. Указанная нелинейность обуславливает нелинейный характер изменения величины горизонтального смещения, и, как следствие, нелинейный характер изменения разницы горизонтального смещения различных фракций.

image008_19_7b4c12a711ba09eec7d380d5417b8e2d Расчетный анализ внутренней аэродинамики пневмосепарирующей машины. image009_18_6a5abfe8199d82b4ca0602068f09a1c9 Расчетный анализ внутренней аэродинамики пневмосепарирующей машины.

Рис.2. Рис.3.

Увеличение времени падения частицы приводит к росту абсолютной величины горизонтального смещения, то есть к увеличению длины проточной части сепарирующей машины. Сократить длину позволяет наличие начальной скорости частицы Wxo, то есть частица должна быть вброшена навстречу сносящему потоку воздуха. На рис.4 представлены траектории движения частиц (обозначения на рис.2) при Wcx=10м/с, Wcy=8м/с, Wyo=1м/с (Wxo=0м/с – рис.4а, Wxo=4м/с – рис.4б). Наличие вертикальных скоростей увеличило разброс фракций до 300мм, а вбрасывание частиц против потока (рис.4б) уменьшило длину проточной части в 2раза.

image010_14_cabb75806bf38f5c3531af84090e301d Расчетный анализ внутренней аэродинамики пневмосепарирующей машины. image011_12_02c8810b3ec1136b6f0b789d1d60ee8d Расчетный анализ внутренней аэродинамики пневмосепарирующей машины.

a) б)

Рис.4.

Рис.5 представляет зависимость разницы горизонтального смещения различных фракций dl от скоростей Wcy и Wxo при Wcx=10м/с, Wyo=0м/с (на рисунке буквами обозначены разности фракций: а - 1120 и 1220 кг/м3, б – 1220 и 1320 кг/м3). Из графика рис.5 следует, что при отсутствии вертикальной составляющей сносящего потока Wcy увеличение скорости Wxo встречного вбрасывания частиц в сносящий поток приводит к росту разницы горизонтального смещения dl, с ростом Wcy влияние Wxo на dl падает. График демонстрирует также нелинейность влияния Wcy на величину dl: при Wcy=5м/с разница горизонтального смещения фракций растет в среднем в 2 раза, а при Wcy=8м/с – в 5 раз.

Таким образом, вертикальная составляющая скорости сносящего потока Wcy>5м/с, а также вбрасывание частиц в сносящий поток с вертикальной и горизонтальными составляющими, направленными вверх и против горизонтальной составляющей сносящего потока соответственно, позволяют существенно (до 5 раз) увеличить разницу горизонтального смещения сепарируемых фракций при длине проточной части, не превышающей 800-1000мм (высота проточной части 800мм).

image012_9_52e7a9eef13e3e24b31b57ebb692b97a Расчетный анализ внутренней аэродинамики пневмосепарирующей машины.

Рис.5.

Недостатком модели (1) – (2) является неучет взаимного действия всех частиц сепарируемого вороха в сносящем потоке, а также пульсационных составляющих всех скоростей. Для решения задачи сепарации вороха в турбулентном сносящем потоке нами использована конечно-элементная модель турбулентного взаимодействия в поле сил гравитации потока воздуха и вороха - потока, содержащего смесь воздуха с частицами различного размера и массы. Модель позволяет варьировать объемную долю частиц в ворохе, моделируя тем самым процесс сепарации различной массовой производительности. В качестве решателя конечно-элементной модели использована программа Fluent [7].

Конечно-элементная модель представляет собой двумерную прямоугольную сеточную область размером 800х800 мм, содержащую 160000 – 180000 четырехугольных элементов со стороной размером 2мм. В качестве частиц использованы сферические гранулы диаметром 4мм. Гранулы подавались во входное отверстие сверху массовым расходом 1,5 – 2,0кг/с. Сносящий поток поступал в расчетную область через границу слева. Для решения использована k-e турбулентная модель течения и эйлерова модель многофазного течения [7]. На рис.6 представлены результаты расчета траектории движения гранул (объемной доли фракции гранул в потоке) для случая свободного падения гранул массовым расходом 1,5 кг/с в поле горизонтального сносящего потока.

image013_29_fa134ffcfd40f924a343f710b664f360 Расчетный анализ внутренней аэродинамики пневмосепарирующей машины.

а) б)

а) – Wcx=10м/с; плотность гранул – 640кг/м3; б) – Wcx=8м/с; плотность гранул – 1320кг/м3; линиями обозначены контуры объемных долей фракции в потоке, точками – траектория, рассчитанная по модели (1) – (2).

Рис.6.

Конечно-элементная модель «строит» траекторию движения массы частиц ниже по сравнению с моделью (1) – (2), при этом данные расчетов по обеим моделям достаточно хорошо отслеживают друг друга. На рис.7 приведены результаты расчета моделей для случая падения гранул трех различных плотностей с равными объемными долями в ворохе общим расходом 1,5 кг/с. Расчет для смеси гранул не привел к существенным отличиям результатов модели (1) - (2) и конечно-элементной модели по сравнению с расчетом для гранул одной плотности.

image014_28_d89aad9c3ea7aaaafab1490360e4971b Расчетный анализ внутренней аэродинамики пневмосепарирующей машины.

image015_6_ff1f27eee0d715f03cca2a408bed325d Расчетный анализ внутренней аэродинамики пневмосепарирующей машины.

Wcx=12м/с; плотность гранул – 1120, 1220, 1320кг/м3; линиями обозначены контуры объемных долей фракции в потоке, точками – траектории, рассчитанные по модели (1) – (2).

Рис.7.

Wcx=8м/с; Wcy=4м/с; Wxo=2м/с; плотность гранул – 1000кг/м3; линии - контуры объемных долей фракции в потоке: 1 – массовый расход гранул 2,5кг/с; 2 – массовый расход гранул 1,25кг/с; 3 – массовый расход гранул 0,15кг/с; 4 - траектория, рассчитанная по модели (1) – (2).

Рис.8.

На рис.8 приведены результаты расчета траектории движения гранул в сносящем потоке, направленном под углом 30º к линии горизонта. Поток гранул подается по горизонтали навстречу сносящему потоку воздуха. На рис.8 видно существенное отличие траекторий, рассчитанных конечно-элементной моделью и моделью (1) – (2). В случае подачи гранул навстречу сносящему потоку траекторию движения потока гранул можно разделить на две зоны: зона встречного движения – от подающего сопла до места разворота потока, и зона спутного движения – от места разворота потока до плоскости падения гранул. В зоне встречного движения гранулы обтекаются сносящим потоком в «решетке» других гранул, что приводит к изменению отрывного режима обтекания (для случая бесконечного потока) на струйный режим обтекания. Результатом изменения режима обтекания является уменьшение коэффициента аэродинамического сопротивления [6] и, как следствие, увеличение горизонтальной протяженности зоны встречного движения. Эффект «решетки» ослабевает при уменьшении массового расхода гранул через сопло (см. кривые 1, 2 и 3 на рис.8). В зоне спутного движения гранул из-за существенного удаления гранул в потоке друг от друга режим обтекания становится близким к обтеканию бесконечным потоком, а траектория движения гранул удовлетворительно описывается моделью (1)–(2) при Сх=0,4. Уменьшение коэффициента аэродинамического сопротивления в зоне встречного движения потока гранул приводит к усилению эффекта разделения фракций – разница горизонтального смещения сепарируемых фракций по результатам расчета конечно-элементной модели больше, чем по модели (1)-(2) при постоянном значении Cx. На рис.9 приведено сравнение траекторий движения гранул различной плотности, рассчитанные по конечно-элементной модели и по модели (1)-(2).

image016_3_df61a2f83fbfd65ae68c6761d29940e3 Расчетный анализ внутренней аэродинамики пневмосепарирующей машины.

Wcx=8м/с; Wcy=4м/с; Wxo=2м/с; массовый расход гранул 1,25кг/с; 1 - плотность гранул 1100кг/м3; 2 - плотность гранул 1000кг/м3; линии - контуры объемных долей фракции в потоке; точки - траектории, рассчитанные по модели (1)–(2).

Рис.9.

Выводы.

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

2. Для решения задач эскизного проектирования проточной части и выбора диапазонов режимных параметров сепарирующей машины для различных видов зерна возможно применение модели (1) – (2).

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

Литература

1. Карпенко А. Н., Халанский В. М. Сельскохозяйственные машины. М: Колос, -1983. 495с.

2. Справочник конструктора сельскохозяйственных машин. Под ред. А. В.Красниченко. М: ГНТИ, 1961. 860с.

3. Электронный источник. Режим доступа: www. petkus. com.

4. Головков А. Н. Пневмосортировальные машины семейства ПСМ. // Центрально-Черноземная государственная зональная МИС. Эл. источник. Режим доступа: http://rmz. menzelinsk. ru/smi/4.html.

5. Электоронный источник. Режим доступа: http://www. agrotech-almaz. com/ production/ general. html.

6. Цаплин С. В., Романов А. Е., Болычев С. А., Давыденко С. В., Тютьмин Д. В. Гидродинамика и газовая динамика Лабораторный практикум // Самарский государственный университет, физический факультет, кафедра физики твердого тела и неравновесных систем. Межкафедральная лаборатория «Взаимодействие излучения с веществом», Издательство «Универс-групп», 2005.

7. Электронный источник. Режим доступа: http://www. fluent. com.

Завалий А. А.

 

Расчетный анализ внутренней аэродинамики пневмосепарирующей машины. - 4.0 out of 5 based on 1 vote