基于弯曲激波理论的全三维乘波体反设计方法

技术领域

本发明涉及临近空间高超声速乘波体领域，尤其涉及基于弯曲激波理论的全三维乘波体反设计方法。

背景技术

临近空间飞行器的发展涉及国家安全与和平利用空间，是目前国际竞相争夺空间技术的焦点之一。以美国、俄罗斯为代表的世界强国都在大力推进各自的高超声速飞行研制计划(Joseph,M.H,James S.M.Richard C.M.,The X-51A Scramjet Engine FlightDemonstration Program,15

虽然在高超声速乘波体研究领域，各项研究已经取得了显著的进展，部件性能也在不断提升。然而，迄今为止，即使吻切乘波体也仍旧采用轴对称基准流场，进而利用流线追踪技术获得对应型面。但是基于轴对称基准流场得到的三维乘波体不存在横向流动，仅有流向流动。目前，科研人员尚未发现有效的方法获得全三维基准流场，进而得到具有全三维流动的三维乘波体。而横向流动对飞行器性能提升有至关重要的作用。与此同时，科研人员普遍采用传统特征线法进行基准流场反设计，不仅编程复杂，而且稳定性差，限制了基本流场选择范围，进而减小了乘波体的几何构造范围。由此可见，目前制约高超声速乘波体性能的问题之一是缺乏一种基于弯曲激波理论的全三维乘波体反设计方法。

发明内容

本发明的目的在于解决现有技术中的上述问题，提供基于弯曲激波理论的全三维乘波体反设计方法。

为达到上述目的，本发明采用如下技术方案：

基于弯曲激波理论的全三维乘波体反设计方法，包括以下步骤：

1)根据设计要求指定全三维基准流场内全三维外流激波，全三维外流激波采用椭圆或非轴对称形状设计；

2)将全三维外流激波离散为一系列参考平面，根据全三维外流激波角、激波曲率以及波后参数，利用弯曲激波理论求解对应的全三维外压缩基本流场；

3)设计全三维乘波体捕获型线并在步骤2)的全三维外压缩基本流场中进行流线追踪得到高超声速全三维乘波体下表面压缩型面，所述型线采用椭圆曲线或非轴对称形状；

4)以下表面压缩型面为基础对高超声速全三维乘波体进行几何构造：根据容积率要求，将乘波体上表面构造为平面或凸面，获得在设计飞行状态下基于弯曲激波理论的全三维乘波体。

所述非轴对称形状包括二次曲线。

步骤2)中，弯曲激波理论的控制方程如下所示：

其中，p为压力，δ为流动角，μ为马赫角，ρ为密度，V为流向速度，w为周向速度，γ为比热比，j为判断因子，s为流线，l为特征线，φ为周向角度，P为沿流线压力的导数，D为气流角沿流线的导数，y为笛卡尔坐标系内纵坐标，a为当地声速。

在步骤3)中，首先在全三维外压缩流场内，利用结尾捕获型线上每一点坐标(x，y)，求解得到不同参考平面内有效流线上每一点的三维坐标(x，y，z)，之后将有效流线上每一点三维坐标组合得到压缩型线，不同参考平面内压缩型线组合构成全三维乘波体下表面压缩型面。

一种基于弯曲激波理论反设计法的全三维乘波体，所述全三维乘波体设有全三维乘波体下表面压缩型面、全三维乘波体上表面型面、全三维乘波体前缘捕获型线；所述全三维乘波体下表面压缩型面于全三维乘波体前缘捕获型线处进入全三维乘波体上表面型面。

所述全三维乘波体上表面型面构造为平面或凸面。

所述全三维乘波体下表面压缩型面为非轴对称曲面。

所述全三维乘波体呈左右对称布置。

相对于现有技术，本发明技术方案取得的有益效果是：

利用本发明生成的基于弯曲激波理论的全三维乘波体同时兼顾三维乘波体横向流动与三维乘波体的高升阻比，并且实现已知激波的全三维流动反设计。所述全三维乘波体为全三维带横向流动的乘波体，增大升力的同时减小外流阻力；且横向流动可帮助低能流向两侧移动，高能流向中心汇聚，提高飞行器性能。此外，弯曲激波理论相较于传统特征线法计算效率更高，精度更高，获得的基本流场范围更广。全三维乘波体通过合理设计横向流动拓宽了三维乘波体的设计范围，从而进一步提高飞行器工作性能。

附图说明

图1为基于弯曲激波理论的全三维乘波体方案基准流场示意图。

图2为弯曲激波理论的求解示意图。

图3为基于弯曲激波理论的全三维乘波体方案仰视示意图。

图4为基于弯曲激波理论的全三维乘波体方案半剖结构示意图。

图5为基于弯曲激波理论的全三维乘波体方案总体结构示意图。

具体实施方式

为了使本发明所要解决的技术问题、技术方案及有益效果更加清楚、明白，以下结合附图和实施例，对本发明做进一步详细说明。

如图1～5所示，基于弯曲激波理论的全三维乘波体反设计方法，包括以下步骤：

1、根据设计要求指定全三维基准流场内全三维外流激波1，全三维外流激波1采用椭圆或其他非轴对称形状设计；其中，由于全三维基准流场的非轴对称性，不同参考平面6内的全三维外流激波1的型线不相同，同理，待求解的不同参考平面6内的压缩型线2也不相同。

2、将全三维外流激波1离散为一系列参考平面，根据全三维外流激波角、激波曲率以及波后参数，利用弯曲激波理论求解对应的全三维外压缩基本流场；

具体地，求解过程在不同参考平面6内进行，如图2所示，根据全三维外流激波1上离散点激波角、激波曲率以及波后参数和全三维基本流场中心线3，利用弯曲激波理论求解得到不同参考平面6内压缩型线7，将所有参考平面6内压缩型线7组合得到对应的全三维外压缩流场；

弯曲激波理论的控制方程如下所示：

其中，p为压力，δ为流动角，μ为马赫角，ρ为密度，V为流向速度，w为周向速度，γ为比热比，j为判断因子，s为流线，l为特征线，φ为周向角度，P为沿流线压力的导数，D为气流角沿流线的导数，y为笛卡尔坐标系内纵坐标，a为当地声速。

3、设计全三维乘波体捕获型线并在步骤2的全三维外压缩基本流场中进行流线追踪得到高超声速全三维乘波体下表面压缩型面，所述型线采用椭圆曲线或二次曲线等非轴对称形状；

具体地，首先在全三维外压缩流场内，利用结尾捕获型线5上每一点坐标(x，y)，求解得到不同参考平面6内有效流线4上每一点的三维坐标(x，y，z)，之后将有效流线4上每一点三维坐标组合得到压缩型线7。不同参考平面6内压缩型线7组合构成全三维乘波体下表面压缩型面9。

4、以下表面压缩型面为基础对高超声速全三维乘波体进行几何构造：将乘波体上表面11构造为平面或凸面。获得在设计飞行状态下基于弯曲激波理论的全三维乘波体。为了保证飞行过程中力矩平衡，全三维乘波体以对称面8呈左右对称布置。

本发明基于弯曲激波理论设计的全三维乘波体由全三维乘波体下表面压缩型面9、全三维乘波体上表面型面11和全三维乘波体前缘捕获型线10组成；且三维乘波体能够实现外部乘波。全三维乘波体下表面压缩型面9于全三维乘波体前缘捕获型线10处进入全三维乘波体上表面压缩型面11。所述全三维乘波体上表面型面11构造为平面或凸面。所述全三维乘波体下表面压缩型面9为非轴对称曲面。

本发明基于弯曲激波理论的全三维乘波体方案在保持三维乘波体高升阻比优点的同时，实现三维乘波体外部全三维流动。通过合理设计全三维基准流场可以利用气流的三维横向流动提高乘波体的工作性能，进而增加飞行器的总体性能。此外，全三维基准流场拓宽乘波体设计范围，为飞行器构造匹配提供更多的选择。