Application de La Logique Floue - Commande de La Température PDF
Application de La Logique Floue - Commande de La Température PDF
Application de La Logique Floue - Commande de La Température PDF
a logique floue a été introduite en 1965 par L.A. Zadeh [1]. Ses principes ont
L été appliqués en 1974 par E.H. Mamdani [2] à la construction d’un premier
contrôleur flou. Mais ce n’est que depuis quelques années que la commande
floue a connu, essentiellement au Japon, un essor remarquable ; elle a été
7 - 1993
Toute reproduction sans autorisation du Centre français d’exploitation du droit de copie est strictement interdite.
© Techniques de l’Ingénieur, traité Informatique industrielle R 7 428 − 1
APPLICATION DE LA LOGIQUE FLOUE : COMMANDE DE LA TEMPÉRATURE D’UN FOUR _______________________________________________________________
1. Principe et réalisation
du contrôleur flou
1.1 Généralités sur les ensembles flous
1.1.1 Définition
Si l’on considère une grandeur physique e, dire que e appartient
à un certain ensemble flou revient à lui attribuer une propriété de
définition imprécise (« linguistique ») : e est grand, petit, voisin de
zéro, etc. La frontière d’un tel ensemble est mal définie, si bien
qu’un élément peut appartenir à la fois à un ensemble et à son
complémentaire. Figure 1 – Fonctions d’appartenance du paramètre e
aux sous-ensembles flous E j
Toute reproduction sans autorisation du Centre français d’exploitation du droit de copie est strictement interdite.
R 7 428 − 2 © Techniques de l’Ingénieur, traité Informatique industrielle
_______________________________________________________________ APPLICATION DE LA LOGIQUE FLOUE : COMMANDE DE LA TEMPÉRATURE D’UN FOUR
Les (2N + 1) ensembles flous E 1 , ..., E j ... ont donc pour bornes : Par exemple, il est clair que, si l’écart est fortement négatif et son
accroissement aussi, le signal de commande doit l’être également.
1 2 1 3 Au contraire, si l’écart est petit et son accroissement aussi, il en sera
Ð ∞,Ð 1 + -------- , Ð 1 , Ð 1 + ------- , Ð 1 + ------- , Ð 1 + -------- , ... de même de la commande. Si l’écart est négatif et petit en module
N N N N
mais son accroissement fortement négatif, ou si l’écart est fortement
et pour « nombres centraux » : négatif mais son accroissement négatif et petit en module, le signal
de commande devra être négatif, mais de module moyen.
jÐ1 Ces considérations nous ont conduits à adopter une table de
v 1 = Ð 1, ..., v j = Ð 1 + ------------ , ..., v 2 N + 1 = 1
N décision « anti-diagonale », résumant les règles choisies [5] ; le
■ On associe à chacun de ces ensembles flous Ej une fonction tableau 1 montre une telle table de décision, en indiquant la
caractéristique de la sortie en fonction de celles de l’écart et de son
d’appartenance µ Ej ( v ) , centrée sur vj , pour laquelle on choisit une accroissement dans le cas où N = 1 ; le tableau 2 généralise au cas
forme triangulaire (figure 3) ; dans le cas général, les fonctions de N quelconque. Ce type de table impose un nombre d’ensembles
d’appartenance sont données par : flous égal à 4N + 1 pour la variable de sortie s.
On choisit, pour les ensembles flous relatifs à la variable de sortie
µ Ej ( v ) = 0 pour v < v j Ð 1 ou v > v j + 1 (2) S1 , S 2 , ..., S4 N + 1 , des fonctions d’appartenance triangulaires
(comme pour les ensembles flous relatifs aux variables d’entrée),
Ð v + vj + 1 définies dans le domaine de variation [– 1, + 1] et centrées sur les
µ Ej ( v ) = ----------------------------- pour v j < v < v j + 1 (3) nombres s k ( 1 < k < 4 N + 1 ) :
vj + 1 Ð vj
1 kÐ1
+ v Ð vj Ð 1 s 1 = Ð 1 , s 2 = – 1 + ---------- , ... , s k = Ð 1 + -------------- , ... , s 4 N + 1 = 1
µ Ej ( v ) = ----------------------------- pour v j Ð 1 < v < v j (4) 2N 2N
vj Ð vj Ð 1
Les règles R adoptées sont donc les suivantes :
Toutefois, pour v < – 1, µ E1 (v ) = 1 et les autres µ Ej ( v ) sont nuls — règle R1 : si l’écart appartient à l’ensemble E 1 et son
et, pour v > 1, µ E2 N + 1 ( v ) = 1 et les autres µ Ej ( v ) sont nuls. accroissement à E 1′ , alors la sortie appartient à l’ensemble flou S1 ;
— règle R k ( 1 < k < 2 N + 1 ) : si l’écart appartient à E k et son
De même, l’intervalle de variation de la deuxième variable (accrois-
sement de l’écart) est ramené à [– 1, + 1] ; on décompose cet inter- accroissement à E 1′ , ou si l’écart appartient à E k – 1 et son
valle en 2 N + 1 ensembles flous E j′ , on définit une variable accroissement à E 2′ , ou si l’écart appartient à E k – 2 et son
normalisée v ’ dont les coefficients d’appartenance µ E ′j ( v ′ ) à ces accroissement à E 3′ ..., alors la sortie appartient à l’ensemble flou Sk ;
2N + 1 ensembles flous sont donnés en fonction de v ’ par les mêmes — règle R k ( 2 N + 1 < k < 4 N + 1 ) : si l’écart appartient à E 2N + 1 et
relations que les µ Ej ( v ) en fonction de v .
son accroissement à E k′ Ð 2 N ou si l’écart appartient à E 2N et son
accroissement à E k + 1 – 2N .... , alors la sortie appartient à l’ensemble
1.2.2 Règles floues flou Sk ;
— règle R 4N + 1 : si l’écart appartient à E 2N + 1 et son accrois-
Ces règles permettent de déterminer le signal de sortie s du sement à E 2′ N + 1 , alors la sortie appartient à l’ensemble flou S 4N + 1 .
contrôleur en fonction des signaux d’entrée e et e ’. Elles relient s
à e et e ’ par des conditions linguistiques prenant en compte
l’expérience ou le savoir-faire acquis par un opérateur humain, ou
traduisant simplement des remarques de bon sens.
négatif
s grande s négative s voisine
négative de zéro On généralise facilement la méthode au calcul du coefficient
d ’ a p p a rt e n a n c e à ch a c u n e d e s r è g l e s R k é n o n c é e s a u
voisin de zéro s négative s voisine s positive paragraphe 1.2.2.
de zéro
µR 1 ( e , e ′, s ) = min µ E1 ( e ), µ E ′ ( e ′ ), µ S1 ( s )
1
= min δ 1 ( e , e ′ ), µ S1 ( s ) (5)
Toute reproduction sans autorisation du Centre français d’exploitation du droit de copie est strictement interdite.
R 7 428 − 4 © Techniques de l’Ingénieur, traité Informatique industrielle
_______________________________________________________________ APPLICATION DE LA LOGIQUE FLOUE : COMMANDE DE LA TEMPÉRATURE D’UN FOUR
Il reste à déterminer, à partir des valeurs envisageables, la vera une description plus détaillée de l’appareillage et de la méthode
valeurs s 0 effective de la sortie du contrôleur à l’instant où on a d’identification dans la référence [7].
mesuré e 0 et e 0′ : c’est l’opération de « défuzzification » (de fuzzy :
flou).
Divers choix sont possibles. Par exemple, on peut prendre pour 2.1 Dispositif expérimental
s 0 le nombre central s k correspondant au plus grand δk (Mamdani),
ou l’abscisse du centre de gravité de la surface limitée par l’ensemble Dans une chaîne de fabrication de pâtisseries industrielles, les
des trapèzes, etc. Pour nos expériences, nous avons choisi la moules où est déposée la pâte passent sur un tapis roulant dans
méthode de Sugeno qui prend pour s 0 le barycentre des s k affectés un four constitué par un long tunnel (20 à 30 m), dans chaque section
des coefficients δk : duquel la température est maintenue à une valeur constante, déter-
4N + 1 minée par le profil de cuisson à respecter. Il ne saurait être question
∑ sk δk d’utiliser un tel four pour la mise au point de produits nouveaux.
1 Les essais sont toujours faits dans un four statique et leur repro-
s0 = ----------------------------
- (10)
4N + 1 ductibilité laisse à désirer, les conditions de chauffage et de refroi-
∑ δk dissement étant réglées manuellement de manière peu précise.
L’objectif de l’essai est de réguler la température d’un four de ce type
Dans nos expériences, le processus à commander ne comporte ou d’obtenir la poursuite d’un profil de température en fonction du
pas d’intégrateur. Pour éviter l’erreur de position, on en introduit temps fixé à l’avance.
un en série avec le contrôleur flou ; la commande u (t ) est donc Ce four est constitué (figure 5) d’une chambre de cuisson en
donnée, pendant la période d’échantillonnage [nT, nT + T ], par : acier, d’épaisseur 2 mm, de 620 mm de largeur, 280 mm de hau-
teur et 880 mm de profondeur, entouré d’un calorifugeage en laine
u (nT ) = u (nT – T ) + Ks 0 (11) de verre capoté de tôle peinte. La pâte est introduite par une porte
où K est un facteur de gain préalablement choisi. en verre qui permet de suivre l’évolution de la cuisson. Les buées
sont évacuées par une cheminée de diamètre 40 mm. Le chauffage
est assuré par 3 résistances électriques couplées en voûte (1 750 W
1.2.6 Résumé de puissance totale sous 220 V) et 3 autres couplées en sole
(2 250 W).
En définitive, les paramètres du correcteur sont les limites des La constante de temps du four est très élevée, de sorte qu’il
domaines de variation [a, b ] et [a’, b ’] de l’écart e et de son serait impossible de suivre les profils de cuisson imposés lorsque
accroissement e ’, le nombre N et le gain K . Son fonctionnement ceux-ci comportent une étape de refroidissement rapide. Il n’est
conduit à : pas possible d’abaisser la température en introduisant de l’air frais,
car on modifierait ainsi le « microclimat » qui se crée autour du
biscuit pendant le chauffage (taux d’humidité en particulier). Nous
— mesurer à chaque instant d’échantillonnage les valeurs e 0 avons donc utilisé une circulation d’eau dans des radiateurs à eau
et e 0′ de e et e ’ ; perdue, constitués de 40 m de tube de cuivre de 10 mm de
— calculer les valeurs normalisées v 0 et v 0′ correspondantes diamètre placés dans le calorifugeage, et 10 m à l’intérieur de la
à l’aide de la relation (1) ; chambre de cuisson entre la tôle et les résistances chauffantes.
— déterminer tous les coefficients d’appartenance µ Ei ( e 0 ) et Dans les essais de commande décrits ci-après, nous avons opéré
µ E j′ ( e 0′ ) ; avec un débit d’eau permanent de 0,6 L /minute. En effet, on
diminue ainsi la constante de temps du système, qui serait énorme
— remplir le tableau à 2N + 1 lignes et 2N + 1 colonnes, sans circulation d’eau (plus de 3 h). De plus, on doit utiliser un
associées aux ensembles Ei , E j′ par les nombres : débit de 4 L /minute pour les phases de refroidissement rapide et,
si l’on passait brusquement d’un débit nul à un débit de 4 L /minute
min µ Ei (e 0 ), µ E j′ ( e 0′ ) ; lorsque le four est chaud, l’eau se transformerait en jet de vapeur,
— prendre, sur chaque « anti-diagonale » telle que i + j soit ce qui hacherait le débit et ne permettrait pas d’obtenir un
constant, le maximum δk des nombres qui s’y trouvent ; fonctionnement reproductible. Le passage d’une valeur du débit à
l’autre se fait par l’ouverture ou la fermeture d’une électrovanne.
(k Ð 1)
— prendre le barycentre s 0 des nombres s k = Ð 1 + ------------------- Le chauffage du four est commandé par des gradateurs de
2N
affectés des coefficients δk (k = 1, ..., 4N + 1) ; puissance à triac (Eurotherm 240 V, 25A) qui hachent le signal
— multiplier ce résultat par le gain K pour obtenir sinusoïdal du réseau. La température est déterminée à l’aide d’un
l’accroissement à donner à la commande [relation (11)]. thermocouple Fe /Cu-Ni chemisé, de diamètre 2 mm. La précision
est de quelques dixièmes de degré de 0 à 40 oC, et de 1 à 1,5 oC
pour des températures supérieures à 100 oC. Le signal de quelques
Remarque : une caractéristique importante d’un tel contrôleur millivolts du thermocouple est transformé en un signal de
est sa non-linéarité. Il est clair en effet que, si le signal e par quelques volts à l’aide d’un convertisseur (type EMO 10 t 30/30GT)
exemple sort de l’intervalle [ a , b ], seul le coefficient qui est, en fait, un amplificateur parfaitement linéaire.
d’appartenance à l’ensemble flou E 1 ou E 2N + 1 est non nul ; dans Les programmes de commande ont été rédigés en Turbo-Pascal
ce cas, la sortie du contrôleur ne dépend plus de la valeur de e et version 4.0. Nous avons utilisé un calculateur IBM. XT286 muni
il en est de même de l’accroissement de la commande. d’une carte comportant des convertisseurs analogique-numérique
et numérique-analogique (carte DAC) dont les entrées-sorties
permettent de lire la température et de commander les gradateurs
et l’électrovanne.
2. Expériences
2.2 Identification du processus
Les expériences décrites dans ce paragraphe ont été réalisées au
Laboratoire d’Automatique et de Contrôle des Processus de l’ISMRa Une identification du système est indispensable lorsqu’on utilise
par MM. CONDETTE et DORLÉANS [6] [7] [8], que les auteurs des algorithmes de correction classiques. Elle peut paraître
remercient pour leur collaboration active et efficace. Le lecteur trou-
Toute reproduction sans autorisation du Centre français d’exploitation du droit de copie est strictement interdite.
© Techniques de l’Ingénieur, traité Informatique industrielle R 7 428 − 5
APPLICATION DE LA LOGIQUE FLOUE : COMMANDE DE LA TEMPÉRATURE D’UN FOUR _______________________________________________________________
Toute reproduction sans autorisation du Centre français d’exploitation du droit de copie est strictement interdite.
R 7 428 − 6 © Techniques de l’Ingénieur, traité Informatique industrielle
_______________________________________________________________ APPLICATION DE LA LOGIQUE FLOUE : COMMANDE DE LA TEMPÉRATURE D’UN FOUR
Les résultats obtenus sont de même très médiocres dans le cas De plus, une expérience a été réalisée récemment dans notre
du deuxième profil de cuisson, constitué de paliers séparés [8]. laboratoire pour le contrôle d’un système non linéaire instable
(asservissement de position d’un ludion dans une cuve remplie
■ L’emploi d’une commande optimale avec intégrateur utilisant un d’eau), très difficile à commander par les techniques traditionnelles
critère quadratique (commande LQI [9]) donne des résultats tout à de l’automatique. Les résultats ont été excellents, le système
fait similaires. Pour obtenir une régulation et une poursuite répondant de manière parfaitement satisfaisante à des consignes
satisfaisantes et moins sensibles aux variations des paramètres du en échelon ou en succession de rampes constituant des triangles.
système, il faut utiliser des algorithmes beaucoup plus sophistiqués
du point de vue mathématique, comme la commande GPC
(generalized predictive control ) [10] qui est nettement moins facile à
mettre en œuvre.
3. Conclusion
2.6 Discussion Les résultats que nous avons décrits dans cet article sont très satis-
faisants. Faut-il en conclure que la commande floue est capable
■ Si l’on ne considère que l’expérience décrite ici, le contrôleur flou d’apporter une solution simple et rapide à tout problème de
utilisé semble présenter des avantages multiples. Il ne nécessite pas commande pouvant se poser dans l’industrie ?
un formalisme mathématique difficile et sa programmation est aisée
Notre pratique de la commande floue reste encore assez limitée,
(il faut toutefois noter qu’il impose un temps de calcul relativement
mais nous pensons cependant pouvoir affirmer que cette conclusion
long, alors que les algorithmes traditionnels – éventuellement plus
serait à la fois prématurée et probablement trop optimiste. Le
élaborés du point de vue mathématique – se traduisent souvent par
contrôleur que nous avons décrit peut être qualifié de « standard »
des programmes très courts). Il peut sembler, à première vue, que le
ou « passe-partout », dans la mesure où il n’a pas été choisi pour
contrôleur flou ne nécessite que très peu de connaissance sur le sys-
s’adapter spécifiquement au problème posé. Nous l’avons d’ailleurs
tème à asservir. Enfin, ce contrôleur est peu sensible aux variations
utilisé pour un système relativement facile à commander (sa réponse
des paramètres du système ; sa robustesse paraît donc bonne,
en boucle ouverte ne comporte pas d’oscillations). Néanmoins,
meilleure par exemple que celle d’un contrôleur LQI testé dans les
même dans ce cas favorable, le réglage du contrôleur a demandé
mêmes conditions [6].
des essais en simulation assez longs. Les simulations que nous
■ En réalité, il ne faut pas dissimuler un certain nombre avons faites ont montré qu’il en était de même dans le cas d’autres
d’inconvénients inhérents à ce type de contrôleur. Les essais réalisés systèmes, ce qui peut nécessiter des tâtonnements avant qu’on
montrent qu’il est imprudent de faire fonctionner un contrôleur flou n’arrive à un fonctionnement dont on ne peut être sûr qu’il soit opti-
sans l’avoir testé au préalable en simulation, ce qui suppose une mal. Notons cependant que les réglages sont beaucoup plus faciles
identification au moins approximative du système. De plus, la pour- si l’on supprime l’intégrateur mis en série avec le contrôleur, mais
suite de l’étude par une série de simulations sur des processus on ne peut éviter alors une erreur de position.
d’ordre 1, 2 ou 3, soumis à différentes perturbations ou comportant Par opposition à ce type de contrôleur, on peut aussi imaginer
des erreurs de modélisation [11], a montré que, si la commande floue d’utiliser un contrôleur « spécifique » , c’est-à-dire plus particu-
est effectivement assez robuste pour des systèmes du premier ordre, lièrement adapté au problème de commande à résoudre. Cette atti-
il n’en est pas de même pour des systèmes du deuxième ordre, tude est plus conforme aux principes même de la logique floue, qui
principalement lorsqu’ils ont des pôles complexes. Le contrôleur flou est destinée à traduire la connaissance et le savoir-faire humains.
utilisé donne alors des résultats assez médiocres. Avec un tel contrôleur, on peut espérer obtenir de bons résultats dans
En fait, une étude récente [12] montre qu’un contrôleur du type de très nombreux problèmes, mais il est clair que le choix de sa struc-
décrit ci-dessus, sans l’intégrateur mis en cascade, a un signal de ture – limitation des ensembles flous pour les variables d’entrée et
sortie très voisin de celui d’un contrôleur proportionnel et dérivé dont de sortie, forme des fonctions d’appartenance, tableau des règles,
les paramètres sont reliés simplement à ceux du contrôleur flou. Avec qui peut ne pas être symétrique, méthode de « défuzzification », etc.
l’intégrateur en cascade, le signal de commande est donc voisin de – et son réglage seront alors empiriques, nécessiteront encore des
celui que fournirait un correcteur proportionnel et intégral. Si l’on tâtonnements parfois longs et délicats.
met l’intégrateur K i /p en parallèle avec le contrôleur flou (figure 10), Enfin, soulignons que l’utilisation d’un intégrateur en parallèle, au
le contrôleur ainsi construit se comporte d’une manière peu dif- lieu d’un intégrateur en série, avec le contrôleur flou « standard »
férente d’un contrôleur PID. On dispose ainsi d’un paramètre sup- que nous avons décrit, ouvre des perspectives particulièrement
plémentaire, ce qui permet de contrôler plus aisément le système. prometteuses.
Les essais réalisés en simulation avec un contrôleur de ce type pour
un système du second ordre instable sont très encourageants [11].
Références bibliographiques
[1] ZADEH (L.A.). – Fuzzy sets. Inform. Control. 8, p. 338-353 (1965). [8] BARRAT (M.), DORLÉANS (P.), VILLENEUVE (G.) et LÉCLUSE (Y.). –
[2] MAMDANI (E.H.). – Application of fuzzy algorithms for control of a Partial state model reference control of an industrial oven. Proc. Instr.
simple dynamic plant. Proc. IEEE, 212, p. 1585-1588 (1974). Mech. Engrs, 206, p. 145-149 (1992).
[3] SUGENO (M.) (éditeur). – Industrial applications of fuzzy control. North [9] FOULARD (C.), GENTIL (S.) et SANDRAZ (J.P.). – Commande et régulation
Holland Pub. Comp. (1992). par calculateur numérique. Eyrolles (1982).
[4] KAUFMANN (A.). – Introduction à la logique floue. R 7032, traité [10] CLARKE (D.W.), MOHTADI (C.) et TUFFS (P.S.). – Generalized predictive
Informatique industrielle, Techniques de l’Ingénieur, oct. 1992. control. Automatica, 23, p. 137-160 (1987).
[5] BUCKLEY (J.J.) et YING (H.). – Fuzzy controller theory : limit theorems for [11] BARRAT (J.P.) et BARRAT (M.). – Perfomances d’un contrôleur flou.
linear fuzzy control rules. Automatica, 25, p. 469-472 (1989). Rapport du Groupe de recherche en Automatique (Pôle Automatisation
intégrée), juin 1992.
[6] CONDETTE (C.), DORLÉANS (P.), BARRAT (J.P.) et BARRAT (M.). –
Application de la commande floue à un four pilote. RAPA, 4, p. 219-231 [12] KAWAJI (S.), MAEDA (T.) et MATSUNAGA (N.). – Fuzzy control using
(1991). knowledge acquired from PD control. Proc. IECON, p. 1549-1554 (1991).
[7] DORLÉANS (P.), BARRAT (J.P.) et LÉCLUSE (Y.). – Commande par [13] GILLE (J.C.). – Systèmes et signaux déterministes. Transformées et
calculateur d’un procédé thermique. RAPA, 2, p. 65-78 (1989). abaques. R 7 010, traité Informatique industrielle, Techniques de
l’Ingénieur, avril 1995.
Toute reproduction sans autorisation du Centre français d’exploitation du droit de copie est strictement interdite.
© Techniques de l’Ingénieur, traité Informatique industrielle R 7 428 − 7
APPLICATION DE LA LOGIQUE FLOUE : COMMANDE DE LA TEMPÉRATURE D’UN FOUR _______________________________________________________________
Toute reproduction sans autorisation du Centre français d’exploitation du droit de copie est strictement interdite.
R 7 428 − 8 © Techniques de l’Ingénieur, traité Informatique industrielle
_______________________________________________________________ APPLICATION DE LA LOGIQUE FLOUE : COMMANDE DE LA TEMPÉRATURE D’UN FOUR
Figure 7 – Courbes expérimentales obtenues avec le four pilote corrigé par le contrôleur flou
Toute reproduction sans autorisation du Centre français d’exploitation du droit de copie est strictement interdite.
© Techniques de l’Ingénieur, traité Informatique industrielle R 7 428 − 9
APPLICATION DE LA LOGIQUE FLOUE : COMMANDE DE LA TEMPÉRATURE D’UN FOUR _______________________________________________________________
Toute reproduction sans autorisation du Centre français d’exploitation du droit de copie est strictement interdite.
R 7 428 − 10 © Techniques de l’Ingénieur, traité Informatique industrielle
_______________________________________________________________ APPLICATION DE LA LOGIQUE FLOUE : COMMANDE DE LA TEMPÉRATURE D’UN FOUR
Figure 9 – Courbes expérimentales obtenues avec le four pilote corrigé par un contrôleur PI numérique
Toute reproduction sans autorisation du Centre français d’exploitation du droit de copie est strictement interdite.
© Techniques de l’Ingénieur, traité Informatique industrielle R 7 428 − 11