Ok, it was only a sign problem. :-)
It seems your conclusion on bad conservation when out of equilibrium is
correct (consistent with your previous findings too).
You'll see that I initialized total energy at 0 at the begining of
compression.
BTW, Q. : how do you initialize to zero a variable you were storing in
plotData()? How to do it in python?
plot.addData(i=O.iter,e22=-triax.strain[1]-refstrain2,...)
With refstrain2=0 initialy, and triggering later
refstrain2=-triax.strain[1].
Bruno
_______________________________________________
Mailing list: https://launchpad.net/~yade-dev
Post to : [email protected]
Unsubscribe : https://launchpad.net/~yade-dev
More help : https://help.launchpad.net/ListHelp
