readme

这是一个模拟第3类模式地震波的matlab脚本。 This a collection of Matlab scripts that solve the antiplane (mode III)

Last update: Thu Sep 29 18:27:00 CEST 2005
Jean Paul Ampuero 	[e-mail]

This is a collection of Matlab scripts that solve the antiplane 
(mode III) earthquake dynamic problem with slip-weakening friction, 
on a 1D fault embedded in a 2D homogeneous elastic unbounded medium.
The problem is formulated as a boundary integral equation
and the elastodynamic kernels are analytically derived in
the spectral domain (spatial wavenumber).
The method is explained e.g. by Morrysey and Geubelle (1997),
and has been improved and extensively used by Nadia Lapusta, 
Alain Cochard, etc.

Time truncation of kernel convolutions is implemented.
Both the "displacement" and the "velocity" formulations of the
interaction term are illustrated, although I most extensively
use the latter.

The friction law used here is slip weakening. The method is
however largely independent of the nature of the friction law.

This directory contains:

  fulldynD.m	"displacement" formulation
  fulldynV.m	"velocity" formulation
  C3.m W3.m	functions used to compute the kernels
  friction.m	linear slip weakening friction law
  ShortPlotDyn.m plots some results (time series)

This scripts are intended to introduce new users to the method.
For serious simulations you should turn to a Fortran 
optimized implementation, such as:
+ BIMAT-PCSI, http://www.sg.geophys.ethz.ch/geodynamics/ampuero/software.html
+ MDSBI, http://www.people.fas.harvard.edu/%7Eedunham/codes.html