Измельчение старого асфальтобетона обычно включает две стадии. На первой стадии, например при измельчении в молотковой дробилке, происходит дробление крупных кусков. На второй стадии происходит тонкое измельчение, которое позволяет выделить битумно-минеральную составляющую.
При тонком измельчении используются как аппараты непрерывного действия, так и аппараты периодического действия. Следует отметить, что для описания аппаратов непрерывного действия могут использоваться уточненные модели аппаратов периодического действия [1].
В данной работе рассматривается процесс тонкого измельчения старого асфальтобетона в электромагнитном поле [2]. Описание процесса измельчения осуществляется с использованием функций измельчения [1]. При этом сами функции измельчения определяются на основе эксперимента.
Процесс порционного измельчения в общем случае может быть описан следующим кинетическим уравнением [3]:
(1)
где F(y, t) - интегральная функция распределения частиц по размерам; S(x) - селективная функция измельчения; B(y, x) - распределительная функция измельчения.
Распределительная функция измельчения должна удовлетворять следующим условиям:
B(y, 0) = 0; B(y, y) = 1. (2)
Для решения интегро-дифференциального уравнения (1) должны быть поставлены начальные и граничные условия, которые определяются из следующих соображений. В начальный момент времени распределение частиц по размерам задается некоторой функцией от аргумента
F(y, 0) = φ(y). (3)
Интегральная функция распределения частиц по размерам удовлетворяет двум граничным условиям, которые вытекают из ее определения
F(0, t) = 0; F(ymax, t) = 1, (4)
где ymax - максимальный размер измельчаемых частиц.
Для решения интегро-дифференциального уравнения (1) относительно F(y, t) будем использовать численный анализ для значений y [ymin, ymax] и t
[0, T] методом сеток [4], сделав следующие разбиения y и t с шагами, которые определяются следующими формулами:
,
(5)
где n, m - соответственно число точек для переменных y и t. Время измельчения T определяется степенью измельчения, которая вычисляется по известной формуле:
(6)
где - средний размер частиц до измельчения;
- средний размер частиц после измельчения.
В соответствии с численным анализом используем приближенные выражения для производных:
(i = 0, 1, 2, ..., n j = 0, 1, 2, ..., m); (7)
(i = 0, 1, 2, ..., n j = 0, 1, 2, ..., m - 1); (8)
(k = 0, 1, 2, ..., n - 1 j = 0, 1, 2, ..., m). (9)
В формуле (9) учитывается тот факт, что в уравнении (1) стоит под знаком интеграла, а интегрирование ведется по переменной x, значения которой определяются индексом k.
Подставляя формулы (7)-(9) в уравнение (1), получим:
(10)
(i = 0, 1, 2, ..., n j = 0, 1, 2, ..., m - 1).
Для удобства введем следующие обозначения:
(11)
В результате получаем следующее уравнение:
(i = 0, 1, 2, ..., n j = 0, 1, 2, ..., m - 1). (12)
Используя начальные условия
Fi,0 = φ(xi). (13)
и граничные условия
F0,j = 0; Fi,m = 1, (14)
стандартным методом сеток [4] находятся значения массива Fi,1, Fi,2, ..., Fi,m (i = 0, 1, 2, ..., n).
Следующим этапом решения задачи является определение функций измельчения Bi,k, Sk. В известной работе [5], актуальность которой не потеряна и в наше время, были сделаны первые попытки экспериментального определения данных функций. Однако из-за сложной техники эксперимента предложенный авторами метод не нашел широкого применения. Тем не менее общие принципы определения функций измельчения, заложенные в работе [5], могут быть использованы при другой технике эксперимента.
В настоящей работе при определении функций измельчения будем исходить из того, что в процессе измельчения каждая фракция, имеющая свои размеры частиц, будет вести себя независимо от других. В соответствии с этим постулатом в ходе экспериментальных исследований изучалось измельчение отдельных фракций, которые были получены путем ситового анализа.
Пусть F* весовая доля материала, крупность которого меньше размера отверстия сита с размером x, который является нижней границей крупности данной фракции. Естественно, что F* должно быть функцией x и времени измельчения. Тогда весовая доля материала, соответствующая приращению размера частиц dx, будет равна:
(15)
Соответственно доля материала, подвергшегося измельчению за единицу времени, будет равна
(16)
Согласно определению, селективная функция S(x) является весовой долей частиц размером x, измельченных за единицу времени. В результате получаем следующее выражение:
(17)
В работе [5] показано, что выражение (17) может быть приведено к более простому и удобному для вычислений виду
(18)
где - средний размер частиц данной фракции;
- весовая доля материала, остающегося на сите с размером, соответствующим нижней границе крупности данной фракции; t - время измельчения.
В соответствии с выражением (15) весовая доля материала, измельченного за время dt, будет равна
(19)
Пусть ∂F*(y, t) является весовой долей измельченного материала за время t, которая приходится на фракцию с нижней границей крупности, соответствующей ситу с размером y. Тогда в соответствии с определением распределительной функции получим
(20)
В работе [5] показано, что выражение (20) может быть приведено к более простому и удобному для вычислений виду
(21)
где - средний размер частиц данной фракции; F*(y, t) - весовая доля измельченного материала, которая приходится на фракцию с нижней границей крупности, соответствующей ситу с размером y.
При выводе формул (18), (21) предполагается, что ни одна частица не подвергается повторному измельчению, т.е., проще говоря, измельчаются только частицы исходной крупности. Как показали эксперименты, это предположение справедливо только для времени измельчения, не превышающего 25 с.
Для определения функций измельчения были проведены эксперименты на лабораторной установке измельчения старого асфальтобетона в электромагнитном поле [6]. В экспериментах изменялись следующие параметры: соотношение объема Vмел, занимаемого мелющими телами, и объема Vизм, занимаемого измельчаемым материалом, определяемое величиной
(22)
напряженность магнитного поля H.
На рис. 1, 2 представлены графики функций измельчения, построенные с использованием формул (18), (21) по экспериментальным данным. Во всех экспериментах соотношение объема Vмел, занимаемого мелющими телами, и объема Vизм, занимаемого измельчаемым материалом, принимало значение δ = 2.
Рис. 1. График селективной функции измельчения старого асфальтобетона в электромагнитном поле
Рис. 2. График распределительной функции B(y,x) измельчения старого асфальтобетона в электромагнитном поле
Анализ полученных экспериментальных данных показал, что в общем случае функции измельчения могут быть описаны выражениями
(23)
(24)
где k, α - некоторые константы.
На рис. 3 представлено сопоставление экспериментальных и расчетных данных, определенных по предложенной модели.
Статистическая проверка соответствия экспериментальных данных по измельчению старого асфальтобетона в ЭМИ распределению, определенному по формуле (10), осуществлялась по критерию Пирсона [7]. Проверка позволяет сделать вывод о соответствии статистических данных распределению, определяемому формулой (10) с доверительной вероятностью, превышающей 0,9.
Рис. 3. Сопоставление расчетных и экспериментальных данных по гранулометрическому составу измельченного материала: H = 45 ка/м; δ = 2; α = 0,013; k = 1,53
Рецензенты:
-
Приоров А.Л., д.т.н., доцент, Ярославский государственный университет им. П.Г. Демидова, г. Ярославль.
-
Басурин В.И., д.ф.-м.н., профессор, Ярославский филиал Московского института инженеров транспорта (университет), г. Ярославль.
Работа поступила в редакцию 20.02.2012.
Библиографическая ссылка
Бадоев В.А., Лупанов А.П., Мурашов А.А. МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ИЗМЕЛЬЧЕНИЯ СТАРОГО АСФАЛЬТОБЕТОНА В ЭЛЕКТРОМАГНИТНОМ ПОЛЕ // Фундаментальные исследования. – 2012. – № 3-3. – С. 624-627;URL: https://fundamental-research.ru/ru/article/view?id=29770 (дата обращения: 17.08.2024).