Changeset 941 for branches/gui/sample/controllers/CSTR_noniso.mso

Timestamp:

Jul 15, 2011, 3:44:52 PM (12 years ago)

Author:

mamuller

Message:

updated samples

File:

• branches/gui/sample/controllers/CSTR_noniso.mso (modified) (2 diffs)

branches/gui/sample/controllers/CSTR_noniso.mso

@@ -21 +21 @@
 
 using "types";
+using "controllers/PIDs";
 
 const_valv as positive(Brief = "Valve Constant", Default=1,Lower=0,Upper=100, Unit='m^2.5/h');
@@ -152 +153 @@
         DAESolver(File = "dassl");
 end
+
+Model CSTR_pid as CSTR
+        
+        VARIABLES
+out     LI as control_signal (Brief="Level Indicator");
+out     TI as control_signal (Brief="Temperature Indicator");
+        
+        EQUATIONS
+        "Level sensor"
+        LI * 'm' = h;
+
+        "Temperature sensor"
+        TI * 'K' = T; 
+        
+end
+
+# Process with controlled CSTR and multiple steady-states
+FlowSheet CSTR_controller
+
+        DEVICES
+        FEED as stream_cstr;
+        CSTR1 as CSTR_pid;
+        PIDL as PID;
+        PIDT as PID;
+ 
+        VARIABLES
+        Lsp as length;
+        Tsp as temperature;
+
+
+        CONNECTIONS
+        FEED to CSTR1.Inlet;
+        CSTR1.LI to PIDL.Input;
+        CSTR1.TI to PIDT.Input;
+        
+        SET
+#       CSTR Parameters
+        CSTR1.R   = 8.3144 * 'kJ/kmol/K';
+        CSTR1.U   = 300 * 'kJ/h/m^2/K';
+        CSTR1.ro  = 1000 * 'kg/m^3';
+        CSTR1.Cp  = 4*'kJ/kg/K';
+        CSTR1.Hr  = -7000 * 'kJ/kmol';
+        CSTR1.Ea  = 6e4 * 'kJ/kmol';
+        CSTR1.ko  = 89 * '1/s';
+        CSTR1.D   = 3.2 * 'm';
+        CSTR1.Cv  = 2.7 * 'm^2.5/h';
+        
+        PIDL.PID_Select = "Ideal_AWBT";
+        PIDT.PID_Select = "Ideal_AWBT";
+        
+#   Level control: PID parameters
+        PIDL.bias=0;
+        PIDL.alpha=0.1;
+        PIDL.Action="Direct";
+        PIDL.gamma=1;
+        PIDL.beta=1;
+        PIDL.Clip="Clipped";
+        PIDL.Mode="Automatic";
+        PIDL.gain=1;
+        PIDL.intTime=2.5*'h';
+        PIDL.derivTime=0*'s';
+        PIDL.tau=1*'s';
+        PIDL.tauSet=1*'s';
+
+#   Temperature control: PID parameters
+        PIDT.bias = 0;
+        PIDT.alpha=0.1;
+        PIDT.Action="Direct";
+        PIDT.gamma=1;
+        PIDT.beta=1;
+        PIDT.Clip="Clipped";
+        PIDT.Mode="Automatic";
+        PIDT.gain=1;
+        PIDT.intTime=2.5*'h';
+        PIDT.derivTime=1*'h';
+        PIDT.tau=1*'s';
+        PIDT.tauSet=1*'s';      
+        
+#       Operating range for control variables
+        PIDL.MaxInput=5;
+        PIDL.MinInput=0;
+        PIDT.MaxInput=700;
+        PIDT.MinInput=230;      
+        PIDT.MaxOutput=700;
+        PIDT.MinOutput=230;     
+
+        EQUATIONS
+
+        "Setpoints"
+        PIDL.SetPoint * 'm' = Lsp;
+        PIDT.SetPoint * 'K' = Tsp;
+
+        "Manipulated Variables"
+        CSTR1.x = PIDL.Output;
+        CSTR1.Tw = PIDT.Output * 'K';
+
+        "Feed Stream"
+        FEED.Ca = 300 * 'kmol/m^3';
+        FEED.F = 3.5 * 'm^3/h';
+
+#   Disturbance
+        if time < 50 * 'h' then 
+                FEED.T = 300 * 'K';
+        else
+                FEED.T = 285 * 'K'; # change to 350 K to saturate controller
+        end
+
+#   Set-point changes
+        if time < 100 * 'h' then
+          Tsp = 630 * 'K';
+        else
+          Tsp = 400 * 'K';
+        end
+
+        if time < 150 * 'h' then
+          Lsp = 1.7 * 'm';
+        else
+          Lsp = 4 * 'm';
+        end
+
+        INITIAL
+        CSTR1.Ca = 50 * 'kmol/m^3';
+        CSTR1.h = 1.7 * 'm';
+        CSTR1.T = 570 * 'K';
+        
+        OPTIONS
+        TimeStep = 1;
+        TimeEnd = 250;
+        TimeUnit = 'h';
+        DAESolver(File = "dassl");
+end
+