← Back to team overview

yade-users team mailing list archive

[Question #694974]: Proper Law2 for Mindlin Physics w/ Poly, Sphere and Facets

 

New question #694974 on Yade:
https://answers.launchpad.net/yade/+question/694974

I am trying to create a simulation with polyhedra, spheres, and facet.  The facets are a stationary geometry generated using ymport.stl().  Materials are:

matSpheres = O.materials.append(ViscElMat(
	young=1.0e8, 
	frictionAngle=radians(fAng), 
	density=2160, 
	poisson=0.3, 
    tc = tcIn,
))

matPoly = PolyhedraMat(
	young=1.0e8, 
	frictionAngle=radians(fAng), 
	density=2160, 
	poisson=0.3,
)

matWalls = O.materials.append(ViscElMat(
	young=1.0e8,
	poisson=0.3,
	frictionAngle=radians(fAng), 
    tc = tcIn,
))

I tried the following engines:

InsertionSortCollider([Bo1_Polyhedra_Aabb(),Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()],verletDist=(2.0e-3),label='collider',ompThreads=1),
	InteractionLoop(
		[Ig2_Sphere_Sphere_ScGeom(), Ig2_Facet_Sphere_ScGeom(), Ig2_Sphere_Polyhedra_ScGeom(), Ig2_Facet_Polyhedra_PolyhedraGeom(), Ig2_Polyhedra_Polyhedra_PolyhedraGeom()],
		[Ip2_FrictMat_FrictMat_MindlinPhys(en=0.3, gamma=0.4, krot=0.3), Ip2_FrictMat_PolyhedraMat_FrictPhys()],
		[Law2_ScGeom_MindlinPhys_Mindlin(includeMoment=True, includeAdhesion=False), Law2_ScGeom_FrictPhys_CundallStrack()],
	),

yadedaily says: "python3.7: /builds/yade-dev/trunk/deb/yadedaily/pkg/dem/FrictPhys.cpp:19: virtual void yade::Ip2_FrictMat_FrictMat_FrictPhys::go(const boost::shared_ptr<yade::Material>&, const boost::shared_ptr<yade::Material>&, const boost::shared_ptr<yade::Interaction>&): Assertion `dynamic_cast<GenericSpheresContact*>(interaction->geom.get())' failed.
Aborted"

yade says: "FATAL /build/yade-KdQLrC/yade-2019.01a/pkg/common/InteractionLoop.cpp:143 action: None of given Law2 functors can handle interaction #1716+4252, types geom:PolyhedraGeom=10 and phys:FrictPhys=3 (LawDispatcher::getFunctor2D returned empty functor)"

so I suppose I need a different Law2 but I have tried a number of them and do not see one that works.  Can Yade simulate a mix of spheres and polys with adhesion between all particles and surface and with rotating resistance on the spheres?  If so, what combination of materials and interactions does one use?

 

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.