FR 2890905 A1 20070323 FR 0509716 A 20050922 2890905 1 L'invention concerne un dispositif de commande d'une suspension d'un véhicule automobile. Un domaine d'application de l'invention concerne les véhicules automobiles ayant une suspension à ressort, une suspension 5 hydropneumatique ou une suspension d'un autre type. Ces suspensions comportent sur chaque roue un amortisseur à loi d'amortissement variable, pouvant être réglée par un actionneur commandé par un calculateur embarqué sur le véhicule. Le calculateur reçoit en entrée des mesures fournies par des 10 capteurs et calcule à partir de celles-ci la ou les grandeurs de commande des actionneurs des amortisseurs. Le calculateur tient compte notamment des accélérations subies par la caisse du véhicule au cours du trajet, telles que par exemple l'accélération modale de pompage dans la direction verticale, l'accélération modale de roulis autour d'un axe longitudinal et l'accélération modale de tangage autour d'un axe transversal. Par intégration, le calculateur calcule les vitesses modales correspondantes de la caisse. On connaît ainsi des dispositifs dans lesquels le calculateur met en oeuvre une régulation pour faire tendre vers zéro la vitesse modale verticale en pompage, la vitesse modale angulaire en roulis et la vitesse modale angulaire en tangage, cette logique étant communément appelée "Skyhook" et devant améliorer le confort des personnes dans la voiture. L'un des problèmes de ces dispositifs est qu'ils ne garantissent pas toujours l'adhérence du véhicule à la route. En effet, dans certaines circonstances, comme les fortes sollicitations, les suspensions commandées par une logique de type Skyhook jouent au détriment de l'effort vertical sur le pneu. Par exemple, sur une bosse, la logique Skyhook commandera d'abord une souplesse de l'amortisseur pour favoriser le confort au détriment de l'effort vertical. Egalement, les sollicitations imposées par le conducteur du véhicule entraînent pas ou peu de remontées d'efforts venus du sol vers la caisse du véhicule, ce qui est ressenti comme peu rassurant ou donne l'impression au conducteur que ses actions n'ont pas de prise sur le véhicule. Par ailleurs, les vitesses de caisse sont considérablement dégradées sous de fortes sollicitations du conducteur, notamment lorsqu'elles sont estimées à l'aide d'accéléromètres. En effet, les accéléromètres verticaux présentent le gros défaut d'être sensibles à leur inclinaison en tangage ou en roulis du fait de la projection de la gravité, mais surtout de mesurer verticalement environ le dixième des sollicitations inertielles longitudinales ou latérale. L'estimation des vitesses de caisse sous de fortes sollicitations du conducteur en est rendue médiocre. Les vitesses de caisse étant médiocres, la logique SkyHook ne peut être convenablement utilisée. L'invention vise à obtenir un dispositif de commande d'une suspension, qui pallie les inconvénients de l'état de la technique et dispose d'une logique autre que la logique Skyhook et faisant suivre à la caisse le profil de la route pour privilégier la tenue de caisse. En particulier, la logique utilisée dans le dispositif de commande suivant l'invention doit traiter au moins un mode de caisse, que sont le pompage, le roulis et le tangage. A cet effet, un premier objet de l'invention est un dispositif de commande d'une suspension d'une caisse de véhicule automobile sur ses roues, comportant un calculateur apte à calculer une grandeur de commande d'un actionneur d'au moins un amortisseur variable de la suspension en fonction d'au moins un effort modal de consigne de la caisse, caractérisé en ce que le calculateur comporte: - un premier moyen de calcul d'au moins une vitesse modale relative de caisse par rapport au plan moyen des roues, 2890905 3 - un deuxième moyen de calcul de l'effort modal de consigne de l'amortisseur, en fonction de la vitesse modale relative de caisse par rapport au plan moyen des roues. Suivant l'invention, le dispositif de commande met en oeuvre une logique de type Roadhook, c'est-à-dire qui suit le profil de la route, cette logique étant également appelée logique de tenue de caisse ou logique de comportement. Le principe de cette logique de tenue de caisse est de faire tendre vers zéro ou de minimiser l'une ou plusieurs des accélérations modales de la caisse par rapport au plan des roues: accélération modale en pompage, accélération modale en roulis, accélération modale en tangage. Un deuxième objet de l'invention est un véhicule automobile comportant une caisse, des roues, une suspension de la caisse sur les roues et un dispositif de commande de la suspension tel que décrit ci-dessus. Un troisième objet de l'invention est un procédé d'obtention d'un véhicule automobile, le véhicule automobile étant muni de roues, d'une caisse, d'une suspension ayant au moins un amortisseur à amortissement variable de la caisse sur les roues, et d'un dispositif de commande de la suspension tel que décrit ci-dessus, le dispositif de commande ayant au moins un calculateur apte à calculer une grandeur de commande d'un actionneur dudit au moins un amortisseur de la suspension, le procédé d'obtention comportant une étape de montage du calculateur sur le véhicule, caractérisé en ce que le procédé d'obtention comporte au moins une étape de programmation du calculateur suivant au moins un programme comportant des instructions de programme mettant en oeuvre les moyens de calcul du dispositif de commande de la suspension. Un quatrième objet de l'invention est un programme informatique de pilotage d'un calculateur, comportant des instructions de programme pour le calcul d'au moins une vitesse modale relative de caisse par rapport au plan moyen des roues, pour le calcul de l'effort modal de consigne de l'amortisseur, en fonction de la vitesse modale relative de caisse par rapport au plan moyen des roues et pour le calcul de la grandeur de commande d'un actionneur d'au moins un amortisseur variable de la suspension en fonction d'au moins l'effort modal de consigne de la caisse, lorsqu'il est mis en oeuvre dans un dispositif de commande de suspension tel que décrit ci-dessus. L'invention sera mieux comprise à la lecture de la description qui va suivre, donnée uniquement à titre d'exemple non limitatif en référence aux dessins annexés, sur lesquels: - la figure 1 est une vue schématique en perspective d'un dispositif 15 de liaison au sol d'un essieu avant d'un véhicule, - la figure 2 est un schéma fonctionnel montrant le dispositif de commande de suspension, - la figure 3 est une vue schématique en perspective de la caisse d'un véhicule, muni de la suspension sur ses roues, - la figure 4 est un synoptique modulaire d'une unité de calcul de vitesses modales du dispositif de commande selon l'invention, - la figure 5 est un synoptique modulaire d'un estimateur prévu dans le dispositif de commande suivant la figure 4, - la figure 6 est un synoptique modulaire d'une unité de calcul 25 d'efforts modaux de type Skyhook, - la figure 7 est un synoptique modulaire d'une unité de calcul d'une masse suspendue à l'avant et à l'arrière, - la figure 8 est un organigramme du procédé de calcul des masses suspendues de l'unité selon la figure 7, - la figure 9 est un synoptique modulaire d'une unité de calcul de niveaux de mouvement et de tressautement de la caisse, - la figure 10 est un synoptique modulaire d'une unité de calcul d'efforts modaux de type Roadhook, - la figure 11 est un synoptique modulaire d'une unité de calcul de termes d'efforts modaux anticipatifs, - la figure 12 est un synoptique modulaire d'une unité de calcul d'efforts de consigne aux roues, comportant l'unité de calcul d'efforts modaux de type Skyhook et l'unité de calcul d'efforts modaux de type Roadhook, - la figure 13 représente des chronogrammes de signaux de détection de sollicitations et d'un coefficient de pondération intermédiaire calculé en fonction de ceux- ci, intervenant dans l'unité de calcul selon la figure 12, - la figure 14 montre des chronogrammes de l'angle du volant au cours d'une mise en virage simple, et d'un coefficient de pondération entre les efforts Skyhook et les efforts Roadhook, intervenant dans l'unité de calcul selon la figure 12, - la figure 15 représente des lois d'amortissement des amortisseurs variables de la suspension, - la figure 16 est un synoptique modulaire d'une unité de calcul 20 d'une loi d'amortissement de consigne en cas de détection d'une percussion, - la figure 17 d'une unité de calcul d'une loi d'amortissement de consigne en cas de détection d'une grande amplitude de mouvement de caisse, - la figure 18 est un schéma en coupe transversale montrant la liaison d'un capteur de débattement à la caisse et à une roue avant ou arrière. Aux figures 1 à 3, le véhicule 1 comporte une caisse 2 montée sur quatre roues, à savoir une roue avant gauche A, une roue avant droite B, 30 une roue arrière droite C et une roue arrière gauche D. Chaque roue A, B, C, D est reliée à la caisse 2 par son propre système S de suspension à ressort R entre des butées, mais pourrait également être une suspension hydropneumatique. Chaque système S de suspension comporte un amortisseur AM muni d'un actionneur M commandé par un calculateur CSS embarqué. Cet actionneur M est par exemple un moteur qui permet de modifier la section de passage de l'huile dans l'amortisseur AM. Il correspond donc à chaque section de passage de l'huile dans l'amortisseur, une loi d'amortissement différente de celui-ci. Ces lois d'amortissement, également appelés états d'amortissement, sont mémorisées sous la forme de courbes, de tableaux de valeurs, de formules mathématiques ou autres. La figure 15 représente ces lois ER d'amortissement, où chaque loi d'amortissement est une courbe prédéterminée de la force exercée par l'amortisseur vers la caisse en fonction de la vitesse VDEB de débattement de cet amortisseur AM, les lois de plus en plus fermes ayant, à vitesse de débattement constante, des forces plus grandes. Les états ER d'amortissement sont par exemple numérotés dans un ordre croissant pour les états d'amortissement de plus en plus fermes, c'est-à-dire correspondant à une force d'amortisseur de plus en plus grande à vitesse VDEB de débattement constante. Ainsi, un état d'amortissement minimum correspond à un état d'amortissement ayant une fermeté minimale, c'est-àdire correspondant à une force d'amortisseur supérieure ou égale à un minimum pour chaque vitesse VDEB de débattement. Le calculateur CSS est relié au réseau CAN du véhicule pour récupérer une grande partie des signaux utiles (vitesse véhicule, régulation ABS, accélérations latérale et longitudinales fournies par le système de freinage, sportivité demandée par le conducteur fournie par une interface avec l'utilisateur (boîtier de servitude intelligent), etc...). Il utilise également ses propres capteurs (liaisons filaires directes avec les capteurs) pour connaître les mouvements de la voiture à chaque instant. Il est enfin relié aux actionneurs dont il assure le pilotage. Le moteur peut être pas à pas, auquel cas l'amortisseur AM possède un nombre déterminé N de lois d'amortissement discrètes ou un moteur à courant continu asservi en position, auquel cas l'amortisseur AM possède une infinité de lois d'amortissement. Par exemple, l'actionneur à moteur pas à pas peut prendre neuf positions stables distinctes, ce qui permet d'obtenir neuf lois d'amortissement, du souple au ferme. En effet, plus la section de passage de l'huile sera petite, plus l'effort d'amortissement sera important, et plus l'amortisseur 'sera ferme. II peut exister des lois stables et des lois instables. Pour les lois stables, il suffit de piloter le moteur pas à pas pour qu'il trouve sa consigne angulaire. Une fois le pilotage terminé, l'actionneur à loi stable reste dans cette position même s'il n'est plus alimenté. A l'inverse, pour les lois instables, le moteur doit être maintenu alimenté pour rester dans cette loi. Par exemple, dans un mode de réalisation, il y a à la fois des lois stables et des lois instables, les lois instables étant par exemple positionnées entre les lois stables consécutives. Par exemple, à la figure 15, on a neuf lois stables et huit lois instables. Dans un autre mode de réalisation, toutes les lois sont stables, par exemple avec 16 lois stables. Chaque actionneur M possède une entrée COM de commande reliée au calculateur CSS pour recevoir de celui-ci une grandeur ER de commande sélectionnant une position de l'actionneur M parmi plusieurs pour imposer une loi d'amortissement prédéterminée, correspondant à cette position. Suivant l'invention, il est prévu un capteur CAP-DEB de débattement sur au moins une des roues A, B, C, D du véhicule, et de préférence sur chaque roue A, B, C, D. Chaque capteur CAP-DEB mesure donc le débattement DEB de sa roue associée par rapport à la caisse 2. Les capteurs CAP-DEB de débattement de roue sont par exemple angulaires et donnent la valeur instantanée de l'angle compris entre l'axe de rotation de la roue et la caisse 2. Par exemple, aux figures 1 et 18, chaque capteur CAP-DEB de débattement comporte une partie fixe CAPF, telle qu'un boîtier, fixée à la caisse 2, et une partie mobile CAPM, reliée à un élément fixé à la roue. Une biellette BIEL de connexion relie la partie mobile CAPM à la partie fixe CAPF et fait tourner un organe MES de mesure angulaire contenu dans la partie fixe CAPF, lorsque la roue monte ou descend par rapport à la caisse 2. La partie mobile CAPM est par exemple fixée sur un élément SUP de support de l'axe AX de rotation de la roue. Cet élément SUP de support est mobile autour d'un axe sensiblement longitudinal SUPL par rapport à la caisse 2. La partie mobile CAPM est fixée sur l'élément SUP de support à distance de son axe SUPL de rotation. Les mesures DEB de débattement des roues A, B, C, D sont envoyées des capteurs CAP-DEB au calculateur CSS, qui comporte des entrées E-DEB correspondantes. Accélérations modales Le calculateur CSS calcule à partir des mesures DEB de débattement des roues l'accélération modale zG en pompage de la caisse, l'accélération modale angulaire â. en roulis et l'accélération modale angulaire W en tangage, selon les formules ci-dessous. (FB+FC FA FD)+CBAD +Ce 8=2 I Io (e lg)(FC + FD) lg(FA + FB) + C cp= où G est le centre de gravité de la caisse 2, zG est l'altitude de G suivant une direction verticale ascendante Z, 0 est l'angle de roulis de la caisse 2 autour d'un axe longitudinal X passant par G et dirigé de l'arrière vers l'avant, cp est l'angle de tangage de la caisse 2 autour d'un axe transversal passant par G et dirigé de la droite vers la gauche, les axes X, Y, Z formant un repère orthonormal. 2G FA+FB+FC+FD M FA, FB, FC, FD sont les forces exercées par les roues respectives A, B, C, D sur la caisse 2 par l'intermédiaire de leur suspension S. v est la voie de la caisse 2, c'est-à-dire la distance entre les roues droite et les roues gauche dans le sens transversal, e est l'empattement du véhicule, 1g est la distance longitudinale entre le centre de gravité G et l'axe transversal des roues avant A et B, M est la masse prédéterminée de la caisse 2 sans occupant du véhicule. le est le moment d'inertie en roulis, et I,p est le moment d'inertie en tangage. CBAD est un couple exercé par la barre anti-devers BAD sur la caisse 2. Co est un couple de roulis, et C,p un couple de tangage. Ci-dessous sont décrits les différents moyens de calcul mettant en oeuvre le procédé de commande suivant l'invention. Le mode de calcul des accélérations modales dans le calculateur CSS est par exemple mis en oeuvre par le module 10 représenté aux figures 4 et 5. Les blocs modulaires décrits aux figures sont mis en oeuvre dans le calculateur CSS par tout moyen automatique approprié, notamment logiciel. Le module 10 comporte un premier moyen CAL de calcul des accélérations modales 2G, 6' et (fi, recevant en entrée les mesures DEB de débattement des roues. Le moyen CAL comporte: - un estimateur 11 du couple CBAD généré par la barre BAD anti-devers, - un estimateur 12 des forces FA, FB, FC, FD exercées par les roues respectives A, B, C, D sur la caisse 2, - un filtre 13 de la mesure de débattement DEB fournie en entrée sur le moyen CAL. Le filtre 13 élimine les basses fréquences de la mesure DEB de débattement fournie par les capteurs CAP-DEB. Ce filtre 13 comporte par exemple un filtre passe-haut ayant une fréquence de coupure basse supérieure ou égale à 0,2 Hz. Le filtre 13 est par exemple réalisé par un filtre passe-bande ayant donc en plus une fréquence de coupure haute, par exemple supérieure ou égale à 8 Hz, qui permet de garder une phase suffisamment constante dans la bande passante. Le débattement de roue filtré DEBF fourni à la sortie du filtre 13 à partir de la mesure DEB de débattement de roue est envoyé à l'entrée de l'estimateur 11, ainsi qu'à une autre entrée de l'estimateur 12. A partir des quatre mesures DEB(A), DEB(B), DEB(C), DEB(D) de débattement fournies par les capteurs CAP-DEB sur les roues respectives A, B, C, D, le filtre 13 fournit quatre mesures de débattement filtrées DEBF(A), DEBF(B), DEBF(C), DEBF(D). Barre anti-devers L'estimateur 11 calcule le couple CBAD de barre antidevers en fonction des valeurs de débattement filtrées DEBF fournies par le filtre 13 20 de la manière suivante: - pour la roue avant gauche: CBAD(A) = (DEBF(A) DEBF(B)).(Kbadav)/v2, - pour la roue avant droite: CBAD(B) _ -CBAD(A), - pour la roue arrière gauche: CBAD(D) = (DEBF(D) DEBF(C)).(Kbadar)/v2, - pour la roue arrière droite: CBAD(C) _ -CBAD(D), Où Kbadav est un paramètre prédéterminé correspondant à la 30 raideur de la barre anti-devers avant BAD, Kbadar est un paramètre prédéterminé correspondant à la raideur de la barre anti-devers arrière, non représentée. Effort de suspension L'estimateur 12 d'effort de suspension comporte une entrée pour les débattements filtrés DEBF, une entrée pour les débattements non filtrés DEB, une entrée pour l'état réel ER de l'actionneur, c'est-à-dire la loi ER d'amortissement qu'il met actuellement en oeuvre, cet état réel et ses changements étant par exemple mémorisés, une entrée DEAV d'effort statique sur les roues avant et une entrée DEAR d'effort statique sur les roues arrière. Ci-dessous, cet estimateur 12 est décrit à la figure 5 à titre d'exemple pour le calcul de l'effort de suspension FA sur la roue avant gauche A. Bien entendu, le calcul est analogue pour les autres efforts FB, FC, FD, en remplaçant ce qui se rapporte spécifiquement à la roue A par les valeurs correspondant à la roue B, C ou D. Dans l'estimateur 12, le débattement DEB(A) mesuré par le capteur CAP-DEB sur la roue A est envoyé à un filtre PB passe-bas limitant la bande passante du débattement DEB(A), suivi d'un module dérivateur DER pour l'obtention de la vitesse de débattement VDEB pour la roue A. Les vitesses VDEB de débattement des roues sont fournies sur une sortie de l'estimateur 12 et du module 10. Un module MFAM de calcul de la force FAM d'amortissement exercée par l'amortisseur AM sur la caisse 2 reçoit en entrée l'état réel ER et la vitesse de débattement VDEB de la roue concernée. Les lois d'amortissement des amortisseurs AM sont par exemple mémorisées à l'avance, ou peuvent être recalculées une fois l'état ER spécifié. Chacune des lois ER d'amortissement permet de calculer ou déterminer la vitesse VDEB de débattement en fonction de la force FAM d'amortissement exercée par l'amortisseur AM et inversement. Le module MFAM détermine à partir de l'état ER la loi d'amortissement actuellement en vigueur pour l'amortisseur AM de la roue A et détermine à partir de la vitesse de débattement VDEB(A) de la roue A pour cette loi sélectionnée, par exemple 2890905 12 par lecture de la courbe de cette loi, la force FAM d'amortissement de la roue A. Un autre module MFSEC de calcul d'une force FSEC de frottement sec de l'amortisseur AM de la roue A reçoit en entrée également la vitesse VDEB de débattement de celle-ci et calcule la force FSEC de frottement sec par la formule suivante: Fsec = (FsAv) . tanh(VDEB/10-) où VDEB est en cm/s, et FsAv est un coefficient de frottement sec des roues avant, ayant été calculé au préalable sur un banc d'essai et est 10 par exemple égal à environ 200 Newton. Ce coefficient de frottement est remplacé par un coefficient FsAr de frottement pour les roues arrière. Estimateur de caractéristiques statiques Un module MAS de calcul d'assiette statique AS reçoit en entrée les débattements DEB des quatre roues A, B, C, D et calcule à partir de ceux-ci l'assiette statique AS, qui représente la position statique d'équilibre de la suspension S lorsque le véhicule est immobile sur un sol horizontal. Ce module MAS calcule une assiette statique avant ASav et une assiette statique arrière ASar. L'assiette statique avant ASav est par exemple calculée comme étant le débattement moyen DEBAVMOY (demi-somme) des débattements DEB des roues avant A, B, filtré par un filtre passe-bas, par exemple de type Butterworth du second ordre, débattement moyen filtré auquel est ensuite rajoutée une constante de décalage d'assiette avant. L'assiette statique arrière ASar est par exemple calculée comme étant le débattement moyen DEBARMOY (demi-somme) des débattements DEB des roues arrières C, D, filtré par un filtre passe-bas, par exemple de type Butterworth du second ordre, débattement moyen filtré auquel est ensuite rajoutée une constante de décalage d'assiette arrière. On suppose que le capteur de débattement CAP-DEB est calibré pour mesurer le débattement par rapport à cette assiette statique AS. Un additionneur AD1 additionne le débattement filtré DEBF-A de la roue A à l'assiette statique AS calculée pour la roue A, c'est-à-dire l'assiette statique avant, pour obtenir la longueur réelle LR du ressort R associée à la roue A. Le module MAS de calcul d'assiette statique AS fait par exemple partie d'un estimateur 20 de caractéristiques statiques représenté à la figure 6, recevant en entrée les débattements DEB des quatre roues A, B, C, D, une pression statique avant et une pression statique arrière dans le cas d'une suspension hydropneumatique, la vitesse VVH du véhicule, une information d'ouvrant 10. La vitesse WH du véhicule est par exemple fournie par un capteur de vitesse ou tout autre moyen de calcul. L'estimateur 20 de caractéristiques statiques comporte: - un moyen de calcul d'une masse dynamique apparente avant MDAAV et d'une masse dynamique apparente arrière MDAAR, en fonction des débattements DEB, - un moyen de calcul d'un biais aérodynamique avant BAAV et d'un 15 biais aérodynamique arrière BAAR, en fonction de la vitesse du véhicule VVH, un moyen de calcul de la masse suspendue MSUS du véhicule et d'une valeur de répartition de masse RMAvAr entre l'avant et l'arrière du véhicule, en fonction de la masse dynamique apparente avant MDAAV, de la masse dynamique apparente arrière MDAAR, du biais aérodynamique avant BAAV et du biais aérodynamique arrière BAAR, - un moyen de calcul du moment d'inertie en roulis le et du moment d'inertie en tangage I en fonction de la masse suspendue MSUS et de la masse suspendue arrière MSUSAR, - un moyen de calcul de la longueur Ig séparant le centre de gravité G de l'axe des roues avant A, B, - un moyen de calcul d'une raideur modale en pompage kZ, d'une raideur modale en tangage kp et d'une raideur modale en roulis ke, en fonction de l'assiette statique AS et de la valeur de répartition de masse RMAvAr entre l'avant et l'arrière. La masse dynamique apparente avant MDAAV est calculée en -calculant le débattement relatif avant, égal au débattement moyen (demi-somme) des débattements DEB des roues avant A, B, auquel est ensuite rajoutée une constante de décalage avant, - extrayant un effort dynamique avant de flexion de ressort EDFAV à 5 partir d'une table ou courbe enregistrée donnant cet effort EDFAV en fonction du débattement relatif avant, calculant la masse dynamique apparente avant MDAAV par la formule: MDAAV = (EDFAV.2/g) + constante avant, Où g est la constante d'accélération de la pesanteur = 9,81 m.s-2. La masse dynamique apparente arrière MDAAR est calculée en -calculant le débattement relatif arrière, égal au débattement moyen (demi-somme) des débattements DEB des roues arrière C, D, auquel est ensuite rajoutée une constante de décalage arrière - extrayant un effort dynamique arrière de flexion de ressort EDFAR à partir d'une table ou courbe enregistrée donnant cet effort EDFAR en fonction du débattement relatif arrière, calculant la masse dynamique apparente arrière MDAAR par la formule: MDAAR = (EDFAR.2/g) + constante arrière. L'effort dynamique de flexion de ressort est nul dans la position d'équilibre du ressort correspondant à sa position statique, le débattement relatif avant étant le débattement par rapport à la position d'équilibre statique, l'extraction se faisant par exemple par interpolation de la table mais pouvant également être effectuée à partir d'une courbe enregistrée de EDFAV, EDFAR. Dans le cas d'une suspension hydropneumatique, la masse MDAAR et la masse MDAAV sont calculées en utilisant la pression statique avant et la pression statique arrière. Le biais aérodynamique avant BAAV, homogène à une masse en kg, est calculé par la formule: BAAV = (CAV.VVH2)/g, Où CAV est un coefficient prédéterminé aérodynamique avant. Le biais aérodynamique arrière BAAR, homogène à une masse en kg, est calculé par la formule: BAAR = (CAR.VVH2)/g, Où CAR est un coefficient prédéterminé aérodynamique arrière. Calcul de la masse suspendue MSUS du véhicule et de la valeur de répartition de masse RMAvAr On calcule d'abord une masse suspendue d'essieu avant MSUSEAV. Pour ce faire, ainsi que cela est représenté aux figures 7 et 8, on filtre lors de l'étape SI la somme (masse dynamique apparente avant MDAAV + biais aérodynamique avant BAAV) par un filtre passe-bas PB1 pour obtenir une masse suspendue d'essieu avant filtrée MSUSEAVF. Puis, on examine - à l'étape S2, si la vitesse VVH du véhicule est comprise entre un seuil bas prédéterminé WH1 et un seuil haut prédéterminé VVH2, - à l'étape S3, si l'information d'ouvrant IO est à fermé ou la vitesse VVH du véhicule est supérieure à un seuil prescrit WH3, - à l'étape S4, si l'écart entre la masse suspendue d'essieu avant filtrée MSUSEAVF(n) et sa valeur MSUSEAVF(n-1) précédemment enregistrée en mémoire est suffisant (supérieur en valeur absolue à un écart prescrit A). Dans le cas où ces conditions sont réalisées, la masse suspendue d'essieu avant MSUSEAV est prise égale à la masse suspendue d'essieu avant filtrée MSUSEAVF et est enregistrée dans la mémoire MEM, à l'étape S5 et dans la position du commutateur logique COMLOG représentée à la figure 7. Dans le cas où l'une, plusieurs de ou toutes ces conditions ne sont pas réalisées, la masse suspendue d'essieu avant MSUSEAV(n) est inchangée et reste égale à la valeur MSUSEAV(n-1) précédemment enregistrée dans la mémoire MEM, à l'étape S6 et dans l'autre position du commutateur logique COMLOG. Puis, à l'étape S7, on calcule une masse suspendue avant MSUSAV en filtrant la masse suspendue d'essieu avant MSUSEAV par un filtre passebas PB2, et éventuellement en saturant les valeurs obtenues par ce filtrage au-dessus d'un seuil haut et au-dessous d'un seuil bas. Les filtres passe-bas PB1 et PB2 sont par exemple d'ordre 1 avec chacun une fréquence de coupure à 0,02 Hz. Le déroulement est analogue pour le calcul de la masse suspendue d'essieuarrière MSUSEAR et de la masse suspendue arrière MSUSAR, en remplaçant MDAAV+BAAV par MDAAR+BAAR et en remplaçant MSUSEAVF par MSUSEARF. La masse suspendue du véhicule MSUS est alors calculée en faisant la somme de la masse suspendue avant MSUSAV et de la masse 15 suspendue arrière MSUSAR MSUS = MSUSAV + MSUSAR La valeur de répartition de masse avant arrière RMAvAr est alors calculée en divisant la masse suspendue avant MSUSAV par la masse suspendue du véhicule MSUS RMAvAr = MSUSAV / MSUS Calcul des moments d'inertie Le moment d'inertie en roulis le est calculé en fonction de la masse suspendue arrière MSUSAR par la formule le = Ay.MSUSAR + By Avec MSUSAR = (1-RMAvAr).MSUS, Où Ay et By sont des paramètres prédéterminés. Le moment d'inertie en tangage I9 est calculé en fonction de la masse suspendue MSUS par la formule 19 = AX.MSUS + BX Où AX et B, sont des paramètres prédéterminés. Calcul de la longueur Ig et des raideurs modales On calcule une raideur de suspension avant kAV et une raideur de suspension arrière kAR. La raideur de suspension avant kAV est obtenue en extrayant de la table ou courbe préenregistrée donnant la raideur de la suspension avant en fonction de l'assiette statique avant, la valeur de raideur avant correspondant à l'assiette statique avant ASav, par exemple par interpolation linéaire. La raideur de suspension arrière kAR est obtenue en extrayant de la 10 table ou courbe préenregistrée donnant la raideur de la suspension arrière en fonction de l'assiette statique arrière, la valeur de raideur arrière correspondant à l'assiette statique arrière ASar, par exemple par interpolation linéaire. La longueur Ig est calculée par la formule suivante: Ig = (1 RMAvAr).e Le module CGI de la figure 4 met en oeuvre ce calcul de la longueur Ig et fait par exemple partie de l'estimateur 20. La raideur modale en pompage kZ, est calculée comme étant la somme de la raideur de suspension avant kAV et de la raideur de 20 suspension arrière kAR kZ=kAV+kAR La raideur modale en tangage kq, est calculée par la formule ki, = kAV.(Ig)2 + kAR.(e-lg)2 La raideur modale en roulis ke est calculée par la formule ke = Kbadav + Kbadar + v2.(kAV+kAR)/4 Calcul des accélérations modales de caisse A la figure 5, un module MLR calcule, à partir d'une table ou courbe enregistrée donnant une force de flexion en fonction de la longueur du ressort R, la force de flexion absolue FLEXABS correspondant à la valeur réelle LR entrée de cette longueur. Cette courbe enregistrée de flexion tient également compte des butées de suspension, qui sont par exemple en caoutchouc et qui exercent un effort plus important sur la caisse lorsque le ressort appuie sur ces butées en fin de course de l'amortisseur AM. En outre, un module MDEA reçoit en entrée l'assiette statique AS et calcule en fonction de celle-ci l'effort statique correspondant DEAV de flexion sur les roues avant et l'effort statique correspondant DEAR de flexion sur les roues arrière. Un soustracteur SOUS soustrait de la force FLEX-ABS de flexion absolue calculée l'effort statique DEAV ou DEAR, à savoir l'effort DEAV dans le cas de la roue avant A, pour obtenir une force FLB de flexion de ressorts et des butées de suspension, correspondant à la force exercée par le ressort R et les butées d'extrémité sur la caisse 2. Un additionneur AD2 additionne la force FAM d'amortissement, la force FSEC de frottement sec et la FLB de flexion de ressorts et des butées de suspension pour obtenir la force FA selon la formule suivante: FA = FAM + FSEC + FLB. Un module CAL-ACC reçoit en entrée le couple CBAD calculé par le module 11, les forces FA, FB, FC, FD de suspension calculées par l'estimateur 12, ainsi que la masse M de la caisse, le moment d'inertie le en roulis et le moment d'inertie Iq, en tangage, qui sont préenregistrés, pour calculer en fonction de ceux-ci les accélérations modales zG, â. et q, en négligeant l'influence des couples Co et C(p, c'est-à-dire en prenant Co = 0 et = 0, dans un mode de réalisation. Dans le perfectionnement décrit ci-dessous, on tient compte des couples Co et Cq, dans le calcul des accélérations modales. Un module CGI de calcul de grandeur inertielle calcule en fonction de M, le et I9 ainsi que d'une valeur entrée de répartition de masse RMAvAr entre l'avant et l'arrière du véhicule, une masse totale MTOT = MREF du véhicule tenant compte d'un chargement normalisé du véhicule, comme par exemple quatre personnes de 67 kg dans l'habitacle du véhicule et 28 kg de bagages dans le coffre arrière, et la longueur Ig séparant le centre de gravité G de l'axe des roues avant A, B, qui est entrée dans le module CAL- ACC. La valeur de répartition de masse RMAvAr est estimée en permanence à l'aide des valeurs de débattement DEB fournies par les capteurs CAP-DEB de débattement et de la comparaison de chacune de ces valeurs à une moyenne calculée des débattements DEB. Un accéléromètre CAP-ACCT est prévu sur le véhicule pour fournir une accélération transversale ACCT à un estimateur 14 du couple de roulis Co, recevant également en entrée la masse totale MTOT et une valeur RECT de recalage de l'accélération transversale ACCT. L'accéléromètre transversal CAP-ACCT est positionné sur le centre de gravité G et non sur le centre CR de roulis. La valeur RECT de recalage de l'accélération transversale est calculée par le module CAL-ACC de la manière suivante: RECT(n) = ACCT(n) - 6' (n-1).(HCdG - hRoulis) Où d'est l'accélération modale de roulis non filtrée, et Où n désigne la valeur de la variable au cycle actuel et (n-1) désigne la valeur de la variable au cycle précédent. L'estimateur 14 calcule le couple de roulis Co par la formule suivante: co = (ACCT - RECT) . (MTOT) . d(G, CR) Où d(G,CR) = HCdG -hRoulis est la distance entre le centre de 20 gravité G et le centre CR de roulis, et est préenregistrée. Un estimateur 15 du couple C,p de tangage reçoit en entrée la longueur Ig, la masse totale MTOT, une accélération longitudinale ACCL fournie par un accéléromètre longitudinal CAPL disposé dans la caisse du véhicule, une information IF de freinage et une valeur RECL de recalage d'accélération longitudinale calculée par le module CAL-ACC. La valeur RECL de recalage d'accélération longitudinale est calculée par le module CAL-ACC de la manière suivante: RECL(n) = ACCL(n) -(fi (n-1). (HCdG) Où (fi est l'accélération modale de tangage non filtrée. L'estimateur 15 calcule le couple C( de tangage selon la formule suivante: cq, = (ACCL - RECL).(MTOT).hG + CfB hG = HCdG représente la hauteur suivant l'axe Z du centre de gravité G par rapport au centre CT de tangage, et est préenregistrée. La composante c(pB du couple c(i, est la composante du couple de tangage dû à l'effet Brouilhet et est calculée en fonction de l'information IF de freinage. Un module 16 de détermination fournit cette information IF de freinage en fonction d'une valeur de pression PMC de maître cylindre, elle-même fournie par un capteur CAP-P de pression du maître cylindre des freins. Les valeurs calculées des couples Co et C9 sont entrées dans le module CAL-ACC qui calcule à partir de celles-ci et de la valeur des autres entrées et fournit en sortie l'accélération zG modale en pompage, l'accélération modale 6. en roulis, l'accélération modale (fi en tangage, et les valeurs de recalage RECT et RECL. L'accélération modale À' en roulis et l'accélération modale (fi en tangage sont envoyées respectivement à deux convertisseurs Cl et C2 de degrés en radians par seconde pour ensuite être fournis avec zG sur une sortie SACC pour les trois accélérations modales non filtrées et de là à une sortie SACC2 du module 10 vers l'extérieur. En outre, ces trois accélérations modales présentes sur la sortie SACC du module 10 sont chacune envoyées à un filtre 17 éliminant les basses fréquences au-dessous d'une fréquence de coupure basse de par exemple 0.1 Hz, 0.2 Hz ou 0.3 Hz. Le filtre 17 présente par exemple en plus de cette composante passe-haut une composante passe-bas pour former un filtre passe-bande. La fréquence de coupure basse du filtre 17 peut être différente selon l'accélération modale EG, 9' ou cp. Les accélérations modales filtrées en sortie du filtre 17 sont ensuite envoyées à un module 18 intégrateur comportant en sortie un filtre passehaut, fournissant en sortie les vitesses modales de caisse estimées, à savoir la vitesse modale G en pompage de caisse, la vitesse modale 6 en roulis de caisse, et la vitesse modale qi en tangage de caisse, sur une sortie du module 10. Ces vitesses modales de caisse, G en pompage, 0 en roulis, qi en tangage sont des vitesses absolues par rapport à un repère galiléen et sont 5 appelées premières vitesses modales de caisse pour la logique Skyhook de confort. Le calculateur CSS calcule ensuite la grandeur ER de commande de l'actionneur M de l'amortisseur AM de la roue A et des autres roues B, C, D en fonction de ces vitesses modales calculées G, Â et (fi, et fournit les grandeurs ER de commande ainsi calculées aux actionneurs M correspondants sur leur entrée COM de commande. Commande de type Skyhook On décrit ci-dessous le calcul d'un gain modal bmod d'amortissement 15 variable et d'un premier effort Fmod modal de consigne de l'amortisseur pour la commande de l'amortisseur de type confort ou Skyhook . Cette logique de type Skyhook utilise les premières vitesses modales absolues de caisse, G en pompage, Â en roulis, cp en tangage produites par le module 10, symbolisées par le signe général Vmod dans ce qui suit. Niveau de mouvement de la caisse et niveau de tressautement de la caisse Un estimateur 24 est prévu pour calculer, en fonction des débattements DEB des roues, un niveau NMC de mouvement de la caisse et un niveau NTC de tressautement de la caisse. A la figure 9, le niveau NMC de mouvement de la caisse et le niveau NTC de tressautement de la caisse sont obtenus dans l'estimateur 24 par: - calcul de la moyenne DEBAVMOY des débattements des roues avant A, B, - filtrage de la moyenne DEBAVMOY des débattements avant par un 30 filtre passe-bande PB3, pour obtenir une grandeur filtrée DEBAVMOYF, 2890905 22 - prise de la valeur absolue de la grandeur filtrée DEBAVMOYF, dans un module redresseur RED, pour obtenir une grandeur redressée DEBAVMOYF I, - maintien des maxima de la grandeur redressée I DEBAVMOYF 1 5 dans un module MMAX de maintien, fournissant le niveau NMC de mouvement de la caisse. Pour le calcul du niveau NMC de mouvement de la caisse, le filtre passebande PB3 est réglé pour laisser passer les fréquences des mouvements de la caisse, qui sont relativement basses. Le filtre passe-bande PB3 des mouvements de la caisse est par exemple réglé de 0,5 à 2,5 Hz et est proche de la fréquence de résonance de la suspension. Il peut par exemple être choisi entre deux pentes pour obtenir un niveau NMC de mouvement atténué et un niveau NMC de mouvement non atténué. Pour le calcul du niveau NTC de tressautement de la caisse, le filtre passe-bande PB3 est réglé pour laisser passer les fréquences des tressautements de la caisse, qui sont relativement grandes. Le filtre passe-bande PB3 des tressautements de la caisse est par exemple réglé avec une fréquence de coupure basse de 3 Hz et une fréquence de coupure haute de 8 Hz ou plus. Il peut par exemple être choisi entre deux pentes pour obtenir un niveau NTC de tressautement atténué et un niveau NTC de tressautement non atténué. Le module MMAX de maintien peut avoir une pente paramétrable de descente et une temporisation paramétrable de maintien des maxima. La temporisation de maintien des maxima est choisie plus courte pour l'obtention du niveau NTC de tressautement de la caisse que pour l'obtention du niveau NMC de mouvement de la caisse. Efforts modaux de consigne Skyhook et gains modaux Un estimateur 21 est prévu pour calculer les gains modaux bmod 30 d'amortissement variable et les premiers efforts modaux Fmod de consigne d'amortissement selon la formule Fmod = - bmod É Vmod Ilyadonc: - un gain modal de pompage bZ pour calculer le premier effort modal de pompage FZ1 = -bZ. 2G - un gain modal de roulis be pour calculer le premier effort modal de roulis Fol = -be.9 un gain modal de tangage bp pour calculer le premier effort modal de tangage Fp1 = -4q Les gains modaux bZ, b9, bp sont variables en fonction des débattements DEB des roues A, B, C, D et sont calculés par l'estimateur 21 en fonction de grandeurs qui ont été calculées. précédemment en fonction de ces débattements DEB des roues A, B, C, D. Les gains modaux bZ, b9, bp peuvent comprendre un ou plusieurs coefficients multiplicatifs, avec par exemple les coefficients multiplicatifs 15 suivants: - un coefficient multiplicatif de référence bZREF, bOREF, byREF, respectivement en pompage, en roulis et en tangage, - un coefficient multiplicatif d'atténuation bZATT, bOATT, b(pATT en pompage, en roulis et en tangage, - un coefficient multiplicatif de recalage bZREC, boREC, b(pREC, respectivement en pompage, en roulis et en tangage, - un coefficient multiplicatif de type de conduite b7TYF, boTyp, b1TYP, respectivement en pompage, en roulis et en tangage. Dans le mode de réalisation représenté à la figure 6, l'estimateur 21 reçoit en entrée les grandeurs suivantes: - le niveau NMC de mouvement de la caisse fourni par l'estimateur - le niveau NTC de tressautement de la caisse fourni par l'estimateur 3 0 24, - la vitesse VVH du véhicule, -les raideurs modales fournies par l'estimateur 24: la raideur modale en pompage kZ, la raideur modale en tangage kt, et la raideur modale en roulis k9, - les vitesses modales Vmod fournies par le module 10: la vitesse modale G en pompage de caisse, la vitesse modale 6 en roulis de caisse, la vitesse modale cp en tangage de caisse, - les moments modaux d'inertie fournis par l'estimateur 20: le moment d'inertie en roulis le et le moment d'inertie en tangage l, - la masse suspendue MSUS fournie par l'estimateur 20, - une information IS de sportivité, pouvant être à un état booléen 0 d'absence de sportivité, ou à un autre état booléen 1 de sportivité, selon l'actionnement d'un bouton correspondant du tableau de bord du véhicule par son conducteur dans une position respectivement de conduite sportive ou d'absence de conduite sportive. Pour chacun des gains modaux bZ, b9, b(r, , le coefficient multiplicatif de référence bZREF, bOREF, bcpREF, respectivement en pompage, en roulis et en tangage, est obtenu en extrayant d'une table ou courbe de référence préenregistrée donnant le coefficient multiplicatif de référence en fonction de la vitesse du véhicule, la valeur bZREF, bOREF, bcpREF du coefficient multiplicatif de référence, qui correspond à la valeur entrée WH de la vitesse du véhicule, par exemple par interpolation linéaire. Pour chacun des gains modaux bZ, b9, b(, le coefficient multiplicatif d'atténuation bZATT, bOATT, b,Au, respectivement en pompage, en roulis et en tangage, est obtenu - en calculant une résistance RZ, R8, Rp, respectivement en pompage, en roulis et en tangage en fonction du niveau NMC de mouvement de la caisse et du niveau NTC de tressautement de la caisse, par la formule: RZ= NTC-/3Z.NMC 2890905 25 Re = NTC -/3e.NMC Rp= NTC-/39.NMC Où /3z, /3e, sont des paramètres préenregistrés permettant de régler le rapport entre les deux niveaux NMC et NTC, ces paramètres /3z, 5 étant par exemple réglés entre 0,5 et 1; - en extrayant d'une table ou courbe préenregistrée donnant le coefficient multiplicatif d'atténuation beATT, bOATT, b(pATT en fonction de la résistance respectivement en pompage, en roulis et en tangage, la valeur bzATT, bOATT, b(pATT du coefficient multiplicatif d'atténuation, qui correspond à la valeur calculée de la résistance Rzi Re, Rp, respectivement en pompage, en roulis et en tangage, par exemple par interpolation linéaire. Le coefficient multiplicatif d'atténuation beATT, bOATT, b(pATT en pompage, en roulis et en tangage est par exemple donné par la formule: beATT = 1 /(1 +az. Rz) bOATT = 1 /(1 +ae. Re) b(Au = 1 /(1 +a(p. R(p) où az, ae, a(p sont des paramètres préenregistrés. La valeur obtenue de bzpn, bOATT, b(pATT n'est par exemple retenue que si la résistance Rz, Re, Rp associée est supérieure à un seuil prescrit. Si la résistance Rzi Re, R(p associée est inférieure ou égale à ce seuil prescrit, on prend 1 comme coefficient multiplicatif d'atténuation beATT, bOATT, b(pATT. Pour chacun des gains modaux bz, be, b(p, le coefficient multiplicatif de recalage bZREC, bOREC, bcpREC, respectivement en pompage, en roulis et en 25 tangage, est obtenu par la formule bzREC - beREC - kZ É MSUS kZREF MREF ke ÉIe 4REF. IBREF bcpREC = 1 k(pnEF ' IcpnEr où kZREF est une raideur de référence en pompage, constante, keREF est une raideur de référence en roulis, constante, kWnEF est une raideur de référence en tangage, constante, IbREF est un moment d'inertie de référence en roulis, constant, IWREF est un moment d'inertie de référence en tangage, constant, kzREF, kOREF, kWREF, MREF, I'REF, IWREF sont des paramètres préenregistrés, correspondant à un chargement normalisé du véhicule, comme par exemple quatre personnes de 67 kg dans l'habitacle du véhicule et 28 kg 10 de bagages dans le coffre arrière. Pour chacun des gains modaux bZ, b9, bW, le coefficient multiplicatif de type de conduite bzTYP, boTyp, bcpTYP, respectivement en pompage, en roulis et en tangage, est égal à un gain préenregistré de sportivité GSz, GS0, GST, si l'information IS de sportivité est à l'état booléen 1 de sportivité et est égal à 1, si l'information IS de sportivité est à l'état booléen 0 d'absence de sportivité. Les gains modaux bZ, b9, bW sont calculés en fonction des coefficients 20 multiplicatifs selon les formules: bz = bzREF É bzATT É bzREC É bzTYP bo bOREF É bOATT É boREC É b9TYP bW = b(pREF É bq,ATT É bcpREC É bcpTYP Le premier effort modal de pompage Fz1, le premier effort modal de roulis F91, le premier effort modal de tangage FW1 sont calculés et sont également appelés efforts modaux de confort ou skyhook . Le premier effort modal de pompage Fz1, le premier effort modal de roulis F91, le kW IW premier effort modal de tangage Fq,l sont fournis sur une sortie de l'estimateur 21. Logique Roadhook Ci-dessous est décrit la logique de type Roadhook, c'està-dire qui suit le profil de la route, cette logique étant également appelée logique de tenue de caisse ou logique de comportement. Le principe de cette logique de tenue de caisse est de faire tendre vers zéro ou de minimiser l'une ou plusieurs des accélérations modales de la caisse par rapport au plan des roues: accélération modale en pompage, accélération modale en roulis, accélération modale en tangage. A la figure 10, le dispositif comporte un estimateur 31 des vitesses modales Vmod2 de la caisse par rapport au plan moyen des roues en fonction des débattements DEB mesurés des roues A, B, C, D. Ces vitesses modales Vmod2 par rapport au plan moyen des roues sont appelées vitesses relatives et comprennent la vitesse modale relative zG2 de caisse en pompage, la vitesse modale relative 02 de caisse en tangage et la vitesse modale relative 62 de caisse en roulis. Cet estimateur 31 des vitesses modales relatives Vmod2 reçoit en entrée: les débattements DEB mesurés sur les roues A, B, C, D, - la voie v, - au moins deux des paramètres suivants: la valeur de répartition de 25 masse RMAvAr entre l'avant et l'arrière, la longueur Ig séparant le centre de gravité G de l'axe des roues avant A, B, et l'empattement e. Les débattements DEB sont d'abord filtrés dans un filtre passe-bas, par exemple d'ordre 2 et de type Butterworth, pour n'obtenir que les débattements basses fréquences et éliminer dans une large mesure les tressautements hautes fréquences. Puis, un dérivateur dérive les débattements DEB ainsi filtrés pour obtenir les vitesses de débattement Roadhook des roues A, B, C, D. Les vitesses modales relatives Vmod2 sont ensuite calculées selon les formules suivantes: - vitesse modale relative de caisse en pompage par rapport au plan moyen des roues: (e -1g) dA + dB Ig de + dD Zcz _ + e 2 e 2 - vitesse modale relative 'de caisse en tangage par rapport au plan moyen des roues: dA+dB dc dD (P2 2e - vitesse modale relative de caisse en roulis par rapport au plan moyen des roues: 92dA dB dc+dD 2v avec dA = vitesse de débattement VDEB de la roue A avant gauche, dB = vitesse de débattement VDEB de la roue B avant droite, = vitesse de débattement VDEB de la roue C arrière droite, dD = vitesse de débattement VDEB de la roue D arrière gauche. Jerk transversal anticipé Un estimateur 32 est prévu pour calculer un jerk transversal anticipé f (dérivée troisième de la coordonnée Y par rapport au temps) à partir de la vitesse WH mesurée du véhicule et de la vitesse de rotation b du volant de direction du véhicule, où 6 est l'angle de rotation mesuré de ce volant, 25 mesuré par tout capteur ou moyen approprié. Cet estimateur 32 reçoit en entrée: - la masse suspendue MSUS, - la valeur de répartition de masse RMAvAr entre l'avant et l'arrière, - la vitesse WH du véhicule, - la vitesse de rotation S du volant de direction. Le jerk transversal anticipé @ est estimé selon la formule: D.S.VVH2 Y _ e(1+K.VVH2) où D est la démultiplication du volant de direction et K est une constante de gain de surbraquage, calculée en fonction de la valeur de répartition de masse RMAvAr entre l'avant et l'arrière et de la masse suspendue MSUS. Le gain K de surbraquage est une grandeur du véhicule, déterminée par des mesures sur le véhicule. Couple moteur anticipé aux roues Un estimateur 40 est prévu pour calculer ce couple moteur anticipé aux roues, désigné par CR. Pour ce faire, on estime le numéro i du rapport embrayé REMBR(i) de la boîte de vitesse du véhicule, qui va par exemple de 1 à 5. On calcule le vitesse WH1 qu'aurait le véhicule à une vitesse de rotation prescrite WMOT1 du moteur, qui ne dépend en position embrayée que du rapport REMBR engagé, selon la formule VVH1 = VVH. wMOT1 / wMOT où wMOT est la vitesse de rotation du moteur à la vitesse WH du véhicule. Par exemple, oMOT1 = 1000 tours / minute. On calcule pour chaque rapport i d'embrayage, les paramètres P; = 0,5. (VVH1(i) + VVHI (i+1)). En comparant WH1 à P; et en retenant la valeur P; la plus proche de WH1, on déduit le rapport i d'embrayage. Le couple moteur anticipé CR aux roues est alors: CR = CM. REMBR(i), avec REMBR(I) = WMOT/wROUE Où REMBR(i) est le rapport d'embrayage ayant le numéro i, CM est le couple moteur déterminé par tout moyen approprié, par exemple par un calculateur du contrôle moteur. (ROUE est la vitesse de rotation des roues. Jerk longitudinal anticipé Un estimateur 33 est prévu pour calculer un jerk longitudinal anticipé X (dérivée troisième de la coordonnée X par rapport au temps) à partir de la dérivée du couple moteur anticipé et de la dérivée PMC de la pression Pm du maître cylindre. Cet estimateur 33 reçoit en entrée: la masse suspendue MSUS, la pression PMC du maître cylindre, le couple moteur anticipé aux roues CR. Le calcul est effectué de la manière suivante. On extrait d'abord d'une table ou courbe préenregistrée donnant un effort de freinage du maître cylindre en fonction de la pression du maître cylindre, la valeur EFR de cet effort de freinage correspondant à la pression PMc du maître cylindre, par exemple par interpolation linéaire. On applique ensuite un filtre passe-bas, par exemple d'ordre 1 de type Butterworth, à cet effort EFR de freinage, et on dérive dans un dérivateur l'effort EFR de freinage ainsi filtré pour obtenir la dérivée EFRF de l'effort EFR filtré. On calcule un effort anticipé du moteur aux roues EMR, égal au couple moteur anticipé aux roues CR, divisé par un rayon moyen Rmoy des roues, prédéterminé et préenregistré. On applique ensuite un filtre passe-bas, par exemple d'ordre 1 de type Butterworth, à cet effort EMR anticipé du moteur aux roues, et on dérive dans un dérivateur l'effort EMR anticipé du moteur ainsi filtré pour obtenir la dérivée ÉMU. de l'effort EMR filtré. Le jerk longitudinal anticipé est alors égal à la somme des dérivées EFRF, , divisée par la masse totale MTOT: 2890905 31 X _ EFRF- + EMMF MTOT Dans cette formule, la masse totale MTOT inclut la masse suspendue MSUS, peut inclure la masse des roues et peut être bornée entre deux seuils. Ces jerks et sont estimés et ne proviennent pas d'une dérivation d'accéléromètres, trop bruités et trop en retard. Termes d'efforts modaux anticipatifs Un module 34 est prévu pour calculer des termes d'efforts modaux 10 anticipatifs, à savoir: - un couple modal anticipatif en tangage, désigné par cq,2ant, - un couple modal anticipatif en roulis, désigné par celant. Il n'est pas calculé d'effort modal anticipatif en pompage, étant donné que seul un effort modal correctif Roadhook agit en pompage, ainsi 15 que cela sera décrit ci-dessous. Dans le mode de réalisation représenté à la figure 11, l'estimateur 34 reçoit en entrée les grandeurs suivantes: - le jerk transversal anticipé f fourni par l'estimateur 32, - le jerk longitudinal anticipé fourni par l'estimateur 33, - la vitesse VVH du véhicule, - les raideurs modales fournies par l'estimateur 24: la raideur modale en pompage kZ, la raideur modale en tangage kI, et la raideur modale en roulis ke, - les vitesses modales relatives Vmod2 par rapport au plan moyen des 25 roues, fournies par le module 31: la vitesse modale relative G2 en pompage de caisse, la vitesse modale relative 02 en roulis de caisse, la vitesse modale relative q52 en tangage de caisse, - les moments modaux d'inertie fournis par l'estimateur 20: le moment d'inertie en roulis le et le moment d'inertie en tangage Iq), 2890905 32 - la masse suspendue MSUS fournie par l'estimateur 20, - l'information IS de sportivité. Ainsi que cela est représenté à la figure 11, chacun de ces termes d'efforts modaux anticipatifs en tangage Ccp2ant et en roulis ce2ant est calculé par traitement respectivement du jerk longitudinal anticipé X et du jerk transversal anticipé Y, pour obtenir un jerk longitudinal anticipé traité XT et un jerk transversal anticipé @T, puis multiplication respectivement par un gain en sollicitation longitudinale ' Gsx pour obtenir le couple modal anticipatif en tangage Ccp2ant et par un gain en sollicitation longitudinale GsY pour obtenir le couple modal anticipatif en roulis ce2ant selon les formules: C(p2ant = GSX É XT ce2ant = GSY É YT Le gain en sollicitation longitudinale Gsx et le gain en sollicitation transversale GsY sont des paramètres prédéterminés de mise au point, déterminés par des essais sur véhicule pour obtenir les bonnes prestations de tenue de caisse sous sollicitation du conducteur. Cette mise en forme est décrite ci-dessous pour le calcul du couple modal anticipatif en tangage, désigné par Ccp2ant à partir du jerk longitudinal anticipé X - passage du jerk longitudinal anticipé X dans un filtre 341 d'annulation des faibles amplitudes, possédant un seuil haut positif SHJL d'activation de jerk longitudinal et un seuil bas négatif SBJL d'activation de jerk longitudinal, pour remplacer par des valeurs nulles les valeurs du jerk longitudinal anticipé X situées entre le seuil haut d'activation de jerk longitudinal SHJL et le seuil bas d'activation de jerk longitudinal SBJL au cours du temps; - passage du jerk longitudinal anticipé X filtré, issu du filtre 341, dans un module 342 de maintien des maxima, pouvant avoir une temporisation paramétrable de maintien des maxima, pour obtenir un jerk filtré et maintenu à ses maxima, désigné parXf,,,aX - passage du jerk Xf,,,ax filtré et maintenu à ses maxima, issu du module 342, dans un module 343 limiteur de pente, limitant en valeur absolue la pente de la descente du jerk Xf.filtré et maintenu à ses maxima, pour obtenir le jerk longitudinal anticipé traité XT, puis multiplication respectivement par le gain en sollicitation longitudinale Gsx pour obtenir le couple modal anticipatif en tangage c92antÉ La temporisation doit être suffisamment longue pour que le terme correctif Roadhook (voir supra) ait le temps de devenir significatif sur action simple (mise en virage simple, freinage ou accélération) et être suffisamment courte pour ne pas perturber le fonctionnement Roadhook et pour ne pas demander d'amortissement inutile. Le passage du jerk transversal anticipé f dans le filtre 341 d'annulation, ayant son seuil haut positif SHJT d'activation de jerk transversal et son seuil bas négatif SBJT d'activation de jerk transversal, puis dans le module 342 de maintien des maxima produit un jerk filtré et maintenu à ses maxima, désigné par @f x, qui est envoyé au module 343 limiteur de pente ayant le gain en sollicitation transversale Gsy pour fournir en sortie le couple modal anticipatif en roulis celant. Les seuils hauts SHJT et SHJL peuvent être égaux et opposés aux seuils bas égaux SBJT et SBJL. Ces seuils sont paramétrables et sont un compromis entre la limitation des actions intempestives et le non traitement des petites sollicitations. De préférence, chacun des seuils SHJT, SHJL, SBJT et SBJL est compris entre 1 et 10 m.s-3. Le fait d'utiliser des termes anticipatifs permet de gagner en temps de réponse pour placer les actionneurs dans le bon état avant que la caisse ait le temps de prendre de la vitesse. Il en résulte une amélioration notable de la tenue de caisse. Termes d'effort modaux correctifs Le module 34 calcule également, en fonction de la vitesse modale relative Vmod2 = ZG2, q 2, 02 par rapport au plan moyen des roues, au moins un deuxième terme d'effort modal correctif F2COR selon la formule générale F2coR = - bmod2 É Vmod2 à savoir: - une deuxième force modale corrective en pompage, désignée par Fz2cor, - un deuxième couple modal correctif en tangage, désigné par C(p2c0r, 10 - un deuxième couple modal correctif en roulis, désigné par Ce2cor selon les formules: Fz2cor = -bZ2 zc2 Ccp2cor = -bcp2 É PZ Ce2cor = -be2 É e2 où bmod2 est un deuxième gain modal correctif, bz2 est un deuxième gain modal correctif de pompage pour calculer la deuxième force modale corrective en pompage Fz2cor, bu est un deuxième gain modal correctif de roulis pour calculer le deuxième couple modal correctif en roulis Ce2cor, 42 est un deuxième gain modal correctif de tangage pour calculer le deuxième couple modal correctif en tangage c1,2corÉ Les deuxièmes gains modaux correctifs bZ2, be2, 42 peuvent comprendre un ou plusieurs coefficients multiplicatifs, avec par exemple les coefficients multiplicatifs suivants: - un deuxième coefficient multiplicatif de référence bzREF2, beREF2, bcpREF2, respectivement en pompage, en roulis et en tangage, - un deuxième coefficient multiplicatif de recalage bzREC2, bOREC2, bTREC2, respectivement en pompage, en roulis et en tangage, - un deuxième coefficient multiplicatif de type de conduite bZTYp2, 30 boTYp2, bcpTYp2, respectivement en pompage, en roulis et en tangage. Pour chacun des deuxièmes gains modaux b,2, be2, bgp2, le deuxième coefficient multiplicatif de référence bzREF2, beREF2, b9REF2, respectivement en pompage, en roulis et en tangage, est obtenu en extrayant d'une deuxième table ou courbe de référence, préenregistrée pour la logique Roadhook, donnant le deuxième coefficient multiplicatif de référence en fonction de la vitesse du véhicule, la valeur bzREF2, bOREF2, bcpREF2 du deuxième coefficient multiplicatif de référence, qui correspond à la valeur entrée VVH de la vitesse du véhicule, par exemple par interpolation linéaire. Pour chacun des deuxièmes gains modaux bz2, be2, 42, le deuxième coefficient multiplicatif de recalage bzREC2, bOREC2, byREC2 est par exemple égal au premier coefficient multiplicatif de recalage bzREC, beREC, bÇREC respectivement en pompage, en roulis et en tangage, décrit ci-dessus: bzREC2 = bzREC, bOREC2 = boREC, b9REC2 = b(pRECÉ Pour chacun des deuxièmes gains modaux bz, be, b1,, le deuxième coefficient multiplicatif de type de conduite bzTYP2, beTYP2, bTTYP2, respectivement en pompage, en roulis et en tangage, est par exemple égal au premier coefficient multiplicatif de type de conduite bzTYP, b0TYP, bgTYP, décrit ci-dessus: bzTYP2 = bzTYP, boTYP2 = boTYP, b(pTYP2 = b(pTYPÉ Les deuxièmes gains modaux correctifs bz2, be2, b1,2 sont calculés en fonction des deuxièmes coefficients multiplicatifs selon les formules: bz2 = bzREF2 É bzREC2 É bzTYP2 be2 = b9REF2 É bOREC2 É boTYP2 bcp2 = bcpREF2 bcpREC2 É bcpTYP2 Efforts modaux Roadhook L'estimateur 34 fait ensuite - la synthèse entre le couple modal anticipatif en tangage cp2ant et le deuxième couple modal correctif en tangage C 2cor, pour obtenir en sortie le deuxième couple ou effort modal c(p2 de tangage, - la synthèse entre le couple modal anticipatif en roulis celant et le 5 deuxième couple modal correctif en roulis ce2cor, pour obtenir en sortie le deuxième couple ou effort modal cet de roulis. La deuxième force modale corrective en pompage, désignée par F22cor est prise en sortie comme deuxième force ou effort modal Fz2 de pompage F22 = Fz2cor. Ces deuxièmes efforts c(2, c02 et F22 sont appelés efforts modaux de comportement ou de tenue de route ou Roadhook . La synthèse est effectuée par choix du terme anticipatif ou du terme correctif en fonction de leurs valeurs, de la manière décrite dans le tableau ci-dessous. Terme Anticipatif -- Petit Grand Terme correctif. Petit Cas 1: Terme correctif Cas 3: Terme anticipatif Grand Cas 2: Terme correctif Cas 4: - Maximum des 2 si de même signe - Terme correctif si signes opposés Pour l'obtention du deuxième effort modal cq)2 de tangage, celui-ci est égal - au deuxième couple modal correctif en tangage ccp2c0r, lorsque la valeur absolue du couple anticipatif en tangage ccp2ant est inférieure ou égale à une première valeur prescrite de tangage V1cp, (cas 1 et 2 dans le tableau, correspondant au terme anticipatif petit), - au couple modal anticipatif en tangage c92ant, lorsque la valeur absolue du couple anticipatif en tangage C(p2ant est supérieure à la première valeur prescrite de tangage V1(p, et lorsque la valeur absolue du couple modal correctif en tangage ccp2cor est inférieure ou égale à une deuxième valeur prescrite de tangage V2(p (cas 3 dans le tableau, correspondant au terme correctif petit et au terme anticipatif grand). Si la valeur absolue du couple anticipatif en tangage Ccp2ant est supérieure à la première valeur prescrite de tangage V1 (p et si la valeur absolue du couple modal correctif en tangage Ccp2cor est supérieure à la deuxième valeur prescrite de tangage V2(p (cas 4 dans le tableau, correspondant au terme correctif grand et au terme anticipatif grand), -si le couple modal correctif en tangage Ccp2cor et le couple anticipatif en tangage Ccp2ant ont le même signe, le deuxième effort modal c92 de tangage est égal à max(I c(p2cor I, I 42ant I).sgn(cç2ant), Où sgn désigne la fonction signe et max la fonction maximum, et - si le couple modal correctif en tangage Ccp2cor et le couple anticipatif en tangage C(p2ant n'ont pas le même signe, le deuxième effort modal ce de tangage est égal au couple modal correctif en tangage Ccp2cor. L'obtention du deuxième effort modal c02 de roulis est analogue à ce qui précède, à partir de ce2cor et co2ant au lieu de ccp2cor et Ccp2ant, avec une première valeur prescrite de roulis V10 au lieu de V1cp, et une deuxième valeur prescrite de roulis V20 au lieu de V2(p. Synthèse entre Skyhook et Roadhook Le premier effort modal de pompage Fz1, le premier effort modal de roulis F91 et le premier effort modal de tangage Fo fournis par l'estimateur 21 (efforts modaux de confort par la logique Skyhook, désignés d'une manière générale par premiers efforts modaux F1 de consigne), ainsi que le deuxième effort modal Fz2 de pompage, le deuxième effort modal CO2 de roulis et le deuxième effort modal c92 de tangage fournis par l'estimateur 34 (efforts modaux de comportements par la logique Roadhook, désignés d'une manière générale par deuxièmes efforts modaux F2 de consigne), sont envoyés à un estimateur 22 d'une force de consigne pour chaque amortisseur, soit pour les roues A, B, C, D, les forces de consigne FA1, FB1, FC1, FD1. Pour chaque mode, l'estimateur 22 pondère le premier effort FI de confort et le deuxième effort F2 de comportement pour calculer l'effort modal F de consigne. L'estimateur 22 calcule: - un effort modal F = Fz de consigne en pompage, en fonction du premier effort de pompage Fz1 de confort, du deuxième effort de pompage Fz2 de comportement et d'un coefficient a de pondération, selon la formule: Fz = a. Fz2 + (1-a) . Fz1 - un effort modal F = F,p de consigne en tangage, en fonction du premier effort modal de tangage Fo de confort, du deuxième effort modal c(p2 de tangage de comportement et du coefficient a de pondération, selon la formule: F = a. cq,2 + (1-a) . Fo - un effort modal F = Fe de consigne en roulis, en fonction du premier effort modal de roulis F91 de confort, du deuxième effort modal c92 de roulis de comportement et du coefficient a de pondération, selon la formule: FB=a. cA2+(1-a). F91 La détermination de ce coefficient a de pondération en fonction de sollicitations détectées est décrite ci-dessous. Le coefficient de pondération est normalement à 0 pour faire suivre aux efforts modaux de consigne les premiers efforts de confort Fz1, Fol et F1 de la logique Skyhook. Accélération longitudinale corrigée L'accélération longitudinale corrigée COR est calculée par un estimateur 25 à partir de l'accélération longitudinale mesurée ACCL, fournie par l'accéléromètre longitudinal CAPL. L'estimateur 25 reçoit en entrée: - la vitesse WH mesurée du véhicule, la masse suspendue MSUS, fournie par l'estimateur 20, -l'accélération longitudinale mesurée ACCL, - la pression PMC du maître cylindre des freins, fournie par le capteur CAP-P, - le couple moteur anticipé aux roues CR, fourni par l'estimateur 40. Le calcul est effectué de la manière suivante. On extrait d'abord de la table ou courbe préenregistrée donnant l'effort de freinage du maître cylindre en fonction de la pression du maître cylindre, la valeur EFR de cet effort de freinage correspondant à la pression PMC du maître cylindre, par exemple par interpolation linéaire. On calcule l'effort anticipé du moteur aux roues EMR, égal au couple moteur anticipé aux roues CR, divisé par un rayon moyen Rmoy des roues, prédéterminé et préenregistré. On calcule un effort de traînée longitudinale ETR en fonction de la 20 vitesse du véhicule WH selon la formule: ETR = COEF. (VVH)2 + DEC Où COEF est un coefficient prédéterminé et préenregistré et DEC est un décalage prédéterminé et préenregistré. L'effort longitudinal total ELT est alors égal à la somme de l'effort 25 EFR de freinage, de l'effort EMR anticipé du moteur aux roues et de l'effort ETR de traînée longitudinale: ELT = EFR + EMR + ETR On calcule la masse totale MTOT, qui inclut la masse suspendue MSUS, peut inclure la masse des roues et peut être bornée entre deux seuils. L'accélération longitudinale anticipée XANT est calculée en divisant l'effort longitudinal total ELT par la masse totale MTOT: XANT=ELT/MTOT L'accélération longitudinale anticipée X. est éventuellement ensuite bornée entre deux seuils. L'accélération longitudinale corrigée XcoR est ensuite calculée en calculant une évolution EVAL de l'accélération longitudinale, égale à l'accélération longitudinale anticipée XANT à laquelle est soustraite l'accélération longitudinale mesurée ACCL: EVAL = X ANT - ACCL - en appliquant à cette évolution EVAL de l'accélération longitudinale un filtre passe-haut PH, par exemple de type Butterworth d'ordre 1, pour obtenir l'évolution longitudinale EVAL filtrée, égale à PH(XANT- ACCL), en additionnant l'évolution longitudinale EVAL filtrée à l'accélération 15 longitudinale mesurée ACCL, pour obtenir l'accélération longitudinale corrigée XcoR XcoR = ACCL + PH(XANT - ACCL) La fréquence de coupure du filtre passe-haut PH permet de régler la vitesse de recalage de l'estimation sur la mesure. Accélération transversale corrigée L'accélération transversale corrigée YcoR est calculée par un estimateur 26 à partir de l'accélération transversale mesurée ACCT, fournie par l'accéléromètre transversal CAPACCT. L'estimateur 26 reçoit en entrée: la masse suspendue MSUS, la valeur de répartition de masse RMAvAr entre l'avant et l'arrière, la vitesse WH du véhicule, - l'angle de rotation 8 du volant de direction, -l'accélération transversale mesurée ACCT. L'accélération transversale anticipée TANT est estimée selon la formule: D.6.VVHZ Y = e(1+K.VVH2) où D est la démultiplication du volant de direction et K est la constante de gain de surbraquage, calculée en fonction de la valeur de répartition de masse RMAvAr entre l'avant et l'arrière et de la masse suspendue MSUS. La constante K de gain de surbraquage est une grandeur du véhicule, déterminée par des mesures sur le véhicule. L'accélération longitudinale anticipée @ANT est éventuellement ensuite bornée entre deux seuils. L'accélération longitudinale corrigée @coR est ensuite calculée en calculant une évolution EVAT de l'accélération transversale, égale 15 à l'accélération transversale anticipée TANT à laquelle est soustraite l'accélération transversale mesurée ACCT: EVAT = @ANT - ACCT - en appliquant à cette évolution EVAT de l'accélération transversale un filtre passe-haut PH2, par exemple de type Butterworth d'ordre 1, pour 20 obtenir l'évolution transversale EVAT filtrée, égale à PH( @ANT - ACCT), - en additionnant l'évolution transversale EVAT filtrée à l'accélération transversale mesurée ACCT, pour obtenir l'accélération transversale corrigée GR YcoR = ACCT + PH2(@ANT - ACCT) La fréquence de coupure du filtre passe-haut PH2 permet de régler la vitesse de recalage de l'estimation sur la mesure. Détection des sollicitations et coefficient de pondération des efforts Skyhook et des efforts Roadhook 2890905 42 A la figure 12, un estimateur 23 calcule le coefficient a de pondération des premiers efforts de confort et des deuxièmes efforts de comportement. L'estimateur 23 reçoit en entrée: - le jerk longitudinal anticipé , fourni par l'estimateur 33, - le jerk transversal anticipé , fourni par l'estimateur 32, - l'accélération longitudinale corrigée XcoR, fournie par l'estimateur - l'accélération transversale corrigée G R, fournie par l'estimateur 26, - l'information IS de sportivité. Par défaut, les premiers efforts de confort FZ1, Fei et Fo de la logique Skyhook sont sélectionnés pour les efforts modaux de consigne, c'est-àdire que le coefficient a de pondération est à O. Les sollicitations sont détectées sur les valeurs prises par ces entrées. Dès qu'une sollicitation est détectée, le coefficient a de pondération passe en tout comportement ou Roadhook, c'est-à-dire à 1, pour sélectionner les deuxièmes efforts FZ2, cet et c92 de comportement comme efforts modaux de consigne. Si l'on détecte une stabilisation au sein d'une sollicitation, typiquement une grande courbe d'autoroute comme à la figure 14, il y a la possibilité de refaire passer le coefficient a de pondération progressivement à 0 en logique Skyhook pour favoriser le confort. Si, au sein de cette stabilisation, on détecte une variation des grandeurs accélérométriques, la répartition repasse immédiatement en tout 'comportement', c'est-à-dire à 1. On crée un signal booléen 'sollicitation conducteur latérale' (SSOLT) et un signal booléen 'sollicitation conducteur longitudinale' (SSOLL) sur franchissement de seuils paramétrables d'accélération corrigée ou de jerk anticipé. Le coefficient de pondération passe à 1 et on réinitialise la temporisation lorsque les évènements suivants sont détectés: - front montant de sollicitation conducteur longitudinale, - front montant de sollicitation conducteur latérale, - dépassement de seuil sur jerk longitudinal en sollicitation conducteur longitudinale, - dépassement de seuil sur variation d'accélération longitudinale en 5 sollicitation conducteur longitudinale, - dépassement de seuil sur jerk latéral en sollicitation conducteur latérale, - dépassement de seuil sur variation d'accélération latérale en sollicitation conducteur latérale. L'estimateur 23 détermine en fonction de l'information IS de sportivité une modulation de seuil MODL longitudinale et une modulation de seuil MODT transversale. Si l'information IS de sportivité est égale à 1, la modulation de seuil MODL longitudinale est égale à une valeur prescrite longitudinale plus petite que 1 et la modulation de seuil MODT transversale est égale à une valeur prescrite transversale plus petite que 1. Si l'information IS de sportivité est égale à 0, la modulation de seuil MODL longitudinale est égale à 1 et la modulation de seuil MODT transversale est égale à 1. On détermine ensuite des signaux de détection de sollicitations: un signal logique longitudinal SSOLL de sollicitation, un deuxième signal logique longitudinal SL2, un troisième signal logique longitudinal SL3, un signal logique transversal SSOLT de sollicitation, un quatrième signal logique transversal ST4 et un cinquième signal logique transversal ST5, de la manière suivante: - si I XcoR I > THAL1. MODL ou X I > THJL1. MODL alors SSOLL = 1, - sinon SSOLL = 0. - si SSOLL = 1 et I X I > THJL2 alors SL2 = 1, - sinon SL2 = 0. - l'accélération longitudinale yL est initialisée à O. - si I XcoR - YL I > THAL2. I yL I alors É YL = XcoR est enregistrée pour le prochain calcul de SL3, É si SSOLL = 1 alors SL3 = 1 et sinon SL3 = 0, 10 - si I XcoR - YL I THAT1. MODT ou I I > THJT1. MODT alors SSOLT = 1, - sinon SSOLT = O. -si SSOLT = 1 et I Y I > THJT2 alors ST4 = 1, - sinon ST4 = O. - l'accélération transversae yT est initialisée à O. - si I YcoR - YT I > THAT2. I yT I alors 25 É yT = GR est enregistrée pour le prochain calcul de ST5, É si SSOLT= 1 alors ST5 = 1 et sinon ST5 = 0, - si 1 GR - YT I 2890905 45 THAL1 est un premier seuil d'accélération longitudinale, THAL2 est un deuxième seuil de variation d'accélération longitudinale, THJL1 et THJL2 sont des premier et deuxième seuils de jerk longitudinal, THAT1 est un premier seuil d'accélération transversale, THAT2 est un deuxième seuil de variation d'accélération transversale, THJT1 et THJT2 sont des premier et deuxième seuils de jerk transversal, ces seuils étant préenregistrés. Les états 1 des signaux de détection correspondent à des états de présence d'une sollicitation et les états 0 correspondent à des états 10 d'absence de sollicitations. Un signal logique SSOL de sollicitation du conducteur est déterminé comme étant égal à 1 si le premier signal logique longitudinal SSOLL de sollicitation est à 1 et/ou si le signal logique transversal SSOLT de sollicitation est à 1 (Opérateur logique OU non exclusif). Un premier signal logique SL1 est rendu égal au signal logique SSOL de sollicitation du conducteur. On détermine, en fonction de l'information IS de sportivité, un temps TMOD de modulation entre les premiers efforts Skyhook et les deuxièmes 20 efforts Roadhook: - si IS = 1 alors le temps TMOD de modulation est égal à TMOD = TPER. MODSPORT, - sinon TMOD = TPER, où TPER est un temps de régime permanent prédéterminé et préenregistré, qui représente le temps de passage de la logique Roadhook vers la logique Skyhook en régime établi, et MODSPORT est un facteur multiplicatif de temps de modulation en cas de choix de conduite sportive, qui est plus grand que 1 et qui est prédéterminé et préenregistré. A la figure 13, représentant des chronogrammes en fonction du 30 temps t, un coefficient aINTER de pondération intermédiaire est ensuite calculé de la manière suivante: - initialisation à 0, (stade S10) - mise à 1 du coefficient aINTER de pondération intermédiaire sur chaque front montant détecté de l'un, plusieurs ou tous les premier, deuxième, troisième, quatrième, cinquième signaux logiques SL1 = SSOL, SL2, SL3, ST4, ST5 (stade S11), - maintien à 1 du coefficient aINTER de pondération intermédiaire pendant un temps mort TMORT prédéterminé et préenregistré après chacun de ces fronts montants détectés, (stade S12), - décroissance du coefficient aINTER de pondération intermédiaire, par exemple linéaire, pendant le temps TMOD de modulation après ce temps 10 mort TMORT (stade S13), jusqu'à 0, - si un nouveau front montant est détecté, le coefficient (INTER de pondération intermédiaire est remis à 1 suivant le stade S11 et le processus S11, S12, S13 décrit ci-dessus est recommencé. On calcule un signal logique limité SSOLLIMIT de sollicitation du conducteur, en filtrant le signal logique SSOL de sollicitation du conducteur dans un limiteur de pente négative, pour qu'il passe de 1 à 0 au minimum en le temps TMOD de modulation. Le coefficient a de pondération est égal au coefficient (INTER de pondération intermédiaire, multiplié par le signal logique limité SSOLLIMIT de 20 sollicitation du conducteur: a = aINTER É SSOLLIMIT La figure 14 montre les chronogrammes de l'angle 6 du volant au cours d'une mise en virage simple, qui provoque le passage du coefficient a de pondération à 1 (Roadhook) au début et à la fin du virage, alors que le coefficient a de pondération est à 0 (Skyhook) avant le virage, après le virage et au milieu du virage. Forces de consigne aux roues On extrait de la table ou courbe préenregistrée donnant le coefficient de répartition d'effort sur l'avant en fonction de la valeur de répartition de masse entre l'avant et l'arrière, la valeur du coefficient CAV de répartition d'effort sur l'avant correspondant à la valeur de répartition de masse RMAvAr entre l'avant et l'arrière, par exemple par interpolation linéaire. Ce coefficient CAV de répartition d'effort sur l'avant est supérieur ou égal à 0 et inférieur ou égal à 1. Un ratio anti-devers RAD, supérieur ou égal à 0 et inférieur ou égal à 1, est calculé en fonction de la vitesse VVH du véhicule. Par exemple, on extrait de la table ou courbe préenregistrée donnant le ratio anti-devers en fonction de la vitesse du véhicule, la valeur du ratio anti-devers RAD correspondant à la vitesse VVH du véhicule, par exemple par interpolation linéaire. L'estimateur 22 calcule les forces de consigne des amortisseurs AM aux roues A, B, C, D en fonction des efforts modaux FZ, Fe et Fq, de consigne, selon les formules suivantes: - la force FA1 de consigne pour la roue avant gauche A: 15 FA1= F.CAV F F.RAD 2 2.e v - la force FB1 de consigne pour la roue avant droite B: = F.CAV F + FB.RAD FB1 2 2.e v la force FC1 de consigne pour la roue arrière droite C: FCl F.(1 CAV)+F +F.(1 RAD) 2 2.e v - la force FD1 de consigne pour la roue arrière gauche D: FDl=F. (1 CAV)+F Fo.(l RAD) 2 2.e v L'estimateur détermine ensuite à partir des forces FA1, FB1, FC1, FD1 de consigne des amortisseurs AM aux roues A, B, C, D, et de la vitesse VDEB de débattement valable respectivement pour ces roues A, B, C, D, la loi ERc =ERcA, ERcB, ERcc, ERcp d'amortissement de consigne que doit prendre l'amortisseur AM de la roue A, B, C, D, par exemple en positionnant le point (VDEB(A) ; FA1) sur le graphique de la figure 15 et en recherchant la loi ER d'amortissement la plus proche. Etats minimaux Un estimateur 27 calcule des états minimaux d'amortissement. Cette fonction permet d'interdire à la suspension d'occuper des états d'amortissement trop souples, en imposant des états ERM minimaux, c'est-à-dire des lois ERM minimales d'amortissement, et ce, en fonction de 4 flux différents en entrée: - la vitesse du véhicule pour l'obtention du premier état minimal ERM1: ce critère est utilisé pour les situations de vie du véhicule à l'arrêt, ou à vitesse très faible (pour descentes de trottoirs par exemple), ou à vitesse très élevée pour la sécurité et le maintien de caisse. - l'accélération longitudinale corrigée pour l'obtention du deuxième état minimal ERM2: ce critère est utilisé pour la sécurité sur sollicitations longitudinales très importantes, au cas où la logique RoadHook ne donnerait pas satisfaction, et pour les situations d'accélérations ou de freinages stabilisés, par opposition aux phases longitudinales transitoires. - l'accélération transversale corrigée pour l'obtention du troisième état minimal ERM3: ce critère est utilisé pour la sécurité sur sollicitations latérales très importantes, au cas où la logique RoadHook ne donnerait pas satisfaction, et pour les situations de virages stabilisés, pendant lesquels la logique de synthèse privilégie la logique SkyHook. - le jerk transversal anticipé pour l'obtention du quatrième état minimal ERM4: ce critère travaille en parallèle avec la logique RoadHook avec termes anticipatifs. Il permet d'assurer une faible prise d'angle en pilotant l'actionneur par anticipation, et également, selon le paramétrage, d'utiliser des états minimaux typés en sur-vireur ou sousvireur afin de jouer sur la vivacité du véhicule lors de la mise en virage. Ces états minimaux sont par exemple calculés séparément pour chaque roue. 2890905 49 Le premier état minimal ERM1 est obtenu en extrayant de la table ou courbe préenregistrée donnant le deuxième état minimal en fonction de la vitesse du véhicule, la valeur du premier état minimal ERM1 correspondant à la vitesse VVH mesurée du véhicule, par exemple par interpolation linéaire. Le premier état minimal peut être calculé séparément pour les roues avant et pour les roues arrière. Le deuxième état minimal ERM2 est obtenu en extrayant de la table ou courbe préenregistrée donnant le deuxième état minimal en fonction de la vitesse du véhicule et de l'accélération longitudinale corrigée, la valeur du deuxième état minimal ERM2 correspondant à la vitesse VVH mesurée du véhicule et à l'accélération longitudinale corrigée XcoR, par exemple par interpolation linéaire. Le troisième état minimal ERM3 est obtenu en extrayant de la table ou courbe préenregistrée donnant le troisième état minimal en fonction de la vitesse du véhicule et de l'accélération transversale corrigée, la valeur du troisième état minimal ERM3 correspondant à la vitesse WH mesurée du véhicule et à l'accélération transversale corrigée YcoR, par exemple par interpolation linéaire. Le quatrième état minimal ERM4 est obtenu en extrayant de la table ou courbe préenregistrée donnant le quatrième état minimal en fonction du jerk transversal anticipé, la valeur du quatrième état minimal ERM4 correspondant au jerk transversal anticipé f, par exemple par interpolation linéaire. Pour chaque roue, l'état minimum ERM global d'amortissement fourni par l'estimateur 27 est alors égal au maximum des états minimaux ERM1, ERM2, ERM3, ERM4. On obtient ainsi un état minimum global ERMA, ERMB, ERMC, ERMD d'amortissement respectivement pour les roues A, B, C, D. Chacune des deux fonctions RoadHook et SkyHook a comme flux principal en entrée l'information des quatre capteurs de débattement. Par exemple, pour un véhicule roulant à moins de 20 km/h sanssollicitation conducteur, la fonction SkyHook demandera un amortissement le plus souple possible car les vitesses modales absolues seront très faibles. Cependant, dans cette situation de vie, le véhicule risque de réaliser des montées ou des descentes de trottoirs, qui sont des sollicitations contraignantes, pour lesquelles il serait préférable que le véhicule se trouve dans un état d'amortissement un peu plus ferme. De même, pour une vitesse du véhicule très élevée (sur autoroute par exemple), sans sollicitations conducteur et sur bonne route, le SkyHook demandera un état souple. Cela peut poser problème pour des grandes vitesses, car l'amortissement risque de devoir passer en très ferme en très peu de temps, ce qui n'est pas possible avec les actionneurs utilisés. D'autre part, la logique RoadHook peut se trouver légèrement en retard par rapport aux sollicitations conducteurs: les efforts anticipatifs estimés par la logique RoadHook ne sont pas en retard, mais pour qu'un passage en loi ferme soit appliqué, il faut que la roue ait déjà pris de la vitesse de débattement. Or, lorsque la roue prend de la vitesse de débattement, il est déjà trop tard. Il faut donc assurer un amortissement suffisamment ferme, indépendamment de la vitesse de débattement des roues, en intégrant des états d'amortissement minimaux sur accélérations longitudinales et latérales, ainsi que sur jerk latéral (en avance par rapport 2 0 aux accélérations). Pour améliorer le confort du véhicule, il est préférable de repasser en logique SkyHook dans les situations de vie virages stabilisés, ou accélérations longitudinales stabilisées. Cela permet de calmer les vitesses absolues de caisse. Cependant, on doit veiller dans ces situations de vie à ne pas trop sous amortir le véhicule car ces situations sont potentiellement dangereuses (virage qui se resserre, dégradation de la chaussée au cours du virage, etc...). On appliquera donc des états minimaux sur les accélérations stabilisées pour permettre d'utiliser la fonction SkyHook de manière sécurisée. Enfin, les états minimaux sur jerk permettent d'ajouter à la synthèse une marge de manoeuvre sur la vivacité et l'agrément de conduite lors des mises en virage. Commande de la loi d'amortissement Un module 28 de commande reçoit en entrée la loi ERcA, ERcB, ERcc, ERcp d'amortissement de consigne, fournie par l'estimateur 22 et l'état minimum global ERMA, ERMB, ERMC, ERMA d'amortissement, fourni par l'estimateur 27, respectivement pour les roues A, B, C, D et calcule à partir de ces états les états ERA, ERB, ERc, ERD de commande des amortisseurs des roues A, B, C, D, en prenant pour chaque roue le maximum de la loi d'amortissement de consigne et de l'état minimum global d'amortissement: ERA = max (ERcA, ERMA) ERB = max (ERcB, ERMB) ERC = max (ERcc, ERMC) ERD = max (ERco, ERMD) Ces états ERA, ERB, ERc, ERD de commande déterminent la loi d'amortissement prise par chaque amortisseur AM et sont les grandeurs ER de commande envoyées sur l'entrée COM de commande vers l'actionneur de chaque amortisseur AM pour chaque roue A, B, C, D. Les états ERA, ERB, ERc, ERD de commande sont en outre envoyés à l'entrée de l'estimateur 12 pour l'état réel ER de l'actionneur. Ci-dessous sont décrites des fonctions supplémentaires, pouvant être prévues dans le dispositif pour calculer les états ERA, ERB, ERc, ERD de commande des amortisseurs des roues A, B, C, D. Prise en compte des percussions Une détection des percussions s'effectue sur les roues avant. Il n'est pas possible d'anticiper l'obstacle. On détectera donc un obstacle lorsque les roues avant le franchiront. La détection d'une percussion s'effectue en surveillant la vitesse de débattement des roues avant du véhicule. La particularité d'une percussion est l'importante vitesse de débattement qu'elle génère au niveau des roues. L'obstacle peut être de faible amplitude (nid de poule peu profond par exemple) mais il est génère un choc, car les roues débattent à grande vitesse. A la figure 16, un estimateur 50 est prévu pour calculer un état de consigne ou loi d'amortissement de consigne ERP en cas de détection 5 d'une percussion. Cet estimateur 50 reçoit en entrée: - les débattements DEB(A), DEB(B) des roues avant A, B, fournis par les capteurs CAP-DEB de débattements, - les vitesses de débattement VDEB(A), VDEB(B) des roues avant A, B, - la vitesse WH mesurée du véhicule, - l'accélération transversale corrigée @cOR, - le coefficient a de pondération des premiers efforts Fz1, Fei et Fq,1 de confort et des deuxièmes efforts FZ2, cet et c1,2 de comportement. La détection et le traitement des percussions s'effectuent indépendamment sur les roues gauche ou droite du véhicule. Si une percussion n'est détectée que sur la roue avant droit, alors le traitement des percussions ne sera activé que sur les roues du côté droit. Si une percussion n'est détectée que sur la roue avant gauche, alors le traitement des percussions ne sera activé que sur les roues du côté gauche. L'estimateur 50 comporte: - un module 51 de détection des percussions en fonction des débattements DEB et des vitesses VDEB de débattement, - un module 52 de calcul d'un retard d'activation et d'un signal d'inhibition de traitement en fonction de la vitesse WH du véhicule, de l'accélération transversale corrigée toRet du coefficient a de pondération. - un module 53 de traitement des percussions du côté gauche, - un module 54 de traitement des percussions du côté droit. Détection des percussions Un seuil SDP de détection de percussion est prédéfini dans le module 51. Lorsque, d'un côté du véhicule, par exemple dans ce qui suit le côté gauche, la vitesse de débattement VDEB(A) de la roue avant est supérieure en valeur absolue au seuil SDP de détection de percussion, un signal logique booléen P de détection probable de percussion est mis à 1, alors que si la vitesse de débattement VDEB(A) de la roue avant est inférieure ou égale en valeur absolue au seuil SDP de détection de percussion, le signal P de détection probable de percussion est à O. Afin d'optimiser le réglage, ce seuil SDP de détection de percussion est paramétré en fonction de la vitesse WH du véhicule. On extrait de la table, courbe ou cartographie préenregistrée donnant le seuil de détection de percussion en fonction de la vitesse du véhicule, la valeur du seuil SDP de détection de percussion correspondant à la vitesse VVH du véhicule, par exemple par interpolation linéaire. A vitesses WH très élevées par exemple, n'importe quel obstacle risque de générer une vitesse de débattement importante. A vitesses de véhicule élevées, il est donc nécessaire d'augmenter le seuil SDP de détection des percussions, afin de ne pas faire de traitement intempestif sur des sollicitations route ne correspondant pas à des percussions proprement dites. Les vitesses de débattement, après un choc, peuvent osciller pendant quelques instants et risquent de dépasser plusieurs fois le seuil SDP à cause d'un seul choc de départ. Une temporisation TEMP déclenchée au premier dépassement du seuil SDP permet alors d'éviter de détecter plusieurs percussions pour un même franchissement d'obstacle. Par exemple, une percussion détectée n'est validée que si elle est détectée pendant plus d'une durée prescrite DDP de détection de percussion, par exemple de 15 millisecondes. Inhibition de la détection de percussion Un signal S = SIDP d'inhibition de détection de percussion est généré comme étant égal à 1 pour inhiber la détection de percussion, lorsque au moins l'un des débattements avant DEB(A), DEB(B) devient inférieur à un premier seuil SDEB1 de butée, ou supérieur à un deuxième seuil SDEB2 de butée, et est égal à 0 sinon. En effet, lors de forts mouvements de caisse, le débattement peut être tel que le train ira au contact de ses butées. L'écrasement des butées génère une vitesse de débattement importante, susceptible d'activer la fonction de traitement des percussions. Si cette fonction est activée dans cette situation de vie, elle imposera des états d'amortissement souple à l'arrière pendant un certain temps. Le problème est que si l'état d'amortissement passe en souple lorsque le train est en contact avec ses butées, les mouvements de caisse ne seront pas du tout freinés, des phénomènes de pompage excessif de l'essieu arrière apparaîtront. On inhibera donc la détection des percussions dans cette situation de vie. Pour cela, on surveille la valeur des débattements de roue. Lorsque ces débattements dépassent le seuil paramétrable SDEB1 ou SDEB2 (qui correspond à la course de débattement possible de la roue, préalablement au contact avec les butées d'attaque ou de détente), la détection des percussions est inhibée. Le module 51 génère un signal W de validation de percussion à partir du signal P de détection probable de percussion, de la manière suivante. Un signal Q de percussion validable et le signal W de validation de percussion sont générés au cycle n de calcul en fonction de leurs valeurs au cycle n-1 précédent et d'un signal T de temporisation TEMP écoulée, calculé en fonction du signal P de détection probable de percussion. Le signal Q de percussion validable est initialisé à 1. Un signal T de temporisation TEMP écoulée est mis à 1, si le signal P de détection probable de percussion est resté à 0 depuis son dernier front descendant pendant un temps supérieur à la temporisation TEMP. Le signal T de temporisation TEMP écoulée est sinon à O. Le signal Q de percussion validable est égal à : Q'=Q.W.T +Q.W.T + Q.W.T + Q.W.T où Q' désigne l'état au cycle suivant, et désigne le complément. Le signal W de validation de percussion est alors mis à 1, signifiant qu'une percussion est bien détectée, lorsque à la fois - le signal P de détection probable de percussion est à 1 pendant un nombre prescrit de cycles consécutifs, par exemple 3 cycles formant la 5 durée DDP, - le signal Q de percussion validable est à 1, - le signal S = SIDP d'inhibition de détection de percussion est à 0, indiquent une non inhibition, - l'accélération transversale corrigée COR est inférieure en valeur 10 absolue à un seuil prescrit SY d'inhibition de l'accélération transversale corrigée: YcoR Lorsque la percussion est détectée sur les roues avant, c'est-à-dire lorsque le signal W de validation de percussion est mis à 1, le module 52 calcule le retard DEL de franchissement des roues arrière par rapport aux roues avant d'une manière générale de la manière suivante: DEL = (e / VVH) TR où TR est un temps prescrit de réaction correspond au temps nécessaire aux actionneurs pour passer en état souple. Si la vitesse WH du véhicule est trop faible (inférieure ou égale à un seuil SWH de vitesse du véhicule) ou si le coefficient a de pondération des premiers efforts Fz1, F81 et F1 de confort et des deuxièmes efforts FZ2, cet et Traitement des roues arrière Dès que la percussion est détectée sur la roue avant gauche, une temporisation est lancée pendant le retard DEL de franchissement dans le module 53 de traitement des roues gauche. A la fin de cette temporisation, un état ERP d'amortissement de consigne souple et prescrit est imposé pour la roue arrière gauche du véhicule pendant une durée de traitement prescrite, pour que la percussion soit convenablement amortie par l'amortisseur de la roue arrière gauche. L'état d'amortissement à appliquer et la durée du traitement sont des données de mise au point paramétrables. Traitement des roues avant Dès que la percussion est détectée sur la roue avant gauche, le traitement sur la roue avant gauche ne peut être qu'un post-traitement. Celui-ci a pour objectif de limiter les ébranlements de train et de freiner les mouvements et rebonds de roue juste après l'obstacle. Le post-traitement des roues avant consiste à imposer un état ERP d'amortissement de consigne ferme et prescrit pendant une durée de posttraitement prescrite. L'état d'amortissement à appliquer et la durée du post-traitement sont des données de mise au point paramétrables. Post-traitement des roues avant et arrière A la fin du traitement des roues arrière, on effectue alors un post- traitement de la percussion, sur les roues avant et sur les roues arrière. Pour freiner les mouvements de la roue dus au passage de l'obstacle, on impose un état ERP d'amortissement de consigne ferme et prescrit pour les roues arrière pendant une durée de post-traitement prescrite. L'état d'amortissement à appliquer et la durée du post-traitement des roues avant et arrière sont des données de mise au point paramétrables. Inhibition du traitement Les modules 53, 54 de traitement des percussions produisent des états imposés ERP d'amortissement des percussions, qui peuvent être prépondérants par rapport aux états ER d'amortissement demandés par les fonctions SkyHook et RoadHook. Dans certaines situations de vie, ces états imposés ERP d'amortissement des percussions peuvent soit dégrader le confort du véhicule, soit s'avérer dangereux pour sa sécurité. C'est pourquoi le traitement des percussions est soumis à d'éventuelles inhibitions. Lorsque le véhicule se déplace sur une route très dégradée, avec des sollicitations à haute fréquence (type route pavée), les vitesses de débattement des roues atteindront des niveaux élevés, susceptibles d'activer la fonction de traitement des percussions. Si cette fonction est activée, elle imposera des états ERP de consigne d'amortissement des percussions, qui seront fermes pendant un temps déterminé sur les quatre roues. Sur une route pavée, ces états ERP d'amortissement fermes généreront de l'inconfort pendant tout le posttraitement. La stratégie idéale sur des pavés ne générant pas de mouvements de caisse est en effet de rester en loi la plus souple possible. On inhibera donc le traitement des percussions dès qu'un nombre déterminé, par exemple trois, percussions seront détectées en un court laps de temps déterminé, par exemple sur le signal W de validation de percussion. L'inhibition résultante aura une durée paramétrable. Un autre cas possible d'inhibition du traitement est les vitesses WH trop faibles du véhicule. D'autre part, lorsque la logique de synthèse de I'AMVAR se trouve en mode comportement , c'est-à-dire lorsque la logique RoadHook est activée et que le coefficient a de pondération est égal à 1 ou est proche de 1, on inhibe également le traitement des percussions (voir supra SINV). Un autre cas d'inhibition du traitement peut être prévu pour la sécurité du véhicule. Lors des fortes sollicitations de la part conducteur, ou lorsque le véhicule est installé en virage stabilisé, imposer un état d'amortissement souple peut s'avérer dangereux pour la tenue de route. Dans ces conditions de roulage, la logique RoadHook optimisant le comportement du véhicule ne doit surtout pas être désactivée par d'autres fonctions. Cela relève de la sécurité des personnes. On surveille donc d'une part l'accélération latérale du véhicule: lorsque celle-ci dépasse un certain seuil paramétrable, on inhibe le traitement des percussions, ainsi que cela est prévu ci-dessus, lorsque l'accélération transversale corrigée YcoR est supérieure ou égale en valeur absolue au seuil prescrit SY d'inhibition de l'accélération transversale corrigée: SY. YCOR Le module 52 génère un signal INHIB d'inhibition de traitement des percussions, égal à 1 pour inhiber le traitement des percussions par les modules 53 et 54, lorsque au moins l'une et/ou l'autre des conditions suivantes est réalisée: - un nombre prédéterminé de percussions, représentées par des front montant du signal W de validation de percussion est détecté en une 20 durée prédéterminée; - le signal d'inhibition SINV pour les vitesses faibles est mis à 1, pour indiquer que la vitesse WH du véhicule est trop faible ou que le coefficient a de pondération des premiers efforts Fz1, F81 et Fo de confort et des deuxièmes efforts FZ2, cet et c92 de comportement est trop grand, pour 25 indiquer que la logique Roadhook est en vigueur, YcoR >_ SY Le retard DEL de franchissement et le signal INHIB d'inhibition de traitement des percussions sont envoyées à deux entrées de chacun des modules 53, 54 de traitement. Chacun des modules 53, 54 comporte également une entrée CLK d'horloge, combinée par un opérateur logique ET avec respectivement l'entrée W(A) du signal W de validation de percussion de la roue A avant gauche et l'entrée W(B) du signal W de validation de percussion de la roue B avant droite, pour indiquer la fréquence de calcul des modules 53 et 54. Une entrée d'horloge est également prévue pour chacun des blocs, estimateurs et modules représentés aux figures. Dans le cas où l'estimateur 50 est prévu, celui-ci fournit les états de consigne ERP en cas de détection d'une percussion, c'est-à-dire pour les roues A, B, C, D, les états de consigne ERPA, ERPB, ERPc, ERPD, à une autre entrée du module 28 de commande. Le module 28 de commande calcule à partir de ces états les états ERA, ERB, ERc, ERD de commande des amortisseurs des roues A, B, C, D, en prenant pour chaque roue le maximum des états ERc, ERP d'amortissement de consigne et de l'état minimum global ERM d'amortissement: ERA = max (ERcA, ERPA, ERMA) ERB = max (ERcB, ERPB, ERMB) ERc = max (ERcc, ERPc, ERMC) ERD = max (ERCD, ERPD, ERMA) Prise en compte des grandes amplitudes de mouvement (logique des grands débattements) Une détection des grands débattements et des grandes vitesses de débattement des roues avant ou des roues arrière est prévue. L'objectif est de détecter au plus tôt les obstacles qui peuvent générer des grandes amplitudes de mouvement de caisse, en marche avant et/ou en marche arrière. La détection de ces situations de vie est prévue pour traiter les obstacles sollicitant simultanément les roues droite et gauche du train avant ou arrière. Ces obstacles peuvent être détectés en compression pour les dos d'âne ou en détente pour les saignées ou creux de taille importante. En marche avant, ce genre d'obstacles va générer sur les roues avant de fortes amplitudes de débattements et de vitesses de débattement. A la figure 17, un estimateur 60 est prévu pour calculer un état de consigne ou loi d'amortissement de consigne ERGD en cas de détection d'une grande amplitude de mouvement de roue. Cet estimateur 60 reçoit en entrée: - les débattements avant DEB(A), DEB(B) des roues avant A, B, et les débattements DEB(C), DEB(D) des roues arrière C, D, qui sont par exemple ceux filtrés DEBF(A), DEBF(B), DEBF(C), DEBF(D), fournis par le filtre 13 à partir des débattements DEB(A), DEB(B), DEB(C), DEB(D) fournis par les capteurs CAP-DEB de débattements, - les vitesses de débattement avant VDEB(A), VDEB(B) des roues avant A, B, et les vitesses de débattement DEB(C), DEB(D) des roues arrière C, D, fournies par le module dérivateur DER, - la vitesse VVH mesurée du véhicule, - le niveau NTC de tressautement de la caisse, fourni par l'estimateur 24. L'estimateur 60 met en oeuvre une logique de détection et de traitement des grandes amplitudes de mouvement et comporte: - un module 61 de détection des grandes amplitudes de mouvement de roue, - un module 62 de validation et d'inhibition de la détection des grandes amplitudes de mouvement de roue, - un module 63 de calcul d'un coefficient x de traitement des grandes 25 amplitudes de mouvement de roue, - un module 64 de calcul de l'état de consigne ou loi d'amortissement de consigne ERGD pour les grandes amplitudes de mouvement de roue. Détection des grandes amplitudes de mouvement de roue Un premier seuil SDGD de détection de grands débattements et un deuxième seuil SVGD de détection de grandes vitesses de débattement sont prédéfinis dans le module 61. Lorsqu'à la fois le débattement DEBF(A) de la roue avant gauche dépasse le premier seuil SDGD de détection de grands débattements, le débattement DEBF(B) de la roue avant droite dépasse le premier seuil SDGD de détection de grands débattements, la vitesse de débattement VDEB(A) de la roue avant gauche dépasse le deuxième seuil SVGD de détection de grandes vitesses de débattement et la vitesse de débattement VDEB(B) de la roue avant droite dépasse le deuxième seuil SVGD de détection de grandes vitesses de débattement, alors un premier signal SDGDAV de détection de grands mouvements avant est mis à 1 pour indiquer la détection d'une grande amplitude de mouvement de roue sur les roues avant. II en est de même pour un deuxième signal SDGDAR de détection de grands mouvements arrière, qui est mis à 1 pour indiquer la détection d'une grande amplitude de mouvement de roue sur les roues arrière, lorsque les quatre conditions de dépassement à la fois sont remplies avec les débattements DEBF(D) et DEBF(C) et les vitesses de débattement VDEB(D) et VDEB(C) pour les roues arrière. Les premier et deuxièmes seuils SDGD et SVGD peuvent être différents pour l'avant et l'arrière. Le dépassement des premier et/ou deuxièmes seuils SDGD, SVGD peut être le passage du débattement et/ou de la vitesse de débattement sous le seuil SDGD, SVGD inférieur, par exemple en détente des amortisseurs, et/ou le passage du débattement et/ou de la vitesse de débattement au-dessus d'un autre seuil SDGD supérieur au seuil SDGD, SVGD inférieur, par exemple en attaque des amortisseurs. Un signal SGD de détection de grands mouvements est mis à 1 pour indiquer la détection d'une grande amplitude de mouvement de roue sur les roues, lorsque le premier signal SDGDAV de détection de grands mouvements avant et/ou le deuxième signal SDGDAR de détection de grands mouvements arrière détecte à 1. Le signal SGD de détection de grands mouvements est fourni par le module 61 de détection au module 62 de validation et d'inhibition. Pour plus de précision et pour éviter des traitements intempestifs, le premier seuil SDGD de détection de grands débattements et le deuxième seuil SVGD de détection de grandes vitesses de débattement sont paramétrés en fonction de la vitesse VVH du véhicule. Par exemple, pour chacun de ces seuils SDGD, SVGD, on extrait de la table, courbe ou cartographie préenregistrée donnant le seuil de détection en fonction de la vitesse du véhicule, la valeur du seuil SDGD, SVGD de détection correspondant à la vitesse VVH du véhicule, par exemple par interpolation linéaire. Inhibition de la détection des grands mouvements de roue Un signal INSGD de validation ou d'inhibition de détection de grande amplitude de mouvement de roues est généré par le module 62, comme étant égal à 0 pour inhiber la détection, lorsque l'une ou plusieurs des conditions suivantes sont réalisées: - le coefficient a de pondération des premiers efforts Fz1, F01 et Fo de confort et des deuxièmes efforts FZ2, c82 et ce de comportement est trop grand, (supérieur à un seuil SCOEFF2 de coefficient de pondération, par exemple nul) pour indiquer que la logique Roadhook est au moins partiellement en vigueur, - le niveau NTC de tressautement est supérieur à un seuil SNTC prescrit de niveau de tressautement. Si aucune des conditions d'inhibition n'est réalisée et si le signal SGD de détection de grands mouvements est 1 pour indiquer la détection d'une grande amplitude de mouvement des roues, le signal INSGD prend l'état 1 de validation de la détection de grande amplitude de mouvement de roues. Dans le premier cas d'inhibition (coefficient a de pondération), sous sollicitations du conducteur, il est plus sûr de laisser la logique Roadhook agir et réagir aux sollicitations routières pour améliorer la tenue de caisse et notamment pour maximiser le contact de la roue avec le sol. Si la logique Roadhook veut transmettre une consigne de passage en états d'amortissement souple, il ne faut pas l'en empêcher. C'est pourquoi on inhibe la détection et le traitement des grandes amplitudes de mouvement lorsque la logique Roadhook est active. Dans le deuxième cas d'inhibition (niveau NTC de tressautement), le traitement des grandes amplitudes de mouvement peut être pénalisant pour le confort vibratoire, car un état d'amortissement trop ferme transférera dans la caisse les irrégularités de la route, et ne filtrera pas donc pas les tressautements et les trépidations engendrées par cette route. C'est pourquoi il est préférable d'inhiber la logique de traitement des grandes amplitudes de mouvement, lorsque la route est dégradée. On utilise une logique de reconnaissance de l'état de la route, basée sur des filtrages passe-bande des débattements DEB. Ainsi que cela est indiqué ci-dessus pour le calcul du niveau NMC de basses fréquences et du niveau NTC de tressautement, un filtrage autour du mode de caisse (autour du Hertz) et un filtrage dans la bande des tressautements (entre 3 et 8 Hz) sont utilisés pour caractériser l'état de la route (bonne route, route à bon revêtement mais générant des mouvements de caisse, route à revêtement dégradé mais plate, route à revêtement dégradé mais générant des mouvements de caisse). Pour l'inhibition, on utilise le niveau de tressautement calculé à partir du filtrage entre 3 et 8 Hz. Le seuil SNTC prescrit de niveau de tressautement est paramétrable. Ainsi, le compromis entre tenue de caisse et confort vibratoire est optimisé. Traitement des grands mouvements de roues L'estimateur 63 calcule, à partir du signal INSGD de validation ou d'inhibition de détection de grande amplitude de mouvement de roues, le coefficient x de traitement des grandes amplitudes de mouvement de roue. Le coefficient x de traitement est une variable supérieure ou égale à 0 et inférieure ou égale à 1. Par défaut, le coefficient x de traitement est à 0. Lorsque le signal INSGD passe de l'état 0 d'inhibition de détection de grande amplitude de mouvement de roues à l'état 1 de validation de grande amplitude de mouvement de roues, le coefficient x de traitement augmente de 0 à 1 avec une pente de montée prescrite, par exemple paramétrable par une première temporisation TEMP1 en entrée du module 63. Le coefficient x de traitement est ensuite maintenu à sa valeur maximum 1 pendant une durée prescrite, par exemple paramétrable par une deuxième temporisation TEMP2 en entrée du module 63, et redescend à 0 avec une pente de descente prescrite, par exemple paramétrable par une troisième temporisation TEMP3 en entrée du module 63. Etats minimaux en cas de détection des grands mouvements de roues Le module 64 reçoit le coefficient x de traitement des grandes amplitudes de mouvement de roue et la vitesse VVH du véhicule et calcule en fonction de ceux-ci la loi d'amortissement de consigne ERGD en cas de détection d'une grande amplitude de mouvement de roue. Le traitement des situations de grande amplitude de mouvement de roues est effectué à l'aide d'états ERGD d'amortissement de consigne minimaux. Les différents paramètres intervenant dans le calcul du coefficient x de traitement permettent de maîtriser exactement l'instant et le tempspendant lequel les états ERGD d'amortissement minimaux seront appliqués par le module 64. Ces états minimaux ERGD sont paramétrables en fonction de la vitesse WH du véhicule pour optimiser le compromis entre tenue de caisse et confort vibratoire quelle que soit la vitesse du véhicule: les états minimaux à utiliser sont par exemple moins élevés à 30 km/h pour le passage des dos d'âne qu'à vitesse plus élevée où une sollicitation de la route créant un grand débattement va nécessiter des états minimaux élevés. Les états minimaux ERGD peuvent également être calculés séparément pour les roues avant et pour les roues arrière. Le calcul des états ERGD d'amortissement de consigne est par exemple effectué de la manière suivante: - un état intermédiaire ERGD-INTER de grande amplitude de mouvement de roues (numéro de loi d'amortissement intermédiaire) est extrait d'une table ou courbe préenregistrée donnant cet état intermédiaire état intermédiaire en fonction de la vitesse du véhicule, la valeur ERGD- INTER de l'état intermédiaire de grande amplitude de mouvement de roues, correspondant à la vitesse WH du véhicule, par exemple par interpolation linéaire. - l'état ERGD d'amortissement de consigne de grande amplitude de mouvement de roues est alors égal à l'état intermédiaire ERGD-INTER d'amortissement, multiplié par le coefficient x de traitement des grandes amplitudes de mouvement de roue, arrondi par exemple au numéro de loi d'amortissement le plus proche. Dans le cas où l'estimateur 60 est prévu, celui-ci fournit les états d'amortissement de consigne ERGD en cas de détection d'une grande amplitude de mouvement de roue, c'est-à-dire pour les roues A, B, C, D, les états de consigne ERPA, ERPB, ERPc, ERPD, à une autre entrée du module 28 de commande. Le module 28 de commande calcule à partir de ces états les états ERGDA, ERGDB, ERGDc, ERGDD de commande des amortisseurs des roues A, B, C, D, en prenant pour chaque roue le maximum des états ERS, ERGD (et éventuellement ERP pour la prise en compte des percussions) d'amortissement de consigne et de l'état minimum global ERM d'amortissement: ERA = max (ERcA, ERGDA, ERMA) ERB = max (ERCB, ERGDB, ERMB) ERc = max (ERcc, ERGDc, ERMC) ERD = max (ERcg, ERGDD, ERMA) 2890905 66 REVENDICATIONS 1. Dispositif de commande d'une suspension d'une caisse (2) de véhicule automobile sur ses roues (A), comportant un calculateur (CSS) apte à calculer une grandeur (ER) de commande d'un actionneur (M) d'au moins un amortisseur (AM) variable de la suspension (S) en fonction d'au moins un effort modal de consigne de la caisse, caractérisé en ce que le calculateur (CSS) comporte: - un premier moyen (31) de calcul d'au moins une vitesse modale (Vmod2) relative de caisse par rapport au plan moyen des roues, - un deuxième moyen (34, 22) de calcul de l'effort (cq)2, cet, Fz2) 10 modal de consigne de l'amortisseur, en fonction de la vitesse modale relative (Vmod2) de caisse par rapport au plan moyen des roues. 2. Dispositif de commande d'une suspension suivant la revendication 1, caractérisé en ce qu'il comprend au moins un capteur (CAP-DEB) de débattement d'une roue (A, B, C, D) par rapport à la caisse (2) du véhicule, le premier moyen (31) de calcul de vitesse modale (Vmod2) relative de caisse étant apte à calculer ladite vitesse modale (Vmod2) relative de caisse par rapport au plan moyen des roues en fonction au moins du débattement (DEB) fourni par le capteur (CAP-DEB) de débattement. 3. Dispositif de commande d'une suspension suivant l'une quelconque des revendications précédentes, caractérisé en ce que le premier moyen (31) de calcul de vitesse modale (Vmod2) relative de caisse comprend un filtre passe-bas du débattement (DEB) fourni par le capteur (CAP-DEB) de débattement et un dérivateur pour l'obtention de ladite vitesse modale (Vmod2) relative de caisse par rapport au plan moyen des roues. 4. Dispositif de commande d'une suspension suivant l'une quelconque des revendications précédentes, caractérisé en ce que le mode est au moins l'un parmi le pompage vertical (z), le roulis (0) et le tangage (O. 5. Dispositif de commande d'une suspension suivant la revendication 4, caractérisé en ce qu'un capteur (CAP-DEB) de débattement est prévu pour chaque roue (A, B, C, D), et le premier moyen (31) de calcul de vitesse modale (Vmod2) relative de caisse est apte à calculer au moins une des vitesses modales relatives (Vmod2) suivantes: - la vitesse modale relative G2 de caisse en pompage par rapport au plan moyen des roues selon la formule: _ (e -1g) dA+dB lg d+dD Zc2 + e 2 e 2 - la vitesse modale relative cp2 de caisse en tangage par rapport au 10 plan moyen des roues selon la formule: dA + dB - dc - dD (P2 - 2e - la vitesse modale relative 02 de caisse en roulis par rapport au plan moyen des roues selon la formule: e2-dA-dB-dC+dD 2v où dA est la vitesse de débattement de la roue (A) avant gauche, dB est la vitesse de débattement de la roue (B) avant droite, di est la vitesse de débattement de la roue (C) arrière droite, dD est la vitesse de débattement de la roue (D) arrière gauche, e est l'empattement prédéterminé du véhicule, v est la voie prédéterminée du véhicule, et Ig est une longueur prédéterminée séparant le centre de gravité (G) du véhicule de l'axe des roues avant (A, B). 6. Dispositif de commande d'une suspension suivant la revendication 5, caractérisé en ce que le calculateur (CSS) comprend un moyen (20) de calcul de la longueur (Ig) séparant le centre de gravité (G) de l'axe des roues avant (A, B) par la formule Ig = (1 RMAvAr).e où RMAvAr est une valeur de répartition de masse entre l'avant et l'arrière. 7. Dispositif de commande d'une suspension suivant la revendication 6, caractérisé en ce qu'il comprend un moyen de mesure de la vitesse (WH) du véhicule, le moyen (20) de calcul de la longueur (Ig) séparant le centre de gravité (G) de l'axe des roues avant (A, B) étant apte à calculer la valeur de répartition de masse (RMAvAr) entre l'avant et l'arrière en fonction des débattements fournis par les capteurs (CAP-DEB) de débattement et de la vitesse (WH) du véhicule. 8. Dispositif de commande d'une suspension suivant l'une quelconque des revendications précédentes, caractérisé en ce que le 15 deuxième moyen (34, 22) de calcul d'effort comporte: un premier moyen (34) de calcul d'un premier terme d'effort modal anticipatif en fonction d'un jerk de caisse déterminé sur le véhicule, un deuxième moyen (34) de calcul d'un deuxième terme d'effort modal correctif en fonction de ladite vitesse modale (Vmod2) relative de caisse par rapport au plan moyen des roues, un troisième moyen (34) de calcul de l'effort modal (cq)2, c82, Fer) de consigne en fonction du premier terme d'effort modal anticipatif et du deuxième terme d'effort modal correctif. 9. Dispositif de commande suivant la revendication 8, caractérisé en ce que le troisième moyen de calcul de l'effort modal (cq)2, Cet, FZ2) de consigne est apte à calculer l'effort modal (cT2, ce2, FZ2) de consigne par sélection - du deuxième terme d'effort modal correctif, lorsque la valeur absolue du premier terme d'effort modal anticipatif est inférieure ou égale à 30 une première valeur prescrite non nulle, - du premier terme d'effort modal anticipatif, lorsque la valeur absolue du premier terme d'effort modal anticipatif est supérieure à la première valeur prescrite et lorsque la valeur absolue du deuxième terme d'effort modal correctif est inférieure ou égale à une deuxième valeur prescrite non nulle, lorsque la valeur absolue du premier terme d'effort modal anticipatif est supérieure à la première valeur prescrite et lorsque la valeur absolue du deuxième terme d'effort modal correctif est supérieure à la deuxième valeur prescrite, - - de celui des premier terme d'effort modal anticipatif et deuxième terme d'effort modal correctif, ayant la valeur absolue maximum, lorsque le premier terme d'effort modal anticipatif et le deuxième terme d'effort modal correctif sont de même signe, - - du deuxième terme d'effort modal correctif lorsque le premier terme d'effort modal anticipatif et le deuxième terme d'effort modal correctif sont de signes opposés. 10. Dispositif de commande d'une suspension suivant l'une quelconque des revendications 8 et 9, caractérisé en ce qu'un capteur (CAP-DEB) de débattement est prévu pour chaque roue (A, B, C, D), le moyen de calcul du deuxième terme d'effort modal correctif est apte à calculer en fonction au moins des mesures de débattement des roues fournies par lesdits capteurs (CAP-DEB) de débattement, un deuxième gain modal correctif bmod2 d'amortissement variable pour le calcul d'un deuxième terme d'effort modal correctif F2COR selon la formule: F2COR = - bmod2 É Vmod2 11. Dispositif de commande d'une suspension suivant la 25 revendication 10, caractérisé en ce que le moyen de calcul du deuxième terme d'effort modal correctif est apte à calculer le deuxième gain modal correctif (bmod2) d'amortissement variable, de manière à ce que ce deuxième gain modal correctif (bmod2) d'amortissement variable comprenne au moins un coefficient multiplicatif de type de conduite (bzTYP, beryp, bo-yp), qui est égal à un gain préenregistré de sportivité, dans le cas où un bouton de commande de conduite sportive sur le tableau de bord du véhicule est dans une position de conduite sportive, et qui est égal à l'unité, dans le cas où le bouton de commande de conduite sportive est dans une position d'absence de conduite sportive. 12. Dispositif de commande d'une suspension suivant l'une quelconque des revendications 10 et 11, caractérisé en ce qu'il comprend un moyen de mesure de la vitesse (VVH) du véhicule, et le moyen de calcul du deuxième terme d'effort modal correctif est apte à calculer le deuxième gain modal correctif (bmod2) d'amortissement variable, de manière à ce que ce deuxième gain modal correctif (bmod2) d'amortissement variable comprenne au moins un coefficient multiplicatif de référence (bZREF, bOREF, bcpREF), calculé en fonction de la vitesse mesurée (WH) du véhicule. 13. Dispositif de commande d'une suspension suivant l'une quelconque des revendications 10 à 12, caractérisé en ce que le moyen de calcul du deuxième terme d'effort modal correctif est apte à calculer le deuxième gain modal correctif (bmod2) d'amortissement variable, de manière à ce que ce deuxième gain modal correctif (bmod2) d'amortissement variable comprenne au moins un coefficient multiplicatif de recalage (bzREC, bOREC, bcpREC), le coefficient multiplicatif de recalage (bZREC) en pompage est égal à _ kZ É MSUS bZREC kZREF MREF ' le coefficient multiplicatif de recalage (bZREC) en roulis est égal à beREC = , I y le "9REF IeRFF le coefficient multiplicatif de recalage (b(pREC) en tangage est égal à b(pREC = k(pREF Iq)REF où kZF est une raideur de référence en pompage, constante, est une raideur de référence en roulis, constante, est une raideur de référence en tangage, constante, IeREF est un moment d'inertie de référence en roulis, constant, kp Iw keREF kcpREF REF est un moment d'inertie de référence en tangage, constant, MREF est une masse de référence, constante, le est un moment d'inertie en roulis, calculé en fonction au moins des débattements, I est un moment d'inertie en tangage calculé en fonction au moins des débattements, MSUS est une masse suspendue du véhicule, calculée en fonction au moins des débattements, kZ est une raideur modale en pompage, calculée en fonction au moins des 10 débattements, ki, est une raideur modale en tangage, calculée en fonction au moins des débattements, ke est une raideur modale en roulis, calculée en fonction au moins des débattements. 14. Dispositif de commande d'une suspension suivant la revendication 7 ou 13, caractérisé en ce qu'il comprend: - un moyen de mesure de la vitesse (WH) du véhicule, - un moyen (20) de calcul d'une assiette statique avant (ASav) et d'une assiette statique arrière (ASar) du véhicule en fonction des mesures (DEBF) de débattement, fournies par les capteurs de débattement, - un moyen (20) de calcul d'une masse dynamique apparente avant (MDAAV) et d'une masse dynamique apparente arrière (MDAAR), en fonction des mesures de débattement (DEB), fournies par les capteurs (CAPDEB) de débattement, - un moyen (20) de calcul d'un biais aérodynamique avant (BAAV) et d'un biais aérodynamique arrière (BAAR), en fonction de la vitesse du véhicule (WH), - un moyen (20) de calcul de la masse suspendue (MSUS) du véhicule et d'une valeur de répartition de masse (RMAvAr) entre l'avant et 30 l'arrière du véhicule, en fonction de la masse dynamique apparente avant 2890905 72 (MDAAV), de la masse dynamique apparente arrière (MDAAR), du biais aérodynamique avant (BAAV) et du biais aérodynamique arrière (BAAR). 15. Dispositif de commande d'une suspension suivant la revendication 14, caractérisé en ce que le moyen (20) de calcul de l'assiette statique avant (ASav) et de l'assiette statique arrière (ASar) du véhicule est apte à calculer l'assiette statique avant (ASav), respectivement arrière (ASar), comme étant le débattement moyen (DEBAVMOY) des débattements (DEB) des roues avant (A, B), respectivement arrière (C, D), filtré par un filtre passe-bas, débattement moyen filtré auquel est ensuite rajoutée une constante de décalage d'assiette avant, respectivement arrière. 16. Dispositif de commande d'une suspension suivant l'une quelconque des revendications 14 et 15, caractérisé en ce que le moyen (20) de calcul des masses dynamiques apparentes avant 15 et arrière (MDAAV, MDAAR) comprend des moyens pour - calculer le débattement relatif avant, respectivement arrière, égal au débattement moyen des débattements (DEB) des roues avant (A, B), respectivement arrière (C, D), auquel est ensuite rajoutée une constante de décalage avant, respectivement arrière, - extraire un effort dynamique de flexion de ressort avant EDFAV, respectivement arrière EDFAR, à partir d'une table ou courbe enregistrée donnant cet effort EDFAV en fonction du débattement relatif avant, - calculer la masse dynamique apparente avant MDAAV, respectivement arrière MDAAR, par la formule: MDAAV = (EDFAV.2/g) + constante avant, MDAAR = (EDFAR.2/g) + constante arrière, où g est la constante d'accélération de la pesanteur = 9,81 m.s 2. 17. Dispositif de commande d'une suspension suivant l'une quelconque des revendications 14 à 16, caractérisé en ce que le moyen (20) de calcul de la masse suspendue (MSUS) du véhicule et de la valeur de répartition de masse (RMAvAr) entre l'avant et l'arrière du véhicule comprend au moins un filtre passe-bas (PB1, PB2) de la somme de la masse dynamique apparente avant (MDAAV), respectivement arrière (MDAAR), et du biais aérodynamique avant (BAAV), respectivement arrière (BAAR), pour obtenir une masse suspendue avant (MSUSAV), respectivement arrière (MSUSAR), la masse suspendue (MSUS) étant égale à la somme de la masse suspendue avant (MSUSAV) et de la masse suspendue arrière (MSUSAR), la valeur de répartition de masse (RMAvAr) entre l'avant et l'arrière du véhicule étant égale à la masse suspendue avant (MSUSAV) divisée par la masse suspendue du véhicule (MSUS). 18. Dispositif de commande d'une suspension suivant la revendication 17, caractérisé en ce qu'il comprend un moyen d'entrée d'une information (10) d'ouvrant, le moyen de calcul de la masse suspendue (MSUS) du véhicule et de la valeur de répartition de masse (RMAvAr) entre l'avant et l'arrière est apte à enregistrer des valeurs successives de la masse suspendue avant (MSUSAV), respectivement arrière (MSUSAR), et à retenir la valeur précédemment enregistrée de la masse suspendue avant (MSUSAV), respectivement arrière (MSUSAR), au lieu de sa valeur suivante, lorsque au moins l'une des conditions suivantes est réalisée: - la vitesse (VVH) du véhicule n'est pas comprise entre un seuil bas prédéterminé (WH1) et un seuil haut prédéterminé (WH2), - l'information d'ouvrant (IO) indique une ouverture de cet ouvrant du véhicule et la vitesse (WH) du véhicule n'est pas supérieure à un seuil prescrit (VVH3), - l'écart entre ladite valeur précédemment enregistrée (MSUSEAVF(n-1)) de la masse suspendue avant (MSUSAV), respectivement arrière (MSUSAR), et sa valeur suivante (MSUSEAVF(n)), est inférieur ou égal en valeur absolue à un écart prescrit (A). 19. Dispositif de commande d'une suspension suivant l'une 30 quelconque des revendications 14 à 18, caractérisé en ce qu'il comprend: -un moyen (20) de calcul du moment d'inertie en roulis (le) et du moment d'inertie en tangage (19) en fonction de la masse suspendue (MSUS) et de la valeur de répartition de masse (RMAvAr) entre l'avant et l'arrière du véhicule, un moyen (20) de calcul de la longueur (Ig) séparant le centre de gravité (G) de l'axe des roues avant (A, B), - un moyen (20) de calcul d'une raideur modale en pompage (kr), d'une raideur modale en tangage (ky) et d'une raideur modale en roulis (ke), en fonction de l'assiette statique (AS) et de la valeur de répartition de masse (RMAvAr) entre l'avant et l'arrière. 20. Dispositif de commande d'une suspension suivant lâ 10 revendication 19, caractérisé en ce que ledit moyen (20) de calcul est apte à calculer: - le moment d'inertie en roulis le par la formule le = Ay.MSUSAR + By Où Ay et By sont des paramètres prédéterminés, MSUS est la masse suspendue du véhicule, RMAvAr est la valeur de répartition de masse entre l'avant et l'arrière, MSUSAR est la masse suspendue arrière du véhicule, égale à MSUSAR = (1-RMArAv).MSUS - le moment d'inertie en tangage I9 par la formule I = AX.MSUS + BX Où AX et BX sont des paramètres prédéterminés. 21. Dispositif de commande d'une suspension suivant l'une quelconque des revendications 19 et 20, caractérisé en ce que le moyen (20) de calcul de la longueur (Ig) séparant le centre de gravité (G) de l'axe des roues avant (A, B) est apte à calculer ladite longueur Ig séparant le centre de gravité (G) de l'axe des roues avant (A, B) par la formule Ig = (1 RMAvAr).e où RMAvAr est la valeur de répartition de masse entre l'avant et l'arrière, e est l'empattement prédéterminé du véhicule. 22. Dispositif de commande d'une suspension suivant l'une quelconque des revendications 19 à 21, caractérisé en ce que le moyen (20) de calcul des raideurs modales est apte à calculer - la raideur modale en pompage kZ, comme étant la somme d'une 5 raideur de suspension avant kAV et d'une raideur de suspension arrière kAR - la raideur modale en tangage Kp par la formule kp = kAV.(Ig)2 + kAR.(e-lg)2 - la raideur modale en roulis ke par la formule ke = Kbadav + Kbadar + v2.(kAV+kAR)/4 où la raideur de suspension avant kAV, respectivement arrière kAR sont obtenues par extraction d'une table ou courbe préenregistrée donnant la raideur de la suspension avant, respectivement arrière, en fonction de l'assiette statique avant, respectivement arrière, de la valeur de raideur avant, respectivement arrière, correspondant à l'assiette statique avant, respectivement arrière calculées, où Kbadav est un paramètre prédéterminé correspondant à la raideur d'une barre anti-devers avant du véhicule, Kbadar est un paramètre prédéterminé correspondant à la raideur 20 d'une barre anti-devers arrière du véhicule. 23. Dispositif de commande suivant l'une quelconque des revendications 8 à 22, lorsqu'elles dépendent au moins de la revendication 8, caractérisé en ce qu'il comprend un moyen (33) de calcul d'un jerk longitudinal anticipé (X) et un moyen (32) de calcul d'un jerk transversal 25 anticipé (@ ), le premier moyen (34) de calcul du premier terme d'effort modal anticipatif étant apte à calculer un couple modal anticipatif en tangage (Ccp2ant) en fonction du jerk longitudinal anticipé (X) et un couple modal anticipatif en roulis (celant) en fonction du jerk transversal anticipé (@ ). 24. Dispositif de commande d'une suspension suivant la revendication 23, caractérisé en ce que le premier moyen (34) de calcul du premier terme d'effort modal anticipatif comprend - un filtre (341) d'annulation des faibles amplitudes du jerk longitudinal anticipé (X), respectivement du jerk transversal anticipé (@ ), - un module (342) de maintien des maxima du jerk longitudinal anticipé (X) filtré, respectivement du jerk transversal anticipé (@) filtré, - un module (343) limiteur de pente du jerk respectif issu du module (342) de maintien des maxima, pour obtenir un jerk longitudinal anticipé (Z) traité, respectivement un jerk transversal anticipé (@T) traité, qui sont multipliés respectivement par un gain en sollicitation longitudinale (Gsx) pour obtenir le couple modal anticipatif en tangage (C(p2ant) et par un gain en sollicitation longitudinale (Gsy) pour obtenir le couple modal anticipatif en roulis (Ce2ant)É 25. Dispositif de commande d'une suspension suivant l'une quelconque des revendications 14 à 22 et l'une quelconque des revendications 23 et 24, caractérisé en ce qu'il comprend un moyen de mesure d'une pression de maître cylindre de freinage, le moyen (33) de calcul du jerk longitudinal anticipé (X) étant apte à 20 calculer le jerk longitudinal anticipé (X) en fonction de la pression de maître cylindre de freinage, d'un couple moteur anticipé aux roues et de la masse suspendue (MSUS). 26. Dispositif de commande d'une suspension suivant la revendication 25, caractérisé en ce que le moyen (33) de calcul du jerk longitudinal anticipé (X) est apte à extraire d'une table ou courbe préenregistrée donnant un effort de freinage du maître cylindre en fonction de la pression du maître cylindre, la valeur (EFR) de cet effort de freinage correspondant à la pression (PME) du maître cylindre, le moyen (33) de calcul du jerk longitudinal anticipé ( ) comprend 30 un filtre passe-bas de l'effort (EFR) de freinage, et un dérivateur de l'effort (EFR) de freinage ainsi filtré pour obtenir la dérivée ÉFRF de l'effort (EFR) de freinage filtré, le moyen (33) de calcul du jerk longitudinal anticipé (X) est apte à calculer un effort anticipé du moteur aux roues (EMR), égal au couple moteur anticipé aux roues (CR) divisé par un rayon moyen (Rmoy) préenregistré des roues, le moyen (33) de calcul du jerk longitudinal anticipé (X) comprend un filtre passe-bas de cet effort (EMR) anticipé du moteur aux roues, et un dérivateur de l'effort (EMR) anticipé du moteur ainsi filtré pour obtenir la dérivée EMRF de l'effort (EMR) du moteur filtré, le jerk longitudinal anticipé X étant égal à X = EFxF + EMnF MTOT où MTOT est la masse totale MTOT du véhicule, prédéterminée. 27. Dispositif de commande d'une suspension suivant la revendication 26, caractérisé en ce qu'il comprend un moyen de mesure de la vitesse (VVH) du véhicule, un moyen de calcul de la vitesse wMOT de rotation du moteur à la vitesse (VVH) du véhicule et un moyen (40) de calcul du couple moteur anticipé aux roues CR selon la formule: CR = CM. REMBR(i), avec REMBR(i) = wMOT/wROUE Où REMBR(i) est le rapport d'embrayage ayant le numéro i, WROUE est la vitesse de rotation des roues, CM est un couple moteur prédéterminé. 28. Dispositif de commande d'une suspension suivant l'une quelconque des revendications 23 à 27, caractérisé en ce que le moyen (32) de calcul d'un jerk transversal anticipé (@) est apte à calculer le jerk transversal anticipé f en fonction de la masse suspendue (MSUS), de la valeur de répartition de masse (RMAvAr) entre l'avant et l'arrière, de la vitesse WH du véhicule, et de la vitesse de rotation ô du volant de direction selon la formule: D.g..VVH2 Y_ e(1 + K.VVH2) où D est la démultiplication prédéterminée du volant de direction du véhicule, e est l'empattement prédéterminé du véhicule et K est une constante de gain de surbraquage, calculée en fonction de la valeur de répartition de masse (RMAvAr) entre l'avant et l'arrière et de la masse suspendue (MSUS). 29. Dispositif de commande d'une suspension suivant la revendication précédentes, caractérisé en ce qu'il comprend au moins pour une roue un capteur (CAP-DEB) de débattement de la roue par rapport à la caisse (2), un moyen (10) de calcul de la vitesse (VDEB) de débattement de l'amortisseur en fonction du débattement de la roue correspondante, un moyen (22) de calcul d'une force (FA1) de consigne de l'amortisseur (AM) en fonction dudit au moins un effort modal de consigne, et un moyen (22) de calcul de la grandeur (ER) de commande de l'actionneur (M) de l'amortisseur en fonction de cette force (FA1) de consigne et de la vitesse (VDEB) de débattement de la roue. 30. Dispositif de commande d'une suspension suivant l'une quelconque des revendications précédentes, caractérisé en ce que la grandeur (ER) de commande est une loi d'amortissement déterminée parmi une pluralité de lois (ER) d'amortissement différentes imposant la force (FA) de l'amortisseur en fonction de sa vitesse (VDEB) de débattement. 31. Véhicule automobile (1) comportant une caisse (2), des roues (A, B, C, D), une suspension (S) de la caisse (2) sur les roues (A, B, C, D) et un dispositif de commande de la suspension (S) suivant l'une quelconque des revendications précédentes. 32. Procédé d'obtention d'un véhicule automobile, le véhicule automobile étant muni de roues, d'une caisse, d'une suspension ayant au moins un amortisseur à amortissement variable de la caisse sur les roues, et d'un dispositif de commande de la suspension, le dispositif de commande ayant au moins un calculateur (CSS) apte à 2890905 79 calculer une grandeùr (ER) de commande d'un actionneur (M) dudit au moins un amortisseur (AM) de la suspension, le procédé d'obtention comportant une étape de montage du calculateur (CSS) sur le véhicule, caractérisé en ce que le procédé d'obtention comporte au moins une étape de programmation du calculateur suivant au moins un programme comportant des instructions de programme mettant en oeuvre les moyens de calcul du dispositif de commande de la suspension suivant l'une quelconque des revendications 1 à 30. 33. Programme informatique de pilotage d'un calculateur (CSS), comportant des instructions de programme pour le calcul d'au moins une vitesse modale (Vmod2) relative de caisse par rapport au plan moyen des roues, pour le calcul de l'effort (c(,2, ce2, Fz2) modal de consigne de l'amortisseur, en fonction de la vitesse modale relative (Vmod2) de caisse par rapport au plan moyen des roues et pour le calcul de la grandeur (ER) de commande d'un actionneur (M) d'au moins un amortisseur (AM) variable de la suspension (S) en fonction d'au moins l'effort modal de consigne de la caisse, lorsqu'il est mis en oeuvre dans un dispositif de commande de suspension selon l'une quelconque des revendications 1 à 30.