E Print 9677
E Print 9677
MÉMOIRE PRÉSENTÉ À
L'UNIVERSITÉ DU QUÉBEC À TROIS-RIVIÈRES
PAR
NOFINIDY HO FANILO ANDRIAMINAHY
AVRIL 2021
Université du Québec à Trois-Rivières
Service de la bibliothèque
Avertissement
Résumé
Les nouvelles normes pour réduire la pollution et les émissions de gaz à effet de serre pour
les opérateurs dans les industries automobiles ont utilisé différentes technologies permettant
solution efficace et fiable pour économiser du carburant et aussi pour diminuer la pollution
de l'environnement. Bien que cela soit une issue prometteuse, l'électrification totale du
batterie qui est la seule source d'énergie des véhicules électriques rechargeables à freinage
externes dont le conducteur, les conditions de la route, les conditions météorologiques, cela
impliquerait de connaitre comment l'énergie est dispersée afin de pouvoir proposer une
Cette étude a pour but d'estimer l'autonomie d'un véhicule électrique en se basant sur
l'apport de chaque force extérieure sur l'efficacité énergétique. Afin d'évaluer l'effet de
chaque amélioration sur l'efficacité énergétique, il faut être en mesure de déterminer avec
précision la consommation énergétique en route. Le défi actuel réside dans l'absence d'un
système de mesure capable de fournir cette information. En effet, les différents capteurs
d'énergie, principalement pour des fins de normes environnementales. Il est donc important
d'estimation correcte des pertes aérodynamiques, des pertes par friction roues/routes et des
pertes dues à la configuration de la route ainsi que la puissance électrique fournie par la
batterie du véhicule. Dans ce projet, nous avons utilisé comme méthodes les moindres carrés
récursifs ou Recursive Least Square (RLS) en anglais et le filtre de Kalman ou Filter Kalman
Les résultats ont montré que le KF donne une bonne performance que le RLS. Le modèle
NARX ou Nonlinear AutoRegressive eXogenous est utilisé pour estimer la puissance. Pour
estimer l'autonomie, on s'est basé sur le même modèle NARX avec d'autres calculs basés
sur des mises à jour de données. Les simulations faites ont montré que le modèle NARX
La validation de ces estimations a été réalisé sur un véhicule électrique à freinage régénératif :
KlA SOUL EV, appartenant à l'Institut de Recherche sur l'Hydrogène (IRH) de l'Université
Mots clés: force aérodynamique, résistance au roulement, moindres carrés récursifs, filtre
Dédicaces
Je dédie ce mémoire :
Ji ma très chère femme Francia Ravalison et mon fils Ethan Andriaminahy pour leur
soutien, leur amour. Merci de toujours croire en moi, vous êtes ma source de motivation.
Ji mes chers parents Laurent et Noro Andriaminahy pour tous leurs sacrifices et leur
Ji mes chères sœurs et leurs maris ainsi que leurs enfants pour leurs encouragements
Remerciements
C'est avec un grand plaisir que je résume en quelques lignes en guise de reconnaissance
et de gratitude envers les personnes qui m'ont aidé à la réalisation de ce mémoire car tout
dans ce travail et de leurs recommandations qui m'ont été très utiles. Je souhaite leur rendre
hommage pour leur amabilité, leurs disponibilités et leurs conseils dans la réalisation de ce
travail.
Je tiens à remercier aussi mes collègues pour leurs conseils et leurs contributions dans ce
travail.
pour leur accueil, cela m'a enrichi de conseils et de discussion me permettant de connaitre
d'autres cultures.
Enfin, je tiens à remercier les membres du jury qui ont pris la peine d'évaluer de travail.
VI
Dédicaces ........ .............. ...................... ............... ......... ...................... ..... ............ ..... .... .... . iv
Remerciements ...... .......... .................. .. ... ............... .... .................................... .... .......... ...... v
Table des matières .... ......... ....... ... ........................ ... .......... ................... ............................. vi
Liste des tableaux .............. .. ........................ .. ......... .... .......... ... ................. ................... ...... x
Liste des abréviations ... .................................. ........ ......... ............ ....... .......... ..... ..... ...... ... xiii
1.1 Mise en contexte ......... .............................. ............. ........... .. ... ........................... ... 1
1.2 Problématique ........ ... ... ....................... ..... ......... .. .... .... ...... .. .... .... ........................ . 3
1.3 Questions relatives au sujet de recherche ......... ..... ................ ..... ............. .. ...... ..... 5
1.4 Objectifs ........... .......... .......... ................... .... ........... ..... .... .... .... ........................ ... . 5
1.6 Organisation du mémoire ........................ .. ......... .... ........... ... .... ................... ..... .... 8
Chapitre 2 - Revue de la littérature .... ......................... .. .............................. .. ......... ..... ... .... . 9
2.1 Type de véhicule se lon la source d'énergie utilisée ....................... ..... ...... .. ...... .... 9
3.1 Architecture d'un véhicule électrique ...... ................................ ... .... ..... .. ....... ...... 28
3.2 Modèle dynamique longitudinal d'un véhicule électrique ........ ............... ............ 28
3.3.2 Estimation de l'énergie consommée ....... ..................... ..... ............... .... .... 33
4.1 Introduction ............................ .... ...... ............. ... ............. ..................... ... ... ......... 35
4.2.3 Estimation avec le filtre de Kalman (KF) ........ .. .. .... .. .................. ........ .... 41
consommée .. .... ...... .......... ............. ....... .... ........ .... ..... .... .... .......... .. ...... .. .. ....... .... 43
4.4 Estimation de l'autonomie ....................................... .. ........... ...... ... ... .. ... ... .. ....... 47
5.1 Validation et simulation ......... .... ... ........ ... ...................................... .................... 51
5.1.2 Expériences et tests ...... ........ .... .. ...... .. .. .. .... ............ .... .. .. .............. .. .. .. .. ... 53
5.2 Résultats des expériences et des simulations .... .. .. ........ .... .. .......................... .... .. 55
roulement .... ............. .... ............ .. ......... ..... ............. .... .. .......... ........ ... .... .. 55
consommée ..... .......... ....... ...... ...... .... ...... ......... .... ...................... .... .......... 60
Annexe A - Fiche technique du véhicule KIA SOUL EV 2017 .................... .... .... .. .. .. .. .. .. 77
IX
Annexe B - Article publié ....... ... ............ ........... .. ................................ ............. .......... ... ... 80
Tableau 2: Caractéristiques du KlA Soul EV 2017 ... .... ....... ........... .... ................... ..... ...... 52
Tableau 4: Récapitulation des tests .... .. .... .. .......... ....... .. .......... ... ........ ..... .................... ..... 54
Tableau 5: MSE et valeurs moyennes estimées .. .... ..... .......... .... .... ...... ... ... .... ... ...... ........... 58
xi
Figure 2: Architecture d'un véhicule à combustion interne .. .................... .............. ... ..... .. . 10
Figure 3: Architecture d'un véhicule électrique rechargeable à batterie ........... ... ... ..... ... .. Il
Figure 4: Architecture d'un véhicule hybride série ................. .......... ............. ...... ...... ..... .. 12
Figure 6: Architecture d'un véhicule hybride série/parallèle .............. .... .... ............. ........ . 13
Figure 7: Structure détaillée d'un véhicule électrique .......... ... ........ .... ....... .. ................... 28
Figure 9: Circuit électrique équivalent de la batterie ............. ...... ... .................. ......... ...... 32
Figure 10: Dispersion des pertes énergétiques d'un véhicule à combustion sur
des trajets: a) en ville, b) sur autoroute .... ............... .... ........ .. ........... .. ..... .. 36
Figure 12: Diagramme d'estimation des paramètres ................ ... ........... .. ................... .... . 38
Figure 13: Rendement du moteur électrique jusqu 'aux roues ................ ... ..... .. ........... .. .... 39
Figure 18: Coefficients estimés Vs valeurs réelles dans les conditions sans
neige: a) coefficient aérodynamique, b) coefficient de résistance
au roulement. .. ... ... ................ ........ ...... ... ......... ................... .... .... ........ ....... 56
Xli
Figure 19: Coefficients estimés Vs valeurs réelles dans les conditions avec
neige : a) coefficient aérodynamique, b) coefficient de résistance
au roulement. ........ .............. ..... .. ......... .... ... .. .... .. .... .. .......... ..... ....... ........... 57
Figure 24: Efficacité de conduite pour un test de conduite dans une zone
scolaire avec une présence de pente (Vmax = 30 km/h) .... .. .............. ... .... .. 64
Figure 26: a) module OBD Link MX, b) logiciel Torque Pro .......... .. ...... ............... ......... .. 85
XIII
AH: Ampère-Heure
AR : Autonomie Restante
CO 2 : Dioxyde de Carbone
cv: Cheval-Vapeur
kg: Kilogramme
kW: Kilowatt
kWh: Kilowattheure
km: kilomètre
xiv
Nm: Newton-Mètre
s : Seconde
Le secteur du transport représente le plus grand défi à relever dans la lutte aux
changements climatiques. Au Québec, près de 40% des émissions de gaz à effet de Serre
(GES) et 25% de la consommation totale d'énergie lui sont attribuables. Les véhicules
énergétique en transport. Le principal gaz à effet de serre émis par les véhicules est le dioxyde
transport (fig. 1). Avec 93 % de son énergie utilisée provenant de produits pétroliers, il
apparaît comme le moins diversifié de tous sur le plan énergétique. En effet, en raison de leur
densité énergétique très élevée (~ 12 kWh/kg), liée à leur état sous forme liquide en condition
carburants pétroliers ont été jusqu'à présent massivement utilisés par les principaux modes
s'affiche comme le deuxième contributeur aux émissions anthropiques de GES, après celui
de la production d'électricité.
Des efforts sont toutefois déployés pour réduire la dépendance du secteur du transport au
• Véhicules égers
Fret
• Maritime
• Aviation
• Industrie
• Habitation s
• Production électrique
• Autres
établi au deuxième rang avec une part de 27,2 %. Il a été suivi du carburéacteur (11,7 %), du
Aujourd'hui, plus d'un milliard de véhicules circulent dans les rues et les industriels
produisent plus de 80 millions de nouveaux véhicules chaque année. Compte tenu de cette
production en masse, les progrès du secteur du transport sont inévitables pour accommoder
les besoins modernes de déplacement avec certaines contraintes comme l'épuisement des
• Efficacité énergétique: le moteur électrique est 3 fois plus efficace que les véhicules
• Émissions de gaz à effet de serre: un véhicule électrique propulsé par batterie ou pile
à combustible ne produit aucun GES local par rapport aux véhicules à moteur à
combustion interne
1.2 Problématique
que cela dépend de plusieurs facteurs que ce soit au niveau du véhicule ou de son
environnement. Les paramètres qui affectent cette consommation peuvent être des inconnus
plusieurs paramètres, tels que le poids du véhicule, l'efficacité du moteur et des paramètres
et l'économie mondiale de carburant dans de tels systèmes est étroitement liée au contrôle du
flux de puissance entre chacune des sources d'énergies et le groupe motopropulseur [6] [7].
L'estimation en ligne des paramètres du véhicule a été un domaine de recherche actif au cours
des dernières décennies [8]. En raison de la différence de taille des véhicules, de la forme de
à l'air varient d'un véhicule à l'autre [9]. Cependant, ces facteurs sont normalement supposés
les mêmes dans la littérature quelle que soit la configuration du véhicule. Le coefficient de
route [10] [11] [12]. De plus, la masse du véhicule et le coefficient de traînée aérodynamique
est attaché à la voiture. Par conséquent, l'estimation des paramètres en ligne est nécessaire
4
car elle réduit l'effort d'étalonnage initial. De plus, il aide le contrôleur du groupe
De plus, pour les véhicules électriques, malgré tous ses avantages cités précédemment,
les véhicules électriques nécessitent une recherche plus profonde car ces véhicules présentent
Ainsi leur acceptation et leur adoption par les clients sur le marché est principalement
entravée par leur densité énergétique limitée, le long temps de recharge et le déploiement
électriques résultant de leur densité d'énergie limitée est considérablement réduite dans les
disponible pour chauffer les accessoires (en raison de la basse température) et surmonter les
enneigées et verglacées. Ces inconvénients des VE entraînent ce que l'on appelle « l'anxiété
d ' autonomie », qui doit être sérieusement combattue par les constructeurs automobiles afin
Par rapport aux problèmes cités précédemment, on peut alors se questionner sur les points
suivants:
consommation d'énergie?
Pour pouvoir être plus précis par rapport à la question précédente, on peut se demander:
des pertes par friction roues/routes et des pertes dues à la configuration de la route
configuration ?
1.4 Objectifs
d'énergie de la famille des véhicules électriques afin que l'on puisse estimer l'autonomie de
manière exacte.
l'autonomie restante permettant de réduire l'anxiété de distances chez les usagers des
véhicules électriques, notre travail de recherche consiste alors à inclure les valeurs des deux
• Établir un modèle de prédiction afin d'avoir une bonne précision sur l'estimation de
1.5 Méthodologie
consommation énergétique ainsi que l'autonomie restante, tout en se basant sur l'estimation
La première étape consiste à effectuer une revue de la littérature afin d ' établir une
Dans la stratégie d'estimation de ces deux paramètres, le modèle dynamique longitudinal est
l'autonomie du véhicule.
7
En premier lieu, afin de choisir les techniques d'estimation les plus adaptés à notre
paramètres est une étape importante pour ce choix. Cela a permis de relever les limites de
certaines méthodes des études antérieures et de développer, dans une deuxième étape, les
algorithmes utiles pour cette étude. Les algorithmes d'estimation développés s' appuient sur
deux approches différentes mais un peu similaires : l'approche des moindre carrés récursifs
(MCR) ou Recursive least square (RLS) et le filtre de Kalman. Une comparaison entre les
deux méthodes estimant les deux paramètres simultanément nous a permis d ' identifier
l'approche la plus appropriée à l' estimation en temps réel. Le filtre de Kalman s'est montré
plus efficace.
Une fois que ces deux paramètres sont estimés, ils seront intégrés au fur et à mesure dans
batterie pour obtenir l'énergie électrique sur un trajet ainsi que l' estimation de l'autonomie
restante. Pour cela, le modèle NARX ou Nonlinear AutoRegressive eXogenous est utilisé
pour ces estimations. Ces estimations sont basées sur deux critères: (i) la configuration du
réseau, et (ii) l'intégration des forces estimées dans le modèle. Ce scenario permettra de
savoir quelles configurations sont adaptés pour ce travail et aussi s'il est nécessaire d'inclure
La troisième étape consiste à valider les méthodes et modèles proposés tout en utilisant
les données réelles et les résultats expérimentaux obtenus sur le véhicule électrique KlA
SOUL EV de l'Institut de Recherche sur 1'Hydrogène (IRH) . Puis, finalement, les résultats
8
seront analysés et discutés tout en tenant compte des cycles de conduites et des itinéraires
définis.
La revue bibliographique est illustrée dans le chapitre 2. Cette partie illustrera les
aussi les facteurs influençant la consommation d'énergie d'un véhicule et aussi les différentes
présentée dans le chapitre 4 ainsi que les estimations de la puissance/énergie fournie par la
estimation.
Le chapitre 6 conclut le mémoire par une discussion sur les résultats et les perspectives.
Chapitre 2 - Revue de la littérature
innovation nécessite une continuité de travail pour l'efficacité énergétique et pour plus de
sécurité. Le but de ce chapitre est de présenter, d'abord, des généralités sur les différentes
technologies en termes de sources d'énergie existantes ainsi que les principaux organes
constitutifs. La deuxième partie sera consacrée aux facteurs qui peuvent influencer la
par de l'essence ou du gazole (diesel). Le moteur produit de l'énergie en brûlant à très haute
mode de fonctionnement de la plupart des véhicules de transport. Les gaz brûlés sont ensuite
Un moteur à combustion interne (MCl) est un type de moteur à combustion dans lequel
Différentielle
Les pistons du moteur descendent et aspirent le carburant et l'air par l'intermédiaire d'une
pompe. Les pistons remontent et le mélange (carburant + air) est comprimé dans les cylindres
du moteur produisant une combustion du mélange avec de l'étincelle. Cette explosion renvoie
les pistons en bas et ces derniers remonteront pour pousser les gaz d'échappement vers
un système de bielles et de vilebrequins, vers l'arbre de transmission qui est connecté avec
Un véhicule électrique rechargeable à batterie est un véhicule cent pour cent électrique
doté d'une seule chaine de traction électrique. Les véhicules électriques possèdent un
chargeur utilisé comme interface entre le pack de batterie et le réseau électrique. Ce pack de
Il
batterie alimente et transmet à l' onduleur le courant continu pour le transformer en courant
Différentielle
d'énergie différentes pour se déplacer. Généralement, on désigne par voiture hybride une
voiture qui recourt à un carburant et à l'électricité pour se mouvoir au moyen de deux moteurs
Ce type de véhicule possédant deux sources dont le moteur thermique n'est pas couplé
directement aux roues. La propulsion est assurée par le moteur électrique, lui qui est alimenté
par une batterie ou/et par un autre système (groupe électrogène) [17].
12
Différentielle
ROllL'
Différentielle
ROllL'
Le type d'hybride parallèle possède deux sources d'énergie pour se propulser: le moteur
thermique et le moteur électrique et qui sont associés aux roues. Les flux énergétiques
13
provenant de ces deux types de sources arrivent en parallèle aux roues pour assurer la traction
simultanée, comme son nom l' indique. Conceptuellement, le moteur thermique et le moteur
électrique sont connectés à la transmission par l'intermédiaire des voies séparées, pennettant
Différentielle
L'alternance entre les deux moteurs s'opère suivant plusieurs modes (un moteur principal
simultanée. Cela assure un usage sélectif du moteur thermique et électrique pour économiser
la consommation [19].
14
véhicule électriq ue
Pour pOUVOIT étudier comment l'énergie est utilisée, il faut savoir les facteurs qUI
influenceront sur la consommation de cette énergie afin de gérer l'énergie nécessaire pour un
déplacement. Dans cette partie, on va présenter des facteurs que l'on ne peut pas négliger
comme le mode de conduite du conducteur avec l'efficacité du moteur ainsi que les forces
résistives telles que la force aérodynamique, la résistance au roulement sans oublier la masse
et les conditions météorologiques qui peuvent faire varier la consommation de l'énergie d'un
véhicule électrique.
dont on joue avec les pédales de frein, accélérateur, etc, affectera sur l'énergie consommée.
de carburant, les émissions de gaz à effet de serre et les accidents [20]. La consommation
d'énergie et les émissions de C02 peuvent être réduites grâce à une variété de stratégies qui
temps, réduisant ainsi les émissions de C02 de la conduite d'un pourcentage équivalent [21].
fournissant des conseils statiques généraux au conducteur (par exemple, n'accélérez pas trop
rapidement, ne réduisez pas la vitesse, etc.) [22]. En allant plus loin, il est possible de fournir
des conseils en temps réel aux conducteurs en fonction de l'évolution du trafic et des
importantes [23]. Une consommation de carburant excessive peut être évitée si une stratégie
de conduite optimale est mise en œuvre en fonction de l'état environnant d'un véhicule et des
contraintes existantes [24]. Selon l'étude faite dans [25], voici des conseils, en deux étapes,
qui pourraient aider sur la façon de conduire afin d'optimiser la consommation d'énergie:
trajet le plus rapide. Sauf exception, il permet d'éviter les détours, empêche de se perdre (et
donc de parcourir des kilomètres pour rien) et fait gagner un temps précieux.
3. Vérifier la pression des pneus chaque mois: rouler avec des pneus sous-gonflés accroît
sensiblement la consommation d'énergie: 2,4% par 0,5 bar de déficit. Une pression trop
Les pneus se dégonflent naturellement, il est donc nécessaire de contrôler leur pression
au moins une fois par mois. Le contrôle doit s'effectuer à froid, c'est-à-dire après moins de
3 kilomètres de roulage, sinon il faut attendre 10 minutes pour laisser refroidir les pneus.
5. Rouler moins vite: rouler autant que possible à vitesse modérée. Rouler moins vite
permet également de ralentir plus facilement, en utilisant le frein moteur autant que possible.
Le frein moteur permet non seulement de réduire sa consommation et l'usure des freins mais
d'attente
de pouvoir maintenir la vitesse la plus stable possible, anticiper le trafic est essentiel.
C'est aussi un conseil essentiel pour la sécurité: on est moins pris au dépourvu si on
regarde le trafic au loin et qu'on anticipe les événements pouvant survenir (piéton qui va
traverser, voiture qui va dépasser ou freiner, feu de signalisation qui passe au rouge ... ).
Conduire comme si on n'avait pas de freins: maintenir le plus possible une vitesse stable,
8. Utiliser la climatisation à bon escient: la climatisation peut être très utile lors de
longues conduites en pleine chaleur mais n'est pas nécessaire la plupart du temps. On évite
les forces qui peuvent influencer la consommation. Cette sous-partie présentera les forces
roulant sur une surface causée par des effets non-élastiques. Selon [26], la résistance au
roulement des pneus est définie comme l'énergie consommée par unité de distance de
déplacement lorsqu'un pneu roule sous une charge. Elle est importante pour le mouvement
17
simulation car le mouvement d'un véhicule au sol est principalement déterminé par les forces
de frottement transmises par les routes via les pneus [28] . D ' après [29], la résistance au
roulement affecte non seulement la consommation d' énergie mais aussi le confort de
conduite pour les passagers ainsi que les couts d' exploitation du véhicule. Dans un contexte
écoénergétiques pourraient jouer un rôle important pour économiser de l'argent, limiter les
D'autre part, d'un point de vue général, l'aérodynamique est la science qui étudie les
phénomènes accompagnant tout mouvement entre un corps et l'air qui le baigne. Dans le
L'aérodynamique des véhicules joue un rôle beaucoup plus important dans les considérations
De même que les deux cités précédemment, la masse est un élément clé dans la gestion
de la consommation. La réduction de la masse est l'un des moyens les plus pratiques dans la
hybrides (HEV) à moteur à combustion interne et les VEH à pile à combustible, par rapport
aux véhicules conventionnels. Les résultats dans [34] ont confirmé que la consommation de
carburant augmente à mesure que la masse augmente et est différente pour différentes
donné en cas de collision. Une réduction de la masse des véhicules de 10% (en supposant
consommation de carburant pour les trains et les camions légers, moyens et lourds, 4% pour
[36] a fait une étude selon laquelle la température affecte l'autonomie d'un véhicule. La
voiture en question est une Nissan Leaf. La consommation d'énergie et l'autonomie varient
provenant de la batterie. Sur la base de ce test et des tests ultérieurs effectués ce jour-là (qui
avait une température moyenne de 14 0 F), l'autonomie estimée du véhicule était de 50 miles.
à diminuer, ce qui augmentait l'autonomie. La période d'essai la plus chaude était en moyenne
de 71 0 F et les tests ont abouti à une autonomie estimée à 91 miles. Tous les tests ont été
Durant le trajet, le véhicule nécessite plus d' énergie pour le confort des passagers comme
la climatisation, la ventilation, le chauffage, etc ... Or, des facteurs tels que la température
trempage thermique avec ces conforts ont été identifiés comme les principaux facteurs des
température et de l'humidité. Ainsi, on peut dire que ces conditions sont non-négligeables en
énergétique du fait que ces deux paramètres peuvent influencer la consommation malS
pourtant cette partie n'est pas vraiment prise en compte dans la littérature.
Plusieurs méthodes ont été rapportées dans la littérature et peuvent être classées comme
la moyenne.
• Recherche d'événements:
Cette catégori e regroupe les méthodes dont les conditions de conduite assurent une
L'étude [39] compare les coefficients de force aérodynamique obtenus par les
expériences respectives afin d'évaluer les effets des coefficients de fo rce aérodynamique sur
critique du vent est examinée en calculant la vitesse critique du vent pour des quantités
aérodynamique uti lisés dans cette étude ont été obtenus à partir d'expériences en soufflerie.
D'autre part, [40] montre une étude sur l'économie de carburant d' un véhicule lourd par
la réduction de la trainée aérodynamique. Les opérations routières réal istes sont prises en
compte en simulant des itinéraires de conduite typiques sur les distributions long-courriers et
urbaines; les variations de poids du véhicu le sont également prises en compte. Les résultats
montrent que les performances de ces dispositifs aérodynamiques dépendent à la fois de leurs
généralement deux fois plus de carburant que ceux conduits dans les zones urbaines. Les
avec des capteurs qui sont principalement disponibles dans les voitures particulières typiques.
mais étendus plus tard pour couvrir les cas de glissement combinés. E n ce qui concerne les
méthodes basées sur la dynamique latérale, plus robustes aux bruits haute fréquence, un
non linéaires robustes a été développé en utilisant la mise en œuvre numérique des conditions
longitudinale, la rigidité longitudinale des pneus est utilisée pour déduire le frottement de la
fonctionnent bien que dans des conditions d'excitation par glissement pur. Pour augmenter la
plage de travail, une logique d'intégration a été développée pour basculer entre les
Dans [42J, ils ont utilisé la méthode Coast down afin de comprendre les améliorations
pourraient être masqués par les essais traditionnels en soufflerie sur sol fixe. L'intégration de
• Calcul de la moyenne:
comme la méthode récursive des moindres carrés ou le fi Itre de Kalman ou des méthodes
d'estimation non linéaires. En effet, certaines études ont présenté une estimation simultanée,
En effet, l'étude dans [8J présente une estimation simultanée des coefficients
récursifs. Des données expérimentales d ' un véhicule lourd ont été utilisées pour tester
22
l'algorithme dans une simul ation. L'estimateur proposé a montré une précision adéq uate dans
la prédiction de la charge routière. De même que dans les études [43] et [44] qui estiment les
De son côté, [9] présente une estimation des deux coefficients: en premier lieu
simultanément, et puis séquentiellement avec deux méthodes en ligne basées sur le filtre de
Kalman : le modèle de base et l'étendu. Les expériences ont montré la difficulté de séparer
la résistance au roulement et la traînée de l'air. Il est montré que des estimations simultanées
des deux sont possibles, mais en pratique une trop grande variation de vitesse est nécessaire
pour obtenir des estimations précises avec les méthodes étudiées. Il est également montré
que lors de l'estimation d'un paramètre à la fois, des estimations précises peuvent être
générées. Cependant, il s'est avéré difficile de baser ces estimations les unes sur les autres,
Basé sur la dynamique longitudinale d ' un véhicule, l'étude [45] a proposé d'estimer les
paramètres qui peuvent varier durant un trajet par simulation y compris les coefficients
aérodynamique et résistance au roulement à l'aide d ' observateurs non linéaires basés sur le
les conclusions suivantes peuvent être tirées: Les limites des estimateurs basés sur la
L'étude dans [46] utilise un réseau de neurones RBF pour compenser l'influence négative
qu'une seule variable, ce qui permet que les centres et les étendues des neurones de la couche
23
cachée puissent être mis hors ligne d'une manière très simple. Les poids des couches de sortie
sont adaptés en ligne par une loi d'adaptation dérivée de l'analyse de stabilité de Lyapunov.
Une telle approche garantit une estimation rapide et robuste de la force de frottement pneu /
route, même dans les cas où des changements larges et brusques des paramètres de frottement
se produisent.
batterie ou l' autonomie restante fournie par le véhicule. Pourtant, cette estimation nécessite
une recherche plus profonde sur le fait que tout se base sur l'énergie restante ou l'énergie
consommée.
VE. L'une des caractéristiques les plus avancées d'un véhicule électrique, par rapport aux
système de freinage régénératif (RBS).L'étude dans [47] propose une approche d'estimation
de l'état de charge utilisant un système neuro-flou basé sur le regroupement soustractif est
avancé en comparaison avec le réseau neuronal à propagation arrière et les réseaux neuronaux
Elman. Les paramètres d'entrée pour modéliser l'approche d'estimation de l'état de charge en
ont montré que le modèle proposé présente une précision suffisante et surpasse les modèles
basés sur les réseaux neuronaux et les réseaux neuronaux Elman. Ainsi, le modèle proposé
sous différents cycles de pilotage montre une avancée remarquable dans l'estimation de l'état
de charge avec un potentiel élevé pour surmonter les inconvénients des méthodes
traditionnelles et fournit donc une approche alternative dans l'estimation de l'état de charge.
Selon [48], pour assurer une gestion fiable de la batterie dans les véhicules électriques
(VE), un estimateur avancé de l'état de charge (SOC) via une méthodologie d'apprentissage
means) basée sur un algorithme génétique est d'abord utilisée pour partitionner les données
d'entraînement échantillonnées dans le test basé sur le cycle de conduite d'une batterie
paramètres antécédents du modèle. L'algorithme des moindres carrés récursifs est ensuite
utilisé pour extraire ses paramètres conséquents. Pour assurer une bonne précision et une
vérifient que l'estimateur proposé présente une précision suffisante et surpasse ceux
consommation d'énergie des véhicules, y compris un cadre de boîte noire qui identifie les
collaborative basée sur la factorisation matricielle. En outre, une étude de cas de prédiction
de la distance à vide pour les véhicules électriques par des données de détection participative
25
est menée et évaluée de manière emplrtque, ce qUi montre que les approches peuvent
Avec l'étude dans [50], des paramètres influençant la consommation d'énergie ont été
utilisés afin de les implémenter dans le modèle dynamique longitudinal d'un véhicule pour
estimer la puissance fournie par la batterie et ainsi la puissance sera dérivée par rapport au
ressource pour les utilisateurs de VE. Cependant, l'ensemble de données dans cette recherche
est limité car il contient des données d'un seul conducteur, un véhicule, qui n'est pas l'un des
véhicules actuellement produits dans le commerce, et comprend une conduite sur autoroute
limitée. L'estimation de l'état de charge est l'un des facteurs les plus critiques pour résoudre
véhicule électrique.
Les travaux dans [51] sont basés sur le modèle VT-CPEM ou Virginia Tech
(vitesse instantanée du véhicule et niveau d'accélération) est utilisé comme variables d'entrée.
Ce modèle complète les simulateurs des véhicules disponibles dans la littérature ouverte en
produisant une erreur moyenne de 5,9% par rapport aux données empiriques. Les résultats
démontrent que dans les cycles de conduite urbaine, il est possible de récupérer plus d'énergie
26
conduite. L'étude a confirmé que l'avantage énergétique des véhicules électriques dans la
conduite urbaine pourrait avoir un impact significatif sur les choix d'itinéraire des gens et
Conclusion:
D'après les études que nous avons cité précédemment dans cette même section, nous
avons pu trouver les lacunes dans les autres études et ainsi pour délimiter les différents axes
cité dans la littérature. Ainsi, nous allons nous concentrer sur l' estimation du coefficient
D'un autre côté, nous allons nous focaliser aussI sur l'amélioration de
l' estimation/prédiction de la puissance fournie par la batterie à partir d'un modèle plus
efficace pour pouvoir donner une meilleure estimation de l'autonomie restant si on connait
est composée d'un réseau d'alimentation alternatif, d'un chargeur de batteries, d'une source
et de la transmission mécanique.
dynamique d'un véhicule et les paramètres physiques de ses composants seront présentés. La
utilisée pour développer les méthodes d' estimation des coefficients aérodynamique et de
résistance au roulement. De plus, un modèle physique de la batterie sera aussi décrit dans
cette partie.
28
Battene de Moteur
servItude TransmISSion électrique
( ) ( i ) (Tri a )
.....
\)
~.
VI
m Co nve rtlsse ur
m- OC/OC
m +
n
~
Chargeur
~. mba rq ué
.0 ( dress. )
C
m
Calculateur CCME
Gestion puissance
résultante exercée sur un objet est toujours égale au produit de la masse de cet objet par son
Chaque force appliquée sur un objet entraîne cet objet à accélérer dans la direction de la force
appliquée. Or, lorsque plusieurs forces sont appliquées sur un objet, il faut déterminer la force
résultante, soit la force équivalente à la somme vectorielle de toutes les forces agissant sur
cet objet [53]. La deuxième loi de Newton se résume par l'application de l'équation suivante:
l.Fext = Ma (1)
Dans le cas des véhicules, voici une figure montrant les forces externes appliquées sur le
longitudinale:
énergi e mécanique nécessaire pour faire tourner les roues. Elle peut être exprimée comme
suit:
30
(3)
Où Pb, 'Y'f et v sont respectivement la puissance fournie par la batterie, le rendement global
déformation instantanée des pneumatiques lors de leur contact avec la route, peut être
(5)
Dans l'équation (4), P est la densité de l'air qui dépend de la température, tandis que A et
aérodynamique. Le paramètre V vent dans l'équation 5 est la vitesse de l'air ambiant qui a une
valeur négative pour le vent arrière et une valeur positive pour le vent de face.
en particulier lorsque l'inclinaison 8 n'est pas négligeable. Elle peut être calculée à partir de
l'équation suivante:
31
La force de freinage (Ffreinage) est liée étroitement à la force d'accélération. Elle permet
chaleur. À partir du théorème de l'énergie cinétique, on peut obtenir cette force avec
l'équation suivante:
_ 1 2
Ffreinage - 2d mv (7)
implémenté en sens inverse, c'est-à-dire que la puissance de sortie de la batterie est calculée
en fonction du profil de vitesse v ainsi que des forces externes constitutives, en utilisant
l'équation 7:
modélisation exacte des batteries qui est difficile à concevoir et difficile à calculer au cours
des simulations.
32
<:1
Em , Rp _ J
r
1
(P,U.I' III
r ,.ch)
fournissent des informations pertinentes sur son fonctionnement, dans lesquelles se trouvent
(9)
(10)
RD = Roo (1 + Ao(1- SOC))
(11)
Ke : constante (V1°C)
Ao : constante (N/A),
La batterie est modélisée par une analogie du circuit équivalent. Ainsi, la tension aux bornes
d'énergie nécessaire pour parcourir les trajets peut être possible avec les informations des
itinéraires avec les destinations finales fixes et les durées des trajets avant les parcours en
l'énergie des batterie (équation (14)) nécessaire à chaque temps d'échantillonnage (t) est
Pm : la puissance mécanique,
Donc, l'énergie jusqu'à la destination finale peut être estimée comme suit :
4.1 Introduction
L'estimation des paramètres dynamiques d'un véhicule est une solution afin de savoir
sur quelle partie du véhicule consomme le plus d'énergie afin d'optimiser sa consommation.
véhicule: la masse ainsi que les conditions de roulages telles que les coefficients
des paramètres les plus influents dans la consommation d'énergie. Estimer ces paramètres
[55J et [56J, l'aérodynamique et la résistance au roulement présentent une grande partie non-
négligeable en termes de pourcentage des pertes sur l'efficacité énergétique selon les figures
ci-après: en premier lieu pour les véhicules à combustion et en deuxième les véhicules
Acceuoire:s
A'lOdy.......qu<
1/ 131.
3".1.
Corb\nm
lOO1~
- 2~~1 TnDJaUSS>OIl
1
fremage
2%
Pene _ _ 1 Pme TfUlSIIUSSlOD
6~. S~.
Figure JO: Dispersion des pertes énergétiques d'un véhicule à combustion sur des
trajets: a) en ville, b) sur autoroute
66
1~nergie
,... Efflcacjt~: 66/100 =66" et 6 /73 91" • «utile
N
p, .r , 111 Il 7
~------------------------------------------~ - . -
Perto, do 10 boltollt Ptrto, du systtrM do
10
Pertes des KCfSSOlffS et
• Uluiairts .ltctll'iUes
.... _._._._ .
13 propul«on ""'trique
1Én~~je
Efflcadté: 74/100 =74 "et 4 93 - 80" • 1 utile
l ' ql 8
1
roulement
D'après la revue de littérature citée dans le chapitre 2 sur l'estimation des coefficients
aérodynamique et de résistance au roulement, la majorité des travaux effectués ont utilisé les
méthodes des moindres carrés récursifs et le filtre de Kalman tout en se basant sur l'équation
puissance mécanique (Pm) pour la traction du véhicule. Dans ce cas, l'énergie des batteries
ACTIONNEUR
:=:====-:::::
PACKS DE
] . ..... ' COVERTISSEUR
BATTERIE
,
,~ . _._ . _ . _ . _ . _ . ! . _ . _ . _ . - _._._ ._.-._.- _._._.- _.,.
Données du véhicule : Estimation :
!, • Puissance de la batterie • Coefficient aérodynamique
• Vitesse • Coefficient de résistance au
roul"""",t
Estimateurs récunifs:
• Model dynamique longitudinal
• Moindres carrés récunifs
.. . Fi1trede ~ 1_
l . _ . _ . _ . __ . __ . _. _ . _ . _ . _ . _ . _ . _ . _ . _ _ _ . _ . _ . _ . _ . _ '
coefficients
1 (18)
[Mav + Mgsin(8)v + 2 d mv 2 ]
Puisque la puissance est appliquée sur l'axe du véhicule au niveau des roues, il est
nécessaire de connaître les pertes de puissance entre le moteur et les roues. Les pertes
d'énergie au niveau du moteur et des roues sont des pertes électriques générées par le moteur
39
Pack de
batterie
TRANSMISSION
'7 re d
suivant:
Donc, il est nécessaire d'estimer le rendement final (T]) du fait que le rendement varie
Les algorithmes d'estimation récursive sont plus utilisés pour les applications en ligne,
embarqués. Plusieurs algorithmes couramment utilisés sont le RLS, le RLS avec facteur
Dans cette partie, pour estimer les coefficients aérodynamique et résistance au roulement,
nous avons opté pour deux méthodes: le RLS et le KF. Ces deux méthodes sont étroitement
liées.
40
4.2.2 Estimation par la méthode des moindres carrés récursifs (MeR ou RLS)
Le RLS est un filtre adaptatif qui estime de manière récursive certains paramètres pour
minimiser une fonction de coût des moindres carrés linéaires pondérés. Cette fonction de
coût est associée aux signaux d'entrée. Le RLS donne des performances de qualité lorsqu'il
fonctionne dans des conditions variant dans le temps [57]. On peut dire que le RLS est un
filtre adaptatif déterministe car il considère que le signal va avoir la même réponse à la même
impulsion d'entrée. Cette méthode se rapproche d'une estimation des moindres carrés avec
une fonction de coût. Cette fonction de coût inclut un facteur d'oubli À qui donne plus
Cela permet de suivre le signal souhaité plus rapidement, en revanche, il est plus
complexe par rapport aux autres comme le LMS, NLMS. Le facteur d'oubli est lié à la vitesse
de convergence; si c'est un, cela signifie que le système a une mémoire infinie.
Dans cette approche, le vecteur de paramètre inconnu est mis à jour de manière itérative
à chaque pas de temps en utilisant les données mesurées d'entrée et de sortie précédentes.
y = c/J8
y = [Mav + Mgsin(8)v]
- Mgcos(8)v -'21 PA (v + Vwind) 2 ] (20)
o= [1] J1 cdF
ALGORITHME:
41
g(n) (22)
k (n) - ------=--:-'-:---
- À+cpT(n)*g(n)
(23)
e(n) = den) - eT (n - 1) * cp(n)
Nous avons divisé en deux étapes le calcul du vecteur de gain variable dans le temps k(n)
pour réduire la complexité du calcul. La première partie utilise la valeur passée de l'équation
de Riccati P(n) multipliée par le vecteur d'entrée. Pour calculer le vecteur de poids w(n) , nous
utilisons e(n) la différence entre le signal souhaité d(n) et le signal estimé, multiplié par le
La valeur de À est reliée au rapport signal bruit (SNR). 8 a une valeur faible pour un SNR
élevé et une valeur élevée pour un SNR faible. La gamme de À est 0 < À < 1.
Le filtre de Kalman peut être utilisé si la dynamique des paramètres est connue.
Il est perçu comme un estimateur optimal, qui fonctionne selon un processus en deux
étapes pour estimer les paramètres d'intérêt. Ce filtre estime d'abord les variables d'état
courant puis les met à jour à la réception de la mesure suivante [58]. La structure du KF
y = [Mav + Mgsin(8)v]
1
[Pb - Mgcos(8)v - Z-vpA(v + V w ind)2]
(26)
X = [1] Il CdF
précédemment estimées, car cela est plus efficace en termes de calcul. Le filtre de Kalman
va estimer les variables d'état avec les mesures pour minimiser l'erreur quadratique. Ce filtre
convient au traitement numérique en temps réel lorsque les données d'entrée présentent du
Prédiction :
Correction :
consommée
Dans cette partie, les réseaux de neurones artificiels ou Artificial NeuralNetwork (ANN)
ont été utilisés pour l'estimation de la puissance/énergie fournie par la batterie. D'après [59],
les réseaux de neurones ont été suggérés comme outils pour la solution de problèmes
d'optimisation combinatoire difficiles. De plus, [60] affirme que l'avantage des modèles ANN
est qu'ils n'ont pas besoin de plus de paramètres ou ne compliquent pas les calculs
Les ANN sont des outils informatiques qui imitent le flux de traitement de l'information
des réseaux de neurones biologiques. Par exemple, les ANN sont constituées de plusieurs
l'espace de sortie, où la relation entre l'entrée et la sortie est suffisamment complexe pour être
poids qui est modifié de manière itérative par l'algorithme ANN afin de minimiser la
différence entre la sortie connue (cible) et la prédiction. Les architectures ANN peuvent être
réseau de neurones récurrents (RNN). Dans le FFNN, les informations sont propagées à
travers le réseau dans une direction de l'entrée vers la sortie, contrairement au RNN où les
sorties précédentes peuvent être réutilisées comme entrées sur le réseau. Le choix d'un réseau
Dans ce travail , nous avons sélectionné l'auto-régression non linéaire avec entrées
véhicule.
Le réseau NARX est une architecture neuronale dynamique couramment utilisée pour la
dans la littérature ont démontré la supériorité du réseau NARX sur les réseaux de neurones
formes en ligne et la conception de contrôleurs. Selon [61], le réseau NARX est une
systèmes dynamiques non linéaires. Le réseau NARX peut être facilement et efficacement
basés sur les réseaux neuronaux, tels que les architectures TDNN et Elman.
NARX appartient à la classe des réseaux récurrents, où les entrées sont constituées des
valeurs passées de la sortie ainsi que de certaines entrées exogènes. De manière analogue au
réseau à anticipation traditionnel, les entrées exogènes peuvent être considérées comme les
caractéristiques susceptibles d'influencer les sorties. L'architecture NARX est si flexible que
les entrées peuvent être soit les valeurs vraies passées de la cible (NARX en boucle ouverte),
soit les valeurs prédites passées des sorties (NARX en boucle fermée). Dans l'étape
d'apprentissage, nous avons utilisé le NARX en boucle ouverte pour accélérer le calcul,
tandis que dans l'étape de test, le NARX en boucle fermée est utilisé afin de faire des
prédictions en temps réel, c'est-à-dire lorsque les données arrivent. Les NARX en boucle
ouverte et en boucle fermée sont présentés dans la figure 14, où u(t) est l'entrée exogène.
45
u(/) u(/)
Focd Focd
Fo rw ard Fo rw ard
Y</) K c t\vo rk
K c two rlc
y(/)
suivante:
Où
Une autre flexibilité du réseau NARX réside dans le fait que la fonction f peut prendre
plusieurs formes. Dans ce travail, f est représenté par la fonction de réseau d'anticipation. Par
exemple, un réseau NARX implémenté avec un réseau FeedForward avec une couche cachée
Comme le montre la figure 15, le réseau NARX implémenté dans ce travail peut être
représenté par les caractéristiques typiques d'un réseau FeedForward (couches cachées,
neurones, etc.) et le délai. Les principales caractéristiques du réseau NARX utilisé dans ce
travail ainsi que leurs valeurs typiques sont présentées dans le tableau 1. Il est à noter que
nous ne présentons dans le tableau 4 que les paramètres typiques utilisés dans ce travail. En
fait, dans nos simulations préliminaires, nous avons testé plusieurs autres configurations afin
1 10
Figure ]5: Diagramme schématique d'un réseau NARX implémenté avec un réseau
à réaction contenant une couche cachée avec] 0 neurones
Le réseau NARX tel que conçu est formé à l'aide de données de séries chronologiques
avec une sortie connue (cible) ainsi que certaines fonctionnalités (exogènes). Dans l'étape de
formation, la fonction de coût utilisée dans ce travail est l'erreur quadratique moyenne.
convergence rapide.
Le réseau NARX est utilisé dans ce travail pour modéliser le profil de pUissance
électrique ainsi que l'autonomie restante. En général, le choix de l'entrée réseau dépend de
leur influence connue sur les données de sortie. Par exemple, dans le cas de la modélisation
Pour l'estimation de l'autonomie restante, les entrées sont la vitesse du véhicule, l'énergie
Les données utilisées dans ce travail sont enregistrées lors de plusieurs tests de conduite
qui ont été réalisés dans de nombreuses conditions afin d'analyser leur effet sur la
divisé en deux parties principales: la formation et les ensembles de données de test. Pour
mettre en évidence la capacité de généralisation du réseau NARX, nous avons utilisé jusqu'à
deux données de tests de conduite comme ensemble d'entraînement, dont 75 % sont utilisés
pour la formation tandis que le pourcentage restant (25%) est utilisé pour la validation. En
tant que tel , le but de l'ensemble de données de validation est d'ajuster le réseau pour l'étape
de prédiction. D'autre part, l'ensemble de données de test se compose des données de trajets
restantes (10 tests de conduite). En termes simples, notre objectif est de montrer que le réseau
Toutes les données utilisées dans ce travail ont été enregistrées à partir du bus CAN du
véhicule.
Généralement, l' autonomie restante (AR) d ' un véhicule peut être estimée en se basant
un examen plus approfondi de cette équation ne montre qu'aucun des termes impliqués n'est
batterie Eres est à peine estimée à partir du modèle de batterie complexe, tandis que
l'effIcacité de conduite future éft est généralement supposée être égale à l'efficacité passée
la plus proche. L'effIcacité de conduite peut être décrite par la relation générale suivante:
__ E(lld) (36)
ett = lld
distance historique passée comme illustré sur la figure 16 montrant le profil de vitesse d'un
100
-..
]
'-"
80
- Q)
:0:s
I1d
:.a
- Q)
60 :04 .:
;>
::s
"'0
Q) 40
Vl
...
Vl
Q)
......
> 20
exemple, lorsque /)"d est égal à la distance de conduite cumulée (di)' l'efficacité de conduite
courte durée, c'est-à-dire /)"d < di , l'efficacité de conduite est appelée efficacité historique.
Intuitivement, on peut affirmer que lorsque le /)"d diminue, l'efficacité de conduite sera
plus instable en raison de la capacité de récupération d'énergie de ,'EY. Ce fait induit un défi
important dans le choix de la valeur appropriée de /)"d afin d'éviter une sous-estimation ou
une surestimation de la plage restante qui, à son tour, peut augmenter l'anxiété de la plage.
50
Pour résoudre ce problème, nous proposons l'utilisation du réseau NARX pour apprendre
directement l'autonomie réelle restante. Dans ce cas, nous utilisons la vitesse du véhicule, la
l'autonomie réelle restante est utilisée comme cible. La véritable autonomie restante est
complète de la batterie. Le réseau formé est ensuite utilisé pour fournir l'estimation en ligne
N
(37)
MSE = ~ l
1
(x t - Xt)2
Où N est le nombre de points de données, xl' état estimé et x l'état de référence. Cette
méthode ne peut donc être utilisée que lorsque la vraie valeur de l'état est connue.
Chapitre 5 - Expérimentations et simulations
Dans ce chapitre, les expérimentations et les résultats seront décrits afin de valider les
L'Institut de Recherche sur l' Hydrogène (IRH) de l' Université du Québec à Trois-
Rivières (UQTR) dans lequel est réalisé ce travail dispose d'un véhicule électrique: KIA
Soul EV 2017 avec une vitesse maximum de 145 kmlh. Tous les tests de conduite présentés
dans ce travail ont été réalisés en utilisant ce véhicule. La figure 17 montre une image du
véhicule d'essai. Dans le tableau 2, nous présentons les paramètres du véhicule tels que
Figure 17: KlA Soul EV 2017 : A) à l 'intérieur de l 'IRH, B) sur asphalte, C) sur
neige
52
Phlyslque
.
Roues 205/60R1692H
Masse 1492 kg
Performances
Type de batterie Polymère de lithium-ion
Capacité de la batterie 45 Ah
Paramètres
A 2.88 m 2
Cd 0.35
P 1.22 kg m- 3
Il 0.018
Nous avons effectué plusieurs tests avec le véhicule KlA Soul EV 2017 pour récupérer
des données afin de valider les algorithmes d'estimations proposés. Les données
expérimentales incluant la température ambiante et les positions du véhicule ont été acquises
logiciel dédié: Torque Pro, avec une fréquence f= 1 Hz. (Voir Annexe C).
Des tests coast-down ont été faits afin de trouver les valeurs du coefficient de résistance
au roulement. Ce test est un procédé fréquemment utilisé et nécessaire pour évaluer les
valeurs expérimentales des forces résistantes et par conséquent les valeurs expérimentales du
neige avec un sol asphalte pour notre cas). Le principe de base de cet essai est le suivant:
accélérer le véhicule à une certaine vitesse, puis le laisser décélérer en vitesse point mort
résistance totale à l'avancement est alors calculée à partir de la deuxième loi de Newton
o
Sans neige -20 e 34.7239 11.2072
13 Vi lle
Pour valider les estimations, on a fait des tests dans différentes places avec différentes
Dans cette partie, on illustre les résultats sur l'estimation des coefficients aérodynamique
proposées : les moindres carrés récursifs (RLS) et le filtre de Kalman (KF). Pour chaque
méthode, on a utilisé les valeurs du coefficient de résistance au roulement obtenu avec les
0.38 1 1
- Réfén.::nc.:e
i1
0.175 •
RI.S
- --- Kalmllil
i
0.37 !
.
gIUM +
ê2 0.J6 ~!
>. ! 1
1l . ;
~ 0 JS5! i..~, .......
~ U.35 i "......-.-.---------.-.-.-. ----_ .... _---
'"'ü i
8 0 .\45 +
j
0.34 ~
:
OlB
.j
j
(a)
OJJ j
o IIKI :no mu -'{}fI 500 NXI
Temps (sec)
56
0.0194 ,
1
i
0.0193 !
i
"!1 O.O)l):! ,~
6 !
~,,0.019 1 i,
~
,!
!"
0019 ,\
E 0.OIX9 ' ,
(1 [1
\J>~:~.::-.~:.::: _~~~=._~~~~~~._~~~=._~~~._. . .-.:'.
r'j-..
,
' ••.•• _._ ••••.•••••• _._ •••••••••• J'"
~ 'Ir
g 0.01 8 i\i
U i.'
!~
O.01R7 +i (b)
1
O.OIU i
o 1110 200 300 400 5(10 600
Temps (sec)
Figure 18: Coefficients estimés Vs valeurs réelles dans les conditions sans neige:
a) coefficient aérodynamique, b) coefficient de résistance au
roulement.
La figure 18 représente les résultats des estimations des coefficients obtenus à partir d'un
test effectué dans des conditions sans neige. Comme on peut le voir sur cette figure, les
coefficients estimés par différentes approches convergent vers les valeurs réelles du
roulement pour une route asphaltée sec (Il = 0,019). La valeur de Il est supposée être la valeur
Ces résultats montrent que les coefficients estimés par toutes les méthodes sont presque dans
la même gamme. Cependant, les valeurs estimées par KF sont plus précises que RLS.
Pareil qu'avec les tests sans neige, les coefficients aérodynamique et résistance au
roulement convergent vers les valeurs obtenues à partir du coast down dans les conditions
... ", -- - -
i
' - - - - -- ---
- Réfêrcncc
i
i -·'· KlS
i
i --- Filrrc de Kalman
04 ......
!!
1
!..l i
, i
r \\.,
-2'
'~ o~rl-··-··~··:::·':::··:.:;·':::·,:.:;·':::·'''''·-''''·,=·,~··-.-----------_-.,-.--.,-.,-.--.--.--.,-.--.--.--.--.,-._-._-
._
-~t .,.,.-.,._.,.,. .,.,. ._.,.- ,.,
1 ----------------_._- ... ----- ..... ---
~- u1,
"1 1
;
(a)
025 '
o '<XI ...
-Rérérence
RLS
- - ' Fihr< de Ka1man
_..._-----
0.0 2~ 2
(b)
0.024
o 100 200 300 400 soo 600
Temps (scc)
Figure 19: Coefficients estimés Vs valeurs réelles dans les conditions avec neige :
a) coefficient aérodynamique, b) coefficient de résistance au
roulement.
La convergence significative des valeurs estimées vers les valeurs réelles illustrées dans
En comparant les résultats des simulations, selon les figures, on remarque que les valeurs
estimées à partir du filtre de Kalman sont plus proches des valeurs réelles que celles du RLS .
58
Trajets Type de MSE-RLS MSE-KF Valeurs- Valeurs- MSE-RLS MSE-KF Valeurs- Valeurs-
route RLS KF RLS KF
1 Sans neige 0.00 Il 2.4205 e-4 0.3626 0.3516 7.1797 e-7 6.4158 e-8 0.018 0.019
2 Sans neige 0.0093 e-3 0.0013 0.3531 0.3550 1.2943 e-6 6.4671 e-7 0.0189 0.019
3 Sans neige 0.0012 1.9941 e-4 0.3633 0.3516 5.9337 e-7 5.7551 e-8 0.0188 0.019
4 Avec neige 3.1002 e-4 5.7619 e-4 0.3407 0.3538 1.1224 e-6 1.1140 e-6 0.0250 0.0250
5 Avec neige 2.5665 e-4 6.6577 e-4 0.3441 0.3537 1.1375 e-6 1.1339 e-4 0.0250 0.0250
6 Avec neige 1.6265 e-4 8.3625 e-4 0.3499 0.3545 7.9022 e-7 7.9412 e-7 0.0250 0.0250
7 Avec neige 2.3276 e-4 5.6984 e-4 0.3453 0.3535 1.0749 e-6 1.0739 e-6 0.0250 0.0250
8 Avec neige 1.4136e-4 9.5087 e-4 0.3493 0.3556 7.0072 e-7 7.0624 e-7 0.0250 0.0250
9 Avec neige 2.1427 e-4 7.9605 e-4 0.3462 0.3545 1.0305 e-6 1.0311 e-6 0.0250 0.0250
10 Avec neige 2.3541 e-4 7.4813 e-4 0.3442 0.3546 1.0430 e-6 1.0435 e-6 0.0250 0.0250
11 Sans neige 0.00\4 5.9646 e-4 0.3517 0.3523 1.2418 e-6 1.8089 e-8 0.0189 0.0190
12 Sans neige 2.2694 e-5 4.9607 e-4 0.3515 0.3524 3.3519 e-8 7.2025 e-8 0.0189 0.0190
13 Sans neige 0.0045 0.0036 0.3458 0.3612 1.0\94e-5 3.1568e-8 0.0\84 0.0\90
14 Sans neige 0.0110 0.0064 0.3782 0.3737 2.6892 e-7 2.529\ e-7 0.0189 0.0180
On remarque qu'il y a plus de fluctuations des valeurs estimées dans le cas de la conduite
réelle probablement en raison des erreurs dans les signaux mesurés ou de petits changements
soudains des conditions routières . De plus, ces résultats indiquent que l'estimation du
coefficient de résistance au roulement est plus précise et plus stable que le coefficient de
59
traînée aérodynamique. Cela peut être dû à la variation de facteurs inconnus tels que le vent
pendant le voyage.
Sur le fait que les résultats obtenus à partir du filtre de Kalman sont plus précis, on
À partir des paramètres estimés, nous avons effectué une analyse de la contribution de
chaque force pour créer une carte énergétique des forces externes appliquées au véhicule. La
figure 20 montre les résultats des tests sur route enneigée et sur route sèche effectués sur le
même itinéraire: en ville .. Les résultats indiquent que la contribution de la force de résistance
au roulement est plus élevée pour la route enneigée que pour les essais sur route sèche. Ce
résultat est tout à fait attendu puisque la présence de neige augmente la résistance au
plus élevée dans le test de conduite sur route que dans l'état de la route en ville,
vitesse, la contribution de la résistance au roulement ne varie pas beaucoup pour le test sur
route sèche
60
1.2
Energie total - Energie totale
~ lO Aérodynamique Aérodynamique
~ Résistance au roulement - Résistance au roulement
u 0.8
.~
0 0.6
l'l0
u G.4
u
j 0.2
0.0
0 100 200 Temps (5) 400 500 600 o 100 200 Temps (5) 500
- EDcrgietotale ~
- Aérodynamique (d)
- Résistance au roulement
Figure 20: Carte énergétique des différents tests de conduite : a) route enneigée en
ville, b) route sèche en ville, c) autoroute sèche (vitesse :s 100 km/h),
d) campagne sèche (vitesse :s 70 km/h)
Dans cette partie, nous allons illustrer les résultats obtenus sur l'estimation de la
puissance fournie par la batterie ainsi que l'énergie consommée. En fait, la puissance est un
paramètre énergétique fondamental des véhicules électriques à partir duquel la plupart des
Nous avons utilisé les valeurs des coefficients obtenues à partir du filtre de Kalman pour
calculer les forces. Dans cette étape, on peut montrer qu'en introduisant les forces, on peut
40
--- &rimée .\.« forces
J5
JO
...
10
-s l 1
(a)
o 100 200 lOO 400 SOO 600
T<-mps (,,:c)
40
lO
20
-10
-20
(b
_lO L - - - - . . l - - - = - - - - - - - - - - L - - - - - - - - - - - - - - -
o 100 200 lOO t\OO
Tcml'" (o;cc)
- - E..... ink.~l\cck.~I~~
1.6
1.4 ~
:=
~
~1.2 -
.g
:;
" l-
B
·î.08
t .
:§
0.6-
0.4 -
0.2
(a)
o
o 100 2011 300 400 soo
Temps (sec)
0.7
~0.6 1
~
.g Il.S
:;
E
~
~ 0.4
':i)
t
jO.3
Il.2
0.1
(b
1 1
100 200 300 4011 SOO
Temps (sec)
Ces figures ci-dessus nous montrent la capacité du NARX à estimer le profil de puissance
instantanée de la batterie. Ces estimations ont été faites en temps réel : aucune donnée du test
n' a été utilisée pendant la phase d'entrainement et les prédictions sont effectuées au fur et à
Dans cette section, nous présentons les analyses sur l'estimation de l'autonomie. Cela
passe par l'évaluation de l'efficacité de conduite. C'est une étape critique dans l'estimation
de l'autonomie restante. Dans cette partie, on remplacera le temps par la distance sur les axes
La figure ci-dessous (figure 23) montre l'efficacité de conduite pour un test de conduite
sur autoroute. On constate que l'efficacité de conduite cumulée est presque constante (5
km/kWh) pour des longs trajets parce que les conditions de conduite comme la vitesse ne
D'autre part, les fluctuations des efficacités d'enregistrement de courte durée s' atténuent
lorsque la distance d'enregistrement augmente. Cette figure est très pertinente et illustre toute
30
•
Efficacité cumul~
~ à jour tous les 1 km
Mise à jour tous les 5 km !"
25 ""
,,,l'
,
,,,
,,,
,,,
.
'I
,
,,
.
1 1
l ,
l ,
L' 1
'\
\ , ;
1'"
'- ... l '
1 _,
5 ... , ...
\
\~ .L --- \ ;'
"
"
5 10 15 20 25 30 35
Distance parcourue (km)
Figure 23: Calcul de l'efficacité de conduite pour un test de conduite sur autoroute
(Vmax = 100 km/h)
64
Les différents rendements de conduite présentés dans cette même figure sont estimés à
partir des informations énergétiques et ne diffèrent que par la distance de mise à jour.
Cependant, une mise à jour tous les 1 km est susceptible d ' augmenter l' anxiété par
rapport à une mise à jour tous les 5 km qui est plus stable. Le fait est que l' estimation de
l' autonomie n ' est pas aussi simple qu' elle semble. De plus, il faut faire attention au choix de
la distance de mise à jour afin de fournir une bonne estimation plus précise de la distance
restante.
Afin de prendre en compte chaque force sur la consommation énergétique, la pente d ' une
route peut affecter l' efficacité de la conduite. La figure 24 illustre cet effet dans les conditions
d ' une route vallonnée dont la pente positive correspond à une descente et une pente négative
8 _ Profile de la pente
Efficacité cumulée
• Mise à jour tous les SOO m
7
0 .04
0 .0 2
'ô'
~
-
\00
'-'
~
• 0 .00
d
~
Il.
-0.0 2
l
Figure 24: Efficacité de conduite pour un test de conduite dans une zone scolaire
avec une présence de pente (Vmax = 30 km/h)
65
pour les 500 premiers mètres et diminue pour les 500 suivants. En effet, pendant les 500
premiers mètres parcourus, le véhicule roule en descente ainsi il consomme moins d'énergie
tandis que quand le véhicule roule sur les 500 derniers, il consomme plus d'énergie pour la
Il est assez compliqué d'estimer l'autonomie restante sans avoir recours à l'efficacité de
conduite mais il est nécessaire de bien choisir les mises àjour. Pour contourner ce problème,
on a utilisé le réseau NARX pour connaitre l'autonomie restante et ce réseau entrainé est
1
Estimée \'Îa NARX 1
1
250 Historique basé sur les données cumulées 1
1
Historique basé sur les mises à jour tous les 1 km 1
1
• Historique basé sur les mises à jour tous les 5 km 1
1
â 200 1
.:oc 1
1
'-"
1
~ 1
1
S14 150 1
CI) 1
'"'
CI) 1
1
·ë ; 1
0
c:
100 ,~ ,
l'
\
1 t
, 1
...
0 \ 1
•
~ ,
... ~
- - - - - __ J _ I
Il
~ 50
'... ,
\,
1
0
5 10 15 20 25 30 35
Distance parcourue (km)
Selon la figure sur l'estimation de l'autonomie sur autoroute (figure 24), on a utilisé deux
calculs de mise àjour avec le réseau NARX pour l'estimation de l'autonomie. On observe
que la portée restante estimée avec le NARX est plus précise que la méthode de calcul sur
l 'historique traditionnelle de mise à jour. Pour les méthodes basées sur la mise à jour tous les
5 km et les données cumulées, fournissent des estimations qui suivent assez bien la tendance
réelle de l' autonomie restante pour une distance parcourue supérieure à 15 km. En effet,
l'efficacité de conduite est presque constante avec ces distances et ainsi l'autonomie peut être
facilement calculée à partir des historiques des données le plus proche. Il convient de noter
que la méthode basée sur l ' historique des données de mise à jour tous les 1 km augmente
l' anxiété de l'autonomie. En définitive, le réseau NARX est très robuste dans l'estimation de
l'autonomie restante et constitue une alternative sérieuse aux approches des constructeurs
pour estimer l'autonomie restante d'un véhicule électrique. De plus, le NARX n'a pas besoin
d'aucune hypothèse sur la distance de mise à jour, ce qui le rend très utile pour l' estimation
Dans ce travail, on a élaboré des algorithmes pour les simulations afin que nous puissions
avoir des résultats satisfaisants. Dans ce cas, il est nécessaire de discuter de ces résultats
obtenus afin d' améliorer un peu plus le travail et de savoir si les objectifs initiaux sont atteints
En effet, l'idée est de savoir comment le véhicule utilise l' énergie qu'il a consommé
permettant ainsi de savoir jusqu'où le véhicule peut rouler avec le reste d'énergie. Les
chapitres l et II nous permettent d'avoir des idées et de présenter le sujet pour mettre en
évidence le défi d'estimer précisément l' énergie consommée d'un véhicule électrique afin de
conducteur. Pour cela, notre objectif principal était de proposer une approche sur l'estimation
de la consommation d'énergie de la famille des véhicules électriques afin que l' on puisse
consommation énergétique, ce projet consistait alors à inclure les valeurs des deux
résistance au roulement dans l'approche de gestion énergétique. Donc, cela nous a permis de
• Établir un modèle de prédiction afin d'avoir une bonne précision sur l'estimation de
Le modèle longitudinal du véhicule utilisé, décrit dans le chapitre III a été simplifié afin
de faire une étude sur l'impact des forces aérodynamique et résistance au roulement sur la
estimation non efficace de ces deux paramètres est intégrée dans l'estimation de l'autonomie,
une valeur résiduelle non fiable est indiquée augmentant l'angoisse de l'autonomie et
accentuer la peur du conducteur face au panne d'énergie. Cependant, en offrant une bonne
estimation, précise et fiable de ces paramètres, il est plus facile pour le conducteur d'avoir
Dans notre cas, on a utilisé deux algorithmes pour l'estimation des coefficients
algorithmes a permis de déterminer ces coefficients. En se basant des résultats obtenus par
rapport à la convergence, le filtre de Kalman a été plus performant par rapport aux moindres
carrés récursifs. L'estimation de ces paramètres nous a permis d'effectuer des analyses de la
D'autre part, un modèle NARX est utilisé pour l'estimation de la puissance ainsi que
l' autonomie restante du véhicule. On a pu remarquer qu'en rajoutant les forces dans le
modèle, on obtient plus de précision par rapport sans les forces . Cependant, on a fait quelques
analyses et calculs sur l'estimation de l'autonomie restante avec des mises à jour de données
permettant de prédire cette autonomie via le modèle NARX. Cela nous a permis d'avoir une
nécessitant aucune hypothèse et nécessite que peu de données de pour entrainer ce qui le rend
Les algorithmes élaborés dans ce projet ont été simulés en utilisant des données réelles
obtenues sur le véhicule KIA SaUL EV de l' Institut de Recherche sur 1'Hydrogène (IRH)
connecté de son OBD 2 avec un aBD Link et le logiciel Torque pour l'acquisition des
travail.
Perspectives:
Les algorithmes élaborés dans cette étude ont été simulés à l'aide des données réelles
obtenues avec le KIA Soul EV de l'IRH. Il est donc intéressant de pouvoir implémenter tous
ces algorithmes pour estimer tous les paramètres qui peuvent influencer la consommation
énergétique d'un véhicule électrique dans le modèle dynamique longitudinal tels que la
roulement dans différentes conditions de conduite comme en été ou les conditions sèches et
aussi pour les conditions hivernales c'est-à-dire à basse température, présence de neige et
D'autre part, il serait aussi intéressant d'implémenter ces algorithmes pour d'autres types
de véhicule comme les voitures à moteur thermique mais surtout les véhicules lourds. Cela
De plus, savoir comment le véhicule dépense de l'énergie est un atout. Il serait aussi très
intéressant de mettre au point une stratégie pour récupérer cette énergie dépensée. Un point
Références
[5] National Academies of Sciences, Engineering, and Medicine, «Review of the 21st
Century Truck Partnership: Third Report,» National Academies Press, 2015.
[9] R. Andersson, «Online estimation of rolling resistance and air drag for heavy duty
vehicles,» 2012.
[11] K. B. Singh and S. Taheri, «Estimation of tire-road friction coefficient and its
application in chassis control systems,» Systems Science & Control Engineering,
vol. 3, pp. 39-61,2015.
[12] Jadranko Matusko; Ivan Petrovié; Nedjeljko Peri, «Neural network based tire/road
friction force estimation,» Engineering Applications ofArtificial Intelligence, vol.
21, n° % 14, p. 442 456, 2008.
[13] Richard Barney Carlson; Henning Lohse-Busch; Jeremy Diez and Jerry Gibbs, «The
Measured Impact ofVehicle Mass on Road Load Forces and Energy Consumption
for a BEV, HEV and ICE Vehicle,» SAE International Journal of Alternative
Po wertrains, 2013.
[14] A. Alagumalai, «InternaI combustion engines: Progress and prospects,» Renewable and
Sustainable Energy Reviews, vol. 38, pp. 561-571,2014.
[17] A. Moretti, «The new electric vehicle "Renault Kangoo",» Electric vehicle symposium
E V's 17, 2000.
[19] L Martellucci; M Santoro; C Boccaletti, «A power train with planetary gear system:
Advantages and a design approach,» Proc. ofEVS, 2001.
[28] L. Li, F.- Y. Wang et Q. Zhou, «Integrated longitudinal and lateral tire/road friction
modeling and monitoring for vehic1e motion control,» IEEE Transactions on
Intelligent Transportation Systems, vol. 7, n° Il , pp. 1-19, 2006.
73
[29] A. M. A. Soliman, «Effect of Road Roughness on the Vehicle Ride Comfort and
Rolling Resistance,» chez 2006 World Congress, Detroit, Michigan, 2006.
[30] J. Barrand et J. Bokar, «Reducing Tire Rolling Resistance to Save,» SAE International,
vol. 1, nO 11, pp. 9-17, 2008.
[33] S. Pagerit, P. Sharer et A. Rousseau, «Fuel Economy Sensitivity to Vehicle Mass for
Advanced Vehicle Powertrains,» chez Advanced Hybrid Vehicle Powertrains
2006, Detroit, Michigan, 2006.
[36] S. Salisbury, «Cold Weather On-Road Testing ofa 2015 Nissan Leaf,» Idaho National
Laboratory, Idaho, 2016.
[41] C. Ahn, H. Peng et E. Tseng, «Robust estimation of road friction coefficient using
lateral and longitudinal vehicle dynamics,» Vehicle System Dynamics, vol. 50,
nO 16, pp. 961-985, 2012.
[43] Y. Chen et 1. Wang, «Adaptive Vehicle Speed Control With Input Injections for
Longitudinal Motion Independent Road Frictional Condition Estimation,» IEEE
Transactions on Vehicular Technology, vol. 60, n° 13, pp . 839 - 848, 201l.
[46] 1. Matusko, 1. Petrovié et N. Perié, «Neural network based tire/road friction force
estimation,» Engineering Applications ofArtificial Intelligence, vol. 21, n° 13, pp.
442-456, 2008.
[47] T. Zahid, K. Xu, W. Li , C. Li et H. Li, «State of charge estimation for electric vehicle
power battery using advanced machine learning algorithm under diversified drive
cycles,» Energy, vol. 162, pp. 871-882, 2018.
75
[49] C.-M. Tseng et C.-K. Chau, «Personalized Prediction ofVehicIe Energy Consumption
Based on Participatory Sensing,» IEEE Transactions on Intelligent
Transportation Systems, vol. 18, nO 111, pp. 3103 - 3113, 2017.
[54] H. K. Fathy, K. Dongsoo et 1. L. Stein, «Online vehicIe mass estimation using recursive
least squares and supervisory data extraction,» chez 2008 American Control
Conference, Seattle, W A, USA, 2008.
[55] Département de l'Énergie des États-Unis, «Tires and Passenger VehicIe Fuel Economy
: Informing Consumers, Improving Performance,» National Research Council,
2006.
[56] Département de l'Énergie des États-Unis, «Where the energy goes: Electric cars».
[61] José Marie P. Menezes Jr. et Guilherme A. Barreto, «Long-term time series prediction
with the NARX network: An empirical evaluation,» Neurocomputing, vol. 71,
nO 16-18, pp. 3335-3343, 2008.
APERÇU
Consommation de carburant :2.0 (Automatique ville), 2.6 (Automatique autoroute)
Moteur: 100% Électrique 81,4kW
Puissance: 109 cv @ 2730 rpm
Boîte de vitesses: Engrenage de réduction à une vitesse
Carrosserie : À hayon
Portières: 5; places: 5
SUSPENSION ET CONDUITE
Barre stabilisatrice avant
Direction à assistance électrique avec direction Flex
Indicateur de basse pression des pneus: Moniteur de pression des pneus (TPMS)
Pneus: P205/60R16
Suspension arrière à barre de torsion
Suspension avant indépendante
78
Abstrncl- The new strin gen t standards laid dow n for resulting in a better fuel economy. Therefore, the fuel
rf'ducing poUutiOD and greenhou.5t ga eDl15S10nS iD eonsumption of the vehicle and the re lea e of C0 2 in the
.r-.. nspo.otation !!cctor compd the ve hiclc mnnufacturcn to use almOsphere signi fieantly depends on lhe resistanee forces of
difftrent ttchoologies to dee,nrbonize tbe v('hicles and reduce the vchicle [9]. The drivi ng resistanee comp rises rolling
th cir cnergy consu mptiun. Powcrtruin contro l n pt imiz.~t i on , resistanee, air drag and mass. Aerodyoamic drag arises fro m
cspcdnlly bnscd o n instantancou nccura le information, bas nn the action of air in the opposite direction to the mo tion of lhe
impurta nt rult.' in imprO\'ing th e rud ecunumy uf "chide.
vehicle [ 10]. It is the mai n source of res istanee and
Howcvcr, its successful implcmcnhltion requins an occunte
proportional to the square of the vehiele ' s specd. At highway
pl.trnm e t~r estiml.tio n method. T he m:lÎn pu rpose orthis :.rticle
condi tions, over 50% of the power ofa car is uscd to overcome
is to propose an Recura te adaptive method for esti mariog tbe
ai r drag. For instance, when the specd of the ve hicle is 50
rolUng resÏ5tance .and auody namic drô.g, which are h vu
important cacton in powertraÎD op timization of a vebic1e. ln tbis km/h, the aerodynamic drag represents less than 50% of the
rt'gard, two onlin!' adaptive metbods, namety l'ecursive l e~ut total resistance while it inereases to more lhan 80% at 130
squares (RLS) and Kalman flller (KF), are formul. lecl for the km/h. Vehicle control systems req ui re lateral and longitudinal
parnmtttrs estimation of tht powtru'ain by UJing th!' conditions (vcloc ities and forces) to control wheel spin,
cxperimc ntal data from a rcal drivi ng cycle 1tnd thelr vehicle yaw raIeS, and lateral slip angles. Amo ng these states,
perfol'mnncts are compn red with an offiine tl'nined al'tifidal the estimation of longi tudinal states makes a major
neural nctwork., (A N1\'). T he cx~ri mc nt al and si mulation contribution to vehic1e stability and traction control [I I J,[ 12].
results show that KF iJ more lICCuflite tb'ID RLS in terms of
prcdicting th e nmd l'ulHnl( rcsistuncc and :lcrodynllmic dnl~
Accurdte real-ti me estimation of veh icle loads, such as
and il performance is very Deôtl' to the developed AN , vehicle mass and paveme nt conditions, aerodynamic drag and
rolling reSiSlanCe, can impro ve safety, effieieney and
Kt!)I,..'ortls- Aerody"amic forc~, ollli,,~ ~.\timatio", mlling perfo rmance [ 13]. The aerodynamie drag estimate has not
rl!sisfan ce, recllrsive itlentification, "Id .icle (Iylltlmie mot/el. been fu lly explored yet in the literature and req uires extens ive
investigations. Somc tesls have becn carri ed ouI in \Vi nd
1. INTROlJuc.-nO, tunne ls to si mul ate the rond conditions for aerodynamie drag
ealculation. Since aerodynamics has become essent ial in
The energy eonsum ption of eonventiona l vehie les depends
aut omotive des ign, the main development tool is the large-
on several paramcters, sueh as vehiele ' s weight, engine
seale wind tunnel. Unfortuna lcly, for a realistie simulation of
emeieney, and dyna mie parameters like rollin)l res istanee and
road conditions. the conventional large-scale wind tunne ls
r
ve hiele aerodyna mi e dral( Il and the ~ l ob.1 fuel savinl( in
have some limitations [1 4]. The ground is fixed relali ve to the
slich systems is thoroughly related to the power flow control
ve hicle, which allows a boundary layer to develop and the
between each eoergy source and the powertrain L2J,[3] . Online
estima tion of vehicle parameters has been an acti ve area of wheels of the development mule not to rotate [1 5]. ln addition,
researeh in rceent dee.des [4]. Duc to the diO"erencc in the size the aerodynamie data obtai ned From di fTerent wi nd tunnels
show a signi fieant degree of di . persion. Not onl y there arc
o fule vehieles, the form of the bodywork and the number of
wheel .xles, the rolling resistanee and air resistance vailles variations in the absolute dmg and the rolling res istance
values, but ulere is also a dispersion due to Ihe change of
vary From one vehicle to another [5]. However, lhese fac tors
config umtion [1 6]. In this regard, an independent and
are normally a sumed the same in the literatu re rega rdless o f
inexpensive method for identifying di lIerences in
tbe vehicle configuration. The ro lling resis~1nce eoem cient
aerodynamic dmg res ulting from various configurations
may vary depending on lhc Iype of li re and the road surface
should be devcloped. Online estimaIio n of vehiele para meters
[6],[7] ,[8]. Moreover, veh icle mass and ae rodynamie drag
has becn an acti ve area of research in recent years
coefficient can signi fi eantly change when a tra iler or anolher
[11 ], [1 7],[ 18]. The model-based approach (using lhe
un powered vehiele is attac hcd to the car. Therefore, onli ne
longitud inal dynamie model of the vehicle and the CA bus
parameler eslimation is necessary in the formulation of
signais) is the preferrcd method for estimating parameters duc
powcrtrain control in commercial vchiclcs si ncc il rcduccs the
initial ca libra tion cO"ort. Furthermo re, it helps tlle powertrain to the minima l additional sensors req uirements. The
estimaI ion o f the.e parameters has been the center of
controll er to quick ly adapt to di ffe rent ve hicle configurat ions,
attenti ons as they remarkably eOlltribute to fuel consumption
{
Fx = Frou + Fa.ra + FRrad • + Face
r;"; MgltCosll + i pA C..... (V)2 + Mgsln6+ Mil (1) 1 H ; lipA(v)2 Mgcos 61
x = lCaero ~lf
Authorized licensed use limilad to: UnIversite du Quebec a Tro is-Rivieres. Oownloaded on Seplember 30.2020 at 14:17'48 UTC from IEEE Xplore Restnctions apply.
82
1
'"
1000 2000 3000 4000 5000
Timo (s)
(b)
100
(a) (b)
KO
Fig. 2. al The tesling vebicle and bl Can-Bus-OBD
ê
~ 60
T ABLE 1. PARAMHERSOFTH E TESTING VEHICLE [22]. "~ 40
Panlmetu uni. Value ~ 20
Wld.h mm 1 747
H<!ah.
C
.,. .. mm
k'
1533
1547
0.375
O ~----------------------~
o 1000 2000
Time(s)
3000 4000 5000
~ ln "Dllalt O.OIR Fig. 4 . The long distance measurcd data profile (a)
Tire size P205170RI5 vehiele specd and (b) engine torque.
p kg/m' 1.200
These measured engine torque and vehicle speed values
2) Experimental tests analyses arc applied to the vchicle model to estimate rolling resistance
Fig. 3 presents the vehicle speed and engine torque values and aerodynamic drag coefficients by using the tbree proposed
during the real world shon distance test. methods.
(a) To test the robustness, each method has bccn applicd using
60 -------------------------, several reallife driving cycles for the long and shon distances
including the highway and urban driving context.
B. Experimental rest/lts
~40 Fig. 5 represents the simulation results obtained from the
~ sbon distance test. As it cao he seen in this figure, the
1l cstimatcd coefficients by difTcrcnt approaches are converging
ÎL20
en to the real values of the acrodynamie drag coefficient (Cu.... =
0.375) [22] and the roUing resistance coefficient for an asphall
road (II = 0.018). The value of ~I is assurned to bc the truc
50 100 ISO 200 250 measured value of the tire rol li ng resistance coefficient on the
Time (s) asphalt). These results iIIustrate thatthe cstimated coefficients
(b) by aU the methods are almost within the same range. However,
100
the estimated values by KF and ANN are more accurale than
RLS.
RO
~ 60
~
~
2- 40
~
20
o ~
o 50 100 150 200 250
Timc (s)
Fig. 3. The shon distance measurcd data profile (a) vehicle
speed and (b) engine torque. Time (s)
Authorized licensed use limited to: Universite du Quebec 8 Trots-Rivieres. Downloaded on September 30,2020 at 14:17:48 ure fr~ IEEE Xplore. Restrictions aPP'Y.
83
(b) TAULe II. ERROR RAT[ AND MEAN Of ESTIMM1!O VAllJ[S DURINO
SHORT OIST A~CE T~"S
ë
~ 0.0181 1 ~
E 1 - Reference ( ·~mdent.ll Aerody .... mlc drllg RuUine rl.~lsbtnc~
~om51
-- KF (long distance test), it is c lear that the RLS pcrfomlance has
cxperienccd sorne impro vements comparcd to the short test.
~E 1
_<1"
The results of the performed sim ulation show that the
1
parameters vary wi th time but nOl significant ly according to
~ 0.375 \"-~---'~---------- the tbree methods on each driving cycle. By linearizing the
~ iden ti fication modcl, th c evaluated mcthods can reasonably
.;: converge loward their nominal va lues. The accuracy of the
1000 2000 3000 4000 5000
Timc(s)
RLS and KF est imations is strongly dependcnt on the torque
va lues. However, the accuracy of the ANN high ly dcpends on
the vehiele spccd and thc cngine torque. Since the
(b)
measuremeOlS of these values are always prone to the
0 .01 8 1 oecurrence of error, they have a signi fi cant efTect on the
- Refe rence
ANN estimation results. The comparison of the results between RLS
ï1
~ O . ()I ~ 05 --Kr
It l..S and KF methods shows that the KF quality of estimations is
better than RLS. Moreover, the convergence oflhe KF method
~ is fas ter than the RLS. Regarding the error mte, the ANN
j U.UIX 1---=-,..,- ""..-------- shows a lower va lue than the other two methods, followed by
.~ ~'i" KF .
"" 0 .01795
.ê ln brief, the accuracy of the estimation depends on the
~ measuring instruments and the configuration of the methods.
ln addition. RLS and KF methods are more suitable fo r real·
0 .01 19 time monitoring than neural nctwork s ince they can easily
o 1000 2000 3000 4000 5000
Timc( O\) adapt to the d rills owing the variation of the conditions.
Fig. 6. Estimatcd coefficients during long distance test. al
aerodynamic drag and bl ro lling resistance .
Authorfzed Ncensed use limited 10: Universite du Quebee a Trois-Rivieres Oownloaded on September 30 ,2020 at 14'17:48 ure from IEEE Xpk>re Restrictions appfy_
84
V . CONCLUS ION [II J A. Vuhidi. A. Stelhnopoulou. and II. Peng. "Recursive ICOlSt squares
with forgctting for online estimation of vehicle mass and rood grade:
ln tbis work, tbree methods, namely RLS, KF, and ANN, theory and cxperiments," Vellicle System Dynomics, vol. 43, pp. 31 -5 5,
are uscd 10 simullancously esti ma te the aerodynamic drag and 2005.
rolling resistance coefficients under normal conditions. RLS r121 Vincent Winslead et llyn V. KolmaJ1o\'S.ky, "Estimation ofRood Grode
and KF conduct lhe eslim.1lion based on Ihe longiludinal and Vehielc Mas! vin Model Prcdieti\'c Control," Proceedings of the
1005 IEEE Conference on Control Applica/ÎOIu," Toronto, Canada,
model oflhe vehicle while ANN is independenl Of lhis mode!.
AugnSl 28-3 1,2005
The eSlimalion results of lhe melhods arc compared wilh
[13] Kidambi , N., Hume, R.• rujii, Y., l'ictron, G. ct n!., "Methods in
respect to the attained precision, convergence speed, and error Vchiele Mas.'\ and Kood Grade E.'\timation", SAE Inl. J. Pa.uerrg. Û lrs ,
rate. The final oulcomes of Ihis paper ind ic'Ie lhal the cn2014
cstimated values of the aerodynamic drag and ro lling [ 14] A. VahKli, M, Dnv.hiniro, A. StcranOJ)Oulou, and P. Huci,
resistance coefficients by these mcthods arc similar to those of "Simultancous mass and lime-v3rying grade estimation for hcavyduty
lhe nominal values. ln fact, the results of the tbree methods vehides." ;" Amer"-'on Control Conferenœ. 2003. Procccdings or the
show an acceptable error. However, KF, as an online 2003, , pp. 4951-4956 vol.6, 2003
[151 Ko~t , H., White, R., and Metz, L ," Road Eva luation of the
estimutor, has been more accu raIe and pcrformed more
Acrodynamic Charnctcristics of Hca vy Trucks", jOllnlt.rl SAE
approprialely under the normal road conditions. Tech"ica/, 2007
Usually, manufaclure gives these values for the nonnal [l6J K. Zarringh.lam, A. KClllci,n. W. Mclck, Il. Khajcpoor, S.-k. Chen,
opcraling condilion (no snow, no ice, sunny days, ClC.). The and N. Moshchuk, "A comparntivc study on identification or vchide
Înct1ial parameLcrs," in 10 / 1 American Conlml Conjf!f'f:llt'f! (ACC), •
proposed melhods can allow 10 eSlnTh1le these parameler under pp. 3599·3604: IEEE, 2012
difTcrenle opcraling including win ter conditions. (17] H. K. Fathy, D. Kang, and J. L Stein. "Onli ne vmide mass estimation
ln add ition, for eleclrieal vchic1e, those methods could bc sing recursive lcast squares and supcrvisory data extraction ," in 2008
AmerÎc:an collirol Cf.Hrjeren,'e, 2008, pp, 1842- 1848: lEEE.
applied to esti mate energy usage and also to predict the
[l8J R. R.jamani, G. Phanomchoeng, D. Piyabongkam, ,nd J. Y. l.ew,
required energy to rcach a des tination. ..Aigorilhrns for Real-Time Estimation of Individual \Vheel Tire- Rond
Friction Coeflicients." MechlltfOnicJ, IEEE/ASME Transllctio"s 0 11,
voL 17, pp. 11 83-1195,2012.
ACKNOWU, DG EMENTS [l9J R. Zarringhalam, A. Rezaei.n. W. Melek. A. Khajepour, S.-k. Chen.
and N. Moshchuk, liA compar.uive study on identification of vebic1e
This work was supported in part FRQNT (Fonds dc incrtial pammctcrs," in 1011 American COlltrol Confere"ce (ACCJ, ,
Recherche du Québec - Nature et technologies) in part by the pp. 3599-3604: IEEE, 20 12
Natural Science and Engineering Roseareh Counc il of r20] O. Trigui, Y. Dubc, S. KclouwanÎ. and K.. Agbossou, "Comparative
Canada. Est imation of Elcctrk: Vchicle Rolling Resistance Cœfficicnt in
WinKT Conditions," in 10/6 IEBI:: Yehit:le Power and P,.o/JlIl,viol1
C{J"Je/'f!"'~ (VPPC) , , pp. 1-6: IbEE, 2016
(21) M. K.andidayeni. A. Macias. A. Amamou, L Boulon. S. Kelouwani.
R EFERENCES and Il. Chaoui, "Ovcrview and bendurusrk. anaIysis of fuel cell
paramctCf'S estimation ror cncrgy management purposcs." Journal of
[ 1} F.. National Academies of Sciences and Medicine. Hel'iew of/he 11.~·1 pmt'ttr s(wrce~·. vol. 380, pp. 92-104. 20 18,
Cenlllly Tnu'/i PClrI"erxltip: n,ira N(1)Ort . National ACl1dcmics Press, [221 Dodge C. liber 2012. FlCffE TECHNIQUE POUR LE CANADA, 2012
20 t5 .
[2J K. Maalcj, S. Kclouwani, K. Agbossou, Y. Dubé, .nd N.flcn.o, "l.ong-
Irip optimal energy planning with online mass estimation for bauery
c lcclric vchiclcs." )o"nlOl o/Po .....erSollrces, vol. 64, no. Il . pp. 4929-
4941 , 2015.[I-5J
[3J K. Ma.lej. S. KelouW'Jni, K. Agbossou .• nd Y. Dubé, "Enhanced fuel
cell hybrid clettric vehicle power sharing method based on fuel cost
and mass cstimation." Jounloi nf Power Smm:c\', vol. 24R, pp. 668-
678, 2014.
[41 D, Zhang. A. (vanco, and Z, Fitipi. "Modd-Based Estimat ion of
Vehicle Acrodynamie Drag and Rolling Resistance," SAE
IntenultioIJal Jollnwl of CQmmen:i,,1 Yehicles, vol , 8. no, 2015.01-
2776, pp. 433-439, 2015.
r51 K. Andcrsson, "Online estimation ofmlling rcsist3ncc and air drag for
hcavy dnly vchiclcs," cd, 2012.
[6J H. Taghal'ijùr, A. Marc/lmi. H. KarÎm Maslak. Ulrd H. Ka/blchanl .
..Anificial Neural Nctwork estimalion of wheel rolling resistance in
c1,y lcam soi!.· Applied Soft Compl/ril/g. \'0/. 13. pp. 3544-355/. 8//
20/3.
Pl K. B. Singh and S. Tahcri, "Estimalion of tire rond friction cocflicicnt
and its upplicution in chassis control systems." Systemç Science &
Co"trol Engineering. vol. 3, pp. 39--61, 2015.
[8] J. Matmko, l. l'etrO\;é . and N. l'crié, "Neural nctwork bascd tire/rond
friction force estimation," Engi"eerilW Applicationv of ArlifidlJ
Intelligence. vol. 2 1, pp. 442 456, 4//2008.
[9J 1. Gibbs. II. Lohs-Busch, R. CBrlsoll. "The Mca5Ul'Cd Impact ofVebic1e
Mas!' on Road Load Forces and F.ncrgyConsumption for 3 BEV. HEV,
and ICF. Vchicl(!", SAF. IIIJenra/ÏfJlwl Jounral (if Ahemul;W!
PO"l!rlm;R\·.2013
[101 Shu-Qing Ynng" Donghoug Ding. "Dr:lg Rcducing Fto"~ in l.:uninar-
Turbulent Tr,msition Region", JUlInllI/ of FlllitL'i E,r,:i"eerinx. 2014
Authorized licensed use limited to: Universite du Quebec a Trois-Rivieres. Downloaded on Seplember 30,2020 al 14:17:48 UTC trom IEEE Xplore . Restrictions apply.
85
(a) (b)
Étape d'installation:
https://github.com/JejuSoul/OBD-PIDs-for-HKMC-EVs
3. Après avoir téléchargé les fichiers, transférer sur l'appareil Android dans le
5. Cliquer sur l'icône avec les 3 points en haut à droite après sur « Add predefined
6. Cliquer sur un fichier et répéter l'opération pour ajouter tous les fichiers voulus.
Étape d'utilisation:
86