search

Mostafa-Nakhaei / PyMastic

Multi-layer Elastic Program in Python
Apache License 2.0
15 stars 7 forks source link

Boundary conditions at Z=0 #4

Open FredDuhart opened 1 year ago

FredDuhart commented 1 year ago

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 commented 1 year ago

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.