We developed a model for a fault in which various areas of the fault phase have different stress-slip constitute laws. The model is conceptually simple, involving nonlinear algebraic equations which can easily be solved by a graphical method of successive iterations. Application is made to the problem of explaining seismicity patterns associated with great earthquakes. The model quantitatively explains phenomena associated with seismic gaps, asperities, and barriers.