New question #676451 on Yade:
https://answers.launchpad.net/yade/+question/676451
Hi everybody!
I'm trying simulate a mono-disperse vibro-fluidized granular media in a
cylinder with a cone-shape base.
I'm using the geom tool to build a facetcylinder+facetcone container and the
harmonicEngine to vibrate it.
Regarding the interactions, i'm using the HertzMindlin physics.
I've done exactly the same simulation with LAMMPS and run both on the same
machine using 4 cores:
once reached a stationary state, simulations with 200 or 300 particles
(gas-like behaviour) have, for the two softwares, almost the same computational
time (YADE is just 20% slower) while in the cases with more than 2000 grains
(liquid-like behaviour) LAMMPS is ten times faster!
This is the code i'm using to create the container and the engines. Is there a
specific tool to speed up YADE simulation in condition of vibro-fluidized
granuar system in the liquid phase namely when the particles are in a vibrating
box but constantly in contact to each other?
Thanks in advance.
Andrea
___________________________________________________
O.materials.append(FrictMat(young=yPMMA, poisson=poissPMMA,
frictionAngle=frictAngle,density=densPMMA, label='lmpPMMA'))
O.materials.append(FrictMat(young=ySteel, poisson=poisSteel,
frictionAngle=frictAngle,density=densSteel, label='lmpSteel'))
coneId=O.bodies.append(geom.facetCone(Vector3(0,0,hcone/2.),rcont,rconelow,hcone,orientation=Quaternion((0,
0, 1),0),wallMask=(False*1+True*2+True*4),material='lmpPMMA'))
cyliId=O.bodies.append(geom.facetCylinder(Vector3(0,0,hcil/2.+hcone),rcont,hcil,orientation=Quaternion((0,
0, 1),0),wallMask=(True*1+False*2+True*4),material='lmpPMMA'))
contenitore=coneId+cyliId
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb(),
],label='collider',verletDist=rball*2),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom(),
],
[Ip2_FrictMat_FrictMat_MindlinPhys(betan=0.4, betas=0.4),
],
[Law2_ScGeom_MindlinPhys_Mindlin(),
],
),
NewtonIntegrator(gravity=(0,0,-9.8),damping=0.0),
]
O.engines = O.engines + [HarmonicMotionEngine(ids = contenitore, A = (0,0,A), f
= (0,0,fr), label='shaker')]
--
You received this question notification because your team yade-users is
an answer contact for Yade.
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
