《基于激波装配法的乘波体设计与分析-陈冰雁.pdf》由会员分享，可在线阅读，更多相关《基于激波装配法的乘波体设计与分析-陈冰雁.pdf（8页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、 第3 5卷 第3期空气动力学学报V o l . 3 5,N o . 3 2 0 1 7年0 6月ACTAAERODYNAMICASINICAJ u n e , 2 0 1 7 文章编号: 0 2 5 8 - 1 8 2 5(2 0 1 7)0 3 - 0 4 2 1 - 0 8基于激波装配法的乘波体设计与分析陈冰雁,刘传振*,纪楚群(中国航天空气动力技术研究院,北京 1 0 0 0 7 4 )摘 要:乘波体外形通过在给定的高超声速激波流场中使用流线追踪法设计得到,传统的乘波体设计一般采用二维楔形或者轴对称锥形流场生成,设计空间受限。为拓展乘波体设计空间,引入激波装配法,数值计算带激波的流场,
2、建立了普适于一般三维流场的乘波体设计方法,突破了传统乘波体设计方法的流场限制。引入“导波体”定义三维流场的生成外形,采用平切导波体流场激波及正向追踪流线的方法设计乘波体外形。通过几类导波体流场生成乘波体的算例,分析了导波体与乘波体之间的关系,包括纵向截面外形、横截面外形以及前缘外形等。分析表明,双锥体在后体产生的激波或膨胀波能有效地改变导波体速度沿纵向的分布特性,利用该特性可在俯视平面形状受约束的情况下实现对乘波体构型纵向稳定性的按需设计;通过改变导波体横截面外形,可以改变速度沿展向的分布特性,从而实现乘波体横截面外形改变;乘波体前缘外形主要与导波体激波在俯视面上的外形相关,通过后体扩张的双锥
3、外形能够生成双前缘乘波体外形。本文指出乘波体与导波体外形之间存在定性对应关系,其关系可作为乘波体设计时导波体选择的参考依据,为扩大设计空间、设计优良乘波体奠定了基础。关键词:激波装配法;三维流场;乘波体;气动外形;数值模拟中图分类号:V 4 1 1 . 4 文献标识码:A doi: 1 0 . 7 6 3 8 / k q d l x x b - 2 0 1 7 . 0 0 3 9收稿日期:2 0 1 7 - 0 3 - 0 6 ; 修订日期:2 0 1 7 - 0 3 - 2 7基金项目:自然科学基金( 1 1 6 7 2 2 8 1 )作者简介:陈冰雁( 1 9 7 7 - ) ,男,广东惠
4、州人,研究员,研究方向:飞行器气动设计. E - m a i l : c h e n _b i n y a n h o t m a i l . c o m通信作者:刘传振* ( 1 9 8 9 - ) ,男,山东德州人,博士,研究方向:气动外形设计. E - m a i l : c h u a n z h e n l i u 1 2 6 . c o m引用格式:陈冰雁,刘传振,纪楚群.基于激波装配法的乘波体设计与分析 J .空气动力学学报, 2 0 1 7 , 3 5 ( 3 ) : 4 2 1 - 4 2 8 .d o i : 1 0 . 7 6 3 8 / k q d l x x b -
5、2 0 1 7 . 0 0 3 9 C H E N B Y , L I U C Z , J I C Q . W a v e r i d e r d e s i g n a n d a n a l y s i s b a s e d o n s h o c k - f i t t i n g m e t h o d J . A c t a A e r o d y n a m i c a S i n i c a , 2 0 1 7 , 3 5 ( 3 ) : 4 2 1 - 4 2 8 .Waveriderdesignandanalysisbasedonshock-fittingmethodC H
6、 E N B i n g y a n , L I U C h u a n z h e n * , J I C h u q u n(ChinaAcademyofAerospaceAerodynamics,Beijing 1 0 0 0 7 4 ,China)Abstract: M o s t c o m m o n l y , a w a v e r i d e r i s d e f i n e d a s a n a e r o d y n a m i c c o n f i g u r a t i o n t h a t i si n v e r s e l y c r e a t e d
7、 f r o m a p r e s c r i b e d h y p e r s o n i c f l o w f i e l d b a s e d o n a p l a n a r o r c o n i c a l s h o c k w a v e .I n t h i s w o r k , t h e p r e s c r i b e d s h o c k w a v e w a s n u m e r i c a l l y c r e a t e d f r o m a n a r b i t r a r y c o n f i g u r a t i o nu s
8、 i n g s h o c k - f i t t i n g m e t h o d . A S h o c k G e n e r a t i n g B o d y ( S G B ) w a s d e f i n e d a s t h e o b j e c t t og e n e r a t e 3 D f l o w f i e l d w i t h s h o c k w a v e f o r w a v e r i d e r c r e a t i o n . T h e c o r r e l a t i o n o f t h e l o n g i t u
9、d i n a ls e c t i o n , c r o s s s e c t i o n , a n d l e a d i n g e d g e f e a t u r e b e t w e e n t h e w a v e r i d e r c o n f i g u r a t i o n a n d i t sc o r r e s p o n d i n g S G B c o n f i g u r a t i o n w e r e a n a l y z e d . T h e a n a l y s i s r e s u l t s s h o w t h
10、a t t h e s h o c k w a v eo r e x p a n s i o n w a v e g e n e r a t e d f r o m t h e a f t - b o d y o f a d o u b l e - c o n e S G B c a n e f f e c t i v e l y a l t e r t h ev e l o c i t y d i s t r i b u t i o n c h a r a c t e r i s t i c s o f t h e S G B f l o w f i e l d . T h i s f e
11、a t u r e c a n b e u s e d t o d e s i g n t h el o n g i t u d i n a l c h a r a c t e r i s t i c s o f a w a v e r i d e r c o n f i g u r a t i o n w h e n i t s p l a n a r p r o f i l e i s r e s t r i c t e d . T h ev e l o c i t y d i s t r i b u t i o n o n c r o s s s e c t i o n c a n b
12、e a l t e r e d b y c h a n g i n g t h e c r o s s s e c t i o n p r o f i l e o f t h eS G B , t h e r e b y a c h i e v i n g t h e c r o s s s e c t i o n p r o f i l e a l t e r a t i o n o f t h e g e n e r a t e d w a v e r i d e r . Aw a v e r i d e r w i t h d o u b l e - s w e e p l e a d
13、i n g e d g e f e a t u r e c a n b e g e n e r a t e d u s i n g a d o u b l e - c o n e w i t he x p a n d i n g a f t - b o d y . Q u a l i t a t i v e c o r r e l a t i o n s e x i s t b e t w e e n t h e w a v e r i d e r a n d t h e S G B , w h i c hp r o v i d e s v a l u a b l e i n s i g h
14、t a s t h e g u i d e l i n e t o c h o o s e a p p r o p r i a t e f l o w f i e l d f o r t h e d e s i g n o f h i g hp e r f o r m a n c e w a v e r i d e r s .Keywords: s h o c k - f i t t i n g m e t h o d ; t h r e e - d i m e n s i o n a l f l o w f i e l d ; w a v e r i d e r ; a e r o d y
15、n a m i cs h a p e ; n u m e r i c a l s i m u l a t i o n万方数据0 引 言乘波体设计是当前国际上高超声速飞行器气动布局研究的重点和热点之一。乘波体构型在设计状态将激波附着于下表面前缘,阻止了气流泄露,具有很高的升阻比。升阻比是反应高超声速飞行器气动特性非常重要的一个参数,升阻比很大程度上决定了高超声速飞行器能够实现哪一种飞行轨道。对于滑翔式和巡航式高超声速飞行器的气动设计,升阻比往往是非常重要的一个技术指标,因为升阻比直接关系到飞行器能够达到的航程、横向机动能力等关键战技术指标。因此,乘波体往往被选择为高超声速飞行器气动布局的基础构型
16、。乘波体的设计要素包括两点:带激波的流场与沿流线追踪生成型面。本文定义生成激波流场的外形为“导波体” ,与在流场中追踪流线得到的外形“乘波体”相对应。传统的乘波体设计方法均可以据此定义导波体,例如楔形流场乘波体 1 的导波体为二维楔形,锥导乘波体 2 - 3 对应圆锥,密切锥方法 4 - 5 的流场是由已知激波形状反推得到,基准流场为锥形流,导波体也可以认为是锥体的一种。基准流场对乘波体的性能有着根本性的影响,选择导波体就是选择基准流场,是乘波体设计过程的重要环节,也是保证乘波体满足气动、装载等设计要求的基础。圆锥的激波流场一般通过求解T a y l o r - M a c c o l l方程
17、数值 6 ,在乘波体设计中应用较多,但这也限制了导波体的选择范围,设计空间不大,同时流场模拟不够准确 7 。为了规避传统乘波体设计方法的局限性,相关学者提出使用三维流场进行乘波体设计 8 - 9 。求解三维流场时,激波曲面的准确分辨是一个难题,M a r c u使用激波捕捉法计算了圆锥流场,并沿流线根据压力梯度确定激波位置,这一方法计算量较大,激波曲面的光滑性难以保证,推广到复杂外形存在困难。激波装配法 1 0 可以准确定位激波曲面,为寻找性能更优的流场进行乘波体设计提供了良好的工具。前期工作中笔者使用三维流线追踪技术 1 1 ,探索了三维流场乘波体快速设计方法,在扩大设计空间方面具有优势。在
18、三维流场乘波体设计中,随着设计空间的扩大,合理选择流场变得越来越重要。如果能找到乘波体与导波体外形间的对应关系,就可以先根据乘波体的要求寻找导波体,再用所得的导波体设计外形,从而提高乘波体设计的灵活性。本文使用激波装配法计算带激波的流场,固定流线追踪( F C T , F l o w C a p t u r eT u b e )初始线为水平直线,投影到激波作为流线追踪的起始点,即“平切”激波,正向追踪流线生成乘波体外形。研究几类导波体生成乘波体的算例,指出乘波体与导波体之间存在一定程度上的对应关系,可作为未来设计中导波体外形选择的参考依据。1 激波装配法带激波的超声速流场通常采用激波捕捉法数值
19、求解。激波捕捉法无需将激波作为未知边界进行求解,因此所得到的激波边界实际上是一个流动量连续急剧变化的狭小区域,给激波的分辨带来了困难,难以用于乘波体设计。与激波捕捉法不同,激波装配法则将流场的外激波作为未知外边界同流场一起求解(图1 ) 。为此引入三维非定常激波的求解方程,包括R a n k i n e - H u g o n i o t激波关系式及描述激波边界处特征波传播特性的特征相容关系式。将上述方程联立数值求解,可求得波后流动量及激波速度。由于外激波边界在求解过程中是移动的,因此求解过程需在动网格中进行。当方程解收敛时,激波速度趋近于零,非定常激波收敛为定常激波,即得到流场外激波的准确外
20、形,所得激波面可直接用于乘波体设计。图1 激波捕捉法及装配法比较Fig.1 Comparisonbetweenshock-capturingmethodandshock-fittingmethod采用超声速E u l e r方程激波装配法数值计算程序求解导波体激波及流场。选择文献 1 2 所给物形作为测试算例,其中Ma= 6 ,= 1 5 ,图2为计算得图2 激波外形计算结果Fig.2 Shocksimulationresultsusingshock-fittingmethod224空 气 动 力 学 学 报 第3 5卷万方数据到的激波,图3和图4分别为激波位置rs和激波斜率rsx与文献 1
21、2 数据的比较,可以看出两者非常吻合。图3 激波位置计算结果比较Fig.3 Comparisonofshocklocations图4 激波斜率计算结果比较Fig.4 Comparisonofshockwaveslope2 导波体2.1 外形基本特征描述方法在描述乘波体和导波体外形时采用如图5所示的坐标系,其中Ux、Uy、Uz表示流场速度分量。图5 坐标系定义Fig.5 Definitionofcoordinatesystem在此坐标系下,乘波体和导波体由三个坐标面上的外形描述: 1 )xy方向上描述纵向截面外形; 2 )yz方向对应横截面外形; 3 )xz方向描述前缘外形。固定流线追踪初始线时
22、,决定乘波体外形基本特征的因素包括:纵向截面外形与导波体流场速度Uy/Ux沿纵向(x向)分布特性相关;横截面与Uy/Ux、Uz/Ux沿展向(z向)的分布特性相关;前缘外形及长宽比与激波形状和流线追踪初始线切割激波的位置等因素相关。因此乘波体外形的基本特征可以通过改变导波体的激波外形及速度分布得到。下面给出两类可改变激波外形及流场速度沿纵向及展向分布的导波体 锥形导波体和双锥导波体,作为重点研究对象。2.2 锥形导波体锥形导波体的外形为尖头,母线为直线,横截面可为任意外形。图6给出了几种典型的横截面形状,其中截面为圆的外形即为锥导乘波体中的圆形尖锥导波体。图6 锥形导波体的横截面外形Fig.6
23、CrosssectionsofconicalSGBs锥形流的流场特性与长度无关,即沿纵向分布的各横截面具有相同的流场特性。改变横截面外形,可以改变锥形导波体流场的速度Uy/Ux和Uz/Ux沿z向的分布特性,进而改变乘波体的横截面外形。2.3 双锥导波体双锥导波体的外形特征是尖头锥形前体、后接扩张或收缩的后锥体,横截面为任意外形。过前体的流场为锥形流场,当后体收缩时,存在沿纵向传播的膨胀波,后体扩张时,存在激波。通过后体出现的膨胀波或激波,可以改变流场速度Uy/Ux、Uz/Ux沿纵向的分布特性,从而改变乘波体纵向截面外形;俯视外形前后锥体的长度、角度决定激波外形,进而决定乘波体的前缘外形;改变横
24、截面外形,可以改变流场速度沿展向的分布特性,从而改变乘波体横截面外形。3 乘波体与导波体外形对应关系3.1 纵向截面外形对应关系分析乘波体外形通过追踪流线生成,纵向截面流场对于乘波体性能影响很大,目前常用的锥导乘波体等使用圆锥流动,流场分布较简单,不存在激波膨胀波等现象,但也限制了设计空间。而增加波系则可以使外形具有更好的特性,比如吕侦军等 1 3 提出了多级压缩锥导乘波体的概念,使用多波系充分发挥前体的预324第3期 陈冰雁等:基于激波装配法的乘波体设计与分析万方数据压缩作用,为进气道的正常工作提供均匀流场。本文使用准确的三维流场进行设计,可以方便地构造流场内激波或膨胀波,提高了设计灵活性及
25、精度。3.1.1 流场内波系对速度分布的影响导波体流场内存在激波或膨胀波时,由于波前波后速度的明显变化,流线穿过这些波时走向会发生明显改变。分析三个典型外形:尖锥、后体收缩的双锥、后体扩张的双锥,考察波系对速度的分布。图7给出了尖锥和后体收缩、扩张的双锥导波体流场内Uy/Ux的分布特性。图8给出了取定导波体流场内的同一位置截面a b ,查看速度分布的特性。可以看到收缩后体产生膨胀波,波后Ux大于波前值,而Uy小于波前值,因此膨胀波波后Uy/Ux小于波前值;扩张双锥内存在激波,激波波后速度分布与膨胀波相反,Uy/Ux大于波前值。( a )尖锥( b )后体收缩膨胀波( c )后体扩张激波图7 双
26、锥导波体流场的Uy/Ux分布特性Fig.7 Uy/Uxdistributioncharacteristicsofdouble-conicalflow根据波前波后速度分布Uy/Ux的变化特性,流线穿过内膨胀波及内激波的趋势可以大体确定,如果双锥导波体后体收缩,所生成的乘波体边界后段也向上收缩;双锥导波体后体扩张,所生成的乘波体边界后段也向下扩张。图8 双锥导波体流场ab截面的Uy/Ux分布特性Fig.8 Uy/Uxdistributioncharacteristicsatabcrosssection3.1.2 具体算例分析考虑三个导波体生成流场进行乘波体设计,设计状态为Ma= 6 ,攻角= 0
27、,其中外形1为8 . 5 半锥角的尖锥外形;外形2为后体收缩的双锥外形( 1 0 半锥角转7 ) ;外形3为后体扩张的双锥外形( 7 半锥角转1 0 ) 。横截面均为圆,图9所示为三个导波体外形流场的Uy/Ux云图。( a )锥形( b )收缩双锥( c )扩张双锥图9 导波体外形的流场速度Uy/Ux云图Fig.9 VelocityUy/UxcontourofSGBflowfield构造乘波体外形,下表面使用追踪流线得到,上表面则考虑扩充容积设计 1 4 ,截面形状用三次函数424空 气 动 力 学 学 报 第3 5卷万方数据表示,形状为“凸” 。在不同乘波体的构造过程中保持长宽比不变,得到的
28、乘波体外形为W 1 、 W 2和W 3 ,对应于图9中的导波体外形1 、 2和3 ,其纵向截面、横截面及前缘外形的比较如图1 0所示。可见三个乘波体的前缘外形和横截面外形均相差不大,主要区别是在纵向截面外形。图10 乘波体W1、W2和W3外形特征比较Fig.10 ConfigurationsofW1,W2andW3使用C F D方法计算乘波体的流场,图1 1给出了三个乘波体外形纵向截面和迎风面表面的压力分布云图。三个外形的压力分布差别很大: W 1的压力分布比较均匀; W 2由于膨胀波的影响后体表面压力明显低于前体; W 3则由于激波的影响后体表面压力明显高于前体。这种不同的压力分布特性对乘波
29、体的纵向稳定特性有着显著的影响。图11 乘波体压力云图对比Fig.11 Pressurecontourofthewaveriders图1 2给出了三个乘波体的升阻及纵向压心特性,即升力CLS(CL为升力系数,S为参考面积) ,升阻比L/D和纵向压心系数Xcp的比较。由图1 1给出的乘波体迎面压力云图可以看出,内激波的存在导致W 3的迎风面高压区较大,迎风面压力高于W 1 ,而膨胀波流场得到的外形迎风面压力低于W 1 ,因此W 3的升力最大, W 2升力最小。与升力特性相对应,在设计点(= 0 )附近, W 3的L/D最大而W 2最小,但是随着攻角增大, W 3的阻力迅速增加,导致其升阻比低于W
30、 1和W 2 。此外,由于W 3迎风面高压区集中在后部,导致压心靠后, W 2迎风面低压区集中在后部,导致压心靠前。( a )CLS( b )L/D( c )Xcp图12 乘波体气动性能比较Fig.12 Aerodynamicpropertiesviaangle-of-attack上述分析可知,双锥体在后体产生的激波或膨胀波能有效地改变导波体流场内速度沿纵向分布的特性,利用该特性可以设计不同纵向截面外形的乘波体,从而实现在俯视平面形状受约束的情况下对乘波体构型纵向稳定性的按需设计,而这种乘波体性能设计的灵活性正是常规乘波体设计方法所难以提供的。本节算例导波体采用了双锥轴对称外形,实际应用中只要
31、求物形在纵向截面的迎风面存在收缩或扩张即可。3.2 横截面外形对应关系分析乘波体的横截面外形与导波体流场内速度524第3期 陈冰雁等:基于激波装配法的乘波体设计与分析万方数据Uy/Ux、Uz/Ux沿展向的分布特性有关,此速度特性可由导波体外形的横截面确定。分析三个导波体外形:尖锥A 0 、双椭圆锥D 1 、双椭圆锥D 2 ,如图1 3所示,这三个导波体外形的纵向截面外形完全一致,横截面不同 A 0横截面为圆, D 1和D 2横截面为两个椭圆的组合, D 2的下部份外形(迎风面)比D 1分布较为平坦。( a ) A 0( b ) D 1( c ) D 2图13 双椭圆锥导波体与尖锥导波体Fig.
32、13 DoubleellipticalconeSGB图1 4给出了导波体A 0 、 D 1和D 2的横截面比较,并给出了迎风面处沿c d位置流场速度Uy/Ux、Uz/Ux沿展向分布的比较。导波体在相同迎风面锥角条件下,横截面迎风面越平,Uy/Ux的负值越大,而且在整个c d段都保持这一特性。Uz/Ux的分布特性有所不同,在c d前段,随着导波体迎风面逐渐变平,Uz/Ux逐渐变小,而在c d后端,Uz/Ux快速增大。上述特性导致乘波体横截面外形及气动特性差别明显。在要求长宽比相同的条件下,由导波体A 0 、 D 1和D 2生成的乘波体W - A 0 、 W - D 1和W - D 2的外形如图1
33、 5所示。三个乘波体的纵向截面外形和前缘外形基本不变,区别主要在横截面外形。由前面给出的速度沿展向分布特性可以得出,导波体迎风面越平,( a )三个导波体横截面及激波形状( b )速度Uy/Ux沿展向cd分布图( c )速度Uz/Ux沿展向cd分布图图14 导波体流场的横截面速度分布Fig.14 VelocityoncrosssectionofdifferentSGBs图15 乘波体三视图比较Fig.15 Waveriderconfigurationswithdifferentcrosssectionprofiles生成的乘波体越平,内侧厚度越大,而展向外缘则较薄。这种横截面特征的对应关系可以
34、在对乘波体外形的迎风面进行扩容设计时获得应用。3.3 前缘外形对应关系分析乘波体前缘由导波体流场的激波曲面决定,由于锥形流导波体的激波外形为锥形或近似锥形,采用F C T为直线平切激波曲面时,其生成的乘波体前缘形状一般为近似抛物线,如图1 6所示。而我们扩大设计空间,则可以得到前缘线具有特殊性质的乘波体外形,比如双后掠外形等,这为我们改善乘波体低速阶段的性能缺陷提供了一种新的思路。如果需要得到双后掠外形的乘波体,则需采用俯视图为双前缘的导波体。如图1 7所示的两个导波体外形, Y 1横截面为双椭圆, Y 2为椭圆,纵向截面外形均为锥形,锥角为8 ;俯视平面外形均为后体扩张的双锥,锥角分别为7
35、和1 8 。624空 气 动 力 学 学 报 第3 5卷万方数据图16 典型类锥导流场生成的乘波体外形Fig.16 Conicalflowwaverider图17 双前缘导波体外形Fig.17 DoublesweptconicalSGB从图1 8看到,导波体Y 1 、 Y 2的激波的俯视图为扩张的双前缘形状,可生成前缘为双后掠的乘波体外形,同时由于Y 1 、 Y 2的横截面不同,生成的乘波体W Y 1的横截面比W Y 2更平(图1 9所示) 。图2 0为乘波体W Y 1和W Y 2气动特性的比较。( a )侧视图 ( b )俯视图( c ) Y 1正视图 ( d ) Y 2正视图( e ) W
36、 Y 1 ( f ) W Y 2图18 双锥导波体激波流场及乘波体外形Fig.18 ShockwavefromdoubleconicalSGBs图19 双后掠乘波体外形Fig.19 Viewsofthedoublesweptwaverider( a )CLS( b )L/D( c )XcpS图20 双后掠乘波体气动特性对比Fig.20 Aerodynamicperformancesofthedoublesweptwaverider724第3期 陈冰雁等:基于激波装配法的乘波体设计与分析万方数据与其他乘波体类似,升力呈线性增长,升阻比在= 4 左右取得最大,同时压心靠后,并随攻角增大压心逐步前移
37、。虽然W Y 1和W Y 2的纵向对称面外形及前缘外形相差不大,但由于W Y 1横截面下表面更平,横截面积比W Y 2大,导致两者气动特性不同。4 结 论通过引入激波装配法建立了普适于一般三维流场的乘波体设计方法,突破了导波体外形的限制。采用平切导波体激波及正向跟踪流线生成乘波体,通过锥形流导波体和双锥导波体生成乘波体的算例分析了导波体与乘波体在纵向截面外形、横截面外形、前缘外形方面的对应关系。分析表明,双锥体在后体产生的激波或膨胀波能有效地改变导波体速度沿纵向分布特性,利用该特性可实现在俯视平面形状受约束的情况下对乘波体构型纵向稳定性的按需设计;通过改变导波体横截面外形,可以改变速度沿展向的
38、分布特性,从而实现乘波体横截面外形改变;乘波体前缘外形主要与导波体激波在俯视面上的外形相关,通过后体扩张的双锥外形能够生成双前缘乘波体外形。参 考 文 献: 1 N o n w e i l e r T R . A e r o d y n a m i c p r o b l e m o f m a n n e d s p a c ev e h i c l e s J . J o u r n a l o f t h e R o y a l A e r o n a u t i c a l S o c i e t y , 1 9 5 9 ,6 3 : 5 2 1 - 5 3 0 . 2 R a s m
39、 u s s e n M L . W a v e r i d e r c o n f i g u r a t i o n s d e r i v e d f r o mi n c l i n e d c i r c u l a r a n d e l l i p t i c c o n e s J . J o u r n a l o f S p a c e c r a f t a n dR o c k e t s , 1 9 8 0 , 1 7 ( 6 ) : 5 3 7 - 5 4 5 . 3 J o n e s J G , M o o r e K C , P i k e J , e t a
40、 l . A m e t h o d f o r d e s i g n i n gl i f t i n g c o n f i g u r a t i o n s f o r h i g h s u p e r s o n i c s p e e d s u s i n ga x i s y m m e t r i c f l o w f i e l d J . A r c h i v e o f A p p l i e d M e c h a n i c s ,1 9 6 8 , 3 7 ( 1 ) : 5 6 - 7 2 . 4 C e n t e r K , S o b i e c
41、z k y H , D o u g h e r t y F . I n t e r a c t i v e d e s i g n o fh y p e r s o n i c w a v e r i d e r g e o m e t r i e s R . A I A A 1 9 9 1 - 1 6 9 7 , 1 9 9 1 . 5 S o b i e c z k y H , D o u g h e r t y F C , J o n e s K . H y p e r s o n i c w a v e r i d e rd e s i g n f r o m g i v e n s
42、h o c k w a v e C / / T h e F i r s t I n t e r n a t i o n a lW a v e r i d e r S y m p o s i u m . M a r y l a n d : U n i v e r s i t y o f M a r y l a n d ,1 9 9 0 . 6 A n d e r s o n J D . F u n d a m e n t a l s o f a e r o d y n a m i c s M . 3 r de d i t i o n . M c G r a w - H i l l C o m p
43、 a n i e s , 2 0 0 1 . 7 M a r k J L , M a i r e L C . S h o c k - b a s e d w a v e r i d e r d e s i g n w i t hp r e s s u r e g r a d i e n t c o r r e c t i o n s a n d c o m p u t a t i o n a l s i m u l a t i o n s J . J o u r n a l o f A i r c r a f t , 2 0 0 5 , 4 2 ( 5 ) : 1 3 5 0 - 1 3 5
44、2 . 8 M a r c u A . L , K o j i r o S . E x p e r i m e n t a l i n v e s t i g a t i o n o f a M a c h3 . 5 w a v e r i d e r d e s i g n e d u s i n g c o m p u t a t i o n a l f l u i d d y n a m i c s J .A I A A J o u r n a l , 2 0 1 5 , 5 3 ( 6 ) : 1 5 9 0 - 1 6 0 1 . 9 L e G G , M a D W , L i
45、Z Y . C o m p u t a t i o n o f h y p e r s o n i cf l o w f i e l d s f o r e l l i p t i c - c o n e - d e r i v e d w a v e r i d e r J . J o u r n a l o fN a n j i n g U n i v e r s i t y o f S c i e n c e a n d T e c h n o l o g y ( N a t u r a lS c i e n c e ) , 2 0 0 6 , 3 0 ( 3 ) : 2 5 7 - 2
46、 6 0 . ( i n C h i n e s e )乐贵高,马大为,李自勇.椭圆锥乘波体高超声速流场数值计算 J .南京理工大学学报(自然科学版) , 2 0 0 6 , 3 0 ( 3 ) : 2 5 7 -2 6 0 . 1 0 M o r e t t i G . T h i r t y - s i x y e a r s o f s h o c k f i t t i n g m e t h o d J .C o m p u t e r s & F l u i d s , 2 0 0 2 , 3 1 : 7 1 9 - 7 2 3 . 1 1 L i u C Z , B a i P
47、, C h e n B Y , e t a l . R a p i d d e s i g n a n d o p t i m i z a t i o no f w a v e r i d e r f r o m 3 D f l o w J . J o u r n a l o f A s t r o n a u t i c s , 2 0 1 6 ,3 7 ( 5 ) : 5 3 5 - 5 4 3 . (i n C h i n e s e )刘传振,白鹏,陈冰雁,等.三维流场乘波体快速设计方法及多目标优化 J .宇航学报, 2 0 1 6 , 3 7 ( 5 ) : 5 3 5 - 5 4
48、3 . 1 2 , . M . . , 1 9 7 0 1 3 L y u Z J , W a n g J F , W u Y Z , e t a l . D e s i g n a n d a n a l y s i s o fm u l t i s t a g e c o m p r e s s i o n c o n e - d e r i v e d w a v e r i d e r c o n f i g u r a t i o n J . J o u r n a l o f A s t r o n a u t i c s , 2 0 1 5 , 3 6 ( 5 ) : 5 1 8 - 5 2 3 . ( i nC h i n e s e )吕侦军,王江峰,伍贻兆,等.多级压缩锥导乘波体设计与分析