New question #698392 on Yade:
https://answers.launchpad.net/yade/+question/698392
Dear YADE developers and users,
I am conducting a simulation where the energy results are important for the
study. So I want to study how the kinetic energy is calculated in YADE. In the
documentation it shows kinetic energy is calculated by E=0.5*m*v^2+0.5*I*w^2,
which is definitely normal and accurate. So I have tried to calculated kinetic
energy by myself, as the particle mass m, velocity v, inertia I and angular
velocity w can be all obtain by using O.bodies[]state.mass,
O.bodies[].state.vel, O.bodies[].state.inertia and O.bodies[].state.angVel. But
when I have set the value of O.cell.velGrad, the E=0.5*m*v^2+0.5*I*w^2
calculated by myself is much larger than utils.kineticEnergy(). And when I
didn't use O.cell.velGrad, the E=0.5*m*v^2+0.5*I*w^2 calculated by myself is
the same with utils.kineticEnergy()
I am thinking why this happened? Why setting the value of O.cell.velGrad will
make E=0.5*m*v^2+0.5*I*w^2 larger than the utils.kineticEnergy()?
My code for calculating the E=0.5*m*v^2+0.5*I*w^2 is:
kineE=0
for i in bodies_id:
m=O.bodies[i].state.mass
v2=numpy.dot(O.bodies[i].state.vel, O.bodies[i].state.vel)
I=O.bodies[i].state.inertia
w2=numpy.dot(O.bodies[i].state.angVel, O.bodies[i].state.angVel)
kineE=kineE+0.5*m*v2+0.5*I*w2
Thank you very much!!
Yuxuan
--
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 : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
