yade-users team mailing list archive

[Xue <question707124@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Answered => Open

Xue is still having a problem:
Hello jan!
I rechecked the shape of the two particles in the above interaction. As you expected, the contact is not a combination of spherical particles + polyhedral particles.
The above model is a coupling model built with oofem and yade, the id=3841 in the error report exceeds the number of particles I imported via yade, so I deduce that this should be a transformed unit (facet?) on the coupling interface. Therefore, the contact type should be polyhedral particles + transformed unit from FEM.
However, I did not figure out where the contact parameters of the interface unit are defined when it is involved in the contact force calculation (in DEM)? Besides, the above reported error only appears when I defined FrictMat in the script. When I define only PolyhedraMat and use only the following engine, the problem did not occur.

O.engines =[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb(),Bo1_Polyhedra_Aabb()],verletDist=.0005),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom(),Ig2_Facet_Polyhedra_PolyhedraGeom(),Ig2_Polyhedra_Polyhedra_PolyhedraGeom(),Ig2_Sphere_Polyhedra_ScGeom()],
[Ip2_FrictMat_FrictMat_FrictPhys(),Ip2_PolyhedraMat_PolyhedraMat_PolyhedraPhys()],
[Law2_ScGeom_FrictPhys_CundallStrack(),Law2_PolyhedraGeom_PolyhedraPhys_Volumetric()],
),
NewtonIntegrator(damping=0.7,gravity=(0,0,-0.3),label='newton'),
]

I apologize for providing wrong information.
xue

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