yade-users team mailing list archive

[jduriez <question201064@xxxxxxxxxxxxxxxxxxxxx>]

Question #201064 on Yade changed:
https://answers.launchpad.net/yade/+question/201064

jduriez proposed the following answer:
Hello Mukesh,

Nice "Yade exercise" ;-)

Your results are consistent with the way Yade computes things (see
http://bazaar.launchpad.net/~yade-
dev/yade/trunk/view/head:/pkg/dem/NewtonIntegrator.cpp, or documentation
surely)

E.g, you have indeed :
position(t=0.4) = position(t=0.3) + speed(t=0.4) * 0.1, and so on for other times
(velocities are computed first, before modifiying positions)

In my opinion, the discrete (in time !) feature of this integration
scheme is the only reason why you do not get the exact solution. Try
with smaller dt, and you should certainly approach the result you
expect.

-- 
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.