yade-users team mailing list archive

[luc scholtes <lscholtes63@xxxxxxxxx>]

Hi Jiao,

Thank you for your email. I think you are right, it should be written

ratio=MomentMax/scalarMoment;

instead of

ratio *= MomentMax/scalarMoment; // to fix the moment to its yielding value

It is always good to have someone behind your shoulders to avoid this kind
of error ;-)

Feel free to commit the change. Otherwise I will try to do it ASAP.

  Luc

2010/11/17 jiaoliqing <jiaoliqing@xxxxxxx>

> Hi Luc Scholtes,
>
> I have read the code in CohesiveFrictionalPM.cpp written by you. I read
> them carefully, so I found there would be a bug in the file. In the part of
> description of moment, this part of code is as followed:
> /*Plastic moment */
>     if(scalarMoment > MomentMax)
>     {
>       Real ratio = 0;
>       ratio *= MomentMax/scalarMoment; // to fix the moment to its yielding
> value
>       moment *= ratio;
>       moment_twist *=  ratio;
>       moment_bending *= ratio;
>     }
> Furthermore, the same code is in CohesionlessMomentRotation.cpp to
> calculate the moment when the moment is over the plastic limit. It is
> distinct bug that is " ratio *= MomentMax/scalarMoment;".  it is written as
> this, then the ratio will still be zero during the simulation. Then the
> moment will be zero in the plastic stage. That is different from the
> algorithms. So, this one should be changed to "ratio =
> MomentMax/scalarMoment;", others don't need to modify. Do you agree with me?
>
> Look forward your reply. thank you!
>
> Regards
> Liqing
>
>
>
> ------------------------------
> 网易163/126邮箱百分百兼容iphone ipad邮件收发<http://help.163.com/10/1108/18/6L03I4G500753VB8.html>
> _______________________________________________
> 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
>
>