Question #703169 on Yade changed:
https://answers.launchpad.net/yade/+question/703169
Status: Answered => Open
Xu is still having a problem:
Hi Bruno,
Thank you for your reply.
> I noticed that you don't take kinetic energy into account, is it negligible?
Kinetic energy has already been considered in MWE(#9) above. Anyway, kinetic
energy is negligible for not conserving energy.
The microscale energy calculated by:
Law=Law2_ScGeom_FrictPhys_CundallStrack(traceEnergy=True)
Estrain=Law.elasticEnergy()
Efriction=Law.plasticDissipation()
Edamp=O.energy['nonviscDamp']
Ekin=utils.kineticEnergy()
Etot=law.elasticEnergy()-Ep_ini+O.energy['nonviscDamp']-Edamp_ini+law.plasticDissipation()+utils.kineticEnergy()-Ekin_ini,
The input work done by external stress:
dw=(0.5*(O.cell.preVelGrad+O.cell.velGrad)*utils.getStress())*O.dt*(O.cell.hSize).determinant()
W=W+dw
The two parts being the same when using sphere particles . However, the
microscale energy ('Etot') is greater than input work when using clump
particles. Other details can be seen in MWE#9.
Cheers,
Xu
--
You received this question notification because your team yade-users is
an answer contact for Yade.
_______________________________________________
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
