Linear and Nonlinear Buckling in FEA
Linear and Nonlinear Buckling in FEA
Linear Buckling
The most basic form of buckling analysis in FEA is linear buckling. This is directly related to the
classic Euler type of calculation. A small displacement of a perturbed shape is assumed in each
element that induces a stress dependent stiffening effect. This adds to the linear static stiffness
in the element. Imagine a guitar string tightened the strings total stiffness goes up and
results in a higher pitch. If the string is slackened the total stiffness goes down, and the pitch
corresponds. The stress dependent stiffness is now subtracting from the linear static stiffness
term. This latter effect causes buckling.
In an assembly of elements in an FEA model there will be a subtle interaction between the
original linear stiffness matrix and the stress dependent stiffness matrix. This is analogous to the
linear stiffness matrix and the mass matrix in a normal modes analysis. The same solution
method is used an eigenvalue extraction. For a linear buckling analysis, this will find what
scaling factors applied to the nominal static load will scale the stress stiffening terms to subtract
sufficiently from the linear static terms to give unstable solutions.
Fig. 2: Typical FEA buckling analysis set up, axially loaded cylinder.
Fig. 2 shows a typical analysis setup. An axial load of a nominal 1KN is applied to the top of a
thin-walled cylinder. The constraint systems are shown. A linear buckling analysis is carried out.
The stress stiffening matrix and the linear static stiffness matrix are calculated in the first linear
static step. In the second step the unstable roots are found using the two matrices in an
eigenvalue solution. The user doesnt have to do anything special here because all buckling
solvers are hard-wired to do the two steps.
The result of the analysis is a table of eigenvalues as shown in Fig. 3 and a set of mode
shapes, or eigenvectors. The first mode shape is shown in Fig. 4.
Fig. 6: Nonlinear buckling showing stress transition just below and at buckling.
The information we get is very useful in design, but it is more of an indicator than a hard
number. We also have to be aware that if we use linear buckling on a structure that is more like
the intermediate category, then we are likely to get a non-conservative over estimate of the
buckling load. We may also find the mode shape transitions very quickly into something very
different.
The boundary condition assumptions for buckling are also critical. However, if the structure can
be categorized as slender and we can show a good margin over the critical linear buckling
load, then in many cases that is sufficient for design.
Fig. 7 shows the higher order modes shapes associated with eigenvalues 1 to 10. Very often
the default in an FEA solver is to just have the first eigenvalue and mode shape. In fact, the
study of the higher modes is useful.
How do we overcome this limitation? Some components and loading will have such a large
natural eccentricity that the solution will find instability. However for our stubborn cylinder we
have to introduce an eccentricity. There are several ways of doing this. All methods can benefit
from our understanding of the linear buckling mode. The nonlinear mode may transition through
this, but it is a good starting point. We can either apply a dummy loading to induce the mode
shape, or actually distort the structure very slightly in favor of the mode shape. The first method
is usually easiest, as any sympathetic load will usually work. Pressures are better than point
loads as they avoid local singularities. If possible a sympathetic pressure can be applied in the
same distribution as the normal displaced mode shape from the linear analysis. It can be
captured as a field function and scaled to suit. It is difficult to assess what level of load to apply,
but it should be a lot smaller than the main axial loading.
In the case of the cylinder, I applied a pressure over one quarter vertical strip to give a net
sideways thrust. This dummy static load was set to give a deflection of 3.5e-4 m a very small
wall deflection. The best way to do this is to ramp up the dummy load as the main loading is
applied, and then to ramp it down to zero until 100% main load is achieved. Alternatively, a prestress load case may be possible, which will not get scaled up as the main loading is increased.
The setup will depend on your FEA program. If all else fails you may have to apply a small
dummy load, which will increase with the main load and be there as a small artifact at the end.
The other alternative is to capture the linear buckling mode shape and apply this back to the
structural mesh as an initial distortion. Again this can use a field function from the mode which is
scaled and used as an offset, or the same thing can be done via an Excel export. Your FEA
solver may well have automated ways of handling this. A word of caution here though: if you
make the imperfections too big you will effectively corrugate the structure and stiffen it up. Trial
and error is required.
Fig. 9 shows the results of the nonlinear buckling analysis. The cylinder configuration and the
level of eccentricity assumed result in a very stable structure that resists buckling until a mode
occurs, similar in nature to the linear mode. There is then a transition to a highly localized mode.
The initial buckling occurs at around 63% of the linear estimate for critical load. The maximum
load level achievable before the structure sees a local collapse is about 82% of the linear critical
load.
The nonlinear method used to track the geometric instability was the Arc Length method. This is
always recommended as it enables the softening response of the structure to be tracked. When
a structure buckles it has less resistance to load. A real-world load controlled test would result in
a runaway collapse. Alternatively, a real-world displacement controlled test would allow us to
monitor the reducing load resistance. The Arc Length method permits this, even though we are
applying load.
Important considerations are: What is the effective critical buckling load level and what is the
post-buckling behavior? I like to use key point plots to investigate this. Fig. 10 shows a shell
element identified as being a good indicator of the changing stress state as instability occurs.
Membrane hoop stress is zero and axial stress is steadily increasing up to the start of instability
at point A. This is a classic linear response. From A to B the hoop stress increases as the
structure distorts. From B to C a second mode occurs with a further instability local to this
element, the axial stress decreases as the structure softens. The overall response in Fig. 9
shows the transition to local dimpling. In fact this region would rapidly exceed yield in practice,
but a linear material model has been used.
10
Fig. 11: Key point plot; radial displacement response at key nodes.
In practice, the model should now be investigated for sensitivity to initial induced eccentricity
and preferably comparing several forms of induced eccentricity. Effects of constraint and loading
implications can be compared to the real-world conditions by experimenting with DOF (degrees
of freedom) and using boundary spring stiffnesses. The load steps can be adjusted to give finer
results closer to initial instability. Plastic behavior could also be investigated in the transition to
the second instability.
Conclusion
Buckling is a critical failure condition for many classes of structure. Accurate estimates of critical
load and response modes are difficult unless a structure falls well into the slender category.
Linear solutions may suit such structures if loads and boundary conditions are carefully
assessed. However for the majority of instability prone structures a full nonlinear analysis is
required. This type of analysis is very sensitive to assumptions on eccentricity and boundary
conditions. A methodology is required that will deal with structural softening. The key point
method is recommended to help identify the onset of instability and subsequent transitional
modes.
11