Метод дискретного элемента


Метод дискретного элемента

Метод дискретного элемента (DEM, от англ. Discrete element method) — это семейство численных методов предназначенных для расчёта движения большого количества частиц, таких как молекулы, песчинки, гравий, галька и прочих гранулированных сред. Метод был первоначально применён Cundall в 1971 для решения задач механики горных пород. Williams, Hocking и Mustoe детализировали теоретические основа метода. В 1985 они показали, что DEM может быть рассмотрен как обобщение метода конечных элементов (МКЭ, FEM). В книге Numerical Modeling in Rock Mechanics, by Pande, G., Beer, G. and Williams, J.R. описано применение этого метода для решения геомеханических задач. Теоретические основы метода и возможности его применения неоднократно рассматривалось на 1-й, 2-й и 3-й Международной Конференции по Методам Дискретного Элемента. Williams, и Bicanic (см. ниже) опубликовали ряд журнальных статей описывающих современные тенденции в области DEM. В книге The Combined Finite-Discrete Element Method, Munjiza детально описано комбинирование Метода Конечного Элемента и Метода Дискретного Элемента.

Этот метод иногда называют молекулярной динамикой (MD), даже когда частицы не являются молекулами. Однако, в противоположность молекулярной динамике, этот метод может быть использован для моделирования частиц с несферичной поверхностью. Различными ответвлениями семейства DEM являются метод отдельных элементов (distinct element method), предложенный Cundall в 1971, обобщенный метод дискретного элемента (generalized discrete element method), предложенный Williams, Hocking и Mustoe в 1985, дискретный деформационный анализ (discontinuous deformation analysis) (DDA) предложенный Shi в 1988, и метод конечных дискретных элементов (finite-discrete element method), предложенный Munjiza и Owen в 2004.

Методы дискретного элемента очень требовательны к вычислительным ресурсам ЭВМ. Это ограничивает размер модели или количество используемых частиц. Прогресс в области вычислительной техники позволяет частично снять это ограничение за счет использования параллельной обработки данных. Альтернативой обработки всех частиц отдельно является обработка данных как сплошной среды. Например, если гранульный поток подобен газу или жидкости, можно использовать вычислительную гидродинамику.

Содержание

Применение

Фундаментальным предположением метода является то, что материал состоит из отдельных, дискретных частиц. Эти частицы могут иметь различные поверхности и свойства. Примеры:

Типичные отрасли промышленности использующие DEM:

  • Горнодобывающая
  • Фармацевтическая
  • Нефтегазовая
  • Сельскохозяйственная
  • Химическая

Основные принципы метода

Моделирование МДЭ начинается c помещения всех частиц в конкретное положение и придания им начальной скорости. Затем силы, воздействующие на каждую частицу, рассчитываются, исходя из начальных данных и соответствующих физических законов.

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

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

На молекулярном уровне, мы можем рассматривать

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

Типичные методы интеграции используемые в методе дискретного элемента:

Дальнодействующие силы

Когда во внимание принимаются дальнодействующие силы (гравитация, сила Кулона), взаимодействия каждой пары частиц необходимо рассчитывать. Число взаимодействий, а следовательно, ресурсоёмкость расчёта, возрастает с увеличением количества частиц квадратично, что не приемлемо для моделей с большим числом частиц. Возможный путь решить эту проблему — объединить некоторые частицы, которые находятся на расстоянии от рассматриваемой частицы, в одну псевдочастицу. Рассмотрим, например, взаимодействие между звездой и отдаленнной галактикой: ошибка, возникающая из-за объединения массы всех звезд в удалённой галактике в одну точку, незначительна. Для того, чтобы определить, какие частицы могут быть объединены в одну псевдочастицу, используются так называемые древесные алгоритмы. Эти алгоритмы распределяют все частицы в виде дерева, квадрадерева в случае двухмерной модели и октадерева в случае трехмерной модели.

Модели в молекулярной динамике делят пространство, в котором происходит моделируемый процесс, на ячейки. Частицы, уходящие через одну сторону ячейки просто вставляются с другой стороны (периодические граничные условия); так же происходит и с силами. Силы перестают приниматься в расчёт после так называемой дистанции отсечения (обычно половина длины ячейки), так что на частицу не воздействует зеркальное расположение той же частицы на другой стороне ячейки. Таким образом, можно увеличивать количество частиц простым копированием ячеек.

Алгоритмы для обработки долговременных сил:

  • Barnes-Hut,
  • быстрый мультипольный метод.

Библиография

  • P.A. Cundall, O.D.L. Strack, A distinct element model for granular assemblies. Geotechnique, 29:47—65, 1979.
  • Williams, J.R., Hocking, G., and Mustoe, G.G.W., "The Theoretical Basis of the Discrete Element Method, " NUMETA 1985, Numerical Methods of Engineering, Theory and Applications, A.A. Balkema, Rotterdam, January 1985
  • Shi, G, Discontinuous deformation analysis — A new numerical model for the statics and dynamics of deformable block structures, 16pp. In 1st U.S. Conf. on Discrete Element Methods, Golden. CSM Press: Golden, CO, 1989.
  • Williams, J.R. and Pentland, A.P., "Superquadric and Modal Dynamics for Discrete Elements in Concurrent Design, " National Science Foundation Sponsored 1st U.S. Conference of Discrete Element Methods, Golden, CO, October 19-20, 1989.
  • Pande, G., Beer, G. and Williams, J.R., Numerical Modeling in Rock Mechanics, John Wiley and Sons, 1990.
  • Kawaguchi, T., Tanaka, T. and Tsuji, Y., Numerical simulation of two-dimensional fluidized beds using the discrete element method (comparison between the two- and three-dimensional models) Powder Technology, 96(2):129—138, 1998.
  • Griebel, Knapek, Zumbusch, Caglar: Numerische Simulation in der Molekulardynamik. Springer, 2004. ISBN 3-540-41856-3.
  • Bicanic, Ninad, Discrete Element Methods in Stein, de Borst, Hughes Encyclopedia of Computational Mechanics, Vol. 1. Wiley, 2004. ISBN 0-470-84699-2.
  • 2nd International Conference on Discrete Element Methods, Editors Williams, J.R. and Mustoe, G.G.W., IESL Press, 1992 ISBN 0-918062-88-8
  • Williams, J.R. and O’Connor, R., Discrete Element Simulation and the Contact Problem, Archives of Computational Methods in Engineering, Vol. 6, 4, 279—304, 1999
  • Ante Munjiza, The Combined Finite-Discrete Element Method Wiley, 2004, ISBN 0-470-84199-0

Программное обеспечение

Открытые источники и некоммерческое программное обеспечение:

  • BALL & TRUBAL (1979—1980) distinct element method (FORTRAN code), originally written by P.Cundall and currently maintained by C.Thornton.
  • SDEC Spherical Discrete Element Code.
  • YADE Yet Another Dynamic Engine, second incarnation of SDEC written from ground-up, GPL license.
  • LIGGGHTS Open Source Discrete Element Method Particle Simulation Code, распространяется и поддерживается Кжиштофом Клоссом

Доступные за дополнительную плату пакеты DEM-программ, включая PFC3D,EDEM и Passage/DEM:

См. также


Wikimedia Foundation. 2010.

Смотреть что такое "Метод дискретного элемента" в других словарях:

  • Метод подвижных клеточных автоматов — Подвижные клеточные автоматы активно меняют своих соседей за счет разрыва существующих связей между автоматами и образования новых связей (моделирование контактного взаимодействи …   Википедия

  • Метод граничных элементов — Метод конечных элементов: триангуляция Метод конечных элементов (МКЭ) численный метод решения задач прикладной механики. Широко используется для решения задач механики деформируемого твёрдого тела, теплообмена, гидродинамики и электромагнитных… …   Википедия

  • Метод конечных объёмов — Метод конечных объёмов  численный метод интегрирования систем дифференциальных уравнений в частных производных. Содержание 1 Описание 1.1 Неформальное 1.2 Математическое …   Википедия

  • Метод конечных элементов — Решение методом конечных элементов двухмерной магнитостатической задачи (линии и цвет означают направление и величину магнитной индукции) …   Википедия

  • МКЭ — Метод конечных элементов: триангуляция Метод конечных элементов (МКЭ) численный метод решения задач прикладной механики. Широко используется для решения задач механики деформируемого твёрдого тела, теплообмена, гидродинамики и электромагнитных… …   Википедия

  • Компьютерное моделирование — краш теста методом конечных элементов. Компьютерная модель (англ. computer model), или численная мод …   Википедия

  • Вычислительная механика — Вычислительная механика  раздел механики сплошных сред, в котором строятся конечномерные модели сплошных сред, используется компьютерное моделирование и численные методы для решения задач механики деформируемого твёрдого тела и механики… …   Википедия

  • МОДЕЛИРОВАНИЕ В СОЦИОЛОГИИ — метод исследования соц. явлений и процессов на их моделях, т. е. опосредствованное изучение соц. ооъектов, в процессе к рого они воспроизводятся в вспомогательной системе (модели), замещающей в познавательном процессе оригинал и позволяющей… …   Российская социологическая энциклопедия

  • Список алгоритмов — Эта страница информационный список. Основная статья: Алгоритм Ниже приводится список алгоритмов, группированный по категориям. Более детальные сведения приводятся в списке структур данных и …   Википедия

  • Программируемые алгоритмы —       Служебный список статей, созданный для координации работ по развитию темы.   Данное предупреждение не устанавл …   Википедия


Поделиться ссылкой на выделенное

Прямая ссылка:
Нажмите правой клавишей мыши и выберите «Копировать ссылку»

We are using cookies for the best presentation of our site. Continuing to use this site, you agree with this.