FR 3093254 A1 20200828 FR 1901751 A 20190221 Description Titre de l'invention : Procédé d'estimation du couple d'une machine électrique synchrone [0001] La présente invention se rapporte de manière générale au domaine de l'électrotechnique et plus précisément au pilotage d'un moteur électrique, trouvant notamment une application avantageuse dans le domaine automobile. [0002] La commande d'un moteur électrique utilise classiquement une boucle de régulation, dans laquelle étant donné un couple de consigne, on en déduit des courants de consigne et des tensions de commande pour l'onduleur alimentant le moteur électrique. Dans certains procédés de régulation, on corrige l'erreur entre le couple réalisé par le moteur électrique et le couple de consigne. [0003] De plus la connaissance du couple réalisé est également nécessaire pour surveiller le bon fonctionnement du moteur électrique, notamment pour interdire par exemple l'utilisation de celui-ci lorsque le couple réalisé ne suit pas le couple de consigne du fait d'un défaut logiciel et risque d'engendrer des accélérations ou freinages intempestifs dangereux en conduite, et de mettre en cause les lois de gestion énergétique pour les véhicules hybrides. [0004] Dans de tels procédés de régulation du couple ou de surveillance du couple, le couple réalisé est soit mesuré soit estimé. La mesure du couple réalisé nécessite un capteur coûteux, une estimation du couple réalisé est donc préférable. Une telle estimation utilise généralement des mesures des courants de phase du moteur électrique, ainsi que des valeurs représentatives du flux magnétique généré par les bobinages du moteur électrique. [0005] Or le flux magnétique généré par le moteur électrique, qui peut être déduit des mesures de courants de phase, ne varie linéairement avec ceux-ci que pour de faibles valeurs d'intensité, et est donc difficile à estimer, notamment dans un moteur électrique qui produit un fort couplage magnétique. Par exemple la construction de cartographies embarquées de composantes du flux magnétique calculées à partir d'un modèle quasi statique non linéaire permet une estimation du couple, mais celle-ci présente une erreur d'estimation importante avec le couple réellement mesuré sur banc. D'autres méthodes permettent d'estimer les composantes du flux magnétique en utilisant les mesures des tensions statoriques, mais l'erreur d'estimation reste importante du fait des non-linéarités de l'onduleur alimentant la machine, qui rendent ces mesures de tension très imparfaites. [0006] Il est également possible, afin de rester dans un domaine linéaire de variation du flux en fonction des mesures des courants de phase, de surdimensionner le moteur 2 électrique par rapport au besoin de couple. Néanmoins cette solution n'est pas souhaitable dans le domaine automobile où on cherche le moins d'encombrement possible du groupe motopropulscur électrique. [0007] Le brevet FR3054391 propose une estimation du couple d'un moteur électrique utilisant un unique réseau de neurones en entrée duquel on applique un vecteur d'intensité comportant les mesures des courants statoriques et rotorique. Du fait de la relation non linéaire qui existe entre le couple d'une part et les courants d'alimentation ainsi que les composantes du flux magnétique d'autre part, l'apprentissage des paramètres du réseau de neurones ne donne pas de bons résultats, et l'erreur d'estimation du couple en sortie est conséquente. [0008] Un des buts de l'invention est de remédier à au moins une partie des inconvénients de la technique antérieure en fournissant un procédé et un estimateur du couple d'une machine électrique synchrone, qui soit précis tout en ne nécessitant pas d'embarquer beaucoup de ressources en mémoire et en calcul, ni d'effectuer en amont un réglage complexe et long tel que celui d'un réseau de neurones. [0009] A cette fin, l'invention propose un procédé d'estimation du couple d'une machine électrique synchrone comportant : [0010] - une étape de réception de mesures des courants statoriques alimentant ladite machine, et d'une valeur mesurée de la température de ladite machine électrique, [0011] - une étape de détermination d'une composante directe desdits courants statoriques, d'une composante en quadrature desdits courants statoriques alimentant ladite machine, et d'une valeur de flux rotorique correspondant à ladite valeur mesurée de température, ladite valeur mesurée de température étant comprise entre une première température prédéfinie et une deuxième température prédéfinie, [0012] - une étape d'obtention d'une première valeur de flux global équivalent cor- respondant à une valeur de flux rotorique à ladite première température de ladite machine électrique, dans une première cartographie donnant ladite première valeur de flux global équivalent en fonction de ladite composante directe desdits courants statoriques et de ladite composante en quadrature desdits courants statoriques, [0013] - une étape d'obtention d'une deuxième valeur de flux global équivalent cor- respondant à une valeur de flux rotorique à ladite deuxième température de ladite machine électrique, dans une deuxième cartographie donnant ladite deuxième valeur de flux global équivalent en fonction de ladite composante directe desdites courants statoriques et de ladite composante en quadrature desdits courants statoriques, [0014] - une étape d'interpolation utilisant ladite première valeur de flux global équivalent et ladite deuxième valeur de flux global équivalent, et fournissant une estimation du flux global équivalent correspondant à ladite valeur de flux rotorique à ladite valeur mesurée de température, 3 [0015] - une étape d'estimation du couple de ladite machine en utilisant la formule suivante : [0016] [Math.1] Tel phiD(Q 3p [0017] où : [0018] - Tq est le couple de ladite machine électrique, [0019] - p est le nombre de paires de pôles de ladite machine électrique, [0020] - Iq est ladite composante en quadrature desdits courants statoriques, [0021] - 4)1-cf) est ladite valeur de flux rotorique à ladite valeur mesurée de température Trot de ladite machine électrique, [0022] - et phiD4f(T,3) est ladite estimation du flux global équivalent issue de ladite étape d'interpolation. [0023] Grâce à l'invention, l'estimation du couple nécessite des cartographies pour seulement un flux global équivalent et non pas pour chaque composante directe et en quadrature du flux magnétique. De plus ces cartographies sont obtenues par la formule [Math 1] ce qui permet d'utiliser des mesures de couple effectuées sur un banc de test. Les résultats sont donc plus précis que dans l'art antérieur. [0024] Avantageusement, le procédé d'estimation du couple selon l'invention comporte une étape préalable d'établissement d'au moins ladite première cartographie et ladite deuxième cartographie, ladite étape préalable comportant les sous-étapes suivantes : [0025] - mesures de couple Tq de la machine électrique réalisées pour différents triplets (Id, Iq, Ti) correspondant à différents points de fonctionnement de ladite machine électrique, Id étant ladite composante directe desdits courants statoriques et T, la température du rotor, 100261 - calcul pour chacun desdits triplets d'une composante de flux (I), indépendante de la température du rotor et définie par (I), = (Ld - Lq)*Id , où Ld est la composante directe de l'inductance statorique et Lq est la composante en quadrature de l'inductance statorique, en utilisant F expression 100271 IMath.21 Tq 3P f(n) 2 [0028] (1.)f(T,) étant la valeur de flux rotorique correspondant à la température T, du rotor, [0029] - construction de ladite première cartographie en additionnant, pour différentes valeurs Id de composantes directes de courants statoriques et Iq de composantes en quadrature de courants statoriques d'un desdits triplets, ladite composante de flux (hl précédemment calculée, et la valeur de flux rotorique à ladite première température, 0030] - construction de ladite deuxième cartographie en additionnant, pour différentes 0 valeurs Id de composantes directes de courants statoriques et Iq de composantes en quadrature de courants statoriques d'un desdits triplets, ladite composante de flux (Dl précédemment calculée, et la valeur de flux rotorique à ladite première température. [0031] Cette façon de construire les cartographies permet de s'affranchir d'une construction à partir de points de fonctionnement de la machine pour lesquels on a une température fixe égale à la première ou la deuxième température prédéterminée, ce qui est difficile à garantir. De plus il est difficile de faire des mesures à des températures extrêmes de machine, à cause de l'endommagement possible de la machine. Il est plus facile de calculer le flux rotoriquc pour une température prédéterminée. [0032] Selon une autre caractéristique avantageuse de l'invention, ladite première tem- pérature et ladite deuxième température correspondent à des températures limites de ladite machine électrique. Cela permet de n'embarquer dans un véhicule dont le groupe motopropulscur utilise le procédé selon l'invention, que deux cartographies du flux global équivalent, économisant ainsi des ressources mémoire. [0033] L'invention concerne aussi un estimateur de couple d'une machine électrique synchrone à aimants permanents, mettant en oeuvre le procédé d'estimation du couple d'une machine électrique selon l'invention. [0034] L'estimateur de couple selon l'invention présente des avantages analogues à ceux du procédé d'estimation du couple selon l'invention, et est par exemple utilisé dans un véhicule électrique ou hybride. [0035] D'autres caractéristiques et avantages apparaîtront à la lecture d'un mode de réa- lisation préféré décrit en référence aux figures dans lesquelles : [0036] [fig. Il représente des étapes du procédé d'estimation du couple d'une machine électrique selon l'invention, dans ce mode de réalisation préféré, [0037] [fig.2] représente une étape préalable à l'utilisation du procédé selon l'invention, dans ce mode de réalisation préféré, et [0038] [fig.3] représente un estimateur de couple selon l'invention utilisant le procédé d'estimation de couple selon l'invention. [0039] Selon un mode préféré de réalisation de l'invention représenté à la figure 1, le procédé d'estimation du couple selon l'invention comporte des étapes El à E6. Il est mis en oeuvre au moins en partie dans un calculateur muni de ressources mémoire, par exemple un calculateur embarqué dans un véhicule électrique ou hybride comportant une machine électrique ME, représentée sur la figure 3. [0040] La machine électrique ME est, dans ce mode de réalisation préféré de l'invention, une machine synchrone à aimants permanents. Cependant le procédé selon l'invention est utilisable sur d'autres types de machines électriques, par exemple la machine synchrone à rotor bobiné où 4f(T) est remplacé par M*If avec M l'inductance mutuelle et T- le courant d'excitation. [0041] L'étape El est la réception de mesures des courants statoriqucs il, i2, i3 alimentant la machine électrique ME, et de la mesure de la température Trot de la machine électrique ME. La température Tm, est prise, dans ce mode de réalisation de l'invention, par un capteur de température au niveau du rotor de la machine électrique ME. En variante la température de la machine est mesurée au niveau du stator et la température du rotor utilisée par le procédé est estimée à partir d'un modèle thermique. Les mesures de courant sont prises par des capteurs de courant au niveau des phases de la machine électrique ME. [0042] L'étape E2 est la détermination d'une composante directe -Id et d'une composante en quadrature 4 des courants statoriqucs de la machine électrique ME, par une transformation des courants des phases il, in et i' mesurés depuis un repère fixe lié au stator, vers le repère de Park, en utilisant la transformée de Park bien connue de l'Homme du métier. En variante la transformée dq0 est également utilisable. [0043] Dans cette étape E2, la valeur (pf(T,m) de la valeur de flux rotorique correspondant à la température Torprécédemment mesurée est également déterminée, en utilisant par exemple une cartographie donnant le flux rotorique en fonction de la température du rotor. En variante, la relation qui lie la température des aimants au rotor au flux rotor est définie par une droite affine cpf(T1) = cl)f(20°C) + A*(T1-20). On détermine la valeur de flux rotorique f(20°C) à 20°C (degré Celsius) et la constante A par différents essais à vide, en imposant différentes températures rotor connues. De préférence on utilise plusieurs droites affines pour plus de précision, par exemple une sur la plage [-20°C, 40°C] et une sur la plage [40°C, 150°C]. Ces plages peuvent bien sûr être modifiées, par exemple la première peut être étendue à [-20°C, 60°C]. [0044] L'étape E3 est l'obtention dans une première cartographie, d'une première valeur de flux global équivalent phiDM)f(T1)) en fonction des composantes directe et en quadrature des courants statoriques déterminées à l'étape E2, cette première valeur de flux global phiD(f(T I)) étant valable pour une valeur de flux rotorique à une première température prédéfinie T,. Plus précisément, sachant que le couple Tq de la machine électrique est donné par la formule : [0045] [Math.3] Tg = 3f .(Laix Phif 'xJq xId) [0046] où : [0047] Ld est la composante directe de l'inductance statmique de la machine électrique ME, [0048] p est le nombre de paires de pôles de la machine. 6 [0049] Phif est le flux des aimants du rotor vu par le stator [0050] et Lq est la composante en quadrature de l'inductance statoriquc de la machine électrique ME, [0051] On peut donc formuler le couple Tq dc la machine comme : [0052] [Math.4] 3p 3p Tg - 2 . ((L d - g) x Id + Phif).1q 100531 où phiD(()4T,t)) est le flux global équivalent pour une valeur du flux rotorique cor- respondant à la température de la machine mesurée Trot. 100541 Dans cette troisième étape E3, le flux global équivalent phiD(4)(1-1)) déterminé ne correspond pas au flux global équivalent phiD(4)4T,t)) correspondant à la température machine mesurée Trot, mais à celui qu'aurait la machine électrique ME si sa température était égale à la première température prédéfinie TI, pour les mêmes valeurs de courants statoriques. 100551 L'étape E4 est de façon similaire, l'obtention dans une deuxième cartographie, d'une deuxième valeur de flux global équivalent phiD(4)4T2)) en fonction des composantes directe et en quadrature des courants statoriques déterminées à l'étape [2, cette deuxième valeur de flux global équivalent phiD(4)(T2)) étant valable pour une valeur de flux rotorique à une deuxième température prédéfinie 112, différente de la première température prédéfinie TI. De préférence T1 et T2 correspondent à des températures limites de la machine électrique ME, de sorte que la température mesurée Tru, soit toujours comprise entre les températures prédéfinies TI et T2. 100561 L'étape E5 est une interpolation linéaire utilisant les deux valeurs de flux global équivalent phiD(4)4T1)) et phiD(4)r(T2)) obtenues aux étapes E3 et E4, pour estimer une valeur de flux global équivalent phiD(4),(Tr')) valable pour le flux rotorique à la température mesurée T'' de la machine électrique ME. L'étape E5 utilise par exemple, pour le calcul de phiD(4)4T1')), la formule : [0057] [Math.5] phin(tpt(Tv',))= phiD(g o,,(1.1))-phiD(cp i(12)) ( T-) - (1- ') "up r(7"7')-9/(T2))+ phiD(pr(T2)) t" [0058] L'étape E6 est l'estimation du couple Tq dc la machine électrique en utilisant la formule : [0059] [Math.1] 3p ( ro t))etiei 2Tq='phiD( [0060] où : [0061] - p est le nombre de paires de pôles de ladite machine électrique, [0062] - Iq est la composante en quadrature (lesdits courants statoriques déterminée à l'étape 7 E2, [0063] - et phiDttf(Tio,)) est l'estimation du flux global équivalent issue de l' étape d'interpolation E5. [0064] En référence à la figure 2, une étape EO préalable à l'utilisation du procédé d'estimation de couple selon l'invention est représentée sous la forme d'étapes E01 à E04, permettant d'obtenir les première et deuxième cartographies donnant respectivement le flux global équivalent phiD(tf(T,)) et le flux global équivalent philD((pf(T2 )) en fonction des composantes directe et en quadrature des courants statoriques de la machine électrique ME. [0065] L'étape E01 est la mesure des valeurs Tq de couple de la machine électrique ME à différents points de fonctionnement de la machine. Cette étape est réalisée sur un banc de tests de la machine électrique ME. Pour chaque point de fonctionnement testé, on mesure et on mémorise dans cette étape E01 le couple Tq réalisé par la machine électrique ME, les valeurs des courants des phases de la machine électrique, et la température Ti du rotor. On associe ainsi pour différentes valeurs de couple Tq, différents triplets (Id, Iq, Ti) comportant la composante directe des courants statoriqucs correspondants, la composante en quadrature des courants statoriqucs correspondants, et la température Ti de la machine correspondante. [0066] L'étape suivante E02 est le calcul, dans un calculateur, pour chacun des triplets (avec le couple associé) mémorisés à l'étape E01, d'une composante de flux 0, indépendante de la température du rotor et définie par (I), = ([AI - Lq)*Id, où Ld est la composante directe de l'inductance statorique et Lq est la composante en quadrature de l'inductance statorique, en utilisant 1' expression : [0067] [Math.2] 14 3 919 [0068] où [0069] -Tq est le couple associé au triplet (Idaq,Ti) du point de fonctionnement considéré, [0070] -p est le nombre de paires de pôles de la machine, [0071] -et 4f(T) est la valeur de flux rotorique correspondant à la température Ti du rotor, obtenue par exemple dans une cartographie donnant le flux rotorique en fonction de la température de la machine. [0072] L'étape E03 est la construction de la première cartographie en additionnant, pour chacun des points de fonctionnement considérés, la composante de flux 4:131 calculée à l'étape E02, avec la valeur de flux rotorique correspondant à la première température prédéfinie TI, elle-même obtenue par exemple dans une cartographie donnant le flux rotorique en fonction de la température de la machine. En revenant à la formule [Math (1) . I( Ti) 8 4], on obtient bien des valeurs de flux global équivalent phiD(4pf(T1)). En associant chaque doublet de composante directe et de composante en quadrature des courants statoriques d'un point de fonctionnement à la valeur de flux global équivalent correspondante obtenue lors de cette étape E03, on obtient la première cartographie. [0073] L'étape E04 est la construction de la deuxième cartographie de manière similaire, en additionnant, pour chacun des points de fonctionnement considérés, la composante de flux (DI calculée à l'étape E02, avec la valeur de flux rotorique correspondant à la deuxième température prédéfinic T,, ce qui donne une valeur de flux global équivalent phiD(cpr(T2)). En associant chaque doublet de composante directe et de composante en quadrature des courants statoriques d'un point de fonctionnement à la valeur de flux global équivalent correspondante obtenue lors de cette étape E04, on obtient la deuxième cartographie. [0074] Il est à noter qu'il est également possible d'obtenir directement des première et deuxième cartographies de flux global équivalent en utilisant les mesures issues de l'étape E01 respectivement à la première température prédéfinie T, et à la deuxième température prédéfinic T2, en utilisant la formule [Math 1]. Néanmoins cela demande de stabiliser la température de la machine électrique ME respectivement aux températures T, et T2 pour différents points de fonctionnement, ce qui est très difficile. [0075] La figure 3 illustre un estimateur de couple EST mettant en oeuvre le procédé d'estimation du couple de la machine électrique ME selon l'invention. L'estimateur EST reçoit en entrée des mesures des courants des phases de la machine électrique ME en étant relié à des capteurs de courant de la machine électrique ME, ainsi qu'une mesure de température en étant relié à un capteur de température du rotor. Il fournit en sortie une estimation Tq du couple de la machine électrique ME. Ce couple Tq est soustrait à une consigne de couple Tref, l'erreur E résultant de cette opération étant envoyée en entrée d'un moyen de commande CMD de la machine électrique ME. Le moyen de commande CMD comporte typiquement un régulateur et un onduleur apte à appliquer des tensions de commande aux enroulements statoriques de la machine électrique ME. [0076] 11 est à noter que bien que dans ce mode de réalisation de l'invention, on utilise seulement deux températures prédéfinies pour réaliser l'interpolation du flux global équivalent à la température mesurée, en variante de réalisation on utilise plus de deux températures prédéfinies et on choisit celles qui sont les plus proches de la température mesurée pour faire l'interpolation du flux global équivalent à la température mesurée. Cette variante permet un résultat plus précis mais au prix de plus de ressources mémoire. De plus l'ordre des étapes n'est pas figé dans la réalisation de l'invention, les étapes E3 et E4 pouvant être effectuées en parallèle ou l'une après l'autre par exemple. 9 [Revendication 1]