Sphère roulante et méthode de l'angle de protection selon la norme CEI 62305

Fondements Physiques des Décharges Atmosphériques et Modèle Électrogéométrique (EGM)

La protection contre les décharges atmosphériques dans les infrastructures complexes repose sur la compréhension rigoureuse du phénomène physique de la rupture diélectrique de l'air à grande échelle. La formation d'un éclair s'initie par un processus d'ionisation graduelle connu sous le nom de traceur descendant échelonné (downward stepped leader), qui transporte une densité de charge considérable depuis le nuage jusqu'à la surface terrestre. À mesure que la tête du traceur descendant s'approche du sol, le champ électrique macroscopique aux pointes des structures élevées excède la rigidité diélectrique critique de l'air impur ou perturbé par les précipitations ( $Ecrit \approx 300 kV/m$ à $3 MV/m$ selon les conditions d'humidité et de pression locale), déclenchant un ou plusieurs traceurs ascendants connecteurs (upward connecting leaders).

Le point dans l'espace tridimensionnel où le traceur descendant et le traceur ascendant s'intersectent définit la phase de saut final ou final jump. À cet instant exact, l'impédance du canal s'effondre et la décharge de retour (return stroke) s'amorce, caractérisée par un front d'onde de courant extrêmement rapide. La distance euclidienne entre la pointe de la structure émettant le traceur ascendant fructueux et la tête du traceur descendant au moment du saut final est définie formellement comme la distance d'impact ou distance de choc ( $r_s$ ).

Le Modèle Électrogéométrique (EGM, Electro-Geometric Model), dérivé à l'origine par Armstrong, Whitehead et Love, et intégré de manière probabiliste dans la série normative IEC 62305, établit une corrélation non linéaire entre l'amplitude de crête du courant de choc ( $I_p$ ) et la distance d'impact ( $r_s$ ). Physiquement, une charge plus importante stockée dans le traceur descendant produit un champ électrique plus sévère à plus grande distance, augmentant la distance à laquelle la rigidité diélectrique de l'espace inter-électrodique succombe. La formulation mathématique standard adoptée pour relier la distance d'impact au courant de crête est :

r_s = A \cdot I_p^b

Où, selon les coefficients empiriques validés de manière opérationnelle par l'IEC 62305 et le CIGRE pour la modélisation des structures terrestres :

r_s = 10 \cdot I_p^{0.65}

Dans cette relation, $r_s$ est exprimé en mètres ( $m$ ) et $I_p$ en kiloampères ( $kA$ ). Le fondement de la conception du système de captage par la sphère roulante réside dans le concept de courant minimal critique ( $Imin$ ). Si un éclair possède un courant inférieur à $Imin$ , sa distance d'impact sera inférieure au rayon de la sphère roulante ( $R$ ). Par conséquent, pour garantir qu'une décharge avec un courant équivalent à $Imin$ soit interceptée avec une certitude géométrique absolue par les dispositifs de captage et non par la structure protégée, le rayon de la sphère roulante doit être fixé rigoureusement en fonction de ce courant minimal :

R = 10 \cdot (Imin)^{0.65} \implies Imin = \left( \frac{R}{10} \right)^{\frac{1}{0.65}}

Par conséquent, un système dimensionné avec un rayon de sphère roulante plus petit offre un niveau de protection supérieur, puisqu'il est capable d'intercepter des éclairs aux impulsions de courant plus faibles (dont les distances de choc sont plus courtes) qui, autrement, pénétreraient le volume d'ombre électrogéométrique.

Paramètres Statistiques et Impulsionnels selon l'IEC 62305-1

La norme IEC 62305-1 classe les niveaux de protection contre la foudre (LPL, Lightning Protection Levels) en quatre catégories (LPL I à LPL IV). Chaque niveau définit un ensemble standardisé de paramètres du courant de foudre, représentant les valeurs limites maximales pour la conception de la résistance mécanique et thermique des composants, et les valeurs limites minimales pour la géométrie d'intersection des dispositifs de captage.

La forme d'onde du courant de retour primaire est modélisée à l'aide de la fonction analytique de double exponentielle ou de l'équation de Heidler, où l'impulsion de première décharge est caractérisée formellement par un rapport de temps de front et de temps de demi-valeur de $10/350\ \mu s$ . L'équation différentielle de l'énergie spécifique transportée par le courant de foudre $i(t)$ à travers une unité de résistance électrique est exprimée par l'intégrale d'action :

\frac{W}{R} = \int0^{\infty} i^2(t) \, dt

Cette énergie spécifique est responsable de l'échauffement adiabatique instantané et des forces électrodynamiques dans les conducteurs de descente. L'élévation de température dans un conducteur de section transversale $A$ soumis à une impulsion de foudre est calculée par l'expression suivante dérivée du bilan thermique adiabatique :

\Delta T = \frac{1}{\alpha} \left[ \exp \left( \frac{\alpha \cdot \rho_0}{C_v \cdot A^2} \int0^{\infty} i^2(t) \, dt \right) - 1 \right]

Où $\alpha$ est le coefficient de température de la résistance électrique ( $K^{-1}$ ), $\rho_0$ est la résistivité à température ambiante ( $\Omega\cdot m$ ), et $C_v$ est la capacité calorifique volumétrique du matériau ( $J /( m ^3\cdot K )$ ).

Paramètre Technique / LPL	LPL I	LPL II	LPL III	LPL IV
Rayon de la Sphère Roulante ( $R$ )	$20 m$	$30 m$	$45 m$	$60 m$
Courant Minimal Intercepté ( $Imin$ )	$3.0 kA$	$5.4 kA$	$10.1 kA$	$15.7 kA$
Courant Maximal de Crête ( $Imax$ , $10/350\ \mu s$ )	$200 kA$	$150 kA$	$100 kA$	$100 kA$
Charge Totale de la Décharge ( $Qflash$ )	$300 C$	$225 C$	$150 C$	$150 C$
Énergie Spécifique ( $W/R$ )	$10.00 MJ /\Omega$	$5.62 MJ /\Omega$	$2.50 MJ /\Omega$	$2.50 MJ /\Omega$
Taille de Maille Réticulaire ( $W \times W$ )	$5 \times 5 m$	$10 \times 10 m$	$15 \times 15 m$	$20 \times 20 m$
Angle de Protection ( $\alpha$ ) pour $h = 20 m$	$25^\circ$	$36^\circ$	$46^\circ$	$54^\circ$
Probabilité d'Interception ( $P_B$ )	$0.99$	$0.98$	$0.95$	$0.80$

Méthode de la Sphère Roulante (Rolling Sphere Method - RSM) : Formulation et Géométrie Tridimensionnelle

La Méthode de la Sphère Roulante est la méthode analytique universelle imposée par la norme IEC 62305-3 pour la détermination des zones protégées sur des structures à géométrie complexe, étant indépendante de la hauteur ou de l'inclinaison de la surface. Géométriquement, elle consiste à faire rouler une sphère imaginaire de rayon $R$ (associé au LPL correspondant) sur et autour de l'ensemble de la structure à protéger, dans toutes les directions possibles, jusqu'à ce qu'elle entre en contact avec le plan du sol ou avec tout objet mis à la terre accessible.

Tous les points de la structure entrant en contact direct avec la surface de la sphère roulante sont considérés comme exposés à des impacts directs de foudre et doivent être équipés d'éléments de captage (pointes de captage, conducteurs aériens ou maillages). Le volume situé sous la surface tangentielle tracée par la sphère entre les points de captage est entièrement protégé dans ce que l'on appelle la zone d'ombre électrogéométrique.

Pénétration de la Sphère entre Dispositifs de Captage

Lorsque la sphère roulante repose sur deux dispositifs de captage parallèles ou isolés (tels que des pointes de captage de hauteur équivalente $h_p$ ou des câbles de garde aériens) séparés par une distance physique $d$ , la sphère pénètre dans l'espace compris entre eux. Pour éviter que la sphère ne touche un quelconque équipement ou élément de la structure situé sous les dispositifs de captage, il est nécessaire de calculer analytiquement la profondeur de pénétration ( $p$ ).

p = R - \sqrt{R^2 - \left( \frac{d}{2} \right)^2}

Cette formule exige intrinsèquement que la distance entre les dispositifs de captage respecte la condition limite de continuité géométrique $d \le 2R$ . Si $d > 2R$ , la sphère passe entre les dispositifs de captage et frappe directement le plan horizontal inférieur ou la structure à protéger.

La hauteur effective de l'enveloppe de protection ( $h_z$ ) au point médian entre deux dispositifs de captage de hauteur $h$ espacés d'une distance $d$ est exprimée par :

h_z = h - p = h - R + \sqrt{R^2 - \left( \frac{d}{2} \right)^2}

Pour des agencements tridimensionnels constitués de quatre pointes de captage disposées aux sommets d'un rectangle de côtés $a$ et $b$ , la distance diagonale représentative est $ddiag = \sqrt{a^2 + b^2}$ . La pénétration tridimensionnelle maximale au centre géométrique du rectangle devient :

p3D = R - \sqrt{R^2 - \frac{a^2 + b^2}{4}}

Condition de non-pénétration destructive pour un objet de hauteur $hobj$ situé dans l'espace intermédiaire :

hobj \le h - p3D \implies hobj \le h - R + \sqrt{R^2 - \frac{a^2 + b^2}{4}}

Méthode de l'Angle de Protection (Protection Angle Method - PAM) et Méthode du Maillage (Mesh Method)

La Méthode de l'Angle de Protection (PAM) est une dérivation simplifiée et bidimensionnelle du Modèle Électrogéométrique. Elle est géométriquement équivalente à la Méthode de la Sphère Roulante uniquement pour des structures simples dotées de limites horizontales simples et de hauteurs restreintes ne dépassant pas le rayon de la sphère roulante correspondante ( $h \le R$ ).

Variation Non Linéaire de l'Angle de Protection

Dans la méthode PAM, le volume protégé par une pointe de captage ou un conducteur aérien adopte la forme d'un cône ou d'un prisme triangulaire dont le sommet coïncide avec l'extrémité du dispositif de captage. L'angle de protection ( $\alpha$ ) n'est pas constant ; il décroît de manière nettement non linéaire avec l'augmentation de la hauteur du dispositif de captage par rapport au plan de référence ( $h$ ). Cette réduction mathématique garantit que le cône formé demeure strictement inscrit à l'intérieur de l'enveloppe circulaire décrite par la sphère roulante de rayon $R$ .

L'équation régissant l'angle de protection à la limite de tangence stricte avec la sphère roulante est :

\alpha(h) = \arcsin \left( 1 - \frac{h}{R} \right)

Le rayon du cercle de base de protection généré au sol ( $rprot$ ) par une pointe de captage à une hauteur $h$ est déterminé analytiquement par :

rprot = h \cdot \tan[\alpha(h)]

En substituant la dépendance implicite de $\alpha$ en fonction de la hauteur géométrique et du rayon de la sphère roulante $R$ :

rprot = h \cdot \tan \left[ \arcsin \left( 1 - \frac{h}{R} \right) \right] = h \cdot \frac{1 - \frac{h}{R}}{\sqrt{1 - \left(1 - \frac{h}{R}\right)^2}} = h \cdot \frac{R - h}{\sqrt{2Rh - h^2}}

Cette expression démontre rigoureusement que lorsque la hauteur du dispositif de captage s'approche du rayon de la sphère roulante ( $h \to R$ ), l'angle de protection tend vers zéro ( $\alpha \to 0^\circ$ ), annulant ainsi l'applicabilité de la méthode de l'angle de protection pour des hauteurs $h > R$ .

Méthode du Maillage (Mesh Method)

Pour les toitures plates de grande extension ou les structures industrielles où l'application de pointes ou de câbles aériens s'avère techniquement infaisable, la norme IEC 62305-3 prescrit la Méthode du Maillage. Un réseau réticulaire de conducteurs disposés sur la surface du toit assure une protection intégrale si trois conditions géométriques infranchissables sont simultanément remplies :

Les conducteurs de captage sont installés sur les bords de la toiture, sur les parties saillantes de la structure et sur les lignes de faîtage.
Les dimensions des mailles (largeur du réseau $W$ ) ne dépassent pas les valeurs limites définies dans le tableau normatif selon le LPL (de $5 \times 5 m$ pour le LPL I jusqu'à $20 \times 20 m$ pour le LPL IV).
Aucune installation métallique située sur la toiture ne dépasse du volume protégé par le réseau conducteur sans l'aide d'un dispositif de captage auxiliaire.

Analyse Forensique des Défaillances des Systèmes de Captage et des Diélectriques Industriels

Les défaillances de la protection contre les décharges atmosphériques trouvent leur origine dans une sous-estimation du LPL, une modélisation électrogéométrique défectueuse ou un calcul insuffisant des distances d'isolement diélectrique. L'analyse forensique de ces défaillances englobe trois mécanismes physiques principaux :

Impacts Latéraux sur les Structures de Grande Hauteur

Dans les structures dont la hauteur dépasse le rayon de la sphère roulante ( $h > R$ ), la sphère de rayon $R$ peut toucher les parois latérales verticales en des points situés à une hauteur supérieure à $hlim = R$ . Le courant des décharges sur les parois latérales est généralement plus faible, mais il provoque le détachement de maçonnerie, la perforation de façades ventilées en aluminium et la destruction des capteurs extérieurs. L'erreur forensique réside dans l'omission d'étendre l'anneau de captage et les conducteurs de descente aux bords et aux parties saillantes du tiers supérieur de la structure dans les bâtiments de grande hauteur (où $h > 60 m$ selon l'IEC 62305-3).

Claquage en Retour (Back-Flashover) par Inductance du Conducteur de Descente

Quand le courant de foudre $i(t)$ circule à travers un conducteur de descente, la chute de tension totale le long du conducteur par rapport à la masse profonde de la terre ne dépend pas uniquement de la résistance de mise à la terre ( $Rpat$ ), mais aussi de l'inductance distribuée du conducteur de descente ( $Lbaj \approx 1\ \mu H/m$ ) en raison de la valeur élevée de la dérivée temporelle du courant ( $di/dt$ ) :

V(t) = i(t) \cdot Rpat + Lbaj \cdot \frac{di(t)}{dt}

Pour un éclair de niveau LPL I de $200 kA$ avec un temps de front de $10\ \mu s$ , la dérivée initiale moyenne est $\frac{di}{dt} = 20 kA /\mu s = 2 \times 10^{10} A/s$ . Dans un conducteur de descente de $20 m$ de longueur ( $Lbaj = 20\ \mu H$ ), la composante inductive génère un pic de tension transitoire impulsionnelle de :

Vind = 20 \times 10^{-6} H \cdot 2 \times 10^{10} A/s = 400 kV

Si la rigidité diélectrique de l'air ( $\sim 3 MV/m$ ) ou des isolants solides séparant le conducteur de descente des tuyauteries ou du câblage interne n'est pas suffisante pour supporter cette différence de potentiel, un arc électrique destructeur (claquage en retour ou back-flashover) se produit, provoquant des incendies, la perforation de conduites de gaz ou l'injection directe de surtensions brutales dans les réseaux basse tension.

Déformation Thermique et Fusion Diélectrique du Conducteur

La contrainte thermique au point d'impact direct entraîne la fusion du métal par vaporisation ponctuelle. La masse de matériau fondu ( $mfund$ ) au point d'entrée du courant continu secondaire ou de l'impulsion est liée à la charge totale $Q$ par le potentiel de chute cathodique/anodique ( $U_{c,a} \approx 10 V$ ) :

mfund = \frac{U_{c,a} \cdot Q}{L_f + C_p(Tfund - Tamb)}

Où $L_f$ est la chaleur latente de fusion et $C_p$ la chaleur spécifique du conducteur. Une conception forensique déficiente omet l'épaisseur minimale des tôles métalliques continues (définie dans l'IEC 62305-3, Tableau 3), ce qui provoque la perforation directe des réservoirs de stockage d'hydrocarbures.

Stratégies de Conception Avancée et Atténuation des Surtensions Transitoires

Pour garantir l'intégrité structurelle et fonctionnelle des installations complexes, la conception du Système de Protection contre la Foudre (SPCR / LPS) doit mettre en œuvre des topologies isolées ou non isolées évaluées rigoureusement par le biais de la distance de séparation ( $s$ ).

Calcul de la Distance de Séparation Normative

La distance de séparation $s$ représente la distance de sécurité air/matériau requise entre le système de captage ou de descente et les composants métalliques ou les lignes électriques internes, afin de prévenir le phénomène de claquage en retour. La formulation générale selon la norme IEC 62305-3 est :

s = k_i \cdot \frac{k_c}{k_m} \cdot l

Où :

$k_i$ : Coefficient qui dépend du Niveau de Protection contre la Foudre (LPL) sélectionné :
$k_i = 0.08 \quad ( LPL I ), \quad k_i = 0.06 \quad ( LPL II ), \quad k_i = 0.04 \quad ( LPL III / IV )$
$k_m$ : Coefficient d'isolement électrique du milieu où l'on évalue la séparation :
$k_m = 1.0 \quad ( Air ), \quad k_m = 0.5 \quad ( Béton / Maçonnerie ), \quad k_m = 0.7 \quad ( Isolant synthétique / PRV )$
$k_c$ : Coefficient de division du courant de foudre entre les conducteurs de descente. Pour un agencement spatial de $N$ conducteurs de descente interconnectés dans un système non isolé :
$k_c = \frac{1}{2N} \quad ou dérivé de la matrice d'impédances de réseau.$
$l$ : Distance géométrique linéaire le long du conducteur de descente ou du dispositif de captage depuis le point où l'on évalue la distance de séparation jusqu'au point d'équipotentialité le plus proche (le point de mise à la terre).

Coordination des Parasurtenseurs (SPD) dans les Zones LPZ

L'atténuation des effets induits électromagnétiquement nécessite le découpage du volume en Zones de Protection contre la Foudre (LPZ, Lightning Protection Zones) selon la norme IEC 62305-4 :

LPZ 0A : Zone exposée aux impacts directs et aux champs électromagnétiques non atténués.
LPZ 0B : Zone protégée contre les impacts directs par la méthode de la sphère roulante, mais exposée aux champs électromagnétiques complets.
LPZ 1 : Zone intérieure où le courant de choc est limité par la division des courants et par la présence de parafoudres de Type 1 (testés avec une onde $10/350\ \mu s$ ).
LPZ 2 : Zone intérieure avec blindage additionnel où les surtensions transitoires sont réduites davantage par des parafoudres de Type 2 (testés avec une onde $8/20\ \mu s$ ).

Le couplage inductif dans les boucles formées par les câbles de signal et d'énergie à l'intérieur de la structure est régi par la loi de Faraday-Lenz. La tension induite dans une boucle de surface $Aboucle$ soumise à un champ magnétique variable $B(t)$ couplé au courant du conducteur de descente $i(t)$ est :

Vind(t) = -\frac{d\Phi}{dt} = -\frac{d}{dt} \int \mathbf{B}(t) \cdot d\mathbf{A} \approx \frac{\mu_0 \cdot Aboucle}{2\pi \cdot d_p} \cdot \frac{di(t)}{dt}

Par conséquent, la réduction de la surface de la boucle ( $Aboucle$ ) au moyen d'un cheminement parallèle torsadé ou d'un blindage coaxial s'avère tout aussi indispensable que l'installation du parafoudre lui-même.

Intégration Technique et Simulation Électromagnétique dans Vexten Suite

Dans le domaine de l'ingénierie électromagnétique moderne, le logiciel Vexten Suite agit comme le moteur central d'analyse déterministe et stochastique pour la modélisation 3D des systèmes de captage contre les décharges atmosphériques. Vexten Suite automatise les méthodologies de calcul spécifiées par la norme IEC 62305 grâce à une intégration analytique entre la distribution spatiale du courant et la réponse en fréquence du réseau diélectrique.

Algorithme de Tracé de Foudre et Calcul Tridimensionnel de la Sphère Roulante

Le moteur de calcul géométrique de Vexten Suite projette numériquement une surface discrète hypersphérique sur la structure importée à partir de formats BIM/CAD (STEP, IFC). Le système résout l'interaction entre le maillage de la structure et la sphère roulante au moyen d'équations polynomiales d'intersection sphère-surface :

\mathbf{P}_{centre} = \mathbf{P}_{capteur} + R \cdot \hat{\mathbf{n}}

Où $\hat{\mathbf{n}}$ est le vecteur normal à la surface de contact dans l'espace tridimensionnel. Le programme effectue un balayage angulaire infinitésimal ( $d\theta, d\phi$ ) en enregistrant tous les points de pénétration où la distance entre la sphère et l'enveloppe de la structure est inférieure à zéro ( $dreal - R < 0$ ). Les points vulnérables sont immédiatement signalés vectoriellement avec leur Niveau de Protection contre la Foudre (LPL) respectif.

Calcul de l'Impédance Impulsionnelle et Intégration avec les Réseaux de Mise à la Terre

Vexten Suite intègre simultanément l'analyse des courants de court-circuit à fréquence industrielle (IEC 60909) et l'analyse transitoire impulsionnelle à haute fréquence de la décharge atmosphérique (IEC 62305-3). Le calcul de l'impédance transitoire du réseau de mise à la terre ( $Z(t)$ ) ne se limite pas à la valeur de résistance en courant continu ( $Rdc$ ), mais prend en compte le comportement inductif des piquets et conducteurs enterrés à travers la modélisation de lignes de transmission distribuées :

Z(t) = \mathcal{F}^{-1} \left\{ \sqrt{\frac{R + j\omega L}{G + j\omega C}} \cdot \coth\left( \sqrt{(R + j\omega L)(G + j\omega C)} \cdot l_e \right) \right\}

Où $l_e$ est la longueur effective critique de l'électrode de mise à la terre sous impulsions à haute fréquence. Si la longueur physique du conducteur enterré dépasse $l_e$ , toute longueur additionnelle ne contribue pas à la réduction de la surtension initiale de la décharge. Vexten Suite calcule la longueur effective critique au moyen de la relation empirique :

l_e = 1.6 \cdot \sqrt{\rho_{sol} \cdot T_1}

Où $\rho_{sol}$ représente la résistivité du terrain ( $\Omega\cdot m$ ) et $T_1$ le temps de front de l'impulsion de courant ( $\mu s$ ).

Flux Intégré de Conception et de Validation dans Vexten Suite

Évaluation Automatique du Risque (IEC 62305-2) : Détermination de la fréquence des impacts annuels ( $N_d$ ) et calcul des pertes tolérables ( $R_1$ à $R_4$ ) pour fixer par calcul le niveau de protection requis (LPL I à IV).
Génération Géométrique de la Zone Protégée : Exécution du module de la Sphère Roulante en 3D. Le logiciel détermine la position précise des pointes de captage, la hauteur des mâts ou le tracé des réseaux maillés sur les toits afin d'éviter les points de collision directe.
Évaluation des Distances de Séparation ( $s$ ) : Cartographie nodale en 3D de la chute de tension dans chaque conducteur de descente et vérification automatique de la règle de séparation isolante par rapport aux chemins de câbles, tuyauteries et gaines de process.
Coordination Dynamique des SPD : Génération des spécifications des parafoudres de Type 1, Type 2 et Type 3 requis à chaque frontière de zone LPZ, garantissant la coordination de l'énergie basée sur la dissipation du paramètre d'énergie spécifique ( $W/R$ ).

Grâce à la rigueur mathématique de l'algorithme de Vexten Suite, le concepteur obtient un modèle électromagnétique déterministe qui élimine la méthode empirique des essais et erreurs, garantissant une protection intégrale contre les transitoires diélectriques, mécaniques et thermiques sévères provoqués par les décharges atmosphériques, conformément aux normes internationales les plus exigeantes de la série IEC 62305.