On présente ici un modèle numérique de fracturation hydraulique dans un milieu tridimensionnel, base sur un couplage entre : une suite d'états d'équilibre élastique, l'écoulement du fluide à l'intérieur de la fracture, et le phénomène de rupture. Pour l'élasticité, une formulation intégro-variationnelle permet de se ramener à un problème bidimensionnel tout en évitant les singularités classiques. la résolution numérique se fait par éléments finis sur un maillage mobile. On résout d'abord le problème d'équilibre. les résultats numériques pour l'ouverture et les facteurs d'intensité de contrainte sont en très bon accord avec des solutions analytiques disponibles pour des géométries simples. puis le couplage et la propagation sont réalisés. les résultats sont compares a ceux obtenus par des modèles plus simples dans l'industrie pétrolière.