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

I am trying to simulate a transfer chute with a mix of  spherical and polyhedra 
particles.  

When I only have spherical particles the following interaction loop does a 
great job of modelling the flow through the chute:

    InteractionLoop(
        [
            Ig2_Sphere_Polyhedra_ScGeom(),
            Ig2_Sphere_Sphere_ScGeom(),
            Ig2_Facet_Sphere_ScGeom(),
        ],
        [
            Ip2_FrictMat_FrictMat_MindlinPhys(en=Restitution, gamma=Stickiness, 
krot=RollRes),
        ],
        [
            Law2_ScGeom_MindlinPhys_Mindlin(includeMoment=True, 
includeAdhesion=True),
        ],
    ),

When I add polyhedra, I get a bunch off errors like:  <FATAL ERROR> 
InteractionLoop:182 virtual void yade::InteractionLoop::action(): None of given 
Law2 functors can handle interaction #167+2972, types geom:PolyhedraGeom=10 and 
phys:FrictPhys=3 (LawDispatcher::getFunctor2D returned empty functor)

To fix this problem I used the following interaction loop:

InteractionLoop(
        [
            Ig2_Sphere_Polyhedra_ScGeom(),
            Ig2_Sphere_Sphere_ScGeom(),
            Ig2_Polyhedra_Polyhedra_PolyhedraGeom(),
            Ig2_Facet_Polyhedra_PolyhedraGeom(),
            Ig2_Facet_Sphere_ScGeom(),
        ],
        [
            Ip2_FrictMat_FrictMat_MindlinPhys(en=Restitution, gamma=Stickiness, 
krot=RollRes),
            Ip2_FrictMat_PolyhedraMat_FrictPhys(),
            Ip2_PolyhedraMat_PolyhedraMat_PolyhedraPhys(),
        ],
        [
            Law2_ScGeom_MindlinPhys_Mindlin(includeMoment=True, 
includeAdhesion=True),
            Law2_ScGeom_FrictPhys_CundallStrack(),
            Law2_PolyhedraGeom_PolyhedraPhys_Volumetric(),
        ],
    ),

The simulation runs fine but neither polys nor spheres do not appear to have 
any coefficient of restitution after I add polys to the simulation.  Is there a 
ways to force the spheres to use the Mindlin physics?  Ideally I'd like the 
polyhedra to use the mindlin physics as well.  Right now it looks like Mindlin 
Law2 and Ip2 are ignored and everything is defaulting to CundallStrack.

Do I need to modify the C++ code to get polys to have a coefficient of 
restitution or am I missing something here?

-- 
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     : yade-users@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-users
More help   : https://help.launchpad.net/ListHelp

Reply via email to