Une implémentation massivement parallèle et efficace de la méthode de champ de phase permettant la simulation de l'initiation et de la propagation de fissures dans des structures quasi-fragiles fortement hétérogènes est proposée. Cette implémentation est capable de modéliser des structures avec plusieurs milliards de degrés de liberté, dans le cadre de la méthode des éléments finis utilisant des maillages structurés tridimensionnels, sur des supercalculateurs. Elle est adaptée pour les machines à mémoire distribuée et utilise l'interface MPI pour la répartition de la charge de travail et la communication de données entre les processus. Ce contexte est illustré par la simulation de la fissuration dans un modèle géométrique réaliste d'un béton léger obtenu par la microtomographie (XR-µCT).