[bfdf16]: / Body / AAUHuman / Leg / MusPar.any

Download this file

610 lines (559 with data), 33.0 kB

// This include file contains all the muscle parameter 
// data of the lower extremity. Please note that it is 
// necessary to calibrate the model. This data is not 
// calibrated!

// A lot of information in this file is obtained from:
// http: // isb.ri.ccf.org/data/delp/Muscle_parameter_table.txt

// Mark de Zee

AnyVar StrengthScaleShank = .StrengthScaling.Shank.StrengthScale;
AnyVar StrengthScaleThigh = .StrengthScaling.Thigh.StrengthScale;
AnyVar StrengthScalePelvis = ....Scaling.StrengthScaling.Pelvis.StrengthScale;
AnyVar FiberLengthScaleShank = .FiberLengthScaling.Shank.FiberLengthScale;
AnyVar FiberLengthScaleThigh = .FiberLengthScaling.Thigh.FiberLengthScale;


AnyMuscleModel3E SoleusPar = {
  F0 = .StrengthScaleShank*.StressIndex*2830; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleShank* 0.04; // Optimum fiber length
  Lt0 = 0.3773608; // Tendon slack length
  Gamma0 = 25*(3.14/180); // Pennation angle in radian: similar to Delp's soleus
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5 ; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E GastrocnemiusPar = {
  F0 = .StrengthScaleShank*.StressIndex*1605; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleShank* 0.06; // Optimum fiber length
  Lt0 = 0.4484168; // Tendon slack length
  Gamma0 = 12*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5 ; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E FlexorDigitorumLongusPar = {
  F0 = .StrengthScaleShank*.StressIndex*310; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleShank* 0.06; // Optimum fiber length
  Lt0 =  0.4288993; // Tendon slack length
  Gamma0 = 7*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5 ; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E FlexorHallucisLongusPar = {
  F0 = .StrengthScaleShank*.StressIndex*320; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleShank* 0.04; // Optimum fiber length
  Lt0 = 0.3915658; // Tendon slack length
  Gamma0 = 10*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5 ; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E TibialisPosteriorPar = {
  F0 = .StrengthScaleShank*.StressIndex*1270; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleShank* 0.03; // Optimum fiber length
  Lt0 =  0.4528486; // Tendon slack length
  Gamma0 = 12*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  // PEFactor = 2.5 ; // Parameter for influence of parallel stiffness
  PEFactor = 2.5 ; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E PeroneusBrevisPar = {
  F0 = .StrengthScaleShank*.StressIndex*350; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleShank* 0.05; // Optimum fiber length
  Lt0 = 0.2403773; // Tendon slack length
  Gamma0 = 5*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5 ; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E TibialisAnteriorPar = {
  F0 = .StrengthScaleShank*.StressIndex*600; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleShank* 0.056; // Optimum fiber length. Delp has 0.098
  Lt0 = 0.2499643; // Tendon slack length
  Gamma0 = 5*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5 ; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E ExtensorDigitorumLongusPar = {
  F0 = .StrengthScaleShank*.StressIndex*340; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleShank* 0.10; // Optimum fiber length
  Lt0 =  0.4399119; // Tendon slack length
  Gamma0 = 8*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5 ; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E ExtensorHallucisLongusPar = {
  F0 = .StrengthScaleShank*.StressIndex*110; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleShank* 0.11; // Optimum fiber length
  Lt0 = 0.2389452; // Tendon slack length
  Gamma0 = 6*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor =2; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E VastusLateralisPar = {
  F0 = .StrengthScaleThigh*.StressIndex*1870; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleThigh* 0.084; // Optimum fiber length
  Lt0 = 0.271776; // Tendon slack length
  Gamma0 = 5*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5 ; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E VastusMedialisPar = {
  F0 = .StrengthScaleThigh*.StressIndex*1295; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleThigh* 0.089; // Optimum fiber length
  Lt0 =  0.3018; // Tendon slack length
  Gamma0 = 5*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5 ; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E VastusIntermediusPar = {
  F0 = .StrengthScaleThigh*.StressIndex*1235; // Maximum force output at optimum fiber length. Delp had 1365 N
  Lf0 = 1.2*.FiberLengthScaleThigh* 0.087; // Optimum fiber length
  Lt0 = 0.2406; // Tendon slack length
  Gamma0 = 3*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5 ; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E RectusFemorisPar = {
  F0 = .StrengthScaleThigh*.StressIndex*780; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleThigh* 0.084; // Optimum fiber length
  Lt0 = 0.5170; // Tendon slack length. Delp has 0.35
  Gamma0 = 5*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5 ; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E SemitendinosusPar = {
  F0 = .StrengthScaleThigh*.StressIndex*330; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleThigh* 0.201; // Optimum fiber length
  Lt0 = 0.3297697; // Tendon slack length
  Gamma0 = 5*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5 ; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E SemimembranosusPar = {
  F0 = .StrengthScaleThigh*.StressIndex*1030; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleThigh* 0.08; // Optimum fiber length
  Lt0 =  0.3879423; // Tendon slack length
  Gamma0 = 15*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5 ; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E BicepsFemorisCaputLongumPar = {
  F0 = .StrengthScaleThigh*.StressIndex*720; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleThigh* 0.11; // Optimum fiber length
  Lt0 = 0.369793; // Tendon slack length
  Gamma0 = 0*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5 ; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E BicepsFemorisCaputBrevePar = {
  F0 = .StrengthScaleThigh*.StressIndex*400; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleThigh* 0.173; // Optimum fiber length
  Lt0 = 0.07296989; // Tendon slack length
  Gamma0 = 23*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5 ; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E SartoriusPar = {
  F0 = .StrengthScaleThigh*.StressIndex*104; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleThigh* 0.579; // Optimum fiber length
  Lt0 =  0.005;// Tendon slack length
  Gamma0 = 0*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // slow twitch Factor for determining shortening speed at optimum length
  K2 = 0; // fast twitch (zero when no info available) determines shortening speed at optimum length
  Fcfast = 0.4; // percentage fast to slow Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5 ; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E GracilisPar = {
  F0 = .StrengthScaleThigh*.StressIndex*108; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleThigh* 0.352; // Optimum fiber length
  Lt0 = 0.1420071; // Tendon slack length
  Gamma0 = 3*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // slow twitch Factor for determining shortening speed at optimum length
  K2 = 0; // fast twitch (zero when no info available) determines shortening speed at optimum length
  Fcfast = 0.4; // percentage fast to slow Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5 ; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E IliopsoasPar = {
  F0 = .StrengthScalePelvis*.StressIndex*1000; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleThigh* 0.1; // Optimum fiber length
  Lt0 =  0.0600; // Tendon slack length
  Gamma0 = 7.5*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5 ; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E GluteusMinimus1Par = {
  F0 = .StrengthScalePelvis*.StressIndex*180; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleThigh* 0.068; // Optimum fiber length
  Lt0 =  0.04375294; // Tendon slack length
  Gamma0 = 10*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5  ; //PEFactor = 2.5 ; //PEFactor = 2.5; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E GluteusMinimus2Par = {
  F0 = .StrengthScalePelvis*.StressIndex*190; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleThigh* 0.056; // Optimum fiber length
  Lt0 = 0.07103281; // Tendon slack length
  Gamma0 = 0*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5  ; //PEFactor = 2.5; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E GluteusMinimus3Par = {
  F0 = .StrengthScalePelvis*.StressIndex*215; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleThigh* 0.038; // Optimum fiber length
  Lt0 = 0.09763591; // Tendon slack length
  Gamma0 = 1*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5  ; //PEFactor = 2.5; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E GluteusMedius1Par = {
  F0 = .StrengthScalePelvis*.StressIndex*550; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleThigh* 0.054; // Optimum fiber length
  Lt0 = 0.0758717; // Tendon slack length
  Gamma0 = 8*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5  ; //PEFactor = 2.5; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E GluteusMedius2Par = {
  F0 = .StrengthScalePelvis*.StressIndex*380; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleThigh* 0.084; // Optimum fiber length
  Lt0 = 0.08978282; // Tendon slack length
  Gamma0 = 0*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5  ; //PEFactor = 2.5; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E GluteusMedius3Par = {
  F0 = .StrengthScalePelvis*.StressIndex*435; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleThigh* 0.065; // Optimum fiber length
  Lt0 = 0.1219481; // Tendon slack length
  Gamma0 = 19*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5  ; //PEFactor = 2.5; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E GluteusMaximus1Par = {
  F0 = .StrengthScalePelvis*.StressIndex*382; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleThigh* 0.142; // Optimum fiber length
  Lt0 =  0.4643645; // 0.125; // Tendon slack length - includes fascia lata length as in TFL
  Gamma0 = 5*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5 ; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E GluteusMaximus2Par = {
  F0 = .StrengthScalePelvis*.StressIndex*546; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleThigh* 0.147; // Optimum fiber length
  Lt0 =0.4471073; // 0.127; // Tendon slack length - includes fascia lata length as in TFL
  Gamma0 = 5*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5 ; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E GluteusMaximus3Par = {
  F0 = .StrengthScalePelvis*.StressIndex*368; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleThigh* 0.144; // Optimum fiber length
  Lt0 = 0.08199298; // Tendon slack length
  Gamma0 = 5*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5 ; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E TensorFasciaeLataePar = {
  F0 = .StrengthScaleThigh*.StressIndex*155; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleThigh* 0.095; // Optimum fiber length
  Lt0 =  0.4512272; // Tendon slack length must include lateral fasciae
  Gamma0 = 3*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5 ; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E PiriformisPar = {
  F0 = .StrengthScalePelvis*.StressIndex*295; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleThigh* 0.026; // Optimum fiber length
  Lt0 = 0.1216; // Tendon slack length
  Gamma0 = 10*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5 ; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E AdductorLongusPar = {
  F0 = .StrengthScaleThigh*.StressIndex*420; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleThigh* 0.138; // Optimum fiber length
  Lt0 =0.05765888; // Tendon slack length
  Gamma0 = 6*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5 ; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E AdductorMagnus1Par = {
  F0 = .StrengthScaleThigh*.StressIndex*345; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleThigh* 0.087; // Optimum fiber length
  Lt0 =  0.03893472; // Tendon slack length
  Gamma0 = 5*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5 ; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E AdductorMagnus2Par = {
  F0 = .StrengthScaleThigh*.StressIndex*310; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleThigh* 0.121; // Optimum fiber length
  Lt0 =  0.05939183; // Tendon slack length
  Gamma0 = 3*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5 ; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E AdductorMagnus3Par = {
  F0 = .StrengthScaleThigh*.StressIndex*445; // Maximum force output at optimum fiber length
  Lf0 = .FiberLengthScaleThigh* 0.131; // Optimum fiber length
  Lt0 =0.2068808; // Tendon slack length
  Gamma0 = 5*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5 ; // Parameter for influence of parallel stiffness
};



AnyMuscleModel3E QuadratusFemorisPar = {
  F0 = .StrengthScaleThigh*.StressIndex*255; // Maximum force output at optimum fiber length.
  Lf0 = .FiberLengthScaleThigh* 0.054; // Optimum fiber length
  Lt0 =0.04295;//0.024; // Tendon slack length
  Gamma0 = 8*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5 ; // Parameter for influence of parallel stiffness

};

AnyMuscleModel3E AbductorBrevisPar = {
  F0 = .StrengthScaleThigh*.StressIndex*285; // Maximum force output at optimum fiber length.
  Lf0 = .FiberLengthScaleThigh* 0.133; // Optimum fiber length
  Lt0 =0.00866;//0.02; // Tendon slack length
  Gamma0 = 0*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor =2.0 ; // Parameter for influence of parallel stiffness

};
AnyMuscleModel3E PictineusPar = {
  F0 = .StrengthScaleThigh*.StressIndex*285; // Maximum force output at optimum fiber length.
  Lf0 = .FiberLengthScaleThigh* 0.133; // Optimum fiber length
  Lt0 =0.005;//0.01; // Tendon slack length
  Gamma0 = 0*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5 ; // Parameter for influence of parallel stiffness

};

AnyMuscleModel3E ObturatoriusExternusPar = {
  F0 = .StrengthScaleThigh*.StressIndex*800; // Maximum force output at optimum fiber length.
  Lf0 = .FiberLengthScaleThigh* 0.049; // Optimum fiber length
  Lt0 =0.05152; // Tendon slack length
  Gamma0 = 5*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5 ; // Parameter for influence of parallel stiffness

};

AnyMuscleModel3E ObturatoriusInternusPar = {
  F0 = .StrengthScaleThigh*.StressIndex*800; // Maximum force output at optimum fiber length.
  Lf0 = .FiberLengthScaleThigh* 0.0474; // Optimum fiber length
  Lt0 =0.1122;//0.08; // Tendon slack length
  Gamma0 = 5*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E GemellusSuperiorPar = {
  F0 = .StrengthScaleThigh*.StressIndex*110*0.5; // Maximum force output at optimum fiber length.
  Lf0 = .FiberLengthScaleThigh* 0.024; // Optimum fiber length
  Lt0 = 0.09978118; // Tendon slack length
  Gamma0 = 10*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5 ; // Parameter for influence of parallel stiffness
};

AnyMuscleModel3E GemellusInferiorPar = {
  F0 = .StrengthScaleThigh*.StressIndex*110*0.5; // Maximum force output at optimum fiber length.
  Lf0 = .FiberLengthScaleThigh* 0.024; // Optimum fiber length
  Lt0 =0.09978; // Tendon slack length
  Gamma0 = 10*(3.14/180); // Pennation angle in radian
  Epsilon0 = 0.053; // Tendon strain at F0
  K1 = 10; // Factor for determining shortening speed at optimum length
  K2 = 0; // Factor for determining shortening speed at optimum length
  Fcfast = 0.4; // Factor for determining shortening speed at optimum length
  Jt = 3.0; // Shape parameter for the tendon stiffness
  Jpe = 3.0; // Shape parameter for the parallel stiffness
  PEFactor = 2.5 ; // Parameter for influence of parallel stiffness
};