[发明专利]一种适用于旋成体外形的来流参数确定方法

权利要求书

1.一种适用于旋成体外形的来流参数确定方法，其特征在于通过下述方式实现：

(1)建立适用于旋成体外形的表面压力近似模型：

首先，确定表面压力计算公式，将表面压力表示为来流动压q_c与压力系数C_pi的乘积与来流静压p_∞的加和；

然后，采用多项式的形式表达所述的压力系数C_pi，多项式因子为飞行攻角α、侧滑角β和来流压力比R；

最后，通过每个表面测点在多组状态下的压力数值进行拟合或者回归或者最小二乘法，得到上述多项式中的系数，得到每个表面测点对应的表面压力近似模型；

(2)获取飞行试验中的旋成体表面测点压力，根据表面测点压力结合上述近似模型进行反算，得到多项式因子，进而得到来流参数；

多项式的阶数至少4阶；

当表面测点位于旋成体的顶点、机体纵向和横向平面内时，所述的多项式表达形式简化对应的侧滑角和或攻角的奇次项，只保留侧滑角和或攻角其偶次项；

当多项式为5阶时，多项式形式如下：

f(α,β,R)＝a₀+a₁α+a₂β+a₃R+a₄α²+a₅β²+a₆R²+a₇αβ+a₈αR+a₉βR+a₁₀α³+a₁₁β³+a₁₂R³+a₁₃α²β+a₁₄αβ²+a₁₅α²R+a₁₆αR²+a₁₇β²R+a₁₈βR²+a₁₉αβR+a₂₀α⁴+a₂₁β⁴+a₂₂R⁴+a₂₃α³β+a₂₄α²β²+a₂₅αβ³+a₂₆α³R+a₂₇α²R²+a₂₈αR³+a₂₉β³R+a₃₀β²R²+a₃₁βR³+a₃₂α²βR+a₃₃αβ²R+a₃₄αβR²+a₃₅α⁵+a₃₆β⁵+a₃₆R⁵+a₃₇α⁴β+a₃₈α³β²+a₃₉α²β³+a₄₀αβ⁴+a₄₁α⁴R+a₄₂α³R²+a₄₃α²R³+a₄₄αR⁴+a₄₅β⁴R+a₄₆β³R²+a₄₆β²R³+a₄₇βR⁴+a₄₈α³βR+a₄₉αβ³R+a₅₀αβR³+a₅₁α²β²R+a₅₂α²βR²+a₅₃αβ²R²

式中，a₀～a₅₃为多项式中的系数；

当表面测点位于旋成体顶点时，多项式形式如下：

f₁(α,β,R)＝a₀+a₃R+a₄α²+a₅β²+a₆R²+a₁₂R³+a₁₅α²R+a₁₇β²R+a₂₀α⁴+a₂₁β⁴+a₂₂R⁴+a₂₄α²β²+a₂₇α²R²+a₃₀β²R²+a₃₆R⁵+a₄₁α⁴R+a₄₃α²R³+a₄₅β⁴R+a₄₆β²R³+a₅₁α²β²R

当表面测点位于旋成体的纵向平面时，多项式形式如下：

f₂(α,β,R)＝a₀+a₁α+a₃R+a₄α²+a₅β²+a₆R²+a₈αR+a₁₀α³+a₁₂R³+a₁₄αβ²+a₁₅α²R+a₁₆αR²+a₁₇β²R+a₂₀α⁴+a₂₁β⁴+a₂₂R⁴+a₂₄α²β²+a₂₆α³R+a₂₇α²R²+a₂₈αR³+a₃₀β²R²+a₃₃αβ²R+a₃₅α⁵+a₃₆R⁵+a₃₈α³β²+a₄₀αβ⁴+a₄₁α⁴R+a₄₂α³R²+a₄₃α²R³+a₄₄αR⁴+a₄₅β⁴R+a₄₆β²R³+a₅₁α²β²R+a₅₃αβ²R²

当表面测点位于旋成体的横向平面时，多项式形式如下：

f₃(α,β,R)＝a₀+a₂β+a₃R+a₄α²+a₅β²+a₆R²+a₉βR+a₁₁β³+a₁₂R³+a₁₃α²β+a₁₅α²R+a₁₇β²R+a₁₈βR²+a₂₀α⁴+a₂₁β⁴+a₂₂R⁴+a₂₄α²β²+a₂₇α²R²+a₂₉β³R+a₃₀β²R²+a₃₁βR³+a₃₂α²βR+a₃₆β⁵+a₃₆R⁵+a₃₇α⁴β+a₃₉α²β³+a₄₁α⁴R+a₄₃α²R³+a₄₅β⁴R+a₄₆β³R²+a₄₆β²R³+a₄₇βR⁴+a₅₁α²β²R+a₅₂α²βR²。