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 PDFInfo
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 239000002243 precursor Substances 0.000 title abstract description 33
- 230000035939 shock Effects 0.000 claims abstract description 53
- 230000007704 transition Effects 0.000 claims description 32
- 230000000630 rising effect Effects 0.000 claims description 3
- 235000013399 edible fruits Nutrition 0.000 claims 3
- 230000010339 dilation Effects 0.000 claims 2
- 230000015572 biosynthetic process Effects 0.000 claims 1
- 238000003786 synthesis reaction Methods 0.000 claims 1
- 230000006835 compression Effects 0.000 description 4
- 238000007906 compression Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 230000007123 defense Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005290 field theory Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000010899 nucleation Methods 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64F—GROUND 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/00—Designing, 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
技术领域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)
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)
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)
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 |
-
2015
- 2015-11-16 CN CN201510785607.XA patent/CN105329462B/en active Active
Patent Citations (4)
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)
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 |