fbpx

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 (P-\Delta); a form of non-linear behaviour that can lead to large magnitude sway deflections in columns. P-\Delta 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, P-\Delta describes the phenomenon whereby an additional or secondary moment is generated in a column due to the combination of axial load (P) and lateral sway (\Delta). 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 P-\Delta 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 P-\Delta concept, consider a column segment of length L subject to an axial force P and undergoing a relative sway \Delta 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.

P-delta effects 1 | DegreeTutors.com
Figure 1. Column segment with length L and axial force P experiencing relative sway \Delta between its ends.

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, \Delta. Again, notice that the axial forces, P are no longer co-linear. As a result, an extra second-order moment is developed at the base of the column,

(1)   \begin{equation*}M_2 = P\Delta\end{equation*}

If \Delta is a non-negligible value, we must also consider the additional secondary sway deflection caused my M_2. 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 P-\Delta 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, \Delta, 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 P-\Delta analysis with a simple example. Consider a simple column with one end fixed and the other free, subject to a vertical load P_V. We will assume that this column experiences a lateral load P_L at B that induces an initial first-order deflection \Delta.

P-delta effects 2 | DegreeTutors.com
Figure 2. Fixed-free column with axial load P_V and lateral load P_L which induces an initial sway deflection \Delta at its free end.

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)   \begin{equation*}P_{cr} = \frac{\pi^2EI}{4L^2}\end{equation*}

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

Let’s assume the column is subject to a compression force of P_V =0.7 P_{cr}. Let’s also assume that the column is subjected to a lateral force magnitude P_L = 0.05P_V.

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

(3)   \begin{equation*}\Delta_1 = \frac{M_BL^2}{3EI} = \frac{P_L L^3}{3EI}\end{equation*}

where M_B=P_L L is the moment generated at the base of the column by the lateral load P_L.

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

P-delta effects 3 | DegreeTutors.com
Figure 3. Second-order bending moment.

At this point we can make a simplification and assume that the second-order moment varies linearly between M_2 = 0 and M_2=P_V\Delta_1, Fig.4.

P-delta effects 4 | DegreeTutors.com
Figure 4. Fixed-free column with first-order, second-order and linearised second-order bending moment diagrams.

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 \Delta_1,

(4)   \begin{equation*}\Delta_{2\_1} = \frac{M_2 L^2}{3EI}\end{equation*}

where \Delta_{2\_1} is the first value of second-order deflection calculated. We refer to this value of \Delta_2 as the first iteration value, \Delta_{2\_1}.

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

The entire process must be repeated yet again with each subsequent value of \Delta_2 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)   \begin{equation*}\Delta_{total} = \Delta_1 + \sum_{i=1}^N \Delta_{2\_i}\end{equation*}

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.

P-delta effects 5 | DegreeTutors.com
Figure 5. Plot showing convergence of lateral column deflection arising from first and second-order effects.

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.

P-delta effects 6 | DegreeTutors.com
Figure 6. Lateral deflection normalised by column length (top), ratio of secondary to primary deflection (bottom).

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.

P-delta effects 7 | DegreeTutors.com
Figure 7. Ratio of secondary to primary deflection for increasing levels 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 P-\Delta 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 P-\Delta 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, P_L to calculate first-order sway deflection, \Delta_1 for the structure.
  2. Manually displace the structural model by applying the nodal displacements to each node in the model according to \Delta_1 nodal deflections.
  3. Apply any vertical loads to the deformed structure (P_V\times \Delta_1) to calculate second-order deflections, \Delta_{2\_1} and note the second-order nodal deflections obtained.
  4. Displace the structure (from its undisplaced position) according to \Delta_{2\_1} nodal deflections.
  5. Apply vertical loads again to the deformed structure (P_V\times \Delta_{2\_1}) to calculate second-order deflections, \Delta_{2\_2}.
  6. Repeat steps 4 and 5 until the second-order deflections observed have become negligibly small.
  7. Calculate the total deflection, \Delta_{total} = \Delta_1 + \sum \Delta_{2\_i}.

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 P-\Delta 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.

Tutorial by:

Dr Seán Carroll

Hi, I’m Seán, the founder and lead tutor at DegreeTutors.com. I hope your found whatever I wrote up there helpful! Feel free to get in touch or follow DegreeTutors on any of the social accounts.