yade-users team mailing list archive

Thread
Date

[Question #228778]: About considered facet's radius in Ig2_Facet_Sphere_ScGeom.

To: yade-users@xxxxxxxxxxxxxxxxxxx
From: Kneib François <question228778@xxxxxxxxxxxxxxxxxxxxx>
Date: Mon, 13 May 2013 13:56:12 -0000
Reply-to: question228778@xxxxxxxxxxxxxxxxxxxxx
Sender: bounces@xxxxxxxxxxxxx

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

Hi, 

I found something weird in the source code into Ig2_Facet_Sphere_ScGeom::go().
Usually when there is a contact between a sphere and a surface (wall, box ...), the contact geometry contains radius1=radius2=rs. (assuming rs is the sphere radius). So the stiffnesses of the contact are computed like if there is two spheres having the same radius.

In Ig2_Box_Sphere_ScGeom::go :
          scm->radius1 = s->radius;
          scm->radius2 = s->radius;

In Ig2_Wall_Sphere_ScGeom::go :
          ws->radius1=ws->radius2=radius;

But into Ig2_Facet_Sphere_ScGeom::go, there is :
          scm->radius1 = 2*sphereRadius;
          scm->radius2 = sphereRadius;

Don't you think this is a mistake ?

-- 
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.