Pertes par Hystérésis et Foucault dans l'Acier CRGO

Thermodynamique Électromagnétique et Microstructure de l'Acier au Silicium CRGO

Le circuit magnétique (noyau) d'un transformateur de puissance agit comme le circuit de couplage de flux fondamental entre les enroulements, opérant sous des régimes cycliques non linéaires d'excitation électromagnétique. Le rendement global et la densité de puissance du transformateur sont conditionnés de manière déterminante par le comportement microstructural et cristallographique du matériau ferromagnétique utilisé. Le standard industriel prédominant pour les transformateurs de moyenne, haute et très/ultra haute tension (HTB/THT/EHV/UHV) est l'acier électrique à grains orientés laminé à froid (Cold-Rolled Grain-Oriented Silicon Steel, CRGO), optimisé par des traitements thermomécaniques afin d'aligner ses axes cristallographiques de facile aimantation selon la direction de laminage.

L'acier CRGO présente une structure cristalline cubique centrée (CC / BCC, réseau $\alpha$-Fe), où l'axe cristallographique $\langle 100 \rangle$ constitue la direction d'énergie minimale d'anisotropie magnétocristalline, tandis que les directions $\langle 110 \rangle$ et $\langle 111 \rangle$ représentent respectivement les axes d'anisotropie intermédiaire et dure. Par une séquence rigoureuse de laminage à froid et de recuits de recristallisation secondaire (inhibée par des précipités de sulfure de manganèse, $MnS$, ou de nitrure d'aluminium, $AlN$), se développe la texture dite de Goss, caractérisée par l'orientation cristallographique $\{110\}\langle 001 \rangle$. Dans cette configuration, les plans $\{110\}$ sont situés parallèlement à la surface de la tôle et les directions $\langle 001 \rangle$ s'orientent parallèlement à la direction longitudinale de laminage (DL / RD, Rolling Direction).

Où $E_a$ est la densité d'énergie d'anisotropie magnétocristalline ($J/m^3$), $K_1 \approx 4{,}8 \times 10^4 \ J/m^3$ et $K_2 \approx 1{,}5 \times 10^4 \ J/m^3$ sont les constantes d'anisotropie du fer à température ambiante ($298 \ K$), et $\alpha_1, \alpha_2, \alpha_3$ représentent les cosinus directeurs du vecteur d'aimantation $\mathbf{M}$ par rapport aux axes de la maille unitaire. La désorientation angulaire moyenne de la texture de Goss par rapport à la direction de laminage définit la qualité de l'acier : dans les aciers CRGO conventionnels, la dispersion angulaire oscille entre $5^\circ$ et $7^\circ$, tandis que dans les aciers à haute perméabilité (Hi-B), elle est réduite à des plages de $2^\circ$ à $3^\circ$.

L'addition de silicium ($Si$) à des concentrations massiques nominales de $3{,}0\%$ à $3{,}4\%$ augmente la résistivité électrique volumique ($\rho$) d'environ $10 \times 10^{-8} \ \Omega\cdot m$ (pour le fer pur) jusqu'à $(45 - 50) \times 10^{-8} \ \Omega\cdot m$. Cet accroissement de résistivité atténue drastiquement le transport libre des porteurs de charge responsables des courants de Foucault. Cependant, des concentrations supérieures à $3{,}5\%$ en masse induisent une fragilité mécanique sévère due à la formation de phases ordonnées intermétalliques $Fe_3Si$ et $FeSi$ ($B2$ et $D0_3$), empêchant les procédés de mise en forme à froid à l'échelle industrielle. De surcroît, le silicium réduit la constante de magnétostriction de saturation longitudinale ($\lambda_{100}$), minimisant ainsi les déformations élastiques induites par le champ magnétique et l'émission acoustique qui en résulte.

La structure en domaines magnétiques (domaines de Weiss) dans le CRGO est composée de domaines principaux orientés à $180^\circ$ séparés par des parois de Bloch, et de domaines de fermeture à $90^\circ$ au voisinage des discontinuités de surface et des inclusions. La largeur caractéristique du domaine ($2d_0$) dépend de l'équilibre entre l'énergie magnétostatique de surface et l'énergie élastique de la paroi de domaine. L'application de revêtements diélectriques sous tension superficielle (comme le revêtement de phosphates et silicates, désigné commercialement sous le nom de Carlite) introduit une traction biaxiale permanente dans la tôle ($\sigma_t \approx 4 - 10 \ MPa$), réduisant la largeur moyenne des domaines à $180^\circ$ et éliminant les domaines secondaires transversaux de fermeture, ce qui diminue radicalement les pertes par courants de Foucault anomaux.

Tripartition Classique et Spectre des Pertes : Le Modèle de Bertotti

Sous excitation magnétodynamique cyclique, la dissipation volumique d'énergie par unité de temps dans le noyau est quantifiée par l'intégrale de contour fermée sur le cycle du cycle d'hystérésis dynamique. La formulation générale des pertes spécifiques de puissance ($P_t$, en $W/kg$) est régie par la théorie statistique des pertes électromagnétiques formulée par Giorgio Bertotti, laquelle dépasse les limitations des formulations empiriques classiques de Steinmetz en décomposant le phénomène en trois mécanismes dissipatifs physiquement découplés :

Où $f$ représente la fréquence fondamentale d'excitation ($Hz$), $B_p$ est la valeur crête de l'induction magnétique ($T$), $\sigma$ est la conductivité électrique du matériau ($S/m$), $d$ est l'épaisseur de la tôle laminée ($m$), $\rho_m$ est la masse volumique de l'acier ($kg/m^3$, typiquement $7650 \ kg/m^3$), $k_h$ est le coefficient d'hystérésis quasi-statique, $\alpha$ est l'exposant de Steinmetz (fréquemment $1{,}6 \le \alpha \le 2{,}0$), et $kexc$ est le paramètre microscopique des pertes en excès.

Pertes par Hystérésis Quasi-Statique

Les pertes par hystérésis statique ($P_h$) trouvent leur origine dans les processus microscopiques de dissipation thermodynamique irréversible lorsque les parois des domaines magnétiques se déplacent à travers un milieu ferromagnétique présentant des défauts de réseau cristallin, tels que des dislocations, des lacunes, des joints de grains, des contraintes résiduelles et des inclusions non métalliques (e.g., précipités de $SiO_2$ ou de $Fe_3C$). Durant le mouvement d'une paroi de Bloch à $180^\circ$, celle-ci se retrouve ancrée mécaniquement sur ces sites de piégeage (pinning sites).

Pour libérer la paroi de domaine du puits de potentiel de piégeage, une augmentation du champ magnétique externe jusqu'à atteindre le champ critique local est nécessaire. Une fois cette barrière énergétique franchie, la paroi avance de manière discontinue à grande vitesse vers le puits de potentiel suivant par un saut microscopique non linéaire, phénomène connu sous le nom d'effet Barkhausen. L'énergie dissipée lors de ce saut est convertie de manière irréversible en phonons (chaleur) au sein du réseau cristallin :

L'aire du cycle d'hystérésis quasi-statique ($w_h$) dépend exclusivement de la structure métallurgique du matériau et demeure indépendante de la fréquence de variation temporelle du flux pour les régimes où le temps de relaxation du domaine est infiniment plus court que la période de l'excitation. La modélisation mathématique de $P_h$ en présence de formes d'onde complexes ou de saturation non linéaire requiert la formulation de l'opérateur continu de Preisach :

Où $\hat{\gamma}_{\alpha\beta}$ est un opérateur hystérétique élémentaire bistable présentant des seuils de basculement $\alpha$ (d'activation) et $\beta$ (de désactivation), et $\mu(\alpha, \beta)$ est la fonction de distribution de densité de Preisach propre au matériau CRGO caractérisé, identifiable expérimentalement par des familles de courbes de renversement de premier ordre (FORC, First-Order Reversal Curves).

Pertes par Courants de Foucault Classiques

Les pertes classiques par courants de Foucault ($Pcl$) dérivent rigoureusement de l'application macroscopique des équations de Maxwell à un conducteur homogène et isotrope continu de géométrie plane finie. Considérons une tôle d'acier d'épaisseur $d$ selon la direction de l'axe $z$ (de $z = -d/2$ à $z = +d/2$), de largeur infinie sur l'axe $y$, et de longueur finie sur l'axe $x$. En admettant une perméabilité magnétique effective $\mu$, une conductivité volumique $\sigma = 1/\rho$, et un vecteur densité de flux magnétique parallèle à la surface : $\mathbf{B}(t) = B_p \sin(\omega t) \hat{\mathbf{x}}$.

Sous l'hypothèse où l'épaisseur $d$ est substantiellement inférieure à l'épaisseur de peau électromagnétique (skin depth) $\delta = \sqrt{2/(\omega \sigma \mu)}$, l'induction magnétique est spatialement uniforme à travers la section transversale de la tôle : $B_x(z,t) \approx B_p \sin(\omega t)$. En intégrant par rapport à $z$ à partir du plan neutre central ($z = 0$, où $E_y = 0$ par symétrie) :

La densité de courant induite $\mathbf{J}(z,t) = \sigma \mathbf{E}(z,t)$ dissipe une puissance instantanée par effet Joule dont la densité volumique locale est $p(z,t) = J_y^2(z,t) / \sigma = \sigma E_y^2(z,t)$. En intégrant sur l'épaisseur totale de la tôle et en effectuant la moyenne temporelle sur une période $T = 2\pi/\omega$ :

En divisant par la masse volumique $\rho_m$, on obtient la formulation spécifique normalisée :

Cette dérivation démontre la dépendance quadratique directe vis-à-vis de l'épaisseur $d$ de la tôle et de la fréquence fondamentale $f$, justifiant la nécessité de réduire les calibres commerciaux de $0{,}35 \ mm$ (M4) jusqu'à $0{,}23 \ mm$ ou $0{,}18 \ mm$ (grades traités par laser à ultra-faibles pertes).

Pertes Anomales ou Pertes en Excès (Excess Losses)

Historiquement, la somme directe des pertes par hystérésis quasi-statique et des pertes classiques de Foucault s'avérait significativement inférieure aux pertes totales mesurées expérimentalement par la méthode du wattmètre sur cadre d'Epstein ($P_{total, mes} > P_h + Pcl$). Cet écart énergétique a été formellement qualifié d'anomalie de pertes.

Le modèle microscopique de Pry et Bean (1958) a démontré que le flux magnétique ne varie pas uniformément dans tout le continuum de la tôle, mais se concentre exclusivement au sein des domaines magnétiques et se transfère par le déplacement localisé des parois de Bloch à $180^\circ$. La vitesse de déplacement d'une paroi de domaine $v_w(t)$ sous excitation sinusoïdale induit un gradient de champ électrique local infiniment plus prononcé au voisinage immédiat de la paroi que celui prédit par la théorie classique du continuum :

Où $L$ représente l'espacement entre parois de domaines et $B_s$ est l'induction de saturation. Le micro-courant de Foucault confiné dans le micro-volume adjacent à la paroi génère un champ de freinage magnétique local. Bertotti a généralisé cette physique en formulant la théorie statistique des Objets Magnétiques Actifs (Magnetic Objects, MO), compris comme des ensembles cohérents de parois de domaines interagissant avec les hétérogénéités locales du milieu :

Où $G = 0{,}1356$ est une constante adimensionnelle associée à la dissipation du champ de Foucault en milieu infini, $S$ est l'aire de la section transversale de la tôle perpendiculaire à la direction du flux magnétique, et $V_0$ est un paramètre intrinsèque du matériau homogène à un champ magnétique ($A/m$), quantifiant l'interaction stochastique entre les objets magnétiques et le champ microscopique de piégeage (internal pinning field distribution).

Grade CRGO (Classification AISI / EN)	Épaisseur Nominale $d$ (mm)	Induction Maximale $B$ à $800 \ A/m$ (T)	Pertes Spécifiques $P_{1{,}7/50}$ (W/kg)	Pertes Spécifiques $P_{1{,}7/60}$ (W/kg)	Facteur de Foisonnement / Empilage (\%)
M4 (AISI 35G155 / EN 10107)	0,35	1,82 - 1,84	1,25 - 1,35	1,64 - 1,77	$\ge 96{,}0$
M3 (AISI 30G130 / EN 10107)	0,30	1,84 - 1,86	1,05 - 1,15	1,38 - 1,51	$\ge 95{,}5$
Hi-B Conventionnel (27Q110)	0,27	1,90 - 1,93	0,95 - 1,05	1,24 - 1,38	$\ge 95{,}0$
Hi-B Haute Perméabilité (23Q090)	0,23	1,92 - 1,95	0,80 - 0,90	1,05 - 1,18	$\ge 94{,}5$
Hi-B Traitement Laser (20Q080-L)	0,20	1,93 - 1,96	0,70 - 0,78	0,92 - 1,02	$\ge 94{,}0$
Métal Amorphe (Fe-Si-B Metglas 2605SA1)	0,025	1,56 - 1,58	0,18 - 0,25	0,24 - 0,33	$\ge 86{,}0$

Effets Non-Linéaires, Distorsion Harmonique et Composante Continue (DC Bias)

Dans les réseaux électriques contemporains, la pénétration massive de l'électronique de puissance (onduleurs centralisés pour parcs photovoltaïques, liaisons HVDC, dispositifs FACTS et variateurs de vitesse) engendre des régimes de magnétisation non sinusoïdaux et des tensions à fort taux de distorsion harmonique ($THD_v$). La réponse non linéaire de la courbe d'aimantation dynamique $B(H)$ transforme ces contraintes en pertes énergétiques supplémentaires sévères.

Excitation Harmonique et Cycles d'Hystérésis Mineurs

Lorsque la forme d'onde de tension appliquée contient des harmoniques de rang supérieur ($h = 3, 5, 7, \dots$), la dérivée temporelle de l'induction magnétique $\partial \mathbf{B}/\partial t$ présente de multiples passages par zéro au cours d'un cycle fondamental. Physiquement, cela force l'inversion transitoire du mouvement des parois de domaines magnétiques avant que la polarisation macroscopique totale ne soit complétée, générant ainsi des cycles d'hystérésis mineurs (minor hysteresis loops) imbriqués au sein de la boucle principale de saturation.

Sous tension purement sinusoïdale, $B_h \propto V_h / (h f_1)$. Si la tension présente une distorsion harmonique caractérisée par des composantes $V_h$, le champ harmonique induit est approximé par $B_{p,h} \approx \frac{V_h}{h \, V_1} B_{p,1}$. En substituant dans l'expression des courants de Foucault classiques :

Cependant, les pertes en excès ne suivent pas ce découplage linéaire. L'interaction multi-fréquentielle accélère la vitesse locale des parois de Bloch selon une intégrale temporelle de la dérivée de flux :

Magnétisation avec Polarisation Continue (DC Bias)

L'injection de courant continu dans les enroulements du transformateur — causée par les courants induits géomagnétiquement (GIC), les topologies d'onduleurs à couplage direct sans transformateur d'isolement ou les défauts asymétriques monopolaires sur les lignes HVDC à retour par la terre — produit une translation du cycle dynamique selon l'axe du champ magnétisant ($Hdc = N Idc / lfe$).

Ce décalage asymétrique repousse le matériau vers le coude de saturation sur une demi-alternance complète (saturation semi-cyclique). Dans la région de saturation ($B > 1{,}9 \ T$ pour le CRGO), la perméabilité différentielle $\mu_d = dB/dH$ s'effondre drastiquement pour atteindre des valeurs proches de la perméabilité du vide ($\mu_0 = 4\pi \times 10^{-7} \ H/m$). Il en résulte des pointes massives de courant magnétisant asymétrique ($imag(t)$ avec des facteurs de crête supérieurs à 10), transformant le champ magnétique confiné au circuit laminé en un flux de fuite tridimensionnel de forte intensité qui traverse les pièces structurelles adjacentes (poutres de serrage, cuve, tirants et écrans magnétiques).

La puissance réactive consommée $Q$ augmente d'un ordre de grandeur, provoquant des chutes de tension sévères aux bornes, l'apparition d'harmoniques pairs et impairs de bas rang (particulièrement $h = 2$ et $h = 3$), et des échauffements thermiques destructifs par courants induits dans les éléments structuraux métalliques.

Analyse Forensique des Défaillances Thermiques et Dégradation Diélectrique Induite par le Circuit Magnétique

Les anomalies électromagnétiques et les dégradations thermiques prenant naissance dans le circuit magnétique constituent un mode de défaillance majeur des transformateurs de grande puissance. Contrairement aux avaries mécaniques des enroulements provoquées par les efforts électrodynamiques radiaux et axiaux de court-circuit, les défaillances du noyau se manifestent par des processus de dégradation électrochimique et thermique progressive sur de longues périodes d'exploitation.

Dégradation de l'Isolation Interlaminaire (Revêtement de Type Carlite)

L'isolation interlaminaire de surface des tôles magnétiques (épaisseur $\sim 1{,}5 - 3{,}0 \ \mu m$, composée d'une couche d'accrochage en forstérite $Mg_2SiO_4$ recouverte d'une couche de phosphate d'aluminium $AlPO_4$ et d'acide phosphorique) est dimensionnée pour supporter des tensions interfeuilles de quelques volts ($1{,}5 \ V$ à $5{,}0 \ V$ efficaces selon la norme CEI 60404-6). Toutefois, sous l'action de contraintes thermiques locales continues ($T > 150^\circ C$), de la fatigue vibratoire d'origine magnétostrictive ($100/120 \ Hz$ et harmoniques acoustiques) ou de contraintes mécaniques excessives de serrage des culasses, ce film diélectrique se fissure et se pulvérise.

La rupture de l'isolation galvanique entre tôles adjacentes crée des micro-soudures par arc ou un contact électrique franc. Dès lors que deux points de contact ou plus s'établissent au sein d'un même paquet de tôles, une boucle fermée en court-circuit de section transversale non négligeable embrasse le flux magnétique principal. La force électromotrice induite dans cette spire de défaut s'exprime par :

La résistance de boucle ohmique $Rloop$ étant de l'ordre du milliohm ($m\Omega$), le courant de circulation $Ifault$ atteint des amplitudes permanentes de plusieurs centaines à plusieurs milliers d'ampères. La puissance dissipée par effet Joule génère des points chauds localisés (hot-spots) dont les températures dépassent rapidement $500^\circ C - 800^\circ C$.

Dégradation de l'Huile Diélectrique et Analyse Forensique par DGA

Le transfert thermique intense entre le point chaud du noyau et le fluide diélectrique environnant (huile minérale naphténique, paraffinique ou esters synthétiques) entraîne une pyrolyse et un craquage thermocatalytique des chaînes hydrocarbonées ($C_n H_{2n+2}$). Les liaisons covalentes carbone-hydrogène ($C-H$, énergie de liaison $\sim 413 \ kJ/mol$) et carbone-carbone ($C-C$, $\sim 348 \ kJ/mol$ ; $C=C$, $\sim 614 \ kJ/mol$ ; $C \equiv C$, $\sim 839 \ kJ/mol$) se scindent en fonction du gradient thermique local.

Conformément aux normes internationales CEI 60599 et IEEE C57.104, le profil des gaz dissous dans l'huile (Dissolved Gas Analysis, DGA) permet de déterminer la sévérité thermique et la nature physique du défaut :

Plage Thermique de Défaut	Gaz Combustible Dominant	Rapports Diagnostiques (CEI 60599 / Duval)	Mécanisme Physique au Cœur du Noyau	Action Corrective Requise
Défaut Thermique $T < 300^\circ C$ (T1)	Méthane ($CH₄$) / Éthane ($C₂H₆$)	$\frac{C₂H₂}{C₂H₄} < 0{,}1$, $0{,}1 < \frac{CH₄}{H₂} < 1$, $\frac{C₂H₄}{C₂H₆} < 1$	Échauffement diffus des culasses, dégradation superficielle des vernis interfoliaires ou des peintures de blocage.	Suivi chromatographique mensuel ; contrôle du spectre harmonique réseau.
Défaut Thermique $300^\circ C < T < 700^\circ C$ (T2)	Éthylène ($C₂H₄$)	$\frac{C₂H₂}{C₂H₄} < 0{,}1$, $\frac{CH₄}{H₂} > 1$, $1 < \frac{C₂H₄}{C₂H₆} < 4$	Spire en court-circuit entre tôles CRGO multiples. Décomposition avancée du revêtement minéral.	Réduction de charge à $70\%$ ; mesure des pertes à vide et test d'isolement du circuit magnétique.
Défaut Thermique $T > 700^\circ C$ (T3)	Éthylène ($C₂H₄$) avec traces d'Acétylène ($C₂H₂$)	$\frac{C₂H₂}{C₂H₄} < 0{,}15$, $\frac{CH₄}{H₂} > 1$, $\frac{C₂H₄}{C₂H₆} > 4$	Fusion localisée du paquet de tôles (core burning) ; court-circuit franc tôle-châssis de serrage.	Déclenchement forcé immédiat. Inspection endoscopique interne et réfection intégrale du noyau.
Décharges Partielles / Étincelage dans le Noyau	Hydrogène ($H₂$) dominant avec $CH₄$	$\frac{CH₄}{H₂} < 0{,}1$, $H₂$ au-delà des seuils normaux (> 100 ppm)	Potentiel flottant du noyau (rupture de la liaison équipotentielle unique), créant des micro-arcs capacitifs vers la cuve.	Vérification de la continuité du circuit de mise à la terre externe du noyau au mégohmmètre sous $2{,}5 \ kV$.

Défaillance par Mises à la Terre Multiples (Boucle de Masse et Courants de Circulation)

Selon les règles de conception normatives, le circuit magnétique d'un transformateur doit être relié à la terre en un point unique au moyen d'une tresse ou lame de cuivre raccordée à une traversée basse tension dédiée sur le couvercle de la cuve. Cette mise à la terre unique évite la montée en potentiel électrostatique flottant induit par couplage capacitif avec les enroulements haute tension ($Vind = C_{w-c}/(C_{w-c} + C_{c-g}) \cdot Vwinding$).

En cas de défaillance de l'isolation des tirants de serrage (tie-rods), d'accumulation de ponts de particules métalliques au fond de la cuve ou d'erreur de montage établissant une seconde liaison accidentelle à la masse à l'autre extrémité du noyau, une boucle de surface importante est formée, interceptant le flux principal :

Cette tension induit un courant de circulation permanent à travers les blindages, la cuve et les tôles de mise à la masse, atteignant couramment $10 \ A$ à plus de $150 \ A$. Ce courant détruit l'isolation des goujons d'assemblage, calcine l'huile adjacente et libère des concentrations massives d'éthylène ($C₂H₄$) et de méthane ($CH₄$), sans qu'aucune protection à maximum de courant ou différentielle ($87T$) ne détecte d'anomalie sur les courants de phase.

Techniques Avancées d'Atténuation dans la Conception et la Fabrication

L'optimisation des performances magnétiques requiert une approche multidisciplinaire combinant l'ingénierie cristallographique, le traitement des microdomaines par rayonnement cohérent et l'optimisation géométrique tridimensionnelle des zones de jonction des culasses.

Assemblages de Culasses de Type Step-Lap à $45^\circ$

L'assemblage géométrique entre les colonnes et les culasses dans les circuits empilés engendre des discontinuités de réluctance. Dans les joints traditionnels droits à recouvrement simple à $90^\circ$ (Butt-Lap joints), le flux magnétique est contraint de passer d'une tôle à l'autre perpendiculairement à la direction de laminage (selon l'axe cristallographique dur $\langle 110 \rangle$), ce qui provoque une réfraction sévère des lignes de champ, une augmentation de la réluctance interfaciale et des saturations locales intenses aux angles de coupe.

La technologie moderne d'empilage fait appel à des coupes d'onglet à $45^\circ$ avec décalage progressif des tôles, dite Step-Lap (configurée couramment sur 5 à 7 échelons par cycle). Cet échelonnement répartit l'entrefer le long d'un gradient axial étendu :

Grâce à un décalage spatial $\Delta x \approx 3{,}5 - 7{,}0 \ mm$ entre couches successives, le flux magnétique ne franchit plus un entrefer unique à forte réluctance, mais se redistribue progressivement sur une section effective équivalente beaucoup plus large :

La géométrie Step-Lap permet de réduire les pertes localisées aux nœuds d'assemblage de $25\%$ à $30\%$, d'abaisser la puissance réactive d'excitation d'environ $40\%$ et de diminuer le bruit acoustique global de $3$ à $6 \ dB(A)$ par rapport aux architectures classiques.

Raffinement des Domaines Magnétiques par Rayonnement Laser (Laser Scribing)

Dans les aciers CRGO à haute perméabilité (Hi-B), la taille des grains cristallographiques est importante ($\sim 10 - 30 \ mm$). Bien que cette structure favorise une orientation de Goss quasi parfaite ($\le 3^\circ$), elle présente l'inconvénient d'accroître l'écartement entre parois de Bloch à $180^\circ$ ($2d_0$), qui évolue proportionnellement à la racine carrée de la dimension du grain :

Un espacement excessif entre parois accélère la vitesse instantanée de déplacement de ces dernières $v_w(t)$, ce qui accroît fortement les pertes anomales par courants de Foucault ($Pexc$).

Le procédé de raffinement des domaines par laser (Domain Refinement ou Laser Scribing) projette un faisceau continu focalisé (laser Nd:YAG ou fibre de forte puissance) transversalement à la direction de laminage selon un pas régulier ($p \approx 4 - 8 \ mm$). L'apport thermique localisé crée, sans ablation superficielle, un profil de contraintes élastiques résiduelles de compression sous la surface en raison du fort gradient thermique de refroidissement :

Ces zones périodiques sous contrainte agissent comme des barrières d'énergie magnéto-élastique qui subdivisent les larges domaines à $180^\circ$, réduisant leur largeur moyenne $2d_0$ de moitié ou des deux tiers. Par conséquent :

Il convient de distinguer les traitements laser sensibles à la température (dont l'effet s'annule lors des recuits de détensionnement à $T > 500^\circ C$ requis sur les circuits enroulés) des méthodes d'indentation mécanique ou d'attaque ionique/plasma dites "recuisables", qui conservent le fractionnement des domaines après traitement thermique à $800^\circ C$.

Implémentation et Modélisation Numérique avec la Suite Vexten

La quantification des pertes magnétiques et leur couplage multiphysique avec le réseau s'effectuent au sein de la Suite Vexten en associant les modules d'écoulement de charge harmonique non linéaire (selon IEEE 519 / IEEE C57.110), de calcul thermique dynamique (CEI 60076-7) et de simulation des transitoires magnétiques et de court-circuit (CEI 60909 / IEEE 141).

Algorithme de Correction par Distorsion Harmonique et Facteur K

En présence de spectres harmoniques élevés générés par des charges non linéaires ou des sources d'énergie renouvelable, le module Vexten PowerFlow & Harmonics réalise une partition dynamique des pertes en recalculant le coefficient de pertes fer d'après la distribution spectrale normalisée :

Le logiciel relie ce coefficient au déclassement thermique (derating) du transformateur via le calcul du Facteur $K$ et du Facteur Harmonique des Pertes de Foucault ($FHL$), régulant le courant admissible maximal afin de garantir que la température du point le plus chaud des enroulements et du circuit magnétique ne dépasse pas la classe thermique de l'isolation :

Procédure de Modélisation Multiphysique dans la Suite Vexten

1. Définition du Spectre d'Excitation Magnétique : Importer dans le module de qualité de l'onde les relevés de mesure sur site ou les spectres harmoniques normalisés de tension ($V_h$) ainsi que la composante continue résiduelle ($Idc$) issue des flux de puissance des onduleurs.
2. Sélection du Matériau Magnétique dans la Base de Données : Définir la nuance spécifique de tôle CRGO (e.g., M3, 23Q090 ou 20Q080-L traitée laser). La Suite Vexten charge automatiquement les paramètres microscopiques : conductivité $\sigma$, constante d'hystérésis $k_h$, épaisseur $d$, masse volumique $\rho_m$ et paramètre de Bertotti $V_0$.
3. Paramétrage de la Géométrie des Jonctions : Renseigner l'architecture d'assemblage des culasses : Butt-Lap $90^\circ$, Step-Lap à 5 échelons ou Step-Lap à 7 échelons. L'algorithme ajuste automatiquement le coefficient d'accroissement de réluctance interfaciale ($\mathcal{K}_{joint}$).
4. Exécution du Calcul Thermique Transitoire Couplé : Le solveur intègre les pertes volumiques totales $Ptotal(x,y,z)$ comme terme source dans le modèle thermo-fluidique conforme à la norme CEI 60076-7, déterminant la distribution tridimensionnelle des températures dans les culasses supérieures, inférieures et les noyaux.
5. Contrôle des Seuils et Diagnostic Forensique Préventif : Les gradients thermiques locaux et les densités de pertes sont confrontés aux seuils critiques de la norme IEEE C57.104. Si la température du noyau dépasse $130^\circ C$ en régime établi, des alertes de déclassement sont générées et le taux prévisionnel d'émanation des gaz de décomposition ($C₂H₄$, $CH₄$, $H₂$) est calculé, permettant l'élaboration de stratégies de maintenance prévisionnelle (CBM) fondées sur la physique de l'appareil.

Grâce à cette chaîne de calcul intégrée, les ingénieurs d'études, d'essais et de maintenance peuvent anticiper avec rigueur le vieillissement thermique du circuit magnétique, écarter les risques de court-circuit destructif entre tôles et sélectionner les matériaux ferromagnétiques les plus performants pour garantir une efficacité énergétique optimale des transformateurs de puissance.