FR 2970337 A1 20120713 FR 1150099 A 20110106 La présente invention est relative à un procédé de caractérisation ultrasonore conjointe de l'épaisseur et du taux de porosité de matériaux ayant des propriétés ultrasonores atténuantes. Cette invention relève du domaine des opérations de contrôle non destructif de pièces constituées de matériaux ayant des propriétés d'atténuation des ultrasons. Ces matériaux sont, par exemple, les matériaux composites ou les matériaux métalliques à gros grains utilisés dans l'aéronautique. Contexte de l'invention et problèmes posés La grande majorité des pièces de construction aéronautique en cours de production sont inspectées par des méthodes ultrasonores permettant de vérifier la qualité des matériaux et structures mis en jeu. Lors de ces opérations d'inspection on souhaite obtenir des informations notamment sur la présence d'hétérogénéités assimilables à des défauts, sur la présence de porosité dans le matériau, sur l'épaisseur géométrique des zones inspectées. On s'intéresse ici particulièrement à la mesure du taux de porosité et de l'épaisseur d'une pièce en matériau composite, ou plus généralement, ayant des propriétés d'atténuation des ultrasons. Le problème est alors de mesurer, avec une incertitude la plus réduite possible, à la fois l'épaisseur du matériau inspecté et son taux de porosité. Dans les procédés d'évaluation de porosité dans un matériau actuellement utilisés, la quantité physique communément mesurée est l'atténuation d'une onde ultrasonore traversant le matériau considéré. Cette atténuation est mesurée sur un écho de référence ayant traversé le matériau, par exemple sur l'écho de fond de pièce comme indiqué sur la fiqure 1. On peut citer notamment, dans cette famille de procédés, un procédé décrit dans la demande de brevet US 2007/0089479. Dans celui-ci, on réalise en premier lieu un échantillon dans un matériau de référence homogène (c'est-à-dire n'étant pas en matériau composite) présentant des inclusions maîtrisées représentatives des bulles présentes dans la matrice du matériau composite de la pièce. On effectue ensuite une mesure d'atténuation ultrasonore sur cet échantillon. Puis on compare le résultat avec le résultat de l'inspection de la pièce en matériau composite, dont on cherche à déterminer le taux de porosité. Un autre procédé de mesure de porosité, permettant d'améliorer la caractérisation de la porosité, en tenant compte notamment des caractéristiques fréquentielles des capteurs utilisés (et donc des longueurs d'ondes mises en jeu par rapport aux dimensions des porosités) a été proposé dans la demande de brevet WO2009130140 (Al) ("Procédé et dispositif de contrôle non-destructif par ultrasons de la porosité d'une pièce en matériau composite") par la demanderesse. Dans ce procédé, la mesure de l'épaisseur de la pièce analysée se fait par mesure du « temps de vol » qui correspond au temps mis par l'onde ultrasonore pour parcourir la distance séparant les deux faces d'un matériau inspecté. Une valeur d'épaisseur est ensuite obtenue en multipliant le temps mesuré par une vitesse de propagation dans le matériau (voir fiqure 2). Plus précisément, dans ce procédé, on établit un ensemble de résultats de référence correspondant à plusieurs échantillons de référence représentatifs de différents niveaux de porosité du matériau composite considéré. Ces résultats de référence sont obtenus par calcul à partir d'un ensemble restreint de données expérimentales et rendent possible l'évaluation d'un taux de porosité pour toute épaisseur de matériau et quel que soit le transducteur ultrasonore utilisé. Ce procédé et les autres procédés actuels souffrent cependant d'incertitudes de mesure élevées. Par ailleurs, ces valeurs sont en général respectées tant que le matériau inspecté est sain (c'est-à-dire exempt de porosité) mais potentiellement dépassées si le matériau devient poreux, puisque la porosité modifie la vitesse de propagation. De telles incertitudes de mesures se révèlent très pénalisantes pour la conception par exemple de pièces d'aéronefs, compte tenu de la prise en compte insuffisante des phénomènes de propagation ultrasonore à l'origine des variabilités observées. Il en résulte la nécessité d'augmenter les marges de sécurité de conception, ce qui conduit à un alourdissement des pièces pour assurer la tenue de contraintes du cahier des charges. Exposé de l'invention La présente invention a donc pour but de réduire cette incertitude de mesure de porosité et d'épaisseur d'un échantillon de matériau. A cet effet, elle vise selon un premier aspect un procédé d'estimation de porosité d'une pièce en matériau ayant des propriétés ultrasonores atténuantes, ledit procédé utilisant une caractérisation ultrasonore conjointe de l'épaisseur et du taux de porosité de ladite pièce. En effet, il a été noté par les inventeurs que la vitesse de propagation ultrasonore varie avec l'épaisseur du matériau ainsi qu'avec le taux de porosité de celui-ci. Dès lors, le fait de considérer une vitesse de propagation constante introduit des incertitudes importantes dans la mesure de l'épaisseur. Les fiqures 3a, 3b et 3c donnent des ordres de grandeur de variation de ces vitesses. Une caractérisation conjointe de l'épaisseur et de la porosité de la pièce permet de résoudre ce problème et de réduire significativement l'incertitude de mesure. Dans une mise en oeuvre préférée, le procédé de caractérisation conjointe d'épaisseur et de porosité, utilise au moins un transducteur ultrasonore, et comporte une étape initiale 100 d'estimation d'une première valeur du taux de porosité de la pièce. Cette estimation est, dans un cas particulier de mise en oeuvre, réalisée par une comparaison d'une valeur d'atténuation par ladite pièce d'une onde ultrasonore émise par chaque transducteur ultrasonore avec au moins une valeur d'atténuation de référence d'une onde ultrasonore de chaque dit transducteur ultrasonore par un échantillon de référence Rnm d'épaisseur em et représentatif d'un taux de porosité Pn connu du matériau composite dans lequel la pièce est fabriquée. Plus particulièrement, chaque valeur d'atténuation de référence est calculée entre : - une fonction d'onde représentative d'une onde ultrasonore émise par chaque transducteur ultrasonore, - et le résultat d'un calcul de filtrage de la fonction d'onde par une fonction de transfert Fnm représentative d'atténuations fréquentielles subies par des ondes ultrasonores traversant l'échantillon de référence Rnm Le procédé comporte une étape 200 itérative dans laquelle, partant de la première valeur estimée du taux de porosité, on réalise successivement des sous-étapes suivantes : 210 - estimation de la vitesse moyenne / du front d'onde généré par le transducteur utilisé sur l'épaisseur É de la pièce avec le taux de porosité P` précédemment estimés, 220 - nouvelle estimation de l'épaisseur É+' déduite de la mesure du temps de vol et de l'évaluation de la vitesse moyenne / calculées à l'étape 210, 230 - nouvelle estimation du taux de porosité P'+', en utilisant la nouvelle estimation d'épaisseur É+' calculée à l'étape 220, dans laquelle la nouvelle estimation du taux de porosité 0+9 est réalisée par une comparaison d'une valeur d'atténuation a par ladite pièce d'une onde ultrasonore émise par chaque transducteur ultrasonore avec au moins une valeur d'atténuation de référence a' d'une onde ultrasonore de chaque dit transducteur ultrasonore par un échantillon de référence R' et représentatif d'un taux de porosité P' connu du matériau composite dans lequel la pièce est fabriquée, 240 - calcul de la différence entre la nouvelle estimation de l'épaisseur et l'estimation précédente d'épaisseur, et de la différence entre la nouvelle estimation de taux de porosité et l'estimation précédente de taux de porosité, Ces sous-étapes 210, 220, 230 étant répétées jusqu'à l'obtention d'un critère de convergence prédéterminé. Avantageusement, le critère de convergence consiste à vérifier que a) la différence entre la nouvelle estimation de l'épaisseur É+' et l'estimation précédente d'épaisseur É, soit inférieure à un seuil prédéterminé, ou que b) la différence entre la nouvelle estimation de taux de porosité 0+9 et l'estimation précédente de taux de porosité P` soit inférieure à un second seuil prédéterminé. Présentation des figures Les caractéristiques et avantages de l'invention seront mieux appréciés grâce à la description qui suit, description qui expose les caractéristiques de l'invention au travers d'un exemple non limitatif d'application. La description s'appuie sur les figures annexées qui représentent : Figure 1 (déjà citée) : une illustration de mesure de l'atténuation ultrasonore pour l'évaluation de porosité : le matériau correspondant au signal rouge a un taux de porosité supérieur à celui du signal bleu Figure 2 (également déjà citée) : une illustration de mesure du temps de 15 vol pour l'évaluation de l'épaisseur. L'épaisseur est égale à : temps de vol aller retour x vitesse / 2 Figures 3a, 3b et 3c (déjà citées) : variation de vitesse des ultrasons en fonction (a) de l'épaisseur (matériau non poreux), (b) du taux de porosité (pour une épaisseur donnée), (c) épaisseur pour différents taux de porosité 20 (P3>P2>P1) Figure 4 : une représentation schématique des éléments du dispositif utilisé dans la mise en oeuvre du procédé, Figure 5a : des courbes caractéristiques d'atténuation (dB/mm) en fonction de la fréquence pour l'ensemble des taux de porosité de la base de 25 données de référence, Figure 5b : des courbes caractéristiques de vitesse de phase (m/s) en fonction de la fréquence pour l'ensemble des taux de porosité de la base de données de référence, Figure 5c : un spectre fréquentiel de capteur, 30 Figure 6 : un schéma de principe de l'algorithme itératif d'évaluation de l'épaisseur et de mise à jour du taux de porosité. Description détaillée d'un mode de réalisation de l'invention La présente invention concerne un procédé de caractérisation conjointe d'épaisseur et de porosité d'une pièce 2 en matériau composite comportant des fibres maintenues dans une matrice, par un contrôle non-destructif par ultrasons de ladite pièce 2. Ce procédé est destiné à être mis en oeuvre sous la forme d'un logiciel exécuté, par exemple, sur un micro-ordinateur 20 de type connu en soi, ou sur tout autre type de calculateur numérique (voir fiqure 4). Ledit micro-ordinateur comporte une interface adaptée à recevoir des informations de capteurs externes. Le micro-ordinateur 20 comporte par ailleurs des moyens 21 de mémorisation d'informations numériques, que celles-ci soient des résultats intermédiaires de calcul ou des informations numériques reçues de capteurs externes. Ces divers éléments sont connus en soi et ne sont donc pas détaillés plus avant ici. Le dispositif utilisé pour le contrôle non-destructif par ultrasons de la pièce en matériau composite 2 comporte, outre les moyens de calcul 20 et les moyens de mémorisation 21, des équipements pour l'inspection ultrasonore, notamment un transducteur ultrasonore 22 et un générateur/ récepteur 23 de signaux, généralement électriques. En ce qui concerne le principe utilisé dans le procédé d'estimation conjointe d'épaisseur et de porosité de matériau, l'invention WO2009130140 de la demanderesse, mentionnée ci-dessus, permet d'évaluer le taux de porosité du matériau en un point donné à partir de la valeur mesurée d'atténuation et en postulant une valeur de vitesse moyenne. On rappelle préliminairement que, dans cette demande WO2009130140, on définit un procédé non-destructif d'évaluation d'un taux de porosité P d'une pièce 2 en matériau composite, au moyen d'au moins un transducteur ultrasonore 22, le taux de porosité P de la pièce 2 étant estimé (dans une étape 65) par une comparaison d'une valeur d'atténuation a par ladite pièce 2 d'une onde ultrasonore émise par chaque transducteur ultrasonore 22 avec au moins une valeur d'atténuation a' de référence d'une onde ultrasonore de chaque dit transducteur ultrasonore 22 par un échantillon de référence R' représentatif d'un taux de porosité P' connu du matériau composite dans lequel la pièce 2 est fabriquée. Dans une étape préalable 64 de l'invention WO2009130140, chaque valeur d'atténuation a' de référence est calculée par comparaison : - d'une fonction d'onde Fo(f) représentative d'une onde ultrasonore émise par chaque transducteur ultrasonore 22, - avec le résultat d'un calcul de filtrage de la fonction d'onde Fo(f) par une fonction de transfert F'm représentative d'atténuations fréquentielles subies par des ondes ultrasonores traversant l'échantillon de référence Rnm La fonction d'onde Fo(f) est déterminée, dans une étape préalable 62 de l'invention WO2009130140, sans la pièce 2 et sans l'échantillon de référence Rnm en mesurant une onde ultrasonore émise par chaque transducteur ultrasonore 22 dans un milieu de couplage, dans des conditions voisines de celles devant être mises en oeuvre pour la pièce 2. Dans l'étape 65 de l'invention WO2009130140, pour améliorer la précision du taux de porosité estimé, ledit taux de porosité de la pièce 2 est estimé en comparant la valeur d'atténuation a par la pièce 2, avec des valeurs d'atténuation a' de référence calculées pour au moins N échantillons de référence Rnm, N étant supérieur ou égal à deux, de niveaux de porosité P' connus différents, 1 n N, chaque échantillon de référence étant caractérisé par une fonction de transfert F'm propre. Les fonctions de transfert F'm sont établies dans une étape 63 de l'invention WO2009130140, après détermination d'une épaisseur E de la pièce, et l'épaisseur em de chaque échantillon de référence R'm est proche de l'épaisseur E de la pièce 2. Dans l'étape 64 de l'invention WO2009130140, les valeurs d'atténuation de référence a' sont calculées, préalablement à la détermination de l'épaisseur de la pièce E, pour M échantillons de référence Rnm, par taux de porosité, d'épaisseurs différentes em (1 m M), M étant supérieur ou égal à deux. Dans l'étape 65 de l'invention WO2009130140, le taux de porosité P estimé est égal au taux de porosité pour lequel la valeur d'atténuation de référence a' la plus proche de la valeur d'atténuation a par la pièce 2 a été calculée, éventuellement la plus proche par valeurs supérieures. Dans l'étape 63 de l'invention WO2009130140, la fonction de transfert Fnm d'un échantillon de référence Rnm de taux de porosité P' et d'épaisseur em est déterminée dans le domaine des fréquences en calculant pour des valeurs de fréquences fdiscrètes l'expression : Fn (f)=exp/ em \ Sn(f , ~adn AS(fl expression dans laquelle Sn(f) est une mesure d'un spectre fréquentiel d'une onde ultrasonore ayant traversé un échantillon de référence réel R' (on note ici à titre de clarification que le terme Rnm désigne dans la présente description l'échantillon Rn lorsqu'il est considéré avec une épaisseur em) d'épaisseur dn, réalisé en matériau composite de même type que celui de la pièce 2 et de taux de porosité Pn, S(f) est une mesure d'un spectre fréquentiel d'une onde ultrasonore de mêmes caractéristiques n'ayant pas traversé ledit échantillon de référence réel Rn, a et A sont des facteurs positifs qui dépendent du mode de mesure desdits spectres fréquentiels S(f) et Sn(f). Le procédé, objet de la présente demande, utilise une étape 200 itérative dans laquelle, partant d'un taux de porosité P' et d'une épaisseur E' de la pièce 2 préalablement estimés (le taux de porosité P' étant par exemple mais non limitativement estimé par la méthode de l'invention WO2009130140) on réalise successivement des sous-étapes suivantes : 210 - évaluation de la vitesse moyenne / du front d'onde généré par le transducteur ultrasonore 22 utilisé sur l'épaisseur É considérée de matériau avec le taux de porosité P` précédemment évalué, 220 - nouvelle estimation de l'épaisseur É+' déduite de la mesure du temps de vol et de l'évaluation de la vitesse moyenne / calculées à l'étape 210, 230 - nouvelle estimation du taux de porosité P'+' en utilisant la nouvelle estimation de l'épaisseur É+' et la procédure d'estimation de l'étape 65 de la demande WO2009130140 (voir plus haut), 240 - calcul de la différence entre la nouvelle estimation de l'épaisseur É+' et l'estimation précédente d'épaisseur E', et de la différence entre la nouvelle estimation de taux de porosité 0+9 et l'estimation précédente de taux de porosité P`, Ces sous-étapes 210, 220, 230 étant répétées jusqu'à ce que a) la différence entre la nouvelle estimation de l'épaisseur É+' et l'estimation précédente d'épaisseur E', soit inférieure à un seuil prédéterminé, ou que b) la différence entre la nouvelle estimation de taux de porosité 0+9 et l'estimation précédente de taux de porosité P` soit inférieure à un second seuil prédéterminé. Ces deux valeurs convergent alors vers des valeurs limites de taux de porosité PL et d'épaisseur EL. Lorsque l'un ou l'autre des critères d'arrêt est vérifié dans l'étape 240, cela signifie que les valeurs d'atténuation a et de vitesse a sont physiquement cohérentes et que l'on a convergé vers la caractérisation d'un matériau « pouvant exister » (les atténuations et porosités sont stables d'une itération sur l'autre). En effet, l'atténuation aet la variation de la vitesse a de phase sont des grandeurs physiques duales. Ce sont deux manifestations ayant pour origine les mêmes phénomènes physiques, et pour lesquels la connaissance de l'un permet théoriquement de déduire l'autre. En particulier, la relation, connue en soi, de Kramers-Krdnig lie ces deux grandeurs (v pour vitesse moyenne du front d'onde, a pour atténuation et co=2*pi*f la fréquence) par la formule : 1 1 _ 2 a(co') dom, v(co) v(coo) n wo co On constate avec cette équation que la vitesse moyenne du front d'onde v, au même titre que l'atténuation a, est une fonction qui varie avec la fréquence co. Cette dépendance de la vitesse moyenne v du front d'onde généré par le transducteur ultrasonore 22, et de l'atténuation a à la fréquence co, est responsable de la variation de la vitesse v en fonction de l'épaisseur traversée. Ces variations de vitesse v proviennent en fait du caractère dispersif (c'est-à-dire d'une variation de la vitesse de phase v en fonction de la fréquence co) de la propagation ultrasonore dans les milieux atténuants. En pratique, le procédé de caractérisation conjointe de porosité P et d'épaisseur E d'une pièce 2 utilise, comme on l'a vu, un traitement logiciel de données d'acquisition ultrasonore. Le traitement comporte la connaissance préliminaire d'un ensemble de données : A/ Des données caractéristiques du matériau constituant la pièce 2: i. Des courbes caractéristiques de la variation de l'atténuation a' (dB/mm) en fonction de la fréquence co pour l'ensemble des taux de porosité P' de la base de données de référence (fiqure 5a). la méthode d'obtention de telles courbes est décrite dans la demande WO2009130140, déjà citée, ii. Des courbes caractéristiques de la variation de la vitesse de phase v' en fonction de la fréquence co pour l'ensemble des taux de porosité Pn de la base de données de référence (fiqure 5b), B/ Des données caractérisant le transducteur ultrasonore 22 et comportant son spectre fréquentiel (fiqure 5c). La méthode d'obtention d'un tel spectre fréquentiel (fonction d'onde Foe est décrite dans l'étape 62 de la demande WO2009130140, déjà citée. Ces données caractéristiques du matériau et du capteur 22 sont, dans le présent exemple non limitatif de mise en oeuvre du procédé, mémorisées dans les moyens de mémorisation 21 du micro-ordinateur 20 Les données d'acquisition, associées à la pièce 2 que l'on souhaite caractériser, sont constituées de : - l'amplitude d'un écho de référence de la pièce considérée, - ainsi que son temps d'arrivée (temps de vol). Pour obtenir ces données, la pièce en matériau composite 2 est, au moins partiellement, immergée dans un milieu de couplage 3. Le transducteur ultrasonore 22 est placé, pour l'évaluation de la porosité, de sorte à émettre des ondes ultrasonores qui ont une incidence sensiblement normale à une face avant 2a de la pièce en matériau composite, en regard dudit transducteur ultrasonore 22. Le transducteur ultrasonore 22 est utilisé en émission et en réception, et mesure des échos dus à des discontinuités d'impédances acoustiques. Un écho, dit « d'entrée », est reçu pour une interface entre la face avant 2a de la pièce 2 et le milieu de couplage 3, et un écho, dit « de fond », est en général reçu pour une interface entre ledit milieu de couplage et une face arrière 2b de ladite pièce, située du côté de la pièce opposé au transducteur ultrasonore 22. Un écho peut également être reçu pour une interface entre le milieu de couplage 3 et une plaque réfléchissante 4, par exemple en verre, qui est placée en vis-à-vis de la face arrière 2b de la pièce 2. L'inspection est alors dite par « double transmission » du fait que les ondes ultrasonores traversent deux fois la pièce 2. Il est également possible d'effectuer une inspection dite par « transmission », en utilisant un premier transducteur ultrasonore 22 pour l'émission du côté de la face avant 2a, et second un transducteur ultrasonore (non illustré sur les figures) pour la réception du côté de la face arrière 2b. Pour évaluer le taux de porosité P de la pièce 2, l'atténuation a par ladite pièce d'une onde ultrasonore émise par le transducteur ultrasonore 22 est mesurée. Par taux de porosité, il faut comprendre une valeur représentative de caractéristiques structurelles ou mécaniques de la pièce induites par la présence de bulles dans la matrice du matériau composite constituant la pièce 2. Le taux de porosité P est évalué sous la forme d'un pourcentage volumique de « vide » (c'est à dire de bulles) ou d'un niveau d'abattement de la résistance mécanique de la pièce. Le procédé de caractérisation conjointe de porosité et d'épaisseur, fournit en sortie : - L'épaisseur EL du matériau inspecté - Le taux de porosité PL du matériau inspecté, par application de l'algorithme décrit schématiquement sur la figure 6: 15 Un exemple de sorties de boucle de l'algorithme itératif, utilisé par le procédé de caractérisation conjointe, est illustré dans le tableau ci-dessous: Itération 0 Itération 1 Itération 2 Itération 3 Vitesse (m/s) 2950 2626 2613 2613 Epaisseur (mm) 11.80 10.50 10.45 10.45 Taux de Porosité (% 2.7 2.8 2.9 2.9 vol.) 20 Il est à noter que l'établissement des courbes caractéristiques de vitesse de phase v en fonction de la fréquence co, pour les taux de porosité P' faisant partie de la base de référence peut se faire soit : - par exploitation directe des mesures, par analyse spectroscopique en analysant la phase des spectres des signaux. 25 - par calcul de la vitesse de phase v en utilisant la relation de Kramers-Krdnig et une vitesse mesurée à basse fréquence. Ce procédé s'adresse notamment aux usines de production effectuant des contrôles de la qualité de pièces en matériaux composites ou métalliques à gros grains. Il permet de réduire significativement, par exemple autour de 20/0 (et 100/0 avec une méthode utilisant le temps de vol seul) les incertitudes de mesure sur la porosité des pièces, et donc autorise une meilleure connaissance prédictive de leurs valeurs de résistance mécanique. Variantes de réalisation La portée de la présente invention ne se limite pas aux détails des formes de réalisation ci-dessus considérées à titre d'exemple, mais s'étend au contraire aux modifications à la portée de l'homme de l'art. On a mentionné dans la description un matériau composite, mais la présente invention s'applique plus généralement à tout matériau ayant des propriétés ultrasonores atténuantes. REVENDICATIONS1. Procédé d'estimation de porosité d'une pièce en matériau ayant des propriétés ultrasonores atténuantes, caractérisé en ce que ledit procédé utilise une caractérisation ultrasonore conjointe de l'épaisseur et du taux de porosité de ladite pièce. 2. Procédé selon la revendication 1, utilisant au moins un transducteur ultrasonore, le procédé comportant une étape initiale 100 d'estimation d'une première valeur du taux de porosité et de l'épaisseur de la pièce, le procédé étant caractérisé en ce qu'il comporte une étape 200 itérative dans laquelle, partant des valeurs estimées précédemment du taux de porosité P' et de l'épaisseur É de la pièce (2), on réalise successivement des sous-étapes suivantes : 210 - estimation de la vitesse moyenne / du front d'onde généré par le transducteur (22) utilisé sur l'épaisseur É de la pièce (2) avec le taux de 15 porosité P`, précédemment estimé, 220 - nouvelle estimation de l'épaisseur É+' déduite de la mesure du temps de vol et de l'évaluation de la vitesse moyenne / calculées à l'étape 210, 230 - nouvelle estimation du taux de porosité P'+', en utilisant la 20 nouvelle estimation d'épaisseur É+' calculée à l'étape 220, dans laquelle la nouvelle estimation du taux de porosité 0+9 est réalisée par une comparaison d'une valeur d'atténuation a par ladite pièce (2) d'une onde ultrasonore émise par chaque transducteur ultrasonore (22) avec au moins une valeur d'atténuation de référence a' d'une onde ultrasonore de chaque dit 25 transducteur ultrasonore (22) par un échantillon de référence R' représentatif d'un taux de porosité P' connu du matériau composite dans lequel la pièce (2) est fabriquée, 240 - calcul de la différence entre la nouvelle estimation de l'épaisseur É+' et l'estimation précédente d'épaisseur É, et de la différence entre la 30 nouvelle estimation de taux de porosité 0+9 et l'estimation précédente P` de taux de porosité,ces sous-étapes 210, 220, 230 étant répétées jusqu'à l'obtention d'un critère de convergence prédéterminé. 3. Procédé selon la revendication 2, caractérisé en ce que l'estimation d'une première valeur du taux de porosité est réalisée par une comparaison d'une valeur d'atténuation par ladite pièce d'une onde ultrasonore émise par chaque transducteur ultrasonore avec au moins une valeur d'atténuation de référence an d'une onde ultrasonore de chaque dit transducteur ultrasonore par un échantillon de référence Rnm d'épaisseur em et représentatif d'un taux de porosité P' connu du matériau composite dans lequel la pièce 2 est fabriquée. 4. Procédé selon la revendication 3, caractérisé en ce que chaque valeur d'atténuation de référence an est calculée par comparaison : - d'une fonction d'onde Fo(f) représentative d'une onde ultrasonore 15 émise par chaque transducteur ultrasonore, - avec le résultat d'un calcul de filtrage de la fonction d'onde Fo(f) par une fonction de transfert Fnm représentative d'atténuations fréquentielles subies par des ondes ultrasonores traversant l'échantillon de référence Rnm, 20 5. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que le critère de convergence consiste à vérifier que a) la différence entre la nouvelle estimation de l'épaisseur É+' et l'estimation précédente d'épaisseur É, soit inférieure à un seuil prédéterminé, ou que b) la 25 différence entre la nouvelle estimation de taux de porosité 0+9 et l'estimation précédente de taux de porosité P` soit inférieure à un second seuil prédéterminé.