Hello Vaclav, You are right. I see the force is nan indeed. Thank you so much. Thank you for so warm help. I will check the contact law carefully. Thank you again.
Best Liqing At 2010-09-20 17:25:06，"Vaclav Smilauer" <[email protected]> wrote: >> I am sorry to tell you maybe there are some problems if the rolling >> resistance is considered in Law2_ScGeom_CohFrictPhys_ElasticPlastic. >> Yade give the error and stop the calculation. I thought maybe it is >> derived from the velositybins. But, I still cannot get the force value >> even if I deleted the all parameters in InsertionSortCollider. >> The mentioned error is described as follow: >> 55217 FATAL >> yade.VelocityBins >> /home/jiaoliqing/yade/pkg/common/DataClass/VelocityBins.cpp:77 setBins: Body >> #12179, velSq=-nan not put in any bin!? >> Aborted >> Could you help me to check it? Thank you so much. >That error is coming almost certainly from the contact law (unhandled >division by zero or such); you can try saving the simulation just before >the problem occurs (it should be always at the same place, if you run >the same simulation) and the sub-stepping through the iteration using >the inspector, you will see that the force will be nan after >InteractionDispatchers. > >Cheers, v. > > >_______________________________________________ >Mailing list: https://launchpad.net/~yade-users >Post to : [email protected] >Unsubscribe : https://launchpad.net/~yade-users >More help : https://help.launchpad.net/ListHelp
_______________________________________________ Mailing list: https://launchpad.net/~yade-users Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
