ТЕСТИРОВАНИЕ КОМПЬЮТЕРНЫХ ПРОГРАММ, ИСПОЛЬЗУЕМЫХ ПРИ МОДЕЛИРОВАНИИ МЕТОДОМ ДИСКРЕТНЫХ ЭЛЕМЕНТОВ

Аннотация


Для исследования динамики сыпучих материалов используется моделирование методом дискретных элементов (Discrete element method - DEM ). Для верификации программной реализации метода, а также для выбора временного шага интегрирования требуется наличие надежных тестовых задач. В настоящей статье представлено тестовое аналитическое решение, полученное в рамках теории удара. Рассматриваются варианты соударения двух сферических частиц со скольжением и вращением. Для аналитического решения была принята гипотеза, что упругая деформация при ударе близка к нулю, пластическая деформация отсутствует, потеря нормальной скорости при ударе описывается коэффициентом восстановления, сила тяжести не учитывается. Начальная скорость движения частиц известна. Получены решения для случая идентичных частиц, а также для частиц разного диаметра. Аналитическое решение используется для верификации собственной программной реализации метода DEM. Проведена серия вычислительных экспериментов с выбором разных параметров соударяющихся частиц и параметров метода. Приводятся графики сравнения тангенциальной, нормальной скорости частиц и скорости вращения частицы после соударения для аналитического и численного решения. Показано, что после удара эти параметры имеют качественное и количественное сходство. Для исследования сходимости метода и для определения оптимального вычислительного шага по времени проведена дополнительная серия вычислительных экспериментов. В качестве параметра, характеризующего точность вычислений, выбран модуль скорости движения части после удара. Показана сходимость реализации метода при уменьшении вычислительного шага по времени. На основании этих данных можно выбрать оптимальный шаг по времени для вычислений.

Полный текст

Для исследования динамики сыпучих материалов [1-4], прессования, гранулирования, измельчения широко используется моделирование методом дискретных элементов (Discrete element method - DEM, или молекулярная динамика) [5-13]. В настоящее время благодаря увеличению компьютерной производительности стало возможно моделировать сложные системы с миллионами частиц в трех измерениях. Программная реализация метода требует наличия надежных тестов для проверки правильности реализации метода. Даже использование готовых программных комплексов требует проверки правильности вычислений при использовании конкретных гипотез и замыкающих соотношений. Целью работы является разработка тестов, позволяющих находить ошибки и верифицировать компьютерные программы, использующие DEM, находить оптимальный шаг расчета по времени, изменять параметры для повышения производительности расчетов при допустимой погрешности. В настоящей работе с применением теории удара [14-17] выведены формулы (1)-(3), необходимые для тестирования соударения двух сферических частиц со скольжением и вращением. С использованием этих формул проведено верифицирование компьютерной программы трехмерного моделирования методом DEM с шарообразными частицами, разработанной авторами [18-26]. Производится сравнение характеристик, полученных в результате использования компьютерной программы DEM, и характеристик, полученных с применением аналитических формул теории удара. Далее используются следующие обозначения: - скорость вращения частицы до удара, рад/с; - скорость вращения частицы после удара, рад/с; - коэффициент трения; k - коэффициент восстановления; - скорость частицы до удара, м/с; - составляющая скорости частицы до удара, нормальная к поверхности второй частицы (поверхности стенки), м/с; - составляющая скорости частицы до удара, касательная к поверхности второй частицы (поверхности стенки), м/с; - составляющая скорости частицы после удара, нормальная к поверхности второй частицы (поверхности стенки), м/с; - составляющая скорости частицы после удара, касательная к поверхности второй частицы (поверхности стенки), м/с; r - радиус частицы, м; - угол между направлением удара и поверхностью второй частицы (или угол между направлением удара и стенкой), °. Рассмотрим моделирование косого удара двух частиц (рис. 1, а). Если две частицы имеют идентичные параметры, то можно упростить задачу, выбрав относительную систему координат так, что две частицы двигаются навстречу друг другу, причем их скорости параллельны оси Y, эта постановка используется при моделировании (рис. 1, б). Если две частицы идентичны, то задача упрощается до удара одной частицы об абсолютно твердую стенку (рис. 1, в), которая является касательной к обеим частицам в момент удара в точке соприкосновения. Задача в такой постановке может быть решена аналитически. Для использования этого решения в качестве теста реализации численного метода стенка заменяется абсолютно неподвижной частицей (рис. 1, г). Таким образом, при моделировании задача может рассматриваться либо как удар двух частиц (рис. 1, б), двигающихся с одинаковой скоростью навстречу, либо как удар частицы о неподвижную частицу (рис. 1, г). Для аналитического решения была приняты гипотеза, что упругая деформация при ударе близка к нулю, пластическая деформация отсутствует, потеря нормальной скорости при ударе описывается коэффициентом восстановления, сила тяжести не учитывается. Начальная скорость движения частицы известна. Трение без скольжения: при (1) Трение со скольжением: (2) где (3) Используем это аналитическое решение для тестирования компьютерной программы. В качестве параметров, по которым будут сравниваться численное и аналитическое решения, используются тангенциальные и нормальные скорости частиц и скорости вращения частицы после соударения. а б в г Рис. 1. Схема расположения частиц для расчета соударения двух частиц одинакового размера Рис. 2. Тангенциальная скорость, нормальная скорость после удара При расчетах использовались параметры: - V0 = 0,1 м/с; - шаг по времени dt = 1/300000 c; - плотность материала ρ = 562 кг/м; - модуль Юнга E = 1000 МПа; - радиусы частиц равны, r = 0,0008 м; - скорость вращения частицы до удара ω0 = 0 м/с. Результаты вычисления при α = 0…90°, k = 1, µ = 0 приведены на рис. 2. Так как µ = 0, вращение после удара отсутствует, ω = 0 (рис. 3). Результаты вычисления при α = 0…90°, k = 1, µ = 0,1 показаны на рис. 4-6. Результаты вычисления при α = 0…90°, µ = 0,1 приведены на рис. 7-9. Коэффициент восстановления в этом тесте зависит от нормальной скорости соударения частиц. Исследователями используются различные зависимости для определения силы реакции от скорости и других параметров. В конечном счете не имеет значения, как изменяется с течением времени сила реакции в направлении нормали при соударении. Для определения скоростей после удара аналитическим способом достаточно знать ударный импульс, который определяется коэффициентом восстановления. В предыдущих тестах использовался коэффициент восстановления равный единице. В этом тесте используется коэффициент восстановления, полученный в процессе моделирования. Таким образом, как в математическом моделировании, так и в аналитическом расчете используются одинаковые коэффициенты восстановления для каждого угла соударения частиц. В настоящей работе проверка коэффициента восстановления не показана, так как его зависимость от скорости соударения, свойств материалов частиц является частным случаем и не имеет общей формулы для применения. Рис. 3. Отношение скорости после удара к скорости до удара Рис. 4. Тангенциальная скорость, нормальная скорость после удара Рис. 5. Отношение скорости после удара к скорости до удара Рис. 6. Скорость вращения после удара, ω Рис. 7. Тангенциальная скорость, нормальная скорость после удара Рис. 8. Отношение скорости после удара к скорости до удара. Скорость Vn (аналитическое решение) не показана, так как равна Vn (dem) Рис. 9. Скорость вращения после удара, ω Схема соударения двух частиц разного размера следующая. Устанавливается стенка в плоскости x = 0. Частицы ударяются о стенку с двух сторон. Точка удара находится в ноле системы координат. (4) В таком случае частицы имеют по модулю равный нормальный и тангенциальный импульс: (5) Частицы двигаются на встречу, ударяются о стенку с одинаковым углом α. Соотношение скоростей после столкновения. (6) Рис. 10. Схема расположения частиц для расчета соударения двух частиц разного размера Рис. 11. Скорости частиц после удара Таким образом, удар двух частиц о стенку в рассматриваемой схеме равнозначен удару двух частиц без стенки. Аналитическое решение используется то же, формулы (1)-(3). Далее по тексту будут добавлены только индексы частиц. При компьютерном моделировании необходимо установить для частиц начальные скорости и координаты согласно схеме (рис. 10) и формулам (4)-(6). Вариант расчета при α = 0…90° град, k = 1, µ = 0, ω01 = 0, ω02 = 0, r1 = 0,0008 м, r2 = 0,0016 м. Погрешность по модулю скорости движения частиц после столкновения между составила между результатами программы DEM и аналитическим решением 2,13 %. Скорости вращения после столкновения равны нолю, так как µ = 0. Вариант расчета при α = 0…90° град, k = 1, µ = 0,1, ω01 = 0, ω02 = 0, r1 = 0,0008 м, r2 = 0,0016 м. По результатам теста погрешность по модулю скорости движения частиц после столкновения составила 2,11 %. Погрешность скорости вращения 1,48 %. Вариант расчета при α = 0…90° град, µ = 0,1, ω01 = 0, ω02 = 0, r1 = 0,0008 м, r2 = 0,0016 м. Для аналитического расчета используется коэффициент восстановления, полученный в программе DEM. По результатам теста погрешность по модулю скорости движения частиц после столкновения составила 0,76 %. Погрешность скорости вращения 1,15 %. Скорости движения, полученные с применением программы DEM, после удара частицы 2 меньше, чем скорости движения частицы 1, в 8 раз, то есть соблюдается формула (7). Углы отражения удара для обеих частиц равны, что подтверждает корректность работы программы DEM. Ниже представлены графики скоростей (рис. 11). Графики скоростей, полученные аналитическим способом, не показаны, так как погрешности менее 1,5 %. Таким образом, получено аналитическое решение для тестирования программ, реализующих DEM. Проведено тестирование собственной программы. Погрешность метода определялась по модулю скорости движения частицы после столкновения. В результате тестирования показано, что относительная погрешность вычислений программы DEM относительно аналитического решения при выбранном шаге по времени составляет не более 2,13 %. Полученное аналитическое решение можно использовать для проверки погрешности программ DEM, входящих в коммерческие пакеты.

Об авторах

Д. В Лобовиков

Пермский национальный исследовательский политехнический университет

А. В Харченко

ООО «Конструктив-СП»

Е. В Матыгуллина

Пермский национальный исследовательский политехнический университет

Список литературы

  1. Лобовиков Д.В., Ханов А.М., Храмов Б.Л. Математическая модель окатывания частиц в барабане, движущемся по планетарной траектории // Химия, технология и промышленная экология неорганических соединений: сб. науч. тр. / Перм. гос. техн. ун-т. - Вып. 5. - Пермь, 2002. - C. 143-153; 134-142.
  2. Dynamics of drag and force distributions for projectile impactin a granular medium / M.P. Ciamarra, A.H. Lara, A.T. Lee, D.I. Goldman, I. Vishik, H.L. Swinney // Phys. Rev. Lett. -2004. - Vol. 92, № 19. - P. 194301.
  3. Jop P., Forterre Y., Pouliquen O. A constitutive law for dense granular flows // Nature. - 2006. - Vol. 441, № 7094. - P. 727-730.
  4. Attractive particle interaction forces and packing density of fine glass powders / E.J.R. Parteli, J. Schmidt, C. Blümel, K.-E. Wirth, W. Peukert, T. Pöschel // Sci. Rep. - 2014. - Vol. 4. - P. 6227.
  5. Buchholtz1 V., Freund J.A., Poschel T. Molecular dynamics of comminution in ball mills // Eur. Phys. J. B. - 2000. - Vol. 16. - P. 169-182.
  6. Coefficient of restitution of colliding viscoelastic spheres / R. Ramırez, T. Poschel, N.V. Brilliantov, T. Schwager // PHYSICAL REVIEW E. - 1999. - Vol. 60 (4). - P. 4465-4472.
  7. Pöschel T., Buchholtz V.Complex flow of granular material in a rotating cylinder // Chaos, Solitons and Fractals. - 1995. - № 4. - P. 1901.
  8. Schwager T., Poschel T. Contact of viscoelastic spheres. “Friction, Arching, Contact Dynamics” // World Scientific. - Singapore, 1997. - P. 293-299.
  9. Volkhard Buchholtz and Thorsten Poschel. A Vectorized algorithm formolecular dynamics of short range interacting particles // International Journal of Modern Physics C. - 1993. - Vol. 4. - P. 1049.
  10. Poschel T., Schwager T.Computational Granular Dynamics // Springer Berlin Heidelberg. - New York, 2005.
  11. Cundall P.A., Strack O.D.L. A discrete numerical model for granular assemblies // Geotechnique, 1979. - Vol. 29. - P. 47.
  12. Hauger W., Schnell W., Gross D. Technische Mechanik. Bd 3: Kinetik 7. - Berlin: Springer, 2002. - 267 p.
  13. Sondergaard R., Chaney K., Brennen C.E. Measurements of solid spheres bouncing off flat plates // ASME J. Appl. Mech. - 1990. - Vol. 57, № 3. - P. 694-699.
  14. Stronge W.J. Impact Mechanics. - Cambridge: Cambridge University Press, 2004. - 280 p.
  15. Бать М.И., Джанелидзе Г.Ю., Кельзон А.С. Теоретическая механика в примерах и задачах. Т. 2: Динамика. - М.: Наука, 1966. - 663 с.
  16. Лобовиков Д.В., Ханов А.М., Храмов Б.Л. Условие адгезии упругопластических сферических тел // Вестник ПГТУ. Аэрокосмическая техника. - 2002. - № 13. - C. 67-71.
  17. Лобовиков Д.В. Влияние изменения параметров на процесс окатывания // Вестник ПГТУ. Механика и технология материалов и конструкций / Перм. гос. техн. ун-т. - Пермь, 2002. - C. 263-272.
  18. Давление в сыпучем материале при гранулировании в планетарном грануляторе / А.М. Ханов, Д.В. Лобовиков, Л.Д. Сиротенко, Е.В. Матыгуллина // Вестник ПГТУ. Проблемы современных материалов и технологий. - 2005. - № 11. - C. 163-169.
  19. Лобовиков Д.В. Образование гранул в планетарном грануляторе // Конструкции из композиционных материалов. - 2006. - Вып. 4. - C. 55-60.
  20. Лобовиков Д.В., Матыгуллина Е.В. Получение композиционных гранулированных материалов в планетарном грануляторе. - Пермь: Изд-во Перм. гос. техн. ун-та, 2008. - 153 с.
  21. Lobovikov D., Hanov A., Hramov B. Peculiarity of sliding of granular material on the surface of a rotating drum // Proceedings of XXX Summer School Advanced Problems in Mechanics 2002. - CПб.: Изд-во Института проблем машиноведения РАН, 2003. - С. 441-446.
  22. Lobovikov D. Dry granulation of powder in a drum. XXXI International Summer School. Conference "Advanced Problems in Mechanics": book of abstracts. - CПб., 2003. - C. 65-66.
  23. Lobovikov D. Experimental data of a granulation in a planetary granulator. XXXII International Summer School. Conference "Advanced Problems in Mechanics": book of abstracts. - CПб., 2004. - C. 68-69.

Статистика

Просмотры

Аннотация - 87

PDF (Russian) - 35

Ссылки

  • Ссылки не определены.

Данный сайт использует cookie-файлы

Продолжая использовать наш сайт, вы даете согласие на обработку файлов cookie, которые обеспечивают правильную работу сайта.

О куки-файлах