自旋飛行器角速率和過載測量數據擬合修正方法

嚴東升，童 偉，唐林卡，賈平會，郭 珂

（北京航天長征飛行器研究所，北京，100076）

自旋飛行器角速率和過載測量數據擬合修正方法

嚴東升，童 偉，唐林卡，賈平會，郭 珂

（北京航天長征飛行器研究所，北京，100076）

為正確分析飛行試驗飛行器迎角和氣動參數偏差等飛行動力學特性，根據具有非零慣性積旋轉飛行器在大氣層外的運動規律，推導飛行器角速率和過載分量的解析表達式。利用測量數據擬合估計相應的系數，給出測量數據系統偏差修正解算公式。計算實例表明，經擬合、系統偏差修正和過載質心換算的測量數據更接近實際飛行狀況，可用于確定飛行器大氣層初始段飛行迎角變化范圍和大氣層內氣動系數辨識等的飛行動力學特性結果分析。試驗數據結果分析實例驗證了該方法的有效性，研究結果對測量數據處理和誤差修正具有現實指導意義。

數據擬合；偏差修正；自旋飛行器

0 引 言

對飛行器角速率、過載參數的實測數據進行分析處理，是進行飛行試驗的飛行器飛行動力學特性結果分析計算的前提條件和重要依據[1]。由于測量傳感器本身存在制造工藝誤差[2]，電量與物理量之間轉換關系存在非線性影響，天地之間環境條件也存在差別，上述因素均會給測量數據帶來系統偏差和隨機偏差。這些偏差的存在直接影響飛行器飛行試驗動力學特性分析的可信度和準確性，嚴重時甚至可能得出不正確的結論。為此，一般需要在分析計算之前對所獲得的角速率、過載測量數據進行預處理，用擬合方法平滑這些參數消除隨機偏差，用已知的運動特性規律進行解算求出系統偏差，并對測量數據進行修正，以獲得更符合實際情況的測量參數數據。

本文以軸對稱旋轉飛行器為研究對象，推導給出了飛行器大氣層外角速率和過載分量的解析表達式。對飛行試驗過載和角速率測量數據進行擬合，修正了測量數據系統偏差，并將過載換算至質心處，獲得了更接近實際飛行狀況的試驗數據，為利用角速率分析大氣層初始段飛行器迎角范圍和利用大氣層內過載換算飛行器氣動特性參數奠定了數據基礎，試驗數據結果分析實例驗證了該方法的有效性。

1 真空段角速率變化特性

根據剛體動力學動量矩定理，飛行器繞質心轉動方程的矢量形式[3]為

式中 I為慣量張量矩陣； Jx， Jy， Jz為轉動慣量；Jxy， Jxz， Jyz為慣性積； 為飛行器相對慣性空間的旋轉角速度；ωx，ωy，ωz為在飛行器體軸坐標系各軸的角速率分量；t為飛行時間；為對質心合力矩的主向量。

絕大多數旋轉對稱體飛行器的慣性積 Jyz項均可以忽略不計，而慣性積 Jxy,Jxz不為零。因此可設Jy=Jz=J，ωx為中等大小，通常Jxy，Jxz，ωy，ωz，ω˙x，ω˙y，ω˙z也為足夠小量，這些小量的平方或小量的乘積都可以忽略不計。

大氣層外的真空飛行段，所有氣動力項均為零，僅存在陀螺效應，簡化上述轉動動力學方程，角速率微分方程可表示為[4]

解式（4）得：

式中 K為振動角頻率；A，B為待定常數；

角速率xω，yω，zω或可寫成式（8）所示的形式：

式中 為相位。

因此，式（8）在 ωz?ωy平面上的曲線形狀是一個圓，其圓心在(ωz=Dz，ωy=Dy)處，半徑為ωT。在真空飛行段，對于具有非零慣性積 Jxy和 Jxz的飛行器，其隨時間變化的測量記錄曲線將顯示出zω圍繞zzDω=作余弦振動，而yω圍繞著yyDω=作正弦振動，常數A，B，Tω可以從t=0給定的初始條件確定。

如果已精確地知道xω，yω和zω，也可以利用式（8）從飛行試驗測量數據求得飛行器慣性積 Jxy和 Jxz。

2 非質心處的附加過載

設飛行器非質心處任意一點 C，在飛行器體軸坐標系中的位置矢量r為

式中 x，y，z為位置分量；i，j，k為飛行器體軸坐標系三軸單位矢量。

由于飛行器的轉動運動，C點相對飛行器質心的附加加速度矢量a記為

式中xa，ya，za為附加加速度分量。

根據矢量導數定義，則有：

非質心處附加過載分量xnΔ，ynΔ，znΔ為：

式中0g為海平面引力加速度。

展開式（12）得到非質心處附加過載的分量形式[5]：

式中xω˙，yω˙，zω˙為角速率對時間導數。

在飛行器體軸坐標系中，安裝于非質心處（x，y，z)的過載傳感器測量得到的過載分量xn，yn和zn，它可以表示為質心處的過載分量xn′，yn′，zn′與附加過載分量xnΔ，ynΔ，znΔ之和：

將式（8）代入式（14），可得到真空段過載分量的解析表達式，如式（15）所示。

3 角速率和過載擬合系數估值及應用

根據最小二乘法原理，使上述目標函數最小，求得待估參數和相應的估值精度為

得到式（6）各參數擬合估計結果后，按式（7）從飛行試驗測量數據中求得xyJ和xzJ。利用式（8）可計算飛行器章動角、進動周期等真空段重要陀螺運動特征參數[5，6]。進一步利用式（6）和姿態角與角速率關系式可以解算飛行器真空段的姿態變化，結合外測數據從而得到進入大氣層初始段附近飛行器迎角及其變化范圍。

4 測量數據系統偏差及修正

式中 ωxi，ωyi，ωzi為對應ti時刻角速率測量結果；nxi，nyi， nzi為對應ti時刻過載分量測量結果。

大氣層內飛行段修正系統偏差后角速率、過載測量數據為

式中et為進入大氣層飛行時間。

消除測量數據系統偏差后，對于大氣層內飛行段過載數據，還需進行傳感器安裝位置修正，根據式（13）將其換算至質心處，用于進行進一步的大氣層飛行段動力學特性結果分析。

5 結果分析

采用本文方法對飛行器測量數據進行處理分析，得到滾轉角速率的擬合結果和俯仰、偏航角速率各系數的擬合值。圖1至圖3給出了測量角速率與擬合結果的對比曲線，兩者相吻合。采用擬合后的角速率表達式，可計算飛行器真空段飛行迎角α，從而得到大氣初始段附近的迎角變化范圍，具體曲線見圖4，符合小迎角飛行的調姿要求。

利用真空段角速率和過載擬合結果，對大氣層內過載進行了系統偏差修正，并換算至質心處。大氣層高空段，由于空氣稀薄，即使迎角較大，其橫、法向過載也應該很小。從圖5、圖6可以看出，經修正后的橫向過載zn、法向過載yn顯然更符合實際飛行情況。

6 結束語

根據飛行器真空段的陀螺運動效應，推導了角速率和過載的解析表達式，為測量數據的擬合分析提供了理論依據。利用飛行試驗測量數據擬合估計了解析式的相應系數，給出了角速率和過載測量數據系統偏差修正解算值。按給定的起始姿態和擬合角速率，計算飛行器真空段的姿態變化，得到了大氣層初始段飛行器迎角變化范圍，為飛行器軌道重建和飛行偏差分析提供初始依據。經過系統偏差修正和傳感器安裝位置影響修正的橫向、法向過載更接近實際飛行狀況，適合用于氣動系數換算等更深入的結果分析。

大氣層外俯仰、偏航角速率和3個過載通常數值較小，采用小量程傳感器并提高測量精度對數據處理分析至關重要。

[1] 王國雄. 彈頭技術[M]. 北京: 宇航出版社, 1993.

[2] 孟晨光等. 光纖通信技術在再入飛行器遙外測系統遠距離測控中的應用[J]. 導彈與航天運載技術, 2012年(增刊1): 84-86

[3] 錢杏芳, 林瑞雄, 趙亞男. 導彈飛行力學[M]. 北京: 北京理工大學出版社, 2011.

[4] 趙漢元. 飛行器再入動力學和制導[M]. 長沙: 國防科技大學出版社,1997.

[5] 賈沛然, 陳克俊, 何力. 遠程火箭彈道學[M]. 長沙: 國防科技大學出版社, 1993.

[6] 李曉軒, 汪群山. 軸對稱再入航天器質心橫偏對滾轉速率的影響[J]. 導彈與航天運載技術, 2012年(增刊1): 28-29.

Method of Angular Rate and Acceleration A Fitting and Revising Measurement Data for Spinning Vehicle

Yan Dong-sheng, Tong Wei, Tang Lin-ka, Jia Ping-hui, Guo Ke
(Beijing Institute of Space Long March Vehicle, Beijing, 100076)

In order to analyze the flight dynamic characteristics of hypersonic vehicles, including angle of attack and aerodynamic parameter error, the analytic expressions for angular rate and acceleration of hypersonic vehicles were deduced base on the motion law of exoatmospheric vehicle possessing nonzero product of inertia. The corresponding coefficients and deviation correction formulas of data measurement system were estimated using the measured data. The proposed method can treat the measured data effectively, and then the actual flight dynamic characteristics of hypersonic vehicles can be analyzed exactly based on the revised data. The computational result shows that the proposed fitting and revising method is effective and the research results are helpful on data processing and error correction.

Data fitting；Error correction；Spinning vehicle

V412.5

A

1004-7182(2017)05-0029-05

10.7654/j.issn.1004-7182.20170507

2017-06-26；

2017-08-16

嚴東升（1963-），男，研究員，主要研究方向為飛行器飛行動力學