CN112948976A - 一种平面前缘升力体前缘线优化方法、系统及升力体 - Google Patents

一种平面前缘升力体前缘线优化方法、系统及升力体 Download PDF

Info

Publication number
CN112948976A
CN112948976A CN202110274578.6A CN202110274578A CN112948976A CN 112948976 A CN112948976 A CN 112948976A CN 202110274578 A CN202110274578 A CN 202110274578A CN 112948976 A CN112948976 A CN 112948976A
Authority
CN
China
Prior art keywords
leading edge
original
quasi
edge line
waverider
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
CN202110274578.6A
Other languages
English (en)
Other versions
CN112948976B (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.)
Institute of Mechanics of CAS
Original Assignee
Institute of Mechanics of CAS
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 Institute of Mechanics of CAS filed Critical Institute of Mechanics of CAS
Priority to CN202110274578.6A priority Critical patent/CN112948976B/zh
Publication of CN112948976A publication Critical patent/CN112948976A/zh
Application granted granted Critical
Publication of CN112948976B publication Critical patent/CN112948976B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B64—AIRCRAFT; AVIATION; COSMONAUTICS
    • B64C—AEROPLANES; HELICOPTERS
    • B64C30/00—Supersonic type aircraft
    • 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
    • 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
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Computational Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Physics (AREA)
  • Manufacturing & Machinery (AREA)
  • Transportation (AREA)
  • Eyeglasses (AREA)

Abstract

本发明公开了一种平面前缘升力体前缘线优化方法、系统及升力体,所述方法包括在计算坐标系下获取原始准乘波体构型的原始前缘线的首部端点以及多个尾部端点；所述计算坐标系为原始准乘波体构型的本体坐标系，所述原始准乘波体构型为原始尖头前缘准乘波体构型，所述原始前缘线为具有拐点且不可导的前缘线；将所述首部端点与其中一个所述尾部端点连接形成直线段；经过位于所述直线段上的所述尾部端点形成与Y轴平行的直线；本申请实施例提供的平面前缘升力体前缘线优化方法，可以将原始尖头准乘波体构型的前缘线在指定平面上投影生成新的前缘线，可以克服原始尖头准乘波体构型外形畸变、不利于工程加工制造的问题。

Description

一种平面前缘升力体前缘线优化方法、系统及升力体
技术领域
本发明涉及乘波构型优化技术领域，特别是涉及一种高超声速的平面前缘升力体前缘线优化方法、系统及升力体。
背景技术
在高超声速条件下，由于常规外形的飞行器在超音速流中前缘大都是脱体激波，激波前后存在的压差使得常规外形的飞行器上的波阻非常大，飞行器将遭受极大的摩阻和波阻，升阻比提升十分困难，会面临难以逾越的“升阻比屏障”。为了解决上述问题，乘波构型的飞行器应运而生。乘波构型也称为乘波体，是一种适宜高超音速飞行的飞行器外形，其所有的前缘都具有附体激波。乘波构型飞行时其前缘平面与激波的上表面重合，就像骑在激波的波面上，依靠激波的压力产生升力。由于乘波体的上表面与自由流面同面，所以不会形成大的压差阻力。因此，乘波体被认为是最有希望打破高超声速“升阻比屏障”的一种新型气动布局。
基于各种基准流场的乘波体生成方法的发展,使得在飞行器设计过程中可以根据不同需求选择不同类型的乘波体构型。但是,乘波体构型一般难以直接满足容积率、配平、稳定性等基本的飞行器设计工程需求。针对这一情况,“准乘波体”构型优化设计方法应运而生，“准乘波体”设计过程为：保留乘波体的前缘线；从前缘线上的点出发，沿每个纵向截面采用可解析表达的曲线确定该截面型线，直到底面截止；保持对称面处型线高度与原乘波体一致；上表面仍然沿用自由流面。最后，再采用遗传算法对曲线进行优化，即可生成满足配平和纵向稳定性的高升阻比准乘波体构型。
根据前缘形状，可将准乘波体分为圆头前缘准乘波体和尖头前缘准乘波体两种形式，其典型外形分别如图2(a-d)和图3(a-d)所示。从图中可见，圆头前缘准乘波体外形比较规则，主要特征是在对称面处前缘线连续可导；而尖头前缘准乘波体外形不规则，主要特征是在对称面处前缘线有拐点、不可导，从图3(a-d)中的侧视图可见，该外形的前缘线高低不平，使得整个外形畸变比较严重。因此，在工程应用中，尖头前缘准乘波体外形不利于加工制造。然而，相比圆头乘波体，尖头乘波体在靠近头部处前缘线后掠角大，有利于降低热流，这一点对高超声速飞行器又尤为关键。
因此，如何克服尖头准乘波体外形畸变严重、不利于工程加工制造的难题，是迫切需要本领域技术人员解决的技术问题。
发明内容
本发明提供了一种平面前缘升力体前缘线优化方法、系统及升力体。
本发明提供了如下方案：
一种平面前缘升力体前缘线优化方法，包括：
所述方法包括：
在计算坐标系下获取原始准乘波体构型的原始前缘线的首部端点以及多个尾部端点；所述计算坐标系为原始准乘波体构型的本体坐标系，所述原始准乘波体构型为原始尖头前缘准乘波体构型，所述原始前缘线为具有拐点且不可导的前缘线；
将所述首部端点与其中一个所述尾部端点连接形成直线段；
经过位于所述直线段上的所述尾部端点形成与Y轴平行的直线；
根据所述直线段以及所述直线形成投影平面；
通过沿Z轴方向将所述原始前缘线向所述投影平面投影，对所述原始前缘线进行优化形成连续可导的目标前缘线。
优选地：在所述目标前缘线选取多个点，从每个点沿每个纵向截面采用相同的曲线方程生成多条目标截面型线且均直到底面截止；
根据所述目标前缘线以及多条所述目标截面型线生成目标准乘波体构型的下表面。
优选地：所述曲线方程为与生成所述原始准乘波体构型相同的曲线方程。
优选地：所述曲线方程的表达式如下：
Figure BDA0002976070870000031
式中，L为给定的乘波体总长，ai方程系数，bi＝0.5+0.1i(1≤i≤11)。
优选地：所述本体坐标系为获得所述原始准乘波体构型的原始坐标系。
优选地：所述原始坐标系的X轴的指向为由原始准乘波体构型的尾部指向头部、Y轴的指向为在水平面上指向左侧、Z轴的指向为在铅锤面内指向下方。
优选地：所述原始准乘波体构型为以配平和纵向稳定性为优化目标获得的尖头前缘准乘波体。
一种平面前缘升力体前缘线优化系统，所述系统包括：
端点获取机构，用于在计算坐标系下获取原始准乘波体构型的原始前缘线的首部端点以及多个尾部端点；所述计算坐标系为原始准乘波体构型的本体坐标系，所述原始准乘波体构型为原始尖头前缘准乘波体构型，所述原始前缘线为具有拐点且不可导的前缘线；
直线段形成机构，用于将所述首部端点与其中一个所述尾部端点连接形成直线段；
直线形成机构，用于经过位于所述直线段上的所述尾部端点形成与Y轴平行的直线；
投影面形成机构，用于根据所述直线段以及所述直线形成投影平面；
目标前缘线形成机构，用于通过沿Z轴方向将所述原始前缘线向所述投影平面投影，对所述原始前缘线进行优化形成连续可导的目标前缘线。
一种平面前缘升力体构型，包括：
前缘线，所述前缘线为将原始前缘线在投影平面上投影获得的连续可导的目标前缘线；
其中，所述原始前缘线为原始尖头前缘准乘波体构型的具有拐点且不可导的前缘线；所述投影平面为通过直线段以及直线形成的平面，所述直线段为所述原始前缘线的首部端点与其中一个尾部端点连接形成的直线段，所述直线为经过位于所述直线段上的所述尾部端点形成与Y轴平行的直线，所述Y轴为原始准乘波体构型的本体坐标系中的Y轴。
优选地：下表面，所述下表面为采用与所述目标前缘线相交的多条截面型线生成的表面；所述截面型线与生成所述原始尖头前缘准乘波体构型的下表面的多条截面型线相同。
优选地：上表面，所述上表面为自由流面。
根据本发明提供的具体实施例，本发明公开了以下技术效果：
本申请实施例提供的平面前缘升力体前缘线优化方法，可以将原始尖头准乘波体构型的前缘线在指定平面上投影生成新的前缘线，可以克服原始尖头准乘波体构型外形畸变、不利于工程加工制造的问题。
另外，在优选地实施方式下，对原始尖头准乘波体构型进行修型，在新前缘基础上，保持下表面生成方法和截面型线不变，从而生成新的平面前缘升力体构型，对升力体的升阻比影响小，既保留了尖头前缘可以降低热流的优点，又克服了其外形畸变严重、不利于工程加工制造的难题，同时还可保持该类构型的高升阻比特性，为高超声速高升阻比乘波类构型的工程化应用提供了一种有效的设计方法。
当然，实施本发明的任一产品并不一定需要同时达到以上所述的所有优点。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种平面前缘升力体前缘线优化方法的流程图；
图2(a-d)是典型圆头前缘准乘波体构型示意图；
图3(a-d)是典型尖头前缘准乘波体构型示意图；
图4(a-b)是典型锥导乘波体生成示意图；
图5是准乘波体下表面生成示意图；
图6是本发明实施例提供的目标前缘线生成过程示意图；
图7(a-d)是本发明实施例提供的根据目标前缘线生成的平面前缘线升力体构型示意图；
图8(a-d)是本发明实施例提供的原始尖头准乘波体构型示意图；
图9是本发明实施例提供的原始尖头准乘波体构型的升阻比线形图；
图10(a-d)是本发明实施例提供的在原始尖头准乘波体构型基础上采用目标前缘线生成的平面前缘升力体构型示意图；
图11是本发明实施例提供的平面前缘升力体构型的升阻比线形图；
图12是本发明实施例提供的一种平面前缘升力体前缘线优化系统的示意图。
图中：前缘线1、原始前缘线11、目标前缘线12、截面型线2、投影平面3、激波4、轴对称体5、准乘波体构型6、基准平面7。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员所获得的所有其他实施例，都属于本发明保护的范围。
实施例
参见图1，为本发明实施例提供的一种平面前缘升力体前缘线优化方法，如图1所示，该方法可以包括：
S101：在计算坐标系下获取原始准乘波体构型的原始前缘线的首部端点以及多个尾部端点；所述计算坐标系为原始准乘波体构型的本体坐标系，所述原始准乘波体构型为原始尖头前缘准乘波体构型，所述原始前缘线为具有拐点且不可导的前缘线；在实际应用中，所述原始准乘波体构型为以配平和纵向稳定性为优化目标获得的尖头前缘准乘波体。该本体坐标系可以根据原始准乘波体构型当前所处的姿态所确定，为了方便进行计算，所述本体坐标系可以为获得所述原始准乘波体构型的原始坐标系。具体的，所述原始坐标系的X轴的指向为由原始准乘波体构型的尾部指向头部、Y轴的指向为在水平面上指向左侧、Z轴的指向为在铅锤面内指向下方。
S102：将所述首部端点与其中一个所述尾部端点连接形成直线段；由于准乘波体具有左右对称的特性，因此在实际应用中，为了简化操作可以仅对左右半侧中一个进行优化，然后镜像生成另一半即可。
S103：经过位于所述直线段上的所述尾部端点形成与Y轴平行的直线；在实际应用中，需要保证直线段与直线具有相同的交点。
S104：根据所述直线段以及所述直线形成投影平面；由于直线段与直线具有相同的交点，因此可以形成一个平面作为投影平面。
S105：通过沿Z轴方向将所述原始前缘线向所述投影平面投影，对所述原始前缘线进行优化形成连续可导的目标前缘线。
本申请实施例提供的方法，可以将原始尖头准乘波体构型的前缘线在指定平面上投影生成新的前缘线，在新前缘基础上，保持下表面生成方法和截面型线不变，从而生成新的平面前缘升力体构型，可以克服原始尖头准乘波体构型外形畸变、不利于工程加工制造的问题。
本申请实施例提供的方法是在基于传统生成准乘波体方法的基础上进行的改进，因此，在目标前缘线生成后，可以采用现有技术中任何一种可以生成准乘波体的方法。例如，本申请实施例可以提供在所述目标前缘线选取多个点，从每个点沿每个纵向截面采用相同的曲线方程生成多条目标截面型线且均直到底面截止；根据所述目标前缘线以及多条所述目标截面型线生成目标准乘波体构型的下表面。采用在目标前缘线上选取点，再由各个点为起点采用相同的曲线方程生成目标截面型线，可以保证生成的各条目标截面型线均与目标前缘线相交。进而保证获得的目标准乘波体构型的下表面的流场在理论上和基准流场的对应部分完全相同，激波可以完全附着在前缘线处，阻止下表面波后高压气体上溢，在小攻角即可获得较高升力。
可以理解的是，该曲线方程可以是根据需要定义的任意的曲线方程，在实际操作时，只需要保证生成各条截面型线时采用的曲线方程相同即可。例如，在一种实现方式下，本申请实施例可以提供所述曲线方程为与生成所述原始准乘波体构型相同的曲线方程。采用与原始准乘波体构型相同的曲线方程生成的下表面，可以保留原始准乘波体构型的下表面的大部分特征，变化较小，便于快速得到满足升阻比要求的升力体。
具体的，所述曲线方程的表达式如下：
Figure BDA0002976070870000071
式中，L为给定的乘波体总长，ai方程系数，bi＝0.5+0.1i(1≤i≤11)。
总之，本申请实施例提供的平面前缘升力体前缘线优化方法，可以对原始尖头准乘波体构型进行修型，生成平面前缘升力体构型，既保留了尖头前缘可以降低热流的优点，又克服了其外形畸变严重、不利于工程加工制造的难题，同时还可保持该类构型的高升阻比特性，为高超声速高升阻比乘波类构型的工程化应用提供了一种有效的设计方法。
下面通过具体实施方式对本申请实施例提供的方法进行详细介绍。
首先，对原始准乘波体构型的具体生成步骤进行详细介绍。
具体设计步骤如下：
步骤一、前缘线定义。
首先，根据锥导乘波体生成前缘线，其生成过程如图4(a-b)所示。给定来流马赫数Ma和锥形激波的激波角βs，可以得到三维激波面(Generating Shock)，在激波面的底面(Base Surface)YZ平面上定义一条任意形式的基准曲线，此处采用如下形式的三次多项式来定义底面上的基准曲线：
Figure BDA0002976070870000081
为了更清晰直观地描述该曲线，设激波圆半径为Rs，基准曲线Z截距为R0，方位角为
Figure BDA0002976070870000082
令曲线与激波圆交点位置处的切线与Y轴的夹角为η，曲线与Z轴交点位置处的切线与Y轴的夹角为λ，并令参数kw＝R0/Rs，化简后可得：
Figure BDA0002976070870000083
这样，给定设计参数kw、
Figure BDA0002976070870000084
η和λ，即可完全确定式(1)的基准曲线形式；再将该曲线沿X轴方向朝激波面投影得到一条交线，即为乘波体的前缘线(Leading Edge)；
步骤二、原始准乘波体构型生成
在上述前缘线上选取若干个点，从每个点，沿每个纵向截面采用相同的曲线形式生成该截面型线，直到底面截止，曲线表达式如下：
Figure BDA0002976070870000085
式中，L为给定的准乘波体总长，ai方程系数，bi＝0.5+0.1i(1≤i≤11)。
根据前缘线和截面型线，即可生成准乘波体下表面，如图5所示。
上表面采用与自由来流方向平行的自由流面。最终生成的原始准乘波体构型如下图3(a-d)所示。
由图3(a-d)中可见，生成的原始尖头前缘准乘波体构型在对称面处前缘线有拐点、不可导，这样使得在工程应用中，尖头前缘准乘波体外形不利于加工制造。
为了解决上述问题，可以采用本申请实施例提供的方法，对原始尖头前缘准乘波体构型的原始前缘线进行优化，具体的实现过程为：
考虑到左右两侧外形的对称性，仅以左半侧说明该构型的设计过程，如图6所示。
①设原始前缘线端点为点A和点B，连结两点得到线段AB；
②过B点作与Y轴平面的直线BC；
③根据AB和BC两条直线确定投影平面；
④将原始前缘线沿Z轴方向朝投影平面投影，即可获得平面前缘线；
⑤在平面前缘线上选取若干点，采用与原始准乘波体构型相同的截面型线生成下表面，上表面同样为自由流面，最终得到的平面前缘升力体构型如图7(a-d)所示。由图7(a-d)中可见，优化后获得的目标前缘线连续可导，有利于在工程应用中的加工制造。
采用本申请实施例提供的方法对原始前缘线进行优化后，即可以获得连续可导的前缘线，同时对于原始尖头前缘准乘波体构型的升阻比影响不大，完全可以满足生成的升力体对于高升阻比的要求。下面对本申请实施例对升阻比影响进行验证。
飞行工况：马赫数20，飞行高度40km。
飞行器长度：5.5米。
生成原始准乘波体构型：
原始前缘线设计参数：kw＝0.717、φ＝33.457°、η＝37.554°和λ＝25.183°。
下表面型线系数(方程3)：
Figure BDA0002976070870000101
采用上述原始前缘线设计参数以及下表面型线系数生成的原始尖头准乘波体构型如图8(a-d)所示，生成的原始尖头前缘准乘波体构型在对称面处前缘线有拐点、不可导。图9给出了该外形在20马赫40km高度的升阻比，在1°攻角升阻比最大，为5.44。
采用申请实施例提供的方法对原始尖头准乘波体构型的前缘线进行优化，并根据优化后获得的目标前缘线保持下表面型线系数以及曲线方程不变生成目标准乘波体构型，如图10所示，优化后获得的目标前缘线连续可导，有利于在工程应用中的加工制造。图11给出了该外形在20马赫40km高度的升阻比，该构型在1°攻角升阻比最大，达到5.5，略高于原始构型。满足升阻比的设计要求。
参见图12，与本申请实施例提供的一种平面前缘升力体前缘线优化方法相对应，如图12所示，本申请实施例还提供了一种平面前缘升力体前缘线优化系统，该系统具体可以包括：
端点获取机构201，用于在计算坐标系下获取原始准乘波体构型的原始前缘线的首部端点以及多个尾部端点；所述计算坐标系为原始准乘波体构型的本体坐标系，所述原始准乘波体构型为原始尖头前缘准乘波体构型，所述原始前缘线为具有拐点且不可导的前缘线；
直线段形成机构202，用于将所述首部端点与其中一个所述尾部端点连接形成直线段；
直线形成机构203，用于经过位于所述直线段上的所述尾部端点形成与Y轴平行的直线；
投影面形成机构204，用于根据所述直线段以及所述直线形成投影平面；
目标前缘线形成机构205，用于通过沿Z轴方向将所述原始前缘线向所述投影平面投影，对所述原始前缘线进行优化形成连续可导的目标前缘线。
与本申请实施例提供的一种平面前缘升力体前缘线优化方法以及系统相对应，本申请实施例还可以提供一种平面前缘升力体构型，包括：
前缘线，所述前缘线为将原始前缘线在投影平面上投影获得的连续可导的目标前缘线；
其中，所述原始前缘线为原始尖头前缘准乘波体构型的具有拐点且不可导的前缘线；所述投影平面为通过直线段以及直线形成的平面，所述直线段为所述原始前缘线的首部端点与其中一个尾部端点连接形成的直线段，所述直线为经过位于所述直线段上的所述尾部端点形成与Y轴平行的直线，所述Y轴为原始准乘波体构型的本体坐标系中的Y轴。
进一步的，在一种实现方式下，还可以包括下表面，所述下表面为采用与所述目标前缘线相交的多条截面型线生成的表面；所述截面型线与生成所述原始尖头前缘准乘波体构型的下表面的多条截面型线相同。
再进一步的，在另一种实现方式下，还可以包括上表面，所述上表面为自由流面。
该平面前缘升力体构型具有设计要求的升阻比，同时由于前缘线经过优化后具有连续可导的优点，因此，有利于实际工程制造。
需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
通过以上的实施方式的描述可知，本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统或系统实施例而言，由于其基本相似于方法实施例，所以描述得比较简单，相关之处参见方法实施例的部分说明即可。以上所描述的系统及系统实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性劳动的情况下，即可以理解并实施。
以上所述仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内所作的任何修改、等同替换、改进等，均包含在本发明的保护范围内。

Claims (11)

1.一种平面前缘升力体前缘线优化方法，其特征在于，所述方法包括：
在计算坐标系下获取原始准乘波体构型的原始前缘线的首部端点以及多个尾部端点；所述计算坐标系为原始准乘波体构型的本体坐标系，所述原始准乘波体构型为原始尖头前缘准乘波体构型，所述原始前缘线为具有拐点且不可导的前缘线；
将所述首部端点与其中一个所述尾部端点连接形成直线段；
经过位于所述直线段上的所述尾部端点形成与Y轴平行的直线；
根据所述直线段以及所述直线形成投影平面；
通过沿Z轴方向将所述原始前缘线向所述投影平面投影，对所述原始前缘线进行优化形成连续可导的目标前缘线。
2.根据权利要求1所述的平面前缘升力体前缘线优化方法，其特征在于，还包括：
在所述目标前缘线选取多个点，从每个点沿每个纵向截面采用相同的曲线方程生成多条目标截面型线且均直到底面截止；
根据所述目标前缘线以及多条所述目标截面型线生成目标准乘波体构型的下表面。
3.根据权利要求2所述的平面前缘升力体前缘线优化方法，其特征在于，所述曲线方程为与生成所述原始准乘波体构型相同的曲线方程。
4.根据权利要求3所述的平面前缘升力体前缘线优化方法，其特征在于，所述曲线方程的表达式如下：
Figure FDA0002976070860000011
式中，L为给定的乘波体总长，ai方程系数，bi＝0.5+0.1i(1≤i≤11)。
5.根据权利要求1所述的平面前缘升力体前缘线优化方法，其特征在于，所述本体坐标系为获得所述原始准乘波体构型的原始坐标系。
6.根据权利要求5所述的平面前缘升力体前缘线优化方法，其特征在于，所述原始坐标系的X轴的指向为由原始准乘波体构型的尾部指向头部、Y轴的指向为在水平面上指向左侧、Z轴的指向为在铅锤面内指向下方。
7.根据权利要求1所述的平面前缘升力体前缘线优化方法，其特征在于，所述原始准乘波体构型为以配平和纵向稳定性为优化目标获得的尖头前缘准乘波体。
8.一种平面前缘升力体前缘线优化系统，其特征在于，所述系统包括：
端点获取机构，用于在计算坐标系下获取原始准乘波体构型的原始前缘线的首部端点以及多个尾部端点；所述计算坐标系为原始准乘波体构型的本体坐标系，所述原始准乘波体构型为原始尖头前缘准乘波体构型，所述原始前缘线为具有拐点且不可导的前缘线；
直线段形成机构，用于将所述首部端点与其中一个所述尾部端点连接形成直线段；
直线形成机构，用于经过位于所述直线段上的所述尾部端点形成与Y轴平行的直线；
投影面形成机构，用于根据所述直线段以及所述直线形成投影平面；
目标前缘线形成机构，用于通过沿Z轴方向将所述原始前缘线向所述投影平面投影，对所述原始前缘线进行优化形成连续可导的目标前缘线。
9.一种平面前缘升力体构型，其特征在于，包括：
前缘线，所述前缘线为将原始前缘线在投影平面上投影获得的连续可导的目标前缘线；
其中，所述原始前缘线为原始尖头前缘准乘波体构型的具有拐点且不可导的前缘线；所述投影平面为通过直线段以及直线形成的平面，所述直线段为所述原始前缘线的首部端点与其中一个尾部端点连接形成的直线段，所述直线为经过位于所述直线段上的所述尾部端点形成与Y轴平行的直线，所述Y轴为原始准乘波体构型的本体坐标系中的Y轴。
10.根据权利要求9所述的平面前缘升力体构型，其特征在于，还包括：
下表面，所述下表面为采用与所述目标前缘线相交的多条截面型线生成的表面；所述截面型线与生成所述原始尖头前缘准乘波体构型的下表面的多条截面型线相同。
11.根据权利要求10所述的平面前缘升力体构型，其特征在于，还包括：
上表面，所述上表面为自由流面。
CN202110274578.6A 2021-03-15 2021-03-15 一种平面前缘升力体前缘线优化方法、系统及升力体 Active CN112948976B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110274578.6A CN112948976B (zh) 2021-03-15 2021-03-15 一种平面前缘升力体前缘线优化方法、系统及升力体

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110274578.6A CN112948976B (zh) 2021-03-15 2021-03-15 一种平面前缘升力体前缘线优化方法、系统及升力体

Publications (2)

Publication Number Publication Date
CN112948976A true CN112948976A (zh) 2021-06-11
CN112948976B CN112948976B (zh) 2024-02-02

Family

ID=76229771

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110274578.6A Active CN112948976B (zh) 2021-03-15 2021-03-15 一种平面前缘升力体前缘线优化方法、系统及升力体

Country Status (1)

Country Link
CN (1) CN112948976B (zh)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105329462A (zh) * 2015-11-16 2016-02-17 中国人民解放军国防科学技术大学 基于可变壁面压力分布规律的吻切流场乘波前体设计方法
CN108100291A (zh) * 2017-12-26 2018-06-01 中国人民解放军国防科技大学 一种给定三维前缘线的吻切乘波体设计方法
CN108304611A (zh) * 2017-12-26 2018-07-20 中国人民解放军国防科技大学 一种给定三维前缘线的锥导乘波体设计方法
CN108595856A (zh) * 2018-04-28 2018-09-28 中国人民解放军国防科技大学 一种“双体”乘波体的参数化设计方法
CN108629108A (zh) * 2018-04-28 2018-10-09 中国人民解放军国防科技大学 一种机身机翼部件可控的参数化乘波体设计方法
CN109250144A (zh) * 2018-09-30 2019-01-22 中国人民解放军国防科技大学 后掠角及上/下反角直接可控的吻切锥乘波体设计方法
CN109279043A (zh) * 2018-10-23 2019-01-29 中国人民解放军国防科技大学 融合低速翼型的冯卡门乘波体设计方法
CN109573093A (zh) * 2018-12-18 2019-04-05 中国人民解放军国防科技大学 融合低速翼型的吻切轴对称冯卡门乘波体设计方法
CN109573092A (zh) * 2018-12-18 2019-04-05 中国人民解放军国防科技大学 吻切轴对称冯卡门乘波体设计方法
US20200283169A1 (en) * 2017-11-09 2020-09-10 China Academy Of Aerospace Aerodynamics Osculating cone theory-based fixed-plane waverider design method

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105329462A (zh) * 2015-11-16 2016-02-17 中国人民解放军国防科学技术大学 基于可变壁面压力分布规律的吻切流场乘波前体设计方法
US20200283169A1 (en) * 2017-11-09 2020-09-10 China Academy Of Aerospace Aerodynamics Osculating cone theory-based fixed-plane waverider design method
CN108100291A (zh) * 2017-12-26 2018-06-01 中国人民解放军国防科技大学 一种给定三维前缘线的吻切乘波体设计方法
CN108304611A (zh) * 2017-12-26 2018-07-20 中国人民解放军国防科技大学 一种给定三维前缘线的锥导乘波体设计方法
CN108595856A (zh) * 2018-04-28 2018-09-28 中国人民解放军国防科技大学 一种“双体”乘波体的参数化设计方法
CN108629108A (zh) * 2018-04-28 2018-10-09 中国人民解放军国防科技大学 一种机身机翼部件可控的参数化乘波体设计方法
CN109250144A (zh) * 2018-09-30 2019-01-22 中国人民解放军国防科技大学 后掠角及上/下反角直接可控的吻切锥乘波体设计方法
CN109279043A (zh) * 2018-10-23 2019-01-29 中国人民解放军国防科技大学 融合低速翼型的冯卡门乘波体设计方法
CN109573093A (zh) * 2018-12-18 2019-04-05 中国人民解放军国防科技大学 融合低速翼型的吻切轴对称冯卡门乘波体设计方法
CN109573092A (zh) * 2018-12-18 2019-04-05 中国人民解放军国防科技大学 吻切轴对称冯卡门乘波体设计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
刘文: "高超声速乘波体气动布局优化及稳定性研究", 《中国优秀博士学位论文全文数据库》, no. 2, pages 031 - 31 *
王晓燕: "基于三维前缘线的乘波体设计方法研究", 《中国优秀硕士学位论文全文数据库》, no. 1, pages 031 - 831 *

Also Published As

Publication number Publication date
CN112948976B (zh) 2024-02-02

Similar Documents

Publication Publication Date Title
CN107180134B (zh) 一种可重复使用的天地往返飞行器外形设计方法
CN107963236B (zh) 基于密切锥理论的定平面乘波体设计方法
EP3484764B1 (en) Airplane wing
CN109279044B (zh) 变马赫数吻切流场乘波体的气动外形设计方法
CN109573093B (zh) 融合低速翼型的吻切轴对称冯卡门乘波体设计方法
CN106126791B (zh) 一种考虑几何不确定性的高超声速机翼气动力/热分析方法
CN104527971B (zh) 一种高超声速飞行器减阻防热反向喷流喷孔设计方法
CN110589010A (zh) 高超声速大装载空间乘波体设计方法
Chen et al. Waverider configuration design with variable shock angle
CN110891857A (zh) 具有至少两个小翼的飞机机翼
CN110104164B (zh) 一种用于跨声速机翼的前加载-吸气组合流动控制方法
CN112948977A (zh) 一种乘波体横航向稳定性优化方法及系统
CN109279043B (zh) 融合低速翼型的冯卡门乘波体设计方法
CN112948976A (zh) 一种平面前缘升力体前缘线优化方法、系统及升力体
Woodward USSAERO computer program development, versions B and C
CN109677630B (zh) 基准流场激波形状可控的强几何约束下的乘波体设计方法
CN112949222B (zh) 一种准乘波体构型横航向稳定性优化方法、系统及升力体
Zhang et al. Supercritical wing design and optimization for transonic civil airplane
CN114611390A (zh) 一种基于全弹道航程的锥导乘波体优化方法及系统
CN110816871A (zh) 一种新的基于锥导法的两级乘波体设计方法
CN114169070A (zh) 一种飞行器的翼型生成方法
Obayashi et al. Transonic wing design by inverse optimization using MOGA
CN112949199B (zh) 一种幂次乘波体的纵向稳定性优化方法及系统
Gan Computational Fluid Dynamics (CFD) of Drag Force for Bullet’s Shape Design
Karimi Kelayeh et al. Aerodynamic assessment of a control strategy based on twist morphing wing in a flying wing aircraft

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