Is it possible to remove the last block from a the lowest layer of a Jenga tower?

4,681

Let's take this in steps.

First - assume there is zero friction between the bottom block and the tower, and between the bottom block and the ground. If I move the block "infinitely quickly" the tower will have had no time to tilt at all, and it will drop vertically (and remain upright).

This tells me there are two things to consider:

  1. the time take to move (if I move slowly, the tower will fall)
  2. the friction between block and tower (with enough friction, I will drag the tower along)

Now I know what equations I have to set up.

First, let's look at the "move slowly, no friction" case. If the tower has height $H$, width $w$ and uniform mass distribution, then it will fall increasingly quickly when we support it off center. I wrote an answer about balancing a pencil a while ago that includes the calculation of the equation of motion of the inverted pendulum. This is complicated by the assumption that the support moves - so let's simplify. If we supported the tower "on the edge" how long would it take to tip over?

In part this depends on the height of the bottom block: for a very thin block, the tower would tilt, then come to rest at a slight angle. For a bigger block, it will tilt until the corner hits the ground, and do so with sufficient speed to continue toppling over. The following diagram shows what we're thinking about:

enter image description here

There is going to be a certain height of block $h$ below which the tower will not topple. Let's start finding that limit...

For mass $m$ initially displaced $w/2$, the total torque is given by

$$\alpha = \tan^{-1}\frac{w}{H}\\ D = \sqrt{H^2 + w^2}\\ \Gamma = m\cdot g \cdot D \sin (\alpha + \theta)$$

And the angular momentum of the rod about the pivot is roughly (for H>>w) $\frac13 m D^2$.

Now the easiest way to calculate whether the tower will continue to topple is by using an energy argument. The tower will topple if it has any speed left when the center of mass is above the corner - that is, when $\theta = \alpha$. If the center of mass is below the initial position when that happens, the tower will fall.

Initially the c of m is at $H/2 + h$; after the corner drops, the maximum height of the center of mass is $\frac{D}{2} = \frac12 \sqrt{H^2 + w^2}$. In other words, the tower is stable as long as

$$H + 2h < \sqrt{H^2 + w^2}$$

if we assume $w$ and $h$ are both much less than $H$, we can rewrite this by expansion:

$$H\left(1 + \frac{2h}{H}\right) \lt H\left(1 + \frac12\left(\frac{w}{H}\right)^2\right)\\ h \lt \frac{w^2}{4H}$$

Note that this analysis was done with the assumption that the block is moved all the way to the edge of the tower: while that would cause the tower to initially start tilting most rapidly, it does restrict the angle it reaches when the corner comes to rest. If we take the other extreme where the block is moved "just off center", then the initial angle of tilt when the corner first hits the ground is $\theta = \tan^{-1}\frac{w}{2H}$. Obviously, if that angle is greater than the internal angle $\alpha$, the tower will continue to fall - but assuming it is not, then the same analysis as above applies. This puts a second limit on $h$:

$\frac{h}{2w} \lt \frac{w}{H}\ h \lt \frac{2 w^2}{H}$$

As you can see, this is in fact a less severe limitation - the chance that the tower will keep moving when it lands sets a tighter limit on how large the block can be.

Now I think we can use this "limiting block height" as a way to estimate the time we get to move a larger block out of the way: if the block is removed before the corner hits the ground, the tower will keep moving at the same rate of rotation as before, while at the same time dropping to the ground. The equations now start to get really messy - especially if we add friction. I think it's time to go to a numerical solution...

The following Python code sample is a very approximate way to approach this. For the particular parameters I chose, you need the block to be pulled out with a velocity of at least 0.36 m/s - not particularly fast. Of course I make all kinds of assumptions here, and I suspect that the physics I modeled was quite sloppy. But you might use this as a starting point for a proper simulation. Draw diagrams, set the equations of motion up properly (as functions instead of inline), and use a better integration scheme (leapfrog, Runge-Kutta) than I used. I am sure that there is at least one error in the calculation - I am not properly accounting for lateral motion of the center of gravity during the pull (force of friction should move entire tower sideways). This will change the solution for very slow velocities - the tower should just move with the bottom block!

So - plenty to do, but I can't spend more time on this now. Thought you might appreciate this nudge in the direction of a solution - but please don't treat this as "correct final answer".

import math
import matplotlib.pyplot as plt
import numpy as np

# very simplistic jenga simulation
H = 0.3  # m
w = 0.05 # m
h = 0.01 # m
D = math.sqrt(H*H+w*w) # diagonal
m = 1    # kg - just to make equations work... drops out everywhere
g = 9.8  # m/s/s
I = m*D*D/3.0

# initial stability criterion:
if (h < w*w/(4*H)):
    print 'stable'
else:
    print 'falls over'

plt.close('all')

#numerical solution - including friction. Loop over range of velocities
for v in np.arange(0.3,0.4,0.01):
 mu = 0.1 # set=0 for frictionless...
 dt = 0.0001
 theta = 0
 omega = 0
 alpha = math.atan(w/H)
 x = 0
 y = h + 0.5*H # center of mass height
 yc = h # corner height
 t=[0]
 a = [0]
 corner=[h]
 vy = 0
 while ((x < w/2) & (yc > 0)):
    beta = theta + alpha
    # torque due to gravity and friction
    # not sure about that second term...
    gamma = 0.5 * m * g * D * math.sin(beta) + m*g*mu*H/2 * (I/(I+m*H*H/4)) 
    dOmega = dt * gamma / I
    omega = omega + dOmega
    dy = omega * math.sqrt(x*x + 0.25*H*H) * math.cos(beta)
    y = y - dy # center of mass height
    yc = h - (x + w/2)* math.sin(theta) # corner height
    theta = theta + omega * dt
    a.append(theta)
    t.append(t[-1]+dt)
    corner.append(yc)
    x = x + v * dt
 # part 2 = free fall
 while(yc > 0):
    theta = theta + omega * dt
    vy = vy + g * dt
    dy = vy * dt
    yc = yc - vy*dt - dt * omega * D/2 * math.sin(alpha-theta)
    t.append(t[-1]+dt)
    corner.append(yc)
    a.append(theta)

 # part 3: does it keep falling:
 while((omega > 0) & (theta < math.pi/2)):
    gamma = m*g*D/2*math.sin(theta-alpha)
    dOmega = dt * gamma / I
    omega = omega + dOmega
    theta = theta + omega*dt
    a.append(theta)
    t.append(t[-1]+dt)
    corner.append(0)

 plt.figure()
 plt.plot(t, a)
 plt.plot(t, corner)
 plt.title('v = %.2f m/s'%v)
 plt.show()

Plotting the height of the corner and the angle of the tower as a function of time, you get different regimes. If you move very slowly, the tower falls over completely (blue line = angle; green line = height of corner above ground):

enter image description here

If you pull very quickly (above 0.37 m/s in this instance - but that depends a lot on the chosen parameters) it will tip a little, then recover (simulation stops when the tower starts turning in the other direction - i.e. when it turned out to be stable):

enter image description here

Finally, in this simulation the tower just goes unstable at 0.36 m/s pull speed: you can see that it teeters, then falls (note the longer time scale):

enter image description here

Share:
4,681

Related videos on Youtube

tlewis3348
Author by

tlewis3348

Updated on June 13, 2020

Comments

  • tlewis3348
    tlewis3348 over 3 years

    A while back, my brother and I had a discussion about the possibility of removing the last block from the lowest layer of a Jenga tower with the layer immediately above it having only one block in it as well (see illustration below). His theory was that if you were fast enough, you could pull it out like a magician pulls a table cloth out from under a table full of china.

                                                                  Illustration

    I protested against that idea by arguing that the tower is so unstable that in order for the block (made out of maple) to be pulled out fast enough, it would end up being damaged.

    I know enough about physics to know that it would be possible to calculate the minimum speed at which something would have to be moving to smack the block hard enough to knock it out from under the tower without knocking the tower over, and I know that if that speed is known, calculating whether the block would get damaged would be relatively simple. However, I don't know enough about physics to actually be able to know how to do the calculation.

    Is there anyone here that would be able to give an explanation of how that speed could be calculated?

    Edit:

    The dimensions of the standard Jenga block are listed here as 1.5 × 2.5 × 7.5 cm, and they are made of maple (I believe it's hard maple, if you want to be precise). Additionally, as the image above shows, I am assuming the shortest tower possible. That is, 17 layers of three blocks each on top of two layers composed of one block each with the extra block resting on top.

    • innisfree
      innisfree over 8 years
      you might consider stability by looking at the potential energy as a function of the angle of tilt, $V(\theta)$, to see whether it is stable wrt small perturbations
    • innisfree
      innisfree over 8 years
      what's your question got to do, specifically, with Lagranign-formalism? I'm removing that tag for now
    • tlewis3348
      tlewis3348 over 8 years
      @innisfree I seem to remember using Lagrangian methods to solve problems similar to this when I was in school. Specifically, I can remember solving a problem where someone is pushing a large cart that hits a rock, and we had to figure out whether the cart would tip over or not (or something like that). I haven't been able to find a copy of that homework, but I remember doing it, and I'm pretty sure we were studying Lagrangian methods when it was assigned.
    • Shin Kim
      Shin Kim over 8 years
      It depends on the height of the tower because the gravitational force exerted on the lowest block is proportional to the mass and the mass is proportional to the height in this case. So it can damage the lowest block due to the strong frictional force if it is high enough, and the lowest block can be pulled out without damage if it's low enough.
    • tlewis3348
      tlewis3348 over 8 years
      @ShinKim According to this, the dimensions of Jenga blocks are 1.5 × 2.5 × 7.5 cm, and they are made of Maple. Also, I'm assuming the shortest tower possible to meet the stipulations mentioned in the question. This tower is shown in the attached image in the post.
  • innisfree
    innisfree over 8 years
    how did you make that professional looking illustration? tikz?
  • Admin
    Admin over 8 years
    Really wonderful +1.
  • tlewis3348
    tlewis3348 over 8 years
    @Floris Thanks! Just a little bit of additional information: According to this, the dimensions of Jenga blocks are 1.5 × 2.5 × 7.5 cm, and are made of Maple. Also, I'm assuming the shortest tower possible to meet the stipulations mentioned in the question. This tower is shown in the attached image in the post.
  • Floris
    Floris over 8 years
    @innisfree - I draw (most of) my diagrams in Powerpoint and take a screen shot (Mac: cmd-shift-4); for the graphs, I use matplotlib.pyplot (at other times I use Excel).
  • Floris
    Floris over 8 years
    @tlewis3348 - thank you. If I have time I will update with those dimensions. Might even try to improve the equations of motion to take account of the horizontal motion... this is fun but time consuming.
  • tlewis3348
    tlewis3348 over 8 years
    @Floris "this is fun but time consuming" <-- I most certainly agree and understand. Ultimately, as you probably guessed from my OP, it's not extremely pressing for me to get the solution, but rather is just an interesting problem that I've been thinking about for a little while.
  • innisfree
    innisfree over 8 years
    that's interesting - powerpoint has come on a lot since I last used it >~ 10 years ago!