Changeset 542 for trunk/sample/miscellaneous/sample_pend.mso

Timestamp:

Jun 21, 2008, 1:01:41 PM (15 years ago)

Author:

Argimiro Resende Secchi

Message:

Adding polar coordinate for pendulum sample.

File:

• trunk/sample/miscellaneous/sample_pend.mso (modified) (2 diffs)

trunk/sample/miscellaneous/sample_pend.mso

@@ -53 +53 @@
         SET
         g     = 9.8 * 'm/s^2';
-        L     = 1 * 'm';
+        L     = 0.9 * 'm';
 
         INITIAL
         "Initial Position x"
-        x = -1 * 'm';
+        x = 0.9 * 'm';
 
         "Initial x Velocity"
         w = 0 * 'm/s';
         
-        #"Initial Position y"  y = 0 * "m";
-        #"Initial y Velocity"  z = 0 * "m/s";
-
         OPTIONS
         TimeStep = 0.1;
-        TimeEnd = 30;
-        Integration = "original"; # original, index0 or index1
+        TimeEnd = 36;
+        Integration = "original"; # original, index0, or index1
 
         NLASolver(
@@ -75 +72 @@
         );
         DAESolver(
-                File = "dasslc",
+                #File = "dasslc",
                 #File = "dassl", # fail when "original" (high-index)
+                File = "mebdf",
+                #File = "pside",
+                #File = "sundials", # fail when "original" (high-index)
+                RelativeAccuracy = 1e-6,
+                AbsoluteAccuracy = 1e-8
+        );
+        SparseAlgebra = true;
+end
+
+
+FlowSheet pend_polar
+        PARAMETERS
+        g     as acceleration (Brief = "Gravity acceleration");
+        L     as length (Brief = "Pendulum cable length");
+        
+        VARIABLES
+        phi   as angle (Brief="Angle");
+        omega as frequency (Brief="Angular velocity", Lower=-100);
+        x     as length_delta(Brief="Position x");
+        y     as length_delta(Brief="Position y");
+        w     as velocity(Brief = "Velocity for x");
+        z     as velocity(Brief = "Velocity for y");
+        T     as Real(Brief = "Tension on cable",Default=10,Unit='1/s^2');
+        
+        EQUATIONS
+        "x Position"
+        x = L*sin(phi);
+
+        "y Position"
+        y = -L*cos(phi);
+
+        "Velocity on x"
+        w = L*omega*cos(phi);
+        
+        "Velocity on y"
+        z = L*omega*sin(phi);
+        
+        "Tension"
+        T = -g*cos(phi)/L-omega^2;
+
+        "Angular velocity"
+        diff(phi) = omega*'rad';
+
+        "Angular acceleration"
+        diff(omega) = -g*sin(phi)/L;
+
+        SET
+        g     = 9.8 * 'm/s^2';
+        L     = 0.9 * 'm';
+
+        INITIAL
+        "Initial Position x"
+        x = 0.9 * 'm';
+
+        "Initial x Velocity"
+        w = 0 * 'm/s';
+        
+        OPTIONS
+        TimeStep = 0.1;
+        TimeEnd = 36;
+        Integration = "original"; # original, index0, or index1
+
+        NLASolver(
+                RelativeAccuracy = 1e-8,
+                AbsoluteAccuracy = 1e-9
+        );
+        DAESolver(
+                #File = "dasslc",
+                File = "dassl", # fail when "original" (high-index)
                 #File = "mebdf",
                 #File = "pside",
                 #File = "sundials", # fail when "original" (high-index)
                 RelativeAccuracy = 1e-6,
-                AbsoluteAccuracy = 1e-8,
-                EventAccuracy = 1e-3
+                AbsoluteAccuracy = 1e-8
         );
-        SparseAlgebra = true; # dense or sparse
+        SparseAlgebra = true;
 end