yade-dev team mailing list archive

[Václav Šmilauer <eudoxos@xxxxxxxx>]

> detectionFactor doesn't enter any equation 
> relative to contact geometry itself or contact law.
OK, I looks like the real cause for that.
> As a consequence, this is perhaps not what you want :
> 
> Ig2_Sphere_Sphere_ScGeom.cpp :
> Real penetrationDepth=s1->radius+s2->radius-distance;
> scm->contactPoint=se31.position+(s1->radius-0.5*penetrationDepth)*normal;
Ah, that's right, that is not what I want.

> Later used in CpmLaw like this :
> applyForceAtContactPoint(BC->normalForce+BC->shearForce,contGeom->contactPoint,I->getId1(),pos1,I->getId2(),pos2);
I will recompute the contact point by hand then.

> I'm also wondering if avoidGranularRatcheting=true is not a problem in 
> your case, since it means rad1/rad2 will be used whatever the distance.
> I added avoidGranularRatcheting attribute in Ig2_Sphere_Sphere_ScGeom to 
> let you choose.
ok, thanks.

> Or maybe you are changing refRadius in the interaction, leading to 
> refRad1+refRad2=initialDistance? In that case it should work even with 
> avoidGranularRatcheting=true.
Good point, I don't do that although I should. I will add it at the
beginning of Law2 when the contact isFresh().

I let you know about the result.

Cheers, v.