CN103234729A - 风洞常规测力试验中气动刚度与气动阻尼的视频测量方法 - Google Patents
风洞常规测力试验中气动刚度与气动阻尼的视频测量方法 Download PDFInfo
- Publication number
- CN103234729A CN103234729A CN2013101193990A CN201310119399A CN103234729A CN 103234729 A CN103234729 A CN 103234729A CN 2013101193990 A CN2013101193990 A CN 2013101193990A CN 201310119399 A CN201310119399 A CN 201310119399A CN 103234729 A CN103234729 A CN 103234729A
- Authority
- CN
- China
- Prior art keywords
- angle
- attack
- ladder
- model
- wind tunnel
- 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
Landscapes
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明公开了一种风洞常规测力试验中气动刚度与气动阻尼的视频测量方法,利用常规测力试验采用的迎角阶梯变化的风洞试验方式,通过迎角视频测量***测得的迎角数据与天平测得的气动力数据,计算出该测力试验参数下,各迎角阶梯的气动阻尼和气动刚度,从而在风洞常规测力试验中,测得给定迎角阶梯气动力的同时,获得了该试验马赫数与雷诺数下的气动阻尼和气动刚度。本发明的积极效果是:实现了在常规测力试验中,采用常规测力试验模型,获得俯仰自由度的气动阻尼和气动刚度,扩展了常规测力试验获得的气动数据;同时,由于本发明采用常规测力试验模型,测得俯仰自由度的气动阻尼和气动刚度,较颤振试验模型或动导数试验模型研制成本低廉得多。
Description
技术领域
本发明涉及一种风洞试验或飞行器地面试验方法,尤其涉及一种风洞常规测力试验中气动刚度与气动阻尼的视频测量方法。
背景技术
颤振是导致飞行器机毁人亡的破坏性振动,故在飞行包线内绝不允许出现颤振现象。
目前,公认的颤振驱动机理分为:阻尼驱动型和刚度驱动型。其中阻尼驱动是指随着飞行速度的增加飞行器***的总阻尼将由正转负,气动负阻尼抵消了结构阻尼,从而使***能量不能得到耗散,引起飞行器结构振动发散,如分离流颤振;刚度驱动型则是指气流带来的气动刚度效应改变了飞行器***扭转振动和弯曲振动的频率,使弯曲频率增大,扭转频率减小,从而在颤振临界点耦合成统一的颤振频率,驱动结构振动发散。
因此,气动阻尼和气动刚度对于飞行器设计至关重要。
尽管计算流体动力学(CFD)已迅速发展,但对分离流与旋涡、激波/边界层、剪切层等复杂流动现象,CFD尚不能建立较精确的物理模型,在相当长的时间内尚难高精度模拟飞行器的气动阻尼和气动刚度。
所以,目前需要通过颤振试验或动导数试验才能得到气动阻尼和气动刚度,上述试验模型除了外形模拟外,对于颤振试验模型还必须刚度模拟和惯性模拟,而动导数试验需要综合考虑弹性铰链的刚度和试验模型的重量及质量分布,故试验模型研制困难,制作周期长,成本高。
发明内容
为了克服现有技术的上述缺点,本发明提供了一种风洞常规测力试验中气动刚度与气动阻尼的视频测量方法,利用常规测力试验采用的迎角阶梯变化的风洞试验方式,通过迎角视频测量***测得的迎角数据与天平测得的气动力数据,计算出该测力试验参数下,各迎角阶梯的气动阻尼和气动刚度,从而在风洞常规测力试验中,测得给定迎角阶梯气动力的同时,获得了该试验马赫数与雷诺数下的气动阻尼和气动刚度。
本发明解决其技术问题所采用的技术方案是:一种风洞常规测力试验中气动刚度与气动阻尼的视频测量方法,包括如下步骤:
步骤一、测量各阶梯上风洞模型的迎角:利用安装在风洞试验段侧壁面上的迎角视频测量***,通过测量风洞模型上的特征点三维坐标,计算出各阶梯风洞模型迎角;
步骤二、计算风洞模型、测力天平与模型支撑机构构成的***在俯仰自由度的一阶振型的频率ω0、弹性铰链常数ks、转动惯量I和结构阻尼cs:
①在风洞未吹风时,对风洞模型施加初始位移或施加初始速度扰动后释放,然后测得模型迎角自由振动数据α(t);
②对α(t)进行傅里叶变换,提取俯仰自由度的一阶振型频率ω0及阻尼比ξ0;
③按如下公式计算振幅对数衰减率δ:
式中Q0为迎角自由振动时初始振幅值,Qn为第n个周期后的迎角振幅值;
④计算弹性铰链常数ks:
对风洞模型施加已知的俯仰力矩Mz,同时测出模型迎角变化量Δα,则
⑤计算俯仰自由度转动惯量I,计算式为:
⑥计算不吹风时的阻尼系数q0和俯仰自由度的结构阻尼cs,计算式为:
cs=2Iq0
步骤三、计算各个阶梯上风洞模型的俯仰气动阻尼和气动刚度:
①测量模型在各个阶梯的迎角数据;
②求解各个阶梯的俯仰气动阻尼和气动刚度:
设迎角视频测量***测得模型在第i个阶梯的迎角数据为αi(t),测力天平在第i个阶梯测得的升力和俯仰力矩分别记为Yi和则
式中Yi为第i个阶梯的升力,ky为升力Yi对应的迎角弹性铰链常数,通过天平校准获得;
式中,ωi为第i个阶梯的迎角数据αi(t)的一阶振型频率。
与现有技术相比,本发明的积极效果是:与传统常规测力试验、颤振试验和动导数试验不同,本发明在常规测力的阶梯式风洞试验方式中,通过迎角视频测量***,测得模型在吹风与未吹风时模型、天平与支杆***俯仰自由度的迎角测量数据,并对迎角测量数据进行傅里叶变换,从俯仰自由度一阶振型的频率变化,计算出风洞吹风状态下(即给定)的气动阻尼和气动刚度,从而实现了在常规测力试验中,采用常规测力试验模型,获得俯仰自由度的气动阻尼和气动刚度,扩展了常规测力试验获得的气动数据;同时,由于本发明采用常规测力试验模型,测得俯仰自由度的气动阻尼和气动刚度,较颤振试验模型或动导数试验模型研制成本低廉得多。
附图说明
本发明将通过例子并参照附图的方式说明,其中:
图1是风洞模型试验段的侧视图;
图2是图1的A-A视图;
图3是模型迎角的计算示意图。
具体实施方式
如图1、图2所示,风洞试验段包括:两端分别固定在风洞试验段上壁面6和风洞试验段下壁面7的模型支撑机构5;风洞模型1固定在测力天平3上、测力天平3固定在模型支撑机构5上,在风洞模型1设置有标记点2;在风洞试验段侧壁面8上安装有迎角视频测量***9。
一种风洞常规测力试验中气动刚度与气动阻尼的视频测量方法,包括以下步骤:
步骤一、测量各阶梯上风洞模型的迎角:
利用安装在风洞试验段侧壁面上的迎角视频测量***9,通过测量风洞模型上的特征点三维坐标,计算出各阶梯风洞模型迎角,计算方法如下:
设模型的旋转中心为风洞坐标系的原点O,风洞模型上的标记点为P,当模型迎角、侧滑角和滚转角都为零时,P在平面XOZ上的投影为当模型迎角运动到第i个阶梯时,P在平面XOZ上的投影为如图3所示,则第i个阶梯模型迎角αi的计算方法如下:
步骤二、计算风洞模型1、测力天平3与模型支撑机构5构成的***在俯仰自由度的一阶振型的频率ω0、弹性铰链常数ks、转动惯量I和结构阻尼cs:
①在风洞未吹风时,对风洞模型1施加初始位移(或施加初始速度扰动后释放),然后运用迎角的视频测量技术,测得模型迎角自由振动数据α(t);
②对α(t)进行傅里叶变换,提取俯仰自由度的一阶振型频率ω0及阻尼比ξ0;
③按如下公式计算振幅对数衰减率δ:
式中Q0为迎角自由振动时初始振幅值,Qn为第n个周期后的迎角振幅值。
④计算弹性铰链常数ks:
对风洞模型施加已知的俯仰力矩Mz,同时测出模型迎角变化量Δα,则
⑤计算俯仰自由度转动惯量I,计算式为:
⑥计算不吹风时的阻尼系数q0和俯仰自由度的结构阻尼cs。计算式为:
cs=2Iq0
步骤三、计算各个阶梯该模型缩比尺度下的俯仰气动阻尼和气动刚度:
①测量模型在各个阶梯的迎角数据:
在风洞常规测力试验过程中,采用迎角阶梯变化的风洞试验方式,是在给定的气流速度(或马赫数Ma)下,预先给定模型迎角变化阶梯,当模型姿态达到给定的迎角后,迎角视频测量***9和测力天平3等进行数据采集,当数据采集完毕后再进入下一个预定的迎角,如此反复直到测完预定的全部迎角的数据为止。
②求解各个阶梯该模型缩比尺度下的俯仰气动阻尼和气动刚度:
式中ky为升力Yi对应的迎角弹性铰链常数,通过天平校准获得。
B)计算第i个阶梯的气动刚度
式中,ωi为第i个阶梯的迎角数据αi(t)的一阶振型频率,对αi(t)进行傅里叶变换,即可得到ωi。
Claims (2)
1.一种风洞常规测力试验中气动刚度与气动阻尼的视频测量方法,其特征在于:包括如下步骤:
步骤一、测量各阶梯上风洞模型的迎角:利用安装在风洞试验段侧壁面上的迎角视频测量***,通过测量风洞模型上的特征点三维坐标,计算出各阶梯风洞模型迎角;
步骤二、计算风洞模型、测力天平与模型支撑机构构成的***在俯仰自由度的一阶振型的频率ω0、弹性铰链常数ks、转动惯量I和结构阻尼cs:
①在风洞未吹风时,对风洞模型施加初始位移或施加初始速度扰动后释放,然后测得模型迎角自由振动数据α(t);
②对α(t)进行傅里叶变换,提取俯仰自由度的一阶振型频率ω0及阻尼比ξ0;
③按如下公式计算振幅对数衰减率δ:
式中Q0为迎角自由振动时初始振幅值,Qn为第n个周期后的迎角振幅值;
④计算弹性铰链常数ks:
对风洞模型施加已知的俯仰力矩Mz,同时测出模型迎角变化量Δα,则
⑤计算俯仰自由度转动惯量I,计算式为:
⑥计算不吹风时的阻尼系数q0和俯仰自由度的结构阻尼cs,计算式为:
cs=2Iq0
步骤三、计算各个阶梯上风洞模型的俯仰气动阻尼和气动刚度:
①测量模型在各个阶梯的迎角数据;
②求解各个阶梯的俯仰气动阻尼和气动刚度:
式中Yi为第i个阶梯的升力,ky为升力Yi对应的迎角弹性铰链常数,通过天平校准获得;
B)计算第i个阶梯的气动刚度
式中,ωi为第i个阶梯的迎角数据αi(t)的一阶振型频率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310119399.0A CN103234729B (zh) | 2013-04-08 | 2013-04-08 | 风洞常规测力试验中气动刚度与气动阻尼的视频测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310119399.0A CN103234729B (zh) | 2013-04-08 | 2013-04-08 | 风洞常规测力试验中气动刚度与气动阻尼的视频测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103234729A true CN103234729A (zh) | 2013-08-07 |
CN103234729B CN103234729B (zh) | 2015-04-08 |
Family
ID=48882780
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310119399.0A Active CN103234729B (zh) | 2013-04-08 | 2013-04-08 | 风洞常规测力试验中气动刚度与气动阻尼的视频测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103234729B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103604579A (zh) * | 2013-11-28 | 2014-02-26 | 陕西飞机工业(集团)有限公司 | 一种模型主动式气动力影像网格试验方法 |
CN105550383A (zh) * | 2014-10-29 | 2016-05-04 | 北京临近空间飞行器***工程研究所 | 一种非定常气动力测量试验***的设计方法 |
CN108195543A (zh) * | 2017-11-29 | 2018-06-22 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种飞机风洞模型颤振吹风试验*** |
CN108204879A (zh) * | 2017-12-29 | 2018-06-26 | 中国航天空气动力技术研究院 | 一种转动惯量的光学测量方法及*** |
CN109632251A (zh) * | 2018-12-25 | 2019-04-16 | 同济大学 | 一种大尺度风洞试验模型同步测力试验装置 |
CN110006616A (zh) * | 2019-04-15 | 2019-07-12 | 中国空气动力研究与发展中心高速空气动力研究所 | 高速风洞模型的准备平台及其使用方法以及制作方法 |
CN110132527A (zh) * | 2019-06-24 | 2019-08-16 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种基于天平信号的风洞试验中模型振动监测方法 |
CN115290295A (zh) * | 2022-10-09 | 2022-11-04 | 中国空气动力研究与发展中心高速空气动力研究所 | 高速风洞级间分离与网格测力试验控制***及应用方法 |
CN115655635A (zh) * | 2022-12-14 | 2023-01-31 | 中国空气动力研究与发展中心高速空气动力研究所 | 用于全机体自由度颤振或阵风试验的两自由度支撑*** |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060096364A1 (en) * | 2004-11-08 | 2006-05-11 | Carlson Mark A | Measurement of coupled aerodynamic stability and damping derivatives in a wind tunnel |
RU2375691C1 (ru) * | 2008-06-06 | 2009-12-10 | Российская Федерация, от имени которой выступает государственный заказчик - Министерство промышленности и торговли Российской Федерации (Минпромторг России) | Динамический успокоитель колебаний аэродинамической модели |
CN101629822A (zh) * | 2009-08-10 | 2010-01-20 | 张征宇 | 振动环境中的多相机动态摄影测量方法 |
CN102175416A (zh) * | 2011-01-28 | 2011-09-07 | 中国空气动力研究与发展中心高速空气动力研究所 | 风洞试验中模型姿态角测量的多相机动态标定方法 |
CN102680201A (zh) * | 2012-05-15 | 2012-09-19 | 空气动力学国家重点实验室 | 基于视频测量的抖振风洞试验方法 |
-
2013
- 2013-04-08 CN CN201310119399.0A patent/CN103234729B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060096364A1 (en) * | 2004-11-08 | 2006-05-11 | Carlson Mark A | Measurement of coupled aerodynamic stability and damping derivatives in a wind tunnel |
RU2375691C1 (ru) * | 2008-06-06 | 2009-12-10 | Российская Федерация, от имени которой выступает государственный заказчик - Министерство промышленности и торговли Российской Федерации (Минпромторг России) | Динамический успокоитель колебаний аэродинамической модели |
CN101629822A (zh) * | 2009-08-10 | 2010-01-20 | 张征宇 | 振动环境中的多相机动态摄影测量方法 |
CN102175416A (zh) * | 2011-01-28 | 2011-09-07 | 中国空气动力研究与发展中心高速空气动力研究所 | 风洞试验中模型姿态角测量的多相机动态标定方法 |
CN102680201A (zh) * | 2012-05-15 | 2012-09-19 | 空气动力学国家重点实验室 | 基于视频测量的抖振风洞试验方法 |
Non-Patent Citations (1)
Title |
---|
张征宇等: "风洞试验中模型迎角的视频测量及精度研究", 《实验流体力学》 * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103604579B (zh) * | 2013-11-28 | 2016-08-17 | 陕西飞机工业(集团)有限公司 | 一种模型主动式气动力影像网格试验方法 |
CN103604579A (zh) * | 2013-11-28 | 2014-02-26 | 陕西飞机工业(集团)有限公司 | 一种模型主动式气动力影像网格试验方法 |
CN105550383A (zh) * | 2014-10-29 | 2016-05-04 | 北京临近空间飞行器***工程研究所 | 一种非定常气动力测量试验***的设计方法 |
CN105550383B (zh) * | 2014-10-29 | 2018-05-18 | 北京临近空间飞行器***工程研究所 | 一种非定常气动力测量试验***的设计方法 |
CN108195543A (zh) * | 2017-11-29 | 2018-06-22 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种飞机风洞模型颤振吹风试验*** |
CN108204879B (zh) * | 2017-12-29 | 2019-11-29 | 中国航天空气动力技术研究院 | 一种转动惯量的光学测量方法及*** |
CN108204879A (zh) * | 2017-12-29 | 2018-06-26 | 中国航天空气动力技术研究院 | 一种转动惯量的光学测量方法及*** |
CN109632251A (zh) * | 2018-12-25 | 2019-04-16 | 同济大学 | 一种大尺度风洞试验模型同步测力试验装置 |
CN109632251B (zh) * | 2018-12-25 | 2019-10-18 | 同济大学 | 一种大尺度风洞试验模型同步测力试验装置 |
CN110006616A (zh) * | 2019-04-15 | 2019-07-12 | 中国空气动力研究与发展中心高速空气动力研究所 | 高速风洞模型的准备平台及其使用方法以及制作方法 |
CN110132527A (zh) * | 2019-06-24 | 2019-08-16 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种基于天平信号的风洞试验中模型振动监测方法 |
CN110132527B (zh) * | 2019-06-24 | 2020-11-20 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种基于天平信号的风洞试验中模型振动监测方法 |
CN115290295A (zh) * | 2022-10-09 | 2022-11-04 | 中国空气动力研究与发展中心高速空气动力研究所 | 高速风洞级间分离与网格测力试验控制***及应用方法 |
CN115290295B (zh) * | 2022-10-09 | 2022-12-23 | 中国空气动力研究与发展中心高速空气动力研究所 | 高速风洞级间分离与网格测力试验控制***及应用方法 |
CN115655635A (zh) * | 2022-12-14 | 2023-01-31 | 中国空气动力研究与发展中心高速空气动力研究所 | 用于全机体自由度颤振或阵风试验的两自由度支撑*** |
CN115655635B (zh) * | 2022-12-14 | 2023-03-07 | 中国空气动力研究与发展中心高速空气动力研究所 | 用于全机体自由度颤振或阵风试验的两自由度支撑*** |
Also Published As
Publication number | Publication date |
---|---|
CN103234729B (zh) | 2015-04-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103234729B (zh) | 风洞常规测力试验中气动刚度与气动阻尼的视频测量方法 | |
CN110928201B (zh) | 一种飞机航电***半物理试验方法及*** | |
CN102680201B (zh) | 基于视频测量的抖振风洞试验方法 | |
Avin et al. | Experimental aeroelastic benchmark of a very flexible wing | |
Sahoo et al. | An accelerometer balance system for measurement of aerodynamic force coefficients over blunt bodies in a hypersonic shock tunnel | |
CN105136423B (zh) | 考虑摩擦力的自由振动动导数试验的数据分析方法 | |
CN103424225B (zh) | 一种测试转动部件动静不平衡量的方法 | |
CN103400035B (zh) | 一种高可信度快速预测飞行器滚转动导数的方法 | |
CN106950853A (zh) | 一种球形贮箱微重力环境下液体晃动的建模方法 | |
CN110160740A (zh) | 一种一体化冲击气动力测量*** | |
Zhao et al. | Structural and aeroelastic design, analysis, and experiments of inflatable airborne wings | |
Saravanan et al. | Aerodynamic force measurement using 3-component accelerometer force balance system in a hypersonic shock tunnel | |
DURHAM et al. | A status report on a model for benchmark active controls testing | |
CN101526352B (zh) | 运动平台上重力方向确定方法 | |
Raza et al. | Applicability of 3D printed fighter aircraft model for subsonic wind tunnel | |
Gillis et al. | Wing-on and Wing-off Longitudinal Characteristics of an Airplane Configuration Having an Thin Unswept Tapered Wing of Aspect Ratio 3, as Obtained from Rocket-propelled Models at Mach Numbers from 0.8 to 1.4 | |
CN106248281A (zh) | 一种平面扭转式微冲量测量方法 | |
Bansal et al. | Stochastic variations in aerodynamic influence coefficients (aics) on flutter prediction of a generic wing | |
Mitchell | Assessment of the accuracy of gust response calculations by comparison with experiments | |
Guglieri et al. | Dynamic stability derivatives evaluation in a low-speed wind tunnel | |
CN116296237B (zh) | 一种低速风洞大型运载火箭竖立风载试验方法 | |
Tyler et al. | Rapid technology focused experimental and computational aerodynamic investigation of a strike tanker | |
Malik | Experimental and computational analysis of a 3d printed wing structure | |
Sanmugadas et al. | Development of a low-cost flutter test bed with an EPS foam model for preliminary wing design | |
Sandford et al. | Transonic Study of Active Flutter Suppression Based on an Energy Concept |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |