Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

глава3

.pdf
Скачиваний:
6
Добавлен:
03.06.2015
Размер:
353.74 Кб
Скачать

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

v = Rб(μ0χтH grad H + g cosϕм),

где μ0 – магнитная постоянная; χт – удельная магнитная восприимчивость тела; H – напряженность магнитного поля; grad H – градиент напряженности магнитного поля.

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

 

 

v2

R − μ

χ

т

H grad H

ϕм

 

б

0

 

 

= arccos

 

 

g

 

,

 

 

 

 

 

 

 

и немагнитных частиц:

 

 

 

 

ϕ

м

= arccosv2 (gR

)1.

 

 

 

 

 

б

 

 

 

gcosφ

V2

 

Rб

φgsinφ

fм g

Rб

Рис. 3.1. Силы, действующие на частицу руды в сепараторах с верхним питанием (трение не учитывается)

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

dmin =

36μhu

 

,

L( fм g(δ −

) δ)

 

 

где h – расстояние от днища лотка до барабана; L – длина зоны разделения; u – скорость подачи материала.

Необходимую длину зоны разделения для извлечения магнитных частиц с

размером большим d определяют по формуле

 

 

L =

18μhu ± u

2h(16μ2h + d 4δ2 )[ fм g(δ −

) δ]

.

 

d 2δ[ fм g(δ − ) δ]

 

 

 

 

 

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

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

11

Сепаратор включает в себя плоскость с трамплином и вращающийся ба-

рабан. Материал, попадая в сепаратор, поступает в верхнюю часть наклонной

плоскости, по которой он начинает движение к нижнему краю. Вследствие раз-

ности коэффициентов трения компонентов происходит расслоение материала, а

на сходе с плоскости образуется веер разделения. Далее частицы попадают на

барабан, где, ввиду различий по характеру отскока частиц и коэффициентам

трения, возможно их различное дальнейшее движение.

 

 

 

 

Составляя модель, выделим следующие фазы движения частиц, описыва-

емые принципиально различными уравнениями движения по наклонной плос-

кости: движение по трамплину, движение после схода с полки (полет), соуда-

рение частицы с движущимся барабаном, движение после соударения. Соуда-

рение и полет после первого контакта могут повторяться многократно. Случаев

удара о барабан может не быть, т. е. не все фазы обязательны.

 

 

Рассмотрим математические модели движения частиц в различных фазах.

Схема, принятая для описания движения частицы обогащаемого материа-

ла в зоне сепарации, приведена на рис. 3.2. Частица начинает движение без

начальной скорости по наклонной плоскости 1 и движется ускоренно из точки 0

до точки А, после чего плавно выходит на криволинейный участок АВ, который

можно для простоты описания считать дугой окружности 2. Направление дей-

ствующих сил при этом изменяется, и движение замедляется; изменяется также

по величине и направлению скорость частицы. Далее частица входит в зону

 

 

 

свободного полета, где движется

0

 

 

под влиянием силы тяжести и силы

β

 

 

сопротивления воздуха, зависящей

1

 

 

A

B

от относительной скорости части-

L

цы. В зависимости от траектории

 

 

2

 

 

движения

возможно

возникнове-

 

 

 

 

 

 

ние двух ситуаций: либо частица

 

 

 

пролетает мимо барабана 3 до вы-

 

 

 

хода из зоны сепарации (т. е. до

 

 

3

соприкосновения с горизонтальной

 

 

поверхностью 4), либо, ударив-

 

 

 

4

 

 

шись о барабан, меняет величину и

Рис. 3.2. Схема движения частицы в бара-

направление скорости и с новыми

начальными

условиями

вновь

банно-полочном фрикционном сепараторе:

начинает

совершать

свободный

1 – наклонная плоскость; 2 – криволинейный

полет. При этом возможно не-

участок; 3 – барабан; 4 – горизонтальная по-

верхность

 

 

сколько соударений частицы с ба-

 

 

 

рабаном, но в конечном итоге про-

цесс обязательно завершится падением на горизонтальную плоскость. Эффек-

тивность разделения будет тем выше, чем больше расстояние между точками

падения на горизонтальной плоскости частиц с различными физическими свой-

ствами. Блок-схема, описывающая процесс, приведена на рис. 3.3. Выполним

декомпозицию модели в соответствии с выделенными фазами движения.

 

12

Вывод координат частицы
конец Рис. 3.3. Блок-схема исследования процесса сепарации
Да
Частица выходит из зоны сепарации?
Соотношение при ударе
Есть контакт с барабаном?
Да
Нет
начало
Ввод варьируемых параметров
Движение по наклонной плоскости
Движение по дуге окружности
Свободный полет частицы

Математическое моделирование этого процесса возможно, если будут составлены уравнения движения частицы на каждом этапе движения. При этом конечные условия каждого предыдущего этапа являют-

ся начальными условиями каждого последующего.

Таким образом, для математического моделиро-

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

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

3.3. Модели процессов разделения

Основные этапы получения уравнений сепарации и сепарационных характеристик

Единый подход к моделированию различных обогатительных процессов предполагает использование однотипных этапов, единого математического аппарата при моделировании различных сепарационных процессов. В соответствии с теорией сепарационных процессов любой продукт обогащения характеризуется двумя видами функций: w(ξ) и β(ξ). Любой процесс обогащения про-

Единый подход к моделированию различных обогатительных процессов с позиции общих закономерностей движения минераль-

ных частиц в рабочих зонах обогатительных аппаратов разработан О. Н. Тихоновым. Данная глава составлена по материалам его работ.

13

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

Любая точка пространства характеризуется функцией фракционного состояния w(ξ, х, у, z, t) и усредненной скоростью движения частиц элементарных функций в точке. Для окрестности любой точки справедлив локальный закон сохранения любой элементарной фракции, выражение которого является первым уравнением сепарации. Второе уравнение сепарации вытекает из баланса статистически усредненных сил, действующих на частицы любой элементарной фракции. Совместное решение уравнений позволяет определить функцию фракционного состояния для любой точки по известным начальным и граничным условиям. Усреднение функции фракционного состояния для стационарного режима по зонам концентрата и хвостов дает средний фракционный состав продуктов. По известному фракционному состоянию в зонах разгрузки и питания wк(ξ) и wхв(ξ), wисх(ξ) могут быть аналитически найдены сепарационные характеристики аппарата по концентрату и хвостам, представляющие собой зависимости вероятности извлечения элементарных фракций материала в продукты от значения физического свойства элементарных фракций и характеризующие степень совершенства аппарата.

Сепарационная характеристика, например, для концентрата, вычисляется

по формуле:

Qкwк(ξ)

 

 

Qкwк(ξ)

 

 

εк(ξ) =

=

,

Qисхwисх(ξ)

Qкwк(ξ)+ Qхwх(ξ)

 

 

 

поскольку в соответствии с уравнением баланса

Qисхwисх(ξ) = Qкwк(ξ)+ Qхwх(ξ),

Qисх, Qк, Qхв – производительность аппарата соответственно по исходному питанию, концентрату и хвостам;

wк(ξ) = Vк1∫ ∫∫ w(ξ, x, y, z)dV ;

Vк

wх(ξ) = Vх1∫ ∫∫ w(ξ, x, y, z)dV.

Vх

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

Составление уравнений сепарации

Как уже отмечалось выше, чтобы полностью охарактеризовать сепарационный процесс, достаточно составить уравнения сепарации и решить их относительно w(ξ, х, у, z) для любой точки. Затем, усредняя интегрированием по зо-

14

нам исходного питания, концентрата и хвостов, получить характеристики

фракционного состояния в зонах wисх(ξ, х, у, z), wк(ξ, х, у, z), wхв(ξ, х, у, z) с учетом начальных и граничных условий при допущении стационарности.

Одно из уравнений сепарации – уравнение типа закона сохранения, второе – баланса сил.

Большинство аппаратов можно представить одномерной моделью, выделив основное направление расслоения частиц под действием сил х, и, следова-

тельно, записать уравнение типа закона сохранения в виде:

 

(mw)

= −

(mwvx )

 

+ W.

 

t

 

x

 

 

 

 

Второе уравнение будет видоизменяться для конкретного аппарата и условий разделения в нем.

Заметим, что все силы отнесены к единице объема частицы с плотностью ρ. Составим уравнение сепарации на примере аппаратов с естественной разде-

ляющей средой переменной плотности (рис. 3.4), где действуют следующие

силы: гравитационная

F грав = ρg ;

среднестатистическая

гравитаци-

онно-архимедова

 

F г-а = − g

 

(x,t), где

 

 

– средняя

ρ

ρ

плотность среды постели, изменяющаяся по пространству x и c и с течением времени t; сопротив-

ления F сопр = −αvx и градиентная

F град = −kw1 wx . Здесь W=0.

t=w/vтр

x

 

xверх

Fграв Fсопр Fград

 

хвосты

Питание

xр

h

 

wисх(ρ)

Fг-а или Fарх

 

концентрат

 

 

 

 

 

 

 

xниж

Рис. 3.4 Схема, используемая при составлении уравнений сепарации для аппаратов с естественной разделяющей средой переменной плотности

Для этого случая уравнения сепарации имеют вид:

gρg

 

(x,t)αvx kw1 w

= 0;

ρ

 

(mw)

= −

(mwvx )

. x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

x

 

 

 

 

При исключении vx(ρ, x, t) при m=const получим уравнение сепарации с

одной неизвестной функцией w(ρ, x, t):

 

 

1 [w(ρ

 

)]

w

 

 

 

2w

 

 

 

 

= D

gα

ρ

 

 

x2

 

 

 

 

;

t

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

ρmax

ρ = ρw(ρ, x,t)dρ,

ρmin

где D=kα-1 – коэффициент макродиффузии, м2/с. Время сепарации t=y/vтр соответствует продвижению материала от загрузки к разгрузке.

15

Решение уравнений сепарации и получение сепарационных характеристик

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

Чтобы определить фракционный состав материалов в зоне разгрузки концентрата и хвостов, требуется ввести граничные условия. В зоне разгрузки материал полностью расслоился и поэтому vx=0, t=tк=wраз/vтр, а функция фракционного состава может рассматриваться как предельная wпред, не зависящая от t

(в зоне разгрузки достигается стационарный режим,

при

котором

wпред( ξ, x ) = w ( ξ, х,t)/t→∞). Для стационарного режима в

зоне

разгрузки

∂w/∂t=0. С учетом граничных условий основное уравнение сепарации решается

относительно w пред( ξ, x ) .

Общее решение w пред( ξ, x ) превращается в конкретное частное с учетом начальных условий, зависящих от фракционного исходного питания:

w ( ξ, x , t 0 )= w исх( ξ) .

Для нахождения w пред( ξ, x 0 ) через заданную функцию w исх( ξ) надо воспользоваться интегральным законом сохранения:

( ) ( )1 xверх ( )

wисх ξ = xверх xниж wпред ξ, x dx,

xниж

где xверх, xниж – координаты верхней и нижней частей слоя разделения. Сепарационную характеристику аппарата εк( ξ) находят по известным

wисх(ξ) и wпред(ξ).

 

 

 

 

 

xниж

 

εк(ξ)=

Qкwк(ξ)

 

wпред dx

 

=

xгр

,

Q

w

(ξ)

xниж

 

исх

исх

 

wпред dx

 

 

 

 

 

 

 

xверх

где Qк/Qисх= (xниж xгр)/(xниж xверх); xгр – координата положения отсекателя концентрата.

Здесь средний фракционный состав концентрата получается усреднени-

ем:

xниж

 

wк(ξ) = (xниж xгр)1 wпред(ξ, x)dx.

(3.1)

xгр

Окончательный вид сепарационной характеристики получают, подстав-

ляя найденное ранее w пред( ξ, x ) в уравнение для εк(ξ).

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

16

Основные уравнения баланса сил и закон сохранения пишутся следующим образом в предположении, что расслоение идет в одномерном простран-

стве:

 

 

 

 

(3.2)

 

Fграв + Fарх + Fсопр + Fград = 0;

 

(mw)

= −

(mwvx )

.

(3.3)

 

 

 

 

t

 

 

x

 

Запишем уравнение (3.2) с использованием таблицы 2.3:

gρgρср αvx k(mw)1grad(mw)= 0.

(3.4)

Если допустить постоянную концентрацию твердого в объеме m=const, то

vx из уравнения (3.4) запишется так:

 

vx = gα-1

(ρρср)k(wα)1grad w.

(3.5)

Подставим (3.5) в (3.3) и запишем одно уравнение сепарации:

w

k

t

=

 

 

α

grad w

g

[(ρ- ρср)w]

 

 

 

 

 

 

;

(3.6)

x

 

x

 

α

 

w= w(ρ, x, t).

Сучетом D=kα-1 и в предположении, что плотность разделения постоянна по объему и во времени и равна средней плотности кусков, уравнение (3.6) перепишем в виде:

w

= D

2w

g(ρρср)α1 w

(3.7)

t

 

x2

 

t

 

или

 

2w

 

 

 

w

= D

A

w

(3.8)

t

x2

t

 

 

 

где D – коэффициент диффузии; А = g(ρ-ρср-1 – коэффициент сноса.

Для любой узкой фракции уравнение (3.8) относится к типу уравнений Фоккера-Планка. Для стационарного режима (∂w/∂t=0) решение предельной задачи при условии полного расслоения в коне аппарата имеет вид:

w

(ρ, x)

 

 

 

= w

(ρ, x

)exp[AD1(x x )].

(3.9)

 

 

 

пред

 

t→∞

пред

0

0

 

 

 

 

 

 

 

 

В соответствии с

интегральным

законом сохранения,

если известно

w ( ρ, x 0 ) для какой-то точки x0 между верхней границей слоя расслоения с координатой xверх и нижней границей слоя с координатой xниж:

wпред(ρ, x0 ) =

 

AD1hw

 

(ρ)exp (AD1x

)

 

 

 

 

 

исх

 

0

 

 

,

(3.10)

exp (AD1x

 

 

)exp (AD1x

 

)

 

 

 

 

 

 

 

 

ниж

верх

 

 

 

 

где h= xнижxверх – толщина слоя, м. Подставляя (3.9) в (3.8),имеем:

 

wпред(ρ, x) =

 

AD1hw

 

(ρ)exp (AD1x)

 

 

 

 

 

 

исх

 

 

.

(3.11)

 

exp (AD1x

 

 

)exp (AD1x

 

 

 

 

 

 

 

)

 

 

 

 

ниж

верх

 

 

 

 

Используя (3.1) и учитывая:

17

 

 

 

xниж

 

wк(ρ) = (xниж xгр)1 wпред(ρ, x)dx;

(3.12)

 

 

 

xгр

 

 

 

 

xниж

 

wисх(ρ) = (xниж xгр)1 wпред(ρ, x)dx,

(3.13)

 

 

 

xгр

 

а также (3.14)

 

 

 

 

 

Qк

=

xгр xверх

,

(3.14)

 

 

 

 

Qисх

 

xверх xниж

 

имеем выражение для сепарационной характеристики:

 

xниж

 

εк(ρ)==

wпред dx

 

xгр

.

xниж

 

wпред dx

 

 

xверх

 

Подставляя (3.11) в (3.15) и интегрируя, получаем:

εк(ρ) =

exp (AD1x )exp (AD1x

)

 

гр

верх

.

exp (AD1x

 

 

)exp (AD1x

 

)

 

ниж

верх

 

 

(3.15)

(3.16)

Если h – высота слоя разделения, а за точку отчета принята координата xгр, то координаты xверх= –0,5h; xниж= 0,5h. Преобразуем выражение (3.16):

 

 

 

 

 

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1exp AD

 

 

 

 

 

 

 

 

 

εк(ρ) =

 

 

 

 

 

 

 

2

 

 

 

.

(3.17)

 

 

h

 

 

 

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

AD

1

 

 

 

 

 

 

 

 

 

 

 

exp AD

 

 

 

exp

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

2

 

 

В последнем уравнении числитель и знаменатель правой части умножим на 2, к числителю прибавим и вычтем exp(AD-1h/2) и проведем ряд преобразований:

 

 

 

 

 

h

 

 

 

 

 

 

 

 

h

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

exp AD

 

 

 

 

 

exp AD

 

 

 

 

 

 

 

 

εк(ρ) = 0,5 0,5

 

 

 

 

2

 

 

 

 

 

 

 

 

2

 

.

 

(3.18)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h

 

 

 

 

 

 

 

h

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

exp AD

 

 

 

 

 

 

 

exp AD

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

2

 

 

 

Поскольку для гиперболических функций справедливо:

 

 

 

x

 

 

 

ch x 1

 

ex + ex 2

 

 

 

th

 

 

 

 

 

=

 

 

 

 

=

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sh x

 

 

ex + ex

 

 

 

 

 

 

2

 

 

 

 

 

 

 

выражение (3.18) окончательно запишем в виде:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 0,5 0,5th[gh(4αD)

(ρρср)]. .

εк(ρ) = 0,5 0,5th AD

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

18

Это закон гиперболического тангенса.

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

Сепарационная характеристика εк(ξ) любого аппарата имеет два важных элемента: положение границы разделения ξгр («рабочей» точки) и крутизну в «рабочей» точке:

ε(ξгр)=

к(ξ)

 

 

 

 

.

 

 

ξ=ξгр

 

 

Уравнение границы разделения εкгр)=0,5.

Крутизна сепарационной характеристики – мера точности разделения сепаратора. Увеличение крутизны означает повышение точности разделения в сепараторе, т. е. уменьшение взаимного засорения фракций. Для идеального се-

паратора ε´(ξгр)→∞.

Крутизна связана со средним вероятным отклонением Еpm, еще одной характеристикой точности разделения в аппарате, – соотношением:

ε´(ξгр)≈0,25/Еpm

Сепарационные характеристики для типовых аппаратов, полученные теоретическим путем с экспериментальным подтверждением, приведены в таблице 3.3, а их графическое изображение дано на рис. 3.5.

Основное использование сепарационных характеристик – расчет и проектирование схем обогащения.

εк(ξ)

а

1

Q1

0,5 Q2>Q1

0

εт(ρ)

1

tн

0,5

0

Q2>Q1

Q1

ξ

г

t≈2tн

t→∞

ρ

εк(ξ)

б

 

εн(d)

 

в

 

1

 

 

 

 

1

 

 

 

0,5

 

 

 

 

0,5

tmin

tmax

 

1 2

 

 

 

 

 

 

t→0

t→∞

 

2

1

 

 

 

 

 

d0

 

0

 

 

 

ξ

0

dгр(t)=var

d

εп(k) д

1 tmax

0,5

tmin

0

k

 

Рис. 3.5. Сепарационные характеристики различных обогатительных аппаратов:

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

19

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]