Электронный научный журнал
Современные проблемы науки и образования
ISSN 2070-7428
"Перечень" ВАК
ИФ РИНЦ = 0,737

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ СМЕШИВАНИЯ СЫПУЧИХ МАТЕРИАЛОВ В ЛОПАСТНЫХ СМЕСИТЕЛЯХ

Смолин Д.О. 1 Дёмин О.В. 1 Першин В.Ф. 1
1 ФГБОУ ВПО «Тамбовский государственный технический университет»
Основываясь на исследованиях движения сыпучих материалов в аппаратах с лопастными рабочими органами, предложена математическая модель смешивания сыпучих материалов в лопастных смесителях, построенная на основе теории цепей Маркова. В работе движение частиц сыпучего материала рассматривается в трёх направлениях: осевом (между рядами рабочих органов), радиальном (в пределах рядов рабочих органов) и окружном (по замкнутому контуру). Основываясь на схеме перемещения частиц материала, предложена схема разбиения смесительной камеры с условием равенства объёмов ячеек. Показаны в общем виде матрицы переходных вероятностей и описан способ идентификации вероятности перехода частиц ключевого компонента в осевом, радиальном, окружном направлениях, а также представлены особенности состояния системы в случае двухвального лопастного смесителя.
цепи Маркова
лопастные смесители
смешивание сыпучих материалов
1. Баранцева Е.А. Процессы смешивания сыпучих материалов: моделирование, оптимизация, расчет / Е.А. Баранцева, В.Е. Мизонов, Ю.В. Хохлова. — Иваново : ГОУ ВПО «Ивановский государственный энергетический университет им. В.И. Ленина», 2008. — 116 с.
2. Конструирование и расчет машин химических производств : учебник для машиностроительных вузов по специальности «Химическое машиностроение и аппаратостроение» / Ю.И. Гусев, И.Н. Карасев, Э.Э. Кольман-Иванов и др. — М. : Машиностроение, 1985. — 408 с.
3. Дёмин О.В. Совершенствование методов расчёта и конструкций лопастных смесителей : дис. … канд. техн. наук. — Тамбов, 2003. — 210 с.
4. Марков А.А. Избранные труды. Теория чисел. Теория вероятности. — М. : АН СССР, 1951. — 717 с.
5. Макаров Ю.И. Аппараты для смешивания сыпучих материалов. — М. : Машиностроение, 1973. — 216 с.
6. Berthiaux H. Application of the theory of Markov chains to model different processes in particle technology / H. Berthiaux, V. Mizonov, V. Zhukov // Powder Technology. — 2005. — Vol. 157, № 1–3. — P. 128–137.

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

Опираясь на результаты предыдущих исследований, предложена математическая модель, в которой перемещение частиц в процессе смешивания представлено в трёх направлениях:

1) осевом (рис. 1а). Представляет собой перемещение частиц материала между рядами рабочих органов вдоль смесительной камеры;

2) радиальном (рис. 1б). Частицы материала перемещаются в пределах рядов рабочих органов смесителя;

3) окружном (рис. 1в). Перемещение частиц материала происходит по замкнутому циркуляционному контуру.

Рисунок 1 – Схема перемещения частиц и разделения на ячейки: (а) – осевое, (б) – радиальное, (в) – окружное.

Смесительная камера разделяется на ячейки равного объёма. В осевом направлении на одинаковых частей, равных числу рядов рабочих органов (рис. 1а). Нумерацию частей производим от одного торца камеры до другого.

Каждую часть, полученную при разбиении в осевом направлении, разобьём последовательно в радиальном направлении на колец (рис. 1б) и в окружном направлении на (рис. 1в) секторов:

(1)

где - угол сектора.

Количество ячеек в каждой полученной части определяется формулой:

(2)

Нумеруются ячейки (рис. 1б) с индексом приращения где - номер кольца.

Радиусы колец (рис. 1б), при условии одинакового объёма ячеек, вычисляются по формуле:

(3)

где - радиус смесительной камеры.

Состояние системы в момент времени (где - номер перехода, - длительность перехода), представим в виде вектора-столбца размером

(4)

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

(5)

где - матрица переходных вероятностей.

Матрица переходных вероятностей с учетом трёх направлений перемещения частиц материала определяется по формуле:

(6)

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

(7)

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

(8)

где - нулевая матрица размером - блочная матрица размером .

(9)

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

Матрица переходных вероятностей при перемещении частиц в радиальном направлении имеет следующий вид:

(10)

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

(11)

где - блочная матрица размером :

(12)

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

Матрица переходных вероятностей при перемещении частиц в осевом направлении имеет вид:

(13)

где - нулевая матрица размером , а - блочная матрица, .

(14)

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

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

(15)

где и - состояния системы при вращении I и II-го валов соответственно (рис. 2), определяемые по формуле (4).

 

Рисунок 2 – Схема перемещения частиц и разделения на ячейки в двухвальном лопастном смесителе.

Состояние системы можно представить в матричной форме:

(16)

где - матрица смешивания, имеющая вид:

, (17)

где - нулевая матрица размером .

- матрица обмена частиц материала при вращении I и II-го валов. Элементы матрицы, кроме элементов , где , , , , равны нулю.

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

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

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

Рецензенты:

Мизонов В.Е., д.т.н., профессор, заведующий кафедрой прикладной математики ФГБОУ ВПО «Ивановский государственный энергетический университет им. В.И. Ленина», г. Иваново.

Лабутин А.Н., д.т.н., профессор, заведующий кафедрой технической кибернетики и автоматики, ФГБОУ ВПО «Ивановский государственный химико-технологический университет», г. Иваново.


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

Смолин Д.О., Дёмин О.В., Першин В.Ф. МАТЕМАТИЧЕСКАЯ МОДЕЛЬ СМЕШИВАНИЯ СЫПУЧИХ МАТЕРИАЛОВ В ЛОПАСТНЫХ СМЕСИТЕЛЯХ // Современные проблемы науки и образования. – 2013. – № 2.;
URL: http://science-education.ru/ru/article/view?id=8703 (дата обращения: 25.05.2019).

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

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