Model Dynamics of Unit Mass Sliding Down y=f(x) Under Gravity

In summary, we can model a mass sliding down a stationary curve as sliding down a stationary inclined plane which is tangential to the curve at the point of contact. The slope of this inclined plane is equal to the derivative of the curve at that point. We can use vector analysis to find the net force exerted on the mass, which includes the force of gravity, the force parallel to the plane, the force perpendicular to the plane, and the normal force. However, there was a flaw in the model due to an angle error in the derivation of the force of friction. This can be corrected by using the multi-parameter version of the tangent function or by finding the normal vector in a different way.
  • #1
LuculentCabal
20
0
For a unit mass sliding down the stationary curve [itex]y = f(x)[/itex] at the point [itex](x, y)[/itex], we can consider our mass to be sliding down a stationary inclined plane which is tangential to the curve at the point [itex](x, y)[/itex]. The slope of this inclined plane is thus [itex]\frac{dy}{dx}(x)[/itex].

For the remainder of this thread, [itex]y'[/itex] means [itex]\frac{dy}{dx}[/itex]

For the purposes of this excersize, I want the coefficient of static friction to be zero.
Let [itex]\mu_k[/itex] be the coefficient of kinetic friction
Let [itex]g[/itex] be the acceleration due to gravity
Let [itex]\vec{F_g} = mg[/itex] be the force exerted on our unit mass due to gravity.
Let [itex]\vec{F_{||}}[/itex] be the force exerted on our mass parallel to the plane.
Let [itex]\vec{F_\bot}[/itex] be the force exerted on our mass perpendicular to and into the plane.
Let [itex]\vec{F_N} = - \vec{F_\bot}[/itex] be the force exerted on our mass perpendicular to and out of the plane .
Let [itex]\vec{F_f} = \mu_k \vec{F_\bot}[/itex] be the force of friction which resists [itex]\vec{F_{||}}[/itex]
Then [itex]\vec{F_Ʃ} = \vec{F_{||}} + \vec{F_f} + \vec{F_\bot} + \vec{F_N}[/itex] is the net force exerted on our mass. This is our objective.
---------------------------------------------------------------------------------------------------------------
From the vector analysis, I have found [itex]\vec{F_Ʃ}[/itex] to be given by:
[itex]\frac{1}{g} \vec{F_Ʃ} = [cos(θ) sin(θ) - \mu_k cos^2(θ)] \ \hat{i} + [sin^2(θ) - \mu_k cos(θ) sin(θ)] \ \hat{j}[/itex]

Where:
[itex]θ = tan^{-1} (y')[/itex]


The model behaves as expected whenever [itex]\mu_k[/itex] = 0.
However, consider the trivial cases 1. and 2. below:

Case 1.
[itex]\ \ \ g = -9.81[/itex] // arbitrarily
[itex]\ \ \ \mu_k = 0.5[/itex] // arbitrarily
[itex]\ \ \ y = x[/itex] // the standard inclined plane

Case 2.
[itex]\ \ \ g = -9.81[/itex] // same as Case 1.
[itex]\ \ \ \mu_k = 0.5 [/itex] // same as Case 1.
[itex]\ \ \ y = -x[/itex] // this is simply observing Case 1. from the other side of the plane

These parameters result in the following dynamics:

Case 1.
[itex]\ \ \ \frac{d^2 x}{d t^2} = -2.4525[/itex]
[itex]\ \ \ \frac{d^2 y}{d t^2} = -2.4525[/itex]

Case 2.
[itex]\ \ \ \frac{d^2 x}{d t^2} = 7.3575[/itex]
[itex]\ \ \ \frac{d^2 y}{d t^2} = -7.3575[/itex]

---------------------------------------------------------------------------------------------------------------

For the previous cases, my intuition told me to expect the following:

a) That the sign on the horizontal acceleration would be negative in Case 1. and positive in Case 2.
b) That the sign on the vertical acceleration would be negative in Case 1. and negative in Case 2.
c) That the magnitude of the horizontal acceleration in Case 1. would be equal to the magnitude of the vertical acceleration in Case 1.
d) That the magnitude of the horizontal acceleration in Case 2. would be equal to the magnitude of the vertical acceleration in Case 2.
e) That the magnitude of horizontal and vertical acceleration in Case 1. would be equal to those of Case 2.

As we can see, the results are in agreement with a), b), c), d). However, hypothesis e) has failed.

Question:
What is the flaw in the model (or hypothesis e) ) that made the results diverge from e) when [itex]\mu_k[/itex] was non-zero (holding all other parameters constant)?

**If it is helpful, I can post my derivation of [itex]\vec{F_Ʃ}[/itex] in a followup.**
 
Physics news on Phys.org
  • #2


In case 2, your friction force has the same direction as the component of gravity parallel to the plane.
 
  • #3


willem2 said:
In case 2, your friction force has the same direction as the component of gravity parallel to the plane.

Hi willem2. Thanks for pointing that out.

In trying to find out why that is the case, I found an angle error in my derrivation of the force of friction. I hope this is the only bug...
 
  • #4


Well it seemed like a good idea, but it appears as though we can't have signed angles
ie [itex]θ ≠ tan^{-1}(y')[/itex] but rather [itex]θ = |tan^{-1}(y')| = tan^{-1}(|y'|)[/itex]
 
  • #5


That's sounds like a classic inverse trig obstacle that we encounter all the time in computer science (I'm a programmer by trade).

Have you considered using the multi-parameter version of the tangent function?

i.e. arctan(x, y) rather than arctan(y/x)

It's used in programming libraries to correct the deficiencies of the traditional arctan. arctan(y/x) essentially loses some of the information of the system. Human beings will often subconsciously reinsert lost information into the system when they use arctan(y/x) by hand, if my memory is correct.

Also, there are other ways of finding the normal vector than using trig. For example, in 2 dimensions you can convert a vector to a complex number (x = real component, y = imaginary component) and then multiply by i in order to rotate the vector by 90 degrees to get the normal vector. Quaternions can achieve a similar function for 3 dimensions I believe, from what I've heard.

Hope that helps. I only glanced over the other parts of your post but these few things caught my eye quick.
 
  • #6


Thank you for taking the time to offer your suggestions WraithGlade.

Unfortunately the problem with the model that I have presented is that displacement instead of distance was used for the purpose of determining theta. It is because of this that I am getting signed angles which are problematic for the sine function being that it is odd. It looks inevitable that fixing this bug is going to introduce some form of absolute value function.

I will have to play around with it now to solve for x(t). However, it looks like Newton-Raphson will still be king at the end of the day!
 
  • #7


LuculentCabal said:
Let [itex]\vec{F_\bot}[/itex] be the force exerted on our mass perpendicular to and into the plane.
You mean, that component of gravity, yes?
Let [itex]\vec{F_N} = - \vec{F_\bot}[/itex] be the force exerted on our mass perpendicular to and out of the plane .
This one is exerted by the wire, yes?
They're only going to be equal and opposite if the acceleration is parallel to the wire.
Where the wire is curved, there will be acceleration normal to the wire.
I note the equation you end up with has no term involving y'', which cannot be right.
 
  • #8


Hello haruspex. Thank you for your review.

Yes, [itex]\vec{F_\bot}[/itex]is due to gravity.

I would be grateful if you (or someone) could explain why our mass would experience acceleration normal to the curve on the points where [itex]\frac{d^2 y}{dx^2} ≠ 0[/itex].

I had imagined that at the point (x,y), our mass would experience acceleration parallel to a tangent at the point (x,y). As our mass moves instantaneously along the curve, our tangent would be instantaneously changing too as it moves along with the mass. Why is this not the case?
 
  • #9


LuculentCabal said:
Hello haruspex. Thank you for your review.

Yes, [itex]\vec{F_\bot}[/itex]is due to gravity.

I would be grateful if you (or someone) could explain why our mass would experience acceleration normal to the curve on the points where [itex]\frac{d^2 y}{dx^2} ≠ 0[/itex].

I had imagined that at the point (x,y), our mass would experience acceleration parallel to a tangent at the point (x,y). As our mass moves instantaneously along the curve, our tangent would be instantaneously changing too as it moves along with the mass. Why is this not the case?

Well, as for arbitrary two dimensional curvilinear motion, the acceleration components are:
[tex]a_{n}=\frac{v^2}{r} [/tex]
[tex]a_{t}=\frac{dv}{dt} [/tex]
The first refers to the acceleration components normal to the tangent and thus the second refers to the acceleration components parallel to the tangent.The r is radius of curvature, which, for an arbitrary two dimensional curve, has a specific formula for it.So here you may see that you have to know the velocity of the mass so as to get the [itex]a_{n}[/itex]
 
  • #10


LuculentCabal said:
I would be grateful if you (or someone) could explain why our mass would experience acceleration normal to the curve on the points where [itex]\frac{d^2 y}{dx^2} ≠ 0[/itex].
Just to add a little to what azureth wrote:
Acceleration in the direction of travel will change the speed but not the direction. To travel along a curve requires a component of acceleration normal to the velocity, a centripetal force.
 
  • #11


Sorry guys, but I have no idea how to find v if all that I know is y(x), g, mu_k, and the initials v0, x0. It seems v cannot be found without first either knowing x(t) or F_x(t) and neither x(t) or F_x(t) can be found without knowing v...
 
  • #12


Because there is friction involved, I think you have to start with the acceleration equations. Otherwise you could have used the energy equation to relate y and v.
 
  • #13


LuculentCabal said:
For a unit mass sliding down the stationary curve [itex]y = f(x)[/itex] at the point [itex](x, y)[/itex], we can consider our mass to be sliding down a stationary inclined plane which is tangential to the curve at the point [itex](x, y)[/itex]. The slope of this inclined plane is thus [itex]\frac{dy}{dx}(x)[/itex].

If by "sliding" you really mean that the mass can move only along the curve, it would be easier to consider the parametric representation of the curve, [itex]x = x(s), y = y(s)[/itex], where [itex]s[/itex] is the length along the curve. Then what you would need to find is [itex]s = s(t)[/itex]. The magnitude of the velocity at any point is [itex]\dot{s} = \frac {ds}{dt}[/itex], so [itex]m\ddot{s} = F_{\tau}[/itex], where the entity on the right is the tangential component of the net force. [itex]F_{\tau} = G_{\tau} - R[/itex], where [itex]G_{\tau}[/itex] is the tangential component of the gravitational force, and [itex]R[/itex] is the magnitude of the force of friction. [itex]R = \mu_kN[/itex], where [itex]N[/itex] is the magnitude of the reaction force. The reaction force is always normal to the trajectory, so [itex] ma_n = N[/itex], where [itex]a_n[/itex] is the normal acceleration. Now, [itex]a_n = v^2\kappa[/itex], where [itex]v[/itex] is the magnitude of the velocity and [itex]\kappa[/itex] is the curvature, so [itex]R = \mu_k m \dot{s}^2 \kappa = \mu_k m \dot{s}^2 \left| x'y'' - y'x'' \right|[/itex]
 
  • #14


Continuing, [itex]G_{\tau} = - mg sin \theta = -mg y'[/itex]. Putting everything together: [tex]m\ddot{s} = -mg y' - \mu_k m \dot{s}^2\left|x'y'' - y'x''\right|[/tex]
Let's see how that works for a unit circle. [itex]x = cos(s), y = sin(s)[/itex]. Then:[tex]\ddot{s} = -g cos(s) - \mu_k m \dot{s}^2[/tex] In the vicinity of [itex]s = \frac {3 \pi } {2}[/itex], [itex] cos (s) \approx s [/itex], and assuming small velocities [itex]\dot{s}^2 \approx 2\dot{s}[/itex], thus [tex]\ddot{s} + 2\mu_k \dot{s} + gs = 0[/tex] which is the well-known equation of damped harmonic oscillations, as is to be expected physically.
 

1. What is the purpose of modeling unit mass sliding down y=f(x) under gravity?

The purpose of this model is to study the motion of a unit mass as it slides down a surface with a given shape and under the influence of gravity. This can help us understand the behavior of objects in real-world scenarios, such as a ball rolling down a hill or a car driving down a curved road.

2. How is this model different from other models of motion?

This model specifically focuses on the motion of a unit mass sliding down a surface with a given shape. Other models may take into account additional factors such as air resistance or friction, but this model simplifies the scenario to only consider the effects of gravity.

3. Can this model be applied to real-world situations?

Yes, this model can be applied to real-world situations where an object is sliding down a surface under the influence of gravity. However, the accuracy of the model may vary depending on the specific scenario and other factors that may affect the motion.

4. What are the variables involved in this model?

The variables involved in this model are the position of the unit mass on the surface (x and y coordinates), the shape of the surface (represented by the function f(x)), the initial velocity of the mass, and the acceleration due to gravity.

5. How can this model be used to make predictions?

By inputting different values for the variables, this model can be used to predict the motion of the unit mass as it slides down the surface. For example, by changing the initial velocity or the shape of the surface, we can see how it affects the speed and position of the mass at different points in time.

Similar threads

Replies
6
Views
2K
  • Introductory Physics Homework Help
Replies
29
Views
985
  • Introductory Physics Homework Help
Replies
25
Views
212
  • Introductory Physics Homework Help
Replies
2
Views
186
  • Calculus and Beyond Homework Help
Replies
3
Views
488
Replies
2
Views
816
  • Introductory Physics Homework Help
Replies
11
Views
169
Replies
1
Views
747
  • Introductory Physics Homework Help
Replies
1
Views
271
  • Introductory Physics Homework Help
Replies
9
Views
859
Back
Top