CN116412991A - 在风洞流场试验中同步模拟飞行高度、飞行速度的方法 - Google Patents

在风洞流场试验中同步模拟飞行高度、飞行速度的方法 Download PDF

Info

Publication number
CN116412991A
CN116412991A CN202310688518.8A CN202310688518A CN116412991A CN 116412991 A CN116412991 A CN 116412991A CN 202310688518 A CN202310688518 A CN 202310688518A CN 116412991 A CN116412991 A CN 116412991A
Authority
CN
China
Prior art keywords
wind tunnel
target
regulating valve
test
pressure regulating
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.)
Granted
Application number
CN202310688518.8A
Other languages
English (en)
Other versions
CN116412991B (zh
Inventor
田嘉懿
凌忠伟
刘为杰
熊能
陈海峰
吴琦
邹岱玘
肖晋
余强
王亮
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
High Speed Aerodynamics Research Institute of China Aerodynamics Research and Development Center
Original Assignee
High Speed Aerodynamics Research Institute of China Aerodynamics Research and Development Center
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by High Speed Aerodynamics Research Institute of China Aerodynamics Research and Development Center filed Critical High Speed Aerodynamics Research Institute of China Aerodynamics Research and Development Center
Priority to CN202310688518.8A priority Critical patent/CN116412991B/zh
Publication of CN116412991A publication Critical patent/CN116412991A/zh
Application granted granted Critical
Publication of CN116412991B publication Critical patent/CN116412991B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
    • G01M9/02—Wind tunnels
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
    • G01M9/06—Measuring arrangements specially adapted for aerodynamic testing
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Physics & Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • General Physics & Mathematics (AREA)
  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)

Abstract

本发明公开了一种在风洞流场试验中同步模拟飞行高度、飞行速度的方法，涉及风洞流场控制领域，包括：S100，将试验模型目标飞行高度转换为风洞试验舱段当地静压目标；S200，根据风洞试验要求，确定风洞试验运行条件参数，进行风洞引射充压启动；S300，进行流场稳定判定，根据判定结果转入S400或进行风洞闭环控制；S400，进行风洞流场测量并判定是否完成所有风洞试验模型姿态调整，根据判定结果重新转入S300或风洞关车。本发明提供一种在风洞流场试验中同步模拟飞行高度、飞行速度的方法，能够同时实现对试验对象飞行所在高度、飞行速度的模拟，提高风洞试验模拟能力和相应试验结果的置信度。

Description

在风洞流场试验中同步模拟飞行高度、飞行速度的方法
技术领域
本发明涉及风洞流场控制领域。更具体地说，本发明涉及一种在风洞流场试验中同步模拟飞行高度、飞行速度的方法。
背景技术
精确的气动试验数据是各型航空、航天飞行器设计研制、技战术性能考核评估、气动力技术研发的基础与依据，也是提高飞行器设计研制水平、缩短研制周期、降低研制成本的关键保障。尽管计算流体力学专业近年来已取得了长足进步，能实现对大型复杂气动问题的准确求解，但是通过风洞试验对全尺寸或缩比试验模型进行吹风试验仍是获取高置信度气动数据不可或缺的关键手段。
传统风洞试验主要关注于对试验对象（模型）飞行速度进行精确模拟，通过调节风洞喷管不同的成型型面等手段以实现这一目的。但对于某些飞行器或部件，同一飞行速度在不同飞行高度下往往具有不同特性，如飞行器发动机、环控系统的高度特性。因此，在风洞试验中实现对飞行器的飞行高度、飞行速度的同步精确模拟对飞行器设计研制和性能考核评估具有十分重要的现实意义，而现有技术中难以在风洞试验中实现对飞行器的飞行高度、飞行速度的同步精确模拟。
发明内容
本发明的一个目的是解决至少上述问题和/或缺陷，并提供至少后面将说明的优点。
为了实现本发明的这些目的和其它优点，提供了一种在风洞流场试验中同步模拟飞行高度、飞行速度的方法，包括：
S100、将试验模型目标飞行高度h转换为风洞试验舱段的当地静压目标
Figure SMS_1
；
S200、根据风洞试验要求，确定相应风洞试验运行条件参数，进行风洞引射充压启动；
S300、根据风洞试验运行试验舱段当地静压
Figure SMS_2
、试验运行马赫数/>
Figure SMS_3
、引射器集气室压力进行流场稳定判定，根据判定结果转入S400，或进行风洞闭环控制；
S400、进行风洞流场测量并判定是否完成所有风洞试验模型姿态调整，根据判定结果重新转入S300，或风洞关车；
其中，在S300中，所述风洞闭环控制包括飞行速度、调压阀开度、飞行高度、风洞引射压力、引射调压阀开度的控制，且所述飞行速度和飞行高度闭环控制被配置为并行控制。
优选的是，在S100中，基于下式将试验模型目标飞行高度转换为风洞试验舱段当地静压目标
Figure SMS_4
：
Figure SMS_5
其中，H表示试验模型目标飞行高度h相对应的地势高度，R 0表示地球半径为6356.766公里，P sl 表示地面标准大气压为101.325kPa。
优选的是，在S200中，所述风洞试验运行条件参数包括：
基于试验模拟目标飞行速度确定的风洞试验运行马赫数目标
Figure SMS_6
，风洞试验舱段当地静压目标/>
Figure SMS_9
，总压初始目标/>
Figure SMS_12
，引射压力初始目标/>
Figure SMS_8
，调压阀预置开度/>
Figure SMS_10
、引射调压阀预置开度/>
Figure SMS_13
，静压偏差阈值/>
Figure SMS_14
，马赫数偏差阈值/>
Figure SMS_7
，引射压力偏差阈值/>
Figure SMS_11
；
所述风洞引射充压启动方式包括：
以开环控制方式将引射调压阀打开至预置开度
Figure SMS_15
，待风洞引射器集气室压力/>
Figure SMS_16
达到引射压力初始目标/>
Figure SMS_17
的设定阈值/>
Figure SMS_18
后，以开环控制方式将调压阀打开至预置开度，待风洞稳定段总压达到总压初始目标的设定阈值，风洞引射充压启动完成。
优选的是，在S300中，所述流场稳定判定的方式为：
将风洞试验舱段当地静压
Figure SMS_19
与当地静压目标/>
Figure SMS_20
的偏差是否小于静压偏差阈值/>
Figure SMS_21
作为判定条件一；
将试验运行马赫数
Figure SMS_22
与马赫数目标/>
Figure SMS_23
的偏差是否小于马赫数偏差阈值/>
Figure SMS_24
作为判定条件二；
将引射器集气室压力
Figure SMS_25
与引射压力目标/>
Figure SMS_26
的偏差是否小于引射压力偏差阈值
Figure SMS_27
作为判定条件三，
其中，在第一次运行时，引射压力目标
Figure SMS_28
为引射压力初始目标/>
Figure SMS_29
；
当且仅当上述三个判定条件均为真时转入S400，否则进行风洞闭环控制。
优选的是，在S300中，飞行速度控制是采用增量式马赫数PID控制器，根据试验运行马赫数
Figure SMS_30
与马赫数目标/>
Figure SMS_31
，更新调压阀开度目标/>
Figure SMS_32
，所述/>
Figure SMS_33
的获取方式为：
设在当前时刻k，试验运行马赫数
Figure SMS_34
与马赫数目标/>
Figure SMS_35
的马赫数偏差/>
Figure SMS_36
为：
Figure SMS_37
则根据
Figure SMS_38
，马赫数PID控制器通过下式计算得到/>
Figure SMS_39
：
Figure SMS_40
其中，
Figure SMS_41
、/>
Figure SMS_42
为增量式马赫数PID控制器控制参数，/>
Figure SMS_43
为上一时刻调压阀开度目标。
优选的是，在S300中，所述调压阀开度随动控制是采用增量式调压阀开度PID控制器，根据计算得到的调压阀开度目标
Figure SMS_44
和测量得到的当前时刻k调压阀开度/>
Figure SMS_45
，更新调压阀开度指令/>
Figure SMS_46
，以控制调压阀开度，所述/>
Figure SMS_47
的获取方式为：
设在当前时刻k，调压阀开度目标
Figure SMS_48
和调压阀开度/>
Figure SMS_49
的开度偏差/>
Figure SMS_50
为：
Figure SMS_51
则根据
Figure SMS_52
，增量式调压阀开度PID控制器通过下式得到/>
Figure SMS_53
：
Figure SMS_54
其中，
Figure SMS_55
、/>
Figure SMS_56
为增量式调压阀开度PID控制器控制参数，/>
Figure SMS_57
为上一时刻调压阀开度指令。
优选的是，在S300中，所述飞行高度控制由静压闭环控制器根据当前时刻k风洞试验舱段当地静压
Figure SMS_58
与当地静压目标/>
Figure SMS_59
，以及上一时刻引射压力目标/>
Figure SMS_60
，更新引射压力目标/>
Figure SMS_61
，包含以下步骤：
S301、在当前时刻k，测量风洞试验舱段当地静压
Figure SMS_62
与当地静压目标/>
Figure SMS_63
，并通过下式得到两者静压偏差/>
Figure SMS_64
：
Figure SMS_65
S302、基于S301中得到的
Figure SMS_66
，通过下式确定引射压力目标补偿量/>
Figure SMS_67
：
Figure SMS_68
其中，函数f( )表示补偿函数，
Figure SMS_69
表示风洞试验舱段当前时刻k风洞试验舱段当地静压/>
Figure SMS_70
与当地静压目标/>
Figure SMS_71
差值的绝对值，sgn( )符号函数为：
Figure SMS_72
S303、基于S302得到的
Figure SMS_73
和上一时刻引射压力目标/>
Figure SMS_74
，通过下式确定引射压力目标/>
Figure SMS_75
为：
Figure SMS_76
。
优选的是，在S300中，所述风洞引射压力随动控制中采用增量式引射压力PID控制器，根据当前时刻k引射压力目标
Figure SMS_77
和当前风洞引射器集气室压力/>
Figure SMS_78
，更新引射调压阀开度目标/>
Figure SMS_79
，所述/>
Figure SMS_80
的获取方式为：
设在当前时刻k，风洞引射器集气室压力
Figure SMS_81
与引射压力目标/>
Figure SMS_82
的引射压力偏差
Figure SMS_83
为：
Figure SMS_84
则根据
Figure SMS_85
，增量式引射压力PID控制器通过下式得到/>
Figure SMS_86
：
Figure SMS_87
其中，
Figure SMS_88
、/>
Figure SMS_89
为增量式引射压力PID控制器控制参数，/>
Figure SMS_90
为上一时刻引射调压阀开度目标。
优选的是，在S300中，所述引射调压阀开度随动控制是采用增量式引射调压阀开度PID控制器，并根据当前时刻k引射调压阀开度目标
Figure SMS_91
和测量得到的引射调压阀开度/>
Figure SMS_92
，确定引射调压阀开度指令/>
Figure SMS_93
，控制引射调压阀开度，所述/>
Figure SMS_94
的获取方式为：
设在当前时刻k，引射调压阀开度目标
Figure SMS_95
与引射调压阀开度/>
Figure SMS_96
的开度偏差/>
Figure SMS_97
为：
Figure SMS_98
则根据
Figure SMS_99
，增量式引射调压阀开度PID控制器通过下式得到/>
Figure SMS_100
：
Figure SMS_101
其中，
Figure SMS_102
、/>
Figure SMS_103
为增量式引射调压阀开度PID控制器控制参数，/>
Figure SMS_104
为上一时刻引射调压阀开度指令。
本发明至少包括以下有益效果：本发明通过将试验模型目标飞行高度转换为风洞试验舱段当地静压目标，进一步通过对风洞试验舱段中来流总压和当地静压的分别调节，来实现风洞试验中对于试验对象（模型）目标飞行高度和目标飞行速度的同步精确模拟，从而进一步提高风洞试验模拟能力和相应试验结果的置信度，满足飞行器风洞试验需求。
本发明的其它优点、目标和特征将部分通过下面的说明体现，部分还将通过对本发明的研究和实践而为本领域的技术人员所理解。
附图说明
图1为本发明的试验流程图；
图2为本发明的风洞马赫数试验曲线；
图3为本发明的风洞引射压力试验曲线；
图4为本发明的风洞试验舱静压试验曲线。
具体实施方式
下面结合附图对本发明做进一步的详细说明，以令本领域技术人员参照说明书文字能够据以实施。
S100、试验模型目标飞行高度h转换为风洞试验舱段的当地静压目标
Figure SMS_105
；
根据标准大气参数表，采用下述公式，将试验模型目标飞行高度h转换为风洞试验舱段当地静压目标
Figure SMS_106
：
Figure SMS_107
其中，H表示试验模型目标飞行高度h相对应的地势高度，R 0表示地球半径为6356.766公里，P sl 表示地面标准大气压为101.325kPa。
S200、根据风洞试验要求，确定风洞试验运行条件。
风洞试验运行条件包括依据试验模拟目标飞行速度确定的风洞试验运行马赫数目标
Figure SMS_109
，风洞试验舱段当地静压目标/>
Figure SMS_113
，总压初始目标/>
Figure SMS_115
，引射压力初始目标/>
Figure SMS_110
，调压阀预置开度/>
Figure SMS_111
、引射调压阀预置开度/>
Figure SMS_114
，静压偏差阈值/>
Figure SMS_116
，马赫数偏差阈值/>
Figure SMS_108
，引射压力偏差阈值/>
Figure SMS_112
。
S300、风洞引射充压启动。
以开环控制方式将引射调压阀打开至预置开度
Figure SMS_117
，待风洞引射器集气室压力/>
Figure SMS_120
达到引射压力初始目标/>
Figure SMS_122
的设定阈值/>
Figure SMS_119
后，以开环控制方式将调压阀打开至预置开度
Figure SMS_121
，待风洞稳定段总压/>
Figure SMS_123
达到总压初始目标/>
Figure SMS_124
的设定阈值/>
Figure SMS_118
，风洞引射充压启动完成。
S400、风洞闭环控制。
为实现对试验对象（模型）飞行高度、飞行速度的精确模拟，风洞开启闭环控制。
S410、流场稳定判定。判定风洞试验舱段当地静压
Figure SMS_125
与当地静压目标/>
Figure SMS_129
的偏差是否小于静压偏差阈值/>
Figure SMS_132
，试验运行马赫数/>
Figure SMS_127
与马赫数目标/>
Figure SMS_128
的偏差是否小于马赫数偏差阈值/>
Figure SMS_131
，以及引射器集气室压力/>
Figure SMS_134
与引射压力目标/>
Figure SMS_126
（当第一次运行时，引射压力目标
Figure SMS_130
为引射压力初始目标/>
Figure SMS_133
）偏差是否小于引射压力偏差阈值/>
Figure SMS_135
。当且仅当上述三个判定条件均为真，则直接跳转至S500，否则转入如下的闭环控制。
S420、飞行速度控制。由马赫数闭环控制器根据试验运行马赫数
Figure SMS_136
与马赫数目标
Figure SMS_137
，以及其他可能控制器补偿量，包含但不限于当前气源压力、大气压等，更新调压阀开度目标/>
Figure SMS_138
。
S430、调压阀开度随动控制。由调压阀开度控制器根据计算得到的调压阀开度目标
Figure SMS_139
和测量得到的当前调压阀开度S，更新调压阀开度指令/>
Figure SMS_140
，控制调压阀开度。
S440、飞行高度控制。由静压闭环控制器根据风洞试验舱段当地静压
Figure SMS_141
与当地静压目标/>
Figure SMS_142
，以及前一时刻引射压力目标/>
Figure SMS_143
，更新引射压力目标/>
Figure SMS_144
，包含以下步骤：
S441、测量当前风洞试验舱段当地静压
Figure SMS_145
与当地静压目标/>
Figure SMS_146
，并计算两者静压偏差:
Figure SMS_147
其中
Figure SMS_148
表示风洞试验舱段当地静压/>
Figure SMS_149
与当地静压目标/>
Figure SMS_150
差值。
S442、根据计算得到的风洞试验舱段当地静压
Figure SMS_151
与当地静压目标/>
Figure SMS_152
差值，通过下式确定引射压力目标补偿量:
Figure SMS_153
其中函数f( )表示设计的补偿函数，
Figure SMS_154
表示风洞试验舱段当地静压/>
Figure SMS_155
与当地静压目标/>
Figure SMS_156
差值的绝对值，符号函数sgn( )为:
Figure SMS_157
S443：根据计算得到的引射压力目标补偿量
Figure SMS_158
和上一时刻引射压力目标/>
Figure SMS_159
，确定为进行高度控制，修正后的引射压力目标/>
Figure SMS_160
为：
Figure SMS_161
S450、风洞引射压力随动控制。由引射压力控制器根据计算得到的修正后引射压力目标
Figure SMS_162
和当前风洞引射器集气室压力/>
Figure SMS_163
，以及其他可能的控制器补偿量，包含但不限定于当前气源压力、大气压等，确定引射调压阀开度目标/>
Figure SMS_164
。
S460、引射调压阀开度随动控制。由引射调压阀开度控制器根据计算得到的引射调压阀开度目标
Figure SMS_165
和测量得到的当前引射调压阀开度/>
Figure SMS_166
，确定引射调压阀开度指令/>
Figure SMS_167
，控制引射调压阀开度。
需要注意的是，在实际的应时，S420至S430以及S440至S460分别为同时并行运行。
S500、风洞流场测量。
试验流场稳定，达到试验要求条件，风洞测量系统采集风洞试验数据。
S510、判定是否完成所有风洞试验模型姿态调整。当判定为假（即没有完成），则跳转至S600；当判定为真（即已完成），则跳转至S700。
S600：调整风洞试验模型姿态，并跳转至S400。
重复S400至S600直至完成所有不同试验模型姿态下风洞试验数据采集，由S510跳转至S700。
S700、风洞关车。
将调压阀开度目标
Figure SMS_168
置为0，待调压阀开度S减小至关车设定目标/>
Figure SMS_169
时，将引射调压阀开度目标/>
Figure SMS_170
置为0，待引射调压阀开度/>
Figure SMS_171
和调压阀开度S均关为0，风洞关车完毕。
实施例：
S100、试验模型目标飞行高度转换为风洞试验舱段当地静压目标；
根据标准大气参数表，采用如下公式，将试验模型目标飞行高度h（8km）转换为风洞试验舱段当地静压目标
Figure SMS_172
为35.6515kPa：
Figure SMS_173
其中，H表示试验模型目标飞行高度h相对应的地势高度，R 0表示地球半径为6356.766公里，P sl 表示地面标准大气压为101.325kPa。
S200、根据风洞试验要求，确定风洞试验运行条件。
风洞试验运行条件包括依据试验模拟目标飞行速度确定的风洞试验运行马赫数目标
Figure SMS_174
为0.7，风洞试验舱段当地静压目标/>
Figure SMS_177
为35.6515 kPa，总压初始目标/>
Figure SMS_180
为50.8054kPa，引射压力初始目标/>
Figure SMS_175
为872 kPa，调压阀预置开度/>
Figure SMS_179
为241mm、引射调压阀预置开度
Figure SMS_181
为434.081mm，静压偏差阈值/>
Figure SMS_182
为0.65 kPa，马赫数偏差阈值/>
Figure SMS_176
为0.01，引射压力偏差阈值/>
Figure SMS_178
为1%。
S300、风洞引射充压启动。
以开环控制方式将引射调压阀打开至预置开度
Figure SMS_185
，待风洞引射器集气室压力/>
Figure SMS_187
达到引射压力初始目标/>
Figure SMS_189
的设定阈值/>
Figure SMS_183
（为96%），而后以开环控制方式将调压阀打开至预置开度/>
Figure SMS_186
，待风洞稳定段总压/>
Figure SMS_188
达到总压初始目标/>
Figure SMS_190
的设定阈值/>
Figure SMS_184
（为95%），风洞引射充压启动完成。
S400、风洞闭环控制。
为实现对试验对象（模型）飞行高度、飞行速度的精确模拟，风洞开启闭环控制。
S410、流场稳定判定。判定风洞试验舱段当地静压
Figure SMS_193
与当地静压目标/>
Figure SMS_194
的偏差是否小于静压偏差阈值/>
Figure SMS_197
，试验运行马赫数/>
Figure SMS_192
与马赫数目标/>
Figure SMS_196
的偏差是否小于马赫数偏差阈值/>
Figure SMS_199
，以及引射器集气室压力/>
Figure SMS_201
与引射压力目标/>
Figure SMS_191
（当第一次运行时，引射压力目标/>
Figure SMS_195
为引射压力初始目标/>
Figure SMS_198
）偏差是否小于引射压力偏差阈值/>
Figure SMS_200
。当且仅当上述三个判定条件均为真，则直接跳转至S500，否则转入如下闭环控制。
S420、飞行速度控制。优选的，采用增量式马赫数PID控制器根据试验运行马赫数
Figure SMS_202
与马赫数目标/>
Figure SMS_203
，更新调压阀开度目标。
当前时刻k试验运行马赫数
Figure SMS_204
与马赫数目标/>
Figure SMS_205
的马赫数偏差/>
Figure SMS_206
为：
Figure SMS_207
根据马赫数偏差
Figure SMS_208
，由马赫数PID控制器计算得到调压阀开度目标/>
Figure SMS_209
为：
Figure SMS_210
其中，
Figure SMS_211
、/>
Figure SMS_212
为增量式马赫数PID控制器控制参数，/>
Figure SMS_213
为上一时刻调压阀开度目标。
S430、调压阀开度随动控制。优选的，采用增量式调压阀开度PID控制器根据计算得到的调压阀开度目标
Figure SMS_214
和测量得到的当前时刻k调压阀开度/>
Figure SMS_215
，更新调压阀开度指令/>
Figure SMS_216
，以控制调压阀开度，控制调压阀开度。
当在当前时刻k，调压阀开度目标
Figure SMS_217
和调压阀开度/>
Figure SMS_218
的开度偏差/>
Figure SMS_219
为：
Figure SMS_220
根据调压阀开度偏差
Figure SMS_221
，由增量式调压阀开度PID控制器根据计算调压阀开度指令/>
Figure SMS_222
为：
Figure SMS_223
其中，
Figure SMS_224
、/>
Figure SMS_225
为增量式调压阀开度PID控制器控制参数，/>
Figure SMS_226
为上一时刻调压阀开度指令。
S440、飞行高度控制。由静压闭环控制器根据当前时刻k风洞试验舱段当地静压
Figure SMS_227
与当地静压目标/>
Figure SMS_228
，以及上一时刻引射压力目标/>
Figure SMS_229
，更新引射压力目标/>
Figure SMS_230
，包含以下步骤：
S441、测量当前时刻k风洞试验舱段当地静压
Figure SMS_231
与当地静压目标/>
Figure SMS_232
，并计算两者静压偏差/>
Figure SMS_233
：
Figure SMS_234
S442、根据计算得到的当前时刻k风洞试验舱段当地静压
Figure SMS_235
与当地静压目标/>
Figure SMS_236
差值/>
Figure SMS_237
，确定引射压力目标补偿量/>
Figure SMS_238
：
Figure SMS_239
其中，函数f( )表示补偿函数，
Figure SMS_240
表示风洞试验舱段当前时刻k风洞试验舱段当地静压/>
Figure SMS_241
与当地静压目标/>
Figure SMS_242
差值的绝对值，sgn( )符号函数为/>
Figure SMS_243
S443、根据计算得到的引射压力目标补偿量
Figure SMS_244
和上一时刻引射压力目标/>
Figure SMS_245
，确定为进行高度控制，修正后的引射压力目标/>
Figure SMS_246
为：
Figure SMS_247
S450、风洞引射压力随动控制。优选的，采用增量式引射压力PID控制器根据当前时刻k引射压力目标
Figure SMS_248
和当前风洞引射器集气室压力/>
Figure SMS_249
，更新引射调压阀开度目标/>
Figure SMS_250
。
当前时刻k，风洞引射器集气室压力
Figure SMS_251
与引射压力目标/>
Figure SMS_252
的引射压力偏差/>
Figure SMS_253
为：
Figure SMS_254
根据引射压力偏差
Figure SMS_255
，由增量式引射压力PID控制器计算得到的引射调压阀开度目标/>
Figure SMS_256
为：
Figure SMS_257
其中，
Figure SMS_258
、/>
Figure SMS_259
为增量式引射压力PID控制器控制参数，/>
Figure SMS_260
为上一时刻引射调压阀开度目标。
S460、引射调压阀开度随动控制。优选的，采用增量式引射调压阀开度PID控制器根据当前时刻k引射调压阀开度目标
Figure SMS_261
和测量得到的引射调压阀开度/>
Figure SMS_262
，确定引射调压阀开度指令/>
Figure SMS_263
，控制引射调压阀开度。
当前时刻当前时刻k，引射调压阀开度目标
Figure SMS_264
与引射调压阀开度/>
Figure SMS_265
的开度偏差
Figure SMS_266
为：
根据引射调压阀开度偏差
Figure SMS_267
，由增量式引射调压阀开度PID控制器计算得到的调压阀开度指令/>
Figure SMS_268
为：
Figure SMS_269
其中，
Figure SMS_270
、/>
Figure SMS_271
为增量式引射调压阀开度PID控制器控制参数，/>
Figure SMS_272
为上一时刻引射调压阀开度指令。
需要注意的是，S420至S430，和S440至S460为同时并行运行。
S500、风洞流场测量。
试验流场稳定，达到试验要求条件，风洞测量系统采集风洞试验数据。
S510、判定是否完成所有风洞试验模型姿态调整。当判定为假（即没有完成），则跳转至S600；当判定为真（即已完成），则跳转至S700。
S600、调整风洞试验模型姿态，并跳转至S400。
重复S400至S600直至完成所有不同试验模型姿态下风洞试验数据采集，由S510跳转至S700。
S700、风洞关车。
将调压阀开度目标
Figure SMS_273
置为0，待调压阀开度S减小至关车设定目标/>
Figure SMS_274
时，将引射调压阀开度目标/>
Figure SMS_275
置为0，待引射调压阀开度/>
Figure SMS_276
和调压阀开度均关为0，风洞关车完毕。
图1给出本发明的飞行高度、飞行速度同步精确模拟方法的试验流程图，在这一方法作用下，图2给出了风洞马赫数试验曲线，可看出在风洞开车后25s左右已将马赫数稳定控制在目标马赫数±0.01稳定带内，并在试验全程保持稳定在这一带内，实现了飞行速度的精确稳定控制；
图3给出了风洞引射压力试验曲线，可看出在风洞开车后12s左右已将引射压力稳定控制在目标引射压力±0. 1%稳定带内，并在试验全程保持稳定这一带内，且通过不断调节引射压力，进一步实现了对风洞试验舱静压的稳定控制。
进一步地，从图4给出的风洞试验舱静压试验曲线可以看出，通过调节引射压力，在风洞开车后15s左右，试验舱静压已稳定在目标试验舱静压±0.65kPa带内，实现了飞行高度的精确稳定控制，其中，图2-图4中箭头端所指区域是对箭头引出端所在区域进行的放大处理。
以上方案只是一种较佳实例的说明，但并不局限于此。在实施本发明时，可以根据使用者需求进行适当的替换和/或修改。
这里说明的设备数量和处理规模是用来简化本发明的说明的。对本发明的应用、修改和变化对本领域的技术人员来说是显而易见的。
尽管本发明的实施方案已公开如上，但其并不仅仅限于说明书和实施方式中所列运用。它完全可以被适用于各种适合本发明的领域。对于熟悉本领域的人员而言，可容易地实现另外的修改。因此在不背离权利要求及等同范围所限定的一般概念下，本发明并不限于特定的细节和这里示出与描述的图例。

Claims (9)

1.一种在风洞流场试验中同步模拟飞行高度、飞行速度的方法，其特征在于，包括：
S100、将试验模型目标飞行高度h转换为风洞试验舱段当地静压目标
Figure QLYQS_1
；
S200、根据风洞试验要求，确定相应风洞试验运行条件参数，进行风洞引射充压启动；
S300、根据风洞试验运行试验舱段当地静压
Figure QLYQS_2
、试验运行马赫数/>
Figure QLYQS_3
、引射器集气室压力进行流场稳定判定，根据判定结果转入S400，或进行风洞闭环控制；
S400、进行风洞流场测量并判定是否完成所有风洞试验模型姿态调整，根据判定结果重新转入S300，或风洞关车；
其中，在S300中，所述风洞闭环控制包括飞行速度、调压阀开度、飞行高度、引射压力、引射调压阀开度的控制，且所述飞行速度和飞行高度闭环控制被配置为并行控制。
2.如权利要求1所述的在风洞流场试验中同步模拟飞行高度、飞行速度的方法，其特征在于，在S100中，基于下式将试验模型目标飞行高度转换为风洞试验舱段当地静压目标
Figure QLYQS_4
：
Figure QLYQS_5
其中，H表示试验模型目标飞行高度h相对应的地势高度，R 0表示地球半径为6356.766公里，P sl 表示地面标准大气压为101.325kPa。
3.如权利要求1所述的在风洞流场试验中同步模拟飞行高度、飞行速度的方法，其特征在于，在S200中，所述风洞试验运行条件参数包括：
基于试验模拟目标飞行速度确定的风洞试验运行马赫数目标
Figure QLYQS_8
，风洞试验舱段当地静压目标/>
Figure QLYQS_9
，总压初始目标/>
Figure QLYQS_11
，引射压力初始目标/>
Figure QLYQS_7
，调压阀预置开度/>
Figure QLYQS_12
、引射调压阀预置开度/>
Figure QLYQS_13
，静压偏差阈值/>
Figure QLYQS_14
，马赫数偏差阈值/>
Figure QLYQS_6
，引射压力偏差阈值/>
Figure QLYQS_10
；
所述风洞引射充压启动方式包括：
以开环控制方式将引射调压阀打开至预置开度
Figure QLYQS_15
，待风洞引射器集气室压力/>
Figure QLYQS_16
达到引射压力初始目标/>
Figure QLYQS_17
的设定阈值/>
Figure QLYQS_18
后，以开环控制方式将调压阀打开至预置开度，待风洞稳定段总压达到总压初始目标的设定阈值，风洞引射充压启动完成。
4.如权利要求1所述的在风洞流场试验中同步模拟飞行高度、飞行速度的方法，其特征在于，在S300中，所述流场稳定判定的方式为：
将风洞试验舱段当地静压
Figure QLYQS_19
与当地静压目标/>
Figure QLYQS_20
的偏差是否小于静压偏差阈值/>
Figure QLYQS_21
作为判定条件一；
将试验运行马赫数
Figure QLYQS_22
与马赫数目标/>
Figure QLYQS_23
的偏差是否小于马赫数偏差阈值/>
Figure QLYQS_24
作为判定条件二；
将引射器集气室压力
Figure QLYQS_25
与引射压力目标/>
Figure QLYQS_26
的偏差是否小于引射压力偏差阈值/>
Figure QLYQS_27
作为判定条件三，其中，在第一次运行时，引射压力目标/>
Figure QLYQS_28
为引射压力初始目标/>
Figure QLYQS_29
；
当且仅当上述三个判定条件均为真时转入S400，否则进行风洞闭环控制。
5.如权利要求1所述的在风洞流场试验中同步模拟飞行高度、飞行速度的方法，其特征在于，在S300中，飞行速度控制是采用增量式马赫数PID控制器，根据试验运行马赫数
Figure QLYQS_30
与马赫数目标/>
Figure QLYQS_31
，更新调压阀开度目标/>
Figure QLYQS_32
，所述/>
Figure QLYQS_33
的获取方式为：
设在当前时刻k，试验运行马赫数
Figure QLYQS_34
与马赫数目标/>
Figure QLYQS_35
的马赫数偏差/>
Figure QLYQS_36
为：
Figure QLYQS_37
则根据
Figure QLYQS_38
，马赫数PID控制器通过下式计算得到/>
Figure QLYQS_39
：
Figure QLYQS_40
其中，
Figure QLYQS_41
、/>
Figure QLYQS_42
为增量式马赫数PID控制器控制参数，/>
Figure QLYQS_43
为上一时刻调压阀开度目标。
6.如权利要求1所述的在风洞流场试验中同步模拟飞行高度、飞行速度的方法，其特征在于，在S300中，所述调压阀开度随动控制是采用增量式调压阀开度PID控制器，根据计算得到的调压阀开度目标
Figure QLYQS_44
和测量得到的当前时刻k调压阀开度/>
Figure QLYQS_45
，更新调压阀开度指令/>
Figure QLYQS_46
，以控制调压阀开度，所述/>
Figure QLYQS_47
的获取方式为：
设在当前时刻k，调压阀开度目标
Figure QLYQS_48
和调压阀开度/>
Figure QLYQS_49
的开度偏差/>
Figure QLYQS_50
为：
Figure QLYQS_51
则根据
Figure QLYQS_52
，增量式调压阀开度PID控制器通过下式得到/>
Figure QLYQS_53
：
Figure QLYQS_54
其中，
Figure QLYQS_55
、/>
Figure QLYQS_56
为增量式调压阀开度PID控制器控制参数，/>
Figure QLYQS_57
为上一时刻调压阀开度指令。
7.如权利要求1所述的在风洞流场试验中同步模拟飞行高度、飞行速度的方法，其特征在于， 在S300中，所述飞行高度控制由静压闭环控制器根据当前时刻k风洞试验舱段当地静压
Figure QLYQS_58
与当地静压目标/>
Figure QLYQS_59
，以及上一时刻引射压力目标/>
Figure QLYQS_60
，更新引射压力目标/>
Figure QLYQS_61
，包含以下步骤：
S301、在当前时刻k，测量风洞试验舱段当地静压
Figure QLYQS_62
与当地静压目标/>
Figure QLYQS_63
，并通过下式得到两者静压偏差/>
Figure QLYQS_64
：
Figure QLYQS_65
S302、基于S301中得到的
Figure QLYQS_66
，通过下式确定引射压力目标补偿量/>
Figure QLYQS_67
：
Figure QLYQS_68
其中，函数f( )表示补偿函数，
Figure QLYQS_69
表示风洞试验舱段当前时刻k风洞试验舱段当地静压/>
Figure QLYQS_70
与当地静压目标/>
Figure QLYQS_71
差值的绝对值，sgn( )符号函数为：
Figure QLYQS_72
S303、基于S302得到的
Figure QLYQS_73
和上一时刻引射压力目标/>
Figure QLYQS_74
，通过下式确定引射压力目标
Figure QLYQS_75
为：
Figure QLYQS_76
。
8.如权利要求1所述的在风洞流场试验中同步模拟飞行高度、飞行速度的方法，其特征在于，在S300中，所述风洞引射压力随动控制中采用增量式引射压力PID控制器，根据当前时刻k引射压力目标
Figure QLYQS_77
和当前风洞引射器集气室压力/>
Figure QLYQS_78
，更新引射调压阀开度目标/>
Figure QLYQS_79
，所述/>
Figure QLYQS_80
的获取方式为：
设在当前时刻k，风洞引射器集气室压力
Figure QLYQS_81
与引射压力目标/>
Figure QLYQS_82
的引射压力偏差/>
Figure QLYQS_83
为：
Figure QLYQS_84
则根据
Figure QLYQS_85
，增量式引射压力PID控制器通过下式得到/>
Figure QLYQS_86
：
Figure QLYQS_87
其中，
Figure QLYQS_88
、/>
Figure QLYQS_89
为增量式引射压力PID控制器控制参数，/>
Figure QLYQS_90
为上一时刻引射调压阀开度目标。
9.如权利要求1所述的在风洞流场试验中同步模拟飞行高度、飞行速度的方法，其特征在于，在S300中，所述引射调压阀开度随动控制是采用增量式引射调压阀开度PID控制器，并根据当前时刻k引射调压阀开度目标
Figure QLYQS_91
和测量得到的引射调压阀开度/>
Figure QLYQS_92
，确定引射调压阀开度指令/>
Figure QLYQS_93
，控制引射调压阀开度，所述/>
Figure QLYQS_94
的获取方式为：
设在当前时刻k，引射调压阀开度目标
Figure QLYQS_95
与引射调压阀开度/>
Figure QLYQS_96
的开度偏差/>
Figure QLYQS_97
为：
Figure QLYQS_98
则根据
Figure QLYQS_99
，增量式引射调压阀开度PID控制器通过下式得到/>
Figure QLYQS_100
：
Figure QLYQS_101
其中，
Figure QLYQS_102
、/>
Figure QLYQS_103
为增量式引射调压阀开度PID控制器控制参数，/>
Figure QLYQS_104
为上一时刻引射调压阀开度指令。
CN202310688518.8A 2023-06-12 2023-06-12 在风洞流场试验中同步模拟飞行高度、飞行速度的方法 Active CN116412991B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310688518.8A CN116412991B (zh) 2023-06-12 2023-06-12 在风洞流场试验中同步模拟飞行高度、飞行速度的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310688518.8A CN116412991B (zh) 2023-06-12 2023-06-12 在风洞流场试验中同步模拟飞行高度、飞行速度的方法

Publications (2)

Publication Number Publication Date
CN116412991A true CN116412991A (zh) 2023-07-11
CN116412991B CN116412991B (zh) 2023-08-22

Family

ID=87054685

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310688518.8A Active CN116412991B (zh) 2023-06-12 2023-06-12 在风洞流场试验中同步模拟飞行高度、飞行速度的方法

Country Status (1)

Country Link
CN (1) CN116412991B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117744250A (zh) * 2024-02-21 2024-03-22 中国空气动力研究与发展中心计算空气动力研究所 一种高超飞行器流场高效复用方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5627311A (en) * 1994-08-23 1997-05-06 National Aerospace Laboratory Of Science & Technology Agency Transportable three-dimensional calibration wind tunnel system, verification method of flight control system and flight simulator using same
RU2005106495A (ru) * 2005-03-09 2006-08-20 Евгений Борисович Пашуков (RU) Способ и устройство для обнаружения вихревого образования над (перед) винтом летательного аппарата (варианты)
CN111024357A (zh) * 2019-12-11 2020-04-17 中国航天空气动力技术研究院 一种大尺寸自由活塞高焓激波风洞模拟飞行环境的方法
CN113280834A (zh) * 2021-05-24 2021-08-20 北京安达维尔民用航空技术有限公司 一种飞机综合备份电子仪表系统
CN114154441A (zh) * 2022-02-10 2022-03-08 中国空气动力研究与发展中心计算空气动力研究所 一种飞行器的环境湍流场生成与模拟计算方法
CN114383801A (zh) * 2021-12-20 2022-04-22 中国航天空气动力技术研究院 一种变体飞行器风洞虚拟飞行试验系统及方法
CN114878133A (zh) * 2022-05-18 2022-08-09 西北工业大学 一种超音速自由射流中的变马赫数试验方法
CN115290287A (zh) * 2022-10-08 2022-11-04 中国空气动力研究与发展中心低速空气动力研究所 一种高空低密度风洞试验系统、方法及风速测量方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5627311A (en) * 1994-08-23 1997-05-06 National Aerospace Laboratory Of Science & Technology Agency Transportable three-dimensional calibration wind tunnel system, verification method of flight control system and flight simulator using same
RU2005106495A (ru) * 2005-03-09 2006-08-20 Евгений Борисович Пашуков (RU) Способ и устройство для обнаружения вихревого образования над (перед) винтом летательного аппарата (варианты)
CN111024357A (zh) * 2019-12-11 2020-04-17 中国航天空气动力技术研究院 一种大尺寸自由活塞高焓激波风洞模拟飞行环境的方法
CN113280834A (zh) * 2021-05-24 2021-08-20 北京安达维尔民用航空技术有限公司 一种飞机综合备份电子仪表系统
CN114383801A (zh) * 2021-12-20 2022-04-22 中国航天空气动力技术研究院 一种变体飞行器风洞虚拟飞行试验系统及方法
CN114154441A (zh) * 2022-02-10 2022-03-08 中国空气动力研究与发展中心计算空气动力研究所 一种飞行器的环境湍流场生成与模拟计算方法
CN114878133A (zh) * 2022-05-18 2022-08-09 西北工业大学 一种超音速自由射流中的变马赫数试验方法
CN115290287A (zh) * 2022-10-08 2022-11-04 中国空气动力研究与发展中心低速空气动力研究所 一种高空低密度风洞试验系统、方法及风速测量方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
伍宗效: "飞行载荷投放模拟系统气动设计与分析", 中国优秀硕士学位论文全文数据库基础科学辑, no. 05, pages 004 - 70 *
李乾 等: "一种临近空间飞行器静/动态气动特性研究", 空气动力学学报, vol. 35, no. 04, pages 504 - 509 *
李云娇;蒋东方;李晓颖;: "一种无人机高度与空速测量系统设计", 计算机测量与控制, no. 05, pages 68 - 71 *
罗金玲;周丹;康宏琳;王济康;: "典型气动问题试验方法研究的综述", 空气动力学学报, no. 05, pages 46 - 55 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117744250A (zh) * 2024-02-21 2024-03-22 中国空气动力研究与发展中心计算空气动力研究所 一种高超飞行器流场高效复用方法
CN117744250B (zh) * 2024-02-21 2024-05-17 中国空气动力研究与发展中心计算空气动力研究所 一种高超飞行器流场高效复用方法

Also Published As

Publication number Publication date
CN116412991B (zh) 2023-08-22

Similar Documents

Publication Publication Date Title
CN106527519B (zh) 一种大型直流暂冲式超声速风洞变速压流场控制方法
CN116412991B (zh) 在风洞流场试验中同步模拟飞行高度、飞行速度的方法
CN114061890B (zh) 一种面向大型开口射流风洞下吹式静压匹配控制方法
CN103630363B (zh) 涡轮发动机高空点火能力的模拟试验方法和试验装置
CN114061891B (zh) 一种面向大型开口射流风洞下吹引射式静压匹配控制方法
CN111665857B (zh) 基于复合智能学习的变体飞行器控制方法
CN107942651B (zh) 一种近空间飞行器控制系统
CN102880055B (zh) 基于预测模型的高超声速飞行器神经网络控制方法
CN114967474B (zh) 一种基于神经网络的通用风洞流场控制方法
CN112327922B (zh) 一种飞翼无人机自主起降综合控制方法
CN111006843B (zh) 一种暂冲式超声速风洞的连续变速压方法
CN108693897A (zh) 引射驱动的闭环回流的暂冲式亚跨声速风洞流场控制方法
Davidson et al. Flight control laws for NASA's Hyper-X research vehicle
CN113008507A (zh) 一种基于暂冲气源的大流量高稳定马赫数风洞快速调节系统及方法
CN113532786A (zh) 一种用于暂冲式风洞的第二喉道控制系统及方法
CN108153323A (zh) 一种高空无人飞行器高精度再入制导方法
CN114200827B (zh) 一种超音速大机动靶标的多约束双通道控制方法
CN107315419B (zh) 一种实现飞艇远距离定点悬停的三段式控制策略
CN113267314A (zh) 一种暂冲式风洞的超声速流场总压控制系统
CN111506113A (zh) 飞行器制导指令计算方法、侧滑角计算方法及制导方法
Sun et al. Accurate modeling and homing control for parafoil delivery system based on wind disturbance rejection
CN109973221A (zh) 超声速进气道与涡扇发动机综合控制方法及装置
CN110231774A (zh) 干扰观测变进气道高超声速飞行器模糊协调控制方法
CN113253606B (zh) 一种校准箱高压供气与真空吸气联合控制系统及方法
CN113200153B (zh) 一种飞行器飞行振动预测方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant