Mathematical Problems in Engineering

Volume 2014, Article ID 231726, 11 pages

http://dx.doi.org/10.1155/2014/231726

## The Static WKB Solution to Catenary Problems with Large Sag and Bending Stiffness

Department of Construction Engineering, National Taiwan University of Science and Technology, Taipei City 10607, Taiwan

Received 30 June 2014; Accepted 11 September 2014; Published 28 September 2014

Academic Editor: Bruno Briseghella

Copyright © 2014 Yuhung Hsu and Chanping Pan. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Large sag with a bending stiffness catenary is a subject that draws attention in the realm of fatigue analysis, estimation of suspension cable sag for bridge cable hoisting, and ocean engineering of the employment of mooring systems. However, the bending stiffness is the cause of boundary layers at the anchorage of cables, thereby finding a solution of the differential equation can be extremely difficult. Previous studies have tackled this problem with the perturbation method; yet, due to the complexity of the matching process and solution finding, the method might not be an ideal solution for engineering applications. Moreover, the finite difference method and the finite element method in numerical analysis can often be ineffective because of inappropriate parameter configuration and the drastic variation of functions in the boundary layers. Therefore, this study proposed a novel bending moment expression of a large sag catenary. The expression was derived from the sag identified using bending moment equations, and a solution was found by applying the WKB method (Wentzel-Kramers-Brillouin method) to overcome the complex problem of boundary layers. Consequently, a simple solution of various mechanical properties in a cable with bending stiffness and large sag could be obtained.

#### 1. Introduction

Slim tension members have been comprehensibly applied to civil and ocean engineering to create cables of cable-stayed and suspension bridges that satisfy traffic demands and mooring systems that satisfy the demands of mining deep-sea petrol and natural gas. Such structural problems could be understood as general cable problems where the mechanical behavior of cables varies significantly with the magnitude of tension. In the study of cables of bridges, researchers often regard the ratio of midpoint sag to span as an appropriate standard for simplifying the parabola theory. When the ratio of midpoint sag to span is 1/8, the horizontal force happens to be the total cable weight. When the ratio of midpoint sag to span is smaller than 1/8 (horizontal force greater than the total cable weight), the resulting parabola is appropriate and correct. However, when the ratio of midpoint sag to span is greater than 1/8 (horizontal force smaller than the total cable weight), a greater error can be found in the parabola [1]. The relatively high tension applied on the cables of a cable-stayed bridge can result in a deflected shape that can be simplified as a parabola because of its relatively small sag and their bending stiffness have only minor effects and are negligible; however, deflected shapes of cables with relatively low tension (main cables of a suspension bridge or marine risers in ocean engineering) can be regarded as a catenary because the large sag increases the influence of the bending stiffness and cannot be neglected. In addition, the bending moment and stress of the anchor points are important parameters in fatigue analysis. Damages to cables caused by fatigue result in significant losses; thus, the understanding of mechanical behavior in cables with bending stiffness demands immediate attention.

The problem, however, of incorporating the bending stiffness effect in cables results from the complex differential equations. Bending stiffness is the cause of boundary layers in anchorages and rapid variations of bending moments occur near the differential regions of the cable anchorage. Burgess [2] utilized the finite difference method to analyze cables with bending stiffness and low tension and indicated the importance of parameter configuration to avoid the loss of effectiveness when employing general finite difference or finite element methods. Reasonable precision should be obtained by adopting a smaller step size and finer grid. Moreover, an approximation algorithm (e.g., perturbation method) should be adopted to find the analytical expression, which cannot be found by an analytical approach because of the complexity of differential equations.

In dealing with the boundary layer problems, the perturbation method consists of two main techniques [3–5]: matched asymptotic expansions and multiple scales method. Most scholars adopted the matched asymptotic expansions to analyze boundary layers in cables with bending stiffness, as M. S. Triantafyllou and G. S. Triantafyllou [6] analyzed the hanging string, Wolfe [7] analyzed the inextensible cable, Irvine [8] analyzed the stay-cables, Stump and Fraser [9] analyzed the high-speed transport of thin-sheet materials, and Stump and van der Heijden [10] analyzed bending and twisting rods. In later research, Denoël and Detournay [11] employed multiple scales method. Furthermore, Denoël and Canor [12] embraced a new perturbation method, patching asymptotics, for analysis. Nonetheless, methods of perturbation are not convenient choices for engineering applications. All perturbation methods require the matching of solutions of cable segments outside the boundary layers and beam segments inside the boundary layers and the matching procedures as well as the solutions are extraordinarily complicated.

To correct the weaknesses of a conventional perturbation method, this study proposed a novel catenary bending moment equation. The core idea originated from the comparison between two catenary models, which consisted of anchorages with the same horizontal distance and vertical elevation. The self-weight and horizontal force between the two models were also identical. The only difference is that one model contains bending stiffness but the other does not. The bending moment solution of the model with a bending stiffness catenary can be expressed by that of the model without a bending stiffness catenary. By winding around the difficulty of a conventional perturbation method, which requires a fourth order differential equation for finding the solution, the proposed equation could directly identify cable sag with the help of bending moment equations.

To satisfy the demands of different cable forces and to verify the correctness of the proposed bending moment equation, two respective second order differential equations of parabola and catenary incorporating bending stiffness were established and discussed. The parabola was expressed by a small sag linear moment-curvature relation, and the solution of the parabola was used to verify the catenary result. The catenary was expressed by a large sag nonlinear moment-curvature relation; moreover, the WKB method (Wentzel-Kramers-Brillouin method is named after physicists Wentzel, Kramers, and Brillouin, who all developed it in 1926.) in the perturbation method [3–5] was adopted to solve the sag differential equations. The solution could resolve the complexity of boundary layers and does not require the matching of cable segments and beam segments as does a general perturbation method. With only one equation, the proposed method could describe the drastic variations in functions inside and outside the boundary layer and find the sag and bending moment effortlessly. Furthermore, this study established the influence of bending stiffness on cables’ static behavior and the standard where the parabola could replace the catenary.

#### 2. Static Theories of Catenaries and Parabolas without Bending Stiffness

##### 2.1. The Catenary Solution

Figure 1 illustrates a catenary problem that considers only tensile forces without the bending moment. Self-weight effects are given and set, and the catenary is located in the first quadrant. The equation of static equilibrium parallel to the -axis is is the tension, is the horizontal tension, and is a constant. The equation of static equilibrium perpendicular to the -axis is is the per-unity-length mass of the cable, is acceleration of gravity, and is the length of the catenary curve. is the horizontal distance between the anchorages. is uniformly distributed along the curved cable, and . Dimensionless parameters are introduced as follows: , , and , where implies the effects of cable self-weight on the equation of static equilibrium. Equation (2) is modified and becomes is the dimensionless sag of the catenary, where the prime denote differentiation with respect to . It can be substituted into the boundary conditions of the hinges on the two ends ( and ), and the solution of (3) is the solution of the catenary. Consider

Consider . The two hinges are at the same elevation level (, ); (4) can be simplified as

##### 2.2. The Parabola Solution

As shown in Figure 2, given that the self-weight distributes uniformly along the chordwise direction, is a constant, and is the chordwise inclination of the suspension cable. The suspension cable length indicates the shortest distance between the two anchorages at the two ends. Equation (3) is modified as

The boundary condition is set to the hinges at the two ends and the solution can be expressed by a parabola. Consider

The two hinges are at the same elevation level (), and (7) is simplified as

Equation (6) must accept the small sag cable assumption. Most previous studies regarded a ratio of midpoint sag to span below 1/8 as small sag cable (, ). On the contrary, a ratio of midpoint sag to span greater than 1/8 is considered as large sag cable (, ) [1].

#### 3. Conventional Solution of the Equation of Static Equilibrium of Catenaries and Parabolas with Bending Stiffness

##### 3.1. Differential Equation of Catenaries with Bending Stiffness

As shown in Figure 3, the bending moment and shear are incorporated into the catenary. The equation of static equilibrium parallel to the -axis is

The horizontal tensile force is a constant. The equation of static equilibrium perpendicular to the -axis is

The equilibrium of cable segments yields

Compilation yields

The moment-curvature relation of a large-sag cable is is flexural rigidity, the elastic modulus, and is the second moment of area. In (13), bending moments applied with compression on the upper part of the cable is regarded as positive. By substituting (13) into (12),

##### 3.2. Solution of Parabolas with Bending Stiffness

Due to the difficulty in finding a solution to (14), previous methods could not solve catenaries consisting of bending stiffness and large sag. In the textbook [1], the most common solution is to simplify (12); it is a classical small sag simplified solution using linear curvature. Assuming a tense cable caused by an enormous tensile force, neglect the nonlinear terms in (13) and simplify it as . Assuming (12) being a parabola, simplify it as

Dimensionless parameter is introduced to show the influence of bending stiffness. Equation (15) is modified as

Boundary conditions of hinges (, , and ) are substituted into the solution to (16), yielding a sag equation that includes hinges and excludes bending moments at the two ends. Consider

The two hinges are at the same elevation level (), and (17) is simplified as

Dimensionless bending moment equation of hinges on a horizontal cable is

Bending moments applied with compression on the upper part of the cable are regarded as positive. The substitution of the fixed-end boundary condition shows that the slope of the hinge pivot is not equal to zero but to the slope of the cable (, , , ). The obtained dimensionless fixed-end bending moment is and are bending moments at the anchorage, having positive values when rotating counterclockwise. The sag equation is

For a horizontal parabola (), (21) can be simplified as

The dimensionless bending moment equation for the fixed-end on a horizontal cable is

Equations (17) and (21) are solutions based on parabolas. This is different from the large sag cable approximation in [7–11].

#### 4. Novel Catenary Model of Static Equilibrium with Bending Stiffness

##### 4.1. Novel Catenary Bending Moment Equation

To overcome the difficulty of finding a solution using (14), this study adopted a novel concept to modify (14). Figure 4 is a comparison model between number 1 and number 2 catenaries. Based on (5), number 1 catenary was set to include the tensile force but not the bending moment. Consider is the per-unity-length mass of number 1 catenary, is the length of number 1 catenary curve, and is the sag of number 1 catenary.

Based on (12), number 2 catenary was introduced to include the shear and bending moment. Consider is the per-unity-length mass of number 2 catenary, is the length of number 2 catenary curve, and is the sag of number 2 catenary.

The mechanical relation existing between number 1 and number 2 catenaries was independent and indirect. The two catenaries had two definitely different shapes, , and had different masses along the chord, . However, they shared the following in common: they had the same horizontal distance and vertical elevation between the two anchorages, the same identical total cable self-weight, and the same horizontal tensile force, .

Mathematically, (24) can be subtracted by (25) as follows:

By integrating (26), and (27) is modified as

Another integration could obtain the bending moment function, . Consider

Equation (29) is a new expression of the cable bending moment equation. The first characteristic allows it to avoid the difficulty of (14), which requires a fourth order differential equation to find the solution; instead, (29) directly finds the sag from the bending moment equation. The second characteristic is that a relation is not needed between sag and of the number 1 and number 2 catenaries. Simply put, (29) is a mathematical equation that expresses the bending moment of number 2 catenary from the sag of number 1 catenary. The two catenaries happen to have the same self-weight and horizontal tensile force. In other words, any form of can be selected to calculate the with bending stiffness.

Constants of integration and are obtained based on the boundary conditions. The two catenaries shared the same anchorage elevation, and . and are the bending moments at the anchorage that show positive value when rotating counterclockwise. Bending moments, when applied with compression on the upper part of the cable, are regarded as positive. The substitution of and results in and . Equation (29) is modified as

The dimensionless bending moment equation is

Equation (30) can also be easily derived by examining the static equilibrium of the free-body diagram of number 1 and number 2 catenaries.

##### 4.2. Verification of the Correctness of the Novel Bending Moment Equation Using the Parabola

This study used the parabola to verify (31). Let (7) be , . By substituting and into (31),

The correctness of (29) can be verified with (which is (21)), a solution obtained by substituting the fixed-end boundary condition.

##### 4.3. Novel Catenary Differential Equation

The moment-curvature relation of the large sag catenary is , and nonlinear terms are not neglected. Parameter is set to represents difference of sag between number 1 and number 2 catenaries caused by bending moments. and . By substituting the results into (31) to obtain the differential equation of large sag catenary,

Given that (4) is , and . In a cable-stayed bridge, 95% of the cables have [13]. , and can be replaced by . Subsequently, (34) becomes the differential equation with a bending moment sag of . Consider

The large horizontal tensile force and small bending stiffness make much smaller as compared to other parameters. The in (35) is multiplied by small parameter to form a paradigmatic boundary layer problem that involves the multiplication of a highest order derivative with small parameter . This implies that an enormous function variation exists in the differential region of the catenary anchorage. Leaving the differential region, function variations quickly come to a mild plane and maintain at a steady state. The division that signifies the rapid and drastic function variation is known as the boundary layer and the differential region at the catenary anchorage is known as the thickness of the boundary layer. Since (35) did not have an analytical solution, this study adopted the WKB method in the perturbation method to find the approximation solution.

#### 5. WKB Catenary Solution with Bending Stiffness

First, the homogeneous solution of (35) was found as follows:

The WKB approximation method developed by Wentzel, Kramers, and Brillouin [2–4] was adopted to find the solution to (36).

In the classical Sturm-Liouville equation [5], . When , the first order approximation is . When , the first order approximation is .

Equation (36) is , and the WKB approximation is

The result of in (37) is is the elliptic integral of the first kind, , and is the imaginary unit. Consider is too small and therefore the nonhomogeneous solution could be obtained by neglecting in (36). Consider

The final result is . Define

The total sag is . Consider

For a horizontal catenary (, ), (41) can be simplified as . The total sag is

A solution could be found by substituting (42) and (44) into the boundary conditions. The advantage of the horizontal catenaries in (44) is their arbitrary abilities to be substituted into any desired boundary conditions. When the catenaries are hinges and do not possess bending moments, . Considering that and , let and simplify the solution to yield

The substitution of (45) into (43) generates the bending moment sag equation for horizontal hinged catenaries as follows:

The substitution of (45) into (44) generates the total sag equation for horizontal hinged catenaries as follows:

Boundary conditions of fixed-ends (, , , ) are substituted into the solution to (47), and the obtained constants of integration are , and the bending moment, , of the left fixed-end is

Substitute (48) and (49) into (43) to obtain the bending moment sag equation for horizontal fixed-end catenaries as follows:

Substitute (48) and (49) into (44) to obtain the total sag equation for horizontal fixed-end catenaries as follows:

Substitute (43) into (31) to obtain the bending moment equation of horizontal fixed-end catenaries as follows: is illustrated in (45) and (48). Without the necessity of a general perturbation method [10, 11] that requires the matching of cable and beam segments, the proposed equations (44) and (52) can describe the drastic function variations within and without the boundary layer with a single equation.

#### 6. Results and Discussion

Figure 5 shows a comparison between the total sag of parabolas and catenaries. The figure indicates that under any and , the WKB solution of catenaries derived from (51) was between (5) and (22). At the same (1.0), the three solutions demonstrated greater differences as reduced. As and , the total sag retrieved by (51) of the cable span midpoint was −0.105. This corresponds to [12]. At the same (), the three solutions demonstrated greater differences as increased. When , (22) shows a good approximation of (51).

Fundamental mechanics show that the difference between (5) and (8) increased with . Therefore, when discussing the influence of on sag, the sag of (5) and (8) should be divided by (51) and (22) to compare with sag caused by the bending moment. The bending moment sag, , corresponds to the solution of (22) that contains a bending stiffness parabola. Consider

First the influence of bending moment sag, , on catenaries is evaluated.

Define the ratio of span midpoint bending moment sag to hinged catenaries corresponding to division of (5) and (46). Equation (5) is the sag of the catenary (), and (46) is the bending moment sag equation for horizontal hinged catenaries (). Consider is the ratio of span midpoint bending moment sag to hinged catenaries.

Under four different ’s (0.5, 1, 2, and 3), Figure 6 views the ratio of span midpoint bending moment sag to catenaries supported by a hinge. The ratio was found to drastically fluctuate between . At all the values, occupied a higher ratio when was smaller, a ratio that amounted to 7.86%; occupied a smaller ratio when was larger. When , the ratio was only 1.7%, and when the ratio was only 0.3%. Bending stiffness had effective influences when in the hinged catenary.

Define the ratio of span midpoint bending moment sag to fixed-end catenaries corresponding to division of (5) and (50). Equation (50) is the bending moment sag equation for horizontal fixed-end catenaries (). Consider is the ratio of span midpoint bending moment sag to fixed-end catenaries. is illustrated in (48). is illustrated in (49).

Figure 7 demonstrates the fixed-end catenary sag ratio occupied by span midpoint bending moment sag . The figure showed a variation trend identical to Figure 6 but with a greatly enhanced value. When , had a 35.9% ratio, had a 6.9% ratio, and had a 3.4% ratio. Obviously whether the boundary condition is hinge or fixed-end results in significant differences. Fixed-end cables should take bending stiffness influences into account.

Define the ratio of span midpoint bending moment sag of parabolas to catenaries corresponding to division of (50) and (53). Equation (53) is the bending moment sag equation for horizontal fixed-end parabolas (). Consider is the ratio of span midpoint bending moment sag of parabolas to catenaries. is illustrated in (48). is illustrated in (49).

The ratio differences between bending stiffness and solutions of parabolas and catenaries were compared using (50) and (53). Figure 8 signifies that the ratio varied with under the four different values (0.5, 1, 2, and 3). The parabola solution was not applicable because a large led to increased differences; on the contrary, a smaller led to reduced differences and a parabola solution could replace a catenary solution. When , the ratio varied between 101.9% and 103.5%. Equation (22) was a good approximation of (51). When , the ratio varied between 107.6% and 114%. When the ratio could be as high as 242.7% and (22) was an inappropriate selection.

Figure 9 shows the comparison between the parabola and catenary fixed-end bending moments. As (20) was based on a small sag linear hypothesis, varied linearly with . However, (49) varied nonlinearly.

Define the ratio of fixed-end bending moments of parabolas to catenaries corresponding to division of (20) and (49). Consider is the ratio of fixed-end bending moments of parabolas to catenaries. is illustrated in (49).

Figure 10 adopts (20) and (49) to evaluate the ratios of influences of bending stiffness on fixed-end bending moments. The figure showed a variation trend extremely similar to Figure 8. When , the ratio varied between 102.5% and 103.5%, and (20) was a good approximation of (49). When , the parabola solution demonstrated that ratio differences increased with , and when , the ratio fluctuated between 110.1% and 114.3%. When , the ratio could be as high as 250%, and the selection of (20) was an inappropriate choice.

Figure 11 shows the distribution of parabola and catenary bending moments. The figure indicates that changed the boundary layer thickness and . At the same (), increased and reduced . This resulted in a reduced thickness of boundary layers and a more drastic variation in bending moment functions. increased the difference between (23) and (52). When , the catenary solution could be replaced by the parabola solution. In addition, the two solutions of showed significant difference in values and shapes and in such a circumstance the parabola was not applicable.

#### 7. Conclusions

(1)This study proposed a novel large sag catenary bending moment expression that allows finding the sag directly from the bending moment equation and finding the WKB catenary solution using the WKB method. The matching of cable and beam segments is no longer necessary, thereby overcoming the complexity of boundary layers. Moreover, with only a single equation, the proposed method could simultaneously describe the drastic function variations inside and outside the boundary layer. The method provides a simple calculation of cables with bending stiffness and large sag and fulfills the engineering needs for tensile cable fatigue stress analysis and estimation of suspension cable sag for bridge hoisting.(2)Statistical analysis revealed a significantly different influence on the behavior of cables with bending stiffness given that the boundary condition was hinges or fixed-ends. Bending stiffness was regarded as effective when for the hinged catenary. However, the fixed-end catenary should always consider the influence of bending stiffness. When catenary , the horizontal tensile force was massive (). The parabola shows a very good approximation of catenary. It could replace the WKB catenary solution. In addition, had the ability to change the boundary layer thickness where bending moments were distributed and fixed-end bending moment. Higher lowered the fixed-end bending moment, reduced the boundary layer, and increased the drastic variation of the bending moment functions.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### References

- H. M. Irvine,
*Cable Structures*, Dover Publications, 1992. - J. J. Burgess, “Bending stiffness in a simulation of undersea cable deployment,”
*International Journal of Offshore and Polar Engineering*, vol. 3, no. 3, pp. 197–204, 1993. View at Google Scholar · View at Scopus - E. J. Hinch,
*Perturbation Methods*, Cambridge University Press, Cambridge, UK, 1991. View at Publisher · View at Google Scholar · View at MathSciNet - J. Kevorkian and J. D. Cole,
*Multiple Scale and Singular Perturbation Methods*, Springer, Berlin, Germany, 1996. View at Publisher · View at Google Scholar · View at MathSciNet - A. H. Nayfeh,
*Introduction to Perturbation Techniques*, John Wiley & Sons, New York, NY, USA, 2011. - M. S. Triantafyllou and G. S. Triantafyllou, “The paradox of the hanging string: an explanation using singular perturbations,”
*Journal of Sound and Vibration*, vol. 148, no. 2, pp. 343–351, 1991. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - P. Wolfe, “The effect of bending stiffness on inextensible cables,”
*International Journal of Engineering Science*, vol. 30, no. 9, pp. 1187–1192, 1992. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - H. M. Irvine, “Local bending stresses in cables,”
*International Journal of Offshore and Polar Engineering*, vol. 3, no. 3, pp. 172–175, 1993. View at Google Scholar · View at Scopus - D. M. Stump and W. B. Fraser, “Bending boundary layers in a moving strip,”
*Nonlinear Dynamics*, vol. 21, no. 1, pp. 55–70, 2000. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - D. M. Stump and G. H. M. van der Heijden, “Matched asymptotic expansions for bent and twisted rods: applications for cable and pipeline laying,”
*Journal of Engineering Mathematics*, vol. 38, no. 1, pp. 13–31, 2000. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - V. Denoël and E. Detournay, “Multiple scales solution for a beam with a small bending stiffness,”
*Journal of Engineering Mechanics*, vol. 136, no. 1, Article ID 006001QEM, pp. 69–77, 2010. View at Publisher · View at Google Scholar · View at Scopus - V. Denoël and T. Canor, “Patching asymptotics solution of a cable with a small bending stiffness,”
*Journal of Structural Engineering*, vol. 139, no. 2, pp. 180–187, 2013. View at Publisher · View at Google Scholar · View at Scopus - A. B. Mehrabi and H. Tabatabai, “Unified finite difference formulation for free vibration of cables,”
*Journal of Structural Engineering*, vol. 124, no. 11, pp. 1313–1322, 1998. View at Publisher · View at Google Scholar · View at Scopus