SWBPIPE · The open manual
Part II · FundamentalsII–1

Part II

Fundamentals of Piping Flexibility

A. Thermal expansion and restraint

A pipe that carries a hot fluid grows. If it is free to grow, it does so without stress. If it is not, the growth is turned into force, and the force into stress in the pipe and load on everything that holds it. Piping flexibility is the study of that conversion, and the design of lines that keep it within bounds.

This part teaches the discipline itself, written to what a correct calculation should give, whatever program makes it. Where SWBPIPE falls short, a mark headed Today links to the status of this edition. How SWBPIPE builds and solves its model is the subject of Part III.

1.Free growth.

A straight pipe of length L, warmed uniformly by ΔT and free at one end, lengthens by

δ = α ΔT L(1)

where α is the coefficient of thermal expansion. For carbon steel α is about 12 × 10−6 /K near ambient and rising with temperature; for austenitic stainless steel it is about 16 × 10−6 /K. A carbon steel run 30 m long, warmed by 150 K (270 °F), grows 54 mm (2.1 in). A stainless line of the same length grows about a third more.

A line of any shape, warmed uniformly and free, keeps its shape and grows in size, so its free end moves by αΔT times the straight line from the fixed point to that end, whatever route the pipe takes. The movement a line must absorb between two anchors depends only on where the anchors are. The loops between them decide how much force it takes.

2.Full restraint.

Hold both ends of the same straight run fixed and warm it. It cannot lengthen, so the anchors must push it back by the whole of its free growth. The mechanical strain is equal and opposite to the thermal strain, and the axial stress is

σ = −E α ΔT = −200 000 MPa × 12 × 10−6 /K × 150 K = −360 MPa(2)

compressive, about 52 ksi, with E the elastic modulus of carbon steel near ambient. That is well above the yield strength of ordinary carbon steel pipe. Notice what is absent from equation (2): the length, the diameter and the wall. A straight run between anchors reaches this stress whatever its size. For a DN200 line (219.1 mm outside diameter, 8.18 mm wall, metal area 5420 mm2) the anchor force is 360 MPa × 5420 mm2 ≈ 1.95 MN (about 440 kip). No anchor, nozzle or pipe rack is built for that. A heavier wall does not help: it raises the force in proportion and leaves the stress where it was.

3.Mean and instantaneous coefficients.

The coefficient of expansion is not constant. The instantaneous coefficient is the slope of the length–temperature curve at one temperature; the mean coefficient is the total strain from a reference temperature divided by the whole change. Growth is the mean coefficient times the change. The instantaneous value rises with temperature, so multiplying its hot value by the whole change overstates the growth. The B31 codes avoid the question by tabulating the total expansion from 21 °C (70 °F) to each temperature, in millimetres per metre (inches per 100 feet). With e(T) that total strain, the growth between any two temperatures is

δ = [e(T2) − e(T1)] L(3)

So a flexibility program should ask for temperatures, installed and operating, and look up the total expansion from them. Entering a coefficient invites the error above, and applying a mean coefficient measured from 21 °C to a change that starts elsewhere gives a second one.

4.Relief by bending.

A straight run cannot absorb its growth in direct stress, so the line must take it in bending. Turn the run through a right angle and the growth of each leg becomes a sideways deflection of the other. A pipe is far less stiff in bending than in direct stress, so the same movement is taken with a small fraction of the force. In an L-bend the short leg bends to absorb the growth of the long one; in a Z-bend the offset leg absorbs the growth of the two parallel legs; a U-loop set into a long run absorbs the growth on either side by the bending of its two legs.

5.The guided cantilever.

A leg that absorbs a sideways movement Δ at one end, with both ends prevented from rotating, bends as a guided cantilever. The force and moment at each end are

F = 12EIΔ/L3,   M = 6EIΔ/L2(4)

with I the second moment of area and L the length of the leg. Divide the moment by the section modulus Z = I/(Do/2), and

σ = M/Z = 3EDoΔ/L2(5)

The force falls with the cube of the leg length and the stress with its square. Doubling a leg divides its stress by four and its end force by eight. The wall has again dropped out of the stress, and the diameter has come in on top: a larger pipe needs a longer leg for the same movement, in proportion to the square root of its diameter.

Take the DN200 line of section 2, 30 m long, whose 54 mm of growth is taken by a leg at right angles to it. The section has I = 30.19 × 106 mm4 and Z = 275.6 × 103 mm3. With E = 200 000 MPa:

Guided cantilever, DN200, Δ = 54 mm
Leg length LEnd force FEnd moment MBending stress M/Z
6 m18.1 kN54.3 kN·m197 MPa
8.4 m6.6 kN27.7 kN·m100 MPa
12 m2.3 kN13.6 kN·m49 MPa

Set the 18 kN of the 6 m leg beside the 1.95 MN of the straight run between anchors. The movement is the same; the shape has cut the force a hundredfold.

The guided cantilever is a sizing tool. It ignores the rotation of the corners and the extra flexibility of the bends, both of which relieve the leg, so it usually overstates the force. It gives the nominal stress in straight pipe and says nothing of the higher stress at the bend, which is the business of section E. Use it to choose a layout; analyse the layout you choose.

6.Flexibility, not strength.

Most structural design is a matter of strength: find the load, then make the member strong enough to carry it. Thermal loading does not work that way. The load is whatever it takes to impose a movement, and it grows with the stiffness of the thing moved. Make the pipe stronger and it pushes harder. The design problem is one of flexibility: to arrange the line so that it takes its movements with forces and stresses that the pipe, its supports and its equipment can bear.

For the student

Equation (5) is worth working backwards. To keep the nominal bending stress in the DN200 leg to 100 MPa, L = √(3EDoΔ/σ) = √(3 × 200 000 × 219.1 × 54/100) mm ≈ 8.4 m. A layout engineer who knows this can place a loop on a plot plan before any analysis is run.

Part II · FundamentalsII–2

B. Primary and secondary loads; pressure

Loads on a pipe are of two kinds, and the whole structure of the piping codes rests on telling them apart. They differ not in size but in what happens when the pipe begins to yield.

1.Primary loads.

Weight and pressure are primary loads. They are forces, and they do not diminish when the pipe deforms. If the stress they cause exceeds what the material can carry, the pipe keeps deforming until it collapses. Primary loads are not self-limiting, and the codes limit them against gross deformation, with the allowable stress at the operating temperature.

2.Secondary loads.

Thermal expansion and anchor movements are secondary loads. They are imposed movements, and the stress they cause is only what it takes to make the pipe conform to them. If it reaches yield somewhere, the pipe yields locally, the movement is accommodated, and the stress stops rising. The load is self-limiting; one application cannot cause collapse. What it can cause is fatigue. Each time the line heats and cools, the stress at its most highly strained points cycles, and enough cycles of a large enough range crack the pipe, most often at a weld, a bend or a branch. So the codes limit secondary stress by its range, with an allowable that depends on the number of cycles.

3.Self-springing and relaxation.

Suppose the calculated hot stress at some point is above yield. On first heating the point yields, and the hot stress settles lower than the elastic calculation said. On cooling, the pipe returns elastically through the full range, and because it was yielded hot it now carries a stress of the opposite sign cold. The line has sprung itself. At high temperature, creep does the same work more slowly: the hot stress relaxes over time and reappears reversed when the line is cold.

After the first cycle or few, a line whose calculated range is within the bounds the codes set cycles elastically between a hot and a cold stress whose difference is the range. This is shakedown. The hot stress and the cold stress are each uncertain, but their difference is not; it is fixed by the movement.

4.Pressure stresses.

Pressure p in a pipe of outside diameter Do and wall t sets up a hoop stress around the wall,

σH ≈ pD/2t(6)

with D the mean diameter for the average stress through a thin wall, or the outside diameter for the conservative form. Wall thickness for pressure is designed against the hoop stress, by the pressure-design rules of the code, before flexibility is considered.

The pressure acting on the closed ends of a run, or on the equivalent area at a bend, is carried by the wall as a longitudinal stress. The exact mean value, the end force pπDi2/4 over the metal area, is

σL = pDi2/(Do2 − Di2)   or, simply and slightly conservatively,   σL = pDo/4t(7)

The simpler form has long been used in the B31 codes; which form applies in a sustained stress calculation is set by the code you work to. For the DN200 line at 4 MPa, the two give 23.8 and 26.8 MPa. The longitudinal pressure stress is a primary stress and belongs in the sustained stress, added to the bending stress of weight.

5.Where the pressure force goes.

In a continuous line the closed-end force is carried by the wall and is in equilibrium with itself. At an elbow, the pressure on the outer part of the bend is balanced by the tension in the wall of the two legs; the elbow is not pushed off its supports. Supports and anchors on a continuous welded line do not see the closed-end pressure force.

The force reaches the supports only where the line is not continuous in tension: at an untied expansion joint, whose bellows cannot carry it and which transmits the pressure times its effective area to the anchors on either side; at an open end or a discharge; at a slip joint or a sleeve coupling. At such points the pressure thrust is a real external load, often the largest on the anchor, and it must be applied there. Elsewhere it must not be.

6.Pressure elongation.

A pressurised pipe free to move changes length. The hoop stress stretches the circumference and, through Poisson’s ratio ν, shortens the pipe; the longitudinal stress lengthens it. The net axial strain is

εL = (σL − νσH)/E ≈ (0.5 − ν)σH/E = 0.2σH/E   for ν = 0.3(8)

taking σL = σH/2 for a thin wall. Programs often call this the Bourdon effect, after the curved gauge tube that straightens under pressure; bends do the same, opening slightly. For the DN200 line at 4 MPa, σH ≈ 52 MPa and the strain is 5.2 × 10−5: 1.6 mm over 30 m, beside 54 mm of thermal growth. It is small in most plant piping, but not on long high-pressure lines or large thin pipe. Applying the whole closed-end force pAi as an axial load, without the Poisson contraction, overstates the elongation two and a half times.

7.A fully restrained, heated, pressurised pipe.

Hold a straight run so that it cannot change length at all, as a long buried line is held by the soil, and warm and pressurise it. The axial strain is zero, so the mechanical strain must cancel the thermal strain: (σL − νσH)/E + αΔT = 0, which gives

σL = νσH − EαΔT(9)

Pressure puts the wall in longitudinal tension and temperature puts it in compression. In the wall stress they act in opposite senses and do not add. For the DN200 line, 0.3 × 52 − 360 ≈ −344 MPa. The restraint force, measured from the free state, is another matter: both pressure elongation and thermal growth would lengthen the free pipe, so both add to the push on whatever holds it. A program must get both right, and it can only do so by treating pressure as a state of the wall with its Poisson strain, not as an axial force.

Today

SWBPIPE applies pressure as an axial force pAi on each pressurised element, without the Poisson term, and leaves pDo/4t out of the straight-pipe stresses. Status.

Part II · FundamentalsII–3

C. Load cases

A flexibility analysis is not one solution but several, each a state the line will be in, and each checked against a different limit. The load cases follow from the distinction of section B: primary loads are checked as they act; secondary loads are checked as a range between states.

1.The basic cases.

A line may have several operating conditions: a normal state, an upset, a steam-out, a start-up, a cold winter shutdown. Write them as temperature states T1 … Tn and pressure states P1 … Pn, with W the weight of pipe, contents, insulation and components, and D any anchor movements that go with a temperature state.

CaseLoads, and what it is for
SustainedW + P1, with spring hangers carrying their design loads. Checked against the sustained allowable (B31.3 clause 320).
OperatingW + T1 + P1 (+ D1), and one such case for each operating state. Gives the hot displacements, the hot support and nozzle loads, and which supports lift off, close or slide. Not itself a code stress case.
ExpansionThe difference between two states: operating less sustained, or one operating state less another. Gives the displacement stress range (section D).
OccasionalWind, earthquake, relief-valve thrust and the like: primary loads of short duration, solved as their own cases, each acting in both senses, and added to the sustained result with the sign that makes the sum larger. Checked against a higher allowable than the sustained one (B31.3 clause 302.3.6, 1.33Sh). Section K returns to them.
HydrotestWater-filled weight and test pressure, with spring hangers locked. Gives support loads and stresses during the test.

2.Expansion as a difference of states.

The expansion case is not a set of loads. It is the change in forces, moments and displacements between two solved states. With one operating condition it is the operating case less the sustained case: what changed when the line went from cold to hot. With several, each pair of states gives a range, and the one that governs is usually that between the two most distant states, for example a hot operating condition and a cold condition below the installed temperature. It is formed directly between those two states, not by adding the stress ranges each makes with the installed condition.

A difference of two states is not an envelope. An envelope takes, for each quantity at each point, the largest value found in any of several cases. The values it collects may come from different cases and do not describe any one state the line is ever in. That is useful for sizing a support, which must carry the worst load in any case. It is wrong for a stress range, which must be formed between two real states, quantity by quantity with its sign.

Today

SWBPIPE has no sustained, expansion or occasional case types; combinations are formed by factors, subtraction and envelopes that you set up. It gives the forces, moments and stress components, not the code stresses. Status.

3.Superposition and its limits.

For a linear system, the solution for several loads acting together is the sum of the solutions for each alone. Solve W, T and P separately and add, and the result is exact. That is what makes the difference OPE − SUS the same as the thermal case alone.

A line with supports that lift off, gaps that close or rests that slide is not linear. Whether a rest carries load depends on everything acting on it. Solved alone, a thermal case has no weight on the rests, so it lets a rest lift at the first upward movement, or holds the pipe down where no rest can; neither is the truth. So each state the line is actually in is solved as a whole, with all its loads together: the operating case as W + T + P in one nonlinear solution, the sustained case likewise. The expansion range is then the difference of those two complete solutions.

The sustained case needs one more thought. Where a rest lifts off when the line is hot, the weight it carried cold is carried by its neighbours for as long as the line is hot. That hot distribution of weight is also a sustained condition and must also meet the sustained limit. It is found by solving the sustained loads in the operating support state, with the lifted supports removed, as well as in the cold state, and taking the worse.

Today

SWBPIPE adds cases solved with nonlinear supports as if the system were linear, without warning. Solve each combined loading as a case of its own. Status.

4.Hydrotest.

A line is tested before it is put in service, usually with water at a pressure above design (for B31.3, at least 1.5 times the design pressure, adjusted for temperature). A gas or steam line designed for a light fluid may never have been supported for a pipe full of water. The test case uses the water-filled weight, the test pressure and the test temperature, with spring hangers locked by their travel stops so that they act as rigid supports. It checks the supports, the stresses under test, and whether temporary supports are needed.

Today

SWBPIPE has no hydrotest case: a hydrotest pressure stops the solve, and hangers cannot be locked. Status.

Part II · FundamentalsII–4

D. The displacement stress range

The rule that the range of thermal stress governs, and not its peak, is the foundation of modern piping flexibility. Rossheim and Markl proposed an allowable stress range in 1940, in “The Significance of, and Suggested Limits for, the Stress in Pipe Lines Due to the Combined Effects of Pressure and Expansion” (Trans. ASME, vol. 62), and the rule rests on the behaviour described in section B.

1.Why the range governs.

A secondary stress cannot cause collapse; its danger is fatigue, which is driven by the cycle. Self-springing and relaxation shift the hot and cold stresses but leave the difference between them, which is fixed by the movement. If the range is kept within the shakedown limit, roughly the sum of the yield strengths at the two extremes of the cycle, the line settles into elastic cycling, and the range and the number of cycles decide its life.

Markl then measured the life. The same 1940 work began a programme of fatigue tests on piping components, reported in “Fatigue Tests of Welding Elbows and Comparable Double-Mitre Bends” (Trans. ASME, vol. 69, 1947) and “Fatigue Tests of Piping Components” (vol. 74, 1952). Components were cycled in bending to failure and compared with the fatigue data for straight commercial pipe containing a butt weld. That comparison gave the stress intensification factors of section E, and the shape of the fatigue curve gave the cyclic reduction below. Markl’s “Piping-Flexibility Analysis” (Trans. ASME, vol. 77, 1955) brought the method together, and it entered ASA B31.1-1955.

2.Forming the range.

The range is formed from the forces and moments, not from stresses. Take the two states that bound the cycle. At each point, subtract the moments of one state from those of the other, component by component, with their signs: the in-plane moment range, the out-of-plane moment range, the torsion range, and the axial force range. Then compute the stress from those ranges.

Never difference a summary or combined stress. A combined stress is built from squares and absolute values; the signs and directions of its parts are gone. Suppose the in-plane moment at an elbow is +10 kN·m hot and −10 kN·m cold. The moment range is 20 kN·m. The combined stress is the same in both states, and differencing it gives a range of zero. A moment that turns from one plane to another does the same thing less obviously.

Where there are several thermal states, the range is formed between each pair that makes a real cycle, and the pair with the largest range is checked against the allowable. Smaller cycles of other sizes are counted into the number of cycles, as described below.

3.The computed range in B31.3.

The current form of B31.3 (clause 319.4.4) computes the displacement stress range as

SE = √[(|Sa| + Sb)2 + (2St)2](10)

where Sa is the axial stress range from the axial force range, intensified by an axial SIF where one applies, and the bending and torsional ranges are

Sb = √[(iiMi)2 + (ioMo)2]/Z,   St = itMt/2Z(11)

with Mi, Mo and Mt the in-plane, out-of-plane and torsional moment ranges, ii, io and it the corresponding SIFs, and Z the section modulus. Earlier editions omitted the axial term and wrote SE = √(Sb2 + 4St2). B31.1 forms the range differently in detail: a single SIF applied to the resultant of all three moment ranges, torsion included.

Flexibility calculations use nominal dimensions (B31.3 clause 319.3.5): the stiffness of the line, and the section modulus for the range, come from the nominal wall, not a wall reduced by corrosion allowance or mill tolerance. The sustained stress, by contrast, is computed on a section reduced by the allowances.

4.The allowable range.

The allowable displacement stress range (B31.3 clause 302.3.5) is

SA = f (1.25Sc + 0.25Sh)(12)

with Sc and Sh the basic allowable stresses at the minimum and maximum metal temperatures of the cycle. Where the sustained stress SL is less than Sh, the unused part may be added, which gives the liberal form

SA = f [1.25(Sc + Sh) − SL](13)

The logic is that of shakedown. The allowable is built to keep the range within about the sum of the cold and hot yield strengths; the sustained stress, which is present at the hot end of the cycle, uses up part of that margin, and what it leaves is available to the range.

The factor f is the cyclic stress range reduction factor. It comes from Markl’s fatigue curve and falls as the number of equivalent full cycles N rises:

f = 6.0 N−0.2,   not more than 1.0(14)

so f is 1.0 up to about 7000 cycles, roughly one a day for twenty years, and falls below it for lines that cycle more often. Cycles of smaller range are counted into N by weighting each with the fifth power of its ratio to the largest range. Recent editions of B31.3 let f rise a little above 1.0 for some lower-strength steels at moderate temperature. B31.1 sets out the same framework, with its stress checks under clause 104.8.

5.The modulus for the range, and hot reactions.

The range is computed with the reference modulus Ea, the modulus at the installed temperature, whatever the operating temperature. This is not an oversight. The allowable range and the SIFs come from fatigue tests at room temperature, their results expressed as elastic stresses at the room-temperature modulus, and the calculated range must be expressed on the same basis.

The forces and moments the line exerts hot are smaller, because the hot metal is less stiff. B31.3 clause 319.5.1 estimates the reactions of a two-anchor line without intermediate restraints from the range R computed at Ea:

Rm = R (1 − 2C/3)(Em/Ea)(15)

with Em the modulus at the extreme (usually hot) temperature and C the cold-spring factor of section G, zero when there is no cold spring. The cold reaction is estimated from the cold spring, or from the self-springing the line will undergo if its range is large, whichever is greater. A general analysis computes the reactions directly, but the principle is the same: the stress range at Ea, the hot loads at the hot modulus.

Today

SWBPIPE takes temperatures pipe by pipe as a change, not as operating states for a line, and solves each case with one set of properties; a thermal case solved at the hot modulus gives a range the codes compute with Ea. Status. It computes no SE or SA. Status.

Part II · FundamentalsII–5

E. Bends, SIFs and branch connections

A pipe bend is more flexible than a straight pipe of the same length, and more highly stressed. Both facts come from the same cause, and both have to be in any flexibility calculation that is to be believed.

1.Ovalization.

Bend a curved tube and its cross-section does not stay round. The longitudinal stresses in the wall, following the curve, have components toward the middle of the section, and they flatten the section into an oval. A flattened section has a smaller second moment of area, so the bend rotates more under a given moment than beam theory predicts, and the wall carries high local bending stresses around its circumference.

Bantlin observed the ovalization in 1910, and found bends much more flexible than the theory of bars predicts. Von Kármán gave the first theory in 1911, for in-plane bending, and his analysis is the beginning of the theory of curved pipe. Vigness extended it to out-of-plane bending in 1943.

2.The flexibility factor.

The effect on stiffness is carried by the flexibility factor k: the ratio of the rotation of the component to that of a straight pipe of the same length under the same moment. A bend with k = 7 bends as seven times its own arc length of straight pipe would. In a calculation, the bend’s contribution to the flexibility of the line is multiplied by k.

For a bend of radius R, with pipe of mean radius r and wall t, the governing parameter is the flexibility characteristic or pipe factor

h = tR/r2(16)

Thin walls and tight bends make h small, and the bend flexible. The classical result for the flexibility factor is the asymptotic closed form of Clark and Reissner (“Bending of Curved Tubes”, Advances in Applied Mechanics, vol. 2, 1951):

k = 1.65/h,   not less than 1(17)

These analyses all give the same flexibility factor for in-plane and out-of-plane bending, and one k for both planes is the classical practice.

3.The stress intensification factor.

The stress intensification factor, or SIF, is defined by fatigue test, not by stress analysis. It is the ratio of the moment that fails straight pipe with a girth butt weld to the moment that fails the component in the same number of cycles. It therefore takes in the local bending of the wall, the geometry and the welds of the component, and it is measured relative to a butt weld, whose SIF is 1.0 by definition. Markl’s fatigue tests gave the classical forms for a bend:

ii = 0.9/h2/3,   io = 0.75/h2/3,   each not less than 1(18)

for in-plane and out-of-plane bending. In-plane bending acts in the plane of the bend, opening or closing it; out-of-plane bending acts across that plane. The two are different because the ovalization, and so the stress, is different, and a program must keep the moments apart so that each can be multiplied by its own factor.

A long-radius elbow in the DN200 line makes the numbers concrete. With R = 304.8 mm, t = 8.18 mm and r = 105.5 mm, h = 0.224, so k ≈ 7.4, ii ≈ 2.44 and io ≈ 2.03. The elbow bends like 7.4 times its arc length of pipe, and fails in fatigue under a moment range less than half of what fails the butt weld beside it.

SIFs multiply moments, and only moments (with, in recent practice, axial force where a factor for it is given). They are never applied to the pressure stress, which has its own rules.

4.Pressure, flanges and the modern source.

Internal pressure resists ovalization. In large, thin-walled bends at high pressure it stiffens the bend and lowers both k and i appreciably, and the codes have long allowed a correction for it. Flanges or other stiff components welded close to the ends of a bend also restrain the ovalization, and lower k and i for the same reason.

For many years the B31 codes carried their flexibility and stress intensification factors in an appendix table. They now take them from ASME B31J, whose 2017 edition gives in-plane, out-of-plane and torsional factors for bends, tees and branch connections, with test methods for components it does not cover. The classical forms above are the right place to learn the subject; the values for design come from the standard your code calls up.

5.Branch connections.

At a branch, three pipes meet. The simplest model joins them rigidly at the point where their centrelines cross. That is stiffer than the truth. The wall of the run pipe deforms locally under the branch’s moments, like a shallow shell, and the branch rotates a little relative to the run. In a thin-walled run this local flexibility can reduce the branch moments considerably; it is modelled as a rotational spring at the surface of the run, with flexibility factors for branch connections.

The branch SIF and the header SIF are applied to the moments in the branch and in the run respectively. Stress on the branch side of a reduced branch is computed with a section modulus belonging to the branch, not the run; B31.3 has long used an effective section modulus Ze for it. The SIF and its section modulus come as a pair from the same source, and must be used together.

Today

A bend drawn in SWBPIPE is analysed as the straight chord between its ends, with no flexibility factor or SIF; the curved element, which takes the k you enter, can be reached only from a model file. Status. The SIF review row scales a summary stress rather than the in-plane and out-of-plane moments. Status. A tee is a rigid junction, with no branch section modulus and no local flexibility. Status.

Part II · FundamentalsII–6

F. The elastic centre and the force method

Before computers, flexibility calculations were made by hand by the force method, and the classical manuals set it out for lines of every shape. The 1941 Kellogg manual, Design of Piping Systems, written for The M. W. Kellogg Company by D. B. Rossheim, A. R. C. Markl, H. V. Wallstrom and E. Slezak, is the classic of that tradition and the model for this one. The method remains the clearest way to see how a line absorbs its growth, and the best check on a program.

1.Release and restore.

Take a line between two anchors, A and B. Release anchor B. The line is now a cantilever from A, statically determinate, and under a uniform temperature rise it simply grows: by section A, B moves by αΔT times the straight line from A to B. Anchor B must put it back. The forces and moments it applies to do so are the redundants: three in a plane (two forces and a moment), six in space.

Let Xj be the redundants and mj(s) the bending moment along the line due to a unit value of each. The movement at B in the direction of redundant i caused by a unit value of redundant j is the flexibility coefficient

fij = ∫ mi mj ds/EI(19)

the integral taken along the whole line, with each bend’s arc length multiplied by its flexibility factor k. Axial and shear deformations are small in ordinary lines and are usually left out; in space, torsion adds a term in ds/GJ. Compatibility at B then reads

∑j fij Xj = −Δi(20)

with Δi the free thermal movement of B relative to A, less any movement imposed on B relative to A. Solve these three or six equations, and every force and moment in the line follows by statics.

2.The elastic centre.

The equations are coupled: each redundant appears in every equation. The coupling can be removed by moving the point at which the redundants act. Think of each element of the line as carrying an elastic weight dw = ds/EI, with a bend’s weight multiplied by k. The elastic centre is the centroid of that weight:

W = ∫ dw,   x̄ = ∫ x dw/W,   ȳ = ∫ y dw/W(21)

Tingey introduced the idea into piping analysis, as the “virtual centre of gravity” of the line.

Attach a rigid arm from B to the elastic centre and let the redundants act there, as a force (Hx, Hy) and a moment M0. Measure x′, y′ from the elastic centre and write Ixx = ∫y′2dw, Iyy = ∫x′2dw and Ixy = ∫x′y′dw, the moments of inertia of the elastic weight. Because the first moments about the centroid vanish, the moment equation separates from the force equations, and for a plane line in uniform expansion, with e = αΔT(rB − rA) the free movement of B,

M0 = 0,   HxIxx − HyIxy = −ex,   HyIyy − HxIxy = −ey(22)

The first result is the classical one: the thermal thrust of a two-anchor line is a single force whose line of action passes through the elastic centre. The bending moment at any point is that force times the point’s distance from its line of action, so the largest moments are at the points farthest from it. Take the force axes along the principal axes of the elastic weight, where Ixy = 0, and the two force equations separate as well; each redundant is then found alone.

In space the same construction holds with six redundants and the moments of inertia of the elastic weight in three dimensions, with torsion entering through ds/GJ. It is lengthy by hand.

3.A worked case, and the stiffness method.

The verification case Expansion loop with a curved bend solves a plane L-bend between two anchors by exactly this method: anchor B released, the flexibility coefficients integrated in closed form with the arc’s share multiplied by k, and the three equations solved together for k = 1, 5, 10 and 20. Moving the redundants to the elastic centre would separate the moment from the forces and give the same reactions. The hand result then checks the solver.

SWBPIPE, like every modern program, solves the same problem by the stiffness method (Part III): unknown displacements at the nodes, rather than unknown forces at a released anchor. The two methods give the same answer for the same model. Its curved bend element is built from the same flexibility integrals, evaluated around the arc and then inverted to a stiffness.

For the student

Sketch a line between two anchors and mark its elastic centre. The thrust passes through it. Its direction follows from equation (22); it lies along the line joining the anchors only when that line is a principal axis of the elastic weight. The elbows farthest from the line of thrust carry the largest moments, and with their SIFs are where to look for the highest stress. Doubling a bend’s k pulls the elastic centre toward it; a reader who sees why understands the method.

Part II · FundamentalsII–7

G. Cold spring and hot and cold reactions

Cold spring is the deliberate pre-straining of a line when it is erected, so that part of its thermal movement is taken up before it is ever heated. A line that will grow is cut short and pulled together at the last joint; a line that will shrink, one operating below the installed temperature, is cut long and pushed together.

1.What it does.

Let the expansion range produce a reaction R at an anchor. Cut the line short by a fraction C of the movement it must absorb, and close the gap at erection. The cold line now carries a reaction of about CR at its anchors, opposite in sense to the hot reaction, and when it is heated the reaction swings through the full range to about (1 − C)R hot. With C = ½, the hot and cold reactions are about equal and half the uncorrected hot value. The usual purpose is to bring down the hot load on a nozzle or anchor that cannot take the full value.

2.No credit in the range.

Cold spring moves both ends of the cycle by the same amount. It changes the hot and cold stresses but not the range between them, and the range is what governs fatigue. So the codes give no credit for cold spring in the displacement stress range (B31.3 clause 319.2.4); the range is computed as if there were none.

3.Two-thirds credit in the reactions.

Cold spring is hard to achieve exactly and hard to verify once the line is closed. The codes therefore credit only part of it in the hot reactions: equation (15) of section D takes two-thirds of C, giving Rm = R(1 − 2C/3)(Em/Ea), while the cold reaction is taken at the full CR. Both are estimates for two-anchor lines; a general analysis applies the cold spring explicitly and applies the same caution to the result.

4.Modelling it.

In a model, cold spring is an imposed relative displacement across the cut: two nodes at the location of the closing weld, a gap between them equal to the cut, and the gap closed as a load. Once closed, the cut is present in every state, cold and hot, so it cancels from the range between them, as the codes require. It is a displacement load, not a primary one, and the sustained case leaves it out. The cold-spring load alone shows the force needed to close the joint at erection.

For the practitioner

A line whose layout depends on cold spring to pass is fragile. The cut must be made to the right length, in the right direction, and closed without overheating or distortion, and none of that can be seen afterwards. Where a small change of route gives the flexibility instead, it is usually the better answer.

Today

SWBPIPE cannot yet apply imposed movements, so cold spring cannot be modelled. Status.

Part II · FundamentalsII–8

H. Supports and restraints in practice

A line is held by supports, and each support acts in some directions and not in others. A model is only as good as its description of what each support actually does: which way it holds, whether it can let go, how much it gives, and whether the pipe slides on it.

1.The kinds of restraint.

SupportWhat it restrains
AnchorAll six: three translations and three rotations. Divides a system into parts that can be analysed apart.
Rest, or +Y supportDownward movement only. It carries weight and lets the pipe lift off. Most rack supports are rests.
GuideMovement across the pipe axis, usually the horizontal direction perpendicular to the pipe, while letting it slide along its axis.
Line stopMovement along the pipe axis, to direct growth toward a loop or away from a nozzle.
Limit stopAxial movement beyond a set clearance, in one or both directions.
Hanger (rod)Downward movement, from above. A rod can pull but not push.

Any of these may have a gap: a clearance through which the pipe moves freely before the restraint takes hold. A gap may be on one side only, as with a rest, or on both sides, as with a guide with clearance on each side of the shoe. The two behave quite differently and must be described separately.

Supports do not always act along global axes. A guide on a run that is not parallel to the plant axes acts perpendicular to the pipe, not along north or east. A program must allow restraints along any direction, and restraints defined relative to the pipe, so that a guide stays a guide whatever way the pipe runs.

Today

SWBPIPE restraints act along global axes only: there are no skewed or pipe-relative restraints and no double-acting gaps or limit stops, and only a spring takes a stiffness. Status.

2.Stiffness.

A support is usually modelled as rigid. A rigid support draws more load than a soft one, so this is conservative for the support itself, but it passes less load to its neighbours than they will see in the plant. A cantilevered bracket, a long rack beam or a vessel shell can move appreciably under load; where it matters, the support is given its real stiffness.

3.Friction.

A pipe that slides on a rest is resisted by friction. The friction belongs to the rest itself and depends on that rest’s own normal reaction N, which changes with the load case. It acts in the plane of the support, against the sliding, and it is Coulomb friction: the pipe sticks until the resultant tangential force reaches μN, then slides against a force of that size opposing the direction of slip. With T1 and T2 the two tangential components,

√(T12 + T22) ≤ μN(23)

The limit is a circle, not a square: a pipe sliding diagonally is resisted by μN, not by μN in each direction. Typical coefficients are about 0.3 for steel on steel and about 0.1 for PTFE slide plates.

Friction usually raises the operating loads on anchors and nozzles, because it resists the growth that the line was laid out to accommodate. It can also lower them, by holding a leg that would otherwise push on a nozzle. So the operating case is commonly run both with and without friction, and the worse result taken for each quantity.

Friction dissipates energy, so its result depends on the path. A static solution assumes the loads were applied in one step from a stated starting state; on cooling the friction reverses, and the line does not return exactly to where it began. Treat friction results as estimates, not as one exact answer.

Today

SWBPIPE models friction as a separate support on one global axis, whose normal force must come from a two-way restraint; it cannot sit on a rest that may lift off, and sliding in two directions follows a square law. Status.

4.Lift-off.

A rest carries weight only while the pipe presses on it. As a line heats, parts of it rise, and a rest there stops carrying load. Its share of the weight goes to its neighbours. The same happens with a rod hanger: a rod can hold a pipe up but not hold it down, and when the line grows upward the rod goes slack.

This is the classic way nozzle loads go wrong. A pump discharge line rises from the nozzle, turns and runs toward the rack, with a rod hanger or a rest near the pump carrying the weight of the valves and the riser. When the line heats, the riser grows upward from the nozzle and lifts the line off that support. The weight it carried now has nowhere to go but the next support along and the pump nozzle. An analysis that treats the support as two-way has it holding the pipe down, which no rest or rod can do; the plant shows a misaligned coupling. The remedy is a spring, which section I takes up.

Lift-off makes the analysis nonlinear. A program finds the support states by iteration: assume a state for each support, solve, check each support against the result, change those that disagree, and solve again until none changes.

Today

SWBPIPE stops this iteration after four passes, which is too few for a long rack line. Status.

Part II · FundamentalsII–9

I. Spring and constant hangers; hanger design

Where a line must be supported at a point that moves vertically between cold and hot, a rigid support will either lift off or hold the pipe against its growth. A spring carries the weight and lets the point move. The price is that a spring’s load changes with its travel, and the change goes somewhere.

1.The procedure.

Spring hangers are designed, not guessed, by a standard procedure in five steps.

Step 1. Restrained weight. Solve the weight case, W, with a rigid vertical support at each hanger location. The load found at each is the hot load: the share of weight the spring is to carry when the line is hot. A spring that carries exactly this in the hot position leaves the rest of the line, and the equipment, with the weight distribution of a rigidly supported line.

Step 2. Travel. Solve the operating case with each hanger removed and its hot load applied as an upward force in its place. Under weight alone this force holds the pipe exactly where the rigid support did, so the movement found at the hanger is the thermal movement: the travel. Removing the hanger, or giving it zero stiffness, is what lets the travel appear.

Step 3. Rate and variability. A spring of rate k changes its load by k times its travel δ. The variability is that change as a fraction of the hot load,

V = k|δ| / Fhot(24)

and it is commonly held to 25 %, following MSS SP-58 practice, or lower near sensitive equipment. Choose a rate low enough to meet the limit, and a spring whose working range contains both the hot and the cold load. If no variable spring will do, use a constant-effort support.

Step 4. Cold load. The load the spring is set to at installation, the cold or installed load, follows from the hot load and the travel. A spring’s load falls as the pipe rises, because the spring extends, and rises as the pipe falls. So for travel measured from cold to hot, positive upward,

Fcold = Fhot + kδup(25)

a point that moves up is set cold at more than its hot load; a point that moves down is set cold at less. It is worth stating in words on every hanger schedule: a sign error here doubles the load change the spring was chosen to limit, and puts it on the equipment.

Step 5. Final analysis. Put each spring into the model as its rate with its cold load applied as a preload, and run every load case again. In the operating case each spring should now carry close to its hot load. In the sustained case it is common practice to carry each spring as a force equal to its hot load, the weight share it was designed for; the difference between cold and hot load is caused by thermal movement, and belongs with the displacement stresses.

2.A worked variability check.

A hanger location has a hot load of 12.0 kN and moves 20 mm upward from cold to hot. For 25 % variability the rate may not exceed 0.25 × 12 000 N / 20 mm = 150 N/mm. Choose a spring of 120 N/mm (about 685 lbf/in). Its variability is 120 × 20 / 12 000 = 20 %, and its cold load is 12.0 + 0.120 × 20 = 14.4 kN. The 2.4 kN it carries cold above its hot load is taken off the neighbouring supports and nozzles when the line is cold, and handed back to them as it heats. The variability limit is a limit on that exchange.

3.Constant-effort supports.

A constant-effort support uses a spring and a lever or cam so that its load stays nearly the same through its whole travel. It is chosen where the travel is too large for a variable spring to stay within the variability limit, or where the load on nearby equipment must not change at all. It is set to the hot load. It carries a fixed force in every case that includes weight, and contributes nothing to a thermal range: its load does not change when the line moves.

Today

SWBPIPE models a spring by its rate alone; the installed load is not applied, and there is no hanger design. Status. A constant-effort hanger’s force is applied in every case, including thermal cases, and counted twice in combinations. Status.

Part II · FundamentalsII–10

J. Equipment nozzles and imposed movements

The piping is often the least of what a flexibility analysis protects. A pump casing distorted by its piping loses alignment and wears its seals; a compressor misaligns; a vessel nozzle cracks at its weld. The loads the line puts on equipment are therefore limited, and often they govern the layout more than any pipe stress.

1.Where the limits come from.

Allowable nozzle loads are not set by the piping codes. They come from the equipment standards and from the equipment vendor:

EquipmentSource of nozzle-load limits
Centrifugal pumpsAPI 610
Centrifugal and axial compressorsAPI 617
Steam turbinesNEMA SM 23; API 611 and API 612
Air-cooled heat exchangersAPI 661
Storage tanksAPI 650, Annex P
Pressure vessels and exchangersLocal-stress methods: WRC 107 (and its successor WRC 537), WRC 297, or finite elements, against the vessel code

The limits are given in the axes the standard defines, usually the nozzle axis and two directions across it, or the equipment’s own axes. Loads must be reported in those axes, as forces and moments with their signs, not as resultants in global axes.

2.Anchor or flexible nozzle.

The simplest model ends the line at the nozzle with an anchor. For pumps, turbines and compressors, whose casings are stiff, that is usually right. A vessel or tank shell is not stiff: it gives under the nozzle load, and a rigid anchor there overstates the load and puts the stress in the wrong place. Nozzle flexibility is modelled with rotational and axial springs at the shell, from WRC 297 for cylindrical vessels or from the stiffness coefficients of API 650 Annex P for tank nozzles, or from finite elements.

3.The equipment grows too.

The equipment is at temperature as well, and its nozzles move. A pump nozzle moves by the growth of the casing from the point where the pump is held, its feet or its centreline supports and the key that locates it, to the nozzle face. A vertical vessel grows upward from the base of its skirt. A tank shell bulges outward under its liquid head and rotates at the nozzle. Each is an imposed displacement at the end of the line, belonging to the operating state that causes it: a vessel at T1 moves its nozzle by one amount, at T2 by another. Like thermal expansion, an anchor movement is a secondary load and part of the displacement stress range.

4.Which cases to check.

Nozzle loads are checked in the operating case, which is usually the largest, but not only there. The sustained case gives the load the equipment carries whenever the line is cold. Where friction is present, the cases with and without it both count. Two pumps in parallel, one running and one spare, are checked with each hot and the other cold, because the cold spare still sits on a line pushed by its hot neighbour. The vendor and the standard decide which cases apply; the analysis must supply them all.

Today

SWBPIPE cannot yet apply imposed movements. Status. It reports each support’s load as a resultant magnitude only, so nozzle loads cannot be read from it. Status. It has no nozzle flexibility and no nozzle-axis reporting. Status.

Part II · FundamentalsII–11

K. Occasional loads, and how results feed a code check

Occasional loads are the last of the loads a static flexibility analysis carries. After them comes the use of all the results: each code check needs particular quantities from particular cases, and the analysis is organised to supply them.

1.Wind.

Wind on a pipe is a drag load per unit length,

w = q Cd De,   q = ½ρV2(26)

where q is the velocity pressure at the height of the pipe, Cd the drag (force) coefficient of a cylinder, which depends on the Reynolds number and the surface, and De the exposed diameter, including insulation. Wind speed rises with height, and the wind standard in use gives the profile. The load acts on the component of the wind perpendicular to each pipe; a run parallel to the wind takes little. Wind is applied from each direction that matters, and in both senses.

2.Static seismic.

The simplest seismic analysis applies a static load equal to a seismic coefficient, in g, times the weight of pipe, contents and insulation. It is applied in each horizontal direction and, where required, vertically, each in both senses, as separate cases. The directional results are then combined by the rule the governing standard gives, such as the square root of the sum of squares or a percentage rule. B31.3 does not require wind and earthquake to be taken as acting together.

3.Relief thrust.

A relief valve discharging to atmosphere pushes back on the line with a force equal to the momentum of the jet plus the excess pressure at the outlet times its area,

F = ṁVe + (pe − pa)Ae(27)

with ṁ the mass flow, and Ve, pe and Ae the velocity, pressure and area at the exit. It arrives suddenly, so as a static load it is multiplied by a dynamic load factor, up to 2. B31.1 sets out a method in a nonmandatory appendix.

Response spectrum and time-history analysis, fluid hammer, slug flow and vibration are dynamic problems, and outside the scope of this part.

Today

SWBPIPE applies wind to the whole span in the direction you give, not resolved onto the pipe, with no height profile; seismic load acts in one sense per case, on all the axes you give at once. Status.

4.From results to a code check.

Each category of code check draws on particular cases. The table sets them out, with what SWBPIPE supplies for each today.

Code categoryInputs come fromSWBPIPE today
Sustained stressThe sustained case, W + P, with springs at their hot loads; also in the operating support state where supports lift off. Moments on the section reduced by allowances, with SIFs as the code gives, plus the longitudinal pressure stress.Forces, moments and stress components per case. Longitudinal pressure stress by hand.
Displacement stress rangeThe difference of two solved states, operating less sustained or one operating state less another, formed on moment and force ranges at Ea, with in-plane and out-of-plane SIFs; the extreme pair governs.Component-by-component ranges between two result states. No SE.
Occasional stressEach occasional case, in both senses, added to the sustained result.Occasional case results, within the limits of the wind and seismic loads.
Hot and cold reactionsOperating and sustained cases, with the hot modulus and any cold spring (section D).Support load as a resultant magnitude only.
Nozzle loadsOperating and sustained cases, with and without friction, and alternate operating states, in nozzle axes.None usable today.
Support loadsThe envelope of sustained, operating, occasional and hydrotest cases at each support.Support load as a resultant magnitude only.

What SWBPIPE supplies in every row is mechanics: forces, moments and stress components. It does not compute a code stress, apply a code allowable, or judge compliance with any code. The code check is made in the program you validate in, as Part V describes.

Today

SWBPIPE has no sustained, expansion or occasional case types and no code stress. Status.

The full theory

This part teaches the discipline; Part III says what the solver does. The derivation of every equation behind it, from the element stiffness to the nonlinear iteration, is being written as the SWBPIPE theory, alongside the source in docs/theory, and will be published here chapter by chapter as Part VII.

Contents · Part I · Part III · The program: swbpipe.com · MIT licence