CN108120581A - 旋转导弹俯仰动导数高速风洞试验装置及方法 - Google Patents
旋转导弹俯仰动导数高速风洞试验装置及方法 Download PDFInfo
- Publication number
- CN108120581A CN108120581A CN201711307469.XA CN201711307469A CN108120581A CN 108120581 A CN108120581 A CN 108120581A CN 201711307469 A CN201711307469 A CN 201711307469A CN 108120581 A CN108120581 A CN 108120581A
- Authority
- CN
- China
- Prior art keywords
- angle
- attack
- moment
- guided missile
- pitching
- 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
Landscapes
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
一种旋转导弹俯仰动导数高速风洞试验装置及方法。该装置包括:高速风洞;支撑机构,支撑机构用于将导弹支撑于高速风洞内,且能够驱动导弹旋转和进行强迫俯仰振动;动导数天平,动导数天平设置于导弹的内部,用于测量导弹的力矩信号;位移元件,位移元件设置于支撑机构上,用于测量导弹的振动角位移信号;处理器,处理器针对攻角序列中的每一个攻角,计算俯仰动导数。在高速风洞中在导弹绕自身轴线旋转的同时,依靠支撑机构驱动导弹做强迫俯仰振动,能够精确简便地计算在高速风洞中导弹旋转时的俯仰动导数。
Description
技术领域
本发明涉及试验空气动力学领域,更具体地,涉及一种旋转导弹俯仰动导数高速风洞试验装置及方法。
背景技术
旋转飞行可以简化导弹的控制***,用一个控制通道实现俯仰和偏航两个方向的控制;可以减小推力偏心、质量偏心、气动偏心等非对称因素对飞行性能的不利影响。但是,旋转飞行也带来了一系列复杂的空气动力学问题,飞行器绕其轴线旋转飞行时,旋转对边界层的剪切效应引起体涡的非对称分离以及引起的边界层转捩区的非对称,从而进一步引起气动特性的变化,尤其是动导数特性的变化。
俯仰动导数是战术飞行武器的重要动态参数,是动态品质分析和控制***设计不可缺少的数据。通过高速风洞试验获取动导数参数是主要研究方法之一。高速风洞试验一般利用振动法测定动导数,振动法通常可分为两种类型:一种是自由振动法,另一种是强迫振动法。实践证明强迫振动法计算工作量小而且准确。
通过风洞试验研究旋转情况下的俯仰动导数参数时,早期的研究人员寄希望于将稳定旋转的气动数据和单自由度动导数试验数据结合起来,估算旋转流场下的振荡动导数,但是这种方法并不可靠。Jaccob Kay曾在12英尺立式风洞中进行了旋转天平加强迫振荡试验,试验数据表明在大攻角范围内,振荡旋转的非定常气动特性与稳定旋转相比有较大差距,体现出更强的非线性。因此,国内外研究机构相继开发各自的试验技术来进行相关的研究,但主要在低速风洞进行研究。旋转导弹的飞行速度一般在高速范围,高速情况下旋转导弹的俯仰动导数特性目前主要通过数值计算进行研究,受到高速风洞载荷、尺寸、堵塞度等因素限制,在高速风洞中进行地面模拟比较罕见。因此,有必要开发一种旋转导弹俯仰动导数高速风洞试验装置及方法。
发明内容
本发明提出了一种旋转导弹俯仰动导数高速风洞试验装置及方法,在高速风洞中在导弹绕自身轴线旋转的同时,依靠支撑机构驱动导弹做强迫俯仰振动,能够精确简便地计算在高速风洞中导弹旋转时的俯仰动导数。
根据本发明的一方面,提出了一种旋转导弹俯仰动导数高速风洞试验装置,包括:高速风洞;支撑机构,所述支撑机构用于将导弹支撑于所述高速风洞内,且能够驱动所述导弹旋转和进行强迫俯仰振动;动导数天平,所述动导数天平设置于所述导弹的内部,用于测量所述导弹的力矩信号;位移元件,所述位移元件设置于所述支撑机构上,用于测量所述导弹的振动角位移信号;处理器,所述处理器针对攻角序列中的每一个攻角,根据以下公式(1)计算所述俯仰动导数:
其中,C为俯仰动导数,为俯仰时差力矩系数,Cmq为俯仰阻尼力矩系数,为气动力矩基波幅值,λ为振动角位移信号和总力矩信号之间的相位差,θ0为角位移的波动幅值,q∞为动压,S为导弹的参考面积,cA为导弹的参考长度,K为减缩频率。
优选地,所述支撑机构为扇片形,其堵塞度小于10%。
优选地,所述处理器通过以下步骤确定气动力矩基波幅值:用总力矩信号减去惯性力矩得到气动力矩变化曲线,其中,所述总力矩信号为在有风情况下,通过所述动导数天平测量获得,所述惯性力矩基于在无风情况下,通过所述动导数天平测量获得的惯性力矩信号得到;根据所述气动力矩变化曲线确定所述气动力矩基波幅值。
优选地,所述处理器通过以下步骤确定每一个攻角对应的惯性力矩:在所述攻角下,以指定的采样间隔测量一个强迫俯仰振动周期内的惯性力矩信号,以所述周期内测量的惯性力矩信号的平均值作为所述攻角对应的惯性力矩。
优选地,所述处理器通过以下公式(2)计算所述减缩频率:
K=ωcA/(2V) (2)
其中,ω为强迫俯仰振动频率,V为风速。
优选地,所述旋转为绕导弹自身轴线的自旋转,所述强迫俯仰振动为绕导弹质心的正弦强迫俯仰振动。
根据本发明的另一方面,提出了一种旋转导弹俯仰动导数高速风洞试验方法,可以包括:在无风情况下,针对攻角序列中的每一个攻角,使导弹在所述攻角下进行旋转和强迫俯仰振动,确定每一个攻角对应的惯性力矩;在有风情况下,针对所述攻角序列中的每一个攻角,使导弹在所述攻角下进行旋转和强迫俯仰振动,测量每一个攻角对应的总力矩信号,并采集每一个攻角对应的振动角位移信号,其中,在有风和无风情况下,所述导弹的旋转速度和强迫俯仰振动的振动角速度和振幅相同;针对攻角序列中的每一个攻角,根据以下公式(1)计算所述俯仰动导数:
其中,C为俯仰动导数,为俯仰时差力矩系数,Cmq为俯仰阻尼力矩系数,为气动力矩基波幅值,λ为振动角位移信号和总力矩信号之间的相位差,θ0为角位移的波动幅值,q∞为动压,S为导弹的参考面积,cA为导弹的参考长度,K为减缩频率。
优选地,通过以下步骤确定气动力矩基波幅值:用所述总力矩信号减去所述惯性力矩得到气动力矩变化曲线;根据所述气动力矩变化曲线确定所述气动力矩基波幅值。
优选地,通过以下步骤确定每一个攻角对应的惯性力矩:在所述攻角下,以指定的采样间隔测量一个强迫俯仰振动周期内的惯性力矩信号,以所述周期内测量的惯性力矩信号的平均值作为所述攻角对应的惯性力矩。
优选地,通过以下公式(2)计算所述减缩频率:
K=ωcA/(2V) (2)
其中,ω为强迫俯仰振动频率,V为风速。
优选地,所述旋转为绕导弹自身轴线的自旋转,所述强迫俯仰振动为绕导弹质心的正弦强迫俯仰振动。
本发明的有益效果在于:在高速风洞中在导弹绕自身轴线旋转的同时,依靠支撑机构驱动旋转导弹做强迫俯仰振动,解决了在高速风洞中导弹旋转时的俯仰动导数测量问题。
本发明的方法和装置具有其它的特性和优点,这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的,或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述,这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述,本发明的上述以及其它目的、特征和优势将变得更加明显,其中,在本发明示例性实施例中,相同的参考标号通常代表相同部件。
图1示出了根据本发明的一个实施例的旋转导弹俯仰动导数高速风洞试验装置的示意图。
图2示出了根据本发明的旋转导弹俯仰动导数高速风洞试验方法的步骤的流程图。
附图标记说明:
1、高速风洞;2、导弹;3、支撑机构;4、动导数天平;5、位移元件。
具体实施方式
下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例,然而应该理解,可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了使本发明更加透彻和完整,并且能够将本发明的范围完整地传达给本领域的技术人员。
图1示出了根据本发明的一个实施例的旋转导弹俯仰动导数高速风洞试验装置的示意图。
旋转导弹俯仰动导数高速风洞试验装置包括:高速风洞1,其马赫数范围为0.3-4.0;支撑机构3,支撑机构3用于将导弹2支撑于高速风洞1内,且能够驱动导弹2旋转和进行强迫俯仰振动;动导数天平4,动导数天平4设置于导弹2的内部,用于测量导弹2的力矩信号;位移元件5,位移元件5设置于支撑机构3上,用于测量导弹2的振动角位移信号;处理器,处理器针对攻角序列中的每一个攻角,根据以下公式(1)计算俯仰动导数:
其中,C为俯仰动导数,为俯仰时差力矩系数,Cmq为俯仰阻尼力矩系数,为气动力矩基波幅值,λ为振动角位移信号和总力矩信号之间的相位差,θ0为角位移的波动幅值,q∞为动压,S为导弹2的参考面积,cA为导弹2的参考长度,K为减缩频率。
具体地,公式(1)中的俯仰动导数为俯仰时差力矩系数与俯仰阻尼力矩系数之和,在实际应用中,将俯仰时差力矩系数与俯仰阻尼力矩系数组合使用,作为俯仰动导数,因此,在本发明中根据公式(1)计算俯仰动导数。
在一个示例中,支撑机构3为扇片形,其堵塞度小于10%,即支撑机构3的迎风面积与高速风洞1的横截面积的比值小于10%,可承受高速风洞1载荷,满足高速风洞1尺寸要求。
在一个示例中,处理器通过以下步骤确定气动力矩基波幅值:用总力矩信号减去惯性力矩得到气动力矩变化曲线,其中,总力矩信号为在有风情况下,通过动导数天平4测量获得,惯性力矩基于在无风情况下,通过动导数天平4测量获得的惯性力矩信号得到;根据气动力矩变化曲线确定气动力矩基波幅值。
在一个示例中,处理器通过以下步骤确定每一个攻角对应的惯性力矩:在攻角下,以指定的采样间隔测量一个强迫俯仰振动周期内的惯性力矩信号,以周期内测量的惯性力矩信号的平均值作为攻角对应的惯性力矩。
在一个示例中,处理器通过以下公式(2)计算减缩频率:
K=ωcA/(2V) (2)
其中,ω为强迫俯仰振动频率,V为风速。
在一个示例中,旋转为绕导弹2自身轴线的自旋转,强迫俯仰振动为绕导弹质心的正弦强迫俯仰振动。
在一个示例中,处理器通过以下公式(3)计算动压:
其中,q∞为动压,M∞为马赫数,p∞为静压,静压由高速风洞测控***测量得到。
具体地,高速风洞1能够提供吹风条件和控制马赫数,通过马赫数可以确定风速,其马赫数范围为0.3-4.0;支撑机构3用于将导弹2支撑于高速风洞1内,且能够驱动导弹2旋转和进行强迫俯仰振动,其中,旋转为绕导弹2自身轴线的自旋转,强迫俯仰振动为正弦强迫俯仰振动,特别地,其是绕导弹质心的正弦强迫俯仰振动;处理器可以为计算机。
图2示出了根据本发明的旋转导弹俯仰动导数高速风洞试验方法的步骤的流程图。
在该实施例中,根据本发明的旋转导弹俯仰动导数高速风洞试验方法可以包括:
步骤101,在无风情况下,针对攻角序列中的每一个攻角,使导弹在该攻角下进行旋转和强迫俯仰振动,确定每一个攻角对应的惯性力矩;
具体地,在无风情况下,首先通过支撑机构将导弹调整至指定的攻角,然后通过支撑机构驱动导弹旋转,并驱动导弹作强迫俯仰振动,通过动导数天平测量相应的惯性力矩信号。然后,将导弹调整至下一攻角,重复前述过程,直到针对攻角序列中的每一个攻角均测量完毕。
步骤102,在有风情况下,针对攻角序列中的每一个攻角,使导弹在攻角下进行旋转和强迫俯仰振动,测量每一个攻角对应的总力矩信号,并采集每一个攻角对应的振动角位移信号,其中,在有风和无风情况下,导弹的旋转速度和强迫俯仰振动的振动角速度和振幅相同;
具体地,启动高速风洞,高速风洞提供吹风条件,并可控制马赫数,通过马赫数可确定高速风洞内的风速,其马赫数范围为0.3-4.0。待高速风洞内的流场稳定到指定的马赫数后,通过支撑机构将导弹调整至指定的攻角,通过支撑机构驱动导弹旋转,并驱动导弹作强迫俯仰振动,通过动导数天平测量相应的总力矩信号,通过位移元件测量导弹的振动角位移信号。然后,将导弹调整至下一攻角,重复前述过程,直到针对攻角序列中的每一个攻角均测量完毕。
步骤103,针对攻角序列中的每一个攻角,根据以下公式(1)计算俯仰动导数:
其中,C为俯仰动导数,为俯仰时差力矩系数,Cmq为俯仰阻尼力矩系数,为气动力矩基波幅值,λ为振动角位移信号和总力矩信号之间的相位差,θ0为角位移的波动幅值,q∞为动压,S为导弹的参考面积,cA为导弹的参考长度,K为减缩频率。
在一个示例中,通过以下步骤确定气动力矩基波幅值:用总力矩信号减去惯性力矩得到气动力矩变化曲线;根据气动力矩变化曲线确定气动力矩基波幅值。
在一个示例中,通过以下步骤确定每一个攻角对应的惯性力矩:在攻角下,以指定的采样间隔测量一个强迫俯仰振动周期内的惯性力矩信号,以周期内测量的惯性力矩信号的平均值作为攻角对应的惯性力矩。
在一个示例中,通过以下公式(2)计算减缩频率:
K=ωcA/(2V) (2)
其中,ω为强迫俯仰振动频率,V为风速。
在一个示例中,旋转为绕导弹自身轴线的自旋转,强迫俯仰振动为绕导弹质心的正弦强迫俯仰振动。
在一个示例中,还可以先对惯性力矩信号、总力矩信号、振动角位移信号进行滤波,再进行相应的计算。
本发明在高速风洞中在旋转导弹绕自身轴线旋转的同时,依靠支撑机构驱动旋转导弹做强迫俯仰振动,解决了在高速风洞中旋转导弹旋转时的俯仰动导数测量问题。
为便于理解本发明实施例的方案及其效果,以下给出三个具体应用示例。本领域技术人员应理解,该示例仅为了便于理解本发明,其任何具体细节并非意在以任何方式限制本发明。
应用示例1
旋转导弹俯仰动导数高速风洞试验装置包括:高速风洞1,提供吹风条件和控制马赫数,通过马赫数确定风速;支撑机构3为扇片形,用于将导弹2支撑于高速风洞1内,且能够驱动导弹2旋转和进行强迫俯仰振动,其中,旋转为绕导弹2自身轴线的自旋转,强迫俯仰振动为绕导弹质心的正弦强迫俯仰振动;动导数天平4,设置于导弹2的内部,用于测量导弹2的力矩信号;位移元件5,设置于支撑机构3上,用于测量导弹2的振动角位移信号;处理器,为计算机,针对攻角序列中的每一个攻角,在无风情况下,通过动导数天平4以指定的采样间隔测量一个强迫俯仰振动周期内的惯性力矩信号,以周期内测量的惯性力矩信号的平均值作为攻角对应的惯性力矩,为在有风情况下,通过动导数天平4测量获得总力矩信号,用总力矩信号减去惯性力矩得到气动力矩变化曲线,根据气动力矩变化曲线确定气动力矩基波幅值,动压通过风速计算获得,导弹2的参考面积与参考长度为导弹2的基本信息,进而根据公式(1)计算俯仰动导数,其中,通过公式(2)计算减缩频率。
应用示例2
旋转导弹俯仰动导数高速风洞试验方法包括以下步骤:
步骤101:在无风情况下,针对攻角序列中的0°,使导弹在0°攻角下进行旋转和强迫俯仰振动,其中,旋转为绕导弹自身轴线的自旋转,转速为300r/min,强迫俯仰振动为绕导弹质心的正弦强迫俯仰振动,振动方程为θ=3°sin(2πt)。在0°攻角下,以指定的采样间隔测量一个强迫俯仰振动周期内的惯性力矩信号,以周期内测量的惯性力矩信号的平均值作为攻角对应的惯性力矩;
步骤102:高速风洞开车,提供吹风条件并控制马赫数为0.6,通过马赫数确定风速,在有风情况下,针对攻角序列中的0°,使导弹在0°攻角下进行旋转和强迫俯仰振动,测量0°攻角对应的总力矩信号,并采集0°攻角对应的振动角位移信号,其中,在有风和无风情况下,导弹的旋转速度均为300r/min,强迫俯仰振动方程均为θ=3°sin(2πt);
步骤103:针对攻角序列中的0°,用总力矩信号减去惯性力矩得到气动力矩变化曲线;根据气动力矩变化曲线确定气动力矩基波幅值,动压通过风速计算获得,导弹的参考面积与参考长度为导弹的基本信息,通过公式(2)计算减缩频率,进而根据公式(1)计算俯仰动导数为-2.03rad-1。
应用示例3
旋转导弹俯仰动导数高速风洞试验方法包括以下步骤:
步骤101:在无风情况下,针对攻角序列中的0°,使导弹在0°攻角下进行旋转和强迫俯仰振动,其中,旋转为绕导弹自身轴线的自旋转,转速为600r/min,强迫俯仰振动为绕导弹质心的正弦强迫俯仰振动,振动方程为θ=3°sin(4πt)。在0°攻角下,以指定的采样间隔测量一个强迫俯仰振动周期内的惯性力矩信号,以周期内测量的惯性力矩信号的平均值作为攻角对应的惯性力矩;
步骤102:高速风洞开车,提供吹风条件并控制马赫数为0.8,通过马赫数确定风速,在有风情况下,针对攻角序列中的0°,使导弹在0°攻角下进行旋转和强迫俯仰振动,测量0°攻角对应的总力矩信号,并采集0°攻角对应的振动角位移信号,其中,在有风和无风情况下,导弹的旋转速度均为600r/min,强迫俯仰振动方程均为θ=3°sin(4πt);
步骤103:针对攻角序列中的0°,用总力矩信号减去惯性力矩得到气动力矩变化曲线;根据气动力矩变化曲线确定气动力矩基波幅值,动压通过风速计算获得,导弹的参考面积与参考长度为导弹的基本信息,通过公式(2)计算减缩频率,进而根据公式(1)计算俯仰动导数为-2.48rad-1。
综上所述,本发明在高速风洞中在旋转导弹绕自身轴线旋转的同时,依靠支撑机构驱动旋转导弹做强迫俯仰振动,解决了在高速风洞中旋转导弹旋转时的俯仰动导数测量问题。
以上已经描述了本发明的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。
Claims (11)
1.一种旋转导弹俯仰动导数高速风洞试验装置,其特征在于,该装置包括:
高速风洞;
支撑机构,所述支撑机构用于将导弹支撑于所述高速风洞内,且能够驱动所述导弹旋转和进行强迫俯仰振动;
动导数天平,所述动导数天平设置于所述导弹的内部,用于测量所述导弹的力矩信号;
位移元件,所述位移元件设置于所述支撑机构上,用于测量所述导弹的振动角位移信号;
处理器,所述处理器针对攻角序列中的每一个攻角,根据以下公式(1)计算所述俯仰动导数:
<mrow>
<mi>C</mi>
<mo>=</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>m</mi>
<mover>
<mi>&alpha;</mi>
<mo>&CenterDot;</mo>
</mover>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>m</mi>
<mi>q</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mover>
<mi>M</mi>
<mo>&OverBar;</mo>
</mover>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&lambda;</mi>
</mrow>
<mrow>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
<msub>
<mi>q</mi>
<mi>&infin;</mi>
</msub>
<msub>
<mi>Sc</mi>
<mi>A</mi>
</msub>
<mi>K</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,C为俯仰动导数,为俯仰时差力矩系数,Cmq为俯仰阻尼力矩系数,为气动力矩基波幅值,λ为振动角位移信号和总力矩信号之间的相位差,θ0为角位移的波动幅值,q∞为动压,S为导弹的参考面积,cA为导弹的参考长度,K为减缩频率。
2.根据权利要求1所述的旋转导弹俯仰动导数高速风洞试验装置,其中,所述支撑机构为扇片形,其堵塞度小于10%。
3.根据权利要求1所述的旋转导弹俯仰动导数高速风洞试验装置,其中,所述处理器通过以下步骤确定气动力矩基波幅值:
用总力矩信号减去惯性力矩得到气动力矩变化曲线,其中,所述总力矩信号为在有风情况下,通过所述动导数天平测量获得,所述惯性力矩基在无风情况下通过所述动导数天平测量获得的惯性力矩信号得到;
根据所述气动力矩变化曲线确定所述气动力矩基波幅值。
4.根据权利要求3所述的旋转导弹俯仰动导数高速风洞试验装置,其中,所述处理器通过以下步骤确定每一个攻角对应的惯性力矩:
在所述攻角下,以指定的采样间隔测量一个强迫俯仰振动周期内的惯性力矩信号,以所述周期内测量的惯性力矩信号的平均值作为所述攻角对应的惯性力矩。
5.根据权利要求1所述的旋转导弹俯仰动导数高速风洞试验装置,其中,所述处理器通过以下公式(2)计算所述减缩频率:
K=ωcA/(2V) (2)
其中,ω为强迫俯仰振动频率,V为风速。
6.根据权利要求1所述的旋转导弹俯仰动导数高速风洞试验装置,其中,所述旋转为绕导弹自身轴线的自旋转,所述强迫俯仰振动为绕导弹质心的正弦强迫俯仰振动。
7.一种旋转导弹俯仰动导数高速风洞试验方法,利用根据权利要求1所述的旋转导弹俯仰动导数高速风洞试验装置,所述方法包括:
在无风情况下,针对攻角序列中的每一个攻角,使导弹在所述攻角下进行旋转和强迫俯仰振动,确定每一个攻角对应的惯性力矩;
在有风情况下,针对所述攻角序列中的每一个攻角,使导弹在所述攻角下进行旋转和强迫俯仰振动,测量每一个攻角对应的总力矩信号,并采集每一个攻角对应的振动角位移信号,其中,在有风和无风情况下,所述导弹的旋转速度和强迫俯仰振动的振动角速度和振幅相同;
针对攻角序列中的每一个攻角,根据以下公式(1)计算所述俯仰动导数:
<mrow>
<mi>C</mi>
<mo>=</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>m</mi>
<mover>
<mi>&alpha;</mi>
<mo>&CenterDot;</mo>
</mover>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>m</mi>
<mi>q</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mover>
<mi>M</mi>
<mo>&OverBar;</mo>
</mover>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&lambda;</mi>
</mrow>
<mrow>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
<msub>
<mi>q</mi>
<mi>&infin;</mi>
</msub>
<msub>
<mi>Sc</mi>
<mi>A</mi>
</msub>
<mi>K</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,C为俯仰动导数,为俯仰时差力矩系数,Cmq为俯仰阻尼力矩系数,为气动力矩基波幅值,λ为振动角位移信号和总力矩信号之间的相位差,θ0为角位移的波动幅值,q∞为动压,S为导弹的参考面积,cA为导弹的参考长度,K为减缩频率。
8.根据权利要求7所述的旋转导弹俯仰动导数高速风洞试验方法,其中,通过以下步骤确定气动力矩基波幅值:
用所述总力矩信号减去所述惯性力矩得到气动力矩变化曲线;
根据所述气动力矩变化曲线确定所述气动力矩基波幅值。
9.根据权利要求7所述的旋转导弹俯仰动导数高速风洞试验方法,其中,通过以下步骤确定每一个攻角对应的惯性力矩:
在所述攻角下,以指定的采样间隔测量一个强迫俯仰振动周期内的惯性力矩信号,以所述周期内测量的惯性力矩信号的平均值作为所述攻角对应的惯性力矩。
10.根据权利要求7所述的旋转导弹俯仰动导数高速风洞试验方法,其中,通过以下公式(2)计算所述减缩频率:
K=ωcA/(2V) (2)
其中,ω为强迫俯仰振动频率,V为风速。
11.根据权利要求7所述的旋转导弹俯仰动导数高速风洞试验方法,其中,所述旋转为绕导弹自身轴线的自旋转,所述强迫俯仰振动为绕导弹质心的正弦强迫俯仰振动。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711307469.XA CN108120581B (zh) | 2017-12-11 | 2017-12-11 | 旋转导弹俯仰动导数高速风洞试验装置及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711307469.XA CN108120581B (zh) | 2017-12-11 | 2017-12-11 | 旋转导弹俯仰动导数高速风洞试验装置及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108120581A true CN108120581A (zh) | 2018-06-05 |
CN108120581B CN108120581B (zh) | 2020-07-28 |
Family
ID=62229726
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711307469.XA Active CN108120581B (zh) | 2017-12-11 | 2017-12-11 | 旋转导弹俯仰动导数高速风洞试验装置及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108120581B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109557901A (zh) * | 2018-11-21 | 2019-04-02 | 中国航空工业集团公司哈尔滨空气动力研究所 | 一种用于风洞旋转流场下模型振荡试验的控制*** |
CN109612680A (zh) * | 2019-01-24 | 2019-04-12 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种可复校核的双位移滚转动导数试验装置 |
CN109632252A (zh) * | 2018-12-27 | 2019-04-16 | 中国航天空气动力技术研究院 | 外式强迫振动动导数试验的振动角位移测量装置及方法 |
CN110044570A (zh) * | 2019-04-16 | 2019-07-23 | 沈阳航空航天大学 | 一种旋成体机身测压实验气动震颤误差修正方法 |
CN110940483A (zh) * | 2019-11-13 | 2020-03-31 | 中国航天空气动力技术研究院 | 用于大长细比飞行器俯仰偏航自由振动动导数试验装置 |
CN110989338A (zh) * | 2019-12-10 | 2020-04-10 | 北京理工大学 | 考虑气动非线性的飞行器旋转稳定控制***及方法 |
CN112067245A (zh) * | 2020-09-07 | 2020-12-11 | 中国航空工业集团公司沈阳空气动力研究所 | 一种高速风洞平移振动动导数试验装置及试验方法 |
CN116952524A (zh) * | 2023-09-18 | 2023-10-27 | 中国航空工业集团公司沈阳空气动力研究所 | 一种高速风洞动导数试验监测方法、电子设备及存储介质 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100132446A1 (en) * | 2008-12-03 | 2010-06-03 | David A Corder | Wind tunnel testing technique |
CN101726401A (zh) * | 2009-12-09 | 2010-06-09 | 中国航空工业第一集团公司沈阳空气动力研究所 | 一种用于俯仰动导数实验的天平测量装置 |
CN102494864A (zh) * | 2011-11-24 | 2012-06-13 | 北京航空航天大学 | 飞行器俯仰运动下偏航/滚转***模拟装置 |
CN104850759A (zh) * | 2015-06-16 | 2015-08-19 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种风洞强迫振动动稳定性导数试验数据处理方法 |
CN105136423A (zh) * | 2015-10-10 | 2015-12-09 | 中国航天空气动力技术研究院 | 考虑摩擦力的自由振动动导数试验的数据分析方法 |
CN105258904A (zh) * | 2015-11-02 | 2016-01-20 | 中国航天空气动力技术研究院 | 一种旋转导弹锥形运动稳定性风洞试验方法 |
CN105527068A (zh) * | 2015-12-29 | 2016-04-27 | 中国航天空气动力技术研究院 | 一种鸭式布局旋转导弹的风洞动态测力试验装置 |
CN205642791U (zh) * | 2015-12-29 | 2016-10-12 | 中国航天空气动力技术研究院 | 一种风洞用旋转导弹的锥形运动模拟装置 |
CN106768816A (zh) * | 2016-12-22 | 2017-05-31 | 中国航空工业集团公司沈阳空气动力研究所 | 一种尾振动的俯仰动导数实验测量装置 |
CN106828819A (zh) * | 2017-01-23 | 2017-06-13 | 哈尔滨工程大学 | 船舶减摇鳍与翼鳍矢量控制方法 |
CN106840574A (zh) * | 2016-12-21 | 2017-06-13 | 中国航天空气动力技术研究院 | 一种用于风洞动导数强迫振动试验的装置 |
-
2017
- 2017-12-11 CN CN201711307469.XA patent/CN108120581B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100132446A1 (en) * | 2008-12-03 | 2010-06-03 | David A Corder | Wind tunnel testing technique |
CN101726401A (zh) * | 2009-12-09 | 2010-06-09 | 中国航空工业第一集团公司沈阳空气动力研究所 | 一种用于俯仰动导数实验的天平测量装置 |
CN102494864A (zh) * | 2011-11-24 | 2012-06-13 | 北京航空航天大学 | 飞行器俯仰运动下偏航/滚转***模拟装置 |
CN104850759A (zh) * | 2015-06-16 | 2015-08-19 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种风洞强迫振动动稳定性导数试验数据处理方法 |
CN105136423A (zh) * | 2015-10-10 | 2015-12-09 | 中国航天空气动力技术研究院 | 考虑摩擦力的自由振动动导数试验的数据分析方法 |
CN105258904A (zh) * | 2015-11-02 | 2016-01-20 | 中国航天空气动力技术研究院 | 一种旋转导弹锥形运动稳定性风洞试验方法 |
CN105527068A (zh) * | 2015-12-29 | 2016-04-27 | 中国航天空气动力技术研究院 | 一种鸭式布局旋转导弹的风洞动态测力试验装置 |
CN205642791U (zh) * | 2015-12-29 | 2016-10-12 | 中国航天空气动力技术研究院 | 一种风洞用旋转导弹的锥形运动模拟装置 |
CN106840574A (zh) * | 2016-12-21 | 2017-06-13 | 中国航天空气动力技术研究院 | 一种用于风洞动导数强迫振动试验的装置 |
CN106768816A (zh) * | 2016-12-22 | 2017-05-31 | 中国航空工业集团公司沈阳空气动力研究所 | 一种尾振动的俯仰动导数实验测量装置 |
CN106828819A (zh) * | 2017-01-23 | 2017-06-13 | 哈尔滨工程大学 | 船舶减摇鳍与翼鳍矢量控制方法 |
Non-Patent Citations (10)
Title |
---|
中国航空工业空气动力研究院: "《航空气动力技术》", 31 December 2013, 航空工业出版社 * |
卢学成 等: "超音速、高超音速飞行器动导数的高效计算方法", 《航空计算技术》 * |
张瑞民 等: "稳定性导数计算方法研究", 《航空工程进展》 * |
李周复: "《风洞试验手册》", 31 January 2015, 航空工业出版社 * |
王明慧等: "《渝黔铁路白沙沱长江大桥建造关键技术》", 31 October 2015, 西南交通出版社 * |
许可法 等: "旋转导弹风洞六自由度自由飞动导数实验研究", 《空气动力学学报》 * |
许晓斌: "《常规高超声速风洞与试验技术》", 31 August 2015, 国防工业出版社 * |
贾区耀: "天空飞行与地面风洞实验动态气动相关性研究", 《实验流体力学》 * |
路波: "《高速风洞测力试验数据处理方法》", 31 July 2014, 国防工业出版社 * |
陈怀瑾: "《防空导弹武器***总体设计和试验》", 31 December 1995, 中国宇航出版社 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109557901A (zh) * | 2018-11-21 | 2019-04-02 | 中国航空工业集团公司哈尔滨空气动力研究所 | 一种用于风洞旋转流场下模型振荡试验的控制*** |
CN109557901B (zh) * | 2018-11-21 | 2021-06-18 | 中国航空工业集团公司哈尔滨空气动力研究所 | 一种用于风洞旋转流场下模型振荡试验的控制*** |
CN109632252A (zh) * | 2018-12-27 | 2019-04-16 | 中国航天空气动力技术研究院 | 外式强迫振动动导数试验的振动角位移测量装置及方法 |
CN109612680A (zh) * | 2019-01-24 | 2019-04-12 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种可复校核的双位移滚转动导数试验装置 |
CN109612680B (zh) * | 2019-01-24 | 2024-01-30 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种可复校核的双位移滚转动导数试验装置 |
CN110044570B (zh) * | 2019-04-16 | 2021-03-16 | 沈阳航空航天大学 | 一种旋成体机身测压实验气动震颤误差修正方法 |
CN110044570A (zh) * | 2019-04-16 | 2019-07-23 | 沈阳航空航天大学 | 一种旋成体机身测压实验气动震颤误差修正方法 |
CN110940483B (zh) * | 2019-11-13 | 2021-12-07 | 中国航天空气动力技术研究院 | 用于大长细比飞行器俯仰偏航自由振动动导数试验装置 |
CN110940483A (zh) * | 2019-11-13 | 2020-03-31 | 中国航天空气动力技术研究院 | 用于大长细比飞行器俯仰偏航自由振动动导数试验装置 |
CN110989338A (zh) * | 2019-12-10 | 2020-04-10 | 北京理工大学 | 考虑气动非线性的飞行器旋转稳定控制***及方法 |
CN112067245A (zh) * | 2020-09-07 | 2020-12-11 | 中国航空工业集团公司沈阳空气动力研究所 | 一种高速风洞平移振动动导数试验装置及试验方法 |
CN116952524A (zh) * | 2023-09-18 | 2023-10-27 | 中国航空工业集团公司沈阳空气动力研究所 | 一种高速风洞动导数试验监测方法、电子设备及存储介质 |
CN116952524B (zh) * | 2023-09-18 | 2023-12-05 | 中国航空工业集团公司沈阳空气动力研究所 | 一种高速风洞动导数试验监测方法、电子设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN108120581B (zh) | 2020-07-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108120581A (zh) | 旋转导弹俯仰动导数高速风洞试验装置及方法 | |
Ozbek et al. | Challenges in testing and monitoring the in-operation vibration characteristics of wind turbines | |
CN103592467B (zh) | 二维超声波风速仪零点在线自校正装置及方法 | |
CN103003565A (zh) | 用于确定风力涡轮机的转子叶片的弯曲角的方法和装置 | |
CN206959778U (zh) | 一种基于滑坡***移和坡度姿态变化的监测*** | |
CN103429852A (zh) | 通过压力监控确定风扇参数 | |
CN105136423B (zh) | 考虑摩擦力的自由振动动导数试验的数据分析方法 | |
CN102636183A (zh) | 基于光纤监测和双轴转台离心机的挠性陀螺二次过载项测试方法 | |
CN102410139A (zh) | 确定转子共同叶片频率的方法和装置及操作风轮机的方法 | |
CN105928833A (zh) | 一种同轴圆筒流变仪流变测试数据的修正方法 | |
CN107367941A (zh) | 基于非线性增益的高超声速飞行器攻角观测方法 | |
CN102778333A (zh) | 一种在大型转动部件上做动平衡测试的方法 | |
CN105478245A (zh) | 基于主轴振动检测的双自由度精密离心机副轴动不平衡量辨识方法 | |
Voisin et al. | Wind tunnel test method to study out-of-service tower crane behaviour in storm winds | |
Hoskoti et al. | Modal analysis of a rotating twisted and tapered Rayleigh beam | |
CN103591919A (zh) | 用于精密离心机静态半径测量的方法与装置 | |
Majhi et al. | Helicopter blade flapping with and without small angle assumption in the presence of dynamic stall | |
RU2344397C2 (ru) | Способ определения демпфирующих свойств моделей самолетов с винтовыми движителями | |
Sytsma et al. | Wind tunnel generated turbulence | |
CN208360517U (zh) | 一种飞行器旋翼组件测试装置 | |
Vermeer et al. | Identification of operational aerofoil state by means of velocity measurements | |
CN106092442B (zh) | Kk轴承法转动惯量测量仪器及其测量方法 | |
Lin et al. | Effects of yaw-roll coupling ratio on lateral-directional aerodynamic characteristics | |
Crowther et al. | Simulation of a spinstabilised sports disc | |
Hohenemser et al. | Model tests on unsteady rotor wake effects. |
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 |