# P-Delta Analysis and Geometric Non-linearity

by Dr Seán Carroll
|
Published: 23rd November 2021
|
Tutorial

In this tutorial, we’ll consider P-Delta analysis (); a form of non-linear behaviour that can lead to large magnitude sway deflections in columns. behaviour occurs in members subject to compression and it presents a particular challenge for columns within tall flexible structures. If you want to download the complete Jupyter Notebook for this tutorial, just follow the link below.

## 1.0 What is the P-Delta effect?

Put simply, describes the phenomenon whereby an additional or secondary moment is generated in a column due to the combination of axial load and lateral sway . This leads to non-linear structural behaviour and can result in lateral deflections far in excess of those arising from lateral loading alone. These so-called second-order deflections will induce additional stresses within the structure that may be significant and require special consideration in design.

The effect is referred to as a geometric non-linearity because it arises as a result of the deformed geometry of the structure. This is in contrast to a material non-linearity such as plastic hinge formation that arises due to the properties of the material.

## 2.0 Secondary moments caused by sway deflections

To flesh out the concept, consider a column segment of length subject to an axial force and undergoing a relative sway between its ends. The sway may have been caused by wind loading or inertia forces due to lateral ground motion. Regardless of how the initial sway deflection came to be, the key point is that the compression forces in the column are no longer co-linear.

We can think of the axial forces, shears and moments shown above as the primary actions on the structure consistent with the first-order sway deflection, . Again, notice that the axial forces, are no longer co-linear. As a result, an extra second-order moment is developed at the base of the column,

(1)

If is a non-negligible value, we must also consider the additional secondary sway deflection caused my . Since the second-order sway deflection further increases the overall sway, leading to yet more secondary sway, we have a feedback loop that could ultimately lead to collapse.

For a tall flexible structure that may undergo significant sway from one storey to the next (storey drift) under lateral base excitation for example, this has catastrophic potential. Only when we are sure that the first-order sway deflection is negligibly small can second-order deflections be ignored.

The effect is a good example of non-linear structural behaviour arising from geometric non-linearity. The typical method of analysis requires iteration to determine the final value that the sway deflection converges on.

This value of sway deflection may be mathematically stable but represent an unsustainable physical configuration for the column. In other words, the second-order deflection may be so high that the column collapses. Ultimately the deciding factor as to whether collapse will occur is the initial sway deflection, , as this sets in motion the feedback loop that delivers the additional second-order deflection.

Since this is a non-linear structural behaviour wherein the inputs to the structural system are not linearly proportional to the outputs (the structure’s behaviour), superposition should not be used in the analysis.

## 3.0 P-Delta Analysis Example: Fixed-Free Column

Next, we can demonstrate the iterative nature of analysis with a simple example. Consider a simple column with one end fixed and the other free, subject to a vertical load . We will assume that this column experiences a lateral load at that induces an initial first-order deflection .

Let’s assume the column has the following parameters,

We can determine the Euler buckling load for this column to work out a sensible upper limit for the axial load. We can use the following formula for a fixed-free column,

(2)

Note that we’re assuming buckling about the minor axis. This is why we selected from the section tables.

Let’s assume the column is subject to a compression force of . Let’s also assume that the column is subjected to a lateral force magnitude .

We can calculate the first-order deflection (hereafter referred to as ) due to the lateral load using he standard formula for cantilever deflection,

(3)

where is the moment generated at the base of the column by the lateral load .

We know that this first-order sway deflection induces an additional second-order moment that has a value of at the base of the column. If we examine the internal moment at some distance from the base of the column, we can see that this moment actually varies between at the tip of the column and at the base of the column, with the moment diagram having the same shape as the first-order deflected shape, , Fig. 3.

At this point we can make a simplification and assume that the second-order moment varies linearly between and , Fig.4.

If we do this, we can quickly determine the additional deflection induced by the second-order moment using the same equation that was used to calculate ,

(4)

where is the first value of second-order deflection calculated. We refer to this value of as the first iteration value, .

Now the iterative nature of the analysis becomes apparent as we need to calculate an additional second-order deflection, induced by this first iteration value, ; in other words, our axial load now generates a further second-order moment, , leading to a second iteration which yields a further second-order deflection, .

The entire process must be repeated yet again with each subsequent value of being added to the preceding value. The iterations continue until the new values of second-order deflection being added become sufficiently small. In this way the second-order deflection converges to its final value. We must then judge the acceptability of the total lateral deflection,

(5)

### 3.1 Writing a function to calculate the secondary deflection

Rather than continuing to perform iterations manually, let’s write a short function to perform the iterations and identify the final second-order deflection.

The function defined above continues to perform iterations to identify second-order deflection values until the value of second-order deflection calculated falls below a threshold percentage of the total column deflection. The function then returns the array of all second-order deflections calculated. Next, we need to call the function and print out the array of deflections calculated.

The total deflection is then easily obtained as the sum of the first-order and all second-order deflections.

Next we can plot the deflection values on each iteration to observer the convergence.

Now that the additional second-order sway deflection has been determined, subsequent analysis can be carried out to determine the acceptability of the stresses induced. The take away message is the need to be aware of second-order effects when considering the behaviour of slender structures that must resist large axial loads.

## 4.0 Exploring the P-Delta parameter space

Now that we have a function to calculate the second-order deflection, we can perform a parameter sweep to get a better understanding of the behaviour of the column within the parameter space. After all, we’ve only considered a single combination of axial and lateral load. It would be nice to plot the lateral deflection as a function of these input parameters. We start by defining a range of axial loads as a function of the column critical load and then a range of lateral loads as a proportion of the axial load.

We can now iterate through these values; looping through the axial load with an outer for loop and then the lateral loads with a nested for loop. For every combination of axial and lateral load we will calculate the maximum lateral deflection and the ratio of second to first-order deflection to get a sense of the degree of non-linearity for each parameter combination.

We can see from the first plot above that for any particular value of axial load, the total deflection increases linearly with increasing lateral load. This is not at all surprising, however we also note that the rate of increase (slope of each line) increases as the axial load increases.

If we plot the ratio of second to first-order deflection for each load combination (second plot) we can see that for each value of axial load the ratio of second to first-order defection remains constant but as the axial load increases linearly, the increase in the deflection ratios is not linear. We can see this more clearly if we simply plot a single value of deflection ratio against each value of axial load.

We can clearly see that the non-linear deflection exceeds the linear deflection as the axial load increases; the column deflection is quickly dominated by the non-linear component of deflection as the axial load increases. Even though in this example we’ve based our analysis on a column that is particularly prone to non-linear deflection, this highlights the importance of not neglecting second-order sway deflections for slender structures that must resists large axial loads.

## 5.0 P-Delta analysis for more complex structures

The example above is a pretty contrived but simple and common example used to explain and illustrate effects. However, the question now arises; how do we apply what we’ve learned to more complex structures?

At this point, the road forks and there are two approaches. In each case we assume that you have a model of your structure that you can perform a stiffness method analysis on (check out these courses for building stiffness method analysis tools).

Iterative approach

The first approach is to implement basically the same solution strategy described above:

1. Apply any lateral load, to calculate first-order sway deflection, for the structure.
2. Manually displace the structural model by applying the nodal displacements to each node in the model according to nodal deflections.
3. Apply any vertical loads to the deformed structure to calculate second-order deflections, and note the second-order nodal deflections obtained.
4. Displace the structure (from its undisplaced position) according to nodal deflections.
5. Apply vertical loads again to the deformed structure to calculate second-order deflections, .
6. Repeat steps 4 and 5 until the second-order deflections observed have become negligibly small.
7. Calculate the total deflection, .

Geometric stiffness matrix approach

A more sophisticated and efficient approach would be to modify the element stiffness matrix to take into account geometric non-linearity. This will result in an element stiffness matrix that exhibits coupling between the axial force and bending moments in the element. Derivation of this stiffness matrix is beyond the scope of this tutorial, but once obtained, this approach will allow the structure to be solved without iteration. Most commercial frame solvers will have the option of taking into account geometric non-linearity and will calculate the appropriate element stiffness matrix.

Well that’s all for this tutorial. I hope you found it helpful and now have a good understanding of effects as a form of geometric non-linearity and how to calculate them in practice. Remember that you can download the complete Jupyter Notebook for this tutorial for free by following this link.

## Dr Seán Carroll

BEng (Hons), MSc, PhD, CEng MIEI, FHEA

Hi, I’m Seán, the founder of DegreeTutors.com. I hope you found this tutorial helpful. After spending 10 years as a university lecturer in structural engineering, I started DegreeTutors.com to help more people understand engineering and get as much enjoyment from studying it as I do. Feel free to get in touch or follow DegreeTutors on any of the social accounts.