发明名称 风洞常规测力试验中气动刚度与气动阻尼的视频测量方法
摘要 本发明公开了一种风洞常规测力试验中气动刚度与气动阻尼的视频测量方法,利用常规测力试验采用的迎角阶梯变化的风洞试验方式,通过迎角视频测量系统测得的迎角数据与天平测得的气动力数据,计算出该测力试验参数下,各迎角阶梯的气动阻尼和气动刚度,从而在风洞常规测力试验中,测得给定迎角阶梯气动力的同时,获得了该试验马赫数与雷诺数下的气动阻尼和气动刚度。本发明的积极效果是:实现了在常规测力试验中,采用常规测力试验模型,获得俯仰自由度的气动阻尼和气动刚度,扩展了常规测力试验获得的气动数据;同时,由于本发明采用常规测力试验模型,测得俯仰自由度的气动阻尼和气动刚度,较颤振试验模型或动导数试验模型研制成本低廉得多。
申请公布号 CN103234729A 申请公布日期 2013.08.07
申请号 CN201310119399.0 申请日期 2013.04.08
申请人 中国空气动力研究与发展中心高速空气动力研究所 发明人 路波;张征宇;王水亮;赵涛;朱龙
分类号 G01M9/06(2006.01)I 主分类号 G01M9/06(2006.01)I
代理机构 成都九鼎天元知识产权代理有限公司 51214 代理人 邓世燕
主权项 1.一种风洞常规测力试验中气动刚度与气动阻尼的视频测量方法,其特征在于:包括如下步骤:步骤一、测量各阶梯上风洞模型的迎角:利用安装在风洞试验段侧壁面上的迎角视频测量系统,通过测量风洞模型上的特征点三维坐标,计算出各阶梯风洞模型迎角;步骤二、计算风洞模型、测力天平与模型支撑机构构成的系统在俯仰自由度的一阶振型的频率ω<sub>0</sub>、弹性铰链常数k<sub>s</sub>、转动惯量I和结构阻尼c<sub>s</sub>:①在风洞未吹风时,对风洞模型施加初始位移或施加初始速度扰动后释放,然后测得模型迎角自由振动数据α(t);②对α(t)进行傅里叶变换,提取俯仰自由度的一阶振型频率ω<sub>0</sub>及阻尼比ξ<sub>0</sub>;③按如下公式计算振幅对数衰减率δ:<img file="FDA00003021894000011.GIF" wi="295" he="135" />式中Q<sub>0</sub>为迎角自由振动时初始振幅值,Q<sub>n</sub>为第n个周期后的迎角振幅值;④计算弹性铰链常数k<sub>s</sub>:对风洞模型施加已知的俯仰力矩M<sub>z</sub>,同时测出模型迎角变化量Δα,则<maths num="0001"><![CDATA[<math><mrow><msub><mi>k</mi><mi>s</mi></msub><mo>=</mo><mfrac><msub><mi>M</mi><mi>z</mi></msub><mi>&Delta;&alpha;</mi></mfrac></mrow></math>]]></maths>⑤计算俯仰自由度转动惯量I,计算式为:<img file="FDA00003021894000013.GIF" wi="154" he="134" />⑥计算不吹风时的阻尼系数q<sub>0</sub>和俯仰自由度的结构阻尼c<sub>s</sub>,计算式为:<maths num="0002"><![CDATA[<math><mrow><msub><mi>q</mi><mn>0</mn></msub><mo>=</mo><mo>-</mo><mfrac><msub><mi>&delta;&omega;</mi><mn>0</mn></msub><mrow><mn>2</mn><mi>&pi;</mi></mrow></mfrac></mrow></math>]]></maths>c<sub>s</sub>=2Iq<sub>0</sub>步骤三、计算各个阶梯上风洞模型的俯仰气动阻尼和气动刚度:①测量模型在各个阶梯的迎角数据;②求解各个阶梯的俯仰气动阻尼和气动刚度:设迎角视频测量系统测得模型在第i个阶梯的迎角数据为α<sup>i</sup>(t),测力天平在第i个阶梯测得的升力和俯仰力矩分别记为Y<sup>i</sup>和<img file="FDA00003021894000021.GIF" wi="104" he="79" />则A)计算俯仰力矩<img file="FDA00003021894000022.GIF" wi="77" he="68" />引起的迎角<img file="FDA00003021894000023.GIF" wi="171" he="89" />计算式为:<maths num="0003"><![CDATA[<math><mrow><msubsup><mi>&alpha;</mi><msub><mi>M</mi><mi>z</mi></msub><mi>i</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msup><mi>&alpha;</mi><mi>i</mi></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>k</mi><mi>y</mi></msub><mo>&CenterDot;</mo><msup><mi>Y</mi><mi>i</mi></msup></mrow></math>]]></maths>式中Yi为第i个阶梯的升力,k<sub>y</sub>为升力Y<sup>i</sup>对应的迎角弹性铰链常数,通过天平校准获得;B)计算第i个阶梯的气动刚度<img file="FDA00003021894000025.GIF" wi="80" he="81" /><maths num="0004"><![CDATA[<math><mrow><msubsup><mi>k</mi><mi>a</mi><mi>i</mi></msubsup><mo>=</mo><msubsup><mi>M</mi><mi>z</mi><mi>i</mi></msubsup><mo>/</mo><mover><mrow><msubsup><mi>&alpha;</mi><msub><mi>M</mi><mi>z</mi></msub><mi>i</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mo>&OverBar;</mo></mover><mo>-</mo><msub><mi>k</mi><mi>s</mi></msub></mrow></math>]]></maths>式中<img file="FDA00003021894000027.GIF" wi="150" he="110" />为<img file="FDA00003021894000028.GIF" wi="146" he="76" />在第i个阶梯的算术平均值;C)计算第i个阶梯的俯仰气动阻尼<img file="FDA00003021894000029.GIF" wi="52" he="72" />与俯仰运动的阻尼系数q<sup>i</sup>,计算式为<maths num="0005"><![CDATA[<math><mrow><msup><mi>q</mi><mi>i</mi></msup><mo>=</mo><msqrt><mfrac><mrow><msub><mi>k</mi><mi>s</mi></msub><mo>+</mo><msubsup><mi>k</mi><mi>a</mi><mi>i</mi></msubsup></mrow><mi>I</mi></mfrac><mo>-</mo><msup><mrow><mo>(</mo><msup><mi>&omega;</mi><mi>i</mi></msup><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><msubsup><mi>c</mi><mi>a</mi><mi>i</mi></msubsup><mo>=</mo><mn>2</mn><mi>I</mi><msup><mi>q</mi><mi>i</mi></msup><mo>-</mo><msub><mi>c</mi><mi>s</mi></msub></mrow></math>]]></maths>式中,ω<sup>i</sup>为第i个阶梯的迎角数据α<sup>i</sup>(t)的一阶振型频率。
地址 621000 四川省绵阳市涪城区剑门路西段278号副2号205研究室