Kinematics of a Holonomic Constraint Rod & Cube System

A loose rod of mass m1 and length l leans against one of the faces of a cube of mass m2 and side length a. The assembly is placed on a horizontal table with one end of the rod touching the table and its other end leaning against the edge of the cube. We set the rod and the center of mass of the cube on the same vertical plane, and then we release the assembly from the rest. For frictionless contacts, we calculate the separation runtime of the rod from the cube as a function of m2/m1 and a/l. This entails forming the equation describing the motion of the system. The equation of motion is analytically unsolvable nonlinear differential equation. Applying a Computer Algebra System, specifically Mathematica [1] [2], we solve the equation numerically. Utilizing the solution, in addition to evaluating the separation runtime, we quantify a list of dynamic quantities, such as the time-dependent interface forces, and, geometric quantities, such as the trajectory of the loose end of the rod. A robust Mathematica code addresses the “what if” scenarios.


Introduction
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.As shown in Figure 1, one end of a uniform rod of length ℓ and mass m 1 touches the horizontal surface at x 1 , its other end leans against the left vertical face of the cube at C(x 2 , y 2 ).The cube has a mass of m 2 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 21 , are also shown.Figure 1.A rod of mass m 1 and length ℓ leans against a cube of mass m 2 and side length a.

Physics of the Problem and Its Solution
The assembly is on a frictionless horizontal surface.

H. Sarafian
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 21 .As depicted, the same value of the force F 21 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 1 .While the rod falls its "circulation radius" shortens hence separation is warranted.In a different scenario if the rod was hinged at x 1 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  and its potential energy with respect to the table top is . 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, where ω equals θ  is the angular velocity of the rod about the cm;

V m gy =
; 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, ( )  .Utilizing the last equation and Equation (2), Lagrangian of the system, e.g.

(
) Equation ( 3) reveals the Lagrangian is independent of x 1 , meaning x 1 is a cyclic variable; = warrants the conservation of the linear momentum of the system along the horizontal direction.Explicitly, momentum of the system 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, ( ) Therefore substituting Equation ( 6) to Euler-Lagrange equation, 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, ( ) Equating this to the value of the energy of the system at rest, E 0 , gives, ( ) ( ) where, (a/ℓ), the first term of the RHS equals sin(θ 0 ) = a/ℓ, is the initial inclination angle of the rod w/horizontal.
Equation ( 9) explicitly shows the dependence of the angular velocity, θ  to the masses and the relevant dimensions.Rearranging Equation (9) its integration gives, ( ) ( ) where we set n = (m 2 /m 1 ).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 12 , accelerates the cube according to, 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, Setting Equation ( 13) zero gives the needed equation whose root is the separation angle, namely, ( ) ( ) ( ) ( ) As expected terms such as, n ≡ m 2 /m 1 and a/ℓ control the value of the separation angle, θ sep , yielding the needed runtime, ( ) ( )

Rod-Cube Contact Force
As the rod slides its inclination angle, θ, with respect to the table top decreases, reduces the contact force between the rod and the cube.To quantify the contact force at any given time, we write, 12 2 2 F m x =  .Acceleration of the cube, 2 x  , uti- lizing Equation ( 11), ( 12) and ( 9), yields, This is the formal solution, however, because θ = θ(t) there is a need to solve the equation of motion for θ(t).The equation of motion is actually Equation (9).
As demonstrated, by arranging Equation ( 17) we calculate the runtime; see Equation ( 10) and (15).Reversing these equations to obtain θ(t) is not possible; one needs to solve Equation (17) directly, this can't be done analytically.In the next section for a chosen set of parameters we present its numeric solution.

Trajectory of the Tip of the Rod
As the contact point of the rod with the table top slides to the left its other end that leans against the mobile cube before breaking away traces a certain curve.It is our interest to plot its trajectory.We offer two solutions.First, a formal analytic solution, and second a straight forward plotting routine.Both solutions hinge upon the know solution of Equation (17).

Analytic solution:
Integrating both sides of Equation (5) w/time yields, ( ) ( ) here we set the initial contact point of the rod with the table at the origin yielding the value of β = cos[arcsin(a/ℓ)].Substituting Equation (18) in the first equation of set of Equation ( 1) after simplification, gives, on the other hand Figure 1 shows, ( ) Squaring both sides of Equation ( 19) and (20) and adding them results, ( ) ( ) where θ(t) is the solution of the equation of motion; Equation (17).For a set of chosen parameters n and ℓ and utilizing the solution of Equation (17) for the same set of parameters along with a parametric plot of Equation ( 22) yields the seek for trajectory.This is shown in the next section.

Results
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 2 /m 1 = 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 3 (θ), 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, 0.159 s a t g = = . As expected the aforementioned elapsed time overvalues the latter, this is because the rotating rod requires a longer time.

Figure 1
Figure 1 depicts the schematic of the problem on hand.
 .At separation instant contact force, F 12 , 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,

Figure 2
Figure 2 depicts the time dependent variation of the inclination angle.The rod begins sliding with the initial inclination angle of arcsin(1/4) = 14.4˚, in 0.108 s at 9.61˚ it separates from the cube.

Figure 3
Figure3shows the maximum value of the contact force between the rod and the cube occurs at the beginning where the inclination angle is at maximum.As the rod falls and rotates the contact force weakens.By the time the rod separates from the cube its value plunges zero.Comparing Figure2and Figure3 reveals

Figure 4 .
Figure 4.The solid line is the trajectory of the tip of the rod as it falls to the table.
Equation (21) deceptively appears as an ellipse; it is not.Coordinates of {x 2 , y 2 } are time-dependent, they are implicit functions of θ(t).Equation (21) can be utilized tracing the tip of the rod; this is left to the interest of the reader.Alternatively, utilizing Equations (19) and (20) are arranged as,