一种风洞跨声速阶梯变攻角试验的马赫数复合控制方法技术

本发明专利技术公开了一种风洞跨声速阶梯变攻角试验的马赫数复合控制方法，涉及风洞流场控制领域，在风洞跨声速阶梯变攻角试验中基于复合控制方式对马赫数进行调节，且所述复合控制包括前馈控制和反馈控制；其中，所述前馈控制被配置为采用基于局部单调的分段三次

【技术实现步骤摘要】
一种风洞跨声速阶梯变攻角试验的马赫数复合控制方法


[0001]本专利技术涉及风洞流场控制领域
。
更具体地说，本专利技术涉及一种风洞跨声速阶梯变攻角试验的马赫数复合控制方法
。

技术介绍

[0002]阶梯变攻角是一种风洞常用的试验方式，单次试验中有若干个试验模型攻角序列，需要在每个攻角状态下测量稳态的气动参数
。
跨声速阶梯变攻角试验，攻角对马赫数控制有比较明显的干扰，并且这种干扰是强非线性的，攻角的变化会造成马赫数的波动
。
[0003]现有的马赫数控制方法主要基于经典的
PID
控制或者在其基础上进行改进得控制方法，但是这类控制方法难以抑制攻角对马赫数的干扰，在阶梯变攻角过程中，马赫数控制误差通常会超过误差带，马赫数再次控稳需数秒时间
。
[0004]风洞是一种能耗巨大的设备，风洞试验通常以秒来计时，为了进一步缩短马赫数控稳时间，降低风洞试验能耗，有研究人员提出了马赫数前馈
‑
反馈复合控制
。
反馈控制仍基于
PID
控制，关键环节在于前馈控制器的设计
。
前馈控制需要根据迎角的变化对马赫数调节装置进行前馈补偿，前馈补偿量与马赫数
、
模型堵塞度
、
攻角
、
模型气动外形等多种因素有关，特别是模型的气动外形难以用特定的函数来表征
。
而现有的前馈控制方法主要是通过统计数据得到前馈补偿量，对于新的试验模型
、
试验工况，则难以保证前馈补偿量的准确性，反而会加剧马赫数的波动，特别是对于新建风洞由于缺乏历史试验数据，上述基于统计数据的前馈控制方法难以实施
。

技术实现思路

[0005]本专利技术的一个目的是解决至少上述问题和
/
或缺陷，并提供至少后面将说明的优点
。
[0006]为了实现本专利技术的这些目的和其它优点，提供了一种风洞跨声速阶梯变攻角试验的马赫数复合控制方法，在风洞跨声速阶梯变攻角试验中基于复合控制方式对马赫数进行调节，且所述复合控制包括前馈控制和反馈控制；其中，所述前馈控制被配置为采用基于局部单调的分段三次
Hermite
插值多项式
MPCHIP
实现，所述反馈控制则采用增量式
PI
控制实现
。
[0007]优选的是，所述前馈控制中
MPCHIP
插值函数的获取方式为：
S1、
基于已知样本点
(x0，
y0)
，
(x1，
y1)
，
……
，
(x
n
，
y
n
)
，建立分段插值函数
f(x)
；
S2、
对样本点中的各内点
、
各端点进行求导处理，以得到分段插值函数
f(x)
中各系数的值
。
[0008]优选的是，在
S1
中，所述分段插值函数
f(x)
在
[x
j
，
x
j+1
]分段区间内，表示为：上式中，
c1、c2、c3、c4为分段插值函数
f(x)
中的各系数；
设是一阶导数，则：由端点
x
j
和
x
j+1
两个点处的值，可以得到：令
x
j
处的差商为为
δ
j
，记为
d
j
，令
x
j+1
‑
x
j
=h
j
，则：上式中的
y
j
为已知样本点
(x0，
y0)
，
(x1，
y1)
，
……
，
(x
n
，
y
n
)
中的第
j
个点
。
[0009]优选的是，在
S2
中，各内点处的导数求取规则为：若第
j
个节点附近的差商
δ
j
‑1和
δ
j
符号相反，或者其中一个为0，则该点的一阶导数为0；若第
j
个节点附近的差商
δ
j
‑1和
δ
j
符号相同，则该点处的一阶导数
d
j
等于左右两个差商的加权平均；令
hs
j
=h
j+1
+h
j
，，，，；则：上式中，若
δ
j
=
δ
j+1
，则
δ
minj
=
δ
maxj
=
δ
j
=
δ
j+1
，若
δ
j
≠
δ
j+1
，则
δ
minj
δ
maxj
=
δ
j
=
δ
j+1
，故：
。
[0010]优选的是，在
S2
中，各端点处的导数求取方式为：采用非中点的三点公式求解端点处的一阶导数值，即采用三个点
(x0，
y0)
，
(x1，
y1)
，
(x2，
y2)
的二次插值多项式进行如下表示：
则，其导数为：将
x=x0，
x2代入上式，并化简得到：则：若
d0和
δ0符号不一致时，则将
d0设为0；若
δ0和
δ1符号不一致，即
x1处是一个极值点，令
d0=3
δ0；若
d
n
和
δ
n
‑1符号不一致，则将
d
n
为0；若
δ
n
‑2和
δ
n
‑1符号不一致，即
x
n
‑1处是一个极值点，令
d
n
=3
δ
n
‑1。
[0011]优选的是，所述前馈控制的实现方式为：
S3
，创建两个等维数的一维数组
A、S
，数组
A
用于存储攻角序列值，数组
S
用于存储流场稳定后每个攻角对应的马赫数调节结构的位置，将元素（
A[1]，
S[1]）
、
（
A[2]，
S[2]）本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.
一种风洞跨声速阶梯变攻角试验的马赫数复合控制方法，其特征在于，在风洞跨声速阶梯变攻角试验中基于复合控制方式对马赫数进行调节，且所述复合控制包括前馈控制和反馈控制；其中，所述前馈控制被配置为采用基于局部单调的分段三次
Hermite
插值多项式
MPCHIP
实现，所述反馈控制则采用增量式
PI
控制实现
。2.
如权利要求1所述的风洞跨声速阶梯变攻角试验的马赫数复合控制方法，其特征在于，所述前馈控制中
MPCHIP
插值函数的获取方式为：
S1、
基于已知样本点
(x0，
y0)
，
(x1，
y1)
，
……
，
(x
n
，
y
n
)
，建立分段插值函数
f(x)
；
S2、
对样本点中的各内点
、
各端点进行求导处理，以得到分段插值函数
f(x)
中各系数的值
。3.
如权利要求2所述的风洞跨声速阶梯变攻角试验的马赫数复合控制方法，其特征在于，在
S1
中，所述分段插值函数
f(x)
在
[x
j
，
x
j+1
]
分段区间内，表示为：上式中，
c1、c2、c3、c4为分段插值函数
f(x)
中的各系数；设是一阶导数，则：由端点
x
j
和
x
j+1
两个点处的值，可以得到：令
x
j
处的差商为为
δ
j
，记为
d
j
，令
x
j+1
‑
x
j
=h
j
，则：上式中的
y
j
为已知样本点
(x0，
y0)
，
(x1，
y1)
，
……
，
(x
n
，
y
n
)
中的第
j
个点
。4.
如权利要求3所述的风洞跨声速阶梯变攻角试验的马赫数复合控制方法，其特征在于，在
S2
中，各内点处的导数求取规则为：若第
j
个节点附近的差商
δ
j
‑1和
δ
j
符号相反，或者其中一个为0，则该点的一阶导数为0；若第
j
个节点附近的差商
δ
j
‑1和
δ
j
符号相同，则该点处的一阶导数
d
j
等于左右两个差商的加权平均；令
hs
j
=h
j+1
+h
j
，，，，；
则：上式中，若
δ
j
=
δ
j+1
，则
δ
minj
=
δ
maxj
=
δ
j
=
δ
j+1
，若
δ
j
≠
δ

【专利技术属性】
技术研发人员：刘为杰，凌忠伟，邓晓曼，李聪建，陈海峰，高川，黄昊宇，羊玺，田嘉懿，吴琦，王亮，肖晋，余强，宗鹏，
申请(专利权)人：中国空气动力研究与发展中心高速空气动力研究所，
类型：发明
国别省市：