I'm working on the code atm. Just two remarks for now :
1) //! stretch ratio ?(n) (http://en.wikipedia.org/wiki/Finite_strain_theory#Stretch_ratio)
This comment is wrong, the strectch ratio is |F.x|/|x|, which is a function of non-diagonal terms. The values returned by getStretchRatio are not stretchRatios, they have no clear meaning and are not invariants of F. As mentionned before, out of diagonal terms have an effect on stretch ratios (this is why I don't know how to work with _size without Hsize).
2) The stress update is not done at each step, this is a problem since the size of the period will be controlled based on outdated values.
I think I'll commit to a branch, I have too many changes in view and I don't want to mess with the trunk.
Bruno -- _______________ Chareyre Bruno Maître de Conférences Grenoble INP Laboratoire 3SR - bureau E145 BP 53 - 38041, Grenoble cedex 9 - France Tél : 33 4 56 52 86 21 Fax : 33 4 76 82 70 43 ________________ _______________________________________________ Mailing list: https://launchpad.net/~yade-dev Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-dev More help : https://help.launchpad.net/ListHelp
