CN108304611B - 一种给定三维前缘线的锥导乘波体设计方法 - Google Patents
一种给定三维前缘线的锥导乘波体设计方法 Download PDFInfo
- Publication number
- CN108304611B CN108304611B CN201711428331.5A CN201711428331A CN108304611B CN 108304611 B CN108304611 B CN 108304611B CN 201711428331 A CN201711428331 A CN 201711428331A CN 108304611 B CN108304611 B CN 108304611B
- Authority
- CN
- China
- Prior art keywords
- point
- waverider
- dimensional
- costa
- shock wave
- 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.)
- Active
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/15—Vehicle, aircraft or watercraft design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Aviation & Aerospace Engineering (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Structures Of Non-Positive Displacement Pumps (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明提供一种给定三维前缘线的锥导乘波体设计方法,解决现有的设计方法不能直接指定三维前缘线的问题。该方法通过给定三维前缘线以及流场参数,根据三维前缘线计算对应的激波出口型线,完成吻切面流场求解,最终生成对应的乘波体。本发明基于锥导理论,在给定三维前缘线和流场参数的情况下设计乘波体,不需要考虑激波型线的影响,直接控制乘波体的三维前缘线,实现了工程上对三维前缘线有要求乘波体的按需设计。
Description
技术领域
本发明属于高超声速飞行器技术领域,特别涉及一种基于给定三维前缘线的锥导乘波体设计方法。
背景技术
乘波体设计是当前国际上高超声速飞行器气动布局研究的重点和热点之一。当前设计方法主要分为两类:一种是给定激波生成体的设计理论,另一种是吻切设计理论(包括吻切锥和吻切轴对称)。前者是在基准流场中沿给定的前缘型线追踪流线获得乘波体下压缩面,该基准流场可以是轴对称或非轴对称。参见李永洲,孙迪,张堃元.前后缘型线同时可控的乘波体设计[J].航空学报,2017,38(1):120153.。
现有的锥导乘波体设计方法通常是给定圆锥激波和乘波体三个典型型线(前缘线水平投影型线、上表面后缘线和下表面后缘线)之一,通过计算流场流线追踪生成乘波体外形。尽管该方法可以同时控制乘波体某一特征型线和激波出口型线,从另一角度说却限制了对乘波体外形的直接控制程度,无法根据工程需求对乘波体外形进行更直接的控制设计。
发明内容
本发明的目的在于提供一种给定三维前缘线的锥导乘波体设计方法,解决现有的设计方法不能直接指定三维前缘线的问题。本设计方法可以根据工程实际要求中给出对三维前缘线,通过该三维前缘线来生成乘波体,对乘波体的外形进行更直接的控制,且无需考虑激波出口型线形状。
为实现上述技术目的,本发明采用的技术方案如下:
一种给定三维前缘线的锥导乘波体设计方法,包括以下步骤:
S1.给定三维前缘线以及流场参数;
其中流场参数包括马赫数Ma,激波角β;乘波体长度L和乘波体宽度W在给定三维前缘线时已经确定。给定的三维前缘线是通过锥导法生成的三维前缘线。
S2根据三维前缘线计算对应的激波出口型线。
首先取三维前缘线其纵向对称面上点作原点,沿乘波体纵向为x方向,垂直方向为y方向,展向方向为z方向,建立乘波体坐标系。
给定三维前缘线,即三维前缘线上在纵向对称面上的点P以及三维前缘线其两端的端点P3、P4是已知的,P点在乘波体底面的投影为P1点。
激波出口型线在纵向对称面上的点为P2点,由公式(1)计算P1点和P2点间的距离:
D=Ltanβ (1)
其中L为P点至乘波体底面的距离,β为激波角,D为P1点和P2点间的距离。
由公式(1)得到P1点和P2点的距离,即可知激波点P2的位置坐标。
对于锥导乘波体,其激波出口型线是圆弧即激波出口型线为圆的一部分,即可根据P2、P3、P4三点确定该三点对应的激波圆的圆心P’点和激波出口型线,并由P’点作平行于x轴的直线,过P2点作一条与该直线夹角为β的直线,两线交点即为圆锥激波顶点。
S3流场求解;
根据S2求得的圆锥激波顶点,通过求解Taylor-Maccoll流场控制方程获得圆锥激波的基准锥的锥半角和对应流场。具体的求解Taylor-Maccoll流场控制方程参见AndersonJ.D.Fundamental of aerodynamics[M].3rd edition.McGraw-Hill Companies,2001.。
对三维前缘线进行离散获得一系列离散点,这里采用等参数离散,具体的离散方法参见毋河海.曲线离散表达的优化研究.全国地图学与gis学术研讨会,2008.。对三维前缘线上的每一离散点进行流线追踪获得各离散点对应的下表面流线,流线追踪方法可参见丁峰.高超声速滑翔-巡航两级乘波设计方法研究[D];国防科学技术大学,2012.。各下表面流线在乘波体底面上的点即各下表面流线对应的下表面后缘线上的点。三维前缘线在乘波体底面的投影即为上表面后缘线。对于三维前缘线上的每一离散点利用自由流线法生成对应的上表面流线,具体方法可参见王庆文.基于吻切理论的两级乘波体设计[D]:[硕士].长沙:国防科学技术大学;2015。
S4生成乘波体外形。
一系列下表面流线平滑连接构成乘波体下表面;一系列上表面流线平滑连接构成乘波体上表面;一系列下表面后缘线上的点平滑连接构成乘波体下表面后缘线;最后对乘波体底部进行封闭,则乘波体设计完成,生成乘波体外形如图4所示。
本发明具有以下技术效果:
本发明基于锥导理论,在给定三维前缘线和流场参数的情况下设计乘波体,不需要考虑激波型线的影响,直接控制乘波体的三维前缘线,实现了工程上对三维前缘线有要求乘波体的按需设计。
附图说明
图1是乘波体部分几何参数示意图,其中图1(a)为乘波体侧视图,图1(b)为乘波体底面后视图;
图2是求解激波型线示意图;
图3是基于锥导理论的乘波体设计原理图;
图4是本发明最终生成的乘波体外形图;
图5为采用传统锥导法生成的乘波体外形图;
图6为激波出口型线对比图;
图7为前缘线水平投影型线对比图;
图8为下表面后缘线对比图。
具体实施方式
以下将结合具体实施例和说明书附图对本发明做进一步详细说明。
首先对待设计的乘波体中的各典型型线以及表面进行介绍,参照图1为乘波体部分几何参数示意图,图1中Ma为来流马赫数,L为乘波体长度,H为乘波体高度,W为乘波体翼展,点1、2、3、4分别为乘波体的三维前缘线、上表面出口型线、下表面出口型线、激波出口型线的中点,点5和点6为乘波体的三维前缘线、上表面出口型线、下表面出口型线和激波出口型线这四条典型型线相交的两个交点,曲线5-1-6为乘波体的三维前缘线,其中点5和点6分别为三维前缘线的两个端点;曲线5-2-6为乘波体的上表面出口型线,其中点5和点6同时也为上表面出口型线的两个端点;曲线5-3-6为乘波体的下表面出口型线,其中点5和点6同时也为下表面出口型线的两个端点;曲线5-4-6为激波出口型线,其中点5和点6同时也为激波出口型线的两个端点;曲线5-1-6和曲线5-2-6所形成的曲面1-5-2-6为乘波体上表面,曲线5-1-6和曲线5-3-6所形成的曲面1-5-3-6为乘波体下表面,曲线5-2-6和曲线5-3-6所在的平面区域5-3-6-2为乘波体底面。
基于给定三维前缘线的锥导乘波体设计方法,包括以下步骤:
S1.给定乘波体的三维前缘线以及流场参数;
采用本发明方法进行乘波体设计时,所需要输入的典型型线仅有三维前缘线,不需要考虑其他激波型线的设计选择。
在实际应用中,设计师是根据乘波体飞行条件、尺寸约束和性能要求等给出满足设计要求的三维前缘线以及流场参数。在本实施例中:给定的三维前缘线是通过锥导法生成的三维前缘线。
其中流场参数包括马赫数Ma,激波角β,乘波体高度H,乘波体长度L和乘波体宽度W在给定三维前缘线时已经确定。
S2根据三维前缘线计算对应的激波出口型线。
首先取三维前缘线其纵向对称面上点作原点,沿乘波体纵向为x方向,垂直方向为y方向,展向方向为z方向,建立乘波体坐标系。
给定三维前缘线,即三维前缘线上在纵向对称面上的点P以及三维前缘线其两端的端点P3、P4是已知的,如图2所示。P点在乘波体底面的投影为P1点。
由公式(1)计算激波出口型线在纵向对称面上的点P2
D=Ltanβ (1)
其中L为P点至乘波体底面的距离,β为激波角,D为P1点和P2点间的距离。
由公式(1)得到P1点和P2点的距离,即可知激波点P2的位置坐标。
对于锥导乘波体,其激波出口型线是圆弧即激波出口型线为圆的一部分,即可根据P2、P3、P4三点确定该三点对应的激波圆的圆心P’点和激波出口型线,并由P’点作平行于x轴的直线,过P2点作直线与该直线的夹角为β,两线交点即为圆锥激波顶点。
S3吻切面流场求解;
参照图3,为基于锥导理论的乘波体设计原理图。根据S2求得的圆锥激波顶点,通过求解Taylor-Maccoll流场控制方程获得圆锥激波的基准锥的锥半角和对应流场。具体的求解Taylor-Maccoll流场控制方程参见Anderson J.D.Fundamental of aerodynamics[M].3rd edition.McGraw-Hill Companies,2001.。
对三维前缘线进行离散获得一系列离散点,这里采用等参数离散,具体的离散方法参见毋河海.曲线离散表达的优化研究.全国地图学与gis学术研讨会,2008.。对三维前缘线上的每一离散点进行流线追踪获得各离散点对应的下表面流线,流线追踪方法可参见丁峰.高超声速滑翔-巡航两级乘波设计方法研究[D];国防科学技术大学,2012.。各下表面流线在乘波体底面上的点即各下表面流线对应的下表面后缘线上的点。三维前缘线在乘波体底面的投影即为上表面后缘线。对于三维前缘线上的每一离散点利用自由流线法生成对应的上表面流线,具体方法可参见王庆文.基于吻切理论的两级乘波体设计[D]:[硕士].长沙:国防科学技术大学;2015。
S4生成乘波体外形。
一系列下表面流线平滑连接构成乘波体下表面;一系列上表面流线平滑连接构成乘波体上表面;一系列下表面后缘线上的点平滑连接构成乘波体下表面后缘线;最后对乘波体底部进行封闭,则乘波体设计完成,生成乘波体外形如图4所示。
为验证本发明的正确性和有效性,将采用本发明方法设计出的乘波体与传统锥导法设计的乘波体进行对比。参照图5为采用传统锥导法生成的参考乘波体外形。本发明采用的三维前缘线采用的即是传统锥导法中的三维前缘线,将利用本发明方法生成的乘波体与利用传统锥导法生成的乘波体进行对比,具体地,将两种方法生成的乘波体上的激波出口型线、前缘线水平投影型线、下表面后缘线进行对比,分别如图6、图7、图8所示,其中黑色线条为采用传统锥导法生成的参考乘波体上的参考型线,黑色点为本发明生成对应型线上的离散点。经对比两种方法生成的乘波体在各典型型线上均吻合,两种方法生成的乘波体是相同的。由此可见,本发明仅通过给定三维前缘线就能够实现乘波体的设计,克服了传统锥导法需要给定多种典型型线,输入条件较为繁杂的缺陷。本发明方法直接基于给定的三维前缘线完成乘波体设计,方法正确且有效。
以上仅是本发明的优选实施方式,本发明的保护范围并不仅局限于上述实施例,凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理前提下的若干改进和润饰,应视为本发明的保护范围。
Claims (2)
1.一种给定三维前缘线的锥导乘波体设计方法,其特征在于,包括以下步骤:
S1.给定三维前缘线以及流场参数;
其中流场参数包括马赫数Ma,激波角β,乘波体长度L和乘波体宽度W在给定三维前缘线时已经确定;
S2根据三维前缘线计算对应的激波出口型线;
首先取三维前缘线其纵向对称面上点作为原点,沿乘波体纵向为x方向,垂直方向为y方向,展向方向为z方向,建立乘波体坐标系;
给定三维前缘线即三维前缘线其在纵向对称面上的点P以及三维前缘线其两端的端点P3、P4是已知的,P点在乘波体底面的投影为P1点;
激波出口型线在纵向对称面上的点为P2点,由公式(1)计算P1点和P2点间的距离:
D=L tanβ (1)
其中L为P点至乘波体底面的距离,β为激波角,D为P1点和P2点间的距离;
由公式(1)得到P1点和P2点间的距离D,即可知激波点P2的位置坐标;
对于锥导乘波体,其激波出口型线是圆弧即激波出口型线为圆的一部分,即可根据P2、P3、P4三点确定该三点对应的激波圆的圆心P’点和激波出口型线;并由P’点作平行于x轴的直线,过P2点作一条与该直线夹角为β的直线,两直线的交点即为圆锥激波顶点;
S3吻切面流场求解;
根据S2求得的圆锥激波顶点,通过求解Taylor-Maccoll流场控制方程获得圆锥激波的基准锥的锥半角和对应流场;
对三维前缘线进行离散获得一系列离散点,对三维前缘线上的每一离散点进行流线追踪获得各离散点对应的下表面流线;各下表面流线在乘波体底面上的点即各下表面流线对应的下表面后缘线上的点;三维前缘线在乘波体底面的投影即为上表面后缘线;对于三维前缘线上的每一离散点利用自由流线法生成各自对应的上表面流线;
S4根据流场求解结果生成乘波体外形;
利用S3中得到的一系列下表面流线平滑连接构成乘波体下表面;利用S3中得到的一系列上表面流线平滑连接构成乘波体上表面;利用S3中得到的一系列下表面后缘线上的点平滑连接构成乘波体下表面后缘线;最后对乘波体底部进行封闭,则乘波体设计完成。
2.根据权利要求1所述的给定三维前缘线的锥导乘波体设计方法,其特征在于,S3中采用等参数离散方法对三维前缘线进行离散。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711428331.5A CN108304611B (zh) | 2017-12-26 | 2017-12-26 | 一种给定三维前缘线的锥导乘波体设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711428331.5A CN108304611B (zh) | 2017-12-26 | 2017-12-26 | 一种给定三维前缘线的锥导乘波体设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108304611A CN108304611A (zh) | 2018-07-20 |
CN108304611B true CN108304611B (zh) | 2019-01-11 |
Family
ID=62870693
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711428331.5A Active CN108304611B (zh) | 2017-12-26 | 2017-12-26 | 一种给定三维前缘线的锥导乘波体设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108304611B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112298599B (zh) * | 2020-11-02 | 2022-05-03 | 厦门大学 | 基于弯曲激波理论的全三维乘波体反设计方法 |
CN112948976B (zh) * | 2021-03-15 | 2024-02-02 | 中国科学院力学研究所 | 一种平面前缘升力体前缘线优化方法、系统及升力体 |
CN112949199B (zh) * | 2021-03-15 | 2024-02-02 | 中国科学院力学研究所 | 一种幂次乘波体的纵向稳定性优化方法及系统 |
CN116611176B (zh) * | 2023-07-21 | 2023-09-15 | 中国空气动力研究与发展中心空天技术研究所 | 一种三维前缘可定制的流线追踪类乘波体设计方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103029830A (zh) * | 2012-12-20 | 2013-04-10 | 中国科学院力学研究所 | 一种双乘波体对拼吸气式高超飞行器前体及其设计方法 |
CN103106295A (zh) * | 2012-12-31 | 2013-05-15 | 中国人民解放军国防科学技术大学 | 一种用于数值模拟的非定常激波生成方法 |
CN103770935A (zh) * | 2013-12-13 | 2014-05-07 | 中国航天空气动力技术研究院 | 乘波体外形设计方法 |
CN104210672A (zh) * | 2014-07-18 | 2014-12-17 | 中国人民解放军国防科学技术大学 | 高超声速乘波机身与进气道一体化设计方法 |
CN104973266A (zh) * | 2015-07-16 | 2015-10-14 | 中国人民解放军国防科学技术大学 | 基于吻切锥理论的滑翔-巡航两级乘波体设计方法 |
CN105151316A (zh) * | 2015-09-24 | 2015-12-16 | 中国人民解放军国防科学技术大学 | 基于变激波角吻切锥理论的滑翔-巡航两级乘波体设计方法 |
CN105667812A (zh) * | 2016-01-29 | 2016-06-15 | 中国人民解放军国防科学技术大学 | 高超声速飞行器前体、进气道和机翼乘波一体化设计方法 |
CN106043737A (zh) * | 2016-06-29 | 2016-10-26 | 中国人民解放军国防科学技术大学 | 一种“等物面‑变马赫数”宽速域乘波飞行器设计方法 |
CN106043738A (zh) * | 2016-06-29 | 2016-10-26 | 中国人民解放军国防科学技术大学 | 一种等激波流场‑变马赫数宽速域乘波飞行器设计方法 |
CN106250597A (zh) * | 2016-07-26 | 2016-12-21 | 厦门大学 | 一种完全流向抽吸的三维内转进气道设计方法 |
-
2017
- 2017-12-26 CN CN201711428331.5A patent/CN108304611B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103029830A (zh) * | 2012-12-20 | 2013-04-10 | 中国科学院力学研究所 | 一种双乘波体对拼吸气式高超飞行器前体及其设计方法 |
CN103106295A (zh) * | 2012-12-31 | 2013-05-15 | 中国人民解放军国防科学技术大学 | 一种用于数值模拟的非定常激波生成方法 |
CN103770935A (zh) * | 2013-12-13 | 2014-05-07 | 中国航天空气动力技术研究院 | 乘波体外形设计方法 |
CN104210672A (zh) * | 2014-07-18 | 2014-12-17 | 中国人民解放军国防科学技术大学 | 高超声速乘波机身与进气道一体化设计方法 |
CN104973266A (zh) * | 2015-07-16 | 2015-10-14 | 中国人民解放军国防科学技术大学 | 基于吻切锥理论的滑翔-巡航两级乘波体设计方法 |
CN105151316A (zh) * | 2015-09-24 | 2015-12-16 | 中国人民解放军国防科学技术大学 | 基于变激波角吻切锥理论的滑翔-巡航两级乘波体设计方法 |
CN105667812A (zh) * | 2016-01-29 | 2016-06-15 | 中国人民解放军国防科学技术大学 | 高超声速飞行器前体、进气道和机翼乘波一体化设计方法 |
CN106043737A (zh) * | 2016-06-29 | 2016-10-26 | 中国人民解放军国防科学技术大学 | 一种“等物面‑变马赫数”宽速域乘波飞行器设计方法 |
CN106043738A (zh) * | 2016-06-29 | 2016-10-26 | 中国人民解放军国防科学技术大学 | 一种等激波流场‑变马赫数宽速域乘波飞行器设计方法 |
CN106250597A (zh) * | 2016-07-26 | 2016-12-21 | 厦门大学 | 一种完全流向抽吸的三维内转进气道设计方法 |
Non-Patent Citations (4)
Title |
---|
"High-temperature gas effects on aerodynamic characteristics of waverider";liu Jun等;《Chinese Journal of Aeronautics》;20141224;第57-65页 |
"Influence of surface pressure distribution of basic flow field on shape and performance of waverider";Feng Ding 等;《Acta Astronautica》;20141211;第62-78页 |
"基于组合控制的高超声速乘波构型飞行器BTT与STT控制方案对比";李健 等;《2007中国控制与决策学术年会论文集》;20071231;第89-92页 |
"锥形流乘波体优化设计研究";耿永兵 等;《航空学报》;20060131;第27卷(第1期);第23-28页 |
Also Published As
Publication number | Publication date |
---|---|
CN108304611A (zh) | 2018-07-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108304611B (zh) | 一种给定三维前缘线的锥导乘波体设计方法 | |
CN108100291B (zh) | 一种给定三维前缘线的吻切乘波体设计方法 | |
WO2019091057A1 (zh) | 基于密切锥理论的定平面乘波体设计方法 | |
CN105059530B (zh) | 一种后掠角可控的尖顶点密切锥乘波体 | |
Rodi | The osculating flowfield method of waverider geometry generation | |
CN106428620B (zh) | 一种大容积高升阻比脊形吻切锥乘波体设计方法 | |
CN105173116A (zh) | 高超声速飞行器密切曲面乘波体设计方法 | |
KR101308348B1 (ko) | 날개-장착된 견인 프로펠러들/로터들에 의해 구동되는 비행체의 날개의 유도 또는 전체 항력을 감소시키는 최적의 날개 평면도들 | |
CN110525679B (zh) | 高超声速嵌入式乘波体设计方法 | |
CN109573092A (zh) | 吻切轴对称冯卡门乘波体设计方法 | |
CN109573093A (zh) | 融合低速翼型的吻切轴对称冯卡门乘波体设计方法 | |
CN106394874A (zh) | 一种双后掠布局的乘波体 | |
Rodi | Vortex lift waverider configurations | |
CN105059531B (zh) | 一种后掠角可控的曲线头部密切锥乘波体 | |
CN108763806B (zh) | 高速列车过特长双线隧道时列车及隧道表面所受静压的数值模拟分析方法 | |
Liang et al. | Detached-eddy simulation of wing-tip vortex in the near field of NACA 0015 airfoil | |
CN103303461B (zh) | 一种飞机天线整流罩 | |
CN112926132A (zh) | 一种考虑三维效应影响下的固定翼翼型气动外形设计方法 | |
CN206155818U (zh) | 一种涡效应增升的乘波体 | |
CN107016199A (zh) | 一种无激波边界层排移鼓包的设计方法 | |
CN108629108B (zh) | 一种机身机翼部件可控的参数化乘波体设计方法 | |
Takashima et al. | Wedge-cone waverider configuration for engine-airframe interaction | |
Rosenfeld et al. | Numerical study of the skin friction on a spheroid at incidence | |
CN203428019U (zh) | 一种飞机天线整流罩 | |
CN105752315A (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 |