FR 2832804 A1 20030530 FR 0115172 A 20011123 PROCEDE DE MESURE ACCELEROMETRIQUE, COMPENSEE THERMIQUEMENT, AVEC AU MOINS UN DISPOSITIF A DEUX RESONATEURS VIBRANTS La présente invention concerne le domaine des mesures accélérométriques, et plus spécifiquement elle concerne des perfectionnements apportés dans le domaine des mesures accélérométriques effectuées avec des dispositifs à deux résonateurs vibrants travaillant respectivement dans des états de contraintes de signes opposés (en général en traction et en compression). Ainsi, il est déjà connu d'utiliser des dispositifs à deux résonateurs vibrants comportant deux poutres vibrantes respectives, parallèles et identiques, travaillant l'une en traction et l'autre en compression, pour détecter une composante d'accélération (voir par exemple les documents FR 00 10675 et FR-A-2 784 752 au nom de la Déposante ; FR-A-2 685 964, FR-A-2 784 752). Comme illustré à titre d'exemple à la figure 1 des dessins ciannexés, un tel dispositif accélérométrique 1 comprend, de façon très schématique, deux accéléromètres miniatures Accu, Acc2 à masses pendulaires non asservies. Chacun d'eux comporte une masse pendulaire (ou masse sismique) M sollicitée vers une position de repos par une liaison avec un socle S qui est constituée d'une poutre ou lame vibrante P (les moyens d'excitation en vibration et de détection de la vibration n'étant pas illustrés) et d'une pièce C formant charnière. Les deux poutres Pi, P2 sont identiques et parallèles et sont disposées en position inversée l'une par rapport à l'autre, de telle sorte que, sous l'action des deux masses sismiques identiques respectives Mi, M2, l'une Pi des poutres travaille en compression et que l'autre P2 des poutres travaille en traction. L'ensemble des deux accéléromètres individuels Accu, Acc2 est enfermé dans un boîtier étanche B dans lequel généralement un vide poussé est établi. Lorsque le dispositif accélérométrique 1 est soumis à une accélération y dirigée suivant l'axe sensible, c'est-à-dire par exemple perpendiculaire à l'étendue des poutres vibrantes Pi, P2, il en résulte une variation des fréquences de résonance des poutres en relation avec l'accélération appliquée et avec la contrainte (compression ou traction) subie par les poutres : fi = fol + kll kl2 Y2+ kl3 Y3 f2 = f02 - k21 &gamma;- k22 &gamma;2 - k23 &gamma;3 où foi et f02 sont les fréquences de résonance des deux poutres respectives Pi et P2 en l'absence d'accélération, et les coefficients k sont des constantes. Ce type de dispositif accélérométrique présente plusieurs intérêts, notamment en raison de son faible encombrement (structure plate dont l'épaisseur peut être inférieure à 1 mm pour des dimensions en plan de l'ordre du centimètre) et de sa fabrication simple et économique par des procédés d'attaque chimique sur une tranche de cristal en général piézoélectrique, mais pouvant être du silicium, avec possibilité de fabrication collective simultanée de plusieurs accéléromètres sur une même tranche de matériau. Toutefois, dans le type de dispositif accélérométrique considéré à double résonateur vibrant, les fréquences de résonance des poutres vibrantes Pi, P2 dépendent non seulement de l'accélération y à laquelle est soumis le dispositif, mais aussi de la température T des poutres. On est donc conduit à compléter les équations précédentes par un polynôme en T : Un calcul simple de la valeur de l'accélération y en fonction de la fréquence mesurée fi et/ou f2 nécessite de connaître la température des poutres Pi, P2. Or il n'est pas possible, en pratique, de disposer un capteur approprié, dans le boîtier B, à proximité immédiate des poutres vibrantes, ni de disposer un capteur thermique au contact même des poutres vibrantes Pi, P2 (ce qui constituerait la solution idéale pour obtenir une valeur exacte de la température) en raison de la perturbation que la présence de ce capteur provoquerait sur la vibration des poutres. Enfin, le recours à un procédé de mesure thermoélectrique serait susceptible de perturber la température de la poutre elle-même. L'invention a donc pour but de proposer un procédé de calcul de la valeur de l'accélération, compensée en température, à partir des mesures des fréquences fi, f2 des deux poutres vibrantes Pi, P2, sans avoir recours à un moyen quelconque de détection thermique disposé au contact des poutres vibrantes. A ces fins, l'invention propose un procédé de mesure accélérométrique par détermination de la valeur, compensée en fonction de la température, d'une accélération à partir de mesures de fréquences faites sur un dispositif accélérométrique à deux résonateurs vibrants de géométries similaires, travaillant respectivement dans des états de contraintes de signes opposés, lequel procédé se caractérise, selon l'invention, en ce qu'il comprend les étapes qui suivent : - le dispositif étant soumis à une composante d'accélération, on détermine les fréquences fi et f2 des vibrations respectives des deux résonateurs vibrants, sachant que les valeurs respectives de ces deux fréquences sont influencées par l'accélération et par la température ; on détermine la somme fi + f2 des deux fréquences mesurées, sachant que cette somme est influencée par la température ; à partir d'une corrélation préétablie de la variation de la quantité fi + f2 en fonction de la température, on relève les valeurs possibles Ti de température correspondant à la valeur de la somme fi + f2 ci-dessus déterminée ; on mesure la valeur (Gmes) d'une grandeur représentative de la température à proximité du dispositif à deux résonateurs vibrants ; à l'aide d'une fonction de transfert H préalablement établie, on détermine une valeur estimée Gestimée = H (Gmes) de la susdite grandeur à partir de ladite valeur mesurée Gmes, cette grandeur estimée étant représentative de la température estimée des résonateurs vibrants ; à partir de la valeur estimée Gestimée de ladite grandeur et des valeurs possibles de température Ti, on détermine, pour chaque valeur Ti, la probabilité P (tel) pour que la valeur de la température des résonateurs soit Ti ; on calcule une valeur pondérée T de la température des résonateurs vibrants T = L P (Ti) x T enfin, à partir d'une et/ou des fréquences mesurées fi et f2 et de la valeur pondérée T de la température, on détermine la valeur compensée en température de l'accélération subie par le dispositif accélérométrique. Grâce aux dispositions de l'invention, on est en mesure de tenir compte, dans le calcul de la valeur de l'accélération y, d'une valeur de température T, ou d'une grandeur représentative de cette température (par exemple une grandeur variant linéairement en fonction de la température) au voisinage des résonateurs-et donc intrinsèquement inexacte puisque n'étant pas la valeur de la température des résonateurs-et à partir d'un choix, basé par un calcul de probabilité, entre plusieurs valeurs possibles de la température. En pratique, les deux résonateurs vibrants ne peuvent pas être rigoureusement identiques du point de vue géométrique : il est alors souhaitable de tenir compte de cette différence géométrique en corrigeant la somme fi + f2 avec un terme d'erreur estimé à partir de l'accélération calculée lors du cycle de calcul précédent, s'il existe. En pratique toujours, la mise en oeuvre du procédé de l'invention avec un polynôme n (Tn) d'ordre n élevé se révèle complexe, en raison du nombre n élevé des solutions mathématiques, c'est-à-dire du nombre n de valeurs possibles de température Ti correspondant à une valeur donnée de la fréquence fi + fs. Il en résulte également un calcul plus long des probabilités P (Ti) et de la valeur T = z P (Ti) x Ti. Une valeur T d'une précision suffisante est obtenue, de façon concrète, avec un polynôme d'ordre 2 : n (T2) ou d'ordre 3 : Il (T3). La relation mathématique entre la fréquence, l'accélération et la température peut s'écrire dans le cas d'un polynôme en T d'ordre 2 et dans le cas d'un polynôme en T d'ordre 3. En pratique, les termes en y2 et y3 sont faibles et peuvent être négligés devant le terme en y. Les deux équations ci-dessus peuvent alors s'écrire sous forme simplifiée respectivement : La courbe correspondant à la fonction f (T) est une parabole ou une cubique, respectivement, dont l'extremum local (situé dans la fourchette des températures de mise en oeuvre) correspond à la température d'inversion ("turn over") To. Ainsi, à toute valeur f de la fréquence sont associable deux valeurs possibles de température Ta et Tb, dont une seule est la température exacte du résonateur vibrant. Dans le cas d'un polynôme du troisième ordre, la troisième solution est aberrante par rapport au domaine de température considéré et peut être rejetée d'entrée. La détermination de la température des résonateurs fait appel à une estimation probabiliste et consiste à associer, à la somme fi + f2, deux valeurs Ta et Tb, à calculer la probabilité P (Ta) pour que la température des résonateurs soit Ta et la probabilité P (Tb) pour que la température des résonateurs soit Tb, puis à calculer une valeur pondérée de la température des résonateurs Dans le cas de mesures accélérométriques et/ou gyroscopiques selon plusieurs axes (en pratique selon deux ou trois axes d'un référentiel) à l'aide de plusieurs dispositifs vibrants accélérométriques et/ou gyroscopiques et/ou gyrométriques respectifs (bloc senseur inertiel), on peut avantageusement prévoir que, selon chaque axe, on détermine une valeur de l'accélération compensée en température selon l'une ou l'autre des modalités exposées plus haut et, pour la détermination de l'une au moins des valeurs de l'accélération compensée en température, on utilise en tant que valeur pondérée de température, la valeur pondérée de température déterminée pour le dispositif vibrant associé à cet axe ainsi qu'au moins une valeur pondérée de température déterminée pour au moins un dispositif vibrant associé à au moins un autre axe. Une telle mise en oeuvre trouve un intérêt tout particulier dans le cas où tous les résonateurs vibrants des dispositifs accélérométriques et/ou gyroscopiques et/ou gyrométriques sont usinés dans un même bloc de cristal. L'invention sera mieux comprise à la lecture de la description détaillée qui suit de certains modes de réalisation donnés à titre d'exemples non limitatifs. Dans cette description, on se réfère aux dessins annexés sur lesquels : - la figure 2 est une représentation très schématique d'un dispositif accélérométrique à deux résonateurs vibrants agencé conformément à l'invention ; - les figures 3A et 3B sont des schémas par blocs illustrant respectivement deux mises en oeuvre préférées du procédé conforme à l'invention ; - la figure 4 est un schéma par blocs illustrant un perfectionnement du procédé de la figure 3A ; - la figure 5 est un schéma bloc illustrant de façon simplifiée un exemple d'agencement de capteur inertiel mettant en oeuvre les dispositions de l'invention ; - la figure 6 est un schéma par blocs, calqué sur celui de la figure 3A, détaillant une partie fonctionnelle du capteur inertiel de la figure 5 ; - la figure 7 est un schéma par blocs illustrant une variante de réalisation d'une portion de la partie fonctionnelle de la figure 6 ; et - la figure 8 est un schéma par blocs, calqué sur celui de la figure 3A, illustrant encore une autre variante de réalisation. L'invention vise à déterminer une valeur, corrigée en fonction de la température, de l'accélération y subie par un dispositif accélérométrique à deux résonateurs vibrants de géométries similaires, travaillant respectivement dans des états de contraintes de signes opposés. Dans l'exemple illustré, les deux résonateurs comprennent respectivement deux poutres vibrantes parallèles et sensiblement identiques travaillant respectivement notamment en traction et en compression, mais pouvant éventuellement être combinés à au moins un autre état de contrainte telle que la flexion et/ou la torsion. Dans l'exemple illustré, l'accélération y ou la composante de l'accélération qui est mesurée est dirigée sensiblement perpendiculairement aux poutres. Les fréquences de vibration des deux poutres peuvent s'écrire, en fonction de la température T et de l'accélération y : Selon l'invention, on détermine les fréquences fi et f2 de vibration des deux poutres vibrantes, puis on calcule la somme fi + f2. A partir d'une corrélation préétablie et mémorisée entre les valeurs de fréquence et les valeurs correspondantes de température f-7 n (Tn), on détermine, pour la fréquence fi + f2, les n valeurs Ti possibles de température. On procède à une évaluation de la température estimée des résonateurs. A cette fin, on mesure la valeur Gmes d'une grandeur représentative de la température à proximité du dispositif à deux résonateurs vibrants ; puis, à l'aide d'une fonction du transfert H préalablement établie, on détermine une valeur estimée Gestimée = H (Gmes) de la susdite grandeur Gmes, cette grandeur estimée étant représentative de la température estimée des résonateurs vibrants. A partir de la valeur estimée Gestimée de ladite grandeur représentative de la température estimée des résonateurs et des valeurs possibles de températures Ti, on détermine les probabilités P (Ti) pour que chaque valeur possible de température Ti soit la valeur exacte de la température des résonateurs. Une façon simple de mettre en oeuvre les dispositions précédentes consiste à effectuer directement une mesure de température, auquel cas la grandeur G représentative de la température est la température ellemême. C'est dans ce contexte qu'est exposé ce qui suit. A proximité immédiate du boîtier du dispositif accélérométrique et soit à l'intérieur, soit plus aisément à l'extérieur de celui-ci, on dispose un capteur de température qui délivre une valeur de température mesurée Tmes. En faisant intervenir une fonction de transfert H prédéterminée par des essais antérieurs, on peut évaluer une température estimée des résonateurs vibrants Testimée = H (Tmes). A partir des valeurs Testimée et Ti, on détermine, pour chaque valeur Ti, la probabilité P (tri) pour que la température des résonateurs soit Ti. Puis on détermine une valeur pondérée T de la température que l'on prend pour valeur de température des deux résonateurs. Ensuite de quoi, connaissant la température T, on est en mesure de calculer la valeur de l'accélération y à partir d'une et/ou des fréquences mesurées fi et f2. La mise en oeuvre d'un polynôme n (Tn) d'ordre n élevé entraîne des calculs complexes et longs, et rend difficile l'établissement de la corrélation f- TI (Tn) en raison du nombre élevé des valeurs possibles Ti associées à chaque valeur de fréquence. Une précision acceptable en pratique peut être obtenue avec un polynôme d'ordre 2 : Il (T2), ou de préférence d'ordre 3 : ri (T3). La fréquence de vibration s'exprime, en fonction de l'accélération et la température, par la relation dans le cas d'un polynôme en T d'ordre 2 ou par la relation dans le cas d'un polynôme en T d'ordre 3. Les termes en y2 et y3 étant faibles en pratique peuvent être négligés devant le terme en y, de sorte que les deux relations précitées peuvent s'écrire sous forme simplifiée, respectivement : Comme représenté schématiquement à la figure 3, on mesure en 10 et 11 respectivement les fréquences fi et f2 des deux résonateurs vibrants Pi et P2 soumises à une accélération y, sachant que les fréquences fi et f2 sont fonction de y et de la température T, cas d'un polynôme d'ordre 2 en T : (où Toi et sont les températures d'inversion ("turn over") respectives), ou bien, . dans le cas d'un polynôme d'ordre 3 en T : On constate que la somme fi + f2 est indépendante de 11 accélération y si ki = k2 et ne dépend que de la température. De plus, la différence fi-fi dépend essentiellement de l'accélération y si les termes de courbure ai et a2, ainsi que les valeurs d'inversion Toi et To2 sont proches les unes des autres, respectivement (les termes en T2 et en T3 s'éliminant). C'est sur l'exploitation des deux valeurs fi + f2 et fi-f qu'est fondé le procédé conforme à l'invention qui est exposé ciaprès en regard de la figure 3A (cas d'un polynôme d'ordre 2 en T) et de la figure 3B (cas d'un polynôme d'ordre 3 en T). A partir des mesures de la fréquence fi en 10 et de la fréquence f2 en 11, on détermine en 12 la somme fi + f2. A partir d'une corrélation préétablie, et tenue en mémoire en 14, entre la fréquence et la température des résonateurs (polynôme en T3 et/ou T2 et T), on relève en 13 les deux valeurs possibles de température Ta, Tb correspondant à la valeur fi + f2 issue des mesures fi et f2. Dans le cas d'un polynôme d'ordre 2 en T (cas illustré à la figure 3A), les deux valeurs recherchées Ta et Tb sont les deux solutions de l'équation du second degré en T. Dans le cas d'un polynôme d'ordre 3 en T (cas illustré à la figure 3B), l'équation du troisième degré en T a trois solutions, dont une cependant est aberrante par rapport au domaine de température considéré et doit être rejetée ; on est ainsi amené à tenir en mémoire en 25 (figure 3B) des valeurs limites, par exemple Tmin ; Tmax, définissant le domaine de température et à comparer en 26 les trois solutions de l'équation du troisième degré avec lesdites valeurs limites pour éliminer la solution aberrante et ne conserver que les deux solutions convenables Ta, Tb. Les deux valeurs de température Ta, Tb sont les deux valeurs possibles pouvant correspondre à la fréquence fi + f2, mais une seule des deux valeurs de température est exacte en pratique. Pour effectuer une approche aussi précise que possible de la température exacte des résonateurs (on suppose les deux résonateurs à la même température), on dispose un capteur thermique CTH à proximité immédiate des résonateurs. Ce capteur thermique peut être disposé à l'intérieur du boîtier B (comme montré en tirets à la figure 2), mais il peut également être disposé (car il s'agit d'une mise en oeuvre plus facile) à l'extérieur de ce boîtier comme illustré en trait plein à la figure 2. Ce capteur thermique CTH fournit en 15 (figures 3A et 3B) une valeur mesurée de température Tmes, laquelle valeur Tmes n'est qu'une valeur approchée, et non exacte, de la température des résonateurs Pi, P2. Préalablement à la mesure, on a établi expérimentalement une fonction de transfert H entre les valeurs exactes de température des résonateurs Pi, P2 situés à l'intérieur du boîtier B et les valeurs de température mesurées, simultanément, par le capteur CTH à l'extérieur du boîtier B. La fonction de transfert H est tenue en mémoire en 16. Pour fixer les idées, la fonction de transfert H peut être une fonction simple de deuxième ou troisième ordre, telle que, avec H (p) = T (p)/TCHT (p) H (p) = 1/ (1 + a p) 2 ou H (p) = 11 [ (1 + al p) (1 + a2 p) (1 + a3 p)] dans lesquelles p est l'opérateur de Laplace et a, (Xi, (Xz, a3 sont des constantes dépendant de la structure du capteur de température et du temps de propagation de la chaleur entre les poutres vibrantes situées à l'intérieur du boîtier B et le capteur de température CTH situé à l'extérieur de ce boîtier. Ainsi, à la température Tmes mesurée à l'extérieur du boîtier B, on est en mesure d'associer en 17, par la fonction de transfert H, une valeur estimée Test de la température des poutres Pi, P2 à l'intérieur du boîtier B. Ensuite, à partir des valeurs Ta, Tb et Test, on effectue en 18 des calculs de probabilité donnant la probabilité P (Ta) pour que la température Ta soit la température des poutres Pi, P2 et la probabilité P (Tb) pour que la température Tb soit la température des poutres Pi, P2. Pour fixer les idées, les probabilités P (Ta) et P (Tb) peuvent être calculées de la façon suivante : avec On pourrait également effectuer les calculs avec les fonctions suivantes : dans lesquels a est une constante. Les formules indiquées ici présentent l'avantage d'être simples et de conduire à des calculs rapides. Toutefois, si nécessaire, l'homme du métier est en mesure d'avoir recours à des formules plus complexes et plus précises, mais aussi entraînant des calculs plus longs. On calcule ensuite, en 19, une valeur T de la température qui est une pondération des valeurs possibles Ta et Tb affectées des coefficients de pondération précédemment calculés P (Ta) et P (Tb) : On définit ainsi une valeur T de température des résonateurs qui est statistiquement la plus vraisemblable. Enfin, à partir d'une des fréquences mesurées fi ou f2, ou à partir des deux fréquences mesurées fi et f2, et de la valeur pondérée T de la température qui vient d'être déterminée en 19, on détermine la valeur compensée en température de l'accélération y subie par le dispositif accélérométrique. De façon pratique, on calcule en 20 la quantité fi-fs. On sait que, dans le cas du polynôme d'ordre 2 en T, l'expression mathématique en est : Autrement dit, cette expression est du type : A partir de la valeur T de la température déterminée en 19, il est alors aisé de calculer, en 21, la valeur y à partir de T et de fi-f2. Dans le cas d'un polynôme d'ordre 3 en T, l'expression mathématique de fi-fez est plus complexe, mais est toutefois déterminable et permet également le calcul de la valeur de y à partir du calcul de la valeur fl-f2. On obtient finalement, en 22, une valeur de l'accélération y compensée en température. La corrélation entre la fréquence et la tempé- 3 2 rature (T3 et/ou T2, T), tenue en mémoire en 14, peut être mémorisée sous forme d'une équation mathématique du second ou du troisième degré en T dont les coefficients sont préétablis par des mesures expérimentales préalables effectuées sur le dispositif accélérométrique 1. Les deux valeurs possibles Ta et Tb sont alors calculées en résolvant cette équation pour chaque valeur donnée de fi + f2. On peut aussi tenir en mémoire en 14 un tableau donnant, pour toutes les valeurs possibles fi + f2, les deux valeurs possibles correspondantes Ta, Tb. De la même manière, la fonction de transfert H tenue en mémoire en 16 peut être mémorisée soit sous forme d'une relation algébrique Test = H (Tmes), soit sous forme d'un tableau de données préétablies donnant, pour chaque valeur Tmes, la valeur associée Test. L'expérience a montré que la courbe traduisant mathématiquement la relation entre la fréquence et la température (T3 et/ou T2, T) était susceptible de variation dans le temps, cette variation se traduisant essentiellement par une dérive de la fréquence pour une 3 2 température donnée, la forme de la courbe f (T3 et/ou T21 T) restant sensiblement inchangée. Pour obtenir un procédé de calcul de l'accélération y qui reste précis dans le temps, il est très souhaitable de prévoir une recalibration comme illustrée sur la figure 4 des dessins annexés. Dans ce but, on ajoute, aux étapes précédemment décrites, les étapes suivantes apparaissant à la figure 4 (laquelle est dérivée plus précisément de la figure 3A). A partir de la valeur mesurée de la température Tmes, on calcule en 22 la dérivée par rapport au temps dTmes/dt. Si on détecte, en 23, que cette dérivée est nulle ou voisine de zéro, on commande en 24 un processus de recalibration RECAL propre à agir sur la fonction de transfert tenue en mémoire en 16, sur le calcul/sélection en 13 des valeurs possibles de température Ta, Tb, et sur le calcul en 21 de la valeur fi-fez = p + q y + r T conduisant, d'une façon simple, à la détermination de l'accélération y. Enfin, on notera que l'accélération y est calculée de façon itérative, par passes de calcul successives, et que chaque passe de calcul prend en compte la valeur de l'accélération trouvée à l'issue de la passe de calcul précédente, en particulier pour la compensation de la somme fi + f2 évoquée plus haut lorsque ki + k2. En pratique, plusieurs dispositifs vibrants peuvent être regroupés pour fournir des informations accélérométriques et/ou gyrométriques ou gyroscopiques selon deux ou trois axes de référence (par exemple bloc senseur inertiel), tous ces dispositifs vibrants faisant individuellement appel à une mesure de température pour la détermination des valeurs compensées des fréquences de vibration. Selon l'invention, la mesure accélérométrique et/ou gyroscopique et/ou gyrométrique compensée en température, effectuée pour chaque dispositif à deux résonateurs vibrants, prend en compte non seulement sa propre valeur de température pondérée, mais aussi au moins certaines températures pondérées prédéterminées pour d'autres dispositifs vibrants de manière à obtenir une valeur de fréquence interpolée en relation avec au moins certaines des températures mesurées. A titre d'exemple, on a représenté très schématiquement, à la figure 5, un capteur inertiel 28 regroupant trois dispositifs vibrants de mesures accélérométriques selon trois axes x, y, z d'un référentiel, chaque dispositif vibrant pouvant être du type illustré à la figure 2 et étant référencé respectivement lc, 1y1 et 1z1, ainsi que trois dispositifs vibrants de mesures gyroscopiques par rapport aux trois axes x, y, z du référentiel, chaque dispositif vibrant pouvant être du type illustré à la figure 2 et étant référencé respectivement 1, 1 et 1. x y Z. A chaque dispositif vibrant précité est associée une unité de traitement d'information qui reçoit les informations fi, f2 et Tmes fournies par chaque dispositif. Comme illustré à la figure 6, chaque unité 27 est agencée (à titre d'exemple) conformément à la représentation donnée à la figure 3A. La figure 6 montre, à titre d'exemple, l'unité 27x1 associée au dispositif vibrant 1, étant entendu que les autres unités associées respectivement aux autres dispositifs vibrants lui sont identiques. Toutes les valeurs pondérées de température déterminées par toutes les unités 27, T, Ty1, Tz1, Tx2, Ty2, Tz2, sont envoyées à un dispositif de corrélation 29 qui, sur la base d'une corrélation préétablie (par exemple calcul de la moyenne), délivre une valeur corrélée de température Tc, qui est ensuite délivrée à chaque unité 27. Dans l'unité 27, la valeur corrélée Tc peut alors être utilisée directement en 21 pour calculer, à partir de fi-fs, la valeur de la composante d'accélération (y à la figure 6). On pourrait aussi envisager, comme illustré à la figure 7 qui ne montre qu'une partie de l'unité 27L que la valeur corrélée de température Tc soit combinée, selon une loi préétablie, en 30 avec la valeur pondérée'N, la valeur résultante étant utilisée pour le calcul de la composante y à partir de la valeur fi-fs. Bien entendu, de nombreuses variantes de corrélation en température peuvent être envisagées, par exemple : corrélation pour les températures des seuls dispositifs vibrants accélérométriques et corrélation pour les températures des seuls dispositifs vibrants gyroscopiques ; corrélation croisée entre les températures des dispositifs accélérométriques et gyroscopiques ; etc.... Un avantage économique peut être obtenu lorsque tous les dispositifs vibrants sont usinés dans un seul et unique bloc de cristal ou"wafer". Les dispositifs vibrants sont alors situés à proximité immédiate les uns des autres et il n'est plus besoin d'affecter un capteur de température à chacun d'eux : selon la configuration d'implantation mutuelle des dispositifs, un seul capteur de température disposé centralement peut suffire, ou à tout le moins deux capteurs de température peuvent être prévus, l'un associé aux dispositifs vibrants de mesures accélérométriques et l'autre aux dispositifs vibrants de mesures gyroscopiques ou gyrométriques. On recueille ainsi un gain de dimension pour le substrat et une économie quant à la capacité de traitement des informations. C'est en particulier dans le cas de la mise en oeuvre de plusieurs dispositifs vibrants localisés de façon adjacente que l'on peut être amené à traiter une grandeur G représentative de la température sans faire appel à la température elle-même. Un tel cas peut se présenter par exemple lorsqu'au dispositif 1 à deux résonateurs vibrants de fréquence fi et f2 est adjoint, de façon adjacente, un dispositif vibrant tel qu'un résonateur gyroscopique du type dit"en bol", dont la fréquence de résonance varie linéairement en fonction de la température et peut s'écrire f = mt + n où m et n sont des constantes. On comprend alors que mesurer la fréquence f de ce dispositif vibrant fournit une indication de température qui, à une constante près, est équivalente à une mesure directe de la température à l'aide d'un capteur de température. De ce fait, le traitement de la grandeur précitée G (T), si celle-ci varie de façon simple en fonction de la température, permet d'assurer le déroulement du procédé de mesure accélérométrique selon l'invention sans qu'il soit nécessaire de mettre en oeuvre un capteur de température à proprement parler. C'est un agencement de ce type qui est illustré à la figure 8 dont le schéma est calqué sur celui de la figure 3A. En 15, on mesure la valeur Gmes d'une grandeur G (T), notamment dans un dispositif vibrant adjacent au dispositif 1 à deux résonateurs vibrants de fréquence fi et f2, ladite grandeur G (T) étant ainsi représentative de la température à proximité dudit dispositif vibrant 1. Puis, à l'aide d'une fonction de transfert H préalablement établie en 16, on détermine une valeur estimée Gestimée = H (Gmes) de la susdite grandeur G (T), la valeur estimée Gestimée étant représentative de la température estimée des résonateurs vibrants de fréquence fi et f2. C'est ensuite cette valeur Gestimée qui est utilisée, en conjonction avec les valeurs Ta, Tb, pour déterminer les probabilités P (Ta) et P (Tb) évoquées plus haut, sans que la valeur de la température estimée Testimée intervienne de façon effective. Une telle solution présente l'intérêt d'éviter de mettre en oeuvre un ou plusieurs capteurs de température, ce qui est avantageux sur le plan économique. En outre, l'élimination d'un ou plusieurs composants dans un bloc senseur inertiel est toujours particulièrement apprécié en raison du gain de place et d'encombrement ainsi obtenu qui va dans le sens de la recherche de blocs senseurs inertiels toujours plus compacts. REVENDICATIONS 1. Procédé de mesure accélérométrique par détermination de la valeur, compensée en fonction de la température, d'une accélération à partir de mesures de fréquences faites sur au moins un dispositif accélérométrique à deux résonateurs vibrants de géométries similaires travaillant respectivement dans des états de contraintes de signes opposés, caractérisé en ce qu'il comprend les étapes qui suivent : - le dispositif étant soumis à une composante d'accélération (y), on détermine les fréquences fi (en 10) et f2 (en 11) des vibrations respectives des deux résonateurs vibrants, sachant que les valeurs respectives de ces deux fréquences sont influencées par l'accélération (y) et par la température (T) ; on détermine (en 12) la somme fi + f2 des deux fréquences mesurées, sachant que cette somme est influencée par la température ; - à partir d'une corrélation préétablie de la variation de la quantité fi + f2 en fonction de la température (en 14), on relève (en 13) les valeurs possibles Ti de température correspondant à la valeur de la somme fi + f2 ci-dessus déterminée ; on mesure la valeur (Gmes) d'une grandeur représentative de la température à proximité du dispositif à deux résonateurs vibrants ; - à l'aide d'une fonction de transfert H préalablement établie, on détermine une valeur estimée Gestimée = H (Gmes) de la susdite grandeur à partir de ladite valeur mesurée (Gmes), cette grandeur estimée étant représen- tative de la température estimée des résonateurs vibrants ; à partir de la valeur estimée (Gestimée) de ladite grandeur et des valeurs possibles de température (Ti), on détermine (en 18), pour chaque valeur possible (Ti), les probabilités P (Ti) pour que chaque valeur possible (Ti) soit la valeur exacte de la température des résonateurs ; on calcule (en 19) une valeur pondérée T de la température des résonateurs vibrants T = E P (Ti) x Ti ; et à partir d'une ou des fréquences mesurées fi et f2 et de la valeur pondérée T de la température, on détermine (en 21) la valeur compensée en température de l'accélération subie par le dispositif accéléro- métrique. 2. Procédé selon la revendication 1, caractérisé en ce que la somme des fréquences (fi + f2) est corrigée, pour tenir compte de différences géométriques entres les deux résonateurs vibrants, avec un terme d'erreur estimé à partir de l'accélération calculée lors du cycle de calcul précédent, s'il existe. 3. Procédé selon la revendication 1 ou 2, caractérisé en ce que l'influence de la température sur la valeur de la vibration de chaque résonateur vibrant est traduite par un polynôme d'ordre 2 : II (T2) ou par un polynôme d'ordre 3 : II (T3), et en ce que, dans la susdite suite d'étapes du procédé, - à partir de la corrélation préétablie (en 14,25, 26) de la variation de la somme des fréquences mesurées (fi + f2) en fonction de la température, on relève (en 13) les deux valeurs possibles de température (Ta, Tb) correspondant à la valeur de la somme des fréquences (fi + f2) ci-dessus déterminée et à partir de la valeur estimée (Gestimée) de ladite grandeur (en 17) et des deux valeurs possibles de température (Ta, Tb), on détermine (en 18) la probabilité (P (Ta) ) pour que la température des résonateurs soit l'une de ces températures (Ta) et la probabilité (P (Tb) ) pour que la température des résonateurs soit l'autre de ces températures (Tb), puis on calcule (en 19) une valeur pondérée de la température T = P (Ta). Ta + P (Tb). Tb. 4. Procédé selon la revendication 3, caractérisé en ce que la corrélation préétablie est une fonction mathématique du second ou du troisième degré, tenue en mémoire (en 14), représentant la variation de la fréquence en fonction de la température établie pour le dispositif accélérométrique considéré à deux résonateurs vibrants. 5. Procédé selon la revendication 3, caractérisé en ce que la corrélation préétablie est une table de correspondance tenue en mémoire (en 14), donnant, pour chaque valeur de la fréquence (fi + f2), les deux valeurs possibles (Ta, Tb) de la température. 6. Procédé selon l'une quelconque des revendications 1 à 5, caractérisé en ce que la fonction de transfert (H) est une fonction mathématique tenue en mémoire (en 16) et établie expérimentalement pour le dispositif accélérométrique considéré à deux poutres vibrantes. 7. Procédé selon l'une quelconque des revendications 1 à 5, caractérisé en ce que la fonction de transfert (H) est une table tenue en mémoire (en 16) et donnant une correspondance préétablie entre chaque valeur mesurée (Gmes) et une valeur estimée (Gestlmée). 8. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'on effectue (en 22 à 24) un recalage sur la fonction de transfert (H), sur la sélection des valeurs possibles (Ti) de température et sur le calcul final lorsqu'on constate une dérive de la courbe de corrélation entre la fréquence et la température. 9. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que les fréquences (fi, f2) des vibrations respectives des deux résonateurs vibrants sont mesurées. 10. Procédé de mesures accélérométriques et/ou gyroscopiques et/ou gyrométriques selon plusieurs axes à l'aide de plusieurs dispositifs vibrants accélérométriques et/ou gyroscopiques et/ou gyrométriques respectifs, caractérisé en ce que, selon chaque axe, on détermine une valeur de l'accélération compensée en température selon l'une quelconque des revendications 1 à 9, et en ce que, pour la détermination de l'une au moins des valeurs de l'accélération compensée en température, on utilise, en tant que valeur pondérée de température, la valeur pondérée de température déterminée pour le dispositif vibrant associé à cet axe ainsi qu'au moins une valeur pondérée de température déterminée pour au moins un dispositif vibrant associé à au moins un autre axe. 11. Procédé de mesure accélérométrique et/ou gyroscopique et/ou gyrométrique selon plusieurs axes selon la revendication 10, caractérisé en ce que tous les résonateurs vibrants des dispositifs accélérométriques et/ou gyroscopiques et/ou gyrométriques sont usinés dans un même bloc de cristal. 12. Procédé selon l'une quelconque des revendications 1 à 11, caractérisé en ce que la grandeur représentative de la température à proximité du dispositif à deux résonateurs vibrants est la température (Tmes) ellemême qui est mesurée (en 15) à proximité du dispositif à deux résonateurs vibrants, et en ce qu'à l'aide de la fonction de transfert H préalablement établie (en 16), on détermine (en 17) une valeur estimée Testimée = H (Tmes) de la température des résonateurs vibrants à partir de la valeur mesurée (Tmes). 13. Procédé selon la revendication 12, caractérisé en ce que la température (Tmes) est mesurée à l'extérieur d'un boîtier du dispositif qui renferme les deux résonateurs vibrants. 14. Procédé selon la revendication 12, caractérisé en ce que la température (Tmes) est mesurée à l'intérieur d'un boîtier du dispositif qui renferme les deux résonateurs vibrants. 15. Procédé selon la revendication 10 ou 11, caractérisé en ce que l'un au moins des dispositifs vibrants a une fréquence de résonance qui varie linéairement en fonction de la température et en ce que c'est cette fréquence dépendant linéairement de la température qui est mesurée et utilisée en tant que grandeur représentative de la température.