Behaviour of Thin-Walled Structures: 14, Spence), Le
Behaviour of Thin-Walled Structures: 14, Spence), Le
f ue ateParatty ig ;
rH ti PPA abtetat ine teaattaesy |4 {
} 14 sh i i As? t i +
Behaviour of
Thin-Walled
Structures
| | Edited by |
| nN picchmed
14, SPENCE}, Le
https://archive.org/details/behaviourofthinw0000unse
BEHAVIOUR OF
THIN-WALLED STRUCTURES
FRONTISPIECE
James McQuistan Harvey obtained his PhD degree from the University of Glasgow
in 1953 for research in the field of cold formed steel sections. After five years spent in
the cold formed steel industry, he returned to academic life in 1958 as a senior
lecturer at the University of Strathclyde (then the Royal College of Science and
Technology), Glasgow. He later became Reader in Mechanics of Materials and in
1969 was appointed the first Trades House of Glasgow Professor of Mechanics of
Materials. He occupied this position until his early retirement in September 1982
and for the last five years was also Deputy Principal of the University of
Strathclyde.
Throughout his academic life Professor Harvey was involved in thin-walled
structures research and played a significant role in the area of cold formed sections
at both national and international level. In the light of this involvement it is fitting
that his retirement should be marked by a Conference devoted to the behaviour of
thin-walled structures.
BEHAVIOUR OF
THIN-WALLED
STRUCTURES
Edited by
J. RHODES
and
J. SPENCE
Department of Mechanics of Materials,
University of Strathclyde, Glasgow, UK
ISBN 0-85334-246-6
All rights reserved. No part of this publication may be reproduced, stored in a retrieval
system, or transmitted in any form or by any means, electronic, mechanical, photocopying,
recording, or otherwise, without the prior written permission of the copyright owner,
Elsevier Applied Science Publishers Ltd, Ripple Road, Barking, Essex, England
This volume contains the papers presented at the Retiral Conference for
Professor James M. Harvey, held at the University of Strathclyde on 29-30
March 1983 and organised by the Department of Mechanics of Materials.
Attendance at the Conference was by invitation and was limited to about 50
delegates, all specialists in the thin-walled structures field.
The Conference was sponsored by:
subtleties were not fully realised have in these years been subjected to
intensive study. As a result, great advances have been made in the areas of
elasto-plastic buckling and post-buckling behaviour, interaction of local
buckling with other forms of buckling, and the effects of imperfections on
the behaviour of plate and shell structures.
Also, during this time the trend towards lighter, thinner structures has
led to the desire to use extremely light gauge metal, for example in roof
sheeting, and has necessitated a substantial amount of research and
appraisal of various aspects of design at loads far beyond local buckling.
The growth in use oflightweight strong materials, such as fibre-reinforced
plastics, has also been a contributory factor towards this need for advances
in the knowledge of the far post-buckling range.
New problems in the field of buckling have also been posed with the
emergence of offshore oil platforms. Problem sources have arisen both in
the platform structures and in the associated subsea pipelines. This has
resulted in significant research input in the area of buckling and collapse of
stiffened cylinders, and on the various problems of pipeline buckling and
buckle propagation in pipelines.
The present papers cover some aspects of all of the above problems.
Twenty papers deal with the general topics of: compression members,
column behaviour, cold formed sections, lightweight structures, shells and
pressure vessels and plates. Of these, nine are devoted to new applied
research within these topics, five present surveys of theoretical, experimen-
tal and industrial aspects of past research and six are concerned with design
implications.
We hope that readers will find the papers interesting and informative.
J. RHODES
J. SPENCE
Contents
MedCLIC AeWe Hee aD c 5 AC, Meets ecb ates ae oss, Sod I ats
P. S. BULSON
Military Vehicles and Engineering Establishment,
Christchurch, Dorset,
and University of Southampton, UK
ABSTRACT
The history of research and analysis to establish the behaviour of plates and
sections reinforced by bulb and lip stiffeners along free edges is reviewed.
Simple design formulae are given to establish the elastic critical stress in
uniform compression. The treatment of ultimate compressive strength is
discussed with particular reference to tests on aluminium alloy extrusions,
and the paper ends with random thoughts on the research/design interface.
1. INTRODUCTION
2. ELASTIC BUCKLING
by Casilev’ in 1914, but his work was not summarised or made widely
available to English-speaking structural engineers until the publication of
Timoshenko’s monograph? on the theory of elastic stability in 1936. It was
also in this book that the results of analytical studies by A.J. Miles of the
University of Michigan were presented in graphical form, and these gave a
first insight into the fundamentals of the problem to many a designer and
research student. Timoshenko also reported that the problem had been
discussed in the early 1930s by Melan? and Rendulic.*
The designers of aircraft had begun to use edge reinforced sections
before this theoretical work was reported, but behaviour in compression
was almost entirely judged from the results of large numbers of strut tests.
Thus we find Sechler and Dunn? in their 1942 book on airplane structural
analysis and design stating that ‘for sections having bulbs no rational
methods of analysis are available, and recourse must be had to testing
enough lengths of the section to determine the coefficients in the short
column equation’. They suggested that design methods which assumed that
the bulb supplied simple support to the edge of the sheet should be used
with caution.
When rational theories became available in the 1950s and 1960s, it was
shown that thin plates in compression having one longitudinal edge simply
supported and the other reinforced by a symmetrically placed stiffener
equal in thickness to the plate behaved as indicated in Fig. |. The critical
stress in end compression, o,,, 1s given by the familiar equation
2 kn?E Ne
Per 19(1 — v2) \b (1)
where k is the buckling coefficient, ¢ is plate thickness, b is plate width, and
E and vy are the elastic modulus and Poisson’s ratio respectively. If the
length of the plate =a, and the width of the reinforcement = b,, Fig. 1
shows the relationship between k and a/b for various values of b, /b. The
analysis assumes that when buckling occurs the simply supported
longitudinal edge remains straight, whereas the edge between the stiffener
and the plate bends outwards. However, if b, is made wide enough, the
section is transformed into three thin component plates, and at buckling
the junction between the plates will remain sensibly straight. In fact, there is
a gradual transition from one mode to the other as b, grows wider, but let
us for the moment concentrate on the outward bending mode, associated
with lips having a relatively small cross-sectional area compared to the
component plate.
4 P. S. Bulson
a/b
Fig. 1. Symmetric lips on hinged flanges.
(a) Long plates simply supported on one longitudinal edge buckle into
only a single longitudinal half wave. This is true whether the other
edge is free, or whether it is supported by a small lip.
(b) As the ratio a/b increases, the value of & for a lipped plate falls
below that for a plate with a free edge. At the value of a/b where the
buckling coefficients are equal, the critical stress of the lip as an
isolated strut constrained to bend in a plane normal to the plate is
equal to that of a plate with a free edge.
If the plate is fixed along one longitudinal edge, and reinforced by a lip
along the other, the picture changes somewhat, and a typical solution for a
lip equal in thickness to the plate is shown in Fig. 2. The curves linking k
and a/b now take the familiar ‘garland’ form, with &k reaching a minimum
value (k,,:,) aS a/b increases. We can make the following observations:
(c) Adding asmalllipcancause k,,;, to fall slightly below its value for a
plate fixed along one edge, free along the other.
Bulbs, Lips and Beads 5
(d) As the width of the lip is increased, k,,,,, approaches its value for a
plate fixed along one edge, simply supported along the other. When
this value is reached the buckling mode changes to several
longitudinal half waves.
In the calculations leading to Figs. 1 and 2 the torsional stiffness ofthe lip
was neglected, and the error in so doing is small. The author®’’ has shown
that it is about 5%. A close approximation to the value of k,,;,, in Fig. 2 can
be obtained by employing the energy method of solution, and examples are
given in ref. 6. The parameters that influence the value of k,,;, for plates
fixed along one edge are EJ/bD and A/bt, where E/ is the flexural rigidity of
the lip bending in its own plane, D is the unit flexural stiffness of the plate,
and A is the cross-sectional area of the lip. The variation ofk,,;, for a range
of values of EJ/bD and A/bt is shown in Fig. 3.
The critical stress of a uniform thickness thin-walled channel or I section
can be expressed in terms of a buckling coefficient applied to the web plate.
6 P. S. Bulson
EI/bD
Fig. 3. Variation of k,,;, with area and stiffness parameters.
byt s 50
0 02 0.4 io ee 10
b,/y
Fig. 4. Buckling coefficients for lipped channel sections of uniform thickness.
Bulbs, Lips and Beads 7
Thus,
wE fave
Gor = (nin) 12(1 p= vy?) (5) (2)
where b,, is the width ofthe web and f its thickness. If the flange width = b,,
it can be shown that for a channel of uniform thickness (k,,;,,), Varies with
b,/b,, according to the dashed curve in Fig. 4. If a symmetric lip, equal in
thickness to the flange, is added at the extremity of each flange, the value of
(K min)w Will increase by a factor («) that is governed by the ratio 5, /b,,. The
exact solution for a lipped channel is also shown in Fig. 4, which is drawn
for a section with 6, /t = 50.
In the range 0 < b, /b,, <0.25, the lip factor («) can be expressed quite
accurately by the approximate relationship
b 2
a=1+40 (Fe) (3)
b,
as indicated in Fig. 5. For lipped channel sections with 6,t= 20, the factor
is
b 2
=1+8(—
a=1+ (5) 4
(4)
iS
ee a ==
(5)
60 \ t
so that, in general,
n?E pave
O oy
he)
( ee
A. Ot 12(1 =) (+)
6
( )
applicable in the range 0 < b, /b, <0.25, 20 < b,,/t < 50.
So far, the discussion has been limited to symmetrical lips, but most
practical sections, whether rolled, folded or extruded, have the reinforce-
ment located on one side of the plate, to form inward or outward turning
lips. A method ofdealing with this was outlined by the author in ref. 8, and
shown to give reasonable agreement with tests by Dwight,” but a simpler
and equally accurate procedure is to replace Jin the expression EJ/bD by an
effective value equal to double the moment ofinertia of the symmetrical lip.
8 P. S. Bulson
3-0
[Keir ty
ZO
ee)
of k,;, can be accurately computed, and the results for a plate having
b/t = 20 are given in Fig. 7. The effect of neglecting torsion is shown by the
dashed curve, and this indicates that at higher values of d/b, k,...min would be
reduced by almost 50%.
Thin-walled channel and I sections reinforced by heavy bulbs have been
examined in refs.6 and 7, and the results for a bulbed channel with
b,,/t = 20 are given in Fig. 8. For uniform sections with 20 < b,/t < 50 and
0 <d/t <6, it is approximately true to use the relationship
ipa
Lac
=1+—|- 9
" 20 (“) 0)
ed _—_—
_—
—
=
<i ——- _
NEGLECTING TORSIONAL
STIFENESS OF BULB
nee ol A/b Or 03
(+)
Fig. 7. Buckling coefficient for bulbed flanges.
Bulbs, Lips and Beads 11
b,/t = 20
be/b
Fo -W
Fig. 8. Buckling coefficients for bulbed channel sections of uniform thickness.
usually in the form of graphs linking the buckling coefficient for one of the
component plates in terms of geometrical ratios of plate widths, thicknesses
and lip dimensions. However, the writers of codes of practice for structures
other than aircraft have, since those days, always considered that the
graphs were unnecessarily complex, particularly for initial design
calculations. At best, a few were relegated to the appendices to interest
those who wished to make academically accurate calculations of critical
stress.
A slight simplification was proposed by Rhodes and Harvey'! in the
early 1970s, who devised critical loads by considering the stability of
component plates with rotational restraint along their longitudinal edges,
and matching these restraints at the common junctions. However, as
Walker’? has pointed out in a book that indicated the sources ofthe design
recommendations for cold formed steel sections in building (British
Standard 449 Addendum No. 1) 1974, the calculation of o,,, and hence
12 P. S. Bulson
2
_ 12(1—v?)
ay? (10)
Using this approach, the design clause for a thin flange reinforced with a lip
of the same thickness as the flange reads as follows:
The value of mb/t to be used is the largest of those obtained by
The expressions for the increase in buckling coefficient to allow for the
flange—lip combination were similar to those given in eqns (7) and (9), but
they were adjusted to allow for the fact that m decreases as k increases, in
proportion to syle. They were also slightly more conservative than eqns (7)
and (9). Further, in order to get a meaningful set of design rules,
it was assumed that the minimum critical stress of a plate simply
supported along one edge, lipped along the other, would increase as the lip
was increased in size, whereas we know from Fig. 1 that the buckling
coefficient for long hinged plates does not increase when the plate is
reinforced by an edge stiffener. The ‘lip factor’ [1 — (b?7/80t7)], which
should strictly be applied to values of m corresponding to plates with edge
restraint along the hinged edge, has been applied to m = 5.1, which is the
value for a plate simply supported along this edge, free on the other. This
notion allows the idea of ‘web’ elements and ‘flange’ elements to be used,
and at the same time takes account of the increase in elastic buckling stress
of the section due to the presence of a lip.
In the revision of CP 118, taking place at the time of writing, there is a
possibility that the critical stress of web or flange elements will be calculated
on the basis of an increase in the plate buckling coefficient, k, using eqns (7)
and (9). The precise nature of the design clauses is still under discussion.
For the flange—lip combination, it is interesting to judge the value of 5, /t
that will produce a buckling coefficient value for the flange equal to that for
a simply supported plate. As an example, consider a channel with equal
flange and web widths. The value of A, is 0.9, and the flange-lip
combination to produce a ‘stiffened’ flange is given by
2 b
0.9 make fu =4 so that= 11.7
40 \ ¢t t
For a flange with b,,/t = 50, the value of 5, /b, is just over 1/5, which agrees
with the figure used in the design of cold formed steel sections. The required
value, however, varies with the 5,./t ratio of the flange, and with b,/b,, for
the section.
386
COMPKESSIVE
MAXIMUM
(N/mm?)
STRECS
443
462
43)
400
370
(N
COMPLESSIVE
(mm?)
MAXIMUM
STRESS
1/3
(N/mm)
STRESS
FAILING
re) 30 6O Le)
SLENDERNESS RATIO
Fig. 11. Strut tests on a lipped I section in high strength aluminium alloy.
when this is influenced by the collapse of the lip in a torsional mode. The
reduction in axial stiffness is between 55 % and 65 %, depending on lip size,
and the reduction in flexural stiffness of the section is very severe—over
90 %. This has a significant effect on the column buckling stress of members
containing such elements. With regard to lipped channels having
asymmetrically stiffened flanges, Rhodes suggests that there is a high
sensitivity to local imperfections, particularly of the lip. He applies his
theory to lipped channels with 5, /t = 50 and symmetrical reinforcement
and shows, for example, that with a geometry given by b,/b, =0.6,
b, /b, =0.4, the reduction in axial stiffness at buckling is about 30%.
The distribution of compressive stresses around a lipped channel
immediately after buckling is given in ref. 18. The non-uniformity results in
a shift of the neutral axis in bending and a decrease in column buckling
load. The ultimate load relationship obtained in this way agrees reasonably
well with experiments conducted by Harvey’? of similar vintage to the
author’s tests on aluminium sections.
For sections having lower values of 6,,/t, the loss of buckling strength is
much less, and in a series of strut tests on lipped I sections, reported by the
author in ref. 13, the use of the maximum local buckling stress as the upper
18 P. S. Bulson
4. DISCUSSION
This paper has not reviewed the problem of lipped or bulbed flanges
forming the compressive material of thin-walled beams subjected to
bending stress; neither has it explored the behaviour of large thin-walled
sections with edge reinforcements welded into position. Both of these
aspects are frequently met in design, and are important subjects.
Nevertheless, within the limited scope of uniform end compression of
folded, rolled or extruded sections, there are enough problems left to keep
research workers locked in discussion for some years.
In looking back three decades or so, the author has been aware of
features of the research/design interface that perhaps have a wider
background. The first is the timeless value of well conducted experiments.
Theories come and go and treatments in codes of practice are subject to
cycles of fashion, but good test results are indestructible. It is fitting, in view
of the purpose of this conference, that the penultimate of the references
quotes the PhD thesis of Professor Harvey, from 31 years ago; and perhaps
it is more than coincidence that the tests in support of the treatment of
lipped aluminium flanges come from the author’s postgraduate work ofthe
same period; and those in support of lipped cold rolled steel channels
quoted in ref. 12 were first recorded in the 1950 PhD thesis of Sir Henry
Chilver. Great emphasis was placed in those days on the importance of
productive experimental work as a vital part of engineering postgraduate
experience. In one of the research investigations mentioned above the
student machined, prepared and tested 377 specimens. There was no
productivity bonus, but a lifetime grasp of the fundamentals and difficulties
of structural testing was guaranteed.
The treatment of bulbs, lips and beads in design codes is a representative
illustration of the difficulties faced by code committees. If the author’s own
experience may be quoted, the static design of thin-walled, lipped
aluminium alloy sections for general structural purposes was carried out in
the 1950s and 1960s by using ref. 2, a collection of aeronautical data sheets,
a text book or two on aircraft structures, and a few reports issued by the
Royal Aircraft Establishment in this country and by the National Agencies
in the USA. As everything was written in a different notation and in
contrasting styles, the designer had to come to grips with the fundamentals
Bulbs, Lips and Beads 19
to be sure that he understood the problem and could evolve his own design
procedures. Moreover, he was using the ‘allowable stresses’ design method
and limiting his needs to a knowledge of elastic critical stress and plasticity
reduction factors. Over the years he became very familiar with the plate
buckling coefficient, effective widths, and formulae like eqn (11).
In 1969 the first Code of Practice for the structural use of aluminium was
published by the British Standards Institution. It was a splendid attempt at
codification, but written in a new language. In the local buckling sections
there was no mention of the conventional plate buckling coefficient, k,
because the equivalent slenderness ratio was used to give a concise method
of using a small family of design curves for local, lateral, torsional and
column buckling. It was simple to use, but not easy for the designer to
retain the grip on fundamentals that he had developed in the preceding
years. The effect of lips and bulbs on buckling strength could be easily
calculated, but without thoroughly understanding the principles of
interaction between a plate element and its edge stiffener.
The 1969 Code will probably be replaced by another in the mid-1980s. It
will no doubt be a most useful document, but the language will have
changed again. The policy now is to use ‘limit state’ design and the
reinforcing lips and bulbs will be attached to ‘compact’, ‘semi-compact’ or
‘non-compact’ sections. The old-fashioned ‘plate buckling coefficient’ is
likely to be hidden away again in new, concise methods of summarising the
relationship between elastic buckling and ultimate strength. Since research
leaders in civil and structural engineering departments of the universities
are still rightly content to use the traditional approach to buckling
problems, as illustrated in refs. 16 and 18, designers find themselves having
to master several approaches to the same subject when they wish to practice
according to the code and at the same time keep abreast of progressive
research. In searching for simplicity we seem to complicate the divergence.
This results in those who can find time having to add to the strain of their
lives by writing books and articles in an attempt to bridge the gap.
The third feature of the research/design interface illustrated by the
subject of this paper was touched upon by Professor Harvey himself in
writing the introduction to the proceedings of a conference held at the
University of Strathclyde four years ago.”° He pointed out that work used
in design recommendations is like the visible part of an iceberg and that
designers should be aware of the vast unseen substructure of fundamental
work beneath. This is important because the lack of a clear understanding
of the physics of structural behaviour can lead to accidents. A quick review
of published work on lipped sections suggests that the properties of his
20 P. S. Bulson
iceberg are about 500 pages of research and analysis to 16 lines in a design
code. In wishing Professor Harvey a long and rewarding retirement we
naturally wonder whether he will find a similar relationship between
practice swings and match-winning iron shots.
REFERENCES
ABSTRACT
NOTATION
1. INTRODUCTION
peeenagy
ON
———
ES FU
f 1
(1)
is less than the 0.8 limit value as adopted by ECSC in its
Recommendations. In this definition, the width b,, of the profile is
measured between the median planes of the opposite walls and 1¢
represents the wall thickness; the yield stress is f, and E is Young’s
modulus.
JA\E
/is the column buckling length and A the cross-sectional area. The value
of parameter y depends on the manufacturing process of the section:
(a) hot finished profiles: Ve So 55
(b) cold finished profiles: y=1.014.
24 R. Maquoi, C. Massonnet
and J. Rondal
= VFy/ocry = (b,/1.9t) y N
Sarl He (N=
le Geen anak ]
= ot heen. = (Wai) vf JE
Manufacturing process
Notation:
== (4) =y (3)
Such a technique has produced sets of optimal ranges (Figs. 2-5). Figures 2
and 3 are devoted to hot finished profiles the yield stresses of which are
equal to 240 and 360 N/mm/ respectively, while Figs. 4 and 5 concern cold
26 R. Maquoi, C. Massonnet and J. Rondal
t(mm)
30
=10m
25
8m
6m
20
1SO_657/XIV
AON
10 Begin
b(mm)
finished profiles with similar values of the yield stress. Calculations leading
to the aforementioned figures show that:
t(mm)
t (mm)
30
E
S
25 E a
E
wo
20
i: os
5S
~—wt
D
15 E =
a
N Oo
vp)
/
be)
=| / /| ioe
| WYa
EEN vise
b(mm)
The ISO range seems rather well calibrated, although the thicknesses
seem to be a little too large for small sections.
The BSC range is less well calibrated; it uses excessively large thickness
values for small sections; for large sections, it does not propose a
sufficiently wide fan of thicknesses, especially towards the large
thicknesses.
The Use of Thin-Walled Hollow Rectangular Sections 29
t(mm)
30
E
25 L is)
Bh ae
i
ao
€
wo
20
S a
a
g
dy
155 Ae =
(je)
“
l= 2-4-6-8-10m b (mm)
0
100 200 300 400
Figures 4 and 5 (cold finished sections) show that the ISO range is here less
well calibrated than for hot finished sections; indeed, the sections are too
thin for the large widths.
One should keep in mind that, in the range of the small sections, various
problems—such as the requirements for the connections, the technological
conditions due to fabrication and corrosion—require that the walls should
not be too thin. This justifies the manufacturing of sections whose wall
thickness exceeds the optimal value computed in axial compression.
30 R. Maquoi, C. Massonnet and J. Rondal
It has been stressed above that, in the vicinity of the optimal thickness,
the change of efficiency on the side of the large thicknesses is substantially
less than on the side of the small thicknesses. This remark enables a simple
extension of the results obtained for square hollow sections to the case of
rectangular ones. One has only to relate the choice of the optimal thickness
to the width of the largest wall and use directly the results obtained for
square tubes. This procedure yields only a very small error because, on the
one hand, the area of the narrow walls is always less than that of the wide
ones and, on the other hand, the efficiency of the narrow walls is only
slightly changed.
t (mm)
30
25
20
657/X1V
1S0
15
b(mm)
From the consideration of Figs. 2—5, it is found that, for the common steel
grades of hollow sections, the thin-walled sections have a rather marginal
position in the optimal ranges. However, Figs. 6 and 7, which are drawn for
a yield stress of 700 N/mm7?, show that the situation becomes quite different
for higher steel grades.
Such a high value of the yield stress could look strange; it must, however,
be stressed that, in Japan, it is presently planned to fabricate welded boxes
t(mm)
30
E
ro)
u
25
20
x
15 ,
10+
Y {=2-4-6-8-10m
VE: 3b(mm)
made from rolled plates whose yield stress is nearly 700 N/mm? or even
greater.°
Progress recently achieved in steel metallurgy should soon enable one to
put on the market, with excellent economic conditions, plates with very
high yield stress that could be used in the fabrication of cold finished tubes.
For hot finished tubes, the situation does not favour such an important
increase of the yield stress; nevertheless, some manufacturers are already
able to roll H profiles the yield stress of which exceeds 500 N/mm?.
Such a viewpoint is also supported by the fact that the detrimental
influence of the residual stresses on the ultimate strength capacity in axial
compression decreases relatively when the yield stress increases. This is due
to the fact that the intensity of residual stresses increases less than in
proportion to the yield stress. Proposals in this way have already been
presented by one of the authors’ but would require a more extended
experimental confirmation.
The authors are convinced that the use of high steel grades may lead to
substantial savings, even when the buckling phenomenon is governing.
4. CONCLUSIONS
For the steel grades that are presently used in the fabrication of hollow
profiles, the thin-walled sections occupy only a marginal place in the
optimal range established in axial compression. Various factors will
apparently lead, in the near future, to an increase of the yield stress of rolled
plates, which constitute the basic material for hollow sections. As a
corollary, an increased use of thin-walled sections may be predicted.
REFERENGES
1. BRAHAM, M., RONDAL, J. and MAssonnetT, C., ‘Large size buckling tests on steel
columns with thin-walled rectangular hollow sections. Comparison with design
methods’, Thin-walled structures, J. Rhodes and A.C. Walker (eds), London,
Granada, 1980, pp. 3-13.
2. BRAHAM, M., GRIMAULT, J. P. and RONDAL, J., “Flambement des profils creux a
parois minces’, Commission des Communautés Européennes, rapport eur 6730 fr,
décembre 1979.
3. BRAHAM, M., GRIMAULT, J. P., MASSONNET, C., Mouty, J. and RoNDAL, J.,
‘Buckling of thin-walled hollow sections. Cases of axially loaded rectangular
sections’, Acier-Stahi-Steel, No. 1, 1980, pp. 30-6.-
The Use of Thin-Walled Hollow Rectangular Sections 33
. RONDAL, J. and MaAquol, R., ‘On the optimum design of square hollow
compression members’, I[UTAM Symposium on Collapse, University College
London, 31 August—3 September 1982.
. RONDAL, J. and MAQuotl, R., Etude d’une gamme optimale de profils creux carreés
et rectangulaires, Annales de Il’Institut Technique de Batiment et des Travaux
Publics, No. 409, novembre 1982.
. Usami, T. and Fukumoto, Y., Local and overall buckling of welded box
columns, Journal of the Structural Division, ASCE, 108, No. ST 3 (March 1982),
525-42.
. Maaquol, R., ‘Some improvements to the buckling design of centrally loaded
columns’, Structural Stability Research Council, Annual Technical Session and
Meeting, New Orleans, 30-31 March, 1982.
¥
tere
Amplitude Modulation of Short-Wave Buckling
Modes
W. T. KOITER
Delft University of Technology, The Netherlands
ABSTRACT
1. INTRODUCTION
localized dimple. The result of the approximate analysis is in this case that
more or less localized imperfections of the type considered here are equally
harmful as periodic imperfections with an amplitude reduced by a factor of
PIER
A different type of amplitude modulation of short-wave buckling modes
may occur in stiffened plate or shell structures due to the interaction of local
panel buckling with overall instability. The possibly disastrous conse-
quences of this interaction in some cases has been generally recognised
since Van der Neut’s pioneering paper.’’ The possibility of an approximate
analysis by means of amplitude modulation is here based on the property
that overall instability is characterised by a relatively large wavelength of
the deformation pattern.'!~ '*
Our discussion of amplitude modulation of short-wave buckling modes
in previous papers and also here is admittedly of a somewhat rough and
approximate character. It is essentially an application of the Rayleigh—Ritz
method, starting from assumed deflection patterns and neglecting ‘small’
terms in the subsequent evaluation of the energy. It may therefore be
worthwhile to note that a more mathematical approach, employing
multiple scaling in an asymptotic analysis, may serve as some justification
of our engineering analysis (cf. e.g. refs. 15, 16, 18).
Wee.
ii 9x eel DW Vy Vy,
ail ah Fe
a5 2W,)
R, 1 o
aye Sy tO Waele
x oy Wied
+40 eat +3(1 —v)(u,, +0, + W,.Wsy ,
Eh (2)
+ 2(1 — v)W,3y |dA
24(1 — v2) {Wy24 + Wyoy + 2VW,..Woyy»yy
This energy functional presupposes conservative external loads, but for the
shallow buckling modes under consideration it is immaterial whether the
loads are so-called dead loads or an external fluid pressure.’ According to
the general theory of elastic stability,*”’ initial imperfections Kw°(x, y),
where w° specifies the shape and x is a non-dimensional measure of the
magnitude, are taken into account to a first approximation by adding to
eqn (2) the integral
1 y
a(t a5 V)U, xy as ne om V Urey + Vy yy a ( an 7] Wsy = 0
: R, R 1 /
where at least one of the wave numbers p and q is large in shallow buckling
modes. The ratios A/C and B/C are obtained from the first pair of eqns (4):
h? v)
——,
4c?R3
(p* +q’)* + (p? +497)? + = (p?+9°)(pPopt+q?ae) (7)
ie
Eh?
|HC A068 + ded) +553c bu | 0)
The post-buckling paths are obtained by putting the partial derivatives
with respect to the non-dimensional amplitudes by and c,, equal to zero.
The result is
bo = —2(1 — A)/3e C= —e §(l—A)/3e (10)
and the deflection in the unstable post-buckling path is given by
8 16
ay =— (1-A) a, = +—(1-A) (15)
FASC } 9c
Let x = 0 denote the equator in a sphere and y/R the longitude along the
equator. Hutchinson’s mode,°*
may give an adequate description near the equator but it fails completely
when the poles are approached. This defect may be removed by the
introduction of an additional modulating factor (sin 6)‘ in eqn (18), and
also in the associated in-plane displacement components, where @ is the
polar angle and & is an integral number not less than 2 and small in
comparison with m. All singularities in middle surface strains and changes
of curvature near the poles which are present in the deflection pattern
eqn (18) disappear by the introduction of the factor (sin 6)*. Our
assumption that k is small in comparison with m implies that all strains and
changes of curvature, obtained from eqn (18) and the associated in-plane
components, are also multiplied approximately by a factor (sin 6)*. The
relative error in the evaluation of the energy is of order k?/m*. We consider
in particular the case that k is large in comparison with unity but small in
comparison with m, say k is of order m'/*. The buckling pattern is now
confined to an equatorial zone. The integrals of the quadratic and cubic
terms in the energy functional now have a factor
|(sin @)7*"*
dé and |(sin 6)3**1d@
(0) 10)
respectively. For large values ofk the ratio of these two integrals is x) hee In
other words, the cubic term in the energy is reduced by a factor J 2/3 in
comparison with the quadratic term. As a consequence the limit load factor
for imperfections described by eqn (16) multiplied by a factor (sin 6)* is
given by
Amplitude Modulation of Short-Wave Buckling Modes 41
E*(x, y).=exp[—4n7
(x? + y7)/R?] (19)
where pu’ is a fairly large number (in order that the modulated periodic
shape decays sufficiently rapidly) but at the same time small in comparison
with m?. The latter condition enables us to treat the modulating factor
again as a constant factor for the purpose of differentiation of the mode.
The behaviour of the shell with the more or less localised imperfections
described above is analysed again by the Rayleigh—Ritz method. The
deflections are assumed in the same shape as the imperfections, in the
circular cylindrical shell in the form
The quadratic and cubic terms in the energy functional in the modulated
modes eqns (20) or (21) are now obtained from the similar terms in the
periodic modes by a multiplication of their integrands by the square and
cube of the modulating factor E*(x, y) respectively. In view of the rapid
decay of the exponential factor the integration may now be carried out
from — © to + «© in both xand y directions. In both directions we arrive at
quadratic integrals of type
The cubic integrals are of the same type with p? replaced by (3/2)? and
their value for n=0 is thus \/2n/y,/3, a factor ./2/3 smaller.
The essential part of the energy expression eqn (9) and its counterpart for
the spherical shell is the expression between the brackets. The single
modification of the expression between the brackets in the case of
modulated modes is that the coefficient of the cubic term is decreased by a
factor 2/3 in comparison with the quadratic term. The critical load factor
remains unity, and the post-buckling paths in the absence of imperfections
are described for the cylindrical shell by (cf. eqn (10))
bb = —(—-A)j/e Cn = £41 —A)/e (23)
and for the spherical shell by (cf. eqn (15))
4 8
a= (1-4) = 45 (1-4) (24)
Imperfections of types eqn (12) or eqn (16), modulated by the factor
eqn (19), now result in a reduced limit load factor A* given by
(1 — A*)? = —4ckd*
for the cylindrical shell, and by
(l= 7"). ScK At
for the spherical shell. More or less localised imperfections of the type
considered here are equally harmful, in the case of an inward dent, as
periodic imperfections of an amplitude reduced by a factor 2/3.
It may be worthwhile to repeat that our present Rayleigh—Ritz type
analysis is quite approximate in character. One of the most doubtful
aspects is perhaps the treatment of the modulating factor E*(x, y) defined
by eqn (19) as a constant factor for the purpose of differentiation of
displacement components. The effect of derivatives of E*(x,y) has been
investigated in ref. 1 in the case of a cylindrical shell under axial
compression. For the value y?/m* =0.1, representing a fairly localised
dent, Gristchak has established that our simple first approximation has an
error of no more than about 11 %.
loads were exhibited most clearly in Van der Neut’s pioneering paper at the
1968 International Congress of Theoretical and Applied Mechanics.!7 In
his analysis the load-shortening relations for buckled flat plate strips with
imperfections were applied to a column built-up from such plate strips. A
mathematically entirely different, but physically nearly equivalent
approach consists of amplitude modulation of the basic short-wave local
mode.*'!!~!3 In this alternative approach the amplitude of the basic local
mode is taken to be a slowly varying function f(x) of the axial coordinate.
The errors introduced by this assumption and its various consequences are
of order /*/L?, where / is the wavelength of the short-wave local mode and L
is the wavelength of the overall mode.
Amplitude modulation of a single basic short-wave local mode yields
quite satisfactory results, as described in the references given above, and
exhibited in the oral presentation, whenever the deflections in the local mode
are largely confined to one side of the neutral axis of the overall mode. This
happens to be the case for a number of well-designed stiffened panels with
rather solid or hat-shaped stiffeners. In such cases the effect of
imperfections in the local mode is quite significant for very small
imperfections but it decays quite rapidly for larger imperfections. The
result is that the reduction in critical load due to such imperfections may be
bounded above by some 10%.
The sensitivity to local imperfections may even be reduced further, if the
basic local mode involves significant deflections on both sides of the neutral
axis for overall buckling. The reason is that in this case the overall mode
tends to increase the magnitude of the deflection in the local mode only on
one side of the neutral axis and to decrease it, approximately equally, on the
other side. An adequate analysis of this case, however, requires that at least
one additional local mode of the same short wavelength should be taken
into accountin order that a combination of the two local modes may
involve dominating deflections on only one side of the neutral axis. The
analysis becomes here more complicated, and we refer to a first attempt to
cover such situations in a previous paper.'*
REFERENCES
Mech., M. Hetenyi and W.G. Vincenti (eds), Berlin, Springer Verlag, 1969,
pp. 389-99.
. POTIER-FERRY, M., ‘Amplitude modulation, phase modulation and locali-
zation of buckling patterns’, Paper presented at IUTAM Symposium on
Collapse: The Buckling of Structures, 1982. To be published by Cambridge
University Press, 1983.
4
J. LOUGHLAN
Cranfield Institute of Technology, UK
and
A. R. UPADHYA
National Aeronautical Laboratory, Bangalore, India
ABSTRACT
NOTATION
1. INTRODUCTION
2. THEORETICAL APPROACH
The axial load P is applied on the axis of symmetry of the section and at a
distance e from its minor principal axis thus causing a moment M about the
web centre line given by
M =P(e+d+6) (2)
By considering a representative section of the column situated between
two local buckling nodes and treating this as a short strut of length 4 which
has its ends compressed according to the displacement system shown in
Fig. 3, we can obtain expressions for the moment M and axial load P in
terms of the displacement parameters u and « of the compression system at
the buckle nodes.
To obtain these, consider first the cross-section shown in Fig. 4 with
initial imperfections wy) and total deflections w. The local buckling
deflections are described by the expression
1X
w; = AY;(y;)
cos He (3)
BUCKLED _FORM_AT
COLUMN CENTRAL
SECTION
P(APPLIED LOAD)
WEB BUCKLES
PRINCIPAL STIFFNESS
DIRECTIONS BUCKLE NODES
Fig. 1. Central portion of column showing local buckles.
FLANGE
b=,
WEB
Fig. 2. Cross-sectional notation of plain channel.
52 J. Loughlan and A. R. Upadhya
+. U(1-%)
where A and Ay are magnitude coefficients for the local deflections and
initial imperfections respectively and Y,(y;) is the locally deflected shape of
the cross-section. It is assumed that the shape Y;();) is the same as the
buckling mode of the locally perfect section and this is obtained very
accurately from the multi-term Rayleigh—Ritz analysis outlined in ref. 10.
Although an infinite number oflocally imperfect shapes are possible for the
section it is convenient in analysis to assume that the imperfections are
similar in form to the buckling mode. This not only leads to simplicity in
analysis but has been shown by Rhodes et a/.'’ to give results for plates
which closely approximate actual plate behaviour.
LOCAL IMPERFECTION Wo
ToTAL DEFLECTION W
To obtain the stresses within the section which comply with the
deformations that take place during loading use is made of von Karman’s
compatibility equation written in the following form to take due account of
imperfections:
ot 4 ot
F, F, ;
(43>); l aay.
ex Ht (ds5); ax?
Ax dy2 + (414); oye
=
62w, l
\? 2, 62, L u =
07Wo.\2
1 ay
07Wo, 07 Wo,
1 1 5
je Ox at (ee Ox MOY, | e
where F; is a middle surface force function for the ith plate such that
0?F, 0° F, 0?F,
Nx), => Ny) See Nxy), = — 6
Substituting for w, and Wo, from eqns (3) and (4) into eqn (5) we can now
solve for the stresses and, by applying the principle of minimum strain
energy as outlined in ref. 10 for the locally perfect column, we may write an
expression relating end displacements at the buckle nodes and out of plane
deflections as
a Ao
C3
(ET). = (12)
where
: = b3 = ( bo G ~) ca
: (441)2 L3 X, \C@y1)2_ 33C C3 2C;
and
ee ee
eer eles
2D OG (aa types (dia
i Ao jeAe Ae
(Perucle =
ame hera ear Ae) 89
Locally Imperfect Plain Channel Columns 55
When A, =0 and A £0 eqns (12), (13) and (14) are identical to those
developed in ref. 10 for a locally perfect column which has undergone local
buckling. It is also of note that when A, #0 the values of ([E/]*)., (d*),
and (Pcruc), tend towards those obtained from the equations developed in
ref. 10 for large values of A. This is due, of course, to the fact that the
equilibrium behaviour of a locally imperfect column tends towards that of
its geometrically perfect counterpart well into the post-buckling range.
When A, = 0 and the effects of local buckling are not present in eqns (12),
(13) and (14) then ([E7]*), becomes (£7), the flexural stiffness of the locally
perfect unbuckled column, (d*), becomes d, the neutral axis position ofthe
locally perfect unbuckled column and, as expected, (Pcryc). =9 since
Pecrucis only associated with a column which has undergone local buckling.
It is evident from eqns (12) and (13) that due to the presence of even a small
value of local imperfection the flexural stiffness and neutral position of the
column cross-section are always different from those of its locally perfect
counterpart, even at zero load, 1.e. when A = Ao.
z= P 5” b?
2 A2 — 42
2¢
Resid iS
: ee CRSb ae 0) xX.
b2 Cab?
Ple+d+6,-—~~ |+ (A?
— A2)| C,b, -—*+-
s"= (a,,)2 X, (4, 1)2X, (16)
E b3 E b,
(4;,)23 reel
Substituting for é«, €é and 6” from eqns (10), (15) and (16) respectively
56 J. Loughlan and A. R. Upadhya
where
NO oe
3x.)
F 2CP , C2X(ay1)2
=a
ees 4©:
_—— —Q _ |NN)
aIs
ae
oa
ne ; aor
mds 2G) G oy i
1 ae —C 2(e+d
é anal 3 2|+* oe (4,1).
na {oe 6
; b, (441)2X,
Since the parameters ({EJ]*),, (d*), and (Peruc). vary with the local
deflection magnitude coefficient A then eqn (11) cannot be used directly to
obtain the load—deflection equilibrium path of the deforming column.
Using eqn (17), however, in conjunction with eqn (11) the equilibrium path
can be obtained quite easily by solving these numerically in a step-by-step
procedure using iterative techniques.
Knowing that the initial state of the system 1s P = 0, A = A, and 6, = 6,
then, for a value of 6, slightly in excess of 0,,, a realistic ‘guess’ at the
equilibrium value of P corresponding to this can be made and hence A can
be obtained numerically from eqn (17). Knowing A the values of ({E/]*)
(d*), and (Pcryc), can now be obtained from eqns (12), (13) and (14) and
by substituting for these in eqn (11) a numerical evaluation of Pcan now be
made for the given 6,. If the value of P obtained from eqn (11) is the same as
that used in eqn (17) then obviously this is the required solution for the
given 6,. If itis not, however, as is more likely, then the value of P obtained
from eqn (11) may be substituted into eqn (17) for the same 6, and the
process repeated until the difference between two consecutive values of P is
zero. The solution will then be complete for the given 6, and to establish the
complete load—deflection equilibrium path of the column the procedure
outlined is repeated for a number of small increments in 6,.
So far in the analysis it has been considered that the locally buckled shape
of the cross-section remains constant with change in load and that it is of
Locally Imperfect Plain Channel Columns 57
the same form as the buckling mode of the locally perfect section. Since the
parameters C, to C, are functions of the locally deflected shape and are
therefore constant for a given shape, the implication of this assumption is
that changes in the values of ([E/]*),, (d*), and (Peruc), occur only with
change in the local deflection magnitude A for any given imperfection
amplitude A,. During loading, however, changes in the locally deflected
shape do in fact occur, caused principally by column bending, and the
coefficients C, to C, are variable which means that changes in the values of
({EJ]*)., (d*), and (Peryuc). depend, significantly, on the current locally
deflected shape. The degree of this significance can be considerable for
certain cross-sectional geometries and an analysis which precludes the
possibility of a changing locally deflected shape can give results which are
seriously misleading.
In this paper local form change is accounted for in an approximate
manner and a study ofits effects on column equilibrium behaviour is made.
The indirect approach considered in analysis first makes use of the readily
available facility of being able to obtain a range of buckled shapes
corresponding to different compression eccentricities from the multi-term
Rayleigh—Ritz analysis described in ref. 10. Using these shapes it is now
possible to obtain a series of solutions describing the load—deflection
equilibrium behaviour of a column, each solution being based on a
different buckled shape and hence on different values of ({E/]*),, (d*), and
(Peruc)e: Lhe lowest envelope of the equilibrium paths obtained in this way
is then considered to be the true solution since for a given column deflection
the solution which gives the lowest axial load is that which prescribes the
least value of the strain energy of the system.
The approach outlined gives results for locally perfect columns which
sufficiently describe the effects of local form change on column behaviour.
The method can also be used, quite easily, to account for these effects for
columns with local imperfections. In the case of the imperfection analysis
presented in this paper, however, the implication of using this approach is
that as well as the locally deflected form changing with load the form of the
local imperfections also changes. Allowing for this fictitious change in the
local imperfection shape with load, however, gives results which are
considered by the authors to be realistically representative of the practical
situation for the case of columns which are sensitive to the combined effects
of local imperfections and local form change.
associated with a variation in the magnitude of the local buckles along the
column. In practice the local buckle magnitudes are more pronounced at
the centre of the column and diminish along the length of the member. This,
of course, is intuitively what one would expect and in fact was evident some
thirty years ago in the work by Bijlaard and Fisher'*:'* who conducted a
series of tests on centroidally loaded I section and square tube members for
which local buckling occurred prior to overall column buckling. In this
work there is clear photographic evidence relating to an I section column
which highlights the proportioning of the local buckle magnitudes along
the length of the column due to overall bending. Due to centroidal loading
and the symmetry of the cross-section, column bending, in this case, is
initiated by overall bifurcation of the locally buckled member and the
equilibrium behaviour after bifurcation is shown in ref. 14 to be fairly
unstable for I sections. In this instance, therefore, the variation of local
deflections along the column does not affect the ultimate carrying
capability of the member but has a direct influence on post-failure
equilibrium behaviour only.
In the case of pin-ended, centroidally loaded, singly symmetric sections
such as lipped and plain channels for which local buckling occurs before
Euler buckling, local deflection variation along the column occurs from the
onset of local buckling and in this case will influence the behaviour of the
column prior to failure. In this paper the approach used considers
essentially the behaviour of a single local buckle along the column and from
this consideration the governing differential equation for the lateral
equilibrium of the locally buckled column is developed (see eqn (20) in
ref. 10). This method of approach is considered to lead to a somewhat
elegant solution to the interactive buckling problem and has been shown to
give results which are in good agreement with experiment. Typical
comparisons showing this agreement are to be found in refs. 6 and 8 for
plain and lipped channel columns respectively.
The use of the differential equation developed in the analysis
automatically introduces a form of local deflection variation along the
column and this can be examined by suitable manipulation of the relevant
equations in the analysis. On examination it has been found that the local
deflection magnitudes do in fact diminish along the column in the presence
of column bending as has been observed in practice.
For geometrically perfect eccentrically loaded columns, depending on
eccentricity magnitude, the analysis can lead to a discontinuity of the
column deflections between the pre- and post-local buckling solutions at
the local buckling load. The reason for this is obtained from an
Locally Imperfect Plain Channel Columns 59
The numerical work involved in the analysis was carried out using the
VAX-11 computer system at Cranfield Institute of Technology and the
main scientific programming language FORTRAN 4 was used to write
the computer program. Although an extensive amount of data is available
from the computer the authors have elected to present results for isotropic
columns and have chosen these to highlight, in their view, the most
significant effects of local form change and local imperfections.
Figure 5 shows load-deflection equilibrium paths and reduced Euler
load curves P*/P,. for a centroidally loaded perfect column with a flange-
to-web ratio of 0.4. The broken curves are based on an unchanging locally
deflected shape with load and the full lines take local form change into
consideration in the approximate or indirect manner described earlier. The
length of the column is such that its local and Euler buckling loads are
simultaneous and due to this the nature of the equilibrium at buckling is
unstable. From the figure it can be seen that the immediate change of
flexural stiffness for this section at local buckling is given by P£/P, ~ 0.71.
If post-local buckling behaviour is now based on the deflected shape
obtained at buckling then the load—deflection equilibrium curve for the
column would tend asymptotically towards the P/P, value of 0.71 as
indicated. The actual behaviour, however, is seen to be quite different from
60 J. Loughlan and A. R. Upadhya
b¢
ie
Ar
DOCHcud,
a7 WC e.pian Didare anOMe nieee a
t =0.0, ‘la =0-4, +-= 100, bw 72h 49
Fig. 5. Effect of local form change on equilibrium behaviour and reduced Euler
load (é = 0.0).
0.8
Gee
doc _-99, Woc
“96-00,
_ Pf
bie-0.4, ¥-100
bw _
+e © LcAD PoinT
10 i; ale bs X NEUTRAL AXIS
08
aes
Dit ell a
oa 13-57
ie
Fe
ad @) 2 4 6 8 10 oe
00¢
t
0.0, Woe
t
=0.0, Bw
bf =0.4, BW
t
=100
Fig. 7. Effect of local form change on equilibrium behaviour (é= —0.1).
10
Bt os , &=0.0
Bw
Woc Pw = 833 , = 0.0
0-57 ‘
N WwW Wo y
Vie |a
Dy
WceMoc-web cent deflection
5 ( We “oc)/t
0 1 2 3 4
Fig. 8. Reduction in flexural stiffness for locally imperfect section.
°9 1 2 3 4 Bam cr Oct
Fig. 9. Shift of neutral axis for locally imperfect section.
64 J. Loughlan and A, R. Upadhya
When local imperfections are present the local deflections begin to grow
from the onset of loading and these are accompanied by a gradual
reduction in the flexural stiffness and a gradual change in the neutral axis
position as shown in Figs. 8 and 9. When the imperfection magnitudes are
small the initial values of flexural stiffness (E/)* and neutral axis position
d*, at zero load, are close to those of the locally unbuckled perfect column,
but the rate of change ofthese values with the growth of local deflections is
high. For larger values of imperfection magnitude, (E/)* and d* change at
a smaller rate with increasing local deflections but their initial values at zero
load are progressively smaller, with increasing imperfection magnitude,
than those of the locally unbuckled perfect column. As the local deflections
increase both (E/)* and d* approach the corresponding values for the
locally buckled perfect column for all imperfection magnitudes. It is seen
from Figs.8 and 9 that for the unloaded column an imperfection
magnitude of w,,/t =0.5 results in a 20 % loss of section flexural stiffness
and a 26 % shift in neutral axis position towards the web centre line. For an
imperfection magnitude of wo,/t = 1.0 the loss in flexural stiffness is 43 °%
and the shift in neutral axis is 57° towards the web centre line.
Figure 10 shows the effect of local imperfections and local form change
0 4 8 12 16 20 éc
1.0
0.9
bt
ey =0-4,
bw
t = 100,
L
bw = 2t49,
a
e =(:0)
Soc_
t =00
5. CONCLUSIONS
The ultimate loads obtained from the solution which takes this into
account have also been shown to be lower.
REFERENCES
. VAN DER NEUuT, A., ‘Mode interaction with stiffened panels’, Proc. [UTAM
Symp. on Buckling of Structures, Harvard University, 17-21 June 1974,
pp. 117-32.
. Korter, W. T. and PIGNATARO, M., ‘An alternative approach to the interaction
between local and overall buckling in stiffened panels’, Proc. 1UTAM Symp. on
Buckling of Structures, Harvard University, 17-21 June 1974, pp. 133-48.
. WALKER, A. C., Interactive buckling of structural components, Science
Progress, 62 (1975) 579.
. Fox, W. C., WALKER, A. C. and RuHopes, J., Buckling of locally imperfect
stiffeners in plates, Proc. of the ASCE, J. of the Eng. Mechs. Div., 103,
No. EMS, October 1977.
. GRAVES SMITH, T. R., The effect of initial imperfections on the strength ofthin-
walled box columns, J. Mech. Sci., 13 (1971) 911-25.
. RHopes, J. and Harvey, J. M., ‘Interaction behaviour of plain channel
columns under concentric or eccentric loading’, Second Int. Colloquium on the
Stability of Steel Structures, Liege, 13-15 April 1977, pp. 439-44.
. LOUGHLAN, J. and Ruopes, J., ‘Interaction buckling of lipped channel
columns’, Stability problems in engineering structures and components, T.H.
Richards and P. Stanley (eds), London, Applied Science Publishers, 1979,
pp. 179-98.
. LOUGHLAN, J. and Ruopes, J., ‘The interactive buckling of lipped channel
columns under concentric or eccentric loading’, International Conference on
Thin-walled Structures, University of Strathclyde, Glasgow, Scotland, 3-6
April 1979.
. Hancock, G. J., Interaction buckling in I-section columns, Proc. of the
ASCE, J. of the Structural Div., 107, No. ST1, January 1981.
. Upapuya, A. R. and LOUGHLAN, J., ‘The effect of mode interaction in
orthotropic fibre reinforced composite plain channel section columns’, In:
Composite Structures, 1. H. Marshall (ed.), London, Applied Science Publishers,
1981.
. Ruopes, J., HARvey, J. M. and Fox, W. C. The load-carrying capacity of
initially imperfect eccentrically loaded plates, /nt. J. Mech. Sci., 17 (1975)
161-75.
. Upapnya, A. R., ‘A study of the buckling behaviour of composite reinforced
metal panels and composite plain channel section columns in compression’,
Ph.D. Thesis, Cranfield Institute of Technology, Cranfield, Bedford, England,
September 1980.
. BIJLAARD, P. P. and FisHeEr, G. P., ‘Interaction of column and local buckling in
compression members’, NACA Technical Note 2640, 1952.
. BULAARD, P. P. and FIsHER, G. P., ‘Column strength of H-sections and square
tubes in post-buckling range of component plates’, NACA Technical Note
2994, 1953.
An Examination of Shakeout in Stiffened Plates
ABSTRACT
NOTATION
Cross-sectional area.
Smallest and largest eigenvalues.
ok
a o>
oO Plate length and width.
Tangential axial rigidity.
S Gy z Damping coefficients for u and w displacements.
Tangential flexural rigidity.
oh
ao
eke™~7) Flat-bar stiffener dimensions.
Modulus of elasticity.
>= Damping factors for u and w displacements.
Column length.
Cross-sectional moment.
TE
Ah Axial load.
69
> Oo P. A. Frieze and E. Drymakis
Subscripts
a End of a time increment.
b Beginning of a time increment.
i Node reference number (Fig. 2).
Z At distance z.
Superscripts
Value non-dimensionalised with respect to value at yield.
/
1. INTRODUCTION
Both the pattern and magnitude of these stresses vary enormously. Their
effect on the response of structures also demonstrates a wide variation
which is, however, probably more dependent on the structural action in the
member concerned than on the source of the stresses themselves.
The effect of hot-rolling is probably most noticeable in structural
sections of thicker proportions. This process tends to produce stresses
which show little variation along the length of the member but significant
changes across the section. The latter arise from differential cooling of the
section following rolling which generally results in tensile stresses being
developed in the thickest portions, such as the flange—web junctions, and
equilibrating compressive stresses elsewhere. These compressive stresses
can reach 50% of yield and since they occur at the tips of flanges have
proven to have a severe weakening influence on such sections subjected to
compression loading.
Cold-rolling of plates tends to produce stresses again uniform along the
length and varying through the thickness. These arise from plastic working
of the material. They are generally not large in magnitude but usually result
in a reduction in stiffness before yield stress is reached under both tensile
and compressive loading. In coupon tests this produces a ‘rounded’ stress—
strain curve. From a practical viewpoint, if the plating is subjected to
further processing, such as cold-bending or welding, the stresses introduced
by these are usually of greater magnitude so that the cold-rolling ones can
be considered as secondary.
Cold-bending particularly of sections can result in significant residual
stresses. These exhibit a ‘zig-zag’ pattern through the height of the section
and in particular circumstances can seriously reduce the strength of
members in flexure.' The magnitude of these stresses and their
consequences can be reduced by the use of double-bending where the
member is over-bent before bending back to the required curvature.
Flame-cutting of plates tends to produce a pattern of stresses similar to
the longitudinal ones which arise from welding, viz. tensile along the line of
the heat input with equilibrating compressive stresses across the remainder
of the plate. As flame-cutting is generally precursory to welding, its residual
stresses tend to be secondary in their effect on structural response.
Welding produces complex patterns of stresses which are manifested in
12 P. A. Frieze and E. Drymakis
two ways. First, there are those directly associated with the weld and the
immediately surrounding parent material, the ‘heat affected zone’. These
are usually tensile in nature and of a magnitude at least equal to the yield
stress of the material (property changes sometimes occur during welding to
produce a material of higher yield stress). Second, there are those which
arise as a consequence of maintaining equilibrium with the weld-associated
tensile stresses. These are obviously compressive in nature, unless the plate
is unrestrained when shrinkage will occur instead, and may be of a
magnitude up to 40% of yield if the welding is predominantly uni-
directional. Considerably higher values may result if orthogonal weld
patterns are adopted. As in the case of hot-rolling of sections, it is the
compressive stresses that generally are of concern from a structural
strength point of view since, in conjunction with compressive loading, they
precipitate either premature buckling or yielding both of which tend to lead
to some form of general or overall instability.
In attempting to deal with residual stresses, designers have normally used
one of the following:
2, -MODEL
procedure. The one used in the present work is finite difference based and
assumes the following:
(i) Plane sections remain plane and normal to the neutral axis;
(ii) Torsional behaviour of the stiffener can be ignored;
(iii) The stress-strain curve for the stiffener is identical to that of the
basic material, here steel;
(iv) The stiffener is initially stress-free;
(v) The axial response of the plating can be determined via a separate
stress-strain curve;
(vi) The plating is initially stress-free although the specification for the
plate from which its stress-strain curve is derived may contain
initial stresses;
(vii) Deflections normal to the longitudinal axis are restricted to the
plane of the stiffener;
(viii) Shear and normal stresses are negligible.
where M and P are the moment and axial load acting on the cross-section
and w, and w are the initial and ‘under load’ deflections normal to the
longitudinal x-axis, and compatibility by
d*w
as (3)
and
du 1/dw\? dwdw,
: = +3 ($) + ax dx 4
where ¢ and ¢ represent the curvature and axial strain in the section, and wis
the axial displacement.
An Examination of Shakeout in Stiffened Plates 75
Pai o,dA
A
=| E,dA +6 |E,zdA
A A
=Ce+G¢ (7)
m= | o,zdA
A
=) B24 +6 |E,z’dA
A A
=Ge+D¢ 7 (8)
where C, G and D are referred to as tangential rigidities. Since E, is stress
dependent, C, Gand D have to be evaluated frequently as loading proceeds.
This process is facilitated by recasting eqns (7) and (8), and (3) and (4) in
incremental form.
The numerical procedure used is referred to as dynamic relaxation or
DR. It is an iterative procedure which uses eqns (1) and (2) to generate out-
of-balance forces until equilibrium of the solution is achieved. These out-
of-balance forces are applied to equations of motion for the displacements
corresponding to the equilibrium conditions, w and uw respectively, as
follows:
du ;
Py ae + C,u = [out-of-balance, eqn (2)] (10)
76 P. A. Frieze and E. Drymakis
| PLATE WIDTH b |
ts
Fig. 1. Beam-column model.
where (') represents differentiation with respect to time ¢, and p and C are
densities and damping coefficients.
Introducing non-dimensional damping factors of the form
Che
K as (11)
Pe
where At is a time increment, and writing eqns (9) and (10) in finite
difference form, the velocity at the end of an increment (represented by
subscript a) can be found in terms of the velocity at the beginning of an
increment (subscript b) by rearranging:
) (alfa),&
4D |P| G, ( dw dwo d?w d?wo
Pwo = (xy? * (Ax)? | 4Ax \“ fax
—
dx l
dx? dx L
| = Nl Ax =
Damping factors are selected to lie just below critical so that oscillations
can occur about the static values. Apart from providing more rapid
convergence than damping just in excess of critical, it provides a visual
means by which convergence can be evaluated in printed output. The sums
of the square of the w and u velocities are also used to provide a second
means of checking convergence.
a OE STRESS LEVEL
Sf,
PLATE WIDTH b
Fig. 3. Idealised welding residual stress distribution.
stiffeners since this was the only form of stiffener incorporated within the
procedure.
10
PRESENT
METHOD
0:8
0°6
0-4
0:2
1-0
PRESENT
METHOD
0:8
0-6
0-4
0-2
4. PARAMETRIC STUDY
The results showed that strength did not vary significantly with aspect
ratio and that an a/b of 0.8 would be most appropriate.
Using this aspect ratio, plates of slenderness 1.04 and 2.07 (b/t = 60 for
mild steel) were analysed in the presence of 25 °% compressive residual stress
and an initial bow of 0.116ft (=0.145ft x 0.8) when loaded first by a
limited amount of tension followed by compression up to failure and into
the post-buckling range. The extent of tensile preload used were strains
equal to 0, 10, 25 and 40% of tensile yield. The results are presented in
Figs. 6 and 7 in the form of average stress-strain curves where both stress
and strain have been non-dimensionalised with respect to their yield values.
82
50
NOISN3IL OVO1-3Yd
%) 49 NIVYLSOT3IA
(
S:0-
P. A. Frieze and E. Drymakis
YA)
advIDAY UIBIIS—Ssd1]
SIAIND
s JOJ
9
passordu09
sajejd YIM snoueA saazZap
Jo uorsua) ‘peoj-aid
9 = ‘p01 49/2 ‘80=
*2 =
“Sly
75Gc
90= O FONT
ZZ NOISN3L QVO1-3ud
S°0
%) JO GIZA (NIVYLS
-0-
An Examination of Shakeout in Stiffened Plates
From the figures it can be seen that the initial loading into tension has a
slope less than unity. This is partly due to the presence of the initial bow
which tends to reduce the stiffness of plates under any in-plane loading and
will be most pronounced for the more slender plates. It is predominantly
due to the edges ofthe plate being initially at tensile yield (see Fig. 3) so that
any tensile loading will render these strips ineffective.
Upon loading in compression, the initial slope appears to be almost
independent of the initial level of tension. This is not unexpected as this is
dictated by the level of initial bow, and for the degree of tension loading
applied little change in this was noted. In-plane yield as manifested by the
rapid change in slope at a compressive strain ratio of about 0.7 is seen to be
affected by the extent of pre-load. Beyond that point, the differences in the
responses are maintained until they converge upon a single post-buckling
curve. From previous studies > this is known to be the post-buckling curve
for the same plate when initially free of any stresses.
This tendency to converge to the initially stress-free result implies that
the strength of stocky plates will be independent of the degree of tension
applied, and that maximum differences will probably occur in the
slenderness range where, for the level of initial bow considered, residual
stresses have their most pronounced influence, i.e. 6 between 1.7 and 1.8.
Comparing Figs.6 and 7 shows that greater differences are obtained
between the results for different pre-loads when 6 =2.07 than when
B = 1.04.
From acompressive strength viewpoint, the influence of tension pre-load
is not marked. For the f =2.07 case, the maximum enhancement in
strength, due to the 40 % initial tensile loading, is some 7%.
1:0
A=0:57
TENSION PRE- LOAD
(% OF YIELD STRAIN )
iS 0:4
TENSION PRE-LOAD
(% OF YIELD STRAIN)
= 0:4
Stith
Fig. 8. Average stress-strain curves for stiffened plates subjected to various
degrees of tension pre-load, $ = 1.04, «=0.4, Ay =//1000.
86 P. A. Frieze and E. Drymakis
d/t,, was set to 10. The ratio of stiffener area to plate area was generally
taken as 0.4.
Typical results in the form of non-dimensional stress—strain curves are
presented in Figs. 8 and 9. For Fig. 8, a plate of slenderness # = 1.04 and an
initial bow of magnitude A, =//1000 was used: the corresponding values
for Fig.9 were 2.07 and //500. For clarity, generally only the results
pertaining to tension pre-loads of 0 and 40 % are illustrated.
Neal 0)
Oe TENSION PRE-LOAD
(% OF YIELD STRAIN)
ee 0-4
31V1d
40-1 4989LL-0 %SZ
L0:2 0 499Lb %SZ
Oe
20
o”7
SENEFEDIES
©e)
/
000b/2
005
7-0 9-0
S-0
NOISN3L OVO1-3Yd
%) 40 CT3IA (NIVYLS
An Examination of Shakeout in Stiffened Plates
5. DESIGN RECOMMENDATIONS
Having shown that residual stress shakeout did not have a significant effect
on compressive strength, the analysis was then used to examine the strength
of plates stiffened by flat-bar stiffeners when subjected to compression only.
For stiffened plates subjected to a longitudinal stress gradient, the strength
of the plates at the end of the panel will sometimes dictate strength. Plate
strength is considered in section 5.1. For panels subjected to small or zero
stress gradients, beam-column strength will be paramount: this is discussed
in section 5.2.
i a/b = 08
g! 55 = 0-116 Bt
On (%)
08 0
----- 10
en 5
0-6 Saree
a aU
0-4
0-2
plate strength with slenderness ratio 6. These can be used for determining
the strength in longitudinal compression of a practical range of plates.
EULER PLATE
a! B 4 Op
nis @ 1:04 o16pt 25%
@ 2:07 ON6Bt 25%
STIFFENER
0:6
0-4
0-2
6. CONCLUSIONS
REFERENCES
7 ae ;
ears —
eae i 3 7
ete C pletene ee yhted haat
aaah ce \ ees 7a om od
tarrite (— “eet Ovencs>
AE
higtc Reon sans RET Aniwr nas tan 8 a
a5 da *) 0 © pot nes wi
HW toje" rA.Yenliey , 2 ine i) Tuite ie
4, 1
oo eer See
AeA treete trian ve
The Cold Formed Industry Scene
R. D. HOMER
Metal Sections Ltd, Oldbury, UK
ABSTRACT
In common with most manufacturing sectors, the cold formed industry scene
is dominated by the recent dramatic contraction of its markets. After
consolidation and retrenchment, the industry is currently seeking ways to
offset reduced and replaced lost markets. A study of past product development
can often reveal helpful leads or pointers and so the paper covers three
decades of such work. Whilst the developments are related to Metal Sections
Ltd, they are intended to represent the industry in general. Factory buildings
with three designs of main frames are described; an aircraft hangar door, a
product from a special variation of the cold formed process which found
unexpected major use in civil engineering and a composite building panel
with a particular example of use are also featured. In conclusion, the paper
deduces two factors valid for consideration. The reader is recommended not
to be put off by size and to include fields previously thought inappropriate.
1. INTRODUCTION
it was realised that this was likely to be a long-term trend and currently the
industry is urgently seeking ways to increase its penetration into existing,
and to find replacements for those lost markets.
A study of past product development can often reveal helpful leads and
pointers for formulating plans of action and so the paper is concerned with
this aspect. Since it is only possible to describe work in which I have been
involved, the paper will review some three decades of work at Metal
Sections Ltd. It is intended, however, that the developments described
should be regarded as indicative of the industry as a whole.
Appropriately, the decade of the 1950s—the period in which my story
begins—found Professor Harvey in industry and he took part in some of
the activities described.
The period immediately after the Second World War was an exciting
time, as finally the fetters of a war economy were being loosened. There
was a grave shortage of many things, even essentials, and a large unfilled
demand for a wide variety of industrial products. The cold formed industry
of those times was enterprising and entrepreneurial, even as it is now, and
willing to push out the frontiers of its markets as far as it could. One
area where a major effort was made was the building industry and it is
mainly, although not exclusively, in this field that I wish to describe a few
developments.
ans AN
<<
“Spence
10’0” (3.05 m) were supported at eaves and apex by Warren girders which
spanned the 40’ 4” between portals (a much greater unobstructed area was
specified for the Header building).
Specified design loads were as follows:
(1) The tender value accepted by the main contractor for the main
warehouse and the Header building’s structural framing was
£193,000.
(2) Atthe time the building was the largest in Europe by a considerable
margin.
(3) The time from the go-ahead to completing the weather proofing
was 14 months which was an all-time record for a building of this
size and which earned the main contractor considerable acclaim.
‘OSNOYIIVM UIPU JO SUOTIIIG “fp “BI
Is™GGy
6S) SOE]GO-Z
ys zLoLE = AAT
+e"“Sr ys"
"SN DEC / = xx>
«
=J
y"GrDE ‘sm ZZ) = “my
w
srcle} = A>
Pd lsoe = XXT
"SMIOS
yINGIG-F = KAT
SIG-O
SN'OS IPE T = YI “WN OSIGSE/ = yay
| =|
ie 2¢
x y id
A A
a dh 6 f
€
ize
ay
‘SUIP[INg JopeoH{ “(pluor—p “B14
sw ges} = XS 'SNIIQQ) = ALY
POOL ai AT aan: OSepy = AAT
SMOSODE|
7 y
of
ra
&
y
BS , T
|
“
=6 22
w
or
a i ==
=)
a] hee J)
t A
aS
Ins v7 =|
“ie
4
100 R. D. Homer
iia
561,
WE.)=
ae
(4) Cold formed sections were used entirely in the construction and the
greatest individual thickness of any section was 7,” (4.76mm),
although the gusset plates joining the member were, of course,
thicker.
(5) In spite of there being some corrosion problems (caused by the
proximity to the ICI chemical plants at Widnes and Runcorn),
Metal Sections made a ?% profit on the project.
Fig. 7. A building at Oldbury, West Midlands, using cold formed section two-
pinned arches.
LOADING SHACKLE
portals made from two cold formed channels welded back to back. Figure 9
shows a three 30’ (9.1 m) span bay x 80’ (24.4m) building. Designs were
made for 40’0” (12.2m) span frame with the same construction though
none were actually made.
Although the foregoing remarks principally concern the main frame,
complete buildings, clad, insulated and erected, were offered and cold
formed section solutions were used for gable framing, wind bracing and
other items of secondary steelwork. In addition, some conventional
triangular trusses were made using similar construction. Figures 10 and 11
show two examples.
Another project undertaken during the same decade was the design of
aircraft hangar doors (Fig. 12). Made entirely from Corten steel cold
formed sections each door measured 32’0” (9.75m) in height, 23’ 6”
(7.16m) long and 12” (0.30m) thick and the specification called for a
tolerance level of +4” (6mm) across diagonals. The doors were erected in a
hangar at London Airport.
104 R. D. Homer
During the 1960s, the most noticeable achievement was the development of
the Z purlin which is, of course, still a major product within the industry.
Early in this decade an interesting special development of the cold formed
process lead to the development of a machine for producing lock seamed
spiral tube. Figure 13 illustrates the product which was produced in
diameters from 3” (76mm) to 5’ (1.5m) and in a range of materials and
thicknesses. Originally, it was thought that the tube would find its main
outlet in high velocity air conducting but it eventually found its largest
volume use as voids and piling casings in civil engineering concrete
constructions (Fig. 14).
Towards the end of the decade when, yet again, there seemed a market
for industrial buildings, a composite polyurethane panel, together with a
comprehensive building system for its use, was developed. This found a use
in many types of buildings, including hospitals, offices, as curtain walling
for concrete buildings and was considered for domestic houses, particularly
for hot climates.
106 R. D. Homer
5. CONCLUDING REMARKS
oh S583 a
Fig. 15. Composite foam filled panel being used in house construction.
The Cold Formed Industry Scene 109
Ba a:
Fig. 17. View of finished house using foam filled panel walls.
several were also working along the same lines—the Brockhouse ‘Clasp’
system, for example, being a particularly successful project.
In conclusion, are there any pointers from the pastto helpus formulate
plans for the 1980 decade? Two pieces of advice seem relevant: one is not to
be put off by size of the project and, second, not to be afraid to look for cold
formed solutions in areas which it was previously thought impossible or
inappropriate.
Acknowledgement
The Author wishes to thank the Directors and Management of Metal
Sections Ltd, for permission to use material and data from their archives.
UMR Research on Cold Formed Steel Structures
WEI-WEN Yu
University of Missouri—Rolla, Rolla, USA
ABSTRACT
During the past 14 years, several research projects related to cold formed
steel structures have been conducted at the University of Missouri—Rolla
(UMR) under the sponsorship of the American Iron and Steel Institute
(AISI). The purpose of these investigations has been to study the structural
strength of cold formed steel members and bolted connections. Some of the
research findings have been used in the AISI Specification for the Design of
Cold Formed Steel Structural Members.
1. INTRODUCTION
Since 1969, eight different research projects have been conducted at the
University of Missouri—Rolla (UMR) under the sponsorship of the Ameri-
can Iron and Steel Institute (AISI) for the purpose of studying the
structural strength of cold formed steel members and connections. The
subjects of these investigations were the compressive strength of perforated
elements, performance of thick sheet steels and steel plates, web strength of
flexural members, effect of bending on web crippling of beams, bolted
connections, design of stainless steel structural members, development of
tentative recommendations on load and resistance factor design of cold
formed steel, and automotive structural components composed of high
strength sheet steels. Some of the research findings have been used to
111
112 Wei-Wen Yu
improve the AISI Specification for the Design of Cold Formed Steel
Structural Members.' The objectives and descriptions of each of these
projects are given below along with a discussion of the research findings and
their application to the design of structural members cold formed to shape
from steel sheet, strip or plate. Details of the research are contained in the
publications referred to in the text.
2. RESEARCH PROJECTS
reduction factors can be used to account for the effects of holes on local
buckling of unstiffened compression elements, shear buckling and web
crippling strength of webs when perforated elements are used.
Since 1973, research concerning the webs for cold formed steel flexural
members has been conducted at the University of Missouri—Rolla. This
study dealt with the structural behavior of cold formed steel beam webs
subjected to bending stress, shear stress, combined bending and shear, web
crippling, and the effect of bending on web crippling load. In addition, this
project has also included an investigation of beam webs reinforced by either
transverse or longitudinal stiffeners. The details of the investigations and
the research findings are presented in refs. 26 through 32. This information
was used extensively in the revision of Sections 2.3.4, 3.4 and 3.5 of the 1980
Edition of the AISI Specification.
In the study of bending strengths of beam webs, the postbuckling
strength of web elements was determined to be a function of the web
slenderness ratio, the bending stress ratio of the web, the flat width-to-
thickness ratio of the compression flange, and the yield point of steel. A
formula for calculating the postbuckling strength factor was developed
from this study. As an alternative method, two formulas were derived from
the test data for reducing the bending capacity of beams with large web
slenderness ratios. In addition, for the range of cross-sectional parameters
covered in this study, ‘effective web depth’ equations were developed for
determining the ultimate bending capacity of cold formed steel
beams.*°:*7:3! These formulas differ from those developed by Bergfelt,**
Hoglund,** and Thomasson.?°
For web crippling and a combination of web crippling and bending,
additional tests of channels and I sections were conducted for end one-
flange loading, end two-flange loading, interior one-flange loading and
interior two-flange loading to supplement the Cornell tests conducted in
the 1940s and 1950s.°°3’ In addition, the publications of Yamaki,*®
White and Cottingham,*? Rockey, El-gaaly and Bagchi,*°~*° Baehre,*’
Ratliff*® and others and the current design criteria used in various
countries have been reviewed in the initial phase of the investigation. Based
on the available test data, modified design formulas were developed for use
in Section 3.5 of the 1980 Edition of the AISI Specification.
2.4. Web Crippling and Combined Web Crippling and Bending of Steel
Decks
The modified design formulas discussed above are based on the test data of
channels, I beams and hat sections having vertical webs with relatively
small R/tand N/t ratios (R being the inside bend radius, ¢ the thickness and
N the length of bearing plate). It has been realized that steel decks often
have inclined webs with large R/t and/or large N/t ratios. In addition, for
UMR Research on Cold Formed Steel Structures ibis
the steel decks used for composite slabs, various types of embossments and
indentations are usually formed in the webs. These deformations may affect
the web crippling strength of steel decks. For this reason, additional studies
were conducted in 1979 at the University of Missouri—Rolla for the purpose
of studying the web crippling strength of steel decks under the
cosponsorship of Steel Deck Institute, American Iron and Steel Institute,
and H.H. Robertson Company. Based on the test data on one-flange
loading, it was found that the modified formulas for web crippling are
adequate for steel decks having flat webs, provided that the decks are
properly connected to support beams. However, the results of continuous
beam tests indicated that the interaction formula for combined bending
and web crippling developed from simple span beam tests of channels and
hat sections is conservative for multi-web steel decks used as continuous
beams if the localized deformation over the interior support is permissible.
This can be attributed to the different postbuckling behaviour of the steel
decks over the interior support and the moment redistribution that occurs
in the continuous beams. Consequently, additional studies have been
recommended for future research in this area.
prepare a design manual for stainless steel, cold formed structural members
under the guidance of the AISI Advisory Group on the Specification for the
Design of Cold Formed Stainless Steel Structural Members (D. S.
Wolford, Chairman). This document was prepared on the basis of the 1968
Edition of the AISI stainless steel specification on annealed grades°’ and
on the research work conducted by Johnson, Wang, Errera, Tang and
Popowich at Cornell University under the direction of Dr George
Winter.°’~°3
The AISI Stainless Steel Design Manual was published in 1974.°* It
contains the following four parts: specification for the design of annealed
and cold-rolled grades of sheet and strip stainless steels, commentary,
illustrative examples and design tables and charts.
economical use of these high strength steels in car bodies are needed by
structural designers and engineers.
In February 1981, the ‘Guide for Preliminary Design of Sheet Steel
Automotive Structural Components’ was issued by the American Iron and
Steel Institute for assisting automotive structural designers to achieve
weight reductions through the efficient utilization of carbon and high
strength steels.°® These design recommendations were based primarily on
the 1968 Edition of the AISI ‘Specification for the Design of Cold Formed
Steel Structural Members’ except for the following major differences with
respect to the AISI Specification written for the design of buildings:°?
In view of the fact that many types of high strength steels with yield
strengths from 80 to 140ksi can be economically used for automotive
structural components, a comprehensive design guide for the use of a
broader range of high strength sheet steels is highly desirable.
Since early 1982, a new research project entitled “Structural Design of
Automotive Structural Components Using High Strength Sheet Steels’ has
been conducted at the University of Missouri—Rolla. The main purpose of
the project is to determine the material characteristics of typical high
strength sheet steels having yield strengths in the range of 80-140 ksi and to
develop the design criteria for the cold formed steel structural members
composed of such high strength steels.
The preliminary study has been completed at the present time. The
structural research will be initiated by a study of the strength of beam webs
and the structural components consisting of curved and flat elements.
The above discussion is related to the UMR research work on cold formed
steel structures. In view of the fact that many research projects on thin-
walled metal construction have been conducted by a large number of
researchers at various institutions and organizations throughout the world,
118 Wei-Wen Yu
4. CONCLUSIONS
Acknowledgement
All of the research projects discussed in this paper have been sponsored by
the American Iron and Steel Institute. The technical guidance provided by
various AISI task groups and the AISI Staff is gratefully acknowledged.
Appreciation is also extended to Steel Deck Institute and H. H. Robertson
Company for cosponsoring the project on steel decks.
UMR Research on Cold Formed Steel Structures 119
REFERENCES
33: BERGFELT, A., Profile’s minces formes a froid, Bulletin Technique De La Suisse
Romande, 99, No. 17 (August 1973) 363-8.
34, HOGLUND, T., ‘Design ofthin plate I girders in shear and bending with special
reference to web buckling’, Bulletin No. 94, Division of Building Statics and
Structural Engineering, Royal Institute of Technology, Sweden, September
1973.
Bo: TuHomasson, P., ‘Livbuckling ens Inverkan pa Burformagan Has Trapets
Profilerad Stalplat’, Nordiske Forskningsdager for Stalkonstruksjoner, Oslo,
Norway, August 1973.
36. WINTER, G. and Pian, R. H. J., Crushing strength of thin steel webs, Cornell
Bulletin, No. 35, Part I, 1946.
37. ZETLIN, Lev, Elastic instability of flat plates subjected to partial edge loads,
Journal of the Structural Division, Proceedings of the American Society ofCivil
Engineers, 81 (September 1955).
38. YAMAKI, N., ‘Buckling of a rectangular plate under locally distributed force
applied on the two opposite edges’, Ist and 2nd report, The Institute of High
Speed Mechanics, Tohoku University of Japan, Vol. 3, 1953.
39. Wuite, R. N. and CottinGHaM, W. S., Stability of plates under partial edge
loadings, Journal of Engineering Mechanics Division, ASCE Proceedings, 88,
No. EMS (October 1962).
40. BaGcul, D. K. and Rockey, K.C.,‘A note on the buckling ofa plate girder web
due to partial edge loadings’, International Association for Bridge and
Structural Engineering, Final Report, September 1968.
41. Rockey, K.C. and BaGcui, D. K., Buckling of plate girder webs under partial
edge loadings, International Journal of Mechanical Science, 12 (1970).
42. Rockey, K. C. and Et-Gaaty, M. A., ‘Ultimate strength of plates when
subjected to in-plane patch loading’, International Association for Bridge and
Structural Engineering, Proceedings of the Seminar on Design of Plate and
Box Girders for Ultimate Strength, 1971.
43. Rockey, K. C., EL-GAALy, M. A. and BAGCHI, D. K., Failure of thin walled
members Gader patch loading, Journal of eral Division, ASCE
Proceedings, 98, No. ST12 (December 1972).
. Rockey, K. C. and Er-Gaaty, M. A., Stability of load bearing trapezoidal
diaphragms, International Association for Bridge and Structural Engineering,
Publications, 32-Il (1972).
45. E.t-GaaLy, M. A. and Rockey, K. C., ‘Ultimate strength of thin walled
members under patch loading and bending’, Proceedings of the Second
Specialty Conference on Cold-Formed Steel Structures, University of
Missouri—Rolla, October 1973.
46. Rockey, K. C. and Et-Gaaty, M. A., ‘Ultimate strength of plates when
subjected to in-plane patch loading’, International Association for Bridge and
Structural Engineering, Colloquium on Design of Plate and Box Girders for
Ultimate Strength, London, England, March 1971; published July 1972.
47. BAEHRE, R., ‘Sheet metal panels for use in building construction—recent
research projects in Sweden’, Proceedings of the Third International Speciality
Conference on Cold-Formed Steel Structures, University of Missouri—Rolla,
24-25 November 1975.
122 Wei-Wen Yu
48. Rat irr, G. D., ‘Interaction of concentrated loads and bending in C-shaped
beams’, Proceedings of the Third Specialty Conference on Cold-Formed Steel
Structures, University of Missouri-Rolla, 24-25 November 1975.
49, Yu, W. W., ‘AISI design criteria for bolted connections’, Proceedings of the
Sixth International Specialty Conference on Cold-Formed Steel Structures,
University of Missouri—Rolla, November 1982, pp. 675-98.
50. WINTER, G., Tests on bolted connections in light gage steel, Journal of the
Structural Division, ASCE Proceedings, 82, No. ST2 (March 1956).
51. WINTER, G., ‘Light gage steel connections with high-strength, high-torqued
bolts’, Publication of the International Association for Bridge and Structural
Engineering, Vol. 16, 1956, p. 513.
52 Porowicn, D. W., ‘Tension capacity of bolted connections in light gage cold-
formed steel’, Thesis, Cornell University, 1969.
Ssh DHALLA, A. K., ERRERA, S. J. and WINTER, G., Connections in thin low-
ductility steels, Journal of the Structural Division, ASCE Proceedings, 97,
No. ST10 (October 1971).
54. CHONG, K. P. and Mattock, R. B., Light gage steel bolted connections
without washers, Journal of the Structural Division, ASCE Proceedings, 101,
No. ST7 (July 1975).
5). GILCHRIST, R. T. and CHONG, K. P., Thin light-gage bolted connections
without washers, Journal of the Structural Division, ASCE Proceedings, 105,
No. ST1 (January 1979).
56. Yu, W. W. and Mossy, R. L., ‘Bolted connections in cold-formed steel
structures’, University of Missouri—Rolla, January 1981.
7, ‘Specification for the design of light gage cold-formed stainless steel structural
members’, American Iron and Steel Institute, 1968 edn.
58. JOHNSON, A. L., ‘The structural performance of austenitic stainless steel
members’, Ph.D. Dissertation, Cornell University, February 1967;
Department of Structural Engineering, Rept. No. 327, Cornell University,
November 1966.
oo: JOHNSON, A. L. and WINTER, G., Behavior of stainless steel columns and
beams, Proceedings, ASCE, 92, No. STS (October 1966) 97-118.
60. Jounson, A. L. and KELSEN, G. A., Stainless steel in structural applications,
Stainless Steel for Architecture, ASTM STP 454, American Society for Testing
and Materials, August 1969.
61. Wana, S. T., ‘Cold-rolled austenitic stainless steel: material properties and
structural performance’, Department of Structural Engineering, Rept. No. 334,
Cornell University, July 1969.
62. WANG, S. T. and Errera, S. J., ‘Behavior of cold-rolled stainless steel
members’, Proceedings of the First Specialty Conference on Cold-Formed Steel
Structures, University of Missouri-Rolla, August 1971.
63. ERRERA, S. J., TANG, B. M. and Popowicn, D. W., ‘Strength of bolted and
welded connections in stainless steel’, Department of Structural Engineering,
Rept. No. 335, Cornell University, August 1970.
64. ‘Stainless Steel Cold-Formed Structural Design Manual’, American Iron and
Steel Institute, 1974.
65. RANG, T. N., GALAMBOs, T. V., YU, W. W. and RAvINDRA, M. K., ‘Load and
resistance factor design of cold-formed steel structural members’, Proceedings
UMR Research on Cold Formed Steel Structures 123
J. RHODES
University of Strathclyde, Glasgow, UK
ABSTRACT
Some of the problems associated with the buckling behaviour of cold formed
steel structural sections are examined. The methods of approach to these
problems used in various design codes are outlined and compared. Effective
widths of compression elements, edge and intermediate stiffener require-
ments for these elements, and some aspects of column and beam behaviour
are considered. Suggestions are made for appropriate design methods to deal
with the problems examined. .
1. INTRODUCTION
In recent years there has been a great deal of effort devoted to the
compilation of new design codes for cold formed steel sections. New AISI!
and Swedish? design codes were introduced in 1980, and substantial work
has been carried out since then on drafting European Recommendations?
for cold formed steel design. At the same time, new national design codes
are in the process of being developed in various countries including the UK.
It is perhaps, therefore, an opportune time to give some thought to the
problems encountered in cold formed steel design and on how these can
best be dealt with. In this paper some impressions of typical aspects of cold
formed section design requirements are presented and some suggestions are
made as to how these aspects can be dealt with in a design code.
125
126 J. Rhodes
In the drafting of a new design code there are factors other than purely
technical details which must be considered. In the case of cold formed steel
design, in which local buckling adds complexity to the various design
problems, one of the most important requirements of a design code is
simplicity of use. There is no merit in providing a designer with
sophisticated analytical methods if he is not able to use them.
Generality of application is also a requirement. Lack of generality is at
present the major failing of the UK code, which users find does not deal
with a number of commonly occurring situations.
If a design code must be simple and general then it cannot always be
accurate in dealing with every situation. In such cases, if simple rules are
inadequate to accurately cover all situations then the range of application
must be known and it must be ensured that the simple rules at least err on
the safe side. Conservatism is therefore a further requirement.
These three requirements are not always easy to satisfy simultaneously,
while at the same time ensuring efficient design. This is especially so when it
is considered that during the lifespan of a code there may be significant
changes in the use of cold formed steel. It is therefore necessary also to
make some attempt at keeping the provisions of a code suitable for dealing
with possible future needs.
With these thoughts in mind, some ideas as to the way 1n which the design
of some elements and members may be dealt with are given in the remainder
of this paper.
is +
' es i} ll
ul II lI
i i Nl
Lee | cov i]
ca) T
Stiffened Element Hu
(ob)
Unstiffened Element
cc) (é
Edoe Stiffened Element aterediote. SiufFenead
Element
Fig. l(a) or on one edge only (unstiffened) as in Fig. 1(b). The elements may
also have lips, or edge stiffeners (c) or intermediate stiffeners (d).
The major factor affecting the behaviour of compression elements is
local buckling. This is taken into account in most design codes using the
concept of ‘effective width’ which admirably describes the effects of local
buckling. Many design codes use the effective width formulation of
Winter° or slight modifications of this, although some codes use other
approaches. The UK specification uses, instead of effective widths, stress
reduction factors (CL factors) based on rigorous analysis by Walker® of a
simply supported square plate with a carefully devised imperfection
parameter which describes imperfections commonly found in cold formed
steel sections.
While the writer feels that an effective width formulation is more suitable
in design than a stress reduction factor approach, the rigorous basis of the
UK specification provides a consistently conservative prediction of plate
strength.
With this in mind, an effective width formulation which gives results in
broad agreement with the UK code predictions, or a little less conservative,
128 J. Rhodes
4-022
ral +14 ([Pau 0,35) | for 2™*
> 0.59 (1)
b Gor Cor
where 8, is the effective width, 5 is the full width and o,,,, is the maximum
compressive stress.
This expression permits evaluation of effective width for any stress level.
The buckling stress o,, is given by
kn’°E py
Per ~ 19(1 — v2) (*) 2)
where the buckling coefficient, k, may be taken as 4 for a stiffened element
and 0.425 for an unstiffened element.
Equation (1) is a little more complex than the effective width expression
of the AISI code, but is a closer fit to the existing UK CL factor curve and
also gives results more consistent with curves derived for heavier structures
in the low b/t range. Figure 2 shows comparison of the effective widths
given by eqn (1) with those of the AISI specification and with the CL factors
of the UK specification. Also shown are the curves of Chatterjee and
Dowling’ and Scheer and Nolke® for bridge structures.
For simplicity in design, eqn (1) may be used for both stiffened and
unstiffened elements using the appropriate & factors. It is known (see ref. 9)
that the behaviour of unstiffened elements can be conservatively predicted
by a formulation set up for stiffened elements (with the appropriate k
factors) so long as the unstiffened elements can be prevented from bending.
However, as these elements have very low bending rigidity in the presence
of local buckling, safeguards should always be built in to take account of
this.
These state simply that to ensure that the stiffened plate behaves as a
stiffened element (rather than as an unstiffened element) the lip width must
be equal to one-fifth of that of the element to be stiffened.
An examination of this problem shows that, while all the specifications
give rules from which the rigidity requirements of edge stiffeners can be
evaluated by consideration of the element in isolation, this is not the case in
reality. Actually the requirements of edge stiffeners are highly dependent
at (<->)
be ‘\ Eq? a1
b
os 0-3
be
b
oe 0-6 PRenax
fe) os J 1-5'Ccr
x AISI Specification
XN
Ok : Eq? ao
we ON
Specification
Oz
_ =
European =
12 Recs. Dn
/
UK
3 ; .
2 Eq“ cad
fe) 20 LO CO 80 x
Fig. 3. Required minimum edge stiffener dimensions.
Future Cold Formed Steel Design Rules 131
by eqn (3) it is not felt desirable to make any reduction to come into line
with the other specifications. It is quite possible that in the future the AISI
code and similar codes will require larger edge stiffeners in this range. The
European Recommendations give larger stiffener requirements than those
of the AISI in this range. These are based on recent research in the USA at
Cornell University’? in which it was shown, on the basis of experiments,
that the AISI requirements in this range are not sufficient. Furthermore, the
tests of ref. 11 were carried out on beam elements in which efforts were
made to ensure that the supported edges of the edge stiffened elements were
kept as close to the fixed condition as possible. There are many situations in
which edge stiffened elements will suffer more severe conditions so that in
the writer’s opinion even the increased rigidities required by this research
will be inadequate in some circumstances.
The requirements given by eqn (4) are shown in Fig. 4, and these are
substantially the same as those of ref.3, but provide a little more
conservatism for sub-elements of high w/t ratios beyond the limit of the test
results.
Although design codes tend to’ give minimum stiffener rigidity
requirements in terms of the Jc, of the stiffener, in actuality stiffeners bend
132 J. Rhodes
log
TA 300
300
about an axis closer to the plate middle surface than to their own neutral
axis. In order to ascertain more accurately the relative merits of different
stiffener shapes it is therefore more realistic to specify the rigidity
requirements in terms of the second moment of area about the plate middle
surface, say J,,;,. Since the ratio of J,,;, to [¢g varies from about 4 to 5 for
common stiffener shapes, a suitable value for required rigidity based on
eqn (4) is
| bk Wie 6
r= 0.2(—] — 5
t* (*) 280 ©)
For comparison purposes Fig. 5 shows the variation of required height, h
of a stiffener of equal height and width / obtained from eqn (5), ref. 3 and
the AISI and CSA’? specifications, with variation in sub-element w/t ratio.
td
12
JO
|=
4. COLUMNS
{20
Ah) Ref Ce) =3%
alo
Failouve ay oN ALSI1
Stress s ™~ fas pe
CN 2,00 > oe ,
/vamn) Cc) R
Suggested ce : oe,
approach ie mage
30
1.6
,No warping restrawmt
Partial Warping
restraint
J+
Full warping
/ restraint
Sa
oO 25 5 Us 10 12°5
5. BEAMS
In the past few years there has been a number of developments in the design
analysis of cold formed beams which should be taken into account in any
new code. A few of these developments are discussed here.
20 25 30 4O
In the case of sections bent about an axis of symmetry and having low
shape factors, the plastic capacity so evaluated will not be a great deal more
than the elastic capacity, but in other cases significant increases in
calculated capacity can be obtained. An indication of the increased
capacity is illustrated in Fig. 9 which plots the ratio of plastic and elastic
moment capacities given by the above suggestions for a hat section beam
with compression flange b/t less than 25,/280/a,.
It is suggested that if a section is compact, as described above, plastic
limit analysis may be used for continuous beams. If the section is not fully
compact then the only advantage which is intended to be obtained using the
plastic capacity given above is the increase in moment capacity at a point,
and limit analysis should not be used for continuous beams.
a P Saal compression
tension
where b,/b is obtained from eqn (1) using the web 5/t ratio and the k factor
corresponding to the compressive stress at buckling of the web element.
This factor for a stiffened element can be accurately and conservatively
calculated using the expression
Curve 3 takes compression flange buckling into account but neglects web
buckling, and thus illustrates the significant effects of web buckling.
Curve 4 shows the predictions of failure using the AISI maximum bending
stress (with safety factors omitted) together with the compression flange
effective width determined on the basis of eqn (1). The predictions of this
approach are very good, although non-conservative, and the sheer
25
2
CNo buckling considered) 3
CWel buckling
els regarded)
Mexe
5
CSuggested
o— Experimental results from approach)
[eXe) Various Seurces given mM
Ref CiG)
wax
oO 2 Ly @ 3 JO Ger
6. CONCLUDING COMMENTS
This paper has only considered a few of the topics which must be taken care
of in future design codes. Many other topics such as beam deflections,
stiffeners for bending elements, the use of enhanced k factors when these
are available, are of importance and should be part of a future design code.
Various aspects of cold formed sections behaviour other than buckling
are also important, for example connection behaviour, and these must also
be included, as well as guidance on the use of testing as a basis for design.
Examination of the various design codes indicates that the current UK
code is deficient in the range of design problems it can deal with. The AISI
specification is impressive and deals with a wide range of aspects of cold
formed section behaviour with simple, accurate analysis methods. This
code forms the basis of many other National Specifications and must be
considered carefully in the setting up of new design rules.
Finally, it may be remarked that both tabular and graphical presentation
of design information are very useful, and should be included where
possible in design codes. However, these should not be used as the sole basis
for design analysis, as has been the case, to some extent, in the past. The
availability and widespread use of cheap computational facilities makes it
imperative that design codes become primarily formula orientated for ease
in computerisation. Indeed, considering the advances in microcomputers
in recent years, and the likely continuation of these advances, one may well
ask the question: should the next generation of design codes, say in the
early 1990s, be written as microcomputer packages?
142 J. Rhodes
REFERENCES
. American Iron and Steel Institute, ‘Specification for the design of cold formed
steel structural members’, 1980.
. National Swedish Committee on Regulations for Steel Structures, ‘Swedish
code for light-gauge metal structures’, English translation, 1982.
. ECCS, ‘European recommendations for the design of profiled sheeting and
sections Part 2—Sections’, 1982.
. British Standards Institution, ‘Specification for the use of cold formed steel
section in building: addendum No. | to BS 449: Part 2: 1969’, 1975.
. WINTER, G., ‘Performance of thin steel compression flanges’, Cornell
University Engineering Experiment Station, Reprint No. 33, 1950.
. WALKER, A. C., The post-buckling behaviour of simply supported square
plates, Aero Quarterly, XX (1969).
. CHATTERJEE, S. and DOWLING, P. J., The design of box girder compression
flanges, Steel plated structures, P.J. Dowling, J.E. Harding and P. A. Frieze
(eds), London, Crosby Lockwood Staples, 1977.
. SCHEER, J. and Novke, H., The background to the future German plate
buckling design rules, Stee/ plated structures, P. J. Dowling, J. E. Harding and
P.A. Frieze (eds), London, Crosby Lockwood Staples, 1977.
. Ruopes, J., Effective widths in plate buckling, Developments in thin-walled
structures—1, J. Rhodes and A.C. Walker (eds), London, Applied Science
Publishers, 1981.
. Ruopes, J., ‘Buckling and failure of edge stiffened plates’, IUTAM Symposium
on Collapse, 1982.
. DESMOND, T. P., PEKOz, T. and WINTER, G., Edge stiffeners for thin-walled
members, J. Structural Div. ASCE, 107, No. ST2 (1981).
. DESMOND, T. P., PEKoz, T. and WINTER, A., Intermediate stiffeners for thin-
walled members, J. Structural Div. ASCE, 107, No. ST4 (1981).
. Canadian Standards Association, ‘CSA Standard $136—1974, cold formed
steel structural members’, 1974.
. Ruopes, J. and Harvey, J. M., ‘Interaction behaviour of plain channel columns
under concentric or eccentric loading’, Preliminary Report, Second
International Colloquium on Stability of Steel Structures, Liege, 1977.
. SCHUSTER, R. M., Cold formed steel design manual, Solid Mechanics Division,
University of Waterloo Press, 1975.
. Ruopes, J., ‘The post buckling behaviour of bending elements’, Proc. 6th
International Specialty Conference on Cold Formed Steel Structures, St Louis,
1982.
. Standards Association of Australia, ‘Rules for the use of cold formed steel in
structures AS1538— 1974’, 1974.
Design of Profiled Steel Sheeting and Decking
ABSTRACT
Figures are given for the increase in production of coated steel strip and the
economic importance of cold rolled steel products in industrial buildings. It
is shown how design methods for profiled steel sheeting have evolved and the
present position of British, European and International Codes is
summarised. Factors influencing the use of profiled sheeting in roofing,
decking and wall cladding are considered. Calculation procedures are given
for the strength and deflection of unstiffened and stiffened profiles, according
to the European Recommendations. It is advantageous for some of the
design assumptions to be based on test data—a procedure only applicable to
mass produced structures—and detailed information is given on the testing
of profiled sheeting. The results of many such tests are correlated with
theory. Finally, some factors influencing future trends in the design of
sheeting, decking and sandwich panels are discussed.
1. INTRODUCTION
Over the last four years, UK steel production in many areas has fallen by
more than 30°. The only area in which there has been increase in
production has been in coated steel strip.
In Europe, production of coated material increased from 284 million
square metres in 1980 to 323.4 million square metres in 1981, and there is no
143
144 E. R. Bryan and J. M. Davies
reason to believe that this trend will not continue. For the years quoted the
increase was 16 °% for coated steel and 9.3 °% for coated aluminium.’
Much of the output of coated strip is used in profiled sheeting and
decking, and there has been a dramatic increase in such usage in recent
years. The trend is most evident in Sweden, where the annual consumption
is approaching 3 square metres per capita” but other European countries
are now experiencing the same phenomenon.
It is of interest too to note the results of a survey of costs of four
industrial buildings made by two of the country’s largest steelwork
contractors.* In three of the buildings, the cost of the cold rolled steel
members and sheeting in the walls and roof exceeded the cost of the hot
rolled steelwork. Yet it is safe to assume that the latter accounted for nearly
all the design effort.
It follows that overall efficient structural design can only be achieved
when full consideration is given to the design of the cladding and its
structural implications. The design of the frame may be more important
from the point of view of safety but less important from the points of view
of serviceability, economics and overall aesthetics.
The design approach to mass produced components such as sheeting and
decking may be very different from that used in the design of individual
structures.* With the former it may well be possible to base the design
assumptions on observed behaviour, but such a procedure would rarely be
economic for an individually designed structure. Moreover, the brevity of
the design process for sheeting is not usually of prime importance since it is
worth spending a considerable time in evolving the most economic shape.
Nevertheless, it is of great help to those who check the design if the method
can be readily followed by the non-specialist.
It is against this background that the design philosophy and procedures
for profiled sheeting and decking are developed in this paper.
2. PRESENT POSITION
For many years after the war, sheeting and decking were designed either on
a semi-experimental basis or according to the principles set out in the AISI
specification for the design of cold formed steel structural members.* In
Britain, the only design guide which has been available is that published by
the Metal Roof Deck Association.® In addition, there are two relevant
Codes of Practice,”’*? but these are orientated towards installation and
performance rather than actual design. Nevertheless the upsurge of interest
Design of Profiled Steel Sheeting and Decking 145
Part 4 Code of practice for design of floors with profiled steel sheeting;
Part 5 Code of practice for design of cold formed sections;
Part 6 Code of practice for design of light gauge sheeting, decking and
cladding;
Part 7 Specification for materials and workmanship: cold formed
sections and sheeting;
Part 9 Code of practice for stressed skin design.
Part 4 has just been published‘? and includes a design method for profiled
steel sheeting which is more advanced than any previous method published
in Britain, and which owes a considerable amount to the Swedish work. In
due course the method will be further revised to bring it in line with the
latest procedures proposed for Part 6.
In Europe, Preliminary European Recommendations for the design of
profiled sheeting have been published by ECCS‘? and the revised version
will shortly be available. The final document will incorporate the European
Recommendations for the testing of profiled metal sheets, published by
ECCS in 1977.'* In addition, the European Recommendations for good
practice in steel cladding and roofing'? specify some structural perfor-
mance criteria. It is expected that these documents, which have benefited
from CIB and US participation will be largely accepted as the ISO
Standard in this field.
3. ROOF SHEETING
Profiled steel sheeting used as the outer waterproof skin of a roof, with the
insulation placed inside, is known as ‘cold roof’ construction. In this
146 E. R. Bryan and J. M. Davies
Fig. 2(b) or (c) are often better able to meet this condition. Work is
presently in progress to define the ‘walkability’ of sheeting more objectively,
since the current position depends far too much on subjective judgements.
Because roof sheeting is usually continuous over several spans, the
maximum span is usually determined by strength rather than deflection.
Wind suction, particularly in local areas such as the eaves, gable or apex, is
likely to be a more important design case than dead plus imposed load, and
due attention must be given to the design of fixings.
If long span sheets are required, decking profiles used in the inverted
position (to improve water shedding) may always be used.
4. ROOF DECKING
be used. The profiles shown in Fig. 3, of depth about 100mm and yield
stress 350 N/mm7?, can be used for spans between 6 and 8 m for UK loading.
Since the permissible deflection of decking is limited to span/250,° the
design criterion is usually that of stiffness rather than strength.
Simple trapezoidal profiles have usually been designed according to the
MRDA code? but neither this document nor alternatives readily available
in the UK envisaged the complicated shapes shown in Fig. 3. Such profiles
can only be calculated according to the Swedish code'! or the European
Recommendations. '? In addition, the maximum trough width requirement
given in the MRDA code was related to the thickness of insulation used at
that time and so it may not be appropriate to the much greater thicknesses
used today.
The most recent development in long span decking profiles is the 206 mm
deep cassette profile, stiffened in both directions, as shown in Fig. 4. It can
span up to 12m and can thus span directly between main frames without
the need for secondary members. An assessment of the design has been
made at the Royal Institute of Technology, Stockholm and at the University
of Salford.'® It was first used in Britain in a large building for the
Warrington Development Corporation. In Sweden it is often used with an
additional metal skin to replace the traditional felt and bitumen finish.
5. WALL CLADDING
For many applications, the profiles used for roof sheeting are also suitable
for vertical cladding. In many respects the loading is similar except that the
requirement for the profile to sustain a concentrated point load no longer
applies.
Since the walls are often the most prominent feature of a building, the
choice of wall cladding is often determined more by aesthetic con-
siderations than structural performance (Fig. 5). In fact, a number of
profiles have been specially developed to give different impressions of scale
and light and shade (Fig. 6). In addition, the sheeting may span vertically or
horizontally or even diagonally (Fig. 7). Curved sheeting around the apex
or eaves of a building may also be used to good effect (Fig. 8). The choice of
colour and texture of wall cladding is probably even more important than
for roof sheeting.
A review of the whole field of sheeting and cladding is given in ref. 19.
6. CALCULATION PROCEDURES
This section describes briefly the analytical model that forms the basis of
the Swedish code’? and the European Recommendations.!? A full account
of the derivation of the design expressions has been given by Hoglund.”°
where @, is the stress at the level of the compression flange in the effective
cross-section.
Using the doubly reduced cross-section, the new section properties are
calculated so that, for first yield in the compression flange as shown in
Fig. 9(b), the ultimate moment of resistance M@, may be found.
In the above calculation it can be shown that, for reasonably symmetrical
profiles, the web is fully effective unless the ratio of profile height to
thickness is above about 70. For shallow trapezoidal profiles, therefore,
with bends of small radius, the calculated bending strength is almost
identical to that given in the MRDA code.
The European Recommendations also make provision for an increase of
bending capacity to be allowed when first yield occurs in the tension flange
as shown in Fig. 9(c). The tension flange is allowed to become plastic and
the ultimate moment of resistance M, is calculated when the compression
flange also reaches yield. If advantage is taken of this concession, a, is
always equal to the yield stress of the material, o,.
For continuous sheeting, where the bending moment changes sign over
the support, the above calculation has to be repeated in the case of
unsymmetrical sheeting. The ultimate moment of resistance over the
support will not generally be the same as that near midspan.
Unlike the MRDA code, the European Recommendations also require
that the support reaction capacity R, is calculated and that account is taken
154 E. R. Bryan and J. M. Davies
0:25
M_ when Ra <0:25,
Mg
M
then Mg < 1:0
ies —0:52
2 en R_ >0:25 <10,
d
M R
then Ma ae Oe eee
6)
10) OFS Yip
1-0 aabe
Fig. 10. Interaction diagram for permissible values of ultimate moment M and
support reaction R at an internal support.
Design of Profiled Steel Sheeting and Decking 155
Y [\Sen=t 5 Pa
NA
be obtained by dividing the member into slices and using either computer
analysis or another appropriate numerical procedure. However, this is
rather impractical for general use and a suitable value of /, for practical
purposes can be obtained by calculating a single value at a bending moment
averaged over the entire length.
k
1 x
' Side 2 A
| : ! :
oe Side 2
Conditions at central support above are
equivalent to conditions under point load: — - Side14
O4L
The tests illustrated in Figs 12-16 were carried out at the University of
Salford during a comprehensive assessment of the decking profiles
concerned. The results confirmed the design principles of the Swedish
160 E. R. Bryan and J. M. Davies
Fig. 17. Interaction between bending moment and support reaction for 1500
trapezoidal profiles.
Design of Profiled Steel Sheeting and Decking 161
The results of Baehre and Fick’s tests as well as the Swedish tests are plotted
in Fig. 17 together with several interaction formulae suggested as lower
bounds. Of these A is the expression given above (and as presently included
in the Preliminary European Recommendations), B is the more
conservative linear expression:
and Cisa parabolic interaction curve favoured by the German DIN 18807.
It is probable that the final European Recommendations will adopt
expression B. In fact, because in practice the values of R/R, are rather
small, the important part of the interaction curves is the upper part where
the differences are quite small.
The present state of the art thus requires elastic or quasi-elastic design of
continuous sheeting based on the calculated moment of resistance at the
internal supports. Figure 17 shows that this can in itself be very
conservative because many of the test points lie well above the lines
representing the alternative design expressions. When the results of tests
on two-span sheeting are compared with the above theory a further degree
of conservatism is often revealed. Most sheeting and decking profiles are
capable of considerable redistribution of bending moment so that the
attainment of the calculated moment of resistance at an internal support is
162 E. R. Bryan and J. M. Davies
9. FUTURE TRENDS
spans. To rectify this deficiency, ECCS Committee TC7 has set up a Study
Group to make proposals for eventual European Recommendations on
Sandwich Panels. It is hoped that this document will serve as a base for all
national codes and will stimulate further the use of sheet steel in building.
REFERENCES
. Makins, N., Coil coating comes into its own, Financial Times, 20 May 1982.
No. BAEHRE, R., ‘Sheet metal panels for use in building construction’, Proceedings
third international specialty conference on cold formed steel structures,
University of Missouri-Rolla, November 1975.
. BRYAN, E. R., ‘European Recommendations for cold formed sheet steel in
building’, Proceedings fifth international specialty conference on cold formed
steel structures, University of Missouri-Rolla, November 1980.
. University of Salford, Conference on light gauge steel, 20 May 1981, p. 34.
. American Iron and Steel Institute, ‘Specification for the design of cold formed
steel structural members’, first edition April 1946, latest edition 3 September
1980.
. Metal Roof Deck Association, ‘Code of design and technical requirements for
light gauge metal roof decks’, first published 1965, third edition 1970.
. British Standards Institution, ‘Code of practice for sheet roof and wall
coverings’, CP 143: Part 10: 1973.
. British Standards Institution, ‘Code of practice for performance and loading
criteria for profiled sheeting in building’, BS 5427: 1976.
. ‘Profiled steel cladding and decking for commercial and industrial buildings’,
London, Constrado, 1980.
. ‘Stahltrapezprofil im Hochbau’, IFBS, Stuttgart, Karl Kramer Verlag, 1980.
. Swedish Institute of Steel Construction, ‘Swedish code for light-gauge metal
structures’, Publication 76, March 1982 (in English).
. British Standards Institution, ‘Structural use of steelwork in building: code of
practice for design of floors with profiled steel sheeting’, BS 5950: Part 4: 1982.
. European Convention for Constructional Steelwork, ‘Preliminary European
Recommendations for the design of profiled sheeting and sections: Part 1—
profiled sheeting’, ECCS-T7-1981, March 1981.
. European Convention for Constructional Steelwork, ‘European Recommenda-
tions for the testing of profiled metal sheets’, ECCS-XVII-77-2E, April
1977.
. European Convention for Constructional Steelwork, ‘Preliminary European
Recommendations for good practice in steel cladding and roofing’, ECCS-T7-
1981, December 1981.
. ANDREW, T., Cladding for life, Building with steel, December 1979. See also
other articles in this issue and in that of July 1979.
. Tarmac Building Products Ltd, Flat roofing—a guide to good practice, RIBA,
1982.
. BRYAN, E. R., Long-span steel decking, The Structural Engineer, 57A, No. 4
(April 1979).
164 E. R. Bryan and J. M. Davies
Sy, Davies, J. M., Profiled steel cladding, The Consulting Engineer (February
1980).
20. HOGLUND, T., ‘Design of trapezoidal sheeting provided with stiffeners in the
flanges and webs’, Swedish Council for Building Research, Document D28:
1980 (in English).
DAN. BAEHRE, R. and Fick, K., ‘Berechnung und Bemessung von Trapezprofilen—
mit Erlatrungen zur DIN 18807’, Berichte der Versuchanstalt fiir Stahl, Holz
und Steine der Universitat Fridericiana, Karlsruhe, 1982.
22) Lawson, R. M. and CRIsINEL, M., ‘Evaluation of the properties of profiled
sheeting’, private communication, 1982.
10
I. H. MARSHALL
Paisley College of Technology, UK
and
W. M. BANKS
University of Strathclyde, Glasgow, UK
ABSTRACT
NOTATION
The buckling of thin, isotropic, homogeneous, flat plates has been the
subject of a vast amount of previously published work. Papers on this
subject have appeared in the technical press for approximately a century
with in excess of two thousand publications presently available. With so
many references available for consultation the researcher could be forgiven
for being not only overawed by the quantity of previous studies but also
slightly bewildered by the diversity of topics considered and the variety of
analytical approaches utilised. Indeed the diversity of past analytical
approaches used in the solution of seemingly identical problems appears to
leave the investigator more than sufficiently armed to solve the particular
problem in hand. Unfortunately such thoughts are often misleading.
From the classical paper by Bryan! almost a century ago, through
notable landmarks? ~* to the present date many analytical techniques have
been successfully employed. Apart from classical closed-form solutions for
specific problems (relatively few in number) approximate techniques
utilising energy, semi-energy, finite difference, finite element, etc. have
found favour.
Although ‘approximate’, such solutions can converge rapidly to yield
realistically accurate buckling and post-buckling characteristics up to two
or three times the critical load. In general the accuracy of these solutions
depends on the number of terms utilised in mathematically modelling the
problem, or the mesh size employed in the case of finite element techniques.
In the case of semi-energy solutions the earliest recorded example used to
determine plate post-buckling characteristics was in 1937 in a paper by
Marguerre.® Essentially the technique involved postulating an expression
to specify the plate lateral deflections (w) thereby allowing an exact solution
to the Von Karman compatibility equation to be obtained. Having
determined an appropriate relationship between plate middle surface
Instability of Composite Structures 167
Providing the in-plane forces N, and N, are constant and the in-plane
shear loading is zero an exact solution to (1) can be found. Typically this
solution is in the form of an eigenvalue (defining the buckling load) and an
eigenvector (defining the mode of buckling).
The Von Karman-type compatibility equation for thin, flat, orthotropic
plates in the presence of large deformations can be written in the form
Ue Ee O*F - 1 O¢F (0?w\? 0?w d?w
Eos Ox oy? ES Gy Shdx oy dx? dy? (2)
Simultaneous solution of eqns (1) and (2) would yield an exact, closed form
solution for the buckling and post-buckling characteristics of the problem
under consideration. However, such solutions are rare and generally
limited to the most straightforward of plate boundary conditions.
Commonly, solutions represent the plate deformations and stresses as
trigonometric Fourier series, thereby producing a set of simultaneous non-
linear algebraic equations equal in number to the chosen number of terms
in the respective series. Mathematically such solutions are ‘exact in the
limit’, 1.e. an ‘exact’ solution only results if an infinite number of terms in the
Fourier series are used. However, since an iterative technique is usually
involved in the solution, numerical approximations are unavoidable.
As an alternative to solving eqns (1) and (2) simultaneously minimum
energy principles can be employed in lieu of equilibrium considerations.
Thus, after satisfaction of compatibility and constitutive requirements is
obtained by solution of eqn (2) the total energy of the system is minimised
according to the semi-energy approach to yield on upper bound solution to
the problem. Providing realistic functions are chosen when defining the
plate lateral deflections a satisfactory solution with acceptable accuracy
will result. The Galerkin technique has also been successfully employed to
obtain an approximate solution to eqn (1) along with the aforementioned
minimum energy principles or occasionally a ‘double’ Galerkin solution to
eqns (1) and (2). In this latter case it is not possible to bound the resulting
solution.
initial curvature specified using a trigonometric series (Fig. 1). Other cases
of hyperbolic, parabolic etc. curvature have received scant attention.
w = AY cos ae (4)
eb
and wy
Fig. 3. Load-end displacement curves for square SS GRP plates for various initial
imperfections.
172 I. H. Marshall and W. M. Banks
Fig. 4. Relationship between post-buckling load and central deflection for square
SS GRP plate with various imperfections.
marked around the buckling load of the plate. A similar set of results to this
figure is obtained by plotting the post-buckling load against the maximum
membrane stress in the plate. It is thus evident that while the magnitude of
the stresses is significantly affected around the buckling load, they are
virtually unaffected at the higher loads.
The same relative effect is observed for the out-of-plane deflections
presented in Fig. 4.
Banks and Harvey**> found that for load out-of-plane deflection
behaviour for orthotropic plates the deviation of experimental results from
those for a perfect plate increased as the aspect ratio of the plates increased.
This was considered to be due mainly to the difficulty in fabricating the
plates and the resulting imperfections obtained. This initial imperfection
will of course occur in practical applications for composite structures. A
typical load out-of-plane deflection curve for a GRP plate is shown in
Fig. 5. The experimental points are seen to tie in with a theoretical
prediction of an imperfection of half the plate thickness.
Stress distributions for this same plate are seen in Fig. 6. Again it can be
seen that significant variation from the perfect plate obtains for the
logitudinal membrane and bending stresses around the buckling load.
Imperfections are obviously an important governing factor on stresses in
this region of loading.
The isotropic problems of Rushton'*:!* have been extended to include
orthotropic material properties in recent publications,*°3° viz. the
snap-buckling of imperfect plates loaded on the convex surface.
Instability of Composite Structures 173
Fer
The latter two terms in eqn (6) are those necessary to define the initially
deflected shape, i.e. Fig. 2. Thus, a solution of eqn (6) along with any of the
aforementioned analytical techniques would -yield the relevant load/
deformation characteristics. The case of constant curvature plates subject
to in-plane compression appears to have received scant attention in the
XO. EXPERIMENTAL
Og8 EXPERIMENTAL
CLAMPED
SUPPORTS
@ LOADING
0 UNLOADING
cl :
SIMPLE SUPPORTS
nz\i1300mm
weh
Fig. 8. Lateral pressure-out of plane deflection relationship for constant
curvature plate.
3. CONCLUSIONS
REFERENCES
. BRYAN, G. H., On the stability of a plane plate under thrusts in its own plane,
with applications to the ‘buckling’ of the sides of a ship, Proc. London Math.
Soc., 22 (1891) 54-67.
. Cox, H. L., ‘The buckling of thin plates in compression’, Tech. Report of the
Aero Committee, 1933-1934, pp. 443-63.
. Levy, S., ‘Bending of rectangular plates with large deflections’, NACA Report
No. 737, 1942, pp. 139-57.
. Korter, W. T., ‘The effective width of infinitely long, flat rectangular plates
under various conditions of edge restraint’, National Luchtvaartlaboratorium,
Rep. S287, December 1943.
. STEIN, M., ‘Loads and deformations of buckled rectangular plates’, NASA TR
R-40, 1959.
. MARGUERRE, K., ‘The apparent width of the plate in compression’", NACA TM
No. 833, July 1937.
. DuNCcAN, W. J., ‘The principles of the Galerkin method’, ARC Tech. Report,
1938.
. YAMAKI, N., Postbuckling behavior of rectangular plates with small initial
curvature loaded in edge compression, J. App/. Mech., Trans. Am. Soc. Mech.
Engrs., 26 (1959) 407-14.
. YAMAKI, N., Postbuckling behavior of rectangular plates with small initial
curvature loaded in edge compression—(continued), J. Appl. Mech., Trans.
Am. Soc. Mech. Engrs., 27 (1960) 335-42.
. YAMAKI, N., Experiments on the postbuckling behavior of square plates loaded
in edge compression, J. Appl. Mech., Trans. Am. Soc. Mech. Engrs., 28 (1961)
238-44.
. Ruopes, J. and Harvey, J. M., Examination of plate post-buckling behavior,
J. Mech. Eng. Sci., ASCE, 103, No. EM3 (1977) 461-78.
. RHopEs, J., HARvEy, J. M. and Fox, W. C., The load-carrying capacity of
initially imperfect eccentrically loaded plates, Int. J. of Mech. Sci., 17 (1975)
161-75.
. RUSHTON, K., Large deflexion of plates with initial curvature, Int. J. Mech.
Sci., 12 (1970) 1037-51.
. RusHTOoN, K., Buckling of laterally loaded plates having initial curvature, Jnt.
J. Mech. Sci., 14 (1972) 667-80.
. Letssa, A. W., ‘Advances and trends in plate buckling research’, Proc. Symp.
on Advances and Trends in Struct. and Solid Mech., Washington, USA,
October 1982.
. LEKHNITSKI, S. G., Anisotropic plates, English translation by S. W. Tsai and T.
Cheron, J. E. Ashton (ed.), Technomic Pub. Co., USA, 1970, pp. 248.
. AMBARTSUMYAN, S. A., Theory of anisotropic plates, English translation by T.
Cheron, J. E. Ashton (ed.), Technomic Pub. Co., USA, 1970.
. ASHTON, J. E. and Wuitney, J. M., Theory of laminated plates, Technomic
Pub. Co., USA, 1970.
19. Jones, R. M., Mechanics of composite materials, Scripta Book Co., USA, 1975.
20. VINSON, J. R. and CHou, T. W., Composite materials and their use in structures,
London, Applied Science Publishers, 1975.
Instability of Composite Structures 177
G. D. GALLETLY
University of Liverpool, UK
ABSTRACT
The background to design proposals currently being considered by the BSI
for two shell buckling problems is discussed in the paper. The problems are
(a) the buckling of torispherical shells subjected to internal pressure and (b)
the buckling of cylindrical shells under the combined action of external
pressure and axial compression. The shells considered are those occurring in
practice rather than theoretically perfect ones, i.e. they are likely to have
initial shape imperfections and contain residual stresses. Some comments on
further work which is needed are also given.
NOTATION
Amplitude of wo.
Perpendicular distance from the axis of revolution to the centre of
curvature of the torus.
Constants—see eqns (7) and (8).
Diameter of the cylinder to which the torispherical shell is
attached.
Modulus of elasticity.
Constants—see eqns (9)-(12).
Length of cylindrical shell (combined loading case).
Circumferential direct stress resultant.
(i) Internal pressure in a torispherical head, or (ii) external pressure
at which buckling occurs under the combined loading.
179
180 G. D. Galletly
1. INTRODUCTION
The pressure vessel ‘Design Methods’ subcommittee of the BSI (i.e. PVE-/5)
deals with many design problems relating to shells. Two of the current
Design Proposals for Shell Buckling Problems 181
Very thin torispherical shells made from stainless steel and having diameter
thickness (D/t) ratios in the range 700-1300 are frequently used in the food
processing and brewing industries; the shells are also used, but to a much
smaller extent, in the nuclear industry, e.g. in the primary tanks of some
fast reactors and in the roof structures of containment buildings in PWRs.
Despite the fact that many of these thin torispherical shells have been
manufactured (and have given satisfactory service), there are no rules
covering their design in BS 5500. This is because there is, at present, a
restriction in Section 3.5 of the Code that D/t should not exceed 500. The
reason for this limit is not entirely clear but, in part, its purpose must be to
minimise the occurrence ofinternal pressure buckling in the knuckle region
of the heads.
During the last ten years there have been many experiments conducted
on internally pressurised torispheres. As-manufactured stainless steel
models up to 3.6m diameter were tested by Kemper! and by Stanley and
Campbell.* Small (0.13 m dia.) machined models of aluminium were tested
by Gill and his co-workers** and by Galletly>® (the latter also tested
some carbon steel models). The last ten years has also seen the emergence of
a powerful elastic-plastic shell of revolution buckling program known as
Bosor 5.’ With this program it has been possible to analyse the internal
pressure buckling problem under discussion, at least for nominally perfect
shells.2~1° Large-deflections need to be taken into account in the
theoretical treatment and the relevant failure mode, for steel shells, is
frequently asymmetric plastic buckling, with buckles occurring around the
circumference (see Fig. 1).
Using Bosor 5 to conduct parametric surveys, on ellipsoidal and
182 G. D. Galletly
Just over 25 years ago, a 15m diameter fluid coker with a torispherical
bottom head failed while undergoing its hydrostatic proof test. A
subsequent elastic analysis of the head by Galletly!’:'® showed that the
direct hoop stress in the knuckle was compressive and exceeded the yield
point of the material at a number of locations. As the design seemed
questionable, it was decided to have a limit analysis of the head carried out.
This was done by Drucker and Shield and the results of the study provided
the stimulus for their subsequent work on the subject.!®:?° In essence, their
results confirmed the predictions of the elastic analysis with regard to stress
levels.
Design Proposals for Shell Buckling Problems 183
a
Ta
Toroidal shell
It may seem strange that the knuckle (or toroidal portion) of the
torisphere should have compressive stresses in it when the loading is
internal pressure. If the pressure vessel is divided into its component parts
(cylinder, spherical cap and toroidal segment) then the first two will move
outwards when subjected to internal pressure whereas the toroidal part will
move inwards. This latter can be shown from the membrane theory of
shells,*? from which one finds for the hoop direct stress resultant, No:
Since a/r is usually greater than 1.5, the ratio R,/R, will be greater than
two for the toroidal segments occurring in practice.
Since the hoop stress resultants are compressive in the toroidal region, it
means that buckling of the knuckle in the circumferential direction is a
possibility.* Whether it will actually occur depends on several factors, e.g.
* See Fig. 8 in ref. 13 for a simple model for demonstrating internal pressure
buckling.
184 G. D. Galletly
= Asymmetric
One 60000 Plastic Buckling
\
200r . x
Fe Asymmetric
& . Elastic Buckling
Fer 50+ <Soe ~
( bf /in* ) = SS
yp= 30000.
0 Ibf/in
500 1000 1500 2000
Fig. 3. Curves ofp,, and p, for a steel torispherical shell and two values of ¢,,.
have shown that the use of the large-deflection theory of shells, and
considering an elastic-plastic analysis, often results in higher collapse
pressures than those given by the limit analysis.?”
Both circumferential buckling (with waxes or wrinkles in the hoop
direction) and the large deformation axisymmetric mode (obtained from a
limit analysis or otherwise) have to be considered in the design of a head.
Which one occurs in practice depends on the geometric characteristics and
the material properties of the head. This point is illustrated in Fig. 3. While
buckling is the main concern in this paper, it should be realised that it is not
the controlling failure mode for the thicker heads. Also, while elastic
buckling may not be of concern for steel heads made from the normal
steels, it is of interest for heads made from plastic materials (PVC, etc.) and
sometimes aluminium.
Design Proposals for Shell Buckling Problems 185
The sequence of events leading to the proposed design rules was more or
less as follows.
(1) The tests of Galletly and Lagae on small machined models were
analysed using the Bosor 5 shell buckling computer program.*’®’?°
The actual stress-strain properties of the shell wall material were
used in the analyses, as well as the meridional variations in
thickness. The wall thickness of the models in the hoop direction
was not constant but, for the D/t = 500 models, the variation was
not large and it was ignored in the analyses.
In general, the correlation between theory and experiment was
reasonably good. However, the definition of when buckling occurs
186 G. D. Galletly
150
Stress
(Ibf /in? )
+10?
100
50
ed
0. le2igea5 10 5 20
Strain (in/,_x103)
(a) The stress-strain curve B in Fig. 4 was used throughout, for all
the heads. Actually, it should only be employed for the toroidal
region of the heads; however, as this region is where buckling
occurs, there should be no great error in using it for the whole
head.
(b) The meridional variations of thickness were taken into
account.
(c) Residual stresses were ignored.
(d) The nominal radii of curvature were used in the analyses.
__ 285[1 — 1250,,/E](r/D)°-*4
Per! Oyp = (Dif *2(RID Ys (5)
When the steels normally used in dished end construction are considered,
an average value of the square-bracketed term in eqn (5) is approximately
0.8. Substituting this value in eqn (5) leads to
230)
(6)
Paley ea
(Dj)! 53(R/D)"
As mentioned earlier, the value of o,,, for austenitic stainless steel is taken as
the 0.2% proof stress of the ‘as-received’ plate material for crown and
segment ends. For cold spun ends, which experience thinning during the
forming operation, a factor of 1.5 is applied to the foregoing 0.2 % proof
stress.
As manufactured torispherical shells contain residual stresses (due to
welding and forming) and shape imperfections. To account for these and
other items, the theoretical p,, values given by eqn (6) have to be divided by
certain factors, C, to bring them into agreement with the experimental
buckling results. The value of C depends on the method used to
manufacture the ends and is as follows:
ee e for cold spun ends (7)
2.3 for crown and segment ends
The design pressure, p, for a cold formed torispherical shell is then found
from
P=P./C (8)
and this equation gives a safety factor of about 1.5.
Design Proposals for Shell Buckling Problems 191
P Ginyes
0.2% proof stress ae (Ditech) 0)
where
oe 100 for crown and segment ends (10)
203 for cold spun ends
and the 0.2 % proof stress is that of the as-received plate. Equation (10) is
useful for determining the design pressure of a torispherical head whose
dimensions are known. Inversion of eqn (10) gives
0.65 R/D Unne4
where
0.050 for crown and segment ends
Kee 0.032 for cold spun ends
(12)
The thickness of the knuckle region can be calculated directly from eqn (11)
if the pressure, etc. are known. The form of eqn (11) is very similar to that
used in the 1982 edition of the French Pressure Vessel Code CODAP??
for this problem. However, CODAP does not appear to cite the source of
its buckling equation nor does it differentiate between spun ends and crown
and segment ends. The Appendix in ref. 16 discusses both of these points.
An equation which is slightly simpler than eqn(6) has also been
suggested by the author:
7. FINAL COMMENT
One of the problems remaining is how the new buckling rules outlined
above should be merged with the curves which are given in Section 3.5 of the
present BS 5500. These curves apply only to ends with t/D > 0.002 and their
origin is a little uncertain. The curves do not seem to have had much
experimental verification but it is believed that they refer to the
axisymmetric large displacement failure mode. It is also believed that the
present curves are safe. Indeed, for cold spun heads made from austenitic
stainless steel, they are probably too safe near t/D = 0.002. It is this facet of
the problem which is being studied currently. It is hoped to arrive at a
satisfactory solution for the transition region around t/D = 0.002 in the
near future.*
This may involve the use of the simple approximate equation for the
axisymmetric plastic collapse of steel torispheres. This is, from ref. 24,
_ 12.6[1 + 2400,,/E](r/D)!-°4
Peo (Di iyk2? (RD yas (14)
Ry, =PIPex
R= ee >
and p=external pressure at which buckling occurs under combined
loading, o, = axial compressive stress at which buckling occurs under
combined loading, p,, = theoretical buckling pressure of a perfect elastic
* For crown and segment heads, the proposed buckling rules merge quite well with
the present BS 5500 rules near t/D = 0.002.
Design Proposals
for Shell Buckling Problems 193
cylindrical shell subjected to external pressure alone, and o,, ..= theoretical
axial buckling stress of a perfect elastic cylindrical shell under axial
compression load alone.
In many cases, there is very little difference between the theoretical
buckling pressures for hydrostatic and lateral pressure and the distinction
will be ignored herein. In this paper, the pR/2t term due to the external
pressure will also be added to o,.
The theoretical interactive buckling curves given by the analysis in ref. 25
reduce to a straight line for reasonably long, simply-supported, perfect
elastic cylinders, 1.e.
R,+R,=1 (16)
For very short perfect elastic cylinders, ref. 27 queries the adequacy of
eqn (16). However, as is shown in ref. 28, the controlling buckling mode for
perfect cylindrical shells made from normal structural steels is the elastic-
plastic one and not the elastic one. The interaction curves for elastic-plastic
buckling are (see ref. 28) curved away from the origin and eqn (16) is a
conservative estimate of them.
When the elastic buckling test results of ref. 26 are plotted in terms of R,
and R,, the plot shown in Fig. 5 ensues. It is quite clear that many of the
points fall below the line R, + R, =1. In the main, this is due to the initial
geometric imperfections in the models which can reduce the theoretical
buckling load considerably, especially for axially compressed cylinders.
Some theoretical results on the effect of axisymmetric imperfections for this
[Ne
Model/LHe R ee
eo} 1 3 |400
[4 2 | 2 [400/526]
e| 3 | 2 |800/3056
[e[ 4 [15«@ li800 [1715|
L
SHS, © | 1 | 800 764
Riy@enl2 viseslzose
[=[9.10] 2 [So0|i906
eg: 4& 4
—1 , 19 e-entnac—|—__|_BR gy
005025040608 10 12,14 16,00
P Ocre
1:0
0:8
re) 06
Tah
0:4
0-2
latter problem are shown in Fig. 6. They were obtained by Koiter (see
ref. 29), assuming the imperfections were axisymmetric and had the shape
shown in Fig. 7. As may be seen from Fig. 6, considerable decreases in
buckling resistance are predicted by this analysis for relatively small
amplitudes of the initial imperfections.
Simple linear, or non-linear, interaction equations have obvious
attractions to designers. Sometimes their form can be derived theoretically
but often they have to be determined empirically. An account of the
philosophy behind the derivation of the latter may be found in ref. 30.
=
wo=A cos q,(R)
q,2182/Rit
Ss. es 6,/Ox,x
02 04 06 08 10“,
x x
Fig. 8. Theoretical elastic buckling interaction curves for perfect and imperfect
cylinders subjected to combined axial load and external radial pressure (R/t = 100,
Z = 200). (a) Simply supported ends, from ref. 31; (b) Clamped ends, from ref. 32.
, Perfect; --—~—, axisymmetric imperfection (A = 0.12).
Rather surprisingly, there have not been many tests conducted on this
interactive buckling problem, particularly in the elastic-plastic region. A
summary of all the tests?°°*’~*! known to the writer, and published
before 1983, is given in Table 1. As may be seen, only a few of the tests have
been conducted on welded steel cylinders and only those tested by Miller
have been on large diameter cylinders.
More tests on this interactive buckling problem are currently under way
in the UK and the USA. However, it will probably be some time before these
results are available. When they are, they will enable any interaction
equations which have been proposed in the meantime to be validated (or
rejected).
If attention is now turned to the test data available on the buckling of
cylinders subjected either to external pressure alone or to axial compression
alone, then the situation is very different. For both of these cases, the
number of tests which have been conducted runs into many hundreds. The
test results for external pressure buckling (from ref. 42) are given in Fig. 9
while those for axial compression (from ref. 43) are shown in Fig. 10. There
is considerable scatter on both these plots, which may be due to initial
imperfections, boundary conditions, residual stresses, non-uniformity of
axial load, etc. However, for design purposes, most national Codes utilise
9qey
T
Burpyong
sysoy,uo [eolpUT[AD
s|fays poyoafqng
07 pauIquioD
[eIxW UOIssoIdWIOD
PUB [BUIOXY SANSSoIg
fay ssoyinp /2pOW dupss [apow ‘ON y/T 1/4 -xoaddy ‘xosddy ainioy dappunog
"9 00x
joraaiou
— Sujuaffiis‘wip fo "q en fo -_ suopuox2 2
(ut) sjapou— ¥) /Aj (uu
Ny) (,wud/
yong Burj rp ‘
97 ‘udIeBUIaA\
TeAW ON WANED car =¢-l —-00r €8'r CL onse[q UseMjog —100'0 —100
uesiojy
pue ayeordnp 008 puess 000 €0'0
9Pl9g dO
SOsEe CUNTUTUIN ON COTO 8~y I =e 69 cs -onselq ‘xoiddy—-v0 -LIl
Aojye (yos) LY onseydSS 61 007
8¢ uesun| Sse[SIXxa]q ON esp'0 8 C798"
=6 00'S — onseg ‘xoiddy
— rs
whe We HO
6€ TONAL “[291§ BUlI—soX
SE9 ¢ =3'O OOS Vie DLE onseyq sury “Oo KEW
popjom sIouayys LOI siouayNs
9¢'0
Op sA0I)pue WITTY -BulI—soX
OTT: = b GaenCGi0) OL 69 OIP onseyq Bury =e OO
udsyLIpIG Aojye ysry) sIOUdNS bro siouayNs
p70
(yy3uans
Ip “‘IOMeM *[291$ ON 0 ce 9 0 sIE S061 c0c ele -onselg “xoiddy -€€0 O1¥
[esagpur poprom
pue 18°0 onseyd dO 880
ITROOW PoAdT]aI-sso3ys
, » porejngyeo
Wo GSQOq sUONepUSUIUIODEI
%) = ‘UONONpel
JO ‘UMOP-YOUY
10]9RJ
JOJ [VIX ‘(WorssaIdUIOD
198 G. D. Galletly
Lower bound of
test results
rer/y,
Lower bound
of
test results
lower-bound curves to the test data and give equations whereby the lower-
bound values may be calculated. These latter are the quantities py, and ox,
which appear in eqn (17) and they will be discussed further in the next
section.
Despite the fact that the various national Code-making bodies must start
with more or less the same test data, their predicted values of py and ox,.are
not the same. The reasons for the differences are discussed in ref. 44 but the
main point is the fact that differences do occur. The magnitudes of the
differences, for axial compression alone and external pressure alone, are
shown in Figs 11 and 12 respectively. As may be seen from Fig. 12 (external
pressure case), there is a large difference between the ASME III N-284
values and those of ECCS/BS 5500 in the elastic buckling region (i.e. the
portion near the origin).
Since the various Codes often predict different values of py and o, ..fora
cylinder with a given geometry, it is clear that the appearance of the same
. PS
~SA ~ SS
0:2F (a) Oyp= 248 MN/m?
[
le
Sak
Sy 2F (b) Oyp= 344MN/m2 Sse,
0 i | al al a a
100 300 600” 100 S00 a econ
R/t R/t
Fig. 11. Predicted minimum buckling stresses according to several codes: axial
compressive load.
200 G. D. Galletly
test results, when plotted in terms of S, and S,, will differ and will depend on
which Code is used. Examples of this point will be given shortly.
The equations used by various Codes to determine py, and ox , are given
in ref. 44 and will not be repeated here.
In essence, eqn (18) has been used by the US aerospace industry for the
elastic buckling of cylindrical shells. It is natural, therefore, to examine it
on the tests listed in Table 1. This will be done in this section and the Codes
used to evaluate py and ox ,will be ASME III Code Case N-284,*° DASt
O134°Dav* andsECGS.**
It is, of course, not necessary to utilise a linear interaction relation in a
buckling design procedure. A part-linear, part-quadratic equation which
includes the effective stress and the yield point of the material has recently
been suggested by Odland.*? The equation is not difficult to use and its
predictions for the tests of Walker er al. were conservative—see ref. 28.
This was also the case when his equation was used for Miller’s tests.
However, in this paper the main interest will be on the linear interaction
relation S,, + S, = 1. To examine its adequacy, all the tests in Table 1 were
plotted in terms of S,, and S, and the latter were calculated according to the
four Codes mentioned earlier. As examples of the results obtained, Figs. 13
and 14 show the plots for the two series of models made from welded
Design Proposals for Shell Buckling Problems 201
Model
S =f1_[07]
500]551 S
Ye [al 2 [o72{soof2u|
ro [#| 3 [os6fsoo]61] 5)
ef 4 [0-18] 500153]
eS
(paz 12
° 0-8
‘4 @ 5
a Spt Sy! Ce
L 1 1 | 4 4 »> 0 wie n ‘nas Fee
TS TER aa 0 04 08 1) 18 20 if
x
(c) ASME Ill N-284 (d) Dnv
Fig. 13. The elastic buckling tests of Miller? on 6.35 m dia. welded steel shells.
Results plotted in terms of S, and S, and according to four different codes.
steel.°?’*! As may be seen for these cases, the straight line S, + S, = 1 was
safe for all four Codes. It was safer for the smaller models of ref. 41 than for
those of ref. 39. Whether this is due to the fact that the smaller models were
stress-relieved or because their ends were clamped is not clear at the
moment. As mentioned before, it is also clear that the same test results may
look safer according to one Code than another (it should be noted that no
safety factors have been introduced into the discussion so far, although
they will need to be considered in a design procedure).
G. D. Galletly
0 0
Oreo eOs 2 tb 7-20 OM O4e Oa 2 16820
(a) ECCS (b) DASt 013 _
0) 0
0) OA OMe Ie Ht AO S 0 04 08 12 16 20.
x
(c) ASME III N-284 (d) Drv %
Fig. 14. The elastic-plastic buckling tests of Walker er a/.*! on 0.32 mdia. welded
and stress-relieved steel shells. Results plotted in tems of S, and S, and according to
four different codes.
Design Proposals for Shell Buckling Problems 203
28 rN
2o[t Stee 6
e@aa
Ole: ae ’ o 4
OS a. & a 8a
e, Ba ee a
1:2 eo a 9 4
io) a8
08 a
-
0:4
0 6
o—»>
DR04 20:8) T2eet he 087 eas OO OL OB VV25 7-682: oO “4 Se
(a) ECCS (b) DASt 013
lo] 1 | 3 [400[3434|
e[ 2} 2 400]1525]
e[ 3) 2 [200|2056,
rel_4 [+5 [800]1715,
Pa ffs,6| 1 |600]764] Sp
Fi) >(nze\|e2e)ges]2035) Opal
20 j=}910]
<2||500]1906] |
a
hor Onn
121° ca aa 4
%9. 4 4
0-8b oq, £
e 4
04 .
0 @®
0
0 04 08 12 16 Sx
(b) DASt 013
reant[276[ra|73]
fa] IV [265] 117 [777 |
ralYe[225]128613 |
rolvs[r86[128[15|
re[wr [265]"9
ra vit2-06[735]25|
re[vit]2-66735]232]
rel[267]200]
tx 1351
WZ
0:8
04 StS yal
> ——s
OQ Oxs Oe WZ. ANS Sy 0. “O08 U8 12" 16 2-0 Sy
(c) ASME III N-284 (d) Dav
Fig. 16. The elastic buckling tests of Mungan*® on 0.45m dia. Plexiglass
cylindrical shells. Results plotted in terms of S, and S, and according to four
different codes.
Design Proposals for Shell Buckling Problems 205
The elastic buckling results of the plastic models (Mylar and Plexiglass)
are also plotted in terms of the S,, S, coordinates in Figs. 15 and 16.
Considering all six series of tests listed in Table 1, the only Code (of the
four listed above) for which all the test points were above the line
S, +S, = 1 was DASt 013.°°°! The recommendations of ASME III Code
Case N-284 were the least conservative of the four Codes. However, this
lack of conservatism will be partly offset by the higher factors of safety
(1.67—2.0) used by the ASME III Code.
The ECCS recommendations give very similar results to BS 5500 for
external pressure alone and it is of interest to know how the ECCS
predictions fare for the combined loading cases. For most of these, the line
S,, + S,,=1 was below the test points. However, as may be seen from
Fig. 16, two of Mungan’s*® tests (under axial compression alone) were
below the line. Since Mungan used a small dynamic disturbance to initiate
buckling in these cases, some researchers believe that his results for axial
compression alone are too low.
It is also instructive to use the 1982 edition of BS 5500 for this interactive
buckling problem. For this purpose, section 3.6 was used with Enquiry
Case 5500/34 (although one is not supposed to do this). The end result
(allowing for the fact that BS 5500 deals with allowable values and not
lower-bound test values) was that all the test points were above the line
S,+S,=1. However, using BS 5500 in this matter gave the most
conservative predictions of the Codes investigated.
Insofar as the choice of a simple linear interactive design procedure is
concerned, one has several options. However, assuming one wishes to
utilise the BSI method for calculating the effects of external pressure, then
the choice reduces to:
If either of the above two options is chosen as a design procedure, one will
have to be careful about the use of safety factors. The ECCS, DASt, etc.
methods are concerned with predicted minimum test values and, therefore,
their predictions have to be divided by 1.5 (the safety factor) to obtain
206 G. D. Galletly
allowable values. The BS 5500 method deals with allowable values directly
and no additional safety factor (SF) is required. Thus one has
S, +S, = 1/SF
Le,
pity Se ae 1 (19)
Px/SF o,,/SF
The quantity p,/SF is the allowable external buckling pressure, which is
determined from Figure 3.6.3 in BS 5500. The quantity o, ,/SF is the
analogous quantity for axial compression—see the curve in Enquiry Case
5500/34.
It is clear from Table 1 that the number of interactive buckling tests which
have been conducted to date on welded steel or aluminium shells of a
reasonable size is very small. Some more results will become available in the
near future when the series being conducted at present, under the auspices
of the American Bureau of Shipping and Conoco, is completed. Other
buckling results, obtained at Liverpool University on smaller models, will
be published later this year.°? This additional information on interactive
buckling will be very useful in evaluating any interaction equations which
may be proposed (the linear eqn(8) or any appropriate non-linear
equations).
Due to the lack of experimental data on welded models of a reasonable
size, any design equations proposed in the interim will have to err on the
conservative side. The simple linear interaction equation, used in
conjunction with the BS 5500 rules, seems to meet this requirement. This is
also the case when the linear equation is used with some other Codes.
REFERENCES
1. Kemper, M. J., ‘Buckling of thin dished ends under internal pressure’, Proc. of
the Symposium on Vessels under Buckling Conditions, I. Mech. E., London,
1972, pp. 23-32.
2. (a) STANLEY, P. and CAMPBELL, T. D., Very thin torispherical pressure vessel
ends under internal pressure: test procedure and typical results, Journal of
Strain Analysis, 16 (1981) 171-86.
Design Proposals for Shell Buckling Problems 207
35). HARDING, J. E., Ring-stiffened cylinders under axial and external pressure
loading, Proc. Institution of Civil Engineers, London, Part 2, 71 (September
1981) 863-78.
36. Reis, A. J., WALKER, A. C. and Virtuoso, F. E., ‘Externally pressurized
cylindrical shells: buckling and collapse under combined loading’, Proc.
International Union of Theoretical and Applied Mechanics (IUTAM)
Symposium on Collapse, University College, London, 1982.
37k Lee, L. H.N., ‘Inelastic buckling of cylindrical shells under axial compression
and external pressure’, Proc. 4th US National Congress on Applied Mechanics,
ASME, 1962, pp. 989-98.
38. MunGan, I., Buckling stress states of cylindrical shells, Proc. American
Society of Civil Engineers, Journal of Structural Division, 100, No. ST11
(November 1974) 2289-306.
39) MILLER, C. D., Summary of buckling tests on fabricated steel cylindrical shells
in USA’, Paper No. 16, Buckling of shells in offshore structures, J. E.
Harding, P. J. Dowling and N. Agelidis (eds.), London, Granada, 1982.
40. Grove, T. and DipRIkSsEN, T., ‘Buckling experiments on 4 large ring-stiffened
cylindrical shells subjected to axial compression and lateral pressure’, DnV
Report No. 77-431, Oslo, Norway, 1977.
41. WALKER, A. C., SEGAL, Y. and MCCALL, S., ‘The buckling of thin-walled ring-
stiffened steel shells’, Proc. Colloquium on Buckling of Shells, University of
Stuttgart, May 1982.
42. KENDRICK, S. B., ‘Collapse of stiffened cylinders under external pressure’, Proc.
Conference on Vessels Under Buckling Conditions, I. Mech. E., London,
December 1972, Paper C190/72.
43. BrusH, D. O. and ALMROTH, B. O., Buckling of bars, plates and shells, New
York, McGraw-Hill, 1975.
44. GALLETLY, G. D., ‘Buckling of fabricated cylinders subjected to compressive
axial loads and/or external pressure—a comparison of several codes’,
University of Liverpool, Department of Mechanical Engineering, Report
A/059/82, March 1982. Also (minus the Appendix) in PVP-57 Pressure vessel
design, G. E. O. Widera (ed.), Bk. No. H00214, June 1982, New York, ASME.
45. ASME (American Society of Mechanical Engineers) Boiler and Pressure Vessel
Code Case N-284, ‘Metal containment shell buckling design methods’, Section
III, Div. 1, Class MC. Supplement No. 2 to Nuclear Code Case Book, 1980.
46. DASt (Deutscher Ausschuss fiir Stahlbau) Richtlinie 013, ‘Beulsicher-
heitsnachweise fiir Schalen’, Cologne, Germany, July 1980.
47. DnvV (Det norske Veritas), ‘Rules for the design, construction and inspection of
offshore structures’, Appendix C, Steel Structures, Oslo, Norway, 1977.
Reprinted 1979.
48. ECCS (European Convention for Constructional Steelwork), “European
recommendations for steel construction, Section 4.6, Buckling of Shells’,
London, The Construction Press, 1981.
49. ODLAND, J., Buckling resistance of unstiffened and stiffened circular cylindrical
shell structures, Norwegian Maritime Research, 6, No. 3 (1978) 2-22.
50. GALLETLY, G. D., ‘Towards a design procedure for the buckling of cylindrical
shells subjected to external pressure and axial compression’, presented to
British Standards Institution, Sub-Committee PVE-/5 (Design Methods),
October 1982 (Document 82/79229).
210 G. D. Galletly
J. G. A. CROLL
University College, London, UK
ABSTRACT
NOTATION
E Modulus of elasticity.
ve Function of Z used in eqn (19).
g Function of Z used in eqn (22).
a| Non-dimensional axial half-wave number (= (j/j,,)).
i Axial half-wave number.
k Coefficient used in eqn (19).
Kwe,kwe, Normalised equivalent load imperfections.
l Cylinder length or ring spacing.
M Moments for fundamental equilibrium.
m Incremental moments due to buckling displacements.
N Stress resultants for fundamental equilibrium.
n Incremental stress resultants due to buckling displacement.
Dp External pressure loading.
r Radius of cylinder.
t Thickness of cylinder wall.
W,, Wy, Total equivalent bending and membrane imperfections.
w Radial displacement—positive inwards.
x Axial coordinate.
Z Cylinder (Batdorf) parameter (=/*(1 — p?)!/?/rt).
Superscripts
Ie Belonging to fundamental equilibrium.
e Relates to ‘equivalent loading imperfection’.
o Relates to initial geometric imperfection.
ae) Location for maximum stress or displacement.
Subscripts
cl Indicates classical critical stress or mode.
C; Critical stress associated with jth mode.
FY First surface yield state.
FP First full section plasticity state.
imp Indicates geometrically imperfect cylinder.
J Relates to jth mode.
Axisymmetric Collapse of Cylinders 213
1. INTRODUCTION
Tests show that cylinders under various combinations of axial and pressure
loading may collapse into one of two distinct deformation modes. For
short, thick cylinders the high bending stresses induced near boundary
constraints combine with the membrane stresses to produce a collapse
characterised by axisymmetric deformation modes. This form of collapse
has been the focus for much recent theoretical analysis.'~ 7 In contrast,
long thin shells are observed to collapse into modes having circum-
ferentially period, asymmetric, modes. Elastic-plastic collapse into both
these modes, shown in Fig. 1, displays an at times considerable degree of
sensitivity to the levels of their respective initial geometric imperfections.
Between these two limiting classes of stocky and slender cylinder there
exists a class of shell for which collapse could take the form of either an
axisymmetric or an asymmetric mode. Which will control buckling collapse
depends critically upon the precise shapes and relative amplitudes of the
axisymmetric and asymmetric initial geometric imperfections.
(a) (b)
Fig. 1. (a) Axisymmetric and (b) asymmetric elastic-plastic collapse modes for
cylinders.
214 JSGZASGroll
Existing design practice does not distinguish between these two quite
distinct forms of collapse. Consequently it fails to emphasise the importance
of controlling both axisymmetric and asymmetric imperfection forms. If
future design of the large cylindrical shell components being incorporated
into many forms of marine structure is to be put on a more rational footing,
it will be necessary to provide separate guidance on strengths against each
of the axisymmetric and asymmetric buckling collapse modes, and to relate
these to clearly delineated tolerance limits on imperfections in each of the
respective buckling modes. To do so will require analytical procedures for
estimating collapse that are reliable, but also sufficiently explicit and simple
that they can be effectively incorporated into design practice.
For periodic buckling such design guidance could conveniently be based
upon the recently described reduced stiffness lower-bound procedures.*~ !°
Analysis of axisymmetric collapse has been concentrated on sophisticated
finite element and finite difference computer codes which, for the numbers
of apparently independent parameters involved in cylinder collapse, are
not immediately suited to the presentation of compact and accessible
design guidance procedures. The following summarises an analytic
alternative which could provide a basis for the presentation of design
guidance against axisymmetric collapse of cylinders. At the very least it
could provide a parametric delineation which would permit numerically
based design guidance to be presented in more compact and universal
forms.
¢,=-9¢
(oy = —He (1)
which, because it involves no radial expansion, is independent of the nature
of end support conditions.
Axisymmetric Collapse of Cylinders 215
w—w® o wr+ws\.
( , ls » eae )sin (3)
izcilig Onna.
where the critical axial stress o,, associated with the jth mode is
a, [(t\? Cin)? i? 1
E -|(;) ay + (7) od oe
This spectrum ofcritical stresses reaches a minimum at the classical critical
stress
~ 0.592Z1/ (6)
i ec
where the Batdorf parameter Z =/?(1 — yv?)'/?/rt. It will assist subsequent
analysis if the critical stresses of eqn (4) are reformulated as
Cena
a3 (s+5) (7
cl
conditions. For the simply supported cylinder these pressure induced radial
expansions give rise to ‘equivalent load imperfection’ of the form"!
ws 2(B — p)
t (yJxBd—p*))"
The total response of eqn (3) is made up from the amplified model
responses as depicted in Fig. 2. Associated with this non-linear buckling
behaviour the elastic membrane and bending stress resultants would be
N,=NF +n,=-—ot
ss et ee be e
Ng= No +g = — pot — 7 (w; — w;) sinjrx/l
Siewaee
Et® jt 2 ei cctael
MS er eae » (=) (w;- w;) sin jax/l (9)
jaoney
Me Ma hee OS os(amy
pet »
jn\* oe:
(=) (w; — wj) sin jrx/l
psoas
This model of non-linear elastic behaviour has been shown to reproduce
that obtained from a full non-linear finite difference solution.’
Axisymmetric Collapse of Cylinders 217
(10)
No —A Ast
rang Wy mt
+ Wh,UW)
where all stresses are normalised with respect to the uniaxial yield stress,
o,, so that A=a/o,, A..=o4/o,, and the total equivalent bending and
membrane imperfections are defined as
W,=
“
» (Ag A) (weJ ;+e)
3 ah (v=
- A) 4 sin j7x/l
(11)
Wu = » EXC Vee eine (“2““1)sin jrx/l
Gay pe A)
ag = — Ad
a, = +2A,
€ Va
a
Ox;
(a)
ie SS
y ee
Va ~
y / = =
Vie
£ me = /t-tn
fo) e fo) /,
Xe a 8i in 7
/ 7
Ox;
(b)
Fig. 3. Stress states used for (a) first surface yield, and (b) first full section
plasticity.
Axisymmetric Collapse of Cylinders 219
5 ZESD8A9 AG a)
[
2= T2849 Ue)
(14)
21
and the critical surface depends upon the relative amplitudes and senses of
initial imperfection and external pressure loading.
Again, solution effectively depends upon just the three parameters A,,, W,
and Wy.
It may be noticed that the equivalent imperfections (W,, W,,) depend
upon A,p. This means that the coefficients a, and b, in eqns (13) and (15)
respectively are themselves implicit functions of A,p. Solution is easily
undertaken as an iterative process, with just a few cycles required to
determine A, and A,;p. The whole process can be easily adapted to
microcomputer applications.
To illustrate the nature of the solutions from eqns (13) and (15), Fig. 4
summarises the behaviour for a cylinder having A,, =¢,./a, =4.0 in the
form of a two-dimensional imperfection sensitivity plot. In all cases shown
the surface or section corresponding with the lowest solution for eqns (13)
and (15) is given at the location indicated in eqn (14). Accompanied by
suitable charts or figures for specifying (W,, W,,), imperfection sensitivity
plots like those shown in Fig. 4 could be used for design purposes. Before
this can be achieved consideration needs to be given to the appropriate
forms of initial geometric imperfection.
Making use of the general forms given for the total bending and membrane
imperfection parameters (W,, W,,), eqn (11), it is possible to consider any
arbitrary form of initial geometric imperfection. Imperfections in the form
of a single axial harmonic j have often been considered. Of all the
wavelengths for imperfections it is that associated with the minimum
critical mode that produces the most adverse reductions in collapse
load.** This form of imperfection is also particularly convenient in the
present analytic formulation,’* and it would through appropriate choice of
imperfection amplitudes be possible to build-up a code specification based
upon the ‘worst shape’.
But of considerably greater direct practical relevance is an imperfection
induced by typical fabrication processes. As illustrated in Fig. 5 these
would arise from any lack-of-fit or shrinkage associated with the welding at
——S SSS,
CEE EN Ee
Ig =2l ring = lack
icone stiffener of fit
ring a
weld P,
l l
Sores | = = —) oo 2
| | ee SESH in
edge fg Wimax edge" W max
shear shear “
|
(a) (b)
Fig. 5. Localised distortions of shell wall caused by (a) weld shrinkage and
(b) the lack of fit of shell plating and ring stiffener, and modelled using fabrication
imperfection given by eqns (16)-(20).
222 J. G. A. Croll
Notice that for Z > 5.698, and on account of eqn (6), the gauge length
1, = 21 =~ 3.46(rt)'/? over which axial imperfections are likely to have their
maximum deviation is close to the 4(rt)/* recommended in many current
codes.'*~'* Choosing the maximum amplitude over the gauge length
Lael to be the prescribed geometric tolerance w,/t, allows c, to be
calibrated as
oO
WwW
t
= (18)
sin jri/l
w.. i 3 ~
iL) oT 2 (20)
ee Z
while the value of k would depend upon the precise nature of the fabrication
process. Design codes usually employ this form of tolerance specification,
Axisymmetric Collapse of Cylinders 223
1:0
f,g 1
f(Z)=——— (Gn, (19))
e Fewoa
Sis ereeg)
7 1
06 giz) = ; cats (eqn (22))
Vi Neljseok ed
0-4
0:2
10 20 30 40
Zz
Fig. 6. Geometric parameters used in eqns (19) and (22) to define fabrication
imperfections.
with DnV adopting k =1/30'* and DASt Ri 13’? and ECCS'* using
k = 1/25. Weighting function f(Z) is taken to be unity in DnV,'? while
ECCS!* recommend f=0.25Z'/? and DASt!?
f= 500Z"/?/1.
Based upon the fabrication imperfection of eqns (16)—(20), and k = 1/30,
Fig. 7 shows typical interaction collapse loads based upon the solution of
eqns (13) and (15).
———
f +/3a woy uba (LZ) te
J/) L-0= wee
XV
n=g
0-=Xpo
3/40
= 5100-0 K
4/0= O€/L (4/4)
3, (2) Z-4 = jy-= og-
(e) J
wor
J. G. A. Croll
Ott Once
oi
kd_ ko
‘SIq°z uosueduiog
jo joexo pur ajewixoidde uoneoliqey ‘suonooysodurt
Burmoys
isry jjny Ayonseyd
se e ainsvour
Jo
OLJOWIUWIASIXe
asdeyjoo
Joy jerxe pue ainssaidpopeoy Ajduuts poysoddns ‘ssapurAo
stays Sutaey(e) 4/7= [0 pue (Q) 4/)= ‘S70
Axisymmetric Collapse of Cylinders 225
‘pluor—L
*31q
diy
(Z)d
(94)
(12) ube wouy 47),
0E/L= (4/4)
ey
uba
wosy
SL00-0= 3/40
$20 =4/1
4/3m
4/0™
[
226 J. G. A. Croll
ae w°
Cc, = sinjni/[]=—— t 9(7) (22)
Er
Hwee gece J
the imperfection profiles based upon eqn (21) differ very little from those
given by eqn (16) for the range of stocky cylinders (those for which Z < 50)
for which the present axisymmetric mode is the most likely form of
collapse. Not surprisingly, little difference can be observed in the first yield,
Ary, or first full plasticity, App, loads based upon these two alternative
fabrication imperfections for the range of shells shown in Fig. 7. If this
form of approximate imperfection was to be used for more slender shells it
would be necessary to recalibrate the tolerance specifications to ensure that
the results are representative of those obtained from the more realistic
fabrication imperfections of eqn (16).
The real advantage of the approximate fabrication imperfection of
eqn (21) over either that of eqn (16) or any of the alternative imperfection
forms is that it may be directly combined with the equivalent load
imperfections. This means that the total equivalent bending and membrane
imperfections of eqn (17) may be rewritten
W, =(B—B)kwe
(23)
Nia B)kwe,
ce JS
> t ea i
la owase
43a4uad}J8ys XV x5
vt i
JIajsedwi
}9yS
asde}]}03
sn30)
ajewixoudde
ubisap ‘sn30)ube (€€) iy i
J/) S2-0=
4/4= O0L
3/Ao
= S$L00-0
/oM
4= (4/J)OE/L
2 (Z)3
Axisymmetric Collapse of Cylinders
31°§ Suiddeyy
jo yoajiod UONORIOIUTWIRISeIp
0} Als JOoJJOdur! WOHOeIN}U! ‘sueISeIp
BuIsnoy) uonewoj
jo ubasuen
(LZ)
227
228 J. G. A. Croll
eles (26)
where f’ = c,x[3(1 — n”)]'/?/2 represents for the geometrically imperfect
system the pressure loading required for the geometric imperfections to be
exactly nullified by the equivalent load imperfections—in much the same
way that the load eccentricity may nullify the out-of-straightness in an
axially loaded column. The constant c, required for the specification of f’ is
obtained from eqn (22) or directly from Fig. 6.
As a result of the above mapping the load interaction diagrams for an
imperfect cylinder (f’ 4 0) may be obtained directly from that of the perfect
cylinder (f’ = 0). All that this entails, as Fig. 8 illustrates, is shearing of the
perfect interaction diagram such that at the appropriate value, f;,,, for the
loaded but imperfect cylinder the collapse load would be the same as that
for the perfect cylinder with £,..,, where
p
aeSo aod oakVicia |
is Ws ee =)
ie ‘lr W=Wr
Pa a4 ie —_ ate
f, :
half ring force te
(qs 2/2 kp Wr) at} = aeons
shear
of
ring — Qs = Ks We = 1/2 kp We
Ws =TWs Ws = Ws
Fig. 9. Interaction between subshell and ring, allowing shell displacements for
flexible rings to be expressed in terms of those for infinitely stiff rings.
cylinder relative to a flexible ring, w¢, is related to the shell displacement for
a radially rigid ring, w,, through
w= nw, (28)
where yn = k,/(k, + 2k,). The ring and shell stiffness coefficients (k,, k,) may
be easily computed in terms of the ring parameter (A,//t) and the Batdorf
parameter Z, where A, is the cross-sectional area of the ring. Having done
so the total equivalent imperfections now become
W, = (n(B— 2) — B’lkws
(29)
Wy = [n(B — 1) — B’Kws,
and the interaction diagram for the imperfect ring supported cylinder once
again may be obtained from a straightforward transformation of the
230 J. G. A. Croll
perfect simply supported case. In the present context the collapse load of
the perfect system at .,,.,, will relate to the imperfect flexible ring, loaded in
such a way that the actual load ratio is f,,,,,ring? when
Otherwise the treatment of the flexible ring follows exactly the same pattern
as the approximate fabrication imperfection of the previous section.
For subshells that occur near end diaphragms a closely related technique
to the above could be developed. However, for these end bays, design based
upon rings having infinite radial stiffness should be conservative.
6. A DESIGN APPROACH
“boiler code:
1-0 a a a es el Eicien ee
Pp/py
/0g
Oy
=
10 20 30 40 50 60
Zi
Fig. 10. Radial pressure loading (f — 00) to induce first surface yield and first full
plasticity. Used to locate points B and C in Fig. 8.
limiting design state. Notice that for internal pressure the equivalent hoop
stresses will be equal but opposite to those under external pressure. With
the help of Fig. 10 the points B and C may be located in Fig. 8.
For an arbitrary inward, w;, and outward, w?, imperfection range,
appropriate values of f’ to locate the points A; and A, may be calculated
using the function g (shown in Fig. 6) to give
Bo=EU+¢3 -= g(Z)
(32)
R= U>cs— 8(Z)
where c, =2[3(1 —w)]'/7/2 ~2.592. Having determined the failure
pressure, Py = — Pc, from Fig. 10, design could then be based upon the
linear interaction formulas
A+c,P =c; (33)
7upyo J. G. A. Croll
(= Fe Pig)
In the linear interactions of eqn (33) the values of 6’ appropriate for inward
and outward initial imperfections would be chosen from eqn (32) and used
to produce the two sets of interaction curves shown in Fig. 8. For design it
would be the lesser of these two solutions that should apply.
7. CONCLUSIONS
REFERENCES
one a
oun &
13
G. DUTHIE
British Gas, Midlands Research Station, Solihull, UK
and
A. §. TOOTH
University of Strathclyde, Glasgow, UK
ABSTRACT
NOTATION
Radius of vessel.
Distance from end of vessel to centre of loading patch.
Half-width of the loaded area in x direction.
Extensional rigidity = Et/(1 — v?).
lave) Differentials with respect to x/a and ¢ respectively.
Modulus of elasticity.
Flexural rigidity = Er°/12(1 — v’).
TrmRHOTSTS
Non-dimensional factor = ¢?7/12a’.
235
236 G. Duthie and A. S. Tooth
1. INTRODUCTION
Fourier expansion methods have been available for many years for the
derivation of displacements and stresses which occur in shells and pressure
vessels. They have been particularly useful for solving local loading
problems and a variety of solutions exist for vessels subjected to point, line
and distributed loads. Those relevant to cylindrical vessels are presented in
refs. 1-12. In general the solutions for the cylindrical vessels are based on
simplified shell equations of the form proposed by Donnell!? and,
because of the limited computer facilities available when the original work
was undertaken, have often only been evaluated to a limited degree of
accuracy. Nevertheless, a variety of results generated by these solutions
have been incorporated into National Pressure Vessel codes, for example in
Local Loads on Cylindrical Vessels 237),
BS 5500,'* and when employed along with the appropriate design code
stress limitations have proved extremely useful and reliable.
{t has, therefore, become normal practice for designers to employ the
most appropriate results presented in the codes and the exception rather
than the rule for designers to attempt to develop rigorous analyses for
problems under consideration. This is not surprising when one considers
the complexity of Fourier expansion methods of solution, the timescales
and computer resources required, with the success of current code design
methods.
However, advances in micro-, mini- and mainframe computer facilities
have made direct application of Fourier expansion local loading computer
models attractive. Indeed such models need no longer be considered merely
as a research approach suitable for use by only a few. Although these
Fourier expansion solutions involve a significant amount of computational
effort they do not consume a great deal of computer core store and could
therefore be mounted and run on relatively small computer facilities.
The authors have used the method extensively in the derivation of the
stresses and displacements which arise when a cylindrical vessel is
supported on saddles.!°~ ‘8 In the process of developing the above analysis
a wide range of particular local load cases were examined to test the validity
of the approach. These cases are considered to have a general interest in
design and for that reason they are presented in this paper. Attention is
confined to the results which arise from one particular series solution for
circular cylindrical vessels subjected to mechanical loads in the radial,
tangential and longitudinal directions. These series are particularly useful
when the loads are applied remote from the vessel ends in such a manner
that the flexibility of the vessel ends plays no significant role in the
behaviour under consideration. A similar approach may be used for other
vessel end conditions using different series. Thermal stress problems may
also be tackled using a similar technique.’?
2. THEORETICAL ANALYSIS
(1 —v) 3 = YW
frat OS) aor — hah +{o- k ——
5 aro}
The positive directions of the mid-surface displacements uw, v and w and the
loading components Px Po and p, indirections x, and zare given in Fig. 1.
An intermediate step in the derivation of the above equations generates
the following expressions for the stress resultants in terms of mid-surface
displacements:
ad gat eeegroa
coe =e aE vie a v) ;
K
M, ma [w+ 6?w + vd?w]
Height of fill
K
OS) (1 — v)[dow
— dv]
OssK [(1
— v)(d* z ow
— d*v)
2
+ dw + 6°w
3
+ vd’ Disksdw]
(2)
Because of approximations made during their derivation, eqns (1) are
linear and can therefore be manipulated and simplified to give the following
partial differential equations relating mid-surface displacements to loading
components:
4
V
{DEXP}w =~ p,a?/D
(1 +) (1 —v)
2 doe REL Wid Oe Ta do
143
can = 2 as?map
{DEXP}u = eel] a
7 [ (i) d® + k(2 = v)d*d? +k(1 + 5?)?6?
1 ==
Gye epmy
ep omme oem
: ”)a4254 a | p,a?/D
1
(DEXP}w = —> [2kd*6 — (2 + v)d?6 + 2kd?6° — 5°] p4a7/D
ie 2
cea) a + 609? + 3¢425°|
v8 a (d? ae Bly
VS Fe)?
Local Loads on Cylindrical Vessels 241
The choice of the above expansion implies that the loading system is
symmetric with respect to the generator passing through ¢ = 0. They also
imply that certain boundary conditions, discussed in detail in 3.3, must
exist at the ends of the vessel. Other expansions can be used to represent
loading conditions which are not of the symmetric form described above.
Substituting these expressions into eqns (2) and (3) leads to the following
solution for the mid-surface displacements and stress resultants:
5 ie0) io6)
zy
u= > » » Z 21m PnmCOS NG COS (AX/a)
n=0 m=0
242 G. Duthie and A. S. Tooth
2 = =
a :
N, =F Z5
1mPn,m COS Ng Sin (Ax/a)
n=0 m=1
ie ae |
a my » » Z6,, mPn,m Sin ng cos (Ax/a)
n=1 m=0
— =
rm)syZ6A 4 mPn.m Sinn cos (Ax/a)
n=1 m=0
M,2=a@7(1—) »)»:Z9
1mPnm Sinnd cos(Ax/a)
n=1 m=0
= «@
ZA, = (4 4 97)7/DEN
==:
|— 1+3
egg anyiicnes = ale ; y in|DEN
244 G. Duthie and A. S. Tooth
+
Cis f= 64° Fa one
k
For each loading component considered the Z4,,,, to Z12,,,, coefficients
can be calculated from the corresponding Z1,,,, Z2,,, and Z3,
coefficients using the relationships
70 an Ze
ZAl nm = —N( =n? —2?)Z 1am + (1 —v)A7Z39on
N Ss ] =A? +n) Zl E =;a | Zo. ni a ;” Zam
Referring back to eqn (5) it is noted that the only unknowns are the loading
terms p,,,,.The method of deriving these is considered at a later stage of this
path.
(b) No rigid body displacement or rigid body rotation of the ends can
occur.
(c) The ends cannot carry applied axial loading.
(d) Generators are free to rotate in a plane normal to the profile.
If the vessel ends conform to the above boundary conditions then the
Fourier expansion of eqns (5) are a complete solution to the problem. No
complementary solution need be added to the particular solution since
eqns (5) satisfy both the governing differential eqns (1) and the boundary
conditions of the problem. However, in practice some deviation from these
conditions is likely to occur; for example, in the case of a vessel with a
torispherical end closure condition (a) is not ideally satisfied. In such cases
it is still possible to use the results of eqns (5) with confidence if the local
loading is applied some distance from the vessel ends.
p= » Po.msin mnx/I
m=1
Multiplying both sides of this equation by sin (m’1x/1) dx d@ and, since the
loading is symmetric about the vertical diameter, integrating over half of
246 G. Duthie and A. S. Tooth
1 / /
MTX
{| asin dxd@ = intYmain ee stl dxd@
0 10)
Noting that
4 dx =U when m 4m’
l
== 1/2 when m =m’
we obtain
2 1
Repeating this procedure for all values of n, m and all loading components
the following set of expressions for p,,,, are obtained.
sf
For tangential loading
4 |
= Po as sinn@dxd®@
0 (n>0Q)
Pan
’ =F1 |
a \p,.Aaxd@ C0)
2 1
ran =z |IE
p.cos —— dx dg (n=0,m>0)
In Jo
9) i
Pan | ie
p,cosngdxd@d (n>0,m=0)
BG Mis
4 1
Ian=F | |p.cos
: lie ite
~~ cosnp dxd¢ (n>0,m>0) (6)
Local Loads on Cylindrical Vessels 247
N, 12M,z Ng 12My6Z
Pas Byaane 7s tab t ie
ete. (see ref. 20):
The above approach is of particular value in dealing with loads
distributed over discrete areas (local loads) and can deal with both
uniformly and non-uniformly distributed loading. The solution is not
intended to apply to point and line loading problems which can be dealt
with more efficiently using other methods.”'8~ '! An acceptable solution for
a point or line loading problem can, however, be obtained by considering
the loading to act over a small but finite area which is, in fact, a reasonable
approximation to what occurs in practice.
lowest generator (f = 0°) and at a distance b from one end of the vessel as
shown in Fig. 2.
This loading system is defined by the expression
P,=p inthe region —8,<@<f, and (b-—c)<x<(b+c)
P,=9 otherwise
Substituting this expression into eqns (6) gives
2 le Se manx
Poni | |Pp, Sim ] dxd@
0 JO
4 1 us
Pram = 7 ||p, sin = cosn@dxd@d
Tt (0)
4 b+c (B,
= — psin oe cosnddxd¢
In l
Spe tres IRD +, (INR ©
= > sin sin i sinnp,
4 ] Bik : es b .
P, = 5 ae E + »ysin np, cos na |Sita ogee sie (8)
m= n=1
ras
Hele
ve
:
ne
EDs wue-o5=9“wUGLjZI=
L i | “5='4 wusoe=4 ;S
9 a=
= Ss is jo
vaquinn Ul
SWJa4 2
E4 S
Salas W=U zy
4 al; G2 eo bey “
ee
ow 3
”
( ,ww/N ) d aunssaud
El —s
|
JeIpey
[ joexg a
ae
0
L-
q <+—X
=) otherwise
It is noted in the above that because of rotational symmetry only one term
in the n series (n = 0) exists and in eqns (5) the terms v, N,, Ny,, Nig, Mgx;
M,,, and Q, all vanish. The remaining terms can be obtained using eqns (5).
The best known classical solution for this problem is that of the infinite
cylinder carrying a uniform line load round a cylindrical profile.?°?* A
small steel cylinder, length /= 508 mm, mean diameter 2a = 229 mm and
thickness t = 1.6mm, was used for the comparisons. In the case of the
Fourier expansion solution a band-width of 0.8 mm was selected. In both
solutions the loading was taken as 1 N/mm of circumferential length.
Table 1 shows the comparison for w, M,, M, and N, with increasing
distance x from the band of loading. The only significant variation occurs
in the M, and M, values at the x = 0 position, no doubt due to the fact that
the classical analysis assumes a line load of zero width. In point of fact if the
modification to include the band-width given in ref. 22 is used the
agreement is within 0.1%. A series of different band-widths (3.2mm,
12.7 mm, 50.8 mm) were examined using both the Fourier expansion and
254 G. Duthie and A. S. Tooth
the infinite cylinder case with the band-width modification. The agreement
is always better than 0.1%.
It was found that, with the exception of the zones near the cylinder ends (i.e.
Local Loads on Cylindrical Vessels 255
4pa_
Pam=—— — acos a)
(sina m=1, 3,5...)
(n=0,
mn
nm Se
mn? a in2 a)
sin ; (n=1,m=1,
A=1,mM=1,2,3,5...)
ehec!
8pa ; :
Pn (sin na COS & — NCOS na sin a) (ies ede
~ mren(n? — 1)
TaN)
diameter,
Za
(MN/m*)
Axial
N,/t
stress
membrane
of
6M,IP
(MN/
stress
bendi
Circumferential
¢ (degrees)
Fig. 5. Comparison with barrel vault solution.
Local Loads on Cylindrical Vessels 257
Load
—o— Bijlaard (ref 3)
A Duthie and Tooth
eS euhy
Mg/P
/
Circumferential
Radial
Moment
6. CONCLUDING COMMENTS
When dealing with a vessel subjected to local loading over small regions
of its surface it is common to assume that each loading component is
uniformly distributed, or in the case of an applied moment is linearly
distributed, over the area considered. This assumption implies that the
loading system, for example a bracket, is flexible enough to take up the
displaced shape of the vessel. Obviously this is not the case when the vessel
is loaded through a very rigid attachment or the loaded area is large like a
saddle type support. However, for cases in which the loaded area is small,
the maximum stresses predicted for the assumed flexible loading patch are
found to agree reasonably well with the stresses occurring in practice
although the maximum values occur at different locations.
Recognising this fact it is possible to use the compendium of p,,,,loading
terms, presented in Appendix A, to analyse the vast majority of loading
cases likely to occur in practice and obtain stress values which are of
genuine significance for design. The ease whereby these can be obtained
258 G. Duthie and A. S. Tooth
Acknowledgements
The authors wish to record their thanks to the United Kingdom Science
Research Council and to the British Gas Corporation for the support given
to this work. This paper is published by permission of British Gas.
REFERENCES
APPENDIX A
ariitt
TE
LE SS SN
Eas
I
: dad
(esos)
1>x30
aa)
“T= gE
ST Cease “d
UZ>PS0
P+)
SxS(9-
O=u
=u
Oe
“eT
|
Local Loads on Cylindrical Vessels
ps0)
zs
ee ete pau) (
u~Soso oy ef
ee” On
261
(panuijuo2)
aiqeL ‘piuos—e
N
\o
ON wossvig
fo SUIpvOTwWasds | PG a
suUlpvo] wajsds =) 00 (asim1ay10 =) (0 40f‘u uw sanjpa
jou (pajonb
I C="
'g- os dS eae eee
7 a= boty)
fom: as
Il eG
6 = =u6 6 © eee
(gustaen)i (
| ee wy o
4
oans(I \SOSCJ—'P)
ee z hire
o=u ee 3
aM
S F 1>x>0 ty be 66 tee =
eels aS)
=
+'o-)So5('9
Cg
eee
—'0-)
-
(‘pu)
[7
urs
gu)
(
= x
soo
Ue uu , ; a Cre
Be
—
S
=S
td J Ei) (
ee sear
a
9 Cc Hy Secuve
'g-— isos
yee (9+q)5xS(9-q)
ig
— usp (Tf) eee, )
ee
'g@So5'¢-
«2 (dea
1)) ia
(9+9)SxS(9-@) Can
Jena
'gsos'g v =u ‘TT‘€
— uIspy (‘gu)
(
a)ae
a
G7 =e Cag= )
”
‘p) (eine Q=uU
(9-9)
<5>
ec da “) re("
(ing)
—"O- \SoS('g (Gao
4
Sxs(0—¢9
OF)) ~ soop.7(pu)urs (‘u)‘)i
u
=
seyee)
lat
Z e
9) S(9—=x V+9)
ayiaIe i
Ol
are (eee
)a- 9)Sx(9+
5(9-q)
SoS (19+)
ja—19) ('g
(
8
soapy“
(9+9—))>x5(0—q=) (“ou)us
) (gu) =u os (
Local Loads on Cylindrical Vessels
—'o-)
('g SoS '-) ("9+ Sa SeeEe BOS
(Ogss ) ei (2a ie
)a-
263
264
“p IIqeL “suIa}sAS BuUIPeO’T [eIpeY ofoqeieg JO} SWI] SUIPVOT
‘ON
wajsds Suppo] 0'0=)
40f‘u uw sanjpa
jou
(pajonb
jig Q=u
_ G+ q)sx>FC-q) € ‘7, =u
F Lar CH ee
en)
uz! eral
( Z Jiga (
(9+9)>xS(9—-@q)
heta Sian
7 asesy-) ( : )ea Cees
(9+¢-)S>xS(9-qG-))
tet Ce
G. Duthie and A. S. Tooth
Teau\., eae
Ter
a (ome
—'f-) +'9
(*G¢-)So5('f
een ee(
N19 e
Os ere
f= ( eae (
as(9-4) SxS9) (9+
eae ras
ores
9-*)V=17
+ 8d ((9-D—-*)¥=7T
+ 899d -b9=€T
d+ EI Sra recsurs X('gu)
uu uu
195 \
uls= pe gua) a sum) UISs=79qua}
soa!
\ Ee ay =F3 ae WAOUS Wel= 91aoa CAISE
d JZ é
i ly 'fd+ (‘gu)us7g=7x
+ soo gu) —('
= us (Igu)5+ uis (‘gu)
gq
Local Loads on Cylindrical Vessels
‘ON
ead wajsds SUuIpvOT fo wos3vIG
“SUIdISAS
SuIpeoT
SUIPCOT
[PIPPY
IVIUIT
SUID,
sqey
JOJ
°s
266
—"))(19 565!o) (f+
lI
4)\.,
“d
Uv(a
i sere Cee PO peu) (
eee
tH Cr Oe
ae eee Ce
( e else
(9+q-))SxS(9-q—-eo
)) gargs Se
p=) (ig o> l= = Cg
mae ee ee
(
(‘ou)sooz7dxZ7=¢d
=rqd = 4a'g =¢q er urs —_AA("gu) (‘pu)soogg*z7=9q
UU : UU
Local Loads on Cylindrical Vessels
— uu 97)uu
— (7 JUU fay amin! pny) —97Z)uu
(1
Lf= uls aS uIS | °y soo pea! cia uIS f +
gum}
uls zit soo a | ny edie
eee
ee ee a n
267
268
sIquyL
‘9 BuipeoyTsway
JOJ [eIpeY BUIPLOT "SUI9]SAS
wopisvig
fo
BUIPDOT
‘ON
Wajsds CERO
ee
peau SSNS
|H="d (~SOS0-)
ea
PODS
ce hee
G. Duthie and A. S. Tooth
ee
=—d
Geers
len
PD LF LP FL a a Ze NN RCE
®y,
u~Soso
A ew="a
)
say iieeiaes
PPLPLE LE A OE CTP
u~Soso
Q—¢g—S45
)) 1) Ot¢
|
SR,
ee (gay)
“d =
1>x>S0
LH = yd = $09)pd
— soo(x 1)°d=7y
+ (soo Psool'd= Ddp
1D EeH
Local Loads on Cylindrical Vessels
US)
» — S09”
(0
zulu
® voz )=(
07]
— urs [(0Z) ana?
uls] dp° quu} IUUL
= zulu
(nu) ns00
— soou (xu)uls[w = pO —_=uls} ——
,¥)u,um
— (1
(
wu ]
1
=¢D (um)/1*dp
269
ORL
“L BuIpeoTSWAY
IOJ [eNUAsUL]ZUIPROT "Sua}SAG
270
‘ON woisvig
fo SUIPDOT wajsas GL
BUuIpDO] Wajsds =) 00 (asimsayjo =) 0.0 40f‘u ul sanjoa
jou (pajonb
DAAAAARAARALS
= u>o>0
Ad RA RRA eta —e\()
I Tay tote
A? ARAARBAAL lg
tee yeneetty >P>u
UT cers
ae
(
S
4
u>o>Q =
: Sx50-9)
Q oo :
oe x-=.
ieee 4
S
u>o>gQ
§
F ee)
ae S‘eT=u
Cee Soke
ie
XS Xe
ee (ee
lees (cad
ee) ae
u~>p>uu>o>g u~>p>u
—'6) So5('
+o) (9g
Ae
io) xs Cheam
(a) UTS
"4 gs0s'y_)@use=Me ae
ea se
eet Sak
th= Ce
sa ie (cae
£ 7Y= [7 uls (aa
us (se)¢g = zg us us(‘ou)
(‘gu) ust—'d
(dz)
]w=pq
Local Loads on Cylindrical Vessels
89 = ag)
7) ['g(u— Djus igenchar ener Ga
3 eal
(u—1)Z [)Z (u+ zu / /
271
°*g aIQUL “suId}skS SUIPeOT [PUIPNIIZUOT JOJ SUID, SUIPCOT
‘ON woisvig
fo 8UuIpDOTwmajsds De bali
Dads
P>0
UTS
i
tea! sera
: Coorer us
Se
I P>0 UZ Tau ‘eT
a ee © ee
oe
UZ
>P>0
ao poe See
x ) ere Ne O=u
c cv
waaPG PEG) uz~sos ai) ere SES
Z) ae eae
[
—9)CI SOS CI+'P)
t
G. Duthie and A. S. Tooth
(
—6—
(9
ee
ia, cancers i@yam =
+'6-)So5('f
(9 +'o-)SO5('9 —'¢-)
ae4(9+
("gu
b)
(“*)
is
ial
(1
——=|V
EV
TW
USS
*7=
uls
vv
'd
ey
usu)
cv
900
soo
=Cy
==
uu Z uu ] ] u
( q—-)SxS
14
G. THOMSON
Ferranti PLC, Edinburgh, UK
ABSTRACT
1. INTRODUCTION
With current design trends in the nuclear and chemical industries towards
higher operating temperatures, the problem of creep in piping systems has
become increasingly important. Furthermore, the necessity for accurate
predictions of the safe working life of components in high capital cost plant
is essential. The fundamental problem in piping structural analyses is to
design a system with sufficient flexibility to contend with thermal loads and
stresses. Smooth pipe bends are one of the essential parts of any pipework
flexibility assessment; being more flexible than an equivalent length of
straight pipe they can provide a significant contribution to the overall pipe
system flexibility.
273
274 G. Thomson
The linear elastic behaviour of thin smooth pipe bends under various
loadings has been extensively examined by many authors in the last seventy
years.’ In 1957, Kachanov? investigated the effect of creep on pipe bends
under in-plane bending. Using a complementary energy method with a
creep power law he derived an upper bound approximate flexibility factor
using a two-term series. Spence later published a series of papers? * giving
more accurate creep flexibility factors employing upper and lower bound
energy analyses with up to five-term series. In his analyses he used a Norton
power law of the form
Zen LEbORY
2.1. Introduction
The analysis presented herein is a combination of the elastic solution of the
present author! together with the method used by Spence* to extend the
elastic theory to the creep regime. Compatible velocity fields are selected
which satisfy the kinematic constraints of the flanged bend problem,
employing trigonometric series. A suitable energy function based on the
In-plane Bending Behaviour of Flanged Pipe Bends Dip
Note that in assumption (b) radius ratio (r/R) itself is not being
neglected.
276 G. Thomson
eC Utena
ig Rlgp~ 880 + 005 0)
) bil OL yee
YeoR\ G6 + ag
. 1fdv 6?w 5
fect) 2
For convenience in the current problem, velocities uv, 0 and w are divided
into two components, the first termed the ‘rigid section velocities’,
associated with the movement of the tube cross-sections relative to the bend
centreline with no change in circularity, and the second are termed
‘distortion velocities’ representing the distortion rates of the initially flat
and circular cross-sections.
Rigid section velocities tig, tg and Wg can be found from U, and V,, the
bend centreline (¢ = +7/2) velocities, i.e.
iif e
Up = U,+),rcos@ je RUF]
tg =V, sing
We =V,cosd (3)
y, is the rotation rate of the bend centreline. Substituting these into the
complete strain-displacement rate eqns(1l) together with the above
assumptions, particularly (e), gives
& = Yop = Ky =0
ees
Ree iaa: : Fae is eU, c
ae (’. 7) (4)
Note that the rotation rate between the ends of the bend at 0= +«/2
153):
The distortion displacement rate component normal to the shell surface,
Wp, can be approximated by the following series which satisfies the relevant
flange and symmetry boundary conditions:
F 7
>=) )Cancos geos? (“2 m=1,3,5..., 9=2,4,6... (6)
m 4
278 G. Thomson
7
r [Cts
R } 00
+)Dai qpcos cu
yale
— gg sin @) cos (mm
3:
DSye
ma (=)
R cos gp cos =)
k, iy >
: ]
nla? = Neos gdcos?(™22)
:
(8
mq
x (i+42 K) rRdg
dd — Mj},
where
n 2n/(n+1)
A=
(5;+ :) (9)
In-plane Bending Behaviour of Flanged Pipe Bends 279
ae 20) JO eae a ow
1 m0
ne «((cosadeos 4 +1 sin gd sin
eos (“=)
ge
+—|— ] cosq@cos
R \ qa
72 i
ak a
mnO\ 12) "+ vi2n
+A by Y Cnt? — 1)cos qd cos? (=) dg dé
Rt & . a @ : a
A= =D, (=) Cut Oe, (=) (11)
3. FLEXIBILITY FACTORS
The end rotation rate for a similar length of straight pipe (Ra), under the
same bending moment, M, is?
_ (Ra\( Me" 1
Ome PNT ig De
|2
0 =4 | (sin p)"*/"d@ (12)
0)
Kral) () 9
‘ (n+1)/n n
(0)
280 G. Thomson
az
PIPE FACTOR (A)
Fig. 2. Comparison of flexibility factors.
which is the end rotation rate of a pipe bend in creep divided by the end
rotation rate of a corresponding straight pipe in creep.
Flexibility factor convergence required about three terms in each velocity
series for A>0.1, «<< 180°, R/r < 10 and n<7. Acceptable convergence
could be obtained with fewer terms inside these ranges.
Figure 2 presents a comparison of the flexibility factors from the present
theory n = 1, with the elastic theory of ref. 8 with y =0.5. Although this
comparison is not strictly valid, since the present theory gives flexibility
rates, its use can be justified by the elastic analogy.'' It can be seen that the
present theory with all its assumptions shows a similar trend to the analysis
of ref. 8 which has more complete displacement series and virtually uses
only the assumptions of thin shell theory for the same problem. However,
the current results are appreciably lower in value, particularly for smaller
radius ratios, R/r, and bend angles, «. Attempts at improving the present
type of analysis for the elastic problem’ have shown that a significant
In-plane Bending Behaviour of Flanged Pipe Bends 281
100
R/r=10
O¢=180"
20
{‘
FACTOR
FLEXIBILITY
1
04 02 05 1 04 02 05
PIPE FACTOR (A)
Fig. 3. Creep flexibility factors for flanged bends with « = 180°.
improvement in the results can be obtained by the removal of the zero direct
shear strain assumption. However, this requires the introduction of a
further independent displacement series for up. Usefully this can have the
same form as used in the current analysis but with independent unknown
coefficients. Introducing these constants in the present analysis is not
difficult but it requires a substantial increase in the solution computation
effort. After this, the biggest improvement can be obtained by use of more
complete displacement series, 1.e. functions which are more capable of
modelling the true displacements. Removal of the zero direct meridional
strain, é, = 0(as used by Karman in 1910) and (r/R) cos @ < | assumptions,
have been shown in the elastic analyses to have little significance. A fully
282 G. Thomson
FACTOR
FLEXIBILITY
—— — SPENCE AG]
01 02 05 1
PIPE FACTOR (A)
Fig. 4. Creep flexibility factors for flanged bends with a=90°.
convergent solution from an analysis of this type may require the use of
about 80 displacement rate coefficients. Numerically minimising a non-
linear problem of this size would be a difficult task, possibly prone to
numerical errors.
Also shown in Figs. 2 and 4 are the results from Spence’s* analysis which
neglected the influence of end constraints. It can be seen that flanges on the
ends of a bend will significantly reduce the potential bend flexibility.
Spence’s results will be more representative of bends terminated by long
tangent pipes.
Figures 3 and 4 present creep flexibility factors of bends with rigid flanges
under in-plane bending for bend angles of 180° and 90°, as determined
In-plane Bending Behaviour of Flanged Pipe Bends 283
00
—— PRESENT THEORY
—-— BOYLE & SPENCE[12
Ss
FACTOR
FLEXIBILITY
i
01 02 05 - 1
PIPE FACTOR ())
Fig. 5. Comparison of present theory with approximate formula of ref. 12.
from the present analysis. From the figures it can be seen that flexibility
increases with reducing pipe factor, A, as was the case in Spence’s analysis,
and increases with increasing radius ratio, R/r, n, and bend angle, «.
Boyle and Spence!’ proposed the following formula:
creep = p" creep
end effects = no re effects
where
p= |Ge etna, Iap effects (l4)
which estimates the flexibility of a bend with end effects under steady state
creep using available data on creep in bends with no end effects and elastic
results for K withn = 1, with and without end constraints. This formula was
284 G. Thomson
4. CONCLUSIONS
Acknowledgements
The author gratefully acknowledges advice and assistance from Professor
J. Spence and Dr J. T. Boyle of the Department of Mechanics of Materials
at the University of Strathclyde.
REFERENCES
ABSTRACT
NOTATION
Area [L’].
Average diameter D = 2R [L].
Modulus of elasticity [KL ~*].*
General notation of normal force [K].
Saha
Theoretical value of the yield force [K].
(xDto,.) of a pipe cross-section.
ss Rocker force (contact force between rocker element and bending
element) [KL *].
* For the sake of simplicity the derived dimension of force is here designated by the
symbol K (K = [MLT] ’).
287
J. W. B. Stark and P. E. de Winter
1. INTRODUCTION
In deep water the loads and deformations that a pipeline can safely
withstand are reduced by external water pressure. The magnitude of such
reduction to a large extent depends upon the pipe parameters and the loads
Plastic Design of Submarine Pipelines 289
and deformations to which the pipe is subjected. The pipe will collapse if the
load carrying capacity or the deformation capacity is exceeded.
This paper presents an approximate method of analysis enabling the
relation between pipe parameters, load combinations and the deformation
capacity to be determined (see Table 1).
Tests have been carried out on model pipes. A brief description will be given
of the test specimen and the test set-up. The results of this test programme
and test results from other sources will be compared with the theoretical
outcome.
2. THEORETICAL MODEL
This part of the model is called the rocker element. The geometrical non-
linear deformation in the circumferential direction can now be determined
very simply with only one parameter: the ovalisation angle f (see
appendix). By virtue of the four hinges the rocker element can be deformed
without any resistance. Therefore, a flexural element, representing the
circumferential stiffness, is also needed.
The stiffness of the pipe cross-section can be modelled either by
rotational springs at the hinges or by a bending element. Professor J.
Witteveen of Delft University proposed the use of the bending element
(Fig. 2) because it allows, in principle, for lower bound solutions.° It is
essential to know that the rocker element is only connected at the hinges to
(a) (b)
Fig. 1. (a) Cross-section after collapse. (b) Cross-section at the onset of collapse.
(c) Theoretical idealisation.
31}DWeUIy
japow UJIM
D
Bulpueq juewe}e
J9HDOI JUBWI]e Buipueq 1UeWej}a
Plastic Design of Submarine Pipelines
“BY«°*7 OneUIOUTY
‘JEpoul
291
292 J. W. B. Stark and P. E. de Winter
3. TEST PROGRAMME
Tests have been carried out in order to verify the results of the theoretical
model. Test results from other sources have gratefully been used when
available and permitted.
In order to minimise the influence of geometrical imperfections some of
the test specimens used at IBBC-TNO have been machined carefully, first
internally and then externally, on a lathe, thus ensuring that the test
specimens had a truly circular cross-sectional shape with constant wall
thickness.
Variations in both wall thickness and diameter were less than +0.02 mm.
The mechanical properties and dimensions of the specimen which have
been tested are presented in Table 2.
All tests specimens were amply instrumented. Specially designed
callipers were attached to all test specimens to measure the change of
diameter during loading.
Plastic Design of Submarine Pipelines 293
S: specially machined test pipe; H: seamless commercial pipe; W: welded pipe; o,:
yield stress; og: failure stress; ¢,: strain beyond which work hardening occurs; é.:
strain at failure.
dead weight
compensation
pressure vessel
test. An additional advantage is that the need for a large pressure vessel is
greatly reduced.
With the test set-up as described above, the theoretical model can be
verified adequately. However, the distribution within the pipe does not
exactly correspond with the actual situation in the sea. In Fig. 5 the test
situation (no endcap) is compared with the actual situation in the sea
(endcap).
For this reason IBBC-TN O is considering constructing a new test facility
in which the specimen is installed inside the pressure vessel. The diameter
will be 1200mm (48”) and the operating pressure of the vessel will be
approximately 250 bar. An additional advantage is that larger specimens
can be tested to a larger curvature.
Plastic Design of Submarine Pipelines 295
test facility
R me
nm elastic
Uv ae)
ov
[P/Pp]
hydrostatic
pressure
Po ovalization e
COLLAPSE
OP RESSURE [P/Pp]
00 02 04 06 08 40 roundness [35]
Fig. 7. Collapse pressure of out-of-round tubes.
296 J. W. B. Stark and P. E. de Winter
4.1. Theory
In this section the theoretical model, described in Section 2 and the
appendix, is used to derive an interaction formula for external pressure in
combination with axial tension. The effects of initial out-of-roundness will
be included. For this loading case the equilibrium equations for the rocker
element as given in the appendix can be reduced as follows:
M=0
F=|o,,dA
8pR?B =4fR (1)
The elastic pressure—ovalisation curve (Fig. 6) can be found by combining
the equations for the rocker and bending elements. The same applies for the
[o)
©
(=)
— Sa
S
wu
= ©
q
oO
S)
je)
@
Oo
©
W225 26
- F/Fp : -10%. - 90%
ss A : Ys = 450 N/mm2 —
: Ys = 700 N/mm2
Oo
2 : Ys = 900 N/mm2 —4
eo}
5.1. Theory
One of the loading cases that occurs in connection with the installation of
subsea pipelines is the load due to external pressure and bending. The pipe
emerges from the pipelaying vessel at a slope (Fig. 9) and has to curve to a
horizontal position in order to settle down on the sea bed. Deformation
associated with this change of direction is imposed upon the pipe. In
connection with this imposed character it is of considerable importance to
know what maximum curvature the pipeline can safely absorb.
(1) the external pressure (P) is kept at a constant level, while the
curvature is increased from zero to the critical value;
(2) a three-dimensional yield surface w for bending moment (M),
external pressure (P) and rocker forces (f) is adopted;
(3) the normality principle is used for calculating the critical curvature;
(4) this analysis is performed on pipes with endcaps (hydrostatic
pressure).
Cals
s Ghana
nfo)
wi -a(f)
The yield surface has a boundary for
Pp
In the analysis given here the pressure P is kept at a constant level and the
curvature is increased from zero until collapse of the pipe occurs. This
means that it will be sufficient to consider a straight section through the
Plastic Design of Submarine Pipelines 299
115 P/Pp
“Zk
a 78 aV3t =1-3(
Fig. 10. Yield surface.
4 M/M>
ee)
t/tp
ftp
Fig. 11. Loading path by constant pressure.
Oo Re K
Se s-3(FY (4)
Py
Plastic Design of Submarine Pipelines 301
For point C, where the critical curvature is attained, the following equation
for plastic collapse curvature can be derived:
(ze) to
re"
jo et
t 4\P
Ker = 4/3 55 ‘
The above formula applies to pipes with pressure on the endcaps. For the
‘no endcap’ situation, an analytical solution is not available.
09 4
0. 2s test results ee
| | r) D/t =25 pe A
07 o Ot =32
8
g
jo)
OoWw
Ooi)
external
CP/Pp]
pressure
——»
oO=
Oo RO 4 6 8 Toe WA SIS)
——+» critical curvature [C/ce]
rigid plastic interaction curve
Fig. 12. Rigid plastic interaction curve.
302 J. W. B. Stark and P. E. de Winter
2 4 6 8 10 12 [C/Ce]
test
Fig. 13. Bending + external pressure: Shell test data vs. theory (formula).
aS ee
6 8 ! 12 [C/Ce]
Fig. 14. Bending + external pressure: Shell test data vs. theory (computer).
Plastic Design of Submarine Pipelines 303
M/Mp
P/Pp
t/fp
Fig. 15. Loading path by constant curvature.
304 J. W. B. Stark and P. E. de Winter
curvature (*/x,)
O 510" 15° *20'" 25°77 30
Critical pressure-curvature combinations
Fig. 16. Critical pressure-curvature combinations.
the loading path is drawn schematically on the yield surface (route O-A-B).
As with other structures, the collapse load here is also path dependent.
The computation of the curvature up to point B is done in the same way
as described in section 5.1 for pressure P = 0 and f=f,. At this point the
curvature is kept at a constant level, while the pressure is increased from
zero until collapse. The path over the yield surface depends on the
equilibrium equation of the rocker element and the normality principle.
With the normality rule the deformation increments on each point of the
yield surface can be established:
Since 0/0M #0, the scaling factor 4 has to be equal to zero. Proceeding
Plastic Design of Submarine Pipelines 305
over the yield surface, while the curvature is kept at a constant level, the
other deformations are also maintained constant.
b= k= A= 0
On the boundary of the yield surface the ultimate pressure is reached
(point C).
The ultimate pressure in point C can be calculated using the
deformations f and x, attained in point B. In Fig. 16 the computed
instability curves for xP and P-« are plotted together. In this particular
case there is considerable influence (D/t = 25, ¢, = 320 N/mm’).
7. CONCLUSIONS
Acknowledgement
This contribution forms a part of the Marine Technological Research
Programme sponsored by the Dutch Government. The work is carried out
jointly by the Institute for Building Materials and Building Structures of
the Netherlands Organisation for Applied Scientific Research (IBBC-
TNO), Protech International and the Universities of Delft and Eindhoven.
The authors wish to thank Nippon Steel Corporation and Shell SIPM for
making their test data available to us.
306 J. W. B. Stark and P. E. de Winter
&q0:dx
SS,
axial deformation due to axial force
APPENDIX
For the kinematic model, three degrees of freedom will be considered (see
ig Wy:
—ovalisation (f)
—axial deformation (¢,0)
—curvature (k)
Rocker element
Geometry (see Fig. 18)
Linearised:
AA =4R’B? 1° variation 6AA = 8R?BdB (14)
Kinematic equations
Strain in axial direction:
Static equations
Equilibrium of the cross-section is expressed using virtual work theory (see
Figs 17):
Jo,.66,,4dA = MOx + Foe,)
+ POAA — 4fR6B (17)
After substitution of the variations derived above, an equation is found
which is valid for random variation of the three degrees of freedom.
Therefore, three equilibrium equations may be written:
M =\{o,,2(1 — B)dA (18)
F=\q.dA (19)
8PR?B = —|o,,kzdA +4fR (20)
The first two equations can be easily recognised. The external bending
moment (M) equals the resulting bending moment caused by axial stresses
in the deformed model. The same applies for the axial force (F). The third
equation shows the connection between both external pressure and bending
with the ovalisation of the tube.
Bending element
The bending element is loaded only by those (opposite) contact forces (/)
which act on the rocker element. Stiffness of the bending element is
nevertheless influenced by the load acting on the rocker element. The
starting point is a linear analysis in which the displacements on the
horizontal and vertical axes are of equal magnitude.
Static equation
Because of symmetry it is sufficient to consider the equilibrium of a quarter
of the cross-section (see Fig. 19):
(ae (21)
Plastic Design of Submarine Pipelines 309
US
iA
Mr
Mr, Ve
ly,
Fig. 19. Equilibrium of the bending element.
For rigid plastic material behaviour, a simplified yield function is used. The
yield function describes the interaction between bi-axial bending and
extension.
On the assumption of three zones in an elementary block of material,
each satisfying the Von Mises yield condition, the following equation can
be derived:
Cl rastal ealaarss te
2 2 2 2 2
m, mM, mM, ny Ny ny
where
m, = 1/4176, ne = 10; n,=—PR
From the tests observed and from the collapse mode in the theoretical
model, it is clear that the change of curvature in the pipe wall material in the
axial direction is negligibly small in comparison to the curvature change in
the circumferential direction.
K, <K,
In the further analysis it is assumed that
K, =0
310 J. W. B. Stark and P. E. de Winter
vah-EJ-QPEII so
The bending moment in the circumferential direction now becomes
The normality rule is applied for deriving the deformation in the plastic
area (see interaction between bending and external pressure). Based on the
relationship presented above, the interaction between bending, external
pressure and axial force can be derived.
REFERENCES
W. J. M. STEEL
Britoil PLC, Glasgow, UK
and
J. SPENCE
University of Strathclyde, Glasgow, UK
ABSTRACT
1. INTRODUCTION
In the early 1970s techniques had been developed for routinely laying
offshore pipelines in water depths up to about 122m. However, the
capability has since been extended as the search for crude oil and natural
gas has moved into deeper water. Examples include a line in 265 m of water
off the California Coast, the planned construction of a line across the
313
314 C. J. Tay, W. J. M. Steel and J. Spence
tensioner
barge“ y
\ stinger anchor cable
\
pipeline
The buckle, once initiated, may propagate along the pipe under the
influence of the pressure alone even into regions free of bending. The
financial and time penalties resulting from such a buckle are severe. It is
therefore of paramount importance to have a comprehensive understand-
ing of the bending/pressure buckling behaviour to ensure the integrity of
offshore pipelines during the laying process.
Previously a great many investigations had been carried out on the
buckling problems of cylinders under different loading conditions. The
contributions of Ades,! Seide and Weingarten,” Lee,? Wilhoit* and
Berdnikov and Nozdrin° reported work being carried out on pure bending
of cylindrical shells. Others such as Stokey er al.,° Gerdeen et al.’ and
Larson® reported work on pipes subjected to combined internal pressure,
bending, axial load and torsion. However, relatively little attention has
been focused on the plastic behaviour of pipes under combined bending
and external pressure.
In recent years experimental work on the buckling of pipes under
combined bending and external pressure has been reported by Johns et al.?
on 51mm diameter pipes of an unspecified grade of steel, and by
Broussard!° on 152mm pipes with a d/t ratio of 24. The pipeline material
was unfortunately not specified and also the detailed development of the
theoretical analysis was not given.
An experimental programme on the bending/pressure behaviour and
buckle propagation has been undertaken by the present authors under the
sponsorship of BP (Trading) Ltd. The test facility for this investigation
comprised a pressure vessel together with suitable instrumentation and the
supporting systems for externally pressurising and/or bending the pipes.
The pressure vessel, which was 1.07 m in diameter and 6.7 m long, allowed
the model pipes (approximately 1/5 scale) to be subjected to external
pressure while being simultaneously subjected to bending by means of two
external 200kN hydraulic rams. The rams were mounted on an external
test rig and acted through glanded nozzles on top of the vessel.
In addition, a number of small scale tests were conducted on mild steel
pipes using the test facility. This background stimulated the theoretical
work outlined herein. The detailed results of the test programme are not
reported here. Theoretical results of the work on some aspects of buckle
propagation and sleeve arrestors is reported separately.'? For interest a
typical buckled cross-section is shown in Fig. 2.
2. THEORETICAL ANALYSIS
€, €z
-€;5 mes
pr P
C= = Uta (1)
318 C. J. Tay, W. J. M. Steel and J. Spence
2.3. Compatibility
From assumption (i), the axial strain may be expressed as
e, =e tee (6)
Inelastic Response of Pipes 319
Stress
:
bi-linear idealisation
tan” E,
\ actual
Strain
de?
Zod,
oO Zz
(7)
Oi (10)
320 C. J. Tay, W. J. M. Steel and J. Spence
where g, is the yield stress of the material. In eqn (10) dé? is the effective
plastic strain increment and may be expressed, in terms of the strains, ina
similar form to eqn (11).
A relationship between effective stress and plastic strain may be obtained
from a simple uniaxial test:
o = g(e") (12)
As stated earlier, the material is assumed to have a bi-linear constitutive
relationship which, from Fig. 5, is
Eae er :
ars GSB,
relVS laoS
edOe
5 eye be
aes
(13)
=é* +e?
Then
ag — 109 — %)
‘ (14)
_ dé
a
d
dee =
E
ue (15)
dep = 22 260
0
+ &)] do
Ey
(16)
Inelastic Response of Pipes 321
eee
jon 4) ] -a|; at -|- és on |
[o, — v(o, +
by (Per eC Es
322 C. J. Tay, W. J. M. Steel and J. Spence
When a = 9, the value of o, is obtained from eqn (11) while o, and a, are
given by eqns (1).
In the case of pure bending, the external pressure is zero and thus
eqn (21) becomes simply
ana(ste)-%
=
1
—+—)]-—
1 oO
(23)
23
M 2
—,— =—JS,sin9dé (25)
REO SH
where
saeoO Ox
The axial strain in the elastic and plastic regions of the pipe, eqns (4a) and
(21) respectively, may also be non-dimensionalised as follows:
3. RESULTS
Typical results obtained from the analysis are shown in Figs. 6-8. As
mentioned above, the material of interest for pipeline evaluation is typically
API 5L Grade X65 which has a nominal yield stress of 450 MN/m?
(65 000 Ibs/in?). The dimensionless parameters in Figs. 6-8 were therefore
chosen relative to this yield stress and a Young’s modulus of 210 MN/m?
(30 x 10° lbs/in?). If these basic material properties apply, Fig. 6 then
shows the effect of varying the tube size when the external pressure is
3.45 MN/m? (500lbs/in*) and the tangent modulus is 3.45 GN/m?
(0.5 x 10° Ibs/in?). This value of the tangent modulus was judged to be
reasonably representative of Grade X65 material. For one particular size
(r/t =15), Fig. 7 illustrates the variation in the moment versus strain
relationship resulting from two different external pressures. For the
324 C. J. Tay, W. J. M. Steel and J. Spence
M Oy
p = 130 Ee 60
9 == 3 & 461-5fia
, ==
O “
1:0
0-75
05
0:25
0 2 4 6 8 10 12
Extreme fibre strain x10"
Fig. 6. Moment vs. extreme fibre strain.
= 461-5
10
0-75
0-5
0-25
0 Z 4 6 8 10 12
Extreme fibre strain x10°
Fig. 7. Moment vs. extreme fibre strain.
4. CONCLUSIONS
v te
petri Dy = 130 lar
ago =Ors a
ul
0-75
05
025
0 2 4 6 8 10 12
Extreme fibre strain x10>
Fig. 8. Moment vs. extreme fibre strain.
REFERENCES
1. Apes, C. S., Bending strength of tubing in the plastic range, J. Aero. Scis., 24
(1957), 605-10.
2. SEIDE, P. and WEINGARTEN, V. I., On the buckling of circular cylindrical shells
under pure bending, Trans ASME, 83, Series E, J. of Applied Mechanics (1961)
112-16.
Inelastic Response of Pipes 327
. Leg, I. H., ‘Buckling and ovaling of pipeline’, M.S. Thesis, Rice University,
Houston, Texas, May 1971.
. Wituoit, J. C. JR, MERwIN, J. E. and Jirsa, J. O., Critical plastic buckling
parameters for pipes in pure bending, ASME Publication paper, 72-pet-29,
1972.
. BERDNIKOV, Y. N. and Nozprin, A. S., Stability of cylindrical shell in the
elasto-plastic region with pure bending (yield theory), Soviet Aeronaut., 17,
No. 4 (1974) 11-17.
. STOKEY, W. F., PETERSON, D. B. and WunpkerR, R. A., Limit loads for tubes
under internal pressure, bending moment, axial force and torsion, Nuclear
Engineering and Design, 4 (1966) 193-201.
. GERDEEN, J. C., SIMORREN, F. A. and Hunter, D. T., Large deflection analysis
of elastic-plastic shells using numerical integration, AJAA Journal, 9, No. 6
(June 1971).
. Larson, L. D., ‘Inelastic response of pressurised tubes under dynamic bending
and torsional loads’, PhD thesis, Mech. Eng. Dept., Carnegie—Mellon
University, 1973, Microfilm order No. 73-22872.
. JOHNS, T. G., MESLOH, R. E., WINEGARDNER, R. and SORENSON, J. E., ‘Inelastic
buckling of pipelines under combined loads’, OTC paper No. 2209, 1975.
. BROUSSARD, D. E., Deepwater pipelaying research under way, Oil and Gas
Journal (June 1975).
. Tay, C. J., ‘Inelastic collapse of pipes under external pressure and bending’,
Doctoral dissertation, University of Strathclyde, June 1979.
. STEEL, W.J.M.and Spence, J., On propagating buckles and their arrest in sub-
sea pipelines, Proc. I. Mech. E. London, 197a (1983) 139-47.
. FRANZEN, W. E. and STokey, W. F. ‘The elastic-plastic behaviour of stainless
steel tubing subjected to bending, pressure and torsion’, 2nd Int. Conf. on
Pressure Vessel Technology, Part 1, Design and Analysis, 1973, pp. 457-67.
. JOHNSON, W. and MELLoR, P. B., Plasticity for mechanical engineers, New
York, Van Nostrand, 1962.
. Hopce, P. G. and ROMANO, F., Deformations of an elastic-plastic cylindrical
shell with linear strain hardening, J. Mech. and Phys. of Solids, 4(1956) 145-61.
ae eee
F Af tay
eye #7 Sognkion
ee a
gee meteooi valeer as ‘ "
on aw a yet vet rea vein
wa i ore ery a¥ aa “9 W wwett..o vw eS
- s ins keg emai DOLE ,crea nies Sipe cath ty amt t 4
; Vets As ey dvinwel Gee cme wryaatSs, 9 a) -ereae oe Yi
4 4
413 maf — Teve6 ama) .7
ied gay: tri, (meri Hew waitbi ptr ethan Jb 73 eet
: la are mit onjbtint hy ees) ae) Vier
i reed Re orean?, A4RSA2st «4 As d,
ae : ik Ai) gO wee) YW iatkh A &
or chien’ a ined A Pig Liha ae pines 2 a ae
Bt, Sow! [A aeety gaaeS reTa ace
ah SOF! epygh re Nye).Lc eet 98)er ord lesen ted
ames tee, in” “tind i
my Ye ities ae ; H
ies ae
oprk Anu ne
; hbhaseaace a me ys che. AnLe aa
Rees
tadoliakg atin SD % gia :
Lites s wy aes ie
: - : =P
® , I el rr. =
She iy & _ om
o_o v
4% + lewes
. by rp
17
ABSTRACT
NOTATION
a Length of plate.
b Breadth of plate.
t Thickness of plate.
E Young’s modulus.
K, Design coefficient for longitudinal direction.
K, Design coefficient for transverse direction.
4 Design coefficient.
B Non-dimensional plate slenderness (=6/t,/09/E).
69, Geometric imperfection amplitude of a single half-wave component.
329
330 A. F. Dier and P. J. Dowling
1. INTRODUCTION
Since the mid-nineteenth century plates have been used extensively in steel
structures as elements participating in the primary load carrying action.
The first solution to the linear plate buckling equation for a simply
supported plate under edge compression was due to Bryan’ in 1891. Since
then many solutions, accounting for different boundary conditions, types
of loading, etc., have been found and a large number of these are collected
together by Bulson.? For stocky plates the elastic critical stress lies above
the yield stress but sensible results can be obtained by using a reduced value
of Young’s modulus in the critical stress expressions. During the 1940s the
deformation (history independent) theory of plasticity was used by workers
such a Bijlaard,* Ilyushin* and Stowell® to evaluate the reduced tangent
modulus. The more correct incremental flow theory of plasticity was first
applied by Handelman and Prager.°®
Further knowledge of elastic plate behaviour was obtained by solving the
large deflection equations of von Karman which can account for the post-
buckling reserve of slender plates. Among the notable early workers in this
area are Levy,’ Coan,® Hu et al.? and, more recently, Yamaki.!° These
‘classical methods’ gave way to less exact and numerical methods such as
Rayleigh-Ritz, finite differences, finite elements, etc. The perturbation
approach'! has been used to find design coefficients for uniaxially and
biaxially compressed plates.'* By assuming a plate fails when the midside
edge stress reaches yield an estimate of the load carrying capacity of the
plate may be found. Much of the information obtained from the above
analyses and other works, including experimental, has been used to derive
effective width expressions for plates in compression. Faulkner!? has given
The Strength of Plates Subjected to Biaxial Forces 331
2. METHOD OF ANALYSIS
3. NUMERICAL RESULTS
he
q/0 = €
SOIl}DJ
We. = ave CS)
A3:x3
yoy uOlyDUD|dxa
JO
SaAUND
Y PpuDG 22S 1X2]
A. F. Dier and P. J. Dowling
ON
oe
poy
JO UOONIWSUOD
JO UOTOVIIUT
"SAINI
‘|
3q
The Strength of Plates Subjected to Biaxial Forces 335
Transversely, the width of the tensile residual stress block is taken as half
that of the longitudinal block for 3:1 plates and for plates in the transverse
strength study. However, in the investigation of square plate biaxial
strength, the widths of the tensile blocks are taken to be the same in both
directions. This makes the x and y directions indistinguishable for a square
plate and thus only approximately half the number of program runs are
needed to generate interaction curves leading to savings in computer costs.
Due to out-of-plane components of the residual stresses the plate is not
initially in equilibrium. Before the plate is loaded it is allowed to find its
equilibrium configuration using the constrained edge condition. It is
generally found that the out-of-plane components magnify the geometric
imperfection and so the initial magnitudes are factored down and a second
equilibrium configuration found. If the resulting magnitude is now close to
the desired value, loading can begin; otherwise the process is repeated until
it does become acceptable.
positions the long edges take up at the end of the first stage, meanwhile
keeping the short edges fixed. The second loading sequence [B] is a reversal
of the above, i.e. transverse loading followed by longitudinal loading to
failure. The stress paths resulting from the two load sequences have been
plotted with those obtained from proportional displacement loading
(Fig. 1).
In case [A], longitudinal loading is applied up to ¢,/é) =0.5 (i.e. that
displacement which would cause a strain of half yield in a perfect plate). At
this stage the classical three half-wave buckling mode has developed and as
a consequence of this the longitudinal stress initially increases with the
application of the transverse load. However, as the plate is still
substantially in the elastic state, it quickly snaps through and takes up the
single half-wave mode of deformation. From this point the longitudinal
stress decreases with increasing transverse load. As can be seen the stress
path grazes the interaction curve produced from the proportional
displacement stress paths. At the point of contact the ¢,:¢, ratio of plate [A]
is approximately 1/4:1 and it is interesting to note that this ratio is the same
as that of the proportionally loaded plate which fails at almost the same
place. However, the absolute magnitudes of ¢, and ¢, are not the same for
the two cases, being 18° more for the non-proportionally loaded plate.
In case [B] transverse loading is applied up to ¢,/é) = 0.5 with the single
half-wave mode predominating. This mode is maintained while longitu-
dinal loading is applied, although there is a tendency for the three half-wave
mode to develop. The maximum longitudinal stress occurs outside the
drawn interaction curve and lies very close to the 1:1/2 (¢,:¢,) ratio stress
path failure point. The interaction curve has not been extended to this point
as the single and three half-wave buckling modes interact to anomalously
strengthen the plate in this region. Small deviations in the initial
imperfections could bring the failure point back towards the interaction
curve as drawn. Here again, the ¢,:e, ratio of plate [B] at collapse is
approximately the same as that of the associated proportionally loaded
plate (¢,:¢, = 1:1/2).
Further evidence for the soundness of the approach will be found in the
following section where, for the special case of loading along the stress axes
(i.e. proportional stress loading), approximately the same strengths are
predicted.
for plates with an aspect ratio of 3 are shown in Fig. 3. For both aspect
ratios, increasing slenderness causes a steady decrease of the proportion of
squash that a platecancarry where at least one of the stress components is
compressive (i.e. positive stress values).
In the absence of transverse loading (a, = 0), approximately the same
load can be carried in longitudinal compression by plates of the same
slenderness but of different aspect ratios. Physically, this may be explained
by considering the 3:1 plate, which buckles into three half-waves, as three
338 A. F. Dier and P. J. Dowling
square plates arranged end to end. When only transverse loading is present
the plates buckle into a single half-wave mode, the longer plates being
weaker due to column-like behaviour in the mid-field of the plate. For
combined loading on 3:1 plates, the transition of the single half-wave and
the three half-wave buckling modes is indicated by sharp changes of
curvature in the interaction curves (Fig. 3). It may be noted that these
transitions do not occur on the radial line predicted by the classical linear
solution.
The effect of aspect ratio on transverse strength is further illustrated in
Fig. 4. These results were generated by loading in the transverse direction,
keeping the short edges constrained. As can be seen, most of the fall-off in
‘ 4/4 (sw)
—~_ OZ
Ss
~
SS
SS
ASS
SS OE
SS
SS
=
~
SS
Ora
a
aS
Ss
ae
The Strength of Plates Subjected to Biaxial Forces
a<
=$ Coon
-<
OlL
wW— =
q/p
4D9dsSD Ol}Du
*BIY«“p OYA
Jo oodse ONIUO ssJoAsuRI] “YuINs
339
340
q/o
=)
bib
JO
L:€
(241d
U}bHu2u4s
|DUIPNz1BUC|
AD
b-s
A. F. Dier and P. J. Dowling
Lie Lia
L:G
OL 4/Q (SW)
SO OL SL Oz Guc O-€ Se ov d
“Bly*¢ OYA
JO SsaUIApUITS
UO ISIOASURI] ‘YIBuas
The Strength of Plates Subjected to Biaxial Forces 341
strength occurs for the lower aspect ratios. Also shown in Fig. 4 are the
limiting strengths for infinitely long plates.
The results of Fig. 4 are replotted in Fig. 5 using the slenderness for the
abscissa. The strength values from the interaction curves for 3:1 plates are
also indicated on this figure. The good agreement between these values
(produced by proportional displacement loading) and the constrained edge
results (proportional stress loading) establishes confidence in the method of
interaction curve construction and in the use of constrained results for the
transverse strength study.
Although Figs. 4 and 5 indicate the sizes of interaction curves for
different aspect ratios, they do not give any indication as to the shape of the
curves. Figure 6 shows the interaction curves for three aspect ratios
° PduUIDIZSUOD
S3|NSDU
x WOU4 UO!}DD42\UI
S$2AUND
—
sbi4zg pud ‘€ $2u-pud
62)(
A. F. Dier and P. J. Dowling
—
Moy
|
ybiu
aBouaao
————
SUOI}I2jJyadw
ele| lor o3|
O
yi
4) bo.
_ | oat 3/q (sw)
suoljo2ajuadwi
SO OL L S O~ Ve Qe Ge Ov 9)
suol}a2syuadw1
BIgLZ poy
Jo suoNdajJoduI
uo ‘yiBuaNs
=q/b
'€
The Strength of Plates Subjected to Biaxial Forces 343
(a/b = 1, 3, 5) for two slendernesses (b/t = 30, 80m.s.). For the longest
plates, a combined single and five half-wave imperfection mode has been
assumed. As can be seen there is a considerable change in the shape of the
interaction curves for aspect ratios between | and 3. However, the curves
for plates of aspect ratios 3 and 5 are similar in shape for any given
slenderness.
Before any design procedure based on theoretical work can be used with
confidence, it must be shown that the theory is in accordance with
experimental results. Some carefully conducted tests carried out at
Cambridge by Bradfield?* have shown that elasto-plastic large deflection
344 A. F. Dier and P. J. Dowling
jDJUaW142dxKyA
ere yuzuoduos
|
( An
¥€) patjddp
¥SD|
c-&G
6-LE
$24 1/9 pel
a= -¢
Y
@
‘
s3jnsau
a (a) oy To
LZiGOm
4
‘syuouttiodxs
q/p= ‘¢
uosuedWOD
YIM
“BI°*g
Me 4 S
=_ 2
Ko
The Strength of Plates Subjected to Biaxial Forces 345
5. DESIGN CONSIDERATIONS
It is usually found that not all four points can be satisfied simultaneously
and a compromise must be sought, especially with regard to the first point.
sanjD~A
WOJ}
PaUIDIGO
UOl}DID4I4zUI!
S2AGND
——
S2UI|
+ Sag OOS SaN|IDA
JO} |:€ $2yDjd
——
P2}D|Oduazu;
Xy
j\D‘'
SoOlyDu
yoadsp
A. F. Dier and P. J. Dowling
that can represent the interaction curves (and which are also valid for
tensile components), they would undoubtedly be too complicated for use in
a design code (see point (i) above). However, an examination of Figs. 2, 3
and 5 suggests an expeditious method to account for tensile components;
that is to use horizontal or vertical lines as appropriate to approximate the
interaction curves in these regions. For square plates (Fig. 2) the ordinates
and abscissae of the lines would be safely given by the values indicated
where the interaction curves cross the stress axes. In some instances for 3:1
and 5:1 plates (Figs. 3 and 6), slightly reduced values would have to be used
to ensure conservatism in the combined tension—-compression regions.
It is widely recognised that for plates of the same slenderness,
approximately the same longitudinal strength is expected, no matter what
the aspect ratios may be (see Fig. 6 or compare Figs. 2 and 3). It is thus
common to use a single design curve cavering all aspect ratios for the
longitudinal strength. In the present case the curve would be that obtained
from a consideration of the longer plates, again to ensure conservatism.
Based on the present numerical results and the above proposed format, the
values of the design coefficients are shown in Fig.9. Some design
coefficients of BS 5400°’ are also indicated and these shall be referred to
later.
An interaction formula is needed for the compressive quadrant to
complete the design curve. Figure 10 shows the normalised interaction
curves (with respect to the K values of Fig. 9) in the quadrant for square and
3:1 plates. The curves do not pass through unity on the axes because the K
values, on which the normalisation is based, are lower than the
corresponding o™/o, values of Fig. 5. The general interaction formula
proposed in ref. 23, which covers other commonly used formulae, is of the
form
although there is no cross product term (i.e. 7 =0). A safe value of 7 for
square plates calculated from the curves of Fig. 10 is 0.45.
Having decided on a value ofn the values of K for stocky plates may now
be found. As an example the K value for the b/t=20 plate may be
considered. The value of o,/a) when o, =, is 0.88 and thus K is given by
or
The Strength of Plates Subjected to Biaxial Forces 349
The one major disadvantage in this approach is that the designer may be
confused as to how a plate can have K values greater than unity implying
strengths greater than squash (see point (iii) above). It is, of course, just a
design artifice for biaxially loaded plates and care in the wording of the
rules should ensure that it is not misapplied. The corresponding
normalised curves for 3:1 plates (Fig. 10) are more widely spaced than
those of the square plates and they are non-symmetric about the diagonal.
Parabolic*® and skew parabolic”? forms have been proposed for which
a = 1 in eqn (1). However, for the sake of simplicity, it is felt that keeping
a = 2 is still satisfactory. A function for 7, dependent on f, can be found to
take account of the different shapes of the interaction curves.
For longer plates, the same values of y as those of 3:1 plates can be used
as the shapes are very similar (Fig. 6), and using, of course, the appropriate
K values. It is suggested that for plates of aspect ratios between 1 and 3, a
linear interpolation of 7 should suffice.
slender plate region is not steep enough and thus gives rise to high
coefficients for transversely loaded long slender plates. It may be noted
that, by using this approach, very long slender plates are out of the range of
the BS 5400 curves.
6. CONCLUSIONS
(1) Interaction curves for square plates and plates of 3:1 aspect ratio
have been presented for biaxial loading.
(2) Although longitudinal strength is not greatly affected by aspect
ratio, the transverse strength is. Curves have been presented
showing the transverse strength of plates of practical aspect ratios.
(3) Adesign method covering the range of possibilities from combined
compression—compression to compression-tension has been
outlined.
(4) Thenew UK bridge code has been shown to be deficient for cases of
combined compression and tension and also needs a re-
examination for transversely loaded plates of large aspect ratios.
REFERENCES
1. BRYAN, G. H., On the stability of a plane plate under thrust in its own plane
with application to the buckling of the side of a ship, Proc. London Math. Soc.,
22 (1891).
2. BuLson, P. S., The stability of flat plates, London, Chatto & Windus, 1970.
3. BIJLAARD, P. P., Some contributions to the theory of elastic and plastic
stability, International Association of Bridge and Structural Engineering, 8
(1947).
4. ILyusuin, A. A., ‘Stability of plates and shells beyond the proportional limit’,
National Advisory Committee for Aeronautics, T.M.116, 1947.
5. STOWELL, E. Z., ‘A unified theory of plastic buckling of columns and plates’,
National Advisory Committee for Aeronautics, T.N.1556, 1948.
6. HANDEELMAN, G. H. and PraGer, W., ‘Plastic buckling of rectangular plate
under edge thrust’, National Advisory Committee for Aeronautics, T.N.1530,
1948.
7. Levy, S., ‘Bending of rectangular plates with large deflections’, National
Advisory Committee for Aeronautics, T.N.846, 1942.
8. Coan, J. M., Large deflection theory for plates with small initial curvature
loaded in edge compression, Journal of Applied Mechanics, 18, No. 2 (1951)
143-51.
352 A. F. Dier and P. J. Dowling
>: Hu, P. C., LuNpquist, E. E. and BATporz?, S. B., ‘Effect of small deviations
from flatness on effective width and buckling of plates in compression’,
National Advisory Committee for Aeronautics, T.N.1124, 1946.
10. YAMAKI, N., Postbuckling behaviour of rectangular plates with small initial
curvature loadedin edge compression, Journal of Applied Mechanics, 26, No. 3
(1959) 407-14. Continued in 27 (1960) 335-42.
ll. THOMPSON, J. M. T. and WALKER, A. C., A non-linear perturbation analysis of
discrete structural systems, International Journal of Solids and Structures, 4
(1968) 757-68.
1 WILLIAMS, D. G. and WALKER, A. C., Explicit solutions for the design of
initially deformed plates subject to compression, Proceedings of the Institution
of Civil Engineers, 59, Part 2 (1975) 763-87.
iley. FAULKNER, D., A review of effective plating for use in the analysis of stiffened
plating in bending and compression, Journal of Ship Research, 19, No. 1 (1975)
1-17.
. Ruopes, J., Effective widths in plate buckling, Developments in thin-walled
structures 1, J. Rhodes and A.C. Walker (eds.), London, Applied Science
Publishers, 1982.
. Dow.inc, P. J., ‘Some approaches to the non-linear analysis of plated
structures’, Symp. on Non-linear Behaviour and Techniques in Structural
Analysis, Transport and Road Research Laboratory Supplementary Report
164 UC, 1974, pp. 86-93.
. CRISFIELD, M. A., ‘Large-deflection elasto-plastic buckling analysis of plates
using finite elements’, Transport and Road Research Laboratory Report LR
5931973.
. HarpinG, J. E., Hopss, R. E. and NEAL, B. G., The elasto-plastic analysis of
imperfect square plates under in-plane loading, Proceedings of the Institution
of Civil Engineers, 63, Part 2 (1977) 137-S8.
. Frieze, P. A., Hospps, R. E. and Dow Lina, P. J., Application of dynamic
relaxation to the large deflection elasto-plastic analysis of plates, Computers
and Structures, 8 (1978) 301-10.
. Littte, G. H., Rapid analysis of plate collapse by live energy minimisation,
International Journal of Mechanical Science, 19, No. 12 (1977) 725-44.
. BRADFIELD, C. D. and CHLADNY, E., ‘A review of the elasto-plastic analysis of
steel plates loaded in in-plane compression’, Cambridge University Report
CUED/D-Struct/TR 77, 1979.
. FRIEZE, P. A., DOWLING, P. J. and Hosss, R. E., Ultimate load behaviour of
plates in compression, Stee/ plated structures, P. J. Dowling, J. E. Harding and
P. A. Frieze (eds.), London, Crosby Lockwood Staples, 1977, pp. 24-50.
Pd, HARDING, J. E., Hosss, R. E. and NEAL, B. G., Ultimate load behaviour of
plates under combined direct and shear in-plane loading, Steel plated
structures, P.J. Dowling, J. E. Harding and P.A. Frieze (eds.), London,
Crosby Lockwood Staples, 1977, pp. 369-403.
2% VALSGARD, S., ‘Ultimate capacity of plates in biaxial in-plane compression’,
Det Norske Veritas Report 78-678, 1978.
24. LittLe, G. H., The collapse of rectangular steel plates under uniaxial
compression, The Structural Engineer, 58B, No. 3 (1980) 45-61.
The Strength of Plates Subjected to Biaxial Forces 353
Jb fe saaeth
He aM plihe ’ 4 i’. : et aed | vx 1 itapayy vy ma ug we
——~ j
yo
; ; oneal ieee
=; lat
| 7 =
| 7
a
( =|
ra.
1 y a ‘YY pormhineg. Vi,
i iA) ‘ ane “Aga aca a ey Mkag y Soe Wy Ww wohilits ;
ut PD, ' Peete “Tie ; vat i i hier anf pie 4)
:
7 x
vent
) P R t -
4a} 7 .
: A od @
18
ABSTRACT
It is shown that simply supported stiffened panels with high aspect ratios can
buckle either locally, in which case the cross-section undergoes considerable
distortion, or globally, in which case the panels undergo an overall dishing.
In shorter panels this dishing is often suppressed with the result that the
panels behave like wide columns because then the effect of the support along
each side on the buckling behaviour is small. By using a so-called buckling
plot obtained by using thefinite strip method it is shown that the change from
local to global buckling can be predicted accurately. Experiments carried
out to verify the application of the finite strip method to stiffened plates
demonstrate that this change can be brought about by changes in the number
of stiffeners or in the depth of the stiffeners and results confirm the
theoretical work.
NOTATION
1. INTRODUCTION
All edges
get
AA q Pe a simply - supported
Oy n stiffeners oe fe
panels are used in the compression flanges of long, narrow box girders, in
the legs of derrick cranes and in box columns. The panels can buckle in at
least two different modes, viz. a local buckling mode and a global buckling
mode (see Fig. 2). Local buckling involves a distortion of the cross-section
with the buckles having a relatively short half-wavelength approximately
equal to the stiffener spacing. For usual cross-sections local buckling
results in the buckling of the plate strips in the well-known chequerboard
pattern with sympathetic twisting of the stiffeners. Global buckling of a
stiffened panel is like that of an orthotropic plate and involves bending of
the stiffeners about their major axis. It is, however, not sufficient to analyse
a panel separately for local and global buckling as these buckling modes
can interact and this leads to a lower buckling stress than both the stresses
for local and global buckling would suggest. This phenomenon is well
known and examples are given by Timoshenko and Gere.’ It is important,
therefore, for a theoretical analysis to analyse the whole cross-section for
buckling in all possible modes and not simply elements of the cross-section,
as for example the plate strips between stiffeners.
The study carried out here is both theoretical and experimental. For the
theoretical analysis the finite strip computer program PLATE, which was
developed by Murray and Thierauf for the preparation of their book of
tables for the design of stiffened steel plates, was used. This program is
similar to the one developed by Stanley.? It is based on a small-deflection
elastic analysis and calculates critical buckling stresses of panels with given
dimensions and boundary conditions. Some results of an experimental
program carried out by Katzer* are used to check and interpret theoretical
results.
A systematic and comprehensive parametric study of stiffened plates is
impossible as there are too many independent parameters. This paper
investigates local and global buckling behaviour and the interaction of
these buckling modes. After briefly describing some of the background of
the theory used (in Section 2) and the experimental program (in Section 3),
358 W. Katzer and N. W. Murray
two studies of the transition from local to global buckling are presented. In
the first, the number of stiffeners n and in the second, the stiffener depth d,
are systematically varied. In Section 6 theoretically and experimentally
obtained buckling stresses are compared and in Section 7 a number of
conclusions which can be drawn from the present study are presented.
2. THEORETICAL ANALYSIS
The stiffened panels described below have been analysed for buckling using
Cheung’s® finite strip method. Details of its specific application to the
analysis of stiffened panels are given by Murray and Thierauf’ and by
Thierauf et al.’
Their method solves for the eigenvalues of a panel with known cross-
section as the length L, of the buckle is increased. The information
contained in their book is the stress at which local buckling occurs, the
corresponding length of the local buckles and the length of panel above
which the panel will buckle in a global manner.
3. EXPERIMENTAL PROGRAMME
were contour maps of the buckled specimen. In every case two nominally
identical specimens were tested in order to check on the reproducibility of
the results.
Experiments on the type of stiffened plates investigated here had earlier
been carried out by Stanley.? Stanley used specimens made from
acrylonitrile butadiene styrene (ABS), a material which has a large elastic
range. His specimens were fabricated by gluing plate strips together and
they were therefore relatively perfect and stress free. In the present testing
programme mild steel plate was used because it was desired to study not
only elastic buckling but also the collapse behaviour of the specimens. All
360 W. Katzer and N. W. Murray
80
60
40
20
20
40
60
80
100
120
140
160
Oy = 171 MPa
+) tension
Stiffeners : SS +1587
+171:0
<a55 +168-°5
[345 +624
+103'5
specimens had initial imperfections and also carried residual stresses owing
to the jigging and welding process. Residual strains were evaluated in some
specimens by measuring the length of the as-welded specimen at certain
grid points of the cross-section, then cutting the specimens into strips in the
longitudinal direction and finally measuring the length of all strips again.
Figure 4 shows the result for a three-stiffener specimen.
In the first study the effect of varying the number of stiffeners on the elastic
buckling behaviour is investigated. One major problem in a theoretical
buckling analysis is that stiffened plates can buckle in many different ways.
From the designer’s point of view, of all possible buckling modes the one
which gives the lowest buckling stress is most important. This buckling
mode is not known beforehand so the theoretical analysis can only be
carried out for an assumed number of buckling half-wavelengths and one
must search systematically for the buckling mode corresponding to the
lowest critical buckling stress. One way of doing this is to assume the
buckling mode is a sine curve with half-wavelength in the longitudinal
direction of L, = L/m with m= 1, 2, 3,.... Another possibility is to vary
the buckling length L, froma very small value in finite steps up to the actual
panel length and calculate the stress o, which is required to maintain the
specimen in equilibrium in the configuration ofa half-sine wave of length
L,. If these values o, are plotted against values of L, a curve or so-called
‘buckling plot’ is obtained. The full line in Fig. 5 shows a typical buckling
Op
(MPa)
TEST C1-23
owe \ SR PER ada EE)
Vomecs ee Bee
' 175 175
400 \ t= 1-97
\
|
|
=950m
TH
5c |
7
al
ish
re)
-_
=
ca |
o |
af
OL, 1000 2000 3000 | (mm)
Fig. 6. Buckling plot for a one-stiffener panel.
500
TEST C2-24
Ween
77777
pen news
777
tes
gee e ‘= a
t= 1:95
300
Og o
200 \8
SS ee ee
panel
of
Length
° Lz fe) } re) 2000 3000 L(
mm)
Fig. 7. Buckling plot for a two-stiffener panel.
Elastic Buckling of Stiffened Steel Plates 363
Ob
(MPa) |
50 |
| Apt
c H TEST C3-27
400
300
Oc. »
Ogata: oo ae
200 | §
| BIN
100 Ll ont
| | ™. EURER
| 2| i> ee
aij — oe
plot for all of the panels investigated here. It shows two distinct minima.
The first minimum at L, is due to local buckling. The corresponding
buckling stress o,, is referred to as the ‘critical stress for local buckling’. The
buckling stress corresponding to L 1s called the ‘critical stress for
global buckling’, o,. The smaller value of o,, and o, is referred to as 04m.
For the panels considered here the second minimum usually occurs at a
length which is greater than the actual panel length L. This buckling plot was
found to be the most satisfactory way of presenting the results of this study.
It is self-evident that ifa panel of length L, can buckle locally at o,, and with
a buckle length equal to L, panels of length 2L,, 3L, and so on will also
buckle at the same critical stress. The dotted lines in Fig. 5 show how this
effect can be included in a buckling plot but since it does not provide
additional information the presence of these additional lobes is tacitly
assumed in the discussion which follows.
Figures 6-8 show the buckling plots for three specimens each with a plate
thickness of approximately 2 mm, a stiffener depth of 40 mm, and one, two
and three stiffeners, respectively. Panels of these cross-sections and a length
of 950 mm have been tested in the laboratory. From the buckling plot of the
one-stiffener specimen one would predict that the panel buckles mainly by
364 W. Katzer and N. W. Murray
(c)
Fig. 9. Contour maps of buckled specimens (¢~ 2mm, d= 40 mm).
Elastic Buckling of Stiffened Steel Plates 365
local buckling, as o,, is significantly smaller than o,. This behaviour was
indeed observed in the laboratory as can be seen from the contour map
shown in Fig. 9(a) and obtained by using the Moire technique. Figures 9(b)
and (c) show the buckling modes observed for the two- and three-stiffener
specimens. As the number of stiffeners is increased above one the b/t ratio
of the plate strips increases and therefore also the resistance against local
buckling. Consequently the first minimum of the buckling plot is shifted
significantly upwards. For the two-stiffener specimen (Fig. 7) o,, is still
slightly smaller than o,. The panel shows an interactive buckling with local
and global modes and this is reflected in the contour map shown in
Fig. 9(b). For the three-stiffener specimen (Fig. 8) o,, is greater than a,. The
contour map (Fig. 9(c)) for this panel shows that there is very little local
buckling and that the main shape indicates global buckling approximating
to that of an overall dishing.
For the sake of comparison the column buckling stress of the specimens,
the Euler stress o,, is also included in all three buckling plots. o, was
calculated using the simple formula o, = n*E/(L/r)* of a pin-ended Euler
strut, thus neglecting the influence of the two lateral supports. It can be seen
that the difference between the Euler stress o, at L and the critical stress for
global buckling is greatest for the one-stiffener specimen. For the three-
stiffener specimen the Euler curve is an asymptote for the part of the
buckling plot which corresponds to global buckling. However, with
increasing L, the difference between o, and a, grows. For large L, the Euler
stress approaches zero whereas for those specimens which had simply
supported sides the global buckling stress reaches a minimum and increases
again with increasing L,. This latter behaviour is similar to that of a very
long simply supported orthotropic plate for it is well known!’ that such a
plate develops long buckles at its critical stress.
Local ee a
10° LGlobal buckling 7
panel length
= 950 mm
predicted to occur first for all three panels, i.e. m = 1. The observed contour
maps from the Moiré technique of the relevant specimens after buckling are
shown in Fig. 13(a—c). The one-stiffener specimens showed severe local
buckling deformations although o, was about 25 % smaller than o,,. The
two- and three-stiffener specimens showed the global buckling behaviour
predicted by theory.
For the one-stiffener specimen and d>20mm the buckling mode
changes to local buckling (Fig. 10, m = 5, thus L, = 190 mm). Up to about
d= 60 mm the curve is fairly horizontal indicating that in this region mainly
local buckling of the plate strips occurs. With further increase in stiffener
depth o,,, decreases and thus the cross-section becomes weaker against
local buckling. The corresponding buckling mode involves both plate and
stiffener deformations. The two- and three-stiffener panels have a lower b/t
ratio and hence a greater resistance against local plate buckling. The graphs
in Figs. 11 and 12 show that the change from the global buckling mode to
the local buckling mode consequently takes place at a greater value of o,,,.
Elastic Buckling of Stiffened Steel Plates 367
One
5x10° NOTE : Logarithmic Scale and False Origin for Ordinate
10°
5x107
12 !126 ‘112
Local Saati
Global eas ode Se
buckling ote aad nae pear
10°
5x10
m=5
+ t=1-5
dt AST sipanel length
104 7707 7777 = 950mm.
eae
80 95 95 80
All panels had initial imperfections and therefore did not buckle suddenly
at their critical load. Table 1 shows a comparison between the theoretically
obtained critical buckling loads and the load range at which local plate
buckles formed during testing. Measurements of residual strain had
shown‘ that the plate strips carry compressive residual stresses o, over most
of their width. o, was found to be approximately 30, 50 and 70 MPa for the
Elastic Buckling of Stiffened Steel Plates 369
(b)
Table 1. Comparison of Observed and Predicted Buckling Loads for Local Plate
Buckling.
a f ; =
AD 13 68 35 55-10
Blagps (a8 : a
GIL din ; ma
Head no : a
1233 3 22 55-15
r Pe = (6.,A, te P,).
4 Re = (Soy 7 0,)A, as Pe
© The first value in this column indicates the approximate load at which first local
plate buckles were observed anywhere in the panel and the second value indicates
the approximate load at which a regular pattern of local plate buckles had formed
over the whole panel.
7. CONCLUSIONS
(i) The test series described here is confined to stiffened panels with
high aspect ratio, with simple rectangular stiffeners and with simple
supports along each longitudinal edge.
372 W. Katzer and N. W. Murray
(il) One aim of the programme was to check that the elastic buckling
stresses and patterns predicted by Murray and Thierauf * who used
the finite strip method could be confirmed by experiments.
Although the testing programme was limited in its scope it can be
concluded that experimental and theoretical results so far obtained
agree quite well. No discrepancies have yet been observed.
(iii) Buckling plots, typically shown in Figs. 6, 7 and 8, show that as the
number of stiffeners increases, the buckling mode may change from
local to global with interactive buckling in the transition region.
Contour maps of the buckling patterns obtained by experiment
(Fig. 14) confirm this.
(iv) Graphs of the theoretical buckling stress against depth of stiffener,
typically shown in Figs. 10, 11 and 12, show that as d increases
from zero, where the panels behave as an isotropic flat plate, up to
an intermediate value the panels behave as orthotropic plates.
Beyond a certain intermediate value the buckling stress decreases
and the local buckling mode is largely governed by torsional effects
in the deep narrow stiffeners. Contour maps of the buckling
patterns obtained by experiment (Figs. 13(b) and 14(b)) confirm
this change from global to local buckling as d increases.
(v) The global buckling is relatively easily suppressed by use of
transverse stiffeners which do not significantly increase the local
buckling stress.
Acknowledgements
The financial help of the Deutsche Forchungsgemeinschaft, Deutscher
Akademischer Austauschdienst and Monash University which have
enabled the authors to pursue this research is gratefully acknowledged.
REFERENCES
1. TIMOSHENKO, P.S. and Gere, J. M., Theory of elastic stability, 2nd edn., Tokyo,
McGraw-Hill and Kogakusha, 1961.
2. Murray, N. W. and THIERAUF, G., Tables for the design and analysis ofstiffened
steel plates, Braunschweig, Wiesbaden, Vieweg, 1981.
3. STANLEY, C.R., ‘The critical buckling stress of plates stiffened with longitudinal
stiffeners’, M.Eng.Sc. Thesis, Dept. of Civil Engg., Adelaide University, 1973.
4. KaTzerR, W., ‘The collapse behaviour of stiffened plate panels of high aspect
ratio under uniaxial compression’, M.Eng.Sci. Thesis, Monash University,
August 1982.
Elastic Buckling of Stiffened Steel Plates 373
OPRF ’
19
A. C. WALKER
University of Surrey, Guildford, UK
ABSTRACT
NOTATION
b Piate width.
Bs Effective width.
c Coefficient.
E Young’s modulus.
Hy, uy Coeihcients.
DP. Critical load.
t Plate thickness.
w Deflection.
Wo Initial deflection.
€ Plate end shortening.
Eo Yield strain.
v Poisson’s ratio.
Average applied stress.
375
376 A. C. Walker
6. Critical stress.
Om Maximum stress.
Fo Yield stress.
oe Ultimate value of average applied stress.
1. INTRODUCTION
Number of ay Wa
papers 40 y
published Vy)
: We
30 y g
10 | Z g
AY AY
SERRE
REE ES RRS
Fig. 1. Papers published on plate buckling.
Plate Buckling Research 377
this is a fairly recent phenomenon and reflects the tremendous growth in the
power of computers and numerical analysis. However, this is not the case;
the great majority of the papers on plate buckling during the past ninety
years has been concerned with theoretical developments. Notwithstanding
the value of these contributions, I feel that the really important
developments in understanding the phenomenon of plate buckling and
developing information for designers have been made when theory and
experiment are brought closely together. In this paper, therefore, the
experimental aspects of research are emphasised.
Indeed it has been one of the great merits of the work carried out in the
UK that the development of analysis and testing facilities has proceeded
hand-in-hand. Ona personal note, I have a feeling of gratitude to Professor
Harvey and his colleagues who impressed this attitude so strongly upon me
when I worked in their department at Strathclyde, or the Royal College of
Science and Technology as it then was.
2. 1850-1930
iy
wsq
his railway then under construction. He had decided to use wrought iron
which could be produced in small sheets and rivetted together to form
larger plate structures. He designed a novel form of structure in which the
train would pass through the centre of the box which itself would rest on
piers and be supported by chains. Figure 2 shows the general arrangement
and Fig. 3 shows the cross-section.
Realising that his structure was in fact quite unique in concept and in
advance of contemporary engineering practice, Stephenson arranged with
Fairbairn to carry out a number of tests to determine the strength of box
beams having various cross-sectional forms. It is important to remember
Plate Buckling Research 379
fig
3. 1930-1950
it is seen that von Karman’s value for the coefficient c would be c= 1.9.
However, as Fig. 5 shows, the value of c is indeed variable and Sechler
0.067 Dural
0:012 Steel
0-050 Aluminum
0-016 Steel
0.052 Steel
ex 0:016
—-Q20 Brass
Loading bar
specimen
Deflection reading
apparatus
0
0000000
qajiooo0oo00e0
Testing machine
Fig. 6. Testing arrangement.
estimated from his test results that the constant should have a value of
almost 2 for very thick plates but reduce significantly for thin plates. Thus
the simple concept of effective width provided a good basis for the
incorporation of test results; however, it is worthwhile to remember the
basis of these tests. They were carried out on single plates loaded axially by
flat plattens. Sometimes the plate edges were rounded to ensure that the
load was introduced in the mid-plane. The longitudinal edges were
supported in vee grooves or between knife-edges and nominally were free to
pull-in as the plate deflected (see Fig. 6, ref. 3, p. 424). It was generally
assumed that these supports gave conditions which were similar to the
theoretical simple support. However, it is by no means certain that simple
support was indeed attained and the interpretation of the load—deflection
curves using Southwell’s method will not be reliable in assessing the efficacy
Plate Buckling Research 383
4. 1950-1965
b,=1.9t
JicLt-048(5) JerI
/—]1-—0.475{—]
oe b
/—
Gin
(4)
in which g,, is the maximum stress on the plate element. In other words,
Winter showed that the effective width defined by eqn (4) is applicable
throughout the elastic post-buckled regime of the plate and not just at the
ultimate load. This formula acted as one of the bases for design for cold
formed sections in the USA until 1968 at which time it was considered that
sufficient experience had been accumulated to change it to
pai [E[i—oais(t)
Os b
[F
Oe.
(5)
or, in terms of the collapse stress,
vb ify
b= =Om
Fig. 7. Correlation between the effective design width formula and test data.
This is shown in Fig. 7, where it is seen that there is fair agreement between
the test results of Winter and the design equation.
In a very interesting recent study Roorda® carried out a critical analysis
of a number of effective width formulae relevant to cold formed section
plate elements and concluded that Winter’s form of relationship is the best
of those proposed in the available literature, but commented that some
revision of the coefficients would result in even better agreement with the
published experimental results.
In Britain about this period there was also active interest in the buckling
behaviour of cold formed sections. Research into local buckling was being
actively pursued at Bristol University and in the Royal College of Science
and Technology in Glasgow. Chilver® proposed a formula for plates as
1/3
ow a6 @ (7)
a) a)
and Kenedi et a/.!° showed that in comparison with the result of tests on
some 200 mild steel sections of different shapes the limits of the scatter of
these tests are shown in Fig. 8; eqn (7) defines the average maximum stress
Plate Buckling Research 385
Equation ( 7)
Equation (8) was also shown to be adequate for mild steel single plates with
different types of longitudinal edge support, which suggested that this
expression could be used for a wide range of thin-walled mild steel sections.
Indeed, eqn (8) formed the basis of the calculations of buckling strength
of cold formed sections and was contained in Addendum No. | to BS 449.
The proposed approach was to calculate the strength of a section as the sum
of the individual strengths of the component plates, assuming these to be
simply supported or free on the longitudinal edges. In a later paper‘! it was
shown that the approach could be in error in the range 0.25 < (a,/a9) < 1.
However, this factor was recognised in a subsequent revision of the
Addendum; that design curve is shown in Fig. 9.
It can be seen that up to that time the mechanics of behaviour of single
plate behaviour was fairly well explored’? '* and that design formulations
were based on non-linear analysis and the empirical interpretations of test
results. There was also a growing interest in the development of more
precise solutions of von Karman’s non-linear equations. A particularly
significant paper was published in 1959 by Yamaki'® in which he used a
multi-term deflection series solution and evaluated the coefficient from the
386 A. C. Walker
0-4
0:2
1 2 3 4 5 6 i
VO9/ Oc
Fig. 9. Comparison of generalised collapse with test results for simply supported
plates in various metals.
Bearing holder
Loading beam
J 1965-1983
This period has been one of tremendous activity in the research of plate
buckling, as evidenced by the large number of papers published. So far as
steel plates are concerned two factors have driven this work: one is the
phenomenal increase in the power of computers and their ready
availability to researchers; the other factor has been the unfortunate
accident to a number of box girder bridges which caused engineers to
review their state of knowledge with regard to the design of these complex
structures.
Taking the analytical development first, there can be no doubt that the
advent of high-speed digital computers has revolutionised our approach to
plate analysis. Although the mechanics of plate behaviour was well
established, the design approaches were essentially empirical. The
numerical analysis developed during the past fifteen years have permitted
the design of plates having a wide variety of boundary and loading
conditions to be put on to a firmer theoretical basis. Plasticity and non-
linear geometrical deformations can both be incorporated into the
numerical analyses and from a knowledge ofthe initial deflections and state
of stress it is possible now to predict, with acceptable accuracy, the
maximum load. The problems with the communication of this information
from researcher to designer is something I will return to later in the paper.
Two schools of research in Britain have been pre-eminent in the
development of computer methods although, ofcourse, individual workers
elsewhere have made significant contributions. The work at Cambridge
University directed by John Dwight led to the first publications'®:?° in
which the important factors of non-linear deformations and plasticity were
combined effectively in a single numerical analysis scheme. Later work at
Cambridge, particularly by Moxham and Little, developed these ideas to
388 A. C. Walker
8 e
(b/t)
¥ (Oo/E) = 2-074
E /Eo
Fig. 12. Load end shortening curves for constrained plates (wot = 0.188).
390 A. C. Walker
UMOYS
3a
‘asplruquiey
|
AjUO‘ssaBuUI) wo2q
NOILOSS
V-V¥
3oq 12>
dwoj> poDoO7uvoly>024y
+++4+44
42 4 ae
104
ysoddnsvawisads NW
Saket
U2dIWDId
$)0H
uawisads
I
"|
teeh]
a6pzZ
Bulpyong
3
=a!
i /
ie
—/ yy
uF /
/ / /
I. |
/
va saro0s6-A
/
++
++
ae
dbo
Joj ayeid
JUsWUIDsUBLIe
ye
x0qJD26 Ze Ul
ssajjou
buls02q disys
Se
¥4JO}OWM2J2S Le
ce Wier
ajpaau
6
ul} Buispaq
ysNJYy
1adD) uaqD
WOJ} DWIY ‘saping VE
% jmaoy Id
‘26pam
AIIP
ANU Ulz/y-|— VA
{UHL
7
ISO,
oe
BAIIQ ca pe
UIDUD
TWMMLUIIIE
—
“BIA “ET
tee
eae
ttt
Aj0\010qQ07
8
‘q poxig SUd|]OU
ysay DOIG
- asoo7
uld Pd}OIO}
Bulsnlpy ul yoisul Jo6ulj)
MdIIS
GudjD BxOKA SPUB}YXD
O}
YUOU
Y xD0Q
JO Jd6ulj)
ayo, payoddns
AQ DBIpzou
42}j}ou Buiy22foudY6nosy, 4s6uij
BulAOWwBxOK0} 442] saysnd
uid
dn adojs
0} diuh UdBWIDDIS
uo4 podwo)}>d
BOpD BSN OM) suid
PuDdD DAOW see UDI}OU
OS DxOA
Ssy>d0U j1yUN yyOQ Suid Giu6
so6ui4
Plate Buckling Research
yabuly pduuid
©}
Buljsoddns
woaq
B14 ‘pluoo—€]
391
392 A. C. Walker
W
Se
a3
<
1 ; Southwell plot
0-10 045
Ah aan seen tye target
(a)
all
Dc = 2850lb ( 12-7 kN)
= 10
E
=
"o
PE
=
6. CONCLUSIONS
REFERENCES
34, VepA, Y. et al., J. Soc. of Naval Architects of Japan, 140 (1976) 20S.
35% KRAGERUP, J., Report No. 143, Series R, Technical University of Denmark,
1982.
36. BRADFIELD, C. D., J. Constr. Steel Research, 1 (1980) 27.
3). Becker, H. et al., Report No. SSC-217, Ship Structure Committee,
Washington, D.C., 1970.
38. BECKER, H. and Cotao, A., Report No. SSC-267, Ship Structure Committee,
Washington, D.C., 1977.
20
ABSTRACT
NOTATION
a Length of plate.
b Width of plate.
b. Effective width.
Cy Om/Fo.2 = 5/5.
é Width of residual tension block at one edge.
399
400 D. S. Mofflin and J. B. Dwight
E Modulus of elasticity.
n Ramberg—Osgood formula index.
t Thickness of plate.
Wo Initial imperfection amplitude.
u, v, w, 8 Displacement parameters.
V Volume of plate.
eye Coordinate axes.
B ./ 69 ./Oc, = (b/t)/1.92 ./o9.,/E for aluminium (v = 0.33).
E Applied strain.
£2 Strain at the 0.2% proof stress.
o Stress or average applied stress.
Bo 0.2% proof stress.
Centre zone stress (simple model).
Elastic critical stress.
Modified elastic critical stress (simple model).
Edge zone stress (simple model).
Limiting value of average applied stress.
Residual compressive stress.
Yield stress of elastic-perfectly plastic materials.
1. INTRODUCTION
Simply Supported
Free to Pull-in
(d)
Fig. 2. Factors affecting plate strength. (a) Initial out-of-flatness; (b) residual
stress; (c) HAZ softening; (d) stress-strain curve.
Buckling of Aluminium Plates in Compression 403
3. STRESS-STRAIN CURVE
oO @ \"
é=—, + 0.002 (=) (1)
90.2
which produces a curve which automatically has the correct strain at 0.2%
proof stress level. The knee-factor nis a measure of the shape of the curve, a
high n indicating a relatively abupt knee and a low na more rounded one
(Fig. 3).
0002 €
Fig. 3. Ramberg—Osgood stress-strain curves.
$80
9/STSOE 100°0
1 9S 0 Ll
9/ST0€ 0$00°0 9IL0
9/SMSOE 100°0
~
Stl
Sbl'0
9/SMA0E $00°0
7
Is]
Sbl'0
IS]
6/VOSOE O-€80S 16 100°0
0 L8 9°°0 0 7$6
O€
Ol
6/VNA0E 1$00°0 18 6880
6/VMA0E 09000 68 1€60
b/9NSOE AL-C809 COC 000'0
$ OVC vv30
0€
NE
b/9NA0E $000
1 VEC c08'0
c0'l
b/9TSOE 1000
~ 87 6LL'0
b/9TEOE 0900°0 LNG 0 cb
b/9MSOE 100°0
0 807 0 cIL
p/OMA0E 700°0
9 SOC 0 COL
1/SASOr W-€80S 000°0
v 0 0€L
OV
91
aor 1/SN
C81
700°0
8
cel
vOL0
£01
1/STd0r
8c
$00°0
$ 1790
1/SM0r 65000 S90
6/VASOr OV Ol 000°0
8 0 vl 0 tv8
16
ag Lil LL
6/VW180r L00°0
¢ 0820
Buckling of Aluminium Plates in Compression
O-€806
0b 6/VM 0900°0 0€L'0
We 99
SC
s/9Nd0r 05000 [p90
Gs
€LS'0
Ls/9NS0r 0 c00
€ 9150
L8I a IST
405
(panuijuod)
406 D. S. Mofflin and J. B. Dwight
aspo
ue) 6690ps0
brs 8eL'0
€19°0 ges'0 plr'0€1s‘0 p8r'0
€7S'0
6SS'0 8LE0oss'0
Srs'0 ple0
Lee‘0
ISI]
0
sooUaP
g
°[00'0
] = ee LY = hes a = Su ce = 67 «1 = — = M
~
g/°M
g/on 710008 0900°0
€s000
ZS00°0 8 LS00°0
1100°0 0500°09500°0
€600°0 6 0S00'0
6900'°0
700°0 500°0 S 700'0
S sajOUAP
~3=
#91 OL 8 bl LI §
aiqeL
‘pluoo—{
y ce
“PlaM
(,wu/A’)
ISIOASUB
700 781 16 767 Z0Z STI 807
SO]OUDP
Aon
W-£805 O-€80¢AL-7809 W-€805 W-£805 W-€805
,
“UW
jDUIWO
1/q 0S 0S (Ns 09 OL i) Q
=
NN SSoUYIY
°0¢'()
“©
‘fou
=
Oo
300
ae
6082-TF
200
508 3~m
100
ee 5083-0
0.005 €
Fig. 4. Comparison between experimental stress-strain curves and Ramberg—
Osgood approximation.
4. TEST PROGRAMME
Referring to Fig. 13a of Chapter 19, the special features are the motor-
driven wedge-straining device and the array of‘fingers’ along the unloaded
edges of the specimen. The role of the fingers, which owed their original
development to Moxham, was to provide a determinate edge condition,
namely simple support with free pull-in, without taking any appreciable
part of the longitudinal load.'°
In each test straining continued without pause, readings of load and
specimen shortening being fed into an X—Y plotter. The duration of a
typical test was about one hour. Load was measured using a conventional
load-cell, and shortening by means of two displacement transducers
located either side of the specimen.
4.2. Specimens
The following list summarises the various specimens tested (see also
Table 1):
The 6082 was all in the fully-heat-treated TF condition. The bulk of the
5083 plates were tested in the as-received M condition (i.e. hot rolled). A
limited amount of the 5083, however, was annealed in the laboratory before
testing (O condition).
All the specimens were given a controlled ‘bump’ (out-of-flatness
imperfection w,) using a hydraulic jack. Springback made it difficult to
control this operation precisely, and Table | lists actual w,-values achieved.
A controlled fusion weld was made down the longitudinal edges of about
half the plates, using a TIG torch mounted on a tractor. The resulting
level of compressive residual stress o, was monitored by taking strain
readings in the middle of each plate with a Weldscan extensometer,!° the
corresponding value of c being then calculated by statics. The results
Buckling of Aluminium Plates in Compression 409
(Table 1) show that the heat-treated alloy 6082 tends to suffer less residual
stress than the 5083, as was predicted by Wong.'’ The extent of the HAZ
softening adjacent to the weld was determined by means of a hardness
survey on two of the 6082 specimens, and was found to correspond roughly
to the width c of the residual tension zone.
A transverse weld bead was laid at midlength on ten of the specimens by
the MIG process, the bead being ground off before testing.
4.3. Results
Table 1 lists failure loads for all the specimens in the form of C,(=a,,/do >).
Typical load-shortening curves appear in Fig. 9.
individual webs or outstands with various edge conditions. It has also been
successfully extended to cover assemblages of plates.
The method uses an energy solution procedure, which provides stable
convergence characteristics and allows unloading to be followed. Here we
describe its application to the specific case of a web simply supported on all
four edges, with the unloaded edges free to pull in and the loaded edges
constrained to stay straight. The plate is assumed to buckle in one half-
wave.
In this application the method has similarity to the earlier analyses of
Moxham! and Little.* The main difference is the splitting up of each
quadrant of the plate into strips, so providing greater flexibility in the
description of the deflected form.
(a)
Fig. 6. Finite-strip model. (a) Typical strip division; (b) definition of parameters;
(c) monitoring points within each strip.
= |[{o} + 7{Ao}]
{Ae} dV
A@g* = change in potential energy of the loading system.
(5) By minimising the energy with respect to the displacement
parameters, using a quasi- Newton minimisation procedure, obtain
the equilibrium deflected shape.
(6) Compute the corresponding stresses and hence the applied load.
This provides one point on the required curve.
out-of-flatness for welded plates will therefore be very slightly greater than
the values quoted.
HAZ softening in welded plates is dealt with, when necessary, by
assigning to the two outer strips (representing the zones of residual tension)
values of ¢) , and appropriate to the weakened material. This implies that
the extent of the HAZ is effectively the same as that of the tension zones, a
convenient assumption to make until more is known about HAZs.
6. SIMPLE ANALYSIS
SSS ee
fee ey
Fig. 7. Simple model. Assumed stress distribution.
Fig. 8. Simple model. o—¢ curves for plate with high b/t. (1) Edge zone; (2) centre
zone; (3) average.
414 D. S. Mofflin and J. B. Dwight
where E£, is the tangent modulus and E, the conventional reduced modulus.
In order to determine the falling branch of the curve, involving the
calculation of the reduced critical stress a/,, it is assumed that material is
fully effective until the strain reaches oo.,/E and ineffective thereafter.
The effect of HAZ weakening adjacent to edge welds is readily allowed
for by splitting each edge zone into two parts. The outer part, representing
the HAZ material, is given a suitably reduced stress-strain curve.
7. DISCUSSION
(1) There is a clear difference between the curves for the two alloys,
with the 5082 tending to keep on rising. This is clearly a function of
the o—é curve (Fig. 4).
(ii) For the 6082 the experimental curves exhibit a much steeper
unloading path than the theoretical ones. This is not a defect of the
theory. It is merely that the test specimens were all of the length 4),
whereas the analyses considered a plate only one buckle long
(a =0.8b or b). For a long plate loaded in isolation only one buckle
is critical. The other buckles unload into it after peak load is passed,
producing a more steeply descending curve when plotted in terms
of average strain on the overall length.'?
Buckling of Aluminium Plates in Compression 417
Fig. 10. Definition of ultimate strength. (1) o—-e curve of material; (2) high-n
material; (3) low-n material.
be more important, the welded strengths being significantly lower than the
non-welded.
The absolute strengths of the 5083 are of course well below those for the
6082, because of its lower proof stress. When, however, the comparison is
made non-dimensionally (in terms of C,) the difference between the two,
although significant, is much reduced. In the low to medium range of f the
tests have C, for the 5083 lying some 10-15 % below that for 6082, although
the theoretical results show a difference of only a few percent.
It should here be noted that, for a given n value, the supposedly non-
dimensional plot of C,, against f is not quite independent of a,,,tending to
shift up and down slightly with og,,.All the theoretical curves presented in
Fig. 9 are based on o, ,= 200 N/mm?. Had the correct value been used
(69.. ~ 290 N/mm? for the 6082, and a, ,~ 180 for much of the 5083) the
effect of C, would have been small, pushing it up very slightly for the 6082
and down for the 5083.
WSSSSEEOS es
Sa cre-s3c---— S
3S
& 1S)i)
nH
=2 s
(2)
(P)
“POP[oM
:PePiem
“W-€80¢
“41-7809
:Ppoplom-uoU
‘pluoo—TT “By
sve oz set
D. S. Mofflin and J. B. Dwight
O = 3/5
100°O = a/¥
aT3M ON ‘SCO'O = a/y ~~}
"INLNAW TYR IXS A
aQTaM ON “100°0O = Q@/9 °
2LINIA
XYOMILL STINIS
dINLS
-——-—
2 ar ee
Oo
420
——-—~ aLINIG d1luLs
FIdWIS ANOGHL ie) z°o= 0oz
‘ =u S
Ss‘
“By
| ‘pluod—]
421
422
e
“Se. °
cae °
Se
= ie e q/¥= $00°O
ee ° 3/2= €
Sse HLIA ONINSLIOS
ASN e
Sa
D. S. Mofflin and J. B. Dwight
“Bi Piuoo—T]
vt
"
=)
Vv
°
"
Buckling of Aluminium Plates in Compression
5083 plates, but does have an appreciable effect for 6082—as one would
expect, due to the HAZ softening which is produced right across the plate.
The reduction on a,, was of the order of 30 % for 6082 plates containing no
other weld, and about 20 °% for plates which already had longitudinal edge
welds. These decreases seem remarkably small, when one considers that
0... for the HAZ material is about half that for the parent metal.
It is seen that for the large imperfection the three curves are virtually
indistinguishable. At the smaller imperfection there is only a slight
difference between the n = «© and n =25 curves, but an appreciable drop
tonn = 3:
The important question is: ‘can a unified C,—f design curve be used for
all materials? The answer appears to be a limited yes, with some exception
made for materials having a very rounded stress-strain curve. It is probable
that cold-reduced light-gauge steel has a g—-e curve comparable to that of
6082-TF in terms of n value, so that a unified plate strength curve could
cover light-gauge as well as hot finished heavy gauge steel. But this
conclusion is tentative at this stage bearing in mind the considerably greater
difference between 6082 and 5083 which has been observed experimentally.
8. CONCLUSIONS
(1) The finite strip theoretical model gives reasonable agreement with
tests.
(2) The simple theory gives remarkably good results, bearing in mind
its relative crudity.
Buckling of Aluminium Plates in Compression 425
(3) The tests showed an appreciably lower strength for 5083 plates in
the low to medium range of b/t, as compared with 6082—when
plotted non-dimensionally (C, vs. B). The theories predicted a
much smaller difference.
(4) Plates with longitudinal edge welds have an appreciably reduced
strength due to residual stress and (for the 6082) HAZ effects.
(5) Transverse welds had negligible effect on the strength of the 5083
plates, but reduced the strength of the 6082 plates by some 20-30 %
due to HAZ softening.
(6) Theoretical comparisons suggest that a common non-dimensional
design curve (C,, vs. £3) might possibly be acceptable for unwelded
heavy gauge steel, light gauge steel and heat treated aluminium—
provided the levels of out-of-flatness are comparable.
Acknowledgements
The work formed part of a Ministry of Defence programme of research on
structural aluminium at Cambridge University Engineering Department,
funded jointly by the Admiralty Marine Technology Establishment and by
the Military Vehicles and Engineering Establishment. A major role was
played by Mr M. J. Free in the Structures Laboratory at Cambridge.
APPENDIX
The displacement functions used in the finite strip method for describing
the deflected shape within each strip were (see Fig. 6):
REFERENCES
COMPOSITE STRUCTURES—2
edited by I. H. Marshall
6 x 9”, xvi + 576 pp., 309 illus.