I R - Shock and Vibration Bulletin: (Part 2 of 5
I R - Shock and Vibration Bulletin: (Part 2 of 5
I R - Shock and Vibration Bulletin: (Part 2 of 5
(Part 2 of 5 Parts)
THE
ir _ SHOCK AND VIBRATION
Vo BULLETIN
Part 2
Ground Motion,
Dynamic Analysis
JANUARY 1972
A Publication of
THE SHOCK AND VIBRATION
INFORMATION CENTER
Naval Research Laboratory, Washington, D.C.
1 DDC
Office of
The Director of Defense
NATIONA TIC b Research and Engineering
INFORM4ATIOtJ SERVICE
Spdngfl.Id, Va. 22151
This
67 document has /(-~
been approved for public release and sale; its distribution is unlimited.
Best
Available
Copy
2_ .
SYMPOSIUM MANAGEMENT
Bulletin Production
.. II for
N ............................
" - o
IF
BOultin 42
(Par -o Prs
THE
SHOCK
. AND VIBRATION
BULLETIN
JANUARY 1972
A Publication of
THE SHOCK AND VIBRATION
INFORMATION CENTER
Naval Research Laboratory, Washington, D.C.
Office of
- The Director of Defense
Research and Engineering
---------------------
CONTENTS,
Ground Motion
-
SINE BEAT VIBRATION TESTING RELATED TO EARTHQUAKE.RESPONSE:SPECTRA .. . 1---
E. G. Fischer, Westinghouse Research-Laboratories; ittjsburh, Pennsylvania
SHOCK INPUT FOR EARTHQUAKE STUDIES USING GROUND MOTION FROM UNDERGROUND
NUCLEAR EXPLOS IONS.......... ........ ...................... ......... 21
D. L. Bernreuter, D. M. Norris, Jr., axd F. J. Tokaz, Lawrence Livermore Laboratory,
University of California, Livermore, California
ROCKING OF A RIGID, UNDERWATER BOTTOM-FOUNDED STRUCTURE SUBJECTED TO
SEISMIC SEAFLOOIR EXCITATIN ...................................... 33
J. G. Hamm2r and H. S. Zwlbel, Naval Civil Engineering Laboratory, Port Hueneme,
California
DEVELOPMENT OF A WAVEFORM SYNTHESIS TECHNIQUE-A SUPPLEMENT TO RESPONSE
SPECTRUM AS A DEFINITION OF SHOCK ENVIRONMENT ...................... 45
R. C. Yang and H. R. Saffell, The Ralph M. Parsons Company, Los Angeles, California
iii
Dynamic Analysis
Invited Papers
SMALL SHIPS-HIGH PERFORMANCE
Rear Admiral H. C. Mason, Commander, Naval Ship Engineering Center, Washington, D.C.
v2
<~.-rA~W.CJ
.'.<fl~ . . . rr,~-.'V~rt
. .- ,--. . . . .
Specifications
SURVEY, OF VIBRATION TEST PROCEDURES IN USE BY THE AIR FORCE
W. B. Yarchb, Air Force Flight-Dynamics Laboratory, Wright-Pitterson Air Force
Base, Ohio
I¢
DEVELOPMENT OF'A ROUGH'ROAD SIMULATOR AND SPECIFICATION FOR TESTING
OF EQUIPMENT TRANSPORTED IN WHEELED-VEHICLES
H. M. Forkiois and E. W. Clements, Naval Research Laboratory, Washington, D.C.
LABORATORY CONTROL OF DYNAMIC VEHICLE TESTING
J. W, Grant, U. S. Army Taik-Au"omftive Command, Warren, Michigan
IMPACT VULNERABILITY OF TANK CAR HEADS
J. C. Shang and J. E. Everett, General American Research Division,
General American Transportation Corporation, Niles, Illinois
Test Control
ON THE PERFORMANCE OF TDM AVERAGERS IN RANDOM VIBRATION TESTS
A. J. Curtis, Hughes Aircraft Company, Culver City, California
A MULTIPLE DRIVER ADMITTANCE TECHNIQUE FOR VIBRATION TESTING OF
- COMPLEX STRUCTURES
S. Smith, Lockheed Missiles & Space Company, Palo Alto Research Laboratory,
Palo Alto, California, and A. A. Woods, Jr., Lockheed Missiles & Space Company,
Sunnyvale, California
vi
1'
'1p
vii
. . . ..o.
PAPERS APPEARING IN PART 4
Isolation and Dmping
TRANSIENTRERONSE OF REAL DISEPATIVE STRUCTURES
"1. Plunkett, UnWfisity of Minnesota, Minneapolis, Minnesota
DYNAMIC RESPONSE OF A RING SPRING
R. L. Eshleman, lIT Research Institute, Chicago, Illinois
SHOCK MOUNTING SYSTEM FOR ELECTRONIC CABINETS
W. D. Delany, Admiralty Surface Weapons Establishment, Portsmouth, U.K.
METHODS OF ATTENUATING PYROTECHNIC SHOCK
S. Barrett and W. J. Kacena, Martin Marietta Corporation, Denver, Colo.'ado
ix
----------
THE EFFECT OF "Q' VARIATIONS IN SHOCK SPECTRUM ANALYSIS
* M. B. McGrath, Martin Marietta Corporation, Denver, Colorado, and W. F. Bangs,
. National Aeronautics and Space Administration, Goddard Space -'Uight Center, Maryland
Structural Analysis
NASTRAN OVERVIEW: DEVELOPMENT, DYNAMICS APPLICATION, MAINTENANCE,
ACCEPTANCE
J. P. Raney, Head, NASTRAN Systems Management Office and D. J. Weidman, Aerospace
Engineer, NASA Langley Research Center, Hampton, Virginia
EXPERIENCE WITH NASTRAN AT THE NAVAL SHIP R&D CENTER AND OTHER
NAVY LABORATORIES
P. Matula, Naval Ship Research & Development Center, Bethesda, Maryland
x
STRUCTURAL DYNAMIC ANALYSIS AND TESTING OF A SPACECRAL T DUAL TRACKING
ANTENNA
D. D. Walters, R. F. Heidenreich, A. A. Woods and B. G. Wrenn, Lockheed Missiles
and Space'Company, Sunnyvale, California
j Ship's Problems
DETERMINATION OF FIXED-BASE NATURAL FREQUENCIES OF A COMPOSITE
STRUCTURE OR SUBSTRUCTURES
C. Ni, R. Sccp, and J. P. Layher, Naval Research Laboratory, Washington, D.C.
ix
GROUND MOTION
E. G. FISCHER
WESTINGHOUSE RESEARCH LABORATORIES
PITTSBURGH, PENNSYLVANIA
Vibration test criteria are developed for evaluating the earthquake re-
sistance and reliability of electrical switchgear, including sensitive
control devices. A series of sine beat vibrations applied at experimentally
determined, natural frequencies of the equipment is potentially more
damaging than the original seismic motion. The test table input can be
related to the floor response spectra as calculated for a particular power
plant structure and location in an active earthquake zone.
TFrequency. u
I
tistical properties are not affected by a trans- the system filter. The vibration response wave
lation of the origin of time. Consequently, by (see Fig. 1b) appears to be a sine wave at a
means of the ergodic hypothesisp the required single frequency but with pulsating amplitude;
assembly-averaging of random data can be hence the so-called sine beat vibration. The
replaced by the more simple task of time- probability density of the instantaneous values
averaging over a single record of long duration, of the filtered response tends to be normal, or
Gaussian (synmetrical)(5].
Because resonant vibration buildup is an
important engineering phenomenon, it is On the other hand, the density distri-
essential to evaluate the frequency content of bution function of the amplitude variation of the
seismic vibration excitation. For this purpose envelope of the random sine wave (see Fig. 1b)
the quantity "power spectral density" (PSD) is tends to be skewed. It. can be expressed
used as a measure of frequency content of explicitly by means of the well-known Rayleigh
random-type functions. (Important frequency distribution curve, which is employd in studies
effects in earthquake free-field accelerographs of cumulative vibration fatigue. It is also
appear to be limited to a range from 1 to 25 Hz.) called the two-dimensional error distribution
with reference to the "random walk" problem,
Figure 2 shows the frequency spectra which was first solved by Lord Rayleigh(6].
chart of PSD for a random broadband excitation Strength failures caused by random ceismic
and a narrow-band response, the latter demon- disturbances involve a relatively ifw load
strating the filter action of a simple oscillator, cycles, but they can work th, iiat. rial in the
Hence, the acceleration-time histories shown elasto-plastic range of cumulative fatigue.
in Figs. la, b have now been characterized by
means of a statistical analysis as plotted in A special case of-tho fluctuating sine
Fig. 2. (Earthquake grcind motion response wave is) of course, the true sine beat vibration
spectra usually appear as broadband excitation, which can be expressed in terms of two rotating
whereas building floor response spectra corres- (acceleration) vectors as follows:
pond to narrow-band quasi-resonance buildup at
one or more natural frequencies of the building.)
LIP
Z ZZ;
[ArI Fl,t;-
l: 0'.1
- -- ,.3,--
-. ,,,g-.j•merit[8, hock response
93.
hh::~
table motions which, in turn,eisi
o develop test
specified
spectra in the mounted equip-
3
However, a distinction~must be made FILTER ACTION OF THE BUILDING
between the broadband seismic iexcitation at STRIICTURE
ground level and the narrow-ltand floor motion
at the various building elevations-the latter Figure 4 illustrates how (a).the broad-
resultk.ng from the filtering action of'the mech- -band seismic excitation at the base of the
anical-soil-structure system. Alsol i 'the buil dtigcan be' (b) flteredW'andgniied by
following discussion an essential distinction the building structure, anl(c) amplified by the
must be made between the customary shock equipment response at the coincident building
response spectra used in design-- and the natural frequency of 6 Hz.
original motion-time history-used as the basis
for the equivalenitsine beat vibration testing. The resulting floor motion shcwn in
Fig. 4b conslsts of various harmonic oscil-
Figure 3 shows (a) a typical mathe- lations depending upon the different paths along
matical model of a nuclear power plant building, which the ground disturbance-has been propa-
and (b) thedynamically equivalent seriis of gated. Essentially, this motion under the
simple oscillators. Only the three lowest equipment -has been filtered at 6'Hz and magni-
modes of vibration with natuial frequencies of fied 2. 0 times by the building structure. In
2, 6 and 22Hz need-be considered, 'since'they other words, it now approximatesfa random
are the only ones calctated within the most sine wave as shown in Fig. 1b, and in-the form
imp6tant earthquake hazard range ibr vibration of an, equivalent sine beat vibratlon ciabe
buildup from about 1 to-25 Hz. (Similar freq- applied to the test machine mounting~plate in
uencies found'in the equipment;model are not order to'evaluate tfie seismic capabilities of
necessarily coincident with those in -the the equipment.
building.)
JIn turn, the test machine will duplicate
It Is appreciated that in complex struc- (or exceed) the resulting equipment response
tures there may be important cross-coupling motion shown'in Fig 4c, which occuriswhen a
effects and fluctuating response motions when natural frequency of the equipment coincides
natural frequencies are close together. with a natural frequency of the building
However, this situation is usually covered by structure. At 6 Hz with.5 p'rcenT damping,
specifying an envelope-type of building floor the magnification of the' floor motion in-the
response spectra. (Similarly, when there is equipment amounts to 5. 5 times. (An equiva-
interaction of two or more equipment modes, lent sine beat vibratio,itest using 5 cycles/beat
then a more severe response condition might at 6 Hz will produce this same magnification in
have to be recognized.) equipment having 5 percent damping.)
4
horizontal motion-time history for anOperl-tting
2:0 Basis Earthquake (OBE). Presumably the
_ Fmodel
_
Damping
includes a soil-structure interaction
factor, and the OBE has been established on
1 .i25% the basis of the seismic history of a specific
1.56 L9- power plant site.
1-15 The filtering action of the building, as
previously illustrated by Fig.. 4, was based
W 1.2 upon the data to be found and interpreted from
__ _
H 'random
Fig. 5. In other words, the ground motion is
with a peak value of 0. 13 g. At an
. 0.8- upper floor level where the equipment might
I I be located, the input motion has been magnffied
11
. f
-V2 times to give apeak value of 0.27 g corre-
sponding to the horizontal asymptote for freq-
* 0.iI uencies above 30 Hz. (Compare Fig, 4b with
Fig. 5). Finally, for 5 percent damping in the
equipment at 6 Hz, which is the coincident
00 1 13 building-equipment natural frequency; the
0.30.5 1 3 to 30 maximum equipment response acceleration
Equlpmenthf'erequency, Hz becomes 1.5 g. (Compare Fig. 4c with Fig. 5&)
5- OBE horizontal floor response spectra Figure 5 also gives floor motion
Fig. 5response spectra (actually the peak acceleration
at equipment location in building measured on the equipment at each of its natu-
ral frequencies) for several values of percent
critical damping, which can be determined for
the equipment from the motion buildup during
the continuous sweep frequency test. The
-29-
5
maximum overall magnification of the ground-
motion into the equipment at 5 percent damping 25
is 1.5 g/0.13 g = 11.5 times, which shows that
the selective filtering action of the building Il
structure can be quite. influential. (For the Steady Stale,
L El Centro earthquake at 5 percent damping-the
shock response spectra magnification is only •
6
Figure 7 also shows the Q-factors for
steady-state vibration resonance (Q = 100/2 g )
and for random (white noise) excitation equal to 50
the square root of Q. The quasi-resonance .
buildups produced by the sine beat vibrations S Cri/ 15
fall'in between the latter two extremes. By 40 - 5 Damping /
Damplin -
comparison, seismic.motion response magni- 6OIl
f-ResConace Sle
fication factors fall soniewhat below the random Tests /
excitation values. The most damaging sections 2
of earthquake oscillographs usually correspond N
to a sine beat vibration excitation of about 3 2-_
cycles/beat at various preferred frequencies 1 15
from about 1 to 10 Hz.
7I
addition, since equipment fragility levels are
Of course, when there are many equip- REFERENCES
ment natural frequencies distributed over the
range from 1 to 25 Hz, the stair-step construc- 1. John A. Blume, et al, Design of Multistory
tion may indicate that additional off-resonance Reinforced Concrete Buildings for Earth-
-tests are unwarranted. On the otherhand, when quake Motions, Portland Cement Assoc.,
there are only one or two obvious equipment 1961.
natural frequencies, or only building structure
natural frequenciei to consider, then they can 2. E. G. Fischer, et al, "Mathematical Model
be used to start the stair-step construction in Analysis for the Dynamic Design of
both directions to establish additional off- Machinery", SESA, Experimental Mech-
resonatice test points. anics, October 1967.
The primary purpose of the previous 3. S. 0. Rice, Bell Sys. Tech. J., 23, 282
development is to avoid unwarranted cumulative (1944) and 24, 46 (1945).
fatigue and wear of equipment being tested for
earthquake resistance. (As a practical matter, 4. J. W. Miles and W. T. Thompson:
typical complex electrical switchgear systems "Statistical Concepts in Vibration," Chapter
usually introduce nonlinear effects in terms of 11, Shock & Vibration Handbook, McGraw-
snubber springs, clearances, friction, cross- Hill, New York, 1961.
coupling, etc. Equipment resonances appeai
heavily damped, but persist over a relatively 5. J. L. Bogdanoff, et al, "Response of a
wide frequency range and they can be excited in Simple Structure to a Random Earthquake-
411 three directions of testing.) Type Disturbance," Bull. SSA, 51, 2 April
1961.
CONCLUSIONS
6. G. S. Mustin, "Theory and Practice of
In general, there are insufficient Cushion Design," Shock & Vibration Info.
strong-motion earthquake accelerograms avail- Center, U.S. Dept. of Defense, 1968.
able for power plant structures, let alone
specific types of equipment, to establish "an 7. J. Penzien, "Design and Research
acceptable seismic risk". Also, it is generally Potential of Two Earthquake Simulator
agreed that an isolated peak acceleration Facilities," Richmond Fie!d Station, Univ-
response is not a reliable indication of damage. ersity of California in Berkeley.
In switchgear equipment, where possible loss
of principle function is a more important con- 8. G. Shipway, OA New Technique for Seismic
sideraton, there does not appear to be any Shock Simulation,)) Wyle Labs; Norco, Cal.
obvious correlation with field service reports. and Huntsville, Ala.
8
DISCUSSION
Mr. Gaynes (Gaynes Testing Laboratories): How quencies, and determine the resonance modes by
-did you-monitor the relays and-swltches to-determine moitoiligl the dtaplaceiment of the specimen. We al-
wheter they were functioning or not functioning? so have gone quite a bit further. We have gone into
the shock specters"., approach by usigrandom trn-
Mr. Fischer: Usually we had an electrical hook- sients In determining the ability of the product to
up to an oscillograph element. The definition of real- withstand vibration. We will be presenting a paper in
function is quite a touchy point. In some applications, the near future on that particular comparison, and I
such as computers for the Safeguard system, any kind think you will find that your feelings are correct. Us-
of relay flutter is not allowed at all, so one practi- bg the sine beat seems tobe a more severe test than
cally has to use solid state circuitry. Again, depend- using random transients, but If the shock spectrum
ing on the application, somebody gives you a defini- can be adequately defned, I think that random tran-
tion of what constitutes a malfunction. There are ac- sients would be v.aperkir i;y to go.
celerometers located on the structure etc, but as far
as malfunctioning, it seems to be the electrical oper- Mr. Fischer: Well it could be more authentic. I
ation that is the most significant. Of course ifthe cir- agree, but I started out by saying: "In the interests
cult breaker pops open it is pretty obvious that you of simplicity we pretty much stuck to the sine beat.,,
are Introuble. We Justify the simplicity by-saying it is a conserva-
tive test. Frankly, I will take testing with simultan-
. Heous input anytime. With actual earthquake records,
Mr. Haag,(MTh Systems Corporation): We have as I believe Dr. Piunkett pointed out tids morning,
been performing several tests in the seismic shock you just think you are pushing something In acertain
area including the sine beat, and we pretty much con- direction. It can always escape at right angles. When
firm your approach. We feel it is a severe test. To you use three dimensional testing you never test
find the resonances, we do pretty much the same as nearly as severely as we are suggesting with sine
required in MIL-STD-167, where you sweep the ire- beat testing.
* -,-- -
'-. t t.W ' V ° ' ii , ' 4, H - o *. , - ,1$---
.o='A f. .. $,
_ _ J
A basic problem faced by these manufacturers is the acceleration time histories of the floor
that only rarely are those persons responsible response could be included directly in equip-
for the design of this type of electrical ment specifications, it has been customary in
equipment also knowledgeable in earthquake earthquake eniineering to describe the equip-
engineering technology. Adding to this prob- ment vibration environment by a response spec-
le is the fact that there are no standard trum. The use of a response spectrum-has
requirements for seismic evaluation, so a particular signLfcance:as& methbdof char-
manufacturer receivesa-different set of ieis- acteriziing a complex ti.ansient vibration in a
mic specifications for each new nuclear power way that is useful lor determiniag h.
jenerating station. While-it is expected that simple structure will respond to that vibra-
the severity of the DBE wiii vary for differ- tion. In particular, the earthquike-response
ent locations in the country, there is also spectrum for a building floor ie the maximum
considerable difference in the type of speci- response of a series of 3ingle-degree-of-'
fications being used to assure compliance with freedom oscillators that are excited by the
seismic requirements. The greatest source of floor motion. Each oscillator has a fixed
confusion from nearly all types of seismic percent of critical damping, but a different
specifications is5the use of response spectra natural frequency, so the series of maximum
as the principal means of describing the earth- responses gives a good representation of the
quake vibration environment. Therefore, a frequency content of the floor motion. The
brief review of the development and implica- floor response spectrum does not describe the
tions of earthquake response spectra for acceleration versus frequency characteristic
earthquake engineering seems justified, since of the floor motion directly,,and this is the
experience has shown that this-type of descrip- principal sourie-of misun4derstanding-by those
tion often confuses even specialists in other who are familiar with the-use of "the.Fourier
1
areas of structural dynamics and vibrations. spectrum or Fourier series to describe cotp ,ex
signals.
EARTHQUAKE'RESPONSE SPECTRA
12
The diagram shown in Fig. 3 of a Fig. 4 shows an acceleration response
Sseries of mechanical oscillators resting on a spectrum for the horizontal floor acceleration
building floor ha proved useful for explain- time history shown in Fig. 2. Most of the
ing the response spectirum. The floor ancel- horizontal vibratory eitergy fron earthquakes
eration time history, '(t), would be predicted occurs at frequencies below 10 to 15 Hz, and
from a structural dynamics analysis of the the largest response of the building occurs at
b-ilding response to the DBE. The equation of its natural frequencies. Therefore, the most
motion for any one of the michanical oscilla- severe vibrations will be experienced by floor-
torn is wounted equipment with natural frequencies
close to those of the building. Equipment
2 which is quite rigid (natural frequencies
n 2nzn(t) nn above about 30 Hz) will follow the floor motion
exactly. Therefore, the acceleration ampli-
where tude which is approached asymptotically at high
natural frequencies on the response spec.crum
zn(t) - X(t) - y(c) (2) is identical to the maximum floor acceleration
n+ that would be observed from the floor" accel-
is the relative displacement between the mass eration time history (see Fig. 2).
and the floor, is the undamped natural fre-
-quency of the nth oscillator, and C is a
selected value of critical damping ratio.
The'ire are many techniques available for cal-
culating the displacement response time his-
tory, z(t), and the maximum absolute value 60
r, rvency. "erttz
NOalwaO
13
Also, the response spectrum amplitudes for the spectrum of horizontal acceleration, with
these motions vary more than the earthquake minor peaks at 1.5, 2.2, and 6.5-Hz. Negli-
response spectrum for changes in damping. gible response was found above 10 Hz. The
For example, the peak response spectrum ampli- vertical acceleration power spectrum showed a
tude of 1.1 g at 8 Hz for 5 percent damping predominant peak at 1.0 Hz, with lower ampli-
shoin in Fig. 4 can be satisfied by using a tude peaks from 7 to 11 Hz and from 17 to 21
0.11 g, constant-amplitude (Q - 1/2C - 10), Hz. Thecs amplitude and frequency reuirements
sinusoidal vibration at the equipment base. led to the selection of a servovalve-controlled
However, if the equipment being evaluated has hydraulic actuator to power the vibration
its lowest resonance at 8 Hz with only 2 per- table.
cent damping (Q = 25)9 the sinusoidal excita-
tion would produce a peak acceleration of 2.75 Since the type of equipment to be
g compared with about 1.6 g'Sa for the floor evaluated by this seismic vibration simulator
earthquake motion--an "overtest" of 70 p'rcent. was not expected to respond to large displace-
The use of sine beat excitation is preferable ment components at low frequencies (below I
to constant-amplitude sinusoidal excitation Hz), one of the first design decisions was to
because the variation in the response spectrum limit the total actuator stroke to 6 inches
amplitude with damping is less, but it is still (peak-to-peak). This was done in order to
different from that of an earthquake motion. minimize the entrapped oil in the hydraulic
This requires making conservative estimates actuator, thereby keeping the oil-column
of the expected equipment damping or measuring stiffness high and maximizing the high-fre-
equipment damping at resonances-- a procedure quency response capability. To keep the
which is both expensive and possibly inaccurate commanded displacement within this 6-inch
because of nonlinearities. limit, a high-pass filter was used In series
with the acceleration comand signal to
The alternative of using the maximum attenuate frequency components below 1.5 Hz.
flooracceleration amplitude to test at all Analog acceleration signals were recorded on
frequencies may also severely overtest equip- FM tape by digital-to-analog conversion of the
ment. Furthermore, equipment with several original digital computer data for use as
resonant frequencies below 30 Hz will not simulator command signals. Command acceler-
respond to single frequency excitation in the ations and the resulting filtcred displaciments
same way it would respond to an earthquake used as a basis for the facility design are
motion, and the large number of evaluations shown in Fig. 5. Peak accelerations of 0.37 g
required might produce fatigue failures that in the horizorital and 0.25 g in the vertical
would never occur during an earthquake, direction were recorded during five seconds
of the floor vibration response data for the
For these reasons, it seemed desir- DBE.
able to develop a more realistic simulation of
the earthquake vibration environment. The The design of the seismic vibration
principal goal was to be able to evaluate facility included consideration of several
equipment efficiently while reducing the tradeoffs necessary to achieve a "cost-effec-
possibities of subjecting it to vibration tive" facility. Simultaneous vibration in
which might be considerably less severe or more vertical and horizontal directions, while
severe than necessary. It would be quite dan- desirable for achieving the most realistic
gerous to "undertest" the equipment, while simulation, was not incorporated because of
overtesting might result in an unnecessary its considerably greater cost and complexity,
increase in the equipment cost. The coupling effects of simultaneous vertical
and iorizontal floor motions during an actual
earthquake ure reduced in importance by their
SEISMIC VIBRATION FACILITY DESIGN considerable difference in both amplitude and
frequency content. It is quite unlikely that
The first task in designing the the maximum accelerations will occur simmlta-
Battelle-Columbus seismic vibration facility neously in time for the two directions.
was to deffne the amplitude and frequency
response require-ents needed to reproduce the Physical Nize of the simulator was
typical floor acceleration time history shown based on an estimate of equipment sizes to be
in Fig. 2. Preliminary procesaing of the evaluated in the foreseeable future. T Is
acceleration data consisted of numerical limit was chosen to beobout 3000 pour .sin
integration to obtain velocity and displacement weight with maximum base dimensions itabout
signals. The horizontal displacement reached 80 by 48 inches.
* a peak amplitude of 12 inches, with a predom-
inant 0.3-11z frequency plus some additional To provide a table wit's minimum
lower frequency component (;he displacement weight and maximum stiffness, o equilateral
did not return to the initial position within triangle of steel I-bears with three angle-
5 seconds). Acceleration power spectral den- bisccting beams In the center was designed.
sity curves for the 5-second signals were The single hydraulic actuat-.r is attached below
generated by Fast Fourier Transform techniques the table centroid, and the actuator can be
to provide some idea of frequency content, rotated 90 degrees to provide motion in either
Large peaks were found at 0.3 and 0.9 Ilzin vertical or horizontal directions. Normal
!4
[ Rqproduced 'fn
: modes of vibration, of the table when actuateit O
lvailable COPY-
at the centroid vere calculated by means of a
15
L *
effects of oil-column resonance [4]: accelera- Fig. 8 shows the control system in
tion feedback, pressure feedback (inconjunc- block-diagram-form with linearized-transfer
tion vith a pressure control servovalve), or functions for the major cimponents. Primary
controlled bypass flow, for example. servo-system feedback is derived from a post-
tion transducer (DCDT) mounted on the actuator.,
Adjustable bypass flow using a nee- The positin feedback signal is subtracted from
die valve and small tube across-the actuator the filtered, double-integrated acceleration
is an effective, yet inexpensive, method of command signal to provide the error-signal to
c~ntrolling the resonance. The effects of the servovalve. A &odifiedlacceleration feed-
introducing bypass flow (based on a linear back signal from a strain-gage accelerometer
estimate of system response) are sketched in mounted on the table is also summed with the
the root loci of Fig. 7. For comparable error signal to improve system response in-the
values of forward-loop gain (the parameter Ka 10 to 20-Hz range.
represents servo amplifier gain in amps per
volt), the system with bypass flow is far more
stable--this is shown by the reduced angle of a
vector from the origin to the gain-dependent .. ,'.
root. One disadvantage of bypass flow is the
lower frequency of the resonant peak for desir-
able levels of damping.
Rod 20/
Sec 1200
OO
X×z .0 Fig. 8 - Seismic vibration facility
Jos control system
1100
SEISMIC VIBRATION FACILITY PERFORANCE
at 5 05 0o Closed-loop frequency respons:
=Rod 'oo .1"' evaluations were conducted with a '1500-pound
sec dummy equipment load on the table at normal
acceleration amplitudes. A constant-amplitude
iloo voltage command to the double integration
0 No bypol flow circuit wgas varied from 1 to 40 Hz, and the
resulting table acceleration was monitored
using a Quan-Tech Model 304TD-wave analyzer to
0 track the fun "uental component of acceleration
with a 1-lz bandwidth filter. The results of
the sinusoidal closed-loop frequency response
evaluation are shown in Fig. 9. A second-order
•Rod resonance with a + 3-dB peak ic evident at 10
See Ht, and the effective system bandwidth (-3 dB
Zo 2oo point) is about 14 Hz for the command signal
1amplitude of 0.2 g.
16 *
. iU
ance by comparing the time histories of the
comand and output signals. However, a com-
parison of the response spectra, Fig. 12, for
0the floor-acceleration command signal and the
o* resulting table acceleration indicates accept-
04. able agreement. The amplitude of the tabl e.
01
2 0 comand ,al acceleration can be increiiid above the 0.37-g
maximum accelerution for the DBE so that. the
response spectrum of the table motion exceeds
o-the
0 required spectrum at all frequencies where
any type of equipment resonance might occur
t (i.e., above 4 to 5 Hz). The increase in amp-
3 litude required to compensate for filter char-
01
oacteristics can be reduced by using an in-
0001 creased rate of filter attenuation (sharper
filter) in the control system or by-preprocess-
00 Equgment Weight- 1500:b ing the original floor-acceleration time his-
05. Mtories with digital computer filter algorithms
o
04 prior to obtaining an analog signal. High-
i ~03 ! ' I actuatori
frequency noise can be reduced by using an
S3 actuator with low-friction seals and extending
requency, orti, the control system bandwidth by improved com-
Fig. 9 -Table acceleration response pensation for the oil-column resonance.
to sinusoidal, constant-
amplitude acceleratiorn comoand
signal
(Al . .. .. . . . . ---
. - - --- --
2.0-
3.0-
17
LI
VIBRATION CHARACTERISTICS
TABLE I
Typical Natural Frequencies of
Electrical Equipment Cabinets
Equip. Natural
Item Equipment Wt., Frequencies, Hz
No. Size, in. lb Horiz. Vert.
i -
:•,1 . , 36 24 90 2350 4.5,
5.0 >25
18
ACKNOWLEDGMNT
* - Constant-amplitude horizontal
sinusoidal vibration at 8 1z The authors are grateful to the
natural frequency. assistance of Mr. Charles Rodman and the other
members of the Mechanical Dynamics Division at
Battelle who have contributed to the design
EQUIPMNT PERFORMANCE and operation of the seismic vibration facili-
ty. They also wish to express their apprecia-
It is encouraging to report that most tion for the cooperation and encouragement of
of the electrical equipment that has been sub- Mr. Julius Tangel of the Public Service Elec-
jected to vibration simulating the DBE for tric and Gas Company.
nuclear power generating stations has performed
(satisfactorily without any significant design
changes. No examples of cabinet structural REFERENCES
1'damage to the main structural members or hold-
down bolts have been observed. The few elec- 1. John A. Blume & Associates, "Summary of
trical failures that have occurred were usual- Current Seismic Design Practice for
ly caused by high-voltage arcing or relays Nuclear Reactor Facilitios", United States
which malfunction. Failures have been ob- Atomic Energy Report TID-25021, Sept. 1967.
served in some meters which are often mounted
on a flexible door panel, but a failure in 2. R.D. Kelly and C. Richman, "Principles and
this type of component seldom effects the Techniques of Shock Data Analysis", SVM-5,
primary functional performance of the equip- The Shock and Vibration Information Center,
ment Naval Research Laboratory, Washington, D.C.,
1969.
19
3. V.H. Larson, "The Control of Acceleration
by Electrohydraulic Shaker Systems", 145
Systems Corp., Technical Bulletin 840.00-1.
4. L.H. Geyer, "Controlled Dampin Through
Dynamic Pressure Feedback",, Moog, Inc.
Technical Bulletin 101.
20
21
?4
explosion (UNE) to excite properly located In addition, no clear relation between EQ
structures to EQ-like ground motion. Several magnitude and peak g-level has been established.
major advantages of using the ground motion Ibis is shown in Figs. 1 and 2. Figure I gives
from UNEs are: (1) the ground motion is avail- a correlation-of g-level vs distance from the
able free as a by-product from nuclear-tests, closest point of observed faulting. The EQ
(2),there is no limit on the size or type of magnitudes range from 5.5 to 8.3. It-should be
structure that could be tested, and (3) true noted that the peak g-level from the San
soil-structure interaction would be achieved. Fernando EQ-was around 1-lg. Figure 2 gives
a correlation by Houcner (111 relating EQ m ag-
The ability to test any type of structure nitude to peak g-level and range. Also shown
and the achievement of true soil-structure . in this figure are the Parkfield EQ (M = 5.7)
interaction are important. For example, under- and the San Fernando EQ (M = 6.3). As can be
ground nuclear reactors and storage containers seen, these clearly do not fit Housner's pro-
can be tested. In fact, for marr types or large jected correlation.
structures this technique of subjecting the full-.
scale structure to ground motion from a UNE The response spectrum has been pre-
seems to be the only practical way to investi- ferred for structural engineering studies of
gate the structure's response to an EQ. strong-motion earthquakes, because it com-
bines both the representation of the exciting
CHARACTERIZATION OF EQ GROUND force and the response calculations. It thus
MOTION lumps together under one representation the
major parameters of interest to the structural
The feasibility of simulating EQ ground engineer. The major disadvantage of the re-
motion with UNE- induced ground motion de- sponse spectrum is that only peak response is
pends on the similarity of the two phenomena. determined. The numbei of near-peak response
The mechanism of energy release is much dif- cycles, which are important for studying post-
ferent in a UNE than in an EQ. There are con- yield behavior of structures, is lost.
siderable similarities and also considerable
differences between the induced ground motions The Fourier spectrum of an input function
depending on the criteria of comparison. The shows directly the significant frequency charac-
choice of criteria will depend upon the purpose teristics of the function, and -from it-the time-
for making the comparison. For example, cri- histcry response of the system can be computed.
teria used by seismologists to compare UNE However, in the study by Jenschke et al. [51 it
and EQ ground motion from the viewpoint of was found that results obtained using the
seismic detection of underground testing [4] are Fourier spectral method were not satisfactory:
quite different from criteria used by structural the sine and cosine Fourier transforms showed
engineers whose basic concern is the response high irregularities of the same order as the
of structures. ground acceleration function, and no correlated
characteristics could be found. These high ir-
In order to show that UNE-induced ground regularities make it difficult to use the sine
motion is similar to that from EQs it is neces- and cosine Fourier spectra for comparison
sary to establish a means of properly charac- purposes. Hudson [6] showed that the Fourier
terizing EQs from the point of view of struc- amplitude spectrum is much more regular and
tural response. Many studies have been made is very closely related to the relative velocity
to develop the best means of characterizing EQ spectrum. It has no advantage over the rela-
ground motion [5-101. However, this has tive velocity spectrum and in fact is inferior in
proved difficult to accomplish, one reason being that the peak levels may be lost; time-history
that only a few records of strong-motion EQs aspects are also lost.
exist (here strong motion is arbitrarily taken to
be a peak g-level greater than 0.1 g). These The power spectral density apprcach is
records show a large variation in pulse shape, attractive in that it allows a probabilistic ap-
time of due'ation, peak g-level, and number of proach to be used. However,. Jenschke
cycles. Some of the complex causes of these et al. [5] found that this approach was inade-
variations are discussed in a paper by quate since ground motions produced by EQs
Trifunac 181. and UNEs are essentially nonstationary phe-
nomena, even for dynamic systems having nat-
The most generally used criteria to ural periods considerably shorter than the du-
characterize EQs for structural studies are ration of the ground motion. This is a serious
(1) peak g-level, (2) response spectra, problem when dealing with the ground motion
(3) Fourier spectra, (4) power spectral density, from UNEs because the duration of shaking is
and (5) time history comparisons. None of quite short-abt.ut 5 to 15 sec as compared to
these criteria are completely satisfactory. 30 sec for the El Centro EQ.
Peak g-level has the advantage of being Several investigators have made attempts
extremely simple to compare. It is also phys- to generate artificial accelerograms using
icclly meaningful. Nevertheless, peak g-level various nonstationary processes. But as
alone is not adequate to characterize EQs (7]. pointed out by Trifunac 181 the models generally
22
1I
Pld.
o,,ld Calif. 27J.one 19(6-M.). iI. Hegnloire, Moto 17A 09 . 195(M-. 7). \ 6 6.\
.C 1940(M -7.1).
2. El Centro, Calif. I8 Maky 22. SonJbt., CostaRica 5sOct. 1950(M 7.7).\
10 3. PugetSound, Wash. 13Apr. 1949(M "7.1). 23. Me,.ico City, Mjsxco I I May 1960(M - 7.2). \, •6
• 5. Longkoch, Calif. IMar.
1 1933(M 6.2). IS. SonJose, CasioRica lB8Nov. 1945 (M -7.0)
6. Kern Sounty,Calif. 21 July 1952(M -7.6). 16. Sishop,Calif. lOApr'. 947 (M. 6.4).
7. El Contre, Calif. 30 D¢c. 1939(Me 6.5). 27. Limo, Peru 17 Ocr. 1966 (M *7.5).\ i
8. Logan, Utah 30 As.g. 1962 (M * 5.7). I8. Ltuyo Say, Aloslr 10 July 1958 (M *7.8).
9. Portland, Oregon 6 Nov. 196?(M -*4.8). 19. PrinceWliam Soun, Alaka 27 Mar. 1964 (M c 8.3).
20.sontlag, oi . et, 1945 (\ ,7,)
0.01 0.1 1 00
Distance -miles
Fig. 1 - Plots of g-level vs distance from closest point of observed faulting. From Cloud and
Perez 010]
used are too simple. This results from a the time histories relative to post-yield be-
lack
strong knowledge
of motion of the actualAtcharacter
accelerograms. of
the present havior of structures.
i time these investigations have not progressed Trifunac (81 suggests using the response
Sfr enough to make use of here, envelope spectrum which is a three-dimensional
plot of magnitude of response of one-degree-of-
Direct
ally not madetime-history comparisons
except to obtain qualitativeareinsight
usu- freedom oscillators vs time and frequency.
Only some accelerograms from the 1040 El
into the nature of the ground motion. Some gen- Centro EQ have been studied in this way; hence
eralized time-history characterization of EQ at the present time this approach-while inter-
ground motion is needed to study the post-yield esting-does not offer a meaningful character-
behavior of structures,. No such characteriza- ization of EQs.
tion exi5ts. Cloud and Perez 1101 suggested
comparing total time the acceleration
y1he was above above discussion summarizes the ap-
a certain g-level Figure 3 shows such a cm- proaches generally suggested to characterize
parison for the Parkfield and El Centro EQs. It EQ ground motion for structural response
is not at all evident from this figure that the studies. It is clear that no adequate single
El Centro EQ was by far the more damaging. characterization exists. Based on this discus-
sion we plan to establish a similarity between
numberOther investigators
of cycles of a given suggest counting
g-level. Both the
this EQ and peak
UNE g-levels
ground motion by showing that
() the from UNEs are within the
behvio
suggestion and ofstrctues
the previous Nosuc chracerza-
one (total time Izaionof23s
range estimated for the strongest EQs. (2) the
above a certain g-level) give some useful infor- response spectra from typical UNEs aresimilar
mation, but they do not properly characterize
tion~~~~ exst.
ugese lod to those of strong EQs. and (3) the duration of
ndPeez11]
gd
1.2 * Figure 4- is a plot of peak g-level versus
range for "typical" nuclear explosive yields of
I0, 100, and1000kt(kilotons). For ranges
less than I DOB the valke plotted is the initial
1.0- NI Pokfield (M=5.7) acceleration. In this regime, the motion is
such that the top several hundred feet of earth
a San Fernando (M=6.3) spalls away from.the lower layers, often caus-
Ing a larger peak'g-level when the spall gap
closes. It would be undesirable to locate test
0.8- _ structures within this spall regime.
It is clear from this figure that the peak
g-level of the strongest EQ can be easily du-
*plicated, or exceeded (if desired), by a UNE.
-- 0.6 See Fig. I for estimated EQ peak g-levels.
Hous
Distance to fault - miles
•0.3
N
Vt
(a) Parkfield, California "
Station 2
.2 0.2
COMPARISONS OF UNE AND EQ GROUND
MOTION "u 0.1
24
10 . t (although several have been proposed). Com-
parison herein will be made relative to the El
Centro (1950, NS), Olympia (1949, N8OE), and
Taft (1952, N21E) spectra; these are the spec-
tra most used for design studies.
*6 40
0
10/
0
20
-
2 3
0 I 2
25
- ,. ~. - ~ "L - 777- 4. *7 4- *.,. * - --.
Parkflstd EQ, Cholame Shandon purposes. These events were chosen since
Array,' Station 1o. 5 N85E they cover a wide range in yield and geological
50 A Sconditions. Figure 7 compares Mississippi (a)
3.3 miles from fault and El Centro. Figure 8 compares Missis-
Peak g-level, 0.46 sippi (b), Aardvark, and Olympia. Figure 9
.40 compares Blanca and Taft, and Fig. 10 com-
pares liallbeak and El Centro.
30-" -A Study of Figs. 7-40-shows that the re-
sponse spectra from typical UNEs are similar
: 20 to those from strong EQs. Furthermore, the
variations between the UNEs and EQs are no
greater than the variations befween the EQs
themselves, as shown in Figs. 5 and 6. Carder
and Cloud [17] also noted this similarity be-
tween EQ spectra and UNE spectra.
0 1 2 3 4 Time-History Comparisons
Figure 11 shows accelerograms from
k00 E m three.UNEs, namely Aardvark, Mississippi,
Parkfield EQ, Cholame and Boxcar (1.2 Mt). (No accelerogram was
- Shandon Array, Station available for Halfbeak; the spectrum was com-
[ 80 - No. 2 N65E - puted from a velocity transducer.) The dura-
C 270 ft from fault tion of strong phase motion (peak g-level
60 Peak g-lvel, 0.5 greater than 0.1 g) is approximately 2 sec for
60- - Aardvark and over 7 'sec for Boxcar. This fig-
ure shows in addition that the duration of strong
"* motion from a UNE increases with larger yield
and also with increased distance from ground
40- zero.
20-,
Figure 12 shows the accelerograms from
the El- Centro, San Fernando, and Parkfield
J
EQs. The duration of strong shaking for these
EQs was Z0, 10, and 1.5 sec, respectively.
The number of near-peak g-level cycles of
ground motion ranged from 3 for the Parkfield
0 1 2 4 EQ to over 10 for El Centro.
Period - ec
The accelerograms of the EQs and
UNEs
shown are quite different. Nevertheless, the
Fig. 6 - Response spectrum from Parkfield EQ duration of strong motion and the number of
(June 1966) at two stations cycles of near-peak g-level for the UNEs fall
TABLE 1
Sources of UNJ, Spectra
Yield Recording Rango Peak
Event (kt) Geology (it) g-Level Ref.
213
I'
-Mississippi UNE (Ikin)
.ii
2 ~ '
iv , . ~
El Centro
EQ (48 km)
0 20 40 60 so 100 12
Undamped natural frequency --rad/sec
Fig. 7 Mississippi UNE v El Centro EQ. Relative velocities normalized by dividing original
velocities by 0.1 g
1.5
0.5 , \...u,,
10 A•
I-j..'
0 2. 40 6 80 t20
0
0 20 40 60 80 100 120
Undamped natural frequency - rod/sec
r'ig. 8 - Mississippi and Aardvark UNEs vs Olympia EQ. Relative velocities normalized by dividing
original velocities by 0.1 g
27
I~~~~ _ _ _ _ _ _ _ _ _I_
1.2 __________
0 80-
Aardvark UNE
.60 -.
* . -Taft EQ (N21 E)
00
*J0.2
~90
- ~80 -Halfieck UNE0
k-5 V- -0.4
840- I ' -Cs.6 1 I1 I
> 0 Af 0 1 2 3 4 5 6
1~0 0.4-
0 1 2 3 0.2-24,0 f
7Natural period -sec 0
-0.2
Fig. 10 -lialfbeak UNE vs El Centro EQ
-0. 4
2 3 4 5 6 7 8 9 10
within the range establisHed by the EQs. It is
therefore noncluded that from the standpoint of The- see
a time history comparison, the ground motion
from any of these UNEs could be considered as Fig. I I - Accelerograms from Aardvark,
being representative of a future major EQ. Mississippi, and Boxcar UNEs.
SUMMARY
The oreoingcomarisns how hatnear-peak g-level cycles for UNEs fall within
7be oreoingcomarisns how hatthe range established by strong EQs.
from the structural engineer's viewpoint,
close-in ground motions from UNEs are sim- From these results it is concluded that
liar to those produced by major EQs. More typical UNEsof generate
specifically: lent to that a possibleground
futuremotion
major equiva-
EQ.
(1) T7he peak g-leveis from UNF~s are Therefore, by predetermining the yield and
location of UNEs (relative to the location of a
compar., -? to tho~se estimated for the strongest test structure), a test structure can be excited
HQs. by grouand motion comparable to that of a future
major EQ.
(2) The upectra !rom typical UNEs are
slnr~ir to tliis,2~r EQs (i.e., the envelope However, to study post-yield behavior of
UCVelOpcir ypical UNE spectra would
.. structures it will be necessary to generate UNE
match an envelope from strong EQs). ground motion with longer pulse durations and
also more cycles of near-peak g-level. This
(3) A time history comparison of ground can be accomplished by sequentially firing
motion from typical IJNEs shows that both the IJINEs or by subjecting test structures to ground
duration of strong motion and the number of motions from UNEs fired over a large time
28
0.4 Parkfield EQ, 6-27-66, 2026 PST-
0.2 Cholme Shandon No. 2
- , _V PA, %A %&-- - I
0 2 6 8 10 12 14 16 18 20 22 24 26 2830
El Centro EQ, 5-18-40, 2037 PST
0.2 -. 1 Strong motion
0 4 6 8 10 12 14 16 18 20 22 24 26 28 30
J
*~1.0-
0.8-
0.6-
4 0.4
0.2 -
-0.2 -
-0.4 --
-0.8
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30
Time - sec
span. Both of these possibilities need much One important fact to note is that lower
further study. yield explosives can be used. The only reason
for using very large yield explosives is to in-
crease the duration of shaking. This is of
In order to properly locate structures a some importance since a large number of lcwer
study must be made to correlate the ground yield explosives are detonated as compared to
nuction
and depthspectra fromofUNEs with yield, geology, the very high-yield explosives. Springer a--'
of burial the explosive. Once this Kinnamannj 181 give a list of U. S, UNEs. This
is accomplished it would be possible to subject list indicates that more than a sufficient num-
test structures to increasingly strong motion ber of UNEs are exploded every year to supply
ranging from elastic response to fi.allysevere the required ground motion for any type of test
damage, program.
29
REFERENCES 4th World Conf. Earthquake Eng., Vol. I,
pp. A2-119-A2-132, 1969
I. F. Tokarz and D. L. Bernreuter, "Com-
parison of Calculated and Measured Re- 11. G. W. Housner, "Intensity of Earthquake3
sponse of a High-Rise Building to Ground Ground Shaking Near the Causative Fault,"
Motions Produced by Underground Nuclear Proc. 3rd World Conf. Earthquake Eng.,
Detonations," Lawrence Livermore Labo- Vol. I, pp. 111-94-111-115, 1965
ratory Rept. UCRL-50977, Dec. 23, 1970
12. D. L. Bernreuter, E. C. Jackson, and
2. D. E. Hudson, "Dynamic Tests of Full A. B. Miller, "Control of the Dynamic
Scale Structures," pp. 127-149 in Earth- Environment Produced by Underground
quake Engineering (R. Wiegel, editor). Nuclear Explosives," in Proc. Symp. on
Prentice-Hall, Englewood Cliffs, N.J., Eng. With Nuclear Explosives (Las Vegas,
1970 Nev., 1970), U. S. Atomic Energy Com-
mission Rept. CONF-700101, Vol. 2,
3. J. Penzien, J. G. Bouwkamp, R. W. pp. 979-993, May 1970
Clough, and D. Rea, "Feasibility Study of
Large-Scale Earthquake Simulator Facil- 13. J. R. Murphy and J. A. Lahoud, "Analysis
ity,' Earthquake Engineering Research of Seismic Peak Amplitudes from Under-
Center, University of California, Berkeley, ground Nuciear Explosions," Bull. Selsm.
Rept. EERC-67-1, Sept. 1967 Soc. Am., Vol. 59, pp. 2325-2341, 1969
30
DISCUSSION
Voice In terms of the duration at certain am- soil conditions in the surrounding areas. For in-
plitudes in those plots that you showed, what was the stance, the acceleration levels in a particular earth-
origin of that.data? quake may vary from place to place depending on the
type of soil or rock encountered. Are you aware of
Mr. Bernreuter: That was from a paper by Cloud any studies that have been done?
and a coauthor in the fourth world earthquake confer-
ence. He was looking at how long the level remained Mr Berreuter: Yes, there have been quite a few
above a certain value; for example,how many seconds studies. Harry Seed, University of Southern Calif-
it remained above 1/2 g or 1/4 g for different earth- ornia has made studies tryiny to relate ampliflca-
quakes. tions of base rock motions to soil structure. There
are also several other papers on the subject.
*yjjThen this was not the actual time he
dwelled at any particular level?
Mr. lUernreuter: No. Mr. Haag: Have they been conclusive inany way?
Voice Was it the time after the initiation of the Mr. Bernreuter: I do not really feel that they
event? have been conclusive. One main reason being simply
the lack of recorded motion. Once the San Fernando
Mr. Bernreuter: Yes. records are studied in detail, one might be able to
understand this better. This is the first time wbhave
Mr. Haag (MTS Systems Corporation): I am in- had so many recording stations for a given earth-
tcrcsted in knowing if any one has done any work In qtake over varied geologies. So there is some hope
relating the earthquake accelerations to the type of that this may be clarified.
31
ROCKING OF A RIGID, UNDSRWATER BOTTOM-FOUNDED
STRUCTURE SUBJECTED TO SEISMIC SEAFLOOR EXCITATION
the shear resistance of the soil prevent sliding of the structure, but the
structure tends to rock rigidly in a manner that compresses the soil beneath
it and may cause lifting off.
INTRODUCTION
34
- --Ti,, ., .-2 .e. .. '- >..,, , = . .... . .. . . i
The shear strength seems to increase almost classical way, and solutions usually used are
linearly with depth below the water-sediment those due to Rankine and to Coulomb (Ref.[7J).
interface. The void ratio apparently does not It would be possible, therefore, to make a
decrease with depth. This might indicate that rough estimate of the horizontal restraint on a
the older, lower deposits have become stronger sunken footing, given the properties of the
through age rather than by consolidation. At soil material and the structural forces and
the soil-water interface, however, one can con- geometry. Lacking specific information about
elude that the material supporting the struc- the sediment properties, it is convenient to
ture will be very soft with a low shearing and assume that sliding is prevented, and to com-
bearing. In fact 89% of the samples reported pute the horizontal restraint required during
in Ref.1] had bearing strengths in the range the response. In a specific location the capa-
0.5 to 2.5 psi. bilities of the bottom material could be com-
pared with what is needed to prevent sliding.
When a negatively buoyant structure is
'placed on the ocean floor, there will be some The thick layer of soft material over-
imnediate settlement due to an almost elastic laying the bedrock will tend to alter the
distortion of the bottom material. Following nature of seismic disturbances occurring in
this, there will be a much slower, and greater, the bedrock as they are transferred up to the
settlement duc to consolidation of the material sediment-water interface. Studies have been
under the additional weight of the structure, made of layered soil masses to show the effect
An equilibrium is eventually reached, and this on a seismic disturbance applied at the lowest
condition is assumed at the time the seismic layer and transmitted by shear from horizontal
disturbance occurs. layer to layer until it reaches the uppermost
layer (Ref.[8]). These studies seem to show
If a footing has settled into the sediment that an input earthquake such as that recorded
it will resist being pulled out because of at El Centro in 1940 will emerge with similar
suction formed beneath it. These break-cut low frequency components but with attenuated
forces are commonly experienced in anchor and high frequency components. Of interest is the
salvage work. It is believed that pullout fact that the peak accelerations of the move-
finally results when a failure mechanism ment do not change appreciably since they seem
occurs in the soil material surrounding the to be associated with the lower frequencies.
footings. It is sometimes considered that the The particular structural size ad type that is
average foundation stress resisting pullout of of interest here responds primarily to the
a footing is comparable to the ultimate bearing lower frequencies, and it is therefore assumed
capacity of the soil material. This was that ignoring the moderating effect of the
assumed in Ref.[2], and was substantiated by layers of sediment would not alter the struc-
actual experience reported in that study. tural response significantly if the bedrock in-
put resembled a large land earthquake. The
It is believed that the rate of loading problem, of course, is that one cannot predict
has an effect on the static resistance of the what the disturbance will be, even in the
soil mass. There are several studies (Ref.[3], bedrock.
[41) that have looked at the dynamic properties
of soil under footings of structures on land. What can be done is to take a record of a
Ref.[51 assumes a homogeneous, isotropic, elas- strong typical earthquake and use it for study
tic medium and derives expressions for a soil purposes as though it were applied to the foot-
resistance function involving a linear combi- ings of the structure. The analytical proe-
nation of the deformation and the velocity of dure established for this can be refined as
deformation. The deformation is multiplied by information concerning sediment properties and
a coefficient K, and the velocity by a coef- depths and bedrock motions becomes available.
ficient C. K and C are functions of the shear-
ing modulus of the soil, its mass density, the
size of the footing, Poisson's ratio of the EQUATIONS
soil, and the frequency of the periodic force
that is loading the soil. Ref.[6] uses a simi- Fig. 3 shows the system of forces acting
lar soil resistance function to study dynamic on the structure when it is displaced so that
field tests on a soil-pile system, with the 0, y and xb are all positive and the motion is
additional assumption that C is also propor- such that the three accelerations 0, y and xb
tional to the depth of the hole. are also positive.
Another possible motion of the footing The center of mass is accelerated horizon-
(assumed to rest on or in a horizontal ocean tally by the combination of two motions: a
floor) is horizontal sliding. Sliding motion motion relative to the base as the rigid struc-
would be resisted by friction and by the shear- ture rotates, and the motion of the base itself.
ing resistance of the sediment material. It is The total horizontal acceleration of the mass
believed that the shearing resistance manifests is thus h0 + xb, and the inertia force is
itself in a way that offers passive resistance
to slipping within the soil along some interior
=
surface outward from the footing. This kind of F 0N(hO + xb)
passive resistance has been studied in a
35
'_'
*y
where H is the mass of the structure and 8 is a with H the mass, a the radius, and Y the proper
multip)ier that incorporates the added mass of coefficient so that YMa 2 is the mass moment of
the water, inertia of the spherical structure about its
own diameter.
The structure is considered to move through
the water with a velocity resulting from the RL and RR describe the static and dynamic
same two motions: A horizontal component of resistance of the soil beneath the footings.
the rotation about the base, plus the base They are assumed to have the following form
motion. The total velocity uf the structure
through twe stationary water is then h6 + ;b,
and the drag force acting on the sphere is RL = RS(L) + t LAL)
10 yHF - OMY
36
? a ~ aaA~a' c ~ ~ a
and the submerged weight Eq. (3) is a set of simultaneous nonlinear
first order differential equations. The ini-
tial conditions at time zero are
37
coefficient is assumed to be unity. The mass numerically equal to one-tenth the slope k of
moment of inertia of the structure about a the static resistance. This static slope was
horizontal diameter of the spherical hull was assumed to be 6,000 lb/in, for both types (b)
assumed to be 2/3 Ma2 with no added mass effect and (c); and the coefficient c was therefore
when the sphere rotates. taken as 600 lb-sec/in. 2 .
Three types of support media were con- The above assumptions imply that the struc-
sidered. The first was an almost rigid rack- ture placed on type (a) material will not settle.
like surface that has the static resistance When the seismic loading occurs, the footings
function shown in Fig. 4(a). The second was a will tend to rock or slide, but will not depress
the supporting medium.
In type (c) material an initial settlement
of
, (Ws4)(l/K) - \ """'oo
4 b"
A 6,000 Wi/n..n Z 6 in.
(a) Rigid bottom material was assumed under each footing. In type (c)
material an initial settlemint of 6 inches was
assumed due to elastic distortion, and an addi-
tional settlement of 18 inc'ies was assumed due
r barltin to consolidation.
uliabearinh
CaacitY
{ Type (a) material was assumed to have un-
limited bearing strength. Types (b) and (c)
were assumed to have an ultimate bearing
strength three times the static bearing pres-
(b) Soft eiato-plastic bottom sure. The pullout strength of type (c) mater-
ial then becomes the force necessary to lift
the footin a distance equal to the consolida-
tion settlement against a resistance of 6,000
lb/in, plus the addition velocity dependent
viscous force.
ultimate
bearin rle
i structure-soil-water system is then
capacity assumed to be excited by a horizontal motion
of the base resembling a strong earthquake on
pullout land. This obviously will not be the same
for, motion as occurs in the bedrock below the sup-
porting material, except perhaps for material
type (a); but indications are that the layers
of sediment will have a small effect on the
(c) Soft llterir exhibiting frequency component that affect this structure
the most. Using the base motion
as though it
Fig. 4 - Static resistance functions were applied to the base of the structure is
assumed for seafloor material therefore not considered unconservative; and
of course one is not sure anyway of the nature
of the movement in the bedrock. It is pos-
soft material that deformed elasticially under sible to model the effects of the sediment
the dead load of the structure but did not con- layers on the motion passing through thcm, but
solidate. When loaded further it followed the this could be done if more specific information
pattern shown in Fig. 4(b). The third type were available about a particular structural
deformed initially under the dead load both by site.
distortion and consolidation. It follows the
pattern of Fig. 4(c). An additional property The input used in this example is an arti-
of type (c) material is that it offers a ten- ficial accelerogram, taken from Ref.[10], that
sile type resistance to having the footing lift is believed to contain all the characteristics
out. of a very strong earthquake, stronger than that
for which recorded data exist. The earthquake
The type (a) bottom material does not Is of 120 seconds total duration and is designed
change with rate of loading so its resistance to represent an upper bound for the ground
does not depend upon A. Types (b) and (c) are motions to be expected in the vicinity of the
assumed to be velocity dependent and to have a causative fault during an earthquake of magni-
resistance function of the form R - kA + cLA. tude 8 or greater.
The coefficient c was asaumed to have a value
38
Ref.[i0] followed the following procedure. having the same support materlal, did not lift
An approximation to white noise was passed off at all. It instead moved horizontally with
through a filter to give the process the desired the ground.
frequency content as determined by the power
spectral density. The resulting stationary In Fig. 6(b) a departure from 6 - 0 does
Gaussian process was then given the desired non- not necessarily mean that lift-off has occurred
stationary properties by multiplying by a suit- because the rotation may be due to vertical
ably chosen envelope. The record was then movement of the footings. Lift-off did occur,
processed and corrected to filter the spurious however, for this case on both sides. This can
long period componets. Response spectra were be seen from Figs. 7(a)and 7(b), which are plots
then calculated and compared to those of re- of the force exerted by the left and right foot-
corded motions. Undesirable frequency compo- ings as a function of time. It is interesting
nents were removed by filtering again. Finally to note that the reaction force can momentarily
the accelerograms were scaled to the appropriate exceed the static bearing capacity becLuse of
intensivities of shaking. Fig. 5 shows the the dynamic effect. Comparing the final por-
resulting simulated records, which are used as tions of the curves in Fig. 6(b) and Fig. 7
the input for this structure. shows that the structure is settling into a
10
- "
t
IN''
. .M 25 A 5M 4 L. . 0.j
. . W
,, tx, V
RESULTS
Fig. 6 shows the displacement xr of the canted position even though the foundation
center of mass of the structure relative to ith reactions are fairly balanced.
base and the angle of rotation in degrees, both
as a function of time for the first 80 or 100 In Fig. 7(c) tie response is less because
seconds of the earthquake. Fig. b(a) As the of the resistance to uplift of the footings.
response if tie supporting material is type (a), In this particular case, none of the footings
Fig. 6(b) is for type (b), and Fig. 6(c) is for broke loose to lift off.
type (c).
In computing the response the assumption
In Fig. 6(a) any departure from 0 - 0 means was made that the structure was restrained from
a footing has lifted off. The maximum lift at horizontal sliding. The- results in Figs. 6 and
0 - 7* is about 3-1/2 feet. The same structure 7 do not include any sliding relative to the
responding in air to the same excitation, and support material. The total horizontal
39
2.0
2.0
I II Iq I I I I I I I
0 0 20 to w0 2 W ~ I of02
6(b) tYr-9b)
1.0
2.0
-2.0L
0 201 0 20 40 50 40 0 70
0 typo Wp
40
0 20 3 0 50 6 0 4
o¢
lo
02.
1.
I I V a I I I
2.0
0
,AA
1O 20 O0
' Z
0t
restraining force required to prevent slid'ng (c) would have a very low scatic shearing
was calculated, however, and the values were as strength, possibly not exceeding 1 psi. It
follows: would appear that types (b) and (c) might de-
velop these forces under the footings assumed
if all the footings were acting. The structure
Table I on the type (a) materiel would have to develop
its horizontal restraint principally through
friction between footings and support material.
Type of Bottom Maximum FS if The force required would be nearly equal to the
Material no Sliding submerged weight of the structure, so it ap-
pears doubtful that the required frictional
(a) 130,000 lb force could develop. Some kind of shear fas-
(b) 111,400 teners might be used to hold the footings
(c) 104,600 against shear. These fasteners would of course
affect the tendency toward lifting off, so the
problem would require re-analysis.
These horizontal forces must be resisted Finally, the peak values of absolute ac-
by friction between the footings and the bottom celeration of the center of mass of the sphere
material, and by the dynamic shearing resis- were computed during the response time in eaci
tance of the material. The horizontal force of the three cases. The peak acceleration was
would probably be shared unequally by the four C.250 g for the structure on type (a) material,
footings if none lifted off; and by perhaps two 0.202 g for type (b), and 0.189 g for type (c).
of them if lift-off occurred. Whether or not The Input acceleration of the base had a peak
horizontal1 resistances of the m.ignitude shown value of about 0.44 g that did not coincide
in 7able I could be developed is;,ot certain, with the peak response-; of the structure.
It is certain that materials of type (b) and
41
F I,
REFERENCES
42
NOTATION
Symbol Meaning
t
of center mass
CD drag coefficient
FD drag force
F% inertia force
.q mass of structure
p mass denqlty
43i
I ,
% --. . '
DISCUSSION
Voice: Did you include any kind of damping to Mr. Zudanz (Franklin Institute), The model
account for the soil radiation? looked to me like you had a weightless type of found-
ation at each support. 'You did not have any coupling
Mr.' Hammer: We have run cases where we cou- between different support points. Is that correct?
pled the damping of the soil with cquivalent struc-
tural damping to get a new equivalent factor, but in
this study we just assume 5% critical damping for Mr..Hammer: No. There were only soil springs
each mode. to ground and not soil springs from support to support.
44
DEVELOPMENT OF A WAVEFORM SYNTHESIS TECHNIQUE -
45
--. -;7T
77 77,~~~ 7 X 7,;
46
Solution of Equations. The solution
of Equation I (Reference 1), substituting A --
£ Equations 2 and 3, is
where
Figure 1. Selected Acceleration
Forcing Function, N - 5
2 - )b t ]
21nt cos n(l
c;2 -cos
t
(11M+l)2 1cos 2 rnt - cos 2 7, -+ 16( m+a ) 1 n ]
lJ-211+4(11_1) 2 (Ifm+1)
+ ~~1-21;
f +6 2Ift T
47
='
Amplitude Coefficients. From the Thus, at each iteration step, the
response spectrum to be matched, M variation of Am is a single valued
numbers of Yk maximum responses corre- function of Yk and the process converges.
sponding to frequencies 1 (k - 1,2,...M) To generate large diagonal elements in
can be selected. If the response spec- [Qkm] , the value of b. for a given Nm
trum is to be matched closely at must be selected so as to produce a near
intermediate frequencies, the choice of resonant response; that is, the maximum
the number and relative 'values of the response at a selected system frequency
frequencies is not completely arbitrary should be dominated by one frequency
nor is the selection of component fre- component of the forcing function and
quencies bm. Procedures for establishing none of the other frequency components
these parameters are discussed later. should contribute substantially to it.
Assume that values have been assigned Amplification Ratio. The amplifi-
to Yk, Zk, and b, and that a trial set of cation ratio of each individual frequency
amplitude coefficients Al, A2,,... based component of the waveform can be deter-
on the desired amplification ratio, have mined from its normalized response
been selected. Either Equation 5, 6, or spectrum (Figure 2). Maximax responses
7 may be searched numerically for the are indicated by the solid lines while
maximum response occurring during the the residual responses are shown by the
time period o Z t V (Tm + 1/ 2lTk) and for dotted and dashed lines.
the time tk at which the maximum response
occurred.
A set of matrix equations can be 00- I N. -
written - L; _TII]I
6.[-MINM SECIRA "- " V 1
URk Q
= X
Lkm] I'm k,m =1,2,3 ... 1 (8 4.0
RSIDLAI.SPCIU.--9---.
48
in equal approximately to the selected To demonstrate the method consider
number of half-cycle oscillations the trapezoidal spectrum shown in Figure
N. 3 and the three collinear points which
are to be matched. The problem is to
select the three frequencies, ;1, ;2,
As an example, assume that it is and ;3 such that both the solution will
desired to match a given spectrum at a converge and, In the region of the spec-
frequency of 10 Hz with a single fre- trum between these frequencies, the
quency component with an amplification spectrum of the synthesized wave will
ratio of 7. The selected component lie close to the line which includes the
should have N n 7 half-cycle oscillations three points.
at a frequency of bm a 10/7 - 1.43 Hz. How-
ever, if the forcing function is to
include components at other frequencies,
convergence of the iteration procedure 0
can be assured only so long as the 0
response of any single frequency compo-
nent is not influenced substantially by C-3
those due to the remaining components. 0
Note also, that the amplification ratio 2
of a single component can only be an odd t
integer and that the minimum amplifica- 0
tion ratio is 3 in order to satisfy the 0
requirement that all motions vanish at LOG FREQUENCY
t = Tm .
49
L
TABLE 1 (Reference 1) Dr (2Nmm) td, (M+1), M
OPTIMUM SYSTEM FREQUENCY RATIOS FOR
DIFFERENT VALUES OF N Note than when Dr ; N, the time-
.. histories of the two frequency components
N 3 5 7 9 11 13 occur consecutively. For match point
frequency ratios shown in Table 1, a
Wi correlation with the desigh spec-
c1.35cose
S=1 1.35 1.6 1.3 trum in the region between these
k_1 frequencies cannot be achieved where the
time delay ratio exceeds N.
Thus, in the example, if II = 5, 12 = 7,
3 = 5,
and N1 and assuming , = 1.0 Hz, ;2 Ranges of time delay ratios for
should be 1.6 liz and w3 = 1.6 x 1.35 = 2.16 values of N from 3 thru 13 which will
Hz. minimize deviations from the design
spectrum in the region between match
Phasing of Waveform Comnonents. point frequencies are shown in Table 2.
Phasing of the waveform components can These ranges were calculated for three
be represented by specifying a delay in frequency components whose frequency
the starting time of each waveform com- ratios wer, ai indicated in Table .1,
ponent. Mathematically, the waveform and whose values of N were equal. For
with phased components can be written other frequency ratios and combinations
similar to Equation 3. Thus of N, of course, other ranges of time
delays might yield equally acceptable
N correlations with the design spectrum.
m=1 TABLE 2
fm= 0 0 < t < tdm RANGE OF ACCEPTABLE TIME DELAY RATIOS
= sin 2i, bm(t-td) sin 2, Hmb(t-t) (10) FOR DIFFERENT VALUES OF N
RANGES OF TIME DELAY RATIOS, Dr
tdm i t <tdm + T
m
50
the time delays are forced to stay with- -
51
Figure 4. As noted, the accuracy is
within +15 percent at frequencies of . ..-. -
lower than 0.3 Hz and better than +10 _
percent at higher frequencies. The
required amplification ratio spectrum
and the amolification ratio spectrum of - - _,
the synthesized waveform are compared
in Figure 5. Acceleration, velocity, .i 1-4t i I-
and displacement time-histories are
shown, in Figures 6 thru 8. 47
SUMMARY
A technique has been developed which -d d
simplifies the synthesis of an oscil- .1i il
latory waveform incorporating specified r
amplification ratios and time delays and - T
matching a given response spectrum. The . . . ... .
nature of the basic synthesized waveform
limits the frequency components which
can be included in the waveform and Figure 7. Velocity-Time Function of
requires that their amplification ratios the Synthesized Waveform
be equal to odd integers greater than
unity.
Desnite these restrictions, however,
the technique is a significant improve-
ment over the older cut-and-try methods
and, in addition, has the important
advantage of oermitting the systematic
variation of amplification ratio and
phase, two parameters of vital importance
to the responses of most practical
systems. ,,
FT
.~-
A-
,--- !..,~ ....
_JaJ.
_
~,. .
IrIT
.... ..
I
-
a. F...
-
,
i
T
....
_i _I
j-
52
TABLE 3
SYSTEM FREQUENCIES, COMPONENT FREQUENCIES, TIME DELAYS, AND
CALCULATED AMPLITUDE COEFFICIENTS OF THE SYNTHESIZED WAVEFORM
REFERENCES
1. R. C. Yang, Modification of the
WAVSYN Computer Program, Document
No. SAF-U2, The Ralph M. Parsons
Company, 30 April 1971
I
THE RESPONSE OF AN ISOLATED FLOOR SLAB-RESULTS
OF AN EXPERIMENT IN EVENT DIAL PACK (U)
J. M.Ferritto
Naval Civil Engineering Laboratory
Port Hueneme, California
This paper outlines atest of a horizontal cylinder covered with an earth berm subjected to
the pressure and drag forces in the 300,psi overpressure region from the detonation of the
500-ton high-explosive shot of Event DIAL PACK. The objective of this test was to obtain
information on the response of an isolated floor slab placed on asoil fill inside the concrete
cylinder. Data from seventeen channels of active instrumentation, composed of a pressure
cell, velocity gages and accelerometers, were recorded. Reduction of the data was made.
$ 55
EXPERIMENTAL PROGRAM Instrumentation
Test Structure Seventeen channels of active instrumentation consisting
of eirht structure velocity gages, eight accelerometers, and
The isolated floor slab test was planned as part of Event one pressure cell were used. The layout of the instrumenta-
DIAL PACK [1], a high explosive field test of 500 tons con- t.ca is shown in Figures 6 and 7. The first letter, A or V,
ducted at the Dcfence Research Establishment Suffield indicates an accelerometer or a velocity gage; the second
(DRES), Ralston, Alberta, Canada. The test structure, letter, V or H,indicates vertical or horizontal orientation.
located at an azimuth of 120 degrees 270 feet from ground The data were conditioned, amplified, and recorded on 32-
zero, was at an anticipated side-on pressure range of 300 psi. track tape recorders located in a bunker 2,700 feet from
The structure location and construction detcils are shown in ground zero. An inter-range instrumentation group timing
Figures 1 and 2. It consisted of a 6-foot-inner-diameter 7- system was recorded on one track of each tape recorder and
inch-thick horizontal right-circular cylinder aligned later used in the data reduction. Timing from an NCEL-
perpendicular to the direction of propagation of the blast designed timing generator was also recorded on one track
wave. The invert of the cylinder was 1-1/2 feet below the of all recorders. A detonation zero pulse provided by DRES
natural grade. The cylinder was covered with 2-1/2 feet of was recorded as received directly onto the last track of each
compacted fill forming a tapered earth berm extending along tape recorder to provide a reference for data located on dif-
the sides of the cylinder and one end at a 3:1 slope. A 2- ferent tape recorders.
foot-thick retaining wall held the other end of the earth The gages were mounted on steel plates cast in or bolted
berm in place; a bolt-on steel closure plate was incorporated to the structure. The pressure ceh was installed in a specially
in the retaining wall to provide access into the cylinder. One designed concrete mount which was cast in the wall with the
and one-half feet of compacted soil was placed inside the heat shield flush with the face of the wall. Passive instrumen.
cylinider, and a 6-inch concrete slab was cast over the soil. tation consisted of tree orthogonally oriented reed gages
A spring-mass system, Figure 3, was installed 4 feet from mounted on the slab and a scratch gage (Figure 7) mounted
the rear end of the slab to simulate equipment on the slab. between the retaining wall and the cylinder.
Figure 4 shows the structure under construction, and Figure
5 shows the completed structure.
:=
" 53' - 2'-,-
berm 2'I_\&
6" - cylinderc
wallZ
Figur 1T 1''6
,_ .....
... _ _T_
_
Nt.SeNoe:Stelrenfrcmetomttendaing.
renoreenomtte
l ndaing.1.. Section, B..-B u ,-
Jectionoam
F3 pie cyinde
I..- -- ' o"
270ft
ground zero I C L
270 ft
53 ft x 42 ft (berm dimensions)
A..
1/4.in.-thlck end plates welded to beams
7"ItF 1/4.ln. end plate"
597I
!P7. -
groujnd zeroj
AVi
andAH
V4ad V
VH57 andVH
AV3 and 58
~AV1 and AH2
AH 4
The Pressure gage Was Set in the outside (aec of the RESULTS AND DISCUSSION
retaining wall to measure thc fre.icld sidc-on overpressure.
Horizontal and vertical accclcromcters wcrc set on thle cvI Observed Damagc
indcr ncar each cnd to measure thc horizontal and vertical
components of tlic cylinder motion.. The -lab instrumnenta- Thc detonation of the 500-ton highi-explosivc chiarge
tion was placed in tlac same cross-sectional planes as thc occurred on 23 July 1970. At D+lI hour project personnel
cylinder instrumentation to evaluate the relative motions returned to 6r.L LN322 site. The fireball Wablackened the
of the slab with respect to the cylindei. A vcrtici'l and two retaining wall and the berm with anlayer of carbon dust.
horizontal acclerMCM' etesre set at the center of the floor Twelve inches of crater ejccta were deposited in front of
slab to record floor accelerations in three orthogonal direc the ground zero end of the retaining wall. This level was
tions. The two vertical velocity gages at eath end were reduced to about 2 inches in)front of the closure plate.
installed near opposite edges of the floc. slab to record IThe maximnumn size of the ejrcta was estimated to be 8
rotational motions of the floor slab. The three orthogonal inches. Approximately 18 inches of the berm on the
reed gages were used to provide records of horizontal and groueMd Meo side were compressed and/or blown tway.
vertical shock spectra. The scratchi gage was used to record
relative motion between the retaining wall and the cylinder.
59
A pattern of major diagonal cracks was noted on the Postshot measurements revealed that the slab rotated
wall. S,,vceal of these cracks were observed on the top sur- clockwise approximately 3 degrees. The permanent relative
face of the v'all and extended completely through the wall. horizontal translation of the slab was 1/2 inch away from
The wall had about 1.degree permanent rotation into the the ground zero side of the cylinder. The center of the slab
berm. Some cracks in the concrete were noted around the remained at the same elevation relative to the cylinder. The
closure plate. At D+I day the closure plate was removed permancat horizontal translation of the cylinder with respect
revealing additional cracking in the wall. These cracks to the retaining wall was 3-5/16 inches away from ground
appeared to go completely through the wall. zero.
The concrete cylinder sections suffered significant
dam. 6 v. A horizontal line of compressive failure was located Reed Gages
at about 30 dcgrees from the top of the cylinder on the lec-
ward side of all the sections. A region of major apalling and The plates from the recd gages were removed, and
tensile cracking was noted at about 90 degrees from the com- measurements of the traces were taken and converted to
pression zone (60 degrees from the top of the cylinder) on displacements. These values, plotted in Figure 9, give the
the windward side of all of the sections (Figure 8). The first shock frequency spectra of the slab.
three sections had been unintentiorally oriented so that the The reed gage consists ot'a number of masses on a rigid
splice in the circular reinforcing steel was located in the bs,. The response of equipment mounted on the slab to a
region of the tension failure. Th6 splice appeared to have given shock can be determined fr3m the measured response
failed, and a section of the pipe was pushed inward approx- of a reed having the same frequency and damping as the
imately 4 inches along this line. The floor slab was littered equipment to the shock. Thus, the shock spectra diagram
with spalled concrete and sand. can be used to determine the shock isolation requirements
of the equipment.
1 , -Active Instrumentation
60
11000Shock
Spectra
0~~~
HorzotalRaia
9W
ElHorizontal Transverse
SVertical
011 10 100100
-4
Frequency (cps)
350-
-70'
Time (msec)
61
Table 1. Peak lnstrumer.tation Values
Cylinder
19.6 146.8 2.3 20.1 0.0 8.5
Vertical accelerometer AVI
Vertical accelerometer AV7 112.8 173.9 7.6 22.4 2.4 5.6
Horizontal accelerometer AH2 90.2 52.6 9.9 5.0 2.3 3.2
Horizontal acceleromctar AH8 72.4 86.0 6.3 11.2 0.4 13.9
Floor Slab
Spring-Mass
Cylinder Displacement
-9
Horizontal Displacement (in.)
Figure 11. Displacement of cylinder.
62
-
i -
1 1
II I I I ' I I l l I
-1 1 3 -1 1 3 5
ground zeroundro
down
0-U
>
-5
Figure 12. Displacement of floor slab, windward side. Figure 13. Displacement of floor slib, leeward side.
Figur:- 11, 12, and 13 indicate tie motion of the inches. The amount of rotation shown is about 3 degrees,
structure was initially downward and away from ground confirming the ficld measurement. From Figure 14, the
zero, chcn downward and toward ground zero. Subsequent rotation of the slab began at about 125 msec and rotated
motion at about D+100 msec was upward and away from the full amount at about 200 mscc.
ground zero. The motion toward ground zero is quite Generally the data obtained in the first few hundred
unusual and unexpected, but is believed to be valid. Pos. milliseconds are very reliable. After this period secondary
sible sources of error such as gage rotation, reversed effects may adversely affect the data. The peak values and
calibration, or an erroneous gage were considered and wave shapes in early time (several hundred milliseconds) are
found not to have been present. This effect was obse'ved quite reliable and are relatively the most important. Accel-
on four independent sets of instrumentation. Ovaling of crometers are usually more reliable than velocity gages when
the cylinder was also dismissed as a possible explanation rotation of the gage is suspected. Vcloc.ity gages experiencing
because the direction of ovaling near the gage location rotation as slight as three degrees may erroneously indicate
was inward rather than outward. The initial downward apparent motion of 100 percent of the anticipated value in
movement was caused by the direct compressive wave. A magnitude. Evidence that this has occurred is noted when
possible explanation of the horizontal movement toward velocities at late times (2 seconds) fail to return to zero.
ground zero is the reaction to the magnitude shear wave Accelerometers are relatively insensitive to rotation and are
which enveloped the berm producing planes of opposing in.;uenced only by the prc luct of the sine of the angle of *
shear forces. These forces gave rise to motion as slippage rotation and the component of acceleration iii the perpcn-
occurred. The reverse movement of the structure was dicular direction. Thus, for small rotations the influence
caused by a combination of rebound and direct induced is negligible.
ground shock arriving about 100 mscc after detonation. Integrations of acceleration data give reasonable
The time to peak positive accelerations of the cylinder indications of wave shape and velocity. l)isplacenents
was about 12 msec after the arrival of the blast wave- how- from integration of accelerations should be capable of giving
ever, the time to peak rositive acceleration of the floor slab an order of magnitude of niovemnent and direction. The inte-
was about 110 insec after the arrival of the blast wave. The gration of data represents a second level of confidence.
magnitude of the floor slab accelerations was substantially
reduced compared to that of tile cylinder. Once set in Shock Isolation
motion, the movement of the floor slab was independent
of the cylinder. The slab appears to have remained station- Usiiig all isolated concrete slab "flo.itinmg' on sand within
ary during the initial movements of the cylinder, amid then a horiLontal cylinder reduced the peak a.celcration of 174g's
moved downward until coming to rest on the soil. on the cylinder to about 2 8 g's on the floor slab. The 22-fps
Figure 14 shows the relative vertical displacenient of peak vertical vclocit) of the cylinder was reduced to about
two velocity gages located on tl,e slab aid separated by 36 iI fps. Table 2 compares tihe peak motions and shows the
63
Ai
shock isolation of the isolated floor slab. Table 2. Comparison of Peak Motions
REFERENCES
6-
VV2
~VVll
1.000 2.000
-2
Time (msec)
64
~I,
DISCUSSION
Voice: Concerning the measurement on the floor about 4,800 degrees of freedom, so it is quite a large
slab, does this refer to the isolated mass sitting on problem. We are attempting first to look at it struc-
top of the springs, or ig the slab sitting on top of the turally, and the time step is made as large as possi-
sol? ble keeping the economics in mind in order to be able
to run the problem.
Mr. Ferritto: All the measurements were made
on the actual concrete floor slab. The spring mass Vcice: Do you solve 4,000 dynamic degrees of
system shown in the slides is simply used to provide freedom?
an equivalent equipment on the floor slab.
Mr. Ferritto: Yes, We are running approximate-
Voice: Is it correct that you were not trying to ly 300 time increments to approximate about 100 mil-
establish the effect of isolation using a layer of soil? liseconds, I believe.
Mr. Ferritto: Basically that is correct. The Voice: What is the highest frequency of the mod-
springs of the mass and the weight of the mass were el you are putting together?
selected to provide the frequency that we would ex-
pect from a typical piece of equipment placed on the Mr. Ferritto: It is a very complex problem be-
slab. It actually represented what we tried to approx- cause the loading is quite complicated. We have a
imate in the model, traveling wave. We have a very complex soil-struc-
ture interaction. The extent to which the berm partic-
Mr. Keen (Bell Telephone Laboratories).: The ipates in the problem still remains an unknown. This
high accelerations ii your last slide uould indicate is one rescn we are using a plane strain, finite ele-
the presence of very high frequency data. You also ment analysis with a quadrilateral element rather
mentioned that you planned to use nonlinear finite el- than another type of modeling procedure. The period
ement analysis techniques. Would this predict any of of the cylinder by itself in compression is about 2
the high frequency phenomena which I assume to be milliseconds. The period of the cylinder in flexure is
present? about 11 milliseconds.
Mr. Ferritto: To this date we have been running Mr. Zudans (Fraklin Institute): I am interested
a structural analysis primarily interested in the op- in your 4,000 dynamic degrees of freedom. Are you
timization of the structure itself without looking at modeling it as a two-dimensional infinitely long type
the isolation characteristics. We have made several of strip?
runs. The time step and other information we have
used in sizing the finite element mesh has been scl- Mr. Ferritto: Yes, it is a plane strain model.
lected, basically, to satisfy the structural character- The quadrilateral element that we are using has 12
istics. The structure, as we are now looking at it, has degrees of freedom and approximately 400 elements.
65
A SHOCK-ISOLATION SYSTEM FOR
22 FEET OF VERTICAL GROUND MOTION*
Experimental data from underground nuclear For maximum utilization and mobility, most
detonations are conditioned and recorded by electronic equipment is housed in truck-trailer
sensitive electronic equipment. Normally, the vans or portable buildings on skids. Most
experiment requires the equipment to be located trailers are conventional highway type, but a
relatively close to the explosive source. In few have been fabricated for higher shock
many cases, ground motion induced by these requirements (7 g vertical). Loaded weights
nuclear explosions requires shock-isolation vary from less than 10,000 lb to 72,000 lbs for
systems capable of supporting heavy instrumen- trailers and up to 140,000 lb for skid
tation vans and isolating them from a very buildings. There is a large variety of
severe three-dimensional dynamic environment, electronic equipment and accessories. Shock
fragility levels vary from 1/2 g to more than
During the past several years we have 50 g. When properly mounted, most equipment
designed, developed, and fielded several falls within our medium fragility range - 6 g
different shock-mitigation systems. Two basic vertical.
systems have been standardized and are discussed
in this paper. These standard designs have The ground motion induced by an underground
allowed us to reduce signal attenuation and to nuclear detonation varies considerably,
reduce very large cable costs by locating depending upon yield of device, geology, and
equipment close to the source. Shock-isolation location of interest. The optimum location,
costs themselves have also been reduced, and with respect to reliability and overall costs
overall reliability has been increased by this for diagnostic and other portable instrumenta-
stat.dardization. tion bunkers located on the surface, regardless
of yield, is safely outside the subsidence
To meet requirements for even more severe crater, but not more than half the depth of
ground motion, we have designed a shock- burial from surface ground zero. In this area
mitigation system to withstand vertical ground the surface usually spalls and follows a
motion up to 32 ft/sec, or about 22 ft ballistic path. Figure 1 is a time history
displacement. The new system is a modification curve of surface ground motion for a large
of one of our standard designs, with the usual event, but the characteristics are typical of
crushable materials replaced by a columnar the spall region. The vertical slapdown
energy absorber. This energy absorber allows a acceleration pulse is usually followed by a
longer stroke without an increase in payload horizontal radial pulse, which can be either
initial height. away from or toward surface ground zero.
Horizontal tangential pulses are also occasion-
In this paper,
the instrumented we discuss
testing the and
program, new the
design,
test ally significant.
results. Work performed under the auspices of the
67
The basic requirement of all shock-
mitigation systems is controlled relative
displacement and force transmission between
the shock input and the package or system to
be isolated. The fundamental classification
of shock and vibration systems is the manner
in which it stores, absorbs, or dissipates
energy. A great many different materials and
methods are used in shock and vibration
isolation. A complete discussion of all
.
AV --,., pare-eters involved in shock mounting items for
underground detonations would be quite long;
therefore, in this paper we will describe only
'dII 0 h....,l a few methods that we have been using.
" :,4 In some cases, accurate final position or
.1 V alignment of the equipment with respect to the
ground is important and must be considered
in the overall design. However, in most cases
position and alignment are iot important. This
paper concerns only these cases. This
variance allows us to design around the three-
--------.....------ %,6 dimensional ground-motion environment by
decoupling the horizontal shock components.
Decoupling is achieved by placing the system
on surfaces with very low friction. Vertical
accelerations are isolated by constant force
vs displacement energy absorbers. Ideally,
these absorbers should have negligible rebound
characteristics. With this condition tha
conventional analysis is relatively simple
o-.(see Ref. 1). Accurate measurement of th6 pay-
load weight and center of gravity and the
"-" appropriate sizing of the energy absorbers
minimize the amount of differential vertical
displacements (tilting) of the system. When
required, excessive horizontal displacements
are controlled by nylon tethers. Energy-
Fig. I Typical time history of ground motion absorber deceleration set values are based on
within the spall region. equipment fragility levels and a structural
amplification factor of 1.75. This factor
accounts for the elasticity of the shock mount
DESIGN METHODS structures,
trailer equipmentDynamic
structure. mounting brackets and the
effects of impact
velocity on energy absorber forces are
The overall problem blends together: (1) considered separately.
definition of input, i.e., ground shock
parameters; (2) fragility level or shock
sensitivity of equipment to the input and We have developed several shock mount
1..,ose (i.e., does it record or transmit systems varying in complexity and cost. The
through shock arrival time?); and (3) design selection of a system for a given event depends
of a fail-safe shock mount system to modify primarily on the estimated maximum vertical
the shock environment when required. ground motion as shown in Table 1. The energy-
absorber system does not necessarily have a
The reliability of a system is strongly factor of safety greater than unity with these
dependent on a foreknowledge of the input maximum estimated input conditions, except for
conditions. Prediction of surface ground the incorporation of fail-safe features. The
motion for our purposes relies heavily on respective nominal ground motion values are
empirical data and methods (see Ref. 1). The much less.
shock-isolation design for a given event is
based on nominal and maximum ground motion
estimates in order to obtain a factor of FOAM AND REUSABLE CRIBBING SYSTEM
safety based on energy. The maximum estimates
are based on maximum credible yield and maximum For many events the estimated maximum
scatter of applicable empirical ground-otion ground motion is less than 10 ft/sec and we use
data. Whenever possible, the design is based the foam and reusable cribbing (F&RC) system
only on peak input parameters since detailed for trailers and vans. It includes permanent
time histories of ground motions are very reusable wood cribbing with an angle iron base
difficult to predict.
68
TABLE
I
Shock-Mitigation System Limits
[Maximum
System
Vertical Ground Motion
Velocity, ft/sec Displacement, ft
69
A primary difference in the system is the type
and location of the ene.rgy absorber (see Fig.
5). All remaining components function the same
as the UGC system.
B
The constant-force columnar energy
due to
Momrent absorber is called a TOR-SHOK and is manufac-
hrnl o f tured by A. R. A. Products, Inc. A single
stage of this device consists of two concentric
tubes with a coil of ductile wire forced
between them. The interference fit between
oC the wireand
sliding andtotheforce
tubestheis sufficient to prevent
wires to rotate.
Rotation of the wires is similar to rotating"
the ring of a torus inside out. The resulting
tensile and compressive strains are in the
F plastic range. Contracting or extending the
I-Veitical force
surfaces. Nylon rope tie downs are also used SECT A-A
to prevent excessive horizontal displacement
and to absorb some of the energy. Colun
70
tubes apart axially forces the wires to rotate ratings and cluster arrays of energy absorbers
many times. The repeated cyclic plastic were tested. At some locations the clusters
straining results in almost constant energy were arranged to impart very high eccentric
absorption per cycle of rotation (or inch of loadings into the columns. The FSGC system set
linear stroke) until eventual fatigue failure. up before drop testing is shown in Fig. 6.
( Before incorporating TOR-SHOKs into a Test instrumentation consisted of acceler-i
shock-mounting system we purchased some for ometers and high-speed movies. Thirteen
static and dynamic tests of individual units accelerometers were placed on the top of the
in both the extension and contraction columns, on the FSGC beams and on the equipment
directions. During the static tests loads were inside the trailer. All channels were recorded
uniform and increased breakaway loads occurred on magnetic tape, with seven selected channels
only a few times. Impact acceleration pulses on m scietigraph , f i cev aluatio n
(three times the average deceleration values)
were measured on the payload in about half of
the extension and compression drop tests. The The drop tests are summarized in Table 2.
axial component of the lateral resonant Tilting occurred during the first three tests,
vibration was quite apparent during the but averaged deceleration strokes agreed with
compression drop tests. In one case, this was test parameters. An erroneous center of
a function of two stages and not the overall gravity location, which was discovered before
length of the column. In all tests drop height, testing began, caused the tilting. After the
stroke, and the average deceleration values
agreed very well. There were no changes in load
during the tests, even though one unit
experienced five complete strokes.
DROP TESTS
A series of full-scale drop tests using a
30,000 lb trailer was conducted. Various load
Fig. 6 Prototype FSGC system before testing.
TABLE 2
Drop-Test Results - Average Values
Deceleration
Drop Velocity Desi!n Conditions Actual
Test Height at Impact Force Stroke Stroke
No. (in.) (ft/sec) (g) (in.) (in.)
71
)-
third test the energy absorbers were reset
(compressed) and redistributed in accordance
with the correct center of gravity location.
Figure 7 shows the trailer after the fourth
test during preparation forthe 15-ft drop of
the fifth test. Fifteen energy absorbers were
used during the fifth test, and three of them
failed. They represented 26.6% of the total _
A AA
A ff~T7~7T~N
A, 72
Fig. 7 Prototype F5GC system ready for Fig. 9 Response of equipment inside trailer
15-ft drop test. during tests 2 and 5. ,
impact-tested to failure. This unit finally pads. These are under the protective polyeth-
failed at a swaged tubing flange. ylene covers at the base of each column in
Fig. 11.
Our comprehensive test program identified
discrepancies in the as-built FSGC system.
These were corrected before using the system ACKNOWLEDGMENT
under actual severe ground-motion conditions.
Normally we desire fail-safe shock-mitigation The authors would like to acknowledge
systems. In actual field applications of this the cooperation and assistance provided by
system, we have installed rigid-foam backup Holmes & Narver, Inc., (Nevada Test Site) in
connection with the F&RC and UGC systems and
EG&G, Inc., (Las Vegas) in the drop test
portion of this program. In particular, the
efforts of Mr. E. Fuller and Mr. H. Montalvo
,,Wer Nbof H&N and Mr. P. Hulhall and Mr. R. Nakanishi
of EG&G are gratefully acknowledged.
: REFERENCE
II
I"
73
DISCUSSION
Mr. Fox: Was that fifteen feet? Mr. Peralta (Bell Telephone
sizing the shock absorbers, do youlaboratories):
do any kind of
In
Mr. Bernreuter: No, I think their stroke was calculation as to how much energy they should be
able to absorb?
designed to be twelve feet capacity. In this case they
stroked 12 feet because they bottomed out. In fact, a
Mr. Bernreuter: Yes we do quite a bit of calcu-
couple of the units failed. We picked the trailer
feet off the ground up 15
and dropped it, so it had a velocity lation based on the weight of the trailer, the expected
level, etc. The formulas are worked out,
of about 28 feet per second when It struck the ground, and we have
performed a compute*, analysis.
74
THE C014PARISON OF THE RESPONSE
Neil E.Mechanics
Johnson and Robert D.
Research, Inc.Galletly
75
COPY
2 23
Coordinate System
A finite element model of the bridge struc- to minimize the number of dynamic degreos of
tural system was developed. The natural fre- freedom for the problem. Note that a much more
quencies, mode shapes and related properties detailed model might be desired if an actual
arc identified. Normal mode techniques were highway bridge were to be analyzed using the
utilized in the problem formulation. The techniques described herein. Such a model
STARDYNE Structural Analysis System and its would probably include a more extensive three-
auxiliary subroutines were used exclusively for dimensional model considering both beam and
the modal vibration analysis and the dynamic plate finite elements and a more detailed de-
response computations. scription of soil/structure interaction includ-
BRIDGE MODEL ing pile foundations, if used. The assumed mass
~ asue.
tion ~ ~ ~SpotclmsaBrid
~ Fgwa scaphrct
ge risisan The isazdYN
distribution summarized in Ta.,le 2.
rgadvlpdb
A mathematical model of a typical 5 ix-span, The roadway superstructure and support columrns
singe
sppot
pdesal hghwa brdgewaswere modeled as a series of beam mrembers. Si.c-
developed to illustrate the analysis techniques tion properties evaluated for these beams are -
described herein. The bridge, wnich is shown
i~
Fi.
n Fg.
Fi. Fg.2
9 m~a:d
is
isann iealzaton
dealzaton ff agiven
aassume in Table 3.as These
concrete section
the basic properties
material with
typical structure and is patterned after typical aprpit adfctosfrrifrigsel
~~~~two-lane
utiled ihighway
the bridge fyeo
prolme designs.
raion. The seldtehniued. prpit oilain
Teodes her o enocn
i.Sch a meogeelte.: l
Using the finite element model described
Thew lnidealized bridgeconsiststh of six 120 foot above, the natural frequencies and normal modes -
long toln spn wih heroadway superstruc- of' the bridge were determined through the use of
ture placed approximately 50 feet above the teSAON tutrlAayi ytmcmue
ground urface. Reinforced concrete construe- thegram.eTNE STA l Anas Semomputer
ens ae assumed. toe supporte bridge
sv pa- and in use by a large number of engineering firm
tpe fodane a s, worldwide through CoNtrol Data Corporation Data
typefoudatins.Centers,
Foundation/soEf interaction was idealized as dy is a large
oamic
analysis (up based
system to 6000onDOF) itatic and
the finite ele-
a series of linear support springs. These ment, normal mode bethodof analysis.
springs have been designed to give reasonable The lowest sixty (60) normal modes were
translational and rotational stiffness proper-dermnd Ths oswihcnb gnral
sri te aSsrued nye ob e
Ta soi classed as either horizontal bending, vertical
spxilaring
r srutne s i eT dex bending or longitudinal are summarized in Table
The example bridge was modeled as a lumped . The generalized
s
weight
h
and modal particpa-
ascmtrilwt
n s atredatrtyia asm onrt
tyia srcue
mass, finite element ap shown in Fig. 1. screi
tion factors correspnding to each o these c
This modtl consists of 45 unconstrained nodes mldes are alsoI givenm. in Table A brief de-
each having 6 degrtes of freedom. Only 2 n es scriptono f the normlr ode ethod an the ters
were selct.ed as m,ss poi locations
dt in order in Tble T is given in Appendix A. e
76
Due to the symmetry properties of the ideal- Table4 were used in the various dynamic response
ized structural system, it should be noted that analyses described below. A modal damping fac-
only a portion of the normal modes listed in tor of .05 was used for all response analyses.
17"-0" TA= I
SUMIARY OF SOIL SPRING CONSTANTS
,4 9.1475 x 109 in lb
10"
S6 x 109 in ib
9.)47,5
R3,AM PROPERTI ES
Supp__rt
Colu-n 23114 ,522, , 1,21l.44
1i3 ].2.
.,, ,l'h.
77
TABLE 4
RESULTS OF THE MODAL ANALYSIS
78
PRESENTATION OF THE METHOD OF ANALYSIS Rearranging Eq. (2), the equations of motion
The earthquake response of structures is de- become
scribed in many sources (5, 4]. A variety of I cl'
N 4I 1xl - Wflyl + (5)
computational tools are available to the struc-
tural analyst. Both transient response tech- Just as the terms on the left hand side of the
niques (5] and shock spectra methods (6] have equations can be considered as the forces acting
achieved wide acceptance in predicting the be- upon the system in the absence of ground motion
havior of structural systems subjected to ground (lyl =jW =(Y' -0), the terms of the right hand
motion. This paper addresses the problem of the side can be viewed as those forces Induced upon
response of structures with large characteristic the system by the ground motions lyf when the
length and considers the case where unattenuated mass points are constrained (lxi =.il = x] =O).
strong ground motion traverses the supports of a This provides a straightforward way for construc-
structure at a uniform velocity. The method is ting [c')and [').
developed In a straightforward manner and is out- For seismic analysis the ground motion is
lined below. often described in terms of ground acceleration.
Consider a system with n degrees of freedom Eq. (3) may be rewritten in the following manner
and r resilient supports to ground. Fig. I shows to accommodate these terms.
a typical system. Assume that the equations of Introducing a new setof auxiliary coordinates
motion are to be written in center-of-mass Jz(, the forces induced by ground motion on the
coordinates and the mass matrix of the system constrained system (right hand side of Eq. (3))
Nt] is diagonal. For simplicity, assume the
damping is small and proportional to the stiff- may be written
ness. This is done for ease of the present K]IZI+ [CilCI =']lyl +['Jlfl ()
analysis; other restricted types of damping
could be considered. Modal techniques will be Substituting the transformation
introduced, modal damping will be utilized and
modal coupling due to damping will be neglected. Izi = [T)JYj (5)
In the absence of ground motion the equations
of motion are: where
r
rmIYJI+CI1 -, [K] XI=0 0) niT] is a transformation of ?onstant
wcoefficients;
where n1
n nIz I is an auxiliary set of coordinates
nNq is the diagonal mass matrix of the into Eq. (4), it is seen that
system;
nn [K][Tflyl + [C][Tfl I. ,Jjyl + ,]iJi (61
n[C) is the damping matrix of the system; [
~,~(
'I[C]l -l c'IJ] +[K]ix!-
PkI+'I'h o (P) (T) -- Kf 1 f'C[c'1 (8)
where The two right hand portions of Eq. (0S) are
r identical expressions since it was previously
n[kI] is a rectangular array of stiffness assumed that
coefficients relating ground attach
nodes to adjacent internal nodes of [C) =o{KI and 'c'] =a .'J (9)
the system;
r Furthermore, the existenc- of iK1 is assured
n , a rectangular array of dampini: since !;,] repre-ents the stiffness %atrixof a
coefficients relating ground attach constrained zyster.
nodes to adiacent internal noles of Next, u.'in F1. (',, F'. (1) be writken
h..ay
tree system;
1 j'zx -. C , i*JK ;W 'KJlr7 + ( 1l )
r ly! d'scribes thp absOlute ivround1 motions
(displacements) nt each of the r hoarrn{ in,'rrr, Fh. (10, .eco.es
supports.
79
N4I'i+[C(ikl-fli)+[K(txl-Izi),o (11) Then
Defining y,(t) (Y(t) t ) I
IeI =xl-f, 0 )]lY'(t){ Y(t- 2
I'v~~~~~~~~~lY
KII: I' cIIi ,(t - tr)
Yi'
Noting that
Fig. h shows the time histories 6ira ground dis-
"I =[T] I= Ill
f[K (1i) turbance for the special case of several evenly
spaced colinear ground points.
Eq. (13) may be written
where 04)
where
(j is an arraj of mode shapes o& the un-
di damped constrained system;
IIin a set of Penerallzerl (modil) co,)rdi-
Fig. 3 - Wave Train Propagation nates.
80
-
Substituting and premultiplyine, by [O]T, Eq. (16) that for the May 18, 1940 El Centro earth-
becomes quake. The vertical component of ground
acceleration equal to 0.6 times the hori-
S- V zontal component (Case 2);
0) . The horizontal component of the ground
°
acceleration in a direction 45 to the
On defining, in the usual manner, the generalized longitudinal direction of the bridgc with
parameters a magnitude equal to that for the May 18,
1940 El Centro earthquake. The vertical
[¢]T N([I = r'-aj component of ground acceleration equal to
I JT(CJ( 4],= 0.6 times the horizontal component (Case3).
RESULTS
J I¢J=X Using the loading conditions described above,
,JI I (' ' = a series of nine independent analyses were per-
formed. Each of these analyses resulted in the
The equations of motion become computation of displacements, velocities and
accelerat.,ns at the various nodes on the I leal-
~JI'i + + (21) ized bridge structure. These results can be
further extended to include equivalent nodal
After solving Eq. (21) for I1, the absolute forces and internal member loads and stresses.
accelerations of the system become To briefly illustrate some the results ob-
tained and to provide a basis for comparing the
' ] = li! - IZ - (TJYJ three methods used, the following representative
data are presented.
and 72) Acceleration responses of a point on the road-
way at the midspan of the bridge are presented
I~x!
4,1 J*.fTHY
Iin Fig. 7 through 12. Figs. 7, 9 and 11 show
81
.50
0215
0.5 -
0 5 10 Time, shc
15 20 25 50
s 82
-V
__..1
". 0
Z .25
flt. 7 - Hid Span'Horitontal Acceleratioo History Fig. 8 Kid Span Horizontol ActeleratlonHistory
thiforn Oroua Mtion Grond lotlon
ov4'rg
-|
0-
0 2 3 4 50 1 3 4 3
Tie*, see Time see
Fig. Span Longitudinal Acceleration History
i -d1 Fig. 10 M
kid Sjpn Longitudical Aceleraton History
UniformGround Mtion Moving GrOund W.tio
,oNc Came3
vvvv 3CM
.3 8.2
.22
7v 0
v. -. 2-
-.2
0 1 2 3 4 5 0 1 2 3 4
Tio'. see
Ti- , e
Figure It Kid SpunVertical. Acceleration History ?Ig. I 1 Mid~f
4 lpnertical" lraln.
Wlifors Crownd Flotion ivll ud tin
-.
1
1- UGM
1.0 A- lMM
o a- SS
. ..o..
.. o
It
.- ".-1-'d ",- o\--'."
",I,,
\--,.,-.0*<. , "0 , 10%" % -.-
< --
, .. .0- --
ot m
I I , I I I , i ,I ..
1 3 5 11
7 13 915 17 19 21 P3 25
Span Location
Fig. 13 - lMhimun Horizontal Acceleration Distribution
(Case I)
195
---
O...
-O ...
<)--
O --<'- -O ... - 0- - -0 - - - - 0O-
--O0
,4
~' 1.0 A- MGM
[3- SS
.94
.5
.. ..---. -- GA . --.-
-- A--.-o-A-. - -
01 I I I ,, a I ; I I 2; 2
1 3 5 7 9 !I1 1 7 9 >
Span iocation
Fig. 1, - 1aximt1. Lonitudinal Acceleration Dintributon
(Case ")
1.0 0- UG'
A- IJM
0 V
43/" .1/ / r%, /1\a-
//
l0.5 / .
o - _
1 3 5 9 11 15 1 1) 1 19 21 25 25
Span Location
Fi . 1) - Maximum Vertical Acceleration Distribution
(Case ")
84
5. Richard Rosen, et al, STARDYNE: User's Manual, Note that modal participation factors are com-
Mechanics Research, Inc. and Control Data puted for each of the threetranslational degrees
Corporation, 1971 of freedom.
6. Richard T. Haelsig, DYNRE IV User's Guide,
Mechanics Research, Inc., 1968
where
K 1,( =j
0~ (2)
where
2
o is the 1.th elgenvalue;
th
(r) is the r eigenvector.
The generalized mass of the structure is given
by
(4)
where
fined to be
,mil'ir
Mr
where
'C ~85 -
p.
i1'
Donald F. Hlaskell
Vulnerability Laboratory
U.S. Army Ballistic Research Laboratories
Aberdeen Proving Ground, Md.
The results of the effort are presented in simple graph and nomogram
format for rapid armor areal weight determinations and mine blast-tank
bottom plate evaluations. This method of analysis will facilitate vuln-
erability assessments and engineering design of armor subjected to mine
blast attack.
: _2
,7 ".,-
88
2 55
Deformation Correlation. Table I lists theoret- action. Iowever, until this analysis is
ical results for plate deformation calculated by improved to bring the theory into closer agree-
use of Eq. 1 as well as the deviation of these ment with test, adequate fracturc predictions
predictions from the test results. As indicated, may be made by using the least squares data fit
the average positive error is 14% and the nega- equation as the fracture criterion. The aveiage
tive error is 18.1% with the overall average absolute error between this least squares fit
error -6.4%. The average absolute error is equation and the plate test data is 13.2%.
16.6%. These test and analytical results are Besides, as indicated by Figure 4, the actual
plotted in Figure 3. In this figure the test tank mine damage data points which were not
data has beei. plotted as reduced deformation used in arriving at the least squares equation
amplitude A defined by: and its 2a limits fall well within the 2a prob-
= 2A a 2 b 2 F hI11/2 ability limits of the least squares fit equation.
Ia
2 2
b Y DISCUSSION
versus mine blast energy flux density Em . The General. The deformation relation, Eq. 1, was
test data is as follows: circles represent 5083 derived from first principles without making use
al-iminum and the squares class 2 rolled hon.o- of test data to develop the equation. First
geneous steel armor. Lach data point represents principles were also used in developing the form
a single test plate or the average of either two of the fracture relation, Eq. 2, although the
or three tests. A least squares fit of the data :xact value of the fracture criterion constant K
is shown as the solid (- ) line with the 2o, was obtained from a least squares fit of test
or 95% probability of occurrence, limits as the data.
dotted (----) lines. The dashed line is the
theoretical curve. As indicated, the theoreti- A relatively simple and straightforward
cal curve is slightly below the least squares means for calculating the thickness of a given
curve. Consequently, plate deformation pre- plate that will fracture when attacked by mine
dicted by the present theory is,on the average, blast is presented by the nomograms of Figures
6.4% lower than would be obtained in practice. S through 7. These nlomograms have been prepared
If a conservati'.e method for predicting deforma- from the fracture thickness equatio)n (Eq. 2).
tion is desire(* tile equation corresponding to Figure 5 i. a nomogram for mean fracture thick-
either of the 2o limits could be employed. For ness bosed on Eq. 2 with K r 1.315 where this
example, use of the upper 2o curve to predict value of K has been determined by a least squares
deformation would yield results that 95% of the fit of the plate fracture data. As indicated
time are higher than would be obtained in actual previously, the average error of this relation
practice. is -0.6% and the average abiolute error is 13.2%.
89
11
IV)
90
PLAN VIEW
OF
TEST PLATE
HOLD-DOWN
t PLATE
FRM
TEST PLATE
MINE
300 f
2o
200- 0 0 01
0L6bI9.oil 0 '4
00
100
010 20 30
E m~ los IN-LB/IN!
Figure 3 -Plate Deformation; 0: 5083 Aluminum Test Data, 0:-Rolled Homogeneous Armor
Steel Test Data,-: Least Squares Fit of Data,--: Theory
91
.. 4).* .~.a* *q00O4
4Jq t .4 -1% N'. SNO) O 4 (1
Nq .4C;-1( 4
+ I. + 1 +I
4 +4 +I + II
0 m wo"l '., ~N
m N'ae m
4J 11 4.+1.44 In .11 .
*m
0
~ e N. -N
-q ~Or )
r,
r tfN
N a) N
0
ID t,
-.4
0 0) r,
1"NL) -too l
It in Go N In -TQO
F-.0
03 ~ ~~
CL-. - -
cl-
CC 9
W 00
4l o '4o
CL
41
i-in
CSA4r-NrtCD(
rt,(1 I I 00 CS92a0
LIN
Ic
W h a a Fy ef Deviation
Test Charge Plate Plate Plate Yield Failure From
Case Weight Thickness Width Length Strength Strain Least Sqs
No Material lb in in in KSI in/in
Average * -0.6%
Error
Average
Absolute = 13.2%
Error
30~
2-
/0 -
20 ,.7e / rI~
. o).4 -
Ri
~
0 ~ 0
¢394
Fqgure 4- Plate Froture. 5083 Aluminum Test Dot 0. Po'led homoi$ ecs Armor Steel T¢il
T
00!0: 0, Ton Test Ooto6 LeOst .u' s rI Ot "ti - - t :yt -
93
S) co00 0 0 In 0 In 0 0
. 11 . . 1 1. y .... T., A ~ .L .
C0C
00
2 cr 0n
0 w 00 q 0 0 <
£0 _j3 W NC
W>-
-an
W
00 0 W0smco
U; %on<-.
94
OD -
No,
0
CL
K 0
11 E
0
o
LL
cm-
00
Ir CL
0 fa i
-j Zt
0 0 t w"ifu- :E
40(4M
00
- N F)
0 - ~~ ~ ~ ~CK) IDI
~ n ( j2
'2 40 0
95~
on InJi
-'I
?. 0.20
E
0
22
LL
w0
00
U,, w0
0-*- 0
Z I
962
,I
THE IMPULSE IMPARTED TO TARGETS BY THE
DETONATION OF LAND MINES
Peter S. Westine
Southwest Research Institute
San Antonio, Texas
All of the analytical studies which have 2) Buried charges demonstrate a greater
been conducted assumed that the loading is an efficiency in transmitting impulse to a
air blast phenomena. Personnel at Cornell target than does a blast wave from a
Aeronautical Laboratories [2] used 13RL data spherical charge in air.
on reflected impulses and pressures from ex-
plosive charges in air to estimate loads on In this paper, the author will demonstrate
tank hull bottoms from mine detcnation. These that the impulse imparted to a target by deto-
air blast data were applied without considering nating a land mine is not an air blast phenome-
ground reflection factors or accounting for the non. The basic premise behind these earlier
effects of burial. In another study conducted studies, chiefly that the loading wab caused by
at the Cleveland Army Tank Automotive Plant, a shock wave propagated through air, was in-
AllIson Division of General Motors 131, the correct. Aerojet General's conclusion that
same BRL data for impulsive loads from air buried charges are more efficient was a cor-
blast were not properly used as in the Cornell rect observation; however, they were
97
incorrect in comparing the loads to air blast spheric density, or velocity of sound in air,
ones. The impulse imparted to a target from can be found in Eq. (1), The specific impulse
a land mine explosion is caused by the momen- imparted to the body being loaded is caused
turn in the explosive products from the charge primarily by the momentum of the explosive
and primarily from the soil encasing the products rather than being caused by a shoc.k
charge. This paper considers only the i- wave propagated through air. Although this
pulse, i.e. the area under the positive rres- equation is Zor a spherical explosive charge in
sure history of the load imparted to a target air, a land mine represents a very similar
from a land mine explosion. Peak pressures circumstance, One can think of the ground
and transient loading histories are not in- which surrounds the mine as a charge sur-
eluded, and cannot be discussed until direct rounded by a weak case of soil. The mass
measurements are made of these effects, term, M, is thus the mass of the mine and an
These limitations will not affect structural equivalent mass of soil. The philosophy rep-
studies of hull plate response or vehicle over- resented by Eq. (1) is applicable to a land mine
turning provided the response of the system explosion because the mass of the engulfed *i-
being stiadied falls within the impulsive loading at standoffs typical of mine attacks is very
realm. T41.4 paper will indicate how specific small relative to the mass of the explosive
impulse is di-tr* ,ted over the surface nf products. Ifwe assume that a 20-lb land mine
several targets, particularly wheels and flat is covered by 3 in. of earth and that the hull of
plates. a structure is 17 in. above the center of the
mine, thin we can demonstrate the validity of
LOADING FROM A MINE BLAST this claim. The density of air under ambient
sea-level atmospheric conditions is 0.0765
The nature of the loading applied to P lb/ft3 and the volume of this air enclosed within
structure auch as the hull of a combat vehicle a sphere 17 in. in diaimeter is 11. 9 ft3 . Thus,
from a land mine explosion involves very com- the weight of an equivalent amount of engulfed
plex phenomena. A shock wave propagated air is 0.0765 times 11.9 or 0. 91 lb, which is
through air is only a minor cause of loading very small relative to 20 lb of high explosive.
and not the principal one. Rudimentary calcu- Actually this ratio of air mass relative to mass
lations indicate that the mass of earth and ex- of explosive products is smaller than 0.91 to
plosive products impinging upon the floor plate 20. 0, because the effective weight of soil (ap-
contribute considerably to the loading. There- proximately 100 Ib) is 5 times greater than the
fore, the loading is a very complex wave form weight of the charge in this illustrative example.
that differs considerably from the classical, Later discussions will show howto estimate the
exponentially-decaying wave that is usually effective maos term, M, in Eq. (1).
associated with blasts.
Jack and Armendt [6] have measured tran-
Baker [51 has shown that very close to the sient pressures at the surface of a plate from
charge the normally-reflected specific im- spherical pentolite charges detonated in a
pulse from a spherical charge detonated in air vacuum. This loading is believed to be some-
equals: what analogous to loads caused by land mines
because the pressures are primarily caused by
I/ 2
(2 ME) the explosive products. A typical pressure
4iRS21 () trace as recorded by Jack and Armendt may be
seen in Fig. I. Observe that this normally-
where M = total mass of explosive and reflected pressure history differs markedly
engulfed air from conventional air blast waves. Two shock
fronts may be seen at A and B in the pressure
history in Fig. I. The rest of the wave has a
S = stand-off distance from center gradual rise time and a rounded shape. The
of charge gradual rise time and rounded shape are
caused by the mass of the explosive products
Sreflected specific impulse, impinging upon the plate and the pressure
This equation is applicable whenever the mass transducer in it. Theoretical studies such as
of the engulfed air is less than approximately those by Lutzky [7] indicate that explosions in
one-tenth the mass of the explosive products. an absolute vacuum would have these charac-
If the explosive charge represented by Eq. (1) teristics. Jack and Armendt feel that the ini-
is encased, the mass, M, represents the mass tial, sharp-rising shock front at point A is an
of the case and explosive. Observe that no Initial reflected air blast wave that would dis-
parameters defining ambient atmospheric con- appear entirely were they to have had a corn-
ditions, such as ratio of specific heats, atmo- plete rather than partial vacuum. They offer
98
* Pressure3
Time
Fig. 1. Normally- Reflected Pressure History Taken
In Vacuum Near a Spherical Pentolite Charge
no opinion as to the cause of the second shock so that the spLCific normally reflected pulse,
front at location B ini Fig. 1. The second shock ipR, can be computed, we will consider how
front may be a secondary shock (,sometimes Eq. (1) can be applied to real targets which
called a "pete" wave by other investigators). possess complex configurations. Eq. (1)
Ordinarily, these shocks-are less intense than allows one to compute the normally reflected
the first shock; however, under these circumn- impulse per unit area at various standoff dis-
stances the wave has an atmosphere of gases tances. This equation does not tell one how
from the explosive products through which it to calculate the total impulse, 1, imparted to
can expand; whereas, the Incident wave has
almost no medium through which it can be con-
ducted. A land mine blast will probably give
a loading somewhat similar to that seen in
a wheel or to a plate where the impulses im-
parted to differential areas on these targets
are not normally reflected. Total impulse
can be deterrilned by knowing the projected
I
Fig. 1; however, the relative magnitude of area of a target and the peak normally reflect-
first shock, explosive products, and second ed impulse imparted to this area :)y using
shock wave pressures will vary considerably. Eq. (2).
One can observe in Fig. I that the impulse
(area under the pressure history) is princi- IR =iR Af (2)
pally caused by the explosive products and
that air shocks produce very little impulse. The term I in Eq. (2) is a shape factor which
Probably in a land mine explosion, the initial is a function of target shape and standoff con-
peak would be larger because of the pre- -ditions. This shape factor will be calculated
sence of an atmos.-here, but the secondary for wheels and rectangular plates, both of
shock would be greatly diminished, because of which are common targets.
reflections being transmitted far into the
earth and because the cylindrical shape of a SHAPE FACTOR FOR A WHEEL
land mine would disperse reflections rather
than focus them as in a spherical charge, A To calculate the shape factor for a wheel,
strong possibility exists that a peak pressure consider encounter conditions as presented
will occur which is a shock rather than being in Fig. 2a. The wheel being considered is
*caused by momentum from soil flung upward located directly over the center of the explo-
by a land mine explosion: however, the vast sive charge. We determine the total vertical
majority of the impulse is caused by momen- impulse imparted to the wheel by considering
tum in soil products. the specific imnpulse imparted to a differential
area located on the rim of the wheel. The
LOADING IMPARTED TO COMPLEX total impulse is then computed by integrating
TARGETS the differential areas around that portion of
Before ,iscsslng wat num,,,ecal valuethe rim being loaded. Fig. 2b shows the 4
Befoe dsct~sin
wha nuericl vluespecific impulse being applied to a differential
should be assigned to the mass, M, in Eq. (1) area, dA.
99
r cote '
rfl-cosO)N rsinO~
L._
(a) (b)
Fig. 2. Wheel Traversing a Land Mine
If the wheel is a thin strip of thickness, h, This impulse may be divided into its vertical
then the differential area, dA, equals: and horizontal components. We are only in-
terested in the vertical component as all hori-
dA = hrd0 (3) zontal components cancel because of sym-
metry. The vertical component, dl, of the
The specific impulse, i, as shown in Fig. 2b impulse equals:
at the location of the differential area equals
after substituting for S in Eq. (1): dl = i cos (* + e) cos e dA (8)
or at
We will assume that only that component of the
specific impulse which is perpendicular to theCo(9 cos 8 -r (9)
differential area imparts momentum to the r -, d
target. This assumption would be invalid were
this to be an air shock; however, it is valid Let us define a nondimensional quantity, 6,
because we are considering momentum from suoh that
moving particles of soil. The impulse which
is applied to the differential area equals: 6 1 +A (10)
100
e arc coo (-) () The differentl:n area, dA, equals:
dA:. dy dx (17)
Because of symmetry, the total impulse may
be obtained by multiplying Eq. (8) by 2.0 and Eq. (1) indiciaes that the specific impulse
integrating the resulting expression from 6 directed at the differential area is:
equal to zero to the limit expressed by Eq. (11).
This integration gives: i (2ME)I/ (18)
2
e 41R
I =fzi coso( e) cosO dA (12)
0 Once again we assume that the component of
this impulse which is tangential to the differ-
Substituting Eqs. (3) and (4) into Eq. (12) and ential area does not load the differential area,
recognizing that 2 rh equals the projected and that all of the impulse imparted to this
area, A, of the wheel yields: differential area is caused by the normal com-
0 ponent of the impt,.se, i, expressed in Eq. (18).
J= (2ME) 1 2 A f
_d2 cos 2 $cos(*+0)cos0 dO The impulse imparted to the differential area
4TS 0 (d+ r-r cose)2 equals:
'A S2
= /X 2 7 $ ( 19 )
But Eq. (13) is Eq. (2) provided the shape
factor, 0, for a wheel equals:
0 One calculates the total impulse imparted to
Sird)2 f cos2# cosO cos(*+9) do the entire plate by substituting into Eq. (20).
(d+ r-r cos 0) 2 X Y
If 4f f aAdA (20
Eq. (14) for the shape factor may be simplified 0 0
further if the appropriate substitutions from
Eqs. (5) and (6) are made into Eq. (14). Making Substituting Eq. (16) into Eq. (18), Eq. (18) into
these substitutions and gathering terms yields: Eq. (19), and Eq. (19) and (17) into Eq. (20)
e (6 Cos
2 -1) cos 0 d9 yields:
0= 62 (15) x Y
f6+1-26os3= f f (2ME) _/4Sdxdy (21)
0 0 4r(x2+y2 5 2 )3 / 2
Although this integral is a very difficult one to
compute in closed form, a numerical solution If one multiplies and divides Eq. (21) by S2 XY
is easily obtained on a computer. The results and recognizes that 4 XY equals the projected
from a numerical investigation are presented area, A. then:
in Fig. 3 where shape factor 0 for a wheel is
presented as a function of 6. 2 s
i (2ME) 1/ 2 A $- r Sdxdv
SHAPE FACTOR FOR RECTANGULAR 4n S2 00 (x 2 -+y2 S2)3/2
PLATES
But Eq. (22) is Eq. (2) provided the shape fac-
To calculate the ahape factor for a rect- tor. 0. for a rectangular plate equals:
angulhr plate, consider a plate with a rect-
angular x-y coordinate system having an origin rS3 X dy
at the center of the plate. Because of the sym- E
J
r f0 x
d3
y 2 + $2)3/2
metry which exists when the charge is buried
directly under the center of the plate, we will
consider the vertical impulse imparted to a After performing this integration, one obtains
differential area located in only one quadrant of for the shape factor of a rectangular plate:
the plate. The plate will have half spans of
X and Y. The slant range, R. from the center ! a (Y/S) (XIS) 1 2
2
of the charge to the differential area equals: (X/S)(Y/S) a iY/S)(IS) 4 1] i
R F/(x24y2+S2) (16)
101
Ye,
4-
~7/4- -- '/
1 2 4 6 8 10 20 40 60 80 100
(I + dIR 3061
102
F -'r
TABLE A I fI '6
Parameters For Determining Impulse d 73)
,1
From A Land Mine Explosion
Eq. (26) defines a three-dimensional space.
Fundamental This space can be reducedi in a two-dimen-
Symbol Parameter Dimensions sional space by squaring the second pi term
and multiplying the result by the first pi term
i specific impulse FT/L 2 to form a new dependent parameter. After
W charge weight FL performing this operation, Eq. (26) may be re-
d depth.of burial L written as Eq. (27).
S standoff distance L
p density of soil FT 2 /L 4 iS2 = f"(, (27)
c seismic velocity L/T d3 d3
103
'~' g
normally reflected total impulse could be The wheels were held by a yoke whose other
easily calculated from Eq. (29). end was attached by means of a torsio.a spring
to a rigid wall. This test configuration is
V 2gh(9 effectively a spring-loa~ed pendulum with a
2 Vwheel as the ball of the pendulum. By mea-
surIng the maximum rotation of the yoke and
Table C contains the test number, depth of wheel system as a result of the detonation of
burial, air gap b~tween ground and bottom of various size explosive charges, the impulse
the mass, and the calculated total impulse for imparted to a wheel could be calculated. The
these tests. The depth of burial was mea- explosive charges used in this test varied
sured from the ground to the top of the charge from 0. 106 lb to 0. 468 lb of pentolite. The
rather than to the center of the chrge. *By as- charges were rectangular parallelepipedes
suming that the 5.0 lb charges were cylinders with a constant surface area of 2" by 21" and a
with diameters of twice their thicknesses, thickness that depended upon the size of the
depths of burial and standoff distances fromthe explosive charge. 1/2" of soil was placed
center of the charge were estimated. Because over the charges. Steel wheels which were 7"
the masses being loaded by the explosive in diameter and 1" thick were in direct con-
charge and soil products are very small, the tact with the ground. On a few occasions the
shape factor § is essentially equal to 1.0 in yoke held only one wheel which was centered
these tests. The final two columns in Table C directly over the mine; however, most ex-
present scaled impulse and scaled charge periments were conducted on a 3-wheel array
weight. In subsequent discussion, these two with a 1. 75" center-to-center spacing between
quantities will be plotted to determine the successive wheels. The middle wheel in any
functional format for Eq. (28). The quantities array was always over the center of the land
(IS 2 )/(A~d 2 ) and W/d 3 are dimensional be- mine. In a 3-wheel array, the impulse was
cause the soil constants have been deleted applied to all 3 wheels. To estimate the im-
from this analysis. Throughout the rest of this pulse imparted to the center wheel only, the
paper, the units for I will be (lb-sec), A will following equation was applied and solved for
be (in 2 ), S will be (ft), d will be (ft), and i the impulse imparted
will be (psi-sec).
TABLE C
I S2- W
Versus - Using Rigid Mass Test Data
A§d3 d
104
data indicate that Eq. (30) proportions the im- A straight line fits the experimental data
pulso appropriately. In developing Eq. (30), we presented in Fig. 4 very accurately. The
assune that the impulse imparted to the out- equation of this line as obtained by a least-
aide wheels is applied parallel to.a line from squares fit to the data in Fig. 4 gives Eq. (31).
the center of the charge to the bottom of the
wheel, and that the upward momentum of out- I W0 . 72 d0 . 84
#A = 1. 725 s2 (31)
side wheels is caused by the vertical compo-
nent of this impulse. Table D contains the
total impulse measured by this test arrange- Eq. (31) is the explicit expression for either
ment, the number of wheels in the test array, Eq. (27) or Eq. (1). If one substitutes
the equivalent impulse imparted to the center 1.4 x 10+6 W for E in Eq. (1) and equates the
wheel only, the charge weight, depth of burial; left hand side of Eqs. (1) and (31) after making
and standoff distance. Before computing the the expressions dimensionally consistent, he
total impulse imparted to a center wheel, one obtains the effective mass term, M. M when
must look up the shape factor for the wheel in expressed in slugs is given by Eq. (32).
Fig. 3. Both 6 and # are listed in Table D
for each test configuration. The dependent M(slugs) = 3.49 W0 44
d1 . 6 5 (32)
quantity, (IS 2 )/(At d3 ), and independent quan-
tity, W/d 3 , are computed from this informa- Usually the reader prrers to think of mass in
tion and listed in the last two columns of units of pounds. Mass in pounds is given by
Table D. Eq. (33).
TABLE D
I S2 W
A Versus - Using Wheel Data
IS 2 W
Total Impulse on AI d3 d3
Test Impulse No. Center d S 6 f (psi-sec
No. (lb-sec) Wheels Wheel (ft) (fit) \ft/ ft3
2 23.5 3 12.1 .0832 .3749 1.285 4.3 98.3 400.
3 32.6 3 16.1 .0930 .3847 1.319 4.0 105.5 353.
4 18.6 3 10.7 .0607 .3524 1.208 5.45 155. 475.
5 27.7 3 14.6 .0778 .3695 1.266 4.55 131.5 425.
6 28.5 3 14.1 .0930 .3847 1.319 4.0 92.6 353.
8 53.3 3 24.2 .117 .4084 1.402 '.25 110.5 261.
9 28.1 3 13.9 .0930 .3847 1.319 4.0 91,2 353.
10 43.2 3 20.1 .110 .4017 1.378 3.4 102. 286.
11 47.1 3 21.4 .117 .4087 1.402 3.25 97.5 261.
13 71.4 3 31.7 .125 .4167 1.430 3.1 130. 240.
14 13.7 1 13.7 .0930 .3847 1.275 4.0 90. 353.
15 15.9 1 15.9 .103 .3947 1.354 3.6 90. 311.
16 23.9 1 23.9 .113 .4047 1.387 3.4 113.5 274.
105
600 ' ' ' ' '' 1 ' '
100 X _
60-
is' 40 r
10 m7
f \ . . , WO-7d o.8+
4-"4 " --"s" -+
1.725
44
2
I Ih II I I ii,1l I1
only two types of experiments are used in the ground. On the other hand, if the charge
developing Fig. 4, a significant number of is located at the surface with no cover over
data points (25 points) from a variety of test it, the impulse begins to be an air blast phe-
conditions is used to develop this curve. The nornenon. Three data points were not includ-
charge weights in the wheel experiments were edt in Table C, because in those rigid mass
varied by a factor of 4. 0. Depths of burial on tests the top of the charge was flush with the
the rigid mass tests ranged from 4" to IZ" of stil'face of the ground, and no soil covered the
soil cover, and the air gap varied from zero charge. Because the center of the charge was
to 16". By combining these facts with the ob- below the surface of the ground, W/d 3 equaled
servation that theae data involve two different a finite value of 2570. Extending the least-
types of targets, one can see that 4, W, d. and squares fit in Fig. 4 would predict that scaled
S have all been varied in Eq. (31). inipulse, (I S)/(d3 ), should equal 490. Using
Eq. (t) as applied to air blast with a ground
SUMMARY reflection factor of 2. 0 predicts that scaled
impulse, (I SZ)/(d 3), would equal 4745. The
Fig. 4 and Eq. (31) do have limits of W/d 3 average of the three tests on the rigid masses
for which they apply. Obviously, if the depth was a scaled impulse, (I S2 )/(d 3 ), of 3227.
of burial is too great, the detonation of a Obviously this test was in a transition range
buried explosive will not disturb the surface of where the impulse was changing from being
106
caused by momentum in soil particles to being The author is indebted to Mr. Bruce
caused by momentum in the explosive products Morris, the Army's technical monitor on this
from the charge. Some finite soil cover is re- project and Mr. Alexander Wenzel, SwRI
quired for Eq. (31) to be valid. Provided project leader, for being given the opportunity
scaled charge weight, W/d 3 , falls between of probing into this problem. In addition,
1 and 1000, the analysis procedure recom- Dr. Wilfred E. Baker of SwRI is hereby
mended in this report should be valid, thanked for reviewing this text and making
several helpful suggestions.
The other major restrictions to this analy-
sis procedure is that the standoff distance to REFERENCES
the target must be sufficiently small so that the
weight of an equivalent sphere of air must be 1. J.Sova, "Summary of Armor Materials and
less than 0. 1 times the weight of the explosive Configuration Tests at Aberdeen Proving
products given by Eq. (33), For a 20 lb charge Ground, "Combat Vehicle Mine Protection
buried 3 in. deep, this observation means that Conference (U), U.S. Army Weapons Con-
the standoff distance must be less than 38 in. mand, 28 June 1967 (Confi( -itial Report).
for this approach to be strictly valid, Pro-
bably this standoff distance can be doubled 2. J.K. Cockrell, R. Anderson, et al., "Phase
without causing serious error. Most targets III Parametric Design/Cost Effectiveness
are much closer to the ground than several Study for a Mechanical Infantry Combat
feet. This final restriction does not appear to Vehicle (MICV), "Cornell Aeronautical
be very restrictive. Labs., Report 6M-2144-H-4, 20 February
1968 (Confidential Report).
In this paper a procedure has been devel-
oped for predicting either the specific impulse 3. A. B. Wenzel, R, C. Young, and C. R.
or the total impulse imparted to any target ex- Russell, "Structural Response and Human
posed to a land mine explosion. We have seen Protection From Land Mines (U), "Allison
that the impulse is not primarily an air blast Division of General Motors Corp., Cleve-
phenomenon. The loading iscaused by mo- land Army Tank-Automotive Plant, TR 3481,
mentum in the enplosive soil products. To June 1968 (Secret Report).
calculate the impulse imparted to complea
targets, one must determine a shape factor 4. W. L. Kincheloe, "Reduction of Blast Ef-
which Isa function of the geometry associated fects, " Final Quarterly Report, 0477-01
with encounter conditions. The shape factor (04)FP, Contract DA-44-009-ENG-4780,
for a wheel may be obtained from Fig. 3 and May 1962.
the shape factor for rectangular objects such
as plates may be calculated using Eq. (24). 5. W. E. Baker, "Prediction and Scaling of
These shape factors are substituted into Eq. Reflected Impulse From Strong Blast
(31) to compute the total impulse imparted to Waves, " Int. Jour. Mech. Sci. , 9, pp. 45-51,
any target. The specific impulse at any loca- (1967).
tion on a target may be estimated by substi-
tuting Eq. (33) into Eq. (1) and taking the appro- 6. W. H. Jack, Jr. . and B. F. Armendt. Jr.,
priate component of the resulting specific "Measurements of Normally Reflected
impulse. Shock Parameters From Explosive Charges
Under Simulated High Altitude Conditions,"
ACKNOWLEDGMENTS BRL Report No. 1280, Aberdeen Proving
Ground, Maryland, April, 1965.
This paper is a direct outgrowth of con-
tract DAAK 02-70-C-0579 between the U.S. 7. M. Lutzky, "Explosions in Vacuum,"
Army Mobility Equipment Research and NOLTR 62-19, White Oaks, Maryland,
Development Center, Fort Belvoir, Virginia, November 1962.
and Southwest Research Institute, San Antonio,
Texas. Under the terms of this contract. 8. P. S. Westine. "Explosive Cratering, ",
Southwest Research Institute is to design, Journal of Terramechanics, Vol. 7, No. 2,
develop, and test a new mine clearing roller 1970, pp. 9-19.
system. A rational design of a new roller
system required that we determine the nature
and magnitude of the loads from a land mine
explosion.
107
[I
CIRCULAR CANTILEVER BEAM
ELASTIC RESPO1JF4 TO AN EXPLOSION
110
p(t) a P ff(t), the time function equation CM = H - m if (H-m) is even
becomes
PO = M- m - if (H-m) is odd.
in(t)+ Co )2gnCt) nft 6
n (12) can also be expressed as
m 0,1,2,..,M e
In (12)
The constant colunn vector (Cn) can be expressed
111
Zl4
So far dmping has not been considered. With the consideration of damping, the nth
However, in the actual case there is a slight mode quasi-static response time function,
damping in the response due mainly to internal g e'¢nunt,
-4(t) can be expressed in an infinite
damping and aerial damping. When damping is 9sIe as
considered, the strain becomes
1 2 2--~t
*(t)en nt c0 + t+ C ~ 2+
tc + c' t
ed(xt) P n l End(x) Gnd(t) (22) gn +C ne ne ne
+... (31)
where End(x) ' En(X)/[-(Cn ) 21
In (22), G (t)is the daired nth mode response where
time funct gn which corresponds to Gn(t) in the m (-rnwn)(
1
case of no damping. n cm E 32n
ne j=o (Mj) CndJ 32)
Gnd(t) • 0nd 6t f(,)e-4nns(t-r)in wnd(t.i)d
By taking terms up to the M'th order (!'Q M) in
(23) (31), gnd(t) e-Cnwnt can be expressed in matrix
0
*d( ) form as
[g d(t) - g d(O) cos wdt - - - n t
gd(t) e = (t)T(Cne (33)
2 2 (0) } .non~hen
cd the strain corresponding to g*d(t).
becomes
And '([gnd2(O)J + 2
(
-I Wnd g9nd0) C -1(x,t) ) • NOT} T F Dn)(A) (36)
nd = tan ( n),d tan (l/tanend Cn(x t) = P't) xn '
112
expressing e*(xt) in polynomial form. For this where 4
case, an error will result due to truncation of
the infinite series when taking terms up to the • )n(w) Rf(w) (45)
N'th order only. However, if N' is taken suf- n)2 -?
ficiently large, the error can be minimized.
RESPONSE IN FREQUENCY DONAIN
(46)
Taking Fourier Transform (FT) of (7) with
zero initial conditions gives Taking IFT of a (w)using the real part only as
given in (44) yfelds
Pn PO 2
an W)Mwnfiw)Jf(w) (37) M - n [g~)# (t)] (47)
n n gn (tan)2% wnn n~
Substituting (37) into FT of (5) yields a(t) and g2(t) are time functions obtained from
Px n((w) and d2CW), respectively.
n i)( P(2) (38) n
2
where F1(w) is FT of impulse response function
sin w t of the beam. This impulse response
functdon is causal, i.e. sin w t is zero for g(t) J 0[Xf(w)6(w-tn) - Xf(w)6(w+wn)].
t cO. Thus
Wn j
cos Wt dw
fn (w)= F{sin wnt u(t) "
a (0 ~ = ~ 2 - f(wn) cos nt (9
)
[6(w + ton 6(w- n)] (39) If the response in the frequency domain
n(w) is known, the response in the time domain
where 6(w) is the delta function, can be obtained from (48) and (49). Also, if
the response in the frequency domain is known,
Forcing function f(t) is also zero for M(w)
can be readily obtained because _2(w) has
t < 0. Let FT of f(t) be delta functions only. Knowing (hw),
Af(w) can
be obtained from (45). Then by taking IFT of
No= f + jf (40) Rf(w), the forcing function f(t) can be obtaine.
If the imaginary part of %n(w) rather than the
real part is used, a parallel analysis will give
Substituting (39) and (40) into (37) will give the same result.
an(w) which has real and imaginary parts.
The domain conversion is made possible
FT of a real causal time function h(t)(71 using the Fast Fourier Transform (FFT). lowever
Is FIT does not have the characteristic property of
FT of a causal time function. Thus the respons,
() + jX( )
M(w) (41) likewise the forcing function, in the frequency
domain obtained from the above analysis and that
where obtained using FFT will be different. Actually,
R( NO cos w obtaining the forcing function directly from
(45) is possible In theory only. In practice,
R) h(t) sin wt dt (42)' the quasi-static response can be obtained using
the characteristics of FFT. The forcing functon
can then the
be quasi-static
obtained from response
the derived
Inverse Fourier Transform (IFT) of 11(w) is between and relation
the forc-
h(t) ;2 o. R(w)cos wt dw ing computing
of function. FFT Is simply an efficient method
the Discrete Fourier Transform
2 RMX sin wt dw (DFT) [8).
(43)
As shown in the Appendix, a sinusoidal
As indicated by (43), the causal tine function function (sine and/or cosine) can be represented
can be obtained using only the real or thefrequency domain. This Is
cainary t d u g oy taccomplished by taking the sampling duration as
imaginary part of FT. an integer number of periods corresponding to
The real part of Gn(w).is the frequency of the function. Impulses on the
p n imaginaiy axis correspond to a sine function and
R (w))
a n (c.)J (44) Impulses on the real axis correspond to a cosine
) } )
n {wn n n function.
113
Cutting (or smoothing) thase impulses in The FT of this response was obtained using
the frequency domain corresponds to removing the 1024 sampling points (N) for the duration of
sinusoidal function in the time domain. This eight times the fundmental period (the same
"cutting impulse technique" can thus be used to sampling duration which was used for the compil-
obtain the quasi-static response of the beam. ation of experimental results for theS ft.- 3 in
The sinusoidal response is removed by making use beam). Since no damping was considered, sampi-
of this domain conversion method. Ing was started at t = 0 (with sample intervzl
being AT'. The result is shown in Fig. 3. (FT
A slightly damped sine function can still is shown for positive w only in all Figures).
be represented by impulses at the corresponding In Fig. 3, an impulse at wl on the imaginary
frequency on the imaginary axis, except that axis corresponds to the sine function and an
there will now be a small side lobe across the impulse at wI on the re-al axis corresponds to
impulses. The sampling duration must be the the cosine function. Since the sampling dura-
same as for the case of no damping. Thus, using tion is not an integer multiple of the period
the imaginary part of FT only (refer to the corresponding to w2 , the response in the frc-
Appendix), the cutting impulse technique can quency domain at w2 is not of impulse form.
still be used to remove the damped sine function.
In this case an error due to the existing small
side lobe will result. However, the error is
quite small as shown In the Appendix. 1500 ..........- REAL
MGNRPPRr
TS nd
n_ or T. nd a
'nd : 'nd in (25) / -
o /
- Soo0
In Fig. 2 the response (strain) of an o0 4 0 600
aluminum cantilever beam of length S ft. and FRoUECY Ia fa
diameter 3 in. (which was also tested in a field
explosion) subjected to a uniformly distributed
exponential forcing function is shown. Here -,CC¢
modes up to the 2nd were taken into account and 499 4
damping was neglected. The units of the respore Fig. 3 Response in Frequency Domain of Circular
are psi, which is the result of dividing the Aluminum Cantilever Beam Subjected to Uniformly
strain by n W.l istributed Exponential Forcing Function e-kt
Dn(x)"
(W1 1 4S.5 rad.sec., k=20, L=5 ft,D:3 in.)
114
sinusoidal response corresponding to w2 can also mately l.S%.
be removed if needed.
4 0.0
,.....
..
... .........
.....
.!
10 1 40*',.OU
+ o. o o
co.+,,.,,<,, 1, ..,,,'++I
o€* '¢u i 200 400 60
- -- I"t011.0PO
o
2. Existing damping effect on the modified the field by bolting then to concrete bases.
Four strain gauges (one facing charge one on the
response. Damping effect on the quasi-static opposite side, and the other t4o similarly
response is very small and the difference located, but higher on the beam) were Installed
due to damping was estimated to be only about on each beam. Two strain gauges were used for
0.3% each strain output. The straip gauges were
connected through bridges to r magnetic tape re-
3. The remaining sinusoidal difference cor- corder which was located in . bunker.
responding to 1. This is duo. to the exist-
ing side lobe in the frequency domain when An overpressure rec,,d obtained from a
damping Is considered. The maximum differ- piezo-electric gauge at lhe expected 12 psi
ence is about 2% (refer to the Appendix). nominal overpressure lo~ation Is shown in Fig.6.
From this record it coa be seen that the actual
The similarity between the exponertial peak overpressure is ipproximately 11 psi and
forcing function used and the quasi-static re- the duration is ahoyo. 0.23 sec.
sponse was examined. The quasi-static response
was smaller than the forcing function by approxi-
115
All A
Predicted(rad/sec) Actual(rad/secl
D WI~t~MACi1YIIMA
- - -o..,,,.. C
I
The beam response from the air blast wave 0 o
A t: &a
. G ©-4
- .o
116
either of the two parts of (50) and (51) in each
case. However, differences will be small be-
cause the slope differences at the intersection
h ~ -QASI-STATIC
------ACtuAL os stATIC RsPt
R[SIPONS[ NO1 If-V94~
points are small.
o PRIsIIIR
DAAG-- Actual drag pressures were obtained from
(50) and (51) using (21). In (21), modes up to
A the 6th were taken. Since slight damping
exists in the actual response, (36) should have
S 10
been
modesused.
higherHowever,
than the because
second damping ratios
were very hard of
to
1[ i iobtain and the damping effect was in fact very
small,(21) was used instead of (36). The drag
os /UIIIJ'I1IM pressures obtained are shown in Figs. 8 and 9
ind can be represented by the following:
ocL0 ,-----.--
01 00 003 0
I
~ 00 -; 1-- --
coo0co
5 ft.- 3 in.beam
•IM, IsEco*os 0 - 35 ms Pd(t) = 1.396-24.986t-S87.10t2+
26127.7!"3 (52)
Fig. 9 Responses of Circular Aluminum Canti- 2( 2
+ lever Beam of Length 2.5 ft. and Diameter 2 in. - Pd(t)
535ms 55ms
Pdt) -5.625-404.82t-7990.08t
-5625-074.8t379.t
are lower than the critical Mach number (which Since the response was mainly of fundamental
Is approximately 0.41). Thus with regard to the mode type, the maximum velocity at midsection
steady state values of the drag coefficients, t was obtained using the maximum displacement at
Reynolds number effect on drag coefficients is midsection and the fundamental frequency. Thus,
expected to be dominant in this Mach number r.ngc the maximum velocities were:
5 ft.- 3 in. beam 1.55 ft/sec
Drag coefficients of a circular cylinder
subjected to an air blast wave generally follow 2.5 ft.- 2 in. boam 1.45 ft/sec
steady state values. However, the following The velocities of the beams are very small com-
118
pared to the air particle velocity during the The forner is symmetric about the folding samp-
blast loading (at t a 0, the air particle vel- ling du.ation tf(k z N)and the latter is anti-
ocity was 504.S ft/sec). Therefore, the effect symmetric about tf.
of the velocity of the beams on Reynolds number
of air flow around the beams can be neglected. Expanding R(x(k)) for a point k* (k I N
For both beams, maximum strains were ob- gives
tained near the quarter point of the first cycle R(x(k*)) * [R{(i)cos WI-I{ (i))sin 2 i]
of each fundamental frequency. Analytically, i-0
maximum strains were computed using the loading _N
configuration of Fig. I together with initial 1 N-i N-i 2-k* 21
conditions. For the drag loading, mean drag kEO iN
coefficients of 0.48 and 0.562 for the S ft.- 3
in. beam and the 2.5 ft.-"2 in.beam, respect- 1EI I(x(k)csN-i sk* 2k
ively, were used. Mean drag coefficients are N k-0 icO
the mean values of the drag coefficient for the
duration of a quarter of the first fundamental I N-I N-1 2k* 2k
period of each beam. Computed and measured N kE0 R(x(k)1ii0 sin -- 'sin -W- i
maximum strains are shown in Table 2.
I N-I N-I 2vk* 2%k
Computed(uin/in) N kWO-
k4p I
S ft.-3 in. 'beam 653
(A.3)
2.5 ft.- 2 in.beaml 490 510 N-i 21 k* 2lk•
In (A.3), lie cos -W--i-cos N izO Ni
TABLE 2
Experimental and Computed Maximum Strains 2r Nl
sin 2k- , N 27Tk*. .21Nk i Ni
sin -ix*sn - and t=0sin
The computed maximum strains show close
agreement with the experimental values. Instead 2rk* 2o1
of using mean drag coefficients with dynamic N are analogous to the integral
pressure of Friedlander decay type as was done forms
above, maximum strains aere also obtained using
the actual drag pressures of (52) and (53). The c s
maximum strains so found were very close to the 0o 2"k*y'cos 2wky dy, Icos 2wk*sin 2wky di
strains determined using mean drag coefficients. s 2ak*ysin 2wky dy and 1sin 2aky.cos 2xkydy
The good agreement between computed maxi-
mum strains and experimentally obtained maximum respectively. Thus R(x(k*)) becomes
strains implies that the obtained drag coeffic-
ients up to at least 4 of the period of each * N-I
fundamental frequency are correct. Also, the 2-Nk=O (k))}kk**
good agreement suggests that the loading con-
figuration used is justified. I R{x(k)N- 6
+N k=O k** kk**
where 6kk,, is the Kronecker delta
k*i
The Fourier Transform pair for continuous
signals can be written in the form k** Is k* and N - k*
X(W) z .ftx(t) e'jwtdt rk** ' I when k** is k*
2 (A.1l) k
x(t) z f' A(l))Jwtdw = -1 when k** is N -
The analogous Discrete FoL.,ier Transform pair In (A.4), the first term resulted from
to (A.1) is ~ 2nk*
N-1 R{X(i)}cos -r - and the second term from
X(i) j 2lI.)0
N-
NI x(k) 2A kkF i-
1
N~i)
N=- 1 x~k) e N (A.2) N1-
1s2 2ak F ,{Rqi)j
the same
N 2
3
x(k) = i X(I) -N result will be obtained. Therefore, it can be
119
t! 2,
about tf. Thus R(x(k)) :an be obtained using
only the real part or the imaginary part of X(i)
If Rfx(k)) is finite up to tf. The same is true NN
for I(x(k)). 40
120
From the above discussion, it is evident
that when slight damping exists, the sinusoidal
a co-
function must be a sine function only or
sine function only to make use of the cutting
impulse technique; while in the case of no dampA
ing, the sinusoidal function does not necessar-
ily have to be a sine or a cosine function only.
REFERENCES
1. W.E. Baker, W.O. Ewing, Jr., J.W. Hanna,
and G.W. Burnewitch, "The Elastic and
Plastic Response of Cantilevers to Air
Blast Loading", Proceedings of the second
U.S. National Congress of Applied Mechanics,
1962, pp. 853-866.
ACKNOWLEDGMENTS
The work reported herein was supported by
the Defence Research Board (Grant No. 1678-09)
and the National Research Council (Grant No.
A-3384). The authors wish to thank Defence
Research Establishment Suffield personnel for
their helpful suggestions and for the oppor-
tunity to participate in event Dial Pack (500
ton TNT field explosion).
121
j MEASUREMENT OF IMPULSE FROM SCALED BURIED EXPLOSIVES
Bruce L. Morris
U.S. Amy Mobility Equipment Research and Development Center
Fort Belvoir, Virginia 22060
I-dimensional anaiysiswais peorfoR-d to determine the physical scaItng
parameters governing the response of wheels to blast loadihg. Hopkinson
scaling was used to determine the proper charge size and location for one-
iquarter scale blast tests. The total energy imported to the test wheels
by the detonation was determined, and the scaled specific impulse was
1clwculated. The test and calculation procedures are descrbed.J
INTRODUCTION IT tLl/2p 1/2 Time scaling
.* The utilization of mine neutralization
hardware require that this equipment operate
under the intense pressure of near-field ex- T2" .L3..
plosive detonations. Since accurate theoretical "0P
knowledge of this explosive-target interaction Soil conditions
is limited, designers have had in the past to 113. cM1/2
resort to full-scale explosive tests to evaluate 3-02P/
and prove their designs. This process is ex-
pensive and time-consuming, so this Center
elected to use scale models to evaluate materi- 7r4- rt
als and configurations for mine clearing roller
wheels capable of withstanding the blast effects Irs
5 9 . Initial conditions
of 30 lbs of explosive. Zp and restraints
/i0=
1 2
IL 1
p1/2 M112
Q
P
Lm=,
I
ditions. These parameters, along with others For replica models, E /7where m and
governing this phenomena, are listed inTable I p denote model and prototype respectively. If
along with their
time (FLT) system.dimensions ina force-length- the same material MM
isused in the model and
prototype wheels, I? = X 3." Ie assume equality
Ten dimensionless products, or Pi terms, of blast pressure, ie,Pm Pp. These con-
I
can be formed from these 13 parameters. There straints are then applied to the above Pi terms
are many techniques for creating this list of to establish the scaling law below.
terms, and no matter which method isused, the
analysis isnot modified as only the algebra is
deleted. Listed below isone set of dimension-
less products or Pi terms.
•; .
A7_
_ ----
_ _ _ _ _ _ - - - '; - " - . i
Table I. Physical Parameters Governing Explosive-Target Interaction
t Time T
2
p Mass density of soil FL- 4T
L Characteristic length L
17 7 Em- AE p scribe1.
f8 m=l a p springs
ykes
kmYl orVeel
fr-0 Tmp 0 cag
Hopkinson has shown that blast pressure is a
function of stand-off distance R and charge
weight W as P=f(R/WI/ 3 ). Thus, if Rm=ARp,
ur= A.JWp to produce equalblast pressures.
124
5l
IMPULSE CALCULATIONS -(M
o0+ -- 1
For computational purposes, the wheel
torsjon-s rtfig system is as shown in Figure 2. 0 + ( K )Q= "( - ) (5)
-
%% i5orig"
9
0max maximum rotation
K spring constant Substituting the boundary conditions into eqs.
(6)and (7)yields
I total impulse
total energy A= Mo
K and B= Ii
mass moment of inertia Sin w tmax and cos w tmax are calculated from
T torque tan w tmax , and the results are substituted into
Ignoring gravitational effects, the total energy eq. (6) to determine Qmx,
absorbed in the test system is given by
125
~!
Westine [2) has developed a method of REFERENCES
calculating the impulse imparted to a target
from a land mine detonation given by [1]. P.S.,Westine, "Explosive Cratering,"
J. Terramichanics, Vol. 7, 11o.2, pp. 9
I= A0 to 19, 1970.
heeII -
where a specifi impulse
total impulse []
[2). by Westine, of
P.S.Detonation "Impulse ImpartedintoSympesium
Land Mines," Targets
I - specific impulse
A - projected area of target 7 in2 for on tn Detection and Neutralization,
these tests March 24-25, 1971, Vol 1 of 2, Fort
0 = shape factor which is a function of Mc Nair, Washington, D.C.
target shape and standoff conditions [3]. W. D. Kennedy, "Explosions and Explosives
Shape factors are determined for the various in Air, " in Effects of Impact and Ex-
configurations tested under this program. These plosions Volume I, Summary Technical
shape factors, together with the total impulse Report, NDRC, Washington, 1946.
as calculated from eq. (4), are used to cal-
culate the specific impulse generated by the
detonations. This impulse is transformed into
scaled specific impulse by dividing by the cube
root of the charge weight. These scaled
specific impulses are presented in Table II
along with the scaled distances and other data
items and are compared to previous extrapolated
data [3] for TNT in Figure 3. The TNT data has
been adjuated to C-4 explosive. Data generated
from these small charges is thus seen to fall
within the limits of data generated using ex-
plosive charges of up to 550 lbs.
Table II. Scaled Specific Impulse
Shot Charge wt. Scaled Dist. Impulse on 1 Wheel 9 Specific Imp Scaled Specific
Number (lbs) (ft/lb" 3 ) (lb-sec) (psi-sec) Impulse
pst-msec
lbl/3
126
2000 - iii
1000 *
800 4+
600 +
140 -
1 11111II
I lI I III I !
.06 .1 .2 .4 .6 1 2 4 6 10 20 40
3
SCALED DISTANCE (FT/LB1/ )
127
DYNAMIC ANALYSIS
THE EFFECTS OF MOMENTUM WHEELS ON THE FREQUENCY RESPONSE
CHARACTERISTICS OF LARGE FLEXIBLE STRUCTURES
INTRODUCTION
The generally accepted method for deter- equations of notion only when modal damping or
raining transfer functions on large systems is damping proportional to either the mass or stiff-
based on calculating the normal modes of the ness matrix is assumed. The introduction of the
system and using them to uncouple the equations momentum wheel terms, however, recouples the
of motion. This method is based on the assump- modal damping matrix, and a second eigenvalue
tion that the modal damping matrix is diagonal. s3lution must be determined. This is done using
There are times, however, when the assumntion of the state vector form of the equations of motion.
a diagonal damping matrix is not valid; for ex-
ample, when a structure has gyros or momentum This meLhod was used to calculate transfer
wheels. The introduction of the momentum term functions for the Skylab Apollo telcscope mount
produces a skew-symmetric, rather than a diago- (ATM) in order to determine the effects of the
nal, damping matrix. In the past, analyses of control moment gyros (CMGs) on the structure's
this type have neglected the effects of the mo- behavior. The importance of this example and
mentum wheels. The purpose of this study is to the reason for choosing it lie in the fact that
include the momentum term and develop a method these same transfer functions are used In de-
to facilitate the calculation of the transfer signing the experiment poitting control system
functions so that their effects on the structure (EPCS) for the ATM. Because of the proximity of
can be evaluated. the COGs (which are used for overall cluster con-
trol) to the ATM, their effect on the EPCS son-
The method makes use of existing mathemati- sors could be significant. Transfer functions
cal techniques to develop the state vector form with and without the C1Gs were computed. The
of the dynamic equations -i motion. From these results are compared in this paper.
equations, complex eigenvalue/eigenvector trans-
formations are then generated and used to un-
couple the equations of motion. Once the equa- THEORETICAL DEVELO.-NT
tions are uncoupled, they are easily solved for
the harmonic steady-state response. The derivation of techniques used to gener-
ate acceleration responses or transfer functions
Undamped normal modes are used to describe for large flexible structures with momentum
the characteristics of the flexible structure, wheels will be developed in four parts:
These modes are sufficient to uncouple the
3. Equations of motion for a flexible of complex numbers yielding magnitude and phase
structure with spinning momentum wheels information.
or gyros; Substituting Eq. 3 and 4 into Eq. 1 yields,
4. Procedure used In the digital computer when the normal modes are ortho.-ormalized on the
program for uncoupling and solving the mass matrix,
response equations.
(j~
- ~] i 2CUwoRJ) io}
1. Equations of Motion for Large Flexible . 1T o (5)
Structures (0 (5)o
since
Response analyses of large flexible struc-
tures generally employ component modal substi- - ol0
tution or modal coupling [Ref 1] techniques to [meqJ
reduce the number of equations-required for the when the modes are orthoncrmalized on the wass
solution. These techniques rely on the ortho- matrix. Solving for 4qoj
gonality of the undamped normal modes of the
structure or substructure to produce diagonal
equivalent mass and stiffness matrices. In ad- lqol (I2r21 2 +I r2~c-fl1)
dition, modal damping is usually assumed ao that
the final equations are uncoupled and readily 10]T f (6)
solved. T{o")
In terms of the generalized or modal coor- Substituting Eq. 6 into Eq. 2 yields a ma-
dinates, the uncoupled equations are of the form trix equation relating responses to forcing
functions in discrete coordinates:
ml()+ [rw1 { . 2a ](q) -1
q P 0eq eq
- iF(t)) (1) {Xo}= [41 - I + i 2;. o
x x y y
130
mI + - In the previous section, the equations of
yV X X notion in modal coordinatei,'Eq. 6, were un-
z I +6 1y - y 1X (8) coupled and could be-readily solved for re-
sponses. Equation 12, however, is fully coupled
in the velocity coefficients. If one pvocaeds
For an axisymetric rotor (Ix -y = Is) with as in section I ari lets
constant angular velocity X about the z axis, int
Eq. 8 becomes (f(t)) = If t e ;
M. -XI x Ox +
6y ( Iz " 1s) X;
(F(t)) - IFo1 e ift ;
y - -Ox
60 - Ao
; (9) (q~t)) -ilol
z 0. for steady-atate response to a harmonic input,
the solution to Eq. 12 would be in the form
Since the inertias about any axes in the x-y T
plane are the same, the coordinates can be fixed - - 02 + i + []
at the center of mass with the z axis coincident [r10 j onj
with the spin axis and the x and y axes nonrota- (13)
tin&. The angular velocities and accelerations, Cr1 [ ] / fFo}.
6 and Oi' respectively, are then descriptive of The solution then involves invertiog a fully
the coordinate system's notion, while the spin coupled matrix that is also dependent on the
rate, A, relates the body motion to the coor- frequency of the forcing function. Since the
dinate system. inversion must be performed for each forcing
frequency, the determination of responses
from
In matrix form, letting y = I z - Io A, Eq. 13 is very time-consuming and costly. How-
~z ever, these equations can be uncoupled by re-
Eq. 9 becomes writing them in state vector form and generating
complex eigenvalues and eigenvectors. The pro-
1. 0n
00
0
00
L
F 000
( 10) c
4.
t b f
Procedure
e
Used to of
System Equations
ul
Uncouple
Motion the Combined
matrix is in-
The acceleration coafficient
dependent of the spin rate and can be included Rewriting Eq. 12 as
in the system inertia matrix when the momentum
wheels are coupled with the elastic structure. {q) + [B) {) + [E] (q) (F)
( (14)
The velocity coefficient matrix, however, .is
dependent on the spin rate and will couple the where
system modal equations. _
B [ + ]T[] [w;
30 -1I
(q)+ ] [,IT
+2W [,F]) ;) I
131
The modal equations of notion can then be ex- -
pressed as hC -3 ItoI 17 till
!
-(i) -[R]-i [H] {z} -JR]- (D) This formulation of the equations of motion
involves the inversion of a forcing frequency-
or dependent matrix that Is diagonal. If responses t
or transfer functions are to be generated for
JU] (Z) - {7) - (J), (17) more than just a few forcing frequencies, it be-
comes much more efficient to generate the com-
where
[Rh
[0
L-.-
,.0
"'-.-i,
] plex eigenvalues and eigenvectors once and to
invert the diagonal matrix r'o - i0.
than to invert the coupled matrix In Eq. 13 for
each time
-0](H
[UB -B
-0 -------
132
I
FIGURE 1. ASSEMBLY,SO0ARARRAYS
ATh RACK, DEPLOYMENT
Figure 1 shows the ATM rack, deployment The design of the EPCS depends on transfer
assembly, and ATM solar arrays. The rack is the functions calculated for various points on the
main support, and houses the structure contain- ATM system. The goal of this task was to de-
ing the experiments; in addition, it houses the termine the effect of the spinning CHGs on the
CMGs, which are used to stabilize the entire EPCS transfer functions. Although the cluster
Skylab cluster in orbit. These three MGs, control system, of which the CHGs are a part,
which have a spin rate of 9300 rpm, are the mo- and the EPCS are independent control systems,
mentum wheels considered in this analysis. the elastic motion of the CMG support structure
Figure 2 shows the orientation of each of the will cause the wheel (gyro) to generate a torque
three CMGs to the Skylab cluster, that may be picked up by the EPCS sensors and
cause the system to respond.
Figure 3 shows the ATM spar, canister, and
gimbal ring assembly (GRA). The spar is the Using existing models of the various ATM
structure on which the various experiment pack- substructures for the vibration analysis led to
ages are mounted. The GRA is used to aim the a total structural model with 1321 degrees of
various photoeraphic experimenLs. The GRA is freedom. Two factors associated with this study
controlled by the EPCS, which has motion sen- dictated that the size of this model be reduced:
sors on the spar. first, the computer cost to obtain an eigenvalue
TY.
,-/ ,.
eTZ TX+Z
133
solution to this large a model was excessive, To perform the transfer function analysis
and second, the current transfer function pro- on the ATKA force input points and acceleration
gram could only accommodate 50 modes due to output points were selected on the basis of their
computer storage limitations. In light of these effect on the EPCS. The input points selected
two factors, we decided to use the constrained were the GRA flex actuator torque motors and the
component mode substitution method [Ref. 3], or rotor imbalance moments on the CMGs. The output
inertial modal coupling method, as it is more points selected were the fine sun sensor and the
commonly called, to determine the modal proper- EPCS rate gyros (both of which are on the spar),
ties. To use this method, the ATH system was the flex actuators, and the CMGs.
broken into three substructures, consisting of
the spar and GRA, the ATM rack, and the canis-
ter. For this analysis, the deployment assembly RESULTS AND CONCLUSIONS
and ATh solar arrays are considered part cf the
rack. The results shown in Fig, 4 thru 7 indicate
that the CMGs on the ATh have a distinct effect
By using this method, the eigenvalue solu- on the ATM transfer functions. In general, the
tion for the final modes and frequencies was effect is random: that is, the amplitude of the
performed on a 288-degree-of-freedom system, transfer function for a given input/output com-
rather than on a 1321-degree-of-freedom system. bination may remain unchanged. increase, or de-
These 288 degrees of freedom represent component crease, depending on the frequency range in
modes of the various substructures selected with question, Prcdictably, the greatest effect
a frequency cutoff criterion. The frequency of occurs when either the input or output point is
the modes obtained ranged from 0.009 to 68 Hz. near the CMGs, as can be seun in the figures.
*E--5
Despun CMGs
- I-I Spun-up CMS.
:III.; ; IIIII I
Iii II
-u -_-_=
FREQUENCY
(HZI i FREQUENCY IHZ
FIG. 4.-FINE SUNSENSOR
VS. FLEX ACTUATOR FIG. 6- SPARCENTER V1. CMG
Z-6
- . UiI/I -4
13
134
Although the results of this study indicate
that including momentum wheels or gyros in the [Ri -
r ..--
I-
j1.0 I
[ = state-vector velocity co-
calculation of transfer functions for structures
containing them can influence the magnitude of
[1.0J 0 efficient matrix
EI
the the
on transfer functions,
magnitude of theirthe
effect
governing
are not clearly
factors [= -----|0 state-vector
cotaefictoi mti paeent
displacement
understood. Limitations on the scope of this 0 -1.0] coefficient matrix
study prevented investigation of some key fac- -
tors, such as the proximity effect exhibited in -[R] [H] - state-vector characteristic
the example and the effect of higher-mode trun- matrix
cation on the validity of the transfer functions [I] - matrix of complex eigen-
near the truncation frequency. vectors
[a J = diagonal matrix of complex
NOMENCLATUREegnaus
- matrx of udanpednormaleigenvalues~f[~
1711 - upper half of [Y)
matrix
[E] 0 - generalized stiffness
135
I -mass inertia about principal axis
a complex eigenvalue
[ I transpose of matrix
-1
- inverse of matrix
[--] - partitioned matrix
- partitioned vector
REFERENCES
ACKNOWLEDGEMENTS
136
T
i'a
Marvin D. Rhodes
NASA Langley Research Center
Hampton, Virginia
137
'p ~~~~"'
~ ~~ ~ I ~ 4 ~4 C~ ~ ~ ~ .k ~ ~ ~ Lt a-C2JSr~s<~
considered. The rigid body/flexible body Interface
was described by space station acceleration in terms
of Induced inertial loads forcing the flexible arrays;
in like manner the flexible array root forces and
moments acted as forcing functions on the space
station. Structural mode descriptiona of the arrays
were required as Input to the simulation; therefore,
a structural analysis of the elastic system was
synthesis analysis of the array structure was
performed by a stiffness matrix method utilizing
equivalent discrete element structural models
representing 600 Inertial degrees of freedom. The
modes utilized in the simulation were chosen on
the basis of a significant percentage of load
participation in Interface force and moment.
Simulations were performed on the orbiting
t
I
required prior to the performance of this simulation, structural system perturbed by Initial attitude
Provision was made In the simulation for closed errors and external forces representing docking.
loop attitude control system dynamics of the space
station and OCS dynamics for the solar arrays. The To demonstrate the adequacy of the method-
latter control system provided the desired orientation elegy which has been simulated, a problem of
of the arrays with the sun by controlling the rotation known solution was selected --- the uniform free-
about the orbit-adjust and seasonal-adjust ax.. free beam planar response to a unit step load
Outputs of the simulation include interaction forces applied at mid-span. The flexible appendage solar
and moments, magnitudes of all motion variables arrays were represented as cantilever uniform
and control parameters as functions of time. beams, having the first five bending modes as
flexible degrees of freedom. When coupled
The formulated simulation has been applied Inertially with the rigid body translatien mode, the
to an extendible solar array structural concept cantilever mode solution yielded results for freq-
and space station which are presently undergoing ueny and loads which compared favorably with the
separate engineering evaluations. Also, candidate exact free-free beam solution.
array orientation control and space station attitude
control systems have been mathematically described CONSIDERED PARAMIETERS
and digitally programmed for this application. Two The presented analysis and corresponding
attitude control systems were provided for the space simulation Is intended to be applicablo to future
station; they are the reaction jet and control moment space stations with controllable solar arrays such
gyro (CMG) systems. The necessary modal as that shown In Figure 1. Structural concepts
II
138
2,of arrays and space stations as shown and associated liroviding sun aignment within a specified'time
control systems are presently undergoing separate after leaving tie earth's shadow. In addition,
engineering developments without regard to a total the OCS mus, meet accuracy requirements
system dynamics criteria for minimizing inter- despite experienced space station disturbances
action loads. One objective of the present analysis andf provide minimum dynamic excitation to the
is to assist In the development of dynamics arrays. Two generic types of OCS drive systems
criteria for each of the cormponent structures from have been considered in the simuiation and ares
analytical results of total system characteristics. he continuous and non-linear drive system,s.'
The continuous-type drive system employs either
In order to account for all of the significant a DC torque motor or a variable-frequeny
dynamic influences upon the space station and synchronous motor as its drive element. A
solar array load interactions, the following block diagram of the cortnuous-drive-OCS model
parameters were considered as basic and are contained In the simulation is shown In Figure 3.
accounted for in the simulation
Conhns.aio "4colr
139
Z Axis C Outer Gkim Spbas
X Aded ot lomon
Odter
Gimbal
Axis
KKsu)
_°
1+
MO Axisl Sc Attitude Anigle
a Commended
0,; - Attitude Artle: late
IVASr K*, Kh Attitude Argle Galn. Rate Gala
Gimbl P v Control LoopDeadbad
Axis
T a RCS Torque L l
I w Sp-e Station Moment of tnertla
ANALYTICAL FORMULATION
Fig. 4. 3 PM Centrot Moment Gyro Array
Space station and solar array motion equations
were formulated together with the Interactive
of 0. 707. The RCS Is primarily used for dynamics provided by the respective attitude and
reference attitude acquisition maneuvers and the orientation control systems. The method given
momentum desaturation of the CMG system. It by Likino [4] has been used as a basis for the
is also an alternate to the CMG for controlling methodology provided in the simulations. The
attitude of the space station. The RCS is depicted in simulation developed from the employed mathe-
the sketch of Figure 5. It is comprised of 4 matical models has been successfully run on the
sets of quad thrusters providing redundant control CDC 6600 computer. The simulation model
about the pitch, roll and yaw body axes. All utilized for the mathematical system is described
maneuvers using the RCS are performed by firing below.
the thrusters in pairs. The RCS model included
in the simulation is depicted in Figure 6 and * The space station and the two arrays are
the corresponding control equations are given In each modeled as interconnected bodies
Reference 3. with each of the arrays permitted
controlled rotations about the spacecraft
In order to account for space station disturb- attachment points. The allowed axes
arices such as crew motions and docking forces, of rotation consist of those parallel to the
provision has also been made In the simulation to space station roll axis and the array
allow for the application of time dependent forcing vane axis.
functions.
* The flexibility of the solar array Is
modeled by means of a truncated set of
+yau Z cantilever modes which Is excited by
the acceleration of the array support.
+ A difference equation technique [5J is
utilized to obtain the modal response.
140
The array driver gear train for the. axis The following notation applies
parallel to the roll axis Is modeled as
an ideal mechanical transformer. The J - lndex.of solar array, J is equal
vare driver axis is directly driven, and to 1 or 2
either or both motions about these axes
may be rigidly constrained. FAj TA " qare
the transient forces and
J torques produced by the flexible
9 The simulation orbit generator uses array dynamics
Lyddane's method (6) for near earth
orbits. The generator is included In T - torque exerted on space station
the simulation to provide a reference for by hinged body along constrained
guidarnco commands. A block diagram axes
representation of the simulation program
is presented in Figure 7 where Impor- mT - total system mass
tant logical switches and function Inter- Tr
connection have been clearly delineated. 0 - Newtonian reference point
141
. n-
, ~r~~i&r
- ~ Va
~ ~ z.~ 's4
.t .~ ra~n~z n~z ~ MudS~aas 3I~....Eg ai!
-777"N
IC
Miajor Cycle Furtiona i
o Appendage Equation Update
o space Station Guidance '113 it
o Solar Array Guidance
Intermediate Step 2 0mi
411. 12 i
r 3 Vwherec l, w12, W13 are the rotational rates about
Dynamica Equationa the ith coordinate frame axis. The rigid body
o CNIG/teaction Jets
o Appendage ynamics scalar equations derived from those presented
o Rigid Bodies Dynamics above reduce to the matrix form shown in Figure
9. The submatrfees (Aij represent the linear
Blk 4 i term coefficients, Vol and ol0represent the rigid
Inisgration Packsage 1ody translational and rotational accelerations
respectively and LAI the rigid array rotational
accelerations relatiA to the space station
Fig. 7. Simulation Flow Chart coordinates for the two unconstrained axes of
rotation. The right hand side of the equations
represent the applied forces, torques (control
Space sttion torques Included) and all non-linear terms.
Controlling torque profiles are computed in the
- CG simulation at designated time increments by
0S Space the appropriate control equations and are used to
SeCO staion force the above matrix equations. The array and
142
V V.
LPVZVL
- t>V-, ~ t.I'V~t4I
- -V
-- -- -- -- M is the mass matrix
analysis.
I'lexillet Arrav
+ 7 The assumption has been made in going from
Pl.~l I 4r Ir Equation 6 to Equation 7 that the motion dependent
i
If- - e'ctot i,.hn of IT t I , matrices which are functions of rigid body rotation
mat particle of flexlii h ass 1'oin rates, are small and have a negligible effect upon
asv,;r pirn. the resulting transformation procedure. Without
corsion o T i t this assumption the simulation would be required
falp pnrice for urml.,rmed to he performed in discrete coupled coordinates
aippervIlsge
- t ior ileflhtlion of iast with resulting manipulations of large order
porticle matrices.
Fig. 10. Flexible Array Geometry The left hand sides of Equations 4 and 7 are
constructed in the simulation from computed
direction cosines, the rigid body inertia tensor,
Substitution of the appropriate direction cosine center of gravity and appendage attachment locations
matrices and consilderailon of the appropriate In the space station coordinate frame and modal
properties resulting from elastic deformation properties of flexible appendages. The latter
gives the following: includes deflection coefficients, frequencies,
damping coefficients and masses, for a chosen
[MNl*q +[K~q -[G~q -[11)4 + 1, (6) nunlber of modes.
where q "
ri21 1
3
2
1
2
2U
2
U1
N Rloth the rigid body andi elastic equations
are solved sequentially employing a finite
N N] difference method. A change to the Integration
U2 U procedure in the simulation is presently being
1
143
4.4
made and incorporates the simultaneous solution exercised independently and the results correlated
of the equations. Excitation of Equation 7 is with known data. At that point, the complete
accomplished by the internal loading acting on each program was checked for continuity and a problem
of the discrete masses resulting from the trans- of known solution was then executed to verify
lational ani.rotational accelerations of the rigid the structural dynamics methodology contained in
bodies. The rigid body equations are in tirn the simulation. The problem ---a free-free beam
forced by the application of the flexible body with zero damping subjected to a concentrated
interaction loads, control system torques and force at mid-span ---was selected because it
external forces. The simulation computes was considered to be a good test for solution
interaction loads from the following definitions, convergence. It also provided information
concerning the accuracy of the analytical approach
F = m Transient and programming techniques. Closed form solu-
J
A 1. 71 force tions for the modal response of free-free beams
subjected to concentrated forces are provided by
T r in O " Transient Leonard (8J. In addition, two other solutions
A [ moment for the response of a free-free beam were obtained
using numerical integration methods. One
solution was obtained by the method provided in
defined by aReference 5 and the other by an independently
derived method using a variable order Adams
F -F + C F integrator. These solutions provided a basis
INTj II
J A
Aj for verification of the simulations.
144
MEQ 1 + M 2 derived and are shown below. The subscript "F"
1(1) EQ (a refers to free-free in these equations.
i I, 2,3... N Fi
2#1 it) 0 F
F EQF 2 o (15)
and
F0 applied step load at mid-span It should be noted that anti-symmetric modes of the
free-free beam were not excited due to the
X(t) = coordinate for rigid body motion positioning of the disturbance force at mid-span.
hEQ = generalized mass of the ITil Modal data for both the free-free and canti-
I cantilever mode lever uniform beams were obtained from stzndArd
reference tables (9] in order to compare the
o = iT1 cantilever mode natural cantilever and free-free beam formulations.
frequency (uncoupled) Numerical data assumed for this comparison were:
MA = 5,0 slugs, Fo = 0. 5 lb. I 1 = 12. 566 rad/sec.
(t) = generalized modal coordinate The evaluation of the numerical coefficients of
for the ITI1 cantilever mode Equations 8-16 was based upon a 25 point
discretized mass representation of both the canti-
= coupled frequency of O.bratior. lever and free-free beamt. The frequencies
of the system, approxlmates the obtained by an orthogonal coordinate transformation
iTil free-free frequency analysis of Equations 8 and 9, for the cantilever
beams, are given in Table 1. These frequencies,
The solutions of Equations 8 and 9, using only along with corresponding free-free beam
the fundamental cantilever mode, are given by the frequencies have been normalized with respect to
following: the frequency of the fundamental free-free beam
mode. The rate of convergence in this frequency
F2 comparison is demonstrated by the successive
X (t) " - + _ cos g t (10) number of modes used. Similarly, the degree of
A A ( A Q1 I
Vcorrelation and convergence in the modal
~amplitude domain can be seer. in Figure 12
which is based upon use of three cantilever modes.
F The comparlson, show that the cantilever beam
-2 formulation can be used to accurately descrbr
(AMEQ -l ) free-free beam modal properties If a sufficient
numiber of cantilever modes are usedl.
145
T'ABLE 1
F'requency Comparlisan of (Uniform Beam) Cantllor~r + Rigid Body Mode Representation
of a Free-Free Uniform Beam
Frequency Ratios: f ,
n I Jrei-free
1.03
'L4o
Sm.0.t abyrn-rI ______.__._ _1_If__1
.6-ree "s
3. ..
FilM. 13Mise313t
Timel.V
.. . 0
146
which Is the method used in the simulation. .3
The comparisons given by Figure 13 show an .2
excellent agreement between the simulation results .1 . 2 .0 .06
and those obtained by an independent method
0 I.r
utilizing the varlable order Adams numerical -. 1 . ... I-
Integration technique. .3
modes. The low frequency waveform and magni- NI]2 Ro.. ts.uise frst rl,,c,-tlvr wed"
tude in both histories is seen to be in good agree- .. S
' ~ ~ ~ ~ l 06f/o,
.08 .,o fIt.,, .,6 .1, .10 .11
'14-.6-e *Fig. 15. Uniform Beam Comparisons of Shear
... __ . istory @1/4 Span for a Unit Step Force Applied
..4 at Mid-Span
simlaltion Ietutt. teg firr ite Carfiler .ode. are attributable to the degree of convergence in
approximatltug the given number of free-free modes
by the same number of cantilever modes.
l08 .o 0. 12. 14. 16. 18. 20. 22.,
The above presented comparisons between
Time.e cantilever an free-free beam response re-
suits demonstrated that the structural dynamics
methodology contained in the simulation was
sufficient for the accurate evaluation of Inter-
action loads.
-. 8 . I t -Il-M 11/4) 9(1) - Z V 41/4, I )
Fig. 14. Uniform Bean Comparisons of Shear Maximum her at Span x. Nlde
listories @ 1/4 Span for a Unit Step Force Applied 1 Mro
at Mid-Span icw00 Ilekrence 8
90 Figure 9
80 Sim latlon
Figure 15 presents a comparison of simulation antlr
C70 Aroximaion
reflults using five cantilever modes with results
obtained using tho Nastran "Direct Transient 40
Itesponse Method" [5]. The free-free .
30
beam representation was discretized Into 40
masse; for use i this method. Modal truncation 20
was not considered in this method; each of the o
0 2 3
F4 5
discr2te masses was allowed two Inertial degrees
of freedom co.'responding to planar bending of l(ilHfd-Dlmly) symmetric hlode Number. I
the beam. In general, good agreement exists
between the frequency content of the two shear Fig. 16. Comparison of Modal Shear Force
histories. Also, the maximum shear force Participation for a Uniform Beam Subjected
given by the hitorles is in excellent agreement. to a Unit Slop Load Appli.d at Mid-Span
1147
i'
•4
TYPICAL SIMULATION RESULTS 3 degrees of freedom. The stiffness of the array
membrane was considered as a function of applied
An extensive analysis of the parameters that tension loading, and the central boom was modeled
influence the structural dynamics of solar arrays as a beam column. A detailed description of the
is planned. This analysis will form the basis of structural modeling and resultant modal data is
a generalized dynamic design criterion for solar reported In Reference 3. Selected frequencies are
array structures. Ithough this analysis has not listed in Table 2 in terms of in-plane and out-of-
been completed, some preliminary results have plane modes.
been obtained.
t
The selections of modes to be used In the
The solar array/space station conceot being simulation was made on the basis of those contri-
evaluated is shown In Figure 17. The space station biting a large percentage of load participation.
is a rather stiff structure (modeled as a rigid Load participation for symmetric modes is
body) which contains 96. 5%of the total mass. evaluated by calculating the shear at the array
The solar array is a large flexible structure attachment point due to a base translational
composed of membrane strips stretched between acceleration. The participation is equal to
the Inner and outer structural support members. (Z mi 0i) 2 /Rn Zmi where mI is the discrete mass,
An extendible boom in the center of each vane 0i is the modal deflection coefficient and Rn
applies a tension load to the membrane strips. is the modal mass. Load participation for anti-
The solar cells and associated interconnects are symmetric modes is evaluated by calculating
cemented to the membrane substrate and generate .the moment at the array attachment point due to
the power required for space station operations, a base rotational acceleration.2 This participation
This array was designed Iy the Lockheed Company function is equal to (Zmo i ri) AZn mr12
under Contract NAS9-11039 and is reported In where r[ is the distance from the vane axis of
Reference 11. the array to the mass point. A typical mode shape
seais shown in Figure 18 for an out-of-plane anti-
Outter Strxtursl SRtdst) symmetric mode. The load participation factor
Support for this mode was 67. 7%. The results of the
. 2! analysis also Indicate that modes with a significant
load participation ;all within the frequency band-
width of the space station attitude control
system. Therefore, coupling of the attitude
.. control system with array modes can be expected.
'y 7
148
TABLE 2
Frequencies and Modal Participation Factors of Selected Roilup Array Modes
Out of Plane
_______ Symmetric Anti symmetric
In-Plane
Symmetric Antisymmetric
2500
1000Foc
Force
1000
I~~
- Accelerationx
- 500
0 0.5 .0 1.5 4 6 0 10
Timne (Second@)
"NC IFig. 19 (a). Simulation Rlesults for an Externally
Applied Docking Force
149
in initial parameters and basic structural data, interaction Force ad Moment Histories
will provide the basis for ihe derivition of a 200
otructural design critetion.
.1 Solar Array Anglar Error Hitorlos
-
0
-. 05 - --
02 4 6 10 10
Time (Seconds)
-3L
0 2 4 6 a 1O
Time (Secors)
150
REFERENCES 7. "Structural Interaction Simulation System",
Technical Report R104, Revision 1, Wolf
i. P. W. Likins and H. K. Bouvier, "Attitude Research and Development Corporation, River-
Control of Non-rigid Spacecraft," J. Astronautics dale, Maryland.
and Aeronautics, May 1971/Vol. 9 No. 5 Pg.
64-71. 8. R. W. Leonard, "On Solutions for the Transient
Response of Beams", NASA Technical Report
2. "Preliminary Synthesis and Simulation of the R-21, 1959.
Selected CMG Attitude Control System," General
Electric'Report EL-506-D, 5 March 1970, 9. D. Young and R. P. Felwar. Jr., "Tables of
General Electric Company, Binghamton, New Characteristic Functions Representing Normal
York. Modes of-Vibrition of A Beam," University of
Texas Publication No. 4913, 1 July 19)49.
3. "Interim Report, The Study of Dynamic Inter-
actions of Solar Arrays with Space Stations and 10. R. L. Bisplinghoff and H. Ashley, Principles of
Development of Array Structural Requirements" Aeroelasticity, pp. 344-350. John Wiley and
Fairchild Industries Report 8581R-1, February Sons, Inc., 1962.
1971, Fairchild Industries, Germantown,
Maryland. 11. "Evaluation of Space Station Solar Array
Technology and Recommended Advanced Develop-
4. P. W. Likins, "Dynamics and Control of ment Program,", First Topical Report LMSC-
Flexible Space Vehicles" Jet Propulsion Labor- A981486, Lockheed Missiles & Space Company,
atory Technical Report 32-1329, Revision 1, December 1970.
Januay 15, 1970.
12. "Stardyne User's Manual", Mechanics Research
5. "The NASTIRAN Theoretical Manual", NASA Inc., Document, Los Angeles, California,
SP-221, Section 11. 3, 1970, Office of Technology January 1971.
Utilization, NASA, Washington, D. C. %
DISCUSSION
Mr. Zudans (Franklin Institute): When you indi- Mr. Mains (Washington University): You said at
cated on one of these diagrams the docking force in- one step that you were using a 600 degree-of-freedom
troduced into the system, how was it introduced rel- system for analysis. Is that correct?
ative to the mass-center of the system? You seem to
be getting the moments, but there were no indications Mr. Weinberger: That is correct.
on the slide as to how this was done.
Mr. Weinberger: In this particular case the Mr. Mains: Did you use a direct integration
docking force was applied at the aft end of the space technique to get those response curves that you
station along the axis of the space station. There showed?
were some out-of-plane motions. These motions
were due primarily to a slight misalignment of the
solar arrays that were active with a linear control
system. This means that even 0.1 of a degree atti- Mr. Weinberger: Yes we did.
tude error would cause some motion of the solar ar-
rays, and hence some torques Into the system. But
the docking force was applied along the axis of the Mr. Mains: How do you have any handle on the
space station. There was no eccentric force or load- meaning, the reliability, of a solution of that size for
Ing. this kind of problem?
151
DISCUSSION
Mr. Weinberger: The 600 degrees of freedom Mr. Weinberger: We are examining this problem
simply referred to the structural analysis model that in connection -with the tension loads that are trans-
was developed. From the structural analysis model, nritted through the boom and so on. This is one of the
the generalized mass and modal coordinates were areas in which we have made certain assumptions in
used for the modes that were selected. In this case linearizing. The validity of these assumptions has not
we used 12 modes, so the 12 elastic modes were been established at this point, although we do have
used for the solar arrays plus the 6 rigid body de- some preliminary test data from Lockheed on the ar-
grees of freedom for the space station. We did not ray of this type. This was a Lockheed array geometry
have a 600 degrees of freedom model for direct into- that we were studing. From what we have been able to
gration. ascertain, we have chosen the model In such a way
that the results agree fairly well with the test data.
XrMjps Did you then do an elgenvalue solu- We get z correlation in that respect. This is the only
tion on the 600 degrees of freedom? way that we can have any confid ence in the lineariz-
ation of the array.
Mr. Weinberger: Yes. We did It to get the fre-
quencies and the mode shapes and the generalized
masses. Mr. Clevensen (Langley liesearch Center): If I
understood you correctly, you used NASTRAN for
Mr. Mains How do you have any handle on the verifying some of your results?
reliability of an elgenvalue solution of that size?
Mr. Weinberger: That's right.
Mr. Wr¢nberger: I think the reliability of struc-
tural models of that size is fairly well documented, Mr. Clevensen: Why could not you have used
for example, in the NASTRAN program and users NASTRAN exclusively and saved considerable work?
manual.
Mr. Weinberger: Of course one of the problems
Mr. Mains: You might be surprised if you would that we looked at was the coupling of the control sys-
check the orthogonality of the vectors sometime. tem. I am not familiar with the NASTRAN program
with regard to the demap instructions and the auxil-
iary useage of NASTRAN, other than for structural
Mr. Zudans: I would like to comment more on and vibration analysis. I assume that you might be
these questions, because 600 degrees of freedom dy- able to code subroutines which represent the orient-
namically today Is nothing. You can handle 3,000 and ation control system and the rigid body mechanics.
there is a perfcct orthogonality. The NASTRAN pro- One thing that comes to mind immediately is the size
gram and many other programs use the invrse iter- of NASTRAN. We have been able to modularize this
ation routine with spectral shapes and it is. very, very program in such a way that it uses much less digital
reliable. However, I wanted to ask a question. How computer time in core than the NASTRAN program
did you account for obvious nonlinear atti'.udes in would use. We felt it was more efficient in develop-
your solar array? It is such a flexible st,ucture that ing our own program, rather than to resort to the
it could not be handled as a linear one. NASTRAN program.
152
PARAMETRICALLY EXCITED COLUMN WITH HYSTERETIC MATERIAL PROPERTIES
. T. Mozer
IBM Corporation
East Fishkill, New York
and
R. M. Evan-Iwanowski, Professor
Syracuse University
Syracuse, New York
INTRODUCTION
153
thermal conductivity. Such situations may pre-
sent the necessity to use materials whose me-
chanical properties are non-Hookean over-all but
the smallest stress levels. Once we leave the
domain of the Hookean solid elastic column, a
large field opens up even to the most casual ob-
server. The number and complexity of the mater-
ial properties, along with their relatively
crude rheological models becomes overwhelming if -co
one attempts to apply them with rigor to the
problem of the parametrically excited column, co c
(an exhaustive discussion is given in Ref. 1).
Several authors have considered material
properties other than the simple Hookean case in
the analysis of the arametrically excited col-
umn. K. K. Stevens [1), for example, solves for
the cases of the Maxwell Element and the Three-
Parameter Model as material properties for the Fig. 2 - DaVidenkov's Model
stationary case. The case of a simple viscous
damper in parallel-with a spring (Voigt Element) cross section at any time will be developed from
is shown in references treated for stationary the stress distributions. On other instances
and nonstationary modes. Later K. K. Stevens expressions for the temporal part of the motion
and R. M. Evan-lwanows'd [2] introduced the corn- of the column are derived.
plex modulus material property representation to
the stationary response of the parametrically Analysis of the system configuration yields
excited column. It is important to note here the partial differential equation
that in this analysis energy dissipation occurs
only due to the bending of the column and not a2 t 2 t a2 U
due to its axial compression. V. V. Bolotin (3] x JR) + Pit)2u+m...ZO (1.2)
considers amplitude dependent damping forces due a xa x a t-
to viscous and dry friction at an end support of
the column, but these forces do not arise from where M is the moment at any cross section and
the material properties of the column itself. M = fA a y d A. We have made the usual Bernoul-
Detailed discussion is also given by Mozer (4]. li-Euler assumptions, and the deflection u is re-
garded as small.
This paper deals with determination of re-
gions of stability and instability as well as The strain dt any point is
the lateral amplitude response of the axially 22 U
excited column whose material properties are of C = C + Y2 c + y K = E + y 2 (1.3)
the pointed hysteretic loop variety. The Davx- a X
denkov model is used to represent this material
property. So far the treatment is for any material proper-
154
-j
the same as the period of P(t).
0
AlI
155
E<< 1Eand
where iswhere
a small parameter such that E > 0,
her)
(lnd 2 i
2n-l(a u)n yn+l d A (2.7)
(a,
X) x2 max
aX
is the dissipation functional depending only on Comparing (2.7) with (2.1) we deduce that
the curvature and material constants. The ar-
rows above t represent the branches correspond- E f [( _U a2 u n -
- P(t)f(t)!gin x+ m i(t) Sin tx- 0 (2.4) where fm(t) is the maximum value of f(t). It
can be shown that the expression + Ey/n{ )yn(
Inorder to make any headway toward solving in (2.9) is always odd about y = 0 regardless of
the value of n. Thus the total integrand of
(2.4) we need to know more about the functional (2.9) is always even. Due to this property we
,. We rewrite the equation for Davidenkov's may integrate the above expression from zero to
-model h/2 and multiply by two instead of integrating
from -h/2 to +h/2. This removes the necessity
-[(co
- )n 2n-I - on)) (2.5) to keep track of the signs of y and c. Thus
n4 h/2 2
__2____ -
n'
and the definition for the moment a X2 ,= ; 2 EYiW"0 2 L2 n [(n-1) Sinn
". ydA ax2 L
and also noting that .nSin n ix (I)n+2 [(-f(t)m f(t)n
a2 u
Y -y . 2n-l (.
1 )n f(t)nm)) y
yn+l d
dy
ax
in the case where the axial strains are neglec- where W indicates the width of the column. In-
tegrating the above, we obtain
ted. We may now substitute the above value for
c from (2.5) and substitute the resulting ex-
pression for the stress into M we get a2 2 E W n+44 2
C 322 2 WnIT ( [(n-1) Sinn n -
ax 1. 1?
2
E 2y ~ U[,2U +
a x, d A + x2 max - n Sinn x [(-f(t)m ;f)) n
- 2
+ 3 u)n
2
yn+l d A (2.6) f(t)nmn] (h)n n (2.10)
ax a X max
We consider the special case where n 2
2 2
=E I -u -:EX [(,2 2U
)M +. u n
a x 2 n a x max2 - =2
ax
156
4form
2 m 6
- f(t)m]} (2.1l) where Es 96L mn
96
Equation (2.11) represents the distributed lat- and the term 2n2p has been replaced by i2n2 p.
eral load due to dissipative material properties. Note that if (2.14) were expanded we would get
We now substitute (2.11) into (2.4) and apply the term + cB 2fmf. Itwould appear at first
Galerkin's method by multiplying the resultant glance thit the coefficient + Z02fm could be
by sin vx/L and integrating, we get combined with n2 on the left hand side of (2.14).
It must be noted, however, that this term is ac-
, L .2 L tually a variable coefficient of f due to the al-
2El 24L f(t) - 2 P(t) f(t) Z- + ternating signs and we are thus justified in
L 1 leaving it on the right hand side of the equa-
+ 2m)it! TEy Wh4s r5 . [ f(t)2 tion. We represent (2.14) in the form
48 L f + n2 f = ;F (f,o,r) (2.15)
]
i2f(t)2 f(t) + f(t) 0 (2.12) where
w2 El
L'2V( We have assumed analysis of the first instabil-
ity region only by using the ang,,lar displace-
is the transverse natural frequency of the col- ment term o/2 in (2.16). The terms a(t) and
umn without axial load and the transverse natu- t(t) are to be determined from the usual rela-
ral frequency of the column loaded by Po is tions,
a = EA1 (r,a, ) (2.17)
157
: 2The region of stability is distinguished from
T 2 S_0 (2 the region of instability by the curve repre-
4 2
-ITO( am -a ) - Sin 2¢u.O (2.21) sented by
V (2.28)
. n B~nam n-- Cos 2p- 0 (2.22) a=n - ).¢v
We may now determine the stationary phbse Stationary response curv for the parame-
angle o from (2.21) a using the value of a from trically excited column whose mwiterial may be
(2.24) represented by the Davidenkov Model with the
special assumption that the axial stresses are
1
[ negligible are plotted in Figs. 6-10. Fig. 6
= Sin- [2A - v and Fig. 7 show the effect of varying the load
2nP(4 f2 + 9) parameter P on the column response, Fig.8 shows
the effect of varying the Davidenkov model pa-
[(2 -f-V2 2 -4
2
V 2+ 9 rameter y, while Fig. 9 shows the change of the
) - V)2 - first instability region with amplitude of re-
sponse. Fig. 10 shows the instability region in
three dimensions.
- P-t] ] (2.25) It is seen from Fig. 6 and Fig. 8 that the
point defined by da/dv = 0 separates the region
for which non-trivial stable solution can exist
The so-called "backbone curve" of the sta- to the right and the region for which no non-
tionary response is that curve lying halfway be- trivial stable solutions can exist to the left,
tween the two values of a obtained from (2.24). n, u, P0 and y being constant. This means that
It is no non-trivial solutions exist for y/2aless than
2___f?(2ov - Eywh4V 6 that at which da/dv = 0. This result has not
a, Z = E wh 6 been frequently observed in the literature, al-
si(4n2+9) 96L m though it can be shown that the response curve
of a parametrically excited column having non-
The boundaries of instability zone are: linear damping of the form dff 2 where d is the
nonlinear damping coefficient and structural
V [111 2p] nonlinear elasticity does close at the backbone
curve in a similar manner to the present case.
determination of
steps in we
known bnownhes,
Applying nwellsable It is further noted that as long as y is posi-
,saly oteratn heof- tive the backbone curve leans toward decreasing
stable or unstable branches, we obtain the fol- .frequencies, characterizing soft systems.
lowing: For d ao/dv > 0 the solution is stable (Negative y implies that the material property
if is such that it generates energy over a cycle).
Pisarenko [6] arrives at dynamic response curves
8, ao for various problems Using the Davidenkov rela-
- V 2 - 4 - > 0 (2.26) tions which have also soft characteristics, but
beyond this, comparison of results loses most of
and for d ao/dv < 0 the solution is stable if its meaning since Pisarenko analyzes systems
near dynamic resonance, and the present work is
a, ao 2 concerned with systems near parametric reson-
2a - v - 2 - 4 V3 <0 (2.27) ance.
158
It is seen from Fig. 9 that the curve de- second frequency to the response of the model
fining the region of instability for zero ampli- does not undully complic2tesubsequent use in
tude of vibration is the same as that for the the governing equations then such work would
perfectly elastic case, The authors believe constitute a valuable 6ontribution. It may also
that this is due to the fact that the axial be possible to incorporate into.this model'the
stresses were considered small. The coupling of nonlinear elastic effects without difficulty.
the bending and the axial stress probably would Such nonlinear elasticity could be obtained from
have provided a more significant dissipation experimental data on the maximum point locus
term for a 1 0 and thus the curve would have curve.
been shifted slightly to the right from u - 0;
v/2n - 1 and would have been rounded somewhat, REFERENCES
similar to the case where linear velocity de-
pendent damping is included. In the case con- 1. K. K. Stevens, "On the Parametric Excitation
sdered the energy dissipation per cycle is pro- of a Viscoelast'ic Column," A.I.A.A, Journal,
portional to y a3 and so the curves in Fig. 9 Vol. 4, No. 12, Dec. 1966.
defining the region of instability shift to the
right and become rounded with increasing ampli- 2. K. K. Stevens, R. M. Evan-Iwanowski, "Para-
tude. metric Resonance of Viscoelastic Columns,"
•The downward shift in the curves defining int. J. Solids Structures, Vol. 5, pp. 755-
765, 1969.
the region of instability with increasing ampli-
tude may be due to the nature of the maximum 3. V. V. Bolotin, The Dynamic Stability of
point locus curve for the Davidenkov model: Elastic Systems, Holden-Day, Inc., San
Francisco, London, Amterdam, 1964.
159
Radius of curvature C1f~j
O4 Anua dslcmetla
oftVayn
6 - Efftaictt Fig. -n-mltd
Fig. - Cangeof
Intabiity i icreasing
160
DISCUSSION
Mr. Zudans (Franklin Inetitute): You surprised Hertz. It was in this general frequency range that we
me with very poor orthogonality. Was that mainly be- were concerned. We had a number of modes which,
cause of the use of experimental modes wihich hadnot through increased growth of the spacecraft, had grad-
been orthogonalized before their usage? ually crept down into the pogo frequency range. We
were also concerned with some of the large amplitude
Mr. Stahle: The orthogonality referred to the modes like the first longitudinal which actually went
abbreviated model test. The criterion that we had set as high as 55 Hertz.
up was that the measured experimental mode would
check within 10 percent -no' the analytical modes. The Mr. Schrantz (Comsat Labs): Did you couple
problem is the very limited amount of instrumenta- your model with the Thor Delta to check out the
tion used on the solar array panels. I think our main responses?
confidenco was gained from the fact that the frequen-
cies matched up very well, and that the main struc-
tural modes agreed relative to the modal admittance
through the base shear. Mr. Stable: This is done by Douglas personnel.
The model that we have been using is the model that
Mr. Prause (Battelle Institute): What are some I presented here. Essentially it Is a modal model us-
of the important frequencies? We saw a Ioc of natural ing modal coupling techniques to marry the spacecraft
frequencies in the presentation but what rre the con- back to the launch vehicle. It follows the basic iner-
trol system frequencles and what are the pogo fre- tial coupling procedures of component synthesis dis-
quencies for this type of space stations? cussed in the literature to couple this analytical rep-
resentation of the spacecraft back to the luanch ve-
Mr. Stable: This is the Earth Resources Tech- hicle.
nology Satellite which is a fairly small, 2,000 pound,
space craft going up on the Thor Delta. The main po- *The paper was presented and discussed by C. V.
go frequency varies somewhere between 17 and 23 Stahle for the authors.
161
DYNAMIC INTERACTION BETWEEN
Mario Paz
Professor, Civil Engineering Department
University of Louisville
Louisville, Kentucky
and
Oscar Mathis
Design Engineer, Rex Chainbelt Inc.
Louisville, Kentucky
A A A A A
C B B a B B
Oe G G G( G NG
x
'/ D GLOBAL
COORDINATES
164
of the elements are numbered consecu- In general, any member of the
tively. truss, although assumed to be ideally
pin connected at its ends, may undergo,
under the action of inertial forces,
El flexural deformation in each of the two
A 5 principal planes in addition to the
2 A extensional deformations along the Ion-
mgitudinal axis. It is assumed that for
each individual element of the system
-1 4 these three deformations are uncoupled;
thus, the dynamic stiffness matrix for
the element of the truss is obtained
L independently for the two flexural de-
formations in each of the principal
planes and for the axial deformations.
The dynamic stiffness matrix for a
uniform pin-jointed bar element shown
in Fig.5 is obtained by solving the
5 corresponding Bernoulli-Euler differen-
tial equation for flexural deformation
and the wave equation for the axial and
by introducing the appropriate boundary
conditions. The dynamic stiffness
6 matrix
is givenforinthe pin-jointed bar element
the appendix.
2 5
ELATOER -4
1~
9 L
EXAMPLE
A dynamic booster conveyor support-
5 ed by a truss type bridge between two
64 buildings is presented to illustrate
the interaction analysis.
Fig. 6 shows the schematic diagram
2- | of the space truss supporting the con-
3veyor. As explained above, the convey-
or is analyzed initially under the
action of the driving force and the
condition of zero displacements at the
supporting points. Then, the reactive
Fig.4 -Isolation spring element forces are applied in performing the
analysis of the truss to obtain the
first approximation for the displace-
SUPPORTING STRUCTURE ments at the points of support of the
conveyor. In the next cycle these
installations of conveyor support- dynamic displacements are imposed as
insionaltres anr complex three external actions on the conveyor. The
di m inple
dinbui
rfuefls or
a ndyst omtr asim
sses to exs
ums. m-- f irst four cycles of the interacting
supporting structure
ple building floor systems. A space effects between are shownandinthe
the conveyor
truss is used in the interaction exam- Table I.
ple presented.
165
Ph, -------- 7: ,-A j -, -t- 4-1--
*I
BWOTTO DIAGONALS
Plan View
BOOSTERSRO-
Elevation View
TABLE I
Interaction Results
(a) Amplitude of truss vibration at conveyor suports fin.)
Suppot Support 2 support suport4
Cycle Supo 3
-t - up
1 y lorizontal
0 Vertical
0 Horizontal
0 Vertical
0 - Horizontal
0 Vertical
0 Horizontal
0
Vertical
0
2 -0.0077 -0.0370 0.0022 -0.0830 0.0082 -0.0998 0.0389 -0.0480
3 -0.0365 -0.0710 -0.0198 -0.1560 -0.0045 -0.1890 0.0600 -0.0770
4 -0.0390 -0.0730 -0.0200 -0.i590 -0.0053 -0.1940 0.0620 -0.0800
5 -0.0390 -0.0730 -0.0200 -0.1580 -0.0052 -0.1940 0.0620 -0.0800
(b) Amplitude of conveyor reactive forces at supports (lbs.)
Support I
Cycle iorizontal Support 2 Support 3 Support 4
Vertical Horizontal Vertical Horizontal Vertical Rorizonta Vertical
1 339 -172 203 -T6 175 -33 139
2 153 -123 162 -111 178 -62 204
3 133 -132 153 -110 173 -48 237
4 136 -132 153 -110 173 j -48 239
166
*'4
NOTES ON COMPUTATIONAL METHOD CONCLUSIONS
A computer program in Fortran IV An iterative method for the analy-
is developed for the analysis of the sis for vibrating conveyors mounted on
space truss, the dynamic booster con- supporting structure has been presented.
veyor, and the initerative procedure The two systems are analyzed separately
as described in this paper. The flow using as boundary conditions the de-
diagram of the computer program follows. flections and forces developed at the
points where the conveyor is supported
Flow Diagram by the structure.
Machine-Structure Interaction
Several testing problems as well
as actual cases of installation of con-
START veyors supported by structural systems
were analyzed by the iterative method.
1(l) READ AND PRINT INPUT DATA This study of the interaction of a vi-
brating conveyor and the supporting
2) SET DEFLECTION AT CONVEYOR structure indicates that the method
SUPPORT POINTS EQUAL TO ZERO presented requires four to eight cycles
4to converge to the final solution.
(3) LOOP I = 1 CALL
CONVEYOR ANALYSIS PROGRAM APPENDIX
(4) SET REACTION ON TRUSS EQUAL The dynamic stiffness matrices for
TO FORCES DETERMINED IN (3) the basic elements of the booster con-
iveyor and for the pin-jointed bar ele-
(5) CALL TRUSS ANALYSIS PROGRAM] ment of the truss shown in Figs. 2
TO DETERMINE DEFLECTION through 5 may be written as follows:
AT CONVEYOR SUPPORTS
Ia.- Isolation Spring
1 (6) IS THERE A SIGNIFICANT CHANGE
IN DEFLECTIONS AT CONVEYOR SUPPORTSI Kt 0 -Kv -Kt 0 Kv
NO YES 0 Ka 0 0 -Ka 0
vSET NEW -K 0 K Kv 0 KR/2
DEFLECTION VALUES
-K t 0 Kv Kt 0 K
rCALL CONVEYOR PROGRAM TO t v t v
IDETERMINE MOTION AND FORCES 0 -K 0 -Ka 0
a
CALL TRUSS POGRAM TO -K 0 Kv 0 KR
R
DETERMINE MOTION AND FORCES
IPRINT OUTPUT MOTION AND FORCES
~where:
STOP
Ka = Axial spring constant
167
A
E Nb iw2
o DaQ-Mit2
Symmetric
o DS DP-J.W 2
-RA0UDaSaHsa
-EABUb 0 0 EABNb-afW 2
o "Da (S+Hsa) D (Ca-Hca) 0 Da2Q-mIW2
-he: Da (Hca-Ca) D(Ha-Sa) 0.aSaHsa DP-Jf 2
where:
E = Modulus of elasticity
Ha = eosh(aL)
A = Cross-sectional area
Hs sinh(aL)
B =(m2/EA)
D -Ea/(l-Ca ca)
w Angular velocity of forcing
frequency S - SIN(BL)
b
m Mass per unit length =
Cb COS(BL)
L a Length of beam
Ub - COSEC (EL)
Nb = COT(BL) Q aHsa+Salca
M i - Concentrated mass at left end of
beam P Salca-CaHsa
Mf = Concentrated mass
at right end of
a = (mw/EI) 1/4 beam
C s= Concentrated
mass moment of
I Cross-sectional moment of inertia inertia at left end of beam
a = COS(aL)
Jf Concentrated mass moment
of
Sa = SIN(aL) inertia at right end of beam
C2 K c +21 M2
s-
SC(Kc-Ks) SSKe+CCKsmW2
Symmetric
-Cx c - i K, (1+e/l 2)-j
sSC(K-K) CK 1C 2
Cse+Sc
K+S2X
SC(K-K)(SKcC 2 s)S1cl
SC(K _- (SI +C KCS SC(K-Ks) S2 Kc.C2Ks
0 0
sc "ejc/12
c s 0
c(cK8)S 0+I s e2 / 2
where:
C = coso
m = mass of booster
S = sinO
j= mass moment of inertia
about
Kc = Spring modulus-compression Point 0 (iq. 3)
ESpring moduluse.
s =Sprig
moulu-she~r ,= , 0 (See Fig. 3)
18
d.- Pin-jointed Bar REFERENCES
E = Modulus of elasticity
L, l,e,h = dimension Mr. Paz: No, sir. The forces coming from the
conveyer to the truss were dynamic forces and the
G = Modulus of elasticity in truss was analyzed as a dynamic problem. The truss
shear members had distributed mass. The paper shows the
dynamic stiffness matrix for a truss which, although
I,yIy, I z = Cross-sectional moment of it is pin-connected at the ends, still has bending due
inertia to the inertia effect.
JiJf = Mass moment of inertia Mr. Zudans: Did you have your truss represent-
with dynamic
systemdegrees degrees of
J'Jxk = Polar moment of inertia
kfreedom and notmass
ed as a lumped only static of freedom?
kaikc,k s = Spring constant
Mr. Paz: Actually this was done by the co-author,
Mimi'm j = Mass but this is not the case presented here. This is the
distributed mass case It has a finite number of de-
grees of freedom because of the masses of the ma-
0 = Angle trix method of structural analysis. But the equation to
determine the dynamic stiffness for each element
w = Angular velocity takes Into account the distributed mass and elasticity.
169
RESPONSE OF A SIMPLY SUPPORTED CIRCULAR PLATE
I
kv 117P- ,
DYNAMIC DISPLACEMENTS
We seek a solution to the governing
equations together with the boundary and
initial conditions in the following
general form
t= j,)T* a,, ,.()
%4 M21 6 (7)
Figure 1 - Typical Circular Plate
in which Rm(r) are the normal functions
arising from the free vibration problem,
Consider the thin circular plate of i.e.;the homogeneous solution. gl(r)
thickness h and radius a, as shown in and g2 (r) are to be selected to satisfy
Fig.l, with its median plane in the r-B the inhomogeneous boundary conditions.
plane and with z denoting the distance fl(t) and f (t)represent the time
from this plane. Within this plate the dependency 8f the boundary conditions (2].
displacement w must satisfy the govern-
ing differential equation 1. Free Vibrations
172
,7
173
Initial Conditions
Making use of the assumed dis- (Ihd, i .. 1
placement form (7) we find that the -1 4))o/'L' T Tr
initial conditions (6) at t-0, for all -,
values of r, become O&
(19)
0 Thus the coeficient may be computed
Te R-'+
according to the equation
The second term of each of these equations
may be expanded in infinite series form,
thereby simplifying these equations: Yr -)fk
Similarly
Am =Em (22)
S* where
Bm EmIf-II L i' L
' r.
The formal solution for w may now be -- IL
written: Finally, using the actual values p=p(t),
9
f =0, f = "Ir/(,-4), gl=l, and g =r0 A(,,L,)
/ (d.)1"z.) I,)
. m) aAd carrying out the integratigns these
become:
,¢ .. . ,-
174
wr
175
early stages of the heating period.
ACKNOWLEDGEMENT
This work has been supported by the
United States Navy under Contracts
N00140-70-C-0019 and N60921-71-C-0197
monitored by the Ntval Applied Science
Laboratory, Brooklyn, New York, and the
Naval Ordnance Laboratory, Silver Spring,
Maryland.
REFERENCES a
I - Boley, B.A. and Weiner, J.g., Theory ,
of Thermal Stresses, Wiley, New"York,
#4 '##*S'-''t'''C. 9
176
iti
1 HiI!~~~~llll
-1IYI
49VINUIX2
fiur 2010 esponse ofT 1/8i ! Thick lat
NYin To ANC71 cal cm Radan Pule
Having ACPeak
InMnit At 0.i
71 cal/cmiRadiant Puls
177I',
A''' PH
tfI
I'~~i 1 il 11i
HavWi A PekintniyAt03le
111
4 Ti! IJil ........ ll! I~......... fil
it H)etr
.1 14 ca/m Ra iant il
Pule 1ni ll
-CC(
-,2t0 -'itM o ~T4C'~ C''
Srs
Ffue6 itiu i t Cne Of /8 i lat !i Atl im 0 0 e
1!!o*n Pulse1
Afe ExpsurAI71 cal/cm Rad
;M wl Ii .
179.
DISCUSSION
Mr. Mercurio (Sperry Gyroscope Company): I Mr. Zudans (Franklin Institute): I would like to
would like to commend you on a very interesting refer to your last slide. Something confused me, and
paper. I noticed that you bad fluxes of 71, and I be- I want to clarify it. You showed the static stressesas
lieve, and 145. being perfectly symmetrical.
Mr. Koch: Yes, these happen to be particular Mr. Koch: No, not quite --- almost.
weapon pulses.
Mr. Zudans: How can that be explained? Was
Mr. Mercurio: What were the associated pres- the heat flux coming to one side of the plate or both
sure loadings? sides of the plate?
Mr. Koch: I think they were 10 psi. They would Mr. Koch: One side of the plate only. It happen-
have wiped out these plates compietely. I had to get ed, I suspect, that the temperature distribution was
rather flexible plates in order to show the dynamics, such that, if you just turned the plate a bit, you get
and I was constrained to use particular pulses by my something that looked symmetrical. It was not actu-
contract, ally, it just looked like it.
Mr. Mercurio: You have not shown the pressure Mr. Zudans: How could the temperature distri-
stresses on the plate. Is that correct? bution be anything near symmetrical If you bad heat
flux from one side only. This is what I do not under-
Mr. Koch: Yes stand.
Mr. Mercurio: Have you done any work in this Mr Koch: If you shifted the distribution curve
area? I would like to know, because some of te com- you would have something that would appear to be
mon materials that we are dealing with get into nearly symmetrical by coincidence. Accidentally this
problems when you apply both the thermal and pres- looked symmetrical, but the actual numerical value
sure loads simultaneously, was such that it was not actually symmetrical.
Mr Koch: We are getting into it right now. In Mr. Yang (University of Maryland): I have two
fact in about a months time I am due to have a report questions. On one of your curves you show various
on this subject. We are also attempting to answer the diameters of 1.5, 2 and 3 meters. Since the rise time
question: "What actually happens if the properties is dependent upon the diameter, I wonder why the rise
vary with temperature, let's say, through the thick- time for a diameter of 3 meters is the slowest, where-
nesis?" The frequency changes and the modulus as the time for 1.5 meters Is in the middle and for 2
changes, so the whole dynamic picture might change. meters is the fastest?
Mr. Koch: I think I may have flashed the curve
Mr. Mercurio: There is also the problem of the too fast. I do not think that actually happened. The 1.5
properties changing over very short time durations meter one was quickest.
where very little irlormation is available.
Mr. Yang: The second question is, lave you
Mr. Koch: I have access to some data which is done any work in the thermal stresses of composite
in a range less than one second -- perhaps in the materials?
millisecond range. I have not looked at it, but I was
promised the data. I agree there is some lack of data Mr, Koch: No, however I hace done some non-
on materials in the short time range. thermal work in the vibration of layered materials.
+ 180
1
WHIRL FLUTTER ANALYSIS OF PROPELLER-NACELLE-PYLON
Yuan-Ning Liu
Naval Ship Research and Development Center
Washingtonr D.C. 20034
181
whirl-flutter phenomenon as
related to
the design of the propeller-rotor k8 Effective blade bending
sys- stiffness constant at
tern for a SEV. Design criteria may be hub-center 1
developed in order to define
a stable
Equivalent mass matrix
propeller-nacelle-pylon system mb Blade mass per unit span-
ae func-
tions of design parameters,
such as
wise length
nacelle mounting stiffness,
nacelle
inerti~s, etc., and blade vibration m Nacelle (include shaft)
mass per unit length
characteristics. The overall design N Number of blades
for SEV must also consider the
vibra-
tion characteristics of other
mcchiner- n Blade designation number
ies and the platform,
(n - 1,2,3, . . . , N)
R Propeller radius
Re(s) Real point of s
II. NOIENCLATJRE
a oa lift
Local r sr Spanwise distance along
curve
blade section slope of a
the blade measured from
hub-center to a blade
ACJ
Equivalent damping
matrix
section
s Laplace transform variable
CdC Sectional blade airfoil
drag and lift coefficients,
t Time
respectively
Blade chord T Kinetic energy
y
c ,c Effective viscous damping U Resultant air velocity on
Ax y constants on nacelle pitch-
ing and yawing motions, a blade element
V Proveller advance speed
respectively
Dissipating fucinx Normalized spanvise dis-4
tance, r/R
D, dL Blade sectional drag an d
lift, respectively 4 41
Xb, s Position vectors for a
E Potential energy blade sectin and a
nacelle se'tion,
h Nacelle length respectively
ea
Instantaneous
I Blade moment of inertia blade sec-
about hub-center tion angle of attack
A'
Im(s) Imaginary part of s ao Initial
at angle
a blade of attack
section
1I11 Nacelle moment of inertid
A including the mass of an Effective flapping angle
pro- i
of the nth blade at its
peller in pitch and yaw,
first bending mode,
respectively
referred to the plane of
rotation
182
~AAA
IA- .~
-~. ~
A-: AAA
A~A __-__________________4_
Nacelle damping ratio For a forward moving propeller-rotor in
183
I.{
FLAPPING
HANGE
. AXISOF LADE
47
• L.6'96AD[ MEJINT
JL
184
Utilizing the Lagrangian equations, the
equations of motions of the dynamic
system shown ia Fig. 2 are given by
dt(~ q aqj C j
Silaerodynamic + ilexternal if +
oill (5
- CCK K
Kd + K d a + K d2a (2
(4b) and can be expressed in the form
d
d0 d1 2
a - Re(s) + jIm(s) (7)
185
_______________ ~ ~ ~ ~ 4 .- -7-T.p r, 771- 7,17- -,.
0.0FORWARD WHIRLING0.
0.00 BACKWARD)
WHIRLING
________0.0
0 0046
~ Is 000
02s1
04 0.1 0.001 .4
Fequency
-. Bending Rati PITCH
FREONCY 0orjir0 u
18j
A digital computer program 112] zero or by omitting time-dependent terms
developed for the Naval Ship Research in Eq. (5). One obtains
and Development Center (NSRDC) was used
to solve Er,. (6).
~187
The propeller-roto'r-system on SR.N4 ,1 one can evaluate the variation of
represents a large high-soced SEV'system, eigenvalues of different configurations
and its structural configuration will and hence define the stability bounda-
probably resemble one of an anticipated ries. The results are shown in Fig. 6.
design. It is believed that design
values for the propeller-nacelle-pylon I
system on a larger and faster vehicle .-- ,
would not vary too much from those of SIATIC VIIGCCOOUNARI$ /
SR.N4. Therefore, a set of typical input CA.-...,._.._.
TABLE 1 0. -- 0
Nacelle and Blade Parameters Used in
Defining Stability Boundaries -
Variable z
N4
h/R 0.83 Fig. 6 - Propelle..-Rotor Whirl Flutter
Stability Boundaries for Various
0.2 Blade Bending Frequencies
" 0.3
'1 3
K 0.0087
The stability boundaries for the
K -0.0216 case of ideally rigid blade and static
d1 divergence were also evaluated and pre-
Kd 0.4 sented in Fig. 6. By examining the
188
'1 -I l yrk -
V. TRANSIENT RESPONSE ANALYSIS VI. COMPARISON OF THEORY WITH AVAILABLE
EXPERIMENTAL INFORMATION ON SURFACE
A physical interpretation of a EFFECT VEHICLES
propeller-rotor configuration for dynamic
stability is the transient response due Since test results and structural
to some initial disturbances. For a sta- information on existing SEV are not
ble system, the amplitude of response readily available, it is very difficult
would damp out in time, but in an unsta- to make any rational precise experimental
ble system the dmplitude increases with evaluation of the applicability of the
time. A special case is the neutral theory to propeller-rotor systems. How-
stable condition, i.e., the response due ever, with the limited structural infor-
to any initial disturbance could be rep- mation available for SR.N2, SR.N3 and
resented by some harmonic functions with SK.5 al- cushion vehicles, the following
constant amplitudes. Numerical evalu- evaluations are made. The propeller-
ations were performed on a PACE nacelle-pylon on SR.N2 is essentially the
Model-2312 active analog computer for same as the one on SR.N3. The values of
five different propeller-rotor config- nacelle pitch and yaw natural frequen-
urations to illustrate the above cdes, and the propeller blade funda-
phenomena. mental beniing frequency onf SR.N2 or
SR.N3 were obtained from Ref. [I]. The
Referring to Fig. 6, use n, 1 to values of the same parameters on SK.5
define the stability boundary and pick were obtci!,ed from a vibration shake
=
five points along the n /,no I line. test performed by NSRDC.*
x y
In practice, this line represents an
The stability boundaries of the
isotropic mounted nacelle. These five
pointspropeller-rotor configurations on SR.N2,
y) SR.N3, and SK.5, shown in Fig. 8, were
0.4, 0.3, 0.2825, 0.25, and 0.2. The obtained from the interpolation of those
first two cases are in the stable shown in Fig. 6. This was based on the
region, the third one is on the neutral assumption that the nondimensional param-
stable boundary, and the last two are in ters used for stability apalysis for
the unstable region. An initial angular SR.N2, SR.N3, and SK.5 were the same as
velocity of 0.01 in the pitch direction those shown in Table 1. This may not
was assumed throughout. The results necessarily be true; however, for the
were expressed through the usage of an purposes of preliminary evaluation, the
X-Y plotter and are shown in Figs. 7a-7e. results should be close e:..ugh to give a
These figures clearly show the backward designer a rough guidance. Nevertheless,
whirling phenomena as well as the degree experimental confirmation of the vklidity
of stability for each configuration. In of the stability boundaries shown in
these plots, the X-axis represents the Fig. 8 as well as those shown in Fig. 6
pitch response and the Y-axis represents is still lacking. The only experimental
the yaw response. Heuce, these figures
actually represent the locus of the * A ]lagcn at NSRDC performed the
motion of a Propeller hub. SK-5 vibration shake test in
March 1971.
189
4C
- 0.
-- - - - - - - - 4 I
$ .1 ~
4-- _ __ ~0
(44
I 14.0 1
I ION SPaY Ma CC
A4 0
__ _ 44 a
-4 u
4w 4.
a 4
v 190
T -and the condition is not serious . . .
STAIC
DYEUN~CII UIS 1 (Ref'. (111, p. 219.)
uSTABLE A
toe ALLY
RIMPKIAN
. ~e . The above quoted statement clearly
NAC PITCHREQUENYRAIO
icurve. As one can see from Fig. 8, the
data point associated with SR.N2 or SR.N3
is very close to the neutral stable line
Fig. 8 - Propeller-Rotor Whirl Flutter
Stability Boundaries for SK.5, SR.N2, and nacelle vibratory motion on SR.N2 and
and SR.N3 Air Cushion Vehicles SR.N3 were observed.
191i
NACELLE FREO'JEOCaz
PITCH q#0.1
MACILL
I YAW
FREQUENCYs
II AA^
Il-
10 to 0 30 4 0 0 i I F -10
0.04
ili
- 0 t0 20 0 0 0 to go
10 to
192
193
important nonClmensional parameters such 6. Niblett, T., "A Graphical Repre-
as nacelle damping, nacelle length, sentation of the Binary Flutter Equa-
tions in Normal Coordinates," Royal
blade inertia, propeller power setting, Aircraft Establishment Tech. Rept. 66001,
Jan 1966
etc.
5. Study the effects on the change of 7. Lytwyn, R. T., "Propeller-Rotor
stability boundaries due to the overall Dynamic Stability," The Boeing Co.,
Vertol Div., Tech. Rept. D8-0L95, 1966
structural dynamic characteristics of
surface effect vehicles. 8. Edenborough, H. K., "Investigation
of Tilt-Rotor VTOL Aircraft Rotor-Pylon
Stability," Journal of Aircraft, Vol. 5,
No. 6, Mar-Apr 1969
DISCUSSION
Mr. Gayman (Jet Propulsion LaboratorX): I ask Mr. Liu: Perhaps I did not make the point very
for a point of clarification. Early in your presenta- clear. We do consider the blade ae an elastic blade
tion you discussed the degrees of freedom you were but restrict it to the first bending mode. I made an
admitting to the problem in reference to blade bend- equlvelant system by considering the blade as rigid,
ing. Did you not mean blade flapping as a ridgid mo- but hinged at the hub. The system also included an
tion? equivalent rotation spring at the hub, and the bending
or flexing frequency was the same frequency as the
Mr. Llu: Yes, I meant the flapping motion of the first bending mode.
blade, but only restricted to the first bending mode.
Mr. Zudans: (rranklin Institute): Because of the
yaw and pitch of the nacelle, the blade plane moves as
a rigid body. Was M'at motion included in your hydro-
Mr. Gayman: That is associated with the oscil- dynamic forces?
lation of the plane ofthe propeller disk, Is it not? The
blades themselves are treated as rigid, are they not? Mr. Liu: Yes.
194
THE DYNAHIC RESPONSE OF STRUCTURES SUBJECTED TO TINE-DEPENDENT
BOUNDARY CONDITIONS USING THE FINITE ELEMENT METHOD
George H. Workman
Battelle, Columbus Laboratories
Columbus, Ohio
Once the general matrix equation has been properly partitioned, then by
standard matrix manipulations, the original mixed boundary value problem
is transformed to a modified force motion problem.
'195
iC,
1
11 J -1 +
toroidal parts. Each part was described by its
major and minor radii and the angular coordi-
nates, tp, so that the exact shape of the bellows
in its neutral position was reproduced quite
Lii ~J accurately. The bellows is of uniform thick-
ness, t - 0.062 inch, and is constructed of
K IK Fsteel sheet -: E a 29 x 106 bin and
(2) (2) poisson's ratio - 0.3.
L;]i J] J
K is
where subscript i is associated with the uncon-
strained displacements and subscript J is
associated with the constrained displacements,. -- '
andIh
Fj = I + [M 3 j + k&W*
[Kj1) [Ai + CKj)IAj1 (4) N'ID
Matrix Equation (4) yields the re- Fig. I - Cross section of bellows used for
active forces at the constrained displacement. vibration analysis (neutral position)
196
I
TABLE I
Natural Frequencies of Bellows Model
4 # :
2197
004~
R4Cor,*.. \h
S, \30
21
204
I 7 7 4 7 24-
4 4 ?
24 0
197
Two different techniques were emp!oy-
ed for the solution of Matrix Equation (3).
to One was the straight numerical integration of
Matrix Equation (3) by a f urth-order single-
step Runge-Kutta method.( 7 The other was the
, modal acceleration method.(5 ) Within the modal
"- ' acceleration method only the first mode was
employed. The static displacement component
~~4 wasn determined utilizing standard static. finite
/3 4 element techniques. By comparing the numerical
:s "0 integration solution and modal acceleration
% /solution with only the first mode employed, it
4
\.14
%
Ithe was clear that, for engineering purposes, only
first vibration mode makes any appreciable
15s contribution to the dynamic solution.
198
Wesley, Reading, Hassachusetts, 1955 2342-2347, 1963.
6. Grafton, P. E., and Strome, D. R., "Analysis 7. Ralston, A., A First Course in Numerical
of Axi-Symmetric Shells by the Direct Stiff- Analysis, pp 191-202, McGraw-Hill, New York,
neas Method", Journal A.l.A.A., Vol. 1, pp 1965.
leo
000
Tim milliseconds
g-0 10 / 39 40 so s
i5
0. I
40-
/ Tkne,nilsecods
"X_' D-ynomlc
-40
Fig. 6
- '!."'' -ia' , -'
..¢.' "".',-'. , , '"-,* **" " '',
Time history of maximum
~d,t "". !'" l" L' ' ''' '' ' ' stress
"'"
at Node 5
'" ... '' *' :Y'rv :/d;t'/"d r
199
ICI
D.-€-
Fig.
7 - Time h~stor fmaiu trs t oe2
Time, miliseconds .
\% \-,-DiNe R1
-% Node 16
-a-Nod 24-
200
DISCUSSION
Mr. Zudans (Franklin Institute): You mentioned Mr. Schrantz (Comsat Laboratory): You said
comparison between modal solutions and you aid that you used conical elements to define the bellows?
the response was entirely in the first mode. Yet the
mode you showed was incompatible with boundrydis- Mr. Workman: Rightl I used a series of conical
placement, so you must have done something in add- elements.
ition to that.
201
I VIBRATION ANALYSIS AND TEST OF THE
This paper presents a unique approach used for the launch vibraticn analysis
of the Earth Resources Technology Satellite (ERTS) and compares the analy-
tical results with experimental measurements. The ERTS is basically a mod-
ification of the Nimbus vehicle with solar arrA7 paddles unchanged. The com-
plex paddle system could best be represented by measured data obtained from
previous Nimbus modal testing. The successful extraction and subsequent re-
coupling on ERTS of the solar array paddle modes from the original Nimbus
experimental mode shapes is given. The analytical model is described along
with its verification by an abbreviated modal test. Good correlation between
test and analysis was evidenced by frequency and mode shape comparisons.
Some areas of discrepancy in the analytical model were uncovered which were
subsequently modified to improve the analytical representation of the spacecraft.
204
so: :.:: .. L AC
ACS
ERTS/ACS
INTERFACE
ASSE4BLT
L/V IhTERFACE
Equation I can be manipulated to derive the stiffness matrix if the modal, spectral and inertial ma-
trices are known (2, 3).
205
2i
[] 1 T [H]J it)I 1#1
gives the identity
!.] []T [H] E [I]
KI K2 XI 0 M O XI
[21 1
X1 ,
K22 X2 [ 22 2
Let = Y4
where: K1 1 is the desired cantilevered paddle stiffness matrix
K21
I 1
K22]
12 ,
-o
[
)2 coordinates being restrained, i.e., spacecraft
] s]
I 1 I [1
[
,T
T
T
2 T
rl
I o
(I
(6)
2'2 ] 0 2J
01
=
F0
N
M2 1
01
M2 j
[1
h
.t2 -;
A
o T
i *1
*
1¢ 1
H
0 A
l
A
1*T~.
2 1 h 0T~ T
2M
2 2
MM2 (7)
T --- (8)
2 H " 2 2 A
[ 1 A] A i 0T i
1 Ml A1 H1 A * 2 [1
122 A 1T MIV2 42 A 02 T 2 0
#
*T N1 * A ~T m
1 *Jy I ht M
1 A (9)~j
solving the above elgenvalue problem yields Y , and transforming
206
'Llak4, 4L.-1-o Y O
... . .. ,-. , , ,= ,.f ........... ,,
-- I',
n ,z
"t I
'NI
I I
I * I
Basically, a simplified Hurty (5,8) attachment primary objective of this test was to obtain
technique was used which entailed free mode fixed-free modal vibration data suitable for
coupling at statically determinate interfaces, evaluating the analytical model of the space- I
The resulting eigenvalue equation form consisted craft for use in the flight loads analysis. This
of a coupled generall -,,d mass matrix and dia- required the measurement of data sufficiently
gonal generalized stiffness matrix with 30 de- detailed to enable the response in the funda-
grees of freedom. Eleven modes of the struc- mental pitch, roll and yaw modes and the lateral
ture and nin, een modes nf the solar array modes in the frequency range of the POGO exci-
were utilized in developing the complete system tation to be evaluated.
modes. The first twenty-four complete space-
craft modes ranging from 13 to 85 hz were sub- Sine sweeps for modal definition were
sequently correlated by vibration testing. performed with low level base excitation at a
rate of 8 minutes per octave. A total of ninety
MODAL TEST DESCRIPTION accelerometer channels were recorded during
the resonant dwells. These were located on
Modal testing of the ERTS full-scale the spacecraft structure in triaxial and biaxial
structural dynamic model was performed as a groups at points corresponding to mass points
4 207
........ ........... ,
- p~il .I
I Illel
aIl l~~* I
4'4
4 J-- I
4 I + 4 1
I A I
I I I
14ACJ4V5W"A.4.
4 208 E
4'V
A constant damping coefficient of g 2 C/Cc
0. 10 was assumed for all modes in the analysis,
I
The extraction technique was successfully dem-
onstrated by the correlation with the modal test. I. Beitch, L., "MASS System - The Computer
Program for General Redundant Structures
The limited quantity of measurements With Vibratory and General Static Loading,
209
TABLE 1.
COMPARISON OF TEST AND CALCULATED NATURAL FREQUENCIES
CALCULATED TEST
NATURAL FREQUENCY NATURAL FREQUENCY
,NU.]BER CPS- MAJOR RESPONSE REGION CPS
a. Frequency Very good Correlation within less than 12% in most modes.
c. Modal Shear Good Primary axis of response less than 16% deviation.
Those in excess lack paddle representation.
210
- - ...- ~ .IVII
General Electric Co., TIS R66FPD172, Sep-
tember 13, 1966.
211
2W,
William F. Hartman
The Johns Hopkins University
Bal timore, Maryland
0.30-
- -
E 9 kg/mme e *
E
E0
0'0
II II A
*0
214
Y4,
0.05
?[ '
I+
ft
i: 0.04 1
,E, Imm
14 kig/mm 2
I2
0.02- /9 kg/mm2
00
O0 ,'"i ,i i
0.01 0.02 Strain 0.04 0.05 0.06
I I I I I
0 0.001 0.002 (StraIn) 2 0.004
Fig. 2. Stress versus strain, o, and squared strain s - from fourth lumbar disc data of
Hirsch[ll]. The near-zero modulus for canine disc[8] isshown as ---.
kind. Furthermore, since neither [3] nor [10] specimens contained half the upper and half the
report the cross-sectional areas of the speci- lower vertebral body and the corresponding
mens, the shift in the curves of Fig. 1 might intervertebral joints. The posterior pedicles
be partially due to the use of inappropriate and facets have been shown to take up approxi-
areas for one or both sets of data. mately 20% of the axial load[l]. Therefore in
calculating the stress on the disc from Hirsch's
Reproducibility and consistency can be data, I use 80% of the load values. The
found inother mechanical properties of spinal longitudinal strain is calculated from the
components. For example, the average vertebral reported lateral bulge of the disc, assuming
fracture stress calculated from the data of uniform incompressible deformation there. This
Perey[7J is .34 kg/mm 2 , which is the same value isa poor assumption for large values of stress
obtained from averaging compressive fracture because Nachemson[13] has shown that the lateral
strengths given in [3]. Of course, itis bulging varies nonlinearly with axial strain and
reasonable to assume that a general material the disc's volume isknown to decrease[3]. How-
description will aoply to vertebrae and discs ever, it should hold approximately for small
only for the averaged data of several specimens values of stress and give at least an estimate
from several bodies. Variations are due to for moderate stress. This p)-codure results in
differences inage, size, sex, disease, injury, the stress-strain curve of Fig. 2. The initial
and specimen preservation. Normal deteriorations nonlinearity is precisely the same as that of
are not yet well documented and certain patholo- Fig. 1. The change in the slope of the stress-
gical disorders certainly go undetected. Never- squared strain plot might correlate with the
theless, the degree of consistency noted above inapplicability of the incompressibility
suggests that gross material characterization assumption. A qualitative correction, based on
isboth sensible and ultimately useful. Nachemson's measurements, would increase the
strains for increasing stress, thereby making
The compressive deformation of interverte- the slope more akin to the initial behavior
bral discs is also reported by Hirsch[l1]. His which agrees with the data of Fig. 1.
215
4 4.
WAVE PROPAGATION INTHE SPINE
A distinguishing feature of the stress- fracture stress -" -"
strain curve for the compression of interverte-
bral discs is its concavity away from the strain 0.30
axis. Such material behavior is typical of many 4
biological materials but here it has interesting
implications regarding stress-wave propagation --
216
stress isapplied sufficiently rapidly, such as
a stress pulse whose length is less than the
average pore size, then the effective cross-
sectional area is reduced and failure will occur
at a lower nominal stress. Thus the shock wave
formation as described here could account for
100 the fact that vertebral fractures occir in-the
thoracic vertebrae during pilot ejection even
[12]. no damage is sustained by the lumbar spine
when
- ~~a
shock input.+(R (]
- [a ll + [ RI - 1 T]
T(4)
The porous structure of the soft cancellous In the above, let [a n * 0.26 and
bon2 of the vertebrae make them inherently c, c = 0.17, 0.05, respectiely. This
weaker in tension than in compression. During orr sornds to a dynamic jump from an tnilial
sufficiently gradual application of compressive stress of .04 kg/mm2 to a valie,.30 kg/mmz.
forces the deforming vertebrae uniformly distri- which is 90% of the fracture stress. The
bute the load effectively over the entire croEs remaining stress is easily accomplished through
section due to the collapse of voids and the the multiple reflections dt the infinite m-
reduction of porosity. However, if compressive pedance. Using the curve of Fig. 4, this gives
217
j
6- 6- Calculated
A Measured (-Hirsch 1955)
E4
Fig. 6. The calculated additional strain increment due to impact is nonlinear with the speed.
The case shown is for a lumbar disc-vertebrae segment prestre-sed 0.005 kgria2 .
an impact speed in the vertebra of 12.2 m/sec, Surg., Vol. 39-A, No. 5, pp. 1135-1164,
which corresponds to a free terrestrial fall Oct., 1957.
from approximately 7.5 meters. Since buttock
and pelvic elasticity have been ignored, this (4] D. C. Keyes and E. L. Compere, "The Normal
value is surprisingly large and certainly does and Pathological Physiology of the Nucleus
not seem to agree with ordinary experience. On Pulposus of the Intervertebral Disc,"
the other hand, impacts which are truly axial J. Bone and Joint Surg., ol. 14, pp. 897-
are seldom achieved and the consideration of 938, 1932.
bending and shear would reduce the injury speed.
(5] S. Fribert, "Low Back and Sciatic Pain
c ile mchspeculative
of what I have discussed will Caused by Intervertebral
be considered because it is based Acta Ortho. Scand., Suppl.Disc
25, Herniation,"
1957.
upon minimal appropriate data, the existence
of finite amplitude shock waves in nonlinear (6] J. E. A. O'Connell, "Protrusions of the
intervertebral discs has been shown to be Lumbar Intervertebral Discs," J. Bone and
plausible and its study is potentially Joint Surgery, Vol. 33-8, pp. 8-30, 1951.
important in understanding dynamically
induced spinal injuries. [7] 0. Perey, "Fracture of the Vertebral End
Plate In the Lumbar Spine," Acta Ortho.
REFERENCES Scand., Suppl. 25, 1957.
(8] Edwin R. Fitzgerald and Alan E. Freeland,
[1) Af. Nachemson, "Lumbar Intradiscal "Viscoelastic Response of intervertebral
Pressure," Acta. Ortho. Scand., Suppl. 43, Disks at Audiofrequencies," to appear in
1960. J. Ned. and Biol. Eng.
(2] M. B. Coventry, "Anatomy of the Interverte- [9] Keith L. Harkolf and Robert F. Steidel,
bral Disk," Clin. Ortho. and Rel. Res., "The Dynamic Characteristics of the Human
Vol. 67, pp. 9-15, 1969. Intervertebral Joint," ASME publication
218
[10) W. J. Virgin, "Experimental Investigations
into the Physical Properties of the
Intervertebral Disc," J. Bone and Joint
Surg., Vol. 33-8, pp.607-611, 1951.
219
219
a't
a'4~I
ACCELERATION RESPONSE OF A BLAST-LOADED PLATE
Lawrence W. Fagel
Bell Telephone Laboratories, Inc.
Whippany, New Jersey
INTRODUCTION MOTIVATION
A commonly used technique for calculating Structures which are designed to withstand
the approximate response of a plate subjected nuclear-weapon effects are often either shell-
to blast loading
equivalenit is to consider the plate
to a one-degree-of-freedom to be
spring- like structures
tective shell-likeorstructurei
are encapsulated within
whibh, for pro-
design
mass system, the frequency of which corre- and analytical purposes, are sometimes re-
sponds to the fundamental frequencyof the plate. garded as composites of plate elements. More-
simply supported plate acted upon by a uni- over, for structural-motion-response studies,
form step-function pressure is analyzed here to it is a usual inherent requirement to consider
determine the validity of this approximation. these plate elements to be externally loaded by
Classical bending equations for plates, that is, blast-induced overpressures* which may in
plane-stress equations, are assumed to be ap- some cases be approximated by a step function
plicable and these lead to a time-dependent, for determining early-time responses. In this
double infinitetrigonometric series solution in regard, a prevalent practice for analysts is to
which each term of the series represents the construct a mathematical model of the structure
transverse vibrational response of a different assuming that an entire panel can be repre-
mode. Comparison of the relative amplitudes sented by a single degree of freedom.t A direct
of each modal response demonstrates that the consequence of this approach is that the calcu-
first mode predominates for the dynamic re- lated peak acceleration response at the center
sponses of displacement and stress; hence the of a square plate element exposed to a step
one-degree-of-freedom approximation is ap- overpressure is
propriate for these quantities. However, the
amplitudes of some of the higher modes of ac-
celeration response are significant compared Peak Peak Overpressure
that these
to those of the first mode, indicatingsolution Acceleration 5Mass per Unit Area of Panel
in-
should not be ignored. Because the
volves an infinite series in time, the maximum (1)
amplitude of response is not obvious from the
solution expression. Acceleration responses
are plotted as functions of time to determine
the contributions associated with frequencies ;Pressures in excess of ambient.
higher than the fundamental. Peak acceleration
response for an undamped plate appears to be tThis practice is recommended in "Design of
about 2.6 times the first mode response i.e., Structures to Resist the Effects of Atomic
what would be calculated using the one-degree- Weapons," U.S. Corps of Engineers Manual
of-freedom analog. EMl110-345-110, 15 March 1957.
Vh w + Eh V4 w =q(x,y,t). (2) ab a 0, J
T
0
g 12(1- 2
v ) (
Therefore,
222
Eqs. (3), (4), and (5) substitute into Eq. (2) to The initial conditions are 0i(0) .
0
form Solving Eq. (8),
i_ (t) sin i
sin 2 2N
g " a b S +W
(9)
Eh
121- (ifr)2 + (ig)2]2 t1i(t) =2(I - Cos coj t)
sin--sinsi a b
0 ~t Eq. (9 susiue intE._3)yild
- 16PoU(t)
Mryy +.U
IL-6I Eh-- ++
2 a b (11)
ij-16 'o Ez +V
with
displacements for evenbeing
all initial conditions values of I Modal
zero. and j the
are trivially equal to zero by Eq. (6a). For odd
values of i and J the equations are in the form Ez 2w 2w
wihaliiilcniiosbigzr.Mdlt
The Laplace transform of Eq. (7) is (l 7), (Lb
~j fr evn vaues f i223
dispacemnts
Salu(0s
for o ) i(+ ijd(S) bi sin irx sin ! y I - co s Wijt) (
are rivallyequl
6a) teoForodd2 tozer byEq.+ L
valus frm
i ad jthea quaionsarein
Ez
-------------------------------------
--------------- W2 - - -
Ez. 2
w where I and j can only be odd values. The rela-
._ ... tive amplitudes of higher-modal values com-
xy = -pared
MY to the respective amplitudes associated
with the fundamental mode are plotted on Fig. 1.
Information from these graphs concurs with the
( 00 O0 philosophy that displacements and stresses re-
2 sulting from mode shapes other than the funda-
1G Po Ez
",2W/g1+v
+and
I mental are negligible
stresses
compared to-displacements
associated with the fundamental
I=1 jl mode. On the other hand, the contributions to
acceleration from some higher modes appear to
be significant. Just how the higher-mode accel-
({I- Cos it) eration amplitudes supplement the acceleration
COS cos jXy . (13) response of the first mode is not obvious from
a b 2 Fig. I nor from Eq. (14); therefore this phe-
a j nomenon will now be more closely examined.
Interest is focused at the center of the plate
The acceleration response is where the single-degree-of-freedom analog had
been presumed to be applicable. At the midpoint
sin Go the acceleration is
16 o a b COS Wijt . 0 0 sin -sin
,r2 6/go16
W/g~ J4-I Wmidpoint 2 W/g L4 i
(14) I=1 1=
In Eqs. (10) through (14) 1 and J are odd nd Cos Wit (15)
#h
ff 3ij\ i + -(. where l and j are oddJ. At t 0, ros wit I1
v 2) \2 2 for all i and J, and Eq. (15) degenerates to
a2
S wi, 00o in
, (first mode) 1 )2 ,
16 PO (g (_)r
Accelerg.tion:
_i 16 P'O IT2 P0
"
w(first mode) iJ' - Wg - =W'
224
~ £
'. 1 7.,7//777/7T. .
. . . . . . . . . . . . . . . . .. . .. . . . . . .
•~~~~~~~~~~~~ " . . . . . . . . .
21 . . . . . . . . . . . . . . . . al . . . . . . . . . . . . . . . .
2$ . . . . . . . . . . . . . . . . 3- . . . . . . ... . . .J . . . .
25S.. . . . . . . . . . . . . . . . 113. . . . . . . . . . . . . . . ..
27 .. . . . . . . . . . . . . . . . T- . .. . . . . . . . . . ..
.S. . .. .. . . . . . . . . . .. .. . . . . .' . . . . . . . . . ..
1
o
• • • • • •I ~O .. . . . . . . . . . . . . . . "
.
IT- 7-.4
5-. .
Ai
, ........ .. . .. . . . . . .........
... . . . . .. . . . . . . . . . . . .
. . . . .. ... . . . .7 . . ' , . . . . .
. . . . I . . . . . . . .
.. . . . . . . . . . .I- .. . . .. . . . . . . . . . . . .
J
1.. . . . . . . . . . . . . . .. . II. . . . . . . . . . . . . . . .
9 ........................................ 1 - j' . . .. . . .
. . . . . . 19 -I . .. . . . . . . . . . . . .
. . . . . . . . . .
I..I ................................................ _j
i2, . . . . . . . . . . . . . .
21 .. . . . . . . . . ...
23. . . . . . . . . . . . . ... . -J . . . . ... . . . . . . ....
225
By-comparing this with the first-mode contri- they are based on damping forces defined'In a
button at t =0 [(16/) Po/(W/g) 1.6 Po/(W/g) manner which preserves the-linear nature of
the modal differential equations of motioi by
which approximately agrees.with Eq. (1)], it ntrodicing a viscous damping term. Damping
appears that the one-mode-only representation is thereby defined so that Eq. (7) is repFaced by
of response Is reasonable and conservative for
the initial value of acceleration. At t> 0 how-
ever, it is not'obvious that the absolute value 2 16' Po (t)
of-the-sum-.of-the contributlonsafromuall middeb
should remain bounded by the amplitude the 4'j + 2(3 wI 4,Ii + W4'=1
;,Iof i Mi-
I Ii 0
i t U~ ,
response of the first mode, and in fact, plots of (17)
normalized acceleration versus normalized
time show a contrary trend. That is, the
equation where 2i w-0 is the modal viscous damping
force,and iJ is the p~rcent-of-critical viscous
m n sin Lr sin Jr damping ratio. When the governing partial
Wmidpoin 6 2 2 differential equation i. definedso that Eq. (17)
Po/(W/g r2 I j results the total response w and its time-
*
P( ) J=l derivatives are summations of the modal values
oij and time-derivatives of OiJ respectively.
cos wijt (I and j are odd) (16) The solution to Eq. (17) in
226
PLATE VIBRATION, 1 MODE CONSIDERED
MI n 1 0 PERCENT DAMPING
S 6.000 -
a.0o
00
.00
~V-
0.000.0.200 0.400 0.600 0.8010 1.000 1,200 1.400 1.600 LOW0 2.000
(TI FUNDAMENTAL PEPIOD)
600
40000 .0 0. .. 0.60 .0 . .0 .0 .0 .0 . .0
(TFUDAETA0PRID
Fig. 2 - oprsnotceeainRsos
Cosierin One and Many Mode
4 2.220
tis i e n e l x ep bn
tI' it It is w i g
INDE[X OF m AND
228
229
-f1
PLATE VIBRATION, 225 MODES CONSIDERED
mx29 n:29 0 PERCENT DAMPING
S6.000
. . . . . . . . . . . . . . . . . .°
-00
0 0.000
Z -4.000 K. .....
....... .
0.000 0.200 0.400 0.600 0.00 1.000 1.200 1.400 1.600 1.600 2.000
(T/FUNDAMENTAL PERIOD)
-. 2 000 -
-...... ...............
inA
: 0.000 AA /
1 - -- I-
0.000 0.200 0.400 0.600 0.800 1.000 1.200 1.400 1.600 1.900 2.000
(T/FUNDAMENTAL PERIOD)
4.00
2.000
S-2.000 .. V* o .
-4.0001 . . . H.... -
-6000 , i
0000 0.200 0.400 0.600 0.800 1.000 1200 1.400 1.600 1.800 2.000
(T/FUNDAMENTAL PERIOD)
230
U )
to NOmsECONSINNRED
9 AIOSS CONSINERI
i i ! I I i I I
0 I 2 3 4 1 6 7 S
PIERCENT OF CRITICAL DAMPIIiN
231
APPENDIX I - REFERENCES
2
APPENDIX II - NOTATION
233
EFFECT OF CORRELATION IN HIGH-INTENSITY
Charles T. Morrow
V , j
Advanced Technology Center, Inc.
Dallas, Texas
(U) A narrow strip (bar or ribbon) is taken as a theoretical test case for
the realism of high-intensity no4 se testing in much the same way that a simple
mechanical resonator is commonly taken as a theoretical test case for more
conventional shock and vibration testing. It is shown that in an effort to
design a realistic test, one must consider tile point-to-point correlation of
the applied field as well as the sound pressure level, even when the damping of
the strip is large enough to prevent significant retu:-n reflections from the
ends. Three types of correlation are investigated in particular--complete cor-
relation at the coincidence angle for the given frequency, and independent ex-
citations at the different antinodal regions. With minor reinterpretation, the
formulas remain applicable when, in addition, an exponential decay of correla-
tion with distance in either direction along the strip is introduced.
(U) At one time, I carried out a simplified (U) That c is of magnitude not greater than
analysis of the infinite strip by assuming that unity msy be proved by squaring the sum and dif-
* excitation took place only at the antinodes, as ference of two random signals A and B, normal-
in Figure 1. This is a convenient simplified ized to unity variance by division by the re-
* model to visualize and to use for inferences in spective standard deviations oA and OB, and
advance of any calculation. For a completely averaging over time
correlated wave at normal incidence, the trans-
verse waves propagating in either direction from
alternate antinodes should be of opposite phase - -
2
and tend to cancel each other. The lower the A 2 2
2 2
damping, the more complete the cancellation. At (A/ I B/OB = A2/O 2AB/OA B + B =
the other extreme, for a completely correlated
wave incident at the coincidence angle so that 2(1 > 0. (2)
its trace velocity equals the velocity of pro- -AB
pagation of transverse waves in the strip, the
transverse waves propagating from the antinodes Hence
should all be directly additive. The lower the
damping, the more antinodes contribute signif-
B i
icantly to the summation. For completely un- IcAB- (3)
correlated excitations exhibiting no dominant
phase angle, the waves propagating from the
antinodes should combine as the square root of (U) In the analysis to follow, an expression
the sum of the squares, producing an inter- analogousof to equation (1) must be derived for
mediate behaviour. the sum an infinite number of random signals,
subject to a phase reversal according to whether
(U) Although the analysis based on this aim- the integer numbering a particular antinode is
plified model does not permit accurate prediction even or odd, and subject to attenuation of the
of the response of the strip to continuously dis- transverse waves on their way to the observation
tributed excitation, it does provide insights point, which will be chosen as the origin.
into the nature and effect of correlation. It
is therefore worth summarizing here as a pre- (U) Let a sinusoidal pressure p ut the i'th
liminary exercise. antinode result in a transverse 4elocity
8
vi = i P"*|
V, pi
236
at, the origin, with a real quantity. Simi-
larly, for the k in antinode, wI = (
Vk = Ok Pk' K5 1 (10)
(U) The square is will merely produce a form more similar to that
of a derivation to come. Equation (6) becomes
(t ) (t )
V f2(t) = IP~ " Pf
p8
2
wf) =K w (f -k
i e-"- IiI+IkI)
pci(f). (12)
Sl
I k PiAf (t) PkAf(t) (T)
(U) For the problem at hand, we will assume preliminary exercise to the analysis of a more
a constant power spectral density wp(f), a realistic model excited continuously as a
function of position. Except in special cases,
phase reversal when the integer i or k is odd, the summations become integrals. As phase
indicating an odd number of half wavelengths shifts in propagation are not limited to 0
spacing from the c;ign, and a transverse wave and 1800, the quad spectrum or imaginary part
propagating in either lirection from each anti- of the cross-power spectrum of the sound pres-
node and decaying exponentially. In short, sure has an effect.
237
(U) In Appendix 1, the derivatior of a gcneral (U) For a completely correlated field incident
equation for the response of the continuous strip at the coincidence angle such that the trace
and expressions for three specific cases is car- velocity equals the velocity of propagation of
ried out in detail. transverse vibration in the strip, the final
# result is
to Equation (8)
(U) The equation corresponding
~is
i ) X2H2W(r) a2+T
cf.
22 ~v p 2
(02+h,2) (8
X(f)
Hwp (f) d a(Igl+lhla)
-2 -1 0 1 2
0 0 0o 0 1 0° I
-2 -1 0 1 2
11- 2yfl
238
i.
terested in ratios of response. It follows CONCLUSION
immediately that
W(f) "2a
(U) The effect of the type of correlation of
vn l-e the correlation of the pressure field on
Rn =- 2 ctructural response is by no means negligible
Wvu(f) 2+3e'a-e
3
, (21) even for infinite structures of typical Q's or
structures long enough so that return reflect-
ions are negligible. Introducing an exponential
is the ratio of normal incidence to uncorrelated decay of correlation along the strip has the
response. The ratio of coincidence to un- same effect on response as increasing the inter-
correlated response is nal attenuation. If the analysis given-here
were extended to cover a strip of finite length,
it would also show a marked difference in re-
W (ff le.2a sponse to the different fields according to
vc • (a (22) whether the number of antinodes is odd or even.
2 2
wVU(f) a2 (a +4V ) 2+3e'o-e'3a (U) It follows that the correlation of the
pressure field as well as its sound pressure
(U) We can find an equivalent Q corresponding level must be considered in establishing a
to the decay constant a in the following way. simulation. As usual, however, realism must
Consider a finite strip starting at g=O and be a compromise limited by practical constaints.
ending at gun, an integer. The return wave When it is not possible to coritr0l the correla-
from the m'th reflection at the far end, after tion closely, it may be desirable to introduce
stoppage of the excitation, is proportional to corrections in the SPL to compensate. Analyses
such as the one given here may be useful in
suggesting the amount of correction to apply.
ACKNOWLEDGMENT
I am indebted to my colleague,
is The time of arrival of the m'th return wave Warren A. Meyer, for verifying the derivation
of Equation (20).
t = mn/f,
so that
mn=ft.
Q = /2u (23)
239
10
20 10 5 0 2
o.1
7 € ". EXCITATION AT NORMAL INCIDENCE
0.1
0 0.2 0.4 0.6 0.8 .
240
- V
APPenIDIX I
ANALYSIS OF AllINFINITE STRIP
+ 2 ft
(U) The analysis to follow will be expendited
e n = cos 21ft + J sin 2nft, (24) by using the complex exponential.
=
2 (25) p P MP (30)
plex multiplier:
e. o
J
j2
ft
e(2tft-0) at the origin. The yXM is taken to be complex,
e e for all possible phase shifts can occur in
propogation Prom various x to the origin, but
it is assumed constant over Af, equal to
cos(2,nft-0) + J sin(21ft-0) Yxfo
(U) The velocity produced by all pxn within
Af is
or
f-I" AvJ(2nf -0)
mAV
mxm
oe 0 J2ift . ej(2,7ft-)
Af
cos(2ft-0) - J sln(2dft-0) (27)
•f x J2nf 2JI x~
e-x ( 2
t~
+Let
J(2nf-0)
0)
cos(2nft-0) ej(2nt- ) + e- J (21ft -
2 (28) "Jexm (33)
Note also that (U) The total complex velocity at the origin
is
d j2nft J2nfr j2 n ft
dt ' (29)
241
/V
F-i1
I * J(2f"-Q-0=)d
j (20 2021 X
f xn xn e dx G G (38)
Af(34)
x=e
X (39)
(U) The total real velocity at the origin Is 0
eyn y (40)
xn Pyn e 0 dy (4
"J" Gxmui
-j (211fro"
Ox0-+xC)]dx,+e-XW E • dy
+x •M (3)A xn Pyn (l
xm~d (35)
P P OS 0-0 yn)
+ sin ( Iy
V f 2 dx f Gx G P PAf
xfm yn xm yn
+ 2f
r-j 2
+ ___ -0y)- As(
f
fe xy mon+ e xw yn xn yn Idy PP
exm yn sInn(mxm-y
n dy h3)
• dx P P
dx G Pxn Pyn
Pxn
2yn
242
e dy
+ eJ [Wpxy(f) - Jqpxy ()]1 1
)k(f xy ) - f - '€
0 (45)
wLPC u(f
PY J .PY COSW(0)dy ,
(5)
Max W (f) Wpx (f) w (f) W f) 1/2 , (19) which is the same as Euation (16).
Cpx(f)
p
*1(f)/
Wpxy
Wpx(fW
(OW y M ,px
f)
Mv/2- A2CWP f ,
ea(Igl+lhl)
" (50)
243
- ,
r
P e-0 d9H~
+f -(.
g+2rgdg
2
)2 2w M
(f
0 Jc(g+h)4
4
-aO) Iffo
f&5h-Jguhdh]
X2H2W (f)
is the
Eqaton(1).p~
which~~5 12n~ .~
poiie2 in 2 2
2 22
A H~w
I (f
''gjIhI)A24 f
-i
we )-ajrhd3 - '-Jdg e~ Boh)r
4o whic istesaea qato 1)
(f)
A2HHVw d e1d U ialcnie
1_gj on il uhta
k4 k+-2 PX
excep
haeefetHs
fo ht 1 (f) w 2
24
X2H2W (f
cos~gj--jhjh~d +2 _c(j~jjhj+J,,j~jjhj
8 gf ed
k- k-
~j - [ r~ii
-, -jh~j)Jv 'j h 2 , ma~ja aS
- 9.chl+jwilhlh
e )22H2W -2a)- 3c -2m +e:c
+f . dik
. , - 2+2e-2t e 2
+j~1
* 2 e +0 +0
2 2 1a 2 2
2 [1- 3 2a
f ) 2H2w (fM
~~2H2w 3
- (f +e~ a
I .1 (. 2 2) l~e-2a
edg dg 2 ejIhIJ
eI dh
kf-ki
2uw k-l which is the same as Equation (20).
2 , 2
2 2
- HW 5 -cg+Jirg
ef~ 2-ah-juh
Mdg f e dl'
,k k4
+ I. fT IF C-ag+J1Tg d +2 -ah-ith dh
2 2 _g e +2~
I a- I
2
2 4- j (k+-i-)
+~~~~[.ky- + e
p +C -A2[J /2 + j k4)]
k-a + C-0 jJ
~2+~.2
12
245