Bridge Engineering by Christian Menn
Bridge Engineering by Christian Menn
Bridge Engineering by Christian Menn
A
B
Life Sciences
Physics
Mathematical and
Physical Sciences
Behavioural and
Social Sciences
Applied Sciences
Computer and
Systems Sciences
Ecological Sciences
Springer-Verlag
Berlin/Heidelberg/New York
E
F
No. 74
Cetin
Y,lmaz
,
Associate Professor
Civil Engineering Department
Middle East Technical University
Ankara, Turkey
s.
Tanvir Wasti
Professor
Civil Engineering Department
Middle East Technical University
Ankara, Turkey
1984
ISBN-13: 978-94-009-6124-1
DO I: 10.1007/978-94-009-6122-7
e-ISBN-13: 978-94-009-6122-7
Distributors for the United States and Canada: K1uwer Boston, Inc., 190 Old Derby
Street, Hingham, MA 02043, USA
Distributors for all other countries: Kluwer Academic Publishers Group, Distribution
Center, P.O. Box 322, 3300 AH Dordrecht, The Netherlands
v
PREFACE
The Proceedings of the NATO Advanced Study Institute on
Analysis and Design of Bridges held at ~eme, lzmir, Turkey from
28 June 1982 to 9 July 1982 are contained in the present volume.
The Advanced Study Institute was attended by 37 lecturers and
participants from 10 different countries. The Organizing Committee
consisted of Professors P. Gtilkan, A. C. Scordelis, S.T. Wasti and
9. Yl.lmaz.
The guidelines set by NATO for the Advanced Study Institute
require it to serve not only as an efficient forum for the
dissemination of available advanced knowledge to a selected group
of qualified people but also as a platform for the exploration of
future research possibilities in the scientific or engineering
areas concerned. The main topics covered by the present Advanced
Study Institute were the mathematical modelling of bridges for
better analysis and the scientific assessment of bridge behaviour
for the introduction of improved design procedures. It has been
our observation that as a result of the range and depth of the
lectures presented and the many informal discussions that took place,
ideas became fissile, the stimulus never flagged and many gaps in
the engineering knowledge of the participants were "bridged".
Here we particularly wish to mention that valuable informal presentations of research work were made during the course of the Institute
by Drs. Friedrich, Karaesmen, Lamas and Parker. Contributing
tangibly to the success of the Advanced Study Institute were the
several opportunities for social get-togethers in which the ladies
a Iso took part.
The editing of the proceedings has entailed a moderate amount
of work in the case of a few manuscripts and only very slight
interference with some others. In general the proceedings faithfully
represent the lectures as transmuted from the spoken to the written
word by the lecturers themselves. Unfortunately Professor B.T. Yenis
contribution was not received in time for publication.
Inspite of the inevitable overlapping, it has been found
suitable and convenient to group the lectures under the following
five main headings :
I.
II.
III.
IV.
V.
Design Considerations
Analysis Methods for Steel Bridges
Concrete Slab and Box Girder Bridges
Dynamics and Vibration
Earthquake and Wind Effects
VI
t;:ETtN YILMAZ
S. TANVIR WAS TI
Edi tors
VII
TABLE OF CONTENTS
Preface
I. DESIGN CONSIDERATIOr.rS
Christian Menn
Prestressing of Curved Bridges
Fritz Leonhardt
Certain Aspects of Bridge Design
23
49
57
H. Roy Evans
Simplified Methods for the Analysis and Design of
Bridges of Cellular Cross-section
95
119
191
Celal Kostem
Lateral Live Load Distribution
Highway Bridges
213
~n
Prestressed Concrete
225
VIII
IV. DYNAMICS AND VIdRATION
Nuri Akka and ~etin YLlrnaz
Fluid-Structure Interaction in Analysis of Bridge piers.
Part I : Theoretical Preliminaries
243
259
~etin
Etienne Zeller
Vibrations and Dynamic Behaviour of Actual Bridges
285
343
Polat Giilkan
Analysis and Design of Bridges for Earthquake
Effects
375
Etienne Zeller
Wind Loads on Bridge Structures
415
List of Participants
449
I. DESIGN CONSIDERATIONS
Christian Menn
Professor fur Baustatik und Konstruktion
Eidgenossische Technische Hochschule
Zurich
SUMMARY
1.
INTRODUCTION
FIG. 1 J
FIG. 2:
~Q
Jr
J
1/2
1/2
f!
T = Q. ~
GK -
O-collapse!
Fig 30
FIG. 3A:
EQUILIBRIUM TORSION
~J~-l/-2-;Q-l-/2-t/.
Fig 3b
FIG. 3B:
CONSTRAINT TORSION
3 GK'I
T=Q 16.(3.EI+GK)
GK -
0 -
T -
K;7
, _-,
....
./
Moments occording to
th~ory
of elasticity
Continous girder
Fig 4
FIG. 4:
MOMENT REDISTRIBUTION
redistributed moments;
Theory of plasticity
all supports
Fig 5
FIG 5:
<p
ARRANGEMENT OF BEARINGS
'Po
=wLl.T+t. cs
Fig 6
FIG 6:
(P)
10
curvature could conceivably be designed without any Joints.
design, although certainly feasible, is seldom used.
2.
Such a
+ qds
0
0
0
Thus,
d2 M
ds 2
dT
Os
1dI'
-(q - FOs)
M
-(F + eq + t)
=-~
curved
continuous
girder,
the
bending moments
are
Fig 7
FIG. 7:
simplified
method
M*
A -125.00
T*
0
-12550 -0.023
- 0.278
2.863
1.042
62.461
1.047
62.50
62.50 -2.344
0
2865
exact
method
-4.167
62.549 -2.339
0
-4.164
Fig8
FIG. 8:
12
ds2
dT
as
-q
_(!i
+ e'q +
r
t)
PRESTRESSING
13
d+dv h+dT
----.. ----....T ~
M
o . d., =ro .ds =r:ti
. ds
dT =- M d.,. = dv
M
r
ds
=.....2L= ___M_ds
2bh
2bhr
rh
2rh
Fig 9
FIG. 9:
t
C = Shear centre
S = Centre of gravity
y
Fig 10
p
FIG. 10:
14
v cy
da
da
-p x (ay .~
~.)
ux - a Z CI'X"
pxa y
daz
da
-p x '[(a-c
(a-c
...::::::yj
y y )'=-:-ux
Z
z ) ax:
AS can be seen, torsion in the concrete cross section may be
caused by tendons in the webs as well as by tendons in the flanges.
In the following, the required arrangement of the prestressing
tendons in order to ensure a specific distribution of the torsion
moments (e.g. Tcp = - Tcg) will be presented, presuming the girder
cross section is a symmetrical one cell hollow box.
3.1
TCp '1 0
_a l
y
a yr
da
-b/2
Vcy
Vcz
Ncy
Mcz
-..
Mcz
ail
-..
Vcy
a zr
_a l
z
Mcy
da~
Vcz
dar
z
CI'X"
= CI'X"
da~
Px
- 2 ' b CI'X"
mtp
= - d:~p = ~
d2 a l
1 b
qpz
~n
15
The distribution of torsion moments in a simply supported
curved girder under dead load is shown in Fig. 11. The equalizing
prestress torsion Tcp = - Tcg requires the following geometry for the
tendons in the web on the outer side of the curve:
da;
2Tc
--ax
Pxb
azl
2.
fTc.dx
and the reverse geometry in the web on the inner side of the curve
(Fig. 12).
This tendon geometry for resistance of torsion may always be
superimposed on the tendon geometry designed for the resistance of
bending. Therefore, in simply supported girders, equalization of
torsion and resistance to dead load bending may be provided with the
same tendons (Fig. 12).
The distribution of torsion moments in a continuous curved
girder under dead load is shown in Fig. 13. This distribution is
clearly different from that in the simply supported girder under the
same load (Fig. 11). Redundant moments due to prestressing designed
to resist bending influence the torsional moments,
generally
resulting in the torsional moment distribution shown in Fig. 14.
In this case the equalizing prestress torsion TCl' = - Tcg
requires the following geometry for the tendons in the web on the
outer side of the curve:
da; _
--ax -
~_~
3.2
Tcp
t-
Mcy
pi
Px
_ps.a s
x z
pi.ai
x Z
....
as
z
-(h-a zi )
Vcy
Vcz
Mcy
Mcz
p~, pi
x
16
Fig 11
FIG 11:
a:
b.
!=
=-
~
c:
~.
~
a:
b:
Fig 12
FIG.
12:
-~~
=:24
17
Fig 13
FIG 13:
Fig 14
FIG. 14:
a;
Fig 15
FIG. 15:
a:
b:
18
da z
---ax
Vcz
Mcz
Vcy = 0
s
-p~. y
pi.ai
x y
da s
-pstfx
x
x
P~tfx
-+
da i
-+
a;
as
y'
da s
da i
tfx, tfx
pS.h.~
x
x
Tcp
. da s
-p .al."tfx
x z
x
dT
~p
px
a~.~y
2 s
l.. d a v i s
pxaz~
azqpy
mtp
3.3
~~
PY
~
da
dx
d2 a
__
z
dx2
cause
the
- qpy.(az-c z }
qpz .(ay - c y )
The deviation forces q; of the tendons are in equilibrium
with the deviation forces
the normal stresses in the concrete.
y
cf6 o'l
qpy
qcy
19
qpy +
~~y
dV c
-?xqpz +
dV cz
dT(V C:i)
ax
dA c
Ac
----a:x-
dV cz
J dT(Vdxcz ) .dA c
----a:x-
Ac
dTc(qE:i)
mtp(qpy) +
ax
dTc(qEZ)
ax
mtp(qpz) +
dT(T )
c
---ax :; f --axdTc
dA c
20
Fig 16
FIG. 16: DEVIATION FORCES qpy' qpz
foTo(-Vcyl
y-
Fig 17
FIG. 17:
21
qpy in C
I~
. 11
Il,~=~ll
dVCY )
dx
S (qPY)
S --
qpz in C
I~
It
11
dVCZ)
S (-
S (qpz)
dx
m,p (qpy)
It
S (m,p)
(~~
--
~I
If" "JI
~I
~I//=/'It
-----.-
m,p (qpz)
If
S (m,p)
Fig 18
It
S(
~~
22
S(mtp)
m,p(qpz)
=-
m,
It
Fig 19
FIG. 19:
S(m,p)
23
PARTIAL
PRESTRESSIL~G
24
25
~ ductile
I~red
by AT
--I ~ deformation
eAT
Fig. I
f.i. curvature
failure
by 1.5 AT
26
safety, if the structural details are well designed.
1.1
max ~L+LL
This definition has the advantage that the engineer can directly
see if the degree which he intends to choose will, for instance,
cover the full dead load moment, ~L' or will cover dead load plus
30% of live load or will cover only 80% of dead load.
Prof. H. Bachmannof Zurich has published charts giving the
necessary amounts of prestressing and reinforcing steel for
rectangular cross sections depending on ttle degree of pres tress ing
(Fig. 3). He further gives the increase of steel stress in the
prestressing steel and the steel stress in reinforcing bars due to
full service load, again depending on the degree of prestressing
(Fig. 4). For a prestressed concrete railroad bridge with I-beams
such charts were calculated in my office (Fig. 5 and 6). In these
27
300
~= res M
Ph
=-0.6}
=-0.4 prestressed with P
=-0.2
0.2%
Fig. 2
0.5%
1.0%
28
0.1
0.2
0.3
80
120
150
The charts of the steel stresses show that the stress amplitudes
remain far below the fatigue limits if for rectangular cross sections
the degree of prestressing is chosen above 0.4 and for I-beams above
0.5 and if the frequently repeated portion of the live load remains
below about 7070 of the dead load. This means that choosing partial
prestressing in most practical cases does not cause any fatigue
problems if the degree of prestressing and the ratio LL:DL are
within these limits. Therefore, even for railroad bridges partial
prestressing can be chosen.
1.3
29
<~<'
~~!;J
Ap
A. = S 460 '550
A.
/'
/'
=S 1500 '1700
.,/
/'
3Ap
A
min. ~ =-'-~ 0.5'.
r Ac ,.'
./
10
pr.str.ssing
d.gr
~~~~--0~.4-----0~.6-----0~.8----~1D~x
Fig. 3 and 4
--
-..""
100
"
dUp
..........
"
",
.................
L--------L~----~L-------~=_----~~~~--~--~~~~.x
0.2
0.4
prestressing
0.6
degree
0.8
1.0
30
Ap = S 1500 /1700
As
= S 450 /500
200
,~g
It>
(")
0.90
prestressing degree
Fig. 5 and 6
.l
200
100
~----~~------~------~--------L------~--~--~~X
0.6
prestressing
degree
0.8
to
31
2.
2.1
Fig. 7
32
s Eeff
y2
1
1/,2
(1)
c E0
l2a'
where :
A
s
I/,c
s c
horizontal span of cable
33
However, is is easy to reduce the sag effect for longer cables by
introducing so-called stiffening ropes as shown in Fig. 9. In this
way the high stiffness of the cable-stayed system can be maintained
for very long spans, up to main spans of 1500 m for railroad and
highway bridges.
180
140
60
O+---~----~--~--~~lc
100
200
300
Fig. 8
Effective modulus of
elasticity gives influence
of the sag of the cable
on its stiffne s s
400
T=Ac G'
stiffening ropes
Fig. 9
34
2.2
35
The ratio between side span 1 and main span has influence
on the stress changes mainly of the back stay cables, which hold
the tower head back to the anchor pier. Live load in the main span
increases these stresses, live load in the side span decreases them,
and in long side spans with 1/ > 0.4, the cables could become
more or less slack. The back stay cables get the largest stress
amplitudes of all cables and these amplitudes must be kept safely
below the fatigue strength of the cables for that part of the live
load which may occur a million times (for instance at highway
bridges 40% of maximal live load).
The ratio 1/ has further influence on the amount of vertical
anchoring forces at the anchor pier. This anchor force decreases
with increasing 1 / . Towers on both sides of the main span are
assumed here. If there is only one tower on one side (unsymmetrical
case) then the main span acts similar as if the span length is
about l.8.
36
Vl
fan shape
harp shape
1i~:~ ~:~=th
J.--/--J.
.:;-3
c
c
a;
Fig. 13
Quantity of cable steel as
a functIOn of relative height
iii
of tuwers
:l
CT21----\l---'<--,--~--_.--
QJ
L.
2
u
.E
37
Longitudinal bending stiffness of towers, which is characteris tic for Morand i' s early stayed beam bridges, ge t ve ry large
moments, which must be carried by the foundation and can easily
double their cost. In order to avoid this, some towers of cablestayed bridges have been built with foot hinges so that the foundations are concentrically loaded (Mannheim and Oberkasse 1 bridges
across the River Rhine) (4, 6).
0)
C3[~-1
-
__ ...J
section 1-1
cl
b)r -, I
i
low
Flg. 14
Two cable planes along
the edges of the deck
with different solutions
for the towers
2.5
38
39
a}
b)
'''CO M"" co
c)
~2t"
:I~
\
~~
ii?~~
Fig. 15 a
<..
Fig. 15 b c
> 12
15 m
m
40
side span
appnooch span
Fig. Hi
2.7
'
anchored
41
Fig.17
Anchorage of HiAm-cables
in the tower and in the concrete edge girder of the
deck
42
2.9
The towers
43
44
~l~
..J...--.330 m --,j.
[4i)
kN/m
6,0
L35m 1 98.3m
10000n"g
train
at
V: 108.8 km/h
system damping
(mm) 0 102030405060
t(s)
-2~
-400
-6:'"
-80.,
Fig. 18
Oscillation diagram
cable-stayed bridge
of a
45
REFERENCES
1.
2.
3.
4.
46
5.
6.
Volke, E. ; Rademacher, C.-H. : NordbrUcke Mannheim Ludwigshafen (Kurt-S chumache r-BrUcke) De r S tah lb au 1973, Heft
4, p. 97, Heft 5, p. l38, Heft 6, p. 161.
7.
8.
9.
10.
11.
12.
13.
14.
15.
49
INTRODUCTION
Highway bridges may be designed in accordance with the allowable stress, load factor, or load and resistance factor design
approach. When these bridges are loaded with the regular design
vehicles the response of the structure will be linearly elastic, or
very close to it. The collapse load levels of these bridges can be
computed by anyone of the accepted methods available in the literature. However, the overloading of the bridge may propel its structural response to a level that is above the proportional limit but
far below the collapse load. In the bridge engineering literature
there exists a major gap in the prediction of the structural response of the bridge at post-proportional limit, and also the determination of the extent of damage to the superstructure. This paper
provides an overview of the inelastic response of highway bridges
composed of steel girders, rolled, built-up or plate girders, and
reinforced concrete deck slab.
The various research programs that are summarized herein have
all employed the finite element displacement method, Derivations
and the development of the mathematical models have been extremely
laborious. Presentation of any formulation would not have increased the engineering contributions of this paper. It is recommended that readers interested in the mathematical derivations of
the employed models, and the verification of these models, should
refer to the references cited. The paper is aimed at bridge
engineers.
1.1
50
The "routine" methods of analysis of beam-slab highway bridgesthe structures which inherently have a very high degree of structural indeterminancy-employ the oversimplified "distribution factor"
approach (9,10). The studies have indicated that even for the
elastic regime of the structure 'this approach is far from being
satisfactory (2,7). Thus, a more detailed scheme needs to be
employed which will account for the structural indeterminancy, as
well as the gradual initiation and spread of material nonlinearities.
The finite element modeling with plate bending elements, with membrane stiffness properties, and the eccentric beam elements would
provide a realistic approach to the elastic simulation of the superstructure. To account for the material nonlinearities, it is essential that both the beam and the plate bending elements be subdivided
into layers, through the depth of the structural member (Fig. 1).
Throughout the analysis with the gradually increased load levels,
using the tangent stiffness formulation, stress and strain, as well
as the damage, history of each layer of each element need to be
monitored. A complicated analysis scheme as such requires the
development of a detailed software package. Program BOVAS (Bridge
Overload Analysis-Steel) has been developed to analyze simple or
continuous steel multigirder highway bridges (4). A number of
benchmark studies have been conducted to verify the accuracy of the
program (2,3,4). It has been found that the program can predict
the full load versus deformation (or damage) history of a given
51
Interface Slip
The slip between the slab and the girders could theoretically
be a major problem, especially at increased overload levels. This
problem could be further compounded due to the recent trend in the
8y
---~
~
z
'w
Fig. 1
............ 8
~.
52
increased use of "partial shear connectors."
Detailed analytical research, and its correlation with the
available test data, have demonstrated that for fully composite
construction the slip is not a progressive phenomenon for increased
overload levels (8).
In the case of the use of partial shear connectors slip does not occur under service loads, not even at load
levels in excess of service loads. For higher load levels, a slip
of a finite amount can take place, which causes a slight "jump" in
the load-deformation curve, corresponding to a slight incremental
loss of stiffness (8). The load deformation curve continues to
climb, with the previous slope, without any further disturbance,
as far as the slip contribution is concerned, until the spread of
extensive damage to the superstructure. Additional finite slip can
also take place. It has been noted that for fully composite or
partial shear connection cases the neglect of the slip will not
result in any major inaccuracy in the prediction of the post-linear
elastic response of the bridge superstructure. The result of this
simplification will be the slight overestimation of the stiffness
of the bridge.
In the case of fully noncomposite construction it is noted that
there still exists a substantial initial composite interaction
between the deck and the girders at design load levels.
For
increased overload levels there will be intermittent slip of undefinable magnitude. The slip, in this case, again is not progressive. For increased live load levels the frictional forces
between the deck and the girders increase, thus arresting the
slip in part.
In any engineering analysis and research the contribution of the slip should preferably not be overemphasized. The
complications arising from its inclusion do not warrant the increase
in accuracy.
3.2
Material Nonlinearities
53
Stability Considerations
The research have indicated that for rolled or built up sections the web buckling, with the exception of at the supports,
could easily be retarded by vertical stiffeners.
In the case of
plate girder construction, especially those with deep and slender
webs, the web buckling takes place almost at the same load level
that some relatively deep cracks, for example, one-fourth the slab
thickness, form in the deck. Regardless of the negative connotation of the web buckling, this is a local constrained phenomenon.
After the web panel buckling, the "truss action" takes over (3).
The post buckling strength of the web is usually more than sufficient
to permit the traverse of infrequent overloads (2,3).
In the case of continuous construction lateral torsional buckling is a possibility in the negative moment region. However, even
though this buckling mode is theoretically possible, due to the
contribution of the wind bracing, and X-bracing or diaphragms, the
lateral torsional buckling becomes an unlikely phenomenon (2).
3.4
Serviceability Limits
54
The key issue for the bridge engineer is to define the extent of
permissible "recoverable damage" to the bridge. In the US bridge
engineering practice, i.e. AASHTO Bridge Maintenance Manual, the
allowable stresses for the steel girders have been set low enough
to prevent the damage to the steel superstructure. It should be
noted that these low stresses are set forth to account for the
adverse effects of the residual stresses and the stress concentrtion at the connection details. The current serviceability limits
require extensive in depth studies in order to relate the possible
post-linear elastic behavior of the bridge components to the bridge
life expectancy as well as the permissible load levels.
EXAMPLE BRIDGE
~W36lt170
(Typ)
Ilj48m.m\ .2540mm
(57')
(100")
48m,!!' I
I .2540mm.
I .2540mm.
I11(57')
000")
(100")
Fig. 2
Cross Section
55
DEFLECTION AT SECTION A-A (In.)
10
20
r---------~r---------~r-~====~~!;:;::~~::l===~1200
5000
4000
-;;;
Co
0
<I
-'
800 :;;;
0
f-
-'
-'
<I
af-
3000
2000
o PrOf>lrom BOVAS
-'
f-
400 a
f-
o
DEFLECTION AT SECTION A-A (mm)
Fig. 3
GENERAL OBSERVATIONS
(2)
56
(3)
(4)
(5)
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
Fisher, T.A., J.W. Fisher, C.N. Kostem and D.R. Mertz. Design
and Retrofit for Fatigue Damage in Web Gap, Fatigue of Steel
and Concrete Structures, IABSE Colloquium Proceedings (Lausanne,
Switzerland, 1982) 535-543.
Hall, J.C. and C.N. Kostem.
Inelastic Analysis of Steel Multigirder Highway Bridges, Fritz Engineering Laboratory Report
No. 435.1 (Lehigh University, Bethlehem, Pennsylvania, USA,
1980).
Hall, J.C. and C.N. Kostem.
Inelastic Overload Analysis of
Continuous Steel Multigirder Highway Bridges by the Finite
Element Method," Fritz Engineering Laboratory Report 432.6
(Lehigh University, Bethlehem, Pennsylvania, USA, 1981).
Heishman, C.A. and C.N. Kostem. Users Manual for Program BOVAS Detailed Version, Fritz Engineering Laboratory Report No. 435.2
(Lehigh University, Bethlehem, Pennsylvania, USA, 1982).
Kostem, C.N. Further Parametric Studies on the Overloading of
Highway Bridges, Fritz Engineering Laboratory Report No. 434.2
(Lehigh University, Bethlehem, Pennsylvania, USA, 1980).
Kostem, C.N. A Parametric Study on the Overloading of Steel
Multigirder Bridges, Fritz Engineering Laboratory Report No.
435.3 (Lehigh University, Bethlehem, Pennsylvania, USA, 1983).
Peterson, W.S. and C.N. Kostem. The Inelastic Analysis of
Beam-Slab Highway Bridge Superstructures, Fritz Engineering
Laboratory Report No. 378B.5 (Lehigh University, Bethlehem,
Pennsylvania, USA, 1975).
Tumminelli, S.C. and C.N. Kostem. Finite Elements for the
Elastic Analysis of Composite Beams and Bridges, Fritz
Engineering Laboratory Report No. 432.3 (Lehigh University,
Bethlehem, Pennsylvania, USA, 1978).
Standard Specifications for Highway Bridges. The American
Association of State Highway and Transportation Officials,
Washington, D.C., USA, 1977.
Manual for Maintenance Inspection of Bridges. The American
Association of State Highway and Transportation Officials,
Washington, D.C., USA, 1978.
57
C. P. Heins, Professor
Institute for Physical Science and Technology and
Civil Engineering Department
University of Maryland, College Park, MD 20742
INTRODUCTION
During the past ten years extensive research studies have been
undertaken in the area of the design and analysis of curved bridge
structures.
It is therefore
ANALYSIS TECHNIQUES
The determination of the induced forces in a horizontally curved
bridge system requires consideration of the interaction of the entire
bridge structure.
58
20
30
3,4
30
40
40
50
or more
2 0.4
~ ~
10.
be neglected.
where:
L[G~/EIW]1/2
2.
Three DOF
b)
Five DOF
3.
4.
(29)
(23)
(23)
(23,24,28)
(23,24,28~33).
59
5.
(8,9,10,15)
The method
Noting that the resultant forces (F) in the flange are not
1
h
(1)
however:
H h
therefore;
T
T
therefore;
where:
!i.
R
H
R
and
(2)
60
EM
(3)
3D
11 6
v1 - 1.3
and
3D
VI
x 2
"6
VI - D + VI
VI
_ 10 D
9
VI
E~
V2
Vl/3
!:! ~
R
D
~
"9= E !:!
R
- 1-11I D
(4)
(5)
E M~
where
= E !:! ~
R
R ~
In general the
D
+ V2 - "6
- E
RD
(6)
Girders
1.0
1.0
1.11
5
1.25
6
1.40
1.56
10
1.72
1.88
2.04
61
degreees of
freedom by two.
The space frame grid method, consists of 6 DOF, and considers
the girder flanges modeled as beams and the webs as a series of
cross and vertical elements.
in
"
EI
''''
w <P
EI
"
+~<P
R
(7)
62
EI
-~Tl"" + (
R
EIX+G~
) n" - EI
9'"
EI
- ~ + m
R2
z
+ GK "
T
(8)
where:
Vertical Deflection
Transverse Rotations
E1x
Bending Stiffeners
G~
Torsional Stiffeners
E1w
Warping Stiffness
= Radius
r2GK~
T
[~'
2t:,2G~
("+EI )+
. [ EI} .2( El +0"'1 ~'7 2.2(EI2+~l
4-f-H1 x
4\R
it-2
n-l
EI
~ + EI
R2
EI
w
R
n+2
R2
R +
r~Iw
-[,t:t.
-R2
n+l
n
=~t:,
} .2~J
2(EI; .o",l]
- 4 -+EI
R2
(9a)
63
El
w
R
-[4~~l + .2(EIX:eR,16~~~ +
w
n - 2
114El
x
6(El) + 2112GK.r + _ _
R2
n
El
2112(E1w + GK.r)
n - 1
-f~~~ + "(EYeR,
-[4( E1W) + 112GK.rJ
n + 1
E1
w
R
m 114
(lOa)
El
w
n + 2
or
[Kl]{n} + [K 2 ]{<P}
114
qy
(9b)
m 114
(lOb)
where the
node~points
~a
a
na
QDa
~n
QDn
and
[K]
is a
diaphragms.
4x4
[K5 ]
(11)
1h
64
~n
4>n
nn
QDn
~b
[K 6 ]
(12)
4>b
QDb
The total diaphragm
4>b
"n"
can be
~J [::J
=
where:
(14)
(15)
and
K ..
~J
= force
and (lOb) represent the general response of curved girder #1, subjected to a
uniform loads qy and mz '
At node
"n"
65
Mnn/n and
QDn/n.
Therefore
(17)
where:
[K 7 ]
(K8 ]
[K9 ]
[K3l , K13 ]
(K lO ]
[K32 , K14 ]
(Kn]
(K 12 ]
[K13 ]
[K4l , K23 ]
[K14 ]
(K42 , K24 ]
{Ij>}
{n}
G:J
,
=
[~: J
If the node at which Eqs. (16) and (17) are being applied does not
correspond to a diaphragm connection, then the
K ..
1J
to zero.
The above referenced equations have been incorporated into a computer
program (31) and utilized in developing design aids, as will now be described.
66
DESIGN AIDS
Ix'~'
and Iw.
The designer,
and bimoment
which induce normal stresses, were known with some degree of confidence.
Then the basic bending equation,
x
I
(18)
BiW
o
Bi
(I
(19)
Iw )
w n
Ie)
0b
and (I Iw ) properties.
w n
67
These techniques
for the bending, Fig. 7, and Fig. 8 lists the impact factor for the
torsional effects.
e=
(L/R).
Bi
MSTATIC(l + IF)
(20)
BiSTATIC(l + IT)
(21)
where
M
Bi
68
be used to establish the induced stresses and thus, proper girder section.
Amplification Factor - K1
All the internal forces and deformations for a single curved and
straight girder have been evaluated using the various computer programs
(33,36).
The ratio of the reactions for these two girders gives the
following,
f(S.C.G.)
f(S.S.G.)
(22)
its analysis gives the following general equations (23) and (24),
Kmoment
K
-bimoment=
0.15(L/R) + 1
(23)
(24)
where
n
R/100,
R > 100'
Radius - Ft.
Distribution Factor - K2
The evaluation of the true distribution of load to each girder,
and this realistic values of internal forces, can be considered by
analyzing the curved girder as a system.
69
The ratio of
f(Sy.C.G.)
f(S.C.G.)
A plot of this ratio vs.
Rand
(25)
These equations are valid for girder systems which contain four, six,
and eight girders spacings of 7', 8', 9', and 10'.
Kmoment
(n
+ 3) (OR4L) + 0.6
K
- .11(R/L)
-bimoment- M
L < 70'
K
- (M - l)(R/L)
-bimoment6
L > 70'
(26)
(27a)
(27b)
where
R > 100'
R/lOO,
L/50,
Radius - Ft.
L > 50'
Reduction Factor K3
Because many bridge structures are continuous, it is desirable to
obtain some factors which can be applied to the simple span data to give
70
f(Sy.C.G.)N
f(Sy.C.G.)
(28)
where
N
Number of spans (2 or 3)
K3
and three span girder systems must all contain equal span lengths with
a maximum given span length of 100.0'.
= 300.0'.
This
71
3.
K2
4.
K3
span.
5.
(2
and
Max. Bimoment Static
6.
(30)
S,
where
1.29
for
7' or 8'
1.57
for
9' or 10'
Rand
This technique
72
curves.
warping stress
S.
(0 )
w
and
where
0b
0b
Cl
sIs
ow
a
ow
(32)
Referenced Value
Design Value
The design charts, Figures 9 through 14, list the induced stresses
and
73
length (K/in.).
in an actual design.
= 100',
These plots
spacing.
The plots of
vs.
(R/L) for L
vs.
(R/L)
(I ), as determined
x
(~).
These curves
x -rx
I
(33)
where
rx
Actual Stiffness
Reference Chart Value given in Figures 10, 12, and 1.
.yot_.
": )
chart value _ t b.
~difi.d.
74
()
a
/j
chart
(- )
cr
(34)
where
"C':)=
(": y
(35)
I
I
w
w
Girder Spacing
The curves have been developed for girder spacing of 8.0', however
spacings up to 10' are acceptable.
Torsional Properties.
If the girder
75
slab thickness, as shown in Figure 12, with the dimensions defined
as follows:
G /G.
c
Beam:
W
n
(37)
c
(38)
Warping Stiffness
I
(39)
Torsional Constant
1
(40)
"3
where
m
G /G
76
Slab:
Steel:
BiW
n
(J
(41)
BiW
n
___
5
(J
(42)
DESIGN SPECIFICATIONS
The AASHTO design criteria (1,2), as applied to steel I and Box
girder bridges, is given in Table 2.
straight and curved elements.
77
~------------------------0
100
Mono bo.
"2:30
0
\()
30
.,-:
1.0
-----------------------------~
Multiple I
ITT
'.0
It-O.S-3
~~
D5~
0.1
0.2
0.1
0.'
O~
Cent-,I angle.'
0.'
OJ
0.1
0.9
Figure 1
78
Fig. 2
lq , lDl
~
0/3
Ie
.I'"
1\
0/6
,.'t:
0/6
I
.,
1\
0/3
Fig. 4
n-l
n+1
~2
,/ \.-I
Fig. 5
~Ol
_l
~
n-
nH
n+2
~""O<"
Fig. 6
79
.t
I
80
".
200.
I"'.
I!I C.S.
JO.
o.
I.
..
1.
zo.
10.
la.
Figure 9
IOt1OUZI!D IWU'IHl mESS
L-lOOO'
12,000. ~It
Ia -
"".
400.
lOO
200
1/L-2,S
100
IlL -
---.--24
-r-I-22.
20.
18
16
l~
.
I
figure 10
81
mas
lQIW.IZID IEHJDC
L .. I25.D'
lID.
UO.
200
uo.
1\;"10
.ao.
...
..
z.
1.
...
zo.
Figure 11
.................... mas
L-WO'
x.." 10,000
m4
5<10
...
lID
IlL" lot
'00.
lO'
'4.
22
20.
II
Ii.
14.
Figure 12
82
1GIIIiIl.mq,~_
L-UDO'
~.U5.OIIl ~Io
F~f:ute
14
Y[SE<:T1O.
IC)DIJUIS
.,.. sp_
...-l
83
m)
La:ilth
l.
2.
1.
10.
50
no.
90.
130.
ISO
10
...
...
I~
.or.
...
figure 15
100.
60
60.
t,..Spm (fcJ
100.
JOQ
250.
x..' ",000
100
..4
1. ... 125
Z50,
lO.
200
;...
lb~
O.
1lO.
.. l2,ooo
~!1..50.
lO.
O.
5.
10.
--60l<don
--IGUdc.
.0
. 0
..'
't
..
OIL
5.
10.
I
I
In'
L .. lSQ'
\
"
I}
f/--ol
bt
x..' W."'"
100
SO.
O.
l3O.
L 100'
;...
100.
140.
;... 200.
I~P~.
JOQ.
1Jll.
OIL
1
Figure 16
..
.,
Figure 17
Three
Span
Span
Two
.70
.70
.70
.60
.75
.65
.35
.35
.45
.40
.75
.75
.90
1.00
Number
of
~ending Moment ~eflection ~otation KSt . Venant ~arping Torsion ~imoment KShear
Spans
K3 - REDUCTION FACTOR FOR MAXIMUM FUNCTION IN TWO AND THREE SPAN BRIDGES
TABLE 1
-1:0-
00
moment)
(pos! ti ve
it
Flange
Compression
Item
t~~~
Straight
TABLE 2
< 4400
--to
(b)
1.
0.6 ( r )
b
[30+8000(0.1-R )
fw
and Pw
0.95
-(-r) [1"""]51
b
--f
/
b
Pw2 =
iff
or
Pw1 =
+ (!:.)(!:.)
R b
(!,.) 2F
0.55 F [1- ~1p P
Y
41T2E
B w
PB
where
Fb
and
t~~
Curved
VI
00
bJ
moment)
Compression
Flange
(negative
Item
I fy
6140
k
Y
< 60 or 13,300
< 60
t~
IF
13300 <
(~)2
fb :'. 0.55Fy
_ .
0 224F [1 -sin _11
Y
2
~~t~
6140 <
fb ~ 0.55 Fy
Straight
TABLE 2
- b/ t IF
( 1330071Ln
y
)
=
0.55 F
IF
y
y
6140 . X
I
f
< 13,300 or 60
Y 0.224F { sin
0.326Fy +
y
t~1F
if
t
IF
13300
11
"2
t 2 x 106 _
57 6(.!.)Z1066
. b
Fb - 57.6(1;)
F
b
k~JJ
113.4(-~)2
f 2
x 106
13,300 - t
(13,300 6140X )
k~
6 = 11_9.0(;V)2
y
where r---~-
FbY-
IF
6140 <
IFy
1-9.2[/1
y
Curved
fv
where X = 1 +4/3(- - 0.15) ~1
Fb
k:'.
0-
00
~:~
10
Stiffener
with
moment)
Compression
Flange
(negative
Item
iK
~ 0.S5F y -
t~
>
B=
~t3w
0.125K
for n=l
where
3 4
$ ={O.07K: for 0>1
Stiffener requirement
IF
y
3580/K
~]
t
y
Xl
Xl
where
Fb
<
Curved
<
(n
1)2
K
fv
0.93 + (1.6 - ({)(y)
1-912(/)
1 (n>1)
0.5SF
IF
~ < 3070/K Xl
0.224F y [1 - Bin!
2 66S0/K -
fb
IF
fb ~ O.SSF y
IFy
~ ~ 3070
TABLE 2
Straight
5.34
1 (n=l)
--.J
00
hJ
with
stiffener
moment)
Compress ion
Flange
(negative
Item
Straight
TABLE 2
t~
O.l{~s)t~
or
Fb
14.4K(;)26 106
= 14.4K(.!.)2 x 10 6 -
Fb is smaller value of
IF"y
w
~
106
66S01i< X2 -30701i< Xl
66S01i< X
2 -
v K
2 t 2
14.4(K.) (~)
f 2
(f)
in
ylf i
S.34J
66S01i< X2
< ."': < 60
+ 0.224F
f 2
1-9.0(f)
y
3Z6F
rLO. y
I
IF
Curved
X2 ~ 1 - 2.13(F~)
y
Where
'b
IF
30701i< X
}] 6
00
00
Stiffener
With
Long1 tud Ina1
and Transverse
Item
> 8t 3 w
sAf
a
~
E
a:
where
I t "- o.10(nH)3w3
Straight
TABLE 2
Curved
'"'
00
ll
Transverse
Sti ffener
With
L.lb= =
Stiffener
r:n:.
Without
Web
Item
fv -
do
0.B7(1-C)
stiffener spacing
F (d /t )2
<J[C+~
2
3
( old)
d0 < 1. 5d
t~
t -
~ < 150
Straight
TABLE 2
::'1.0
0.02
F (d /t )2
Y
+ (dId )2 J
0.B7(1-C) J
,d---:L
11 +\-o/d)
C = 2.2 x lOB [1
[C
1. 5d
Fb
- J
<
< 1.0
fv - 3
%::.
Curved
If d O/R
IIf
Same
)g
with
Transverse
Stiffener
Web
Item
d
25 (--%)2 _ 20 > 5.0
rl t 3
> 0
- 10.92 J
Stiffener Criteria
Straight
TABLE 2
< 2600
!?
n;
d2
0.95flt
(/
~d
= 1.0+
~)2
Curved
0.78
- 0.78) 4
)
1775
Z
2.
d
1.0 for dO
[25(
t3 J
> d
Stiffener Criteria
0.78~-f~1.0
92
REFERENCES
1.
2.
3.
4.
5.
C.P. Heins, "Box Girder Design," AISC Engineering Journal, Vol. IS,
No.4, December 1978.
6.
C.P. Heins, L.J. Hua, "Proportioning Box Girder Bridges," ASCE St. Div.
Jour., Vol. 106, No. STll, November 1980.
7.
C.P. Heins, D.H. Hall, "Designers Guide to Steel Box Girder Bridges,"
Bethlehem Steel Corp., Bethlehe~, PA, 1981.
8.
9.
10.
C.P. Heins, J.C. Olenick, "Diaphragms for Curved Box Beam Bridges,"
ASCE St. Div. Jour., Vol. 101, No. STlO, October 1975.
11.
12.
13.
14.
C.P. Heins,et aI., "Curved Steel Box Girder Bridges: State of Art,"
ASCE St. Div. Jour., Vol. 104, No. STll, November 1978.
15.
C.P. Heins, W.H. Lee, "Curved Box Girder Field Test, " ASCE St.
Div. Jour. , Vol. 107, No. ST2, February 1981.
93
16.
17.
18.
Heins, C.P., Huo, C., "Truck and Wind Load Effects on Girder Slab
Bridge With/Without Wind Bracing," AISI Project 308, Effects of
Bracing on I Girder Bridges, Inst. for Phy. Sci. and Tech.,
Tech. Note BN-988-A, University of Maryland, College Park, MD,
April 1980.
19.
20.
21.
22.
Heins, C.P., Kuo, T.C., "Live Load Distribution of Simple Span Steel
I-Beam Composite Highway Bridge at Ultimate Load," Department of
Civil Engineering. University of Maryland, April 1973.
23.
Wilson, L., Bathe. K.J., Peterson, F.E., Dovey, H., "User's Guide
to Program SAP," University of California. Department of Civil
Engineering, November 1972.
24.
25.
26.
94
27.
28.
29.
30.
31.
32.
33.
34.
35.
Vashi, K.M., Heins, C.P., "Impact Factors for Curved Highway Bridges".
Civil Engineering Rept. No. 32, Spetember 1969.
36.
37.
38.
39.
95
H.R. Evans
SYNOPSIS
The naner describes a number of simplified methods that have
been developed for the analysis of cellular bridge decks; the
behaviour in the linear elastic and in the non-linear nost-buckling
range is considered. Each method assumes an idealized structural
behaviour to provide a simplified solution that may be carried out
on a microcomputer. The methods have been develoned specifically
for use during the nreliminary stages of design when reneated
analyses may be required to establish ontimum nroportions. The
paper includes some typical comparisons of the results given by
the proposed methods with the results of other established methods.
1
INTRODUCTION
In recent years there has been a great increase in the computing power that is available to the bridge designer and this has
made it ~ossible for powerful techniques, requiring extensive
computations, to be applied to the analysis of bridge structures.
Of major importance in this context is the finite element method (1)
which, by virtue of its adantability, enables a wide range of
bridge loading, geometry and sunnort conditions to be taken into
account in analysis.
Despite the increase in available computing Dower, the finite
element method must still be regarded as a rather exnensive method
of analysis. Certainly, at the initial design stage, when the
designer is endeavouring to ontimize the prODortions of his structure by repeated analysis, the costs involved in using the method
are unaccentable. Furthermore, since most modern design codes are
based on limit state considerations, it is essential to be able to
96
97
98
nodal sections
transverse
frame action
~
H [
interaction
between frame ;
and
:
'I
plate systemsl=t:
:
;:
I
,I
,
longitudinal
plate action
FIG. 1
'"
M = 6 EIlI
12
FIG. 2
100
II
Cal
IDDlce
(a)
(b)
(c)
FIG. 4
______other methods
G nodal section
101
'?
ie
10
~...
-- ---
1000
so
.-.- _
040
"
III
>
20
10
!II
.=
_ _ - -'>- -
It"
,-
;'
v'
2'0
1ii
radius
span
1'8
I'"
--
III
....A
.c.
e 1.4~
III
-I/,- _
inner edge
"
2
1
...
outer edge
"JO
"2
,0
g 0'&
CII..Itu edge:.
:;; 0'"
:l
"'=
0'4
WO.2
'l.
....J
radius
span
radius
-span
l'
.e
V60
...u
.2
..
<;:::
til
'0
-0u
'"
'50
'40
.~
20
outer edge:
'10
>
FIG. 5
102
2.2
103
stress distribution
due to shear lag
FIG. 6
(a)
single-cell box
spine-beam girder
double-cell box
FIG. 7
104
112
Finite element
Bar Simulation
span/width
FIG. 8
f-!
015
hoi
Interm_cHace
suppo"
48
Two-span girder
FIG. 9
II
105
[5t...=:',
II'
II'
d:::::J~
==",,==",==m
Spu _ 12m
Span_<tOm
FIG. 10
FIG. 11
GRILLAGE IDEALIZATION
106
107
108
"'-.,,"--r<'-::..-
~_~+-1"--(,
/!:::~===:?'
FIG. 12(a)
W
v b_
, ,'l,.<:9
ij
,I,
~n
! .
I
1"-
1-2.em
,,-1.~ em
I'
*\.2~1'2em
,_
25P
60
I.
2'5P
40
20
25P
o
3
5
0
10
2
- - FIn;le elemenl melhod
o
Grillage me1hod
FIG. 12(b)
LONGITUDINAL
STRESSES (N/mlh
3
4
FIG. 12 (c)
VERTICAL DEFLECTIONS
Crrun)
109
~O~------~------~~'-----~~------7.40
Web opening: 'III
(a)
CKS)
~O~------~~------~20'-------~~------~40
Web_nlng:'III
Cb)
FIG. 13
CKr)
110
"
~ .......
-
1'()
~,
2-0
,----,
1IJ
-,-i
4-0
-\..
-" "-
........
C---
(l)
lH
S-O
(l)
--- --------j
x --_
3-0
-,-
................
---
-, _
I.
'--l
lor
&0
,----,
200
-Z
- - - __
_
J
...
. "~-~~
'-'
,S-o
1IJ
<l)
1IJ
1IJ
<l)
/'
;...
10-0
C/J
/'
--
/-/
,.
../
5-0 -- --
()o()~~-----------------~~~------------------~~~
6-0
"-0
3-0
tl ------
2-0
Jf",.,.""'''
'"
... ..-
"
"
x
Ala
,"
_-8
80-0
- - GRllEJ re8ulll
o
Finite element reautll
60-0
400
2()O
'-a
(b)
10
24-cell girder
--------,,.,,,O,...--------.iiOI<-------"3""h---
O-O O
!:-
(a)
OpenIng: '1ft
16-cell girder
FIG. IS
III
112
113
Flange edgeltr...: N/mm2
8~0~____~4~0____~80~____~1~20~__~1~80~__~2~OO~__~2~40T-___
:li6000
Can".1 CSellexlon: mm
Cen"al dolitxlon 01 tllo ,truclU," with thlckor wOb,
.round tho porlphe'Y
Flange Idoe Itrell a' the centre 01 the Itructure with
thlckor wobi around the poripho'Y
-I
-I
1-
sr.
:t
2-
bendinG
Shed I
~t
strcs~
..ess
rIItio
u/o v
...
0.~2
0.86
0.84
1.00
0.57
1.00
0.78
0.67
85 kN
1.00
l/T v
1.00
0.~8
1.00
1.00
0.78
1.00
0.00
0.34
1.00
10
0.60
0.99
stage 3
80 kN
ratIo
.:;
2-
stage 2
Ii
1.00
0.~7
12
1.00
0.00
114
comparisons that have been carried out to date have shown encouraging agreement but the author is currently conducting a series of
tests on large-scale models to provide further data.
The extension of the grillage approach has been carried out
without adding greatly to the complexity of the solution process.
Consequently, even in its non-linear form, it may still be regarded
as a convenient design aid.
4
CONCLUSION
A number of different methods of analysis for cellular bridge
structures have been described in this paper. In each method, an
idealized structural behaviour is assumed so as to provide a simplified solution procedure, making each method suitable for use during
the preliminary stages of design when repeated analyses may be
required to establish optimum proportions. Some typical comparisons
with results obtained from other methods have been presented to
illustrate the accuracy that can be achieved by the simplified
approaches; such comparisons are presented in more detail in the
listed references. All the methods presented are suitable for
solution on the types of microcomputers that are becoming increasingly available to the design engineer.
REFERENCES
1.
2.
3.
4.
5.
6.
7.
8.
115
9.
119
Alex C. Scordelis
Professor of civil Engineering
University of California
Berkeley, California USA
SUMMARY
A review of the analytical methods and computer programs, which
have been developed at the University of California at Berkeley for
the linear and nonlinear analysis of reinforced and prestressed
concrete box girder bridges, is presented. Short descriptions of
15 different computer programs based on one of the following methods
of analysis are given:
(1) folded plate elasticity method; (2)
finite strip method; (3)
finite element method; (4)
finite segment
method; and (5)
nonlinear analysis method. These various programs
provide the capability for the detailed analysis of single and
multicell box girder bridges of one or more spans, which. are straight,
curved, skew or of arbitrary geometry in plan. General loadings and
arbitrary boundary conditions may be treated.
The linear analysis programs assume the analytical model to be
a linear uncracked homogeneous concrete system. The nonlinear
analysis programs utilize analytical models which include the
concrete, reinforcing steel and prestressing steel and can predict
the structural response of bridges, under initial and time dependent
load histories, through their elastic, cracking, inelastic and
ultimate ranges.
Five representative numerical examples are presented to
illustrate some of the solutions which can be obtained using the
Berkeley computer programs. Results obtained are compared with
experimental results and those from other analytical solutions.
120
1.
INTRODUCTION
121
FIG. 1
Ny
FIG. 2
Nxy
Nx
MXY~MXY
FIG. 3
, Rtr
~
I
FIG. 4
(JOINT
.... \'j"j
(i=-ijj
\~J
122
approach has been to: (1) study the available literature; (2) develop
analytical methods and general computer programs; (3) perform
experimental studies on elastic models to verify the analytical
methods developed if deemed necessary; (4) make analytical parameter
studies; (5) test large scale reinforced concrete models or prototypes; and (6) develop recommended design procedures.
Initially,
reinforced concrete bridges and linear methods of analysis were
studied, but in recent years since 1977 the emphasis in this
research has shifted to prestressed concrete bridges and to nonlinear
methods of analysis. The results of this continuing research program
at the University of California have been described in a series of
research reports and technical papers which are listed as references
in Appendix 2, under A and B respectively.
The purpose of the present paper is to review the analytical
methods and computer programs, which have been developed at Berkeley
for box girder bridges, and to present numerical results from
several examples. Both linear and nonlinear analytical models and
methods will be described in the body of the paper, while in
Appendix 1 short descriptions of 15 different computer programs for
box girder bridges, which are available, will be described. Detailed information on these computer programs and their theoretical
basis can be found in the published research reports and technical
papers listed as references in Appendix 2.
2.
123
3.
124
computations, has been devel~ped for a single harmonic, it can be
reused for any harmonic, and thus the method is ideally suited to
the application of a digital computer.
The analysis for each harmonic load has the advantage that for
straight or curved circular bridges such loads will produce displacements of the same variation and vice versa and thus a single
characteristic value may be used to describe any force or displacement pattern. For example, the displacement pattern for the nth
harmonic:
mrx
rex)
ro sin
reel
ro sin
e-; for
mr8
(1)
(2)
eo
2.
125
3.
4.
126
(a)
----1
-(b)
(c)
Primary Structure
~OLDED
~LATES
RIGIO LINKS
~+F------'~4r ~
GIRDER AXIS
SUPPORT
BENT
(d)
(e)
127
the folded plate structure is analyzed for the
given external applied load using a direct stiffness
harmonic analysis as described above. A displacement vector 00 (24 x 1 in Fig. 5) is found for this
case, which defines the displacements at the points
where the redundants are to act.
2.
3.
4.
(3a)
(3b)
or
00
FX
(3c)
where
(4)
x
5.
F -I
(3c)
00
(5)
128
4.
129
L 20
82
.. TwVzz
Uz
-I
"
FIG. 6&
EDGE I
UI
y..... y
"'1
X,U,X
III
w,
81
F
&l
I,
Z,W
= K 10
Y= y/b
J(
I
I
.-~IIIIIIIII~
'
COS nTTi!
[-"-1-...:"
....
~
FIG. 6b
130
FIG. 7a
FIG. 7b
131
4.2
132
rve d Br idg e
Ci rc ul ar ly Cu
(a)
,/
,.
"
(h)
.......
_----
Ri ng Pl ate
El~ent
,
,r,. (C)
Co nic al El~ent
FIG . B
(d)
'-;.~-.,
me nt
Cy lin dr ica l Ele
E ELEMENTS
GIR DE R BR IDG
CURVED BOX
133
(3) the standard friction formula is used to calculate the prestressing tendon force Fi existing in each tendon segment.
F.
l
(6)
5.
Both the folded plate elasticity method and the finite strip
method make use of the direct stiffness harmonic analysis procedure described earlier. Advantages of this procedure are:
1.
2.
3.
4.
2.
134
SECTIONS
TENDON POINT
.......
PORTION (J)
10
FIG. 9
SEG~'ENTS
PORTION (1+ 1)
2 SEGMENTS
SEGMENT:
,,,
,
CD
TENDON POINT:
TENDON ANGLE CHANGE:
0.
FIG.
lOa
FIG. lOb
135
i,j
FIG. 11a
FIG. llb
= STRIP EDGES
= TENDON POINT
136
6.
Many multicell box girder bridges used in separation and interchange structures fall into this category. The computer program
CELL was written especially for this type of bridge in which the
bridge is made up of top and bottom horizontal deck slabs and
vertical longitudinal webs and transverse diaphragms (Fig. 14).
The deck slabs are idealized by quadrilateral elements having
a total of 5 degrees of freedom (DOF) per node, 3 translations and
2 rotations (Fig. 13a, b). The in-plane action of the quadrilateral elements is represented by the plane stress mixed model
Q8Dll having 2 translational DOF at each external corner node and
3 internal DOF (Fig. l3a). The mixed model is constructed using
separate expansions for the displacement and strain fields.
The
variations of the u and V components of the displacement field are
approximated by the standard bi-linear expansion for the 8 corner
node DOF and by bi-quadratic expansions for 2 of the internal DOF.
137
FIG. 12
J,Y
YI.k+
8Yll~
~Y.Y=CONST.
L -_ _ _....;Z"--__
la)
J,V
-t=i-
8z 1
u"-tv
8., w,
Id)
1,8x
~I""
u~
u = a [ 1-)
a
1.
I,U
FIG. 13
r:
.,U
"~~
&
.,8x
::I
f"2
J,Y
IG)
Ib)
8z 1
'-----~_~
Ie)
U,
WI~
.,u
x2
21-)
a
~~ .L
I f)
I~)J
a
21~)3
a
Zl
31~) 2 + 1
a
138
PLAN VIEW
E!
ELEVATION
TOP DECK
II
II
,I
F=RTICAL WEB
BOTTON DECK/
SECTION
FIG.
14
A- A
139
The third internal DOF is used to enforce a constant shear-strain
variation over the entire element, which produces a more flexible
and better element. After the element stiffness is formed the 3
internal DOF are eliminated by an internal static condensation
process. The quadrilateral plate bending element Q19 (Fig. 13b)
used for the deck slabs is a compatible element made up of four
subtriangles, each of which has 11 DOF associated with full cubic
expansions of the w-displacement field and an enforced linear
variation of the normal slope along one edge. In combining the
four sub-elements, a quadrilateral with 19 DOF is obtained. However,
the 7 internal DOF are eliminated by static condensation leaving the
essential 3 DOF at each corner node, 2 rotations and a translation
(Fig.13b).
The veritcal webs of the bridge are idealized by special
rectangular spar elements having a total of 5 DOF at each corner
node, 3 translations and 2 rotations. A single element over the
entire depth of the bridge can be used to capture the essential
behavior of the web. The in-plane action of these elements is
represented by the model QUSP12 (Fig. 13c). A bilinear expansion
for u and v is associated with the 2 translational DOF at each
node and a cubic variation in the x-direction of V is defined by
the rotation 6 z i = dV/dX at each node. The plate bending is
represented by a simple one-way bending element ONEW having 2 DOF
at each node, a rotation and a translation (Fig. 13d). Orthotropic
plate properties and arbitrary loadings and boundary conditions
can be treated in CELL. Automatic element and coordinate generation options which minimize the required input data, have made
CELL a widely used program.
6.2
140
each node is also assigned a rotational DOF (Fig. 13e). This is
defined as the average rotation about the element z-axis. The
internal element displacements u and v are assumed to vary linearly
with the 2 translational DOF and as beam functions (Fig. 13f) with
the rotational DOF. The assumed nodal rotations introduce small
angular discontinuities so that the element is not fully compatible.
The plate bending model used
to that described in the preceding
program. The stiffness matrix for
frame elements can be found in any
7.
141
3.
8.
1.
2.
3.
4.
142
Each plate segment is assumed to obey the ordinary theory of
folded plates. Thus slab action is determined by the behavior of
one way slab strips spanning between longitudinal joints, so that
torsional and longitudinal slab moments are neglected.
In-plane
membrane action is defined by elementary beam theory applied to
individual finite segments which results in a linear valuation of
longitudinal stresses over the width of each segment between longitudinal joints.
Advantages of the finite segment method are:
1.
2.
9.
1.
2.
NONLINEAR ANALYSIS
143
144
(a)
(b)
Element
---11lZlZl'8IID.~
---I'---'+'-~--...REFERENCE
PLANE rl~~~~PREST,RESSING ~
STEEL
(e)
(d)
(e)
FIG. 15
145
146
(I)
-cr
1.--' ]
(I)
LLJ
FIG. 16a
('1)
('1)
\T
'--IT-'
(4)
.-T"T""I---' .-r---.-, ,-
(~)
LT
Cross-S~ctlons
~)
-,;;::::r-
ElevatIons
~)
z-.----zz
G)
('5)
'J""
,~--.--.
i..-~=--l"7
Plan
.j>.
.....
148
-----';1:
a
Box Cross Sectlon
rC.
!&
\;
-x"
Plan
FIG. l(,b
Lc
149
150
.t
I
---
",~
, 2 '
"
~]
TRANSFORMATION TO GLOBAL COORDINATES
DIRlCTIOllS OF
1'IWISVERSE SECTIONS
--- LOCAL
-GLOBAL
(al
-?- t:t]
r t13i': ~ t 3
B
ELElIENT DOF' S
(bl
2"L 4
\J
(c)
FIG. 17
d'-----
'2
.~
I
---- LOCAL
_
GLOBAL
3 2
Beam Element
SHEAR
PANEL
TOP
C;
(d)
3-D GRILLAGE KlDEL AND STRUCTlJRAL ELEMENTS FOR COMPUTER PROGRAM NOBOX
151
10.0
NUMERICAL EXAMPLES
152
A.2. In general, the numerical results from the computer programs
have been verified by one or more of the following ways:
(1) comparison with experimental results from physical model studies;
(2) comparison with numerical results from other analytical models;
and (3) performing statical checks.
Five representative numerical examples will be presented here
to illustrate some of the solutions which can be obtained using the
Berkeley computer programs. Additional numerical examples may be
found in the references given in Appendix A2.
10.1
0 Jill
~ L.-
"
-:"'-";===,R4
WEB NO. 1 :,
-,;-..;::c::::-:..-.-.-:.:.-.-.-.'ii/P-::::o;::::::::;::.-:.__-.-_-.::.:.:::_-.:..;
b)
Rl
r-~ I - - - - - . - - - - - - - - -
RJ
FIG. 18
d)
L _ _ _ _ _ _-1
0 AT 3" , '2"
0 49
I/B
--
-:.-i
/,
No-f--~-"- ' 1. /
F"'
110 DIAPHRAGM
I ~,'
'" I-
S~AN
e)
10"
--------l
60"
j~Fi~~~~:~:i~~tl~;~:;l.
c)
lLJCJ~B6pr~
---.---------.------WEB NO.
Ci:~~~-:~~~~~~~~:~::;~~~~H[~::~~~~~~~::~:~;::?-END DIAPHRA~
-'
_R2
SPAN L - - - - - - 1
3"SKEW ANGLE
1"
J _
a)
I,
VI
W
154
0.60
a)
= 45,
-0.90
_ _ EXPERIMENT
-o.eo.
--<>_
0-.0..
in)
I>
.'
"0
..0"
= 35.5
J~"
THEo.RY
All
~
-0.70.
Span
~-O-
0.-0-
.0'.-
-0..40
-0.30
-0.20.
~
-0..10.
0.
w
&
w
&
= 116
0..10.
0..20
0..30
o.Ao.
O.~O
0.60.
0..70
0.80
0..90
FIG. 19
b)
in, Span
= 60
in)
155
distribution that' is a function of span, skew angle or curvature;
and (3) the transverse distribution of My at midspan, for a given
load position, is essentially independent of span, skew angle or
curvature.
Typical theoretical and experimental results for the transverse
distribution of Nx at midspan for three midspan load positions are
presented in Fig. 19 for skew Model 3A having a skew angle = 45
and a span = 35.5 in. (902 mm) Fig. IBa) and for curved Model lA
having a radius of lIB in. and a span = 60 in. (1524 mm) (Fig. IBc).
Note the good agreement between theoretical and experimental values.
It can also be seen that the distribution curves for the two bridges
are eseentially parallel, as indicated in conclusion (2) above,
indicating that the departure from a uniform transverse distribution
is primarily a function of the local plate bending moments in the
bridge.
10.2
156
A X
II
WESTI
Y 0
lEAST
~======~F=======~========~=F======~
SPAN I
SPAN
n
3600
GI~ER
A X
-r
~
I
~) r, f-t
1442
~DIAPHRAGM+358
\300 30e.
1500
1
1
1500
T
I
358'
-r
I
1442
}
~~~==~~~~~~~~~~
_~~~~__~~~__~~~~~~__~~__~~~__~~___ ~
--~-----+~~-----4~----~-~!~
04MENSIONS
FIG. 20a
EXAMPLE 2 -
IN FT (1 FT .0305 m )
35/8"
~.I
FIG. 20b
5'-0"
II
i
L.
"
2'-6114"
FT.
IN.
=25.4 mm
=0.305 m
L,
l' -8"
l' -6"
ROUND COLUMN
iOiOI
12' - O
4(@2'-6 7/8"=10-3/1/2"
v.
-"
158
B21, B24, B25]. Extensive correlative studies of experimental
vAlues for reactions, deflections, internal forces and moments with
theoretical values predicted by various computer programs are
presented in the above references.
Linear elastic analyses were carried out for the straight
bridge using FINPLA2 and for the curved and skew bridges using
CELL. Only a typical comparison of theoretical results for the
three bridge models [A23] is given in Table 1 for midspan deflections
and in Table 2 for longitudinal membrane forces Nx at Section A
near midspan, for midspan point loads of 100 kips (445 kN)
The transverse distribution of deflections at midspan sections
X and Y given in Table 1 indicates that for point loads on the
longitudinal centerline at 3X + 3Y the skew bridge deflections are
generally smaller than the corresponding deflections for the
straight and the curved bridge. For a point load on the acute
side of the span the skew bridge deflections at the loaded point
are higher than for the straight or curved bridge, while for
loads o~ the obtuse side of the bridge the opposite is the case.
From the transverse distribution of longitudinal membrane
forces Nx at Section A given in Table 2, it can be observed that
for loads along center girder 3, the membrane forces in the
longitudinal direction in the skew bridge are of the same order
of magnitude as for the straight and the curved bridges. It is
interesting to note that higher membrane forces exist at Section
A on the acute side of the span than on the obtuse side for the
skew bridge. For loads along girder 1 the membrane forces Nx at
Section A are considerably higher for the skew bridge than the
straight or the curved bridge models, while for loads along girder
5 the membrane forces at Section A of the skew bridge are much
smaller than in the corresponding straight and curved bridge cases.
This emphasizes the sensitivity of the skew bridge to transverse
load position as compared to the straight or curved bridge.
Nonlinear analyses were also carried out for the straight and
curved bridge models using NOBOX [A28, B27, B28). For both box
girder bridge models a longitudinal collapse load analysis was
performed first in order to obtain an upper bound on the ultimate
collapse load and an indication of the ultimate collapse mechanism.
A complete nonlinear analysis was then done for both examples and
results were compared with the experimental findings. Finally
eccentric midspan point loads were applied in both cases as
opposed to the more uniform midspan loadings used for the
experimental ultimate load test, and again longitudinal collapse
load analyses and complete nonlinear analyses were performed.
Only typical results for load vs. midspan deflections are
presented in Fig. 21 and 22 for the curved bridge model. Fig. 21
BRIDGE
MODEL
2
3
GIRDERS
SECTION X
4
2
3
GIRDERS
SECTION Y
4
STRAIGHT 5.86 5.96 6.10 5.96 5.86 -1.98 -1.97 -1.97 -1. 97 -1.98
CURVED 5.25 5.61 6.01 6.13 6.31 -1.84 -1.91 -1.97 -2.06 -2.16
6.81 6.27 5.77 4.66 3.58 -.99 -1.48 -1.99 -2.51 -3.03
SKEW
STRAIGHT 4.17 4.94 5.86 6.97 8.30 -2.55 -2.28 -1.99 -1.71 -1.44
CURVED 4.06 5.11 6.31 7.73 9.36 -2.55 -2.37 -2.18 -2.02 -1.87
2.59 3.01 3.58 4.30 5.90 -1.24 -1.12 -1.02 -.91 -.78
SKEW
3X
5X
IIFT~O.305ml
STRAIGHT 8.29 6.97 5.85 4.93 4.17 -1.43 -1.71 -1.99 -2.28 -2.55
CURVED 7.14 6.11 5.25 4.59 4.06 -1.16 -1.49 -1.81 -2.15 -2.50
SKEW 12.46 9.18 6.80 4.60 2.56 -.76 -1.88 -3.02 -4.1B -5.37
IX
STRAIGHT 1.61 2.68 3.88 5.28 6.86 1.36 2.25 3.50 5.30 8.00
5X + 5Y CURVED 1.56 2.78 4.15 5.72 7.49 1.27 2.30 3.69 5.63 8.46
SKEW -2.82 -1.18 .55 2.42 5.12 1.10 ~.16 5.55 8.50 12.33
STRAIGHT 3.86 3.98 4.13 3.98 3.86 3.50 3.89 4.90 3.89 3.50
3X + 3V CURVED 3.44 3.72 4.04 4.07 4.12 3.15 3.67 4.81 3.95 3.71
SKEW
3.78 3.77 3.78 3.16 2.56 2.36 3.06 4.36 3.65 3.55
AX
YD
YD
,<
II I II I;
AX
/1 II I II
~I
STRAIGHT 6.86 5.27 3.88 3.59 1.61 8.00 5.30 3.50 2.25 1.36
1X + lY CURVED 5.79 4.61 3.41 2.40 1.50 4.77 4.77 3.18 2.01 1.33
SKEW 11.69 8.30 5.81 3.99 1.32 5.75 2.63 0.33 -1.48 -3.04
LOAD
CASE
TABLE 1
'-0
UI
5X+ 5Y
3X+ 3Y
1X+ 1Y
GIRDER
12
3
5
7
9
11
WEB
STRAIGHT -23.2 -21.6 -20.7 -22.1 -24.1 -26.2 -28.1 15.6 17.6 21. 6 27.0 32.3
CURVED -24.3'-19.7 -18.0 -20.7 -23.8 -25.9 -37.3 18.1 18.6 21. 5 25.6 29.3
SKEW
-6.5 -5.2 -5.6 -5.3 -12.7 -26.0 -34.2 4.0 5.5 5.9 15.1 28.6
PLATE
STRAIGHT
CURVED
SKEW
STRAIGHT
CURVED
SKEW
WEB
STRAIGHT -21.6 -22.1 -23.0 -23.8 -23.0 -22.1 -21.6 21.7 22.8 23.1 22.8 21.7
CURVED -27.2 -24.2 -20.1 -23.0 -21.8 -23.1 -24.7 22.8 23.3' 23.5 23.0 21.2
SKEW -30.1 -25.4 -24.0 -20.4 -18.7 -18.6 -19.2 26.4 25.7 21. 2 20.1 17.4
10
BOTTOM NODES
1FT :O.305m
PLATE
II
STRAIGHT
CURVED
SKEW
(!)
WEB
TOP NOOES
G>
STRAIGHT -2B.1 -26.2 -24.1 -22.1 -20.7 -21.6 -23.2 32.3 27.0 21.6 17.6 15.6
CURVED 41.6 -32.3 -24.7 -20.8 -18.0 -19.4 -22.5 32.5 26.7 21.8 18.6 17.2
SKEW -42.2 -37.2 -34.2 -32.5 -25.7 -28.0 -33.4 34.8 38.1 34.6 27.6 26.1
<D
81012
PLATE
24
-1:1:1:1:1-
NODE I
LOAD
BRIDGE
CASE LOCATION MODEL
TABLE 2
161
162
(KIP)
180
.... ,.;=g
160
:z
""
CVI
m
~
140
C-
120
c>:
:z
CVI
'"
w 100
C-
80
""-'
0
EXPERIMENT
'"w
NOBOX
( KIP)
176 K
173 K
160
140
120
C-
c>:
0
60
D
W
180
-'
EXPERIMENT
NOBOX
Cl
u.J
-'
-'
C-
CC-
C-
c>:
c(
__~__~--~--~
0.5 1.0 1.5 2.0 2.5
DEFLECTION AT IX (IN.)
O~--~--
FIG. 21
(KJ P)
mDSPAN
LO.~D 1BO
mDSPM
( KIP)
LO.~D 180
160
.~
"
w
140
'-
l!,:..
120
100
LOAD AT
2X,3X,4X
160
140
120
100
LOAD AT
5X
LOAD AT
5X
LOAD AT
lX
20
2
MIDSPAN DEFLECTION AT IX
FIG. 22
00
(IN.)
"1IDSPAN DEFLECTION AT 5X
(I N. )
163
10
0.667'
-1t
20
30
40
: O. 542'
II
50
60
1K1P
0.548',
~}.O'
'>
70 (FT.)
432000 KSF
r,
188000 KSF
9.0,L
t--9.0'
(!)
GIRDER <D
/
'-- (!)
FOR SECTION 7
(a)
SECTION 1
GIRDER
il'
CELL
BOX
(b)
FIG. 23
Girder Deflections
164
-4
-2
CELL
~--- ...
BOX
2
4
O'I~
C'.
>
M
0 CELL
SECTION
GIRDEll ,
(b)
FIG. 24
BOX
165
A large difference exists, however, in the computer time
necessary for the solution by the grillage model of NOBOX as compared
to the finite element model of CELL. CELL requires about six times
the central processor time required by NOBOX.
In addition, it
should be noted that NOBOX yields directly design quantities such
as girder moments and shears, which have to be calculated in an
additional step from the internal nodal forces when CELL is used.
10.4
Example 4
Single Cell, Curved, Two Span, Prestressed Box
Girder Model
166
FIG. 25"
~1 f1+U"11111111
~tmml
.
CI>
8 0 $rI
~
; ..
:n-P"
CD
c i
FIG. 25b
~I~
~ ~
!:l P:\ ~
N
~ r:
rII
S~ ~ ~
.'S r,; i h =,
jII
18~
~ !:!
!!!
l!
~ ~
!!? ~
..
lIS
N
~, ~, ~ , ~, ~ ~ ~
~ r;, ,..
\!l !!! ~ ~ ~
...r-ill
"S 'S,
r;;
~ ~
y.;, .'
(VI
Po
167
r._'7~
(a)
",
'P
(b)
(e)
FIG. 26
168
- - - EXPERII1ENT
-1.0
Go --oi) CURDIP
A,MEGAN
cr (MPa)
(a)
ELEMENT
33
y! / p,,-'p
I
0.4 cr (MPa)
0.3
pI, ""-
0.2
0.1
2224 N
-0.1
-0.2
-0.3
-0.4
(b)
Transverse (Radial) Stresses (MPa) at OUter and Inner Surface of Top Slab
FIG. 27
169
concrete box girder bridge (Fig. 28) using the computer program
PCFRAME [A29, B23] in order to determine its response to time
dependent effects and to increasing truck overloads up to failure.
The cross section (Fig. 28), with a two lane roadway width, is
assumed constant over the entire length of the bridge.
It is
assumed that simple support conditions exist at the ends of the
structure and vertical bearing supports exist at the two interior
supports.
The cast-in-place, post-tensioned bridge is designed using
California highway design criteria. One prestressing (PS) posttensioned tendon in each vertical web is assumed to have the
profile shown in Fig. 28. Two bridges are analyzed for comparison,
one without and one with additional longitudinal mild reinforcing
steel (RS), uniformly distributed over the cross-section. This
steel, not required for strength, is provided for construction
purposes.
For PCFRAME analysis the cross section is discretized into
ten concrete layers, two for each of the top and bottom slabs and
six for the webs. Where included the mild RS is modeled by three
steel layers, one in each of the top and bottom slabs and one at
the midheight of the web.
The time dependent analysis traces the response of the bridge
under dead load only from its initial prestressing at 28 days
through a period of 10,000 days (27 years). Time dependent
variations of deflections, curvatures, moments, and stresses as
well as prestressing forces are traced. Only typical results for
the vertical deflection (Fig. 29a) and the prestressing tendon
stress (Fig. 29b) versus time are presented here.
Notably evident are the smaller deflections for the bridge
with, as compared to without, mild RS (Fig. 29a). The ratios of
the dead load deflection after 10,000 days to the initial deflection
at 28 days are about 3.0 and 4.2 for the two cases. As expected
the prestressing tendon stress decreases with time (Fig. 29b) due
to creep and shrinkage of the concrete and relaxation of the stress
in the prestressing steel. Most of the prestress loss occurs within
the first 100 days. The analysis by PCFRAME predicts considerably
smaller losses than code values, probably because the latter are a
conservative estimate for design.
In order to study the response of the bridge to increasinq
load levels up to ultimate failure, a nonlinear material and time
dependent analysis was performed with PCFRAME in which the bridge
was loaded to failure using increasing levels of truck live loads.
The truck live loads were applied at two different ages after the
initial application of dead load plus prestress, first at an early
age of 30 days and second at a much later age of 10,000 days (27
48.8
I!I
AI_
II
6
1 m
-A.71'T"
48 8
1.1..
J..
10.4
I!I
..
_I
FIG. 28
.ld:.
5.5m
ro'"Pf:~:": )0"
~
I
25 em
l~oor~
..
f'
::;
171
14
E
u
...>3
5o
--'
UJ
I-
c::
...
l-
e<
.L
UJ
:>
::0....
:a:
2
O~----__----~~----~
o~----~----~----~~
10
100
1000
10000
10
TIME IN DAYS
(al
Vertical Displacement
~v
(bl
VS. T~e
'"
I-
u.
10000
IN !lAYS
a vs. Time
Tendon Stress
6
a
0-
""-
1000
Tl~'E
FIG. 29
0-
100
0-
;:;:-5
e<
~ 4
~
o 3
c::
C3
--'
--'
.JL
:x
!lL!
6V
O~~--~--~~--~~
o 0.2 0.4 0.6 0.8 1.0 1.2
DISPLACEt~ENT
(al
BV, em
Vertical Displacement
vs. OVerload P/PO
FIG. 30
6v
12
14
16
TENDON STRESS
(bl
18
If, MPa
Tendon Stress fI
VS. OVerload PlPo
20
172
years) .
The truck live load consisted of an overload vehicle requiring
a special permit to be used on California highways. Total length
of the vehicle is 108 ft (33m), total weight Po= 107 tons (952 kN),
which is carried on 13 axles. In the overload analysis one truck
only is positioned in the middle of the center span of the bridge,
and the structural load vector representing Po is incremented until
ultimate failure occurs. Only the bridge without mild RS was
analyzed since its contribution to the ultimate capacity is small.
In general, overload response of the bridge at the two loading
ages proved to be quite similar. Typical results are presented for
the vertical deflection (Fig. 30a) for both loading ages and for the
prestressing tendon stress (Fig. 30b) loaded only at an age of 30
days, since loading at 10,000 days results in similar curves, but
shifted slightly.
The load-deflection curves (Fig. 30a) show that the structure
is slightly stiffer when loaded at 30 days than when loaded at
10,000 days. This is a result of the reduction in prestressing
force, due to time dependent effects, and the decrease in the
cracking load of the bridge. The large overload capacity of the
bridge is evident with ultimate failure occuring under dead load
plus 7.2 times the live load due to one vehicle.
It can also be
seen in Fig. 30a that the ultimate load does not vary significantly
with age at loading.
Fig. 30b plots tendon stresses at two critical locations
along the girder under increasing vehicle overload. Response is
shown only for loading at the age of 30 days.
The initiation of
cracking at the critical locations is readily identified in this
figure. After dead load plus three times the load of one vehicle,
cracking is initiated at the center of the main span, resulting in
a large increase in tendon stress at that location. After dead load
plus five times the load of one vehicle, cracking is initiated over
the interior support, resulting in a similarly large increase in
tendon stress at that location, and a further increase in tendon
stress at center span due to a redistribution of moment. Tendon
stress continues to increase under increasing overload until ultimate
failure of ,the bridge occurs due to rupture of the tendons at the
center of the main span.
11.0
173
12.0
ACKNOWLEDGEMENTS
174
APPENDIX 1.
175
A.
176
ers. The structures must be simply supported at its two ends.
Each plate element is idealized by a number of longitudinal finite
strips in which the properties of the longitudinal and transverse
stiffeners are distributed uniformly over the area of each strip
and are accounted for in the analysis. The finite strip method
is used to determine the strip stiffness. The displacement
patterns are assumed to vary as harmonics longitudinally.
In the
transverse direction, a linear variation of the in-plane displacements and a cubic variation of the normal displacements are chosen.
A direct stiffness harmonic analysis is used to analyze the
assembled structure.
CURSTR Program (1970)
[A20]
[A26]
177
[A27]
178
coordinate generation options minimize the required input data.
FINPLA Program (1967)
[A12,B6, B16]
179
generations options minimize the required input data.
D.
[A2, B4]
180
before, at and anytime after the transfer of prestress.
SEGAN Program (1978)
[A21, B19]
181
failure and collapse of the structure.
and failure loads can be determined.
Both collapse
mech~nisms
182
APPENDIX 2.
183
Research Reports
1.
2.
3.
4.
Willam, K. J., and Scordelis, A.C., "Analysis of Orthotropic Folded Plates with Eccentric Stiffeners,"
Structural Engineering and Structural Mechanics Report
No. SESM 70-2, University of California, Berkeley,
February 1970 (PB 191 051).
5.
6.
7.
8.
9.
184
Girder Bridge Model, Volumn I, "Structural Engineering and
Structural Mechanics Report No. UC SESM 71-5, University
of California, Berkeley, April 1971 (PB 199 187).
10.
11.
12.
13.
14.
15.
16.
17.
18.
185
"Structural Behavior of a Curved Two Span Reinforced
Concrete Box Girder Bridge Model, Volume II," structural
Engineering and Structural Mechanics Report No. UC SESM
74-6, University of California, Berkeley, September 1974
(PB 242 524/AS).
19.
20.
21.
22.
23.
24.
25.
26.
186
for Curved Prestressed Box Girder Bridges," Structures
and Materials Research Report No. UC SESM 80-10, Department
of Civil Engineering, University of California, Berkeley,
December 1980.
B.
27.
28.
29.
Technical Papers
1.
2.
3.
4.
5.
187
Vienna, September 1970.
6.
7.
8.
9.
10.
11.
12.
13.
14.
188
15.
16.
17.
18.
19.
20.
21.
22.
23.
24.
189
25.
26.
27.
28.
191
S. Tanvir Wasti
Alex C. Scordelis
SUMMARY
Comparisons of theoretical and experimental values of selected
structural quantities, such as reactions, deflections and moments
are presented for three large scale (1:2.8), two span, four cell,
reinforced concrete box girder bridge models tested under dead load,
working loads, and overloads, including loading to failure. All
three models had the same cross-sectional and longitudinal span
dimensions, but in plan one was straight, one was curved and one
was skew. The models were 72 ft (21 m) long along the longitudinal
centerline, 12 ft (3.7 m) wide and 1.71 ft (0.52 m) deep.
Amounts of steel reinforcement used in the straight and curved
model were similar, but in the skew model the amount was considerably less.
General agreement was found between theoretical and experimental results. The structural response of the straight and curved
models to point loads at midspan sections was similar, however that
of the skew bridge was different and highly dependent on the transverse position of the applied point load. All three bridge models
exhibited excellent behavior with high overload capacities before
ultimate failure.
192
1.
INTRODUCTION
2.
193
9.8 sq. in. (6323 mm 2 ) in the bottom slabs and 68 No.4 rebars of
area 13.4 sq. in. (8645 mm 2 ) in the top slab over the center bent.
The reductions represent smaller design moments for skew bridges
because of the observed decrease in dead load moments in skew
bridges [7]. It should, however, be pointed out that the main
longitudinal reinforcement quantities comprise only a part of the
total steel in box girder bridges because of the additional requirements for transverse steel and web reinforcement.
The elevation and plan views of the three bridge models are
given in Fig. 1 with basic dimensions and designations of significant longitudinal and transverse sections. The common transverse
section is shown in Fig. 2. In each model A, B, C and D represent
the instrumented sections, Wand E the West and East abutment
sections and X and Y the midspan sections. Z is the central right
section showing the orientation of the center bent diaphragm for
the straight and curved models, whereas the skew center bent section
along the diaphragm for the skew model is designated as T. Span I
for each of the models possesses a midspan diaphragm which, for the
skew model, is staggered as shown in Fig. 1. In the straight and
curved models it may be observed that the distances between sections
A and X and sections D and Yare 3.58 ft (1.09 m) and that sections
Band C are in each case 3 ft (0.91 m) on either side of section Z.
For the skew model it was decided to divide the sections in mUltiples
of 2.57 ft (0.78 m). In each case it was considered preferable to
locate the instrumented sections A and D close to but not at the
midspan loaded sections.
For all models, locations were designated by the intersection
of the girder and the associated transverse (right or skew) section.
Thus IX refers to the location where girder 1 and the midspan section
X cut each other. For the skew bridge model, the acute side and
obtuse side refer to the bridge portions on either side of girder 3,
these portions are skew-symmetrically placed in the two spans of
the bridge.
3.
EXPERIMENTAL PROGRAM
194
A X
II
WESTI
Y 0
~~~~~~==~==~
SPAN I
SPAN
36.00
A X
~I
..- DIAPHRAGM
1442
358
Y 0
r-
t---o300 30P
1500
358
1500
III
N
1442
r-6
oil
GIRDER
J'
~l
~1543 ~
G>;
.!
2'/\45.
X \" DIAPHRAGM
2S
=---- ~
VY/
ltj
1800
1//
// /
SCHEMATIC
SKEW
BOX
REPRESENTATION
GIRDER
Y D
I
1543
III
//
!@
oil
=0.305 m )
OF ELEVATION
BRIDGE MODELS
~ '"
r-
1800
///
DIMENSIONS IN FT (1 FT
FIG. 1
) }}
CURVED
AND
=,0' - 311ft
'2' - O
M
35/8M
FIG. 2
5' -cf
TYPICAL
SECTION
L, "
" _8 H
ROUND COLUMN
" -6"
OF BOX
... i
GIRDER
2' - 6 114
BRIDGE
MODELS
= 25.4mm
=0.305 m
IN.
FT.
~~i--~--------~--------~----~~~--~--~~~
iOiOI
4 (@ 2'-6 7/8
\0
v.
196
276, 345 and 414 MPa) at sections of maximum positive and negative
moment. After the removal of the conditioning loads, point loads
were applied in each phase in several combinations, to produce
total stresses in the tensile reinforcement of each girder (considered separately) of the order of the working stresses, i.e.,
24 and 30 ksi (166 and 207 MPa). The conditioning loads represented
the successive deterioration of the bridge models under the effects
of overloading.
Details of the very extensive loading program for each model,
incorporating different support restraints, truck and construction
vehicle loads and also loading to failure shall not be dealt with
herein, but may be studied in the above mentioned references.
4.
THEORETICAL ANALYSES
5.
5.1
197
for the straight, curved and skew bridge models. The load cases
treated here will be restricted for lack of space mainly to the
following:
IX, 3X, 5X, IX + lY, 3X + 3Y and 5X + 5Y. The experimental values will comprise results recorded for these point loadings after the 30 ksi (207 MPa) conditioning load level. Because
the analytical models do not allow for cracking, comparisons of
deflections will be limited to experimental results only.
In order to enable an assessment of the response of the three
bridge models to overload stresses, some results of reactions,
deflections and moments at the 60 ksi (414 MPa) stress level will
be included in Section" 6.
5.2
Reactions
198
TABLE'
COMPARISON OF REACTIONS FOR STRAIGHT. CURVED AHO SKEW BRIDGE MODELS
~~~TE
\Mw
'
ME
\:_4iY=fIt2~r'
t
0
D
LOAD
WEST END
E
L
CASE
Mw
TW
39.2
38.0
MF
TF
RE
-355
69.0
96
- 42
-8.2
- :35
-300
69.0
135
- 66
-8.3
-1 <'2
41.2
-283
64.5
110
- 83
-5.7
- 99
41.2
-280
68.1
24
- 78
-7.1
- 96
32.5
193
-193
72.9
117
-196
-5.4
125
-12S
33.6
181
-181
7ol.n
133
-216
-6.0
126
-L.. o
39.3
68.8
86
-S.l
38.1
70.3
135
-8.3
-2
38.8
57
69.1
122
-7.9
- 17
40.0
67
6~.7
103
-7.2
- 20
42.S
-93
93
64.4
91
- 69
-6.9
35
- 35
42.2
-78
78
69.1
107
- 65
-8.1
10
- 10
39.2
335
6~.O
96
42
-8.2
95
37.9
306
69.5
128
72
-8.3
116
36.2
395
74.1
133
82
-10.2
63
36.4
400
76.0
145
74
-9.7
.)5
52.6
-390
390
55.9
65
51
-8.5
-73
73
57.0 -373
373
62.3
91
68
-9.1
-83
83
CU
ST
CU
SK
ST
5X
CU
SK
1 KIP
EAST END
RF
SK
3X
CENTER FOOTING
RW
ST
lX
= 4.448 kN
, FT .0.305 m
ME
1FT-KIP =1.356kN-m
TE
199
TABLE 1 (CONT)
COMPARISON OF REACTIONS
REACTIONS
0
0
LOAD
T
ST
1X
+
1Y
3X
+
3Y
CU
SK
SX
+
SY
1 KIP
END
RW
MW
TW
RF
MF
TF
31. J
-432
133.0
- 84
3' . J
-_ 31
30.1
-374
1 3~. 2
14
-140
30 <1
-<l5 ~~
Re
ME
Te
-332
129.3
-168
35.3
35.5
-381
132.5
-54
-149
33.2
- 346
24.1
264
-264
12H.7
54
-247
47 3
517
-317
24.5
244
-244
132.3
58
-277
47 2
5~S
-5)5
?'.~2
.-I...
31. 2
137.7
31.2
29.7
30.5
30.9
"
138.4
39
138.3
-3
30.8
39
31. 3
37
141.4
72
30.1
35
35.6
-118
118
128.9
-2
35.5
117
-117
36.0
- 41
99
137.6
35.9
85
- 85
31. 0
432
138.0
84
31.0
431
30.3
389
140.2
-2
144
31. 5
468
26.0
458
147.9
166
26.1
4,8
28.7
461
151.4
-31
162
27 .1
<117
47.1
-515
515
128.9
-52
249
24.0
-265
-265
47 .8
-475
475
140.1
-49
300
24.0
-277
!. 77
CU
SK
EAST
T
ST
KIPS)
35.3
T
ST
FT
MODelS
CU
SK
AND
CENTER FOOTING
WEST END
CASE
( KIPS
BRIDGE
= 4448 kN
1 FT
=0.305 m
200
moving towards the acute side of the span. Note that there is no
end moment for the simply supported straight or curved bridges and
that the end moment for the skew bridge is only introduced by the
skewness of the simply supported end abutment.
The torque at the end abutments of the skew bridge has the
same magnitude but opposite sign as the end moment and is smaller
for loads on the acute side but larger for loads on the obtuse side
of the skew span than in corresponding straight or curved bridge
cases.
5.3
Deflections
201
: I
-I
I-~-S
,-3
I~
,-
0.0
90
I
180
I
210
I
45.0
36.0
I
54.0
I
63.0
I
-0.5
1.0
STRAIGHT
CURVED
o SKEW
-05
LOAD AT IV
1.0
2.0
LOAD AT IX .. IV
10
2.0
1 IN ... 25.4 mm
FIG 3
DISTANCE IN FT
GIRDER
W'
R3
3 I"
VERTICAL
1 FT.
= o 3048 m
LOADS AT 1X.W. 1X .. 1Y
72.0
IE
202
II-~-II
3-3
,-
y-'
,,
II-
GIRDER 3
aD
11.0
I
X
/\ 3
DISTANCE IN FT.
18.0
27.0
36.0
45.D
54.0
63.0
72.0
IE
W~I----~I------~I------+------+------+------+------+-----~
-0.5
ao~
__~lO~A~D~A~T~3~X~__________~~____________________~~_
1.0
STRAIGHT
CURVED
o SKEW
2.0
-0.5
ao~~----~~--~~----~~~--------------------~
1.0
2.0
2.0
1 IN . 25.4 mm
FIG. 4
VERTICAl
1 FT. 0.3048 m
203
Ia-I
,-
-I
IY I-,a
I-~'I
a'3
,
I-
3 <<
~, 3
GIRDER S
OD
Wi
DISTANCE IN FT
11.0
I
11.0
I
'n.0
36.0
I
45.0
I
54D
I
63.0
I
720
IE
-S.O
STRAIGHT
CURVED
o SKEW
2.0
-O.S
0.ok~~~2:;=::~.,...---------,
LOAD AT SY
1.0
2.0
1 IN . 2S.4
AG. 5
VERTICAL
mm
1 FT 0.3048 m
204
1.5 times the analytically predicted values.
5.4
Internal Moments
205
TABLE
THEORETICAL
(LOAOS
ANO
EXPERIMENTAL
NORMALIZEO
TOTAL
MOMENTS
AT
1X
3X
5X
1X+1Y
3X+3Y
5X+5Y
INTERNAL
LOAD CASE
A ANO 0
SECTIONS
TO 100 KIPS)
BRIDGE
MODEL
SECTION A
SECTION
THEORY
EXPER.
THEORY
STRAIGHT
565
524
-118
-125
CURVED
557
515
- 97
- 75
SKEW
693
657
-20~
-183
STRAIGHT
561
535
-116
-123
CURVED
569
538
-116
-124
SKEW
561
556
-131
-123
STRAIGHT
565
542
-118
-127
CURVED
582
586
-139
-143
SKEW
420
407
- 58
- 62
STRAIGHT
448
404
443
461
CURVED
458
427
457
431
SKEW
632
597
212
268
STRAIGHT
446
406
448
443
CURVED
453
467
452
399
SKEW
430
433
430
437
STRAIGHT
448
429
443
475
CURVED
444
440
445
467
SKEW
211
273
634
576
1 KIP. 4.448 kN
1 FT -KIP.1.356 kN-m
EXPER
206
6.
6.1
Moments
The point load of about 20 kips (89 kN) at location 5Y has been
taken as a basis of comparison for the assessment of changes in the
transverse distribution of the moments at sections A and D after the
successive application of conditioning loads resulting in nominal
maximum tensile stresses in the steel reinforcement at Sections X
and Y of 24 ksi to 60 ksi (166 MPa to 414 MPa).
In Fig. 7a the
theoretical distribution in percentage of the total section moment
to each girder for a point load at 5Y is given for all three bridge
models. This distribution does not change with an increase in
stress level because of the assumption of linear elastic behavior
in the analytical models. The experimentally measured percentages
in Fig. 7b indicate fluctuati~ns of up to 4% in the girder distributions, but comparisons with theoretical percentage predictions
continue to be reasonable.
....Z
-'
'"
g
....
z
~
Ie
....
FIG.6
14
16
18
1 KIP. 4.44BkN
1 IN .254mm
THEORY
0.3
DEFLECTION
0.2
AT 5Y
0.4
CONDITIONING LOADS
0.5
0.6
FOR A POINT
STRAIGHT
CURVED
SKEW
10D
10
30 KSI
~~
20
10
STRAIGHT
CURVED
SKEW
I
8.0
l f..v/
12
6D
60 KSI ~
12
14
16
18
~O
2~
-.]
208
I-
3-1
1-
1-5
1~~-53
I I II 1~
AX
3~
,~
YD
AX
YD
~s::~
SECT.
A
SECT.
D
(0) THEORY
EI STRAIGHT
oCURVED
fjSKEW
~
SECT
A
...
~'"
~",,,,
1 1
"5!~
1 1 1
00 . . . .
;~~
SECT
(b) EXPERIMENT
FIG. 7
00 ... '"
~~'"
LOADS
W9
7.
CONCLUSIONS
210
S.
REFERENCES
1.
2.
3.
4.
5.
6.
7.
S.
9.
211
10.
11.
213
INTRODUCTION
Design Loading
The design loading employed in major bridges follow the HS
214
20-44 design vehicle (7). For simplified design purposes this loading vehicle is replaced by a lane loading and a concentrated loading. The research reported herein have investigated the practicality of this design vehicle. This vehicle is a pragmatic load.
However, it is known that this loading should be subject to careful
scrutiny in view of the changes in the frequently encountered
vehicular axle combinations and weights (4). The research employed
the design vehicle, rather than the simplified equivalent lane
loadings. Thus, the results reflect the intent of the AASHTO
Specifications more realistically (7).
1.2
Distribution Factor
The research have indicated that for simple span right prestressed concrete I-beam bridges the design parameters that govern
the load distribution are the beam spacing, as it is already
recognized by the current Specifications, span length of the bridge,
and width of the bridge. (Spacing/span length) ratio is critical
in the definition of the effective participatory area of a given
beam. This ratio also facilitates the consideration of the extent
of the single or bi-directional bending of the deck slab. Furthermore, the width of the bridge also needs to be considered. The
215
width can be taken into account either as is, or through the specification of the number of beams and the beam spacing.
It is also
noted that depending upon whether the load distribution formulae
is to be used for old bridges as a rating tool, or for the design
of new bridges, the number of loaded traffic lanes also needs to
be considered as a secondary factor.
In view of the many independent design parameters that need
to be investigated in assessing their full contribution to the
load distribution, the conduct of a test program using full scale
bridges became a highly improbable proposition. However, an analytical research program that can include the contribution of all
design parameters was expeditiously executed.
In the reported
research it has been found that an analytical experiment can be
designed and executed using a finite element software package.
In
excess of 300 right bridges were analyzed having span lengths from
about 40 ft. (12.20 m) to 140 ft. (42.70 m), and with the number
of beams varying from 3 to 15.
Since the computation of the distribution factors require the definition of the mid-span transverse
influence line for the definition of loaded lanes, each bridge had
to be analyzed for ten different loading conditions.
Through the
full computerization of the input and output options of the finite
element analysis, and automatic file transfers from program to
program, it was possible to conduct the extensive parametric
investigation with great ease, after the development of the required
computer programs (1).
All the theoretical analysis results were benchmarked against
available test results of full scale bridges and the tests conducted in the laboratories under controlled conditions for scaled
down bridges (4). An acceptable agreement, even with the relatively
coarse "production tool" version of the finite element mesh, was
observed, the maximum discrepancy being less than 5% (1).
2.1
216
64'-0
.span
Beam Spacing
Design Lanes
FJ-o"
2-12'-0 Wide
1.60-
4-Beam Bridae
5-Beam Brid(le
S/!5.5
cz::
1.20 I-
oz
1= 0.80 I-
be d
abed
abed
abed
abed
;:)
II)
a:
I-
~ 0.40-
Key
Fig. 1
2.2
A
B
B
A
24 fI. Wid. BridQ.
I AI
Diaphragms
217
Certain bridge specifications require that if the span length
is less than some set length then mid-span diaphragms are to be
used. However, if the span length is larger than this set value,
instead of midspan diaphragm, the use of diaphragms at the third
points is specified. The research have showed that for a real
bridge, with an adverse loading condition as would be commonly
encountered in the traffic stream, the diaphragms at midspan lessen
the portion of the load carried by the beam immediately under the
vehicle (1,6). The effects of the diaphragms are imperceptible
for the beams sufficiently away from the vehicle. Furthermore, at
the maximum moment section of the bridge, i.e. mid-span, more
uniform load distribution can be achieved through the use of midspan diaphragms. Third-point diaphragms do not have any justifiable contribution. Therefore, it is recommended that if any
diaphragm is to be used it should be at the mid-span. Also, if
more diaphragms are to be used due to any reason that may be put
forth, then the diaphragms should be located at mid-span and at
quarter points. These conclusions can easily be observed in
Figure 2.
2.3
Skew
218
50.0.-----------------.
40.0
~- 30.0
f!?
z
lIJ
U
Li:
Ib
20.0
o Na Diaphragms
A Diaphragms At L/3
IZ
+ Diaphragms At
~ 10.0
o
L/4
X Diaphragms At L/!5
:E
o Diaphragms At Y2
.0
-5.0 L.....!,..------:!:---~:_--::_-___;!:_-~
Roadway Width
Ec
36'-0"
Fig. 2
J01
J01
Span
Beem Spacing
g
71'-0"
7'-2"
219
where C1 and C2 are constants, sand L are beam spacing and the
span length respectively, and 0 is the skew angle.
The formulae
for the load distribution for prestressed concrete spread boxbeam bridges have a slightly more complex form and the formulae
used for interior and the exterior beams are different (2).
------0-_
OJ2
-a...
0.11
..... .......
0.10
...
Beam C
en 0.09
ffi
0.08
0.07
...80.0&
It
Beam B
Beam A
0.01
45
30
f=o.1
Fig. 3
f-o.
~aO
EI
.
t=0.5
-W- '0.035
220
2.2
20
a::
1.8
u
~
1.6
14
0
I-
t=
::l
In
ii:
IUI
...
1.2
1.0
Nt.-4
Na=7
NL- 3
0.8
06
1/10
1/8 1/1
1f6
1/4
S/L
Fig. 4
221
CONCLUSIONS
Simple span prestressed concrete bridges with in-situ reinforced concrete deck and with prestressed concrete I-beams or
spread box-beams are one of the most commonly used types of bridge
2.2
NB=5
2.0
II:
I-
1.8
1.6
1.4
LL
i=
::l
'!!
1.2
I~
1.0
NB= 6
--;;
NB=8
i"O'
NB=9
NB= 10
.../. )U
0.8
0.6
NB=7
II:
0
-I
1/3
S/L
Fig. 5
222
1.
2.
3.
4.
223
5.
6.
7.
225
Christian Menn
Professor fur Baustatik und Konstruktion
Eidgenossische Technische Hochschule
Zurich
SUMMARY
The applications, typical dimensions, structural behavior,
analysis methods, and ultimate strength considerations of skew slab
bridges are presented.
The effects of cracking and inelastic behavior on the
deflections and servicabili ty of skew slab bridges are reviewed.
The need for prestressing in order to control behavior and provide
servicability is discussed.
Theories of elasticity and plasticity are discussed with
respect to their applicability to the design of skew slab bridges.
It is shown that an inherent stress condition can be superimposed on
the computed ultimate section forces with no effect on ultimate
strength.
Finally, the relationships necessary for
slab bridge for ultimate strength are presented.
designing a
skew
226
1.
INTRODUCTION
227
Fig. 1
FIG. 1:
~.:::::::::::::::::;::::::::::::::::::::::::::
6:::::::::::::::::-::::::::::::::::::::::::::::::::::::1
It)o:c.to::f~AJ
1.0.0.:0.0:.0]
Fig. 2
FIG. 2:
[m]
40
30
20}:::===========
<p = 30'
<p = 60'
<p = 90'
10
Fig. 3
FIG. 3:
228
m2
-------
- - - - me
y
6~O~
lo"b = 0.63
J0.3;;J
t/m2
~x
mt
a) narrow plate
lo/b= 1.67
.,
b) wide plate
narrow
plate
m1
ml
0.42l
()(
wide
plote
m1
ml
()(
1.000
0.164
290
1.060
0.016
22.0
1048
-0.104
18.3
-0.096
16.8
0.448
-0480
51.0
1.088
0.272
-0.304
475
ql~18
ql:/8
[0]
ql~/8
ql~/8
[0]
c) points of design
Fig 4
FIG. 4:
229
Consequently, l4/h3 is a deflection index which gi ves an
indication of the extent of the deflection. An 80 cm plate with a
span of 20 m exhibits a deflection value of
2.
Theory of Plasticity
for
ultimate
strength
is
normally
formulated
as
<
associated safety factors:
dead load
live load
characteristic concrete strength
characteristic steel strength
'14
F IG . 5 :
114
AND
OF DESIGN
TIC POINTS
IS
R
TE
C
A
R
A
CH
ENT ZONES
REINFORCEM
LI2
--------
_ _ -1. ~
b /6
<::)
IN
tv
231
iN
dx
stress
conditions
or
avy
avx
for plates: - - +
dx
ay
Such an inherent
or
are
d2 M
dx2
32 mx
defined
in
principle
as
a2m
+ 2.~ +
ax dy
ax2
a2 m
from
support
232
a)
support displacement
I I
I I I
J1
.Jf
b)
I I qo
Fig 6
FIG 6:
233
where F d,E
or:
where F d,i
2.3
redundant
moments
represent
part
of
the
entire
234
i(Z----------#~B~7
R
~
L_._B
A
h4=f-'-'-'-'-'-'-f
'I
l = 2b
.............
_----
my
m x y .
Section
A- A
Section
Fig 7
_e_e_e_e_.
FIG 7:
Fd
B-B
235
prestressing moment
where IDap
mzp
~:
- p'e
represents the redundant moment due to prestress.
The value chosen for the width over which a prestressing tendon
is considered effective is not of great importance. More important
is the necessity for using the same width (for example b p '" 2' h)
throughout the design calculation.
In the formula for the design ultimate forces, the redundant
moments due to prestressing (an inherent stress condition) are
generally included with a factor of 1.0 (Fig. 8). Taking into
account a possible moment redistribution, however, this factor may
actually vary from about 0.8 to 1.2.
3.
for
narrow
simply
236
+=~l==~~~';;;;;";';;;;;========l+
L = 2b
-----
............ _.-
-.- -
Section
8- 8
~
-'-'
Section
A- A
Fig 8
----1.8 m
g+p
_._-_._. 1.8 m
g+p
FIG. B:
+ 1.0 m
zp
237
A-
...A/~=============/47
k ' - / - - - - - - - - - - - - - - - -...// /
/
/
-c
_I
A-A
B-B
F'-'- _._._._._._.-. .,
Hypothetical arrangement
of prestressing for the
calculation of mzp '
h
AlBIC
Fig 9
FIG. 9:
238
chosen reinforcement directions.
W
x'
y'
xy
W
11
the
ultimate
internal
m~
m*
11
*
2 COl + nyocos
*
2 ~1
mxosin
m* = 21 (1Il_* i;;ll
:t
0
~)
~y sin
2 oC0 1
osin2<Ol +
*
* - mi;;)O(m
ull - mil)
~~
> 0
>0
Negative moments:
> 0
239
FIG 10:
243
~etin Y~lmaz
INTRODUCTION
244
(2)
where
p
p(t) = density of the
u = u(X, Y, Z) = (ul, u2'
velocity of the fluid
Fi= body force at a point
0ij= dynamic stress in the
fluid
u 3)
at a point (X, Y, Z)
(X, Y, Z) of the fluid
fluid at a point (X, Y, Z)
245
Finally,
D( ) _ ~+
at
1ft -
a()
ax.-
uj
(3)
au.
ClO
~J
( 4)
Pat=~+Fi
J
-pO ..
~J
21l(e ..
~J
- 3" ekko ij )
(5)
KEkko ..
q
2G(E .. - -3 Ekko .. )
~J
~J
(6)
KE:kkO ij
(7)
where h
-+
h(X,Y,Z,t)
246
cf(V h
curlcurl h)
(11)
where
(12)
is the speed of sound in the fluid.
Equations (10) and (11) are both forms of the wave equation
for which various methods of solutipn are well established. Equation
(10) is a scalar equation in one unknown, p, whose solution yields
the dynamic pressure field in the fluid domain. Equation (11) is a
vector equation in one, two, or three unknowns governing the timevarying displacements of the fluid medium.
Recall that equations (10) and (11) have been derived under
the assumption that the acoustic fluid is compressible. This implies
that the speed of sound in the fluid, cf, is finite. If the fluid
is assumed to be incompressible, which is a reasonable assumption
for water, the speed of sound becomes infinitely large. In this
case, the response of the fluid is governed by either one of the
following equations
o
V2 h +
curlcurl
(13)
h=
( 14)
247
3.
I'"
r--.
Fluid
~IC
Fig. 1
Consider a single-degree-of-freedom mass-spring system interacting with an acoustic fluid medium as shown in Fig. 1. The hydrodynamic pressure acting on the face of the mass is denoted by p.
Assuming that this face has a unit area p is also equal to the
hydrodynamic force acting on the mass. The displacement of the mass
m from its position at rest is denoted by x. The spring constant is
k. The mass is assumed to be at a distance L from the origin of the
y axis fixed in space.
The equation of motion governing the displacement x(t) of this
idealized structure is
mX
kx
(15)
= -p
a2 q,
ay
1"
(16)
;r=Z<P
cf
-p~
at y = L,
(17)
248
that the velocities of the mass and the fluid particles at the
interface must be equal so that no gap can occur between the two.
This condition is stated mathematically as follows :
at y = L
(18)
The assumption that displacements are small has been used in equation
(18). From equation (16) it is seen that we need two boundary
conditions on in y. Equation (18) provides one of them. The second
one is obtained from the fact that, for propagating waves, the
radiation condition must be satisfied at y = 00. This implies that
as y + 00, must go to zero. The definition of the problem will be
completed by assigning the initial conditions which are taken as
x
0,
= V,
0,
at
0,
(19)
in which V is a constant.
The problem is now well-defined and the solution can be obtained
by various means. The Laplace Transform method used in the present
work yields the following
x=
~ e -a.t SinQt
in which
Q = w
i _a~
w=1f
m'
(20)
(21a)
(2lb)
(2lc)
249
in vacuo is the following: According to Chopra (1981), damping
ratios less than or equal to 0.2 constitute a range which includes
most structures in vacuo. Thus, damping effects on the frequency and
period of such structures are negligible. On the other hand, the
same can not be said for a mass-spring system in a fluid medium.
The a./w for tne interacting system is rIOt necessarily less than
0.2. It is commonly larger than 0.2 and it may even be larger than
1.0. Recall that, for p=lOOO kg/m J , c f =1500 m/sec and w=200 rad/sec,
a.
pc f
3750
m
ThU1l, for a mass m > 3750 kg the a./w ratio will be less than 1.0.
For smaller masses it will be even larger than 1.0. The gE'neralization stated by Chopra (1981) is, therefore, not applicable for
damping caused by the presence of the fluid. The difference between
structural damping and fluid damping can be explained as follows :
Recall that
t,;
and
= 2rrw
a.
w
(22)
which are the damping ratios for a damped structure in vacuo and a
structurally undamped structure in fluid. t,; depends on-cne-ratio
between c and m wnich are both related to the properties of~
same structure. On tne other hand, a./w depends on the ratio between
the quantities which are related to the structure and the surrounding
fluid. Thus, it is, in principle, possible to Come up with a./w
ratios which are much different from corresponding ~ ratios. In case
a./w turns out to be greater than 1.0, tne solution becomes
x
witn
IT
/a.~
L
W
(23)
_ 1
mx +
kx
(24)
iF;
1
cf
=0.
(25)
250
I
(26)
y = L,
y = L.
(27)
rnX + pCf~ +
kx = -rnx B - pcfx B
LS
(21:3)
Wt
(29)
irlt
(30)
substituting equations (29) and (30) into (28) we obtain the real
part of x and it is
2 2
pC f 2
2
P cf
rl) --2] Cos~lt + liJf W Sinrlt}
m
(31)
Real x
2
t2
[ (w
Here,
W =
Nm as
~
m
before.
Nmv cO',ldder tne r('3por.se of a single-degree-of-freedom massspring-damper system in vacuo subjected to a base motion. The
governing equation is (Chopra, 1981)
mx
cx
kx
-mxB
(32)
A comparison of equations (28) and (32) shows that, in fluidstructure interaction problems with base motion, base velocity x
enters the equation in addition to the base acceleration xB. Thi~
is another important difference between the responses of a structure
in vacuo and of that in a fluid medium. The damping force proportional
to tne base velocity is added to the damping force proportional to
tne relative velocity of the mass. Having Ui.sc1lssed seme significant
aspects of fluid-structure interaction using a single-degree-offreedom system, we can now consider a slightly more complicated and
more realistic model which is studied below.
251
4.
d2 <Jl
1 Cl<Jl
Cl 2 <Jl
r or
~8-L
-+--,,-+-,,-~
or 2
0 -
1 Cl <Jl
------.,.~ L
c 2 at
(33)
- - t -__ Vg
Fig. 2
subm~
long rigid
cylinder
at
r = a,
(34a)
252
3 = 0
3r
as
r -+
(34b)
00.
-iwt
(35)
From (34a) :
3
ar
=
v g CosO e
-iwt
at
(36)
Accordingly,
(r,O,t) = l(r) CosO e
-iwt
(37)
+ .!:.
+~
CP' _ 1 <p
r12l
r
21
=- 0
(38)
cf
(39)
in wr.ich .\ is a constant to be detE'cn:ined fro1!\ ("I.e boundary condition at r = a and Hfl) is the Hankel function of the first kind
and fi~5t order. The general solution of (33) is, therefore,
(1)
<P(r,8,t)=A Hi
Wr
(-) CosO
-iwt
(40)
cf
The boundary condition (36) yields
A = v / {~
g dr
H(l) (wr)}
1
cf
(41)
r = a
- ipw{}
a is
r=a
}t;;~-wa/c
f
- - - - - - - - - " - - l.:os8 ,
d
(1)
{d~ Hi
in which
S=
wr/c f
(t;;)}'=wa/c
'"
( 43)
253
= H(l) (s) -.!:.. H(l) (0
l o s 1
~ H(l) (0
ds
( 44)
cf
{Ho
(1)
(wa/c f ) - wa HI
------ Case
( 45)
(wa/c f )}
Call the terms on the right hand side of (45) S, except for the
Case term. Thus,
( 46)
Pa = S Cos8
H is
P
T
4aHSf
11/Z
a
= TIaHS
cosZe de
( 47)
Finally,
iTIaHpc
{H(l) (wa/c f
a
~
Hl(l) (wa/c f )}
wa
( 48)
cf
( 49)
H (1) (wa)
a
c
f
For
exp{i(~ - ~) }
4
cf
-iTT/Z
-1
(50)
254
PT = naHpc V
f
(51)
Assuming that the highest frequency that is encountered in earthquakes is about 10 cps; thus, for (19) to be valid, one must have
2nfa/c f
Tnis yields,
cf
a 21Tf
or
24 m.
This yields the added mass results because (na 2Hp) is the mass of
the fluid displaced by the cylinder. For this assumption,
(53)
l.n 'AIhich
H( 1 ) (wa/c )
1
(54)
H~l) (z)
= J)z)
i Y)z)
(55)
255
we obtain
J 1 -I- i Yl
/\ =
------,-
(J
(56)
iY ) - ~(J
o
z 1
l.n which z = wa/c f and the argument of t,le Bessel functions has
been d:opped. Equation (56) can be brought to the following form
Jl
A=
i Y1
(57)
-----,-----------,,-
(J
~
z
J \
l'
iCY
~
z
Y )
+
y
is
(5tl)
in which
(59a)
(59b)
(59c)
PT
iTIaHPC f Vg
(a
+
y
is)
(60)
PT = TIaHpc f vg(Sinwt
Let
~ =
TIaHpc
=
iCoswt)
(ex
+y
is)
(61)
v. Thus
(62)
Finally,
(63)
in which
.
= -~( O,Sl.nwt
y
-(O,Coswt
y
- SCoswt) ,
(64a)
(64b)
SSinwt) .
l'londimensionalize PT as follows :
256
na Hpg
(65)
'
which
.
-to ( aS~nwt
(67 a)
SCoswt),
XY
to
-- (aCoswt + SSinwt) .
XY
(67b)
i Re(PT) i
X~{aSin(Arctan
iIm(PT)i =
~{aCos(Arctan(~ + SSin(Arctan(~}
(- F
- SCos(Arctan(-
B-} ,
(68)
(69)
A computer program has been coded for the calculation of the above
results. The numerical results obtained from the approximate solution (52) and the exact solution (68) are compared in Fig. 3. It is
seen that, for the range of frequencies realistic for earthquakes,
the low frequency (or the added mass) approximation yields results
very close to the exact solution.
5.
The effect of the surrounding fluid on the response of a singledegree-of-freedom spring-mass system is studied analytically. An
analytical method of solution is used also for the steady state
response analysis of a submerged rigid cylinder. The flexibility
of the cylinder has not been considered to keep the analytical
solution feasible. The solution of more realistic problems, in
wnich, for instance, flexibility of the pier is considered, generally
requires numerical techniques which will be the topic of Part II of
the present paper.
Acknowledgement: This work has been supported by Yliksel Proje
architect-engineers for the Karakaya Dam Bridge project. The
authors are indebted to Mr. Erdal Erkakan, project manager.
A.~.,
257
~
~
15
"'.>'"
.~
.3
lit
cb
'"
tl
0
)(
""
;I
"
~
~
...J
'"
""'
~
.....
""
'!I.!I.
"0
0
c:
.~
Ill."""
258
REFERENCES
1. Abramowitz, M., and Stegun, I.A., (editors), Handbook of
Mathematical Functions, Dover Publ., New York, 1965.
2. Chopra, A., Dynamics of Structures, A Primer, EERI, California,
1981.
3. Lamb, H., Hydrodynamics, Dover publ., 6th edn., New York, 1945.
259
l~UMERICAL
i~uri
Akka
1.
INTRODUCTION
260
approximation to the wave equation using pressure as the dependent
field variable. This meti10d was used to compute the natural
frequencies of submerged structures in finite fluid regions and in
unbounded fluid regions in which the fluid-structure interaction
effects are negligible beyond finite distances from tne structure.
This method was extended by Hunt (1974) and Zarda (1976) to applications in wi1ich the far field effects are important by a procedure
which matches finite element solutions at finite external fluid
boundaries with known analytical solutions for infinite fluid Cases,
A similar formulation is presented in (Nastran, 1976) and implemented
in the finite element, structural analysis computer program, NASTRAN.
In this Case tne finite element equations were derived from variational metnods using a form of the wave equation in which the fluid
pressure is the dependent field variable. Pressure analog methods
nave also been developed which use either "lumped" (Everstine et
al., 197) or consistent (Schroeder and Marcus, 1974,1976) formulations equivalent to the previous finite element pressure methods.
The pressure analog methods are derived from an analogy between
fluid pressure and one of the Cartesian displacement components in
finite element formulation of classic linear elasticity. These
methods have been implemented using readily available finite element,
structural analysis programs such as iiastran.
Further investigations have produced structurally compatible,
acoustic fluid finite elements for linear two-dimensional analysis
based on pressure as the dependent variable (Nahavandi et al. 1975;
Marcus, 1973). Bathe and Hahn (1979) have implemented two-and
three-dimensional, acoustic fluid elements in the structural
analysis program ADIHA which is capable of nonlinear and large
deformation analyses. These methods require the addition of separate
acoustic fluid finite elements to the element libraries of existing
finite element, structural analysis programs.
In all of tne preceding finite element methods, the boundary
conditions at the fluid-structure interfaces are satisfied such that
the acoustic-fluid finite elements are compatible with structural
elements. The coupling of the fluid-structure interaction is
achieved directly throug;1 the assembly of the system mass, damping,
and stiffness matrices. The generality and versatility of the
finite element method in purely structural applications is thus
extended to the acoustic fluid-structure application with relative
ease.
This present work represents an extension of techniques
developed in (Akay et al., 1977; Akka et al., 1979) for the
implementation of two-dimensional and axisy~tric, acoustic fluid
finite elements in the structural analysis, finite element program
SAPIV. The finite element formulation is derived from the application of the method of weighted residuals using Galerkins criterion
to a form of the wave equation based on fluid displacement as the
independent field variable. Implementation of acoustic fluid finite
261
pp -
->-
KIJ
->-
IJp =
(1)
We define the sum of the terms on the left hand side of equation
(1) to be the gene ral ized force,
J
->-
->-
= pp - KV Vp
(2)
262
-+
-+
-+--+
-+
( 4)
'lp dxdydz
-+
-+
(5)
-+
-+
-+
Ids KopV'p . n ds
where tne vector ri is tne outward unit normal to tile surface s. Tne
first and second integrals in equation (5) may be rewritten using
tne ident it ies
(3t)
dP 2
..
dp
-+
-+
-+-+
= 2'lp
.V'(op)
(7)
1
2
o{/~(2 P (p)
1-+
+ 2 KV'p
-
-+
V'p) dxdydz}
-+
-+
fdS Kop(V'p).n ds
(d)
(9)
263
-;.
-;.
( 10)
by tne equation
2
d V _
dP dp. 1- J
K .
1-J
( 11)
A ma5S matrix 1m. 1 1-S formed from tile internal kinetic energy
1-J
T = f rI
2
p(p) dxdydz
(12 )
by the equation
2
m .
1-J
(13)
( 14)
acl.
f.
1-
acl.
1-,,, + __
1-N +
f "sKN(
a
i ax I X
Cly y
aN.
az1-N)d
Z
S
(15)
264
K'V P -
~
dt
;0
L(p)
;0
(16)
->-
Ii '"
m
l: NiPi
i;ol
(18)
(20)
265
51
(e)
N. (e)
1
a2 (e)
a2 (e)
ax
ay
'i
a2 (e)
(e)
at
dxdydz Ce ) = 0,
i = 1,2, .. ,r
(21)
wnere r is the number of nodes for the given element. We now write
(22)
Upon substitution of equation (22) into equation (21) and integrating
the Laplacian terms by parts we obtain (dropping the (e) superscripts)
t
ClN.
K(_1
aN.
"N
I~I +_1
ax ax
ay
eM.
,,-.T
I~I
ay
(23)
k .. =
1J
I"
f. =
1
51
"
p N .i~. dxdydz
(25)
1 J
a,L
ai~.
dX
ax
1 _J
KC_
al~.
aN. aN.
_1
ay
aN.
_J
dy
aN. aN.
_1
dZ
_J) d d
az
x ydz
(26)
aN.
(27)
266
3.
FORMULATIOi~
-+
-+
h = s
on s
(28)
V.h
(29)
-+
h Cx,y,z,t)
y
IN (x,y ,z)h
h/x,y,z,t)
ILL (x ,y ,z) h
1 l
0
yl
0
(t)
Zl
t)
r)
J
(30)
~quations
Cl 2h
,?h
Cl 2 h
Cl 2 h
Clx
ClxCly
ClxClz'
Clt
a2 h
Cl 2 h
Cly
CiyCiz
j{C--;' + __
x + __x)_p--;, =
Cl 2 h
K(--y
ClxCly
+ ~+ --y)
'1
Cl 2 n
Cl 2
+ __
z +
dXClz
ClyClz
K( __z_
,lh
(31)
-p----f = 0
(32)
Cl 2 h
Cit
Ci 2h
---T-)
-P---T=0
Ciz
Cit
C33)
267
Applying the method of weighted residuals using Galerkin IS criterion
to equation (31) yields (dropping the (e) superscripts)
,?h
3Zh
32h
:lh
3x3y
dXdZ
IS-Pi IK(~ + __
x + __
X)_ p~1
3x
dxdydz
(34)
dt
dZh
__
x) dxdyd z
dXdZ
d
dh x
= I nl~ ~ (-" ~ dX
Clx
ah
ah
+ ---..:. + ~)
ely
(35)
dxdydz
dZ
We now subst i tute for ilx from equat ions (30) and inte grate
(35) by parts to obtain
IN. ( I ClN I +
dS
dX
I dN I
dy
_ I
I dN I){h }n ds
dZ
x x
~.
r;:
_...2( I 31~
dX
dX
"
+ I dN
dy
+ I dl~
3z
I ) {h } dxdydz
x
( 36)
IN. (I eM
dS
dX
+ I Cli~ I + I ,h~
dy
- I
3H.
r;:
I ) {h }n
3z
"'"
ds
'"
'"
~
I~I
+ I~I + I~I
dy
dX
dy
dZ
) {h }dxdydz
y
(37)
and
I",
oS
d. ( I ~N I
~
oX
+ I oy
~N I + I aZ
~i~ I
- I
3N.
r;:
'"
) {h}n
Z
-~
(I~I
+
dZ
dX
ds
IClNI +
dy
I~Nzl
a
){h }dxdydz
Z
(38)
268
f as
~ K.t~.L
3N
I"I
aX
3N I + I,,)
3N I {h}.n
+ I"oy
oZ
ds
3N;
3N;
dN;
I ~N I
~N
~N
dX
3y
dZ
dX
3y
dZ
){h}
(39)
( 41)
~ N. {p }ds
as
L
n
IK LJ.. 1
fll
wuere
IB L1
3N.
L
3X
0
31~
.
L
dY
3i~
.
L
3N.
L
dxdydz
( 42)
o1
( 43)
0
3N.
L
8Z
3Y
8X
az
3H.
L
3N.
L
3Y
ax
3N.
L
8Z
3N.
L
269
and
IC I
( 44)
Defining
1m.l.J1
= f~p
( 45)
N.N. dxdydz
l. J
l.J
{f.}
(46)
l.
Icl
(l+\!) (l-2v)
I-\!
\!
I-v
\!
I-v
(1-2v) /2
(1-2v)/2
Sym.
( 47)
(1-2v) /2
270
4.
k~
APPLICATION
k~D
CONCLUSIONS
271
fOr the present problem. It was also decided tnat the Taft record
of 1952 and the EI-Centro record of 1940 modified with respect to
tne appropriate peak acceleration should be employed in the studies
of the superstructure and its interaction with the fluid medium.
Tne use of two different earthquake records allows one to study the
effect of the frequency content of an earthquake on the response of
the system. In Fig. 5 the north-south component of the EI-Centro
record of 1940 which is modified for 0.16g is given. The Taft record
~s not snown but it is also modified for 0.16g.
Tne maximum responses of the pier model subjected to the
modified Taft earthquake and to the modified EI-Centro earthquake
are illustrated in Figs. 6a and 6b, respectively. In these figures,
maximum displacement, maximum bending stress and maximum shear
stress are plotted along the length of the pier for both the in
vacuo and submerged cases. The effect of the fluid interaction-is
clearly seen in these figures. For all the cases considered, the
presence of tne fluid amplifies the response of the structure. The
maximum responses resulting from the modified EI-Centro record are
generally higher than those resulting from the modified Taft record.
The displacement at the top and the mOment at the base of the
pier subjected to the modified EI-Centro record are given in Figs.
7 and 0, respectively, as functions of time. The corresponding
responses for the modified Taft record are given in Fig. 9 and 10.
The amplifying effect of the presence of water is clearly noted in
these figures. Respective comparisons of Figs. 7 and 9 and Figs. 8
and 10 reveal the effect of the frequency content of the earthquake
record on the structural response.
The numerical results given in this section are for a specific
pier model only. The results indi~ate that the presence of water
interacting with the structure affects the response of the latter
significantly. The maximum responses are, in sOme cases, almost
doubled. However, it should be emphasized that these conclusions
are valid for the specific pier considered only. For another pier
with different cross sectional dimensions and/or height the effect
of water may turn out to be not so much significant. Thus it is
suggested that each case be studied on its own and the conclusions
tnereof should not be generalized.
ACKNOWLEDGMENTS
The work described in this paper was performed in conjunction
witn consultancy services provided to Yliksel Proje A.$., architectengineers for the project. The encouragement of Mr. Erdal Erktakan,
Project Manager, is appreciated. The sections on the finite element
formulation are essentially based on a report by Mr. J.G. Bennett of
tne University of Michigan submitted to N .Akka as part of the
requirements for a graduate course.
272
6i7.719
~
.:1.619
Sl.....-
IL
..........
698.084
~
~.284
~!I.OOO
693.000
.~
~.
A-A
~.OOO
I. ' I f
'f/lfff
Ifill
C - C
FIG. 1 ELEVATION AND CROSS SECTION OF THE PIER
273
t-.
ACTUAL
EQUIVALENT
:sz
..
WATER LEVER
I:
I:
I
II
II
III
II
II
II
II
II
II
11
II
-1
I
I
I
I
I
--I
44(
EQUIVALENT PIER
I
I
I
I
I
I
I
I
I
I
I
I
I
m.
z
40
37
34
31
28
25
22
19
16
13
10
13
12
11
10
41
38
35
32
29
26
23
20
17
14
11
Is
:;
42
39
36
33
30
gr
24
21
18
B.
26
25
24
23
22
21
20
19
18
17
16
15
14
5. 3 m.
~
L
51 m.(16732 3 It)
..'\7
'"
FIG. 3
ACTUAL PIER
II
II 1"II
II
I"
II
II
II
II
II
" I"II
II
II j II
II
II
II
II
II
II
II III
II
II
II
II
II
II
II
II
II
II
II
II
II
II
II
II
II
11
8 It )
+>
-.)
40
37
34
28
..!1!
..!.
..!Q
...1
--!
13
12
11
10
11
21
24
33
30
.!!
J1
....
41
26
!I
r
r
I
,~
MOOE 1
r
,,
MODE 2
0>-- -
IN VACUO
M00E3
-0
FIG. 4
1
2
CASE
FREQUENCY (CPS)
MODE 1 MOllE 2
MODE 3
10.03 T :51.401 118.40
7.98
40.28
91.:54
0.10 020
WATER LEVEL
MIlE
OUTSIDE
ZERO
I ZERO
_MAL
NORMAL
42~__~4=~========~==~~~~F=~~========~==~==~~F=~-=====__
3825 ~
3$
23
129
22
26
23
32
19
20
20
17
18
17
, 14
15
16
14
....,
IV
vo
276
....
o
o
a:
o
~
a:
40l
37
~l
31
28!
2S!
22~
19!
16;
~!
10~
7~
4!
13
12
11
10
19
21
./
42~,
FIG. 60
3e
J11
41
1\'
"
21 2 42 4 848 1060
BENDING STRESS (k~/e",)
"I
.........
IN VACUO
WITH WATER
1.40 :1 80
g.60
SHEAR STRESS (k~/em')
~-----o
0----00
II,
II
24
II
3O
27
24
/
/
.
.
:
.I
~
I
-I!
21JI
18 II
IS
12
..j
..!
38 25
3:1
32
23
29
22
26
23
18'
17
20
20
17
14
11
16
15
141
1 :I
tv
...,
...,
4' ~
3[
3~
3!
2!
2~
16!
.!Q!
13
12
11
10
26
25
24
23
22
21
20
19
18
17
16
gil
---I
42
!!
24
3!
18
~.IO
: :
7.6~ 10.2
DISPLACEMENT ( cm )
2.~~
iI
~JI
I
I
~J r
225
675
1I.2~
8ENOING STRESS ( kg/cm~ )
0---0
0 - _ _ -0
IN VACUO
WITH WATER
3.77
11.31
SHEAR STRESS (kg/c",)
.,
II
'I
41
38
32
29
26
23
20
17
14
11
1 ~ 141
15
2
4~
00
--.J
- 9
- 6
- 3
31
i:1
'
FIG. 7
t .,
,
:~
:Ii'1':i1'l
II
"II'
" I',I
"
,III
(, ,t
em.
[,
WITH WATER
(%
5 DAMPING)
IN VACUO
hrne(sec)
DISPLACEMENT
tv
.....
\0
30.000\
- 5Q<XXj
_40000\
_ 20.000 \
_ 10.,
100001
20.000\
300001
40 000\
50=\
FIG. 8
"
I.
,,
{%
5 DAMPING}
IN VACUO
WITH WA1'ER
MOMENT Tm.
_ 6
- 3
'3
f'I
I
I
I.
"
;~
I
I'
em,
FIG, 9
:,
',I
~~
"
I',
\I
TAFT
",'
,>;11'
i'
, ,
RECORD
l >il\
'
::
! \ :: \
:\
'.
Ii
'
'
!.: \:if
' ,I,
' I
"
,
hi
::\1
'
"
,!'
.I,
'I'
(0/0 5 DAMPINGl
'/11TH WATER
IN VACUO
' '
i~
'
'
>,'
I 11~I\i, ' ,
",
',j'
;1,
, ,
'llme(s
: V' I
'~\
"
"
' \ , ~ ,{
I: , :
____ -
"!
\.
,! ,: I
,:".', i
>"
"
',21
I
!,
__
~'
\~"
"
:' \
I'
! 'I
(i.
I:
eel
00
'
- :i0.000
- 40.000
- 30.000
-20.000
-Io,ooof
10.000
20000
\!
'
\11
'II
~
I
II
I I
rill)'"
I
if'J
It
il
I'
,I
I,"
I
I,
\1
V:: :1~ \Y
,i "
I' I: ::
'I
1t-.;lt{III::'~
:11:
. I' , If. ,
I
I I I I i Iill, I"'
I
I
II I 'i~'i\
II!
, " L
I,,,
'\'
I
I'PII
,! ' I
1 "
I \
IiI
II
'!
(0/0
~ OAMPING I
IN VACUO
WliH WATER
I'
If
rV ~ I
,Ill I
II' " I I
II
J,
I
"
, ' I II I
I '/
~
I
\~\"
II I, ~
11111",
1:1 " II II JI
, '
t
I II II /1" I
,
~
II /1 II
I
I
!
f I 11
I I ' \ /1
I
J
D
MOMENT AT BASE OF TO PIER SUBJECTED TO THE MODIFIE
RECORD
TAFT
:
I
II
I I
I , II II
I ,\ ,
' Ii I
I,
:'/1
jl"11
'I'
,
'i'!
1\
::I!I~:
\,
Tm.
'II!~I"
'II
FIG. 10
';
: 1\
MOMENT
l
oo,oooL,j,~~~'\1
30,000
40,000
50000
I
\
\}
'~i
,I
l\7Ttime (SeCI
'.
t..>
co
t..>
283
REFERENCES
1.
2.
3.
Bathe, K.J., and Hahn, W.F., "On Transient Analysis of FluidStructure Systems ," Computers and Structures, Vol. 10, 1979,
pp. 303-391.
4.
5.
Bedross ian, Bedros and Di Maggio, Frank L., "Acoust ic Approximations in Fluid-Shell Interactions ," Journal Engr. Mech. Div.,
A.S.C.E., Vol. 93, No. EM3, June, 1972, pp. 731-743.
o.
7.
d.
9.
Glilkan, P., Erdik, M., Akka, N., Y~lmaz, <;;., and Qner, M.,
"Seismic Analysis of a Multispan Railway Bridge: A Case Study,"
Proc. 7th European Conference on Earthquake Engineering,
September 20-25,1982, Athens, Greece. Vol. 6, pp. 163-170.
10.
284
11.
12.
13.
Marcus, M. S., "A Finite Elemont Method Appl ied to ti1e Vibration
of Subm2rged Plates," Journal of Ship Research, Vol. 22, No.2,
June 1978, pp. 94-99.
14.
15.
Aeronaut~cs
10.
17.
lS.
19.
285
E. ZELLER
Associate Professor
Department of Civil Engineering - Institut Universitaire
de Technologie - Universite STRASBOURG-III (FRANCE)
INTRODUCTION
Bridges have been built since ancient times. As
crossing and linking needs increased and developed,bridges have
multiplied and improved.
In the beginning, the know-how was closely related to
previous experience. From the point-of-view of materials used, as
well as the knowledge of structural behaviour and an understanding
of building techniques, bridges sum up the science and know-how
of an era.
If we observe the evolution of the most important
structures in history, we will be able to follow up the series of
progressive innovations and improvements. Moreover, it can clearly
be seen that those bridges nowadays "old", were modern and even
ahead of their time since the most advanced scientific and
technical elements had to be worked out and then assembled.
Nowadays, industry regularly provides high specification
materials, and present technology supplies powerful means for a
quick setting up and construction. But undoubtedly the most
outstanding progress in the last twenty years is due to the
elaboration of mathematical modelling of bridges, thus allowing
be,..tter analysis.
If the hand solution procedures have permitted predicting fairly accurately the deflection and the stresses of
structures understatic or quite static loads, when it comes to
deal with complex structures or carry out a dynamic analysis, the
286
Chauss~es,
LCPC,
287
288
v.
CONCLUSIONS
I. METHOD OF MEASUREMENT AND ANALYSIS
PROCES~
289
Excitation Method
J F (t)
A---
A [tr-
A>
rp
LL__
:QAtLet go Load
-.
~
Vehicle
+ Obstacle
\:0
Shaki ng machi ne
~QQQ---
~~~
-,
-----
tt
-- ----
.t
'~
~
h .
~~:
,
14fV -
4A--- ~
a.
rI.
t)l{V+v-
t~.
Wind Excitation
Fi g. 1
t~
I
IAt~.(\,~
TV 'PV"
290
types of gauges
291
~"
accelerometer
n $&i .mometar
~-----
==:c--
.~..
magnetic
Cl
II
:/~~~':
:frequencies
moment
di.splaceL'1ent
~
Ji
stress str..-.in
stra:l.n gagee
/'
:frequencies
di.!lplac6mt'nt
#.
e1ectronic theo
dolite
"trackercoil or
capnc:l.ty induc-
tion
"
~F
displacemont.Creq~eDcies
s tntic.d)"1l.:u,1iC
opticnl measure
displacell1\)n.t
laser
hologram
de:f'ormation
.~-~liGht
.,x-~
stat1c + dynamic
!'requencies
f'requencies
weighbridge
axle load
s'trnin sage
tr,,:fNc ..,tn'tisticlf
quartz
IIUltpletic 0011s
speod. intern\l:J
Fig. 2
~ j
~~
Fig. 3
I
I
IF
292
Frequency analysis
For the treatment of vibrational characteristics
different types of spectrum analysis (Narrow Band) are now
popular. The common mathematical basis is the Fourier Transform.
The discrete equivalent of Fourier's integral, called Discrete
Fourier Transform D.F.T. can be written in two equations:
n-l
_ 1
f(n) exp (_j 2~nk)
forward transform
I:
F(k) -"N
n=O
N-1
inverse transform
f(x)
I: F(k) exp (j 2~nk)
k=O
The Fast Fourier Transform F.F.T. is more efficient in
digital computing ; th~s algorithm reduces the number of complex
multiplication from N to N lo92N.
Analysis bandwidth T
Spectrum or
S~gnal.
period~oally
~t ~i
....J.-7'
Fig. 4
A typical bandwidth choice would be approximately 1/5 of the
bandwidth BAff = iE ; TE is the length of the pulse.
Nevertheless it is possible to increase Trep to obtain a better
likelihood of encountering problems with dynamic range and crest
factor. The value of bandwidth Beff can in most cases only be
fixed by analysis of the data (If we gradually vary the bandwidth,
a stage is reached where the results remain unchanged. Only if
we want to express the results as an "energy spectral density"
the bandwidth must be less than the value of the transient itself
B :< Beff .
Window length and incremental step
It is always desirable to make the step somewhat shorter
than T so that there is a certain amount of redundant information in the succession of analysis
293
Fi g. 5
Window type
The choice of the window type is also important.
The Hanning function as well as the rectangular function are
easily generated and give good results.
o t)'"'-~~_f~_":f.I.-1X :1.
-to ~:i'." I ! ! ' , ~~ . l.,
20 :
i:::: I I i II
:;0
i Li.
bo
~H(:
10 ;:;;:':.
III
! III
;~ ~l ii~i; i: !
dB
Frequency
.:
Att.nuat~on
j.'-
20 dB
Decade
60 dB / Decado
~
1-1--\
-{I\i----Corrtperi.on
or
F1nt woiehtin~
(rectansular)
~ctlon
Hanning
(cooine)
Fig. 6
294
At
.r
HEASUllli
Fig. 7
attenuation
Consequently the above mentioned report becomes much
more favorable.
,
s
.AMPLITUDE
Accelerltion spectnn CIIftCnte bridg. :
30 II $'ngl. ,pin SUbJKted to rood
tl'lff!e
Fig. 8
'1
II
FREQUENCIES
295
hmpl1.tude
~\,~
-
A'
A.
<111
~.'"j ~ k;-"~'<O
III
A.
t:
Fig. 9
Autocorrelation function, Convolution, Transfer function
This function is obtained by multiplying a function
fx(t) by the function fx(t + T) as follows:
Rxx (T)
1 JT/2
lim f
T-+-oo
fx(t) fx(t+T)dt
- T/2
f(t) and
h(t). The
296
+OO
9 (t)
-00
f (t)
1-1----.-..-
~_--,g (t)
input
~I
h(t)
out put
and H (f)
(f)
= F (f) .
of h (t)
H
and
then
(f)
fB
297
298
=a
factor of elastic
sensibility which
depends on the length
L of the element
0,4
a
= 'l-+"--;0",""2-;-L
1.5
----SO
100
BRIDGE SPAN m
Fi g. 10
factor similar to a
mechanical admittance ;
dynamic load
(vehicles) ;
2. Current research
Traffic conditions as well as higher stresses in
material as lighter bridge decks lead to greater stress changes.
Therefore the vibrational behavior and the risk of fatigue
damage must be taken into account.
In order to study the relation between traffic load
and stress in the different bridges components, the European Coal
and Steel Community sponsored a large investigation program.
Some interesting research on existing bridges has already been
undertaken by various laboratories.
a) Road-traffic loads on a steel bridge
A single carriage way bridge with only two lanes, built
with a simple structural system has been chosen and studied by
Mr. J. DE BACK and MR. M.H. KOLSTEIN, Delft University. The
Netherlands.
The bridge
A movable bascule bridge crossing a river has been
chosen : "the Haagsche Schouw Bri dge" on Hi ghway 44 Den Haag -
299
Amsterdam.
A
B
C
D
E
k- _1._ _ _ _ _
7.40 m _ _ _ _
9.50 ..
Fi g. 11
Measuring experiment
Strain gages were applied to a cross girder (bottom
flange) and to the stringer located under the traffic lanes.
----,"l)OO~
DRI VE DIRECTION
FAST LANE
"
"
Fi g. 12
Weightbridges and magnetic coils were installed a distance from
the bridge in order to measure axle load, vehicle speed spacing,
and to give all statistical traffic data.
All required data for the traffic load as well as the
resulting loading effects are recorded and processed with a
minicomputer. All the electronic equipment is situated in a
measuring lorry in order to produce immediately 6 histograms of
300
and the stress at the measuring point was recorded. The calculated
stresses are about 30% higher than the measured values.
- Dynamic response of the bridge to moving vehicles
- Moving vehicles transmit to the bridge-deck a local dynamic
action which fluctuates around the static load. The direct
impact action can be defined by an impact factor: maximum
force between wheel and road divided by the corresponding
static load; the factor for trucks and trailers can vary from
1. 2 to 2.1
- Moving vehicles generate vlbrations in the bridge structure,
and the resulting stress fluctuations at a measured point also
depend on the general oscillation movement of the bridge even
if the vehicle is far away.
This can be seen in the signal obtained by the passage of a truck
and trailer over the bridge. This analogous recording enables us
to clearly see the 4 maxima in spite of the superposition of
301
Fi g. 14
Bridges subjected to all road traffic vehicles
The real problem in dynamic requirements of
structures is the passage of vehicles at different speeds, lateral
position, load, spacing.
Even though the axle-load of the traffic and stresses
are measured simultaneously at any given moment, because of the
random superposition principle, the mechanism of loading actionloading effects must be statistically studied. The vehicles
composing the traffic sample have been grouped according to
position and number of axles in 21 groups.
The main factor in the traffic statistics are vehicle
group load, spacing, velocity and lateral position.
In regard to stress, the number of level crossings
is of major importance. These measurements, carried out on the
cross girder, clearly show that the slow lane bears the
strongest stress.
Fig. 15
10 01
"'ll/Nl:l
302
Counting methods
For a given composition of traffic, the experimental
stresses in the structure are analysed for a chosen time period.
Two counting methods can be used in order to work out the stresshistogram. The classical levelcrossing method and the rainflowcounting method.
In the rainflow counting method, the collection of the
stress variation with time (vertically represented) is similar to
the different roof surfaces on which the rain flows. Each half
fluctuation of the signal corresponds to a roof surface which
leads the rainflow from a maximum to a minimum and inversely.
o, - -_ _ _ _ _ STRESS
Fig. 16
o
1'>
l'
TI C
303
.,
"
Fig. 17
Each vehicle type has a statistical frequency for
crossing the bridge in the slow or fast lane in a given period of
time, and it contributes to an increase in fatigue. In order to
estimate with precision the damage potential, a new concept of
prime importance has been used.
Characteristic axle-load
based on the equivalence of the fatigue damage potential
between
[characteristic axle-load]N times = N axles of the axle load
spectrum
304
For each vehicle, this value can be calculated. For a bridge having a defined structure and a given distribution of traffic the
potential damage by each vehicle type can be estimated for a
given duration.
characteristic axle-load
Obviously the sum of the potential damage induced by the
different vehicle types increases with time. The time needed for
this cumulative damage to reach the value 1 is nothing more than
the theoretical durability of the bridge under the considered
traffic data.
The application of the notion of characteristic axleload is extremely viable. It shows that for this bridge, 91% of
the accumulated damage in the slow lane is caused by only 5% of
all traffic.
In particular, it points out the most damaging types of
vehicles. Several other bridges have been studied by the same
laboratory. Even though, the structure of the bridge and the
traffic data are very different, the results are similar. The
above power relationship made it possible to calculate the damage
potential of a load range class or a stress range class with the
Palgren Miner damage rule.
For a given vehicle type the characteristic axle-load
can considerably differ for two different bridges, but also-for
each bridge type, the fatigue life expectancy can be obtained.
This datum is fundamental for future bridge codes.
b) Railway bridges
Some interesting research on bridges required for traincrossing has been realised by M. M.R. Pfeiffer at the Frauenhofer
Institute of Darmstadt (W. Germany).
Axle-loads are measured on the rail itself due to an
ingenious strain gage system which only delivers an electric
tension when the wheel is situated between points 1 and 2. At the
same time this system gives the count and 2 of these systems give
the speed.
Fi g. 18
Rail + wheel
strain gages
305
The axle-loads are much more homogeneous and the axleload spectrum, obtained by the passage of a passenger train is
shown by Fig. 19.
Fig. 19
log
dB
20
Hz
10
10 1
10 2
Main girder
longitudinal stringer
70
STRESS
N/ ...2
~.
60
LONGITUDINAL BEAll
- - - - TRANSVERSAL GIRDER
_ _ MAIN GIRDER
_.
-- .--._.
~J -~~))
-10
Fi g. 20
~-----=r~'::' --,-~
10- 1
10
10 1
10 2
10 3
Log
~f
104
306
Conclusion
For railway as for road traffic subjected bridges,
the dynamic response can be written in terms of stress cycles.
The bridge design method, in order to take into account the
potential fatigue damage of structural :.elements is inseparable
from the stress spectrum which must be converted into a design
spectrum. Finally for a given traffic, it is the design spectrum
which improves the fatigue life expectancy calculations.
This type of research on bridges subjected to traffic
clearly shows the insufficiency of the classical notion of
dynamic factor or homogeneous stress for ensuring saf~ty
This research introduces a more efficient approach: it
means attempting to level the cumulated damage for all points of
the structure. With this concept all elements of a bridge
structure have the same life time expectancy.
III. DYNAMIC TEST ON ACTUAL BRIDGES
1. Dynamic test with vertical impulse: Auscultation of full
concrete slab bridge
It concerns a two-lane road-bridge, with short span.
This full slab bridge (cast-in-place) lies on 4 elastomer
bearings with a dead load of 120 metrical tons. The principal
geometrical characteristics are given in fig. 21 and 22.
The dynamic test was made with a mechanical system.
It consists in dropping a 50 Kg-steel plate from a hei9ht of
0.7 m on springs.
The dynamic characteristics have been obtained
through the classical method: accelerometer, recorder, analysis.
Cross section
Fig. 21
9.00
307
/,
-;
/:
?
/
,..
/,
"
'/
/
A:
2 x 3 x 2 12 pilln
, 10 <II
L. 570 ~
: .. terttght 3 ell
ISptNil t top lIyer 9 CII
C :
Fi g. 22
The two main frequencies very clearly appeared, and
whatever the conditions of measures might be, they hQve the
following value.
2,95 Hz (calculated 3)
Fundamental mode frequency
"
10,3)
11,0 Hz (
Second
"
Third
c
23
Hz (
23)
The damping obtained by logarithmic decrement has also
given a constant value, similar for the 2 modes 0 = 0,01.
As in most concrete slab bridges, geometrical imperfections, errors of concrete quantity (mass), errors of steelbar
positioning (inertia) ... give experimental frequencies slightly
different to the theoretical ones. The amplitude of the vibrations
depends on the position of the measuring position M. Fig. 23
shows the response for 2 characteristic positions of M : in the
middle and at 1/3 of the bridge span.
It would be extremely interesting to redo in identical
conditions this test in order to follow the mechanical behavior
through the years.
308
Acceleration
M
,
tt,
i},-t
'"'--:la
20
Fi g. 23
10
Aropl itude
Fi g. 24
309
* The
* the
. resonance method
A constant sinusoidal force with a progressive
variable frequency is applied. For each frequency N. the
1
corresponding amplitude Al is measured. For the resonance,
frequency N2 ' the maximum of amplitude reaches Ar
If 2 different frequencies NI and N2 (so that Nl < N2 < N2 )
produce the same amplitude A = A = Ar
1
2
-= the damping can be
/2
obtained:
Nr
N2-N 1
Q surtension factor
1
'V
=Q
= 2Q
310
AMPLITUDE RATIO A
AD
AMPLITUDE
1
...
.,...:
1[2
0.5
...
...
"
".
FREQUENCY
...
~REQUENCY
N1Nr N2
No
resonance
Impulse spectrum
Cross section:
I~
Fig. 25
311
DAMPING
resonance
steel girder
1.19
1.23
109.
decrement
2.35
.. wi th nonnal
1.19
1.15
2.60
2.84
1.18
1.15
2.50
2.50
mea::.ured
concrete deck
concrete deck
computed
2.24
Fig. 26
312
;I
Damping
8 %
stee 1 qi rder
'-- . ....c
..mm
I
6
Amplitude
..
............. p
.... ....
3
2
.... - ....
2.5 mm
. . . ...
'
10
1'3
mm
Fig. 27
IV. DYNAMIC TEST IN LABORATORY
The dynamic behaviour of a large number of bridge
elements can be usefully studied inthe laboratory with dynamic tests
to assess durability.
1. Accumulative damage
If n1 n2 .. , are the number of stress cycles applied
to the material at amplitudes of stress 01 02 .. The o-N Miner
relationship give us the linear accumulative damage criterion D
in the discrete form
313
'-:I"-::!---:---;----,;+
"1 nj
number of cycles"
stress spectrum
stress histogram
c~mulative
Fi g. 28
Fig. 29
Fig. 30
damage
Fig. 31
constant amplitude
cycling
in adviln,=e' fixed
bloc. eyel'!!..
l'"3ndcm proces$
c./cl1n']
314
315
Fi g. 33
Fig. 32
Phase
II
Phase
III
Cumulative damage
ratio 0i
i l FAILURf
O.S
O.S
Fig. 34
1 N
cycle ratl0
2. Damping
Themodal analyses of a bridge furnish the frequencies
and mode from the mass and stiffness matrixes without necessitating other hypotheses. On the contrary the damping value has to
be introduced arbitrarily. This value is based on an experimental
damping behaviour in similar bridges, generally small. The errors
316
-f ~
+ -fFORCE
FORCING S Y S T E M .
AMPLITUDE
Fi g. 35
Fig. 36
317
n=3
n 2
n I
n 0
'--_ _ _ _ _ _ _ _ tillle
Amp itu
Fig. 37
Fi g. 38
Considering now the influence of the prestressing stress
on the value of 0, tests on the prestressed beams show clearly
that from a certain stress level 0p stabilisation of 0 is
obtained,Fig. 39.
6 S
-1.5
Fi g. 39
3.5
2.5
~~~-~~~~-----------+----~--~
100
Variation of
..
0pdaN/cm 2
O"p
318
FAILURE
"-
n 0
I!
Fig. 40
II!
n I
AMPLITUDE MAX.
319
10%.
320
40
)J",i/e
..\
Fig. 41
Characteristics such as rail-wheel creep coefficient, stiffness
of the vertical and horizontal wagon suspension, hysteresis
damping must be taken into account.
321
~
._./
,.
"
Lat.. ,.al
.s
Fig. 42
J"ne ..,."
322
Fig. 43
running gradient
track grad i ent
for V: 100 km/hr
323
Speed
.........
\""".j.":":
....;'
...........
__ 4
65 kmthr
AmI,
'Ih,. ,.
0 5 0 '00 15021>0250
Dynamic amplification
- running gradient
.... track gradient
no
0.030 9
250
0.162 9
Fig. 44
Fig. 45
324
400 m train
6 m
5
--r---.,~---
4
3.l---~'"
,,
o
\
\
lBridge Type 2
~------~----- -...
,- - - -/
./
.'
Bridge Type 3
"
300m
1l00m
Fi g. 46
""
""' "-
",
""
300m
'.
""
325
~/.~
,
~ia
4oo~
-f~vtZ.ncitt!)
('t>""we..'.! ,"~'>IJrc.! Damp/~
0.64'''
~683
.t
0.827
0.4'(.1
0.11"
0.1'7
1.3/0
/.361
1./U4
1.3.B
0133
0.100
o.4'J?
ad'6;t
o.(.r! 0120
0.199 0.8'-10
1.3'3
1.34'3
U9J
1.300
a.f.].!
a100
0161
O.sr-o
0.of1e:.
ao~!>
.z
"3
o.,ro
o.no
aJ{'1
1./5'
(J,
~-
1.21/
'----'----
a.ol"
---
Dalnp;"..!!
c.ot,s"
0.0.288
aal4~J
I. aU 1.233
- - ._------
Fig. 47
*"F1,.
road
Ci;~:ar
type)
*Te-tr()n-
circular type
48
326
~ ~'
-- !r
~
/
I a..
1.".-
---
l"
* ORE
* SETRA
Fig. 49
[JJ1D" ;...,.,..-1
I~
7:
'l~
A_
of a bearing
327
Horizontal forces
All the dimensional variation of the bridge (temperature, creep
... ) with norizontal forces applied produce distortion in the oearing and a shearing stress a H is developed.
~
.
II
Ii
"___ 2.~
f~
Fig. 50
Rotation
In addition an elastomer plate can be a subjected to a
rotation a and the corresponding stress is T . When the
linear or angular displacements are too large, a the bearing has
to slide so that overly larger efforts can be avoided on the
bridge deck as on the foundation ; the friction coefficient on
a PTFE plate is experimentally given (approximately a few percent).
Dimensioning
Whatever the types of actions applied to the elastomer
bearings the shear-stress is limited T = TN + TH + Ta < 5G
as well as normal ~)stress aM < 15 MPa
For a railway bridge the values T = 1,5 TN + 1,5 Tare
a
recommended because of the fatigue phenomenon resulting from the
great variation of the vertical ,loads.
By their elasto-plastic adaptation, the bearing allows
a good distribution of loads even with small differential
movements of the foundation; after an initial static deformation,
a supplementary adaptation must be taken into consideration.
Dynamic Behaviour
The visco elastic properties interfere at the time of
the transmission of dynamic charges. However,in case of road
traffic the dynamic components remain limited, even if there is
a small vertical oscillation around the average position, the
amplitudes remain low compared to the dynamic bridge deflection.
The influence on proper frequencies of the bridge consequently
remains insignificant. In any case the damping effect on the
bridge structure will be limited and can not easily be separated
328
Fig. 51
PERIOIl(
Ir--.--.----,.--__
on-"
0.24
:c:
AA'
~~----.----
.,
II
lr ~-~.
-.
I ~
O:'~UJ ___,j---
o.7'l
..
"'OJC
- ,.--,
1~QQ.~
172 m
320 bars
21 bars
Density = 1.B5
Mod. E = 21.0 x 106 KN/m2
(+)
1 110 tons
"28 (-) =
"28
4J 37.5
350 bars
28 bars
Mod. E = 39.3 x 106 KN/m2
"28 (-)
"28 (+) =
_6~~---1
144 m
prefab. elem
U gh:t concrtte
Traditional concrete
prefab. elemt
Fig. 52
__ ~ ___ _
-L.
00
-L
n0
330
331
I
~!L
i~ ~,,=o--I--::=*I:=x
:::::~-=~-=t-A---:====---A-::::I"
.42
-----,I
'
.60
.67
Fi g. 53
4
.85
1.03
Z 1.26
I--~~~---"""-<-
1.38
- -- -
- - ----
--~~=---
1.60
332
PERHAN t
LOAD
TYPE OF BEARINGS
~
c:
tons
SOLUTION II
SOLUTION I
::.:
0
+
+----
ident,cal
;,. .s-,
sS-'
0
1
2
3
4
5
identical I
Neop. blocks
00 x 900
9[20+3]
--L,",I;>
--
Ix
4090
Fi g. 54
identical I
2450
identical I
2500
4 Neop
block~
500 x 500
520
-t--+-
'2910
300
identical I
r~
11,
l1y
151,:'11
WOq
II
177B
bJ ':$3
Ii 116
.7-:5
2305
6472
40
II
7324
21,7
189?
175
Si34
2b67
4~O4
2B~
181
l6,
_I
IS"
11
I
~b91
l411
4419
190
:bl
II
6458
2276
4810
337
356
IS GIl
4432
:39
JS"
I} ,"(,1
27
ZO 173
J~:';
11
MuD
~1~!J
3701
7609
II
3825
2475
1691
2~6
110
.....
I
II
2430
19 79;
2418
1 338
Fx : kN
Fig. 55
"x
kll.m
I":.
333
IA,-,",
~fi/61"11C.
ACTION/F"I!~
1
t N-
N+ ton
T-
~
II
N.
1807
1433
1471
16~
203
409
453
:;85
163
142
H':
191
1632
1376
114
liZ
t NT-
Fi g. 56
334
values.
The dynamic behaviour of the girder-foundations unit is
inseparable from the mechanical characteristics of the bearings,
especially in case of seismic action.
The dynamic analysis of the mechanical system with
solution I and II makes an advised choice of the bearings possible.
Safety is increased by limitation of the forces in the whole
structure.
2. Effect of profile imperfection of the bridge subjected to
moving loads
The dynamic effects due to traffic are generally
introduced in calculations with a dynamic amplification coefficient.
This coefficient is not always the consequence of a
dynamic approach of mechanical response or of experimental
measures.
These coefficients depend on different parameters, the
different national code values show significant differences
because they are not based on the same parameters.
Most of the codes give a global dynamic coefficient
which allowsthe obtention from a static response of the so-called
dynami c moment and defl ection of the bridge.
A more subtle approach has been conducted by Mf.Palamas
(L.C.P.C. Paris). Mr. Palamas shows in fact, that the
coefficient has to be defined for each section of the bridge. The
dynamic majoration coefficient of the vertical deflection (Cd)
and of the moments (Cm) for a given section then depends not only
on the mechanical characteristics of a bridge and of the moving
load, but is also greatly influenced by the profile imperfections
and the vehicle speed.
Calculation of the bridge-moving load system was
realised by a computer program "Palmares", the numerical application was done for a classical concrete bridge and one vehicle.
The essential characteristics are given in Fig. 57.
When the profile is without irregularities and when the vehicle
arrives at the entrance of the bridge without vertical oscillation,
the force exerced is constant and Cd at midspan, scarcely
exceeds the unit value (Fig. 58).
335
llechanical charaterist1cs :
Bndge ..pan 34. or 11
2 ..
Yehtcule mass
+---~dL---~--~Q-----r
----~~z=-==-s=~~~==-====~~
"17" 16.7 HZ
1. 9.92 x 10 10 "/ra2
RlgtcJUy
t/~
31.7 t
suspension stiffness
suspension dal!ping
vertical period
veloclt,yV
To 0,37 sec
o< V<
100 KlD/hr
Fig. 57
1t,1).
S'"
JIIc"UC~.1II'
"tn/,.,
11.55
L. ]q",
Fig. 58
Fig. 59
Fig. 60
,.1
336
~~
rUD
I
Fig. 61
$~",
tit,."""...
L"
-I
f-.V= 1001,,,,/),,._
Fig. 62
Fig. 64
./
Fig. 63
Fig. 65
337
338
CON C L U S ION
The problem of the dynamic behaviour of a bridge is
complex. Different methods of measurements and analysis process
have been presented. The control of the bridge response under
traffic demands important infrastructures. Whereas the dynamic
tests, such as vertical impulse or sinusoidal excitation which
can be repeated in the same conditions, are the indicators for
survey and maintenance.
This paper shows that on the one hand different bridge
types must be tested in situ in order to establish the diagnostic,
and on the other hand, the dynamic behaviour of new elements and
material must be evaluated in advance by a more scientific
approach through laboratory tests and mathematical simulation.
These 2 aspects are complimentary. In the stage of bridge design
the life time prediction must be obtained with an acceptable level
of threshold. That means that the whole structure must resist
simultaneously two types of phenomena :
- mechanical loading, which introduces stress cycles with the consequence of fatigue and cumulative damage
- physical and chemical environmental effects and their interaction with the first one.
The mechanical characteristics can be considered as
indicators of the vibrational behaviour so that the durability of
a bridge or bridge element is nothing more than the maintenance
of performance in time. The value of frequencies and damping of
the different modes are with the stress the main indicator. Of
course their evolution depends on the environmental characteristics such as temperature effect, chemical action, corrosion,
time.
In laboratories,these characteristics can be more or less
simulated but the correlation with the in situ behaviour is not
easy when the time factor is of great importance.
In fact, concerning the environmental agents which
can attain the mechanical performance, the life time cannot be
correctly estimated. Exterior conditions, bridge type, load
cycles ... vary to such a great degree that present experience
is almost always founded on numerous spot observations made in
situ.
339
REFERENCES
CHAPTER I
- Bendat J.S., Piersol A.G. - Random data: analysis and measurements procedures. Wiley Interscience, New York, 1971.
- Guillard Y. - Vibration dues au trafic routier, degradation de
l'environnement. These 3e cycle. Fac. Science, Universite Le
Mans, 1978.
- Zeller E. - Comportement vibratoire d'une structure metallique
test impulsion et excitation sinusoldale. Etude Juin 1982.
CHAPTER II
- De Back J., Ko I stei n M. H., Van Leeuwen J. L. - r1easurements and
interpretation of traffic loading in bridges. Lewin Inst.
Univ. Technology Civil Eng., DELFT, 1981.
- R. Pfeiffer, J.M. Zaschel - Ermittlung der Verkehrslasten und
der Verkersbeanspruchungen von Stahlbrucken als Grundlage fur
kunftige Bemessungsregeln, LBF Fraunhofer Institut, Darmstadt
FhG, Berichte n 4, 1979.
CHAPTER III
- J. Caracilli - Constatation sur ouvrages metalliques en
construction et en service. Etude Paris 1980.
- L. Ferron, E. Zeller - Comportement vibratoire d'un plancher
industriel, du calcul aux tests dynamiques. Construction
Metallique nO 4 1982.
340
CHAPTER IV
- H.P. Lieurade - Essais de fatigue sous sollicitation d'amplitude variable, RFP 255, IRSID 1978.
- C.P. Heins - Bridge response and damage, J. Washington Ac.
Science 1979, 69 n 1.
- M. Gregoire - Interpretation des resultats d'analyse statistique de sollicitations relevees en service, CETIM. 1961
- B.T. Yen - Minimizing fatigue and fracture in steel bridges,
AS ME 1979 Conference, Washington.
- J. Caracilli, A. Aloustaz - Influence de la precontrainte
sur 1 'amortissement des vibrations transversales des poutres
en Beton Precontraint, Bulletin Liaison Labo P.C., Mai-Juin
1978.
- J.L. Van Leeuwen, AJM Siemens - Miner's rule with respect to
plain concrete, Heron vol. 24, 1979.
- H. Van Koten - Structural damping, TNO CUR report n 75, 1977.
- Y. Yasoshima, Y. Matsumoto, T. Nishioka - Studies on the
running stability of railway vehicles on suspension bridges.
Journal of Fac. of Eng. UnivcrsiLY of Tokyo, n 1, 1981.
CHAPTER V
- Appareils d'appui en elastomere frette, Bulletin technique
n 4, 1974.
- Palamas J. - Imperfections du profil d'un pont et vibrations
sous trafic, recherche effectuee au LCPC Paris, These 1982.
- Y. Guillard, E. Zeller - Initiation Dynamique des Structures,
constructions parasismiques; mesures sur ossatures industrielles
et ponts sous trafic routier. Cours Ecole d'Ingenieurs
Geophysiciens - Institut de Physique du Globe, Strasbourg.
v.
343
C. P. Heins, Professor
Institute for Physical Science and Technology and
Civil Engineering
University of Maryland, College Park, MD. 20742
INTRODUCTION
With the
~dvent
Fernando Earthquake and more recently the 1978 Santa Barbara Earthquake
(1), bridge structures in the USA have undergone considerable destructive
forces.
The earthquakes have caused the bridge professionals to reassess the
design techniques that have been applied, up until that time, for seismic
design.
The present 1977 AASHTO bridge code (2), as related to seismic
design, suggests an equivalent static force method for simple structures
and when the structure is complex, for example curved bridges, a
computer based response spectrum or dynamic analysis should be considered.
In this present 1977 AASHTO code, most engineers would utilize the
seismic coefficient method (SCM), however, the utilization of the SCM,
may give erroneous results when designing under seismic conditions (3).
3M
Because of such errors and from experience gained from recent earthquake
the 1975 AASHTO code was reassessed in 1977 by the ATC (4) research group.
The work of this council, developed a new specification (5), which will be
given in part herein, relative to straight girder bridges.
Although the ATC-6 design guide (5) has advanced the state of seismic
design of highway bridges, it was recognized (6) that additional research
was necessary to determine the seismic response of bridges which are
influenced by geometry and curvature.
resulting in a proposed design guide for curved steel box girder bridges.
These results will be given in detail in the following section.
STRAIGHT HIGHWAY BRIDGES
The analysis techniques, as recommended in the proposed design guide (6),
are:
1.
2.
in which both methods can employ simple analytical techniques not requiring
dynamic computer programs (10).
The (MMSM) requires modeling of the structure as a lumped mass or
distributed mass continuum and application of a given response spectrum
or time history.
345
shapes and eigen vectors and actions requires use of a computer program (1).
Such application is generally warrented when the structure has an unusual
configuration.
The details to be presented herein will only be relative to the (SMSM)
technique and direct design criteria.
General - The development of the (SMSM) technique requires utilization
of the generalized displacement function concept.
Eq.
(2)
where;
Amplitude terms
~n
displacement functions
l\(x, t)
where:
~(x)
Z(t)
Z(t)
Amplitude of motion
~(x)
Shape function
(3)
The generalized displacement function as given by Eq. (3), will now be used
to develope design oriented equations.
Formulation of SDOF Equations - Consider a SDOF system, as shown in Fig. 1,
subjected to a ground motion.
346
mli
where: {f1
(4)
= cl:J.
fD
(5)
fS = kI:J.
The virtual work method will now be used to obtain the general equations
of motion; and is given by;
oW =
where: JrF
ol:J.
=0
(6)
Applied Forces
lLol:J.
IF.
= Virtual
Displacement
Z(t)
+ /).
(t)
(7)
fD
c[1jJ(x)Z(t)]
fS
k[1jJ(x)Z(t)]
1jJ(x)oZ(t)
(8)
347
= f:m(x)[~(X)Z(t) +
JL
~g(t)]~(x)dx6Z(t)
+ oc(x)[w(x) Z(t)]W(x)dx6Z(t)
Since
6~
f:k(X)[W(X)Z(t)]~(X)dX6Z(t)
~(x)
6z(t)
s=
[w(x)2]dx
f:m(X) [W(x)]dx
(9)
J:C(X)
c*
k*
= J:k(X)
[~(x)
+ c*
Z(t)
] dx
therefore;
y Z(t)
+ k*
Z(t)
= -!J. g (t)
or
Z(t) + (c*)
y
Z(t) + (k*)
y
Z(t) =
..
-~
(t)-
(10)
Let:
~
c*
2m (x)w
=-;:---::;~
~ (t)~
g
Y
(11)
348
l1(t) = t:"(t)
where:
and
III
III
= JtX (T)exp[-~(t-T)sinlll(t-T)dx
og
t:"(t)
l/J(x)
Z(t)
l/J(x)~
!let)
III
(12)
Local Forces:
Neglecting damping:
= m(x)
fl(x,t)
i;(x,t)
b.
III
A cos
-Ii
Substituting
..
III
A sin lIJt
b.
= A sinwt
gives
11=-1Il11
2
therefore
or
fS(x,t)
fS(x.t)
= m(x)
= -fS(x,t)
[m(x)
(13)
349
Maximum effects:
Displacements
Using the response spectrum, the spectral velocity, which by definition
is a maximum of the response integral or
~(x,t)max
1
s
w v
n:uS
sd
mS
therefore:
~
max
and
S
cs
.. t::max
c s g/w
(14)
Forces
a-
fs(x,t)
fs(x,t)
m(x)w(x)yw~(t)
=
m(x)
w(x)~
m(x)
. w(x)~
y
Sa
wOJ
However
Sa = g
{
m(x) = w(x)
and
350
=~
y s w(x) ~(x)
fS(x)
(15)
where:
Cs
= Elastic
Natural Frequency:
The natural frequency can be determined by equating external and
internal work.
f Fll
WE
where:
= Po
2
JL t.(x)dx
0
po. a
2"
f:t.(X)dX
2
00
~ --
max:
max:
~ 2dx
w(x)(t.)
s
max:
= !!LL
fo
L
w(x)(t. ) dx
2g
max:
P
o
2"
2
2g
_~
a -
max:
and
00
21T
(16)
Po
351
The coefficient
where:
A
= Acceleration
Coefficient (Fig. 3)
General Coefficients A, S:
The coefficient A and S given in Eq. (17) are prescribed as follows:
Acceleration Coefficient A. - The acceleration coefficient A, is given
in Fig. 3, as applied to the United States.
Site Coefficient S. - The site coefficient approximates the effects of the
site conditions on the elastic response coefficient and is;
(S) Coefficient
Soil
I
1.0
Profile
Type
II
III
1.2
1.5
foll~'s:
352
Stiff soil conditions where the soil depth is less than 200 ft
(61 m) and the soil types overlying rock are stable deposits of
sands, gravels, or stiff clays.
Po
~(x)
STEP 2
Calculate parameters
a
y
where
J ~(x)dx
J m(x) w(x)dx
J m(x)~(x) 2dx
a,
or
or
structure.
8, and y
f w(x)~(x)dx
J w(x)~(x) 2dx
w(x)
353
STEP 3
Calculate natural period;
unit load
STEP 4
Calculate equivalent static earthquake loading:
f s (x).
=~
w(x)$(x)
y s
(force/unit length)
STEP 5
Apply loading fs(x) to structure and determine resulting member
faces and displacements.
354
modeling of the structure such that the support conditions can represent the
physical restraints.
~.
and
(Iw,Wn,Wa) can be
computed. they were not considered in this study. as it has been shown
(12). that with proper bracing of the box girder, warping and distortional
effects are negligible.
355
to determine the effective earthquake effects, using the response spectrum curves.
The general procedure using this method will now be described.
The support restraints, to be imposed on the bridge model, can be identified
as releases in the computer model.
Equivalent Dynamic Analysis.The natural frequency response of a single degree of freedom system can be
predicted by the following;
w
=1 ~
(18)
where
k
spring constant
tained therein are then used to determine the induced dynamic force
M . ~.S.
356
where
~.S.
curves.
If the system is subjected to angular accelerations
e,
dynamic torque is
M
where:
(20)
laRS
pr 2 dA
eRS
1J
steel
J<x2+y2)dA +p
concrete
J(X 2+/)dA
Natural Frequency.As indicated by Equation (18 ), the natural frequency of a SDOF system is
given as a function of the spring stiffness (k) and the spring mass (m).
If the
(21 )
where:
~
= induced
Therefore by determining the response of a given box girder bridge, when subjected
to a unit load, an equivalent spring constant can be obtained.
This constant,
k , and the total mass of the bridge, will then permit evaluation of the natural
frequency as given by Eq. (18).
In the instance of curved structures, the dynamic action can occur in three
The
357
resulting maximum displacement. induced by these unit loads. will then give the
correspond1ng equivalent stiffnesses
kx k_". k
y
z
The corresponding natural frequencies can then be determined by applying Eq. (18).
BRIDGE STUDIES
Typical Sections.In order to develop a simplified design technique. the response of
various curved box girder bridges must be examined.
which have been used in previous studies (13). were utilized in this
parametric study.
system were considered because this is most typical of the curved box
girder structures (12).
The basic span length configurations that were examined, are shown
in Fig. 9, where L=50' to 150' and N=l,2.
k = 0 ~ 2 x 10 3 k/in; K
x
y
co
0.5 x 10 3 k/in.
The equivalent seismic resonse of the single, two, three, and four span
structures were examined.
for all bridge spans and their corresponding induced accelerations were
then obtained (7).
given by
and
wt )
358
have been plotted as a function of (i) and the span length (7), for
the single span structure.
accelerations for the continuous spans and the Single span however have been
plotted as a function
F
F, where
(22)
DESIGN CRITERIA
Trends.The seismic design of continuous curved box girders were related to
the response of single span curved girders.
2.2 (~)2+0.0llL+0.45
(23)
-0.016(L)+4.7
(24)
3.8
(25)
Ax
....L
Y
A
z
359
(26)
at = 1. 5
(longitudinal)
F
= -0.02
where:
F
J...
(L)
+ 3.75 +
= 0.00125
(27)
1.0
(vertical)
(28)
(transverse)
F
(29)
-0.005L + 1.5
F t (torsion)
Ft
(30)
= -0.075
where:
NS
(NS) + 1.15
= number
of spans (2, 3, or 4)
Design Approach.The equivalent seismic design of curved box girder bridges will
incorporate the primarily developed equations and the effective peak
acceleration map (Kp )' given in AASHTO.
is of the form:
EQ n
= Fn
EQ n
= Fn
A
n
Translation
(31)
K
p
Rotation
(32)
360
where:
Fn
An
= single
span acceleration
= wI = total
weight of structure
gravity
EQn
for the specific direction n, the continuity factor F and single span
acceleration A are given by:
n
Fx
for R
<
600'
>
600'
Fz
Az
3.8
Ft
At =
361
EXAMPLE
In order to illustrate the application of these proposed equations,
consider a two span curved box structure of spans 100'-100' and
radius R = 600'
x
Ax= 2.2
(*
)2 + 0.011 L
+ 0.45
100
2.2 ( 600
)2 + 0.011 (100) + 0.45
1.61
A
J....
A
z
Az = 3.8
At = 1.5
the continuity factors can now be computed as:
F
x
F
(K
0.00125 R,
600')
362
F
J...
Fy
Fz =
Ft
= 442.15
Kips
= Fx
W
A
(g)
(2.5) (1.61 g)
4.025 W
EQ
=F.A
Y
Y
=
(3.1) (0.9 g)
2.79 W
w
g
= (1.0) (3.8
g)
w
g
363
3.8 W
EQt
Ft. At
the equivalent induced uniform seismic forces qx' qy' qz' qt' are therefore:
E~
4.025 x 442.15
0.7415 kips/in
qx
-TX'"""L
2 x 100 x 12
EQ
2.79 x 442.15
Y
q
0.5140 kips/in
2 x 100 x 12
y
2 x L
EQz
3.8 x 442.15
q
0.7 kips/in
2 x 100 x 12
z
lXL
EQt
1.5 x 2273.63
1.4210 kip-in/in
qt
2 x 100 x 12
2 x L
Evaluation of the reactions of curved box girder bridge, using a
straight girder gives:
i)
x direction
I' E E E E (
qx
Il
.
2
Rx 1
qx
Rx2
Rx 3
0
0
(2L)
EQx
1779.65 kips
364
ii)
Y direction
rf--~L----ft--~L~----f
= 0.375q y L
Ryl
= 0.375qy L
1. 25qyL
0.375q y L
0.375 (0.5140) (100 x 12)
Ry
231.3 kips
1.25qy L
Ry3
iii)
Ry1
231.3 kips
z direction
)~+~+~~~~+~~+~~+~~~~~~+~~~++++++++++
~l
~2
' _ _L=--_--+'
Rzl
Rz2
0.3759, L
z
Rz1
~3
=
=
Rz3
= 1.25q'L
,Z
.075q zL
0.375 (0.7) (100 x 12)
315 kips
Rz2
1. 25q z L
1.25 (0.7) (100 x 12)
= 0.375q L
z
365
1050 kips
Rz3
Rz1
= 315
kips
iv) Torsion
0.375QtL
Rt 1
= 1. 25Qt L
0.37SQtL
0.375 (1.4210) (100 x 12)
639.45 kip - in
Rt 2
1. 25QtL
Rt 1
= 639.45
kip
in
=d
2131.5
48
d = 48".
44.41 kips
comparing these results with the results obtained from the (SFS) computer
program gives:
366
~
Direction
No.
Eqs
2
program
E.qs
program
Eqs
program
x direction
Rx
1779.65
1298
223.9
y direction
Ry
231.3
250
771-
775.3
231.3
250
z direction
Rz
315
417.5
967.7
315
417.5
Torsion
639.45
869.8
639.45
869.8
* Ry x
d= T
1050
44.41
55,5*
367
P(t)
or
grd motion
m.6 g
Fig. 1
o(x)
J~~----,~ ~
rt t f t
f f f t f t r~
Po
Fig. 2
368
F'GURE 3.
,/(x)
~I I f I ! 1 t
-- - -I"'(X)
~x
(a)
Po
r r
H"'(X)
Po _ _ x
(b)
Fig.
roo
...,
lJ>
'1
"..
()Q
>
...>-
!:
0.25
'" ,"
3.6
\: ,"
FREQUENCY( CPS )
< ""","
33
,,',<'"
10
i 1000
1000
01
\OX, ""X"
...,
.....
0'
(1)
"...
()Q
''', """,,,"
'"
,,','"
"'"
'"
< ,.lOOCI
>!:
u
~
5
025
2.:>
9
FREQUENCY (CPS)
l3
1000
01
_100~~~~100
tOOO,
-..J
....
371
0.002r---t-~--~~r-~+-~~~~~~~~-4~~
0.001
L....__J.__ _~~~__J._~:.......J~_~_.l_~~~~
0.1
0.2
0.5
1.0
2.0
FREQUENCY (CPS)
MEAN
Figure 7
372
(b)
(0)
F
Fig. B
Fig. g
373
REFERENCES
1.
2.
3.
4.
5.
"Seismic Design Guidelines for Highway Bridges:, ATC-6, Palo Alto, Calif.
October, 1981.
6.
7.
8.
9.
13. Heins, C. P., Sahin, M. A., "Natural Frequency of Curved Box Girder
Bridges", ASCE Structural Journal, Vol. 105, No. ST12, Dec., 1979.
375
Polat Giilkan
Faculty of Engineering, Hacettepe University, Ankara
ABSTRACT
The emergence of procedures related to increasingly more
precise earthquake structural analysis has necessitated the
implementation of similar procedures regarding the definition of
the seismic input. It is evident that bridges, which form a
significant groLlp within lifeline structures, should also be given
progressively greater attention during design and construction if
earthquake induced loads are among the primary load conditions.
Following a review of damage to bridges during past earthquakes, ani outline of methods on dynamic structural analysis is
presented. A survey of code requirements in effect in several
countries is presented next to emphasize the manner in which
theoretical considerations are reflected in regulatory documents.
This is followed by a thorough discussion of procedures for
estimating earthquake ground motions. The report is concluded with
the details of the earthquake analysis of a long multispan railway
bridge planned to cross over a large darn reservoir.
376
1.
INTRODUCTION
377
detailed according to code requirements, and this process will
result in perfectly acceptable structures. For the more important
cases, however, the procedures described in the following should
be considered, and code provisions should only be used for checking
the results of the more detailed dynamic analyses or for "calibrating"
them.
The rest of this paper is devoted to several considerations.
First, observations based on the past seismic behavior of bridges
and damage patterns will be described. This will be followed by a
discussion of the code requirements in effect in several countries
regarding bridge design for seismic effects. Considerable attention
will be paid to the methodology for developing vibratory ground
motion criteria at a specific site in a format suitable for use in
analysis. In the last section, an example will be provided for a
railway bridge currently under study and sample results will be
presented.
2.
378
of backfills resulting from compaction is often observed.
3. Seismic damages due to vibration effects are less common
than damages due to other effects; however they may occur
hi taller bridge structures which lack sufficient strength.
stiffness, and strength. A commonly encountered form of
failure is the insufficient length of bearing supports,
hinges and other devices for preventing spans from falling
off.
4. To minimize damages, bridge structures should be designed
with proper recognition of the stability and bearing
capacities of foundation soils, force-deformation and energy
absorption characteristics of substructure, superstructure
and linking elements, the dynamic nature of structural
response, and the dynamic characteri~tics of all forces
acting on the complete soil-structure system.
3.
DYNAMIC ANALYSIS
+ I ci
{~}
IKI
. {x} = {R(t)}
( 1)
(2)
379
necessary to solve the corresponding eigenproblem.
3.1
2
{y} + diagI2w.I;:I.{y} + diaglwl{y} = II {R(t)}
(4)
1.
1.
1.
In Eq. (4) it is assumed that the damping matrix lei also satisfies
the modal orthogonality condition
{cp.}T lei {cp.} =0
1.
i '" j
i
(5 )
380
algorithms (4).
3.3
,x
where x_
and
(6)
yr
where y
r
rx
W r yr
= r rx + r ry + r rz
(7)
r ry
_{cjl } T
r rz
-{cjl } T
r
r
IMI
IMI
i MI
gx
x
x
( 8)
gy
gz
h(max)} = {cp}
rx
r
{cp}T
r
IMI
{I}
x
Sd (w ,I; )
x r r
(9)
{cjl }
{cjl } T
{cjl }
{cjl } T
r
r
IMI
IMI
~n
the
{I }
y
Sdy (w r ,I; r )
(9a)
{I }
z
S dz (w r ,I; r )
(9b)
(Ill)
381
aW
a.gM
(11)
(12)
In this equation,
the basic coefficient
factor related to regional seismicity (i.e. zone factor)
factor related to the importance of the structure, which
ideally contains its vulnerability
factor that considers the vibratory characteristics of
the structure, most typically its fundamental period
factor related to structural type (ductility)
soil factor
For superstructures with uniform distribution of mass over the
382
Cl.g -
3{'T
with
T)- O.ss
(13)
where
F
F C
(14)
New Zealand
Similar to the U.S. requirements the basic lateral force is
specified by
CI.
= C F
(15)
383
(.)
0.2
Zone
.!!
;;:::
Zone B
(.)
...
0.1
Zone
en
u
"iii
0
0.4
0.8
1.2
1.6
Figure 1.
384
2
3
Zone A
Zone B
Zone C
0.24
0.18
0.15
0.13
0.12
0.10
0.09
0.20
0.17
Category
Minimum Value
of F
1.0
0.85
0.7
Note
California
These criteria specify the following value for the lateral
force coefficient
ARS
= -Z-
(16)
In Eq. (16)
A
385
01
analytical studies
5.1
Determination of Seismicity
386
Steps
ESTIMATE SEISMICITY
ON BASIS OF
EXISTING
REGIONAL DATA
USER DEFINES
ACCEPTABLE
LEVEL OF RISK
USER PROVI DES
DAMAGE OR DESIGN
CRITERIA
No
Yes
ESTI .. ATE LOCAL SOIL AMPLIFICATION EFFECTS
No
USER PROVIDES
ANALYSIS
CRITERIA
Figure 2.
387
Date of occurrence
Epicentral intensity, and where possible, magnitude
Epicenter and hypocenter locations
Epicentral maps showing the position of all reported earthquakes within an 80 km radius of the site, and, ideally,
all earthquakes with intensity larger than V (or magnitude
larger than about 4.2)
Correlation of epicenter location with tectonic features
Frequency of occurrence for earthquakes in terms of
magnitudes or intensities.
This task is often undertaken within the format of an earthquake catalog compilation. Catalogs of felt or recorded earthquakes
do not in general cover a long enough time interval to permit
meaningful extrapolations into the future except on a board regional
scale (13).
Earthquake recurrence law is of the inverse logarithmic form
log N(M)
a - bM
( 17)
388
One rationale for classifying fault activity uses the available
geologic, seismological and historical data to specify the fault
activity descriptions. Faults are considered to be significant if
they satis fy one of the following conditions : (i) faults cross ing
the site vicinity are capable of fracturing during the lifetime of
the proposed structure, (ii) faults near the site are recognized as
being capable of generating large damaging earthquakes. Current
criteria indicate that a fault should be considered capable of
permanent surface displacement if movement occurred at least once
during the last 35 000 years or more than once during the last
500 000 years. Empirical relationships relating fault length and
earthquake magnitude have been derived from observations
5.3
~n(R+
30)
(18)
389
0.7
0.6
~
.
...'"...
..'"
0
;::
0.5
-'
u
u
u
0
...
II:
Q
II)
0.2
0.1
5
10
50
100
Figure 3.
390
5.4
391
w. (RD)
w (PSRV)
(19)
=
w2 (RD)
(20)
392
(ern! s )
(2la)
log Av
(cm/s )
(2lb)
log Vh
(cm/s)
(2lc)
log V
v
(cm/s)
(2ld)
log Dh
( ern)
(2le)
log D
v
-1.13+0.24 IMM
(ern)
(2lf)
-'
II:
......
...~
>=
...
!E.. ~
UII:
II:
II)
10
100
1000
...
t;;o
...~ .,
I!:'
II:
'"
!d
I
I
I
I
I
"~I
.....
-I
Figure 4.
i
U
...>g
...
>=
z
:II
....ffi...
I!:'
.,
II:
8~
It
MERCALLI
INTENSITY
I~
10
100
--.--r--T
o Vertical component
Horizontal component
ffi
.,...-'
:II
~
,.:
...z
...
...
...
t;;
:!
I-
.,
10
100
~~
I
I
It
I
I
r~t
Il I
\C
....
....
394
395
Table 3. Aplification Factors in the Newmark-Hall
Procedure (20)
Amp 1 ification Factor
~ercent
lmping
0
1
2
5
10
20
Acceleration
(PSAA)
Velocity
(PSRV)
Displacement
(RD)
6.4
5.2
4.3
2.6
1.5
1.2
4.0
3.2
2.8
1.9
1.3
1.1
2.5
2.0
2.0
1.8
1.2
1.0
Displacement
A(33 Hz)
B(9 Hz)
C(2.5 Hz)
1.0
1.0
1.0
1.0
1.0
4.96
3.54
2.61
2.27
1.90
5.95
4.25
3.13
2.72
2.28
D(u.25 Hz)
3.20
2.50
2.05
1.88
1. 70
Note that the factors given for the construction of the vertical
spectra are to be applied to the horizontal components of
acceleration and displacement.
5.5.2 Site Dependent Spectra. The data base used for determining
the site independent spectrum shapes can be reexamined with regard
396
PERIOD, IN SECONDS
5.0
1.0
0.5
0.2
0.02
0.1
0.01
...
0
&oJ
'"
...'"
&oJ
'"
II::
&oJ
....&oJ
~
....
Z
...
&oJ
;!!;
>-
...!::
0
...J
&oJ
.
>
...J
.........
II::
....
'"
10
0.3L-____L -____
0.1
0.2
~~
0.5
20
FREQUENCY, IN HERTZ
Figure 5.
50
100
397
398
PERIOD. IN SECONDS
10
100
0.1
0.01
....
..........
Rock
0
<.>
II:
a::
..
...
;::
...z
t-
0.5
<.>
'!
>-
100
!:
<.>
0
...>
-'
...-'
II:
.......o
t-
10
EXPLANATION
E._blo mIG. pl
0.1
'lGftdord dl.latla. :
Saft-mldlum .oil.
o..p collllla.ll...on.
Stiff .oil
dip.,,,,
Rock
Revul.tory Guidi I 60
0.5~____~~~~______~~-L
0.01
____~__~~________~~
'0.1
10
FREQUENCY. IN HERTZ
Figure 6.
100
399
made.
When generalized parametric curves are judged to be
insufficiently precise, it might be desirable to utilize finite
element or similar procedures (24,25). The degree of accuracy with
which the physical characteristics of the earth materials at a site
are known should be evaluated critically before embarking upon a
costly site response analysis investigation.
6.
A CASE
Figure 7.
An
401
EL. 711.75
.--_----,-!L
.:r
11
-,1
II
I'
II
Ir-
I:
II
II
II
II
'I
693.00
II
=11
II
II
II
II
II
II
II
0
0
CO>
0
0
....
320
75 165
800
( b)
Section
Dimensions in cm
I
II
I 642.00
"""~"'~
(a)
Figure 8.
Elevation
402
than the end product of a broad optLmLzation effort aimed at
reducing the number of foundation excavations in difficult conditions and a more economical use of materials.
6.2
403
36'
37'
39'
38'
40'
4"
4"
40'
~.
"ERlURUM
---..............
Foult
Probable foult
Normal fault
~ Overthrust
_
Strike slip fault
-=:::;-
Figure 9.
404
405
100
~ (Hard
.. .....................
soils)
............
-.
...
......
,,"
.~"'z.:-..:::::...
.-~~ . . .".
"""
,',"
',.~.....
Rock
,.~
10
'".\>
.
',' ...
~\'.
\"
.
u
0
>
1.0
0.1
10
1.0
Frequency in
Figure 10.
Hz
406
comparison purposes two design spectra for "hard soils" and "deep
cohensionless" soils reported in (23) scaled to the appropriate
maximum acceleration are shown in Fig. 10. It is noted that the
general dynamic features of the soil columns are reflected in
these curves.
6.3
co~sidered
407
darn result from the action of a certain volume of water attached
to the surface (26). For
flexible structures interaction effects
are more complex but they are again approximated by an extension
of the added mass concept.
Earthquake induced vibrations of submerged structures can be
studied satisfactorily by assuming that wave action is negligible
and the velocity of the structure relative to the surrounding water
is low. Under these conditions the analysis can be performed by
adding to the mass of the structure the mass of
a certain volume of water which gives a new virtual mass but otherwise treating the structure as if it stood in vacuo. Classical
normal modes do not exist for the piers when hydrodynamic effects
are included, but their displacements in terms of their modes of
vibration in vacuo may still be used to approximate the actual
displacements.
Analysis of the response of structure interacting with fluids
is best carried out in the frequency domain. The response to unit
harmonic ground acceleration .
x.. ()
t
= Re {e iwt}
(22 )
may be expressed by
Y (t)
n
Re{Y (w)e
iwt
(23)
yew) = L(w)
(24)
~n
-:;;z
+ -322P
3z
32 P
= 2-2
(25 )
dt
408
409
0.5
direction
I -
'I - direction
0.4
:!!:::
u
.,
..,
.,
0.3
CIt
0.2
CD
0.1
o~----~------~----~------~----~
0 3 5
Period in I
Figure 11.
410
Sample Results
The design elastic forces for which the cross sections shown
in Fig. 8 were designed were derived from the spectra shown in
Fig. 11. It is of interest to note that good agreement was obtained
between the "enhanced" design spectrum applied to the structure
assumed to be standing in air and the reduced spectrum excluding
hydrodynamic forces applied to the piers with appropriately modified
mass distribution derived from the added mass concept to account
for the submerged portions. Although eight of the lowest modest were
considered the fundamental mode was by far the most dominant, and
results based upon it would have been quite satisfactory.
A sample of typical results is given in Fig. 12 for a 70 m
tall pier in terms of displacements, snear forces and b,ending moments
in the direction transverse to the bridge axis. Reinforcement
requirements were basically determined on the basis of working
stress procedures although the ultimate strength theory was used
in checking the strength against critical modes of failure. Contrary
to most design codes the computations were performed for the case
when half of the design train load was effective at t,le same time
as the design earthquake. This explains the finite value of ,the
shear force at the top of the pier in Fig. 12; it stems from the
mass of the truss and the train lumped at the top of the pier.
Because of these additively conservative assumptions and because
of the rather strong ground motion considered for design substantial
overturning moments were obtained for which the pile foundations
need to be designed.
7.
CONCLUSIONS
411
o
(a) Olsplacement In em
( b) Shear Force in t
30000
(c) Overturning Moment in tom
Figure 12.
412
REFERENCES
1.
2.
Iwasaki, T., J. Penzien, and R. W. Clough, "Literature SurveySeismic Effects on Highway Bridges," Earthquake Engineering
Research Center Report No. EERC 72-11, University of California,
Be rkeley, Novembe r, 1972.
3.
4.
5.
Wilson, E.L., A. der Kiureghian, and E.P. Bayo, "A Replacement for
the SRSS Method in Se ismic Analysis ," International Journal on
Earthquake Engineering and Structural Dynamics, Vol.9, pp. 187194, 19tH.
6.
"Norm fur die Belastungsannahmen und die Uberwachung der Buaten ,"
SIA 162, Art. 22, Schweizerischen Ingenieur und Architekten
Verein, Zurich, 1970.
~n
Finite
413
7.
d.
9.
!~ew
10.
11.
12.
13.
Glilkan, P., and M.S. Ylicemen, "Seismic Risk Analysis for Nuclear
Power Plants ," Journal of Pure and Applied Sciences, Middle East
Technical University, Ankara, April, 1977.
14.
15.
Schnabel, P.B., and H.B. Seed, "Accelerations in Rock for Earthquakes in the Western United States ," Bulletin of the
Seismological Society of America, Vol. 62, pp. 501-516, 1973.
16.
17.
414
18.
19.
20.
21.
22.
23.
24.
25.
26.
27.
2d.
Chopra, A.K., and N .M. Newmark, "AnaLysis ," in Design of Earthquake Resistant Structures, edited by E. Rosenblueth, Wiley,
198v.
415
E. ZELLER
Associate Professor
Department of Civil Engineering - Institut Universitaire
de Technologie - Universite STRASBOURG-III (FRANCE)
Introduction
In the history of construction. men have always had to
deal with the effects of wind. Any time a structure or construction
was unable to resist it, the wind always made a ruin of it. Each
time, man was given a lesson by the confrontation of his structural
concepts with reality.
Thus for ages, the knowledge of wind effects has
remained inseparable from pratical experience. A structure was
considered stable if it could resist ~ind ; thus, in many languages,
transversal stability was expressed in a similar fashion ...
"wi nd-braci ng", "contreventement", "controvento" ...
Certainly from the end of the last century, in dealing
with certain bold projects (the Eiffel Tower, skyscrapers, the
first long span bridges) research had to be done. However there
remained a world of difference between knowledge of structural
wind effects and the otherwise advanced science of aerodynamics
(aviation).
And so, in 1940, when the Tacoma bridge collapsed. other
similar bridges had already met the same fate. But in the case of
the Tacoma bridge~ this occurred during daylight hours. The
resonance effect was observable for a long time, and was even
recorded on film, as is well known. This served as a starting
point in the scientific investigation of the aeroelastic phenomenon.
416
Nowadays, in constructions or structures whose structural characteristics (dimensions, sections, mass and stiffness
distribution, damping ... ) stay in a certain defined range, national wind codes, although sometimes presenting differences, bring
ample security to classic constructions.
In the setting of the European Convention for Constructional Steelwork, a great amount of work has been achieved by
some twenty experts from ten nations This group has just
brought together the main strong points of various wind codes,
the results of wind-tunnel work and in situ measurement, such as
constructive experiments. Because the author has had the privilege of participating in the Technical Committee on Wind Effects
since 1968, the analytical approach in the beginning of his
report, chapter I, reflects the results of this work-group.
Nowadays, aeroelasticity has made notable strides, but
in direct consequence of the bridge sections being made of
sharp edges, the position and intensity of shedding vortices have
not been revealed by a faithful mathematical model.
That is why in great bridges with low frequencies, wind
tunnel study on complete or partial models remains the tool number
one for guaranteeing the required serviceability and security in
design.
When the construction of a bridge has been decided, the
theoretical analysis in relation to the parameters (length,
structure type, load, material ... ) leads inevitably to a solution
in which the two following aspects are sufficient to determine the
behaviour in the wind :
Point 1 - mechanical characteristics of the bridge:
rigidity, frequency, damping ...
Point 2 - a geometrical section determining the wind
flow conditions
When the corresponding aerolastic model leads us to a
structure whose response remains inferior to allowed limit values
in all cases of possible excitation, behaviour will be acceptable.
In the opposite case, when excitation conditions (different wind directions, speed, wind turbulence) produce flexural or
torsional instability of the bridge, the solution must be
modified by changing the frequencies or the section, until the
modified model yields a satisfying aerodynamic response.
The author will show the influence that the two above
aspects can have on actual bridges in chapter II of this report.
417
418
CONTENTS
Introduction
I.
2.
3.
4.
I I.
Wi nd tunne 1 approach
1.
2.
3.
III.
IV.
419
f..[:
where
a(n)
6(n) = Arc tg
[1/:
Jf
o
n of its
+ 6(n) ) dn
[1;
phase angle
(+
420
~+
$
:z
0
IU
t..I
a:
Q
\I}
....
01
II
~
~
-t.
~~
""d
S
~o
>
...
"
>
s:!
421
(m/s)2
2
1
l/n
Fig. 2
Wind spectrum measured by van den Hoven
In the aerodynamicsof structures the velocity fluctuation
of turbulent wind need to choose a base of time T of the order
of minutes and then the expression of wind speed becomes :
V = Vm + Vf
where
.!fv
Vm = T
Vf
dt
422
m,t,z
U-
(1 -
.!)
t ]
rL loge
(1 - P)
liN]
kt Vm,50,z with kt
t
t
200
100
YEARS 50
20
10
5
2
0.7
/
0.8
0.9
Fig. 3
/
,/
1.1
1.2
423
Wind direction:
As a general rule. the design speed Vmt z is
determined taking no account of direction. But in all cases the
wind is quasi-parallel to the ground surface.
Turbulence intensity :
The statistical analysis of the anemometer records shows
that the standard deviation 0v of the speed fluctuation is
f
= Vm, t ,z
f
-V--
m.t,z
Sv
(n)
where
S (n)
vf
x
= power
1800 n
m.t.10
424
iI,t,z = k log ~
e Zo
m,t,10
where
VII
k, z00 b and
at
z.o .. 2.5
~=
Fig. 4
425
Remark :
Vm,t,z does not take into account some phenomena such as
tornadoes, local wind, local topographic effect (escarpments,
hills, valleys).
The vertical gradient of wind velocity depends on the
terrain roughness of the five conventional roughness classes,
defined as:
sea or lake shore with more than 5 km of fetch
over the water
II flat terrain with some isolated obstacles (trees,
buildings, ... )
III rural areas with low buildings, trees, hedges
IV urban, industrial or forest areas
V
Roughness class
logartthmlc
law
II
III
IV
0,005
0,05
0,3
2,5
0,16
0,19
0,23
0,26
0,29
{'o
(m)
= _-..:..:K_ _
k log ~
e
zo
II
III
IV
0,16
0,19
0,22
0,23
0,23
426
(1)
Vf
r--~
~
;JY'
.l
(2)
Fig. 5
This variation can be described by a correlation
function p expressing the simu_ltaneity of occurrence -of gusts
at two points (1 and 2 ) :
427
=e
The parameters
relationship
42 ( z ) 0,25
20
"2
428
WIND
Fig. 6
When a slender structure oscillates in bending perpendicular to the wind, a periodic variation of pressure can be
produced and an amplification of vertical oscillation can lead
to aerodynamic instability if critical wind speed is reached.
-+
V wlnd
~~'d
V relative vorl ge
Fi g. 7
The torsion of the bridge section modifies the angle of
incidence; the resultant pressure variation can produce
instability in torsion.
WIND
Fig. 8
429
1T
oK
CL ~ S
Co
II
stable interaction
Fig.
aeroelastic instability
9
430
431
<:
Vd eSlgn
c
for a constant cross section
V = 1,63 0
c
with
n1
(J
n1~
x b
(J
Example
(J
=11-
2,5
_--I
2n ~ ~ 2 n ~ if ~ small
logarithmic decrement = ~
(~
% of critical damping)
the damping
. structural damping
Os
0.01
<
Os
<
0.06
0a
432
fl ex
rL 1
d 1 .1:1 tors .
1 fl ex
alon~
v'-)
c
r-__________~~d;~~~;-------------Ir---------433
Bri dge secti on
--cP-
-d
---i1r-----------~d
--~I 0.1-0.2
Lr-f
'I
0.2
0.1-0.2
0.2
_N
'20.5 d
- D-
0.3
0.3
0.3
0.15-0.Sd
I--
- c:::J
0.05
0.1-0.15
0.15 - 0.5d
...
-'
<.
0.13
0.11
.,
t-
-I
"1'e
-0) - {3d
......
:--
Fi g. 10
0.05
0.5
0.1
0.06-0.13
0.7
0.1-0.2
0.12
e:
I
---1
434
Remark 1
Remark 2
Remark 3
435
II.
436
2. Energy evolution/cycle
Fig. Ii
Fig. 12
A
B
C
D
E
F
mm fig
@'
I/"
10
30
40
mm
...
AMPLl TU DE
o
10
20
30
nvn
AMPLITUDE
/--_.
Vd
TORSION
= 9,6
= 3,9
-Vd
m:
Al1PLlTUDE
0tors
0nex
--- ----
HER
km/h
FLUTIER
245 kili/h
Fig. 13 - Dampino and Amrl i tude of the aeroelastic rrodel (for steel solution)
20
IJ
165 ~m/h
1,8 %
195 kmth
:;..-"'FLUHER
0tars = 3,7 %
'desiqn
10flex
Q(=+~.
--I---~--==~
20J.
30T J
42 ....
40
50
~LJ
Dynamic pressure
= 6,5 %
= 14,1 %
AMPLITUDE
0. ___
0flex
) DESIGN VdSPEED
.j:..
w
.....,
438
be~aviour
439
0+0
a
s
O2
0a+ 0s+ 0d
a = 0 0 s = 6,5 %
Fi g. 14
440
AMPLITUDE
em
10
1
I
BASIC BRIDGE
2,--1_~
SOFFIT PLATE
Fi g. 15
441
The initial damping 6 1 became increased with the temporary damping device 6d . With 02 the pseudo periodical movement
had now an amplitude ~ 3 cm .
Definitive solution :
The investigation on a model in a wind tunnel gave
clearly the same mechanical behaviour as observed on the bridge.
After exploring different shapes for the bridge as well as the
use of perfored girders (permeability modifies wind flow favorably) the definitive solution consisted in a simple external modification of the bridge deck : adding soffit plating. The effectiveness of different dimensions of fairing is illustrated in
Figure 15.
Importance of small details:
Furthermore, such insignificant details for a bridge designer, such as handrails, can considerably modify the amplitude of
the response when blocked by snow. The sectional model response of
the unmodified Long's Creek bridge allowed to reproduce this clearly.
Fig. 16.
Amplitude
on
18
16
!~
'I'
II:I I,
II
Ilt
12.
10
J \
:
I
I
Handrail clear
I
I
I
't
f \
I
2-
10
20
30
40
50
Fi g. 16
V
km/hr
I~B2m
19ot m
Fi g. 17
}V
630 m
,::
CO"f~2
0.52 m
L
.lQ
Am ~
Edge
i.:T]
tN
443
I CONF.
NUMUR
WIND
ANeLi
VORTEX SKtDDIN'C
VELOCITY, MPH
~IAXIMml
A~IPLITUDE
0
LOIo.'ER
UPPER
LI~lT
LIMIT
TORSION
179
171
+4
In
122
70
+6
29.0
43.7
+8
25.2
41.7
0.409
0.101
0.469
100
77
60
157
122
148
+6
26.4
37.3
0.261
+8
24.5
51.1
0.139
-8
26.3
44.3
1. 24
0.438
+4
+6
26.3
45.7
+8
26.3
41.1
106
147
110
+4
-6
1
VELOCITY
DECREASING
+2
-6
-8
VELOCITY
INCREASINC
0.335
0.468
80.2
In
152
III
102
155
153
132
126
99
88
-6
114
125
-8
124
117
Fig. 18
Limit
servic
Wind velocity
V
Fig. 19
11 y t he
For the brldge
des Meules,1/ curve I shows schematlca
solution of flexural and curve II that of torsional movement with
the wind speed. In fact,each response type is a combination of
response I and II. In order to improve structural safety,the
range V flexion V torsion has to be large ; bending and
cr
cr
torsion or coupled flutter have to be situated largely over the
design speed.
445
Bendi~g
I'llIft
incidence
= + 4
vm = 8 m/s
Cl
V = Vm + Vf = 17,3 m/s
2.0
,,
,0
,,
\'\
11
IS
Damping
Fig. 20
~6
3. Aerodynamic stability of a
brid~e
under erection
447
REFERENCES
CHAPTER I
E.C.C.S. - Recommendations for the calculation of wind effects
on buildings and structures, 1978.
SCRUTON C., NEWBERRY C. - On the estimation of wind loads for
Buildings and structural design, 1963.
SFINTESCO D., ZELLER E. - Prise en compte des effects du vent
sur les constructions, 1975.
VAN KOTEN - Windbelasting, TNO, 1974.
WYATT T.A. - A note on the significance of Cross-Modal
Correlations on the Evaluation of quasi-static
response to wind, 1974.
CHAPTER II
DAVENPORT A.G.
CHAPTER II I
BORGES J.F., CASTANHETA M., BORGES A.R.J. - Design criteria for
wind loads on statistical Bases,
Proceedings
Wind effects on Buildings and structures, Tokyo,
1971.
WARDLAW R.L. - Approaches to the suppression of wind induced
vibrations on structures, 1979.
Laboratory report on the aerodynamic stability of
the Pasco-Kennewick intercity bridge, 1974.
449
LIST OF PARTICIPANTS
Mr. K. Ak1ll1lar (Turkey)
Dr. N.
Akk~
(Turkey)
450
YargL~oglu
(Turkey)