yade-users team mailing list archive

["Sergei D." <sj2001@xxxxxxxxx>]

25.03.2010 18:57, chiara modenese:
> I am testing my script with the class you committed but there is something
> wrong. In the contact law I used to "play" with the damping in the shear
> direction I coded something as in the attached file. Please look at that (I
> changed the terminology but it is clear what it does) so you can check if
> you are doing well.
> I do not commit my law simply because it is quite a mess and I was playing
> with it to clarify this issue around the damping. Please update your law
> following the procedure as I suggested.
> 
> Thanks, Chia
> 


Sorry, but I'm confused... You are wrote in the ViscoGlobal.cpp:

> ///////////////////////////////////////// MOHR-COULOMB LAW
>	Real maxFs =
phys->normalForce.SquaredLength()*std::pow(phys->tangensOfFrictionAngle,2);
>	if (shearElastic_cum.SquaredLength() > maxFs)
>	{
>	Real ratio = Mathr::Sqrt(maxFs)/shearElastic_cum.Length();
>	shearElastic_cum *= ratio;
>	phys->shearForce = shearElastic_cum;
>	}
>	else
>	{phys->shearForce = shearElastic_cum + shearViscous;} // ad

So, you have only elastic friction for dynamic, and (elastic+viscous)
for static friction. Is'n it?


-- 
Best regards,
Sergei D.