CN104932565B - 一种高精度暂冲型引射式跨声速风洞流场控制结构 - Google Patents
一种高精度暂冲型引射式跨声速风洞流场控制结构 Download PDFInfo
- Publication number
- CN104932565B CN104932565B CN201510300756.2A CN201510300756A CN104932565B CN 104932565 B CN104932565 B CN 104932565B CN 201510300756 A CN201510300756 A CN 201510300756A CN 104932565 B CN104932565 B CN 104932565B
- Authority
- CN
- China
- Prior art keywords
- pressure
- wind tunnel
- main
- flow field
- control
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000003068 static effect Effects 0.000 claims abstract description 41
- 238000004364 calculation method Methods 0.000 claims abstract description 30
- 238000000034 method Methods 0.000 claims abstract description 30
- 239000011159 matrix material Substances 0.000 claims abstract description 26
- 238000012360 testing method Methods 0.000 claims abstract description 22
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 9
- 238000002347 injection Methods 0.000 claims description 24
- 239000007924 injection Substances 0.000 claims description 24
- 230000008569 process Effects 0.000 claims description 14
- 230000001105 regulatory effect Effects 0.000 claims description 13
- 239000013598 vector Substances 0.000 claims description 12
- 238000011010 flushing procedure Methods 0.000 claims description 9
- 206010034719 Personality change Diseases 0.000 claims description 6
- 230000008878 coupling Effects 0.000 claims description 6
- 238000010168 coupling process Methods 0.000 claims description 6
- 238000005859 coupling reaction Methods 0.000 claims description 6
- 230000004044 response Effects 0.000 claims description 6
- 230000001052 transient effect Effects 0.000 claims description 4
- 238000012937 correction Methods 0.000 claims description 3
- 239000012530 fluid Substances 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 claims description 3
- 239000000243 solution Substances 0.000 claims description 3
- 230000033228 biological regulation Effects 0.000 claims description 2
- 238000010276 construction Methods 0.000 claims description 2
- 230000008859 change Effects 0.000 abstract description 4
- 239000000203 mixture Substances 0.000 abstract 1
- 230000000694 effects Effects 0.000 description 6
- 230000000903 blocking effect Effects 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000007664 blowing Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005265 energy consumption Methods 0.000 description 2
- 230000001276 controlling effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000002706 hydrostatic effect Effects 0.000 description 1
- 238000002513 implantation Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Landscapes
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明涉及一种高精度暂冲型引射式跨声速风洞流场控制结构及方法。该控制结构将风洞流场控制分解为主引压力控制和总压静压控制两部分,由设定值计算模型、主引控制器、流场控制器三个模块构成。在主引压力控制部分,采用主调压阀控制主引压力。在总压静压控制部分,通过主排气阀和栅指控制总压和静压。设定值计算模型由一维变截面管道非定常欧拉方程组构成,通过计算流体力学软件FLUENT对风洞模型求解主引压力设定值和静压设定值。使用分段增量式PID控制算法构成主引控制器,调节主引射器出口处的主引压力。以模型姿态前馈动态矩阵控制方法为核心构成流场控制器,调节稳定段的总压和试验段的静压。该控制方法实时性好,控制精度高,试验阶段马赫数的变化控制到了0.001的范围。
Description
技术领域
本发明涉及一种暂冲型引射式跨声速风洞流场的高精度控制结构及方法,流场马赫数控制精度可达0.001。
背景技术
风洞是一种产生可控均匀气流的管状试验装置,是研究空气流动规律、航空航天飞行器及其它物体气动特性的重要地面模拟平台。风洞试验作为获取飞行器气动力数据、评估分析布局设计和预测飞行性能的主要手段,其试验精准度直接关系到飞行器布局方案的优化选择和对实际飞行性能的准确预测,从而关系到飞行器的研制质量和研制水平。
作为进行空气动力学试验的主要平台,暂冲型引射式跨声速风洞需要解决流场高精度控制以及降低能耗等问题。大飞机(我国C919大型客机、ARJ飞机等)等飞行器的研制对风洞流场马赫数的控制精度提出了严苛的要求,试验过程中流场马赫数的精度(准确度和精确度)都必须达到0.001,而作为大飞机这类我国自主在研项目唯一可依赖的跨声速气动力试验平台(2.4m暂冲式跨声速风洞)的马赫数控制精度仅为0.002~0.003。另外,风洞属于高能耗的设备。以国内2.4m暂冲式跨声速风洞为例,运行时间一般为45~90秒,试验车次超过4000次/年。按照一次试验60秒进行计算,能源动力成本高达上万元。
暂冲型引射式风洞主要有以下控制难点:①扰动频繁,主要体现在两个方面:设定值变化和模型攻角变化。空气动力学试验要求在试验阶段马赫数需要阶梯或连续变化,从而带来设定值阶跃扰动;攻角阶梯或连续变化也会带来扰动。所以,如何能够快速准确地消除扰动成为主要解决的控制问题。②时变,暂冲型引射式风洞以储气罐中的气体作为气源,由于气源容量的限制,气源压力随时间逐渐下降,导致风洞***的模型发生变化。如何能够尽量地消除气源压力对风洞流场的影响也是需要解决的控制问题。③时滞,2.4米风洞的主体结构是一个全长近200米的环形管道,控制量和被控量在实际物理位置上存在一定距离,这导致控制量的变化不能直接作用于被控量,例如主排气阀与稳定段之间的距离有近60米,使得主排气阀对总压的作用存在滞后。④耦合严重,在风洞***中,各控制量和被控量之间存在严重的耦合关系。如在使用主排气阀调节风洞流场时,总压和马赫数都会受到影响;同理,在使用栅指调节风洞流场时,总压和马赫数也都会发生变化。因此,在设计风洞流场控制器时,必须同时考虑主排气阀、栅指对总压和马赫数的影响。⑤非线性严重,暂冲型引射式风洞以总压和马赫数表征风洞流场品质,传统的控制***都采用总压和马赫数作为被控参数。但是,由于马赫数不能直接测量,根据亚音速一维等熵流理论,根据马赫数由总压和静压经公式计算得出,但这种关系是非线性的,不利于流场精确控制。
发明内容
本发明的目的是针对暂冲型引射式跨声速风洞流场中存在的控制问题,建立了一种带有解耦、前馈功能的高精度风洞流场控制结构。
一种高精度暂冲型引射式跨声速风洞流场控制结构,由主引压力控制和总压、静压控制两部分组成,包括设定值计算模型、主引控制器、流场控制器;所述设定值计算模型的输出端一端输出到主引控制器,主引控制器输出到主调压阀子***用于调压,调压后由主引射器输出,且输出的数据传输回主引控制器;所述设定值计算模型的输出端另一端输出两路到流场控制器,流场控制器,输出两路,一路输出到主排气阀子***,一路输出到栅指***,主排气阀子***和栅指***控制风洞流场分别输出总压和静压,且输出的总压和静压数据传输到流场控制器。
设定值计算模型由一维变截面管道非定常欧拉方程组构成,该模型用于根据总压设定值和马赫数设定值计算主引压力设定值和静压设定值。
所述的设定值计算模型的建立步骤如下:
步骤1.建立风洞进气处一维变截面管道非定常欧拉方程组。
利用一维等熵流理论计算主引射器喷嘴处的气体压力,密度和速度,在此基础上建立主调压阀子***处的质量方程、动量方程和能量方程。
步骤2.建立风洞排气处一维变截面管道非定常欧拉方程组。
利用一维等熵流理论计算主排气阀处的气体压力,密度和速度,在此基础上建立主排气阀子***处的质量方程、动量方程和能量方程。
步骤3.建立风洞气源、集气室气动力方程,建立风洞模型。
所述设定值计算模型的计算过程为:使用计算流体力学软件FLUENT对风洞模型求解。采用非稳态隐式求解,粘性设置为无粘。压力速度耦合采用SIMPLE算法,模型输入量为入口边界总压、静压以及出口边界静压,在稳定段和试验段添加平面作为监测面以获得所需输出量。模型中压力入口边界参数通过瞬态Profile导入到FLUENT中,软件根据时间步的推进而读取对应数据作为压力入口边界的参数值。压力出口参数设置为恒定压力(当地大气压),动量方程的离散格式选择二阶迎风格式以减少数值耗散的影响,从而计算出主引压力设定值和静压设定值。
使用分段增量式PID控制算法构成主引控制器,调节主引射器出口处的主引压力。所述的分段增量式PID控制算法将控制量分成不同的阶段,根据不同的偏差采用不同的调节参数,在大偏差时加快调节速度,小偏差时减少超调。
以模型姿态前馈动态矩阵控制方法为核心构成流场控制器,调节稳定段的总压和试验段的静压。所述的模型姿态前馈动态矩阵控制方法的计算步骤如下:
步骤1.根据风洞阶跃响应建立模型向量;
步骤2.构造一步动态矩阵和多步动态矩阵A;
步骤3.计算控制矩阵D;
步骤4.计算k-1时刻的被控量一步预测值
步骤5.采集数据。采集k时刻的总压、静压y1(k)、y2(k);
步骤6.利用控制偏差补偿k时刻的被控量一步预测值
步骤7.利用模型姿态对被控量一步预测值进行前馈补偿修;
按步骤1构造模型姿态变化对总压及模型姿态变化对静压的模型向量,记为b1、b2,其中,bi=[bi(1),bi(2),…,bi(N)]T,i=1,2;
按步骤2构造一步扰动动态矩阵具体如下:然后计算扰动前馈补偿量记k时刻的模型姿态变化量为△v(k),则扰动前馈补偿量为:
在步骤6的基础上,将扰动前馈补偿量叠加到完成对的修正,公式如下:
步骤8.计算k时刻的被控量初始预测值
步骤9.提取预测时域P内的被控量初始预测值
步骤10.获取未来P个时刻的总压、静压设定值序列w(k);
步骤11.计算k时刻的由主排气阀、栅指构成的控制向量u(k),在将k时刻的控制量施加到风洞后,转到步骤4准备计算下一时刻的控制量。
综上所述,由于采用了上述技术方案,本发明的有益效果是:
提出一种新的控制结构,通过引入主引压力控制回路,减少气源压力下降对风洞流场的影响,有效解决了时变问题。
将被控参数从总压、马赫数转换为总压、静压、主引压力,降低***的非线性度,减轻了被控参数之间的耦合,有利于提高控制性能。
提出一种适合风洞流场控制的模型姿态前馈动态矩阵控制方法,实现了风洞流场的抑制模型姿态扰动和优化解耦控制。
风洞流场控制精度得到提高,如马赫数控制精度提高至0.001,总压控制精度提高至0.1%,满足了大飞机的实验要求。
风洞流场控制精度的提高,使得实验过程中数据测量等待时间减少,从而加快了实验的进程,缩短实验时间,降低了能耗。
附图说明
本发明将通过例子并参照附图的方式说明,其中:
图1是暂冲型引射式跨声速风洞控制***图;
图2是本发明的暂冲型引射式跨声速风洞控制结构图。
具体实施方式
本发明以国内2.4m暂冲型引射式跨声速风洞为例,说明所发明的控制结构和方法的实施方式。
2.4m暂冲型引射式跨声速风洞流场控制***结构如图1所示,虚线框中为广义风洞对象。为消除风洞试验过程中气源压力的降低对风洞流场的影响,引入主引压力控制回路,将风洞流场控制分解为主引压力控制和总压静压控制两部分,如图2所示。该控制结构由设定值计算模型、主引控制器、流场控制器三个模块构成。在主引压力控制部分,采用主调压阀控制主引压力。在总压静压控制部分,通过主排气阀和栅指控制总压和静压。具体实施步骤如下:
步骤一、建立设定值计算模型
在使用图2所示的控制结构调节风洞之前,首先应建立设定值计算模型,建立过程为:
步骤1.建立风洞进气处一维变截面管道非定常欧拉方程组
在主调压阀子***处进气对风洞空气环流的影响相当于给风洞空气环流注入质量、动量和能量。由于储气罐压力与风洞压力之比总高于临界值,即主调压阀阀门喉道处(最小面积处)始终处于音速状态,因此相应的源项可表示为:
注入的质量:
注入的动量:S2=S1Ue+PeAe
注入的能量:
其中,Ae是主引射器所有喷嘴总面积,Pe,ρe,Ue分别是主引射器喷嘴处的压力,密度和速度,A*为主调压阀门的喉道截面积。由面积比Ae/A*,用下述一维等熵流理论可计算出Pe,ρe,Ue。
由上式可计算出M数,再由下面三个公式可计算出Pe,ρe,Ue。
其中,P0和ρ0为集气室内的总压与密度。
步骤2.建立风洞排气处一维变截面管道非定常欧拉方程组
在主排气阀子***排气对风洞空气环流的影响相当于从风洞空气环流中抽取质量、动量和能量。相应源项表达式为:
其中,P0e和T0e是风洞管道上排气处的总压和总温,是排气***喉道面积。上式只有当排气***喉道面积处达音速时才适用。风洞中压力与排气***出口处的压力比未达到临界值时,可按亚音速一维等熵流理论公式计算,具体计算如下:
S1=ρbUbAb
其中,Pb,ρb,Ub分别是排气***出口处的压力、密度和速度。ρe是风洞管道上排气处的密度,由于排气方向与风洞气流方向垂直,所以在风洞气流方向上动量损失为零,即
S2=0
上式中Ue是风洞管道排气处的速度,Pe是风洞管道排气处的静压,并且该式也只适用于排气阀门喉道面积处为音速时的情况,对亚临界情况,仍然要用亚音速一维等熵流公式计算出ρb、Pb和Ub,计算出S1,进而利用上式计算出S3。
步骤3.建立风洞气源、集气室气动力方程,建立风洞模型
假设:
(Ⅰ)气源、集气室容积内气体参数均匀分布。
(Ⅱ)气源、集气室容积内各质点速度U=0。
对容积V内的气体应用积分形式气体动力学方程组:
应用假设(1)、(2)对气源边界得:
Us=0
其中,为流出气源的气体质量流率,值为ρsUsAs,Vs为气源容积,Ps,ρs,Us,As分别为与气源结合处得送气管道截面上的压力、密度、速度和截面积。应用假设(1)、(2)对集气室边界得:
Um=0
其中,分别为流入、流出集气室的气体质量流率,值分别为ρCUCAC和ρBUBAB,Vm为集气室容积,XC为与集气室结合处的送气管道截面上的参数,XB为与集气室结合处的喷管截面上的参数。
步骤二、计算主引压力设定值和静压设定值
建立好风洞模型后,根据设定值模型计算主引压力设定值和静压设定值,其计算过程为:使用计算流体力学软件FLUENT对风洞模型求解。采用非稳态隐式求解,粘性设置为无粘。压力速度耦合采用SIMPLE算法,模型输入量为入口边界总压、静压以及出口边界静压,在稳定段和试验段添加平面作为监测面以获得所需输出量。模型中压力入口边界参数通过瞬态Profile导入到FLUENT中,软件根据时间步的推进而读取对应数据作为压力入口边界的参数值。压力出口参数设置为恒定压力(当地大气压),动量方程的离散格式选择二阶迎风格式以减少数值耗散的影响,从而计算出主引压力设定值和静压设定值。
步骤三、求取控制率
在每个控制周期(假定控制周期为k)内,求取主引射器、主排气阀和栅指控制率的过程如下:
首先,使用分段增量式PID控制算法求取主调压阀控制率,根据不同的偏差采用不同的调节参数,在大偏差时加快调节速度,小偏差时减少超调。增量式PID控制算法如下:
其中,e为控制偏差,Kp为比例系数,Ti为积分时间常数,Td为微分时间常数。
使用模型姿态前馈动态矩阵控制方法求取主排气阀和栅指的控制率,计算步骤如下:
步骤1.根据风洞阶跃响应建立模型向量。首先对采集到的阶跃响应数据单位化,然后利用输入输出之间的单位阶跃响应分别构造主排气阀对总压、主排气阀对静压、栅指对总压和栅指对静压的模型向量,记为a11、a21、a12、a22。下面以a11为例说明模型向量的构造过程:令a11=[a11(1)…a11(N)]T,其中a11(i)为主排气阀对总压的单位阶跃响应在iT时刻的采样值,i=1,2,…N11,T为采样周期,N11为建模时域。为了便于计算,将风洞变量的建模时域统一选取为N。按照相似的方法,可以构造a21、a12、a22。
步骤2.构造一步动态矩阵和多步动态矩阵A。动态矩阵可以依据模型向量构造,具体如下:
其中,i=1,2;j=1,2;P为预测时域;M为控制时域。P表示预测未来P个时刻的被控量,M表示所预测的控制量是根据未来M个时刻的控制增量计算的。
步骤3.计算控制矩阵D。
D=(ATQA+R)-1ATQ (5)
其中,Q为误差权矩阵,R为误差权矩阵。
步骤4.计算k-1时刻的被控量一步预测值表示:当k-1时刻控制有一增量△u(k-1)时,被控量未来N个时刻的预测值。计算公式如下:
其中,表示:在k-1时刻,假定控制作用不变时对未来N个时刻的被控量初始预测值;
步骤5.采集数据。采集k时刻的总压、静压y1(k)、y2(k)。
步骤6.利用控制偏差补偿k时刻的被控量一步预测值
其中,e(k)为k时刻的预测误差,e(k)为k时刻被控量真实值与k时刻被控量预测值的差值,H为误差权矩阵。
步骤7.利用模型姿态对被控量一步预测值进行前馈补偿修。
按步骤1构造模型姿态变化对总压及模型姿态变化对静压的模型向量,记为b1、b2。其中,bi=[bi(1),bi(2),…,bi(N)]T,i=1,2。按步骤2构造一步扰动动态矩阵具体如下:然后计算扰动前馈补偿量记k时刻的模型姿态变化量为△v(k),则扰动前馈补偿量为:
在步骤6的基础上,将扰动前馈补偿量△yd(k)叠加到完成对的修正,公式如下:
步骤8.计算k时刻的被控量初始预测值计算公式如下:
其中,S0为移位矩阵,
步骤9.提取预测时域P内的被控量初始预测值计算公式如下:
其中,I0为移位矩阵,
步骤10.获取未来P个时刻的总压、静压设定值序列w(k)。
w(k)=[sP0(k),sPS(k)]T (12)
其中,
sP0(k)=[sP0(k+1),…,sP0(k+P)]T为总压设定值,
sPS(k)=[sPS(k+1),…,sPS(k+P)]T为静压设定值。
步骤11.由主排气阀、栅指构成的控制向量u(k)。
u(k)=u(k-1)+△u(k) (14)
在将k时刻的控制量施加到风洞后,转到步骤4准备计算下一时刻的控制量。
试验效果:在国内2.4m暂冲型引射式跨声速风洞上进行控制效果的测试。考察在吹风实验过程中的马赫数和总压的控制精度。风洞试验的不同工况对应不同的马赫数、总压设定值,设定值计算模型求出的静压和主引压力的设定值如表1。
表1不同工况对应的被控变量设定值
分别选择两种阻塞度不同的飞行器模型进行控制效果对比试验,模型A的阻塞度为1%,模型B的阻塞度为0.8%。实验过程中,飞行器模型从-4°到10°阶梯变攻角,攻角变化速度2°/s,阶梯4°。选择四种吹风工况进行对比,采用控制过程中被控量的平均值和标准差两个指标衡量风洞控制***的准度和精度。控制效果如表2所示。
表2实验阶段应用验证控制性能对比表
表2显示,基于所发明的控制结构及方法,风洞试验过程中马赫数控制精度提高至0.001,总压控制精度提高至0.1%。
本发明并不局限于前述的具体实施方式。本发明扩展到任何在本说明书中披露的新特征或任何新的组合,以及披露的任一新的方法或过程的步骤或任何新的组合。
Claims (7)
1.一种高精度暂冲型引射式跨声速风洞流场控制结构,其特征在于由主引压力控制和总压静压控制两部分组成,包括设定值计算模型、主引控制器、流场控制器;所述设定值计算模型的输出端一端输出到主引控制器,主引控制器输出到主调压阀子***用于调压,调压后由主引射器输出,且输出的数据传输回主引控制器;所述设定值计算模型的输出端另一端输出两路到流场控制器,流场控制器输出两路,一路输出到主排气阀子***,一路输出到栅指***,主排气阀子***和栅指***控制风洞流场分别输出总压和静压,且输出的总压和静压数据传输回流场控制器;
所述流场控制器采用模型姿态前馈动态矩阵控制方法为核心,其具体计算方法如下:
步骤1:根据风洞阶跃响应建立模型向量;
步骤2:构造一步动态矩阵和多步动态矩阵A;
步骤3:计算控制矩阵D;
步骤4:计算k-1时刻的被控量一步预测值
步骤5:采集k时刻的总压、静压y1(k)、y2(k);
步骤6:利用控制偏差补偿k时刻的被控量一步预测值
步骤7:利用模型姿态对被控量一步预测值进行前馈补偿修正得到
步骤8:计算k时刻的被控量初始预测值
步骤9:提取预测时域P内的被控量初始预测值
步骤10:获取未来P个时刻的总压、静压设定值序列w(k);
步骤11:计算k时刻的由主排气阀、栅指构成的控制向量u(k), 在将k时刻的控制量施加到风洞后,转到步骤4准备计算下一时刻的控制量。
2.根据权利要求1所述的一种高精度暂冲型引射式跨声速风洞流场控制结构,其特征在于所述主引控制器使用分段增量式PID控制算法构成,调节主引射器出口处的主引压力。
3.根据权利要求1所述的一种高精度暂冲型引射式跨声速风洞流场控制结构,其特征在于:
按步骤1构造模型姿态变化对总压及模型姿态变化对静压的模型向量,记为b1、b2,其中,bi=[bi(1),bi(2),…,bi(N)]T,i=1,2;
按步骤2构造一步扰动动态矩阵具体如下:然后计算扰动前馈补偿量记k时刻的模型姿态变化量为△v(k),则扰动前馈补偿量为:
4.根据权利要求1所述的一种高精度暂冲型引射式跨声速风洞流场控制结构,其特征在于:在步骤6的基础上,将扰动前馈补偿量叠加到完成对的修正,公式如:
5.根据权利要求1所述的一种高精度暂冲型引射式跨声速风洞流场控制结构,其特征在于所述设定值计算模型的构建步骤如下:
建立风洞进气处一维变截面管道非定常欧拉方程组;利用一维等熵流理论计算主引射器喷嘴处的气体压力,密度和速度,在此基础上建立主调压阀子***处的质量方程、动量方程和能量方程;
建立风洞排气处一维变截面管道非定常欧拉方程组;利用一维等熵流理论计算主排气阀处的气体压力,密度和速度,在此基础上 建立主排气阀子***处的质量方程、动量方程和能量方程;
建立风洞气源、集气室气动力方程,建立风洞模型。
6.根据权利要求5所述的一种高精度暂冲型引射式跨声速风洞流场控制结构,其特征在于所述设定值计算模型的计算过程为:使用计算流体力学软件FLUENT对风洞模型求解;采用非稳态隐式求解,粘性设置为无粘;压力速度耦合采用SIMPLE算法,模型输入量为入口边界总压、静压以及出口边界静压,在稳定段和试验段添加平面作为监测面以获得所需输出量;模型中压力入口边界参数通过瞬态Profile导入到FLUENT中,软件根据时间步的推进而读取对应数据作为压力入口边界的参数值;压力出口参数设置为恒定压力,动量方程的离散格式选择二阶迎风格式以减少数值耗散的影响,从而计算出主引压力设定值和静压设定值。
7.根据权利要求6所述的一种高精度暂冲型引射式跨声速风洞流场控制结构,其特征在于所述恒定压力为风洞所在处的大气压力。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510300756.2A CN104932565B (zh) | 2015-06-04 | 2015-06-04 | 一种高精度暂冲型引射式跨声速风洞流场控制结构 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510300756.2A CN104932565B (zh) | 2015-06-04 | 2015-06-04 | 一种高精度暂冲型引射式跨声速风洞流场控制结构 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104932565A CN104932565A (zh) | 2015-09-23 |
CN104932565B true CN104932565B (zh) | 2017-06-06 |
Family
ID=54119769
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510300756.2A Active CN104932565B (zh) | 2015-06-04 | 2015-06-04 | 一种高精度暂冲型引射式跨声速风洞流场控制结构 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104932565B (zh) |
Families Citing this family (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105424311B (zh) * | 2015-11-10 | 2017-12-15 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种大长细比带尾舵旋成体模型风洞测力试验方法 |
CN105388764B (zh) * | 2015-12-15 | 2017-12-08 | 重庆科技学院 | 基于动态矩阵前馈预测的电液伺服pid控制方法及*** |
CN106840582B (zh) * | 2016-12-23 | 2019-04-09 | 中国航天空气动力技术研究院 | 一种脉冲风洞运行控制***及脉冲风洞运行控制方法 |
CN106919747B (zh) * | 2017-02-24 | 2020-12-18 | 中国航天空气动力技术研究院 | 一种飞行器舱内压力预测方法 |
CN108572054B (zh) * | 2017-03-13 | 2020-01-07 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种暂冲式风洞模拟试验方法及装置 |
CN106979852B (zh) * | 2017-03-24 | 2019-02-26 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种适用于全尺寸内埋武器舱高速风洞试验的流场控制方法 |
CN107103123A (zh) * | 2017-04-06 | 2017-08-29 | 大连大学 | 基于Bagging‑FSE算法的风洞马赫数混合模型建模方法 |
CN107272412B (zh) * | 2017-07-13 | 2020-06-09 | 中国航天空气动力技术研究院 | 一种暂冲式风洞流场控制的辩识方法 |
CN108225712B (zh) * | 2017-12-29 | 2021-07-27 | 中国空气动力研究与发展中心设备设计及测试技术研究所 | 连续式跨声速风洞喷管段的半柔壁喉块与柔板控制方法 |
CN109115447A (zh) * | 2018-08-01 | 2019-01-01 | 东北大学 | 超声速风洞空气流场模型非线性阶次辨识方法 |
CN108828954B (zh) * | 2018-08-15 | 2021-11-02 | 苏州佐竹冷热控制技术有限公司 | 气候风洞自适应预测控制***及其控制方法 |
CN110161841A (zh) * | 2019-06-05 | 2019-08-23 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种适用于暂冲式跨声速风洞的前馈-模糊pid控制方法 |
CN111006843B (zh) * | 2019-12-04 | 2021-09-07 | 中国航天空气动力技术研究院 | 一种暂冲式超声速风洞的连续变速压方法 |
CN112462813B (zh) * | 2020-12-01 | 2022-10-25 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种基于虚拟阀位解算的暂冲式高速风洞调压阀安全控制装置及控制方法 |
CN113268824B (zh) * | 2021-04-23 | 2023-07-11 | 上海卫星工程研究所 | 基于温控阀内流场的流体回路数值模拟方法、***及介质 |
CN113375897B (zh) * | 2021-05-26 | 2022-09-20 | 北京航空航天大学 | 一种面向暂冲式风洞***的建模仿真方法 |
CN113532786B (zh) * | 2021-06-10 | 2024-05-03 | 中国航天空气动力技术研究院 | 一种用于暂冲式风洞的第二喉道控制***及方法 |
CN114061891B (zh) * | 2022-01-18 | 2022-03-29 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种面向大型开***流风洞下吹引射式静压匹配控制方法 |
CN114061890B (zh) * | 2022-01-18 | 2022-03-29 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种面向大型开***流风洞下吹式静压匹配控制方法 |
CN114166459B (zh) * | 2022-02-11 | 2022-04-19 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种风洞群中压调度方法 |
CN114185265B (zh) * | 2022-02-15 | 2022-04-15 | 中国空气动力研究与发展中心高速空气动力研究所 | 大型开***流风洞超声速定总压连续变马赫数控制方法 |
CN114185266B (zh) * | 2022-02-15 | 2022-04-22 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种适用于暂冲式风洞的总压复合控制方法 |
CN116499700B (zh) * | 2023-06-26 | 2023-09-01 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种风洞主引射压力分段控制方法及*** |
CN116499701B (zh) * | 2023-06-26 | 2023-09-01 | 中国航空工业集团公司沈阳空气动力研究所 | 暂冲式跨声速风洞马赫数控制方法、电子设备及存储介质 |
CN116754176B (zh) * | 2023-08-22 | 2023-10-24 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种暂冲式高速风洞气源压力在线精确估计方法 |
CN117687446B (zh) * | 2024-02-04 | 2024-04-16 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种暂冲式风洞主引射器压力控制方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11237304A (ja) * | 1998-02-23 | 1999-08-31 | Kawasaki Heavy Ind Ltd | 吹出し式風洞の制御装置 |
US7401505B1 (en) * | 2007-01-09 | 2008-07-22 | The United States Of America As Represented By The Secretary Of The Navy | Low cost wind tunnel for supersonic and hypersonic aerothermal testing |
CN102539107A (zh) * | 2012-02-17 | 2012-07-04 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种实现风洞试验信号精确同步的方法 |
CN103123504A (zh) * | 2012-12-18 | 2013-05-29 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种适用于暂冲式跨声速风洞的变速压流场控制方法 |
CN103135624A (zh) * | 2012-12-19 | 2013-06-05 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种带引射功能的暂冲式超声速风洞控制方法 |
CN104048807A (zh) * | 2014-03-21 | 2014-09-17 | 西北工业大学 | 连续式跨声速风洞实验段流场变雷诺数控制方法 |
-
2015
- 2015-06-04 CN CN201510300756.2A patent/CN104932565B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11237304A (ja) * | 1998-02-23 | 1999-08-31 | Kawasaki Heavy Ind Ltd | 吹出し式風洞の制御装置 |
US7401505B1 (en) * | 2007-01-09 | 2008-07-22 | The United States Of America As Represented By The Secretary Of The Navy | Low cost wind tunnel for supersonic and hypersonic aerothermal testing |
CN102539107A (zh) * | 2012-02-17 | 2012-07-04 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种实现风洞试验信号精确同步的方法 |
CN103123504A (zh) * | 2012-12-18 | 2013-05-29 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种适用于暂冲式跨声速风洞的变速压流场控制方法 |
CN103135624A (zh) * | 2012-12-19 | 2013-06-05 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种带引射功能的暂冲式超声速风洞控制方法 |
CN104048807A (zh) * | 2014-03-21 | 2014-09-17 | 西北工业大学 | 连续式跨声速风洞实验段流场变雷诺数控制方法 |
Non-Patent Citations (2)
Title |
---|
2.4m ×2.4m 跨声速风洞多变量控制;李尚春 等;《贵州科学》;20021231;第20卷(第4期);正文第140-142页 * |
暂冲式风洞马赫数模糊自适应控制***研究;洪志伟;《中国优秀硕士学位论文全文数据库工程科技II辑》;20150515;正文第7-9页,图2.5 * |
Also Published As
Publication number | Publication date |
---|---|
CN104932565A (zh) | 2015-09-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104932565B (zh) | 一种高精度暂冲型引射式跨声速风洞流场控制结构 | |
WO2019237320A1 (zh) | 一种基于模型预测的航空发动机在线优化及多变量控制设计方法 | |
CN108762089B (zh) | 一种基于模型预测的航空发动机在线优化及多变量控制设计方法 | |
CN106647253B (zh) | 航空发动机分布式控制***多性能鲁棒跟踪控制方法 | |
CN108304684B (zh) | 一种火箭发动机尾喷射流仿真方法及*** | |
CN108108562B (zh) | 一种基于高斯分布的风电机组尾流的解析建模方法 | |
Cunningham et al. | Modern wing flutter analysis by computational fluid dynamics methods | |
CN111006843B (zh) | 一种暂冲式超声速风洞的连续变速压方法 | |
CN106679925A (zh) | 一种微质量射流流量高精度控制装置和控制方法 | |
CN101887267A (zh) | 风洞马赫数控制器 | |
CN113267314A (zh) | 一种暂冲式风洞的超声速流场总压控制*** | |
CN110866336B (zh) | 一种单孔超音速氧枪喷嘴设计优化*** | |
CN116736723A (zh) | 一种气动热影响下的弹性高超声速飞行器建模和模糊自适应滑模控制方法 | |
Chung et al. | Real time analytical linearization of turbofan engine model | |
McKernan et al. | A linear state-space representation of plane Poiseuille flow for control design: a tutorial | |
Starshak et al. | Method of characteristics design of high-temperature wind tunnel nozzles with vibrational relaxation | |
Prasannakumar et al. | Design and power calculation of HLFC suction system for a subsonic short-range aircraft | |
Yu et al. | Perturbation analysis and control of Mach number 2.4-meter transonic wind tunnel | |
Mikhail et al. | Computation of a supersonic flow past an axisymmetric nozzle boattail with jet exhaust | |
Fang et al. | Design optimization of unsteady airfoils with continuous adjoint method | |
Konig et al. | Transonic Lattice Boltzmann Simulations of the NASA-CRM in the European Transonic Windtunnel | |
Thormann et al. | Efficient aerodynamic derivative calculation in three-dimensional transonic flow | |
Pralits | Optimal design of natural and hybrid laminar flow control on wings | |
CN107489531B (zh) | 基于半积分与增益自适应的高超声速发动机供油规律设计方法 | |
Harish et al. | Simulation of control of pressure regulating valve in high speed wind tunnels PID control in LabVIEW |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |