CN105329462A - Changeable wall surface pressure distribution rule-based osculating flow field ride precursor design method - Google Patents

Changeable wall surface pressure distribution rule-based osculating flow field ride precursor design method Download PDF

Info

Publication number
CN105329462A
CN105329462A CN201510785607.XA CN201510785607A CN105329462A CN 105329462 A CN105329462 A CN 105329462A CN 201510785607 A CN201510785607 A CN 201510785607A CN 105329462 A CN105329462 A CN 105329462A
Authority
CN
China
Prior art keywords
point
waverider
flow field
curve
osculating
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201510785607.XA
Other languages
Chinese (zh)
Other versions
CN105329462B (en
Inventor
丁峰
柳军
刘珍
黄伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN201510785607.XA priority Critical patent/CN105329462B/en
Publication of CN105329462A publication Critical patent/CN105329462A/en
Application granted granted Critical
Publication of CN105329462B publication Critical patent/CN105329462B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64FGROUND OR AIRCRAFT-CARRIER-DECK INSTALLATIONS SPECIALLY ADAPTED FOR USE IN CONNECTION WITH AIRCRAFT; DESIGNING, MANUFACTURING, ASSEMBLING, CLEANING, MAINTAINING OR REPAIRING AIRCRAFT, NOT OTHERWISE PROVIDED FOR; HANDLING, TRANSPORTING, TESTING OR INSPECTING AIRCRAFT COMPONENTS, NOT OTHERWISE PROVIDED FOR
    • B64F5/00Designing, manufacturing, assembling, cleaning, maintaining or repairing aircraft, not otherwise provided for; Handling, transporting, testing or inspecting aircraft components, not otherwise provided for

Landscapes

  • Engineering & Computer Science (AREA)
  • Manufacturing & Machinery (AREA)
  • Transportation (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Structures Of Non-Positive Displacement Pumps (AREA)

Abstract

一种基于可变壁面压力分布规律的吻切流场乘波前体设计方法,步骤如下:S1,给定三种壁面压力分布规律不同曲线,将三曲线分别作为回转体母线,设计三种不同的回转体,分别求解绕这三种回转体的超声速轴对称流场,并将求解得到的流场作为三种壁面压力分布规律不同的基准流场;S2,根据乘波前体的性能要求,设计激波出口型线,应用吻切原理在每个吻切平面进行流场求解,所有的吻切平面流场组合成三维基准流场;S3,给定乘波前体前缘线在底部的投影曲线,根据激波出口型线,在S2的每个吻切平面内使用流线追踪法,生成吻切流场乘波前体气动构型。该方法设计的乘波前体既能够满足进气道入口高压升比的要求,又能够充分发挥乘波前体的高升阻比特性。

A design method for waveriding precursors of kissing flow field based on variable wall pressure distribution. rotators, respectively solve the supersonic axisymmetric flow fields around these three rotators, and use the obtained flow fields as the reference flow fields with different wall pressure distribution laws; S2, according to the performance requirements of the waverider precursor, Design the shock wave outlet profile, apply the kissing principle to solve the flow field at each kissing plane, and combine the flow fields of all the kissing planes into a three-dimensional reference flow field; The projected curve, according to the shock wave exit profile, uses the streamline tracing method in each kissing plane of S2 to generate the aerodynamic configuration of the waveriding precursor of the kissing flow field. The waverider precursor designed by this method can not only meet the requirement of high lift-to-lift ratio at the inlet of the inlet, but also fully utilize the high lift-to-drag ratio characteristic of the waverider precursor.

Description

基于可变壁面压力分布规律的吻切流场乘波前体设计方法Design method of waverider precursor for kiss-cut flow field based on variable wall pressure distribution

技术领域technical field

本发明涉及高超声速飞行器气动外形设计的技术领域,具体涉及一种基于可变壁面压力分布规律的吻切流场乘波前体设计方法。The invention relates to the technical field of aerodynamic shape design of a hypersonic vehicle, in particular to a design method for a waverider precursor of a kiss-cut flow field based on a variable wall surface pressure distribution law.

背景技术Background technique

吸气式高超声速飞行器是指飞行马赫数大于5、以吸气式发动机或其组合发动机为主要动力、能在大气层和跨大气层中远程飞行的飞行器,其应用形式包括高超声速巡航导弹、高超声速有人/无人飞机和空天飞机等多种飞行器。An air-breathing hypersonic vehicle refers to an aircraft with a flight Mach number greater than 5, powered by an air-breathing engine or a combined engine, and capable of long-distance flight in the atmosphere and across the atmosphere. Its application forms include hypersonic cruise missiles, hypersonic Manned/unmanned aircraft and aerospace aircraft and other aircraft.

高超声速飞行器的气动外形,主要有轴对称构型、升力体构型和乘波体构型三大类,其中乘波体构型利用激波压缩原理(乘波原理)实现了在高超声速飞行条件下高升阻比的气动要求,除了如何兼顾前缘气动热防护与气动性能的问题以外,对该构型的研究已趋于成熟。The aerodynamic shape of hypersonic vehicles mainly includes three types: axisymmetric configuration, lift body configuration and waverider configuration. Among them, the waverider configuration uses the shock wave compression principle (waverider principle) In addition to the aerodynamic requirements of high lift-to-drag ratio under high conditions, in addition to the problem of how to balance the aerodynamic thermal protection and aerodynamic performance of the leading edge, the research on this configuration has become mature.

乘波体构型设计方法的核心是基准流场的选择,不同的基准流场对应不同的设计方法,已有的方法有以下几种:楔导法、锥导法、楔锥法及基于超声速轴对称流动的吻切类生成方法。其中,吻切类设计方法主要包括吻切锥设计理论、吻切轴对称设计理论以及吻切流场设计理论。吻切锥理论使得乘波体底部横截面的激波线形状不再局限于圆弧或直线,而是可以根据进气道唇口外形进行合理设计。吻切轴对称理论提出吻切平面内的基准流场可以不再局限于锥形流场,而是可以根据设计需要选用合适的轴对称基准流场。应用吻切锥和吻切轴对称理论可以得到横向流动更均匀的乘波体构型,并能够有效地改善乘波前体(进气道入口)流场品质,有利于机身-进气道一体化设计。吻切流场理论提出每个吻切平面内的基准流场不再局限于同一个轴对称流场,而是可以根据设计需要在每个吻切平面内选用不同的轴对称基准流场。The core of the waverider configuration design method is the selection of the reference flow field. Different reference flow fields correspond to different design methods. The existing methods are as follows: wedge guidance method, cone guidance method, wedge cone method and supersonic-based Kiss-cut class generation method for axisymmetric flow. Among them, kiss-cut design methods mainly include kiss-cut cone design theory, kiss-cut axisymmetric design theory and kiss-cut flow field design theory. The kiss-cut cone theory makes the shock line shape of the bottom cross-section of the waverider no longer limited to circular arc or straight line, but can be reasonably designed according to the shape of the inlet lip. The kissing axisymmetric theory proposes that the reference flow field in the kissing plane can no longer be limited to the conical flow field, but an appropriate axisymmetric reference flow field can be selected according to the design requirements. The waverider configuration with more uniform lateral flow can be obtained by applying the kiss-cut cone and kiss-cut axisymmetric theory, and can effectively improve the flow field quality of the waverider precursor (intake inlet), which is beneficial to the fuselage-inlet integrated design. The kissing flow field theory proposes that the reference flow field in each kissing plane is no longer limited to the same axisymmetric flow field, but can choose different axisymmetric reference flow fields in each kissing plane according to design requirements.

2011年,贺旭照等人在文献《密切曲面锥乘波体——设计方法与性能分析》(力学学报)中提出了一种采用具有直线激波和等熵压缩波系的曲面锥作为基准流场的设计方法,有效克服了传统吻切锥乘波体压缩量不足及容积率偏小的缺点。但是,该方法存在一定的缺陷,即由于压升比的提高,导致乘波前体高升阻比特性得不到充分发挥。In 2011, He Xuzhao and others proposed a method using a curved cone with a straight line shock wave and an isentropic compression wave system as the reference flow field The design method effectively overcomes the shortcomings of insufficient compression and small volume ratio of traditional kiss-cut cone waveriders. However, this method has certain defects, that is, due to the increase of the pressure-to-lift ratio, the high lift-to-drag ratio of the waverider precursor cannot be fully utilized.

发明内容Contents of the invention

本发明提供一种基于可变壁面压力分布规律的吻切流场乘波前体设计方法,解决现有乘波前体只能按壁面压力恒定分布规律进行设计的不足。基于这一方法设计的乘波前体既能够满足进气道入口高压升比的要求,又能够充分发挥乘波前体的高升阻比特性。The invention provides a design method for waveriders in kiss-cut flow field based on the variable wall pressure distribution law, which solves the problem that the existing waverider body can only be designed according to the constant wall pressure distribution law. The waverider precursor designed based on this method can not only meet the requirement of high lift-to-lift ratio at the inlet of the intake port, but also give full play to the high lift-to-drag ratio characteristics of the waverider precursor.

为解决上述问题,本发明采用如下的技术方案:In order to solve the above problems, the present invention adopts the following technical solutions:

一种基于可变壁面压力分布规律的吻切流场乘波前体设计方法,包括以下步骤:A method for designing a waverider precursor in a kiss-cut flow field based on variable wall pressure distribution, comprising the following steps:

S1、给定三种壁面压力分布规律不同曲线,分别为壁面压力升高的曲线、壁面压力恒定的曲线及壁面压力降低的曲线,将三曲线分别作为回转体母线,设计三种不同的回转体,分别求解绕这三种回转体的超声速轴对称流场,并将求解得到的流场作为三种壁面压力分布规律不同的基准流场;S1. Given three different curves of wall pressure distribution law, which are the curve of rising wall pressure, the curve of constant wall pressure and the curve of decreasing wall pressure, respectively, the three curves are used as the busbar of the rotating body, and three different rotating bodies are designed. , respectively solve the supersonic axisymmetric flow fields around these three rotating bodies, and use the obtained flow fields as the reference flow fields with different wall pressure distribution laws;

S2、根据乘波前体的性能要求,将激波出口型线分为三段进行设计,并保证设计的型线在每段曲线的连接处曲率连续,即曲线的二阶导数连续,然后根据所设计的激波出口型线,应用吻切原理在每个吻切平面进行流场求解,所有的吻切平面流场组合成三维基准流场;S2. According to the performance requirements of the waverider precursor, divide the shock wave exit profile into three sections for design, and ensure that the designed profile has continuous curvature at the junction of each section of the curve, that is, the second derivative of the curve is continuous, and then according to For the designed shock wave outlet profile, the kissing principle is applied to solve the flow field in each kissing plane, and all the kissing plane flow fields are combined into a three-dimensional reference flow field;

S3、给定乘波前体前缘线在底部的投影曲线,根据步骤S2的激波出口型线,在步骤S2的每个吻切平面内使用流线追踪法,生成吻切流场乘波前体气动构型。S3. Given the projection curve of the leading edge line of the waverider front at the bottom, according to the shock wave exit profile in step S2, use the streamline tracing method in each kissing plane in step S2 to generate the kissing flow field waveriders Precursor aerodynamic configuration.

本发明采用吻切流场理论为基础,给定相同的来流条件,根据乘波体展向不同位置性能的要求,在进气道入口位置(即乘波体展向中间部位)选取壁面压力升高的基准流场,保证进气道入口高压升比特性;在乘波体展向两端选取壁面压力降低的基准流场,使得乘波前体高升阻比的特性得到充分发挥;在上述两个位置之间的部分选取壁面压力恒定的基准流场,作为过渡。The present invention adopts the kiss-cut flow field theory as the basis, and given the same incoming flow conditions, according to the performance requirements of different positions of the waverider, the wall pressure is selected at the inlet position of the air inlet (that is, the middle part of the waverider). The elevated reference flow field ensures the high-pressure lift ratio characteristics at the entrance of the inlet; the reference flow field with reduced wall pressure is selected at both ends of the waverider spanwise, so that the characteristics of the high lift-to-drag ratio of the waverider precursor can be fully utilized; in the above In the part between the two positions, the reference flow field with constant wall pressure is selected as the transition.

本发明的有益技术效果:Beneficial technical effect of the present invention:

本发明提供了一种新的乘波体设计方法,即,基于可变壁面压力分布规律的吻切流场乘波前体设计方法,这一方法克服了传统吻切锥乘波体只采取一种壁面压力分布规律的基准流场进行设计的不足,该方法可根据乘波体展向不同位置性能的要求,在进气道入口位置(即乘波体展向中间部位)选取壁面压力升高的基准流场,保证进气道入口高压升比特性;在乘波体展向两端选取壁面压力降低的基准流场,使得乘波前体高升阻比的特性得到充分发挥;The present invention provides a new waverider design method, that is, a waverider design method based on the variable wall pressure distribution law of the kissing flow field waverider. This method can be used to select the wall pressure increase at the inlet position of the inlet (that is, the middle part of the waverider span) according to the performance requirements of different positions of the waverider spanwise. The reference flow field ensures the high-pressure lift ratio characteristics of the inlet inlet; the reference flow field with reduced wall pressure is selected at both ends of the waverider spanwise, so that the characteristics of the high lift-to-drag ratio of the waverider precursor can be fully utilized;

附图说明Description of drawings

图1为三种不同壁面压力分布规律的基准流场剖面图;Fig. 1 is the cross-sectional view of the reference flow field of three different wall pressure distribution laws;

其中,直线1-3-11表示由母线1-2生成的圆锥在设计状态下产生的激波;2-11表示由2点发出的左行特征线;1-2-5、1-2-6、1-2-7分别表示壁面压力升高、壁面压力恒定、壁面压力降低的三条回转体母线;曲线3-4-8表示在母线1-2-5生成的回转体产生的流场中流线追踪获得的流线;曲线3-4-9表示在母线1-2-6生成的回转体产生的流场中流线追踪获得的流线;曲线3-4-10表示在母线1-2-7生成的回转体产生的流场中流线追踪获得的流线;12表示壁面压力升高回转体母线在点5处的倾斜角;13表示壁面压力恒定回转体母线在点6处的倾斜角;14表示壁面压力降低回转体母线在点7处的倾斜角;15表示基准流场的总长;16表示基准流场的总宽,即激波出口半径;17表示回转体的初始压缩角,即母线1-2-6和x轴的夹角;Among them, the straight line 1-3-11 represents the shock wave generated by the cone generated by the bus 1-2 in the design state; 2-11 represents the left-hand characteristic line emitted by 2 points; 1-2-5, 1-2- 6. 1-2-7 represent the three gyratory generatrixes with increasing wall pressure, constant wall pressure and decreasing wall pressure respectively; Curve 3-4-8 represents the flow field generated by the gyrator generated by busbar 1-2-5 The streamline obtained by streamline tracing; curve 3-4-9 represents the streamline obtained by streamline tracing in the flow field generated by the rotator generated by bus 1-2-6; curve 3-4-10 represents the streamline obtained by the bus 1-2-6 2-7 The streamlines obtained by tracing the streamlines in the flow field generated by the gyrator; 12 represents the inclination angle of the generatrix of the gyrator at point 5 with increasing wall pressure; 13 represents the inclination angle of the generatrix of the gyrator with constant wall pressure at point 6 Inclination angle; 14 indicates the inclination angle of the generatrix of the rotary body at point 7 when the wall pressure decreases; 15 indicates the total length of the reference flow field; 16 indicates the total width of the reference flow field, that is, the radius of the shock wave exit; 17 indicates the initial compression angle of the rotary body , that is, the angle between the bus 1-2-6 and the x-axis;

图2为激波出口截面示意图;Fig. 2 is a schematic diagram of a section of a shock wave outlet;

其中,18为乘波前体前缘线在激波出口截面的投影线,也为乘波前体上表面后缘线;19为乘波前体下表面后缘线;20为激波出口型线;21为点25对应的吻切锥;22为点25对应的吻切锥顶点在激波出口截面的投影点;23点为点22和点25连线和乘波前体上表面后缘线的交点;点24为点22和点25连线和乘波前体下表面后缘线的交点;25为激波出口型线上的任意点;26为过点25的吻切平面;曲线27-28、28-29、29-29'、28'-29'和27'-28'构成激波出口型线20;30为过点25的曲率圆;Among them, 18 is the projection line of the front edge line of the waverider body on the section of the shock wave exit, which is also the trailing edge line of the upper surface of the waverider body; 19 is the trailing edge line of the lower surface of the waverider body; 20 is the shock wave exit type line; 21 is the kissing cone corresponding to point 25; 22 is the projection point of the apex of the kissing cone corresponding to point 25 on the shock wave exit section; point 23 is the line connecting point 22 and point 25 and the trailing edge of the upper surface of the waverider precursor The intersection point of the line; point 24 is the intersection point of the line connecting point 22 and point 25 and the trailing edge line of the lower surface of the waverider precursor; 25 is any point on the shock wave exit profile line; 26 is the kissing plane passing through point 25; the curve 27-28, 28-29, 29-29', 28'-29' and 27'-28' constitute the shock wave exit profile 20; 30 is the curvature circle passing through point 25;

图3为吻切面内流线追踪示意图;Figure 3 is a schematic diagram of streamline tracing in the kissing plane;

其中,31为点25对应的吻切锥顶点;32为点25对应的前缘点;Among them, 31 is the apex of the kissing cone corresponding to point 25; 32 is the leading edge point corresponding to point 25;

具体实施方式detailed description

本发明提供一种可变壁面压力分布规律的吻切流场乘波前体设计方法,包括以下步骤:The invention provides a method for designing a waverider precursor of a kiss-cut flow field with a variable wall pressure distribution law, which includes the following steps:

步骤S1、给定三种壁面压力分布规律不同曲线,分别为壁面压力升高的曲线、壁面压力恒定的曲线及壁面压力降低的曲线,将三曲线分别作为回转体母线,设计三种不同的回转体,分别求解绕这三种回转体的超声速轴对称流场,并将求解得到的流场作为三种壁面压力分布规律不同的基准流场。Step S1, given three different curves of wall pressure distribution law, which are the curve of wall pressure increase, the curve of wall pressure constant and the curve of wall pressure decrease, respectively, the three curves are used as the generatrix of the rotary body, and three different rotary curves are designed. The supersonic axisymmetric flow fields around these three rotating bodies are respectively solved, and the obtained flow fields are used as the reference flow fields with different wall pressure distribution laws.

(1)给定直线段1-2-6,直线段的一端点为点1,直线段的另一端点为点6,点2为直线段中间的一点,点1位于x轴上,直线段1-2-6与x轴的夹角为17,夹角17等于直线段在点6处的夹角13,以直线段1-2-6作为回转体母线,使其绕着x轴360°旋转形成直圆锥,所成直圆锥即为壁面压力恒定的回转体;(1) Given a straight line segment 1-2-6, one end point of the straight line segment is point 1, the other end point of the straight line segment is point 6, point 2 is a point in the middle of the straight line segment, point 1 is located on the x-axis, and the straight line segment The included angle between 1-2-6 and the x-axis is 17, and the included angle 17 is equal to the included angle 13 of the straight line segment at point 6, and the straight line segment 1-2-6 is used as the generatrix of the revolving body, so that it circles the x-axis 360° Rotate to form a straight cone, which is a body of revolution with constant wall pressure;

(2)在点6的正上方给定一点5,在点2和点5之间设计一外扩张曲线2-5,本发明外扩张曲线2-5为抛物线形式的曲线,曲线2-5和直线1-2在点2处一阶导数连续,并结合2点和5点的坐标,利用公式(1)能够唯一确定曲线2-5,曲线2-5在点5处的切线与x轴夹角为12,夹角12大于夹角17,由曲线1-2-5绕着x轴360°旋转形成外扩张式曲面锥,所成外扩张式曲面锥即为壁面压力升高的回转体;(2) a point 5 is given directly above the point 6, and an outer expansion curve 2-5 is designed between the point 2 and the point 5, and the outer expansion curve 2-5 of the present invention is a curve of a parabolic form, and the curve 2-5 and The first derivative of straight line 1-2 at point 2 is continuous, and combined with the coordinates of point 2 and point 5, the formula (1) can be used to uniquely determine curve 2-5, the tangent of curve 2-5 at point 5 and the x-axis The angle is 12, and the included angle 12 is greater than the included angle 17. The curve 1-2-5 rotates 360° around the x-axis to form an outwardly expanding curved surface cone, and the formed outwardly expanding curved surface cone is a body of revolution with increased wall pressure;

y=ax2+bx+c(1)y=ax 2 +bx+c(1)

(3)在点6正下方给定一点7,在点2和点7之间设计一内膨胀曲线2-7,本发明内膨胀曲线2-7为抛物线形式的曲线,曲线2-7和直线1-2在点2处一阶导数连续,并结合2点和7点的坐标,利用公式(1)能够唯一确定曲线2-7,曲线2-7在7点处的切线与x轴夹角为14,夹角14小于夹角17;由曲线1-2-7绕着x轴360°旋转形成内膨胀式曲面锥,所成内膨胀式曲面锥即为壁面压力降低的回转体;(3) a point 7 is given directly below point 6, and an internal expansion curve 2-7 is designed between point 2 and point 7. The internal expansion curve 2-7 of the present invention is a curve of parabolic form, and curve 2-7 and straight line 1-2 The first-order derivative at point 2 is continuous, combined with the coordinates of 2 points and 7 points, using formula (1) can uniquely determine curve 2-7, the angle between the tangent line of curve 2-7 at point 7 and the x-axis is 14, and the included angle 14 is smaller than the included angle 17; the curve 1-2-7 rotates around the x-axis 360° to form an internal expansion curved surface cone, and the formed internal expansion curved surface cone is a body of revolution with reduced wall pressure;

(4)给定来流条件(马赫数、静压和静温),采用有旋特征线方法,求解绕这三种回转体的超声速轴对称流场,得到三种基准流场,分别定义为:壁面压力恒定的基准流场、壁面压力升高的基准流场和壁面压力降低的基准流场。有旋特征线方法为本领域的公知技术,具体可参见“《气体动力学》,M.J.左克罗,J.D.霍夫曼,国防工业出版社,1984年,p138-195”。(4) Given the incoming flow conditions (Mach number, static pressure, and static temperature), the method of rotating characteristic lines is used to solve the supersonic axisymmetric flow field around these three rotating bodies, and three reference flow fields are obtained, which are defined as : The reference flow field with constant wall pressure, the reference flow field with increased wall pressure and the reference flow field with reduced wall pressure. The spin characteristic line method is a well-known technology in the art. For details, please refer to ""Gas Dynamics", M.J. Zocrow, J.D. Hoffman, National Defense Industry Press, 1984, p138-195".

其中,在壁面压力升高的基准流场中进行流线追踪得到的乘波体构型具有高压升比,在壁面压力升高的流场中追踪的流线对应图1中曲线3-4-8;在壁面压力降低的基准流场中进行流线追踪得到的乘波体构型具有高升阻比,在壁面压力降低的流场中追踪的流线对应图1中的曲线3-4-10。Among them, the waverider configuration obtained by streamline tracing in the reference flow field with increased wall pressure has a high lift ratio, and the streamlines traced in the flow field with increased wall pressure correspond to curve 3-4- in Figure 1 8. The waverider configuration obtained by streamline tracing in the reference flow field with reduced wall pressure has a high lift-to-drag ratio, and the streamline traced in the flow field with reduced wall pressure corresponds to curve 3-4-10 in Figure 1 .

(注:压升比=乘波前体底部横截面压力/来流压力)(Note: Pressure rise ratio = cross-sectional pressure at the bottom of the waverider precursor / incoming flow pressure)

步骤S2、应用吻切原理将步骤S1中得到的三个轴对称基准流场组装成三维基准流场;Step S2, assembling the three axisymmetric reference flow fields obtained in step S1 into a three-dimensional reference flow field by applying the kiss cut principle;

首先,设计激波出口型线。First, design the shock wave exit profile.

激波出口型线以其中线呈左右对称,将激波出口型线分为五段,分别为乘波前体左段27-28、乘波前体左过渡段28-29、进气道入口位置段29-29'、乘波前体右过渡段28'-29'、乘波前体右段27'-28',其中乘波前体左段27-28和乘波前体右段27'-28'是呈左右对称的,乘波前体左过渡段28-29和乘波前体右过渡段28'-29'是呈左右对称的。将乘波前体左段27-28和乘波前体右段27'-28'统称为乘波前体两端段,将乘波前体左过渡段28-29和乘波前体右过渡段28'-29'统称为乘波前体过渡段,因此设计激波出口型线时只需设计乘波前体两端段27-28和27'-28'、乘波前体过渡段28-29和28'-29'、进气道入口位置段29-29'即可;The shock wave exit profile is symmetrical about the center line, and the shock wave exit profile is divided into five sections, namely, the left section 27-28 of the waverider precursor, the left transition section 28-29 of the waverider precursor, and the inlet of the inlet The position section 29-29', the right transition section 28'-29' of the waverider body, the right section 27'-28' of the waverider body, of which the left section 27-28 of the waverider body and the right section 27 of the waverider body '-28' is left-right symmetrical, and the left transition section 28-29 of the waveriding precursor and the right transition section 28'-29' of the waveriding precursor are left-right symmetrical. The left section 27-28 of the waveriding front body and the right section 27'-28' of the waveriding front body are collectively referred to as both end sections of the waveriding front body, and the left transition section 28-29 of the waveriding front body and the right transition section of the waveriding front body Sections 28'-29' are collectively referred to as the transition section of the waverider precursor, so when designing the shock wave exit profile, it is only necessary to design the two end sections 27-28 and 27'-28' of the waverider precursor, and the transition section 28 of the waverider precursor. -29 and 28'-29', the inlet position section of the air inlet 29-29' is sufficient;

如图2所示,将激波出口型线分为三段进行设计,这三段分别是乘波前体两端段27-28和27'-28'、乘波前体过渡段28-29和28'-29'、进气道入口位置段29-29'。给定乘波前体左段、乘波前体左过渡段、进气道入口位置段、乘波前体右过渡段和乘波前体右段之间连接点的坐标,即给定点27、点28、点29、点27'、点28'和点29'的坐标。将进气道入口位置段29-29'设计成水平直线,波前体两端段27-28和27'-28'与乘波前体过渡段28-29和28'-29'均设计为四次曲线,各段曲线连接点处左右两侧二阶导数相同以保证段与段之间曲率连续,即设计的曲线在点28、29、28'、29'处二阶导数连续。利用四次曲线方程即式(2),由连接点的坐标和连接点处二阶导数连续可以唯一确定四次曲线方程中的系数,进而确定波前体两端段和乘波前体过渡段对应的四次曲线。As shown in Figure 2, the shock wave exit profile is divided into three sections for design, these three sections are the two end sections 27-28 and 27'-28' of the waverider precursor, and the transition section 28-29 of the waverider precursor. And 28'-29', air inlet inlet position section 29-29'. Given the coordinates of the connection points between the left section of the waverider front body, the left transition section of the waverider front body, the inlet position section, the right transition section of the waverider body and the right section of the waverider front body, that is, the given points 27, The coordinates of point 28, point 29, point 27', point 28' and point 29'. The inlet position section 29-29' of the air inlet is designed as a horizontal straight line, the two end sections 27-28 and 27'-28' of the wave front body and the transition sections 28-29 and 28'-29' of the waverider body are designed as For quartic curves, the second order derivatives at the left and right sides of the connecting points of each section of the curve are the same to ensure that the curvature between sections is continuous, that is, the designed curve has continuous second order derivatives at points 28, 29, 28', and 29'. Using the quartic curve equation (2), the coefficients in the quartic curve equation can be uniquely determined from the coordinates of the connection point and the continuity of the second order derivative at the connection point, and then the two ends of the wave front body and the transition section of the wave-riding body can be determined The corresponding quartic curve.

y=ax4+bx2+c(2)y=ax 4 +bx 2 +c(2)

其中,设计乘波体两端段27-28和27'-28'的基准流场为S1中定义的壁面压力降低的基准流场,用来提高乘波前体升阻比;设计进气道入口位置段29-29'采用S1中定义的壁面压力分布升高的基准流场,以保证进气道入口高压升比的要求;乘波前体过渡段28-29和28'-29'采用S1中定义的壁面压力恒定的基准流场;Among them, the design reference flow field of the two end sections 27-28 and 27'-28' of the waverider is the reference flow field of the wall pressure reduction defined in S1, which is used to increase the lift-to-drag ratio of the waverider precursor; The inlet position section 29-29' adopts the reference flow field with increased wall pressure distribution defined in S1 to ensure the high pressure lift ratio requirement at the entrance of the inlet; the waverider precursor transition sections 28-29 and 28'-29' adopt The reference flow field with constant wall pressure defined in S1;

如图2所示,对激波出口型线20进行离散,在激波出口型线上每5mm取一个点,可以保证不同点产生的流线能够形成光滑曲面;在激波出口型线离散后所取得的离散点中任意取一点25,得到过该点的曲率圆30,曲率圆30即为点25对应的吻切曲面锥在出口截面的吻切圆,该吻切曲面锥的轴线与X轴平行,过点25和吻切曲面锥轴线的平面为吻切平面26;点22为曲率圆30的圆心;As shown in Figure 2, the shock wave exit profile 20 is discretized, and a point is taken every 5 mm on the shock wave exit profile, which can ensure that the streamlines generated at different points can form a smooth surface; after the shock wave exit profile is discretized Randomly take a point 25 among the discrete points obtained, and obtain the curvature circle 30 passing through this point, and the curvature circle 30 is the kissing circle of the kissing surface cone corresponding to point 25 in the exit section, the axis of the kissing curved surface cone is in line with X The axes are parallel, and the plane passing through the point 25 and the cone axis of the kissing curved surface is the kissing plane 26; the point 22 is the center of the curvature circle 30;

如果点25位于乘波前体过渡段,则选择步骤S1中壁面压力恒定的基准流场作为过点25的吻切面内流场;如果点25位于乘波体两端段,则选择步骤S1中壁面压力降低的基准流场作为过点25的吻切面内流场;如果点25位于进气道入口位置段,则选择步骤S1中壁面压力分布升高的基准流场作为过点25的吻切面内流场;以此规律对激波出口型线离散后所取的每一个点根据起各自所处的激波段来选择对应壁面压力分布的基准流场作为该点对应的吻切流场。图2所示的点25位于乘波前体左过渡段28-29,因此,选择步骤S1中壁面压力恒定的基准流场作为过点25的吻切面内流场。If point 25 is located at the transition section of the waverider body, select the reference flow field with constant wall pressure in step S1 as the flow field in the kissing plane passing through point 25; if point 25 is located at both ends of the waverider, select step S1 The reference flow field with reduced wall pressure is taken as the flow field in the kissing plane passing through point 25; if point 25 is located at the entrance of the inlet, the reference flow field with increased wall pressure distribution in step S1 is selected as the kissing plane passing through point 25 Inner flow field; each point taken after discretizing the shock wave outlet profile according to this rule selects the reference flow field corresponding to the wall pressure distribution as the kiss cut flow field corresponding to the point according to the respective shock wave segment. The point 25 shown in FIG. 2 is located at the left transition section 28-29 of the waverider body. Therefore, the reference flow field with constant wall pressure in step S1 is selected as the flow field in the kissing plane passing through the point 25.

步骤S3、给定乘波前体前缘线在底部的投影曲线,根据步骤S2的激波出口型线,在步骤S2的每个吻切平面内使用流线追踪法,生成吻切流场乘波前体气动构型。Step S3, given the projection curve of the leading edge line of the waverider precursor at the bottom, according to the shock wave outlet profile in step S2, use the streamline tracing method in each kissing plane in step S2 to generate the kissing flow field product Aerodynamic configuration of wave precursors.

(1)首先给定乘波前体前缘线在底部的投影曲线,即乘波体上表面后缘线。从步骤S2设计的激波出口型线上等间距地取出足够密的离散点,一般每5mm取一个点,可以保证不同点产生的流线能够形成光滑曲面,具体可参见丁峰.高超声速滑翔-巡航两级乘波设计方法研究[D].长沙:国防科学与技术大学(硕士).2012。(1) Firstly, the projection curve of the leading edge line of the waverider body at the bottom is given, that is, the trailing edge line of the upper surface of the waverider body. Take sufficiently dense discrete points at equal intervals from the shock wave exit profile designed in step S2. Generally, one point is taken every 5 mm to ensure that the streamlines generated by different points can form a smooth surface. For details, see Ding Feng. Hypersonic Gliding - Research on two-stage waveriding design method for cruising [D]. Changsha: University of National Defense Science and Technology (Master). 2012.

(2)如图2和图3所示,在激波出口型线离散后所取得的离散点中任意取一点25,由点25获得过点25的曲率圆30,过25点的曲率圆即为25点对应的吻切锥激波在出口截面的吻切圆,吻切锥的轴线平行于x轴(即自由流线)。点22为过点25的曲率圆30的圆心,也为点25对应的吻切锥顶点31(图3)在激波出口截面的投影点。25点、22点和31点构成过点25的吻切面26。点25和点22的连线交乘波体上表面后缘线于点23,由点23做平行于x轴的直线交吻切锥激波于点32,点32也为过点25的吻切面内对应的前缘点。从前缘点32向后进行流线追踪至激波出口截面,得到后缘点24,前缘点与后缘点连成的曲线32-24即为点25对应的下表面流线,将点32和点23用直线连接,获得点25对应的上表面流线;对激波出口型线离散后所取得的所有离散点采用同样的方式获得各自吻切面内的前缘点、乘波体下表面流线、乘波体上表面流线和下表面后缘线上的点。(2) As shown in Figure 2 and Figure 3, a point 25 is arbitrarily selected among the discrete points obtained after the shock wave exit profile is discrete, and the curvature circle 30 passing through point 25 is obtained from point 25, and the curvature circle passing through point 25 is is the kissing circle of the kissing cone shock corresponding to point 25 at the exit section, and the axis of the kissing cone is parallel to the x-axis (ie, the free streamline). The point 22 is the center of the curvature circle 30 passing through the point 25, and is also the projection point of the vertex 31 ( FIG. 3 ) of the kissing cone corresponding to the point 25 on the shock wave exit section. 25 o'clock, 22 o'clock and 31 o'clock constitute the kiss section 26 passing through point 25. The line connecting point 25 and point 22 intersects the rear edge line of the upper surface of the waverider at point 23, and a straight line parallel to the x-axis is drawn from point 23 to meet the tangent cone shock at point 32, and point 32 is also the kiss passing through point 25 Corresponding front edge point in the tangential plane. Tracing the streamline from the leading edge point 32 backward to the shock wave exit section, the trailing edge point 24 is obtained, and the curve 32-24 formed by connecting the leading edge point and the trailing edge point is the lower surface streamline corresponding to the point 25. Point 32 Connect with point 23 with a straight line to obtain the upper surface streamline corresponding to point 25; for all the discrete points obtained after the shock wave exit profile is discretized, use the same method to obtain the leading edge points in the respective kissing planes and the lower surface of the waverider Streamlines, streamlines on the upper surface of the waverider, and points on the trailing edge line on the lower surface.

(3)一系列前缘点平滑连接构成乘波体的前缘线;一系列上表面流线构成乘波体上表面;一系列下表面流线构成乘波体下表面;一系列后缘线上的点平滑连接构成乘波体下表面后缘线;由前缘线、后缘线及乘波体下表面构成乘波前体气动构型。(3) A series of leading edge points are smoothly connected to form the leading edge line of the waverider; a series of upper surface streamlines form the upper surface of the waverider; a series of lower surface streamlines form the lower surface of the waverider; a series of trailing edge lines The points above are connected smoothly to form the trailing edge line of the lower surface of the waverider; the aerodynamic configuration of the waverider front is formed by the leading edge line, the trailing edge line and the lower surface of the waverider.

Claims (5)

1., based on an osculating flow field waverider forebody derived method of designing for the variable wall surface pressure regularity of distribution, it is characterized in that, comprise the following steps:
S1, the different curve of given three kinds of wall pressure distribution rules, be respectively the curve that the constant curve of curve, wall pressure that wall pressure raises and wall pressure reduce, using three curves as gyro-rotor bus, design three kinds of different gyro-rotors, solve the supersonic speed axisymmetric flow field around these three kinds of gyro-rotors respectively, and will the flow field that obtains be solved as the different benchmark flow field of three kinds of wall pressure distribution rules;
S2, performance requriements according to waverider forebody derived, shock wave being exported molded line is divided into three sections to design, and ensure the junction continual curvature of molded line at every section of curve of design, namely the second derivative of curve is continuous, then according to designed shock wave outlet molded line, application osculating principle carries out flow field calculation at each osculating plane, all osculating plane group of flow fields synthesis three-dimensional references flow fields;
S3, the given drop shadow curve of waverider forebody derived costa in bottom, according to the shock wave outlet molded line of step S2, use streamlined impeller method, generate osculating flow field waverider forebody derived aerodynamic configuration in each osculating plane of step S2.
2. the osculating flow field waverider forebody derived method of designing based on the variable wall surface pressure regularity of distribution according to claim 1, it is characterized in that, the concrete grammar of step S1 is:
The given linear portion 1-2-6 of S1.1, one end points of linear portion is point 1, another end points of linear portion is point 6, point 2 is a bit in the middle of linear portion, and point 1 is positioned in x-axis, and the angle of linear portion 1-2-6 and x-axis is 17, angle 17 equals the angle 13 of linear portion at point 6 place, using linear portion 1-2-6 as gyro-rotor bus, make its around x-axis 360 ° rotate formed straight circular cone, become straight circular cone to be the constant gyro-rotor of wall pressure;
S1.2 is directly over point 6 given 1: 5, design one outer dilation curve 2-5 between point 2 and point 5, curve 2-5 and straight line 1-2 is continuous in point 2 place first derivative, and in conjunction with the coordinate of 2 and 5, can uniquely determine curve 2-5, the tangent line of curve 2-5 at point 5 place and x-axis angle are 12, and angle 12 is greater than angle 17, rotated to be formed around x-axis 360 ° by curve 1-2-5 and extend out a formula curved surface and bore, institute becomes outer expanding curved surface cone to be the gyro-rotor of wall pressure rising;
S1.3 is immediately below point 6 given 1: 7, expansion curve 2-7 in design one between point 2 and point 7, curve 2-7 and straight line 1-2 is continuous in point 2 place first derivative, and in conjunction with the coordinate of 2 and 7, can uniquely determine curve 2-7, curve 2-7 is 14 at 7 tangent lines located and x-axis angle, and angle 14 is less than angle 17; By curve 1-2-7 around x-axis 360 ° rotate is formed inner intumescence type curved surface bore, institute become inner intumescence type curved surface cone be wall pressure reduction gyro-rotor;
The given inlet flow conditions of S1.4, employing has revolves characteristic line method, solve the supersonic speed axisymmetric flow field around these three kinds of gyro-rotors, obtain three kinds of benchmark flow fields, be defined as respectively: the benchmark flow field that the benchmark flow field that wall pressure is constant, wall pressure raise and the benchmark flow field that wall pressure reduces.
3. the osculating flow field waverider forebody derived method of designing based on the variable wall surface pressure regularity of distribution according to claim 2, is characterized in that, outer dilation curve 2-5 and interior expansion curve 2-7 is the curve of parabolic.
4. the osculating flow field waverider forebody derived method of designing based on the variable wall surface pressure regularity of distribution according to claim 3, it is characterized in that, the concrete grammar of described S2 is:
S2.1 designs shock wave outlet molded line;
Shock wave outlet molded line is symmetrical with its center line, shock wave is exported molded line and be divided into five sections, be respectively left section of waverider forebody derived, the left transition phase of waverider forebody derived, inlet mouth position section, the right transition phase of waverider forebody derived, right section of waverider forebody derived, wherein left section of waverider forebody derived and right section of waverider forebody derived are symmetrical, the left transition phase of waverider forebody derived and the right transition phase of waverider forebody derived are symmetrical, left for waverider forebody derived section and right section of waverider forebody derived are referred to as waverider forebody derived two ends section, left for waverider forebody derived transition phase and the right transition phase of waverider forebody derived are referred to as waverider forebody derived transition phase;
Left section of given waverider forebody derived, the left transition phase of waverider forebody derived, inlet mouth position section, the coordinate of point of connection between the right transition phase of waverider forebody derived and right section of waverider forebody derived, inlet mouth position section is designed to horizontal linear, wavefront body two ends section and waverider forebody derived transition phase are designed to quartic curve, each section of curve point of connection place left and right sides second derivative is identical to ensure continual curvature between section and section, and the point of connection namely between waverider forebody derived two ends section and waverider forebody derived transition phase, point of connection place second derivative between waverider forebody derived transition phase and inlet mouth position section are continuous;
Utilize quartic curve equation, uniquely can be determined the coefficient in quartic curve equation by the coordinate of point of connection and point of connection second derivative continuously, and then determine wavefront body two ends section and quartic curve corresponding to waverider forebody derived transition phase; Wherein quartic curve equation is as follows:
y=ax 4+bx 2+c
Wherein, the benchmark flow field designing Waverider two ends section is the benchmark flow field that the wall pressure defined in S1 reduces; The benchmark flow field that design inlet mouth position section adopts the wall pressure distribution defined in S1 to raise; The benchmark flow field that the wall pressure defined in waverider forebody derived transition phase employing S1 is constant;
S2.2 obtains three-dimensional references flow field;
Carry out discrete to shock wave outlet molded line, on shock wave outlet molded line, every 5mm gets a point, can ensure that the streamline that difference produces can form smooth surface; 1: 25 is got arbitrarily in the discrete rear acquired discrete point of shock wave outlet molded line, the circle of curvature 30 of this point must be, the circle of curvature 30 is the osculating circle of osculating curved surface cone in outlet of a little 25 correspondences, the axis of this osculating curved surface cone is parallel with X-axis, cross a little 25 and the plane of osculating curved surface axis of cone line be osculating plane 26, point 22 is the center of circle of the circle of curvature 30;
As fruit dot 25 is positioned at waverider forebody derived transition phase, then the benchmark flow field that in selection step S1, wall pressure is constant is as the osculating face flow field crossing point 25; As fruit dot 25 is positioned at Waverider two ends section, then select the benchmark flow field of wall pressure reduction in step S1 as the osculating face flow field crossing point 25; As fruit dot 25 is positioned at inlet mouth position section, then select the benchmark flow field of wall pressure distribution rising in step S1 as the osculating face flow field crossing point 25; Select the benchmark flow field of corresponding wall pressure distribution as the osculating flow field of this some correspondence to discrete rear each the got point of shock wave outlet molded line according to shock wave section residing separately using this rule.
5. the osculating flow field waverider forebody derived method of designing based on the variable wall surface pressure regularity of distribution according to claim 3, it is characterized in that, the concrete grammar of described S3 is:
The first given drop shadow curve of waverider forebody derived costa in bottom of S3.1, i.e. Waverider upper surface trailing edge line; Take out enough close discrete point with exporting the first-class spacing of molded line from the shock wave of step S2 design, on shock wave outlet molded line, every 5mm gets a point, to ensure that the streamline that difference produces can form smooth surface;
S3.2 gets arbitrarily 1: 25 in the discrete rear acquired discrete point of shock wave outlet molded line, the circle of curvature 30 of a little 25 was obtained by point 25, the circle of curvature crossing at 25 is the osculating circle of osculating cone shock wave in outlet of a little 25 correspondences, the axis being parallel of osculating cone is in x-axis, point 22 was the center of circle of the circle of curvature 30 of a little 25, and point 22 is also the subpoint of osculating conic node 31 in shock wave outlet of point 25 correspondences; Point 25, point 22 and point 31 formed the osculating face 26 of a little 25; The line of point 25 and point 22 hands over Waverider upper surface trailing edge line in point 23, and the straight line being made to be parallel to x-axis by point 23 hands over osculating to bore shock wave in point 32, point 32 be also a little 25 osculating face in the leading edge point of correspondence; Carry out streamlined impeller backward to shock wave outlet from leading edge point 32, obtain trailing edge point 24, leading edge point and the curve 32-24 that trailing edge point is linked to be are a little 25 corresponding lower surface streamlines, will put 32 and be connected with straight line with point 23, and obtain the upper surface streamline putting 25 correspondences; The discrete rear acquired all discrete points of shock wave outlet molded line are adopted to the point on the leading edge point obtained in the same way in respective osculating face, Waverider lower surface streamline, Waverider upper surface streamline and lower surface trailing edge line;
The a series of leading edge point smooth connection of S3.3 forms the costa of Waverider; A series of upper surface streamline forms Waverider upper surface; A series of lower surface streamline forms Waverider lower surface; Point smooth connection on a series of trailing edge line forms Waverider lower surface trailing edge line; Waverider forebody derived aerodynamic configuration is formed by costa, trailing edge line and Waverider lower surface.
CN201510785607.XA 2015-11-16 2015-11-16 Osculating flow field based on variable wall surface pressure regularity of distribution waverider forebody derived method for designing Active CN105329462B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510785607.XA CN105329462B (en) 2015-11-16 2015-11-16 Osculating flow field based on variable wall surface pressure regularity of distribution waverider forebody derived method for designing

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510785607.XA CN105329462B (en) 2015-11-16 2015-11-16 Osculating flow field based on variable wall surface pressure regularity of distribution waverider forebody derived method for designing

Publications (2)

Publication Number Publication Date
CN105329462A true CN105329462A (en) 2016-02-17
CN105329462B CN105329462B (en) 2016-08-24

Family

ID=55280293

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510785607.XA Active CN105329462B (en) 2015-11-16 2015-11-16 Osculating flow field based on variable wall surface pressure regularity of distribution waverider forebody derived method for designing

Country Status (1)

Country Link
CN (1) CN105329462B (en)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106428620A (en) * 2016-10-31 2017-02-22 中国人民解放军国防科学技术大学 Design method of ridge-shaped osculating-cone wave-rider with large capacity and high lift-to-drag ratio
CN107016199A (en) * 2017-04-13 2017-08-04 中国人民解放军国防科学技术大学 It is a kind of that the design method for moving bulge is arranged without shock-boundary
CN108019279A (en) * 2017-12-07 2018-05-11 中国人民解放军国防科技大学 A Design Method for Hypersonic Inlet
CN109279043A (en) * 2018-10-23 2019-01-29 中国人民解放军国防科技大学 Design method of Von Karman waverider with low-speed airfoil
CN109598062A (en) * 2018-12-04 2019-04-09 中国人民解放军国防科技大学 Design method of osculating flow field waverider with variable wall surface pressure distribution rule
CN110414016A (en) * 2018-04-27 2019-11-05 中国航天科工飞航技术研究院(中国航天海鹰机电技术研究院) The Waverider geometry parameterization design method and system of ultrahigh speed pipeline transportation tool
CN112765725A (en) * 2020-12-30 2021-05-07 四川京航天程科技发展有限公司 Analytic Riemann resolving method for multi-dimensional Euler equation
CN112949199A (en) * 2021-03-15 2021-06-11 中国科学院力学研究所 Method and system for optimizing longitudinal stability of power waverider
CN112948976A (en) * 2021-03-15 2021-06-11 中国科学院力学研究所 Plane front edge lifting body front edge line optimization method and system and lifting body
CN114313253A (en) * 2022-03-03 2022-04-12 中国空气动力研究与发展中心计算空气动力研究所 Aerodynamic layout and design method of high lift-drag ratio air-breathing hypersonic aircraft
CN114852366A (en) * 2022-05-30 2022-08-05 厦门大学 A Design Method for Dual Waverider Configurations with Simultaneous Controllable Shock Shape and Pressure Distribution

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103770935A (en) * 2013-12-13 2014-05-07 中国航天空气动力技术研究院 Wave rider appearance designing method
CN104192302A (en) * 2014-07-18 2014-12-10 中国人民解放军国防科学技术大学 Waverider designing method based on reference flow field of revolution body of cuspidal Von Karman curve
CN104724281A (en) * 2015-02-13 2015-06-24 中国科学院力学研究所 Combined front-edge wave rider design method and combined front-edge wave rider
CN104973266A (en) * 2015-07-16 2015-10-14 中国人民解放军国防科学技术大学 Gliding-cruising two-stage wave rider design method based on osculating cone principle

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103770935A (en) * 2013-12-13 2014-05-07 中国航天空气动力技术研究院 Wave rider appearance designing method
CN104192302A (en) * 2014-07-18 2014-12-10 中国人民解放军国防科学技术大学 Waverider designing method based on reference flow field of revolution body of cuspidal Von Karman curve
CN104724281A (en) * 2015-02-13 2015-06-24 中国科学院力学研究所 Combined front-edge wave rider design method and combined front-edge wave rider
CN104973266A (en) * 2015-07-16 2015-10-14 中国人民解放军国防科学技术大学 Gliding-cruising two-stage wave rider design method based on osculating cone principle

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106428620B (en) * 2016-10-31 2017-07-21 中国人民解放军国防科学技术大学 A kind of big volume high lift-drag ratio ridged osculating cone Waverider design method
CN106428620A (en) * 2016-10-31 2017-02-22 中国人民解放军国防科学技术大学 Design method of ridge-shaped osculating-cone wave-rider with large capacity and high lift-to-drag ratio
CN107016199A (en) * 2017-04-13 2017-08-04 中国人民解放军国防科学技术大学 It is a kind of that the design method for moving bulge is arranged without shock-boundary
CN108019279A (en) * 2017-12-07 2018-05-11 中国人民解放军国防科技大学 A Design Method for Hypersonic Inlet
CN110414016A (en) * 2018-04-27 2019-11-05 中国航天科工飞航技术研究院(中国航天海鹰机电技术研究院) The Waverider geometry parameterization design method and system of ultrahigh speed pipeline transportation tool
CN109279043A (en) * 2018-10-23 2019-01-29 中国人民解放军国防科技大学 Design method of Von Karman waverider with low-speed airfoil
CN109598062B (en) * 2018-12-04 2022-12-02 中国人民解放军国防科技大学 Design method of osculating flow field waverider with variable wall surface pressure distribution rule
CN109598062A (en) * 2018-12-04 2019-04-09 中国人民解放军国防科技大学 Design method of osculating flow field waverider with variable wall surface pressure distribution rule
CN112765725A (en) * 2020-12-30 2021-05-07 四川京航天程科技发展有限公司 Analytic Riemann resolving method for multi-dimensional Euler equation
CN112765725B (en) * 2020-12-30 2023-04-07 四川京航天程科技发展有限公司 Analytic Riemann resolving method for multi-dimensional Euler equation
CN112949199A (en) * 2021-03-15 2021-06-11 中国科学院力学研究所 Method and system for optimizing longitudinal stability of power waverider
CN112948976A (en) * 2021-03-15 2021-06-11 中国科学院力学研究所 Plane front edge lifting body front edge line optimization method and system and lifting body
CN112949199B (en) * 2021-03-15 2024-02-02 中国科学院力学研究所 Method and system for optimizing longitudinal stability of power wave body
CN112948976B (en) * 2021-03-15 2024-02-02 中国科学院力学研究所 Plane leading edge lifting body leading edge line optimization method, system and lifting body
CN114313253A (en) * 2022-03-03 2022-04-12 中国空气动力研究与发展中心计算空气动力研究所 Aerodynamic layout and design method of high lift-drag ratio air-breathing hypersonic aircraft
CN114852366A (en) * 2022-05-30 2022-08-05 厦门大学 A Design Method for Dual Waverider Configurations with Simultaneous Controllable Shock Shape and Pressure Distribution
CN114852366B (en) * 2022-05-30 2024-07-05 厦门大学 Design method for shock wave shape and pressure distribution simultaneously controllable double waverider configuration

Also Published As

Publication number Publication date
CN105329462B (en) 2016-08-24

Similar Documents

Publication Publication Date Title
CN105329462B (en) Osculating flow field based on variable wall surface pressure regularity of distribution waverider forebody derived method for designing
CN104210672B (en) Hypersonic rider fuselage and inlet channel integrated design method
CN105697150B (en) A Design Method of Bulky Inlet with Consideration of Cone Guide and Close Waverider
CN106005475B (en) Hypersonic inside and outside flow integrated full Waverider aircraft method for designing
CN104192302B (en) Based on the Waverider method of designing around tip Feng karman curve gyro-rotor benchmark flow field
CN109250144B (en) Design method of kiss-cut cone waverider with direct controllable sweep angle and dihedral angle
CN105151316B (en) Design method of gliding-and-cruising two-stage waverider based on variable shock wave angles and osculating cone principle
CN104973266B (en) Based on glide-two-stage of the cruising Waverider method of designing of osculating cone theazy
CN109573093B (en) Design method of kiss-cut axisymmetric Von Karman waverider with low-speed airfoil
CN105667812B (en) Integrated design method of hypersonic vehicle precursor, inlet and wing waveriders
CN107672821B (en) Variable Mach number parallel wide-speed-range waverider aircraft design method based on osculating cone theory
CN109573092B (en) Design method of kiss-cut axisymmetric von Karman waverider
US20070181743A1 (en) Method for streamline traced external compression inlet
CN107514311A (en) Based on rotatable air intake duct/waverider forebody derived integrated design method in precursor shock wave
CN110304267B (en) Hypersonic aircraft design method and system
CN107310748B (en) Become the Design of Aerodynamic Configuration method of Angle of Shock Waves osculating flow field Waverider
CN115056998A (en) A design method for longitudinal segmented and graded compression of a cone-guided waverider precursor
CN110450963A (en) Method and system for integrated design of hypersonic aircraft body and inward turning inlet
CN108100291A (en) A Kiss Cut Waverider Design Method with a Given 3D Front Line
CN108595856A (en) Parametric design method of 'binary' waverider
CN104912667A (en) Design method of hypersonic speed internal-contraction air inlet channel carried out in steps
CN109279043B (en) Von Karman waverider design method fused with low-speed airfoil
CN111003196A (en) Full-wave-rider aircraft and design method and system thereof
CN107016199B (en) A Design Method of Shock-Free Boundary Layer Displacement Bulge
CN113942651A (en) A Novel Flight Control Device for SACCON Type Aircraft

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant