CoordinateLimitForce - nonLinear region

Provide easy-to-use, extensible software for modeling, simulating, controlling, and analyzing the neuromusculoskeletal system.
POST REPLY
User avatar
Sina Porsa
Posts: 99
Joined: Thu Feb 03, 2011 7:21 pm

CoordinateLimitForce - nonLinear region

Post by Sina Porsa » Mon Feb 25, 2013 9:00 pm

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

User avatar
Michael Sherman
Posts: 807
Joined: Fri Apr 01, 2005 6:05 pm

Re: CoordinateLimitForce - nonLinear region

Post by Michael Sherman » Mon Feb 25, 2013 10:08 pm

Hi, Sina. The basic quintic step function is step(x)=10x^3-15x^4+6x^5, valid only for x between 0 and 1:
untitled.png
untitled.png (3.85 KiB) Viewed 534 times
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

User avatar
Sina Porsa
Posts: 99
Joined: Thu Feb 03, 2011 7:21 pm

Re: CoordinateLimitForce - nonLinear region

Post by Sina Porsa » Mon Feb 25, 2013 10:14 pm

Thank you Sherm

POST REPLY