Bead sliding on rotating rod - Lagrange Mechanics
The position of the bead is given by two coordinates : phi and rho - angle and radius in polar coordinates.
The 2x 2nd order Equations of motion are derived with Lagrange2 formalism.
In order to be solved numerically by matlab buildin ODE solvers the equations have to be
linearized to 4x 1st order ODE's as follows:
w - angular frequency
y0=[0.01 0 0 w]; % [ r dr/dt phi dphi/dt ] initial conditions at t=0
tspan = [0 60];
f=@(t,y)[y(2);w^2 * y(1) ; y(4);0];
[t,y]=ode45(f,tspan,y0); % call ode45 solver
The .zip file contains a mp4-video of the animation.
Cite As
Lucas Tassilo Scharbrodt (2024). Bead sliding on rotating rod - Lagrange Mechanics (https://www.mathworks.com/matlabcentral/fileexchange/67999-bead-sliding-on-rotating-rod-lagrange-mechanics), MATLAB Central File Exchange. Retrieved .
MATLAB Release Compatibility
Platform Compatibility
Windows macOS LinuxCategories
- Sciences > Physics > General Physics >
- Engineering > Mechanical Engineering > Statics and Dynamics >
Tags
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!Discover Live Editor
Create scripts with code, output, and formatted text in a single executable document.
bead on rotatig rod/
Version | Published | Release Notes | |
---|---|---|---|
1.0.0.0 |