CN114018531A - 一种高超声速风洞连续变总压的试验数据处理方法 - Google Patents
一种高超声速风洞连续变总压的试验数据处理方法 Download PDFInfo
- Publication number
- CN114018531A CN114018531A CN202111316804.9A CN202111316804A CN114018531A CN 114018531 A CN114018531 A CN 114018531A CN 202111316804 A CN202111316804 A CN 202111316804A CN 114018531 A CN114018531 A CN 114018531A
- Authority
- CN
- China
- Prior art keywords
- total pressure
- alpha
- wind tunnel
- hypersonic wind
- fitting
- 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
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
- G01M9/065—Measuring arrangements specially adapted for aerodynamic testing dealing with flow
-
- 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/08—Aerodynamic models
-
- 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
本发明公开了一种高超声速风洞连续变总压的试验数据处理方法。本发明的高超声速风洞连续变总压的试验数据处理方法包括以下步骤:a.对高超声速风洞不同总压条件下的流场马赫数进行修正;b.对高超声速风洞连续变总压的试验数据进行降噪和去干扰处理;c.对飞行器模型进行真实姿态角修正。本发明的高超声速风洞连续变总压的试验数据处理方法能较为完整、准确地获得飞行器模型气动力随雷诺数的变化规律。
Description
技术领域
本发明属于高超声速风洞试验技术领域,具体涉及一种高超声速风洞连续变总压的试验数据处理方法。
背景技术
当前高超声速技术飞速发展,各种复杂外形的高超声速飞行器不断出现。这些高超声速飞行器设计研制离不开气动数据的支撑,高超声速风洞试验是获取其气动数据的有效手段之一。
目前,在研的高超声速飞行器的气动性能对雷诺数变化十分敏感,因此,需要通过地面试验模拟高超声速飞行器随雷诺数连续变化的规律,以便为高超声速飞行器设计和控制系统设计提供输入。高超声速风洞连续变总压试验能够获取高超声速飞行器随雷诺数连续变化的规律。但是,相比于常规高超声速风洞试验数据处理,高超声速风洞连续变总压试验数据处理需要解决以下几个问题:
1.不同总压条件下,高超声速流场的真实马赫数估值不同。通常,流场的马赫数是根据流场校测结果给定的,但由于试验成本的限制,往往仅开展了典型总压状态的流场马赫数校准。但是,对于高超声速风洞连续变总压试验,流场总压状态是连续变化的,无法直接应用流场校测数据得到总压连续变化状态下的流场马赫数,需要根据流场校测数据对流场马赫数进行修正;
2.相比于常规高超声速风洞测力试验,高超声速风洞连续变总压测力试验数据存在更多的噪声和干扰,因此需要对高超声速风洞连续变总压测力试验数据进行降噪和去干扰处理;
3.不同总压条件下,飞行器模型受到的载荷会不同,使得应变天平和天平支杆的弹性变形不同,从而使得飞行器模型的姿态角发生变化。如果不修正模型姿态角,则获取的飞行器模型气动力系数随总压变化规律不准确,因此需要对飞行器模型姿态角进行修正。
当前,亟需发展一种高超声速风洞连续变总压的试验数据处理方法。
发明内容
本发明所要解决的技术问题是提供一种高超声速风洞连续变总压的试验数据处理方法。
本发明的高超声速风洞连续变总压的试验数据处理方法,包括以下步骤:
a.对高超声速风洞不同总压条件下的流场马赫数进行修正;
a1.统计高超声速风洞试验段的模型安装位置的不同总压状态下流场校测数据,并绘制成表格;
a2.将步骤a1的表格,取横坐标为流场总压,纵坐标为流场马赫数,绘制成散点图,对散点图进行多项式拟合,得到马赫数和总压的对应关系曲线;
a3.通过步骤a2获得的马赫数和总压的对应关系曲线,查找任一总压值对应的流场实际马赫数;
b.对高超声速风洞连续变总压的试验数据进行降噪和去干扰处理;
b1.原始数据低通滤波;
b2.采用小波分析方法对试验数据进行降噪和扣除干扰处理;
c.对飞行器模型进行真实姿态角修正;
c1.根据不同总压条件下,飞行器模型所承受的载荷,得到天平和支杆的弹性角Δα,进而得到飞行器模型的实际姿态角α实际=α名义+Δα;
c2.根据固定总压条件下,飞行器模型气动力系数随攻角变化曲线,采用二次多项式进行拟合,得到气动力系数在变总压试验名义攻角α名义附近的气动力系数与攻角的拟合曲线f(C~α);
c3.通过步骤c1的拟合曲线f(C~α),得到名义攻角α名义的气动力系数C拟合,α名义和真实攻角α真实的气动力系数C拟合,α真实;
c4.得到不同总压条件下,模型在名义攻角α名义的真实气动数据:
C试验,α真实-C试验,α名义=C拟合,α真实-C拟合,α名义
即:
C试验,α真实=C试验,α名义+(C拟合,α真实-C拟合,α名义)
进一步地,所述的步骤b1的低通滤波的截止频率为3Hz。
进一步地,所述的步骤b1的低通滤波方法为高超声速风洞采集系统硬件滤波或者数字滤波。
本发明的高超声速风洞连续变总压的试验数据处理方法能较为完整、准确地获得飞行器模型气动力随雷诺数的变化规律。
附图说明
图1为本发明的高超声速风洞连续变总压的试验数据处理方法流程图;
图2为Φ1米高超声速风洞流场校测获得的P0~Ma曲线;
图3为高超声速风洞连续变总压试验获得的原始试验数据(法向力CN);
图4为高超声速风洞连续变总压试验的原始试验数据与进行小波分析降噪和去干扰处理后的试验数据对比(法向力CN);
图5为小波分析降噪和去干扰处理后的试验数据在小波分析后,姿态角修正前后的数据对比(法向力CN)。
具体实施方式
下面结合附图和实施例详细说明本发明。
如图1所示,本发明的高超声速风洞连续变总压的试验数据处理方法,包括以下步骤:
a.对高超声速风洞不同总压条件下的流场马赫数进行修正;
a1.统计高超声速风洞试验段的模型安装位置的不同总压状态下流场校测数据,并绘制成表格;
a2.将步骤a1的表格,取横坐标为流场总压,纵坐标为流场马赫数,绘制成散点图,对散点图进行多项式拟合,得到马赫数和总压的对应关系曲线;
a3.通过步骤a2获得的马赫数和总压的对应关系曲线,查找任一总压值对应的流场实际马赫数;
b.对高超声速风洞连续变总压的试验数据进行降噪和去干扰处理;
b1.原始数据低通滤波;
b2.采用小波分析方法对试验数据进行降噪和扣除干扰处理;
c.对飞行器模型进行真实姿态角修正;
c1.根据不同总压条件下,飞行器模型所承受的载荷,得到天平和支杆的弹性角Δα,进而得到飞行器模型的实际姿态角α实际=α名义+Δα;
c2.根据固定总压条件下,飞行器模型气动力系数随攻角变化曲线,采用二次多项式进行拟合,得到气动力系数在变总压试验名义攻角α名义附近的气动力系数与攻角的拟合曲线f(C~α);
c3.通过步骤c1的拟合曲线f(C~α),得到名义攻角α名义的气动力系数C拟合,α名义和真实攻角α真实的气动力系数C拟合,α真实;
c4.得到不同总压条件下,模型在名义攻角α名义的真实气动数据:
C试验,α真实-C试验,α名义=C拟合,α真实-C拟合,α名义
即:
C试验,α真实=C试验,α名义+(C拟合,α真实-C拟合,α名义)
进一步地,所述的步骤b1的低通滤波的截止频率为3Hz。
进一步地,所述的步骤b1的低通滤波方法为高超声速风洞采集系统硬件滤波或者数字滤波。
实施例1
本实施例的飞行器模型为升力体模型,高超声速风洞为CARDC的Φ1米高超声速风洞,流场校测数据见表1,获得的M~P0曲线见图2。进行高超声速风洞连续变总压测力试验。
本实施例的低通滤波的截止频率取3Hz;选用Φ1米高超声速风洞的高精度放大器PFI28000的巴特沃思低通滤波器进行硬件滤波;获得的CN~P0原始数据见图3,进行小波分析降噪和去干扰处理后的CN~P0对比曲线见图4,姿态角修正前后的CN~P0曲线见图5。
从图5可知,试验数据进行小波分析降噪和去干扰处理后,数据干扰噪声明显降低。未进行姿态角修正之前,图中显示法向力系数CN随总压增大而增大;进行姿态角修正之后,图中显示法向力系数CN随总压增大变化较小,量值基本保持不变。即:未进行姿态角修正之前,法向力系数CN随总压变化规律是不正确的,有必要进行姿态角修正。
尽管本发明的实施方案已公开如上,但其并不仅仅限于说明书和实施方式中所列运用,它完全可以被适用于各种适合本发明的领域。对于熟悉本领域的人员而言,在不脱离本发明原理的前提下,可容易地实现另外的改进和润饰,因此在不背离权利要求及等同范围所限定的一般概念下,本发明并不限于特定的细节和这里示出与描述的图例。
表1
Claims (3)
1.一种高超声速风洞连续变总压的试验数据处理方法,其特征在于,包括以下步骤:
a.对高超声速风洞不同总压条件下的流场马赫数进行修正;
a1.统计高超声速风洞试验段的模型安装位置的不同总压状态下流场校测数据,并绘制成表格;
a2.将步骤a1的表格,取横坐标为流场总压,纵坐标为流场马赫数,绘制成散点图,对散点图进行多项式拟合,得到马赫数和总压的对应关系曲线;
a3.通过步骤a2获得的马赫数和总压的对应关系曲线,查找任一总压值对应的流场实际马赫数;
b.对高超声速风洞连续变总压的试验数据进行降噪和去干扰处理;
b1.原始数据低通滤波;
b2.采用小波分析方法对试验数据进行降噪和扣除干扰处理;
c.对飞行器模型进行真实姿态角修正;
c1.根据不同总压条件下,飞行器模型所承受的载荷,得到天平和支杆的弹性角Δα,进而得到飞行器模型的实际姿态角α实际=α名义+Δα;
c2.根据固定总压条件下,飞行器模型气动力系数随攻角变化曲线,采用二次多项式进行拟合,得到气动力系数在变总压试验名义攻角α名义附近的气动力系数与攻角的拟合曲线f(C~α);
c3.通过步骤c1的拟合曲线f(C~α),得到名义攻角α名义的气动力系数C拟合,α名义和真实攻角α真实的气动力系数C拟合,α真实;
c4.得到不同总压条件下,模型在名义攻角α名义的真实气动数据:
C试验,α真实-C试验,α名义=C拟合,α真实-C拟合,α名义
即:
C试验,α真实=C试验,α名义+(C拟合,α真实-C拟合,α名义)
2.根据权利要求1所述的高超声速风洞连续变总压的试验数据处理方法,其特征在于,所述的步骤b1的低通滤波的截止频率为3Hz。
3.根据权利要求1所述的高超声速风洞连续变总压的试验数据处理方法,其特征在于,所述的步骤b1的低通滤波方法为高超声速风洞采集系统硬件滤波或者数字滤波。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111316804.9A CN114018531B (zh) | 2021-11-09 | 2021-11-09 | 一种高超声速风洞连续变总压的试验数据处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111316804.9A CN114018531B (zh) | 2021-11-09 | 2021-11-09 | 一种高超声速风洞连续变总压的试验数据处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114018531A true CN114018531A (zh) | 2022-02-08 |
CN114018531B CN114018531B (zh) | 2023-03-21 |
Family
ID=80062653
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111316804.9A Active CN114018531B (zh) | 2021-11-09 | 2021-11-09 | 一种高超声速风洞连续变总压的试验数据处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114018531B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115615655A (zh) * | 2022-11-15 | 2023-01-17 | 中国航空工业集团公司沈阳空气动力研究所 | 一种连续式跨声速风洞试验数据干扰修正方法 |
CN116929703A (zh) * | 2023-09-18 | 2023-10-24 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种考虑堵塞效应的低温风洞马赫数确定方法及应用 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2942038A1 (fr) * | 2009-02-09 | 2010-08-13 | Airbus France | Procede de determination d'une valeur rigide d'un coefficient aerodynamique |
CN104048807A (zh) * | 2014-03-21 | 2014-09-17 | 西北工业大学 | 连续式跨声速风洞实验段流场变雷诺数控制方法 |
RU2547473C1 (ru) * | 2013-12-12 | 2015-04-10 | Федеральное государственное унитарное предприятие "Центральный аэрогидродинамический институт имени профессора Н.Е. Жуковского" (ФГУП "ЦАГИ") | Рабочая часть аэродинамической трубы |
CN105716826A (zh) * | 2016-02-18 | 2016-06-29 | 江西洪都航空工业集团有限责任公司 | 一种零升阻力系数的雷诺数效应修正方法 |
CN105973563A (zh) * | 2016-05-12 | 2016-09-28 | 中国科学院力学研究所 | 一种强激光与超声速风洞联合破坏效应的试验系统及方法 |
JP2017049014A (ja) * | 2015-08-31 | 2017-03-09 | 公益財団法人鉄道総合技術研究所 | 騒音測定方法及び騒音算出装置 |
CN110160734A (zh) * | 2019-05-23 | 2019-08-23 | 华中科技大学 | 基于多孔介质的风洞喷管、主动降噪装置及方法 |
CN110763420A (zh) * | 2019-10-23 | 2020-02-07 | 中国空气动力研究与发展中心高速空气动力研究所 | 低温风洞流场自抗扰控制方法 |
CN111351628A (zh) * | 2020-03-27 | 2020-06-30 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种高超声速风洞连续变攻角试验数据处理方法 |
CN111459175A (zh) * | 2019-12-18 | 2020-07-28 | 北京航空航天大学 | 一种基于l1自适应控制器的尾座式无人机轨迹跟踪容错控制方法 |
CN112304563A (zh) * | 2020-10-30 | 2021-02-02 | 中国空气动力研究与发展中心超高速空气动力研究所 | 研究转捩对高超声速飞行器气动特性影响的风洞试验方法 |
CN112525474A (zh) * | 2020-12-22 | 2021-03-19 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种获得跨声速风洞堵塞干扰因子的方法 |
CN113267314A (zh) * | 2021-04-25 | 2021-08-17 | 中国航天空气动力技术研究院 | 一种暂冲式风洞的超声速流场总压控制系统 |
-
2021
- 2021-11-09 CN CN202111316804.9A patent/CN114018531B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2942038A1 (fr) * | 2009-02-09 | 2010-08-13 | Airbus France | Procede de determination d'une valeur rigide d'un coefficient aerodynamique |
RU2547473C1 (ru) * | 2013-12-12 | 2015-04-10 | Федеральное государственное унитарное предприятие "Центральный аэрогидродинамический институт имени профессора Н.Е. Жуковского" (ФГУП "ЦАГИ") | Рабочая часть аэродинамической трубы |
CN104048807A (zh) * | 2014-03-21 | 2014-09-17 | 西北工业大学 | 连续式跨声速风洞实验段流场变雷诺数控制方法 |
JP2017049014A (ja) * | 2015-08-31 | 2017-03-09 | 公益財団法人鉄道総合技術研究所 | 騒音測定方法及び騒音算出装置 |
CN105716826A (zh) * | 2016-02-18 | 2016-06-29 | 江西洪都航空工业集团有限责任公司 | 一种零升阻力系数的雷诺数效应修正方法 |
CN105973563A (zh) * | 2016-05-12 | 2016-09-28 | 中国科学院力学研究所 | 一种强激光与超声速风洞联合破坏效应的试验系统及方法 |
CN110160734A (zh) * | 2019-05-23 | 2019-08-23 | 华中科技大学 | 基于多孔介质的风洞喷管、主动降噪装置及方法 |
CN110763420A (zh) * | 2019-10-23 | 2020-02-07 | 中国空气动力研究与发展中心高速空气动力研究所 | 低温风洞流场自抗扰控制方法 |
CN111459175A (zh) * | 2019-12-18 | 2020-07-28 | 北京航空航天大学 | 一种基于l1自适应控制器的尾座式无人机轨迹跟踪容错控制方法 |
CN111351628A (zh) * | 2020-03-27 | 2020-06-30 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种高超声速风洞连续变攻角试验数据处理方法 |
CN112304563A (zh) * | 2020-10-30 | 2021-02-02 | 中国空气动力研究与发展中心超高速空气动力研究所 | 研究转捩对高超声速飞行器气动特性影响的风洞试验方法 |
CN112525474A (zh) * | 2020-12-22 | 2021-03-19 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种获得跨声速风洞堵塞干扰因子的方法 |
CN113267314A (zh) * | 2021-04-25 | 2021-08-17 | 中国航天空气动力技术研究院 | 一种暂冲式风洞的超声速流场总压控制系统 |
Non-Patent Citations (4)
Title |
---|
NAGAI S.: "Validity of Statistical Uncertainty in Comparison of Winged Vehicle Force Data at Hypersonic Wind Tunnels" * |
QIU, Y: "Modeling the mean wind loads on cylindrical roofs with consideration of the Reynolds number effect in uniform flow with low turbulence" * |
谢飞: "高超声速风洞试验模型底压测量方法研究" * |
谢飞等: "高超声速风洞变雷诺数试验技术研究" * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115615655A (zh) * | 2022-11-15 | 2023-01-17 | 中国航空工业集团公司沈阳空气动力研究所 | 一种连续式跨声速风洞试验数据干扰修正方法 |
CN115615655B (zh) * | 2022-11-15 | 2023-03-31 | 中国航空工业集团公司沈阳空气动力研究所 | 一种连续式跨声速风洞试验数据干扰修正方法 |
CN116929703A (zh) * | 2023-09-18 | 2023-10-24 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种考虑堵塞效应的低温风洞马赫数确定方法及应用 |
CN116929703B (zh) * | 2023-09-18 | 2023-11-21 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种考虑堵塞效应的低温风洞马赫数确定方法及应用 |
Also Published As
Publication number | Publication date |
---|---|
CN114018531B (zh) | 2023-03-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114018531B (zh) | 一种高超声速风洞连续变总压的试验数据处理方法 | |
EP2711678A2 (en) | Correction of pressure signals measured during supersonic wind tunnel testing | |
Pisasale et al. | A novel method for extending the calibration range of five-hole probe for highly three-dimensional flows | |
CN108255786B (zh) | 一种称重结果的干扰补偿计算方法及系统 | |
CN108318491B (zh) | 一种基于频谱曲率分析的织物缺陷检测方法 | |
CN117093828B (zh) | 一种液压缸缸体密封性检测方法及系统 | |
CN108168680B (zh) | 一种动态称重滤波方法及其系统 | |
CN112213092A (zh) | 一种采用惯性法测试拱桥吊杆内力增量的测量方法 | |
CN113848963B (zh) | 一种飞行控制系统的控制律参数设计方法 | |
CN113029512B (zh) | 用于暂冲式跨声速风洞的空气流场模型校正控制方法 | |
Xue et al. | Effects of tube system and data correction for fluctuating pressure test in wind tunnel | |
CN112762962A (zh) | 基于温度滞后模型的微机电系统加速度计零偏补偿方法 | |
CN117330278A (zh) | 用于全桥半导体天平的温度补偿信号调理方法 | |
CN116295740A (zh) | 一种基于自适应变分模态分解的信号去噪重构方法 | |
CN109471192B (zh) | 一种全自动重力测试仪高精度动态数据处理方法 | |
CN110031087B (zh) | 一种用于测量飞机舱内噪声的方法 | |
CN106525223A (zh) | 一种齿轮传动装置异响的下线检测方法 | |
CN117147022A (zh) | 力传感器非线性补偿方法及系统 | |
CN112461432A (zh) | 一种实现大气数据测量长期稳定的在线修正方法 | |
CN102297745B (zh) | 基于gm(0,2)模型的多维力传感器的静态解耦方法 | |
CN114813096B (zh) | 一种用于高超声速风洞的多层烧结网试验选型方法 | |
CN106873363B (zh) | 一种飞行器迎角信号的建模方法 | |
CN103076035B (zh) | 一种基于双支持向量机的传感器测量方法 | |
CN104331862A (zh) | 一种并联型分数阶零相位滤波器及其滤波方法 | |
CN114993606A (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 |