Modeling Superellipsoids using LIGGGHTS

Submitted by MahmoudY on Wed, 03/26/2025 - 18:15

Hi All,

I am working on simulating superellipsoids particles free falling into a cylinderical tank. I am aware that spehere in the code are easily modeled. However, I am not sure how to model shapes other than shperes, such as superellipsoids. I am studying the packing configuration of different particles shapes using DEM solver. I would appreciate your suggestions and comments on this problem.

Thanks in advance !

Mishal | Thu, 03/27/2025 - 14:17

Have a look at the link below. The package will essentially take a .stl file and fill the model with spheres to give you the required shape, the data is printed to a .txt file which can be read by LIGGGHTS. Have a look at the multisphere section in the readme files. I have added the Python code below which will convert the txt output from Clumps into the format for LIGGGHTS (sorry the indenting did not come out in the post).

https://pypi.org/project/clump-python/

with open(output, 'r') as F:
L = 0
for L, line in enumerate(F):
pass
L += 1

with open(output, 'r') as F:
data = F.readlines()

with open(output, 'w') as F:
for line in data:
line = line.replace(",", " ")
F.write(line)