26

I have a paper which lies on a flat surface. The paper is fixed on one side and the opposite side can slide in the direction of the opposing side. As side end slides toward the other, a "bump" forms. I want to know what the solution is to the shape of this bump. There must be some standard solution for this case.

EDIT: Luboš Motl did provide a really nice answer for the case when the paper is assumed to lie flat against the surface at the ends (also with the assumption that the bump is small). I also found this interesting paper on a similar topic.

Krumelur
  • 363

3 Answers3

16

A very natural and interesting question.

If the paper stays nearly flat, we may use the linearized approximation. I will neglect the gravitational potential energy because it's probably much smaller than the bending energy (note that the shape of the bent paper is almost the same when the desk is vertical as when it is horizontal).

The paper doesn't want to bend much, so the energy contains a term that punishes the second derivatives of $y$ (the curvature) $$E = K \int_0^L {\rm d}x\,(y'')^2$$ where $y'\equiv dy/dx$, and so on. We want to keep the length of the paper fixed - it's prescribed by the way how you attach it at the end points. In the linear approximation, the length is a linear function of $$D = \int_0^L {\rm d} x\, (y')^2 $$ which may be seen, if you need it, by a Taylor expansion of the exact expression for the length of the graph, $\int (1+y^{\prime 2})^{1/2}{\rm d}x$. Add it with a Lagrange multiplier, requiring $$\delta (E+\lambda D) = 0.$$ I think that the resulting equations say $$ y'''' = \frac{\lambda}{K} y''.$$ Note that in the Euler-Lagrange equations, all the primes "clump". If you write $Y=y''$, it says that the second derivative of $Y$ is proportional to $Y$ itself. The physically relevant solution is $$ Y = Y_0 \cos (2\pi x n / L) $$ where only $n=1$ is possible if there is a table underneath the paper. Consequently, $$ y = y_0 [1-\cos (2\pi x / L)],$$ too. I added the term $1$ as an integration constant (while the other is zero) to guarantee that $y=0$ and $y'=0$ at the boundaries.

So what you get in practice is one wave of a cosine, from one minimum to the next one.

A more detailed discussion would be needed to eliminate the other cosine-like function with the same frequency, the sine, as well as linear functions, and to discuss whether the exponentially increasing/decreasing solutions may ever be relevant. Well, it wouldn't be too difficult. A more general solution for the same (cosine-like) sign of $\lambda$ would be $$ y = y_0 [1-\cos (\phi_0+ 2\pi x N / L)] + Ax^2+B,$$ and the four conditions $y=0$, $y'=0$ at $x=0,L$ as well as the fixed length of the paper would imply $N=1$, $\phi_0=0$, $A=0$, $B=0$, as well as the right normalization $y_0$. As suggested above, the conditions would have other solutions where $N=2,3,4\dots$ but these solutions wouldn't satisfy $y\geq 0$ for $0\leq x \leq L$, so the paper couldn't arrange itself above the table. If there were a hole in the table, these higher harmonics (several periods of the wave containing) solutions would be possible - but I guess that they would be unstable.

I am not sure whether I would be able to solve it analytically if the angle $y'$ were not infinitesimal. It seems clear that the exact solution for significant angles is not just a cosine: the bent paper tends to resemble a balloon - for which $y(x)$ is not even single-valued - when there's too much redundant paper in the middle.

On the other hand, if you sharply bend the paper at $x=0,L$ so that $y'$ may be arbitrary at these two points, the paper in between will be bent as an arc of a circle - see another answer below - and this result is accurate even if $y'$ is of order one. Note that in the linearized approximation, the arc gives $y$ being a quadratic function of $x$ so that $y''=0$ still solves the fourth-order equation above.

Luboš Motl
  • 179,018
  • Hey, great to see that my intuition still works well 20 years later... ;-) – ysap May 16 '11 at 16:43
  • 1
    Dear ysap, excellent to hear about it. But is it really intuition, or just plain good old memory that remembers someone else's intuition or, more likely, calculation? ;-) – Luboš Motl May 16 '11 at 16:44
  • Well, I guess it's both. Honestly, this specific course was one of my least favorites (to say the least) while studying Aerospace Eng. What I did vaguely remembered is that the this problem is modeled by a simple diff. eqn., and that the shape of the buckled beam is an eigenfunction of that eqn., which is harmonic. Testing it with an actual sheet of paper, I assured myself that it is a full cycle of the harmony (due to the table B.C., as you described), and not just half of it. – ysap May 16 '11 at 16:52
  • @Lubos Pardon my ignorance, but I did not understand the statement "The paper doesn't want to bend much, so the energy contains a term that punishes the second derivatives of y (the curvature) " – yayu May 16 '11 at 16:58
  • Dear @yayu, the $(y')^2$ term vanishes if the paper is straight. If it is curved at a given point $x$, the term is nonzero. It's quadratic because you may really imagine that there are springs on both sides of the paper, trying to keep is straight, and the $(y')^2$ term generalizes the $kx^2/2$ potential energy of a spring. – Luboš Motl May 16 '11 at 19:09
  • Excellent answer! Thank you. The answer is clearly observable, but I wanted to know why. – Krumelur May 17 '11 at 06:19
  • @Luboš: A general solution of $ y'''' = a y''$ is $ y(x) = \frac{e^{-\sqrt{a} x} (c_1 e^{2 \sqrt{a} x}+c_2)}{a}+c_4 x+c_3$, i.e. not a trigonometric. – Martin Gales May 17 '11 at 07:37
  • Dear @Martin, $a$ in your equation has to be negative, so that its square root is imaginary, so that you have the complex exponential, and it has to be combined to the trigonometric functions - by having $c_1=c_2^*$ - for the result to be real. By the way, the positive $a$ solutions - which are exponential (and which I have surely mentioned in my answer as well) - may be shown to be unstable. They require the opposite sign of $\lambda$ which is translated to the shape's being a local maximum, rather than minimum, of the energy. – Luboš Motl May 17 '11 at 07:56
  • Correction from "by the way": By the way, I have mentioned the exponentially increasing and decreasing solutions in my answer as well and I don't know how to make them realize in practice. The pressure between the end points of the paper should probably have the opposite sign in this case, because of the opposite sign of $\lambda$, and I don't think it's possible. So I think that only the trigonometric plus quadratic solutions are relevant. – Luboš Motl May 17 '11 at 08:02
  • @Luboš - you mentioned that the solutions for N>1 imply that y(x)<0 at some point. I fail to see why. Even with higher frequencies, the cosine is bounded between -1 and 1, so adding the constant 1 guarantees y(x)>=0. Am I missing something? – ysap May 17 '11 at 22:53
  • Neglecting gravity, would it not be easier to use the coordinate system of the paper itself? – Krumelur May 18 '11 at 06:29
2

You gave some boundary contraints, but not enough to solve this problem without some assumptions on the part of the answerer. The missing piece is: For each end; is the end...

1.) constrained to lie flat against the table, or

2.) allowed to pivot up.

If the answer is 1 for either or both ends, the answer is non-trivial.

If the answer is 2 for BOTH ends, and the weight of the paper is non-negligible, then, again the answer is non-trivial.

If the answer is 2 for BOTH ends -AND- the weight of the paper is negligible, then the paper will bend in the arc of an almost perfect circle. With the distace between two ends = paper height, the circle radius is infinity. With the distance between the ends = 0, the radius is paper height divided by (2*Pi); and the arc is the complete circle.


EDIT, based on fearlesscoder's, Lubos' and Ysap's comments:

If the paper is constrained to lie in the plain of the table at the ends (dy/dx = 0), then the answer could indeed be (or approximate) a complete cycle of a cosine wave (with an offset in the vertical axis); as long as the ends are constrained to not approach closer than L/8 to each other, where L is the length of the sheet. Not sure what to do as the ends of the paper get really close to one another: the cosine thingy breaks down in that region.

It would not, in any case (where the ends are constrained to lie flat), resemble a circle, nor a segment of a circle, nor a quadratic nor a Gaussian curve (which is the first thing that popped to mind when I realized exactly what boundary conditions the OP meant).

I would make the assumption that the resistance of the sheet to bending is proportional to the radius of curvature (((d^2)x)/(dy^2) = ((d^2)y)/(dx^2)). I would make the peak of the bump at x=0. I would have to put in a lot of thought getting to the point where I can find out if y=cos(x) is a sol'n to that problem.

Vintage
  • 4,991
  • 2
  • 16
  • 15
  • 1
    I think that the effect of the paper weight can be eliminated (at least, reduced to a great extent) by considering a vertical setup, where the paper is bound to a vertical surface and is pressed from both sides. Other than that, real life solutions are never trivial, but simple, and pretty accurate models do exist frequently. @Lubos' answer develops such a model and solution that show the shape to be harmonic, and not circular. – ysap May 16 '11 at 17:58
  • Dear @Vintage, I agree that if one is allowed to sharply bend the paper at $x=0$ and $x=L$, then the paper in between will make an arc of a circle, or - in my approximation - $y=ax^2+bx+c$. This is a sort of trivial solution. I think that the OP meant 1) that the paper is never sharply bent and it continues as flat plane on both end points. – Luboš Motl May 16 '11 at 20:35
  • 1
    never mind the hand waving, how about proving this mathematically. – Larry Harson May 17 '11 at 11:05
  • @ user2146 - No, thank you. As I said, the sol'n is non-trivial. – Vintage May 17 '11 at 17:06
  • Sorry about the lack of boundary constraints, your 1.) was the one I had in mind, though the other one is interesting as well. – Krumelur May 17 '11 at 17:17
1

I'm quite sure that the idealized shape is one cycle of a cosine. This is due to the differential equations governing the buckling of thin structures (beams).

However, there has been many, many years since I studied this kind of stuff, so I don't have reference at hand.

ysap
  • 149
  • What if you have a shape that does not support bending moments (zero thickness paper)? It will still form a shape due to the conservation of mass, but the EI component will not be there... Interesting! – John Alexiou May 16 '11 at 17:54
  • @ja72 - I am not sure, but I think that the simple model assumes zero thickness. – ysap May 16 '11 at 18:01