CN107766673B - 一种半径可控的参数化三维前缘钝化设计方法 - Google Patents

一种半径可控的参数化三维前缘钝化设计方法 Download PDF

Info

Publication number
CN107766673B
CN107766673B CN201711102724.7A CN201711102724A CN107766673B CN 107766673 B CN107766673 B CN 107766673B CN 201711102724 A CN201711102724 A CN 201711102724A CN 107766673 B CN107766673 B CN 107766673B
Authority
CN
China
Prior art keywords
leading edge
curve
waverider
passivation
point
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
Application number
CN201711102724.7A
Other languages
English (en)
Other versions
CN107766673A (zh
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 CN201711102724.7A priority Critical patent/CN107766673B/zh
Publication of CN107766673A publication Critical patent/CN107766673A/zh
Application granted granted Critical
Publication of CN107766673B publication Critical patent/CN107766673B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/15—Vehicle, aircraft or watercraft design
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Automation & Control Theory (AREA)
  • Numerical Control (AREA)
  • Prostheses (AREA)

Abstract

一种半径可控的参数化三维前缘钝化设计方法，首先给定初始乘波体外形，然后对初始乘波体外形中的尖锐前缘进行钝化，包括以下步骤：先将乘波体上表面抬升一定高度，接着获取抬升后的乘波体前缘点P2的插值条件，给出抬升后的钝化前缘曲线上不同横向站位处的钝化半径，设置抬升后的乘波体前缘点P2以及乘波体的上表面沿法向抬升高度h；最后利用3次Bézier曲线构造钝化前缘曲线，然后放样得到钝化前缘曲面。本发明在保证钝化前缘曲线与上、下表面曲线光滑连接的同时，还可以根据热防护需求指定各前缘点处的钝化半径，有效地实现了热防护与气动外形的多学科耦合设计。

Description

一种半径可控的参数化三维前缘钝化设计方法
技术领域
本发明涉及高超声速飞行器技术领域，尤其涉及一种半径可控的参数化三维前缘钝化设计方法。
背景技术
在高超声速飞行时，空气受到剧烈压缩而产生的气动加热使得飞行器面临的巨大的热防护压力，因此高超声速飞行器在工程应用中必须要对存在的尖锐前缘进行钝化处理。在旋成体、升力体、翼身融合体等常见高超声速飞行器气动构型之中，乘波体具有升阻比大、设计过程简单等优点，是一种应用前景广阔的高超声速飞行器气动构型方案。但乘波体的设计原理决定了其外形必然带有尖锐前缘，给热防护和加工制造带来一定困难，在设计时有必要对其前缘进行钝化处理。
目前，虽然已有相应三维钝化方法可使钝化前缘与飞行器表面光滑连接，但尚不能指定各前缘点处的钝化半径，故在设计飞行器气动外形时无法根据热防护需求所提出的前缘约束对钝化前缘进行针对性设计。
发明内容
针对现有技术存在的上述缺陷，本发明目的在于提供一种半径可控的参数化三维前缘钝化设计方法，其能够解决现有的前缘钝化方法不能兼顾钝化前缘与飞行器上下表面光滑连接与指定各前缘点钝化半径的问题。本方法使钝化前缘能在满足与飞行器上下表面光滑连接且钝化前缘形状沿横向连续变化的同时，还可根据热防护需求来设计各前缘点处不同的钝化半径，有效地实现了热防护与气动外形的多学科耦合设计。
为实现本发明之目的，采用以下技术方案予以实现：
一种半径可控的参数化三维前缘钝化设计方法，首先给定初始乘波体外形，然后对初始乘波体外形中的尖锐前缘进行钝化。给定的初始乘波体外形，既可是吻切锥乘波体，但不局限于该乘波体或锥导乘波体等其他类型的乘波体外形，还可用于其他带有尖锐前缘的飞行器外形。也就是本发明接下来重点介绍的钝化方法可以应用于带有尖锐前缘的乘波体或者是带有尖锐前缘的飞行器，实现乘波体或飞行器的三维前缘钝化。
具体地，对尖锐前缘进行钝化，包括以下步骤：
S1将乘波体上表面抬升一定高度；
给定初始乘波体外形，即已知乘波体未钝化时的原始前缘点为P0和未钝化前的乘波体纵向截面曲线；未钝化前的乘波体纵向截面曲线由前缘点、上表面曲线和下表面曲线连接而成。
设对乘波体的上表面沿法向抬升高度h，抬升后的乘波体上表面端点为P1，抬升后的乘波体下表面端点为乘波体未钝化前的原始前缘点P0。抬升后的乘波体纵向截面曲线由上表面曲线、下表面曲线以及钝化前缘曲线连接而成，其中上表面曲线和下表面曲线即是未钝化前的乘波体纵向截面曲线上的上表面曲线和下表面曲线，因此是已知的，钝化前缘曲线是待求解的；抬升后的乘波体上表面端点P1是上表面曲线上的点，抬升后的乘波体下表面端点P0是下表面曲线上的点。
S2获取抬升后的乘波体前缘点P2的插值条件；
抬升后的乘波体前缘点P2的插值条件为抬升后的乘波体下表面端点P0处的切矢v0、曲率k0和抬升后的乘波体上表面端点P1处切矢v1和曲率k1。
S3给出抬升后的钝化前缘曲线上不同横向站位处的钝化半径。
S4设置抬升后的乘波体前缘点P2以及乘波体的上表面沿法向抬升高度h。
S5求出所有钝化前缘曲线，然后放样得到钝化前缘曲面。
本发明S2中，对上表面曲线进行均匀离散，获得多个离散点，获得的多个离散点中包括抬升后的乘波体上表面端点P1，对多个离散点进行差分，能够获得抬升后的乘波体上表面端点P1的1阶和2阶导矢，进而获得抬升后的乘波体上表面端点P1处的切矢v1和曲率k1；
同样的，对下表面曲线进行均匀离散，获得多个离散点，获得的多个离散点中包括抬升后的乘波体上表面端点P1，对多个离散点进行差分，能够获得抬升后的乘波体下表面端点P0的1阶和2阶导矢，进而获得抬升后的乘波体下表面端点P0处的切矢v0和曲率k0。
本发明S3中，未钝化前的乘波体纵向截面曲线中的前缘曲线已知，故前缘曲线上各横向站位处的前缘后掠角已知，前缘曲线上各横向站位均一一对应钝化前缘曲线上的横向站位，钝化前缘曲线上各横向站位处的当地前缘后掠角与前缘曲线上与之对应的横向站位处的当地前缘后掠角看作相同，即钝化前缘曲线上各横向站位处的当地前缘后掠角λ已知。
设定抬升后的乘波体前缘点P2处的前缘钝化半径，由预先设定的抬升后的乘波体前缘点P2处的前缘钝化半径和钝化前缘曲线上各横向站位处的当地前缘后掠角λ来确定抬升后的钝化前缘曲线上各横向站位处的前缘钝化半径，如式(3)所示：
rswept＝rsymmetry(cosλ)2.2 (3)
其中，rswept为钝化前缘曲线上各横向站位处的前缘钝化半径；rsymmetry为预先设定的抬升后的乘波体前缘点P2处的前缘钝化半径。
本发明S4中，抬升后的乘波体前缘点P2和乘波体的上表面沿法向抬升高度h的设置应满足以下两个要求：
1)抬升后的乘波体前缘点P2和乘波体的上表面沿法向抬升高度h的设置能使乘波体钝化前缘曲线能够成功生成；
2)抬升后的乘波体前缘点P2的曲率大于乘波体钝化前缘曲线上各横向站位处的曲率。
进一步地，本发明将乘波体的上表面沿法向抬升高度h设置为当地前缘钝化半径的1～2倍。
抬升后的乘波体前缘点P2设置在上表面曲线以及下表面曲线的延长线所构成夹角的角平分线上，抬升后的乘波体前缘点P2距P0点的水平距离l可由式(4)得出：
l＝rswept/cosλ (4)
其中，rswept为当地前缘钝化半径。
本发明S5中，根据点P0和P1以及S2获取的抬升后的乘波体前缘点P2的插值条件求出插值曲线即钝化前缘曲线；将求解钝化前缘曲线的问题看作插值于两个端点的GHI问题，利用3次Bézier曲线构造钝化前缘曲线，求出所有的钝化前缘曲线，选取其中同时满足S4中要求的钝化前缘曲线进行放样，从而得到所需的钝化前缘曲面。具体地，在本发明中是以P0、P2为端点利用一段3次Bézier曲线构造点P0和P2间的钝化前缘曲线，以P2、P1为端点利用一段3次Bézier曲线构造点P2和P1间的钝化前缘曲线，完整的钝化前缘曲线即由点P0和P2间的钝化前缘曲线和点P2和P1间的钝化前缘曲线连接而成，求出所有的钝化前缘曲线，选取其中同时满足S4中要求的钝化前缘曲线进行放样，从而得到所需的钝化前缘曲面。
3次Bézier曲线共需要4个控制点来进行曲线的求解，因此首先根据已知的两个控制点(即点P0和P2或者点P2和P1)求得插值曲线中的另外两个控制点，然后根据已知的两个控制点以及求解得到的另外两个控制点，即4个控制点来进行插值曲线的求解。
下面介绍如何以P0、P2为端点利用一段3次Bézier曲线构造点P0和P2间的钝化前缘曲线，方法如下：
已知的两个控制点P0和P2(即两个端点)与待求的两个控制点P3和P4，这4个控制点间的关系可由下式表示
其中：δ0和δ1为设定的中间参数，为待求值；v0和v2分别为端点P0和P2处的切矢。
两个端点P0和P2处的曲率为
由S2获取的抬升后的乘波体前缘点P1的插值条件可以得到关于δ0和δ1的方程组
其中d是两个端点P0和P2之间的位移矢量。
通过求解方程组(7)可以得到未知参数δ0和δ1的值，从而计算得到其余两个控制点P3和P4的坐标，进而得到一条曲率连续的插值曲线，即利用3次Bézier曲线构造的点P0和P2间的初步钝化前缘曲线。
方程组(7)的解可能有0至3个解。由于乘波体钝化前缘曲线为三维空间曲线，钝化前缘曲线上各个横向站位处的端点插值条件不一，可能导致方程无解或者多解。如要将该方法应用于三维外形上，则对于方程组解的讨论将是主要障碍之一。此处对插值条件进行简化，即不要求钝化前缘曲线与上表面曲线、下表面曲线间的曲率连续，将钝化前缘曲线与上表面曲线、下表面曲线连接处的曲率值设置为0。
方程组(7)可简化为
原方程组(7)变为一个一元一次方程与一个二元二次方程的组合，不难看出，方程组(8)必然有且只有一个解。通过对抬升后的乘波体前缘点P2的插值条件的调整，以降低钝化前缘曲线与上表面曲线、下表面曲线连接的光滑性为代价，简化了Bézier曲线的求解过程，有利于保持不同横向站位处钝化前缘曲线之间的连续性。
以P2、P1为端点利用一段3次Bézier曲线构造点P2和P1间的钝化前缘曲线的方法与以P0、P2为端点利用一段3次Bézier曲线构造点P0和P2间的钝化前缘曲线的方法完全相同，在此不再赘述。
本发明的有益技术效果是：
本发明基于增加材料的前缘钝化思路和基于Bézier曲线和有理Bézier曲线的前缘钝化方法的半径可控的参数化三维前缘钝化设计方法，在保证钝化前缘曲线与上、下表面曲线光滑连接的同时，还可以根据热防护需求指定各前缘点处的钝化半径，有效地实现了热防护与气动外形的多学科耦合设计。
综上所述，该方法主要有以下几个特点：
1)乘波体前缘钝化过程实现了参数化，有利于程序实现和优化设计，减少了飞行器设计的工作量；
2)钝化前缘曲面与乘波体表面光滑连接，可以避免飞行器表面不光滑对性能所带来的不利影响；
3)前缘点以及前缘点处的钝化半径可以直接指定，使得气动外形设计与热防护系统设计可以实现有效衔接。
附图说明
图1是实施例给定的乘波体外形形体示意图；
图2是移除材料和增加材料两种前缘钝化思路对比示意图，上部分为乘波体尖锐前缘的侧视图和后视图，中间部分为移除材料钝化前缘示意图，下面部分是增加材料钝化前缘示意图；
图3是乘波体上表面“抬升”前后的乘波体纵向截面示意图；
图4是乘波体上表面抬升后的新前缘点P2设置示意图，其中l为钝化前后前缘点的水平距离；
图5是三次Bézier曲线示意图，其中：P0、P1、P2、P3为设计插值曲线的四个控制点，v0、k0分别是P0点处的切矢和曲率，相应地，v2、k2分别是P2点处的切矢和曲率，d为P0和P2之间的的位移矢量；
图6是利用2段3次Bézier曲线构造钝化前缘曲线示意图，P0、P2和P2、P1分别是2段3次Bézier曲线的端点，v、k则分别是各点处的切矢和曲率；
图7是脊形吻切锥乘波体等轴测视图及头部钝化细节图；
图8是有粘和无粘升力系数随攻角变化曲线；
图9是有粘和无粘阻力系数随攻角变化曲线；
图10是有粘和无粘俯仰力矩系数随攻角变化曲线；
图11是有粘和无粘升阻比数随攻角变化曲线。
具体实施方式
下面将结合本发明实施例图中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，做进一步详细说明，但本发明的实施方式不仅限于此。
一种半径可控的参数化三维前缘钝化设计方法，其首先给定乘波体外形。本实施例以脊形吻切锥乘波体为例进行前缘钝化设计，参照图1，为本实施例给定的乘波体外形形体示意图。特别强调，本发明所提出的钝化方法，既可用于吻切锥乘波体，但不局限于该乘波体或锥导乘波体等其他类型的乘波体，还可用于其他带有尖锐前缘的飞行器外形，本实施例仅以吻切锥乘波体为例进行说明。
现有的各种乘波体前缘钝化方法在材料的增减方面存在差异：Steven提出了通过移除材料和增加材料两种途径进行前缘钝化的思路，参见Stevens D R.Practicalconsiderations in waverider applications[M].AIAAAircraft Design SystemsMeeting.1992.。如图2所示，其中移除材料的方法被一些学者在研究中所采用。考虑到乘波体的前缘较薄，移除材料的方法可能对升力面的几何形状会产生较大影响，进而使得升阻比和容积减小。然而实际上采用增加材料的方法对乘波体容积、升力面和流场结构的影响较小；本发明采用增加材料的钝化思路，沿流向设计各个横向站位上钝化前缘曲线。
具体的，对给定的乘波体外形中的尖锐前缘进行钝化，包括以下步骤：
S1将乘波体上表面抬升一定高度
给定初始乘波体外形，即已知乘波体未钝化时的原始前缘点为P0和未钝化前的乘波体纵向截面曲线；未钝化前的乘波体纵向截面曲线由前缘点、上表面曲线和下表面曲线连接而成。
设对乘波体的上表面沿法向抬升高度h，抬升后的乘波体上表面端点为P1，抬升后的乘波体下表面端点为乘波体未钝化前的原始前缘点P0。抬升后的乘波体纵向截面曲线由上表面曲线、下表面曲线以及钝化前缘曲线连接而成，其中上表面曲线和下表面曲线即是未钝化前的乘波体纵向截面曲线上的上表面曲线和下表面曲线，因此是已知的，钝化前缘曲线是待求解的；抬升后的乘波体上表面端点P1是上表面曲线上的点，抬升后的乘波体下表面端点P0是下表面曲线上的点。参照图3是乘波体上表面“抬升”前后的乘波体纵向截面示意图。
S2获取抬升后的乘波体前缘点P2的插值条件
抬升后的乘波体前缘点P2的插值条件为抬升后的乘波体下表面端点P0处的切矢v0、曲率k0和抬升后的乘波体上表面端点P1处切矢v1和曲率k1。
对上表面曲线进行均匀离散，获得多个离散点，获得的多个离散点中包括抬升后的乘波体上表面端点P1，对多个离散点进行差分，能够获得抬升后的乘波体上表面端点P1的1阶和2阶导矢，进而获得抬升后的乘波体上表面端点P1处的切矢v1和曲率k1；
同样的，对下表面曲线进行均匀离散，获得多个离散点，获得的多个离散点中包括抬升后的乘波体上表面端点P1，对多个离散点进行差分，能够获得抬升后的乘波体下表面端点P0的1阶和2阶导矢，进而获得抬升后的乘波体下表面端点P0处的切矢v0、曲率k0；
其中：具体的差分方法参见杨炯,宁涛,席平.前缘点曲率可控的曲率连续前缘几何设计[J].计算机辅助设计与图形学学报,2016,28(7):1195-200.以抬升后的乘波体上表面端点P1为例，具体的差分方法为：由于上表面曲线已知，所以可以将上表面曲线离散为多个点P1、P2、P3、P4、P5…Pn，Pi＝(xi，yi)，根据导矢的定义，点Pi对应的一阶导矢可由点位置的差分计算得到:
一阶差分：δf(xi+1)＝f(xi+2)-f(xi)
一阶导矢：ui+1＝(2h，yi+2-yi)，其中h为相邻离散点间的x轴向距离
二阶差分：δ2f(xi+1)＝f(xi+2)+f(xi)-2f(xi+1)
二阶导矢：u′(xi+1)＝(2h，yi+2+yi-2yi+1)
即可得到插值条件。
S3给出抬升后的钝化前缘曲线上不同横向站位处的钝化半径；
抬升后的乘波体前缘点P2(乘波体纵向对称面处的驻点)是热防护压力最大的部位，该处的前缘钝化半径应由热防护材料和飞行工况等条件确定，是预先设定好的。
未钝化前的乘波体纵向截面曲线中的前缘曲线已知，故前缘曲线上各横向站位处的前缘后掠角已知，前缘曲线上各横向站位均一一对应钝化前缘曲线上的横向站位，钝化前缘曲线上各横向站位处的当地前缘后掠角与前缘曲线上与之对应的横向站位处的当地前缘后掠角相同，即钝化前缘曲线上各横向站位处的当地前缘后掠角λ已知。
由预先设定的抬升后的乘波体前缘点P2处的前缘钝化半径和钝化前缘曲线上各横向站位处的当地前缘后掠角λ来确定钝化前缘曲线上各横向站位处(即不同横向坐标下前缘曲线上的点)的前缘钝化半径，如式(3)所示。
rswept＝rsymmetry(cosλ)2.2 (3)
其中，rswept为钝化前缘曲线上各横向站位处的前缘钝化半径(由公式(3)得到)，rsymmetry为抬升后的乘波体前缘点P2处的前缘钝化半径。
S4设置抬升后的乘波体前缘点P2以及乘波体的上表面沿法向抬升高度h；
抬升后的乘波体前缘点P2和乘波体的上表面沿法向抬升高度h的设置应满足以下两个要求：
1)参数设置能使钝化前缘曲线能够成功生成；2)抬升后的乘波体前缘点P2的曲率大于乘波体钝化前缘曲线上各横向站位处的曲率。
进一步地，本发明优选将乘波体的上表面沿法向抬升高度h设置为当地前缘钝化半径的1～2倍即可。
将抬升后的乘波体前缘点P2设置在上表面曲线以及下表面曲线的延长线所构成夹角的角平分线上，如图4所示，其距P0点的水平距离l可由式(4)得出，
l＝rswept/cosλ (4)
S5求出所有钝化前缘曲线，然后放样得到钝化前缘曲面；
根据点P0和P1以及S2获取的抬升后的乘波体前缘点P2的插值条件求出插值曲线即钝化前缘曲线，将求解钝化前缘曲线的问题看作插值于两个端点的GHI问题，以P0、P2为端点利用一段3次Bézier曲线构造点P0和P2间的钝化前缘曲线，以P2、P1为端点利用一段3次Bézier曲线构造点P2和P1间的钝化前缘曲线，完整的钝化前缘曲线即由点P0和P2间的钝化前缘曲线和点P2和P1间的钝化前缘曲线连接而成，求出所有的钝化前缘曲线，选取其中同时满足S4中要求的钝化前缘曲线进行放样，从而得到所需的钝化前缘曲面。
根据端点条件构造插值曲线的问题，可以看作插值于两个端点的GHI(GeometricHermite Interpolation，几何厄米特插值)问题。该问题由de Boor首次提出，参见Boor CD,K,Sabin M.High accuracy geometric Hermite interpolation[J].ComputerAided Geometric Design,1987,4(4):269-78.，其利用3次Bézier曲线插值于给定条件的平面上的两个点，从而解决了GHI问题。应用于GHI问题中的曲线类型有多种，其中Bézier曲线有着良好的凸包性和较简单的计算方法，在实际应用中有一定的优势。
3次Bézier曲线共需要4个控制点来进行曲线的求解，因此首先根据已知的两个控制点(即点P0和P2或者点P2和P1)求得插值曲线中的另外两个控制点，然后根据已知的两个控制点以及求解得到的另外两个控制点，即4个控制点来进行插值曲线的求解,如图5所示。
下面介绍如何以P0、P2为端点利用一段3次Bézier曲线构造点P0和P2间的钝化前缘曲线，方法如下：
已知的两个控制点P0和P2(即两个端点)与待求的两个控制点P3和P4，这4个控制点间的关系可由下式表示
其中：δ0和δ1为设定的中间参数，为待求值；v0和v2分别为端点P0和P2处的切矢。
两个端点P0和P2处的曲率为
由S2获取的抬升后的乘波体前缘点P2的插值条件可以得到关于δ0和δ1的方程组
其中d是两个端点P0和P2之间的位移矢量。
通过求解方程组(7)可以得到未知参数δ0和δ1的值，从而计算得到其余两个控制点P3和P4的坐标，进而得到一条曲率连续的插值曲线，即利用3次Bézier曲线构造的点P0和P2间的初步钝化前缘曲线。
方程组(7)的解可能有0至3个解。由于乘波体钝化前缘曲线为三维空间曲线，钝化前缘曲线上各个横向站位处的端点插值条件不一，可能导致方程无解或者多解。如要将该方法应用于三维外形上，则对于方程组解的讨论将是主要障碍之一。此处对插值条件进行简化，即不要求钝化前缘曲线与上表面曲线、下表面曲线间的曲率连续，将钝化前缘曲线与上表面曲线、下表面曲线连接处的曲率值设置为0。
方程组(7)可简化为
原方程组(7)变为一个一元一次方程与一个二元二次方程的组合，不难看出，方程组(8)必然有且只有一个解。通过对抬升后的乘波体前缘点P2的插值条件的调整，以降低钝化前缘曲线与上表面曲线、下表面曲线连接的光滑性为代价，简化了Bézier曲线的求解过程，有利于保持不同横向站位处钝化前缘曲线之间的连续性。
以P2、P1为端点利用一段3次Bézier曲线构造点P2和P1间的钝化前缘曲线的方法与以P0、P2为端点利用一段3次Bézier曲线构造点P0和P2间的钝化前缘曲线的方法完全相同，在此不再赘述。
本实施例中，利用2段3次Bézier曲线构造钝化前缘曲线，具体如图6所示。图6是利用2段3次Bézier曲线构造钝化前缘曲线示意图，P0、P2和P2、P1分别是2段3次Bézier曲线的端点，v、k则分别是各点处的切矢和曲率。钝化前缘曲线由两段曲线连接而成，即点P0和P2间的钝化前缘曲线和点P2和P1间的钝化前缘曲线。同样的将求解钝化前缘曲线的问题看作插值于两个端点的GHI问题，以P0、P2为端点利用一段3次Bézier曲线构造点P0和P2间的钝化前缘曲线，以P2、P1为端点利用一段3次Bézier曲线构造点P2和P1间的钝化前缘曲线，求出所有的钝化前缘曲线。实际生成前缘钝化曲面时，较难保证所有生成的钝化前缘曲线均能满足前述两点要求，因此选取其中同时满足S4中所有要求的钝化前缘曲线进行放样，从而得到所需的钝化前缘曲面。本实施例最终得到钝化后的吻切锥乘波体外形如图7所示。
对钝化后的乘波体体进行气动性能分析，由于乘波体为面对称外形，在只研究俯仰特性时可取一半流场进行数值模型，计算域包括了乘波体的底部。计算网格为非结构网格，为提高对边界层的预测精度，在壁面附近划分了棱柱网格以提高边界层的预测精度。流动控制方程为可压流的三维雷诺平均Navier-Stokes方程(RANS)，在计算中采用二阶迎风格式进行空间离散，时间项为隐格式。采用RNG-两方程湍流模型和增强壁面函数模拟湍流，壁面边界条件设置为绝热无滑移壁面。在马赫数Ma＝8.0和高度H＝31km时，气动性能随攻角变化规律如图8～11所示。
不难看出，升力系数和俯仰力矩系数与攻角基本上呈线性关系，且值随着攻角的增大而增大。粘性对升力系数和俯仰力矩系数的影响不大，在飞行器设计初期可通过无粘计算快速获取升力特性和俯仰特性。由于高超声速流动具有强粘性的特点，乘波体阻力系数的有粘和无粘计算结果存在较大的差异。乘波体的最大升阻比出现在8°攻角，与一般研究认为的最大升阻比出现在小攻角的结论不一致。可能的原因是本文采用了脊形的设计，上表面与一般乘波体所用的自由流面设计有较大区别，进而导致了升阻比和俯仰特性存在一定的差异。
以上包含了本发明优选实施例的说明，这是为了详细说明本发明的技术特征，并不是想要将发明内容限制在实施例所描述的具体形式中，依据本发明内容主旨进行的其他修改和变型也受本专利保护。本发明内容的主旨是由权利要求书所界定，而非由实施例的具体描述所界定。

Claims (5)

1.一种半径可控的参数化三维前缘钝化设计方法，其特征在于：首先给定初始乘波体外形，然后对初始乘波体外形中的尖锐前缘进行钝化，包括以下步骤：
S1将乘波体上表面抬升一定高度；
给定初始乘波体外形，即已知乘波体未钝化时的原始前缘点为P0和未钝化前的乘波体纵向截面曲线；未钝化前的乘波体纵向截面曲线由前缘点、上表面曲线和下表面曲线连接而成；
设对乘波体的上表面沿法向抬升高度h，抬升后的乘波体上表面端点为P1，抬升后的乘波体下表面端点为乘波体未钝化前的原始前缘点P0；抬升后的乘波体纵向截面曲线由上表面曲线、下表面曲线以及钝化前缘曲线连接而成，其中上表面曲线和下表面曲线即是未钝化前的乘波体纵向截面曲线上的上表面曲线和下表面曲线，因此是已知的，钝化前缘曲线是待求解的；抬升后的乘波体上表面端点P1是上表面曲线上的点，抬升后的乘波体下表面端点P0是下表面曲线上的点；
S2获取抬升后的乘波体前缘点P2的插值条件；
抬升后的乘波体前缘点P2的插值条件为抬升后的乘波体下表面端点P0处的切矢v0、曲率k0和抬升后的乘波体上表面端点P1处切矢v1和曲率k1；
S3给出抬升后的钝化前缘曲线上不同横向站位处的钝化半径；
未钝化前的乘波体前缘曲线已知，故前缘曲线上各横向站位处的前缘后掠角已知，前缘曲线上各横向站位均一一对应钝化前缘曲线上的横向站位，钝化前缘曲线上各横向站位处的当地前缘后掠角与前缘曲线上与之对应的横向站位处的当地前缘后掠角看作相同，即钝化前缘曲线上各横向站位处的当地前缘后掠角λ已知；
设定抬升后的乘波体前缘点P2处的前缘钝化半径，由预先设定的抬升后的乘波体前缘点P2处的前缘钝化半径和钝化前缘曲线上各横向站位处的当地前缘后掠角λ来确定抬升后的钝化前缘曲线上各横向站位处的前缘钝化半径，如式(3)所示：
rswept＝rsymmetry(cosλ)2.2 (3)
其中，rswept为钝化前缘曲线上各横向站位处的前缘钝化半径；rsymmetry为预先设定的抬升后的乘波体前缘点P2处的前缘钝化半径；
S4设置抬升后的乘波体前缘点P2以及乘波体的上表面沿法向抬升高度h应满足以下两个要求：
1)抬升后的乘波体前缘点P2和乘波体的上表面沿法向抬升高度h的设置能使乘波体钝化前缘曲线能够成功生成；
2)抬升后的乘波体前缘点P2的曲率大于乘波体钝化前缘曲线上各横向站位处的曲率；
S5求出所有钝化前缘曲线，然后放样得到钝化前缘曲面；
根据点P0和P1以及S2获取的抬升后的乘波体前缘点P2的插值条件求出插值曲线即钝化前缘曲线；将求解钝化前缘曲线的问题看作插值于两个端点的GHI问题，利用3次Bézier曲线构造钝化前缘曲线，求出所有的钝化前缘曲线，选取其中同时满足S4中要求的钝化前缘曲线进行放样，从而得到所需的钝化前缘曲面，方法如下：
以P0、P2为端点利用一段3次Bézier曲线构造点P0和P2间的纯化前缘曲线]首先将点P0和P2作为已知的两个控制点求得插值曲线中的另外两个控制点，然后根据已知的两个控制点以及求解得到的另外两个控制点即4个控制点来进行插值曲线的求解，求得的插值曲线即P0和P2间的钝化前缘曲线；
以P2、P1为端点利用一段3次Bézier曲线构造点P2和P1间的钝化前缘曲线，首先将点P2和P1作为已知的两个控制点求得插值曲线中的另外两个控制点，然后根据已知的两个控制点以及求解得到的另外两个控制点即4个控制点来进行插值曲线的求解，求得的插值曲线即P2和P1间的钝化前缘曲线；
完整的钝化前缘曲线即由点P0和P2间的钝化前缘曲线和点P2和P1间的钝化前缘曲线连接而成，求出所有的钝化前缘曲线，选取其中同时满足S4中要求的钝化前缘曲线进行放样，从而得到所需的钝化前缘曲面。
2.根据权利要求1所述的半径可控的参数化三维前缘钝化设计方法，其特征在于：S2中，对上表面曲线进行均匀离散，获得多个离散点，获得的多个离散点中包括抬升后的乘波体上表面端点P1，对多个离散点进行差分，能够获得抬升后的乘波体上表面端点P1的1阶和2阶导矢，进而获得抬升后的乘波体上表面端点P1处的切矢v1和曲率k1；
同样的，对下表面曲线进行均匀离散，获得多个离散点，获得的多个离散点中包括抬升后的乘波体上表面端点P1，对多个离散点进行差分，能够获得抬升后的乘波体下表面端点P0的1阶和2阶导矢，进而获得抬升后的乘波体下表面端点P0处的切矢v0和曲率k0。
3.根据权利要求2所述的半径可控的参数化三维前缘钝化设计方法，其特征在于：S4中，乘波体的上表面沿法向抬升高度h设置为当地前缘钝化半径的1～2倍；
抬升后的乘波体前缘点P2设置在上表面曲线以及下表面曲线的延长线所构成夹角的角平分线上，抬升后的乘波体前缘点P2距P0点的水平距离l可由式(4)得出：
l＝rswept/cosλ (4)
其中，rswept为当地前缘钝化半径。
4.根据权利要求1所述的半径可控的参数化三维前缘钝化设计方法，其特征在于：S5中，以P0、P2为端点利用一段3次Bézier曲线构造点P0和P2间的钝化前缘曲线，方法如下：
已知的两个控制点P0和P2与待求的两个控制点P3和P4，这4个控制点间的关系可由下式表示
其中：δ0和δ1为设定的中间参数，为待求值；v0和v2分别为端点P0和P2处的切矢；
点P0和P2处的曲率为
由S2获取的抬升后的乘波体前缘点P2的插值条件可以得到关于δ0和δ1的方程组
其中d是点P0和P2之间的位移矢量；
通过求解方程组(7)可以得到未知参数δ0和δ1的值，从而计算得到其余两个控制点P3和P4的坐标，进而得到一条曲率连续的插值曲线，即利用3次Bézier曲线构造的点P0和P2间的初步钝化前缘曲线；
将钝化前缘曲线与上表面曲线、下表面曲线连接处的曲率值设置为0；方程组(7)简化为
通过求解方程组(8)得到最终的P2和P1间的钝化前缘曲线；
按照与以P2、P1为端点利用一段3次Bézier曲线构造点P2和P1间的钝化前缘曲线相同的方法构造点P0和P2间的钝化前缘曲线。
5.根据权利要求1所述的半径可控的参数化三维前缘钝化设计方法，其特征在于：给定的初始乘波体外形是带有尖锐前缘的乘波体外形或者是带有尖锐前缘的飞行器外形。
CN201711102724.7A 2017-11-10 2017-11-10 一种半径可控的参数化三维前缘钝化设计方法 Active CN107766673B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711102724.7A CN107766673B (zh) 2017-11-10 2017-11-10 一种半径可控的参数化三维前缘钝化设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711102724.7A CN107766673B (zh) 2017-11-10 2017-11-10 一种半径可控的参数化三维前缘钝化设计方法

Publications (2)

Publication Number Publication Date
CN107766673A CN107766673A (zh) 2018-03-06
CN107766673B true CN107766673B (zh) 2018-07-31

Family

ID=61273317

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711102724.7A Active CN107766673B (zh) 2017-11-10 2017-11-10 一种半径可控的参数化三维前缘钝化设计方法

Country Status (1)

Country Link
CN (1) CN107766673B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109443690B (zh) * 2018-12-10 2020-10-09 中国科学院力学研究所 一种高焓风洞喷管型线的光滑化定型方法
CN111173619B (zh) * 2019-12-20 2024-09-20 中国科学技术大学 内转式进气道及其v字形钝前缘

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7431242B1 (en) * 2005-11-01 2008-10-07 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration (Nasa) Re-entry vehicle shape for enhanced performance
CN103049597A (zh) * 2012-11-27 2013-04-17 中国人民解放军国防科学技术大学 变尺度边缘钝化乘波构型及其钝化方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2952454B1 (fr) * 2009-11-06 2011-12-23 Airbus Operations Sas Procede et outil de simulation du comportement aerodynamique d'un element aerodynamique d'un aeronef presentant un angle de fleche variable

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7431242B1 (en) * 2005-11-01 2008-10-07 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration (Nasa) Re-entry vehicle shape for enhanced performance
CN103049597A (zh) * 2012-11-27 2013-04-17 中国人民解放军国防科学技术大学 变尺度边缘钝化乘波构型及其钝化方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Effectiveness of side force models for flow simulations downstream of vortex generators;Liesbeth Florentie等;《AIAA Journal》;20170430;第55卷(第4期);第1373-1384页 *
Performance and regression rate characteristics of 5-kN swirling-oxidizer-flow-type hybrid rocket engine;Takashi Sakurai等;《Journal of propulsion and power》;20170830;第33卷(第4期);第891-901页 *
乘波构型的钝化方法及其对性能影响研究;刘济民等;《宇航学报》;20110530;第32卷(第5期);第966-974页 *
乘波构型钝化方法分析及性能研究;陈小庆等;《空气动力学学报》;20170228;第35卷(第1期);第108-114页 *
前缘点曲率可控的曲率连续前缘几何设计;杨炯等;《计算机辅助设计与图形学学报》;20160730;第28卷(第7期);第1196-1200页 *

Also Published As

Publication number Publication date
CN107766673A (zh) 2018-03-06

Similar Documents

Publication Publication Date Title
CN107963236B (zh) 基于密切锥理论的定平面乘波体设计方法
Lanfrit Best practice guidelines for handling Automotive External Aerodynamics with FLUENT
Eastwood et al. Toward designing with three-dimensional bumps for lift/drag improvement and buffet alleviation
CN105913483A (zh) 一种三维交叉道路模型生成的方法及装置
CN107766673B (zh) 一种半径可控的参数化三维前缘钝化设计方法
KR102616901B1 (ko) 넓은 속도 범위 극초음속기의 공기역학적 구성 설계 방법 및 시스템
CN104063746B (zh) 一种基于旅行商问题求解器的曲面精加工刀路生成方法
CN110414168A (zh) 基于与前机身耦合优化的高超声速隔离段设计方法及系统
WO2014071884A1 (zh) 单曲风挡机头一体化设计方法
CN110210185B (zh) 一种高超声速隔离段优化设计方法及系统
CN113536458B (zh) 一种基于类别形状函数变换的可变形机翼的快速几何建模方法
Marcum et al. Unstructured mesh generation using advancing layers and metric-based transition for viscous flowfields
CN103676786B (zh) 一种基于加速度原理的曲线光顺方法
CN103413019B (zh) 一种不规则轮廓双曲率外形面整体壁板离散方法
CN110414016A (zh) 超高速管道运输工具的乘波体外形参数化设计方法及系统
CN102930586A (zh) 一种基于线性旋转不变微分坐标的可交互几何变形方法
CN108804791A (zh) 一种适用于埋入式进气道布局的飞行器参数化方法
Wissink et al. Application of 3D Strand Solver to Rotorcraft Hover
Ito et al. Solution adaptive mesh generation using feature-aligned embedded surface meshes
Shaw et al. Study of dimple effect on aerodynamic drag characteristics of a car
Oxyzoglou et al. Effects of heaving motion on the aerodynamic performance of a double element wing in ground effect
CN112606996A (zh) 一种乘波体结构
Kuzmin et al. Airfoils admitting anomalous behavior of lift coefficient in descending transonic flight
Kashyap et al. Computational analysis of flap Camber and ground clearance in double-element inverted airfoils
Kuwahara et al. Direct simulation of a flow around a subsonic airfoil

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