Design Guide 26 Design of Blast Resistant Structures
Design Guide 26 Design of Blast Resistant Structures
Design Guide 26 Design of Blast Resistant Structures
Design of Blast
Resistant Structures
Exterior Gauges
Pressure
Impulse
Impulse (psi-msec)
Pressure (psi)
Time (msec)
26Steel Design Guide
Design of Blast
Resistant Structures
RAMON GILSANZ, Lead Author
Gilsanz Murray Steficek LLP
New York, New York & Los Angeles, California
RONALD HAMBURGER
Simpson Gumpertz & Heger, Inc.
San Francisco, California
DARRELL BARKER
ABS Consulting
San Antonio, Texas
JOSEPH L. SMITH
Applied Research Associates, Inc.
Vicksburg, Mississippi
AHMAD RAHIMIAN
WSP Cantor Seinuk
New York, New York
by
The information presented in this publication has been prepared in accordance with recognized
engineering principles. While it is believed to be accurate, this information should not be used
or relied upon for any specific application without competent professional examination and
verification of its accuracy, suitability and applicability by a licensed professional
engineer, designer or architect. The publication of the material contained herein is not
intended as a representation or warranty on the part of the American Institute of Steel
Construction or of any other person named herein, that this information is suitable for any general
or particular use or of freedom from infringement of any patent or patents. Anyone making use of
this information assumes all liability arising from such use.
Caution must be exercised when relying upon other specifications and codes developed by other
bodies and incorporated by reference herein since such material may be modified or amended
from time to time subsequent to the printing of this edition. The Institute bears no responsibility
for such material other than to refer to it and incorporate it by reference at the time of the initial
publication of this edition.
Ronald Hamburger, P.E., S.E., is a senior principal at Simpson Gumpertz & Heger, Inc. He chairs the
ASCE 7 Committee and the AISC Connection Prequalification Review Panel, and is a member of AISC Task
Committee 9 on Seismic Design. He also chaired the NCSEA Ad Hoc Joint Industry Advisory Committee
that developed the structural integrity provisions contained in the International Building Code.
Darrell Barker, P.E., is vice president for Extreme Loads and Structural Risk at ABS Consulting. He is a
member of ASCE, ACI, ASIS, ASME and PGCI. He contributed to the Handbook for Blast Resistant Design
of Buildings and to ASCE/SEI 59-11, Blast Protection of Buildings, and CSA 850-12, Design and Assess-
ment of Buildings Subjected to Blast Loads.
Joseph L. Smith, PSP, is a director and senior vice president of Applied Research Associates, Inc. He is an
active member of ASCE, SAME, SARMA, PGCI and the ASIS International participating on several profes-
sional committees. He is member of the Department of Homeland Security Explosives Standards Working
Group and past Co-Chair of the Explosives Hardening and Mitigation Subgroup, and also contributed to
ASCE/SEI 59-11, Blast Protection of Buildings, and the Handbook for Blast Resistant Design of Buildings.
Ahmad Rahimian, Ph.D., P.E., S.E., is chief executive of WSP Cantor Seinuk in New York. He is an active
member of ASCE, AISC and chairs ACI committee 375. He was a member of the Department of Buildings
panel that wrote the New York City structural integrity code recommendations.
Acknowledgments
The authors appreciate the input and ongoing guidance during the development of this design guide provided
by the following:
David Holgado, ABS Consulting Anders Carlson, Gilsanz Murray Steficek LLP
James T. Brokaw, Applied Research Associates, Inc. David Chlebus, Gilsanz Murray Steficek LLP
Larry M. Bryant, Applied Research Associates, Inc. Thibaut Dehove, Gilsanz Murray Steficek LLP
Kenneth W. Herrle, Applied Research Associates, Inc. Karim Ezzeldin, Gilsanz Murray Steficek LLP
Charles C. Ellison, Applied Research Associates, Inc. Wenjun Guo, Gilsanz Murray Steficek LLP
J. Mikhael Erekson, Applied Research Associates, Inc. Eugene Kim, Gilsanz Murray Steficek LLP
Karl Rubenacker, Gilsanz Murray Steficek LLP Arturo Montalva, Gilsanz Murray Steficek LLP
Brett Benowitz, Gilsanz Murray Steficek LLP Andrew Sparn, Gilsanz Murray Steficek LLP
Preface
This Design Guide provides guidance for the design of blast resistant structures and progressive collapse
mitigation. Background information and some basic principles are reviewed, as well as the presentation of
design examples. The goal of this Design Guide is to provide enough information for a structural engineer to
effectively interact with a security or blast consultant.
i
ii
TABLE OF CONTENTS
CHAPTER 1 INTRODUCTION. . . . . . . . . . . . . . . .1 CHAPTER 4 STRUCTURAL RESPONSE
TO BLAST LOADS. . . . . . . . . . . . . . . . . . . . . . . . 35
1.1 HISTORY OF INCIDENTS . . . . . . . . . . . . . . .1
1.1.1 Blast Incidents . . . . . . . . . . . . . . . . . . .1 4.1 REPRESENTATION OF BLAST LOADING. . .36
1.1.2 Progressive Collapse Incidents. . . . . . . . .2 4.2 SINGLE DEGREE OF FREEDOM SYSTEMS . .36
1.2 CHARACTERISTICS OF 4.3 BLAST RESPONSE OF ELASTIC SINGLE
BLAST EFFECTS . . . . . . . . . . . . . . . . . . . . .3 DEGREE OF FREEDOM SYSTEMS. . . . . . . .37
1.3 BLAST EFFECTS VERSUS 4.3.1 Time-History Analysis . . . . . . . . . . . . . 37
SEISMIC EFFECTS . . . . . . . . . . . . . . . . . . . .4 4.3.2 Graphical Solution . . . . . . . . . . . . . . .38
Example 4.1—Determination of the Peak
CHAPTER 2 BLAST LOADS . . . . . . . . . . . . . . . . .7 Dynamic Force and Displacement . . . . . . . . . . 39
4.3.3 Energy Solution . . . . . . . . . . . . . . . . . 40
2.1 EXPLOSION PARAMETERS . . . . . . . . . . . . .7
4.4 ANALYSIS OF NONLINEAR SINGLE
2.2 EXPLOSIVE THREAT SCENARIOS. . . . . . . . .8
DEGREE OF FREEDOM RESPONSE . . . . . . .40
2.3 BLAST PHENOMENA . . . . . . . . . . . . . . . . . .8
4.4.1 Time-History Methods . . . . . . . . . . . . . 41
2.3.1 Key Parameters. . . . . . . . . . . . . . . . . . .8
4.4.2 Graphical Solutions . . . . . . . . . . . . . . . 41
2.4 BLAST LOAD PREDICTION . . . . . . . . . . . . 10
Example 4.2—Determination of
2.4.1 Empirical Relationships . . . . . . . . . . . . 10
Ductility Demand . . . . . . . . . . . . . . . . . . . . . 41
2.4.2 External Loads . . . . . . . . . . . . . . . . . . 10
4.4.3 Energy Methods . . . . . . . . . . . . . . . . . 42
2.4.3 Internal Loads . . . . . . . . . . . . . . . . . . 12
4.5 MULTIPLE DEGREE OF FREEDOM
2.4.4 Analytical Methods . . . . . . . . . . . . . . . 12
STRUCTURES . . . . . . . . . . . . . . . . . . . . . .42
2.5 LOADS ON STRUCTURES . . . . . . . . . . . . . . 12
4.6 SOFTWARE . . . . . . . . . . . . . . . . . . . . . . . .45
2.5.1 Equivalent Load Shapes . . . . . . . . . . . . 13
2.5.2 Drag Loads . . . . . . . . . . . . . . . . . . . . 14
CHAPTER 5 BLAST RESISTANT DESIGN OF
2.5.3 Nonreflected Surface Loads . . . . . . . . . 14
STRUCTURAL SYSTEMS . . . . . . . . . . . . . . . . . .47
2.5.4 Shielding and Reflection . . . . . . . . . . . 14
2.5.5 Net Lateral Loads . . . . . . . . . . . . . . . . 14 5.1 ENERGY METHOD . . . . . . . . . . . . . . . . . . . 47
2.5.6 Negative Phase . . . . . . . . . . . . . . . . . . 15 5.2 SIMPLIFICATIONS BASED ON DYNAMIC
2.5.7 Interior Loads Due to Leakage. . . . . . . . 15 PARAMETERS . . . . . . . . . . . . . . . . . . . . . .50
2.6 RESOURCES . . . . . . . . . . . . . . . . . . . . . . . 15 5.3 DESIGN EXAMPLES . . . . . . . . . . . . . . . . . . 51
2.7 DESIGN EXAMPLE. . . . . . . . . . . . . . . . . . . 16 Example 5.1—Blast Resistance
Example 2.1—Preliminary Evaluation of Blast of a One-Story Building . . . . . . . . . . . . . . . . . 51
Resistance of a One-Story Building . . . . . . . . . 16 Example 5.2—Blast Resistance
of a Three-Story Building . . . . . . . . . . . . . . . . 59
CHAPTER 3 DESIGN CRITERIA
FOR BUILDINGS . . . . . . . . . . . . . . . . . . . . . . . . . 31 CHAPTER 6 BLAST RESISTANT ANALYSIS
AND DESIGN OF STRUCTURAL MEMBERS . . . .69
3.1 THREAT ASSESSMENT METHODS . . . . . . . 31
3.1.1 DOJ Report . . . . . . . . . . . . . . . . . . . . 31 6.1 MATERIAL PROPERTIES OF STEEL
3.1.2 GSA Security Criteria . . . . . . . . . . . . . 31 FOR BLAST DESIGN. . . . . . . . . . . . . . . . . . 69
3.1.3 ISC Security Design Criteria . . . . . . . . 31 6.1.1 Strength Increase Factor (SIF) . . . . . . . .69
3.1.4 Unified Facilities Criteria . . . . . . . . . . . 32 6.1.2 Dynamic Increase Factor (DIF) . . . . . . . 69
3.1.5 Department of State Criteria . . . . . . . . . 32 6.1.3 Dynamic Design Stress . . . . . . . . . . . . 69
3.1.6 Additional Criteria . . . . . . . . . . . . . . . 33 6.2 DESIGN CRITERIA FOR BLAST DESIGN . . .71
3.2 GOOD PRACTICE . . . . . . . . . . . . . . . . . . . 33 6.2.1 Load Combinations . . . . . . . . . . . . . . . 71
3.2.1 Exterior Considerations . . . . . . . . . . . . 33 6.2.2 Ultimate Strength . . . . . . . . . . . . . . . . 71
3.2.2 Interior Considerations. . . . . . . . . . . . . 34 6.2.3 Deformation Criteria . . . . . . . . . . . . . . 71
iii
6.3 FAILURE MODES . . . . . . . . . . . . . . . . . . . . 73 CHAPTER 8 RESISTANCE TO PROGRESSIVE
6.3.1 Breaching . . . . . . . . . . . . . . . . . . . . . 73 COLLAPSE . . . . . . . . . . . . . . . . . . . . . . . . . . . 117
6.3.2 Tension . . . . . . . . . . . . . . . . . . . . . . . 74
8.1 OVERVIEW . . . . . . . . . . . . . . . . . . . . . . . 117
6.3.3 Compression . . . . . . . . . . . . . . . . . . . 74
8.1.1 Progressive Collapse Definition . . . . . . 117
6.3.4 Shear . . . . . . . . . . . . . . . . . . . . . . . . 74
8.1.2 Brief Explanation of the
6.3.5 Flexure . . . . . . . . . . . . . . . . . . . . . . . 75
Design/Analysis Problem . . . . . . . . . . 117
6.3.6 Combined Forces . . . . . . . . . . . . . . . . 76
8.1.3 Basic Concepts . . . . . . . . . . . . . . . . . 117
6.4 DESIGN EXAMPLES . . . . . . . . . . . . . . . . . . 77
8.2 ANALYSIS AND DESIGN CODES
Example 6.1—Design of Structural Elements
AND GUIDELINES . . . . . . . . . . . . . . . . . . 118
Subject to Indirect Blast Loading . . . . . . . . . . . 77
8.2.1 Introduction . . . . . . . . . . . . . . . . . . . 118
Example 6.2—Design of Structural Elements
8.2.2 U.S. General Services Administration
Subject to Direct Blast Loading: Façade Girt
Guidelines . . . . . . . . . . . . . . . . . . . . 118
and Column . . . . . . . . . . . . . . . . . . . . . . . . . 87
8.2.3 Department of Defense Criteria . . . . . . 119
Example 6.3—Design of Structural Elements
8.2.4 British Standards . . . . . . . . . . . . . . . 121
Subject to Direct Blast Loading: Composite
8.2.5 Eurocode. . . . . . . . . . . . . . . . . . . . . 122
Roof Beam . . . . . . . . . . . . . . . . . . . . . . . . 102
8.3 ANALYTICAL APPROACHES TO
PROGRESSIVE COLLAPSE . . . . . . . . . . . . 122
CHAPTER 7 DESIGN OF CONNECTIONS FOR
8.3.1 Analysis Concepts. . . . . . . . . . . . . . . 122
BLAST RESISTANT STRUCTURES . . . . . . . . . . 113
8.3.2 Nonlinear Static Pushover Analysis:
7.1 GENERAL CONSIDERATIONS. . . . . . . . . . 113 Energy Balance Approach. . . . . . . . . . 125
7.2 DESIGN RESPONSIBILITY . . . . . . . . . . . . 113 8.3.3 Nonlinear Dynamic Analysis:
7.3 CONNECTION DUCTILITY . . . . . . . . . . . . 113 Time-History Approach . . . . . . . . . . . 127
7.4 CONNECTION STRENGTH . . . . . . . . . . . . 113 8.4 RECOMMENDATIONS . . . . . . . . . . . . . . . 129
7.4.1 Required Strength . . . . . . . . . . . . . . . 113 8.4.1 Prescriptive Recommendations . . . . . . 129
7.4.2 Available Strength. . . . . . . . . . . . . . . 114 8.4.2 General Design Recommendations . . . . 130
7.5 BOLTED CONNECTIONS . . . . . . . . . . . . . 114 8.4.3 Analytical Design Recommendations . . 131
7.5.1 Shear Connections . . . . . . . . . . . . . . 114 8.5 DESIGN EXAMPLE. . . . . . . . . . . . . . . . . . 131
7.5.2 Tension Connections . . . . . . . . . . . . . 114 Example 8.1—Analysis of Structural System
7.6 WELDED CONNECTIONS . . . . . . . . . . . . . 114 with Removal of an Interior Column . . . . . . . . 131
7.6.1 Filler Metals . . . . . . . . . . . . . . . . . . 114 8.6 EXAMPLE SUMMARY . . . . . . . . . . . . . . . 159
7.6.2 Quality Assurance . . . . . . . . . . . . . . . 114
7.6.3 Tension Applications . . . . . . . . . . . . . 114 SYMBOLS . . . . . . . . . . . . . . . . . . . . . . . . . . . 163
7.6.4 Flexural Applications. . . . . . . . . . . . . 115
REFERENCES . . . . . . . . . . . . . . . . . . . . . . . . . . 165
7.7 BRACING AND MOMENT-RESISTING
CONNECTIONS . . . . . . . . . . . . . . . . . . . . 115
iv
Chapter 1
Introduction
The purpose of this guide is to disseminate knowledge of Chapter 8 addresses basic progressive collapse concepts.
blast resistance and progressive collapse mitigation to the Progressive collapse design is independent of blast design
structural engineering community, presenting basic the- because progressive collapse may be caused by other pos-
ory with design examples so engineers and architects can sible events such as fire, accident, impact, etc. Examples
achieve simple and effective designs. demonstrating the determination of the structural response
Presently, security consultants with the assistance of the to progressive collapse are included.
owner evaluate the particular vulnerabilities of a given facil- The guide addresses only the behavior of structural steel
ity and determine the appropriate and acceptable level of under blast loading. It does not cover doors, windows, or any
security risk. The risk assessment study determines the loca- other structural material.
tion and the size of the explosive threat. The blast consul-
tants then calculate the blast pressures and review the design 1.1 HISTORY OF INCIDENTS
produced by the engineer of record. If the design is found
In years past, blast resistant design was typically only used
to be insufficient, the blast consultant recommends upgrad-
for facilities that either housed (or were in close proximity
ing the design and these revisions are incorporated into the
to) explosive material or were known as potential targets for
construction drawings. It is advisable to involve the security
attack. Munitions plants and storage facilities, strategic mili-
consultant and blast consultant as early as possible in the
tary and government facilities, and natural gas and petro-
planning and design process.
leum refineries are a few examples of facilities that might
There is enough information provided in this guide to
have been designed specifically to resist blasts. However,
allow practicing structural engineers with a background
the threat of bombings has increased in recent years. The
in structural dynamics to interact with blast consultants to
incidents described in the following are closely associated
produce effective designs. The engineer of record can then
with the evolution of the different security design criteria
proceed with the structural design based on the blast pres-
described in Chapter 3.
sures given by the blast consultant. As it is with any unusual
design, a peer review is a good idea and it is suggested that
1.1.1 Blast Incidents
the final design be reviewed by a qualified blast consultant
with experience in the design of blast resistant structures. While numerous bombing events have occurred throughout
This guide is divided into the following chapters: the world, a small number of these events over the past three
Chapter 2 addresses external blast explosions and is decades has had the largest impact on how the U.S. prepares
focused on the shock wave—not on fragment or projectile for, and responds to, such events.
loading. The chapter does not cover the loads generated by a Notable events include:
large blast in close proximity to the structure. • April 18, 1983—A suicide car bomber attacked the
Chapter 3 addresses the evolution of documents related to U.S. Embassy in Beirut, Lebanon, killing 63 people,
the design of buildings for blast loading and provides guid- 17 of whom were Americans.
ance on the relevant factors in protective building design.
• October 23, 1983—The U.S. Marine barracks in Bei-
Chapter 4 addresses methods of dynamic analysis, sim-
rut, Lebanon, were attacked by a suicide truck bomb
plifying multiple degrees of freedom into single degree of
killing 241 American military personnel.
freedom systems, and determining the dynamic response to
defined loads. It also explains the use of general structural • December 1983—Suicide truck bombers attacked
engineering software to solve simple multiple degree of free- the U.S. and French embassies in Kuwait killing 5
dom problems. and injuring 86.
Chapter 5 addresses the overall response of a building’s • September 20, 1984—The annex of the U.S. embassy
structural system to blast loading. in Beirut, Lebanon, was attacked with a truck bomb
Chapter 6 addresses member design, failure modes killing 24 and injuring the ambassador.
and design criteria including breaching, shear failure and
• December 21, 1988—A terrorist bomb destroyed
bending.
Pan Am Flight 103 over Lockerbie, Scotland, killing
Chapter 7 addresses steel connection design for blast
270 people.
loading.
Fig. 2-4. Relationship of reflected and free-field (side-on) pressure-time histories (Pr = reflected pressure; Pso = side-on pressure).
Fig. 2-5. Positive phase parameters for surface burst TNT explosions (DOD, 2008).
Given:
Evaluate the one-story steel building shown in Figure 2-8, Figure 2-9 and Figure 2-10 for a blast caused by a charge with W =
500 lb located just above ground level at the location shown in Figure 2-8. The building is 50 ft by 70 ft in plan and has a 15-ft
story height. The lateral force resisting system consists of rigid frames in the long dimension and braced frames on the exterior
walls in the short dimension. The roof is metal decking over structural steel purlins with a 0% slope.
Scaled Distance
For a stand-off distance of R = 50 ft and a TNT equivalent charge weight of 500 lb, the scaled distance is:
R
Z=
3
W
50.0 ft
=
3
500 lb
3
= 6.30 ft/lb
Since the explosive is located just above ground, the charge is considered to be a hemispherical surface burst explosion.
Fig. 2-11. Positive phase shock parameters for hemispherical TNT explosion on the surface at sea level (DOD, 2008).
Because the area under the two curves for equivalent reflected pressure in Figure 2-12 and Figure 2-13 must be equal, with the
area representing the impulse, the equivalent load duration is always less than the actual duration.
2 Ir
te,r =
Pr
2(246 psi ms)
=
79.5 psi
= 6.19 ms
2 I so
t e,so =
Pso
2(96.0 psi ms)
=
24.9 psi
= 7.71 ms
Figure 2-12 and Figure 2-13 show the predicted pressure time-history loading function and the equivalent triangular loading
function for reflected pressure and side-on pressure.
Front Wall
The typical blast load on a front wall is shown in Figure 2-14.
Side Walls
The typical blast load on a side wall is shown in Figure 2-16 as a plot of pressure versus time.
The blast load on the side wall is computed near the front corner of the building. For simplicity, the blast parameters (Pso, Iso,
qo, U) are computed using the scaled distance Z computed at the front wall on the centerline. The sidewall blast parameters are
given in Figure 2-13. A more detailed calculation would consider the reduction of blast pressure and impulse over the wall mov-
ing toward the rear.
Roof
Similar to the side walls, the blast load on the roof is computed near the corner of the building. For simplicity, the blast param-
eters (Pso, Iso, qo, U) are computed using the scaled distance Z computed at the front wall on the centerline. Thus, the roof param-
eters are the same as those for the sidewall and shown in Figure 2-13.
R
Z=
W3
120 ft
=
( 500 lb ) 3
= 15.1 ft/lb 3
Fig. 2-16. Typical side wall blast loading. Fig. 2-17. Typical rear wall blast loading.
Based on these parameters, the side-on blast pressure loads at the top of the rear wall are as shown in Figure 2-19.
Computation of the rear wall load is completed by evaluating the rise time and total duration as the blast wave sweeps down the
wall. The combined load is shown in Figure 2-20. In this figure, the equivalent positive phase load duration is used. The negative
phase of the load is not shown in this example. The negative phase is often ignored for simplicity and this is typically a conser-
vative approach. The span of the element parallel to the traveling wave is the building height, L1 = 15 ft. The free-field pressure
at the top of the wall is 4.60 psi with a time of arrival, ta = 66.0 ms. Thus, the time to peak pressure, t2, is the rise time plus the
time of arrival:
L
t2 = 1 + t a
U
15.0 ft
= + 66.0 ms
1.26 ft/ms
= 77.9 ms
The time at the end of the blast load (positive phase), tf, is the time to peak pressure plus the side-on load duration. In the follow-
ing calculation, the equivalent linear decay duration is used:
t f = t2 + te,so
= 77.9 ms + 19.1 ms
= 97.0 ms
Frame Loads
To establish the lateral response of the framing system, blast loads on the front wall and on the rear wall may be computed sepa-
rately to get the net combined load. The time of arrival of the loads to the walls must be taken into account. The results from the
computation of the loading for the frame along gridline B are shown in Figure 2-23.
The positive phase loading on the front (reflected) wall is often conservatively used alone for design for lateral response rather
than the net loading calculation results shown in this example.
Rear Wall
The pressure load applied to the rear wall is the side-on pressure applied to the rear wall which was shown in Figure 2-19. The
interaction of the pressure loading applied in each face of the frame is shown in Figure 2-22. Figure 2-23 shows the superposi-
tion to scale of the blast loads for the front and rear walls, accounting for time phasing. Note that positive pressure for the rear
wall is shown as a negative value in the graph because this load is in the opposite direction to that applied to the front wall load.
5
t= 77.9 ms, P= 4.6 psi
3
Pressure (psi)
t= 97.0 ms
t= 66 ms
-1
-2
0 20 40 60 80 100 120 140 160 180 200 220
Time (ms)
90
P= 79.5 psi
80
70
60
Pressure (psi)
50
40
30
ta = 15.6 ms
20
10
t= 21.8 ms
0
0 20 40 60 80 100 120 140 160 180 200
Time (ms)
Rear Wall
Front Wall
Charge
50ft 70ft
Frame B - Elevation
Direction of pressure
applied against wall
80
70
60
50
Pressure (psi)
40
Pressure load from Front wall
Pressure load from Rear wall
30
20
10
-10
0 20 40 60 80 100 120 140 160 180 200 220
Time (ms)
Fig. 4-1. Typical triangular loading function used to represent Fig. 4-2. Triangular load function with separate build-up and
blast pressure loading on structures. tail-down phases.
This is more commonly expressed in the form of cyclic fre- This single degree of freedom equation can be solved by
quency, f, in units of cycles per second or Hz, where f is numerical integration, sometimes performed using a spread-
given by: sheet program and algorithms available in a number of texts
on structural dynamics, such as Wang (1967). More com-
1 K monly, any of several structural analysis software programs
f = (4-6)
2π m can be used to perform this analysis directly. Regardless, it is
important to select an appropriate time increment over which
Rather than using natural frequency to characterize dynamic the integration is performed. The time increment needs to be
structural behavior, structural engineers often find it more short enough to provide an accurate solution yet long enough
Fig. 4-3. Dynamic load factor and triangular impulse loading for elastic response.
Given:
Consider a structure with a stiffness of 1,000 kip/in. and a natural period of vibration of 0.5 s. It is subjected to a simple triangular
impulsive loading, like that of Figure 4-1, with a peak force on the structure of 10,000 kips and a duration, td, of 0.05 s. Using
Figure 4-3 or Figure 4-4, as appropriate, determine the peak displacement and force in the structure.
Solution:
Step 1: Compute the force and displacement under statically applied peak loading:
T 0.5 s
=
td 0.05 s
= 10
Step 3: Enter Figure 4-3 for the pulse of Figure 4-1 with T/td = 10, the DLF is:
DLF = 0.280
Fig. 4-4. Dynamic load factor, triangular impulse loading, loading and unloading phases for elastic response.
The velocity imparted to the structure by the impulse pro- Fmax Δ max KΔ 2max
vides it with kinetic energy (WK). As the mass moves away WS,el = = (4-14)
2 2
from its at-rest position, it will strain the structure, dissipat-
ing the kinetic energy into stored strain energy (WS,el). When The strain energy that accumulates in the inelastic structure
the stored strain energy identically equals the imparted is:
kinetic energy, the structure will have reached its maximum
response. This condition is solved as follows: Fyield Δ yield
WS,inelastic = + Fyield ( Δ max − Δ )
WS,el = WK 2
= Fyield Δ yield ( μ − 1 2 ) (4-15)
Kxmax2 mv i2
=
2 2
m I Substituting Fyield = KΔyield and KΔmax2/ 2 = WS,el into Equa-
x max = vi = (4-11)
K Km tion 4-15, and rearranging, it can be shown that at maximum
displacement, Δmax, the strain energy stored by the elastic-
K plastic structure is related to the strain energy stored by the
Fmax = Kx max = I (4-12) elastic structure as follows:
m
2 (μ − 1 2 )
WS,elastic-plastic = WS,el (4-16)
where, xmax and Fmax are, respectively, the peak displacement μ2
and resisting force in the structure. This energy solution is
not exact, as it neglects the initial movement of the struc- At a ductility of 2, the elastic-plastic structure will store only
ture while the impulse is being applied. For structures in 75% of the strain energy stored by the elastic structure. At a
which the ratio of the structural period to the load duration is ductility of 4, the elastic-plastic structure will store only 44%
greater than about 10, this effect is negligible, and the energy of the strain energy stored by an elastic structure. Therefore,
solution is sufficiently accurate.
Given:
Determine the ductility demand for the structure of Example 4.1 if Fyield is 1,400 kips. From Example 4.1, the structure has a
stiffness of 1,000 kip/in. and a natural period of vibration of 0.5 s. It is subjected to a simple triangular impulsive loading, like
that of Figure 4-1, with a peak force on the structure of 10,000 kips and a duration, td, of 0.05 s. Find the peak displacement and
force in the structure, as well as the strength ratio and ductility demand.
Solution:
Step 1: Find the displacement and force in the structure if it remains elastic. Referring to Example 4.1, the maximum force for
elastic response is Fmax = 2,800 kips and the maximum displacement for elastic response is Δmax = 2.80 in.
Fig. 4-5. Force-displacement diagram for (a) elastic response and (b) elastic-plastic response.
Step 3: Find the ductility demand from Figure 4-6, at a strength ratio of 2.0:
μ ≈ 2.50
4.4.3 Energy Methods As with the elastic energy methods, this solution provides
sufficient accuracy when the ratio of structural period to
Energy methods for nonlinear structures are similar to those
impulse duration is 10 or greater.
for linear structures. However, the term for strain energy is
modified to account for the nonlinear behavior. Assuming
4.5 MULTIPLE DEGREE OF FREEDOM
elastic-plastic behavior, the ductility demand, μ, is deter-
STRUCTURES
mined as follows based on previous derivations:
Real structures do not have mass concentrated at a single
WS,elastic-plastic = WK point, but instead have mass distributed throughout the struc-
ture. Such real structures tend to have multiple modes of
mv 2i I2
K Δ 2yield ( μ − 1 2) = = vibration, each characterized by a unique shape and a unique
2 2m natural frequency or period. Figure 4-7 illustrates this con-
1⎛ ⎞ cept with modes of vibration for a structure with mass con-
I2
μ = ⎜1 + ⎟ (4-18) centrated at three points. Such structures are termed multiple
2 ⎜ KmΔ 2yield ⎟ degree of freedom (MDOF) structures.
⎝ ⎠
Load-mass
Mass factor, KM
factor, KLM
Con- Con-
Load cen- cen- Maximum Spring Dynamic
Loading Strain Factor, trated Uniform trated Uniform Resistance, Constant, Reaction,
Diagram Range KL Mass* Mass Mass* Mass Rm K V
8Mp 384EI
Elastic 0.64 — 0.50 — 0.78 0.39R + 0.11F
L 5L3
8Mp
Plastic 0.50 — 0.33 — 0.66 0 0.38Rm + 0.12F
L
4Mp 48EI
Elastic 1.0 1.0 0.49 1.0 0.49 0.78R – 0.28F
L L3
4Mp
Plastic 1.0 1.0 0.33 1.0 0.33 0 0.75Rm – 0.25F
L
6Mp 56.4EI
Elastic 0.87 0.76 0.52 0.87 0.60 0.525R – 0.025F
L L3
6Mp
Plastic 1.0 1.0 0.56 1.0 0.56 0 0.52Rm – 0.02F
L
*Equal parts of the concentrated mass are lumped at each concentrated load.
Source: Design of Structures to Resist the Effects of Atomic Weapons, U.S. Army Corps of Engineers, EM 1110-345-415 (USACE, 1957).
Fig. 4-8. Complex framing in which one flexural element is supported by other flexural elements.
ALGOR X X X X
ANSYS X X X X
CSI-ETABS X X X X
CSI-PERFORM X X X X
CSI-SAP X X X X
LARSA X X X X
LARSA 4D X X X X
LS-DYNA X X X X
NASTRAN X X X X
RISA-3D X X X
STAAD X X X X
(beam) differs by a factor of 2 from the period of the second- 4.6 SOFTWARE
ary element (girder), they can be treated as individual SDOF A number of common structural engineering software prod-
structures. One degree of freedom is the primary element ucts can be used to assist in the evaluation of structures and
(beam) analyzed with the blast load and unyielding supports structural elements for blast loads. Table 4-2 lists some of
and the other degree of freedom is the secondary element the more popular structural analysis software found in
(girder) loaded with the reactions from the primary element design offices and their ability to perform the various types
analysis. If the periods are within a factor of 2 of each other, of analyses described previously, as applicable to blast
more exact methods of analysis must be used. response evaluation. The types of elements available and the
flexibility of each of these software packages vary consider-
ably. Engineers should consult the software documentation,
available from the licensors, prior to using the software for
specific applications.
I 1 ⎛ I2 ⎞
Δ el = (5-6) Δ pl = ⎜⎜ + WS,el,max ⎟⎟ (5-9)
me K Fyield ⎝ 2me ⎠
where
For linear elastic-perfectly plastic behavior (Figure 5-2), the
Fyield = force that would cause the structure to yield
strain energy is:
In accordance with UFC 3-340-02 Table 5-8 (DOD, 2008),
I2 for blast response designed to avoid imminent collapse, the
WS,pl = Fyield Δ pl − W S,el,max = (5-7)
2m e deflection criteria used for frames limits the interstory drift
to H/25 (where H is the height between stories), and the
The maximum elastic strain energy occurs when the system maximum member end rotation (measured from the chord
yields:
1 2
WS,el = K Δ el
2
DISPLACEMENT
Fyield
WS,el =
1
Fyield × Δel ,max WS,pl = Fyield Δ pl − WS,el ,max
2
∆ el,max ∆ pl
DISPLACEMENT
F = KΔ (5-16)
π Fpeak
Feq,el = = k DRF Fpeak (5-17)
T td
Given:
In this example, the one-story building from Chapter 2, shown in Figure 5-4, is initially designed to support a 25 psf wind load
and to deflect less than H/400 under the design wind load. It is then designed for a 500-lb TNT equivalent blast load located 50 ft
away from the short façade. The blast load is assumed uniform on the front façade and any blast pressure on the rear façade is
ignored. The W-shapes are ASTM A992 steel and the rod braces are ASTM A572 Grade 50 steel. The gravity loads used in this
example are a 30 psf roof dead load and 40 psf façade dead load.
Simplifying assumptions include that the response of the building is predominantly in the first mode of vibration and that the load
path from the façade subjected to the blast, back to the lateral force resisting system, is completely rigid. This second assumption
is conservative because it implies that there is no energy dissipation along the load path.
Hand calculations are performed to determine the behavior of the system based on the equivalent blast pressure and the lateral
capacity of the building. An elastic-perfectly plastic equivalent system is used to obtain the plastic deflection of the system.
Secondly, a computer analysis is used to determine the behavior of the structure under blast load. A nonlinear time-history
dynamic analysis is performed for the triangular blast load shown in Figure 5-5 (see Chapter 2 for details). Nonlinear material
properties are modeled for the diagonal rods through the use of elastic-plastic axial hinges (see Chapter 6 for details).
For this example, SAP2000 was used. There are many other commercially available software packages that would perform this
analysis as well. See Chapter 4 for a discussion of available computer programs; in particular, Table 4-2 shows these packages
and their capabilities, including nonlinear dynamic analyses.
The results will show that the system is able to safely withstand the blast load if footings and connections are designed for the
overload. Design of columns, beams and connections will be introduced in Chapters 6 and 7.
The building dimensions are shown in Figure 5-4. The load used in the example is shown in Figure 5-5.
Solution:
Each of the braced frames takes only half of this load, hence the load modeled is 4.70 kips. Figure 5-6 shows the design and the
load used for wind design. Beams are modeled as rigid elements to capture the diaphragm behavior. All members are assumed
to be pinned at their ends.
Under wind load, the deflection is H/900 = 0.200 in. This is less than the H/400 limit set in the design statement. This system has
a stiffness of 4.70 kips/0.200 in. = 23.5 kip/in. The tension in the rod under wind loading is determined as follows:
where
α = tan −1 (15.0 ft 35.0 ft )
= 23.2°
The yield strength should be increased by a factor of 1.30 as discussed in Chapter 6, Equation 6-5. Hence, the available tensile
strength of the rod due to yielding is (tensile rupture will not control):
Tmax = 1.30 Fy Ag
⎡ ⎛ 3 4 in. ⎞2 ⎤
= 1.30 ( 50 ksi ) ⎢ π ⎜ ⎟ ⎥
⎢⎣ ⎝ 2 ⎠ ⎥⎦
= 28.7 kips
For this maximum tension, the maximum lateral load that the system can carry is:
⎛ T max ⎞
Fyield = ⎜ ⎟ ( applied wind load )
⎝ T wind ⎠
⎛ 28.7 kips ⎞
=⎜ ⎟ (4.70 kips)
⎝ 5.11 kips ⎠
= 26.4 kips
Considering a 30 psf roof dead load and a 40 psf façade dead load using the tributary area of the façade assuming half is carried
by the roof, the total weight involved in the movement of the bracing on one side is:
we = roof dead load + façade dead load
( 30.0 psf )( 50.0 ft ) ( 70.0 ft ) /2 ( 40.0 psf )(50.0 ft + 70.0 ft) (15.0 ft/2 )
= +
1, 000 ft/kip 1, 000 ft/kip
= 88.5 kips
According to Equation 5-10, the fundamental period of this single degree of freedom system is:
me
T = 2π
K
we
= 2π
gK
88.5 kips
= 2π
(386 in./s2 )(23.6 kip/in.)
= 0.619 s
The period of 0.619 s is approximately 100 times the duration of the load, te = 0.00620 s, as given in Figure 5-5; therefore, impul-
sive behavior is assumed.
The peak blast force applied per brace, Fpeak, correlating to the peak pressure (Pr from Figure 5-5), is half of the total blast pres-
sure, and is determined as follows:
Fpeak = (peak blast pressure)(tributary facade area)
⎡ (50.0 ft/2)(12 in./ft)(15.0 ft/2)(12 in./ft) ⎤
= 79.5 psi ⎢ ⎥
⎣ 1, 000 lb/kip ⎦
= 2,150 kips
Note that using only half of the story height (15.0 ft /2) is a departure from UFC 3-340-02 (DOD, 2008), which conservatively
uses the full façade area in calculating this force. This guide will follow the more conventional method used in wind design.
Based on Equation 5-12, the impulse is:
Fpeak t d
I=
2
( 2,150 kips ) ( 6.20 × 10− 3 s )
=
2
= 6.67 kip-s
= 67.7 kipss
The equivalent blast pressure is 14 times greater than the wind pressure. The maximum lateral force that the system can carry,
determined previously, is Fyield = 26.4 kips. The equivalent blast load is 67.7 kips, approximately 2.5 times larger than the actual
elastic capacity. Hence, the assumption that the structure remains elastic is not appropriate and plastic considerations are neces-
sary for design.
Plastic Deflection
Considering an elastic-plastic behavior of the building with the lateral capacity previously defined, Equation 5-9 gives the plastic
displacement as:
1 ⎛ I2 ⎞
Δ pl = ⎜⎜ + WS,el,max ⎟⎟
Fyield (5-9)
⎝ 2me ⎠
where
I = 6.67 kip-s (as determined previously)
me = we /g
From Equation 4-14, the maximum elastic energy is:
Fyield Δ el,max
WS,el,max =
2
( 26. 4 kips )(1.12 in. )
=
2
= 14.8 kip-in.
1 ⎛ I2 ⎞
Δ pl = ⎜⎜ + WS,el,max ⎟⎟ (5-9)
Fyield ⎝ 2me ⎠
1 ⎡ ( 6.67 kip-s )2 ⎤
= ⎢ + 14.8 kip-in.⎥
26.4 kips ⎢ 2 ( 88.5 kips 386 in./s ) ⎥
⎣ ⎦
= 4.24 in.
The structure as designed for wind is able to resist the blast pressure with an acceptable deflection and ductility ratio provided
that the connection has enough strength to allow the member to reach its yield strength. Chapter 6 will introduce resistance and
ductility criteria to design and check the structural elements.
Computer Analysis
The structure shown in Figure 5-6 is modeled for computer analysis (using SAP 2000) with both rods included to capture the
behavior of the structure on rebound. The rods are modeled as tension-only elements with nonlinear axial hinges (see Chapter 6
for further discussion). The beams are modeled as rigid to approximate rigid diaphragm behavior and a rigid path to the lateral
system. As discussed earlier, this is a conservative assumption that neglects energy dissipation.
The maximum tensile strength of the rod was computed previously to be 28.7 kips. It is assumed to have no compressive strength.
The total weight of the system was computed as 88.5 kips. This weight is distributed linearly over the beams in the frame.
The blast load applied on the structure is given in Figure 5-7. A general material nonlinear time-history analysis was performed.
The structure was analyzed for the effect of the blast load alone, excluding dead load. From the computer modal analysis, the
period of the first elastic mode is 0.620 s. This is the same as the value obtained by the hand calculations performed previously.
However, the period of the blast response (T ≈ 1 s) is much longer than the calculated elastic period of 0.619 s due to plastic
behavior.
Time dependent displacement results are shown in Figure 5-8. The maximum displacement is 4.24 in., which is less than H/25 =
7.20 in. This is the same value obtained previously. The structure deforms plastically. The displacement past the maximum elastic
displacement is permanent plastic deformation of Rod 1.
The axial forces in the rods are shown in Figure 5-9. Note that Rod 2, which was initially inactive, is activated in the rebound of
the structure. The energy absorbed by Rod 1 includes both an elastic and a plastic component. The energy absorbed during the
rebound by Rod 2 is elastic and is equal to the elastic energy released by Rod 1 as it unloads. As shown in Figure 5-9, from 0.430
to 0.700 s neither rod is active. Due to the permanent plastic deformation shown in Figure 5-8, both rods are in compression at
this time. Rod 2 is compressed because the structure has not reached the original elastic equilibrium point. Rod 1 has been per-
manently elongated and thus is compressed before reaching the original elastic equilibrium point.
Figure 5-10 shows the axial force in Column 1 and Column 2. Figure 5-11 shows the reactions at the supports. For comparison,
the elastic reactions under wind load are ±2 kips. The foundation reaction under blast load is 12 kips, 6 times greater than the
wind reaction, but with less than 0.6 s of duration. Note that the foundations are subject to both downward and upward load. In
Chapter 6, the columns are designed to remain elastic to avoid the failure of the structure.
The ductility in the system, μ, can be computed as the ratio between the plastic and the elastic yield displacement as determined
previously. Chapter 6 defines criteria to classify the blast behavior of the structural system based on ductility ratios.
Fig. 5-8. Displacement results for the movement of the roof under blast load.
Given:
The three-story building shown in Figure 5-12 is designed for blast loading. The lateral system of the building is formed by a
chevron braced frame on each of the 150-ft sides of the building. This building is designed for 25 psf wind load with a deflection
limited to H/400. The resulting structural design is shown in Figure 5-13. The blast load affects the 120-ft-wide façade and has
the same triangular time-history blast pressure used in Example 5.1. The HSS diagonal braces are ASTM A500 Grade C material.
All of the floors are modeled as rigid diaphragms. This neglects the axial load in the beams that can affect the beam design.
The façade and the load path to the lateral system are assumed rigid. As described earlier, this is a conservative assumption as it
neglects the energy absorbed by the façade and the load path.
Lumping the building’s mass at each floor, the mass moving with each frame per floor is:
Roof = 450 kips (50 psf)
3rd Floor = 900 kips (100 psf)
2nd Floor = 900 kips (100 psf)
Solution:
T = 2π
∑ wi Δ2i
g ∑ fi Δ i
( 450 kips )( 0.120 in.)2 + ( 900 kips )( 0.100 in. )2 + ( 900 kips )( 0.0600 in.) 2
= 2π
( 386 in./s ) ⎡⎣( 9.00 kips )( 0.120 in.) + (18.0 kips )( 0.100 in.) + (18.0 kips )( 0.0600 in.) ⎤⎦
= 0.695 s
Using the same impulse time history function as Example 5.1, using Figure 5-5, the total peak load in this three-story braced
frame is:
This force is equivalent to a uniform pressure on the façade for each braced frame of:
As with the previous example, the equivalent blast pressure is much greater than the design wind pressure of 25 psf. The elastic
yield point is identified as the total load that makes the first diagonal fail in compression. In Figure 5-13, the force in the bottom-
most diagonal is 32 kips and the total wind shear is 45 kips. Assuming a blast compressive strength for a 17-ft-long HSS6×6×¼
of 165 kips (see Chapter 6 for discussion of dynamic element strength), the lateral strength is:
The equivalent blast force is 2.5 times larger than the actual elastic strength of the system; therefore, a plastic analysis is required.
As the braces are behaving plastically, design the system as a special concentrically braced frame (SCBF) as defined in the AISC
Seismic Provisions for Structural Steel Buildings (AISC, 2010b). Therefore, as an SCBF, the post-buckling and post-yielding
strength of the braces is explicitly modeled and considered through the use of nonlinear axial hinges in the next section. Note
that, however, the structure is not subject to high seismic loading.
Computer Analysis
The structure shown in Figure 5-13 is modeled as a multiple degree of freedom system in a structural software package (SAP2000)
to carry out the nonlinear time-history analysis. The strength and ductility properties for the different elements are calculated in
Chapter 6. Columns are assumed to remain elastic and are designed as such in Chapter 6.
Using the dynamic strength of the members to be defined in Chapter 6, the HSS6×6×¼ diagonals have 340 kips of tensile
strength and 165 kips of compressive strength (see Example 6.1). The W12×35 beams have 262 kip-ft of flexural strength. The
diagonals are modeled in the computer analysis with nonlinear axial hinges at their centers and the beams with nonlinear moment
hinges at their midspan. The hinge properties are in accordance with FEMA 356 (FEMA, 2000b). See Example 6.1 for the deriva-
tion of these hinge properties. Member ends are assumed pinned.
Figure 5-14 shows the time-history load applied for the nonlinear time-history analysis. The total peak pressure is distributed to
the different floors by tributary façade area. For modeling purposes, a leaning column is added to the side of the building and the
floor mass is lumped at these nodes, as shown in Figure 5-14. This will prevent high frequency axial modes from developing in
the results.
Figure 5-15 shows the nonlinear time-history displacement at each floor of the building. The structure goes from an original
elastic period of 0.695 s to an apparent period of roughly 3 s.
The maximum deflections are shown in Table 5-1, where Δ is the total deflection, δ is the interstory deflection, and H is the
story-to-story height. The deflection between the 1st (ground) floor and 2nd floor is the most critical. With δ = H/46 < H/25 and
a maximum end rotation of θ = 1.2° < 2.0° (based on a maximum beam deflection of 3 in.), it meets the deflection and rotation
criteria, as stated in Section 5.1. Hence, this structure, as designed for wind, safely withstands the blast load.
Figure 5-16, Figure 5-17, Figure 5-18 and Figure 5-19 show the response of the different elements that either start to yield or
completely fail. Figure 5-20 shows the deflection of the second floor beam at midspan. Figure 5-21 shows a key elevation locating
these elements. Figure 5-16 shows that the first floor compressed diagonal fails in compression. After this diagonal loses all of
its strength, the other elements remain elastic or start to yield without considerable ductility. Hence, the permanent deformation
The purpose of this chapter is to define the structural proper- be used to model and analyze the structures used as exam-
ties used for nonlinear dynamic models, to understand the ples in this chapter. Table 4-2 shows a list of some of this
analysis results, and to classify or design the structural ele- software.
ment for the desired criteria. Failure modes will be explained
to determine the mechanism of collapse for the structural 6.1 MATERIAL PROPERTIES OF STEEL
element. Deformation criteria will be established to classify FOR BLAST DESIGN
the response of the element. The primary failure modes used
to design structural elements are breaching, tension-com- 6.1.1 Strength Increase Factor (SIF)
pression, bending, shear, and axial-bending interaction. The
For steel grades of 50 ksi or less, the average yield stress of
secondary modes of failure such as brittle failure caused by
steels currently produced is approximately 10% larger than
local stress concentrations, welding, low temperatures, etc.
the stress specified by ASTM; therefore, for blast design the
are beyond the scope of this guide.
specified minimum yield stress should be multiplied by a
Member failure is defined through support rotation and
strength increase factor, SIF, of 1.10. For higher grades this
ductility ratio. This definition is intended for the single
average is smaller than 5%; therefore, no factor is used on
degree of freedom systems (SDOF) and simplified multiple
those grades. Ultimate strength is not factored in any case.
degree of freedom systems (MDOF). Strain based failure
criteria may also be justified if strains are calculated using
6.1.2 Dynamic Increase Factor (DIF)
finite element methods that characterize the material prop-
erties, the details of construction, and many of the pos- Steel mechanical properties vary with the time rate of strain.
sible failure mechanisms. As used in this guide, ductility is As compared with the static values normally used in design,
defined as the ratio between the maximum deflection and the properties vary for dynamic loading as follows:
the maximum elastic deflection. This parameter is smaller • The yield point increases substantially.
than one if the behavior is elastic and larger than one if the
• The ultimate tensile strength increases slightly.
behavior is plastic.
The blast design of several building elements is exempli- • Modulus of elasticity does not vary and the elonga-
fied in this chapter. The examples follow the load path from tion at rupture either remains constant or is slightly
the façade to the lateral load resisting system. Examples reduced.
include the design of a façade girt, façade column, perimeter The factor used to modify the static stress due to dynamic
beam, and some elements of the lateral load resisting sys- load is the dynamic increase factor, DIF. These factors are
tems designed in Chapter 5. defined in Design of Blast Resistant Buildings in Petrochem-
This chapter does not consider the energy dissipated along ical Facilities (ASCE, 2010b), UFC 3-340-02, and PDC
the load path and assumes that the elements are absorbing TR-06-01 (USACE, 2008).
all of the energy, which is a conservative assumption. The The values summarized in Table 6-1 are based on an aver-
designer should give attention to the reactions and the con- age strain rate of 0.10 in./in./s which is characteristic of low
nections along the load path from one member to another as pressure explosions. These values are appropriate for most
these will be required to transfer the full energy of the sys- conventional explosive load environments. Higher values of
tem. It will be shown that on the order of 80% of the energy strain rate give larger values of DIF. UFC 3-340-02 provides
from a blast load can be dissipated along the load path. values of DIF for different average strain rates.
Blast mitigation design should be integrated with the
overall structural design, not left to some later stage in 6.1.3 Dynamic Design Stress
design, as it may increase the stiffness or the mass of the
Based on the expected ductility ratio and/or the damage
structure which would affect the response of the structure to
allowed in the structural element, different dynamic yield
other loads.
points are defined by UFC 3-340-02. If the ductility ratio
Various programming tools or commercial software can
is smaller than or equal to 10, the dynamic design stress is:
DEFLECTION
DESIGN SECTION
MODULUS S (S+ Z)/2 Z
DESIGN STRESS
< fdy fdy fdy +(fdu-fdy )/4
as a function of the damage in the building. UFC 3-340-02 members, as the connections between members may have
also classifies the response as a function of the protection different criteria. The ductility limits are linked to a given
provided by the structural elements. This guide uses a low mode of response. A flexural ductility is different from shear
level of protection that implies a high response because the or tension.
intent is to provide a design that avoids imminent collapse There are several other sources for response criteria in
but allows substantial damage. addition to the NYC Building Code. It is up to the designer
For different structural elements, this guide follows the to determine which criteria are most applicable and should
response criteria shown in Table 6-2, which conforms to be used. Tables from UFC 3-340-02, Design of Blast Resis-
Table 1626.9.3 of the New York City (NYC) Building Code tant Buildings in Petrochemical Facilities (ASCE, 2010b),
(NYCBC, 2008). The rotation criteria in Table 6-2 refer to and FEMA 356 (FEMA, 2000b) are included in this guide
support rotations. The criteria are defined for the behavior as Table 6-3, Table 6-4 and Table 6-5. The 2010 edition
of a single element. These criteria apply to the design of of Design of Blast Resistant Buildings in Petrochemical
Facilities (ASCE, 2010b) also has medium response and (15.0 ft)(12 in./ft)(tan10°) = 31.7 in.
high response criteria in addition to the low response criteria
A ductility of 1 implies a deflection of 4 in. Thus, a ductility
included in Table 6-3. It is recommended that these sources
of 20, again taken from Table 6-2, would require a deflection
be consulted for more information in determining design
of 80 in. Obviously, this level of ductility is not achievable
criteria.
and the support rotation controls. These, and values for the
In Table 6-5, ΔT refers to the axial deformation at expected
other criteria sources, are tabulated in Table 6-6.
yielding load. It is important to note that FEMA is the only
Other criteria exist, such as U.S. Army Corps of Engineers
reference to have criteria for a tension member. This guide
Protective Design Center (USACE PDC) criteria, which are
will therefore use these criteria for tension elements. Because
more restrictive and are usually used in commercial and
this guide is concerned with a low level of security and pre-
governmental buildings. In Chapter 5, the maximum lateral
vention of collapse, it is recommended that the collapse
displacement was limited to H/25 and a maximum plastic
prevention (CP) values be used. For information on distin-
rotation of 2° was imposed to the overall behavior of the
guishing between primary and secondary members, see the
lateral system.
discussion in FEMA 356 (FEMA, 2000b). FEMA 356 also
includes acceptance criteria for other element types.
6.3 FAILURE MODES
To help put these criteria in perspective, examine the val-
ues for a beam with a 30-ft-long span that starts to yield at 6.3.1 Breaching
midspan when the deflection is 4 in. A support rotation of For blast loads in contact (or in very close proximity) to an
10°, as indicated in Table 6-2 as the maximum rotation for element, the element will be breached before responding in
steel beams, implies a deflection of: a flexural manner due to the high pressure produced by the
explosion. If the scaled distance as defined in Chapter 2 is Axial hinge properties for these elements are elastic-plastic.
below 2, it is possible for the element to breach before the If necessary, for computational stability, a hardening slope
overall response of the structural element starts. of 0.1% can be used in the plastic area. Examples of these
For very close charges, temperature and the shock wave hinges will be shown later in this chapter.
are important, in addition to the airblast. An explosion in
direct contact, with a wall will interact directly with it and 6.3.3 Compression
will induce a shock wave inside the wall. The speed and
The available strength of compression elements, based on
magnitude of the shock waves can cause the wall to crack
AISC Specification Chapter E is governed by the following
internally. For example, when the shock wave reaches the
equation:
opposite face of a concrete wall, a section of concrete may
scab or separate from the wall due to the energy of the shock ϕPn = ϕFcr Ag (6-10)
wave exceeding the tensile strength of the material.
For the preliminary design of the slab thickness, this where
Design Guide uses an experimental graphic found in UFC Ag = gross area of the element
3-340-02 and shown in Figure 6-3. The input for this graph is Fcr = critical stress
the thickness (T, feet), the stand-off (R, feet) and the charge ϕ = 0.90
(W, equivalent pounds of TNT). Depending on the perfor-
For blast loading:
mance requirements, the slab may be designed to allow a
Fcr = critical stress determined in accordance with AISC
local breach in the bay closest to the blast, yet withstand the
Specification Chapter E substituting ƒds for Fy
blast in adjacent bays. Typically, the stand-off and charge
ϕ = 1.00, the resistance factor defined in Section 6.2.2.1
weight is specified in the data provided by the blast consul-
tant. Other breaching curves can be found in UFC 3-340-01 For a tension-compression hinge, the recommendations of
(DOD, 2002), however, this document is export controlled FEMA 356 (FEMA, 2000b) can be used, if no other proper
and not available to the public. criteria are defined.
Fig. 6-5. Design moment for beams with moderate ductility (DOD, 2008).
The elastic-plastic behavior for each element used in Chapter 5 is defined in the following. The elements designed in this section
all carry blast load, but are not subject to direct pressure from the blast. Example 6.2 will examine the effects of elements directly
subject to blast pressure. The criteria used in this example to define the strength of the elements is based on UFC 3-340-02,
unless otherwise noted. For element design, assume that there is no energy dissipation along the load path. A more accurate and
less conservative procedure is to use the effect of the dynamic reactions of one member on another member along the load path.
When using member reactions to load adjacent members, natural frequencies should be compared. If the natural frequencies are
close, a simultaneous solution is required to account for the interaction between the two members. If the period of the primary
element (i.e., beam) is at least twice the period of the secondary element (i.e., girder), they can be treated as individual single
degree of freedom (SDOF) structures on unyielding supports. If not, a multiple degree of freedom (MDOF) solution of the same
system is required similar to what was discussed in Section 4.5.
The tension rods are ASTM A572 Grade 50 steel, the W-shapes are ASTM A992 steel, and the HSS shapes are ASTM A500
Grade C.
Solution:
From Table 2-4 of the AISC Steel Construction Manual (AISC, 2011a), hereafter referred to as the AISC Manual, the material
properties are:
ASTM A572 Grade 50
Fy = 50 ksi
Fu = 65 ksi
ASTM A992
Fy = 50 ksi
Fu = 65 ksi
ASTM A500 Grade C
Fy = 50 ksi
Fu = 62 ksi
From AISC Manual Table 1-1 and Table 1-12, the geometric properties are:
HSS6×6×¼
Ag = 5.24 in.2
r = 2.34 in.
W12×53
Ag = 15.6 in.2
ry = 2.48 in.
W12×35
Ix = 285 in.2
Sx = 45.6 in.3
Zx = 51.2 in.3
For a ductility ratio μ > 10, the dynamic design stress, defined in Equation 6-2, is:
⎛ fdu − fdy ⎞
fds = fdy + ⎜ ⎟⎟ (6-2)
⎜ 4
⎝ ⎠
⎡1.05 ( 65 ksi ) − 60.5 kssi ⎤
= 60.5 ksi + ⎢ ⎥
⎣ 4 ⎦
= 62.4 ksi
where
Comparing these dynamic design stresses to the simplified value given by Equation 6-5:
Note that the simplified value is slightly less conservative. The remainder of the calculations throughout this chapter will utilize
the simplified expression of Equation 6-5.
The dynamic design shear stress defined in Equation 6-4, for a ductility ratio μ = 10, is:
The available tensile strength of the rod is determined from Equation 6-9, as follows:
ϕPn = ϕFy Ag (6-9)
ϕPn = ϕ fds Ag
⎡ ⎛ w in. ⎞2 ⎤
= 1.00 ( 65.0 ksi ) ⎢ π ⎜ ⎟ ⎥
⎢⎣ ⎝ 2 ⎠ ⎥⎦
= 28.7 kips
The hinge properties used in Example 5.1 for this rod follows the elastic-plastic curve shown in Figure 6-6. A slight slope is
included in the plastic region for computational convergence purposes.
The maximum horizontal displacement obtained in Example 5.1, Figure 5-8, is 4.24 in. This is less than the horizontal deflection
limit of H/25 = 7.20 in. for frame structures. The maximum rotation at the base of the column for this displacement is:
⎛Δ⎞
α = tan −1 ⎜ ⎟
⎝H⎠
⎡ 4.24 in. ⎤⎥
= tan −1 ⎢
⎢⎣ (15.0 ft) (12 in./ft)⎥⎦
= 1.34° < 2°
For small displacement formulation, the elongation in the rod can be determined. First, the angle of the diagonal is determined as:
Then, the elongation, assuming the top of the column moves horizontally is:
Δ L = Δ cos θ
= ( 4.24 in.) ( cos 23.2°)
= 3.90 in.
fds L
Δ el =
E
( .0 ksi )( 38.0 ft )(12 in./ft )
65
=
29, 000 ksi
= 1.02 in..
3.86 in.
μ=
1.02 in.
= 3.78
This is the same value as the one obtained for the whole lateral system in Example 5.1. From Table 6-5, the allowable defor-
mation is 9ΔT. This is analogous to a response ratio of μ = 9. With a ductility of 3.78, the demand is less than the capacity and
therefore acceptable.
2π 2 E
Cc =
fds
This can be compared to the AISC Specification Section E3 limit for inelastic behavior of:
KL E
≤ 4.71
r Fy
For blast loading, Fy = fds = 65.0 ksi and ϕ = 1.00, and the limit is:
Note that this limit is significantly less than the suggested maximum slenderness of 200 given in the User Note in Section E2 of
the AISC Specification. Assuming K = 1.0, the slenderness of this element is:
π2 E
Fe = (Spec. Eq. E3-4)
( KL / r )2
π 2 ( 29, 000 ksi )
=
(87.2 )2
= 37.6 kssi
⎛ Fy ⎞
⎜ F ⎟
Fcr = ⎜0.658 e ⎟ Fy (Spec. Eq. E3-2)
⎝ ⎠
For blast loading, Fy = fds = 65.0 ksi, and the critical stress is:
⎛ fds ⎞
⎜ ⎟
Fcr = ⎜ 0.658 Fe ⎟ fds
⎜ ⎟
⎝ ⎠
⎛ 65.0 ksi ⎞
= ⎜ 0.658 37.6 ksi ⎟ ( 65.0 ksi )
⎜ ⎟
⎝ ⎠
= 31.5 ksi
Note that the buckling stress used here corresponds to the AISC Specification. This differs from UFC 3-340-02 which uses a
buckling stress corresponding to the 1989 Specification (AISC, 1989).
Hence, for blast loading, with ϕ = 1.00, the available compressive strength of the diagonal brace is:
(
= 1.00 ( 31.5 ksi ) 5.24 in.2 )
= 165 kips
The available tensile yielding strength of the diagonal brace, from AISC Specification Section D2, is:
ϕPn = ϕFy Ag (6-9)
For blast loading, Fy = fds = 65.0 ksi and ϕ = 1.00, and the available tensile strength is:
ϕPn = ϕ fds Ag
(
= 1.00 ( 65.0 ksi ) 5.24 in.2 )
= 341 kips
With these tensile and compressive capacities, the tension-compression hinge properties used in Example 5.2 are derived as put
forth in Figure 6-7, Table 5-6 and Table 5-7 of FEMA 356 (FEMA, 2000b). These hinge properties are used to model plastic
hinges in the braces to allow for a nonlinear plastic analysis of the structure, as shown in Chapter 5.
From Table 5-7 of FEMA 356, the HSS section has the modeling parameters shown in Table 6-7 for tension and compression. ΔT
refers to the axial deformation at expected tensile yielding load, and Δc is the axial deformation at expected buckling load. Based
on E = 29,000 ksi, its relationship to stress and strain, i.e., ε = σ/E and the definition of strain, i.e., ε = ΔL/L, the values of ΔT and
Δc can be determined. Setting these two equations for ε equal to each other and solving for ΔL results in:
σL
ΔL =
E
For tension:
σL
ΔT =
E
( 65.0 ksi )(17.0 ft )(12 in./ft )
=
29, 000 ksi
= 0.457 in.
σL
Δc =
E
( 31.5 ksi )(17.0 ft )(12 in./ft )
=
29, 000 ksi
= 0.222
Using the allowable strain-hardening slope of 1.5% of the elastic slope for tension produces a maximum stress of approximately
1.15 of the yield stress (Fell et al., 2006). For compression, to avoid computational instabilities, use a 0.1% hardening slope. With
this information, the hinge is compiled as Figure 6-8.
Results from Section 5.3, Example 5.2, indicate that only one diagonal at the first floor fails in compression while the others only
started yielding. The structure remains stable despite the failure of one brace in compression due to its redundancy.
2π 2 E
Cc =
fds
As determined previously, for blast loading, Fy = fds = 65.0 ksi, and the limit becomes:
The buckling stress is defined in Equation 6-14 and for this case is:
π2E
Fe = (Spec. Eq. E3-4)
( KL / r )2
π 2 ( 29, 000 ksi )
=
(72.6 in.) 2
= 54.3 kssi
For blast loading, Fy = fds = 65.0 ksi, and the critical stress is:
⎛ fds ⎞
⎜ ⎟
Fcr = ⎜ 0.658 Fe ⎟ fds
⎜ ⎟
⎝ ⎠
⎛ 65.0 ksi ⎞
= ⎜ 0.658 54.3 ksi ⎟ 65.0 ksi
⎜ ⎟
⎝ ⎠
= 39.4 ksi
For blast loading, ϕ = 1.00, and the available compressive strength is:
ϕPn = ϕ Ag Fcr
ϕPn = ϕFy Ag
For blast loading, Fy = fds = 65.0 ksi and ϕ = 1.00, and the available tensile strength is:
ϕ Pn = ϕ fds Ag
The maximum axial compressive load in Column 1 due to the blast load is 11.5 kips, as shown in Figure 5-10. Considering a
17.5 ft by 12.5 ft tributary area, a dead load of 50 psf, a live load of 30 psf, and the load combination given by Equation 6-7, the
total column load is:
Pu = 1.0D + 0.25L + 1.0B
= 1.0 (17.5 ft )(12.5 ft )( 0.050 ksf ) + 0.25 (17.5 ft )(12.5 ft ) ( 0.030 ksf ) + 1.0 (11.5 kips )
= 24.1 kips < ϕPn
The maximum axial compressive load in Column 2 due to the blast load is 11.5 kips, also shown in Figure 5-10. Considering a
35 ft × 12.5 ft tributary area, the total axial compressive load for this column is:
Pu = 1.0D + 0.25L + 1.0B
= 1.0 ( 35 ft )(12.5 ft )( 0.050 ksf ) + 0.25 ( 35 ft )(12.5 ft ) ( 0.030 ksf ) + 1.0 (11.5 kips )
= 36.7 kips < ϕPn
The maximum compression is 24.1 kips in Column 1 and 36.7 kips in Column 2. These are significantly below the buckling
load of the column (ϕPn = 615 kips). The maximum stresses in the columns are 1.5 ksi and 2.4 ksi, respectively, therefore, the
columns remain elastic.
⎛S+Z ⎞
M pʹ = fds ⎜ ⎟ (6-12)
⎝ 2 ⎠
⎛ 65.0 ksi ⎞ ⎛ 45.6 in. + 51.2 in. ⎞
3 3
=⎜ ⎟ ⎜⎜ ⎟⎟
⎝ 12 in./ft ⎠ ⎝ 2 ⎠
= 262 kip-ft
where fds was determined previously as 65.0 ksi for blast loading using the simplified Equation 6-5.
M ʹp L2
Δ=
12 EI
( 262 kip-ft )(12 in./ft )( 24 ft )2 (12 in./ft )2
=
12 ( 29, 000 ksi ) ( 285 in.4 )
= 2.63 in.
⎛ Δ ⎞
θ = tan − 1 ⎜ ⎟
⎝L/2⎠
⎡ 2.63 in. ⎤
= tan − 1 ⎢ ⎥
⎢⎣ ( ft )( in./ft ) ⎥⎦
24 12 / 2
= 1.05°
Note that the hinge rotation is double this support rotation (see Figure 6-2). These parameters define an elastic-perfectly plastic
moment-rotation curve. Since many programs have convergence problems with a perfectly plastic zone, a sloped line should be
introduced in this plastic region. This slope shown in Figure 6-9 is based on fds determined from Figure 6-1 and Equation 6-2 for
ductility greater than 10. Therefore, the dynamic design stress is:
fdu − fdy
fds = fdy + (6-2)
4
where
and
Mult = fds Z x
(
θult = min 10°, 20 θ )
= min (10°, 21°)
= 10°
Again, note that hinge rotation is double this support rotation. The moment-rotation diagram for the beam hinge is given in
Figure 6-9.
Figure 5-20 shows the time-history deflection at midspan for the second-floor beam. The maximum deflection results in the first
cycle with a value of 3 in. For this deflection the ductility is smaller than the maximum value defined in Table 6-2 for beams:
3.00 in.
μ=
2.63 in.
= 1.14 < 20
Note, also, that the ductility ratio is smaller than 3, indicating that the initial assumption was valid. If, conversely, the ratio was
larger than 3 this process would have to be repeated with the correct Mp.
Example 6.2—Design of Structural Elements Subject to Direct Blast Loading: Façade Girt and Column
Given:
In the previous examples, the element behaviors were defined and included in the structural models used in Examples 5.1 and
5.2. The results from these examples were used in this chapter to check the adequacy of the members to support the loads defined
in Chapter 5. These particular elements were not directly exposed to blast loads. In this example, the elements are designed and
analyzed based on the blast load applied directly to them. Elements designed in Example 5.1 are simplified into an SDOF model
and will be redesigned according to the requirements defined in this chapter. All steel is ASTM A992 material. Specifically, the
following elements are designed:
(a) Façade Girt Design: Design an 8-in.-deep section. This element has been designed for wind as a MC8×20 of ASTM A992
material, with a deflection limitation of L/260.
An introduction to the MDOF-to-SDOF simplification method was presented in Chapter 4. Figure 6-10 provides an overview of
this method.
Note that the sample procedure equations given in Figure 6-10 are based on a simply supported beam with load and mass uni-
formly distributed. Equations for other boundary conditions were defined in Chapter 4. The same transformation as shown in
Figure 6-10 can be performed by multiplying only the mass by the load-mass factor (Biggs, 1964). Figure 6-11 summarizes the
SDOF solution.
KM
MSDOF = MK LM = M (6-14)
KL
This can be seen by starting from the simple force equilibrium equation, and applying the transformation shown in Figure 6-10.
The simple force equilibrium equations are:
F = ku + Mu (6-15)
KL F = kKL u + MKM
u (6-16)
The KLM approach is simpler because it only uses one transformation factor and is standard practice in blast analysis/design.
Here, the load factor, KL, and the mass factor, KM, are used because they have a more physical interpretation.
The connection should be detailed for this high rotational capacity. Connection design will be discussed in Chapter 7.
In this example, the process shown in Figure 6-10 and Figure 6-11 will be followed. Begin by converting the mass and load from
the MDOF system to an equivalent SDOF system that is easy to solve. Within this equivalent SDOF system, find the necessary
yield force from the total energy and maximum deflection in the system. Next, convert this yield force back to the MDOF system
and find the maximum moment. With this moment, design a section with the required plastic section modulus using the AISC
Specification.
The self-weight of the façade is 40 psf. Therefore, the weight of the façade supported by the system is:
The self-weight is important for the calculation of the mass involved in the movement of the system. The girt is included in
the dynamic and modal calculations; however, it provides only lateral resistance to the wind loading and supports only its own
weight.
From Example 2.1, the blast lasts for a duration of 6.19 ms and peaks at a pressure of 79.5 psi. The load associated with this blast
pressure spread about the tributary area of the girt is given as:
( )
Fpeak = ( 79.5 psi ) 144 in.2 ft 2 ( 5.00 ft )( 25.0 ft ) (1,000 lb/kip )
= 1, 430 kips
Based on the SDOF simplification used for mass and loads uniformly distributed in the plastic range, the load and the stiffness
are multiplied by the load factor, KL = 0.50, and the mass is multiplied by the mass factor, KM = 0.33, as found in Table 4-1.
Therefore, the load and mass parameters used in the discrete system are:
wSDOF = KM w
= 0.33 ( 5.00 kips )
= 1.65 kips
The equivalent impulse due to the 6.19 ms blast in this SDOF system is determined from Equation 5-12 as follows:
Fpeak,SDOF td
I SDOF =
2
( 715 kips ) (6.19 × 10− 3 s)
=
2
= 2.21 kip-s
The total energy produced by the blast load in the SDOF system is determined from Equation 5-3:
2
I SDOF
WP,SDOF =
2 m SDOF
( 2.21 kip-s )2
=
⎛ 1.65 kips ⎞
2⎜ ⎟
⎝ 386 in./s2 ⎠
= 571 kip-in.
To limit the maximum displacement (determined previously) to comply with the support rotation criteria, the SDOF yield
force is:
WP,SDOF
R yield,SDOF =
Δ max
571 kip-in.
=
26.0 in.
= 22.0 kips
R yield,SDOF
Ryield =
KL
22.0 kips
=
0.50
= 44.0 kips
The maximum moment can be found using the equation for the maximum resistance given in Table 4-1:
Assuming the full plastic moment is developed, the minimum plastic section modulus can be determined from Equation 6-13,
with fds = 1.30Fy, as follows:
Mp
Z min =
fds
(138 kip-ft )(12 in./ft )
=
1.30 ( 50 ksi )
= 25.5 in.3
Based on this preliminary design, there is no MC8 strong enough to support the blast. There are several possible modifications to
improve the behavior of the system: increase the excited mass, increase the strength-stiffness of the system, or decrease the blast
load by integrating a variable blast pressure that is a function of the distance to the charge along the girt. In this example, increase
the steel section to an ASTM A992 W8×28.
For the W8×28, from AISC Manual Table 1-1, the geometric properties are:
d = 8.06 in.
tw = 0.285 in.
Zx = 27.2 in.
Ix = 98.0 in.4
The plastic moment is:
M p = fds Z (6-13)
=
(
1.30 ( 50 ksi ) 27.2 in.3 )
12 in./ft
= 147 kip-ft
From AISC Manual Table 3-23, using the equation for maximum deflection for a uniformly distributed load on a simply sup-
ported beam, the elastic displacement is:
2
5 Mpl L
Δ el =
48 EI
⎡ ⎤
5 ⎢ (147 kip-ft )(12 in. ft )( 25.0 ft ) (12 in. ft )
2 2
= ⎥
48 ⎢ ( 29, 000 ksi ) (98.0 in.4 ) ⎥
⎢⎣ ⎥⎦
= 5.82 in.
The strength parameters to use in the dynamic calculation are determined in the following. From Table 4-1, the maximum resis-
tance is:
Using the load factor, KL, from Table 4-1, the SDOF yield force is:
R yield,SDOF = KL R yield
= 0.50 ( 47.0 kips )
= 23.5 kips
R yield
K=
Δ el
47.0 kips
=
5.82 in.
= 8.08 kip/in.
Using the load factor, KL, from Table 4-1, the SDOF stiffness is:
K SDOF = KL K
= 0.50 ( 8.08 kip/in. )
= 4.04 kip/in.
The mass and load do not change from the previous calculations.
The period of the system, based on the SDOF mass determined previously, is:
m SDOF
T = 2π
K SDOF
1.65 kips
= 2π
(386 in./s2 ) ( 4.04 kip/in.)
= 0.204 s
The structural period (0.204 s) is more than 10 times longer than the load duration (0.00619 s), hence the assumption of impulsive
load is correct.
Figure 6-13 shows the displacement computed using SDOF software. As can be seen from this plot, no damping was introduced
into the system. The maximum deflection is 28.4 in. > 26.0 in., hence the rotation criterion is not met by a slight margin. The
elastic displacement for this beam is 5.82 in., therefore the ductility ratio for this system is:
28.4 in.
μ=
5.82 in.
= 4.88 < 20
where
F = load applied at the time of maximum response
The dynamic available shear strength, from Equation 6-11, is:
ϕVn = ϕfdvAw
where
fdv = 0.55fds (6-4)
= 0.55(SIF)(DIF)fds
Aw = dtw
Fig. 6-13. SDOF displacement for façade girt in Example 6.2— determined by computer modeling.
ϕVn = ϕ fdv Aw
= 1.00 ( 0.55 )(1.10 )(1.19 )( 50 ksi )( 8.06 in. ) ( 0.285 in.)
= 82.7 kips > 18.2 kips
Hence, the section can support the shear. The connection should be designed for this capacity.
From Example 2.1, the blast lasts for a duration of 6.19 ms and peaks at a pressure of 79.5 psi. The load associated with this blast
pressure spread about the tributary area of the column is given as:
Following a similar procedure to the girt design example previously solved, convert the MDOF system to an equivalent SDOF
system to solve. Based on the SDOF simplification used for mass and loads uniformly distributed in the plastic range, the load
and the stiffness are multiplied by the load factor, KL = 0.50, and the mass is multiplied by the mass factor, KM = 0.33, as found
in Table 4-1. Therefore, the load and mass parameters used in the discrete system are:
The equivalent impulse due to the 6.19 ms blast in this SDOF system is determined from Equation 5-12 as follows:
Fpeak,SDOF td
I SDOF =
2
( 2,150 kips )( 6.19 × 10− 3 s)
=
2
= 6.65 kip-s
The total energy produced by the blast load in the SDOF system is determined from Equation 5-3:
2
I SDOF
WP,SDOF =
2 mSDOF
( 6.65 kip-s )2
=
⎛ 4.95 kips ⎞
2⎜ ⎟
⎝ 386 in./s2 ⎠
= 1, 720 kip-in.
For the W12×53, assuming a uniformly distributed load, the following properties define the structural behavior:
Mp = fds Z (6-13)
=
(
1.3 ( 50 ksi ) 77.9 in.3 )
12 in./ft
= 422 kip-ft
The strength parameters to use in the dynamic calculation are determined in the following. From Table 4-1, the maximum resis-
tance is:
8 M pl
R yield =
L
8 ( 422 kip-ft )
=
15.0 ft
= 225 kips
Using the load factor, KL, from Table 4-1, the SDOF yield force is:
R yield,SDOF = K L R yield
= 0.50 ( 225 kips )
= 113 kips
Ryield
K=
Δ el
225 kips
=
1.39 in.
= 162 kip/in.
Using the load factor, KL, from Table 4-1, the SDOF stiffness is:
K SDOF = K L K
= 0.50 (162 kip/in. )
= 81.0 kip/in.
mSDOF
T = 2π
K SDOF
4.95 kips
= 2π
(386 in./s2 ) (81.0 kip/s )
= 0.0791 s
R 2yield,SDOF
WS,el,max =
2 K SDOF
(113 kips )2
=
2 ( 81.0 kip/in. )
= 78.8 kip-in. < 1, 720 kip-in.
The maximum elastic energy is smaller than the energy induced by the impulse. Therefore, this element will achieve plastic
behavior. As stated earlier, the intent is for columns to remain elastic. For the tributary blast load directly applied to the column,
there is not an economical solution for this column to remain elastic. But the preliminary assumption of rigid behavior of the girt
is highly conservative; the maximum load that this element is carrying comes from the reaction in the girt, not the blast pressure
on the tributary area of the column. Assuming this reaction is static, the system to solve is defined in Figure 6-15, where the 23.5-
kip end reactions from the girt are shown, calculated as half of Ryield. Note that this is different than the dynamic end reaction of
18.2 kips. The 23.5-kip load is used for redundancy and to be conservative.
The maximum bending moment for this configuration is:
The factored or required flexural strength is 1.0MB = 1.0(2,820 kip-in.) = 2,820 kip-in.
The axial load on the column, based on a roof dead load of 30 psf, is:
⎛ 50 ft ⎞ ⎛ 70 ft ⎞
PD = ( 30.0 psf ) ⎜ ⎟⎜ ⎟ 1, 000 lb/kip
⎝ 2 ⎠⎝ 4 ⎠
= 13.1 kipss
The factored or required compressive strength is 1.0PD = 1.0(13.1 kips) = 13.1 kips based on Equation 6-7.
Next, the axial compressive and flexural available strength of this element are determined and checked using the combined com-
pression and flexure interaction equations. For buckling about the weak axis, the column is assumed not to be braced by the girts.
The buckling length is Ky L = 1.0(15.0 ft) = 15.0 ft, therefore, the slenderness is:
For buckling about the strong axis, the column is also unbraced. The buckling length is KxL = 1.0(15.0 ft) = 15.0 ft; therefore,
the slenderness is:
2π 2E
Cc =
fds
This compares to the AISC Specification limit for inelastic behavior of:
KL E
< 4.71
r Fy
π2E
Fe = (Spec. Eq. E3-4)
( KL / r )2
π 2 ( 29, 000 ksi )
=
(72.6 in.) 2
= 54.3 kssi
⎛ Fy ⎞
⎜ F ⎟ (Spec. Eq. E3-2)
Fcr = ⎜0.658 e ⎟ Fy
⎝ ⎠
From Equation 6-5, Fy = fds = 1.30(50 ksi) = 65 ksi, and the critical stress is:
⎛ 65 ksi ⎞
⎜
Fcr = 0.658 54.3 ksi ⎟ ( 65 ksi )
⎜ ⎟
⎝ ⎠
= 39.4 ksi
Determine the available flexural strength, assuming that this beam is laterally unsupported because the girts cannot carry any
axial load after the blast. Therefore, the unbraced flexural strength is determined as follows from AISC Specification Section F2,
with Fy = fds :
The unbraced length of the column, Lb = 15 ft, falls between Lp and Lr. Therefore, the available flexural strength is given by:
⎡ ⎛ Lb − Lp ⎞⎤
⎢⎣
(
ϕ Mn = ϕ Cb ⎢ Mp − Mp − 0.7 fds Sx ⎜ ) ⎟ ⎥ ≤ ϕMp
⎜ Lr − L p ⎟⎥
(from Spec. Eq. F2-2)
⎝ ⎠⎦
⎧⎪ ⎞ ⎤ ⎛ 15.0 ft − 7.68 ft ⎞ ⎪⎫
⎡
(
⎛
= 1.00 (1.0 ) ⎨422 kip-ft − ⎢ 422 kip-ft − 0.7 ( 65 ksi ) 70.6 in.3 ⎜
⎩⎪ ⎣
1
) ⎟⎥ ⎜ ⎟⎬
⎝ 12 in./ft ⎠ ⎦ ⎝ 23.2 ft − 7.68 ft ⎠ ⎭⎪
≤ 1.00 ( 422 kip-ft )
= 349 kip-ft ≤ 422 kip-ft
We can compare the moment from the girt reactions to the maximum moment if subjected to full blast pressure as an indication
of the energy dissipated by the façade. The maximum moment in the column due to the girt reactions was found above to be 235
kip-ft. Determine the moment due to the full blast pressure, as follows. The natural period of the column is:
me
T = 2π (5-10)
K
15.0 kips
= 2π
(386 in./s )162 kip/in.
2
= 0.0973 s
π
k DRF = (5-15)
T td
π
=
0.0973 s 6.19 × 10 −3 in.
= 0.200
(
⎡ ( 79.5 psi ) 144 in.2 /ft 2
= 0.200 ⎢
) ⎤⎥ ( 25 ft )
⎢ 1, 000 lb/kip ⎥
⎢⎣ ⎥⎦
= 57.2 kip/ft
qeq L2
Mmax =
8
( 57.2 kip/ft )(15.0 ft )2
=
8
= 1,610 kip-ft
⎛ 235 kip-ft ⎞
⎜ ⎟ = 14.6%
⎝ 1, 610 kip-ft ⎠
Example 6.3—Design of Structural Elements Subject to Direct Blast Loading: Composite Roof Beam
Given:
The composite beam shown in Figure 6-16 will be designed for the roof blast load defined in Figure 6-17. Note that this analy-
sis is not based on the previous building. The structure consists of a 25-ft-long composite beam, consisting of a W14×22 with
a 5½-in.-thick slab with a 3-in. metal deck. The beams are spaced at 6 ft on-center. The steel material is ASTM A992 and the
concrete is normal weight with a specified compressive strength of 3 ksi. For blast design, the strength of concrete in compres-
sion is multiplied by 1.12 (dynamic increase factor from UFC 3-340-02). The composite beam is assumed to be fully braced at
both flanges.
As shown in Figure 6-16, the effective flange width is 6 ft. There is 2.5 in. of concrete above the metal deck. The center of this
concrete area is 11 in. from the center of the steel beam.
Solution:
From AISC Manual Table 2-4, the material properties are:
ASTM A992
Fy = 50 ksi
Fu = 65 ksi
Ec = 57, 000 f cʹ
⎛ 1 ksi ⎞
= 57, 000 3, 000 psi ⎜ ⎟
⎝ 1, 000 psi ⎠
= 3,120 ksi
Therefore:
Es
n=
Ec
29, 000 ksi
=
3,120 ksi
= 9.29
( 6.00 ft )(12 in./ft )
From this, the 6 ft width of the composite section becomes = 7.75 in. in the transformed uniform section.
9.29
The elastic neutral axis of the composite section is found by taking the first moments of area about the top of the concrete deck
as follows:
= 4.01 in.
where y is the distance to the elastic neutral axis from the top of the concrete deck. Then, taking the second moments of area
about the elastic neutral axis:
( ⎛
)
I steel = 199 in.3 + 6.49 in.2 ⎜ 11.0 in. +
⎝
2.50 in.
2
⎞
− 4.01 in. ⎟
⎠
= 640 in.4
h E
= 53.3 ≤ 3.76
tw fds
E 29, 000 ksi
3.76 = 3.76
fds 1.3 ( 50 ksi )
= 79.4
Therefore, because h/tw = 53.3 < 79.4, the available flexural strength is determined from the plastic stress distribution on the
composite section, and Mn = Mp. To find Mp, the concrete compression area is first found by force balance using the dynamic
strength of concrete, f ′dc = 1.12(3.00 ksi) = 3.36 ksi and fds = 1.30(50 ksi) = 65 ksi from Equation 6-5. The depth of the compres-
sion block is:
As fds
a=
0.85 fdcʹ beff
=
(6.49 in. ) (65 ksi )
2
Therefore:
⎛d a⎞
Mp = As fds ⎜ + t − ⎟
⎝ 2 2⎠
(
= 6.49 in.2 ) (65 ksi ) ⎛⎜⎝ 13.72 in. + 5.50 in.. − 2.052 in. ⎞⎟⎠ ⎛⎜⎝ 121 ftin. ⎞⎟⎠
= 398 kip-ft
ϕMp,steel = ϕ Z fdy
( )
⎡ 33.2 in.3 (1.30 )( 50 ksi ) ⎤
= 1.00 ⎢ ⎥
⎢ 12 in./ft ⎥
⎢⎣ ⎥⎦
= 180 kip-ft
The peak load for the 12.2 psi peak blast pressure is:
Based on the SDOF simplification used for mass and loads uniformly distributed in the plastic range, the load and the stiffness
are multiplied by the load factor, KL = 0.50, and the mass is multiplied by the mass factor, KM = 0.33, as found in Table 4-1.
Therefore, the parameters used in the discrete system are determined as defined in Figure 6-9. The equivalent load is:
From AISC Manual Table 3-23, the maximum elastic deflection for a uniformly distributed load is:
5 M p L2
Δ el =
48 EI tr
5 ( 398 kip-ft )(12 in./ft )( 25 ft ) (12 in./ft )
2 2
=
48 ( 29, 000 ksi )( 798 in.4 )
= 1.93 in.
From Table 4-1, the maximum resistance to use in the dynamic calculation is:
8Mp
R yield =
L
8 ( 398 kip-ft )
=
25 ft
= 127 kips
For the design of this beam, the existing dead load is applied simultaneously with the blast load and will reduce the beam strength.
From Table 4-1 and Figure 6-9, for the SDOF system, the equivalent load is:
R yield
K=
Δ el
127 kips
=
1.93 in.
= 65.8 kip/in.
mSDOF
T = 2π
K SDOF
wSDOF
= 2π
gK SDOF
2.97 kips
= 2π
(386 in./s2 ) ( 32.9 kiip/in.)
= 0.0961 s
For this example, the beam period is less than 10 times the load duration, 0.0961 s < 10(0.0144 s) = 0.144 s; therefore, the impulse
formulation cannot be used. For the rebound in the upward direction, where Mp,steel was determined previously, the force is deter-
mined from the expression in Table 4-1:
8 M p,steel
R yield,steel =
L
8 (180 kip-ft )
=
25 ft
= 57.6 kips
Again, the dead load is applied simultaneously with the blast load and therefore, for the rebound, the available strength is
increased.
R Rebound,SDOF
K Rebound,SDOF =
Δ el,steel
33.3 kips
=
3.51 in.
= 9.49 kip/in.
m SDOF
T = 2π
KRebound,SDOF
2.97 kips
= 2π
(386 in./s2 ) ( 9.49 kip/in. )
= 0.179 s
For the rebound, the beam period is greater than 10 times the load duration, 0.179 s > 10(0.0144 s) = 0.144 s; therefore, the
impulse formulation can be used for the rebound response.
Graphical Solution
This system can be solved graphically using the approach presented in Chapter 4. Using Figure 4-4, the dynamic amplification
factor is obtained based on the ratio between the structural period and the load duration:
T 0.0961 s
=
td 0.0144 s
= 6.67
From Figure 4-4, the dynamic load factor, DLF, is found to be 0.4. Dynamic load applied to the beam is:
This is less than the strength of the beam, Ryield,reduced = 118 kips, hence the system remains elastic in the first cycle.
⎢ ⎥
⎣ 2 ⎦
=
⎛ 2.97 kips ⎞
2⎜ ⎟
⎝ 386 in./s2 ⎠
= 58.7 kip-in.
The plastic displacement due to the rebound is determined as follows, where WRebound,el,max is determined from Equation 5-8:
Δ Rebound,pl =
(WImpulse,SDOF + WRebound,el,max ) (from Eq. 5-9)
R Rebound,SDOF
3.52 in.
μ=
3.51 in.
= 1.00 < 20
which is acceptable.
Computer Calculations
As an alternative, a software program may be used to evaluate the composite beam for blast loading, using an SDOF simplifica-
tion. The parameters for the load used in the analysis were previously calculated:
The load applied follows the time history shown in Figure 6-17. The yield force and stiffness were obtained previously as:
For this particular example, SAP2000 was used. However, any of the software packages or programming tools mentioned in
Chapter 4 would suffice. To solve, an SDOF model was constructed of two tension-only elements with the appropriate stiffness,
as shown in Figure 6-18. Plastic hinges according to FEMA 356 (FEMA, 2000b) are introduced in both elements to account for
any plastic behavior. These are derived in the same manner as those used in previous examples and in the modeling of Chapter 5.
The mass and peak blast forces were both applied to the node at which the tension elements connect.
Note that these values include the deflection from the dead load.
Hence, the ductility is:
Δ max
μ=
Δ el
1.98 in.
=
1.93 in.
= 1.03
2.87 in.
μ=
3.51 in.
= 0.818
Note that, as the composite action is not utilized in the rebound, the beam is much more flexible in this direction. As a result, the
beam deflects considerably more upwards (during the rebound) than it does downwards.
For the SDOF solution, the element just starts to yield as a composite beam but not during the rebound. Figure 6-20 shows the
force resultant from the SDOF system. Note that this plot, again, includes the dead load.
The maximum shear force corresponds to the maximum end reaction. Using the equation from Table 4-1 again, with F = 0 at the
time of maximum response, we get:
The available shear strength for blast loading is determined in accordance with Section 6.3.4:
ϕVn = ϕ fdv Aw
= 1.00 ⎡⎣0.55 (1.30 )( 50 ksi ) ⎤⎦ ( 13.7 in. )( 0.230 in. )
= 113 kips > 48.3 kips
where
fdv = dynamic design stress for shear defined in Section 6.1.3
Aw = area of the web
Hence, the section can support the shear. The connection should be designed for this available strength.
Py 2mg
The system must have twice the design load to
arrest the collapse.
h 0.5L
Assume the spring designed for deflection:
G Design L h mg
240 120 k
h 120 mg
k
Assume the spring ductility:
Py
GUlt 10 G y 10 k
Solve for required strength:
Py
2mg (10 ) mg (120 mg )
Py k k
10 Py k
Py 4.6 mg
8.4 RECOMMENDATIONS There are currently two sets of building code requirements
intended to address structural integrity—the New York City
The following recommendations are grouped into prescrip- Building Code (NYCBC, 2008) and the International Build-
tive detailing, general design and analytical methods. The ing Code (ICC, 2012). The New York City Building Code
prescriptive detailing recommendations are applicable to all requires all bolted connections to have a minimum of two
of the members of the structure and are intended to increase bolts and that bolted connections of all columns, beams,
the toughness and redundancy without additional analysis. braces, and other structural elements that are part of the lat-
This is similar in concept to detailing requirements in high eral load resisting system be designed as bearing connections
seismic zones. These requirements are applied to all designs, with pretensioned bolts or as slip-critical connections. All
without regard to the actual forces. The benefit derived from end connections of beams and girders must have a minimum
prescriptive detailing has not been quantified. Future test- available axial tensile strength equal to the larger of the pro-
ing could quantify the benefit of these prescriptive detailing vided vertical shear strength of the connections at either end,
requirements. The general design recommendations should but not less than 10 kips. Elements and their connections that
apply to the building as a whole and should make the engi- brace compression elements should have an available axial
neer aware of the significance of bay size, key elements, tensile strength of at least 2% of the required strength of the
column location, beam size, and metal deck slabs. The ana- compression element being braced but not less than 10 kips.
lytical design and analysis recommendations are applicable
Given:
Figure 8-11 gives the geometry of a floor system with a 24-ft by 24-ft bay size (Model A). The composite floor consists of a
52-in. composite slab consisting of 3-in. metal deck and 22-in. normal weight concrete. The wire mesh is 6×6 W1.4×W1.4
WWF (Fy = 60 ksi) and the metal deck is 20 gage (Fy = 33 ksi) with the ribs perpendicular to the beams. The composite beams
are ASTM A992 W12×19 spaced 8 ft on center, with (24) w-in.-diameter steel headed stud anchors. The beam shear connections
consist of a w-in.-thick ASTM A36 single plate with three w-in.-diameter ASTM A325 bolts. The composite girders are ASTM
A992 W16×31, with (27) w-in.-diameter steel headed stud anchors. The girder shear connections consist of a 4-in.-thick ASTM
A36 single plate with four w-in.-diameter ASTM A325 bolts. The specified compressive strength of the concrete is 3.5 ksi, and
the modulus of elasticity is:
( ) (33 )
1.5
Ec = 145 lb/ft 3 3, 500 psi 1, 000 lb/kip
= 3, 410 ksi
The floor framing and slab are continuous on all sides of the framing plan shown. The floor system will be analyzed for the
removal of the center column. At collapse initiation, the uniform service load is 87.5 psf (dead load = 50 psf, sustained dead
load = 25 psf, 25% of live load = 12.5 psf).
For the 36-ft by 36-ft floor system (Model B), the geometry, distributed loading (87.5 psf), and materials are equivalent to those
of the 24-ft by 24-ft model; however, the bays are 36-ft square, with 12-ft spacing between beams. There is also more reinforce-
ment and the girders, beams shear studs, and connections are slightly more robust. The wire mesh is 6×6 W2.9×W2.9 WWF
(Fy = 60 ksi). The composite beams are ASTM A992 W12×50 spaced 12 ft on center, with (36) 1-in.-diameter steel headed stud
anchors. The beam shear connections consist of a 4-in.-thick ASTM A36 single plate with (3) w-in.-diameter ASTM A325 bolts.
The composite girders are ASTM A992 W18×97, with (42) 1-in.-diameter steel headed stud anchors. The girder shear connec-
tions consist of a 4-in.-thick ASTM A36 single plate with four w-in.-diameter ASTM A325 bolts.
60 ksi
ε yield =
29, 000 ksi
= 0.00207
Solution:
From AISC Steel Construction Manual (AISC, 2011a), Table 2-4, the material properties are:
ASTM A36
Fy = 36 ksi
Fu = 58 ksi
ASTM A992
Fy = 50 ksi
Fu = 65 ksi
From AISC Manual Table 1-1, the geometric properties are:
W12×19
Ag = 5.57 in.2
W16×31
Ag = 9.13 in.2
W12×50
Ag = 14.6 in.2
W18×97
Ag = 28.5 in.2
Nonlinear geometry and nonlinear material properties are necessary to capture load redistribution for the vertical collapse pro-
gression scenarios. To look at the axial forces developed through catenary action in the beams, large displacements must be
considered. Nonlinear material properties are included for all failure scenarios. Discrete plastic hinges are used in the steel frame
elements to model these effects based on FEMA 356 (FEMA, 2000b). The axial hinges are simply a modeling concept that allows
the beams to yield axially. For this example, the only hinge modeled is the axial tensile yielding capacity of the beams. Figure
8-12 shows this generic hinge property. From AISC Specification Section D2(a), with ϕ = 1.00, the available axial tensile yield-
ing strength is:
For Model A
W12×19
φ Pn = φ Fy Ag
For Model B
W12×50
ϕ Pn = ϕ Fy Ag
(
= 1.00 ( 50 ksi ) 28.5 in.2 )
= 1,430 kips
Procedure
The procedure to obtain the pushover curve for this structure is based on the following assumptions:
• The structure is symmetric; no horizontal displacement is expected at the center node.
• The beams act as truss elements; no bending deflection is assumed. Therefore, the deflection between the girder and the
beam is directly related.
• The material behavior of the structure is elastic perfectly-plastic.
⎛Δ⎞
θ = tan − 1 ⎜ ⎟ (8-9)
⎝L⎠
Δ
L +ΔL = (8-10)
sin θ
EA
N= Δ L ≤ AFy (8-12)
L
6. And the vertical force associated with the vertical deflection for only one beam is:
F = N sin θ (8-13)
The pushover curve for one beam can be plotted as F versus Δ. For more beams framing into the same node, the forces for each
one can be added. For beams framing into the girder, the vertical deflection is associated with the control displacement, Δ. For
this problem there are two beams framing into the girder and the vertical displacement associated with each beam is qΔ and 3Δ
as shown in Figure 8-16. The sum of the vertical forces from the two beams framing into the girder and the vertical force from
the girder give the total vertical force. This approach assumes that the deflections along the girder vary linearly with the distance
to the point of maximum deflection. If the bending deformation of the girder is included, the solution is more accurate. When
bending is not considered, as in this example, the result is conservative.
Capacity Curve
The theory presented in the previous section and Section 8.3.2 is used to compute the pushover and capacity curves for both
models. Figure 8-17 and Figure 8-19 show the pushover curve and the capacity curve for Model A and Model B, respectively,
and Figure 8-18 and Figure 8-20 show the tensile force for Model A and Model B, respectively. The area below the pushover
curve is the energy that the structure can absorb. The area below the pushover curve divided by the corresponding displacement
yields the capacity curve of the structure.
Fig. 8-17. Pushover curve, capacity curve and load for Model A.
Parametric Study
A parametric study of Model A is presented in this section. The area of the elements is varied from 0.5 to 2 times the initial
area. The element capacities remain proportional to the element areas. All other parameters in the system remain constant. The
pushover and capacity curves for each case have been determined. The results are shown in Figure 8-21. Note that the behavior
of this system is nonlinear for each load due to geometric and material nonlinearities. Therefore, twice the deflection does not
result in twice the force.
As beam size is increased, vertical deflection decreases. The whole assembly has to absorb the same amount of energy, but with
a larger beam the connections must be able to carry more of the load because the beam absorbs less energy in bending. As beam
size decreases, vertical deflection increases and thereby the axial force in the beams and connections is reduced. However, the
larger the deflection the more longitudinal and rotational ductility must be absorbed by the beams and connections.
Table 8-1 shows the axial load in the different elements in Model A. Note that the stiffer the system, the bigger the axial load
needed to be supported by the connection. These values are higher than the axial capacity from a regular shear connection which
would require significant strengthening. A more economical approach may be to account for the tensile strength developed in the
slab reinforcing, which will be checked in part (b) of this example. Based on this parametric study, elements with greater axial
stiffness result in smaller vertical deflections, greater axial force, and smaller ductility demand in the connection. The resulting
connection is strong and expensive. To reduce the cost of the connection, the engineer should consider using as small a member
as possible.
Fig. 8-19. Pushover curve, capacity curve and load for Model B.
Figure 8-22 presents the vertical deflection versus time results from the time-history analysis. For Model A, the maximum
deflection is 22.5 in. Comparing the nonlinear dynamic approach (time-history method) to the nonlinear static approach (energy
balance/steel pushover method), it is seen that both methods produce similar results. The energy balance method gives a maxi-
mum deflection of 22.0 in. (Figure 8-17). Based on the time-history analysis, the W16×31 and W12×19-1 framed to the removed
column have yielded (Figure 8-23). Beam W12×19-2 has a tensile force of 227 kips and beam W12×19-3 has a tensile force of
59 kips (Figure 8-23). These axial forces are slightly higher than those calculated by way of the energy balance method (Figure
8-18). It is important to note that the nonlinear dynamic approach is substantially more time consuming and complex than the
simplified pushover method. In both models the axial forces in the beams and girders are substantial at large deflections. Connec-
tions sufficient to develop these axial forces, while exhibiting the necessary rotational ductility, would be necessary.
Compression Ring
The boundary conditions modeled in these problems assume pinned connections at the edge of the collapsed span. This assump-
tion requires that the reaction be distributed to the slab using continuity in the reinforcement or compression in the slab (Figure
8-24). For composite floor systems with sufficient horizontal restraint (i.e., interior bays) this is a reasonable assumption. Based
on testing by Allam et al. (2000), simply supported slabs will self-equilibrate at large deflections through tension in the central
regions and a compression ring along the perimeter zones. The amount of catenary action to resist collapse is dependent on the
stiffness of the anchorage. The following examples will verify that the compression ring induced in the slab provides sufficient
anchorage. For a more accurate assessment, a finite element analysis or other comparable analysis could be used.
FComp
w=
0.85 fconc t
323 kips
=
0.85(3.5 ksi)(2.50 in.)
= 43.4 in.
The compressive force causes axial deformation in the compression ring which results in a radial displacement at the perimeter
beam-column joints (Figure 8-26):
To check the slippage of the steel headed stud anchors and the capacity of each stud, refer to Grant et al. (1977) and Easterling
et al. (1993). According to Easterling et al. (1993), the nominal strength of each steel headed stud anchor is 24 kips and the slip-
page is 8 in. To develop the full strength of the beam, 19 shear studs are required. Because the design is conservatively using
more than one steel headed stud anchor per foot (27 studs on the girder and 24 on the beam), the actual slip should be less than
that calculated.
The total displacement due to the flexibility of the compression ring (0.224 in.) plus the slippage of the steel headed stud anchors
(8 in.) is 0.349 in. Using the procedure discussed previously, the extended length of the beam in Model A due to the central
vertical deflection of 22.0 in. is:
Adding the 0.349 in. to the elongated length as shown in Figure 8-27, the total vertical deflection is:
The resulting vertical deflection, ΔTOT, is 20% greater than the calculated deflection, Δ, when complete anchorage was assumed.
This derivation is conservative and the additional displacement is within an acceptable range.
Fig. 8-27. Additional deflection due to anchorage flexibility and stud slippage.
Conclusions
In this part of Example 8.1, the energy method was used to develop the pushover and capacity curves. The results were compared
to the results of a time-history analysis and found to be similar. This example has explored the concept of catenary action of the
steel beams and girders and has determined that the axial demand on the connections, in many instances, requires the develop-
ment of the axial yield strength of the beams and girders. Without test data or complex analytical models, it is difficult to ascertain
if the connection design is adequate. By design, connections that achieve the tension capacity of the steel member, while main-
taining a level of longitudinal and rotational ductility, become ductile moment connections that resist bending behavior in addi-
tion to tensile forces. Increasing the shear tab thickness and/or directly welding the beam web to the column will also enhance
the connection strength and ductility as shown by Khandelwal and El-Tawil (2007). Through a parametric study, the idea was
introduced that smaller steel members increase the deflections and decrease the demands on the connections.
In addition, the idea was introduced that it is superior to rely on the concrete slab rather than to rely on the adjacent structural
bays to provide the axial restraint for the catenary forces. The slab can develop a compression ring that produces the required
beam restraint. The tension in the member can be transferred through the steel headed stud anchors to the concrete, creating a
compression ring around the perimeter of the bay that contains the removed column. In this case, the connections do not develop
the full tensile capacity of the member. However, because the transfer of the beam tensile force to the slab is critical, the spac-
ing of the steel headed stud anchors should not exceed one stud anchor per foot of beam. A simplified method for checking the
strength of the compression ring in the concrete was presented. If the concrete compression ring is not considered, perimeter
connection forces become quite large.
When designing to resist progressive collapse, it would generally not be economical for the engineer to exclusively consider
catenary action of the structural steel. In this example, cable action or catenary action in the steel beams redistributed vertical
loads and assisted in reaching a new equilibrium. However, the axial forces necessary to redistribute the vertical loads are large.
Inclusion of the flexural behavior of the steel connections, or of the composite steel-slab action at connections, will give a more
economical design. Moment connections resist forces in bending, thus generating a plastic hinge mechanism that dissipates
energy and resists the collapse of the structural bay. The analytical concepts introduced in this example are also valid for hinges
in bending. If the engineer wants to use moment connections, the procedure will be similar, with replacement of the axial hinges
with moment-rotation hinges. Further still, the engineer may take advantage of membrane action of the slab, considering the
tensile strength of the steel reinforcing within the slab and the compression capacity of the concrete in the slab to achieve a more
economical design. This will be further explored in Part b of this example.
where
w = predicted distributed load
3L x ε x
δ= (8-15)
2 sin 6ε x
The complete load versus central deflection response can be obtained by using the following solution procedure:
1. Choose a value of εx.
2. Calculate εy.
3. Determine Tx and Ty corresponding to εx and εy using the stress-strain relationship.
4. Calculate the load, w.
5. Calculate the defection, δ.
See Table 8-2 for the calculations.
Model Results
The parameters used here are those defined in Part (a) of this example. Following the procedure derived from Mitchell and Cook
(1984) from the previous section, the pushover curve and capacity curve are plotted and shown in Figure 8-28 for Model A and
Figure 8-29 for Model B. The capacity curve is obtained by the same procedure utilized in Part (a) of this example.
Considering only the strength of the steel reinforcement in the slab, the total vertical deflection is 48.5 in. for Model A (wire mesh
reinforcement 6×6-W1.4×W1.4 WWF which is the equivalent of No. 3 bars spaced at 48 in. on center) and 54.5 in. for Model B
(wire mesh reinforcement 6×6-W2.9×W2.9 WWF which is approximately the equivalent of No. 3 bars spaced at 24 in. on cen-
ter). The deflections calculated for Model A and Model B are equivalent to rotations of 9.6° and 7.2°, respectively.
Increasing the area of reinforcement reduces the deflection in the slab. For example, if the reinforcement in Model A is increased
to 6×6-W2.9×W2.9 WWF, the deflection is reduced to 26.5 in. or 5.26°. Similarly, if the reinforcement in Model B is increased
to 6×6-W5×W5 WWF, the deflection is reduced to 36.3 in. or 4.8°. Figure 8-28 and Figure 8-29 show the pushover and capacity
curves for Model A reinforced with 0.028 in.2/ft and Model B reinforced with 0.058 in.2/ft, respectively. Table 8-3 shows the
deflection and corresponding rotation for various areas of steel slab reinforcement for both models.
Conclusions
In this part of Example 8.1, it was shown that slab structures could develop secondary load carrying mechanisms and exhibit a
degree of membrane action. The deflections in Model A and Model B for this formulation are within reason when considering
only the strength of the slab reinforcement. The deflection calculations neglected any contribution of the steel in the metal deck.
If continuity of the metal deck can be guaranteed, the capacity of the system increases. Additionally, only the tension capacity
was checked because in compression there is concrete in addition to the reinforcing providing strength; thus, tension controls.
The performance of this system can be further improved by combining the capacity curve of the steel and the slab. Increasing the
slab reinforcement reduces the impact of the collapse on the beam connections and reduces the rotation and axial load demand.
In this example, the slab’s resistance to collapse was investigated with respect to the individual bay. Horizontal restraint was not
considered because the tests by Mitchell and Cook (1984) show that, with large deformations, the response of a simply supported
slab is essentially the same as a fully restrained slab. This result is due to the ability of the slab to form its own in-plane edge
Fig. 8-28. Pushover and capacity curve for the slab for Model A.
Boundary Conditions
The following boundary conditions are assumed:
• Pin supports at the base of the columns (Restraints = Ux, Uy and Uz).
• Horizontal restraints at the top of the columns (Restraints = Ux and Uy).
Fig. 8-29. Pushover and capacity curve for the slab for Model B.
Loads
The total service load of 87.5 psf is applied uniformly as point loads at the nodes of the shell elements.
Fig. 8-30. CSI Perform-Collapse model for the composite floor system.
Composite Beams
The composite beams were modeled as elastic steel beams with moment releases and nonlinear axial connections at the ends. The
beams are ASTM A992 W12×19 sections.
Composite Girders
The composite girders were modeled as elastic steel beams with moment releases and nonlinear axial connections at the ends.
The girders are ASTM A992 W16×31 sections with Fy = 50 ksi.
where the edge distance is 1.25 in. The bolt holes were assumed to be standard holes (m in.). Tearout is shown to control and
the stress-strain relationship of the single-plate tear-out was assumed as shown in Figure 8-34.
Deformed Shape
The deformed shapes for both phases are shown in Figure 8-38 and Figure 8-39. The deformation increased almost linearly with
the load during phase 1. During phase 2, the deformation at the center point, where the column was removed, kept increasing
gradually with the load to 2.7 in., and then it increased suddenly from 2.7 in. at 42% of the load to 14.8 in. at 60% of the load.
After the sudden increase, the deformation increased gradually again until it reached a maximum of 20 in. at the full static load.
The deformation continued to increase gradually until an energy balance was reached at 200% of the static load and a deforma-
tion of 33.6 in. The pattern of displacement is shown in Figure 8-40.
Fig. 8-38. Deformed shape for phase 1 (gravity load with column in place).
Fig. 8-42. Yielded steel layers at the energy balance stage (200% static load).
Fig. 8-43. Steel tension strain in composite slab: steel layers start yielding at 44% of the load.