# Protective and Blast Resistive Design of Posttensioned Box Girders Using Computational Geometry

Authors: Majid Aleyaasin
Journal: Advances in Civil Engineering
 (2018)

---

## Abstract

The optimal tendon profile and its associated duct geometry for posttensioned box girders are investigated. A computational algorithm has been developed to determine an ideal shape for the tendon and duct. The algorithm is based on a diagonal cross point in the Magnel quadrilateral and uses computational geometry instead of graphical drawing. Thereafter, an ideal parabolic and linear tendon profiles can be calculated from which the covering duct could be shaped. To check the optimality of the results, an automatic examination of Magnel diagrams in various cross sections is incorporated in the algorithm. This enables a unique prestress level to be selected that suits all sections. Then, the mideccentricity of the two crossing points with a common prestress line is chosen as a design eccentricity in each cross section. The optimal duct shape is determined based on such automatic inspections. In a numerical example, the linear, parabolic, and optimal duct geometries are compared and drawn. It is concluded that both linear and parabolic duct shapes can be very close to an optimal shape. In a numerical example, an optimal, box girder with linear open access-type tendons is designed that can withstand extra blast load when explosion occurs.


---

## Body

## 1. Introduction

Steel is used in concrete bridges as a reinforcement material. Therefore, the protection against corrosion for increasing lifetime of the bridge is an important issue and remains an active field of research [1]. Advanced methods for protection needs accessibility inside the box girder as described in [2]. A suitable open access type to the prestressed tendons is open linear duct type as shown in Figure 1 [3]. In this scheme, tendons are accessible easily and can be monitored continuously.Figure 1
(a) Externally posttensioned: inside of box girder cross section. (b) Externally posttensioned beam: girder outside web.
(a)
(b)It is obvious that the accessibility to the steel tendon as shown in Figure1 is superior to other types of design. This will be useful particularly when inspection for corrosion and crack protection is required. The corrosion protection in precast bridges is a formidable task as explained in [4]. Therefore, open access posttensioned segments similar to Figure 1 are corrosion and also crack protective.Recently, protection of the bridges against explosions due to, for example, car bombs is studied seriously, and one of the damages is shown in Figure2. The permanent deformation is caused by immediate microcracks just after explosion [5]. The author believes that the open access design shown in Figure 1 can also be blast protective because it can be optimised to take the extra blast load as it is shown in a numerical example in this paper. Moreover, the possible cracks after explosions can be inspected much easier.Figure 2
Permanent deformation of box girder after explosion.Regarding the extra blast load, the bomb power or explosive amount (equivalent TNT mass) and also the distance between the explosion centre and the bridge are required. Then, we need to find the scaled distance from the following equation:(1)Z=RMTNT1/3.Then, the resulting overpressure (in kPa) due to the explosion [6] is(2)p=912.88Z−1.46635.Then, the extra blast load isWex=pb (b is the width of the bridge).This paper provides a new method that is suitable for the optimal design of an accessible linear-type, posttensioned tendon and its duct. The method is based on the traditional Magnel diagram [7], but its novelty is the combination with computational geometry, by which an algebraic computerised algorithm can be developed.There are two types of prestressed concrete members. In pretensioned type, the strands are pretensioned before casting, and then it is essential that the prestress level in various cross sections should be different. In order to achieve this objective, debonding technology is adopted by which some strands are partially debonded from concrete [7]. This article is not related to pretensioned beams since the majority of box girders are built via posttensioned members.However, in posttensioned beams and members, the prestress level in all cross sections is similar. The tensioning is planned after casting; therefore, the tendons are not straight and can be covered in a curved duct [8]. There are various methods for determination of the ideal tendon and duct profile; for example, the criteria in [8] are different from the one in [9]. In order to obtain an accurate profile, “design curves” are introduced in [10] by which in any cross section the suitable eccentricity is determined by the design curve and finally the profile could be built.Since an ideal tendon profile determination is a design problem, it should be noted that the designers have been using a graphical design method which was developed by Magnel [11] in late forties. Since then, the Magnel design diagram has been revisited in [12], and the edges of the Magnel quadrilateral are specified as criteria for choosing an appropriate eccentricity and prestress level. Recently, Calladine [13] has used the mapping technique and converted the Magnel quadrilateral to a parallelogram by which the design procedure becomes easier. Moreover, Stanton [14] introduced a simple sign convention in order to draw the Magnel diagram by one equation only.In this article, it is shown that the ideal tendon profile in the posttensioned beam could be determined by automated computer inspection of the Magnel diagrams, in the various cross sections without drawing any quadrilateral on paper. By using computational geometry, first the coordinates of the edges of the Magnel quadrilateral is calculated. Then, the coordinates of the crossing points of the two diagonals of the quadrilateral in the cross section in which maximum service bending moment occurs is determined. Similar procedure is repeated for two more cross sections, and thereby a parabolic or linear curve could be computed, knowing three eccentricities across the beam span. Further investigation is carried out to show the accuracy of this type of design.In a numerical example for a long span box girder bridge, by examining the Magnel quadrilateral in many more cross sections and considering a prestress level that suits all the cross sections, the corresponding eccentricities are computed and a smooth curve is fitted to all the points. When comparing this optimal duct shape with the parabolic and linear ones, it does not show any significant difference.The whole drawing process and inspection and selection are algorithmic and can be done via one computer program. It is shown that we can determine the ideal parabolic and linear duct shapes and also check their closeness to optimal duct shape, all in one algorithm. It is concluded that both parabolic and linear duct geometries introduced in this paper are simple and optimal. Therefore, we can choose an accessible and easily protectable linear tendons and duct shape which is also an optimal design. The method is based on automatic inspection of the Magnel diagram and particularly can be useful for checking the Magnel polygon in the cross section where the deviator is located and sharp changes (discontinuities) in prestress occur.


## 2. Magnel Design Diagram

The top fibres of a beam in transfer (before service) are in tension since the beam planned to be in hogging, and if the allowable tensile transfer stress isftt≤0, then the stress constraint would be(3)σ=MTZT+PA+PeZT≥ftt.In (3), σ>ftt while σ<ftt because in the transfer state, the top fibres are in tension (− sign by convention) and then ftt≤0 also Zt>0, P>0, Mt>0, and e<0. Manipulating (3) leads to the following inequality:(4)1P≥efttZt−Mt+ZtAZtftt−AMt.The bottom fibres of a beam in transfer, when it is in hogging, are in compression, and if the allowable compressive transfer stress isftc≥0, then the stress constraint would be(5)σ=MtZb+PA+PeZb≤ftc.In (5), σ<ftc while σ<ftc because in the transfer state, bottom fibres are in compression (+ sign by convention) and then ftc>0 also Zb<0, P>0, Mt>0, and e<0. Manipulating (5) leads to the following inequality:(6)1P≥eftcZb−Mt+ZtAZbftc−AMt.The top fibres of a beam in service obviously are in compression because the beam is in sagging, and if the allowable compressive service stress isfsc≥0, then the stress constraint would be(7)σ=MsZt+αPA+αPeZt≤fsc.In (7), σ<fsc while σ<fsc because in the service state, top fibres are in compression (+ sign by convention); so fsc>0 also Zt>0, P>0, Ms>0, and e<0. Parameter 0≤α≤1 represents the minimum conserved prestress factor in service since during service prestress level would be lost [8] and 1−α expresses this maximum loss. Now, by manipulating (7), the following inequality would be obtained:(8)1P≥αefscZt−Ms+αZtAZtfsc−AMs.The bottom fibres of a beam in service (sagging) are in tension, and if the allowable tensile service stress isfst≤0, then the stress constraint would be(9)σ=MsZb+αPA+αPeZb≤fst.In (9), 0≤α≤1 appears again because of service condition. σ>fst while σ<fst because in the service state, bottom fibres are in tension (− sign by convention), so fst≤0 also Zb<0, P>0, Ms>0, and e<0. By manipulating (9), the following inequality would be obtained:(10)1P≤αefstZb−Ms+αZbAZbfst−AMs.The four stress constraints (3), (5), (7), and (9) lead to the inequalities (4), (6), (8), and (10) in a coordinate system in which e is the horizontal axis and P−1 is the vertical axis. The borders of these four inequalities are designated by four lines, and when they cross each other, they form a quadrilateral in which the design parameters e and P−1 are chosen from inside of this quadrilateral. Overall design diagram was first introduced by Magnel [11] and has been used by designers ever since [7, 8]. Obviously, in various types of prestress losses and for this reason for design purposes, the minimum conservation (maximum loss) should be considered.In Section3, we will show that, by symbolic computation of the four edges in the Magnel quadrilateral, we can seek inside the region to find an appropriate design point for each cross section. First, parabolic and linear tendon profiles will be derived by using three cross sections. Thereafter, an algorithm for an ideal tendon profile will be explained by using several cross sections. In a numerical example, it will be shown that parabolic and linear profiles are very close to the ideal profile.


## 3. New Algorithm for Computation of Tendon Profile

The borderlines of the constraints forftt in (4) and ftc in (6) cross each other in the edge designated by coordinates ett−tc and P−1tt−tc as follows:(11)ett−tc=Ztftt−MtZb−Zbftc−MtZtZbftc−MtA−Ztftt−MtA,P−1tt−tc=ett−tcA−ZtZtftt−MtA.The borderlines of the constraints forfst in (10) and fsc in (8) cross each other in the edge designated by coordinates est−sc and P−1st−sc as follows:(12)est−sc=Zbfst−MsZt−Ztfsc−MsZbZtfsc−MsA−Zbfst−MsA,P−1st−sc=est−scA+ZbαZbfst−MsA.The borderlines of the constraints forfst in (10) and ftt in (4) cross each other in the edge designated by coordinates est−tt and P−1st−tt as follows:(13)est−tt=Zbfst−MsZt−Ztftt−MtαZbZtftt−MtαA−Zbfst−MsA,P−1st−tt=est−ttA−ZtZtftt−MtA.The borderlines of the constraints forfsc in (8) and ftc in (6) cross each other in the edge designated by coordinates esc−tc and P−1sc−tc as follows:(14)est−tc=Ztfsc−MsZb−Zbftc−MtαZtZbftc−MtαA−Ztfsc−MsA,P−1sc−tc=esc−tcA+ZtZbftc−MtA.The edge with the coordinatesest−tt and P−1st−tt is opposite to the edge with the coordinates esc−tc and P−1sc−tc. Connecting these two edges provides a diagonal with the following slope:(15)msc−tcst−tt=P−1st−tt−P−1sc−tcest−tt−esc−tc,and the corresponding line can be expressed by the following equation:(16)P−1−P−1st−tt=e−est−ttmsc−tcst−tt.The edge with the coordinatesett−tc and P−1tt−tc is opposite to the edge with the coordinates est−sc and P−1st−sc. Connecting these two edges provides a diagonal with the following slope:(17)mst−sctt−tc=P−1tt−tc−P−1st−scett−tc−est−sc,and corresponding line can be expressed by the following equation:(18)P−1−P−1tt−tc=e−ett−tcmst−sctt−tc.Intersection of the lines (16) and (18) represents the centre of the quadrilateral which possesses the eccentricity eDM and can be calculated by(19)eDM=P−1tt−tc−P−1st−tt+mst−tcst−ttest−tt−mst−sctt−tcett−tcmst−tcst−tt−mst−sctt−tc.Similarly, we can calculateeDM in two other cross sections to obtain parabolic and linear curves. For example, in a simply supported beam with span L, the maximum bending moment occurs in middle span 0.5 L, the other two cross sections could be at the supporting ends, and their corresponding eccentricity eDC can be calculated from (19) by substituting Mt=0 and Ms=0, that is,(20)eDC=eDMMs=Mt=0.Since atx=0.5L, the eccentricity is e0.5L=eDM and in the supports eL=e0=eDC, and then an ideal parabolic tendon profile could be computed by(21)ex=4eDM−eDCxL−x2L2+eDC.And for linear tendon profile similar to Figure 1, we have(22)ex=2eDM−eDCxL+eDC,0≤x≤L2,(23)ex=2eDM−eDC1−xL+eDC,L2≤x≤L.Another ideal tendon profile can also be achieved, if the Magnel diagram would be drawn in many cross sectionsn>>3, across the span. This is a formidable task for a designer, but can be programmed easily by a computer. Then, for every Magnel quadrilateral, we have a pair of diagonals with equations(24)P−1−Pi−1st−tt=e−eist−ttmisc−tcst−tt,i=1,2,…,n,(25)P−1−Pi−1tt−tc=e−eitt−tcmist−sctt−tc,i=1,2,…,n.By simple visual inspection of all diagrams in one figure, it is very easy to check which prestress level suits all the quadrilaterals in the figure. This inspection is also possible without any drawing and can be done by intelligent computer searching. When a design prestress level was selected and namedP−1D, by its substitution into (24) and (25), two candidates for the eccentricity in each cross section can be found via the following formulas:(26)eiC1=P−1D−Pi−1st−ttmisc−tcst−tt+eist−tt,i=1,2,…,n,(27)eiC2=P−1D−Pi−1tt−tcmist−sctt−tc+eitt−tc,i=1,2,…,n.To be on the conservative side, the average of the above could be considered as the design eccentricity for each cross section, that is,(28)eiD=0.5eiC1+eiC2,i=1,2,…,n.


## 4. Numerical Example

The numerical example is taken from [12]. The prestressed beam, simply supported at the ends with the cross section A=0.3045m2, has the span L = 18.3 m, the width b = 1.22 m, the height h = 0.61 m, and the section modulus Zt=0.0352m3 and Zb=−0.03869m3; the limit stresses are ftt=−1.34MPa, ftc=17.38MPa, fsc=18.6MPa, and fst=−3.21MPa; the dead weight is Wt=8.76kN/m; and 500 kg of TNT exploded 200 m away from the bridge while traffic load Wtr=5.5kN/m in service condition should be considered. The prestress loss is expressed by the parameter α=0.8. We need to find the optimal duct shape and the eccentricities as shown in Figure 1 and displayed in [3].The scaled distance of this explosion will be considered as a part of service load:(29)Z=RMTNT1/3=2005001/3=25.1984m/kg1/3.The resulting overpressure (in kPa) due to the explosion is(30)p=912.88Z−1.46635=912.88×25.1984−1.46635=8.045kPa.This produces overpressure load as a result of explosion:(31)Wex=pb=8.045×1.22=9.815kN/m,Wa=Wex+Wtr=9.815+5.5=15.32kN/m.The maximum transfer bending moment occurring at a midspan is(32)Mt=WtL28=8.76×18.328=366.1kN⋅m.The maximum transfer bending moment occurring at a midspan is(33)Mt=Wt+WaL28=8.76+15.32×18.328=1006.8kN⋅m.The Magnel diagram in the maximum bending moment cross section is drawn in Figure3; the quadrilateral occupies a small region in the coordinate system, and for the designer, it is difficult to spot the centre, and the quadrilateral has been magnified in the figure.Figure 3
Magnel design diagram in the midspan cross section.In Figure4, the design diagram is drawn at the ends, in which both the transfer moment Mt=0 and service moment Ms=0, we can see that the relevant quadrilateral occupies much larger region in the coordinate system plane. Although the quadrilateral losses its convexity, the centre of it is still an appropriate design point. One can see that the design point will be outside the quadrilateral, but its eccentricity coordinate is in the acceptable range. The designer finally will choose a prestress level which suits all cross sections.Figure 4
Magnel design diagram in the supporting ends.Based on (19) and (20), parabolic and linear tendons are designed, and the position of the strand, relative to the neutral axis in each cross section, is shown in Figure 5. The neutral axis represents zero eccentricity and is shown by dash ∗ line; sign × represents top and bottom fibres. The linear strand position is shown by the solid line, and the point line shows the parabolic strand.Figure 5
Position of the tendon, neutral axis, and top and bottom fibres across the span.In order to see if these types of tendon profiles are accurate enough, we will draw the Magnel diagram in nine cross sections each withL/8 apart. Because of symmetry, we need only five Magnel diagrams, which all are drawn in Figure 6. Although we have plotted each diagram in separate figure, it is difficult to compare all quadrilaterals. However, it is easy to find out that a prestress level equal to the average of prestresses at edges of the quadrilateral of the midspan is an appropriate criterion. The horizontal line corresponding to this prestress level intersects all Magnel quadrilaterals.Figure 6
Magnel diagrams for the cross sections: (a) (L/2), (b) (L/4, 3L/4), (c) (L/8, 7L/8), (d) (3L/8, 5L/8), and (e) (ends 0, L).
(a)
(b)
(c)
(d)
(e)This enabled (26)–(28) to be used, and in each cross section, eiD could be computed. Using eiD, which resulted in nine points, by which a symmetric curve is produced and the tendon profile is displayed in Figure 7 by the point line, one can see that the simple linear profile in Figures 1(a) and 5 and the optimal profile in Figure 7 are very close, and in fact the parabolic and linear shape tendons can be ideal ones, if we use the expressions in (21)–(23) in this paper.Figure 7
Ideal tendon profile based on the Magnel diagram of 9 cross sections eachL/8 apart.


## 5. Conclusions

In this article, the ideal tendon profile in the posttensioned beams is investigated by multiple automatic examinations of the Magnel design diagrams. An algorithm was introduced by which one can calculate ideal parabolic and linear tendon profiles. It is shown that ideal parabolic and linear shape tendons and ducts, calculated by the formulas in this paper, are very close to the ideal and optimal tendon profile. It is recommended that both linear and parabolic ducts can be used by designers because they do not have any disadvantage. Moreover, the design is based on simple formulas without engaging with graphical work. It can also be implemented in design of continuous girders, assuming that the moments in intermediate pier supports are found by static calculations. Moreover, the blast load as a result of car bomb and so on can be considered and added to the service load. In this case, the relevant eccentricity and prestress level can be found as shown in the numerical example. This avoids the permanent deformation (Figure2) as a result of explosion near the bridge.As far as that author is aware, there is not any design method which is particularly suitable for design of accessible linear-type tendons and ducts in posttensioned beams. Instead, designers use traditional methods and increase the factor of safety to compensate lack of accuracy in modelling. Therefore, the method in this paper can be used for future design applications in which bridge design should be optimised to take an extra blast load as a result of explosion due to car bomb and so on. Moreover, in future by installing sensors in the open access box girder, compression conservation can be monitored and the validity of formulas in this paper can also be checked.



---

*Source: 4932987-2018-02-22.xml*