Summary

Blast Quantification Using Hopkinson Pressure Bars

Published: July 05, 2016
doi:

Summary

This protocol details the use of Hopkinson pressure bars to measure reflected blast loading from near-field explosive events. It is capable of interpolating a pressure-time history at any point on a reflective boundary and as such can be used to fully characterize the spatial and temporal variations in loading produced.

Abstract

Near-field blast load measurement presents an issue to many sensor types as they must endure very aggressive environments and be able to measure pressures up to many hundreds of megapascals. In this respect the simplicity of the Hopkinson pressure bar has a major advantage in that while the measurement end of the Hopkinson bar can endure and be exposed to harsh conditions, the strain gauge mounted to the bar can be affixed some distance away. This allows protective housings to be utilized which protect the strain gauge but do not interfere with the measurement acquisition. The use of an array of pressure bars allows the pressure-time histories at discrete known points to be measured. This article also describes the interpolation routine used to derive pressure-time histories at un-instrumented locations on the plane of interest. Currently the technique has been used to measure loading from high explosives in free air and buried shallowly in various soils.

Introduction

Characterizing the output of explosive charges has many benefits, both military (defending against buried improvised explosive devices in current conflict zones) and civilian (designing structural components). In recent times this topic has received considerable attention. Much of the knowledge gathered has aimed at the quantification of the output from charges to enable the design of more effective protective structures. The main issue here is that if the measurements made are not of high fidelity then the mechanisms of load transfer in these explosive events remain unclear. This in turn leads to problems validating numerical models which rely on these measurements for validation.

The term near-field is used to describe blasts with scaled distances, Z, less than ~1 m/kg1/3, where Z = R/W1/3, R is the distance from the center of the explosive, and W is the charge mass expressed as an equivalent mass of TNT. In this regime the loading is typically characterized by extremely high magnitude, highly spatial and temporally non-uniform loads. Robust instrumentation is hence required to measure the extreme pressures associated with near-field loading. At scaled distances Z < 0.4 m/kg1/3, direct measurements of the blast parameters are either non-existent or very few1 and the semi-empirical predictive data for this range is based almost entirely on parametric studies. This involves using the semi-empirical predictions given by Kingery and Bulmash2, which is outside of the author's intended scope. Whilst tools based on these predictions3,4 allow for excellent first-order estimations of loading they do not fully capture the mechanics of near-field events, which are the focus of the current research.

Near-field blast measurements have in recent times focused on quantifying the output from buried charges. The methodologies employed vary from assessing the deformation caused to a structural target5-7 to direct global impulse measurement8-13. These methods provide valuable information for the validation of protective system designs but are not capable of fully investigating the mechanics of load transfer. Testing can be done at both laboratory scales (1/10 full scale), or at near to full scale (> 1/4), with pragmatic reasons such as controlling burial depth or ensuring no inherent shape of the shock front is generated by the use of detonators rather than bare charges14. With buried charges the soil conditions need to be highly controlled to guarantee the repeatability of the testing15.

Independent of the whether the charge is placed in free air or is buried, the most fundamental issue in measuring the resulting blast is ensuring the validity of measurements made by the instrumentation deployed. In the designed test apparatus16 a fixed 'rigid' target plate is used to shield the Hopkinson pressure bars17 (HPBs) whilst at the same time ensuring that the ends of the bars can only record the fully reflected pressures. The authors have previously shown that measurement of reflected pressure from a rigid target is more accurate and repeatable than incident, or 'free-field' measurements18-20. The geometry of this plate is such that any pressure relief generated by clearing or flow around the target edge21 would be negligible. This new test apparatus has been constructed at 1/4 scale. At this scale tight control over the burial conditions and the explosives can be ensured, with the full scale charge size of 5 kg scaled down to 78 g, at a burial depth of 25 mm.

Protocol

1. Rigid Reaction Frame

  1. Determine scaled distance at which testing will take place using Equation 1, where R is the distance from the center of the explosive, and W is the charge mass expressed as an equivalent mass of TNT.
    Z = R/W1/3     (1)
  2. Calculate approximate maximum impulse this arrangement will generate via numerical modelling (see Appendix A) or specific tools such as ConWep3.
    Note: The use of ConWep3 is only valid for free air blast, if an estimation of the pressures generated from buried charges is required the more advanced numerical modelling is required.
  3. Check the estimated loading from the modelling will not generate in-plane displacements of more than 0.5 mm in the target plate.
  4. Increase the loading calculated by a factor of 10 to account for inaccuracies in the modelling and to add flexibility for future testing.
  5. Design a rigid reaction frame to be able to resist the maximum loading calculated16. In an Engineering department, perform these calculations in house; otherwise seek the services of a Structural Engineer.
    1. Procure rigid reaction frames, contract a specialist contractor to fabricate and install the frames to the designs of the structural engineer.
  6. Procure target plate, contract a specialist steel fabricator.
    Note that the plate will need to be mounted on load cells (if used) and that holes for the HPBs (designed in section 3) will need to be drilled through the plate before mounting.

Figure 1
Figure 1. Schematic of the test frame. (A) Overall arrangement, (B) plan of target plate, (C) close-up view of target plate. The Hopkinson pressure bars are hung from the bar assembly receiver so that they sit flush with the face of the target plate. This allows the fully reflected pressure acting on the target plate to be recorded. Please click here to view a larger version of this figure.

2. Load Cell Design

  1. Procure or fabricate load cells (if used). These can either be off-the-shelf universal (compression/tension) strain-gauge canister models or built in-house using sections of thick wall mild steel tubing welded to mounting plates with strain gauges affixed in a Wheatstone bridge formation as shown in Figure 2.
  2. If the load cells have been fabricated in-house, send them to an external contractor for calibration.

Figure 2
Figure 2. Diagram of the in-house fabricated load cells. (A) Side elevation, (B) end elevation. The dark grey cylinder is a thick wall steel tube which strains under loading. This strain is recorded using a single strain gauge as no rotation is experienced during the loading. From the calibration of the load cell the strain can be related back to the stress applied. Please click here to view a larger version of this figure.

3. Hopkinson Pressure Bar Design

  1. Determine the duration of recording, Equation 9, required to capture the full loading from the blast. The minimum duration required is the time taken in the numerical model (section 1.2) for the pressure to return to zero, after the initial pressure spike. Here, use 1.2 msec.
  2. Decide on the material of choice for the HPBs. This affects the elastic wave speed, Equation 10, in the bar which is given by Equation 11 where Equation 12 is the Young's modulus and Equation 13 is the density. For measuring a high pressure shock, use stiff materials such as steel; where as if a weaker shock is expected, use less stiff materials such as a magnesium alloy or even nylon.
  3. Choose the position on the HPB that the strain gauge will be positioned, being as close as possible to the loaded face of the HPB to minimize dispersion. In the current set-up the thickness of the target plate and the manoeuvrability required to fit the bars in place meant that the gauges could only be installed 250 mm from the loaded face.
  4. Calculate the HPB length required using Equation 14, where Equation 15 is the distance from the loaded face of the HPB to the strain gauge and Equation 16 (3.25 m).
  5. Determine required HPB radius to have sufficient bandwidth to capture the event using: Equation 17 kHz, where Equation 18 is the HPB radius in mm22,23 (5 mm).
  6. Decide on the spatial resolution required to capture the distribution of pressure across the plate. This is generally as close as possible while maintaining the structural integrity of the target plate. In the current work, use 25 mm.
  7. Drill holes in the target plate to mount the HPBs (this can be part of the fabrication process). A close fit is required without the HPBs being in contact with the plate. Here, use 0.5 mm tolerance with 17 holes being drilled in a cross shape (Figure 1b).
  8. Procure the HPBs (17), making sure to have the distal ends threaded to allow for suspension in the bar assembly receiver (Figure 3A).

4. Experimental Setup & Data Acquisition

Note: With the reaction frame, target plate, load cells and HPBs designed and fabricated, assembly can begin as shown in Figure 1, and designed in protocol section 1.

  1. Attach semiconductor strain gauges to HPBs (Figure 3B) and load cells using cyanoacrylate, being careful to ensure continuity of earth through all cabling. An example of the Wheatstone bridge used for the HPBs is shown in Figure 3C.
    1. Verify all earth cables are attached to ensure continuity of earth. Well earthed test apparatus will improve signal quality notably.
  2. Ensure wiring is sufficiently long to make sure the oscilloscope is locatable in a blast free area (shielded wiring should be used which has sufficient signal bandwidth).
  3. Fit the target plate to the rigid reaction frame, using the optional load cells if present (Figure 1C).
  4. Hang HBPs from the bar assembly receiver, passing the loaded end through the correct hole in the target plate. Hang the HPBs freely from a nut screwed onto the threaded distal end of the HPB.
  5. Ensure bars are vertical using a spirit level (adjusting the receiver accordingly).
  6. Check the faces of the HPBs are level with the target plate, adjusting the nut accordingly.
  7. Set the trim on the variable resistor in the conditioning circuit (Figure 3C) to keep voltage within the limits of oscilloscope during testing. Do this through trial and error aiming to set the out of balance for each channel as seen on the digital readout on the amplifier boxes to zero.
  8. Connect the amplified gauge output to a suitable digital oscilloscope. Configure to have a sampling frequency (1.56 MHz), recording duration (28.7 msec) with a pre-trigger duration of 3.3 msec.
    1. Set the recording to trigger when the voltage in the break wire channel (which is itself wired into the oscilloscope) exceeds an 'out-window'. Record voltage for each gauge connected (22 in total, 17 HPBs, 4 load cells and the break wire) and time.

Figure 3
Figure 3. (A) Diagram of a HPB fitted into the target plate, (B) section through HPB at gauge location, (C) example Wheatstone bridge circuit. Two strain gauges are used in the Wheatstone bridge so that and bending of the Hopkinson bar is cancelled out. Please click here to view a larger version of this figure.

5. Explosive preparation

  1. Decide on the explosive charge mass and stand-off to be used in the tests (100 g PE4 at 75 mm).
  2. Decide whether the charges are to be detonated in free air or within another medium (soil, water etc.). For free air tests a spherical charge shape is normally utilized whereas with buried charges the standard is a 3:1 squat cylinder24,25.
  3. For free air tests:
    1. Suspend the charge below the target plate at the correct stand-off (75 mm). Achieve this with a thin timber strip or by placing the charge on a sheet of polythene.
    2. Place the charge co-axially with the measurement array to ensure valid readings.
    3. For free air tests use an electrical detonator, with the detonator being placed half way into the charge from the base. Do this at the last moment before firing and when the range has already been made safe.
  4. For buried tests:
    1. Fabricate a suitable container for the medium. For soils, the current testing uses 1/4 scale containers23.
    2. Decide upon the soil type to be utilized and the geotechnical conditions: moisture content and dry density of the soil, see ref.15 for more details.
    3. Decide on the burial depth to use in the testing. This is usually 100 mm in a full scale test, as the current tests are done at ¼ scale this means a 25 mm burial depth.
    4. Mix the soil thoroughly using a suitably sized construction mixer to achieve the target moisture content. For sands the mixing time required is 10 min.
      1. Check the moisture content of the mix by removing a small amount and weigh it to calculate the total mass, Equation 19. Dry the removed soil and re-weigh to calculate the mass of water, Equation 20. Geotechnical moisture contents are specified in terms of gravimetric moisture content, Equation 21.
      2. If the moisture content is within tolerance continue, otherwise remix the soil. A tolerance of ±0.05-0.1% has been achieved in the current work.
    5. Weigh the empty soil container and calculate the volume to enable calculation of the soil density once full (step 5.4.7).
    6. Compact the soil in layers, thin enough to guarantee the target density, ensuring that the mass of soil entering the container is known. For Leighton Buzzard Sand15 this is done in two layers.
    7. Once the container is full, check that the density of the soil within is in tolerance (±0.2%). The target dry density in all tests with Leighton Buzzard Sand was 1.6 Mg/m3. Calculate dry density, using Equation 22, where ρd is the dry density, M is the total mass of soil added to the container, V is the volume of the soil container and w is the moisture content.
    8. Excavate a small hole ≈50 mm to allow the charge to be placed with the top surface at the correct burial depth (25 mm).
    9. Place a non-electrical detonator into the base of the charge, and excavate a suitable channel to the side of the container to ensure the top surface of the container is uninterrupted once the soil is replaced.
    10. Place charge and detonator into the excavated hole, checking the burial depth is correct. Back fill the hole with the excavated material.

6. Firing sequence

Note: there is a small amount of overlap with protocol section 5 due to the nature of the testing. The firing sequence should aim to minimize risk and should only be conducted by suitably trained staff.

  1. For free air tests:
    1. Arrange charge support below the target plate at the correct stand-off (75 mm).
    2. Close the range. Deploy sentries to ensure range is clear during firing.
    3. Place charge on the support co-axial to the instrumentation. Attach the break wire to the detonator, and place the detonator in the charge.
  2. For buried tests:
    1. Place soil container so that the charge is placed co-axial to HPB array.
    2. Close the range. Deploy sentries to ensure range is clear during firing.
    3. Connect the break wire, ensuring it is wrapped around the periphery of charge (this gives a more repeatable time of detonation in buried charges).
  3. Move to firing point and confirm instrumentation is running.
  4. Supply power to the break wire. Check with sentries it is safe to proceed with the firing.
  5. Initiate explosives. Make the test area safe.
  6. Download and back up data.
  7. Re-open test range.

7. Numerical interpolation for a 1D HPB array

  1. Import the data from the raw data files into Matlab.
  2. Time-shift all data in the radial direction so that the peak pressure for each bar arrives at the same time as the peak pressure of the central bar using Equation 2 (Figure 4B).
    Equation 23     (2)
  3. Interpolate the pressure at any radial distance from Figure 4B.
  4. Plot the arrival times (Equation 24) used to align the peak pressures and fit a cubic equation through the data (Figure 4C).
  5. Time-shift the interpolated data to fit the arrival times, generating a continuous shock front (Figure 4D).
  6. Repeat for each individual set of test data.

Figure 4
Figure 4. Interpolation sequence for 1D HPB array. (A) Original data, (B) time-shifted data, (C) shock front arrival times, and (D) final interpolated pressure time data16. The discrete nature of the pressure time histories can clearly be seen in (A) with there being no continuity between the peak pressures at each of the five gauge locations. When aligned by peak pressure as in (B) the interpolation of pressure at any radial distance (assuming the same arrival time) is possible. By recording the time shift required to align the peak pressures the arrival time of the shock front can be calculated as shown in (C). This then allows the arrival time and pressure time history to be calculated for any radial distance be interpolation of pressure from (B) and time from (C) giving the final interpolated pressure as seen in (D). Please click here to view a larger version of this figure.

8. Numerical interpolation for a 2D HPB array

Note: The code used to run the interpolation in Matlab has been provided along with an example results file which will be referred to in this section.

  1. Import the data from the raw data files into Matlab. For the example test data, double click on the test_data.mat file, and then click 'Finish' in the Import Wizard.
  2. Open the interpolation2d.m Matlab script.
  3. Define a regular grid over which the interpolation will run by changing the mesh. Ensure this is the same resolution as the mesh in any future numerical modelling26,27. This is set in the '%mesh details' section of the code.
  4. Run the interpolation2d.m Matlab script. Note the following steps are implemented in the code and are listed here for clarity.
    1. Time-shift all HPB pressure traces by Equation 24(Equation 2). Original data is shown for Equation 25 mm in Figure 5B, with the same data time-shifted in Figure 5C.
      Note: The time shift is required to allow the interpolation routine to successfully locate the shock front at any given time. This essentially involves aligning the data for each radial array so all the maximum pressures align.
    2. Calculate the radius, Equation 26, and angle, Equation 27 for a given point of interest on the grid, as shown in Figure 5A.
    3. Apply the 1D interpolation to the two HPB arrays closest to the point of interest for the current radius Equation 26 (for Equation 28 the interpolation would use the Equation 29 and Equation 30 arrays).
    4. Interpolate linearly between the 2 pressures based on Equation 31 (again for a Equation 28 the weighting would be 50% of the Equation 29 and 50% of the Equation 30 array calculated pressures).
    5. Calculate the instantaneous load by multiplying the interpolated pressure by the grid spacing (area) to give the load.
    6. Multiply the load by the time step of the sampling to obtain the instantaneous impulse.
    7. Repeat for all locations and times (summing the instantaneous impulse to give the total impulse).
    8. Time-shift the pressure time history for each location based on cubic interpolation of the shock arrival time (Figure 5D).

Figure 5
Figure 5. Interpolation sequence for 2D HPB array. (A) Sign conventions used, (B) original data Equation 35mm, (C) time-shifted data Equation 36 mm, and (D) arrival times for each radial direction16. For a 2D array of bars the pressure time history at any point is dependent on both radial distance and which quadrant the point of interest is located. If the blast were perfectly symmetric then the pressures in (B) would form vertical lines as shown in (C). In (B) it can be seen that the shock front is reaches the 50 mm location on Equation 30 axis first.
Please click here to view a larger version of this figure.

Representative Results

An effectively rigid reaction frame needs to be provided. In the current testing a total imparted impulse of several hundred Newton-seconds needs to be resisted with minimal deflection. An illustration of the rigid reaction frame used is given in Figure 1. In each frame a 50 mm steel 'acceptor' plate has been cast into the base of the cross beams. Whilst not explicitly required, this allows for easy fixing of the load cells / target plate and provides added protection to the face of the concrete beam. The closest scaled distance currently conducted has been 0.15 m/kg1/3.

The current frame has been tested up to 500 Ns, and has 500 mm square columns with a 750 mm deep, 500 mm wide beam spanning the columns as shown in shown in Figure 1A. The critical element in the design is the target plate which is 100 mm thick mild steel, this was estimated to deform 0.3 mm when resisting a 100 g spherical free air blast at 75 mm stand-off (using the LS-DYNA28 routine load_blast4). The construction of the frames was performed by a specialist concrete contractor who provided site equipment and formwork. The factors applied in the design stage very much depend upon the nature of the testing and whether any further safety factors will be applied by the structural engineer. A safety factor of 10 was used in the current work.

Indications of the equipment used have been given in the main protocol sections where appropriate. To provide representative results a single test with 17 HPBs configured in a 2D array as shown in Figure 1B was conducted. In the current work the bars utilized are 3.25 m long with a radius of 5 mm, with the strain gauge being attached 0.25 m from the loaded face as shown in Figure 3A. The spacing of the HPBs in the target was chosen to be 25 mm as shown in Figure 1B.

The test conducted used a 78 g PE4 3:1 squat cylinder buried at 28 mm in saturated Leighton Buzzard sand15. The sand had a bulk density of 1.99 Mg/m3 and moisture content of 24.77%. The stand-off between the finished soil surface and the target plate was 140 mm.

Once the test had been conducted the individual pressure-time histories were run through a simple 5 point moving average smoothing algorithm to remove any high frequency noise from the data. It was noted when collating the data that the 75 and 100 mm bars in the Equation 32 array had not recorded the data properly. This was likely due to the glue of the strain gauge de-bonding from the HPB giving false readings. To compensate for this the data from the 75 and 100 mm Equation 30 bars were used instead. The data from each of the 4 radial arrays are plotted in Figure 6, with the central (0 mm) HPB being common to all plots. In the saturated soil a very clear shock front is seen, with the pressure decaying slowly with radial distance.

The recorded pressure-time histories were then run through the 2D interpolation routine, with the zone of interest being set as a 200×200 mm square (-100 to 100 mm). This zone was subdivided into a 5 mm grid, which was deemed fine enough to capture the shock front propagation across the target plate accurately. Plots of the interpolated pressure acting on the target plate at selected times are shown in Figure 7. The ≈ 20 msec delay in the arrival of the shock front is the time taken for the shock wave to cover the distance between the charge and the target plate. The asymmetric nature of the loading shown in the recorded data (Figure 6) can be seen clearly in the and Equation 32 axes. This is especially clear at Equation 34 msec.

Figure 6
Figure 6.  Recorded pressure-time histories for a single test with a 2D HPB array. (A)Equation 29 (BEquation 30 (C)Equation 37 (D) Equation 38. This figure shows the processed trace for each bar location. The black central bar trace is common to all plots to indicate the arrival of the shock front. The non-continuous nature of the shock front is again clearly visible as there is little overlap between the peak stress zone for each bar. The lower pressure in the 25 mm Equation 38 bar has an interesting effect on shape of the shock front as plotted in Figure 7. Please click here to view a larger version of this figure.

Figure 7
Figure 7. Interpolated pressure distribution plots at specified times16. The horseshoe shape of the shock front can be seen in the t = 0.22-0.23 plots. This is caused by the dip in pressure seen in the 25 mm Equation 38 bar shown in Figure 6. By 0.3 sec after detonation the shock front is almost symmetric along all axes. Please click here to view a larger version of this figure.

Discussion

Using the protocol outlined above the authors have shown that it is possible to get high fidelity measurements of the highly varying loading from an explosive charge, using an array of Hopkinson pressure bars. Using the interpolation routine outlined the discrete pressure-time histories can be transformed into a continuous shock front which is usable directly as the loading function in numerical modelling or as validation data for the output of such models.

When using buried charges the methodology used for preparing the soil containers outlined in protocol section 5 must be checked to ensure sufficient compaction energy is provided to reach the target density. If the target density is not reached then the lift height should be reduced to increase the effectiveness of the compaction. From previous research it has been seen that uniform soil types provide more repeatable test data than tests conducted with well-graded soils15. Tests with soils at full saturation are also prepared with a slightly different methodology as outlined in the authors' previous work15.

For all tests with explosive charges it has been shown in previous studies16,29,30 that the location of the detonator in near-field tests is critical to producing repeatable tests which are free from signal noise. In this respect the detonator should always be placed into the bottom of the charge (furthest from the target plate) so that any fragments from the detonating assembly will not strike the HPBs ahead of the main shock front.

Whilst every attempt is made to ensure the testing is as robust as possible, data loss does still occur. This is usually due to the strain gauges partially de-bonding from the HPBs which can be a particular problem in cold weather (the current apparatus is set up in an unheated building). Great care must also be taken when setting up the break wire as this not only allows the recording of the detonation time but also provides the signal that triggers the data recording. A loss of this signal or error in setup can cause all data from a test to be lost. In regard to data management, the data from the tests are immediately duplicated from the recording computer to a USB drive to ensure no data is lost once the testing is complete.

In the current testing the load cells attach the target plate to the rigid reaction plate and are used to measure the total impulse imparted to the target plate (as the HPBs only cover a limited area). If quantification of only the localized loading (and not global data) is required then the target plate can be mounted directly to the rigid reaction frame.

With the HPB pressure-time histories being only applicable to known points on the target plate an interpolation routine is required to evaluate the pressure-time history for any point on the target plate, and hence to calculate the total recorded impulse.

If only a single radial array has been utilized in the testing, interpolation is still possible by assuming the point HPB loadings to be indicative of the loading at any polar rotation at the same radius on the target plate. Time shifting is also required to interpolate across the discontinuous data (Figure 4A).

The main advantage of using a 2D HPB array is the ability to capture asymmetry in the pressure-time histories. This requires a more complex interpolation routine. In principal this theory can be applied to any number of radial arrays. In the current research this has been limited to four arrays (Equation 29, Equation 37, Equation 30,Equation 38 ) from 0 to 100 mm with the central HPB being common to all (Figure 5A). A total of 17 HPBs have been used in each test.

The interpolation routine in the form presented here assumes that for each pressure-time history there is a single well-defined pressure spike which corresponds with the arrival of the shock front. It can be seen in Figure 6 that for all the bars this is a good assumption. In certain test conditions however this assumption may not be valid and so thought should be given on how best to align the pressure-time histories to allow for the most representative interpolation of pressure.

Modifications can be easily made to cater for different scaled distances (Z) in the current protocol by moving the charge further away from the target plate. Care however must be taken if the scaled distance is brought down below 0.15 to ensure that the loading will not damage the face of the HPBs. The shape of explosive and explosive type may also be changed, with the caveat that the initial modelling done to validate the experimental design will have to be checked.

Divulgations

The authors have nothing to disclose.

Acknowledgements

The authors wish to thank the Defence Science and Technology Laboratory for funding the published work.

Materials

Load Cell RDP RSL0960 This is only indicative, the exact load cell should be able to resolve the required loading
Steel target plate / HPBs Garratts  N/A Fabricated to order
Strain gauge Kyowa KSP-2-120-E4 To use with steel HPBs
Cyanoacrylate Kyowa CC-33-A Check with manufacturer depending on mar material to be used
Digital Oscilloscope TiePie HS4 16-bit Handyscopes  6 used in parallel in current testing
Leighton Buzzard sand Garside sands Garside 14/25 Uniform silica sand 

References

  1. Esparza, E. Blast measurements and equivalency for spherical charges at small scaled distances. Int. J. Impact Eng. 4 (1), 23-40 (1986).
  2. Kingery, C. N., Bulmash, G. ARBRL-TR-02555. Airblast parameters from TNT spherical air burst and hemispherical surface burst. , (1984).
  3. Hyde, D. W. . Conventional weapons program (ConWep). , (1991).
  4. Randers-Pehrson, G., Bannister, K. A. ARL-TR-1310. Airblast loading model for DYNA2D and DYNA3D. , (1997).
  5. Neuberger, A., Peles, S., Rittel, D. Scaling the response of circular plates subjected to large and close-range spherical explosions. Part II: Buried charges. Int. J. Impact Eng. 34 (5), 874-882 (2007).
  6. Xu, S., et al. An inverse approach for pressure load identification. Int. J. Impact Eng. 37 (7), 865-877 (2010).
  7. Pickering, E. G., Chung Kim Yuen, S., Nurick, G. N., Haw, P. The response of quadrangular plates to buried charges. Int. J. Impact Eng. 49, 103-114 (2012).
  8. Bergeron, D. M., Trembley, J. E. Canadian research to characterize mine blast output. , (2000).
  9. Hlady, S. L. Effect of soil parameters on landmine blast. , (2004).
  10. Fourney, W. L., Leiste, U., Bonenberger, R., Goodings, D. J. Mechanism of loading on plates due to explosive detonation. Int. J. on Blasting and Fragmentation. 9 (4), 205-217 (2005).
  11. Anderson, C. E., Behner, T., Weiss, C. E. Mine blast loading experiments. Int. J. Impact Eng. 38 (8-9), 697-706 (2011).
  12. Fox, D. M., et al. The response of small scale rigid targets to shallow buried explosive detonations. Int. J. Impact Eng. 38 (11), 882-891 (2011).
  13. Ehrgott, J. Q., Rhett, R. G., Akers, S. A., Rickman, D. D. Design and fabrication of an impulse measurement device to quantify the blast environment from a near-surface detonation in soil. Experimental Techniques. 35 (3), 51-62 (2011).
  14. Pope, D. J., Tyas, A. Use of hydrocode modelling techniques to predict loading parameters from free air hemispherical explosive charges. , (2002).
  15. Clarke, S. D., et al. Repeatability of buried charge testing. , (2014).
  16. Clarke, S. D., et al. A large scale experimental approach to the measurement of spatially and temporally localised loading from the detonation of shallow-buried explosives. Meas Sci Technol. 26, 015001 (2015).
  17. Hopkinson, B. A Method of Measuring the Pressure Produced in the Detonation of High Explosives or by the Impact of Bullets. Philos. Trans. R. Soc. (London) A. 213, 437-456 (1914).
  18. Rigby, S. E., Tyas, A., Fay, S. D., Clarke, S. D., Warren, J. A. Validation of semi-empirical blast pressure predictions for far field explosions – is there inherent variability in blast wave parameters?. , (2014).
  19. Rigby, S. E., Tyas, A., Bennett, T., Clarke, S. D., Fay, S. D. The negative phase of the blast load. Int. J. of Protective Structures. 5 (1), 1-20 (2014).
  20. Rigby, S. E., Fay, S. D., Tyas, A., Warren, J. A., Clarke, S. D. Angle of incidence effects on far-field positive and negative phase blast parameters. Int. J. of Protective Structures. 6 (1), 23-42 (2015).
  21. Tyas, A., Warren, J., Bennett, T., Fay, S. Prediction of clearing effects in far-field blast loading of finite targets. Shock Waves. 21 (2), 111-119 (2011).
  22. Tyas, A., Watson, A. J. A study of the effect of spatial variation of load in the pressure bar. Meas Sci Technol. 11 (11), 1539-1551 (2000).
  23. Tyas, A., Watson, A. J. An investigation of frequency domain dispersion correction of pressure bar signals. Int. J. Impact Eng. 25 (1), 87-101 (2001).
  24. NATO Standardisation Agency. Procedures for evaluating the protection level of logistic and light armoured vehicles. Allied Engineering Publication (AEP) 55. 2, (2011).
  25. Elgy, I. D., et al. . UK ministry of defence technical authority instructions for testing the protection level of vehicles against buried blast mines. , (2014).
  26. Clarke, S. D., et al. Finite element simulation of plates under non-uniform blast loads using a point-load method: Buried explosives. , (2015).
  27. Rigby, S. E., et al. Finite element simulation of plates under non-uniform blast loads using a point-load method: Blast wave clearing. , (2015).
  28. Hallquist, J. O. . LS-DYNA theory manual. , (2006).
  29. Fay, S. D., et al. Capturing the spatial and temporal variations in impulse from shallow buried charges. , (2013).
  30. Fay, S. D., et al. Measuring the spatial and temporal pressure variation from buried charges. , (2014).

Play Video

Citer Cet Article
Clarke, S. D., Fay, S. D., Rigby, S. E., Tyas, A., Warren, J. A., Reay, J. J., Fuller, B. J., Gant, M. T. A., Elgy, I. D. Blast Quantification Using Hopkinson Pressure Bars. J. Vis. Exp. (113), e53412, doi:10.3791/53412 (2016).

View Video