PeriDEM
0.1.0
|
Implementation of high-fidelity model of granular media that combines the advantages of peridynamics and discrete element method (DEM). The model has following advantages over existing mechanical models for granular media:
For more details about the model and results, we refer to the paper:
Prashant K. Jha, Prathamesh S. Desai, Debdeep Bhattacharya, Robert P Lipton (2020). Peridynamics-based discrete element method (PeriDEM) model of granular systems involving breakage of arbitrarily shaped particles. Journal of the Mechanics and Physics of Solids, Volume 151, 2021. Doi https://doi.org/10.1016/j.jmps.2021.104376.
Download pdf here.
Some example cases. For more details, look at the problem_setup.py
file or input_0.yaml
file in examples.
To create example input files, the python script is provided. Python script allows parameterization of various modeling, geometrical parameters, and creating .geo
files for gmsh
and particle locations file. Typically, the input files consists of:
input.yaml
– the main instruction file for PeriDEM
with details about material models, particle geometries, time step, etcparticle_locations.csv
– this file provides location and other details of the individual particles. Each row in the file consists ofi
– zone id that particle belongs tox
– x-coordinate of the center of the particle. Next two columns are similarly for y
and z
coordinatesr
– radius of the particleo
– orientation in radians. This is used to give particle (particle mesh) a rotationmesh.msh
– mesh file for the reference particle or wall. For example, in compressive test example, there are total four mesh files: one each for circular and hexagon shaped particle and one each for fixed and mobile wallCircular without damping | Circular without damping |
Different materials | Different radius | Different radius different material |
Concave particles |
Setup for this test consists of 502 circular and hexagonal shaped particles of varrying radius and orientation within a rectangle container. The top wall of the container is moving downward with prescribed speed resulting in compression of the particle system. The quantity of interest is the compressive strength of media; it is expected that the moving wall register reaction for with increasing penetration, however, after a certain amount of compression of media, the damage will initiate in individual particles especially those connected by force chain resulting in yielding of the media. For more details, we refer to Jha et al 2021
Compressive test setup |
Top: Plot of reaction force per unit area on top wall. Bottom: Particle state at four times. Color shows the damage at nodes. Damage 1 or above indicates the presence of broken bonds in the neighborhood of a node. |
Compressive test simulation |
Core dependencies for building the executible:
Dependencies for running the examples:
If all the libraries are installed in global space, commands for building the code is as simple as
If libraries such as hpx, pcl are installed in custom paths, we will write
input.yaml
and all other files such as .msh
file for particle and wall and particle locations file are created and their filenames with paths are provided in the input.yaml
file, we will run the problem (using 4 threads)Some shell scripts are provided that may help in building the code on ubuntu 18.04, 20.04, and mac.
We also provide the docker image of the peridem code. The docker image is built using the three layers
Dockerfile associated to all layers listed above can be found in build scripts. Image corresponding to layer 2 above can be downloaded from the dockerhub:
docker pull prashjha/u1804-pd
docker pull prashjha/u2004-pd