Estimation et simulation de position quasi-SMO basées sur la force électromotrice en retour de ligne
Principe d’estimation de la position de détection de passage à zéro par champs électromagnétiques
Dans la technologie de détection des champs électromagnétiques des moteurs à courant continu sans balais, il existe une différence de phase d’environ 30° degrés électriques entre le temps de déclenchement réel de la commutation et le point de passage à zéro des champs électromagnétiques arrière. Lorsque le point de passage à zéro est directement utilisé comme référence de commutation, cette caractéristique de retard inhérente entraîne une augmentation de l’ondulation du couple. De plus, la méthode de reconstruction du point neutre comporte des erreurs de mise en œuvre physique, et il y a un léger décalage entre le point neutre virtuel reconstruit et le point neutre théorique, ce qui introduit une distorsion supplémentaire du signal de détection. Afin d’améliorer la précision de la détection de position, la stratégie de détection différentielle de la contre-CEM de ligne est couramment utilisée dans la pratique de l’ingénierie : par le calcul en temps réel de la différence instantanée entre la contre-CEM des deux enroulements de phase sous tension, l’influence des fluctuations de tension de mode commun peut être efficacement éliminée, ce qui permet d’obtenir une capture plus précise du point de passage à zéro et une estimation de la position du rotor. La relation entre la contre-CEM de ligne et la CEM opposée est illustrée dans la formule, de sorte que la forme d’onde de la contre-CEM de ligne de l’enroulement du stator dans des conditions idéales est illustrée à la figure 4.1.
![]()
Figure 4.1 Formes d’onde de contre-CEM et de CEM opposées dans des conditions idéales
Les points de commutation et les points de passage par le zéro de la figure 4.1 correspondent aux étiquettes de la figure 2.6. Les points de passage à zéro de la ligne BLDCM sont Q1, Q2, Q3, Q4, Q5 et Q6. Comme on peut le voir sur la Figure 4.1, le point de passage par le zéro de la ligne contre-CEM est égal au point de temps de commutation du moteur. Par conséquent, en détectant les six points de passage à zéro de la ligne triphasée back-EMF, le contrôle de commutation du BLDCM peut être réalisé.
Principe de contrôle SMO
Le contrôle en mode glissant (SMC), en tant que méthode de contrôle intelligente non linéaire émergente, présente la caractéristique d’une commutation discontinue des lois de contrôle. L’ajustement adaptatif est réalisé en reconstruisant dynamiquement la structure de contrôle du système : le système est composé de plusieurs sous-structures, chacune correspondant à un mode de fonctionnement spécifique. Le contrôleur commute différents modes de fonctionnement en temps réel en fonction de l’état actuel. Ce mécanisme de commutation permet essentiellement l’évolution dynamique de la structure du système en modifiant fréquemment la fonction de commande. Lorsque le système répond à des conditions spécifiques, ce mécanisme de commutation peut reconstruire la trajectoire d’état, forçant les variables du système à osciller de manière à haute fréquence et de faible amplitude le long d’une surface géométrique prédéfinie, appelée surface de glissement, formant un mode de glissement unique. Par rapport aux méthodes de contrôle traditionnelles, cette méthode ajuste dynamiquement la structure de contrôle tout en supprimant les vibrations, en équilibrant la vitesse de réponse du système et la précision en régime permanent. Son avantage réside dans le fait qu’il permet de rompre avec la dépendance à l’égard de modèles de système précis grâce à un mécanisme de commutation non linéaire, de maintenir la robustesse du contrôle dans des conditions de travail complexes, et il est particulièrement adapté aux systèmes fortement non linéaires avec des perturbations externes.
Dans la théorie du contrôle en mode glissant, la surface de commutation sert de limite de l’espace d’état, divisant l’ensemble du domaine d’état en deux régions complémentaires. La figure 4.2 montre intuitivement cette caractéristique de partitionnement dynamique, où la surface de commutation n’est pas seulement une ligne de démarcation mathématique, mais porte également des informations clés sur les transitions d’état du système. Lorsque la trajectoire se déplace à travers cette surface, trois points caractéristiques d’état typiques apparaissent : la trajectoire peut pénétrer directement la surface de commutation, telle que le point A ; il peut également être dispersé des deux côtés avec la surface de commutation comme référence, comme le point B ; ou il peut être attiré par une zone spécifique de la surface de commutation, formant une agrégation stable, telle que le point C.
Dans l’analyse dynamique de SMC, comme le montre la figure 4.2, lorsque la trajectoire passe à travers la surface de commutation, trois modèles de comportement typiques peuvent se produire : Du point de vue de l’analyse du mécanisme de contrôle, le point de pénétration régulier A et le point de départ B ne reflètent que le processus de transition d’état du système ; tandis que le point de terminaison C a une importance de contrôle importante et doit faire l’objet d’une analyse.
Lorsqu’il y a plusieurs terminaisons avec des caractéristiques similaires dans le système, ces points forment une région d’attraction dynamique, et le mouvement de trajectoire dans cette région a des caractéristiques de convergence significatives : tout point de mouvement approchant de cette région sera forcé dans cette région, formant un mode de glissement stable. Définissant la région spéciale formée par plusieurs points de terminaison comme la région du mode glissant, son existence permet au système d’obtenir un contrôle robuste du suivi de trajectoire en reconstruisant dynamiquement la surface de commutation.
Figure 4.2 Schéma caractéristique de trois points sur la surface de commutation Figure 4.3 Deux étages de mouvement du système SMC
Dans le cadre de la théorie du contrôle en mode glissant, les caractéristiques dynamiques du système doivent répondre à de multiples contraintes pour obtenir l’effet de contrôle attendu. Tout d’abord, le mode glissant doit exister, de sorte que la trajectoire de l’état du système soit strictement contrainte dans la structure géométrique prédéfinie de la surface du mode glissant. Deuxièmement, le système doit répondre à la condition d’accessibilité pour s’assurer que toutes les trajectoires de mouvement convergent vers la surface de commutation préconçue dans un temps limité, ce qui est la contrainte fondamentale pour atteindre l’objectif de contrôle. Troisièmement, en construisant la condition de stabilité de Lyapunov, la stabilité asymptotique du mouvement en mode glissant est garantie, et la contrainte stable sur le processus dynamique du système est obtenue par la caractéristique définie négative de la fonction d’énergie. Quatrièmement, grâce à la contrainte globale des conditions ci-dessus, le système peut réaliser une double optimisation de la qualité de la réponse dynamique et de la précision du contrôle, formant une caractéristique de contrôle qui combine robustesse et rapidité.
Dans le contrôle en mode glissant, la sélection de la loi d’atteinte a un impact important sur les caractéristiques de réponse dynamique du système. En fonction des différences dans les caractéristiques dynamiques de l’objet de contrôle réel, diverses formes de lois de portée sont généralement utilisées pour optimiser l’effet de contrôle. Par exemple, la loi d’atteinte exponentielle réalise une convergence rapide grâce aux caractéristiques de décroissance exponentielle, mais peut être accompagnée de bavardage ; tandis que la loi d’atteinte de puissance permet un contrôle équilibré de la vitesse de convergence et de la précision à l’état stationnaire en ajustant l’indice de paramètre ; La loi d’atteinte adaptative améliore encore la capacité d’adaptation robuste du système aux perturbations variables dans le temps grâce à l’ajustement des paramètres en temps réel. Cette stratégie de conception différenciée permet au contrôleur en mode glissant de répondre de manière flexible aux exigences dynamiques des différents systèmes physiques, tout en optimisant les indicateurs de performance globaux et en garantissant la précision du contrôle.
Le contrôle du mode glissant construit une surface de commutation avec les caractéristiques dynamiques souhaitées, conduisant la trajectoire de l’état du système à s’approcher progressivement et finalement à se stabiliser sur cette surface géométrique à partir de la position initiale, formant une caractéristique unique de mouvement du mode coulissant. Ce mécanisme de contrôle permet d’atteindre la stabilité globale de l’état du système en un temps limité, et son processus de contrôle peut être divisé en deux étapes : lorsque la trajectoire est dans le segment AB, le système est dans l’étape du mode d’approche et les variables d’état s’approchent progressivement de la surface de glissement prédéfinie le long d’une trajectoire non linéaire ; lorsque la trajectoire atteint le segment BC, le système entre dans l’étape du mode glissant, et le mouvement d’état est strictement contraint dans la structure géométrique de la surface du mode glissant, et son processus de changement de caractéristique dynamique est illustré à la figure 4.2.
SMO, en tant que système d’estimation dynamique, construit un modèle d’estimation des variables d’état à l’aide de variables mesurables externes telles que les signaux d’entrée et de sortie, résolvant efficacement le problème selon lequel les capteurs traditionnels ne peuvent pas mesurer directement des paramètres tels que la vitesse du moteur, la position du rotor et la force électromotrice induite. Dans le domaine du contrôle moteur, cette technologie est devenue un outil important pour améliorer les performances de réponse dynamique du système.
Principe de fonctionnement du quasi-SMO
Dans les systèmes de commande de moteur, en raison des caractéristiques physiques des enroulements du stator et de l’effet de décalage inhérent aux dispositifs de commutation électronique de puissance, des phénomènes évidents d’inertie et de retard se produiront en fonctionnement réel. Ces caractéristiques dynamiques non idéales conduisent à un modèle de mouvement complexe de la trajectoire d’état du système près de la surface de glissement : dans l’étape du mode d’approche, le système ne peut pas atteindre complètement un mouvement de glissement idéal, mais forme une trajectoire en dents de scie près de la surface de commutation, accompagnée d’oscillations à cycle limite. Bien que ce phénomène de broutage ne puisse pas être complètement éliminé, l’amplitude de broutage peut toujours être efficacement supprimée en optimisant la conception de la surface en mode glissant et la configuration des paramètres de la loi de contrôle. Cependant, un broutage continu réduira non seulement considérablement la vitesse de réponse dynamique et la précision en régime permanent du système, mais entraînera également une perte d’énergie supplémentaire, voire une instabilité du système.
Afin de résoudre ce problème, les chercheurs ont proposé dans les années 1980 la théorie du contrôle en mode quasi-glissant, qui fournit une nouvelle solution pour la stabilité du système et l’optimisation des performances dynamiques en introduisant une fonction de saturation au lieu de la fonction de signe traditionnelle. Par rapport au problème de discontinuité de la fonction de signe, la fonction de saturation évite efficacement le broutage à haute fréquence dans le processus de commutation du système grâce à des caractéristiques non linéaires lisses, tout en maintenant des capacités d’estimation précise de l’état du système. Cette amélioration atténue non seulement les fluctuations drastiques des paramètres dans l’étape du mode d’atteinte de l’observateur traditionnel en mode glissant, mais renforce également la robustesse du système de contrôle, fournissant une base théorique et une voie de mise en œuvre technique pour le contrôle précis de processus industriels complexes.
L’expression de la fonction de saturation couramment utilisée est
Où : k est la constante de gain ; ±ε est la limite de la couche limite ; |ε| est une constante, représentant l’épaisseur de la couche limite.
Le schéma de coordonnées correspondant est illustré à la figure 4.4. Dans la stratégie SMC, le mécanisme de contrôle de la couche limite permet d’optimiser les performances de manière dynamique grâce à l’ajustement régional. Lorsque l’état du système est situé à l’extérieur de la couche limite, une méthode de contrôle par relais est utilisée pour faire se déplacer le système vers la surface de glissement ; Lorsque l’état entre dans la couche limite, passez en mode de contrôle linéaire pour supprimer les vibrations à haute fréquence. Cette méthode de contrôle garantit les performances dynamiques de l’étage d’approche et la convergence en douceur de la trajectoire dans l’étage stationnaire en ajustant les paramètres de gain de contrôle linéaire et d’épaisseur de la couche limite.
Figure 4.4 Fonction de saturation
Conception de Quasi-SMO
L’équation de tension de ligne triphasée du moteur est la suivante :
Dans la formule : uab est la tension de ligne ; iab est le courant de ligne.
La relation triphasée entre la ligne et la force électromotrice du moteur est la suivante :
La formule s’écrit sous forme matricielle comme suit :
Définissez iab dans la formule en tant que variable d’état, uab en tant qu’entrée système et eab en tant que sortie système. L’équation d’état du moteur peut être dérivée de la formule comme suit :
L’expression de la fonction de saturation f(x) utilisée dans la partie SMO de cet article est la suivante :
Simulink Simulation Modeling Analysis
Le modèle de simulation global du BLDCM quasi-SMO sans capteur est construit dans Simulink dans MATLAB, comme le montre la Figure 4.5. Parmi eux, la valeur initiale de la source contrôlée en courant continu est de 5 V et le tube de commutation de l’onduleur triphasé adopte l’IGBT, comme le montre la figure 3.14. Comme dans le chapitre précédent, le module de corps de moteur fourni avec Simulink, à savoir le module de moteur synchrone à aimants permanents (PMSM), est utilisé, et sa force électromotrice arrière configurée est remplacée par une onde trapézoïdale, qui est BLDCM. Les paramètres de son module sont définis comme suit : résistance du stator 2,875 Ω, inductance du stator 0,0085H, moment d’inertie 0,0008 kg.m2, coefficient d’amortissement 0N.m.s, nombre de paires de pôles du rotor 2, frottement statique 0 N.m. Définissez les variables du moteur dans l’interface Modélisation/Paramètres du modèle/Propriétés du modèle/Callbacks/InitFcn, comme illustré à la Figure 4.6.
![]()
![]()
Figure 4.5 Modèle de simulation global du BLDCM sans capteur quasi-SMO
Figure 4.6 Propriétés du modèle
Le module de sortie et d’acquisition du signal est illustré à la Figure 4.7.
Figure 4.7 Module de sortie et d’acquisition de signal
Le module d’observation des champs électromagnétiques arrière de l’OMM est illustré à la figure 4.8. La limite supérieure du réglage des paramètres du module de limite (Saturation) est 1 et la limite inférieure est -1. Le module de filtre passe-bas de SMO est illustré à la Figure 4.9. L’observation de la contre-CEM de la ligne SMO et les courbes de forme d’onde de la réplique réelle des phases B et C coïncident presque, ce qui indique que la conception de la SMO est correcte, comme le montre la Figure 4.10.
Figure 4.8 Module d’observation de la FOM arrière de l’OMS
Figure 4.9 Module de filtre passe-bas dans SMO
Figure 4.10 Observation de la contre-CEM de la ligne SMO et rétro-CEM de la ligne réelle des phases BC
Dans le système de commande de moteur à courant continu sans balais sans capteur, la surveillance en temps réel de la vitesse du moteur est la condition essentielle du réglage de la boucle de vitesse. Selon la relation entre la vitesse du moteur et la force électromotrice maximale de la ligne, la vitesse du rotor peut être obtenue :
Dans la formule : n est la vitesse du rotor du moteur, Em est la valeur maximale de la ligne contre-CEM et Ke est le coefficient de la contre-CEM.
La ligne back-EMF change périodiquement dans une forme d’onde trapézoïdale, et sa valeur maximale subit un saut tous les 60 degrés électriques, et chaque phase a une valeur maximale pendant la conduction. Utilisez le module ABS de Simulink pour convertir le signal contre-EMF d’origine en une valeur absolue afin d’éliminer les interférences de signe causées par la différence de phase ; Ensuite, le comparateur produit la valeur maximale ; enfin, la vitesse est calculée en divisant la valeur maximale par le double du coefficient de rétro-CEM. Le module de mesure de vitesse est illustré à la Figure 4.11. Afin de réduire les interférences de signal haute fréquence dans la mesure de la vitesse, le filtre passe-bas illustré à la Figure 4.11 est utilisé.
Figure 4.11 Module de mesure de la vitesse Figure 4.12 Module de filtre passe-bas
Comme dans le chapitre précédent, le module de régulation de vitesse utilise la méthode de séparation intégrale et la limitation intégrale pour limiter l’amplitude cumulative du terme intégral et réduire l’impact du dépassement intégral. Comme le montre la figure 4.13. Le régulateur de courant est illustré à la Figure 4.14, où le coefficient proportionnel est fixé à 5, le coefficient intégral est 1 et la limite supérieure de sortie est limitée à 1 et la limite inférieure est -1.
Figure 4.13 Module régulateur de vitesse
Figure 4.14 Module régulateur de courant
Le module logique de commutation est illustré à la Figure 4.15. Étant donné que la bobine du rotor BLDCM est connectée en Y, la ligne de rotor arrière-CEM mène la CEM opposée de 30°, et parce que le point de passage à zéro de la CEM opposée précède le point de changement de pas du signal Hall de 30°, elle mène d’un total de 60°, il est donc nécessaire de retarder la commutation d’un battement. Cette fonction est écrite dans le module Function de MATLAB, et le code est illustré à la Figure 4.16, qui réalise le retard du battement de commutation de un cinq quatre six deux trois à trois un cinq quatre six deux, comme le montre la Figure 4.16.
Fig. 4.15 Module logique de commutation
Fig. 4.16 Code de la fonction MATLAB Fig. 4.17 Diagramme du retard de la fonction MATLAB
Comme dans le chapitre précédent, la commande HPWM_LON illustrée à la Fig. 3.2 (c) est adoptée, et le module de commande PWM est illustré à la Fig. 4.18. Le module de séquence répétitive est utilisé pour générer une onde porteuse triangulaire de 10 kHz, avec des valeurs de temps définies sur [0 1/20000 2/20000] et des valeurs de sortie définies sur [0 1 0].
Fig. 4.18 Module de commande PWM
Analyse des résultats de simulation
Le couple électromagnétique réel et les formes d’onde de couple électromagnétique souhaitées sont illustrés à la Fig. 4.19. Selon les résultats de la simulation, lorsque le couple électromagnétique souhaité change brusquement, le couple électromagnétique réel peut également suivre le changement. Cependant, en raison de l’utilisation de la commutation en six étapes, il y a une pulsation évidente.
Fig. 4.19 Couple électromagnétique réel et formes d’onde de couple électromagnétique souhaitées
La forme d’onde de la vitesse globale est illustrée à la Fig. 4.20. Selon les résultats de la simulation, lorsque la vitesse souhaitée augmente soudainement (0,4 s) ou diminue (0,7 s), la vitesse de sortie peut être suivie efficacement en peu de temps, et elle peut également maintenir rapidement la vitesse lorsque la charge augmente soudainement (0,15 s et 0,5 s) ou diminue (0,85 s). Comme le montre la forme d’onde de vitesse locale de la Fig. 4.21, l’erreur de vitesse est concentrée entre -2,5 tr/min et +5 tr/min. La forme d’onde de la vitesse de démarrage est illustrée à la Fig. 4.21. Un démarrage en trois étapes est adopté : pré-positionnement du rotor de 0 à 0,01 s, positionnement secondaire de 0,01 s à 0,02 s, passage à l’accélération haute tension et début de l’utilisation du contrôle PWM à partir de 0,02 s, et la vitesse oscille légèrement autour de la vitesse donnée à partir de 0,03 s, complétant le démarrage du moteur.
Si la vitesse souhaitée est modifiée à 1600 tr/min, la forme d’onde de vitesse est illustrée à la Fig. 4.23. Selon les résultats de la simulation, la vitesse est stable, mais la fluctuation de vitesse augmente, et l’erreur de vitesse est concentrée entre -10 tr/min et +15 tr/min.
Fig. 4.20 Forme d’onde de vitesse globale
Fig. 4.21 Forme d’onde de vitesse locale Fig. 4.22 Forme d’onde de démarrage de vitesse
Fig. 4.23 Forme d’onde de vitesse avec vitesse initiale donnée de 1600 tr/min
Le courant de phase A du stator BLDCM est illustré à la Fig. 4.24. Les formes d’onde de courant triphasé sont similaires, les trois phases ABC étant en retard de 120° en séquence, et oscillant toutes autour de 0. L’amplitude de changement de courant est positivement corrélée avec le couple de charge.
Fig. 4.24 Forme d’onde du courant de phase A du stator
Résumé du chapitre
Ce chapitre présente d’abord le principe de détection du passage à zéro de la force électromotrice arrière pour l’estimation de la position, puis introduit le principe de contrôle SMO, puis modélise et analyse le principe de contrôle quasi-SMO à l’aide d’une fonction par morceaux pour remplacer la fonction de signe dans Simulink de MATLAB, en analysant principalement le couple de charge, la vitesse du rotor et le courant de phase du stator pour vérifier la précision de la théorie.