Analysis and Tests of A Cylindrical Shell Roof Model. Proc. ASCE
Analysis and Tests of A Cylindrical Shell Roof Model. Proc. ASCE
Lehigh Preserve
Fritz Laboratory Reports Civil and Environmental Engineering
1954
B. G. Johnston
Recommended Citation
Thurlimann, B. and Johnston, B. G., "Analysis and tests of a cylindrical shell roof model., Proc. ASCE Separate No. 434, 1954, Reprint
No. 95 (54-2)" (1954). Fritz Laboratory Reports. Paper 1469.
http://preserve.lehigh.edu/engr-civil-environmental-fritz-lab-reports/1469
This Technical Report is brought to you for free and open access by the Civil and Environmental Engineering at Lehigh Preserve. It has been accepted
for inclusion in Fritz Laboratory Reports by an authorized administrator of Lehigh Preserve. For more information, please contact
[email protected].
.' . ~ .. . .
Al'TALYSIS A~'D 'IESTS OF' A
CYL!--NDRICAL SHELL ROOF MODEL
by
and
Bruce G. Johnston"C--:l-
Synopsis:
2 13K -2-
ip the Uni ted Sta tes. The Ice Arena in Hershey, ?ennsylv ani a, 1936
( 6) *..~, wa s the f:1r st c:v lindri cel she 11 roof in the Uni ted Sta tes.
In the following years and especially during and after the war,
~r ell4 M pie
Lateral horizontal forces, caused) ~ )by .the wind action on X
the front and back door of large span cylindrical shell roofs,
(Fig.l) can reach a considerable magnitude. Their influence on
the structure can no longer be neglected as secondary. A difference
of tempera tUre inside and outside of the bu1lding may produce rela-
the problen s was made and the results were checked experimenta lly
* The Project was started Jan. ~, 1949 and terminated Feb. 15, 1951.
2l3K --4-
b<N<f lx=o - ) N1 dx
(Ring) (Cylinder)
the stresses in the ribs are obviously identical for both structures
under an equal load systen. Once the width b, called the effective
wid th of the cylinder, is known, the calculation of the rib stresses
becomes a problem of tne "Simple Beam Thea ry" (cross sections rems in
plene) 6
-
* Only a very ahort outline of this problem can be given here. For
s complete study references (15) and (16) should be consulted.
r ..'
de~lved
expressions
-~ "
I indeflnitiGly
.
special caae (
on either
.
side of a rib..
(3f-:: 00 In Fig .. 4>
This 113 a
of the author's
...'
,,'
2l3K 5
Where: b
a
= effective' width
\
"k- n~
8 ,
where n depends on the variation o~ the stresses in circumfer-
entia1 direction (n is' the number of complete waves made by e
harmonic function ,around the circumference of the cylinder, see
Fig.4). For n=O, the stress distribution has axial s~etry
(e.g. constant radial 'line load around the rib) Bnd the effective
width has its ne ximum value. The parameterA depends on the
length of the overhang ~ of the are 11 Bnd on the shell cons tan tx.
/ i: 4{3 * (3)
{Sll
the left side. This may be safely done, if the factor I d > 20<4,
where d is t~ distance of the end from the rib (see Ref. (15)
p. 25).
minimum.
The action of the wind on the front end back door of cylin-
drioal shell roofs (Fig. -1.) produces considerable horizontal force.
1.) Problem:
213K 7
~~x __
-
~it
1..
(7 ).
Of.)'
By integra ting Eq. (7) end c m sidering the fact that the shear
force N is zero at the center of the spen, due to s'YJIlmetry
xC;f
conditions, N,x becomes
N yx =-h p l'
- -'
*The erect1. on of cylindrical s1::811 roof s is mede in un1 ts in order
to re-use the form work Elnd furthermore to , provide diletatioD .1oints o
,
*~fThis assumption holds certe1 nly for shell roofs whose width is'
relatively small to the span L.
.i.
213K 8
the load ,P and subjected to' shear forces NX <1'snd N ~x; and di,rec·t,
forces Nx only. ActuallY ·thel!!e Imaginary fp'rce: El arenpt ~cting;. t\,
Introducing shear forQ88 -::N xr s<:tlng ~long' the, connecting I,i,ne, (.:) .'
of the ::ell .and the' ri ba th}ese ima~lnary forces N.xp are ellm~ \( ~"X'
instad. The interaction between the. rib and 'the shell
~. . is ts ken
. - - - -' - .-
, I
geometrical c ondi ti 0 ns fur ,n1 shing the two equa t ions far the
J (10)
CP - (,p.
J - J
•
K •
JK ~ (JK)~eBS
, 'K S (¢K)meas -}:MK
(J X)mea 8 and (¢K)meas ere the mea aured horizontal d1 apIa cement
and rotation of the abutment respectively. Kbeing the coefficient
of e18s~ie
restraint, the product -tM , where M is the restraining
K K
moment, represents the elastic end rotation of the rib. *
In general, the shspe of the center line of the rib ani the
and Me. For the special CElse of a circular center line and e con-
stant I, the calculetion of Hand M by direct integration is ahown~
c c
Extension easilyma y be made to other cases, e.g. parabolic shape
of the rib and I = Ie , or to s numeri cs I solu ti on.
cos~
For the rlb shown in Fig. 8 the normal force No end the be!lding
moment 1.1 ' due to the sl'E-a r load Nx~' may be obtained by integrati~g
0
the contribution of t~ di stri buted shear loed" Ifris the angle at
which No end Mo are determined, let Nx~oadu.>represent the sh:la~
loe d applied over an infinitesimal rib length adW a t any angle t.J,.) 'f
between 0 and r:t. 'Then: (Fig. 8b)
NQ ., - xr • co. ('1'-"" dW
Substituting. from Eqc (9) - Nxces -!-P1' and performing the in-
J 21.
tegrat10n
\ The total normal force N and the to ta 1 bending moment M are hence
(F'ig. 8e)
. N =- No - Hc co~ '! (14)
system is expressed as the work d one by ·-the virtual .loed P' : 1 * :::
(
d k
~jM'EI
tMd~S +fN EA'Nds
Where: ds = Sed t:f
The influence of the normal force N on the displacement is extremely
small and therefore, disregarded. Assuming, EI is constant, end
·S~e eny textbook about Elastic Strain E Methods
213K 11 .
J[-
•
becomes: f
J;.: k "e {eo0'f - eO"'/kJ &0 t He".(1-eos1 lt Me J"od r
Replacing M~ by (13) and integrating gives
(18) and (20) are two equations with the 2 unknowns He and M •
c k dl
am ¢k are given by (10) or (11) respectively. Once the r·edundante
are known the normal force N and the bendi ng moment M for any angle ~
are given by (14) and (15).' The stresses in the ribs are determined
by the usual formula for combined direct force N and the bending
moment M
"..._ N i
u - ~ rMz
213K 12
A and I are the cross sectional area snd the moment of inertia
of the effeotive section,- cmsisting of the rib ss s web and the
In the first case the NCf - . forces of the shell will be taken by
the support and the effective width 1s constant down to the edge-
member (Fig. 9a). If there 1s no support at all, the effective
width must reduce to zero at thesprlnging line. (f'ig.9b). The
calcula ted by means of Eq. (21), once th e ncr rna 1 force Nand. the
S &" bh(\A
the present case the boundary condi.ticns for the shell at the ribs·
x .: 0 and x = 21, are:
x =~: M
X
e 0
X 15 0 S = given
Mx being the bending moment of the shell in axial direction per
unit width, the conditionM x • 0 presupposes the neglecting of, the
the forces end moments being rapidly damped out with increasing
secti on, the normal force N and the bending moment M for any engle l'
along the rib, the deflect! on of the rib is readily determined 'by
the work equat ion. Putting B dummy load P t =1 at th('J point for
which the deflection is (lestred, the work of this unit load due to
r
deflection under cons ideration, or
d= MtMds
EI
"'fNtNgs
EA
2l3K l<t
where Nt Bnd M' are due to P' ~ 1, end Mend N due to the actual
load Sj"stem.
6.) Numerical Example:
The foregoing theory is now applied to the analytical solution
of a model of a shell roof construction. }ilig. 10 g1 ves the dimensions
of the model approximately 1/30 of that of the actual structure
(e.g. hangar at Rapid City, S.D.). Structural steel was chosen as
materiel for the model, the dimensions (shell thickness about 1/8")
leeving practically no other choice. The following list gives the
principal dimensions and properties:
List:
Outer ribs: Height r :: 2.113"
Thickness t e 0 0505"
Angle 'fk 81
0.5866 radians
Middle rib: Height r .:' 2.113"
(Used la ter on)
'I'hickness t :: 0.625"
Shell: Radius a = 108 ft
Thickness h - 00118"
Distance 1 l~ 12"
Load: '* p :: 1~0.62 Ib/in
coeff1 c1 ent
1 0 3161 s 0 .. 3687
1/ sh'
To know the number n, the stress distribution along the ri.b
"
f1 nal. results, Fig. 12 shows the stress in the. ribs. The stress
Vt (vA has the same variation) in the lower fiber has approximately,
tb. e fo I'm. of a cosine-funct:l. on, wi th a helf-wave length of ¢: 0035 .
to 0.40 radians o !he number n becomes
n = 1f_"~_~~_/::::-'
half-wave length
A,:: n~ s 0 0 26
Ent~rlng the chert of Fig. 4 one may read out the coeff:i.cie- nt K:
For:/! =0
1...: 0 0 26
section are calcu18ted next. Fig. lIe gives all pertinent 1n-
format ion:
Area A =1.227 in 2
Distances z z
u - 10186
~
in
zL = 0.926 in
ze s 0.867 in
Radius of effective section 8
e = 108.926 in.
1n terest ing
Ratio I (effective section) ='0.536 g 1.35
I (rib on11) 0.397
All values are now available to solve the two equations (18)
and (20) for the two redundants He and Me. The solution 1s made
for the e~sumpt1on of a rigid" restraint (condi ti ons (10) ) and an
elastic ,restraint with meesured displecement and rotation of the
are found
a 0) Rigid restra 1 nt: H
c = -1068 1b
hI
c ::
3156 + in - 1b
Mc = 4371 + in - Ib
The calculation of the direct forces Ntf and the cross bending
moment Mx iI( the shell fo llowed the procedure briefly outlined on
po 13. 'Ihe r,esu It s are pre sen ted in Fig. 13.
two level bars chflcked rotations (accuracy 0.0002 rad.). The actual
load! ng co~s iated of 10 equally spaced horizontal single loads as
,
between the theoretical. loading, assumed uniformly distributed, end
the actual loading can cause only local differences in behevior o
I
I
The stresses in the ribs ccmputed from the strain gage 'readings
I
ere shown 'in Fig: 12. Note the agreement between test and theory
if the analysis is made with arbitrary allowance for the elastic
80) Conclusions:
thrust I1t
and the em moment M • Introducing the virtual lQsd H'k
k
:: 1 (Fig. 16b) and the virtual moment M'k 1 (Fig. 16b), the hori- =
zontal disPlacementdk end the end rotatlcn ¢k respectively can be
calculeted by means of tl:e work equation.
The actusl normal force N and bending moment M are (see Fig.
16s) :
N ~
- Hc C08 f ( 22)
Due to H' .- I
1
(F lg. 16b)
k -
.N' :- - cos 'f (24)
,"p
-... -Se (cos 1- COSfk) ( 25)
""
Sk - M'Mds+
EI
W'Nds
FA
Taking Mt k l: 1 an d
Nt =0 (27 )
Mt =- 1 (28)
B. ¢k
ee
~Mk (':!k + i) + Hk8e (t.fk cOB)'k, - Sb~k) (29),
Where: K = g lC
ae
•
S-Forces ere calcula ted by "Eq. (4), once the stress (jA a long
the conre cting line of rib am shell is known (Eq. (21) ).
Knowing the boundary conditiC'ns for the shil, which are
,
Rib x - 0 l ~x =0
Rib x : 22:),
Rib x =1. ~
OX
= 0,
*See p. 12 for further explanation.
2l3K 21
dis pIa cen;tent dk a1?-d, end rotation ¢k. The effective width of
the ribs must be ca1cula ted· first. In order to know the number
n ::.
005866
ra -... 0.14
n ~i
The coeffic lentIL for the middle and outer ribs ( / being
(2) Out er
rib:jlf: 0] K--
0 .. 38
k: 0
. b :: K' feh l:J 1.,36 in
Distances Z u ~ s -1.382 in
. zL e O~73l in
zA = 0.672 in
Radius of. effect:1.ve section
4
Moment of inertia I :: 0.921 in
Ratio I (effective section): 0 920 :: 1084 0
I ( r1 b only) 0.491
Outer Ribs~ (as on p. 16)
Area 1022'7 in 2
Distances ·z ~a -1 .. 186 in
t(~ \ .
zt= 00926 in
zA §
00867 in
Radius of effective section 8
e - 108 .. 926 in
MornEll t of inertie I -= 0 536 in 4
0
dk ~ =003300 in
'" c
~k 70'753xlO-4 ra di ans
Taking the coefficient of elastic restraint
2l3K 23
Ie:: 0005 8e
EI
and inserting in Eqo (26) end (29) the proper values, the
Mk • -1'7850 in-lb
Ou tar Rib: Hk :. -n46 Ib
~k -10350 in-lb
~
6 6) Con c lu s ions:
CII
-r--
+ c2 I z
I
= -c l rI -+ c
2 z (31 )
exactly the hori zontal d1 sple cern en t Jk alr-eedy tree ted under'
III. Be replacing in Eqo (26) end (29).
](32)
the horizontal thrust Hk and the end moment Mk due to s temp-
erature r1 se L1 to are thus determined.
213K 26
been developed.
The ribs, reinforcing the shell at regular intervals l ,
increase the buckling stiffness of the shell panel in between.
Furthermore, they are l~dispensible for carrying a one-sided
live load (snow), acting on the structure. To provide B
Area Ai = 1.,9608
Momol .1 1 =0 0 9204
essentially as predicted.
/" ~0;~~~ ~~~~:;~,.
Mariy tests not reported herein were made during the course
..
-
of the investigation but, space does not permit their lnclusiono ' , .
Similar procedures of analye1sl' based on the use of ·the "effective
",
1'or e.ach of. the s 19n1ticant terms or the J:i'ourier Series, and the actual
stress at any point could be determined by superpositiono
------------
',.' ..
;: .
. , .
.~: :;.'~.: -.'t"''!'t.-., .;.
"
.''1.
~ .....,.,
~.:"~ . . ""_, ;._.-. __
~i.::.
"~
,,,.,
ACl\NOWLEDGME'N T
. .
O. Gruenwald, W. A. Renner and P. Rongved, all of Roberts end
Schaefer Company, during several meetings are sincerely 8cknow-
ledged.
LIS'l' OF" REFERENCES
.I I I
I I: !! I I ,!,
I' : I I II;
T- Beam Fig. 2
Cylinder Ring
------ ---'--
- - -I-- - - -
f----------
- X ----110 f - I- - - X
Fig. 3
1.40
;--In finitely long cylinde'r
I. 20' I-----~
.60 f3l=O.2
Semi-infinite
.20 I
I. Effective width b = KM I
0 .....- -........- - -........- - - - - - - -.....- - - -
----.. A
0. .50 1.00 1.50 2.00 2.50
l
--f~-}~--:--I of..
(3 =
= n~
1]1
, I--I--l S = S"cos n<P
, Fig. 4
0) Theoretical LoadH19
. p (lbjin)
t ttttttt ~ t t t tt~ ~!~ t tt t ~ tt t tt+t tttt
I ! ~
I I
-- :.....-..J
-Fig.5-
-Fig 6 -
• • '1--
I .
c.) ~----
I /I
I / I
II /
°e ,---
I
M'= I I
I t "....\ . ) P~1
~--f-- --+---~
r<.; J i
-4--
Shell Shell
IEffect,ve Effective
W'dt'h Width
R,b Rib
-Fig,g
SHELL ROOF MODEL
.• - • - <-.
ALL DIMENSIONS IN
INCHES
I I
, I
59.78 59.78
C
B
12.0 12.0·
2.0
j~O.625 ~1t-O.505
SECTION a-a
11
1"=4
0.20
1.0 0.1\8
~--r--r--
2.5 3.5
SECTION A-A
3.5 3.5
\.0
1.5 1.5
10.0
3.0 .~-~------=~-------+-------,.------.
12.0 12.0
I.. 3.0 '
..I
FIG.IO·
a.) Outer Ri b b.) Middle Rib
0.505"
Effective section of
the ri bs
-10000 1/ ,,_- 1--/ I '-: --v/· \
-15 0 OOt--l~'----t---+---+---+--__+--t-\----+-1- - + t - - t - : - - - (~--I---I---+---+----f---I----+--+---+---'H"
Edge -member:.-r-~ ,
1-+-\-+---+---+--...----,--r------ 1
----4
1
~==I===,=F===;I;:::===;to=~-~d3!!::;====1 ~-I I
5000-~
~
( 0 - (0 _
en (\J
tow -100 _ ~_ ~-
~ 00
en -
r<)
en_en __ ~ -+--+-- ~
ro to to
-- en -
to
en -
ro
en -
-
0 0 _ 00_
(J) I'-
ro- I'- -
ro
(0 -
-
~ (0
I'-~
-
(/)
o ~,.,'
~ ~X q..
'
~I
~I
~~ (\1
I'
~.
'"
='
q
,
+-.
~
~ = I'- ~ ~,
~
.
.
~
",'
~.
o'
~• x
~ LQ
FIBER-STRESS
<fsMf>
AT UPPER
AND LOWER
...
EDGES OF
..
. RIBS.
FIG. 12
12" ------------1
MIDDLE RIB
/ / /
<t :ep=O o
o ~~~-:------,----.,......v----:::::O--X-===+===+==./=+>-=--t'J:::::;~::::;;;;;;;:~I=X=~-X-,----JO
\ 0·
~--t---t---/-+V~--+---t----t--t---t---+--4-
-400f-\-\--+-V-. 4
/' f'/'
-800 X ~ y - -8
V~, 0 ... / '
"- -r----
-1200 I-------'--.. -~ -.-._-.-._---'-_---'-_---'-_---'-_----l-_-'-_ t (62) -
o
-12
.. 400
\ /V
/"-t---+---t--'---J7-!'--t---t---t---t----+---+--1 -4
~\-t-I-7"
,//
)<~ y/
_.
- 8 00 I-/---j~--+---+---:.-"./~---+----+---+----+--_+__--+--+-~
-8
/ "......... ---../
- -- . c:
"-..c
-"
c:
-1200 I------'------L- -12
-
I
c:
.-
--e-
..c
z
o
~,,~~
~
Col. 2: ep:0.1739
0 0 _--
~--t--I
0
·1)( I Ix
0
I
0
x
a ~
.c
.....
z
(/) -400 ~--~-- FORCES Nep8 MOMENT Mx in SHELL 1-----1 -4 w
w ~
'"
0::
o
lL.
- 800 I---+---+---+--_+__
THEORETICAL Nep
- - - - THEORETICAL Mx -8
0
~
I"-
...... --
o
-.....:
EXPERIMENTAL Mxt - - - 1 8
I
I 4
a
z
w
al
x
o v/. "'------ ..... - - -
1 1--_
a
Col: 3: ep = 02898
0
~ (-6.4) 16
I 12
~~I
I
0;-a-
1200 0 1'-',,: I
1
12
,
1
r- ' " '
80 a / I
0 "I :
8
/
/ iI '"' " ;'"
I "-
3: 40 A ! I I " "- 4
0 Y .~.
3: ." I x
i~
0.......... " 0
0
0 ll "'-, x
"TI
",
Z
::0
0 a Co I. 5:
.J..
'fI =0.4637
- )( r- a
-t ",
Gl (J)
3:
)(
1(-10.1)
0' Z
Z -e-
.. eft l>
:I: Z
ITI 0
r
.r
Shell-;::=========i=X=========r;.]2L
Edge-member ,
o ad¢
I Theoretical Deflection
x=o~"--
I
at center (E =29.0 x 106 )
. ./" ~"
~/ ~
x=21.~_
/d""': ,.,.,._:..---'
~
. - - __ + __---- ~
/'
h- /:--"
~.
. -".l'-""'--
FIG. 14
Actual 0Y5tern:
. a)
d K
= horizontal . dtspl~cement
1K :; end rotation
.-
.!
~. to •
30000 '
x
,---- ------ -------- ---- - - --- - ----- ,---- -- - - --- - - - - - - ... - - -------
-,
1 ----- 1= - --J
25000 1-1------+------ ------;-t----------e--- --- GL Theoretical stress in lower fiber ---i------------- -1-- ----------~----_H
X Experimental stress in upper fiber} - x
20000
r\.
I ''\.. :~ f-----c-----1--r 0 Experimental stress in lower f'lber. Outer, rib x = 0" and middle rib --.v-
_ . __ _
'\ I: . 6 Experimentat stress in upper fiber} Outer rib x= 21,
15000 ~ Experimental stress in lowerfiber--rr-----j----+--~/'_t_--__t__l
- 1___ ~ i---if- - 0--_
1-1----- •
I ~__ J /
I --
5000
a
--------··1
"" 'v///'-"--/
'~//--~~TER RIBS
I
t
i ------ " '/
-',~__/'-/---+---I--------j-
-5000
-10000
(/)
0-
c: -15000
w
w
Cf)
Cf) 30000 fR
0
• w I
0::
I- 25000
Cf)
20000
15000
10000
5000
-5000
-10000
-15000
12" -----------~.~
1500
1000 ~ x ~
20
k* --
'\
I--
-- /
-. :><v*./
0
0 1-8-_-' - - 0
500 " 10
~
/' 0 o
~
8
'II
o ""- *
I--. ~ I ''""~
0
o
t. ep=O \,
-500 -10
\0
0
1\ -20
\
\ -30
\
\
-40
"-
~
x
V- m
~- ---1 - _ - v
-. ~
/ JO
fTI
z
~
......... ~ t--_~ 0
// ~--- ~-o-
v/ 0
i'z t'-- ~ ill'....
'", 0
Z'
G>
3:
4> = 0.1739 " 0
'", , -10 . 3:
m
z
.;
° \ -20
CJ)
\
~
THEORETICAL 3:
\ )(
- - - - - THEORETICAL \
X -30
x
0
EXPERIMENTAL Nep
EXPERIMENTAL M x
-~
•
CT
........
~
,...- 10
-S'
_ ....
10
z
-t
CJ)
~
fTI
CJ)
0 --
lL-----
ep =0.4057
+0(5.53
)(
~ .1-- .- --- I- l- .- 3- .-
-- le;
-L
~
- ,, .
o
3:
)(
» -500 -10
z0
:z
CJ) 20
I /1
fTI
r /
r / 10
ep =0.4637
0 ",
+0 (7.20~
l-- x.
_"3: ~ l-o-
v/ " o
~ ---
"-
cO
(»
-500
V
..........
./
V
X
f--.-
-- 1-0-- '>---
Ix--r---
>I. -10
-1000 -20
J