Motivation of suggested investigation stems from the fact that the proposed assembly is composed of a sliding cube that acts as a point-like object with only linear kinematics in contrast to the rod that in addition to the former possesses rotational degrees of freedom. The length of the rod adds additional features to the physics of the problem in contrast to the movement of the pair point-like objects in the two-body problem reported in [3] . The leaning rod against the cube constitutes a holonomic constraint reducing the number of degrees of freedom. These two features, namely, 1) rotation of the rod and 2) its constraint movement sparked the investigation. The problem appears to be a 2D two-body problem with numerous degrees of freedom however; the restricted motion of the rod practically reduces to a problem with only one degree of freedom. In so many words, the kinematic information about the entire system comes about from the solution of a single nonlinear differential equation with its dependent variable being the inclination angle of the rod w/floor. The equation of motion is analytically unsolvable, proving the effective usefulness of a Computer Algebra System, specifically Mathematica. Utilizing the numeric solution in hand, we quantify a list of kinematic information primary as noted the separation runtime of the rod-cube system. Additional quantities such as speed, acceleration and contact forces of the bodies are also quantified as well. This article is composed of three sections. In addition to Introduction in Section 2 applying Lagrangian and Hamiltonian methods, we develop the needed formulation conducive to equation of motion. This section includes sub-sections: Rod-Cube contact force and Trajectory of the tip of the rod. Section 3 is the Results. In this section, we utilize Mathematica numeric utilities solving the needed equations conducive to quantifying kinematic and dynamic quantities of interest. We close the report with concluding remarks suggesting ideas furthering the investigation.

Figure 1 depicts the schematic of the problem on hand.

As shown in Figure 1, one end of a uniform rod of length ℓ and mass m₁ touches the horizontal surface at x₁, its other end leans against the left vertical face of the cube at C(x₂, y₂). The cube has a mass of m₂ with side length a. The rod and the center of mass (cm) of the cube are on a vertical plane. Relevant forces acting on the rod namely the contact with the floor, N, the weight, W, and the one with the cube, F₂₁, are also shown.

Releasing the assembly from rest sets the system in motion. This is because for the frictionless surfaces there is an unbalanced horizontal contact force acting on the cube, not shown in Figure 1; it acts opposite to shown F₂₁. As depicted, the same value of the force F₂₁ acts on the rod pushing it to the left setting it in motion. However, as the rod slides its center of mass falls and the rod rotates about the sliding contact point x₁. While the rod falls its “circulation radius” shortens hence separation is warranted. In a different scenario if the rod was hinged at x₁ so that circulation radius would have stayed constant separation would have occurred. For either scenario one speculates the separation runtime would depend to the masses in motion and respective dimensions of the bodies. It is one of the objectives of this investigation to quantify separation runtime. As we were developing the formulation of the problem also we quantify a list of interesting quantities. The proposed problem is different from the “classic” version [4] [5] where the rod leans against a vertical wall.

To bypass the dynamic quantities, i.e. forces, for time being we formulate the problem following Lagrangian method. Lagrangian of the system is composed of two pieces, one for the cube and the other one for the rod. Kinetic energy of the

cube is T c u b e = 1 2 m 2 x ˙ 2 2 and its potential energy with respect to the table top is

V c u b e = m 2 g a 2 . Kinetic energy of the rod in addition to a 2D motion of the cm includes its rotational energy about the cm. These two terms are,

T r o d = 1 2 m 1 v c m 2 + 1 2 I c m ω 2

where ω equals θ ˙ is the angular velocity of the rod about the cm; I c m = 1 12 m 2 l 2 . Its potential energy is V r o d = m 1 g y 1 c m ; g is the gravity constant; according to tradition dotted quantities are the derivatives w/time. As shown in Figure 1, it is clear that the angular velocity of the rod about its cm is the same as its angular velocity w/contact point; θ is the inclination angle of the rod w/horizontal.

Conveniently the origin is set as shown. Accordingly,

{ x 2 = x 1 + l cos ( θ ) x 1 c m = x 1 + 1 2 l cos ( θ ) y 1 c m = 1 2 l sin ( θ ) (1)

their respective velocities are,

{ x ˙ 2 = x ˙ 1 − l θ ˙ sin ( θ ) x ˙ 1 c m = x ˙ 1 − 1 2 l θ ˙ sin ( θ ) y ˙ 1 c m = 1 2 l θ ˙ cos ( θ ) (2)

conducive to, v c m 2 = x ˙ 2 1 c m + y ˙ 2 1 c m .

Utilizing the last equation and Equation (2), Lagrangian of the system, e.g. L = T r o d + T c u b e − ( V r o d + V c u b e ) gives,

L = 1 2 m 1 [ x ˙ 1 2 + ( l 2 ) 2 θ ˙ 2 − l x ˙ 1 θ ˙ sin ( θ ) ] + 1 2 I c m θ ˙ 2     + 1 2 m 2 [ x ˙ 1 − l θ ˙ sin ( θ ) ] 2 − m 1 g l 2 sin ( θ ) − m 2 g a 2 (3)

Equation (3) reveals the Lagrangian is independent of x₁, meaning x₁ is a cyclic variable; ∂ L ∂ x 1 = 0 . Because ∂ L ∂ x ˙ 1 ≡ P x 1 according to Euler-Lagrange equation d d t P x 1 = 0 warrants the conservation of the linear momentum of the system

along the horizontal direction. Explicitly, momentum of the system P x 1 is,

P x 1 ≡ ∂ L ∂ x ˙ 1 = m 1 [ x ˙ 1 − 1 2 l θ ˙ sin ( θ ) ] + m 2 [ x ˙ 1 − l θ ˙ sin ( θ ) ] = 0 (4)

The RHS of Equation (4) is set zero in accordance the initial rest state of the system. Accordingly, Equation (4) sets a relationship between the two variables, namely,

x ˙ 1 = ( 1 2 m 1 + m 2 m 1 + m 2 ) l θ ˙ sin ( θ ) (5)

Substituting Equation (5) in (3), yields,

L ( θ , θ ˙ ) = 1 2 m 1 l 2 [ − 1 4 m 1 m 1 + m 2 sin 2 ( θ ) + 1 3 ] θ ˙ 2 − m 1 g l 2 sin ( θ ) − m 2 g a 2 (6)

Therefore substituting Equation (6) to Euler-Lagrange equation, d d t ∂ L ∂ θ ˙ − ∂ L ∂ θ = 0

yields the equation of motion. Equation of motion is a challenging second order ODE. Although it can be solved numerically, however, the primary objective is not to identify θ = θ(t), but rather t = t(θ). Therefore, we make a strategy change, namely, because there are no dissipative forces, conservation of energy is warrant we consider the Hamiltonian of the system, namely, H = T + V ≡ Energy, this yields,

H ( θ , θ ˙ ) ≡ E = 1 2 m 1 l 2 [ − 1 4 m 1 m 1 + m 2 sin 2 ( θ ) + 1 3 ] θ ˙ 2 + m 1 g l 2 sin ( θ ) + m 2 g a 2 (7)

Equating this to the value of the energy of the system at rest, E₀, gives,

1 2 m 1 l 2 [ − 1 4 m 1 m 1 + m 2 sin 2 ( θ ) + 1 3 ] θ ˙ 2 + m 1 g l 2 sin ( θ ) + m 2 g a 2 = E 0 ≡ m 1 g l 2 ( a l ) + m 2 g a 2 (8)

where, (a/ℓ), the first term of the RHS equals sin(θ₀) = a/ℓ, is the initial inclination angle of the rod w/horizontal.

Equation (8) simplifies,

θ ˙ 2 = 4 g l a l − sin ( θ ) 4 3 − ( 1 1 + m 2 m 1 ) sin 2 ( θ ) (9)

Equation (9) explicitly shows the dependence of the angular velocity, θ ˙ to the masses and the relevant dimensions. Rearranging Equation (9) its integration gives,

t = 1 2 l g ∫ 0 arcsin ( a l ) 4 3 − 1 1 + n sin 2 ( θ ) a l − sin ( θ ) d θ (10)

where we set n = (m₂/m₁). The lower limit of integration indicates the runtime of the rod to the table top. Our primary objective is to evaluate the separation runtime; hence, the lower limit of the integration should be the separation angle. This angle comes about from dynamics aspect of the problem.

The driving force acting on the cube, F₁₂, accelerates the cube according to, F 12 = m 2 x ¨ 2 . At separation instant contact force, F₁₂, vanishes result an equation conducive to the separation angle. First we calculate the needed acceleration. Applying the first equation of set Equation (2) and equation (5) gives,

x ¨ 2 = − l 2 ( 1 + n ) [ θ ¨ sin ( θ ) + θ ˙ 2 cos ( θ ) ] (11)

On the other hand differentiating Equation (9) w/time gives,

θ ¨ = g 2 l { − 4 cos ( θ ) 4 3 − 1 1 + n sin 2 ( θ ) + 1 1 + n ( 1 − 4 sin ( θ ) ) sin ( 2 θ ) [ 4 3 − 1 1 + n sin 2 ( θ ) ] 2 } (12)

And finally substituting Equations (9) and (12) in Equation (11) and simplification yields the needed equation. The output is a lengthy fraction, its numerator is,

− 6 g l ( 1 + n ) cos ( θ ) [ − 8 a l ( 1 + n ) + 12 ( 1 + n ) sin ( θ ) − 3 sin 3 ( θ ) ] = numerator (13)

Setting Equation (13) zero gives the needed equation whose root is the separation angle, namely,

3 sin 3 ( θ ) − 12 ( 1 + n ) sin ( θ ) + 8 a l ( 1 + n ) = 0 (14)

As expected terms such as, n ≡ m₂/m₁ and a/ℓ control the value of the separation angle, θ_sep, yielding the needed runtime,

t = 1 2 l g ∫ θ s e p arcsin ( a l ) 4 3 − 1 1 + n sin 2 ( θ ) a l − sin ( θ ) d θ (15)

It is important emphasizing that entire calculation presented in this section is carried out applying one of the most powerful Computer Algebra System, Mathematica [1] . Equation (14) is generic; its root(s) is the separation angle. The input parameters of this equation may be any desired values, for instance for a practical case we set n ≡ m₂/m₁ = 3, and a/ℓ = 1/4, making the combination a/ℓ (1 + n) = 1, simplifies the equation. With this choice of parameters it yields a unique solution, namely, θ = 9.61˚. In other words the rod from its initial position θ(t = 0) = 14.47˚ falls only 4.86˚ and separates. Equation (14) is a cubic trigonometric equation, its cubic term, 3sin³(θ), can be dropped reducing it to a linear trig equation, 6sin(θ) − 1 = 0, with a root of, θ = 9.59˚. This neat approach underestimating the exact solution only by 0.2%!

Utilizing the separation angle and applying Equation (15) we evaluate the associated runtime by numeric integration. For a 1.0 m rod this gives, t_runtime = 0.108 s. This is compared to the runtime when we assume the rod falls all the way to the table, t = 0.187 s; Equation (10). The latter is compared to the runtime

of a point-like mass at the cm of the rod, t = a g = 0.159   s . As expected the

aforementioned elapsed time overvalues the latter, this is because the rotating rod requires a longer time.

“Sliding Ladder” is a classic physics problem. Its solution is readily available. For instance, the author referred to reference [4] . This particular reference is cited because its free animated solution is also available as well [5] . Our investigation is a modified version of the “Sliding Ladder” problem. Here we replace the stationary wall with a massive mobile cube.

The author searched for the analysis of the proposed problem in vain ending up crafting its solution in this note. Applying Mathematica, the equation of motion, Equation (17), is solved numerically. Utilizing the solution, we determined time-dependent quantities such as runtime, contact force, and trajectory of the tip of the rod. We have also determined additional information such as separation angle of the rod and the cube. Our investigation embodies long-hand analytic formation conducive to parameter-dependent equations including: side length of the cube and its mass, length of the rod and its mass. Applying Mathematica numeric solver and its powerful graphic features as an example, we presented variety of output for a specific set of parameters. Presented information in this report challenges the interested reader to explore new features applying a different set of parameters. Reference [6] embodies Mathematica codes and information needed to produce the graphs embedded in this article.

Sarafian, H. (2018) Kinematics of a Holonomic Constraint Rod & Cube System. World Journal of Mechanics, 8, 227-235. https://doi.org/10.4236/wjm.2018.86018