# Pyroshock Acceleration Field Reconstruction in Temporal and Spectral Domains Based on Laser Shock Scanning and Iterative Decomposition and Synthesis Considering Stop Band Effects

Authors: Yong-Woon Kim; Jae-Kyeong Jang; Jung-Ryul Lee
Journal: Shock and Vibration
 (2017)

---

## Abstract

Pyrotechnic devices are used to separate substructures from main structures. Pyroshock can cause failure in electronic components that are sensitive to high frequency shock. Most of the existing methods to analyze pyroshock have limitations for high frequency simulations and are only available for simulation of point explosive-induced pyroshock. To solve the problem of existing methods, we developed a laser shock-based pyroshock reconstruction algorithm covering high frequency range that can predict linear explosive-induced pyroshock, as well as point explosive-induced ones. The developed algorithm reconstructs pyroshock from laser shock test in both temporal and spectral domains using an iterative signal decomposition and synthesis method. In the signal decomposition and synthesis process, unremoved signals in the stopbands occurred and were compensated by iteration to improve the results. At the end of this paper, various types of pyroshock were processed through the proposed method. Pyroshock wave propagation images and shock response spectrum images were presented as a result. To verify the algorithm, we compared the obtained result with a real pyroshock. The time domain signal was reconstructed with an averaged peak to peak acceleration difference of 20.21%, and the shock response spectrum was reconstructed with an average mean acceleration difference of 25.86%.


---

## Body

## 1. Introduction

Pyrotechnic devices are used to separate substructures from main structures in such situations as vehicle launches, rocket-stage separations, missile launches, and pilot ejections [1–3]. Pyrotechnics have advantageous features, such as high energy per unit weight and volume [4]. These features have ensured that pyrotechnic devices are widely used in aerospace structures. On the other hand, pyrotechnics have a critical safety problem. During operation, high acceleration and high frequency structural shock, called pyroshock, always occur. Pyroshock rarely influences structural members, but its high frequency characteristic can easily cause failure in electronic components [5]. In a NASA technical memorandum, Laurence described a survey compiled over a 23-year period, which included 84 serious component or system failures. Out of the failures, 12 occurred in-flight with fully developed and qualified hardware [6]. Therefore, it is necessary to evaluate the pyroshock transmitted to the structure and to predict the damage of the electronic equipment.Previous research implemented pyroshock using real pyro devices that downscaled structures. This technique had the advantage of being capable of simulating pyroshock in the most reasonable fashion. However, since the actual pyro device was used, the structures were damaged and could not be used again. Therefore, it was necessary to manufacture the respective structures separately for all the experiments, so that the economic burden and time were inevitable [7–9]. Because of problems with cost and safety, pyroshock imitation techniques using mechanical impact instead of pyro devices were studied. This technique simulates pyroshock by causing mechanical collisions using pendulum motion. These tests yielded excellent simulation results in the mid-field and far-field, but simulation results were not good in the near-field [10–12].In addition, previous studies that simulated pyroshock using numerical analysis were reviewed. Among them, recent research into Hydro-code-based pyroshock simulation has been reported as suitable for simulating highly dynamic phenomena such as pyroshock [13–15]. However, most of these numerical techniques required accurate knowledge of all parameters of the structure to obtain reliable simulation results, and computational time and cost limitations emerged when applied to complex structures [16, 17].Recently, laser shock has been used to reconstruct pyroshock in [18]. The pyroshock reconstruction method, using laser shock, had the advantage that it could confirm propagation of pyroshock in real structure nondestructively. Since shock generating and measuring processes are distinct, there are characteristic differences between laser shock and pyroshock, even though both are elastic wave in solid medium. In [18], the reconstructed results has used an algorithm that only adjusts the peak value to compensate for this characteristic difference. This resulted in a below average reconstruction response. In [19], Chong has reconstructed pyroshock by statistically predicting how the amplitude of the pyroshock changes with distance from the shock source and provided the predicted information to the laser shock. However, since the pyroshock not only is dependent on the distance from the explosion source but also depends on the geometry of the structure, the algorithm developed in [19] was able to reconstruct only on the simple specimens. Therefore, to reconstruct pyroshock using the laser shock, an innovative algorithm that can appropriately compensate for the characteristic difference between the two shocks is needed.In this paper, we present a novel algorithm to reconstruct the waveform and spectrum response of pyroshock and laser shock for both linear and point-wise pyrotechnical devices to overcome the limitations of existing laser shock-based methods. The proposed algorithm compensates for the difference in characteristics of laser shock and pyroshock by using the iterative signal decomposition and synthesis method, which achieved the gain needed to compensate for the difference in characteristics between laser shock and pyroshock based on shock response spectrum (SRS). The gain was applied to laser shock signal decomposition via Butterworth filter in each frequency band. In the laser shock decomposition process, due to the inevitable characteristics of the band-pass filter, the decomposed signal in the passband included the unwanted signals that in-flowed from stopbands. This is an interference problem that was compensated by iteratively removing the unwanted signal to improve the reconstruction result.The reconstruction results are presented using the proposed pyroshock wave propagation imaging (PWPI) and shock response spectrum imaging (SRSI) algorithms. We can validate propagation characteristics such as propagation path and reflection waves of pyroshocks based on PWPI results in the time-space domain and determine how much the structure is affected by the pyroshock for each frequency using the SRSI results. To verify the proposed algorithm, the comparison results with the real pyroshock are presented in terms of similarity between the real pyroshock and reconstructed pyroshock.


## 2. Simultaneous Temporal and Spectral Pyroshock Reconstruction Algorithm

Figure1 shows the flow chart of the pyroshock reconstruction algorithm developed in this paper. The algorithm is divided into 5 steps. In the first step, real pyroshock signals are measured using LDV sensors. In the second step, laser shock signals generated by Q-switch laser are acquired through a PZT sensor. In the third step, the shock response spectrum (SRS) of the pyroshock is reconstructed by compensating the characteristic difference between the laser shock and pyroshock. In the fourth step, the signal decomposition and synthesis method is used to reconstruct the time domain signal of pyroshock by giving the information of the reconstructed SRS to laser shock data. In the fifth step, PWPI and SRSI are generated using the reconstructed pyroshock data. The details of each step are explained in the following sections.Figure 1
Flow chart of pyroshock reconstruction algorithm based on laser shock.


### 2.1. Real Pyroshock Measurement at Multipoints

The pyroshock measurement is performed through a shock measurement system based on multiple laser Doppler vibrometers (LDVs, Polytech Inc.) as shown in Figure2. First, the target structure and the pyro devices are installed and the area to be analyzed is specified. The point at which you want to acquire the pyroshock is determined and the coordinates of that point are indicated by (xs,ys). Finally, we set up the DAQ configuration through the GUI platform and wait for the pyro device to explode. Then, when the pyro device is exploded, pyroshock signals are collected from the LDV sensor according to the DAQ condition set through the GUI platform. The pyroshock signals collected from the sensors are stored in the PC as digital voltage form. The pyroshock signal is stored in the digital voltage form and is denoted by VD(xs,ys,t). It is then calibrated to the pyroshock signal P(xs,ys,t) of the acceleration unit through a conversion process using (1).(1)Pxs,ys,t=dVDxs,ys,t×CF×10AdB/20×CFmrdt,where VD(xs,ys,t) is pyroshock signal of digital voltage form [V]. P(xs,ys,t) is pyroshock signal of acceleration unit [g]. CF is calibration factor of LDV sensor. AdB is attenuation factor of attenuator. CFmr is measurement unit of LDV sensor [m/s/V].Figure 2
Pyroshock measurement process using a shock measurement system based on multiple laser Doppler vibrometers (LDVs).


### 2.2. Laser Shock Scanning Measurement

Through the pyroshock measurement process, we can observe the characteristics of the pyroshock. However, since the number of LDV sensors is limited, pyroshock can be detected only at a few limited points. To confirm the propagation of pyroshock for the whole structure, we reconstruct pyroshock using laser shock. First, the same structure used in the pyroshock measurement is prepared to obtain the laser shock signal. Then, PZT sensors (R-Cast M204A, Fuji ceramics Inc.) are installed in place of the pyrotechnical device because the reciprocal setup of the laser shock scanning system (Mobile PZT-UPI, X-NDT Inc.) designates the sensor locations shock source locations. Once the structure and sensors are installed, the laser scanning process proceeds as shown in Figure3. The GUI platform, shown in Figure 3(a), sets the DAQ conditions and transmits commands to the Q-switch laser and the laser mirror scanner. The laser mirror scanner moves the laser along the scan path. The laser beam moves through the laser mirror scanner and generates the laser shock at the position where it will stay, as shown in Figure 3(b). The laser shock is generated as shown in Figure 3(c). It is collected through a PZT sensor and subsequently stored in a PC. When the laser scanning process is completed, four-dimensional laser shock data are obtained as shown in Figure 3(d). The obtained laser shock data are denoted by L(x,y,t). If the pyroshock is not induced by a single point source but is induced by a combined explosion such as linear explosive induced pyroshock, the wave superposition principle is used under the assumption that the structural response to the pyroshock is linear. Many PZT sensors are placed at the site of the combined explosion, and laser shock data are acquired from each sensor. The laser scanning data obtained from each sensor are denoted as Lscan(x,y,t,n). To simulate pyroshock induced by combined explosion sources, Lscan(x,y,t,n) is phase arrayed and superimposed according to the explosion mechanism. The superimposed result is denoted as L(x,y,t). In L(x,y,t), the data at the same point as the point at which the pyroshock was measured is designated L(xs,ys,t).Figure 3
Laser shock measurement process scheme; (a) PC with GUI platform. (b) Thermoelastic mechanism. (c) Laser shock time domain signal. (d) Laser scanning data.


### 2.3. Shock Response Spectrum Reconstruction

To reconstruct the pyroshock using laser shock, it is necessary to compensate the characteristic differences between laser shock and pyroshock. The compensation procedure is based on the SRS. The SRS is the most widely used function to quantify pyroshock and is calculated based on the time domain signal. It applies the time domain signal as a base excitation to an array of single degree of freedom (SDOF) systems [21], where the independent variable is the natural frequency. Therefore, the calculation is performed on multiple independent SDOF systems. Any arbitrary set of natural frequencies can be used for the SRS calculation. However, the typical method is based on a proportional bandwidth such as 1/12 octave. Therefore, in this paper, the SRS is obtained at 1/12 octave intervals. The detailed procedure for calculating the SRS is described in [22]. Using the obtained P(xs,ys,t) and L(xs,ys,t), we calculate the SRS and define the results as SRSP(xs,ys,fn) and SRSL(xs,ys,fn), respectively. The ratio of the two SRSs is defined as GSRSxs,ys,fn.(2)GSRSxs,ys,fn=SRSPxs,ys,fnSRSLxs,ys,fn.In the system,GSRS(xs,ys,fn), the value that best compensates the difference between the laser shock and the pyroshock on average, for the whole set of training points, is designated GSRS,O(fn). By multiplying the system GSRS,O(fn) by SRSL(x,y,fn), the SRS reconstruction process is completed. The value obtained through the SRS reconstruction process is defined as SRSref(x,y,fn). (3)SRSrefx,y,fn=SRSLx,y,fn∗GSRS,Ofn.The similarity between any two SRSs is represented by MAD. The MAD can be obtained from [20](4)MAD%=1z∑fn=f1fzSRS1fn-SRS2fnSRS1fn∗100%.


### 2.4. Time Domain Signal Reconstruction

The time domain signal of the pyroshock is reconstructed by substitutingSRSref(x,y,fn) to the time domain signal of the laser shock. The reconstruction of time domain signal is based on signal decomposition and synthesis. Signal decomposition is conducted via a 2nd order Butterworth filter with a zero-phase filter, which is a type of signal processing filter designed to have a very flat frequency response in the passband. Figure 4 shows the results of signal decomposition and synthesis for arbitrary shock using Butterworth filter. The black signal represents the original signal and the blue signal represents the reconstructed signal after signal decomposition and synthesis using a Butterworth filter, which was suitable for the signal decomposition and synthesis since the original signal and the reconstructed signal were almost identical. Figure 5(a) shows a test algorithm that uses the signal decomposition and synthesis method. The gain needed to compensate for the difference between pyroshock and laser shock was obtained in the SRS domain. Then the gain was applied to a laser shock signal decomposed by Butterworth filter as shown in Figure 5(a). Figure 5(b) shows the result of the test algorithm. The blue signal shows reconstructed pyroshock and the black signal shows the real pyroshock. As shown in Figure 5(b), the reconstruction was not performed well because of the filter’s problem. Since the Butterworth filter is not ideal filter, decomposed signal contains unwanted signals. For example, Figure 5(c) shows the laser shock signal decomposed between 100 Hz and 106 Hz, via the Butterworth filter. It could be seen that unremoved high frequency components in the stopbands were mixed. These unremoved signals in the stopbands cause interference in the passband signals and ruin reconstruction results. To solve this stopband effect in band-pass filtering, the iteration method was applied as shown in Figure 6(a). As a result, we could obtain significantly improved results as shown in Figure 6(b). Finally, we designed the reconstruction algorithm of time domain signal as shown in Figure 7 by applying SRSrefx,y,fn obtained in Section 2.3 to the algorithm of Figure 6(a). First, we calculate the ratio of SRSL(x,y,fn) and SRSref(x,y,fn) in each frequency domain as shown in Figure 7 and store it as a gain value. Next, the laser shock signal is decomposed by each frequency range using a Butterworth filter. The decomposed laser shock signal is denoted by LBW(x,y,t,fn) and the previously stored gain result is applied. After applying the gain, Ptrainedx,y,t was obtained by synthesizing LBW(x,y,t,fn) again as (5). Finally, the interference problem of the unremoved signal in the laser shock decomposition process is compensated for by iteration. The compensation process is based on (5); the gain is updated by feedback of the stop band effect as shown in Figure 7. This process is repeated until the MAD value between SRSref(x,y,fn) and SRStrained(x,y,fn) is less than 20%. If SRStrained(x,y,fn) is sufficiently similar to SRSref(x,y,fn), Ptrainedx,y,t and SRStrained(x,y,fn) are stored as PRx,y,t and SRSR(x,y,fn).(5)PRx,y,t=∑fn=100Fs/2LBWx,y,t,fnSRSrefx,y,fnSRSLx,y,fn∏SRSrefx,y,fnSRStrainedx,y,fn.Figure 4
Signal decomposition and synthesis using Butterworth band-pass filter.Figure 5
(a) Test algorithm. (b) Reconstruction results of the test algorithm. (c) The problems of test algorithm.
(a)
(b)
(c)Figure 6
(a) Test algorithm 2. (b) Reconstruction results of the test algorithm 2.
(a)
(b)Figure 7
Reconstruction algorithm of pyroshock time domain signal.


### 2.5. Pyroshock Wave Propagation Image and Shock Response Spectrum Image

When the amplitude ofPRx,y,t, obtained from the previous process is plotted in color scale over time, the PWPI is obtained as shown in Figure 8(a). In addition, when the amplitude of SRSR(x,y,fn) is imaged in color scale as frequency increases, the SRSI is obtained as shown in Figure 8(b).Figure 8
Pyroshock propagation evaluation; (a) PWPI; (b) SRSI.
(a)
(b)


## 2.1. Real Pyroshock Measurement at Multipoints

The pyroshock measurement is performed through a shock measurement system based on multiple laser Doppler vibrometers (LDVs, Polytech Inc.) as shown in Figure2. First, the target structure and the pyro devices are installed and the area to be analyzed is specified. The point at which you want to acquire the pyroshock is determined and the coordinates of that point are indicated by (xs,ys). Finally, we set up the DAQ configuration through the GUI platform and wait for the pyro device to explode. Then, when the pyro device is exploded, pyroshock signals are collected from the LDV sensor according to the DAQ condition set through the GUI platform. The pyroshock signals collected from the sensors are stored in the PC as digital voltage form. The pyroshock signal is stored in the digital voltage form and is denoted by VD(xs,ys,t). It is then calibrated to the pyroshock signal P(xs,ys,t) of the acceleration unit through a conversion process using (1).(1)Pxs,ys,t=dVDxs,ys,t×CF×10AdB/20×CFmrdt,where VD(xs,ys,t) is pyroshock signal of digital voltage form [V]. P(xs,ys,t) is pyroshock signal of acceleration unit [g]. CF is calibration factor of LDV sensor. AdB is attenuation factor of attenuator. CFmr is measurement unit of LDV sensor [m/s/V].Figure 2
Pyroshock measurement process using a shock measurement system based on multiple laser Doppler vibrometers (LDVs).


## 2.2. Laser Shock Scanning Measurement

Through the pyroshock measurement process, we can observe the characteristics of the pyroshock. However, since the number of LDV sensors is limited, pyroshock can be detected only at a few limited points. To confirm the propagation of pyroshock for the whole structure, we reconstruct pyroshock using laser shock. First, the same structure used in the pyroshock measurement is prepared to obtain the laser shock signal. Then, PZT sensors (R-Cast M204A, Fuji ceramics Inc.) are installed in place of the pyrotechnical device because the reciprocal setup of the laser shock scanning system (Mobile PZT-UPI, X-NDT Inc.) designates the sensor locations shock source locations. Once the structure and sensors are installed, the laser scanning process proceeds as shown in Figure3. The GUI platform, shown in Figure 3(a), sets the DAQ conditions and transmits commands to the Q-switch laser and the laser mirror scanner. The laser mirror scanner moves the laser along the scan path. The laser beam moves through the laser mirror scanner and generates the laser shock at the position where it will stay, as shown in Figure 3(b). The laser shock is generated as shown in Figure 3(c). It is collected through a PZT sensor and subsequently stored in a PC. When the laser scanning process is completed, four-dimensional laser shock data are obtained as shown in Figure 3(d). The obtained laser shock data are denoted by L(x,y,t). If the pyroshock is not induced by a single point source but is induced by a combined explosion such as linear explosive induced pyroshock, the wave superposition principle is used under the assumption that the structural response to the pyroshock is linear. Many PZT sensors are placed at the site of the combined explosion, and laser shock data are acquired from each sensor. The laser scanning data obtained from each sensor are denoted as Lscan(x,y,t,n). To simulate pyroshock induced by combined explosion sources, Lscan(x,y,t,n) is phase arrayed and superimposed according to the explosion mechanism. The superimposed result is denoted as L(x,y,t). In L(x,y,t), the data at the same point as the point at which the pyroshock was measured is designated L(xs,ys,t).Figure 3
Laser shock measurement process scheme; (a) PC with GUI platform. (b) Thermoelastic mechanism. (c) Laser shock time domain signal. (d) Laser scanning data.


## 2.3. Shock Response Spectrum Reconstruction

To reconstruct the pyroshock using laser shock, it is necessary to compensate the characteristic differences between laser shock and pyroshock. The compensation procedure is based on the SRS. The SRS is the most widely used function to quantify pyroshock and is calculated based on the time domain signal. It applies the time domain signal as a base excitation to an array of single degree of freedom (SDOF) systems [21], where the independent variable is the natural frequency. Therefore, the calculation is performed on multiple independent SDOF systems. Any arbitrary set of natural frequencies can be used for the SRS calculation. However, the typical method is based on a proportional bandwidth such as 1/12 octave. Therefore, in this paper, the SRS is obtained at 1/12 octave intervals. The detailed procedure for calculating the SRS is described in [22]. Using the obtained P(xs,ys,t) and L(xs,ys,t), we calculate the SRS and define the results as SRSP(xs,ys,fn) and SRSL(xs,ys,fn), respectively. The ratio of the two SRSs is defined as GSRSxs,ys,fn.(2)GSRSxs,ys,fn=SRSPxs,ys,fnSRSLxs,ys,fn.In the system,GSRS(xs,ys,fn), the value that best compensates the difference between the laser shock and the pyroshock on average, for the whole set of training points, is designated GSRS,O(fn). By multiplying the system GSRS,O(fn) by SRSL(x,y,fn), the SRS reconstruction process is completed. The value obtained through the SRS reconstruction process is defined as SRSref(x,y,fn). (3)SRSrefx,y,fn=SRSLx,y,fn∗GSRS,Ofn.The similarity between any two SRSs is represented by MAD. The MAD can be obtained from [20](4)MAD%=1z∑fn=f1fzSRS1fn-SRS2fnSRS1fn∗100%.


## 2.4. Time Domain Signal Reconstruction

The time domain signal of the pyroshock is reconstructed by substitutingSRSref(x,y,fn) to the time domain signal of the laser shock. The reconstruction of time domain signal is based on signal decomposition and synthesis. Signal decomposition is conducted via a 2nd order Butterworth filter with a zero-phase filter, which is a type of signal processing filter designed to have a very flat frequency response in the passband. Figure 4 shows the results of signal decomposition and synthesis for arbitrary shock using Butterworth filter. The black signal represents the original signal and the blue signal represents the reconstructed signal after signal decomposition and synthesis using a Butterworth filter, which was suitable for the signal decomposition and synthesis since the original signal and the reconstructed signal were almost identical. Figure 5(a) shows a test algorithm that uses the signal decomposition and synthesis method. The gain needed to compensate for the difference between pyroshock and laser shock was obtained in the SRS domain. Then the gain was applied to a laser shock signal decomposed by Butterworth filter as shown in Figure 5(a). Figure 5(b) shows the result of the test algorithm. The blue signal shows reconstructed pyroshock and the black signal shows the real pyroshock. As shown in Figure 5(b), the reconstruction was not performed well because of the filter’s problem. Since the Butterworth filter is not ideal filter, decomposed signal contains unwanted signals. For example, Figure 5(c) shows the laser shock signal decomposed between 100 Hz and 106 Hz, via the Butterworth filter. It could be seen that unremoved high frequency components in the stopbands were mixed. These unremoved signals in the stopbands cause interference in the passband signals and ruin reconstruction results. To solve this stopband effect in band-pass filtering, the iteration method was applied as shown in Figure 6(a). As a result, we could obtain significantly improved results as shown in Figure 6(b). Finally, we designed the reconstruction algorithm of time domain signal as shown in Figure 7 by applying SRSrefx,y,fn obtained in Section 2.3 to the algorithm of Figure 6(a). First, we calculate the ratio of SRSL(x,y,fn) and SRSref(x,y,fn) in each frequency domain as shown in Figure 7 and store it as a gain value. Next, the laser shock signal is decomposed by each frequency range using a Butterworth filter. The decomposed laser shock signal is denoted by LBW(x,y,t,fn) and the previously stored gain result is applied. After applying the gain, Ptrainedx,y,t was obtained by synthesizing LBW(x,y,t,fn) again as (5). Finally, the interference problem of the unremoved signal in the laser shock decomposition process is compensated for by iteration. The compensation process is based on (5); the gain is updated by feedback of the stop band effect as shown in Figure 7. This process is repeated until the MAD value between SRSref(x,y,fn) and SRStrained(x,y,fn) is less than 20%. If SRStrained(x,y,fn) is sufficiently similar to SRSref(x,y,fn), Ptrainedx,y,t and SRStrained(x,y,fn) are stored as PRx,y,t and SRSR(x,y,fn).(5)PRx,y,t=∑fn=100Fs/2LBWx,y,t,fnSRSrefx,y,fnSRSLx,y,fn∏SRSrefx,y,fnSRStrainedx,y,fn.Figure 4
Signal decomposition and synthesis using Butterworth band-pass filter.Figure 5
(a) Test algorithm. (b) Reconstruction results of the test algorithm. (c) The problems of test algorithm.
(a)
(b)
(c)Figure 6
(a) Test algorithm 2. (b) Reconstruction results of the test algorithm 2.
(a)
(b)Figure 7
Reconstruction algorithm of pyroshock time domain signal.


## 2.5. Pyroshock Wave Propagation Image and Shock Response Spectrum Image

When the amplitude ofPRx,y,t, obtained from the previous process is plotted in color scale over time, the PWPI is obtained as shown in Figure 8(a). In addition, when the amplitude of SRSR(x,y,fn) is imaged in color scale as frequency increases, the SRSI is obtained as shown in Figure 8(b).Figure 8
Pyroshock propagation evaluation; (a) PWPI; (b) SRSI.
(a)
(b)


## 3. Pyroshock Reconstruction and Verification

### 3.1. Reconstruction of Linear Explosive Induced Pyroshock in CFRP Specimen

In this section, pyroshock data in ref [16] was secured to verify and compare the performance of the developed algorithm (Table 1). The pyroshock data were obtained by exploding the linear pyrotechnical device in the quasi-isotropic CFRP sandwich panel. The laser shock was obtained for the same CFRP specimen, and the PWPI and SRSI were obtained using the developed algorithm. The verification of the algorithm was conducted by comparing the real pyroshock with the reconstructed pyroshock. We describe each step of Figure 1 in detail to help understand the reconstruction process.Table 1
Laser shock measurement condition for reconstruction of linear explosive induced pyroshock.
Shock source
F
s (MHz)
Filter (kHz)
Number of samples
PRF (Hz)
Scan area (mm)
Laser shock
1.25
0.1~100
600
1,000
Width = 640Height = 640Interval = 2


#### 3.1.1. Linear Explosive Induced Pyroshock Measurement in CFRP Specimen

The pyroshock was generated by exploding a 230-mm linear pyrotechnical device on a CFRP specimen as shown in Figure9(a). Three LDV sensors were used to measure pyroshock at each measurement point indicated in Figure 9(b). The measured pyroshock was converted into a time domain signal in acceleration as shown in Figure 9(c) through the conversion process and stored in the PC.Figure 9
Linear explosive induced pyroshock measurement on the CFRP specimen. (a) Experimental setup. (b) Sensing point on the test specimen. (c) Time domain signal of linear explosive induced pyroshock.
(a)
(b)
(c)


#### 3.1.2. Laser Shock Scanning Measurement in CFRP Specimen

The same specimen used in pyroshock measurement was installed as shown in Figure10(a), and the laser scanning process was performed along the scanning path as shown in Figure 10(b). When the laser scanning process was finished, the laser scanning data of the 4-dimensional shape was obtained as shown in Figure 10(c). To simulate 230 mm-long linear pyrotechnical device, eleven PZT sensors were installed at intervals of 23 mm, as shown in Figure 11(a). The pyrotechnical device used in this paper started to explode from the right and had a principal propagation direction of 19.6 degree [18], so we phase arrayed the laser scanning data and synthesized it as shown in Figure 11(b). The phased-array laser scanning data were used as Lx,y,t. Figure 11(c) showed the extraction of laser shock data at the same point as the pyroshock measurement point.Figure 10
(a) Laser shock measurement setup. (b) Laser scanning process and laser shock time domain signal at an arbitrary point. (c) Four-dimensional laser scanning data structure.
(a)
(b)
(c)Figure 11
(a) Snapshots of laser scanning data at 65μs obtained by each sensor. (b) Snapshot of synthesized laser scanning data at 100 μs. (c) Training point laser shock, L(xs,ys,t).


#### 3.1.3. Shock Response Spectrum Reconstruction of Linear Explosive Induced Pyroshock

After completion of the laser shock and pyroshock data acquisition, the SRS reconstruction process was performed as described in Section2.3. We used the measurement points 1 and 2, as shown in Figure 12(a), as the training points to reconstruct the SRS. Measurement point 3 was used as the verification point to verify the reconstruction results. In other words, the real pyroshock signals obtained in the only measurement points 1 and 2 were used to train the algorithm and the algorithm provided the reconstruction results and then the results were compared with the real pyroshock at the measurement point 3. The reconstruction results at the training points are shown in Figure 12(b). The black signal represents the SRS of the real pyroshock and the blue signal represents the SRS obtained through the reconstruction process. At the point 1, the SRS was trained with a MAD of 4.34%. At the point 2, the SRS was trained with a MAD of 20.80%.Figure 12
SRS reconstruction on the training points. (a) CFRP specimen. (b) Reconstruction results.
(a)
(b)


#### 3.1.4. Time Domain Signal Reconstruction of Linear Explosive Induced Pyroshock

The time domain signal of pyroshock was reconstructed by substituting the information fromSRSref(x,y,fn) to Lx,y,t via the method described in Section 2.4. The reconstruction results of measurement points 1 and 2 used as training points are shown in Figures 13 and 14. The blue signal represents the reconstruction result and the black signal represents the real pyroshock. At measurement point 1, the time domain signal was reconstructed with a PAD of 8.96% and the SRS was reconstructed with a MAD of 22.04%. At the measurement point 2, the time domain signal was reconstructed with a PAD of 19.43% and the SRS was reconstructed with a MAD of 21.84%. Since we confirmed that the two training points were reconstructed to good similar levels, we obtained PRx,y,t and SRSR(x,y,fn) by reconstructing pyroshock at all points using GSRS,O(fn) obtained in the training process.Figure 13
Time domain signal reconstruction at point 1, training point.Figure 14
Time domain signal reconstruction at point 2, training point.


#### 3.1.5. Pyroshock Wave Propagation Image and Shock Response Spectrum Image of Linear Explosive Induced Pyroshock

The obtainedPRx,y,t was imaged in color scale at each point along the time and PWPI was obtained as shown in Figure 15(a). Figure 15(a) shows a snapshot of the PWPI at 467 μs. We can confirm propagation characteristics such as propagation path and pyroshock reflection wave via PWPI. Additionally, PWPI allowed us to see how much of each point was affected by shock waves over time. When the amplitude of SRSR(x,y,fn) was imaged in color scale as frequency increases, SRSI was obtained. Through the SRSI, we could ascertain how much the structure was affected by the pyroshock for each frequency. For example, in Figure 15(b), we could see how much each point was affected by the pyroshock with 10 kHz frequency components. We expect to be able to use this SRSI to avoid dangerous regions for the electronic components which are sensitive to certain frequency ranges.Figure 15
Linear explosive induced pyroshock propagation evaluation; (a) PWPI at 467μs; (b) SRSI at 10 kHz.
(a)
(b)


#### 3.1.6. Verification of Linear Explosive Induced Pyroshock Reconstruction Result

To verify the obtained PWPI and SRSI, we compared the reconstruction result with the real pyroshock at the measurement point 3 set as the verification point. The comparison result is shown in Figure16. The black signal is the real pyroshock, and the blue signal is the reconstruction result. At the verification point, the time domain signal was reconstructed with a PAD of 19.83%, and the SRS was reconstructed with a MAD of 21.61%. Using these reconstruction results, we could confirm the attenuation characteristics of pyroshock over time and distance. Figure 17(b) could be obtained by plotting SRSR(x,y,fn) along the red line shown in Figure 17(a). As shown in Figure 17(b), it could be seen that the SRS attenuates exponentially with distance from the explosion source. Figure 18 could be obtained by plotting PRx,y,t along the red line shown in Figure 17(a). As shown in Figure 18, it could be seen that the time domain signal attenuates exponentially with time and distance from the explosion source.Figure 16
Time domain signal reconstruction at point 3, verification point.Figure 17
(a) CFRP specimen with red line. (b) Reconstructed SRS along the red line.
(a)
(b)Figure 18
Reconstructed time domain signal along the red line.Reconstruction using measurement points 1 and 2 as training points and measurement point 3 as verification point confirmed that pyroshock was reconstructed with averaged PAD of 16.07% and MAD of 21.83%. To verify more closely the developed algorithm, this process was repeatedly conducted by changing the training and verification points. Figure19 shows the results of reconstruction using measurement points 1 and 3 as training points and measurement point 2 as verification point. Figure 20 shows the results of reconstruction using measurement points 2 and 3 as training points and measurement point 1 as verification points. All three verifications were reconstructed at similar levels in PAD and MAD.Figure 19
Time domain signal reconstruction result in using point 2 as a verification point.Figure 20
Time domain signal reconstruction result in using point 1 as a verification point.


### 3.2. Reconstruction of Point Explosive Induced Pyroshock in Disk Type Specimen

#### 3.2.1. Point Explosive-Induced Pyroshock Measurement in Disk Type Specimen

In this section, we analyzed the propagation of the point explosive induced pyroshock in specimens of steel 4340 material. The specimen had a thickness of 27 mm and a diameter of 300 mm. The pyroshock measurement setup is shown in Figure21. Pyroshock was induced by a 1/2-inch explosive bolt as shown in Figure 21 and measured by three LDVs. The pyroshock was measured at 10 mm, 50 mm, and 90 mm from the explosion source, and detailed DAQ conditions are shown in Table 2.Table 2
Measurement condition of point explosive induced pyroshock on the disk type specimen.
Shock source
CF
F
s (MHz)
Filter (kHz)
Number of samples
C
F
m
r (mm/s/V)
A
d
B (dB)
Pyroshock
C
F
1 = 1.42CF2 = 1.42CF3 = 1.42
1
0.1~100
30,000
1,000
A
d
B
1 = 9AdB2 = 6AdB3 = 6Figure 21
Measurement system of point explosive-induced pyroshock and the disk type specimen.


#### 3.2.2. Laser Shock Measurement in Disk Type Specimen

Laser shock data were obtained for the same specimen used in pyroshock measurement. The laser shock generated by Q-switch laser was collected by a PZT sensor. The scanning area was set to include pyroshock measurement points as shown in Figure22. Detailed DAQ conditions are shown in Table 3.Table 3
Measurement condition of laser shock for reconstruction of pyroshock on the disk type specimen.
Shock source
F
s (MHz)
Filter (kHz)
Number of samples
PRF (Hz)
Scan area (mm)
Laser shock
1
0.1~100
512
50
Width = 136Height = 104Interval = 1Figure 22
Scanning measurement system of laser shock for reconstruction of explosive induced pyroshock on the disk type specimen.


#### 3.2.3. Reconstruction of Pyroshock in Disk Type Specimen

The pyroshock reconstruction process was performed using two arbitrary points of the acquired signals, measurement point 2 and point 3 as the training points. The reconstruction results at the training points are shown in Figure23. The blue signal represents the reconstruction result and the black signal represents the real pyroshock signal. The reconstruction results showed that the time domain signal at the measurement point 2 was reconstructed with the PAD of 4.88% and the SRS was reconstructed with the MAD of 29.50%. At measurement point 3, the time domain signal was reconstructed with the PAD of 1.62% and the SRS was reconstructed at a similar level of MAD 26.57%. Figure 24 shows PWPI and SRSI. Figure 24(a) shows a snapshot of PWPI at 60 μs, and Figure 24(b) shows a snapshot of SRSI at 40 kHz.Figure 23
Reconstruction results of pyroshock on the disk type specimen at the training points.Figure 24
Pyroshock propagation evaluation on the disk type specimen; (a) PWPI at 60μs. (b) SRSI at 40 kHz.
(a)
(b)To verify the algorithm, we compared the reconstruction result and the real pyroshock at the measurement point 1 set as the verification point as shown in Figure25. Because of comparison at the verification point, the time domain signal was reconstructed with the PAD of 20.59% and the SRS was reconstructed with the MAD of 30.11%.Figure 25
Reconstruction results of pyroshock on the disk type specimen at the verification point.


### 3.3. Comparison of Developed Algorithm and Previous Algorithms

Reconstruction of pyroshock using laser shock was first attempted in [18]. In [18], SRS reconstruction results have not been accurate because they tried to compensate for the difference in characteristics between pyroshock and laser shock by simply matching the peak to peak acceleration in the time domain. In [19, 20], these researchers have developed an algorithm to compensate for the difference in characteristics of laser shock and pyroshock by statistically predicting the amplitude of pyroshock with distance from the explosion source. They have achieved accurate results by reconstructing point explosive-induced pyroshock on an aluminum plate. However, since the pyroshock is affected not only by the distance from the explosion source but also by the geometry of the structure, the algorithm proposed by [19, 20] did not have any mechanism needed to predict pyroshock in complex structures. Figure 26 shows the reconstruction results of the algorithm developed in this paper compared with the previously developed algorithm in [18–20]. Figure 26 shows the reconstruction results of the point explosive induced pyroshock in disk type specimen through each system. Reconstruction was performed using only point 2 and point 3 in Figure 22; point 1 was used as a verification point to verify the algorithm, and Figure 26 shows the reconstruction results at verification point. In each result, the black signal represents the real pyroshock, the blue signal represents the reconstruction result, and the red signal represents the laser shock before gaining. Figure 26(a) shows the reconstruction result using the algorithm presented in [18]. From the results of the SRS reconstruction, only the peak value was adjusted, and the SRS shape was maintained as it is. This indicates that the reconstruction result will be not good if the SRS shape of laser shock differs from the SRS shape of pyroshock. On the other hand, in the case of the algorithm presented in this paper, since the characteristic difference is corrected for each frequency range, the reconstruction result was excellent even though the SRS shape of laser shock differs from the SRS shape of pyroshock as shown in Figure 26(c). Figure 26(b) shows the reconstruction result using the algorithm presented in [19, 20]. We can see that the reconstruction results were not good because of the problem where the algorithm proposed in [19, 20] assumes that the amplitude of the pyroshock is only dependent on the distance from the explosion source. However, the pyroshock not only was dependent on the distance from the explosion source but also depends on the geometry of the structure. Thus, the algorithm presented in [19, 20] cannot be successfully used for complex and local structures which do not allow a propagation modelling based exponential attenuation. On the other hand, in the case of the algorithm presented in this paper, it can be confirmed that various types of pyroshock were excellently reconstructed as shown in Table 4.Table 4
Reconstruction results and comparison of the proposed algorithm and previous algorithms.
Developed algorithm
Ref [18]
Ref [19, 20]
Training point
Verification point
Training point
Verification point
Training point
Verification point
MAD (%)
PAD (%)
MAD (%)
PAD (%)
MAD (%)
PAD (%)
MAD (%)
PAD (%)
MAD (%)
PAD (%)
MAD (%)
PAD (%)
CFRP-linear
21.94
14.02
21.61
19.83
31.87
20.55
34.51
19.10
Could not be reconstructed
Disk-point
28.04
3.25
30.11
20.59
45.65
7.30
52.64
18.03
36.40
45.30
63.22
48.33
Average
24.99
8.64
25.86
20.21
38.76
13.93
43.58
18.57
36.4
45.3
63.22
48.33Figure 26
Comparison of reconstruction results of pyroshock on the disk type specimen at the verification point. (a) Ref [18]. (b) Ref [19, 20]. (c) The algorithm proposed in this paper.
(a)
(b)
(c)


## 3.1. Reconstruction of Linear Explosive Induced Pyroshock in CFRP Specimen

In this section, pyroshock data in ref [16] was secured to verify and compare the performance of the developed algorithm (Table 1). The pyroshock data were obtained by exploding the linear pyrotechnical device in the quasi-isotropic CFRP sandwich panel. The laser shock was obtained for the same CFRP specimen, and the PWPI and SRSI were obtained using the developed algorithm. The verification of the algorithm was conducted by comparing the real pyroshock with the reconstructed pyroshock. We describe each step of Figure 1 in detail to help understand the reconstruction process.Table 1
Laser shock measurement condition for reconstruction of linear explosive induced pyroshock.
Shock source
F
s (MHz)
Filter (kHz)
Number of samples
PRF (Hz)
Scan area (mm)
Laser shock
1.25
0.1~100
600
1,000
Width = 640Height = 640Interval = 2


### 3.1.1. Linear Explosive Induced Pyroshock Measurement in CFRP Specimen

The pyroshock was generated by exploding a 230-mm linear pyrotechnical device on a CFRP specimen as shown in Figure9(a). Three LDV sensors were used to measure pyroshock at each measurement point indicated in Figure 9(b). The measured pyroshock was converted into a time domain signal in acceleration as shown in Figure 9(c) through the conversion process and stored in the PC.Figure 9
Linear explosive induced pyroshock measurement on the CFRP specimen. (a) Experimental setup. (b) Sensing point on the test specimen. (c) Time domain signal of linear explosive induced pyroshock.
(a)
(b)
(c)


### 3.1.2. Laser Shock Scanning Measurement in CFRP Specimen

The same specimen used in pyroshock measurement was installed as shown in Figure10(a), and the laser scanning process was performed along the scanning path as shown in Figure 10(b). When the laser scanning process was finished, the laser scanning data of the 4-dimensional shape was obtained as shown in Figure 10(c). To simulate 230 mm-long linear pyrotechnical device, eleven PZT sensors were installed at intervals of 23 mm, as shown in Figure 11(a). The pyrotechnical device used in this paper started to explode from the right and had a principal propagation direction of 19.6 degree [18], so we phase arrayed the laser scanning data and synthesized it as shown in Figure 11(b). The phased-array laser scanning data were used as Lx,y,t. Figure 11(c) showed the extraction of laser shock data at the same point as the pyroshock measurement point.Figure 10
(a) Laser shock measurement setup. (b) Laser scanning process and laser shock time domain signal at an arbitrary point. (c) Four-dimensional laser scanning data structure.
(a)
(b)
(c)Figure 11
(a) Snapshots of laser scanning data at 65μs obtained by each sensor. (b) Snapshot of synthesized laser scanning data at 100 μs. (c) Training point laser shock, L(xs,ys,t).


### 3.1.3. Shock Response Spectrum Reconstruction of Linear Explosive Induced Pyroshock

After completion of the laser shock and pyroshock data acquisition, the SRS reconstruction process was performed as described in Section2.3. We used the measurement points 1 and 2, as shown in Figure 12(a), as the training points to reconstruct the SRS. Measurement point 3 was used as the verification point to verify the reconstruction results. In other words, the real pyroshock signals obtained in the only measurement points 1 and 2 were used to train the algorithm and the algorithm provided the reconstruction results and then the results were compared with the real pyroshock at the measurement point 3. The reconstruction results at the training points are shown in Figure 12(b). The black signal represents the SRS of the real pyroshock and the blue signal represents the SRS obtained through the reconstruction process. At the point 1, the SRS was trained with a MAD of 4.34%. At the point 2, the SRS was trained with a MAD of 20.80%.Figure 12
SRS reconstruction on the training points. (a) CFRP specimen. (b) Reconstruction results.
(a)
(b)


### 3.1.4. Time Domain Signal Reconstruction of Linear Explosive Induced Pyroshock

The time domain signal of pyroshock was reconstructed by substituting the information fromSRSref(x,y,fn) to Lx,y,t via the method described in Section 2.4. The reconstruction results of measurement points 1 and 2 used as training points are shown in Figures 13 and 14. The blue signal represents the reconstruction result and the black signal represents the real pyroshock. At measurement point 1, the time domain signal was reconstructed with a PAD of 8.96% and the SRS was reconstructed with a MAD of 22.04%. At the measurement point 2, the time domain signal was reconstructed with a PAD of 19.43% and the SRS was reconstructed with a MAD of 21.84%. Since we confirmed that the two training points were reconstructed to good similar levels, we obtained PRx,y,t and SRSR(x,y,fn) by reconstructing pyroshock at all points using GSRS,O(fn) obtained in the training process.Figure 13
Time domain signal reconstruction at point 1, training point.Figure 14
Time domain signal reconstruction at point 2, training point.


### 3.1.5. Pyroshock Wave Propagation Image and Shock Response Spectrum Image of Linear Explosive Induced Pyroshock

The obtainedPRx,y,t was imaged in color scale at each point along the time and PWPI was obtained as shown in Figure 15(a). Figure 15(a) shows a snapshot of the PWPI at 467 μs. We can confirm propagation characteristics such as propagation path and pyroshock reflection wave via PWPI. Additionally, PWPI allowed us to see how much of each point was affected by shock waves over time. When the amplitude of SRSR(x,y,fn) was imaged in color scale as frequency increases, SRSI was obtained. Through the SRSI, we could ascertain how much the structure was affected by the pyroshock for each frequency. For example, in Figure 15(b), we could see how much each point was affected by the pyroshock with 10 kHz frequency components. We expect to be able to use this SRSI to avoid dangerous regions for the electronic components which are sensitive to certain frequency ranges.Figure 15
Linear explosive induced pyroshock propagation evaluation; (a) PWPI at 467μs; (b) SRSI at 10 kHz.
(a)
(b)


### 3.1.6. Verification of Linear Explosive Induced Pyroshock Reconstruction Result

To verify the obtained PWPI and SRSI, we compared the reconstruction result with the real pyroshock at the measurement point 3 set as the verification point. The comparison result is shown in Figure16. The black signal is the real pyroshock, and the blue signal is the reconstruction result. At the verification point, the time domain signal was reconstructed with a PAD of 19.83%, and the SRS was reconstructed with a MAD of 21.61%. Using these reconstruction results, we could confirm the attenuation characteristics of pyroshock over time and distance. Figure 17(b) could be obtained by plotting SRSR(x,y,fn) along the red line shown in Figure 17(a). As shown in Figure 17(b), it could be seen that the SRS attenuates exponentially with distance from the explosion source. Figure 18 could be obtained by plotting PRx,y,t along the red line shown in Figure 17(a). As shown in Figure 18, it could be seen that the time domain signal attenuates exponentially with time and distance from the explosion source.Figure 16
Time domain signal reconstruction at point 3, verification point.Figure 17
(a) CFRP specimen with red line. (b) Reconstructed SRS along the red line.
(a)
(b)Figure 18
Reconstructed time domain signal along the red line.Reconstruction using measurement points 1 and 2 as training points and measurement point 3 as verification point confirmed that pyroshock was reconstructed with averaged PAD of 16.07% and MAD of 21.83%. To verify more closely the developed algorithm, this process was repeatedly conducted by changing the training and verification points. Figure19 shows the results of reconstruction using measurement points 1 and 3 as training points and measurement point 2 as verification point. Figure 20 shows the results of reconstruction using measurement points 2 and 3 as training points and measurement point 1 as verification points. All three verifications were reconstructed at similar levels in PAD and MAD.Figure 19
Time domain signal reconstruction result in using point 2 as a verification point.Figure 20
Time domain signal reconstruction result in using point 1 as a verification point.


## 3.1.1. Linear Explosive Induced Pyroshock Measurement in CFRP Specimen

The pyroshock was generated by exploding a 230-mm linear pyrotechnical device on a CFRP specimen as shown in Figure9(a). Three LDV sensors were used to measure pyroshock at each measurement point indicated in Figure 9(b). The measured pyroshock was converted into a time domain signal in acceleration as shown in Figure 9(c) through the conversion process and stored in the PC.Figure 9
Linear explosive induced pyroshock measurement on the CFRP specimen. (a) Experimental setup. (b) Sensing point on the test specimen. (c) Time domain signal of linear explosive induced pyroshock.
(a)
(b)
(c)


## 3.1.2. Laser Shock Scanning Measurement in CFRP Specimen

The same specimen used in pyroshock measurement was installed as shown in Figure10(a), and the laser scanning process was performed along the scanning path as shown in Figure 10(b). When the laser scanning process was finished, the laser scanning data of the 4-dimensional shape was obtained as shown in Figure 10(c). To simulate 230 mm-long linear pyrotechnical device, eleven PZT sensors were installed at intervals of 23 mm, as shown in Figure 11(a). The pyrotechnical device used in this paper started to explode from the right and had a principal propagation direction of 19.6 degree [18], so we phase arrayed the laser scanning data and synthesized it as shown in Figure 11(b). The phased-array laser scanning data were used as Lx,y,t. Figure 11(c) showed the extraction of laser shock data at the same point as the pyroshock measurement point.Figure 10
(a) Laser shock measurement setup. (b) Laser scanning process and laser shock time domain signal at an arbitrary point. (c) Four-dimensional laser scanning data structure.
(a)
(b)
(c)Figure 11
(a) Snapshots of laser scanning data at 65μs obtained by each sensor. (b) Snapshot of synthesized laser scanning data at 100 μs. (c) Training point laser shock, L(xs,ys,t).


## 3.1.3. Shock Response Spectrum Reconstruction of Linear Explosive Induced Pyroshock

After completion of the laser shock and pyroshock data acquisition, the SRS reconstruction process was performed as described in Section2.3. We used the measurement points 1 and 2, as shown in Figure 12(a), as the training points to reconstruct the SRS. Measurement point 3 was used as the verification point to verify the reconstruction results. In other words, the real pyroshock signals obtained in the only measurement points 1 and 2 were used to train the algorithm and the algorithm provided the reconstruction results and then the results were compared with the real pyroshock at the measurement point 3. The reconstruction results at the training points are shown in Figure 12(b). The black signal represents the SRS of the real pyroshock and the blue signal represents the SRS obtained through the reconstruction process. At the point 1, the SRS was trained with a MAD of 4.34%. At the point 2, the SRS was trained with a MAD of 20.80%.Figure 12
SRS reconstruction on the training points. (a) CFRP specimen. (b) Reconstruction results.
(a)
(b)


## 3.1.4. Time Domain Signal Reconstruction of Linear Explosive Induced Pyroshock

The time domain signal of pyroshock was reconstructed by substituting the information fromSRSref(x,y,fn) to Lx,y,t via the method described in Section 2.4. The reconstruction results of measurement points 1 and 2 used as training points are shown in Figures 13 and 14. The blue signal represents the reconstruction result and the black signal represents the real pyroshock. At measurement point 1, the time domain signal was reconstructed with a PAD of 8.96% and the SRS was reconstructed with a MAD of 22.04%. At the measurement point 2, the time domain signal was reconstructed with a PAD of 19.43% and the SRS was reconstructed with a MAD of 21.84%. Since we confirmed that the two training points were reconstructed to good similar levels, we obtained PRx,y,t and SRSR(x,y,fn) by reconstructing pyroshock at all points using GSRS,O(fn) obtained in the training process.Figure 13
Time domain signal reconstruction at point 1, training point.Figure 14
Time domain signal reconstruction at point 2, training point.


## 3.1.5. Pyroshock Wave Propagation Image and Shock Response Spectrum Image of Linear Explosive Induced Pyroshock

The obtainedPRx,y,t was imaged in color scale at each point along the time and PWPI was obtained as shown in Figure 15(a). Figure 15(a) shows a snapshot of the PWPI at 467 μs. We can confirm propagation characteristics such as propagation path and pyroshock reflection wave via PWPI. Additionally, PWPI allowed us to see how much of each point was affected by shock waves over time. When the amplitude of SRSR(x,y,fn) was imaged in color scale as frequency increases, SRSI was obtained. Through the SRSI, we could ascertain how much the structure was affected by the pyroshock for each frequency. For example, in Figure 15(b), we could see how much each point was affected by the pyroshock with 10 kHz frequency components. We expect to be able to use this SRSI to avoid dangerous regions for the electronic components which are sensitive to certain frequency ranges.Figure 15
Linear explosive induced pyroshock propagation evaluation; (a) PWPI at 467μs; (b) SRSI at 10 kHz.
(a)
(b)


## 3.1.6. Verification of Linear Explosive Induced Pyroshock Reconstruction Result

To verify the obtained PWPI and SRSI, we compared the reconstruction result with the real pyroshock at the measurement point 3 set as the verification point. The comparison result is shown in Figure16. The black signal is the real pyroshock, and the blue signal is the reconstruction result. At the verification point, the time domain signal was reconstructed with a PAD of 19.83%, and the SRS was reconstructed with a MAD of 21.61%. Using these reconstruction results, we could confirm the attenuation characteristics of pyroshock over time and distance. Figure 17(b) could be obtained by plotting SRSR(x,y,fn) along the red line shown in Figure 17(a). As shown in Figure 17(b), it could be seen that the SRS attenuates exponentially with distance from the explosion source. Figure 18 could be obtained by plotting PRx,y,t along the red line shown in Figure 17(a). As shown in Figure 18, it could be seen that the time domain signal attenuates exponentially with time and distance from the explosion source.Figure 16
Time domain signal reconstruction at point 3, verification point.Figure 17
(a) CFRP specimen with red line. (b) Reconstructed SRS along the red line.
(a)
(b)Figure 18
Reconstructed time domain signal along the red line.Reconstruction using measurement points 1 and 2 as training points and measurement point 3 as verification point confirmed that pyroshock was reconstructed with averaged PAD of 16.07% and MAD of 21.83%. To verify more closely the developed algorithm, this process was repeatedly conducted by changing the training and verification points. Figure19 shows the results of reconstruction using measurement points 1 and 3 as training points and measurement point 2 as verification point. Figure 20 shows the results of reconstruction using measurement points 2 and 3 as training points and measurement point 1 as verification points. All three verifications were reconstructed at similar levels in PAD and MAD.Figure 19
Time domain signal reconstruction result in using point 2 as a verification point.Figure 20
Time domain signal reconstruction result in using point 1 as a verification point.


## 3.2. Reconstruction of Point Explosive Induced Pyroshock in Disk Type Specimen

### 3.2.1. Point Explosive-Induced Pyroshock Measurement in Disk Type Specimen

In this section, we analyzed the propagation of the point explosive induced pyroshock in specimens of steel 4340 material. The specimen had a thickness of 27 mm and a diameter of 300 mm. The pyroshock measurement setup is shown in Figure21. Pyroshock was induced by a 1/2-inch explosive bolt as shown in Figure 21 and measured by three LDVs. The pyroshock was measured at 10 mm, 50 mm, and 90 mm from the explosion source, and detailed DAQ conditions are shown in Table 2.Table 2
Measurement condition of point explosive induced pyroshock on the disk type specimen.
Shock source
CF
F
s (MHz)
Filter (kHz)
Number of samples
C
F
m
r (mm/s/V)
A
d
B (dB)
Pyroshock
C
F
1 = 1.42CF2 = 1.42CF3 = 1.42
1
0.1~100
30,000
1,000
A
d
B
1 = 9AdB2 = 6AdB3 = 6Figure 21
Measurement system of point explosive-induced pyroshock and the disk type specimen.


### 3.2.2. Laser Shock Measurement in Disk Type Specimen

Laser shock data were obtained for the same specimen used in pyroshock measurement. The laser shock generated by Q-switch laser was collected by a PZT sensor. The scanning area was set to include pyroshock measurement points as shown in Figure22. Detailed DAQ conditions are shown in Table 3.Table 3
Measurement condition of laser shock for reconstruction of pyroshock on the disk type specimen.
Shock source
F
s (MHz)
Filter (kHz)
Number of samples
PRF (Hz)
Scan area (mm)
Laser shock
1
0.1~100
512
50
Width = 136Height = 104Interval = 1Figure 22
Scanning measurement system of laser shock for reconstruction of explosive induced pyroshock on the disk type specimen.


### 3.2.3. Reconstruction of Pyroshock in Disk Type Specimen

The pyroshock reconstruction process was performed using two arbitrary points of the acquired signals, measurement point 2 and point 3 as the training points. The reconstruction results at the training points are shown in Figure23. The blue signal represents the reconstruction result and the black signal represents the real pyroshock signal. The reconstruction results showed that the time domain signal at the measurement point 2 was reconstructed with the PAD of 4.88% and the SRS was reconstructed with the MAD of 29.50%. At measurement point 3, the time domain signal was reconstructed with the PAD of 1.62% and the SRS was reconstructed at a similar level of MAD 26.57%. Figure 24 shows PWPI and SRSI. Figure 24(a) shows a snapshot of PWPI at 60 μs, and Figure 24(b) shows a snapshot of SRSI at 40 kHz.Figure 23
Reconstruction results of pyroshock on the disk type specimen at the training points.Figure 24
Pyroshock propagation evaluation on the disk type specimen; (a) PWPI at 60μs. (b) SRSI at 40 kHz.
(a)
(b)To verify the algorithm, we compared the reconstruction result and the real pyroshock at the measurement point 1 set as the verification point as shown in Figure25. Because of comparison at the verification point, the time domain signal was reconstructed with the PAD of 20.59% and the SRS was reconstructed with the MAD of 30.11%.Figure 25
Reconstruction results of pyroshock on the disk type specimen at the verification point.


## 3.2.1. Point Explosive-Induced Pyroshock Measurement in Disk Type Specimen

In this section, we analyzed the propagation of the point explosive induced pyroshock in specimens of steel 4340 material. The specimen had a thickness of 27 mm and a diameter of 300 mm. The pyroshock measurement setup is shown in Figure21. Pyroshock was induced by a 1/2-inch explosive bolt as shown in Figure 21 and measured by three LDVs. The pyroshock was measured at 10 mm, 50 mm, and 90 mm from the explosion source, and detailed DAQ conditions are shown in Table 2.Table 2
Measurement condition of point explosive induced pyroshock on the disk type specimen.
Shock source
CF
F
s (MHz)
Filter (kHz)
Number of samples
C
F
m
r (mm/s/V)
A
d
B (dB)
Pyroshock
C
F
1 = 1.42CF2 = 1.42CF3 = 1.42
1
0.1~100
30,000
1,000
A
d
B
1 = 9AdB2 = 6AdB3 = 6Figure 21
Measurement system of point explosive-induced pyroshock and the disk type specimen.


## 3.2.2. Laser Shock Measurement in Disk Type Specimen

Laser shock data were obtained for the same specimen used in pyroshock measurement. The laser shock generated by Q-switch laser was collected by a PZT sensor. The scanning area was set to include pyroshock measurement points as shown in Figure22. Detailed DAQ conditions are shown in Table 3.Table 3
Measurement condition of laser shock for reconstruction of pyroshock on the disk type specimen.
Shock source
F
s (MHz)
Filter (kHz)
Number of samples
PRF (Hz)
Scan area (mm)
Laser shock
1
0.1~100
512
50
Width = 136Height = 104Interval = 1Figure 22
Scanning measurement system of laser shock for reconstruction of explosive induced pyroshock on the disk type specimen.


## 3.2.3. Reconstruction of Pyroshock in Disk Type Specimen

The pyroshock reconstruction process was performed using two arbitrary points of the acquired signals, measurement point 2 and point 3 as the training points. The reconstruction results at the training points are shown in Figure23. The blue signal represents the reconstruction result and the black signal represents the real pyroshock signal. The reconstruction results showed that the time domain signal at the measurement point 2 was reconstructed with the PAD of 4.88% and the SRS was reconstructed with the MAD of 29.50%. At measurement point 3, the time domain signal was reconstructed with the PAD of 1.62% and the SRS was reconstructed at a similar level of MAD 26.57%. Figure 24 shows PWPI and SRSI. Figure 24(a) shows a snapshot of PWPI at 60 μs, and Figure 24(b) shows a snapshot of SRSI at 40 kHz.Figure 23
Reconstruction results of pyroshock on the disk type specimen at the training points.Figure 24
Pyroshock propagation evaluation on the disk type specimen; (a) PWPI at 60μs. (b) SRSI at 40 kHz.
(a)
(b)To verify the algorithm, we compared the reconstruction result and the real pyroshock at the measurement point 1 set as the verification point as shown in Figure25. Because of comparison at the verification point, the time domain signal was reconstructed with the PAD of 20.59% and the SRS was reconstructed with the MAD of 30.11%.Figure 25
Reconstruction results of pyroshock on the disk type specimen at the verification point.


## 3.3. Comparison of Developed Algorithm and Previous Algorithms

Reconstruction of pyroshock using laser shock was first attempted in [18]. In [18], SRS reconstruction results have not been accurate because they tried to compensate for the difference in characteristics between pyroshock and laser shock by simply matching the peak to peak acceleration in the time domain. In [19, 20], these researchers have developed an algorithm to compensate for the difference in characteristics of laser shock and pyroshock by statistically predicting the amplitude of pyroshock with distance from the explosion source. They have achieved accurate results by reconstructing point explosive-induced pyroshock on an aluminum plate. However, since the pyroshock is affected not only by the distance from the explosion source but also by the geometry of the structure, the algorithm proposed by [19, 20] did not have any mechanism needed to predict pyroshock in complex structures. Figure 26 shows the reconstruction results of the algorithm developed in this paper compared with the previously developed algorithm in [18–20]. Figure 26 shows the reconstruction results of the point explosive induced pyroshock in disk type specimen through each system. Reconstruction was performed using only point 2 and point 3 in Figure 22; point 1 was used as a verification point to verify the algorithm, and Figure 26 shows the reconstruction results at verification point. In each result, the black signal represents the real pyroshock, the blue signal represents the reconstruction result, and the red signal represents the laser shock before gaining. Figure 26(a) shows the reconstruction result using the algorithm presented in [18]. From the results of the SRS reconstruction, only the peak value was adjusted, and the SRS shape was maintained as it is. This indicates that the reconstruction result will be not good if the SRS shape of laser shock differs from the SRS shape of pyroshock. On the other hand, in the case of the algorithm presented in this paper, since the characteristic difference is corrected for each frequency range, the reconstruction result was excellent even though the SRS shape of laser shock differs from the SRS shape of pyroshock as shown in Figure 26(c). Figure 26(b) shows the reconstruction result using the algorithm presented in [19, 20]. We can see that the reconstruction results were not good because of the problem where the algorithm proposed in [19, 20] assumes that the amplitude of the pyroshock is only dependent on the distance from the explosion source. However, the pyroshock not only was dependent on the distance from the explosion source but also depends on the geometry of the structure. Thus, the algorithm presented in [19, 20] cannot be successfully used for complex and local structures which do not allow a propagation modelling based exponential attenuation. On the other hand, in the case of the algorithm presented in this paper, it can be confirmed that various types of pyroshock were excellently reconstructed as shown in Table 4.Table 4
Reconstruction results and comparison of the proposed algorithm and previous algorithms.
Developed algorithm
Ref [18]
Ref [19, 20]
Training point
Verification point
Training point
Verification point
Training point
Verification point
MAD (%)
PAD (%)
MAD (%)
PAD (%)
MAD (%)
PAD (%)
MAD (%)
PAD (%)
MAD (%)
PAD (%)
MAD (%)
PAD (%)
CFRP-linear
21.94
14.02
21.61
19.83
31.87
20.55
34.51
19.10
Could not be reconstructed
Disk-point
28.04
3.25
30.11
20.59
45.65
7.30
52.64
18.03
36.40
45.30
63.22
48.33
Average
24.99
8.64
25.86
20.21
38.76
13.93
43.58
18.57
36.4
45.3
63.22
48.33Figure 26
Comparison of reconstruction results of pyroshock on the disk type specimen at the verification point. (a) Ref [18]. (b) Ref [19, 20]. (c) The algorithm proposed in this paper.
(a)
(b)
(c)


## 4. Conclusions

In this paper, we developed a pyroshock reconstruction algorithm using laser shock. The developed system used a laser scanning technique to acquire laser shock and reconstructed it with pyroshock using iterative signal decomposition and synthesis method. This algorithm obtained the gain to compensate the difference in characteristics between laser shock and pyroshock based on SRS. The gain was applied to laser shock signal decomposed by Butterworth filter in each frequency band. The stopband effect in band decomposition was compensated by iteration to improve the reconstruction result. We reconstructed linear explosive induced pyroshock using the developed algorithm as well as point explosive. PWPI and SRSI were obtained through the reconstruction. We could confirm propagation characteristics such as propagation path and reflection wave of pyroshock through PWPI. Through the SRSI, we could ascertain how much the structure was affected by the pyroshock for each frequency. To verify the algorithm, we compared the reconstruction result with the real pyroshock. At the verification point, the time domain signal was reconstructed with an averaged PAD of 20.21% and the SRS was reconstructed with an average MAD of 25.86%. This was a remarkably improved result compared to the developed pyroshock reconstruction system in [18–20]. Additionally, according to [21], the repeatability of the pyroshock obtained through the same pyro device in the same structure had a MAD of 19.6%. Considering this fact, the pyroshock reconstruction algorithm developed in this paper provided sufficiently reliable reconstruction results. The developed system can be used to determine the most suitable location for installing electronic equipment with a specific resonance frequency. In addition, it can be used to check the propagation path of pyroshock and effectively reduce the pyroshock transmitted to the structure. To reduce the economic burden caused by the pyro device, the developed system can be applied to the evaluation of the multi-explosive induced pyroshock using only the single pyro device test by applying the wave superposition principle. In addition, it is expected that the iterative signal decomposition and synthesis method proposed in this paper can be used widely because it effectively removes the stop band effect caused by the signal decomposition process using Butterworth filter.



---

*Source: 8351791-2017-10-26.xml*