變質心高超聲速再入飛行器建模與特性分析

樊博璇 陳桂明 林洪濤

變質心高超聲速再入飛行器建模與特性分析

樊博璇1陳桂明1林洪濤2

（1. 火箭軍工程大學，西安 710025；2. 中國人民解放軍31004部隊，北京 100091）

建立了一維變質心滾控高超聲速再入飛行器的完整動力學模型。對變質心飛行器彈道及姿態動力學進行研究，明確了變質心滾控飛行器的控制機理與特點，同時針對其通道間的耦合效應以及執行機構參數設置對控制能力的影響分析，提出了飛行器結構設計和控制系統設計的相關要求，為高超聲速再入飛行器的工程實踐提供理論基礎。

高超聲速再入飛行器；變質心控制；動力學

1 引言

變質心技術是機動變軌突防的一種新形式，相比其他傳統控制技術（如空氣舵控制、發動機控制等）具有低成本、高效率等優勢，可有效提高高速飛行器的突防能力[1～3]。

美國是最早開展再入飛行器變質心控制研究的國家之一[4～6]。1997年，俄羅斯研制的白楊M導彈也采用變質心控制技術成功實現了彈道修正[7]。國外對變質心控制已有較為深入的研究，并進入工程實踐階段。國內對變質心控制技術的研究起步較晚，西北工業大學是較早涉及這方面研究的單位，易彥等[8]利用拉格朗日法對變質心旋轉飛行器建立了六自由度動力學方程，并對活動質量體的運動所帶來的姿態角變化進行了仿真。周軍等[9]研究了三維變質心非交互控制的可解性問題，給出三種變質心非交互控制實現性類別及分類判斷方法。陳升澤等[10]采用反演設計方法，引入虛擬控制量，設計了一種針對雙活動質量塊變質心飛行器的滑模變結構姿態控制器。高長生[11]利用 Newton-Euler 法建立了攜帶三個活動質量體的飛行器動力學模型，并在滾動穩定作為約束條件下對俯仰-偏航通道動力學進行了分析。劉智陶等[12]針對單質量塊滾控式變質心飛行器的欠驅動問題，設計了一種自抗擾控制器進行滾偏耦合控制。目前國內關于變質心控制技術仍處于理論研究階段，主要集中于多質量塊的配平控制分析及控制器設計，而對一維單質量塊滾控式變質心飛行器的機理分析、控制性能問題的研究較少，對于高超聲速再入飛行器變質心控制的完整建模和系統分析缺乏系統研究。

本文以一維變質心滾控高超聲速再入飛行器為研究對象，基于牛頓力學建立了包含質量塊運動模型的完整動力學模型；然后通過理論分析與仿真，研究變質心高超聲速再入飛行器的飛行路徑及姿態動力學，明確變質心控制機理與特點；最后結合飛行器質心運動方程，分析了各方案下飛行器機動能力及變質心機構控制能力，為高超聲速再入飛行器變質心控制機理的分析與制導控制系統的設計提供一定的理論基礎。

2 飛行器動力學建模

本文研究對象如圖1所示，整個飛行器系統主要包括飛行器殼體和一套變質心機構，變質心機構由可移動質量塊（實線部分）和不可移動部分（虛線部分）組成，O為飛行器殼體質心，O為可移動質量塊質心。飛行器采用固定配平攻角外形，通過單軸滑動飛行器內部質量塊形成的氣動力矩調節飛行器姿態角，使飛行器按一定的彈道軌跡運動。

圖1 一維變質心滾控飛行器示意圖

飛行器建模時所用的坐標系及符號定義如下：

在建立一維變質心再入飛行器動力學模型時，只考慮再入段飛行器質點系所受的重力、空氣動力和空氣動力矩。

2.1 殼體質心平動動力學方程

在慣性系下，飛行器各部分質心矢徑和飛行器系統總質心矢徑有關系式：

則式（3）可改寫為：

式（6）在彈道系中表示為：

a. 相對加速度

式（7）等號左邊可表示為：

b. 地球引力

再入系中，根據式（4），彈道系下的地球引力可表達為：

c. 空氣動力

d. 牽連加速度

e. 哥氏加速度

f. 變質心機構引起的附加慣性力

將以上各項代入式（7），可得到彈道系下一維變質心飛行器的質心動力學方程：

2.2 繞殼體質心轉動動力學方程

根據相對動量矩定理可得：

a. 殼體相對動量矩對時間的導數

b. 質量塊相對動量矩對時間的導數

c. 不可動部分相對動量矩對時間的導數

d. 合外力對殼體質心的力矩

e. 牽連慣性力矩

牽連慣性力矩可表示為：

將上述各項代入式（16），可得：

設質量塊和不可動部分初始安裝位置質心重合，根據矢量求導法則和矢量叉乘運算法則，將力和力矩表達式代入式（22），化簡可得繞殼體質心轉動動力學方程的矢量形式為：

3 飛行器特性分析

由理論力學可知，當質量塊相對殼體運動時，所產生的慣性力會影響質心運動。由于變質心機構質心的偏移，不僅使飛行器轉動慣量發生變化，也引入了慣性力矩的干擾，從而影響飛行器的姿態運動。因此，一維變質心飛行器具有更加復雜的動力學方程，有必要分析其模型特性。

3.1 飛行器動力學分析

3.1.1 軌道動力學分析

一維變質心再入飛行器所受外力包括氣動力、地球引力以及由變質心機構運動引入的慣性力。

a. 氣動力

一維變質心再入飛行器所受氣動力大小與飛行器動壓及氣動系數成正比。由于高空段空氣密度較小，因此飛行器動壓較小。同時，飛行器俯仰與偏航通道依靠自身靜穩定度保持穩定，因此其氣動系數也較小，此時飛行器所受的空氣動力不大。隨著再入段飛行器高度的高速下降，空氣密度迅速增大，動壓與空氣動力也急劇增長，成為飛行器運動的最主要影響因素。

圖3 飛行器氣動力與各向過載變化曲線

對飛行器法向力和軸向力進行仿真，仿真曲線如圖3所示。可以看到，隨著高度下降，飛行器動壓增大，各向空氣動力隨之增大。在40km以上的高空，軸向力與法向力最大值約3000N，隨著高度下降，各向氣動力增加，軸向力與法向力最大量級為105。

b. 地球引力

c. 變質心機構慣性力

通過以上分析可知，對于一維變質心再入飛行器的彈道，在高空段主要受地球引力影響，在40km以下，主要受氣動力影響。而變質心機構慣性力則對彈道影響較小，可通過力量級比較對彈道運動適當簡化。

3.1.2 姿態動力學分析

由動力學模型可知，飛行器姿態運動中所受力矩主要包括氣動力矩、陀螺力矩、變質心機構控制力矩以及慣性力矩。此外，由于飛行器系統轉動慣量與慣性主軸的改變，還將引入各通道姿態運動的耦合。下面逐一進行仿真分析。

a. 氣動力矩

圖4 飛行器滾轉力矩與俯仰力矩隨高度變化曲線

飛行器氣動力矩大小主要取決于飛行器氣動角和動壓。以飛行器殼體質心為參考點，從圖4的仿真曲線可以看到，隨著飛行器高度下降，飛行器俯仰力矩可達103量級，因此氣動力矩是影響飛行器姿態的主要因素。

b. 陀螺力矩

c. 控制力矩

可以看出，控制力矩包括兩部分，一是由質量塊初始安裝位置引起的附加氣動力矩，可對飛行器的氣動力矩進行平衡；二是由質量塊主動移動引起的控制力矩，可對飛行器姿態進行主動調節。假設變質心機構初始位置，結合飛行器飛行彈道，仿真得到附加力矩曲線如圖5所示。

通過曲線可知，附加氣動力矩是影響飛行器姿態的主要因素之一，對飛行器起著姿態控制作用，且隨著動壓的增大而增大。由于控制耦合的存在，質量塊運動會為偏航通道引入附加氣動力矩，其量級與滾轉通道控制力矩相當。由于飛行器偏航通道轉動慣量約為滾轉通道轉動慣量的10倍，且在實際飛行控制中，質量塊并非持續做大幅運動，因此偏航通道在側向靜穩定性下可保證振蕩收斂。

d.慣性力矩

質量塊的運動將導致飛行器受到變質心機構慣性力的作用，當變質心機構與飛行器質心存在位置偏差時，慣性力將以慣性力矩的形式被引入飛行器繞質心運動中，慣性力矩表達式為：

e.慣量耦合

由于質量塊運動將導致飛行器總質心偏移，飛行器總轉動慣量將發生變化，而轉動慣量的變化將會進一步加強各通道間的運動耦合與控制耦合。質量塊運動對系統轉動慣量的影響主要體現在兩個方面，一是改變了慣性主軸轉動慣量大小；二是改變了慣性主軸位置，產生了慣量積。

由此可見，質量塊的運動使主慣量軸上的轉動慣量增大，各慣量軸上變化范圍不超過5%。質量塊運動主要增大了滾轉與偏航通道以及滾轉與俯仰通道間的慣量積，而對偏航與俯仰通道的影響較小。

3.2 飛行器控制能力分析

由3.1.2c分析可知，控制力矩表達式的第一部分為機構安裝位置引入的附加氣動力矩，主要用來平衡飛行器氣動力矩；第二部分為控制力矩，用來主動調節飛行器姿態。控制力矩的大小直接決定了飛行器滾控響應的快速性。

結合飛行器姿態方程，俯仰與偏航通道仍采用“瞬時平衡”假設，可得控制力矩作用下飛行器滾轉角加速度表達式為：

將上式作為一維變質心滾控飛行器控制能力的衡量指標，對其進行仿真分析，仿真曲線如圖6所示。

仿真曲線表明，由于40km以上的動壓較小，飛行器最大滾動控制力矩較小。但隨著高度下降，動壓增大，飛行器最大滾控力矩與最大滾轉角加速度明顯增大，其量級為10。因此，變質心控制的優越性主要體現在低空（40km以下）段。

4 結束語

a. 變質心控制機構處于飛行器內部，因而無需特殊解決控制機構的燒蝕問題，同時也減少了某些控制參數的時變性；

b. 飛行器俯仰與偏航姿態均為振蕩收斂狀態，其收斂點為各通道配平角，對于側滑即為零側滑角；

c. 飛行器彈道主要受氣動力與地球引力影響，變質心機構運動產生的慣性力對其影響較小；

d. 隨著高度下降，飛行器最大滾控力矩與最大滾轉角加速度明顯增大。因此，變質心控制的優越性主要體現在低空段；

e. 利用高速氣動力和氣動力矩進行姿態控制，只要設計合理，在同樣姿態控制要求下，可最大限度地節約能源，提高飛行器有效載荷的攜帶能力；

f. 只要飛行器的質心有很小的變化，就能改變飛行器靜穩定度，產生足夠的過載改變飛行器的飛行軌跡，為控制系統提供了一種方便有效的作動方式，大大提高了飛行器的機動能力。

1 Petsopoulo T, Regan F J. A moving-mass roll control system for a fixed-trim re-entry vehicle [C]. The 32nd Aerospace Science Meeting and Exhibit, Reno, US, January, 1994: 10～13

2 Petsopoulos T, Regan F J. Moving-mass roll control system for fixed-trim re-entry vehicle [J]. Journal of Spacecraft and Rockets, 1996, 33(1): 54～61

3 高長生，李君龍，荊武興，等. 導彈質量矩控制技術發展綜述[J]. 宇航學報，2010，31(2)：307～314

4 Anon. Generic aerocapture atmospheric entry study, volume 1 (final report) [R]. US: General Electric Co., October 1980

5 Regan F J, Kavetsky R A. Add-on controller for ballistic reentry vehicle [J]. IEEE Transactions on Automatic Control, 1984, 12(6): 869～880

6 White J E, Robinett R D. Principal axis misalignment control for deconing of spinning spacecraft [J]. Journal of Guidance, Control and Dynamics, 1994, 17(4): 823～830

7 Norris R, Kristensen H. Russian nuclear force[J]. Bulletin of the Atomic Scientists, 2005, 61(2): 70～72

8 易彥，周鳳岐，周軍. 基于變質心控制導彈的運動分析[J]. 航天控制，2000(3)：1～5

9 周軍，王霄婷，林鵬. 三維變質心非交互控制的可實現性分析[J]. 宇航學報，2013，34(3)：1000～1328

10 陳升澤，張旋，解春雷，等. 基于反演設計的變質心飛行器滑模變結構控制器設計[J]. 導彈與航天運載技術，2018，362(4)：65～69

11 高長生，陳爾康，荊武興. 高超聲速飛行器機動規避軌跡優[J]. 哈爾濱工業大學學報，2017，49(4)：16～21

12 劉智陶，李澗青，高長生. 欠驅動變質心飛行器的滾偏耦合自抗擾控制[J/OL]. 北京航空航天大學學報：1-10[2020-11-11]. https://doi.org/10.13700/j.bh.1001-5965.2020.0206.

Modeling and Characteristic Analysis of Hypersonic Reentry Vehicle with Moving Mass Control

Fan Boxuan1Chen Guiming1Lin Hongtao2

(1. Rocket Force University of Engineering, Xi’an 710025；2. Unit 31004 of the People’s Liberation Army, Beijing 100091)

This paper establishes a complete dynamic model of a one-dimensional variable-centre-rolling hypersonic reentry vehicle. Research on the trajectory and attitude dynamics of the variable-cent-of-mass aircraft, and clarify the control mechanism and characteristics of the variable-cent-of-mass roll control aircraft. At the same time, the structural design of the aircraft is proposed for the analysis of the coupling effect between the channels and the influence of the actuator parameter settings on the control ability. The related requirements of control system design provide a theoretical basis for the engineering practice of hypersonic reentry aircraft.

hypersonic reentry vehicle；moving mass control；dynamics

樊博璇（1987），博士，兵器科學與技術專業；研究方向：武器系統運用與保障工程。

2020-09-23