5 - Buckling of Columns
5 - Buckling of Columns
𝑀 2·𝛽 ·𝜃 … 𝐸𝑞𝑢. 1
F.B.D
Since the angle 𝜃 is a small quantity, the lateral displacement of point B is 𝜃 L/2.
Therefore, the equation of equilibrium is obtained by summing moments about
point B for bar BC.
Second, consider bar BC as a free body and note that it is subjected to the action
of the axial forces P and the moment M in the spring. The moment 𝑀 is equal
to the rotational stiffness 𝛽 times the angle of rotation 2𝜃 of the spring; thus,
𝜃𝐿
𝑀 𝑃 0 … Equ. 2
2
𝑃𝐿 4·𝛽
2·𝛽 𝜃 0 … Equ. 3 𝑃 … Equ. 4
2 𝐿
At higher values of the load, the column is unstable and may collapse by buckling,
that is, by excessive bending. For the ideal case, the column will be in equilibrium in
the straight position even when the axial force P is greater than the critical load.
However, since the equilibrium is unstable, the smallest imaginable disturbance will
cause the column to deflect sideways. Once that happens, the deflections will
immediately increase, and the column will fail by buckling.
Differential Equation for Column Buckling
𝐸𝐼𝑣 𝑀 …(Equ. 1)
Where:
𝑀 𝑃𝑣 0 or 𝑀 𝑃𝑣 …(Equ. 2)
𝐸𝐼𝑣 𝑃𝑣 0 …(Equ. 3)
Solution of the Differential Equation
𝑘 or 𝑘 …(Equ. 4)
𝑣 𝑘 𝑣 0 …(Equ. 5)
𝐶 0, so 𝑣 𝐶 𝑠𝑖𝑛 · 𝑘𝑥 … 𝐸𝑞𝑢. 7
𝐶 𝑠𝑖𝑛 · 𝑘𝐿 0 … Equ. 8
Case 2:
The second possibility for satisfying Equation 8 is given by the equation
known as the buckling equation:
𝑠𝑖𝑛 · 𝑘𝐿 0 … 𝐸𝑞𝑢. 9
𝑛𝜋𝑥
𝑣 𝐶 𝑠𝑖𝑛 · 𝑘𝑥 𝐶 𝑠𝑖𝑛 𝑛 1, 2, 3, … … 𝐸𝑞𝑢. 12
𝐿
2.1 Critical Loads
The lowest critical load for a column with pinned ends is obtained when 𝑛 1
𝜋 𝐸𝐼
𝑃 … 𝐸𝑞𝑢. 13
𝐿
𝜋𝑥
𝑣 𝐶 𝑠𝑖𝑛 … 𝐸𝑞𝑢. 14
𝐿
• The constant 𝐶 represents the deflection at the midpoint of the column and may have
any small value , either +ve or –ve.
• Therefore, the part of the load‐deflection diagram corresponding to P is a horizontal
straight line.
• Thus, the deflection at the critical load is undefined, although it must remain small for our
equations to be valid. Above the bifurcation point B, the equilibrium is unstable, and
below point B it is stable.
• Buckling of a pinned‐end column in the first mode is called the fundamental case of
column buckling.
• The type of buckling described in this section is called Euler buckling, and the
critical load for an ideal elastic column is often called the Euler load.
• The famous mathematician Leonhard Euler (1707–1783), generally recognized as
the greatest mathematician of all time, was the first person to investigate the
buckling of a slender column and determine its critical load.
• Taking higher values of the index 𝑛 in Equations 11 and 12 leads to an infinite
number of critical loads and corresponding mode shapes. The mode shape for 𝑛
2 has two half‐waves. The corresponding critical load is four times Larger than the
critical load for the fundamental case.
• The magnitudes of the critical loads are proportional to the square of 𝑛, and the
number of half‐waves in the buckled shape is equal to 𝑛.
• Buckled shapes for the higher modes are often of no practical interest because the
column buckles when the axial load P reaches its lowest critical value.
• The only way to obtain modes of buckling higher than the first is to provide lateral
support of the column at intermediate points, such as at the midpoint of the
column
• The flexural rigidity can be increased by using a “stiffer” material (that is, a
material with larger modulus of elasticity E)
• Reducing the wall thickness of a tubular member and increasing its lateral
dimensions (while keeping the cross‐sectional area constant) also increases the
critical load because the moment of inertia is increased.
2.2 Critical Stress
After finding the critical load for a column, calculate the corresponding critical
stress by dividing the load by the cross‐sectional area. For the fundamental case
of buckling, the critical stress
𝑃 𝜋 𝐸𝐼
𝜎 … 𝐸𝑞𝑢. 15
𝐴 𝐴𝐿
In which 𝐼 is the moment of inertia for the principal axis about which buckling occurs. This
equation can be rewritten in a more useful form by introducing the notion
𝐼
𝑟 . . . 𝐸𝑞𝑢. 16
𝐴
in which 𝑟 is the radius of gyration of the cross section in the plane of bending. The
equation for critical stress becomes
𝜋 𝐸
𝜎
𝐿
𝑟
𝐿
𝑆𝑙𝑒𝑛𝑑𝑒𝑟𝑛𝑒𝑠𝑠 𝑟𝑎𝑡𝑖𝑜
𝑟
The slenderness ration only depends on the dimensions of the column.
A column that is long and slender will have a high slenderness ratio and therefore a low critical
stress.
A column that is short and stubby will have a low slenderness ratio and will buckle at a high stress
Typical values for slenderness ratio for actual columns are between 30 and 150.
The critical stress is the average compressive stress on the cross section at the instant the load
reaches its critical value
A graph of this stress as a function of slenderness ratio is obtained and a curve known as Euler’s
curve.
The curve is valid only when the critical stress is less than the proportional limit, because the
equations were derived using Hooke’s law.
2.3 Effects of large Deflections,
Imperfections, and Inelastic Behavior
• The equations for critical loads were derived from ideal
columns.
• As a consequence, the magnitudes of the small
deflections at buckling were undefined.
• When 𝑃 𝑃 , the column may have any small
deflection, which is a condition represented by the
horizontal line “A” in the load – deflection diagram
• The theory for ideal columns is limited to small
deflections as second derivative 𝑣 for the curvature
was used.
• A more exact analysis, based upon the exact expression
for curvature shows that there is no indefiniteness in
the magnitudes of the deflections at buckling.
• Instead, for an ideal, linearly elastic column, the load‐
deflection diagram goes upwards in accord with curve
“B”.
• Thus, after a linearly elastic column begins to buckle, an
increasing load is required to cause an increase in the
deflections.
Suppose a column is not constructed perfectly, the column might have an
imperfection in the form of a small initial curvature, so that the unloaded
column is not constructed perfectly straight.
The larger the imperfections, the further curve “C” moves to the right,
away for the vertical line.
By comparing lines “A, B and C”, note that for practical purposes the critical
load represents the maximum load‐carrying capacity of an elastic column,
because large deflections are not acceptable in most applications.
When the stresses exceed the proportional limit and the material
no longer follows Hooke’s law, the load‐deflection diagram is
unchanged up to the level of load at which the proportional limit
is reached
Then the curve for inelastic behavior “D” departs from the elastic
curve, continues upward, reaches a maximum, and turns
downward.
A column shaped as shown in the figure on the right will have a larger
critical load than a prismatic column made from the same volume of
material.
• Consider a prismatic column with pinned ends that is
free to buckle in any lateral direction.
• Assume the column has a solid cross section, which
cross section gives the largest critical load.
Buckling of a column with pinned ends is usually considered as the most basic case of buckling.
However, in practice, engineers encounter many other end conditions such as fixed ends, free ends, and
elastic supports.
The critical loads for columns with various kinds of support conditions can be determined from the
differential equation of the deflection curve by following the same procedure used to analyse the pinned‐
end column.
Procedure:
1. Obtain an expression for the bending moment in the column by
assuming the column to be in the buckled state.
2. Set up the differential equation of the deflection curve, sing the bending‐
moment equation
3. Solve the equation and obtain its general solution, which contains two
constants of integration plus any other unknown quantities.
4. Apply boundary conditions pertaining to the deflection 𝑣 and the slope
𝑣 and obtain a set of simultaneous equations.
5. Solve the equations to obtain the critical load and the deflected shape of
the buckled column
3.1. Column Fixed at the base and free at the Top
Consider an ideal column that is fixed at the base, free at the top, and subjected
to an axial load “P”. The bending moment at distance 𝑥 from the base is:
𝑀 𝑃 𝛿 𝑣 … 𝑎
Where 𝛿 is the deflection at the free end of the column. The differential equation
of the deflection curve the becomes
𝐸𝐼𝑣 𝑀 𝑃 𝛿 𝑣 … 𝑏
Where 𝐼 is the moment of inertia for buckling in the x‐y plane. Using the notation
𝑘 , equation (b) rearranges into the form
𝑣 𝑘 𝑣 𝑘 𝛿 … 𝑐
𝑣 𝐶 sin 𝑘𝑥 𝐶 𝑐𝑜𝑠 𝑘𝑥 … 𝑑
𝑣 𝛿 … 𝑒
When 𝑣 is substituted into the left –hand side of the differential equation, it
produces the right‐hand side, the term 𝑘 𝛿.
𝑣 𝐶 sin 𝑘𝑥 𝐶 𝑐𝑜𝑠 𝑘𝑥 𝛿 … 𝑓
𝑣 0 0 𝑣 0 0
𝐶 𝛿 … 𝑔
To apply the second condition, first differentiate equ. (f) to obtain the slope:
𝑣 𝐶 𝑘 𝑐𝑜𝑠 𝑘𝑥 𝐶 𝑘 𝑠𝑖𝑛 𝑘𝑥
Now substitute the expressions for 𝐶 and 𝐶 into the general solution of equ. (f),
we obtain the equation of the deflection curve for the buckled column:
𝑣 𝛿 1 cos 𝑘𝑥 … ℎ
The third boundary condition applies to the upper end of the column, where the
deflection 𝑣 𝛿.
𝑣 𝐿 𝛿
𝛿 cos 𝑘𝐿 0 … 𝑖
If 𝛿 0, there is no deflection and you have the trivial solution – the column
remains straight and buckling does not occur.
𝑛𝜋
𝑘𝐿 𝑛 1, 3, 5, … 𝑘
2
The expression 𝑘 is used to obtain the critical load formula:
𝑛 𝜋 𝐸𝐼
𝑃 𝑛 1, 3, 5, … … 𝑙
4𝐿
𝑛𝜋𝑥
𝑣 𝛿 1 𝑐𝑜𝑠 𝑛 1, 3, 5, … 𝑚
2𝐿
𝜋 𝐸𝐼
𝑃 … 𝑛
4𝐿
𝜋𝑥
𝑣 𝛿 1 𝑐𝑜𝑠 … 𝑜
2𝐿
Taking higher values of the index “n” leads to an infinite number of critical loads
from equ. (l). The corresponding buckled mode shapes have additional waves in
them
3.2. Effective Lengths of Columns
• The critical loads for columns with various support conditions
can be related to the critical load of a pinned‐end column
through the concept of an effective length.
• Consider a deflected shape of a column fixed at the base and
free at the top.
• This column buckles in a curve that is 1/4 of a complete sine
wave.
• If the deflection curve is extended, it becomes 1/2 of a
complete sine wave, which is the deflection curve for a pinned‐
end column.
• The effective length 𝐿 for any column is the length of the
equivalent pinned‐end column.
• Thus for a fixed‐free column, the effective length is:
𝐿 2𝐿 … 1
𝜋 𝐸𝐼
𝑃 … 2
𝐿
• If the effective length of the column is known (no matter how complex the end
conditions may be), you can substitute into the preceding equation (2) and
determine the critical load.
𝐿 𝐾𝐿 … 3
Where 𝐿 is the actual length of the column. Thus the critical load is:
𝜋 𝐸𝐼
𝑃 … 4
𝐾𝐿
The factor K =2 for a column fixed at the base and free at the top
K = 1 for a pinned‐end column
3.3. Columns with Both Ends Fixed Against Rotation
𝐿
𝐿 … 𝑖
2
4𝜋 𝐸𝐼
𝑃 … 𝑖𝑖
𝐿
3.4. Column Fixed at the Base and Pinned at the Top
When a column buckles, a reactive moment 𝑀 develops at the base, as there can be no rotation at that point.
From the equilibrium of the entire column, there must be horizontal reactions R at each end such that
𝑀 𝑅𝐿
The bending moment inn the buckled column, at distance 𝑥 from the base is,
𝑀 𝑀 𝑃𝑣 𝑅𝑥 𝑃𝑣 𝑅 𝐿 𝑥 … 𝑖
𝐸𝐼𝑣 𝑀 𝑃𝑣 𝑅 𝐿 𝑥 … 𝑖𝑖
𝑅
𝑣 𝑘 𝑣 𝐿 𝑥 … 𝑖𝑖𝑖
𝐸𝐼
The general solution of this equation is
𝑅
𝑣 𝐶 sin 𝑘𝑥 𝐶 cos 𝑘𝑥 𝐿 𝑥 … 𝑖𝑣
𝑃
The solution contains three unknown quantities (𝐶 , 𝐶 , and R). The following three boundary conditions are
required:
𝑣 0 0 𝑣 0 0 𝑣 𝐿 0
𝑅𝐿 𝑅
𝐶 0 𝐶𝑘 0 𝐶 𝑡𝑎𝑛 𝑘𝐿 𝐶 0 … 𝑣
𝑃 𝑃
To obtain the solution for buckling, equ. (v) is solved in a more general manner, by eliminating R from thr forst two
equations, which yields
𝐶 𝑘𝐿 𝐶 0 𝑜𝑟 𝐶 𝐶 𝑘𝐿 … 𝑣𝑖
Substitute 𝐶 expression into equ. (v) (last expression), we obtain the buckling equation
𝑘𝐿 𝑡𝑎𝑛𝑘𝐿 . . 𝑣𝑖𝑖
The buckling equation is a transcendental equation, it cannot be solved explicitly. Nevertheless, the values of 𝑘𝐿 that
satisfy the equation can be determined numerically by using computer program for find the roots id equations. The
smallest nonzero value of 𝑘𝐿 that satisfies equ. (vii) is
𝑘𝐿 4.4934
The corresponding critical load is
20.19𝐸𝐼 2.046𝜋 𝐸𝐼
𝑃 𝑣𝑖𝑖
𝐿 𝐿
Which is higher than the critical load for a column with pinned ends and lower than
the critical load for a column with fixed ends.
The effective length of the column is (distance from the pinned end of the column
to the point of inflection in the buckled shape.
The equation of the buckled mode shape is obtained substituting equ. (v) in the
general solution equ. (iv):
𝑣 𝐶 sin 𝑘𝑥 𝑘𝐿 cos 𝑘𝑥 𝑘 𝐿 𝑥 … 𝑖𝑥
Inn which
4.4934
𝑘
𝐿
Critical loads, effective lengths, and effective‐length factors for ideal columns
Limitation: The Euler buckling theory used in this section is only valid if the
column is perfectly straight before the load is applied, the column and its
supports have no imperfections, and the column is made of a linearly elastic
material that follows Hooke’s las.
4. COLUMNS WITH ECCENTRIC AXIAL LOADS
When a column is compressed by loads 𝑃 that are applied with a small
eccentricity 𝑒 measured from the axis of the column.
Each eccentric axial load is equivalent to a centric load 𝑃 and a couple of
moment 𝑀 𝑃𝑒
This moment exists from the instant the load is first applied; therefore, the
column begins to deflect at the onset of loading.
The deflection then becomes steadily larger as the load increases.
𝑀 𝑀 𝑃 𝑣 𝑃𝑒 𝑃𝑣 𝑎
𝐸𝐼𝑣 𝑀 𝑃𝑒 𝑃𝑣 … 𝑏
or
𝑣 𝑘 𝑣 𝑘 𝑒 … 𝑐
𝑣 𝐶 sin 𝑘𝑥 𝐶 𝑐𝑜𝑠 𝑘𝑥 𝑒 … 𝑑
In which 𝐶 and 𝐶 are constants of integration in the homogeneous solution and 𝑒 is the
particular solution.
The boundary conditions for determining the constants 𝐶 and 𝐶 are obtained from
the deflections at the ends of the column:
𝑣 0 0 𝑣 𝐿 0
These conditions yield
𝑒 1 cos 𝑘𝐿 𝑘𝐿
𝐶 𝑒 𝐶 𝑒𝑡𝑎𝑛
sin 𝑘𝐿 2
𝑘𝐿
𝑣 𝑒 tan sin 𝑘𝑥 cos 𝑘𝑥 1 … 𝑒
2
4.1. Maximum Deflection
The maximum deflection 𝛿 produced by the eccentric loads occurs at the mid‐point of
the column and is obtained by setting 𝑥 in equ. (e).
𝑘𝐿
𝛿 𝑒 sec 1 … 𝑓
2
The above equation can be rewritten in a different form by replacing 𝑘 with its equivalent
in terms of critical load:
𝑃 𝑃𝜋 𝜋 𝑃
𝑘 … 𝑔
𝐸𝐼 𝑃 𝐿 𝐿 𝑃
𝑀 𝑃 𝑒 𝛿 .. 𝑗
𝑘𝐿 𝜋 𝑃
𝑀 𝑃 · 𝑒 · 𝑠𝑒𝑐 𝑃 · 𝑒 · 𝑠𝑒𝑐 .. 𝑘
2 2 𝑃
The manner in which 𝑀 varies as a function of the axial load 𝑃 is shown in the figure below
Acting at this cross section are the compressive force P and the bending moment 𝑀 .
‐The stresses due to the force P are equal to , where A is the cross‐sectional area of the
column.
𝑃 𝑀 ·𝑐
𝜎 … 𝑖
𝐴 𝐼
𝜋 𝑃
𝑀 𝑃 · 𝑒 · 𝑠𝑒𝑐
2 𝑃
Since 𝑃 𝜋 for a pinned‐end column, and since 𝐼 𝐴𝑟 , where 𝑟 is the radius of gyration
in the plane of bending, the preceding equation becomes
𝐿 𝑃
𝑀 𝑃 · 𝑒 · 𝑠𝑒𝑐 … 𝑖𝑖
2𝑟 𝐸𝐴
Substitute equ. (i) into (ii) we obtain the maximum compressive stress formula, also known as
the secant formula;
𝑃 𝑃·𝑒·𝑐 𝐿 𝑃
𝜎 𝑠𝑒𝑐
𝐴 𝐼 2𝑟 𝐸𝐴
Or
𝑃 𝑒·𝑐 𝐿 𝑃
𝜎 1 𝑠𝑒𝑐 … 𝑖𝑖𝑖
𝐴 𝑟 2𝑟 𝐸𝐴
The eccentricity ratio is a measure of the eccentricity of the load as compared to the
dimensions of the cross section.
Its numerical value depends upon the position of the load, but typical values are in the range
from 0 – 3 with the most case value being less than 1.
𝑒𝑐
𝐸𝑐𝑐𝑒𝑛𝑡𝑟𝑖𝑐𝑖𝑡𝑦 𝑟𝑎𝑡𝑖𝑜 … 𝑖𝑣
𝑟
‐THE END ‐