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

   Для моделирования механического поведения твердых частиц сыпучих материалов широко применяется метод частиц или метод дискретных элементов (МДЭ, 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.

При копировании материалов ссылка на сайт www.sunspire.ru обязательна. Также, вы можете использовать библиографическую ссылку на учебное пособие:

 

"Белов, В.В. Компьютерная реализация решения научно-технических и образовательных задач: учебное пособие / В.В. Белов, И.В. Образцов, В.К. Иванов, Е.Н. Коноплев // Тверь: ТвГТУ, 2015. 108 с."

Официальная группа ВК:

Сайты-партнеры:

Центр научно-образовательных электронных ресурсов Тверского государственного технического университета:

http://cdokp.tstu.tver.ru

 

Официальный сайт Тверского государственного технического университета:

http://www.tstu.tver.ru

Русскоязычный форум по языку программирования Dark Basic Professional:

http://area.mediahouse.ru/

Открытия, Изобретения, Новые технические разработки:

http://www.belashov.info

Электронный магазин 3D моделей CGTrader:

http://www.cgtrader.com

Рейтинг@Mail.ru
Яндекс.Метрика
Топ Разработка игр
Bourabai Research - Технологии XXI века