CN108414185A - 对称飞行器风洞试验数据零点误差处理方法 - Google Patents
对称飞行器风洞试验数据零点误差处理方法 Download PDFInfo
- Publication number
- CN108414185A CN108414185A CN201810129236.3A CN201810129236A CN108414185A CN 108414185 A CN108414185 A CN 108414185A CN 201810129236 A CN201810129236 A CN 201810129236A CN 108414185 A CN108414185 A CN 108414185A
- Authority
- CN
- China
- Prior art keywords
- wind tunnel
- zero
- tunnel test
- test data
- positive
- 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
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/06—Measuring arrangements specially adapted for aerodynamic testing
Landscapes
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明提供了一种对称飞行器风洞试验数据零点误差处理方法:将具有对称性的基本飞行姿态选定为风洞试验数据零点的标准状态;确定合成攻角为零时飞行器横截面内空气动力合力和合力矩的方向和大小;在旋转体轴系下对飞行器横截面内合力和合力矩进行分解,获得相应的力和力矩的新值;通过代数平均,获得轴向力和滚转力矩零点的新值;获得新旧值的差量;将随合成攻角变化的旧数据加上前述差量获得新数据,实现曲线的整体平移,以供配套使用。本发明从零点状态是同一个的物理实际出发,创新提出零点力和力矩归一化的思路,建立一种对称飞行器风洞试验数据零点误差处理方法,可全面考虑风洞试验模型安装误差和风洞流场均匀度误差,提高了风洞试验数据的精度,为制导、控制和自动驾驶仪的飞行仿真和特性评估解决了零点误差的干扰问题。
Description
技术领域
本发明涉及试验数据误差处理领域,特别涉及一种对称飞行器风洞试验数据零点误差处理方法。
背景技术
飞行器所受空气动力和力矩主要有计算和风洞试验两种获取手段。任何试验都存在试验误差,风洞试验也不例外。风洞试验从设计到实施已经有一套完备的模拟准则和多个相关的标准、规范来指导工作,同时其测量系统的误差也已形成成熟的修正方法。但是,对于风洞试验模型安装误差和风洞流场均匀度误差等受人为和时间发展的误差来源还没有形成规范的修正方法,不同气动设计师采用不同的方法处理,特别是合成攻角为零时的零点误差。
风洞试验获得的飞行器空气动力和力矩,最终以各种形式提供给制导、控制和自动驾驶仪进行飞行仿真和特性评估。风洞试验误差的存在,导致零点时虽然物理状态一样,但结果却不一样的现象,在仿真中引起零点振荡,易给出系统不稳定的错误结论。
现有技术中,一般统一将所有情况零点的值强行置为零,且不整体平移数据,改变了零点附件物理量对合成攻角的导数,或者不进行零点误差的处理。对上述虚假的零点振荡现象没有本质改善,仍然存在零点误差的干扰问题。
发明内容
针对现有技术中的缺陷,本发明的目的是提供一种对称飞行器风洞试验数据零点误差处理方法。
根据本发明提供的一种对称飞行器风洞试验数据零点误差处理方法,包括如下步骤:
步骤S1、将飞行器的基本飞行姿态选定为风洞试验数据零点的标准状态,其中,飞行器的基本飞行姿态至少具有左右对称性;
步骤S2、利用标准状态的对称性,确定合成攻角为零时飞行器横截面内的空气动力合力、空气动力合力矩的方向和大小;
步骤S3、在旋转体轴系OXYZ下,对飞行器OYZ坐标面内的空气动力合力、空气动力合力矩进行分解,获得相应滚转角下零点力和零点力矩四分量的新值;
步骤S4、通过代数平均,获得轴向力和滚转力矩零点的新值;
步骤S5、根据步骤S3和步骤S4获得的新值,获得所有零点力和零点力矩新旧值的差量;其中,步骤S3获得的新值是指相应滚转角下零点力和零点力矩四分量的新值,步骤S4获得的新值是指轴向力和滚转力矩零点的新值;
步骤S6、将随合成攻角变化的一套旧数据加上前述差量,获得一套新数据,实现曲线的整体平移,以供配套使用。
优选地,由滚转角决定旋转体轴系OXYZ相对于飞行器机体的方位,在旋转体轴系的OYZ坐标面内,不考虑轴向力和滚转力矩,飞行器左右对称的部分因产生相互抵消的侧向力和侧向力矩而进行置零处理,飞行器在OYZ坐标面的空气动力合力空气动力合力矩的方向和大小即分别为标准状态的法向力、俯仰力矩的方向和大小,公式如下:
其中
旋转体轴系OXYZ:风洞试验时,飞行器不动,空气来流,
原点O:飞行器质心,
X轴正向:沿体轴指向来流,单位矢量
Y轴正向:来流速度在飞行器横截面投影的方向,单位矢量
Z轴正向:按照右手法则确定,单位矢量
α:合成攻角,来流与X轴的夹角,取正值;
FZ:侧向力,正方向为Z轴正向;
MY:质心处偏航力矩,正方向为Y轴正向;
FY:法向力,正方向为Y轴正向;
MZ:质心处俯仰力矩,正方向为Z轴正向;
下标s:标准状态;
下标n:新值。
优选地,所述在旋转体轴系下,对飞行器OYZ坐标面内的空气动力合力空气动力合力矩进行分解,获得以下力和力矩的新值:
其中
Φ:滚转角,4#舵与Y轴的夹角,逆时针为正。
优选地,通过代数平均,获得轴向力和滚转力矩零点的新值,公式如下:
其中
FX:轴向力,正方向为X轴负向;
MX:滚转力矩,正方向为X轴正向;
m:试验时滚转角的总数;
下标l:试验时滚转角的序号,取值为1到m的正整数。
优选地,通过减法运算,获得零点所有力和力矩分量新、旧值的差量:
其中
前缀Δ:差量;
下标o:旧值。
优选地,在所述步骤6中,将随合成攻角变化的一套旧数据加上前述差量,获得一套新数据,实现曲线的整体平移,以供配套使用:
优选地,所述标准状态,还具有上下对称性。
优选地,所述标准状态,还具有轴对称性。
与现有技术相比,本发明具有如下的有益效果:
本发明从零点状态是同一个的物理实际出发,创新提出零点力和力矩归一化的思路,可全面考虑风洞试验模型安装误差和风洞流场均匀度误差,提高了风洞试验数据的精度,为制导、控制和自动驾驶仪的飞行仿真和特性评估解决了零点误差的干扰问题。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1是本发明实施例中对称飞行器风洞试验数据零点误差处理方法的坐标系和滚转角定义的示意图,为后视图。
图中:弹体10,尾舵20。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变化和改进。这些都属于本发明的保护范围。
针对物理状态一样但结果却不一样的现象,本发明的目的在于提供一种对称飞行器风洞试验数据零点误差处理方法,抑制仿真中虚假系统不稳定的出现,以解决零点误差的干扰问题,提高风洞试验数据的精度。
本发明一种对称飞行器风洞试验数据零点误差处理方法包括以下步骤:
(1)飞行器的基本飞行姿态一般具有对称性,将其选定为风洞试验数据零点的标准状态;
(2)利用标准状态的对称性,确定合成攻角为零时飞行器横截面内的空气动力合力和合力矩的方向和大小(不考虑轴向力和滚转力矩);
(3)根据需要,在旋转体轴系OXYZ下,对飞行器OYZ坐标面内的合力和合力矩进行分解,获得相应滚转角下零点力和力矩的分量新值;
(4)通过代数平均,获得轴向力和滚转力矩零点的新值;
(5)通过减法运算,获得零点所有力和力矩分量新、旧值的差量;
(6)将随合成攻角变化的一套旧数据加上前述差量,获得一套新数据,实现曲线的整体平移,以供配套使用。
现以一具体实施例详细说明本发明一种对称飞行器风洞试验数据零点误差处理方法:
图1所示,是对称飞行器风洞试验数据零点误差处理方法的坐标系和舵偏的定义,包括旋转体轴系OXYZ、机体10和尾舵20,其中有4个尾舵,分别为1#舵、2#舵、3#舵、4#舵,即一号舵、二号舵、三号舵、四号舵。
本实施例中涉及的定义如下:
坐标系
旋转体轴系OXYZ:风洞试验时,飞行器不动,空气来流,
原点O:飞行器质心,
X轴正向:沿体轴指向来流,单位矢量
Y轴正向:来流速度在飞行器横截面投影的方向,单位矢量
Z轴正向:按照右手法则确定,单位矢量
符号
α:合成攻角,来流与X轴的夹角,取正值;
Φ:滚转角,4#舵与Y轴的夹角,逆时针为正;
FY:法向力,正方向为Y轴正向;
Mz:质心处俯仰力矩,正方向为Z轴正向;
FZ:侧向力,正方向为Z轴正向;
MY:质心处偏航力矩,正方向为Y轴正向;
FX:轴向力,正方向为X轴负向;
MX:滚转力矩,正方向为X轴正向;
m:试验时滚转角的总数;
前缀
Δ:某量的差量;
下标
l:试验时某个滚转角的序号,取值为1到m的正整数;
n:某量的新值;
o:某量的旧值;
s:标准状态。
本实施例的一种对称飞行器风洞试验数据零点误差处理方法包括以下步骤:
(1)飞行器的基本飞行姿态一般至少具有左右对称性,将其选定为风洞试验数据零点的标准状态;
(2)利用标准状态的对称性,确定合成攻角为零时飞行器OYZ坐标面内的空气动力合力和合力矩的方向和大小;
不考虑轴向力和滚转力矩时,飞行器左右对称的部分产生相互抵消的侧向力和力矩,进行置零处理,飞行器OYZ坐标面内的合力和合力矩的方向和大小即为标准状态的法向力和俯仰力矩的方向和大小,公式如下:
(3)根据需要,在旋转体轴系下对飞行器OYZ坐标面的合力和合力矩进行分解,获得以下力和力矩的新分量,分解公式如下:
(4)通过代数平均,获得轴向力和滚转力矩零点的新值,公式如下:
(5)通过减法运算,获得零点所有力和力矩分量新、旧值的差量;
(6)将随合成攻角变化的一套旧数据加上前述差量,获得一套新数据,实现曲线的整体平移,以供配套使用:
本发明一种对称飞行器风洞试验数据零点误差处理方法,从零点状态是同一个的物理实际出发,创新提出零点力和力矩归一化的思路,可全面考虑风洞试验模型安装误差和风洞流场均匀度误差,提高了风洞试验数据的精度,为制导、控制和自动驾驶仪的飞行仿真和特性评估解决了零点误差的干扰问题。本发明可应用于基本飞行状态是对称的所有飞行器的风洞试验数据零点误差处理上。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变化或修改,这并不影响本发明的实质内容。在不冲突的情况下,本申请的实施例和实施例中的特征可以任意相互组合。
Claims (10)
1.一种对称飞行器风洞试验数据零点误差处理方法,其特征在于,包括如下步骤:
步骤S1、将飞行器的基本飞行姿态选定为风洞试验数据零点的标准状态,其中,飞行器的基本飞行姿态至少具有左右对称性;
步骤S2、利用标准状态的对称性,确定合成攻角为零时飞行器横截面内的空气动力合力、空气动力合力矩的方向和大小;
步骤S3、在旋转体轴系OXYZ下,对飞行器OYZ坐标面内的空气动力合力、空气动力合力矩进行分解,获得相应滚转角下零点力和零点力矩四分量的新值;
步骤S4、通过代数平均,获得轴向力和滚转力矩零点的新值;
步骤S5、根据步骤S3和步骤S4获得的新值,获得所有零点力和零点力矩新旧值的差量;
步骤S6、将随合成攻角变化的一套旧数据加上所述差量,获得一套新数据,实现曲线的整体平移,以供配套使用。
2.根据权利要求1所述的对称飞行器风洞试验数据零点误差处理方法,其特征在于,由滚转角决定旋转体轴系OXYZ相对于飞行器机体的方位,在旋转体轴系的OYZ坐标面内,不考虑轴向力和滚转力矩,飞行器左右对称的部分因产生相互抵消的侧向力和侧向力矩而进行置零处理,飞行器在OYZ坐标面内的空气动力合力空气动力合力矩的方向和大小即分别为标准状态的法向力、俯仰力矩的方向和大小,公式如下:
α=0
α=0
α=0
α=0
其中
旋转体轴系OXYZ:风洞试验时,飞行器不动,空气来流,
原点O:飞行器质心,
X轴正向:沿体轴指向来流,单位矢量
Y轴正向:来流速度在飞行器横截面投影的方向,单位矢量
Z轴正向:按照右手法则确定,单位矢量
α:合成攻角,来流与X轴的夹角,取正值;
FZ:侧向力,正方向为Z轴正向;
MY:质心处偏航力矩,正方向为Y轴正向;
FY:法向力,正方向为Y轴正向;
MZ:质心处俯仰力矩,正方向为Z轴正向;
下标s:标准状态;
下标n:新值。
3.根据权利要求1所述的对称飞行器风洞试验数据零点误差处理方法,其特征在于,所述在旋转体轴系下,对飞行器OYZ坐标面内的空气动力合力空气动力合力矩进行分解,获得以下力和力矩的新值:
α=0
α=0
α=0
α=0
其中
Φ:滚转角,4#舵与Y轴的夹角,逆时针为正。
4.根据权利要求1所述的对称飞行器风洞试验数据零点误差处理方法,其特征在于,通过代数平均,获得轴向力和滚转力矩零点的新值,公式如下:
其中
FX:轴向力,正方向为X轴负向;
MX:滚转力矩,正方向为X轴正向;
m:试验时滚转角的总数;
下标l:试验时滚转角的序号,取值为1到m的正整数。
5.根据权利要求1所述的对称飞行器风洞试验数据零点误差处理方法,其特征在于,通过减法运算,获得零点所有力和力矩分量新、旧值的差量:
α=0
α=0
α=0
α=0
α=0
α=0
其中
前缀Δ:差量;
下标o:旧值。
6.根据权利要求1所述的对称飞行器风洞试验数据零点误差处理方法,其特征在于,在所述步骤6中,将随合成攻角变化的一套旧数据加上前述差量,获得一套新数据,实现曲线的整体平移,以供配套使用:
7.根据权利要求1所述的对称飞行器风洞试验数据零点误差处理方法,其特征在于,不同滚转角下零点的飞行器构型一致。
8.根据权利要求1所述的对称飞行器风洞试验数据零点误差处理方法,其特征在于,零点是指合成攻角为零的状态点。
9.根据权利要求1所述的对称飞行器风洞试验数据零点误差处理方法,其特征在于,所述标准状态,还具有上下对称性。
10.根据权利要求9所述的对称飞行器风洞试验数据零点误差处理方法,其特征在于,所述标准状态,还具有轴对称性。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810129236.3A CN108414185B (zh) | 2018-02-08 | 2018-02-08 | 对称飞行器风洞试验数据零点误差处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810129236.3A CN108414185B (zh) | 2018-02-08 | 2018-02-08 | 对称飞行器风洞试验数据零点误差处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108414185A true CN108414185A (zh) | 2018-08-17 |
CN108414185B CN108414185B (zh) | 2020-01-31 |
Family
ID=63128147
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810129236.3A Active CN108414185B (zh) | 2018-02-08 | 2018-02-08 | 对称飞行器风洞试验数据零点误差处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108414185B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112484955A (zh) * | 2020-12-07 | 2021-03-12 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种飞行器滚转控制的风洞仿真方法 |
CN114264404A (zh) * | 2021-12-28 | 2022-04-01 | 中国航空工业集团公司西安飞机设计研究所 | 一种结构强度试验用测力传感器的调零方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004347389A (ja) * | 2003-05-21 | 2004-12-09 | Central Glass Co Ltd | 風圧測定方法 |
CN102095430A (zh) * | 2010-11-18 | 2011-06-15 | 合肥工业大学 | 基于阶跃响应的传感器动态误差频域修正技术 |
CN103335814A (zh) * | 2013-06-09 | 2013-10-02 | 电子科技大学 | 风洞中实验模型倾角测量误差数据修正系统及修正方法 |
KR20140127647A (ko) * | 2013-04-25 | 2014-11-04 | (주)수도프리미엄엔지니어링 | 가상 풍동의 운전 방법 |
CN105136422A (zh) * | 2015-09-10 | 2015-12-09 | 中国航天空气动力技术研究院 | 风洞试验中修正飞行器模型侧滑弹性角的方法 |
CN105571811A (zh) * | 2015-12-22 | 2016-05-11 | 中国航天空气动力技术研究院 | 测量风洞实验中飞行器实际攻角值的方法 |
CN105628325A (zh) * | 2014-10-29 | 2016-06-01 | 北京临近空间飞行器系统工程研究所 | 一种锥型面气动压力场实时高精度获取方法 |
CN105716826A (zh) * | 2016-02-18 | 2016-06-29 | 江西洪都航空工业集团有限责任公司 | 一种零升阻力系数的雷诺数效应修正方法 |
-
2018
- 2018-02-08 CN CN201810129236.3A patent/CN108414185B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004347389A (ja) * | 2003-05-21 | 2004-12-09 | Central Glass Co Ltd | 風圧測定方法 |
CN102095430A (zh) * | 2010-11-18 | 2011-06-15 | 合肥工业大学 | 基于阶跃响应的传感器动态误差频域修正技术 |
KR20140127647A (ko) * | 2013-04-25 | 2014-11-04 | (주)수도프리미엄엔지니어링 | 가상 풍동의 운전 방법 |
CN103335814A (zh) * | 2013-06-09 | 2013-10-02 | 电子科技大学 | 风洞中实验模型倾角测量误差数据修正系统及修正方法 |
CN105628325A (zh) * | 2014-10-29 | 2016-06-01 | 北京临近空间飞行器系统工程研究所 | 一种锥型面气动压力场实时高精度获取方法 |
CN105136422A (zh) * | 2015-09-10 | 2015-12-09 | 中国航天空气动力技术研究院 | 风洞试验中修正飞行器模型侧滑弹性角的方法 |
CN105571811A (zh) * | 2015-12-22 | 2016-05-11 | 中国航天空气动力技术研究院 | 测量风洞实验中飞行器实际攻角值的方法 |
CN105716826A (zh) * | 2016-02-18 | 2016-06-29 | 江西洪都航空工业集团有限责任公司 | 一种零升阻力系数的雷诺数效应修正方法 |
Non-Patent Citations (3)
Title |
---|
W. GAWRONSKI 等: "Antenna Mean Wind Torques: A Comparison of Field and Wind-Tunnel Data", 《IEEE ANTENNAS AND PROPAGATION MAGAZINE》 * |
吕治国 等: "激波风洞测力试验不确定度分析", 《江汉大学学报(自然科学版)》 * |
恽起麟: "风洞实验数据误差分析", 《气动实验与测量控制》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112484955A (zh) * | 2020-12-07 | 2021-03-12 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种飞行器滚转控制的风洞仿真方法 |
CN114264404A (zh) * | 2021-12-28 | 2022-04-01 | 中国航空工业集团公司西安飞机设计研究所 | 一种结构强度试验用测力传感器的调零方法 |
CN114264404B (zh) * | 2021-12-28 | 2024-02-23 | 中国航空工业集团公司西安飞机设计研究所 | 一种结构强度试验用测力传感器的调零方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108414185B (zh) | 2020-01-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109324636B (zh) | 基于二阶一致性和自抗扰的多四旋翼主从式协同编队控制方法 | |
CN106681348B (zh) | 考虑全捷联导引头视场约束的制导控制一体化设计方法 | |
CN104960674B (zh) | 一种运动目标的指向跟踪控制方法 | |
CN103926931B (zh) | 轴对称高速飞行器运动特征综合识别方法 | |
CN106228014B (zh) | 一种导弹气动系数的获取方法 | |
CN109612676B (zh) | 基于飞行试验数据的气动参数反算方法 | |
CN107844128B (zh) | 一种基于复合比例导引的高超声速飞行器巡航段制导方法 | |
CN109635494A (zh) | 一种飞行试验与地面仿真气动力数据综合建模方法 | |
CN108414185A (zh) | 对称飞行器风洞试验数据零点误差处理方法 | |
Bao et al. | Integrated guidance and control for hypersonic morphing missile based on variable span auxiliary control | |
CN108180910A (zh) | 一种基于气动参数不确定的飞行器快速高精度制导方法 | |
CN115639830B (zh) | 一种空地智能体协同编队控制系统及其编队控制方法 | |
CN109063391A (zh) | 旋转条件下的动导数计算检测方法及动导数风洞试验方法 | |
CN109724626A (zh) | 一种针对极区传递对准动态挠曲杆臂效应的模型补偿方法 | |
CN113253617A (zh) | 用于四旋翼无人机的在线自适应控制方法 | |
CN115366109A (zh) | 一种旋翼飞行机械臂复合分层抗干扰方法 | |
CN114721266B (zh) | 飞机舵面结构性缺失故障情况下的自适应重构控制方法 | |
CN106774385B (zh) | 一种采用自适应变结构的飞艇定点悬停控制方法 | |
CN106570242B (zh) | 低空大动压整体式整流罩高速分离流固耦合仿真方法 | |
CN108427428A (zh) | 基于改进迭代算法的自适应滑模变结构航天器姿态控制方法 | |
CN110414110B (zh) | 一种用于飞行失速状态下的飞机受力仿真方法 | |
Charters | The linearized equations of motion underlying the dynamic stability of aircraft, spinning projectiles, and symmetrical missiles | |
CN104613984B (zh) | 临近空间飞行器传递对准模型不确定性的鲁棒滤波方法 | |
Bogos et al. | Similarity criteria for “full” and “scale” aircraft on the lateral stabilty analysis | |
CN113110428A (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 |