← Back to team overview

yade-users team mailing list archive

Re: Linear law + moment rotation

 

On 23 November 2010 19:09, Bruno Chareyre <bruno.chareyre@xxxxxxxxxxx>wrote:

> On 23/11/10 19:37, chiara modenese wrote:
> >  I want to assign a
> > treshold plastic value for the bending moment (as Luc CR). May I modify
> it?
> > I think those options should be there in the code, anyway.
>
> Sure! Make the threshold very high by default please, so it will not modify
> existing stuff.
> Something like 100π*stiffness should be enough...
>
Ok, I will.

>
> In order to nullify twist moment, you'll need to add a "ktw" attribute in
> CohFrictMat and
> CohFrictPhys, as currently the same "kr" is used for both bending and
> twisting. Then you
> can, set it to a value in Ip2_2xCohFrictMat_CohFrictPhys.cpp, and replace
> kr->ktw in
> CohesiveFrictionalContactLaw.cpp:90.
> That way, you'll be able to set kr=something and ktw=0.
>
Oh yes, I was already thinking to split it into two separate attributes. It
is the easiest and less intrusive thing to do.
Thanks for suggestions!

>
> By the way, this law is used in the preprocessor CohesiveTriaxialTest.
>
Oh yes, actually I found it after a while :-)
cheers, Chiara

>
> B.
>
> _______________________________________________
> Mailing list: https://launchpad.net/~yade-users<https://launchpad.net/%7Eyade-users>
> Post to     : yade-users@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~yade-users<https://launchpad.net/%7Eyade-users>
> More help   : https://help.launchpad.net/ListHelp
>

References