Abstract: | Earthquakes are recognized as resulting from a stick–slip frictional instability along faults. Based on the node‐to‐point contact element strategy (an arbitrarily shaped contact element strategy applied with the static‐explicit algorithm for modelling non‐linear frictional contact problems proposed by authors), a finite element code for modelling the 3‐D non‐linear friction contact between deformable bodies has been developed and extended here to analyse the non‐linear stick–slip frictional instability between deformable rocks with a rate‐ and state‐dependent friction law. A typical fault bend model is taken as an application example to be analysed here. The variations of the normal contact force, the frictional force, the transition of stick–slip instable state and the related relative slip velocity along the fault between the deformable rocks and the stress evolution in the total bodies during the different stages are investigated, respectively. The calculated results demonstrate the usefulness of this code for simulating the non‐linear frictional instability between deformable rocks. Copyright © 2003 John Wiley & Sons, Ltd. |