Column buckling analysis

Last time, we looked at the solutions for six elementary column buckling problems. In this post, we’ll outline the method for getting those solutions and give a couple of detailed examples.

In all of these problems, the column is assumed to be straight and prismatic, which means that the cross-section doesn’t change along its length. The minimum moment of inertia of the cross-section is I, and the column buckles in the direction associated with I.1 The column has a compressive load P acting at either end, and the load is aligned with the centerline of the column. The material is linearly elastic with a Young’s modulus of E, and there is no yielding of the material at any point.

With these assumptions the governing differential equation of the column is

EId 4ydx 4+Pd 2ydx 2=0

where x is the coordinate along the centerline of the column and y is the (small) lateral deflection of the column when it buckles. This is a fourth-order linear differential equation with constant coefficients, which is great for us because it’s one of the few types of differential equation that can be solved analytically.

Here’s the general form of the solution,

y=Asinλx+Bcosλx+Cx+D

where λ is a constant that combines the coefficients,

λ=PEI

and A, B, C, and D are the four constants of integration that appear when solving a fourth-order differential equation. We’ll use the boundary conditions (how the column is supported at either end) to solve for these.

Simple-simple

The column with simple supports at either end is the fundamental problem of column buckling, and it’s pretty easy to work out its solution.

Simple-simple column shown horizontally

(I’m showing the column horizontally to take up less room on the screen. Because we’re not considering the weight of the column itself, the orientation doesn’t matter.)

We’ll start with the conditions at x=0. Here, the column is prevented from moving laterally, so

y(0)=B+D=0

or

D=B

The column is free to rotate at this end, which means there’s no bending moment. Since the bending moment is proportional to the second derivative of the deflection,

y (0)=λ 2B=0

where I’m using the prime notation for derivatives. This means B=0 and therefore D=0, which simplifies things nicely.

Now we look at the other end, x=L, where again the lateral defection is zero,

y(L)=AsinλL+CL=0

and the bending moment is zero,

y (L)=Aλ 2sinλL=0

Since λ isn’t zero, the second equation means

AsinλL=0

and combining that with the first equation tells us that C=0.

So we’ve eliminated all the constants of integration except A. The equation above means either

A=0

or

sinλL=0

The first of these is the solution when the column is straight and y=0 throughout its length. It is a solution, but it’s trivial—we know the column can be straight if the load is small.2 We’re interested in the nontrivial solution for which the column isn’t straight. That’s the buckling solution.

So when does sinλL=0? When λL=nπ, for any integer n. The buckling load will be when P is the smallest positive value that satisfies this equation. In other words,

(λL) cr=P crEIL=π

or

P cr=π 2EIL 2

which is the answer we saw in the last post for the simple-simple column.

What about A? Well, A is undetermined. Remember at the top of the post where I snuck the parenthetical “(small)” into the description of y and didn’t explain it? Here’s where that comes in. Buckling deflections really aren’t small, and to solve for them properly we’d have to start with a nonlinear differential equation. But that doesn’t mean what we’ve done is useless. The value of P cr is correct, and what we’ve done tells us that the deflected shape of the buckled column,

y=Asinλx=AsinπxL

is one-half of a sine wave. This is useful because it helps us visualize the buckling.

Fixed-simple

This is the problem that inspired the previous post. It won’t take as long because we can use some of the simple-simple results to speed things up.

Fixed-simple column shown horizontally

Again, we start with the general solution,

y=Asinλx+Bcosλx+Cx+D

The boundary conditions at x=0 are the same as before,3 so we know that

B=D=0

The boundary conditions at x=L are partly the same and partly different. The deflection is zero,

y(L)=AsinλL+CL=0

and the slope is zero,

y (L)=AλcosλL+C=0

The slope condition means

C=AλcosλL

Substituting this into the deflection boundary condition gives

AsinλLAλLcosλL=0

or

A(sinλLλLcosλL)=0

Once again, there’s the trivial solution, A=0, and the nontrivial solution,

sinλLλLcosλL=0

Divide this through by cosλL and move one term to the other side of the equation to get

tanλL=λL

As we saw in the previous post, the smallest positive solution is

λL=4.4934

and therefore

P cr=4.4934 2EIL 2=20.19EIL 2

Once again, we can’t solve for A, but we can see that the buckled shape is

y=A(sin4.49xL4.49xLcos4.49xL)

which isn’t any sort of standard trig function, but it is what I plotted for the fourth column in this image:

Buckled column shapes


All the other column problems in the image above can be solved the same way. The boundary conditions are different, but the solution process is the same. Notice that all the columns other than the fixed-simple buckle into shapes that are portions of sine waves. As we showed above, the simple-simple column is half of a sine wave. As for the others:

The effective length is based on the number of half sine waves the buckled shape has.

This simple visualization of the buckled shapes allows you to work out the buckling loads for all four of these columns based on the simple-simple column. The only other thing you have to remember is 4.49.


  1. If this sentence is confusing, this post may help. 

  2. I’m not going to get into determining the stability of the straight solution for large loads. That’s beyond this post. 

  3. Because I cleverly put the simple end at the left. You learn a few things after 40 years.