Brevet FR2934907A1 — publié le 20100212

MÉTHODE DE CRÉATION D'UN MODÈLE DE DYNAMIQUE VIBRATOIRE D'UN ENSEMBLE PNEUMATIQUE ROUE CAVITE Domaine de l'invention [0001] L'invention concerne une méthode de création d'un modèle de prévision de la dynamique d'un ensemble Pneumatique/Roue/Cavité. [0002] Ce modèle vise d'une part à étudier le comportement vibratoire de l'ensemble pneumatique/roue/cavité, d'autre part à être intégré dans un modèle d'assemblage avec un système de liaison au sol d'un véhicule pour étudier le comportement vibratoire du Io véhicule ou d'un sous-ensemble du véhicule incluant l'ensemble-monté. Arrière plan technologique [0003] Dans ce qui suit : [0004] � Pneumatique � désigne l'élément de liaison entre le sol et le site d'accrochage à la roue. 15 [0005] � Roue � désigne l'élément de liaison entre le pneumatique et l'attache de roue du véhicule. [0006] � Cavité � désigne le volume de fluide, en général mais pas exclusivement de l'air, contenu entre le pneumatique et la roue. [0007] � Ensemble Monté � désigne le système constitué du pneumatique, de la roue et 20 de la cavité fluide. [0008] En dynamique vibratoire, un système mécanique se met classiquement sous la forme : KTOTOUTO +Crorouro + MTOTOUTO = fro [0009] Dans cette équation, l'indice TO fait référence au nombre total de degrés de liberté 25 de la structure, généralement élevé dans le cas de calculs utilisant la technique des éléments finis ; u,o et fTo sont respectivement les vecteurs des variations de -2- déplacements et de forces généralisées ; Kroro , Croro et MTOTO sont les matrices de rigidité, d'amortissement et de masse de la structure. [0010] Les techniques de sous-structuration dynamique sont bien connues de l'homme de l'art. Elles permettent de réduire l'information mécanique contenue dans les modèles éléments finis, au juste nécessaire pour représenter le comportement vibratoire de la structure. On passe classiquement de modèles à plusieurs centaines de milliers de degrés de liberté à moins d'une centaine. Ces techniques sont pratiquées couramment dans l'industrie. Elles sont particulièrement intéressantes pour assembler des sous-systèmes complexes issus de plusieurs partenaires. L'assemblage des systèmes se fait à travers leurs interfaces, le reste de la structure interne étant décrit par la liste des modes propres. L'application au pneu est particulièrement adaptée (malgré la présence de matériaux viscoélastiques) dans la mesure où d'une part la densité modale reste relativement faible (30 modes sont nécessaires pour couvrir la gamme de fréquences 0-200 Hz) et d'autre part les interfaces sont physiquement évidentes et parfaitement définies (centre-roue et aire de contact). [0011] En référence de l'état de l'art, on peut citer l'article de D. Tsihlas, T. Lacroix, B. Clayton, � A Comparison of Two Sub-Structuring Techniques for Representing the Modal Properties of Tires �, Tire Science and Technology, TSTCA, Vol. 29, N°1, Jan-March 2001, pp. 23-43. [0012] En particulier, on retient de cet article l'écriture finale du modèle sous-structuré selon la technique dite Craig-Bampton : K SSO û +C ssn ll + M ssn = f [0013] Avec f le vecteur efforts généralisés s'appliquant au système et i le vecteur déplacement généralisé décrivant les mouvements du système. [0014] L'originalité de cette formulation tient dans le choix des degrés de liberté ir . [0015] On décompose en effet ir en deux grandeurs distinctes, û = où sont les coordonnées modales, se référant aux intensités des modes propres de la structure. Typiquement pour un pneumatique, la taille de ce vecteur, correspondant au nombre de modes retenus, vaut entre 30 et 100. sont les degrés de liberté géométriques des interfaces externes du système. Par définition de ces interfaces, elles regroupent l'ensemble des noeuds géométriques susceptibles d'être soumis à des efforts extérieurs au système. [00161 Dans le cas du pneumatique, il est courant de retenir en tant qu'interface externe le centre-roue, usuellement mais pas exclusivement considéré comme un corps rigide ; et une aire de contact, dont le nombre de degrés de liberté retenu dépend de la complexité de l'interaction entre le sol et le pneumatique dont l'homme du métier souhaite rendre compte avec le modèle. [00171 Le grand avantage de cette écriture réside dans la condensation de l'information vibratoire sur les seules coordonnées modales. C'est grâce à cette technique que l'on passe classiquement de modèles à plusieurs dizaines de milliers de degrés de liberté à des modèles réduits à moins d'une centaine de degrés de liberté, tout en conservant une représentation fidèle du comportement vibratoire de la structure. [00181 Par la suite, on appellera � modèle � l'ensemble des trois matrices K sSD , CSSD et ms" respectivement matrices de rigidité, d'amortissement visqueux et de masse du système. [00191 Pour préciser les écritures qui seront utilisées par la suite, il faut noter que ces matrices peuvent prendre la forme suivante : KSSD = Kn K/E CssD = CH CIE et MSSD _ MI/ MIE _KEI KEE_' CEI CEE _ME/ MEE_ [00201 Les indices I se réfèrent aux coordonnées modales et les indices E se réfèrent aux degrés de liberté géométriques des interfaces externes. [00211 Classiquement, on adjoint à ce modèle deux termes (D,N_, et A,N_E (notés par la suite cl) et A) qui permettent de reconstituer la géométrie du système initial selon l'équation : 3 a, uE uT0 - ~UEX OEx / 1 X E l{ u I les degrés de liberté géométriques de la partie de la structure initiale que l'on souhaite représenter ; et u1 sont les degrés de liberté représentés par les modes et lEx les degrés de liberté d'interface dans l'écriture sous-structurée. [0022] Usuellement, ce type de modèle est utilisé pour des calculs dans le domaine fréquentiel sous l'hypothèse que les degrés de liberté sont de la forme u(t) = U(r))e'" avec w, pulsation d'excitation du système. On peut alors mettre les matrices du modèle sous la forme : K = K(w) = KSSD + iCSS°w MSS°w2 [0023] Par la suite, K ou les matrices KSS° CSS° et MSSD constituent indifféremment le modèle. [0024] On notera pour l'état de l'art comme pour l'invention que toutes les grandeurs constitutives du modèle sont complexes. [0025] Par construction, le modèle tel que décrit dans le document en référence prend en 15 compte la dynamique du pneumatique entre l'aire de contact et le centre-roue. [0026] Il omet toutefois de prendre en compte les dynamiques de la roue et de la cavité fluide qui se trouvent pourtant incluses entre l'aire de contact et le centre de roue. [0027] L'homme de métier constate alors par comparaison entre le calcul et la mesure que le modèle de pneumatique seul ne représente pas toute l'activité vibratoire dès 150 Hz 20 pour des valeurs usuelles de pneumatique et de roue. [0028] Pour introduire le comportement vibratoire de la cavité fluide dans la modélisation, une technique peut alors consister à résoudre dans le domaine temporel les équation d'interaction entre le fluide et la structure en imposant comme contrainte au fluide de suivre à tout moment le déplacement de la paroi du solide en contact tel que 25 décrit dans la demande de brevet US20070137290. Cette technique, certes intéressante, est toutefois très coûteuse en temps de simulation et capacité de machine de calcul. Par mu uro avec -5- ailleurs, elle n'autorise pas aisément l'interaction avec un modèle du comportement vibratoire du véhicule. Enfin elle n'intègre pas le comportement vibratoire de la roue dont il a été vu qu'il participait également à l'activité vibratoire. [0029] En alternative, W.B. Clayton, R. Saint-Cyr, dans l'article � Incorporation of the Tire Air Cavity Resonance into The Modal Tire Model �, Proceedings IMAC XVI, January 1998, proposent trois types de formulations pour introduire la dynamique du fluide interne de la cavité dans le modèle de pneumatique. Les modèles déduits de ces formulations sont in fine d'une forme proche voire identique de celle du modèle de pneumatique seul. Les méthodes décrites dans cet article présentent donc les avantages d'utilisation du modèle de pneumatique initial tout en ajoutant une représentation fidèle du comportement vibratoire du fluide interne à la cavité. Comme la description de l'invention s'appuiera pour partie sur les écritures utilisées dans cet article, on en rappelle ici les grandes lignes. [0030] Selon cet article, on peut mettre en équation le comportement vibratoire du fluide interne de la cavité pour aboutir à une formulation généralisée du modèle de cavité sous la forme de deux grandeurs : d'une part la liste S2? des pulsations propres de la pression dans le fluide, d'autre part, la grandeur notée (K'Ff) dans l'article, que l'on notera I'fs par la suite, qui représente le couplage entre les modes du fluide et la structure mécanique à l'interface entre le fluide et la structure. [0031] Selon l'état de l'art, en reprenant les écritures rappelées plus haut pour le modèle de la structure, on écrit alors le modèle couplé fluide/structure sous la forme : - PC K = TFfs. AFIS K 1E K 1E K EE [0032] Malgré des avantages évidents, la démarche présentée dans cet article n'intègre toujours pas le comportement vibratoire de la roue dont il a été vu qu'il participait 25 également à l'activité vibratoire. [0033] En alternative intéressante, le document JP2005161945 A propose également de s'appuyer pour partie sur une formulation de type Craig-Bampton du pneumatique et dans -6- une certaine mesure de la roue. Selon les revendications de ce document, on peut alors assembler les deux termes pour mener une simulation dans le domaine temporel. La technique proposée pour assembler les modèles de pneumatique et de roue présente toutefois l'inconvénient d'imposer la compatibilité des maillages à l'interface entre les deux structures. Ceci suppose en général une forte coordination entre les parties responsables de chacun des modèles. On peut imaginer même si ce n'est pas revendiqué, que cette approche permet d'améliorer la qualité de la simulation. Il manque en outre le comportement vibratoire de la cavité fluide. Description de l'invention [0034] L'objet de l'invention est une méthode de création d'un modèle de prévision de la dynamique vibratoire d'un ensemble monté constitué d'un pneumatique monté sur une roue et comportant un fluide dans la cavité interne définie par le pneumatique et la roue ou � cavité �. [0035] Cette méthode comporte les étapes suivantes : (1) définir des noeuds d'interface communs au pneumatique et à la roue ; (2) créer un modèle sous-structuré dynamique du pneumatique condensé aux ù P noeuds d'interface communs (K ) ; (3) créer un modèle sous-structuré dynamique de la roue condensé aux noeuds ù d'interface communs (K R ) ; (4) créer un modèle de l'activité vibratoire de la cavité (Q?, Ffs'', rfsR) ; et (5) assembler le modèle de la cavité aux modèles des sous-structures mécaniques ùPRC pneumatique et roue (K ). [0036] Le modèle résultant a l'avantage de décrire de façon très satisfaisante le comportement vibratoire de l'ensemble monté incluant notamment les modes du pneumatique, les modes de roue, les modes de la cavité fluide et leurs éventuels couplages sur une bande de fréquence d'au moins 0-500 Hz. Cette bande de fréquence est particulièrement utile pour la prévision du comportement vibro-acoustique d'un véhicule équipé des ensembles montés modélisés. -7- [0037] Avantageusement, la méthode comporte comme dernière étape, une étape de sous- - PRC structuration du modèle assemblé (K ) pour obtenir un modèle sous-structuré ù dynamique de l'ensemble monté K Fine! . [0038] Cette étape optionnelle permet d'augmenter l'efficacité du modèle en se ramenant 5 à un format strictement identique à celui du modèle de pneumatique seul. [0039] Pour définir des noeuds d'interface communs au pneumatique et à la roue, on effectue de préférence les étapes suivantes : - on établit un maillage éléments finis du pneumatique ; - on définit dans le maillage du pneumatique les noeuds (X) d'interfaces physiques du t o pneumatique et de la roue ; - on établit un maillage éléments finis de la roue ; - on définit dans le maillage de la roue les noeuds (Y) d'interfaces physiques de la roue et du pneumatique ; et on définit des points de condensation (NJ) communs au pneumatique et à la roue 15 avec pour chacun desdits points de condensation, sa position géométrique et sa relation de condensation entre ledit point de condensation et le ou les noeuds associés (X) d'interfaces physiques pneumatique ainsi que le ou les noeuds associés (Y) d'interfaces physiques roue. [0040] Cette méthode a l'avantage d'autoriser en particulier l'utilisation de maillages 20 incompatibles dans les calculs de roues et de pneumatiques qui peuvent de ce fait être travaillés séparément. En respectant cette procédure de définition des interfaces, les positions géométriques des noeuds du maillage du pneumatique et du maillage de la roue peuvent être différents, mais les modèles créés aux étapes ultérieures seront automatiquement compatibles. 25 [0041] De préférence, on décrit la géométrie du pneumatique et de la roue dans un nombre donné (NP) de plans méridiens ; et dans chaque plan méridien et pour chaque zone de contact pneumatique et roue, on assujettit tous les noeuds du contact pneumatique et roue à un déplacement de corps rigide. Le nombre de plans méridiens est avantageusement compris entre 25 et 100. -8- [0042] Pour créer un modèle sous-structuré dynamique du pneumatique condensé aux noeuds d'interface communs, on construit avantageusement une matrice sous-structurée -P (K ) en séparant les composantes : - représentant le comportement vibratoire des modes internes (K11") ; - représentant les noeuds (NJ) d'interfaces communs avec la roue (K,vjNJP) ; et - représentant les noeuds d'interfaces externes différents des noeuds communs avec la P roue (KEE ) ; pour obtenir : -P -P -P Ku KINJ KJE -P -P -P -P K = K NJ~ K NJNJ K NJE -P -P -P KE! K ENJ K EE [0043] Le modèle résultant KR décrit le comportement vibratoire du pneumatique et est compatible avec un couplage par les interfaces avec une roue. [0044] Pour pouvoir coupler le modèle de pneumatique avec le modèle de cavité fluide, on calcule également les termes cDP etAPqui permettent de reconstituer la géométrie du système initial du pneumatique selon l'équation : uPC =~PaP +AP i,' avec u" ' les degrés de liberté du maillage initial en contact avec la cavité fluide selon l'écriture rappelée précédemment. [0045] Pour créer un modèle sous-structuré dynamique de la roue condensé aux noeuds R d'interface communs, avantageusement, on construit une matrice sous-structurée (K) en 20 séparant les composantes : - représentant le comportement vibratoire des modes internes (KiR) ; représentant les noeuds (NJ) d'interfaces communs avec le pneumatique (KNJNJR ) ; et représentant les noeuds d'interfaces externes différents des noeuds communs avec le pneumatique (KEER) ; 25 pour obtenir : -9- -R ùR ùR K 11 K INJ K JE ùR ùR ùR ùR K = K Nui K NJNJ K NJE ùR ùR ùR KE! KENJ K EE [0046] Le modèle résultant KR décrit le comportement vibratoire de la roue et est compatible avec un couplage par les interfaces avec un pneumatique. [0047] Pour pouvoir coupler le modèle de roue avec le modèle de cavité fluide, on 5 calcule également les termes (DR et AR qui permettent de reconstituer la géométrie du système initial de la roue selon l'équation : uRC = Ra/R+ARU ER avec uRC les degrés de liberté du maillage initial en contact avec la cavité fluide selon l'écriture rappelée précédemment. 10 [0048] Les modèles sous-structurés dynamiques du pneumatique et/ou de la roue peuvent utiliser une formulation à interface fixe dite de Craig-Bampton. Ils peuvent aussi utiliser une formulation à interface libre dite de Mac-Neal. Dans ce dernier cas, on adjoint aux matrices KP respectivement KR deux matrices (I)cR1, ACRE respectivement (DcRR, AcRR selon la théorie bien connue de l'homme de l'art. 15 [0049] Avantageusement, pour créer un modèle de l'activité vibratoire du fluide de la cavité interne : - on établit un maillage élément fini de la cavité ; - on définit les noeuds (NCP) d'interfaces physiques entre la cavité et le pneumatique ; - on définit les noeuds (NCR) d'interfaces physiques entre la cavité et la roue 20 - on résout par éléments finis les trois équations suivantes : (1) P=pco2O (2) 1 , Pn=un pw- (3) P=ù /30 avec : 25 P le champ de pression dans le fluide ; -10- 0 le champ de déplacement dans le fluide ; u le champ de déplacement de la structure ; n le vecteur normal à la surface de contact fluide / structure ; /3 le module de compressibilité du fluide ; et p la densité du fluide ; pour obtenir la liste des pulsations propres de la cavité (n?), un terme de couplage entre la cavité et le pneumatique (Ffsp) et un terme de couplage entre la cavité et la roue (FfSR). [0050] Avantageusement, en considérant dans le maillage de la cavité les noeuds (NCR) d'interfaces physiques de la roue et de la cavité et dans le maillage de la roue les noeuds (NRc) d'interfaces physiques de la roue et de la cavité, on définit des points de condensation (NJ') communs à la cavité et à la roue avec pour chacun desdits points de condensation, sa position géométrique et sa relation de condensation entre ledit point de condensation et le ou les noeuds associés (NCR) d'interfaces physiques cavité ainsi que le ou les noeuds associés (NRc) d'interfaces physiques roue. [0051] Selon un mode de réalisation préféretiel, on décrit la géométrie de la cavité et de la roue dans un nombre donné (NP) de plans méridiens ; et dans chaque plan méridien, on regroupe les noeuds de la roue en contact entre la roue et la cavité en un nombre donné de segments (NS), chaque segment étant assujetti à un déplacement de corps rigide et étant représenté par un noeud de condensation. Le nombre de plans méridiens peut être compris entre 25 et 100. Préférentiellement, le nombre de segments (NS) est compris entre 3 et 5. [0052] Avantageusement, on assemble le modèle de la cavité aux modèles des sous-structures mécaniques pneumatique et roue, et on obtient : PFP cD RFR 0 A P r fs ARrR- ùP ùP 0 KlvJ K!E 0 [ RFR, 0 ùR K 1E K NJ! K NJ! K NJNJ + K NJNJ K NJE K N'JE Ke! 0 KENJ KÉE 0 K E! KENJ 0 KÉE K PRC ù -11- [0053] Avantageusement, pour obtenir un modèle sous-structuré dynamique de ùFinal l'ensemble monté (K ) : - on considère comme noeuds externes les noeuds d'interface externes des sous-structures initiales pneumatique et roue ; on considère comme noeuds internes les coordonnées modales des sous-structures initiales pneumatique et roue, augmentées des coordonnées modales de la cavité et des noeuds communs d'interface des sous-structures initiales pneumatique et roue ; et on obtient : KùFanal = ù Final ù Final Ku KIE ù Final ù Final K El K EE [0054] Cette dernière étape permet d'augmenter l'efficacité du modèle en se ramenant à un format strictement identique à celui du modèle de pneumatique seul. Description brève des dessins [0055] Tous les détails de réalisation sont donnés dans la description qui suit, complétée 15 par les figures 1 à 6, dans lesquelles : - la figure 1 présente un organigramme synthétique de la méthode de création d'un modèle de dynamique vibratoire d'un ensemble monté ; - la figure 2 illustre l'étape 1 de définition des noeuds d'interface communs au pneumatique et à la roue ; 20 - la figure 3 illustre l'étape 4 de création d'un modèle de l'activité vibratoire du fluide de la cavité ; - la figure 4 présente un exemple d'assemblage pneumatique et roue par maillages compatibles ; la figure 5 présente un exemple d'assemblage pneumatique et roue par maillages 25 incompatibles ; et - la figure 6 présente un exemple de calcul du transfert vibratoire entre un déplacement vertical imposé dans l'aire de contact et un effort vertical relevé au centre de roue pour un modèle du pneumatique seul, un modèle pneumatique et roue et un modèle pneumatique roue et cavité. -12- Description détaillée de modes de réalisation de l'invention [0056] La figure 1 présente un organigramme général de la méthode objet de l'invention. [0057] Cette méthode comprend six étapes dans sa forme la plus complète. Ces étapes sont à réaliser successivement ou non selon le cas. Étape n°1 : Définition des interfaces communes [0058] Pour préparer le modèle de pneumatique et le modèle de roue nécessaires à la réalisation du modèle d'ensemble monté, il convient de suivre les étapes indiquées à la figure 2 respectivement pour le responsable du pneumatique et le responsable de la roue. [0059] En comparaison avec l'état de l'art, représenté en Figure 4, il n'est pas nécessaire selon l'invention que la position géométrique des points d'interfaces communs coïncident avec la position géométrique des noeuds des modèles initiaux. Ceci autorise en particulier l'utilisation de maillages incompatibles dans les calculs de roues et de pneumatiques qui peuvent de ce fait être travaillés séparément. [0060] En pratique, l'innovation permettant d'accéder à cet avantage consiste à appliquer 15 à la zone d'accrochage entre le pneumatique et la roue des techniques de condensation d'interface. [0061] La technique de condensation consiste à retenir sur le maillage initial un certain nombre de noeuds, le plus souvent mais pas forcément, les noeuds qui seront effectivement compris dans la zone de contact entre le pneumatique et la roue. La technique consiste 20 alors à imposer à l'ensemble (ou à des sous-ensemble) des noeuds retenus un mouvement d'ensemble piloté par le mouvement d'un (ou plusieurs) point(s), dit(s) de condensation. Une illustration est donnée en Figure 5. Une fonction préalablement définie, par exemple une relation linéaire, permet de relier le mouvement des noeuds géométriques de l'objet initial au mouvement du (ou des) point(s) de condensation. Selon l'invention, l'astuce 25 consiste à imposer l'identité des noeuds de condensation pour les interfaces communes entre le pneumatique et la roue. [0062] Pour fonctionner correctement, cette technique suppose qu'à l'étape 1, tel qu'indiqué à la figure 2, les responsables du pneumatique et de la roue se soient entendus -13- sur une définition commune de NJ points d'interfaces de condensation. En particulier, la définition de chacun des NJ points d'interfaces communs devra préciser sa position géométrique et la relation de condensation entre ses degrés de liberté et les degrés de liberté des noeuds physiques du maillage initial. [0063] En respectant cette procédure de définition des interfaces, la position géométrique des noeuds des maillages du pneumatique et de la roue peuvent être différents mais les modèles créés aux étapes 2 et 3 seront automatiquement compatibles lors de l'assemblage de l'étape 5. [0064] Pour utiliser ces techniques de condensations l'astuce consiste, avantageusement, à condenser sur un même point de condensation le maximum de noeuds physiques possible sans dégrader la représentation du comportement vibratoire de la structure dans le domaine fréquentiel requis. Dans le cas du couplage entre le pneumatique et la roue, préférentiellement : (1) on décrit la géométrie selon NP plans contenant l'axe de rotation du pneumatique et de la roue. Ces plans sont dits plans méridiens. Il faut choisir NP de façon à bien représenter les ordres de déformées circonférentielles du pneumatique et de la roue dans la bande fréquentielle qui nous intéresse, en pratique 0-500 Hz. NP peut être compris entre 25 et 100. Classiquement, à partir de 30 plans méridiens, l'ordre 3 de la déformée circonférentielle du crochet de jante est bien représenté, ce qui couvre usuellement la gamme 0-500 Hz ; et (2) dans chaque plan méridien et pour chaque zone de contact pneumatique et roue, on assujettit tous les noeuds du contact pneumatique et roue à un déplacement de corps rigide comme illustré sur la figure 5. Un tel choix de condensation autorise une représentation fidèle de la structure couplée jusqu'à la bande de fréquence où des modes locaux de déformations des crochets de jantes apparaissent, c'est-à-dire pour un ensemble-monté tourisme usuel, au-delà de 500 Hz. [0065] A la figure 4, les maillages de la roue et du pneumatique sont compatibles, c'est-30 à-dire que à chaque noeud d'interface physique du maillage de la roue correspond un noeud - 14 - physique du maillage du pneumatique. Les deux maillages doivent être établis avec une forte concertation entre les deux responsables de la roue et du pneumatique. [0066] En revanche, à la figure 5, on assujettit tous les noeuds des deux interfaces physiques gauche et droite à un déplacement de corps rigide. Tous ces noeuds sont 5 condensés sur deux noeuds virtuels disposés près du centre roue, un à gauche et un à droite. [0067] Pour ce cas particulier, il y aura donc 2xNP noeuds d'interface pour chaque structure, ce qui est notablement moins que le nombre de noeuds géométriques des maillage initiaux, la réduction allant usuellement d'un facteur 5 à 20. La seule contrainte de calcul pour chacune des parties sera d'avoir le même découpage sectoriel de sa structure et de 10 déclarer les noeuds de condensation aux mêmes coordonnées. Étape n°2 : Construction du modèle de pneumatique : [0068] Le modèle dynamique du pneumatique est alors construit selon l'état de l'art rappelé plus haut de la sous structuration dynamique soit à interface fixe, également désignée par Craig-Bampton, soit à interface libre, également désignée par Mac-Neal. 15 Selon l'invention, par écart au modèle usuel de pneumatique seul, décrit dans l'état de l'art, on construit alors la matrice sous-structurée de telle sorte à bien séparer les composantes représentant respectivement le comportement vibratoire des modes internes, des NJ noeuds d'interface commun avec la roue définis lors de l'étape 1 et les noeuds d'interface externe différents des noeuds communs avec la roue. 20 [0069] Le résultat de cette étape consiste en trois matrices K, C et M décrivant le comportement vibratoire du pneumatique et compatibles avec un couplage par les interfaces avec une roue. Dans le cas d'un calcul de type Mac-Neal, le modèle comportera une matrice supplémentaire pour décrire les déformées modales des interfaces. Par construction, le modèle prend par exemple la forme 25 -P -P -P -P K = KK/NJ -P K/E -P K NJNJ -P K NJ/ K NJE -P -P -P KEl K ENJ K EE 20 -15- [0070] Classiquement, on adjoint à ce modèle les deux matrices JP et AP qui permettent de reconstituer la géométrie du système initial selon l'équation définie plus haut dans l'état de l'art. L'astuce est ici de calculer ces grandeurs sur la frontière entre le pneumatique et la cavité fluide pour permettre ensuite le couplage entre modes de structure et mode de fluide. [0071] L'exposant P désigne ici les grandeurs se rapportant à la structure pneumatique. Étape n°3 : Construction du modèle de roue : [0072] Le modèle dynamique de la roue est également construit selon l'état de l'art rappelé plus haut de la sous structuration dynamique soit à interface fixe, également désignée par Craig-Bampton, soit à interface libre, également désignée par Mac-Neal. Selon l'invention, par écart au modèle usuel de pneumatique seul, décrit dans l'état de l'art, on construit alors la matrice sous-structurée de telle sorte à bien séparer les composantes représentant respectivement le comportement vibratoire des modes internes, des NJ noeuds d'interface commun avec le pneumatique définis à l'étape 1 et les noeuds d'interface externe différents des noeuds communs avec le pneumatique. [0073] Le résultat de cette étape consiste en trois matrices K, C et M permettant de décrire le comportement vibratoire de la roue et compatible avec un couplage par les interfaces avec un pneumatique. Dans le cas d'un calcul de type Mac-Neal, le modèle comportera une matrice supplémentaire pour décrire les déformées modales des interfaces. Par construction, le modèle prend par exemple la forme : -R -R -R Ku KlNJ KJE R K = K Ni! K NJNJ K NJE -R -R -R KEr KENJ K EE [0074] Classiquement, on adjoint à ce modèle les deux matrices (DR et AR qui permettent de reconstituer la géométrie du système initial selon l'équation définie plus haut dans l'état de l'art. L'astuce est ici de calculer ces grandeurs sur la frontière entre la roue et la cavité fluide pour permettre ensuite le couplage entre modes de structure et mode de fluide. 25 [0075] L'exposant R désigne ici les grandeurs se rapportant à la structure roue. - 16- Étape n°4 : Construction du modèle de cavité : [0076] La construction de modèle de cavité reprend les principes de formulation décrits dans l'état de l'art en référence [Clayton & All]. [00771 Comme indiqué dans l'article, la démarche de calcul repose sur 3 équations 5 fondamentales : (1) L'équation du mouvement dans le fluide : P=pco2O ; (2) La continuité du mouvement à l'interface fluide/structure : (3) L'équation constitutive du fluide : P=ù /30 . [0078] Cette dernière équation rend compte du fait qu'en acoustique, les variations de 10 pression sont suffisamment rapides pour être considérées comme adiabatiques. [0079] Dans ces équations, les variables sont P le champ de pression dans le fluide 0 le champ de déplacement dans le fluide u le champ de déplacement de la structure 15 n le vecteur normal à la surface de contact fluide / structure /3 le module de compressibilité (Bulk modulus) du fluide p la densité du fluide. [0080] En combinant ces trois équations selon les règles décrites dans l'article en référence, l'homme de l'art peut calculer par éléments finis les grandeurs nécessaires pour 20 prendre en compte l'interaction entre l'activité vibratoire de la cavité fluide et l'activité vibratoire des structures mécaniques à sa frontière. En particulier, il accède à deux grandeurs utilisées par la suite Ç242 les pulsations propres de la pression dans le fluide, et Ffs terme descriptif du couplage entre les modes du fluide et la structure mécanique en contact. Ces deux grandeurs sont parfaitement connues lorsque la géométrie de la cavité est définie 25 et que les caractéristiques du fluide, module de compressibilité (Bulk modulus) /3 et densité p, sont choisies. Leurs règles d'utilisation sont décrites dans l'état de l'art. [0081] En nouveauté par rapport à cet état de l'art, l'invention consiste à séparer l'interface entre le fluide et la structure en plusieurs interfaces correspondant d'une part 1 Pii=un ; et p0)2 -17- aux noeuds en contact entre le fluide et la structure du pneumatique et d'autre part entre les noeuds en contact entre le fluide et la structure de la roue comme décrit à la Figure 3. [0082] De ce fait, selon l'invention, par comparaison avec les deux grandeurs initiales 12? et Ffs, le modèle de la cavité contient finalement trois grandeurs : d'une part il conserve la définition de la liste Q? des pulsations propres de la pression dans le fluide, d'autre part, le terme Ffs est séparé en 2 grandeurs FfsP et FfSR correspondant respectivement aux termes de couplage du fluide avec le pneumatique et la roue. [0083] Toujours en nouveauté par rapport à cet état de l'art, l'invention propose également d'autoriser une incompatibilité entre le maillage à l'interface du fluide et de la structure en retenant NC noeuds physiques sur le maillage du fluide à l'interface entre le fluide et la structure. Ces noeuds sont répartis en NCP à l'interface avec le pneumatique et NCR à l'interface avec la roue. Pour la structure, on peut alors définir une relation linéaire entre le déplacement de NCP (resp. NCR) noeuds d'interfaces et les déplacements des P (resp. R) noeuds physiques du pneumatique (resp. de la roue). Cette relation peut par exemple s'écrire : uNCR = DRuR et uNCP = DRuR [0084] La formulation du couplage évolue alors simplement en modifiant les termes de couplage Ffs1' et FfSR qui deviennent respectivement FffPDP et FfSRDR. Par souci de simplification, on gardera pour la suite une écriture sous la forme générale FfsP et FfSR qui couvrira par convention les cas de maillages compatibles ou incompatibles. [0085] Pour utiliser cette technique de condensation, avantageusement, on condense sur un même point de condensation le maximum de noeuds physiques possible sans dégrader la représentation du comportement vibratoire de la structure et du fluide. Sur la gamme [0-500] Hz, le fluide présentera au maximum des modes d'ordre 2 dans le cas usuel d'un pneumatique tourisme gonflé à l'air. En revanche, sur cette gamme de fréquences, le pneumatique présente une densité modale importante et imposera que le nombre de noeuds de condensation du couplage entre le pneumatique et la cavité soit proche du nombre de noeuds du maillage initial. - 18- [0086] En revanche, la densité modale pour la roue est bien plus faible, de l'ordre d'une dizaine de modes entre 0 et 500 Hz. Il devient alors astucieux de condenser les interfaces de couplage entre la roue et la cavité. Pour optimiser au mieux le ratio entre taille des modèles et représentation physique, de préférence : (1) on décrit la géométrie selon NP plans contenant l'axe de rotation de la roue et de la cavité. Ces plans sont dits plans méridiens. Il faut choisir NP de façon à bien représenter les ordres de déformées circonférentielles de la roue et de la cavité dans la bande fréquentielle qui nous intéresse, en pratique 0-500 Hz. Le nombre NP est usuellement compris entre 25 et 100. Classiquement, à partir de 30 plans méridiens, l'ordre 3 de la déformée circonférentielle de la roue et de la cavité est bien représenté, ce qui couvre correctement la gamme 0-500 Hz. (2) dans chaque plan méridien, on regroupe les noeuds de la jante de la roue en contact entre la roue et la cavité en un nombre donné NS de segments, chaque segment étant assujetti à un déplacement de corps rigide et étant représenté par un noeud de condensation. Le nombre de segments peut usuellement être compris entre 3 et 5. Un tel choix de condensation autorise une représentation fidèle de la structure couplée jusqu'à la bande de fréquence où des modes locaux de déformations de la jante de la roue apparaissent, c'est-à-dire pour un ensemble-monté tourisme usuel, au-delà de 500 Hz. Les modes de cavité intéressants dans cette bande de fréquence sont eux aussi bien représentés. [0087] Pour ce cas particulier, il y aura donc NPxNS noeuds d'interface entre la cavité et la roue, ce qui est notablement moins que le nombre de noeuds géométriques des maillage initiaux, la réduction allant usuellement d'un facteur 5 à 20. La seule contrainte de calcul pour chacune des parties sera d'avoir le même découpage sectoriel de sa structure et de 25 déclarer les noeuds de condensation aux mêmes coordonnées. Étape n°5 : Assemblage des modèles pneumatique roue et cavité [0088] Selon l'invention, avec la nouveauté de séparer l'interface fluide/structure en sous-interfaces correspondant respectivement au pneumatique et à la roue dans les équations du couplage entre fluide et structure, le modèle assemblé 30 pneumatique/roue/cavité s'écrit comme suit : -19- K PRC = [~PrP]1 IS [(D R TA1 Rs1 0 [APFp' [A RF PrP ARrsR ,,. ~. o As. , K O K/NJ K E 0 0 K, KINJ 0 KE ùP ùR ùP ùR ùP ùR K NJ/ K NJ/ K NJNJ + K NJNJ K NJE K 1JE KEr 0 KENJ KP 0 O K E/ K ENJ O K EE (DPFP (DRFR [0089] Ce modèle est constitué à partir de l'ensemble des grandeurs calculées lors des étapes précédentes. L'avantage d'assembler pneumatique, roue et cavité est illustré sur les courbes de transfert en Figure 6. [0090] La Figure 6 présente un exemple de calcul du transfert vibratoire entre un déplacement vertical imposé dans l'aire de contact et un effort vertical relevé au centre de roue avec et sans présence de l'activité vibratoire de la roue, et de la cavité. La figure a correspond au modèle de pneumatique seul, la courbe b au modèle pneumatique et roue et la courbe c au modèle complet, pneumatique-roue-cavité. On constate que dans la plage de fréquences 0-190 Hz environ, les trois modèles donnent des résultats proches. En revanche, la prise en compte de la roue est sensible notamment entre 200 et 260 Hz, et seul le modèle complet associant les activités vibratoires du pneumatique, de la roue et de la cavité permet de mettre en évidence un pic de transfert situé vers 220 Hz. [0091] Un tel modèle est toutefois notoirement plus volumineux que le modèle initial de 15 pneumatique, notamment du fait de l'existence des degrés de liberté communs au pneumatique et à la roue. Étape n°6 : Sous structuration [0092] Pour accéder à l'avantage d'utilisation du modèle pneumatique seul, une opération de sous structuration peut alors être menée sur la base des matrices obtenues à 20 l'étape précédente, en séparant les degrés de liberté interne du système sur lesquelles ne s'appliquera pas de forces externe, des degrés de liberté externes qui pourront être soumis à des efforts externes. Une telle technique est en particulier décrite dans [Clayton & All] pour le cas pneumatique et cavité. 2934907 - 20 - [0093] Par rapport à cet état de l'art la nouveauté consiste à retenir la partition suivante : les noeuds externes sont les noeuds d'interface externes des sous structures initiales pneumatiques et roue et les noeuds internes sont les coordonnées modales des sous structures initiales pneumatique et roue, augmentés des coordonnées modale de la cavité 5 fluide et augmentés des noeuds communs d'interface des sous structures initiales pneumatique et roue. [0094] Selon l'invention, avec ces définitions d'interface, les techniques usuelles de sous-structuration peuvent alors être pratiquées car les noeuds d'interface commun au pneu et à la roue ne seront pas soumis à des efforts externes. Les techniques de sous-structuration 10 conduisent à la création de trois matrices K, M et C, éventuellement augmenté de la matrice des déformées modales des interfaces restantes dans le cas d'une sous-structuration de type Mac-Neal, de format strictement analogue à celles représentant l'activité vibratoire du pneumatique seul, et donc présentant les avantages d'utilisation de ces dernières, mais intégrant fidèlement le comportement vibratoire du pneumatique, de la roue et du fluide de la cavité interne. [0095] Le modèle final reprend donc la forme du modèle de pneumatique seul KùFinal = ù Final ù Final Ku KIE ù Final ù Final KEI KEE [0096] La taille d'un tel modèle est alors le nombre de modes retenus par l'utilisateur, typiquement entre 50 et 100 pour représenter la dynamique de l'ensemble monté de 0 à 500 Hz, augmenté du nombre de degrés de liberté des sous-structures pneumatique et roue, typiquement quelques dizaines de degrés de liberté dans l'aire de contact et au centre de la roue. [0097] De façon évidente, il est possible d'omettre le modèle de roue ou le modèle de cavité dans la démarche pour créer un modèle partiel de l'activité vibratoire au centre-roue. [0098] L'invention n'est pas limitée aux exemples décrits et représentés et diverses modifications peuvent y être apportées sans sortir de son cadre défini par les revendications annexées. REVENDICATIONS1. Procédé de création d'un modèle de prévision de la dynamique vibratoire d'un ensemble monté constitué d'un pneumatique monté sur une roue et comportant un fluide dans la cavité interne définie par le pneumatique et la roue ou � cavité �, comportant les étapes suivantes : (1) établir un maillage éléments finis du pneumatique, établir un maillage éléments finis de la roue et définir des noeuds d'interface communs au pneumatique et à la roue ; (2) créer un modèle sous-structuré dynamique du pneumatique condensé aux 1 o noeuds d'interface communs (K P ) ; (3) créer un modèle sous-structuré dynamique de la roue condensé aux noeuds d'interface communs (KR) ; (4) créer un modèle de l'activité vibratoire du fluide de la cavité (SVP, rfSP, FfSR) ; et (5) assembler le modèle de la cavité aux modèles des sous-structures mécaniques 15 pneumatique et roue (KRRC ) 2. Procédé selon la revendication 1, comportant comme dernière étape, une étape de sous-structuration du modèle assemblé (KPRC) pour obtenir un modèle sous-structuré ù dynamique de l'ensemble monté (KFinal ). 20 3. Procédé selon l'une quelconque des revendications 1 et 2, dans laquelle, pour définir des nœuds d'interface communs au pneumatique et à la roue : on définit dans le maillage du pneumatique les noeuds (X) d'interfaces physiques du pneumatique et de la roue ; 25 - on définit dans le maillage de la roue les noeuds (Y) d'interfaces physiques de la roue et du pneumatique ; et on définit des points de condensation (NJ) communs au pneumatique et à la roue avec pour chacun desdits points de condensation, sa position géométrique et sa relation de condensation entre ledit point de condensation et le ou les noeuds associés 30 (X) d'interfaces physiques pneumatique ainsi que le ou les noeuds associés (Y) d'interfaces physiques roue. 2934907 - 22 - 4. Procédé selon la revendication 3, dans laquelle : on décrit la géométrie du pneumatique et de la roue dans un nombre donné (NP) de plans méridiens ; et 5 dans chaque plan méridien et pour chaque zone de contact pneumatique et roue, on assujettit tous les noeuds du contact pneumatique et roue à un déplacement de corps rigide. 5. Procédé selon la revendication 4, dans laquelle le nombre de plans méridiens est 10 compris entre 25 et 100. 6. Procédé selon l'une des revendications 1 à 5, dans laquelle, pour créer un modèle sous-structuré dynamique du pneumatique condensé aux noeuds d'interface communs, on construit une matrice sous-structurée (KP) en séparant les composantes : 15 - représentant le comportement vibratoire des modes internes (Kit") ; représentant les noeuds (NJ) d'interfaces communs avec la roue (KNJNJP ) ; et représentant les noeuds d'interfaces externes différents des noeuds communs avec la - P roue (Kcc ) ; pour obtenir : KP KN] KIR -P -P -P -P 20 K = KN/1 KNJNJ KNIE -P -P -P Ker KEN/ KLE 7. Procédé selon l'une quelconque des revendications précédentes, dans laquelle, pour créer un modèle sous-structuré dynamique de la roue condensé aux noeuds d'interface communs, on construit la matrice sous-structurée (KR) en séparant les composantes : 25 - représentant le comportement vibratoire des modes internes ( K h i ? ) ; représentant les noeuds (NJ) d'interfaces communs avec la roue (KNJNJR ) ; et représentant les noeuds d'interfaces externes différents des noeuds communs avec la roue (KEER) ;-23-pour obtenir : -R ùR ùR K11 KJNJ K1E ùR ùR ùR ùR K = K NJI K NJNJ K NJE ùR ùR ùR KEI K ENJ K EE 8. Procédé selon l'une des revendications 6 et 7, dans laquelle le modèle sous-structuré 5 dynamique utilise une formulation à interface fixe dite de Craig-Bampton. 9. Procédé selon l'une des revendications 6 et 7, dans laquelle le modèle sous-structuré dynamique utilise un formulation à interface libre dite de Mac-Neal et dans laquelle on adjoint aux matrices (KP ; KR) deux matrices (OCRP , ACRP ; (DCRR , AcRR) pour 1 0 permettre de reconstituer la géométrie du système initial. 10. Procédé selon l'une quelconque des revendications précédentes, dans laquelle, pour créer un modèle de l'activité vibratoire du fluide de la cavité interne : on établit un maillage élément fini de la cavité ; 15 - on définit les noeuds (NCP) d'interfaces physiques entre la cavité et le pneumatique ; on définit les noeuds (NCR) d'interfaces physiques entre la cavité et la roue ; on résout par éléments finis les trois équations suivantes : (1) P=pv2e (2) 1 2 Pn=un pw 20 (3) P=û X39 avec : P le champ de pression dans le fluide ; B le champ de déplacement dans le fluide ; u le champ de déplacement de la structure ; 25 h le vecteur normal à la surface de contact fluide / structure ; fi le module de compressibilité (Bulk modulus) du fluide ; et p la densité du fluide ; 2934907 - 24 - pour obtenir la liste des pulsations propres de la cavité (û?), un terme de couplage entre la cavité et le pneumatique (I'f?) et un terme de couplage entre la cavité et la roue (FfsR). 11. Procédé selon la revendication 10, dans laquelle, en considérant dans le maillage de la cavité les noeuds (NCR) d'interfaces physiques de la roue et de la cavité et dans le maillage de la roue les noeuds (NRc) d'interfaces physiques de la roue et de la cavité, on définit des points de condensation (NJ') communs à la cavité et à la roue avec pour chacun desdits points de condensation, sa position géométrique et sa relation de condensation entre ledit point de condensation et le ou les noeuds associés (NCR) d'interfaces physiques cavité ainsi que le ou les noeuds associés (NRc) d'interfaces physiques roue. 12. Procédé selon la revendication 11, dans laquelle : - on décrit la géométrie de la cavité et de la roue dans un nombre donné (NP) de plans méridiens ; et - dans chaque plan méridien, on regroupe les noeuds de la roue en contact entre la roue et la cavité en un nombre donné de segments (NS), chaque segment étant assujetti à un déplacement de corps rigide et étant représenté par un noeud de condensation. 13. Procédé selon la revendication 12, dans laquelle le nombre de plans méridiens est 20 compris entre 25 et 100. 14. Procédé selon l'une des revendications 12 et 13, dans laquelle le nombre de segments (NS) est compris entre 3 et 5. 25 15. Procédé selon l'une quelconque des revendications précédentes, dans laquelle, on assemble le modèle de la cavité aux modèles des sous-structures mécaniques pneumatique et roue, et on obtient :-25- Q~ f Prfs ~Rrfs 0 [API ]' K 1 0 K ENI K P O [ARrR] 0 KE1 KENJ 0 KEE 16. Procédé selon l'une quelconque des revendications 2 à 15, dans laquelle pour obtenir un modèle sous-structuré dynamique de l'ensemble monté (KFina} on considère comme noeuds externes les noeuds d'interface externes des sous-structures initiales pneumatique et roue ; on considère comme noeuds internes les coordonnées modales des sous-structures initiales pneumatique et roue, augmentées des coordonnées modale de la cavité et des noeuds communs d'interface des sous-structures initiales pneumatique et roue ; lo on réalise une sous-structuration dynamique ; et- on obtient : APFP ARr.R ù fs fs {cDPrfs K, 0 KINJ K e 0 KPRC = L~Rrfs 0 K KINJ 0 Km 0 K NJI K NJI K NJNJ + K NJNJ K NJE K NJE K Fanal = --Final ù Final ù K11 KrE ùFinal ùFinal KEI KEE