Boundary conditions at Z=0

[FredDuhart]

Hi Mostafa,

First of all, thank you very much for this great job.
Your program works very well, but i noticed one little issue. When z = 0 and x <= a then vertical stress should be equal to q.
But we can see in example 'Stress_Validation.py' that Stress_Z[0] = 107.18604438890017 and not 100 (q value).

Have you some explanation for this issue?

Thanks!

Frédéric

[FredDuhart]

I answer to my self :
A way to resolve this problem is to make 2 two succesive iteration (for example iteration =25 and iteration = 26) and the good result is the mean of the 2 results.