CN108280264A - 一种跨声速喷管设计方法 - Google Patents
一种跨声速喷管设计方法 Download PDFInfo
- Publication number
- CN108280264A CN108280264A CN201711445425.3A CN201711445425A CN108280264A CN 108280264 A CN108280264 A CN 108280264A CN 201711445425 A CN201711445425 A CN 201711445425A CN 108280264 A CN108280264 A CN 108280264A
- Authority
- CN
- China
- Prior art keywords
- section
- curve
- nozzle
- points
- wall
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/14—Pipes
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Lining And Supports For Tunnels (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
一种跨声速喷管设计方法,涉及风洞试验领域;包括如下步骤:步骤(一)、初设跨声速喷管外形;步骤(二)、构建喷管膨胀段初始外壁曲线;步骤(三)、对初始的膨胀段的外壁曲线中的圆锥曲线段进行优化;步骤(四)、对初始的膨胀段的外壁曲线中的三次曲线段进行优化;步骤(五)、得到优化后的膨胀段外壁曲线;步骤(六)、得到收缩段的外壁曲线和最终跨声速喷管外壁曲线;本发明减少了喷管膨胀段的长度,缩减空间,较少加工成本;且气流流过喷管到达喷管出口时,流场品质满足国军标先进指标。
Description
技术领域
本发明涉及一种风洞试验领域,特别是一种跨声速喷管设计方法。
背景技术
在跨声速范围内,流动是比较复杂的。流场中既有亚声速区域,又有超声速区域,而且流场经常不稳定,流动十分复杂。跨声速风洞必须解决试验段马赫数接近1时所出现的各种问题,技术比较复杂。在跨声速范围内,影响模型气动特性的主要相似准则是马赫数,所以跨声速喷管的出口流场品质非常重要。
目前,扩声速喷管主要利用抽气机除去驻室压力,是在声速截面前排除一部分流量,使其在试验段产生所需要的跨声速气流。为了防止气流拥塞,气动模型需要设计较小;由于需要抽气机吸除驻室压力,对抽气机的功率要求很大,导致较高的试验成本。
发明内容
本发明的目的在于克服现有技术的上述不足,提供一种跨声速喷管设计方法,减少了喷管膨胀段的长度,缩减空间,较少加工成本;且气流流过喷管到达喷管出口时,流场品质满足国军标先进指标。
本发明的上述目的是通过如下技术方案予以实现的:
一种跨声速喷管设计方法,包括如下步骤:
步骤(一)、初设跨声速喷管外形
包括收缩段和膨胀段;收缩段和膨胀段沿轴向首尾连接;连接处为喉道;喉道对应喷管外壁的位置为T点;且收缩段轴向远离膨胀段的一端为喷管入口端;膨胀段轴向远离收缩段的一端为喷管出口端;所述膨胀段的外壁包括三次曲线段、圆锥曲线段和Bézier曲线段;三次曲线段、圆锥曲线段和Bézier曲线段依次连接,形成膨胀段的外壁曲线;
步骤(二)、根据预先设定的喷管出口端直径和喷管喉道的位置,利用Bézier曲线r(u)和伯恩斯坦基函数Jn,i构建喷管膨胀段初始外壁曲线;
步骤(三)、对初始的膨胀段的外壁曲线中的圆锥曲线段进行优化;
建立跨声速喷管坐标系oxy;求得步骤(二)中Bézier曲线r(u)一阶导数最大点,即为圆锥曲线段和Bézier曲线段的连接点D点;得到圆锥曲线段和Bézier曲线段的连接点D点的坐标(xD,yD);以D点的坐标(xD,yD)为基点,根据Bézier曲线段在D点的一阶导数即为圆锥曲线段的斜率,根据斜率;生成优化后的圆锥曲线段;
步骤(四)、对初始的膨胀段的外壁曲线中的三次曲线段进行优化;
步骤(五)、依次连接步骤(四)得到的优化后的三次曲线段、步骤(三)得到的优化后的圆锥曲线段和步骤(二)得到的Bézier曲线段;得到优化后的膨胀段外壁曲线;
步骤(六)、根据步骤(四)得到的优化后的三次曲线段,获得T点下游的斜率;令T点下游的斜率与T点上游斜率相同;根据预设的喷管入口直径;得到收缩段的外壁曲线;连接收缩段和优化后的膨胀段外壁曲线,得到最终跨声速喷管外壁曲线。
在上述的一种跨声速喷管设计方法,所述的步骤(二)中,利用Bézier曲线r(u)和伯恩斯坦基函数Jn,i构建喷管膨胀段初始外壁曲线的方法为:
Bézier曲线r(u)采用特征多边形顶点的位置矢量与伯恩斯坦基函数的线性组合表达:
式中,n为Bézier曲线的次数,6≤n≤8,且n为正整数;
i为Bézier曲线征多项形顶点的序号,0≤i≤n;
u为Bézier曲线参数,0≤u≤1;
Vi为特征多边形顶点的位置矢量;
Jn,i为伯恩斯坦基函数;
式中,为组合数;
得到初始的膨胀段的外壁曲线。
在上述的一种跨声速喷管设计方法,所述的步骤(三)中,跨声速喷管坐标系oxy的建立方法为:以T点为坐标原点,x轴正方向为沿收缩段指向膨胀段方向;y轴正方向为竖直向上方向。
在上述的一种跨声速喷管设计方法,所述的步骤(四)中,对三次曲线段进行优化的具体方法为:
根据步骤(三)得到的圆锥曲线段,选取圆锥曲线段的端C点,得到C点的坐标(xC,yC);得到T点的坐标(xT,yT);以C点坐标(xC,yC)为起始点,以T点坐标(xT,yT)为终点;构建三次曲线段的方程,生成优化后的三次曲线段。
在上述的一种跨声速喷管设计方法,所述的步骤(四)中,所述三次曲线段的方程为:
式中,其中tanβC为C点下游处斜率,C点下游处斜率与圆锥曲线段的斜率相同。
本发明与现有技术相比具有如下优点:
(1)本发明采用步骤1中的Bézier曲线,能够减少喷管膨胀段的长度,缩减空间,较少加工成本;
(2)本发明采用步骤4的方法,实现了喷管收缩段和膨胀段连接处斜率连续,气流流过喷管收缩段和膨胀段时,气流无分离;
(3)本发明步骤5生成的喷管膨胀段曲线,使得气流流过喷管到达喷管出口时,流场品质满足国军标先进指标。
附图说明
图1为本发明跨声速喷管型面示意图;
图2为本发明跨声速喷管设计流程图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步详细的描述:
如图1所示为跨声速喷管型面示意图,由图可知,一种跨声速喷管设计方法,包括如下步骤:
步骤(一)、初设跨声速喷管外形
包括收缩段1和膨胀段2;收缩段1和膨胀段2沿轴向首尾连接;连接处为喉道;如图2所示为跨声速喷管设计流程图,由图可知,喉道对应喷管外壁的位置为T点;且收缩段1轴向远离膨胀段2的一端为喷管入口端;膨胀段2轴向远离收缩段1的一端为喷管出口端;所述膨胀段2的外壁包括三次曲线段2-1、圆锥曲线段2-2和Bézier曲线段2-3;三次曲线段2-1、圆锥曲线段2-2和Bézier曲线段2-3依次连接,形成膨胀段2的外壁曲线;
步骤(二)、根据预先设定的喷管出口端直径和喷管喉道的位置,利用Bézier曲线r(u)和伯恩斯坦基函数Jn,i构建喷管膨胀段2初始外壁曲线;只要给出给定特征多边形的顶点矢量,就能构造一条Bézier曲线。构造的Bézier曲线需要一阶、二阶导数连续,Bézier曲线按照一阶导数的单调性和二阶导数的正负分为两部分:第一部分一阶导数单调增加,二阶导数大于零;第二部分一阶导数单调递减,二阶导数小于零。一阶导数最大D点为Bézier两部分曲线的分界点。
利用Bézier曲线r(u)和伯恩斯坦基函数Jn,i构建喷管膨胀段2初始外壁曲线的方法为:
Bézier曲线r(u)采用特征多边形顶点的位置矢量与伯恩斯坦基函数的线性组合表达:
式中,n为Bézier曲线的次数,6≤n≤8,且n为正整数;
i为Bézier曲线征多项形顶点的序号,0≤i≤n;
u为Bézier曲线参数,0≤u≤1;
Vi为特征多边形顶点的位置矢量;
Jn,i为伯恩斯坦基函数;
式中,为组合数;
得到初始的膨胀段2的外壁曲线。
步骤(三)、对初始的膨胀段2的外壁曲线中的圆锥曲线段2-2进行优化;
建立跨声速喷管坐标系oxy;跨声速喷管坐标系oxy的建立方法为:以T点为坐标原点,x轴正方向为沿收缩段1指向膨胀段2方向;y轴正方向为竖直向上方向。求得步骤(二)中Bézier曲线r(u)一阶导数最大点,即为圆锥曲线段2-2和Bézier曲线段2-3的连接点D点;得到圆锥曲线段2-2和Bézier曲线段2-3的连接点D点的坐标(xD,yD);以D点的坐标(xD,yD)为基点,根据Bézier曲线段2-3在D点的一阶导数即为圆锥曲线段2-2的斜率,根据斜率;生成优化后的圆锥曲线段2-2;
步骤(四)、对初始的膨胀段2的外壁曲线中的三次曲线段2-1进行优化;
对三次曲线段2-1进行优化的具体方法为:
根据步骤(三)得到的圆锥曲线段2-2,选取圆锥曲线段2-2的端C点,得到C点的坐标(xC,yC);得到T点的坐标(xT,yT);以C点坐标(xC,yC)为起始点,以T点坐标(xT,yT)为终点;构建三次曲线段2-1的方程,生成优化后的三次曲线段2-1。
三次曲线段2-1的方程为:
式中,其中tanβC为C点下游处斜率,C点下游处斜率与圆锥曲线段2-2的斜率相同。
步骤(五)、依次连接步骤(四)得到的优化后的三次曲线段2-1、步骤(三)得到的优化后的圆锥曲线段2-2和步骤(二)得到的Bézier曲线段2-3;得到优化后的膨胀段2外壁曲线;
步骤(六)、根据步骤(四)得到的优化后的三次曲线段2-1,获得T点下游的斜率;令T点下游的斜率与T点上游斜率相同;根据预设的喷管入口直径;得到收缩段1的外壁曲线;连接收缩段1和优化后的膨胀段2外壁曲线,得到最终跨声速喷管外壁曲线。
通过数值计算,不断迭代,优化D点和C点的坐标取值,实时对喷管出口流场进行测量,当喷管出口流场品质达到国军标优秀指标时,迭代终止,完成跨声速喷管的型线的设计。
本发明喷管膨胀段圆锥曲线区域2-2在有些跨声速马赫数条件下可以取消,此时喷管膨胀段由喷管膨胀段三次曲线2-1和喷管膨胀段Bézier曲线2-3两部分组成,此种方式获得的喷管出口流场也满足国军标先进指标。
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
Claims (5)
1.一种跨声速喷管设计方法,其特征在于:包括如下步骤:
步骤(一)、初设跨声速喷管外形
包括收缩段(1)和膨胀段(2);收缩段(1)和膨胀段(2)沿轴向首尾连接;连接处为喉道;喉道对应喷管外壁的位置为T点;且收缩段(1)轴向远离膨胀段(2)的一端为喷管入口端;膨胀段(2)轴向远离收缩段(1)的一端为喷管出口端;所述膨胀段(2)的外壁包括三次曲线段(2-1)、圆锥曲线段(2-2)和Bézier曲线段(2-3);三次曲线段(2-1)、圆锥曲线段(2-2)和Bézier曲线段(2-3)依次连接,形成膨胀段(2)的外壁曲线;
步骤(二)、根据预先设定的喷管出口端直径和喷管喉道的位置,利用Bézier曲线r(u)和伯恩斯坦基函数Jn,i构建喷管膨胀段(2)初始外壁曲线;
步骤(三)、对初始的膨胀段(2)的外壁曲线中的圆锥曲线段(2-2)进行优化;
建立跨声速喷管坐标系oxy;求得步骤(二)中Bézier曲线r(u)一阶导数最大点,即为圆锥曲线段(2-2)和Bézier曲线段(2-3)的连接点D点;得到圆锥曲线段(2-2)和Bézier曲线段(2-3)的连接点D点的坐标(xD,yD);以D点的坐标(xD,yD)为基点,根据Bézier曲线段(2-3)在D点的一阶导数即为圆锥曲线段(2-2)的斜率,根据斜率;生成优化后的圆锥曲线段(2-2);
步骤(四)、对初始的膨胀段(2)的外壁曲线中的三次曲线段(2-1)进行优化;
步骤(五)、依次连接步骤(四)得到的优化后的三次曲线段(2-1)、步骤(三)得到的优化后的圆锥曲线段(2-2)和步骤(二)得到的Bézier曲线段(2-3);得到优化后的膨胀段(2)外壁曲线;
步骤(六)、根据步骤(四)得到的优化后的三次曲线段(2-1),获得T点下游的斜率;令T点下游的斜率与T点上游斜率相同;根据预设的喷管入口直径;得到收缩段(1)的外壁曲线;连接收缩段(1)和优化后的膨胀段(2)外壁曲线,得到最终跨声速喷管外壁曲线。
2.根据权利要求1所述的一种跨声速喷管设计方法,其特征在于:所述的步骤(二)中,利用Bézier曲线r(u)和伯恩斯坦基函数Jn,i构建喷管膨胀段(2)初始外壁曲线的方法为:
Bézier曲线r(u)采用特征多边形顶点的位置矢量与伯恩斯坦基函数的线性组合表达:
式中,n为Bézier曲线的次数,6≤n≤8,且n为正整数;
i为Bézier曲线征多项形顶点的序号,0≤i≤n;
u为Bézier曲线参数,0≤u≤1;
Vi为特征多边形顶点的位置矢量;
Jn,i为伯恩斯坦基函数;
式中,为组合数;
得到初始的膨胀段(2)的外壁曲线。
3.根据权利要求1所述的一种跨声速喷管设计方法,其特征在于:所述的步骤(三)中,跨声速喷管坐标系oxy的建立方法为:以T点为坐标原点,x轴正方向为沿收缩段(1)指向膨胀段(2)方向;y轴正方向为竖直向上方向。
4.根据权利要求1所述的一种跨声速喷管设计方法,其特征在于:所述的步骤(四)中,对三次曲线段(2-1)进行优化的具体方法为:
根据步骤(三)得到的圆锥曲线段(2-2),选取圆锥曲线段(2-2)的端C点,得到C点的坐标(xC,yC);得到T点的坐标(xT,yT);以C点坐标(xC,yC)为起始点,以T点坐标(xT,yT)为终点;构建三次曲线段(2-1)的方程,生成优化后的三次曲线段(2-1)。
5.根据权利要求4所述的一种跨声速喷管设计方法,其特征在于:所述的步骤(四)中,所述三次曲线段(2-1)的方程为:
式中,其中tanβC为C点下游处斜率,C点下游处斜率与圆锥曲线段(2-2)的斜率相同。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711445425.3A CN108280264B (zh) | 2017-12-27 | 2017-12-27 | 一种跨声速喷管设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711445425.3A CN108280264B (zh) | 2017-12-27 | 2017-12-27 | 一种跨声速喷管设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108280264A true CN108280264A (zh) | 2018-07-13 |
CN108280264B CN108280264B (zh) | 2021-04-13 |
Family
ID=62802353
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711445425.3A Active CN108280264B (zh) | 2017-12-27 | 2017-12-27 | 一种跨声速喷管设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108280264B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114329847A (zh) * | 2022-01-06 | 2022-04-12 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种“双钟型”收缩扩张喷管标模设计方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101477710A (zh) * | 2009-01-20 | 2009-07-08 | 北京航空航天大学 | 一种高超声速飞行器机身推进一体化外形建模方法 |
EP2433716A1 (en) * | 2010-09-22 | 2012-03-28 | Hexagon Technology Center GmbH | Surface spraying device with a nozzle control mechanism and a corresponding method |
CN103291495A (zh) * | 2013-05-21 | 2013-09-11 | 南京航空航天大学 | 超声速/高超声速飞行器发动机过膨胀喷管旁路式装置 |
CN104316287A (zh) * | 2014-10-24 | 2015-01-28 | 中国人民解放军国防科学技术大学 | 二维变马赫数喷管及使用该喷管的超声速变马赫数风洞 |
KR101541072B1 (ko) * | 2014-04-30 | 2015-07-31 | 조선대학교산학협력단 | 표층가열 재생공법 아스팔트 포장장치용 가열판의 최적 설계 방법 |
CN105676640A (zh) * | 2016-01-13 | 2016-06-15 | 南京航空航天大学 | 基于贝塞尔曲线的涡扇发动机加速过程控制律设计方法 |
CN106813890A (zh) * | 2016-12-19 | 2017-06-09 | 中国航天空气动力技术研究院 | 一种用于风洞跨声速试验前室超压的预判方法 |
CN106837549A (zh) * | 2017-02-06 | 2017-06-13 | 厦门大学 | 内并联式高超声速双通道进气道的设计方法 |
-
2017
- 2017-12-27 CN CN201711445425.3A patent/CN108280264B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101477710A (zh) * | 2009-01-20 | 2009-07-08 | 北京航空航天大学 | 一种高超声速飞行器机身推进一体化外形建模方法 |
EP2433716A1 (en) * | 2010-09-22 | 2012-03-28 | Hexagon Technology Center GmbH | Surface spraying device with a nozzle control mechanism and a corresponding method |
CN103291495A (zh) * | 2013-05-21 | 2013-09-11 | 南京航空航天大学 | 超声速/高超声速飞行器发动机过膨胀喷管旁路式装置 |
KR101541072B1 (ko) * | 2014-04-30 | 2015-07-31 | 조선대학교산학협력단 | 표층가열 재생공법 아스팔트 포장장치용 가열판의 최적 설계 방법 |
CN104316287A (zh) * | 2014-10-24 | 2015-01-28 | 中国人民解放军国防科学技术大学 | 二维变马赫数喷管及使用该喷管的超声速变马赫数风洞 |
CN105676640A (zh) * | 2016-01-13 | 2016-06-15 | 南京航空航天大学 | 基于贝塞尔曲线的涡扇发动机加速过程控制律设计方法 |
CN106813890A (zh) * | 2016-12-19 | 2017-06-09 | 中国航天空气动力技术研究院 | 一种用于风洞跨声速试验前室超压的预判方法 |
CN106837549A (zh) * | 2017-02-06 | 2017-06-13 | 厦门大学 | 内并联式高超声速双通道进气道的设计方法 |
Non-Patent Citations (2)
Title |
---|
宋红超: "基于离散型伴随方法的单边膨胀喷管优化设计研究", 《工程热物理学报》 * |
黄炳修: "高超音速风洞轴对称喷管设计方法的比较", 《第八届全国实验流体力学学术会议论文集》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114329847A (zh) * | 2022-01-06 | 2022-04-12 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种“双钟型”收缩扩张喷管标模设计方法 |
CN114329847B (zh) * | 2022-01-06 | 2023-03-21 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种“双钟型”收缩扩张喷管标模设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108280264B (zh) | 2021-04-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108195544B (zh) | 一种脉冲型风洞串列喷管 | |
CN108846224B (zh) | 一种超声速流道设计方法及装置 | |
CN107089341B (zh) | 与飞行器一体化的高超声速进气道外压缩面设计方法 | |
CN109927917A (zh) | 一种超声速飞行器内转式乘波前体进气道一体化设计方法 | |
CN108999725B (zh) | 一种带双钟型引射套管的引射喷管 | |
CN102998084B (zh) | 高超声速喷管结构 | |
CN204495535U (zh) | 1.2米量级风洞用动力模拟引射器 | |
CN208310917U (zh) | 一种解决s弯喷管与涡扇发动机匹配问题的转接段结构 | |
CN105043711B (zh) | 一种兼容多喷管的风洞扩压器及风洞扩压方法 | |
CN108757217B (zh) | 一种双钟型膨胀偏转喷管 | |
CN110284994A (zh) | 一种基于喉道偏移式气动矢量喷管的并联式推力矢量排气系统 | |
CN109723570A (zh) | 具有非对称后体型面的喉道偏移式气动矢量喷管 | |
CN102434315A (zh) | 旁路式双喉道无源矢量喷管 | |
CN112035952B (zh) | 一种模拟飞行器外流的引射喷管实验装置的设计方法 | |
CN105138787B (zh) | 基于特征线追踪的超声速流场设计方法 | |
CN107120210A (zh) | 一种超音速喷管 | |
CN105329462A (zh) | 基于可变壁面压力分布规律的吻切流场乘波前体设计方法 | |
CN107013334A (zh) | 一种双燃烧室超燃冲压发动机进气道及进气控制方法 | |
CN105205220B (zh) | 一种高超声速内转式进气道的内通道设计方法 | |
CN108280264A (zh) | 一种跨声速喷管设计方法 | |
CN103678774B (zh) | 考虑进口参数非均匀的超声速推力喷管设计方法 | |
CN110633522A (zh) | 一种基于最大推力理论的超声速推力喷管反设计方法 | |
CN106837601B (zh) | 带有侧向膨胀的喉道偏移式气动矢量喷管 | |
CN108563896A (zh) | 一种提高火箭发动机喷管性能的扩张段型面设计方法 | |
CN105464838B (zh) | 用于被动推力导向和羽流偏转的方法和装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |