Physics:Euler's critical load

From HandWiki
Short description: Formula to quantify column buckling under a given load
Fig. 1: Critical stress vs slenderness ratio for steel, for E = 200 GPa, yield strength = 240 MPa.

Euler's critical load is the compressive load at which a slender column will suddenly bend or buckle. It is given by the formula:[1]

[math]\displaystyle{ P_{cr} = \frac{\pi^2 EI}{(KL)^2} }[/math]

where

  • [math]\displaystyle{ P_{cr} }[/math], Euler's critical load (longitudinal compression load on column),
  • [math]\displaystyle{ E }[/math], Young's modulus of the column material,
  • [math]\displaystyle{ I }[/math], second moment of area of the cross section of the column (area moment of inertia),
  • [math]\displaystyle{ L }[/math], unsupported length of column,
  • [math]\displaystyle{ K }[/math], column effective length factor

This formula was derived in 1744 by the Switzerland mathematician Leonhard Euler.[2] The column will remain straight for loads less than the critical load. The critical load is the greatest load that will not cause lateral deflection (buckling). For loads greater than the critical load, the column will deflect laterally. The critical load puts the column in a state of unstable equilibrium. A load beyond the critical load causes the column to fail by buckling. As the load is increased beyond the critical load the lateral deflections increase, until it may fail in other modes such as yielding of the material. Loading of columns beyond the critical load are not addressed in this article.

Around 1900, J. B. Johnson showed that at low slenderness ratios an alternative formula should be used.

Assumptions of the model

Fig. 2: Column effective length factors for Euler's critical load. In practical design, it is recommended to increase the factors as shown above.

The following assumptions are made while deriving Euler's formula:[3]

  1. The material of the column is homogeneous and isotropic.
  2. The compressive load on the column is axial only.
  3. The column is free from initial stress.
  4. The weight of the column is neglected.
  5. The column is initially straight (no eccentricity of the axial load).
  6. Pin joints are friction-less (no moment constraint) and fixed ends are rigid (no rotation deflection).
  7. The cross-section of the column is uniform throughout its length.
  8. The direct stress is very small as compared to the bending stress (the material is compressed only within the elastic range of strains).
  9. The length of the column is very large as compared to the cross-sectional dimensions of the column.
  10. The column fails only by buckling. This is true if the compressive stress in the column does not exceed the yield strength [math]\displaystyle{ \sigma_y }[/math] (see figure 1): [math]\displaystyle{ \sigma = \frac{P_{cr}}{A} = \frac{\pi^2E}{(L_e/r)^2} \lt \sigma_y }[/math] where:
    • [math]\displaystyle{ {L_e}/{r} }[/math] is the slenderness ratio,
    • [math]\displaystyle{ L_e = KL }[/math] is the effective length,
    • [math]\displaystyle{ r = \sqrt{I / A} }[/math]is the radius of gyration,
    • [math]\displaystyle{ I }[/math] is the second moment of area (area moment of inertia),
    • [math]\displaystyle{ A }[/math] is the area cross section.

For slender columns, the critical buckling stress is usually lower than the yield stress. In contrast, a stocky column can have a critical buckling stress higher than the yield, i.e. it yields prior to buckling.

Mathematical derivation

Pin ended column

The following model applies to columns simply supported at each end ([math]\displaystyle{ K = 1 }[/math]).

Firstly, we will put attention to the fact there are no reactions in the hinged ends, so we also have no shear force in any cross-section of the column. The reason for no reactions can be obtained from symmetry (so the reactions should be in the same direction) and from moment equilibrium (so the reactions should be in opposite directions).

Using the free body diagram in the right side of figure 3, and making a summation of moments about point x: [math]\displaystyle{ \Sigma M = 0 \Rightarrow M(x) + Pw = 0 }[/math] where w is the lateral deflection.

According to Euler–Bernoulli beam theory, the deflection of a beam is related with its bending moment by: [math]\displaystyle{ M = -EI\frac{d^2 w}{dx^2}. }[/math]

Fig. 3: Pin ended column under the effect of Buckling load

so: [math]\displaystyle{ EI\frac{d^2w}{dx^2} + Pw = 0 }[/math]

Let [math]\displaystyle{ \lambda^2 = \frac{P}{EI} }[/math], so: [math]\displaystyle{ \frac{d^2w}{dx^2} + \lambda^2 w = 0 }[/math]

We get a classical homogeneous second-order ordinary differential equation.

The general solutions of this equation is: [math]\displaystyle{ w(x) = A \cos(\lambda x) + B \sin(\lambda x) }[/math], where [math]\displaystyle{ A }[/math] and [math]\displaystyle{ B }[/math] are constants to be determined by boundary conditions, which are:

  • Left end pinned: [math]\displaystyle{ w(0) = 0 \rightarrow A = 0 }[/math]
  • Right end pinned: [math]\displaystyle{ w(\ell) = 0 \rightarrow B \sin(\lambda \ell) = 0 }[/math]
Fig. 4: First three modes of buckling loads

If [math]\displaystyle{ B = 0 }[/math], no bending moment exists and we get the trivial solution of [math]\displaystyle{ w(x) = 0 }[/math].

However, from the other solution [math]\displaystyle{ \sin(\lambda \ell) = 0 }[/math] we get [math]\displaystyle{ \lambda_n \ell = n\pi }[/math], for [math]\displaystyle{ n = 0, 1, 2, \ldots }[/math]

Together with [math]\displaystyle{ \lambda^2 = \frac{P}{EI} }[/math] as defined before, the various critical loads are: [math]\displaystyle{ P_{n} = \frac{n^2 \pi^2 EI}{\ell^2} \; , \quad \text{ for } n = 0, 1, 2, \ldots }[/math] and depending upon the value of [math]\displaystyle{ n }[/math], different buckling modes are produced[4] as shown in figure 4. The load and mode for n=0 is the nonbuckled mode.

Theoretically, any buckling mode is possible, but in the case of a slowly applied load only the first modal shape is likely to be produced.

The critical load of Euler for a pin ended column is therefore: [math]\displaystyle{ P_{cr} = \frac{\pi^2 EI}{\ell^2} }[/math] and the obtained shape of the buckled column in the first mode is: [math]\displaystyle{ w(x) = B \sin \left({\pi \over \ell} x\right) . }[/math]

General approach

Fig. 5: forces and moments acting on a column.

The differential equation of the axis of a beam[5] is: [math]\displaystyle{ \frac{d^4 w}{dx^4} + \frac{P}{EI}\frac{d^2 w}{dx^2} = \frac{q}{EI} }[/math]

For a column with axial load only, the lateral load [math]\displaystyle{ q(x) }[/math] vanishes and substituting [math]\displaystyle{ \lambda^2 = \frac{P}{EI} }[/math], we get: [math]\displaystyle{ \frac{d^4 w}{dx^4} + \lambda^2\frac{d^2 w}{dx^2} = 0 }[/math]

This is a homogeneous fourth-order differential equation and its general solution is [math]\displaystyle{ w(x) = A\sin(\lambda x) + B\cos(\lambda x) + Cx + D }[/math]

The four constants [math]\displaystyle{ A, B, C, D }[/math] are determined by the boundary conditions (end constraints) on [math]\displaystyle{ w(x) }[/math], at each end. There are three cases:

  1. Pinned end:
    [math]\displaystyle{ w = 0 }[/math] and [math]\displaystyle{ M = 0 \rightarrow {d^2w \over dx^2} = 0 }[/math]
  2. Fixed end:
    [math]\displaystyle{ w = 0 }[/math] and [math]\displaystyle{ {dw \over dx} = 0 }[/math]
  3. Free end:
    [math]\displaystyle{ M = 0 \rightarrow {d^2w \over dx^2} = 0 }[/math] and [math]\displaystyle{ V = 0 \rightarrow {d^3w \over dx^3} + \lambda^2{dw \over dx} = 0 }[/math]

For each combination of these boundary conditions, an eigenvalue problem is obtained. Solving those, we get the values of Euler's critical load for each one of the cases presented in Figure 2.

See also

References

  1. "Column Buckling | MechaniCalc". https://mechanicalc.com/reference/column-buckling. 
  2. Euler, Leonard (1744) (in Latin). Methodus inveniendi lineas curvas maximi minimive proprietate gaudentes sive solutio problematis isoperimetrici latissimo sensu accepti. Geneva, Switzerland: Marc Michel Bousquet et Cie.. pp. 267-268. https://archive.org/details/methodusinvenie00eule/page/266/mode/2up.  From pp. 267-268: "37. Quae ante de specie prima sunt annotata inservire possunt viribus columnarum dijudicandis. […] contra vero si pondus P fuerit majus, columna incurvationi resistere non poterit." (37. Those [things] which have been noted before about the first type can serve in judging the strength of columns. So let the column AB be placed vertically on the base A, bearing the load P. For if the column is already arranged in such a way that it cannot slide [away] from the load P, [then] if [the load] will have been too great, there will be nothing else to fear except the bending of the column; in this case, therefore, the column may be regarded as endowed with elasticity. Therefore let the absolute elasticity of the column = Ekk, and its height AB = 2f = a; and [in] §25 above [p. 261] we have seen that the required bending force on this column or the minimum [load that is required to bend this column] = ππ Ekk/4ff = ππ Ekk/aa. Thus, unless the load P being carried is greater than Ekk/aa, absolutely no bending will have to be feared; but on the other hand, if the load P will have been greater, the column will not be able to resist bending.)
  3. "Twelve Viva Questions on Columns and Struts" (in en). 2015-03-28. http://engineering.myindialist.com/2015/twelve-viva-questions-on-columns-and-struts/. 
  4. "Buckling of Columns". http://web.aeromech.usyd.edu.au/AMME2301/Documents/Chapter09.pdf. 
  5. Timoshenko, S. P.; Gere, J. M. (1961). Theory of Elastic Stability (2 ed.). McGraw-Hill.