MuscleFunctionFactory.h

Go to the documentation of this file.

#ifndef MUSCLEFUNCTIONFACTORY_H_
#define MUSCLEFUNCTIONFACTORY_H_
/* -------------------------------------------------------------------------- *
 *   OpenSim:  SmoothSegmentedFunctionFactory.h                               *
 * -------------------------------------------------------------------------- *
 * The OpenSim API is a toolkit for musculoskeletal modeling and simulation.  *
 * See http://opensim.stanford.edu and the NOTICE file for more information.  *
 * OpenSim is developed at Stanford University and supported by the US        *
 * National Institutes of Health (U54 GM072970, R24 HD065690) and by DARPA    *
 * through the Warrior Web program.                                           *
 *                                                                            *
 * Copyright (c) 2005-2012 Stanford University and the Authors                *
 * Author(s): Matthew Millard                                                 *
 *                                                                            *
 * Licensed under the Apache License, Version 2.0 (the "License"); you may    *
 * not use this file except in compliance with the License. You may obtain a  *
 * copy of the License at http://www.apache.org/licenses/LICENSE-2.0.         *
 *                                                                            *
 * Unless required by applicable law or agreed to in writing, software        *
 * distributed under the License is distributed on an "AS IS" BASIS,          *
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.   *
 * See the License for the specific language governing permissions and        *
 * limitations under the License.                                             *
 * -------------------------------------------------------------------------- */
/*
  Update:
  This is a port of the original code so that it will work with
  the multibody code RBDL written by Martin Felis.
  
  Author:
  Matthew Millard
  
  Date:
  Nov 2015
 
*/
#include "../geometry/SmoothSegmentedFunction.h"
#include "../geometry/SegmentedQuinticBezierToolkit.h"
 
#include <cstdio>
#include <iostream>
#include <fstream>
#include <cmath>
 
 
 
namespace RigidBodyDynamics {
  namespace Addons {
  namespace Muscle{
 
class MuscleFunctionFactory
{
 
 
  public:
 
   // friend class SmoothSegmentedFunction;
 
 
  static void createFiberActiveForceLengthCurve(
                  double lce0, 
                  double lce1, 
                  double lce2, 
                  double lce3, 
                  double minActiveForceLengthValue, 
                  double plateauSlope, 
                  double curviness, 
                  const std::string& curveName,
                  RigidBodyDynamics::Addons::Geometry
                    ::SmoothSegmentedFunction&
                    smoothSegmentedFunctionToUpdate);
 
  static void createFiberForceVelocityCurve(
                  double fmaxE, 
                  double dydxC, 
                  double dydxNearC, 
                  double dydxIso, 
                  double dydxE, 
                  double dydxNearE,
                  double concCurviness, 
                  double eccCurviness,
                  const std::string& curveName,
                  RigidBodyDynamics::Addons::Geometry::
                    SmoothSegmentedFunction&
                    smoothSegmentedFunctionToUpdate);
 
  static void createFiberForceVelocityInverseCurve(
                  double fmaxE, 
                  double dydxC, 
                  double dydxNearC, 
                  double dydxIso, 
                  double dydxE, 
                  double dydxNearE,
                  double concCurviness, 
                  double eccCurviness, 
                  const std::string& muscleName,
                  RigidBodyDynamics::Addons::Geometry::
                    SmoothSegmentedFunction&
                    smoothSegmentedFunctionToUpdate);
 
  static void createFiberCompressiveForcePennationCurve(
                double phi0, 
                double kiso, 
                double curviness,   
                const std::string& curveName,
                RigidBodyDynamics::Addons::Geometry::
                  SmoothSegmentedFunction&
                  smoothSegmentedFunctionToUpdate);
 
  static void createFiberCompressiveForceCosPennationCurve(
                double cosPhi0, 
                double kiso, 
                double curviness, 
                const std::string& curveName,
                RigidBodyDynamics::Addons::Geometry::
                  SmoothSegmentedFunction&
                  smoothSegmentedFunctionToUpdate);
 
 
  static void createFiberCompressiveForceLengthCurve(
                double l0, 
                double kiso, 
                double curviness,  
                const std::string& curveName,
                RigidBodyDynamics::Addons::Geometry::
                  SmoothSegmentedFunction&
                  smoothSegmentedFunctionToUpdate);
 
  static void createFiberForceLengthCurve(
                 double eZero, 
                 double eIso,
                 double kLow, 
                 double kIso,
                 double curviness,
                 const std::string& curveName,
                 RigidBodyDynamics::Addons::Geometry::
                  SmoothSegmentedFunction&
                  smoothSegmentedFunctionToUpdate);
 
  static void createTendonForceLengthCurve(double eIso, 
                double kIso,
                double fToe, 
                double curviness,
                const std::string& curveName,
                RigidBodyDynamics::Addons::Geometry::
                  SmoothSegmentedFunction&
                  smoothSegmentedFunctionToUpdate);
 
 
  
 
};
 
}
}
}
 
#endif //MUSCLEFUNCTIONFACTORY_H_