← Back to team overview

yade-users team mailing list archive

Re: [Question #242646]: Changing The Integration Algorithm

 

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

    Status: Answered => Open

burak er is still having a problem:
>Thank you very much for contributing. Do you plan to extend to
rotational DOFs?

>You mention advantages of the adaptative scheme over the symplectic
Euler scheme (NewtonIntegrator). Do you have >benchmarks or paper about
this specific point?

Yes, I think extending also for rotational coordinates. Currently, it is
possible to add rotational coordinates easily but as you know it is a
bit different than the translational coordinates(rotations are not
vector quantities!). I have found a paper here[1] that discusses the
integration of the rotational equations of motion. Fortunately, it
requires only a little modification.

About your question, I do not have benchmark or paper about this. But my
point of view is; although, symplectic Euler scheme might be fine for
the energy preserving simulation of discrete elements, but  I think, its
implementation,  namely, the NewtonIntegrator does not provide well
usability. The most important part is the lack of adaptive integration
where the step size is unpredictable for my cases.

>If you want to share your work, please, commit it in your branch or
>just attach patches here, so we will be able to review it and upload
>to the development branch.

I plan to commit it after I finish the implementation and check
correctness in a few days.

burak

[1]drugg.fgg.uni-lj.si/3947/1/P2011AES.pdf

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