Comments (5)
Hello Xue,
The convex optimization problem does have worse conditioning with a larger tau, making it more difficult to solve. You could try changing the fsolve settings to allow more iterations and have a tighter cutoff tolerance. Parameter stepping the tension is also a good option; when tau = [45, 0, 0, 0] is too hard on the solver, you could often have better luck solving tau = [15, 0, 0, 0], then using that solution as the initial guess for tau = [30, 0, 0, 0], then tau = [45, 0, 0, 0].
What tension is it failing to solve at? Unfortunately I do not have a MATLAB license to run 05_Tendon_Robot.m, but for the C++ version the solver fails somewhere between tau = [150, 0, 0, 0]N and tau = [200, 0, 0, 0]N. The scenario is already not very realistic at tau = [150, 0, 0, 0]N; the robot doubles back on itself.
Hopefully the MATLAB script has similar behavior with tighter solver tolerances!
from continuumrobotexamples.
Hello John,
Iterative solving the problem that has a large tau is a good idea. However, it does not work for a dynamic problem, for example, a large force is suddenly applied to the tender driven manipulator. It is very difficult to find a correct initial guess, especially, the time step is large. But to avoid "catastrophic cancellation" we cannot make the time step pretty small. What do you think of this problem? Do you think this is a drawback compared with the conventional way of solving PDEs, first spatial discretization then time discretization? I haven't tried but I feel the conventional way will be able to handle the case that a large force is suddenly applied, or the force has a large sudden change.
Best,
Jiefeng
from continuumrobotexamples.
Yes, using an implicit method for time integration is great for large time steps, but may cause problems at small time steps. Part of the issue is that shooting methods can have problems with poor conditioning. In an upcoming paper, we couldn't simulate concentric tube snapping with the shooting method, but we could use a finite difference scheme with the implicit time approach. The conventional way can work, but the time step needs to be small enough to meet the CFL condition, which can be quite limiting.
from continuumrobotexamples.
Thanks! Does 'poor conditioning' mean it is quite difficult to get the correct initial guess? The problem that I currently work on has a complicated shape (that has several possible static configurations that are close), a very small time step has to be chosen to make sure the initial guess is close enough. So I have to switch to the conventional method, but as you said CFL condition is limiting. I hope your new method can be helpful and look forward to reading your upcoming paper!
Best,
Jiefeng
from continuumrobotexamples.
Right, elements of the gradient descent Jacobian go to infinity if the guess is too far off, but with a close enough guess you can have a manageable slope.
from continuumrobotexamples.
Related Issues (14)
- Do you have a matlab code for Tendon Static Case HOT 3
- typos in the documentation CosseratRodIvpSolution.pdf HOT 2
- Question about TendonRobotStatics.m HOT 3
- about the fsolve() function in TendonRobotStatics.m HOT 1
- Questions regarding the formulation of the static equation in 01_Statics/02_CSG_Inverse_Kinematics/CSG_Inverse_Kinematics.pdf HOT 7
- Concern about the initial value of the internal forces in RodIVP.m HOT 2
- Concern about the initial value of the internal forces in RodIVP.m
- External force on an arbitry point. HOT 7
- Non-crucial Derivation Error HOT 1
- Rod with an Archimedean spiral reference shape HOT 4
- Frictional model for tendon driven continuum robot HOT 4
- New Blender versions break compatability HOT 1
- Build error of project โ01_Statics\06_Tendon_Robot_Variants\Cpp_3\Tendon_Displacement_Controlโ HOT 7
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
๐ Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. ๐๐๐
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google โค๏ธ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from continuumrobotexamples.