CA3159341C - Procede de commande d'un actionneur de systeme mecanique a frottements imbriques - Google Patents
Procede de commande d'un actionneur de systeme mecanique a frottements imbriques Download PDFInfo
- Publication number
- CA3159341C CA3159341C CA3159341A CA3159341A CA3159341C CA 3159341 C CA3159341 C CA 3159341C CA 3159341 A CA3159341 A CA 3159341A CA 3159341 A CA3159341 A CA 3159341A CA 3159341 C CA3159341 C CA 3159341C
- Authority
- CA
- Canada
- Prior art keywords
- friction
- actuator
- algorithm
- mechanical system
- controlling
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Program-control systems
- G05B19/02—Program-control systems electric
- G05B19/18—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of program data in numerical form
- G05B19/404—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of program data in numerical form characterised by control arrangements for compensation, e.g. for backlash, overshoot, tool offset, tool wear, temperature, machine construction errors, load, inertia
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Program-control systems
- G05B19/02—Program-control systems electric
- G05B19/418—Total factory control, i.e. centrally controlling a plurality of machines, e.g. direct or distributed numerical control [DNC], flexible manufacturing systems [FMS], integrated manufacturing systems [IMS] or computer integrated manufacturing [CIM]
- G05B19/41885—Total factory control, i.e. centrally controlling a plurality of machines, e.g. direct or distributed numerical control [DNC], flexible manufacturing systems [FMS], integrated manufacturing systems [IMS] or computer integrated manufacturing [CIM] characterised by modeling, simulation of the manufacturing system
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/30—Nc systems
- G05B2219/32—Operator till task planning
- G05B2219/32091—Algorithm, genetic algorithm, evolution strategy
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/30—Nc systems
- G05B2219/41—Servomotor, servo controller till figures
- G05B2219/41154—Friction, compensation for friction
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/30—Nc systems
- G05B2219/41—Servomotor, servo controller till figures
- G05B2219/41161—Adaptive friction compensation
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Human Computer Interaction (AREA)
- General Engineering & Computer Science (AREA)
- Quality & Reliability (AREA)
- Feedback Control In General (AREA)
- Control Of Position Or Direction (AREA)
- Control Of Electric Motors In General (AREA)
Abstract
Procédé de commande d'un actionneur électrique d'un système mécanique comprenant une pluralité de zones de contact imbriquées, le procédé comprenant les étapes de : - acquérir des données sur le système mécanique qui comprend un nombre de zones de contact imbriquées, - élaborer une modélisation du système, à partir desdites données et d'un nombre de modèles de Lugre mis en parallèle égal au nombre de zones de contact imbriquées, et déterminer des paramètres de la modélisation et une structure de compensation des frottements modélisés, - inclure la structure de compensation dans une loi de commande de l'actionneur, - commander l'actionneur au moyen de la loi de commande.
Description
PROCEDE DE COMMANDE D'UN ACTIONNEUR DE SYSTEME MECANIQUE
A FROTTEMENTS IMBRIQUES
La présente invention concerne le domaine de la commande des actionneurs électriques dans des systèmes mécaniques comprenant un élément mobile par rapport à un élément fixe.
ARRIERE PLAN DE L'INVENTION
De nombreux systèmes mécaniques incorporent un ou plusieurs actionneurs électriques, comme des actionneurs électromécaniques, qui sont reliés à une source d'alimentation en énergie électrique et sont pilotés pour déplacer l'élément mobile du système mécanique de manière très précise par rapport à l'élément fixe. Le pilotage de chaque actionneur est réalisé à partir d'une loi de commande mettant par exemple en relation des paramètres d'alimentation de l'actionneur avec des positions ou des vitesses de l'élément mobile.
Lorsque deux pièces ayant des surfaces en contact doivent frotter l'une contre l'autre pour que l'élément mobile puisse se déplacer, il est nécessaire de piloter l'actionneur pour que celui-ci produise un effort suffisant pour vaincre ces frottements de manière à ce que le déplacement se fasse sans à-coup. A cette fin, la loi de commande comprend une structure de compensation permettant d'alimenter l'actionneur pour qu'il vainque les frottements et déplace l'élément mobile de manière régulière.
Il existe plusieurs modèles permettant une modélisation des frottements. Ces modèles sont basés sur le comportement d'aspérités présentes sur les surfaces en contact et ces aspérités sont modélisées sous la forme de lames élastiques de telle manière que les lames élastiques d'une des surfaces en contact interagissent avec des lames rigides de l'autre des surfaces en contact.
A FROTTEMENTS IMBRIQUES
La présente invention concerne le domaine de la commande des actionneurs électriques dans des systèmes mécaniques comprenant un élément mobile par rapport à un élément fixe.
ARRIERE PLAN DE L'INVENTION
De nombreux systèmes mécaniques incorporent un ou plusieurs actionneurs électriques, comme des actionneurs électromécaniques, qui sont reliés à une source d'alimentation en énergie électrique et sont pilotés pour déplacer l'élément mobile du système mécanique de manière très précise par rapport à l'élément fixe. Le pilotage de chaque actionneur est réalisé à partir d'une loi de commande mettant par exemple en relation des paramètres d'alimentation de l'actionneur avec des positions ou des vitesses de l'élément mobile.
Lorsque deux pièces ayant des surfaces en contact doivent frotter l'une contre l'autre pour que l'élément mobile puisse se déplacer, il est nécessaire de piloter l'actionneur pour que celui-ci produise un effort suffisant pour vaincre ces frottements de manière à ce que le déplacement se fasse sans à-coup. A cette fin, la loi de commande comprend une structure de compensation permettant d'alimenter l'actionneur pour qu'il vainque les frottements et déplace l'élément mobile de manière régulière.
Il existe plusieurs modèles permettant une modélisation des frottements. Ces modèles sont basés sur le comportement d'aspérités présentes sur les surfaces en contact et ces aspérités sont modélisées sous la forme de lames élastiques de telle manière que les lames élastiques d'une des surfaces en contact interagissent avec des lames rigides de l'autre des surfaces en contact.
2 Le modèle de LuGre est utilisable pour modéliser les frottements lorsque deux pièces seulement sont en contact et permet d'en tirer une structure de compensation stable des frottements comme indiqué dans le document De Wit et al., "A new model for control of systems with friction", pages 419 à 425, n 3, volume 40, IEEE Transactions on automatic control, 1995, IEEE. On a représenté à la figure 1 une pièce P reposant sur une surface S et se déplaçant à une vitesse v avec un frottement F s'opposant au déplacement de la pièce P, ainsi que sa représentation conformément au modèle de LuGre. Le modèle de LuGre s'écrit alors sous la forme :
VO ¨ 0z(t) + 1(t) + 2w(t) 1(0(01 i(t) = WM
g(COM)z(t) g(w) = rc + (rs ¨ rc)e s) w (70 où :
- Pf(t) est le couple de frottement total, ¨ Go est la raideur des lames élastiques, - Gi est le coefficient d'amortissement du mouvement des lames élastiques, - o2 est le coefficient de frottement visqueux, - rc est le couple de Coulomb, ¨ rs est le couple de frottement de Stribeck ou statique, - cos est la vitesse angulaire de Stribeck, - co est la vitesse angulaire relative.
Il existe d'autres modèles de frottements plus précis que le modèle de LuGre mais plus complexes dans leur structure (et donc dans leur mise en uvre) tels que le modèle de Leuven (décrit dans le document Lampaert et
VO ¨ 0z(t) + 1(t) + 2w(t) 1(0(01 i(t) = WM
g(COM)z(t) g(w) = rc + (rs ¨ rc)e s) w (70 où :
- Pf(t) est le couple de frottement total, ¨ Go est la raideur des lames élastiques, - Gi est le coefficient d'amortissement du mouvement des lames élastiques, - o2 est le coefficient de frottement visqueux, - rc est le couple de Coulomb, ¨ rs est le couple de frottement de Stribeck ou statique, - cos est la vitesse angulaire de Stribeck, - co est la vitesse angulaire relative.
Il existe d'autres modèles de frottements plus précis que le modèle de LuGre mais plus complexes dans leur structure (et donc dans leur mise en uvre) tels que le modèle de Leuven (décrit dans le document Lampaert et
3 al., "Modification of the Leuven integrated friction model structure", pages 683 à 687, n 4, volume 47, IEEE
transactions on Automatic Control, 2002, IEEE) ou le modèle GMS ou Generalized Maxwell Slip (Lampaert et al., "A generalized Maxwell-slip friction model appropriate for control purposes", pages 1170 à 1177, volume 4, 2003 IEEE International Workshop on Workload Characterization, 2003, IEEE). De par leur complexité, ces modèles permettent aussi de modéliser des architectures mécaniques plus complexes. Ceci est d'autant plus vrai lorsque plusieurs surfaces de contact sont imbriquées les unes dans les autres comme dans le cas d'un arbre pivotant dans une première douille pivotant elle-même dans une deuxième douille qui pivote elle aussi dans un palier avec éventuellement un joint d'étanchéité frottant sur l'une ou l'autres des pièces en mouvement : il n'est pas possible de déterminer la contribution des différentes surfaces de contact dans le frottement global existant entre l'arbre et le palier.
Cependant, ces modèles ne permettent pas de compenser de manière stable les frottements ainsi modélisés.
OBJET DE L'INVENTION
L'invention a notamment pour but de proposer un moyen pour améliorer la précision et la finesse de la commande d'un actionneur d'un système mécanique ayant plusieurs surfaces de contact.
RESUME DE L'INVENTION
A cet effet, on prévoit, selon l'invention un procédé de commande d'un actionneur électrique d'un système mécanique comprenant une pluralité de zones de contact imbriquées, le procédé comprenant les étapes de :
- acquérir des données sur le système mécanique qui comprend un nombre de zones de contact imbriquées, - élaborer une modélisation du système, à partir desdites données et d'un nombre de modèles de Lugre égal au nombre
transactions on Automatic Control, 2002, IEEE) ou le modèle GMS ou Generalized Maxwell Slip (Lampaert et al., "A generalized Maxwell-slip friction model appropriate for control purposes", pages 1170 à 1177, volume 4, 2003 IEEE International Workshop on Workload Characterization, 2003, IEEE). De par leur complexité, ces modèles permettent aussi de modéliser des architectures mécaniques plus complexes. Ceci est d'autant plus vrai lorsque plusieurs surfaces de contact sont imbriquées les unes dans les autres comme dans le cas d'un arbre pivotant dans une première douille pivotant elle-même dans une deuxième douille qui pivote elle aussi dans un palier avec éventuellement un joint d'étanchéité frottant sur l'une ou l'autres des pièces en mouvement : il n'est pas possible de déterminer la contribution des différentes surfaces de contact dans le frottement global existant entre l'arbre et le palier.
Cependant, ces modèles ne permettent pas de compenser de manière stable les frottements ainsi modélisés.
OBJET DE L'INVENTION
L'invention a notamment pour but de proposer un moyen pour améliorer la précision et la finesse de la commande d'un actionneur d'un système mécanique ayant plusieurs surfaces de contact.
RESUME DE L'INVENTION
A cet effet, on prévoit, selon l'invention un procédé de commande d'un actionneur électrique d'un système mécanique comprenant une pluralité de zones de contact imbriquées, le procédé comprenant les étapes de :
- acquérir des données sur le système mécanique qui comprend un nombre de zones de contact imbriquées, - élaborer une modélisation du système, à partir desdites données et d'un nombre de modèles de Lugre égal au nombre
4 PCT/EP2020/079953 de zones de contact imbriquées, et déterminer des paramètres de la modélisation et une structure de compensation des frottements modélisés, - inclure la structure de compensation dans une loi de commande de l'actionneur, - commander l'actionneur au moyen de la loi de commande.
Le procédé de l'invention permet de modéliser tout type de frottement multiple et de compenser efficacement les frottements ainsi modélisés pour obtenir une loi de commande particulièrement précise et stabilisante.
D'autres caractéristiques et avantages de l'invention ressortiront à la lecture de la description qui suit d'un mode de réalisation particulier et non limitatif de l'invention.
BREVE DESCRIPTION DES DESSINS
Il sera fait référence aux dessins annexés, parmi lesquels :
La figure 1 est une représentation d'une situation de frottement et sa modélisation par un modèle de LuGre classique ;
La figure 2 est une vue schématique représentant un système mécanique en relation avec lequel peut être mis en uvre le procédé selon l'invention ;
La figure 3 est une vue schématique du dispositif de commande ;
La figure 4 est une représentation schématique du compensateur du dispositif de commande ;
La figure 5 est une représentation d'une boucle d'asservissement utilisée pour la mise en uvre du procédé de l'invention.
DESCRIPTION DETAILLEE DE L'INVENTION
En référence à la figure 2, le système mécanique utilisé
pour illustrer le fonctionnement de l'invention comprend quatre pièces en frottement relatif comme par exemple un rotor d'un moteur électrique tournant dans un stator et lié en rotation à un arbre mené qui tourne dans un palier de guidage en rotation. Dans un tel agencement, on trouve
Le procédé de l'invention permet de modéliser tout type de frottement multiple et de compenser efficacement les frottements ainsi modélisés pour obtenir une loi de commande particulièrement précise et stabilisante.
D'autres caractéristiques et avantages de l'invention ressortiront à la lecture de la description qui suit d'un mode de réalisation particulier et non limitatif de l'invention.
BREVE DESCRIPTION DES DESSINS
Il sera fait référence aux dessins annexés, parmi lesquels :
La figure 1 est une représentation d'une situation de frottement et sa modélisation par un modèle de LuGre classique ;
La figure 2 est une vue schématique représentant un système mécanique en relation avec lequel peut être mis en uvre le procédé selon l'invention ;
La figure 3 est une vue schématique du dispositif de commande ;
La figure 4 est une représentation schématique du compensateur du dispositif de commande ;
La figure 5 est une représentation d'une boucle d'asservissement utilisée pour la mise en uvre du procédé de l'invention.
DESCRIPTION DETAILLEE DE L'INVENTION
En référence à la figure 2, le système mécanique utilisé
pour illustrer le fonctionnement de l'invention comprend quatre pièces en frottement relatif comme par exemple un rotor d'un moteur électrique tournant dans un stator et lié en rotation à un arbre mené qui tourne dans un palier de guidage en rotation. Dans un tel agencement, on trouve
5 au moins trois zones de frottements qui s'additionnent et qu'on peut représenter par quatre pièces 1, 2, 3, 4 coaxiales les unes aux autres et en rotation les unes par rapport aux autres, les quatre pièces étant en contact les unes avec les autres :
- la pièce 1 a une surface externe Si' en contact avec une surface interne S2 de la pièce 2, - la pièce 2 a une surface externe S2' en contact avec une surface interne S3 de la pièce 3, - la pièce 3 a une surface externe S3' en contact avec la surface interne S4 de la pièce 4.
Le système comprend donc trois zones de frottement :
S1'/S2, 52'/53, 53'/54. On dit que ces zones de frottement sont imbriquées.
Le procédé de l'invention est agencé pour commander le moteur électrique qui constitue l'actionneur électrique du système mécanique. Le procédé de commande comprend les étapes de :
- acquérir des données sur le système mécanique qui comprend un nombre M de zones de contact imbriquées, - élaborer une modélisation, dite M-LuGre, du système, à
partir desdites données et d'un nombre M de modèles de Lugre, mis en parallèle, égal au nombre de zones de contact imbriquées, en utilisant un algorithme d'optimisation stochastique pour déterminer des paramètres de la modélisation et déterminer une structure de compensation des frottements dans les
- la pièce 1 a une surface externe Si' en contact avec une surface interne S2 de la pièce 2, - la pièce 2 a une surface externe S2' en contact avec une surface interne S3 de la pièce 3, - la pièce 3 a une surface externe S3' en contact avec la surface interne S4 de la pièce 4.
Le système comprend donc trois zones de frottement :
S1'/S2, 52'/53, 53'/54. On dit que ces zones de frottement sont imbriquées.
Le procédé de l'invention est agencé pour commander le moteur électrique qui constitue l'actionneur électrique du système mécanique. Le procédé de commande comprend les étapes de :
- acquérir des données sur le système mécanique qui comprend un nombre M de zones de contact imbriquées, - élaborer une modélisation, dite M-LuGre, du système, à
partir desdites données et d'un nombre M de modèles de Lugre, mis en parallèle, égal au nombre de zones de contact imbriquées, en utilisant un algorithme d'optimisation stochastique pour déterminer des paramètres de la modélisation et déterminer une structure de compensation des frottements dans les
6 zones de contact imbriquées, - inclure la structure de compensation dans une loi de commande de l'actionneur, - commander l'actionneur au moyen de la loi de commande.
Le modèle M-LuGre met en parallèle M modèles de LuGre i pour calculer le couple de frottement total Pf(t) comme suit :
(Go,izi(t) + (0) + 62 W (t) i=1 IW(01 (t) = w(t) _______________________________________ i(co (0)z i(t) g .) r + (rs,i r ,i)e cos,l g i(W) OE0,i dans lequel, pour chaque modèle i :
- Go,i est la raideur des lames élastiques représentant les aspérités en frottement, est le coefficient d'amortissement du mouvement des lames, ¨ G2 est le coefficient de frottement visqueux (on notera que, par factorisation, on évite de devoir introduire G2, ) - zi est la déformation moyenne de la lame, - rc,i est le couple de Coulomb, ¨ rs,i est le couple de frottement statique, - cos,i est la vitesse angulaire de Stribeck, - co est la vitesse angulaire relative mesurée du moteur.
La structure de compensation, de type FEEDBACK, est basée sur une estimation du modèle précédent et a la forme suivante :
Le modèle M-LuGre met en parallèle M modèles de LuGre i pour calculer le couple de frottement total Pf(t) comme suit :
(Go,izi(t) + (0) + 62 W (t) i=1 IW(01 (t) = w(t) _______________________________________ i(co (0)z i(t) g .) r + (rs,i r ,i)e cos,l g i(W) OE0,i dans lequel, pour chaque modèle i :
- Go,i est la raideur des lames élastiques représentant les aspérités en frottement, est le coefficient d'amortissement du mouvement des lames, ¨ G2 est le coefficient de frottement visqueux (on notera que, par factorisation, on évite de devoir introduire G2, ) - zi est la déformation moyenne de la lame, - rc,i est le couple de Coulomb, ¨ rs,i est le couple de frottement statique, - cos,i est la vitesse angulaire de Stribeck, - co est la vitesse angulaire relative mesurée du moteur.
La structure de compensation, de type FEEDBACK, est basée sur une estimation du modèle précédent et a la forme suivante :
7 '\ 2 F + (F . ¨ F)e `a) =
g1(o) = C,1 C,i 6 o, le)1 =
z "E.
dt g() E=Oc Or.
M ( ckµ
Ff = + al,i¨ + a20) 0,1 i=1 dt dans laquelle :
- s est l'erreur entre la consigne d'entrée ou signal de référence (position ou vitesse) et le signal de sortie obtenu Or (position ou vitesse selon la consigne entrée) ;
- les variables ks,i sont au nombre de M et sont déterminées pour améliorer la robustesse du compensateur au regard de modifications des conditions extérieures comme la température ou des évolutions du système lui-même, - les autres paramètres et variables ont été mentionnés précédemment en liaison avec le modèle.
L'ajustement du compensateur comprend l'étape d'identifier tous les paramètres de chacun des modèles de LuGre i de la modélisation M-LuGre et l'étape de déterminer les variables ks,i, ks,1,1 assurant un fonctionnement optimal et stabilisant du compensateur.
Les paramètres du modèle M-LuGre sont définis en utilisant un algorithme d'optimisation stochastique, ici un algorithme à évolution différentielle. Cet algorithme est décrit dans le document P. Feyel et al., "LuGre Friction Model Identification and Compensator Tuning
g1(o) = C,1 C,i 6 o, le)1 =
z "E.
dt g() E=Oc Or.
M ( ckµ
Ff = + al,i¨ + a20) 0,1 i=1 dt dans laquelle :
- s est l'erreur entre la consigne d'entrée ou signal de référence (position ou vitesse) et le signal de sortie obtenu Or (position ou vitesse selon la consigne entrée) ;
- les variables ks,i sont au nombre de M et sont déterminées pour améliorer la robustesse du compensateur au regard de modifications des conditions extérieures comme la température ou des évolutions du système lui-même, - les autres paramètres et variables ont été mentionnés précédemment en liaison avec le modèle.
L'ajustement du compensateur comprend l'étape d'identifier tous les paramètres de chacun des modèles de LuGre i de la modélisation M-LuGre et l'étape de déterminer les variables ks,i, ks,1,1 assurant un fonctionnement optimal et stabilisant du compensateur.
Les paramètres du modèle M-LuGre sont définis en utilisant un algorithme d'optimisation stochastique, ici un algorithme à évolution différentielle. Cet algorithme est décrit dans le document P. Feyel et al., "LuGre Friction Model Identification and Compensator Tuning
8 Using a Differential Evolution Algorithm", pages 85 à 91, 2013 IEEE Symposium on Differential Evolution (SDE), 2013, IEEE. De préférence, l'algorithme est utilisé
plusieurs fois successivement pour affiner la détermination des paramètres de la modélisation. La détermination des paramètres de réglage de l'algorithme est connue en elle-même. On pourra utiliser les paramètres de réglage standards de l'algorithme à
évolution différentielle tels que ceux mentionnés dans le document R. Storn et K. Price, Differential Evolution -A Simple and Efficient Heuristic for Global Optimization over Continuous Spaces , pages 341-359, Journal of Global Optimization 11, Kluwer Academic Publisher, 1997.
Les données du système mécanique comprennent ici les spécifications de celui-ci et des résultats d'essais ou simulations menés sur une boucle d'asservissement comme celle représentée à la figure 5 dans le cas d'un système actionné par un actionneur de type moteur à courant continu déplaçant une charge supposée être une inertie pure J. On comprend que la boucle d'asservissement reçoit en entrée la consigne Oc et fournit la sortie Or. A partir de la consigne Oc sont déterminés le produit Js2 qui est le couple nécessaire pour mettre en mouvement la charge d'un angle û, et l'erreur s entre la consigne Oc et la sortie Or. L'erreur s est utilisée dans le correcteur K(s) qui fournit la commande correspondant au couple de Coulomb rc duquel est soustrait le frottement total rs() extrait de la sortie. Le résultat de cette soustraction est divisé par le produit Js2 pour obtenir la sortie Or.
Les données sont obtenues pour différents scénarios (correspondant à différentes conditions d'utilisation,
plusieurs fois successivement pour affiner la détermination des paramètres de la modélisation. La détermination des paramètres de réglage de l'algorithme est connue en elle-même. On pourra utiliser les paramètres de réglage standards de l'algorithme à
évolution différentielle tels que ceux mentionnés dans le document R. Storn et K. Price, Differential Evolution -A Simple and Efficient Heuristic for Global Optimization over Continuous Spaces , pages 341-359, Journal of Global Optimization 11, Kluwer Academic Publisher, 1997.
Les données du système mécanique comprennent ici les spécifications de celui-ci et des résultats d'essais ou simulations menés sur une boucle d'asservissement comme celle représentée à la figure 5 dans le cas d'un système actionné par un actionneur de type moteur à courant continu déplaçant une charge supposée être une inertie pure J. On comprend que la boucle d'asservissement reçoit en entrée la consigne Oc et fournit la sortie Or. A partir de la consigne Oc sont déterminés le produit Js2 qui est le couple nécessaire pour mettre en mouvement la charge d'un angle û, et l'erreur s entre la consigne Oc et la sortie Or. L'erreur s est utilisée dans le correcteur K(s) qui fournit la commande correspondant au couple de Coulomb rc duquel est soustrait le frottement total rs() extrait de la sortie. Le résultat de cette soustraction est divisé par le produit Js2 pour obtenir la sortie Or.
Les données sont obtenues pour différents scénarios (correspondant à différentes conditions d'utilisation,
9 notamment températures, et différents paramètres d'alimentation pour chacune de ces conditions d'utilisation) et comprennent par exemple : paramètres (tension, intensité) d'alimentation de l'actionneur électrique, vitesse de déplacement et effort fournis par l'actionneur en fonction des paramètres d'alimentation, vitesse de déplacement et effort à obtenir en sortie du système mécanique_ les signaux étant échantillonnés.
On utilise ensuite l'algorithme d'évolution différentielle avec la fonction de coût suivante, pour une température Tc :
NSk , = r I __ 1(0 Pi- 0 , =
k=1 s,k 1=1 dans laquelle :
- Nm est le nombre de scénarios, ¨ Ns,k est le nombre d'échantillons, - O, est la sortie du modèle simulé avec le vecteur de décision (Go,i, G2, rc,i, r,i, cos,i) avec i variant de 1 à M de sorte qu'un vecteur de (5M+1) paramètres doit être ajusté.
Quand les paramètres du modèle M-LuGre ont été déterminés pour les différents scénarios, il est possible de mettre en place une simulation sur la base de la boucle d'asservissement de la figure 4 afin de déterminer la commande et l'erreur échantillonnée pour un jeu donné de variables ks,i, ks,m.
On utilise à cette fin à
nouveau l'algorithme d'évolution différentielle mais avec la fonction de coût suivante :
- I(Errork2N+ ,2 pi) , k=1 N
s,k PI(Fc,k Fc,k P-11)2 1=2 dans laquelle :
- NT est le nombre de scénarios testés, 5 - Ns,k est le nombre d'échantillons, - r,k est la commande, - a et 13 sont des nombres réels positifs choisis pour harmoniser les ordres de grandeur ou donner plus d'importance à l'erreur, la commande ou la variation de
On utilise ensuite l'algorithme d'évolution différentielle avec la fonction de coût suivante, pour une température Tc :
NSk , = r I __ 1(0 Pi- 0 , =
k=1 s,k 1=1 dans laquelle :
- Nm est le nombre de scénarios, ¨ Ns,k est le nombre d'échantillons, - O, est la sortie du modèle simulé avec le vecteur de décision (Go,i, G2, rc,i, r,i, cos,i) avec i variant de 1 à M de sorte qu'un vecteur de (5M+1) paramètres doit être ajusté.
Quand les paramètres du modèle M-LuGre ont été déterminés pour les différents scénarios, il est possible de mettre en place une simulation sur la base de la boucle d'asservissement de la figure 4 afin de déterminer la commande et l'erreur échantillonnée pour un jeu donné de variables ks,i, ks,m.
On utilise à cette fin à
nouveau l'algorithme d'évolution différentielle mais avec la fonction de coût suivante :
- I(Errork2N+ ,2 pi) , k=1 N
s,k PI(Fc,k Fc,k P-11)2 1=2 dans laquelle :
- NT est le nombre de scénarios testés, 5 - Ns,k est le nombre d'échantillons, - r,k est la commande, - a et 13 sont des nombres réels positifs choisis pour harmoniser les ordres de grandeur ou donner plus d'importance à l'erreur, la commande ou la variation de
10 la commande. Par exemple : si a et 13 sont très inférieurs à 1, l'erreur est amplifiée ; si 13 est nul, la dérivée de la commande est élevée.
On a illustré aux figures 3 et 4 un compensateur de type FEEDBACK pour un système actionné par un actionneur de type moteur à courant continu déplaçant une charge supposée être une inertie pure J.
On comprend que le compensateur reçoit en entrée la consigne Oc à partir de laquelle sont déterminés le produit Js2 et l'erreur s entre la consigne Oc et la sortie Or. Le compensateur a ici une structure de compensation connue en elle-même (De Wit et al., "A new model for control of systems with friction", pages 419 à
425, n 3, volume 40, IEEE Transactions on automatic control, 1995, IEEE). L'erreur s est utilisée dans le correcteur K(s) et dans le compensateur PM dont les sorties sont ajoutées au produit Js2 pour obtenir le couple de Coulomb Pc duquel est soustrait le frottement total Ps() résultant de la modélisation. Le résultat est divisé par le produit Js2 pour obtenir la sortie Or.
On a illustré aux figures 3 et 4 un compensateur de type FEEDBACK pour un système actionné par un actionneur de type moteur à courant continu déplaçant une charge supposée être une inertie pure J.
On comprend que le compensateur reçoit en entrée la consigne Oc à partir de laquelle sont déterminés le produit Js2 et l'erreur s entre la consigne Oc et la sortie Or. Le compensateur a ici une structure de compensation connue en elle-même (De Wit et al., "A new model for control of systems with friction", pages 419 à
425, n 3, volume 40, IEEE Transactions on automatic control, 1995, IEEE). L'erreur s est utilisée dans le correcteur K(s) et dans le compensateur PM dont les sorties sont ajoutées au produit Js2 pour obtenir le couple de Coulomb Pc duquel est soustrait le frottement total Ps() résultant de la modélisation. Le résultat est divisé par le produit Js2 pour obtenir la sortie Or.
11 Il en résulte une loi de commande incorporant une compensation des frottements qui permet de piloter l'actionneur de manière efficace, régulière et stabilisante car la structure de compensation du LuGre original est conservée.
Bien entendu, l'invention n'est pas limitée au mode de réalisation décrit mais englobe toute variante entrant dans le champ de l'invention telle que définie par les revendications.
En particulier, le procédé de l'invention est utilisable pour la commande d'un actionneur de n'importe quel système mécanique. Un tel système mécanique est destiné
au déplacement précis d'un élément mobile par exemple pour le pointage d'un organe optique ou optronique dans une direction donnée. Le système mécanique peut comprendre des pièces en rotation et/ou en translation.
Par modèle de LuGre, on entend le modèle de LuGre tel que défini initialement, ou ses dérivés comme les modèles de Dahl et Coulomb.
La structure de compensation peut être d'un autre type que le type FEEDBACK et par exemple de type FEEDFORWARD.
Un autre algorithme stochastique ou un algorithme autre qu'un algorithme stochastique peut être utilisé. On peut par exemple utiliser l'algorithme du simplexe, l'algorithme dit pattern search , l'algorithme dit subgradient method , l'algorithme dit coordinate descent ...
L'utilisation d'un algorithme d'optimisation globale sans gradient (ou derivative free optimization algorithm dans la littérature) facilite la mise en uvre du procédé
mais n'est pas obligatoire.
Bien entendu, l'invention n'est pas limitée au mode de réalisation décrit mais englobe toute variante entrant dans le champ de l'invention telle que définie par les revendications.
En particulier, le procédé de l'invention est utilisable pour la commande d'un actionneur de n'importe quel système mécanique. Un tel système mécanique est destiné
au déplacement précis d'un élément mobile par exemple pour le pointage d'un organe optique ou optronique dans une direction donnée. Le système mécanique peut comprendre des pièces en rotation et/ou en translation.
Par modèle de LuGre, on entend le modèle de LuGre tel que défini initialement, ou ses dérivés comme les modèles de Dahl et Coulomb.
La structure de compensation peut être d'un autre type que le type FEEDBACK et par exemple de type FEEDFORWARD.
Un autre algorithme stochastique ou un algorithme autre qu'un algorithme stochastique peut être utilisé. On peut par exemple utiliser l'algorithme du simplexe, l'algorithme dit pattern search , l'algorithme dit subgradient method , l'algorithme dit coordinate descent ...
L'utilisation d'un algorithme d'optimisation globale sans gradient (ou derivative free optimization algorithm dans la littérature) facilite la mise en uvre du procédé
mais n'est pas obligatoire.
Claims (7)
1. Procédé de commande d'un actionneur électrique d'un système mécanique comprenant une pluralité de zones de contact imbriquées, le procédé comprenant les étapes de :
- acquérir des données sur le système mécanique qui comprend un nombre de zones de contact imbriquées, - élaborer une modélisation du système, à partir desdites données et d'un nombre de modèles de LuGre égal au nombre de zones de contact imbriquées, et déterminer des paramètres de la modélisation et une structure de compensation des frottements modélisés, - inclure la structure de compensation dans une loi de commande de l'actionneur, et - commander l'actionneur au moyen de la loi de commande.
- acquérir des données sur le système mécanique qui comprend un nombre de zones de contact imbriquées, - élaborer une modélisation du système, à partir desdites données et d'un nombre de modèles de LuGre égal au nombre de zones de contact imbriquées, et déterminer des paramètres de la modélisation et une structure de compensation des frottements modélisés, - inclure la structure de compensation dans une loi de commande de l'actionneur, et - commander l'actionneur au moyen de la loi de commande.
2. Procédé selon la revendication 1, dans lequel les paramètres de la modélisation sont déterminés en utilisant un algorithme d'optimisation globale sans gradient.
3. Procédé selon la revendication 2, dans lequel l'algorithme d'optimisation globale sans gradient est un algorithme stochastique.
4. Procédé selon la revendication 3, dans lequel l'algorithme stochastique est un algorithme à évolution différentielle.
5. Procédé selon l'une quelconque des revendications 1 à 4, dans lequel les paramètres de la modélisation propre à chacun des modèles de LuGre i sont :
- raideur des lames représentant le frottement, - coefficient d'amortissement du mouvement des lames, - coefficient de frottement visqueux (.sigma.2), - couple de frottement de Coulomb (.GAMMA.c,i), - couple de frottement statique (.GAMMA.S,i) - vitesse angulaire de Stribeck (.omega.S,i), et - vitesse relative globale du moteur (.omega.).
- raideur des lames représentant le frottement, - coefficient d'amortissement du mouvement des lames, - coefficient de frottement visqueux (.sigma.2), - couple de frottement de Coulomb (.GAMMA.c,i), - couple de frottement statique (.GAMMA.S,i) - vitesse angulaire de Stribeck (.omega.S,i), et - vitesse relative globale du moteur (.omega.).
6. Procédé selon l'une quelconque des revendications 1 à 5, dans lequel la structure de compensation est de type FEEDBACK.
7. Procédé selon l'une quelconque des revendications 1 à 5, dans lequel la structure de compensation est de type FEEDFORWARD.
Date Reçue/Date Received 2022-04-27
Date Reçue/Date Received 2022-04-27
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FRFR1912230 | 2019-10-30 | ||
| FR1912230A FR3102866B1 (fr) | 2019-10-30 | 2019-10-30 | Procédé de commande d’un actionneur de système mécanique à frottements imbriqués |
| PCT/EP2020/079953 WO2021083814A1 (fr) | 2019-10-30 | 2020-10-23 | Procede de commande d'un actionneur de systeme mecanique a frottements imbriques |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CA3159341A1 CA3159341A1 (fr) | 2021-05-06 |
| CA3159341C true CA3159341C (fr) | 2023-01-03 |
Family
ID=69743349
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CA3159341A Active CA3159341C (fr) | 2019-10-30 | 2020-10-23 | Procede de commande d'un actionneur de systeme mecanique a frottements imbriques |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US11803171B2 (fr) |
| CN (1) | CN114641738B (fr) |
| CA (1) | CA3159341C (fr) |
| FR (1) | FR3102866B1 (fr) |
| WO (1) | WO2021083814A1 (fr) |
Family Cites Families (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5233852A (en) * | 1992-04-15 | 1993-08-10 | Aluminum Company Of America | Mill actuator reference adaptation for speed changes |
| US6505703B2 (en) * | 2001-03-14 | 2003-01-14 | Visteon Global Technologies, Inc. | Vehicle steering system control |
| AU2008202736B2 (en) * | 2002-03-21 | 2011-09-15 | Alexandria Nominees Pty Ltd | Incremental material urging system |
| US6901320B2 (en) * | 2003-01-15 | 2005-05-31 | Visteon Global Technologies, Inc. | Friction compensation in a vehicle steering system |
| US7276877B2 (en) * | 2003-07-10 | 2007-10-02 | Honeywell International Inc. | Sensorless control method and apparatus for a motor drive system |
| JP4920612B2 (ja) * | 2008-02-07 | 2012-04-18 | 株式会社ハーモニック・ドライブ・システムズ | アクチュエータの角度伝達誤差補償方法 |
| CA2779013C (fr) * | 2009-11-09 | 2018-03-13 | Advanced Mechanical Technology, Inc. | Systeme et procede pour une simulation de mouvement de joint |
| JP5453606B2 (ja) * | 2010-04-09 | 2014-03-26 | 株式会社ハーモニック・ドライブ・システムズ | 波動歯車装置を備えたアクチュエータの位置決め制御方法 |
| JP5574228B2 (ja) * | 2010-04-09 | 2014-08-20 | 株式会社ハーモニック・ドライブ・システムズ | 波動歯車減速機の温度変化に伴う摩擦特性変動を考慮したアクチュエータの適応型摩擦補償法 |
| US8427093B2 (en) * | 2010-07-02 | 2013-04-23 | Woodward Hrt, Inc. | Controller for actuation system employing Kalman estimator incorporating effect of system structural stiffness |
| CN103398154A (zh) * | 2013-08-10 | 2013-11-20 | 韶能集团韶关宏大齿轮有限公司 | 液压驱动式离合器与制动器结构的行星式自动换挡变速器 |
| FR3016010B1 (fr) * | 2013-12-30 | 2016-02-12 | Chassis Brakes Int Bv | Dispositif de frein a tambour incluant un frein de stationnement fonctionnant en mode duo servo, vehicule et procedes d'assemblage associes |
| US20170261529A1 (en) * | 2014-06-26 | 2017-09-14 | Hiwin Technologies Corp. | Method for identifying friction parameter for linear module |
| US9702187B2 (en) * | 2015-02-13 | 2017-07-11 | Hunter Douglas Inc. | Covering for an architectural opening having nested tubes |
| EP3144774A1 (fr) * | 2015-09-15 | 2017-03-22 | Thomson Licensing | Procédés et appareil de transmission d'une information d'angle de rotation à un ensemble d'actionneurs associés à une surface |
| DE102019002644A1 (de) * | 2018-04-17 | 2019-10-17 | Fanuc Corporation | Steuerung und Steuerverfahren |
-
2019
- 2019-10-30 FR FR1912230A patent/FR3102866B1/fr active Active
-
2020
- 2020-10-23 WO PCT/EP2020/079953 patent/WO2021083814A1/fr not_active Ceased
- 2020-10-23 CN CN202080076810.9A patent/CN114641738B/zh active Active
- 2020-10-23 US US17/773,315 patent/US11803171B2/en active Active
- 2020-10-23 CA CA3159341A patent/CA3159341C/fr active Active
Also Published As
| Publication number | Publication date |
|---|---|
| FR3102866A1 (fr) | 2021-05-07 |
| CN114641738A (zh) | 2022-06-17 |
| FR3102866B1 (fr) | 2021-10-01 |
| US20220382242A1 (en) | 2022-12-01 |
| WO2021083814A1 (fr) | 2021-05-06 |
| CA3159341A1 (fr) | 2021-05-06 |
| CN114641738B (zh) | 2024-06-28 |
| US11803171B2 (en) | 2023-10-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP3472676B1 (fr) | Procédé de compensation des couples de coriolis, centrifuges et de gravité dans un simulateur de mouvements et système à simulateur de mouvements | |
| EP3446067B1 (fr) | Procédé d'estimation du jeu d'un actionneur électro-mécanique | |
| EP2956276A1 (fr) | Procede de detection amelioree de collision d'un robot avec son environnement, systeme et produit programme d'ordinateur mettant en uvre le procede | |
| FR3012112A1 (fr) | Procede de surveillance de fonctionnement d'un dispositif de pilotage d'aeronef et dispositif de pilotage d'aeronef ainsi surveille | |
| EP3221757A1 (fr) | Procédé et système de compensation adaptative de frottements secs | |
| EP4054804B1 (fr) | Articulation à trois degrés de liberté pour un robot et procédé de commande correspondant | |
| FR3040220A1 (fr) | Systeme de regulation d'un parametre regule | |
| WO2015193592A1 (fr) | Estimateur de couple pour double embrayage | |
| CA3159341C (fr) | Procede de commande d'un actionneur de systeme mecanique a frottements imbriques | |
| CA2255996C (fr) | Actionneur de gouverne de commande de vol d'aeronef | |
| FR2900360A1 (fr) | Procede et dispositif de reglage de parametres de fonctionnement d'un robot, programme et support d'enregistrement pour ce procede | |
| WO2009065818A1 (fr) | Actionneur a transfert de moment cinetique pour le controle d'attitude d'un engin spatial | |
| FR2960808A1 (fr) | Procede d'identification de frottements dans une articulation de bras de robot ou de bras manipulateur, et procede de compensation de couple faisant application | |
| EP3391161A1 (fr) | Procédé de surveillance d'un système d'actionnement électromécanique | |
| WO2015040238A1 (fr) | Procédé et dispositif de réglage d'un système mécanique actionné | |
| WO2022129566A1 (fr) | Robot de co-manipulation à loi de commande en effort mixte offrant une forte sensibilité de l'effecteur et permettant l'interaction avec le corps du robot | |
| EP3661702B1 (fr) | Procede d'assistance d'au moins un mouvement d'un utilisateur et dispositif correspondant | |
| JP7058868B2 (ja) | 旋回台の制御装置およびその設定方法 | |
| WO2022254162A1 (fr) | Système de commande de compensateur de vol pour aéronef à retour haptique | |
| Zeng et al. | Optimal Controller Design and Simulation Analysis of Inertially Stabilized Platform for Airborne Remote Sensing | |
| FR2882118A1 (fr) | Procede et systeme d'asservissement d'un embrayage pilote et procede de fabrication du systeme | |
| FR2855282A1 (fr) | Dispositif et procede de commande d'un retour de force a appliquer a au moins une manette d'une interface motorisee |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| EEER | Examination request |
Effective date: 20220427 |
|
| EEER | Examination request |
Effective date: 20220427 |
|
| EEER | Examination request |
Effective date: 20220427 |
|
| EEER | Examination request |
Effective date: 20220427 |
|
| EEER | Examination request |
Effective date: 20220427 |
|
| EEER | Examination request |
Effective date: 20220427 |
|
| MPN | Maintenance fee for patent paid |
Free format text: FEE DESCRIPTION TEXT: MF (PATENT, 4TH ANNIV.) - STANDARD Year of fee payment: 4 |
|
| U00 | Fee paid |
Free format text: ST27 STATUS EVENT CODE: A-4-4-U10-U00-U101 (AS PROVIDED BY THE NATIONAL OFFICE); EVENT TEXT: MAINTENANCE REQUEST RECEIVED Effective date: 20240930 |
|
| U11 | Full renewal or maintenance fee paid |
Free format text: ST27 STATUS EVENT CODE: A-4-4-U10-U11-U102 (AS PROVIDED BY THE NATIONAL OFFICE); EVENT TEXT: MAINTENANCE FEE PAYMENT DETERMINED COMPLIANT Effective date: 20241202 |