Solving The Equations for a Spinning Top in Python

In this video I simulate and animate a 3D spinning top by solving Lagrange's equations in python, and animating the spinning top in Matlab.

Finding Moment of Inertia Tensor for Any Object:

Code:

Discord:

[kashu]

this is so cool, your videos are seriously valuable keep it up man

[cambridgebreaths]

Man I like your tutorials. Want to binge watching soon. Wish me luck. All the best for your channel.

[ricardoferes]

Hey man, I am a mechanical engineer from Brazil, and this content is awesome.

[chachachachang]

So, is this like a gaming channel for solving differential equations? As an instructor for ODEs I really appreciate that you are stepping up the challenge for all of us. Great job!

[andresderudder]

Hi Mr, I study physics in UBA (Argentina) and I'm new in the channel, your videos are great!

[nestorcaticha]

Very nice lecture, informative and clear. A minor comment is that omega $\vec \omega $ is the angular velocity, not the angular momentum, given by $I \vec \omega$

[hussainhassan]

when I used this
sols = smp.solve([LE1, LE2, LE3], (the_dd, phi_dd, psi_dd),
simplify=False, rational=False)
and
sols[the_dd]
it continues to run for more than 1 hour and it did not stop!
thanks

[taybabusra]

Watching your video regularly, can you please make video on epicycle drawing details on python

[jcmschott]

Really nice stuff, Mr P Solver! I was just wondering if it was possible to implement this animation in Blender, maybe by exporting the data containing the solution and key-framing a spinning top modeled in Blender.
Another question: have you ever wondered about using Blender and all its shiny, glossy and nodes-editor features to visualize your data?
Once more, thank you so much for this huge inspiration and all the nice content!

[frankkoslowski]

Same problem as with the 3D-Pendulum.
The expressions sympy.solve([LE1, LE2, LE3], (the_dd, phi_dd, psi_dd), simplify=False, rational=False)
is expected to tackle are soo long that the process of finding a symbolic solution was shown to take seemingly forever.
Are you using a professional version of Python implementing precompiled code capable of streamlining the process perhaps?
You are also doing something rather unorthodox here: sols[the_dd]
namely, indexing with a symbolic expression of the 2nd derivative of Theta, rather than an Integer commonly used for slicing. 😭

[NaneRulz]

Great content. How about using fenics as pde silver?

[felixlupp]

It doesn't work for me I get TypeError: can't convert expression to float in the odeint function. anyone know how to fix this?

[aqueeqazam]

You are a carbon copy of elon musk man. I am from India and I am great fund of your tutorial.