Overall Cost Optimization of Prestressed Concrete Bridge Using Genetic Algorithm
Overall Cost Optimization of Prestressed Concrete Bridge Using Genetic Algorithm
DOI 10.1007/s12205-013-0355-4
www.springer.com/12205
···································································································································································································································
Abstract
Overall cost optimization of Prestressed Concrete (PC) bridges is investigated in this study. The purpose of this study is to
determine the optimum span number and optimum cross-sectional properties of multi-span bridges. Considered bridge superstructure
is constituted by adjacent simply supported pretensioned prestressed I-girders. Also, considered bridge substructure is constituted by
single-column piers and rectangular spread footings. Span number, cross-section dimensions of prestressed girders and the area of
Prestressing Steel (PS) are considered as design variables. PC girders, piers and footings are designed according to AASHTO
Standard Specifications for Highway Bridges. A modified hybrid Genetic Algorithm (GA) is used for the optimum design. Working
stress, ultimate strength, ductility limits, deflection, and geometry constraints are considered. Total cost of the bridge is taken as
optimality criterion. A computer program is coded to perform optimum design and numerical examples from the application are
designed. One of these examples is used in this paper. It is concluded that GA can be effectively used in the overall cost optimization
of PC bridges.
Keywords: prestressed concrete, genetic algorithm, optimization, bridge
···································································································································································································································
o
*Assistant Professor, Dept. of Civil Engineering, Namik Kemal University, 59860 Çorlu, Tekirdag, Turkey (Corresponding Author, E-mail: zay-
[email protected])
**Professor, Dept. of Civil Engineering, Yildiz Technical University, 34210 Esenler, I· stanbul, Turkey (E-mail: [email protected])
− 769 −
Zekeriya Ayd i n and Yusuf Ayvaz
according to shape of the valley. Therefore, shape of the valley is The bridge girders are considered as adjacent. So, slab thickness
considered as a design parameter in this study. The other and slab reinforcement are not taken as a design variable.
important design parameters are the total length of the bridge, the Similarly, cross section properties of the abutments are also not
width of the bridge, live load, material properties and unit prices taken as design variables. Therefore, the costs of the slab and
of the materials. abutments are not considered in the calculation of objective
The purpose of this study is to find out the minimum total cost function.
of the bridge. The total cost of the bridge includes the cost of the
PC girders, the cost of the transverse beams, the cost of the piers, 2.2 Design Variables and Design Parameters
and the cost of the foundations. Optimum span number, optimum Design variables considered in this study are given below:
PC girder number in bridge cross-section, optimum shape of the Number of piers (np)
●
PC girders and optimum number of the PS tendons are Number of PC girders (ng)
●
investigated to minimize the total cost. The total of 33 constraints Height of PC I-girder (h)
●
is considered according to AASHTO Standard Specifications for Bottom flange width (bb)
●
Optimum span number, girder number and shape of girders are Top flange thickness (tt)
●
determined in one step. It is very difficult to solve such a Bottom flange thickness (tb)
●
complicated problem with usual optimization techniques. GA is Top sloped thickness (pt)
●
used to solve this optimum design problem. A computer Bottom sloped thickness (pb)
●
program is coded in BASIC to perform the optimum design The number of PS tendons in each girder (nps)
●
materials, loads, soil bearing capacity, site altitudes, and the other q
g31 = ----------- – 1 (8)
necessary parameters. qallow
Where q is the maximum stress for the soil under the spread
2.3 Constraints
footing, and qallow is the allowable stress of the soil. Ultimate
Constraints for bridge superstructure are investigated according
shear strength constraints and ultimate flexural strength constraints
to AASHTO Standard Specifications (AASHTO, 2002). The
for the footings are given below, respectively.
total of 28 constraints is considered for superstructure. These
constraints can be categorized into six groups as given below: Vu, f
g32 = ----------
-–1 (9)
1) Flexural working stress constraints for girders φ Vc, f
2) Ultimate flexural strength constraints for girders
Mu, f
3) Shear constraints (working stress and ultimate strength) for g33 = -----------
-–1 (10)
φ Mn, f
girders
4) Deflection constraints Where Vu,f, and Mu,f are the factored shear force and moment at
5) Ductility constraints the footings, respectively, Vn,f and Mn,f are the nominal shear
6) Geometrical constraints strength and the nominal moment strength of the footings,
respectively. In the equations given above, φ is the strength
Detailed information about these constrains can be found in the
reduction factor. The value of φ is taken as 0.70 for the Eq. (6),
studies realized by Aydin and Ayvaz (2010), and Aydin (2006).
0.85 for the Eq. (7) and 9, 0.90 for the Eq. (10).
The total of 5 constraints for bridge substructure is also considered
according to AASHTO Standard Specifications (AASHTO, 2002).
3. Optimum Design with Genetic Algorithm
These constraints are given below in addition to the total of 28
constrains mentioned above:
GA is an optimization method that uses the principles of
1) Ultimate flexural strength constraints for piers
Darwin’s natural selection theory (Holland, 1975; Goldberg,
2) Ultimate shear strength constraints for piers
1989). GAs use a group of initial solutions and evaluate them to
3) Soil bearing capacity constraint
reach an optimum solution. These initial solutions are selected as
4) Ultimate shear strength constraints for footings
randomly. Solutions are categorized according to their fitness.
5) Ultimate flexural strength constraints for footings
The fitness of each solution is determined according to its
Cross-section dimensions of piers are not taken as design closeness to the optimum solution. Then, the genetic process
variables. Instead of this, cross-section of pier is selected begins. Reproduction, crossover, and mutation operators are
proportionally with those of application project. Selected cross- applied to the population at each iteration of genetic process.
section is checked according to ultimate flexural strength constraints This process continues until the required convergence is
o
and ultimate shear strength constraints. The cross-section of pier achieved (Goldberg, 1989; Ayvaz and Aydin, 2000; Dalog lu and
is increased proportionally in the case of inadequacy. Ultimate Aydin, 1999).
flexural strength constraints for piers are calculated as: GAs are widely used to solve structural engineering problems.
Several of such studies are realized by Rajeev and Krishnamoorthy
( Pu, Mu )
g29 = ------------------------
-–1 (6) (1998), Jenkins (1992), Rajan (1995), Pezeshk et al. (2000),
ϕ ( Pn, Mn )
Sarma and Adeli (2000), Saka (2007), Kumaran et al. (2010),
where Pu and Mu are the factored axial load and moment at the Degertekin et al. (2008) and Kameshki and Saka (2001).
piers, Pn and Mn are the nominal axial load strength and nominal In this study, a modified hybrid GA is used for the optimum
moment strength of the piers, respectively. Ultimate shear design. Three kinds of crossover and mutation operator are used.
strength constraints for piers can be written as: The major modification in GA is the use of controlled mutation.
The best solutions are reached at less iteration if controlled
Vu
g30 = ----------------------
-–1 (7) mutation is used.
φ ( Vc + Vs )
where Vu is the factored shear force at the pier section, Vc and Vs 3.1 Penalized Objective Function
are the shear strengths of piers provided by the concrete and the GAs are suitable for unconstrained optimization problems.
shear reinforcement, respectively. Therefore, a penalized objective function is used to convert the
Dimensions of spread footing are also not taken as design constrained problem to an unconstrained problem. The penalized
variables. Horizontal dimensions of footings are determined objective function is constituted by adding a penalty function to
according to soil bearing capacity constraint. The thickness of the objective function. So, the penalized objective function can
foundation is determined according to shear strength constraints. be expressed as:
Then, reinforcement of foundation is designed according to
Φ = CT ⋅ [ 1 + K ⋅ P ] (11)
ultimate flexural strength constraints. Thus, soil bearing capacity
constraint can be calculated as: where K is a penalty coefficient that depends on the kind of the
problem. It is important to determine suitable value of K to reach penalized objective function among the solution strings in each
the best solutions in a short time. P is the penalty function, and generation, respectively. A fitness factor is described to decide
can be given as: which solutions are removed from generation. The fitness factor
n for ith solution string is calculated as:
P= ∑ pi (12)
Fc, i = Fi ⁄ Fort (15)
i=1
where n is the total number of constraints, pi is the violation where Fort is the average of the fitness values. Thus, the fittest
factor of the ith constraint and can be calculated as: solutions are placed twice instead of the weakest solutions in the
mating pool by reproduction operator. The solutions in the
pi = gi for gi > 0 , pi = 0 for gi ≤ 0 (13)
mating pool are known as parent solutions (Goldberg, 1989;
where gi is the normalized value of ith constraint. Aydin, 2006).
are selected randomly. The mutation rate is generally determined Table 2. Unit Prices
to be between 0.1% and 1%. Three kinds of the mutation Production Unit Unit price ($)
operator are used in this study. They are standard mutation, PC (uppc) m3 900
controlled mutation, and hybrid mutation (Goldberg, 1989; PS (upps) ton 9500
Aydin, 2006). Detailed information about mutation operator can Cast-in-place concrete for bridges (upbc) m3 220
be found in the studies realized by Aydin and Ayvaz (2010), and Ordinary cast-in-place concrete (upc) m3 80
Aydin (2006). Lean concrete (uplc) m3 80
Reinforcing steel (ups) ton 1500
4. Numerical Example Formwork (upf) m2 35
Excavation (upex) m3 42
Several numerical examples are designed to demonstrate
efficiency of the program. One of these examples is given here.
Selected bridge has been previously designed by Emay variable values. There are 16 possible values for the 2nd, 3rd, and
International Eng. Cons. Co. (Emay, 1998) and constructed in 4th design variables, so that a four-character string must be used
Turkey. Design parameters of the bridge for this example are to demonstrate each of these design variable values. Similarly,
taken as the same with application project and given in Table 1. there are 32 possible values for the 10th design variable, so that a
The unit prices of the productions are taken according to five character string must be used for this design variable.
Turkish unit price criteria and given in Table 2. The site Therefore, each of the solution strings must have 35 characters.
altitudes of the valley are also considered as the same with In other words, 35 number characters (“1” or “0”) must be
those of application project. determined randomly to constitute each initial solution. Design
Design variables values sets for ten design variables are given variable values sets are taken to be the same with those of Aydin
in Table 3. There are 8 possible values for the 1st, 5th, 6th, 7th, 8th, and Ayvaz (2010).
and 9th design variables, as shown in Table 3, so that a three- Several of the optimum design parameters for GA are
character string must be used to demonstrate each of these design determined as follows for this example.
Fig. 4. Cross Section of the Girder: (a) For the Application Project, Fig. 7. Longitudinal Section of the Bridge: (a) For the Application
(b) For this Study (Dimensions in Centimeters) Project, (b) For this Study (Dimensions in Meters)
Table 6. Footing Dimensions for the Application Project and for this Study
Dimension
Solution 1st footing 2nd footing 3rd footing 4th footing 5th footing 6th footing
(m)
B 12.00 12.50 12.50 12.50 12.50 11.00
This study L 15.00 19.00 19.00 18.00 15.00 13.00
H 2.45 2.45 2.45 2.45 2.45 2.45
B 13.00 13.20 13.20 13.00 12.00
Application project L 16.50 20.40 20.40 16.50 15.00
H 2.50 2.50 2.50 2.50 2.50
Fig. 4-Fig. 6, and Fig. 7, respectively. Footing dimensions for the n: Total number of constraints
application project and for the optimum design are given in Table 6. P: Penalty function
pi: Violation factor of the ith constraint
5. Conclusions Pn: Nominal axial load strength of the piers
Pu: Factored axial load at the piers
1. Cost optimization of PC I-girder bridges is realized by using q: Maximum stress for the soil under the spread footing
GA. It is concluded that the program coded in this study can qallow: Allowable stress of the soil
be efficiently used in the optimizations of PC bridges. upc: Unit price of ordinary concrete
2. The optimum number of piers, the optimum number of PC upbc: Unit price of cast-in-place concrete for bridges
girders and optimum shape of PC girder are found out at the upex: Unit price of excavation
same time. Although it is very difficult to solve such a com- upf: Unit price of formwork
plicated problem with usual optimization techniques, it is uplc: Unit price of lean concrete
solved easily by GA. uppc: Unit price of precast concrete
3. The solutions reached by GA are realistic and constructible ups: Unit price of reinforcing steel
without modifying because of using discrete design vari- upps: Unit price of prestressing steel
ables and probable design variable values sets. Vc: Shear strength of piers provided by the concrete
4. Shape of the site is considered, and heights of the piers are Vn,f: Nominal shear strength of the footings
determined according to site altitude. Therefore, optimum Vs: Shear strength of piers provided by the shear
number of piers is determined as more realistic. reinforcement
5. The solution obtained by the proposed GA is up to %12.6 Vu: Factored shear force at the pier section
more economical than the application project. Vu,f,: Factored shear force at the footings
6. In this study, various crossover and mutation types are used. φ: Strength reduction factor
The best solutions are obtained by using the uniform cross- Φ: Penalized objective function
over and the controlled mutation. Φmax: The largest value of the penalized objective function in
7. Constraints can easily be adapted to the problem using the each generation
penalty function and design variable value sets. Φmin: The smallest value of the penalized objective function in
each generation
Notations
References
Cf: Cost of foundations
AASHTO (2002). Standard specifications for highway bridges, American
Cg: Cost of girders
Association of State Highway and Transportation Officials, Washington,
Cp: Cost of piers D.C.
CT: Total cost (Objective function) Aydin, Z. (2006). Optimum design of prestressed concrete bridge
Ctb: Cost of transverse beams girders using genetic algorithm, PhD Thesis, Karadeniz Technical
Fc,i: Fitness factor for the ith solution string University, Turkey.
Fi: Fitness value of the ith solution Aydin, Z. and Ayvaz, Y. (2010). “Optimum topology and shape design
Fort: Average of the fitness values of prestressed concrete bridge girders using a genetic algorithm.”
gi: Normalized value of ith constraint Struct. Multidiscip. Optim., Vol. 41, No. 1, pp. 151-162.
Ayvaz, Y. and Aydin, Z. (2000). “Optimum design of trusses using a
K: Penalty coefficient
genetic algorithm.” Proceedings of the Second International Conference
Mn: Nominal moment strength of the piers on Engineering using Metaphors from Nature, Leuven, Belgium, pp.
Mn,f: Nominal moment strength of the footings 159-168.
Mu: Factored moment at the piers Cohn, M. Z. and Lounis, Z. (1994). “Optimal design of structural
Mu,f: Factored moment at the footings concrete bridge systems.” J. Struct. Eng., Vol. 120, No. 9, pp. 2653-
2674. concrete bridge girder systems.” PCI J., Vol. 38, No. 4, pp. 60-78.
o
Daloglu, A. and Aydin, Z. (1999). “Kafes sistemlerin uygulamaya Lounis, Z., Mirza, M. S., and Cohn, M. Z. (1997). “Segmental and
yönelik optimum tasarimi.” Mühendislik Bilimleri Dergisi, Vol. 5, conventional precast prestressed concrete I-bridge girders.” J.
No. 1, pp. 951-957 (in Turkish). Bridge. Eng., Vol. 2, No. 3, pp. 73-82.
Degertekin, S. O., Saka, M. P., and Hayalioglu, M. S. (2008). “Optimal Pezeshk, S., Camp, C. V., and Clem, D. (2000). “Design of nonlinear
load and resistance factor design of geometrically nonlinear steel framed structures using genetic optimization.” J. Struct. Eng., Vol.
space frames via tabu search and genetic algorithm.” Eng. Struct., 126, No. 3, pp. 382-389.
Vol. 30, No. 1, pp. 197-205. Rajan, S. D. (1995). “Sizing, shapes and topology design optimization
Emay International Eng. Cons. Co. (1998). Uzunçayir baraji varyanti of trusses using genetic algorithm.” J. Struct. Eng., Vol. 121, No. 10,
yolu dinar köprüsü uygulama projesi hesaplari, I·stanbul (in pp. 1480-1486.
Turkish). Rajeev, S. and Krishnamoorthy, C. S. (1998). “Genetic algorithm-based
Goldberg, D. E. (1989). Genetic algorithms in search, optimization and methodology for design of reinforced concrete frames.” Comput.-
machine learning, Addison-Wesley Publishing Company, Inc., New Aided. Civ. Infrastruct. Eng., Vol. 13, No. 1, pp. 63-74.
York, N. Y. Saka, M. P. (2007). “Optimum topological design of geometrically
Holland, J. H. (1975). Adaptation in natural and artificial systems, nonlinear single layer latticed domes using coupled genetic algorithm.”
University of Michigan Press, Ann Arbor, Mich. Comput. Struct., Vol. 85, Nos. 21-22, pp. 1635-1646.
Jenkins, W. M. (1992). “Plane frame optimum design environment Sarma, K. C. and Adeli, H. (2000). “Fuzzy discrete multicriteria cost
based on genetic algorithm.” J. Struct. Eng., Vol. 118, No. 11, pp. optimization of steel structures.” J. Struct. Eng., Vol. 126, No. 11,
3103-3112. pp. 1339-1347.
Jones, H. L. (1985). “Minimum cost prestressed concrete beam design.” Sirca, G. F. and Adeli, H. (2005). “Cost optimization of prestressed
J. Struct. Eng., Vol. 111, No. 11, pp. 2464-2478. concrete bridges.” J. Struct. Eng., Vol. 131, No. 3, pp. 380-388.
Kameshki, E. S. and Saka, M. P. (2001). “Optimal design of nonlinear Totres, G. G. B., Brotchie, J. F., and Cornell, C. A. (1966). “A program
steel frames with semi-rigid connections using a genetic algorithm.” for the optimum design of prestressed concrete highway bridges.”
Comput. Struct., Vol. 79, No. 17, pp. 1593-1604. PCI J., Vol. 11, No. 3, pp. 63-71.
Kumaran, S. S., Muthukumaran, S., and Vinodh, S. (2010). “Optimization Yu, C. H., Das Gupta, N. C., and Paul, H. (1986). “Optimization of
of friction welding of tube to tube plate using an external tool.” prestressed concrete bridge girders.” Eng. Optim., Vol. 10, No. 1, pp.
Struct. Multidiscip. Optim., Vol. 42, No. 3, pp. 449-457. 13-24.
Lounis, Z. and Cohn, M.Z. (1993). “Optimization of precast prestressed