> 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. _______________________________________________ Mailing list: https://launchpad.net/~yade-dev Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-dev More help : https://help.launchpad.net/ListHelp
