Recently, I wanted to do a simple test of the stress distribution within the lithosphere. I give the viscosity in advance, I give the force by Neumann boundary conditions, and I look at the effect of viscosity on the force component.
The results were as I expected: higher stresses at the top of the viscous upper mantle and lower stresses at the bottom of the lower crust.
But at the same time I have some questions:
- I found that the simulation results did not change during the simulation time.
- In the crust, the obvious viscosity difference between the bottom of the upper crust and the top of the lower crust should also produce obvious changes in stress, but the results show that the stress in the crust is gradually reduced.
Does anyone have a similar problem or give me some advice?
gravity_infstrain.cfg (2.4 KB)
pylithapp.cfg (7.6 KB)