geostatic stresses

[dm88]

Dear Code_Bright users,

I would like to simulate geostatic stresses in an oedometer sample (0.056m x 0.02 m).
What I expect at the bottom of the model, is a vertical (total) stress equal to: Syy=(2700 kg/m3 x 9.81 m/s2 x 0.02 m)= 0.52974 kPa=0.00052974 MPa, where 2700 is the solid phase density.

However, code_bright calculates at the bottom of the model a Syy equal to 0.000368 MPa, which is like adopting a weight per unit volume equal to 18 kN/m3, instead of 26 kN/m3 (which I want, considering a density of 2700 kg/m3)

-How could I solve this issue?
-Why Syy profile along Y is not perfectly linear, but change slope at the end of the profile?

Thank you in advance for your reply.

Kind regards,

Dario

[code_bright]

Dear Dario,

Note that the density you are using corresponds only to the solid phase, thus the total density of the porous media will be calculated taking into account the porosity and the degree of saturation.

Regards,

Alfonso