張小青，王 馳，李金輝，羅欣宇，于瀛潔

（1. 上海大學(xué) 精密機(jī)械工程系, 上海200444；2. 同濟(jì)大學(xué)浙江學(xué)院 電子與信息工程系, 嘉興314000；3. 近地面探測(cè)技術(shù)重點(diǎn)實(shí)驗(yàn)室, 無錫214035）

1 Introduction

As a kind of low-cost defense weapon in wartime, landmines have been widely used in previous wars. Landmines have always been a severe threat to the people’s live of those mined countries or regions[1]. Therefore, eliminating hidden dangers caused by landmines has become a global problem that all countries in the world face. Landmine detection remains a worldwide problem to this day, especially for the detection of buried plastic landmines.The difference in electrical characteristics between plastic-cased landmines and the surrounding soil is relatively small, making it difficult for commonly used metal-cased landmine detectors to detect plastic-cased landmines. From the 1940s to the 1950s, a large number of plastic-cased landmines began to appear when the metal shells of landmines were gradually replaced by plastic shells[2]. Due to the limitation of these detection mechanisms, it is difficult to distinguish whether the buried objects are plastic-cased landmines or bricks, rocks, and other interfering objects. Using UAV-based Optical Data Fusion, landmine detection probability can be preliminarily estimated[3]. The methods of detecting the chemical composition of explosives such as hyperspectral images[4], Neutron challenge[5], and Raman spectroscopy[6]developed in recent years have strong identification capabilities, but the system is complicated and expensive, and the detected signal is extremely weak, leading to the problems relating to long detection (signal accumulation) time and a high misdetection rate. It is possible to overcome the difficulty of detecting plastic-cased landmines containing very low metal contents by improving the sensitivity of metal landmine detectors.However, many shrapnel, shell and other metal interferences in a minefield can cause high false positives. Biological mine detection methods, such as the passive method using Honeybee-based on-site explosive sampling, can only estimate the total explosive load in a certain area but cannot locate the mine’s location[7].

Acoustic-to-seismic landmine detection technology based on the unique mechanical properties of landmines and the principle of acoustic-to-seismic coupling shows good application prospects, especially in the safe and effective detection of plastic-cased landmines. The acoustic-to-seismic coupling[8-10]refers to when sound waves propagating in the air are incident on the ground. In addition to most of the energy reflected into the air by the ground surface, a small proportion is coupled to propagate underground due to the influence of soil porosity. This forms seismic waves with different compositions such as transverse waves and fast and slow longitudinal waves. The air cavity and fuze in the landmine structure make its acoustic compliance (the deformation caused by per unit stress,commonly known as flexibility) much greater than the acoustic compliance of the surrounding soil. The"soil-landmine" system can produce the equivalent“Mass-spring" resonance phenomenon under the influence of seismic waves. Sabatier J. M and Donskoyet al.[11-13]successively established linear and nonlinear resonance models for acoustic-to-seismic landmine detection; Zagraiet al.[14-15]studied the multi-mode mechanism of landmine vibration response; Albertset al.[16]studied the impact of landmine burial depth on resonance frequency. These studies have indicated the feasibility of acoustic-toseismic landmine detection based on the mechanical properties of landmines and the principle of acoustic-to-seismic coupling. However, due to the low efficiency of acoustic-to-seismic coupling, the ground surface vibration excited by sound waves is still very weak even when the landmine resonates.Measuring the characteristic signals of landmines accurately and quickly has always been a key problem that limits the research of acoustic-to-seismic landmine detection systems.

Laser interferometric vibration detection technology has a good application prospect in acousticto-seismic landmine detection due to its of highly precise and non-contact mechanism. Qiukun Zhanget al.[17]studied that the vibration measurement technology of high-performance optical coherence velocimeter can achieve nanometer-level accuracy; Xiang and Sabatier[18-19]used a single-beam laser Doppler vibrometer for acoustic-to-seismic landmine detection; Wang Chi research group[20-21]studied the landmine’s mode shape based on laser self-mixing vibrometer; Rajeshet al.[22]used the ultrasonic Doppler vibrometer array for landmine detection, which increased the vibration measurement rate to a certain extent, but the array had the disadvantages in its bulky size and high cost. In this paper, based on the multi-model resonance phenomenon of a landmine’s upper casing under the sound wave excitation, laser speckle interferometric vibration measurement technology based on a time average is used to study the relationship between the mode shapes of the landmine’s upper casing and the Bessel signals, providing theoretical evidence for realizing the rapid scanning technology of acoustic-optics landmine detection.

2 Analysis of mode shapes of mine

According to literatures [14-15], the main vibration mode characteristics of landmines are mainly determined by the upper surface of landmines, i.e., the landmine’s upper casing. To facilitate the analysis of the multi-mode vibration phenomena of landmines, the characteristics of the three-dimensional vibration mode shapes of landmines are simplified to the characteristics of the two-dimensional vibration mode shape of the plane structures where the burying depth is assumed to be 0 to remove the restriction of the soil above the landmine. As shown in Fig. 1, the landmine’s upper casing is equivalent to a cylindrical thin circular plate structure in polar coordinatesOrθ, with a radius ofaand a thickness ofh. Two continuous boundary constrained springs are arranged on the boundary of the plate structure atr=a, namely linear displacement springkand rotation constrained springK. For a circular planar structure, the (m,n) mode can be used to represent the vibration mode.

Fig. 1 Analytical object and model of landmine’s multi-modal vibration characteristics. (a) Anti-tank plastic landmine; (b)equivalent cylindrical thin circular plate

In a certain 1storder state, wheremrepresents the number of pitch diameters of the structure, andnrepresents the number of pitch circles of the structure. Theoretically, the point with a value of 0 in the mode shape is called a node, i.e., the intersection of the mode shape and the undeformed point in the original structure. The pitch diameter refers to the diameter with a value of 0 in the vibration mode,i.e., the lines composed of nodes. The pitch circle refers to the circle with a value of 0 in the vibration mode. As shown in Fig. 2 (Color online), when a landmine is excited by an external vibration with a certain frequency, the surface of its upper casing will produce a unique multi-modal vibration phenomenon. In Fig. 2, green represents pitch diameter,red represents pitch circle.

Fig. 2 Contour diagrams of landmine mode shapes. (a) (0,1) mode; (b) (0,2) mode; (c) (1,1) mode

The mode shapes of landmines are thoroughly analyzed according to literatures [14-15]. When the landmine’s upper casing is subjected to external forces, the vibration equation of the thin circular plate can be expressed as:

whereTmnis a time function, andWmn(r,θ) is the mode shape function. Based on the mathematical model of the thin circular plate structure, its vibration function can be written as:

where Jnand Inare the first kind and second kind Bessel functions, andais the radius of the thin circular plate. Let the modal shape function equal to zero, and the following equation can be obtained:

According to equation (4), combined with the values of the points whose amplitudes are 0, whenr=r1,r2,r3,···rm, the value of the mode shape function is 0, andmis the number of pitch circles. When θ=θ1,θ2,θ3···θn, the value of the mode shape function is 0, andnis the number of pitch diameters.CmnandAmnare coefficients that depend onmandn, andλmnis theathpositive root of them-order Bessel function Jn=0.

For elastic boundaries, the boundary conditions are:

MrandVrare the bending moments at the corresponding positions. According to the plate and shell vibration theory, the following expressions apply:

According to equations (7) and (8), equations(5) and (6) can be expressed in terms ofW(a,θ).Then the parametersCmnandAmncan be obtained by substituting equation (3) into equations (5) and (6).

The parameterAmnof the mode shape function can be obtained by the following equation:

where δ is the Kronecker number, and whenm=p,δmp= 1; whenm ≠ p,δmp=0.mMis the mass of the thin circular plate. According to the equation (9),Amncan be expressed as:

The mode shape functionWmn(r,θ) can be obtained by substitutingCmn,Amnandλmninto equation(3). The complex "soil-mine" system can be described by a relatively simple and direct mathematical equation, which provides a theoretical basis for the analysis of the mode shapes of the landmine’s upper casing and its application in the acoustic-toseismic landmine detection system.

3 Bessel fringe mapping method

Electronic Speckle-Shearing Patten Interferometry (ESSPI) is a new technique for measuring displacement derivative developed after electronic speckle interference and is widely used in the field of non-destructive testing[23-25]. For the measurement of the deformation caused by the vibration of the object, laser shearing speckle interferometry is often combined with the time-average method of a CCD camera[26-27]. Consequently, the use of high-energy sound wave pulses or vibration exciters to excite the vibration energy close to the landmine location will cause abnormal changes in the ground surface vibration under the resonant movement of the landmine (Fig. 3). The laser beam is projected to the location where the landmine is buried, and the highresolution CCD camera system is used to record the speckle interference fringes containing the surface vibration information in real time. The interference fringes are analyzed and processed to obtain small vibration signals on the ground surface, and then the existence of buried landmines can be determined.

Fig. 3 Vibration deformation schematic diagram of the soil surface with and without landmine under the excitation of sound waves

The speckle interference fringe obtained by the time-average method of a CCD camera is modulated by the zero-order Bessel function, that is,Bessel fringe[28]. Thus, studying the mapping relationship between the mode shapes of the landmine’s upper casing and the Bessel fringes can further study the information such as the type and size of the landmine, as well as the rapid identification method. The following is the analysis of mapping relationship between the function of the mode shape of the landmine ’s upper casing and the Bessel fringes.

To analyze the mapping relationship between the mode shape of the mine cover and the Bessel fringe, the optical path of laser speckle interferometry measurement on the mode shapes of the landmine’s upper casing is shown in Fig. 4(Color online). Assuming that the shearing direction is in thexdirection. LaserS(xs,ys,zs) irradiates on the measured landmine ’s upper casing. The two pointsP1(x,y,z) andP2(x+δx,y,z) on the surface of the landmine’s upper casing spaced apart byδxare imaged and interfered with each other at the same pointC(xc,yc,zc) on the CCD camera after passing through the shearing device, whereδxis the shear amount in thexdirection. After the landmine’s upper casing is excited by sound waves, pointP1(x,y,z) moves to(x+u,y+v,z+w) and pointP2(x+δx,y,z) moves to(x+δx+u+δu,y+v+δv,z+w+δw). The distance betweenandis(δu,δv,δw).

Fig. 4 The optical path of laser speckle interferometry measurement on the mode shapes of the landmine’s upper casing

Fig. 5 Vibration diagram under (0,1) mode. (a) Three-dimensional vibration mode diagram; (b) contour map

Fig. 6 Vibration diagram under (0,2) mode. (a) Three-dimensional vibration mode diagram; (b) contour map

Fig. 7 Vibration displacement gradient change under the (0,1) mode. (a) Three-dimensional diagram of the displacement gradient change; (b) contour map

Fig. 8 Vibration displacement gradient change under the (0,2) mode. (a) Three-dimensional diagram of displacement gradient change; (b) contour map

According to our preliminary research and analysis[28], it is assumed that the intensity of speckle interference obtained by the time averaging method based on a CCD camera isI(x,y), and its expression is as follows：

In the formula, 2Re[g1(x,y)(x,y)] is random speckle noise, J0is the zero-order Bessel function of the first kind, and Δ φ(x,y) is the phase difference of the interference laser due to the vibration deformation of the landmine’s upper casing.

In order to facilitate the correlation between the shear speckle interferometry theory and the multimodal vibration characteristics of mines, the vibration mode function is transformed from the polar coordinate system to the Cartesian coordinate systemxOy, that is:x=rcosθ，y=rsinθ .

The equation of the mode shape function of the landmine’s upper casing in the Cartesian coordinate system can be obtained as equation (13) by substituting equation (12) into equation (3).

According to literature [29], the relationship between the pure phase change caused by the outoff-plane displacement of the landmine’s upper casing and the derivative of the out-off-plane displacement is as follows:

In the formula (11), since 1 -J0(Δφ(x,y)) is a value that changes in a period, the resulting speckle interferogram is a light and dark speckle fringe pattern produced by 1 -J0(Δφ(x,y)) modulation, that is Bessel fringe. It is possible to calculate the phase change caused by the vibration of the out-off-plane deformation through observing the obtained Bessel fringe. Then, the amount of change in the gradient of the out-off-plane displacement of the landmine’s upper casing can be obtained. As a result, the vibration deformation of the landmine’s upper casing can be determined. Where there is a constant shear amount, the Bessel fringe order increases as the displacement gradient is larger.

4 Numerical analysis and discussion

The landmine’s upper casing is equivalent to a thin circular plate structure with elastic supports and it has the mode functionWmn(x,y). The numerical analysis is based on the mapping relationship between the mode shapes of landmine and Bessel fringes, which is formula (13).

The simulation parameter setting is shown in Table 1. The radius of the landmine’s upper casinga=13.5 cm, Poisson’s ratio υ=0.33, Young’s modulusE=17×1019Pa, and laser wavelength λ=658 nm. To clearly display the fringe series in the Bessel fringe pattern, the shear amount is set to 6 mm, and the Cartesian coordinate system is established with the center of the thin circular plate as the origin. The shearing direction of the laser speckle interference is in thexdirection. Assuming that there is a circle and 0 pitch diameters in the thin circular plate, i.e., the (0,1) order mode. The simulation results are shown in Fig 5, 7, 9(a) (Color online), while the mode simulation results for the (0,2)order mode are shown in Figs. 6, 8, 9(b) (Color online).

Tab. 1 Simulation parameter setting

Figs. 5 and 6 show the (0,1) and (0,2) order landmine vibration patterns, respectively. The vibration amplitude of the (0,1) order mode decreases from the center to the surroundings with a maximum value of 2.157 3 mm at the center of the mine.Due to the existence of two pitch circles in the (0,2)order mode, the vibration amplitude changes suddenly at the first pitch circle line then increases from there to the center. Its amplitude reaches a maximum of 2.809 9 mm at the center of the mine, which is an increase of 0.652 6 mm compared to the maximum amplitude of the (0,1) order mode.

Figs. 7 and 8 show the (0,1) and (0,2) order modes’ vibration displacement gradient changes, respectively. Since the shear direction of the laser speckle interference is set to thexdirection, the vibration displacement gradient is obtained by calculating the first-order partial derivative of the landmine vibration mode function with respect tox.Since the mode shapes of landmines are evenly symmetric abouty=0, the displacement gradient change of the (0,1) order mode shape obtained by calculating the first-order partial derivative with respect toxis oddly symmetric abouty=0. The maximum absolute value of the gradient change is 0.228 5. It is obvious from the contour map in Fig. 7(b) that aty=0, the displacement gradient change is always zero regardless of any change inx. The two displacement gradients with odd symmetry increase or decrease from the point where ofx=0 cm andy=±7.752 cm to the surroundings, respectively. As shown in Fig. 8, the first section circular vibration mode of the (0,2) order is evenly symmetrical abouty=0, and its displacement gradient change is oddly symmetrical at that same point. The maximum absolute value of the gradient change is 0.706 4, an increase of approximately 0.477 9 compared to the maximum absolute displacement gradient of the(0,1) order mode. Similarly, aty=0, the displacement gradient change is constantly 0 for any change ofx. The two displacement gradients with odd symmetry increase or decrease from the point with the coordinate wherex=0 cm andy=±4.377 cm to the surroundings, respectively. Compared with the center position of the (0,1) order, theycoordinate of the(0,2) order is offset by 3.375 cm toward the center.

Fig.9 shows the mode shapes of the landmine based on speckle interference fringes of the (0,1)and (0,2) order modes, respectively. It can be seen from Fig. 7 that the displacement gradient generated before and after the interference laser vibration in the (0,1) order mode is oddly symmetrical abouty=0. Aty= 0, with the change ofx, the displacement gradient change is always 0. Therefore, Fig. 9(a)shows symmetrical butterfly-shaped interference fringes with dark fringes aty=0. Since the maximum displacement gradient is 0.228 5 , the Bessel fringe order is about 2. Similarly, the Bessel fringe at the first pitch circle in Fig.9(b) is symmetrically butterfly-shaped wherey=0 with the dark fringe appearing aty=0. Since the maximum displacement gradient is 0.706 4, which is approximately 3 times the maximum displacement gradient of the (0, 1) order mode, the Bessel fringe order is about 6.

Fig. 9 Mine modal shapes based on speckle interference fringes. (a) Under the (0,1) mode; (b) under the (0,2) mode

5 Experiment verification analysis

According to the literatures [14-15], landmines will produce different resonance phenomena in different soil environments, different buried depths,different frequencies, and different decibel levels of sound wave excitation, leading to different degrees of vibration and deformation of the soil surface above buried landmines, i.e., different displacement gradient changes. The order of Bessel fringe series is related to the displacement gradient and the shear amount. Therefore, when using laser speckle interferometry to carry out mine detection experiments,the order of interference fringes such as the (0,1) order mode will change based on the actual measurement situation, but the overall shape should be symmetrically butterfly-shaped interference fringes.

For the purpose of carrying out experimental verification analysis, a landmine detection experimental system based on laser speckle measurement was built up as shown in Fig. 10. The sound source excitation system composed of a function signal generator, a speaker and a power amplifier were used to excite high-intensity low-frequency sound waves. The laser speckle interference detection system was used for speckle interference detection of surface vibrations. The decibel meter was used to measure the sound pressure level of the sound wave.The interference laser wavelength was 658 nm. The type-69 plastic case coach mine, the type-72 antipersonnel coach mine and the Brick were taken as the sample to be tested. The samples were buried in dry or wet sand at a depth of 2 cm. The shear amountδxwas set to 12.122 mm in thexdirection in all experiments, except for the experiment of the type-69 plastic case coach mine buried in wet sand wherein the shear amount was 10.935 mm. The frequencies of the sound were 110 Hz for the type-69 plastic case coach mine and 145 Hz for the type-72 anti-personnel coach mine and brick. The test results are shown in Fig. 11～15 (Color online).

Fig. 10 Acousto-optic mine detection experimental system based on laser speckle measurement

Fig. 11 (a) The type-69 plastic case coach mine in dry sand. Bessel fringes obtained with different excitation sound frequencies and different sound pressure levels. (b) 110 Hz, 100 dB; (c) 110 Hz, 95 dB; (d) 110 Hz, 90 dB

Fig. 12 Bessel fringes of the type-69 plastic case coach mine in wet sand with different excitation sound frequencies and different sound pressure levels. (a) 110 Hz, 100 dB; (b) 110 Hz, 95 dB; (c) 110 Hz, 90 dB

Fig. 13 (a) Type 72 anti-personnel coach mine, and it’s Bessel fringes obtained with different excitation sound frequencies and different sound pressure levels. (b) 145 Hz, 100 dB; (c) 145 Hz, 98 dB; (d) 145 Hz, 95 dB

Fig. 14 (a) Brick in dry sand and (b) it’s speckle interferometer detection results obtained with the excitation sound frequency of 145 Hz and the sound pressure level of 100 dB

Fig. 15 Vibration amplitude of the type-69 plastic case coach mine in dry sand with the excitation sound frequency of 110 Hz and the sound pressure level of 100 dB. (a) Three-dimensional vibration; (b)contour map

As shown in Fig.11, the dotted circle indicates the mine burial area, and no Bessel fringe appears on the ground surface of the area where no landmine is buried. The Bessel fringes are symmetrical butterfly-shaped, which is quantitatively consistent with the theoretical analysis in Section 4. In addition, as the decibel value of the sound pressure level decreases, the deformation amplitude of the ground vibration decreases, the magnitude of change in the displacement gradient decreases, and the number of Bessel fringe levels decreases. The shape and change trend of the Bessel fringes in Fig.12 are consistent with those in Fig.11. In Fig.13, the fringe profile is gossip-like, and the number of Bessel fringe levels also decreases as the decibel value of the sound pressure level decreases. There is no fringe in the area where the brick was buried simultaneously(Fig.14). Therefore, Bessel fringes can reflect the mode shapes of landmines and the displacement gradient changes, which shows the feasibility of using laser speckle interference technology for coupled acoustic-to-seismic landmine detection. The detailed experimental steps and results can be seen in our previous research (Ref. [28]). A three-dimensional vibration mode diagram and contour map are recovered for the type-69 plastic case coach mine in dry sand after performing phase unwrapping(Fig.15). The black dotted line area is the likely location of the landmine. From Fig.15, we can see that the acoustic-to-seismic coupling vibration deformation occurs in the buried landmine area, which is consistent with the theory in Section 3. However,the experimental results only qualitatively verify the relationship between the modal vibration mode and the laser speckle interference signal, the mapping relationship between the mode shape and the Bessel fringe signal under the complex environment is the focus of the next research work.

6 Conclusion

In this paper, the mapping relationship of laser speckle interference fringes is preliminarily established by analyzing landmine mode shapes. The derived mapping relationship shows that the different mode shapes of landmines correspond to the unique Bessel fringes. The Bessel fringes of two modes are simulated. It is verified by simulations that the displacement gradient change caused by mine vibration can generate Bessel fringes. The Bessel fringe series reflects the maximum displacement gradient generated by the vibration deformation. The maximum displacement gradient is increased by 0.477 9 after changing from the (0,1) order mode to the(0,2) order mode, and the Bessel fringe order is increased from 2 to 6. The Bessel fringe shape reflects the entire vibration model. Due to the increase in pitch circle number from the order (0,1)mode to the order (0,2) mode, the Bessel fringe pattern changes from a simple pair of butterfly shapes to a combination of a symmetrical butterfly shape on the first pitch circle, and a circular looped stripe on the second pitch circle. Furthermore, experiments were carried out. Both the numerical and experimental results illustrate the theoretical conclusions, providing theoretical evidence for realizing the rapid scanning technology of acoustic-optics landmine detection. Thus, the Bessel fringe database corresponding to the mode shapes of different mines in different environments can be established,and it is feasible to realize the rapid scanning technology of acoustic-optics landmine detection.