混压式超声速、高超声速进气道不起动振荡频率预测方法
混压式超声速、高超声速进气道不起动振荡频率预测方法
【技术领域】
[0001] 本发明涉及混压式超声速、高超声速进气道，尤其是其不起动振荡频率的预测方 法。
【背景技术】
[0002] 超声速、高超声速进气道是吸气式高速推进系统的关键部件之一，其包括外压式、 内压式和混压式三类。其中，混压式进气道由于能在保证较高内流气动性能的同时，减小进 气道的气动阻力、轮廓尺寸和结构重量，故在超声速、高超声速进气道的设计中得到了广泛 应用。然而，由于内压缩的引入，使得混压式进气道容易陷入不起动状态。
[0003] 不起动状态是超声速、高超声速进气道的一类非正常工作状态。通常，若混压式超 声速、高超声速进气道的流量捕获特性未因其内部流态的改变而受到影响，称进气道为起 动，否则为不起动。当进气道处于不起动状态时，不仅总压恢复系数和流量系数急剧下降， 使得推进系统不能产生推力，且极可能伴随着波系以及通道压强的剧烈振荡，带来破坏力 极强的周期性力载荷和热载荷，造成结构破坏，并使飞行器在非定常起动力的作用下边的 更加难以控制。为此，有必要揭示混压式进气道振荡流态的产生机理，建立振荡频率的预测 方法，进而为发动机和飞行器的结构强度设计提供参考。
[0004] 自Oswatitsch于1944年首次观测到超声速进气道的振荡流态（即喘振）以来， 国内外有大量的研宄工作致力于这一现象的描述、预测和控制。通过对外压式超声速进气 道喘振的仿真和实验研宄，各国学者发现：扰动的声波反馈在喘振信号传递环中起着重要 作用，故喘振的主特征频率以及各次要特征频率往往与进气道腔体的各阶声学谐振频率相 关。然而，混压式进气道与外压式进气道的不起动流态存在显著差异，前者为超声速溢流， 后者为亚声速溢流，使得其不稳定机理和振荡特性都有所不同。为此，若依旧采用腔体声学 振荡模型对混压式进气道的不起动振荡频率进行估算，必然会带来较大的误差。并且，近年 来混压式进气道在飞行器上的应用趋于广泛，并已衍生发展出了多种不同的具体形式，如 二维形式、轴对称形式、三维侧压形式、内乘波形式以及组合形式等。因此，发展适用于混压 式超声速、高超声速进气道不起动振荡频率的估算方法显得十分重要。
【发明内容】
[0005] 本发明提供一种适用于混压式超声速、高超声速进气道的不起动振荡频率预测方 法，能够实现对不同来流状态和几何参数下进气道不起动振荡频率的快速、准确预测，并且 适用于各种形式的混压式进气道。
[0006] 为达到上述目的，本发明的混压式超声速、高超声速进气道不起动振荡频率预测 方法可采用如下技术方案：
[0007] 混压式超声速、高超声速进气道不起动振荡频率估算方法包括以下步骤：
[0008] (1)、将一个混压式进气道不起动振荡周期分为两个阶段：进气道腔体内高压气体 积蓄阶段和口外不起动波系运动阶段；
[0009] (2)、在已知具体的来流条件和进气道前体几何参数下，通过已有的斜激波理论得 到进气道进口前的气流参数，并换算得到进气道的捕获流量；
[0010] (3)、通过假设进气道处于临界工作状态，在已知进口气体参数和进气道通道面积 变化的基础上，结合正激波理论和一维定常管流理论计算得到该进气道的极限反压；临界 工作状态指结尾激波停于进气道喉道处时的状态；
[0011] (4)、由上述极限反压和已知的来流总温换算出最大储气密度，再乘以进气道内通 道腔体容积得到对应气体积蓄阶段的最大储气量〇!〇 ;
[0012] (5)、由进气道的极限出口反压、出口喉道面积以及来流总温获得该进气道的最大 出流能力；
[0013] (6)、依据进气道通流时的腔体储气量可知不起动时初始腔体已有50% -60%的 填充；并且此刻的进气道出口基本一直处于极限反压状态，因此取最大出流流量的95%作 为平均出流流量；
[0014] (7)、基于上述步骤（6)的假设计算得到的腔体储气量增益以及进出口流量差，可 计算得到进气道腔体内高压气体的积蓄时间；
[0015] (8)、在口外不起动波系运动阶段中，来流总温对应的滞止声速（〇；)即为该阶段 波系运动至上游的最大速度，由于该值远高于其在离开通道的初始速度，因此选取该滞止 声速的一半为该波系运动阶段的平均速度；
[0016] (9)、由于该振荡阶段的高度非定常性，其振荡幅度覆盖整个进气道前体；进而在 已知该阶段的波系运动距离和运动速度的基础上，可得到口外不起动波系运动时间；其中 波系运动距离为进气道前体压缩面长度1^的两倍；运动速度为来流总温对应滞止声速的一 半；
[0017] (10)、最后将进气道腔体内高压气体积蓄时间和口外不起动波系运动时间叠加， 即可得到对应的进气道不起动振荡周期，从而换算得到相应的振荡频率。
[0018] 本发明的混压式超声速、高超声速进气道不起动振荡频率预测方法的一个关键在 于，通过将进气道内通道结尾激波系的复杂非稳定前传过程转化为对其内腔体储气量变 化，从而可以对其进行定量分析，实现该阶段时长的预测。其基本原理在于：不起动过程中 通道内结尾激波系的前传诱因是进气道的进出口流量不平衡，捕获进入内通道的气流由于 受限于进气道出口的出流能力，不能及时地排出，进而在内通道下游形成低速高压区；随着 时间的推进，该高压区将逐渐向上游蔓延，流场则表现为结尾激波系的前传；为此，可以通 过估算该阶段进气道腔体内储气量的变化值和进出口流量差，得到该阶段的持续时间。本 方法规避了对复杂的内通道激波系运动速度的求解和积分，从气体填充量的角度定量获得 了该阶段的持续时间，且适用于任何形状的进气道内腔体。
【附图说明】
[0019] 图1是进气道处于临界工作状态下的示意图
[0020] 图2是进气道在下游壅塞作用下的高压气体堆积过程示意图
[0021] 图3是进气道在下游壅塞作用下的口外波系运动过程示意图
【具体实施方式】
[0022] 本发明公开了一种混压式超声速、高超声速进气道不起动振荡频率预测方法。
[0023] 请参阅图1至图3所示，结合一种二元混压式高超声速进气道的实例，下面对采用 本发明预测该进气道不起动振荡频率的详细实施步骤进行叙述。
[0024] (1)、首先将一个混压式进气道不起动振荡周期分为两个阶段：进气道内腔体高压 气体的积蓄阶段和口外不起动波系运动阶段。这里需要说明的是，任意的混压式超声速、高 超声速进气道不起动振荡周期均可划分为进气道腔体内高压气体的积蓄和口外不起动波 系运动两个阶段，其原理在于对于该类进气道的不起动诱因均是进出口流量的不平衡，且 其口外波系均存在大振幅的往复运动，因而在振荡中必将出现上述两个阶段。接下来，通过 分别计算两个阶段的时长t2得到不起动振荡周期的长度t。
[0025] (2)、如图1所示，该进气道为一典型的二元混压式进气道，已知来流参数 (M。，VmI；，Pm P (!)以及前体一级压缩面1的压缩角si、二级压缩面18的压缩角S2，来流通 过前体两道斜激波2进入到内通道进口 3,通过斜激波理论可计算得到内通道进口 3上的气 流参数（M3，v3，T3，p3,P3)，并且结合内通道进口 3处的横截面积A3，可计算得到内通道进口 的捕获流量电=P.rvvA;。其中，M、v、T、p、P分别为气流的马赫数、速度、温度、压力和密度， S:为一级压缩面1与水平面的夹角，S2为二级压缩面2与一级压缩面1的夹角。上述计 算中，已知每道斜激波的波前参数时，可以采用以下关系式组计算获得波后参数：
[0026] 由
解得0
[0027] 而后
[0033] (3)、再经过唇罩4处斜激波5 的增压，同样通过斜激波理论可以计算得到斜激波5 后的气流参数（M5,v5,T5,p5,P5)。假设此刻进气道处于临界工作状态，结尾激波6停留在隔 离段7的上游喉道截面8。进而可通过正激波理论得到此刻的波后参数（M8,v8,T8,p8,P8)。上述计算中，正激波波后的参数采用以下关系式组计算获得
[0039] (4)、然后，依据一维定常管流理论以及出口截面9的面积A9、喉道截面8的面积 A8，得到出口处临界状态下的气流参数（M9，T9，p9,p9)。具体如下：
[0043] 进而计算
[0046] 09=1)9八1?.1'9)，其中1? = 287.4 1〇^.1()
[0047] (5)、根据步骤（3)得到的进气道出口临界反压p9，以及已知的气体来流总温 可通过气体状态方程计算出对应的临界储气密度P9。进一步，通过乘上已知的进气道内通 道的体积V。（包括隔离段7以及燃烧室扩张段10在内的腔体体积），即可得到对应的气体 储气阶段最大储气量％=P9XVC。
[0048] (6)、如图2所示，在进气道下游设置堵锥11构成了不起动诱因，将在进气道内通 道下游形成流速极低的高压区13和不断前传的结尾激波14。由于该高压区内的压力分布 均匀，且取值为步骤（3)中得到的出口临界反压&，因此可通该压力值p9、内通道出口喉道 12的流通面积A12获得此时进气道的最大出流能力"、，如下
[0050] (7)、选取内通道腔体的初始气流填充率为50%，进气道出口的平均流量 为95%?岣2。为此，该储气过程中腔体内的气体变化量为50% ?nii，且进出口流量 差为电-95%'办2,进而两者求商可得到该进气道腔体内高压气体积蓄阶段的时长
[0051] (8)、如图3所示，当结尾激波在下游高压推动下越过唇罩前缘15后，将演变为口 外运动波系16向压缩面17上游传去。在高压的推动下，其所能达到的最大速度为来流总 温对应的滞止声速Q；。由于该速度远高于该波系的初始运动速度，因此认为整个运动过程 的平均速度为〇；/2。
[0052] (9)、若该前体压缩面17的长度为U，则整个口外波系的运动行程为2U，再依据步 骤（7)得到的平均运动速度〇；/2,可马上得到该进气道口外不起动波系运动阶段的时长t2 =4IVC:。
[0053] (10)、通过叠加步骤（6)和步骤（8)得到的振荡过程两个阶段的时长&和12，可 得到整个的振荡周期长度t=h+t2，进而换算得到相应的振荡频率f= 1/t。
[0054] 另外，本发明的具体实现方法和途径很多，以上所述仅是本发明的优选实施方式。 该估算方法所适用的混压式超声速、高超声速进气道包括二维形式、轴对称形式、三维侧压 形式、内转形式以及任意组合形式的混压式进气道。应当指出，对于本技术领域的普通技术 人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也 应视为本发明的保护范围。本实施例中未明确的各组成部分均可用现有技术和理论加以实 现。
【主权项】
1. 混压式超声速、高超声速进气道不起动振荡频率预测方法，包括以下步骤： (1) 、将一个混压式进气道不起动振荡周期分为两个阶段：进气道腔体内高压气体积蓄 阶段和口外不起动波系运动阶段； (2) 、在已知具体的来流条件和进气道前体几何参数下，通过已有的斜激波理论得到进 气道进口前的气流参数，并换算得到进气道的捕获流量； (3) 、通过假设进气道处于临界工作状态，在已知进口气体参数和进气道通道面积变化 的基础上，结合正激波理论和一维定常管流理论计算得到该进气道的极限反压；临界工作 状态指结尾激波停于进气道喉道处时的状态； (4) 、由上述极限反压和已知的来流总温换算出最大储气密度，再乘以进气道内通道腔 体容积得到对应气体积蓄阶段的最大储气量（Hl1); (5) 、由进气道的极限出口反压、出口喉道面积以及来流总温获得该进气道的最大出流 能力； (6) 、依据进气道通流时的腔体储气量可知不起动时初始腔体已有50% -60%的填充； 并且此刻的进气道出口基本一直处于极限反压状态，因此取最大出流流量的95%作为平均 出流流量； (7) 、基于上述步骤（6)的假设计算得到的腔体储气量增益以及进出口流量差，可计算 得到进气道腔体内高压气体的积蓄时间； (8) 、在口外不起动波系运动阶段中，来流总温对应的滞止声速（〇；)即为该阶段波系 运动至上游的最大速度，由于该值远高于其在离开通道的初始速度，因此选取该滞止声速 的一半为该波系运动阶段的平均速度； (9) 、由于该振荡阶段的高度非定常性，其振荡幅度覆盖整个进气道前体；进而在已知 该阶段的波系运动距离和运动速度的基础上，可得到口外不起动波系运动时间；其中波系 运动距离为进气道前体压缩面长度1^的两倍；运动速度为来流总温对应滞止声速的一半； (10) 、最后将进气道腔体内高压气体积蓄时间和口外不起动波系运动时间叠加，即可 得到对应的进气道不起动振荡周期，从而换算得到相应的振荡频率。2. 根据权利要求1所述的混压式超声速、高超声速进气道不起动振荡频率估算方法， 其特征在于：步骤（2)中，已知具体的来流条件和进气道前体几何参数（Mtl, ^Tci, P(l，Ptl); 通过斜激波理论可计算得到内通道进口上的气流参数（M3, v3, T3, p3, P3)，并且结合内通道 进口处的横截面积（A3)，计算得到内通道进口的捕获流量电=PWVA.;;其中，m、v、T、p、P 分别为气流的马赫数、速度、温度、压力和密度。3. 根据权利要求2所述的混压式超声速、高超声速进气道不起动振荡频率估算方法， 其特征在于：步骤（3)中，进气道的极限反压（p9)的计算方法为：其中P8为喉道截面处的压力，M8为喉道截面处的马赫数，M9为出口处临界状态下的马 赫数，k取1.4。4. 根据权利要求3所述的混压式超声速、高超声速进气道不起动振荡频率估算方法， 其特征在于：步骤（5)中，进气道的最大出流能力咚2为：A12为内通道出口喉道的流通面积。5. 根据权利要求4所述的混压式超声速、高超声速进气道不起动振荡频率估算方法， 其特征在于：进气道腔体内高压气体的积蓄时间为： tx = 50% · Dil / (m3 - 95%' mr) 〇6. 根据权利要求5所述的混压式超声速、高超声速进气道不起动振荡频率估算方法， 其特征在于：进气道口外不起动波系运动时间t2= 41^/0^。7. 根据权利要求1至5中任一项所述的混压式超声速、高超声速进气道不起动振荡估 算方法，其特征在于：该估算方法所适用的混压式超声速、高超声速进气道包括二维形式、 轴对称形式、三维侧压形式、内转形式以及任意组合形式的混压式进气道。
【专利摘要】本发明公开了一种混压式超声速、高超声速进气道不起动振荡频率预测方法，通过将单个振荡周期划分为进气道内腔体高压气体积蓄和口外不起动波系运动两个阶段，并依次通过将进气道通道内结尾激波系的前传过程转化为对其腔体储气量变化的定量分析，以及将来流总温对应的滞止声速作为不起动波系口外运动最高速度的方法，可对混压式进气道不起动振荡频率进行快速、准确的预估。该方法从振荡机理出发简化了振荡模型，且对进气道的几何参数无特定要求，因此该方法对不同形式的混压式进气道不起动振荡频率估算具有良好的通用性，且同时适用于混压式超声速进气道和高超声速进气道。
【IPC分类】G06F19/00
【公开号】CN104899418
【申请号】CN201510202664
【发明人】谭慧俊, 张启帆, 黄河峡, 陈昊, 孙姝, 宁乐
【申请人】南京航空航天大学
【公开日】2015年9月9日
【申请日】2015年4月24日
【技术领域】
[0001] 本发明涉及混压式超声速、高超声速进气道，尤其是其不起动振荡频率的预测方 法。
【背景技术】
[0002] 超声速、高超声速进气道是吸气式高速推进系统的关键部件之一，其包括外压式、 内压式和混压式三类。其中，混压式进气道由于能在保证较高内流气动性能的同时，减小进 气道的气动阻力、轮廓尺寸和结构重量，故在超声速、高超声速进气道的设计中得到了广泛 应用。然而，由于内压缩的引入，使得混压式进气道容易陷入不起动状态。
[0003] 不起动状态是超声速、高超声速进气道的一类非正常工作状态。通常，若混压式超 声速、高超声速进气道的流量捕获特性未因其内部流态的改变而受到影响，称进气道为起 动，否则为不起动。当进气道处于不起动状态时，不仅总压恢复系数和流量系数急剧下降， 使得推进系统不能产生推力，且极可能伴随着波系以及通道压强的剧烈振荡，带来破坏力 极强的周期性力载荷和热载荷，造成结构破坏，并使飞行器在非定常起动力的作用下边的 更加难以控制。为此，有必要揭示混压式进气道振荡流态的产生机理，建立振荡频率的预测 方法，进而为发动机和飞行器的结构强度设计提供参考。
[0004] 自Oswatitsch于1944年首次观测到超声速进气道的振荡流态（即喘振）以来， 国内外有大量的研宄工作致力于这一现象的描述、预测和控制。通过对外压式超声速进气 道喘振的仿真和实验研宄，各国学者发现：扰动的声波反馈在喘振信号传递环中起着重要 作用，故喘振的主特征频率以及各次要特征频率往往与进气道腔体的各阶声学谐振频率相 关。然而，混压式进气道与外压式进气道的不起动流态存在显著差异，前者为超声速溢流， 后者为亚声速溢流，使得其不稳定机理和振荡特性都有所不同。为此，若依旧采用腔体声学 振荡模型对混压式进气道的不起动振荡频率进行估算，必然会带来较大的误差。并且，近年 来混压式进气道在飞行器上的应用趋于广泛，并已衍生发展出了多种不同的具体形式，如 二维形式、轴对称形式、三维侧压形式、内乘波形式以及组合形式等。因此，发展适用于混压 式超声速、高超声速进气道不起动振荡频率的估算方法显得十分重要。
【发明内容】
[0005] 本发明提供一种适用于混压式超声速、高超声速进气道的不起动振荡频率预测方 法，能够实现对不同来流状态和几何参数下进气道不起动振荡频率的快速、准确预测，并且 适用于各种形式的混压式进气道。
[0006] 为达到上述目的，本发明的混压式超声速、高超声速进气道不起动振荡频率预测 方法可采用如下技术方案：
[0007] 混压式超声速、高超声速进气道不起动振荡频率估算方法包括以下步骤：
[0008] (1)、将一个混压式进气道不起动振荡周期分为两个阶段：进气道腔体内高压气体 积蓄阶段和口外不起动波系运动阶段；
[0009] (2)、在已知具体的来流条件和进气道前体几何参数下，通过已有的斜激波理论得 到进气道进口前的气流参数，并换算得到进气道的捕获流量；
[0010] (3)、通过假设进气道处于临界工作状态，在已知进口气体参数和进气道通道面积 变化的基础上，结合正激波理论和一维定常管流理论计算得到该进气道的极限反压；临界 工作状态指结尾激波停于进气道喉道处时的状态；
[0011] (4)、由上述极限反压和已知的来流总温换算出最大储气密度，再乘以进气道内通 道腔体容积得到对应气体积蓄阶段的最大储气量〇!〇 ;
[0012] (5)、由进气道的极限出口反压、出口喉道面积以及来流总温获得该进气道的最大 出流能力；
[0013] (6)、依据进气道通流时的腔体储气量可知不起动时初始腔体已有50% -60%的 填充；并且此刻的进气道出口基本一直处于极限反压状态，因此取最大出流流量的95%作 为平均出流流量；
[0014] (7)、基于上述步骤（6)的假设计算得到的腔体储气量增益以及进出口流量差，可 计算得到进气道腔体内高压气体的积蓄时间；
[0015] (8)、在口外不起动波系运动阶段中，来流总温对应的滞止声速（〇；)即为该阶段 波系运动至上游的最大速度，由于该值远高于其在离开通道的初始速度，因此选取该滞止 声速的一半为该波系运动阶段的平均速度；
[0016] (9)、由于该振荡阶段的高度非定常性，其振荡幅度覆盖整个进气道前体；进而在 已知该阶段的波系运动距离和运动速度的基础上，可得到口外不起动波系运动时间；其中 波系运动距离为进气道前体压缩面长度1^的两倍；运动速度为来流总温对应滞止声速的一 半；
[0017] (10)、最后将进气道腔体内高压气体积蓄时间和口外不起动波系运动时间叠加， 即可得到对应的进气道不起动振荡周期，从而换算得到相应的振荡频率。
[0018] 本发明的混压式超声速、高超声速进气道不起动振荡频率预测方法的一个关键在 于，通过将进气道内通道结尾激波系的复杂非稳定前传过程转化为对其内腔体储气量变 化，从而可以对其进行定量分析，实现该阶段时长的预测。其基本原理在于：不起动过程中 通道内结尾激波系的前传诱因是进气道的进出口流量不平衡，捕获进入内通道的气流由于 受限于进气道出口的出流能力，不能及时地排出，进而在内通道下游形成低速高压区；随着 时间的推进，该高压区将逐渐向上游蔓延，流场则表现为结尾激波系的前传；为此，可以通 过估算该阶段进气道腔体内储气量的变化值和进出口流量差，得到该阶段的持续时间。本 方法规避了对复杂的内通道激波系运动速度的求解和积分，从气体填充量的角度定量获得 了该阶段的持续时间，且适用于任何形状的进气道内腔体。
【附图说明】
[0019] 图1是进气道处于临界工作状态下的示意图
[0020] 图2是进气道在下游壅塞作用下的高压气体堆积过程示意图
[0021] 图3是进气道在下游壅塞作用下的口外波系运动过程示意图
【具体实施方式】
[0022] 本发明公开了一种混压式超声速、高超声速进气道不起动振荡频率预测方法。
[0023] 请参阅图1至图3所示，结合一种二元混压式高超声速进气道的实例，下面对采用 本发明预测该进气道不起动振荡频率的详细实施步骤进行叙述。
[0024] (1)、首先将一个混压式进气道不起动振荡周期分为两个阶段：进气道内腔体高压 气体的积蓄阶段和口外不起动波系运动阶段。这里需要说明的是，任意的混压式超声速、高 超声速进气道不起动振荡周期均可划分为进气道腔体内高压气体的积蓄和口外不起动波 系运动两个阶段，其原理在于对于该类进气道的不起动诱因均是进出口流量的不平衡，且 其口外波系均存在大振幅的往复运动，因而在振荡中必将出现上述两个阶段。接下来，通过 分别计算两个阶段的时长t2得到不起动振荡周期的长度t。
[0025] (2)、如图1所示，该进气道为一典型的二元混压式进气道，已知来流参数 (M。，VmI；，Pm P (!)以及前体一级压缩面1的压缩角si、二级压缩面18的压缩角S2，来流通 过前体两道斜激波2进入到内通道进口 3,通过斜激波理论可计算得到内通道进口 3上的气 流参数（M3，v3，T3，p3,P3)，并且结合内通道进口 3处的横截面积A3，可计算得到内通道进口 的捕获流量电=P.rvvA;。其中，M、v、T、p、P分别为气流的马赫数、速度、温度、压力和密度， S:为一级压缩面1与水平面的夹角，S2为二级压缩面2与一级压缩面1的夹角。上述计 算中，已知每道斜激波的波前参数时，可以采用以下关系式组计算获得波后参数：
[0026] 由
解得0
[0027] 而后
[0033] (3)、再经过唇罩4处斜激波5 的增压，同样通过斜激波理论可以计算得到斜激波5 后的气流参数（M5,v5,T5,p5,P5)。假设此刻进气道处于临界工作状态，结尾激波6停留在隔 离段7的上游喉道截面8。进而可通过正激波理论得到此刻的波后参数（M8,v8,T8,p8,P8)。上述计算中，正激波波后的参数采用以下关系式组计算获得
[0039] (4)、然后，依据一维定常管流理论以及出口截面9的面积A9、喉道截面8的面积 A8，得到出口处临界状态下的气流参数（M9，T9，p9,p9)。具体如下：
[0043] 进而计算
[0046] 09=1)9八1?.1'9)，其中1? = 287.4 1〇^.1()
[0047] (5)、根据步骤（3)得到的进气道出口临界反压p9，以及已知的气体来流总温 可通过气体状态方程计算出对应的临界储气密度P9。进一步，通过乘上已知的进气道内通 道的体积V。（包括隔离段7以及燃烧室扩张段10在内的腔体体积），即可得到对应的气体 储气阶段最大储气量％=P9XVC。
[0048] (6)、如图2所示，在进气道下游设置堵锥11构成了不起动诱因，将在进气道内通 道下游形成流速极低的高压区13和不断前传的结尾激波14。由于该高压区内的压力分布 均匀，且取值为步骤（3)中得到的出口临界反压&，因此可通该压力值p9、内通道出口喉道 12的流通面积A12获得此时进气道的最大出流能力"、，如下
[0050] (7)、选取内通道腔体的初始气流填充率为50%，进气道出口的平均流量 为95%?岣2。为此，该储气过程中腔体内的气体变化量为50% ?nii，且进出口流量 差为电-95%'办2,进而两者求商可得到该进气道腔体内高压气体积蓄阶段的时长
[0051] (8)、如图3所示，当结尾激波在下游高压推动下越过唇罩前缘15后，将演变为口 外运动波系16向压缩面17上游传去。在高压的推动下，其所能达到的最大速度为来流总 温对应的滞止声速Q；。由于该速度远高于该波系的初始运动速度，因此认为整个运动过程 的平均速度为〇；/2。
[0052] (9)、若该前体压缩面17的长度为U，则整个口外波系的运动行程为2U，再依据步 骤（7)得到的平均运动速度〇；/2,可马上得到该进气道口外不起动波系运动阶段的时长t2 =4IVC:。
[0053] (10)、通过叠加步骤（6)和步骤（8)得到的振荡过程两个阶段的时长&和12，可 得到整个的振荡周期长度t=h+t2，进而换算得到相应的振荡频率f= 1/t。
[0054] 另外，本发明的具体实现方法和途径很多，以上所述仅是本发明的优选实施方式。 该估算方法所适用的混压式超声速、高超声速进气道包括二维形式、轴对称形式、三维侧压 形式、内转形式以及任意组合形式的混压式进气道。应当指出，对于本技术领域的普通技术 人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也 应视为本发明的保护范围。本实施例中未明确的各组成部分均可用现有技术和理论加以实 现。
【主权项】
1. 混压式超声速、高超声速进气道不起动振荡频率预测方法，包括以下步骤： (1) 、将一个混压式进气道不起动振荡周期分为两个阶段：进气道腔体内高压气体积蓄 阶段和口外不起动波系运动阶段； (2) 、在已知具体的来流条件和进气道前体几何参数下，通过已有的斜激波理论得到进 气道进口前的气流参数，并换算得到进气道的捕获流量； (3) 、通过假设进气道处于临界工作状态，在已知进口气体参数和进气道通道面积变化 的基础上，结合正激波理论和一维定常管流理论计算得到该进气道的极限反压；临界工作 状态指结尾激波停于进气道喉道处时的状态； (4) 、由上述极限反压和已知的来流总温换算出最大储气密度，再乘以进气道内通道腔 体容积得到对应气体积蓄阶段的最大储气量（Hl1); (5) 、由进气道的极限出口反压、出口喉道面积以及来流总温获得该进气道的最大出流 能力； (6) 、依据进气道通流时的腔体储气量可知不起动时初始腔体已有50% -60%的填充； 并且此刻的进气道出口基本一直处于极限反压状态，因此取最大出流流量的95%作为平均 出流流量； (7) 、基于上述步骤（6)的假设计算得到的腔体储气量增益以及进出口流量差，可计算 得到进气道腔体内高压气体的积蓄时间； (8) 、在口外不起动波系运动阶段中，来流总温对应的滞止声速（〇；)即为该阶段波系 运动至上游的最大速度，由于该值远高于其在离开通道的初始速度，因此选取该滞止声速 的一半为该波系运动阶段的平均速度； (9) 、由于该振荡阶段的高度非定常性，其振荡幅度覆盖整个进气道前体；进而在已知 该阶段的波系运动距离和运动速度的基础上，可得到口外不起动波系运动时间；其中波系 运动距离为进气道前体压缩面长度1^的两倍；运动速度为来流总温对应滞止声速的一半； (10) 、最后将进气道腔体内高压气体积蓄时间和口外不起动波系运动时间叠加，即可 得到对应的进气道不起动振荡周期，从而换算得到相应的振荡频率。2. 根据权利要求1所述的混压式超声速、高超声速进气道不起动振荡频率估算方法， 其特征在于：步骤（2)中，已知具体的来流条件和进气道前体几何参数（Mtl, ^Tci, P(l，Ptl); 通过斜激波理论可计算得到内通道进口上的气流参数（M3, v3, T3, p3, P3)，并且结合内通道 进口处的横截面积（A3)，计算得到内通道进口的捕获流量电=PWVA.;;其中，m、v、T、p、P 分别为气流的马赫数、速度、温度、压力和密度。3. 根据权利要求2所述的混压式超声速、高超声速进气道不起动振荡频率估算方法， 其特征在于：步骤（3)中，进气道的极限反压（p9)的计算方法为：其中P8为喉道截面处的压力，M8为喉道截面处的马赫数，M9为出口处临界状态下的马 赫数，k取1.4。4. 根据权利要求3所述的混压式超声速、高超声速进气道不起动振荡频率估算方法， 其特征在于：步骤（5)中，进气道的最大出流能力咚2为：A12为内通道出口喉道的流通面积。5. 根据权利要求4所述的混压式超声速、高超声速进气道不起动振荡频率估算方法， 其特征在于：进气道腔体内高压气体的积蓄时间为： tx = 50% · Dil / (m3 - 95%' mr) 〇6. 根据权利要求5所述的混压式超声速、高超声速进气道不起动振荡频率估算方法， 其特征在于：进气道口外不起动波系运动时间t2= 41^/0^。7. 根据权利要求1至5中任一项所述的混压式超声速、高超声速进气道不起动振荡估 算方法，其特征在于：该估算方法所适用的混压式超声速、高超声速进气道包括二维形式、 轴对称形式、三维侧压形式、内转形式以及任意组合形式的混压式进气道。
【专利摘要】本发明公开了一种混压式超声速、高超声速进气道不起动振荡频率预测方法，通过将单个振荡周期划分为进气道内腔体高压气体积蓄和口外不起动波系运动两个阶段，并依次通过将进气道通道内结尾激波系的前传过程转化为对其腔体储气量变化的定量分析，以及将来流总温对应的滞止声速作为不起动波系口外运动最高速度的方法，可对混压式进气道不起动振荡频率进行快速、准确的预估。该方法从振荡机理出发简化了振荡模型，且对进气道的几何参数无特定要求，因此该方法对不同形式的混压式进气道不起动振荡频率估算具有良好的通用性，且同时适用于混压式超声速进气道和高超声速进气道。
【IPC分类】G06F19/00
【公开号】CN104899418
【申请号】CN201510202664
【发明人】谭慧俊, 张启帆, 黄河峡, 陈昊, 孙姝, 宁乐
【申请人】南京航空航天大学
【公开日】2015年9月9日
【申请日】2015年4月24日
最新回复(0)