Cette thèse concerne la modélisation d'ouvrages maçonnés tridimensionnels en grand appareil, précontraints ou comportant un grand nombre de blocs. Ces structures sont des constructions en gros blocs de pierre, éventuellement liés par des joints de mortier, dont les géométries tridimensionnelles ne permettent généralement pas la définition d'un volume élémentaire représentatif. Le comportement mécanique de ces édifices ou structures, peut être étudié par des méthodes qui considèrent le comportement des blocs et leurs interactions avec des blocs voisins, par frottement sec ou en présence de mortier cohésif. Nous avons donc choisi de les modéliser par la méthode des éléments discrets Non Smooth Contact Dynamics. Les monuments étudiés sont des structures en pierres massives, simple appareil régulier (mur) ou irrégulier (pont Julien) et structures précontraintes (escalier, coupôle). D'une part, un nouvel algorithme de détection des contacts a été développé pour gérer les contacts surfaciques, spécifiques aux maçonneries. Diverses applications sur des cas réels, et des comparaisons avec des expériences issues de la littérature, ont montré la robustesse de l'algorithme. Les résultats obtenus par l'analyse des descentes de charges sont cohérents avec les caractères physiques et géométriques des ouvrages considérés. D'autre part, les liaisons de contact conduisent souvent à un problème statiquement indéterminé au niveau local. Nous avons proposé une méthode de résolution qui traite simulatnément tous les noeuds d'une liaison plane. L'étude sur des cas académiques a permis de tester la pertinence de la méthode. Une comparaison a été faite avec la méthode de résolution existante, implémentée dans le code LMGC90.