The Intriguing Stresses in Pipe Bends

January 31, 2023

For many structural engineers, beam theory is a popular analysis tool. Using the equations can be beneficial when considering structural behavior, as they are easy to apply and provide useful results. However, it’s also known that due to the simplicity and convenience, beam theory is applied even when some underlying assumption may not hold up too well. This blog post investigates one such case, where there are severe pitfalls of using beam theory, and the true structural behavior is surprisingly different.

A Thought Experiment

Take a look at the pipe model shown below. It consists of a 90° pipe bend with adjacent straight legs and a moderate-to-thin wall thickness. The material is assumed to be isotropic and linear elastic. Now, imagine bending this pipe in the bend plane as indicated with the force arrows. While doing so, we encourage you to consider two questions:

  1. How does the pipe bend deform?
  2. Where do the highest stresses occur?

The geometry of a curve pipe (bottom) and its straight legs loaded in pure bending (top).
A curved pipe section with straight legs loaded in pure bending.

A view of the COMSOL Multiphysics UI showing the Model Builder with the Boundary Load 1 node highlighted and the corresponding Settings window with the Boundary Selection, Force, and Traction Field sections expanded.
Applying a moment on a boundary results in a varying surface traction. Such a load can be easily applied in the COMSOL Multiphysics® software as of version 6.1 using the Resultant load type in the Boundary Load feature.

Do you have answers in mind? If so, play the following animations to view the von Mises stress and the deformation for a 3D solid model loaded with a varying bending moment.

 

The von Mises stress distribution and deformation of a curved pipe loaded with a varying bending moment.

 

A detailed view of the von Mises stress distribution and deformation of the cross section at the center of the bend for a varying bending moment.

The pipe is slender with a constant cross section, so it would seem like a natural choice to treat such a structure as a beam in a simplified analysis. The bending moment is the only load acting on the structure, and it’s thus constant for any given section along the entire beam axis (both in the straight and curved pipe segments). Beam theory predicts that the axial stress has its peak values at the inner-most and outer-most fibers in the bend plane. But, defying beam theory, the opposite happens!

As shown in the animations, the von Mises stress distribution reaches its peak at the top and bottom of the pipe. The maximum stress even occurs on the inside of the pipe. The cross section of the bend also deforms significantly, and, more specifically, it ovalizes with the major axis either being oriented in the bend plane or perpendicular to it depending on the direction of the bending moment.

Analyzing Pipes with Beam Theory

Pipe bends are common in piping systems, which typically transport liquids or gas, often under high pressure. One place where you may find a lot of pipes are oil tankers. The labyrinthine piping systems can look pretty fascinating.

A close-up view of a red piping system on an oil/chemical tanker.
Piping systems on an oil/chemical tanker. Image licensed under CC BY-SA 3.0 via Wikimedia Commons.

Many piping standards (or codes) used for industrial applications are based on beam theory when it comes to the structural analysis. But, as we have already discovered, pipe bends generally do not behave like beams. When digging in piping standards, you will in fact find a lot of information dedicated to pipe bends. In particular, piping standards recommend to apply correction factors to the stiffness and stresses for curved pipe sections. (Ref. 1) In simplified terms, the bending stiffness from elementary beam theory, E\,I, should be divided by a so-called flexibility factor, or k-factor. Similarly, computed stresses, \sigma, should be multiplied by a stress intensification factor (SIF), or i-factor, when comparing with allowed stress values. In other words, the ordinary beam theory may still be applied, but bending stiffness and stresses should be corrected.

Piping standards provide estimates for these correction factors, which depend primarily on the geometry of the pipe bend. In some cases, it’s also recommended to apply additional corrections to account for the effect of the pressure difference, \Delta p, between internal and ambient pressure. This is because an internal overpressure tends to counteract the ovalization of the pipe cross section, which has a stiffening effect. Typical definitions for the k- and i-factors for a pipe bend subjected to in-plane bending are shown below (Ref. 2).

Stiffness correction Stress correction Flexibility characteristic
Correction:E \, I \rightarrow \, \frac{ E \, I }{ k }
 
k-factor:k = \textrm{max}\left(\frac{1.65}{h}; ~1 \right)
Correction: \sigma \rightarrow i \, \sigma
 
i-factor: i = \textrm{max}\left(\frac{0.9}{h^{2/3}}; ~1 \right)
A schematic of the flexibility characteristic.   h = \frac{R_\textrm{b} t}{r_\textrm{m}^2}

The flexibility characteristic, h, is a dimensionless geometric measure that takes into account the wall thickness, t; the mean pipe radius, r_\textrm{m}; and the centerline bend radius, R_\textrm{b}. The estimates for the correction factors are based on experimental data, or derived from higher-fidelity finite element analyses using solid or shell models. The correction factors can take rather large values, meaning that the deviation from ordinary beam theory is significant. It’s therefore crucial to consider the possible stress increase and stiffness reduction to get conservative results. The table below shows the correction factors for a range of pipe geometries. The wall thickness, t, and bend radius, R_\textrm{b}, are related to the outer pipe radius, R_\textrm{o}.

Geometry An image of the pipe geometry in column 2, according to its measurements in the table. An image of the pipe geometry in column 3, according to its measurements in the table. An image of the pipe geometry in column 4, according to its measurements in the table. An image of the pipe geometry in column 5, according to its measurements in the table. An image of the pipe geometry in column 6, according to its measurements in the table. An image of the pipe geometry in column 7, according to its measurements in the table.
\frac{R_\textrm{b}}{R_\textrm{o}} 4 4 4 2 2 2
\frac{t}{R_\textrm{o}} 20% 10% 2% 20% 10% 2%
k 1.7 3.7 20.2 3.3 7.4 40.4
i 1 1.5 4.7 1.4 2.5 7.6

The correction factors provided by piping standards have a long history, dating back as far as the 1950s. While these corrections are still used, they are only phenomenological models and do not explain what really happens in a pipe bend. As shown in the stress plots above, the true behavior of a pipe bend is a fully 3D problem. But before we further investigate the full 3D case, let’s recall some important relations from elementary beam theory.

Curved Beams

The classical beam equations are derived assuming that the beam axis is initially straight, and therefore the equations are, strictly speaking, only valid for straight beams. If the beam equations are instead derived from an initially curved beam axis, the definition of the strain along the axial direction differs compared to a straight beam, and subsequently the axial stress in the cross section differs. The assumed kinematics for the initially straight and curved beams are compared in the table below. The undeformed geometry is shown in gray, and the deformed geometry in red. The last row shows the analytic expressions for the axial stress resulting from an applied bending moment, M.

  Straight beam Curved beam
Kinematics for
pure bending
A schematic of the straight beam. A schematic of the curved beam.
Strain \varepsilon_{xx} = \frac{(\varrho_\textrm{m} + y) \, \textrm{d} \Theta – \textrm{d}x}{\textrm{d}x}\approx y \frac{\textrm{d} \Theta}{\textrm{d} x} \varepsilon_{xx}= \frac{\varrho \, (\textrm{d}\Theta + \textrm{d}\Phi) – r \, \textrm{d} \Theta}{r \, \textrm{d} \Theta}
Area moment
of inertia
I = \int_A \, y^2 \, \textrm{d} \bar{A} J = \int_A \, \frac{y^2}{1 + y/r_\textrm{c}} \, \textrm{d} \bar{A}
Stress \boxed{\sigma_{xx} (y) = \frac{M}{I} y} \boxed{\sigma_{xx} (y) = \frac{M}{A \, r_\textrm{c}} + \frac{M}{J} \, \frac{y}{1 + y/r_\textrm{c}}}

In the expressions above, \sigma_{xx} represents the stress in the beam’s axial direction. The local coordinate y has its origin at the centroid of the cross section. For the curved beam, r_\textrm{c} is the radius of curvature of the centroidal axis. Note that the curved beam relations recover the straight beam relations when the radius r_\textrm{c} goes to infinity.

According to the straight-beam assumption, the axial stress varies linearly across the section, whereas in a curved beam the stress distribution becomes nonlinear — even with the assumption that the cross section remains plane during deformation. For a straight beam, the neutral axis (where \sigma_{xx} = 0) coincides with the centroid (y = 0). For a curved beam, however, the neutral axis and centroid no longer coincide, as shown in the following example.

A figure showing the axial stress on the cut through of a beam with a solid, circular profile.
Cut through a beam with a solid, circular profile. The applied bending moments are indicated at the beam ends. The resulting axial stress, \sigma_{xx}, as computed with solid elements, ranges from tension (red) to compression (blue). In the curved section, \sigma_{xx} has a nonlinear distribution.

A 1D plot with normal stress on the x-axis and local coordinate y on the y-axis.
A comparison of the axial stress distribution, as computed with a full 3D solid model (blue), beam theory assuming a straight beam axis (red), and beam theory assuming a curved beam axis (green).

Even though this particular beam is by no means slender, curved beam theory matches a full 3D solid solution well. The stress comparison also emphasizes that straight beam theory predicts nonconservative results for highly curved beams. In this case, the maximum stress is underestimated by around 40 percent! However, for beams with a weak curvature, the use of straight beam theory is often still justifiable. As a rule of thumb, classical beam theory is often sufficient when the ratio of mean curvature radius to beam height is >10.

What is Different in Pipes?

We have now seen that a curved beam with a solid, circular cross section is well described with beam theory. So, the reason why certain pipe geometries behave so differently must lie in the wall thickness. As before, the following animations show a 90° pipe bend loaded in bending. In this case, the bending moment is kept constant, while the wall thickness is reduced.

 

The von Mises stress distribution (normalized) and principal stresses in a curved pipe with varying wall thickness loaded in pure bending.

 

Detailed view of the von Mises stress distribution (normalized) and cross section deformation in the center of the bend for a varying wall thickness. The inner-bend radius is on the left side.

For extremely thick walls (t/R_\textrm{o} > 35\,\%), the pipe still behaves as expected from beam theory. In fact, the von Mises stress is practically completely determined by the stress state discussed for the solid beam, that is \sigma_\textrm{mises} \approx | \sigma_{xx} |, where \sigma_{xx} is the stress parallel to the beam axis. The applied moment causes tensional stresses at the inner radius and compressional stresses at the outer radius. The fact that the stresses are slightly higher in the curved section compared to the straight legs is also fully consistent with curved beam theory as seen above.

\frac{t}{R_\textrm{o}} = 75\,\% \frac{t}{R_\textrm{o}} = 35\,\%
A close-up view of the pipe bend showing the Von Mises stress and principal stresses for a wall thickness of 75 percent. A close-up view of the pipe bend showing the Von Mises stress and principal stresses for a wall thickness of 35 percent.
\frac{t}{R_\textrm{o}} = 20\,\% \frac{t}{R_\textrm{o}} = 10\,\%
A close-up view of the pipe bend showing the Von Mises stress and principal stresses for a wall thickness of 20 percent. A close-up view of the pipe bend showing the Von Mises stress and principal stresses for a wall thickness of 10 percent.
\frac{t}{R_\textrm{o}} = 5\,\% \frac{t}{R_\textrm{o}} = 2.5\,\%
A close-up view of the pipe bend showing the Von Mises stress and principal stresses for a wall thickness of 5 percent. A close-up view of the pipe bend showing the Von Mises stress and principal stresses for a wall thickness of 2.5 percent.

View of the pipe bend showing the von Mises stress (normalized) and principal stresses for different wall thicknesses.

At a relative thickness of t/R_\textrm{o} \approx 35\,\% (which for any real application would be considered very thick), the stress distribution starts to change quite significantly. Additional perpendicular tensional stresses superpose the beam solution at the inner- and outer-bend radii. Simultaneously, the top and bottom of the pipe show compressional stresses. These additional circumferential stresses arise due to the ovalization of the cross section. Ordinary beam theory explicitly ignores such cross-sectional deformations, and it does fall short in capturing its effect.

Two side-by-side plots of the normalized strains in axial (left) and circumferential (right) directions at the center of a thick-wall bend.
Comparison of normalized strains in axial and circumferential directions at the center of a thick-walled bend. The ovalization (highly exaggerated) causes circumferential strains that have higher peak values than the axial strains.

Smaller wall thicknesses allow for more pronounced ovalization shapes, and this deviation from the initial circular profile becomes most extreme at the top and bottom. The cross-sectional deformations may be considered as local “bending deformations”, which cause strains in the circumferential direction of the pipe wall. At the inner-wall radius, where the local deformation is compared with a smaller, initial arc length, the strains (and thus stresses) reach their highest value.

Modeling Pipes in COMSOL

The COMSOL Multiphysics® software provides various features for analyzing pipes at different fidelity levels. Let’s take a look at some examples.

The Pipe Mechanics Interface

The Pipe Mechanics interface is a special version of the Beam interface. Its enhanced functionality makes it possible to take fluid properties into account and easily set up multiphysics couplings to include fluid and thermal effects.

The COMSOL Multiphysics UI showing the Model Builder with the Fluid-Pipe Interaction node selected, the corresponding Settings window, and a pipe mechanics model in both the Graphics and Plot windows.
A pipe mechanics model coupling structural, fluid, and thermal effects.

The Pipe Mechanics interface also has a built-in Bend feature that corrects stiffnesses and stresses by providing k-factors and SIFs, as discussed above.

A screenshot of the Settings window showing the Sketch and Correction Factors sections of the Bend feature.
The Settings window of the Bend feature, showing the stiffness and stress correction factors.

Submodeling Parts

While large parts of piping systems are generally well approximated with beam theory, there may be components for which higher-fidelity modeling using solid or shell elements is more appropriate. For such cases, it’s important to set up a correct kinematic coupling between displacement and rotational degrees of freedom (DOFs). A consistent coupling can be easily added using the Structure-Pipe Connection multiphysics coupling.

The COMSOL Multiphysics UI showing the Model Builder with the Structure-Pipe Connection node selected, the corresponding Settings window, and a structure pipe connection model in both the Graphics and Plot windows.
Coupling the Pipe Mechanics interface with the Solid Mechanics or Shell interface using the Structure-Pipe Connection multiphysics coupling.

Geometry Parts

COMSOL Multiphysics also provides a number of predefined geometry parts for 3D piping analyses, for both structural and fluid analyses. These built-in, fully parameterized parts can be found in the Part Libraries. The parts are optimized to be easily meshed, and they contain helpful predefined selections and work planes to easily connect components and build complex piping systems.

A close-up view of the many built-in geometry parts available in the Part Libraries, such as bend.
A section of a piping system created from built-in geometry parts.

The Part Libraries (left) provide many built-in geometry parts, such as fully parameterized piping parts. At right is an example of part of a piping system created from built-in geometry parts.

Further Reading

Learn more about the use of modeling and simulation for pipe design in these COMSOL blog posts:

References

  1. E.A. Wais and E.C. Rodabaugh, “Background of SIFs and Stress Indices for Moment Loadings of Piping Components”, United States, 2005; https://www.osti.gov/biblio/841246
  2. “Stress Intensification Factors (i-Factors), Flexibility Factors (k-Factors), and Their Determination for Metallic Piping Components,” ASME, 2017; https://www.asme.org/codes-standards/find-codes-standards/b31j-stress-intensification-factors-flexibility-factors-determination-metallic-piping-components/2017/drm-enabled-pdf

Comments (1)

Leave a Comment
Log In | Registration
Loading...
Mervin Blankenship
Mervin Blankenship
August 11, 2023

thanks for info

EXPLORE COMSOL BLOG