Rotating Spring-Mass system ODE45 function

5 次查看(过去 30 天)
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 个评论
Jan
Jan 2021-3-19
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?

请先登录,再进行评论。

回答(1 个)

William Rose
William Rose 2021-3-28
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.

类别

Help CenterFile Exchange 中查找有关 Ordinary Differential Equations 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by