Метод дискретных элементов

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

Гранулированный материал представляется совокупностью из N обособленных упругих частиц сферической формы радиуса Ri (рисунок 1.а). Согласно законам классической Ньютоновской механики, движение каждой i-й частицы полностью определяется координатами ее центра тяжести xi и углом поворота θi вокруг центра тяжести как жесткого целого (i=1,…,N).

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

(1)

(2)

где t – время, mi – масса частицы, Ii – момент инерции, векторы Fпредставляют собой сумму сил, действующих на контактах i-й и j-й частиц (включая силу тяжести):

(3)

Вектор Mi – момент сил Fij относительно центра i-й частицы:

(4)

Поверхностные силы Fij состоят из сил отталкивания Fn,ij и сил трения Ft,ij. Сила отталкивания возникает между частицами при условии δij > 0 (рисунок 1.а) и направлена по нормали nij в направлении центра i-й частицы (рисунок 1.б). Для ее определения выберем вязкоупругую модель соударения:
где Fn,ij – упругая; Fn,ij – вязкая составляющие.

Согласно Дж. Герцу, упругая часть силы Fn,ij равна:

(5)

где vi – коэффициент Пуассона, Ei – модуль упругости частицы. В случае соударения частицы и прямолинейной границы, в формуле (5) достаточно положить один из радиусов равным бесконечности.

Рисунок 1  Задание параметров для описания контакта, δij – величина перекрытия частиц:

геометрия модели (а); силы, действующие на контакте (б)

Вязкая составляющая силы отталкивания определяется из соотношения:

(6)

где Mij – приведенная масса частиц, un,ij – проекция относительной скорости точки соударения на ось nijγn – коэффициент демпфирования, оказывающий основное влияние на коэффициент восстановления скорости после удара.

Сила трения Ft,ij направлена против движения i-й частицы относительно j-й, а ее величина определяется соотношением:

(7)

где ut,ij – проекция скорости точки контакта Ci относительно скорости точки Cj на ось tijφt – угол контактного трения между частицами.

Таким образом, система дифференциальных уравнений второго порядка (1–2) относительно неизвестных xiθi и соотношения (3–7) полностью определяют движение и соударение совокупности упругих частиц, моделирующих дискретную среду [1].

Метод дискретных элементов реализуется как на плоскости (рисунок 2), так и в трехмерной системе координат (рисунок 3).

Рисунок 2  Реализация метода дискретных элементов на плоскости

Рисунок 3  Модель плотной упаковки сфер по методу дискретных элементов: 10000 частиц (а); 30000 частиц (б)

Наиболее известны следующие программы, реализующие метод дискретных элементов: Chute Maven (Hustrulid Technologies Inc.), PFC2D и PFC3D, EDEM (DEM Solutions Ltd.), GROMOS-96, ELFEN, MIMES, PASSAGE®/DEM [2].


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

[1]  Клишин, C.B. Применение метода дискретных элементов при анализе гравитационного движения гранулированного материала в сходящемся канале / C.B. Клишин // Горный информационно-аналитический бюллетень.(научно-технический журнал), 2009, №12, С.273-277.

[2]  Арсентьев, В.А. Методы динамики частиц и дискретных элементов как инструмент исследования и оптимизации процессов переработки природных и техногенных материалов / В.А. Арсентьев, И.И. Блехман, Л.И. Блехман, Л.А. Вайсберг, К.С. Иванов, А.М. Кривцов // Обогащение руд. 2010. № 1. С. 30–35.