FR 2897938 A1 20070831 FR 0650637 A 20060224 Procédé d'estimation de l'adhérence d'une surface [1] L'invention se rapporte à un procédé d'estimation de l'adhérence d'une surface et plus particulièrement de vérifier ses caractéristiques structurelles et de mouillure. [2] Dans le but de sécuriser les conditions de roulage, il est un besoin du constructeur de pouvoir prévoir des réactions de ses véhicules automobiles en toutes conditions. Cependant, il est difficile de mettre au point un véhicule automobile pour toutes les situations de roulage et n'est au final qu'un io compromis entre les différentes situations instrumentées. Cela est dû notamment au fait qu'il est difficile de vérifier l'adhérence réelle du véhicule car la structure de la route ou sa mouillure peuvent considérablement varier. Par conséquent, il n'est pas possible de sécuriser les conditions de roulage 15 finement lors de la mise au point du véhicule. [3] Le but de la présente invention est de pallier tout ou partie les inconvénients cités précédemment en proposant un procédé capable de déterminer les caractéristiques d'adhérence de la route. 20 [4] A cet effet, l'invention se rapporte à un procédé d'estimation de la mouillure d'une surface comprenant l'étape consistant à établir un tableau de valeurs d'intensité des ondes électromagnétiques polarisées émises puis réfléchies sur ladite surface sous le même angle selon un 25 échantillonnage prédéterminé permettant de déterminer les caractéristiques réflectométriques de ladite surface caractérisé en ce qu'il comprend l'étape consistant à quantifier la dissymétrie de la distribution des valeurs d'intensité par rapport à une valeur de référence afin de la corréler aux 30 phénomènes de diffusion de la surface. [5] Conformément à d'autres caractéristiques avantageuses de l'invention : - la valeur de référence est la valeur moyenne arithmétique des intensités ; - l'étape de quantification consiste à calculer le moment centré d'ordre 3 de la valeur de référence à partir de la relation : 1 n-1 63 =ùL(UP(i)ùVF) n Z=o où 63 est le moment centré d'ordre 3, n est le cardinal des valeurs d'intensité, UP(i) est la valeur d'intensité du point i de l'échantillonnage et VF est la valeur de référence ; io - Le procédé comprend l'étape consistant à comparer la valeur de ladite dissymétrie à une valeur seuil pour déterminer ou non la présence d'eau sur la surface. - l'étape de comparaison est suivie d'une étape d'alerte si la valeur de ladite dissymétrie est supérieure à une valeur 15 seuil permettant de signaler la présence d'eau sur la surface. -la valeur seuil est comprise entre 5 et 8. [6] L'invention se rapporte également à un procédé d'estimation de la structure d'une surface comprenant l'étape 20 consistant à établir un tableau de valeurs de la somme des intensités des ondes polarisées et dépolarisées réfléchies sur ladite surface sous le même angle que l'émission selon un échantillonnage prédéterminé permettant de déterminer les caractéristiques réflectométriques de ladite surface 25 caractérisé en ce qu'il comprend l'étape consistant à déterminer le spectre énergétique du signal de la somme des ondes mesurées par bande de fréquence à partir des valeurs dudit tableau afin de déterminer la distribution entre la microstructure et la macrostructure de la surface. 30 [7] Conformément à d'autres caractéristiques avantageuses de l'invention : - l'étape de détermination du spectre est réalisée par calcul en utilisant la relation : D(UT) = DSP(UT)ô-1 où D(UT) est le spectre énergétique de UT en fonction de la fréquence, DSP est la fonction mathématique de densité spectrale de puissance obtenue à partir de la transformée de Fourier discrète du signal de UT et n est le nombre d'échantillons de UT; - l'étape de détermination des spectres est suivie d'une étape de détermination de l'énergie basse fréquence par calcul en utilisant la relation : ,z EBF = ù E (D(UT)(i)) N Z=,,, io où EBF est l'énergie basse fréquence de UT, DSP est la fonction mathématique de densité spectrale de puissance de UT et N est le nombre d'échantillons entre les basses fréquences v,.et v2; - l'étape de détermination de l'énergie basse fréquence 15 comporte également une étape de détermination de l'énergie haute fréquence par calcul en utilisant la relation : V4 EHF = ù E (D(UT)(i)) Ni=3 ,3 où EHF est l'énergie basse fréquence de UT, DSP est la 20 fonction mathématique de densité spectrale de puissance de UT et N est le nombre d'échantillons entre les hautes fréquences v3.et v4; - l'étape de détermination de l'énergie haute fréquence est suivie d'une étape de détermination de la répartition entre 25 la microstructure et la macrostructure de la surface par calcul en utilisant la relation : EBF Tl EHF où , représente la répartition de l'énergie, EBF est l'énergie basse fréquence et EHF est l'énergie haute 30 fréquence ; - les fréquences v,; v2; v3 et v4 sont respectivement égales à 0,01; 0,1; 1 et 5 mm-' [8] L'invention se rapporte également à un procédé d'estimation du facteur de polarisation d'une surface comprenant l'étape consistant à établir un tableau de valeurs d'intensité des ondes, respectivement polarisées et de la somme des ondes polarisées et dépolarisées, réfléchies sur ladite surface sous le même angle que l'émission selon un échantillonnage prédéterminé permettant de déterminer les caractéristiques réflectométriques de ladite surface caractérisé en ce qu'il comprend l'étape consistant à io déterminer la valeur du signal polarisé normalisé en utilisant la relation : m n- MPN = 1 ; Up(i) n-m m z 1 z m+l LUT(j) z où MPN est la moyenne du signal polarisé, Up(i) est la valeur d'intensité du point i des ondes polarisées, UT(j) est la valeur 15 d'intensité du point j de la somme des ondes polarisées et dépolarisées, n est le cardinal des valeurs d'intensité et m est un intervalle centré sur un des n points. [9] Avantageusement selon l'invention la valeur d'intervalle est constante et sensiblement compris entre 10 et 40 valeurs. 20 [10] A cet effet, l'invention se rapporte à un procédé d'estimation du facteur de dépolarisation d'une surface comprenant l'étape consistant à établir un tableau de valeurs d'intensité des ondes, respectivement dépolarisées et de la somme des ondes polarisées et dépolarisées, réfléchies sur 25 ladite surface sous le même angle que l'émission selon un échantillonnage prédéterminé permettant de déterminer les caractéristiques réflectométriques de ladite surface caractérisé en ce qu'il comprend l'étape consistant à déterminer la valeur du signal dépolarisé normalisé en 30 utilisant la relation m n- MDN = 1 ; U (i n-m Lm +m z 1 z m+1 LUT(j) z où MDN est la moyenne du signal polarisé, UD(i) est la valeur d'intensité du point i des ondes dépolarisées (Yrd), UT(i) est la valeur d'intensité du point j de la somme des ondes polarisées (yu) et dépolarisées (Yrd), n est le cardinal des valeurs d'intensité, m est un intervalle centré sur un des n points. [1 1 ] Avantageusement selon l'invention la valeur d'intervalle est constante et sensiblement compris entre 10 et 40 valeurs. [12] Enfin, l'invention se rapporte également à un procédé io d'estimation de l'adhérence d'une surface caractérisé en ce qu'il comporte les étapes consistant à calculer la mouillure selon un procédé d'estimation conforme à l'une des variantes ci-dessus, calculer la structure selon un procédé d'estimation conforme à l'une des variantes ci-dessus, calculer l'adhérence 15 par rapport aux calculs de mouillure et de structure en utilisant la relation : . = ki . + kz • ri2 + k3 où ,u représente l'adhérence de la surface, k1, k2 et k3 sont des valeurs constantes de corrélation de l'adhérence, 6-3 20 représente la mouillure de la surface, , représente la structure de la surface. [13] Conformément à d'autres caractéristiques avantageuses de l'invention : -le procédé comporte les étapes consistant à calculer 25 facteur de polarisation selon un procédé d'estimation conforme à l'une des variantes ci-dessus, calculer facteur de dépolarisation selon un procédé d'estimation conforme à l'une des variantes ci-dessus et l'étape de calcul l'adhérence est également effectuée avec les calculs des 30 facteurs de polarisation et de dépolarisation en utilisant la relation : =ki 63 +k2 .112 +k4.MPN 2+k5.MDN 2+k3 où ,u représente l'adhérence de la surface, k1, k2, k4, k5 et k3 sont des valeurs constantes de corrélation de l'adhérence, 6-3 représente la mouillure de la surface, i représente la structure de la surface, MPN représente le facteur de polarisation et MDN représente le facteur de dépolarisation. - les valeurs constantes de corrélation sont calculées par itération successives à partir de mesures d'adhérence sur io un véhicule automobile. - lesdites mesures d'adhérence (uni) sont obtenues par la relation H g Mar+M y L où y est la décélération obtenue lors du freinage d'un 15 véhicule, M est la masse totale du véhicule, Mav est la masse sur l'essieu avant, Mar est la masse sur l'essieu arrière, H est la hauteur du centre de gravité du véhicule, L est l'empattement du véhicule et g est l'accélération de la pesanteur. 20 -la mesure de décélération est effectuée en même temps que lesdites valeurs d'intensité. - la mesure de décélération est établie dans ledit tableau en même temps que lesdites valeurs d'intensité. - le procédé comporte l'étape finale de réinitialisation dudit 25 tableau. [14] D'autres particularités et avantages ressortiront clairement de la description qui en est faite ci-après, à titre indicatif et nullement limitatif, en référence aux dessins annexés, dans lesquels : 30 - la figure 1 est un organigramme du procédé selon l'invention et ; - la figure 2 est une représentation schématique du système selon l'invention. [15] Comme illustré à la figure 2, on peut voir un système de 35 caractérisation généralement annoté 1 d'une surface S. Il M•y-0,01•g•Mar comporte principalement un système d'émission 3 d'ondes électromagnétiques et un système de détection 5 d'ondes électromagnétiques et un support d'enregistrement (non représenté). Le principe utilisé par le système de caractérisation 1 consiste à envoyer des ondes électromagnétiques yep sur la surface S puis récupérer les ondes électromagnétiques réfléchies yrp, yYd sur la surface S. Une des premières caractéristiques selon l'invention réside en ce que, préférentiellement, les directions d'émissions E et de io détection D forme un angle a identique dans le sens trigonométrique par rapport à la normale N. Selon l'invention, l'angle a est constant et est compris entre 15 et 25 . Préférentiellement, a est égale à 20 . [16] Le système d'émission 3 comporte principalement une 15 source d'ondes électromagnétiques, un dispositif de collimation et un dispositif de focalisation. Tous ces moyens sont réunis en une seule référence 7. Préférentiellement, la source d'ondes électromagnétiques est une source laser pour permettre de rayonner selon une puissance suffisante sur la 20 surface S. Le dispositif de collimation est utilisé pour émettre un rayonnement parallèle. Enfin, le dispositif de focalisation est utilisé pour borner la section du rayonnement yep à un diamètre sensiblement égale à 400 ,um. [17] Selon l'invention, le système d'émission 3 rayonne, de 25 façon continue, sur la surface S des ondes électromagnétiques polarisées yep selon une longueur d'onde îs constante. Préférentiellement, selon l'invention, la longueur d'onde est égale à 980 nm car ce type de laser est très répandu dans le commerce. 30 [18] Le système de détection 5 comporte principalement un dispositif séparateur de polarisation 9, un dispositif de filtrage 11 et un dispositif de déplacement 13. Le dispositif de filtrage 11 consiste essentiellement en un filtre optique qui limite le bruit engendré par les rayonnements extérieurs au système 35 d'émission 3. Il est constitué de moyens lenticulaires teintés capables de limiter les longueurs d'onde des rayonnements qui le traversent. Préférentiellement selon l'invention, le dispositif de filtrage 11 est capable de laisser passer la longueur d'onde îs du système d'émission 3 avec une tolérance de 50 nm. Dans l'exemple cité ci-dessus, cela correspondrait à une plage de valeurs comprise entre 930 et 1030 nm. [19] Le dispositif séparateur de polarisation 9 est utilisé en aval du dispositif de filtrage 11 pour dissocier les ondes réfléchies y, selon qu'elles soient polarisées yYp ou dépolarisées yrd. Pour détecter et mesurer ces rayonnements io dépolarisés et polarisés yrd, yYp, on utilise, respectivement, des moyens de détection 15 et 17. Ces moyens sont, selon l'invention, des transducteurs du type des photodiodes capables de quantifier des intensités lumineuses en un signal électrique sensiblement proportionnel (UD, Up). Le système de 15 détection 5 peut ainsi détecter les ondes réfléchies polarisées yYp et/ou dépolarisées yrd provenant de la surface S. Comme visible à la figure 2, on peut voir que ces ondes ne peuvent provenir cependant jusqu'aux moyens de détection 15 et 17 que si elles restent dans le cylindre de détection C centré sur 20 la direction D. Par conséquent, les autres ondes électromagnétiques ne seront pas détectées. [20] Préférentiellement selon l'invention, le système de détection 5 comporte également un dispositif de déplacement 13. Selon l'invention, le dispositif de déplacement 13 consiste 25 en un véhicule automobile. Cela permet lorsqu'il se déplace selon une direction X d'émettre des rayonnements selon la direction E et de détecter ces rayonnements selon la direction D, sensiblement équivalente à E, provenant de la portion de la surface S quelques instants plus tard. Bien entendu, la 30 distance, entre les systèmes d'émission 3 et de détection 5 d'une part, et la surface S d'autre part, est liée la garde au sol du véhicule automobile. [21] Le principe d'utiliser un véhicule automobile comme dispositif de déplacement 13 a nécessité des aménagements 35 quant à la structure du système de caractérisation 1. Ainsi, il comporte spécifiquement un dispositif de variation de fréquence d'échantillonnage 19 et des moyens de protection 21. [22] Les moyens de protection 21 comportent essentiellement un hublot 23 incliné selon le même angle que les directions E et D par rapport à la normale N à la surface S. Ces moyens sont utilisés pour éviter les projections venant de la surface S qui peut être une route ou un chemin empruntés par le véhicule automobile. Les moyens de protection 21 sont inclinés vers l'avant du véhicule automobile afin de minimiser io les éclaboussures sur la section intérieure utile transparente et maximiser leur contact par rapport au bord distal inférieur opaque. L'opacité est préférée afin de diminuer d'avantage le bruit optique, c'est-à-dire les rayonnements ne venant pas du système d'émission 3. On comprend alors que le rayon de 15 la section minimale transparente du hublot est liée à l'éloignement entre les directions E et D. [23] Le dispositif de variation de fréquence d'échantillonnage 19 est monté à même les moyens de détection 15 et 17. Il est utilisé pour que, en fonction de la vitesse Cv de déplacement 20 du véhicule automobile, le pas d'échantillonnage Al reste constant. Préférentiellement selon l'invention, la fréquence F est calculée à partir de la relation : F=C Al où 25 -F est la fréquence d'échantillonnage d'émission ; - C, est la vitesse du véhicule ; - Al est le pas d'échantillonnage. [24] Préférentiellement selon l'invention, le pas d'échantillonnage Al est compris entre 25 et 100 ,um. De plus, 30 de manière préférée, le système de caractérisation 1 comporte un support d'enregistrement capable de stocker des valeurs détectées sur une distance d'échantillonnage de la surface S sensiblement compris entre 0,5 et 2 m le long de la direction X. i0 [25] Nous allons maintenant expliquer le procédé de caractérisation de la surface S généralement annoté 31 à la figure 1. le procédé de caractérisation 31 comporte une étape d'émission 33 d'ondes électromagnétiques et une étape de détection 35. L'étape d'émission 33 consiste à émettre du véhicule des ondes électromagnétiques yep polarisées. Cela est utile notamment pour les comparées à celles qui seront réfléchies sur la surface S. Comme expliqué lors de l'explication du système de caractérisation 1, l'étape io d'émission 33 n'est pas subordonnée à un quelconque pilotage, elle se borne à émettre en continue des ondes électromagnétiques yep sous l'angle a par rapport à la normale N. L'étape peut cependant être asservie à une commande de mise en fonctionnement de la fonction de caractérisation ou du 15 démarrage du véhicule automobile. [26] L'étape de détection 35 comporte principalement des étapes de mesure 37, de variation de fréquence 39, d'enregistrement 41 et de comparaison 43. Comme expliqué lors de l'explication du système de caractérisation 1, l'étape 20 de mesure 37 des ondes électromagnétiques réfléchies sur la surface S permet de détecter celles qui sont polarisées yYp et/ou dépolarisées yrd. La mesure est réalisée selon le même angle a par rapport à la surface S que celles émises yep afin de déterminer les caractéristiques réflectométriques de la surface 25 S. L'étape de mesure 37 relève ainsi les intensités des ondes électromagnétiques réfléchies yrd, yYp qui sont quantifiées respectivement au moyen de la tension des signaux UD, Up provenant des moyens de détection 15 et 17. [27] Selon une variante de l'invention, on mesure également la 30 décélération y, lors de l'étape 37, du véhicule automobile si le véhicule est en phase de freinage. 2897938 Il [28] Comme expliqué lors de l'explication du système de caractérisation 1, le fait que le dispositif de déplacement 13 soit le véhicule automobile entraîne un pas d'échantillonnage Al qui varie suivant sa vitesse Cv. Il a donc été imaginé que 5 l'étape de mesure 37 soit asservie à une fréquence F d'échantillonnage variable pour effectuer une mesure à distance régulière Al , par exemple 50 ,um, le long de l'axe X sur la surface S. [29] Préférentiellement selon l'invention, la fréquence F est io calculée lors de l'étape 39 à partir de la relation : F=C Al où - F est la fréquence d'échantillonnage d'émission ; - C, est la vitesse du véhicule ; 15 -Al est le pas d'échantillonnage. [32] A la suite de l'étape de mesure 37, il est prévue une étape d'enregistrement 41 dans un tableau des mesures UD et/ou UP et/ou leur somme UT (c'est-à-dire UD+UP) et/ou la décélération y du véhicule automobile en fonction d'une désignation de la 20 mesure. Préférentiellement, toutes les valeurs UD, UT, UT et y sont enregistrées dans le tableau lors de l'étape 41 sur le support d'enregistrement du système de caractérisation 1. Cependant, dans un but de simplifier le procédé 31, une ou plusieurs valeur(s) parmi UD, UT, UT et y peuvent ne pas être 25 enregistrées. [33] Préférentiellement, la désignation de la mesure est un chiffre entier qui est incrémenté à chaque étape de mesure 37. Cependant, elle peut correspondre à la distance parcourue le long de l'axe X depuis la première étape de mesure 37. 30 [32] A la suite de l'étape d'enregistrement 41, il est prévu une étape de comparaison 43 par rapport à ladite désignation de mesure. En effet selon l'invention, l'étape de mesure 37 et d'enregistrement ne sont répétées que sur une distance d'échantillonnage prédéfinie, par exemple 1 m. Dès lors suivant l'une des méthodes de désignation de la mesure expliquées ci-dessus, on utilisera respectivement le cardinal des mesures présentes dans le tableau ou directement le quantième de la désignation. Lors de l'étape 43, si la désignation de la mesure est en dessous d'un nombre seuil (c'est-à-dire indique que l'on est en dessous de 1 m d'échantillonnage) alors le procédé retourne à l'étape de mesure 37 sinon le procédé continue. [33] Bien entendu, il peut être prévu une interface homme û Io machine afin de pouvoir moduler le pas constant d'échantillonnage Al ou la distance constante de caractérisation le long de l'axe X. Ils peuvent alors être compris respectivement entre 25 et 100 ,um et entre 0,5 et 2 m. [34] Nous allons maintenant expliquer le procédé d'estimation 15 de la mouillure 51 de la surface S généralement annoté 51 à la figure 1. Le procédé d'estimation de la mouillure 51 utilise les valeurs du tableau établit à l'aide du procédé de caractérisation 31 expliqué ci-dessus. Avantageusement selon l'invention pour déterminer la mouillure (c'est-à-dire détecter 20 si la surface S est humide), le procédé d'estimation 51 utilise uniquement les valeurs Up dudit tableau, c'est-à-dire les valeurs d'intensités des ondes électromagnétiques réfléchies de la surface S qui sont polarisées yrp. [35] Le procédé d'estimation 51 vise notamment à vérifier la 25 distribution des valeurs Up en fonction de la distance d'échantillonnage. Le but de cette routine statistique est de vérifier qu'à l'état sec, aucune direction de réflexion sur la surface S n'est privilégiée. Cela signifie que la réflexion des ondes électromagnétiques est dite diffuse à l'opposé de la 30 réflexion dite spéculaire dans le cas de présence d'eau qui suit les lois de l'optique géométrique (également assimilée à la réflexion sur un miroir). [36] En effet dans le cas d'une réflexion diffuse, les valeurs Up sont toutes sensiblement identiques tout au long de 35 l'échantillonnage. En revanche dans le cas de présence d'eau, la réflexion dite spéculaire a tendance à privilégier des directions de réflexions qui vont évoluées au cours de l'échantillonnage. C'est pourquoi, préférentiellement, le procédé de caractérisation 31 effectue une distance minimale d'échantillonnage comme par exemple 1 m. [37] Selon le procédé d'estimation 51, on vérifie la dissymétrie des valeurs Up par rapport à une valeur de référence VF prédéfinie afin de déterminer si les valeurs présentes dans le tableau offre une répartition écartée ou sensiblement centrée. Préférentiellement selon l'invention, la valeur de référence VF est calculée lors de l'étape 53 et correspond à la moyenne arithmétique des valeurs Up, ci-après nommé Up . [38] A la suite de l'étape 53, on calcule lors de l'étape 55 le moment centré d'ordre 3. En effet comme expliqué ci-dessus, lors d'une réflexion diffuse, les valeurs Up ont tendance à être sensiblement égales. En prenant la valeur moyenne Up puis en vérifiant la dispersion de la répartition Up par rapport à cette valeur moyenne, on trouve que ladite dispersion est faible. Préférentiellement, le moment centré d'ordre 3 o de la valeur de référence est calculé à partir de la relation : -1 _ 63 =' L(UP(i)ùUP) n i=o où : - 63 est le moment centré d'ordre 3 ; - n est le cardinal des valeurs d'intensité ; -Up(i) est la valeur d'intensité du point i de l'échantillonnage ; - Up est la valeur de référence. [39] Selon une variante de l'invention dans laquelle on cherche uniquement à détecter l'humidité sur une surface S, le procédé d'estimation de la mouillure 51 comporte, à la suite de l'étape 55, l'étape 57 de comparaison et d'alerte de la valeur a' obtenue à une valeur seuil 6s qui est représentative de la valeur empirique au-dessus de laquelle il est certain qu'il y a une présence d'eau sur la surface. L'étape 57 prévient alors le 35 conducteur du véhicule qu'un fluide est présent sur la route. [40] Bien entendu, la valeur seuil 6s est à déterminer préalablement, cependant, à l'aide d'essais réalisés sur nombre de bitumes, la valeur 6s constante est généralement comprise entre 5 et 8. Préférentiellement selon l'invention, elle est égale à 6. [41] Nous allons maintenant expliquer le procédé d'estimation de la structure i de la surface S généralement annoté 61 à la figure 1. Le procédé d'estimation de la structure 61 utilise les valeurs du tableau établit à l'aide du procédé de caractérisation 31 expliqué ci-dessus. Avantageusement selon l'invention pour déterminer la structure (c'est-à-dire détecter le rapport entre microstructure et macrostructure sur la surface S), le procédé d'estimation 61 utilise uniquement les valeurs UT dudit tableau, c'est-à-dire la somme des valeurs d'intensités UP+UD des ondes électromagnétiques réfléchies polarisées yu et dépolarisées Yrd. En effet, il est possible à l'aide de ces valeurs de vérifier quelle en est la répartition. [42] Le procédé d'estimation 61 comporte une première étape 63 consistant à déterminer le spectre énergétique du signal de la somme des ondes mesurées (UT) par bande de fréquence spatiale. Cela permet de se faire une représentation par fréquence de la répartition du signal pour les faire correspondre aux critères de micro- et macrostructure. [43] Selon l'invention, l'étape 63 de détermination du spectre est réalisée par calcul en utilisant la relation : D(UT) = DSP(UT)ô-1 où - D(UT) est le spectre énergétique de UT en fonction de la fréquence ; -DSP est la fonction mathématique de densité spectrale de puissance obtenue à partir de la transformée de Fourier discrète du signal de UT ; - n est le nombre d'échantillons de UT. [44] Préférentiellement selon l'invention, il est décidé que la microstructure représente des grains dont le diamètre est sensiblement inférieur ou égal à 380 ,um et, pour la macrostructure, sensiblement supérieur ou égal à 500 ,um. Préférentiellement par rapport à la densité spectrale de puissance calculée D(UT), cela revient à étudier, d'une part pour la microstructure, les bandes de fréquence v3.et v4 sensiblement entre 1 et 5 mm-' et, d'autre part pour la macrostructure, les bandes de fréquence vl.et v2 sensiblement entre 0,01 et 0,1 mm-1. [45] C'est pourquoi, dans une deuxième étape 65, à partir de la densité spectrale de puissance calculée D(UT), on discrimine les deux plages de fréquences, c'est-à-dire vl-v2 et v3-v4. pour déterminer respectivement l'énergie basse fréquence EBF et l'énergie haute fréquence EHF. [46] Préférentiellement selon l'invention, l'énergie basse fréquence est obtenue par calcul en utilisant la relation ,z EBF = ù E (D(UT)(i)) N Z=,,, où - EBF est l'énergie basse fréquence de UT ; - DSP est la fonction mathématique de densité spectrale de puissance de UT ; - N est le nombre d'échantillons entre les basses fréquences vl.et V2. [47] Par conséquent, l'énergie haute fréquence est obtenue consécutivement par un calcul similaire en utilisant la relation : 1 V4 E H F = ù E (D(UT )(i)) N où - EHF est l'énergie basse fréquence de UT ; - DSP est la fonction mathématique de densité spectrale de puissance de UT ; - N est le nombre d'échantillons entre les hautes fréquences v3.et V4. [48] On remarque que, selon l'invention, on réalise une moyenne arithmétique sur la plage de fréquence considérée pour obtenir les valeurs d'énergie EBF et EHF. [49] A la suite de l'étape 65, on estime la structure Il par une étape 67 de détermination de la répartition entre la microstructure et la macrostructure de la surface. Préférentiellement, l'étape 67 est réalisée par calcul en utilisant la relation _ EBF EHF où: - ~ représente la répartition de l'énergie ; - EBF est l'énergie basse fréquence ; - EHF est l'énergie haute fréquence. [53] L'étape 67 permet donc de se faire une idée de la répartition qui est, au travers d'essais sur nombre de bitumes, systématiquement majoritairement sous forme demicrostructure, c'est-à-dire que q, le rapport entre EBF et EHF, est toujours inférieur à 1. [54] Nous allons maintenant expliquer le procédé d'estimation de la polarisation et de la dépolarisation de la surface S respectivement annoté 71 et 75 à la figure 1. Les facteurs de polarisation MPN et de dépolarisation MDN sont liés en ce qu'ils utilisent le même calcul mais à partir des valeurs différentes d'intensité du tableau obtenu, par exemple, par le procédé de caractérisation 31. En effet, pour la détermination du facteur de polarisation MPN, on utilise les valeurs d'intensité Up des ondes électromagnétiques réfléchies polarisées yu et de la somme UT de toutes les ondes électromagnétiques réfléchies yr. On comprend alors que pour la détermination du facteur de dépolarisation MDN, on utilise les valeurs d'intensité UD des ondes électromagnétiques réfléchies dépolarisées Yrd et de la somme UT de toutes les ondes électromagnétiques réfléchies yr. [52] Préférentiellement selon l'invention, le procédé d'estimation du facteur de polarisation 71 comporte une étape principale 73 consistant à déterminer la valeur du signal polarisé normalisé, c'est-à-dire MPN, en utilisant la relation : ( où m n- 1 2 MPN = n-m 2u (i Z 1 m+1 EUTU) 2 - MPN est la moyenne du signal polarisé normalisé ; - UP(i) est la valeur d'intensité du point i des ondes polarisées yu ; io - UT(j) est la valeur d'intensité du point j de la somme des ondes polarisées yu et dépolarisées yrd - n est le cardinal des valeurs d'intensité - m est un intervalle centré sur un des n points. [53] Comme expliqué ci-dessus, de façon similaire le procédé 15 d'estimation du facteur de dépolarisation 75 comporte une étape principale 77 consistant à déterminer la valeur du signal dépolarisé normalisé, c'est-à-dire MDN, en utilisant la relation : m n- 1 2 MDN = n ù m 2UD (i m 2 1 m+1 EUTU) 2 20 Oë - MDN est la moyenne du signal polarisé ; - UD(i) est la valeur d'intensité du point i des ondes dépolarisées yrd ; - UT(j) est la valeur d'intensité du point j de la somme des 25 ondes polarisées yu et dépolarisées yrd - n est le cardinal des valeurs d'intensité - m est un intervalle centré sur un des n points. [54] Préférentiellement selon l'invention, les intervalles m sont égaux à 20 points, cependant ils peuvent être compris entre 10 et 40 points. Cela signifie que l'on effectue une double moyenne arithmétique pour obtenir MPN ou MDN (sur m/2 autour du point n pour UT et sur les points n pour respectivement Up et (ID). [55] Nous allons maintenant expliquer le procédé d'estimation de l'adhérence ,u de la surface S généralement annoté 81 à la figure 1. Le procédé 81 représente une synthèse des procédés io d'estimation de la mouillure 51 et/ou de la structure 61 et/ou du facteur de polarisation 71 et/ou du facteur de dépolarisation. Il consiste essentiellement à corréler les différentes valeurs calculées par rapport à une adhérence mesurée ,um pour estimer a posteriori l'adhérence ,u de la 15 surface S échantillonnée. [56] Préférentiellement selon l'invention, l'adhérence ,u, lors de l'étape 87, est obtenue en utilisant la relation : =ki•63+k2.i 2+k4•MPN2+k5•MDN2+k3 où 20 -,u représente l'adhérence de la surface ; - k1 k2, k3, k4 et k5 sont des valeurs constantes de corrélation de l'adhérence ; - 63 représente la mouillure de la surface ; - , représente la structure de la surface ; 25 -MPN représente le facteur de polarisation ; - MDN représente le facteur de dépolarisation. [61] On remarque que le facteur de mouillure 6-3 est privilégié par rapport aux autres valeurs calculées. En effet, les autres valeurs, qui sont toujours inférieures à 1, sont toutes montées 30 au carré afin de minimiser leur influence. C'est un parti pris en faveur de la mouillure qui a été jugée comme le facteur le plus important pour l'adhérence, cependant, il peut être envisagé de laisser toutes les valeurs à la même puissance. [62] Selon l'invention, les valeurs constantes de corrélation kl, 35 k2, k3, k4 et k5 sont à calculer, lors de l'étape 85, afin que le procédé d'estimation 81 puisse déterminer l'adhérence ,u. Les valeurs de corrélations sont obtenues par itérations successives par rapport à une mesure de l'adhérence ,um effectuée lors de l'étape 83. [59] Préférentiellement selon l'invention, l'étape de mesure 83 de l'adhérence ,um est obtenue à partie de la relation : II M•yù0,01.g•Mar H g .Mar +M .y L où - y est la décélération obtenue lors du freinage d'un véhicule ; - M est la masse totale du véhicule ; - Mav est la masse sur l'essieu avant ; - Mar est la masse sur l'essieu arrière ; - H est la hauteur du centre de gravité du véhicule ; -L est l'empattement du véhicule ; - g est l'accélération de la pesanteur. [60] On comprend que le procédé d'estimation 81 nécessite une liaison avec le procédé de caractérisation 31 afin de récupérer la valeur y de décélération obtenue lors du freinage du véhicule. Selon une première variante, après chaque étape 37 de mesure, la valeur y est transmise à l'étape 83 qui mesure l'adhérence réelle. Selon une seconde variante, lors de l'étape 41, la valeur y est enregistrée sur ledit tableau en même temps que les valeurs d'intensités Up et/ou UD et/ou UT. L'étape 83 n'a alors plus qu'à la récupérer. [61] Il faut comprendre de ce qui expliqué ci-dessus qu'il n'est pas possible de réaliser des mesures d'adhérence ,um en permanence mais uniquement en phase de freinage et selon des conditions particulières. C'est pourquoi, selon l'invention, les étapes 83 et 85 sont préférentiellement réalisée en sortie de ligne de montage du véhicule automobile afin d'obtenir des valeurs k1, k2, k3, k4 et k5 parfaitement calibrées. Un exemple de calibration fait correspondre à kl; k2; k3; k4 et k5, les valeurs numériques respectives -0,005; 9,8; 1,14; -0,97 et -0,45. [62] Selon l'invention à la suite de l'étape de calcul 87 du procédé d'estimation 81, une étape finale 89 de réinitialisation du tableau de valeur est effectuée. L'étape 89 consiste essentiellement en un effacement des valeurs présentes dans le tableau afin d'effectuer un nouvel échantillonnage à l'aide du procédé de caractérisation 31. [63] Bien entendu, la présente invention ne se limite pas à l'exemple illustré mais est susceptible de diverses variantes et modifications qui apparaîtront à l'homme de l'art. En io particulier, le système de détection 5 peut être non pas monté sur un véhicule automobile mais à un endroit fixe. Dans ce cas, le dispositif de déplacement 13 peut alors consister en des moyens du type crémaillère û pignon sur lesquels sont montés les systèmes d'émission 3 et de détection 5 et dont le 15 rayon de courbure de la crémaillère correspond à la distance entre la surface S et ledit système de détection. [64] Il peut également être prévu que selon une variante de l'invention, les valeurs constantes de corrélation k1, k2, k3, k4 et k5 sont préenregistrées et utilisée directement par le procédé 20 d'estimation 81 pour déterminer l'adhérence ,u sans avoir à effectuer de mesure de décélération lors de freinage. [65] Selon une autre variante, l'étape de calcul 87 de l'adhérence ,u n'utilise pas tous les facteurs et se contente d'utiliser la relation 25 =k 63 où -,u représente l'adhérence de la surface ; - k1 k2 et k3 sont des valeurs constantes de corrélation de l'adhérence ; 30 -63 représente la mouillure de la surface ; - q représente la structure de la surface ; 20 REVENDICATIONS 1. Procédé d'estimation de la mouillure (51) d'une surface (S) comprenant l'étape (31) consistant à établir un tableau de valeurs d'intensité (UP) des ondes électromagnétiques polarisées émises (yu) puis réfléchies (yu) sur ladite surface sous le même angle (a) selon un échantillonnage prédéterminé permettant de déterminer les caractéristiques réflectométriques de ladite surface caractérisé en ce qu'il comprend l'étape (55) consistant à quantifier la dissymétrie de la distribution des valeurs d'intensité (UP) par rapport à une valeur de référence (VF) afin de la corréler aux phénomènes de diffusion de la surface. 2. Procédé (51) selon la revendication 1, caractérisé en ce que la valeur de référence (VF) est la valeur moyenne arithmétique des intensités (Up ). 3. Procédé d'estimation selon la revendication 1 ou 2, caractérisé en ce que l'étape (55) de quantification consiste à calculer le moment centré d'ordre 3 (63) de la valeur de référence à partir de la relation : n-1 6 3 = 1 L (UP (i) ù VF ) i=o où : - 63 est le moment centré d'ordre 3 ; - n est le cardinal des valeurs d'intensité ; - UP(i) est la valeur d'intensité du point i de 25 l'échantillonnage ; - VF est la valeur de référence. 4. Procédé (51) selon la revendication 2 ou 3, caractérisé en ce qu'il comprend l'étape (57) consistant à comparer la valeur (63) de ladite dissymétrie à une valeur seuil (6s) pour 30 déterminer ou non la présence d'eau sur la surface. 5. Procédé (51) selon la revendication 4, caractérisé en ce que l'étape (57) de comparaison est suivie d'une alerte si la valeur (63) de ladite dissymétrie est supérieure à une valeur seuil (as) permettant de signaler la présence d'eau sur la surface. 6. Procédé (51) selon la revendication 5, caractérisé en ce que la valeur seuil (as) est comprise entre 5 et 8. 7. Procédé d'estimation (61) de la structure (ri) d'une surface (S) comprenant l'étape (31) consistant à établir un tableau de io valeurs de la somme (UT) des intensités des ondes polarisées (yu) et dépolarisées (Yrd) réfléchies sur ladite surface sous le même angle (a) que l'émission selon un échantillonnage prédéterminé permettant de déterminer les caractéristiques réflectométriques de ladite surface caractérisé en ce qu'il 15 comprend l'étape (63) consistant à déterminer le spectre énergétique du signal de la somme des ondes mesurées (UT) par bande de fréquence à partir des valeurs dudit tableau afin de déterminer la distribution entre la microstructure et la macrostructure de la surface. 20 8. Procédé (61) selon la revendication 7, caractérisé en ce que l'étape (63) de détermination du spectre est réalisée par calcul en utilisant la relation : D(UT) = DSP(UT)ô-1 où 25 -D(UT) est le spectre énergétique de UT en fonction de la fréquence ; -DSP est la fonction mathématique de densité spectrale de puissance obtenue à partir de la transformée de Fourier discrète du signal de UT ; 30 -n est le nombre d'échantillons de UT. 9. Procédé (61) selon la revendication 8, caractérisé en ce que l'étape (63) de détermination des spectres est suivie d'une étape (65) de détermination de l'énergie basse fréquence par calcul en utilisant la relation,z EBF = ù E (D(UT )(i)) Ni,, où - EBF est l'énergie basse fréquence de UT ; - DSP est la fonction mathématique de densité spectrale de puissance de UT ; - N est le nombre d'échantillons entre les basses fréquences vl.et v2. 10. Procédé (61) selon la revendication 9, caractérisé en ce que l'étape (65) de détermination de l'énergie basse fréquence io comporte également une étape de détermination de l'énergie haute fréquence par calcul en utilisant la relation : 1 V4 EHF = ù E (D(UT)(i)) N. où - EHF est l'énergie basse fréquence de UT ; 15 -DSP est la fonction mathématique de densité spectrale de puissance de UT ; - N est le nombre d'échantillons entre les hautes fréquences v3.et V4. 11 . Procédé (61) selon la revendication 10, caractérisé en ce 20 que l'étape (65) de détermination de l'énergie haute et basse fréquence est suivie d'une étape (67) de détermination de la répartition entre la microstructure et la macrostructure de la surface par calcul en utilisant la relation EBF EHF 25 où : - il représente la répartition de l'énergie ; - EBF est l'énergie basse fréquence ; - EHF est l'énergie haute fréquence. 12. Procédé (61) selon l'une des revendications 9 à 11 , 3o caractérisé en ce que les fréquences vl.et v2 sont respectivement égales à 0,01 mm-' et 0,1 mm-1. 13. Procédé (61) selon l'une des revendications 10 à 12, caractérisé en ce que les fréquences v3.et v4 sont respectivement égales à 1 mm-' et 5 mm-1. 14. Procédé d'estimation (71) du facteur de polarisation (MPN) d'une surface (S) comprenant l'étape (31) consistant à établir un tableau de valeurs d'intensité (Up, UT) des ondes, respectivement polarisées (yu) et de la somme des ondes polarisées (yu) et dépolarisées (yrd), réfléchies sur ladite surface sous le même angle (a) que l'émission selon un échantillonnage prédéterminé permettant de déterminer les caractéristiques réflectométriques de ladite surface caractérisé en ce qu'il comprend l'étape (73) consistant à déterminer la valeur du signal polarisé normalisé (MPN) en utilisant la relation : où m n- 1 2 MPN = n-m 2u (i) Z 1 m+1 EUT Ci) 2 - MPN est la moyenne du signal polarisé ; - UP(i) est la valeur d'intensité du point i des ondes polarisées (yu) ; -U TU) est la valeur d'intensité du point j de la somme des ondes polarisées (yu) et dépolarisées (yrd) ; - n est le cardinal des valeurs d'intensité ; - m est un intervalle centré sur un des n points. 15. Procédé (71) selon la revendication 14, caractérisé en ce 25 que la valeur d'intervalle (m) est constante et sensiblement compris entre 10 et 40 valeurs. 16. Procédé d'estimation (75) du facteur de dépolarisation (MDN) d'une surface comprenant l'étape (31) consistant à établir un tableau de valeurs d'intensité (UD, UT) des ondes, respectivement dépolarisées (yrd) et de la somme des ondes polarisées (yu) et dépolarisées (yrd), réfléchies sur ladite surface sous le même angle que l'émission selon un échantillonnage prédéterminé permettant de déterminer les caractéristiques réflectométriques de ladite surface caractérisé en ce qu'il comprend l'étape (77) consistant à io déterminer la valeur du signal dépolarisé normalisé (MDN) en utilisant la relation : où m n- 1 2 MDN = n -m 2UD (i m Z 1 m+1 LU (j) 2 - MDN est la moyenne du signal polarisé ; 15 - UD(i) est la valeur d'intensité du point i des ondes dépolarisées (yrd) ; - UT(j) est la valeur d'intensité du point j de la somme des ondes polarisées (yu) et dépolarisées (yrd) ; - n est le cardinal des valeurs d'intensité ; 20 -m est un intervalle centré sur un des n points. 17. Procédé (75) selon la revendication 16, caractérisé en ce que la valeur d'intervalle (m) est constante et sensiblement compris entre 10 et 40 valeurs. 18. Procédé d'estimation (81) de l'adhérence Cu) d'une surface 25 caractérisé en ce qu'il comporte les étapes consistant à : - calculer la mouillure (a3) selon un procédé d'estimation (51) conforme à l'une des revendications 3 à 6 ; - calculer la structure (ii) selon un procédé d'estimation (61) conforme à l'une des revendications 11 à 13 ; 30 -calculer (87) l'adhérence Cu) par rapport aux calculs de mouillure (a3) et de structure (ii) en utilisant la relation : =k .63 +k2+k3 où - ,u représente l'adhérence de la surface - k1 k2 et k3 sont des valeurs constantes de corrélation de l'adhérence ; - 63 représente la mouillure de la surface ; - q représente la structure de la surface. 19. Procédé (81) selon la revendication 18, caractérisé en ce qu'il comporte les étapes consistant à : io -calculer facteur de polarisation (MPN) selon un procédé d'estimation (71) conforme à l'une des revendications 14 à 15 - calculer facteur de dépolarisation (MDN) selon un procédé d'estimation (75) conforme à l'une des revendications 16 à 15 17 -et en ce que l'étape (87) de calcul l'adhérence Cu) est également effectuée avec les calculs des facteurs de polarisation (MPN) et de dépolarisation (MDN) en utilisant la relation : 20 =k,•G3+k2 r~2+k4•MPN2+k5•MDN2 où - ,u représente l'adhérence de la surface - k1 k2, k4, k5 et k3 sont des valeurs constantes de corrélation de l'adhérence ; 25 -63 représente la mouillure de la surface ; - q représente la structure de la surface - MPN représente le facteur de polarisation ; - MDN représente le facteur de dépolarisation. 20. Procédé (81) selon la revendication 18 ou 19, caractérisé 30 en ce que les valeurs constantes de corrélation (k1, k2, k4, k5 et k3) sont calculées (85) par itération successives à partir de mesures (83) d'adhérence (,um) sur un véhicule automobile. 21 . Procédé (81) selon la revendication 20, caractérisé en ce que lesdites mesures (83) d'adhérence (,um) sont obtenues par 35 la relation : M•yù0,01. g•Mar où - y est la décélération obtenue lors du freinage d'un véhicule ; -M est la masse totale du véhicule ; - Mav est la masse sur l'essieu avant ; - Mar est la masse sur l'essieu arrière ; - H est la hauteur du centre de gravité du véhicule ; - L est l'empattement du véhicule ; io -g est l'accélération de la pesanteur. 22. Procédé (81) selon la revendication 21, caractérisé en ce que la mesure de décélération (y) est effectuée en même temps (37) que lesdites valeurs d'intensité (Up, UD, UT). 23. Procédé (81) selon la revendication 21 ou 22, caractérisé 15 en ce que la mesure de décélération (y) est établie (41) dans ledit tableau en même temps que lesdites valeurs d'intensité (Up, UD, UT). 24. Procédé (81) selon l'une des revendications précédentes, caractérisé en ce qu'il comporte l'étape (89) finale de 20 réinitialisation dudit tableau. H O Mar +M y L