← Back to team overview

yade-users team mailing list archive

Re: [Question #696253]: Polyhedra-Sphere-Facet Coefficient of Restitution

 

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

    Status: Open => Answered

Jan Stránský proposed the following answer:
Hello,

If you want all interactions to be of the same type, there is nothing
easier :-)

Something like this could work:
Spheres - FrictMat
Polyhedrons - FrictMat
Facet - PolyedraMat
( there is currently no Ig2_Facet_Polyhedron_ScGeom, although it should not be difficult to implement it, then you could get rid of PolyhedraMat entirely)
and
InteractionLoop(
        [
            Ig2_Sphere_Polyhedra_ScGeom(),
            Ig2_Sphere_Sphere_ScGeom(),
            Ig2_Polyhedra_Polyhedra_ScGeom(),
            Ig2_Facet_Polyhedra_PolyhedraGeom(),
            Ig2_Facet_Sphere_ScGeom(),
        ],
        [
            Ip2_FrictMat_FrictMat_MindlinPhys(en=Restitution, gamma=Stickiness, krot=RollRes),
            Ip2_FrictMat_PolyhedraMat_FrictPhys(),
        ],
        [
            Law2_ScGeom_MindlinPhys_Mindlin(includeMoment=True, includeAdhesion=True),
            Law2_PolyhedraGeom_PolyhedraPhys_Volumetric(),
        ],
    ),

Have a try and let us know

> Ig2_Polyhedra_Polyhedra_ScGeom creates really strange behavior

ScGeom between Polyhedrons is an easy "hack", not too much tested etc.
If you can isolate the problem in a MWE (e.g. 2 specific polyhedrons with given initial velocities), please open a new question or new issue.

> sometimes crashes so that did not work.

Please be more specific (ideally again in an separate issue or question)

Cheers
Jan

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