于峰，趙志博，劉文一

高空風場對某飛行器級間分離影響研究

于峰1，趙志博2，劉文一1

（1.91550部隊41分隊，遼寧 大連 116023；2. 91049部隊，山東 青島 266102）

研究高空風場對飛行器飛行時級間分離的影響。建立飛行器分離模型，仿真計算飛行器分離時下面級的橫移距離、兩個關鍵部件之間的間隙。在最低飛行條件風速34 m/s時，下面級橫移距離和兩個關鍵間隙均滿足設計條件。在仿真計算風速40 m/s條件下，下面級橫移距離、兩個關鍵部件間隙均滿足設計要求，因此最低飛行條件風速可放寬至40 m/s。在最大風速45 m/s時，三個值均超出了設計值，有可能會對分離造成影響。高空風對飛行器級間分離存在一定影響，隨著高空風風速的增加，面級的橫移距離和兩個關鍵部件之間鍵間隙隨之增加。

高空風；分離；橫移；間隙

飛行器飛行時必須要經(jīng)歷屬于四維空間的高空風場，因此飛行器設計時既要考慮最大風，也要考慮風切變，既要考慮其空間變化，也要考慮其時間變化[1]。飛行器在級間分離時，一方面風使下面級分離體發(fā)生了相對運動，另一方面風從飛行器頭部到尾部有大小變化[2]。這些因素對飛行器的飛行和分離會產(chǎn)生較大影響，因此有必要對其進行分析，以避免高空風場對飛行器正常飛行造成影響。

1 高空風場模型

高空風對飛行器的影響主要是氣動載荷，它取決于動壓和俯仰角。高空風對飛行器有四種物理作用：風載荷、風切變、風振和陣風干擾。高空風往往是以兩種以上的效應作用于飛行器，從而對飛行器的影響更大[3]。

1.1 氣動載荷

氣動載荷可以用動壓和俯仰角的乘積來表示。當風運動遇到障礙物時，單位體積的動能施加于障礙物上的壓力為動壓。當風以一定速度被滯止時，壓力增大，在等熵條件下，將全部轉化為動能[4]。動壓的表達式為[5]：

飛行器飛行時，1為飛行速度，2為風速，=v+v，代入式（1）得：

當1≥2時，2可以忽略。

1.2 風載荷

風載荷指風吹過障礙物時，障礙物在垂直于風速方向上單位面積所受到的壓力。風壓是兩部分作用之和：一是障礙物迎風面的動壓，二是障礙物背風側的壓降。壓降的表達式為：

式中：為結構常數(shù)，對于圓柱形障礙物，=0.3。將動壓和壓降作用加起來，得到風壓的表達式[6]：

1.3 風切變

風切變?yōu)槟骋环较蛏先魏物L分量和風矢量的局地變化，通常用表示風矢量，、、分別表示、、方向的風分量[7]：

式中：、、分別為指向、、的單位矢量。飛行器飛行受到的風切變的影響，可通過建立在飛行器上的坐標系來表達。參照拉格朗日座標，飛行器飛行時的風速變化可由式（6）給出：

是,的函數(shù)，式（7）又可變?yōu)椋?/p>

式中：表示矢量的仰角[8]。

將式（8）代入式（7），整理后得到[9]：

1.4 風振

當風吹過飛行器時，在飛行器的背風側會產(chǎn)生卡曼渦旋。卡曼渦旋會使飛行器在垂直于風的方向上發(fā)生振動，這種現(xiàn)象為風振[10]。這種振動往往會使飛行器發(fā)生結構性破壞，當風振與飛行器發(fā)生共振時產(chǎn)生的破壞力更大[11]。

當飛行器在高空風場中飛行時，可將其看作是一個圓柱，其自振周期為[12]：

式中：為飛行器長度；為飛行器外徑。

2 仿真模型

2.1 物理模型

一個完整的仿真模型包括幾何模型、約束類型、初始運動值和受力等[13]。某飛行器上面級發(fā)動機噴管套在下面級前連接裙里，其模型如圖1所示。深色部分為上面級，淺色部分為下面級，同時標出了兩個關

圖1 仿真模型

鍵件的位置、高空風方向和下面級橫移方向。

對于飛行器的上面級和下面級，都有6個拉格朗日方程及其對應的約束方程[14]：

將方程（13）和方程（14）改寫成矩陣型：

式中：0表示零矩陣；為廣義力矩陣；為約束方程矩陣[16]。

在廣義坐標下，動能可以用式（16）表示：

式（15）與式（16）合并為：

方程（18）、（19）均為代數(shù)-微分方程。在進行多體動力學分析時，可將拉格朗日方程寫為一階微分方程形式，且引入[18]：

得到：

2.2 初始運動

飛行器有一個初始速度，計算時需要設置初始速度在參考坐標系上的分量值[19]，建立的坐標系如圖2所示。

圖2 坐標系

飛行器在級間分離時刻，其速度矢量為：

2.3 內(nèi)力載荷

飛行器分離時，上、下面級之間存在約束，因此上、下面級之間成對出現(xiàn)大小相等、方向相反的內(nèi)力[20]。對飛行器分離不做功的內(nèi)力為保守內(nèi)力，對飛行器分離做功的內(nèi)力稱為非保守內(nèi)力。飛行器分離，受到外來載荷或非保守內(nèi)力作用時，飛行器的上、下面級就會產(chǎn)生加速度，飛行器分離模型就會處于動態(tài)變化中[21]。

3 仿真計算及結果分析

3.1 計算條件

計算模型采用三維六自由度剛體運動學和動力學模型，考慮的受力工況主要有：下面級和上面級主發(fā)動機推力、高空風場產(chǎn)生的氣動載荷、分離時上面級和下面級的速度、質(zhì)量特性等[22]，見表1。飛行器飛行時具有一定攻角，計算時只考慮風載荷垂直加載在飛行器迎風面的分量[23]。

表1 級間分離時飛行器載荷

Tab.1 Load of aircraft during stage separation

3.2 計算結果

仿真計算了飛行器上下面級之間分離后，不同風速條件下，下面級的橫移距離、關鍵間隙1、2的距離。

3.2.1 下面級橫移距離

飛行器上下面級之間分離后，風速分別為27、30、34、38、40、45 m/s時，下面級的橫移距離如圖3所示。

3.2.2 關鍵間隙1距離

飛行器上下面級之間分離后，風速分別為27、30、34、38、40、45 m/s時，關鍵間隙1距離如圖4所示。

3.2.3關鍵間隙2距離

飛行器上下面級之間分離后，風速分別為27、30、34、38、40、45 m/s時，關鍵間隙2距離如圖5所示。

圖3 不同風速時下面級的橫移距離

圖4 不同風速時關鍵間隙1距離

圖5 不同風速時關鍵間隙2距離

3.3 結果分析

仿真計算結果表明，高空風對飛行器分離存在一定影響。隨著高空風風速的增加，飛行器下面級的橫移距離隨之增加，關鍵間隙值1和2隨之減小，見表2。設計值要求分離時下面級最大橫移距離為31 mm，關鍵間隙1的最小值為69 mm，關鍵間隙2的最小值為152 mm。

表2 風速增加與下面級橫移距離、兩個關鍵間隙值之間的關系

Tab.2 The relationship among the increased wind speed, sway distance of the lower stage and gaps of two key points

4 結論

為了研究高空風場對飛行器飛行時級間分離的影響，建立了飛行器分離模型，仿真計算了飛行器分離時下面級的橫移距離和兩個關鍵部位之間的間隙，得出以下結論。

1）高空風對上面級和下面級分離存在一定影響，隨著高空風風速的增加，下面級的橫移隨之增加，使危險點之間的間隙隨之減小，會對級間分離和飛行器正常飛行帶來影響。

2）在最低飛行條件高空風速34 m/s時，下面級橫移距離為29 mm，兩個關鍵間隙分別為73 mm和182 mm，均滿足設計條件。

3）在計算風速40 m/s條件下，下面級橫移距離、兩個關鍵部件間隙均滿足設計要求，因此最低飛行條件可放寬至40 m/s。

4）在計算最大風速45 m/s條件下，關鍵間隙1減小至63 mm，關鍵間隙2減小至140 mm。這兩個間隙均偏小，有可能會對飛行器級間分離和正常飛行造成影響。

[1] 李臣明, 張微, 韓子鵬. 100 km以下風場對遠程彈箭運動的影響研究[J]. 兵工學報, 2006, 28(10): 1170-1173. LI Chen-ming, ZHANG Wei, HAN Zi-peng. Effect of Wind Field at High Altitude Under 100 km on Long-range Shells and Rockets[J]. Acta Armamentarn, 2006, 28(10): 1170-1173.

[2] KELLY J M, PIETER G B, BANDU N P, et al. Overview of TRANSONIC to Hypersonic Stage Separation Tool Development for Multi Stage to Orbit Concepts[R]. AIAA 2004-2595, 2004.

[3] 美國國家海洋和大氣局?國家航宇局和美國空軍部.標準大氣(美國1976)[M]. 北京: 科學出版社, 1982. National Ocean and Atmosphere Bureau of America, NationalSpace Navigation Bureau and America Air Force Office. Standard Atmosphere(US,1976)[M]. Beijing: Science Pros,1982.

[4] 馬瑞平, 廖懷哲. 中國地區(qū)20~80 km高空風的一些特征[J]. 空間科學學報, 1999, 19(4): 334-341. MA Rui-ping, LIAO Huai-zhe. The Characteristics of Winds at Height of 20~80 km in the Chinese Area[J]. Chinese Journal of Space Science, 1999, 19(4): 334-341.

[5] GUNGLE R L, BROSIER W S, LEONARD H W. An Experimental Technique for the Investigation of Tipoff Associated with Stage Separation of Multistage Rocket Vehicles[R]. NASATND-1030, 1962.

[6] 劉文一, 楊春, 楊培源. 固體發(fā)動機金屬材料和復合材料殼體水下力學特性分析[J]. 航天制造技術, 2017, 205(5): 47-49. LIU Wen-yi, YANG Chun, YANG Pei-yuan. Mechanical Property Evaluation of Solid Rocket Motor Case of Steel Material and Composite Material under Water[J]. Aerospace Manufacturing Technology, 2017, 205(5): 47-49.

[7] DE CASTRO L N, VON ZUBEN F J. The Clonal Selection Algorithm with Engineering Applications[C]// Proceedings of GECCO`00, Workshop on Artificial Immune Systems and Their Applications. Las Vegas, USA, 2000.

[8] KECHRIOTIS G I, MANOLAKOS E S. Hopfield Neural Network Implementation of Tile Optimal CDMA Multiuser Detector[J]. IEEE Transactions on Neural Networks(S1045-9227), l996, 7(1): 131-141.

[9] 劉文一, 焦冀光. 某飛行器復合材料薄壁加筋結構艙段穩(wěn)定性分析[J]. 艦船電子境工程, 2016, 36(2): 110-113. LIU Wen-yi, JIAO Ji-guang. Stability Analysis of a Aircraft`s Thin-walled Reinforced Composite Material Cabin[J]. Ship Electronic Engineering, 2016, 36(2): 110-113.

[10] 牛躍聽, 穆希輝, 姜志保. 某型火箭彈控制艙環(huán)境適應性研究[J]. 裝備環(huán)境工程, 2014, 11(1): 86-92. NIU Yue-ting, MU Xi-hui, JIANG Zhi-bao. Research on Environmental Adaptability of a Guided Rocket[J]. Equipment Environmental Engineering, 2014, 11(1): 86-92.

[11] 劉文一, 焦冀光. 固體發(fā)動機裝藥熱安全性數(shù)值分析[J]. 裝備環(huán)境工程, 2016, 13(2): 130-131. LIU Wen-yi, JIAO Ji-guang. Numerical Analysis on the Thermal Safety of Solid Rocket Motor Propellant[J]. Equipment Environmental Engineering, 2016, 13(2): 130-131.

[12] 趙保平, 孫建亮, 蔡駿文, 等. 航天動力學環(huán)境的最新進展與技術展望[J]. 裝備環(huán)境工程, 2015, 12(3): 8-14. ZHAO Bao-ping, SUN Jian-liang, CAI Jun-wen, et al. Recent Development and Prospect of Astrodynamic Environmental Technology[J]. Equipment Environmental Engineering, 2015, 12(3): 8-14.

[13] 李堯. 飛機溫度環(huán)境適應性要求分析和確定技術探討[J]. 裝備環(huán)境工程, 2008, 5(6): 60-64. LI Yao. Analysis and Confirmation Techniques for Requirements of Temperature Environmental Worthiness of Aircraft [J]. Equipment Environmental Engineering, 2008, 5(6): 60-64.

[14] 劉文一, 李玉龍, 吳訓濤. 基于分布水動壓力的飛行器水下載荷計算方法研究[J]. 導彈與航天運載技術, 2019, 336(1): 20-23. LIU Wen-yi, LI Yu-long, WU Xun-tao. Calculation method of Underwater Load of Aircraft Based on Distributed Hydrodynamic Pressure[J]. Missiles and Space Vehicles, 2019, 336(1): 20-23.

[15] 劉章龍, 趙徐成, 胡濤. 基于高原環(huán)境的保障裝備適應性技術研究[J]. 裝備環(huán)境工程, 2016, 13(2): 35-36. LIU Zhang-long, ZHAO Xu-cheng, HU Tao. Adaptive Technology of Support Equipment Based on Plateau Environment[J]. Equipment Environment Engineering, 2016, 2016, 13(2): 35-36.

[16] SENOCAK I, SHYY W. A Pressure-based Method for Turbulent Cavitating Flow Computations[J]. Journal of Computational Physics, 2002, 176(2): 363-383.

[17] EDWARDS J R, LIOU M S. Simulation of Two-phase Flows Using Low-diffusion Shock-capturing Schemes[C]// 44th AIAA Aerospace Sciences Meeting and Exhibit. Reno, Nevada, 2006.

[18] LEE C, BYUN D. Cavitation Flow Analysis of Axisymmetric Bodies Moving in the Water[J]. ICCSA LNCS, 2006, 3984: 537-545.

[19] CHEN Y, LU C J, GUO J H. Modeling and Simulation of Cavitating Flow over Underwater Vehicle with Large Angle of Attack [J]. Journal of Ship Mechanics, 2010, 14(12): 1319-1330.

[20] 李強兵. 基于L-M法的某型UUV系統(tǒng)可靠性評估及應用[J]. 機械, 2015, 42(5): 32-33. LI Qiang-bing. Comprehensive Evaluation of System Reliability and Application of a Certain UUV Based on L-M[J]. Machinery, 2015, 42(5): 32-33.

[21] 李玉杰, 石顏青, 羅振兵. 水滴結冰結霜及合成雙射流的防霜防冰實驗[J]. 氣體物理, 2017, 2(5): 37-45. LI Yu-jie, SHI Yan-qing, LUO Zhen-bing. An Experimental Investigation on the Process of Droplet Icing/Forsting and Anti-frosting /Anti-icing Using Dual Synthetic Jet[J]. Physics of Gases, 2017, 2(5): 37-45.

[22] HAYASHI Y, AOKI A, ADACHI S, et al. Study of Frost Properties Correlating with Frost Formation Types[J]. Heat Transfer, 1997(99): 239-245.

[23] FEUILLEBOIS A, LASEK A, CREISMEAS P, et al. Freezing of a Subcooled Liquid Droplet[J]. Colloid Sci, 1995(169): 90-102.

Effect of High Altitude Wind Field on Stage Separation of an Aircraft

YU Feng1, ZHAO Zhi-bo2, LIU Wen-yi1

(1. Unit 41 of PLA 91550, Dalian 116023, China; 2. PLA 91049, Qingdao 266102, China)

The work aims to research the effect of wind field at high altitude on the stage separation of aircraft during flying. A stage separation model of aircraft was established. The sway distance of the lower stage and the gaps of two key points of aircraft after its stage separation were simulated. The sway distance of the lower stage and the gaps of two key points of aircraft all met the design conditions and requirements when the minimum wind speed and the simulated one were respectively 34 m/s and 40 m/s. Therefore, the minimum wind speed for flying could be relaxed to 40 m/s. When the maximum wind speed was 45 m/s, the sway distance and the gaps above all exceeded the design values, so the stage separation would be affected by it. It proves that the stage separation of aircraft will be affected by high altitude wind in some way, and the sway distance of the lower stage and the gaps of two key points of aircraft will be increased with the increase of high altitude wind speed.KEY WORDS: high altitude wind; separation; sway; gap

2019-09-06；

2019-10-02

10.7643/ issn.1672-9242.2020.03.018

V112.71；O358

A

1672-9242(2020)03-0108-06

2019-09-06；

2019-10-02

于峰（1967—），男，碩士，高級工程師，主要研究方向為武器裝備鑒定。

YU Feng (1987—), Male, Master, Senior engineer, Research focus: weapons and equipment identification.