Wolfram Language Paclet Repository
Community-contributed installable additions to the Wolfram Language
MalcolmWoodruff`TerraPercolatio`
PorePressureField  |
|
| | ||||
▪
Below the water table the pore pressure is the FEM value u = 9810 (h - y) Pa; above it, it is capped hydrostatic suction (negative), with the cap set in metres of head by "SuctionCap".
▪
Returns ["OutsideSoil"] outside the soil body, so mapping it over a grid and dropping missing values gives a pore-pressure field ready to export, e.g. for a limit-equilibrium slope-stability program such as TerraLapsus.
▪
For a confined solution (no phreatic surface) it returns the plain FEM pore pressure everywhere on the mesh.
▪
The following options can be given:
"SuctionCap" | 1. | Positive metres of head at which negative pore pressure (suction) above the water table is clipped; the default of 1 m is about 9.81 kPa. |
Examples
(2)
Basic Examples
(1)
Get the pore pressure at a point below the water table:
In[1]:=
damSpec=Association["Profile"{{0.`,0.`},{52.`,0.`},{28.`,12.`},{24.`,12.`}},"UpstreamLevel"10.`,"DownstreamLevel"0.`,"kx"1.`*^-6,"ky"1.`*^-6];damSol=
[damSpec,"Points"16];u=
[damSol];u[26.`,6.`]
 | SeepageUnconfined |
| PorePressureField |
Out[1]=
23.5794
Options
(1)
""