Changeset 1118 for branches/maintenance/2.2.2/Modelica/Mechanics/Rotational.mo

Timestamp:

05/04/08 21:39:16 (4 months ago)

Author:

otter

Message:

Enhanced efficiency of friction models of Rotational and Translational library in some situations (only of peak > 1)

Files:

• branches/maintenance/2.2.2/Modelica/Mechanics/Rotational.mo (modified) (4 diffs)

branches/maintenance/2.2.2/Modelica/Mechanics/Rotational.mo

@@ -724 +724 @@
       Modelica.Mechanics.Rotational.Torque torque 
                                annotation (extent=[-65, -5; -55, 5]);
-      Modelica.Mechanics.Rotational.Clutch clutch1(peak=1.1, fn_max=20) 
+      Modelica.Mechanics.Rotational.Clutch clutch1(          fn_max=20, peak=1.1)
         annotation (extent=[-25, -5; -15, 5]);
       Modelica.Blocks.Sources.Sine sin1(amplitude=10, freqHz=5) 
@@ -732 +732 @@
       Modelica.Mechanics.Rotational.Inertia J2(J=1, initType=Modelica.Mechanics.Rotational.Types.Init.InitialState) 
                                  annotation (extent=[-5, -5; 5, 5]);
-      Modelica.Mechanics.Rotational.Clutch clutch2(peak=1.1, fn_max=20) 
+      Modelica.Mechanics.Rotational.Clutch clutch2(          fn_max=20, peak=1.1)
         annotation (extent=[15, -5; 25, 5]);
       Modelica.Mechanics.Rotational.Inertia J3(J=1, initType=Modelica.Mechanics.Rotational.Types.Init.InitialState) 
                                  annotation (extent=[35, -5; 45, 5]);
-      Modelica.Mechanics.Rotational.Clutch clutch3(peak=1.1, fn_max=20) 
+      Modelica.Mechanics.Rotational.Clutch clutch3(          fn_max=20, peak=1.1)
         annotation (extent=[55, -5; 65, 5]);
       Modelica.Mechanics.Rotational.Inertia J4(J=1, initType=Modelica.Mechanics.Rotational.Types.Init.InitialState) 
@@ -1817 +1817 @@
       constant SI.AngularAcceleration unitAngularAcceleration = 1;
       constant SI.Torque unitTorque = 1;
-      
+      Real sa0 = (tau0_max + (tau0_max - tau0))/unitTorque 
+        "Value of sa when start of forward sliding at w=0";
     equation 
       /* Friction characteristic
@@ -1824 +1825 @@
       special code shall be generated) */
       
-      startForward = pre(mode) == Stuck and (sa > tau0_max/unitTorque or pre(startForward)
-         and sa > tau0/unitTorque) or pre(mode) == Backward and w_relfric > w_small or 
-        initial() and (w_relfric > 0);
-      startBackward = pre(mode) == Stuck and (sa < -tau0_max/unitTorque or pre(
-        startBackward) and sa < -tau0/unitTorque) or pre(mode) == Forward and w_relfric <
-        -w_small or initial() and (w_relfric < 0);
-      locked = not free and not (pre(mode) == Forward or startForward or pre(
-        mode) == Backward or startBackward);
+      startForward = pre(mode) == Stuck and sa > tau0_max/unitTorque or 
+                     pre(mode) == Backward and w_relfric > w_small or 
+                     initial() and (w_relfric > 0);
+      startBackward = pre(mode) == Stuck and sa < -tau0_max/unitTorque or 
+                      pre(mode) == Forward and w_relfric < -w_small or 
+                      initial() and (w_relfric < 0);
+      locked = not free and not (pre(mode) == Forward or startForward or pre(mode) == Backward or startBackward);
       
       a_relfric = unitAngularAcceleration *
-              (if locked then 0 else if free then sa else if startForward then 
-              sa - tau0/unitTorque else if startBackward then sa + tau0/unitTorque else if pre(mode)
-         == Forward then sa - tau0/unitTorque else sa + tau0/unitTorque);
+              (if locked then 0 else if free then sa else 
+              if startForward then sa - sa0 else if startBackward then sa + sa0 else 
+              if pre(mode)== Forward then sa - sa0 else sa + sa0);
       
       /* Friction torque has to be defined in a subclass. Example for a clutch: