Table Of ContentConcurrent Computational and
Dimensional Analyses of Design of
Vehicle Floor-Plates for Landmine-Blast
Survivability
M. Grujicic & B. A. Cheeseman
Journal of Materials Engineering and
Performance
ISSN 1059-9495
Volume 23
Number 1
J. of Materi Eng and Perform (2014)
23:1-12
DOI 10.1007/s11665-013-0637-5
1 23
Report Documentation Page Form Approved
OMB No. 0704-0188
Public reporting burden for the collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, gathering and
maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this collection of information,
including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington
VA 22202-4302. Respondents should be aware that notwithstanding any other provision of law, no person shall be subject to a penalty for failing to comply with a collection of information if it
does not display a currently valid OMB control number.
1. REPORT DATE 3. DATES COVERED
2014 2. REPORT TYPE 00-00-2014 to 00-00-2014
4. TITLE AND SUBTITLE 5a. CONTRACT NUMBER
Concurrent Computational and Dimensional Analyses of Design of
5b. GRANT NUMBER
Vehicle Floor-Plates for Landmine-Blast Survivability
5c. PROGRAM ELEMENT NUMBER
6. AUTHOR(S) 5d. PROJECT NUMBER
5e. TASK NUMBER
5f. WORK UNIT NUMBER
7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 8. PERFORMING ORGANIZATION
Clemson University,Department of Mechanical Engineering,241 REPORT NUMBER
Engineering Innovation Building,Clemson,SC,29634-0921
9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 10. SPONSOR/MONITOR’S ACRONYM(S)
11. SPONSOR/MONITOR’S REPORT
NUMBER(S)
12. DISTRIBUTION/AVAILABILITY STATEMENT
Approved for public release; distribution unlimited
13. SUPPLEMENTARY NOTES
14. ABSTRACT
15. SUBJECT TERMS
16. SECURITY CLASSIFICATION OF: 17. LIMITATION OF 18. NUMBER 19a. NAME OF
ABSTRACT OF PAGES RESPONSIBLE PERSON
a. REPORT b. ABSTRACT c. THIS PAGE Same as 14
unclassified unclassified unclassified Report (SAR)
Standard Form 298 (Rev. 8-98)
Prescribed by ANSI Std Z39-18
Your article is protected by copyright and
all rights are held exclusively by ASM
International. This e-offprint is for personal
use only and shall not be self-archived in
electronic repositories. If you wish to self-
archive your article, please use the accepted
manuscript version for posting on your own
website. You may further deposit the accepted
manuscript version in any repository,
provided it is only made publicly available 12
months after official publication or later and
provided acknowledgement is given to the
original source of publication and a link is
inserted to the published article on Springer's
website. The link must be accompanied by
the following text: "The final publication is
available at link.springer.com”.
1 23
Author's personal copy
JMEPEG (2014)23:1–12 (cid:2)ASM International
DOI: 10.1007/s11665-013-0637-5 1059-9495/$19.00
Concurrent Computational and Dimensional Analyses
of Design of Vehicle Floor-Plates for Landmine-Blast
Survivability
M.Grujicic andB.A. Cheeseman
(Submitted January 3,2013;published online August9,2013)
Development of military vehicles capable of surviving landmine blast is seldom done using full-scale pro-
totype testing because of the associated prohibitively-high cost, the destructive nature of testing, and the
requirementsformajorlarge-scaleexperimental-testfacilitiesandalargecrewofengineerscommittedto
thetask.Instead,testsofsmall-scalemodelsaregenerallyemployedandthemodel-basedresultsarescaled
up to the full-size vehicle. In these scale-up efforts, various dimensional analyses are used whose estab-
lishment and validation requires major experimental testing efforts and different-scale models. In the
present work, an approach is proposed within which concurrent and interactive applications of the com-
putationalanalyses(oflandminedetonationandtheinteractionofdetonationproductsandsoilejectawith
the vehicle hull-floor) and the corresponding dimensional analysis are utilized. It is argued that this
approachcanguidethedesignofmilitary-vehiclehull-floorswhichprovidetherequiredlevelofprotection
tothevehicleoccupantsunderlandmineblastattackwithoutintroducingunnecessarilyhighweighttothe
vehicle.Tovalidatethisapproach,acombinedEulerian/Lagrangianformulationforlandminedetonation
and the interaction of detonation products and soil ejecta with the vehicle hull-floor (developed in our
previous work) has been utilized along with the experimental results pertaining to small-scale model and
full-scale vehicletesting.
investigations of whole-vehicle performance under various in-
Keywords dimensional analysis, landmine detonation, vehicle
service (Ref 1) and threat (Ref 2-4) conditions. However, the
hull-floordesign
true utility of these computational analyses remains question-
able considering that critical shortcomings of these analyses,
are either not fully identified or understood. For example,
material models used in these analyses need substantial
1. Introduction improvements to take into account the contribution of various
microstructural effects (at various length scales), the effects of
componentsmanufacturingandassemblyaswellastheeffects
Recent efforts of the U. S. Army have been aimed at
of realistic environmental, thermal, mechanical, and dynamic
becoming more mobile, deployable, and sustainable while
loading conditions on material behavior (Ref 5). In addition,
maintaining or surpassing the current levels of lethality and
many issues associated with the interactions and contacts
survivability. Current battlefield vehicles have reached in
between the vehicle and its surrounding (including fluid/liquid
excess of 70 tons due to ever increasing lethality of ballistic
interactions) are still not fully resolved. This situation requires
threatswhichhinderstheirabilitytobereadilytransportedand
thatextensiveexperimentaltestingandvalidationeffortsstillbe
sustained. Therefore, a number of research and development
employed during development of new or retro fitting of the
programs are under way to engineer light-weight, highly
existing military vehicles. In other words, computational
mobile, transportable, and lethal battlefield vehicles with a
analyses are presently capable of reducing somewhat the
target weight under 20 tons. To attain these goals, significant
new-vehicle development time and of lowering the extent of
advances are needed in the areas of light-weight materials and
(but not fully eliminating a need for) the experimental testing
structuresaswellasintheareasofcomputer-aideddesignand
and validation efforts.
experimental testing/validation techniques for these vehicles.
The main issue with the experimental testing/validation
Overthelasttwodecades,theuseofcomputationalanalyses
efforts is that they should ideally involve full-scale prototypes
in the design of military vehicles has continuously increased.
ofthevehicles,whichisquiteexpensiveandrequiresmajortest
Efficient computational codes and fast, massive, and inexpen-
facilities and a large crew of engineers committed to the task.
sive computational resources enable today computational
Ontheother hand,itisquite advantageous to carry outall the
experimental testing and validation using small-scale models
M.Grujicic,DepartmentofMechanicalEngineering,ClemsonUniversity,
and to then scale-up the model results to the full-scale
241 Engineering Innovation Building, Clemson, SC 29634-0921;
prototype, using the appropriate dimensional analysis (defined
and B.A. Cheeseman, Army Research Laboratory—Survivability
Materials Branch, Aberdeen Proving Ground, MD 21005-5069. later). In trying to attain this goal, computational analyses can
Contacte-mail:[email protected]. play a critical role. That is, the computational analyses can be
Journal of Materials Engineeringand Performance Volume23(1) January 2014—1
Author's personal copy
firstvalidatedexperimentallyusingasinglesmall-scalemodel. cause a sequence of failures throughout the entire vehicle
Then,thecomputationalanalysescanbeemployedtoestablish through components/sub-system connections and interfaces
the basic fundamental parameters and their inter-relations and, (e.g., via fuel lines running along the floor, floor-bolted seats,
in turn to develop the appropriate dimensional analysis. This ammunition storage racks, power-train lines, etc.). Tradition-
could potentially lead to significant reductions in the small ally, the floor-rupture problem is solved through the use of
scale model testing efforts (e.g., different-scale models would thicker floor-plates, stronger, tougher, and usually heavier
not have to be employed to establish/validate the dimensional material and through the utilization of applique armor. This
analysis). Lastly, both the computational and experimental approach is generally associated with unnecessary oversized/
modelscale-uppredictionscanbeused,withmoreconfidence, overweight vehicle-hull floors and, is being abandoned now-
in the design/development of full-scale prototypes. adaysduetoaforementionedmilitary(cid:2)srequirementsforlighter
A review of the public-domain literature carried out as part vehicles. Combined application of the computational analyses
of the present work revealed that the landmine-detonation ofbothsmall-scalemodelsandfull-scaleprototypesofvehicles
related research activities can be broadly divided into three (or their hulls) subjected to landmine blast and a dimensional
maincategories:(a)shockandblastwavemechanics/dynamics analysis for model-results scale-up, proposed in the present
including landmine-detonation phenomena and large-deforma- work, is a way to reduce the weight of the vehicles while
tion/high-deformationrateconstitutivemodelsfortheattendant ensuring the required confidence level for its survivability and
materials (high explosive, air, soil, etc.) (Ref 6); (b) the survivability of its occupants.
kinematic and structural response of the target to blast loading Theorganizationofthepaperisasfollows:Insection 2.1,a
including the role of target design and use of blast attenuation brief overview is provided of the Buckingham P theorem
materials(Ref7);and(c)vulnerabilityofhumanbeingstopost- (Ref 9). The dimensional analysis for military-vehicle floor-
detonation phenomena such as high blast pressures, spall plates subjected to landmine-detonation attack originally pro-
fragments, and large vertical and lateral accelerations (Ref 2). posed by Wenzel and Hennessey (Ref 8) is overviewed in
The present work falls into the categories (a) and (b) of the section 2.2. The combined Eulerian/Lagrangian computational
research listed above since it emphasizes both the role of procedure used in the present work is described in section 2.
material behavior under transient-dynamic loading conditions Themainresultsobtainedinthepresentworkarepresentedand
as well as the kinematic and structural responses of the target discussedinsection 3,whilethekeyconclusionsresultedfrom
structure to detonation of landmines shallow buried in soil. the present study are summarized in section 4.
The main objective of the present work is to extend the
scaling analysis of Wenzel and Hennessey (Ref 8) for military
vehicle-hullfloor-platessubjectedtoimpulseandhigh-pressure 2. Computational Procedure
loads associated with detonation of a landmine buried in soil
under the hull body. The dimensional analysis of Wenzel and
2.1 Buckingham P Theorem (Ref 9)
Hennessey (Ref 8) is based on the so-called Buckingham
Theorem(Ref9),aprocedurecommonlyusedinthesituations As discussed earlier, the Buckingham P theorem (also
inwhichfull-scalestructures/prototypesaredesignedandsized known as the ‘‘similitude theory’’) (Ref 9) is employed in the
basedontheresultsobtainedusingsmall-scalemodeltests.To presentworktoaddresstheproblemofcomputational-analysis
extend the dimensional analysis of Wenzel and Hennessey validation via small-scale model testing and the problem of
(Ref 8), the combined Eulerian/Lagrangian computational scale-up ofthe model results to the full-scale prototype. In the
model and analysis of detonation of landmines shallow-buried remainder of this section, a brief overview is provided of the
in soil and of the interactions between detonation products similitude theory and of its utility and limitations.
and soil ejecta with the target structure developed by Grujicic The similitude theory (also known as the ‘‘dimensional
et al. (Ref 10) is used. analysis’’) is essentially a technique/procedure for describing
Themainroleofthemilitary-vehiclehull-flooristoprotect thebehaviorofasystemintermsofanumberofdimensionless
the vehicle occupants in the case of detonation of a landmine parameters (generally referred to as the P terms). This
buried in soil underneath the vehicle belly. Complex interre- procedure typically involves the following main steps:
lationships generally exist between the type and the extent of
vehicle-occupants injuries and the extent of vehicle damage (a) Identification of all physical parameters which govern
resulting from under-the-vehicle landmine blast. The most the behavior of the system under investigation. This is
serious vehicle-occupants injuries and their fatalities are most the most critical step in the dimensional analysis since it
frequently incurred as a consequence of hull-floor rupture. In requires a good understanding of the underlying physics
addition to generating floor-material fragmentation, floor rup- of the problem at hand (even in the case when the func-
ture enables for the ingress of flame and toxic (detonation- tional forms of the governing equations are not known).
product) gases and can lead to on-board fires and explosions. (b) Specification of the units for the physical parameters
Whilethevehicleunderlandmineattackisalwayssubjectedto defined in (a) in terms of the fundamental physical
high‘‘G’’)forces,thisaspectofthethreatisusuallyconsidered dimensions (typically force, length, and time).
assecondarywhenfloorrupturetakesplace.Ontheotherhand, (c) Identification of a subset of parameters among the physi-
when landmine-detonation results in only hull-floor bulging, cal parameters defined in (a) with the number of param-
the vehicle is then propelled upward and off the ground, to an eters in this subset being equal to the number of
extent which scales directly with size of the landmine and fundamental physical dimensions defined in (b). The
inversely with vehicle weight. The accompanying dynamic parameters in this subset will be used for converting the
shock can cause vehicle occupants and on-board instruments remaining parameters in (a) into a set of dimensionless
andweaponstobethrownaboutwithinthevehicleinterior.In parameters (i.e., the P terms). It should be noted that
addition, large hull-floor deflections and deformations can while the number of the resulting P terms is fully
2—Volume23(1) January 2014 Journal of MaterialsEngineering andPerformance
Author's personal copy
defined (equal to the number of physical parameters in identified eleven key physical parameters which govern the
(a) minus the number of fundamental dimensions in (b)), responseoftheplatesubjectedtolandmine-detonationloading.
the functional form of the P terms is not unique. It is Theseparametersarelistedandexplainedandtheirfundamen-
generally advantageous to define the P terms in such a tal dimensions are identified in Table 1.
way that they have a clear physical meaning. However,
regardless of the choice of the P terms, the final impli- P Terms. Next, three of the physical parameters listed in
cations/findings offered by the similitude theory is inde- Table 1 (floor-plate thickness, h, floor-plate material yield
pendent of this choice (provided, all the governing strength, r, and the same material density, q) are used for
physical parameters are correctly identified in (a)). constructing eight P terms (Table 2).
(d) Identification of the functional relationships between It should be noted that the rigid body motion of the floor-
each of the P terms and all other P terms. These rela- plate was neglected due to the fact that it did not modify
tions are generally determined using experimental or landmine-detonation-induced loading and, hence, did not
computational approaches or the combination of the two modify the structural response of the floor-plate. The distrib-
approaches. uted inertial effects, on the other hand, were included through
the selection of the material density and the floor-plate
The power of the similitude theory is that the functional thickness.Itshouldbealsonotedthatnoexplicitconsideration
relationsofferedaregeneralized,i.e.,theeffectofgeometrical, oftheeffectofchangeinthestrainrateswasconsideredwitha
kinematic,ambient,loading,andresponsescalesiseliminated. justification that the strain rates do not appreciably change
Consequently, the similitude theory states that, if all the between the full-scale prototype and its replica-model. How-
pertinent physical parameters were correctly identified in (a) ever,itshouldberecognizedthattheeffectofhighstrainrates
and if all the P terms were kept invariant (i.e. unchanged) was included implicitly, via the use of the dynamic material
between a small-scale model and its full-scale prototype, then yield strength, r.
theresultsobtainedusingasmall-scalemodelcouldbedirectly
used to predict the response of the corresponding full-scale
counterpart.
2.2 Brief Overview of the Wenzel and Hennessey Work
(Ref 8)
As mentioned earlier, one of the main objectives of the
present work is to further advance the dimensional analysis
presented by Wenzel and Hennessey (Ref 8). In the remainder
of the section, a brief overview is provided of this work.
The work of Wenzel and Hennessey (Ref 8) deals with
scalingfloor-platesformilitaryvehiclesexpectedtobetargeted
by landmine blast. The main objective of the work was to
develop and validate a dimensional analysis for the vehicle
Fig.1 Aschematic of theprototype landmine-blasttest set-upused
floor-plate sizing problem, so that the results of small-scale
intheworkofWenzelandHennessey(Ref8)
replica-model can be used to determine the minimal required
floor-plate thickness which will ensure that, under the specific
landmine-detonation scenario, the vehicle floor will not under-
goruptureorextensivedynamicdeflection.Theterm‘‘replica-
Table 1 Key physical parameters controlling structural
model’’ will be defined in greater detail later in this section.
response of the vehicle-floor plates
Untilthen,thistermwillbeusedtodenoteasmall-scalemodel
whichissimilar(geometrically,constitutively,andwithrespect Fundamental
to extent of properly scaled loads and structural responses) to Symbol Parameter dimensions
the corresponding full-scale prototype.
The starting point in the work of Wenzel and Hennessey P Pressureduetolandmine F/L2
(Ref 8) was to take advantage of the fact that floor-plate detonation
deflections do not significantly modify the (landmine-detona- I Specificimpulseduetolandmine FT/L2
detonation
tion-induced)dynamicloading.Towardthatend,theproblemat
h Floor-platethickness L
hand was divided (decoupled) in two separate problems: (a) a
h Arrayofother(non-dimensional) …
floor-plate structural response problem and (b) the problem of i
floor-platedimensions
analyzingpressuresandimpulsesresultingfromdetonationofa r Floor-platematerialyieldstrength F/L2
landmine (shallow buried in soil underneath the vehicle-hull r Arrayofother(non-dimensional) …
i
bottom). strength parameters for floor-plate
2.2.1 Structural Response. Physical Parameters. By material
carrying out a careful and detailed engineering analysis of the e Floor-platematerialequivalent …
structural response of a thick metal plate (a model for the strain
vehicle-hullbottom)supportedonwoodenblocks(toobtainthe q Floor-platematerialdensity FT2/L4
d Floor-platedeflection L
required level of ground clearance), and clamped along its
a Floor-plateacceleration L/T2
edges (to simulate the effect of the surrounding/supported
t Postlandmine-detonationtime T
vehicle-framestructure),Fig. 1,WenzelandHennessey(Ref8)
Journal of Materials Engineering andPerformance Volume23(1) January 2014—3
Author's personal copy
Table 2 P terms controlling structural response of the Table 4 Parametersfordefiningloadingfromalandmine
vehicle-floor plates
Fundamental
Similarity Symbol Parameter dimensions
Pterm Physicalmeaning type
P Pressureduetolandminedetona- F/L2
(cid:3) =h Arrayofnon-thickness Geometrical tion
1 i
(non-dimensional) I Specificimpulseduetolandmine FT/L2
floor-platedimensions detonation
(cid:3) =r Arrayofnon-yield-strength Constitutive E Landminedetonation-releaseden- FL
2 i
(non-dimensional) ergy
parametersforfloor-platematerial r Circular-diskshapelandmine …
(cid:3) =e Floor-platematerialequivalent Constitutive, thickness
3
strain response r Arrayofnon-thickness(non- …
i
(cid:3) =d/h Thickness-scaledfloor-plate Response dimensional)landminedimensions
4
deflection d Landminedepth-of-burial L
(cid:3) =aqh/r Non-dimensionalfloor-plate Response R Landmine/floor-platestandoffdis- L
5
acceleration tance
(cid:3) =(tr1/2)/ Non-dimensionaltime Response P Ambientatmospheric-pressure F/L2
6 0
(hq1/2) a Speedofsoundinair L/T
0
(cid:3) =I/ Non-dimensionalimpulse Loading c Ratioofspecificheatsforair …
7 0
(r1/2q1/2h) duetolandmine q Densityofsoil FT2/L4
detonation c Seismicvelocityofsoil L/T
(cid:3) =P/r Non-dimensionalpressure Loading g Gravitationalacceleration L/T2
8
duetolandmine
detonation
Table 3 Replica-to-prototypescalingforstructuralresponse
2.2.2 Landmine-Detonation-InducedImpulsesandPres-
ofthevehicle-floorplates
sures. Wenzel and Hennessey (Ref 8) also carried a detailed
engineering analysis of detonation of landmine shallow-buried
Symbol Parameter Scalefactor
in soil. In order to quantify the resulting time-dependant
h Arrayofnon-thickness k impulses and pressures (input loading parameters in the
i
(non-dimensional)floor-plate structural response analysis), they included the effect of the
dimensions landmine shape and size, type of explosive, soil properties as
q Floor-platematerialdensity 1.0 well as the properties of the surrounding (ambient) air. These
r Floor-platematerialyieldstrength 1.0 parameters along with their fundamental dimensions are listed
e Floor-platematerialequivalent 1.0 in Table 4.
strain
d Floor-platedeflection k
Landmine Related Parameters. The landmine was considered
a Floor-plateacceleration 1/k
to be of a circular-disk shape and to have a size defined by
t Postlandmine-detonationtime k
I Specificimpulseduetolandmine k thickness, r, and an array of non-dimensional parameters, ri
detonation (i = 1, in the present case and r1 represents thickness-formal-
P Pressureduetolandminedetona- 1.0 izedlandminediameter).Thetypeofexplosiveinthelandmine
tion (of a givensize) isthen expressed by thetotal energyreleased
during landmine detonation, E. The location of the landmine
Scaling. As mentioned earlier, in accordance with the simili- with respect to the model/prototype is defined in terms of the
tudetheory,fortheresultsobtainedusingsmall-scalemodelto depth-of-burial, d (a distance between the top surface of the
be applicable at the prototype scale, the P terms must be landmine and the soil/air interface) and standoff distance,
unchanged between the model and the prototype. In order to R (distance between the top surface of the landmine and the
ensure this condition, the model and the prototype must be bottom of the test plate).
similar geometrically, (i.e., P is unchanged), constitutively
1
(i.e.,P andP termsareunchanged),inresponse(i.e.,P -P Ambient-AirParameters. Toaccountfortheeffectsofairblast,
2 3 3 6
terms are unchanged) and with respect to loading (P and P three ambient-air parameters are considered: ambient pressure,
7 8
terms are unchanged). Hence, if the model is obtained by P , sound speed in air, a , and a ratio of the constant-pressure
o 0
scalingdowntheprototypebyafactork(whilemaintainingall and constant-volume specific heats in air, c .
o
theanglesunchanged),themodelandtheprototypearemadeof
the same material (i.e., r and q are unchanged), specific Soil Parameters. To account for the effect of soil (into which
i
impulse due to landmine detonation, I, is scaled by a factor k the landmine is buried) Wenzel and Hennessey (Ref 8),
while pressure due to landmine detonation, P, is kept un- assumed that detonation of shallow-buried landmine can be
changed, then e, d, a, and t responses of the model are scaled relatedtotheproblemofundergroundlandmineexplosionand
withrespecttotheirprototypecounterpartsinaccordancewith theresultingcraterformation.Consequently,theyassumedthat
the values given in Column 3 of Table 3. When the model is two soil-specific parameters should be included: the (initial)
scaled down in accordance with these conditions, then the soildensity,q,anditsseismicvelocity,c.Inaddition,toinclude
model is referred to as the replica-model. the effect of gravitational forces of the soil surrounding the
4—Volume23(1) January 2014 Journal of MaterialsEngineering andPerformance
Author's personal copy
Table 5 P terms for determining loading from landmine betweenthesmall-scalemodelandthefull-scaleprototype.The
question then arises as to how to scale the landmine/air/soil-
Similarity related physical parameters listed in Table 4 to achieve this
Pterm Physicalmeaning type condition. In addition, the scaling for the pressure and for the
impulse has already been fixed by the structural-response
(cid:3) =r Arrayofnon-thickness Geometrical
1 i analysis,Table 3andthisscalingcannotbechangedwithinthe
(non-dimensional)
present landmine-detonation-analysis.
landminedimensions
If r, d and, R are scaled by k then P and P terms are
(cid:3) =d/R Scaleddepth-of-burial Geometrical 1 3
2 invariant between the model and the prototype. If the same
(cid:3) =r/R Inverseofscaledstandoffdistance Geometrical
3
(cid:3) =c Ratioofspecificheatsforair Atmospheric ambient air and soil is used in both the model and the
4 0
(cid:3)5=a0/c Airsoundspeedtosoilseismic Kinematic prototype,thenP4andP5areinvariant.Also,theinvarianceof
speedratio the P and P terms requires for E to be scaled as k3. This is
6 7
(cid:3) =P R3/E Sach(cid:2)sNumber Ambient reasonable since all the spatial dimensions are scaled in k and
6 0
(cid:3)7=qc2R3/E ModifiedSach(cid:2)sNumber Soil E depends on the landmine volume.
(cid:3)(cid:3)8==gPRR/3c/2E NFroonu-ddeimNeunmsiboenralpressuredueto LGoraavdiitnagtional Since E scales with k3 and R scales with k, then from P9
9 and P it is seen that P remains fixed while I is scaled as k
landmine 10
betweenthemodelandtheprototype.Asmentionedearlier,this
detonation
is exactly the scaling for specific impulse and pressure which
(cid:3) =IcR2/E Non-dimensionalimpulsedueto Loading
10 was required by the structural-response analysis.
landmine
detonation The invariance of the P8 term, requires that gravitational
accelerationbescaledas1/k.Inotherwords,forexample,fora
one-quarter model, model experiments would have to be done
under conditions corresponding to a four times higher gravi-
Table 6 Replica-to-prototypescalingforblast-loadrelated
tationalaccelerationthantheonenormallypresentattheearth.
parameters
Wenzel and Hennessey (Ref 8) recognized that conducting
Scale small-scale experiments under hyper-gravitational conditions
Symbol Parameter factor would be quite challenging. Consequently, they invoked the
studyofWestine(Ref11)whichshowedthatfailingtoaccount
r,R,d Circular-diskshapelandmine k for the gravitational effects during underground landmine
thickness, explosionformodelinarangebetween1/2and1/8,resultedin
landmine/floor-platestandoffdistance,
relatively small changes in the resulting crater size. Hence,
landminedepth-of-burial
Wenzel and Hennessey (Ref 8) argued that even though the
c,a Seismicvelocityofsoil,speedof 1.0
0 invarianceoftheFroudeNumber(P termwhichisdefinedas
soundinair 8
a ratio of the gravitational and inertial effects) cannot be
c Ratioofspecificheatsforair 1.0
0
q Densityofsoil 1.0 generallysatisfied,thesmall-scalemodelresultscouldbeused,
E Landmine-detonation-released k3 withsignificant confidence,forfull-scale prototypedesignand
energy sizing.
P0,P Ambientatmospheric-pressure, 1.0 2.2.3 Validation of the Dimensional Analysis. To vali-
pressureduetolandminedetonation date their dimensional analysis presented above, Wenzel and
I Specificimpulseduetolandmine k Hennessey (Ref 8) carried out a series of one-half and one-
detonation
quartermodelexperimentsoffloor-plateresponsetolandmine-
detonation loads, and compared the results with available full-
scaleprototypeexperimentaldata.Theyfoundthattheone-half
andone-quartermodelresults(presentedindetailsinsection 3)
landmine, gravitational acceleration, g, is also defined as a
are fully mutually consistent which constituted the proof for
physical parameter affecting the temporal evolution of land-
validity of the similitude theory they developed. On the other
mine-detonation-induced impulses and pressure.
hand, the full-scale results could not be used for dimensional-
analyses validation, due to the lack of knowledge of the
Resulting Landmine-Detonation Loading Parameters. The attendant soil parameters. Nevertheless, the full-scale results
discussion presented above introduced 11 landmine, air and revealed that soil parameter play important role in the
soil-related parameters. The remaining two parameters are landmine-detonation-induced loading.
associated with the landmine-detonation-induced specific
impulse, I, and the associated (peak) pressure, P.
2.3 Computational Problem and Solution Method
P Terms. Three of the parameters listed in Table 4 (i.e., soil Thebasicformulationofthecomputationalproblemdealing
density, q, soil seismic velocity, c, and the total detonation- with interactions between the detonation products, shell frag-
released energy, E) are used for non-dimensionalization of the ments, and soil ejecta (all resulting from the explosion of a
remaining parameters. The resulting 10 P terms are listed in shallow-buried landmine) and the target model/prototype hull-
Table 5. bottomplateispresentedinthissection.Allthecalculationswere
carriedoutusingANSYS/Autodyn(Ref12),ageneralpurpose
Scaling. Asmentionedearlier,inorderfortheresultsobtained transientnonlineardynamicsprogram.Thecomputationalmod-
usingasmall-scalemodeltobeapplicabletothecorresponding elingofthisinteractioninvolvedtwodistinctsteps:(a)geomet-
full-scaleprototype,allPtermsmustbemaintainedunchanged ricalmodelingofthemodel/prototypealongwiththeadjoining
Journal of Materials Engineering andPerformance Volume23(1) January 2014—5
Author's personal copy
mine, air, and soil regions, and (b) the associated transient
nonlineardynamicsanalysisoftheimpulseloading(momentum
transfer)fromthedetonationproducts,shellfragments,andsoil
ejectatothetarget.Thepart(b)ofthisanalysiswasperformed
usingamodifiedversionofthetechniquedevelopedbyFairlie
andBergeron(Ref13)whichcouplesamulti-materialEulerian
meshtotwoLagrangianmeshes.TheEulerianmeshcontained
initiallyaComposition-Bsolid-circular-diskshapeburiedland-
mine (and, after explosion, the resulting high-pressure, high-
internal energy-density gaseous detonation products) and the
(initiallystationary,atmospheric-pressure)air(locatedabovethe
soil).TheEulerianmeshwasconstructedintermsofeightnode
cuboidalcells.OneoftheLagrangianmesheswasusedtomodel
thesoil,whiletheothertorepresentthefloor-plate(alongwiththe
adjoininghull-supportedframestructure).Thesoilandthetarget
werebothdiscretizedusingeight-nodesolidelements(Table 6).
An advantage was taken of the inherent symmetry of the
model.Inotherwords,twomutually-orthogonalverticalplanes
of symmetry were placed along the axis of the target plate as
well as along the axis of the air, landmine, and soil regions
whichenabledonlyaquarterofthecomputationalmodeltobe
analyzed.Representativequartersymmetricmodelsforvarious
computational domainsused inthe present study are shownin
Fig. 2.ItshouldbenotedthatthelowerportionoftheEulerian
domain contains the landmine while the rest of the lower
portion of the Eulerian domain is occupied by the Lagrangian
soil domain. Likewise, the upper portion of the Eulerian
domainwhichextendsabovethesoilcontainsinitiallyairandis
partially occupied by the Lagrangian floor-plate domain.
Air and Composition-B detonation-products are treated as
fluidsandmodeledusingonlytheappropriateequationofstate,
EOS [i.e., the gamma EOS (Ref 12) for air and the JWL EOS
(Ref 12) for Composition-B]. The floor-plate was treated as
being made of Rolled Homogenized Armor (RHA) and
modeled using a linear equation of state, the Johnson Cook
strength and failure models and an instantaneous total equiv-
alenterosionstrainalgorithm.Detailsofthesemodels/relations
and their parameterizations can be found in Ref 10. Soil was
represented using our CU-ARL sand material model
Fig.2 Computational sub-domains used in the present analysis of
(Ref 14-20). The Lagrangian soil elements were allowed to
structural response of vehicle-hull floor-plate subjected to landmine-
erode at an incremental geometric strain of 2.0. Upon the
blast dynamic loading: (a) Euler sub-domain and (b) Lagrange sub-
erosion of an element, the resulting free nodes are allowed to domain
retain their mass/velocity (and, thus, momentum) and to
continue to interact with the soil and the target plate. Several gage points were defined within the landmine, soil,
Interactions between the gaseous fluids (i.e., air and air, and floor-plate which allowed monitoring of the quantities
detonation products) with the target plate were accounted for such as pressure, velocity, and (in the case of the Lagrange
through the use of the appropriate Euler/Lagrange coupling domains) the vertical displacements.
option within ANSYS/Autodyn (Ref 12). Likewise, the soil/ At the beginning of the simulation, all the Lagrange and
target-plate interactions were modeled through the use of the Eulerdomainswereactivatedandthelandminedetonated.The
appropriate Lagrange/Lagrange coupling option. (circular-disk shape) landmine was detonated over its entire
The‘‘flowout’’boundaryconditionswereappliedtoallthe bottom face at the beginning of the simulation. A standard
freefaces(thefaceswhichdonotrepresentinterfacesbetween meshsensitivityanalysiswascarriedout(theresultsnotshown
the different domains) of the Euler domain except for the face for brevity) to ensure that the results obtained are not very
associated with the vertical symmetry planes. To reduce the sensitive to the size of the cells/elements used.
effect of reflection of the shock waves at the outer surfaces of
the Lagrange soil domain, ‘‘transmit’’ boundary conditions
were appliedtoallthefreefaces ofthisdomain except forthe
3. Results and Discussion
faces associated with the vertical symmetry planes and the
upper face which defines the soil/air interface. The transmit
boundary conditions enable propagation of the pressure waves To comply with the experimental procedure of Wenzel and
across the boundaries without reflection mimicking wave Hennessey(Ref8),allthecomputationalanalysescarriedoutin
propagation in an infinitely-large soil domain (Ref 12). the present work involved a single (full-scale) size of the
6—Volume23(1) January 2014 Journal of MaterialsEngineering andPerformance
Author's personal copy
vehicle-hull floor-plate with L9W9H dimensions of 2.13 m (Ref 8) is provided in the subsequent sections. It should be
by 1.52 m by 0.0254 m. Following the present dimensional recalledthattwosetsofPtermswereintroducedinthepresent
analysis,thesedimensionswerereducedbyafactorof2and4 dimensional analysis. One set associated with structural
in the case of one-half and one-quarter models, respectively. response of the floor-plate while the other corresponding to
The initial soil density and its seismic velocity (sound speed) various aspects of detonation of a landmine buried in soil. To
were selected to match their experimental counterparts in the enable clear differentiation between the two sets of P terms
work of Wenzel and Hennessey (Ref 8). Also, the same superscriptsSRandBLareaddedtotherespectivePtermsto
explosive-type(CompositionB)andfloor-platematerial(Rolled denote their structural-response/blast-loading nature.
Homogenized Armor, RHA) were used to match those used in
Ref 8.Initially,twovaluesofthestandoff-distance normalized
3.1 Floor-Plate Deflection Versus Explosive-Charge Energy
ground clearance (G.C./R) of 0.8889 and 0.8235 (correspond-
ing to ground-clearance values of 14¢¢ = 0.3556 m and In this section, the experimental results from Ref 8 and the
24¢¢ = 0.6096 m, respectively and a fixed/single value of the present computational results pertaining to the dependence of
depth-of-burialof3¢¢ = 0.0762 m)wereselected,inagreement thickness-scaled floor-plate average displacement (PSR, in the
4
withthoseusedbyWenzelandHennessey(Ref8).Inorderto structural-responseanalysis)onthescaledchargeenergy(PBL,
7
obtain further insight into the dimensional analysis of land- of the detonation-induced impulse/pressure analysis) are pre-
mine-blastinduceddeformationofthevehicle-hull-floor-plates, sented and discussed.
additional ground-clearance values were also used in the 3.1.1 ExaminationofthePBLTermasDefinedbyWenzel
7
present computational analyses. The explosive-charge sizes and Hennessey (Ref 8). In Fig. 4(a), the PSR versus PBL
4 7
used (at the full scale) at different ground-clearance values are results from the work of Wenzel and Hennessey (Ref 8) are
listed in Table 7. presented for the two aforementioned levels of normalized
Anexampleofthepressurefieldinairduringdetonationof ground clearance. In addition to the measured data points, the
a landmine buried in soil underneath a vehicle-hull floor-plate corresponding power-law best-fit curves are also shown. The
is shown in Fig. 3. A more detail representation and the corresponding computational counterpart results obtained in the
discussionoftheresultsobtainedinthepresentworkandtheir presentworkaredisplayedinFig. 4(b).Furthermore,toassistan
experimental counterparts obtained by Wenzel and Hennessey experiment-with-calculationcomparison,thebest-fitcurvesfrom
Fig. 4(a) are reproduced in Fig.4(b). A simple examination/
comparison of the results displayed in Fig.4(a) and (b)
Table 7 Testmatrixpertainingtotheselectionofexplosive-
reveals that
chargesizes,depth-of-burialandfloor-platedimensionsat
differentfull-scaleground-clearancelevels
(a) At each level of the ground clearance, the corresponding
Full-scalegroundclearance,m full-scale, half-scale, and quarter-scale model computa-
tional results are practically identical. This finding sug-
Parameter 0.356 0.483 0.609 gests that the choice of the two P terms in Fig. 4(b),
was appropriate/justified. Similar conclusion can be
Composition-Bexplosive 1.81 4.54
drawn based on the experimental results of Wenzel and
charge,kg 3.18 9.07
Hennessey (Ref 8), Fig. 4(a). Nevertheless, it should be
4.54 18.14
noted that all the materials (Composition-B, soil, air and
9.07 …
RHA) involved in the present dimensional analysis were
13.61 …
Depth-of-burial,m 0.076 kept unchanged. Also, at a given level of scaled ground
Standoffdistance,m 0.432 0.559 0.686 clearance, P7BL value in Fig. 4(a) and (b) was changed
Floor-plate,L9W9H 2.13491.52490.025 by merely changing the value of explosive-charge
dimensions,m energy, E (via changing the charge thickness, r). Conse-
quently, it is not clear at this point if the PSR versus
4
PBL relations suggested in Fig. 4(a) and (b) will be
7
valid if properties of the attendant materials and explo-
sive-charge shape are changed. This point will be
addressedingreater detail in the remainder of this section.
(b) At each scaled ground-clearance level, the corresponding
best-fit curves for the experimental and computational
data differ by less than 10%. This finding suggests that
there is a reasonably good agreement between the exper-
imental test transient blast-loading conditions used and
the floor-plate response observed in the work of Wenzel
and Hennessey (Ref 8) and the present transient nonlin-
ear dynamic computational modeling and analysis proce-
dure of the same problem. This finding, and the fact
that ground-clearance level is seen to have a significant
effect on the PSR versus PBL relationship, suggested
4 7
the need for additional computational analyses at other
values of ground clearance. Since the two aforemen-
Fig.3 An example of the air-pressure field during detonation of a tioned ground-clearance values of 14 and 24 inches
mineshallow-buriedinsoilunderafloor-platetarget define a fairly realistic range of this quantity in the case
Journal of Materials Engineering andPerformance Volume23(1) January 2014—7