3D Pendulum

A point mass m = 1 m = 1 is attached to the free end of a mass-less rigid rod of length L = 1 L = 1 . The other end of the rod is attached to the origin with a ball joint, allowing the free end to be positioned at an arbitrary point in 3 D 3D space (constrained by the length of the rod, or course). The gravitational acceleration is 10 10 in the negative z z direction.

Initially, the rod is rotating in the x y xy plane with angular speed ω = 7 \omega = 7 (radians per second) with respect to the z z axis. As the rod moves, what is the minimum z z coordinate reached by its free end?

Bonus: Plot the trajectory


The answer is -0.357.

This section requires Javascript.
You are seeing this because something didn't load right. We suggest you, (a) try refreshing the page, (b) enabling javascript if it is disabled on your browser and, finally, (c) loading the non-javascript version of this page . We're sorry about the hassle.

3 solutions

Let at time t t , the angle that the rod makes with vertical be α α . Then the angular velocity of the point mass is given by ω = k sin 2 α \omega=\dfrac{k}{\sin^2 α} , where k k is a constant (since the line of centrifugal force intersects the z z -axis and the line of weight force is parallel to that axis, the moments of both about that axis are zero).

Let β = d α d t β=\dfrac{dα}{dt} . Then β d β d α = k 2 cos α sin 3 α g sin α β\dfrac{dβ}{dα}=\dfrac{k^2\cos α}{\sin^3 α}-g\sin α

β 2 = 2 g cos α k 2 sin 2 α \implies β^2=2g\cos α-\dfrac{k^2}{\sin^2 α} . Here I assumed that the rod was horizontal initially, so that k = 7 k=7 .

When α α is minimum, β = 0 2 g sin 2 α = 49 cos α cos α = 4 + ( 49 2 g ) 2 49 2 g 2 0.35633646 β=0\implies 2g\sin^2 α=49\cos α\implies \cos α=\dfrac{\sqrt {4+(\frac{49}{2g})^2}-\frac{49}{2g}}{2}\approx 0.35633646 .

Hence the minimum z z coordinate is 0.3563 \boxed {-0.3563} .

Karan Chatrath
Apr 29, 2020

Let the pendulum be parameterized in spherical coordinates such that:

A = θ A = \theta B = ϕ B = \phi

Coordinates of the bob at an arbitrary instant are:

x = sin A cos B x = \sin{A} \cos{B} y = sin A sin B y = \sin{A} \sin{B} z = cos A z = \cos{A}

The kinetic energy and potential of the rod is:

T = 1 2 ( x ˙ 2 + y ˙ 2 + z ˙ 2 ) \mathcal{T} = \frac{1}{2}\left(\dot{x}^2+\dot{y}^2+\dot{z}^2\right) V = 10 cos A \mathcal{V} = 10\cos{A}

Using the Lagrangian approach or the energy conservation principle, the equations of motion can be derived which are (working left out):

A ¨ = ( sin A cos A ) B ˙ 2 + 10 sin A \ddot{A} = (\sin{A}\cos{A}) \dot{B}^2 + 10\sin{A} B ¨ = 2 B ˙ A ˙ cot A \ddot{B} = -2\dot{B}\dot{A} \cot{A}

Initial conditions are: A ( 0 ) = π 2 A(0) = \frac{\pi}{2} B ( 0 ) = 0 B(0) = 0 A ˙ ( 0 ) = 0 \dot{A}(0) = 0 B ˙ ( 0 ) = 7 \dot{B}(0) = 7

Numerical integration does the rest beyond this point. However, I am presenting an analytical approach similar to that of @Alak Bhattacharya as follows. Consider the second equation of motion

B ¨ = 2 B ˙ A ˙ cot A \ddot{B} = -2\dot{B}\dot{A} \cot{A} B ¨ = d B ˙ d t = d B ˙ d A d A d t = 2 B ˙ A ˙ cot A \implies \ddot{B} = \frac{d\dot{B}}{dt}= \frac{d\dot{B}}{dA} \frac{dA}{dt}= -2\dot{B}\dot{A} \cot{A} d B ˙ d A = 2 B ˙ cot A \implies \frac{d\dot{B}}{dA} = -2\dot{B} \cot{A}

Separating variables and integrating gives:

ln B ˙ = 2 ln ( sin A ) + ln K \ln{\dot{B}} = -2\ln\left(\sin{A}\right) + \ln{K}

Where K K is an integration constant. SImplifying gives:

B ˙ = K sin 2 A \dot{B} = \frac{K}{\sin^2{A}}

Applying initial condition yields that K = 7 K=7 . Thus:

B ˙ = 7 sin 2 A \dot{B} = \frac{7}{\sin^2{A}}

Using the above result in the first equation of motion:

A ¨ = ( sin A cos A ) B ˙ 2 + 10 sin A \ddot{A} = (\sin{A}\cos{A}) \dot{B}^2 + 10\sin{A}

A ˙ d A ˙ d A = sin A cos A ( 49 sin 4 A ) + 10 sin A \implies \dot{A}\frac{d \dot{A}}{dA} = \sin{A}\cos{A} \left(\frac{49}{\sin^4{A}}\right) + 10\sin{A} A ˙ d A ˙ = ( ( 49 cos A sin 3 A ) + 10 sin A ) d A \implies \dot{A} \ d\dot{A} = \left(\left(\frac{49\cos{A}}{\sin^3{A}}\right) + 10\sin{A}\right)dA

Integrating both sides again:

A ˙ 2 = 49 sin 2 A 20 cos A + C \dot{A}^2 = -\frac{49}{\sin^2{A}} - 20\cos{A} + C

Applying initial conditions again:

A ˙ 2 = 49 sin 2 A 20 cos A + 49 \dot{A}^2 = -\frac{49}{\sin^2{A}} - 20\cos{A} + 49

Finally, consider:

z = cos A z = \cos{A}

For z z to be minimum:

z ˙ = sin A A ˙ = 0 \dot{z} = -\sin{A} \ \dot{A} = 0

A ˙ = 0 \implies \dot{A} = 0 49 sin 2 A 20 cos A + 49 = 0 \implies -\frac{49}{\sin^2{A}} - 20\cos{A} + 49 = 0

This leads to a quadratic equation in cos A \cos{A} . Solving and realising that A > π / 2 A> \pi/2 since the bob moves downward results in:

cos A 0.3563 \cos{A} \approx -0.3563

Hence the minimum Z-coordinate is:

z m i n = 0.3563 z_{min} = -0.3563

Missed the bit about the trajectory. Here it is:

One can see the formation of an intricate pattern in 3D space. I have given three fundamental views here.

Karan Chatrath - 1 year, 1 month ago

A 3D view:

The system is simulated for 10 seconds here.

Karan Chatrath - 1 year, 1 month ago
Steven Chase
Apr 29, 2020

Thanks to @Alak Bhattacharya and @Karan Chatrath for the detailed solutions. I will simply add some plots of the trajectory in the x z xz plane. The plots are for 10 10 seconds, 30 30 seconds, and 100 100 seconds. The trajectory is essentially confined to a "band" on the surface of a sphere.

@Steven Chase Nice problem. Post more problems in pendulum

A Former Brilliant Member - 1 year, 1 month ago

Log in to reply

@Steven Chase How can we parametrise ellipsoid and write its d A dA ??

A Former Brilliant Member - 1 year, 1 month ago

0 pending reports

×

Problem Loading...

Note Loading...

Set Loading...