Краткий обзор пневматических устройств для транспортирования сыпучих смесей и процессов внутри них Текст научной статьи по специальности «Физика»

ТЕХНИЧЕСКИЕ НАУКИ

УДК 621.6.04

КРАТКИЙ ОБЗОР ПНЕВМАТИЧЕСКИХ УСТРОЙСТВ ДЛЯ ТРАНСПОРТИРОВАНИЯ СЫПУЧИХ СМЕСЕЙ И ПРОЦЕССОВ ВНУТРИ НИХ

Афанасьев А. И., Потапов В. Я., Костюк П. А., Макаров В. А.

В работе рассматриваются основы расчета пневматических устройств для транспортирования различных сыпучих смесей. Составлена математическая модель данного устройства. Математически обоснованы геометрические параметры устройства и его основные узлы. Приведена схема камерного питателя пневмотранспортной установки и описан принцип его действия. Данные математические зависимости могут использоваться для расчетов пневмотранспортных установок применимо к различным отраслям промышленности.

Ключевые слова: пневмотранспорт; математическая модель; сыпучая смесь; транспортирование; питатель.

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

Транспортные трубопроводы могут быть проложены во многих местах с незначительной занимаемой площадью пространства.

Существенным недостатком пневмотранспортных установок является сравнительно большой расход энергии на единицу транспортного материала. Однако в целом ряде случаев при применении пневматического транспорта общие затраты будут меньше, чем при других видах транспортирования. Кроме того, очень часто пневматический способ транспортирования является частью технологического процесса какого-либо производства, где он становится незаменимым (например, в пневмогазовых сушилках торфобрикетного производства или в пневмоуборочных машинах) [1].

Пневмотранспортная установка состоит из питающего устройства, трубопровода, отделителя и воздуходувной машины. Основным ее элементом является камерный питатель (рис. 1), содержащий приемную камеру 1 с шарнирно смонтированной на ней крышкой 2. В приемную камеру 1 помещена гибкая (мягкая) емкость (тара) 3 с сыпучим материалом в закрытом состоянии. Разгрузочный трубопровод 4 смонтирован на крышке 2. Нагнетательный трубопровод 5 содержит патрубок 6 с нагнетательным соплом 7 с возможностью свободного перемещения вдоль него. Лезвие ножевого устройства 8 направлено вниз для разрезания верхней части гибкой тары 3. Для отвода излишка воздуха по воздуховоду 9 дно 10 приемной камеры 1 снабжено патрубком 11. Для удаления гибкой тары 3 из камерного питателя приемная камера 1 снабжена диском 12 с юбкой 13. Установка диска 12 на определенном расстоянии от дна 10 камеры 1 обеспечивается винтовым устройством 14. Замкнутое пространство 15 приемной камеры 1 сообщается с разгрузочным трубопроводом 4 через вентиль 16 и с нагнетательным трубопроводом 5 посредством вентиля 17 воздуховода. Вентиль 18 транспортного трубопровода предназначен для перекрытия и подачи сжатого воздуха к нагнетательному соплу 7. Диск 12 имеет возможность свободного пере-

28

Известия Уральского государственного горного университета

мещения по вертикали до стопорного устройства 19. Диск 12 снабжен отверстием 20, перекрытым обратным клапаном 21 [2].

Благодаря разности давлений, создаваемых воздуходувной машиной, в трубопроводе пневмосистемы возникает движение воздуш-

Рис. 1. Схема камерного питателя пневмотранспортной установки

ного потока. Материал, предназначенный для транспортирования, с помощью питающего устройства вводится в пневмосистему, где он попадает в воздушный поток [1].

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

Часто в таких устройствах для создания псевдоожижения осуществляется продувка слоя частиц воздухом через пористую сетчатую перегородку. В результате прохода воздуха через слои материала, из-за разности давлений перед слоем Р1 и над ним Р2, объем материала насыщается воздухом и он начинает увеличиваться в объеме, оставаясь неподвижным, при этом частицы материала начинают отделяться друг от друга. Частицы не покидают образовавшийся «псевдоожиженный» слой, поскольку их вес уравновешивается скоростным напором воздуха, т. е. его скорость равна скорости витания частиц v В данном состоянии перепад давления АР =

Р1 - Р2 = const и не зависит от скорости потока воздуха v которая будет равна критической скорости «псевдоожижения» йкр.

Критическая скорость аэросмеси йкр, м/с:

% = (1)

где с - коэффициент, зависящий от крупности материала, склонности к слеживанию и наличию влаги, с = 0,25-0,4; - весовая концентрация материала; ; - плотность соответственно материала и воздуха, кг/м3; g - ускорение свободного падения, м/с2; d - диаметр пневмопровода, м.

Необходимая скорость v потока, м/с:

^ = (1,05-1,1Кр , (2)

Сопротивление P пневмопровода при движении аэросмеси, Па:

Рп = ( 1 + 0,4ц)Рв , (3)

где - потеря давления при движении чистого воздуха, Па.

Потеря давления ЛР Па в поднимающихся коротких участках пневмопровода

ЛРг = ( 1 + ^>Рв gH, (4)

где H - разница геодезических отметок, м (высота подъема).

Полное падение давления ЛР, Па в пневмотранспортной системе

ЛР=Рп + ЛРг+ ЛРд,„ + ЛРз, (5 )

где ЛР...= vпрв/2 - динамические потери

давления, Па; ЛРз - потери давления в загрузочном устройстве, обычно ЛРз = 1 Па. Расход воздуха Q м3/мин:

Qc = бмРм^РРвХ (6)

где Qм - расход материала, транспортируемого по пневмопроводу, м3/мин.

Необходимая подача воздухоподающей машины Q м3/мин:

бв = kQ0, (7)

где k = 1,1-1,15 - коэффициент, учитывающий потери воздуха в пневмопроводе [3].

Расход и скорость воздуха влияют только на степень возрастания объема материала.

№ 3(39), 2015

29

Транспортировка материала может начинаться при условии, когда vn > v и равна скорости уноса v [4].

Во многих работах [1, 3, 4] установлено, что скорость транспортирования материала зависит от скорости витания частиц.

Аналитически скорость витания частиц материала в потоке воздуха определяется из условия равенства силы сопротивления F, действующей на частицу, которая движется относительно потока воздуха и силы тяжести G.

При числах Рейнольдса, характеризующих турбулентное обтекание частиц, имеет место квадратичный закон сопротивлений, т. е.

2

F = cS -2^, (8)

2

где c - коэффициент сопротивления при обтекании частиц потоком; S - площадь сечения частицы, перпендикулярного направлению потока, м2; - плотность воздуха, кг-с2/м4.

Если тело имеет форму шара, то ^ nd п nd

G =-----р g и S =-----, тогда, подставляя

6 1 м с? я у

4

значение S в формулу (8) и приравнивая F и

G, получим формулу для определения скорости витания v , м/с:

вит’

v = v =

в вит

4 £Фм

ЗСрв

где Ym - объемный вес частиц материала.

Коэффициент сопротивления С зависит от многих факторов и определяется экспериментальным путем. Для частиц шаровой формы при Re < 1 коэффициент С может быть определен по уравнению

24 С=—.

R

На практике значения для частиц различных материалов определяют непосредственно экспериментальным путем в аэродинамической (витательной) трубе [1].

Объем камерного питателя может быть оценен по формуле

Q

кп

М.,

“Рм ’

где Gм - транспортная производительность установки с камерным питателем, м3/ч; t -среднее время для загрузки 1м3 объема питателя материалом; и - количество питателей, работающих в цикле загрузки; - плотность материала, кг/м3.

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

Движение частицы материала в потоке газа. На частицу, находящуюся в потоке газа, действуют силы давления (напора) и тангенциальные (касательные) силы. Силы давления действуют перпендикулярно отдельным элементарным площадкам поверхности частиц (рис. 2) и заданы зависимостью

dP = pdS, (9)

где dP - элементарная сила; p - давление газа в рассматриваемом месте; dS - элементарная площадка.

Рис. 2. Силы, действующие на частицы в потоке газа

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

dT = zdS,

где dT - элементарная касательная сила; т -

30

Известия Уральского государственного горного университета

касательное напряжение текущего газа в данном месте поверхности частицы.

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

Р = fs dp = Ss 'jdp* + dpt + dp?. (10)

f = Js dT = fs JdT? + dTy + dTif (11)

где P - равнодействующая сила давления; T - равнодействующая тангенциальных сил; S - поверхность частицы.

Конечное (равнодействующее) воздействие газа U = T + P, веса частицы бч и пассивных сил Nр (ударные силы) можно заменить равнодействующей R действующей в центре тяжести частицы с равнодействующим моментом M, действующим вокруг оси, проходящей через центр тяжести. Равнодействующая в центре тяжести частицы вызывает последовательное продвижение ее и равнодействующего момента вращения частицы вокруг оси, проходящей через ее центр тяжести. Отсюда видно, что относительное обтекание, заданное разностью векторов абсолютных скоростей текущего газа и движения частицы, изменяется в каждый момент, а этим изменяется как равнодействующая, так и равнодействующий момент. Равнодействующую U можно разделить на две составляющие: в направлении обтекания частицы и в направлении, перпендикулярном к направлению обтекания частицы. Составляющая в направлении обтекания частицы называется сопротивлением О, а сила, перпендикулярная к направлению обтекания частицы, называется подъемной силой W. Таким образом, движение частицы весьма сложно и его нельзя математически выразить в общем виде. Движение одной из них направлено по оси трубопровода, а две других взаимно перпен-

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

Движение частицы в неподвижном газе под действием собственного веса. Частица в неподвижной (спокойной) газовой среде изображена на рис. 3. Абсолютная скорость течения газа c равна нулю. Под влиянием веса Q4 частица начинает двигаться в направлении действия веса со скоростью v. Относительная скорость обтекания w = —v равна абсолютной скорости движения частицы. Для того чтобы можно было при этом изучать весь процесс как одноразмерный, мы пренебрежем воздействиями, которые могли бы вызвать движение перпендикулярно направлению обтекания, т. е. подъемными силами.

Тогда конечное движение частицы будет прямолинейным в направлении вектора веса и будет зависеть только от составляющей в направлении обтекания, т. е. от сопротивления частицы О. На частицу действует в направлении движения составляющая гравитационной силы (веса), а против направления движения - сопротивление среды и статическая подъемная сила W обусловленная законом Архимеда. Уравнение движения будет иметь форму

бч - Wc - о=Q а, (12)

g

т. е. ускорение частицы а в каждом моменте

№ 3(39), 2015

31

зависит от мгновенной величины отдельных сил. Рассмотрим эти величины. Вес частицы постоянен и его можно определить из уравнения

0ч = ^чУ м ,

т. е. он зависит от объема частицы и от удельного веса материала частицы. Статическая подъемная сила частицы задана зависимостью

W = бчУр.

Эта подъемная сила зависит от объема частицы и удельного веса среды уг, в которой находится частица. Величина статической подъемной силы в газовой среде более чем в 1000 раз меньше веса частицы, поскольку отношение удельных весов материала частиц Y

и газовой среды — >1000. Поэтому для пра-

Уг

ктических рассуждений о движении частицы в газовой среде можно при пневматическом транспортировании пренебречь действием статической подъемной силы, не нанося этим никакого ущерба точности. Необходимо обратить внимание на то, что в отношении других сред (например, воды при гидравлическом транспортировании) уже нельзя пренебречь воздействием статической подъемной силы, так как удельный вес воды почти в 1000 раз больше удельных весов газов. Этот факт позволяет упростить уравнение движения до выражения

0ч = о = 0 a. (13)

g

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

2

P = P0- —^-tSz, (14)

где P - давление в данной точке поверхности частиц; Р0 - давление в спокойной среде; с -скорость обтекания частицы в данном месте поверхности частицы; - потеря, считая с места, где измерено давление P0.

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

т = 3 —, (15)

dy

где т - тангенциальное напряжение внутри среды; dc/dy - градиент скорости в направлении, перпендикулярном направлению скорости; 3 - коэффициент динамической вязкости в кг-с/м2.

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

Распределение тангенциальных сил вдоль поверхности частицы образует сопротивление трения тела. Оба вида сопротивления воз-

32

Известия Уральского государственного горного университета

действуют друг на друга, и обычно их нельзя разделять. Имеются определенные формы тел, у которых преобладает один из этих двух видов сопротивлений. У очень тонкой пластинки, расположенной перпендикулярно направлению обтекания, преобладает сопротивление формы, и наоборот, если эта пластинка расположена в направлении обтекания, то в ней преобладает сопротивление трения.

В общем виде сопротивление частицы задано зависимостью

O=^ч^уг

2 g

(16)

где £, - коэффициент сопротивления; F4 -аэродинамическое сечение частицы; ю - относительная скорость обтекания; уг - удельный вес газовой среды; g - ускорение силы тяжести (9,81 м/с2).

Аэродинамическое сечение F - это прямоугольная проекция частицы в направлении относительного обтекания. У частиц сферической формы аэродинамическим сечением является круг, диаметр которого равен диаметру частицы, а у частиц кубической формы, стенка которых перпендикулярна направлению относительного обтекания, аэродинамическим сечением является квадрат со стороной, равной грани куба. У частицы кубической формы, диагональ которой расположена в направлении обтекания, аэродинамическим сечением является правильный шестиугольник со стороной 0,8165е, где е -грань куба.

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

юе

(17)

юе

Re =—,

где е - характеристический размер частицы;

v - кинематическая вязкость среды, заданная отношением $/Q где $ - динамическая вязкость, а - удельная масса газовой среды.

У газов динамическая вязкость возрастает с температурой. Значения динамической вязкости для воздуха (чаще всего используемого в качестве транспортирующего газа) в зависимости от температуры приведены на диаграмме (рис. 4). Влияние давления на динамическую вязкость незначительно, и практически им можно пренебречь. Чаще всего применяют к воздуху при 20 °С показатель динамической вязкости 1,84-10-6 кг-с/м2 и кинематической вязкости 15 м2/с.

106'

Рис. 4. Зависимость динамической вязкости воздуха от температуры [9]

В основном различают две группы форм частиц, одна из которых имеет коэффициенты сопротивления, не зависящие от числа Рейнольдса, тогда как вторая группа имеет коэффициенты сопротивления, зависящие от этого числа. У первой группы отрыв граничного слоя обусловлен однозначно влиянием ребра обтекания. Формы тел или частиц отличаются в основном тем, что вплоть до грани обтекания сечение тела либо возрастает, либо остается постоянным. На грани обтекания сечение тела резко изменяется. Коэффициенты сопротивления некоторых форм тел первой группы приведены в табл. 1. Вторая группа включает в округленные тела, не имеющие острых ребер обтекания, и у которых сечение в направлении течения уменьшается не сразу. Здесь место отрыва граничного слоя обусловлено условиями обтекания, и поэтому коэффициент сопротивления зависит от числа

Рейнольдса. Наиболее тщательно исследован коэффициент сопротивления шара в зависи-

№ 3(39), 2015

33

мости от числа Рейнольдса. Этот коэффициент сопротивления приведен на диаграмме (рис. 5) и в табл. 2. Для некоторых других форм эта зависимость указана ориентировочно на диаграмме (рис. 6).

35 ------------------------------------

Рис. 5. Зависимость £ от Re у шара; зависимость £Re2 и Re/£ от Re [10]

Величина сопротивления для обоих этих случаев задана следующими зависимостями:

2

O = (18)

2g

где £ = const

или

£=f

' юеуг ^

V % у

(19)

Зная вязкость и удельный вес газа, скорость обтекания частиц, их форму и размер, можно по приведенным выше зависимостям определить величину сопротивления. Окончательным уравнением движения после подстановок будет

Q, -SF,T-rr = — “■ (20)

2g g

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

v2 Q

Йч - £F, — Y г = Q a. (2 1 )

2g g

Таблица 1

Коэффициенты сопротивления некоторых форм тел первой группы

£ Форма тела

2,01 a/b=x>

1,40 a/b=18

1,33 Полый полушар

1,29 a/b=10

1,22 Диафрагма (прямоугольная плита для различных a/b)

1,19 a/b=4

1,17 Сплошной полушар

1,15 a/b=2

1,11 Круглая плита

1,10 a/b=1

0,99 1/d=7

0,91 1/d=1

0,87 1/d=4

0,85 1/d=2

0,51 Угол при вершине 60°

0,40 Сплошной полушар

0,34 Полый полушар

0,34 Угол при вершине 30°

Однако по мере возрастания скорости непрерывно возрастает второй член левой части уравнения (21) a < g и ускорение уменьшается до достижения равновесия, когда a = 0, вытекающее из уравнения

v2

Йч- F—Уг=о. (22)

2£

10

0 5 10 15 20 25

-2 ------------------■----------------—---------

Рис. 6. Зависимость £ от Re для частиц, имеющих форму валика и эллипсоида [11]

Если бы на частицы не действовало сопротивление среды (£ = 0), то из этого уравнения вытекала бы известная зависимость a = g, подразумевающая равномерно ускоренное движение частицы. При сопротивлении среды (£ Ф 0) в начальный момент, когда частица находится еще в покое (v = 0), сопротивление ее также равно нулю, а ускорение a = g.

Поскольку все показатели, кроме v, в этом уравнении постоянны, получается, что при равновесном состоянии частица достигает некоторой конечной постоянной скорости, которую мы называем скоростью падения v Величина ее равна_____________

v

п

2gQ4 ^;Ум £ F Уг

(23)

34

Известия Уральского государственного горного университета

Для частиц с постоянным коэффициентом 2 определить эту скорость просто, поскольку все величины в уравнении заданы однозначно. Труднее определить скорость падения тех частиц, у которых коэффициент 2 зависит от числа Рейнольдса, т. е. и от vH . Ориентировочно определить скорость падения можно по диаграмме на рис. 5 или 6, узнав коэффициент 2, и по уравнению (23) проверить vH. Если значение vH окажется больше, чем выбранное нами, то нужно заново выбрать так, чтобы коэффициент был большим, чем у ранее выбранной величины; если же vu окажется меньше, то нужно поступить наоборот.

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

Характеристическим размером является диаметр шара d. Вес шаровой частицы равен

nd3

Аэродинамическое сечение частицы

F = — d2 F = 4 d .

Подставим эти зависимости в уравнения (17) и (23):

v d

Re = -

v

4g A

(24)

(25)

“ 3 2Yr

Исключив vH из обоих уравнений, получим после преобразования

2Re2 = -

.А Аа=А А

3 v2 уг 3g v Исключив d, получим

2 YмYг.

Re

2 :

'•y 3

A VlIi 4g v YM

3 2

v Y

п t г

4g2 a Y„

(26)

(27)

v

Далее используем зависимости

= = M Re и d=^ =®g^. (28)

d Yг vH Y^h

Таблица 2

Коэффициенты сопротивления некоторых форм тел

Re 2 2Re2 Re / 2

110-2 2,4-103 2,4-104 4,17-10'6

1-10-1 2,4Д02 2,4^10° 4,17Д0'4

2 1,2 4,8 1,6710-3

3 8,0101 7,2 3,75

5 4,95 1,24101 1,0110-2

7 3,65 1,79 1,92

1-100 2,65 2,65 3,77

2 1,44 5,76 1,3910-1

3 1,04 9,37 2,88

5 6,9100 1,73-102 7,25

7 5,4 2,65 1,296Д00

1-101 4,1 4,10 2,44

2 2,55 1,02103 7,84

3 2,0 1,8 1,50101

5 1,5 3,75 3,33

7 1,27 6,23 5,55

1Д02 1,07 1,07-104 9,35

2 7,7Д04 3,08 2,59102

3 6,5 5,85 4,61

5 5,5 1,38105 9,09

7 5,0 2,45 1,40103

1103 4,6 4,60 2,17

2 4,2 1,68106 4,76

3 4,0 3,60 7,50

5 3,85 9,60 1,299 104

7 3,9 1,91107 1,785

1Д04 4,05 4,05 2,469

2 4,5 1,8108 4,444

3 4,7 4,23 6,383

5 4,9 1,23109 1,020 105

7 5,0 2,45 1,400

1105 4,8 4,8 2,083

2 4,2 1,681010 4,762

3 2,0 1,80 1,500 106

4 8,410-2 1,34 4,761

6 1,010-1 3,6 6,000

1Д06 1,3 1,3-10и 7,692

3 2,0 1,81012 1,500 107

Кроме того, нам необходимо знать зависимость значений 2Re2 и Re / 2 от числа Рей-

№ 3(39), 2015

35

нольдса. Эти значения можно вывести из диаграммы £ = f(Re), также они приведены на рис. 5 и в табл. 2.

Порядок расчета следующий. Если известен диаметр шаровой частицы d, то из уравнения (26) мы определяем значение £Re2. Из диаграммы на рис. 5 или из табл. 1 определяем число Рейнольдса, а из первой зависимости (28) - vn. Наоборот, если известна скорость падения шаровой частицы, то определяем из выражения (27) значение , из диаграммы на рис. 5 или из табл. 1 - значение числа Рейнольдса, а по второй зависимости (28) определяем d, т. е. диаметр шаровой частицы. Таким образом, графоаналитическим методом можно определить скорость падения частиц округлой формы.

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

£

Re”0

(29)

Наиболее известно уравнение Стокса, где K0= 24, а ”0 =1, т. е. £ = 24/Re. Это выражение действительно для чисел Рейнольдса до 0,2, но часто его используют и в диапазоне до 2.

Подставив коэффициент в уравнение (25), получим

v

п

(30)

Однако правильность этой зависимости ограничена условием

VI 73 Re -Ш2 g (31)

УмУг

или v- - 3 ^е2Ум %2 18y2 , (32)

где Re имеет значение 0,2 или соответственно 2.

До сих пор мы рассматривали установившееся движение частицы в спокойной среде. Рассмотрим теперь детально движение частицы в спокойной среде при неравновесном состоянии. При постоянном коэффициенте сопротивления £ и начальных условиях x = 0, v = 0 зависимость пути x от скорости v частицы задано зависимостью

x = -

2

5“

2 v

v vH У

( 3 3)

или же

x=-

v

2 g

mv

(34)

Зависимость значений m1 изображена на диаграмме (рис. 7). Из этой диаграммы видно, что теоретически частица достигает скорости падения лишь на бесконечно длинном пути и что 992% скорости она приобретает на пути x = 0,2 V“.. С помощью диаграммы можно определить скорость движения частицы в любой точке ее пути.

Рис. 7. Диаграмма m1 = -fly/v^

До сих пор мы исходили из постоянного коэффициента сопротивления для выведения кинематических величин при движении частицы. Это справедливо только для частиц, имеющих острую грань обтекания. Для округленных частиц, у которых коэффициент £ зависит от числа Рейнольдса, интегрирование дифференциальных уравнений затруднительно, так как есть еще одна переменная, а именно коэффициент £, зависящий от числа Рейнольдса. Эта зависимость выражена в виде таблицы или графика (рис. 5 и табл. 2). Аналитические выражения, выведенные различными авторами, относящиеся к определённым областям чисел Рейнольдса, сложны, и

36

Известия Уральского государственного горного университета

решение дифференциальных уравнений движения в этом случае затруднительно. Очень простое решение возможно при использовании зависимости Стокса для коэффициента сопротивления £ = шаровой частицы.

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

К=,

const

------е —

£ Yr

(35)

Для всех подобных частиц граненой формы выражение постоянно, поэтому

vH = const

е ь.

Yr

(36)

Скорость падения возрастает с корнем удельного веса материала частиц и с размером частиц. Влияние среды сказывается только ее удельным весом Yr. У одной и той же частицы произведение v^YT = const. Скорость падения частицы обратно пропорциональна квадратному корню из удельного веса среды. Если состояние среды задано,

Yг

Р

RT ’

(37)

то уравнением состояния является скорость падения частицы, обратно пропорциональная квадратному корню из давления и прямо пропорциональна квадратному корню из температуры среды. Скорость падения частицы не зависит в этом случае от вязкости среды. Для определения скорости падения округленных частиц, где коэффициент £ = f(Re), справедливо уравнение

vH = f (Re) const e — = V ^

г

f v у е п t г

1 )

const e—.

Yr

(38)

В этом случае зависимость скорости падения от параметров среды более сложная.

Чтобы показать хотя бы в общих чертах зависимость скорости падения сферической частицы от параметров, заменим зависимость , указанную в диаграмме на рис. 5, для отдельных участков зависимостью (29), т. е.

£

А

Re”

о

Логарифмированием этой зависимости получим lg = lg K0— lg Re. Поскольку зависимость £ от Re представлена в логарифмических координатах, то выражение представляет в этом случае прямую, точка пересечения которой с осью ординат (для числа Рейнольдса Re = 1) является значением К0, а ”0 - тангенс угла ф, образуемого прямой с отрицательной осью абсцисс (см. рис. 8). На диаграмме одновременно показана зависимость значений и от числа Рейнольдса для отдельных точек кривой £ = f (Re).

Для определенного диапазона чисел Рейнольдса можно выбрать значения п0 и К0 в качестве средних. Тогда для скорости падения шаровой частицы будет справедливо выражение

v

п

-A Re”o d Ь-.

3К0 Yг

После подстановки вместо Re лучим

v

2-п,

п

4

3

g

d1 + По _Y^_ К09”0 у^”0 '

(39)

v dy

п t г

9

по-

(40)

Для области действия закона Стокса

1 d2

(”о =1, Ко = 24) справедливо vn = —— Ум.

18 v

№ 3(39), 2015

37

В этой области скорость падения пропорциональна квадрату диаметра частицы, удельному весу материала частицы и обратно пропорциональна вязкости среды. Поскольку вязкость газа возрастает с температурой, то одновременно скорость падения обратно пропорциональна температуре среды. Для области 2 = const («0 = 0, K0 = 0,47)

4 d ум

v = - g-------—

п V 3 0,47 уг

лютного давления среды. Здесь скорость падения не зависит от вязкости среды. Понятно, что в отношении зависимостей между обеими областями имеется большое различие.

Для любого n0 справедливо

(

v=

П

2 - «0

f gл

V 3K0 J

1 + «0 _ X d 2 - «0 Д «

Yr J 1

- 2y 2 - «0 I м

«

0

X

а это выражение тождественно выражению для 2 = const.

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

Для переходной области между обоими крайними случаями заменим зависимость прямой, взятой по диаграмме на рис. 8.

Определим значения К0 и «0, а изменение скорости падения определим из выражения

«2

1 - «0

( ^ 2 - «0 f

v

Y г2 J

1 + «0 d2 '2 - «

V di J

Г 2 V^2 J

f Ym2

. Yмl J

2 - «0

где v , задано уравнением (25).

n

X

БИБЛИОГРАФИЧЕСКИМ СПИСОК

1. Урбан Я. Пневматический транспорт / под ред. Л. М. Шведова. М.: Машиностроение, 1967. 256 с.

2. Энергосберегающее оборудование для транспортировки сыпучих материалов: Исследование, разработка, производство / под ред. И. Г. Южакова. Екатеринбург: УГТУ-УПИ, 2007. 317 с.

3. Транспортировка и складирование порошкообразных строительных материалов / И. П. Малевич [и др.]. М.: Стройиздат, 1984. 184 с.

4. Давыдов С. Я. Новое энергосберегающее печное, транспортное и складское оборудование в производстве цемента: учеб. пособие. Екатеринбург: УГТУ-УПИ, 2004. 170 с.

5. Кузнецов Ю. М. Пневмотранспорт: теория и практика. Екатеринбург: УрО РАН, 2005. 61 с.

6. Камерный питатель нагнетательной пневмотранспортной установки / С. Я. Давыдов [и др.]. Пат. № 1437320 RU. Б.И., 1988. №42.

7. Камерный питатель нагнетательной пневмотранспортной установки / С. Я. Давыдов [и др.]. Пат. №1676970 RU. Б.И., 1997. № 19.

8. Камерный питатель нагнетательной пневмотранспортной установки / С. Я. Давыдов [и др.]. Пат. № 2083458 RU. Б.И., 1997. № 19.

9. Идельчик И. Е. Справочник по гидравлическим сопротивлениям. М.: Рипол Классик, 2013. 468 с.

10. Успенский В. А. Пневматический транспорт. Свердловск: Металлургиздат, 1959. 232 с.

11. Справочник авиационного инженера / В. Г. Александров [и др.]. М.: Изд-во «Транспорт», 1973. 354 с.

Поступила в редакцию 2 сентября 2015 г.

Афанасьев Анатолий Ильич - доктор технических наук, профессор, заведующий кафедрой горнопромышленного транспорта. 620144, г. Екатеринбург, ул. Куйбышева, 30, Уральский государственный горный университет.

Потапов Валентин Яковлевич - доктор технических наук, профессор кафедры горной механики. 620144, г. Екатеринбург, ул. Куйбышева, 30, Уральский государственный горный университет.

Костюк Петр Андреевич - аспирант кафедры горной механики. 620144, г. Екатеринбург, ул. Куйбышева, 30, Уральский государственный горный университет. E-mail: [email protected] Макаров Владимир Александрович - генеральный директор. 620028, г. Екатеринбург, ул. Татищева, 94, ЗАО Центр диагностики и экспертизы «Цветметналадка».

38

Известия Уральского государственного горного университета