simtk.org

Hi everyone.
I am using CoordinateLimitForce to represent ligaments in my model. I know that this force has a constant stiffness except through the transition region which uses a non-linear function to change the stiffness smoothly from zero to the specified stiffness.
I was wondering if I can know the exact non-linear equation which defines the stiffness in this region?
It seems that it's a 5th degree polynomial. It would be great if I could know the constant coefficients.
Thanks,
Sina

Hi, Sina. The basic quintic step function is step(x)=10x^3-15x^4+6x^5, valid only for x between 0 and 1: However, there are many twists and turns in CoordinateLimitForce on the way to producing a force, such as how damping is calculated. These are internal details and are not part of the definition, so are subject to change in later versions. However, if you are interested you can look in the source code, which is available on the OpenSim 3.0 downloads page. The file you want to look in is OpenSim\Simulation\Model\CoordinateLimitForce.cpp. Forces are calculated in method calcLimitForce().

The step function itself is implemented using Simbody's SimTK::Function::Step() class which handles scaling and limiting the function and is documented here.

Regards,
Sherm

Thank you Sherm