Hi Alexander,

the only problem I see is a wrong shape of the load value returned by linear_tension() - replace:

val = nm.concatenate((fx,fy,fz),axis=2)

with

val = nm.concatenate((fx,fy,fz),axis=1)

After this change both the attached examples work for me.

What is you sfepy version? I strongly suggest updating to the latest version (release or git), as that would have warned you about the wrong shape.

r.

On 10/21/2016 10:49 AM, Alexander Stuhl wrote:

Hi,

I am trying to run a programm concerning bending of a beam by terminal couples in a - more or less - complicated geometry. But SfePy calculated that there is no displacement in the entire volume. Then I reduced - step by step - the complexity of the programm in order to find the error why there is no displacement, but I cannot find it. Finally I made a simple programm concerning bending and compared it to a modified example ( linear_elasticity/linear_elastic_tractions.py ). My programm still not shows a displacement of the volume whereas the modified example does. These programms are identical regarding the code except the orientation of the beam.

probe.mesh == block.mesh: Cylinder from x=-5 to x=5. Hexahedra. proby.py is a modified version of linear_elasticity/linear_elastic_tractions.py. torquefg.mesh written in Gmsh: Same cylinder, but from z=-5 to z=5. Tetrahedra

I do not know the difference. Could you tell me where the fundamental error is?

Best regards Alexander