Solving Linear Systems for Multibody Systems
4 views (last 30 days)
Tiago Carvalho on 21 Jul 2022
I'm currently coding a Multibody Foward Dynamcis Simulator in MATLAB, for my master thesis, and I am having issues with solving linear systems due to low RCOND values, ill-condition of the matrices.
This problem is expected in this type of problems since each column entry in a line corresponds to the elements of a joint equation for two bodies resulting in a highly sparse matrix and consequently ill-conditioned matrices.
What my algorithm does is it calculates the initial accelerations of bodies in order for them to be integrated through an ode solver for position and velocity. My problems arise in solving the linear system for the initial accelerations that will be integrated. The system I have to solve is the following:
Where: M is a mass matrix, Phiq is a Jacobian, alpha/omega/mu are scalars, Phi, gamma and upsilon are vectors.
For my algorithm I tried using pinv(A)*B or lsqr(A,B), but I am not able to get the correct results (only get constant values from initial time until the last integration). On the other hand I get good results by using mldivide or \, but i get the following message:
Warning: Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND = 1.394430e-19.
I wanted to know if there is anyway to avoid this issue, since I know it can cause significant noise in my results, I usually use pinv for this but it seems to not be working this time.
Thank you for your time and attention.
(Attached files are the left and right hand side of the linear problem)
Bjorn Gustavsson on 22 Jul 2022
As best I can interpret your flow-chart it seems that it ought to be "reasonably straightforward" to follow Torsten's advice by converting the equation for into one for both and . If I get it right it should be something like:
function Mout = modified_massmatrix(M)
Mout = [eye(size(M,1)),zeros(size(M));
or if the mass-matrix is a function of t, q and :
function Mout = modified_massmatrix(t,qqdot,M)
Mnow = M(t,qqdot);
Mout = [eye(size(Mnow,1)),zeros(size(Mnow));
Then you should be good to go with the ODE-integrating functions.