Научный журнал
Современные наукоемкие технологии
ISSN 1812-7320
"Перечень" ВАК
ИФ РИНЦ = 0,940

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ПРОЦЕССОВ КЛАССИФИКАЦИИ ЗЕРНИСТЫХ МАТЕРИАЛОВ НА СИТАХ

Ахмадиев Ф.Г. 1 Гиззятов Р.Ф. 1 Назипов И.Т. 1
1 ФГАОУ ВО «Казанский государственный архитектурно-строительный университет»
Рассмотрен процесс выделения целевых продуктов из зернистого материала на ситах с учетом стохастической природы подобных процессов, в качестве случайного процесса выбрана концентрация числа частиц выделяемых продуктов на поверхности сит (1/м). Для изучения распределения случайного процесса вдоль сит классификатора построена система стохастических дифференциальных уравнений. Коэффициенты дифференциальных уравнений найдены в зависимости от скорости вибрационного движения сыпучей среды по поверхности сит классификатора и вероятности просеивания частиц в отверстия ячеек, которые определяют количество прохода рассматриваемых частиц через ячейку за единицу времени, тем самым характеризуют скорость изменения случайного процесса. Получено решение системы стохастических дифференциальных уравнений и определены числовые характеристики случайного процесса. В общем случае уравнение Колмогорова – Фоккера – Планка относительно плотности вероятности случайного процесса для первого сита с помощью специально подобранной замены приведено к простейшему уравнению диффузии и решение уравнения записано в виде интегральной свертки. При аппроксимации случайного процесса белым шумом по найденным его числовым характеристикам записано приближенное решение для плотности распределения и вероятностей перехода, которые являются фундаментальными решениями задачи Коши при дельтаобразных начальных условиях. Приведены результаты расчетов распределения среднего значения случайного процесса для выбранных фракций вдоль первого и второго сита многоярусного классификатора.
математическая модель
разделение
зернистый материал
вероятность просеивания
случайный процесс
1. Вайсберг Л.А., Иванов К.С., Мельников А.Е. Совершенствование подходов к математическому моделированию процессов вибрационного грохочения // Обогащение руд. 2013. № 2. С. 21–27.
2. Непомнящий Е.А. Кинетика некоторых процессов переработки дисперсных материалов // Теоретические основы химической технологии. 1973. Т. 7. № 5. С. 754–763.
3. Ахмадиев Ф.Г., Гиззятов Р.Ф., Назипов И.Т. Гидрогазодинамика и кинетика разделения дисперсных сред на ситовых классификаторах // Инженерно-физический журнал. 2017. Т. 90. № 5. С. 1135–1145.
4. Ахмадиев Ф.Г., Гиззятов Р.Ф. Моделирование разделения зернистых материалов на многокаскадных классификаторах на основе теории случайных процессов // Теоретические основы химической технологии. 2018. Т. 52. № 3. С. 306–317.
5. Огурцов В.А., Горохова Е.Р., Огурцов А.В. Моделирование движения частиц при виброгрохочении на основе теории цепей Маркова // Строительство и реконструкция. 2011. № 5(37). С. 85–88.
6. Долгунин В.Н., Иванов О.О., Уколов А.А. Кинетика сегрегации частиц различной шероховатости и упругости при быстром гравитационном течении зернистой среды // Теоретические основы химической технологии. 2009. Т. 43. № 2. С. 199–207.
7. Ахмадиев Ф.Г., Гиззятов Р.Ф., Киямов Х.Г. Математическое моделирование процесса тонкослойного разделения зернистых материалов по размерам на ситовых классификаторах // Теоретические основы химической технологии. 2013. Т. 47. № 3. С. 309–317.
8. Казаков В.А. Введение в теорию марковских процессов и некоторые радиотехнические задачи. М.: Сов. радио, 1973. 232 с.

Для выделения целевых продуктов из зернистого материала используют различные конструкции аппаратов, однако наиболее эффективным из них является ситовое оборудование. Математическое моделирование происходящих при этом процессов является основой для технологического расчета ситовых классификаторов и оптимального управления процессами. В работе [1] и в ряде других работ этих авторов изучено движение частиц в колеблющейся среде, рассмотрены различные модели вибрационного движения и получены зависимости для средней скорости, скорости сегрегации. Следует также отметить, что работы этих авторов были одними из первых в области вибрационного разделения сыпучих сред. В работе [2] процесс разделения дисперсных сред на ситовых классификаторах рассмотрен как диффузионный процесс и изменение концентрации числа проходовых частиц по толщине слоя зернистого материала в зависимости от времени в направлении к вибрирующей поверхности изучено с позиций марковских процессов и описано уравнением Колмогорова – Фоккера – Планка (К. – Ф. – П.). В работах [3, 4] процесс выделения целевых продуктов на многокаскадных классификаторах также изучен с использованием теории случайных процессов, в частности в работе [4] для описания процесса перемещения частиц по поверхностям сит вследствие просеивания применяется теория Пуассоновских процессов. Изучению сегрегации частиц в виброожиженном слое на основе марковских цепей посвящена работа [5]. В работе [6] приведены основные принципы организации технологических процессов с управляемыми сегрегированными потоками, которые возникают при переработке зернистых материалов. Таким образом, для изучения процесса выделения целевых продуктов из зернистого материала можно использовать разные подходы, однако с учетом случайного характера процесса в целом наиболее предпочтительным является стохастический подход.

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

Математическое моделирование

Изучим процесс выделения целевых продуктов из зернистого материала на ситовых классификаторах с учетом стохастической природы подобных процессов. Диаметр частиц разделяемого зернистого материала рассматривается как одномерная случайная величина и для описания дисперсности зернистого материала задается плотность распределения частиц по размерам G(d). Тогда долю частиц выделяемого продукта в исходном материале можно найти по формуле

ahmad01.wmf.

В качестве случайного процесса ahmad02.wmf рассмотрим концентрацию числа частиц целевого продукта (1/м) размерами dj на поверхности i-го сита в точке x от начала координат в момент времени t. Тогда кинетические уравнения для определенного выше случайного процесса в общем виде можно записать [3, 7]:

ahmad03.wmf

Здесь кинетические функции Fi характеризуют скорость просеивания на i-м ярусе, вектор ahmad04.wmf задает свойства материала, n – количество сит классификатора, ηi(t) – дельта-коррелированные во времени случайные функции, которые в каждый момент времени характеризуют отклонение концентрации числа частиц от ее среднего значения на i-м ярусе (белый шум), с характеристиками: ahmad05.wmf и ahmad06.wmf.

Здесь βi – интенсивность, ahmad07.wmf – спектральная плотность белого шума. Особенностью этих уравнений, которая позволяет назвать их стохастическими, является наличие воздействия в виде белого шума.

Число частиц выбранной фракции на поверхности i-го сита в любой момент времени ahmad08.wmf определяется по выражению

ahmad09.wmf,

где Li – длина сита, Vi cp – средняя скорость вибрационного движения. Отклонение числа частиц ahmad10.wmf от среднего значения в любой момент времени связано с вероятностью просеивания частиц в ячейки сита, которая рассматривается как случайное событие. Эта вероятность зависит от размеров и формы самих частиц, гранулометрического состава разделяемого материала, условий стесненности и других факторов, а также от относительной скорости вибрационного движения материала. В области отверстия ячейки относительная скорость частицы может принимать любые значения от нуля до максимального. Поэтому зависимость вероятности просеивания от относительной скорости рассматривается как случайное событие с нормальным распределением [3]. В каждый фиксированный момент времени ahmad11.wmf из-за различной вероятности просеивания в ячейки сит частиц выбранной фракции их число на поверхности сита ahmad12.wmf является случайной величиной, т.е. число частиц ahmad13.wmf рассматривается как случайный процесс. Аппроксимация случайного процесса Ni белым шумом возможна, так как время корреляции случайного процесса много меньше среднего времени пребывания выбранных частиц на поверхности сита. С учетом свойств белого шума [8] процесс Ni является марковским процессом. Для его изучения можно использовать математический аппарат теории марковских процессов и плотность распределения ahmad14.wmf найти из решения системы уравнений Колмогорова – Фоккера – Планка (К. – Ф. – П.) [8]:

ahmad15.wmf.

Здесь приняты обозначения: ahmad16.wmf коэффициенты диффузии. Решение уравнений (К. – Ф. – П.) при больших значениях n представляет собой сложную задачу. Поэтому для изучения случайного процесса Ni построим систему стохастических дифференциальных уравнений, которые описывают его изменение вдоль сит классификатора в зависимости от времени.

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

ahmad17.wmf (1)

Коэффициенты кинетических уравнений вычисляются по зависимости [2]:

ahmad18.wmf (2)

Коэффициенты αi определяют количество прохода рассматриваемых частиц через ячейку за секунду, тем самым характеризуют скорость изменения случайного процесса. Здесь приняты обозначения: P – вероятность просеивания частиц в ячейку сита, 2ai – шаг i-го сита. Вероятность просеивания в ячейку рассматривается как сложное событие, состоящее из геометрической вероятности и вероятности, которая зависит от относительной скорости вибрационного движения. Вычисление вероятности просеивания в ячейку рассмотрено в работе [3]. Вибрационное движение сыпучих сред достаточно подробно изучено в работе [1]. Вычисление средней скорости вибрационного движения и амплитуды относительной скорости, которая используется при вычислении вероятности просеивания в ячейку, приведены в работе [3]. Рассмотрим решение дифференциальных уравнений (1) при следующих начальных и граничных условиях:

ahmad19.wmf при ahmad20.wmf и ahmad21.wmf при ahmad22.wmf (3)

Условия (3) определяют подачу разделяемого материала только на начало верхнего яруса многоярусного классификатора, ahmad23.wmf – число выбранных частиц (1/м), которые поступают на начало первого сверху сита. Используя замену ahmad24.wmf уравнения (1) с учетом условий (3) могут быть решены сведением к обыкновенным дифференциальным уравнениям [7]. Например, общее решение дифференциального уравнения для верхнего яруса можно записать в виде

ahmad25.wmf

где ahmad26.wmf определяет среднее значение случайного процесса ahmad27.wmf.

На основании решений системы (1) для каждого яруса можно определить среднее значение, корреляционную функцию и дисперсию процесса Ni. Среднее значение случайного процесса ahmad28.wmf для верхнего яруса с учетом условий (3) имеет вид

ahmad29.wmf.

Также с учетом свойств белого шума [8] можно найти дисперсию процесса:

ahmad30.wmfahmad31.wmf.

Решение дифференциального уравнения относительно среднего для второго сита из системы (1), с учетом решения для первого и условий (3), имеет вид

ahmad32.wmf

В итоге данное решение можно записать в виде

ahmad33.wmf

где ahmad34.wmf и ahmad35.wmf.

Решение дифференциального уравнения относительно среднего значения случайного процесса для третьего сита из системы (1), с учетом решений для первого и второго, а также начальных и граничных условий (3), имеет вид

ahmad36.wmf

Полученное решение для третьего сита можно представить в виде

ahmad37.wmf

ahmad38.wmf

ahmad39.wmf

где

ahmad40.wmf,

ahmad41.wmf,

ahmad42.wmf.

Используя метод индукции, можно показать, что решение дифференциального уравнения относительно среднего значения случайного процесса для i-го сита имеет вид

ahmad43.wmf (4)

где ahmad44.wmf

В общем случае для определения плотности распределения случайного процесса решение системы уравнений (К. – Ф. – П.) проводят численными методами. При аппроксимации белым шумом случайный процесс Ni является нормальным [8]. Поэтому, зная числовые характеристики случайного процесса, можно записать приближенное решение для плотности распределения и вероятностей перехода, которые являются фундаментальными решениями задачи Коши при дельтаобразных начальных условиях:

ahmad45.wmf

ahmad46.wmf

В общем случае N1 в начальном сечении является случайной величиной с плотностью вероятности ahmad47.wmf. Тогда уравнение (К. – Ф. – П.) для первого сита, с помощью замены

ahmad48.wmf ahmad49.wmf (5)

приводится к простейшему уравнению диффузии ahmad50.wmf и решение полученного уравнения записывается в виде интегральной свертки:

ahmad51.wmf (6)

Тогда общее решение уравнения (К. – Ф. – П.) для первого сита имеет вид

ahmad52.wmf (7)

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

ahmad1a.wmf

ahmad1b.wmf

Изменение среднего значения случайного процесса (концентрации частиц, 1/м) целевых продуктов вдоль первого и второго сита многоярусного классификатора: α11 = 5,78E-3; α21 = 8,42E-2; α31 = 2,57E-1; α41 = 9,08E-1; α12 = 3,66E-4; α22 = 2,77E-3; α32 = 7,18E-1; α42 = 6,93E-1 (сек-1); Vcp = 0,05 м/сек; размеры фракций: 1 – (0,8÷0,9)·10-3; 2 – (0,7÷0,8)·10-3; 3 – (0,6÷0,7)·10-3; 4 – (0,5÷0,6)·10-3 (м)

Заключение

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


Библиографическая ссылка

Ахмадиев Ф.Г., Гиззятов Р.Ф., Назипов И.Т. МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ПРОЦЕССОВ КЛАССИФИКАЦИИ ЗЕРНИСТЫХ МАТЕРИАЛОВ НА СИТАХ // Современные наукоемкие технологии. – 2019. – № 12-1. – С. 30-35;
URL: https://top-technologies.ru/ru/article/view?id=37828 (дата обращения: 18.07.2024).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1,674