← Back to team overview

yade-dev team mailing list archive

Re: Triaxial total energy

 

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

References