Summary: We propose a new discretization scheme for the elastic stress field. Local ruptures give rise to a stress redistribution which can be represented by double couples. The model is applied to earthquake simulation. Computational efficiency is achieved by the use of lattice Green functions. The model allows inclusion of phenomenological features such as annealing and static fatigue.