Hello, I was thinking about implementation of stressTensorOfPeriodicCell functions (as discussed at https://blueprints.launchpad.net/yade/+spec/stressgetters ), and I would suggest this: - create separate functions for normal and shear contribution (as it is in utils.normalShearStressTensors) - other functions can be combined from this function (total stress vector, symmetrized stress vector etc.)
for the implementation applicable to both Dem3dofGeom and SCGeom, it would be nice if the IGeom class had giveLength() method. It would be in GenericSphereContact class and each inherited class would overload it in its specific way. Would you agree with this approach? I was also thinking about computation of (secant) stiffness matrix of periodic cell, the formulas for estimation (even though pretty good according to my recent tests) are very similar as for stress computation. The only problem is that current normal and shear stiffness of the interaction is needed (and in general it can differ from IPhys.kn or IPhys.ks as each material model treat them in different way). The solution would be very similar to the stress one, i.e. each IPhys class would have giveCurrentKn() and giveCurrentKs() methods (which would return kn and ks by default). Again, would you agree with such implementation? Thanks for any feedback and have a nice day Jan _______________________________________________ Mailing list: https://launchpad.net/~yade-dev Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-dev More help : https://help.launchpad.net/ListHelp
