yade-users team mailing list archive

[Bruno Chareyre <bruno.chareyre@xxxxxxxxxxx>]

Ok, thanks. I confirmed that with tests in the meantime. I'm a bit 
surprised, I though the plastic part was removed from total displacement 
(or ref point moved, which is the same). Consequence : it needs to use 
the incremental form (currently bugged with periodic, working on that 
now) or adapt the total one.

FYI : the volumic collapse of the periodic cell we saw in last days 
appears randomly. I got alternatively correct/stupid results with 
exactly the same script. Still investigating that.

Cheers.

Bruno


Václav Šmilauer a écrit :
> > I'm trying to understand strange results in periodic triax (sudden 
> > increase of stress with an increased friction).
> > Quick question : what happens in the total formulation 
> > (Ig2_Dem3D+Ip2_Frict+LawDem3_FrictPhys) when the friction coefficient is 
> > modified?
> > Could it modify the computed force instantly somehow?
> >     
> If you modify friction of the interaction (FrictPhys::ks), you will get
> immediately different shear force (as total u_T will be still the same):
>
> pkg/dem/Engine/GlobalEngine/ElasticContactLaw.cpp:105
>
> 	Real maxFsSq=phys->normalForce.SquaredLength()*pow(phys->tangensOfFrictionAngle,2);
> 	Vector3r trialFs=phys->ks*geom->displacementT();
>
> v
>
>
>
>
> _______________________________________________
> Mailing list: https://launchpad.net/~yade-users
> Post to     : yade-users@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~yade-users
> More help   : https://help.launchpad.net/ListHelp
>
>   


--
_______________
Bruno Chareyre
Associate Professor
Grenoble INP
Lab. 3SR
BP 53 - 38041, Grenoble cedex 9 - France
Tél : 33 4 56 52 86 21
Fax : 33 4 76 82 70 43
________________