Rotating Spring-Mass system ODE45 function

6 views (last 30 days)
Tyler
Tyler on 19 Mar 2021
Answered: William Rose on 28 Mar 2021
I need to plot the motion of a mass sliding in a slot in a circular disk(gravity into page) while connecting to a spring centered at the center of the disk. The plot is supposed to show x vs y rather than the x or y components vs time. Does anyone know the proper x and y equations of motion that would apply here?
  1 Comment
Jan
Jan on 19 Mar 2021
The description as text is not clear enough. Did you made a drawing already?
Creating the equations on motion is not a Matlab problem, so it is off-topic in this forum.
Do you want to solve this by a DAE with depending coordinates or an ODE with independent coordinates?

Sign in to comment.

Answers (1)

William Rose
William Rose on 28 Mar 2021
Find x(t) and y(t) then plot x vs y.
You did not mention friction so we assume there is none. So the equations of motion could be written
d2x/dt2 = f_x(x,y)/M
d2y/dt2 = f_y(x,y)/M
The forces include the inward spring force and the centrifugal force in this rotating frame of reference.
The spring force and centrifugal force are radial. Therefore, rather than trying to develop the equations above, I would express both forces in terms of r, and solve the differential equation for r(t).
d2r/dt2 = f(r)/M
where f(r) = f_spring(r) + f_centrifugal(r)
You did not say if the angular velocity is constant. The angular velocity versus time will determine what the centrifugal force looks like and will determine theta versus time. Once you know r(t) and theta(t), convert to x(t) and y(t) in the standard way. If the angular velocity is constant then it is pretty easy. If the initial angle is not specified, you might as well assume angle theta=0 at t=0.
That should be enough hints to get you going.

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!