发明名称 考虑摩擦力的自由振动动导数试验的数据分析方法
摘要 本发明涉及考虑摩擦力的自由振动动导数试验的数据分析方法,可从俯仰和偏航自由振动动稳定性风洞试验的试验数据中准确计算出飞行器的动稳定性导数,特别适用于系统摩擦力不可忽略的试验系统。本发明的技术方案是:首先在飞行器的振动动力学方程中添加摩擦力矩项,而后利用参数辨识技术从试验数据中辨识出相关气动参数,最后根据静态气动力与模型重力间关系建立筛选准则,获得最终的动稳定性导数数据。
申请公布号 CN105136423A 申请公布日期 2015.12.09
申请号 CN201510651344.3 申请日期 2015.10.10
申请人 中国航天空气动力技术研究院 发明人 张玉石;赵俊波;梁彬;付增良;高清
分类号 G01M9/02(2006.01)I;G06F19/00(2011.01)I 主分类号 G01M9/02(2006.01)I
代理机构 北京远大卓悦知识产权代理事务所(普通合伙) 11369 代理人 史霞
主权项 一种考虑摩擦力的自由振动动导数试验的数据分析方法,包含动力学分析、气动参数估计方法和气动参数筛选准则,其特征在于,包括:步骤一、通过自由振动动稳定性导数风洞试验,测量飞行器姿态角,其中姿态角为飞行器的俯仰角或偏航角;步骤二、以系统摩擦力矩和气动力矩为标准,当摩擦力矩与气动力矩同量级或大于气动力矩时,采用模型(1)~(4)得到俯仰稳定性气动导数或偏航稳定性导数;<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>I</mi><mover><mi>&theta;</mi><mo>&CenterDot;&CenterDot;</mo></mover><mo>=</mo><msub><mi>M</mi><mi>A</mi></msub><mo>+</mo><msub><mi>M</mi><mi>d</mi></msub><mo>+</mo><msub><mi>M</mi><mi>f</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000818427420000011.GIF" wi="734" he="92" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>M</mi><mi>A</mi></msub><mo>=</mo><msub><mi>QSc</mi><mi>A</mi></msub><mo>&lsqb;</mo><msub><mi>C</mi><mrow><mi>M</mi><mi>&theta;</mi></mrow></msub><mi>&Delta;</mi><mi>&theta;</mi><mo>+</mo><msub><mi>C</mi><mrow><mi>M</mi><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover></mrow></msub><mfrac><msub><mi>c</mi><mi>A</mi></msub><msub><mi>V</mi><mi>&infin;</mi></msub></mfrac><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mo>&rsqb;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000818427420000012.GIF" wi="944" he="161" /></maths>M<sub>d</sub>=‑K<sub>d</sub>Δθ           (3)<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>M</mi><mi>f</mi></msub><mo>=</mo><mo>-</mo><msub><mi>K</mi><mi>f</mi></msub><mi>s</mi><mi>g</mi><mi>n</mi><mrow><mo>(</mo><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000818427420000015.GIF" wi="951" he="104" /></maths>当姿态角为俯仰角时,θ为飞行器的通过自由振动动稳定性导数风洞试验测得的俯仰角,S为飞行器模型的俯仰方向参考面积,I为飞行器俯仰转动惯量,M<sub>A</sub>为飞行器的俯仰气动力矩,M<sub>d</sub>为自由振动动导数天平弹性铰链的恢复力矩,M<sub>f</sub>为自由振动动导数系统的摩擦力矩,Q为风洞气流的动压,c<sub>A</sub>为飞行器俯仰方向的参考长度,C<sub>Mθ</sub>为俯仰静稳定性气动导数,Δθ为飞行器相对于平衡点的俯仰角变化量,[QSc<sub>A</sub>C<sub>Mθ</sub>Δθ]为气动恢复力矩,<img file="FDA0000818427420000013.GIF" wi="317" he="145" />为气动阻尼力矩,V<sub>∞</sub>为风洞气流速度,K<sub>d</sub>为天平铰链的恢复力矩系数,K<sub>f</sub>为系统摩擦力矩系数,<img file="FDA0000818427420000014.GIF" wi="100" he="78" />为俯仰动稳定性导数;当姿态角为偏航角时,θ为飞行器的通过自由振动动稳定性导数风洞试验测得的偏航角,S为飞行器模型的偏航方向参考面积,I为飞行器偏航转动惯量,M<sub>A</sub>为飞行器的偏航气动力矩,M<sub>d</sub>为自由振动动导数天平弹性铰链的恢复力矩,M<sub>f</sub>为自由振动动导数系统的摩擦力矩,Q为风洞气流的动压,c<sub>A</sub>为飞行器偏航方向的参考长度,C<sub>Mθ</sub>为偏航静稳定性气动导数,Δθ为飞行器相对于平衡点的偏航角变化量,[QSc<sub>A</sub>C<sub>Mθ</sub>Δθ]为气动恢复力矩,<img file="FDA0000818427420000021.GIF" wi="310" he="154" />为气动阻尼力矩,V<sub>∞</sub>为风洞气流速度,K<sub>d</sub>为天平铰链的恢复力矩系数,K<sub>f</sub>为系统摩擦力矩系数,<img file="FDA0000818427420000022.GIF" wi="99" he="77" />为偏航动稳定性导数。
地址 100074 北京市丰台区云岗西路17号