Modélisation des courants de Foucault : Au-delà du circuit RL

La modélisation de l’impédance d’un haut-parleur au-delà de quelques centaines de Hz est un défi classique en ingénierie acoustique. Si la résistance ($R_e$) et l’inductance ($L_e$) de base suffisent pour le bas du spectre, les pertes par courants de Foucault (eddy currents) dans les pièces polaires transforment la réponse en un système complexe à mémoire, rendant les modèles simples obsolètes.

1. L’héritage : Réseaux de Cauer et limites empiriques

Traditionnellement, l’industrie s’appuie sur le modèle de Cauer ($R_e + L_e + R_2/L_2 + \dots$). En empilant des mailles $RL$ en échelle, on tente d’approximer la nature fréquentielle de la diffusion magnétique par une analogie électrique.

Bien que standard, cette approche se heurte à deux obstacles majeurs dès que l’on vise la haute précision :

• La Rigidité Numérique : Plus on cherche à couvrir de décades de fréquences, plus le système devient mathématiquement “raide” (stiff), rendant les simulations temporelles instables ou extrêmement gourmandes en ressources.
• L’Absence de Localisation : Il n’existe aucune correspondance spatiale entre les composants du réseau et le fer. La physique est “étalée”, rendant la modélisation de la saturation magnétique purement empirique et déconnectée de la géométrie réelle.

2. L’élégance Mathématique : Dérivées Fractionnaires

Pour pallier le manque de précision de Cauer, l’utilisation d’une impédance fractionnaire $Z(s) = Ls^\alpha$ (où $0 < \alpha < 1$) offre une alternative séduisante. Elle capture la pente semi-inductive avec une économie de paramètres remarquable. La formulation de Grünwald-Letnikov (GL) permet de traduire ce concept en une somme pondérée, agissant comme un filtre FIR à mémoire longue.

L’approche Podlubny

Pour extraire les poids $w_i^\alpha$ de ce filtre, la méthode de Podlubny (1999) via une IFFT est l’outil de référence. En manipulant la fonction génératrice $(1-z)^\alpha$, on obtient une précision numérique inégalée pour des approximations d’ordres supérieurs.

% FIR weight calculation via IFFT (Munroe Method)
% Extract from PhD Thesis Appendix
nPoints = 1024;
theta = linspace(0, 2*pi, nPoints);
z = exp(-1i * theta);
alpha = 0.5; % Fractional derivative order

% First order approximation (Euler)
% f2 and f3 for higher orders can also be substituted here
f = (1 - z).^alpha; 

% Transform to time domain to extract weights
F = real(ifft(f));
nMemory = 512; % FIR filter length
w = F(1:nMemory);

Le Mur de la Complexité

Cependant, l’élégance se heurte rapidement au réel. Pour maintenir une mémoire de seulement 1 seconde à un taux d’échantillonnage de 96 kHz, le processeur doit effectuer 96 000 multiplications par échantillon.

Le problème devient critique en mouvement : lorsque la bobine se déplace ($x$), l’ordre $\alpha$ change. Recalculer instantanément des dizaines de milliers de coefficients via IFFT dépasse les capacités des DSP les plus puissants. Nous arrivons ici à une limite : le modèle simule le comportement, mais ignore la structure.

3. La Solution Physique : Method of Lines (MoL)

Pour franchir ce mur, Munroe Research privilégie la Method of Lines (MoL). Ici, nous ne cherchons plus à mimer une courbe, mais à discrétiser directement les équations de Maxwell au sein de la géométrie. Pour une pièce polaire axisymétrique, l’entrée du système est le champ magnétique $H_z$ (direction axiale), généré par la bobine, qui se diffuse radialement ($r$) vers l’intérieur du matériau selon l’équation :

De la PDE à l’implémentation FPGA

En découpant le rayon du pôle en anneaux concentriques, on transforme cette équation aux dérivées partielles (PDE) en un système d’équations différentielles ordinaires (ODE).

Discrétisation radiale

Vue spatiale du pôle (n=0 à la surface).

En utilisant un schéma de différence centrale, l’équation de diffusion se discrétise pour chaque point $n$ du rayon :

En regroupant les termes, on obtient une structure de calcul itérative parfaitement adaptée au matériel :

Le flux est injecté à la surface extérieure ($n=0$). À l’axe central ($r=0$, soit $n=N$), nous appliquons une condition de symétrie (dérivée nulle) pour garantir la continuité physique du flux.

Architecture Haute Performance sur Kintex-7

Pour une exécution en temps réel à 96 kHz, nous avons conçu une architecture sur Kintex-7 utilisant le multiplexage temporel. Cette approche permet de simuler 128 points spatiaux en parallèle avec une efficacité de ressources optimale, là où un processeur séquentiel s’effondrerait.

Architecture FPGA

Implémentation multiplexée sur Kintex-7.

Pourquoi c’est une révolution pour le contrôle

Contrairement aux approches précédentes, la MoL est un modèle “White Box”. Sur Kintex-7, nous accédons en temps réel à la densité de flux ($B$) en chaque point du fer. Cela permet d’injecter nativement l’hystérésis et la saturation locale. Si l’extrémité du noyau sature, le modèle réagit instantanément, offrant un contrôle du transducteur d’une précision inégalée, fidèle aux principes physiques fondamentaux.

Synthèse comparative

Conclusion

La Method of Lines sur FPGA représente l’avenir de l’ingénierie acoustique de précision. En remplaçant l’analogie par la simulation physique directe, nous éliminons les compromis entre complexité de calcul et fidélité. Plus de détails techniques sont disponibles dans nos publications JAES et IEEE ci-dessous.

En savoir plus : Découvrez l’implémentation complète de la Method of Lines sur Kintex-7 dans notre étude de cas dédiée.

Références et Publications