4.5

4.5 | 2 ratings Rate this file 16 Downloads (last 30 days) File Size: 3.48 KB File ID: #25943
image thumbnail

Mogi: point source in elastic half-space

by

 

24 Nov 2009 (Updated )

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   Please login to tag files.
Please login to add a comment or rating.
Comments and Ratings (2)
12 Oct 2011 Alan Rendon  
17 Oct 2010 Forrest Brett  
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