Design of Structures For Missile Impacts
Design of Structures For Missile Impacts
Design of Structures For Missile Impacts
Revision 2
if
SEPTEMBER 1974
TOPICAL REPORT
DESIGN OF STRUCTURES
FOR MISSILE IMPACT
VIASTER
I
{I
DISCLAIMER
This report was prepared as an account of work sponsored by an
agency of the United States Government. Neither the United States
Government nor any agency Thereof, nor any of their employees,
makes any warranty, express or implied, or assumes any legal
liability or responsibility for the accuracy, completeness, or
usefulness of any information, apparatus, product, or process
disclosed, or represents that its use would not infringe privately
owned rights. Reference herein to any specific commercial product,
process, or service by trade name, trademark, manufacturer, or
otherwise does not necessarily constitute or imply its endorsement,
recommendation, or favoring by the United States Government or any
agency thereof. The views and opinions of authors expressed herein
do not necessarily state or reflect those of the United States
Government or any agency thereof.
DISCLAIMER
Portions of this document may be illegible in
electronic image products. Images are produced
from the best available original document.
TOPICAL REPORT
BC-TOP-9-A
Revision 2
PREPARED BY:
R. B. Linderman
J. V. Rotz
G. C. K. Yeh
APPROVED BY:
W. A. Brandes
H. W. Wahl
BC-TOP-9-A
Rev. 2
14
CAVEAT:
UNITED STATES
20545
Mr.
John V. Morowski
-2
NO
4 1974
of
perforation formulae assume that the missile strikes the target normal
to the surface, and the axis of the missile is assumed parallel to
the line of flight.
'he penetration is the Modified Petry equation, while that for perforation
ified
to allow its use for concrete strength other than 3000 psi
The wall
by the thicknesses
However the
governing requirement for an overall structural response design cohsideration is that the maximum deflection of the target shall be
limited so as not to impair the function of other safety related
equipment.
the target
i.e., determination
for
Sample
If properly
utilized in nuclear power plant structural design work, the procedures and criteria contained in the report should provide
b)
c)
BC-TOP-9-A
Rev. 2
ABSTRACT
This report contains methods and procedures for evaluating the effects of
missile impact on structures. A means to evaluate the change of velocity
of a missile passing through a liquid is also included. Missile impact
effects on structures are evaluated in terms of local damage (penetration,
Empirical formulae
perforation, and spalling) and structural response.
Structural dynamic principles are used
are used to evaluate local effects.
to evaluate structural response.
ACKNOWLEDGMENT
This document is
contributors.
The following is a chronological account of major participants contributing to the development of this document:
Revision 0 (issued October, 1972) was prepared by M. Fakhari,
B. Linderman, J. Rotz and M. Suarez and approved by A. J. Bingaman
(Gaithersburg Office Chief Civil Engineer) and D. W. Halligan
(Power and Industrial Division Chief Civil Engineer).
Revision 1 (issued July, 1973) was prepared by R. B. Linderman,
M. Fakhari, J. V. Rotz, E. Thomas, G. A. Tuveson, and G. C. K. Yeh;
and approved by W. A. Brandes (Los Angeles Power Division, Chief
Civil Engineer) and L. G. Hinkelman (Thermal Power Organization,
Chief Civil Engineer); Technical Consultant, N. M. Newmark.
Revision 2 (issued September, 1974) was prepared by R. B. Linderman,
J. V. Rotz and G. C. K. Yeh; and approved by W. A. Brandes (Los
Angeles Power Division, Chief Civil Engineer).and H. W. Wahl (Thermal
Power Organization, Chief Civil Engineer); Technical Consultant,
N. M. Newmark.
BC-TOP-9-A
Rev. 2
CONTENTS
Title
Section
1.
1.1
1.2
1.3
1.4
2.
2.1
2.1.1
2.1.2
2.1.3
2.2
2.3
2.3.1
2.3.2
3.
INTRODUCTION
Page
1-1
1-1
1-1
1-2
1-2
General
Approach
Missile Characteristics
Target Characteristics
2-1
LOCAL EFFECTS
Reinforced Concrete Targets
Penetration
Perforation
Spalling
Steel Targets
Multiple Element Barriers
Reinforced Concrete Barrier
2-1
2-1
2-2
2-3
2-3
2-4
2-5
2-5
Steel Barrier
STRUCTURAL RESPONSE TO MISSILE IMPACT LOAD
3-1
3.1
General
3-1
3.2
3-1
3.3
3.3.1
3.3.2
3.3.3
3.4
3.5
3.5.1
3.5.2
3.5.3
3.5.4
3-2
3-2
3-3
3-3
3-5
3-7
3-7
3-7
3-8
3-9
4.
4.1
4.2
4.3
DESIG'N GUIDELINES
4-1
4
xi
4-1
4-1
4-3
BC-TOP-9-A
Rev. 2
CONTENTS
Title
Section
5.
5.1
5.2
5.2.1
5.2.2
5.2.3
(Cont)
Page
5-1
SPECIAL PROBLEMS
Force-Time History for Automobile Crash
Penetration of a Missile Through a Liquid
Liquid Depth is Less Than or Equal to Missile
Length
Liquid Depth is Greater Than Missile Length (H > L)
Definitions of Notations
5-1
5-1
5-2
5-2
5-3
0
xii
BC-TOP-9-A
Rev. 2
LIST OF APPENDICES
Appendix
Title
Page
APPENDIX A
A-I
APPENDIX B
Glossary
B-i
APPENDIX C
C-I
APPENDIX D
Derivations
D-1
APPENDIX E
Sample Applications
E-1
APPENDIX F
F-I
4
xiii
BC-TOP-9-A
Rev. 2
LIST OF TABLES
Table
Title
Page
4-1
4-2
4-4
Beams
4-3
4-5
4-6
4-4
4-7
5-1
C-1
C-2
Incompressible Flow
Concrete Penetration,
Spalling Formulas
Perforation in
5-4
Perforation,
Steel Formulas
i
xv
and
C-5
C-8
BC-TOP-9-A
Rev. 2
LIST OF ILLUSTRATIONS
Figure
2-1
2-2
(K p)
2-8
2-9
3-12
3-2
Energy-Displacement Functions -
3-13
3-3
Energy-Displacement
With Other Loads
Impact Combined
2-4
4-1
Coefficients
Functions -
3-14
Sections
2-6
3-1
21
for
2-7
21
Page
2-3
Title
4-8
5-1
Penetration of a Missile in
C-I
C-2
thru
C-22
Deleted
xvii
a Liquid
5-5
C-9
BC-TOP-9-A
Rev.
Section I
INTRODUCTION
1.1
GENERAL
response.
Missiles may be generated by an event that is
(such as com-
1.2
APPROACH
Verify that the maximum deflection does not impair the function
of other safety related systems.
1-1
BC-TOP-9-A
Rev. 2
1.3
MISSILE CHARACTERISTICS
1.4
TARGET CHARACTERISTICS
1-2
BC-TOP-9-A
Rev. 2
Section 2
LOCAL EFFECTS
2.1
2.1.1
PENETRATION
X = 12 K
Logl
215,000
(2-1)
where
X = Depth of missile penetration into concrete element of infinite
thickness (inches)
Note: Usually this equation expresses the depth of penetration in feet; however, for this document it has been
modified to express it
K
A
p
p
Vs
= Penetration coefficient
in inches.
for reinforced concrete
Missile weight
Projected frontal area of missile
This formula is
(psf)
When the element has a finite thickness the depth of penetration is:
X L +e
X,
2-1
(t> 2X)
(2-2)
BC-TOP-9-A
Rev. 2
where
XI = Depth of penetration of missile into a concrete element of finite
thickness (inches).
e = Base of Napierian Logarithms
t = Thickness of concrete element
2 1
Penetrations
(inches)
examples of missiles are shown in
PERFORATION
W8 (UVs
D1.
\'331(
)(2-3)
where
T = Thickness of concrete element to be just perforated
W = Weight of missiles
(inches)
(lb)
D = Diameter of missiles
(inches)
Note:
For irregularly shaped missiles, an equivalent
diameter is used.
The equivalent diameter is taken as
the diameter of a circle with an area equal to the circumscribed contact, or projected frontal area, of the
non-cylindrical missile.
V
c
Thi!. formula is
BRL,
formula.
t p = 1.25T ; T + 10 (in
2-2
inches)
(2-4)
BC-TOP-9-A
Rev. 2
2.1.3
12
SPALLING
Spalling
ment may
estimate
that the
of concrete from the side opposite the contact surface of the eleoccur even if the missile will not perforate the element.
For an
of the thickness that will just start spalling, it is recommended
following equation be used:
Ts = 2T
(2-5)
where
T
2.2
1.25 Ts
T 5 + 10(in inches)
(2-6)
STEEL TARGETS
Steel targets, such as pipes and mechanical equipment vessels, may be perforated by a missile.
Sometimes, protruding elements of a missile may
puncture a steel target when the entire missile does not perforate or pass
through the target.
The minimum contact area of a missile protrusion is
used to calculate puncture thickness and the projected area of the entire
missile is used to calculate perforation thickness.
The BRL Formula is shown below, modified by setting a material constant
K = 1 and solving directly for steel plate thickness, T, which will just
be perforated by the missile,
MVs2/
T=
L22
672D
2-3
(2-7)
BC-TOP-9-A
Rev. 2
where
T = Steel plate thickness to just perforate (inches).
M = Mass of the Missile (lb sec 2 /ft)
W = Weight of the Missile (lb)
Vs = Striking Velocity of the Missile Normal to Target Surface (ft/sec)
D = Diameter of the Missile (in.)
Note: For irregularly shaped missiles the
diameter is used.
The equivalent diameter
diameter of a circle with an area equal to
scribed contact, or projected frontal area
non-cylindrical missile.
equivalent
is taken as the
the circumof the
2.3
(2-8)
It may be desirable to construct a missile barrier of several thinner elements, instead of one thick element.
Analysis of missile barriers composed
of several elements involves determining the residual velocity (Vr) after
perforation of one element and using this value for the striking velocity
(Vs) on the next element.
The following formula is used to determine the
residual velocity, Vr (see Appendix C)
r
= 2 _2 ) 1/2
V-V
For((Vs
V)
For (V
V)
(2-9)
Vr =0
where
21
2-4
BC-TOP-9-A
Rev. 2
is
2.3.1
106]
Vr
(2-10)
STEEL BARRIER
V=
1.5]l/
[r
(inches)
4
2-5
(2-11)
BC-TOP-9-A
Rev. 2
0
0.0050
z
UU-
uJ
0
U
0
Ul
z~l
u.
a-
7,000
28 - DAY COMPRESSIVE STRENGTH OF CONCRETE
Figure 2-1
VALUES OF PENETRATION COEFFICIENT (K ) FOR REINFORCED CONCRETE
(Reference 14)
2-6
0I
24
W=12
22
Uj
A=l10
W
A
COMPRESSIVE.STRENGTH OF
CONCRETE = 3000 PSI
u)
20
z
x
Lu
NOTE:
18
I--
LU
oil
1. Win pounds
2. Ain inches 2
3. Equation 2-1 uses
2
Ap in lb/ft
010_Z_
16
0
U
C.,
14
ww
be
I-J
12
Aw
LL
10i
w
z
Z
ZIN
5<00
I-a.
LU
0
z
0
w/t
A
w--=
A
LU
zLU
4..
.
...
W=
___
a.
'0
50
io6
300
200
400
5(
30
I-t
28
I-
LL.
IL
COMPRESSIVE STRENGTH OF
CONCRETE = 3000 PSI
26
= 12 -
__
I
0
24
1in04
D_____
poundsD_
NOTE: W in ihW
22
D in inches
20
Dr
18
w =Z
/006
1000
16
LU
0.
D' 8
14
I--
12
10
b-
D1,.
-4W0
8
LU
Dwa
0)
2.
Wl-
4
2
0
50
100
200
300
MISSILE VELOCITY (FT/SEC)
Figure 2-3
PERFORATION OF REINFORCED CONCRETE FOR
VARIOUS MISSILES (BALLISTIC RESEARCH LAB)
400
500
0.
LEGEND:
40
20
80-
z
d
30
U.'
15
30
60-
z
-j
-j
cc
C.3
a.
U.
C2
C07
20 =
10
I.
LU
ZC.3
0
40LU
c-
02
LL
10
clO
I--
LU
(4
20-
lie
Lu
I-
0.
0
0
50 100
200
300
400
500 0
50 100
200
300
400
500
0-
wc.~
<I
0
'.0
BC-TOP-9-A
Rev. 2
Section 3
STRUCTURAL RESPONSE TO MISSILE IMPACT LOAD
3.1
GENERAL
When a missile strikes a target, large forces develop at the missile target interface, which decelerate the missile and accelerate the target.
If the interface forcing function is known, (experimentally determined),
the target structure can be modeled mathemat:Lcally and conventional numerical techniques can be used to predict structural response.
For most cases,
the forcing function is not known, and a rational method involving energy
balance techniques is used to estimate structural response.
This involves
using the strain energy of the target at maximum response to balance the
residual kinetic energy of the target (or target-missile combination)
resulting from missile impact.
For investigation purposes, it is convenient to model the event as a missile
of mass, Mm, and striking velocity, Vs, impacting a spring-backed target
mass, Me.
The spring may be linear, bilinear, or non-linear, depending on
the target structure resistance-displacement function. Since the actual
coupled mass varies during impact, an estimated average effective target
mass, Me is used to evaluate inertia effects during impact.
The impact may be either elastic or plastic, depending on whether or not
significant energy losses are sustained during impact.
These losses are
associated with inelastic deformations, local damage in the impact zone,
etc.
3.2
Since the duration of impact is very short, (usually less than a few milliseconds), the target mass displacement and the corresponding spring force
are also very small. Neglecting the spring force effect during impact,
(a slight conservatism), the velocities of the missile and target after
impact are calculated from the following relationships:
V
m
VT"
(M - eM)
M +M
m
SV M ((+e)
M+M
m
3-1
(3-2)
12
BC-TOP-9-A
Rev. 2
Vm - Missile velocity after impact
VT = Target velocity after impact
V
M = Mass of missile
m
Me =
e = Coefficient of restitution
3.3
3.3.1
ELASTIC IMPACT
Equations (3-1) and (3-2)(12)* show that the velocity of the missile after
For this
impact is opposite to that of the target if Mm is less than eMe.
case, the strain energy, Es, of the responding target spring required to
diminish the target mass velocity to zero (maximum target response) is
numerically equal to the kinetic energy of the target mass at the end of
the impact duration.
MVT2
Es
(3-3)
21
= --
the impact is
is;
2
2M M V
EEs
(M+mes M)
(m
(3-4)
Me
Referring again to equations (3-1) and (3-2), the velocity of the missile after
impact is in the same direction as that of the target if Mm is greater than
eMe.
In this case, the target spring decelerates the target mass, allowing
the missile to overtake the target, which results in multiple impact.
If
the impact is
missile through a series of impacts and absorb all the initial kinetic
*References are in appendix F.
3-2
BC-TOP-9-A
Rev. 2
M V2
E
3.3.2
(3-5)
12
PLASTIC IMPACT
reduces to zero
For a plastic collision, the coefficient of restitution
(e = 0) and the missile and target masses attain the same velocity at the
end of impact duration.
If the impact is of short duration, the target
displacement and corresponding spring force effect during impact are small,
The strain energy required to stop
and can be conservatively neglected.
the target-missile combination is then the sum of the kinetic energy of the
missile and the target masses at the end of the duration of impact.
MV22
MV~
M2
Es =
From equations
M
2 M+
eT2
(3-6)
=V
MV
m s
Mm + Me
e~
(3-7)
(3-7)
into equation
(3-6),
M2V2
2+M)e5
Es 's22M(Mmmem+sM
e
3.3.3
(3-8)
F(t)
R(x) = Mee"
F(t)
Force-time function
R(x)
K = Acceleration
x
of target mass
3-3
BC-TOP-9-A
Rev. 2
Numerical methods are usually used for solution of equation 3-9 which is
solved for the maximum value of displacement xm. The target strain energy
is then;
2
R(x)dx
Jo
Rdt
t,
is;
- R(x)I
LF(t)
(t
212
2l~,ft)]
= M
E (t)
E(t)
or
dt
Me
0o
is
then
EIt
F(t) -
R(x)]
dt 2(3-10)
Equation (3-10) shows that deletion of the R(x) term will result in a
conservative overestimate of E(t).
If R(x)<< F(t) during impact, t, the
inaccuracy is usually negligible.
For this condition, the kinetic energy
of the target mass at time tI is conservatively estimated as;
12
t 1
tl
E
fF(t)
(3-11)
ti
The applied impulse,
I,
is
dt
2e
by definition,
curve.
eti
I
F(t) dt
3-4
0@
BC-TOP-9-A
Rev.
If
2
(1/2)MeV
e T(3
I2
2M
Et
(3-12)
the velocity of the missile approaches that of the target at the end of
time, ti, equal to the duration of impact. The strain energy of the target
required to stop the missile-target combination is then;
2
MV
2 --+ 2M
Es =-
(3-13)
2
Vm =VT
and,
V2
m
12
M2
e
(Mm + Me) 1
E
3.4
e
2M2
e
(3-14)
The effective target mass during impact varies from a low value at initial
contact and generally increases to an upper limit during or at the end of
the impact duration. Due to the complex phenomenology associated with
missile impact, no general analytical solution is available to evaluate the
effective coupled mass on a continuous time basis. The average effective
mass can, however, be estimated, utilizing the results of impact tests on
reinforced concrete beams( 7 ) wherein the measured maximum structural
response was used to back-calculate the average mass during impact.
3-5
BC-TOP-9-A
Rev. 2
Based on these data, the following formulae provide a lower limit estimate
of Me (which results in an upper limit estimate of kinetic energy after
impact).
For concrete beams:
By cT
21
Me =(DxX + 2T)
fif
B _(Dy + 2T)]
[if
B >(Dy + 2T)]
ycT
Me = (D
+ 2T)(D
7'
+ 2T)
(3-15)
ycT
c
(Dy + T)
(3-16)
(Dx + 2d) Hx
(3-17)
(3-18)
D
Y
Yc
Ys
Acceleration of gravity
3-6
BC-TOP-9-A
Rev. 2
3.5
3.5.1
GENERAL PROCEDURES
The strain energy, Es, required to stop the target (or missile-target
combination), is determined from the relationships in sections 3.2 and 3.3.
The resistance-displacement function, R(x), for a concentrated load at the
area of impact is determined from the target structure physical configuration and material properties.
The estimated maximum target response is determined by equating the available target strain energy to the required strain energy and solving for the
maximum displacement x . (See Figure 3-1.)
3.5.2
kx
(3-19)
Xm =
If other loads are present on the target structure which will act concurrently
with missile impact loads, the maximum combined displacement is determined
as follows:
Then
x
- X +X'
o
Since
2
J2Es] /
Let
x' - Displacement due to missile impact (See Figure 3-1)
x
3-7
BC-TOP-9-A
Rev. 2
that
21it follows
Xm
3.5.3
(3-20)
o~k
For elasto-plastic
R(x) = kx,
(O<XxXe)
wher e
xe = Yield displacement
R = Plastic resistance.
m
Then
Es = Rm(x
or
X
E
x
_I + e
R
2
m
is
(3-21)
xm
x
e
E
xr x-R
(3-22)
If other loads are present on the target structure which will act concurrent
with missile impact loads, the maximum combined displacement is determined
as follows:
Let
xt = Xe - x
3-8
BC-TOP-9-A
Rev. 2
xe = yield displacement
Es
K--2
+ kx'
or
= x
x!
T2
-
xm
Substituting x'
(X -x)
in
X m_
5e
(x2
x2-
(3-23)
Es
=
X0o
2IX0+
(x
(3-24)
The values of Pr should be less than the allowable ductility ratios w given
in section 4.
3.5.4
Ee T
R(x) dx
3-9
(3-25)
BC-TOP-9-A
Rev. 2
When no other concurrent loads are acting, the maximum displacement occurs
at the value of x where Ee is equal to Es.
The correct value of xm is therefore the value of x, which will satisfy the following relationship:
=f
R(x) dx
shown in
(3-26)
figure 3-2.
When other loads are acting concurrent with missile impact loading,
correct value of x. will satisfy the following relationship:
Es
R(x)dx - Ro (x
x
= equivalent static
x )
the
(3-27)
is
safety factor
Fs
Fs
s
0.5 if
R(x)
is
0.25 if R(x) is
analysis)
(such as by failure
Ef =f
R(x) dx
(3-29)
xf = displacement at failure
i
3-10
BC-TOP-9-A
Rev. 2
Ef =
R(x) dx - R
and
R
Ra
r
R f = resistance at failure
3-11
(xf
x0
(3-30)
RESPONSE
RESISTANCE DISPLACEMENT
AVAILABLE STRAIN
ENERGY WITHOUT
AVAILABLE STRAIN
ENERGY WITH
FUNCTION
OTHER LOADING
OTHER LOADING
mc0
0
RI
ELASTIC
R - kx
ELASTO-PLASTIC
I-
kx
Rm
xe
xe
xe
xx
Km
x0
Km
xe
Xm
NONLINEAR
'mT
"O
xf
XQ
SHADED AREA (STRAIN ENERGY) MUST EQUAL E. IFROM SECTIONS 3.2 AND 3.3)
Figure 3-1
RESISTANCE-DISPLACEMENT FUNCTIONS WITH ASSOCIATED STRUCTURAL
RESPONSE WITH AND WITHOUT THE EFFECT OF OTHER LOADS
Km
xf
BC-TOP-9-A
Rev. 2
'U
a',
Xnm
DISPLACEMENT X
Figure 3-2
ENERGY-DISPLACEMENT FUNCTIONSIMPACT LOADS ONLY
3-13
xf
BC-TOP-9-A
Rev. 2
E@- f Rtm)dx
0
.Jca
E
mu
2
mu
Ea
xm
DISPLACEMENT X
Figure 3-3
ENERGY-DISPLACEMENT FUNCTIONS IMPACT COMBINED WITH OTHER LOADS
0
3-14
BC-TOP-9-A
Rev. 2
Section 4
DESIGN GUIDELINES
4.1
ing capacity.
The dynamic capacity of the structural elements must be
based on material dynamic strength properties which are obtained by applying
strength value:
a dynamic increase factor (DIF) to the static
fdyn= (DIF)
fStat
(4-1)
where
2
= specified
static
strength value
table 4-1.
DESIGN PARAMETERS
4.2
Ia =
.1
(Ig
+ Ic)
=11(2
d3)
(4-2)
where
1g -
4-1
BC-TOP-9-A
Rev. 2
b = width of concrete section
F = coefficient for moment of inertia of cracked section with
tension reinforcing only.
(See figure 4-1.)
t = concrete thickness
d
where
A
1.4 F ',
0.25 f'
Y
1. 4 f
f--
bd
0.25 f'
- A'
dI
4-2
(4-4a)
--
fY
BC-TOP-9-A
Rev. 2
where
ft = compression strength of concrete
c
A'
4.3
of concrete
The maximum allowable ductility ratios for concrete and steel members are
However,
4-3
BC-TOP-9-A
Rev. 2
Table 4-1
DYNAMIC INCREASE FACTOR
(DIF)
(From Ref. 19)
I.
DIF
Compression
1.25
1.0
Bond
1.0
Reinforcing Steel
Tension & Compression For 40ksi yield strength steel
1.2
1.0
1.0
It.
Structural Steel
Flexure,
1.2
1.0
2 160ksi
1.0
Shear
4-4
BC-TOP-9-A
Rev. 2
Table 4-2
RESISTANCE-YIELD DISPLACEMENT
VALUES FOR BEAMS
DESCRIPTION
YIELD
DISPLACEMENT
RESISTANCE
CANTILEVER
(il
I'LLF
3E1
12
RL 3
4E
12
SIMPLY SUPPORTED
(2)
R
Xe
JIM
Ft -
L.
FIXED SUPPORTS
IM)
4(M+ + 94)
-e
RL 3
(4)
MULTI-SPAN
2
't (M+ + Mal
__
U.
L
L
I L47 I -L47
Where
M+
u
MU
-e
001 RL
O.
El
4-5
BC-TOP-9-A
Rev. 2
Table 4-3
RESISTANCE YIELD DISPLACEMENT
VALUES FOR SLABS
AND PLATES
21
DESCRIPTION
(1)
YIELD
DISPLACEMENT (33)
RESISTANCE (29,30,31.32)
R-
JR
1R.2
12EI
(,_V2,
b/a
n
121
1.0
1.1
1.2
1.4
1.6
1.8
2.0
.1390
.1518
.1624
.1781
.1884
.1944
.1981
3.0
.2029
0
.2031
POISSON'S RATIO
THICKNESS (in)
2
= MODULUS OF ELASTICITY (Ibrin )
4
MOMENT OF INERTIA PER UNIT WIDTH (in /in)
Mu
21
2"r (Mt + M
Ra2
12EI
-a
I1.,V2 )
b/a
1.0
1.2
1.4
1.6
1.8
2.0
.0671
.0776
.0830
.0854
.0864
.0866
0.0171
4-6
BC-TOP-9-A
Rev. 2
Table 4-4
DUCTILITY RATIOS
(From Reference 28)
Max.
Allowable Value of
Reinforced Concrete
Flexure
10
Beams
0.
1
0p-.
10
Slabs
0.10
ip-p
30
Compression
1.3
A'
As
p'
is
Steel Elements
Members proportioned to preclude lateral
and local buckling
Flexure,
20
Steel Columns
Proportioned to preclude elastic buckling
1.3
e
0.5 eu
y
e u = ultimate strain
ey = yield strain
4-7
a,
BC-TOP-9-A
Rev. 2
1.0
U-
I-z
10"
ULU
C3
10-2
10-2
1.0
10"1
RATIO pn
As
P = bd,
p'=-'-
A's
K!
2
F =-3 __ + pn (1-K)
2n-1
n
.!L
d
nd=-Es
(pn)- - (K --Ai2
0. 10,
K = -m + (m2 + 2q)y
+ 2n-
q = pn (1 + 0.19 P-)
m = pn (1 +
p~
Figure 4-1
COEFFICIENTS FOR MOMENT OF INERTIA
OF CRACKED SECTIONS
4L-8
BC-TOP-9-A
Rev. 2
SECTION 5
SPECIAL PROBLEMS
5.1
s
F(t) -
sin 20t,
(5-1)
(t
where
5.2
5-1
BC-TOP-9-A
Rev. 2
Z (X) = g/a + bA0(1-2ax)/2a 2 + e- 2 ax(V0 2 -g/a-bAo/2a 2)
V0 2L+
2
g(e 2aL
(0 <x<L) (5-1)
-1)/a
(m-1)
(5-2)
> L)
Notations used above are defined at the end of this section. Missile penetration in a liquid can be catagorized by the following cases:
5.2.1
5.2.1.1
If
Z1 (x) is
If
Z (X) is
5.2.2
[1I
(H)]
1/2
(5-3)
5.2.2.1
Z2 (x) is
L (Z 2 (L) < 0)
If
L (Z 2 (L)
If
Z 2 (x)
If
Z2 (x)
is
>
0)
There are two
H (Z 2 (H)
<
0)
It will penetrate a
0, and then float to the
V .
Z2 (H)l 1/2
5-2
(5-4)
BC-TOP-9-A
Rev. 2
In case the missile shape does not have a uniform cross-sectional
area,
refer to equations (D-23) and (D-36) in Appendix D.2 for more general
solutions.
5.2.3
DEFINITIONS OF NOTATIONS
a - yA0 CD/2W
(5-5)
(5-6)
yg/W
g = gravitational acceleration
(g - 32.17 ft/sec2 at sea level)
W = weight of missile
Y - weight density of liquid
terminal velocity
0
5-3
(5-8)
BC-TOP-9-A
Rev. 2
Table 5-1
FLOW (25)
DRAG COEFFICIENT FOR VARIOUSLY SHAPED BODIES IN INCOMPRESSIBLE
Form of Body
L/d
Circular disk
Tandem disks,
L d -
spacing
diameter
Rectangular plate,
>10_3
1.12
>103
1.12
0.93
1.04
1.54
>103
1.16
5
20
1.20
1.50
o1.95
CD
1
2
L - length
d - width
length
diameter
to flow)
d = diameter
21
4
7
0.87
0.99
105
0.63
0.74
20
0.90
1.20
5
Streamlined foil (1 : 3 airplane strut)
1.12
0.91
0.85
L = length
>103
1
2
0.35
0.33
>4 x 104
0.07
>103
1.33
0.34
>5 x 10
L = span
d - chord
2j
21
21
Hemisphere:
Hollow upstream
Hollow downstream
Sphere
Ellipsoid (1 : 2,
105
major axis II to flow)
0.5
>3 x 105
0.20
>2 x 105
0.07
>2 x 105
0.05
0
5-4
BC-TOP-9-A
Rev. 2
MISSILE
Tv=v-0
LIQUID SURFACE
tO0
TTV12-
-a -
-1
MISSILE
c. .
H
V
.x
Fb - W717m
2
F - Wa v /1
Ff
"-
-"V
~'N
Figure 5-1
PENETRATION OF A MISSILE IN A LIQUID
5-5
BC-TOP-9-A
Rev. 2
APPENDIX A
CROSS REFERENCE LISTING TO AEC STANDARD SAR FORMAT
This appendix shows the cross reference between sections of AEC's Standard
SAR format and the sections of this topical report.
BC-TOP-9
3.5.4
A-i
BC-TOP-9-A
Rev. 2
APPENDIX B
GLOSSARY
B.1
PENETRATION
B.2
It
is
PERFORATION
SPALLING OF CONCRETE
Spalling is the peeling off of the back face of the target opposite to the
face of impact.
B.4
DUCTILITY RATIO
B.5
B.6
ELASTIC IMPACT
B.7
PLASTIC IMPACT
B-1
BC-TOP-9-A
Rev. 2
APPENDIX C
REVIEW OF EXISTING FORMULAS
C.l
The most common formulas used in determining the local effects of a missile
on a target, such as penetration, perforation, and spalling for missiles
striking either a concrete or steel target, are given in tables C-1 and C-2.
These tables include equations C-l through C-11. These are the current
state-of-the-art formulas on impact analysis, which consists primarily of
empirical methods based on experiments conducted for specific and limited
applications. Generally, the experiments were conducted for the Government
using missiles, such as bombs and bullets, and having velocities above
1000 ft/sec. Current impact analysis assumes that the missile impinges
the target normal to the surface. The effects of the oblique angle of
striking at various velocities are illustrated in figure C-1.
It can be
seen that assuming normal striking of the target is conservative, since a
small deviation from a normal impact decreases the depth of penetration
considerably.
The Army Corps of Engineers and National Defense Research Committee
equations (table C-l) for penetration, perforation, and spalling have a
term, which depends only on the diameter of the missile. However, this
term provides overly conservative results when a low velocity and large
diameter missile is considered. For example: as V. - 0 the penetration
approaches 0.5D; perforation approaches (1.8)D; and spalling approaches
(2.8)D, which is not realistic.
Experimental data with velocities below 500 ft/sec are just beginning to
develop, with the emphasis on the effect of impact on the target. Some
experiments have been completed with missile velocities in the range of
interest. However, the tests were not necessarily conducted for target
information.( 2 1 ) Therefore, available pertinent data are limited.
The modified Petry formula has had the widest application for determining
the penetration of a mislile into concrete targets and is adopted for use
at the present time. It was developed by the Poncelet theory, provides
estimate of penetration, and has functioned best in the velocity range of
interest. Also, conservatism is built into this approach because of the
following:
A.
The angle of striking the target has a large effect if the angle
is greater than 200.
A normal angle of strike is assumed.
C-1
12
BC-TOP-9-A
Rev. 2
B.
C.
area of impact,
and
Even though tde modified Petry formula was developed in 1910, the material
coefficient for penetration, K., has been revised by experiments and is
reported by Amirikian(1 4 ) and shown in figure 2--l.
The BRL formula for perforation of concrete targets is
used.
It
is
/Mv 2 2/3
T .( 2)
672D
The Stanford Equation (table C-2) has the following defined limits of
applicability:
0.1 < T/D < 0.8,
0.002 < T/L < 0.05,
C-2
(C-12)
BC-TOP-9-A
Rev. 2
T =
0.045
-.
0.047 -
ws
D~mW
(C-13) 1 2
where,
2
wV
2g
C.2
Equation (2-9) assumed the residual kinetic energy of the missile after
perforation (Er) is the difference between the kinetic energy of the missile
before impact (Ek) and the energy required to perforate the steel (EP)
2
Er
r
where
M - mass of the missile
-E
MV
r
2
lb-sec2
c-3
MV 2
2
KV 2
2
(C-14)
I2
BC-TOP-9-A
Rev. 2
Solving for V
This equation(27) neglects the mass of the plug which may be punched out of
the target, which would be very small for a steel target; for a concrete
target, the concrete would fracture and not act: in conjunction with the
missile mass.
VP can be obtained from equations (2-3)
will be the velocity to just perforate,
target, t, is used.
C-4
Table C-i
CONCRETE PENETRATION,
PERFORATION,
Remarks
Formula
Identification
A.
X-
Equation
No.
12K Alog
0
p p10
I +
v2
8
215,000
C-i
Modified Petry
(Refs. 13, 14, 15)
c' D2
X-
(Refs.
(1
0215
B.
Modified Petry
-2)]Depth
of penetration for
ness.
XI4 X when t - 3X
\1.5
W D0"215]C"
X- 282
Research Conmittee
(Refs. 13, 16, 17)
Ammnn & Whitney
(Refs. 18, 19).
-4
1 + e -4
0
tU'
282 N W D
D2
"c
"
12
1i000 j
0 2
1.8
Vs
1000
C-3
T = 2X
X is obtained from
Equation (C-1)
oe~
Table C-i
'.0
Equation
B.
Army Corps of Engineers
Remarks
Formula
Identification
T - 1.35D + 1.24X
obtainted from
C-5
Equation (C-2)
No.
X is obtained from
Equation (C-2)
T - 1.23D + 1.07X
C-6
v 1.33
10
a%'
Ballistic Research
Laboratories
(Modified)
T = 7.8
7"8
w
D1
427
4 &c D
C.
O
"100
~1.8
= 2.28D + 1.13X
X is
obtained from
C-8
Equation (C-2)
0.
For V'
C - 3000 psi
C-7
J 13
1.33
(Ref. 13)
v)
X is obtained from
Equation (C-2)
C-9
*0
Table C-1
CONCRETE PENETRATION,
W
Vs
D =
A
p
X =
S=
t
ft
K
p
PERFORATION,
N =
nI
T
Ts
NOTE:
(in.)
2
0
Table C-2
;00t
0
Identification
Formula
Remarks
No.
0.5 MV2
Ballistic Research Lab
(Refs. 2, 3, 13)
3/2
17,400 K2 D
Stanford Research
Institute
46,500
C-10
2
W
,000 T + 1,500&- T)
(Ref. 20)
C-11
a1*~JU
K -
usually
1,)
S -
ultimate tensile strength of the target minus the tensile stress in the steel (psi)
W -
W5
(See Ref.
20)
BC-TOP-9-A
Rev.
On
VS=
V=
12533
V2=
to
SF77777[2
2533
f525
figal
........
6,263
V4=
71,90
45
140
201
PC
..
'1482 Ric,
1474
Ric.
F777 7777
F7=
Vs:
'724
W$t Ri.
Figure C-i
TYPICAL CRATER PROFILES
C-9
BC-TOP-9-A
Rev.
9APPENDIX
D
DERIVATIONS
D.1
EQUATION
(5-1)
(D-l)
12.5g x
where
-R -
deceleration
(ft/sec
x -
g -
(ft)
)
give the relation
wm .
F =--M
(D-2)
R= 12. 5 W x
m
where
(ib)
Wm - weight of automobile
Equation (D-1) is the motion for an undamped linear oscillator with a unit
iero
Its solution with initial
mass and a spring constant equal to 12.5g.
deformation is
x - C sin (12.5g)I/2
(D-3)
W V2
Em
A
g
+-E
x
-. - F
L
2 max max
where
Vs -
D-1
(D-4)
BC-TOP-9-A
Rev. 2
1/2
1 ~)12V 8a
CM
(D-5)
into equation (D-2)
gives
Vs
sin ( 1 2 . 5 g)1/2 t
i (2
(D-6)
- 0.625 V W sin 20 t
a m
This is a sine wave of frequency w
and then rapidly reducing to zero. Thus under the condition of plastic collision (i.e., missile and target acquire same velocity after impact) the
duration of the impact force is
21
2Tr/w
21
= 0.314 sec.
t = 0.0785 sec.,
At
x
- 3.82 ft
max
gF ave/W
= 25g
/ 1
i
1/2
-)
(78.17) = 3.91 ft
x
2max
12.5
According to the forcing function equation (D-6) the average deceleration
(average over distance, not over time) for Test No. 505-IW is
gF ave/W - gFmax/2W
- (0.
62
5)(78.17)g/
- 24.42g
D-2
BC-TOP-9-A
Rev. 2
D.2
THROUGH A LIQUID
- F(D7
d
(D-7)
where
W - Weight of missile
g - gravitational acceleration
x - depth of missile c.g. below the initial c.g. as shown in figure 5-1
t - time after initial contact of missile with liquid
Fb - buoyant force
Fd
drag force
0 and x
Fb
y rf
0
A(xl) dxl
D-3
Yf(x),
(0 s x s L)
(D-8)
BC-TOP-9-A
Rev. 2
where
a constant
Fb = Wy/ym, (x >L)
(D-9)
where
YM M weight density of missile
The drag force is
Fd
Y AmCDV 2/2g
(D-lO)
where
A - maximum horizontal cross-sectional area of missile
m
Table 5-1 from reference 25 lists some typical values of CD for variously
shaped bodies in incompressible fluid flow. Reference 24 (page 35) presents
some CD values for a family of nose shapes. Other references on fluid
mechanics can also be consulted.
D-4
BC-TOP-9-A
Rev. 2
Substituting equations (D-8), (D-9) and (D-10) into equation (D-7) results
in the following two forms of the equation of motion and solutions:
A.
ForO <x<L
R + ax
g -
+ bf(x) -
0 ,
(D-12)
(0 < x < L)
where
and f(x) is
This is
given in
a nonlinear,
a = yAmC D/2W
(D-13)
b = yg/W
(D-14)
equation (D-8).
second order,
nonhomogeneous,
ordinary differ-
.2
y(x) Then if
2
(D-15)
a v
y'(x)
Equation (D-12)
2*(*)'-
where c is
(D-16)
2 R
becomes
12
(D-17)
2g - 2bf(x)
y'(x) + 2ay(x) -
y(x)
2A M/x
order, nonhomogeneous,
and has the solution
ordinary differen-
(D-18)
ef2adx
D-5
e2ax
(D-19)
BC-TOP-9-A
Rev. 2
Substituting equation (D-19)
y(x) _ v 2 me- 2 ax
=g/a-2be
(D-20)
-2ax
G(x) +c c
(O<x<L)
where
"(x) fe"
2
f(x) dx
S e2
A(x,) dx<I dx
It
(D-21
g/a + 2bG(O)
(D-22)
yW
V02
g/a
(D-23)
+ 2b [G(O)
G(x)]1
(0 < x < L)
f 2 ax
G(L) -
[I
1 x A(x
1 .)dx 1] dx}xL(-4
. V2 2 + gY/Yma + e-2aL
-2
{V0 2
- g/a
(D-25)
+ 2b [G(0) -
G(L)]I
D-6
BC-TOP-9-A
Rev. 2
where V1 is
V 2 . & -Y/Y)
(D-26)
12
Consider the special case of i missile with uniform horizontal
cross-sectional area A0.. Thei i A(x 1 ) - A0 . Equation (D-21) gives
G(x) -
je2ax (tx
A6 dx ) dx - A f xe2aXdx
(D-27)
Ae2ax (Zax-l)/4a2
(0 < x < L)
from which
G(O) - -Ao/4a
(D-28)
and
2
G(L)
%e2aL (2aL-l)/4a
(D-29)
- g/a
(D-30)
For x >L
+ 2 + gy/y
Max
This is
-g-
0 ,
(x >.L)
(D-31)
f(x) - alymb
(x ) L)
,-2&x
D-7
(D-32)
gives
(x > L)
(D-33)
12
BC-TOP-9-A
Rev. 2
The integration constant k can be determined by the condition that
at x - L, v
(V 1 2 _ V22 ) e 2 aL
(D-34)
=- [V2
1V22 + e
v -
2a
[2b (G(O)
, (x > L)
(e
eaL
gives
G(L)) +
11/2
+ g
(D-35)
(_L)
(D-36)
(x
(
/al l1/2
In the special case of a missile with uniform horizontal crosssectional area A0 equations (D-28) and (D-29) are substituted into
equation (D-36) to give
v-
V22+
2
ax [bAO (e aL (1
2aL)
1) /2a 2
'1/2z
2, (x > L)
(D-37)
Y/Ym-l)/ll
+ V0 2 + g (e 2L
At x - H, when the missile strikes the target (See figure 5-1) the
velocity V is
x replaced by H.
D-8
or equation (D-37)
with
BC-TOP-9-A
Rev. 2
APPENDIX E
SAMPLE APPLICATIONS
E.1
CONCRETE (PENETRATION,
A 4-inch x 12-inch wooden plank, weighing 108 pounds, strikes at 300 mph
(440 fps) in a normal head-on collision with a reinforced concrete
(f'c " 3000 psi) wall. The plank has a 48 square inch cross-sectional
area with the equivalent diameter of 7.8 inches.
E.1.1
PENETRATION
Penetration is
X
For 3000 psi concrete K
V2)
+ ____s
215000
12 K
Ap Log1 0 (1
p
and
A
324 psf
48/144
Then
+4402
21'5""6 ) -= 3.77 in.
When the thickness of a wall is less than 3 x 3.77 = 11.3 in., the depth of
penetration is given by equation (2-2):
For example,
4(-L
2x
X-
a)]x
1 + e
-.
E-1
BC-TOP-9-A
Rev. 2
PERFORATION
E.1.2
Vs \1. 3
T =
1 440 .3
WloO/0
108
427
77.8
- 7.01 in.
SPALL
E.1.3
- 2 T
E.2
STEEL TARGETS
Given:
Question:
Solution:
Then
10
]20
2/3
'
2 x 32. 2 (200)2
T =
672(1)
and t
1.25 x 0.5 -
0.625 inches.
E-2
- 0.5 inches
BC-TOP-9-A
Rev. 2
E.3
STRUCTURAL RESPONSE
(12)(13)
12g
13/g
the required
target strain
m s(2001x12)2
EEs =
m2 +
M1
eiW7~f
(2 +1_0
13
_-+
- 32,440 in.-lb
g g
8Ifdy
8421.7)450,000)(1.2)
12
(10 x 12)(6)
.~~29~0
Ld
14,467 lb
and
RR LL3
x
-n
48EI
3
(14,467) (10 x 12) 3
"
0.80 in.
where the value of the moment of inertia, I, for the beam cross-section is
taken from reference 8, and modulus of elasticity E - 30 x 106 and dynamic
increase factor DIF - 1.2 (table 4-1) have been used.
According to figure 3-1 the maximum strain energy for purely elastic
structural response is
E
--I R
e 2 m e
(14,467)
2(4
E-3
BC-TOP-9-A
Rev. 2
Es
E
a
a m
Pr
1 _ (08)1,
32,440
I
+ 0.5 - 3.30
Since, according to table 4-4 the allowable ductility ratio for a steel
beam under lateral loads is 20, this beam can withstand the postulated
In case other
missile impact if no other loads are acting simultaneously.
loads are present as missile impacts and remain in effect throughout the
20 structural response the required ductility ratio should be evaluated by
2 equation (3-24) instead of equation (3-22).
E.4
-,E2(32.17)
(ll)j/
26.6 ft/sec.
V0 d . (26.6) (7)
0.93 x 10-T
R V
. 10772.0
Since L/d - 17/7 - 2.43 the drag coefficient is, according to table 5-1 for
the case of circular cylinder with axis parallel to flow and with R >. 103,
CD
0.854
rd 2 /r
n(7) 2/4 -
38.5 ft 2 .
gives
aCDA.
2W
(62.4)(0.854)(38.5) .. 0.0051 ft
2 (2 x 105)
E-4
BC-TOP-9-A
Rev. 2
gives
b -12
- (62.4)(32.17)
2 x 10 5
W
0.010 ft-
sec-1
2 x 105
fi A- =(38.
m 7m
(38.57(10
2 5) (17)
According to equation (5-8)
- Y/ym)/a
V2 - [g (1
Since H > L,
305.6 lb/ft
/2
[(32.17) (1 - 62.4/305.6)/0.0051]/
= 70.9 ft/sec.
2
Z2(L)2 - V 2+e+ 2 e-2aL(bA0
a 1'-O [e aL
al (1
(1 - 2aL)
+Vo
+a& ( a2
aL Y/ 1m
a
)1
5027 + (0.8408)
[-
12)
1i
4193]
0
E-5
(po0)
(38.5)
2 (0.0050 2
[e
1 7 3 4 (1 - 0.1734)
1502 > 0,
BC-TOP-9-A
Rev. 2
the value of Z2 (H) should be calculated:
+ g (e 2 aLy/Ym-l)/a}
(70.9)2 + e-2(0"005l)(37)(-4193)-
V a1/2
[2g
(h + H)T
[ 2 (32.17)(11 + 37)71
55.6 ft/sec
E-6
BC-TOP-9-A
Rev. 2
APPENDIX F
REFERENCES AND BIBLIOGRAPHY
F.1
REFERENCES
1.
Rotz, J. V., Yeh, G. C. K., Bertwell, W., Tornado and Extreme Wind
Design Criteria for Nuclear Power Plants, Topical Report, BC-TOP-3,
Revision 3, Bechtel Power Corporation, August 1974.
2.
Russell, C. R.,
3.
4.
5.
Goldsmith, W.,
6.
John M. Biggs,
pp. 202-244.
7.
8.
9.
10.
11.
12.
13.
14.
Gwaltney, R. C., Missile Generation and Protection in Li~ht-WaterCooled Power Reactor Plants, ORNL NSIC-22, Oak Ridge National
,Laboratory, Oak Ridge, Tennessee, for the U. S. Atomic Energy
Commiselon, September 1968.
A. Amirikian, Design of Protective Structures, Report NP-3726,
Bureau of Yards and Docks, Department of the Navy, August 1950.
F-1
BC-TOP-9-A
Rev. 2
15.
16.
17.
18.
19.
20.
21.
22.
23.
24.
Maryland,
21
Army Corps of
1960.
Basic Mechanics of Fluids, John Wiley and
25.
26.
27.
28.
F-2
BC-TOP-9-A
Rev. 2
29.
30.
31.
32.
33.
34.
35.
36.
37.
Rao, N. R. N., Lohrmann, M., Tall, L., Effect of Strain Rate On The
Yield Stress Of Structural Steels, Fritz Laboratory Reprint No. 293,
Lehigh University Institute of Research, (also Journal of Materials,
Vol. 1, No. 1, American Society for Testing and Materials,
March 1966).
38.
39.
40.
Albritton, G. E., Response of Deep Reinforced and Unreinforced Concrete Slabs to Static and Dynamic Loading, ASCE National Meeting on
Structural Engineering, September 30 - October 04, 1968.
41.
42.
F-3
S.,
BC-TOP-9-A
Rev. 2
43.
44.
F-4
BC-TOP-9-A
Rev. 2
F.2
BIBLIOGRAPHY
I.
2.
3.
No.
41,
of Steel
1961.
August 1972.
5.
6.
Horger, 0. J., Metals Engineering Design. American Society of Mechanical Engineers Handbook, McGraw-Hill, New York, 1953.
7.
8.
9.
10.
11.
Rinder, R., Saffian, L. W., Wachtell, S., Cohen, E., Dobbs, N.,
Manual for Design of Protective Structures Used in Exposive Processing and Storage Facilities TR3808, Picatinny Arsenal, Dover,
New Jersey/Amman and Whitney, New York, New York, November 1968.
12.
BN-TOP-2,
Bechtel Corp.,
4.
Polyteknish Forening,
Rev 1.
F-5
12