I think I have something to share. How would you feel about co-authoring a paper to this conference - mainly to further expose SfePy but also to check that I'm not applying it incorrectly?

Using the strain rate dependency example you provided I recoded an example looking at the failure of a cylindrical specimen under compression. The cylinder is made up of elements that I treat individually or discretely, checking the strain rates in each and using a relationship established in the lab, I calculate the compressive strengths in the elements based on these strain rates. If the stress in an element exceeds its strength then the element fails and I assign a low stiffness to that element. The code I'm using is available here:

Would you mind looking and commenting on this code. I'm a bit skeptical about the avgqp function which averages the strains at the quad points. Note that I use a normal distribution to vary the moduli of the elements initially.  
I've attached the mesh: cylgeo.vtk.

What I aim to do with the paper is demonstrate the influence of specimen geometry on strength. For simple compression and tension tests (on cylinders) this is trivial but with indirect tensile tests, for example, there is both compression and tension failure, and I hope to show using FEM that the strengths of materials determined using these tests are significantly influenced by the geometry of the specimen.     

The force-displacement curves currently generated by the code look promising. I'm not sure the analysis can be refined by using smaller elements?


On Fri, Jan 14, 2011 at 3:19 AM, Robert Cimrman <> wrote:
I see :) Good luck then!


On Thu, 13 Jan 2011, Andre Smit wrote:

No plans currently - still working on an analysis. Maybe if I can get it
done before the end of the month :)

On Thu, Jan 13, 2011 at 2:43 AM, Robert Cimrman <>
     Hi Andre,

     are you planning to go there?


On Mon, 10 Jan 2011, Andre Smit wrote:



You received this message because you are subscribed to the Google Groups "sfepy-devel" group.
To post to this group, send email to
To unsubscribe from this group, send email to
For more options, visit this group at