《降落伞拉直阶段轨迹及拉直力计算》由会员分享，可在线阅读，更多相关《降落伞拉直阶段轨迹及拉直力计算（15页珍藏版）》请在金锄头文库上搜索。
1、1合理假设为简化计算，假设：1）拉伞过程中,引导伞与回收物运动轨迹为一条直线，且物伞系统作平面运动。)不考虑风的影响，物伞系统没有升力。3)不考虑伞绳的弹性。4)引导伞、物体和拉直中的伞系统微元质量作为三个质点解决。2动力学建模对于回收物及已拉出的主伞系统在回收物的坐标上的运动方程 ()其中,分别为回收物和主伞系统已拉出部分的气动阻力;为伞绳内的张力,最大值即为拉直力。对于拉直中的主伞系统微元质量可写成如下的动量方程 (2)其中，为伞系统微元质量拉出速度。 (3)为拉出的微元质量运动速度,在不考虑伞绳弹性状况下,。固然在方程(2)中忽视了主伞系统微元的自身重力。为伞系统拉出长度。为引导伞速度。
2、因此,（2)式可改写为 (4)其中,为拉动的伞系统单位长度质量。对于引导伞及伞袋在回收物坐标上的运动方程 （5）其中，为未拉出的主伞系统质量，,分别为引导伞及伞袋的阻力。为了消除()式中，将()式代入(1）式可得 (6）由于不考虑伞绳弹性，上式可改写为 ()又知 (8)拉直阶段的轨迹计算可以根据式(），(5)，()和(8）六个一阶微分方程组,已时间为计算步长,采用龙格库塔法算出各个时间内的个参数。对于先拉伞绳法 (9)其中,为伞绳及伞衣总质量，为伞衣底边沿伞绳方向的单位长度的质量，为伞衣底边宽度，为伞绳长度,为伞系统全长,为伞衣名义直径。3计算求解 MATAB中内置求解微分方程的ode函数,现
3、选用ode45进行求解计算。计算输入参数如下:ro1.23; 空气密度g.8; %重力加速度CAys=1.5; %引导伞阻力系数CAd=.2; %伞衣套阻力系数CAw=1； %物体阻力系数CAe=06; %已拉出物体阻力系数Qs=05h（)2*CAs； 引导伞气动阻力Q=0.5rho*x()2*CAd； %伞衣套气动阻力Qw=.5rho*x(5)2*CAw; %物体的气动阻力Q=0.5*hox(5)CAe; %已拉出伞系统的气动阻力Mw=60； %物体质量Mys=； %引导伞质量(涉及伞衣、套伞包)Ms=0.6; %伞绳的质量密度Ms=5; %伞边的质量密度b=0.5; %伞边的宽度s=6;
4、伞绳总质量My10; 伞衣总质量L=10; %伞绳全长Lxt=1405； %伞系统全长D=8; %伞衣名义直径Fs=50; %拉出阻力4成果图 Xd-时间历程图图2 Yd-时间历程图图角度时间历程图图4Vs-时间历程图图5 V-时间历程图图-时间历程图图7 FL-时间历程图图 轨迹曲线图5成果优化在前面的讨论计算中,假设伞绳为非弹性体，无伸长。显然，伞绳假设为非弹性体与实际状况出入较大。然而,由于在拉直过程中,除了伞绳拉直瞬间拉直力较大外，在拉动伞绳时伞绳内张力较小,因此,伞绳假设为非弹性体对计算拉直阶段轨迹来说影响很小,但对于拉直力的计算影响较大,因而在拉直力的计算中必须要考虑伞绳的弹性。假
5、设伞绳为弹性的无质量的弹簧，则弹性纵波沿伞绳的传播速度 (10)其中,为伞绳的理论弹性模量,为单根伞绳的单位长度质量。现假设传播速度为常数，则伞绳张力传播是一种典型的一维波动方程 (11)其中，为伞绳任意横剖面的位移，为伞绳剖面的坐标。为了求解式(1），必须拟定伞绳拉直瞬间的初始条件和边界条件。初始条件可以写成 (12）边界条件 (3)假设伞绳的弹性变形符合虎克定律,则 (1)由初始条件和边界条件可以解出式(1),得 （15）将式(15)对微分，并代入式（14)，可得 （16)将式(4)与式（16)联立求解，可以求得拉出速度 （) 其中，将式()代入式(16),可以求得伞绳拉直力 (1)6拉直
6、力计算 选用伞绳材料为凯夫拉29，弹性模量为62Pa,伞绳数量为0,其他参数参照不考虑伞绳弹性时的计算输入输出。运用MALAB可以简便快捷的求解计算成果,则。求解程序见附录。附录A主函数clercl,=de4(aero，0:0.0:1,0;0;15;10；100;0);输入初值%参数rho=1.93; %空气密度g=98; %重力加速度CAys1.5; 引导伞阻力系数Cd=0.2； %伞衣套阻力系数CAw=； %物体阻力系数CA=0; %已拉出物体阻力系数Qys=0.rho*(4)2*CAys; %引导伞气动阻力Qd=0.x(4)*C； %伞衣套气动阻力Qw0.5*ro*x（)2*CAw; %
7、物体的气动阻力Q=0.5*rhox(5)2*CAe; 已拉出伞系统的气动阻力Mw=60； 物体质量My=5; 引导伞质量（涉及伞衣、套伞包)sh=0.; 伞绳的质量密度Msy1=5; %伞边的质量密度=00; %伞边的宽度Msh; %伞绳总质量Ms=1; 伞衣总质量Lsh=1; %伞绳全长=4.05; %伞系统全长o=8; %伞衣名义直径sh0; %拉出阻力mlngth()；Me=zeros(m,1);Dml=zeros(,1）;F=zer(m,1）；ifx(6）=sh Me=hx(6)；l (6)Lx Me=Msh1*Lh+Ms*b+(MsMs*)*(1*(Lt-x（6）Do)2);endM
8、MshMy+Msy*b-e; %未拉出主伞的质量 i=:m;if x(,6)Lsh Dml(）=h1;elsei (i,6）=Lh+ Dl（i）=1;elseif(i,6)=L Dml(i)=-1.25*Mh*x（i，)+17.5Msh1; nedk=mnu(拉直阶段各变量时间历程曲线：,Xd-时间历程图, d时间历程图,角度时间历程图,Vys-时间历程图,w-时间历程图,L-时间历程图,F-时间历程图，轨迹曲线图）；switch k case 1 figre(1); plo(t,x(:,）); xlbl(t/s)； ylael(X/m); ite (X-时间图）; cas 2 figur(2
9、); pl（,x(:,); xlel(/s)； lel（Yd/）; tile(Yd-时间图）; case 3 igure(3); plo(,x（:,3）; xlabe（t/s）; ylbel(角度)； til(角度-时间图) ase4 figure(); lot(,(:,4）; xlabel(t/s); yel(Vys/m); ttl(Vy-时间图) cae 5 figue(5) plot（,x(:,5）)； xlabel(t/s); ylael(m/); itle(Vw-时间图) cse gur(6) plot(t,x(:,6)）; abel(t/); bel(Lm)； title(L-时间图) case 7 for i1:m; F()=Dm()*（(i,5)-x(i,4)+Fsh; e fgure（)； plo(t,FL); xlabel(ts); ylbl(L/N）；
