FR 3122298 A1 20221028 FR 2104317 A 20210426 La présente invention concerne un procédé de détermination du couple (Te) d’une machine électrique, dans lequel on mesure (MES) les courants et les tensions, et on détermine la vitesse de rotation électrique du rotor de la machine. Puis, au moyen d’un modèle dynamique (MOD) du flux magnétique et d’un filtre de Kalman sans parfum, on détermine le flux ( ) puis le couple (Te) de la machine électrique. Figure 1 à publier Procédé de détermination du couple d’une machine électrique La présente invention concerne le domaine du contrôle des machines électriques, en particulier le contrôle des machines électriques synchrones à pôles saillants. Ces machines électriques trouvent une application notamment dans le domaine des véhicules automobiles. Classiquement, une machine électrique comprend un rotor (partie mobile) et un stator (partie fixe). Le rotor est habituellement logé à l'intérieur du stator. Généralement, le stator est de forme annulaire et est logé à l'intérieur d'un support tubulaire pour y être fixé. Le stator comprend des générateurs de flux magnétique, généralement des bobinages électriques. Ces bobinages sont alimentés par une pluralité (classiquement trois) phases électriques, afin de générer un champ magnétique tournant. De plus, selon le type de machine électrique, notamment pour les machines électriques synchrones à pôles saillants, le rotor peut comprendre des aimants permanents. Lors du fonctionnement d'une telle machine électrique, les bobinages du stator sont parcourus par un courant électrique pour générer le champ magnétique nécessaire à l'entraînement en rotation du rotor. Dans le domaine de machines électriques, la machine synchro-réluctante assistée d’aimants permanents est aujourd’hui de plus en plus développée dans de nombreuses applications industrielles en raison de son facteur de puissance, de son efficacité, de sa compacité et de sa capacité à fonctionner à haute vitesse. Grâce à l’aimant permanent dans la construction du rotor, le facteur de puissance est amélioré et les pertes du stator sont considérablement réduites par rapport aux machines synchro-réluctantes. De plus, le couple de réluctance créé dans une telle machine limite le besoin en aimant permanent, qui reste encore onéreux, et rend donc ce type de machine moins chère à produire. Dans plusieurs applications industrielles, l’information précise du couple est indispensable pour les machines électrique que ce soit pour le contrôle du couple, notamment avec les stratégies MTPA ou MTPV (respectivement de l’anglais « Maximum Torque per Ampere » pouvant être traduit par couple maximum par ampère, et « Maximum Torque per Voltage » pouvant être traduit par couple maximum par volt), pour la validation au banc moteur, pour le diagnostic et la détection des fautes, et aussi pour la surveillance et la sécurité. Les informations de couple peuvent être obtenues au moyen de mesures par des capteurs de couples, mais ces capteurs sont souvent coûteux, contiennent beaucoup de bruit, et sont difficilement utilisables en temps réel. Le document « P. Sue, D. Wilson, L. Farr, and A. Kretschmar, “High precision torque measurement on a rotating load coupling for power generation operations,” in 2012 IEEE International Instrumentation and Measurement Technology Conference Proceedings. IEEE, 2012, pp. 518–523. » décrit une telle détermination du couple au moyen d’un capteur. Une autre possibilité pour déterminer le couple d’une machine électrique consiste en l’utilisation d’algorithmes qui donnent une estimation du couple à partir d’autres mesures telles que les courants, les tensions et la vitesse de rotation du rotor de la machine électrique. Par exemple, le document « K. C. Yeo, G. Heins, and F. De Boer, “Comparison of torque estimators for pmsm,” in 2008 Australasian Universities Power Engineering Conference. IEEE, 2008, pp. 1–6. » décrit une estimation du couple au moyen d’une équation générale de couple, en supposant que les inductances et le flux d’aimant permanent sont tous connus et constants. Cette approche est simple à appliquer, mais son manque de robustesse est un inconvénient majeur en raison de son hypothèse non réaliste (inductance et flux magnétiques connus et constants). Une autre solution décrite dans le document « J. Xu, S. Panda, Y. Pan, T. Lee, and B. Lam, “Improved pmsm pulsating torque minimization with iterative learning and sliding mode observer,” in 2000 26th Annual Conference of the IEEE Industrial Electronics Society. IECON 2000. 2000 IEEE International Conference on Industrial Electronics, Control and Instrumentation. 21st Century Technologies, vol. 3. IEEE, 2000, pp. 1931–1936. » concerne une technique de mode glissant pour faire l’estimation du couple. Toutefois, cette méthode ne concerne que les machines électrique non saillantes, et les inductances de la machine électrique doivent être connues avec précision a priori, ce qui n’est pas toujours facile en pratique. L’utilisation d’un observateur de flux permet également la détermination du flux. Le document « M. Koteich, “Flux estimation algorithms for electric drives: a comparative study,” in 2016 3rd International Conference on Renewable Energies for Developing Countries (REDEC). IEEE, 2016, pp. 1–6. » illustre cette méthode. Un certain nombre d’observateurs de flux sont principalement développés à partir d’un modèle de tension qui comprend l’intégration de signaux de force électromotrice et l’estimation de flux à travers le modèle de courants. En général, la force électromotrice intégrée est traitée par un filtre passe-haut, alors que l’estimation du flux passe à travers un filtre passe-bas. La bonne qualité de l’estimation du flux dépend du choix de la fréquence de coupure, en particulier pour les faibles vitesses de rotation de la machine électrique. Toutefois, un filtre qui peut rapidement atténuer le décalage constant pour des vitesses moyennes et élevées peut être imprécis à faible vitesse et inversement. En outre, la plupart des observateurs du flux en général ont besoin d’informations a priori précises sur les paramètres de la machine électrique tels que la résistance du stator, les inductances et le flux d’aimant permanent. Toutefois, ces paramètres ne sont généralement pas connus exactement dans la réalité, notamment le flux d’aimant permanent en raison du changement de température du rotor et les inductances en raison de saturation magnétique. Par conséquent, il est nécessaire de concevoir un observateur de flux qui prend en compte des incertitudes de paramètres de la machine électrique et aussi des incertitudes des mesures de capteurs. L’invention a pour but de déterminer de manière précise et robuste le couple d’une machine électrique. Pour cela, la présente invention concerne un procédé de détermination du couple d’une machine électrique, dans lequel on mesure les courants et les tensions, et on détermine la vitesse de rotation électrique du rotor de la machine. Puis, au moyen d’un modèle dynamique du flux magnétique et d’un filtre de Kalman sans parfum, on détermine le flux puis le couple de la machine électrique. L’utilisation d’un filtre de Kalman sans parfum permet de prendre en compte les non linéarités du modèle dynamique du flux magnétique, assurant ainsi une efficacité et une précision de l’observateur sur toute la plage de vitesse de rotation du rotor de la machine électrique. La présente invention concerne en outre un procédé et un système de contrôle d’une machine électrique, qui met en œuvre le procédé de détermination du couple de la machine électrique. L’invention concerne un procédé de détermination du couple d’une machine électrique, ladite machine électrique comprenant un rotor, un stator, ledit stator comprenant des bobinages connectés à une pluralité de phases électriques. Pour ce procédé, on met en œuvre les étapes suivantes : On mesure un courant et une tension dans lesdites phases de ladite machine électrique ; On détermine une vitesse de rotation électrique dudit rotor, notamment en fonction d’une vitesse de rotation mécanique dudit rotor ; On construit un modèle dynamique dudit flux magnétique de ladite machine électrique, ledit modèle dynamique dudit flux magnétique étant une équation différentielle dudit flux magnétique qui relie le flux magnétique au courant et à la tension desdites phases de la machine électrique, et à la vitesse de rotation électrique dudit rotor ; On détermine ledit flux magnétique en appliquant un filtre de Kalman sans parfum adaptatif audit modèle dynamique dudit flux magnétique, ledit modèle dynamique de flux magnétique étant appliqué aux mesures de courant et de tension et à la vitesse de rotation électrique déterminée, et la matrice de covariance du bruit du filtre de Kalman sans parfum étant dépendante de ladite vitesse de rotation électrique du rotor ; et On détermine le couple de ladite machine électrique aux moyens desdites mesures de courant et au moyen dudit flux magnétique déterminé. Selon un mode de réalisation de l’invention, on détermine ladite vitesse de rotation électrique ω e dudit rotor par une formule du type ω e =Nω avec N le nombre de paires de pôles de la machine électrique et ω la vitesse de rotation mécanique dudit rotor. Conformément à une mise en œuvre, ledit modèle dynamique de flux magnétique s’écrit : avec R la résistance de la machine électrique, ω e la vitesse de rotation électrique, le flux magnétique direct, le flux magnétique en quadrature, la dérivée par rapport au temps du flux magnétique direct, la dérivée par rapport au temps du flux magnétique en quadrature, I d le courant direct, I q le courant en quadrature, V d la tension directe, V q la tension en quadrature. Avantageusement, on détermine ledit couple de la machine électrique au moyen de la formule : avec le flux magnétique direct, le flux magnétique en quadrature, I d le courant direct, I q le courant en quadrature et N le nombre de paires de pôles de ladite machine électrique. Selon un aspect, on applique ledit filtre de Kalman sans parfum adaptatif à l’équation d’état avec correspondant à l’état non mesurable, avec le flux magnétique direct, le flux magnétique en quadrature, correspondant à l’entrée de ladite équation d’état, avec R la résistance de la machine électrique, I d le courant direct, I q le courant en quadrature, V d la tension directe, V q la tension en quadrature, et avec ω e la vitesse de rotation électrique. Avantageusement, une fonction qui relie les courants et le flux magnétique est une fonction non linéaire obtenue par une méthode des éléments finis. Selon un option de réalisation, on détermine ledit flux magnétique au moyen des étapes suivantes : on initialise k=0, le vecteur d’état = et la matrice de covariance ; à chaque instant k différent de 0, on acquiert lesdites mesures, et on en déduit le vecteur de sortie y(k) ; et à chaque instant k différent de 0, on détermine ledit flux magnétique au moyen des équations suivantes : avec K le gain du filtre de Kalman sans parfum adaptatif, P xy une covariance croisée dudit état et desdites mesures, P yy la covariance prédite des mesures, m y la moyenne prédite de la sortie, I d le courant direct, I q le courant en quadrature. Conformément à une mise en œuvre, ladite machine électrique est une machine électrique synchrone à pôles saillants, de préférence une machine synchro-réluctante assistée d’aimants permanents. En outre, l’invention concerne un procédé de contrôle d’une machine électrique, pour lequel on met en œuvre les étapes suivantes : On détermine le couple de ladite machine électrique au moyen du procédé de détermination dudit couple selon l’une des caractéristiques précédentes ; et On contrôle ladite machine électrique au moyen dudit couple déterminé. De plus, l’invention concerne un système de contrôle d’une machine électrique comprenant une machine électrique et un contrôleur pour mettre en œuvre le procédé de contrôle selon l’une des caractéristiques précédentes. D'autres caractéristiques et avantages du procédé selon l'invention, apparaîtront à la lecture de la description ci-après d'exemples non limitatifs de réalisations, en se référant aux figures annexées et décrites ci-après. Liste des figures La illustre les étapes du procédé selon un mode de réalisation de l’invention. La illustre les étapes du procédé selon un deuxième mode de réalisation de l’invention. La illustre les courbes du courant, de la vitesse de rotation, d’angle de courant et de la tension en fonction du temps, pour un premier exemple d’application du procédé selon un mode de réalisation de l’invention. La illustre des courbes du couple en fonction du temps pour l’exemple de la , une première courbe correspondant au couple mesuré, et la deuxième courbe correspondant au couple déterminé par le procédé selon un mode de réalisation de l’invention. La illustre les courbes de courant, de la vitesse de rotation, de déphasage de courant, et de la tension en fonction du temps, pour un deuxième exemple d’application du procédé selon un mode de réalisation de l’invention. La illustre des courbes du couple en fonction du temps pour l’exemple de la , une première courbe correspondant au couple mesuré, et la deuxième courbe correspondant au couple déterminé par le procédé selon un mode de réalisation de l’invention. La présente invention concerne un procédé de détermination, en temps réel, du couple d’une machine électrique. La machine électrique comprend un rotor et un stator, ce dernier étant équipé de bobinages connectés à plusieurs phases électriques, par exemple à trois phases électriques pour générer un champ magnétique permettant la rotation du rotor. Selon un aspect de l’invention, la machine électrique peut être une machine électrique synchrone à pôles saillants. En effet, le procédé est particulièrement adapté à ce type de machine, d’une part, car le modèle dynamique de la machine électrique est bien représentatif de ce type de machine électrique, et car la détermination du couple permet le contrôle d’une telle machine électrique, notamment par un contrôle direct du couple. De préférence, la machine électrique peut être une machine synchro-réluctante assistée d’aimants permanents. En effet, ce type de machine électrique nécessite une information précise de couple, notamment pour leur contrôle, par exemple au moyen d’une stratégie de contrôle de type MTPA et MTPV, pour leur diagnostic et pour leur surveillance. La décrit, schématiquement et de manière non limitative, les étapes du procédé selon un mode de réalisation de l’invention. Le procédé de détermination du flux magnétique comprend les étapes suivantes : Mesures des courants et des tensions (MES) Détermination de la vitesse de rotation électrique (VIT) Construction du modèle dynamique de flux magnétique (MOD) Application du filtre de Kalman sans parfum adaptatif (KAL) Détermination du couple de la machine électrique (DET) Les étapes 1) à 3) sont indépendantes, et peuvent être réalisées dans cet ordre, dans un ordre différent ou simultanément. L’étape 4) permet de déterminer le flux magnétique de la machine électrique, et l’étape 5) permet de détermine le couple Te de la machine électrique. Les étapes du procédé peuvent être mises en œuvre par des moyens informatiques, notamment au moyen d’un contrôleur de la machine électrique. Ces étapes seront détaillées dans la suite de la description. L’invention concerne également un procédé de contrôle d’une machine électrique. Le procédé de contrôle met en œuvre le procédé de détermination du couple selon l’invention et une étape supplémentaire de contrôle de la machine électrique en fonction du couple déterminé. La décrit, schématiquement et de manière non limitative, les étapes du procédé selon un deuxième mode de réalisation de l’invention. Dans ce cas, le procédé de contrôle de la machine électrique comprend les étapes suivantes : Mesures des courants et des tensions (MES) Détermination de la vitesse de rotation électrique (VIT) Construction du modèle dynamique de flux magnétique (MOD) Application du filtre de Kalman sans parfum adaptatif (KAL) Détermination du couple de la machine électrique (DET) Contrôle de la machine électrique (CTRL) Les étapes 1) à 3) sont indépendantes, et peuvent être réalisées dans cet ordre, dans un ordre différent ou simultanément. L’étape 4) permet de déterminer le flux magnétique de la machine électrique, et l’étape 5) permet de détermine le couple Te de la machine électrique. Les étapes du procédé peuvent être mises en œuvre par des moyens informatiques, notamment au moyen d’un contrôleur de la machine électrique. Ces étapes seront détaillées dans la suite de la description. Dans la suite de la description, les dérivées par rapport au temps sont indiquées par un point. Les notations indexées par la mention (directe) ou (quadrature ou « quadratique ») signifient que les grandeurs sont exprimées dans le repère de Park (repère tournant lié au rotor). En outre, les valeurs à l’état initial sont indiquées avec un 0 (t ou k = 0). 1) Mesure des courants et des tensions Lors de cette étape, on mesure les courants et les tensions dans les phases de la machine électrique. Ces mesures peuvent être réalisées par des capteurs de tension et de courant instrumentant la machine électrique, de manière à mesurer les courants et les tensions dans les phases de la machine électrique. Avantageusement, les mesures dans les phases peuvent être transformées en mesures dans le repère de Park (repère tournant lié au rotor), de manière à déterminer le courant direct, le courant en quadrature, la tension directe, et la tension en quadrature. 2) Détermination de la vitesse de rotation électrique du rotor Lors de cette étape, il s’agit de déterminer la vitesse de rotation électrique du rotor. Selon un mode de réalisation de l’invention, la vitesse de rotation électrique ω e du rotor peut être déterminée à partir de la vitesse de rotation mécanique du rotor, au moyen de la formule : , avec N le nombre de paires de pôles de la machine électrique et ω la vitesse de rotation mécanique dudit rotor. Conformément à une mise en œuvre de ce mode de réalisation, la vitesse de rotation mécanique du rotor peut être estimée, par toute méthode connue de l’homme du métier. Par exemple, la vitesse de rotation mécanique peut être estimée à partir d’une méthode de type boucle à phase asservie PLL (de l’anglais « phase-locked loop »). En variante, la méthode d’estimation de la vitesse de rotation mécanique peut être conforme à celle décrite dans la demande de brevet FR 2 984 637. Alternativement, la vitesse de rotation mécanique du rotor peut être mesurée au moyen d’un capteur de vitesse placé sur la machine électrique. En variante, la vitesse de rotation électrique peut être déterminée directement. 3) Construction du modèle dynamique du flux magnétique Lors de cette étape, on construit un modèle dynamique du flux magnétique de la machine électrique. Le modèle dynamique du flux magnétique électrique est une équation différentielle qui relie le flux magnétique au courant et la tension des phases de la machine électrique, et à la vitesse de rotation électrique du rotor. Le modèle est dit dynamique car il est fonction de la vitesse de rotation du rotor. Le modèle dynamique du flux magnétique est une représentation d’état de la machine électrique. On rappelle qu’en théorie des systèmes (et en automatique), une représentation d’état permet de modéliser un système dynamique sous une forme matricielle, en utilisant des variables d’état. Cette représentation peut être linéaire ou non, continue ou discrète. La représentation permet de déterminer l’état interne et les sorties du système à n’importe quel instant futur si l’on connaît l’état à l’instant initial et le comportement des variables d’entrée qui influent sur le système. La dynamique du flux magnétique peut être exprimée par l’équation différentielle suivante : Avec R la résistance de la machine électrique (pouvant être mesurée au moyen d’un capteur), ω e la vitesse de rotation électrique, le flux magnétique direct, le flux magnétique en quadrature, la dérivée par rapport au temps du flux magnétique direct, la dérivée par rapport au temps du flux magnétique en quadrature, I d le courant direct, I q le courant en quadrature, V d la tension directe, V q la tension en quadrature. On peut réécrire cette équation sous une forme vectorielle compacte : ou de la manière suivante : avec x(t) l’état non mesurable défini par u(t) est l’entrée définie au moyen des mesures par A est la matrice dépendante de la vitesse de rotation électrique du rotor et définie par Ainsi, l’équation d’état du flux peut décrite par un système linéaire variant dans le temps, Toutefois, l’équation de sortie du flux est non linéaire, étant donné que les courants dépendent du flux. Cette dépendance peut être écrite sous la forme : où f d et f q sont des fonctions, respectivement direct et en quadrature. Selon un mode de réalisation, les fonctions f d et f q peuvent être déterminées par une méthode des éléments finis FEMM (de l’anglais « Finite Element Method Magnetics »). Il en résulte qu’il s’agit d’un problème d’estimation non linéaire, étant donné que l’équation de sortie est non linéaire. 4) Application du filtre de Kalman sans parfum adaptatif Lors de cette étape, on détermine le flux magnétique de la machine électrique. Pour cela, on applique un filtre de Kalman sans parfum (UKF de l’anglais « unscented Kalman filter ») adaptatif pour le modèle construit à l’étape 3), appliqué aux mesures de tensions et de courants obtenues à l’étape 1), et à la vitesse de rotation électrique du rotor obtenue à l’étape 2). L’application du filtre Kalman permet d’obtenir un observateur d’état. Le filtre de Kalman sans parfum est un algorithme de filtrage qui utilise un modèle de système pour estimer l’état caché actuel d’un système, puis corrige l’estimation à l’aide des mesures de capteur disponibles. La philosophie de l’UKF diffère du filtre de Kalman étendu en ce sens qu’elle utilise la transformation sans parfum pour approcher directement la moyenne et la covariance de la distribution cible. Le filtre de Kalman sans parfum peut comporter les étapes de prédiction d’état et de correction des mesures, ces deux étapes étant précédées d’une étape préalable pour le calcul des « points sigma ». Les points sigma sont un ensemble d’échantillons calculés de manière à pouvoir propager de manière exacte les informations de moyenne et de variance dans l’espace d’une fonction non linéaire. Selon l’invention, la matrice de covariance du bruit du filtre de Kalman sans parfum est adaptée en fonction de la vitesse de rotation électrique du rotor. Cette adaptation permet de rendre l’observateur performant (précis et robuste) sur toute la plage de vitesse de rotation de la machine électrique. De manière avantageuse, le modèle d’état peut s’écrire : Pour ce mode de réalisation, on peut noter : est l’estimation de x(k) à partir des mesures du temps k-1. est l’estimation de x(k) à partir des mesures du temps k. est la variance d’erreur à partir des mesures du temps k-1. est la variance d’erreur à partir des mesures du temps k. Conformément à une mise en œuvre de l’invention, on peut déterminer le flux magnétique de la machine électrique au moyen des étapes suivantes : on initialise k=0, le vecteur d’état = et l’état de la matrice de covariance, ; à chaque instant k, on acquiert les mesures u(k) de la machine électrique ; et à chaque instant k, on détermine le flux magnétique au moyen des étapes détaillées ci-dessous. Pour l’estimation d’état donnée et l’estimation de la variance d’erreur donnée , il y a deux étapes dans l’application du filtre de Kalman sans parfum : i) prédiction et ii) correction. On note Q(ω) la variance de bruit du système, qui dépend de la vitesse de rotation du rotor de la machine électrique. A l’instant k, pour l’étape de prédiction, on peut écrire : Après l’étape de prédiction, la distribution de x(k) peut être donnée comme une distribution gaussienne pouvant s’écrire : Avec Les points sigma associés à la moyenne m x et à la matrice P x peuvent être calculés comme suit : Avec n=2, et Si est une racine carrée de Px. Les points sigma se propagent dans le modèle de mesure sous la forme suivante, pour tout i compris entre 0 et 2n : les fonctions f d et f q pouvant être celles déterminées par une méthode des éléments finis FEMM. L’étape suivante peut consister à calculer la moyenne prédite m y , la covariance prédite de la mesure P yy et la covariance croisée de l’état de mesure P xy de la manière suivante : avec et les poids de la moyenne prédite et de la covariance prédite, et la matrice de covariance du bruit de mesure. Ensuite, le gain du filtre de Kalman, l’estimation de l’état et la matrice de covariance au temps k peuvent être exprimés de la forme suivante : L’application de ces étapes permet de déterminer l’état x, et par conséquent, le flux magnétique. 5) Détermination du couple de la machine électrique Lors de cette étape, on détermine le couple de la machine électrique au moyen des mesures de courant de l’étape 1) et du flux magnétique déterminé à l’étape 3). Conformément à une mise en œuvre de l’invention, on peut déterminer le couple T e de la machine électrique par une formule du type : avec le flux magnétique direct, le flux magnétique en quadrature, I d le courant direct, I q le courant en quadrature et N le nombre de paires de pôles de la machine électrique. Alternativement, cette étape peut mettre en œuvre tout modèle du couple qui relie le couple et le flux magnétique. Ce modèle peut dépendre également de la vitesse de rotation du rotor de la machine électrique et/ou des mesures de tension. 6) Contrôle de la machine électrique Cette étape est facultative. L'invention concerne également un procédé de contrôle, en temps réel, d'une machine électrique, notamment une machine électrique synchrone, pour lequel on réalise les étapes suivantes : on détermine le couple de la machine électrique à l'aide du procédé (étapes 1) à 5)) décrit ci-dessus ; et on contrôle le couple de ladite machine synchrone en fonction du couple déterminé. Pour cette étape, on peut utiliser tout moyen classique de contrôle de la machine électrique, qui prend en en compte le couple de la machine électrique. Par exemple, le contrôle de la machine électrique peut être basée sur une méthode de contrôle direct du couple (en anglais « direct torque control method ») performante, et particulièrement adaptée pour les machines électriques synchrones à pôles saillants. Il peut s’agir par exemple d’une stratégie de contrôle de type MTPA ou MTPV (respectivement de l’anglais « Maximum Torque per Ampere » pouvant être traduit par couple maximum par ampère, et « Maximum Torque per Voltage » pouvant être traduit par couple maximum par volt). En outre, l'invention concerne un système de contrôle d'une machine électrique, notamment une machine électrique synchrone, adapté à appliquer le procédé tel que décrit ci-dessus. Un tel système de contrôle de machine électrique peut comprendre des moyens de contrôle de la machine électrique comportant des moyens de détermination du couple de la machine électrique et des moyens de contrôle du couple de la machine électrique. Les moyens de détermination du couple déterminent le couple de la machine électrique à partir des mesures de courants et de tensions. Il s'agit des courants et des tensions de chacune des trois phases de la machine électrique. Les moyens de contrôle du couple appliquent des tensions aux bornes de la machine électrique en fonction du couple afin d'assurer une consigne de couple pour la machine électrique. Avantageusement, le système de contrôle peut être un contrôleur comprenant des moyens informatiques. Ce procédé et ce système de contrôle, peuvent être utilisés pour une machine électrique embarquée à bord d'un véhicule, notamment à bord d'un véhicule automobile électrique ou hybride. Toutefois, le système de commande décrit n'est pas limité à cette application et convient pour toutes les applications des machines électriques, y compris les applications stationnaires. Comme il va de soi, l’invention ne se limite pas aux seules formes de réalisation du procédé, décrits ci-dessus à titre d’exemple, elle embrasse au contraire toutes les variantes de réalisation. Exemples d’application Les caractéristiques et avantages du procédé selon l'invention apparaîtront plus clairement à la lecture des exemples d'application ci-après. Les deux exemples concernent des tests expérimentaux pour une machine électrique synchro-réluctante assistée d’aimants permanents selon la conception décrite dans la demande de brevet FR3084536 (WO2020/020580). Pour ces deux exemples, on compare le couple de la machine électrique mesuré au couple de la machine électrique déterminé par le procédé selon l’invention. La machine électrique utilisé possède les caractéristiques suivantes : nombre de paires de pôles : 4 Résistance : 0,0065 Ohm Nombre de phases : 3 Vitesse nominale de la machine électrique : 9000 tr/min Tension de bus : 350 V Fréquence de découpage : 12000 Hz. Pour le premier exemple, le courant de référence est fixé à 400 A en mode de contrôle de courant (en d’autres termes l’amplitude du courant est imposée, et l’angle de courant est automatiquement calculé pour maximiser le couple), et on fait varier la vitesse de rotation de la machine électrique de 500 tr/min à 9000 tr/min. La illustre pour cet exemple : En haut à gauche, la courbe de l’amplitude du courant I abs en A en fonction du temps T en s, qui correspond au courant imposé, En haut à droite, la courbe de la vitesse de rotation mécanique du rotor ω en tr/min en fonction du T en s, qui est imposée, En bas à gauche, la courbe de l’angle de courant δ i en ° en fonction du temps T en s, cette valeur mesurée étant celle qui permet de maximiser le couple, En bas à droite, la courbe de l’amplitude de tension V abs en V en fonction du temps T en s, qui est une valeur mesurée. Au moyen de ces mesures, on applique le procédé selon l’invention et on détermine le couple de la machine électrique. La illustre pour cet exemple le couple Te en Nm, avec une première courbe correspondant au couple mesuré MES et une deuxième courbe correspondant au couple déterminé par le procédé selon l’invention EST. On remarque que les deux courbes sont quasiment superposées. Par conséquent, le procédé selon l’invention permet de déterminer le couple de la machine électrique de manière précise et robuste sur toute la plage de vitesse de rotation. Pour le deuxième exemple, la vitesse de référence est fixée à 500 tr/min, et on fait varier l’amplitude de courant de la machine électrique de 100 A à 500 A. Pour cet essai, la machine électrique n’est pas en zone de défluxage. La illustre pour cet exemple : En haut à gauche, la courbe de l’amplitude du courant I abs en A en fonction du temps T en s, qui correspond au courant imposé, En haut à droite, la courbe de la vitesse de rotation du rotor mécanique ω en tr/min en fonction du T en s, qui est imposée, En bas à gauche, la courbe de l’angle de courant δ i en ° en fonction du temps T en s, cette valeur mesurée étant celle qui permet de maximiser le couple, En bas à droite, la courbe de l’amplitude de tension V abs en V en fonction du temps T en s, qui est une valeur mesurée. Au moyen de ces mesures, on applique le procédé selon l’invention et on détermine le couple de la machine électrique. La illustre pour cet exemple le couple Te en Nm, avec une première courbe correspondant au couple mesuré MES et une deuxième courbe correspondant au couple déterminé par le procédé selon l’invention EST. On remarque que les deux courbes sont quasiment superposées. Par conséquent, le procédé selon l’invention permet de déterminer le couple de la machine électrique de manière précise et robuste sur toute la plage de variation du courant. Procédé de détermination du couple (Te) d’une machine électrique, ladite machine électrique comprenant un rotor, un stator, ledit stator comprenant des bobinages connectés à une pluralité de phases électriques, caractérisé en ce qu’on met en œuvre les étapes suivantes : On mesure (MES) un courant et une tension dans lesdites phases de ladite machine électrique ; On détermine une vitesse de rotation électrique dudit rotor, notamment en fonction d’une vitesse de rotation mécanique dudit rotor ; On construit un modèle dynamique (MOD) dudit flux magnétique de ladite machine électrique, ledit modèle dynamique (MOD) dudit flux magnétique étant une équation différentielle dudit flux magnétique qui relie le flux magnétique au courant et à la tension desdites phases de la machine électrique, et à la vitesse de rotation électrique dudit rotor ; On détermine ledit flux magnétique en appliquant un filtre de Kalman sans parfum adaptatif (KAL) audit modèle dynamique (MOD) dudit flux magnétique, ledit modèle dynamique de flux magnétique étant appliqué aux mesures (MES) de courant et de tension et à la vitesse de rotation électrique déterminée, et la matrice de covariance du bruit du filtre de Kalman sans parfum étant dépendante de ladite vitesse de rotation électrique du rotor ; et On détermine le couple (Te) de ladite machine électrique aux moyens desdites mesures (MES) de courant et au moyen dudit flux magnétique déterminé. Procédé de détermination du couple d’une machine électrique selon la revendication 1, dans lequel on détermine ladite vitesse de rotation électrique ω e dudit rotor par une formule du type ω e =Nω avec N le nombre de paires de pôles de la machine électrique et ω la vitesse de rotation mécanique dudit rotor. Procédé de détermination du couple d’une machine électrique selon l’une des revendications précédentes, dans lequel ledit modèle dynamique (MOD) de flux magnétique s’écrit : avec R la résistance de la machine électrique, ω e la vitesse de rotation électrique, le flux magnétique direct, le flux magnétique en quadrature, la dérivée par rapport au temps du flux magnétique direct, la dérivée par rapport au temps du flux magnétique en quadrature, I d le courant direct, I q le courant en quadrature, V d la tension directe, V q la tension en quadrature. Procédé de détermination du couple d’une machine électrique selon l’une des revendications précédentes, dans lequel on détermine ledit couple (Te) de la machine électrique au moyen de la formule : avec le flux magnétique direct, le flux magnétique en quadrature, I d le courant direct, I q le courant en quadrature et N le nombre de paires de pôles de ladite machine électrique. Procédé de détermination du couple d’une machine électrique selon l’une des revendications précédentes, dans lequel on applique ledit filtre de Kalman sans parfum adaptatif à l’équation d’état avec correspondant à l’état non mesurable, avec le flux magnétique direct, le flux magnétique en quadrature, correspondant à l’entrée de ladite équation d’état, avec R la résistance de la machine électrique, I d le courant direct, I q le courant en quadrature, V d la tension directe, V q la tension en quadrature, et avec ω e la vitesse de rotation électrique. Procédé de détermination du couple d’une machine électrique selon la revendication 5, dans lequel une fonction qui relie les courants et le flux magnétique est une fonction non linéaire obtenue par une méthode des éléments finis. Procédé de détermination du couple d’une machine électrique selon l’une des revendications précédentes, dans lequel on détermine ledit flux magnétique au moyen des étapes suivantes : on initialise k=0, le vecteur d’état = et la matrice de covariance ; à chaque instant k différent de 0, on acquiert lesdites mesures, et on en déduit le vecteur de sortie y(k) ; et à chaque instant k différent de 0, on détermine ledit flux magnétique au moyen des équations suivantes : avec K le gain du filtre de Kalman sans parfum adaptatif, P xy une covariance croisée dudit état et desdites mesures, P yy la covariance prédite des mesures, m y la moyenne prédite de la sortie, I d le courant direct, I q le courant en quadrature. Procédé de détermination du couple d’une machine électrique selon l’une des revendications précédentes, dans lequel ladite machine électrique est une machine électrique synchrone à pôles saillants, de préférence une machine synchro-réluctante assistée d’aimants permanents. Procédé de contrôle d’une machine électrique, caractérisé en ce qu’on met en œuvre les étapes suivantes : On détermine le couple de ladite machine électrique au moyen du procédé de détermination dudit couple selon l’une des revendications précédentes ; et On contrôle (CTRL) ladite machine électrique au moyen dudit couple déterminé. Système de contrôle d’une machine électrique comprenant une machine électrique et un contrôleur pour mettre en œuvre le procédé de contrôle selon la revendication 9.