|
| Millard2012EquilibriumMuscle () |
| Default constructor. More...
|
|
| Millard2012EquilibriumMuscle (const std::string &aName, double aMaxIsometricForce, double aOptimalFiberLength, double aTendonSlackLength, double aPennationAngle) |
| Constructs a functional muscle using default curves and activation model parameters. More...
|
|
bool | getUseFiberDamping () const |
|
double | getFiberDamping () const |
|
double | getDefaultActivation () const |
|
double | getDefaultFiberLength () const |
|
double | getActivationTimeConstant () const |
|
double | getDeactivationTimeConstant () const |
|
double | getMinimumActivation () const |
|
const ActiveForceLengthCurve & | getActiveForceLengthCurve () const |
|
const ForceVelocityCurve & | getForceVelocityCurve () const |
|
const FiberForceLengthCurve & | getFiberForceLengthCurve () const |
|
const TendonForceLengthCurve & | getTendonForceLengthCurve () const |
|
const MuscleFixedWidthPennationModel & | getPennationModel () const |
|
const MuscleFirstOrderActivationDynamicModel & | getActivationModel () const |
|
double | getMinimumFiberLength () const |
|
double | getMinimumFiberLengthAlongTendon () const |
|
double | getTendonForceMultiplier (SimTK::State &s) const |
|
double | getFiberStiffnessAlongTendon (const SimTK::State &s) const |
|
double | getFiberVelocity (const SimTK::State &s) const |
|
double | getActivationDerivative (const SimTK::State &s) const |
|
double | getPassiveFiberElasticForce (const SimTK::State &s) const |
| get the portion of the passive fiber force generated by the elastic element only (N) More...
|
|
double | getPassiveFiberElasticForceAlongTendon (const SimTK::State &s) const |
| get the portion of the passive fiber force generated by the elastic element only, projected onto the tendon direction (N) More...
|
|
double | getPassiveFiberDampingForce (const SimTK::State &s) const |
| get the portion of the passive fiber force generated by the damping element only (N) More...
|
|
double | getPassiveFiberDampingForceAlongTendon (const SimTK::State &s) const |
| get the portion of the passive fiber force generated by the damping element only, projected onto the tendon direction (N) More...
|
|
void | setMuscleConfiguration (bool ignoreTendonCompliance, bool ignoreActivationDynamics, double dampingCoefficient) |
|
void | setFiberDamping (double dampingCoefficient) |
|
void | setDefaultActivation (double activation) |
|
void | setActivation (SimTK::State &s, double activation) const override |
|
void | setDefaultFiberLength (double fiberLength) |
|
void | setActivationTimeConstant (double activationTimeConstant) |
|
void | setDeactivationTimeConstant (double deactivationTimeConstant) |
|
void | setMinimumActivation (double minimumActivation) |
|
void | setActiveForceLengthCurve (ActiveForceLengthCurve &aActiveForceLengthCurve) |
|
void | setForceVelocityCurve (ForceVelocityCurve &aForceVelocityCurve) |
|
void | setFiberForceLengthCurve (FiberForceLengthCurve &aFiberForceLengthCurve) |
|
void | setTendonForceLengthCurve (TendonForceLengthCurve &aTendonForceLengthCurve) |
|
void | setFiberLength (SimTK::State &s, double fiberLength) const |
|
double | computeActuation (const SimTK::State &s) const override final |
|
void | computeInitialFiberEquilibrium (SimTK::State &s) const override |
| Computes the fiber length such that the fiber and tendon are developing the same force, distributing the velocity of the entire musculotendon actuator between the fiber and tendon according to their relative stiffnesses. More...
|
|
void | computeFiberEquilibrium (SimTK::State &s, bool solveForVelocity=false) const |
| Computes the fiber length such that the fiber and tendon are developing the same force, either assuming muscle-tendon velocity as provided by the state or zero as designated by the useZeroVelocity flag. More...
|
|
void | extendPostScale (const SimTK::State &s, const ScaleSet &scaleSet) override |
| Adjust the properties of the muscle after the model has been scaled. More...
|
|
|
const double & | get_fiber_damping () const |
| Get the value of the fiber_damping property. More...
|
|
double & | upd_fiber_damping () |
| Get a writable reference to the fiber_damping property. More...
|
|
void | set_fiber_damping (const double &value) |
| Set the value of the fiber_damping property. More...
|
|
const double & | get_default_activation () const |
| Get the value of the default_activation property. More...
|
|
double & | upd_default_activation () |
| Get a writable reference to the default_activation property. More...
|
|
void | set_default_activation (const double &value) |
| Set the value of the default_activation property. More...
|
|
const double & | get_default_fiber_length () const |
| Get the value of the default_fiber_length property. More...
|
|
double & | upd_default_fiber_length () |
| Get a writable reference to the default_fiber_length property. More...
|
|
void | set_default_fiber_length (const double &value) |
| Set the value of the default_fiber_length property. More...
|
|
const double & | get_activation_time_constant () const |
| Get the value of the activation_time_constant property. More...
|
|
double & | upd_activation_time_constant () |
| Get a writable reference to the activation_time_constant property. More...
|
|
void | set_activation_time_constant (const double &value) |
| Set the value of the activation_time_constant property. More...
|
|
const double & | get_deactivation_time_constant () const |
| Get the value of the deactivation_time_constant property. More...
|
|
double & | upd_deactivation_time_constant () |
| Get a writable reference to the deactivation_time_constant property. More...
|
|
void | set_deactivation_time_constant (const double &value) |
| Set the value of the deactivation_time_constant property. More...
|
|
const double & | get_minimum_activation () const |
| Get the value of the minimum_activation property. More...
|
|
double & | upd_minimum_activation () |
| Get a writable reference to the minimum_activation property. More...
|
|
void | set_minimum_activation (const double &value) |
| Set the value of the minimum_activation property. More...
|
|
const double & | get_maximum_pennation_angle () const |
| Get the value of the maximum_pennation_angle property. More...
|
|
double & | upd_maximum_pennation_angle () |
| Get a writable reference to the maximum_pennation_angle property. More...
|
|
void | set_maximum_pennation_angle (const double &value) |
| Set the value of the maximum_pennation_angle property. More...
|
|
const ActiveForceLengthCurve & | get_ActiveForceLengthCurve () const |
| Get the value of the ActiveForceLengthCurve property. More...
|
|
ActiveForceLengthCurve & | upd_ActiveForceLengthCurve () |
| Get a writable reference to the ActiveForceLengthCurve property. More...
|
|
void | set_ActiveForceLengthCurve (const ActiveForceLengthCurve &value) |
| Set the value of the ActiveForceLengthCurve property. More...
|
|
const ForceVelocityCurve & | get_ForceVelocityCurve () const |
| Get the value of the ForceVelocityCurve property. More...
|
|
ForceVelocityCurve & | upd_ForceVelocityCurve () |
| Get a writable reference to the ForceVelocityCurve property. More...
|
|
void | set_ForceVelocityCurve (const ForceVelocityCurve &value) |
| Set the value of the ForceVelocityCurve property. More...
|
|
const FiberForceLengthCurve & | get_FiberForceLengthCurve () const |
| Get the value of the FiberForceLengthCurve property. More...
|
|
FiberForceLengthCurve & | upd_FiberForceLengthCurve () |
| Get a writable reference to the FiberForceLengthCurve property. More...
|
|
void | set_FiberForceLengthCurve (const FiberForceLengthCurve &value) |
| Set the value of the FiberForceLengthCurve property. More...
|
|
const TendonForceLengthCurve & | get_TendonForceLengthCurve () const |
| Get the value of the TendonForceLengthCurve property. More...
|
|
TendonForceLengthCurve & | upd_TendonForceLengthCurve () |
| Get a writable reference to the TendonForceLengthCurve property. More...
|
|
void | set_TendonForceLengthCurve (const TendonForceLengthCurve &value) |
| Set the value of the TendonForceLengthCurve property. More...
|
|
const double & | get_max_isometric_force () const |
| Get the value of the max_isometric_force property. More...
|
|
double & | upd_max_isometric_force () |
| Get a writable reference to the max_isometric_force property. More...
|
|
void | set_max_isometric_force (const double &value) |
| Set the value of the max_isometric_force property. More...
|
|
const double & | get_optimal_fiber_length () const |
| Get the value of the optimal_fiber_length property. More...
|
|
double & | upd_optimal_fiber_length () |
| Get a writable reference to the optimal_fiber_length property. More...
|
|
void | set_optimal_fiber_length (const double &value) |
| Set the value of the optimal_fiber_length property. More...
|
|
const double & | get_tendon_slack_length () const |
| Get the value of the tendon_slack_length property. More...
|
|
double & | upd_tendon_slack_length () |
| Get a writable reference to the tendon_slack_length property. More...
|
|
void | set_tendon_slack_length (const double &value) |
| Set the value of the tendon_slack_length property. More...
|
|
const double & | get_pennation_angle_at_optimal () const |
| Get the value of the pennation_angle_at_optimal property. More...
|
|
double & | upd_pennation_angle_at_optimal () |
| Get a writable reference to the pennation_angle_at_optimal property. More...
|
|
void | set_pennation_angle_at_optimal (const double &value) |
| Set the value of the pennation_angle_at_optimal property. More...
|
|
const double & | get_max_contraction_velocity () const |
| Get the value of the max_contraction_velocity property. More...
|
|
double & | upd_max_contraction_velocity () |
| Get a writable reference to the max_contraction_velocity property. More...
|
|
void | set_max_contraction_velocity (const double &value) |
| Set the value of the max_contraction_velocity property. More...
|
|
const bool & | get_ignore_tendon_compliance () const |
| Get the value of the ignore_tendon_compliance property. More...
|
|
bool & | upd_ignore_tendon_compliance () |
| Get a writable reference to the ignore_tendon_compliance property. More...
|
|
void | set_ignore_tendon_compliance (const bool &value) |
| Set the value of the ignore_tendon_compliance property. More...
|
|
const bool & | get_ignore_activation_dynamics () const |
| Get the value of the ignore_activation_dynamics property. More...
|
|
bool & | upd_ignore_activation_dynamics () |
| Get a writable reference to the ignore_activation_dynamics property. More...
|
|
void | set_ignore_activation_dynamics (const bool &value) |
| Set the value of the ignore_activation_dynamics property. More...
|
|
| Muscle () |
| Default constructor. More...
|
|
double | getMaxIsometricForce () const |
| get/set the maximum isometric force (in N) that the fibers can generate More...
|
|
void | setMaxIsometricForce (double maxIsometricForce) |
|
double | getOptimalFiberLength () const |
| get/set the optimal length (in m) of the muscle fibers (lumped as a single fiber) More...
|
|
void | setOptimalFiberLength (double optimalFiberLength) |
|
double | getTendonSlackLength () const |
| get/set the resting (slack) length (in m) of the tendon that is in series with the muscle fiber More...
|
|
void | setTendonSlackLength (double tendonSlackLength) |
|
double | getPennationAngleAtOptimalFiberLength () const |
| get/set the angle (in radians) between fibers at their optimal fiber length and the tendon More...
|
|
void | setPennationAngleAtOptimalFiberLength (double pennationAngle) |
|
double | getMaxContractionVelocity () const |
| get/set the maximum contraction velocity of the fibers, in optimal fiber-lengths per second More...
|
|
void | setMaxContractionVelocity (double maxContractionVelocity) |
|
bool | getIgnoreTendonCompliance (const SimTK::State &s) const |
| Get/set Modeling (runtime) option to ignore tendon compliance when computing muscle dynamics. More...
|
|
void | setIgnoreTendonCompliance (SimTK::State &s, bool ignore) const |
|
bool | getIgnoreActivationDynamics (const SimTK::State &s) const |
| Get/set Modeling (runtime) option to ignore activation dynamics when computing muscle dynamics. More...
|
|
void | setIgnoreActivationDynamics (SimTK::State &s, bool ignore) const |
|
virtual double | getActivation (const SimTK::State &s) const |
| get the activation level of the muscle, which modulates the active force of the muscle and has a normalized (0 to 1) value Note: method remains virtual to permit override by deprecated muscles. More...
|
|
double | getFiberLength (const SimTK::State &s) const |
| get the current working fiber length (m) for the muscle More...
|
|
double | getPennationAngle (const SimTK::State &s) const |
| get the current pennation angle (radians) between the fiber and tendon at the current fiber length More...
|
|
double | getCosPennationAngle (const SimTK::State &s) const |
| get the cosine of the current pennation angle (radians) between the fiber and tendon at the current fiber length More...
|
|
double | getTendonLength (const SimTK::State &s) const |
| get the current tendon length (m) given the current joint angles and fiber length More...
|
|
double | getNormalizedFiberLength (const SimTK::State &s) const |
| get the current normalized fiber length (fiber_length/optimal_fiber_length) More...
|
|
double | getFiberLengthAlongTendon (const SimTK::State &s) const |
| get the current fiber length (m) projected (*cos(pennationAngle)) onto the tendon direction More...
|
|
double | getTendonStrain (const SimTK::State &s) const |
| get the current tendon strain (delta_l/tendon_slack_length is dimensionless) More...
|
|
double | getFiberPotentialEnergy (const SimTK::State &s) const |
| the potential energy (J) stored in the fiber due to its parallel elastic element More...
|
|
double | getTendonPotentialEnergy (const SimTK::State &s) const |
| the potential energy (J) stored in the tendon More...
|
|
double | getMusclePotentialEnergy (const SimTK::State &s) const |
| the total potential energy (J) stored in the muscle More...
|
|
double | getPassiveForceMultiplier (const SimTK::State &s) const |
| get the passive fiber (parallel elastic element) force multiplier More...
|
|
double | getActiveForceLengthMultiplier (const SimTK::State &s) const |
| get the active fiber (contractile element) force multiplier due to current fiber length More...
|
|
double | getFiberVelocity (const SimTK::State &s) const |
| get current fiber velocity (m/s) positive is lengthening More...
|
|
double | getNormalizedFiberVelocity (const SimTK::State &s) const |
| get normalized fiber velocity. More...
|
|
double | getFiberVelocityAlongTendon (const SimTK::State &s) const |
| get the current fiber velocity (m/s) projected onto the tendon direction More...
|
|
double | getPennationAngularVelocity (const SimTK::State &s) const |
| get pennation angular velocity (radians/s) More...
|
|
double | getTendonVelocity (const SimTK::State &s) const |
| get the tendon velocity (m/s) positive is lengthening More...
|
|
double | getForceVelocityMultiplier (const SimTK::State &s) const |
| get the dimensionless multiplier resulting from the fiber's force-velocity curve More...
|
|
double | getFiberForce (const SimTK::State &s) const |
| get the current fiber force (N) applied to the tendon More...
|
|
double | getFiberForceAlongTendon (const SimTK::State &s) const |
| get the force of the fiber (N/m) along the direction of the tendon More...
|
|
double | getActiveFiberForce (const SimTK::State &s) const |
| get the current active fiber force (N) due to activation*force_length*force_velocity relationships More...
|
|
double | getPassiveFiberForce (const SimTK::State &s) const |
| get the total force applied by all passive elements in the fiber (N) More...
|
|
double | getActiveFiberForceAlongTendon (const SimTK::State &s) const |
| get the current active fiber force (N) projected onto the tendon direction More...
|
|
double | getPassiveFiberForceAlongTendon (const SimTK::State &s) const |
| get the total force applied by all passive elements in the fiber (N) projected onto the tendon direction More...
|
|
double | getTendonForce (const SimTK::State &s) const |
| get the current tendon force (N) applied to bones More...
|
|
double | getFiberStiffness (const SimTK::State &s) const |
| get the current fiber stiffness (N/m) defined as the partial derivative of fiber force with respect to fiber length More...
|
|
double | getFiberStiffnessAlongTendon (const SimTK::State &s) const |
| get the stiffness of the fiber (N/m) along the direction of the tendon, that is the partial derivative of the fiber force along the tendon with respect to small changes in fiber length along the tendon More...
|
|
double | getTendonStiffness (const SimTK::State &s) const |
| get the current tendon stiffness (N/m) defined as the partial derivative of tendon force with respect to tendon length More...
|
|
double | getMuscleStiffness (const SimTK::State &s) const |
| get the current muscle stiffness (N/m) defined as the partial derivative of muscle force with respect to muscle length More...
|
|
double | getFiberActivePower (const SimTK::State &s) const |
| get the current active fiber power (W) More...
|
|
double | getFiberPassivePower (const SimTK::State &s) const |
| get the current passive fiber power (W) More...
|
|
double | getTendonPower (const SimTK::State &s) const |
| get the current tendon power (W) More...
|
|
double | getMusclePower (const SimTK::State &s) const |
| get the current muscle power (W) More...
|
|
double | getStress (const SimTK::State &s) const override |
| get the stress in the muscle (part of the Actuator interface as well) More...
|
|
void | setExcitation (SimTK::State &s, double excitation) const |
| set the excitation (control) for this muscle. More...
|
|
double | getExcitation (const SimTK::State &s) const |
|
void | computeEquilibrium (SimTK::State &s) const override final |
| Find and set the equilibrium state of the muscle (if any) More...
|
|
| PathActuator () |
|
GeometryPath & | updGeometryPath () |
|
const GeometryPath & | getGeometryPath () const |
|
bool | hasGeometryPath () const override |
| Return a flag indicating whether the Force is applied along a Path. More...
|
|
void | setOptimalForce (double aOptimalForce) |
|
double | getOptimalForce () const override |
|
virtual double | getLength (const SimTK::State &s) const |
|
virtual double | getLengtheningSpeed (const SimTK::State &s) const |
|
double | getPower (const SimTK::State &s) const override |
|
double | getStress (const SimTK::State &s) const override |
|
void | addNewPathPoint (const std::string &proposedName, const PhysicalFrame &aBody, const SimTK::Vec3 &aPositionOnBody) |
| Note that this function does not maintain the State and so should be used only before a valid State is created. More...
|
|
double | computeActuation (const SimTK::State &s) const override |
|
virtual double | computeMomentArm (const SimTK::State &s, Coordinate &aCoord) const |
|
const GeometryPath & | get_GeometryPath () const |
| Get the value of the GeometryPath property. More...
|
|
GeometryPath & | upd_GeometryPath () |
| Get a writable reference to the GeometryPath property. More...
|
|
void | set_GeometryPath (const GeometryPath &value) |
| Set the value of the GeometryPath property. More...
|
|
const double & | get_optimal_force () const |
| Get the value of the optimal_force property. More...
|
|
double & | upd_optimal_force () |
| Get a writable reference to the optimal_force property. More...
|
|
void | set_optimal_force (const double &value) |
| Set the value of the optimal_force property. More...
|
|
PathActuator * | clone () const override |
| Create a new heap-allocated copy of the concrete object to which this Object refers. More...
|
|
const std::string & | getConcreteClassName () const override |
| Returns the class name of the concrete Object-derived class of the actual object referenced by this Object, as a string. More...
|
|
| ScalarActuator () |
|
virtual double | getControl (const SimTK::State &s) const |
| Convenience method to get control given scalar (double) valued control. More...
|
|
int | numControls () const override |
|
virtual void | setActuation (const SimTK::State &s, double aActuation) const |
|
virtual double | getActuation (const SimTK::State &s) const |
|
virtual void | setSpeed (const SimTK::State &s, double aspeed) const |
|
virtual double | getSpeed (const SimTK::State &s) const |
|
void | setMinControl (const double &aMinControl) |
| Methods to manage the bounds on ScalarActuator's control. More...
|
|
double | getMinControl () const |
|
void | setMaxControl (const double &aMaxControl) |
|
double | getMaxControl () const |
|
void | overrideActuation (SimTK::State &s, bool flag) const |
| Enable/disable a ScalarActuator's override actuation. More...
|
|
bool | isActuationOverridden (const SimTK::State &s) const |
| return ScalarActuator's override status More...
|
|
void | setOverrideActuation (SimTK::State &s, double value) const |
| set the actuation value used when the override is true More...
|
|
double | getOverrideActuation (const SimTK::State &s) const |
| return override actuation More...
|
|
const double & | get_min_control () const |
| Get the value of the min_control property. More...
|
|
double & | upd_min_control () |
| Get a writable reference to the min_control property. More...
|
|
void | set_min_control (const double &value) |
| Set the value of the min_control property. More...
|
|
const double & | get_max_control () const |
| Get the value of the max_control property. More...
|
|
double & | upd_max_control () |
| Get a writable reference to the max_control property. More...
|
|
void | set_max_control (const double &value) |
| Set the value of the max_control property. More...
|
|
| Actuator () |
|
virtual const SimTK::Vector | getDefaultControls () |
| Actuator default controls are zero. More...
|
|
virtual const SimTK::VectorView_< double > | getControls (const SimTK::State &s) const |
|
virtual void | getControls (const SimTK::Vector &modelControls, SimTK::Vector &actuatorControls) const |
| Convenience methods for getting, setting and adding to actuator controls from/into the model controls. More...
|
|
virtual void | setControls (const SimTK::Vector &actuatorControls, SimTK::Vector &modelControls) const |
| set actuator controls subvector into the right slot in the system-wide model controls More...
|
|
virtual void | addInControls (const SimTK::Vector &actuatorControls, SimTK::Vector &modelControls) const |
| add actuator controls to the values already occupying the slot in the system-wide model controls More...
|
|
virtual bool | shouldBeParallelized () const |
| Tell SimBody to parallelize this force. More...
|
|
bool | appliesForce (const SimTK::State &s) const |
| Return if the Force is applied (or enabled) or not. More...
|
|
void | setAppliesForce (SimTK::State &s, bool applyForce) const |
| Set whether or not the Force is applied. More...
|
|
const bool & | get_appliesForce () const |
| Get the value of the appliesForce property. More...
|
|
bool & | upd_appliesForce () |
| Get a writable reference to the appliesForce property. More...
|
|
void | set_appliesForce (const bool &value) |
| Set the value of the appliesForce property. More...
|
|
| ModelComponent () |
| Default constructor. More...
|
|
| ModelComponent (const std::string &aFileName, bool aUpdateFromXMLNode=true) |
| Construct ModelComponent from an XML file. More...
|
|
| ModelComponent (SimTK::Xml::Element &aNode) |
| Construct ModelComponent from a specific node in an XML document. More...
|
|
virtual | ~ModelComponent () |
| Destructor is virtual to allow concrete model component cleanup. More...
|
|
void | connectToModel (Model &model) |
| Connect this ModelComponent to its aggregate- a Model. More...
|
|
const Model & | getModel () const |
| Get a const reference to the Model this component is part of. More...
|
|
Model & | updModel () |
| Get a modifiable reference to the Model this component is part of. More...
|
|
bool | hasModel () const |
| Does this ModelComponent have a Model associated with it? More...
|
|
void | preScale (const SimTK::State &s, const ScaleSet &scaleSet) |
| Perform any computations that must occur before ModelComponent::scale() is invoked on all ModelComponents in the Model. More...
|
|
void | scale (const SimTK::State &s, const ScaleSet &scaleSet) |
| Scale the ModelComponent. More...
|
|
void | postScale (const SimTK::State &s, const ScaleSet &scaleSet) |
| Perform any computations that must occur after ModelComponent::scale() has been invoked on all ModelComponents in the Model. More...
|
|
| Component () |
| Default constructor. More...
|
|
| Component (const std::string &aFileName, bool aUpdateFromXMLNode=true) |
| Construct Component from an XML file. More...
|
|
| Component (SimTK::Xml::Element &aNode) |
| Construct Component from a specific node in an XML document. More...
|
|
| Component (const Component &)=default |
| Use default copy constructor and assignment operator. More...
|
|
Component & | operator= (const Component &)=default |
|
virtual | ~Component ()=default |
| Destructor is virtual to allow concrete Component to cleanup. More...
|
|
virtual void | generateDecorations (bool fixed, const ModelDisplayHints &hints, const SimTK::State &state, SimTK::Array_< SimTK::DecorativeGeometry > &appendToThis) const |
| Optional method for generating arbitrary display geometry that reflects this Component at the specified state. More...
|
|
const SimTK::MultibodySystem & | getSystem () const |
| Get the underlying MultibodySystem that this component is connected to. More...
|
|
bool | hasSystem () const |
| Check if this component has an underlying MultibodySystem. More...
|
|
bool | isComponentInOwnershipTree (const Component *component) const |
| Does the provided component already exist anywhere in the ownership tree (not just subcomponents of this component)? More...
|
|
void | addComponent (Component *subcomponent) |
|
Add a Component (as a subcomponent) of this component. More...
|
|
template<typename T = Component> |
ComponentList< const T > | getComponentList () const |
| Get an iterator through the underlying subcomponents that this component is composed of. More...
|
|
template<typename T = Component> |
ComponentList< T > | updComponentList () |
| Similar to getComponentList(), except the resulting list allows one to modify the components. More...
|
|
template<typename T = Component> |
unsigned | countNumComponents () const |
| Uses getComponentList<T>() to count the number of underlying subcomponents of the specified type. More...
|
|
std::string | getAbsolutePathString () const |
| Get the complete (absolute) pathname for this Component to its ancestral Component, which is the root of the tree to which this Component belongs. More...
|
|
ComponentPath | getAbsolutePath () const |
| Return a ComponentPath of the absolute path of this Component. More...
|
|
std::string | getRelativePathString (const Component &wrt) const |
| Get the relative path of this Component with respect to another Component, as a string. More...
|
|
ComponentPath | getRelativePath (const Component &wrt) const |
| Get the relative path of this Component with respect to another Component. More...
|
|
bool | hasComponent (const std::string &pathname) const |
| Query if there is a component (of any type) at the specified path name. More...
|
|
template<class C = Component> |
bool | hasComponent (const std::string &pathname) const |
| Query if there is a component of a given type at the specified path name. More...
|
|
template<class C = Component> |
const C & | getComponent (const std::string &pathname) const |
| Get a unique subcomponent of this Component by its path name and type 'C'. More...
|
|
template<class C = Component> |
const C & | getComponent (const ComponentPath &pathname) const |
|
const Component & | getComponent (const std::string &pathname) const |
| Similar to the templatized getComponent(), except this returns the component as the generic Component type. More...
|
|
template<class C = Component> |
C & | updComponent (const std::string &name) |
| Get a writable reference to a subcomponent. More...
|
|
template<class C = Component> |
C & | updComponent (const ComponentPath &name) |
|
Component & | updComponent (const std::string &pathname) |
| Similar to the templatized updComponent(), except this returns the component as the generic Component type. More...
|
|
unsigned | printComponentsMatching (const std::string &substring) const |
| Print a list to the console of all components whose absolute path name contains the given string. More...
|
|
int | getNumStateVariables () const |
| Get the number of "continuous" state variables maintained by the Component and its subcomponents. More...
|
|
Array< std::string > | getStateVariableNames () const |
| Get the names of "continuous" state variables maintained by the Component and its subcomponents. More...
|
|
template<class C = Component> |
const C * | findComponent (const ComponentPath &pathToFind) const |
| Find a Component to which this Component is an ancestor—in other words, a Component that is directly owned by this Component or is owned by one of its sub-components, sub-sub-components, etc. More...
|
|
template<class C = Component> |
const C * | findComponent (const std::string &pathToFind) const |
| Same as findComponent(const ComponentPath&), but accepting a string (a path or just a name) as input. More...
|
|
const StateVariable * | traverseToStateVariable (const std::string &pathName) const |
| Get a StateVariable anywhere in the Component tree, given a StateVariable path. More...
|
|
void | finalizeFromProperties () |
| Define a Component's internal data members and structure according to its properties. More...
|
|
void | finalizeConnections (Component &root) |
| Satisfy the Component's connections specified by its Sockets and Inputs. More...
|
|
void | clearConnections () |
| Disconnect/clear this Component from its aggregate component. More...
|
|
void | addToSystem (SimTK::MultibodySystem &system) const |
| Have the Component add itself to the underlying computational System. More...
|
|
void | initStateFromProperties (SimTK::State &state) const |
| Initialize Component's state variable values from its properties. More...
|
|
void | setPropertiesFromState (const SimTK::State &state) |
| Set Component's properties given a state. More...
|
|
int | getNumSockets () const |
| Get the number of Sockets in this Component. More...
|
|
std::vector< std::string > | getSocketNames () |
| Collect and return the names of the sockets in this component. More...
|
|
template<typename T > |
const T & | getConnectee (const std::string &name) const |
| Get the "connectee" object that the Component's Socket is bound to. More...
|
|
const Object & | getConnectee (const std::string &name) const |
| Get the connectee as an Object. More...
|
|
const AbstractSocket & | getSocket (const std::string &name) const |
| Get an AbstractSocket for the given socket name. More...
|
|
AbstractSocket & | updSocket (const std::string &name) |
| Get a writable reference to the AbstractSocket for the given socket name. More...
|
|
template<typename T > |
const Socket< T > & | getSocket (const std::string &name) const |
| Get a const reference to the concrete Socket provided by this Component by name. More...
|
|
template<typename T > |
Socket< T > & | updSocket (const std::string &name) |
| Get a writable reference to the concrete Socket provided by this Component by name. More...
|
|
void | printSubcomponentInfo () const |
| List all subcomponents by name and recurse into these components to list their subcomponents, and so on. More...
|
|
void | printSocketInfo () const |
| List all the Sockets of this component and whether or not they are connected. More...
|
|
void | printInputInfo () const |
| List all the inputs of this component and whether or not they are connected. More...
|
|
template<typename C > |
void | printSubcomponentInfo () const |
|
void | printOutputInfo (const bool includeDescendants=true) const |
| Print outputs of this component and optionally, those of all subcomponents. More...
|
|
int | getNumInputs () const |
| Access the number of Inputs that this component has. More...
|
|
int | getNumOutputs () const |
| Access the number of Outputs that this component has. More...
|
|
std::vector< std::string > | getInputNames () const |
| Collect and return the names of Inputs in this component as an std::vector. More...
|
|
std::vector< std::string > | getOutputNames () const |
| Collect and return the names of Outputs in this component as an std::vector. More...
|
|
const AbstractInput & | getInput (const std::string &name) const |
| Get an Input provided by this Component by name. More...
|
|
AbstractInput & | updInput (const std::string &name) |
| Get a writable reference to an Input provided by this Component by name. More...
|
|
template<typename T > |
const Input< T > & | getInput (const std::string &name) const |
| Get a concrete Input that you can directly ask for its values. More...
|
|
const AbstractOutput & | getOutput (const std::string &name) const |
| Get the Output provided by this Component by name. More...
|
|
AbstractOutput & | updOutput (const std::string &name) |
| Get a writable reference to an Output provided by this Component by name. More...
|
|
SimTK::IteratorRange< OutputConstIterator > | getOutputs () const |
| Iterate through all Outputs of this component. More...
|
|
int | getModelingOption (const SimTK::State &state, const std::string &name) const |
| Get a ModelingOption flag for this Component by name. More...
|
|
void | setModelingOption (SimTK::State &state, const std::string &name, int flag) const |
| Set the value of a ModelingOption flag for this Component. More...
|
|
template<typename T > |
const T & | getInputValue (const SimTK::State &state, const std::string &name) const |
| Get the Input value that this component is dependent on. More...
|
|
template<typename T > |
const T & | getOutputValue (const SimTK::State &state, const std::string &name) const |
| Get the Output value provided by this Component by name. More...
|
|
double | getStateVariableValue (const SimTK::State &state, const std::string &name) const |
| Get the value of a state variable allocated by this Component. More...
|
|
void | setStateVariableValue (SimTK::State &state, const std::string &name, double value) const |
| Set the value of a state variable allocated by this Component by name. More...
|
|
SimTK::Vector | getStateVariableValues (const SimTK::State &state) const |
| Get all values of the state variables allocated by this Component. More...
|
|
void | setStateVariableValues (SimTK::State &state, const SimTK::Vector &values) const |
| Set all values of the state variables allocated by this Component. More...
|
|
double | getStateVariableDerivativeValue (const SimTK::State &state, const std::string &name) const |
| Get the value of a state variable derivative computed by this Component. More...
|
|
double | getDiscreteVariableValue (const SimTK::State &state, const std::string &name) const |
| Get the value of a discrete variable allocated by this Component by name. More...
|
|
void | setDiscreteVariableValue (SimTK::State &state, const std::string &name, double value) const |
| Set the value of a discrete variable allocated by this Component by name. More...
|
|
template<class T > |
SimTK::CacheEntryIndex | getCacheVariableIndex (const CacheVariable< T > &cv) const |
| Get the index of a Component's cache variable in the Subsystem for allocations. More...
|
|
SimTK::CacheEntryIndex | getCacheVariableIndex (const std::string &name) const |
| Get the index of a Component's cache variable in the Subsystem for allocations. More...
|
|
template<class T > |
const T & | getCacheVariableValue (const SimTK::State &state, const std::string &name) const |
| Get the value of a cache variable allocated by this Component by name. More...
|
|
template<class T > |
const T & | getCacheVariableValue (const SimTK::State &state, const CacheVariable< T > &cv) const |
| Get the value of a cache variable allocated by this Component. More...
|
|
template<typename T > |
void | setCacheVariableValue (const SimTK::State &state, const std::string &k, T value) const |
| Set the value of a cache variable, identified by name , to a new value and mark the cache variable as valid. More...
|
|
template<typename T > |
void | setCacheVariableValue (const SimTK::State &state, const CacheVariable< T > &cv, T value) const |
| Set the value of a cache variable to a new value and mark the cache variable as valid. More...
|
|
template<typename T > |
T & | updCacheVariableValue (const SimTK::State &state, const std::string &name) const |
| Returns a mutable reference to the value of a cache variable identified by name . More...
|
|
template<typename T > |
T & | updCacheVariableValue (const SimTK::State &state, const CacheVariable< T > &cv) const |
| Returns a mutable reference to the value of a cache variable. More...
|
|
bool | isCacheVariableValid (const SimTK::State &state, const std::string &name) const |
| Returns true if the cache variable, identified by name , is valid. More...
|
|
template<class T > |
bool | isCacheVariableValid (const SimTK::State &state, const CacheVariable< T > &cv) const |
| Returns true if the cache variable is valid. More...
|
|
void | markCacheVariableValid (const SimTK::State &state, const std::string &name) const |
| Marks the value of a cache variable, identified by name , as valid. More...
|
|
template<typename T > |
void | markCacheVariableValid (const SimTK::State &state, const CacheVariable< T > &cv) const |
| Marks the value of a cache variable as valid. More...
|
|
void | markCacheVariableInvalid (const SimTK::State &state, const std::string &name) const |
| Marks the value of a cache variable, identified by name , as invalid. More...
|
|
template<class T > |
void | markCacheVariableInvalid (const SimTK::State &state, const CacheVariable< T > &cv) const |
| Marks the value of a cache variable as invalid. More...
|
|
const Component & | getOwner () const |
| Access the owner of this Component. More...
|
|
bool | hasOwner () const |
| (For advanced users) Check if this Component has an owner. More...
|
|
const Component & | getRoot () const |
| Obtain the root Component, which is this component if it is orphaned. More...
|
|
virtual | ~Object () |
| Virtual destructor for cleanup. More...
|
|
bool | isEqualTo (const Object &aObject) const |
| Equality operator wrapper for use from languages not supporting operator overloading. More...
|
|
Object & | operator= (const Object &aObject) |
| Copy assignment copies he base class fields, including the properties. More...
|
|
virtual bool | operator== (const Object &aObject) const |
| Determine if two objects are equal. More...
|
|
virtual bool | operator< (const Object &aObject) const |
| Provide an ordering for objects so they can be put in sorted containers. More...
|
|
void | setName (const std::string &name) |
| Set the name of the Object. More...
|
|
const std::string & | getName () const |
| Get the name of this Object. More...
|
|
void | setDescription (const std::string &description) |
| Set description, a one-liner summary. More...
|
|
const std::string & | getDescription () const |
| Get description, a one-liner summary. More...
|
|
const std::string & | getAuthors () const |
| Get Authors of this Object. More...
|
|
void | setAuthors (const std::string &authors) |
| Set Authors of this object. More...
|
|
const std::string & | getReferences () const |
| Get references or publications to cite if using this object. More...
|
|
void | setReferences (const std::string &references) |
| Set references or publications to cite if using this object. More...
|
|
int | getNumProperties () const |
| Determine how many properties are stored with this Object. More...
|
|
const AbstractProperty & | getPropertyByIndex (int propertyIndex) const |
| Get a const reference to a property by its index number, returned as an AbstractProperty. More...
|
|
AbstractProperty & | updPropertyByIndex (int propertyIndex) |
| Get a writable reference to a property by its index number, returned as an AbstractProperty. More...
|
|
bool | hasProperty (const std::string &name) const |
| Return true if this Object has a property of any type with the given name, which must not be empty. More...
|
|
const AbstractProperty & | getPropertyByName (const std::string &name) const |
| Get a const reference to a property by its name, returned as an AbstractProperty. More...
|
|
AbstractProperty & | updPropertyByName (const std::string &name) |
| Get a writable reference to a property by its name, returned as an AbstractProperty. More...
|
|
template<class T > |
bool | hasProperty () const |
| Return true if this Object contains an unnamed, one-object property that contains objects of the given template type T. More...
|
|
template<class T > |
const Property< T > & | getProperty (const PropertyIndex &index) const |
| Get property of known type Property<T> as a const reference; the property must be present and have the right type. More...
|
|
template<class T > |
Property< T > & | updProperty (const PropertyIndex &index) |
| Get property of known type Property<T> as a writable reference; the property must be present and have the right type. More...
|
|
bool | isObjectUpToDateWithProperties () const |
| Returns true if no property's value has changed since the last time setObjectIsUpToDateWithProperties() was called. More...
|
|
void | readObjectFromXMLNodeOrFile (SimTK::Xml::Element &objectElement, int versionNumber) |
| We're given an XML element from which we are to populate this Object. More...
|
|
void | updateXMLNode (SimTK::Xml::Element &parent, const AbstractProperty *prop=nullptr) const |
| Serialize this object into the XML node that represents it. More...
|
|
bool | getInlined () const |
| Inlined means an in-memory Object that is not associated with an XMLDocument. More...
|
|
void | setInlined (bool aInlined, const std::string &aFileName="") |
| Mark this as inlined or not and optionally provide a file name to associate with the new XMLDocument for the non-inline case. More...
|
|
std::string | getDocumentFileName () const |
| If there is a document associated with this object then return the file name maintained by the document. More...
|
|
int | getDocumentFileVersion () const |
| If there is a document associated with this object then return its version number. More...
|
|
void | setAllPropertiesUseDefault (bool aUseDefault) |
|
bool | print (const std::string &fileName) const |
| Write this Object into an XML file of the given name; conventionally the suffix to use is ".osim". More...
|
|
std::string | dump () const |
| dump the XML representation of this Object into an std::string and return it. More...
|
|
virtual bool | isA (const char *type) const |
| The default implementation returns true only if the supplied string is "Object"; each Object-derived class overrides this to match its own class name. More...
|
|
const std::string & | toString () const |
| Wrapper to be used on Java side to display objects in tree; this returns just the object's name. More...
|
|
PropertySet & | getPropertySet () |
| OBSOLETE: Get a reference to the PropertySet maintained by the Object. More...
|
|
const PropertySet & | getPropertySet () const |
|
This class implements a configurable equilibrium muscle model, as described in Millard et al. (2013).
An equilibrium model assumes that the forces generated by the fiber and tendon are equal:
\[ f_{ISO}\Big(\mathbf{a}(t) \mathbf{f}_L(\hat{l}_{CE}) \mathbf{f}_V(\hat{v}_{CE}) + \mathbf{f}_{PE}(\hat{l}_{CE}) + \beta \hat{v}_{CE}\Big) \cos \phi - f_{ISO}\mathbf{f}_{SE}(\hat{l}_{T}) = 0 \]
This model can be simulated in several configurations by adjusting three flags:
- ignore_tendon_compliance: set to true to make the tendon rigid. This assumption is usually reasonable for short tendons, and can result in a performance improvement by eliminating high-frequency dynamics and removing the fiber length from the state vector.
- ignore_activation_dynamics: set to true to use the excitation input as the activation signal. This results in faster simulations by reducing the size of the state vector.
- fiber_damping: set to a value greater than 0.001 to include fiber damping in the model. The addition of damping reduces simulation time while allowing the muscle model to be more physiological (it can have an activation of zero, its active-force-length curve can go to zero, and its force-velocity curve can be asymptotic).
Elastic Tendon, No Fiber Damping
The most typical configuration used in the literature is to simulate a muscle with an elastic tendon, full fiber dynamics, and activation dynamics. The resulting formulation suffers from three singularities: \(\mathbf{a}(t) \rightarrow 0\), \(\phi \rightarrow 90^\circ\), and \( \mathbf{f}_L(\hat{l}_{CE}) \rightarrow 0 \). These situations are all handled in this model to ensure that it does not produce singularities and does not result in intolerably long simulation times.
Numerical singularities arise from the manner in which the equilibrium equation is rearranged to yield an ordinary differential equation (ODE). The above equation is rearranged to isolate \( \mathbf{f}_V(\hat{v}_{CE}) \). We then invert to solve for \( \hat{v}_{CE} \), which is then numerically integrated during a simulation:
\[ \hat{v}_{CE} = \mathbf{f}_V ^{-1} \Big( \frac{ ( \mathbf{f}_{SE}(\hat{l}_{T}) ) / \cos \phi - \mathbf{f}_{PE}(\hat{l}_{CE}) } { \mathbf{a}(t) \mathbf{f}_L(\hat{l}_{CE})} \Big) \]
The above equation becomes numerically stiff when terms in the denominator approach zero (when \(\mathbf{a}(t) \rightarrow 0\), \(\phi \rightarrow 90^\circ\), or \( \mathbf{f}_L(\hat{l}_{CE}) \rightarrow 0 \)) or, additionally, when the slope of \(\mathbf{f}_V ^{-1}\) is steep (which occurs at fiber velocities close to the maximum concentric and maximum eccentric fiber velocities).
Singularities can be managed by ensuring that the muscle model is always activated ( \(\mathbf{a}(t) > 0\)), the fiber will stop contracting when a pennation angle of 90 degrees is approached ( \(\phi < 90^\circ\)), and the fiber will also stop contracting as its length approaches a lower bound ( \( \hat{l}_{CE} > lowerbound\)), which is typically around half the fiber's resting length (to ensure \( \mathbf{f}_L(\hat{l}_{CE}) > 0 \)). The fiber is prevented from reaching unphysiological lengths or its maximum pennation angle using a unilateral constraint. Additionally, the force-velocity curve is modified so that it is invertible.
When an elastic tendon without fiber damping is selected, the minimum active-force-length value is set to 0.1, the minimum permissible activation is set to 0.01, and the maximum permissible pennation angle is set to acos(0.1) or 84.3 degrees. This is done as a convenience for the user to prevent the model from taking an unreasonable amount of time to simulate.
(Rigid Tendon) or (Elastic Tendon with Fiber Damping)
Neither of these formulations has any singularities. The lower bound of the active-force-length curve can be zero (min( \( \mathbf{f}_L(\hat{l}_{CE})) = 0 \)), activation can be zero (i.e., the muscle can be turned off completely), and the force-velocity curve need not be invertible.
The rigid tendon formulation removes the singularities by ignoring the elasticity of the tendon. This assumption is reasonable for many muscles, but it is up to the user to determine whether this assumption is valid.
The formulation that uses an elastic tendon with fiber damping removes singularities by solving the equilibrium equation with Newton's method. This is possible because the partial derivative of the equilibrium equation with respect to fiber velocity is always positive if \( \beta > 0\) and, thus, Newton's method can find a solution to the equilibrium equation.
When either of these singularity-free formulations is selected, the minimum active-force-length value and the minimum permissible activation are set to zero. This is done as a convenience for the user, as these changes make the results of the model more realistic yet incur no performance penalty. The maximum pennation angle is left as acos(0.1) or 84.3 degrees, as allowing higher pennation angles results in an increasingly stiff fiber velocity state as pennation angle increases.
Usage
This object should be updated through the set methods provided.
Example
double maxIsometricForce = 5000;
double optimalFiberLength = 0.025;
double tendonSlackLength = 0.25;
double pennationAngle = 0.5;
bool ignoreTendonCompliance = false;
bool ignoreActivationDynamics = false;
double dampingCoefficient = 0.001;
maxIsometricForce,
optimalFiberLength,
tendonSlackLength,
pennationAngle);
myMuscle.setMuscleConfiguration(ignoreTendonCompliance,
ignoreActivationDynamics,
dampingCoefficient);
Please refer to the doxygen for more information on the properties that are objects themselves (MuscleFixedWidthPennationModel, ActiveForceLengthCurve, FiberForceLengthCurve, TendonForceLengthCurve, and ForceVelocityInverseCurve).
Reference
Millard, M., Uchida, T., Seth, A., Delp, S.L. (2013) Flexing computational muscle: modeling and simulation of musculotendon dynamics. ASME Journal of Biomechanical Engineering 135(2):021005. http://dx.doi.org/10.1115/1.4023390.
- Author
- Matt Millard
-
Tom Uchida
-
Ajay Seth