一种针对超声速民机的旋成体静音锥快速设计方法

1.本发明属于超声速民机声爆抑制技术领域，具体涉及一种针对超声速民机的旋成体静音锥快速设计方法。


背景技术：

2.声爆是制约超声速民机投入商业运营的瓶颈。为解决声爆问题，国内外研究人员从上世纪60年代开展大量的研究。其中，本世纪初美国湾流公司howe提出机头静音锥技术(专利号us6698684)，是一种行之有效的声爆抑制技术。机头静音锥技术是通过在超声速民机头部添加多段细长杆组合而成的装置，可实现声爆抑制，其抑制声爆的机理为：利用这些细长杆产生的多道弱激波代替传统飞机头部产生的一道强激波，且在传播到地面的过程中不发生汇聚。这样不仅降低地面波形机头激波的强度，还能够大幅增加头激波的上升时间。
3.美国nasa(national aeronautics and space administration)于2002年在4
×
4整体式风洞中对howe设计的静音锥开展风洞试验，并将测量结果与cfd(computational fluid dynamic)计算得到的声爆近场波形进行对比，从而验证数值模拟计算结果的有效性。2006年又将可伸缩的机头静音锥技术应用到f-15b飞机上进行飞行试验，从而验证机头静音锥技术的有效性。
4.然而，传统方式中，设计有效降低声爆强度的静音锥往往采用高可信度计算与优化设计相结合的手段，耗时十分长，且严重依赖于计算条件。


技术实现要素：

5.针对现有技术存在的缺陷，本发明提供一种针对超声速民机的旋成体静音锥快速设计方法，可有效解决上述问题。
6.本发明采用的技术方案如下：
7.本发明提供一种针对超声速民机的旋成体静音锥快速设计方法，包括以下步骤：
8.步骤1，选择超声速民机基准构型；所述超声速民机基准构型的机身长度为l；
9.步骤2，在设计马赫数ma和设计攻角α的巡航条件下，对超声速民机基准构型进行数值模拟，提取得到超声速民机基准构型正下方距离为h位置的超声速巡航状态下的基准声爆近场波形；
10.其中，在基准声爆近场波形中，基准声爆近场波形起点为o
start
，基准声爆近场波形终点为o
end
，机身头部(s1)产生的激波在基准声爆近场波形中对应位置为点b1；
11.步骤3，确定激波最优峰值p
s_lb
和膨胀波最优长度l
e_lb
：
12.步骤3.1，自超声速民机基准构型的机身头部(s1)开始，沿飞机轴向查找到满足以下条件的机身截面(s2)：机身截面(s2)的直径d0＝0.5％l～2.5％l；
13.步骤3.2，确定机身截面(s2)产生的激波在步骤2得到的基准声爆近场波形中的对应位置，表示为点b2；
14.步骤3.3，在步骤2得到的基准声爆近场波形中，截去起点o
start
到点b1以及点b1到
点b2之间的基准声爆近场波形前段b1b2；保留点b2到终点o
end
之间的基准声爆近场波形后段b2o
end
；
15.步骤3.4，在起点o
start
和点b2之间，依次选择三个点，分别为：点a1、点a2和点a3；
16.确定点a1、点a2和点a3的初始位置后，采用下面优化设计过程，不断调节点a1、点a2和点a3的位置，使点a1到点a2的连接直线形成的激波a1a2的峰值ps，为激波最优峰值p
s_lb
；使点a2到点a3的连接直线形成的膨胀波a2a3的长度le，为膨胀波最优长度l
e_lb
；
17.其中，激波a1a2的峰值ps，为点a1和点a2之间的高度差；膨胀波a2a3的长度le，为点a2到点a3的水平距离；
18.优化设计过程为：
19.步骤3.4.1，确定约束条件为：
20.条件1：激波a1a2的峰值ps＝0.1*基准声爆近场波形最大超压值～1.5*基准声爆近场波形最大超压值；
21.条件2：激波a1a2的斜率＝0.008～0.012；
22.条件3：膨胀波a2a3的长度le＝0.1l～0.5l；
23.条件4：点a3到点a1连线的斜率＝-0.00132～-0.00124；
24.条件5：激波a1a2的斜率＝点b2到点a3连线的斜率；
25.步骤3.4.2，在满足步骤3.4.1确定的约束条件下，给定点a1、点a2和点a3的初始位置，并在基准声爆近场波形后段b2o
end
的前方，依次连接起点o
start
、点a1、点a2、点a3和点b2，从而得到添加激波和膨胀波的声爆近场波形；
26.结合声爆远场传播方法，以添加激波和膨胀波的声爆近场波形作为输入，使其传播到远场，并求解得到远场声爆强度；
27.判断得到的远场声爆强度是否达到优化目标最低值；如果满足，则停止优化，点a1、点a2和点a3的当前位置，即为最优位置，进而得到激波最优峰值p
s_lb
和膨胀波最优长度l
e_lb
；如果不满足，则改变a1、点a2和点a3的位置，实质为改变激波a1a2的峰值ps和膨胀波a2a3的长度le，重复进行迭代优化；
28.步骤4，对旋成体静音锥外形进行优化设计，得到最优旋成体静音锥外形：
29.步骤4.1，设定旋成体静音锥外形的设计变量为：圆锥段长度l1、圆柱段长度l2和圆柱段半径r；
30.步骤4.2，在满足以下约束条件下，给定圆锥段长度l1、圆柱段长度l2和圆柱段半径r的初始值，从而确定初始旋成体静音锥外形；
31.s.t.l1∈[0.01l,0.08l]
[0032]
l2∈[0.01l,0.08l]
[0033]
r∈[0.05l,0.25l]
[0034]
步骤4.3，根据修正的声爆线化理论，获得初始旋成体静音锥外形的体积等效截面积分布，并对初始旋成体静音锥外形的声爆近场波形进行预测，获得初始旋成体静音锥外形的声爆近场波形；
[0035]
从初始旋成体静音锥外形的声爆近场波形中，提取到初始旋成体静音锥外形对应的激波峰值ps和膨胀波长度le；
[0036]
步骤4.4，判断步骤4.3得到的激波峰值ps和膨胀波长度le是否满足以下关系式：
[0037]
min f＝w1·
|p
s_lb-ps| w2·
|l
e_lb-le|
[0038]
其中：
[0039]
w1和w2为权重系数，且满足w1 w2＝1.0；
[0040]
如果满足，则步骤4.2给定的圆锥段长度l1、圆柱段长度l2和圆柱段半径r，即为最终优化得到的旋成体静音锥外形最优参数，进而得到最优旋成体静音锥外形，然后执行步骤5；如果不满足，在满足约束条件下，调整圆锥段长度l1、圆柱段长度l2和圆柱段半径r，返回步骤4.3；
[0041]
步骤5，将最优旋成体静音锥外形的几何参数表示为：圆锥段最优长度l1(best)、圆柱段最优长度l2(best)和圆柱段最优半径r(best)；
[0042]
自超声速民机基准构型的机身头部(s1)开始，沿飞机轴向查找到半径为圆柱段最优半径r(best)的机身截面(s3)，截去机身头部(s1)到机身截面(s3)的一段机身；然后，将最优旋成体静音锥外形添加到超声速民机基准构型的机身截面(s3)的前面，得到超声速民机旋成体静音锥构型。
[0043]
优选的，步骤4.3具体为：
[0044]
步骤4.3.1，给定初始旋成体静音锥外形的几何参数为：圆锥段长度l1、圆柱段长度l2和圆柱段半径r；
[0045]
采用下式，获得初始旋成体静音锥外形的体积等效截面积分布sv(x)：
[0046][0047]
其中：
[0048]
x
11
＝l
1-r/tan(μ-α)
[0049]
x
12
＝l1 r/tan(μ-α)
[0050]
x
13
＝l1 l
2-r/tan(μ-α)
[0051]f11
(x)＝(x-l1)
·
tan(μ-α)
[0052]f12
(x)＝2arccos(f
11
(x))
[0053]f13
(x)＝π-arccos(f
11
(x))
[0054]
其中：
[0055]
α：设计攻角；
[0056]
μ：步骤2的设计马赫数ma对应的马赫角；
[0057]
β：圆锥段锥角；并且满足：β＜μ-α；
[0058]
x
11
：旋成体静音锥第一轴向位置；
[0059]
x
12
：旋成体静音锥第二轴向位置；
[0060]
x
13
：旋成体静音锥第三轴向位置；
[0061]
x
11
、x
12
和x
13
的确定方法为：静音锥圆锥段与静音锥圆柱段的交界面为圆锥-圆柱交界面，圆锥-圆柱交界面的下轮廓与飞机对称面相交于圆锥-圆柱交界面下端点p1，圆锥-圆柱交界面的上轮廓与飞机对称面相交于圆锥-圆柱交界面上端点p2，静音锥圆柱段下轮
廓与飞机对称面相交于圆柱-机身交界面下端点p3；
[0062]
过圆锥-圆柱交界面下端点p1，作与马赫平面平行的平面，与旋成体静音锥轴线相交于旋成体静音锥第一轴向位置x
11
；
[0063]
过圆锥-圆柱交界面上端点p2，作与马赫平面平行的平面，与旋成体静音锥轴线相交于旋成体静音锥第二轴向位置x
12
；
[0064]
过圆柱-机身交界面下端点p3，作与马赫平面平行的平面，与旋成体静音锥轴线相交于旋成体静音锥第三轴向位置x
13
；
[0065]f11
(x)：第一中间量参数；
[0066]f12
(x)：第二中间量参数；
[0067]f13
(x)：第三中间量参数；
[0068]
x：以初始旋成体静音锥外形的锥尖为原点，过旋成体静音锥轴线作x轴，旋成体静音锥轴线的不同位置点的x轴坐标；
[0069]
步骤4.3.2，根据初始旋成体静音锥外形的体积等效截面积分布sv(x)，计算声爆f函数f(τ)：
[0070][0071]
其中：
[0072]
τ：在初始旋成体静音锥外形为巡航条件下时，在其正下方距离为h位置建立与初始旋成体静音锥外形轴线平行的τ轴，初始旋成体静音锥外形的锥尖在τ轴的投影为τ轴原点；τ轴的不同位置点具有τ轴坐标；
[0073]
步骤4.3.3，采用下式，计算初始旋成体静音锥外形的线化的声爆近场波形：
[0074][0075]
其中：
[0076]
p
inf
：自由来流压强；
[0077]
dp：超声速民机扰动压强与自由来流压强差；
[0078]
γ：大气比热比；
[0079]
b：普朗特-格劳厄特系数，计算公式为：
[0080]
步骤4.3.4，使用面积平衡法或burgers-hayers方法对步骤4.3.3得到的初始旋成体静音锥外形的线化的声爆近场波形进行非线性修正，得到最终的初始旋成体静音锥外形的声爆近场波形。
[0081]
本发明提供的一种针对超声速民机的旋成体静音锥快速设计方法具有以下优点：
[0082]
本发明采用修正的声爆线化理论和优化设计方法，考虑到修正的声爆线化理论预测近场波形简便，以及优化设计过程中设计变量数少的优点，能够直接以远场声爆强度为目标，快速的设计出添加旋成体静音锥的超声速民机外形，有效抑制超声速民机巡航阶段产生严重声爆的问题，更好地满足超声速民机飞行时对声爆强度的要求。
附图说明
[0083]
图1为本发明提供的针对超声速民机的旋成体静音锥快速设计方法的流程示意图；
[0084]
图2为基准声爆近场波形提取过程的原理图；
[0085]
图3为超声速民机基准构型以及机身截面s2位置的示意图；
[0086]
图4为基准声爆近场波形与添加激波和膨胀波的声爆近场波形示意图；
[0087]
图5为优化设计激波最优峰值和膨胀波最优长度过程中远场声爆强度随添加样本点数量的收敛示意图；
[0088]
图6为使用马赫平面截旋成体静音锥的示意图；
[0089]
图7为超声速民机旋成体静音锥构型的外形示意图；
[0090]
图8为选取的超声速民机基准构型几何外形图；
[0091]
图9为本发明最终设计得到的超声速民机静音锥构型图；
[0092]
图10为静音锥构型与超声速民机基准构型的声爆近场波形对比示意图；
[0093]
图11为静音锥构型与超声速民机基准构型的声爆远场波形对比示意图。
具体实施方式
[0094]
为了使本发明所解决的技术问题、技术方案及有益效果更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
[0095]
本发明提供一种针对超声速民机的旋成体静音锥快速设计方法，在降低声爆强度的条件下，有效解决静音锥技术设计过程中耗时长的问题。具体的，本发明采用修正的声爆线化理论和优化设计方法，考虑到修正的声爆线化理论预测近场波形简便，以及优化设计过程中设计变量数少的优点，能够直接以远场声爆强度为目标，快速的设计出添加旋成体静音锥的超声速民机外形，有效抑制超声速民机巡航阶段产生严重声爆的问题，更好地满足超声速民机飞行时对声爆强度的要求。
[0096]
参考图1，本发明提供一种针对超声速民机的旋成体静音锥快速设计方法，包括以下步骤：
[0097]
步骤1，选择超声速民机基准构型；所述超声速民机基准构型的机身长度为l；
[0098]
步骤2，在设计马赫数ma和设计攻角α的巡航条件下，对超声速民机基准构型进行数值模拟，提取得到超声速民机基准构型正下方距离为h位置的超声速巡航状态下的基准声爆近场波形；
[0099]
如图2所示，为基准声爆近场波形提取过程的原理图；在图2中，在基准声爆近场波形中，基准声爆近场波形起点为o
start
，基准声爆近场波形终点为o
end
，机身头部s1产生的激波在基准声爆近场波形中对应位置为点b1；
[0100]
步骤3，确定激波最优峰值p
s_lb
和膨胀波最优长度l
e_lb
：
[0101]
本步骤主要思路为：结合声爆远场传播和优化设计方法，将添加的激波峰值和膨胀波长度作为设计变量，以远场声爆强度为目标，对这两个参数进行优化设计，获得声爆强度最小时对应的激波峰值和膨胀波长度，即为激波最优峰值p
s_lb
和膨胀波最优长度l
e_lb
。其中，声爆远场传播可以采用基于波形参数法的低可信度预测方法或基于广义burgers方程
的高可信度预测方法。
[0102]
步骤3.1，自超声速民机基准构型的机身头部s1开始，沿飞机轴向查找到满足以下条件的机身截面s2：机身截面s2的直径d0＝0.5％l～2.5％l；
[0103]
如图3所示，为超声速民机基准构型以及机身截面s2位置的示意图；在图3中，1、2、3和4，分别代表：机身、机翼、发动机和尾翼。其中，机身长度，是指从机身头部s1到尾部的长度。
[0104]
步骤3.2，确定机身截面s2产生的激波在步骤2得到的基准声爆近场波形中的对应位置，表示为点b2；
[0105]
步骤3.3，在步骤2得到的基准声爆近场波形中，截去起点o
start
到点b1以及点b1到点b2之间的基准声爆近场波形前段b1b2；保留点b2到终点o
end
之间的基准声爆近场波形后段b2o
end
；具体原理，如图4所示，为基准声爆近场波形的放大原理图。
[0106]
步骤3.4，在起点o
start
和点b2之间，依次选择三个点，分别为：点a1、点a2和点a3；
[0107]
确定点a1、点a2和点a3的初始位置后，采用下面优化设计过程，不断调节点a1、点a2和点a3的位置，使点a1到点a2的连接直线形成的激波a1a2的峰值ps，为激波最优峰值p
s_lb
；使点a2到点a3的连接直线形成的膨胀波a2a3的长度le，为膨胀波最优长度l
e_lb
；
[0108]
其中，激波a1a2的峰值ps，为点a1和点a2之间的高度差；膨胀波a2a3的长度le，为点a2到点a3的水平距离；
[0109]
优化设计过程为：
[0110]
步骤3.4.1，确定约束条件为：
[0111]
条件1：激波a1a2的峰值ps＝0.1*基准声爆近场波形最大超压值～1.5*基准声爆近场波形最大超压值；
[0112]
条件2：激波a1a2的斜率＝0.008～0.012；
[0113]
条件3：膨胀波a2a3的长度le＝0.1l～0.5l；
[0114]
条件4：点a3到点a1连线的斜率＝-0.00132～-0.00124；
[0115]
条件5：激波a1a2的斜率＝点b2到点a3连线的斜率；
[0116]
步骤3.4.2，在满足步骤3.4.1确定的约束条件下，给定点a1、点a2和点a3的初始位置，并在基准声爆近场波形后段b2o
end
的前方，依次连接起点o
start
、点a1、点a2、点a3和点b2，从而得到添加激波和膨胀波的声爆近场波形；
[0117]
结合声爆远场传播方法，以添加激波和膨胀波的声爆近场波形作为输入，使其传播到远场，并求解得到远场声爆强度；
[0118]
判断得到的远场声爆强度是否达到优化目标最低值；如果满足，则停止优化，点a1、点a2和点a3的当前位置，即为最优位置，进而得到激波最优峰值p
s_lb
和膨胀波最优长度l
e_lb
；如果不满足，则改变a1、点a2和点a3的位置，实质为改变激波a1a2的峰值ps和膨胀波a2a3的长度le，重复进行迭代优化；
[0119]
如图5所示，为优化过程中远场声爆强度随添加样本点数量的收敛示意图。
[0120]
步骤4，对旋成体静音锥外形进行优化设计，得到最优旋成体静音锥外形：
[0121]
步骤4.1，设定旋成体静音锥外形的设计变量为：圆锥段长度l1、圆柱段长度l2和圆柱段半径r；
[0122]
步骤4.2，在满足以下约束条件下，给定圆锥段长度l1、圆柱段长度l2和圆柱段半径
r的初始值，从而确定初始旋成体静音锥外形；
[0123]
s.t.l1∈[0.01l,0.08l]
[0124]
l2∈[0.01l,0.08l]
[0125]
r∈[0.05l,0.25l]
[0126]
步骤4.3，根据修正的声爆线化理论，获得初始旋成体静音锥外形的体积等效截面积分布，并对初始旋成体静音锥外形的声爆近场波形进行预测，获得初始旋成体静音锥外形的声爆近场波形；
[0127]
从初始旋成体静音锥外形的声爆近场波形中，提取到初始旋成体静音锥外形对应的激波峰值ps和膨胀波长度le；
[0128]
本步骤中，具体采用以下方法，获得最终的初始旋成体静音锥外形的声爆近场波形：
[0129]
步骤4.3.1，给定初始旋成体静音锥外形的几何参数为：圆锥段长度l1、圆柱段长度l2和圆柱段半径r；
[0130]
采用下式，结合图6，获得初始旋成体静音锥外形的体积等效截面积分布sv(x)：其中，图6为使用马赫平面截旋成体静音锥的示意图；在图6中，5代表静音锥圆锥段；6代表静音锥圆柱段。
[0131][0132]
其中：
[0133]
x
11
＝l
1-r/tan(μ-α)
[0134]
x
12
＝l1 r/tan(μ-α)
[0135]
x
13
＝l1 l
2-r/tan(μ-α)
[0136]f11
(x)＝(x-l1)
·
tan(μ-α)
[0137]f12
(x)＝2arccos(f
11
(x))
[0138]f13
(x)＝π-arccos(f
11
(x))
[0139]
其中：
[0140]
α：设计攻角；
[0141]
μ：步骤2的设计马赫数ma对应的马赫角；
[0142]
β：圆锥段锥角；并且满足：β＜μ-α；
[0143]
x
11
：旋成体静音锥第一轴向位置；
[0144]
x
12
：旋成体静音锥第二轴向位置；
[0145]
x
13
：旋成体静音锥第三轴向位置；
[0146]
x
11
、x
12
和x
13
的确定方法为：静音锥圆锥段与静音锥圆柱段的交界面为圆锥-圆柱交界面，圆锥-圆柱交界面的下轮廓与飞机对称面相交于圆锥-圆柱交界面下端点p1，圆锥-圆柱交界面的上轮廓与飞机对称面相交于圆锥-圆柱交界面上端点p2，静音锥圆柱段下轮廓与飞机对称面相交于圆柱-机身交界面下端点p3；
[0147]
过圆锥-圆柱交界面下端点p1，作与马赫平面平行的平面，与旋成体静音锥轴线相交于旋成体静音锥第一轴向位置x
11
；
[0148]
过圆锥-圆柱交界面上端点p2，作与马赫平面平行的平面，与旋成体静音锥轴线相交于旋成体静音锥第二轴向位置x
12
；
[0149]
过圆柱-机身交界面下端点p3，作与马赫平面平行的平面，与旋成体静音锥轴线相交于旋成体静音锥第三轴向位置x
13
；
[0150]f11
(x)：第一中间量参数；
[0151]f12
(x)：第二中间量参数；
[0152]f13
(x)：第三中间量参数；
[0153]
x：以初始旋成体静音锥外形的锥尖为原点，过旋成体静音锥轴线作x轴，旋成体静音锥轴线的不同位置点的x轴坐标；
[0154]
步骤4.3.2，根据初始旋成体静音锥外形的体积等效截面积分布sv(x)，计算声爆f函数f(τ)：
[0155][0156]
其中：
[0157]
τ：在初始旋成体静音锥外形为巡航条件下时，在其正下方距离为h位置建立与初始旋成体静音锥外形轴线平行的τ轴，初始旋成体静音锥外形的锥尖在τ轴的投影为τ轴原点；τ轴的不同位置点具有τ轴坐标；
[0158]
步骤4.3.3，采用下式，计算初始旋成体静音锥外形的线化的声爆近场波形：
[0159][0160]
其中：
[0161]
p
inf
：自由来流压强；
[0162]
dp：超声速民机扰动压强与自由来流压强差；
[0163]
γ：大气比热比；
[0164]
b：普朗特-格劳厄特系数，计算公式为：
[0165]
步骤4.3.4，使用面积平衡法或burgers-hayers方法对步骤4.3.3得到的初始旋成体静音锥外形的线化的声爆近场波形进行非线性修正，得到最终的初始旋成体静音锥外形的声爆近场波形。
[0166]
步骤4.4，判断步骤4.3得到的激波峰值ps和膨胀波长度le是否满足以下关系式：
[0167]
min f＝w1·
|p
s_lb-ps| w2·
|l
e_lb-le|
[0168]
其中：
[0169]
w1和w2为权重系数，且满足w1 w2＝1.0；
[0170]
如果满足，则步骤4.2给定的圆锥段长度l1、圆柱段长度l2和圆柱段半径r，即为最终优化得到的旋成体静音锥外形最优参数，进而得到最优旋成体静音锥外形，然后执行步骤5；如果不满足，在满足约束条件下，调整圆锥段长度l1、圆柱段长度l2和圆柱段半径r，返
回步骤4.3；
[0171]
步骤5，将最优旋成体静音锥外形的几何参数表示为：圆锥段最优长度l1(best)、圆柱段最优长度l2(best)和圆柱段最优半径r(best)；
[0172]
自超声速民机基准构型的机身头部s1开始，沿飞机轴向查找到半径为圆柱段最优半径r(best)的机身截面s3，截去机身头部s1到机身截面s3的一段机身；然后，将最优旋成体静音锥外形添加到超声速民机基准构型的机身截面s3的前面，得到超声速民机旋成体静音锥构型。如图7所示，为超声速民机旋成体静音锥构型的外形示意图。
[0173]
本发明使用修正的声爆线化理论和优化设计方法，结合修正的声爆线化理论预测近场波形的简便特性，以及优化设计过程中设计变量数少的优点，能够直接以远场声爆强度为目标，在有效抑制严重的声爆问题的条件下，快速的设计出添加静音锥的超声速民机外形。
[0174]
下面以一个具体应用实例进一步说明本发明提出的一种针对超声速民机的旋成体静音锥快速设计方法：
[0175]
本实施例中，选取一个典型超声速民机基准构型，在该超声速民机基准构型上实现单级旋成体静音锥的快速设计，并保证添加静音锥的超声速民机构型能够有效抑制远场声爆的强度。
[0176]
步骤1，选取的超声速民机基准构型几何外形如图8所示，该超声速民机基准构型的机身长度l为71米。
[0177]
针对该构型，进行后续对旋成体静音锥快速设计的过程。
[0178]
步骤2，使用计算流体力学方法，使用ns方程对超声速民机基准构型的空间流场进行数值模拟。计算状态：马赫数ma为2.0，攻角α为5.12度。
[0179]
提取超声速民机基准构型正下方3倍机身长度处(h/l＝3.0)的基准声爆近场波形。其中，h为超声速民机正下方距离飞机飞行轨迹的距离。
[0180]
步骤3，为了说明添加设计静音锥对超声速飞机民机声爆强度的抑制效果，对基准构型的远场声爆强度进行评估。首先，基于广义burgers方程对声爆远场波形进行求解。其中，超声速飞机民机的飞行高度为18.592km，传播大气环境为标准大气。然后，以感觉声压级(pldb)作为远场声爆强度的评价指标，计算得到基准构型远场声爆强度为98.90分贝。
[0181]
步骤3，确定激波最优峰值p
s_lb
和膨胀波最优长度l
e_lb
：
[0182]
步骤3.1，超声速民机基准构型的机身头部s1开始，沿飞机轴向查找到满足以下条件的机身截面s2：机身截面s2的直径d0＝0.5％l～2.5％l；
[0183]
步骤3.2，确定机身截面s2产生的激波在基准声爆近场波形中的对应位置，表示为点b2；
[0184]
步骤3.3，在基准声爆近场波形中，截去起点o
start
到点b1以及点b1到点b2之间的基准声爆近场波形前段b1b2；保留点b2到终点o
end
之间的基准声爆近场波形后段b2o
end
；
[0185]
步骤3.4，在起点o
start
和点b2之间，依次选择三个点，分别为：点a1、点a2和点a3；
[0186]
确定点a1、点a2和点a3的初始位置后，采用优化设计算法，不断调节点a1、点a2和点a3的位置，使点a1到点a2的连接直线形成的激波a1a2的峰值ps，为激波最优峰值p
s_lb
；使点a2到点a3的连接直线形成的膨胀波a2a3的长度le，为膨胀波最优长度l
e_lb
；
[0187]
优化设计过程为：
[0188]
结合声爆远场传播和优化设计方法，将添加的激波a1a2的峰值ps和膨胀波a2a3的长度le作为设计变量，以远场声爆的感觉声压级(pldb)为目标，采用基于代理模型的优化算法对上述两个参数进行优化设计。
[0189]
其中，在将声爆波形传播到远场的过程中，将添加了激波和膨胀波的声爆近场波形作为输入，基于广义burgers方程对声爆远场波形进行求解。超声速民机的飞行高度为18.592km，传播大气环境为标准大气。优化设计中初始样本点由拉丁超立方方法产生，采用的加点准则为ei msp lcb pi mse，优化迭代步数为80步(每次迭代添加5个样本点)，计算收敛示意图如图5所示。
[0190]
优化结果表明，当远场声爆强度最低，即相应的感觉声压级为95.41分贝时，最优点处激波峰值p
s_lb
和膨胀波长度l
e_lb
分别为0.0099和12.45，为激波最优峰值p
s_lb
和膨胀波最优长度l
e_lb
。
[0191]
步骤4，对旋成体静音锥外形进行优化设计，得到最优旋成体静音锥外形：
[0192]
步骤4.1，设定旋成体静音锥外形的设计变量为：圆锥段长度l1、圆柱段长度l2和圆柱段半径r；
[0193]
步骤4.2，在满足约束条件下，给定圆锥段长度l1、圆柱段长度l2和圆柱段半径r的初始值，分别为：圆锥段长度l1＝2.0m，圆柱段长度l2＝4.0m，圆柱段半径r＝0.55m；
[0194]
步骤4.3，根据修正的声爆线化理论，获得初始旋成体静音锥外形的体积等效截面积分布的解析表达式，如下式所示：
[0195][0196]
根据修正的声爆线化理论对初始旋成体静音锥外形的线化的声爆近场波形进行预测。根据预测得到的声爆近场波形，提取静音锥对应部分的激波峰值ps＝0.00725和膨胀波长度le＝13.45。
[0197]
步骤4.4，判断得到的激波峰值ps和膨胀波长度le，是否与步骤3得到的激波最优峰值p
s_lb
和膨胀波最优长度l
e_lb
相同。该算例的优化数学模型为：
[0198]
min f＝0.8|p
s_lb-ps| 0.2|l
e_lb-le|
[0199]
w.r.t l1,l2,r
[0200]
s.t.l1∈[1.0,5.0]
[0201]
l2∈[1.0,3.0]
[0202]
r∈[0.4,0.8]
[0203]
优化设计中初始样本点由拉丁超立方方法产生，采用的加点准则为ei msp lcb pi mse。优化结果表明当参数(l1,l2,r)取值(3.5m,4.8m,0.7m)时，目标函数达到最小。
[0204]
步骤5，将最优旋成体静音锥外形的几何参数表示为：圆锥段最优长度l1(best)、圆柱段最优长度l2(best)和圆柱段最优半径r(best)；
[0205]
根据得到的静音锥参数，使用如下方法在超声速民机基准构型头部添加静音锥外
形：
[0206]
自超声速民机基准构型的机身头部s1开始，沿着飞机轴向查找到半径为0.7m的机身截面s3，截去机身头部s1到机身截面s3的一段机身；然后，在机身截面s3前面，依次生成长度为4.8m、半径0.7m的圆柱段，以及长度为3.5m的圆锥段，获得超声速民机静音锥构型，如图9所示。
[0207]
使用计算流体力学方法，使用ns方程对超声速民机基准构型的空间流场进行数值模拟。计算状态：马赫数ma为2.0，攻角α为5.12度。提取超声速民机基准构型正下方3倍机身长度处(h/l＝3.0)的基准声爆近场波形，并基于广义burgers方程对声爆远场波形进行求解，并基于本发明方法，得到旋成体静音锥构型。
[0208]
对超声速民机基准构型和本发明设计得到的旋成体静音锥构型进行声爆比较。如图10所示，为静音锥构型与超声速民机基准构型的声爆近场波形对比示意图；如图11所示，为静音锥构型与超声速民机基准构型的声爆远场波形对比示意图；其中，7为超声速民机基准构型声爆近场波形；8为超声速民机旋成体静音锥构型声爆近场波形；9为超声速民机基准构型声爆远场波形；10为超声速民机旋成体静音锥构型声爆远场波形。
[0209]
从图11可以看出，添加的静音锥使得在基准构型机头激波前又增加了一道弱激波，增加远场波形的上升时间，因此降低地面声爆强度，由基准构型的98.90分贝降低至95.65分贝。
[0210]
因此，本发明提供的针对超声速民机的旋成体静音锥快速设计方法，具有以下特点：
[0211]
本发明结合修正的声爆线化理论预测近场波形的简便特性，以及优化设计过程中设计变量数少的优点，能够直接以远场声爆强度为目标，在有效抑制声爆强度的条件下，快速的设计出添加旋成体静音锥的超声速民机外形。
[0212]
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视本发明的保护范围。

技术特征：
1.一种针对超声速民机的旋成体静音锥快速设计方法，其特征在于，包括以下步骤：步骤1，选择超声速民机基准构型；所述超声速民机基准构型的机身长度为l；步骤2，在设计马赫数ma和设计攻角α的巡航条件下，对超声速民机基准构型进行数值模拟，提取得到超声速民机基准构型正下方距离为h位置的超声速巡航状态下的基准声爆近场波形；其中，在基准声爆近场波形中，基准声爆近场波形起点为o
start
，基准声爆近场波形终点为o
end
，机身头部(s1)产生的激波在基准声爆近场波形中对应位置为点b1；步骤3，确定激波最优峰值p
s_lb
和膨胀波最优长度l
e_lb
：步骤3.1，自超声速民机基准构型的机身头部(s1)开始，沿飞机轴向查找到满足以下条件的机身截面(s2)：机身截面(s2)的直径d0＝0.5％l～2.5％l；步骤3.2，确定机身截面(s2)产生的激波在步骤2得到的基准声爆近场波形中的对应位置，表示为点b2；步骤3.3，在步骤2得到的基准声爆近场波形中，截去起点o
start
到点b1以及点b1到点b2之间的基准声爆近场波形前段b1b2；保留点b2到终点o
end
之间的基准声爆近场波形后段b2o
end
；步骤3.4，在起点o
start
和点b2之间，依次选择三个点，分别为：点a1、点a2和点a3；确定点a1、点a2和点a3的初始位置后，采用下面优化设计过程，不断调节点a1、点a2和点a3的位置，使点a1到点a2的连接直线形成的激波a1a2的峰值p
s
，为激波最优峰值p
s_lb
；使点a2到点a3的连接直线形成的膨胀波a2a3的长度l
e
，为膨胀波最优长度l
e_lb
；其中，激波a
1 a2的峰值p
s
，为点a1和点a2之间的高度差；膨胀波a
2 a3的长度l
e
，为点a2到点a3的水平距离；优化设计过程为：步骤3.4.1，确定约束条件为：条件1：激波a
1 a2的峰值p
s
＝0.1*基准声爆近场波形最大超压值～1.5*基准声爆近场波形最大超压值；条件2：激波a
1 a2的斜率＝0.008～0.012；条件3：膨胀波a
2 a3的长度l
e
＝0.1l～0.5l；条件4：点a3到点a1连线的斜率＝-0.00132～-0.00124；条件5：激波a
1 a2的斜率＝点b2到点a3连线的斜率；步骤3.4.2，在满足步骤3.4.1确定的约束条件下，给定点a1、点a2和点a3的初始位置，并在基准声爆近场波形后段b2o
end
的前方，依次连接起点o
start
、点a1、点a2、点a3和点b2，从而得到添加激波和膨胀波的声爆近场波形；结合声爆远场传播方法，以添加激波和膨胀波的声爆近场波形作为输入，使其传播到远场，并求解得到远场声爆强度；判断得到的远场声爆强度是否达到优化目标最低值；如果满足，则停止优化，点a1、点a2和点a3的当前位置，即为最优位置，进而得到激波最优峰值p
s_lb
和膨胀波最优长度l
e_lb
；如果不满足，则改变a1、点a2和点a3的位置，实质为改变激波a
1 a2的峰值p
s
和膨胀波a
2 a3的长度l
e
，重复进行迭代优化；步骤4，对旋成体静音锥外形进行优化设计，得到最优旋成体静音锥外形：步骤4.1，设定旋成体静音锥外形的设计变量为：圆锥段长度l1、圆柱段长度l2和圆柱段
半径r；步骤4.2，在满足以下约束条件下，给定圆锥段长度l1、圆柱段长度l2和圆柱段半径r的初始值，从而确定初始旋成体静音锥外形；s.t.l1∈[0.01l,0.08l]l2∈[0.01l,0.08l]r∈[0.05l,0.25l]步骤4.3，根据修正的声爆线化理论，获得初始旋成体静音锥外形的体积等效截面积分布，并对初始旋成体静音锥外形的声爆近场波形进行预测，获得初始旋成体静音锥外形的声爆近场波形；从初始旋成体静音锥外形的声爆近场波形中，提取到初始旋成体静音锥外形对应的激波峰值p
s
和膨胀波长度l
e
；步骤4.4，判断步骤4.3得到的激波峰值p
s
和膨胀波长度l
e
是否满足以下关系式：min f＝w1·
|p
s_lb-p
s
| w2·
|l
e_lb-l
e
|其中：w1和w2为权重系数，且满足w1 w2＝1.0；如果满足，则步骤4.2给定的圆锥段长度l1、圆柱段长度l2和圆柱段半径r，即为最终优化得到的旋成体静音锥外形最优参数，进而得到最优旋成体静音锥外形，然后执行步骤5；如果不满足，在满足约束条件下，调整圆锥段长度l1、圆柱段长度l2和圆柱段半径r，返回步骤4.3；步骤5，将最优旋成体静音锥外形的几何参数表示为：圆锥段最优长度l1(best)、圆柱段最优长度l2(best)和圆柱段最优半径r(best)；自超声速民机基准构型的机身头部(s1)开始，沿飞机轴向查找到半径为圆柱段最优半径r(best)的机身截面(s3)，截去机身头部(s1)到机身截面(s3)的一段机身；然后，将最优旋成体静音锥外形添加到超声速民机基准构型的机身截面(s3)的前面，得到超声速民机旋成体静音锥构型。2.根据权利要求1所述的针对超声速民机的旋成体静音锥快速设计方法，其特征在于，步骤4.3具体为：步骤4.3.1，给定初始旋成体静音锥外形的几何参数为：圆锥段长度l1、圆柱段长度l2和圆柱段半径r；采用下式，获得初始旋成体静音锥外形的体积等效截面积分布sv(x)：其中：x
11
＝l
1-r/tan(μ-α)x
12
＝l1 r/tan(μ-α)x
13
＝l1 l
2-r/tan(μ-α)
f
11
(x)＝(x-l1)
·
tan(μ-α)f
12
(x)＝2arccos(f
11
(x))f
13
(x)＝π-arccos(f
11
(x))其中：α：设计攻角；μ：步骤2的设计马赫数ma对应的马赫角；β：圆锥段锥角；并且满足：β＜μ-α；x
11
：旋成体静音锥第一轴向位置；x
12
：旋成体静音锥第二轴向位置；x
13
：旋成体静音锥第三轴向位置；x
11
、x
12
和x
13
的确定方法为：静音锥圆锥段与静音锥圆柱段的交界面为圆锥-圆柱交界面，圆锥-圆柱交界面的下轮廓与飞机对称面相交于圆锥-圆柱交界面下端点p1，圆锥-圆柱交界面的上轮廓与飞机对称面相交于圆锥-圆柱交界面上端点p2，静音锥圆柱段下轮廓与飞机对称面相交于圆柱-机身交界面下端点p3；过圆锥-圆柱交界面下端点p1，作与马赫平面平行的平面，与旋成体静音锥轴线相交于旋成体静音锥第一轴向位置x
11
；过圆锥-圆柱交界面上端点p2，作与马赫平面平行的平面，与旋成体静音锥轴线相交于旋成体静音锥第二轴向位置x
12
；过圆柱-机身交界面下端点p3，作与马赫平面平行的平面，与旋成体静音锥轴线相交于旋成体静音锥第三轴向位置x
13
；f
11
(x)：第一中间量参数；f
12
(x)：第二中间量参数；f
13
(x)：第三中间量参数；x：以初始旋成体静音锥外形的锥尖为原点，过旋成体静音锥轴线作x轴，旋成体静音锥轴线的不同位置点的x轴坐标；步骤4.3.2，根据初始旋成体静音锥外形的体积等效截面积分布sv(x)，计算声爆f函数f(τ)：其中：τ：在初始旋成体静音锥外形为巡航条件下时，在其正下方距离为h位置建立与初始旋成体静音锥外形轴线平行的τ轴，初始旋成体静音锥外形的锥尖在τ轴的投影为τ轴原点；τ轴的不同位置点具有τ轴坐标；步骤4.3.3，采用下式，计算初始旋成体静音锥外形的线化的声爆近场波形：其中：p
inf
：自由来流压强；
dp：超声速民机扰动压强与自由来流压强差；γ：大气比热比；b：普朗特-格劳厄特系数，计算公式为：步骤4.3.4，使用面积平衡法或burgers-hayers方法对步骤4.3.3得到的初始旋成体静音锥外形的线化的声爆近场波形进行非线性修正，得到最终的初始旋成体静音锥外形的声爆近场波形。

技术总结
本发明提供一种针对超声速民机的旋成体静音锥快速设计方法，包括以下步骤：选择超声速民机基准构型，提取基准声爆近场波形；截去基准声爆近场波形B1B2，添加激波A1A2与膨胀波A2A3，确定激波最优峰值和膨胀波最优长度；结合修正的声爆线化理论，对旋成体静音锥外形进行优化设计；将最优旋成体静音锥外形添加到超声速民机基准构型机身截面S3之前，得到超声速民机旋成体静音锥构型。本发明利用修正的声爆线化理论快速预测近场波形，以及优化设计变量数少的优点，能够直接以声爆强度为目标，快速设计出添加旋成体静音锥的超声速民机外形，有效抑制超声速民机巡航阶段产生严重声爆问题，更好地满足超声速民机飞行时对声爆强度的要求。求。求。


技术研发人员：张力文 韩忠华 宋文萍 乔建领 丁玉临 宋科 张科施
受保护的技术使用者：西北工业大学
技术研发日：2022.03.30
技术公布日：2022/5/25