Un bloc de PertDP représente un tuyau hydraulique avec des pertes de charge. Ce composant représente une perte de charge hydraulique où la perte de pression est directement proportionnelle au débit, état laminaire. Ce bloc a une direction; c.-à-d. la direction positive est quand le fluide entre par le port d'entrée (le noir). Ceci signifie que ce bloc suppose que le débit est positif si le fluide découle du port noir vers le port blanc. Les paramètres importants de ce composant sont la longueur, le diamètre du tuyau, les altitudes des ports d'entrée et de sortie, et quelques autres coefficients thermo-hydrauliques.
Parametres du tuyau | |
Longueur du tube : L (m) | 10 |
Diamètre interne du tube : D (m) | 0.2 |
Coefficient de perte de charge-frottement(S.U) : lambda | 0.03 |
Altitude entrée tuyauterie : z1 (m) | 0 |
Altitude sortie tuyauterie : z2 (m) | 0 |
Si ![]() |
0 |
class PerteDP parameter Real L "Longueur du tube"; parameter Real D "Diametre interne du tube"; parameter Real lambda "Coefficient de perte de charge-frottement (s.u.)"; parameter Real z1 "Altitude entree tuyauterie"; parameter Real z2 "Altitude sortie tuyauterie"; parameter Real p_rho "Si > 0, masse volumique imposée du fluide"; // parameter Integer mode "Région IF97. 1:liquide - 2:vapeur - 4:saturation - 0:calcul automatique"; parameter Real g = 9.80665 "Accélération de la pesanteur"; parameter Real pi = 3.1415927; parameter Real eps = 1.e-0 "Limite inf. de la fonction carrée"; Real khi "Coefficient de perte de charge hydraulique"; Real deltaPf "Perte de charge par frottement"; Real deltaP "Perte de charge totale"; Real Q (start=500) "Débit massique"; Real rho (start=998) "Masse volumique"; Real Tm (start=290) "Température moyenne"; Real Pm (start=1.e5) "Pression moyenne"; Real Hm (start=100000) "Enthalpie spécifique moyenne"; Real region (start=1) "Numéro de région IF97"; Real ThermoCar; public PortPHQ1 C1 ; PortPHQ2 C2 ; equation C1.P - C2.P = deltaP; C2.Q = C1.Q; C2.H = C1.H; Q = C1.Q; /* Inversions de débit */ 0 = if (Q > 0.0) then (C1.H - C1.Hm) else (C2.H - C2.Hm); /* Calcul des pertes de charges */ deltaP = deltaPf + rho*g*(z2 - z1); ThermoCar = if (abs(Q) > eps) then Q*abs(Q) else Q*eps; deltaPf = 8*khi*ThermoCar/(pi^2*D^4*rho); /* Formule de Darcy-Weisbach (Idel'cik p. 55). On suppose qu'on est en régime quadratique. */ khi = lambda*L/D; /* Calcul des propriétés thermodynamiques */ Pm = (C1.P + C2.P)/2; Hm = (C1.H + C2.H)/2; Tm = Hm/4187 + 273.15; region = 1; rho = 998; end PerteDP;