基于MATLAB降落伞拉直过程性能教程分析.docVIP

防护与救生技术 降落伞拉直过程性能分析 姓名：WXH 班级： 学号： 学院：能源与动力学院 一、拉直阶段假设 为简化计算，假设： 拉伞过程中，引导伞、物体运动轨迹为一条直线，物-----伞系统作平面运动。 不考虑风的影响，物-----伞系统没有升力。 在拉直过程中，伞绳为非弹性体，无伸长。 引导伞、物体和拉直中的伞系统微元质量dm作为三个质点处理。 此次仿真采用倒拉法进行性能仿真。 拉直阶段计算内容 采用MATLAB软件编写程序，利用已知方程组推导出拉直过程中各个参数随时间的变化，并利用MATLAB输出曲线图像，再利用曲线图像对整个过程进行过程性能分析。 基本方程组： 一共有七个方程，其中六个为微分方程，均为变量对时间t的导数，故采用MATLAB程序编写简介迅速，且输出图像曲线简单明了。 编写程序确定参数 1、编写程序 MATLAB中有专门解决一阶微分方程的ode45（龙格-库塔函数），ode函数一共有其中，在这里我采用了ode45函数，可以自适应变步长的求解方法，从而使计算速度很快，并且计算精度较高。 2、确定参数 由于之前从来没有了解过有关降落伞的知识，所以在这一方面感觉很欠缺，没有什么概念，因此在确定参数的时候浪费了很多时间，参数选取的不准备会使图像输出有很大的出入，因此选参数的时候应该仔细小心。在这次的作业中，对于主伞，我采取了主伞伞衣和主伞伞绳质量相等即比例为1:1的参数，从而确定了主伞质量，伞绳质量密度，伞衣质量密度的分布等。因此，此次仿真都是基于伞衣质量=伞绳质量的条件下进行仿真。其余参数经过请教老师也基本确定。 p=1.2; %空气密度 g=9.81; %重力加速度 CAys=1.5; %引导伞阻力系数 CAd=0.2; %伞衣套阻力系数 CAw=1; %物体阻力系数 CAe=0.6; %已拉出物体阻力系数 Qys=0.5*p*x(4)^2*CAys; %引导伞气动阻力 Qd=0.5*p*x(4)^2*CAd; %伞衣套气动阻力 Qw=0.5*p*x(5)^2*CAw; %物体的气动阻力 Qe=0.5*p*x(5)^2*CAe; %已拉出伞系统的气动阻力 Mw=60; %物体质量 Mys=5; %引导伞质量(包括伞衣、套伞包) Msh1=0.6; %伞绳的质量密度 Msy1=5; %伞边的质量密度 b=0.05; %伞边的宽度 Msh=6; %伞绳总质量 Msy=6; %伞衣总质量 Lsh=10; %伞绳全长 Lxt=14.05; %伞系统全长 Do=8; %伞衣名义直径 Fsh=50; 输出图像 图像分析 由图像可以看出，除了拉直的图像外，其余图像都是平滑的曲线，Xd和Yd的图像随时间的变化均增加，Xd之所以为负值是建立坐标系的关系，但它们的数值都是逐渐变大的。 整个过程时间很短，在这次仿真中，以上图像都是以80m/s的开伞速度下进行仿真，从L-t图像曲线可以看出，大概只需要0.65s左右就可以拉出整个主伞。 Vw和Vys的