FR 3046018 A1 20170623 FR 1562824 A 20151218 1 PROCEDE D'OPTIMISATION DE CHAUFFAGE PAR INDUCTION. La présente invention est relative à un procédé d'optimisation de chauffage par induction et notamment à un procédé d'optimisation de chauffage par induction mis en oeuvre dans un dispositif comprenant des inducteurs couplés magnétiquement. Par couplage magnétique, on entend que les inducteurs produisent entre eux des inductions mutuelles. Plus particulièrement, l'invention se rapporte à un procédé d'optimisation de chauffage par induction mis en oeuvre dans un dispositif de chauffage d'une pièce, mais également à un dispositif de chauffage par induction et à un ensemble de chauffage. Il est connu de l'art antérieur différentes solutions visant à chauffer par induction une pièce. Les techniques les plus classiques de chauffage par induction connues mettent en oeuvre des configurations qui sont satisfaisantes lorsque les pièces à chauffer sont toujours de même nature et de mêmes dimensions. Toutefois, de telles techniques de chauffage ne permettent pas de contrôler de façon fine des profils de chauffage sur des pièces ayant des formes et des dimensions différentes et/ou complexes et/ou variables au cours 20 du temps 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 25 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 30 délicat, sans compter qu'il peut avoir des répercussions néfastes sur les générateurs de fréquence, par exemple une casse de composants. Une solution pour pallier ces problèmes est par exemple de minimiser les interactions magnétiques mutuelles entre les inducteurs, soit en 3046018 2 les éloignant soit en les blindant avec des culasses magnétiques pour canaliser le flux soit en les occultant avec des écrans magnétiques dans lesquels vont se développer des courants induits. Toutefois, cette solution ne permet pas d'obtenir un chauffage de pièces avec une précision et une finesse importante. 5 Elle nécessite de plus l'ajout de pièces supplémentaires volumineuses et coûteuses (culasses et écrans) ainsi que les dispositifs de motorisation et de contrôle pour assurer leur déplacement afin d'adapter le profil de densité de puissance à la pièce à chauffer. Le contrôle de la répartition du champ magnétique et de la densité 10 de puissance injectée dans la pièce peuvent ainsi par exemple être réalisés par l'utilisation d'actionneurs mécaniques. Une utilisation complexe de formes d'inducteurs peut également être envisagée pour effectuer ce contrôle de la répartition du champ magnétique et de la densité de puissance injectée dans une pièce à chauffer par induction. Toutefois, de tels moyens entrainent un 15 temps de réponse particulièrement long, notamment lorsque, dans le processus industriel, un changement des formes des pièces et/ou des dimensions des pièces à chauffer est réalisé. Ainsi, de telles lignes de production ne permettent pas de s'adapter de façon rapide à un fonctionnement continu ni à un changement de position ou de format des pièces à chauffer, et d'adapter en fonction de ce changement le profil de température désiré. De plus, des pertes supplémentaires sont induites dans les organes mécaniques qui peuvent s'échauffer et les forces électromagnétiques en présence peuvent provoquer des vibrations sur ces pièces. La publication W02011048316 de la demanderesse décrit un exemple d'un procédé tel que précédemment décrit, qui remédie aux désavantages précités. Plus précisément, cette publication divulgue un procédé de chauffage dans lequel le chauffage par induction est mis en oeuvre dans un dispositif comprenant des inducteurs couplés magnétiquement et qui permet de s'affranchir de toutes pièces mécaniques de réglage de sorte à éviter les pertes induites dans ces pièces mécaniques et de sorte à simplifier la mise en oeuvre du chauffage de la pièce. En l'espèce le procédé décrit dans la publication W02011048316 a permis de prendre en compte les nombreux couplages, d'une part entre les différents inducteurs et d'autre part entre les inducteurs et 3046018 3 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. Le procédé décrit dans W02011048316 permet également l'ajustement du chauffage à différents profils de température recherchés en temps réel, en agissant sur la commande 5 d'onduleurs alimentant les inducteurs et sans avoir besoin d'ajuster la structure des inducteurs. Le procédé décrit dans la demande W02011048316 apporte satisfaction. Toutefois, dans un souci de s'adapter à de nouvelles configurations d'inducteurs qui conduisent à des couplages plus forts entre la charge et les inducteurs d'une part et afin d'optimiser l'efficacité énergétique 10 d'autre part, il existe le besoin de développer un nouvel outil qui permet une meilleure adaptabilité. En effet, des couplages plus forts entre la charge et les inducteurs peuvent entrainer des fortes perturbations des courants dans les inducteurs alimentés avec des onduleurs de tension ou de courant. Ces fortes perturbations viennent perturber la commande des onduleurs et des inducteurs 15 et par conséquence la qualité du chauffage de la pièce. La présente invention a notamment pour but de pallier ces inconvénients. A cet effet, l'invention concerne un procédé d'optimisation de chauffage par induction mis en oeuvre dans un dispositif de chauffage d'une 20 pièce, le dispositif comprenant des inducteurs couplés magnétiquement, chaque inducteur étant alimenté par un onduleur de tension ou de courant qui lui est propre et associé à un condensateur pour former un circuit oscillant, lesdits circuits oscillants possédant sensiblement la même fréquence de résonance, chaque onduleur étant commandé par une unité de commande de 25 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 chauffée, ledit procédé comprenant les étapes suivantes : 30 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 ; 3046018 4 b. on calcule des courants cibles appropriées pour injecter dans ladite pièce ledit profil de densité de puissance de référence par : - la connaissance de fonctions images vectorielles représentatives des relations liant les densités de courant créées par les 5 inducteurs aux courants traversant les inducteurs de sorte à déterminer les valeurs de l'amplitude des courants, et de la phase des courants des inducteurs, - la connaissance d'une matrice d'impédances déterminée par la connaissance des relations électromagnétiques liant lesdits 10 inducteurs entre eux et à ladite pièce, - la valeur des condensateurs des circuits oscillants, - l'identification et la détermination des amplitudes et de la phase des harmoniques du courant délivré par chaque onduleur de sorte à contrôler le fondamental et les harmoniques du courant dans 15 l'alimentation d'un inducteur par son onduleur associé, 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 20 onduleurs de façon à corriger les courants traversant les inducteurs pour la fréquence fondamentale ainsi que pour les harmoniques. Grâce à ces dispositions, on vient commander les harmoniques des courants aux bornes des inducteurs de sorte à améliorer les performances des inducteurs. En effet, la maîtrise des harmoniques du courant aux bornes des 25 inducteurs permet d'envoyer aux bornes des inducteurs un signal moins perturbé et donc avec moins de distorsions. En outre, le procédé selon l'invention permet une plus grande adaptabilité. Dans des modes de réalisation préférés de l'invention, on peut éventuellement avoir recours en outre à l'une et/ou à l'autre des dispositions 30 suivantes. Selon une réalisation, les harmoniques du courant délivré par chaque onduleur sont optimisés de sorte à améliorer la qualité du chauffage par induction. 3046018 5 Selon une réalisation, on calcule la capacité électrique des condensateurs pour compenser l'énergie réactive portée par les fondamentaux des courants des inducteurs en fonction de la connaissance des relations électromagnétiques liant lesdits inducteurs entre eux. 5 Selon une réalisation, lors de l'étape b), on utilise en outre des fonctions images vectorielles représentatives des relations liant les densités de courant créées par les inducteurs aux courants traversant les inducteurs pour chaque harmonique considéré du courant de sorte à obtenir les valeurs des amplitudes du courant et de la phase des courants des inducteurs. 10 Selon une réalisation, lors de l'étape b) les valeurs cibles sont calculées par la connaissance de la matrice d'impédances déterminée par la connaissance des relations électromagnétiques liant lesdits inducteurs entre eux et à ladite pièce pour le fondamental du courant et la connaissance des matrices d'impédances pour chaque harmonique h considéré. 15 Selon une réalisation, on identifie et on détermine les amplitudes et les phases des harmoniques de rang i, avec i allant de 2 à n, le nombre n étant un entier naturel à déterminer. Par exemple, les harmoniques de rang 2, 3, 4, 5, 6, 7, ..., 10 peuvent être considérées. Dans des variantes de réalisation, les harmoniques de rang 3, 5 et 7 peuvent être considérées. . 20 Selon une réalisation, on détermine les capacités des condensateurs et on associe une matrice d'impédance à un vecteur des capacités. Selon une réalisation, une modulation de largeur d'impulsions est réalisée par insertion d'angles supplémentaires afin de supprimer les harmoniques de rang supérieur ou égal à deux des courants traversant les 25 inducteurs. Selon une réalisation, des harmoniques de rang supérieur à deux des courants traversant les inducteurs sont supprimés. Selon une réalisation, le fondamental du courant délivré par chaque onduleur est contrôlé. 3 0 Selon une réalisation, un signal carré peut être utilisé. La présente invention concerne également un dispositif de chauffage par induction comprenant : - des inducteurs couplés magnétiquement, chaque inducteur étant associé à un 3046018 6 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 5 l'amplitude et la phase du fondamental et des harmoniques du courant traversant l'inducteur correspondant ; et comprenant en outre : - des moyens de détermination du courant fondamental et des harmoniques des courants traversant les inducteurs ainsi que des moyens de détermination 10 d'un profil de température effectif d'une pièce 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 15 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 et leurs harmoniques atteignent des valeurs cibles appropriées pour injecter dans ladite pièce ledit profil de densité de puissance 2 0 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 25 pour commander les onduleurs de façon à corriger les courants traversant les inducteurs. Selon une réalisation, les inducteurs présentent des culasses magnétiques de canalisation du flux magnétique. La présente invention concerne enfin un ensemble de chauffage par 30 induction comprenant un dispositif de chauffage par induction tel que précédemment décrit et une pièce chauffée par le dispositif de chauffage par l'intermédiaire du procédé de chauffage décrit plus haut. Selon une réalisation, la pièce chauffée est métallique. Selon une 3046018 7 réalisation, la pièce chauffée comprend un métal ferromagnétique. Toutefois, dans des variantes de réalisation, la pièce chauffée n'est pas forcément métallique. Par exemple la pièce chauffée peut être uniquement conductrice de l'électricité et par suite réceptive au chauffage par induction. 5 D'autres caractéristiques et avantages de l'invention apparaîtront au cours de la description suivante d'une de ses formes de réalisation, donnée à titre d'exemple non limitatif, en regard des dessins joints. Sur les dessins : 10 - les figures 1 A à lE représentent schématiquement différents exemples d'ensembles de chauffage par induction pour lesquelles le procédé de chauffage selon l'invention peut être mis en oeuvre ; - la figure 2 représente schématiquement un exemple d'inducteur comportant une culasse magnétique, 15 - la figure 3 représente schématiquement un premier mode de réalisation d'un dispositif de chauffage par induction mettant en oeuvre un procédé de chauffage, dans lequel les harmoniques des courants sont supprimés de sorte à réduire les perturbations du courant dans les inducteurs ; - la figure 4 représente un signal carré permettant l'élimination de 20 deux harmoniques et le contrôle du fondamental pour la mise en oeuvre du procédé selon le premier mode de réalisation - la figure 5 représente schématiquement un dispositif de chauffage par induction mettant en oeuvre un procédé de chauffage selon un deuxième mode de réalisation, dans lequel les harmoniques des courants sont utilisés 25 afin d'améliorer la qualité du chauffage de la pièce et d'augmenter le rendement global du système ; - la figure 6 représente un logigramme d'optimisation mis en oeuvre dans le deuxième mode de réalisation du procédé selon l'invention ; - la figure 7 représente schématiquement la superposition des effets 30 des harmoniques à partir des courants inducteurs dans le procédé d'optimisation selon le premier mode de réalisation - la figure 8 représente schématiquement la superposition des effets des harmoniques à partir des courants onduleurs dans le procédé 3046018 8 d'optimisation selon le deuxième mode de réalisation. Sur les différentes figures, les mêmes références désignent des éléments identiques ou similaires. Les figures 1A à 1E représentent des exemples non limitatifs de 5 dispositifs de chauffage dans lesquelles le procédé selon l'invention peut être mis en oeuvre. Par exemple, la figure 1A concerne une configuration de disque métallique 10 chauffé par flux transverse à l'aide de trois couples de bobines jumelles Ind1, Ind2, Ind3. La figure 1B représente un dispositif de chauffage par induction similaire à celui de la figure 1A pour le chauffage d'une tôle 11 qui 10 est mobile lors du procédé de chauffage par induction. La figure 1C représente schématiquement un deuxième exemple de dispositif de chauffage par induction, appliqué au chauffage d'une barre métallique 12 qui est mobile au cours du procédé de chauffage par induction. La figure 1D représente schématiquement un troisième exemple de dispositif de chauffage par 15 induction, appliqué au chauffage d'une tôle 13 qui est entrainée en translation par exemple lors du chauffage par induction. La figure 1E représente schématiquement un quatrième exemple de dispositif de chauffage par induction, appliqué au chauffage d'une tôle 13 qui peut être également entrainée en translation par exemple lors du chauffage par induction. 20 La pièce chauffée par le dispositif de chauffage par induction D est une pièce 10, 11, 12, 13 métallique ou non. Plus précisément la pièce chauffée est conductrice de l'électricité et donc réceptive au chauffage par induction. La pièce chauffée peut être, comme représenté sur les figures 1A à 1E une tôle comprenant une épaisseur sensiblement constante ou non et sensiblement 25 plane, une barre ou un tube s'étendant selon une direction longitudinale avec une section ouverte ou fermée, constante ou variable. Toutefois, dans des variantes de réalisation la pièce à chauffer peut présenter des formes plus complexes, par exemple avec des portions coudées. La pièce est chauffée par zones. 30 Comme représenté sur les figures 1A à 1 E, le dispositif de chauffage D comprend plusieurs inducteurs Ind1, Ind2, Ind3. Le dispositif de chauffage peut comprendre deux, trois, quatre, cinq, six ou plus inducteurs Ind1, Ind2, ..Indp. Dans la description qui suit, il est considéré que le dispositif de 3046018 9 chauffage comprend p inducteurs Ind1, .. Indp, avec p un nombre entier strictement supérieur à 1, par exemple avec p égale à 3. Chaque inducteur lind1, ...Indp est traversé par un courant li, ..lp (visible sur les figures 3 et 5). Le courant li, ..lp traversant les inducteurs lind1, 5 ...Indp peut être décomposé par séries de Fourier en un courant fondamental 11,1, _Io et une pluralité de courants harmoniques I1,h...10-,(avec h la fréquence harmonique ou l'ordre d'harmonique considéré). Les courants harmoniques sont dus au mode de commande des onduleurs (notamment aux courants carrés générés par les onduleurs) et/ou à la présence de dispositifs 10 magnétiques non linéaires (culasses magnétiques, pièces chauffées ferromagnétiques) dans le réseau du dispositif de chauffage. On parlera notamment d'harmoniques d'ordre h, avec h pouvant être un chiffre pair ou impair. Dans une variante de réalisation, on peut considérer les harmoniques à la fois paires et impaires. Dans une autre variante, seules les harmoniques 15 d'ordre impair peuvent être considérées. En général une harmonique d'ordre 1 correspond à un fondamental. Aussi, dans le reste de la description, sauf si spécifiquement mentionné, le terme « fondamental » correspond à l'harmonique d'ordre 1 et le terme « harmoniques » correspond à des harmoniques d'ordre strictement supérieur à 1. Tel que schématiquement 20 représenté sur les figures 3 et 5, chaque inducteur Ind1,...Indp peut être associé à un condensateur C1, ...Cp connecté en parallèle dans le cas d'onduleurs de courant. Toutefois, la présente invention peut aussi être mise en oeuvre avec des onduleurs de tension. Dans le cas d'onduleurs de tension, les condensateurs sont en série avec les inducteurs. 25 Chaque couple inducteur-condensateur forme un circuit oscillant OC1, ..00p. Les circuits oscillants OC1, ..00p peuvent être chacun sensiblement à la même fréquence de résonnance. Les circuits oscillants OC1, ..00p, et par conséquent les inducteurs Ind1, ...Indp sont alimentés par des onduleurs 01,...0p. Plus précisément, chaque circuit oscillant OC1, ..00p est 30 alimenté par un onduleur O1,...,Op correspondant. La pluralité d'onduleurs est en l'occurrence reliée à une unique source 1 continue. Les onduleurs 01,...,Op peuvent être soit des onduleurs de courant soit des onduleurs de tension. Chaque onduleur O1,...,Op peut comprendre des interrupteurs 3046018 10 électroniques bidirectionnels. Comme visible sur les figures 3 et 5, chaque onduleur O1,...,Op est commandé par une unité de commande, également appelée modulateur Mo-, avec p le nombre d'onduleurs et h le nombre d'harmoniques qui sont contrôlées. Plus précisément, à chaque onduleur peut 5 correspondre une pluralité de modulateurs. Ainsi, le premier onduleur 01 est commandé par une pluralité de modulateurs M1,1, - - -Mi,h. De façon similaire, le deuxième modulateur est commandé par une pluralité de modulateurs M2,1, ...M2,h. Il en est de même pour l'ensemble des onduleurs 01, ...Op du dispositif. Chaque modulateur conçoit des ordres de commande des 10 interrupteurs sous la forme d'impulsions, et le décalage dans le temps de ces ordres permet par exemple de faire varier l'amplitude A1,1, et la phase cp 1, 1, ...(p p,h du courant traversant l'inducteur correspondant pour le courant fondamental 111, ...1p,1 et pour les harmoniques du courant 11,h...1p,h que l'on souhaite considérer. La variation de l'amplitude du courant considérée en sortie 15 de chaque onduleur 01,...,Op s'effectue par modulation de largeur d'impulsions et notamment en introduisant un angle de décalage sur le signal généré par le modulateur commandant le courant considéré (fondamental du courant ou harmonique h du courant) de l'onduleur 01, ...Op. Le décalage peut aussi bien se faire en retard ou en avance par rapport au signal de commande de 20 l'onduleur d'un autre inducteur pris comme référence. Afin de déterminer et contrôler en temps réel la densité de puissance à injecter dans la pièce chauffée, des moyens de détermination des paramètres d'amplitude et de phase des courants traversant les inducteurs peuvent être prévus. Ainsi, il est possible de corriger la commande des onduleurs en temps 25 réel par les modulateurs M1,1, ...,Mp,h. Ces moyens de détermination des paramètres d'amplitude et de phase des courants traversant les inducteurs fournissent ces paramètres à des unités comparatrices , Ep. Les moyens de détermination des paramètres d'amplitude et de phase sont par exemple notamment des transformateurs de 30 courant disposés chacun en série avec un inducteur. Toutefois, dans des variantes de réalisation d'autres moyens peuvent être envisagés. En outre, des moyens de détermination d'un profil de température effectif 8. _ 1 mes, e2mes, - - - enmes de la pièce à chauffer sont prévus. Par exemples 3046018 11 des thermocouples peuvent être envisagés et disposés sur une pluralité n de zones chauffées et les températures 8. -lmes, e2mes, - - -enmes de chaque zone peuvent être ainsi mesurées. En variante, des caméras thermiques peuvent être utilisées. La mesure de température peut également être réalisée de façon 5 indirecte par calculs par exemple si les zones de la pièce où la température est à mesurer sont difficiles d'accès ou confinées. Le profil de température peut être déterminé en continu pendant la chauffe et être régulièrement ou de façon continue comparé à un profil de température de référence elref, e2ref, ---enref correspondant au profil final de 10 chauffe souhaité pour la pièce. Le profil final de chauffe souhaité pour la pièce à chauffer peut être préalablement enregistré par un manipulateur ou autre dans une mémoire. La comparaison entre le profil de chauffe de référence Dprefl , - - -Dprefn (aussi appelé profil de chauffe théorique) et le profil de chauffe mesuré Dmesl, 15 ... Dmesn est effectuée par un comparateur. Le résultat est traité par un calculateur. Le calculateur peut par exemple 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 20 puissance de référence. Le calculateur calcule, à partir des équations connues de la chaleur et ou à partir d'équations simplifiées, le profil de densité de puissance de référence Dpref1, ...Dprefn que le dispositif de chauffage doit injecter dans la pièce pour atteindre le profil de température de référence elref, e2ref, - - -enref- 25 A partir du profil de référence D prefl, - - -Dprefn que le dispositif de chauffage doit injecter dans la pièce pour atteindre le profil de température de référence elref, e2ref, ...enref, un calculateur établit des courants I cibles ---- ref1,1, ...lrefp,h que doivent délivrer les onduleurs afin que les courants aux bornes des inducteurs atteignent des valeurs cibles appropriés pour injecter dans la pièce à 30 chauffer le profil de densité de puissance de référence. Ces valeurs cibles appropriées Iref1,1, ---Irefp,h vont être déterminées en prenant en compte non seulement le courant fondamental (aussi appelé harmonique d'ordre un) du courant délivré par l'onduleur mais également les harmoniques du courant. En 3046018 12 effet, la prise en compte uniquement du fondamental du courant ne permet pas d'obtenir un profil de chauffe fin, à cause de fortes distorsions. Toutefois, la prise en compte des harmoniques du courant assure un réglage plus précis. Des unités comparatrices comparent les paramètres des courants mesurés 5 pour le fondamental et pour les harmoniques aux valeurs cibles Irefi,i, ...lrefp,h et détermine les écarts entre ces deux grandeurs. Des unités de traitement CORR1,1, ..., CORRp,h des paramètres d'amplitude et de phase de ces courants de correction (pour le fondamental de chacun des courants et pour les harmoniques considérés de chacun des courants) génèrent des instructions de 10 correction envoyées aux modulateurs pour commander les onduleurs de façon à corriger les amplitudes et le déphasage pour le fondamental et pour les harmoniques considérés des courants traversant les inducteurs. Plus particulièrement, la prise en compte des harmoniques du courant, ou encore la maitrise des harmoniques du courant, pour améliorer ou 15 affiner le profil de chauffe de la pièce à chauffer par le dispositif de chauffage peut être réalisée de deux manières. La première manière est représentée schématiquement sur la figure 3 tandis que la deuxième manière est représentée schématiquement sur la figure 5. De façon générale, les calculs utilisés et grandeurs mises en jeux 20 pour réaliser le procédé d'optimisation de chauffage par induction se basent notamment sur les équations connues et simplifiées ou non de la chaleur (non reproduites ici) et utilisent également des fonctions fk,h dites fonctions images des densités de courant en fonction de la position des points de la pièce chauffée par le dispositif de chauffage D. Les fonctions images fk,h donnent les 25 répartitions de puissance des induits par des inducteurs de courant unitaire pour chaque fréquence fondamental et éventuellement pour chaque harmonique pour obtenir les valeurs des courants en amplitude. Ces fonctions images sont définies notamment dans la publication W02011048316 pour le fondamental. 3 0 Selon un premier mode de réalisation, les courants cibles sont calculés en cherchant à éliminer les harmoniques du courant. Plus précisément, le principe de ce premier mode de réalisation consiste en la détermination des formes des signaux carrés qui permettent 3046018 13 d'éliminer les amplitudes des harmoniques du courant Hn de la série de Fourier (ou en d'autres termes, pour chaque inducteur et pour chaque harmonique Ath, ...Ap,h, avec h strictement supérieur à un). La détermination des formes des signaux carrés suppose l'insertion d'angles supplémentaires lors de la 5 correction du courant qui correspond à une modulation de largeur d'impulsion. La décomposition en série de Fourier de la tension par exemple aux bornes des inducteurs peut s'écrire : 00 V =IHnsin(nCO.t + ÇOn) n=1 Avec : 10 V : Tension aux bornes des inducteurs ; H, : amplitude des harmoniques çan : phase des harmoniques w : pulsation du fondamental Le réglage des amplitudes Hn des harmoniques implique qu'il y ait autant 15 d'angles sur le quart de période que d'harmoniques à contrôler lorsque l'on considère la source de tension ou de courant (définie ci-avant) constante. Les angles sont déterminés en résolvant l'équation suivante : E m 117, = 4 1(- 1)k+1 cos(nak ) nt- k=1 Avec : 20 E : la tension aux bornes des onduleurs ; al min Et a3 - a2 > min 5 peut être imposée avec min un angle qui par exemple peut être de 5°. A partir de la connaissance des condensateurs et autres paramètres fixes ou fixés et des équations ci-dessus, il est possible de calculer les angles de commande pour commander les onduleurs du dispositif de sorte que les 10 harmoniques de rang 3 et 5 (H3, H5) soient nulles. En particulier, la matrice d'impédances du système pour la configuration de chauffage envisagée, afin de refléter l'état magnétique du système pour une géométrie donnée doit être déterminée. En particulier cette détermination est possible, en temps réel, sans connaissance préalable sauf le nombre de lignes et de colonnes de la matrice d'impédances. Les équations 15 qui suivent sont données, à titre d'exemple, pour un dispositif de chauffage comprenant trois inducteurs. Toutefois, dans des variantes de réalisation, les dispositifs de chauffage peuvent comprendre plus de trois inducteurs. L'équation matricielle générale du système, en ne considérant que le courant fondamental, peut s'écrire : 20 V = Z./ Avec : V : Tensions sinusoïdales aux bornes des inducteurs ; I : Courants dans les enroulements des inducteurs 25 Z : matrice d'impédance du système pour le courant fondamental. Dans le cas présent, pour le courant fondamental, la matrice peut s'écrire sous la forme : 3046018 16 zii((o) Z12 \(w) Z13 ((co Z =z,io) Z22 (w) Z23 (w) z,i(co) Z32 (w) Z33 (w) Ou encore : R11 + .1410 R12 ± /L12° R13 ± /L13° Z = R21 ± /L21W R22 + j122 W R23 + IL23 W R31 + ./L31W R32 + /L32 W R33 + jL33 W_ Avec : 5 L' : inductance propre de chaque inducteur Lmn = Lnm : inductances mutuelles entre inducteurs R' : résistances propre de chaque inducteur Rmn = Rnm : résistances équivalente due aux courants induits. Avec la connaissance des relations électromagnétiques entre les 10 bobines et la pièce à chauffer, et notamment des fonction fk,1 définies ci- dessus pour le fondamental, il est possible de procéder au calcul des courants de référence du courant fondamental à injecter dans chacune des bobines afin d'obtenir le chauffage désiré. Ce premier mode de réalisation peut être en particulier appliqué pour le 15 chauffage de pièces métalliques non magnétiques (en d'autres termes des pièces ne comportant pas de matériau magnétiques). Toutefois, des pièces non métalliques mais conductrices d'électricité peuvent également être envisagées, comme mentionné précédemment. En outre, les inducteurs ne 20 présentent pas de culasses magnétiques. La mise en oeuvre du procédé selon le premier mode de réalisation, à savoir la suppression des harmoniques, peut avoir lieu selon les quatre étapes suivantes. enmes Dans une première étape, les températures effectives 8. -Unes, e2mes, mesurées de la pièce sont comparées avec les valeurs de température de référence 8. _lref, e2ref, ---enref, de sorte à déterminer la carte des densités de 3046018 17 5 puissance à injecter dans la pièce. Les températures effectives 8. _1 mes, e2mes, - - - enmes peuvent être mesurées, comme précédemment décrit par une caméra thermique ou par des thermocouples ou encore par d'autres moyens de mesure de température d'une pièce connus. Dans une deuxième étape, on utilise les fonctions de forme fk,1 donnant les répartitions de puissance induits par des inducteurs de courant unitaire pour obtenir les valeurs des courants en amplitude (11,i, ..., lp,1) et en phase (y2,1, ..., (pp,i) des inducteurs Ind1 à Indp. La connaissance de la matrice 10 d'impédance Z1 et des condensateurs permet de déterminer les courants de référence du courant fondamental au niveau des onduleurs 01 à Op. Dans une troisième étape, il est déterminé, pour chaque onduleur 01, ... Op, les courants harmoniques lond_k,h et leur phase pour chaque harmonique h, avec h strictement supérieur à 1. A titre d'exemple non limitatif, 15 h peut être égal à trois et cinq, ou plus et/ou avec des harmoniques paires. Dans une quatrième étape, la comparaison entre les courants de référence du fondamental des inducteurs et ceux réellement mesurés permet de calculer les informations de correction des ordres de pilotage des différents onduleurs pour la fréquence fondamentale (ordres Mk,1) ainsi que pour les 20 harmoniques (ordres Mk,h). On effectue ainsi autant de calculs de correction qu'il y a d'harmoniques à contrôler. Le nombre d'harmoniques à contrôler peut varier en fonction de la finesse souhaitée du profil de température. En l'espèce, le contrôle (et donc la suppression) des harmoniques d'ordre 3 et 5 permet 25 d'affiner le profil de température de façon particulièrement satisfaisante. Ainsi, on compense les harmoniques, par exemple d'ordres 3 et 5 par réglage des angles de commande des onduleurs. Comme visible sur la figure 3, le profil de référence D prefl , - - - D prefn est associé d'une part aux fonctions de forme fk,1 donnant les répartitions de 30 puissance induits par des inducteurs de courant unitaire et d'autre part à la matrice d'impédances Z1,..., Zn. Les capacités des condensateurs C1, ...Cp 3046018 18 sont calculées et on obtient ainsi les courants cibles à injecter dans la pièce. Les fondamentaux des courants cibles sont comparés aux courants mesurés puis envoyés vers des unités de traitement CORR1,1, ..., CORRp,h. Les courants harmoniques sont également envoyés vers des unités de traitement 5 CORRi,i, ..., CORRp,h correspondantes. Les unités de traitement envoient les données aux modulateurs M1,1, ...Mth qui commandent les onduleurs 01, ..Op. Dans un deuxième mode de réalisation, les harmoniques sont contrôlées et utilisées pour affiner le profil de température et améliorer la qualité du chauffage et le rendement global du système. Ce mode de 10 réalisation peut être appliqué à des onduleurs de courant ou de tension. En outre, ce mode de réalisation permet notamment le chauffage par induction de pièces métalliques en matériaux comportant des métaux ferromagnétiques, comme par exemple, à titre non limitatif, fontes et aciers. Le procédé permet également d'être utilisé avec des dispositifs comprenant des inducteurs qui 15 présentent des culasses magnétiques 20 de canalisation du flux magnétique, tel que représenté sur la figure 2. En effet, les culasses magnétiques et/ou les métaux ferromagnétiques génèrent également des courants harmoniques dans les inducteurs qui peuvent venir perturber le signal de commande et l'identification et la commande des harmoniques permettent de gérer ces 20 signaux perturbateurs. Comme précédemment mentionné, l'onduleur peut être un onduleur de courant ou un onduleur de tension. Par exemple, l'onduleur peut être un onduleur de tension ou de courant trois phase, et par suite avec trois inducteurs. Toutefois, un onduleur de tension ou de courant n phases, par suite 25 avec un dispositif comportant n inducteurs pourrait être également utilisé. Par analogie à l'écriture de la densité de courant totale induite dans la matière, il est possible d'écrire le courant total induit par les composantes de rang m (On considère notamment toutes les harmoniques contrôlées, impaires ou paires). La résultante des contributions de toutes les harmoniques est donnée par 30 l'équation suivante : 3046018 _1 19 i=1 i=1 Avec : ' : 'Harmonique de rang m du courant de l'inducteur i et fR : composantes réelles et imaginaires des fonctions de distribution des 5 courants induits dans la matière. Le vecteur des courants des inducteurs correspondant aux harmoniques de rang m est noté . = p11 La démarche sous-jacente de ce deuxième mode de réalisation consiste 10 à déterminer les courants " des inducteurs, ainsi que leurs phases respectives gin . Les amplitudes des courants et leurs phases sont calculées pour l'atteinte d'un profil de température de référence go: _lref, e2ref, - - -enref- Le calcul des courants des inducteurs passe par l'écriture de la fonction 15 à optimiser qui consiste en la superposition des termes à chaque fréquence multiple de la fréquence principale FP. La fréquence principale FP est par exemple de 1500 Hertz. De façon non limitative, il est possible par exemple de se limiter à la connaissance des fonctions de distribution des courants induits pour les harmoniques de rang 1 (correspondant au fondamental), de rang 3, de 20 rang 5 et de rang 7. Toutefois, la présente invention est envisagée avec la connaissance des fonctions de distribution des courants induits pour les harmoniques de rang 1 (fondamental), de rang 2 et jusqu'au rang n. Les courants des trois inducteurs sont décomposés en parties réelles et imaginaires pour chacune des fréquences multiples de la fréquence principale 25 FP (correspondant au fondamental) selon l'équation suivante (dans l'équation suivante, à titre d'exemple non limitatif le nombre d'harmonique m est limité à 1, 3, 5 et 7, toutefois l'homme du métier comprendra que le nombre d'harmoniques peut être supérieur à 7, notamment peut être n, avec n un entier naturel quelconque) : 3046018 20 Avec : Où : : partie réelle du courant induit par les harmoniques de rang m des courants des inducteurs )1f/ : partie imaginaire du courant induit par les harmoniques de rang m des 10 courants des inducteurs. La puissance à injecter caractérisant la rapidité et le profil de chauffage, est imposée en calculant la densité de courant de consigne correspond à l'équation suivante (dans l'équation suivante, à titre d'exemple non limitatif le 15 nombre d'harmonique m est limité à 7, toutefois l'homme du métier comprendra que le nombre d'harmoniques peut être supérieur à 7, notamment peut être n, avec n un entier naturel quelconque et n'est pas limité à la valeur 7) : 3 z=1 5 FIS 20 Avec : : densité de puissance de référence aux abscisses r, p : résistivité du matériau supposée uniforme en tout point considéré de la pièce chauffée. La fonction à optimiser est ainsi la suivante : itetfr - Fer. = 25 3046018 21 Cette fonction est à optimiser pour le fondamental et chaque harmonique h d'ordre 2 à n. La solution X = (Im1, (pm1, Imn, (pmn) est calculée pour l'obtention d'une distribution fixée par une valeur de référence, des courants induits en tout point 5 de la pièce chauffée selon des étapes d'optimisation, notamment selon les étapes d'optimisation représentées par le logigramme de la figure 6. En particulier, selon le logigramme représenté sur la figure 6, on assigne différentes valeurs à X de sorte à comparer pour ces différentes valeurs les fonctions F afin de trouver la solution optimale X = (Im1, (pm1, Imn, (pmn) pour 10 laquelle la fonction F, en valeur absolue, est la plus importante. En d'autres termes, par récurrence on calcule la valeur de X qui permet d'optimiser la fonction F. Plus précisément, le logigramme représenté sur la figure 6 permet de générer N vecteurs parmi lesquels est choisi celui qui offre la plus faible erreur du courant total induit par rapport à la valeur de référence. 15 On s'intéresse alors aux profils homogènes de densité obtenus dans différentes conditions. Bien que le principe soit générale, les figures 7 et 8 représentent un exemple limité aux harmoniques impaires de rang 3, 5 et 7 de superposition des effets harmoniques à partir des courants inducteurs. La première colonne 20 représente les fonctions de distribution des courants induits, la deuxième colonne la superposition des effets des harmoniques et la troisième colonne donne le calcul des températures. Dans un premier cas, les signaux peuvent avoir la même origine. Dans ce cas, il y a un démarrage simultané de tous les signaux à l'instant t=0. Cela 25 revient à considérer pour variables, seulement les amplitudes des harmoniques et imposer un zéro aux phases. Les contraintes imposées peuvent alors être écrites (dans l'équation suivante, à titre d'exemple non limitatif le nombre d'harmoniques m est limité à 1, 3, 5 et 7, toutefois l'homme du métier comprendra que le nombre d'harmoniques peut être supérieur à 7, notamment 30 peut être n, avec n un entier naturel quelconque) : 3046018 22 m = Ï,1, Dans un deuxième cas, les phases des signaux sont libres. Les contraintes imposées peuvent alors être écrites (dans l'équation suivante, à titre d'exemple non limitatif le nombre d'harmonique m est limité à 1, 3, 5 et 7, 5 toutefois l'homme du métier comprendra que le nombre d'harmoniques peut être supérieur à 7, notamment peut être n, avec n un entier naturel quelconque) : x = 10 Dans un troisième cas, les phases des signaux sont liées et référencées au courant de l'onduleur 1 à chaque fréquence. Ce cas consiste à imposer les phases v171= et à écrire (dans l'équation suivante, à titre d'exemple non limitatif le nombre d'harmoniques m est limité à 1, 3, 5 et 7, toutefois l'homme du métier comprendra que le nombre d'harmoniques peut être supérieur à 7, 15 notamment peut être n, avec n un entier naturel quelconque) : 5: 1000 A partir des données qui précèdent, il est possible de déterminer par 20 calcul que l'adjonction d'harmoniques améliore la précision du profil de chauffe à condition de régler leurs déphasages respectifs, ce qui complexifie l'optimisation. Le troisième cas notamment avec l'insertion de déphasages entre les fondamentaux et entre fondamentaux et harmoniques permet d'améliorer sensiblement la précision théorique du profil de chauffe. Toutefois, 25 le troisième cas induit aussi une augmentation des pertes par commutation des onduleurs. La figure 8 représente une vue générale des étapes utilisées dans le deuxième mode de réalisation pour les harmoniques 2 à n du courant et pour le 3046018 23 fondamental. La deuxième colonne regroupe les matrices d'impédance, la troisième colonne les fonctions de distribution des courants induits, la quatrième colonne la superposition des effets des harmoniques et la dernière colonne le calcul des températures. 5 En partant de la connaissance des harmoniques des courants dans chaque inducteur, il est possible de remonter aux courants des onduleurs. Les composantes harmoniques des courants des onduleurs permettent de calculer les angles de commande pour le courant fondamental, pour les harmoniques de 2 à n. On écrit : Avec : courant de l'onduleur modulé par largeur d'impulsions 15 : terme de couplage de l'inducteur i et k à la fréquence m * f, f : la fréquence du fondamental, 271 : impédance propre de l'inducteur i en présence du matériau à chauffer Dans une première étape et pour un onduleur de courant, les 20 condensateurs placés en parallèle (comme représenté sur la figure 5) avec les inducteurs sont calculés. Bien entendu, dans le cas d'onduleurs de tension, les condensateurs sont en série et non en parallèle. Ces condensateurs sont calculés pour une compensation globale de l'énergie réactive portée par les fondamentaux des courants des inducteurs. Il faut à ce niveau, prendre en 25 compte l'ensemble des échanges existants entre les inducteurs. L'écriture matricielle où l'on retrouve les termes, z f','J et 41, caractérise ces échanges à la fréquence principale (en l'occurrence par exemple 1500 Hz) et permet de calculer les capacités des condensateurs. Les tensions 17,hisont calculées par l'équation rond j 10 hl _ 30 3046018 24 Les condensateurs sont calculés par l'équation : Les courants des onduleurs sont calculés à partir des courants inducteurs par la relation qui fait intervenir la tension de l'onduleur suivante : 5 = I * 1- Avec : m= 1 à n 10 Le calcul des angles de commande suppose la résolution sur chaque phase du système d'équations (dans l'équation suivante, à titre d'exemple non limitatif le nombre d'harmoniques est limité à 1, 3, 5 et 7, toutefois l'homme du métier comprendra que le nombre d'harmoniques peut être supérieur à 7, notamment peut être n, avec n un entier naturel quelconque) : coc.i 1 4 = = P3 = P5 15 e4) = p7 Les valeurs pi sont caractéristiques des amplitudes des harmoniques et liées au courant de la source qui constitue également une variable. Le paramètre p est calculé par l'équation qui met en relation les harmoniques des courants des onduleurs à obtenir, en fonction de la source de 20 courant continue qui les alimente : Pm = Avec : i : 'indice de l'onduleur, Is : valeur de la source de courant, 25 m : rang de l'harmonique. L'équation ci-dessus s'applique à la variante de réalisation avec des onduleurs de courant. Le même principe s'applique pour des onduleurs de tension. 3046018 25 Les équations qui précèdent peuvent être implémentées et mises en oeuvre dans le procédé d'optimisation de chauffage selon le deuxième mode de réalisation qui comporte les étapes suivantes. Dans une première étape, on compare les températures mesurées avec 5 les valeurs de référence et on détermine la carte des densités de puissance à injecter dans la pièce. Dans une deuxième étape, on utilise les fonctions de forme fk,h donnant les répartitions de puissance induites par des inducteurs de courant unitaire et pour chaque fréquence fondamentale ou harmonique pour obtenir les valeurs 10 des courants en amplitude (11,1, ..., lp,1), (11,h, ..., lp,h) et en phase (cp2,1, ..., (pp,1), (cp1,h, ..., (pp,h) des inducteurs Indl à Indp. Dans une troisième étape, la connaissance de la matrice d'impédances Z1 et des condensateurs permet de déterminer les courants de référence du courant fondamental au niveau des onduleurs 01 à Op. La connaissance des 15 matrices d'impédances Zh, avec h égale à 2, 3, 4...n (3, 5, 7 par exemple), et des condensateurs permettent de déterminer les courants de référence du courant harmonique de rang h au niveau des onduleurs 01 à Op. Dans une quatrième étape, la comparaison entre les courants de référence des inducteurs et ceux réellement mesurés (obtenue par 20 décomposition en série de Fourier) permet de calculer les informations de correction des ordres de pilotage des différents onduleurs pour la fréquence fondamentale (ordres Mk,1 ) ainsi que pour les harmoniques (ordres Mk,h). La boucle de rétroactions du deuxième mode de réalisation, illustrée sur 25 la figure 5, est en partie semblable à celle du premier mode de réalisation. Le profil de référence Dpref1, ...Dprefn est associé d'une part aux fonctions de forme fk,1 donnant les répartitions de puissance induites par des inducteurs de courant unitaire et d'autre part à la matrice d'impédances Z1,..., Zh. Les capacités des condensateurs C1, ...Cp sont calculées et on obtient ainsi les 30 courants cibles à injecter dans la pièce. Les fondamentaux des courants cibles sont comparés aux courants mesurés puis envoyés vers des unités de 3046018 26 traitement CORR1,1, ..., CORRo. Les courants harmoniques sont également envoyés vers des unités de traitement CORR1,h, ..., CORRp,r, correspondantes. Les unités de traitement envoient les données aux modulateurs M1,1, ... Mp,h 5 qui commandent les onduleurs 01, ..Op. En outre, pour les courants harmoniques, on associe aux profils de référence D pre, -.Dprefil les fonctions de forme fk,r, donnant les répartitions de puissance induits par des inducteurs de courant unitaire pour les harmoniques. A partir de la matrice d'impédances Zh et du calcul du vecteur des capacités C, on obtient pour les harmoniques les 10 courants cibles lp,h. Ces valeurs sont envoyées vers les modulateurs Mp,h. Ainsi, la stratégie de pilotage des onduleurs par modulation de largeur d'impulsions permet de définir complètement les ordres de pilotage par la fourniture de deux paramètres seulement pour chaque harmonique qu'on souhaite contrôler. L'ensemble de chauffage par induction comprenant un dispositif de 15 chauffage par induction tel que précédemment mentionné et une pièce forme un système permettant de s'adapter aux contraintes industrielles de façon rapide et en continu. REVENDICATIONS1. Procédé d'optimisation de chauffage par induction mis en oeuvre dans un dispositif de chauffage (D) d'une pièce (10), le dispositif comprenant des inducteurs (Ind1, Ind2, ..., Indp) couplés magnétiquement, chaque inducteur étant alimenté par un onduleur de tension ou de courant (01, 02, ..., Op) qui lui est propre et associé à un condensateur (Ci, C2, ...Cp) pour former un circuit oscillant, lesdits circuits oscillants possédant sensiblement la même fréquence de résonance, chaque onduleur (01, 02, ..., Op) étant commandé par une unité de commande (1V11,1,---,Mp,h) de façon à faire varier l'amplitude (A1,1, ...Ap,n) et la phase 41,1 , ...(p p,h) du courant (I1,...1p) traversant l'inducteur correspondant, le dispositif (D) comprenant en outre des moyens de détermination dudit courant (11, ...1p) ainsi que des moyens de détermination d'un profil de température effectif (A ,-imes, e2mes, - - -enmes) de ladite pièce chauffée (10), ledit procédé comprenant les étapes suivantes : a. on compare ledit profil de température effectif ue ,,-imes, e2mes, - - - enmes) à un profil de température de référence (A ,-iref, e2ref, - - - enref); et on calcule un profil de densité de puissance de référence (Dpref1; Dpref2, - - - Dprefn) que le dispositif de chauffage doit injecter dans ladite pièce pour atteindre ledit profil de température de référence (8. , -lref; e2ref, - - - enref) b. on calcule des courants cibles (Iref1; Iref2; - - .1refp) appropriées pour injecter dans ladite pièce ledit profil de densité de puissance de référence (Dpref 1; Dpref 2, - - - Dprefn) par : - la connaissance de fonctions images vectorielles (fk,l; - - -f k,h) représentatives des relations liant les densités de courant créées par les inducteurs aux courants traversant les inducteurs de sorte à déterminer les valeurs de l'amplitude des courants (11, 1, -,1p, h),et de la phase des courants 41,1, ...(p p, h) des inducteurs (Ind1, Ind2, ..., Indp), - la connaissance d'une matrice d'impédances (Z1, Zn) déterminée par la connaissance des relations électromagnétiques liant lesdits inducteurs entre eux et à ladite pièce, - la valeur des condensateurs (Ci, ...Cp) des circuits oscillants, 3046018 28 - l'identification et la détermination des amplitudes et de la phase des harmoniques du courant (Iond_k, h) délivré par chaque onduleur (01, 02, ..., Op) de sorte à contrôler le fondamental et les harmoniques du courant dans l'alimentation d'un inducteur (Ind1, Ind2, .. Ando) par son 5 onduleur (01, 02, ..., Op) associé, 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 10 onduleurs de façon à corriger les courants traversant les inducteurs pour la fréquence fondamentale ainsi que pour les harmoniques. 2. Procédé selon la revendication 1, dans lequel les harmoniques du courant délivré par chaque onduleur (01, 02, ..., Op) sont optimisés de sorte à améliorer la qualité du chauffage par induction, et dans lequel, on 15 calcule la capacité électrique (C1, ...Cp) des condensateurs pour compenser l'énergie réactive portée par les fondamentaux des courants des inducteurs (11,i, ..10) en fonction de la connaissance des relations électromagnétiques liant lesdits inducteurs entre eux. 3. Procédé selon la revendication 1 ou 2, dans lequel lors de 20 l'étape b), on utilise en outre des fonctions images vectorielles (fk, h) représentatives des relations liant les densités de courant créées par les inducteurs aux courants traversant les inducteurs pour chaque harmonique considérée du courant de sorte à obtenir les valeurs des amplitudes du courant (li, h, -,Ip, h) et de la phase des courants (y) 1, n, ...(p p, h) des inducteurs (lind1, 25 lind2, ..., lindp). 4. Procédé selon l'une des revendications 1 à 3, dans lequel lors de l'étape b) les valeurs cibles (Ip,href) sont calculées par la connaissance de la matrice d'impédances (Zi) déterminée par la connaissance des relations électromagnétiques liant lesdits inducteurs (Ind1, Ind2, ..., Indp) entre eux et à 30 ladite pièce (10) pour le fondamental du courant et la connaissance des matrices d'impédance (Zn) pour chaque harmonique que l'on souhaite contrôler. 5. Procédé selon l'une quelconque des revendications 3046018 29 précédentes, dans lequel on détermine les capacités des condensateurs (C1, C2, ...Cp) et on associe la matrice d'impédances (Z1, Zh), à un vecteur des capacités. 6. Procédé selon la revendication 1, dans lequel une modulation de 5 largeur d'impulsions est réalisée par insertion d'angles supplémentaires afin de supprimer des harmoniques de rang supérieur ou égale à deux des courants traversant les inducteurs (li, ...Io, dans lequel des harmoniques de rang supérieur ou égal à 2 sont supprimées par un signal carré, et dans lequel le fondamental du courant délivré par chaque onduleur (01, ...Op) est contrôlé 10 par le même signal carré. 7. Dispositif de chauffage par induction comprenant: - des inducteurs (Ind1, Ind2, ..., Indp) couplés magnétiquement, chaque inducteur (Ind1, Ind2, ..., Indp) étant associé à un condensateur (Ci, C2, ...Cp) pour former un circuit oscillant, lesdits circuits oscillants possédant au moins 15 approximativement la même fréquence de résonance; - des onduleurs (01, 02, ..., Op) alimentant chacun un inducteur (Ind1, Ind2, ..., Indp) qui lui est propre, chaque onduleur (01, 02, ..., Op) étant commandé par une unité de commande (Mp,h) de façon à faire varier l'amplitude et la phase du fondamental et des harmoniques du courant traversant l'inducteur 20 correspondant ; comprenant en outre : - des moyens de détermination du courant fondamental et des harmoniques des courants traversant les inducteurs Ind1, Ind2, ..., Indp) ainsi que des moyens de détermination d'un profil de température effectif d'une pièce 25 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 30 de température de référence; - des moyens de calcul, basés sur la connaissance d'une matrice des impédances (Z1, 4), de courants cibles (Irefi, Iref2, ...lrefp ) que doivent délivrer les onduleurs (01, ...Op) afin que les courants des inducteurs et leurs 3046018 30 harmoniques 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 à 5 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. 8. Dispositif de chauffage par induction selon la revendication 7, 10 dans lequel les inducteurs (Ind1, Ind2, ..., Indp) présentent des culasses magnétiques de canalisation du flux magnétique. 9. Ensemble de chauffage par induction comprenant un dispositif de chauffage (D) par induction selon la revendication 7 ou 8 et une pièce (10) chauffée par le dispositif de chauffage (D) par le procédé selon l'une 15 quelconque des revendications 1 à 7, ladite pièce chauffée étant conductrice d'électricité. 10. Ensemble de chauffage selon la revendication 9, dans lequel la pièce (10) est métallique. 11. Ensemble selon la revendication 10, dans lequel la pièce 20 chauffée comprend un métal ferromagnétique.