Calcul du Temps de Réponse d’un Système Asservi
📝 Situation Industrielle
Vous êtes ingénieur au sein du Service Technique Central d'un site sidérurgique majeur (type ArcelorMittal ou Tata Steel), plus précisément affecté au laminoir à froid. Cette installation colossale, longue de plusieurs centaines de mètres, réduit l'épaisseur de bobines d'acier (coils) par écrasement successif entre des cylindres de travail.
La qualité du produit fini (tôle pour carrosserie automobile) exige une régularité d'épaisseur absolue : la tolérance est de +/- 2 microns à une vitesse de défilement de 1200 m/min. Pour atteindre cette performance, les cages de laminage sont équipées d'un système de serrage hydraulique ultra-rapide appelé HAGC (Hydraulic Automatic Gauge Control). Ces vérins, capables de développer 2000 tonnes de poussée, doivent corriger en temps réel les moindres variations de dureté de la bande entrante.
Le Problème : Suite au remplacement des cylindres d'appui par des modèles plus rigides (et donc plus lourds), les opérateurs ont constaté des "pompages" (oscillations) lors des phases d'accélération. Le système semble avoir perdu en stabilité.
Votre Direction Technique vous mandate pour vérifier la fréquence propre hydraulique de la nouvelle configuration. Vous devez déterminer si l'augmentation de la masse mobile a abaissé la fréquence de résonance à un niveau critique, compromettant le temps de réponse du système (cible < 30 ms).
"Les pressions de service atteignent 280 bars. Ne jamais intervenir sur les accumulateurs ou les flexibles sans avoir consigné l'installation et vérifié la décharge complète de la pression résiduelle."
L'étude dynamique repose sur une caractérisation précise des trois sous-systèmes : le fluide de transmission, l'actionneur (géométrie) et la charge mécanique (inertie). Chaque donnée a été vérifiée par le service Métrologie.
📚 Référentiels Applicables
Les calculs doivent se conformer aux standards internationaux pour garantir la validité des résultats auprès du constructeur de la servo-valve.
ISO 10770-1 : Méthode d'essai des régulateurs de débit Théorie des Systèmes Asservis LinéairesLe comportement "ressort" de l'huile est le facteur prédominant de la résonance. Le module de compressibilité (\(\beta\)) n'est pas une constante fixe : il chute drastiquement en présence d'air dissous ou de flexibles souples. Nous utilisons ici un module "effectif" (\(\beta_{\text{eff}}\)) mesuré sur banc, bien inférieur au module théorique (1.8 GPa).
| Fluide utilisé | Huile Minérale ISO VG 46 |
| Module de Compressibilité Effectif (\(\beta_{\text{eff}}\)) | 14 000 bar (soit 1.4 GPa) |
| Pression d'alimentation (\(P_s\)) | 280 bar |
Le vérin HAGC est un vérin "plongeur" ou double effet à simple tige massive. La raideur hydraulique est inversement proportionnelle au volume d'huile emprisonné. Le "Volume Mort Total" inclut le volume dans les chambres du vérin PLUS le volume dans les tuyauteries entre la servo-valve et le vérin. C'est ce volume total qui se comprime.
| Type | Vérin Double Effet Symétrique |
| Diamètre Piston (\(D_p\)) | 320 mm |
| Course Totale | 50 mm (Position de travail à mi-course) |
| Volume Mort Total (Chambres + Tuyaux) (\(V_t\)) | 5.8 Litres |
C'est la modification majeure du projet. Les nouveaux cylindres d'appui (Back-up Rolls) sont plus denses. La "Masse Mobile" représente tout ce qui bouge lors de la régulation : la tige du vérin, les empoises (paliers), et les cylindres eux-mêmes. L'amortissement visqueux (\(\xi\)) est estimé d'après les frottements des joints et des glissières.
| Masse Mobile Totale (\(M\)) | 12 500 kg (Cylindres + Empoises) |
| Coefficient d'amortissement estimé (\(\xi\)) | 0.55 (Système sous-amorti) |
| Exigence Client (Temps de réponse \(t_{\text{r}, 5\%}\)) | < 30 ms |
E. Protocole de Résolution
Pour valider la conformité dynamique du système, nous allons procéder par une approche analytique en "boucle ouverte" focalisée sur la fréquence propre hydraulique.
Calcul de la Surface Active
Détermination de la section du vérin soumise à la pression, base de tous les calculs de force et de raideur.
Détermination de la Raideur Hydraulique
Calcul de l'effet "ressort" de la colonne d'huile (\(K_h\)) en fonction du volume emprisonné et du module de compressibilité.
Pulsation Propre Hydraulique
Identification de la fréquence de résonance naturelle (\(\omega_0\)) du système Masse-Ressort.
Estimation du Temps de Réponse
Déduction du temps de réponse indiciel à 5% (\(t_{\text{r}, 5\%}\)) à partir des abaques du second ordre.
Calcul du Temps de Réponse d’un Système Asservi
🎯 Objectif Scientifique
Dans un système hydraulique de puissance, la surface du piston est l'interface critique de conversion d'énergie. Elle transforme l'énergie potentielle du fluide (Pression) en énergie mécanique (Force). Mais au-delà de la force, cette surface joue un rôle capital dans la rigidité du système : plus la surface est grande, plus la colonne d'huile est large, ce qui influence directement sa compressibilité. Notre premier objectif est donc de déterminer avec précision la section active en unités SI (mètres carrés) pour garantir la cohérence dimensionnelle des calculs de raideur à venir.
📚 Référentiel Mathématique
Géométrie Euclidienne (Aire du Disque) Système International (SI)Le diamètre du vérin est fourni en millimètres (\(320 \text{ mm}\)), ce qui est la norme en dessin industriel. Cependant, les formules de dynamique des fluides (notamment celles impliquant des Pascals et des Joules) exigent impérativement l'usage du mètre. Une erreur classique consiste à calculer en mm² puis à convertir, ce qui induit souvent des erreurs d'ordre de grandeur \(10^6\). Pour sécuriser le calcul, la stratégie la plus robuste est de convertir le diamètre en mètres dès le départ, avant même d'élever au carré.
Le principe de Pascal établit que la pression se transmet intégralement dans un fluide incompressible. La force générée par un vérin est le produit de cette pression par la surface sur laquelle elle s'applique : \(F = P \times S\). Pour un vérin cylindrique, cette surface est un disque parfait. La précision de ce calcul conditionne directement l'estimation de la force de laminage et, comme nous le verrons, de la raideur \(K_h\).
Étape 1 : Données d'Entrée
| Paramètre | Symbole | Valeur Brute |
|---|---|---|
| Diamètre Piston | \(D_p\) | 320 mm |
Pour passer de mm à m, déplacez la virgule de 3 rangs vers la gauche. \(320 \text{ mm} = 0.320 \text{ m}\). Évitez absolument d'utiliser \(S = \pi \cdot R^2\) si on vous donne le diamètre, car diviser le diamètre par 2 introduit une étape de calcul supplémentaire source d'erreur et d'arrondi. La formule avec \(D^2/4\) est plus directe.
Étape 2 : Application Numérique Détaillée
1. Conversion PréalableOn convertit le diamètre en unité SI de base.
On élève le diamètre au carré avant de multiplier par Pi.
On applique la formule complète.
Ce résultat de 0.08 m² (soit environ 800 cm²) est typique pour un vérin de laminoir. C'est une surface considérable, équivalente à une feuille A4 environ, capable de générer des forces colossales sous 280 bars.
La surface active de 0.0804 m² est la pierre angulaire de notre modélisation. Elle va maintenant être élevée au carré dans le calcul de raideur, ce qui signifie que toute imprécision ici sera amplifiée. La valeur obtenue confirme que nous sommes en présence d'un actionneur de grande puissance.
Vérifions l'ordre de grandeur. Un carré de 30cm de côté a une aire de \(0.3 \times 0.3 = 0.09 \text{ m}^2\). Notre disque de 32cm de diamètre doit être légèrement plus petit que le carré qui l'entoure (\(32 \times 32 \approx 0.1 \text{ m}^2\)). Le résultat de 0.08 m² est donc parfaitement cohérent.
Ne jamais confondre Rayon et Diamètre dans la formule. Si vous utilisez \(\pi R^2\), n'oubliez pas de diviser D par 2 d'abord. Ici, l'erreur classique est d'oublier le carré ou la division par 4.
🎯 Objectif Scientifique
C'est le cœur du problème dynamique. Nous devons quantifier l'effet "ressort" de la colonne d'huile. Bien que les liquides soient souvent considérés comme incompressibles en statique, ils se comportent comme des ressorts très raides sous haute pression et haute dynamique. Nous cherchons à calculer \(K_h\) (en Newtons par mètre), qui représente la résistance du fluide à l'écrasement. Plus cette valeur est élevée, plus le système sera réactif et précis.
📚 Référentiel Physique
Loi de Hooke (Généralisée) Module de Compressibilité (Bulk Modulus)La raideur d'un vérin hydraulique varie selon la position du piston. Le pire cas (la raideur la plus faible, donc la plus critique pour la stabilité) se produit lorsque le volume d'huile comprimé est maximal. Pour un vérin double effet symétrique (ou considéré comme tel en asservissement autour du point milieu), la raideur dépend du module de compressibilité \(\beta_{\text{eff}}\), de la section \(S\) et du volume total \(V_t\).
Attention : Nous devons utiliser le \(\beta_{\text{eff}}\) (effectif) et non théorique. La présence de flexibles et de micro-bulles d'air divise souvent la raideur réelle par 2 ou 3 par rapport à la théorie. Ici, la valeur de 14 000 bar est une donnée expérimentale fiable qu'il faut respecter.
La raideur \(K\) se définit comme \(K = \frac{Force}{Deplacement}\). Pour une colonne de fluide, la variation de volume \(\Delta V\) sous une variation de pression \(\Delta P\) est régie par le module \(\beta_{\text{eff}}\) : \(\Delta P = -\beta_{\text{eff}} \frac{\Delta V}{V}\).
En traduisant cela en termes de force (\(F = P \cdot S\)) et de déplacement (\(\Delta V = S \cdot \Delta x\)), on obtient la raideur d'une colonne simple : \(K = \frac{\beta_{\text{eff}} S^2}{V}\).
Pour un vérin asservi bloqué par sa valve (deux chambres en opposition de pression), les deux "ressorts" d'huile (chambre avant et arrière) agissent en parallèle hydrauliquement par rapport au mouvement de la charge. Pour un piston centré, cela introduit un facteur 4 dans la formule.
La raideur globale \(K_h\) est donnée par :
Où \(V_t\) est le volume total d'huile emprisonné (chambres + tuyauteries) en [m³].
Étape 1 : Données d'Entrée
| Paramètre | Valeur | Unité SI requise |
|---|---|---|
| Surface Active (\(S\)) | 0.08042 | m² |
| Module Compressibilité (\(\beta_{\text{eff}}\)) | 14 000 bar | Pa (Pascal) |
| Volume Mort Total (\(V_t\)) | 5.8 Litres | m³ |
Ne vous lancez jamais dans le calcul avec des bars et des litres !
1 bar = \(10^5\) Pa. Donc \(14\,000 \text{ bar} = 1.4 \times 10^9 \text{ Pa}\) (ou 1.4 GPa).
1 Litre = \(1 \text{ dm}^3 = 10^{-3} \text{ m}^3\). Donc \(5.8 \text{ L} = 0.0058 \text{ m}^3\).
Étape 2 : Calculs Détaillés
1. Préparation des valeursOn convertit d'abord toutes les grandeurs.
La surface joue un rôle quadratique, on la calcule séparément.
On calcule le terme de "force" équivalente au numérateur.
Le résultat est colossal : \(6.24 \times 10^9 \text{ N/m}\). Cela signifie que le système est extrêmement rigide. Pour comprimer ce vérin de seulement 1mm, il faudrait appliquer une force de 6.24 millions de Newtons (624 tonnes) ! C'est cette rigidité colossale qui permet de laminer l'acier avec précision.
Nous avons déterminé la constante de raideur de notre "ressort liquide". Avec \(6.24 \text{ GN/m}\), nous disposons d'un actionneur très raide, capable de transmettre les efforts sans s'écraser. Cette valeur élevée est prometteuse pour obtenir une fréquence propre élevée.
Les raideurs hydrauliques sont toujours très élevées. Un ordre de grandeur en \(10^9 \text{ N/m}\) (Giga-Newton/mètre) est standard pour de la grosse hydraulique industrielle. Un résultat en \(10^6\) ou \(10^{12}\) aurait signalé une erreur d'unité (mm vs m).
L'erreur fatale ici est d'oublier de convertir les litres en m³. Si vous divisez par 5.8 au lieu de 0.0058, votre raideur sera 1000 fois trop faible, et tout le dimensionnement s'effondre.
🎯 Objectif Scientifique
Nous connaissons la raideur du "ressort" d'huile (\(K_h\)) et la masse qu'il doit déplacer (\(M\)). Tout système Masse-Ressort possède une fréquence de résonance naturelle. Si l'on essaie de piloter le système plus vite que cette fréquence, il ne suivra plus et se mettra à osciller ou à atténuer le signal. Calculer cette fréquence, c'est trouver la "limite de vitesse" physique de notre actionneur.
📚 Référentiel Physique
Dynamique des Structures Oscillateur Harmonique SimpleLa formule est universelle en mécanique : \(\omega = \sqrt{K/M}\). Plus la raideur est élevée, plus le système est réactif (fréquence haute). Plus la masse est lourde, plus le système est lent (fréquence basse). Dans notre laminoir, nous avons une masse énorme (12.5 tonnes) mais une raideur gigantesque. Lequel va l'emporter ? Le calcul nous le dira.
Tout système masse-ressort non amorti oscille naturellement à une pulsation \(\omega_0\) (en radians/seconde). Cette pulsation ne dépend que des caractéristiques physiques intrinsèques (K et M), et non de la commande.
La relation fondamentale est \(\omega_0^2 = \frac{K}{M}\).
La fréquence \(f_0\) (en Hertz, ou cycles par seconde) est simplement une autre manière d'exprimer cette vitesse angulaire : \(f_0 = \frac{\omega_0}{2\pi}\).
La pulsation naturelle \(\omega_0\) (en radians par seconde) est :
Avec \(K_h\) en [N/m] et \(M\) en [kg]. Le résultat \(\omega_0\) est en [rad/s].
Étape 1 : Données d'Entrée
| Paramètre | Valeur |
|---|---|
| Raideur Calculée (\(K_h\)) | \(6.244 \times 10^9 \text{ N/m}\) |
| Masse Mobile (\(M\)) | \(12\,500 \text{ kg}\) |
Vérifiez toujours vos unités sous la racine carrée : \(\frac{\text{N/m}}{\text{kg}} = \frac{(\text{kg} \cdot \text{m} / \text{s}^2) / \text{m}}{\text{kg}} = \frac{1}{\text{s}^2}\). La racine carrée donne bien des \(\text{s}^{-1}\) (radians/seconde). C'est homogène.
Étape 2 : Calculs Détaillés
1. Calcul du rapport K/MLa fréquence \(f_0\) est plus parlante pour les ingénieurs de contrôle.
Nous obtenons une fréquence propre naturelle de 112.5 Hz.
Dans le monde de l'hydraulique industrielle lourde, c'est une valeur excellente. La plupart des vérins standards oscillent entre 10 et 40 Hz. Le fait d'être au-dessus de 100 Hz indique que la rigidité hydraulique domine très largement l'inertie de la masse, malgré l'augmentation de poids des cylindres.
Avec une fréquence propre de 112.5 Hz, notre système mécanique a un potentiel dynamique très élevé. Cela signifie qu'il est capable, en théorie, de suivre des signaux variant jusqu'à cette fréquence avant d'entrer en résonance. C'est une base très saine pour un asservissement performant.
Une fréquence de 112 Hz correspond à une période de \(T = 1/112 \approx 9 \text{ ms}\). Cela suggère que le système est capable de réagir très vite (dans la gamme de la dizaine de millisecondes), ce qui est bon signe pour l'objectif final.
Attention, cette fréquence est la fréquence propre hydraulique. La fréquence propre structurelle du bâti de la cage (élasticité des montants) est souvent plus basse (autour de 60-80 Hz). Le maillon faible réel pourrait être la structure, et non l'hydraulique.
🎯 Objectif Scientifique
C'est la conclusion et la validation finale de notre mission. Le client a une exigence claire : le système doit répondre en moins de 30 ms (critère \(t_{\text{r}, 5\%}\)). Ce temps correspond à la durée nécessaire pour que, suite à un ordre brusque (échelon), la position réelle du vérin atteigne la consigne et y reste avec une erreur inférieure à 5%. Nous devons lier les paramètres abstraits calculés précédemment (\(\omega_0, \xi\)) à cette grandeur temporelle concrète.
📚 Référentiel Physique
Automatique (Systèmes du 2nd Ordre) Abaques de Réponse IndicielleUn axe hydraulique se comporte typiquement comme un système du second ordre (masse-ressort-amortisseur). Son temps de réponse dépend de deux facteurs :
1. La "vitesse" intrinsèque, donnée par \(\omega_0\).
2. L'amortissement \(\xi\).
Avec \(\xi = 0.55\), le système est "sous-amorti". Il est vif mais va dépasser la consigne (overshoot) avant de se stabiliser. Ce dépassement "perd du temps" car on doit attendre que les oscillations rentrent dans la bande des 5%. Pour estimer ce temps sans résoudre l'équation différentielle complexe, on utilise des abaques ou des formules d'approximation empiriques valables pour \(0.4 < \xi < 0.8\).
En automatique, on définit le "temps de réponse réduit" \(TR_{\text{red}} = t_r \cdot \omega_0\). C'est une valeur sans dimension qui ne dépend que de l'amortissement \(\xi\).
Pour \(\xi \approx 0.7\) (amortissement optimal), \(TR_{\text{red}} \approx 3\).
Pour \(\xi = 0.55\), le dépassement est plus fort, ce qui peut légèrement augmenter cette valeur si le premier pic sort des 5%, mais l'approximation \(TR_{\text{red}} \approx 3\) reste une excellente base pour un pré-dimensionnement rapide, souvent utilisée côté sécurité (car un \(\xi\) plus faible rend la montée plus rapide au détriment de la stabilité).
Pour estimer le temps de réponse à 5% d'un système du 2nd ordre :
Cette formule est une enveloppe couramment admise en ingénierie hydraulique.
Étape 1 : Données d'Entrée
| Paramètre | Valeur |
|---|---|
| Pulsation Propre (\(\omega_0\)) | \(706.76 \text{ rad/s}\) |
| Amortissement (\(\xi\)) | 0.55 |
Si vous n'avez pas de calculatrice scientifique, retenez que pour un bon système hydraulique (\(\xi \approx 0.7\)), le temps de réponse est à peu près égal à la moitié de la période propre : \(t_r \approx T_0 / 2 = 1 / (2 \cdot f_0)\).
Étape 2 : Calculs Détaillés
1. Calcul du DénominateurOn multiplie l'amortissement par la pulsation.
On applique la formule approximative.
Le calcul théorique donne 7.7 ms. C'est extrêmement bas par rapport à la cible de 30 ms. Cela montre que l'axe hydraulique "nu" (vérin + huile + masse) est très rapide.
Le temps de 7.7 ms ne prend en compte que la dynamique mécanique. Dans la réalité, il faut ajouter le temps de réponse de la servo-valve elle-même (qui met environ 10 ms à s'ouvrir). En règle générale, on additionne quadratiquement les temps de réponse : \(T_{\text{total}} = \sqrt{T_{\text{hydro}}^2 + T_{\text{valve}}^2} \approx \sqrt{7.7^2 + 10^2} \approx 12.6 \text{ ms}\). Même avec cette correction, nous restons très confortablement sous la barre des 30 ms.
Le temps de réponse calculé (7.7 ms) est légèrement inférieur à la période de l'oscillation naturelle (9 ms). C'est typique d'un système sous-amorti qui monte très vite ("temps de montée" rapide) avant d'osciller.
Attention, ce calcul suppose que la pression d'alimentation (280 bars) est suffisante pour accélérer la masse à cette vitesse (saturation en débit/effort). Si la servo-valve est trop petite (saturation en débit), le temps réel sera beaucoup plus long. Ce calcul valide la dynamique fréquentielle, pas le dimensionnement en puissance.
Participer à la discussion