← Back to team overview

yade-users team mailing list archive

Re: [Question #703169]: Energy conservation of clump

 

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.