Cette thèse aborde le problème de la rupture localisée dans les roches, qui caractérise un grand nombre d'applications dans le domaine du génie civil, tels que la rupture du barrage, effondrement des fondations, la stabilité des excavations ou les tunnels, les glissements de terrain et les éboulements. Le risque de rupture localisée devrait être mieux appréhendé pour mieux l'éviter. La rupture localisée dans les roches est généralement caractérisée par une une rupture soudaine et quasi-fragile sans avertissement sous forme de grandes déformations et visibles avant la rupture elle-même. Cela se produit également sous l'influence des hétérogénéités matériels, des fissures existantes et d'autres défauts initiaux.Les trois nouveaux modèles numériques, intégrant les mécanismes de ruptures localisées, l'hétérogénéité de la roche et des fissures existantes et d'autres défauts, sont présentés dans cette thèse. Le premier modèle propose une représentation en 2D de roche composite à deux phases, où la phase solide représente la roche intacte et la faible phase indique les défauts initiaux. Le deuxième modèle représente l'extension du modèle précédent vers un espace 3D, où est considéré un ensemble complet de mécanismes de ruptures en 3D. Les propriétés hétérogènes sont considérées ici par une distribution aléatoire en accord avec la variation statistique de Gausse. Ce modèle est également utilisé pour l'analyse de la roche intacte par spécimens possédant des écarts de formes géométriques qui influencent la résistance à la compression uni-axiale. Le troisième modèle est un modèle en 2D, traitant l’interaction volumétrique entre un fluide et la structure sous l’influence de l’écoulement du fluide à travers le milieu de la roche poreuse.L'approche des lattices discrètes est choisie pour construire un cadre général pour les trois modèles, où les cellules de Voronoï représentent les grains de roche tenus ensemble par les poutres de Timoshenko comme des liens de cohésion. La cinématique améliorée est caractérisée par l'approche intégrée des discontinuités comme un supplément à la cinématique standard de liens cohérents. Cela sert pour la propagation de la macro fissure dans tous les modes de ruptures et de leurs combinaisons, entre les grains de la roche. La formation de la zone du processus de rupture suivie par des microfissures coalescentes, précédant la rupture localisée, est aussi considérée dans les modèles. L’écoulement du fluide est régi par la loi de Darcy, tandis que les conditions de couplage obéissent à la théorie de poroélasticité de Biot.Les résultats des modèles numériques ont été vérifiés par des exemples de la littérature dans le cas des modèles en 2D. Le modèle en 3D a été validé suite aux résultats expérimentaux effectués sur 90 échantillons de roches, où sont considérées de légères déviations géométriques des spécimens.La présentation de ces modèles, ainsi que leurs aspects de mise en œuvre sont présentés en détail. L’approche avec une discontinuité intrinsèque et le caractère local des améliorations nécessaires à la simulation des discontinuités de déplacement orientent vers la condensation statique des degrés de liberté améliorés sont efficacement intégrés dans l’architecture des éléments finis.