張孝云,江小華
(南京理工大學機械工程學院,南京 210094)
航空子母彈藥在打擊敵方縱深大面積多目標領(lǐng)域中有著明顯優(yōu)勢。為提高子彈破片殺傷效果,子彈引信須嚴格控制空炸高度[1-2]。考慮子彈引信為脈沖激光定高引信的情況,子彈在空中穩(wěn)速下落時,受氣流影響一般呈傾斜降落姿態(tài),導致激光脈沖光路與目標地面發(fā)生傾斜。由于激光具有高定向性的特點,激光光路傾斜時,在同等高度上引起激光引信的光程加長,導致實際空炸高度低于所設(shè)置的最佳空炸高度。
為修正因彈體傾斜所帶來的空炸高度偏差,文中提出利用彈體傾角修正空炸高度的方法。文獻[3]應用液體表面在重力作用下平衡時可形成絕對水平面的原理,利用激光束在液體表面或漂浮物表面上一次或多次反射折射實現(xiàn)光學位移放大,雖然能夠?qū)崿F(xiàn)高精度傾角檢測,但測量裝置相對于引信狹小空間而言,體積過大、結(jié)構(gòu)復雜,很難應用于常規(guī)彈藥引信中。文中采用具有微小體積的慣性力加速度傳感芯片,利用重力在傳感軸方向的分量與傳感軸傾角的關(guān)系來進行定高修正。
在器件的選擇上要考慮功能、體積、功耗、穩(wěn)定性、集成度、成本等方面的要求。引信的空間小,要求器件及其外圍的電路連接盡量簡單。選擇功耗低、集成度高、功能強、操作控制簡單且成本適中的MEMS慣性力加速度傳感器芯片是較為可行的。
經(jīng)比較,選用意法半導體公司新一代MEMS數(shù)字式三軸線性加速度傳感器LIS331DL。該傳感芯片采用LGA封裝,3mm×3mm×1 mm的微小體積,完全適應引信狹小的空間;抗過載能力強,可承受高達10000g的沖擊;量程可編程選擇為 ±2/±8g;內(nèi)置ADC采樣頻率可編程為100/400Hz,以適應微控器讀取速率;有I2C/SPI兩種外部總線接口,方便與微控器的連接。
設(shè)脈沖激光定高引信測得距目標地面的斜距離值為L,彈體實際高度為H,彈軸傾角為θ,則:
慣性力加速度傳感器按傳感軸的數(shù)量分為單軸、雙軸和三軸。以地面慣性參考系為定系,彈體坐標系為動系。單軸修正模式下,令Z軸沿彈軸方向(激光定高引信光軸方向)布置,θ為彈軸偏離重力方向的角度,如圖1所示。
重力加速度值記為g,傳感軸Z的理論輸出記為Az(θ)。
彈軸傾角的余弦值為定高修正值,記:
靈敏度:
圖1 單軸修正模式示意圖
LIS331DL的傳感軸輸出值寄存器是8位的,輸出數(shù)據(jù)有正負之分,表示方式為二進制補碼。在1g作用下,Z軸最大輸出絕對值是58,可以認為58代表1g輸入時,傳感軸對應的輸出值,對應分辨率為0.017g。
在單軸修正模式下,Z軸與彈軸方向一致,Z軸傾角代表彈軸傾角,以Z軸輸出值與0.017的乘積作為定高修正值。由式(4)可見,在0°~90°范圍內(nèi),傳感軸輸出值對彈軸傾角的靈敏度是逐漸增大的,在0°~20°范圍內(nèi),彈軸傾角平均每變化5.35°,傳感軸寄存器才有1個單位的變化,而在70°~90°范圍內(nèi)時,傾角平均每變化1.01°,就可以引起傳感軸寄存器1個單位的變化。
由式(1)知,當彈軸傾角較大時,不修正會引起很大的測量誤差;由式(4)知,單軸修正模式適合彈軸傾角較大的情況。
雙軸修正模式下,Z軸仍沿彈軸方向布置,X、Y 軸所構(gòu)成的平面P垂直Z軸方向,如圖2所示。
圖2 雙軸修正模式示意圖
X、Y兩軸偏離重力方向的角度分別記為 α、β ;理論輸出值分別記為 Ax(θ)、Ay(θ)。
地面定坐標系記為Ouvr,彈體坐標系記為Oxyz,從Ouvr變換到Oxyz,可以認為坐標系Ouvr先繞v軸傾斜θu,再繞u軸傾斜θv,根據(jù)線性空間變換理論可知:
進一步計算得:
重力加速度在坐標系Ouvr中的向量表示記為:(ugvgrg)= ( 0 0g),由式(8)坐標系變換關(guān)系可得重力加速度在坐標系Oxyz中的向量表示為:
可以驗證向量(xgygzg)的模正好是g,其中:
可以證明平面P相對于水平面的傾角δ滿足關(guān)系式:
由式(5)~式(6)、式(9)~式(12)得:
彈軸Z方向垂直平面P,由幾何關(guān)系知,平面P相對于水平面的傾角δ大小等于彈軸Z偏離重力方向的傾角θ。
與單軸修正模式不同,雙軸修正模式是通過X、Y軸偏離重力方向的角度來間接反映彈軸的傾角。由式(13)可見,定高修正值與X、Y兩軸的輸出值之間是非線性關(guān)系。當彈軸傾角較小時,X、Y兩軸的輸出值較小,對傾角的靈敏度高;當彈軸傾角較大時,X、Y雙軸的輸出值較大,對傾角的靈敏度低。隨著彈軸傾角的增大,平均每1°的變化所帶來的計算誤差在逐漸增大,在0°~20°范圍內(nèi)平均每1°誤差所對應的定高修正誤差是0.319%,在20°~45°范圍內(nèi),該值升到1.533%,約是前者的5倍。因此,為了減小定高修正誤差,較大傾角時,采用對大傾角靈敏度高的單軸修正模式;較小傾角時,采用對小傾角靈敏度高的雙軸修正模式。
LIS331DL的3個傳感軸在結(jié)構(gòu)上相互垂直的特點,恰好滿足小傾角和大傾角范圍內(nèi)對不同修正模式的要求。具體布置時,Z軸沿彈軸方向,X、Y軸構(gòu)成的平面P與彈軸方向垂直。
當子彈處于傾斜平衡狀態(tài)時,文中所分析的修正模式可以較好的修正因子彈傾斜所造成的定高誤差。當子彈出現(xiàn)加速或減速運動時,加速度傳感器3個軸的矢量和不再是1g,在計算修正值前,先計算出3個軸矢量和的模,并以此模作為修正基準,這樣可有效減小子彈非平衡狀態(tài)下的修正值誤差。例如當子彈加速運動時,可減小過高的修正值;減速運動時,可提高過低的修正值。
此外,LIS331DL進行A/D轉(zhuǎn)換時存在一定的轉(zhuǎn)換誤差,由于溫度變化造成芯片零漂移也會帶來一定的誤差,假定傳感軸輸出值由確定性成分f(t)和隨機起伏的誤差成分e(t)組成。離散化采樣時,表達形式可寫成:
對非平穩(wěn)數(shù)據(jù){yk},在適當小的區(qū)間上視為接近平穩(wěn)而做局部平均,以減小{ek}所造成的隨機起伏。沿全長N個數(shù)據(jù)的逐一小區(qū)間上不斷的局部平均,最終可使傳感軸輸出結(jié)果{fk}較平滑,濾掉頻繁起伏的隨機誤差[4]。算式為:
式中:ωi為權(quán)系數(shù)為小于m的任一正整數(shù),且p+q+1=m。
取m=7,p=q=3,即所采用的滑動平均法中的等權(quán)中心平均法,LIS331DLAD采樣頻率設(shè)置為400Hz。
在ZTS-3三軸轉(zhuǎn)臺上,對所設(shè)計的定高修正模塊進行單、雙軸修正值輸出測試。以模塊傾角模擬彈軸傾角。Z軸輸出值與0.017的乘積為單軸模式下的定高修正值;調(diào)用式(13)可根據(jù)X、Y軸的輸出值計算出雙軸模式下的定高修正值。
模塊傾角與轉(zhuǎn)臺傾角之間由于模塊實際安裝過程中不能完全對平的影響存在固有角度偏差。為確定此角度偏差,令轉(zhuǎn)臺處于 5°、10°、15°、20°、25°、30°共6個不同的狀態(tài),由于雙軸修正模式對此范圍內(nèi)的傾角變化靈敏度高,利用式(13)求出對應這6個狀態(tài)的模塊傾角,試驗時輸出計算值分別為:11.60°、16.27°、21.03°、25.86°、30.46°、35.58°,按照誤差平方和最小原則,計算出模塊相對于轉(zhuǎn)臺的理論固有偏角近似為6°。
以6°作為模塊相對于轉(zhuǎn)臺的固有安裝偏角,ZTS-3 三軸轉(zhuǎn)臺的角度誤差是0.02°,轉(zhuǎn)臺在 0°~80°范圍內(nèi),每5°進行一次測試。測試數(shù)據(jù)如表1所列。
表1 試驗測試數(shù)據(jù)表
表1中:φ為轉(zhuǎn)臺傾角(°);θ為模塊傾角(°);T_val為對應模塊傾角的理論定高修正值;D_val為實際測試給出的雙軸模式修正值;S_val為實際測試給出的單軸模式修正值;D_er為雙軸模式修正值相對修正誤差(%);S_er為單軸模式修正值相對修正誤差(%)。
從表1中可以看出:
1)隨著模塊傾角的增大,雙軸修正模式下,定高修正值相對誤差D_er逐漸增大,61°以下,誤差不超過10%;單軸修正模式下,其定高修正值相對誤差S_er起伏變化不大,在0~75°范圍內(nèi),維持在9% ~12%之間。
2)0~61°范圍內(nèi),雙軸模式修正精度要高于單軸模式,傾角越小,對比越明顯,61°以上時,情況相反。
根據(jù)以上結(jié)論,實際應用中,先計算雙軸模式下定高修正值對應的傾角計算值,后與61°比較,比較結(jié)果作為定高修正模式選擇的依據(jù),小于等于61°選用雙軸修正模式;大于61°選用單軸修正模式。
文中提出一種利用彈體傾角對激光定高引信進行定高修正的方法。建立了單、雙軸兩種修正模式,為減小定高修正值相對誤差,通過對兩種模式特點的分析,得出小傾角時選擇雙軸模式、大傾角時選擇單軸模式的結(jié)論,實驗結(jié)果表明,分段采用單、雙軸兩種修正模式的修正方法,可以校準傳感器與彈體軸的安裝角度誤差,并能夠以較高的精度修正因彈體傾斜所造成的激光引信定高誤差。
[1]袁平,張合,陳炳林.實現(xiàn)軟殺傷彈藥指令引信技術(shù)的途徑[J].南京理工大學學報,2004,12(6):571-574.
[2]崔占忠,宋世和,徐立新.近炸引信原理[M].3版.北京:北京理工大學出版社,2009.
[3]樂開端,曹建安.高精度激光傾角測量技術(shù)研究[J].光子學報,2004,33(9):1108-1110.
[4]林洪樺.動態(tài)測試數(shù)據(jù)處理[M].北京:北京理工大學出版社,1995.