This is a classical example to demonstrate the effect of hydromechanical coupling in a poroelastic medium. For more details we refer to a textbook , in which also the analytical solution is derived. As domain we consider a sphere, by symmetry we need to simulate only an octant.
The boundary conditions of hydraulics are atmospheric pressure on the sphere surface and impermeable elsewhere. The boundary conditions of mechanics are an overburden (Neumann) applied as step load on the sphere surface at initial time $t=0$. The remaining sides are fixed in normal direction (Dirichlet).
The material is isotropic, linear elastic. Solid and fluid are assumed to be incompressible. In its initial state the sphere is not deformed and there is ambient pressure everywhere. A sudden load increase on the surface is instantly transferred on the pore pressure, whereas the solid needs time to deform, until it carries the load. Since the pore fluid is squeezed out of the outer layers first, they act like a tightening belt and consequently the pressure in the center rises, it may even exceed the value of the applied load. Finally the pore pressure approaches to ambient pressure.
All parameters are concluded in the following tables.
|Young’s modulus (bulk)||$10\cdot 10^6$||Pa|
|Poisson ratio (bulk)||$0.1$||-|
|Solid density||$2.5\cdot 10^3$||kg/m$^3$|
|Finite Elements||$8741$||Taylor-Hood tetrahedral elements|
|Coupling scheme parameter||$0.7774$||-|
As predicted, the pressure in the center exceeds the applied load and then levels out.
Pressure at center of sphere
For more details about the staggered scheme we refer to the user guide - conventions.
Verruijt, A. (2009): An introduction to soil dynamics. Springer Science and Business Media, DOI:https://doi.org/10.1007/978-90-481-3441-0 https://link.springer.com/book/10.1007/978-90-481-3441-0