1 Graded Problems
Total Page:16
File Type:pdf, Size:1020Kb
PHY 5246: Theoretical Dynamics, Fall 2015 Assignment # 9, Solutions 1 Graded Problems Problem 1 (1.a) In order to find the equation of motion of the triangle, we need to write the Lagrangian, with generalized coor- dinate θ . The potential energy is going to be based on y { } the location of the center of mass, and we find that via l/2 l/2 ρ l/2 0 − x yCM (θ =0) = 2 dx dyy m 0 √3(l/2 x) Z Z− − 2 m l/2 1 l 2 = dx 3 x CM 1 √3 −m 0 2 2 − ! 2 l 2 l Z 2 8 3 1 l l/2 θ = x √3l2 2 3 2 − ! 0 3 4 l l √3l = = , − 2 −√3l2 8 −2√3 where the density is ρ. This is as expected by symmetry considerations and by the fact that the CM has to be at the geometrical center of the triangle, i.e. at one third of it’s height, 1 √3 l l = . −3 2 −2√3 Note the negative sign is because the triangle lies below the x-axis. So, for generalized θ, ρ yCM = cos θ, −2√3 and the potential energy is l U = mgyCM = mg cos θ. − 2√3 To calculate the kinetic energy, 1 T = T (about 0) = I θ˙2, rot 2 3 we need first to find the principal moment of inertia about the axis of rotation, which is an axis perpendicular to the plane of the triangle, through 0. This is the only moment of inertia we need, since ˙ ~ω = θˆe3 and therefore, 1 1 1 T = ~ωT ˆI ~ω = θ˙2I = θ˙2I . 2 · · 2 33 2 3 l/2 0 2 2 I3 = 2ρ dy(x + y ) 0 √3(l/2 x) Z Z− − l/2 0 l/2 0 = 2ρ dxx2 dy +2ρ dx dyy2 0 √3(l/2 x) 0 √3(l/2 x) Z Z− − Z Z− − l/2 l l/2 1 l 3 = 2ρ dxx2√3 x +2ρ 3√3 x Z0 2 − ! Z0 3 2 − ! 1 l l3 1 l4 1 l4 = 2ρ√3 + 3 2 8 − 3 16 4 16! m 1 l4 ml2 = 2√3 = . 1 √3 2 3 16 6 2 2 l Therefore 1 ml2 T = θ˙2, 2 6 ! and the Lagrangian is 1 ml2 l L = T V = θ˙2 + mg cos θ. − 2 6 ! 2√3 The equation of motion for our generalized coordinate is d ∂L ∂L ml2 l =0 θ¨ + mg sin θ =0 dt ∂θ˙ − ∂θ → 6 2√3 g θ¨ + √3 sin θ =0. ⇒ l Under small oscillations about θ = 0 this is the equation for a simple harmonic oscillator with frequency g ω = √3 . r l (1.b) For this situation, the θ = 0 position of the CM is y (θ =0)= l . Therefore CM − √3 y l l yCM = cos θ U = mg cos θ. −√3 ⇒ − √3 We need to calculate the moment of inertia I3 with re- l/2 l/2 spect to the new center of rotation (and axis through − x it). l/2 √3x − 2 2 I3 = 2ρ dx dx(x + y ) 0 √3/2l Z Z− l/2 √3x l/2 √3x 2 − − 2 CM = 2ρ dxx 3 dy +2ρ dx 3 dyy θ 0 √ l 0 √ l Z Z− 2 Z Z− 2 l/2 l l/2 1 l3 = 2ρ√3 dxx2 x + + dx 3 x3 + √3 0 2 0 3 8 l "Z − ! Z − !# − 2 1 l4 1 1 l4 l4 = 2ρ√3 + l4 + −4 16 3 − 4 16 16! m l4 3+2+6 5 = 2√3 − = ml2. 1 √3 2 16 6 12 2 2 l Thus the kinetic energy and Lagrangian are 1 5 T = ml2 θ˙2 2 12 1 5 l L = T V = ml2 θ˙2 + mg cos θ, − 2 12 √3 and the equation of motion is d ∂L ∂L 5 l =0 ml2θ¨ + mg sin θ =0 dt ∂θ˙ − ∂θ → 12 √3 g 12 1 θ¨ + sin θ =0. ⇒ l 5 √3 Thus the frequency of the small oscillation is 12 g ω = . s5√3 l Problem 2 (a) y y y = l cos θ CM √2 θ l The cube rotates about the side perpendicular to the plane of the figure, through 0. We are therefore referring to the case in which the frame has axes parallel to the side of the cube and origin at one corner. In figure (1), the cube is rotating about the back side, along the y-axis. The 2 2 corresponding moment of inertia is I22 = 3 ml (also see p. 409-410 in Goldstein). We can find the angular velocity of the cube when one face strikes the plane, by imposing conservation of energy: l Ei = Ui = mg √2 l 1 2 2 2 Ef = Uf + Tf = mg + ml ωf 2 2 3 l l 1 2 2 Ei = Ef mg = mg + ml ω . ⇒ √2 2 3 f 1 1 1 2 3 g g = lωf ωf = (√2 1) . ⇒ √2 − 2! 3 ⇒ s2 − l z O y x Figure 1: (b) If the cube can slide, there is no more fixed point. So, we have to calculate the kinetic energy as CM CM T = Ttrans + Trot 1 2 1 1 2 2 = mvCM + ml ωf , 2 2 6 where: l l vCM =y ˙CM = sin θθ˙ = θ,˙ −√2 −2 circ and when θ = 45 , the cube hits the plane. With ω(θ = 45◦)= ωf , we have 1 l2 1 5 T = m ω2 + ml2ω2 = lω2, 2 4 f 12 f 24 f Or 1 1 5 2 12 g g = lωf ωf = (√2 1) . √2 − 2! 24 ⇒ s 5 − l We could also have calculated the kinetic energy as CM CM T = Ttrans + Trot 1 2 1 CM 2 1 1 1 2 2 = mvCM + I22 ωf = + l ωf 2 2 2 2 6 1 = l2ω2 = T origin. 3 f rot 2 ˙ l CM ml Now use v = Rθ to find vCM = √2 ωf , and on p. 422 of Goldstein we can find I22 = 6 . We thus obtain the same answer as above. Problem 3 (Goldstein 5.17) Take M as the mass of the cone, (x1′ , x2′ , x3′ ) as the principal axes going through the CM, which 3 is at (0, 0, a = 4 h) in the body system. x3′ is the cone axis, x2′ is perpendicular to the plane going through x3′ and the segment OA, and x1′ is orthogonal to the plane of x2′ and x3′ . As stated in the problem, R is the radius, h is the height and α is the half-angle of the cone. Due to the condition of no-slipping, the angular velocity of the cone, ~ω, is directed along the instantaneous axis of rotation which is the line of contact between the surface of the cone and the surface on which the cone is rolling (i.e. along the segment OA, see figure). x3 x′1 O a x2 α θ CM x R ′3 x1 x′2 A ~ω (a) The kinetic energy (T ) can be calculated either using the center of mass as origin and taking as reference frame the system of principal axes described above, or using a system of axes parallel to that one but with origin at the apex of the cone O, which is a fixed point for the body. In the first case, the kinetic energy is given by the sum of two terms: the kinetic energy of the center of mass (i.e. the kinetic energy of a pointlike object of mass M located at the CM of the cone and moving as CM moves) and the kinetic energy of the body as rotating about the center of mass (expressed in terms of the principal moments of inertia and the component of the angular velocity in the principal axes frame or body frame), CM CM T = Ttrans + Trot kinetic energy 1 1 T CM = = Mv2 = Ma2 cos2 αθ˙2. trans of the CM 2 CM 2 ˙ In the above equation we have used that vCM = a cos αθ as it rotates about ˆx3. This implies that ω (angular velocity about OA) and θ˙ are related. Thus we also have vCM = a sin αω. Combining these two relations we find v ω = CM = θ˙ cot α. (1) a sin α The components of ~ω in the body frame are cos2 α ~ω =(ω sin α, 0,ω cos α)=(θ˙ cos α, 0, θ˙ ). sin α This follows from our choice of x2′ and because of the relation (1). Therefore: 1 1 T CM = I ω2 + I ω2 rot 2 1 1 2 3 3 2 I = I = 3 M R2 + h = 1 2 20 4 I = 3MR2. n 3 10 2 4 1 3 2 h 2 ˙2 1 3 2 cos α ˙2 = M R + cos αθ + MR 2 θ 2 20 4 ! 2 10 sin α 2 4 CM CM 1 2 2 ˙2 1 3 2 h 2 ˙2 1 3 2 cos α ˙2 T = Ttrans + Trot = Ma cos αθ + M R + cos αθ + MR 2 θ 2 2 20 4 ! 2 10 sin α 9 3 1 3 = Mh2 cos2 αθ˙2 + Mh2 tan2 α + cos2 αθ˙2 + Mh2 cos2 αθ˙2. 32 40 4 20 3 To get this last line we used a = 4 h, R = h tan α, and the next line uses the trigonometric identity, 2 1 tan α = 2 1. cos α − 9 3 3 3 3 T = Mh2 cos2 αθ˙2 Mh2θ˙2 Mh2 cos2 αθ˙2 + Mh2 cos2 αθ˙2 32 40 − 40 4 20 3 = Mh2θ˙2(1+5cos2 α).