Changeset 190 for branches/newlanguage/sample/controllers/CSTR_noniso_pid.mso

Timestamp:

Mar 7, 2007, 3:09:13 PM (17 years ago)

Author:

Rafael de Pelegrini Soares

Message:

Adapted PID models for the new language (the usage of switcher is still pending)

File:

• branches/newlanguage/sample/controllers/CSTR_noniso_pid.mso (modified) (7 diffs)

branches/newlanguage/sample/controllers/CSTR_noniso_pid.mso

@@ -22 +22 @@
 using "controllers/PIDs";
 
-const_valv as positive(Brief = "Valve Constant", Default=1,Lower=0,Upper=100, Unit="m^2.5/h");
+const_valv as positive(Brief = "Valve Constant", Default=1,Lower=0,Upper=100, Unit='m^2.5/h');
 
 Model corrente
@@ -34 +34 @@
 
         PARAMETERS
-        ko  as frequency                (Unit="1/h");
+        ko  as frequency                (DisplayUnit='1/h');
         A   as area;
         At  as area;    
-        Ea  as energy_mol               (Unit="kJ/kmol");
-        R   as Real                     (Unit="kJ/mol/K");
-        ro  as dens_mol                 (Unit="kmol/m^3");
-        Cp  as cp_mol                   (Unit="kJ/kmol/K");
-        U   as heat_trans_coeff (Unit="kW/m^2/K");
-        Hr  as heat_reaction    (Unit="kJ/kmol");
+        Ea  as energy_mol               (DisplayUnit='kJ/kmol');
+        R   as Real                     (Unit='kJ/mol/K');
+        ro  as dens_mol                 (DisplayUnit='kmol/m^3');
+        Cp  as cp_mol                   (DisplayUnit='kJ/kmol/K');
+        U   as heat_trans_coeff (DisplayUnit='kW/m^2/K');
+        Hr  as heat_reaction    (DisplayUnit='kJ/kmol');
         
         VARIABLES
@@ -53 +53 @@
         tau      as time_h;
         rA           as reaction_mol;
-        k        as frequency   (Unit="1/h"); 
+        k        as frequency   (DisplayUnit='1/h');
 in  Inlet    as corrente;
 out Outlet   as corrente;
@@ -94 +94 @@
         DEVICES
         FEED as corrente;
-        CSTR as CSTR;
+        CSTR1 as CSTR;
         PIDL as PID_Ideal_AWBT;
         PIDT as PID_Ideal_AWBT;
@@ -107 +107 @@
         
         CONNECTIONS
-        FEED to CSTR.Inlet;
+        FEED to CSTR1.Inlet;
         
         SET
 #       Parâmetros do CSTR"
-        CSTR.R   = 8.3144 * "kJ/kmol/K";
-        CSTR.U   = 300 * "kJ/h/m^2/K";
-        CSTR.ro  = 55.56 * "kmol/m^3";
-        CSTR.Cp  = 70*"kJ/kmol/K";
-        CSTR.Hr  = -7000 * "kJ/kmol";
-        CSTR.Ea  = 6e4 * "kJ/kmol";
-        CSTR.ko  = 89 * "1/s";
-        CSTR.A   = 8 * "m^2";
-        CSTR.At  = 25 * "m^2";
+        CSTR1.R   = 8.3144 * 'kJ/kmol/K';
+        CSTR1.U   = 300 * 'kJ/h/m^2/K';
+        CSTR1.ro  = 55.56 * 'kmol/m^3';
+        CSTR1.Cp  = 70*'kJ/kmol/K';
+        CSTR1.Hr  = -7000 * 'kJ/kmol';
+        CSTR1.Ea  = 6e4 * 'kJ/kmol';
+        CSTR1.ko  = 89 * '1/s';
+        CSTR1.A   = 8 * 'm^2';
+        CSTR1.At  = 25 * 'm^2';
         
         EQUATIONS
 
         "Equações do PID para controle de nível"
-        L_ad*(Lmax-Lmin)=CSTR.h-Lmin;
+        L_ad*(Lmax-Lmin)=CSTR1.h-Lmin;
         PIDL.Ports.input=L_ad;
 
         "Equações do PID para controle de temperatura"
-        T_ad*(Tmax-Tmin)=CSTR.T-Tmin;
+        T_ad*(Tmax-Tmin)=CSTR1.T-Tmin;
         PIDT.Ports.input=T_ad;
         
         "Variáveis manipulada"
-        CSTR.Cv  = 2.2136 * "m^2.5/h" * (1 - PIDL.Ports.output);
-        CSTR.Tw  = PIDT.Ports.output*(Tmax-Tmin)+Tmin;
+        CSTR1.Cv  = 2.2136 * 'm^2.5/h' * (1 - PIDL.Ports.output);
+        CSTR1.Tw  = PIDT.Ports.output*(Tmax-Tmin)+Tmin;
         
         #distúrbio regulatório
         if time<1.6e5 then      
-                FEED.T = 300 * "K";
+                FEED.T = 300 * 'K';
         else
-                FEED.T = 285 * "K";
+                FEED.T = 285 * 'K';
         end
 
@@ -152 +152 @@
         PIDL.Options.autoMan=0;
         PIDL.Parameters.gain=20;
-        PIDL.Parameters.intTime=5*"h";
-        PIDL.Parameters.derivTime=0*"s";
+        PIDL.Parameters.intTime=5*'h';
+        PIDL.Parameters.derivTime=0*'s';
         PIDL.Ports.setPoint=0.55;
-        PIDL.Parameters.tau=1*"s";
-        PIDL.Parameters.tauSet=1*"s";
+        PIDL.Parameters.tau=1*'s';
+        PIDL.Parameters.tauSet=1*'s';
         
         PIDT.Parameters.bias = 0;
@@ -166 +166 @@
         PIDT.Options.autoMan=0;
         PIDT.Parameters.gain=40;
-        PIDT.Parameters.intTime=5*"h";
-        PIDT.Parameters.derivTime=1*"h";
+        PIDT.Parameters.intTime=5*'h';
+        PIDT.Parameters.derivTime=1*'h';
         PIDT.Ports.setPoint=0.85;
-        PIDT.Parameters.tau=1*"s";
-        PIDT.Parameters.tauSet=1*"s";   
+        PIDT.Parameters.tau=1*'s';
+        PIDT.Parameters.tauSet=1*'s';   
         
         "Valores limites para normalizações"
-        Lmax=5*"m";
-        Lmin=0*"m";
-        Tmax=700*"K";
-        Tmin=230*"K";   
+        Lmax=5*'m';
+        Lmin=0*'m';
+        Tmax=700*'K';
+        Tmin=230*'K';   
         
         "Variáveís da corrente de alimentação"
-        FEED.Ca = 300 * "kmol/m^3";
-        FEED.F = 3.5 * "m^3/h";
+        FEED.Ca = 300 * 'kmol/m^3';
+        FEED.F = 3.5 * 'm^3/h';
 
         INITIAL
-        CSTR.Ca = 50 * "kmol/m^3";
-        CSTR.h = 2.5 * "m";
-        CSTR.T = 650 * "K";
+        CSTR1.Ca = 50 * 'kmol/m^3';
+        CSTR1.h = 2.5 * 'm';
+        CSTR1.T = 650 * 'K';
         
         OPTIONS
-        time = [0:0.1:1 1:1:100] * "h";
+        TimeStep = 0.1;
+        TimeEnd = 100;
+        TimeUnit = 'h';
 end