高超声速边界问题求解

最新推荐文章于 2024-10-07 12:07:10 发布

qq_42706301

最新推荐文章于 2024-10-07 12:07:10 发布

阅读量567

点赞数

文章标签：数学建模

本文链接：https://blog.csdn.net/qq_42706301/article/details/109457770

版权

aircraft 专栏收录该内容

1 篇文章 2 订阅

订阅专栏

高超声速滑翔式飞行器边界条件求解

参考文献《高超声速滑翔飞行器运动特性分析及弹道跟踪预报方法研究》

一、方程介绍

在这里插入图片描述

代入约束条件得到最终的数学模型为：
在这里插入图片描述

二、参数选取

%% the initial value
m = 907;            
Sm = 0.4839;        
Kn = 2.582e-7;     
Qmax = 6000;       
qmax = 45000;       
nmax = 4;           
alpha_max = 20;    
g0 = 9.8066;      
rho0 = 1.225;      
beta = 1/7110;    
R = 6371000;

气动系数自行下载该文献进行查找

三、平衡滑翔边界的迭代求解

在这里插入图片描述
其他三个方程均可以通过V表示H，因此很容易求解；

由于上式无法通过V表示高度H，因此采用迭代求解，而相关隐式方程较难写成迭代格式，写成迭代格式之后判断收敛性也是一个较难解决的难题。

在这里，本文采用穷举法，通过H的穷举找到满足条件的H，从而解决了给定一个V，求得一个H。具体如下：

for h=20000:1:80000
        con = exp(-beta*h);
        L=CL*rho0*Sm*v^2*con/(2*m);                         % 加速度
        test_pingheng = (g0*R^2)/((R+h)^2)-v^2/(R+h)-L;     % 式2.10取等号
        if abs(test_pingheng)<0.001
            v_end(1,i)=v;
            h_pingheng(1,i)=h;
            break;
        end
    end