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 <jiaoliq...@163.com> > 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@lists.launchpad.net > Unsubscribe : > https://launchpad.net/~yade-users<https://launchpad.net/%7Eyade-users> > More help : https://help.launchpad.net/ListHelp > >
_______________________________________________ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp