Analyse SFRA pour la Détection de Déformation de Bobinage

𝗗𝗘𝗙𝗢𝗥𝗠𝗔𝗖𝗜𝗢𝗡 𝗜𝗡𝗩𝗜𝗦𝗜𝗕𝗟𝗘 𝗗𝗘 𝗗𝗘𝗩𝗔𝗡𝗔𝗗𝗢𝗦: 𝗘𝗟 𝗥𝗜𝗘𝗦𝗚𝗢 𝗦𝗜𝗟𝗘𝗡𝗖𝗜𝗢𝗦𝗢 𝗧𝗥𝗔𝗦 𝗨𝗡 𝗖𝗢𝗥𝗧𝗢𝗖𝗜𝗥𝗖𝗨𝗜𝗧𝗢 Un transform

Ing. Francisco Ramírez

Fondements Physico-Mathématiques de l'Analyse de Réponse en Fréquence par Balayage (SFRA)

L'Analyse de Réponse en Fréquence par Balayage (SFRA, pour son sigle en anglais Sweep Frequency Response Analysis) constitue la méthodologie non destructive de plus haute résolution pour évaluer l'intégrité mécanique et hydraulique-structurelle de l'ensemble noyau-enroulements dans les transformateurs de puissance. Le principe fondamental de cette méthode réside dans le fait qu'un transformateur de puissance constitue un réseau passif complexe à paramètres distribués RLC (résistance, inductance, capacité et conductance), dont la fonction de transfert spectrale est bijective par rapport à sa géométrie physique et à son hyperstructure diélectrique.

Aux moyennes et hautes fréquences, la géométrie de l'enroulement, les distances diélectriques entre galettes (ou spires), l'intervalle d'isolement par rapport à la cuve et la position relative entre les enroulements haute tension, basse tension et de réglage déterminent une empreinte digitale spectrale unique. Toute altération dimensionnelle imperceptible à l'échelle macroscopique — induite par des forces électrodynamiques de court-circuit, des séismes ou des contraintes mécaniques lors du transport — modifie la distribution locale des inductances et des capacités, provoquant un déplacement univoque des pôles et des zéros de la fonction de transfert du système.

Modèle à Paramètres Distribués des Enroulements

Afin de modéliser analytiquement le comportement spectral d'un enroulement, celui-ci est subdivisé en éléments infinitésimaux ou discrets au sein d'un réseau en échelle de type Ladder Network. Chaque section élémentaire dxdx le long de l'axe axial de l'enroulement possède ses propres inductances, des inductances mutuelles par rapport aux sections adjacentes et non adjacentes, des capacités série entre spires/galettes, des capacités par rapport à la terre (noyau et cuve), ainsi que des pertes diélectriques et ohmiques associées.

Les équations aux dérivées partielles qui régissent la propagation de l'onde de tension v(x,t)v(x,t) et de courant i(x,t)i(x,t) le long de l'enroulement correspondent à la généralisation des Équations des Télégraphistes dans le domaine fréquentiel spectral :

V(x,s)x=(R(x,s)+sL(x,s))I(x,s)0LwsM(x,ξ)I(ξ,s)dξ-\frac{\partial V(x, s)}{\partial x} = \left( R(x, s) + s L(x, s) \right) I(x, s) - \int0^{Lw} s M(x, \xi) I(\xi, s) \, d\xi
I(x,s)x=(G(x,s)+sCg(x))V(x,s)x(sCs(x)V(x,s)x)-\frac{\partial I(x, s)}{\partial x} = \left( G(x, s) + s Cg(x) \right) V(x, s) - \frac{\partial}{\partial x} \left( s Cs(x) \frac{\partial V(x, s)}{\partial x} \right)

Où :

  • s=σ+jω=j2πfs = \sigma + j\omega = j2\pi f représente la fréquence complexe de Laplace.
  • R(x,s)R(x, s) est la résistance distribuée dépendante de la fréquence en raison de l'effet de peau (skin effect) et de l'effet de proximité dans les conducteurs (RfR \propto \sqrt{f}).
  • L(x,s)L(x, s) est l'inductance propre distribuée par unité de longueur.
  • M(x,ξ)M(x, \xi) est le noyau de l'opérateur intégral d'inductance mutuelle entre la position axiale xx et ξ\xi.
  • Cg(x)Cg(x) est la capacité shunt ou distribuée à la terre (entre l'enroulement et le noyau magnétique mis à la terre ou la paroi de la cuve).
  • Cs(x)Cs(x) est la capacité série distribuée (capacité entre spires consécutives et capacité entre galettes adjacentes).
  • G(x,s)G(x, s) est la conductance transversale de l'isolement diélectrique cellulose/huile (G=ωCgtanδG = \omega Cg \tan\delta).
  • LwLw représente la longueur totale équivalente du conducteur de l'enroulement.

Formulation de la Fonction de Transfert

L'essai SFRA injecte un signal sinusoïdal d'amplitude constante Vin(f)Vin(f) (typiquement 10 V crête à crête) à une borne de l'enroulement avec un balayage logarithmique de fréquence s'étendant usuellement de 10Hz10 Hz à 2MHz2 MHz, tout en mesurant la réponse en tension Vout(f)Vout(f) à une autre borne, les deux canaux étant fermés sur une impédance de référence standard de 50 Ω50\ \Omega.

La réponse est généralement exprimée sur une échelle logarithmique en décibels représentant l'atténuation ou le gain de tension :

H(f)dB=20log10Vout(f)Vin(f)H(f)_{ dB } = 20 \log10 \left| \frac{Vout(f)}{Vin(f)} \right|

Alternativement, en termes de réponse en phase spectrale ϕ(f)\phi(f) :

ϕ(f)=arg(Vout(f)Vin(f))=arctan(Im(H(f))Re(H(f)))\phi(f) = \arg\left( \frac{Vout(f)}{Vin(f)} \right) = \arctan \left( \frac{ Im (H(f))}{ Re (H(f))} \right)

Dans le domaine opérationnel de Laplace, l'impédance d'entrée Zin(s)Zin(s) ou la fonction de transfert de tension H(s)H(s) peut être exprimée de manière rigoureuse par un développement en pôles (pjp_j) et zéros (ziz_i) :

H(s)=Ki=1m(szi)j=1n(spj)=Ksm+am1sm1++a1s+a0sn+bn1sn1++b1s+b0(nm)H(s) = K \cdot \frac{\prod_{i=1}^{m} (s - z_i)}{\prod_{j=1}^{n} (s - p_j)} = K \cdot \frac{s^m + a_{m-1}s^{m-1} + \dots + a_1 s + a_0}{s^n + b_{n-1}s^{n-1} + \dots + b_1 s + b_0} \quad (n \ge m)

Les fréquences de résonance (maxima locaux, pôles) et d'antirésonance (minima locaux, zéros) sont régies par les racines du polynôme caractéristique. Pour une section élémentaire d'enroulement discret d'ordre NN, les fréquences résonantes primaires ωk\omega_k sont paramétrées analytiquement par :

ωkkπLeq(Cg+k2π2Cs)\omega_k \approx \frac{k \pi}{\sqrt{Leq \left( Cg + k^2 \pi^2 Cs \right)}}

Cette relation mathématique démontre de manière incontestable l'extrême sensibilité de la méthode : une altération géométrique millimétrique modifiant l'écartement entre galettes alterne la capacité série CsC_s, ce qui déplace directement la fréquence propre ωk\omega_k dans les zones de moyenne et haute fréquence du spectre.

Segmentation par Bandes de Fréquence et Diagnostic Sous-Structurel

Le spectre complet de la fonction de transfert SFRA (10Hz10 Hz à 2MHz2 MHz) ne répond pas de manière homogène à l'ensemble des composants physiques du transformateur. En raison des comportements différenciés de la réactance inductive (XL=2πfLX_L = 2\pi f L) et de la réactance capacitive (XC=12πfCX_C = \frac{1}{2\pi f C}), la réponse fréquentielle est analysée selon quatre bandes fonctionnelles distinctes, chacune étant dominée par un sous-ensemble structurel spécifique de l'appareil.

Bande de Très Basse et Basse Fréquence (< 2 kHz) : Domaine Magnétisant du Noyau

À des fréquences inférieures à 2kHz2 kHz, la réactance capacitive des isolements est extrêmement élevée (XCX_C \to \infty), se comportant pratiquement comme un circuit ouvert. En conséquence, le courant injecté circule principalement à travers l'inductance magnétisante non linéaire du noyau magnétique (LmL_m) et la résistance en courant continu des enroulements (RdcRdc).

La fonction de transfert dans cette zone est dominée par une première antirésonance profonde suivie d'une résonance correspondant au circuit équivalent constitué par l'inductance magnétisante du noyau en parallèle avec la capacité équivalente par rapport à la terre :

fcore=12πLm(Cg+Cinter)f_{ core } = \frac{1}{2\pi \sqrt{L_m \left( C_g + Cinter \right)}}

Tout phénomène modifiant la perméabilité magnétique du noyau μr\mu_r, la continuité des tôles d'acier au silicium, l'état de la mise à la terre du noyau ou la présence d'un magnétisme rémanent après un défaut entraînera une variation massive de l'amplitude et de la position spectrale dans cette bande basse.

Bande de Moyenne-Basse Fréquence (2 kHz à 100 kHz) : Interaction Inter-Enroulements et Structures Globales

Dans cette plage spectrale, la réactance de l'inductance magnétisante augmente considérablement et la capacité inter-enroulements (CinterCinter) commence à offrir un chemin de faible impédance. La signature dans cette bande reflète l'interaction mutuelle entre les enroulements concentriques de Haute Tension (HT), de Basse Tension (BT) et les enroulements de Réglage (Taps).

Le comportement est régi par l'inductance de fuite en court-circuit (LkL_k) et la capacité inter-enroulements CinterCinter. Les variations dans cette zone révèlent une inclinaison globale des enroulements, des déplacements axiaux globaux de l'ensemble du bloc d'enroulements ou des désordres au niveau des connexions internes et des changeurs de prises en charge (OLTC).

Bande de Moyenne-Haute Fréquence (100 kHz à 500 kHz) : Structure Propre de l'Enroulement et Géométrie Axiale/Radiale

À ces fréquences, l'inductance de fuite LkL_k et les capacités série entre galettes ou entre spires (CsC_s) forment un réseau multiphasique complexe de pôles et de zéros. L'acier du noyau magnétique agit pratiquement comme un écran au flux magnétique en raison des courants de Foucault induits à la surface des tôles, ce qui réduit la perméabilité effective du noyau à μr1\mu_r \approx 1.

Par conséquent, la réponse dépend de manière quasi exclusive de la géométrie interne propre à l'enroulement. Les déplacements locaux de galettes, la déformation par flambage radial (radial buckling), l'inclinaison des spires ou galettes (tilting) ou le tassement des cales de pressage axial provoquent des variations directes de la capacité locale inter-galettes Cs(x)C_s(x) et de l'inductance de fuite incrémentale ΔLk\Delta L_k, se traduisant par un glissement horizontal (en fréquence) et vertical (en atténuation) de multiples pics résonants.

Bande de Haute Fréquence (500 kHz à > 2 MHz) : Câblage de Sortie, Connexions et Traversées Isolantes

À l'extrémité supérieure du spectre (> 500 kHz), l'atténuation est dominée par les capacités parasites de plus faible valeur, les conducteurs de liaison interne (sorties/leads), les traversées isolantes (bushings) et l'inductance propre des barres de raccordement.

Toute modification observable dans cette zone qui ne s'accompagne pas de variations aux fréquences inférieures est généralement liée à des différences de mise à la terre de l'équipement de mesure, à une dégradation des connexions des traversées, à une variation de la position physique des liaisons de descente du transformateur, ou à des défauts sur les écrans de blindage électrostatique.

Analyse Forensique des Modes de Défaillance Mécanique et Transformations du Spectre

Les forces électrodynamiques ou forces de Lorentz s'exerçant sur les conducteurs d'un transformateur lors d'un événement de court-circuit externe sont proportionnelles au carré du courant instantané (FIsc2F \propto Isc^2). Ces forces possèdent deux composantes vectorielles fondamentales : la force radiale (FrF_r) et la force axiale (FzF_z). L'action de ces forces produit des déformations plastiques irréversibles dans la structure mécanique de l'enroulement, qui se traduisent par des empreintes spectrales SFRA extrêmement spécifiques.

Flambage Radial (Radial Buckling)

Les forces radiales dirigées vers l'intérieur agissent en comprimant mécaniquement l'enroulement interne (souvent l'enroulement de basse tension entourant directement le noyau). Si la force radiale dépasse la limite critique de déformation plastique ou de flambage du conducteur en cuivre, le phénomène connu sous le nom de Buckling se produit.

On distingue deux typologies mécaniques principales :

  • Flambage Libre (Free Buckling) : L'enroulement se déforme en adoptant un profil lobulaire sans point fixe, formant des crêtes et des creux sur toute sa circonférence.
  • Flambage Forcé (Forced Buckling) : L'enroulement s'effondre vers l'intérieur entre les barrettes de calage ou cales axiales de la colonne du noyau, générant des déformations périodiques écrasées.

Effet sur le Spectre SFRA : La déformation par flambage réduit la distance physique radiale entre l'enroulement interne et le noyau, augmentant considérablement la capacité par rapport à la terre (CgC_g \uparrow). De plus, la symétrie de l'espace inter-enroulements est altérée, modifiant la capacité mutuelle (CinterCinter). Sur le tracé SFRA, ce défaut se diagnostique par un décalage systématique vers les basses fréquences des pics de résonance dans la bande de 100 kHz à 500 kHz, accompagné d'une modification d'amplitude des creux d'antirésonance.

Déplacement Axial et Désalignement (Axial Displacement & Tilting)

Les forces axiales proviennent de l'asymétrie du flux magnétique de fuite aux extrémités supérieure et inférieure de l'enroulement. Si les enroulements HT et BT ne sont pas parfaitement alignés magnétiquement, une composante de force axiale destructrice apparaît, tendant à pousser l'enroulement intérieur vers le bas et l'enroulement extérieur vers le haut.

Si la pression axiale exercée excède la précharge du système de serrage mécanique (anneaux de maintien et cales axiales), deux modes de défaillance surviennent :

  1. Déplacement Axial Global : L'enroulement complet glisse axialement le long du noyau.
  2. Inclinaison des Galettes (Disk Tilting) : Les galettes d'enroulement adjacentes basculent tangentiellement ou se déforment en "zigzag" sous l'effet de l'effondrement des blocs d'isolement inter-galettes.

Effet sur le Spectre SFRA : Le déplacement axial modifie drastiquement le couplage magnétique mutuel et induit une variation de l'inductance de fuite en court-circuit (LkL_k \uparrow). Sur la signature SFRA, le déplacement axial génère une altération marquée dans la bande de moyenne-basse fréquence (2 kHz - 100 kHz), où l'on observe un glissement notable du premier creux de résonance inter-enroulements vers la gauche (fréquences plus basses) ainsi que des variations d'amplitude supérieures à 3-6 dB sur les pics intermédiaires.

Court-Circuit Entre Spires ou Claquage Diélectrique Inter-Galettes

Une défaillance diélectrique découlant de la dégradation du papier isolant Kraft ou d'une surtension transitoire de manœuvre ou de foudre peut court-circuiter deux ou plusieurs spires consécutives. Un court-circuit entre spires annule le flux magnétique au sein de la boucle en court-circuit en vertu de la loi de Lenz, induisant un courant circulant interne massif qui masque et annule l'inductance de ce segment d'enroulement.

Effet sur le Spectre SFRA : L'annulation de l'inductance équivalente provoque une perte dramatique de la capacité de stockage d'énergie magnétique. Dans la réponse SFRA, la signature aux basses et moyennes fréquences perd complètement les pics résonants caractéristiques de l'enroulement affecté, transformant la courbe de fonction de transfert en un profil plat à forte atténuation dans les bandes basse et moyenne, similaire à la réponse d'un circuit fortement amorti ou en court-circuit permanent.

Algorithmes de Comparaison Statistique et Critères Normatifs

L'interprétation de l'essai SFRA repose principalement sur la comparaison d'empreintes spectrales au travers de trois méthodologies standards :

  • Comparaison par Rapport à l'État Initial (Baseline / Fingerprint vs. Actual) : Comparaison directe avec l'empreinte de référence mesurée en usine ou avant la mise en service. C'est la méthode garantissant la plus haute sensibilité et certitude diagnostique.
  • Comparaison Phase à Phase (Phase-to-Phase) : Comparaison entre les phases symétriques d'une même machine (ex. Phase A vs. Phase C) mesurées dans la même configuration.
  • Comparaison avec un Transformateur Jumeau (Sister Transformer) : Comparaison avec une unité de conception mécaniquement identique du même constructeur.

Pour éliminer la subjectivité de l'inspection visuelle des courbes de Bode, les normes internationales telles que l'IEEE C57.149, l'IEC 60076-18 et les recommandations du CIGRE WG A2.26 préconisent l'utilisation de métriques analytiques avancées, au premier rang desquelles le Coefficient de Corrélation Croisée (RxyRxy ou CCCC) et la Distance Absolue du Spectre.

Coefficient de Corrélation Croisée (RxyRxy)

Étant donné deux vecteurs de réponse en fréquence exprimés en décibels, Xk=Hbase(fk)X_k = Hbase(f_k) et Yk=Hactual(fk)Y_k = Hactual(f_k), évalués de manière discrète sur un nombre NN de points au sein de sous-bandes de fréquence prédéfinies :

Rxy=k=1N(XkXˉ)(YkYˉ)k=1N(XkXˉ)2k=1N(YkYˉ)2Rxy = \frac{\sum_{k=1}^{N} \left( X_k - \bar{X} \right) \left( Y_k - \bar{Y} \right)}{\sqrt{\sum_{k=1}^{N} \left( X_k - \bar{X} \right)^2 \cdot \sum_{k=1}^{N} \left( Y_k - \bar{Y} \right)^2}}

Xˉ\bar{X} et Yˉ\bar{Y} correspondent aux moyennes arithmétiques des amplitudes de réponse dans la sous-bande évaluée :

Xˉ=1Nk=1NXk,Yˉ=1Nk=1NYk\bar{X} = \frac{1}{N} \sum_{k=1}^{N} X_k, \quad \bar{Y} = \frac{1}{N} \sum_{k=1}^{N} Y_k

Métrique de Distance Absolue du Spectre (DabsDabs)

Pour quantifier les écarts d'amplitude sans qu'ils ne soient excessivement déformés par de purs déphasages ou des décalages spectraux de pics étroits, on utilise la Distance Absolue Spectrale Moyenne :

Dabs(X,Y)=1Nk=1NXkYkDabs(X, Y) = \frac{1}{N} \sum_{k=1}^{N} \left| X_k - Y_k \right|

Critères de Décision Diagnostique

Conformément au guide DL/T 911-2016 et aux normes IEEE/IEC, les seuils de sévérité structurés pour le coefficient de corrélation croisée dans la bande critique de l'enroulement (100kHz500kHz100 kHz - 500 kHz) sont classés selon l'échelle rigoureuse suivante :

EˊtatDiagnostique={Normal/Intact,siRxy0.998DeˊformationLeˊgeˋre/Marginale,si0.990Rxy<0.998DeˊformationModeˊreˊe,si0.950Rxy<0.990DeˊformationSeˊveˋre/DeˊfaillanceMeˊcanique,siRxy<0.950État Diagnostique = \begin{cases} Normal / Intact , & si Rxy \ge 0.998 \\ Déformation Légère / Marginale , & si 0.990 \le Rxy < 0.998 \\ Déformation Modérée , & si 0.950 \le Rxy < 0.990 \\ Déformation Sévère / Défaillance Mécanique , & si Rxy < 0.950 \end{cases}

Matrice Comparative des Phénoménologies, Critères Normatifs et Impact Diagnostique

Le tableau suivant regroupe de manière exhaustive les signatures spectrales, les variations des paramètres électriques équivalents, les normes applicables, les seuils métriques et les conséquences opérationnelles pour les modes de défaillance détectables par SFRA :

Mode de Défaillance Mécanique / Structurelle Bande de Fréquence Affectée Variation Paramétrique RLC Équivalente Comportement de la Signature Spectrale SFRA Critère Métrique Normatif (IEC 60076-18 / IEEE C57.149) Impact Diélectrique ou Opérationnel Critique
Flambage Radial (Radial Buckling) Moyenne-Haute Fréquence
(100 kHz – 500 kHz)
CgC_g \uparrow (10% - 30%)
CinterCinter \downarrow
LkvariableL_k variable
Déplacement horizontal des pics résonants vers les basses fréquences. Variation d'amplitude des pics. Rxy<0.95Rxy < 0.95 en bande moyenne-haute.
Dabs>3.5dBDabs > 3.5 dB.
Réduction des distances d'isolement dans l'huile. Claquage diélectrique inter-enroulements par décharges partielles.
Inclinaison / Déplacement Axial des Galettes Moyenne-Basse à Moyenne-Haute
(20 kHz – 300 kHz)
CsC_s \uparrow ou \downarrow
LkL_k \uparrow (1% - 5%)
MalteˊreˊeM altérée
Modification sévère de l'amplitude des pics/creux. Apparition de nouvelles paires de pôles et zéros parasites. Rxy<0.98Rxy < 0.98 en bande moyenne.
ΔLk>2%\Delta L_k > 2\% lors de l'essai de réactance de fuite.
Cisaillement des cales d'isolation clés. Perte totale de la précharge de serrage mécanique axial.
Perte de Mise à la Terre du Noyau Basse Fréquence
(10 Hz – 1 kHz)
CgflottanteC_g flottante
G(parasite)G \to \infty (parasite)
Disparition du creux caractéristique à basse fréquence. Décalage vertical du tracé vers le haut. Déviation extrême avec Rxy<0.90Rxy < 0.90 dans la sous-bande basse (< 1 kHz). Génération de décharges destructives par potentiel flottant. Gazeification massive de l'huile diélectrique (DGA).
Magnétisme Rémanent dans le Noyau Très Basse Fréquence
(10 Hz – 500 Hz)
μr(saturationlocale)\mu_r \downarrow (saturation locale)
LmL_m \downarrow
Déplacement du premier pic de résonance vers des fréquences plus hautes. La signature redevient normale après démagnétisation. Rxy<0.95Rxy < 0.95 en bande < 1 kHz, mais normalisation absolue après cycle de démagnétisation DC. Courants d'appel (inrush currents) sévères lors de la mise sous tension. Enclenchement intempestif des protections différentielles (87T).
Court-Circuit Entre Spires Basse, Moyenne et Haute
(10 Hz – 500 kHz)
Lm0L_m \to 0
LkdrastiqueL_k \downarrow drastique
R(pertes)R \uparrow (pertes)
Effondrement total de la signature spectrale. La courbe perd ses pics et devient une réponse platement atténuée. Rxy<0.80Rxy < 0.80 simultanément dans plusieurs bandes. Inductance magnétisante détruite. Défaillance catastrophique immédiate. Déclenchement de la protection Buchholz (63) et surintensité instantanée (50).
Déplacement des Connexions et Inclinaison des Sorties (Leads) Haute Fréquence
(500 kHz – 2 MHz)
LleadalteˊreˊeLlead altérée
CstrayalteˊreˊeCstray altérée
Glissement des crêtes uniquement dans la zone supérieure à 500 kHz. Stabilité absolue à basse et moyenne fréquence. Rxy<0.95Rxy < 0.95 exclusivement pour f>500kHzf > 500 kHz. Sous-bandes inférieures maintenues à Rxy>0.998Rxy > 0.998. Risque faible de défaillance immédiate, sauf réduction possible des distances diélectriques par rapport à la cuve.

Stratégies d'Atténuation Mécanique et Modélisation Quantitative des Forces Électrodynamiques

La prévention des dommages structurels détectés par SFRA nécessite une modélisation rigoureuse des forces électrodynamiques s'exerçant sur la géométrie des enroulements lors de défauts de court-circuit triphasés ou biphasés à la terre. La conception doit garantir que les contraintes dynamiques ne dépassent pas les limites élastiques du cuivre ou des alliages de cuivre dopé à l'argent (Cu-Ag).

Modélisation des Forces Radiales de Court-Circuit

La densité de force radiale par unité de longueur FrF_r appliquée à un enroulement circulaire soumis à un courant crête maximal instantané de court-circuit dissymétrique ipeakipeak est calculée par :

Fr=μ0(Nipeak)2πDm2heqF_r = \frac{\mu_0 \cdot (N \cdot ipeak)^2 \cdot \pi \cdot Dm}{2 \cdot heq}

Où :

  • μ0=4π×107H/m\mu_0 = 4\pi \times 10^{-7} H/m est la perméabilité magnétique du vide.
  • NN est le nombre total de spires de l'enroulement.
  • ipeakipeak est le courant crête de court-circuit (incluant la composante continue d'asymétrie Kpeak2IscKpeak \cdot \sqrt{2} Isc).
  • DmD_m est le diamètre moyen de l'enroulement.
  • heqheq est la hauteur magnétique équivalente de l'enroulement en tenant compte des champs de fuite aux extrémités.

La contrainte mécanique de traction/compression moyenne (contrainte tangentielle, Hoop Stress, σθ\sigma_{\theta}) induite sur la section transversale du conducteur AcA_c s'exprime analytiquement par :

σθ=Fr2πNAc=μ0Nipeak2Dm4heqAc\sigma_{\theta} = \frac{F_r}{2\pi \cdot N \cdot A_c} = \frac{\mu_0 \cdot N \cdot ipeak^2 \cdot D_m}{4 \cdot heq \cdot A_c}

Pour éviter tout flambage permanent, la condition de conception stricte impose que :

σθσadmissible=σ0.2FS\sigma_{\theta} \le \sigma_{ admissible } = \frac{\sigma_{0.2}}{ FS }

σ0.2\sigma_{0.2} est la limite d'élasticité conventionnelle à 0,2 % de déformation du matériau conducteur (typiquement 120220MPa120 - 220 MPa pour le cuivre écroui) et FSFS est le facteur de sécurité statique/dynamique (FS1.5FS \ge 1.5).

Calcul des Forces Axiales de Compression et de Serrage

La force axiale totale de compression FzF_z agissant pour tasser l'enroulement vers son centre magnétique moyen est décrite par l'intégration du produit du flux radial de fuite Br(z)B_r(z) et du courant total :

Fz(z)=2π0zDmNdiripeakBr(ξ)dξF_z(z) = 2\pi \int0^{z} D_m \cdot Ndir \cdot ipeak \cdot B_r(\xi) \, d\xi

La contrainte d'écrasement axial sur les cales d'isolement en cellulose ne doit pas dépasser la résistance à la compression dynamique à chaud du carton comprimé haute densité (Pressboard, typiquement σcomp4070MPa\sigma_{\text{comp}} \le 40 - 70 MPa).

Techniques Constructives Modernes d'Atténuation

  • Conducteurs Transposés Continus Collés à la Résine Époxy (Bonded CTC) : Utilisation de câbles CTC cuits avec un vernis thermo-adhérent au stade B (B-stage). Cela lie mécaniquement tous me brins individuel en cuivre en une poutre rigide monolithique, augmentant le module de section et multipliant la résistance au flambage radial par un facteur de 3 à 5 par rapport aux conducteurs individuels.
  • Serrage Axial Dynamique sous Pression par Ressorts ou Vérins Hydrauliques : Mise en place de dispositifs d'application de précharge constante pour compenser automatiquement le fluage (creep) du papier et du bois stratifié d'isolation tout au long de la durée de vie de l'appareil, garantissant que l'enroulement ne perde jamais la force de maintien axiale minimale requise (Fprecarga>Fz, dinaˊmico Fprecarga > F_{\text{z, dinámico }}).

Intégration Avancée avec Vexten Suite pour l'Analyse SFRA et l'Intégrité Structurelle

Au sein de l'écosystème analytique de Vexten Suite, le module spécialisé de diagnostic SFRA agit comme le moteur central pour l'expertise forensique et la prédiction de la santé structurelle des transformateurs de puissance. Le logiciel intègre le traitement analytique des données spectrales avec des simulations déterministes du réseau électrique.

Correlación de Défauts de Court-Circuit selon l'IEC 60909 et l'IEEE 141

Le moteur de calcul de court-circuit de Vexten Suite (module Vexten Fault Engine) simule les transitoires électromagnétiques de défaut dans les réseaux de transport et de distribution selon les méthodologies IEC 60909 et IEEE 141. Lors de la survenance d'un défaut à proximité de la sous-station, Vexten Suite calcule avec précision le courant de crête asymétrique instantané (ipi_p) qui a traversé le transformateur :

ip=κ2Iki_p = \kappa \cdot \sqrt{2} \cdot I_k''

κ\kappa est le facteur de crête fonction du rapport R/XR/X du réseau au point de défaut et IkI_k'' est le courant de court-circuit initial symétrique.

À partir de ce courant calculé, Vexten Suite détermine la cartographie tridimensionnelle des contraintes mécaniques dynamiques (Fr,FzF_r, F_z) subies par l'enroulement. Si les limites de fatigue accumulée du conducteur sont dépassées, le système alerte automatiquement l'ingénieur de maintenance et émet une ordre d'intervention prioritaire pour la réalisation d'un essai SFRA in situ.

Algorithme Vexten d'Empreinte Spectrale et Décomposition Automatique en Sous-Bandes

Le module logiciel Vexten SFRA Analyzer automatise l'importation des fichiers de tracés normalisés XML/ASCII (formats CIGRE / IEEE C57.149). Il exécute de façon native les processus mathématiques avancés suivants :

  1. Alignement Spectral par Spline Cubique : Interpolation logarithmique de la réponse pour harmoniser des jeux de données fréquentiels discontinus issus d'équipements de test différents (ex. Omicron, Doble, Megger).
  2. Décomposition Adaptative en Quatre Sous-Bandes Standards : Segmentation automatique du vecteur spectral dans les régions physiques du noyau, de l'interaction inter-enroulements, de la structure propre des enroulements et des connexions.
  3. Matrice de Coefficient de Correlación Croisée Vectorielle et Métrique SDM (Spectrum Distance Metrics) : Exécution parallèle d'analyses statistiques multivariables entre l'empreinte originale de référence et les mesures récentes de terrain.
CCVexten=[RxyBajaRxyMedBajaRxyMedAltaRxyAlta]    MoteurdInfeˊrenceetDeˊtectionduBuckling/Tilting\mathbf{CC}_{ Vexten } = \begin{bmatrix} Rxy^{ Baja } \\ Rxy^{ Med-Baja } \\ Rxy^{ Med-Alta } \\ Rxy^{ Alta } \end{bmatrix} \implies Moteur d'Inférence et Détection du Buckling / Tilting

Si la métrique calculée par Vexten Suite dans la bande moyenne-haute passe en dessous de 0.980.98, le module déclenche des alertes précoces et croise ces résultats avec l'essai de capacité et facteur de dissipation de l'isolement (C1/C2) ainsi que la mesure de l'inductance de fuite en court-circuit (Short-Circuit Impedance test selon l'IEC 60076-5), isolant ainsi le mode exact de défaillance et estimant le pourcentage de vie utile mécanique restante du transformateur.