Guide utilisateur — Maintenance Turbines LM2500

Introduction

🛠️ À propos de cette application

Cette application est un outil d'aide à la décision pour la maintenance des turbines à gaz de type LM2500. Elle regroupe cinq modules complémentaires qui couvrent l'ensemble du cycle de vie d'une décision de maintenance : du planning prévisionnel à la justification économique des actions retenues.

Module	Question à laquelle il répond	Méthode
Planning turbines	Quand sont mes prochaines révisions L2 / L3 ?	Calcul par paliers horaires
Vue turbogénérateurs	Quel est l'état de ma flotte aujourd'hui ?	Tableau de bord temps réel
RCM II	Quelle maintenance appliquer sur chaque mode de défaillance ?	Reliability Centered Maintenance
Analyse Weibull	Quel est le comportement de défaillance de mon équipement ?	Régression Weibull sur données de retour d'expérience
Weibull Thermique	Comment la température modifie-t-elle la durée de vie ?	Modèles Arrhenius (turbines) et Montsinger (bobinages)
BCR & Isorisque	Quelle action de maintenance est la plus rentable ?	Benefit-to-Cost Ratio + courbe €/RU
Compression / Débitmètre	Quelle puissance pour comprimer ce gaz, et quel débit dans ma conduite ?	Propriétés gaz réel (CoolProp), méthode de Schultz, ISO 5167

Section 1

🟥 La matrice de risque

La matrice de risque est l'outil central de l'application. Elle permet d'évaluer et de visualiser le niveau de risque associé à chaque mode de défaillance, en combinant deux dimensions indépendantes : la probabilité d'occurrence et la conséquence en cas de défaillance.

📏 Définition de P et de C

Chaque axe est gradué de 1 à 6. Ce choix (plutôt que 1-5 ou 1-10) donne une matrice 6×6 avec 36 cases, suffisamment fine pour différencier les risques sans créer une fausse précision.

Niveau P	Label	Fréquence annuelle	Exemple LM2500
1	Quasi impossible	< 0,01 /an (1 fois en 100 ans)	Rupture du carter de turbine
2	Très rare	0,01–0,1 /an (1 fois en 10–100 ans)	Rupture d'aube en titane
3	Rare	0,1–0,5 /an (1 fois en 2–10 ans)	Fuite huile sur palier principal
4	Occasionnel	0,5–2 /an (1 à 2 fois par an)	Encrassement filtre huile
5	Fréquent	2–12 /an (mensuel à bimensuel)	Alarme vibration capteur
6	Très fréquent	> 12 /an (hebdomadaire)	Encrassement filtre air compresseur

Niveau C	Label	Impact opérationnel	Ordre de grandeur coût
1	Négligeable	Aucun arrêt, aucune blessure	< 1 k€
2	Mineure	Arrêt < 4 h	1–10 k€
3	Modérée	Arrêt 4–24 h, blessure légère possible	10–100 k€
4	Sérieuse	Arrêt 1–7 j, blessure grave possible	100 k€–1 M€
5	Critique	Arrêt > 7 j, invalidité possible	> 1 M€
6	Catastrophique	Perte équipement, décès possible	Perte totale

🗺️ Les 4 zones de risque

Le score de risque est simplement le produit R = P × C. Il va de 1 (minimum, P=1 et C=1) à 36 (maximum, P=6 et C=6). On définit 4 zones à partir de seuils :

R = P × C

Faible : R ∈ [1, 6] → RTF acceptable, surveillance légère
Modéré : R ∈ [7, 14] → tâche préventive souhaitable
Élevé : R ∈ [15, 24] → tâche préventive obligatoire
Critique : R ∈ [25, 36] → action immédiate requise

⚠️ Attention aux diagonales : deux modes de défaillance avec le même score peuvent nécessiter des stratégies très différentes. Un score de 12 peut venir de (P=2, C=6) — rare mais catastrophique — ou de (P=6, C=2) — très fréquent mais bénin. La stratégie de maintenance sera radicalement différente. Lisez toujours P et C séparément avant de regarder le score global.

💡 Exemple : compresseur axial LM2500

✈️ Cas réel — Rupture d'aube compresseur

Mode de défaillance : Rupture d'aube due à la fatigue matière ou à l'ingestion d'un corps étranger (FOD — Foreign Object Damage).

Évaluation :
— Probabilité : 2 — Très rare (environ 1 fois en 20 ans par turbine)
— Conséquence : 6 — Catastrophique (arrêt immédiat, risque d'incendie, dommages en cascade)

Score : R = 2 × 6 = 12 — Zone Modérée.
Mais ne vous fiez pas au score seul : la conséquence catastrophique (C=6) impose une tâche de type surveillance sur condition (TD) malgré la probabilité faible. Une inspection boroscopique trimestrielle est justifiée.

Section 2

🔧 RCM II — Reliability Centered Maintenance

La RCM (Maintenance Centrée sur la Fiabilité) est une méthode structurée issue de l'aviation civile (Nowlan & Heap, 1978, United Airlines) et formalisée en RCM II par John Moubray. Son objectif : pour chaque mode de défaillance, identifier la tâche de maintenance la plus efficace compte tenu du comportement de défaillance et des conséquences.

🧠 Les 7 questions fondamentales de la RCM

1

Quelles sont les fonctions de l'équipement ?
Ex : "Comprimer l'air de 1 bar à 23 bar avec un débit de 86 kg/s"

2

Comment peut-il ne plus remplir sa fonction ?
Ex : "Chute progressive du taux de compression en dessous de 20 bar"

3

Qu'est-ce qui cause cette défaillance fonctionnelle ?
Ex : "Encrassement des aubes par dépôts salins"

4

Que se passe-t-il quand la défaillance survient ?
Ex : "Augmentation de la consommation de carburant, perte de puissance"

5

En quoi cela a-t-il de l'importance ?
Ex : "Impact opérationnel majeur, risque de déclenchement sur survitesse"

6

Que peut-on faire pour prévenir ou anticiper la défaillance ?
Ex : "Surveillance hebdomadaire de la pression différentielle compresseur"

7

Que faire si aucune tâche préventive n'est applicable ?
Ex : "Accepter la panne (RTF) avec stock de pièces de rechange"

📋 Les 5 tâches RCM

Code	Nom complet	Principe	Quand l'utiliser
TD	Tâche sur condition	Surveiller un paramètre jusqu'à un seuil d'alarme	Quand la dégradation est détectable avant la panne (courbe P-F)
TRP	Restauration préventive	Révision complète ou partielle à intervalle fixe	Quand il existe un âge limite ET que la restauration remet à neuf
TRE	Remplacement préventif	Remplacement systématique à intervalle fixe	Quand il existe un âge limite ET que la pièce n'est pas réparable
RTF	Run-To-Failure	Aucune action — on attend la panne	Quand le coût de la maintenance dépasse le coût de la panne, ou risque faible
RDT	Recherche de défaillance	Test périodique pour révéler une défaillance cachée	Pour les équipements de protection dont la défaillance n'est pas évidente

🌳 Logique de sélection de la tâche

La sélection de la bonne tâche suit une logique séquentielle. Les questions se posent dans cet ordre, et la première réponse positive détermine la tâche :

1. La défaillance est-elle évidente pour l'opérateur ? → Non → Considérer RDT (recherche de défaillance cachée)

2. La dégradation est-elle détectable AVANT la défaillance ? → Oui → TD (surveillance sur condition)

3. Existe-t-il un âge limite au-delà duquel le taux de panne augmente ? → Oui + pièce réparable → TRP → Oui + pièce non réparable → TRE

4. Aucune des tâches ci-dessus n'est applicable ou rentable → RTF (si risque acceptable) → Action de conception si risque critique

💡 Exemple — Filtre huile LM2500

Mode de défaillance : Colmatage du filtre huile principal.

Analyse : La dégradation (colmatage) est progressive et détectable via la mesure de pression différentielle. Il n'y a pas d'âge limite strict car le colmatage dépend de la qualité de l'huile et des conditions d'exploitation.

Tâche retenue : TD — Mesure quotidienne de la pression différentielle. Remplacement du filtre si ΔP > 2,5 bar. Simple, efficace, pas de remplacement systématique coûteux.

Section 3

📊 Analyse de Weibull

L'analyse de Weibull est une méthode statistique qui permet de modéliser le comportement de défaillance d'un équipement à partir de données de retour d'expérience (durées de vie observées). Elle répond à la question : "Mon équipement a-t-il un âge limite ?"

📐 La loi de Weibull à 2 paramètres

La loi de Weibull est caractérisée par deux paramètres : β (bêta) le paramètre de forme, et η (êta) le paramètre d'échelle (durée caractéristique, en heures ou cycles).

// Fonction de défaillance cumulée F(t) F(t) = 1 − exp(−(t/η)^β)

// Fiabilité (probabilité de survie) S(t) S(t) = 1 − F(t) = exp(−(t/η)^β)

// Taux de défaillance instantané h(t) h(t) = (β/η) × (t/η)^(β−1)

// Durée de vie caractéristique η À t = η : F(η) = 1 − e^(−1) ≈ 63,2% des équipements ont défailli

🔎 Comment interpréter le paramètre β (bêta)

Le paramètre β est le plus important à lire en premier. Il vous dit pourquoi les équipements tombent en panne :

Valeur de β	Régime	Taux de panne	Cause typique	Implication maintenance
β < 1	Infantile (mortalité juvénile)	Décroissant avec l'âge	Défauts de fabrication, mauvais montage, rodage insuffisant	TRP/TRE contre-productives — préférer TD + contrôle qualité
β ≈ 1	Aléatoire (vie utile)	Constant (loi exponentielle)	Causes externes, surcharges, FOD, erreurs humaines	TRP/TRE inefficaces — seules TD ou RTF sont pertinentes
β > 1	Usure (vieillissement)	Croissant avec l'âge	Fatigue, corrosion, abrasion, usure mécanique	TRP/TRE justifiées — il existe un âge optimal de remplacement

💡 Exemple — Roulements de palier LM2500

Données : 12 roulements remplacés au cours des 5 dernières années. Durées de fonctionnement avant remplacement : 18 000h, 22 000h, 19 500h, 21 000h, 24 000h, 17 500h, 23 000h, 20 000h, 25 000h, 18 500h, 22 500h, 21 500h.

Résultat Weibull : β = 3,2 — η = 21 500h.
Interprétation : β > 1 confirme un mécanisme d'usure progressive. Le taux de panne augmente fortement après 21 500h. À 63,2% des équipements ont déjà défailli à 21 500h.

Décision : TRP justifiée. Remplacement préventif à 18 000h (avant le genou de la courbe) pour maintenir une fiabilité > 90%.

Section 3b

🌡️ Weibull avec correction thermique

La loi de Weibull standard suppose que les conditions d'exploitation sont constantes. En pratique, la température joue un rôle majeur : un équipement qui tourne plus chaud que prévu vieillit beaucoup plus vite. Le module Weibull Thermique permet de quantifier cet effet en recalculant le paramètre η en fonction de la température, tout en conservant β inchangé.

Deux modèles couvrent les cas les plus courants en maintenance industrielle :

🔥 Modèle 1 — Arrhenius-Weibull

Le modèle d'Arrhenius repose sur un principe physico-chimique : la chaleur accélère les réactions de dégradation. Il est utilisé pour les composants mécaniques et thermiques soumis à des variations de charge (turbines, injecteurs, paliers).

// η varie avec la température absolue T (en Kelvin = °C + 273.15) η(T) = C × exp( Ea / (k × T) )

// Paramètres C = constante matériau / composant (calibrée sur les données) Ea = énergie d'activation (eV) ← résistance du matériau à la chaleur k = constante de Boltzmann = 8.617 × 10⁻⁵ eV/K T = température en Kelvin (°C + 273.15)

// Facteur d'accélération entre deux températures AF = exp( Ea/k × (1/T_use − 1/T_stress) )

// Lien avec η : si AF = 3 alors la machine vieillit 3× plus vite à T_stress η(T_stress) = η(T_use) / AF

Ea (eV)	Mécanisme de dégradation typique	Composants concernés
0.3 – 0.5	Dégradation mécanique légère, fatigue	Pièces structurelles acier, engrenages
0.6 – 0.8	Oxydation, corrosion à chaud	Aubes turbine, injecteurs, tuyauteries
0.9 – 1.2	Dégradation de polymères, joints, revêtements	Joints d'étanchéité, câblage, flexibles
1.2 – 1.4	Défauts semiconducteurs, contamination	Électronique de commande, capteurs

⚠️ Attention : La relation est exponentielle, pas linéaire. +50 °C ne double pas le vieillissement : selon Ea, il peut le multiplier par 3, 5 ou 10. C'est pourquoi les petites dérives de température d'exploitation ont des conséquences bien plus grandes qu'on ne l'intuit.

⚡ Modèle 2 — Montsinger-Weibull (bobinages moteurs)

La loi de Montsinger (1930) est une règle empirique fondamentale en électrotechnique, normalisée dans la norme IEC 60085. Elle s'applique aux isolants électriques des bobinages de moteurs et générateurs.

// Durée de vie de l'isolant en fonction de la température L(T) = L₀ × exp( −α × (T − T₀) )

// Paramètres L₀ = durée de vie nominale à T₀ (heures) T₀ = température nominale de la classe d'isolation (°C) α = coefficient de dégradation = ln(2)/10 ≈ 0.0693 (règle des 10°C) T = température réelle du bobinage (°C)

// Règle de Montsinger — la plus simple à retenir +10 °C → durée de vie divisée par 2 −10 °C → durée de vie multipliée par 2

Classe d'isolation	T_max (°C)	Échauffement admissible	Application typique
Classe A	105 °C	60 K	Moteurs anciens, bobinages coton / soie
Classe E	120 °C	75 K	Moteurs standards
Classe B	130 °C	80 K	Moteurs industriels courants ← le plus fréquent
Classe F	155 °C	105 K	Moteurs à variateur de vitesse, applications sévères
Classe H	180 °C	125 K	Moteurs haute performance, traction

💡 Exemples chiffrés

🔥 Exemple 1 — Turbine LM2500 (Arrhenius)

Données connues : β = 2,5 et η = 8 000 h estimés sur des données d'injecteurs à une température d'échappement de 450 °C. Ea = 0,7 eV (oxydation).

Question : La turbine doit tourner à 520 °C lors d'une surcharge. Quel est l'impact sur la durée de vie ?

Calcul :
AF = exp( 0,7 / 8,617×10⁻⁵ × (1/723,15 − 1/793,15) ) ≈ 2,70
η à 520 °C = 8 000 / 2,70 ≈ 2 968 h (au lieu de 8 000 h)
MTTF passe de 7 098 h → 2 634 h · B10 passe de 3 252 h → 1 207 h

Conclusion : une surcharge de +70 °C réduit l'espérance de vie de 63 %. La prochaine inspection préventive doit être avancée en conséquence.

⚡ Exemple 2 — Bobinage moteur classe B (Montsinger)

Données connues : Moteur asynchrone classe B (T_max = 130 °C). Durée de vie nominale estimée à 25 000 h. β = 2,0.

Question : Une surcharge fait monter le bobinage à 150 °C. Quelle est la nouvelle durée de vie ?

Calcul :
ΔT = 150 − 130 = +20 °C → 2 fois la règle des 10°C
Facteur = (1/2)² = 0,25
η à 150 °C = 25 000 × 0,25 = 6 250 h (au lieu de 25 000 h)
MTTF passe de 22 156 h → 5 539 h

Conclusion : +20 °C divise la durée de vie par 4. Un moteur prévu pour durer 25 000 h ne tiendra que 6 250 h si la surcharge est permanente. Surveiller la température bobinage est ici aussi important que surveiller les vibrations.

🧭 Quel modèle choisir ?

Critère	Arrhenius-Weibull	Montsinger-Weibull
Application principale	Composants mécaniques / thermiques	Isolants électriques
Exemples	Turbines, injecteurs, paliers, joints	Bobinages moteur, transformateurs
Paramètre clé	Ea — énergie d'activation (eV)	α ou règle des 10°C (IEC 60085)
Fondement	Physico-chimique (universel)	Empirique (électrotechnique)
Température en entrée	Kelvin obligatoire (°C + 273,15)	Celsius (écart ΔT)

Section 4

💰 BCR — Benefit-to-Cost Ratio

Une fois que la RCM a identifié les actions de maintenance possibles, il faut les prioriser et les justifier économiquement. Le BCR mesure le retour sur investissement d'une action de maintenance : combien de pertes évite-t-on pour chaque euro investi ?

📐 Formule du BCR

// Formule générale BCR = Bénéfice_annuel / Coût_action_annuel

// Calcul du bénéfice (réduction de risque × coût d'une panne) Bénéfice = (ΔRisque / 36) × Coût_panne

// Réduction de risque (sur une échelle de 0 à 36) ΔRisque = (P_avant × C_avant) − (P_après × C_après)

// Interprétation du BCR BCR < 1 → L'action coûte plus qu'elle ne rapporte → Rejeter BCR 1–2 → Rentabilité limite → Étudier BCR 2–5 → Rentable → Retenir BCR > 5 → Très rentable → Prioriser

⚙️ Les deux leviers d'une action de maintenance

Une action de maintenance peut réduire le risque de deux façons indépendantes — sur la probabilité, sur la conséquence, ou sur les deux. C'est une nuance importante que beaucoup d'analyses simplifient à tort.

Levier	Définition	Exemples d'actions
⬇️ Réduire P	Diminuer la probabilité que la défaillance survienne	Inspection vibratoire (détecte avant la panne), remplacement préventif d'usure, nettoyage préventif
⬇️ Réduire C	Diminuer la gravité des conséquences si la défaillance survient quand même	Stock de pièces de rechange (réduit l'arrêt de 7 jours à 4h), formation des opérateurs (réaction plus rapide), équipement de secours
⬇️ Réduire P + C	Agir simultanément sur les deux dimensions	Redondance d'équipement (bascule automatique), modernisation technologique

💡 Exemple chiffré — Système de lubrification LM2500

Situation : Pompe huile principale — risque de pression basse.
P = 3 (Rare), C = 5 (Critique) → Score = 15, zone Élevée.
Coût estimé d'une panne : 300 000 €.

Action A — Inspection hebdomadaire de la pression :
Coût annuel : 3 000 €. Réduit P de 3 → 1 (détection précoce). C reste à 5.
ΔRisque = (3×5) − (1×5) = 15 − 5 = 10 points
Bénéfice = (10/36) × 300 000 = 83 333 €/an
BCR = 83 333 / 3 000 = BCR = 27,8 🚀

Action B — Stock pompe de rechange :
Coût annuel : 8 000 € (amortissement). Réduit C de 5 → 2 (réparation en 4h). P reste à 3.
ΔRisque = (3×5) − (3×2) = 15 − 6 = 9 points
Bénéfice = (9/36) × 300 000 = 75 000 €/an
BCR = 75 000 / 8 000 = BCR = 9,4 ✅

Conclusion : Les deux actions sont rentables. L'inspection (A) est plus efficace par euro investi. Idéalement, on combine les deux : l'inspection détecte, le stock de pièces limite les pertes si la panne survient quand même.

Section 5

📈 Courbe €/RU — Coût par Risk Unit

Le BCR permet de savoir si une action est rentable. La courbe €/RU répond à une question différente : "Dans quel ordre doit-on investir ?" et "Jusqu'où vaut-il la peine d'aller ?".

📐 Formule du €/RU

// Coût par unité de risque réduite €/RU = Coût_action / ΔRisque

// Exemple : action qui coûte 5 000€ et réduit le risque de 10 points €/RU = 5 000 / 10 = 500 €/RU

// Interprétation €/RU bas → action très efficace par point de risque éliminé → à prioriser €/RU élevé → action chère pour peu de résultat → à budgétiser avec précaution

📉 Le coude de la courbe — seuil d'arbitrage

Quand on trie toutes les actions par €/RU croissant et qu'on les empile, on obtient une courbe en escalier. Cette courbe monte progressivement : les premières actions sont peu chères par RU, les suivantes de plus en plus chères.

Le coude est le point où la pente de la courbe change brutalement — c'est-à-dire l'action à partir de laquelle chaque point de risque en moins devient significativement plus cher. C'est votre seuil d'arbitrage budgétaire : les actions avant le coude sont à retenir en priorité ; les actions après le coude doivent être soumises à une analyse coût-bénéfice plus fine.

💡 Exemple chiffré — Compresseur axial LM2500

Situation de départ : P = 4, C = 5, score = 20 (Élevé). Coût d'une panne estimé : 500 000 €.

Trois actions sont envisagées :

Action	Coût/an	P après	C après	ΔRisque	€/RU	BCR
Inspection thermographie IR	5 000 €	2	5	10 pts	500 €/RU	BCR 27,8
Stock pièce de rechange	15 000 €	4	2	12 pts	1 250 €/RU	BCR 5,6
Redondance équipement	80 000 €	1	2	18 pts	4 444 €/RU	BCR 3,1

Lecture de la courbe en escalier :
— Marche 1 (inspection IR) : 10 points de risque éliminés à 500 €/RU
— Marche 2 (stock pièce) : 12 points supplémentaires à 1 250 €/RU
— Marche 3 (redondance) : 18 points supplémentaires à 4 444 €/RU

Coude détecté entre la marche 2 et la marche 3 : le prix par point de risque est multiplié par 3,5. Si le budget est contraint, on retient l'inspection IR + le stock de pièces (40 points de risque éliminés, ΔRisque cumulé = 22 pts sur 36 possibles, score final = 20 − 22 = −2 → score plancher à 2, zone Faible).

⚠️ Attention : Le €/RU ne remplace pas le BCR, il le complète. Une action peut avoir un excellent €/RU (peu chère par RU) mais un BCR < 1 si le coût de la panne est très faible. Utilisez toujours les deux indicateurs ensemble.

Section 5b

🗜️ Compression & Débitmètre

Ce module sort de la logique « risque / maintenance » des sections précédentes : il s'agit d'un outil de thermodynamique appliquée. Il répond à deux questions concrètes du métier : « quelle puissance faut-il pour comprimer ce gaz ? » et « quel débit circule dans ma conduite ? ». Les propriétés des gaz réels sont calculées par la bibliothèque CoolProp (équations d'état de référence), à partir d'un mélange que vous composez librement.

🧪 Gaz parfait, gaz réel et facteur Z

La loi des gaz parfaits (P·v = r·T) suppose que les molécules n'interagissent pas. À basse pression c'est une bonne approximation, mais dès qu'on comprime, les molécules se rapprochent et le gaz s'écarte du comportement idéal. On mesure cet écart par le facteur de compressibilité Z :

// Équation d'état réelle (Z corrige la loi des gaz parfaits) P · v = Z · r · T

// r = constante spécifique du mélange = R / M Z = 1 → gaz parfait Z < 1 → gaz plus compressible que parfait (cas courant en compression) Z > 1 → gaz moins compressible (très hautes pressions)

Pour un gaz naturel typique vers 10 bar, Z vaut environ 0,98 : l'écart au gaz parfait est déjà de 2 %. À 40–100 bar il devient bien plus marqué. Ignorer Z conduit à des erreurs de plusieurs pour-cent sur la puissance — d'où la nécessité d'une méthode « gaz réel » comme celle de Schultz.

📐 La méthode de Schultz (ASME PTC-10)

Une compression réelle n'est ni parfaitement isotherme ni adiabatique réversible : c'est une transformation polytropique, avec des pertes. Le rendement polytropique η_p mesure la part « utile » du travail fourni :

// Définition du rendement polytropique (le long du trajet réel) η_p = (travail réversible) / (travail réel) = (v·dP) / dh

// Hauteur (head) = énergie par kg de gaz, en J/kg Hauteur polytropique H_p = η_p · (h₂ − h₁) Hauteur isentropique H_s = h₂ₛ − h₁ (compression idéale réversible)

// Rendement isentropique (autre façon de noter la performance) η_s = H_s / (h₂ − h₁)

// Puissance absorbée par le gaz P_gaz = débit massique · (h₂ − h₁)

Le problème : sur un gaz réel, les formules polytropiques simples (établies pour le gaz parfait) sont fausses. R. Schultz (1962) a proposé une correction, aujourd'hui au cœur de la norme d'essai ASME PTC-10. Elle repose sur deux coefficients de compressibilité X et Y qui mesurent la sensibilité du volume du gaz à la température et à la pression :

// Coefficients de Schultz (nuls/unitaires pour un gaz parfait) X = (T/v)·(∂v/∂T)_P − 1 → vaut 0 pour un gaz parfait Y = −(P/v)·(∂v/∂P)_T → vaut 1 pour un gaz parfait

// Exposant polytropique réel (déduit des volumes d'entrée/sortie) n = ln(P₂/P₁) / ln(v₁/v₂)

// Facteur correcteur de Schultz f (≈ 1, cale la formule sur la réalité) f = H_s / [ (n_s/(n_s−1)) · (P₂·v₂ₛ − P₁·v₁) ]

// Hauteur polytropique corrigée (forme de Schultz) H_p = f · (n/(n−1)) · (P₂·v₂ − P₁·v₁)

⚠️ Pourquoi c'est important : plus le gaz s'éloigne du gaz parfait (Z loin de 1, X loin de 0, Y loin de 1), plus l'écart entre le calcul gaz réel et le calcul gaz parfait grandit. Sur les hydrocarbures lourds ou à haute pression, négliger Schultz fausse le dimensionnement du compresseur et de son moteur d'entraînement.

💡 Exemple — compression d'un gaz naturel

🗜️ Cas — gaz naturel 10 → 40 bar

Mélange (molaire) : 90 % méthane, 5 % éthane, 3 % CO₂, 2 % azote. Aspiration à 20 °C / 10 bar abs, refoulement 40 bar abs, η_p = 0,78, débit 5 kg/s.

Propriétés à l'aspiration : masse molaire 17,8 g/mol, Z = 0,979, exposant isentropique k = 1,33, coefficients de Schultz X = 0,076 et Y = 1,02 (déjà loin des valeurs « gaz parfait » 0 et 1).

Résultats du calcul gaz réel :
— Température de refoulement : 155,6 °C (idéal isentropique : 125,8 °C)
— Z passe de 0,979 (entrée) à 0,988 (sortie) · exposant polytropique n = 1,39
— Facteur de Schultz f = 1,0018 · rendement isentropique η_s = 0,746
— Hauteur polytropique H_p = 227,8 kJ/kg · puissance 1 460 kW

Comparaison gaz parfait : la même hauteur calculée en supposant Z = 1 donnerait 231,7 kJ/kg, soit un écart de ≈ 1,7 % rien qu'à 40 bar. L'écart croît vite avec la pression : c'est tout l'intérêt de la correction de Schultz.

❄️ Compression multi-étages avec refroidissement

Pour les forts taux de compression, on comprime en plusieurs étages en refroidissant le gaz entre chaque étage (réfrigérant intermédiaire). Deux bénéfices : on limite la température de refoulement, et on réduit la puissance totale (comprimer du gaz froid coûte moins cher). La répartition optimale donne le même taux de compression à chaque étage :

// Taux de compression par étage (réparti également) r = (P₂/P₁)^(1/N) // N = nombre d'étages

// Entre deux étages : on refroidit le gaz (souvent jusqu'à la T d'aspiration) Chaleur évacuée = débit · (h_chaud − h_refroidi)

❄️ Cas — gaz naturel 1 → 40 bar (η_p = 0,78, 5 kg/s)

Le même mélange, comprimé d'1 à 40 bar (taux total 40), avec refroidissement à 20 °C entre les étages :

Étages	Taux / étage	Puissance totale	T refoulement	Gain vs 1 étage
1 (sans refroid.)	40,0	5 327 kW	422 °C	—
2	6,32	4 185 kW	205 °C	−21 %
3	3,42	3 828 kW	139 °C	−28 %

Lecture : en passant de 1 à 3 étages, la puissance chute de 28 % et la température de refoulement passe de 422 °C (inacceptable pour les joints et l'huile) à 139 °C. En contrepartie, il faut évacuer la chaleur dans les réfrigérants (≈ 2 600 kW pour la version 3 étages).

📏 Débitmètre — conversions & diaphragme ISO 5167

Un débit de gaz peut s'exprimer de trois façons, qu'il faut savoir convertir. Seul le débit massique (kg/s, kg/h) est conservé quelles que soient T et P ; les débits volumiques dépendent des conditions :

// Débit volumique RÉEL (aux conditions T, P de la conduite) q_v réel = débit massique / ρ(T, P)

// Débit volumique NORMAL : volume ramené à 0 °C et 1,01325 bar q_v normal (Nm³/h) = débit massique / ρ_normale

Le deuxième outil calcule le débit à partir d'un diaphragme (un disque percé qui crée une perte de charge Δp mesurable), selon la norme ISO 5167. C'est le principe des débitmètres déprimogènes :

// Débit massique à travers un diaphragme (ISO 5167) q_m = (C / √(1−β⁴)) · ε · (π/4)·d² · √(2·Δp·ρ₁)

// β = rapport des diamètres ; C = coeff. de décharge ; ε = expansibilité β = d / D // d = Ø orifice, D = Ø conduite C ≈ 0,6 // équation de Reader-Harris/Gallagher (dépend du Reynolds) ε ≤ 1 // correction de compressibilité du gaz

📏 Cas — diaphragme sur conduite de gaz naturel

Données : conduite D = 100 mm, orifice d = 50 mm (β = 0,5), gaz à 20 °C / 10 bar abs, perte de charge mesurée Δp = 200 mbar, prises dans l'angle.

Résultat : coefficient de décharge C = 0,604, expansibilité ε = 0,994, nombre de Reynolds ≈ 7,5×10⁵ (régime turbulent, domaine valide).
Débit obtenu : 2 396 kg/h, soit 3 005 Nm³/h ou 321 m³/h réels.

Vérification croisée : ces trois valeurs sont bien cohérentes entre elles via les masses volumiques réelle (7,47 kg/m³) et normale (0,797 kg/Nm³) — exactement ce que donne l'outil de conversion de débit.

🧭 Les 4 modes de l'onglet

Mode	On saisit	On obtient
Propriétés	Mélange + T + P	Z, masse volumique, Cp/Cv, k, vitesse du son, X et Y
Compression	P₁, P₂, débit + soit η_p, soit la T₂ mesurée	L'autre des deux (T₂ ou η_p), hauteurs, puissance, f, n
Multi-étages	P₁, P₂, η_p, débit, nb d'étages, T de refroid.	Détail par étage, puissance totale, chaleur évacuée, gain
Débitmètre	Une valeur de débit, ou un diaphragme (D, d, Δp)	Tous les débits (kg, Nm³, m³ réel), C, ε, Reynolds

Section 6

🔄 Flux de travail recommandé

Les cinq modules de l'application s'utilisent dans un ordre logique. Voici le flux de travail recommandé pour une analyse complète :

1

Planning turbines — Identifiez les turbines dont la révision approche. Les révisions L2 (25 000h) et L3 (50 000h) sont des points naturels pour réévaluer les stratégies de maintenance.

2

RCM II — Pour la turbine concernée, créez ou mettez à jour une analyse RCM. Listez les fonctions, les modes de défaillance, évaluez P et C pour chaque mode, et sélectionnez la tâche de maintenance appropriée.

3

Analyse Weibull — Pour les modes de défaillance avec historique de pannes, lancez une analyse Weibull. Si β > 1, vous avez un âge limite → TRP/TRE justifiée. Si β ≈ 1, TD ou RTF.

3b

Weibull Thermique — Si la température d'exploitation varie (surcharge turbine, surchauffe bobinage), utilisez ce module pour recalculer η à la température réelle. Comparez les courbes F(t) pour quantifier l'impact sur les intervalles de maintenance.

4

BCR & Isorisque — Pour les modes à risque élevé ou critique, entrez les actions candidates identifiées en RCM et calculez leur BCR. Visualisez le déplacement dans la matrice de risque.

5

Courbe €/RU — Lisez le seuil d'arbitrage pour prioriser vos actions selon le budget disponible. Retenir systématiquement les actions avant le coude de la courbe.

🎯 Règle de décision synthétique

Zone de risque	BCR > 2	BCR 1–2	BCR < 1
Critique (25–36)	✅ Agir immédiatement	✅ Agir + optimiser	⚠️ Repenser l'action
Élevé (15–24)	✅ Retenir	⚠️ Étudier le budget	❌ Rejeter ou RTF partiel
Modéré (7–14)	✅ Retenir si budget ok	❌ RTF + surveillance	❌ RTF
Faible (1–6)	✅ Si très peu cher	❌ RTF	❌ RTF

Section 7

📖 Lexique

Arrhenius (modèle)

Modèle physico-chimique décrivant l'accélération de la dégradation avec la température. La durée caractéristique η varie selon η(T) = C × exp(Ea / (k×T)). Utilisé pour les composants mécaniques et thermiques (turbines, injecteurs). Référence : MIL-HDBK-217, IEC 61709.

BCR

Benefit-to-Cost Ratio. Ratio bénéfice sur coût d'une action de maintenance. BCR > 1 = rentable.

Censure

En analyse Weibull, une observation censurée est une durée de fonctionnement sans panne (équipement encore en service ou retiré pour une autre raison). Elle doit être incluse dans l'analyse pour ne pas biaiser les résultats.

Courbe P-F

Courbe Potential Failure to Functional Failure. Représente la dégradation progressive d'un équipement. L'intervalle P-F est la fenêtre pendant laquelle la défaillance est détectable avant qu'elle ne devienne fonctionnelle.

€/RU

Euro par Risk Unit. Coût d'une action de maintenance par point de risque éliminé. Permet de comparer des actions qui réduisent des quantités différentes de risque.

Ea (énergie d'activation)

Paramètre clé du modèle d'Arrhenius, exprimé en électronvolts (eV). Représente la "résistance" d'un matériau à la dégradation thermique. Plus Ea est élevé, plus la durée de vie est sensible à la température. Valeurs typiques : 0,3–0,5 eV (mécanique) à 1,2–1,4 eV (semiconducteurs).

Facteur d'accélération (AF)

Ratio entre la durée de vie à la température d'utilisation et celle à une température de stress. AF = 3 signifie que l'équipement vieillit 3× plus vite à la température de stress. Formule : AF = exp(Ea/k × (1/T_use − 1/T_stress)).

FOD

Foreign Object Damage. Dommage causé par l'ingestion d'un corps étranger dans le compresseur. Cause principale de rupture d'aube en phase aléatoire (β ≈ 1).

Isorisque

Courbe reliant tous les couples (P, C) ayant le même score de risque. Équation : P × C = constante (hyperbole). Utilisée pour visualiser le déplacement dans la matrice de risque suite à une action de maintenance.

Montsinger (loi de)

Règle empirique (1930) normalisée IEC 60085 : chaque augmentation de 10 °C divise la durée de vie d'un isolant électrique par 2 (et inversement). Formalisée par L(T) = L₀ × exp(−α × ΔT), avec α = ln(2)/10 ≈ 0,069. Applicable aux bobinages moteurs, transformateurs et câbles.

L2 / L3

Paliers de révision des turbines LM2500. L2 = révision intermédiaire à 25 000h (durée ≈ 30 jours). L3 = révision majeure à 50 000h (durée ≈ 200 jours).

Mode de défaillance

Manière précise dont un équipement cesse d'assurer sa fonction. Un équipement peut avoir plusieurs fonctions, chacune avec plusieurs modes de défaillance potentiels.

RCM

Reliability Centered Maintenance. Méthode structurée de sélection des tâches de maintenance basée sur les conséquences des défaillances plutôt que sur l'âge des équipements. Développée pour l'aviation, aujourd'hui appliquée à l'industrie.

RTF

Run-To-Failure. Décision délibérée de ne pas faire de maintenance préventive et d'attendre la panne. Valide économiquement quand le coût de la maintenance dépasse le coût de la panne et que le risque est faible.

Weibull β (bêta)

Paramètre de forme de