Brevet FR2951606A1 — publié le 20110422

Procédé de chauffage par induction mis en oeuvre dans un dispositif comprenant des inducteurs couplés ma2nétiquement La présente invention se rapporte à un procédé de chauffage par induction mis en oeuvre dans un dispositif de chauffage d'une pièce métallique telle qu'une tôle ou une barre, le dispositif comprenant des inducteurs couplés magnétiquement. Par couplage magnétique, on entend que les inducteurs produisent entre eux des inductions mutuelles. Les techniques les plus classiques de chauffage par induction mettent en oeuvre des configurations qui sont satisfaisantes lorsque les pièces à chauffer sont toujours de même nature et de mêmes dimensions. Mais l'industrie exige de plus en plus de souplesse et de productivité. Il est demandé aux lignes de production de s'adapter en fonctionnement continu au changement de la position ou du format des pièces à chauffer, et d'adapter en fonction de ce changement le profil de température désiré. Des technologies connues permettent d'avoir un contrôle du chauffage par zone de puissance injectée, mais le contrôle du profil de température dans les zones chauffées reste lié à la conception géométrique des bobines et à leur mode d'alimentation, principalement par la variation d'amplitude des courants que l'on y injecte. La détermination de ces courants et la régulation qui en découle est fortement tributaire du couplage magnétique existant entre les bobines du fait des inductions mutuelles, chaque bobine alimentée ayant une influence sur toutes les autres. Le couplage magnétique rend le contrôle du profil de température de la pièce chauffée extrêmement délicat, sans compter qu'il peut avoir des répercutions néfastes sur les générateurs de fréquence, par exemple une casse de composants. La demande de brevet japonais JP2009032498A décrit un système pour chauffer une pièce métallique par des bobines d'induction couplées magnétiquement. Ce système utilise en série avec chaque bobine à partir de la deuxième une paire de thyristors antiparallèles pilotés par un circuit de contrôle permettant de contrôler la tension appliquée à la bobine à partir de données mesurées rendant compte des inductions mutuelles subies par la bobine. Ce -2- système présente néanmoins des inconvénients, et en particulier il permet de contrôler seulement la puissance moyenne produite par chaque bobine sans pouvoir contrôler précisément le profil de température généré par les bobines dans la pièce chauffée. La présente invention vise à résoudre ces inconvénients et à procurer un procédé de chauffage prenant en compte les nombreux couplages, d'une part entre les différents inducteurs et d'autre part entre les inducteurs et la pièce à chauffer, pour permettre de contrôler avec une bonne précision le profil de température généré par les inducteurs. A cet effet, l'invention a pour objet un procédé de chauffage par induction mis en oeuvre dans un dispositif de chauffage d'une pièce métallique, le dispositif comprenant des inducteurs couplés magnétiquement, chaque inducteur étant alimenté par un onduleur qui lui est propre et associé à un condensateur pour former un circuit oscillant, lesdits circuits oscillants possédant au moins approximativement la même fréquence de résonance, chaque onduleur étant commandé par une unité de commande de façon à faire varier l'amplitude et la phase du courant traversant l'inducteur correspondant, le dispositif comprenant en outre des moyens de détermination dudit courant ainsi que des moyens de détermination d'un profil de température effectif de ladite pièce métallique, ledit procédé comprenant les étapes suivantes : a) on compare ledit profil de température effectif à un profil de température de référence, et on calcule un profil de densité de puissance de référence que le dispositif de chauffage doit injecter dans ladite pièce pour atteindre ledit profil de température de référence ; b) à partir d'une matrice d'impédances déterminée par la connaissance des relations électromagnétiques liant lesdits inducteurs entre eux et à ladite pièce et par la connaissance de fonctions images vectorielles représentatives des relations liant les densités de courant créées par les inducteurs aux courants traversant les inducteurs, on calcule des courants cibles que doivent délivrer les onduleurs afin que les courants des inducteurs atteignent des valeurs cibles appropriées pour injecter dans ladite pièce ledit profil de densité de puissance de référence ; -3- c) on détermine les courants traversant les inducteurs pour les comparer aux dites valeurs cibles et déterminer des écarts de courants à corriger, et on envoie aux dites unités de commande des instructions de correction en fonction desdits écarts de courants afin de commander les onduleurs de façon à corriger les courants traversant les inducteurs. Grâce à ces dispositions, on obtient un pilotage précis du profil de température appliqué à la pièce chauffée, ce qui est idéal pour chauffer avec un même dispositif plusieurs pièces de tailles et de natures différentes. Dans des modes de réalisation préférés d'un procédé de chauffage selon l'invention, on a recours notamment à l'une ou l'autre des dispositions suivantes : on détermine les capacités desdits condensateurs, et on associe ladite matrice d'impédances à un vecteur des capacités ; on détermine une valeur initiale de ladite matrice d'impédances pour une température moyenne initiale donnée desdits inducteurs et de ladite pièce, puis on détermine à intervalles variables ou périodiques la matrice d'impédances modifiée pour au moins une valeur augmentée de ladite température moyenne, et on utilise ladite matrice d'impédances modifiée pour recalculer les dites valeurs cibles ; après avoir effectué successivement les étapes (a) et (b) on effectue au moins une fois l'étape (c) pour diminuer lesdits écarts de courants à corriger, puis on réitère au moins une fois les étapes (a), (b) et (c) en actualisant ledit profil de température effectif par des mesures de température en différentes zones chauffées de la pièce ; pour la détermination par calcul des dites valeurs cibles à l'étape (b), grâce à la connaissance desdites fonctions images vectorielles, on calcule des fonctions images des densités de puissance selon les caractéristiques spatiales des zones de la pièce dans lesquelles lesdites densités de puissance sont injectées, et on calcule un vecteur optimisé des courants cibles à déterminer en minimisant la différence entre chacune desdites fonctions images des densités de puissance et une fonction densité de puissance de référence correspondant audit profil de densité de puissance de référence. L'invention a aussi pour objet un dispositif de chauffage par induction comprenant : -4- des inducteurs couplés magnétiquement, chaque inducteur étant associé à un condensateur pour former un circuit oscillant, lesdits circuits oscillants possédant au moins approximativement la même fréquence de résonance ; des onduleurs alimentant chacun un inducteur qui lui est propre, chaque onduleur étant commandé par une unité de commande de façon à faire varier l'amplitude et la phase du courant traversant l'inducteur correspondant ; caractérisé en ce qu'il comprend en outre : des moyens de détermination des courants traversant les inducteurs ainsi que des moyens de détermination d'un profil de température effectif d'une pièce métallique chauffée par le dispositif ; des moyens de comparaison dudit profil de température effectif par rapport à un profil de température de référence ; des moyens de calcul d'un profil de densité de puissance de référence que le dispositif de chauffage doit injecter dans ladite pièce pour atteindre ledit profil de température de référence ; des moyens de calcul, basés sur la connaissance d'une matrice des impédances, de courants cibles que doivent délivrer les onduleurs afin que les courants des inducteurs atteignent des valeurs cibles appropriées pour injecter dans ladite pièce ledit profil de densité de puissance de référence ; des moyens de comparaison des courants traversant les inducteurs par rapport aux dites valeurs cibles, aptes à déterminer des écarts de courants à corriger, et des moyens de traitement desdits écarts de courants aptes à générer des instructions de correction envoyées aux dites unités de commande pour commander les onduleurs de façon à corriger les courants traversant les inducteurs. Dans des modes de réalisation préférés d'un dispositif de chauffage selon l'invention, on a recours notamment à l'une ou l'autre des dispositions suivantes : les onduleurs sont alimentés par une même alimentation source de courant ou source de tension, et lesdits moyens de comparaison desdits courants déterminés traversant les inducteurs comprennent des unités comparatrices recevant chacune des paramètres déterminés d'un courant traversant un inducteur et des paramètres des valeurs cibles correspondantes et étant chacune reliée à une unité de traitement desdits écarts de courants, une desdites unités comparatrices -5- recevant en outre des paramètres représentatifs de ce que délivre ladite alimentation et son unité de traitement associée étant adaptée pour générer des instructions de régulation envoyées à ladite alimentation de façon à modifier le courant ou la tension qu' elle délivre. D'autres caractéristiques et avantages ressortent de la description qui va suivre d'exemples non limitatifs de modes de réalisation, en référence aux figures dans lesquelles : La figure 1 représente schématiquement un premier exemple de dispositif de chauffage par induction dans lequel le procédé de chauffage selon l'invention peut être mis en oeuvre, appliqué au chauffage d'un disque métallique fixe. La figure 2 représente schématiquement une modélisation du système à trois inductances couplées de la figure 1, vu de l'alimentation. La figure 3 représente schématiquement le dispositif de chauffage par induction de la figure 1, appliqué au chauffage d'une tôle que l'on déplace. La figure 4 représente schématiquement un deuxième exemple de dispositif de chauffage par induction, appliqué au chauffage d'une barre métallique que l'on déplace. La figure 5 représente schématiquement un troisième exemple de dispositif de chauffage par induction, appliqué au chauffage d'une tôle que l'on déplace. La figure 6 représente schématiquement un quatrième exemple de dispositif de chauffage par induction, appliqué au chauffage d'une tôle que l'on déplace. La figure 7 représente schématiquement une fonction image de la densité de puissance calculée à partir d'un vecteur optimisé des courants permettant de minimiser la différence entre ladite fonction et une fonction référence de densité de puissance. La figure 8 représente schématiquement un premier mode de réalisation d'un dispositif de chauffage par induction selon l'invention dans lequel l'alimentation des onduleurs est une source de courant. La figure 9 représente schématiquement un second mode de réalisation d'un dispositif de chauffage par induction selon l'invention dans lequel l'alimentation des onduleurs est une source de tension. -6- Sur la figure 1, le dispositif de chauffage montré en exemple concerne une configuration de disque métallique amagnétique chauffé par flux transverse à l'aide de trois couples de bobines jumelles, ce qui a l'avantage de garder l'aspect axisymétrique du problème. Afin de s'assurer de la symétrie de l'ensemble du système, chaque bobine placée d'un côté du disque est connectée en série avec sa bobine jumelle de l'autre face pour former un seul inducteur. De cette façon, le système est invariant par rotation. En outre, afin de travailler avec l'hypothèse de linéarité, on considérera que les matériaux électromagnétiques du système ont une perméabilité constante et unitaire. Chaque inducteur est alimenté par un onduleur qui lui est propre de type série (onduleur de tension) ou de type parallèle (onduleur de courant). Sur la figure 2, la modélisation du système sous forme d'inductances couplées permet de représenter les différentes interactions existantes. Cette modélisation permet également l'étude de l'alimentation électrique des inducteurs et le calcul des valeurs des courants qu'il faut injecter. Il est nécessaire de déterminer la matrice d'impédances du système pour chaque configuration de chauffage envisagé, afin de refléter l'état magnétique et électrique du système pour une géométrie donnée. La dimension N de la matrice est donnée par le nombre d'inducteurs, ici N=3. La matrice d'impédances doit être complète pour tenir compte de tous les effets de couplage. La détermination de cette matrice pouvant être complexe, plusieurs moyens analytiques ou numériques, ou des mesures en ligne et en continu par injection de signaux particuliers, peuvent être mis en oeuvre. Ainsi modélisée, l'équation générale du système peut s'écrire : V=Z.I V : Tensions sinusoïdales aux bornes des inducteurs ; : Courants dans les enroulements des inducteurs; Z : Matrice d'impédances du système. Dans notre cas, la matrice Zpeut s'écrire sous la forme : 2951606 -7- z11(w) z12 (w) z13 (w) z21 (w) z22 (w) z23 (w) Z31 (w) z32 (w) z33 (w) ou encore : R11+jL11w R12+jL12w R13+i-13w R21 + jL21w R22 + jL22w R23 + jL23w R31 + jL31 w R32 + I L32 CO R33 + jL33 CO Umm : représente l'inductance propre de chaque inducteur ; Lmn = Lnm : représente les inductances mutuelles entre inducteurs ; Rmm : représente les résistances propres de chaque inducteur ; Rmn ù Rem : représente les résistances équivalentes dues aux courants induits. Avec la connaissance des relations électromagnétiques entre les bobines et la pièce à chauffer, il est possible de procéder au calcul des courants à injecter dans chacune des bobines afin d'obtenir le chauffage désiré. Il est à noter que différentes configurations ou méthodes classiques de calcul essaient de minimiser les termes de couplage non-diagonaux afin de s'affranchir des problèmes liés aux interactions entre les bobines. De plus, pour de nombreux cas où les couplages sont faibles, les résistances propres de chaque inducteur sont souvent grandes devant les résistances équivalentes dues aux courants induits. Les méthodes classiques utilisent ainsi une matrice simplifiée, c'est-à-dire non complète, qui conserve seulement les termes diagonaux. Ceci implique une régulation simplifiée du chauffage, mais au détriment du contrôle précis du profil de température et de la flexibilité de l'installation, en particulier dans la zone située sous les bobines. Au contraire, la présente invention prend en compte la matrice d'impédances complète du système afin d'améliorer la détermination des courants à injecter dans les bobines et donc améliorer le contrôle du profil de température de la pièce chauffée. Dans l'exemple décrit, nous avons trois inducteurs alimentés par trois sources de courants différentes. La détermination des courants à injecter dans chaque z= Z= -8- bobine revient à déterminer cinq variables inconnues, la phase du courant dans l'inducteur Indl servant de référence et n'étant donc pas une inconnue. En effet, pour une tôle donnée constituant la pièce à chauffer, les inconnues sont : • 11 : Valeur efficace du courant dans l'inducteur Indl, lequel courant est pris comme référence de phase ; • /2 et 92 : Valeur efficace du courant dans l'inducteur Ind2, et déphasage de ce courant par rapport à 11 ; • /3 et 93 : Valeur efficace du courant dans l'inducteur Ind3, et déphasage de ce courant par rapport à 11 . Le vecteur des inconnues peut alors s'écrire : r X =~1 3 /2 3 q)2 3 13 3 q)3 } (1) Il n'est pas possible de déterminer facilement ces inconnues par les méthodes de résolution habituelles. En effet, à l'exception de cas très simples, la formulation analytique liant les données géométriques, les courants électriques dans les inducteurs, la répartition spatiale du champ électromagnétique et la densité de puissance en tous points est quasi impossible avec autant de variables. Les logiciels classiques de calcul de champ basés sur des techniques numériques de découpage du domaine d'études en mailles élémentaires permettent de connaître la répartition du champ magnétique, et par conséquent de calculer les densités de puissance dans les pièces conductrices en fonction des courants injectés dans les inducteurs. Dans notre cas, une problématique inverse se pose, puisqu'il s'agit de savoir s'il existe une ou plusieurs valeurs du vecteur x permettant d'obtenir un profil de densité de puissance voulu dans la pièce. Par l'application de l'équation de la chaleur, il est bien connu que la densité de puissance injectée Dp dans une pièce conductrice donne une bonne image du comportement thermique du produit chauffé. Par exemple, dans le cas d'une chauffe statique où la vitesse de déplacement du matériau traité est nulle, la connaissance de la température instantanée T du matériau traité nécessite -9- classiquement la résolution temporelle d'une forme simplifiée de l'équation de la chaleur: p.Cp âT =div(À.gradT)+Dp : représente la masse volumique ; : représente la capacité thermique massique ; À : représente la conductivité thermique. La résolution de cette équation implique une intégration en temps réel, ce qui n'est pas très difficile. De plus, dans le cas d'un chauffage « flash », c'est-à-dire si le temps de chauffe est petit de telle sorte que l'on peut négliger la diffusion thermique de la chaleur au sein du matériau pendant cette durée, l'expression se simplifie encore de la sorte : pCpat=Dp (2) Nous obtenons donc une expression simplifiée classique permettant de relier la densité de puissance injectée Dp et l'élévation de la température. Ainsi, à partir du profil thermique souhaité pour la pièce chauffée, on obtient le profil de densité de puissance recherché. Dans l'exemple en référence à la figure 1, le système est invariant suivant l'axe de révolution du disque en tôle et dans l'épaisseur de la tôle. On prend donc en compte une seule dimension du disque, à savoir la direction radiale de la zone considérée du disque. Pour la détermination du vecteur x des inconnues, on sait que la densité de puissance suivant le rayon de la zone considérée se calcule par l'équation suivante : Dp(r, x) = , soit : Dp(r, x) = 1 (JR (r, x) + J; (r, x)) (3) Q o où 6 représente la conductivité électrique, J représente le vecteur densité de courant défini sur le rayon r dans la pièce, JR (r,x) et JI (r,x) représentant les composantes réelle et imaginaire de ce vecteur en fonction du rayon de la zone considérée. P Cp - 10 - Le système pris en exemple est complètement linéaire, c'est-à-dire en particulier sans matériaux ferromagnétiques ni hystérésis. Nous pouvons donc appliquer le théorème de superposition des sources pour chacune des alimentations des trois inducteurs. On notera qu'un principe similaire peut être mis en oeuvre dans un système non linéaire. Nous obtenons ainsi des fonctions images des densités de courant en fonction du rayon r de la zone annulaire considérée du disque chauffé, chaque fonction image fk étant représentative de la relation liant la densité de courant Jk(r), créée par un inducteur, au courant ik alimentant cet inducteur. Ces fonctions images sont vectorielles et ont des composantes réelle et imaginaire définies de la manière suivante : fkR(r)= JkR(r) 'k fkl(r)= Jkl(r) Ik Au final, dans notre exemple à trois inducteurs, le calcul vectoriel de la densité totale de courant induit dans la zone annulaire de rayon r du disque peut s'exprimer ainsi : 3 J(r, x) = 1 (fkR (r) + jfkl (r)),Ik k=1 ,avec j2=-1, soit: 3 (r, x) rfkR(rr) l+ rr)).( llIkR + ) = ( Ifk, ( ~Ikl ) k=1 d'où Ce qui peut s'écrire aussi J(r, x) = JR (r, x) + p, (r, x) (4) On obtient donc une relation entre le vecteur densité de courant induit dans la zone considérée de la pièce et les vecteurs des courants dans les inducteurs. Avec d'une part la matrice d'impédances liant les grandeurs électriques entre inducteurs, et d'autre part les fonctions images des densités de courant dans la -11- pièce, nous disposons ainsi de toutes les informations nécessaires au calcul du vecteur des inconnues x à partir d'un profil de densité de puissance déterminé. On notera que l'on peut aussi faire intervenir dans ce calcul le vecteur des condensateurs, c'est-à-dire le vecteur des capacités des circuits oscillants, puisque ces capacités ne sont généralement pas rigoureusement égales du fait des tolérances de fabrication et qu'elles peuvent en outre dériver quelque peu. Pour le calcul, on pourra utiliser des logiciels de résolution des équations aux dérivées partielles, avec diverses techniques numériques possibles telles que les éléments finis, les différences finies, les volumes finis, les intégrales de frontière, les éléments de circuits partiels, ou tout autre technique du même genre. Cette méthode a été décrite pour un exemple donné d'un système magnétiquement couplé relativement simple, mais est néanmoins transposable à tout système plus complexe et non symétrique. Le nombre de bobines n'est pas limité, et diverses formes et configurations des bobines ou des pièces à chauffer sont envisageables, comme dans les exemples visibles sur les figures 3 à 6. Une fois que la fonction image de la densité de courant est déterminée, la fonction image de la densité de puissance DD(r, x) est déterminée par les relations des équations (3) et (4) ci-dessus. Il est avantageux en outre d'optimiser par le calcul le vecteur des inconnues x. Le problème d'optimisation consiste à calculer un vecteur optimisé x permettant de minimiser la différence entre la fonction image de la densité de puissance et une fonction densité de puissance de référence D pre' (r) qui correspond à un profil de densité de puissance de référence que l'on cherche à injecter dans le disque métallique. Cette fonction densité de puissance de référence prend par exemple une valeur constante si nous recherchons une homogénéité de température sur le disque. Il est cependant envisageable d'avoir une fonction non constante afin d'obtenir des profils particuliers de chauffe. Avec l'appareillage de la figure 1, la demanderesse a réalisé des essais avec différentes fonctions densités de puissances de références correspondant par exemple à des profils sinusoïdaux ou triangulaires dans la direction radiale du disque, et les résultats sont très satisfaisants. - 12 - L'optimisation consiste donc à minimiser la fonction g(r, x) = Dp(r, x) ù Dpre' (r tout en fixant des limites hautes et basses X" et xi' sur les inconnues recherchées. Cela nous permet d'éliminer entre autres les solutions aberrantes ou qui n'ont pas de réalité physique. La formulation du r problème d'optimisation revient donc à minimiser g (r, X) avec x xä } et HI x;E x;8,x; J i=1,...,n Après résolution du problème, nous obtenons un vecteur optimisé x contenant toutes les amplitudes des vecteurs des courants dans les inducteurs et leurs phases respectives, pour le disque métallique donné. Un des résultats pour un exemple de disque de 650 mm de diamètre, avec une référence de densité de puissance DpNef égale à lOMW/m3, donne un écart relatif maximal de 3% sur la fonction image de la densité de puissance Dp(r, x) comme représenté sur la figure 7. Cette méthode de résolution peut facilement être élargie afin de prendre en compte plusieurs dimensions d'un disque, par exemple trois si outre le rayon on tient compte de la position angulaire et de l'épaisseur de matière de la zone considérée, tout en tenant compte également de l'égalité de la compensation réactive nécessaire aux bornes de chaque bobine pour que les trois circuits oscillants oscillent à des fréquences très voisines. Nous passerions ainsi d'un vecteur à cinq inconnues à un vecteur à dix-huit inconnues, sans changer le système physique. La méthode expliquée ci-dessus pour la détermination du vecteur optimisé x est avantageusement utilisée dans le procédé de chauffage par induction selon l'invention, ce procédé pouvant être mis en oeuvre en particulier dans l'un ou l'autre des dispositifs de chauffage représentés sur les figures 8 et 9. Sur la figure 8 est représenté schématiquement un premier mode de réalisation d'un dispositif de chauffage par induction selon l'invention, dans lequel l'alimentation 1 des onduleurs est une source de courant. - 13 - Le dispositif de chauffage comprend des inducteurs Indl, Ind2,..., Indp, couplés magnétiquement, chaque inducteur étant alimenté par un onduleur de courant 01, 02,..., Op, qui lui est propre et associé à un condensateur CI, C2,..., Cp, pour former un circuit oscillant OC1, 0C2,..., OCp. Les onduleurs de courant sont mis en série avec l'alimentation 1. Chaque onduleur comprend généralement des interrupteurs électroniques bidirectionnels, et est commandé par une unité de commande aussi appelée modulateur Ml, M2,..., Mp. De façon connue en soi, chaque modulateur conçoit des ordres de commande des interrupteurs sous la forme d'impulsions, et le décalage dans le temps de ces ordres permet de faire varier l'amplitude AI, A2,..., Ap, et la phase (pi, (p2,..., (pp, du courant Ip, traversant l'inducteur correspondant. Les circuits oscillants possédent au moins approximativement la même fréquence de résonance, ce qui permet de maximiser l'efficacité de l'induction puisque les inducteurs travaillent sensiblement à cette fréquence, et permet aussi de diminuer les pertes dans les onduleurs. Des moyens de détermination des paramètres d'amplitude et de phase des courants Ip, des inducteurs, non représentés sur la figure, sont prévus pour fournir ces paramètres à des unités comparatrices Ci, C2,..., C. Ces moyens de détermination peuvent consister par exemple en des transformateurs de courant disposés chacun en série avec un inducteur, mais d'autres moyens sont envisageables. On pourrait par exemple mesurer le courant actif fourni par l'onduleur au circuit oscillant, et calculer le courant dans l'inducteur à l'aide des paramètres d'inductance et de capacité. On prévoit en outre des moyens de détermination d'un profil de température effectif de la pièce métallique chauffée 10, non représentés sur la figure, en disposant par exemple des thermocouples sur un nombre n de zones chauffées et en relevant les températures Ai mes, 02 mes, ..., An mes, mesurées. On peut aussi déterminer ces températures en utilisant une caméra thermique, ou encore procéder par des calculs à partir des courants induits si par exemple des zones chauffées sont trop confinées pour une mesure directe. Le profil de température effectif est par exemple déterminé en continu pendant la chauffe et est régulièrement comparé à un profil de température de 2951606 - 14 - référence 01 ref, 02 ref, ..., 8n ref, correspondant au profil final de chauffe souhaité pour la pièce et préalablement entré dans une mémoire. Cette comparaison est effectuée par un comparateur 2, qui peut intégrer ladite mémoire. Le résultat est traité par un calculateur qui, à partir d'une équation déduite de l'équation de la chaleur et éventuellement simplifiée comme l'équation (2) précédente, calcule le profil de densité de puissance de référence Del, Dpre 2 Dpre n que le dispositif de chauffage doit injecter dans la pièce pour atteindre le profil de température de référence. Le calculateur peut consister en une mémoire dans laquelle est rentré un tableau de profils précalculés de densité de puissance de référence correspondants à différents profils de température effectifs pour une ou plusieurs configurations de pièces et un ou plusieurs profils de densité de puissance de référence. Un calculateur établit des courants cibles que doivent délivrer les onduleurs afin que les courants des inducteurs atteignent des valeurs cibles appropriées Ii ref, 12 ref,..., Ip ref, pour injecter dans la pièce le profil de densité de puissance de référence. Ce calcul utilise la matrice d'impédances Z avec les fonctions images vectorielles fk et préférablement le vecteur des capacités des circuits oscillants, définis précédemment. Les unités comparatrices El, g2,..., Ep comparent les paramètres des courants mesurés ou calculés Ii mes, 12 mes,..., Ip mes, des inducteurs aux valeurs cibles Ii ref, 12 ref,..., Ipref, et déterminent les écarts de courants 8h corr, 812 corn..., 8Ip corn à corriger, aussi appelés courants de correction. Des unités de traitement CORR1, CORR2,..., CORRp, des paramètres d'amplitude et de phase ce ces courants de correction génèrent des instructions de correction envoyées aux modulateurs pour commander les onduleurs de façon à corriger les courants traversant les inducteurs. Avantageusement, on peut déterminer une valeur initiale Z;,,; de la matrice d'impédances Z pour une température moyenne initiale Oin, donnée des inducteurs et de la pièce à chauffer, puis déterminer à intervalles variables ou périodiques la matrice d'impédances modifiée Zmod(0) pour au moins une valeur augmentée Omod de la température moyenne 0, et on utilise la matrice d'impédances modifiée pour recalculer les courants cibles. Dans le cas d'intervalles variables d'échantillonnage, le calcul des courants cibles peut s'effectuer à chaque fois que - 15 - la température moyenne A mesurée atteint sensiblement une nouvelle valeur augmentée Omod parmi une série de valeurs prédéterminées. Avantageusement, l'onduleur de courant alimentant l'inducteur de plus faible impédance, par exemple la bobine Indl dans l'exemple de la figure 1, est choisi comme onduleur de référence commandé en pleine onde, puisque le courant dans cet inducteur, plus grand que celui dans les autres inducteurs, est pris préférablement comme référence de phase. Sur la figure 8, le courant Ii étant aussi pris comme référence de phase, il est avantageux que l'unité comparatrice correspondante Ei reçoive les paramètres du courant le mes délivré par l'alimentation 1. De cette façon, l'unité de traitement CORR1 associée sera adaptée pour générer des instructions de régulation envoyées à l'alimentation 1 via un modulateur de pilotage M'l, de façon à modifier le courant délivré par l'onduleur 01 au circuit oscillant OC1, ce qui permet de contrôler l'amplitude de ce courant et donc de modifier l'amplitude du courant Ii dans l'inducteur Indl. Pour chauffer une pièce métallique avec le dispositif de chauffage décrit ci-dessus, on utilise le procédé comprenant les étapes suivantes: a) on compare le profil de température effectif de la pièce au profil prédéterminé de température de référence, et on calcule le profil de densité de puissance de référence que le dispositif doit injecter dans la pièce pour atteindre le profil de température de référence ; b) à partir de la matrice d'impédances Z du système, préférablement associée au vecteur des capacités des circuits oscillants, et par la connaissance des fonctions images vectorielles fk, on calcule les courants cibles que doivent délivrer les onduleurs afin que les courants des inducteurs atteignent les valeurs cibles appropriées pour injecter dans la pièce le profil de densité de puissance de référence ; c) on détermine par mesure ou par calcul les courants traversant les inducteurs pour les comparer aux valeurs cibles de ces courants et déterminer les écarts de courants à corriger, et on envoie aux modulateurs les instructions de correction afin de commander les onduleurs de façon à corriger les courants. - 16 - Avantageusement, après avoir effectué successivement les étapes (a) et (b), on effectue au moins une fois l'étape (c) pour diminuer les écarts de courants à corriger, puis on réitère au moins une fois les étapes (a), (b) et (c) en actualisant le profil de température effectif par des mesures de température en différentes zones chauffées de la pièce. Sur la figure 9 est représenté schématiquement un deuxième mode de réalisation d'un dispositif de chauffage par induction selon l'invention, dans lequel l'alimentation 1 des onduleurs est une source de tension. Le dispositif de chauffage est analogue à celui du premier mode de réalisation de la figure 8, mais les onduleurs de courant sont mis en parallèle avec la source de tension. Ce mode de réalisation présente certains avantages, notamment celui de diminuer les pertes par conduction dans les onduleurs. Par contre, le paramètre de courant le caic représentatif du courant que délivre l'alimentation 1 à l'onduleur 01 doit être calculé à partir de la tension d'alimentation à l'aide d'une matrice d'impédances Z'. REVENDICATIONS1. Procédé de chauffage par induction mis en oeuvre dans un dispositif de chauffage d'une pièce métallique, le dispositif comprenant des inducteurs (Indl, Ind2,..., Indp) couplés magnétiquement, chaque inducteur étant alimenté par un onduleur (01, 02,..., Op) qui lui est propre et associé à un condensateur (Cl, C2,..., Cp) pour former un circuit oscillant (OC1, 0C2,..., OCp), lesdits circuits oscillants possédant au moins approximativement la même fréquence de résonance, chaque onduleur étant commandé par une unité de commande (Ml, M2,..., Mp) de façon à faire varier l'amplitude (Ai, A2,..., Ap) et la phase (spi, (p2,..., (pp) du courant (Ii, I2,..., Ip) traversant l'inducteur correspondant, le dispositif comprenant en outre des moyens de détermination dudit courant (Ii, I2,..., Ip) ainsi que des moyens de détermination d'un profil de température effectif (A1 mes, 02 mes, ..., On mes) de ladite pièce métallique, ledit procédé comprenant les étapes suivantes : a) on compare ledit profil de température effectif (Ai mes, 02 mes, ..., An mes) à un profil de température de référence (A1 ref, 02 ref, ..., Onref), et on calcule un profil de densité de puissance de référence (Dprefl Dpre 2 Dpre n) que le dispositif de chauffage doit injecter dans ladite pièce pour atteindre ledit profil de température de référence ; b) à partir d'une matrice d'impédances (Z) déterminée par la connaissance des relations électromagnétiques liant lesdits inducteurs entre eux et à ladite pièce et par la connaissance de fonctions images vectorielles (fk) représentatives des relations liant les densités de courant créées par les inducteurs aux courants (Il, I2,..., Ip) traversant les inducteurs, on calcule des courants cibles que doivent délivrer les onduleurs afin que les courants des inducteurs atteignent des valeurs cibles (Il ref, 12 ref,..., Ip ref) appropriées pour injecter dans ladite pièce ledit profil de densité de puissance de référence (Dprefl Dpre 2 Dpren ) ; c) on détermine les courants (Il mes, 12 mes,..., Ip mes) traversant les inducteurs pour les comparer aux dites valeurs cibles (Il ref, 12 ref,..., Ip ref) et déterminer des écarts de courants (8I1 corn, 812 corn..., 8Ip corn) à corriger, et on envoie aux dites unités de commande (Ml, M2,..., Mp) des instructions de- 18 - correction en fonction desdits écarts de courants afin de commander les onduleurs de façon à corriger les courants traversant les inducteurs. 2. Procédé de chauffage selon la revendication 1, dans lequel on détermine les capacités desdits condensateurs (Cl, C2,..., Cp) et on associe ladite matrice d'impédances (Z) à un vecteur (C) des capacités. 3. Procédé de chauffage selon la revendication 1 ou 2, dans lequel on détermine une valeur initiale (Z;,,;) de ladite matrice d'impédances (Z) pour une température moyenne initiale (Oin,) donnée desdits inducteurs et de ladite pièce, puis on détermine à intervalles variables ou périodiques la matrice d'impédances modifiée (Zmod(0)) pour au moins une valeur augmentée (Omod) de ladite température moyenne, et on utilise ladite matrice d'impédances modifiée pour recalculer les dites valeurs cibles (Ii ref, 12 ref,..., Ip ref). 4. Procédé de chauffage selon l'une quelconque des revendications 1 à 3, dans lequel après avoir effectué successivement les étapes (a) et (b) on effectue au moins une fois l'étape (c) pour diminuer lesdits écarts de courants (8Ii corn, 8I2 aorr,..., 8Ip eorr) à corriger, puis on réitère au moins une fois les étapes (a), (b) et (c) en actualisant ledit profil de température effectif (Ai mes, 02 mes, ..., 8n mes) par des mesures de température en différentes zones chauffées de la pièce. 5. Procédé de chauffage selon l'une quelconque des revendications 1 à 4, dans lequel pour la détermination par calcul des dites valeurs cibles (Ii ref, 12 ref,..., Ip ref) à l'étape (b), grâce à la connaissance desdites fonctions images vectorielles (fk) on calcule des fonctions images des densités de puissance (Dp(r,x)) selon les caractéristiques spatiales (r) des zones de la pièce dans lesquelles lesdites densités de puissance sont injectées, et on calcule un vecteur optimisé (x) des courants cibles à déterminer en minimisant la différence entre chacune desdites fonctions images des densités de puissance (Dp(r,x)) et une fonction densité de puissance de référence (DpYef(r)) correspondant audit profil de densité de puissance de référence (Dprefi Dpre i ..., Dpren, ). 2951606 - 19 - 6. Dispositif de chauffage par induction comprenant : des inducteurs (Indl, Ind2,..., Indp) couplés magnétiquement, chaque inducteur étant associé à un condensateur (Cl, C2,..., Cp) pour former un circuit oscillant (OC1, 0C2,..., OCp), lesdits circuits oscillants possédant au moins approximativement la même fréquence de résonance ; des onduleurs (01, 02,..., Op) alimentant chacun un inducteur (Indl, Ind2,..., Indp) qui lui est propre, chaque onduleur étant commandé par une unité de commande (Ml, M2,..., Mp) de façon à faire varier l'amplitude (Ai, A2,..., Ap) et la phase (spi, (p2,..., (pp) du courant (Ii, I2,..., Ip) traversant l'inducteur correspondant ; caractérisé en ce qu'il comprend en outre : des moyens de détermination des courants (Ii, I2,..., Ip) traversant les inducteurs ainsi que des moyens de détermination d'un profil de température effectif (Ai mes, 02 mes, ..., an mes) d'une pièce métallique chauffée par le dispositif ; des moyens de comparaison dudit profil de température effectif (Ai mes, 02 mes, ..., an mes) par rapport à un profil de température de référence (Ai ref, 02 ref, ..., en ref) ; des moyens de calcul d'un profil de densité de puissance de référence (Dprefi, Dpre 2, , Dpre n) que le dispositif de chauffage doit injecter dans ladite pièce pour atteindre ledit profil de température de référence ; des moyens de calcul, basés sur la connaissance d'une matrice des impédances (Z), de courants cibles que doivent délivrer les onduleurs afin que les courants des inducteurs atteignent des valeurs cibles (Ii ref, 12 ref,..., Ip ref) appropriées pour injecter dans ladite pièce ledit profil de densité de puissance de référence (Dprefi Dpre 2 Dpren ) ; des moyens de comparaison (Ei, E2,..., Ep) des courants (Ii mes, 12 mes,..., Ip 111es) traversant les inducteurs par rapport aux dites valeurs cibles (Ii ref, 12 ref,..., Ip ref), aptes à déterminer des écarts de courants (8Ii cOrr, 8I2 corn..., 8Ip cOrr) à corriger, et des moyens de traitement (CORR1, CORR2,..., CORRp) desdits écarts de courants aptes à générer des instructions de correction envoyées aux dites- 20 - unités de commande (Ml, M2,..., Mp) pour commander les onduleurs de façon à corriger les courants traversant les inducteurs. 7. Dispositif de chauffage par induction selon la revendication 6, dans lequel les onduleurs (01, 02,..., Op) sont alimentés par une même alimentation (1) source de courant ou source de tension, et dans lequel lesdits moyens de comparaison desdits courants déterminés (Ii mes, I2 mes,..., Ip mes) traversant les inducteurs comprennent des unités comparatrices (Ci, g 2,ù, g p) recevant chacune des paramètres déterminés (Ai, (pi ; A2, (p2 ;... ; Ap , (pp) d'un courant traversant un inducteur (Ii mes, I2 mes,..., Ip mes) et des paramètres des valeurs cibles correspondantes (Ii ref, 12 ref,..., Ip ref) et étant chacune reliée à une unité de traitement (CORR1, CORR2,..., CORRp) desdits écarts de courants, une (Ci) desdites unités comparatrices recevant en outre des paramètres (Ie mes, le cala) représentatifs de ce que délivre ladite alimentation (1) et son unité de traitement associée (CORR1) étant adaptée pour générer des instructions de régulation envoyées à ladite alimentation (1) de façon à modifier le courant ou la tension qu'elle délivre.