Procédé de prédiction de la propagation d’une fissure par fatigue dans une pièce mécanique soumise à un champ de contraintes externe (9011) et modélisée numériquement entre un instant initial t 0 et un instant final t f avec un pas de Δt, comprenant les étapes suivantes : une étape de maillage à l’instant t 0 ; une étape de détermination (802) d’un champ de contraintes internes (9032) à appliquer aux lèvres de la fissure, par calcul du champ de contraintes (9022) aux mêmes points du maillage pour la structure non fissurée (902) lorsque le champ de contraintes externes (9011) lui est appliqué, le champ de contrainte interne (9032) de la pièce réelle fissurée (903) étant alors égal et de signe opposé au champ de contrainte interne (9022) de la pièce non fissurée équivalente (902) ; caractérisé en ce que l’étape de maillage comprend les sous-étapes suivantes : Maillage (801) d’une pièce non fissurée (902) identique à la pièce fissurée (901) ; Maillage numérique (800) de la fissure ; Insertion du maillage de fissure (800) dans le maillage (801) de la pièce non fissurée (902), de sorte à obtenir un maillage numérique (804) de la pièce fissurée (901) ; Le procédé comprenant en outre une étape de réduction de modèle (803) de la structure réelle à un modèle modal simplifié (805), la réduction de modèle consistant à exprimer le champ de contrainte à chaque nœud du maillage fissuré (804) sous forme d’une combinaison linéaire de champs de contraintes modaux (σ k ) associés chacun respectivement à un mode structurel k de la pièce, pour un nombre de modes prédéterminé ; une étape de calcul de facteurs d’intensité des contraintes modaux (806) le long d’un front de fissure pour le modèle modal simplifié (805) ; une étape de reconstitution (807) de l’évolution temporelle des FIC à chaque nœud du front de la fissure à partir des FIC modaux obtenus à l’étape (806) de calcul ; une étape d’estimation (808) de la valeur d’une avancée de fissure à partir des évolutions temporelles des FIC obtenues à l’étape (807) de reconstitution et au moyen d’une loi de propagation de fissure, l’étape d’estimation (808) comprenant la détermination de la direction d’avancée de chaque nœud du front de fissure par calcul de l’angle qui maximise la vitesse de propagation de ce nœud ; une étape (809) d’obtention d’un maillage fissuré propagé reconstitué à partir du maillage (804) de la pièce fissurée et de l’estimation (808) de la valeur d’avancée de fissure ; les étapes b. à g. étant répétées pour chaque pas de la propagation modélisée, de l’instant t 0 + Δ t jusqu’à l’instant final t f , le procédé étant mis en œuvre par une unité de traitement de données (12) d’un système (10), l’unité de traitement comprenant des moyens de calcul (14). Procédé selon la revendication précédente, comprenant en outre : une étape (810) d’établissement de formulaires consistant à établir, pour chaque pas de propagation, des tableaux de correspondance entre des facteurs d’intensité de contraintes calculés et des dimensions de fissure ; une étape (811) de calcul d’une durée de vie restante de la pièce mécanique sur la base des formulaires établis. Procédé selon l’une quelconque des revendications précédentes, la réduction de modèle (803) étant une décomposition orthogonale propre (POD). Procédé selon l’une quelconque des revendications précédentes, la loi de propagation de fissure utilisée dans l’estimation (808) de la valeur d’avancée de fissure étant une loi de Paris. Procédé selon la revendication précédente, la loi de propagation de fissure utilisée dans l’estimation (808) de la valeur d’avancée de fissure étant une loi de Paris à laquelle une correction est appliquée. Procédé selon l’une quelconque des revendications précédentes, l’étape d’estimation (808) de la valeur d’avancée de fissure comprenant l’imposition, à chaque pas de propagation, d’une valeur de propagation maximale a max de la fissure de façon à réduire le nombre de pas de propagation calculés pour cette étape. Procédé selon l’une quelconque des revendications précédentes, le maillage de la structure fissurée étant décomposé en : une zone remaillée (32), à proximité de la fissure, pour laquelle les étapes d’obtention d’un maillage numérique (800) de la fissure et d’obtention d’un maillage numérique (804) de la pièce fissurée sont renouvelées à chaque pas ; une zone non-remaillée (31), en dehors de la zone remaillée (32), pour laquelle le maillage est déterminé à l’étape d’obtention d’un maillage numérique (801) de la pièce non fissurée et est maintenu pour l’ensemble des pas de propagation de la fissure ; et en ce que les étapes de calcul de facteurs d’intensité de contrainte (806) le long du front de fissure pour chaque mode structurel du modèle simplifié, d’obtention (809) d’un maillage fissuré propagé et le cas échéant de reconstitution (807) du facteur d’intensité des contraintes à tous les instants de la modélisation ne sont effectuées que pour les degrés de liberté de la pièce correspondant à la zone remaillée (32). Procédé de maintenance d’une pièce d’aéronef entre une première série (M n ) de missions en vol et une deuxième série (M n+1 ) de missions en vol, comprenant les étapes suivantes : Inspection (110) de la pièce afin de déterminer si une anomalie susceptible de causer l’apparition d’une fissure est apparue dans la pièce suite à la première série (M n ) de missions en vol Si oui, mise en œuvre du procédé selon l’une quelconque des revendications 1 à 7 afin de déterminer une durée de vie de la pièce ; Sinon, remise en service (112) de l’aéronef pour la deuxième série de missions en vol ; Procédé de maintenance selon la revendication précédente, dans lequel l’étape a.i. de détermination d’une durée de vie est suivie d’une comparaison avec une durée de vie de la pièce non fissurée et/ou de comparaison avec une durée de vie minimum nécessaire à une certification de la pièce. Procédé d’établissement d’une base de données pour une pièce mécanique, reliant des données d’entrée à des données de sortie ; les données d’entrée comprenant des tailles et/ou géométries de fissures, des positionnements de fissures dans la pièce et des contraintes externes à la pièce ; les données de sortie comprenant des facteurs d’intensité des contraintes et/ou des durées de vie de la pièce, comprenant les étapes suivantes : Mise en œuvre du procédé selon l’une quelconque des revendications 1 à 9 de sorte à obtenir des facteurs d’intensité des contraintes et/ou une durée de vie de la pièce mécanique pour les conditions d’entrée ; Enregistrement, dans la base de données, d’une correspondance entre les données d’entrée et les données de sortie. Système comprenant une unité de traitement (12) comprenant des moyens de calcul (14) et une mémoire (16), l’unité étant configurée pour mettre en œuvre un procédé selon l’une quelconque des revendications précédentes. Produit programme d’ordinateur configuré pour être mis en œuvre par un système selon la revendication précédente, et comprenant des instructions pour provoquer la mise en œuvre d’un procédé selon l’une quelconque des revendications 1 à 10.