Solving a Cycloid Equation: Finding θ(t)

  • Thread starter thomas91
  • Start date
42 replies · 7K views · last reply
thomas91
Messages
23
Reaction score
0
Missing homework template due to originally being posted in technical forum
Hello! Nice to meet you. I'm studying a Bachelor of Science and I have a doubt:

A bead slides without friction on a frictionless wire in the shape of a cycloid with equations

x = a(θ-sin θ) y = a(1+cos θ)

where θ's range is 0 to 2π and a>0

a) Find the equation of motion:

And I resolved:
Equation: (1 - cos(θ)) θ'' + (1/2)(θ')² sin(θ) - (g / (2a)) sin(θ) = 0

I checked in a physics book that the solution is correct

b) Make the change of variable u= cos (θ/2) and find a linear equation for u(t). Solve this equation and find the solution θ(t) which in the case of a bead at t = 0 falls from rest from the corresponding point θ= π/2

I don't know how to make the change. Can anybody help me? Thanks and nice to meet you.
 
Exactly what in the change of variables are you having problems with? It is a matter of reexpressing the differential equation in terms of u instead of theta.
 
Orodruin said:
Exactly what in the change of variables are you having problems with? It is a matter of reexpressing the differential equation in terms of u instead of theta.

Thanks for answering! I do not understand how I have to make the change of variable, I'm trying to do since Wednesday and not even know where to start the change of variable
 
It is just a normal change of variables. You insert the new variable in place of the old one using the relation between them.
 
Orodruin said:
It is just a normal change of variables. You insert the new variable in place of the old one using the relation between them.

I tried doing a simple change of variables. But I obtain a non linear equation..

Equation: u'' + 1/(u + 1) (u')^2 + g/(4a) = 0, written as a second-order nonlinear differential equation in u(t).
 
Please show your work and how you got there. Your non-linear term should not be there. The other terms are correct.
 
Orodruin said:
Please show your work and how you got there. Your non-linear term should not be there. The other terms are correct.

Equation: (1 - cos(θ)) θ'' + (1/2)(θ')^2 sin(θ) - (g / (2a)) sin(θ) = 0, with θ'' and θ' the second and first derivatives of θ.


u= cos (θ/2) then:

θ = 2 arccos u

then:

Equation θ′ = −(2 × u′ / √(1 − u²)), with the fraction's numerator u′ and denominator the square root of 1 minus u squared.


and:
Equation θ'' = -2 times a bracket containing u'' divided by the square root of (1 - u^2), plus u' divided by (1 - u^2) raised to the 3/2 power.


Then substitute:

Equation: (1 - cos(2 cos^-1(u))) θ'' + (1/2) sin(2 cos^-1(u)) ((θ')^2 - g/a) = 0


Now:
Trigonometric identity written as cos 2a = cos²a − sen²a.
and
Trigonometric identity written as sen 2a = 2 sen a cos a, with the double-angle formula for sine.
and substitute:

Equation: (1 - (cos^2(cos^-1(u)) - sin^2(cos^-1(u)))) θ'' + (1/2) × 2 (sin(cos^-1(u)) cos(cos^-1(u))) ((θ')^2 - g/a) = 0


We now that cos(arccos (a)) = a and sen(arccos(a) =
Square root of (1 minus a squared), written as sqrt(1 - a^2).


Equation: 2(1 - u^2)θ'' + u√(1 - u^2)((θ')^2 - g/a) = 0


Using first and second derivates and simplifying:

Equation: -4(u'' sqrt(1 - u^2)) + 4(u - 1)(u')^2 × 1/sqrt(1 - u^2) - (g/a) sqrt(1 - u^2) = 0


Then multiplied by
Fraction with numerator 1 and denominator the square root of 1 minus u squared, written as 1 over sqrt(1 - u^2).
and dividing by (-4) we obtain:

proxy.php?image=http%3A%2F%2Fi.imgur.com%2FJHSWflM.png


Which isn't linear.What is your opinion? Thanks!
 
It is much easier to help you if you stop pasting images and instead use the forum LaTeX to type in your equations (they become quotable).

Your second derivative is wrong, the second term needs to contain ##(u')^2##. Also, it becomes much much easier if you do not invert the relation ##u = \cos(\theta/2)## but instead differentiate it as it stands.
 
Orodruin said:
It is much easier to help you if you stop pasting images and instead use the forum LaTeX to type in your equations (they become quotable).

Your second derivative is wrong, the second term needs to contain ##(u')^2##. Also, it becomes much much easier if you do not invert the relation ##u = \cos(\theta/2)## but instead differentiate it as it stands.

¿##(u')^2##? I don't see this.. I don't understand what you mean. Sorry.
Thanks Orodruin!
 
Now I think that the second derivative is :
Equation θ'' = -2 ( u'' / sqrt(1 - u^2) + u u' / (1 - u^2)^(3/2) ), with u'' and u' as primed terms and θ'' on the left.

But it doesn't contains ##(u')^2##
 
thomas91 said:
Now I think that the second derivative is :
Equation θ'' = -2 times the sum of two terms in parentheses: u''/√(1 − u²) plus u u'/(1 − u²)^(3/2).

But it doesn't contains ##(u')^2##
It does, you did it wrong again. What is the derivative of
$$
\frac{1}{\sqrt{1-u^2}}?
$$
 
thomas91 said:
##uu' ((1-u^2))^{-3/2}##
It's that right?
So what is then the derivative of
$$
\frac{u'}{\sqrt{1-u^2}}?
$$
 
Orodruin said:
So what is then the derivative of
$$
\frac{u'}{\sqrt{1-u^2}}?
$$

Mathematical expression: u'' divided by the square root of (1 - u^2), plus (u')^2 divided by (1 - u^2) raised to the power 3/2, all enclosed in large parentheses.

It's ok?

Thanks!
 
What happened to the ##u## coming out of the derivative of ##1/\sqrt{1-u^2}##? But again, the best way forward is not to solve for ##\theta## before differentiating, it is to directly differentiate ##u = \cos(\theta/2)##.
 
Orodruin said:
What happened to the ##u## coming out of the derivative of ##1/\sqrt{1-u^2}##? But again, the best way forward is not to solve for ##\theta## before differentiating, it is to directly differentiate ##u = \cos(\theta/2)##.

Hello! If you do not mind, I prefer to finish the exercise by this method and later solve by the method you suggest

Then, the second derivative is:
Mathematical expression: -2 times the sum of u'' divided by the square root of (1 - u^2), plus u(u')^2 divided by (1 - u^2) raised to the 3/2 power, all in parentheses.


Replacing first and second derivative:

Handwritten-style equation set equal to zero involving u, u', u'' and g/a, with fractional powers and a square root term of 1 - u^2 in the denominator.


And simplifying:

Handwritten equation: (-4u'')√(1-u²) - u√(1-u) × g/a = 0
Simplifying again:

Differential equation u'' + u(g/(4a)) = 0, written in plain text with the second derivative of u plus u times g over 4a equal to zero.
What do you think? Thanks and have a nice sunday
 
Looks fine apart from a square disappearing in your middle step (the end result is correct).
 
Orodruin said:
Ah yes, and the ##a## should be in the denominator in the end result ...

Sorry! It was a mistake...
The linear equation is:
Differential equation u'' + u(8/(4a)) = 0, written as a second-order linear equation with constant coefficient 8/(4a).


The characteristics polynomial is:
Characteristic equation λ² + g/(4a) = 0, written as an equation with λ squared plus g over 4a equal to zero.


With complex roots:

Equation: lambda = -i sqrt(g) / (2 sqrt(a)), written as a complex root of the characteristic polynomial.


Equation lambda = i times the square root of g divided by 2a, written as a fraction with a square root in the numerator and 2a in the denominator.


The general solution:

Handwritten solution u = A cos((1/2) sqrt(g/a) t) + B sin((1/2) sqrt(g/a) t).


Equation: cos(θ/2) = A cos((1/2)√(g/a) t) + B sin((1/2)√(g/a) t), the general solution of a second-order linear ODE with cosine and sine terms in t.


We need to know the particular solution. True?Thanks!
 
Edit: The equation is homogeneous. No particular solution. But then we can not replace, I can not find the solution
 
What do you mean by "cannot replace"? The homogenous equation has constants you have to adapt to the initial conditions. That will give you the solution you want.
 
Orodruin said:
What do you mean by "cannot replace"? The homogenous equation has constants you have to adapt to the initial conditions. That will give you the solution you want.
Hello!

Equation: cos(θ/2) = A cos((1/2) sqrt(g/a) t) + B sin((1/2) sqrt(g/a) t)


I think this:
a) Find the solution θ(t) which in the case of a bead at t = 0 falls from rest from the corresponding point θ= π/2

Handwritten equation: cos(π/2) = A cos(0) + B sin(0), with the left side cos of the fraction π over 2 inside a bracket.


Then:
Handwritten equation: 1/√2 = A


What do you think?
 
Orodruin said:
This settles A. You also need to find B.

But ... B is out of the equation because sin (0) = 0
Then B may be any value. Yes?

Thanks!
 
For that particular value of t it does not affect the position. But this is a second order differential equation so you need two conditions (and you have one more!).
 
Orodruin said:
For that particular value of t it does not affect the position. But this is a second order differential equation so you need two conditions (and you have one more!).

Hello! I can't find the other condition. Can you give me a clue?

Thanks
 
Can I set two arbitrary (Between 0 to 2π ) differentes values of θ and solve a equation system for t and B?
 
thomas91 said:
Can you give me a clue?
thomas91 said:
b) Make the change of variable u= cos (θ/2) and find a linear equation for u(t). Solve this equation and find the solution θ(t) which in the case of a bead at t = 0 falls from rest from the corresponding point θ= π/2

My highlighting.
 
Orodruin said:
My highlighting.

Does the rest position mean that the velocity mean that θ' = 0 for t=0?

I tried to derivate θ and solve B for t=0 and I obtain B=0...
 
thomas91 said:
Does the rest position mean that the velocity mean that θ' = 0 for t=0?
Yes.

thomas91 said:
I tried to derivate θ and solve B for t=0 and I obtain B=0...
So now you can insert this into your solution and solve for theta to get ...?