盧玉斌 程永生 孫遠(yuǎn)程

(1.中國工程物理研究院電子工程研究所 四川綿陽 621900；2. 西南科技大學(xué)制造過程測試技術(shù)教育部重點實驗室 四川綿陽 621010)

?

侵徹混凝土類目標(biāo)過載數(shù)據(jù)的處理原則研究

盧玉斌1,2程永生1孫遠(yuǎn)程1

(1.中國工程物理研究院電子工程研究所 四川綿陽 621900；2. 西南科技大學(xué)制造過程測試技術(shù)教育部重點實驗室 四川綿陽 621010)

針對目前尚無統(tǒng)一、有效的侵徹過載數(shù)據(jù)處理方法這一問題，基于侵徹過載信號特征，提出處理過載信號時濾波截止頻率的選擇原則。該原則即按10/T(T為侵徹過載的持續(xù)時間，單位為s)計算得到的頻率，并結(jié)合實測過載的頻譜分析結(jié)果，確定濾波截止頻率獲取彈體的剛體侵徹過載。按照此原則處理實測過載數(shù)據(jù)，所得剛體過載的一次積分或二次積分分別與實測的撞擊速度和侵徹深度保持了較好的一致性，相對誤差小于10%，驗證了該原則的有效性。

信號處理 侵徹過載 混凝土類目標(biāo) 截止頻率

侵徹過載指彈體侵徹過程中產(chǎn)生的減加速度信號。研究侵徹過載特性的重要內(nèi)容之一是確定剛體過載，準(zhǔn)確的剛體過載為精確有效的計算侵徹深度、穿透目標(biāo)硬層數(shù)或空穴數(shù)等參數(shù)提供可靠的數(shù)據(jù)支持。彈體在侵徹靶板過程中，除了彈體自身的剛體過載外，還有作用于彈體上的各種復(fù)雜的振動信號，加上所疊加的其它信號。實測侵徹過載信號是主要包含兩種信號成分：一是在侵徹靶板過程中彈體遇到的阻力所形成的減加速度信號，由無數(shù)的小擾動疊加形成的脈沖信號；二是在侵徹過程中彈體產(chǎn)生的振動信號，包括彈體橫向振動和彈體縱向振動[1]。

引信所處動態(tài)環(huán)境的復(fù)雜性決定了侵徹過載信號是含有許多瞬變成分的非平穩(wěn)時變信號。同時在頻域上侵徹過程表現(xiàn)為具有無限帶寬的連續(xù)頻譜，因此在侵徹測量中應(yīng)盡量避免使用濾波器。但實際的侵徹過載往往表現(xiàn)為復(fù)雜的時域波形，這就需要采用低通濾波器濾除機械、電的干擾和噪聲以及由加速度計的諧振等引起的高頻噪聲。疊加在剛體過載信號上的這些噪聲信號致使無法準(zhǔn)確確定剛體過載。因此，在侵徹過載測量中，通常選用二階巴特沃斯型低通濾波器[2]。

但在濾除高頻噪聲的同時，濾波器也會造成輸出波形的畸變，以致可能影響所確定剛體過載的準(zhǔn)確性。且隨著濾波器截止頻率的降低，所確定的剛體過載的誤差會變大。為此，在分析實測侵徹數(shù)據(jù)特征的基礎(chǔ)上，本文提出選擇侵徹過載信號濾波截止頻率的原則，并結(jié)合實例驗證該原則的有效性。

1 侵徹混凝土類目標(biāo)過載信號特征分析

彈體侵徹混凝土類目標(biāo)是一個復(fù)雜的物理過程，由于彈體在侵徹過程中并非理想剛體，而是一個存在著復(fù)雜的應(yīng)力波傳播的結(jié)構(gòu)體，侵徹過載曲線全程呈振蕩型。侵徹通常包括開坑階段和成孔階段(即穩(wěn)定侵徹階段)。在開坑階段，由于彈體存在強烈的沖擊振蕩，在時域曲線中過載信號表現(xiàn)為前沿具有高幅值窄脈沖。對于抗壓強度較低的混凝土類目標(biāo)，若彈體能夠侵入或穿透目標(biāo)，其成孔階段的過載信號表現(xiàn)為恒定值。如圖1(a)的曲線是71 mm試驗彈垂直侵徹抗壓強度為18 MPa素混凝土靶的一發(fā)實測侵徹過載波形，彈速為569 m/s[3]；圖1(b)為彈體以632 m/s的速度垂直侵徹設(shè)計抗壓強度為35 MPa的鋼筋混凝土靶的實測侵徹過載曲線[4]。對于抗壓強度較高的混凝土類目標(biāo)，當(dāng)彈體的侵徹速度不足以使其侵入目標(biāo)內(nèi)部時，在成孔階段因靶板表面崩裂而不能形成穩(wěn)定的侵徹過程，在時域曲線中過載信號表現(xiàn)為逐漸衰減的趨勢。如圖2(a)所示為70 mm射彈以556 m/s的速度侵徹抗壓強度為48 MPa的素混凝土圓筒靶的侵徹過載曲線[5]；圖2(b)為彈體以1043.9 m/s的速度垂直侵徹抗壓強度為38 MPa的素混凝土靶的實測過載曲線[6]。

圖1 彈體侵徹和混凝土靶實測侵徹過載曲線Fig.1 Penetration overload curves of 18 MPa plain concrete and 35 MPa reinforced concrete

圖2 彈體侵徹和素混凝土靶實測侵徹過載曲線Fig.2 Penetration overload curves of 48 MPa and 38 Mpa plain concretes

通常，在撞擊瞬間會產(chǎn)生很強的應(yīng)力波，彈體內(nèi)各點的加速度實際上是應(yīng)力波在彈體內(nèi)傳播的綜合表現(xiàn)。在撞擊初期，應(yīng)力波將導(dǎo)致彈體各點加速度的第一個峰值；而應(yīng)力波在侵徹過程中的影響主要是在過載曲線產(chǎn)生振蕩。對于加速度傳感器，因其固有頻率有限，通常很難真實響應(yīng)彈體內(nèi)由應(yīng)力波所產(chǎn)生的第一個加速度峰值。

對圖1(a)和圖2(a)的信號進(jìn)行頻譜分析，結(jié)果見圖3?？梢钥闯鰪楏w侵徹混凝土類目標(biāo)的過載信號的頻譜主瓣寬度在5 kHz以內(nèi)，即彈體的質(zhì)心(剛體)加速度的頻譜分量主要集中在5 kHz以內(nèi)，主瓣后的第一個頻譜峰值對應(yīng)于彈體的一階振動(軸向振動)，該頻率分量主要取決于彈體長度[5]。

圖3 圖1(a)和圖2(a)所示實測侵徹過載頻譜分析Fig.3 Spectral analysis of Fig.1(a) and Fig.2(a)

因此，對于混凝土類目標(biāo)，侵徹過載信號的基本特征為在時域上表現(xiàn)為單峰性，且在200～300 μs量級[7]的高幅值脈沖后持續(xù)ms量級的低幅值加速度，彈體的剛體過載頻譜分量主要集中在5 kHz以內(nèi)，結(jié)構(gòu)響應(yīng)的頻率分量主要集中在5～10 kHz[5]。

2 選擇侵徹過載信號濾波截止頻率的原則

在侵徹撞擊的初始時刻，由于高幅值窄脈沖的作用，彈體的高階模態(tài)被激發(fā)，在過載信號中表現(xiàn)為較高的頻率分量。而在研究侵徹彈體強度、設(shè)計侵徹引信和反侵徹防護(hù)時，有用信號是剛體過載。要獲得彈體的剛體過載，必須對實測過載數(shù)據(jù)進(jìn)行處理，而處理的關(guān)鍵是如何確定數(shù)字濾波器的截止頻率。

對實測侵徹過載信號進(jìn)行頻譜分析(如圖3)發(fā)現(xiàn)，信號的頻譜分量存在許多頻譜峰值點。從理論上來說，這些峰值點應(yīng)該與彈體軸向振動相關(guān)的模態(tài)頻率相對應(yīng)，但實彈侵徹過程中彈體的約束條件比采用有限元軟件進(jìn)行模態(tài)分析時施加的邊界約束條件更復(fù)雜，所以實際的頻譜峰值點與各階模態(tài)頻率之間有一定的偏差(如文獻(xiàn)[3]和文獻(xiàn)[5])。

在信號與系統(tǒng)和數(shù)字信號處理原則中，周期信號的幅度譜是線譜，而非周期信號的幅度譜為連續(xù)譜，信號所占頻率是一個頻率范圍，即具有一定的頻率寬度。因此，截止頻率一般可取為 10/T(T為侵徹過載的持續(xù)時間，單位為s) Hz[2, 8]。如果波形中疊加的高頻噪聲過于嚴(yán)重，為有效濾掉非剛體過載信號，可適當(dāng)降低濾波截止頻率，選取一階模態(tài)頻寬的起點，用這個頻率點進(jìn)行低通濾波可將非剛體過載的高頻信號濾掉，同時滿足沖擊加速度測量的不確定度要求。

因此，本文提出侵徹過載信號濾波截止頻率的選擇原則為在參照10/T的截止頻率的基礎(chǔ)上，結(jié)合實測過載信號的一階模態(tài)頻寬起點，保證濾波前后過載信號的積分速度在一定的誤差范圍內(nèi)。根據(jù)此原則確定的濾波器截止頻率通常稱為臨界截止頻率，按此頻率進(jìn)行濾波獲得的加速度即為彈體的剛體過載。

3 實測過載信號的數(shù)據(jù)處理

圖2(a)所示曲線為988型高g值壓電加速度傳感器測得的70 mm射彈以556 m/s的速度侵徹抗壓強度為48 MPa的素混凝土圓筒靶的過載、經(jīng)修正零點漂移后的曲線[5]。從圖中可知，T=2.8 ms，所以10/T為3.6 kHz。圖3(a)為圖2(a)所示過載曲線的頻譜。從圖中可以看出一階模態(tài)頻寬起點為2.5 kHz。對過載曲線以截止頻率2.5 kHz進(jìn)行濾波，濾波后得到的剛體過載曲線如圖4所示。剛體過載曲線經(jīng)一次積分和二次積分得到的速度和位移分別為566 m/s和0.563 m，與實測得到的速度556 m/s和侵徹深度0.601 m分別相差1.8%和6.3%。

圖4 2.5 kHz濾波后的過載曲線Fig.4 Overload curve after 2.5 kHz filtering

經(jīng)數(shù)值計算可知，當(dāng)截止頻率小于2.5 kHz時，積分后得到的速度會隨著截止頻率的降低而減小。

圖1(a)所示曲線是直徑71 mm的試驗彈以569 m/s的彈速侵徹抗壓強度為18 MPa素混凝土靶的一發(fā)實測侵徹過載波形[3]，由圖可知10/T為2.5 kHz。圖5是它的頻譜分析中2～5 kHz范圍的展開圖，可見一階模態(tài)頻寬起點處對應(yīng)的頻率為3.6 kHz，用此頻率進(jìn)行低通濾波可將非剛體高頻信號全部濾掉。 圖6(a)為截卵形頭部空心彈體以135 m/s速度侵徹混凝土靶的過載曲線[8]，由圖可知其脈沖持續(xù)時間為2.4 ms，得到10/T為4.2 kHz。圖6(b)為選取相應(yīng)的截止頻率得到的濾波后的過載曲線，積分得到的長桿彈侵徹混凝土前的速度與實測速度之間的誤差為6.5%。

圖5 頻譜分析局部展開圖Fig.5 Frequency of single-order modal analysis

圖6 侵徹混凝土過載時間曲線Fig.6 Overload time curves of penetration concrete before and after filtering

圖7為直徑71 mm的彈體以338 m/s的速度侵徹抗壓強度為38.6 MPa的疊層混凝土靶的過載響應(yīng)曲線[9]。由圖可知其脈沖持續(xù)時間為5 ms，得到10/T為2 kHz。選取相應(yīng)的截止頻率得到濾波后的過載曲線也示于圖7，積分得到侵徹混凝土前的速度與實測速度之間的誤差為2.66%。

圖7 侵徹混凝土疊層靶原始和濾波后的過載時間曲線Fig.7 Overload curves of penetration concrete laminated plate before and after filtering

圖8 侵徹混凝土靶過載曲線Fig.8 Overload curves in 3 states

圖8(a)所示曲線為B&K公司的高g值加速度傳感器8309測得的某型號彈以566 m/s的速度侵徹混凝土靶的過載、經(jīng)修正零點漂移后的曲線[10]。從圖中可知，T=2.5 ms，所以10/T為4 kHz。圖8(b)為侵徹過載的頻譜曲線，可見一階模態(tài)頻寬起點為3.2 kHz。對過載曲線以截止頻率3.2 kHz進(jìn)行濾波，濾波后得到的剛體過載曲線如圖8(c)所示。比較濾波后的信號及其積分與空腔膨脹理論預(yù)測值，二者基本吻合。

4 結(jié)論

本文提出的選擇濾波截止頻率的原則為：在參照10/THz的截止頻率的基礎(chǔ)上，結(jié)合實測過載信號的一階模態(tài)頻寬起點，滿足侵徹過載測量的不確定度要求，保證濾波前后過載信號的一次或二次積分值不變或在一定的誤差范圍內(nèi)。數(shù)據(jù)處理結(jié)果表明此原則是有效的，經(jīng)處理后得到的剛體過載信號完整地反映了彈體侵徹混凝土類目標(biāo)過程中過載的變化規(guī)律，彈體剛體過載的一次積分或二次積分與相應(yīng)實測結(jié)果吻合較好，相對誤差均小于10%，滿足侵徹測試的工程要求。

[1] 郝慧艷, 李曉峰, 劉明杰. 侵徹靶板過程中剛體過載特征提取方法[J]. 火力與指揮控制, 2012, 37(12):42-44, 47.

[2] 曹亦慶, 李善明, 王卓, 等. 沖擊加速度峰值測量中低通濾波器截止頻率的影響[J]. 科技導(dǎo)報, 2008, 26(16):81-84.

[3] 黃家蓉, 劉瑞朝, 何翔, 等. 侵徹過載測試信號的數(shù)據(jù)處理方法[J]. 爆炸與沖擊, 2009, 29(5):555-560.

[4] 周寧, 任輝啟, 沈兆武, 等. 侵徹鋼筋混凝土過程中彈丸過載特性的實驗研究[J]. 實驗力學(xué), 2006, 21(5): 572-578.

[5] 范錦彪, 祖靜, 徐鵬, 等. 彈丸侵徹混凝土目標(biāo)減加速度信號的處理原則[J]. 探測與控制學(xué)報, 2012, 34(4):1-5, 9.

[6] FRANCO R J, PLATZBECKER M R. Miniature penetrator (MINPEN) acceleration recorder development test[R]. Sandia National Laboratories, 1998.

[7] 張建新. 侵徹引信炸點控制理論及試驗研究[D]. 南京:南京理工大學(xué), 2012.

[8] 王華軍, 王燕, 李新娥. 空心彈侵徹混凝土過載特性的實驗研究[J]. 測試技術(shù)學(xué)報, 2012, 26(2):150-153.

[9] BOOKER P M, CARGILE J D, KISTLER B L, et al. Investigation on the response of segmented concrete targets to projectile impacts[J].International Journal of Impact Engineering,2009,36(7):926-939.

[10] 徐鵬, 祖靜, 范錦彪. 高速動能彈侵徹硬目標(biāo)加速度測試技術(shù)研究[J]. 振動與沖擊, 2007, 26(11):118-122.

Investigation on Penetration Deceleration Signal Processing Method for Penetrating Targets of Concrete-like Materials

LU Yu-`bin1,2, CHENG Yong-sheng1, SUN Yuan-cheng1

(1.InstituteofElectronicEngineering,CAEP,Mianyang621900,Sichuan,China;2.KeyLaboratoryofTestingTechnologyforManufacturingProcess,MinistryofEducation,SouthwestUniversityofScienceandTechnology,Mianyang621010,Sichuan,China)

Until now, there is no standard and effective processing method for the penetration deceleration data. Based on the analysis of the characteristic of the deceleration signal during the projectile penetration into targets of concrete-like materials, the selection principle of the cut-off frequency is put forward. This principle is combining the cut-off frequency of 10/T (where T is the duration of penetration deceleration, with the unit of s) with the spectrum analysis results of the measured deceleration data, to determine the cut-off frequency of the digital low-pass filter used for extracting rigid-body deceleration. According to the above principle, the penetration deceleration data are processed and the rigid-body deceleration data are obtained. The single or double integration of rigid-body deceleration is consistent with the measured penetrating velocity or penetration depth and the relative error is less than 10%.

Signal processing; Penetration deceleration; Concrete-like target; Cut-off frequency

2014-10-28

NSAF聯(lián)合基金(U1430110)；中國博士后科學(xué)基金(2012M511943)。

盧玉斌(1980—)，男，博士，研究方向為沖擊動力學(xué)。E-mail: luyubin@swust.edu.cn

O39

A

1671-8755(2015)01-0075-05