Changeset 1003 for trunk/sample/controllers

Timestamp:

Jun 11, 2017, 2:37:15 PM (6 years ago)

Author:

Argimiro Resende Secchi

Message:

Fixed sintaxe problem of old version.

File:

• trunk/sample/controllers/CSTR_noniso_pid.mso (modified) (4 diffs)

trunk/sample/controllers/CSTR_noniso_pid.mso

@@ -105 +105 @@
 # Process with controlled CSTR and multiple steady-states
 FlowSheet CSTR_controller
+        PARAMETERS
+        Lmin as length;
+        Lmax as length;
+        Tmin as temperature;
+        Tmax as temperature;
 
         DEVICES
@@ -113 +118 @@
  
         VARIABLES
-        L_ad as Real;
-        Lmin as length;
-        Lmax as length;
-        T_ad as Real;
-        Tmin as temperature;
-        Tmax as temperature;
         Lsp as length;
         Tsp as temperature;
+out     TI      as control_signal;
+out     LI      as control_signal;
         
         CONNECTIONS
         FEED to CSTR1.Inlet;
-        
+        LI       to PIDL.Input;
+        TI       to PIDT.Input;
+
         SET
 #       CSTR Parameters
@@ -136 +139 @@
         CSTR1.D   = 3.2 * 'm';
         CSTR1.Cv  = 2.7 * 'm^2.5/h';
-        
-        PIDL.PID_Select = "Ideal_AWBT";
-        PIDT.PID_Select = "Ideal_AWBT";
-        
-        EQUATIONS
-
-        "Dimensionless level to connect PID"
-        L_ad*(Lmax-Lmin)=CSTR1.h-Lmin;
-        PIDL.Ports.input=L_ad;
-
-        "Dimensionless temperature to connect PID"
-        T_ad*(Tmax-Tmin)=CSTR1.T-Tmin;
-        PIDT.Ports.input=T_ad;
-        
-        "Manipulated Variables"
-        CSTR1.x  = PIDL.Ports.output;
-        CSTR1.Tw  = PIDT.Ports.output*(Tmax-Tmin)+Tmin;
-        
-#   Level control: PID parameters
-        PIDL.Parameters.bias=0;
-        PIDL.Parameters.alpha=0.1;
-        PIDL.Options.action=-1;
-        PIDL.Parameters.gamma=1;
-        PIDL.Parameters.beta=1;
-        PIDL.Options.clip=1;
-        PIDL.Options.autoMan=0;
-        PIDL.Parameters.gain=1;
-        PIDL.Parameters.intTime=2.5*'h';
-        PIDL.Parameters.derivTime=0*'s';
-        PIDL.Ports.setPoint=(Lsp - Lmin)/(Lmax - Lmin);
-        PIDL.Parameters.tau=1*'s';
-        PIDL.Parameters.tauSet=1*'s';
-
-#   Temperature control: PID parameters
-        PIDT.Parameters.bias = 0;
-        PIDT.Parameters.alpha=0.1;
-        PIDT.Options.action=1;
-        PIDT.Parameters.gamma=1;
-        PIDT.Parameters.beta=1;
-        PIDT.Options.clip=1;
-        PIDT.Options.autoMan=0;
-        PIDT.Parameters.gain=1;
-        PIDT.Parameters.intTime=2.5*'h';
-        PIDT.Parameters.derivTime=1*'h';
-        PIDT.Ports.setPoint=(Tsp - Tmin)/(Tmax - Tmin);
-        PIDT.Parameters.tau=1*'s';
-        PIDT.Parameters.tauSet=1*'s';   
-        
-        "Operating range for control variables"
+
+#       Operating range for control variables
         Lmax=5*'m';
         Lmin=0*'m';
         Tmax=700*'K';
-        Tmin=230*'K';   
+        Tmin=230*'K';
+
+        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';
+        PIDL.MinInput=Lmin/'m';
+        PIDL.MaxInput=Lmax/'m';
+
+#   Temperature control: PID parameters
+        PIDT.bias = 0;
+        PIDT.alpha=0.1;
+        PIDT.Action="Reverse";
+        PIDT.gamma=1;
+        PIDT.beta=1;
+        PIDT.Clip="Clipped";
+        PIDT.Mode="Automatic";
+        PIDT.gain=1;
+#       PIDT.gain=15; # for a setpoint of 523K or a faster response
+        PIDT.intTime=2.5*'h';
+        PIDT.derivTime=1*'h';
+        PIDT.tau=1*'s';
+        PIDT.tauSet=1*'s';
+        PIDT.MinInput=Tmin/'K';
+        PIDT.MaxInput=Tmax/'K';
+        PIDT.MinOutput=Tmin/'K';
+        PIDT.MaxOutput=Tmax/'K';
+        
+        EQUATIONS
+
+        "Dimensionless level to connect PID"
+        LI=CSTR1.h/'m';
+
+        "Dimensionless temperature to connect PID"
+        TI=CSTR1.T/'K';
+
+        "Manipulated Variables"
+        CSTR1.x  = PIDL.Output;
+        CSTR1.Tw  = PIDT.Output*'K';
+
+        "Level setpoint"
+        PIDL.SetPoint=Lsp/'m';
+
+        "Temperature setpoint"
+        PIDT.SetPoint=Tsp/'K';
 
         "Feed Stream"
@@ -206 +218 @@
         else
           Tsp = 400 * 'K';
+#         Tsp = 523 * 'K'; # near unstable steady-state (change also the controller's gain)
         end