Featured Image 10

The Role of Catenary Equations in Sag and Tension Calculations

The catenary equation is the backbone of transmission line design — it's the actual shape a conductor takes under its own weight and tension, and it drives every clearance check, sag template, and tower spotting decision that follows. Every shortcut used later, from the parabola to the ruling span, is a simplification of this curve — and knowing where it comes from is what tells you when those shortcuts are safe to use. This post derives it fully, from force balance to final formula.

This post follows the derivation in Kiessling et al., Overhead Power Lines, §14.1–14.2.Notation is converted to this site's convention throughout: S = span, D = sag, L = conductor length, w = weight per unit length, H = horizontal tension.

The catenary equation appears throughout this site's sag-tension calculators and worked examplesThis post derives it directly, starting from a force balance on an infinitesimal element of conductor and working through every intermediate step, then continues to the sag, maximum sag, and conductor length formulas that follow from it.

Notation used in this post
SymbolMeaning
HHorizontal component of conductor tension (constant along the span)
VVertical component of conductor tension (varies along the span)
wConductor weight per unit length
x, yHorizontal and vertical coordinates
dLInfinitesimal length of conductor
SSpan length (horizontal distance between supports)
DSag (vertical drop below the line joining the two supports)
LTotal conductor length between supports
x₀, y₀Constants of integration (coordinates of the origin, fixed in Part 3)
Part 1: Force balance on a conductor element

Let's isolate an infinitesimal segment of conductor, length dL, hanging under its own weight. Three forces act on it: weight w·dL acting downward, tension (H, V) at one end, and tension (H+dH, V+dV) at the other end.

Post Media32
Infinitesimal Conductor Length.
Vertical equilibrium
\[ w\,dL = dV \](1)

Arc length in terms of the coordinates:

\[ dL = \sqrt{dx^2+dy^2} = dx\sqrt{1+(dy/dx)^2} \](2)

Substituting (2) into (1) and dividing by dx:

\[ \frac{dV}{dx} = w\sqrt{1+(dy/dx)^2} \](3)
Horizontal equilibrium

No horizontal load acts on the element:

\[ dH = 0 \](4)

Equation (4) integrates to H = const. The horizontal tension is the same at every point along the span.

Tangent condition

Tension in a flexible cable is always directed along the tangent to the curve, so the ratio of vertical to horizontal components equals the local slope:

\[ V = H\frac{dy}{dx} \](5)

Differentiate (5) with respect to x. Since H is constant, dV/dx = H·d²y/dx². Substitute into (3):

\[ \frac{d^2y}{dx^2} = \frac{w}{H}\sqrt{1+\left(\frac{dy}{dx}\right)^2} \](6)

This is the governing differential equation of the sag curve.

Part 2: Solving the differential equation

Equation (6) is separable, but it takes a specific trick to get there. Start by dividing both sides by √(1+(y')²):

\[ \frac{y''}{\sqrt{1+(y')^2}} = \frac{w}{H} \](6a)

Multiply both sides by y'·dx:

\[ \frac{y''\,y'}{\sqrt{1+(y')^2}}\,dx = \frac{w}{H}\,y'\,dx \](6b)

The left side is the derivative of √(1+(y')²), verified directly by the chain rule:

\[ \frac{d}{dx}\sqrt{1+(y')^2} = \frac{y'\,y''}{\sqrt{1+(y')^2}} \](6c)

Since y'·dx = dy, equation (6b) is equivalent to:

\[ d\left[\sqrt{1+(y')^2}\right] = \frac{w}{H}\,dy \](6d)

Both sides are exact differentials. Integrating:

\[ \sqrt{1+(y')^2} = \frac{w}{H}\,y + C \](6e)

C is undetermined at this stage; no boundary condition has been applied. Write C = -(w/H)y_0, defining a new constant y_0 := -CH/w. This is a valid algebraic substitution since w/H ≠ 0. Equation (6e) becomes:

\[ \sqrt{1+(y')^2} = \frac{w}{H}(y-y_0) \](7)
Second integration

Equation (7) is separable. Substituting p = (w/H)(y-y_0), dy = (H/w)dp:

\[ \frac{H}{w}\cdot\frac{dp}{\sqrt{p^2-1}} = dx \](8)

Integrating (the left side is the definition of cosh⁻¹), and introducing a second constant x₀:

\[ \frac{H}{w}\cosh^{-1}(p) = x - x_0 \](9)

Solving for y:

\[ y = y_0 + \frac{H}{w}\cosh\left[\frac{w}{H}(x-x_0)\right] \](10)

Two constants, x₀ and y₀, matching the order of the differential equation. This is the general solution for any coordinate origin.

Part 3: Locating the vertex

Let's find the vertex. It's the point of zero slope:

\[ \left.\frac{dy}{dx}\right|_{\text{vertex}} = 0 \](11)

Differentiating (10) gives y' = sinh[(w/H)(x-x₀)], which is zero exactly at x = x₀. The vertex is located at x = x₀. Substituting into (10):

\[ y_{\text{vertex}} = y_0 + \frac{H}{w} \](12)

Equation (12) holds for any choice of origin. Setting x₀ = 0 and y₀ = 0 places the origin directly below the vertex. Equation (10) reduces to:

\[ y = \frac{H}{w}\cosh\left(\frac{w\,x}{H}\right) \](13)

By (12), with y₀ = 0, the vertex sits at height H/w above this origin, not at y = 0. The quantity c = H/w is the catenary parameter; the standard form y = c·cosh(x/c) has its vertex at height c, consistent with (12).

Finding the vertex
Location of the catenary vertex.
Part 4: Small-slope (parabolic) approximation

Now for the shortcut most engineers actually reach for. For typical span lengths and sags in overhead line design, the slope dy/dx stays small enough that (dy/dx)² ≪ 1. Under this condition, (6) reduces to:

\[ \frac{d^2y}{dx^2} \approx \frac{w}{H} \](14)

Integrating once, and applying the vertex condition dy/dx = 0 at x = 0:

\[ \frac{dy}{dx} = \frac{w}{H}x \](15)

Integrating again, with the parabola's own vertex placed at y = 0:

\[ y = \frac{w}{2H}x^2 \](16)

Equation (16) is the parabolic approximation to the sag curve, measured as height above the vertex. It was derived independently of the catenary solution. The connection between the two is established by expanding cosh in equation (13) as a Taylor series around x = 0:

\[ \frac{H}{w}\cosh\left(\frac{wx}{H}\right) = \frac{H}{w} + \frac{w\,x^2}{2H} + \frac{w^3x^4}{24H^3} + \dots \](17)

The first term, H/w, is the vertex height from (12). Subtracting it — that is, measuring height above the vertex, the same reference used in (16) — leaves:

\[ y_{\text{catenary}} - \frac{H}{w} = \frac{w\,x^2}{2H} + \frac{w^3x^4}{24H^3} + \dots \](18)

The leading term of (18) is identical to (16). The parabola is the small-slope limit of the catenary; the terms omitted in (14) are exactly the higher-order terms of this series. Their relative size, and therefore the accuracy of the approximation, depends on wS/H — larger for longer spans or slacker conductors, negligible for typical transmission spans.

Part 5: Sag relative to the line joining the supports

Here's a subtlety worth catching early: equation (13) gives height relative to the vertex, not sag. Sag is the vertical drop below the straight line connecting the two support attachment points, A and B.

Sag at level span
The maximum sag is located exactly half the span for level supports.

Let support A sit at horizontal position x_A (measured from the vertex), support B at x_B = x_A + S, with height difference h between them. The straight line from A to B has height (h/S)(x-x_A) above A at position x. Sag D is the gap between that line and the curve:

\[ D = \frac{h}{S}(x-x_A) + \frac{H}{w}\left[\cosh\left(\frac{w\,x_A}{H}\right)-\cosh\left(\frac{w\,x}{H}\right)\right] \](19)

For a level span (h = 0), the vertex sits at midspan, x_A = -S/2. Substituting into (19):

\[ D = \frac{H}{w}\left[\cosh\left(\frac{wS}{2H}\right)-\cosh\left(\frac{w\,x}{H}\right)\right] \](20)

Maximum sag occurs at the vertex (x = 0):

\[ D_{max} = \frac{H}{w}\left[\cosh\left(\frac{wS}{2H}\right)-1\right] \](21)

Expanding (21) as a series produces the same structure as (17)–(18): the leading term is wS²/(8H), the parabolic sag formula, with higher-order correction terms from the exact catenary.

Part 6: Conductor length

One more quantity worth having in closed form: the actual length of conductor between the two supports.

\[ L = \int_{x_A}^{x_B}\sqrt{1+(y')^2}\,dx \](22)

From (7) with the vertex-origin convention, y' = sinh(wx/H), so 1+(y')² = cosh²(wx/H). The square root is exact:

\[ L = \int_{x_A}^{x_B}\cosh\left(\frac{w\,x}{H}\right)dx = \frac{H}{w}\left[\sinh\left(\frac{w\,x_B}{H}\right)-\sinh\left(\frac{w\,x_A}{H}\right)\right] \](23)

For a level span:

\[ L = \frac{2H}{w}\sinh\left(\frac{wS}{2H}\right) \](24)
Part 7: Inclined spans — exact vertex location

Let's go back to Part 5's vertex position, given only implicitly through equation (16), and solve it in closed form. Using the identity cosh(p) − cosh(q) = 2·sinh((p+q)/2)·sinh((p−q)/2) applied to equation (16):

\[ x_A = \frac{H}{w}\sinh^{-1}\left[\frac{h}{2\frac{H}{w}\sinh\left(\frac{wS}{2H}\right)}\right] - \frac{S}{2} \](25)

Equation (25) is exact — not an approximation. Verified numerically against equation (16) across multiple test cases, including negative h (support B lower than A).

For extreme combinations of a short, steep span and low tension, equation (25) can place the vertex outside the physical span (x_A < -S/2 or beyond support B). When this happens, the conductor's lowest point on the span is the lower support itself, not a true vertex — the curve is monotonically rising or falling across the whole span. Equation (25) is still mathematically correct in this case; it just describes a point that lies outside the interval actually built.

This is not a corner case worth skipping. It comes up directly in dead-end or angle structures on steep terrain, and in short jumper spans at substations — exactly the spans where a sag template built assuming a mid-span vertex would place the clearance check at the wrong point entirely. Before running a clearance check on a short, steep span, confirm equation (25) actually falls within [-S/2, S/2] first.

Sag at unlevel span
The maximum sag is not located in the vertex for unequal level supports.
Location of maximum sag

For a level span, maximum sag occurs at the vertex. For an inclined span, it does not — sag is measured relative to the straight line joining the supports, and that line itself is tilted. Differentiating the sag equation (19) with respect to x and setting it to zero:

\[ \frac{dD}{dx} = \frac{h}{S} - \sinh\left(\frac{wx}{H}\right) = 0 \](26)

Solved directly, measuring x from the vertex (not from support A):

\[ x_c = \frac{H}{w}\sinh^{-1}\left(\frac{h}{S}\right) \](27)

Since sinh⁻¹ is an odd, strictly increasing function, x_c has the same sign as h. Maximum sag sits on whichever side of the vertex is closer to the higher support. Verified numerically for positive, negative, and zero h: equation (27) matches the numerically located maximum in every case, to within numerical precision.

This direction result is established in Hatibovic (2014) using support-origin coordinates, requiring a multi-step inequality argument. In the vertex-origin coordinates used throughout this post, the same conclusion follows directly from the sign and monotonicity of sinh⁻¹ in equation (27) — one property of one function, rather than a chain of inequalities.

For hillside spans, this means the governing clearance point is not automatically at midspan or at the vertex — it shifts toward the higher structure. On a sloped profile, checking clearance only at the vertex or only at midspan can miss the actual worst-case point. This is the same effect worked numerically in Problem 4 of the solved-problems post, using the parabolic approximation; equation (27) gives the exact location for the same physical situation.

Part 8: Tension at any point

One last useful quantity: total conductor tension, T = √(H² + V²). Using V = H·sinh(wx/H) from equation (5) and its integrated form:

\[ T = H\sqrt{1+\sinh^2\left(\frac{wx}{H}\right)} = H\cosh\left(\frac{wx}{H}\right) \](28)

Comparing to equation (13), H·cosh(wx/H) = w·y, so:

\[ T = w\,y \](29)

Total tension at any point equals the conductor weight per unit length times the height above the reference line — the same reference line established in Part 3, sitting H/w below the vertex. Verified numerically against the direct calculation of (28) at several positions along the curve.

This gives tension at the support attachment points directly, without separately computing V first — relevant for checking conductor tensile stress against its rated strength at the point where tension is actually highest, which is always at the higher-elevation support, not at the vertex.

If you go looking at how other references handle this, you'll find Hatibovic (2014) sets his origin at the lower support instead of the vertex, which suits direct clearance calculations but means carrying the vertex coordinates as separate quantities throughout his derivation. Deakin (2019) uses the same vertex-origin coordinates as this post and arrives at equation (29) independently, written as T = λgy. Haukaas's UBC lecture notes take yet another route, working with the tangent angle directly instead of dy/dx, and land in the same place.
Summary of equations
QuantityFormulaEq.
Governing differential equation\( \dfrac{d^2y}{dx^2} = \dfrac{w}{H}\sqrt{1+(dy/dx)^2} \)(6)
General solution (any origin)\( y = y_0 + \dfrac{H}{w}\cosh\left[\dfrac{w}{H}(x-x_0)\right] \)(10)
Vertex height above the origin\( y_{\text{vertex}} = y_0 + \dfrac{H}{w} \)(12)
Catenary equation, vertex-origin form\( y = \dfrac{H}{w}\cosh\left(\dfrac{wx}{H}\right) \)(13)
Parabolic approximation (small slope)\( y \approx \dfrac{w}{2H}x^2 \)(16)
Sag at any point, unequal-elevation supports\( D = \dfrac{h}{S}(x-x_A) + \dfrac{H}{w}\left[\cosh\left(\dfrac{wx_A}{H}\right)-\cosh\left(\dfrac{wx}{H}\right)\right] \)(19)
Sag at any point, level span\( D = \dfrac{H}{w}\left[\cosh\left(\dfrac{wS}{2H}\right)-\cosh\left(\dfrac{wx}{H}\right)\right] \)(20)
Maximum sag, level span\( D_{max} = \dfrac{H}{w}\left[\cosh\left(\dfrac{wS}{2H}\right)-1\right] \)(21)
Conductor length, general\( L = \dfrac{H}{w}\left[\sinh\left(\dfrac{wx_B}{H}\right)-\sinh\left(\dfrac{wx_A}{H}\right)\right] \)(23)
Conductor length, level span\( L = \dfrac{2H}{w}\sinh\left(\dfrac{wS}{2H}\right) \)(24)
Vertex location, inclined span (exact)\( x_A = \dfrac{H}{w}\sinh^{-1}\left[\dfrac{h}{2\frac{H}{w}\sinh(wS/2H)}\right] - \dfrac{S}{2} \)(25)
Location of maximum sag, from vertex\( x_c = \dfrac{H}{w}\sinh^{-1}\left(\dfrac{h}{S}\right) \)(27)
Tension at any point\( T = w\,y \)(29)
This post derives the exact catenary relations. The parabolic approximation introduced in Part 4 is the small-slope limit of these results — the next post derives the parabola's own sag, maximum sag, and length formulas directly (not as a series truncation of the catenary), and quantifies the error against the results above.
Summary

Equations (1)–(29) derive from one physical fact — constant horizontal tension along the span, equation (4) — combined with calculus. The parabolic approximation in Part 4 is not a separate model; it is the leading term of the exact catenary, verified directly by series expansion in equations (17)–(18). Parts 7 and 8 extend the level-span results to inclined spans and to tension at any point, verified numerically against independent literature. The only place a choice is made rather than derived is the placement of the coordinate origin in Part 3, and that choice is stated explicitly rather than assumed.

References:
Kiessling, F.; Nefzger, P.; Kaintzyk, U.; Nefzger, F.: Overhead Power Lines: Planning, Design, Construction. Berlin, Springer-Verlag, Chapter 14.
Hatibovic, A.: "Derivation of Equations for Conductor and Sag Curves of an Overhead Line Based on a Given Catenary Constant." Periodica Polytechnica Electrical Engineering and Computer Science, 58/1 (2014), 23–27.
Deakin, R.: "Catenary Curve." Dunsborough, WA, Australia, 2019.
Haukaas, T.: "Cables and the Catenary Shape." University of British Columbia lecture notes.

Similar Posts

Leave a Reply

Your email address will not be published. Required fields are marked *