yade-users team mailing list archive

[Andrew Jennings <question696253@xxxxxxxxxxxxxxxxxxxxx>]

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.