Linearized Buckling Analysis

A linearized buckling analysis can be used for estimating the critical load at which a structure becomes unstable. This is a predefined study type that consists of two study steps: An initial step in which a unit load is applied to the structure, and a second step in which an eigenvalue problem is solved for the critical buckling load.

The idea behind this type of analysis can be described in the following way:

Consider the equation system to be solved for a stationary load f,

Here the total stiffness matrix, K, has been split into a linear part, KL, and a nonlinear contribution, KNL.

In a first order approximation, KNL is proportional to the stress in the structure and thus to the external load. So if the linear problem is solved first for an arbitrary initial load level f0,

(2-14)

then the nonlinear problem can be approximated as

where λ is called the load multiplier.

An instability is reached when this system of equations becomes singular so that the displacements tend to infinity. The value of the load at which this instability occurs can be determined by, in a second study step, solving an eigenvalue problem for the load multiplier λ.

COMSOL reports a critical load factor, which is the value of λ at which the structure becomes unstable. The corresponding deformation is the shape of the structure in its buckled state.

The level of the initial load used is immaterial since a linear problem is solved. If the initial load actually was larger than the buckling load, then the critical value of λ is smaller than 1. It is also possible that the computed value of λ is negative. This signifies that a reversed load will give the critical case.

Geometric Nonlinearity

Sometimes, the preload case requires a geometrically nonlinear analysis in order to produce the correct state. This means that it is no longer solved using the linear set of equations given by Equation 2-14, but rather

(2-15)

The assumption for the buckling analysis is still that KNL is proportional to the external load, even though this may be disputable for a strongly nonlinear case. KNL is based on the stresses, which must be computed in the same way for both cases, that is under the same assumption about geometric nonlinearity. The effect is that the stiffness matrix at the linearization point includes the nonlinear part from Equation 2-15, and the eigenvalue problem is reformulated as


	If the computation of the predeformation requires a geometrically nonlinear analysis, then geometric nonlinearity must be used also in the Linear Buckling study step. In this case, it must be assumed the critical load factor λ is significantly larger than 1.


	Be aware that for some structures, the true buckling load can be significantly smaller that what is computed using a linearized analysis. This phenomenon is sometimes called imperfection sensitivity. Small deviations from the theoretical geometrical shape can then have a large impact on the actual buckling load. This is especially important for curved shells.


	For a structure that exhibits axial symmetry in the geometry, constraints, and loads, the critical buckling mode shape can still be nonaxisymmetric. A full 3D model should always be used when computing buckling loads.

•

Studies and Solvers and Linear Buckling in the COMSOL Multiphysics Reference Manual

•

Linear Buckling in the theory section of the Structural Mechanics Module User’s Guide

•

Bracket — Linear Buckling Analysis: Application Library path

•

Linear Buckling Analysis of a Truss Tower: Application Library path