Mogi: point source in elastic half-space
by François Beauducel
24 Nov 2009
(Updated 08 Nov 2011)
Computes surface displacements, tilt and strain due to point source at depth
|
Watch this File
|
| File Information |
| Description |
The Mogi [1958] model allows to calculate analytic solution for surface deformation due to a point source in an elastic half-space. This model is widely used to simulate ground deformation produced by local perturbation like volcanic magma chamber. It computes displacements, tilt and strain in a polar space, due to a volume variation at depth [Anderson, 1936] or an isotropic pressure variation in a spherical source [Mogi, 1958].
The proposed Matlab script is a literal transcription of Mogi's simple equations, extended to non-isotropic medium (Poisson's ratio different from 0.25). All parameters can be vectorized. See help for syntax, and script comments for details. |
| MATLAB release |
MATLAB 5.2 (R10)
|
|
Tags for This File
|
| Everyone's Tags |
|
| Tags I've Applied |
|
| Add New Tags |
Please login to tag files.
|
| Updates |
| 14 Dec 2009 |
Possibility to specify a Poisson's ratio different from 0.25 for non-isotropic medium. Previous input parameters are still allowed. |
| 05 Jan 2010 |
IMPORTANT: corrects an error when specifying elastic parameters E and nu.
Adds some checks on input arguments.
Adds a scheme explaining the model parameters. |
| 08 Nov 2011 |
Bug correction: radial strain value (optional er output argument) was wrong... Now corrected. |
|
Contact us at files@mathworks.com
