# Self-buckling

A column can buckle due to its own weight with no other direct forces acting on it, in a failure mode called self-buckling. In conventional column buckling problems, the self-weight is often neglected since it is assumed to be small when compared to the applied axial loads. However, when the weight of the column is significant compared to its buckling strength, it is important to take self-buckling into account.
Elastic buckling of a "heavy" column i.e., column buckling under its own weight, was first investigated by Greenhill in 1881. He found that a free-standing, vertical column, with density $\rho$, Young's modulus $E$, and cross-sectional area $A$, will buckle under its own weight if its height exceeds a certain critical value:
    $l_{\text{max}}\approx \left(7.8373\,{\frac {EI}{\rho gA}}\right)^{\frac {1}{3}}$
where $g$ is the acceleration due to gravity, and $I$ is the second moment of area of the beam cross section.
One interesting example for the use of the equation was suggested by Greenhill in his paper. He estimated the maximal height of a pine tree, and found it cannot grow over 300 feet (90 m) tall. This length sets the maximum height for trees on earth if we assume the trees to be prismatic and the branches are neglected.

## Mathematical derivation

Suppose a uniform column fixed in a vertical direction at its lowest point, and carried to a height $l$, in which the vertical position becomes unstable and flexure begins. There is a body force $q$ per unit length $q=\rho gA$, where $A$ is the cross-sectional area of the column, $g$ is the acceleration due to gravity and $\rho$ is its mass density.
The column is slightly curved under its own weight, so the curve $w(x)$ describes the deflection of the beam in the $y$ direction at some position $x$. Looking at any point on the column, we can write the moment equilibrium:
    $M=-\int _{0}^{x}q(y-w)\mathrm {d} x$
where the right-hand side of the equation is the moment of the weight of BP about P.
According to Euler–Bernoulli beam theory:
    $M=-EI{\mathrm {d} ^{2}w \over \mathrm {d} x^{2}}$
Where $E$ is the Young's modulus of elasticity of the substance, $I$ is the second moment of area.
Therefore, the differential equation of the central line of BP is:$$$$
    $EI{\mathrm {d} ^{2}w \over \mathrm {d} x^{2}}=q\int _{0}^{x}(y-w)\mathrm {d} x$
Differentiating with respect to x, we get
    ${\begin{aligned}EI{\mathrm {d} ^{3}w \over \mathrm {d} x^{3}}&=q\int _{0}^{x}\left(-{\mathrm {d} w \over \mathrm {d} x}\right)\mathrm {d} x\\&=-qx{\mathrm {d} w \over \mathrm {d} x}\end{aligned}}$
We get that the governing equation is the third order linear differential equation with a variable coefficient. The way to solve the problem is to use new variables $n,z,k$ and $r$:
    $k^{2}={4 \over 9}{q \over EI},\quad r^{2}=x^{3},\quad {\sqrt {x}}z={\mathrm {d} w \over \mathrm {d} x},\quad n^{2}={1 \over 3^{2}}$
Then, the equation transforms to the Bessel equation
    $r^{2}{\mathrm {d} ^{2}z \over \mathrm {d} r^{2}}+r{\mathrm {d} z \over \mathrm {d} r}+\left(k^{2}r^{2}-n^{2}\right)z=0$
The solution of the transformed equation is $z=AJ_{\frac {1}{3}}(kr)+BJ_{-{\frac {1}{3}}}(kr)$
Where $J_{n}$ is the Bessel function of the first kind. Then, the solution of the original equation is:
    ${\mathrm {d} w \over \mathrm {d} x}={\sqrt {x}}\left(AJ_{\frac {1}{3}}\left(kx^{\frac {3}{2}}\right)+BJ_{-{\frac {1}{3}},1}\left(kx^{\frac {3}{2}}\right)\right)$
Now, we will use the boundary conditions:
* No moment at $x=0\rightarrow {\mathrm {d} ^{2}w \over \mathrm {d} x^{2}}(x=0)=0\rightarrow A=0$
* Fixed at $x=l\rightarrow w(l)=0\rightarrow J_{-{\frac {1}{3}}}\left(kl^{\frac {3}{2}}\right)=0$
From the second boundary condition, we get that the critical length in which a vertical column will buckle under its own weight is:
    $l_{\text{max}}=\left({\frac {j_{{\frac {1}{3}},1}}{k}}\right)^{\frac {2}{3}}=\left({\frac {9{j_{{\frac {1}{3}},1}}^{2}}{4}}{\frac {EI}{q}}\right)^{\frac {1}{3}}$
Using $j_{{\frac {1}{3}},1}\approx 1.86635$ , the first zero of the Bessel function of the first kind of order $-{\frac {1}{3}}$, $l_{\text{max}}$ can be approximated to:
    $l_{\text{max}}\approx \left(7.8373\,{\frac {EI}{\rho gA}}\right)^{\frac {1}{3}}$

## Euler's mistake

The column under its own weight was considered by Euler in three famous papers (1778a, 1778b, 1778c). In his first paper, Euler (1778a) concluded that the column simply supported under its own weight would never lose its stability. In his second paper on this topic Euler (1778b) described his previous result as paradoxical and suspicious (see Panovko and Gubanova (1965); Nicolai, (1955); Todhunter and Pierson (1866) on this topic). In the next, third in series, paper, Euler (1778c) found that he had made a conceptual mistake and the “infinite buckling load” conclusion was proved to be wrong. Unfortunately, however, he made a numerical mistake and instead of the first eigenvalue, he calculated a second one. Correct solutions were derived by Dinnik (1912), 132 years later, as well as Willers (1941), Engelhardt (1954) and Frich-Fay (1966). Numerical solution with arbitrary accuracy was given by Eisenberger (1991).
222 years after Euler's mistake in 1778, Elishakoff revisited this problem and derived closed-form solutions for self-buckling problems.