基于龙格库塔算法对降落伞空投的研究概要.docx
- 文档编号:10803895
- 上传时间:2023-02-23
- 格式:DOCX
- 页数:15
- 大小:210.32KB
基于龙格库塔算法对降落伞空投的研究概要.docx
《基于龙格库塔算法对降落伞空投的研究概要.docx》由会员分享,可在线阅读,更多相关《基于龙格库塔算法对降落伞空投的研究概要.docx(15页珍藏版)》请在冰豆网上搜索。
基于龙格库塔算法对降落伞空投的研究概要
基于龙格2库塔算法对降落伞空投的研究
王景志,陈红林
(西北工业大学电子信息学院,陕西西安 710129
摘 要:
以降落伞空投系统中的降落伞拉直和充气过程以及物伞系统的稳定降落过程为研究目标,采用动力学分析法建立了物伞空投系统完整的运动模型,提出了多吊带系统分析和物伞系统夹角变化的新方法。
根据已有的降落伞性能计算模型和物伞系统运动模型,采用龙格2库塔法对降落伞空投全过程进行了仿真。
通过参数调试来分析空投高度﹑降落伞的名义面积等因素对降落伞空投系统的影响,指出了实际系统的改进方向,为降落伞空投系统的设计提供了理论依据。
关键词:
降落伞空投;运动模型;龙格2库塔法;计算机仿真
中图分类号:
TN911.7;TP391.9 文献标识码:
A 文章编号:
10042373X(20101420124203
ResearchforParachuteAirdropBasedonRunge2KuttaAlgorithm
WANGJing2zhi,CHENHong2lin
(DepartmentofElectronicsandInformation,NorthwesternPolytechnicUniversity,Xi’an710129,China
Abstract:
Focusedontheparachutedeploymentandinflation,andthetrajectoryofcargoparachutesysteminthepara2chuteairdropsystem,thetrajectorymodelofcargoparachutesystemwasbuiltwiththedynamicsmethod.Theanalyticmethodofthemulti2slingsystemandanglebetweenthecargoandparachutesystemisproposed.TheprocessoftheparachuteairdropwassimulatedwithRunge2Kuttaalgorithmonthebasisofparachuteperformancemodelandcargeparachutesystemtrajectorymodel.Withtheanalysisoftheeffectofeveralfactorsconcernedwiththeairdropsuchasheightofairdropandpar2achuteresistancecoefficientontheparachuteairdropsystem.Theimprovementdirectionoftheavailableactualsystemsispointedoutwhichprovidesatheoreticbaseforthedesignofcargoparachuteairdropsystem.
Keywords:
parachuteairdrop;trajectorymodel;Runge2Kuttaalgorithm;computerizedsimulation
收稿日期:
2010202225
2008年我国大型飞机项目正式立项,大型飞机项
目是一个国家工业、科技水平和综合实力的集中体现,对增强中国的综合国力、科技实力和国际竞争力,使中国早日实现现代化具有极为重要的意义。
由于载机的载重、货舱尺寸和接口关系不同,与其相适应的各型空投系统就有很大不同。
因此,在载机研发的同时,需要发展空投系统和精确空投技术。
空投是指利用降落伞将人员和物资投送到指定地点的一种技术,是部队投送弹药装备等物资补给和兵力部署的重要手段。
空投系统的研制必须全面考察系统运动特性以保障空投过程的安全。
一般研究手段包括投放试验、风洞试验和数字仿真等。
国内外对于各种降落伞空投系统[1]开展了很多研究,建立了不同复杂程度的计算模型,但是没有对降落伞的开伞以及附加质量等问题进行分析研究,这些是精确空投系统仿真必须考虑的关键问题。
为此,本文将降落伞空投过程进行了全面的系统过程建模和仿真分析。
降落伞性能计算方法包括基于试验的分析和基于计算结构力学、计算流体力学及流固耦合的数值模拟。
但是由于算法和计算能力的限制,数值模拟方法在现阶
段仅限于研究部分特殊问题。
因此本文的降落伞性能计算采用由试验获取的经验公式。
1 降落伞空投工作过程
本文分析的空投系统采用自动开启降落伞[223]的空投方法。
在装有货物的平台离机瞬间,启动自动开伞装置,降落伞开始充气。
当降落伞充气到一定程度时,收口绳受力,使降落伞第一次充气呈“灯泡”状,以减小降落伞的开伞动载,延时数秒后,切割器将收口绳切断。
降落伞完全张开并稳定下降直至着陆。
2 降落伞充满时间计算[4]
假设:
在充气过程中,空气密度不变;在充气过程中,伞衣阻力系数不变;在充气过程中,伞衣有效透气量不变。
(1降落伞收口充气阶段。
由实验数据得,降落伞的充满时间与拉直速度成反比。
同时,降落伞越大,充满时间越长。
由此可以得出降落伞充满时间的—般经验公式:
科学计算与信息处理王景志等:
基于龙格2库塔算法对降落伞空投的研究
tm1=(kD0/vn
L
(1
式中:
vL为拉直速度,即空投时飞机的速度;D0为伞衣名义直径;k,n为由测量数据推算所得的系数,对于平面圆形伞k=8,n=0.9。
(2降落伞保持收口阶段。
这个阶段的时间tm2,通常在时序中预先设置好,本文中tm2=2s。
(3降落伞第二次充气阶段。
这个阶段的充气时间为:
tm3=(kD0/v
(2
式中:
v为解除收口时系统的速度。
故降落伞充满时间tm为:
tm=tm1+tm2+tm3
(3
3 降落伞附加质量[4]
把降落伞的附加质量分成两部分进行计算:
伞衣内含质量Mn;伞衣表观质量Mb。
(1伞衣内含质量Mn即伞衣内所包含的空气质量,计算公式为:
Mn=ρHV=
2ρHD30
3
π21.058-(T-1.3121.62(4
ρH=ρ0(1-H/443304.265
(5式中:
V为伞衣内的容积,ρH为在H高度时空气的密
度,ρ0=1.225kg/m3
T=t/tm,D0为伞衣名义直径。
(2伞衣表观质量Mb,它只在做不稳定运动时才
体现出来,计算公式为:
Mb=0.25πD3
π3T3
2
ρHT=ρHD30
4
π2
T52(6
故物伞系统的总质量M为:
M=Mn+Mb+Ms+Mw
(7
式中:
Ms为降落伞本身的质量;Mw为货物的质量。
4 物伞系统计算模型[526]
在地面坐标系OXY中建立物伞系统的运动数学模型。
原点O选取在货物平台离机瞬间的货物平台质心在地面上的投影位置;OY轴沿重力方向指向下;OX轴垂直于OY轴指向物伞系统运动方向。
假设条件:
降落伞和货物之间的连接吊带考虑成刚性轻质杆;降落伞是轴对称的,伞衣完全充满后具有固定的形状;由于降落伞运动的非定常性而引起的气动力或力矩的增量部分用附加质量来表示;只考虑伞衣产生的气动力,并且认为伞衣的气动压心与其几何中心相重合。
本文分析的空投系统采用四吊带的连接方式。
如图1所示,图中货物质心为O;连接绳上端与降落伞系统连接点为A;下端点B与4根分吊带Oi(
i=1,2,3,4相连;连接绳AB的延长线与货物相交于C点,物伞系统夹角(降落伞伞轴线
与货物轴线之间夹角为δ。
图1 物伞系统示意图
利用C点的位置来判断各吊带的受力形式以及物伞系统夹角δ的变化:
当C点落在OOi(i=1,2,3,4上时,吊带BOi受力,其他吊带松弛;
当C点落在三角形OOiOj以及边OiOj(i,j=1,2,3,4上时,吊带BOi,BOi受力,其他吊带松弛;
当C点和O点重合时,所有吊带均受力,δ=0。
物伞系统在空投过程中受2个力的作用:
重力G,
降落伞的阻力F。
G=Mg
(8
在伞衣充气和拉直阶段降落伞的阻力F为:
F=CD12ρHv
2π4D0t
tm2
(9
式中:
CD为阻力系数;v为物伞系统速度。
当t=tm时伞
衣完全充满,物伞系统进入稳定降落阶段。
在稳定降落阶段降落伞的阻力F为:
F=CD12
ρHv2
A0(10
式中:
A0为降落伞名义面积:
A0=
π4
D2
0。
物伞系统的运动方程为:
f(X,t=dX
dt
=v0cosθ0-∫
t
0Fcosθ
M
dtf(Y,t=dYdt=v0
sinθ0
+
∫t
Mg-Fsinθ
Mdtf(δ
t=d
δdt=δ0
-∫t
MgLAO
cosθIZ
dt
(11
θ=arctan[f(Y,t/f(X,t](12
式中:
相对质心转矩IZ;LAO为A点到O点的直线距离。
式(11为物伞系统的基本运动方程。
可以看出这是一些非线性的运动微分方程。
在通常情况下,这组微分方程只能通过数值计算求解,所得的结果为一系列离散的数据。
5 系统仿真结果及分析
根据前面建立的降落伞性能计算模型和物伞系统运动模型,应用Matlab仿真技术[7],采用龙格2库塔算
5
21《现代电子技术》2010年第14期总第325期
计算机应用技术
法对物伞系统在伞衣充气和拉直阶段与稳定降落阶段全过程进行仿真[8]。
5.1 龙格2库塔算法介绍[9]
假定微分方程:
dy
dx
=f(x,y(13
及初始条件:
y(x0=y0
(14
其迭代格公式为:
yn+1=yn+(1/6(K1+2K2+2K3+K4
K1=hf(xn,yn
K2=hf(xn+h/2,yn+K1/2K3=hf(xn+h/2,yn+K2/2K4=hf(xn+h,yn+K3
(15
龙格2库塔法的主要优点是计算精度高,式(15的
截断误差为O(h5,仿真时步长选取:
h=0.01。
5.2 仿真流程与结果
设参数及初始条件:
Mw=2000kg;v0=200m/s;θ0=90°
;LA0=10m;CD=0.873,初始充气时间t1=0;稳定降落段初始时间t2=0。
系统仿真流程如图2所示,仿真结果如图3~图5所示
。
图2
系统仿真流程图
图3
空投高度和降落伞名义面积对空投距离的影响
图4
空投高度和降落伞名义面积对空投时间的影响
图5 物伞系统夹角随空投时间变化
6 结 语
由仿真结果得:
空投距离随的空投高度增加而增加,相对应高度的增加量来看,高度的增加对投放距离的影响明显,空投高度对空投时间影响不大,因此空投高度的增加可以实现远距离空投。
在选择降落伞时,应在条件许可的情况下,应选择名义面积较小的降落伞,这样不仅能够增加空投距离,还能减少货物在空中降落所需的时间,较大程度的保证运输机和货物的安全。
物伞系统夹角随空投时间近似呈阻尼变化,最后系统稳定降落。
龙格2库塔算法计算精度较高,能满足通常的计算要求;每次计算yn+1,只用到前一步的计算结果yn。
因此,在已知y0的条件下,就可以自动地进行计算,并且计算过程中改变步长不受限制。
参 考 文 献
[1]STRONGTed.DispellingtheMythsandFearsofManned
VehicleAirdrops[C].18thAIAAAerodynamicDeceleratorSystemsTechnologyConferenceandSeminar,2005.[2]BAGDONOVICHB,DESABRAISKJ,BENNEYR.Over2
viewoftheprecisionairdropimpromentfour2powerslongtermtechnologyproject[C].17thAIAAAerody2namicDeceleratorSystemsTechnologyConferenceandSeminar,2003.
[3]KEPeng,YANGChun2xin,YANGXue2song.Extraction
PhaseSimulationofCargoAirdropSystem[J].ChineseJournalofAeronautics,2006,19(4:
3152321.
(下转第130页
科学计算与信息处理王景志等:
基于龙格2库塔算法对降落伞空投的研究
图9 PID控制器硬件仿真
分析:
由于机器人的3个全向轮所承受的负载不一
样,即在相同的占空比的PWM下,3个轮子的实际速度并不相同,这就使得三轮速度不可能准确合成机器人的速度,进而影响机器人的控制轨迹。
根据图1所示的机器人1号轮和3号轮负载相当,2号轮子承受的负载较大,没有加入PID控制器时,前后运动虽然在一定范围内近似直线,但是机器人运行的速度达不到预期设定的速度,左右运动轨迹就是一个圆,而且设定的机器人左右移动速度大小还决定了机器人是朝顺时针方向还是逆时针方向转圈。
加入PID控制算法后,轮子的速度得到校正,机器人能够以预期设定的速度前后左右运动,特别是左右运动在一定范围内近似为直线,不再是圆圈。
可见PID闭环控制算法明显提高了机器人的控制性能。
5 结 语
针对目前常见的以DSP为核心实现足球机器人底层运动控制系统的方案,提出了一种采用FPGA实现三轮全方位移动足球机器人的底层运动控制系统的方法。
通过在三轮足球机器人上的应用实践,发现这种采用FPGA实现的方案有很好的实时性,精确度较高,而且由于FPGA本身的引脚多特点,其可扩展性较强,比如可以通过串口配置数字罗盘等外围信息传感器等其他传感器,同时,本设计对于研究多电机的机器人运动控制系统的实现方案有重要的参考价值和实用价值。
另外,由于全向轮的随动性较强,且易打滑,在实行精确控制的时候方向容易受到影响,而且PID闭环控制算法反应时间较长,参数还需要更多时间的调试,在以后的研究中,我们将研究更为精确的控制算法,实现对机器人的精确控制。
参 考 文 献
[1]TUKuo2yang,LUOShi2chao.Designandimplementationof
omni2directionalsoccerrobotsforroboCup[J].IEEEInterna2tionalConferenceonSystem,ManandCybernetics,2006,10:
200022005.
[2]黄永贤.全自主足球机器人运动控制系统的研制[D].广州:
广东工业大学,2008.
[3]冷春涛,曹其新,黄彦文,等.基于SoC的全方位移动机器人
运动控制系统设计与实现[C].2005年中国智能自动化会议论文集,2005.
[4]聂晓璐.全方位小型足球机器人运动性能研究[D].天津:
天
津大学,2006.
[5]CHARITKHUANChirot,BHURIPANYOJanjai,CHOO2
MUANGRerngwut.FPGAimplementationofclosed2loopcontrolsystemforsmall2sizedroboCup[J].IEEEConferenceonRobotics,AutomationandMechatronics,2006,6:
126.[6]罗真,曹其新,陈卫东.中型组自主式足球机器人平台设计
回顾[J].机器人,2003,25(4:
3782384.
[7]李拥军,杨文淑.光电编码器测速算法的IP核实现[J].长春
理工大学学报:
自然科学版,2008,31(3:
35237.
[8]李永新.RoboCup小型组足球机器人结构设计与分析[D].
北京:
中国科技大学,2003.
[9]杨文超,张荣标,张春艳.基于FPGA的PID控制器的研究
[J].微计算机信息,2005,21(11:
1352137.
[10]ChakravarthyNarashiman,XIAOJi2zhong.FPGA2basedcon2
trolsystemforminiaturerobots[J].Proceedingsofthe2006IEEE/RSJInternationalConferenceonIntelligentRobotsandSystem,2005,10:
339923404.
[11]刘波.精通VerilogHDL语言编程[M].北京:
电子工业出
版社,2007.
(上接第126页
[4]王利荣.降落伞理论与应用[M].北京:
宇航出版社,1997.[5]龚超,宋玉琴.无人机自动起飞地面滑跑技术的研究与仿真
[J].现代电子技术,2008,31(22:
1362138.
[6]李士明,曹凯.车辆导航定位仿真试验系统研究[J].现代电
子技术,2009,32(5:
1912194.
[7]吴训忠,胥勇军,任培林.一种基于Matlab的虚拟仪器仿真
和应用新途径[J].现代电子技术,2005,28(12:
57259.
[8]王昊鹏,刘兵.飞行器虚拟现实仿真研究[J].现代电子技
术,2008,31(13:
46248.
[9]高晓光.航空军用飞行器导论[M].西安:
西北工业大学出
版社,2002.
[10]王扬.弹射跳伞全过程仿真软件:
仿真支撑软件PROSIMS
应用示例[J].计算机仿真,2003,20(9:
125.
[11]朱勐,庄毅.降落伞投放实验仿真系统的研究与实现[J].
计算机仿真,2004,20(9:
24227.
作者简介:
王景志 男,1984年出生,吉林农安人,硕士研究生。
主要研究领域为复杂系统建模、仿真与评估。
陈红林 女,1951年出生,江苏南通人,教授,硕士研究生导师。
主要研究领域为复杂系统建模、仿真与评估。
科学计算与信息处理王全州等:
中型组机器人运动控制系统的FPGA设计
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 龙格库塔 算法 降落伞 空投 研究 概要