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