# Performance Study of Impact-Point Estimation Systems Using Two Types of Acoustic Sensor Arrangement Methods

Authors: Shin Hur; EunJoong Lee; Jongseong Won; Duck-Gyu Lee; Kyungjun Song; Kyusik Park
Journal: Journal of Sensors
 (2019)

---

## Abstract

In this study, a shooting experiment is carried out with an acoustic sensor array that can detect shock wave at high-sound pressures and withstand the shock wave of the bullet. Delta- and H-type of acoustic sensor units are constructed. The measured values of the temporal arrival differences estimated based on the distance between the acoustic sensors are then compared with the theoretical calculations according to the types of bullets and shooting distance. In addition, an impact-point estimation system (IPES) is fabricated, which comprises delta- and H-type sensor units, a target, signal processing device, and an impact point estimation algorithm. The impact points on the target caused by the shooting experiment are compared with those calculated using the impact point estimation algorithm. Finally, the mean absolute errorMAE between the actual and the calculated impact-point coordinates is calculated, and the performances of the delta- and H-type IPES are evaluated.


---

## Body

## 1. Introduction

When a bullet is fired from a gun owing to the occurrence of the deflagration inside it, the air around the bullet is compressed and the bullet generates a cone-shaped shock wave. The generated blast wave then propagates at the sonic velocity in the radial direction. As the shock wave is generated in a clear N-shape within a short time interval, the position of the bullet can be identified based on the time difference of the arrival (TDoA) of the shock wave. Acoustic sensors capable of measuring shock wave signals must withstand sound pressures greater than 140 dB. They must also exhibit uniform detection performance regardless of temperature and humidity changes. In conventional shooting exercises, whenever a shooter shoots, the person’s visual acuity or the assisting telescope confirm the location at which the bullet has penetrated the target plate. However, it is difficult and time consuming to determine the location of an impact point because the shooter uses the naked eye or the telescope to determine the location of the impact point. Therefore, various countries, including the US, UK, Switzerland, Sweden, and others, have developed impact-point estimation systems (IPES), such as the location of miss and hit (LOMAH) system, to support shooting training programs. The IPES can measure the impact of a bullet and display the bullet’s impact point in real time. The system allows shooters to identify the point of impact with a display at a nearby location, which can shorten the time of shooting training and improve the safety of the training program. Training commanders can also manage training information and perform real-time assessments to summarize the shooting results and monitor the progress of trainees. A typical IPES consists of a shock-wave detection unit based on an acoustic a sensor array, a signal processing device, an impact-point estimation algorithm, a computer, and a monitor. In general, the shock-wave detection unit possesses signal processing functions and detects the shock wave generated by the bullet. To accurately estimate the impact point, the placement method and number of acoustic sensors must be determined, and an effective impact-point estimation algorithm must be used. Previous studies related to this current investigation included theoretical and experimental studies of the shock wave and gunshot of the bullet generated by shooting, and the study of the time delay between the generation of the shock wave and the gunshot [1, 2]. Furthermore, theoretical and experimental studies have been published that allow the detection of the position of a sniper with the use of fire detectors worn by soldiers [3, 4]. However, no prior studies compare the impact point performance with different bullet types and sensor arrangement methods using the IPES.This study selects an acoustic sensor that is capable of sensing the shock wave of a bullet and compares the calculated and experimental values of TDoA of the shock wave caused by the distances between acoustic sensors. Furthermore, the IPES are fabricated with delta- and H-type of acoustic sensor units and the impact-point location of bullets on the target are obtained from an actual shooting experiment and compared with those calculated using the impact-point estimation algorithm. The performances of the IPES are evaluated by comparing the errors of the impact-point coordinates between the shooting experiment and the impact-point estimation calculations.


### 1.1. Theory of Shock Wave

The speed of an object traveling faster than the speed of sound in air can be expressed using the Mach number,M=V/c, as the ratio of the speed of sound c and the velocity of the object V. Figure 1 shows the shock wave generated by a bullet traveling at a velocity higher than the speed of sound from the compression of the air. This shock wave has a conical shape with an angle θ relative to the path of the traveling object, and the wave propagates in the vertical direction of the shock wave at the speed of sound [5].Figure 1
Generation of a cone-shaped shock wave for a bullet traveling with a velocity V.In this case, the propagation angleθ of the shock wave is equal to arcsin1/M. The shock wave measured through an acoustic sensor has a clear N-shape. The time interval T from the highest to the lowest points of the N-shape is expressed as follows [1]:
(1)T=1.82dcMxl1/4where d is the diameter of the bullet, x is the vertical distance between the trajectory of the bullet and the acoustic sensor, and l is the length of the bullet. The speed of sound c is expressed according to 331.5+0.6×t, indicating that the temperature t has a significant impact on the formation of the shock wave. Figure 2 shows a shock wave with a typical N-shape. The time interval T of the shock wave was measured as the peak-to-peak interval of the measured voltage shock wave signal.Figure 2
N-shape form of shock wave signal.Figure3 shows a schematic for the calculation of the shock and blast wave arrival times. The distance between each acoustic sensor and the impact point can be expressed as x, the travel distance from the shock wave generated in the path of the bullet to the acoustic sensor as x/cosθ, and the travel time as x/c·cosθ. In this case, the distance from the time of shooting to the shock wave detection time by the acoustic sensor is expressed as L−xtanθ using the shooting distance L. The corresponding travel time is L−xtanθ/V. The sum of the travel time of the bullet from the time of shooting and the travel time of the shock wave that reaches the acoustic sensor is determined as L−xtanθ/V+x/c·cosθ. Using this expression, the shock wave arrival times at the six sensors can be calculated for an arbitrary shooting distance. The travel distance of the blast wave from the shooting point to the acoustic sensor can be expressed as x2+L2, and the corresponding travel time is expressed as x2+L2/c. When the shooting takes place at a distance of 100m, the shock wave arrival time is 138ms at the temperature of 20°C, the bullet velocity is 920m/s, the bullet diameter is 7.62mm, the bullet length is 15mm, and the vertical distance x between the bullet path and the sensor is 1m. The time difference of arrival (TDoA) between the shock wave and the blast wave measured by the acoustic sensor, i.e., L−xtanθ/V+x/c·cosθ−x2+L2/c, increases as the shooting distance increases [1].Figure 3
Schematic explaining the calculation of the shock and blast wave arrival times.


## 1.1. Theory of Shock Wave

The speed of an object traveling faster than the speed of sound in air can be expressed using the Mach number,M=V/c, as the ratio of the speed of sound c and the velocity of the object V. Figure 1 shows the shock wave generated by a bullet traveling at a velocity higher than the speed of sound from the compression of the air. This shock wave has a conical shape with an angle θ relative to the path of the traveling object, and the wave propagates in the vertical direction of the shock wave at the speed of sound [5].Figure 1
Generation of a cone-shaped shock wave for a bullet traveling with a velocity V.In this case, the propagation angleθ of the shock wave is equal to arcsin1/M. The shock wave measured through an acoustic sensor has a clear N-shape. The time interval T from the highest to the lowest points of the N-shape is expressed as follows [1]:
(1)T=1.82dcMxl1/4where d is the diameter of the bullet, x is the vertical distance between the trajectory of the bullet and the acoustic sensor, and l is the length of the bullet. The speed of sound c is expressed according to 331.5+0.6×t, indicating that the temperature t has a significant impact on the formation of the shock wave. Figure 2 shows a shock wave with a typical N-shape. The time interval T of the shock wave was measured as the peak-to-peak interval of the measured voltage shock wave signal.Figure 2
N-shape form of shock wave signal.Figure3 shows a schematic for the calculation of the shock and blast wave arrival times. The distance between each acoustic sensor and the impact point can be expressed as x, the travel distance from the shock wave generated in the path of the bullet to the acoustic sensor as x/cosθ, and the travel time as x/c·cosθ. In this case, the distance from the time of shooting to the shock wave detection time by the acoustic sensor is expressed as L−xtanθ using the shooting distance L. The corresponding travel time is L−xtanθ/V. The sum of the travel time of the bullet from the time of shooting and the travel time of the shock wave that reaches the acoustic sensor is determined as L−xtanθ/V+x/c·cosθ. Using this expression, the shock wave arrival times at the six sensors can be calculated for an arbitrary shooting distance. The travel distance of the blast wave from the shooting point to the acoustic sensor can be expressed as x2+L2, and the corresponding travel time is expressed as x2+L2/c. When the shooting takes place at a distance of 100m, the shock wave arrival time is 138ms at the temperature of 20°C, the bullet velocity is 920m/s, the bullet diameter is 7.62mm, the bullet length is 15mm, and the vertical distance x between the bullet path and the sensor is 1m. The time difference of arrival (TDoA) between the shock wave and the blast wave measured by the acoustic sensor, i.e., L−xtanθ/V+x/c·cosθ−x2+L2/c, increases as the shooting distance increases [1].Figure 3
Schematic explaining the calculation of the shock and blast wave arrival times.


## 2. Results and Discussion

### 2.1. Experimental Measurements

A typical muzzle explosion produces shock waves with a sound pressure level (SPL) ≥ 140 dB. Thus, acoustic sensors used for measuring these shock waves must be able to withstand these sound pressures [6]. Therefore, several fire tests were conducted on various acoustic sensors. Based on these tests, the acoustic sensor that was capable of withstanding shock waves was selected. As shown in Figure 4(a), the structure of the acoustic sensor is composed of a protection grid, an electret condenser microphone (ECM), a preamplifier, housing, and an SMB connector. Figure 4(b) shows the acoustic sensor (HJ06, SISPIA Co., Korea). Each of the acoustic sensors should also have uniform sensitivity and phase characteristics to predict accurately the impact point of the bullet. The sensitivity and phase characteristics of the acoustic sensors were measured in an anechoic chamber according to the standard measurement method. The sensitivity of the acoustic sensors was approximately -46.0 dB in the frequency band of 100 Hz to 10 kHz, and the phase was approximately within ±3° at 1.0 kHz. Figure 5 shows the sensitivity and phase measurement outcomes from the six acoustic sensors, thus indicating the uniform sensitivity and phase characteristics with respect to frequency. The performance specifications of the acoustic sensors are listed in Table 1. The IPES based on the delta- and H-type sensor units was constructed with the use of the six acoustic sensors with uniform sensitivity and phase value.Figure 4
Acoustic sensor. (a) Schematic of structure and constituent components. (b) Photograph of acoustic sensor.
(a)
(b)Figure 5
Sensitivity (solid lines) and phase (dotted lines) measurements obtained from the acoustic sensors.Table 1
Measurement results of the acoustic sensor.
Sensitivity
-46.0 dB relative to 1 V/pa ± 3 dB
Frequency response
± 2 dB, 100 Hz to 10 kHz
Phase response
± 3°, 100 Hz to 3 kHz± 5°, 3 kHz to 5 kHzTo measure the shock wave of a bullet and to test the performance of the IPES, an experiment was performed in an actual shooting test site. Bullets with diameters equal to5.56mm and 7.62mm were used. Prior to the experiment, the bullet velocity was measured at 50 m, 100 m, 150 m, and 300 m from the shooting point with a speed measuring device (Doppler radar SL-528PE, Denmark) based on the Doppler effects at the shooting position. Figure 6 shows the bullet velocity as a function of distance for bullets with diameters equal to 5.56 mm and 7.62mm. For locations near the shooting distance, the 5.56mm bullet with a diameter of 5.56mm is faster than that the bullet with a diameter of 7.62mm. At the shooting distance of 260 m, the 7.62mm bullet with the diameter of 7.62mm is faster.Figure 6
Bullet velocity measurement setup. (a) Doppler radar equipment. (b) Bullet velocity measurements as a function of the shooting distance.
(a)
(b)Figure7 shows a schematic of the experimental setup for the shock wave detection units. Figures 7(a) and 7(b) show the delta- and H-type of acoustic sensor units with six acoustic sensors in each of them. In addition, Figure 7(c) shows the data processing setup that contains a signal conditioner (PCB-480E09), a data acquisition system (DAQ, NI USB-6366), LabVIEW (NI Signal Express) software, and a recording PC. In addition, a wireless communication device (MAP5020) operating based on Wi-Fi technology and control PC is used to collect and control the data acquisition of generated shock waves at long distances. The sampling rate of the shock wave measurement is 200 kHz, and the data transmission rate is 20 kHz. The shooting distances are 100m and 200m, and the bullet diameters are 5.56mm and 7.62mm.Figure 7
Schematic of the acoustic sensor units. (a) Delta-type sensor unit. (b) H-type sensor unit. (c) Data processing setup.
(a)
(b)
(c)Figure8 shows the test setup of the IPES for shock wave measurement and impact-point identification. Based on the actual shooting and the six acoustic sensors placed in the delta-type and H-type sensor units, shock waves are collected at the shooting distances of 100m and 200m. Figures 9(a) and 9(b) show the shock and blast waves measured at the shooting distances of 100m and 200m, respectively. These results correspond to the shooting experimental results of the 7.62mm bullet with the delta type sensor unit. These outcomes show that the blast waves arrive after the arrival of the shock waves. Figures 9(c) and 9(d) show the shock waves measured by the six acoustic sensors at the shooting distances of 100 m and 200m, respectively, using the same conditions. Herein, S1–S6 represent the six acoustic sensors whose arrangements are shown in Figure 7.Figure 8
Test setup of the impact-point estimation system (IPES). (a) H-type sensor unit. (b) Delta-type sensor unit. (c) Data processing setup.Figure 9
Results of the shooting experiment with a delta-type sensor unit and a 7.62 mm bullet. (a, b) Shock waves and blast waves measured according to the shooting distance. (c, d) Shock waves the showing time differences among the acoustic sensors.
(a)
(b)
(c)
(d)The arrival time of the blast wave can be simply calculated asx2+L2/c. Figure 10 shows the experimental and calculated time delays between the shock and blast waves according to the shooting distance measured from the delta-type sensor unit. Based on this figure, it is confirmed that the time delay of the shock wave increases as the shooting distance increases, and the error range of the experimental values is very small. However, the difference between the experimental and calculated values increases as the shooting distance increases. A maximum difference of 7.0% is obtained for a distance of 200 m. This appears to be caused because the theoretical equation cannot consider the nonlinear characteristics of the flight of the bullet, such as the drag, which occurs in long-distance flights as the shooting distance increases.Figure 10
Experimental and calculated time delays between the shock and blast waves.Use of the point of impact at the instant the bullet penetrated the target in the actual shooting test, allows the calculation of the actual coordinates from the center of the arrayed sensor. To estimate the coordinates of the actual impact point, the time delay of the shock wave arriving at each sensor must be measured. In Figure7, the delta-type sensor unit contains two separate groups of sensors at the left and the right sides, while the H-type sensor unit has two separate groups at the front and rear. Each group comprises three sensors each. Based on this configuration, the arrival times among the sensors can be calculated. Figures 11(a) and 11(b) show the differences between the experimental and calculated TDoA values of the shock wave with bullet diameters of 5.56 mm and 7.62 mm according to the shooting distance in the cases of delta-and H-type sensor units. The TDoA value refers to the difference in the arrival time of the shock wave with respect to the arrival times at other acoustic sensors, and it is based on the signal that arrives first among all the acoustic sensor arrays. The delta-type sensor unit exhibits a deviation difference of ~2.1%, while the H-type sensor unit yields a deviation difference of 14.2%. This indicates that the delta-type sensor unit yields a higher accuracy in the estimation of the arrival time of the shock wave.Figure 11
Differences between the experimental and calculated values of the TDoA between the sensors according to the shooting distance (100 m, 200 m). (a) Delta- and (b) H-type sensor units. The error bars represent the standard deviations from three measurements obtained with two types of bullets.
(a)
(b)


### 2.2. Impact-Point Estimation Algorithm

The impact-point estimation algorithm proposed by Levanon [7] defined the time of arrival (ToA) of the bullet shock wave at each microphone using simultaneous equations and estimated the impact point using least-squares estimation, such as the Gauss–Newton (GN) algorithm.Figure12 shows the shock wave occurrence time point B, microphone position Sk, actual impact point H, and horizontal incidence angle α, when the bullet penetrates the target. If the ToA at which the shock wave that occurs at B and arrives at Sk is defined as tk, tk can be expressed according to Equation (2).
(2)tk=toffset−1Vx−xksinα−ykcosα+x−xkcosα+yksinα2+h−hk2×1C2−1V2where toffset is the common delay time for all the microphones and is expressed as toffset=−y0/V, c is the speed of sound, V is the velocity of the bullet, and y0 is the position of the signal that reaches the microphone. Nonlinear equations are defined in the number of microphones using the arrival time and position of each microphone. The equation listed above contains six variables. Specifically, x and h are impact points, and α is the horizontal incidence angle. Therefore, to calculate the simultaneous equations, six or more equations, i.e., six or more microphone signals are required. In this study, signals were acquired from acoustic sensors using delta- and H-type sensor units.Figure 12
Geometry of target.The method proposed by Levanon [7] uses the ToA method. Compared with the TDoA method, this method produces significantly different operation results when there are ToA measurement errors [8], including the instances at which time synchronization problems occur. Therefore, the impact points of the delta- and H-type sensor units were estimated using the TDoA-based impact-point estimation algorithm proposed by Won and Park [9].Figure13 shows a flow chart of the process used for the impact-point estimation algorithm. The algorithm designates initial values (impact point, speed of sound, and bullet velocity), performs iterative operations using the GN algorithm, and completes the iterative operations and displays the results when the variation of the operation result decreases below the set value. H denotes the partial differential of χ, and R(i) denotes the difference between the measured time difference and the calculated time difference. Δp represents the amount of change in x and y at each iteration.Figure 13
Flow chart of the impact-point estimation algorithm.Even though20 rounds were fired in the shooting experiment, 13 datasets were obtained because seven rounds missed the target. Table 2 lists the shooting conditions and measured coordinates as well as the coordinates calculated using TDoA for the 13 datasets acquired from the shooting experiment. Furthermore, the data for the shooting distance of 200 m included only three datasets. At this time, the x-axis coordinate represents the horizontal coordinate and the z-axis coordinate represents the vertical coordinate. The center point of the target used in the experiment is converted into the coordinates (0, 900).Table 2
The measured and the calculated impact-point coordinates.
Test times
Shooting distance (m)
Bullet diameter (mm)
Measured coordinates (mm)
Calculated Coordinates (mm), Delta-type
Calculated coordinates (mm), H-type
X axis
Z axis
X axis
Z axis
X axis
Z axis
1
100
7.62
60.0
843.0
68.5
844.5
48.7
818.5
2
−68.0
1050.0
−63.5
1048.6
−76.4
1030.3
3
126.0
669.0
128.3
667.0
111.8
648.9
4
−36.0
648.0
−33.0
649.3
−43.6
623.3
5
87.0
1063.0
95.6
1060.0
77.2
1032.4
6
5.56
35.0
965.0
39.6
962.3
No data
No data
7
−47.0
1058.0
−38.4
1048.2
−53.5
1023.0
8
118.0
712.0
112.4
709.4
105.1
689.5
9
−20.0
625.0
−10.2
624.1
−22.6
596.9
10
93.0
922.0
101.4
908.6
77.0
887.6
11
200
7.62
113.0
765.0
116.1
770.7
99.4
743.2
12
231.0
1072.0
228.2
1070.5
208.9
1032.0
13
5.56
55.0
709.0
52.7
712.2
41.9
686.4Figure14 shows the measured impact-point coordinates as well as the calculated impact-point coordinates using the delta- and H-type sensor units. Compared to those estimated for the delta-type sensor unit, the calculated coordinates estimated for the H-type sensor unit are slightly biased to the direction at the bottom–left. In addition, the H-type sensor unit has difficulties the determination of the measured impact point in the upward–downward directions because all six acoustic sensors are present on the same plane, as shown in Figure 15. As the impact-point estimation algorithm performs operations by designating the target center as the initial value, an actual bullet that passes through the bottom of the acoustic sensors can be determined as being positioned above the acoustic sensors. To address this problem, it seems necessary to place additional acoustic sensors or to adjust the positions of the existing acoustic sensors.Figure 14
The measured and calculated impact-point coordinates.Figure 15
Schematic shows a potential problem with the same signal in an H-type sensor unit.Conversely, the delta-type sensor unit is capable of determining the impact point in the upward–downward directions, and thus avoids the problem identified above. The mean absolute errorMAE between the measured coordinates and the calculated coordinates of impact-point position is estimated with Equation (3).
(3)MAE=∑i=1nmi−cinwhere m is the measured coordinate, c is the calculated coordinate and n is the number of data points.Table3 shows the MAE values of the delta- and H-type sensor units for different bullet diameters. In the case of the shooting distance of 100 m, the MAE value of the delta-type is relatively smaller than the MAE value of the H-type. Accordingly, it is considered that the delta-type of impact prediction performance is better. In the case of the H-type sensor unit, the reason why the MAE value of the Z axis is larger than the MAE value of the X axis is attributed to the fact that the problem arises because the sensor is arranged in the X axis as mentioned above. Furthermore, given that the MAE value of the bullet with a diameter of 7.62 mm is relatively smaller than that for the bullet with a diameter of 5.56 mm, it is considered that the impact prediction performance of the bullet with a diameter of 7.62 mm is superior. In the case of the shooting distance of 200 m, the MAE value of the delta-type is also relatively smaller than the MAE value of the H-type. However, given that the MAE value is calculated using only three shooting datasets, the reliability of the MAE value cannot be confirmed.Table 3
Mean absolute error according to bullet type.
Shooting
Bullet
Delta-type
H-type
Distance(m)
Diameter(mm)
X axis MAE (mm)
Z axis MAE (mm)
X axis MAE (mm)
Z axis MAE (mm)
100
5.56
7.40
5.88
9.5
30.00
7.62
5.38
1.84
10.26
23.92
200
5.56
2.30
3.2
13.1
22.60
7.62
2.95
3.6
17.85
30.95


## 2.1. Experimental Measurements

A typical muzzle explosion produces shock waves with a sound pressure level (SPL) ≥ 140 dB. Thus, acoustic sensors used for measuring these shock waves must be able to withstand these sound pressures [6]. Therefore, several fire tests were conducted on various acoustic sensors. Based on these tests, the acoustic sensor that was capable of withstanding shock waves was selected. As shown in Figure 4(a), the structure of the acoustic sensor is composed of a protection grid, an electret condenser microphone (ECM), a preamplifier, housing, and an SMB connector. Figure 4(b) shows the acoustic sensor (HJ06, SISPIA Co., Korea). Each of the acoustic sensors should also have uniform sensitivity and phase characteristics to predict accurately the impact point of the bullet. The sensitivity and phase characteristics of the acoustic sensors were measured in an anechoic chamber according to the standard measurement method. The sensitivity of the acoustic sensors was approximately -46.0 dB in the frequency band of 100 Hz to 10 kHz, and the phase was approximately within ±3° at 1.0 kHz. Figure 5 shows the sensitivity and phase measurement outcomes from the six acoustic sensors, thus indicating the uniform sensitivity and phase characteristics with respect to frequency. The performance specifications of the acoustic sensors are listed in Table 1. The IPES based on the delta- and H-type sensor units was constructed with the use of the six acoustic sensors with uniform sensitivity and phase value.Figure 4
Acoustic sensor. (a) Schematic of structure and constituent components. (b) Photograph of acoustic sensor.
(a)
(b)Figure 5
Sensitivity (solid lines) and phase (dotted lines) measurements obtained from the acoustic sensors.Table 1
Measurement results of the acoustic sensor.
Sensitivity
-46.0 dB relative to 1 V/pa ± 3 dB
Frequency response
± 2 dB, 100 Hz to 10 kHz
Phase response
± 3°, 100 Hz to 3 kHz± 5°, 3 kHz to 5 kHzTo measure the shock wave of a bullet and to test the performance of the IPES, an experiment was performed in an actual shooting test site. Bullets with diameters equal to5.56mm and 7.62mm were used. Prior to the experiment, the bullet velocity was measured at 50 m, 100 m, 150 m, and 300 m from the shooting point with a speed measuring device (Doppler radar SL-528PE, Denmark) based on the Doppler effects at the shooting position. Figure 6 shows the bullet velocity as a function of distance for bullets with diameters equal to 5.56 mm and 7.62mm. For locations near the shooting distance, the 5.56mm bullet with a diameter of 5.56mm is faster than that the bullet with a diameter of 7.62mm. At the shooting distance of 260 m, the 7.62mm bullet with the diameter of 7.62mm is faster.Figure 6
Bullet velocity measurement setup. (a) Doppler radar equipment. (b) Bullet velocity measurements as a function of the shooting distance.
(a)
(b)Figure7 shows a schematic of the experimental setup for the shock wave detection units. Figures 7(a) and 7(b) show the delta- and H-type of acoustic sensor units with six acoustic sensors in each of them. In addition, Figure 7(c) shows the data processing setup that contains a signal conditioner (PCB-480E09), a data acquisition system (DAQ, NI USB-6366), LabVIEW (NI Signal Express) software, and a recording PC. In addition, a wireless communication device (MAP5020) operating based on Wi-Fi technology and control PC is used to collect and control the data acquisition of generated shock waves at long distances. The sampling rate of the shock wave measurement is 200 kHz, and the data transmission rate is 20 kHz. The shooting distances are 100m and 200m, and the bullet diameters are 5.56mm and 7.62mm.Figure 7
Schematic of the acoustic sensor units. (a) Delta-type sensor unit. (b) H-type sensor unit. (c) Data processing setup.
(a)
(b)
(c)Figure8 shows the test setup of the IPES for shock wave measurement and impact-point identification. Based on the actual shooting and the six acoustic sensors placed in the delta-type and H-type sensor units, shock waves are collected at the shooting distances of 100m and 200m. Figures 9(a) and 9(b) show the shock and blast waves measured at the shooting distances of 100m and 200m, respectively. These results correspond to the shooting experimental results of the 7.62mm bullet with the delta type sensor unit. These outcomes show that the blast waves arrive after the arrival of the shock waves. Figures 9(c) and 9(d) show the shock waves measured by the six acoustic sensors at the shooting distances of 100 m and 200m, respectively, using the same conditions. Herein, S1–S6 represent the six acoustic sensors whose arrangements are shown in Figure 7.Figure 8
Test setup of the impact-point estimation system (IPES). (a) H-type sensor unit. (b) Delta-type sensor unit. (c) Data processing setup.Figure 9
Results of the shooting experiment with a delta-type sensor unit and a 7.62 mm bullet. (a, b) Shock waves and blast waves measured according to the shooting distance. (c, d) Shock waves the showing time differences among the acoustic sensors.
(a)
(b)
(c)
(d)The arrival time of the blast wave can be simply calculated asx2+L2/c. Figure 10 shows the experimental and calculated time delays between the shock and blast waves according to the shooting distance measured from the delta-type sensor unit. Based on this figure, it is confirmed that the time delay of the shock wave increases as the shooting distance increases, and the error range of the experimental values is very small. However, the difference between the experimental and calculated values increases as the shooting distance increases. A maximum difference of 7.0% is obtained for a distance of 200 m. This appears to be caused because the theoretical equation cannot consider the nonlinear characteristics of the flight of the bullet, such as the drag, which occurs in long-distance flights as the shooting distance increases.Figure 10
Experimental and calculated time delays between the shock and blast waves.Use of the point of impact at the instant the bullet penetrated the target in the actual shooting test, allows the calculation of the actual coordinates from the center of the arrayed sensor. To estimate the coordinates of the actual impact point, the time delay of the shock wave arriving at each sensor must be measured. In Figure7, the delta-type sensor unit contains two separate groups of sensors at the left and the right sides, while the H-type sensor unit has two separate groups at the front and rear. Each group comprises three sensors each. Based on this configuration, the arrival times among the sensors can be calculated. Figures 11(a) and 11(b) show the differences between the experimental and calculated TDoA values of the shock wave with bullet diameters of 5.56 mm and 7.62 mm according to the shooting distance in the cases of delta-and H-type sensor units. The TDoA value refers to the difference in the arrival time of the shock wave with respect to the arrival times at other acoustic sensors, and it is based on the signal that arrives first among all the acoustic sensor arrays. The delta-type sensor unit exhibits a deviation difference of ~2.1%, while the H-type sensor unit yields a deviation difference of 14.2%. This indicates that the delta-type sensor unit yields a higher accuracy in the estimation of the arrival time of the shock wave.Figure 11
Differences between the experimental and calculated values of the TDoA between the sensors according to the shooting distance (100 m, 200 m). (a) Delta- and (b) H-type sensor units. The error bars represent the standard deviations from three measurements obtained with two types of bullets.
(a)
(b)


## 2.2. Impact-Point Estimation Algorithm

The impact-point estimation algorithm proposed by Levanon [7] defined the time of arrival (ToA) of the bullet shock wave at each microphone using simultaneous equations and estimated the impact point using least-squares estimation, such as the Gauss–Newton (GN) algorithm.Figure12 shows the shock wave occurrence time point B, microphone position Sk, actual impact point H, and horizontal incidence angle α, when the bullet penetrates the target. If the ToA at which the shock wave that occurs at B and arrives at Sk is defined as tk, tk can be expressed according to Equation (2).
(2)tk=toffset−1Vx−xksinα−ykcosα+x−xkcosα+yksinα2+h−hk2×1C2−1V2where toffset is the common delay time for all the microphones and is expressed as toffset=−y0/V, c is the speed of sound, V is the velocity of the bullet, and y0 is the position of the signal that reaches the microphone. Nonlinear equations are defined in the number of microphones using the arrival time and position of each microphone. The equation listed above contains six variables. Specifically, x and h are impact points, and α is the horizontal incidence angle. Therefore, to calculate the simultaneous equations, six or more equations, i.e., six or more microphone signals are required. In this study, signals were acquired from acoustic sensors using delta- and H-type sensor units.Figure 12
Geometry of target.The method proposed by Levanon [7] uses the ToA method. Compared with the TDoA method, this method produces significantly different operation results when there are ToA measurement errors [8], including the instances at which time synchronization problems occur. Therefore, the impact points of the delta- and H-type sensor units were estimated using the TDoA-based impact-point estimation algorithm proposed by Won and Park [9].Figure13 shows a flow chart of the process used for the impact-point estimation algorithm. The algorithm designates initial values (impact point, speed of sound, and bullet velocity), performs iterative operations using the GN algorithm, and completes the iterative operations and displays the results when the variation of the operation result decreases below the set value. H denotes the partial differential of χ, and R(i) denotes the difference between the measured time difference and the calculated time difference. Δp represents the amount of change in x and y at each iteration.Figure 13
Flow chart of the impact-point estimation algorithm.Even though20 rounds were fired in the shooting experiment, 13 datasets were obtained because seven rounds missed the target. Table 2 lists the shooting conditions and measured coordinates as well as the coordinates calculated using TDoA for the 13 datasets acquired from the shooting experiment. Furthermore, the data for the shooting distance of 200 m included only three datasets. At this time, the x-axis coordinate represents the horizontal coordinate and the z-axis coordinate represents the vertical coordinate. The center point of the target used in the experiment is converted into the coordinates (0, 900).Table 2
The measured and the calculated impact-point coordinates.
Test times
Shooting distance (m)
Bullet diameter (mm)
Measured coordinates (mm)
Calculated Coordinates (mm), Delta-type
Calculated coordinates (mm), H-type
X axis
Z axis
X axis
Z axis
X axis
Z axis
1
100
7.62
60.0
843.0
68.5
844.5
48.7
818.5
2
−68.0
1050.0
−63.5
1048.6
−76.4
1030.3
3
126.0
669.0
128.3
667.0
111.8
648.9
4
−36.0
648.0
−33.0
649.3
−43.6
623.3
5
87.0
1063.0
95.6
1060.0
77.2
1032.4
6
5.56
35.0
965.0
39.6
962.3
No data
No data
7
−47.0
1058.0
−38.4
1048.2
−53.5
1023.0
8
118.0
712.0
112.4
709.4
105.1
689.5
9
−20.0
625.0
−10.2
624.1
−22.6
596.9
10
93.0
922.0
101.4
908.6
77.0
887.6
11
200
7.62
113.0
765.0
116.1
770.7
99.4
743.2
12
231.0
1072.0
228.2
1070.5
208.9
1032.0
13
5.56
55.0
709.0
52.7
712.2
41.9
686.4Figure14 shows the measured impact-point coordinates as well as the calculated impact-point coordinates using the delta- and H-type sensor units. Compared to those estimated for the delta-type sensor unit, the calculated coordinates estimated for the H-type sensor unit are slightly biased to the direction at the bottom–left. In addition, the H-type sensor unit has difficulties the determination of the measured impact point in the upward–downward directions because all six acoustic sensors are present on the same plane, as shown in Figure 15. As the impact-point estimation algorithm performs operations by designating the target center as the initial value, an actual bullet that passes through the bottom of the acoustic sensors can be determined as being positioned above the acoustic sensors. To address this problem, it seems necessary to place additional acoustic sensors or to adjust the positions of the existing acoustic sensors.Figure 14
The measured and calculated impact-point coordinates.Figure 15
Schematic shows a potential problem with the same signal in an H-type sensor unit.Conversely, the delta-type sensor unit is capable of determining the impact point in the upward–downward directions, and thus avoids the problem identified above. The mean absolute errorMAE between the measured coordinates and the calculated coordinates of impact-point position is estimated with Equation (3).
(3)MAE=∑i=1nmi−cinwhere m is the measured coordinate, c is the calculated coordinate and n is the number of data points.Table3 shows the MAE values of the delta- and H-type sensor units for different bullet diameters. In the case of the shooting distance of 100 m, the MAE value of the delta-type is relatively smaller than the MAE value of the H-type. Accordingly, it is considered that the delta-type of impact prediction performance is better. In the case of the H-type sensor unit, the reason why the MAE value of the Z axis is larger than the MAE value of the X axis is attributed to the fact that the problem arises because the sensor is arranged in the X axis as mentioned above. Furthermore, given that the MAE value of the bullet with a diameter of 7.62 mm is relatively smaller than that for the bullet with a diameter of 5.56 mm, it is considered that the impact prediction performance of the bullet with a diameter of 7.62 mm is superior. In the case of the shooting distance of 200 m, the MAE value of the delta-type is also relatively smaller than the MAE value of the H-type. However, given that the MAE value is calculated using only three shooting datasets, the reliability of the MAE value cannot be confirmed.Table 3
Mean absolute error according to bullet type.
Shooting
Bullet
Delta-type
H-type
Distance(m)
Diameter(mm)
X axis MAE (mm)
Z axis MAE (mm)
X axis MAE (mm)
Z axis MAE (mm)
100
5.56
7.40
5.88
9.5
30.00
7.62
5.38
1.84
10.26
23.92
200
5.56
2.30
3.2
13.1
22.60
7.62
2.95
3.6
17.85
30.95


## 3. Conclusions

This study presents the first performance evaluation outcomes on the accuracy of the IPES according to two types of acoustic sensor arrangement methods. First, the acoustic sensor suitable for measuring shock waves was selected, and the various acoustic sensors were then tested. Subsequently, acoustic sensors with uniform sensitivity and phase characteristics were used in the shooting experiment, and it was confirmed that they had almost identical response characteristics for shock wave measurements. An IPES was constructed with the use of delta-type and H-type sensor units with the selected acoustic sensors, and shock waves were collected by performing shooting experiments at the shooting distances of100m and 200m, and with two types of bullets with the diameters of 5.56mm and 7.62mm. By comparing the experimental and calculated TDoA values of the shock wave for delta-type and H-type sensor units, the TDoA values of the delta type sensor unit exhibited a low difference deviation of 2.1%. The coordinates of the impact point were obtained from actual shooting tests and the impact-point estimation algorithm based on TDoA datasets for the delta- and H-type sensor units. Because the MAE values of the delta-type sensor unit were relatively smaller than those for the H-type sensor unit, it was considered that the performance of the delta-type impact-point estimation was better than that of the H-type unit. It is expected that this study will contribute to the improvement of the accuracy of the IPES based on the acoustic sensor arrangement method.



---

*Source: 1758765-2019-11-11.xml*