Hello,
I am trying to model the Earth crust. My model is 50 km * 30 km and in 2D (plane strain). I am using the solid-mechanics physics with linear elastic material. Fora first test, I chose E = 40*10^9 Pa, poisson = 0.25, density = 2800 kg/m^3. The left and right boundaries are roller, and the lower boundary is fixed (no displacement). Only the top boundary is free.
There aren't any inital values (no displacement, no velocity). Finally, I have a body load (the pressure; force per unit volume) : -rho*g for value y.
After calculation with Comsol, my problem is:
- The crust isn't static, there is a displacement. Why is that so?
- The stress distribution due to gravity calculed by Comsol differs from the pressure rho*g*y: there is a non-zero deviatoric stress. Why is there a deviatoric stress since the stress distribution should be lithostatic?
Thank you for your help!
Alexis
I am trying to model the Earth crust. My model is 50 km * 30 km and in 2D (plane strain). I am using the solid-mechanics physics with linear elastic material. Fora first test, I chose E = 40*10^9 Pa, poisson = 0.25, density = 2800 kg/m^3. The left and right boundaries are roller, and the lower boundary is fixed (no displacement). Only the top boundary is free.
There aren't any inital values (no displacement, no velocity). Finally, I have a body load (the pressure; force per unit volume) : -rho*g for value y.
After calculation with Comsol, my problem is:
- The crust isn't static, there is a displacement. Why is that so?
- The stress distribution due to gravity calculed by Comsol differs from the pressure rho*g*y: there is a non-zero deviatoric stress. Why is there a deviatoric stress since the stress distribution should be lithostatic?
Thank you for your help!
Alexis