Xcos Modelica libraries: encode Hydraulics.mo in UTF-8 39/18739/2
Paul Bignier [Fri, 2 Dec 2016 15:46:05 +0000 (16:46 +0100)]
 * Like Electrical.mo. OpenModelica needs UTF-8 library files

Change-Id: Ida90caa60762b5c696aab9241247ad33ead5f3bd

scilab/modules/scicos_blocks/macros/Hydraulics/Hydraulics.mo

index 4ad11fc..564c3e9 100644 (file)
@@ -1,39 +1,39 @@
-connector PortPHQ1 
+connector PortPHQ1
 //  extends ThermoLib.Connecteurs.PortPHQ;
   Real P "Pression du fluide";
-  Real H "Enthalpie spécifique du fluide";
-  Real Q "Débit massique du fluide";
-  Real Hm "Enthalpie spécifique de mélange";
+  Real H "Enthalpie spécifique du fluide";
+  Real Q "Débit massique du fluide";
+  Real Hm "Enthalpie spécifique de mélange";
 
-//  input Boolean a=true "Pseudo-variable pour la vérification du sens des connexions";
-//  output Boolean b "Pseudo-variable pour la vérification du sens des connexions";
+//  input Boolean a=true "Pseudo-variable pour la vérification du sens des connexions";
+//  output Boolean b "Pseudo-variable pour la vérification du sens des connexions";
 end PortPHQ1;
-connector PortPHQ2 
+connector PortPHQ2
 //  extends ThermoLib.Connecteurs.PortPHQ;
   Real P "Pression du fluide";
-  Real H "Enthalpie spécifique du fluide";
-  Real Q "Débit massique du fluide";
-  Real Hm "Enthalpie spécifique de mélange";
+  Real H "Enthalpie spécifique du fluide";
+  Real Q "Débit massique du fluide";
+  Real Hm "Enthalpie spécifique de mélange";
 
-//  input Boolean a "Pseudo-variable pour la vérification du sens des connexions";
-//  output Boolean b=true "Pseudo-variable pour la vérification du sens des connexions";
+//  input Boolean a "Pseudo-variable pour la vérification du sens des connexions";
+//  output Boolean b=true "Pseudo-variable pour la vérification du sens des connexions";
 
 end PortPHQ2;
 model Source
   parameter Real P0=300000 "Pression de la source";
-  parameter Real T0=290 "Température de la source";
-  parameter Real H0=100000 "Enthalpie spécifique de la source";
-  parameter Real option_temperature=1 "1:température fixée - 2:enthalpie fixée";
-   
+  parameter Real T0=290 "Température de la source";
+  parameter Real H0=100000 "Enthalpie spécifique de la source";
+  parameter Real option_temperature=1 "1:température fixée - 2:enthalpie fixée";
+
   Real Pm "Pression moyenne";
-  Real Q "Débit";
-  Real Tm "Température moyenne";
-  Real Hm "Enthalpie spécifique moyenne";
+  Real Q "Débit";
+  Real Tm "Température moyenne";
+  Real Hm "Enthalpie spécifique moyenne";
 
-public 
+public
   PortPHQ2 C ;
-equation 
-  
+equation
+
   C.P = Pm;
   C.Q = Q;
   C.Hm = Hm;
@@ -41,25 +41,25 @@ equation
   Pm = P0;
 
   Hm = H0;
-  
+
   Tm = if (option_temperature <= 1) then T0 else Hm/4187 + 273.15;
-  
+
 end Source;
 model Puits
   parameter Real P0=100000 "Pression du puits";
-  parameter Real T0=290"Température du puits";
-  parameter Real H0=100000 "Enthalpie spécifique du puits";
-  parameter Real option_temperature=1 "1:température fixée - 2:enthalpie fixée";
+  parameter Real T0=290"Température du puits";
+  parameter Real H0=100000 "Enthalpie spécifique du puits";
+  parameter Real option_temperature=1 "1:température fixée - 2:enthalpie fixée";
 
   Real Pm "Pression moyenne";
-  Real Q "Débit";
-  Real Tm "Température moyenne";
-  Real Hm "Enthalpie spécifique moyenne";
+  Real Q "Débit";
+  Real Tm "Température moyenne";
+  Real Hm "Enthalpie spécifique moyenne";
 
-public 
+public
   PortPHQ1 C ;
-equation 
-  
+equation
+
   C.P = Pm;
   C.Q = Q;
   C.Hm = Hm;
@@ -67,34 +67,34 @@ equation
   Pm = P0;
   Hm= H0;
 
-  
+
   Tm = if (option_temperature <= 1) then T0 else Hm/4187 + 273.15;
 
 end Puits;
-model PerteDP 
+model 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 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";
+  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 Q (start=500) "Débit massique";
   Real rho (start=998) "Masse volumique";
-  Real Tm (start=290) "Température moyenne";
+  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 Hm (start=100000) "Enthalpie spécifique moyenne";
+  Real region (start=1) "Numéro de région IF97";
   Real ThermoCar;
 
-public 
+public
   PortPHQ1 C1 ;
   PortPHQ2 C2 ;
 
@@ -103,62 +103,62 @@ equation
   C1.P - C2.P = deltaP;
   C2.Q = C1.Q;
   C2.H = C1.H;
-  
+
   Q = C1.Q;
 
-  /* Inversions de débit */
+  /* 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.0*D^4.0*rho);
-  
-  
-    /* Formule de Darcy-Weisbach (Idel'cik p. 55). On suppose qu'on est en régime quadratique. */
+
+
+    /* 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 */
+
+  /* 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;
-model Bache 
+model Bache
   parameter Real Patm "Pression dans le ciel de la bache";
   parameter Real A "Section de la bache";
-  parameter Real ze1 "Altitude du piquage d'entrée 1";
-  parameter Real ze2 "Altitude du piquage d'entrée 2";
+  parameter Real ze1 "Altitude du piquage d'entrée 1";
+  parameter Real ze2 "Altitude du piquage d'entrée 2";
   parameter Real zs1 "Altitude du piquage de sortie 1";
   parameter Real zs2 "Altitude du piquage de sortie 2";
   parameter Real z0 "Altitude initiale du fluide";
   parameter Real T0 "Temperature initiale du fluide";
-  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 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 e=0.25 "Ouverture des piquages";
   parameter Real k=1   "Coefficient de perte de charge des piquages";
-  parameter Real eps=1.e-0 "Limite inf. de la fonction carrée";
-  parameter Real zmin=1.e-6     "Hauteur minimum d'eau pour l'équation  de l'énergie";
+  parameter Real eps=1.e-0 "Limite inf. de la fonction carrée";
+  parameter Real zmin=1.e-6     "Hauteur minimum d'eau pour l'équation  de l'énergie";
 
   Real z(start=z0) "Hauteur du fluide";
-  Real Tm(start=T0) "Température moyenne";
+  Real Tm(start=T0) "Température moyenne";
   Real Pm(start=1.e5) "Pression moyenne";
-  Real Hm(start=100000) "Enthalpie spécifique moyenne";
+  Real Hm(start=100000) "Enthalpie spécifique moyenne";
   Real rho(start=998) "Masse volumique";
   Real BQ "Bilan de masse";
-  Real BH "Bilan d'énergie";
+  Real BH "Bilan d'énergie";
   Real Oe1;
   Real Oe2;
   Real Os1;
   Real Os2;
-  Real region(start=1) "Numéro de région IF97";
+  Real region(start=1) "Numéro de région IF97";
   Real Max_e1;
   Real Max_e2;
   Real Max_s1;
@@ -168,18 +168,18 @@ model Bache
 
   Real yNiveau ;
 
-public 
+public
 //  SortieReelle yNiveau ;
   PortPHQ1 Ce1 ;
   PortPHQ2 Cs2 ;
-//  PortTher Cth ; /* Port non utilisé dans le modèle désiré, non implanté sur le bloc SCICOS */
+//  PortTher Cth ; /* Port non utilisé dans le modèle désiré, non implanté sur le bloc SCICOS */
   PortPHQ1 Ce2 ;
   PortPHQ2 Cs1 ;
 
-//initial equation 
+//initial equation
 
 //  Hm = 100000;
-  
+
 equation
 
   Ce2.Q = 0;
@@ -205,7 +205,7 @@ equation
 
   Max_s1 = if (z - zs1 < 0.0) then 0.0 else (z - zs1);
   Cs1.P = Patm + rho*g*Max_s1;
-  Os1 = 0; 
+  Os1 = 0;
 
   Max_s2 = if (z - zs2 < 0.0) then 0.0 else (z - zs2);
   ThermoCar2 =if (abs(Cs2.Q) > eps) then Cs2.Q*abs(Cs2.Q) else Cs2.Q*eps;
@@ -220,9 +220,9 @@ equation
   /* Bilan de masse */
   BQ = Ce1.Q + Ce2.Q - Cs1.Q - Cs2.Q;
   rho*A*der(z) = BQ;
-  
-  /* Bilan d'énergie */
-  /* le port thermique n'est pas utilisé ici  "+ Cth.W" (on a W=0)*/
+
+  /* Bilan d'énergie */
+  /* le port thermique n'est pas utilisé ici  "+ Cth.W" (on a W=0)*/
   BH = Ce1.Q*(Ce1.H - Hm) + Ce2.Q*(Ce2.H - Hm) - Cs1.Q*(Cs1.H - Hm) - Cs2.Q*(Cs2.H - Hm);
 
   der(Hm) = if (z > zmin) then BH/(rho*A*z) else 0.0;
@@ -231,84 +231,84 @@ equation
   Ce2.Hm = Hm;
   Cs1.Hm = Hm;
   Cs2.Hm = Hm;
-  
-//  Cth.T = Tm; /* le port thermique n'est pas utilisé ici*/
-  
+
+//  Cth.T = Tm; /* le port thermique n'est pas utilisé ici*/
+
   /* Capteur de niveau */
   yNiveau = z;
-  
-  /* Calcul des propriétés thermodynamiques */
+
+  /* Calcul des propriétés thermodynamiques */
   Pm = Patm + rho*g*z/2;
-  
+
   Tm = Hm/4187 + 273.15;
   region = 1.0;
-  
+
   rho = 998.0;
 
 end Bache;
-model VanneReglante 
+model VanneReglante
   parameter Real Cvmax "CV max";
-  parameter Real p_rho "Si > 0, masse volumique imposée du fluide";
-//  parameter Integer mode=0 "Région IF97. 1:liquide - 2:vapeur - 4:saturation - 0:calcul automatique";
+  parameter Real p_rho "Si > 0, masse volumique imposée du fluide";
+//  parameter Integer mode=0 "Région IF97. 1:liquide - 2:vapeur - 4:saturation - 0:calcul automatique";
 
-  parameter Real eps=1.e-0 "Limite inf. de la fonction carrée";
+  parameter Real eps=1.e-0 "Limite inf. de la fonction carrée";
   Real Cv "Cv";
-  Real Q(start=500) "Débit massique";
-  Real deltaP "Perte de charge singulière";
+  Real Q(start=500) "Débit massique";
+  Real deltaP "Perte de charge singulière";
   Real rho(start=998) "Masse volumique";
-  Real Tm(start=290) "Température moyenne";
+  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 Hm(start=100000) "Enthalpie spécifique moyenne";
+  Real region(start=1) "Numéro de région IF97";
   Real ThermoCar;
 
   Real Ouv ;
 
-public 
+public
 //  EntreeReelle Ouv ;
   PortPHQ1 C1 ;
   PortPHQ2 C2 ;
 
-equation 
-  
+equation
+
   C1.H = C2.H;
   C1.Q = C2.Q;
-  
+
   Q = C1.Q;
   deltaP = C1.P - C2.P;
-  
-  /* Inversions de débit */
+
+  /* Inversions de débit */
   0 = if (Q > 0.0) then C1.H - C1.Hm else C2.H - C2.Hm;
-  
+
   /* Calcul de la perte de charge */
 
   ThermoCar = if (abs(Q) > eps) then Q*abs(Q) else Q*eps;
   deltaP*Cv*abs(Cv) = 1.733e12*ThermoCar/rho^2;
-  
+
   /* Calcul du Cv en fonction de l'ouverture */
   Cv = Ouv*Cvmax;
-  
-  /* Calcul des propriétés thermodynamiques */
+
+  /* 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 VanneReglante;
 
-model Flowmeter "Capteur de débit"
+model Flowmeter "Capteur de débit"
   parameter Real Qini=1 ;
 
-  Real Q (start=10) "Débit massique";
+  Real Q (start=10) "Débit massique";
   Real Mesure;
 
   PortPHQ1 C1;
   PortPHQ2 C2;
 
-equation 
+equation
   Mesure = Q;
 
   C1.P = C2.P;