Question #234129 on Yade changed:
https://answers.launchpad.net/yade/+question/234129
Status: Answered => Open
lingran is still having a problem:
Hi Bruno,
Thanks for your comments. You are right and I will try my best to add
law2.bendingElastEnergy() in c++ code.
But I still have several doubts.
1. like Jan said, Why c++ is one time step earlier than python?
2. When use Law2_ScGeom6D_CohFrictPhys_CohesionMoment() and CohFrictMat,
why do we need to set both alphaKr and alphaKtw to a non-zero value, otherwise
the simulation would crush?
3. Energy is not balanced even without considering friction and moments,
that means: utils.kineticEnergy+ Emgh(sum of
(i.state.mass*(i.state.pos[1]*9.81 ))+
law2.normElastEnergy()+law2.shearElastEnergy()+O.energy[‘nonviscDamp’]!=constant..
The error increases with the increasing number of particles.
Thanks.
Best regards.
Lingran
--
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
