Question #295301 on Yade changed:
https://answers.launchpad.net/yade/+question/295301
Jan Stránský proposed the following answer:
>
>
>
> if(asymmetry!=1){
> for(size_t i=0; i<negIds.size(); i++){
> negCoords[i]-=dAX;
> axisVel(negIds[i]) = -dAX/scene->dt; // update
> current position
>
> }
>
> }
>
>
there is no computation of displacement here.. just each **coordinate**
(from 0 to negIds.size()) is updated by displacement dAX, so .. Really try
to understand what the code does before being inspired, otherwise it might
be "dangerous" :-)
dAX (or maximum 2*dAX) is the actual change in length of the specimen. If
you compute external work like this (multiplying it by number of fixed
particles) and it equals the elastic energy, there must be a problem
somewhere..
Jan
--
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 : yade-users@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
