Added code snippets to Muscle Modeling tutorial to enable testing (#47)

* Added code snippets to Muscle Modeling tutorial to enable testing

* Added code snippets to Muscle Modeling lesson 5+6+7

---------

Co-authored-by: Jacob Hilmar Adamsen <jha@anybodytech.com>

Author: JacobAdamsen

A_Getting_started_anyscript/index.md

@@ -1,6 +1,7 @@
 ::: {rst-class} break
 :::
 
+(How_to_write_AnyScript)=
 # How to write AnyScript
 
 This tutorial introduces the AnyScript modeling language, which is used

Muscle_modeling/Downloads/MuscleDemo.4.any

@@ -15,7 +15,7 @@ Main = {
     };  // Global reference frame
 
     // Define one simple segment
-AnySeg Arm = {
+    AnySeg Arm = {
       r0 = {0.500000, 0.000000, 0.000000};
       Mass = 1.000000;
       Jii = {0.100000, 1.000000, 1.000000}*0.03;
@@ -43,7 +43,7 @@ AnySeg Arm = {
       DriverPos = {-10*pi/180};
       DriverVel = {80*pi/180};
       AnyRevoluteJoint &Jnt = .Jnt;
-      Reaction.Type = {0};
+      Reaction.Type = {Off};
     };
     AnyMuscleModel SimpleModel = 
     {
@@ -57,10 +57,10 @@ AnySeg Arm = {
       AnyRefFrame &Ins = .Arm.M1Insertion;
       AnyDrawMuscle drw = 
       {
-        RGB = {1, 0, 0};
+        //RGB = {1, 0, 0};
         //Opacity = 1;
         //DrawOnOff = 1;
-        Bulging = 1;
+        Bulging = 2;
         ColorScale = 1;
         RGBColorScale = {0, 0, 1};
         MaxStress = 250000;

Muscle_modeling/Downloads/MuscleDemo.5-2.any

@@ -80,7 +80,7 @@ Main = {
 
     AnyMuscleModel2ELin Model2 = {
       F0 = 200;
-      Lfbar = 0.3;
+      Lfbar = 0.2;
       Lt0 = 0.3;
       Epsilonbar = 0.05;
       V0 = -0.3;
@@ -95,7 +95,7 @@ Main = {
       AnyRefFrame &Ins = .Arm.M1Insertion;
       SPLine.StringMesh = 20;
       AnyDrawMuscle drw = {
-        Bulging = 0;
+        Bulging = 2;
         ColorScale = 1;
         MaxStress = 250000;
       };
@@ -109,7 +109,7 @@ Main = {
       SPLine.StringMesh = 20;
       SPLine.InitWrapPosVectors = {{-0.2, -0.2, 0},{-0.05,-0.2, 0}};
       AnyDrawMuscle drw = {
-        Bulging = 0;
+        Bulging = 2;
         ColorScale = 1;        
         MaxStress = 250000;
       };

Muscle_modeling/Downloads/MuscleDemo.5.any

@@ -78,7 +78,7 @@ Main = {
       AnyRefFrame &Ins = .Arm.M1Insertion;
       SPLine.StringMesh = 20;
       AnyDrawMuscle drw = {
-        Bulging = 4;
+        Bulging = 2;
         ColorScale = 1;
         MaxStress = 250000;
       };
@@ -92,7 +92,7 @@ Main = {
       SPLine.StringMesh = 20;
       SPLine.InitWrapPosVectors = {{-0.2, -0.2, 0},{-0.05,-0.2, 0}};
       AnyDrawMuscle drw = {
-        Bulging = 4;
+        Bulging = 2;
         ColorScale = 1;        
         MaxStress = 250000;
       };

Muscle_modeling/Downloads/MuscleDemo.6.any

@@ -102,8 +102,8 @@ Main = {
   // The study: Operations to be performed on the model
   AnyBodyStudy ArmStudy = {
     AnyFolder &Model = .ArmModel;
-//    InverseDynamics.Criterion.Type = MR_MinMaxStrict;
-InverseDynamics.Criterion.Type = MR_Quadratic;
+    //InverseDynamics.Criterion.Type = MR_MinMaxStrict;
+    InverseDynamics.Criterion.Type = MR_Quadratic;
     Gravity = {0.0, -9.81, 0.0};
   };
 

Muscle_modeling/Snippets/lesson1/snip.Muscles.main-1.any

@@ -0,0 +1,49 @@
+//# BEGIN SNIPPET 1
+// This is a very simple model for demonstration of muscle modeling
+Main = {
+
+  AnyFolder MyModel = {
+
+    // Global Reference Frame
+    AnyFixedRefFrame GlobalRef = {
+      AnyDrawRefFrame drw = {
+        RGB = {1,0,0};
+      };
+    };  // Global reference frame
+
+    // Define one simple segment
+    AnySeg Arm = {
+      r0 = {0.500000, 0.000000, 0.000000};
+      Mass = 1.000000;
+      Jii = {0.100000, 1.000000, 1.000000}*0.03;
+      AnyRefNode Jnt = {
+        sRel = {-0.5, 0.0, 0};
+      };      
+      AnyDrawSeg drw = {};
+    };
+    
+    // Attach the segment to ground by a revolute joint
+    AnyRevoluteJoint Jnt = {
+      AnyRefFrame &ref1 = .GlobalRef;
+      AnyRefFrame &ref2 = .Arm.Jnt;
+      Axis = z;
+    };
+
+    // Drive the revolute joint at constant velocity
+    AnyKinEqSimpleDriver Drv = {
+      DriverPos = {-10*pi/180};
+      DriverVel = {40*pi/180};
+      AnyRevoluteJoint &Jnt = .Jnt;
+      Reaction.Type = {Off};
+    };
+    
+  }; // MyModel
+
+  // The study: Operations to be performed on the model
+  AnyBodyStudy MyStudy = {
+    AnyFolder &Model = .MyModel;
+    InverseDynamics.Criterion.Type = MR_MinMaxStrict;
+    Gravity = {0.0, -9.81, 0.0};
+  };
+};  // Main
+//# END SNIPPET 1

Muscle_modeling/Snippets/lesson1/snip.Muscles.main-2.any

@@ -0,0 +1,58 @@
+//expect_errors = ["<ObjectName>", "Model loading skipped"]
+
+// This is a very simple model for demonstration of muscle modeling
+Main = {
+
+  AnyFolder MyModel = {
+
+    // Global Reference Frame
+    AnyFixedRefFrame GlobalRef = {
+      AnyDrawRefFrame drw = {
+        RGB = {1,0,0};
+      };
+    };  // Global reference frame
+
+    // Define one simple segment
+    AnySeg Arm = {
+      r0 = {0.500000, 0.000000, 0.000000};
+      Mass = 1.000000;
+      Jii = {0.100000, 1.000000, 1.000000}*0.03;
+      AnyRefNode Jnt = {
+        sRel = {-0.5, 0.0, 0};
+      };      
+      AnyDrawSeg drw = {};
+    };
+    
+    // Attach the segment to ground by a revolute joint
+    AnyRevoluteJoint Jnt = {
+      AnyRefFrame &ref1 = .GlobalRef;
+      AnyRefFrame &ref2 = .Arm.Jnt;
+      Axis = z;
+    };
+
+//# BEGIN SNIPPET 1
+// Drive the revolute joint at constant velocity
+AnyKinEqSimpleDriver Drv = {
+  DriverPos = {-10*pi/180};
+  DriverVel = {40*pi/180};
+  AnyRevoluteJoint &Jnt = .Jnt;
+  Reaction.Type = {Off};
+};
+
+§AnyMuscleModel <ObjectName> = 
+{
+  F0 = 0.0;
+  //Lf0 = 0.0;
+  //Vol0 = 0.0;
+};§
+//# END SNIPPET 1
+    
+  }; // MyModel
+
+  // The study: Operations to be performed on the model
+  AnyBodyStudy MyStudy = {
+    AnyFolder &Model = .MyModel;
+    InverseDynamics.Criterion.Type = MR_MinMaxStrict;
+    Gravity = {0.0, -9.81, 0.0};
+  };
+};  // Main

Muscle_modeling/Snippets/lesson1/snip.Muscles.main-3.any

@@ -0,0 +1,74 @@
+//expect_errors = ["<ObjectName>", "Model loading skipped"]
+
+// This is a very simple model for demonstration of muscle modeling
+Main = {
+
+  AnyFolder MyModel = {
+
+    // Global Reference Frame
+    AnyFixedRefFrame GlobalRef = {
+      AnyDrawRefFrame drw = {
+        RGB = {1,0,0};
+      };
+    };  // Global reference frame
+
+    // Define one simple segment
+    AnySeg Arm = {
+      r0 = {0.500000, 0.000000, 0.000000};
+      Mass = 1.000000;
+      Jii = {0.100000, 1.000000, 1.000000}*0.03;
+      AnyRefNode Jnt = {
+        sRel = {-0.5, 0.0, 0};
+      };      
+      AnyDrawSeg drw = {};
+    };
+    
+    // Attach the segment to ground by a revolute joint
+    AnyRevoluteJoint Jnt = {
+      AnyRefFrame &ref1 = .GlobalRef;
+      AnyRefFrame &ref2 = .Arm.Jnt;
+      Axis = z;
+    };
+
+    // Drive the revolute joint at constant velocity
+    AnyKinEqSimpleDriver Drv = {
+      DriverPos = {-10*pi/180};
+      DriverVel = {40*pi/180};
+      AnyRevoluteJoint &Jnt = .Jnt;
+      Reaction.Type = {Off};
+    };
+//# BEGIN SNIPPET 1
+AnyMuscleModel §SimpleModel§ = {
+  F0 = §100§;
+  //Lf0 = 0.0;
+  //Vol0 = 0.0;
+};
+//# END SNIPPET 1
+
+//# BEGIN SNIPPET 2
+// End of AnyMuscleModel
+
+§AnyViaPointMuscle <ObjectName> = 
+{
+  //viewForce.Visible = Off;
+  //MetabModel = Null;
+  //FatigueModel = Null;
+  //MuscleModel = Null;
+  //viewMuscle.Visible = Off;
+  AnyMuscleModel &<Insert name0> = <Insert object reference (or full object definition)>;
+  AnyRefFrame &<Insert name0> = <Insert object reference (or full object definition)>;
+  AnyRefFrame &<Insert name1> = <Insert object reference (or full object definition)>;
+  //AnyRefFrame &<Insert name2> = <Insert object reference (or full object definition)>;
+  //AnyRefFrame &<Insert name3> = <Insert object reference (or full object definition)>; You can make any number of AnyRefFrame objects!
+};§
+//# END SNIPPET 2
+    
+  }; // MyModel
+
+  // The study: Operations to be performed on the model
+  AnyBodyStudy MyStudy = {
+    AnyFolder &Model = .MyModel;
+    InverseDynamics.Criterion.Type = MR_MinMaxStrict;
+    Gravity = {0.0, -9.81, 0.0};
+  };
+};  // Main

Muscle_modeling/Snippets/lesson1/snip.Muscles.main-4.any

@@ -0,0 +1,70 @@
+//expect_errors = ["<Insert", "Model loading skipped"]
+
+// This is a very simple model for demonstration of muscle modeling
+Main = {
+
+  AnyFolder MyModel = {
+
+//# BEGIN SNIPPET 2
+// Global Reference Frame
+AnyFixedRefFrame GlobalRef = {
+  AnyDrawRefFrame drw = {
+    RGB = {1,0,0};
+  };
+  §AnyRefNode M1Origin = {
+    sRel = {0.0, 0.1, 0};
+  };§
+};  // Global reference frame
+
+// Define one simple segment
+AnySeg Arm = {
+  r0 = {0.500000, 0.000000, 0.000000};
+  Mass = 1.000000;
+  Jii = {0.100000, 1.000000, 1.000000}*0.03;
+  AnyRefNode Jnt = {
+    sRel = {-0.5, 0.0, 0};
+  };
+  §AnyRefNode M1Insertion = {
+    sRel = {0.0, 0.1, 0};
+  };§    
+  AnyDrawSeg drw = {};
+};
+//# END SNIPPET 2
+    
+    // Attach the segment to ground by a revolute joint
+    AnyRevoluteJoint Jnt = {
+      AnyRefFrame &ref1 = .GlobalRef;
+      AnyRefFrame &ref2 = .Arm.Jnt;
+      Axis = z;
+    };
+
+    // Drive the revolute joint at constant velocity
+    AnyKinEqSimpleDriver Drv = {
+      DriverPos = {-10*pi/180};
+      DriverVel = {40*pi/180};
+      AnyRevoluteJoint &Jnt = .Jnt;
+      Reaction.Type = {Off};
+    };
+//# BEGIN SNIPPET 1
+AnyMuscleModel SimpleModel = {
+  F0 = 100;
+  //Lf0 = 0.0;
+  //Vol0 = 0.0;
+};
+
+AnyMuscleViaPoint §Muscle1§ = {
+  AnyMuscleModel &§Model§ = §.SimpleModel§;
+  AnyRefFrame &<Insert name0> = <Insert object reference (or full object definition)>;
+  AnyRefFrame &<Insert name1> = <Insert object reference (or full object definition)>;
+};
+//# END SNIPPET 1
+    
+  }; // MyModel
+
+  // The study: Operations to be performed on the model
+  AnyBodyStudy MyStudy = {
+    AnyFolder &Model = .MyModel;
+    InverseDynamics.Criterion.Type = MR_MinMaxStrict;
+    Gravity = {0.0, -9.81, 0.0};
+  };
+};  // Main