[Getfem-users] problem with dirichlet boundary conditions

[Eugen Wintersberger]

Hi there
 I try to solve a very simple problem with getfem: a bended bar. 
In fact I consider a bar fixed on one side and exerted to a force (in
negative z direction) on the other one. Clearly on the fixed end I have
Dirichlet boundary conditios with \vec{u}=0 and some force on the other
side. Therefore, I expect that the maximum value in uz would be zero
(due to the Dirichlet boundary condition).
However, in the solution I get the following:
...
solve done!
get displacement field
ux min: -1.090228e-08, max: 1.090042e-08
uy min: -4.643458e-07, max: 4.642960e-07
uz min: -5.047865e-06, max: 1.396889e-09
...
Obviously, the maximum in uz direction is not zero. Attached to this
mail you will find the small python program performing the calculation. 
Is there any obvious mistake I did in the program?

best regards 
 Eugen Wintersberger

PS: I'm quite new to getfem and to FEM in general ;) so I hope my
question is not too stupid ;))

simplebar.py
Description: Text Data