Brevet FR2910662A1 — publié le 20080627

PROCEDE POUR SELECTIONNER UNE CONFIGURATION OPMIMISEE DE PNEUMATIQUES POUR UN VEHICULE TERRESTRE TEL QU'UNE VOITURE OU UNE MOTO DE COMPETITION. L'invention concerne, de façon générale, les techniques liées à l'équipement en pneumatiques des véhicules à moteurs. Plus précisément, l'invention concerne un procédé pour sélectionner, parmi un ensemble de configurations possibles, une configuration optimisée de pneumatiques pour équiper un véhicule destiné à courir sur un circuit en suivant une trajectoire prédéterminée, orientée dans un sens unique de parcours et incluant des virages. Cette sélection, qui est traditionnellement réalisée sur la base d'essais et de résultats d'expériences, représente aujourd'hui des investissements considérables en termes de temps et d'argent. Néanmoins, même dans le cas où ces investissements peuvent être réalisés, ce mode de sélection présente l'inconvénient qu'il ne peut être mis en oeuvre qu'après fabrication des pneus utilisés pour les essais. La présente invention a précisément pour but de proposer un procédé de sélection exempt de cette contrainte. A cette fin, le procédé de l'invention, par ailleurs conforme à la définition qu'en donne le préambule ci-dessus, est essentiellement caractérisé en que qu'il comprend au moins les étapes consistant à : -développer un modèle physique du véhicule équipé 30 d'une première configuration de pneus; 2910662 2 - simuler le parcours de la trajectoire par le véhicule modélisé, en optimisant la vitesse de ce véhicule; -mémoriser le temps de parcours mis par le véhicule 5 pour parcourir au moins une fois la trajectoire; - comparer le temps de parcours à au moins un temps de référence et produire un résultat de comparaison; et -sélectionner ou non la première configuration comme configuration optimisée en fonction du résultat de 10 comparaison. Dans un mode de réalisation préféré de l'invention, le modèle du véhicule inclut un modèle thermomécanique des pneus équipant ce véhicule. Ce procédé peut notamment comprendre les étapes 15 consistant à : (a) recenser un ensemble de paramètres physiques tendant à limiter la vitesse adoptée par le véhicule sur la trajectoire, cet ensemble comprenant des paramètres constants et des paramètres dynamiques, ces paramètres 20 étant liés aux différents pneus possibles, au châssis du véhicule et/ou aux conditions de roulement, et incluant une image numérique de la trajectoire, constituée par des coordonnées numérisées de points successifs de cette trajectoire; 25 (b) développer un modèle de pneu décrivant le comportement physique de chaque pneu et faisant intervenir un premier sous-ensemble de l'ensemble des paramètres physiques; (c) développer un modèle de véhicule faisant intervenir à la fois le modèle de pneu et un deuxième sous-ensemble de l'ensemble des paramètres physiques, et décrivant le comportement physique du véhicule en 2910662 3 fonction au moins de sollicitations des pneus, ces sollicitations comprenant des glissements et des dérives; (d) découper l'image numérique de la trajectoire en tronçons adjacents se succédant dans le sens de parcours 5 et dont chacun comprend un ensemble d'au moins trois points de trajectoire incluant un point initial et un point terminal de ce tronçon se succédant dans cet ordre dans le sens de parcours, le point terminal étant localisé dans une entrée de virage; 10 (e) associer respectivement aux points initial et terminal de chaque tronçon des valeurs limites de la vitesse du véhicule; (f) sélectionner, en tant que configuration active, une première configuration de pneus; 15 (g) sélectionner, en tant que tronçon actif, un premier tronçon de trajectoire; (h) simuler, au moyen du modèle de véhicule, le parcours d'une première partie au moins du tronçon actif de trajectoire, depuis son point initial, par le véhicule 20 équipé de la configuration active de pneus, en optimisant les sollicitations des pneus pour maximiser l'accélération subie par le véhicule à partir de la vitesse limite associée au point initial de ce tronçon actif; 25 (i) mémoriser les vitesses du véhicule aux différents points de la première partie du tronçon actif et les temps de parcours entre ces différents points; (j) sélectionner, en tant que section active du tronçon actif, la portion de trajectoire définie entre un 30 premier et un deuxième points d'une paire de points du tronçon actif se succédant dans cet ordre dans le sens de 2910662 4 parcours, le deuxième point de cette paire étant constitué par le point terminal du tronçon actif; (k) sélectionner, en tant que vitesse cible, la vitesse limite associée au deuxième point de la section 5 active; (1) simuler, au moyen du modèle de véhicule, le parcours de la section active du tronçon actif par le véhicule équipé de la configuration active de pneus, en optimisant les sollicitations des pneus et la vitesse 10 adoptée par le véhicule au premier point de la section active pour atteindre le deuxième point de la section active à la vitesse cible, tout en minimisant le temps de parcours sur cette section active; (m) mémoriser le temps de parcours de la section 15 active par le véhicule; mémoriser, en tant que vitesse limite, la vitesse du véhicule au premier point de la section active; et mémoriser, en tant qu'élément d'une deuxième partie du tronçon actif, la section active parcourue; 20 (n) vérifier si la vitesse limite du véhicule au premier point de la section active dépasse ou non la vitesse, mémorisée à l'étape (i), qui est adoptée en ce même point par le véhicule lors du parcours de la première partie du tronçon actif; 25 (o) dans la négative, sélectionner, en tant que nouvelle section active du tronçon actif, une portion de trajectoire située en amont, par rapport au sens de parcours, de la section active chronologiquement précédente, et définie entre un premier et un deuxième 30 points d'une paire de points du tronçon actif se succédant dans cet ordre dans le sens de parcours, le deuxième point de cette paire étant constitué par le 2910662 5 premier point de la section active chronologiquement précédente, et répéter les étapes (k) à (n); (p) dans l'affirmative, mémoriser, en tant que temps de parcours total du tronçon actif, la somme du temps de 5 parcours total sur la deuxième partie de ce tronçon actif et des temps de parcours mémorisés sur la fraction de la première partie de ce tronçon qui ne recouvre pas la deuxième partie; (q) sélectionner, en tant que nouveau tronçon actif, 10 le tronçon de trajectoire qui succède au tronçon précédemment actif dans le sens de parcours; (r) répéter les étapes (h) à (q) jusqu'à parcourir au moins une fois la trajectoire et obtenir un temps de parcours cumulé et stabilisé de la trajectoire par le 15 véhicule. Il est possible de prévoir que le temps de référence soit obtenu en appliquant le procédé de l'invention à une autre configuration de pneus, auquel cas de procédé comprend en outre les étapes consistant à : 20 (s) sélectionner, en tant que nouvelle configuration active, une deuxième configuration de pneus; (t) répéter les étapes (g) à (r); (u) comparer les temps de parcours cumulés et stabilisés de la trajectoire, respectivement obtenus avec 25 les première et deuxième configurations de pneus; et (v) sélectionner, en tant que configuration optimisée, la configuration de pneus conduisant au plus court temps de parcours cumulé et stabilisé de la trajectoire, le temps de parcours obtenu avec chaque 30 configuration de pneus servant ainsi de temps de référence pour le temps de parcours obtenu avec l'autre configuration de pneus. 2910662 6 Dans le cas où le modèle de pneu est un modèle thermomécanique, il est possible de prévoir que l'étape (j) soit précédée d'une étape (w0) consistant à sélectionner, en tant que conditions thermiques actives, 5 des conditions thermiques prédéterminées, que l'étape (1) soit réalisée en appliquant aux pneus les conditions thermiques actives, que l'étape (m) inclue l'opération élémentaire supplémentaire qui consiste à mémoriser les sollicitations des pneus sur chaque section active, et 10 que l'étape (p) soit immédiatement précédée d'une boucle d'itération comprenant les étapes consistant à : (wi) en cas de succès de la vérification réalisée à l'étape (n), simuler au moyen du modèle de véhicule le parcours par ce véhicule équipé de la configuration 15 active de pneus, section active après section active dans le sens de parcours, de la deuxième partie du tronçon actif, en appliquant aux pneus sur chaque section active les sollicitations mémorisées à l'étape (m) pour cette section active, et en déduire de nouvelles conditions 20 thermiques d'utilisation des pneus sur la deuxième partie du tronçon actif; (w2) vérifier si l'écart entre les nouvelles conditions thermiques et les conditions thermiques actives s'inscrit ou non dans une limite de convergence; 25 et (w3) dans la négative, mémoriser, en tant que conditions thermiques actives, les nouvelles conditions thermiques et répéter les étapes (j) à (w2). L'étape (d) peut être mise en oeuvre en fixant le 30 point terminal de chaque tronçon au maximum de courbure d'un virage, et le procédé de l'invention peut comprendre une opération supplémentaire, mise en oeuvre à chaque 2910662 7 nouveau parcours de la trajectoire et consistant à déplacer le point terminal de chaque tronçon vers l'amont, en direction d'une zone de moindre glissement des pneus. 5 Le procédé de l'invention peut être optimisé en prévoyant que le modèle de pneu développé à l'étape (b) et mis en oeuvre à chaque étape de simulation (h, 1, wl) décrive le comportement physique de chaque pneu équipant le véhicule en situation de roulement sur le sol, avec 10 lequel la bande de roulement du pneu présente une aire de contact incluant au moins une zone de contact adhérent et au moins une zone de contact glissant, que l'étape (b) comprenne au moins des opérations de modélisation, que chaque étape de simulation (h, 1, wl) comprenne des 15 opérations d'affectation de valeurs numériques et des opérations de résolution, que les opérations de modélisation consistent au moins à établir, par application de lois physiques connues et / ou construites par expérimentation spécifique, et en tant que premier 20 modèle, un modèle des efforts longitudinaux, des efforts transversaux, et d'un couple d'auto-alignement, transmis par le pneu entre le sol et le véhicule, en fonction de paramètres physiques spécifiques, incluant un coefficient d'adhérence et un module de cisaillement de la gomme du 25 pneu, et en fonction de paramètres dynamiques liés aux conditions physiques de roulement et d'utilisation du pneu, le couple d'auto-alignement étant lié à l'intensité des efforts longitudinaux et transversaux et à leur distribution dans l'aire de contact, en prévoyant que les 30 opérations d'affectation consistent au moins à attribuer des valeurs numériques aux paramètres dynamiques et aux paramètres spécifiques, et que =Les opérations de 2910662 8 résolution consistent au moins à déduire, en utilisant au moins le premier modèle et les valeurs attribuées aux paramètres dynamiques et aux paramètres physiques spécifiques, des valeurs des efforts longitudinaux, des 5 efforts transversaux, et du couple d'auto-alignement, en prévoyant que les opérations de modélisation comprennent aussi l'établissement, en tant que second modèle, d'un modèle local d'échauffement exprimant des variations d'une température de contact de la bande de roulement 10 avec le sol depuis l'entrée jusqu'à la sortie de l'aire de contact par suite du contact et du glissement de la bande de roulement avec le sol, et l'établissement, en tant que troisième modèle, d'un modèle global d'échauffement et de flux thermique, ce troisième modèle 15 exprimant des variations, sur une période d'au moins un tour de roue, d'une température périphérique de la bande de roulement et d'une température interne du pneu en fonction de valeurs antérieurement connues ou estimées des températures périphérique et interne, d'un 20 coefficient de conduction thermique de la bande de roulement, et de phénomènes à composante thermodynamique tels que déformations internes subies par le pneu, échanges thermiques entre le pneu et son environnement, et glissement de la bande de roulement sur le sol, et en 25 prévoyant qu'au moins le coefficient d'adhérence et le module de cisaillement interviennent dans le premier modèle comme des variables respectivement fonctions de la température périphérique et de la température interne, que la valeur de la température périphérique, obtenue par 30 des opérations de résolution afférentes au troisième modèle, soit utilisée par des opérations de résolution afférentes au deuxième modèle pour prendre en compte la 9 2910662 dépendance en température du coefficient d'adhérence, et que la valeur de la température interne, obtenue par des opérations de résolution afférentes au troisième modèle, soit utilisée par des opérations de résolution afférentes 5 au premier modèle pour prendre en compte la dépendance en température du module de cisaillement de la gomme du pneu. Il est avantageux de prévoir que le premier modèle développé à l'étape (b) comprenne des équations liées à 10 des conditions d'équilibre des forces élémentaires de cisaillement et de glissement de la bande de roulement dans l'aire de contact, que chaque étape de simulation (h, 1, wl) comprenne au moins une phase itérative de calcul constituée d'une succession de cycles de calcul, 15 et que chaque cycle de calcul comprenne au moins une opération de résolution afférente au premier modèle et une opération de résolution afférente au deuxième modèle. Le procédé de l'invention, dans un mode préféré de réalisation, peut être mis en oeuvre en prévoyant dans 20 l'étape (b) que le premier modèle soit établi en considérant que l'aire de contact comprend une unique zone de contact adhérent et une unique zone de contact glissant séparées l'un de l'autre par un point de passage, que le premier modèle prenne la forme d'un 25 système d'équations exprimées au moins en fonction des paramètres dynamiques, des paramètres spécifiques, et de l'abscisse du point de passage, que chaque phase itérative soit dédiée à des phénomènes apparaissant pendant un intervalle de temps élémentaire correspondant, 30 et que chaque phase itérative soit mise en oeuvre pour résoudre, par approximations successives et au moins à partir de valeurs antérieurement connues ou estimées de 2910662 10 l'abscisse du point de passage, des efforts transversaux, et du couple d'auto-alignement, de nouvelles valeurs de l'abscisse du point de passage, des efforts transversaux, et du couple d'auto-alignement qui résolvent le système 5 d'équations du premier modèle pour les valeurs attribuées aux paramètres dynamiques et aux paramètres spécifiques, ce dont il résulte que les efforts longitudinaux, les efforts transversaux, et le couple d'auto-alignement du pneu peuvent être calculés en temps réel pendant le 10 roulement du véhicule. Dans ce cas, chaque nouveau cycle de calcul de chaque phase itérative comprend de préférence les opérations consistant au moins à : - calculer une nouvelle valeur provisoire de l'abscisse 15 du point de passage à partir des équations d'équilibre des forces élémentaires et des valeurs antérieurement connues ou estimées des efforts transversaux et du couple d'auto-alignement; -calculer, à partir de la nouvelle valeur provisoire de 20 l'abscisse du point de passage et des équations liant les efforts transversaux et le couple d'auto-alignement aux paramètres dynamiques, aux paramètres spécifiques, et à l'abscisse du point de passage, de nouvelles valeurs des efforts transversaux et du couple d'auto-alignement 25 utilisables pour un éventuel cycle de calcul ultérieur; - interrompre conditionnellement ladite phase itérative au moins lorsque l'écart entre la nouvelle valeur provisoire de l'abscisse du point de passage et la valeur antérieurement connue ou estimée de cette abscisse est 30 inférieur à une limite de précision prédéterminée; et - à l'interruption de ladite phase itérative, affecter aux efforts transversaux et au couple d'auto-alignement, 2910662 11 en tant que valeurs pour cette phase, les nouvelles valeurs des efforts transversaux et du couple d'auto-alignement obtenues au dernier cycle de calcul. En outre, chaque étape de simulation (h, 1, wi) du 5 véhicule intervenant dans le procédé de l'invention comprend alors avantageusement une opération mise en oeuvre après la fin de chaque phase itérative et consistant à mettre au moins à jour les paramètres dynamiques pour tenir compte des évolutions subies par 10 ces paramètres pendant le temps d'exécution de la phase itérative, et à engager une nouvelle phase itérative. Par ailleurs, les opérations de résolution afférentes au troisième modèle sont de préférence effectuées en dehors de chaque phase itérative. 15 Chaque phase itérative peut être précédée d'une phase préparatoire au cours de laquelle sont calculées des grandeurs contingentes, incluant les dimensions de l'aire de contact, en fonction des valeurs attribuées aux paramètres dynamiques et aux paramètres spécifiques, 20 chaque grandeur contingente étant exploitée dans ladite phase itérative avec une valeur qui lui a été attribuée lors de la phase préparatoire. Pour modéliser le véhicule dans son ensemble, il est possible de prévoir que chacun des premier, deuxième et 25 troisième modèles soit exploité pour chaque pneu et associé à un modèle dynamique de châssis, que le modèle de châssis fournisse aux premier, deuxième et troisième modèles, pour chaque pneu, les valeurs de certains au moins des paramètres dynamiques, et que le modèle de 30 châssis exploite, pour chaque pneu, les valeurs des efforts longitudinaux, des efforts transversaux, et du 2910662 12 couple d'auto-alignement obtenues par la mise en oeuvre des premier, deuxième et troisième modèles. D'autres caractéristiques et avantages de l'invention ressortiront clairement de la description qui 5 en est faite ci-après, à titre indicatif et nullement limitatif, en référence aux dessins annexés, dans lesquels : - la figure 1 est une vue schématique d'un exemple de trajectoire possible, constitué par le circuit de 10 compétition automobile de Jerez; - la figure 2 est un organigramme illustrant l'organisation opératoire globale du procédé de l'invention; - la figure 3 est une vue schématique agrandie d'une 15 partie du circuit illustré à la figure 1; - la figure 4 est un diagramme représentant, en fonction de la distance parcourue par le véhicule et repérée en abscisse, la vitesse adoptée par ce véhicule et repérée en ordonnées, d'une part pour la première 20 partie d'un tronçon actif, et d'autre part pour la deuxième partie de ce tronçon; - la figure 5 est une vue schématique en élévation d'une roue équipée d'un pneu auquel est appliqué un modèle de pneu utilisable dans le procédé de l'invention; 25 -la figure 6 est une vue schématique de dessus agrandie de l'aire de contact du pneu de la figure 5 avec le sol, le pneu étant en situation de dérive par rapport au vecteur vitesse (axe X); - la figure 7 est une vue schématique de face du 30 pneu de la figure 5, le pneu étant en situation de carrossage ; 13 2910662 - la figure 8 est une vue schématique partielle et en coupe radiale d'un pneu représenté dans deux situations différentes de sollicitations; - la figure 9 est une vue schématique en élévation 5 d'un pneu représenté dans deux situations différentes de sollicitations; -la figure 10 est une vue schématique de dessus d'un pneu subissant une torsion; - la figure 11 est un schéma récapitulant les 10 différentes grandeurs intervenant dans la mise en oeuvre du procédé de l'invention; -la figure 12 est une vue schématique de dessus d'un pneu subissant diverses sollicitations; - la figure 13 est une vue de dessus agrandie de 15 l'aire de contact d'un pneu subissant diverses sollicitations, et sur laquelle apparaissent la trace de points et de trajectoires caractéristiques; - la figure 14 représente, de façon tridimensionnelle, la loi liant le coefficient 20 d'adhérence à la température de contact, et à la vitesse de glissement d'un pneu, pour une pression de contact donnée; - la figure 15 représente la loi liant à la température le module de cisaillement du mélange 25 caoutchouc de la bande de roulement d'un pneu, pour un niveau donné de sollicitation en fréquence et déformation; - la figure 16 est un schéma récapitulant les phénomènes thermiques pris en compte dans le modèle 30 thermique global du pneu; 14 2910662 - la figure 17 est un organigramme illustrant la structure opérationnelle du modèle de pneu préféré utilisable dans le procédé de l'invention; - les figures 18a et 18b sont des organigrammes 5 illustrant la structure opérationnelle d'un autre modèle de pneu préféré utilisable dans le procédé de l'invention; et - la figure 19 est un schéma illustrant le fonctionnement d'un véhicule bicycle en régime permanent. 10 Comme annoncé précédemment, l'invention concerne un procédé permettant de sélectionner, parmi un ensemble de configurations possibles, une configuration optimisée de pneumatiques pour équiper un véhicule. L'invention s'applique en particulier à l'équipement 15 en pneus d'un véhicule de compétition destiné à courir sur un circuit tel qu'illustré à la figure 1 en suivant une trajectoire prédéterminée, cette trajectoire étant orientée dans un sens unique de parcours SU et incluant des virages. 20 Ce procédé requiert le développement d'un modèle physique du véhicule équipé d'une configuration de pneus, impliquant une modélisation de ces pneumatiques et permettant une simulation du véhicule équipé de ces pneus. 25 Une fois disponibles les modèles des pneumatiques et du véhicule, le procédé de l'invention réalise une simulation du parcours de la trajectoire par le véhicule modélisé, dans laquelle la vitesse du véhicule est optimisée. 30 Le temps de parcours mis par le véhicule pour parcourir au moins une fois la trajectoire est alors mémorisé et comparé à un temps de référence. 2910662 15 Il est ainsi possible, en fonction du résultat de cette comparaison, d'évaluer par rapport à une référence la configuration de pneus utilisée dans la simulation, et donc de classer cette configuration comme étant ou non 5 optimisée par rapport à la référence. Cette référence peut elle-même être constituée par une configuration de pneus équipant un véhicule simulé, de sorte que deux configurations différentes de pneus peuvent servir de référence l'une par rapport à l'autre. 10 Dans le mode de réalisation préféré de l'invention, le modèle utilisé pour les pneus est un modèle thermomécanique, la simulation pouvant ainsi atteindre un haut degré de réalisme. Un exemple détaillé de mise en oeuvre du procédé de 15 l'invention est donné ci-après dans le cas le plus complet, dans lequel la référence est constituée par une configuration de pneus ayant fait l'objet d'une simulation, et dans lequel le modèle de pneus est un modèle thermomécanique. 20 Un tel procédé comprend typiquement les étapes suivantes (figure 2). (a) L'étape (a) consiste à recenser l'ensemble des paramètres physiques ayant une influence, notamment dans le sens d'une réduction, sur la vitesse que peut adopter 25 le véhicule sur la trajectoire. Cet ensemble de paramètres comprend des paramètres liés aux différents pneus possibles, au châssis du véhicule et/ou aux conditions de roulement, et inclut des paramètres constants et des paramètres dynamiques. En particulier, 30 le sous-ensemble des paramètres constants liés aux conditions de roulement comprend une image numérique de la trajectoire, constituée par des coordonnées numérisées 2910662 16 de points successifs de cette trajectoire, tels que les points P11 ou P23 visibles sur la figure 3. (b) L'étape (b) consiste à développer un modèle de pneu décrivant le comportement physique de chaque pneu en 5 fonction des paramètres physiques pertinents, et notamment en fonction de la température. Un contenu possible de l'étape (b) est décrit en détail dans la section intitulée "Modélisation et simulation des pneumatiques". 10 (c) L'étape (c) consiste à développer un modèle de véhicule équipé de pneumatiques. Un contenu possible de l'étape (c) est décrit en détail dans la section intitulée "Modélisation et simulation du véhicule". (d) L'étape (d) consiste à découper l'image 15 numérique de la trajectoire en tronçons adjacents tels que TRi, TR2, etc. (figures 1 et 3), ces tronçons se succédant dans le sens de parcours SU de la trajectoire du véhicule. Chaque tronçon comprend plusieurs points de trajectoire régulièrement espacés les uns des autres, par 20 exemple avec un pas de 5 mètres. Chaque tronçon, supposé orienté dans le sens de parcours SU, commence en un point initial, tel que P11 pour TRI et P21 pour TR2, et se termine par un point terminal tel que P17 pour TR1 et P25 pour TR2, le point terminal et le point initial de deux 25 tronçons adjacents, par exemple P17 et P21 des tronçons TRI et TR2, étant confondus. Chaque tronçon est choisi de manière à présenter au plus une sortie de virage et une entrée d'un autre virage. Le point terminal de chaque tronçon est ainsi localisé dans une entrée de virage, et 30 par exemple choisi, par défaut, au point de courbure maximal ou à mi-parcours du virage. 2910662 17 (e) L'étape (e) consiste à associer respectivement au point initial et au point terminal de chaque tronçon des valeurs limites de la vitesse du véhicule, telles que VL11 et VL12. La vitesse limite VL12 au point terminal 5 P17 d'un tronçon TR1 et la vitesse limite VL12 au point initial P21 du tronçon suivant TR2 sont évidemment fixées à la même valeur dans la mesure où ces points sont confondus. Ces valeurs limites sont par exemple choisies ou estimées sur la base de vitesses concrètement mesurées 10 en compétition. (f) L'étape (f) consiste à sélectionner, en tant que configuration active, une première configuration de pneus destinée à être testée. (g) L'étape (g) consiste à sélectionner, en tant que 15 tronçon actif, un premier tronçon de trajectoire pour commencer la simulation, par exemple le tronçon TRI (figures 1 et 3). (h) L'étape (h) consiste à simuler, au moyen du modèle de véhicule développé à l'étape (c), le parcours 20 du tronçon actif TRI ou d'une première partie au moins Pl TA de ce tronçon actif (figure 4), depuis son point initial P11, par le véhicule équipé de la configuration active de pneus. En fait, cette simulation est réalisée en optimisant les sollicitations des pneus de manière que 25 le véhicule accélère au maximum à partir de la vitesse limite Vil associée au point initial P11 de ce tronçon actif, en l'occurrence du tronçon TRI. Dans la mesure où elle fait intervenir des relations non linéaires, cette optimisation doit être réalisée au moyen d'un outil 30 approprié, tel que l'algorithme SQP de programmation quadratique séquentielle, qui est par exemple disponible dans les outils du progiciel vendu sous la dénomination 2910662 18 "Matlab 7". Bien que la simulation réalisée à l'étape (h) puisse être effectuée sur toute la longueur du tronçon actif, en l'occurrence TRI, elle n'a pas besoin de l'être, ce qui justifie l'indication du fait que cette 5 simulation est réalisée "sur une première partie au moins P1 TA du tronçon actif". En effet, la simulation réalisée à l'étape (h) est conçue pour simuler le déplacement du véhicule en pleine accélération. Or, dans la mesure où chaque tronçon - à l'exception éventuelle du dernier - se 10 termine par un virage, le véhicule estcontraint d'abandonner son régime d'accélération maximale avant la fin de chaque tronçon. (i) L'étape (i), qui est mise en oeuvre au fur et à mesure du déroulement de l'étape (h), consiste à 15 mémoriser les vitesses adoptées par le véhicule aux différents points, tels que P11, P12, P13, etc. de la première partie P1_TA du tronçon actif, en l'occurrence TRi, et les temps de parcours du véhicule entre ces différents points. Le résultat obtenu à cette étape est 20 équivalent à celui qu'illustre le diagramme de la figure 4, qui représente la vitesse du véhicule au passage des différents points du tronçon, les temps de parcours étant définis par les rapports de la distance par la vitesse. (w0) L'étape (w0) consiste à sélectionner, en tant 25 que conditions thermiques actives, des conditions thermiques prédéterminées. Cette étape, qui sert à préparer une boucle d'itération permettant le calcul par approximations successives des conditions de freinage du véhicule à la fin du tronçon actif, peut être mise en 30 oeuvre en utilisant par défaut les températures de pneus préalablement mesurées ou estimées sur des véhicules réels aux différents points du tronçon actif. 19 2910662 (j) L'étape (j) consiste à sélectionner, en tant que section active du tronçon actif, la portion de trajectoire définie entre l'avant-dernier point et le point terminal du tronçon actif, c'est-à-dire la section 5 entre les points P16 et P17 lorsque le tronçon actif est TRI (figure 3). (k) L'étape (k) consiste à sélectionner, en tant que vitesse cible, la vitesse limite associée au deuxième point de la section active, c'est-à-dire la vitesse VL12 10 au point P17 lorsque la section active est celle qui s'étend entre P16 et P17. (1) L'étape (1) consiste à simuler, au moyen du modèle de véhicule, le parcours de la section active (en l'occurrence de P16 à P17) du tronçon actif (en 15 l'occurrence TRI) par le véhicule équipé de la configuration active de pneus. Cette simulation est réalisée en appliquant aux pneus les conditions thermiques actives et en optimisant les sollicitations des pneus et la vitesse V:L6 du véhicule au premier point 20 P16 de la section active pour atteindre le deuxième point P17 de la section active à la vitesse cible VL12, tout en minimisant le temps de parcours sur cette section active. De même que l'étape (h), l'étape (1) fait donc intervenir des outils d'optimisation avancés quoique bien connus de 25 l'homme du métier et commercialement disponibles. En revanche, contrairement à l'étape (h), l'étape (1) est destinée à décrire le déplacement du véhicule en phase de décélération. (m) l'étape (m) consiste à effectuer des stockages 30 de données dans une mémoire informatique, et notamment à mémoriser le temps de parcours de la section active (en l'occurrence P16 à P17) par le véhicule, à mémoriser, en 2910662 20 tant que vitesse limite, la vitesse V16 du véhicule au premier point P16 de la section active, à mémoriser, en tant qu'élément d'une deuxième partie P2_TA du tronçon actif (en l'occurrence TRI), la section active P16 à P17 5 qui vient d'être parcourue, et à mémoriser les sollicitations des pneus sur cette section active. La référence ici faite à une "deuxième partie P2_TA" du tronçon actif est justifiée par le souci de distinguer cette deuxième partie, sur laquelle le véhicule est 10 simulé en décélération, de la première partie Pl_TA sur laquelle le véhicule a été simulé en accélération. Par ailleurs, comme le montre la figure 4 de façon symbolique par des flèches, la première partie P1_TA est simulée en progressant depuis le point initial du tronçon actif, 15 alors que la deuxième partie P2 _TA est simulée en progressant depuis le point terminal du tronçon actif. (n) L'étape (n) consiste à vérifier si la vitesse limite V16 du véhicule au premier point P16 de la section active dépasse ou non la vitesse, mémorisée à l'étape 20 (i), qui est adoptée en ce même point P16 par le véhicule lors du parcours de la première partie Pl TA du tronçon actif. De façon plus intuitive, cette étape consiste à tester si la deuxième partie P2_TA, telle qu'elle a été parcourue jusqu'à présent, est ou non suffisante pour se 25 raccorder à la première partie Pl TA dont les caractéristiques ont été mémorisées à l'étape (i). (o) L'étape (o), qui est seulement mise en oeuvre en cas de réponse négative au test de l'étape (n), consiste d'abord à sélectionner, en tant que nouvelle section 30 active du tronçon actif, la portion de trajectoire définie entre deux points consécutifs du tronçon actif et située immédiatement en amont de la section active 2910662 21 chronologiquement précédente, c'est-à-dire en l'occurrence la section s'étendant de P15 à P16, puis à répéter les étapes (k) à (n) en les appliquant à la nouvelle section active P15 à P16. 5 (wl) L'étape (wl), qui est seulement mise en oeuvre en cas de réponse affirmative au test de l'étape (n), consiste à simuler, au moyen du modèle de véhicule, le parcours par ce véhicule équipé de la configuration active de pneus, section active après section active dans 10 le sens de parcours SU, de la deuxième partie P2_TA du tronçon actif, en appliquant aux pneus sur chaque section active les sollicitations mémorisées à l'étape (m) pour cette section active, et en déduire de nouvelles conditions thermiques d'utilisation des pneus sur cette 15 deuxième partie P2_TA du tronçon actif. De façon plus intuitive, l'étape (wl) consiste donc à parcourir dans le sens réel SU la partie P2_TA du tronçon actif sur laquelle le véhicule est en décélération, en rejouant les sollicitations préalablement mémorisées des pneus pour en 20 déduire des conditions thermiques de fonctionnement des pneus plus adaptées au cas particulier que les conditions thermiques choisies par défaut à l'étape (w0) et précédemment utilisées à l'étape (1) en tant que conditions thermiques actives initiales. 25 (w2) L'étape (w2) consiste à vérifier si l'écart entre les nouvelles conditions thermiques obtenues à l'étape (wl) et les conditions thermiques actives s'inscrit ou non dans une limite de convergence. En d'autres termes, l'étape (w2) teste si les nouvelles 30 conditions thermiques obtenues sont acceptables, ou s'il est nécessaire de les affiner encore davantage. 22 2910662 (w3) L'étape (w3), qui est seulement mise en oeuvre si les nouvelles conditions thermiques ne sont pas encore acceptables, consiste à mémoriser, en tant que conditions thermiques actives, les nouvelles conditions thermiques 5 obtenues à l'étape (wi) et à répéter les étapes (j) à (w2) ; (p) L'étape (p), qui est seulement mise en oeuvre si les nouvelles conditions thermiques ont suffisamment convergé et sont donc acceptables, consiste à mémoriser, 10 en tant que temps de parcours total du tronçon actif (en l'occurrence toujours TR1), la somme du temps de parcours total sur la deuxième partie P2 TA de ce tronçon actif et des temps de parcours mémorisés sur la fraction de la première partie Pi TA de ce tronçon qui ne recouvre pas 15 la deuxième partie, c'est--à-dire sur la fraction de P1_TA située à gauche du point de rencontre de P2_TA et Pl TA sur la figure 4. (q) L'étape (q) consiste à sélectionner, en tant que nouveau tronçon actif, le tronçon de trajectoire, en 20 l'occurrence TR2, qui succède au tronçon précédemment actif TRI dans le sens de parcours SU. (r) L'étape (r) consiste à répéter les étapes (h) à (q) jusqu'à parcourir au moins une fois la trajectoire et à obtenir un temps de parcours cumulé et stabilisé de la 25 trajectoire par le véhicule. (s) L'étape (s) consiste à sélectionner, en tant que nouvelle configuration active, une deuxième configuration de pneus. (t) L'étape (t) consiste à répéter les étapes (g) à 30 (r) en utilisant la nouvelle configuration de pneus. (u) L'étape (u) consiste à comparer les temps de parcours cumulés et stabilisés de la trajectoire, 23 2910662 respectivement obtenus avec les première et deuxième configurations de pneus. (v) Et l'étape (u) consiste à sélectionner, en tant que configuration optimisée, la configuration de pneus 5 conduisant au temps de parcours cumulé et stabilisé de la trajectoire le plus court, le temps de parcours obtenu avec chaque configuration de pneus servant ainsi de temps de référence pour le temps de parcours obtenu avec l'autre configuration de pneus. 10 Comme indiqué précédemment, la position du point terminal de chaque tronçon peut être choisie à l'étape (d) de façon conventionnelle, ce point étant par exemple fixé au maximum de courbure du virage final de ce tronçon. 15 Néanmoins, cette position conventionnelle est rarement optimale, les pilotes ayant tendance à remettre le véhicule en accélération en virage avant le point de courbure maximal. Pour optimiser la simulation du véhicule, le 20 procédé de l'invention peut comprendre une opération supplémentaire, qui est mise en oeuvre à chaque nouveau parcours de la trajectoire, et qui consiste à déplacer le point terminal de chaque tronçon vers l'amont, en direction d'une zone de moindre glissement des pneus. 25 MODELISATION ET SIMULATION DES PNEUMATIQUES La présente section décrit un mode possible de 30 réalisation de l'étape (b) et de chaque étape du procédé de l'invention qui implique la simulation du comportement 24 2910662 physique de chaque pneu équipant le véhicule en situation de roulement sur le sol. Dans une situation de roulement d'un pneu au sol, la bande de roulement du pneu présente avec le sol une 5 aire de contact formée d'au moins une zone de contact adhérent et d'au moins une zone de contact glissant. La finalité de l'étape (b) est de développer les outils permettant de fournir, en fonction d'un certain nombre de grandeurs d'influence, les valeurs de divers 10 efforts qui sont transmis par chaque pneu entre le sol et le véhicule. La présente section de description fera usage des conventions d'écriture ci-après, qui seront spécifiques à la modélisation des pneumatiques, et. dont une partie est 15 déjà communément utilisée dans des modèles de pneu existants connus de l'homme du métier. 6 Angle de dérive b' Angle de dérive au niveau des nappes sommet a~ Angle de mise à plat (plysteer) Y Angle d'inclinaison (carrossage) r Taux de glissement longitudinal Fx Effort longitudinal Fy Effort transversal Fz Charge de la roue V Vitesse des points du sol WX Vitesse des points du sommet suivant son axe Mz Couple d'auto-alignement Ny Composante du couple d'auto-alignement liée à l'effort Fy Nx Composante du couple d'auto-alignement liée à la distribution des efforts Fx dans la largeur de l'aire de contact 2910662 25 RL Rigidité latérale du pneu RLL Rigidité longitudinale du pneu kT Rigidité de torsion du pneu S2 Souplesse à la flexion sur chant du bloc sommet p Coefficient d'adhérence dynamique entre la gomme du pneu et le sol : p(p,Vg,Ts) Po Coefficient d'adhérence statique Lx Longueur caractéristique de l'aire de contact Ly Largeur de l'aire de contact ent Taux d'entaillement de la bande de roulement AssX Coefficient d'assouplissement longitudinal de la bande de roulement AssY Coefficient d'assouplissement latéral de la bande de roulement G* Module de cisaillement du mélange caoutchouc de la bande de roulement du pneu eue, Epaisseur de la bande de roulement a Demi-longueur de l'aire de contact b Abscisse de début de glissement dans l'aire de contact Ti Profil de température interne de la gomme, entre la surface du pneu et les nappes sommet Ts Température de surface moyenne de la bande de roulement Tc Température de contact à l'interface entre la gomme et le sol Vg Vitesse de glissement entre le pneu et le sol XN, Coordonnées des points du sommet du pneu (zone interne YN de la bande de roulement, à la base du pneu, à la verticale de l'aire de contact) Xk,Yk Coordonnées des points de la bande de roulement au niveau de l'interface avec le sol Les figures 5 à 7 précisent le repère utilisé comme référence. Ce repère est défini par : 5 0 : le point origine du repère au centre de l'aire de contact; 2910662 26 OX : l'axe parallèle au vecteur vitesse; OY : l'axe perpendiculaire à OX, parallèle au plan du sol quel que soit le carrossage; Dans ce repère, la convention de signe impose que : 5 pour z>0, un effort longitudinal est généré dans la direction de l'axe OX; pour 5>0, un effort latéral est généré dans la direction de l'axe OY, et pour y>0, un couple dit d'auto-alignement Mz négatif est généré, provoquant une poussée latérale négative (c'est-à-dire 10 orientée dans la direction opposée à OY). Les efforts transmis par le pneu entre le sol et le véhicule comprennent des efforts longitudinaux Fx, des efforts transversaux Fy, et un couple d'auto-alignement Mz, qui est lié à l'intensité des efforts longitudinaux 15 et transversaux et à leur distribution dans l'aire de contact. Les grandeurs d'influence comprennent typiquement d'une part des paramètres dynamiques, c'est-à-dire variables au moins en fonction du temps, et liés aux 20 conditions physiques de roulement et d'utilisation du pneu, et d'autre part des paramètres physiques qui sont des paramètres spécifiques au pneu considéré et de valeur constante dans des conditions de référence déterminées. Les paramètres dynamiques comprennent l'angle de 25 dérive, le taux de glissement, l'angle de carrossage, la charge, la vitesse, la pression de gonflage, les températures air et sol, les températures initiales du pneumatique, et le temps. Les paramètres spécifiques comprennent les 30 dimensions de l'aire de contact (longueur, largeur, coefficient de forme), le profil de pression p(x) le long de l'aire de contact, les rigidités longitudinale Kx et 2910662 27 transversale Ky de la bande de roulement, les rigidités de la structure du pneu, à savoir la rigidité latérale RL, la rigidité longitudinale RLL, la rigidité radiale Rrr, et la rigidité de torsion kT, la rigidité l/S2 du 5 bloc sommet, une loi d'adhérence ,u du couple gomme/sol, des paramètres de transfert de longueur entre reliefs du pneu ("ribs"), et sont reliées aux paramètres dynamiques par des relations construites par expérimentation spécifique dont il est possible de rendre compte par les 10 expressions données ci-après à titre d'exemple. La rigidité latérale correspond à un déport dy de l'aire de contact (cf. figure 8) par rapport au plan de roue sous l'effet d'un effort latéral : 15 RL == RLO + RLp p où RLO [N/m] représente la partie structurelle et RLP [N/m/bars] la partie pneumatique, p étant la pression exprimée en bars. 20 La rigidité longitudinale correspond à un déport dx de l'aire de contact (cf. figure 9) suivant l'axe longitudinal de la roue en présence d'un effort longitudinal Fx : 25 RLL _: RLLO + RLLP p où RLLO [N/m] représente la partie structurelle et RLLP [N/m/bar] la partie pneumatique, p étant la pression exprimée en bars. 5 2910662 28 Le développement du couple d'auto-alignement Mz entraîne une torsion d'un angle 08 de l'enveloppe autour de l'axe Z par rapport au plan de jante (cf. figure 10) : X18=Mz kr La rigidité de torsion de l'enveloppe contient une composante structurelle kTo [N.m/rad] et une composante kTz [m/rad] qui traduit la variation de la rigidité de 10 torsion avec la charge par exemple de la forme : k r = (kro + krz Fz) 11 f' L'angle de dérive réel 8' pris par l'aire de 15 contact s'exprime en fonction de l'angle de dérive au niveau de l'axe de la roue de la manière suivante : 8'=B+Mz kr 20 La déformation du sommet peut: être modélisée par une loi du second degré avec une courbure au centre de l'aire de contact qui s'exprime par exemple de la manière suivante : p=S2Fy 25 où S2 est un paramètre représentatif de la souplesse à la flexion sur chant. La rigidité radiale relie la charge Fz à la flèche du sommet par rapport à la jante. Elle dépend de la pression et se décompose en deux termes : un terme 2910662 29 structurel RRO [N/m] qui correspond à la rigidité radiale du pneumatique à pression nulle, et un terme pneumatique RRP [N/m/bar] . RR = RRO + RRp p 5 La longueur de l'aire de contact est définie par : Lx=a,l Pc +bFZ cette formulation permettant de prendre en compte l'effet de la charge et de la pression de gonflage. 10 La largeur de l'aire de contact est définie par : Ly=Ly,+2*Lye où Ly, est la largeur des reliefs au centre du pneu et 15 Ly, la largeur des reliefs aux épaules, calculée par la formule ci-dessus : _ i F \z Lye = c arctan d z ù e La surface effective de l'aire de contact est définie comme le produit de la largeur par la longueur, 20 pondéré par l'entaillement et un coefficient de forme : SADC = C forme ent Lx Ly] Le coefficient de forme Cforme prend en compte la 25 variation de forme de l'aire de contact en fonction de la charge. 30 2910662 Lors de l'écrasement de la structure du pneu, le sommet prend une inclinaison a, qui est une grandeur caractéristique de la structure étudiée. Le profil de pression entre l'entrée et la sortie 5 de l'aire de contact est déterminé de la manière suivante . 2n+1 FZ x 1 zn P(x)- 2n SADC 1-(Lx/2 ~ Cette répartition de pression est telle que : Lx/2 Fz Jp(x)dx= -Lx/2 SADC A faible charge, le profil est plutôt de type parabolique (n = 1). A forte charge, le profil de pression est quasi uniforme. De préférence, n est un réel variant linéairement 15 avec la longueur de l'aire de contact Lx. Pour éviter qu'aux faibles charges n ne devienne trop faible (voire négatif), n est borné inférieurement à 1 de la façon suivante : n = max(1, na Lx + nb) . 20 Le modèle de pneu mis en oeuvre dans l'invention utilise de préférence au moins un modèle mécanique (premier modèle) couplé à un modèle thermique local (deuxième modèle) et à un modèle thermique global 25 (troisième modèle). Chacun de ces modèles est établi au cours d'une phase préliminaire de la modélisation du pneu et prend la forme d'un système d'équations. 10 2910662 31 Ces modèles peuvent être établis de plusieurs manières, et notamment en utilisant des lois physiques connues de l'homme du métier, ou des lois plus ou moins approchées construites par expérimentation spécifique, de 5 sorte que les équations par lesquelles s'expriment ces modèles peuvent prendre de multiples formes. Ces modèles se caractérisent donc essentiellement par leurs grandeurs d'entrée, leurs grandeurs de sortie, et le fait qu'ils prennent chacun la forme d'un système 10 d'équations rendant compte de phénomènes physiques observables en fonction de propriétés physiques quantifiables. Les grandeurs de sortie du modèle mécanique de pneu sont les efforts longitudinaux Fx, les efforts latéraux 15 Fy, et le couple d'auto-alignement Mz. La grandeur de sortie du modèle thermique local est la température de la bande de roulement le long de l'aire de contact. Et les grandeurs de sortie du modèle thermique 20 global du pneu sont la température moyenne périphérique ou de surface Ts de la bande de roulement, et le profil de la température interne Ti de cette bande dans le sens de son épaisseur. Les grandeurs d'entrée et de sortie de l'ensemble 25 de ces modèles sont recensées sur la figure 11. Le modèle mécanique de pneu préféré se caractérise en outre par le fait qu'il est établi en considérant d'une part que l'aire de contact comprend deux zones (figure 12), à savoir une zone de contact adhérent où 30 l'effort est contrôlé par le cisaillement de la bande de roulement à déplacement imposé et une zone de contact glissant où l'effort est contrôlé par le coefficient de 32 2910662 frottement entre la gomme et le sol, et en considérant d'autre part qu'il existe un point unique N d'abscisse "b" marquant le passage entre la zone de contact adhérent et la zone de contact glissant. 5 Les équations sont formulées selon ce principe, ce qui permet d'obtenir des expressions rapides à résoudre. Un exemple de modèle mécanique avantageux est donné ci-après. Dans cet exemple, la modélisation du fonctionnement 10 de l'aire de contact est basée sur une approche de type poil de brosse avec une première phase de cisaillement en entrée de l'aire de contact et une seconde phase de glissement. Il est supposé que ces deux phases sont distinctes, uniques et connexes, et qu'il n'y 15 a aucun mécanisme parasite de reprise de cisaillement lors de la phase de glissement. Tous les développements qui suivent font l'hypothèse que les angles de dérive restent modérés (inférieurs à 20 degrés environ), de sorte que 20 l'approximation tan(S) S est valide et sera faite systématiquement. Kx et Ky désignent les rigidités de la bande de roulement du pneu, qui peuvent être reliés au module de la gomme et aux caractéristiques de la sculpture de 25 préférence selon la relation suivante : KX = h scre AssX +h sc G * ent K _ ,, - hscre +h AssY s` G*ent 2910662 33 hscre étant l'épaisseur de la sculpture, et hsc l'épaisseur de la sous-couche, telles que eK.M = hSCYe +hsc. La figure 13 est un schéma fonctionnel de l'aire de contact. Le segment NK définit un élément de la bande de 5 roulement ("poil de brosse"). N est le point situé au niveau du sommet et K est le point de la bande de roulement au niveau du sol. Le point d'abscisse b représente la transition entre le contact avec adhérence et le contact glissant. 10 En entrée de l'aire de contact (voir figure 13), l'élément de gomme de la bande de roulement du pneu n'est pas cisaillé (XN = XK) . Le cisaillement de la gomme a en fait plusieurs origines : une mise en dérive de la roue d'un angle 8, 15 une mise en carrossage d'un angle y, et une différence entre la vitesse du point N du sommet et la vitesse de défilement des points du sol. Sous l'hypothèse d'une déformation uniforme dans l'épaisseur de la bande de roulement, l'effort 20 élémentaire généré par le cisaillement d'un l'élément de cette bande peut s'écrire : dFx =Kx(XK -XN)dS 25 dFy = Ky(YK -YN)dS où dS est la surface élémentaire de l'élément NK de la bande de roulement. L'équation de la trajectoire des points du sommet 30 est approximée par la relation suivante : 15 34 2910662 1 2 FY YN=BXN--2S2FYXN +R L expression dans laquelle 8' est l'angle de dérive du sommet, qui diffère de l'angle de dérive du fait de la 5 torsion de la structure du pneu, et qui répond à l'équation : 8'= 8+al+ Mz kr En supposant que la relation YK(a)=YN(a) soit satisfaite 10 en entrée de l'aire de contact, il vient (équation 1) . %K -YN =8'(a-.XN)- 2S2Fy(a2 -XN2) En posant : V la vitesse des points du sol, Wx la vitesse des points du sommet suivant son axe, et WX V = z V 20 l'expression XK - XN devient (équation 2) . -1- 1 K -XN =(a-XN)+z Par définition, r correspond au taux de glissement 25 longitudinal. Les composantes de la vitesse de glissement sont données par : 2910662 Vgy = Dans la partie glissante de l'aire de contact, les 5 efforts élémentaires sont générés par le frottement entre la gomme et le sol et la direction des efforts est colinéaire au vecteur cisaillement, ce qui se traduit par : JVgx d(X dtXN)=Wx V d(YK -YN) = (S'-S2FYXN)Wx dt dFx =Y (XK -XN) dFY K (YK ù YN ) avec r, = Kx KY 10 20 En notantf=S'-~S2FY(a+XN), il vient : dFx z dFY -rK (1+),6 15 Les efforts élémentaires dans la zone de frottement s'écrivent : dFx = rKr ,u(p,Tg,T)pdS j(rKr)2 +(1+2)2)32 (1 + T)(3 dFY = ,u(p, Vg, T )p dS \/(rKi)2 +(l+2)2/32 L'abscisse b du point N à partir duquel est supposé se produire le glissement du pneu sur le sol correspond à l'équilibre entre les efforts élémentaires de 36 2910662 cisaillement et d'adhérence, cet équilibre étant représenté par l'équation 3 : [Kx(XK -X1 )r [KY(YK - TN)2 .. =[u0pz où po est le coefficient d'adhérence statique exprimé au point d'abscisse b. Bien qu'il puisse a priori exister dans l'aire de 10 contact plusieurs points de transition entre une zone d'adhérence et une zone glissante, le modèle mécanique de pneu utilisé dans le mode de réalisation préféré de l'invention fait avantageusement l'hypothèse de l'unicité de ce point de transition. En d'autres termes, dès 15 qu'apparaît un glissement dans l'aire de contact, ce glissement est supposé perdurer jusqu'à la sortie de cette aire de contact. Les équations représentatives du bilan des efforts sous réserve de cette hypothèse d'unicité du point de 20 passage sont données ci-après. Il est toutefois possible d'en donner une version plus générale, correspondant au cas où plusieurs points de passage existeraient dans l'aire de contact. En outre, le mode de réalisation illustré aux 25 figures 18a et 18b, qui utilise une discrétisation de la bande de roulement dans l'aire de contact et qui sera discuté ultérieurement, ne fait aucune hypothèse a priori sur le nombre de points de passage dans l'aire de contact. 30 Les efforts appliqués au centre de la roue équipée du pneu considéré sont obtenus en intégrant les efforts 5 5 10 37 2910662 élémentaires se développant à la surface de l'aire de contact : Fx = f f Kx(XK ùXN)dS+ f f 2 rKr 2 2 p(p,Vg,T)pdS -Ly/2b -Ly/2-a\J(rKZ) +(1+Z) fi Ly/2 a Ly/2 b FY = f f Ky(YK ùYN)dS+ f (ù (1+z)/j ,u(p,Vg,T)pdS -Ly/2 b -Ly/2-a~rKi) +(1+i) z 2 L'intégration conduit respectivement aux équations 4 et 5 suivantes : b Fx =KXLy[I z (aùb)2]+Lyent f- 2 rKr 2 ,u(p,Vg,T)pdx 21+r -auJ(r,r) +(1+z) J3 b Fy =KyLy~8'(aùb)2 -526 Y-(2a3 -3a2b+b3) +Lyent f 1+i)~ y(p,Vg,T)pdx _ -a-\j(rKZ)2 +(1+)2/62 Le couple d'auto-alignement Mz a deux contributions, à savoir un moment NY lié à l'effort Fy dont le centre de poussée est décalé par rapport au 15 centre de l'aire de contact, et un couple Nx lié à la distribution des efforts Fx dans la largeur de l'aire de contact. En général, le couple Nx est un couple de rappel, sauf dans des cas particuliers de fort couple moteur. 20 Dans le même cadre d'hypothèses que précédemment, le moment NY peut être calculé directement par l'équation Ly/2 a Ly/2 b 6: 38 2910662 6 +LyentJ 1+r ,u(p, Vg,T )pxdx -a,j(rKr)2 +(1+r)2p2 Le couple Nx est engendré par une distribution non uniforme des efforts Fx dans la largeur de l'aire de 5 contact, qui tend à s'amplifier lors de la mise en trapèze de l'aire de contact sous l'effet du déport ou du carrossage. Dans une approche de modélisation à une seule bande de relief du pneu, la distribution des efforts FX dans la largeur de l'aire de contact n'est pas 10 directement accessible. Ainsi, le couple Nx est modélisé par une relation ad hoc dont une formulation mathématique est donnée ci-après à titre indicatif en fonction du moment Ny et du carrossage (équation 7) . 15 Nx = fil (Fz)Ny + Q2 (Fz, v)y Un exemple de modèles thermiques local et global est donné ci-après. Le modèle thermique local prend essentiellement en 20 compte les phénomènes thermiques liés au contact de la bande de roulement du pneu avec le sol dans l'aire de contact et au glissement relatif de la bande de roulement dans une partie de l'aire de contact. Le modèle thermique global prend en compte 25 l'ensemble des phénomènes d'échauffement et des échanges thermiques du pneu sur au moins un tour de roue. La formulation des efforts globaux du pneumatique est basée sur une décomposition en effort de cisaillement de la bande de roulement et en effort de frottement. 30 L'effort de frottement est fonction du coefficient d'adhérence p entre la gomme et le sol, ce coefficient NY = KY Ly 6 8'(a3 -3ab2 +2b3)- S2F Y (a2 ùb2)2 2910662 39 dépendant de la pression, de la vitesse de glissement et de la température de contact. La température de contact dans l'aire de contact 5 est modélisée de la façon suivante (modèle thermique local). Lors du passage de la gomme dans l'aire de contact, la température de contact évolue en fonction de la conduction thermique et du frottement entre la gomme et 10 le sol. La température dans le contact peut être calculée de différentes manières connues de l'homme de l'art, par exemple à l'aide de méthodes de discrétisation par différences finies. Le procédé décrit ci-après optimise le temps de 15 calcul tout en conduisant à uneprécision largement suffisante. En faisant l'hypothèse de deux matériaux semiinfinis avec une température homogène (Ts pour la gomme et Tsol pour le sol), la température de surface, lorsque 20 les deux massifs sont brusquement mis en contact parfait, s'écrit . 'se g +T o = sor esor T = e +e g s01 25 où eg et esol sont respectivement l'effusivité de la gomme et du sol. Lorsqu'il y a glissement entre la gomme et le sol, le flux de frottement ÇPF génère une élévation de 30 température de surface, qui lorsque le flux est constant, s'exprime par : 2910662 Tc(t)=T~+2a9, t , eg 2r où a est le coefficient de partage qui détermine la proportion de flux qui pénètre dans la gomme. Dans le cas 5 d'un contact parfait, ce coefficient s'exprime par : a = eg eg + esor Pour la loi d'adhérence gomme sol, dont la figure 14 10 illustre la dépendance vis-à-vis des paramètres Vg, et Tc, plusieurs formulations mathématiques évidentes pour l'homme de l'art permettent de reproduire l'évolution avec la température, vitesse et pression à partir d'expérimentations spécifiques. 15 A titre d'exemple, on peut utiliser la formule : i Tc \\- P1 +(f[2 ù,ul)exp ùa2log2 avec T = To + a21og10 (Vg / Vo)/[al -- loglo (Vg / Vo )] , où / 11 , f12T0, a, al, a2 i el, e2, e3, V0 sont des constantes du modèle 20 Comme le montre la figure 14, le coefficient d'adhérence f.[ a une évolution complexe avec la température et la vitesse de glissement : aux faibles températures, ce coefficient augmente avec la température; aux fortes températures, c'est l'inverse. Le 25 coefficient u passe donc par un maximum avec la température. Ce maximum est différent selon la valeur de la vitesse de glissement. Plus la vitesse de glissement ,u(Tc,Vg,P~) ù (ei.PceZ +e3 41 2910662 est élevée, plus le maximum de ce coefficient est obtenu à des températures élevées. Le modèle thermique global calcule le profil de température dans l'épaisseur de la gomme en moyenne par 5 tour de roue et dans la largeur de la bande de roulement. Ce modèle permet d'obtenir la température Ti de l'intérieur de la bande de roulement, qui détermine la rigidité G*(Ti), ainsi que la température de surface (ou périphérique) Ts de la bande de roulement en entrée 10 d'aire de contact, qui est utilisée pour le calcul thermique dans l'aire de contact (modèle local). Un exemple de loi liant la rigidité à la température est illustré à la figure 15. En fait, cette loi est propre à chaque matériau utilisé et dépend de la 15 formulation des mélanges constitutifs de la gomme du pneu. D'une manière générale, lorsque la température du mélange augmente, sa rigidité diminue. Le modèle thermique global prend en compte les 20 mécanismes suivants : -la conduction dans la gomme; - l'échauffement par le frottement entre la gomme et le sol; - l'échauffement lié aux pertes dans la gomme; et 25 - le refroidissement par conduction avec le sol et convection avec l'air. La figure 16 rappelle de façon schématique l'ensemble de ces mécanismes. En faisant l'hypothèse d'une homogénéité de 30 température dans la largeur de la bande de roulement et sur un tour de roue, il est possible d'obtenir l'équation 42 2910662 unidimensionnelle de la chaleur dans le système de coordonnées polaires liée à la roue : aT 2 â2T at pcp axe pcp 5 T représente le profil de température dans l'épaisseur de la bande de roulement [ K]; 2 est la conductivité thermique de la gomme [W/m/ K]; p est la masse volumique [kg/m3] ; cp est la chaleur massique de la gomme [J/kg/ K]; 10 q est le terme de production de chaleur due aux pertes dans la gomme en [W/m3] ; x représente la direction radiale (i.e. dans l'épaisseur de la bande de roulement); 15 Effusivité : e=.Ppcp z Diffusivité : a= _ -2 [m2/s] pcp e Les conditions aux limites sont cependant différentes selon qu'on considère la surface extérieure 20 du pneu ou l'interface entre la bande de roulement et les nappes du sommet du pneu. Dans le premier cas, la surface du pneu voit ses conditions aux limites changer sur un tour de roue : hors de l'aire de contact, il existe un flux de surface due à 25 la convection avec l'air environnant; dans l'aire de contact, il y a un flux de surface lié à la conduction avec le sol et au frottement entre la gomme et le sol. A l'interface gomme/(sol+air), la condition aux limites de où 43 2910662 flux imposé à la surface du pneu peut s'écrire d'une manière formelle comme suit : -2-=~p [W/m2] 5 où ç est un flux de surface qui reste à expliciter. Les échanges par rayonnement thermique peuvent être négligés. En revanche, au niveau de l'interface entre la 10 bande de roulement et les nappes du sommet du pneu, l'hypothèse d'un flux nul peut être faite (condition adiabatique). Le calcul du terme q [W/m3] de production de chaleur due aux pertes dans la gomme peut être réalisé de 15 la façon suivante. Lorsque la gomme passe dans l'aire de contact, elle subit des déformations par compression et cisaillement qui sont source de chaleur. La puissance interne dissipée dans la gomme sur un tour de roue est calculée comme le 20 produit, par la fonction pertes P, de l'énergie fournie lors du passage dans l'aire de contact Wf et de la fréquence de tour de roue : q=PWff La densité d'énergie de déformation élastique subie 25 par la gomme dans l'aire de contact est explicitée en fonction des efforts longitudinaux, transversaux et de la charge du pneu, ce qui permet d'obtenir la formulation finale . 44 2910662 _ 1 V FX2+F.2+FZ2/3 4 P(Wf'T) 2Gk 2rrRo S,aDC2 où P(WW,T) est la fonction pertes, qui prend en compte le point de fonctionnement de la gomme en température et en 5 amplitude de sollicitation et qui peut être caractérisée par expérimentation spécifique. La température T correspond à la température caractéristique du mélange et détermine entre autres le niveau de pertes et de module. Sachant que les lois 10 pertes et module sont par exemple mesurées à une fréquence de 10Hz, la température T est en fait une température équivalente au sens de la loi WLF, afin d'avoir une estimation des pertes et module pour des fréquences de sollicitation différentes : 15 T =T.. +C21og1o(f /10)/[Cl -1og10(f /10)], où T1 est la température interne du mélange issue du calcul thermique global, et où f = V /(27cRo) est la 20 fréquence de rotation. Comme le comprendra aisément l'homme du métier à la lecture de la présente description, la même relation est utilisée pour relier le module de cisaillement de la gomme à la température interne Ti, ceci afin de prendre 25 en compte les mécanismes de rigidification du mélange lorsque la fréquence du tour de roue augmente. Pour le calcul du flux de conduction avec le sol, la bande de roulement et le sol peuvent être assimilés à deux murs semi-infinis mis en contact pendant 2910662 l'intervalle de temps Lad,. En supposant le contact parfait, le flux de conduction s'écrit : egesol 1 9,,,,d = -2 - (Ts ù Tsoi eg + esai \/tadcit 5 où Ts est la température de surface de la gomme. Le calcul du flux de convection avec l'air est rendu plus difficile par le fait que l'échange thermique avec l'air dépend fortement de la nature de l'écoulement 10 d'air autour du pneu. D'une manière générale, la modélisation des échanges convectifs est basée sur des formulations semi-empiriques. Dans le cas particulier du pneumatique, la formule suivante peut être utilisée : 15 'Pconv Cair 4.126(2R0 )ù0.2 770.805 (Ts ù Tair ) où Cair est une constante qui prend en compte l'effet de la convection forcée. Le calcul du flux de frottement doit rendre compte 20 du fait que le frottement constitué par le glissement de la gomme sur le sol est source de production de chaleur. En toute rigueur, la dissipation d'énergie se fait dans le massif de gomme sur des épaisseurs inférieures au millimètre. Dans le cas d'un contact sec, l'approximation 25 peut être faite que l'énergie est dissipée en extrême surface et qu'elle est modélisée par un flux de frottement. Le flux de frottement moyen dans la zone de glissement de l'aire de contact s'écrit comme suit : 46 2910662 ~frot = a Vg Fp (1- ppa)LxLy ent où a est un coefficient de partage de flux entre la gomme et le sol; une valeur a=l signifie que la totalité 5 du flux de frottement est orienté vers la gomme; une valeur a-0 signifie que la totalité du flux de frottement est orienté vers le sol; Fp est la composante de l'effort induit par le frottement entre la gomme et le sol; Vg est la vitesse de 10 glissement; et ppa est la proportion de points adhérents dans l'aire de contact. Le flux thermique moyen à la surface de la bande de roulement est défini comme la moyenne des différents flux, pondérée par la durée caractéristique pendant 15 laquelle ces flux sont effectifs lors d'un tour de roue, ce que traduit la relation : q) frot (1 - PPaltadc + V Gond tadc -F- ç0convtHadc _ tadc + tHadc 20 où tadc correspond à la durée de passage d'un élément de la bande de roulement dans l'aire de contact; où tHadc est la durée de passage d'un élément de la bande de roulement hors de l'aire de contact; et où (1-ppa) talc est la durée pendant laquelle un élément de la bande de roulement 25 glisse dans l'aire de contact. La figure 17 illustre de façon globale la mise en oeuvre opérationnelle d'un procédé conforme au modèle de pneu préféré utilisable dans le procédé de l'invention. 47 2910662 En amont de sa mise en oeuvre opérationnelle, la réalisation de ce modèle de pneu comprend une phase préliminaire de modélisation au cours de laquelle sont établis le modèle mécanique ou premier modèle, le modèle 5 thermique local ou deuxième modèle, et le modèle thermique global tour de roue, ou troisième modèle. Chaque modèle est établi par application de lois physiques connues et / ou construites par expérimentation spécifique, et prend la forme d'un système d'équations 10 telles que celles présentées ci-dessus pour chaque modèle. Pour rappel, le modèle mécanique fournit les expressions des efforts longitudinaux et des efforts transversaux qui sont transmis par le pneu entre le sol 15 et le véhicule, du couple d'auto-alignement qui est lié à l'intensité de ces efforts et à leur distribution dans l'aire de contact, et de l'équilibre des forces élémentaires de cisaillement et de glissement du pneu au point N, supposé unique dans le mode de réalisation 20 préféré de l'invention, de passage entre les zones de contact adhérent et glissant, ces expressions étant données en fonction des paramètres dynamiques qui sont liés aux conditions physiques de roulement et d'utilisation du pneu, des paramètres physiques 25 spécifiques du pneu, et de l'abscisse b du point de passage. Le modèle thermique local fournit l'expression des variations de la température de contact de la bande de roulement avec le sol depuis l'entrée jusqu'à la sortie 30 de l'aire de contact, ces variations dépendant notamment de la température périphérique de la bande de roulement, 48 2910662 de la température du sol, et du glissement de la bande de roulement sur le sol. Le modèle thermique global fournit le profil de température dans l'épaisseur de la bande de roulement et 5 les variations, sur une période de tour de roue, de la température périphérique de la bande de roulement et de la température interne du pneu, en fonction notamment des valeurs antérieurement connues ou estimées des températures périphérique et interne, du coefficient de 10 conduction thermique de la bande de roulement, et de phénomènes à composante thermodynamique tels que les déformations internes subies par le pneu, les échanges thermiques entre le pneu et son environnement, et le glissement de la bande de roulement sur le sol. 15 Par ailleurs, le modèle mécanique prend en compte le coefficient d'adhérence et le module de cisaillement sous forme de variables respectivement fonctions de la température périphérique et de la température interne. La phase d'exploitation du modèle de pneu, qui 20 s'appuie sur les résultats de la phase de modélisation, comprend des opérations d'affectation de valeurs numériques et des opérations de résolution. Les opérations d'affectation consistent essentiellement à attribuer des valeurs numériques aux 25 paramètres dynamiques et aux paramètres spécifiques, tandis que les opérations de résolution consistent essentiellement à déduire, en utilisant notamment le modèle mécanique et les valeurs attribuées aux paramètres dynamiques et aux paramètres physiques spécifiques, les 30 valeurs des efforts longitudinaux Fx, des efforts transversaux Fy, et du couple d'auto-alignement Mz. 49 2910662 La valeur de la température périphérique Ts, qui est obtenue par des opérations de résolution afférentes au modèle thermique global, est utilisée par des opérations de résolution afférentes au modèle thermique 5 local pour prendre en compte la dépendance en température du coefficient d'adhérence. D'autre part, la valeur de la température interne Ti, qui est également obtenue par des opérations de résolution afférentes au modèle thermique global, est 10 utilisée par des opérations de résolution afférentes au modèle mécanique pour prendre en compte la dépendance en température du module de cisaillement de la gomme du pneu. Le modèle mécanique inclut de préférence des 15 équations liées aux conditions d'équilibre des forces élémentaires de cisaillement et de glissement de la bande de roulement dans l'aire de contact, et sur la base desquelles s'effectue le couplage entre le modèle mécanique et le modèle thermique local du pneu. 20 En pratique, la mise en ouvre du modèle de pneu comprend avantageusement une phase itérative de calcul constituée d'une succession de cycles de calcul, chaque cycle de calcul comprenant au moins une opération de résolution afférente au modèle mécanique et une opération 25 de résolution afférente au modèle thermique local. La phase d'exploitation de ce modèle de pneu, qui fait suite à la phase de modélisation et qui est illustrée à la figure 17 pour le mode de réalisation le plus avantageux, comprend d'abord une opération 1 30 consistant à initialiser un indice compteur n utilisé pour la mesure du temps. 2910662 A l'opération 2, des valeurs initiales Tso et Tio sont attribuées à la température périphérique de la bande de roulement et à la température interne du pneu, par exemple en faisant l'hypothèse que le pneu est 5 initialement en équilibre thermique avec l'air ambiant. L'opération 3 incrémente le temps compté par le chronomètre d'un intervalle de temps correspondant au moins à celui qui est nécessaire pour la réalisation des calculs ultérieurs tels que décrits ci-après. 10 Des valeurs préalablement mesurées ou mémorisées sont ensuite attribuées aux paramètres dynamiques (opération 4) et aux paramètres spécifiques (opération 5). L'opération 6 consiste à initialiser un indice 15 compteur k utilisé pour compter le nombre de cycles de calcul successifs effectués à l'intérieur d'une boucle d'itération qui sera décrite ci-après. L'opération 7 constitue une phase préparatoire destinée notamment à permettre le calcul de grandeurs 20 contingentes dont la valeur peut être considérée comme constante pour les différents cycles de calcul d'une même boucle d'itération, ce qui permet d'éviter l'exécution répétée de ces calculs à chaque cycle de calcul d'une même boucle d'itération. 25 En particulier, :La phase préparatoire 7 est utilisée pour calculer les dimensions Lx et Ly de l'aire de contact, sa surface SAx, le profil de pression p (x) le long de l'aire de contact, de même que les rigidités Kx et Ky de la bande de roulement en fonction de la 30 température interne Ti à l'instant antérieur n-1, c'est-à-dire Tin-l. 51 2910662 Des valeurs estimées Fyo, Mzo et bo sont par ailleurs attribuées aux efforts latéraux Fy, au couple d'auto-alignement Mz, et à l'abscisse b du point N de passage entre la zone de contact adhérent et la zone de 5 contact glissant. Lorsque les paramètres d'entrée varient peu en fonction du temps, les valeurs estimées Fyo, Mzo et bo peuvent être constituées par les valeurs calculées à l'instant précédent. 10 Dans le cas contraire, l'abscisse bo de début de glissement est déterminée en supposant que le profil de pression dans la longueur de l'aire de contact est parabolique et en négligeant la rigidité de torsion et la rigidité du sommet. 15 Dans ce cas, l'équation de glissement (équation 3) a une solution analytique : _ -2 r +[Ky(S+a1)]2 -1 bo = a i 4 a2 3 ,uoFz /LY\ 20 Connaissant bo, l'effort Fyo et le couple Mm sont alors calculés à partir des équations 5 à 7. Il faut en outre vérifier que la position initiale obtenue pour le point de passage satisfait la contrainte sur le signe du cisaillement latéral de la bande de 25 roulement, à savoir (YK -YN)8' > 0 . Si tel n'est pas le cas, la solution envisagée n'a pas de sens physique. Il faut alors imposer YK ùYN = 0 , et la valeur initiale Mm du couple d'auto-alignement est fixée à O. 52 2910662 L'opération 8 incrémente l'indice compteur k et permet d'aborder un premier ou un nouveau cycle de calcul (étapes 9 et 10) de la phase itérative. Cette phase itérative permet de résoudre par 5 approximations successives, et à partir des valeurs antérieurement connues ou estimées de l'abscisse b du point de passage, des efforts transversaux Fy, et du couple d'auto-alignement Mz, de nouvelles valeurs de ces grandeurs b, Fy et Mz qui résolvent les équations 1 à 7 10 précédemment présentées, pour les valeurs attribuées aux paramètres dynamiques et aux paramètres spécifiques. Ces équations sont formulées en faisant apparaître l'abscisse b du point de passage entre la zone de contact adhérent et la zone de contact glissant, les efforts dans 15 l'aire de contact étant décomposés en deux contributions, à savoir un effort de cisaillement qui dépend des rigidités de la gomme, du sommet et de la carcasse du pneu, et un effort de frottement qui dépend de la loi de frottement. 20 L'abscisse b est calculée à partir des équations 1 à 3 et des valeurs de Fy et Mz estimées à l'itération précédente. Il s'agit d'une équation scalaire, dont la solution est bornée (- a S b 5 a). Le calcul de l'abscisse b est par exemple réalisé par combinaison 25 d'une bissection et d'une sécante. Si le déplacement proposé par la méthode sécante sort des bornes inférieures et extérieures, le modèle de pneu bascule sur un pas de bissection. Plusieurs solutions de l'abscisse b étant a priori 30 possibles, la solution retenue est celle qui satisfait la condition (YK -YN)8' > O. 53 2910662 b b Les intégrales f ,u(P(x),T,Vg)P(x)dx et f ,u(P(x),T,Vg)P(x)xdx -a -a liées au frottement sur le sol de la gomme de la bande de roulement sont par exemple calculées en utilisant les formules de quadrature de Gauss. 5 Pour le calcul des résidus du système constitué par les équations 4, 5, et 6+7, et le calcul de la convergence, il s'agit de résoudre le système non linéaire de n équations en n inconnues noté de manière formelle : F(x)=O. 10 Bien que plusieurs procédés itératifs soient possibles, le procédé optimal semble être constitué par le procédé itératif mixte de type Newton-Raphson/Broyden connu de l'homme de l'art. A chaque cycle de calcul d'ordre k, l'étape 9 15 comprend notamment une opération consistant à calculer une nouvelle valeur provisoire bk de l'abscisse du point de passage à partir des équations 1 à 3 d'équilibre des forces élémentaires et des valeurs antérieurement connues ou estimées FYk_1 et MZk_1 des efforts transversaux et du 20 couple d'auto-alignement. Dans le cas du modèle de pneu préféré, dans lequel le modèle mécanique est enrichi par le modèle thermique local et prend notamment en compte l'influence de l'évolution de la température de contact sur la valeur du 25 coefficient de frottement de la gomme du pneu sur le sol, l'étape 9 comprend également, pour chaque point de la zone de contact glissant, le calcul de la température de contact Tc et du coefficient de frottement en fonction de la valeur provisoire nouvellement calculée bk de 30 l'abscisse du point de passage, de la température périphérique Tsn_1 telle que connue à l'instant précédent, 54 2910662 et de la température du sol Tso1, la température périphérique Tsri._1 du pneu étant utilisée pour le calcul de la température de contact en entrée de la surface de contact. 5 Enfin, l'étape 9 comprend l'opération qui consiste à calculer, à partir de la nouvelle valeur provisoire bk de l'abscisse du point de passage, des valeurs des rigidités Kx et Ky de :La bande de roulement pour la température interne Tin-1 telle qu'elle est connue à 10 l'instant antérieur, de la valeur du coefficient de frottement/1, et des équations 1 à 7, de nouvelles valeurs Fyk et MZk des efforts transversaux et du couple d'auto-alignement, utilisables pour un éventuel cycle de calcul ultérieur. 15 L'étape 10 consiste à tester la convergence de la phase itérative. Par exemple, si les écarts respectifs entre d'une part les nouvelles valeurs bk, Fyk et MZk et d'autre part les valeurs provisoires bk_1, FYk_1 et MZk-1 issues du cycle 20 de calcul précédent sont inférieurs à des limites respectives telles que e, alors la phase itérative est interrompue. Dans le cas contraire, cette phase itérative est poursuivie par rebouclage en amont de l'opération 8. A l'interruption de la phase itérative, le modèle 25 thermique global tour de roue est utilisé (étape 11) pour calculer les nouvelles valeurs à jour Tsn et Tin des températures périphérique et interne en tenant compte des évolutions subies par ces températures, depuis la fin de la phase itérative précédente, sous l'effet de l'ensemble 30 des phénomènes à composante thermodynamique, tels que les déformations internes subies par le pneu, les échanges 2910662 thermiques entre le pneu et son environnement (air, sol), et le glissement de la bande de roulement sur le sol. On rappelle que la température Ts est la température de surface moyenne, à la fois dans la largeur 5 et dans la circonférence du pneu, de la bande de roulement, et que le calcul des températures Ti et Ts s'appuie sur une modélisation unidirectionnelle dans l'épaisseur de la bande de roulement.. Les équations du modèle thermique global sont 10 résolues par une approche classique de méthode par différences finies en utilisant un maillage spatial dans l'épaisseur de la gomme et une méthode de résolution dans le temps de type Runge-Kutta d'ordre 2. L'opération 12, mise en oeuvre à la sortie de la 15 phase itérative, consiste à délivrer, pour l'instant n, les valeurs des efforts longitudinaux et transversaux Fxn, Fyn, du couple d'auto-alignement Mzn, de la température interne Tin du pneu, et de la température périphérique Tsn de la bande de roulement. 20 Le procédé est alors rebouclé juste en amont de l'opération 3 d'incrémentation du chronomètre, et avant l'opération 4 de mise à jour des paramètres dynamiques qui permet de tenir compte des évolutions subies par ces paramètres pendant le temps d'exécution de la phase 25 itérative qui vient de se terminer. La nouvelle valeur Tin de la température interne du pneu sera utilisée pour l'opération 5 de mise à jour des paramètres spécifiques, ou au cours de la phase préparatoire 7, pour en déduire la nouvelle valeur de la 30 rigidité G* du mélange constitutif de la gomme du pneu, qui conditionne les valeurs des rigidités Kx et Ky. 56 2910662 Par ailleurs, la nouvelle valeur Tsn de la température périphérique de la bande de roulement sera utilisée lors de l'étape 9 pour le calcul de la température de contact Tc en entrée de la bande de 5 roulement. On comprend donc que le couplage entre le processus de détermination des efforts et le processus de détermination des températures intervient à deux niveaux, à savoir d'une part par la prise en compte du fait que la 10 température moyenne Ti de la bande de roulement influe sur la rigidité G* du mélange et donc sur les rigidités Kx et Ky de la bande de roulement, et d'autre part par la prise en compte du fait que la température périphérique Ts de la bande de roulement dans l'aire de contact influe 15 sur le coefficient d'adhérence entre la gomme et le sol. Le procédé de simulation de pneu qui vient d'être décrit est notamment applicable à la simulation en temps réel du comportement dynamique d'un véhicule doté d'un châssis et de plusieurs pneus en situation de roulement 20 sur le sol. Ce modèle de pneu peut alors par exemple être mis en oeuvre sur un ordinateur portable ayant les caractéristiques suivantes PC Windows XP, Pentium 4 3.6Ghz, 2Go de mémoire les performances du procédé de 25 simulation d'un pneu seul en situation de roulage dérive et en présence de glissement sont de 9000 calculs complets à la seconde, où chaque calcul correspond à la détermination des efforts longitudinaux et transversaux, du couple d'auto-alignement et des températures de 30 surface et interne du pneu à un instant donné en fonction des paramètres dynamiques et spécifiques. 57 2910662 Dans les applications véhicule, chacun des modèles mécanique, thermique local et thermique global, ou au moins le premier d'entre eux, est exploité pour chaque pneu et associé à un modèle dynamique de châssis dont un 5 exemple sera fourni ultérieurement. Chacun des modèles de pneu coopère avec ce modèle de châssis d'une part pour recevoir du modèle de châssis les valeurs des paramètres dynamiques ou de certains d'entre eux au moins, et d'autre part pour permettre au 10 modèle de châssis d'exploiter, pour chaque pneu, les valeurs des efforts longitudinaux, des efforts transversaux, et du couple d'auto-alignement obtenues par la mise en oeuvre des modèles de pneu. Les figures 18a et 18b illustrent, pour un 15 intervalle de temps n déterminé, la phase d'exploitation d'un autre modèle de simulation de pneu qui, bien que beaucoup moins économique que le précédent en calculs, offre lui aussi un grand réalisme. Dans cet autre modèle de pneu, l'aire de contact 20 est discrétisée par exemple sur N zones successives de sa longueur et l'équilibre des forces élémentaires de cisaillement et de glissement de la bande de roulement dans l'aire de contact (ADC) est examiné sur chaque élément de surface. 25 Plus précisément, chaque cycle de chaque phase itérative comprend N opérations successives d'analyse des conditions d'équilibre des forces élémentaires de cisaillement et de glissement sur les N surfaces élémentaires respectives différentes de l'aire de 30 contact, et chaque cycle de la phase itérative est interrompu lorsque l'ensemble des surfaces élémentaires considérées lors de ce cycle recouvre l'aire de contact. 58 2910662 Les opérations 2', et 4' à 8' de ce modèle sont respectivement les mêmes que les opérations 2 et 4 à 8 du modèle de la figure 17. L'opération 13 consiste à initialiser un indice 5 compteur i utilisé pour compter les opérations d'analyse des conditions d'équilibre des forces élémentaires sur les différentes surfaces élémentaires. L'opération 14 constitue une phase préparatoire destinée à fixer à zéro les valeurs des efforts 10 élémentaires transversaux et longitudinaux âM O et â540 s'exerçant en amont de la première surface élémentaire de l'aire de contact, et à fixer la valeur initiale de la température de contact Tc. L'opération 15 incrémente l'indice compteur i et 15 permet d'aborder l'analyse des efforts élémentaires dans une première ou une nouvelle surface élémentaire de l'aire de contact, les opérations 16 et 17 étant consacrées au calcul de ces efforts. Le test 18 vise à déterminer si les efforts 20 élémentaires de cisaillement de la bande de roulement sont ou non supérieurs à l'adhérence. Dans la négative, l'opération 19 fixe à zéro l'échauffement qui aurait été lié au glissement de la bande de contact sur le sol. 25 En revanche, en cas de glissement, l'opération 20 détermine avec précision les conditions de glissement et calcule le flux de chaleurlié à ce glissement conformément au modèle thermique local. Qu'il y ait ou non glissement, le test 21 est 30 effectué pour déterminer si la surface élémentaire considérée est ou non la dernière des N surfaces 2910662 59 élémentaires de la bande de roulement dans l'aire de contact discrétisée. Dans la négative, la surface suivante est analysée après retour à l'opération 15. 5 Sinon, l'opération 22 est effectuée pour calculer les efforts longitudinaux et transversaux ainsi que le couple d'auto-alignement par sommation des efforts élémentaires et des couples élémentaires obtenus pour les différentes surfaces élémentaires. 10 Le test 10', qui est équivalent au test 10 de la figure 17, conduit soit à retourner à l'opération 8' pour poursuivre la phase itérative, soit à fournir les valeurs des efforts et du couple d'auto-alignement pour le cycle de calcul considéré et à mettre en oeuvre les opérations 15 d'affectation et de résolution pour le modèle thermique global, telles que recensées de façon directement lisible sur la figure 18b. MODELISATION ET SIMULATION DU VEHICULE 20 La présente section décrit un mode possible de réalisation de l'étape (c) et de chaque étape du procédé de l'invention qui implique la simulation du comportement physique du véhicule en utilisant un modèle de châssis 25 complété par un modèle de pneu, ce dernier étant par exemple constitué par le modèle décrit dans la section précédente. La présente section de description fera usage des conventions d'écriture supplémentaires ci-après. 30 Tr Traînée aérodynamique g% Taux de glissement (Slip Ratio) av Angle volant 30 35 2910662 P Pression de gonflage d'un pneu (vecteur à quatre composantes) yr Angle de lacet (rotation autour de l'axe vertical du véhicule) 5 çt Vitesse de lacet (rotation autour de l'axe vertical du véhicule) t Temps x Position du véhicule (vecteur à deux composantes) 10 v Vitesse du véhicule (idem) a Accélération du véhicule (idem) R Rayon de courbure de la trajectoire m Masse du véhicule [Kg] hCDG Hauteur du centre de gravité [m) 15 hcDR Hauteur du centre de roulis [m] E Empattement du véhicule [m], avec E = et + e2 et Distance entre le centre de gravité et l'essieu avant [m] e2 Distance entre le centre de gravité et l'essieu 20 arrière [m] Voie Voie [m]. A priori, la voie n'est pas la même à l'arrière et à l'avant. Il faut distinguer : voieAv : Il [m] VOieAR 12 [m] 25 Par ailleurs, quand il sera nécessaire de se référer à une grandeur par roue, les notations suivantes seront utilisées . l'indice 1 désigne la roue avant gauche (par exemple Fxl) ; l'indice 2 désigne la roue avant droite (par exemple Fx2) ; l'indice 3 désigne la roue arrière gauche (par exemple Fx3) ; l'indice 4 désigne la roue arrière droite (par exemple Fx4). En réalité, la simulation d'un véhicule au sens de la présente description inclut implicitement un modèle du pilote en plus d'un modèle de châssis du véhicule et d'un 40 modèle de pneu. Cependant, le pilote est supposé suivre une trajectoire imposée, de sorte que les commandes appliquées par le modèle de pilote au véhicule sont 2910662 61 limitées aux sollicitations appliquées aux pneus, c'est-à-dire aux glissements et aux dérives. Un véhicule de compétition tel qu'une Formule 1 est simulé par un modèle de châssis et de pneus de type 5 quadricycle, alors qu'une moto est simulée par un modèle de type bicycle. Cependant, le modèle quadricycle est typiquement obtenu par extension du modèle bicycle, auquel sont rajoutées des relations décrivant le transfert de charge 10 latéral en virage (il y a donc une charge par roue et non plus par essieu) et des relations décrivant le carrossage par roue (inclinaison). La plupart des données du modèle véhicule peuvent être obtenues directement ou indirectement à partir de 15 mesures effectuées sur des véhicules réels. Le fonctionnement en régime permanent d'un véhicule bicycle, bien connu de l'homme du métier, est rapidement rappelé ci-après en référence à la figure 19. A un instant t, le bicycle suit une trajectoire 20 circulaire de rayon de courbure R, avec une vitesse V. Une accélération centrifuge s'exerce sur le centre de gravité du véhicule. Pour prendre le virage, le véhicule va devoir compenser cette accélération avec des efforts générés par les essieux avant et arrière. Ces 25 efforts seront produits par la mise en dérive des pneus avant et arrière. Les vecteurs vitesses du centre de gravité G et des centres des roues sont représentés par les flèches de direction tangente à la trajectoire circulaire de centre 30 0 et de rayon R. La distance entre le centre de gravité G et l'essieu avant (A) est notée 11 ; la distance entre le centre de 10 62 2910662 gravité G et l'essieu arrière (B) est notée 12 ; la dérive avant est notée 8l; la dérive arrière est notée 82 ; la dérive de caisse, angle entre le vecteur vitesse en G et l'orientation du véhicule est notée 8G ; l'angle 5 de braquage de l'essieu avant est noté fi. Dans ces conditions, l'angle de dérive de caisse se calcule par : 8G=82R Et l'angle de braquage de l'essieu avant est donné par : Dans le modèle quadricycle, le carrossage par roue peut être modélisé (de façon purement heuristique) par 15 une fonction affine de l'accélération transversale aT suivant la formule : y = kyl aT +k y2, dans laquelle les paramètres k;4etky2peuvent être obtenus à partir de mesures sur des véhicules réels, le signe de 20 l'accélération transversal dépendant de la localisation droite ou gauche de la concavité du virage. Dans la mesure où la charge aérodynamique (ici notée Fz par dérogation aux notations choisies) qui est appliquée sur la voiture évolue avec la vitesse, elle est 25 prise en compte sous la forme d'un polynôme d'ordre 4 dépendant de cette vitesse, tel que : Fz = kz0 + kzly + kz2v2 + kz3v3 +kz4v4, ce polynôme étant établi par approximation de relevés expérimentaux. 2910662 63 La charge totale se répartit ensuite sur les essieux avant et arrière suivant la balance aérodynamique B : Fz = FzAV + FzAR FzAf_=B FZ la balance aérodynamique étant supposée être une fonction 5 affine de la vitesse, donc de la forme : B=kBo+kBly Il est d'usage de s'intéresser aux valeurs de la balance aérodynamique à 150 km/h et 250 km/h. Si les Fzi sont les charges par roue mesurées sur 10 les différentes roues du véhicule, et si m la masse (éventuellement corrigée de la consommation de carburant), alors la charge aérodynamique totale est donnée par : 4 1Fz; -9.81 m i=l 15 Il convient également de prendre en compte la traînée, qui est une force aérodynamique longitudinale s'opposant à l'avancée du véhicule, particulièrement importante à haute vitesse. Cette force est en fait 20 retranchée à la somme des efforts Fxi dans l'équation d'équilibre longitudinal du véhicule. La traînée est par exemple modélisée par un polynôme d'ordre 2 du type : Tr = kTro +kTr]V+kTr2Vz . 25 Le moteur du véhicule peut, quant à lui, être simplement modélisé par une puissance maximum disponible : environ 850 ch. avec un moteur V10, et 650 ch. avec un moteur V8 et par exemple un unique régime moteur. 64 2910662 En conséquence, la puissance délivrée par l'essieu arrière est limité par la puissance maximum exprimée en chevaux, à savoir : (Fx3 + Fx4) v REVENDICATIONS . 1. Procédé pour sélectionner, parmi un ensemble de configurations possibles, une configuration optimisée de pneumatiques pour équiper un véhicule destiné à courir sur un circuit en suivant une trajectoire prédéterminée, orientée dans un sens unique de parcours et incluant des virages, ce procédé comprenant au moins les étapes consistant à : -développer un modèle physique du véhicule équipé d'une première configuration de pneus; - simuler le parcours de la trajectoire par le véhicule modélisé, en optimisant la vitesse de ce véhicule; - mémoriser le temps de parcours mis par le véhicule pour parcourir au moins une fois la trajectoire; - comparer le temps de parcours à au moins un temps de référence et produire un résultat de comparaison; et - sélectionner ou non la première configuration comme configuration optimisée en fonction du résultat de comparaison. 2. Procédé suivant la revendication 1, caractérisé en ce que le modèle du véhicule inclut un modèle thermomécanique des pneus équipant ce véhicule. 3. Procédé suivant l'une quelconque des revendications précédentes, caractérisé en ce qu'il comprend les étapes consistant à : (a) recenser un ensemble de paramètres physiques tendant à limiter la vitesse adoptée par le véhicule sur la trajectoire, cet ensemble comprenant des paramètres constants et des paramètres dynamiques, ces paramètres étant liés aux différents pneus possibles, au châssis du 68 2910662 véhicule et/ou aux conditions de roulement, et incluant une image numérique de la trajectoire, constituée par des coordonnées numérisées de points successifs de cette trajectoire; 5 (b) développer un modèle de pneu décrivant le comportement physique de chaque pneu et faisant intervenir un premier sous-ensemble de l'ensemble des paramètres physiques; (c) développer un modèle de véhicule faisant 10 intervenir à la fois le modèle de pneu et un deuxième sous-ensemble de l'ensemble des paramètres physiques, et décrivant le comportement physique du véhicule en fonction au moins de sollicitations des pneus, ces sollicitations comprenant des glissements et des dérives; 15 (d) découper l'image numérique de la trajectoire en tronçons adjacents se succédant dans le sens de parcours et dont chacun comprend un ensemble d'au moins trois points de trajectoire incluant un point initial et un point terminal de ce tronçon se succédant dans cet ordre 20 dans le sens de parcours, le point terminal étant localisé dans une entrée de virage; (e) associer respectivement aux points initial et terminal de chaque tronçon des valeurs limites de la vitesse du véhicule; 25 (f) sélectionner, en tant que configuration active, une première configuration de pneus; (g) sélectionner, en tant que tronçon actif, un premier tronçon de trajectoire; (h) simuler, au moyen du modèle de véhicule, le 30 parcours d'une première partie au moins du tronçon actif de trajectoire, depuis son point initial, par le véhicule équipé de la configuration active de pneus, en optimisant 69 2910662 les sollicitations des pneus pour maximiser l'accélération subie par le véhicule à partir de la vitesse limite associée au point initial de ce tronçon actif; 5 (i) mémoriser les vitesses du véhicule aux différents points de la première partie du tronçon actif et les temps de parcours entre ces différents points; (j) sélectionner, en tant que section active du tronçon actif, la portion de trajectoire définie entre un 10 premier et un deuxième points d'une paire de points du tronçon actif se succédant dans cet ordre dans le sens de parcours, le deuxième point de cette paire étant constitué par le point terminal du tronçon actif; (k) sélectionner, en tant que vitesse cible, la 15 vitesse limite associée au deuxième point de la section active; (1) simuler, au moyen du modèle de véhicule, le parcours de la section active du tronçon actif par le véhicule équipé de la configuration active de pneus, en 20 optimisant les sollicitations des pneus et la vitesse adoptée par le véhicule au premier point de la section active pour atteindre le deuxième point de la section active à la vitesse cible, tout en minimisant le temps de parcours sur cette section active; 25 (m) mémoriser le temps de parcours de la section active par le véhicule; mémoriser, en tant que vitesse limite, la vitesse du véhicule au premier point de la section active; et mémoriser, en tant qu'élément d'une deuxième partie du tronçon actif, la section active 30 parcourue; (n) vérifier si la vitesse limite du véhicule au premier point de la section active dépasse ou non la 70 2910662 vitesse, mémorisée à l'étape (i), qui est adoptée en ce même point par le véhicule lors du parcours de la première partie du tronçon actif; (o) dans la négative, sélectionner, en tant que 5 nouvelle section active du tronçon actif, une portion de trajectoire située en amont, par rapport au sens de parcours, de la section active chronologiquement précédente, et définie entre un premier et un deuxième points d'une paire de points du tronçon actif se 10 succédant dans cet ordre dans le sens de parcours, le deuxième point de cette paire étant constitué par le premier point de la section active chronologiquement précédente, et répéter les étapes (k) à (n); (p) dans l'affirmative, mémoriser, en tant que temps 15 de parcours total du tronçon actif, la somme du temps de parcours total sur la deuxième partie de ce tronçon actif et des temps de parcours mémorisés sur la fraction de la première partie de ce tronçon qui ne recouvre pas la deuxième partie; 20 (q) sélectionner, en tant que nouveau tronçon actif, le tronçon de trajectoire qui succède au tronçon précédemment actif dans le sens de parcours; (r) répéter les étapes (h) à (q) jusqu'à parcourir au moins une fois la trajectoire et obtenir un temps de 25 parcours cumulé et stabilisé de la trajectoire par le véhicule. 4. Procédé suivant la revendication 3, caractérisé en ce qu'il comprend en outre les étapes consistant à (s) sélectionner, en tant que nouvelle configuration 30 active, une deuxième configuration de pneus; (t) répéter les étapes (g) à (r); 2910662 71 (u) comparer les temps de parcours cumulés et stabilisés de la trajectoire, respectivement obtenus avec les première et deuxième configurations de pneus; et (v) sélectionner, en tant que configuration 5 optimisée, la configuration de pneus conduisant au plus court temps de parcours cumulé et stabilisé de la trajectoire, le temps de parcours obtenu avec chaque configuration de pneus servant ainsi de temps de référence pour le temps de parcours obtenu avec l'autre 10 configuration de pneus. 5. Procédé suivant l'une quelconque des revendications 3 et 4, caractérisé en ce que le modèle de pneu est un modèle thermomécanique, en ce que l'étape (j) est précédée d'une étape (w0) consistant à sélectionner, 15 en tant que conditions thermiques actives, des conditions thermiques prédéterminées, en ce que l'étape (1) est réalisée en appliquant aux pneus les conditions thermiques actives, en ce que l'étape (m) inclut l'opération élémentaire supplémentaire qui consiste à 20 mémoriser les sollicitations des pneus sur chaque section active, et en ce que l'étape (p) est immédiatement précédée d'une boucle d'itération comprenant les étapes consistant à : (wl) en cas de succès de la vérification réalisée à 25 l'étape (n), simuler au moyen du modèle de véhicule le parcours par ce véhicule équipé de la configuration active de pneus, section active après section active dans le sens de parcours, de la deuxième partie du tronçon actif, en appliquant aux pneus sur chaque section active 30 les sollicitations mémorisées à l'étape (m) pour cette section active, et en déduire de nouvelles conditions 72 2910662 thermiques d'utilisation des pneus sur la deuxième partie du tronçon actif; (w2) vérifier si l'écart entre les nouvelles conditions thermiques et les conditions thermiques 5 actives s'inscrit ou non dans une limite de convergence; et (w3) dans la négative, mémoriser, en tant que conditions thermiques actives, les nouvelles conditions thermiques et répéter les étapes (j) à (w2). 10 6. Procédé suivant l'une quelconque des revendications précédentes combinée à la revendication 3, caractérisé en ce que l'étape (d) est mise en oeuvre en fixant le point terminal de chaque tronçon au maximum de courbure d'un virage. 15 7. Procédé suivant l'une quelconque des revendications précédentes combinée à la revendication 3, caractérisé en ce qu'il comprend une opération supplémentaire, mise en oeuvre à chaque nouveau parcours de la trajectoire et consistant à déplacer le point 20 terminal de chaque tronçon vers l'amont, en direction d'une zone de moindre glissement des pneus. 8. Procédé suivant l'une quelconque des revendications précédentes combinée à la revendication 3, caractérisé en ce que le modèle de pneu développé à 25 l'étape (b) décrit le comportement physique de chaque pneu équipant le véhicule en situation de roulement sur le sol, avec lequel la bande de roulement du pneu présente une aire de contact incluant au moins une zone de contact adhérent et au moins une zone de contact 30 glissant, en ce que l'étape (b) comprend au moins des opérations de modélisation, en ce que chaque étape de simulation (h, 1, wl) comprend des opérations 73 2910662 d'affectation de valeurs numériques et des opérations de résolution, les opérations de modélisation consistant au moins à établir, par application de lois physiques connues et / ou construites par expérimentation 5 spécifique, et en tant que premier modèle, un modèle des efforts longitudinaux (Fx), des efforts transversaux (Fy), et d'un couple d'auto-alignement (Mz), transmis par le pneu entre le sol et le véhicule, en fonction de paramètres physiques spécifiques, incluant un coefficient 10 d'adhérence et un module de cisaillement de la gomme du pneu, et en fonction de paramètres dynamiques liés aux conditions physiques de roulement et d'utilisation du pneu, le couple d'auto-alignement étant lié à l'intensité des efforts longitudinaux et transversaux et à leur 15 distribution dans l'aire de contact, les opérations d'affectation consistant au moins à attribuer (4) des valeurs numériques aux paramètres dynamiques et aux paramètres spécifiques (5), et les opérations de résolution consistant au moins à déduire, en utilisant au 20 moins le premier modèle et les valeurs attribuées aux paramètres dynamiques et aux paramètres physiques spécifiques, des valeurs des efforts longitudinaux (Fx), des efforts transversaux (Fy), et du couple d'auto- alignement (Mz), en ce que les opérations de modélisation 25 comprennent aussi l'établissement, en tant que second modèle, d'un modèle local d'échauffement exprimant des variations d'une température de contact de la bande de roulement avec le sol depuis l'entrée jusqu'à la sortie de l'aire de contact par suite du contact et du 30 glissement de la bande de roulement avec le sol, et l'établissement, en tant que troisième modèle, d'un modèle global d'échauffement et de flux thermique, ce 74 2910662 troisième modèle exprimant des variations, sur une période d'au moins un tour de roue, d'une température périphérique de la bande de roulement et d'une température interne du pneu en fonction de valeurs 5 antérieurement connues ou estimées des températures périphérique et interne, d'un coefficient de conduction thermique de la bande de roulement, et de phénomènes à composante thermodynamique tels que déformations internes subies par le pneu, échanges thermiques entre le pneu et 10 son environnement, et glissement de la bande de roulement sur le sol, en ce qu'au moins le coefficient d'adhérence et le module de cisaillement interviennent dans le premier modèle comme des variables respectivement fonctions de la température périphérique et de la 15 température interne, en ce que la valeur de la température périphérique (Ts), obtenue par des opérations de résolution afférentes au troisième modèle, est utilisée par des opérations de résolution afférentes au deuxième modèle pour prendre en compte la dépendance en 20 température du coefficient d'adhérence, et en ce que la valeur de la température interne (Ti), obtenue par des opérations de résolution afférentes au troisième modèle, est utilisée par des opérations de résolution afférentes au premier modèle pour prendre en compte la dépendance en 25 température du module de cisaillement de la gomme du pneu. 9. Procédé de simulation suivant la revendication 8, caractérisé en ce que le premier modèle inclut des équations liées à des conditions d'équilibre des forces 30 élémentaires de cisaillement et de glissement de la bande de roulement dans l'aire de contact, en ce que ce procédé comprend au moins une phase itérative de calcul 75 2910662 constituée d'une succession de cycles de calcul, et en ce que chaque cycle de calcul comprend au moins une opération de résolution afférente au premier modèle et une opération de résolution afférente au deuxième modèle. 5 10. Procédé de simulation suivant la revendication 9, caractérisé en ce que le premier modèle est établi en considérant que l'aire de contact comprend une unique zone de contact adhérent et une unique zone de contact glissant séparées l'une de l'autre par un point de 10 passage, en ce que le premier modèle prend la forme d'un système d'équations exprimées au moins en fonction des paramètres dynamiques, des paramètres spécifiques, et de l'abscisse (b) du point de passage, en ce que chaque phase itérative est dédiée à des phénomènes apparaissant 15 pendant un intervalle de temps élémentaire correspondant, et en ce que chaque phase itérative est mise en oeuvre pour résoudre, par approximations successives et au moins à partir de valeurs antérieurement connues ou estimées de l'abscisse (bo) du point de passage, des efforts 20 transversaux (Fyo), et du couple d'auto-alignement (Mzo), de nouvelles valeurs de l'abscisse (b) du point de passage, des efforts transversaux (Fy), et du couple d'auto-alignement (Mz) qui résolvent le système d'équations du premier modèle pour les valeurs attribuées 25 aux paramètres dynamiques et aux paramètres spécifiques, ce dont il résulte que les efforts longitudinaux, les efforts transversaux, et le couple d'auto-alignement du pneu peuvent être calculés en temps réel pendant le roulement du véhicule. 30 11. Procédé de simulation suivant la revendication 10, caractérisé en ce que chaque nouveau cycle de calcul 76 2910662 de chaque phase itérative comprend les opérations (9) consistant au moins à : - calculer une nouvelle valeur provisoire de l'abscisse du point de passage à partir des équations d'équilibre 5 des forces élémentaires et des valeurs antérieurement connues ou estimées des efforts transversaux et du couple d'auto-alignement; - calculer, à partir de la nouvelle valeur provisoire de l'abscisse du point de passage et des équations liant les 10 efforts transversaux et le couple d'auto-alignement aux paramètres dynamiques, aux paramètres spécifiques, et à l'abscisse du point de passage, de nouvelles valeurs des efforts transversaux et du couple d'auto-alignement utilisables pour un éventuel cycle de calcul ultérieur; 15 - interrompre conditionnellement (10) ladite phase itérative au moins lorsque l'écart entre la nouvelle valeur provisoire de l'abscisse du point de passage et la valeur antérieurement connue ou estimée de cette abscisse est inférieur à une limite de précision prédéterminée; et 20 - à l'interruption de ladite phase itérative, affecter (12) aux efforts longitudinaux et transversaux et au couple d'auto-alignement, en tant que valeurs pour cette phase, les nouvelles valeurs des efforts longitudinaux et transversaux et du couple d'auto-alignement obtenues au 25 dernier cycle de calcul. 12. Procédé de simulation suivant la revendication 10 ou 11, caractérisé en ce qu'il comprend en outre une opération de mise en oeuvre après la fin de chaque phase itérative et consistant à mettre au moins à jour (4) les 30 paramètres dynamiques pour tenir compte des évolutions subies par ces paramètres pendant le temps d'exécution de 77 2910662 la phase itérative, et à engager une nouvelle phase itérative (8-10). 13. Procédé de simulation suivant l'une quelconque des revendications 8 à 12, caractérisé en ce que les 5 opérations de résolution afférentes au troisième modèle sont effectuées en dehors de chaque phase itérative. 14. Procédé suivant l'une quelconque des revendications 8 à 13, caractérisé en ce que chacun des premier, deuxième et troisième modèles est exploité pour 10 chaque pneu et associé à un modèle dynamique de châssis, en ce que le modèle de châssis fournit auxdits premier, deuxième et troisième modèles, pour chaque pneu, les valeurs de certains au moins des paramètres dynamiques, et en ce que le modèle de châssis exploite, pour chaque 15 pneu, les valeurs des efforts longitudinaux, des efforts transversaux, et du couple d'auto-alignement obtenues par la mise en oeuvre desdits premier, deuxième et troisième modèles. 20