> so, how to indentify which particle should be calculated and update > first , and which should be later? > (because it is a dynamic process, so just use the order of ID is > not a way to solve this) The trick is that you have constant positions during the whole timestep, you calculate forces from all contacts (order doesn't matter, elements don't move during calculation) and then, before next timestep begins, there is LeapFrogPositionIntegrator, that calculates displacements based on velocity, inertia, applied forces, ... and moves elements to the new position. Obviously, your timestep must be small enough so that you get convergence.
HTH, Vaclav _______________________________________________ Yade-users mailing list [email protected] https://lists.berlios.de/mailman/listinfo/yade-users
