涡环旋转伞系统减速导旋效率研究

马晓冬, 郭 锐, 刘荣忠, 胡志鹏, 吕胜涛

(1. 南京理工大学机械工程学院, 智能弹药技术国防重点学科实验室，江苏 南京 210094; 2. 哈尔滨建成集团有限公司, 黑龙江 哈尔滨 150030; 3. 中国兵器工业第203研究所, 陕西 西安 710065)

涡环旋转伞系统减速导旋效率研究

马晓冬1,2, 郭 锐1,*, 刘荣忠1, 胡志鹏3, 吕胜涛1

(1. 南京理工大学机械工程学院, 智能弹药技术国防重点学科实验室，江苏 南京 210094; 2. 哈尔滨建成集团有限公司, 黑龙江 哈尔滨 150030; 3. 中国兵器工业第203研究所, 陕西 西安 710065)

为了解涡环旋转伞的减速导旋特性，提高旋转伞系统的减速导旋效率，以满足某些载物工作状态对转速、落速的苛刻需求，设计一种由涡环旋转伞和圆形减速伞构成的组合伞系统。制作单伞系统和组合伞系统的试验模型，分别进行伞塔试验。基于伞塔试验的高速摄影和图像处理技术，提取伞物系统的弹道数据点，结合弹道模型，用最小二乘法拟合时间弹道数据，得到单伞系统和组合伞系统的极限速度和阻力特征。基于姿态存储测量法，获得伞物系统的转速变化规律。对单伞系统和组合伞系统的弹道规律及减速导旋效率等数据进行对比分析，提出可有效提高旋转伞系统转速落速比的方法。结果表明：对高速录像进行图像处理可简单快速有效地提取伞物系统的极限速度和阻力特征；单具涡环旋转伞可提供更大扭矩，使载物旋转加速度更大，其导旋效率高；组合伞可提供更大阻力，使载物极限落速更低，其减速效率高；组合伞系统的转速落速比更大，整体性能优于单伞系统。

涡环旋转伞；组合伞；伞塔试验；减速导旋效率；转速落速比

0 引 言

涡环旋转伞是一种常见的旋转降落伞。伞衣的高速旋转，使得带涡环旋转伞的物伞系统在下降时具有良好的稳定性。此外，它还具有阻力系数大、开伞动载小、易维护等优点，也被广泛应用于兵器的弹道控制、飞机的着陆减速和控制、空降空投、航天飞行器的安全回收技术以及外太空的探险和着陆等领域[1-2]。

在工程界和学术界，国内外众多研究者对平面圆形或锥形等结构轴对称伞系统进行了大量的理论研究和探索，但对旋转降落伞系统缺乏比较系统的研究。旋转伞系统的稳定落速和稳定转速是两个重要指标参数，直接决定其工作性能，稳定状态下的转速落速比是评价旋转伞系统的一个重要指标。关于旋转降落伞的已有研究多涉及旋转伞-载物系统的动力学建模[3-10]，对弹道特性和动力学特性进行仿真分析，很少有关于旋转伞系统减速导旋效率的讨论。特别地，文献[11]通过理论分析和伞塔投放试验方法研究了涡环旋转伞的开伞充气稳定性及其减速导旋运动特性，获得了伞衣幅结构及伞绳连接方式对涡环旋转伞减速导旋运动的影响规律，但未指出提高转速落速比的有效措施。

为进一步了解涡环旋转伞的减速导旋特性，提高旋转伞系统的转速落速比，设计了一种由涡环旋转伞和圆形减速伞构成的组合伞系统。针对伞塔投放试验数据，利用图像处理技术和姿态存储测量技术，得到单伞系统和组合伞系统的落速和转速变化，进一步对比分析二者的减速导旋效率，得到提高涡环旋转伞系统转速落速比的有效措施。

1 理论模型

本文研究的涡环旋转伞系统主要由5部分组成，包括伞衣幅、伞绳、伞盘、旋转接头和载物，如图1所示。涡环旋转伞的结构复杂，一般由4片伞衣幅和多根不同长度的伞绳连接而成。伞衣幅上不对称的开口设计可以实现在充满的情况下，形成一定的凸面和倾斜度，在不对称空气动力和力矩的作用下实现旋转性能，伞盘通过旋转接头带动载物旋转。

1.1 涡环旋转伞系统减速运动模型

试验求解涡环旋转伞系统极限落速和阻力特征过程中，假设涡环旋转伞充满后形状不再变化，忽略姿态变化，将其简化为仅受空气阻力和重力的单刚体运动，则其三自由度运动方程为:

式中，v=[vxvyvz]T为系统相对于地面的速度，x为系统空间位置向量，g为重力加速度，mp和md分别为涡环旋转伞和载物的质量，ρ为空气密度，(CpSp+CdSd)为系统的阻力特征，Cp和Cd分别为涡环旋转伞和载物的阻力系数，Sp和Sd分别为二者的迎风面积。 只考虑系统竖直方向的运动，即vy=vz=0。当系统达到平衡状态时，dv/dt=0，由式(1)得系统极限速度

由文献[9]可知由极限速度表示的系统运动方程为

式中，x为系统下落距离，t为飞行时间。

由式(3)知，若已知系统飞行时间向量T[t1，t2，…，tn]及对应高度向量X[X1，X2，…，Xn]，则可通过拟合的方法计算得到伞塔试验中涡环旋转伞系统的极限速度，再通过式(2)计算得到阻力特征。

1.2 涡环旋转伞系统转动运动方程

假设伞轴与铅垂轴平行，且系统仅绕伞轴旋转。涡环旋转伞的转动方程为:

式中，ωb为涡环旋转伞的稳定转速。

2 伞塔试验与数据处理

2.1 试验环境与条件

伞物系统试验模型从高塔顶部100 m处投放，初始落速和转速均为零。试验用Red lake dk100高速录像机，拍摄频率1000 Hz，图片像素640×480，距投放点地面投影100 m。标志物为铅垂放置且有固定尺寸刻度的绳索，如图3和图4所示，通过测量底部标志物N点的实际高度确定各标志物的高度。为准确记录伞物系统下落时间和转速变化，在载物上装有姿态参数存储测量仪。

2.2 基于图像提取落速

伞物系统自由飞行时间较短，天空背景几乎不变，利用Surendra算法[13]提取高速摄影图片中的天空背景。选取伞物系统在视场中的高速摄影图片序列，与天空背景差分并去噪后得到伞物系统运动图像序列[14-16]，最终得到伞物系统的上下左右像素极值。本文关注铅垂方向的运动，故可由载物下缘像素坐标来确定伞物系统的位置和运动。

为得到伞弹系统下落过程中的实际高度，可在高速录像拍摄的背景中设置两个与地面高度确定的标尺M和N，如图3所示，其中YM和YN为图片中标尺M和N的铅垂像素坐标，HM和HN分别为二者距地面的实际距离，YP为载物下缘的铅垂像素坐标。伞物系统的实际高度为：

高速录像图片序列中伞物系统的实际高度组成高度向量H。设伞物系统初始高度为H0，则其下落距离向量X=H0I-H。

设高速录像的频率为f，则拍摄第n张图片至第n+i张图片的飞行时间为((n+i)-n)/f。设第一张图像拍摄于伞物系统抛出后的t1秒，则第n张图片中伞物系统的飞行时间为t1+n/f。根据选取的高速录像图片序列，可确定对应的伞物系统飞行时间向量T。

已得到原始数据X和T，且二者之间满足式(3)，其中V为待定系数，则通过最小二乘法[17]可拟合出伞物系统的极限速度V，进而求出其阻力特征。

初始时间t1未知，可通过逐渐增大试取t1值，由拟合误差判断是否达到要求。当误差小于设定值时，认为所取t1值与试验中伞物系统下落到图像拍摄时刻吻合。

2.3 基于姿态参数存储测量仪提取转速

姿态参数存储测量仪采用三轴加速度计和三轴磁阻传感器，获得3路加速度分量和3路地磁分量，如图5所示，采样频率500 Hz。试验结束回收测量仪，基于MATLAB调用存储数据进行数据处理，可得到试验模型在下降过程中的转速变化和飞行时间。具体实现过程如下：

(1) 利用三轴转台给定载物固定转速，对姿态参数存储测量仪进行标定；

(2) 读入地磁分量数据;

(3) 鼠标选取起始(Start falling)至终点(Touching the ground)数据段；

(4) 利用采样频率将时间坐标转化为秒，步骤(3)中截取数据段的时长即系统的飞行时间，利用测量仪标定曲线将地磁分量数据转化为转速；

(5) 对得到的转速曲线光滑处理。

3 结果分析与验证

3.1 结果分析

处理伞塔试验数据，得到系统极限速度、阻力特征及稳定转速，如表1所示。由于组合伞的涡环旋转伞和减速伞均具有减速作用，因此组合伞-载物系统的极限速度较小，阻力特征较大，即比单伞系统具有更好的减速效率。伞衣迎风面积远大于载物迎风面积，若忽略载物对系统阻力的影响，则可得到涡环旋转伞和组合伞的阻力系数，分别为1.06和1.47。组合伞系统的稳定转速稍小于单伞系统，但转速落速比明显大于单伞系统。

表1 伞物系统主要运动数据Table 1 Main motional data of parachute-payload system

图6为伞物系统转速变化曲线。在导旋初期，单伞系统的旋转加速度较大，在较短的时间内达到最大转速，随后转速小范围波动，认为达到稳定状态。组合伞系统的旋转加速度较小，在下落过程中，转速平稳地逐渐增大，在即将落地时转速才达到最大值。旋转加速度的大小主要由涡环旋转伞的大小决定。伞衣幅较大的涡环旋转伞展开充满后，可为载物提供较大的扭矩，载物一定的情况下，其旋转加速度更大，所以单伞系统比组合伞系统具有更好的导旋效率。

3.2 结果验证

为验证结果的准确性，同时直观地进行比较，将图像处理、数据拟合得到的阻力特征代入伞物系统三自由度运动方程(1)，得到伞物系统自由下落过程中的高度变化曲线和速度变化曲线，如图7和图8所示。图7得到：伞物系统试验飞行高度变化与理论计算结果吻合良好，组合伞系统的滞空时间明显长于单伞系统。图8得到：少数实时平均速度大小出现跳动变化，这是由于图片像素较低，载物轮廓比较模糊，每次进行图像提取载物下缘位置不完全一致而产生的误差。但此时系统落速不再明显变化，可忽略这些少数数据点。由平均速度法测得的速度值与理论速度曲线吻合良好。

从图7可看出，单伞系统和组合伞系统从100 m高下落到地面分别用时约9 s和10.5 s，而从如图5所示的姿态原始数据中可知二者实际飞行时间约为10 s和9 s。理论计算和伞塔试验结果差别较大，但二者并不矛盾。理论计算假设伞物系统一直保持形状不变，即初始状态下伞即处于充满状态；而伞塔试验中下落初期，伞要进行充气，且受到横风的影响，受随机因素影响相对较大。处理得到的伞塔试验数据主要是伞物系统稳定状态下的数据，与理想状态下的理论计算结果吻合，所以本文的数据处理方法是可行的。

4 结 论

1) 通过Surendra方法提取伞物系统背景图像，对高速录像拍摄图像进行运动目标提取，获得伞物系统下落的时间、高度向量，结合运动模型，利用最小二乘拟合方法求解自由下落状态的极限速度和阻力特征。该方法快速简单有效。

2) 单伞系统的伞衣幅面积大，展开后可为载物提供更大的扭矩，使载物具有更大的旋转加速度，在较短的时间内转速达到最大，即单伞系统具有更好的导旋效率。

3) 组合伞系统的涡环旋转伞和圆形减速伞均具有良好的减速效果，故其阻力特征明显大于单伞系统，极限落速小于单伞系统，即组合伞系统具有更好的减速效率；本文设计的涡环旋转伞和组合伞的阻力系数分别约为1.06和1.47。

4) 稳定状态下，组合伞系统的转速落速比大于单伞系统，组合伞系统的工作性能优于单伞系统。

[1]Wang L R. The theory and application of parachute[M]. Peking: Aerospace press, 1997. (in Chinese)王利荣. 降落伞理论与应用[M]. 北京: 宇航出版社, 1997[2]Ewing E G, Bixbu H W, Knacke T W. Recovery systems design guide[R]. AD A070251, Dec, 1978.

[3]Pillasch D W, Shen Y C. Parachute/submunition system coupled dynamics[R]. AIAA-84-0875, 1984.

[4]Doherr K F. Nine-degree-of-freedom simulation of rotating parachute systems[J]. Journal of Aircraft, 1992, 29(5): 774-781.

[5]Guo R, Liu R Z, Huang F F, et al. Study on the decelerating and rotating dynamics model for a certain rotating parachute system[C]//Third International Conference on Information and Computing, 2010: 216-218.

[6]Liu S P, Han Z P. A drag model for parachute-projectile systems of sensor-fuzed munitions[J]. Acta Armamentarii, 1997, 18(3): 221-225. (in Chinese)刘世平, 韩子鹏. 末敏弹旋转伞-弹系统阻力模型研究[J]. 兵工学报, 1997, 18(3): 221-225.

[7]Tang Q G, Zhang Q B, Zhang X J, et al. Nine-degree-of-freedom model of bomb-parachute system[J]. Acta Armamentarii, 2007, 28(4): 449-452. (in Chinese)唐乾刚, 张青斌, 张晓今, 等. 伞-弹系统九自由度动力学模型[J]. 兵工学报, 2007, 28(4): 449-452.

[8]Tang Q G, Wang Y, Zhang Q B, et al. Application of dynamics and kinematics of parachute-bomb system in target identification for a target sensitive projectile[J]. Acta Armamentarii, 2007, 28(7): 796-799. (in Chinese)唐乾刚, 王昱, 张青斌, 等. 伞-弹动力学及运动学在末敏弹目标识别中的应用[J]. 兵工学报, 2007, 28(7): 796-799.

[9]Guo R, Liu R Z. Dynamics model and simulation of rigid and flexible coupling system for terminal-sensitive submunition[J]. Acta Armamentarii, 2007, 28(1): 10-14. (in Chinese)郭锐, 刘荣忠. 末敏弹刚柔耦合系统动力学模型及仿真[J]. 兵工学报, 2007, 28(1): 10-14.

[10]Yin K G. Study on the motion characteristic of terminal-sensitive submunition[D]. Nanjing University of Science and Technology, 2008. (in Chinese)殷克功. 末敏子弹运动特性分析研究[D]. 南京理工大学, 2008.

[11]Guo R, Liu R Z, Hu Z P, et al. Study on the inflated stability and motional characteristics of vortex ring parachute canopy[J]. Acta Aerodynamica Sinica, 2013, 31(6): 733-738. (in Chinese)郭锐, 刘荣忠, 胡志鹏, 等. 涡环旋转伞开伞稳定性及减速导旋运动特性研究[J]. 空气动力学学报, 2013, 31(6): 733-738.

[12]Zhang Y J. Advanced dynamics[M]. Nanjing: Nanjing University of Science and Technology Press, 2003. 张延教. 高等动力学[M]. 南京: 南京理工大学出版社, 2003.

[13]Wang Z Q, Liu F Q. Comparison of adaptive extract algorithm of object scene[J]. Computer Engineer, 2008, 34(23): 220-223. (in Chinese)王正勤, 刘富强. 自适应背景崎岖算法的比较[J]. 计算机工程, 2008, 34(23): 220-223.

[14]Fan Y, Guo C G. A method for detection and recognition of moving objects in sequential images[J]. Journal of Dalian University of Technology, 2004, 44(1): 122-126. 范盈, 郭成安. 一种运动图像的检测与识别技术[J]. 大连理工大学学报, 2004, 44(1): 122-126.

[15]Liu G X, Shao M L, Liu X H, et al. Video moving object auto-extraction in real scene[J]. Acta Optica Sinica, 2006, 26(8): 1150-1155. (in Chinese)

刘贵喜, 绍明礼, 刘先红, 等. 真实场景下视频运动目标自动提取方法[J]. 光学学报, 2006, 26(8): 1150-1155.

[16]Hu Z P, Liu R Z, Guo R, et al. TSP drag coefficient calculation based on projectile image processing[J]. Journal of Detection & Control, 2012, 34(3): 8-12. (in Chinese)胡志鹏, 刘荣忠, 郭锐, 等. 基于弹丸图像的末敏弹阻力系数计算方法[J]. 探测与控制学报, 2012, 34(3): 8-12.

[17]Wan J W, Xiong D L, Huang F K. Robust estimation for the coefficient of resistanceCD[J]. Acta Armamentarii, 1991(1): 21-28. (in Chinese)万建伟, 熊栋梁, 皇甫堪. 弹丸阻力系数CD的Robust估计[J]. 兵工学报, 1991(1): 21-28.

Study on decelerating and spinning efficiency of vortex ring parachute system

Ma Xiaodong1,2, Guo Rui1,*, Liu Rongzhong1, Hu Zhipeng3, Lyu Shengtao1

(1.SchoolofMechanicalEngineering,ZNDYMinisterialKeyLaboratory,NanjingUniversityofScience&Technology,Nanjing210094,China；2.HarbinJianchengGroupLimitedCompany,Harbin150030,China;3.No.203ResearchInstituteofChinaOrdnanceIndustries,Xi’an710065,China)

To understand the decelerating and spinning characteristics of vortex ring parachute, increase the decelerating and spinning efficiency of rotating parachute system, and satisfy the harsh demand of limiting falling speed and rotating speed for some payloads at specific working state, a compound parachute system composed of a circle drag parachute and a vortex ring parachute was designed. The test models of single parachute system and compound parachute system were made out, and then two tower tests were carried out separately. Based on high speed photography of parachute tower test and image processing technology, the trajectory data points of the parachute-payload system were extracted. Combining with the trajectory model, the time and trajectory data was fitted by least square method, the limiting falling speed and the resistance characteristics of the system were obtained. Using attitude memory-measuring technology, the rotating speed changing law of parachute-payload system was gained. The trajectory laws and the decelerating and spinning efficiency of the vortex ring parachute and the compound parachute system were analyzed contrastively, and then the method to increase the ratio of rotating speed and limiting falling speed availably was proposed. The results show that through processing the pictures shot by high-speed camera, the limiting falling speed and resistance characteristics of the parachute-payload system can be extracted quickly and efficiently. The single vortex ring parachute could provide bigger torque so that the rotating accelerating of the payload is higher, which means that the spinning efficiency is better. The compound parachute could bigger resistance so that the limiting falling speed of the payload is lower, which means that the decelerating efficiency is better. The radio of rotating speed and falling speed of compound parachute system is bigger than that of single parachute system, which means that the whole performance of compound parachute system is better.

vortex ring parachute; compound parachute; parachute tower test; decelerating and spinning efficiency; ratio of rotating speed and falling speed

0258-1825(2017)01-0057-06

2015-01-05;

2015-02-27

国家自然科学基金(11372136); 国家自然科学基金-青年科学项目(11102088); 江苏省研究生培养创新计划项目(CXLX12-0210)

马晓冬(1988-),男,黑龙江宝清人,工程师,博士,研究方向:弹药的灵巧化与智能化. E-mail: bqnj6222007@126.com

郭锐*(1980-),男,山东海阳人,副教授,博士,研究方向:弹药的灵巧化与智能化. E-mail: guoruid@163.com

马晓冬, 郭锐, 刘荣忠, 等. 涡环旋转伞系统减速导旋效率研究[J]. 空气动力学学报, 2017, 35(1): 57-61.

10.7638/kqdlxxb-2015.0003 Ma X D, Guo R, Liu R Z, et al. Study on decelerating and spinning efficiency of vortex ring parachute system[J]. Acta Aerodynamica Sinica, 2017, 35(1): 57-61.

V211.73

A doi: 10.7638/kqdlxxb-2015.0003