Nous modelisons les processus de glissement instable, en tenant compte de la geometrie souvent complexe des systemes de failles et des lois de frottement observees en laboratoire. Nous assimilons les failles a des interfaces de discontinuite dans un milieu elastique parfait, la croute, et utilisons une loi de frottement dependant du glissement deduite d'experiences qui montrent que la resistance du materiau diminue du seuil statique au seuil dynamique, proportionnellement au glissement. <br /><br />Nous proposons deux methodes numeriques. La premiere a pour but de simuler l'evolution temporelle spontanee d'un reseau de faille, soumis initialement a un champ de contraintes donne et auquel nous appliquons une perturbation initiale en vitesse. Elle utilise un schema de type Newmark en temps, et une discretisation spatiale en elements finis avec decomposition de domaine. Elle se revele capable de capturer efficacement les instabilites du glissement, et en particulier la phase d'initiation, qui precede la propagation de la rupture dynamique et qui se caracterise par une forme auto-similaire et une croissance exponentielle du glissement au cours du temps. Des experiences numeriques montrent que l'interaction, sur des segments de faille presentant un recouvrement significatif, se manifeste par l'existence de « zones d'ombre » dans lesquelles les contraintes sont dechargees et le glissement inhibe. En cas de recouvrement important de deux segments de faille, on observe une dissymetrie des profils de glissement, correspondant a la disparition de la singularite de contraintes a l'une des pointes de faille.<br /><br />Le deuxieme schema numerique realise l'analyse spectrale non lineaire du probleme de l'initiation « pseudo-linearise » autour de la position d'équilibre dans laquelle le système de failles est initialement au seuil de resistance statique. La non-linearite du probleme provient de la prise en compte des zones d'ombre dont on ne connait pas la geometrie a priori. Cette analyse permet de trouver le mode qui porte la signature de l'initiation, c'est-a-dire la forme auto-similaire remarquee plus haut. La version statique de cette analyse modale fournit un critere de stabilite des reseaux de failles, c'est-a-dire la valeur limite du taux d'affaiblissement au-dela de laquelle un episode de glissement donnera lieu a un evenement sismique.<br /><br />Nous faisons l'hypothese que le mode non lineaire statique, qui caracterise un comportement en affaiblissement a la limite de la stabilite, peut etre utilise pour decrire le glissement cumule a l'echelle tectonique sur un reseau de failles normales particulier en Afar, dont nous connaissons les glissements mesures en surface. Nous montrons qu'un choix judicieux du profil d'affaiblissement « equivalent » a l'echelle tectonique permet un bon accord entre le glissement observe et le mode statique. Nous en tirons des conclusions en termes d'interaction, mais aussi de propagation et/ou branchement des segments de failles. Enfin, nous decrivons deux autres applications : <br /><br />1) l'influence de la fracturation secondaire (endommagement) sur la forme du glissement<br />2) les parametres geometriques favorables a l'apparition d'une zone de relai entre deux segments de faille se propageant l'un vers l'autre.