Physics Library
 An open source physics library
Encyclopedia | Forums | Docs | Random |  
Login
create new user
Username:
Password:
forget your password?
Main Menu
Sections

Meta

Talkback

Downloads

Information
[parent] Calculus of Variations: Beltrami Problems and Catenary-Type Examples (Example)

Calculus of Variations: Beltrami Problems and Catenary-Type Examples

CV06 derived the Beltrami identity for an autonomous first-order functional

       ∫ b
J[y] =    F (y,y′)dx,
        a
(1)

namely

|--------------|
F  − y′Fy′ = C.|
----------------
(2)

The importance of this identity is computational: it replaces the usual second-order Euler–Lagrange equation by a first-order relation. In favorable cases that relation can be solved directly for yand then integrated by a single quadrature. This companion set develops that skill through increasingly physical examples, culminating in the catenary and catenoid and ending with a preview of the brachistochrone [123].

PIC

Figure. The practical Beltrami workflow. Autonomy in the independent variable produces a first integral, which is then solved for the slope and reduced to a quadrature. boundary data determine the integration constants only after the reduction is complete.

1 How to use this set

Attempt every exercise before reading Part II. For each problem, use the same sequence:

  1. identify the integrand F and verify that Fx = 0;
  2. compute Fy carefully;
  3. form F yFy before simplifying;
  4. set the result equal to a constant;
  5. solve the first-order relation for yor dx∕dy;
  6. integrate once more, keeping track of branches and constants;
  7. apply endpoint or symmetry data only after the general stationary family has been obtained; and
  8. remember that satisfying Beltrami proves stationarity under the stated hypotheses, not automatically global minimality.

Part I: Exercises

Exercise 1: recognize when Beltrami applies

For each integrand below, state whether the Beltrami identity can be used immediately. If it can, compute the corresponding first integral.

  1. F = ∘ ------′2-
  1 + (y).
  2. F = 12(y)2 + V (y).
  3. F = x2(y)2 + y2.
  4. F = ey∘1--+-(y′)2-.
  5. F = (1 + x)(y)2.

For the cases where Beltrami does not apply, identify the explicit dependence that prevents its direct use.

Exercise 2: a master weighted-length family

Consider

       ∫ b       ---------
J [y] =    Φ (y)∘ 1 + (y′)2dx,     Φ(y) > 0.
        a
(3)

  1. Use Beltrami to show that every sufficiently smooth stationary curve satisfies
       Φ(y)
∘----------= C.
  1 + (y′)2
    (4)

  2. Solve for (y)2.
  3. Derive the quadrature
                ∫
               ----C-dy------
x −  x0 = ±    ∘ Φ-(y)2 −-C2-.
    (5)

  4. Explain why the condition Φ(y) ≥|C| appears automatically.

This exercise is the common algebraic core of several later examples.

Exercise 3: shortest path in the plane revisited

For the arc-length functional

       ∫ x1∘  ---------
L[y] =        1 + (y ′)2dx,
        x0
(6)

use Beltrami, rather than the full Euler–Lagrange equation, to show that every smooth stationary graph has constant slope and is therefore a straight line. Then impose

y (x0) = y0,    y(x1) = y1,
(7)

and obtain the explicit stationary curve.

Exercise 4: the catenary from an augmented chain functional

A uniform flexible chain in a vertical plane has gravitational potential energy proportional to

∫                ∘ ---------
   yds,     ds =   1 + (y′)2dx.
(8)

Its total length is fixed. Introducing a constant multiplier λ for that length constraint gives the augmented integrand

            ∘ ---------
F  = (y + λ)  1 + (y′)2.
(9)

For this exercise, take the augmented functional as given; CV09 will derive the variational multiplier rule systematically.

  1. Apply Beltrami and show that
       y + λ
∘-------′-2 = a,
  1 + (y )
    (10)

    where a > 0 is a constant.

  2. Solve for yand separate variables.
  3. Integrate to obtain
    |--------------------------|
|            ( x − b)      |
y (x ) = acosh  ------  − λ.|
-----------------a----------
    (11)

  4. Show that the lowest point occurs at x = b and has horizontal tangent.
  5. Explain which constants are fixed by geometry and which constant is associated with the length constraint.

PIC

Figure. A catenary-type stationary profile. The parameter b locates the lowest point, a controls the curvature scale, and an additive vertical shift is absorbed by the multiplier constant in the augmented formulation.

Exercise 5: why a shallow catenary looks parabolic

For the symmetric catenary

             (  )
y(x) = acosh  x-  − a,
              a
(12)

use the Taylor expansion of cosh z to show that, for |x|≪ a,

        x2-  -x4--
y (x) = 2a + 24a3 +  ⋅⋅⋅ .
(13)

Hence derive the leading parabolic approximation

|----------|
|       x2 |
y(x ) ≈ --.|
--------2a--
(14)

Estimate the first neglected correction term and explain why a hanging cable with small sag can appear almost parabolic even though its exact ideal-chain shape is a catenary.

PIC

Figure. A symmetric catenary and its small-sag parabolic approximation. The two agree near the lowest point because cosh z = 1 + z22 + O(z4).

Exercise 6: the catenoid from minimum surface area

A surface of revolution is formed by rotating the graph y(x) > 0 about the x-axis. Its area is

          ∫
            x1 ∘  ------′2-
A [y] = 2π     y   1 + (y ) dx.
           x0
(15)

  1. Ignore the constant factor 2π and apply Beltrami.
  2. Show that
         y
∘----------=  a.
  1 + (y′)2
    (16)

  3. Integrate the first-order equation and obtain
    |------------(------)---|
|              x − b    |
y(x ) = a cosh ------ . |
-----------------a-------
    (17)

  4. Explain why the same hyperbolic cosine appears in both the hanging chain and minimum-surface problems even though the physical functionals are different.

PIC

Figure. The generating curve of a catenoid. Rotating the catenary-shaped profile about the horizontal axis produces the classical minimal surface of revolution.

Exercise 7: autonomous mechanics and the energy integral

Let the independent variable be time t and write qt = dq∕dt. Consider the action

       ∫ t [               ]
          1 1-     2
S[q] =  t   2 m (qt) −  V(q)  dt.
        0
(18)

  1. Apply the Beltrami identity with F = L(q,qt).
  2. Show that the result can be written
    |---------------------|
|1m (q)2 + V (q) = E. |
-2----t---------------|
    (19)

  3. Solve for dt∕dq and obtain the quadrature
              ∘ ---∫
             m-  q -----dξ-----
t − t0 = ±    2     ∘E--−--V-(ξ-).
                q0
    (20)

  4. Explain why turning points satisfy V (q) = E.

Exercise 8: brachistochrone preview by Beltrami

Let y measure vertical distance downward from the starting point. Conservation of mechanical energy gives speed

     ∘ ----
v =    2gy.
(21)

The travel time along a graph y(x) is therefore

       ∫ ∘  ---------
         ---1 +-(y′)2
T[y] =      √2gy--   dx.
(22)

  1. Apply Beltrami and show that the first integral is equivalent to
    |-------′-2------|
y-(1-+-(y-)-) =-2a-
    (23)

    for some positive constant a.

  2. Introduce the parameter
    y = a(1 − cos𝜃)
    (24)

    and show that

      ′      𝜃-
y  = cot 2.
    (25)

  3. Derive
    |------------------------------------------|
-x −-x0-=-a(𝜃-−-sin-𝜃),----y-=--a(1 −-cos𝜃),-
    (26)

    which is a cycloid.

  4. Explain why this exercise is only a preview: CV17 will address the historical problem, endpoint geometry, and full interpretation in detail.

PIC

Figure. The Beltrami reduction of the brachistochrone leads naturally to a cycloidal parameterization. The curve is shown only as a preview of the full CV17 analysis.

Part II: Complete Worked Solutions

Solution 1: recognize when Beltrami applies

Beltrami requires that the integrand have no explicit dependence on the independent variable x.

(a) For

    ∘  ---------
F =    1 + (y ′)2,
(27)

we have Fx = 0. Also

            ′
  ′   ∘---y------
Fy =    1 + (y′)2.
(28)

Thus

F yFy = ∘ ------′2-
  1 + (y ) ---(y′)2---
∘  ------′2-
   1 + (y ) (29)
= ∘----1-----
  1 + (y′)2 = C. (30)

(b) For

     1
F  = --(y ′)2 + V(y ),
     2
(31)

again Fx = 0. Since Fy = y,

F −  y′F  ′ = V (y) − 1(y′)2 = C.
        y           2
(32)

(c) The integrand

F = x2(y ′)2 + y2
(33)

contains x explicitly, so Beltrami does not apply directly.

(d) For

      ∘  ---------
F = ey   1 + (y ′)2,
(34)

there is no explicit x dependence. Therefore

|----------------|
|    ey          |
|∘----------= C. |
---1-+-(y′)2-------
(35)

(e) The factor (1 + x) is explicit x dependence, so the Beltrami identity cannot be replaced by a constant. The more general du Bois–Reymond identity from CV06 must be used instead.

Solution 2: a master weighted-length family

Let

               ∘ ---------
F(y,y ′) = Φ (y )  1 + (y′)2.
(36)

Because Fx = 0, Beltrami gives

F  − y′Fy′ = C.
(37)

First compute

          -----y′----
Fy′ = Φ(y)∘  -----′-2.
             1 + (y )
(38)

Therefore

F yFy = Φ(y)∘  ---------
   1 + (y′)2 Φ(y)   (y′)2
∘--------′2-
   1 + (y ) (39)
= ---Φ-(y-)---
∘ ------′-2
  1 + (y ). (40)

Hence

|----------------|
|   Φ(y)         |
|∘-------′2-= C. |
---1-+-(y)--------
(41)

Squaring gives

                2
      ′2   Φ-(y)-
1 + (y )  =   C2  ,
(42)

so

|------------2-----2-|
|(y′)2 =  Φ(y)-−--C--.|
-------------C2------|
(43)

On an interval where a consistent branch is chosen,

dx           C
---=  ± ∘-----2-----2.
dy        Φ(y)  − C
(44)

Integrating,

|-----------∫----------------|
|              ----C-dy------|
x −  x0 = ±    ∘ Φ-(y)2 −-C2-.
------------------------------
(45)

For the square root to remain real,

Φ(y)2 − C2 ≥  0,
(46)

or

|----------|
Φ (y) ≥ |C||
------------
(47)

because Φ > 0. This restriction is not imposed separately; it is encoded in the first integral itself.

Solution 3: shortest path in the plane revisited

Here

    ∘  ------′2-
F =    1 + (y ) .
(48)

From Solution 1,

    1
∘----------= C.
  1 + (y′)2
(49)

Therefore

(y′)2 = 1--− 1,
        C2
(50)

which is constant. Choosing one continuous branch,

 ′
y =  m,
(51)

where m is constant. Hence

y = mx  + c.
(52)

Apply the two endpoint conditions:

y0 = mx0 +  c,    y1 = mx1  + c.
(53)

Subtracting,

     y1-−-y0-
m  = x1 − x0 .
(54)

Thus

|------------y-−--y----------|
|y(x) = y0 + -1----0(x − x0).|
-------------x1 −-x0----------
(55)

Beltrami has recovered the same stationary line as the full Euler–Lagrange calculation, but with one fewer differentiation step.

Solution 4: the catenary from an augmented chain functional

Take

            ∘ ---------
F  = (y + λ)  1 + (y′)2.
(56)

There is no explicit x dependence. Compute

             -----y′----
Fy ′ = (y + λ)∘1-+--(y′)2.
(57)

Then

F yFy = (y + λ)( ∘ ---------        ′2   )
    1 + (y′)2 − ∘--(y-)-----
                  1 + (y′)2 (58)
= ∘--y +-λ---
  1 + (y′)2. (59)

Set the constant equal to a positive parameter a:

|----------------|
|---y +-λ---     |
|∘1--+-(y′)2 = a.|
------------------
(60)

Rearrange:

           (y + λ)2
1 + (y′)2 =----2---,
              a
(61)

so

               2   2
(y′)2 =  (y-+-λ-)-−-a--.
             a2
(62)

Invert the derivative on a monotone branch:

dx-=  ± ∘------a-------.
dy        (y + λ)2 − a2
(63)

Let

    y + λ
u = --a---,    dy =  adu.
(64)

Then

x − b      ∫    du
------= ±    √--------=  ± arcoshu.
  a            u2 − 1
(65)

Because cosh is even, both branches combine into

y + λ        ( x − b)
------ = cosh  ------ .
  a              a
(66)

Therefore

|------------(------)------|
|              x-−-b-      |
y (x ) = acosh    a     − λ.|
----------------------------
(67)

Differentiate:

            (      )
 ′           x-−--b
y (x ) = sinh    a     .
(68)

Thus

 ′
y (b) = 0,
(69)

and since cosh z 1, the point x = b is the lowest point of the curve for this sign convention.

The constant b locates the horizontal position of the lowest point. The constant a sets the curvature scale. The multiplier λ enters as a vertical shift in this augmented form and is ultimately determined together with a and b by the endpoint geometry and the prescribed total chain length. The justification for introducing λ as a variational multiplier is the subject of CV09 [12].

Solution 5: why a shallow catenary looks parabolic

Use

             z2   z4
cosh z = 1 + ---+ ---+ ⋅⋅⋅ .
             2!   4!
(70)

With z = x∕a,

y(x) = a[     (x-)    ]
 cosh  a   − 1 (71)
= a[   2      4       ]
  x---+ -x---+ ⋅⋅ ⋅
  2a2   24a4 (72)
= x2
---
2a +  x4
---3-
24a + ⋅⋅⋅. (73)

Therefore, when |x|∕a 1,

|----------|
|       x2 |
y(x ) ≈ --.|
--------2a--
(74)

The leading neglected term is

 x4
----3.
24a
(75)

Relative to the quadratic term, its size is approximately

 4     3       2
x-∕(24a-)-=  -x--.
 x2∕(2a)     12a2
(76)

Thus if |x|∕a = 0.3, for example, the first correction is only about

0.32
 12 =  0.0075,
(77)

or less than one percent of the quadratic term. This is why shallow catenaries are visually difficult to distinguish from parabolas over a limited span.

Solution 6: the catenoid from minimum surface area

Ignoring the constant factor 2π, the integrand is

      ∘ ---------
F =  y  1 + (y′)2.
(78)

Beltrami gives

F  − y′Fy′ = C.
(79)

Since

        ----y′-----
Fy ′ = y ∘1-+-(y′)2,
(80)

we obtain

|----------------|
|-----y-----     |
|∘ ------′-2 = a,|
---1-+-(y-)-------
(81)

where a > 0. Rearranging,

         2    2
(y′)2 =  y-−-a--.
          a2
(82)

Thus

dx-     ----a-----
dy =  ± ∘ -2----2.
          y −  a
(83)

Integrating exactly as in the catenary calculation gives

|------------(------)---|
y(x ) = a cosh x-−-b- . |
|                a      |
-------------------------
(84)

The repeated hyperbolic cosine is not an accident. Both problems reduce to an integrand of the general weighted-length form

          ∘ ------′2-
F  = Φ (y)  1 + (y) ,
(85)

with a weight linear in y after a vertical shift. The physical meanings are different: the catenary comes from gravitational potential energy with a length constraint, while the catenoid comes from surface area. The algebraic structure of the Beltrami first integral is nevertheless the same [32].

Solution 7: autonomous mechanics and the energy integral

Let

               1-     2
F  = L (q, qt) = 2 m (qt) −  V (q).
(86)

Because L has no explicit time dependence, Beltrami gives

L −  q L  =  C.
      t qt
(87)

Now

Lqt = mqt.
(88)

Therefore

L qt Lqt = 1
--
2m(qt)2 V (q) m(q t)2 (89)
= [               ]
 1m (q )2 + V(q)
 2    t. (90)

Writing C = E gives

|---------------------|
|1m (q)2 + V (q) = E. |
-2----t---------------|
(91)

Solve for qt:

       ∘ --------------
         -2
qt = ±   m  [E  − V (q)].
(92)

Invert:

        ∘ ---
-dt = ±   m- ∘----1-----.
dq         2   E −  V(q)
(93)

Hence

|--------------------------------|
|          ∘ -m-∫ q      dξ      |
|t − t0 = ±    --    ∘-----------.|
--------------2--q0---E-−--V-(ξ-)-|
(94)

At a turning point, qt = 0. Therefore the energy equation requires

|----------|
-V(q)-=-E.--
(95)

This is the familiar mechanical energy integral, obtained here as a direct Beltrami first integral of an autonomous action [45].

Solution 8: brachistochrone preview by Beltrami

The integrand is

          ∘1 -+-(y′)2-
F(y,y ′) = ---√-------.
               2gy
(96)

There is no explicit x dependence. Compute

               ′
  ′   ------∘y---------
Fy  = √2gy--  1 + (y′)2.
(97)

Then

F yFy = √-1---
  2gy[∘  ---------       ′ 2   ]
    1 + (y ′)2 − ∘--(y-)----
                 1 + (y′)2 (98)
= √-----∘-1--------
  2gy   1 + (y′)2 = C. (99)

Square and absorb the positive constants into a new parameter a:

|------------------|
|y(1 + (y′)2) = 2a. |
-------------------
(100)

Now set

                        2 𝜃-
y = a(1 − cos𝜃 ) = 2a sin 2 .
(101)

The first integral gives

            2a       1
1 + (y′)2 = ---= ---2-----.
            y    sin (𝜃∕2)
(102)

Therefore

(y′)2 = cot2 𝜃,
            2
(103)

and on the descending branch

|----------|
y ′ = cot 𝜃.
---------2--
(104)

Differentiate the parameterization of y:

dy-
d𝜃 = a sin 𝜃.
(105)

Since

 ′   dy∕d𝜃-
y =  dx∕d𝜃 ,
(106)

we have

dx-
d 𝜃 = -a-sin-𝜃--
cot(𝜃∕2) (107)
= a sin 𝜃 tan 𝜃-
2 (108)
= a(1 cos 𝜃). (109)

Integrating,

|---------------------|
x − x0 = a (𝜃 − sin𝜃 ). |
-----------------------
(110)

Together with

|----------------|
|y = a(1 − cos𝜃),|
------------------
(111)

this is the parametric equation of a cycloid. The appearance of the cycloid is one of the classical achievements of the early calculus of variations. CV17 will return to the brachistochrone with its historical development and complete variational interpretation [13].

What these problems should teach

The most important lesson is not that several famous curves can be memorized. It is that one structural observation,

Fx = 0,
(112)

changes the solution strategy. Rather than expand the full Euler–Lagrange equation into a second-order ODE, first form

F  − y′F ′ = C.
        y
(113)

For weighted-length integrands of the form

          ∘ ------′2-
F  = Φ (y)  1 + (y) ,
(114)

this immediately becomes

   Φ(y)
∘----------= C,
  1 + (y′)2
(115)

which often exposes the geometry of the problem before any difficult integration begins.

Common mistakes

  • Using Beltrami when F contains x explicitly. The constant first integral requires Fx = 0.
  • Computing Fy incorrectly. In Φ(y)∘1--+-(y′)2, the factor Φ(y) is held constant when taking the partial derivative with respect to y.
  • Dropping the square-root domain condition. Solving for ycan introduce a requirement such as Φ(y)2 C2.
  • Forgetting the branch sign. A first-order relation usually gives y= ±f(y). A complete smooth curve may switch monotone branches at a turning point.
  • Treating the multiplier λ in the catenary as arbitrary decoration. It enforces the fixed-length constraint; CV09 provides the theorem justifying it.
  • Calling every hyperbolic-cosine graph a hanging chain. The same analytic profile also generates a catenoid, but the underlying variational functional is different.
  • Assuming a first integral proves a minimum. Beltrami supplies a necessary stationarity relation. Classification is a separate issue.

Summary

For an autonomous functional

       ∫
                ′
J [y ] =   F (y,y) dx,
(116)

stationarity implies the Beltrami identity

|-----′--------|
F--−-y-Fy′ =-C.-
(117)

For the common weighted-length family

          ∘ ---------
F  = Φ (y)  1 + (y′)2,
(118)

this reduces to

|----------------|
|---Φ(y)----     |
|∘1--+-(y′)2-= C. |
------------------
(119)

The same algebraic structure produces a straight line for ordinary planar arc length, a hyperbolic cosine for the catenary and catenoid, an energy integral in autonomous mechanics, and the first-order cycloidal relation in the brachistochrone. CV06E2 next emphasizes the complementary special case of cyclic dependent variables and conserved conjugate momenta.

References

[1]   I. M. Gelfand and S. V. Fomin, Calculus of Variations, Dover Publications, 2000.

[2]   Bruce van Brunt, The Calculus of Variations, Springer, 2004.

[3]   Robert Weinstock, Calculus of Variations with Applications to Physics and Engineering, Dover Publications, 1974.

[4]   Cornelius Lanczos, The Variational Principles of Mechanics, 4th ed., Dover Publications, 1986.

[5]   Herbert Goldstein, Charles Poole, and John Safko, Classical Mechanics, 3rd ed., Addison Wesley, 2002.


"Calculus of Variations: Beltrami Problems and Catenary-Type Examples" is owned by bloftin.
(view preamble)
View style:
Other names:  CV06E1
Keywords:  calculus of variations, Beltrami identity, first integral, catenary, catenoid, autonomous integrand, weighted arc length, quadrature, brachistochrone, energy integral, Euler-Lagrange equation, worked exercises

This object's parent.

Cross-references: mechanics, theorem, domain, position, square, speed, parameter, energy, graph, algebraic, boundary, catenary, relation, identity, CV06

This is version 1 of Calculus of Variations: Beltrami Problems and Catenary-Type Examples, born on 2026-09-13.
Object id is 1203, canonical name is CalculusOfVariationsBeltramiProblemsAndCatenaryTypeExamples.
Accessed 10 times total.

Classification:
Physics Classification02.30.Xx (Calculus of variations)
 02.30.Hq (Ordinary differential equations)
 45.20.Jj (Lagrangian and Hamiltonian mechanics)
Pending Errata and Addenda
None.
Discussion
Style: Expand: Order:

No messages.

Interact
rate | post | correct | update request | add example | add (any)