JORQUERA
JORQUERA
JORQUERA
2 (2016) 000-000
DOI: http://dx.doi.org/10.12989/scs.2016.20.1.000
(Received April 22, 2015, Revised July 14, 2015, Accepted October 20, 2015)
Abstract. Despite of the growing number of built examples, the analysis of non-symmetrical cable stayed bridges
has not received considerable attention from the researchers. In fact, the effects of the main design parameters in the
structural behavior of these bridges are not addressed in detail in the literature. To fill this gap, this paper studies the
structural response of a number of non-symmetrical cable-stayed bridges. With this aim, a parametric analysis is
performed to evaluate the effect of each of the main design parameters (the ratio between the main and the back span
length, the pylon, the deck and backstay stiffnesses, the pylon inclination, and the stay configuration) of this kind of
bridges. Furthermore, the role of the geometrical nonlinearity and the steel consumption in stays are evaluated.
Keywords: non-symmetrical cable-stayed bridges; parametric study; design parameters; steel bridge; design
guidelines
1. Introduction
The rapid bloom of cable-stayed bridges has been propitiated by the development of erection
techniques and improvements of construction materials (Svensson 2012). According to many
authors (see Agrawal 1997), this bridge type is the most economical one for bridge spans ranging
between 150 and 360 m. A major concern in the structural efficiency of a cable-stayed bridge and,
therefore, in its cost, lies in the symmetry of its longitudinal configuration, i.e., when the side span
is roughly the 40% of the main span for a three span bridge. For this reason, structurally speaking,
designers prefer symmetric designs rather than non-symmetric ones, which are limited to those
cases where the symmetric solution is discouraged for site conditions. Nevertheless, economic
boom of the last decade has deflected attention from economic cost to iconic aesthetical
appearance increasing significantly designer‟s freedom (Jorquera-Lucerga 2013). This structural
freedom has enabled the proliferation of a number of non-symmetric cable-stayed bridges that
might even go without basic structural elements such as backstays (as the Alamillo Bridge in Spain,
Casas and Aparicio 1998).
In cable-stayed bridges, all the resistant mechanisms are related, to a greater or lesser extent, to
the tensile forces of the stay cables in service. Since current cable-stayed bridges include a highly
indeterminate structural system, no unique solution for pre-tensioning stay cable forces exists
(Marchetti and Lecinq 1999). In fact, these represent a design parameter that can be tailored to
achieve an effective design for the bridge (such as full use of material strength, achievement of
target geometry or minimization of structural cost).
The stay cable forces are usually defined to achieve a certain designer‟s criterion in service, in a
stage known as the Objective Service Stage or OSS (see Lozano-Galant et al. 2012a, b). The
criteria proposed in the literature to estimate the stay forces in the OSS are reviewed by many
authors (see Guan 2000, Chen et al. 2000, Hassan et al. 2012, Lozano-Galant et al. 2013).
According to them, some of the most important criteria are as follows: (1) The Pendulus Rule (see
Virlogeux 1994, SETRA 2001). This is a simplified method to carry out a first approximation of
the stay cable forces in nearly horizontal decks. This criterion assumes that the tension of each stay
cable is close to the tension obtained when the deck is hinged at each stay anchorage. In this way,
the cable forces can be estimated as the projection in the stay direction of the averaged load
introduced into the two adjacent deck segments. (2) Minimization of the Creep Effect Criterion
(see Scotti 2003). This criterion defines the stay cable forces in the OSS to minimize the creep
effect in concrete structures with vertical stay cables. In order to analyze the creep effects, this
method applies the age-adjusted effective modulus method (Bazant 1972). The calculation of the
stay forces is performed by solving a system of equations that relates the deflections of the deck
when unitary forces are applied into the stay cables and when the permanent load is applied into
the structure. (3) Minimal Bending Energy Criterion (Du 1989). This criterion is based on the
minimization of the bending energy of the structure. Additionally, this criterion can be used to
minimize the effects of the time-dependent phenomena Lozano-Galant and Turmo (2014) and the
staggered erection of the superstructure Lozano-Galant et al. (2014, 2015). (4) Rigidly Supported
Continuous Beam Criterion (Manterola et al. 2006, Chen et al. 2000, Gimsing 1997). This
criterion assumes that the long-term behavior of a cable-stayed bridge corresponds with that of a
fictitious rigidly continuous beam. This beam is defined by removing the stay cables and adding
fictitious bearings at the bridge deck anchorages. In this way, the stay cable forces can be obtained
by projecting the vertical reactions of the corresponding fictitious supports into the stay cable
direction. (5) Zero Displacement Criterion (Lazar et al. 1972, Wang et al. 1993). This criterion
defines the stay cable forces to achieve zero deflections at certain control points of the structure. A
common criterion consists of defining zero vertical deflection at the deck-stay connection and zero
horizontal deflection at the top of the pylon in structures with backstays. To achieve this goal
(Wang et al. 1993) proposed the use of an iterative process in which the prestressing stay cable
forces are successively updated to minimize the deflections at the control points. Some other
methods such as the B-spline curve criterion (Hassam et al. 2012, Hassam 2013) are also based, to
a great extent, on the analysis of the bridge deflections. The Minimization of the Sum of Squares
Method criterion defines the stay cable forces by an optimization analysis in which the sum of the
squares of the vertical deflections along the bridge deck is minimized. On the other hand, the B-
spline curves criterion proposes the optimization, by mean of genetic algorithms, of the post-
tensioning functions to achieve minimum deck deflections. These post-tensioning functions are
defined by B-spline curves. (6) The Unit Load Method (Janjic et al. 2002, 2003). This criterion is
based on a linear system of equations that includes a degree of freedom for each stay cable force.
This system relates the bending moments at some control points for two types of load cases: unit
prestressing loads at each stay cable and the target load of the superstructure. This system of
Structural behavior of non-symmetrical steel cable-stayed bridges
equations can be directly solved by mean of an influence matrix after defining the desired moment
distribution to be achieved. A similar approach was used by Guan in terms of deflections. (7)
Optimization Criterion (Baldomir et al. 2010, Negrao and Simoes 1997, Simoes and Negrao 2000).
In this criterion, the prestressing stay cable forces are defined by the minimization of a scalar
objective function. Different tendencies are used to define these objective functions. Some
designers base their objective functions on the structural efficiency while some others base it on
the economy of the structure. To reduce the computation time of the optimization process, Lute et
al. (2009) proposed the use of genetic algorithms.
All criteria presented above might be applied to symmetric cable-stayed structures.
Nevertheless, not all of them are adequate to non-symmetric ones. This is the case of the Pendulus
rule and the Rigidly Continuous Beam Criterion as they provide no information of the backstay
forces. For this reason, sometimes the designers have to look for alternative methods. This was the
case of the design of the Alamillo Bridge (Casas and Aparicio 1998). In this structure the stay
cables forces in the main span were defined to correspond with those of an equivalent continuous
beam. Nevertheless, the weight of the inclined pylon was defined to assure that the resultant of the
stay forces remained inside the pylon. The lack of studies referring to non-symmetrical cable-
stayed bridges is not limited to the backstay forces as the effects of the main design parameters
(such as the ratio between the main and the back span length, the pylon and the deck stiffnesses,
the pylon inclination or the stay configuration) are not conveniently addressed neither.
To fill the detected gaps, this paper studies the structural behavior of highly non-symmetrical
two span cable-stayed bridges made of steel. With this aim, a number of parametric analyses
focused on each of the main design parameters (ratio between the main and the back span length,
the pylon and the deck stiffnesses, the pylon inclination of the stay configuration) are presented.
Furthermore, the effect of the geometrical nonlinearity in these structures is not significant.
This paper is organized as follows: In Section 2, the layout and the load cases in the case of
study are presented. In Section 3, the calculation and the assumptions of the study are presented. In
Section 4, the results of the parametric analyses are discussed. Finally, a set of conclusions are
drawn in Section 5.
2. Cases of study
The cases of study consist of a cable-stayed bridge whose geometry is modified to analyze the
effects of different non-symmetrical configurations. In this section, the configurations and load
cases of the analyzed non-symmetric bridges are presented.
2.1 Configurations
The main parameters that define the geometry of a two span non-symmetric cable-stayed
bridge may be defined as follows: (1) Back span, L1, (2) Main span, L2, (3) Pylon height over the
deck, H1, (4) Pylon height below the deck, H2, (5) Pylon inclination, α, (6) Backstay inclination, α1,
(7) Lower main span stay inclination, α2, (8) Cable arrangement, in which the stay cables are
distributed throughout the pylon. The studied arrangements are Harp (parallel cables) and Semi-
Fan (cables distributed throughout a certain length h, near the top of the pylon), (9) Anchorage
separation in the pylon back span, h1, (10) Anchorage separation in the pylon main span, h2, (11)
Number of stay cables in the back span, n1, (12) Number of stay cables in the main span, n2, (13)
Juan J. Jorquera-Lucerga, Jose A. Lozano-Galant and Jose Turmo
Fig. 1 Main parameters of a non-symmetrical cable-stayed bridge with different stay configurations
Separation of the backstay, s0, (14) Anchorage separation in the deck back span, s1, (15)
Anchorage separation in the deck main span, s2. (16) Separation of the deck last stay, s3. All these
parameters are summarized in Fig. 1.
The assumptions considered in this study might be summarized as follows: (1) Pylon-deck
connection: As considered in many of actual bridges, a configuration pioneered by Severins
Bridge (1960) over the River Rhine in Cologne: see e.g., Zárate-Brazo Largo cable-stayed bridge
(Svensson 2012) no pylon-deck connection is considered. (2) Valid range of stay inclinations: To
assure an adequate efficiency of the backstays, α1 and α2 must remain into a certain range of
validity. In this study, this range has been assumed between 30° and 60°, although smaller values
outside this range could also be considered (e.g., Menn (1990) considered minimum stay
inclinations of 22°). (3) Fixed parameters in non-symmetrical configurations: To homogenize the
characteristics of the analyzed bridges, all the examples share the following fixed parameters L2 =
100 m, s0 = s1 = s2 = 6 m, s3 = 7 m, H2 = 15 m and n2 = 15. (4) Parameters defining the non-
symmetrical configuration: The non-symmetrical configurations are defined in terms of the fixed
parameters and fixing α, α1 and h. Firstly, the pylon top can be obtained by intersecting the last
main span stay and the pylon. In this way, H2 can be calculated in terms of α1 and α as presented in
Eq. (1). This height can be used to calculate L1 as presented in Eq. (2). In order to satisfy s0 = 6 m
and S1 = 6 m with an integer number of stay cables, the admissible L1 are restricted to 10 m, 21 m,
33 m, 51 m, 63 m, 75 m, 81 m and 99 m. For each of these values of L1, n1 can be estimated from
Eq. (3). Finally, h1 can be obtained in terms of h and n1 as presented in Eq. (4). The calculation of
h2 can be carried out by substituting n1 to n2 in Eq. (4). (5) The effects of the staggered erection of
the superstructure are neglected. (6) Steel relaxation is neglected. The simulation of this
phenomena can be found in Gunaydin et al. (2014) and Atmaca and Ates (2012).
tan 𝛼1 ·𝐿2
𝐻2 = (1)
1 − tan 𝛼 ·tan 𝛼1
1
𝐿1 = 𝐻2 tan 𝛼 + < 𝐿2 (2)
tan 𝛼1
Structural behavior of non-symmetrical steel cable-stayed bridges
𝐿 − 𝑆3
𝑁2 = + 1 − 𝑁1 (3)
𝑆1
ℎ
ℎ1 = (4)
𝑁1
The valid cases that satisfy stay cables efficiency condition for L1=100m are located into the
shadowed region of Fig. 2(A) in terms of L1, α and α1. To reduce the number of bridge geometries
available for the parametric analysis only some α (-20°, -10°, 0°, 10° and 20°) are studied, where a
negative angle represents a pylon leaning backwards and a positive angle represents a pylon
leaning forward the main span. The agreement between the studied α and the admissible L1 results
in the 23 analyzed cases summarized in Fig. 2(A). Each of these geometries is analyzed with the
two alternative pylon and deck sections presented in Fig. 2(B). Differences between both sections
(Stiff and Flexible) refer to their stiffness. The sections stiff and flexible sections in the steel pylon
are 4000×2000 and 2000×2000 mm with 30mm thickness, respectively. In the case of the deck, the
Stiff and Flexible sections are 6000×2000 and 10000×500 mm with 30 mm thickness. It is
remarked that all analyzed pylons and decks are assumed to be made of steel. Therefore, time-
dependent phenomena that occur in concrete bridges are not considered. The mechanical
Fig. 2 Cases analyzed in the parametric study. (a) Valid and analyzed cases in the parametric study; and
(b) Analyzed pylon and deck cross-sections
Juan J. Jorquera-Lucerga, Jose A. Lozano-Galant and Jose Turmo
properties of both the deck and the pylon are summarized in Table 1. The steel Young‟s modulus
is 210000 MPa for the deck and the pylon and 195000 MPa for the stay cables. The ultimate stay
strength, fGUTS, is 1860 MPa. These values are obtained from (Eurocode 2005). The areas of the
stay cables are not included in this table as they are calculated for each cable in each model.
The configurations of the studied examples vary to analyze the influence of the following
variables: (1) Pylon inclination α (ranging from -20° to 20°), (2) Deck cross section (Stiff or
Flexible in Fig. 2(B)), (3) Pylon cross section (Stiff or Flexible in Fig. 2(B)), and (4) Stay cable
configuration (Semi-Fan or Harp in Fig. 1(A)). In addition to these analyses, a study of the effects
of the geometrical nonlinearity is presented. All these variations results in a number of 368
analyzed structures.
In this section after summarizing the calculation procedure, the main simulation assumptions
are presented.
3.1 Calculations
Fig. 3 Envelope of bending moments in service. The envelope proposed in Walther et al. (1999) is
presented by a dotted line and the one used in this paper by a continuous one
solid dash, with the bending moment envelope obtained when every live load location is
considered, in red solid dash. This figure shows that both envelopes present the same stresses in
the pylon and in the proximities of the deck abutments (where the maximum deck bending
moments are found). Nevertheless, this is not the case at the inner part of the deck where the
accuracy of the Walther‟s envelope is lost (as higher bending moments are obtained by the red
envelope). In view of these results, and to assure an accurate simulation of the whole bridge, the
stress and deflection envelopes calculated in this paper have been obtained considering all possible
live load locations. To do so, the analyzed FEMs include a number of load cases equal to the
number of deck beam elements.
3.2 Assumptions
(1) Geometry of the FEMs and load cases. The geometry of the FEMs and the load cases
analyzed are presented in Section 2.1 and 2.2, respectively.
(2) Definition of all stay cable cross sections but the backstay one. The stay cable cross
sections must fulfill two criteria simultaneously. The first one is the limitation of the
maximum stresses. To prevent fatigue failure, the stay cables are designed to avoid
exceeding the 45% of the ultimate stay strength, fGUTS under the characteristic combination
(see SETRA 2001). The second one is the limitation of increment of stresses. To prevent
fatigue failure, the maximum increment of stress in stay cables due to live loads, 𝛥𝜎max ,
must be controlled. In this paper, this value has been fixed to 200 MPa. The amount of live
loads that must be considered for fatigue checking are specified in codes. In this paper, the
fatigue load considered is equal to the 100% of the characteristic live load. The area of
every stay cable, Ac, can be expressed as
𝑃𝐿 + 𝐿𝐿 ·𝑠 𝐿𝐿 ·𝑠
𝐴𝑐 = max , (5)
sin 𝛼 ·0.45 ·𝑓𝐺𝑈𝑇𝑆 sin 𝛼 ·𝛥𝜎𝑚𝑎𝑥
Juan J. Jorquera-Lucerga, Jose A. Lozano-Galant and Jose Turmo
where PL and LL are the Permanent and Live loads, s is the cable spacing and α is the
cable inclination defined in Fig. 1. This procedure is based on the method of the tributary
area proposed by Svensson (2012).
(3) Definition of stay cable forces in all stays but the backstays due to permanent loads.
According to Svensson (2012), the Rigidly Supported Continuous Beam Criterion has
been used to define the axial load of the stay cables due to permanent loads.
(4) Definition of backstay area and force due to permanent loads. The definition of the area
and force of the backstay is not as straightforward as in the case of the rest of the stay
cables. The difficulty in defining these elements comes from the fact that backstay force
and stiffness are strongly linked with bridge forces and deflections. This problem has to be
solved by an iterative process, in which each iteration represents a different FEM of the
whole bridge with a certain backstay area. In this model, the axial force of the backstay is
defined to set the bending moment at the bottom of the pylon to zero when the permanent
loads and the axial forces in the rest of stays calculated according to the preceding
paragraph are applied. Once defined the axial force of the backstay for permanent loads,
the variation of stresses, 𝛥𝜎𝑚𝑎𝑥 , and the maximum tensile stresses in the backstay are
calculated using a FEM model. If any of these values exceeds the design thresholds, a new
area of the backstay is defined. This new area is introduced in the FEM that is used in the
following iteration.
This section studies the effects of the main design parameters of the non-symmetrical cable-
stayed bridges in the structural behavior of their structural elements (stays, pylon and deck). The
analyzed design parameters are: the back span length, the pylon inclination, the pylon stiffness, the
deck stiffness, the backstay stiffness, the stay arrangement and the geometrical nonlinearity, After
presenting the results of these analyses, they are discussed in detail.
In this section, the results of the analyzed design parameters are presented. It is remarked that
the geometry and the stay areas of all the proposed models are defined according to the procedures
described in Sections 1 and 3, respectively. The Stiff and Flexible sections for the pylon and deck
are described in Fig. 2.
the envelopes of vertical deflections in the deck, vD, (Fig. 5(B)), the envelopes of bending
moments in the pylon, MP (Fig. 5(C)) and the envelopes of horizontal deflection in the pylon, uP,
(Fig. 5(D)).
Structural behavior of non-symmetrical steel cable-stayed bridges
symmetrical cable-stayed bridges with a vertical pylon (α = 0°) and six different back span lengths
(L1 = 33 m, 51 m, 63 m, 75 m, 81 m and 99 m). With this aim, the effects of two different pylon
stiffnesses (Stiff and Flexible) are analyzed. The stiffness of the deck is assumed as Flexible.
Fig. 7 summarizes the area of the cables, Ac, (Fig. 7(A)), the stress variations in cables due to
live loads, Δσc, (Fig. 7(B)), the envelope of bending moments in the deck, MD (Fig. 7(C)), the
envelope of vertical deflections in the deck, vD, (Fig. 7(D)), the envelope of bending moments in
the pylon, MP, (Fig. 7(E)), and the envelope of horizontal deflections in the pylon, uP, (Fig. 7(F)).
In this figure continuous lines indicate the results for the Stiff pylon and the dotted ones the results
of the Flexible pylon.
are compared with the reference values (γBS = 1) in Fig. 9. This figure summarizes the envelope of
bending moments in the deck, MD, (Fig. 9(A)), the envelope of vertical deflections in the deck, vD,
(Fig. 9.B), the envelope of bending moments in the pylon, MP, (Fig. 9(C)), and the envelope of
horizontal deflections in the pylon, uP, (Fig. 9(D)).
cable-stayed bridges. The analyzed structures include a vertical pylon (α = 0°) and six different
back span lengths (L1 = 33 m, 51 m, 63 m, 75 m, 81 m and 99 m).
Fig. 11 summarizes the cable areas, Ac, (Fig. 11(A)), the envelope of increments of axial
stresses in the stays, Δσc, (Fig. 11(B)), the envelope of bending moments in the deck, MD, (Fig.
11(C)), the envelope of vertical deflections in the deck, vD, (Fig. 11(D)), the envelope of bending
moments in the pylon, MP, (Fig. 11(E)), and the envelope of horizontal deflections in the pylon, uP,
(Fig. 11(F)). In this figure, continuous and dotted lines are used to identify the results of the
analyses with and without P-δ effects, respectively.
4.2 Results
In this section, the results of the parametric analysis presented in the preceding section are
discussed in detail.
values of bending moments, especially at its bottom, and higher deflections, at its top, appear
when L1 decreases. The lack of symmetry of the bending moment diagram at the bottom of the
pylon is a consequence of the lack of symmetry of the span distribution. This behavior is
illustrated in the models with α = -20°, where the maximum sagging moments in the deck vary
from 11 MNm for L1 = 10 m to 6 MNm for L1 = 99 m, and in the pylon, these values vary from 42
MNm for L1 = 10 m to 21 MNm for L1 = 99 m.
4.3 Discussion
since when heavy vehicles are not considered the area of the cables are underestimated. For these
reasons, when only distributed loads are considered, it seems reasonably to reduce the percentage
of the cable strength from 45% to 35% or a similar value. On the other hand, the oscillation of stay
stresses criterion might be the determining criterion in light cable-stayed bridges. In these
structures, the ratio of the live load considered in the oscillation is defined by the structural code.
This ratio is usually lower than the 100% assumed in the presented analyses.
12(D), where the weight of the backstay cables, the rest of the cables and of all the cables are
compared. Fig. 12(D) can be expressed in terms of kg/m2 as presented in Fig. 12(E). In this figure
the minimum steel weight is obtained when L2 is near 40 m (L2/L1 = 0.40).
For all the pylon inclinations the minimum value is obtained for α = -20°, although the
differences are not significant. In symmetrical bridges the optimal pylon inclination is α = 10°. On
the other hand, when 0.55 ≤ L2/L1 ≤ 0.75 the optimal pylon corresponds with α = 0°. It is remarked
that differences between the different cases are not significant enough to advice one or another
pylon inclination. This is not the case in structures with a short L2 span as in this structures the
steel weight is significantly reduced when α < 0°.
Fig. 12 Results backstay areas (A); variation of stay forces due to live loads (B); variation of stay
stresses due to live loads (C); steel consumption (D); and steel in kg/m2 (E)
Structural behavior of non-symmetrical steel cable-stayed bridges
in Fig. 12(D) and 12(E). A reference value might be considered as L2/1000 for the frequent live
load [IAP 11], which corresponds to the Load Model 4 (5 kN/m2) defined at Eurocode 0, factored
by ψ = 0,4. Nevertheless, this ratio might be reduced in agreement with the client. For example, in
the Sunniberg Bridge this value was as low as L2/400 (Marchetti and Lecinq 1999).
5. Conclusions
This paper presents a parametric study addressing the structural response of asymmetrical
cable-stayed bridges. In this analysis, the effects of the main design parameters (the pylon
inclination, the back span length, the stiffness of the deck, pylon and backstay and the stay
configuration) are compared. This study includes the role of the geometrical nonlinearity. Despite
of the growing number of built examples, this structural type has deserved little attention from
researchers.
The parametric analyses carried out show that: (1) In the vast majority of analyzed bridges, the
stress oscillation criterion stands as the most adverse criterion to define the backstay area. On the
contrary, the rest of the stays are usually defined by limiting its axial strength. (2) The higher the
pylon inclination towards the main span, the higher the backstay area. The minimal backstay area
is obtained when the back span length represents a ratio between the 20 and 50% of the main span
length. (3) In most analyzed structures, variations in steel weight between structures with different
pylon inclinations are not significant enough. Nevertheless, this is not the case in structures with a
short backspan. In these structures, the steel weight is significantly reduced when the pylon is lent
towards the backspan. (4) The backspan length plays an important role in the structural response of
asymmetrical bridges. The shorter this length the worse the structural response in the main span is.
(5) Leaning the pylon towards the back span for a certain back span length, worsen the structural
behavior of the main span. (6) Deck deflections are not reduced in the same proportion when the
deck bending stiffness is increased. For this reason, the use of stiffer decks is not advised. (7)
Increasing the backstay stiffness is more efficient to improve the deck behavior than increasing the
bending stiffness of the pylon. (8) Generally speaking, the harp configuration is less efficient than
the semi-fan one. (9) The geometrical nonlinearity does not play an important role in the structural
response of the asymmetrical bridges analyzed in this study.
As a final conclusion, generally speaking, asymmetrical bridges present a less efficient
structural behavior, since the disadvantages of symmetrical bridges are enhanced, especially for
very short back spans. As showed above, factors like the area and inclination of the backstay cable,
the cable arrangement (harp or semi-fan), the stiffness and inclination of the pylon, and the
stiffness of the deck, are of major importance. With no doubts, this kind of bridges need further
studies in order to, for example, establishing criteria to design the area of the backstay cables.
These and other problems will be addressed in the near future by the authors.
Acknowledgments
The authors would like to thank the former students D. Javier Sánchez Mateos Enrique, D.
Mario Pérez Sánchez y D. Alberto de la Torre Palencia for contributing to the initial studies of this
paper. The authors thank the Spanish Ministerio de Economia y Competitividad and the FEDER
funds for the funding provided through the research grant BIA2013-47290-R directed by Jose
Turmo.
References
Agrawal, T.P. (1997), “Cable-stayed bridges-parametric study”, J. Bridge Eng., 2(2), 61-67.
Atmaca, B. and Ates, S. (2012), “Construction stage analysis of three-dimensional cable-stayed bridges”,
Structural behavior of non-symmetrical steel cable-stayed bridges
CC