Experimental Measurement of Specific Impulse Distribution and Transient Deformation of Plates Subjected to Near-Field Explosive Blasts

The shock wave generated from a high explosive detonation can cause significant damage to any objects that it encounters, particularly those objects located close to the source of the explosion. Understanding blast wave development and accurately quantifying its effect on structural systems remains a considerable challenge to the scientific community. This paper presents a comprehensive experimental study into the loading acting on, and subsequent deformation of, targets subjected to near-field explosive detonations. Two experimental test series were conducted at the University of Sheffield (UoS), UK, and the University of Cape Town (UCT), South Africa, where blast load distributions using Hopkinson pressure bars and dynamic target deflections using digital image correlation were measured respectively. It is shown through conservation of momentum and Hopkinson-Cranz scaling that initial plate velocity profiles are directly proportional to the imparted impulse distribution, and that spatial variations in loading as a result of surface instabilities in the expanding detonation product cloud are significant enough to influence the transient displacement profile of a blast loaded plate.


Introduction
When a blast wave interacts with a structure located close to the source of the explosive, the resulting transient blast load is extremely high magnitude, short duration, and highly spatially non-uniform over the face of the target. A considerable challenge to the blast protection community is to understand and control the risks associated with accidental or malicious blast events, and design adequate protective systems to mitigate their effects. Clearly, this

Literature Review
A detailed history of early air-blast experimentation, particularly from the second world war onwards, can be found in Esparza [1]. The data from many of these early studies would later be compiled to form the now wellestablished Kingery and Bulmash (KB) semi-empirical blast predictive method [2]. Close to the charge, direct measurements of the blast parameters were either "nonexistent or very few" [1] and the near-field semi-empirical predictive data were inferred from non-direct measurements [3] or rudimentary numerical analyses [4]. Whilst KB blast parameter predictions have been shown to be accurate for far-field, geometrically simple scenarios [5], there is currently a lack of definitive experimental validation data in the extreme near-field (a few radii from the centre of the charge). As such, there remains a great level of uncertainty on the exact form and magnitude of the loading imparted to a target located close to the source of an explosion, as highlighted in the modelling work of Shin et. al [6] where differences of >400% were reported between the numerical and predicted reflected pressure at short distances from the explosive.
Since the development of the KB predictive method there have been efforts to experimentally measure the blast load in extreme near-field conditions. One method involves measuring the residual momentum of -and hence inferring the impulse imparted to -small, rigid metal plugs embedded within a larger target surface [7,8], however this approach cannot be used to give a temporal description of the blast load. An alternative, developed in 1914 by Bertram Hopkinson, is the apparatus now known as the Hopkinson pressure bar (HPB) [9], consisting of a length of cylindrical bar which propagates an elastic stress pulse along its axis to be recorded by sensitive equipment situated a safe distance from the loaded end. Whilst it is now more commonly used in its 'split' form for high strain-rate material testing [10], the HPB is still a valuable tool for measuring highmagnitude, short-duration loading [11][12][13][14][15][16]. HPBs are used in this study at UoS to record the spatial and temporal distribution of loading acting on a rigid target located close to an explosive.
There is a large body of experimental work on the dynamic response of plates subjected to large dynamic loading arising from blast and impact events. Typically the primary metric used to assess plate performance under extreme loading is some measure of peak dynamic or residual deformation at the centre of the plate [17][18][19][20][21][22][23][24][25][26][27]. Whilst this information may enable researchers to develop relationships between basic parameters, it is often more desirable for validation of numerical modelling approaches to have a more detailed description of the transient displacement profile of the plate. Digital image correlation (DIC) was first introduced in 1983 [28] and has recently been used to more accurately study fullfield transient deformations of blast loaded plates [29][30][31][32][33][34][35][36][37][38][39][40][41][42].
It is clear that whilst the topics of loading and deformation have been studied more extensively in recent years, our knowledge of exactly how the blast load develops on the loaded face and subsequently drives the displacement of a target plate is still inhibited by a lack of repeatable and well-controlled experimental data. There is a need, therefore, to present a combined study into the detailed measurement of blast loading and target deformation from extremely near-field blast events.

Hopkinson-Cranz Scaling
Hopkinson-Cranz scaling [43,44] is based on fundamentals of geometrical similarity and can be used to relate the blast parameters from one set of physical experiments (or simulations) to another set conducted at a different scale. Two explosions can be expected to give identical blast waves 1 at distances which are proportional to the cube-root of the respective energy release [45], i.e. the blast pressure profile at a distance of R from an explosive mass W will be similar to the blast pressure at a distance of KR from a mass of K 3 W . The scaled distance of a blast event, therefore, is defined as the distance from the blast divided by the cuberoot of the charge mass, i.e. Z = R/ 3 √ W . It can be seen that the two examples described previously will have an identical scaled distance, Pressures and velocities do not change between scales, however arrival times, durations and impulses are scaled by the length scale factor, K. The Kingery and Bulmash predictive method [2], for example, utilises a form of Hopkinson-Cranz scaling to present scaled distance-parameter relationships, where scaled and de-scaled pressure parameters are identical, and arrival times, durations and impulses are scaled by the cube-root of the charge mass. In this study, the principles of Hopkinson-Cranz scaling are used to compare the results from blast experiments conducted at different length scales but with identical scaled geometries.

'Characterisation of Blast Loading' Apparatus
Blast load distributions were measured in experiments conducted at the University of Sheffield Blast & Impact 1 Assuming identical charge shapes and ambient conditions   [46]. The CoBL apparatus, shown in Fig. 1(a), comprises a pair of stiff, massive, fibre and bar reinforced concrete frames spaced 1 m apart, with each frame comprising two 500 mm square columns with a 750 mm deep, 500 mm wide concrete beam spanning horizontally between the two columns. A 100 mm thick steel target plate is underslung from the soffits of the horizontal beams and acts as a nominally rigid boundary to reflect the shock wave and detonation products impinging on the target after detonation of an explosive some distance beneath the centre of the plate. The plate is 1.4 m in diameter to negate the effect of blast wave clearing around the target edge [47].
The target plate is drilled through its thickness to allow 10 mm diameter, 3.25 m long EN24(T) steel HPBs to be mounted and set with their loaded faces flush with the underside of the target plate. A total of 17 bars were used; one central bar and four bars located at each radial offset of 25, 50, 75 and 100 mm from the plate centre, with the bar naming convention in this article following the coordinate axes shown in Fig. 1(b).
Kyowa KSP-2-120-E4 semi-conductor strain gauges were mounted in pairs on the perimeter of each HPB at 250 mm from the loaded face, in a Wheatstone-bridge circuit to neglect any bending effects and to ensure that only the axial strain component was recorded. Strain data were recorded using 14-Bit digital oscilloscopes at a sample rate of 3.125 MHz and were triggered via a voltage drop in a separate breakwire channel.

Digital Image Correlation Pendulum
Plate deformations and global impulse were measured in tests conducted at BISRU at the University of Cape Town. The blast pendulum described in Ref. [49], was recently modified to include stereo high speed video (HSV) capabilities [40]. The pendulum consists of an I-beam suspended from four cables with a clamping frame located at the front of the pendulum to constrain the test plates, and added counterbalance at the rear of the pendulum (Fig. 2). A mild steel shroud covers the pendulum to protect the cameras and lighting system housed within. Part of the web of the I-beam had been removed so as to not inhibit the cameras' view (and illumination) of the rear of the test plate. The cameras were positioned appropriately for the focal distance of the lenses selected, and were mounted on an extruded aluminium rail system which was sat on nylon blocks to reduce the transmission of vibration to the cameras.
The test plates were made from 400×400×3 mm thick Domex 355MC (7830 kg/m 3 ), a high-strength, hot rolled, low alloy steel [40], with an exposed circular area of  300 mm diameter when loaded in the clamping frame. The plates were fully clamped around the perimeter, with constrained zero in-plane and out-of-plane displacement and constrained rotations at the boundary.

Transient plate deformation measurements
Digital image correlation (DIC) [28] was used to measure the transient deformation profile of the plates using two IDT vision NR4S3 high speed cameras. The cameras were separated at an included angle of approximately 30 • and were fitted with 35 mm fixed focus lenses. A field-ofview of 1024×76 pixels was selected in order to cover a full-width strip along the centreline of the plates (an approximate region of 300×25 mm). HSV data were recorded at 30,000 fps with the exposure set to 31 μs. The desired illumination was achieved with the use of two custom-built LED lights, focussed on the rear of the test plate. Light diffusers were used to ensure uniform illumination across the field-of-view of the cameras. The rear surface of each plate was first sprayed with a thin layer of white primer before a black speckle pattern was added, generating a random pattern for use in the DIC software. The rear surface of each plate was thoroughly cleaned, abraded, and degreased with acetone before painting to allow for better adhesion of the paint to the plate.
Dantec Dynamics Istra 4D DIC software was used to process the HSV data. Prior to each test a calibration target was filmed and the data was used to generate intrinsic and extrinsic calibration values for each test, ensuring consistent accuracy across the entire test series. These calibration values were then imported into the DIC software enabling full-field out-of-plane transient displacement profiles of the plate to be generated post-test. Further information on the DIC method used in the current testing is available in Ref. [40]. Data were extracted from a gauge line, defined in the software using a straight line drawn through the centreline of the test plate. The transient centreline displacements were then exported into MatLab for postprocessing. Any rigid body motion of the cameras and framing system was automatically removed from the processed results.

Global impulse measurements
The motion of the pendulum was recorded using a CP35MHT80 laser displacement sensor (50 μm resolution) which was focussed on the rear of the pendulum. Displacement data were recorded by a digital oscilloscope which was triggered by the initial movement of the pendulum and recorded at a sample rate of 100 kHz, which was found to be sufficient to fully capture the swing of the pendulum post detonation.
The impulse imparted to the test plate, I , is given as where m is the mass of the pendulum, and v is the initial velocity of the pendulum, calculated from where T is the natural period of the system, taken as the average period under forced vibration measured separately before each test series, and β is the damping coefficient which can be calculated from knowledge of the peak forwards displacement (in the direction of the blast), y 1 , and peak backwards displacement (towards the source of the blast), y 2 , The values of peak forwards and backwards displacement are determined for each test from a 5 th -order polynomial curve fit to the recorded pendulum motion, with separate curves defined for forwards and backwards displacement cycles for greater accuracy.

Explosive Charge Geometries and Test Plan
Seventeen tests were conducted in total: seven tests using the CoBL apparatus, as detailed in the UoS test programme   Table 1; and ten tests using the digital image correlation pendulum, as detailed in the UCT test programme in Table 2.

Charge masses
PE4 explosive charges were used throughout. In tests 1-3 at UoS, 100 g (±0.05 g) spherical PE4 explosive charges were hand-pressed into bespoke 3D-printed moulds, ensuring a consistent and regular spherical shape. In tests 4-7 at UoS, 78 g PE4 charges were formed into cylinders with diameter:height ratio of 3:1 by hand-pressing them directly into 3 mm thick PVC containers with internal dimensions exactly matching the required dimensions of the charges. The charges in test 4-6 were kept in the PVC containers (consistent with previous testing [48]), and the charge in test 7 was removed from the container to provide comparative data on the container's influence on blast load development (see "Effect of charge casing").
In the UCT tests, 50 g (±0.05 g) PE4 chargers were handrolled into spheres (tests 8, 9, and 12-14), or hand-pressed into cylindrical moulds (tests 10, 11, and 15-17) with diameter:height ratio of 3:1 to match the geometry of the UoS tests. Each charge was measured with callipers after being shaped to ensure charge size consistency. Five tests were performed for each charge configuration: two commissioning tests without DIC, and three tests with DIC (Table 2).
Charge geometries and stand-off distances for the UCT tests were designed to be similar to UoS tests using Hopkinson-Cranz scaling ("Hopkinson-Cranz Scaling").

Charge placement
The charges in tests 1-3 were suspended directly under the centre of the target plate on a 'drumskin' comprising a glassfibre weave fabric (density 25 g/m 2 ) held taut in a steel ring, set on adjustable struts mounted in the base of the test arena. The charges in tests 4-7 were placed on a small timber prop, sat inside an empty steel container to match the geometry of previous testing with buried explosives [48,50]. The timber prop was machined to give the required stand-off between the charge surface and the target, and was 3 mm longer in test 7 to account for omission of the PVC case. In all UoS tests the charges were aligned with the plate centre using an alignment laser, and the clear stand-off distance from the charge surface to the plate (55.4 mm tests 1-3, 168.0 mm tests 4-7) was measured and confirmed to be within tolerance (±0.5 mm).
The UoS charges were detonated using Nitronel MS 25 non-electronic shock-tube detonators (700 mg PETN) inserted through the bottom face of the charge to a depth marked on the detonator corresponding to the charge radius (tests 1-3) or half charge height (tests 4-7). A breakwire was wrapped around the detonator to give an accurate reading of the time of detonation and hence serves as a trigger for the HPB channels.
In the UCT tests, the charges were affixed to a thin polystyrene bridge which spanned 250 mm onto polystyrene legs which were machined to the exact stand-off distance plus charge height, as the rear face of the explosives was attached to the polystyrene bridge so as to not obstruct the propagation of the blast wave between the charge and the target. The bridge legs were attached to the target plate using double-sided tape and were positioned according to markers drawn on the plate using a template to ensure consistent charge placement. M2A3 electric detonators were inserted through the polystyrene bridge into the centre of the charges from the rear, to ensure repeatable detonator placement which was consistent with the UoS tests (detonators were inserted to a depth equal to the charge radius in tests 8,9, and 12-14, and a depth equal to half charge height in tests 10,11, and 15-17). As with the UoS tests, the required insertion depth was marked on the detonator prior to placement to ensure consistency. A breakwire was formed using a thin strip of foil attached to the back of the polystyrene bridge which was used to trigger the cameras. Bare charges were used throughout, based on the negligible influence of the cylindrical charge casing as demonstrated in "Effect of charge casing".

Direct Load Measurements
Spherical charge tests to reach the gauge location from the loaded face of the HPB. Aside from oscillations in the data caused by Pochhammer-Chree dispersion [51], the blast pressures appear to resemble the characteristic 'Friedlander' waveform [52], in that there is a sudden rise to peak pressure followed by an exponential decay back to ambient conditions, with recorded positive phase durations of ∼0.05-0.07 ms. As the blast event is located within the extreme near-field regime, negative phase effects are either negligible or non-existent [53]. Dispersion effects can be seen to 'round off' the leading edge of the stress pulse, with the peak pressure in the central (0 mm) bar occurring on the second oscillation, rather than with arrival of the stress pulse as is the case with all other bars. Some higher frequency components of the signal (which travel at a lower velocity relative to the lower frequency components [54]) can be seen to arrive towards the end of the positive phase. Whilst this demonstrates that there may be significant high frequency features associated with the sharp rise of the pressure pulse, current frequency-domain dispersion correction methods are limited to frequencies of 250 kHz [55] for the bar radii used in this study. Transients of < 4μs therefore cannot yet be resolved and details of the blast load immediately following arrival of the shock wave cannot be extracted from the current measurement technique. The specific impulse, determined from cumulative temporal integration of the pressure signals, is unaffected by dispersion, and decays from an epicentral value of 3.87 MPa.ms to a value of 0.82 MPa.ms at a radial offset of 100 mm (approximately 4 charge radii) from the plate centre.  Figure 4 shows the compiled peak pressure and peak specific impulse for tests 1-3. Generally there is a good degree of test-to-test repeatability for this charge size and stand-off distance, with all values of peak pressure and peak specific impulse occupying a narrow band either side of the mean distribution.
Cylindrical charge tests Figure 5 shows example pressure-time and specific impulsetime histories from test 5, again with 50 μs subtracted from the time datum. Compiled peak pressure and peak specific impulse for tests 4-6 are shown in Fig. 6.
Generally, the compiled blast parameters for the cylindrical tests demonstrate a larger degree of variability when compared with the spherical tests. For example, the +y 25 mm bar in test 5 recorded a peak pressure of 243 MPa, which is 58% greater than the mean peak pressure for all bars at 25 mm radial ordinate. The peak specific impulse recorded at the +y 25 mm bar in test 5 was 3.63 MPa.ms, which is 29% greater than the mean peak specific impulse for all bars at 25 mm radial ordinate. This variability is evidence of the growth of local irregularities due to Rayleigh-Taylor [56,57] and Richtmyer-Meshkov [58,59] instabilities on the surface of the expanding detonation product cloud. It is hypothesised that the blast wave from the (extremely near-field) spherical charges in tests 1-3 impinges on the target before these instabilities have formed to any significant degree, whereas in the cylindrical tests (4-6), where there is a greater physical and scaled distance between the charge and the target, these instabilities become more prominent [15]. Figure 7 shows the compiled peak specific impulses (and mean distribution) from test 7 compared to the mean distribution from tests 4-6. Whilst the epicentral impulse for test 7 appears to be slightly below the average from tests 4-6, it is worth noting that this value comes from only one data point and hence some degree of experimental spread should be expected, given the growth of instabilities mentioned previously. There is very little difference between the test 4-6 impulse distribution and test 7 impulse distribution at distances greater than 25 mm from the target centre. Therefore it can be said that, for this geometry, the influence of a 3 mm thick PVC charge case is negligible.

Interpolated impulse distribution
A spline interpolant was fitted to the values of mean peak specific impulse at each radial ordinate for clarity of presentation and to aid future numerical model validation. Three conditions were applied to the fit: zero gradient at the plate centre; zero gradient and zero impulse at an arbitrary large radial offset from the plate centre; and non-negative peak specific impulse at any radial ordinate. These conditions ensured the spline interpolant was physically valid, i.e. radially symmetrical and monotonically decreasing with increasing distance from the plate centre. Figure 8 shows fitted interpolants for both charge arrangements. Integrating these curves with respect to area gives   10(c) show the displacement profile of the plates at the time of peak deflection. The DIC software was unable to determine any displacements between 6.6-6.9 ms after detonation in test 15 due to debris obscuring the cameras' view of the test plate. However, this data drop-out occurs after the peak deflection is reached and hence the results are unaffected.

Dynamic deformation
Generally the spherical peak deflection profiles and histories are all in good agreement, with peak deflections occurring near the centre for all tests. The cylindrical tests appear less consistent, with test 15 demonstrating a considerably lower peak deflection (which also occurs some 15 mm from the plate centre), and a marked difference   Fig. 11 Residual plate profiles and peak dynamic deflections for UCT spherical tests in the deflected shape at peak deflection. The observed 15 mm offset is significantly outside the tolerance on charge geometry and placement, and hence this discrepancy is likely due to an inherent feature of the explosive itself, rather than experimental error. This will be discussed further in "Discussion".

Imparted impulse and residual deformation
The deformed plates were scanned using a NextEngine    Tables 3 and 4 along with the imparted impulse determined from the recorded pendulum motion. Bonorchis and Nurick [60] demonstrated, through numerical analysis, that for a similar test arrangement to that desrcibed in this article, the impulse acting on the exposed area of the plate was typically ∼67% of the total impulse recorded by the UCT blast pendulum. This factor was found to be insensitive to scaled distance, and hence is used in this article to estimate the total impulse acting on the pate. The mean recorded impulses are 66.0 Ns and 67.1 Ns for the spherical and cylindrical shots respectively. That the impulses are comparable despite the cylindrical charges being places three times further (clear distance) from the target plate demonstrates the increased directionality caused by the charge shape; an effect recently discussed by Nurick et. al [61]. The spread of experimental impulses is slightly lower for the spherical charges than for the cylindrical charges; all spherical recorded impulses are within ±3.5% and all cylindrical recorded impulses within ±6.0% of the respective means. It cannot be said for certain whether this increased variability for the cylinders is significant owing to a small dataset (only five global impulses per configuration), however it is clear that the variability in global impulse uptake is considerably less than the variability in specific impulse seen in "Direct Load Measurements".
The mean values of peak residual and peak dynamic deflections are 15.86 mm and 20.54 mm for the spheres and 17.91 mm and 21.96 mm for the cylinders. The peak residual and peak dynamic deflections recorded in test 15 are considerably lower than all other cylindrical tests, despite the global impulse closely matching the mean value from the test series. It is suggested that this discrepancy is caused by non-coaxial impingement of the shock wave and detonation products, likely due to instabilities emerging from the detonation product/air interface, which gives reason for the lower magnitude and non-central peak (dynamic and residual) deformation. This is justified with reference to Fig. 6(b) where the maximum peak specific impulse at 25 mm from the plate centre in test 5 was 21% greater than the mean central bar peak specific impulse and 16% greater than the peak central bar specific impulse for that test. It is possible that a similar situation occurred in test 15 and the maximum peak specific impulse was acting some distance from the plate centre. Figure 13 shows the displacement-time history along the centreline of the plate, recorded from test 12. Figure 14 shows the velocity-time history along the centreline of the plate for the same test, assuming a linear gradient between displacement profiles. It can be seen that peak displacement *denotes impulse acting on the exposed area of the plate has been estimated by factoring the total impulse by 0.67, after [60]  is reached approximately 0.73 ms after detonation, whereas peak plate velocity is reached approximately 0.06 ms after detonation. At this time, the plate has displaced a relatively small amount compared with the peak value, and deformation is concentrated towards the central ∼100 mm diameter region of the plate. This suggests that the work done by the blast load is small and hence the system is responding largely in the impulsive domain. Two sets of flexural waves can be observed in the velocity-time history in Fig. 14: one set beginning shortly after arrival, propagating inwards from the edges of the plate; and one set beginning shortly after peak velocity is reached, propagating outwards from the centre of the plate. This behaviour is consistent with flexural waves generated by uniformly distributed and central point loads respectively, and are indicative of the specific impulse distribution seen in Fig. 8(a). Similar behaviour was seen for the cylindrical charge tests, however the results are not repeated here for brevity.

Plate velocity profiles
The inward travelling flexural waves can be seen to cross at the plate centre at approximately 0.39 ms after detonation, indicating a propagation velocity of approximately 400-450 mm/ms. After crossing, the magnitude of the waves appears to diminish and they are subsequently difficult to distinguish at the current plotting resolution. The outward travelling flexural waves reach the edge of the plates at approximately 0.4-0.45 ms after detonation where they are reflected inwards and cross at the plate centre at approximately 0.68 ms after detonation. This repeated crossing of flexural waves causes the high frequency oscillatory behaviour seen in Figs. 9(a) and 10(a).
Peak velocity profiles were extracted from postprocessing of the DIC data and are shown in Figs. 15

Scale Factors
Using Hopkinson-Cranz scaling [43,44], we can say that where x, i, and W are distance from the centre of the plate, specific impulse, and charge mass, and the subscripts UoS and UCT refer to quantities used in the University of Sheffield and University of Cape Town tests respectively. In order to compare the UCT results to the UoS results, therefore, we must multiply distances and specific impulses by the cube-root of the mass scale factor: 3 √ 100/50 = 1.26 for the spherical tests and 3 √ 78/50 = 1.16 for the cylindrical tests respectively. Furthermore, since total impulse is a product of specific impulse and area, total impulses can be compared by dividing or multiplying by the mass scale factor: 2.00 for spherical tests and 1.56 for cylindrical tests.

Total and Specific Impulse Comparison
An imparted impulse results in an equivalent change in momentum. Under impulsive loading conditions, therefore, Peak specific impulse (MPa.ms) Test where x is distance and i is specific impulse as defined previously, and ρ and t are density and thickness of the plate: 7830 kg/m 3 and 3.00 mm respectively. The average factored impulses for the UCT spherical and cylindrical charge tests are 44.2 and 45.0 Ns respectively (Tables 3 and 4). When expressed at UoS scale, by multiplying by the respective mass scale factors, this gives 88.4 and 70.2 Ns, which is in very good agreement with the values of 83.9 and 71.3 Ns given in "Interpolated impulse distribution". Figure 17 shows the inferred specific impulses from the UCT spherical tests calculated from equation (6), expressed at UoS scale using equation (5), compared to the specific impulses directly measured from the UoS spherical tests. Figure 18 shows the same for the cylindrical tests. Here, markers for the UoS test results have been duplicated either side of the plate centre for comparative purposes.
For the spherical tests (Fig. 17), the inferred specific impulses closely match the measured values. This confirms that Hopkinson-Cranz scaling is valid for the scales used in the current testing. While the DIC velocity measurements do not have the same temporal resolution as the UoS Hopkinson bar results, the inferred specific impulses correlate very well and show the same form and magnitude as the directly measured load. For the cases studied in this paper, as there is little difference between the loading measured on a rigid reflecting surface (UoS tests) and a deformable one (UCT tests), it can be said that fluidstructure-interaction [62][63][64] has negligible influence on the load developed on the reflecting surface.  The directly-measured cylindrical specific impulses demonstrated considerably more spread than those in the spherical tests. This increased variability is also borne-out in the inferred specific impulse distributions from the DIC measurements (Fig. 18). Generally, the inferred specific impulses for the cylinder tests lie within the range of recorded specific impulses, and are in particularly good agreement between ±25-50 mm from the plate centre. When viewed with this dataset, test 15 appears less of an outlier and better fits within the general range of recorded data. In almost all cases, the inferred specific impulse begins to drop below the recorded values at ±100 mm from the UoS-scale plate centre. This is due to the fact that the loading is not truly impulsive, and when considering initial velocities there is a non-negligible delay between the loading being applied in the plate centre and the loading being applied at the periphery of the plate. The velocity profiles in this article are extracted from one particular instant in time (when the maximum velocity is acting at the centre of the plate), rather than the time-invariant peak velocity at any given point on the plate, however the influence that this has on the inferred impulse distribution appears to be minimal.

Summary and Conclusions
This paper presents the results from two experimental test series conducted at different scales, and investigated the loading distribution and subsequent dynamic response of plates subjected to near-field explosions. Direct load measurements were performed at the University of Sheffield (UoS), UK, and plate deformation measurements were preformed at the Blast Impact & Survivability Unit (BISRU) at the University of Cape Town (UCT), South Africa. In each series, two different charge configurations were tested: spherical and cylindrical (3:1 diameter:height) PE4 charges. The spherical and cylindrical charges were 100 g and 78 g respectively in the UoS tests, and 50 g for both spheres and cylinders in the UCT tests. Direct load measurements were performed on an effectively rigid target plate, whereas deformation measurements were performed using 3 mm thick Domex 355MC plates, with fully clamped boundaries and an effective exposed circular area with a 300 mm diameter.
A significant increase in the directionality of loading was observed for the cylindrical charges. The imparted impulses from spherical and cylindrical charges were similar, despite the cylinders being three times further (clear scaled distance) from the target plate. Whilst the total impulses integrated over the instrumented area (UoS) and entire plate surface (UCT) were highly repeatable for both spherical and cylindrical charges, the specific load distributions demonstrated a considerably higher level of variability in the cylindrical tests. This is hypothesised to be due to a greater relative distance from the charge to the target and hence greater propensity for Rayleigh-Taylor and Richtmyer-Meshkov instabilities to develop.
The plate deformation measurements from the cylindrical charges also demonstrated a higher level of variability relative to the spherical charge plate deformation measurements. This suggests that spatial variations in loading as a result of surface instabilities in the expanding detonation product cloud are significant enough to influence the transient displacement profile of a blast loaded plate.
Initial velocity distributions were calculated from UCT DIC data. Conservation of momentum was used to infer the distribution of specific impulse, and Hopkinson-Cranz scaling was used to express the UoS and UCT tests at the same scale. The inferred impulses from the UCT tests showed excellent agreement with the directly-measured specific impulse distributions from the UoS tests.

Data access statement
The data presented in this publication can be obtained on request by contacting sam.rigby@sheffield.ac.uk.