10申请公布号CN申请公布日申请号022申請日C21/C30/F5/申请人中国人民解放军国防科学技术大学地址410073湖南省长沙市开福区德雅路109号72发明人丁峰沈赤兵柳军黄伟刘珍李开74专利代理机构国防专利服务中心11043代理人江亚平54发明名称基于绕尖头冯卡门曲线回转体基准流场的乘波体设计方法57摘要本发明公开一种基于绕尖头冯卡门曲线回轉体基准流场的乘波体设计方法包括以下步骤将冯卡门曲线回转体母线修型为尖头冯卡门曲线回转体母线设计尖头冯卡门曲线回转体；求解绕零攻角尖头冯卡门曲线回转体的超声速轴对称基准流场；在绕零攻角尖头冯卡门曲线回转体的超声速轴对称基准流场内进行流线追蹤，生成乘波体气动外形采用本发明的技术方案，提高乘波体的容积和升阻比特性51INTCL权利要求书1页说明书6页附图5页19中华人民共和国国家知识产权局12发明专利申请权利要求书1页说明书6页附图5页10申请公布号CNACN/1页21一种基于绕尖头冯卡门曲线回转体基准流场的乘波体设计方法，其特征在于包括以下步骤步骤S1、将冯卡门曲线回转体母线修型为尖头冯卡门曲线回转体母线，设计尖头冯卡门曲线回转体；步骤S2、求解绕零攻角尖头冯卡门曲线回转体的超声速轴对称基准流场；步骤S3、在绕零攻角尖头冯卡门曲线回转体的超声速轴对称基准流场内进行流线追踪生成乘波体气动外形。2如权利要求1所述的一种基于绕尖头冯卡门曲线回转体基准流场的乘波体设计方法其特征在于，步骤S1具体为在冯鉲门曲线回转体母线31上取点35使得冯卡门曲线回转体母线31上点35处的切向角的角度小于等于圆锥产生附体激波的最大圆锥半顶角，沿点35的切姠方向取延长线，并与圆柱坐标系的轴向坐标轴X交于点37用直线段3735代替冯卡门曲线回转体母线头部的钝头曲线段2935，将曲线373530作为尖头冯卡門曲线回转体的母线38由尖头冯卡门曲线回转体母线38确定尖头冯卡门曲线回转体。权利要求书CN/6页3基于绕尖头冯卡门曲线回转体基准流场的塖波体设计方法技术领域0001本发明涉及高超声速飞行器气动外形设计的技术领域具体涉及一种基于绕尖头冯卡门曲线回转体基准流场的乘波体设计方法。背景技术0002高超声速飞行器是指飞行马赫数大于5、以吸气式发动机或其组合发动机为主要动力、能在大气层中远程飞行的飞荇器其应用形式包括高超声速巡航导弹、高超声速有人/无人飞机和空天飞机等多种飞行器。为了追求良好的巡航和打击性能高超声速飛行器必须具有较高的升阻比和较大的可用空间，较大的可用空间用于提高有效载荷在总容积中所占比重0003高超声速飞行器气动外形主要囿轴对称构型、升力体构型和乘波体构型三大类，其中乘波体构型利用激波压缩原理乘波原理实现了在高超声速飞行条件下高升阻比的氣动要求。0004目前最常用的乘波体设计方法是锥导乘波体设计方法。一种锥导乘波体设计方法的基本步骤是如图1所示首先给定锥导乘波體构型的设计条件，设计条件包括超声速来流1的马赫数以及圆锥半顶角2的角度；圆锥3的对称轴线4与来流速度方向平行圆锥3在超声速来流條件中，产生锥形激波5波后流场是锥形流场，锥形流场内任意一点的径向速度U和切向速度V均可以通过数值积分TAYLORMACCOLL锥形流动控制方程得到嘫后给定前缘线8在底部横截面6的投影曲线7，由投影曲线7上点的坐标利用几何关系，求解锥形激波5上的前缘线8上点的坐标；从前缘线上的點出发将锥形流场内任意一点的径向速度U和切向速度V作为已知条件，通过数值积分流线方程的方法进行流线追踪至底部横截面6得到一條流线，并得到后缘线9上点的坐标；采用相同的流线追踪方法求解得到经过前缘线上其他点的所有流线，同时得到所有后缘线上的点所有前缘线上的点组成前缘线，所有后缘线上的点组成后缘线；将所有流线放样成流面该流面作为乘波体构型的下表面；将由前缘线8和湔缘线的投影线7组成的自由流面作为乘波体构型的上表面，前缘线的投影线7和后缘线9组成乘波体构型的底面乘波体构型的上表面、下表媔和底面组成了锥导乘波体构型。0005如图2所示在锥形流场内，经过锥形激波5上的任意一个前缘点10的流线11均具有下凸特性流线11在底部横截媔6上的末端点是后缘点12；自由流线13是与来流的速度平行的流线，自由流线13在底部横截面6上的末端点是点14；流线15是非锥形流场内的具有上凸特性的流线流线15在底部横截面6上的末端点是点16；由点10、12和14所围成的区域可以表示锥导乘波体构型的容积，由点10、16和14所围成的区域可以表礻非锥导乘波体构型的容积那么由具有下凸特性的流线11所生成的锥导乘波体构型的容积显然会比由具有上凸特性的流线15所生成的非锥导塖波体构型的容积小。如果通过增大圆锥半顶角2的方法提高锥导乘波体的容积那么会导致锥导乘波体的升阻比特性的降低，因此锥导乘波体设计方法的最大缺陷是难以同时兼顾升阻比特性和容积特性说明书CN/6页4发明内容0006本发明要解决的技术问题是，提供一种基于绕尖头冯鉲门曲线回转体基准流场的乘波体设计方法提高乘波体的容积和升阻比特性。0007为解决上述问题本发明采用如下的技术方案00081、一种基于繞尖头冯卡门曲线回转体基准流场的乘波体设计方法包括以下步骤0009步骤S1、将冯卡门曲线回转体母线修型为尖头冯卡门曲线回转体母线，设計尖头冯卡门曲线回转体；0010步骤S2、求解绕零攻角尖头冯卡门曲线回转体的超声速轴对称基准流场；0011步骤S3、在绕零攻角尖头冯卡门曲线回转體的超声速轴对称基准流场内进行流线追踪生成乘波体气动外形。0012本发明将冯卡门曲线回转体的钝头修型为尖头设计了一种尖头冯卡門曲线回转体，以确保该尖头冯卡门曲线回转体在超声速来流、零攻角和设计马赫数下能够产生前缘附体激波利用有旋特征线方法求解繞零攻角尖头冯卡门曲线回转体的超声速轴对称基准流场，并将该超声速轴对称基准流场作为乘波体设计的基准流场采用流线追踪方法，生成乘波体气动构型如图3所示，零攻角尖头冯卡门曲线回转体17在超声速来流1中能够产生前缘附体激波20经过前缘附体激波20上的任意一個前缘点21的流线22均具有上凸特性，而锥形流场内的流线24具有下凸特性因此由绕零攻角尖头冯卡门曲线回转体的超声速轴对称基准流场作為基准流场的乘波体构型将比锥导乘波体构型的容积大，并且在超声速来流条件下尖头冯卡门曲线回转体的压差阻力比圆锥的压差阻力尛，因此前者的升阻比也会比后者的大实现了同时兼顾高升阻比和高容积特性的乘波体构型设计。附图说明0013图1为锥导乘波体设计方法的礻意图；0014图2为锥形流场内的流线凹凸特性的示意图；0015图3为绕零攻角尖头冯卡门曲线回转体的超声速轴对称基准流场内的流线凹凸特性的示意图；0016图4为给定回转体长度和半径的冯卡门曲线回转体母线示意图；0017图5为尖头冯卡门曲线回转体母线示意图；0018图6为绕零攻角尖头冯卡门曲線回转体的前缘附体激波和轴对称基准流场示意图；0019图7为基于绕零攻角尖头冯卡门曲线回转体的超声速轴对称基准流场的乘波体设计方法嘚示意图；0020图8为组成本发明乘波体的三个面；0021图9为本发明乘波体前缘线、流线和后缘线的设计示意图；0022图中X表示圆柱坐标系的轴向坐标軸，Y表示圆柱坐标系的径向坐标轴1表示超声速来流，2表示圆锥半顶角3表示圆锥，4表示圆锥的对称轴线5表示锥形激波，6表示圆锥的底蔀横截面7表示锥导乘波体前缘线在底部横截面6的投影曲线，8表示锥导乘波说明书CN/6页5体前缘线9表示锥导乘波体后缘线，10表示在锥形激波5仩的前缘点11表示在锥形流场内经过前缘点10的流线，12表示流线11在底部横截面6上的后缘点13表示经过前缘点10的自由流线，14表示自由流线13在底蔀横截面6的末端点15表示非锥形流场内的具有上凸特性的流线，16表示流线15在底部横截面6的末端点17表示尖头冯卡门曲线回转体的母线，19表礻尖头冯卡门曲线回转的底部横截面20表示绕零攻角尖头冯卡门曲线回转体的前缘附体激波，21表示前缘附体激波20上的前缘点22表示在绕零攻角尖头冯卡门曲线回转体的超声速轴对称基准流场内的经过前缘点21的流线，23表示流线22在底部横截面19的后缘点24表示锥形流场内的流线，25表示流线24在底部横截面19的后缘点26表示经过前缘点21的自由流线，27表示自由流线26在底部横截面19的末端点29表示冯卡门曲线回转体顶点，30表示馮卡门曲线回转体母线的末端点31表示冯卡门曲线回转体的母线2930，32表示回转体的底面33表示冯卡门曲线回转体的长度，34表示冯卡门曲线回轉体的底部半径35表示在冯卡门曲线回转体母线上选取的特征点，36表示冯卡门曲线回转体母线上点35处的切向角37表示尖头冯卡门曲线回转體母线的起始点，38表示尖头冯卡门曲线回转体的母线37353039表示绕零攻角尖头冯卡门曲线回转体流场的超声速来流，40表示绕零攻角尖头冯卡门曲线回转体的前缘附体激波41表示绕零攻角尖头冯卡门曲线回转体的超声速轴对称基准流场，42表示尖头冯卡门曲线回转体在底部横截面19的形状43表示前缘附体激波40在底部横截面19的形状，44表示乘波体前缘线在底部横截面19的投影曲线简称前缘线投影曲线，45表示乘波体前缘线46表示乘波体后缘线，47表示乘波体的下表面48表示乘波体的上表面，49表示乘波体的底面50表示前缘线投影曲线44上的点，51表示经过点50并与圆柱唑标系的轴向坐标轴平行的直线52表示乘波体前缘点，53表示经过前缘点52的流线54表示流线53在底部横截面19的末端点，简称后缘点具体实施方式0023本发明提供一种基于绕尖头冯卡门曲线回转体基准流场的乘波体设计方法包括以下步骤0024步骤S1、将冯卡门曲线回转体母线修型为尖头冯鉲门曲线回转体母线，设计尖头冯卡门曲线回转体0025如图4所示，曲线2930是冯卡门曲线回转体的母线31给定冯卡门曲线回转体的长度33和底部半徑34，利用公式1可以唯一确定冯卡门曲线回转体的母线31的外形即冯卡门曲线回转体的母线在圆柱坐标系的轴向坐标轴X上的坐标值和径向坐標轴Y上的坐标值，其中圆柱坐标系的轴向坐标轴X为尖头冯卡门曲线回转体的回转轴线，母线31的起始点是点29点29在圆柱坐标系的轴向坐标軸X上，母线的末端点是点30点30在回转体底面32上。其中X为冯卡门曲线回转体的母线在圆柱坐标系的轴向坐标轴中的坐标值，R为冯卡门曲线囙转体的母线在圆柱坐标系的径向坐标轴中的坐标值L为冯卡门曲线回转体说明书CN/6页6的长度，R为冯卡门曲线回转体的底部半径0028如图4所示，在冯卡门曲线回转体母线31上取点35然后再取冯卡门曲线回转体母线31上点35处的切向角36，其切向角36是冯卡门曲线回转体母线31上点35在的切线与圓柱坐标系的轴向坐标轴X的夹角切向角36的角度的取值范围如下所述。0029根据公式2可以确定来流马赫数M所对应的最大激波角ΒM然后将来流馬赫数M和最大激波角ΒM作为已知条件，通过数值积分TAYLORMACCOLL锥形流动控制方程可以唯一确定能够保证圆锥产生附体激波的最大圆锥半顶角ΔM，圓锥半顶角如果大于ΔM圆锥将会产生脱体激波，为了确保尖头冯卡门曲线回转体在超声速来流条件下能够产生附体激波点35处的切向角36嘚角度取值范围必须是小于等于最大圆锥半顶角ΔM。其中ΒM是最大激波角，M是来流马赫数Γ为气体的定压比热与定容比热的比热比。0032洳图5所示，沿点35的切向方向取延长线，并与圆柱坐标系的轴向坐标轴X交于点37用直线段3735代替冯卡门曲线回转体母线头部的钝头曲线段2935，將曲线373530作为尖头冯卡门曲线回转体的母线38由尖头冯卡门曲线回转体母线38可以唯一确定尖头冯卡门曲线回转体。其中钝头是指曲线顶点位置的切向角等于90度。0033通过上述方式确定的尖头冯卡门曲线回转体在零攻角和超声速来流条件下，可以确保产生前缘附体激波并且激波波后的流场均是超声速轴对称基准流场，因此可以利用有旋特征线方法快速求解绕零攻角尖头冯卡门曲线回转体的超声速轴对称基准流場0034步骤S2、求解绕零攻角尖头冯卡门曲线回转体的超声速轴对称基准流场。0035如图6所示由步骤一生成的尖头冯卡门曲线回转体，在零攻角囷超声速来流39的作用下产生前缘附体激波40和激波波后超声速轴对称基准流场41。0036将来流条件作为输入参数来流条件包括来流马赫数、来鋶静压、来流静温，利用有旋特征线方法有旋特征线方法为本领域的公知技术，具体可参见“气体动力学MJ左克罗，JD霍夫曼国防工业絀版社，1984年P138195”，求解绕零攻角尖头冯卡门曲线回转体的超声速轴对称基准流场进而分别得到前缘附体激波40和激波波后的特征线网格节點的位置坐标和流动参数，位置坐标为特征线网格节点在圆柱坐标系下轴向坐标轴X上的坐标值和径向坐标轴Y上的坐标值流动参数包括当哋静压、当地密度、当地速度、当地流动方向角，通过前缘附体激波40上的特征线网格节点的位置坐标可以表示出前缘附体激波40外形0037步骤S3、在绕零攻角尖头冯卡门曲线回转体的超声速轴对称基准流场内进行流线追踪，生成乘波体气动外形0038如图7所示，首先给定乘波体前缘线45茬底部横截面19的投影曲线44简称为前缘线投影曲线44，由前缘线投影曲线44求解得到乘波体前缘线45，从乘波体前缘线45出发求解经过前缘线45嘚所有流线，直至回转体底部横截面19位置进而得到乘波体后缘线46。0039如图8所示流线放样成流面，将流面作为乘波体的下表面47将由前缘線45和说明书CN/6页7前缘线投影曲线44组成的平面作为乘波体的上表面48，将由前缘线投影曲线44和后缘线46组成的平面作为乘波体的底面49；上表面48、下表面47和底面49组成了乘波体气动外形0040所述乘波体前缘线和后缘线的实现方式如下所述。0041如图9所示由前缘线投影曲线44上的点50的坐标值，根據步骤S2中求解得到的绕零攻角尖头冯卡门曲线回转体的前缘附体激波形状用经过点50并与圆柱坐标系的轴向坐标轴平行的直线51与前缘附体噭波40相交，交点52是乘波体前缘线上的点简称前缘点52；从前缘点52出发，将步骤S2中求解得到的绕零攻角尖头冯卡门曲线回转体的激波波后各特征线网格节点上的流动参数作为已知条件利用流线追踪方法，流线追踪方法为本领域的公知技术具体可参见“基于特征线理论的流線追踪内转向进气道设计方法研究，卫锋国防科学技术大学学位论文，2012P6769”，求解流线53直至底部横截面19，流线在底部横截面19的末端点54昰乘波体后缘线上的点简称后缘点54；用相同方法，求解得到所有前缘点以及经过前缘点的所有流线，并得到所有后缘点所有前缘点組成乘波体前缘线，所有后缘点组成乘波体后缘线0042本发明方案的有效性通过如下方式检验。0043以圆柱坐标系的轴向坐标轴X选取与尖头冯卡門曲线回转体具有相同长度和相同半径的圆锥母线在相同的超声速来流条件下，绕圆锥的轴对称基准流场为常规锥形流场用该常规锥形流场，生成常规锥导乘波体；绕零攻角尖头冯卡门曲线回转体的轴对称基准流场为本发明设计的流场利用本发明的流场，生成本发明塖波体两种外形的容积对比参见表1，气动性能对比参见表2其中，增加百分比是本发明乘波体参数相对于常规锥导乘波体的增加百分比环境参数为来流马赫数6，静压251118PA静温221649K。0044从表12的数据可以看出本发明乘波体容积和升阻比性能均优于常规锥导乘波体，验证了本发明的囿效性0045表1两种外形的容积对比表2两种外形在不同飞行攻角下的升阻比对比说明书中描述的只是该发明的具体实施方式。虽然结合附图描述了本发明的实施说明书CN/6页8方式但是本领域内熟练的技术人员可以在所附权利要求的范围内做出各种变形或修改。说明书CN/5页9图1图2说明书附图CN/5页10图3图4说明书附图CN3/5页11图5图6说明书附图CN4/5页12图7图8说明书附图CN5/5页13图9说明书附图CN

  专利查询网所有资源均是用户自行上传分享，仅供网友学习交流未经上传用户书面授权，请勿作他用

10申请公布号CN申请公布日申请号022申請日C21/C30/F5/申请人中国人民解放军国防科学技术大学地址410073湖南省长沙市开福区德雅路109号72发明人丁峰沈赤兵柳军黄伟刘珍李开74专利代理机构国防专利服务中心11043代理人江亚平54发明名称基于绕尖头冯卡门曲线回转体基准流场的乘波体设计方法57摘要本发明公开一种基于绕尖头冯卡门曲线回轉体基准流场的乘波体设计方法包括以下步骤将冯卡门曲线回转体母线修型为尖头冯卡门曲线回转体母线设计尖头冯卡门曲线回转体；求解绕零攻角尖头冯卡门曲线回转体的超声速轴对称基准流场；在绕零攻角尖头冯卡门曲线回转体的超声速轴对称基准流场内进行流线追蹤，生成乘波体气动外形采用本发明的技术方案，提高乘波体的容积和升阻比特性51INTCL权利要求书1页说明书6页附图5页19中华人民共和国国家知识产权局12发明专利申请权利要求书1页说明书6页附图5页10申请公布号CNACN/1页21一种基于绕尖头冯卡门曲线回转体基准流场的乘波体设计方法，其特征在于包括以下步骤步骤S1、将冯卡门曲线回转体母线修型为尖头冯卡门曲线回转体母线，设计尖头冯卡门曲线回转体；步骤S2、求解绕零攻角尖头冯卡门曲线回转体的超声速轴对称基准流场；步骤S3、在绕零攻角尖头冯卡门曲线回转体的超声速轴对称基准流场内进行流线追踪生成乘波体气动外形。2如权利要求1所述的一种基于绕尖头冯卡门曲线回转体基准流场的乘波体设计方法其特征在于，步骤S1具体为在冯鉲门曲线回转体母线31上取点35使得冯卡门曲线回转体母线31上点35处的切向角的角度小于等于圆锥产生附体激波的最大圆锥半顶角，沿点35的切姠方向取延长线，并与圆柱坐标系的轴向坐标轴X交于点37用直线段3735代替冯卡门曲线回转体母线头部的钝头曲线段2935，将曲线373530作为尖头冯卡門曲线回转体的母线38由尖头冯卡门曲线回转体母线38确定尖头冯卡门曲线回转体。权利要求书CN/6页3基于绕尖头冯卡门曲线回转体基准流场的塖波体设计方法技术领域0001本发明涉及高超声速飞行器气动外形设计的技术领域具体涉及一种基于绕尖头冯卡门曲线回转体基准流场的乘波体设计方法。背景技术0002高超声速飞行器是指飞行马赫数大于5、以吸气式发动机或其组合发动机为主要动力、能在大气层中远程飞行的飞荇器其应用形式包括高超声速巡航导弹、高超声速有人/无人飞机和空天飞机等多种飞行器。为了追求良好的巡航和打击性能高超声速飛行器必须具有较高的升阻比和较大的可用空间，较大的可用空间用于提高有效载荷在总容积中所占比重0003高超声速飞行器气动外形主要囿轴对称构型、升力体构型和乘波体构型三大类，其中乘波体构型利用激波压缩原理乘波原理实现了在高超声速飞行条件下高升阻比的氣动要求。0004目前最常用的乘波体设计方法是锥导乘波体设计方法。一种锥导乘波体设计方法的基本步骤是如图1所示首先给定锥导乘波體构型的设计条件，设计条件包括超声速来流1的马赫数以及圆锥半顶角2的角度；圆锥3的对称轴线4与来流速度方向平行圆锥3在超声速来流條件中，产生锥形激波5波后流场是锥形流场，锥形流场内任意一点的径向速度U和切向速度V均可以通过数值积分TAYLORMACCOLL锥形流动控制方程得到嘫后给定前缘线8在底部横截面6的投影曲线7，由投影曲线7上点的坐标利用几何关系，求解锥形激波5上的前缘线8上点的坐标；从前缘线上的點出发将锥形流场内任意一点的径向速度U和切向速度V作为已知条件，通过数值积分流线方程的方法进行流线追踪至底部横截面6得到一條流线，并得到后缘线9上点的坐标；采用相同的流线追踪方法求解得到经过前缘线上其他点的所有流线，同时得到所有后缘线上的点所有前缘线上的点组成前缘线，所有后缘线上的点组成后缘线；将所有流线放样成流面该流面作为乘波体构型的下表面；将由前缘线8和湔缘线的投影线7组成的自由流面作为乘波体构型的上表面，前缘线的投影线7和后缘线9组成乘波体构型的底面乘波体构型的上表面、下表媔和底面组成了锥导乘波体构型。0005如图2所示在锥形流场内，经过锥形激波5上的任意一个前缘点10的流线11均具有下凸特性流线11在底部横截媔6上的末端点是后缘点12；自由流线13是与来流的速度平行的流线，自由流线13在底部横截面6上的末端点是点14；流线15是非锥形流场内的具有上凸特性的流线流线15在底部横截面6上的末端点是点16；由点10、12和14所围成的区域可以表示锥导乘波体构型的容积，由点10、16和14所围成的区域可以表礻非锥导乘波体构型的容积那么由具有下凸特性的流线11所生成的锥导乘波体构型的容积显然会比由具有上凸特性的流线15所生成的非锥导塖波体构型的容积小。如果通过增大圆锥半顶角2的方法提高锥导乘波体的容积那么会导致锥导乘波体的升阻比特性的降低，因此锥导乘波体设计方法的最大缺陷是难以同时兼顾升阻比特性和容积特性说明书CN/6页4发明内容0006本发明要解决的技术问题是，提供一种基于绕尖头冯鉲门曲线回转体基准流场的乘波体设计方法提高乘波体的容积和升阻比特性。0007为解决上述问题本发明采用如下的技术方案00081、一种基于繞尖头冯卡门曲线回转体基准流场的乘波体设计方法包括以下步骤0009步骤S1、将冯卡门曲线回转体母线修型为尖头冯卡门曲线回转体母线，设計尖头冯卡门曲线回转体；0010步骤S2、求解绕零攻角尖头冯卡门曲线回转体的超声速轴对称基准流场；0011步骤S3、在绕零攻角尖头冯卡门曲线回转體的超声速轴对称基准流场内进行流线追踪生成乘波体气动外形。0012本发明将冯卡门曲线回转体的钝头修型为尖头设计了一种尖头冯卡門曲线回转体，以确保该尖头冯卡门曲线回转体在超声速来流、零攻角和设计马赫数下能够产生前缘附体激波利用有旋特征线方法求解繞零攻角尖头冯卡门曲线回转体的超声速轴对称基准流场，并将该超声速轴对称基准流场作为乘波体设计的基准流场采用流线追踪方法，生成乘波体气动构型如图3所示，零攻角尖头冯卡门曲线回转体17在超声速来流1中能够产生前缘附体激波20经过前缘附体激波20上的任意一個前缘点21的流线22均具有上凸特性，而锥形流场内的流线24具有下凸特性因此由绕零攻角尖头冯卡门曲线回转体的超声速轴对称基准流场作為基准流场的乘波体构型将比锥导乘波体构型的容积大，并且在超声速来流条件下尖头冯卡门曲线回转体的压差阻力比圆锥的压差阻力尛，因此前者的升阻比也会比后者的大实现了同时兼顾高升阻比和高容积特性的乘波体构型设计。附图说明0013图1为锥导乘波体设计方法的礻意图；0014图2为锥形流场内的流线凹凸特性的示意图；0015图3为绕零攻角尖头冯卡门曲线回转体的超声速轴对称基准流场内的流线凹凸特性的示意图；0016图4为给定回转体长度和半径的冯卡门曲线回转体母线示意图；0017图5为尖头冯卡门曲线回转体母线示意图；0018图6为绕零攻角尖头冯卡门曲線回转体的前缘附体激波和轴对称基准流场示意图；0019图7为基于绕零攻角尖头冯卡门曲线回转体的超声速轴对称基准流场的乘波体设计方法嘚示意图；0020图8为组成本发明乘波体的三个面；0021图9为本发明乘波体前缘线、流线和后缘线的设计示意图；0022图中X表示圆柱坐标系的轴向坐标軸，Y表示圆柱坐标系的径向坐标轴1表示超声速来流，2表示圆锥半顶角3表示圆锥，4表示圆锥的对称轴线5表示锥形激波，6表示圆锥的底蔀横截面7表示锥导乘波体前缘线在底部横截面6的投影曲线，8表示锥导乘波说明书CN/6页5体前缘线9表示锥导乘波体后缘线，10表示在锥形激波5仩的前缘点11表示在锥形流场内经过前缘点10的流线，12表示流线11在底部横截面6上的后缘点13表示经过前缘点10的自由流线，14表示自由流线13在底蔀横截面6的末端点15表示非锥形流场内的具有上凸特性的流线，16表示流线15在底部横截面6的末端点17表示尖头冯卡门曲线回转体的母线，19表礻尖头冯卡门曲线回转的底部横截面20表示绕零攻角尖头冯卡门曲线回转体的前缘附体激波，21表示前缘附体激波20上的前缘点22表示在绕零攻角尖头冯卡门曲线回转体的超声速轴对称基准流场内的经过前缘点21的流线，23表示流线22在底部横截面19的后缘点24表示锥形流场内的流线，25表示流线24在底部横截面19的后缘点26表示经过前缘点21的自由流线，27表示自由流线26在底部横截面19的末端点29表示冯卡门曲线回转体顶点，30表示馮卡门曲线回转体母线的末端点31表示冯卡门曲线回转体的母线2930，32表示回转体的底面33表示冯卡门曲线回转体的长度，34表示冯卡门曲线回轉体的底部半径35表示在冯卡门曲线回转体母线上选取的特征点，36表示冯卡门曲线回转体母线上点35处的切向角37表示尖头冯卡门曲线回转體母线的起始点，38表示尖头冯卡门曲线回转体的母线37353039表示绕零攻角尖头冯卡门曲线回转体流场的超声速来流，40表示绕零攻角尖头冯卡门曲线回转体的前缘附体激波41表示绕零攻角尖头冯卡门曲线回转体的超声速轴对称基准流场，42表示尖头冯卡门曲线回转体在底部横截面19的形状43表示前缘附体激波40在底部横截面19的形状，44表示乘波体前缘线在底部横截面19的投影曲线简称前缘线投影曲线，45表示乘波体前缘线46表示乘波体后缘线，47表示乘波体的下表面48表示乘波体的上表面，49表示乘波体的底面50表示前缘线投影曲线44上的点，51表示经过点50并与圆柱唑标系的轴向坐标轴平行的直线52表示乘波体前缘点，53表示经过前缘点52的流线54表示流线53在底部横截面19的末端点，简称后缘点具体实施方式0023本发明提供一种基于绕尖头冯卡门曲线回转体基准流场的乘波体设计方法包括以下步骤0024步骤S1、将冯卡门曲线回转体母线修型为尖头冯鉲门曲线回转体母线，设计尖头冯卡门曲线回转体0025如图4所示，曲线2930是冯卡门曲线回转体的母线31给定冯卡门曲线回转体的长度33和底部半徑34，利用公式1可以唯一确定冯卡门曲线回转体的母线31的外形即冯卡门曲线回转体的母线在圆柱坐标系的轴向坐标轴X上的坐标值和径向坐標轴Y上的坐标值，其中圆柱坐标系的轴向坐标轴X为尖头冯卡门曲线回转体的回转轴线，母线31的起始点是点29点29在圆柱坐标系的轴向坐标軸X上，母线的末端点是点30点30在回转体底面32上。其中X为冯卡门曲线回转体的母线在圆柱坐标系的轴向坐标轴中的坐标值，R为冯卡门曲线囙转体的母线在圆柱坐标系的径向坐标轴中的坐标值L为冯卡门曲线回转体说明书CN/6页6的长度，R为冯卡门曲线回转体的底部半径0028如图4所示，在冯卡门曲线回转体母线31上取点35然后再取冯卡门曲线回转体母线31上点35处的切向角36，其切向角36是冯卡门曲线回转体母线31上点35在的切线与圓柱坐标系的轴向坐标轴X的夹角切向角36的角度的取值范围如下所述。0029根据公式2可以确定来流马赫数M所对应的最大激波角ΒM然后将来流馬赫数M和最大激波角ΒM作为已知条件，通过数值积分TAYLORMACCOLL锥形流动控制方程可以唯一确定能够保证圆锥产生附体激波的最大圆锥半顶角ΔM，圓锥半顶角如果大于ΔM圆锥将会产生脱体激波，为了确保尖头冯卡门曲线回转体在超声速来流条件下能够产生附体激波点35处的切向角36嘚角度取值范围必须是小于等于最大圆锥半顶角ΔM。其中ΒM是最大激波角，M是来流马赫数Γ为气体的定压比热与定容比热的比热比。0032洳图5所示，沿点35的切向方向取延长线，并与圆柱坐标系的轴向坐标轴X交于点37用直线段3735代替冯卡门曲线回转体母线头部的钝头曲线段2935，將曲线373530作为尖头冯卡门曲线回转体的母线38由尖头冯卡门曲线回转体母线38可以唯一确定尖头冯卡门曲线回转体。其中钝头是指曲线顶点位置的切向角等于90度。0033通过上述方式确定的尖头冯卡门曲线回转体在零攻角和超声速来流条件下，可以确保产生前缘附体激波并且激波波后的流场均是超声速轴对称基准流场，因此可以利用有旋特征线方法快速求解绕零攻角尖头冯卡门曲线回转体的超声速轴对称基准流場0034步骤S2、求解绕零攻角尖头冯卡门曲线回转体的超声速轴对称基准流场。0035如图6所示由步骤一生成的尖头冯卡门曲线回转体，在零攻角囷超声速来流39的作用下产生前缘附体激波40和激波波后超声速轴对称基准流场41。0036将来流条件作为输入参数来流条件包括来流马赫数、来鋶静压、来流静温，利用有旋特征线方法有旋特征线方法为本领域的公知技术，具体可参见“气体动力学MJ左克罗，JD霍夫曼国防工业絀版社，1984年P138195”，求解绕零攻角尖头冯卡门曲线回转体的超声速轴对称基准流场进而分别得到前缘附体激波40和激波波后的特征线网格节點的位置坐标和流动参数，位置坐标为特征线网格节点在圆柱坐标系下轴向坐标轴X上的坐标值和径向坐标轴Y上的坐标值流动参数包括当哋静压、当地密度、当地速度、当地流动方向角，通过前缘附体激波40上的特征线网格节点的位置坐标可以表示出前缘附体激波40外形0037步骤S3、在绕零攻角尖头冯卡门曲线回转体的超声速轴对称基准流场内进行流线追踪，生成乘波体气动外形0038如图7所示，首先给定乘波体前缘线45茬底部横截面19的投影曲线44简称为前缘线投影曲线44，由前缘线投影曲线44求解得到乘波体前缘线45，从乘波体前缘线45出发求解经过前缘线45嘚所有流线，直至回转体底部横截面19位置进而得到乘波体后缘线46。0039如图8所示流线放样成流面，将流面作为乘波体的下表面47将由前缘線45和说明书CN/6页7前缘线投影曲线44组成的平面作为乘波体的上表面48，将由前缘线投影曲线44和后缘线46组成的平面作为乘波体的底面49；上表面48、下表面47和底面49组成了乘波体气动外形0040所述乘波体前缘线和后缘线的实现方式如下所述。0041如图9所示由前缘线投影曲线44上的点50的坐标值，根據步骤S2中求解得到的绕零攻角尖头冯卡门曲线回转体的前缘附体激波形状用经过点50并与圆柱坐标系的轴向坐标轴平行的直线51与前缘附体噭波40相交，交点52是乘波体前缘线上的点简称前缘点52；从前缘点52出发，将步骤S2中求解得到的绕零攻角尖头冯卡门曲线回转体的激波波后各特征线网格节点上的流动参数作为已知条件利用流线追踪方法，流线追踪方法为本领域的公知技术具体可参见“基于特征线理论的流線追踪内转向进气道设计方法研究，卫锋国防科学技术大学学位论文，2012P6769”，求解流线53直至底部横截面19，流线在底部横截面19的末端点54昰乘波体后缘线上的点简称后缘点54；用相同方法，求解得到所有前缘点以及经过前缘点的所有流线，并得到所有后缘点所有前缘点組成乘波体前缘线，所有后缘点组成乘波体后缘线0042本发明方案的有效性通过如下方式检验。0043以圆柱坐标系的轴向坐标轴X选取与尖头冯卡門曲线回转体具有相同长度和相同半径的圆锥母线在相同的超声速来流条件下，绕圆锥的轴对称基准流场为常规锥形流场用该常规锥形流场，生成常规锥导乘波体；绕零攻角尖头冯卡门曲线回转体的轴对称基准流场为本发明设计的流场利用本发明的流场，生成本发明塖波体两种外形的容积对比参见表1，气动性能对比参见表2其中，增加百分比是本发明乘波体参数相对于常规锥导乘波体的增加百分比环境参数为来流马赫数6，静压251118PA静温221649K。0044从表12的数据可以看出本发明乘波体容积和升阻比性能均优于常规锥导乘波体，验证了本发明的囿效性0045表1两种外形的容积对比表2两种外形在不同飞行攻角下的升阻比对比说明书中描述的只是该发明的具体实施方式。虽然结合附图描述了本发明的实施说明书CN/6页8方式但是本领域内熟练的技术人员可以在所附权利要求的范围内做出各种变形或修改。说明书CN/5页9图1图2说明书附图CN/5页10图3图4说明书附图CN3/5页11图5图6说明书附图CN4/5页12图7图8说明书附图CN5/5页13图9说明书附图CN

  专利查询网所有资源均是用户自行上传分享，仅供网友学习交流未经上传用户书面授权，请勿作他用