Grillage Analogy C.S.surana R.agrawal - For Word
Grillage Analogy C.S.surana R.agrawal - For Word
Grillage Analogy C.S.surana R.agrawal - For Word
C.S. Surana
IL Agrawal
-
Narosa
Z-
1-1-Toic
-L)%%
,
C.S. Surana
Department of Civil Engineering
Indian Institute of Technology,
Delhi Nevi Delhi-110 016, India
R. Agrawal
Department of Civil Engineering
Institute of Technology, Banaras Hindu University
Varanasi, India
.r
Copyright C 1998 Narosa Publishing House
NAROSA PUBLISHING HOUSE
6 Community Centre, Panchsheel Park, New Delhi 110 017
35-36 Greams Road, Thousand Lights, Madras 600 006 PIA
306 Shiv Centre, D.B.C. Sector 17, K.U. Bazar P.O., New Bombay 400 705 252 2
2F-2G Shivam Chambers, 53 Syed Amir Ali Avenue, Calcutta 700 019
3 Henrietta Street, Covent Garden, London WC2E 8LU, UK
ISBN 81-7319-153-0
_ To
Our Parents
for having shown us the way
Our Families
for their constant support
A
4
V
Preface
Bridge design and construction all over the 'world has undergone remarkable
changes in the past two decades. The increased demand for complex
roadway alignments, advances in construction technology and availability of
computing power for bridge design, are some of the factors for these
developments.
Over the years, a number of methods of analysis of highway bridge
decks have been- evolved. The methods range from the simplified hand
method like Courbon's or graphical methods like Hendry-Jaeger and
Morice & Little etc. to highly sophisticated methods like finite element,
finite strip etc. The former are conservative and the latter which require
fairly complex computer programs and larger computational facilities, are
prone to errors of idealization and interpretation of results.
Grillage analogy method, which is well-established and computer-
oriented, bridges the gap between the two. The present book describes
bridge deck analysis by grillage analogy. The method is versatile in nature
and can be applied to a variety of bridge decks having both simple as well
as complex configurations with ease and confidence. Analysis of bridge
decks employing grillage analogy is possible on commonly available PCs
while retaining the accuracy and versatility of other refined methods that
l y usually require larger computational facilities. A considerable saving in
time in the analysis of the bridge is achieved and the method also provides
a 'feel' of the bridge, behaviour to the designer. Although, the accuracy of
any method of analysis for a particular structure is difficult to predict, the
method of grillage analogy is found to be fairly accurate when compared
to methods like FEM.
The book is mainly intended for professionals and students. Consultants
and researchers who are confronted with the problem of analysis for bridge
design and those who wish to specialize in the subject, will also find it
viii Preface
A I L
Preface ix
is so designed that loads other than Indian standard loadings, can also be
easily incorporated. The listing of the program is given in Appendix II.
A number of worked out examples of different types of slab, T-beam and
box-girder bridges are given in Chapter 7 to explain the use of the program.
The example are chosen from actual life-size bridges and the solutions are
obtained for IRC loadings and also for user's specified loadings. Input and
selected output modules are given for the convenience of the user.
Exhaustive and relevant references are included after each chapter for
the benefit of the readers.
A diskette containing two programs of the Appendices (in a ready to
use form) can be ordered from the authors through the publisher. The
diskette also contains the input data and exhaustive force responses/output
of all the worked out examples given in the book.
Although the general methods and concepts postulated by Lightfoot, West,
Hambly, Jaeger, Bakht and others are further developed and subsequently
expanded, we are deeply indebted to them and to many others whose
writings, teachings or personal help have shaped our thinking and approach to
the subject matter.
The financial assistance and other help rendered by the Curriculum
Development Cell of I.I.T. Delhi towards writing the book is gratefully
acknowledged. The sabbatical leaves granted to the authors by I.I.T. Delhi
and I.T. BHU, Varanasi, respectively for this joint venture are also thank-
fully acknowledged. We are indebted to our research scholars who helped
in developing and checking the computer programs and to our numerous
students for assistance in worked examples and proof checking. The works
of preparing drawings and typing of the manuscript was done by different
persons at different times. We acknowledge their help.
Any suggestions for further improvement from the readers would be
very much appreciated.
C.S.SURANA
RAMJIAGRAWAL
Contents
Preface
1. Introduction
L1 General 1
1.2 Recent Trends in Analysis and Design of Bridges 2
1.2.1 Structural Systems 2
1.2.2 Computer-Aided Methods of Analysis 3
1.23 Design Methodology 4
1.2.4 Modern Construction Techniques 4
1.3 Structural Forms of Bridge Decks 5
1.4 Form of Construction 6
1.4.1 Slab Bridge 6
1.4.2 Slab-on-Girders Bridge 9
1.4.3 Box-Girder Bridge 12
1.5 Plan Geometry or Planforms 15
1.6 Support Configurations 16
1.7 Bridge Loadings 16
1.7.1 Loading Requirements 18
1.7.2 Dead Loads 18
1.7.3 Live Loads 18
1.7.4 Impact Loads 25
1.7.5 Footway, Kerb, Railing and Parapet Live Loads 26
1.8 Comments on Loading Standards 28
1.9 Organisation of the Text 30
References 33
Introduction
1.1 , GENERAL
Bridge construction has been one of the important engagements of mankind
from the earliest days. Bridges are one of the most challenging of all civil
engineering works. It has always fired the imagination of people as they
seem to lead to hitherto uncharted territory.
Bridge construction tqday has achieved a world-wide level of importance.
The nuiribers and sizes of bridges have continuously increased in the last
fifty years. Man's increasing mobility through railway and motorised
transport has caused such complex forms of bridges to be built, which had
seemed unrealistic earlier. To cope-up with this demand, tremendous efforts
all over the world in the form of active research in analysis, design and
construction of bridges is continuing.
Over the years, a number of methods of analysis of bridge superstruc-
tures have been evolved and are being used. Courbon's method, Hendry-
Jaeger method and Morice and Little method are some of the methods
which have been in use since long, and, are still popular, as they are found
to be easy, amenable to design graphs and also reasonably accurate for
bridge decks of simple configurations. But these methods are being gradu-
ally replaced where computer facilities are available or more accurate
analysis is desired or the cross-section and/or layouts of the bridge decks
are complex.
Following the advent of digital computers, computer-aided methods like
Finite Element, Finite Difference, Finite Strip have been developed and are
in use to analyse intricate forms of skew, curved, bifurcated and arbitrary
shapes of bridges having usual support conditions and cross-sections. But
these methods are highly numerical and always carry a heavy cost-penalty_
Grillage Analogy is probably one of the most popular computer-aided
methods for analysing bridge decks. The method consists of representing the
actual decking system of the bridge by an equivalent grillage of beams. The
dispersed bending and torsional stiffnesses of the decking system are
assumed, for the purpose of analysis, to be concentrated in these beams. The
2 Grillage Analogy in Bridge Deck Analysis
stiffnesses of the beams are chosen so that the prototype bridge deck and
the equivalent grillage of beams, are subjected to identical deformations
under loading. The actual deck loading is replaced by an equivalent nodal
loading. The method is applicable to bridge decks with simple as well as
complex configurations with almost the same ease and confidence. The
method is easy to comprehend and use. The analysis is relatively inexpen-
sive and has been proved to be reliably accurate for a wide variety of
bridges. The grillage representation helps in giving the designer a feel of
the structural behaviour of the bridge and the manner in which the loading
is distributed and ,eventually taken to the supports.
The book essentially deals with the Grillage Analogy method and its
applications to variety of bridge decksboth simple and complex. But it is
also thought relevant to introduce in brief other existing methods of analysis
to the readers (Chapter 2).
In order to apprise the readers with the developments in bridge engineer-
ing, the present chapter discusses The recent trends in analysis and design of
bridge decks. The structural forms of decks based on types of construction,
planforms and support conditions have been outlined. The loading standards
used in India for concrete highway bridges, are given in some detail. Indian
Roads Congress (IRC) loadings are compared with the standard loadings
adopted by some other developed countries.
___________________________ E
/
a) Solid Slab
______________ E
a
b) Solid Stab
____________________________ fl
1 7:7000007 - -
c) Voided Stab
Figure LI Slab Bridges
stressed concrete slabs of depth upto 1.2 m are cheaper. For moderate skew
crossings having spans of 15 to 25 m, this type of deck with longitudinal
prestressing is useful but for highly skewed crossings, reinforced concrete -
decks are preferred for ease of construction. If the voided section is found
inadequate in shear, it should be kept solid near supports (refer voided slab
example of Chapter 7).
In R.C. slab bridges, span-depth ratio of 15 for simple spans and 20 to 25
for continuous spans are usually adopted for both solid and voided slabs.
For cast-in-situ, prestressed concrete voided slab bridges, this ratio is nearly
30. In precast prestressed voided slabs, the ratio ranges between 25 and 30.
The deck slab overhang, designated as 'a' in Figs. 1.1b and 1.1c may be
provided to produce the desirable aesthetic effect and also to reduce the
dead load and the width of sub-structure.
In many countries, standard precast prestressed beams are employed for
short- and medium-span bridges. These standard beams are closely
positioned across the width of the bridge and in-situ concrete is poured to
give transverse connection in order to create a slab-type deck. This form of
deck is described as Pseudo-slabs. Such type of two-stage casting is also
referred to as Contiguous Construction and the standard beams are termed as
Contiguous Beams.
8 Grillage Analogy in Bridge Deck Analysis
Shear key 7
a) inverted T-Beams Bond rod
. i . . J. ..
. . . :
r .
. . --. . .c
.
r. 1 .
..1
\
, _,., .
1... ...- '..:.
. . ; :
.....
:. Z.:. . % .
Il
..
Void
b) Box -Beams
with in-situ concrete to form an integral slab. Transverse bond rods are
provided just above the bottom flange passing through pre-formed holes in
the precast beams. This type of construction can be used upto 20 m span.
Another standard group of beams, known as M-Beams have been developed
by Cement and Concrete Association (C&CA), U.K., in collaboration with
the Ministry of Transport, U.K. (Fig. 1.3) for short- and medium-span
bridges.
Pre-tensioned cast-in-situ Pseudo-slabs are usually adopted when the
erection of formwork presents no difficulty. The advantage of this form of
construction is that the structure is monolithic and the stress-distribution in
the slab can easily be evaluated. On the other hand, the precast zire.crs v:ith
in-situ concrete filled up, are preferred when there is difficulty in supporting
Introduction 9
b) T-Beams
Figure 1.3 C & CA, U.K., Standard M-Beams
1
1 1
I
I b
1 I
Cross ben
h
Figure 1.4 T-Beam Bridge
simple spans and 12 to 15 for continuous spans. Higher ratios are possible,
but riding qualities are affected by creep characteristics of concrete. The
girder spacing 'h' (Fig. 1.4) may vary as justified by comparing the cost of
corresponding slab thickness. The usual range of spacing 'h' is between 2 to 3
-
m for these bridges. The stem-width 'b' is about 300 mm. This stem width is
increased to 'B' at the bottom, forming a bulb to accommodate a large
number of reinforcement bars there. This '13' may be kept between 500 to
600 mm. The stem width '6' is increased to '13' in the end-region, to take care
of large shears occuring there.
Slab-on-girders bridge also includes prestressed concrete bridges. Majority
of prestressed concrete bridges, constructed in India are of post-tensioned
type. The bridge decks, with post-tensioned girders suitable for simply
supported construction, may ..have either fully cast-in-situ slab and girders or
deck slab with precast prestressed girders alongwith cross-beams, assembled
together and transversely prestressed (Fig. 1.5a,b). Such types of construction
is convenient for the span range of 20 to 30 m. The span-depth ratio is
usually kept as 20 for simple spans and 25 for continuous spans for
prestressed concrete girder bridges. The girder spacing 'h' (Fig. 1.5a) is
normally kept between 2.0 and 4.5 m. The stem width 'b' should preferably
be a minimum of 300 mm to facilitate prestressing of tendons. The deck slab
overhang 'a' should be provided as required to produce the desirable aesthetic
effect and to reduce transverse moments.
Decks with composite construction are also popular for short- to medium-
span bridges. Composite construction refers to the use of structural elements
made of two materials in combination in such a way that they act
Introduction
11.
Cross ,ocb
beam
a 'I
h _I, h
11111111111111Nal111111enIF
Precast girder
Precast cross- beam
girder
Precast
Cast-in-
situ slab
Steel girder
p
t w
S
V./
4 h ___ I
(b) Multi - Celt Box-Girder
T
5
I.-- a
(c) Single-Celt Trapezoidal Box-Girder
SI b w o -
a 4. b h
+- b
Figure 1.8 Multi-Cell Spine Box-Girders
The span-depth ratio for R.C. box-girder bridges are generally adopted as
16 for simple spans and 18 for continuous spans. For prestressed cast-in-situ
concrete box girders, this ratio ranges from 20 to 25 depending upon support
conditions. In case of precast prestressed box-girder bridges, the span-depth
ratio is taken between 18 and 20 whereas for spine or spaced box-girders, the
ratio lies usually between 18 and 22. In precast, prestressed multi-cell box-
beams, the ratio can be as high as 25 to 30 [9].
The spacing 'h' of box-girders usually lies between 2.0 and 3.5 m. The web
of box-girder superstructure should have a minimum thickness `ki: of 200
mm. It is often useful to increase the thickness of webs near supports to
provide adequate concrete section for shear resistance. Precast box-beams
ordinarily have a width '17' of 1.0 to 1.2 m and height 's' in the range of 0.6 to
1.2 m (Figs. 1.7 and 1.8). The bottom slab thickness 'd,' is kept
Introduction IS
approximately 1/20 of clear span between webs but it should not be less
than 150 mm and may be increased near continuous supports.
1.5 PLAN GEOMETRY OR PLANFORMS
The horizontal and vertical alignments of a bridge are governed by the
geometries of the highway, roadway or channel it crosses. A bridge may
either be right or skew, straight or curved or any combination thereof.
Typical geometrical planforms of bridge decks are illustrated in Fig. 1.9.
Suspended span
(c) Balance
d Cantilever
Arrangemen
t
fi
ft 11
-
(d) Continuous Arrangement
Figure 1.10 Support Configurations of Bridge Decks
dinal distance between two successive trains is 18.5 m. The minimum centre
line distance of the wheel-line from the edge of the Kerb works out to 400
mm. The configuration of the load as well as the position of each wheel is
given in Figs. 1.11(a) and 1.11(b).
1200 1200
1-- 8300 f--4800 H-4800 --1
X
4- + -F. + + + +
Direction of Motion
1800
+ +
(Left Most
Front Wheel
H81-320044-1.300 +3000+3000+3000-
1 (b) Plan
1800
(c) Section
on x x
(411. Dimension are in mm)
Figure 1.11 (a) 1RC Class A and Class B Loadings
Class A-two laneS loading consists of two class A -single lane trains placed
side by side at specified minimum clearance. Class A loading is
20 Grillage Analogy in Bridge Deck Analysis
Clear Roadway
3600
_____ 7200 _____
0) Elevation
4- 1-1-++ ++4- + + -r
Suggested Equivalent
Concentrated Loadings o
un
O
Left Most
Front WheeNt- 1-4--++++++1
[4-9@360z3240 ______
(20 toads 3.5t =70t on Two Tracks)
(ii) Nan
35t 35t
2900
(iii) Cross-Section
(All Dimensions are in mm)
Figure 1.12 (a) IRC Class AA Tracked Vehicle
iv) IRC Class 70R Loading
This is the revised version of Class AA loading and consists of tracked
.and wheel loadings. The loading is detailed in Appendix I of IRC Code
[191 The minimum clearance between the road face of the kerb and the
outer edge of the track or wheel is same as for Class AA loading. The
spacing between successive vehicles is 30 m. This spacing is measured
from the rear most point of ground contact of the leading vehicle to the
forward most point of ground contact of the following vehicle.
22 Grillage Analogy in Bridge Deck Analysis
1 1 1 ___ I
375t 6.25 t 6.25 t 3.75 t
(1) Cross Section
-I-
+Mos +
t---
Lett -
10,2
Front Wheel
Ill\
+
-I- + 4. + I-
k 600 4*
__ 1000 4-
600 --1
2200 ___________
) Plan
Carriage Way Minimum
Width Value of C
Single Lane Bridges
3800 and 300
above
Multi L ane B ridg es
- 5500 600
5500 1200
All Dimensions are in mm
Figure 1.12 (b) IRC Class AA Wheeled Vehicle
70R track loading, as before, weighs 70 tonnes (Fig. 1.13). The track
dimensions are slightly different than those of Class AA track loading. For
design purposes, wherever required, each strip loading could be idealised into
a suitable number of point loads say 8 or 10.
70R wheel loading is of two types: (1) 70R Bogie loading weighing 40
tonnes through two axles each weighing 20 tonnes (Fig. 1.14a) and (2) 70R
train loading weighing 100 tonnes through seven axles, one axle of 8 tonnes,
two axles of 12 tonnes each and four axles of 17 tonnes each (Fig. 1.14b).
Introduction 23
4570
____ 7920 _____
(a) Elevation
t ,+ + + + + + + + -i-
vSuggested Equivalent I
Concentrated Loadings c,
4 + 4 +
Left Most / 1-P-96)457= 4113.1
35T. 35T
14 2900
An axle may have four or eight wheels on it. There are two, four wheel
arrangements and one, eight wheel arrangement leading to three alternate
wheel arrangements termed as Col. 'I', Col. 'm' and Col. 'n' arrangements
[19]. AU axles will have the same arrangement of wheels at a time and all
wheels on an axle will have equal loads. The two alternate four wheel
arrangements namely Col. '1' and Col. 'm' are given in Figs. 1.14(a, b). The
eight wheel arrangement namely Col. 'n' is not found critical and is not given
nere. However, if required, details of Col. 'n' arrangement could be obtained
from IRC Code [19].
24 Grillage Analogy in Bridge Deck Analysis
I- 4880
20t 20t
(i) Elevation
" --"F
4
0
+
I
0 JCI
CO
ci
-I-
Left Mos +t
Front Wheel (ii) Plan
1-e--
n nb
NI
(iii) Cross-section
Note: (i) Min. ro' same as for Class A A Loading
(ii) Max. 'w' =410
(iii) Either a=450 & b =1480 (CoI.1')
or a=795 & b = 790 (CoL`m1);
so that (2a4b) = 2380
Figure 1.14 (a) IRC Class 70R Bogie Loading
IRC Code [19] also gives in Appendix I certain other types of track and
wheeled loadings. These are lighter than Class 70R tracked and wheeled
loadings discussed above. These are to be adopted if a user specifies these for
the bridge.
For detailed loading Standards and their specifications, the reader is ad-
vised to consult the relevant IRC Codes of Practice mentioned under refer-
ences [18, 19].
The above IRC loads are nlared on the. hriricre clerk anr1 mnveri tong;til_ dinally as
well as transversely in small increments to cover the entire bridge
Introduction 25
1
a -
O T
3960
213 137073050 7'713 0 -
1f
1 t 12t 1 t 1 7 1 17t
t P
(i) Elevation
+.
44 +
++++ +4-
00.4co
4 4
(ii) Plan
(Left Most
Front Wheel
Note : Cross-Section and distances are as for
Bogie Loading shown in Fig. 1.14(a)
(All dimension are in mm)
Figure 1.14 (b) IRC Class 70R Wheel Loading
deck. One Class A or Class B loading can be put on every lane of the
roadway of a bridge. For multi-lane bridges, one lane of Class AA or Class
70 R per two lanes of the carriageway is allowed as an alternate to Class A
loading.
the bridge span. This percentage can be calculated using suitable formulae
[19] or could be directly read from Figs. 1.15 (a,b) for both steel and concrete
bridges for different types of IRC loadings.
55 54.5
50
74 . 45
L= Span in m
Lo 40
35
900
30 44 13.5+L
0
25 Steel
Q. 20 450
E 6 +L Concrete
15.4
15
10 8.8
5 I 1 I
0 5 10 15 20 25 30 35 40 45 50
Span (m)
Class A and Class B Loadings
Figure 1.15 (a) Impact Percentage Curves
25
O .
4.
20 0 15
u10
"" 5 0 5 10 15 20
30 35 40 45 50
(a)
25
Span (m)
Class AA and Class 70R Wheeled Loading
.7.30
25
025C
2:20
10 STEEL
40 45 50
.7.30
L.= Span in m
025C
Concrete & Steel
4
(
2:20
215 J.
4-
10
U CON C. 8.6
o. 5
F 0 5 10 15 20 25 .30 35
Span (m)
(b) Class AA and Class 70R Tracked Loading
Figure 1.15 (b) Impact Percentage Curves
= ( p 2 6 0 + 4800)(16.5
L A 15 )
where P` = 400 kg/m2 or 500 kg/m2 as the case may be based on
(i) above.
P = the live load in kg/m2-
L= the effective span of the main girder in m.
W = width of the footway in m.
(iv) . Each part of the footway shall be capable of carrying a wheel load
of 4 tonnes, which shall be deemed to include impact, distributed
over a contact area. 300 mm in diameter; the permissible working
stresses shall be increased by 25 per cent to meet this provision.
(v) The railings or parapets shall be designed to resist a lateral horizontal
force and a vertical force each of 150 kg/m applied simultaneously at
the top of the railing or parapet. These forces need not be considered in
the design of the main structural members if footpaths are provided. In
cases where footpaths are not provided, the effect of these forces shall
be considered in the design of the structural system supporting the
railings and the footpath upto the face of the footpath kerb only.
64
56
E 24
-43 16
c
gel 8
0 20 40 60 80 100
span(m)
(a) Maximum bending moment for two lanes
280 Legend :
240 New
zealand Japan
200
Shear force 00
---France
160 --.---West Germany
--IRC (India)
120
---A ASHO group -
80 -)BS (HA group)
----Sweden
40
0 20 40 60 80 100
Span (m)
(b) Maximum shear force for two lanes
Figure 1.16 Comparison of Highway Bridge Loadings
REFERENCES
1. BAKHT, B. and JAEGER, L.G., "Bridge Analysis Simplified", McGraw Hill Book
Co., New York, 1985.
2. CUSENS, A.R. and PAMA, R.P., "Bridge Deck Analysis", John Wiley and Sons,
1975.
3. HAMBLY, E.C., "Bridge Deck Behaviour", Chapman and Hall, London, 1976.
4. HARVEY, DAVID [., "Recent Trends in Short and Medium Span Highway
Bridges in the United Kingdom", Proc. Intl. Conf. on Short and Medium Span
Bridges, Vol. I, 1982.
5. KULKA, E and LIN, T.Y., "Comparative Studies of Medium Span Box-Girder
Bridges with Other Precast Systems", Proc. Intl. Conf. on Short and Medium
Span Bridges, Vol. I, 1982.
6. LEE, D.J., "Progress in Bridge Engineering", Developments in Structural
Engineering, Proc. Fourth Rail Bridge Centenary Conference, Vol. I, 1990.
7. MATHIVAT, IE., "Recent Developments in Prestressed Concrete Bridges", FLP
Notes, 1988.
8. Pormuswamv, S., "Bridge Engineering", Tata McGraw-Hill Publishing Co.
Ltd., New Delhi, 1986.
9. RAINA, V.K., "Concrete Bridge Practice-Construction, Maintenance and
Rehabilitation", Tata McGraw-Hill Publishing. Co. Ltd., New Delhi, 1988.
10. ROWE, R.E. and Somenviu.e, G., "Research on Slab Type and Spine Beam
Bridge", Proc. Intl. Conf. on Developments in. Bridge Design and
Cnnstruction. Crosby Lockwood & Son Ltd., 1971.
11. STEINMAN, D.B. and WATSON, S.R., "Bridges and their Builders", Dover
Publications, New York, 1957.
12. THOMAS, P.K., "A Comparative Study of Highway Bridge Loadings in Different
Countries", U.K. Transport and Road Research Laboratory-Supplementary Re-
port 135 UC, 1975.
13. VicroR, D.J., "Essentials of Bridge Engineering", Oxford & IBH Publishing
Co., New Delhi, 1980.
14. WIELAND, M., "Modern Bridge Engineering in Structural Concrete", Proc. of
Asia-Pacific Conference on Road, Highway and Bridge Maintenance and
Rehabilitation, Bangkok, 1987.
15. WrreFoirr, H., "Building Bridges", Bentan Verlog GmbH, Germany, 1984.
16. "Specifications for Highway Bridges", Pt. I and Pt. III, Japan Road
Association, 1984.
17. "Bridge Rules Specifying the Loads for Designing the Super and Sub-
structures of Bridges and for Assessing the Strength of Existing Bridges",
Govt. of. India Publication, Ministry of Railways, New Delhi, 1977.
18. IRC:5-1985, "Standard Specifications and Code of Practice for Road Bridges,
Section I-General Features of Design", Indian Roads Congress, New Delhi,
I99A.
34 Grillage Analogy in Bridge Deck Analysis
19. IRC:6-I987, "Standard Specifications and Code of Practice for Bridges, Section
HLoads and Stresses", Indian Roads Congress, New Delhi, 1990.
20. IRC:18-1985. "Design Criteria for Prestressed Concrete Road Bridges (Post-
tensicined Concrete)", Indian Roads Congress. New Delhi, 1990.
21. IS:1343-1981, "Indian Standard Code of Practice for Prestressed Concrete",
Bureau of Indian Standards, New. Delhi, 1981.
Chapter 2
Methods of Bridge
Deck Analysis
2.1 INTRODUCTION
Extensive research into the behaviour of bridge decks under loading had
been carried out pver the past five decades and a number of methods of
bridge deck analysis were evolved from time to time. Prior to the general
use of the computer-aided analysis, hand computation methods and charts
based on some approximations and idealizations, provided convenient
methods of load distribution. These were reasonably accurate for design
purposes. However, with the advent of digital computers, many computer -
aided methods have been developed and are in use although some of these
methods are highly numerical and expensive.
The present chapter aims at, giving a broad idea about the various impor-
tant methods of bridge deck analysis. The applicability of the different
methods in relation to bridge type, plan geometry and support conditions is
also presented in a tabular form.
of Inertia
I .0-1 1
;
Courbon's Method * *
Orthotropic Plate * * * * *
Theory
Finite Difference * * *
Method
Analysis
Grillage Analogy * * * * * * * * * * *
Folded Plate * * * *
Method
Finite Element * * * * * * * * *
Method
Finite Strip Method * * * * * * * *
38 Grillage Analogy in Bridge Deck Analysis
will be maximum in the exterior girder on the side of the eccentric load (or
c.g. of loads, if there is a system of concentrated loads) and minimum on
the other exterior girder.
The load R1 is given by
PI (PLedEL)
R. =L "
EI S E1 , z 4
or,
(2.1)
R.=--P4(1--L1 edij
Eli cf
where I. = moment of inertia of ith longitudinal girder P
= total live load
e = eccentricity of live load (or- in case of multiple loads,
distance of c.g. of loads from centroid of moment of inertias) d.
= distance of ith girder from centroid of moment of inertia
For a N-girder bridge with. all girders
having same moment of inertia, the above
formula reduces to (2.2)
D2'd + 2 H 4 4
+D, d W
P(x,y) (2.3)
d xw 4 d xd 2 y2 1../ Y d y4
where 21-/ = (Dxy, DD + DI + D2)
Dx and D - are the equivalent flexural rigidities and Dxy and D are the
equivalent torsional rigidities per unit width in longitudinal and transverse
directions respectively. Di and D2 are the coupling rigidities per unit width
arising due to Poisson's ratio effects.
The flexural and torsional rigidities have significant influence on the
load distribution. Their effect is considered through two dimensionless
characterising parameters namely flexural parameter 9 and torsional para-
meter a as given below : b Ds 2a
(2.4)
Dy. 30.25
+ D3w+ Di+ D2
and a = _________________________________ (2.5)
2 11 (px
where, 2a and 2b represent span and width of the equivalent plate respec-
tively. The coupling rigidities D and D2 are small and can be neglected
without introducing any significant error.
The dimensionless distribution coefficients, which are given for nine stan-
dard reference points and load positions across the bridge width, are plotted
40 Grillage Analogy in Bridge Deck Analysis
or tabulated against values of O. The charts or tables are given for two values
of a, namely a = 0 and a= 1. Values of the coefficient Ka for any interme-
diate value of a are obtained by the following interpolation function,
section are computed as usual and pose no problem. The torsional rigidity
of the section is evaluated based on St. Venant's method by taking the
summation of the torsional rigidity of the components forming the section
and from this torsional rigidities D and D of the section are computed.
yx
However, the orthotropic plate idealization does not always represent the
physical behaviour of the bridge decks. In bridge with few girders, say 3,
the bending moments obtained are subject to errors, especially if the
bridge is wide and load occupies only a fraction of the width. Also, the
transverse moments which are complex combination of bending between
girders and bending due to non-uniform girder deflections, can not be
accurately obtained.
Also, the method suffers from the drawback of having to assume the deck
to be uniform throughout and also the design curves involve a certain degree
of approximation in use. As only the first harmonic component is used in
load distribution, the method is not recommended to be used- to estimate
shear. Further, the method can not handle skew bridges.
'-aw
Methods of Bridge Deck Analysis 41
difficult to apply. The finite difference method is the answer for such complex
boundary conditions. The method is versatile in nature and has wide
applications:-
The finite difference technique had been used to advantage, first by
Neilsen [34], later by Westergaard [49] and applied to bridge decks by
Naruoka and Ohmura [33], Heins and Hails [17], Robinson [38] and many
others.
In this method of analysis, the deck is notionally divided into grids of
arbitrary mesh size and the deflection values at the grid points are treated as
unknown quantities. The usual governing differential equation of an
orthotropic plate is considered in the finite difference method. The differen-
tial equation and accompanying boundary conditions are expressed in terms
of these unknown deflections. The resulting sets of linear simultaneous
equations are then solved for these unknown deflections. Finally, moments
and shear forces are determined from the known deflection pattern. The
curved deflection profile of the deck is approximated by a series of straight
lines and, naturally, accurate results can be expected only if fine grids are
used.
Finite difference equations for various boundary conditiOns like
simple, fixed, free or a combination of free and simple supports can be
written down for each case. However, a fixed edge condition is not treated
very accurately by this method except when very fine grids are used. The
simultaneous equations formulated from these grid points are solved on
digital computers which have matrix packages specially designed for these
type of problems.
The presence of interconnecting beams below the slab, present a special
boundary situation. A better representation of slab and beam interaction can
be found by treating the two as separate structural elements which are made
compatible by satisfying a set of boundary conditions. The interaction of
beams and the top deck is based on the assumptions that beam and deck slab
have identical deflections, the beams have no torsional rigidity and there is no
horizontal force between the beams and the slab.
The simplicity of the trigonometric solution for simply supported right
bridge decks tends to be lost when extended to the problem of skew decks.
However, the finite difference method has been used extensively for skew
bridges also where the grid may be taken along the orthogonal coordinates as
shown by Jensen [22], Robinson [38] or along skew coordinates as suggested
by Naruoka and Ohmura [33].
Most of the solutions and published values of deflections and moments for
skew necks roicr oni:yr to isotropic plates; the mesh is court4.... and acc
42 Grillage Analogy in Bridge Deck Analysis
racy is doubtful. Naruoka and Ohmura [33] neglected Poisson's ratio. They
dealt with torsionless plate (H = 0) as well as isotropic plate (H = Dx= Dy).
They had difficulty in setting up satisfactory interpolation equations between
these two limiting cases.
.favour [21] produced influence surfaces for deflections and moments for a
wide range of skew orthotropic slabs for both uniform and concentrated loads
making various over-simplifying assumptions. Schleicher- and Wegener [41]
published tables of deflections and stresses for continuous isotropic skew
slabs under uniform loading. Ghali [13] used finite difference equations at
different segments of a_ two-girder skew bridge to determine influence
coefficients for deflection; shear and moment assuming that the applied loads
act directly on girders leading to erroneous conclusions. Robinson [38]
assumed that a concentrated load might reasonably be distributed over one
mesh quadrilateral. Paterson [35] adopted this assumption in developing an
ALGOL computer program based on standard skew mesh over the skew
orthotropic plates.
Although, finite difference method, applied to skew bridge decks, has been
able to solve a large number of bridge structures but still it has certain
inherent difficulties in its use like adoption of fine mesh work, proper as-
sessment of singularities around the obtuse corner, deterioration in conver-
gence with increasing skew 'angles and in some cases, non-compatibility in
boundary conditions. The number of parameters involved in skew bridges are
such that the preparation of design curves to cover a realistic range of
loadings, skew angles and degree of orthotropy, does not appear to be a
practical proposition.
The finite difference method has also been extended to bridge decks curved
in plan by Heins and Hails [17] based on governing equation of orthotropic
plate in polar coordinates, neglecting Poisson's ratio. However, the efforts
needed to analyse a curved deck by this method is cumbersome
and is generally not recommended.
.-
2.6 METHOD OF HARMONIC ANALYSIS
In harmonic analysis, the applied load is broken into a number of harmonic ,
components, each consisting of a distributed load parallel to the longitudinal
axis of the structure and with intensity varying as pure sine-wave as shown in
Fig. 2.2.
11 It X
Mathematically, P(x) a sin ___ (2.7)
,t=3 L
Methods of Bridge Deck Analysis 43
12 (LI3 Eir a
(2.8)
= 7r4 71 El
'= 2 h CJ (2.9)
The extreme value of /3 for torsionally stiff bridges has been assumed as
infinity. Design graphs between a and p for 13 = cc have also been
plotted similar to graphs for torsionally weak bridges namely, 13 = 0.
Coefficients for any intermediate value of /3 may be obtained by using
following inter' polation function;
11(3+16 ____________________________________________________
ffrc7) (2.10)
A
'""--"= '
.hirhods of Bridge Deck Analysis 47
Space frame programs have been used in the analysis of box-girder bridge
decks. However, the simulation of boxes by space frames is not capable of
predicting local effects and the method has proved expensive in use.
Auxiliary
Nodal Line Strips
Nodal Lines
Figure 2.3 Nodal and Auxiliary Nodal Lines in Finite Strips
The continuum is divided into strips by fictitious lines called Nodal Lines
and these strips are assumed to be connected to one another along discrete
number of nodal lines which coincide with longitudinal boundaries of the
Methods of Bridge Deck Analysis 51
REFERENCES
L AGRAWAL, R., "Analysis and Design of Interconnected Skew Girder Bridges",
Ph.D. Thesis, Indian Institute of Technology, New Delhi, 1975.
2. BARES, R. and MAssozorET, C., -"Analysis of Beams and Grids and OrthOtropic
Plates by the Guyon-Massonnet-Bares Method", Crossby Lockwood, London, 1968.
3. BROWN, T.G. and GuAu, A., "Finite Strip Analysis of Skew Slabs", Proc. McGill-
EIC Conference in Finite Element Method in Civil Enginering, 1972.
4. CHEUNG, M.S., CHEUNG, Y.K. and CHAU, A., "Analysis of Slab and Girder
Bridges by the Finite Strip Method", Building Science, Vol. 5, 1970.-
5. CHEUNG,IK., "The Finite Strip Method in the Analysis of Elastic Plates with Two
Simply Supported Ends", Proc. ICE, 40, 1968.
6. CHEUNG, Y.K., "Analysis of Box-Girder Bridges by Finite Strip Method", Proc.
2nd Intl. Symposium on Concrete Bridge Design Chicago, ACI Publications, SP 26,
1969.
7. CHEUNG, Y.K., and CHEUNG, M.S.,,"Analysis of Curved Box-Girder Bridge by
Finite Strip Method," Publication, IABSE Vol. 31/1, 1971.
8. CHEUNG, Y.K., "Finite Strip Method in Structural Analysis", Pergamon Press,
Oxford, England, 1976.
9. CLOUGH, R.W., "The Finite Element in Plane Stress Analysis", Proc. 2nd ASCE
Conf. on Electronic Computation, Pittsburg, Pa., 1960.
10. COURBON, J., "Application de la Resistance des Materiaux au Calculdes Pants",
Dunod, Paris, 1950.
11. CUSENS, A.R. and PAMA, R.P., "Bridge Deck Analysis", John Wiley, London, 1975.
12. DESAI, C.S. and ABEL, J.F., "Introduction to Finite Element Method", Von
Nostrand Reinhold, New York, 1972. 1
Methods of Bridge Deck Analysis 53
13. GHAI I, A., "Designs of Simply Supported Skew Concrete Girder Bridges", Proc.
Intl. Symposium on Concrete Bridge Design, SP 26 Toronto, 1969.
14. -Gm A., "Grillage Analysis. Notes for Course on Bridge Deck Analysis", Civil
Engg. Dept., University of Dundee, 1972.
IS. GOLDBERG, J.E. and LEVE, ELL., "Theory of Prismatic Folded Plate Structures",
Publication IABSE, VoI. 17, 1957.
16. GUYON, Y., :`Calcul des Pants Larges a Poutres Multiples Solidarisees Par des
Entretoises", Annales des Ponts et Chausees, No. 24, 1946.
17. HEINS, C.P. and HAns, R.L., "Behaviour of Stiffened Curve Plate Model", J.
Struct. Div., ASCE, Vol. 95, ST11, 1969.
18. HENDRY, A.W. and JAEGER, L.G., "The Analysis of Grid Frameworks and
Related Structures", Prentice-Hall, Englewood Chatto & Windus, London, 1958.
19. JAEGER, L.G., BAiurr, B., and SURANA, C.S., "Application and Analysis of
Three-Girder Skew Bridges", Proc. Second Intl. Colloquium on Concrete in
Developing Countries, Bombay, 1988.
20. JAEGER, L.G., and &waif, B., "The Grillage Analogy in Bridge Analysis", Cana-
dian Journal of Civil Engineering, 9(2), 1982.
21. JAVOUR, J., Ikew Slabs and Gridwork Bridges", Bratislava Czechoslovakia:
Slovenske Uydavataletsvo Technickij Litratury, 1967.
22. JENSEN, V.P., "Analyses of Skew Slabs", Bulletin Series No. 332 University of
Illinois, Illinois, 1941.
23. JOHNSON, C.D. and LEE, T., "Long Non-prismatic Folded Plate Structures", J.
Struct. Div., ASCE, Vol. 94, 1968.
24. KREYSZIG, E., "Advanced Engineering Mathematics", John Wiley New York, 1962.
25. LIGEar-oar, E. and SAWKO, F., "Structural Frame Analysis by Electronic
Computer: Grid Frameworks Resolved by Generalised Slope Deflection,"
Engineering, 187, 1959.
26. Loo, Y.C. and CUSENS, A.R., "The Finite Strip Method in Bridge Engineering",
Viewpoint Publication, Cement and Concrete Association, London, 1978.
27. MAISEL, B.I., "Review of Literature Related to the Analysis and Design of Thin
Walled Beams", Technical Report 42:440, Cement and Concrete Association,
London, 1970.
28. MARTIN, H.C., "Introduction to Matrix Method of Structural Analysis," McGraw-
Hall Book Co., New York, 1966.
29. MARTIN, H.C. and CAREY, G.F., "Introduction to Finite Element Analysis",
McGraw-Hall, New York, 1973.
30. MASSONNET, C., "Methode de Calcul des Ponts a Poutres Multiples Tenant
Compte de Leur Resistance a la Torsion", Publication, IABSE, No. 10, 1950.
31. MEYER, C., "Analysis and Design of Curved Box-Girder Bridges", Report SESM
70-22, University of California at Berkeley. 1970.
32. MORICE, P.B., LITTLE, G. and RowE, R.E., "Design Curves for the Effects of
Con-Pui-micatIon uri I a, Lcmctit and Concrete Associaton,
London, 1956.
54 Grillage Analogy in Bridge Deck Analysis
33. NARUOKA, M. and OHMURA, I-I,. "On the Analysis of a Skew Girder Bridge by
the Theory of Orthotropic Parallelogram Plates", Proc. 1SBSE, No. 19, 1959.
34. NEILSEN, N.J., "Bestemmelese of Spaedinger i Piader Ved Anvendelse of
Differensligninger", Doctoral Dissertation, College of Engineering,
Copenhagen, 1918.
35. PATERSON, D.K,W., "Load Distribution in Skew Orthotropic Plates", Ph.D.
Thesis, University of DUndee, 1970.
36. PoWELL, G.I-1 and OGDEN, D.W., "Analysis of Orthotropic Steel Plate Bridge
Decks", Proc. Struct. Div., ASCE, ST5, V.95, 1969.
37. PaAsAn, J., "Modified Hormonics Method for Analysis and Design of Skew
Girder Bridges", Ph.D. Thesis, Indian Institute of Technology, New Delhi, 1982.
38. ROBINSON, K.E., "Behaviour of Simply Supported Skew Slabs Under Concentrated
Loads", Research Report No.8, Cement and Concrete Association, London, 1959.
39. ROWE, R.E., "Concrete Bridge Design", C.R. Books Ltd., London, 1962.
40. Sawx.o, F., "Computer Analysis of Grillages, Curved in Plan", Publication,
IABSE, 1967.
41. SCHLEICHER, C. and WEGENER, B., "Continous Skew Slabs: Tables for Statical
Analysis", Verlog fur Bauwesen, Berlin, 1968.
42. SCORDELIS, A.C., "Analysis of Simply Supported Box-Girder Bridges", Report
SESM 66-17, Dept. of Civil Engg., UniVersity of California at Berkeley, 1966.
43. SCORDELIS, A.C., "Analysis of Continous Box-Girder Bridges", Report SESM
6725, University of California at Berkeley, 1967.
44. SCORDELIS, A.C., "Analytical Solutions for Box-Girder Bridges", Proc. Intl. Conf.
on Developments in Bridge Design and Construction, Crosby Lockwood, London,
1971.
45. SURANA, C.S., "Interconnected Skew Bridge Girders", Ph.D.Thesis, University
of Edinburgh, U.K., 1968.
46. TURNER, M.J., CLOUGH, R.W., MARTIN, H.C. and.Torp, L.J., "Stiffness and
Deflection Analysis of Complex Structures", J, Aero Sci., 23, 1956.
47. WEST, R., "The Use of a Grillage Analogy for the Analysis of Slab and Pseudo-
Slab Bridge Decks", Research Report 21, Cement and Concrete Association,
London, 1973.
48. WEST, R., "Recommendations on the Use of Grillage Analysis for Slab and
Pseudo-Slab Bridge Decks", C&CA/CIRIA. Cement and Concrete
Association, London, 1973.
49. WESTERGAARD, H.M., "Formulas for the Design of Rectangular Floor Slabs and
Supporting Girders", Proc. ACI, 22, 1926.
50. ZIENKIEWICZ, 0.C., "The Finite Element Method in Engineering Science",
McGraw-Hill, London, 1971.
51. IRC:21-1987, "Standard Specifications and Code of Practice for Road Bridge:
Section III-Cement Concrete (Plain and Reinforced)", Indian Roads Congress,
New Delhi, 1991.
Chapter 3
3.1 INTRODUCTION
The bridge deck structure may be considered as an assembly of structural
members connected together at discrete nodes forming a grid. The
deformations and forces at nodes are inter-related by corresponding
stiffnesses. In order to satisfy 'the equilibrium and compatibility conditions
at each node, a large number- of simultaneous equations will result and the
manual solutions of these may be prohibitive. But using the matrix method
of structural analysis as a primary approach, it becomes possible to obtain
the computer-oriented solution.
It is assumed that readers have sufficient exposure to matrix methods of
structural analysis. But to maintain the continuity of discussion in succeeding
chapters, a brief presentation of stiffness method and its formulation for
bridge deck analysis by Grillage Analogy is included. A simple but general
computer program, illustrating different steps involved in matrix formulation
needed, is also presented. A sample grid is analysed as an example.
83
0
2
(a)
'31
L (b)
Figure 3.1 Global Axes and Degrees of Freedom
Zm
Figure 3.2 Local Axes and Degrees of Freedom
In the above figure, the single headed arrows denote translations and the
double headed arrows represent rotations and are taken as positive in the
directions shown following right hand screw rules.
3 6
E1, CJ
Xm k 5
Zm (a) L (b)
Yrn 1E1
_6E1 L2
pp-XM
-_ 2E1
Zm 12E1 L3
L3
(I) Unit 'Translation
Along Zm at End
Ym
G7
k 'Xm
1.0
Zm
2, L
E
1 r. __ Xm
4E1Z m ' 1 . 0
6E1 6E1
L2
(iii) Unit Rotation About Ym at End
Figure 3.4 (a) Determination of Member Stiffnesses Corresponding to Degrees
of Freedom at End j
earlier. All vectors are drawn in positive senses, but in cases where the
restraint actions are actually negative, a minus sign precedes the expression
for the stiffness coefficient. In general, it is possible for the member to
undergo any one or more of the six displacements shown in Fig. 3.4.
All the member stiffnesses shown in the figures are derived by determin -
ing the values of restraint actions required to hold the distorted member in
equilibrium. The stiffness matrix [K] for a grid member as above, therefore,
is of the order 6 x 6 and each column in the matrix represents the actions
caused by corresponding unit displacement. The member stiffness 'matrix
thus obtained is symmetric aria is snown in equation 3.3.
60 Grillage Analogy in Bridge Deck Analysis
Ym
6E1 6E1
_ Xm
J 0
k
12(1
12E1
L3
(iv) Unit Translation Along Zm at End k
vm
G GI
GAKt X m
1.0
Zm
6E1
L 2 6 E 1
L2
(vi) Unit Rotation About Ym at End k
Figure 3.4 (b) Determination of Member Stiffnesses Corresponding to Degrees of
Freedom at End k
12 El 6E! 12 El 6 El
0 0 2
L2 L3 L
GJ _GJ
0 0 0 0
6E1 L 4 EI 6 El L 2 El
[K.] L 2 L L (3.3)
12 El 0 6 El 11E1 0 6 El.
L3 L3 L2
L2 GJ
0 GJ 0 0
6 EI L 2 EI 6 El L 4E1
L2 L L2 L
Ym Y
Y Xm
m5
ems
l 0 0
= 0 cos y sin y (3.4)
_0 sin y cosy_
cosy siny
0 sin y cos y I
I
tyesm3
Wf
tn4
0 cosy siny ,,
W m]
10 0I
I 0 siny cosy_ O',6
or {c1} = [T] (3.5)
where {dm} and {d.'} are local and global member deformation vectors
respectively and [T] is a 6 x 6 transformation matrix derived from the rota-
tion matrix fill such that
(3.6)
12E1
R
;n:
L3.
6E1 GJ C2 4E1
34
".n2 s CS
Y?
mi3
6E1 GJ4E1"L LGJ52+4E1C2
12E1 6E7 6E1 12E1
L' L2 L2 c L3
M",5 6E1 GJ2 L 2E1 2 GJ 2E1 6E1 GI 4E1 2
Les L _ +-S L L2 L 2 L
S-
M 4 +1 76 L
6E1_ a...ILI 2 6E1 GJ 4E1 GJ 2 4E1
2
L2 L L L L L2 TCS---ECS TS
+r-C
(3.12)
where s = sin y, c = cos y
K1 0 K13 0 Kik 0 0 0 0 0
K22 0 K25 K.7) 0 0 0 0
K23
K33 K34 K35 K36 )0 0 0 0
Kh4 /(45 Kt.6 K47 KI.8 0 0
K 5 5 K 5 6 K 57 K 5$ 0 0
K66
1(67
K68
1(77 1(78K69
K79
Symmetric K88 K89 1(81
stored are further reduced by the fact that the lower triangle part of the
matrix can be obtained from the upper triangle part and vice-versa. So only
half band of the matrix need be stored and rest of the elements can be
deduced from this, whenever required. The half band-width of stiffness
matrix of a structure can be obtained as
Half Band-Width (Degrees of freedom per node) x (maximum difference
of numbers of connected nodes f 1)
Storage requirements can be further reduced by using Skyline Technique for
assembly. As shown in Fig. 3.6, Skyline is an envelope drawn in matrix nose
zero Ci erfierit S which do not have any non -zero elements
o
where [Km] etc. are partitioned submatrices and {pp} and {DR} are free and
[IC p
START
/ I N P U T : G R I D D
n o d e s , N o . o f e
C o o r d i n a t e s o f
M e m b e r p r o p e r t
c o n n e c t i v i t i e s ,
E T A I L S N o . o f
l e m e n t s ,
n o d e s ,
i e s a n d
H a l f b a n d - w i d
i
TypeINPUT: SUPPORT
of support, No. of DETAILS
supports,
Supported
in case of flexible supportconstant
nodes, Spring
Generate load
vector
4r
4.
Print
nodal deformations
4.
DO 300
Loop for solving for
all members
Compute
member end forces
Print
member end forces
3.8 EXAMPLE
111- ar.d nut of inaluies of the above program have been
illustrated with the help of a simple example of a skeletal grid shown in Fig. 3.8.
70 Grillage Analogy in Bridge Deck Analysis
Simply
supported
Notes:
(i) Load of 10,000 kg at node11
(ii)Dimensions in cm and
E,Gin kg/cm 2
Member Forces
Member Shear Force Bending Moment Torsion
No. End 1 End 2
1 0.145468E+04 0.205935E+05 - 311530E+06 0/02739E+05
2 0.114928E+04 0.318142E+06 -.547997E+06 0.231074E+05
3 0.470264E+03 0.563450E+06 -.657503E+06 0.193761E+05
4 - .466629E+03 0.680321E+06 -.586995E+06 0.497521E+04
5 -.121693E+04 0.599063E+06 -.355676E+06 - .519601E+04
6 -.158519E+04 0.355656E+06 -.386176E+05 -.732661E+04
7 0.196012E+04 - .721619E+05 -.319862E-F06 0361352E-1-01
8 0.263379E+04 0.297518E+06 - .824277E+06 0.107522E+05
9 0.406311E+04 0.80306R.E+k"..,:.; - .1615.i.;21-07 .71.7 i :13-;
Contd.
Stiffness Method Applied to Grillage Analysis 73
10 - .406311E+04 0.161569E+07 -.803068E+06 0.221355E+05
11 -.263379E+04 0.824277E+06 - .297518E+06 0.107522E45
12 - .196012E+04 0.319862E+06 0.721619E+05 0.361352E+01
13 0.158519E+04 0.386176E+05 -.355656E+06 -.732661E+04
14 0.121693E+04 0.355676E+06 - .599063E46 -.519601E-F04
15 0.466629E+03 0.586995E+06 - .680321E+06 0.497521E+04
16 -.470264E+03 0.657503E+06 -.563450E+06 0.193761E+05
17 - .114928E+04 0.547997E+06 - .318142E+06 0231074E+05
18 - .145468E+04 0.311530E+06 - .205935E+05 0.202739E+05
19 - .943386E+02 -.226007E+03 0.335797E+05 -.288976E+05
20 0.305410E+03 -.267163E+04 -.105307E+06 -.667882E+04
21 0.679011E+03 -.135654E-1-05 - .226501E+06 -.828866E+04
22 u.y.sotwit+Os -.263173E+05. - .304925E+06 - .595134E+04
23 0.750304E+03 -.157253E+05 -.249547E-1-06 - .134105E44
24 0.368260E+03 .149277E+04 - .128707E+06 0.152035E4.04
25 -.425352E+02 0.324875E+05 - .174490E+05 0.221261E+05
26 0.425352E+02 0:174490E+05 -.324875E+05 0.221261E+05
27 - .368260E+03 0.128707E+06 0.149277E+04 0.152035E+04
28 -.750304E+03 0.249547E+06 0.157253E+05 -.134105E+04
29 - .936893E+03 0.304925E+06 0.263173E+05 - .595134E+04
30 - .679011E+03 0.226501E+06 0.135654E+05 - .828866E+04
31 - .305410E+03 0.105307E+06 0.267163E+04 - .667882E+04
32 0.943386E+02 -.335797E+05 0.226007E+03 - .288976E+05
REFERENCES
1. BATHE, K.J., "Finite Element Procedures in Engineering Analysis", Prentice-
Hall, USA, 1982.
2. BEAUFAIT, F.W., RowAN, W.H., HOODLEY, P.C. and HACKETT, R.M.,
"Computer Methods of Structural Analysis", Prentice-Hall, USA, 1970.
3. CHAPRA, S.C. and CANALE, R.P., "Numerical Methods for Engineers" McGraw-
Hill, New York, 1989.
4. CUSENS, A.R. and PAMA, R.P., "Bridge Deck Analysis", John Wiley, 1975.
5. GERE, J.M. and WEAVER, W., "Analysis of Framed Structures", Van Nostrand, 1965.
6. JENKINS, W.M., "Matrix and Digital Computer Methods in Structural
Analysis", McGraw-Hill, New York, 1969.
7. RUBINSTEIN, M.F., "Matrix Computer Analysis of Structures", Prentice-Hall,
USA, 1966.
8. SURANA, C.S., "Grillage Analogy for Analysis of Super-Structures of Bridges",
Proceedings of the Seminar on Modern Trends in Analysis and Design, B.H.U.,
Varanasi, February 1984.
9. SURANA, C.S. and HUMAR, J.L., "Beam and Slab Bridges with Small Skews",
Canadian Journal of Civil Engineering, Vol. 11. No. I, 1984.
Chapter 4
Transformation of Bridge
Deck into Equivalent Grillage
4J INTRODUCTION
In recent years, the Grillage Analogy Method, which is a computer-oriented
technique, is increasingly being used in the analysis and design of bridges.
The method is also suitable in cases where bridge exhibits complicating
features such as heavy skew, edge stiffening and. isolated supports. The use
of computer facilitates the investigation of several load cases in shortest
possible time. The method is versatile in nature and the contribution of kerb
beams and the effect of differential sinking of girder ends over yielding
bearings (such as neoprene bearing) can also be taken into account and large
variety of bridge decks can be analysed with sufficient practical accuracy.
Furthermore, the grillage representation is conducive to give the designer a
`feel' for the structural behaviour of the bridge and the manner in which the
bridge loadings are distributed and eventually taken to the supports.
The method consists of 'converting' the bridge deck structure into a
network of rigidly connected beams at discrete nodes i.e. idealizing the
bridge by an equivalent grillage. The deformations at the two ends of a
beam element are related to the bending and torsional moments through
their bending and torsional stiffnesses.
These moments are written in terms of the end-deformations employing
slope-deflection and torsional rotation-moment equations. The shear force
in the beam is also related to the bending moment at the two ends of the
beam and can again be written in terms of the end-deformations of the
beam. The shear and moment in all the beam elements meeting at a node
and fixed end reactions, if any, at the node, are summed-up and three basic
statical equilibrium equations at each node namely EF = 0, EMx = 0 and
EMy = 0 are satisfied.
Transformation of Bridge Deck into Equivalent Grillage 75 ,
In general, a grid having 'n' nodes will have '3n' nodal deformations
and `3n' equilibrium equations relating to these. Back substitution in the
slope-deflection and torsional rotation-moment equations will give the
bending and torsional moments at the two ends of each beam element.
Shear forces are computed from bending moments and external loads.
When a bridge deck is analysed by the method of Grillage Analogy, there
are essentially five steps to be followed for obtaining design responses:
(i) Idealization Of physical deck into equivalent grillage
(ii) Evaluation of equivalent elastic inertias of members of grillage
(iii) Application and transfer of loads to various nodes of grillage
(iv) Determination of force responses and design envelopes and
(v) Interpretation of results
'The first two steps of grillage analogy are discussed in this Chapter. The
other remaining steps are dealt with in the subsequent chapters.
yL2) t3
x
Ly 24 (1 y2)
2 3
r
ly
2
-j 24 (1 7 )
76 Grillage Analogy in Bridge Deck Analysis
.(E Ly (1-3y)) t3
G 24(1y2 )
.(E L, (1.-3y) t3
2
G 24(1y )
= 7(121 +14)13 _________________________________________________ t3-
/ (4.1)
LyJ24(1_72)
where y= Poisson's ratio and E, G Elastic moduli.
Ly
Y
Torsioniess Diagonal
Members
3
t Lx 3
/
=
Ix = 24 Y
24
WAMOM
MOO_Off
MMAM rAM
(a) Nan of Idealized
Grillage Accounting
for Poisson's Ratio
b) Skew bridge
Skewness has considerable effect on the behaviour of the deck and criti-
cal design stresses. Skew decks are associated with special characteristics
like hogging moment and large reactions near the obtuse corner and small
reactions and possible uplift at the acute corner. Moreover there is consid-
erable torsion of deck. Thus special attention is required while laying out
the grid lines for a skew bridge.
Decks with skew angles less than 15 can usually be handled as right
decks. Same guidelines as discussed above, for right bridges, are applicable
for such decks with or without footpaths and do not need any further elabo-
ration for setting out the grid lines. However, bridges having skew angles
more than 15 pose problems in regards to the positioning and orientation of
longitudinal and transverse grid lines and are discussed here.
A skew deck can be analysed with grillage having either a parallelogram
mesh as shown in Fig. 4.5(a) or orthogonal meshes as in Figs. 4.5(b) and 4
.5(c).
Transformation of Bridge Deck into Equivalent Grillage .81
-
r -
r-r-T-7 - - 1 XmL
I 1 I I I I I
I 1 1 I I I
I L --1 A_ i I 1
r T T ---
I I I I1 1 11
I _1_11 ! I
I 1 I II I
_1. _____ -- Grid-Lines
I I 1- I I
i I 1 I I 1
L - i-- - -I- --1- - 4-- -4- -
_I
I I ' 1 1 I
I- T - TI - i-I- 1 - - - 4 - - 1 I
I 11_1 I I
I 1
I Ii , ! I
I 1 1 I
,L__4_,_1-1_I_J__ -I X= 0
I I I I I I- I
______
I (a) 11Plan
i of Solid
11 1 1 1
Slab
1
1
LL __i __..,_ _ _1_ L _ _I
-.4 b
1D
I_________________________________________________________________________________1-
(b) Section on A-A
Figure 4.3 Solid Slab Grillage Lay-out without Footpath
While the parallelogram mesh (Fig. 4.5a) is convenient for low skew
angles, it is not appropriate for angles of skew greater than 15 because it has
-.no members close to the direction of dominating structural action. For
bridges with larger skew, say. greater than 15, a parallelogram mesh as Fig.
4.5(a) will result in an over-estimated maximum deflections and moments,
the amount increasing with angle of skew. The quantity of reinforcement in
such grids is likely to be excessive and uneconomical. An orthogonal grid
lay-out as shown in Figs. 4.5(b) or 4.5(c) will be more realistic.
In skew bridges, the direction of principal bending moment across the
Farancl to skcw span at edge to near normal to
.$
411
rI - - r- - r - - -T - - - 7 - I-- -1 X1.
1 1 I I tI I t 1
I 1 I
1I 11 tI I
I -- i --- 4. - - - 4 - - - -I- - - --I
I
I I 1
I
I 1
I I I
1 I
i
r- --II- -- -1-- --1- -- --11----F - - H A
I I .1
I I I 1 I . 1 I I I
1
x
1 1 1
1 I I I
r-- - i I l i t i- t 1 i
1 ii I 1 I
t 1 1 I
11 I. I ___I
t 1 1 i1 t
i 1 I t I I
I l 1
1
t1
r-- - -I1--.-r - -- -t- --- ir - -- -I- - 1
I I I I I I
I 1 Ii I I
1 t
i
t__ __ i__ __i_ ___L _ _ _ I _ _ _ _ L _ _ _I
Footpath
u ' ( 1 4 . 1 4 1
i s I , T
(b) Section on A. A
abutment in the centre. This is illustrated in Fig_ 4.6. While adopting the
orthogonal arrangement of grillage, two different cases arise, depending
upon the skew angle and span-width ratio as shown in Fig. 4.7.
Case (a) refers to the situation where skew region is small and right
Nyc m re-,,r,-.s.rt. this rasp when length AE < EG. The
recommended limiting condition is L sin 2, .S B/2 (Fig. 4.7a).
E. 1
Transformation of Bridge Deck into Equivalent Grillage 83
Diaphragm Beam
PP-
INN AI
I
1IP"
XIII1111LS
I "
SiIN111111
(a) Skew or Paretlelogram (b) Mesh Orthogonal to
Mesh Span
II II III II II 6..
1IIIIIIIIIIII Edge Beam
IIII II I II 111 k
(c) Mesh Orthogonal to Support ,
Figure 4.5 Grillages for Skew Decks
Case (b) represents the situation where skew region is large and right
region is small i.e. AE > EG. The recommended limiting condition here may
be L sin A > B/2 (Fig. 4.7b).
The above two cases will decide the orientation of longitudinals with
respect to supports. In case (a), the main spanning will follow the direction
normal to support and the longitudinal grid lines will be taken normal to
support lines and transverse grid lines will be perpendicular to them as shown
in Fig. 4.8 (a).
As against this, in case (b). the main spanning will fnllnw ctepw rur,:ction
and hence longitudinal grid lines will be taken parallel to bridge axis and
84 Grillage Analogy in Bridge Deck Analysis
Right Region
Skew span
B
(b) La rge Skew Region
( L sin h > 6/2 )
Figure 4.7 Cases of Small and Large Skew Regions
Case (a)
The spacing of longitudinal grid lines may be different in skew region and
right region. The spacing of these longitudinal grid lines is to be neither
greater than the spacing of transverse grid lines nor greater than three times
the thickness of slab. In the right region, .one 11/4-)nezlit4difia1 caLli jakcn
along ED and GF and one or more between ED and GF. In skew region,
Transformation of Bridge Deck into Equivalent Grillage 85
F C
X-L
-Longitudinds
_ Transversals
14 . / * 4 .
/ Transversals
/
/ / Longitudinals
A /
/ 1'4 / /
A
(b) Longitudinal Grid Lines
Parallel to Bridge Axis
Figure 4.8 Grid Lines in Skew Slab Decks
the longitudinal grid lines may originate from the nodes where the transverse
lines meet the longitudinal grid lines. However, if the number of longitudinals
in skew region become too many, alternate or even lesser longitudinal lines
may suffice. As far as possible, the longitudinal grid lines should be
equidistant in skew region and right region each.
Seven transverse grid lines for skew decks of span upto 10 m and nine
grid lines for spans above 10 m may be used. Out of these transverse grid
lines, one line each be located at ends of the span joining the centre line of
bearings on each abutment and one at the centre of span. The remaining
lines can be set in between these. As far as possible, the grid lines be kept
equidistant.
86 Grillage Analogy in Bridge Deck Analysis
Case (b)
The lay-out of lbrigitudinal grid lines will be similar to as in the case of
right bridge discussed earlier. The actual numbers of transverse grid lines
will be decided so that the spacing of these lines is not greater than 1.5 to 2
times the spacing of longitudinal grid lines. Referring to Fig. 4.8(b), one
grid line each is set along GQ and DR and one at midway between GQ and
DR. Others may be set between these if necessary. In triangular portions
AGQ and DCR, transverse grid lines may be taken from the end of each
longitudinal. However, if the number of transversals in this region
becomes too large then alternate or even lesser transverse lines may also
suffice. Again, in this case also, the transverse grid lines be placed
equidistant, as far as possible.
In skew slab deck with foot path, one extra longitudinal grid line is
taken for the footpath at its centre at each end of the deck width as in the
case of right deck.
Sometimes a thick slab deck with thin cantilever and connecting slabs, as
shown in Mg./4.9 is also encountered in practice. The thin cantilever slab at
the edge, supports the kerb, parapet etc. and the middle slab is meant for
road divider or verge. The longitudinal grid lines can be placed as in Case
(a) or as in Case (b). Location of the longitudinal grid lines as shown, assign
the sectional properties of the transverse grid members better between 3 and
4 in Fig. 4.9(a) or between 4 and 6 in Fig. 4.9(b). It may be also noted that in
Fig. 4.9(b) longitudinal grid lines through 1, 5 and 9 will not have zero end
deflections.
Thin Cantilever
(b)
Figure 4.9 Thick Slab Deck with Thin Cantilever and Connecting Slab
AM .
Partipet
\ m
1 1
I I j I I I - 0
3. Pseudo-Slib Bridge
As discussed in Chapter 1, contiguous beams spaced closely with in-situ
concrete comes under the category of pseudo-slabs. Usually standard pre-
cast beam sections are used in the construction. In such a situation, the
number of closely spaced beams are large. A grillage arrangement with
longitudinal grid lines coincident with every physical beam will lead to
many grid linpc which ,Pc....*.1 arid unmanageable. Therefore, it
is proper to represent more than one physical beam by a longitudinal grid line.
88 Grillage Analogy in Bridge Deck Analysis
8 @ 1= 8 M
V
(a) Cross Section O
I
3 I f
- . (b) Longitudi nal Gri d D
S
Lines'
.1F 5
I
. 1 I
.
I-- - I-- -- - , --4- 1 ----I
I 1
1 1 i
I I I I
I i I
I I I II
I I I I
I, 1 1
I I I
I
I 1
l- - +---- i-- - - + - --- -i Iii;
I I i 1 c.4
1 1 1 1e
1 1 1 1-
1 1 I
----- 1 -.- --1---- - 1 ----1
1 1 1 1
I I I 1
I I I I
11-- -I-- -- ---
1 I -1 i
I i I- 1 - --
1 I
1 1
I 1 I I
I I 1 I I
1--- -i- - --
i E. i
k--
r to n
--Jk-- - -1
Figure 4.11 Grillage Geometry of Inverted T-Beam with in-situ Concrete Slab-Pseudo Slab
Bridge
Transformation of Bridge Deck into Equivalent Grillage 89
Voids
i
1-
EiftiOrD14ffiritrIfffe
(a) Box-Beam Deck
I I -0 1
(c) Cross-Section
(a) Plan
iiiiiiiii
1101111111
11
( d) G ril l a ge Lay - out
The transverse grid lines are also placed at abutments joining the centres of
bearings. A minimum of seven transverse grid lines are recommended,
including end grid lines. It is advisable to align the transverse grid lines
normal to the longitudinal lines wherever cross-girders do not exit. It should
also be noted that the transverse grid lines are extended upto the extreme
longitudinal grid lines.
In skew bridges with small skew angle say less than 15 and with no
intermediate diaphragms, the transverse grid lines are kept parallel to the
support lines as shown in Fig. 4.15(a). Additional transverse grid lines are
provided in between these support lines in such a way that their spacing does
not exceed twice the spacing of longitudinal lines, as in the case of right
bridges, discussed above.
IS
In skew bridges with higher skew angle, the transverse grid lines are
set along abutments PQ and WV and also along PR, ST and UV initially
as shown in Fig. 4.15(b). Then extra grid lines are inserted in between PR
and UV. When the span-width ratio and the skew angle are such that the
skew region QR is larger than the spacing of transverse grid lines in right
region VR, additional transverse members parallel to interior transverse
grid lines starting from the interior nodes are to be placed.
FQ
(
Z.
y,y ,R
Lb ;is./
>1 7T ,
(a) Deck with Small Skew (b) Deck with Large Skew
A>15
Figure 4.15 Grillage Arrangement in Skew T-Beam Bridge without Cross Girders
I-
.0
- -- ___________ elipt
Longitudinal grid lines are usually placed coincident with webs of the
actual structure (Fig. 4.16). If the deck has sloping end webs (Fig. 4.17), the
grillage simulation is not so precise and engineering judgement must be used
to position longitudinal members. However, one grid line may be taken at the
junction of the inclined web with slab as shown. Normally, additional
longitudinal grid lines are located along the edges of the side cantilevers with
nominal stiffnesses for the convenience of analysis. Additional longitudinal
grid lines are to be adopted for bridges with footpaths at their centres as done
in the case of slab and T-beam bridges, discussed earlier.
The transverse medium consisting of top and bottom slabs only (with no
diaphragms), is represented by equally spaced transverse grid lines along
94 Grillage Analogy in Bridge Deck Analysis
the span. The spacing and number of grid lines are similar to as adopted for
slab-on-girders bridge. If the deck is having diaphragms, the transverse grid
lines are placed along each diaphragm including at supports. Additional grid
lines representing the top and bottom slabs are placed in between the dia-
phragms, if needed, to meet the minimum:requirements of transverse grid
lines. A closer spacing of transverse grid lines will result in more continuous
structural behaviour and will provide greater details of forces etc.
For skew box-girder bridges, the procedure to.be followed regarding the
setting of longitudinal and transverse grid lines will be the same as in the case
of slab-on-girders bridge, discussed earlier.
Spaced box-girders, also referred to as Spine box-girders, connected by
top slab only, form a distinct class of decks requiring a special approach.
Figure 4.18 represents a four-cell spaced box-girder bridge.
In deciding the grillage layout for this type of deck, the main problem lies
in defining the stiffness in the transverse direction. In the transverse direction,
the deck is alternately very stiff over the box-beams and very flexible between
the box-beams. It is not possible to replace the transverse media by grillage
beams of uniform stiffness, as is usually possible with other types of decks. In
this type of deck, each physical beam is replaced by grid lines positioned at
the webs as shown in Fig. 4.18(b). This takes into account the effect of abrupt
change of transverse section more correctly. The other de- tails of layout of
transverse grid lines are similar to those of slab-on-girders bridge. The plan of
the arrangement of grid lines is shown in Fig. 4.18(c).
Sometimes, the width of the cell is small in comparison with their spac-
ing and walls are relatively thick preventing distortion of the cell (Fig. 4.19).
For such sections of bridges, the longitudinal grid lines are better placed
coincident with centre lines of boxes with additional nominal grid lines
running in between the boxes as shown in Fig. 4.19 (b). In the transverse
direction, the grid lines will have uniform stiffness corresponding to the
deck slab and other details of layout will be similar to those of slab-on-
girders bridge, discussed earlier.
FF
Transformation of Bridge Deck into Equivalent Grillage 95
1 --4----- t i - -- t - i - -- 4--
i
I I 1 I I 1 ,I 1
_ 4 . . . .
t-
I 1 T1
1_1_ I , I I 1 1
I-
1
1
T -7-
--1
r 1--7-
I -- - - I j I I I
II t - t--- - - F 4
--
I I I 1,
4
1 I I
1-- 1- - --- I - I - -- i
-1-
- - I- -1 A--k--
A k A i 7 1c --- i--- 1--
bi
I
i i
er,
(bi- b2) -
\,.../
1-..-b4,-.-1 1--b.-,-
(a) Cross- Section of Deck with Narrow .-1Box-Beams
1 2
to 3
( b4/2
(d) Nominal Stab
Grid Line
3 b3 d3
= (4.5)
10 (b2 + d2 )
In the case of a thin rectangle where b >
5d, the J value is more accurately given by
bd3
J= (4.6)
3
If the cross-section has reentrant
corners, J is very much less than that given by equation 4.4 above [13]. In
such cases, the value of J is obtained by notionally sub-dividing the section
into rectangular shapes without having reentrant corners and summing the
values of J of these elements. The value of J of a sub-divided portion with
notional cuts on two opposite faces' is to be computed as if sub-division is a
part of wide thin strip for which J = bd3I 3. Figure 4.20 shows a T-section
with reentrant corners and its sub-division. Thus, if J values of the portions
1, 2, 3 and 4 are designated as J , J , J3 and J4 respectively, then,
1 2
(1 3)
J 2 =- 1 (b 2
Jr = - " di
' i3 ) 3
3 b33 3 b: d34
J
4 2 2
J 3 10 (1 4 10 (b 4 d4)
d41:
1I d i
1. ______
b1
ro2
(4.7)
It may be noted that the value of J of the portion of deck slab forming the
flange is' to be halVed to take into account its continuity in the other direction
[12]. Widths b3 and b4 of segments 3 and 4 are so adjusted that areas b3xd3
and b4xd4 are same as original areas of the respective segments.
While notionally sub-dividing the section, it is worth remembering
Prandtl's membrane analogy. It is shown that the torsional stiffness of a
cross-sectional shape is proportional to the volume under an inflated bubble
stretched across a hole of the same shape. Thus, care is to be taken so as to
obtain the largest numerical value of J of the section as a whole. This is
achieved by choosing the elements so that they maximize the volume under
their bubbles [11, 12]. Sub-dividing the cross-section arbitrarily into rect-
angles and not trying to maximize the. volume under the bubble will result in
lower value of J.
As an alternative to Equations 4.5 and 4.6 above, expression given by
Timoshenko and Goodier [12] can be used for rectangular sections viz. J =
yrbd3 where ty depends on the ratio of sides b and d (b > d) of the rectangle.
The values of tir for different ratios of b and d are given in Table 4.1_
TABLE 4.1
bld 1 1.2 1.5 2 2.5 3 4 5 10 00
.0 .141 0 .166 0 .19 6 0 .229 0 .249 0 .26 3 0 .281 0 .291 0 .312 0333
i 3d =
12 t
If the deck has thin cantilever or intermediate slab strips as in Fig. 4.9, the
longitudinal grid lines can be placed as in Fig. 4.9(a) or as in Fig. 4.9(b). In
Fig. 4.9(a) the flexural inertias of all members are calculated about the deck
neutral axis. However if the grid lines are placed as in Fig. 4.9(b), the thin
slabs above members 1,5 and 9 act primarily as flanges to members 2, 4, 6
and 8 respectively. Consequently the inertias of members 1, 5 and 9 are
calculated about centroid of the thin slab, while for members 2, 4, 6 and 8
inertias are calculated with the flanges as in Fig. 4.9(a) but with small
inertias of members 1, 5 and 9 deducted. In the transverse direction, the thin
slab would bend about its own centroid only so that the thin slab depth is
used for members 1-2, 4-5, 5-6 and 8-9 while the thick slab depth is used for
members 2-3, 34, 6-7 and 7-8.
For a voided slab deck such as in Fig. 4.21, the inertias of longitudinal
grid lines are calculated for shaded section about neutral axis. Transversely,
fa
Neutral
xis
Comparing the values of i and j per unit width for solid slab, it will be
seen that
j = 2i (4.9)
There is no simple and rigorous rule for calculating j for voided slabs
and the above relationship between i and j can be used with advantage to
compute j for voided slabs with sufficient accuracy.
Further, the torques in two orthogonal directions of an orthotropio slab are
of equal magnitudes. Consequently, the transverse and. longitudinal grillage
members should have identical torsional inertia per unit width of the deck. It
is suggested that, for such slabs, the torsional inertia of the transverse and
longitudinal grillage members, per unit width, be computed as,
j = 2 ix i (4.10)
where i and i are the longitudinal and transverse flexural inertias per unit
width of slab respectively.
Alcn if hac hppn nhcerveri that the. simulation of grillage and physical deck is
.
of. edge member [6]. However, for simplicity, full width of edge member may be used.
In pseudo-slab construction, usually precast prestressed concrete beams are closely
placed with in-situ reinforced concrete forming the deck. In such cases, the in-situ
concrete slab has lower strength and stiffness than. the prestressed concrete beams so that
it has a modular ratio m = 0.8, compared
to prestressed concrete. Different transformed
cracked section inertias are used in the
two directions because the transverse
reinforcement is usually lighter.
It is to be noted that in a skew gril-
lage, the slab width which is to be
taken for computing the inertias, is the
one normal to the grid line. Referring
to Fig. 4.22, the widths of the slab to be
considered with the grid lines A and B
B of the skew system will be U cos a. and
Figure 4.22 Slab widths for V cos A. respectively. )
inertia in a skew grillage layout
43.4 Flexural
and Torsional Inertias of
Grillage Members: Slab-on-Girders Deck
Slab-on-girders bridge decks consist of a number of beams spanning
longitudinally between abutments with a thin slab spanning transversely
across the top. T-beam bridges are the common examples 'tinder this cat-
egory. The be___ ms may be cast monolithically with the slab or the precast
beams with in-situ slab may be used. The decks may be with or without
intermediate and/or end diaphragms.
The thin slabs can be thought of as flanges of I or T-beams. When such I or
T-beams bend, the flanges are subjected to flexural stresses. An element of
the flange away from the rib or stem of the beam has less stress than the one
directly over the rib due to shearing deformations of the flange. Shear
deformation relieves some amount of compressive stress in more distant
elements. This phenomenon is known as shear lag. The variation of corn:
pressive stress across the width of flange is accounted for by considering the
effective width of flange. The effective width may be smaller than the actual
width and is considered to be uniformly stressed. The effective width of flange
is determined based on the concept of constant compressive stress over the
entire effective width such that the total compressive force carried by the
effective flange is the same as that by the actual flange width with variable
compressive stress (Fig. 4.23). Tnis effective width has been found to depend,
besides other things, on the span and the relative thickness of slab.
102 Grillage Analogy in Bridge Deck Analysis
HEffective widthH
Variation of-' -1_..-- Assumed Uniform
' Ftexural. ' >--'r----'-'------- Stress over the
----
S t r e s s
Effective width
1. In case of T-beams
(i) One-fourth the effective span of the beam
(ii) The distance between the centres of the ribs of the beams
(iii) The breadth of the rib plus twelve times the thickness of the slab
2. In case of L-beams
(i). One-tenth the effective span of the beam
(ii) The breadth of the rib plus one-half the clear distance between the
ribs
(iii) The breadth of the rib plus six times the thickness of slab
The flexural inertia of each grillage member is calculated about its cen-
troid. Often the centroids of interior and edge member sections are located at
different levels. The effect of this is ignored as the error involved is
insignificant.
Once the effective width of slab acting with the beam is decided, the
deck is conceptually divided into number of T or L-beams as the case may
be. Some portion of the slab may be left over between the flanges of
adjacent beams in either directions. In the longitudinal direction, it is suf-
ficient to consider the effective flange width of T, L or composite sections,
in order to account for the effects of shear l ag and irrn^r- th 1-ft
Transformation of Bridge Deck into Equivalent Grillage 103
over slab. However, in the transverse direction, the left over slab should
be considered by introducing additional grid lines at the centre of each
left over slab portion.
Sometimes, for the purpose of improving the simulation, it is desirable
to place additional grillage members of nominal stiffnesses between grid
lines representing the beams. The sectional properties of these additional
members are calculated in a similar manner as outlined in Section 4.3.3
for the deck of Fig. 4.9.
The sectional properties of grid lines representing the slab only, are
calculated in the usual way i.e. I = bd3/12 and J = bd316.
If the construction materials have different properties in the longitudinal
and transverse directions, care must be taken to apply correction for this.
For example, in a reinforced concrete slab on precast prestressed concrete
beams or on steel beams (Fig. 4.24), the inertia of the beam element (1 or
J) is multiplied by the ratio of modulus of elasticities of beam Eb and slab
E: materials to convert it into the inertia of slab material. For example the
total torsional' inertia in terms of the slab material will be given by
Effective Flange Width
1 --*-1 J = Jl + (J2 + J3 + J4) Eb
_________ b E s
-d1
d2 J1, J2, J3 and J4 are calculated as below
=12-(13-121 4)
3b 3 d 3
j 22
2
10 (i4 -
13
- Figure 4.24 Torsional Inertia in
Precast Beam with Cast-in-Situ Slab
J3=3b3d3
3b3 d3
J 44
4
10 (hi + 41)
Following the notations given in Fig. 4.25, the moment of inertia of transverse grillage
member per unit width, it, is calculated about the common centroid and is given by
= 1 + t2 d2(4.11) since ti d1 = t2 d2 and (di +
d2) = d, we get
2
3 b3d3 b.?d?
J= (4.13)
10 (b2 + d2) (0 +4)
J =4 A2 (4.14)
f ds
t
where 'A' is the area bounded by the centre line of the closed cross-section
and 's' and 't' refer to the length and thickness of different units respec-
tively.
Referring to the Fig. 4.27, the equation 4.14 can be written as
4A2 4A2
J= (4.15)
(E (Si S2 2S 3
t t t
l 2 3
______________________________________ 4
-7
F - - - 2
Figure 4.27 Thin-Walled Box-Section
106 Grillage Analogy in Bridge Deck Analysis
I-
t2
(a) Equivalent One -Cell Cellular Deck
J 4 b2d2 (4.16)
(b b 2d)
t2 t3
The term 2d/t3 is small compared to other two terms in the denominator of
equation 4.16 and can be ignored. The equation 4.16 can then be re-written as
J 4 bd2ti (4.17)
t 2)
L
Transformation of Bridge Deck into Equivalent Grillage 107
(b)
Comparing equations 4.12 and 4.18, it is seen that] =. 2i. This is the same
relationship as was obtained for solid slab bridges earlier and can also be conveniently
used for computing the values of J of the transverse members of multi-cell boxes of
closed section without diaphragms.
For obtaining J values of longitudinal members, cell-widths in cross-sections are
considered and for J values of transverse members, spacings of transverse members are
to be taken.
If j is the torsional inertia per unit_width, referring to Fig. 4.29, the values of J for
longitudinal members 1, 2 and 3, will be given by
T
a
(b) Plan
Figure 4.29 Evaluation of J for Longitudinal and Transverse Members of Multi-
Cell Box-Girder Bridge
= 1/2j h and J, = j h
and for transverse members 4 and 5. .1 vaiiips will
J4= = ja
= O M
4 1 6 1 1 . 1 1 1 1 . -
'1 2 3 4
1
the value of the flanged box section. For the transverse grillage beam be-
tween members 1-2 (Fig. 4.30b), the value of I will be based on top and
Transformation of Bridge Deck into Equivalent Grillage 109
bottom slabs computed about their combined centroid and the value of I for
the beam between 2-3 will be based on top slab about its own centre of
gravity. The values of J for beams between 1-2 and 2-3 can be determined
from the usual relation i.e. j = 2i.
Thus, it is seen that the actual deck of the bridge can be idealised into
suitable grillage with appropriate equivalent inertias allotted 'to its each
member. Other steps in grillage analogy i.e. the application of loads and its
transfer to various nodes of grillage, analysis, stress response,
interpretation of results etc. are dealt with in the next chapter.
REFERENCES
1. BAKHT, B. and JAEGER, LG., "Simplified Analysis for Slab-on-Girder
Bridges", Bridge and Structural Engineer 12(4), 1982, INDIA.
2. BAKHT, B., JAEGER, L.G. and CHEUNG, M.S., "Cellular and Voided Slab
Bridges", Journal of Structural Division, ASCE (ST9), 1981.
3. BAxnr, B. and JAEGER, L.G., "Bridge Analysis Simplified", McGraw Hill,
New York, 1985.
4. BAKHT, B.s JAEGER, L.G., CHEUNG, M.S. and MuFrr, A.A., 'Me State-of-the
Art in Analysis of Cellular and Voided Slab Bridges", Canadian Journal of
Civil Engineering 8(3), 1981.
5. CUSENS, A.R. and PAMA, R.P., "Bridge Deck Analysis", Wiley, London, 1975.
6. HAMBLY, E.C., "Bridge Deck Behaviour", Chapman and HaIl Ltd., London, 1976.
7. JAEGER, L.G. and BAKHT, B., "Bridge Analysis by Micro-computer", McGraw-
Hill, New York, 1990.
8. JAEGER, L.G. and BAxtrr, B., 'The Grillage Analogy in Bridge Analysis",
Canadian Journal of Civil Engineering 9(2), 1982.
9. LtarrFocrr, E. and SAWKO, F., "Structural Frame Analysis by Electronic Com-
puter: Grid Frameworks Resolved by Generalised Slope Deflection",
Engineering, 187, 1959.
Id. MAISEL, B.I., "Review of Literature Related to the Analysis and Design of Thin
Walled Beams", -Technical Report 42.440, Cement and Concrete Association,
London, July 1970.
1 I. OMEN, J.T., "Mechanics of Elastic Structures", McGraw Hill, New York, 1967.
12. TIMOSHENKO, S. and GOODIER, J.N., "Theory of Elasticity", McGraw-Hill, New
York, 1951.
13. WEST, R., "Recommendations and the Use of Grillage Analysis for Slab and
Pseudo-Slab Bridge Decks", Cement and Concrete Association, CIRIA,
London; 1973.
14. WEST, R., "The Use of Grillage Analogy for the Analysis of Slab and Pseudo-
Slab Bridge Decks", Research Report 21, Cement and Concrete Association,
London, 1973.
110 Grillage Analogy in Bridge De& Analysis
15. YET-IRANI, Al. and HusATN, M.H., "A Gridwork Frame Analogy for Plates in
Flexure", Journal of Engineering Mechanics Division, ASCE, June 1965.
16. IRC: 21-1987, "Standard Specifications & Code of Practice for Road Bridges,
Cement Concrete (Plain and Reinforced), Section III", The Indian Roads Congress,
1988.
-
Chapter 5
Application
Application of Loads,
Analysis, Force Responses
and Their Interpretations
5.1' INTRODUCTION
The bridge deck has been transformed into an equivalent grillage consisting
of longitudinal and transverse grid members such that the idealized grillage is
very close to the physical deck. Each member of the grillage is allotted
flexural and torsional inertias which are equivalent to the corresponding
physical properties of the bridge deck. The longitudinal and transverse grid
lines form a mesh having number of nodes. The bridge is mainly subjected to
vertical loads comprising dead, live and impact loads. Grillage analysis
requires that these applied loads be transformed into equivalent loads at
nodes. This is done by finding equivalent nodal loads in terms of either
vertical load only or alternately vertical load, bending moment and torsional
moment.
This Chapter discusses different types of loads, identification of panels in
which the wheel loads of a vehicular loading system fall and transfer of loads
to nodes of grillage in the form of equivalent nodal loads. The analysis of
grillage is then carried out and response envelopes and the interpretation of
results are discussed. Both right and skew decks have been dealt with. Local
effects which are to be accounted for in arriving at the design responses, are
also discussed.
I I ' I I
I I I
I I
1I
I 1
Ii
Longitudinal Grid Line
2c
-L- Tributary Area
2c
Tributary Area
b'p+p 2a pqr,p, 2a pi
of longitudinal grid lines. The equivalent vertical load along each grid line is
computed in the form of vertical u.d.l. based on its tributary area. For grid .
line 'A' the loading is non-central and hence the equivalent load will consists
of a vertical -and torsional u.d.l., whereas for grid line `B' the loading is
central and hence the equivalent load will be a vertical load only.
The self weight of cross-beams and diaphragms needs further consider-
ations. These beams, located at specific intervals, are actually small discrete
loads on the longitudinal girders. However, for simplicity of computation,
the total weight of all the cross-beams per span should be calculated and
equally divided in the form of distributed loads to various longitudinal
members of the grillage. The dead weight of railings, kerbs, footpaths etc.
is lumped on the edge longitudinal grid lines.
b) Live Load
The main live loading on highway bridges is of the vehicles moving on
it. Indian Roads Congress (IRC) recommends different types of standard
114 Grillage Analogy in Bridge Deck Analysis
c) Impact Load
Another major loading on the bridge superstructure is due to the vibrations
. Application of Loads, Analysis, Force Responses and Interpretations 115
caused when the vehicle is moving over the bridge. This is considered
through impact loading. IRC gives impact load as a perCentage of live
load. As per IRC Code, impact load varies with type of live loading, span
length of the bridge and whether it is a steel or a concrete bridge. The
impact load can be calculated using formulae or could be directly read
from ready-to-use graphical plot (Fig. 1.15). The impact load, so
evaluated, is - directly added to the corresponding live load.
+P +P77
+ P6 + P6 Upper
Triangular
Panel
+P5 +P5
t
Ly
Lower
+P4 -1-P4
Triangular
+P3 +P3 Panels
+P
. 2 +P2
+121 -1-P1 Width --1
(a) Rectangular Panels (b) Triangular and Rectangular
in Right Bridge Panels in Skew Bridge
Parattelogrum
panels
&Y
known, the type of panetand its location in the grillage are determined.
The load in.the panel is now to be transferred to the contiguous nodes of
the panel in the form of equivalent loads.
A panel may have more than one wheel loads. In such a situation,. each
load is transferred to the nodes independently and the effects are summed
up algebraically at each node due to all wheel loads falling inside the
panel. Similarly, the nodes will receive loads froti other panels around it
and also any load which may come on it directly.
w L 2 ) m (WL WL
2 12
W/unit Length
o. (Z)
nrYerrel 3
1-*-- Li ___________ L
}M2 V
Figure 5.4 Transfer of Dead Load on the Nodes
Similarly, the torsion T per unit length due to the transverse eccentricity of
loads can also be distributed on the nodes treating the concerned longitudinal
grid lines as fixed.
As a sign convention, downward vertical force has been assumed as
positive. For moments, right hand screw rule has been followed.
consists of a uniformly distributed load on its area over part or full length.
Its magnitude etc. are discussed in Section 1.7 of Chapter 1. The footpath
live load may be distributed to the longitudinal grid lines situated in its
Vicinity according to the tributary area of each grid line, like dead load,
discussed above. The transfer of this live load to the nodes of the grillage
will also be done in the same manner as dead load. The bridge structure is
analysed for footpath live load on one side of the deck or ontoth sides and
over part length or full length and the resulting force responses are com-
bined with force responses due to other live loadings, if the combination
results in increase of responses.
ad ac
P P i = Lx L,, P
L., L,
P= P P= be bd P
i (5.2)
- 1,, L, 4L, Ly
For a lower triangular panel with load P as shown in Fig: 5.6, the equiva-
lent nodal vertical loads are obtained by first distributing thz, load on lon-
gitudinal 1-2 and end transversal 1-3 at edges E and F respectively in a
direction parallel to the transversal 2-3. The equivalent vertical loads at E and
F and on nodes 1, 2 and 3 are given by
(s (I)
Pr= P PF =s P
120 Grillage Analogy in Bridge Deck Analysis
Ly
X
2 3
y
P F
1 i 4,
Hb a
Figure 5.5 Statical Distribution of Loads in Rectangular Panel
Ly
d
b
3
TM) (
C
)
L.
1
(A)
Hence
C pE C p L:113
F =
Lx L, L.
P (4a L td)
=- 6 - 1 L PE = - p
L.,. 4.
P3 =-- La PF =---d , F
(5.3)
4.
Application of Lo'ads, Analysis, Force Responses and Interpretations 121
The above expressions for nodal forces can also be alternately obtained by
dividing the panel load P in ratios of the equivalent corresponding areas of
triangles formed and the total area of the triangular element, i.e.
area of A BPC
P 1=
area ofd ABC
area oft APC
P2
area of A ABC
By substituting in terms of the dimensions shown in the figure, the nodal loads
P1, P2 and P3 obtained from equation 5.4 will be same as given by equation 5.3
above.
The equivalent nodal loads in upper triangular elements (Fig. 5.3b) will
,be same in magnitude and direction as in lower triangular element given
above. In panels having parallelogram shape and containing the wheel load
P (Fig. 5.7), the equivalent nodal vertical loads at E and F and on nodes I,
2, 3 and 4 are similarly evaluated and given in equation 5.5.
P
E =. 2- P P F P
Y
= ad ac
_P
p
P
P
P
1
2 LXLy
= be p
P4 bd p
(5.5)
3 4
L , L y
I
.
b a
Ly
Figure 5.7 Statical Distribution of Load in Parallelogram Panel
illustrated and the resulting equations for vertical shear and two moments at
the nodes of the panels in each case are derived.
Vector forces on the diagrams are shown in the positive direction of axes
but proper signs are incorporated in the expressions given: Negative sign
preceding the equation will indicate that the force is in opposite direction to
the vector drawn on the diagrams. Also, the equations pertaining to fixed end
distribution of a point load, a bending moment and a torsional moment for a
beam are given below which will be useful in subsequent discussions.
a) Fixed End Distribution of a Point Load P for a beam AB (Fig. 5.8a)
b 2 (3a + b) D D-
a2 (3b + a) P
I 5
I
A =
ab2 a2 b
M= M= P (5.6)
A L2 L2
M
A 61 I . OmB
1____________
4 _______ a w1-4
L
b .).
,.
PA i
1 P
B
MA rm MB,Ai
t's 11/4 B
a _________________________________________________ I - b ______.1
a
_________________________L .i
Pt: PB
(b) Bending Moment Distribution
(c) Torsional Moment Distribution
Sign Convention : Moment Force Torsion
TA a /6- lAl b
I
Figure 5.8 Fired Ended Beam: Distribution of Forces at the
L Dy Ly
124 Grillage Analogy in Bridge Deck Analysis
M y2 M Y7
MX2 11 ;P 11 MX3
A
t - MXF
MX E
P
F
1 4
M
e )E------ba
l
___1fhP4
_____________ LY
Myi
Figure 5.9 Load Distribution in Rectangular Panel
a2 bn ab2M
p
(5.9)
=
2r
y
'14 L AC
L2 4
Again treating the member 1-2 as fixed at its ends, PE will induce at node
1, vertical force P1 and moment and /1//..x will cause moment Mxi in the
directions shown and their values can be written as
d2 (3 c + d)
P
= Lg PE
x
ar Cd2
Application of Loads, Analysis, Force Responses and Interpretations 125
Similarly for other nodes of the panel, the vertical force and the two
moments can be evaluated as follows:
At node 2
C2 (3d +c)
P2 ' 3 PE
Lx
2 d
M =MXE M =-- r (5.11)
E
Lx L2
At node 3
g (3d+ c)
3 = L 3 F
c2d
M r ri (5.12)
13 L 2F
At node 4
d2 (3c +d)
P4 L3 F
c d2
Figure 5.10 shows a lower triangular panel of the grid acted upon by a load P. The load P is
first distributed to members 1-2 and 1-3, parallel to the transversal 2-3 at E and F respectively.
The vertical forces PE and PF and associated moments MXE and MXF at E and F are given by
d
PE=(s d) 2 (2d s) P P F 2 (3 s 2 d ) P
s3 S3
(S d) 2 d p d2 (s d) p
M XE 2 M XF 4
(5.14)
Now, PE will cause vertical forces P. and P2 and moments
M yi and Mn and the moment MXE will induce moments Alm
and Mr at nodes 1 and 2
126 Grillage Analogy in Bridge Deck Analysis
6\14'\
- e:
\
M
Y2
.4___ Nwr
Mx2 LY _______________________________________________________________________________________
1.R2
C
X
Lx
a
1(A)
.4,45\ Myl
Figure 5.10 Load Distribution in Lower Triangular Panel
P =C2 (3a + c)
L3 E
c
MX1= MXE AlYL = c2 a P (5 .15)
x
At node 2 due to 3
a2 (3 c + a)
P 2 or P E t = LS
=7J
a
M o r M X S =LM X E
s
2
c
At node 1 due to PF
C2 (3a+ c) c2 a
Pim= L3 PF
M =
"I) sin a PF (5.17)
P3(1) = a2 (3c + a) a2 c
PF
M3(1) LPF (5.18)
I sin a
Now MXF acting in the direction of X-axis is resolved into components MXF
sin -a and MXF cos a in the directions parallel and perpendicular to member 1-
3 respectively. MXF sin a will cause torsional moments M1(3) and M3(3) at
nodes 1 and 3 and MXF cos a will cause bending moments M1(.2) and vertical
force Pim at node 1 and bending moment M3(2) and vertical force P3(2) at node 3.
The expressions for these force resultants are given below
a
M1(3) = L MXF sing M3(3) L = MXF sin a
c (2a c)
M1(2) = L MXF cosa
!
a (2c 1 a)
(5.19)
M MXF cosa
3(2) =
:"1 .
Now at node 1, the total vertical force PA, moment about x-axis MxA and
moment about y-axis M, will be,
128 Grillage Analagy.in Bridge Deck Analysis
PA = PI + P1(1) + PI(2)
c2 (3a+ c) c2 (3a+ c) 6ac
PF M xF sina cos a
L3 P+x L x
3
L x
L 2
L 2F
Ls2 Lx
(5 . 2 0 )
Similarly at node 3, the total forces Pc, Mxc and M 1 will be,
P = P + PC 3(1) 30.1
a 2 (3c+ a) , bac
rF + MxF sin a cos a L3x
Mxc = M3(1)cos a + M3(2) cos a +.M3(3)sin a
a2 c a (2c a)
= r2 FcoLa costa MxF + a MxF sine a
"L L! Lx
cr La(2ca)
= +
2
M xF sin a cos a + a MxF sin a cos a (5.21)
x
Similar expressions can be derived for upper triangular panel also (Fig.
5.11). The distribution of vertical forces, bending and torsional moments on
nodes of the panel are shown.
The resultant vertical forces, at joints A, B and C will remain same as in the
case of lower triangular panel. Also, the bending and torsional moments at the
three nodes will be the same in magnitude but opposite in direction. (c)
Parallelogram Panel
The load in a parallelogram panel can be distributed to its nodes in the
similar manner as in the case of rectangular or triangular panels discussed
- - .
limxi
My1
__ y
41-
;'"s
&-
ty
e 4.ar Mx
MXF
4}h>si>
CV') 17,1
/4.
1
P
3 `ry 4N"' .\ 3(C)
4'62 b d 2 (B)
2
_
_
_
__
_
__
_
_
__
_
__
_
_
__
_
__
_
_
__
_
__
_
_
__
_
__
_
_
__
_
__
_
_
__
_
__
_
_
__
_
__
_
_
__
_
__
_
_
__
_
__
_
_
__
_
__
_
_
__
_
__
_
_
__
_
__
_
_
__
_
__
_
_
__
_
__
_
_
__
_
__
_
_
__
_
__
_
_
__
_
__
_
_
__
_
__
_
_
__
_
__
_
_
__
_
__
_
_
__
_
__
_
_
__
_
__
_
_
__
_
__
_
_
__
_
__
_
_
__
_
__
_
_
__
_
__
_
_
__
_
__
_
_
__
_
__
_
_
__
_
__
_
_
__
_
__
_
_
__
_
__
_
Ly_________________________________________________________________"4 M
Y
11M x2
a2 (3 b + a) b 2 (3a+ b)
E L3 PF Ls
P
2
M =a b P sec a b2a
MF 2 See a
E L2.1. L
=
a2 b
P b2 a
tilXE 2 MXF =
L L2
130 Grillage Analogy in Bridge Deck Analysis
P3 mx3
0
Y4
tlx2
P2
my 2. MF MXF
d 0 YF_
PF
Lx
PO
ticlx4
P.4
1
01
p H-- b .i* a 1-1 x
Ly
2b ab2
11.11.E=-- P tan a m= P tan a . (5.22)
L2 YF L
d
xk1L E
c d2 d (2 c - d)
PE YE (5.23)
Lt L;
Application of Loads, Analysis, Force Responses and Interpretations 131
P2 c2 (3d + c) 6cd
P E 3+ Y M
L 3 E 1, E
x
L a .
C 119
M XE
x ' '
e2 d c (2d
2 c) mYE
Y 2
=L 2x E L
(5.24)
(3d+c)
C2 6cd
P3 = L 3 I F + L3 MYF
Lx
Mx3 L i v i x"F
x
Loading Axle
p
1"1"-- YMAX
c
_FH44,R- rncn-751---1
Left Most Right Most
Grid Line B _ Grid Line
(i)
++++Cross-section
X=L ++++ 1
XM 1 1 t+
0 CI)
B 4T*-
0 0
;41
(0,0)
+ 4I+
+ + + +
0 0
()Lett Most Wheel of Leading Axle 0 0
Right Most Wheel of Rear Most Axle 0 CD
(ii) Right Bri dge Plan
Figure 5.13 (a) Initial and Final Positions of Loading to Cover Entire
Right Deck
The total number of positions which a live load may occupy on the deck and for
which the bridge will be analysed. are:
Application of Loads, Analysis, Force Responses and Interpretations 137
F"' 1MAX
Loading Axle
c -0 1 C
4141-p-n'
X=L
XMAX
(0,0)
u
I
Figure 5.13 (b) Initial and Final Positions of Loading to Cover Entire Skew
Deck
138 Grillage Analogy in Bridge Deck Analysis
various live loads on the bridge deck through prescribed steps along the span
and across the carriageway width.
Loading (1) Dead load
Loading (2) Class A Two lane
Loading (3) - Class 70R Train : Col 1
Loading (4) Class 70R Train : Col m
Loading (5) Class 70R Bogie : Col 1
Loading (6) Class 70R Bogie : Col m
Loading (7) Class 70R Track
Loading (8) Class A Single lane
Loading (9) Specified by user
Loading (1) represents dead load. Loading (2) to Loading (9) represent
various types of live load.
It is built-in the program that live load results are given inclusive of dead
load results. If only live Ibad results are desired, then in the Input data, dead
load can be given as zero.
.... ..,
_ l
mertt
\Griltage Output
Only one value of the shear force for a member of the grillage is obtained
from the output and the same may be used in design as such. Similarly,
maximum reactions printed, are taken as design values for reactions at
supported nodes.
The torque per unit width in a true orthotropic slab is same in orthogonal
directions, however, it is often different when read from 2TilIage output.
The torque per unit width at any point should be taken as the average of
these two outputs.
Application of Laads, Analysis, Force Responses and Interpretations 141
The design of the section in bending should be based upon the principal
moments mi and in,. Referring to Fig. 5.15, the principal moments mi, m2 and
their deviation at the nodal point of interest is obtained from the values of
bending moments per unit width mx and my, in X and Y directions and the
corresponding averaged torsional moment per unit width mx,employing the
formulae
2
nix+ my {(mx my )
m1.2 2
- + 2 2 +m 2
(5.27)
ttan2 0 = m xy
mx my
m1
I
2
My
Figure 5.15 Computation of Principal Moments and Their Deviation from
X-Y Axes
(M) 70/,)4
(My) = (M)0 + 7(Mx),3 (5.28)
where (M )0and (M)0 are the responses obtained by grillage analysis for
which Poisson's ratio = 0 and (Mx) and (My) are the relevant corrected
moments for a finite value of Poisson's ratio.
Due to dispersion of load, the actual area of application of load is larger
than the area of contact of the wheel with the slab. If this application area is
larger than the grillage mesh, the load can be assumed to be sufficiently
dispersed for the grillage to reproduce the distribution of moments through-
out the slab. No further modification of moments is necessary. On the other
hand, if the area of application of the load is small compared to the grillage
mesh, no worthwhile information will be obtained about the local high values
under the load, though the grillage distributed moment field will simulate
that in the deck. The additional moments due to high local curvatures can be
obtained for the area of slab within the grillage mesh from influence charts
due to Pucher [6].
4
r
large jumps in moments at the nodes. This is because of the transfer of the
torsional moment in the transverse members at each joint to bending mo-
ments and shear forces in the longitudinal member [3]. The true design
bending moment is obtained by averaging the bending moments on the two
sides of each joint. The top and bottom slab stresses are calculated from these
average moments.
The transverse bending moment in the grillage member is equivalent to
the opposed transverse compression of the top slab and tension of the bottom
slab (or vice-versa) due to transverse flexure without Tlistortion. In
narrow decks, it is usually very small compared to the longitudinal bending
moment (except in the diaphragms) but in wide decks it can be large spe-
cially near skew supports. A grillage output transverse moment diagram has
a saw-tooth shape like the longitudinal moment diagram and the top and
bottom slab stresses are calculated from the average moments.
The slab bending moments are derived from the shear force in the trans-
verse grillage members. The fraction of this shear force carried by each of
the top and bottom slabs of the cell are assumed to be proportional to the
flexural stiffnesses of the slabs. Assuming also that points of contraflexure
lie midway between the webs, the moment at each end of a slab is simply
the shear force it carries multiplied by half the distance between webs. The
transverse slab moment in the cantilever can be taken directly from the
grillage output since this member is not representing the cell. Local moments
are again obtained using influence charts of Pucher t6]. The total design
moment in cellular deck is obtained by adding the slab moments of the cell,
cantilever moment and the local moments.
The torsion-shear flow in the slabs must be calculated from the average
torque per unit width of transverse and longitudinal grillage members.
The grillage output shear force represents the total shear force in each web
of the deck.
Grillage analysis ignores the effects of Poisson's ratio on the interaction of
longitudinal and transverse moments. This neglect of Poisson's ratio
introduces little error in narrow decks while in wide decks with little stiffness
against cell-distortion, the calculated transverse moments can be con-
siderably in error, if moments are small and the Poisson's ratio is significant.
However, since the concrete has relatively low Poisson's ratio of ap-
proximately 0.15, its neglect may not result in any appreciable error.
Reviewing the discussions contained in Chapters 4 and 5, taken together,
on the grillage analogy method, it can be seen that the procedure is straight
forward and is also amenable to an enaine.er who may not he well conversant
with higher mathematics and numerical techniques. The setting out of
Application of Loads, Analysis, Force Responses and Interpretations 145
the grid lines and the evaluation of elastic properties of the members can be
handled manually following the guidelines provided above. The determi-
nation of the dead loads etc. is simple. The remaining steps in grillage
analogy i.e. transfer of loads to nodes, formulation and inversion of matrix,
solution of equations, evaluation of design responses etc. are done on com-
puter using: a suitable software package which are not difficult to find
now-a-days. Preparation of input data is also simple (as will be shown later)
and can be done mechanically even without understanding the computer
program fully However, the interpretations of results and its modifications
due to local effects pertinent to a particular type of deck may present some
difficulty and will need some understanding of the structural behaviour of
the deck 'in question.
A general grillage analysis computer program developed for dead load
and various types of live loads is discussed in the next chapter.
REFERENCES
1. BAKHT, B. and JAEGER, L.G., "Bridge Analysis Simplified", McGraw Hill, New
York, 1985.
2. GERE, M.J. and WEAVER, W., "Matrix Analysis of Framed Struture", D. .Van
Nostrand Co., USA, 1965.
3. HAMBLY, E.C., "Bridge Deck Behaviour", Chapman and Hall, London, 1976.
4. JAEGER, L.G. and Bmarr, B., "The Grillage Analogy in Bridge Analysis", Canadian
Journal of Civil Engineering, 9(2), 1982.
5. JAEGER, L.G. and BAKHT B., "On the Analysis of Slab Bridges by Grillage Anal-
ogy", Proc. Annual Meeting of the Canadian Society for Civil Engineering,
Montreal, May 1987.
6. PUCHER. A., "Influence Surfaces of Elastic Plates", Springer Verlag, Wien and
New York, 1964.
7. RUSCH. E.H.H. and HERGENROEDER, A., "Influence Surfaces for Moments in Skew
Slabs", Werner-Verlag, Dusseldorf, 1969.
8. WEST. R., "Recommendations on the Use of Grillage Analysis for Slab and
Pseudo-Slab Bridge Decks". Cement and Concrete Association, London, 1973.
9. WEST. R., "The Use of Grillage Analogy for the Analysis of Slab and
Pseudo-Slab Bridge Decks", Research Report 21, Cement and Concrete
Association, 1973.
10. IRC 6-1987, "Specifications and Code of Practice for Road Bridges-Loads and
Stresses", Section II, Cement Concrete (Plain and Reinforced), Indian Roads
Congress. New Delhi - 1990.
11. IRC 21-1987, "Specifications and Code of Practice for Road Bridges". Section
III, Cement Concrete (Plain and Reinforced), Indian Roads Congress. New
Delhi, 1988.
-- -
Chapter 6
Computer Program
6.1 INTRODUCTION
A Computer Program based on grillage analogy is developed in FORTRAN
for the elastic analysis of bridge decks, covering right and skew layouts. The
program is compatible with MS-DOS environment in its present form but it
can also be compiled to run on mini or main-frame computer. The program
is modular and includes different subroutines to execute specific operations.
Additional features and modifications can be easily incorporated for a future
modular expansion. Sufficient number of explanatory comment statements
have been introduced for better and quicker understanding of the program.
However, the program is written with the premise that the reader has some
prior exposure to computer programming. But if one does not want to go
through the drill of understanding the program modules completely, he can
still prepare the input data and can analyze the bridge.
This Chapter deals with the important features of the program. The
Program Manual consisting of listing of variables, sign conventions, main
and subroutine segments and flow charts is given. The Users Manual
comprises data input and output formats. Limitations and scope of the
program are also discussed.
The listing of the program is given in Appendix II at the end of the book.
A diskette with executable file of the program compatible to IBM PC along
with input data and result files for all the worked examples illustrated in
Chapter 7, is available and can be ordered separately.
6.3.1 Variables
Some of the variables used in the program for input, output and also
those used in more than one subroutines are listed below:
148 Grillage Analogy in Bridge Deck Analysis
Left-most Longitudinal
Grid line
(
Mx 'ex
r __ My ,6y
Z (downward) Z (downward)
(a) Coordinate System (b) Sign Convention Figure
6.1 Coordinate System and Sign Convention
Call "COMPARE2" to Obtain Final Design Responses by Comparing all Load Types
STOP
X.4
.L. 1:1317
Computer Program 153
(START)
INPUT
Grid type, Skew angle, No. of members,
No. of nodes, No. of types of groups of
members, E, G, No. of longitudinal grid lines,
Spacing of longitudinal grid lines, X-
Projections of bridge beyond support lines,
Node nos. of one representative element of
the group, Total no. of elements in the group,
X and Y lengths of the element, I and J of the element,
Nodes nos. of remaining elements of the group.
Bearing
Stiffness
INPUT
/ I N P U T
of kerb from left most
longitudinal grid line.
Call Subroutine
"DEAD LOAD"
YES
INPUT INPUT
Vertical and Torsional udl Vertical load
loadings on longitudinal grid at each node due
lines due to dead loads to dead load
Computer Program 155
Assemble global
stiffness matrix
in half band-
width form
Decomposes the
Stiffness matrix Call Subroutine
into upper & lower 4r-- "DECOMP"
triangular matrices
/ INPUT
"CH3"
YES
STOP)
NO
156 Grillage Analogy in Bridge Deck Analysis
INPUT
TotalTotal no. of live load
cases to be analysed
DO 1000
Loop for solving all
live load cases
I N P U T
Load case type
"LCASE"
NO
INPUT
I n i t i a l & f i n a l X a n d Y
c o o r d i n a t e s o f l e f t m o s t f r o n t
w h e e l o f l o a d i n g , I n c r e m e n t s i n X
a n d Y c o o r d i n a t e s .
Calculates impact
Call Stibroutine factor according
"IMPACT" to type of loading
and span.
1Call Subroutine
"PAIDRB"
A
r YES
:han _____________________________________________
NO
Triangular Rectangular
Are
all wheel
Next loads analysed
0 4
w h e e l t ___ NO and load vectors
added
Spi
Inflates the deformation
vector for rigid bearing
case with zeros at restrained
degrees of freedom.
4,
Calculates member end
forces using member
stiffness matrices and
deformation vectors.
4
Calculates support reactions
by using equilibrium
conditions at the supported
nodes.
Are
all load
NO positions for
urrent load type
exhausted
YES
INPUT
IK2
160 Grillage Analogy in Bridge Deck Analysis
NO Is
IK2=I
YES
1' Prints the maximum
Call Subroutine values of member end
"WRITEl" forces; deformations,
reactions etc. for the
current load type.
Are all
NO load
Next
load case exhaus
L
17.. 7, ,
. 1 1 1 1 1 1 . m w . - - s - - - - . : . . . _ , -
162 . Grillage Analogy in Bridge Deck Analysis
In load case 9, the *user specifies the details of loading, namely, the
number of wheels and their coordinates. The wheels are numbered
beginning from left most wheel of leading axle, successive wheels
on the same axle in the positive direction of Y, then the wheels on
next following axle and so on.
8. Subroutine WLOAD
This subroUtine describes magnitudes of load on all wheels for each
type of live loading. The minimum alliArable distance of center of
left most wheel from Kerb is set in. All the wheels are taken as point
loads acting at the centers of wheel contact- area. Class 70R Track
loading is described as ten point loads per track, the load on each
point being 3.5 tonnes at spacing of 457 mm as shown in Fig. 1.13.
For User's Specified loading, the' oad on each wheel inclusive of
impact and minimum kerb distance from center of left most front
wheel are given as input.
9 . Subroutines PAIDRB and PAID1
When grid type is skew and the transversals are parallel to the supports
forming parallelogram panels or when the grid is right i.e. skew angle is
zero, the subroutine PAIDRB is called. On the other hand, when the
skew angle is not zero but the transversals are orthogonal to the
longitudinals, the program subroutine PAID1 is called. These subrou-
tines use the spacing between longitudinal grid lines and the lengths of
elements taken on longitudinals, to calculate the boundaries of the
panels with respect to node 1 (specified earlier) and then compare the
coordinates of wheel load with these panel boundaries to identify the
panels containing the wheels. These subroutines are further discussed
below.
a) Subroutine PAIDRB
For each wheel, the Y-coordinate, measured from the origin i.e. from
the left most longitudinal grid line, is checked. If the Y-coordinate is
less than minimum permissible distance, defined by the distance
(Akerb+SCLMIN), the subroutine terminates with position identifica-
tion character as YL. Also, the subroutine terminates if the Y-coordi-
nate exceeds Y mar (Table 5.1) with position identification character
YG. Similarly, if the X-coordinate is such that the wheel is not on the
span, the subroutine terminates with position identification character
as XL or XG as the case be. The panel is now identified as 'P.' i.e.
rectangular or parallelogram. By using the coordinates of wheels, the
panel containing the wheel is identified by its contiguous nodes MNI,
Computer Program 163
y MN11
MN1 MN 22 MN2
The dimensions of the panel and the position of the wheel on the panel
with respect to MN1 are calculated and stored. Distinction between
lower and upper triangular panels as shown in Fig. 6.4(c, d) are also
made using position identification characters as 'XG' or 'XL'. If the load
is out of the permissible range in Y-direction (nearer to Kerbs than
permissible), the program stops. If a particular wheel - load is
longitudinally out of span then it is ignored.
10.: Subroutines LDISTR and LDISTT
These subroutines distribute the load from the panel onto the nodes in
the form of vertical loads and moments in X and Y directions. The
subroutine LDISTR is used for distributing the load to the nodes of
rectangular as well as parallelogram panels while the. subroutine
LDISTT distributes the load to the nodes of triangular panels. The
final expressions for nodal load vectors for rectangular and parallelo-
gram panels are given by equations 5.10 to 5.13 and 5.23 to 5.26
respectively and for triangular panels are given by equations 5.16 to
5.21. These have been used in the program with appropriate signs.
11. Subroutine SOLVE
This subroutine uses the previously decomposed stiffness matrix and
the load vector to solve the set of simultaneous equations by backward
and forward substitution procedure and gives the deformations corre-
sponding to three degrees of freedom. If the support type is rigid, the
deformation vector is 're-inflated' by zeroes at the support nodes. This
subroutine is called directly if the supports are flexible. In case of rigid
supports, the load vector is first reduced in size by eliminating elements
corresponding to degrees of restraints at supported nodes and then the
subroutine SOLVE is called.
12. Subroutines COMPARE1 and COMPARE2
These subroutines are used to produce the output for the maximum
live load design responses.
a) Subroutine COMPAREI
The deflections, slopes, moments, shears and torsions at the nodes and
the reactions at the support nodes are compared for all placings of the
same live loading and the critical values of these force responses are
retained along with the corresponding load positions. Along with a
critical response, wherever relevant, the corresponding values of other
related responses, are also retained. For example, along with maximum
shear force corresponding torsional moment, along with maximum
torsional moment corresponding shear forcc,
Computer Program 165
a check. This is followed by the analysis results. The analysis results consists
of force responses and deformation responses.
The following is the sequence in which the output result is obtained:
I. Grid details as per input data after suitably titling these.
2. Dead load details as per input data.
3. A list of different types of live loads for which the analysis can be.
done and reference numbers of these live loads.
4. Responses under dead load as follows:
(a) Shear forces and torsional moments in each beam element.
(b) Sagging bending moment at each node on its right and on its left.
(c) Hogging bending moment at each node on its right and on its left.
(d) Rotations about X and Y axes and vertical deflections at each node.
(e) Reactions, rotations about X and Y axes and vertical deflections at
each support. Deflections at supports will occur when grid is
supported on yielding supports.
5. Responses under individual live load together with impact load and
dead load as follows:
(a) Span for impact factor calculation and impact factor taken for the live
load.
(b) Type of live load through its reference number, initial and final posi-
tions of reference wheel through its X and Y coordinates and step
lengths in X direction and Y direction respectively. Reference wheel is
the left most wheel on leading axle of the loading system.
(c) Shear force and torsional moment for each beam element
(i) maximum shear force, corresponding torsional moment and posi-
tion of reference wheel;
(ii) maximum torsional moment, corresponding shear force and
position of reference wheel;
(d) Maximum sagging bending moment at each node:
(i) at right of the node and corresponding position of reference
wheel;
(ii) at left of the node and corresponding position of reference wheel;
(e) Same as at (d) above but for maximum hogging bending moment.
(f) Maximum rotations about X and Y axes and maximum deflection at
each node.
rom.ar.tinri :et earth cttrwsnrt with Correspondingdeflection,
REFERENCES
1. GERE. M.J. and WEAveR, W., "Matrix Analysis of Framed Structure', D. Van
Nostrand Co., USA, 1965.
170 Grillage Analogy in Bridge Deck Analysis
2. NAYAR, K.K., RAGHUPAIHI, M., SEETHARAMULU, K. and &RANA, C.S.,
"Computer Aided Bridge AnalysisA Software Based on Grillage Analogy",
International Conference on Bridges and Flyovers, Hyderabad, February 1991.
3. SHARMA, K.G., &R ANA, C.S. and AGRAWAL, A.K., "Automated Analysis of
T-Beam Bridges Using Grid Method", Proceedings of Second International Con-
ference on Computer Aided Analysis and Design in Civil Engineering, U.O.R.,
Roorkee, January 1985.
4. "Ontario Highway Bridge Design Code", Ministry of Transportation and Com-
munications, Highway Engineering Design, Downsview, Ontario, Canada, 1983.
5. "Standard Specifications for Highway Bridges", American Association of State
Highway and Transportatiion Officials (AASHTO), Washington, D.C., 1977.
6. IRC 6-1987, "Standard Specifications and Code of Practice for Road Bridges -
Loads and Stresses", Section II, The Indian Roads Congress, New Delhi, 1990.
Chapter 7
Illustrative Examples
7.1 INTRODUCTION
This Chapter aims at illustrating the use of the grillage analogy method
discussed in earlier chapters through worked examples. A set of six examples
have been chosen covering different types oDbridge superstructures.
Converting an actual bridge deck to an equivalent grid is not without
pitfalls and needs special attention. This conversion involves setting of longi-
tudinal and transverse grid lines and assigning flexural and torsional inertia
values to various grid members. This has been done for all the bridge
examples presented in this Chapter along with brief discussions. It is ex-
pected that these examples and the detailed recommendations given in
Chapter 4 together, will enable a prospective user of the grillage method to
convert a bridge superstructure to an equivalent grid with requisite accuracy
and confidence.
The computer program 'GABS' discussed in Chapter 6 is a versatile
program. These examples are used to explain the use of the aboVe
program when the loading is IRC loading or a user-specified loading or a
stationary loading.
A user of the program at some stage would like to know the sensitivity
of a solution to the size of the grid mesh, i.e. how will the response result
alter, if instead of a normal sized mesh, as recommended in Chapter 4, a
finer or coarser sized mesh is adopted. Example 1 of a right slab bridge is
devoted to investigate this.
End-reactions are non-uniform for a skew bridge even under a
uniformly distributed loading on the deck, such as 'self load'. Example 2
is of a skew slab bridge and is chosen to illustrate this variation, when the
slab is supported on bearings either rigid or flexible.
Computation of Inertia for a voided slab is not straight forward.
Example 3 of a voided slab bridge illustrates this computatiori, as well as
the handling of 'non-uniform dead loading' along the longitudinal grid
172 Grillage Analogy in Bridge Deck Analysis
In Example 4, a T-beam bridge has been analysed for dead load as .well as
live load specified by the Ontario Highway Bridge Design (OHBD) Code of
Canada.
In example 5, a skew T-beam bridge is analysed for different IRC
loadings and the envelope diagrams. for B.M and S.F are presented.
A box-girder bridge develops torsional and distorsional Warping under
loading. The bending stresses are also affected by shear lag. All these effects
are ignored in the program "GAB" presented in this book.
In example 6, a box-girder bridge has been analysed by 'GABS' as
well as by Finite Element Method and the results are compared.
72 ILLUSTRATIVE' EXAMPLES
Examples 1, 2 and 3 pertain to a right slab, skew slab and a voided slab
bridge deck respectively. Analysis,of slab-on-girders bridges both right and
skew have been illustrated in examples 4 and 5 whereas example 6 dis-
cusses a box-girder bridge.
All the bridges are simply supported and have 7500 mm carriageway
width. The bridges are of reinforced concrete. The Young's modulus of
elasticity and Shear modulus of elasticity are taken as 2.0 and 0.87 tlmm2,
respectively. Wherever dead loads are considered, density of reinforced
concrete including density of wearing coat is taken as 2.40 dm'. In the
input data, the loading type is designated by a number which refers to the
type as given. in Section 5.5.3 of Chapter 5.
COmplete Input and Output files of all the six examples of.this Chapter
and one example of Chapter 3, with the two programs 'GABS' and
'GRID!, (in a ready to use form) is available on a diskette. The diskette
can be ordered from the authors through the publisher.
7.2.1 Example 1: Right Slab Bridge
A two-lane right slab bridge (Fig. 7.1) is chosen for this example. The
equivalant grid is shown in Fig. 7.2 and is referred to as normal mesh in
further discussions. It consists of seven longitudinal and seven transverse
grid lines. Recommendations given in Chapter 4 are followed in arriving at
the grid pattern. The bridge is analysed for four different types of IRC live
loadings alongwith corresponding impact factors. The IRC live loadings
chosen are (i) Class A-2 lane loading, (ii) Class 70R Train Col. loading,
(iii) Class 70R Bogie Col. '1' loading and (iv) Class 70R Track loading.
These loadings are moved on the bridge one at a time in suitably chosen
intervaIsboth longitudinally as well as transversely so that the load
traverses the entire length and width of the deck. For this example, an
o b .
P
Wearing Coat .4600
470lieroo
0
7500
L
300
9000
300
interval of 500 mm has been chosen for longitudinal movement of all types
of loadings. In transverse direction, the intervals are so chosen that the load
traverses the full deck width in 5 or 6 steps.
The input data is given in Table 7.1 and selective output is given in Table
7.2. It may be noted that maximum longitudinal bending moment is obtained
under Class 70R Track loading. The maximum bearing reaction'per unit
width for exterior grid lines is under ClassA.-2 Jane loading and for
utidr.-.1 701-: Train Col. T loading. This shows that for
te r
700
_0
7 000
7500
H 60 1350 =8100
300
42 49 35N,
7
fi 13 20 27 34 41 - he
5 12 19 26 33 40 47 0 o
co
%0
IL
11 18 25 32 39 46 ti+
3 10 17 24 31 38 45
2 9 16 23 30 37 44
1 8 15 22 29 36 43
a a
Outer Outer
Transversal
N-10504
T Longitud
Figure 7.2 Normal Grid Mesh and Sections of Various Grid Members
Illustrative Examples 175
f a particular bridge any one particular type of IRC live loading does not give
all the critical responses.
300 300
1 2 12 1500 0 6.48E+10 -12.97E+10
2 3 3 4 4 5 5 6 6 7
43 44 44 45 45 46 46 47 47 48 48 49
8 9 3 0
1500 0 3.86E+10 7.72E+10
9 10 10,
11 II 12 12 13 13 14
15 16 16
17 17 18 .18 19 19 20 20 21
22 23 23
24 24 25 25 26 26 27 27 28
29 30 30
31 31 32 32 33 33 34 34 35
36 37 37
38 38 39 39 40 40 41 41 42
1 8 1 2
0 1350 3_0E+10 6.0E+10
8 15 15
22 22 29 29 36 36 43
7 14 14
2 9 30 21 21 28 28 35 35 42 42 49
9 16 16 0 1350 4.29E+10 8.58E+10
3 10 10 23 23 30 30 37 37 44
4 11 11 17 17 24 24 31 31 38 38 45
5 12 12 18 18 25 25 32 32 39 39 46
6 13 13 19 19 26 26 33 33 40 40 47
14 3 1 20 20 27 27 34 34 41 41 48
300.
1 7 8 14 15 21 22 28 29 35 36 42 43 49
0 0
0
1 0 0 0 0 0 0 0 0 0 0 0
1
0
0
4
2
Y Y
.
90
0
84 49 4; 2.0 0.87 7 1350 1350 1350 1350 1350 1350
0
0 .1705. 500. 385. 22400. 4015.
Contd.
. .`
0
5
0 1705. 500. 385. 10220. 4015.
1
0
7
0 1920. 500. 440. 13113. 4120.
1
0
1
0
S
This example is further used to study the effect of size of the mesh
formed by grid lines on various force responses. Figure 7.3 shows a coarser
mesh for the same bridge where the number of longitudinal grid lines have
been reduced from 7 to 5 but the number of transverse grid lines have been
kept the same. Figure 7.4 shows a finer mesh for the same bridge where the
number of transverse grid lines have been increased from 7 to 11 but the
number of longitudinal grid lines are kept same as in normal sized mesh of
Fig. 7.2. Both these grids are analysed for the same four types of IRC live
loadings as above, keeping the longitudinal and transverse intervals for
various IRC loadings same as in the analysis of grid of Fig. 7.2. The input
(b) Maximum bearing reactions and corresponding bearing reactions per unit width
Reference Exterior Grid Line Central Grid Line
Grid Load Reaction R per unit Mean Load Reaction R per unit Mean
Type R width (r) r Type R width (r) r
(t) (t/m) (t/m) (t) (t/m) (t/m)
Normal 2 7.23 7.42 3 17.50 13.01
Coarse 2 9.59 7.30 7.41 3 26.02 12.tb i2.91
Fine 2 7.32 7.51 3 17.38 12.87
Illustrative Examples 177
4.6001.r-
700 1000
$
4 @ 2025 8100
1 .,,,7 s' v Lc, V J aS
-r-
13 20. 27 . 34
6
12 : 19 26 33
5
0
0
0
4 11 18 " 25 32 171
10 17 "24 31 eU
0
3
9 16 23 30
2
1 8 15 22 29
A A A A P
Inner Longitudinal Outer Longitudina L
1=5.79 E10
6e 1 35 0 = 81 0 0
11 v 22 33 44 a 55 v 66a 77
10 21 32 43 54 65 76
.
9 20 31 42 53 64 75
19 30 41 52 63 74
70 X 900..9000
7 18 29 40 51 62 73
6 17 28 39 50 61 72
5 16 27 38 49 60 71
15 26 37 48 59 70
3 14 25 36 47 58 69
2. 13 24 35 46 57 68
1 12 23 34 45 56 67
a
-
1 = 2. 5 7 E 10 . 4 7 50 1. - - T 1=214E10
700 Jf. 4 29E10
70 0 J = 5 1 5 E 1 0
900 k --pl.
3001.r--
10490
7
0 9 24 5
(a) PLan at mid-thickness of Slab and EquIvaLent grid
50
Q..
. 40
Rigid Bearings
oc 30
Flexible
20
cc 10
0 9 18 27 36 45
Bearing Node Number
TABLE 7.5: Bearing Reactions Under Dead Load and Dead & Live Loads
Combine
Bearing Dead Load Alone Dead and Live Loads Combine
Node Reaction (t) Min. Reaction (t) Max. Reaction (t)
Flexible Rigid Flexible Rigid Flexible Rigid
Bearing Bearing Bearing Bearing Bearing Bearing
9 or 37 44 51 44 49 68 74
18 or 28 31 23 31 22 55 55
27 or 19 28 29 28 29 52 61
36 or 10 24 27 24 27 42 47
45 or 1 13 11 11 9 20 16
Illustrative Examples 183
7.2.3 Example 3: Voided Slab Bridge
A voided slab bridge is chosen for this example. The voids are running
through the central part of the span. In support regions where shear is
large, solid section has been provided. Figure 7.7 shows the bridge details
and the proposed grid arrangement is shown in Fig. 7.8.
The example is chosen to explain (a) how inertia values are to be com-
puted for such voided sections and (b) how non-uniform dead load along
the longitudinal grid line is to be dealt with.
The inertia values are computed on the principle that slab has same
Liu8 , 8
supportest 4
1 Mis.mmiM
14441 .
I l
r - , ._,
in Fin r : !---. 7 '
ri
I 1 1 I 1 iI I1 11
il I
II 1 1 I
I I 1 I II
I I I I II : I I
III I
1ii I I I N.
1I I II :I
C
I Ig 0 -
iI1
I , , I I 1 ec. Q
1
_11.11i r 1 8
iilriii Hi!!I I N,
it
1 11 1;111 i I
I II Ii ' 1 I 1 1 II I
I i. I II I Ill!! 1 I 1 1 I
LJLJLJULJLJUII I ..
co o '
0
90
0
7500 9C11,_ISuf:port
_
F()11(3 01 0 010 0 0
H- 7 x 960 u 6720 .4 1-.
480 480
90 1830 1920 1920 1830 90
Figure 7.7 Voided Slab Bridge
184 Grillage Analogy in Bridge Deck Analysis
7
V
v5
14 21 28 6
13 20 27 3
0 0
0 N
- 2 11
5 0
12 19 26
0
N
(0
3
4
11 18 25 .
30
1
10 17 24
9 16 23
8 15 22
1 A 2 9
i+-1830-441 2 x 19 2 0 1830-44
= 3840
. M . 1 1 1 1 . Figure 7.8 Grid Relating to Voided Slab Bridge
inertia per unit width in the two mutually perpendicular directions. This is
true only if void depth is upto 60 per cent of the overall depth. For voided
section, the moment of inertia per unit width of the slab is computed from
the section showing the void/voids and the same is assumed per unit width
in the orthogonal direction. Torsional inertia .1 of the section is computed
from moment of inertia 1 adopting the formula .1 = 21 for all slab sections.
Table 7.6 gives the sections and their I values. Moment of Inertia values of
members belonging to groups 3 to 10 have been calculated on the basis of I
per unit width (i) of the members of groups I and 2.
The dead load alone the longitudinal grid line is non-uniform as in the end
regions the section consists of solid slab and in the central region the section
consists of voided slab. Such non-uniform loading is dealt in by
1..-:-,,''''''"Ff7e-.72.14.-,1=.0M-;.1';'.:'-,
i,---7411111111.11c12
Section
x 10 - 4
(D = IOW)
--I 160
800 '(i = 0.8333)
--04 960
p 736 (i 17...7671)
1 9 2 0
1600 kg 6001
1400
.8333 x 1400 = 1167
I.- 2000
(.8333 + .7671) 1000
66001 0 = 1600
_L
1.10001.-10001
2000 -.4
2 do-- 2-3
3 Interior long. 8-9
grid line
4 9-10
7500
210 70mm Wearing coat
-
,1
150
C1-1-1-*--3000
--
1 1"-
1- 600
r100.3 3000
51
A B C D Ci e t A'
0.11950 -1..s.-- 1. @1500: 6000 ----1.-11050 I-
-(a) Cross- section at Mid-span
I
10500
E G
(b) Longitudinal section
Figure 7.9 Right T-Beam Bridge
hogging bending moment over interior support. The span length of exterior
transverse beam will be less than centre to centre distance between adjacent
bearings. It may be assumed as 0.8 times the distance between centres of
adjacent bearings, that is 2.40 m. The flange widths of individual transverse
beams applying the relevant Clause will be then as shown in Fig. 7.10. We
now proceed to idealize the bridge deck into grid structure.
Let us take transverse grid lines first. Referring to Fig. 7.9 (b), adopt grid
lines along centres of transverse beams. These are then at E, G and J. The
effective 'slab widths covered by these transverse T-beams still leave out
large widths of slabs between them. To take into account these remaining
O
O
cV
0
1 MO ,
o
O I
+m, OFR
-
0
on
.0-o
(12
slab strips adopt transverse grid lines between them at F and H giving five
transverse grid lines in half length of the bridge or nine transverse grid lines
in full length of the bridge. The widths of the slab strips represented by these
are shown in Fig. 7.10. The assumed locations of all the nine grid lines
are shown in Fig. 7.11. Now, I and T values to these grid members corre-
sponding to their sections are to be assigned.
In longitudinal direction, referring to Fig. 7.9 (a), grid lines along the
centres of three longitudinal beams are adopted. These are at B, D and B'. I
and J values based on corresponding T-sections of these beams are as-
signed to these grid lines. The spacing of longitudinal beams is 3000 mm
N/RiNsam.,_ MUM
18 36
27 45 63
62
61
a
6 0 0
V
II
in
N
59 _ii
50
eV
co
58
57
56
11
10 55
401500 = WOO I
which is more than spacing of transverse grid lines which is 2625 MILL But
the spacing of longitudinal grid lines should be less than spacing of trans-
verse grid lines for the width to span plan dimension of this bridge. Hence
two additional grid lines at C and C' between the longitudinals are intro-
duced. The live load wheel while traversing over the carriageway may come
on the cantilever slab i.e. on the left of B or right of B' and hence additional
longitudinal grid lines, one each, say at centres of each of the kerb slab at A
and A' be adopted. The inertia value of entire transverse section is covered by
die loligatudinal beams at B, D and B' and hence the grid lines at locations A,
C, C' and A' may be assigned zero inertia values or alternately
Illustrative Examples 191
these grid lines could be assigned small nominal inertia values and corre-
spondingly inertia values of beams at B, D and B' be reduced. The former
option is adopted herein. The location of all seven longitudinal grid lines
is .shown in Fig. 7.11.
The bridge is analysed for dead load and live load. The dead load is due
to self load and is assumed to be acting along the seven longitudinal grid
lines as per corresponding tributary widths. The live load is .due; to OHBD
Truck load and is chosen to demonstrate the use of the program! for a user
specified live loading. The OHBD truck loading is used for the design of
highway bridges in Ontario State of Canada and has been taken from
Ontario Highway Bridge Design Code [2]. The impact factor for the loading
is taken as 0.3. The minimum kerb distance for the loading is taken as 600
mm. The truck loading is given in Fig 7.12. The truck is made to traverse
the length and breadth of the entire carriageway. The lengths of the steps in
longitudinal and transverse directions are kept as 1500 mm and 750 mm.
respectively.
The input data is given in Table 7.9. The longitudinal design bending
moment along exterior and interior beams at one-eighth span point along
6 14 14 20 16 Axle load
3 7 710 8 Wheel load
I
250 250
-
.
4 .0
I
. 0 0
- 0 0
GO 0 st
2SCRE3 U
7
L
00
Wheel 1200
Reference T
)7/2fC2 ____ M
1- *Ti 81 44 -1
3600 6000 7200 11
i
1- 180Q0 .. I
with corresponding positions of the truck is given in Table 7.10. The truck
position is indicated through X-Y coordinates of reference wheel, which is
the left most wheel of the leading axle. The X -Y coordinates are given with
respect to node 1 (Fig. 7.11).
KERB
3@ 5500- 16500
Skew span =2070
(a) Plan
111 1
11
I-
450
li w
e
k
rs
a
cross beam
60
A0C0
H 1475 hi-441250r 475
91
1238
_ 30
78
65 77 in00
52 64 76 89 rc
*4'
z
9 51 63 75 88
26 38 50 62 74 87
13 2750
25 37 61 86
85
. 12 _ 84
24 36 48 6072
11 83 0
2 3 3 5 47 59 71 0
0
10 82
22 34 46 58 7 0 II
0
U)
(--
9 81 csa
21 33 45 576 9 x -
4
8 80
X 0 3 2 44 566 8 1512
7
19 31 43 55 1238
6 6 .79
_____y1
18 30 42
17 29
4
3
2
1 0*
(0,0) 1 .
1--,-4x1250= -`.1 1-.-
1475 5000 1475
200
'150
100
.50
co
0
Bending Moment Envelope for Exterior Beam
7 no
100
oi 50
0
Bending Moment Envelope for Central Beam
30
10
0
Shear
Shear Force Envelope for Interior Beam
0246 8 10 12 14 16 18 20'
Span (m)
Figure 7.16 Bending Moment and Shear Force Envelopes for Beams
overall envelope values for bending moments and shear forces under live
loadings for exterior and interior beams are presented in Fig. 7.16. In pre-
paring these envelope diagrams the response values at polar symmetric
points have been compared and wherever these were different, the larger
one has been adopted. The shear force envelope diagrams show sudden
change in values at locations of cross beams. The envelope diagrams for
shear forces also show maximum reaction values at the supports. The reac-
tions at the two ends of the exterior beam differ considerably. The reaction
is minimum at the acute angled corner and is maximum at the obtuse
angled corner. This is expected in a skew bridge.
4
2
0 625 1000 280 47000 2025
1
0
3
0 1630 1000 385 42000 3940
1
0
5
0 1630 1000 385 29000 3940
1
0
7
0 1845 1000 440 32000 4045
1
0
1
Illustrative Examples 199
3 5 0 0
1. 4176 --oi
(a) Cross Section
7
30 )35
0 4- 225
est
30000
50
450 (c) Longitudinal Section
- 5
slabs beyond B and B'. The moment of inertia of the cross section of the
bridge about a common axis is computed and this is divided equally among
the two longitudinals at B and B'. Similarly torsional inertia of the closed
trapezoidal section is computed using equation 4.15 and one half of this is
assigned to each of the longitudinals at B and B'. Thus it is assumed that in
longitudinal direction, the entire inertia is Concentrated along grid lines at B
and B'. The remaining four longitudinals which are located at A, A',
C and C' each are assigned zero inertia values. Referring to Fig. 7.17 (c,
d), nine transverse grid lines as shown have been assumed. The flexural and
torsional inertia values of these transverse members are computed employing
equation 4.12. The end transverse members also have diaphragms and hence
moment of inertia of diaphragms are added while computing moment of
inertia of end transverse members.
The bridge is analysed for a stationery position of Class 70R Col. 'm'
train loading. The equivalent grid and position of live load on it is shown
in Fig. 7.18. Table 7.12 gives the grid input data. The maximum
longitudinal bending moment at mid-span was found to be 350.0 t.m
giving a bending compressive stress of 1.76 MPa and bending tensile
stress of 3.36 MPa.
To compare the results of aboye grillage analysis, the bridge is alter-
nately analysed by Finite Element Method (FEM). Standard package SAP
IV is used for the analysis. Rectangular plate element with six degrees of
freedom at each node is chosen. Boundary elements are introduced at sup-
port locations for support reactions. Figure 7.19 gives the_ structural
idealisation chosen for the FEM analysis. The maximum bending
compressive and tensile stresses at mid span were found to be 1.78 MPa
and 3.20 MPa respectively. The comparison is quite close and acceptable.
For more detailed information on the comparison, readers may refer to
reference [1].
In the above examples, wherever necessary, the sections adopted for
the grillage .elements are given. The inertia values of these sections can be
taken from the respective grid Input data given for each example.
Y
90
0
93 54 8 2.0 0.87 6 2725.0 737.0 4176.0 737.0 2725.0 450 450
1 2 32 3750. 0. 0. 0.
2 3 3 4 4 5 56 67 7 8 8 9
Contd.
' Illustrative Examples 201
19 20 20 21 21 22 22 23 23 24 24 25 25 26 26 27
28 29 29 30 30 31 31 32 32 33 33 34 34 .35 35 36
46 47 47 48 48 49 49 50 50 51 51 52 52 53 53 54
10 11 16 3750.0 0.0 136.76E+10 203.4E+10
11 12 12 13 13 14 14 15 15 16 16 17 17 18
3.7 38 38 39 39 40 40 41 41 42 42 43 43 44 44 45
t 10 4 0.0 2725.0 0.38E+10 0 .76E+10
37 46 9 18 45 54
10 19 0.0 737.0
4 92.86E+10 185.72E+10
28 37 18 27 36 45
19 28 2 0.0 4176.0 92.86E-F10 185.72E+10
27 36
2 11 14 0.0 2725.0 0.62E+10 1.23E+10
38 47 3 12-39 48 4 13 40 49 5 14 41 50
6 15 42 51 7 16 43 52 8 ,17 4453
11 20 14 0.0 737.0 135.04E+10 270.08E+10
29 38 12 21 30 39 13 22 31 40 14 23 32 41
15 24 33 42 16 25 34 43 17 26 35 44
20 29 7 0.0 4176.0 135.04E+10 270.08E+10
21 30 22 31 23 32 24 33 25 34 26 35
4 3 19 27 28 36
1800.
1
0. 0. 0. 0.0. 0. 0. 0. 0. 0. 0. 0.
0
I
1
0
0
1
3
23980. 3205. 1. 1. 23980. 3205.
1
0
S
202 Grillage Analogy in Bridge Deck Analysis
795 790
f3205 795 e
-0
positionsof
7
2
9 3
6 w
h
els
1 I I - 53
11101111 43 2
St
12t
51
8@3750m.30000
12t
F
42
MINI
14
cc;
171
17t -
Yir
4 17
ao
40 17t
N
3
30
121
12 39
17
212
0
1075
9
11725.4 I.-3500 --I 27251
1075
. c
H
0 In on 0 71:7 0 IA 0
01 0 to
.111 -gr-
ati
te r co 0 -t
O
U
1 In
5
O
g
2
0
0
0
3
0
0
0
0
I0 O
w %.t w tp O c-I
-1
1 2 725 1180188411E183 2725 I 1917 11917) 184.1. ["1--
1144 1146
REFERENCES
1. ALOK 13uowt.ucK, "To study the Behaviour of Box-Girder Bridges under Live
Load by Analysing the Bridge Using Several Methods and Comparing the Re-
sults", An M.Tech. Project, Civil Engineering Department, I.I.T. Delhi. 1990.
2. IRC 21-1987, "Standard Specifications and Code of Practice for Road Bridges",
Section III, The Indian Roads Congress, New Delhi. 1991.
3. "Ontario Highway Bridge Design Code", Ministry of Transportation and Com-
munications, Highway Engineering Design. Downsview. Ontario, Canada, 1983.
9
r t
Appendix I
Listing of Programs
Grid'*
C PROGRAM 'GRID'
C THIS PROGRAM CAN ANALYZE A GRID UNDER APPLIED NODAL
C LOADS. THE CONFIGURATION OF THE GRID COULD BE RIGHT,
C SKEW OR QUADILATERAL.
IMPLICIT DOUBLEPRECISION (A-H2O-Z)
COM MON/BX I /X(100), Y(I 00), XI(I00), XJ(100), SM (6,6),
NCD (6), NN(50) COMMON/BX2/GSM(200,50)
COMMON/BX3/NCN1(200),NCN2(200),NSN(20),SUPSTF(20)
COMMON/BX4/P(300),D(300)
COMMON/BX5/GDM(6),DM(6),PM(6),TM(6,6)
COMMON/BX6/LNODE(200),NLT(200),XLOAD(200)
OPEN(UN1T=3,FILE='GRID.IN')
OPEN(UNIT=4,FILE='GRID.OUT)
C COORDINATES OF NODES
READ (3,*)NGROUP
DO 20 I=I, NGROUP
READ (3 ,*)N1 ,N2,NINC, AX,AY,XINC,YINC
DO 15 11=N1,N2,NINC
X(I1)7AX
Y(I1) =AY
AX=AX+XINC
AY =AY+Y1NC
15 CONTINUE
20 CONTINUE
C MEMBER CONNECTIVITIES
READ (3,*)NGROUP
DO 30 1=1, NGROUP
READ(3,*)N1,N2,N1NC,NC I ,NC2,NC11NC, NC21NC
DO 25 11=NI,N2,NINC
NCN1(I1)=NCI
NCN2(I1)=NC2
NC 1 =NCI +NC IINC
NC2=NC2+NC2INC
25 CONTINUE
30 CONTINUE
C MEMBER PROPERTIES
READ (3,*)E,G
READ (3,*)NGROUP
Appendix I 207
DO 40 I=1,NGROUP
READ(_S,*)ALAJ,N,(NN(I1),I1 =1,N)
DO 35 I1=1,N
12=NN(11)
XI(12)=AI
XJ(12)=AJ
35 CO N TI NU E
4 0 C O N TI N U E
SM(3,2)=-X3*C*S+X4*C*S
SM(3,3)=X3*C*C-I-X4*S*S
SM(4,1)=-X1
SM(4,2)=-X2*S
SM(4,3)=X2*C
SM(4,4)=X1
SM(5,1)=X2*S
SM(5,2)=X3/2.*S*S-X4*C*C
SM (5 ,3) =-X3/2.*C*S-X4*C*S
SM(5,4)=-X2*S
SM(5,5)=X3*S*S-i-X4*C*C
SM(6,1)=:-X2*C
SM(6,2)=-X3/2.*C*S-X4*C*S
SM(6,3)=X.3/2.*C*C-X4*S*S
SM(6,4) =X2*C
SM(6,5)=-X3*C*S+X4*C*S
SM(6,6)---X3*C*C+X4*S*S
DO 80 J=1,5
DO 80 K=(.I+1),6
80 SM(J,K)=SM(K,J)
C 'SUPPORT CONDITION
C READ SUPPORT TYPES- 0 FOR FLEXIBLE, 1 FOR RIGID
READ (3,*)NST
IF(NST.EQ.0) GC) TO 100
READ (3,*)NSUP
READ(3,*)(NSN(I),I = 1 ,NSUP)
DO 95 I=1,NSUP
IDOF = 3*(NSN(I)-1)+ 1
GSM(IDOF,1)=GSM(IDOF, 1)*1.0E + 6
95 CONTINUE
GO TO 110
100 READ (3,*)NSUP
READ(3,*) (NSN(1),SUPSTF(0,1=1,NSUP)
DO 105 I=1,NSUP
IDOF=3*(I-1)+1
ncivr(1nrw i)=GSM(1DOF,I)+SUPSTF(I)
105 CONTINUE
Appendix 1 209
C LOAD VECTOR
. READ (3,1)NLOADS
DO 120 I=1,NLOADS
READ(3,*)LNODE(I),NLT(I),XLOAD(I)
. 120 CONTINUE
DO 125 I=1,NLOADS
IDOF=3*(LNODE(1)-1)+NLT(1)
P(IDOF) =XLOAD(I)
125 CONTINUE
C MEMBER FORCES
WRITE(4,*)
WRITE(4,I'MEMBER FORCES'
WRITE(4;145)
145 FORMAT(1X,'MEM.NO.',IX,'S.FORCE',7X,'B.MOMENT',12X,'TORSION')
WRITE(4,150)
150 FORMAT(26X,'END 1',9X,'END 2')
DO 195 I= LNELEMS
TM(3,2)=-S
1:
TM(3,3)=C DO
160 IA=4,6 II 1.
=1A-3
DO 160 IB=-4,6
1
11=1B-3
160 TM(IA,113)=TM(11,1.1) 1
C SUBROUTINE DECOMP:
C DECOMPOSITION.OF STIFFNESS MATRIX SUING
C CHOLESICEY'S METHOD.
C N ---->TOTAL NO. OF EQUATIONS
C IBW = >HALF BAND wrai
OF THE MATRIX
C EXIT = >INDICATING VARIABLE
C S = <STIFFNESS MATRIX TO BE DECOMPOSED
SUBROUTINE DECOMP(N,IBW,EXIT)
IMPLICIT DOUBLEPRECISION (A-H2O-Z)
COMMON/BX2/S(200,50)
WRITE(*,*) DECOMPOSITION STARTS'
DO 260 I=1,N.
WR1TE(s,*)I
IP=N-I+ 1
IF (IBW.LE.IP) GOTO 200
GOTO 205
200 IP=IBW
205 DO 260 LI P
SUM =S(I,J)
IQ=IBW-J
IF ((I-1).LE.IQ)GOTO 210
GOTO 215
210 IQ =I-1
215 IF (IQ.EQ.0) GOTO 230
GOTO 220
220 DO 225 K =1,IQ
225 SUM =SUM-S(I-K,1+K)*S(I-K,J+K)
230 IF (I.NE.1) GOTO 235
GOTO 240
235 S(I,J) =SUM/TEMP
GOTO 260 -
C37 IF (SUM.LE.0.0) GOTO 245
240 GOTO 250
245 EXIT=0.0
RETURN
250 EXIT =1.0
255 TEMP=DSQRT(SUM)
S(I,J) =TEMP
260 CONTINUE
RETURN
END
C SUBROUTINE SOLD':
C SOLUTION OF STIFFNESS MATRIX DECOMPOSED BY
212 Grillage Analogy 0: Bridge Deck Analysis
C CHOLESKEY 'S METHOD FOR GIVEN LOAD 5
D7 = -Pl*Pl*Pa- Fr2*P2*P3
D8 =-P1*P2*(PS -P4) D9=
P I *P 1 *P3 -In *P2*F4
M1=3*M(I,1)
M2 =3*M(I,2)
M3 =M1-1
M4 =M1-2
M5=M2-1
M6 =M2-2
MXI = D l*B(Mrh) +D2*B(M3) + D3*(B(M1)-B(M2)) + D7*B(M6) + D8*B(M5)
MY1=D2*B(146.4)+D4*B(M3)+D5*(B(M1)-B(M2))+D8*B(N16)+D9*B(M5)
MX2 =D7*B(M4) + D8*B(M3) + D3*(B(M1)-B(M2)) +D1*B(1V16) +D2*B(M5)
MY2 =D8*B(145.4) +D9*B(M3)+D5*(B(M1)-B(M2))+D2*B(M6)+D4*B(M5)
FZ =D3*(B(M--4)+B(M6))+D5*(B(M3)+B(M5))+D6*(B(M1)-B(M2))
FZ1 =F5
J=I-IBW+ I
IF ((I +1).LE.JBW) GOTO 265
GOTO 265
265 J=1
DO 270 K=1,1-1
270 SUM = SUM-S(K, (I-K+ 1))*D(K)
275 D (I) = SU/s4/S(I ,I)
280 CONTINUE DO
305 I =- I ,N
J=N-I+IBW IF (J.
GT. GOTO 285
GOTO 290
285 J=N
290 SUM =D(N-1+ I)
IF (I, EQ.1) GOTO 300
DO 295 K=(N-I+2),1
'295 SUM =SUM-S(N-I+1,K-N +D*D(K)
300 D(N-I+1)=SUM/S(N-I+1,1)
305 CONTINUE
RETURN
END
Appendix II
Listing of Program
Gabs*
PROGRAM 'GABS'
C THIS PROGRAM CAN ANALYSE A GRID UNDER GENERALISED DECK
C LOADING.THE CONFIGURATION OF THE GRID COULD BE RIGHT ,
C SKEW OR TRIANGULAR.
IMPLICIT real*8"(A-H2O-Z)
INTEGER GG,DD,WW,R,T,EE,EE1,EE2,T1,R3,PATYPE,GRIDTYPE
real*8 L,MP1:,MP2,MQ1,MQ2,MX1,MX2,MY1,MY2,E1,E2,
+MAXPT(6,4;500),MAXB(4,500),MAXRT(4, 140),MAXPT2(6,2:4,500),
+MAXB2(2:4;500),MAXRT2(2:4,100),MAXRTN(4,100),CORSPT(2,500),
+MAXR2N(2:4,100),CORSPTI(2,500),CORSPT2(2,500),CORK(100),
+CORX1(100),CORY1(100),CORIX(100),COR1Y(160),COR2X(100),
+COR1X1(100),COR1Y1(100),COR2X1(100),COR2Y1(100),
+ coRRx(2,500), CORRY(2 ,500),CORRX1 (2,500), CORRX2(2,500),
+ CORY(100);COR2Y(100), CORRY1(2,500),CORRY2(2,500),
.+CORD(100),CORD1(100),CORD2(100)
CHARACTER*1 CH1,CH3
CHARACTER*2 POS, INPFILE*12, OUTFILE*12
COMMON/BLOCK2/S(1,50000),U(50000),P(1000,7),PT(1000,7),W(100),
+D(100),B(1000),BD(1000),RT(100),FEL(500,6),RTN(100)
COMMON/BLO CK3/XC (100), YC(100),P W(100),PI ,SANG,PT1 (1000 ,2),
+DLB(1000),XLLG(0:100),YLLG(0:100) -
DIMENSION M(1000,2), ITT(100),JDD(100),JDD1(100),
+PM(1000,1)
DIMENSION DLNG(100),BL(1000)
COMMON/BLOCK6/ SCLMIN,XLL,YLL,AKERB
COMMON/BLOCK8/0MAXPT(6,4,1000), oMAX13(4,1000), OMAXRT(4,100),
+0MAXRN(4,100)
COMMON/XX1/AIMP
CHARACTER*20 FILE1, FILE2
WRITE(*,*) 'Enter the file name for Input '
READ(*,5) FILEI
5 FORMAT(A20)
WRITE(*,*) 'Enter the file name for Output
READ(*,5) FILE2
OPEN(UNIT=3,FILE:=F1LE1)
OPEN(UNIT =4 ,FILE =FILE2)
C OPEN(UNIT=3,FILE='*SRC.EX2')
C OPEN(UNITL-4,FILE='*SRC.GRID3OU')
MOVE=O
XL1 =0.0
YLI =0.0
XINT=0.0
YINT=0.0
DO 10 i=1.6
DO 10 J=1,4
Appendix II 215
DO 10 K=1,500
10 MAXPT(I,J,K)=0.0
DO 15 1=1,4
DO 15 3=1,500
15 MAXB(I,J)=0.0
DO 20 I=1,4
DO 20 3=1,100
MAXRTN(I,J)=0
20 MAXRT(I,J)=0
P1=4*ATAN(1.0)
25 READ(3;230) CHI
IF (CHI.EQ.'N') GO TO 315
READ(3,*) GRIDTYPE
C GRIDTYPE CAN EITHER BE SKEW(ANY VALUE OTHER. THAN 90) OR
C ORT H.(90)
READ (3,*) ANG
READ (3,*) MM,N,DD,E,G,NLG,(DLNG(I),I=1,NLG-1)
READ (3,*) E1,E2
SANG = ANG*PI/180
WRITE (4,235)
WRITE (4,240) ANG
WRITE (4,245) MM,N,DD,E,G
WRITE(4,290)NLG,(DLNG(I),I=1,NLG-1)
WRITE(4,295)
WW =0
1=0
DO 45 GG=1, DD
), I= I+ 1
READ (3,*) M(I,1),M (1,2),T,P(I,1),P(1,2),P(1,3),P(I,4)
WRITE (4,250) GG
WRITE (4,255)T,(P(1,1C2),K2=1,4)
PM(I, 1) = P(I , 1).
L=SQRT(P(I,1)**2+P(I,2)**2)
P(I,1) =P(I,
p(1,2)=P(1,2)/L.
P(I,3)=2*E*P(1,3)/L
H =P(1,3)
P(I,4)=G*P(I,4)/L
P(I,5) = 6*H/(L**2)
P(1,6)=2*H
P(I,7) =3 *HiL
EE=I
EE1=EE+1
EE2=EE+T-1
DO 30 II I =EE1,EE2
30 PM(II I, I) = PM(EE,1)
READ (3,*) (M(IAX, 1), M(IAX ,2),IAX =EE1,EE2) Nil
1+T 1
216 Grillage Analogy in Bridge Deck Analysis
Appendix II 217
CALL SORT(IX,JDD,JDD1)
CALL MODSTIF(WW,IX,JDD1,S,NROWS)
ELSE
DO 70 1=1,1X
NE=WW*(JDD(I)-1)+ I
S(1,NE)=S(1,NE)-1-SUFK
70 CONTINUE
ENDIF
CALL DECOMP(NROWS,WW,S,U)
C COMPUTATIONS OF LOAD MATRIX STARTS
75 READ (3,230) CH3
IF (CH3.EQ.'S') GO TO 315
NL=N
N9=3*N+WW
READ (3,*) 'NOLCASE
WRITE(4,80)
80 FORMAWING(1) = DEAD LOAD ONLY'/
+ ' LDING(2) = CLASS A - TWO LANE'/
+ ' LDING(3) = CLASS 70-R TRAIN:COL -L'/
+ ' LDING(4) = CLASS 70-R TRAIN:COL -M'/
LDING(5) = CLASS 70-R BOGIE:COL -1.1/
LDING(6) = CLASS 70-R-BOGIE:COL -M'/
LDING(7) = CLASS 70-R TRACK'/
+ ' LDING(8) = CLASS A - SINGLE LANE'/
+ ' LDING(9) = SPECIFIED BY USER '/)
DO 225 IMM1=1,NOLCASE
READ(3,*)LCASE
IF(LCASE.EQ.I)THEN
WRITE(*,*)'DEAD LOAD ONLY'
WRITE(4,*)'DEAD LOAD RESULTS'
DO 85 JA3=1,3*NL
85 B(JA3)=DLB(JA3)
GOTO 135
ENDIF
READ(3,*)XL,YL,XINC,YINC,XSTOP,YSTOP
MODIN1=0
MODIN2=0
XFIN = XL + REAL(INT((XSTOP-XL)/XINC))*XINC
YFIN = YL REAL(INT((YSTOP-YL)/YINC))*YINC
WRITE(4:*)'
WRITE(4,310)LCASE
WRITE(4,*)' ------------------------
CALL IMPACT (LCASE, SPAN,SANG)
WRITE(4,'(10X,A)TTHE RESPONSE IS INCLUSIVE OF DEAD LOAD RESULTS'
WRITE (4,300)XL, YL, XINC, YINC,XFIN, YFIN,AIMP
MOVE =0
IM1 =0
YL1 =YL-YINC
218 :Grillage Analogy in Bridge Deck Analysis
90 XL1 =XL-XINC
YL1 =YL1+YINC
95 XL1 =XLI + XINC
IF (YLLGT.YSTOP)THEN
GO TO 200
ENDIF
IF (XLLGT.XSTOP)THEN
GO TO 90
ENDIF
IF(MODIN1. EQ.1) THEN
XL1 =XL1-1.0
MODIN1 =0
ENDIF
IF(XLI.EQ.0) THEN
XL1=1.0
MODINI =1
ENDIF
IF(MODIN2.EQ.1) THEN
XL1 =XL1 +1.0
MODIN2=0
ENDIF
IF(XLLEQ.SPAN) THEN
XL1=XL1-1.0
MODIN2=.1
ENDIF
C WRITE(*,*)'LOAD CASE ',LCASE,' LEFT FRONT WHEEL AT ',XL1,
C +' AND ',YL1
IM1=IMI+1
CALL LLOAD (LCASE, XLI, YL1, XC, YC, NW; NWPA,IM1)
CALL WLOAD (LCASE, PW,SCLMIN, NW, IM1)
IXL=0
IXG=0
DO 100 JA2=1,3*NL
100 B(JA2)=0.0
DO 110 JA2=1,MM
DO 105 JA3=1,6
105 FEL(IA2,1A3)=0.0
110 CONTINUE
DO 120 I5P=1,NW
IF ((GRIDTYPE.NE.90).bR.(SANG.EQ.0.0)) THEN
CALL PAIDRB(XLLG,YLLG,XC(15P),YC(I5P),SANG,IXL,IXG,A10ERB
, SCLMIN, SPAN, POS, PATYPE, MN1, MN 2, MN11 , MN22, X LL, YLL, XPL, YPL
+ ,NTG,NLG)
ELSE
CALL PAID1(M,PM,N,NLG,DLNG,MN1,MN2,MN11.MN22,PATYPE,
+ XC(I5P),YC(15P),XLL,YLL,XPL,YPL,SANG,IXG,IXL,POS,
+ AKERB,SCLMIN)
ENDIF
Appendix II 219
IF(POS.EQ.'YOGO TO 90
IF(POS.EQ.'YG')GO TO 90
IF(PATYPE.EQ.0)G0 TO 120
IF(PATYPE.EQ.3)THEN
CALL LDISTT (MN1,MN2,MN22,PW(15P),XPL,YPL,XLL,YLL,BL,POS,NL)
ELSEIF(PATYPE.EQ.4)THEN
IF (GRIDTYPE.EQ.90)THEN
SANG! =0.0
CALL LDISTR (MN1,MN2,MN1I,MN22,PW(I5P),SANG1,13L,XPL,YPL,XLL,
+ Yi,L,POS,NL)
ELSE
SANGI=SANG
CALL LDISTR (MN1,MN2,MN11,
+ MN22,PW(15P),SANG1,BL,XPL,YPL,XLL,YLL,POS,NL)
ENDIF
ENDIF
DO 115 JOG=I,3*NL
115 B(JOG)=B(JOG)+BL(JOG)
120 CONTINUE
IF(IXG.EQ.NW)GO TO 90
IF(LXL.EQ.NW)G0 TO 95
MOVE=MOVE+1
DO 125 IFI=1,3*N
125 B(IFI)=B(IFI)*AIMP
DO 130 JA3=1,3*NL
130 B(JA3)=B(JA3)+DLB(IA3)
C COMPUTATION OF LOAD MATRIX ENDS
135 DO 140 I=1,NROWS1
BD(I)=-B(1)
140 CONTINUE
IF (JRR.EQ.3) THEN
DO 150 1=1X,1,-1
J=JDDI(I)
DO 145 IP= J,NROWS1
B(IP)=B(IP+1)
145 CONTINUE
150 CONTINUE
CALL SOLVE (NROWS,WW,U,B,B)
DO 155 I=NROWS+1,NROWSI
155 B(I)=0
DO 165 I=1,IX
J=.113D1 (I)
DO 160 IP =NROWSI,J,-1
B(IP + I) =B(IP)
160 CONTINUE
B(IP +1)=0
165 CONTINUE
LT C."0
220 Grillage Analogy in Bridge Deck Analysis
CALL SOLVE(NROW S1 , WW , U, B , B)
ENDIF
C COMPUTATIONS OF ELEMENT FORCES STARTS
DO 170 I =1,MM
P1 =P(I,1)
P2 =P(I,2)
P3 =P(I ,3)
P4 =P(1,4)
P5 =P(1,5)
P6 = P(I, 6)
P7 =P(I,7)
D1 =-P1*P1*P4 +P2*P2*P6
D2 =P1*P2*(P4-P6)
D3 = P2*P7
D4 =P2*P2*P4+P 1 *Pl*P6
D5 = -P 1 *P7
D6 = P5
D7 = -Pl*P 1*P4 +P2*P2*P3
D8 =-Pl*P2*(P3 +P4)
D9 =P1*P1*133-P2*P2*P4
M =3*M(I, I)
M2 =3*M(I,2)
M3 =M1-1
M4 =M1-2
M5 =M2-1
M6 = M2-2
MX1 =D1*B(M4)+D2*B(M3)+D3*(B(M1)-B(M2))+D7*B(M6)+D8*B(M5)
MY1 = D2*B(M4) + D4*B(M3) + D5*(B (M1)-B (M2)) + D8*B(M6) + D9*B (M5)
MX2 = D7493(M4) + D 8*B (M3) + D3*(13(M1)-B(M2)) + D 1*B (M6) + D2*B(M5)
MY2 = D8*B (M4) + D9*13(M3) + D5*(B (M1)-B(M2)) + D2*B(M6) + D4*B (M5)
FZ = D3*(B (M4) + B(M6)) + D5*(B (M3) + B (M5 ))+136*(B(M1)-B (M2))
FZ1= FZ
MQ 1= (-MX1)*P2 +(MY1)*P1
MQ2 = (-MX2)*P2 + (MY2)*P1
MP1 = (MX1)*P1 +(MY1)*P2
MP2 = (MX2)*P1 +(MY2)*P2
PT(I,3) =FZ1
PT(I,4) = MP 1
PT(L5) = MQ I
PT(I, 6) =MQ2
170 CONTINUE
C COMPUTATIONS OF ELEMENT FORCES ENDS
DO 175 II =1,1X
RT(11 ) =0
175 CONTINUE
DO 185 11=1,1X
DO 180 12=1,MM
IFWTT(11 ) . EQ. Ivi(12, 1)).OR.(11-f(11).EQ.M(12,2))) THEN
Appendix 11 221
IF (JTT(II).EQ.M(I2,I)) THEN
RT(I1) =RT(I1) +PT(12,3)
ELSE
RT(I1) =RT(I1)-PT(I2,3)
ENDIF
ENDIF
180 CONTINUE
185 CONTINUE
DO 190 II =1,IX
NP1=3*JTT(II)
RT(I1)=RT(11)-BD(NP1)
190 CONTINUE
DO 195 JP8=1,2
DO 195 IP8=1,MM
KP8 =JP8 +4
195 PTI(IP8,11)8)=PT(IP8,KP8)
CALL COMPARE1 (PT,B,RT,MAXPT,MAXB,MAXRT,MAXRTN,MM,N,IX,
+MOVE,PT1,XL1,YLI,LCASE,CORSPT,CORX,CORY,CORXI,CORY1,
+ CORRX, CORRY, CORD I ,JTT)
IF(LCASE.EQ.1) GOTO 200
GO TO 95
200 READ (3,*)IIC2
IF (IK2.EQ.I)CALL WRITE1(B,MAX7F,MAXB,MAXRT,MAXRTN,MM,
+N,IX,MOVE,M,JTT,LCASE,CORSPT,CORX,CORY,CORX1,CORY1,
+CORRX,CORRY,CORDI)
IF(NOLCASE.EQ.1) GO TO 75
DO 205 IL1=1,4
IL2 =11,1+2
DO 205 IL3 =1,MM
PT (IL3,112)=MAXPT(IL1,1,IL3)
DO 205 1=2,4
MAXPT2(IL I ,I ,IL3) =MAXPT(IL1,I,IL3)
CORSPT2(1,IL3)= CORSPT(1,IL3)
CORSPT2(2,IL3) =CORSPT(2,IL3)
205 CONTINUE
DO 210 1=1,N
DO 210 .1=3*1-2,3*1
B(J)=MAX13(1,J)
DO 210 K=2,4
MAXBAK,J)=MAXB(K,J)
CORRX2(1,I)=CORRX(1,I)
CORRY2(1,I)=CORRY(1,I)
CORRX2(2,I)=CORRX(2,I)
CORRY2(2,I) =CORRY(2,I)
210 CONTINUE
DO 215 1=1,1X
RT(I) =MAXRT(I ,I)
RTN(I)=MAXICIN ,
222 Grillage Analogy in Bridge Deck Analysis
DO 215 J=2,4
MAXRT2(J,I) =MAXRT(J,1) 2'
MAXR2N(.1,1)=MAXRTN(J,I) 2
COR2X(D=CORX(I) 2
COR2Y(I)=CORY(I) 2
COR2X1(I)= CORX I (I)
COR2Y1 (I) = CORY1 (I)
CORD2(I)=CORD1(I)
215 CONTINUE
DO 220 JP8=1,2
KP8=JP8+4
DO 220 108=1,MM
PTI(IP8,JP8) = MAXPT(ICP8, 1, IP8)
DO 220 1=2,4
MAXPT2(KP8,I,IP8)=MAXPT(ICP8,LIP8)
220 CONTINUE
CALL COMPARE2(PT,B,RT,RTN,OMAXPT,OMAXB,OMAXRT,OMAXRN,MM,
+N,IX,IMM1,PT1,MAXPT2,
+MAXB2,MAXRT2,MAXR2N,JTT,COR2X,COR2Y,COR2X1,COR2Y1,
+COR1X,CORIY,COR1X1,CORIYI,CORSPT1,CORSPT2,
+CORRX1,CORRX2,CORRY1,CORRY2,CORD,CORD2)
225 CONTINUE
WRITE(4,80)
LCASE =0
READ(3,*) 11C3
IF (IK3.NE.1) GO TO 75
CALL WRITE2(0MAXPT,OMAX13,0MAXRT,OMAXRN,MM,N,IX,MOVE,M,JTT,
+ LCASE,CORIX,CORlY,CORIX1 , CORlY I ,CORSPT1,CORRX1,CORRYI,CORD)
GO TO 75
230 FORMAT (1A1)
235 FORMAT(2(/),30X,'BRIDGE DETAILS'/29X,16('*')//)
240 FORMAT(10X,'SKEW ANGLE = ',F9.3, ' DEGREES')
245 FORMAT(10X,'NO. OF ELEMENTS= ',13/,10X, 'NO. OF NODES=',
+ 13/,10X,'NO. OF ELEMENT GR.OUPS-=',12/,
+ 10X, 'YOUNGS MODULUS OF ELASTICITY=',
+ E14.7, ' T/SQ.MMVIOX,'SHEAR MODULUS OF ELASTICITY = ',E14.7,
+ ' T/SQ.MMV)
250 FORMAT(//15X,'GROUP',I2/15X,7('-'))
255 FORMAT(20X,'NO. OF ELEMENTS IN THIS GROUP = ',12/20X,
+ 'LENGTH OF ELEMENT IN X-DIRECTION =',F7.0,' MM'/20X,
+ 'LENGTH OF ELEMENT IN Y-DIRECTION =',F7.0,' MM'/20X,
+ 'MOMENT OF INERTIA OF ELEMENT =',E14.7,' MM4'/20X,
+ 'TORSIONAL INERTIA OF ELEMENT =', E14.7,' MM4') 260
FORMAT(/20X,'ELEMENTS IN THIS GROUP:'//7(20X,5(12,'-',
+ I2,3X)/))
265 FORMAT (/15X,'BOUNDARY CONDITIONS-SUPPORT NODES'
+ /15X, 33 ('-'))
270 FORMAT (20X,'NO. OF SUPPORT NODES 12120X,
Appendix Il 223
S(1,M3+1)=S(1,M3+1)+Pl*Pl*P4+P2*P2*P6
S(1 ,M3 +2) =S(1,M3 + 2) +P I *P2*(P4-P6)
S(1,M3 +3) =S( I ,M3 +3)-P2*P7
S(1,M2+1)=S(I,M2+1)+P2*P2*P4+Pl*P1*P6
S(1,M2+2)=S(1,M2+2)+P1*P7
S(I,M1+1)=S(1,M1+1)+P5
C ABOVE STATMENTS SET UP DIAGONAL ELEMENTS
U=3*IABS(M(I,2)-M(I,1))
R=3*M(I,1)
RI =(R-1)*'WW+U
R2 =(R-2)*W1V+U
R3 =(R-3)*WW +U
S(1,R3+1)=S(1,R3+1)-Pl*P1*P4+P2*P2*P3
S(1,R3+2)=-S(1,R3+2)-Pl*P2*(P4+P3) S(1
,R3 +3) = S(1,R3 +3)-P2*P7
S(1,R2)=S(1,R2)-P1*P2*(P4+P3)
S(1,R2 + I) = S(1,R2+ 1)-P2*P2*P4 +Pl*Pl*P3
S(1,R2+2)=S(1,R2+2)+P1*P7
S(1,R1-1)=S(I,R1-1)+P2*P7
S(1,R1)=S(1,R1)-PI*P7
S(1,R1+ I)=S(1,R1+1)-P5
325 CONTINUE
C ABOVE STATMENT SET UP OFF DIAGONAL ELEMENTS
RETURN
END
c***********************************************************************
SUBROUTINE INODE (M,I1,I2,16)
IMPLICIT real*8(A-H2O-Z)
DIMENSION M(1000,2)
IF (12.LT.I1) THEN
WRITE(*,*)*ERROR'
READ(*,*)0
ELSE
K1=0
330 K1 =K1+1
IF (M(K1,1).NE.II.OR.M(K1,2).NE.I2) GO TO 330
I6=K1
ENDIF
RETURN
END
c***********************************************************************
SUBROUTINE PAID1(M,PM,N,NLG,DLNG,MN1,Mi',12,MN11,MN22,
+ PT,XL,YL,XLL,YLL,XPL,YPL,SANG,IXG,IXL,POS,AKERB,SCLMIN)
IMPLICIT real*8(A-H2O-Z)
INTEGER PT
CHARACTER*2 POS
DIMENSION M(1000,2),PM(1000,1),DLNG(100),XLLG(100),SLNG(I00)
POS='NN'
Appendix-11 225
X=0
DO 335 I=1,NLG-1
SLNG(I)=0
X=X+DLNG(I)
335 SLNG(1)=X -
SPAN =0
DO 340 !IFF =1 ,(N/NLG-1)
XLLG(IFF) =0
CALL INODE(M,IFF,IFF+1,IN1)
SPAN=SPAN+PM(IN1,1)
XLLG(IFF)=SPAN
340 CONTINUE
IF (YL.LT.(AKERB+SCLMIN))GO TO 470
IF(YL.GT.(SLNG(NLG-I)-(AKERE+SCLMIN)))GO TO 475
IF((XL.LT.(SLNG(NLG-I)*TAN(SANG))).AND.(YL.GT.XL/
+ TAN(SANG))) GO TO 480
IF ((XL.GT.SPAN).AND.(YL.LT.((XL-SPAN)/TAN(SANG))))
+ GO TO 485
XLEN=0.0
DO 345 II =1,NLG-1
CALL INODE(M,I1,II + I ,LI)
345 XLEN=XLEN+PM(L1,1)
IF ((XL.GT.XLEN).AND.(XL.LE.SPAN))GO TO 350
IF (XL.LT.XLEN) GO TO 380
IF (XL.GT.SPAN) GO TO 425
3 5 0 PT = 4
DO 355 II =NLG,N/NLG-I
CALL INODE (M,I1J1+1,L1)
XLEN=XLEN+PM(L1,1)
IF (XL.LE.XLEN) GO TO 360
355 CONTINUE
360 IF(YL.GT.0)GOTO 365
MN 1=11
MN 2=11+1
MN 11 = II +(N/NLG-1)
MN 22=11+N/NLG
XLL=XL-(XLEN-PM(LI,1))
YLL=0.0
XPL=PM(L1,1)
YPL=DLNG(1)
RETURN
365 DO 370 12=1,NLG-1
IF ((SLNG(12)-YL).GE.0.0)G0 TO 375
370 CONTINUE
375 MN 1=11+ (I2-1)*(N/NLG- )
MN 2=iviNi+1
MN 1I=MN1+(N/NLG-1)
MN 22=MN11+1
226 Grillage Analogy in Bridge Deck Analysis
XLL=XL-(XLEN-PM(L1,1))
IF(12.NE.1)YLL=YL-SLNG(I2-1)
IF(12.EQ.1)YLL=YL
XPL=PM(L1,1)
YPL=DLNG(12)
RETURN
380 POS = 'XL'
DO 390 13= 1,NLG-1
IF(13.NE.1)GOTO 385
IFaL.LE.XLLG(1)).AND.(YL.LE.XL/TAN(SANG))) GOTO 395
GOTO 390
385 IF ((XL.LE.XLLG(I3)).AND.(YL.LE.XL/TAN(SANG)).AND.(YL.GE.SLNG
+ (13-1))) GO TO 395
390 CONTINUE
GO TO 400
395 PT=3
MN 1 =1 ,+(I3-1)*N/NLG
MN 2=MN1+I
MN 22=MN1+N/NLG
MN 11=0
CALL INODE(M,I3,13 +1,131)
XPL=PM(I3I,1)
YPL=DLNG(13)
IF(13.NE.1)XLL=XL-XLLG(I3-1)
IF(I3.NE. DYLL = YL,SLNG(13-1)
IF(I3.EQ.1)XLL=XL
IF(I3.EQ.1)YLL=YL
RETURN
400 PT =4
DO 405 14 =2,NLG-1
IF (XL.LE.XLLG(I4)) GO TO 410
405 CONTINUE
410 DO 415 15=1,NLG-1
IF ((YL-SLNG(15)).LE.0) GO TO 420
415 CONTINUE
420 MN I =14+(15-1)*(N/NLG-1)
MN 2=MN1+1
MN 11 =MN1 +N/NLG- I
MN 22 = MNI 1 + 1
CALL INODE(M,I4,14+1,I3I)
XPL=PM(I3I,I)
YPL=DLNG(I5)
XLL = X L-XL.TL-G(I4-1)
IF(15.NE.1)YLL =YL-SLNG(15-1)
IF(I5.EQ.1)YLL=YL
RETURN
425 POS = 'XG'
XXL= (SPAN +XLLG(NLG-1))-XL
Appendix II 227
YYL=SLNG(NLG-1)-YL
DO 435 16=1,NLG-1
IF(I6.NE.1)GOTO 430
IF((XXL.LEALLG(1)).AND.(YYLLE.XXL/TAN(SANG)))GOTO 440
GOTO 435
430 IF ((XXL.LE.XLLG(16)).AND.(YYLLEJO(L/TAN(SANG)).
+ AND.(YYLGE.SLNG(16-1))) GO TO 440
435 CONTINUE
GO TO 445
440 PT =3
MN I =(146)*N/NLG+N
MN 2=MN1-1
MN 22=MN1-NINLG
MN 11=0
CALL INODE(M,I6,16+1,131)
XPL = PM(I3I, 1)
YPL=DLNG(NLG-16)
IF(I6.NE.1)XL,L=XXL-XLLG(16-1)
IF(I6 .NE. DYLL =YYL-SLNG(I6-1)
IE(16.EQ.1)XLL=XXL
IF(I6. EQ. 1)YLL=YYL
RETURN
445 PT=4
DO 450 17=-2,NLG-1
IF (XXL.LEILLG(17)) GO TO 455
450 CONTINUE
455 DO 460 I8=1,NLG-1
IF ((YYL-SLNG(18)).LE.0) GO TO 465
460 CONTINUE
465 MN22= N-(18-1)*N/NLG-(I7-18)
MN11=MN22-1
MN2=MN22-(N/NLG-1)
MN1=MN2-1
CALL INODE(M,I7,17+1,131)
XPL=PM(I3I,1)
YPL=DLNG(I8)
XLL=XLLG(I7)-XXL
IF(18.NE.1)YLL =SLNG(18)-YYL
RETURN
470 POS='YL'
RETURN
475 POS='YG'
RETURN
480 IXL =Del,4-1
GO TO 490
485 IXG =IXG+ I
490 PT =0
RETURN
228 Grillage Analogy in Bridge Deck Analysis
END
c**********************************************************************
SUBROUTINE PAIDRB(XLLG,YLLG,XL,YL,SANG,IXL,IXG,AKERB,SCLMIN,
SPAN,POS,PT,MN1,MN2,MN11,MN22,XLL,YLL,XPL,YPL,NTG,NLG)
IMPLICIT real*8(A-H2O-Z)
DIMENSION YLLG(0:100),XLLG(0:100)
INTEGER PT
CHARACTER*2 POS
POS=
PT=0
IF (YL.GT.(AKERB+SCLMIN)).THEN
IF (YL.LT.(YLLG(NLG-1)-(AKERB+SCLMIN))) THEN
IF (XL.GT.(YL*TAN(SANG))) THEN
IF X(XL-SPAN).LT.(YL*TAN(SANG))) THEN
PT =4
IF (YL.LT.0.0) THEN
11=1
495 IF (XL.LT.XLLG(I1)) THEN
MN1=I1
MN2=I1 +1
MN11=MN1+NTG
MN22=MN11+1
XLL=XL-XLLG(I1-1)
YLL = 0
XPL =XLLG(I1)-XLLG(11-1)
YPL=YLLG(1)
ELSE
II =I1 +1
GO TO 495
ENDIF
ELSE IF (YL.GT.YLLG(NLG-1)) THEN
11=1
G=XL-YLLG(NLG-1)*TAN(SANG)
500 IF (G.LT.XLLG(I1)) THEN
MN II = II +NTG*(NLG-1)
MN22 =MN11 +1
MN1=MN11-NTG
MN2 =MN1 +1
XLL=G-XLLG(I1-1)
YLL=YLLG(NLG-1)-YLLG(NLG-2)
XPL=XLLG(I1)-XLLG(I1-1)
YPL=YLL
ELSE
I1 =II +1
GO TO 500
ENDIF
ELSE
11=1
Appendix 11 229
1 SCLMIN=400 PW(I)=1.35
PW(2)=1.35 PW(3)=135
PW(4)=1.35 PW(5)=5.70
PW(8)=5.70 PW(6)=5.70
PW(7)=5.70 PW(9)=3.40
PW(10)=3.40
PW(11)=3.40
PW(12)=3.40
PW(13)=3.40
PW(14)=3.40
PW(15)=3.40
PW(16)=3.40 IF
(LCASE.EQ.2)GO TO 2
RETURN
2 DO 520 1=17,32
520 PW (I)=PW(I-16)
RETURN .
3 SCLMIN=1405
DO 525 1=1,4
PW (1)=2.0
PW (I+4)=3.0
525 PW (1+8)=3.0
DO 530 1=13,28
530 PW (I)=4.25
RETURN
5 SCLMIN=1405
DO 535 1=1,8
535 PW (I)=5.0
RETURN
7 SCLMIN=1620
DO 540 1=1,20
540 PW (1)=3.50
RETURN
9 IF (MOVE.NE.1) GO TO 550
READ (3,*) SCLMIN1
READ (3,*)(PW1(K),K=1,NW)
WRITE(4,545)(PW1(K),K=1,NW)
545 FORMAT('LOAD ON WHEEL' ,10E10.2)
550 SCLMIN =SCLMINI
DO 555 K=1,NW
PW (K) =PWI(K)
555 CONTINUE
10 RETURN
Appendix 11 231
END
c********************************************************i**************
SUBROUTINE LLOAD (LCASE,XT,YT,XW,YW,NW,NWPA,MOVE)
IMPLICIT reals8(A-H2O-Z)
r. DIMENSION XW(100),YIN(100),XW1(100),YW1(100)
-DO 565 11=1,28
XW (II)=0.0
565 YW (11)=0.0
X=XT
Y=YT
GO TO (10,1,3,4,5,6,7,1,9),LCASE
IF AIL=1
RETURN
1 NW =16
NWPA =2
DO 570 11 =1,8'
YW (I1*2-1)=Y
570 YW (11 41)=Y+1800
XW (1)=X
XW (2)=XW(1)
XW (3) =XW(2)-1100
XW (4)=XW(3)
XW (5) =XW(4)-3200
XW (6)=XW(5)
XW (7)=XW(6)-1200.
XW (8)=XW(7)
XW (9)=XW(8)-4300
XW (10)=XW(9)
XW (11)=XW(10)-3000
XW (12)=XW(11)
XW (13)=XW(12)-3000
XW (14)=XW(13)
XV(I5)=XW(14)-3000
XV, (16)=XW(15)
IF (LCASE.EQ.2) GO TO 2
RETURN
2 NW =32
NWPA =4
DO 575 1=17,32
575 XW (1)=XW(I-16)
DO 580 1=9,16
YW (1*2-1)=YW((I-8)*2)+1700
580 YW (I*2)=YW(I*2-1)+1800
RETURN
3 NW=28
NWPA
=4 K-7
585 DO 590 1=1,K
232 Grillage Analogy in Bridge Deck Analysis
YW(I*4-3)=Y
YW(I*4-2) =Y +450
YW(I*4-1)--= Y+ 1930
590 YW(I*4) =Y +2380
IF (LCASE.EQ.5)GO TO
615 595 DO 600 I=1,4
XW(I) =X
XW(I+4) =XW(I)-3960
XW(I + 8) = XW(I +4)-1520
VW (I+ 12) =XW(I +8)-2130
XW(I+16)=XW(I+12)-1370
XW(I+20)=XW(I+16)-3050
600 XW(I +24) = XW(I +20)-1370
RETURN
4 NW=28
NWPA=4
L=7
605 DO 610 I=1,1,
YW(I*4-3)=Y
YW(I*4-2)=Y+795
YW(I*4-1) = Y + 1585
YW(I*4)=Y +2380
610 CONTINUE
IF (LCASE.EQ.6)GO TO 615
GO TO 595
5 NW =8
NWPA=4
K=2
GO TO 585
615 DO 620 I=1,4
XW(I) =X
620 XW(I+4)=X-1220
RETURN
6 NW=8
NWPA=4
L =.2
GO TO 605
7 NW=20
NWPA=2
DO 625 1=1,10
YW(I*2-1)=Y
625 YVV(I *2) =Y +2060
DO 630 1=1,10
XW(I*2-1)=X-((I-1)*457)
630 XW(I*2)=XW(I*2-1)
RETURN
!F. N!'r 1 )r:Ct TO Fall
READ (3,*)NW I
Appendix 1.1 233
IF(NW1.EQ.1)G0 TO 640
DO 635 12P=2,NW1+1
IF (12P.EQ.NW1+1)GOTO 635
READ (3,*)XW1(I2P),YW1(I2P)
635 CONTINUE
640 NW=NW1
XW(I)=XT
YW(1)=YT
DO 645 I1P=2,NW
XW(I1P) =X1V(1)+ XW 1(I1P)
YW(IIP)=YW(1)+YWI(I1P)
645 CONTINUE
I0 RETURN
END
SUBROUTINE DECOMP(NROWS,IWW,A,U)
IMPLICIT real*8(A-H,G-Z)
DIMENSION A(1,50000), U(50000)
DO 650 I=1,NROWS*IWW
650 U(I)=0.0
DO 670 I=I,NROWS
IP = (I-1)*IWW+ 1
SUM=0.0
DO 655 K=1, I-1
IF ((I+1).GT.(K+IWW)) GO TO 655
SUM= SUM + (U((K-1)*IWW+ I-K+1))**2-
655 CONTINUE
U(IP)=SQRT(A(1,IP)- SUM)
DO 665 J= 2,IWW
SUM=0.0
DO 660 K= 1,1-1
IF ((I+J).GT.(K+IWW)) GO TO 660
SUM= SUM + (U((K-1)*IWW +/-K +I))*(U((K-1)*IWW+I-K +1))
660 CONTINUE
U(IP+.1-1)=(A(1,IP+J-1)-SUM)/U(IP)
665 CONTINUE
670 CONTINUE
RETURN
END
c***********************************************************************
SUBROUTINE SOLVE(NROWS,IWW,U,F,D)
IMPLICIT real*8(A-H2O-Z)
DIMENSION U(50000),F(1000),D(I000),X(1000)
X(1)=F(1)/U(1)
DO 680 I =2,NROWS
SUM=0
DO-675 J=1, I-I
IF ((I-J).GE.IWW) GO "ID t,-75
234 Grillage Analogy in Bridge Deck Analysis
SUBROUTINE SORT(IX,JDD,JDD1)
IMPLICIT real*8(A-H2O-Z)
DIMENSION JDD(100),JDD1(100)
DO 720 I=1 ,IX
720 JDD1(I) =JDD(I)
IP=IX+1
Appendix II 235
DO 730 1=1,a-1
IP=IP-1
IK = 1
DO 725 J=1,IP
IF (JDD1(IK).GT.JDD1(J)) GO TO 725
IK
725 CONTINUE
EMP =.113D1(IP)
IDD1(IP)=./DD1(1K)
IDD1(1K) =EMP
730 CONTINUE
RETURN
END
c***************m*******************************************************
SUBROUTINE LDISTR (MN1,M N2,MN11, M N22, PLOAD , SANG,
+ BL,XPL,DLNG,XLL,YLL,POS,NL)
IMPLICIT real*8(A-H, 0-Z)
DIMENSION BL(I000)-
CHARACTER*2 POS
DO 735 I=1 ,3* NL
735 BL(I)= 0
C =XLL
B =YLL/COS(SANG)
D =XPL-XLL
A =(DLNG-YLL)/COS(SANG)
X =XPL
Y =DLNG/COS(SANG)
PE =PLOAD*A*A*(3*B +A)/Y**3
PF =PLOAD*B*B*(3*A+B)/Y**3
BMF =-PLOAD*A*B*B/Y/Y*COS(SANG)
BME=PLOAD*A*A*B/Y/Y*COS(SANG)
BMFI =-PLOAD*A*B*B/Y/Y*SIN(SANG)
BME1 = PLOAD *A*A*B/Y/Y*SIN(SANG)
BL(3*MN I )=PE*D*D*(3*C + D)/X**3-6*BMEI*D*C/X**3
BL(3*M N2) = PE*C*C*(3*D + C)/X**3 +6*BMEI*D*ca**3
BL(3*MN11)= PF*D*D*(3*C + D)/X**3-6*BMF1*D*C/X**3
BL(3*MN22)=PF*C*C*(3*D +C)/X**3 +6*BMF1*D*C/X**3
BL(3*MN1-1)=-(PE*C*D*D/X/X-BME1*D*(2*C-D)/X**2)
BL(3*MN2-1)= -(-PE*C*C*D/X/X-BMEl*C*(2*D-C)/X**2)
BL(3*MN11-1) --(PF*C*D*D/X/X-BMFI*D*(2*C-D)/X**2)
BL(3*MN22-1) =-(-PF*C*C*D/X/X-BMFI*C*(2*D-C)/X**2)
BL(3*MNI-2) =BME*D/X
BL(3*MN2-2) =BME*CIX
BL(3*MN I1 -2) =BMF*D/X
BL(3*M N22-2) = Bne*ca
RETURN
END
c***********************************************************************
236 Grillage Analogy in Bridge Deck Anidysis
DO 770 J=1,NLG
IF (I.EQ.N/NLG*J)G0 TO 785
IF (I.EQ.N/NLG*J+1)G0 TO 785
770 CONTINUE
DO 775 J=1,NLG
IF (I.LT.N/NLG*J)G0 TO 780
775 CONTINUE
780 CALL INODE(M,I-1,I,I1)
CALL INODE(M,I,I+1,12)
AL= (PM(I1, 1) +PM(I2, 1))/2.
ALl = PM(I1,1)/2.0
AL2 = PM(I2,1)/2.0
B(3*1)=B(3*I) +K(J*2-1)*AL
B(3*I-1)=B(3*I-1)-K(J*2-1)*(-PM(I1,1)**2+PM(12,I)**2)/12
B(3*I-2)=B(3*I-2)+K(J*2)*AL
785 CONTINUE
DO 790 I=1,NLG
B(3*((1-1)*N/NLG+I))=El*K(2*I-1)+B(3*((I-1)*N/NLG+1))
790 B(3*N*IINLG)=E2*K(2*I-I)+B(3*N*I/NLG)
ELSE
READ(3,*)(B(3*I),I=1,N)
WRITE(4,8I0)
WRITE(4,815)(I,B(3*1),I=1,N)
ENDIF
795 FORMAT(/5X,'DEAD LOAD ALONG LONGITUDINAL GRID LINESV5X,39c*'),
+ I/5X,'LONG. GRID',5X,'VERTICAL LOAD',5X,'TORS1ONAL MOMENT'/5X,
+ 'LINE NO.',13X,'TIMM',14X,'T MMTMW/60('*')/)
800 FORIvIAT(8X,I3,9X,E11A,9X,E11.4)
805 FORMAT(!/'THE END PROJECTIONS OVER SUPPORT LINES ARE ',E1I.4,
+' & A,' MM'/' RESPECTIVELY.'//)
810 FORMAT(/5X,'DEAD LOAD AS SUPPLIED,DIRECTLY ON NODES,VERTICAL'
+ ,'. LOAD ONLY',/5X,30('*')//5X,'NODE NO.' ,10X,'LOAD(T)'
+ /30('*'),/)
815 FORMAT((6X,I3,I3X,E11.4))
RETURN
END
c***********************************************************************
SUBROUTINE COMPARE1(PT,B,RT,MAXPT,MAX13,MAX.RT,MAX RTN,MM,N
+ ,IX,IMI,PT1,X,Y
+ ,LCASE,CORSPT,CORX,CORY,CORX1,CORY1,CORRX,CORRY,CORD1
+ ,JTT)
IMPLICIT REAL*8(A-F1,0-Z)
DIMENSION PT(1000,7),B(1000),RT(100),PT1(1000,2),JTT(100)
REAL *8 MAXPT(6,4,500),MAXB(4,500),MAXRT(4,100),MAXRTN(4,100)
+ ,CORSPT(2,500),CORX(100),CORY(100),CORX1(100),CORY1(100)
+ C ORRX(2,500). CORRY(2 .500) CORD I
(100) IF (1M1.GT.1)G0 TO 850
Appendix II 239
DO 820 I=1,4
1 =I +2
DO 820 K=I,MM
MAXPT(J, 1,K) = PT(K, J)
MAXPT(I,2,K)=LCASE
MAXPT(L3,K) = X
MAXPT(I,4,K)=Y
CORSPT(I,K) =PT(K,4)
CORSPT(2,K)=PT(K,3)
820 CONTINUE
KL
DO 825 I=1,N
N11 =3*I-2
N13 =3*I .
CORRX(1,I) -=13(N11 +1)
CORRX(2,I)=B(N13)
CORRY(1,I) =B(N 11)
CORRY(2,I)=B(N13)
DO 825 J=N11,N13
MAXB(1,J)=B(J)
MAXB(2,J)=LCASE
MAXB(3 ,J) =X
MAX13(4,J)=Y
825 CONTINUE
), DO 830 I=1,IX
MAXRT(I ,I) =RT(I)
MAXRT(2,I) =LCASE
MAXRT(3,I) =X
MAXRT(4,I) =Y
CORXI (I) = B(JTT(I)*3 -2)
CORYI (I) =B(JTT(I)*3-1)
CORD1(1) = B(ITT(I)*3)
MAXRTN(1,I) =RT(I)
MAXRTN(2,I)=LCASE
MAXRTN(3,I) =X
MAXRTN(4,I) =Y
WRITE(*,*)LJTT(I),B(JTT(I)*3-2),B(JTT(1)*3-1)
CORX(I)=B(JTT(1)*3-2)
CORY(I)=B(JTT(I)*3-1)
830 CONTINUE
DO 835 1=5,6
DO 835 K=1,MM
J=1-4
MAXPT(I, 1,K) = PT1(K,J)
MAXPT(1, 2,K) = LCASE
MAXPT(1,3,K)= X
MAXPT(I,4,K)=Y
835 CONTINUE
240 Grillage Analogy in Bridge Deck Analysis
DO 845 1=3,6
DO 845 K=1,MM
IF (1.GT.4)G0 TO 840
IF (MAXPT(I,1,K).LT.0.0)THEN
MAXPT(I,1,K)=0
MAXPT(I,2,K)=0
MAXPT(I,3,K)=0
MAXPT(I,4,K) =0
ENDIF
GO TO 845
840 IF (MAXPT(I,1,K).GT.0.0)THEN
MAXPT(L1,K)=0
MAXPT(I,2,K)=0
MAXPT(I,3,K)=0
MAXPT(I,4,K) =0
ENDIF
845 CONTINUE
RETURN
850 DO 855 1=1,2
J=I+2
DO 855 K=1,MM
IF (ABS(MAXPT(I,1,K)).GE.ABS(PT(K,J)))GO TO 855
MAXPT(I, 1,K) = PT(K,J)
MAXPT(I,2,K) =LCASE
MAXPT(1,3,K)=X
MAXPT(I,4,K) =Y
CORSPT(I,K)=PT(K,54)
855 CONTINUE
DO 860 1=3,4
J=1+2
DO 860 K=1,MM
IF (MAXPT(I,1,K).GE.PT(K,J))GO TO 860
MAXPT(1,1,K)=PT(K,J)
MAXPT(I,2,K) =LCASE
MAXPT(1,3,K)=X
MAXPT(1,4,K)=Y
860 CONTINUE
DO 865 1=5,6
DO 865 K=1,MM
J=I-4
IF (MAXPT(1,1,K).LT.PT1(K,J))G0 TO 865
M A'XPT(I , 1 ,K) = PT1(K,J)
MAXPT(1,2,K)=LCASE
MAXPT(1,3,K)=X
MAXPT(1,4,K)=Y
865 CONTINUE
KL1 = I
DO 870 I=1,N
Appendix II 241
N 1 =3*1-2
IF (ABS(MAX13(1,N11)).GE.ABS(B(N11)))G0 TO 870
MAXI3(1 ,N 1) =B(NI1)
MAXI3(2,N11)=LCASE
MAX13(3,N11)=X
MAX13(4,N11)=Y
CORRX(1 ,I) = B(N1 I +1)
CORRX(2,I)=B(N11+2)
870 CONTINUE
KL1 =1
DO 875 I=1,N
N12=3*1-1
IF (ABS(MAXB(1,N12)).GE.ABS(B(N12)))G0 TO 875
MAXB(1,N12) =-13(N12)
MAX13(2,N12)=LCASE
MAXB(3,N12)=X
MAX13(4,N12)=Y
CORRY(1,I)=B(N12-1.)
CORRY(2,1)=B(N12+1)
875 CONTINUE
DO 880 I=1,N
N13 =3*I
IF (ABS(MAX13(1,N13)).GE.ABS(B(N12)))G0 TO 880
MAXB(1,N13)=B(N13)
MAXB(2,N13)=LCASE
MAXB(3,N13)=X
MAXB(4,N13)=Y
880 CONTINUE
DO 885 I=1,IX
IF(MAXRT(1,1).GT.RT(I))G0 TO 885
MAXRT(1,I)=RT(I)
MAXRT(2,I) =LCASE
MAXRT(3 ,I) =X
MAXRT(4,I)=Y
CORX1(I) = B (JTT(I)*3 -2)
CORY1(I) = B (JTT(I) *3-1)
CORD1(I)=B(ITT(I)*3)
885 CONTINUE
DO 890 I=1,IX
IF (MAXRTN(1,I).LT.RT(I)) GO TO 890
MAXRTN(1,I)=RT(I)
MAXRTN(2,I)=LCASE
MAXRTN(3,I)=X
MAXRTN(4,I)=Y
WRIT E(*,*)i,yrixo, B(JTT(1)*3-2),B(JTT(1)*3-1)
CORX(I)=B(ITT(I)*3-2)
CORY(I)=B(JTT(1)*3-1)
890 CONTINUE
242 Grillage Analogy in Bridge Deck Analysis
RETURN
END
c***********************************************************************
SUBROUTINE COMPARE2(PT,B,RT,RTN,MAXPT,MAXB,MAXRT,MAXRTN,
+ MM,N,IX,IM1,PT1,
+ MAX PT2 ,MAXB2,MAXRT2;MAXR2N,JTT ,COR2X,COR2Y,C OR2XI ,COR2Y1 ,
+ COR1X COR I Y, COR1XI ,CORlYI,CORSPT1,CORSPT2,
+ CORRRX I , CORRX2 ,CORRY I ,CORRY2 , CORD,CORD2)
IMPLICIT REAL*8(A-H2O-Z) 90
DIMENSION PT(1000,7),B(1 000), RT(100),RTN(100),PT I (1000, 2),
JTT(100)
REAL*8 MAXPT(6,4,500),MAXB(4,500),MAXRT(4,100)
+ ,MAXPT2(6,2:4,500)
+ ,MAXB2(2:4,500),MAXRT2(2:4,100),MAXRTN(4,100),MAXR2N(2:4,100)
+ ,COR2X(100), COR2Y(100),COR2X1(100) ,COR2Y1( I00) ,CORSPT1 (2 ,500)
+ ,CORSPT2 (2,500),COR1X( I00),CORI Y(100),COR1X1(100),COR1Y1 (100)
+, CORRX1(2,500), CORRX2(2,500) ,CORRY1(2,500) ,CORRY2 (2,500) 9
+, CORD( I00),CORD2(100)
IF (Im1.NE.1)G0 TO 925
DO 895 I =1,4
J=I+2
DO 895 K = 1,MM
MAXPT(I, I ,K) = PT(K,J)
MAXPT(I,2,K)=MAXPT2(1,2,K)
MAXPT(I,3,K) = M AXPT2 (I ,3,K)
MAXPT(I,4,K) =MAXPT2(I,4,K)
CORSPT1 (1 ,K) = CORSPT2 (1,K)
CORSPT I (2 ,K) = CORSPT2(2,K)
895 CONTINUE
DO 900 I=1,N
NI I =3*I-2
N13 =3*1
DO 900 J=N11,N13
MAXB(1,J) = B(J)
MAXB(2,J) = M AXB2(2 ,J)
MAXB(3,J) MAXB2(3 ,J)
MAX B(4, J) = M AXB2 (4, J)
CORRX1(1 ,I) =CORRX2(1,I)
CORRY I (1 ,I) =CORRY2(I,I)
C ORRX I (2 , I) = CORRX2 (2,I)
C ORR Y I (2,D= CORRY2(2,I)
900 CONTINUE
DO 905 I = 1,IX
MAXRT(1,1) = RT (I)
M AXRT(2 , I) = M AXRT2 (2,1)
M AXRT (3 , I) = MAXRT2 (3,1)
M AXRT (4, I) = MAXRT2(4,I)
COR1XI(I)=COR2X1(I)
Appendix II 243
CORI Y f(I) =COR2Y1(1)
CORD(I)=CORD2(I)
MAXRTN(1,I) =RTN(I)
MAXRTN(2,I) =MAXR2N(2, I)
MAXRTN(3,I)=MAXR2N(3,I)
MAXRTN(4,I)=MAXR2N(4,1)
COR1X(I)=COR2X(I)
CORIY(I)=COR2Y(I)
905 CONTINUE
DO 910 1=5,6
DO 910 K=1,MM
J=I-4
MAXPT(L1,K)=PT1(K,J)
MAXPT(I,2,K) =MAXPT2(I,2,K)
MAXPT(1,3,K)=MAXPT2(1,3,K)
MAXPT(L4,K)=MAXPT2(I,4,K)
910 CONTINUE
DO 920 I=3,6
DO 920 K=1,MM
IF (I.GT.4)G0 TO 915
IF (MAXPT(I,1,K).LT.0.0)THEN
MAXPT(I,1,K)=0
MAXPT(I,2,K)=0
MAXPT(I,3,K)=0
MAXPT(1,4,K)=0
ENDIF
GO TO 920
915 IF (MAXPT(1,1,K).GT.0.0)THEN
MAXPT(I,1,K)=0
MAXPT(I,2,K) =0
MAXPT(1,3,K)=0
MAXPT(I,4,K)=0
ENDIF
920 CONTINUE
RETURN
925 DO 930 I=1,2
.1= I + 2
DO 930 K=I,MM
IF (ABS(MAXPT(I,I,K)).GT.ABS(PT(K,J)))GO TO 930
MAXPT(I,1,K)=PT(K,J)
MAXPT(1,2,10=MAXPT2(1,2,K)
MAXPT(I,3,K)=MAXPT2(I,3,K)
MAXPT(I,4,K)=MAXPT2(I,4,K)
CORSPT1(I,K)= CORSPT2(I,K)
930 CONTINUE
DO 935 1=3,4
J=1+2
DO 935 K=i,iviNi
244 Grillage Analogy in Bridge Deck Analysis4
IF (MAXPT(I,1,K).GT.PT(K,D)G0 TO 935
MAXPT(I,1,K)=PT(K,J)
MAXPT(I,2,IC) =MAXPT2 (I, 2, K)
MAXPT(1,3,K)=MAXPT2(I,3,K)
MAXPT(I,4,K)=MAXPT2(I,4,K)
935 CONTINUE
DO 940 1=5,6
DO 940 K=1,MM
J=I4
IF (MAXPT(I,1,K).LT.PT1(K,i))GO TO 940
MAXPT(I, I ,K) =PT1(K,J)
MAXPT(1,2,K)=MAXPT2(1,2,K)
MAXPT(1,3,K) =MAXPT2(1,3 ,K)
MAXPT(1,4,K)=MAXPT2(I,4,K)
940 CONTINUE
KL1=1
DO 945 I=1,N
N11=3*1-2
IF (ABS(MAXB(1,N11)).GT.ABS(B(N11)))G0 TO 945
MAX'S(1 ,N11)= B(N I 1)
MAXB(2, N11) = MAXB2(2, N11)
MAXB(3,N11)=MAXB2(3,N11)
MAXB(4,N1I)=MAX132(4,N11)
CORRX1(1,1)=CORRX2(1,I)
CORRX1 (2 ,I) CORRX2(2 ,I)
945 CONTINUE
DO 950 1=1,N
N12=3*1-1
IF (ABS(MAXB(1,NI2)).GT.ABS(B(N12)))GOTO 950
MAXB(1,N12)=B(N12)
MAXB (2, N12) = MAXB2(2 , N 12)
MAXB(3,N12)=MAX132(3,N12)
MAXB(4,N12) = MAXB2(4,N 12)
CORRY I (1 ,1) = CORRY2(1, I)
CORRY I (2,I) = CORRY2(2, I)
950 CONTINUE
DO 955 I=1,N
N13 =3*I
IF (ABS(MAXB(1,N13)).GT.ABS(B(N13)))GOTO 955
MAXB(1,N13)=B(N13)
MAXB(2,N13)=MAXB2(2,N13)
IVIAXB(3 , N13) =MAXB2(3 , N13)
MAXB (4 , N13) = MAXB2(4, N13)
955 CONTINUE
DO 960 I=1,IX
IF(MAXRT(1,I).GT.RT(I))GO TO 960
MAXIIT(11)=RT(T 1
MAXRT(2,I)=MAXRT2.(2,I)
Appendix11 245
MAXRT(3,I)=MAXRT2(3,I)
MAXRT(4,I)=MAXRT2(4,I)
COR1XI(I)=COR2X1(1)
COR1Y1(I)=COR2Y1(1)
CORD(I) =CORD2(I)
960 CONTINUE
DO 965 I=1,1X
IF(MAXRTN(1,11.1-T.RTN(1))00, TO 965
MAXRTN(1,I)=RTN(I)
MAXRTN(2,I)=MAXR2N(2,I)
MAXRTN(3,I) =MAXR2N(3,I)
MAXRTN(4,I)=MAXR2N(4,i)
COR1X(I)=COR2X(I)
COR1Y(1)=COR2Y(1)
965 CONTINUE
' RETURN
END
c***************************************************************i*******
SUBROUTINE WRITE2 (MAXPT,MAXB,MAXRT,MAXRTN,MM,N,IX,MOVE,M,
+ JTT,LCASE,COR1X,CORlY,CORIX1,COR1Y1,CORSPT1,CORRX1,CORRY1
+ ,CORD)
IMPLICIT REAL*8(A-H2O-Z)
DIMENSION M(1000,2);JTT(100)
real*8 MAXPT(6,4,500),MAXB(4,500),MAXRT(4,100),MAXRTN(4,100)
+ ,CORIX(100),COR1Y(100),COR1X1(100),CORIY1(100),CORSPTI(2,500)
+,CORRXI(2,500),CORRY1(2,500),CORD(100)
WRITE (4,1025)
DO 970 I=1,MM
WRITE (4,1030)I,M(L1),M(L2),(MAXPT(K,1,0,CORSPT1(K,1),
+(MAXPT(K,J,1),J= 2,4),K = 1 ,2)
970 CONTINUE
WRITE (4,1015)
WRITE (4,*)'MAX SAGGING MOMENT ON THE ELEMENTS'
WRITE (4,1040)
DO 975 I=1,MM
975 WRITE (4,1035)I,M(I,I),M(I,2),(MAXPT(5,7,1),J= I,4),(MAXPT(4,5,1)
+ ,J=1,4)
WRITE (4,1015)
WRITE (4,*)'MAX HOGGING MOMENT ON THE ELEMENTS'
WRITE (4,1040)
DO 980 I= I,MM
980 WRITE .(4,1035)I,M(I,1),M(1,2),(MAXPT(3,J,I),J = 1,4), (MAXPT(6,1,1)
+ ,J = 1,4)
READ (3,*) IK4
IF (IK4.NE.1) GO TO 990
WRITE (4,11345)LNAME,LNAME,LNAME
DO 985 I=1,N
N11=3*1-2
246 Grillage Analogy in Bridge Deck Analysis
N12 =31-1
N13 = 3*I
WRITE (4,1050)1,1, (MAX.13(J,N11),J =1,2),(MAX13(1,NI2),J= 1,2),
+ (MAXB(J,N13),J = 1,2)
985 CONTINUE
990 CONTINUE
WRITE (4,1060)
DO 995 I=1,LX
WRITE (4,1065)LITT(1),-MAX.RTN(1,I),COR1X(I),CORlY(I)
+ ,MAXB(I,3*JTT(I)),(MAXRTN(J,I),J =2;4)
995 CONTINUE
WRITE (4,1020)
WRITE (4,1055)
DO 1000 I=1,IX
WRITE (4,1065) LITT(1),-MAXRT(1,1),COR1X1(1),CORIY1(1),CORD(1)
+ ,(MAXRT(I,I),J=2,4)
1000 CONTINUE
WRITE (4,1075)
DO 1005 I=1,IX
N1I=3*.ITT(1)-2
WRITE (4,1085)LITT(1),MAXB(1,N11),CORRX1(1,JTT(1)),
+ C ORRX1(2, ITT(I)),(MAXB(J, N11), I =2,4)
1005 CONTINUE
WRITE (4,1080)
DO 1010 I=1,IX
N12=34TIT(1)-1
WRITE (4,1085)I,JTT(I),MAXB(1,N12),CORRY1(1,ITT(I)),
+CORRY1(2,Ms(I)),(MAXB(J,N12),J=2,4)
1010 CONTINUE
WRITE (4,1020)
1015 FORMAT(/80('*'))
1020 FORMAT(/80('*'))
1025 FORMAT(//,'MAX. SHEAR AND TORSION ON THE ELEMENTS'/
+110('*')/'SNO ELEMENT MAX SHEAR TORSION'
+,' LCASE X Y '
+ ,' MAX TORSION SHEAR',
+ ' LCASE X Y'/19X,'T TMM',12X,'MM MM',5X,
+ TMM T',14X,'MM MM'/110('*')///)
1030 FORMAT(I3,3X,I2, I2,2(2X, E11.4, IX, E11.4,
1X +,F3.0,2X,F6.0,2X,F6.0))
1035 FORMAT(13,3X,I2,'-',12,2(2X,E11.4,1X,F3.0,2X,F6.0,2X,F6.0))
1040 FORMAT(80('*')PSNO ELEMENT AT END1 LCASE X Y '
+ ,' AT END2 LCASE X Y'/18X,'TMM',13X,'MM MM',5X,
+ 'TMM',14X,'MM MM'/80('*')/) 1045
FORMAT(99('*')/'SNO',2X,'JOINT',3X,
+ 'MAX X-ROTATION',1X,1A8,2X,'MAX Y-ROTATION',1X,1A8,2X,
+ 'MAX Z-DEFLECTION',1X,1A8,/22X,'RAD'.24X,'RAD',27X,
'MM',/99('*')/)
Appendix II 247
1050 FORMAT(I3.4X,13,2(4X,E13.6,5X,F6.0),5X,E13.6,6X.
+F6.0)
1055 FORMAT(/,'MINIMUM REACTION ON SUPERSTRUCTURE/K('*')
+PSNO JOINT REACTN COR ROT-X COR ROT-Y
+ CO R D EFL LC A S E X r/
+ 80('*')/)
1060 FORMAT(//,`MAXIMUM REACTION ON SUPERSTRUCTURE/80(n
+PSNO JOINT REACTN COR ROT-X COR ROT-Y COR DEFL
, ' LC A S E X Y '/
+ 80('*')/)
1065 FORMAT(I3,2X,I3,2X,4(E11A,1X),F3.0,2X,F6.0,2X,F6.0)
1070 FORMAT(13,2X,13,2X,3(E11.4,1X),F3.0,2X,F6,0,2X.F6.0)
1075 FORMAT(80(w),/,' SNO JOINT MAX ROT X COR ROT Y '
+ COR DEFL LCASE X Y 'PUNITS IN T AND mm'/
+ 80(P')/)
1080 FORMAT(80C*),/,' SNO JOINT MAX ROT Y COR ROT X '
+,' COR DEFL LCASE X Y '/
+ 80('*')/)
1085 FORIvIAT(I3,2X,I3,3(3X,E11.4),2X,F3.0,2X,F6.0,2X,F6.0)
1090 FORMAT(I3,2X,I3,2X,E11.4,6X,F3.0,2X,F6.0,2X,F6.0)
RETURN
END
c**********************************************************************
SUBROUTINE WRITE1(13,MAXPT,MAXB,IvIAXRT,MAXRTN,MM,N,IX,MOVE,M,
TIT,LCASE,CORSPT,CORX,CORY,CORX1,CORYLCORRX,CORRY,CORD1)
ti
Subject Index
Analysis, 35, 111, 131-139 Clearance, 20, 22-24
Arbitrary planforms, 16 Codes of practices (see bridge codes)
Assumptions, 36, 56 Composite construction, 10
Axes, 42, 46 Coarse mesh, 41, 171
Axle loading, 19 Contributory area (see tributory area)
Computer aided methods, 1, 3
Band-width, 51, 64 Computer programs
C Beam GRID, 32, 66-69, 205-212
contiguous beam, 7, 87 GABS, 32, 146-169, 171-203,
effective flange width, 96, 101 213-250
V equilibrium equations, 74 Contiguous beams (see beams)
spaced beam, 14 Contiguous nodes, 115, 117
spine beam, 14 Courbon's method, 1, 36-38
Bearing stiffness, 133 Curved deck, 42
Bearings Curved geometry, 16
neoprene, 16, 32, 133
rocker, 16 Data
roller, 16 input, 150, 165
Bending moment diagram, 140, 143 output, 150, 165
Box-girder bridges, 6, 12-15, 92-94, Dead Ioad, 18, 112, 118
103-109, 143-146 Degrees of freedom, 46, 56-58
Bridge classification, 5 Design curves, 1, 39, 43, 44
Bridge codes, 5, 17, 18 Design envelopes, 31, 138
Bridge deck analysis (see analysis) Diaphragms, 10, 13
Bureau of Indian Standards (BIS), 17 Direct stiffness method, 31, 55-73, 132._.-
Displacement method, 56
Cable stayed bridges. 3, 6 Distortion, 13, 93, 108, 144
CAD/CAM systems, 4 Ductility-ratio, 4
Cantilever bridges, 16
Cantilever construction, 5 Cell Edge grid lines, 87, 89
distortion, 13, 93, 108, 144 Effective width of flange, 96, 101
Cellular deck, 6, 51, 103 Elastic properties, 94-109
(`am not , Anri nr.rvisrdaCci no,
(.1 Envelope diagrams, 134, 138, 196
CA), 8 Equilibrium equations, 74
Cho lesky' s factorisation Equivalent loads (see transfer of
method, 31, 69, 133 loads)
252 Subject Index
Examples, 69-73, 171-203
Expert system, 4
External prestressing, 3
Gauss-Elimination procedure,
31, 66, 133
Gauss-Seidel method, 65
Grade separation, 12
Grillage examples
box-girder deck, 143, 199-203
grid deck, 69-73
skew deck, 69-73, 179-182
slab deck, 172-187
slab-on-girders deck, 187-198
voided slab deck, 183-187
Grillage section properties
box-girder-deck, 103-109
cellular deck (see box-girder deck)
skew deck, 101
slab deck, 99-101
slab-on-girders deck, 101-103
Launching, 5
Limit state concept, 4
Lines of strength, 78
Live load
AASHO, 28, 147
IRC, 18-25, 28-30
OHBD, 147, 191
Loading (see dead load, live load,
impact load)
Load distribution, 35, 43, 99
Loading standards, 17, 24, 28, 30
Local effects, 111
Longitudinal grid lines, 31, 78
Neutral axis. 99
Nodal deformations, 66, 131, 133
Nodal loads, 31, 69, 117
Re-entrant corners, 97
Responses, 31, 32, 111, 119, 134, 147
Result output, 150, 165
Rigid frame bridges, 6
Rigid supports, 16, 65, 133
Saint-Venant, 40, 96
Segmental construction, 5
Shear deformation, 13, 92, 101 Vehicular live loading (see also live
Shear lag, 13, 92, 101, 172 load), 27, 31, 111, 115, 119
Shear modulus, 70, 76, 166 Voided slab (see also slab bridges),
Sign convention, 56, 150 6, 13, 87, 99, 146, 182
Skeletal structures, 46, 56
Skew angle, 80 Warping, 172
Skew decks
finite difference method, 41, 42 Young's modulus, 70, 76, 166
finite element method, 45
thccry, '!"
Subject Index 253
Skew grid, 31, 42, 69-73, 80
Skyline technique, 64
Slab bridges, 6-9; 79-89, 99-101, 139-
142
Slab-on-girders bridges, 6, 9-12, 90-
92, 101-103, 143
Slope deflection, 74, 75
Solid slab, 79, 80, 90, 99, 146
Spaced box-girders, 94 Space
frame analysis, 46 Spine
beam bridges, 14 Spine box-
girders, 14, 94, 146
Stiffness method (see direct stiffness
method)
Subroutines (see computer programs)
Support conditions, 2, 5, 31, 35, 132
Supports, 6, 13, 15, 69