For building a reasonable finite element geometric model, a method is proposed to determine the strain-free configuration of the multispan cable. The geometric conditions (the end conditions and the unstretched length conditions) as constraints for the configuration of multispan cable are given. Additionally, asymptotic static equilibrium conditions are given for determining the asymptotic shape of the multispan cable. By solving these constraint equations, a set of parameters are determined and then the strain-free configuration of multispan cable is determined. The method reported in this paper provides a technique for building reasonable FEA geometric model of multispan cables. Finally, a three-span cable is taken as example to illustrate the effectiveness of the method, and the computed results are validated via the software ADINA.
1. Introduction
The cable structures are widely used in the engineering, such as bridge, building, and power transmission system [1, 2]. In [1], a new force method is proposed for analysing the dynamic behaviour of flat-sag cable structures. The accepted dynamic model of such cables reduces to a partial differential equation and an integral equation. The support reaction forces are considered as excitations, allowing D’Alembert’s solution to be used. In this way, single or multispan cables have been developed in the form of a single-degree-of-freedom system in terms of the additional dynamic tension. Finally, two examples are presented to illustrate the accuracy of the proposed force method for single and multispan cable systems subjected to harmonic forces. In [2], a three-dimensional modeling procedure is proposed for cable-stayed bridges with rubber, steel, and lead energy dissipation devices. In this study, the cable structure can be simplified as a four-node isoparametric cable element. The multispan cables and towers are the basic structures in the power transmission system. Under the wind load, the multispan cable is apt to generate large displacements because of its high flexibility and large span. So the multispan cable exhibits strong geometrical nonlinearity. Usually, the dynamic behavior of such cable structures is analyzed on the platform of the finite element analysis (FEA) software. To build a reasonable FEA model, two alternative configurations of the cable structures need to be determined in the modeling process; one is the initial static equilibrium configuration and the other is the assumed strain-free configuration. As a prestressed state, the static equilibrium configuration of a cable is the starting point of the subsequent dynamic analysis, whereas the strain-free configuration of a cable is only an assumed state which does not exist in the real world due to the ubiquitous gravity and it is the starting point for static and dynamic analyses. If the FEA geometric model is built based upon the static equilibrium configuration of the cable, the prestrain or prestress in each element caused by gravity should be specified one element by one element; this is always cumbersome. The gravity load is included in the equilibrium configuration. The FEA geometric model contains two gravity loads in this way, so that the model does not match the actual cable structure. The gravity load is not included in the strain-free configuration. Alternately, if the geometric model is built based on the strain-free configuration, the prestrain or prestress in each element needs not to be specified in advance and they are set up naturally in the nonlinear static analysis via the FEA. After that, the gravity load applied in the model; the system achieves equilibrium. Therefore, the assumed strain-free configuration is more preferable to the static equilibrium configuration for building the FEA geometric model of multispan cable.
For the static and dynamic analysis of cable structures, Irvine [3] systematically summarized the classical achievements reached up to 1981 and the static equilibrium shape of a cable is described as an elastic catenary. For the modeling of cable structures, in [4], a method for modeling cable supported bridges for nonlinear finite element analysis is presented in this paper. A two-node catenary cable element, derived using exact analytical expressions for the elastic catenary, is proposed for the modeling of cables. Based upon the elastic catenary, Such et al. [5] proposed a method for determining the static equilibrium configuration of arbitrary three-dimensional cable structures subject to gravity and point loads. Srinil et al. [6, 7] formulated a system of nonlinear partial differential equations describing the large amplitude three-dimensional free vibrations of inclined sagged elastic cables. Ai and Imai [8] proposed an iterative scheme to find the static equilibrium shape of beam-cable mixed system. In [9], based on exact analytical expressions of elastic catenary, Thai and Kim presented a catenary cable element for the nonlinear static and dynamic analysis of cable structures. In [10], Vu et al. presented a spatial catenary cable element for the nonlinear analysis of cable-supported structures and proposed an algorithm for form-finding of cable-supported structures. Rienstra [11] derived the partial differential equations and boundary conditions satisfied by the static equilibrium configuration of multiple spans coupled via suspension strings. In [12], Impollonia et al. obtained the deformed shape of the elastic cable in closed form for the cases of uniformly distributed load and multiple-point forces. Greco and Cuomo [13] proposed a method for obtaining an exact configuration of slack cable nets by using the exact expressions of the equilibrium derived from the equation of the catenary.
All the above papers focused on the determination of static equilibrium configuration of the cable structures. However, the present paper concentrates on determining the strain-free or unstretched configuration of multispan cable. Firstly, based upon the elastic catenary theory, a set of geometric boundary conditions, cable length condition, and static equilibrium conditions are derived for determining the configuration of the multispan cables when they are approaching asymptotically to the strain-free state. These constraint conditions constitute a system of nonlinear algebraic or transcendental function equations. Secondly, these nonlinear equations are solved and one finds a set of parameters which determine the unstretched configuration of the multispan cable. Finally, an example of three-span cable is given to illustrate the method.
2. The Governing Equations of Multispan Cable
First we assume that (1) the bending stiffness of cable is negligible; (2) the material of cable is of Hookian type; that is, its constitutive relation is linear elastic; and (3) only the self-weight as the static force is applied to the cable structure. Under these assumptions, the governing equation for the cable in its static equilibrium state is as follows [6]:(1)EAε01+y0′2′=0,EAε0y0′1+y0′2′+ωc1+y0′21+ε0=0.Here, the “′” denotes the derivative with respect to x, ε0 denotes the tensional strain in the cable when it is in static equilibrium state, E is Young’s modulus, A is the sectional area of the cable and assumed to be constant, y0 represents the static equilibrium shape of the cable, and ωc is the cable weight per unit length in the unstretched state.
The tensional strain ε0 is related to the horizontal tensional force TH as follows:(2)ε0=THEA1+y0′2.According to (1), the final static equilibrium governing equation of the cable can be rewritten in the following form:(3)1+THEA1+y0′2y0′′+ωcTH1+y0′2=0.This is the so-called elastic catenary equation. Being different from the traditional catenary equation, it takes account of the extensibility of the cable. Under the dynamic wind load, the multispan cable usually vibrates nonlinearly due to its large displacement, so the finite element analysis is preferable to investigate the static and dynamic behavior. To build the FEA model, firstly the geometric model of multispan cable is needed. Therefore, one of the initial configurations, the strain-free state or the static equilibrium state, should be determined. As stated in the Introduction, in the following, we only discuss how to determine the strain-free configuration of the multispan cable.
3. Determination of the Strain-Free Configuration
Determining the strain-free configuration of the cable based on the data of static equilibrium configuration can be viewed as an inverse problem in nonlinear structural mechanics. The unstretched state of cable is an assumed state and it can be viewed as an asymptotic state that corresponds to EA→∞ (which means the cable is inextensible), and ωc/TH→-a (note that ωc is negative), so the equation for cable in unstretched state is as follows:(4)y′′-a1+y′2=0.This is the rigid catenary equation. Its solution can be found as(5)yx=1acoshax+c1+c2.
For the three-span cable, as illustrated in Figure 1, the shape of the left span cable in its unstretched state can be expressed as(6)y1x=1a1cosha1x+c11+c12.Similarly, the shape of the middle span cable can be expressed as(7)y2x=1a2cosha2x+c21+c22and the shape of the right span cable is expressed as(8)y3x=1a3cosha3x+c31+c32.
The configuration of three-span cable.
There are 9 unknown parameters appeared in the above expressions. Determination of the strain-free configuration of the three-span cable is equivalent to determining these unknown parameters, so the following constraint equations can be built.
For the left span cable,(9)End conditions:y1xx=0=H1,y1xx=la=ha,L1-la2+H2+r1-ha2=r12.Here, the meaning of the first and second equations of (9) is obvious. In engineering, the rod AP1 (insulator) is very stiff so it can be reasonably assumed as rigid; therefore, the third equation of (9) holds too. Figure 1 can be referred to for the meaning of the symbols L1, H1, H2, la, ha, and r1 in the above equations. In addition, the following condition about the length of the original unstretched cable can be built:(10)Length condition:∫0la1+y1′2dx=s01.Here, s01 represents the length of the left span cable in its unstretched state and can be expressed via the data of the static equilibrium configuration.
It should be pointed out that, in some cases (e.g., in the example given later), the original length s01 of the cable in its unstretched state may be less than the straight distance AP1, so, in order to build a reasonable FEA geometric model, the position of the hanging point P1′ must be determined.
Note that the function y1(x) is in the form of cosh; the above integral actually can be given explicitly. To find the length s01, the following equations can be constructed (see [3], page 18):(11)THs01EA+THs01W1sinh-1V1TH-sinh-1V1-W1TH=L1,W1s01EAV1W1-12+THs01W11+V1TH21/2-1+V1-W1TH21/2=H1-H2,where W1(=ωcs01) is the self-weight of the left span cable and TH and V1 are the horizontal and vertical reaction forces at the support point A, respectively. From the above two equations, the unknown V1 and s01 can be solved.
Being similar to (9), for the middle and right span cables, the following constraint equations can also be constructed.
For the middle span cable,(12)End conditions:y2xx=la=ha,y2xx=L1+L2+L3-lb=hb,L3-lb2+H3+r2-hb2=r22,Length condition:∫laL1+L2+L3-lb1+y2′2dx=s02.For the right span cable,(13)End conditions:y3xx=L1+L2+L3-lb=hb,y3xx=L1+L2+L3=H4,(14)Length condition:∫L1+L2+L3-lbL1+L2+L31+y3′2dx=s03.Again, Figure 1 can be referred to for the meaning of the symbols H3, H4, L2, L3, lb, hb, and r2 in the above equations. s02 and s03 represent, respectively, the lengths of middle and right span cables in their unstretched state.
Similar to (11), the determining equations for s02 are as follows:(15)THs02EA+THs02W2sinh-1V2TH-sinh-1V2-W2TH=L2,W2s02EAV2W2-12+THs02W21+V2TH21/2-1+V2-W2TH21/2=H3-H2and the equations for s03 are as follows:(16)THs03EA+THs03W3sinh-1V3TH-sinh-1V3-W3TH=L3,W3s03EAV3W3-12+THs03W31+V3TH21/2-1+V3-W3TH21/2=H3-H4.Besides (9), (10), and (12)–(14), two other equations describing the equilibrium conditions of forces in the limiting state from P1→P1′ and P2→P2′ are given as follows:(17)At pointP:T1+T2·τ=0.Referring to Figure 2, the unnormalized directional vector n can be written in its coordinate-component form as(18)n=la-L1,ha-H2+l1and the unnormalized directional vector τ can also be given in its coordinate-component form as(19)τ=ha-H2+l1,L1-la.The tensional forces T1 and T2 in the left span and middle span cables can be written in their component form as(20)T1=-T111+y1′2,y1′1+y1′2x=XP1′,T2=-T211+y2′2,y2′1+y2′2x=XP1′.On the other hand, in the limit case (ωc→0), the following relations hold:(21)ωcT1/1+y1′2=-a1,ωcT2/1+y2′2=-a2.Substituting these relations into (13), one finally gets(22)1a2-1a1H2+l1-ha+la-L1y2′a2-y1′a1x=la=0.Similar to (22), at point C the following relation holds:(23)1a3-1a2H3+l2-hb+L3-lby3′a3-y2′a2x=L1+L2+L3-lb=0.Theoretically, the strain-free configurations of cable structures with any number of spans can be determined in the same way as described above. On the other hand, there are infinity of many strain-free configurations for the multispan cables; each of them satisfies the geometric constraints, that is, the end conditions and the length conditions. By imposing the force equilibrium conditions, we obtain actually the asymptotic strain-free configuration of the cable. In a sense, such an asymptotic configuration (from the static equilibrium state to the strain-free state) is optimal, because it is most similar to the static equilibrium shape of the cable, so if one performs the FEA starting from this configuration, the convergence is expected to be most rapid.
The forces at the suspension point.
For determining the configuration of the three-span cable in its unstretched state, totally, there are 13 unknowns including ai, cij (i=1⋯3,j=1,2) and ha, hb, la, and lb and 13 equations including (9)-(10), (12)–(14), and (22)-(23). Because of the complex nonlinearity, these equations can only be solved numerically. In the next section, an engineering example is given as illustration.
4. An Example
To illustrate the proposed method for determining the strain-free configuration of multispan cable, in this section a real three-span transmission line system is taken as example; its data are listed as shown in Table 1.
The parameters of three-span cable.
Parameter
Value
H1
40
H2
206
H3
210
H4
45
L1
580
L2
1910
L3
570
TH
91378
E
6.3e+10
A
6.6655e-4
ρ
4524.1
ωc
29.0668
Substituting these data into the above equations and then using the Newton-Raphson method, one can find a set of solutions as shown in Table 2.
The solution of three-span cable.
Parameter
Value
a1
0.2591e-3
a2
0.2791e-3
a3
0.2590e-3
c11
802.4259
c12
−3903.9257
c21
−1527.5928
c22
−3502.7267
c31
3879.8891
c32
−3903.5631
ha
206.0680
hb
210.0671
la
578.8362
lb
1912.3203
The unstretched lengths of the three span cables are found to be s01=602.7090463, s02=1935.110206, and s03=592.8032497, respectively. Knowing the values of these parameters, the asymptotic configuration of the three-span cable in its strain-free state can be described by (6)–(8).
Via the above computed results, the original unstretched configuration of the three-span cable can be easily plotted as shown in Figure 3. Starting from this configuration, the geometric model for FEA can be easily built, and the corresponding static equilibrium configuration can be further determined by using FEA software, for example, the ADINA. The static equilibrium configuration via the ADINA is also shown in the same figure. The largest difference in height between the unstretched and the static equilibrium configurations is about 17.75 m.
The unstretched and the static equilibrium configurations of the three-span cable.
On the other hand, to see clearly the difference between the original unstretched and the static equilibrium configurations, a zoom-in of the local configurations of the three-span cable near the left insulator is shown in Figure 4.
The magnification of local configuration near the left insulator.
5. Conclusions
In this work, a method is presented for determining the asymptotic strain-free configuration of the multispan cable. The constraints of the geometric conditions and the asymptotic static equilibrium conditions are derived, and the determining equations for the asymptotic strain-free configuration are obtained. A set of parameters which describe the shape of the multispan cable are determined. The proposed method provides a way to build reasonable FEA geometric model of multispan cables. At last, a three-span cable system is taken as example to illustrate the effectiveness of the method; the results are validated via the ADINA.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
MaX.ButterworthJ. W.A force method model for dynamic analysis of flat-sag cable structures200916662363510.3233/SAV-2009-04932-s2.0-74349090266AliH. M.Abdel-GhaffarA. M.Seismic passive control of cable-stayed bridges19952425927210.1155/1995/918721IrvineH. M.1981Boston, Mass, USAMIT PressKaroumiR.Some modeling aspects in the nonlinear finite element analysis of cable supported bridges199971439741210.1016/S0045-7949(98)00244-22-s2.0-0032664745SuchM.Jimenez-OctavioJ. R.CarniceroA.Lopez-GarciaO.An approach based on the catenary equation to deal with static analysis of three dimensional cable structures20093192162217010.1016/j.engstruct.2009.03.0182-s2.0-68549133413SrinilN.RegaG.ChucheepsakulS.Large amplitude three-dimensional free vibrations of inclined sagged elastic cables200333212915410.1023/a:10260192229972-s2.0-0142153908SrinilN.RegaG.ChucheepsakulS.Three-dimensional non-linear coupling and dynamic tension in the large-amplitude free vibrations of arbitrarily sagged cables20042693–582385210.1016/s0022-460x(03)00137-82-s2.0-0345870011AiM.ImaiH.A shape-finding analysis of suspended structures on the displacement-method equilibrium2004211274110.2208/jsceseee.21.27s2-s2.0-8344259871ThaiH.-T.KimS.-E.Nonlinear static and dynamic analysis of cable structures201147323724610.1016/j.finel.2010.10.0052-s2.0-78449243488VuT.-V.LeeH.-E.BuiQ.-T.Nonlinear analysis of cable-supported structures with a spatial catenary cable element20124355836052-s2.0-84866865427RienstraW.ManleyJ.McKeeS.OwensD. R.Nonlinear free vibrations of coupled spans of suspended cablesProceedings of the 3rd European Conference on Mathematics in Industry1988133144ImpolloniaN.RicciardiG.SaittaF.Statics of elastic cables under 3D point forces20114891268127610.1016/j.ijsolstr.2011.01.0072-s2.0-79952315440GrecoL.CuomoM.On the force density method for slack cable nets201249131526154010.1016/j.ijsolstr.2012.02.0312-s2.0-84860722735