CN109063391A - 旋转条件下的动导数计算检测方法及动导数风洞试验方法 - Google Patents
旋转条件下的动导数计算检测方法及动导数风洞试验方法 Download PDFInfo
- Publication number
- CN109063391A CN109063391A CN201811160514.8A CN201811160514A CN109063391A CN 109063391 A CN109063391 A CN 109063391A CN 201811160514 A CN201811160514 A CN 201811160514A CN 109063391 A CN109063391 A CN 109063391A
- Authority
- CN
- China
- Prior art keywords
- derivative
- indicate
- indicates
- axis
- dynamic
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Algebra (AREA)
- Fluid Mechanics (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Operations Research (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明提供了一种旋转条件下的动导数计算检测方法及动导数风洞试验方法,包括:流场解算步骤:根据需要计算的动导数计算公式,选择相应的确定运动规律,根据所述运动规律,对旋转复合振荡的耦合运动进行流场解算,获得计算结果;动导数数据计算步骤:根据所述计算结果,利用动导数计算公式及相应的运动规律,计算获得周期平均动导数数据。本发明首次设计了考虑旋转角速度的旋转复合振荡的耦合运动,巧妙利用Fourier展开和三角函数系的正交性,建立了一种旋转条件下的动导数计算方法,本发明全面考虑旋转非定常效应,提高了数学模型精度,为解决旋转飞行器关键飞行问题提供了技术基础。
Description
技术领域
本发明涉及动态气动特性领域,具体地,涉及一种旋转条件下的动导数计算检测方法及动导数风洞试验方法。尤其是,涉及一种旋转条件下的动导数计算检测方法。
背景技术
飞行器根据控制模式和用途的不同,经常全程或部分采用绕自身纵轴旋转的飞行方式,转速从几转/秒到上万转/秒。工程实际中,为了向制导控制系统描述飞行器的气动特性,需要对此进行建模。建模时一般将气动力分为定常和非定常两部分,非定常部分由动导数和相应变量的乘积和表述。由于计算手段和认识所限,一般绕自身纵轴旋转的飞行器,在计算动导数时并不计及旋转的影响,而是采用与非旋转飞行器同样的方法获得,这与实际情况是不符合的。动导数的计算偏差会给飞行器动态特性的评估和控制系统的设计带来误差,在关键点的精度降低甚至会导致飞行目的无法达成的严重后果。
发明内容
针对现有技术中的缺陷,本发明的目的是提供一种旋转条件下的动导数计算检测方法及动导数风洞试验方法。
根据本发明提供的一种旋转条件下的动导数计算检测方法,包括:
流场解算步骤:根据需要计算的动导数计算公式,选择相应的确定运动规律,根据所述运动规律,对旋转复合振荡的耦合运动进行流场解算,获得计算结果;
动导数数据计算步骤:根据所述计算结果,利用动导数计算公式及相应的运动规律,计算获得周期平均动导数数据。
优选地,所述计算结果包括以下任一个或任多个结果:密度、压强、速度。
优选地,所述动导数计算公式及相应的运动规律包括:
-周期平均纵向洗流时差导数计算公式及沿z轴的小振幅振荡运动规律;
-周期平均俯仰阻尼导数计算公式及小振幅俯仰振荡运动规律;或者
-周期平均滚转阻尼组合导数计算公式及小振幅滚转振荡运动规律。
优选地,周期平均纵向洗流时差导数计算公式及沿z轴的小振幅振荡运动规律包括:
沿z轴的小振幅振荡运动规律为:
其中
z表示气流坐标轴系oxyz的z坐标;气流坐标轴系oxyz为:原点o:质心;x轴正向:沿飞行速度方向;z轴正向:在飞行器对称面内垂直于ox轴指向下方;y轴正向:按照右手法则确定y轴正向;
z(t)表示t时刻的z坐标值;
t表示时间;
z0表示振荡运动的振幅;
ω1表示振荡运动的圆频率;
表示振荡运动的初始相位;
周期平均纵向洗流时差导数计算公式:
其中
表示的周期平均值,即周期平均纵向洗流时差导数;
表示纵向洗流时差导数,指Cm对的偏导数;
Cm表示绕气流坐标轴系oxyz的oy轴的力矩系数;
My表示绕气流坐标轴系oxyz的oy轴的力矩,即俯仰力矩;
ρ表示空气密度;
V表示飞行器速度;
Sr表示参考面积;
Lr表示参考长度;
表示无量纲的合成攻角变化率,;
表示合成攻角变化率,即合成攻角对时间的导数;
表示减缩频率;
π表示圆周率;
cos αΦ表示αΦ的余弦值;
αΦ表示合成攻角,即纵轴与来流速度的夹角,仅取正;
下标Φ表示一个空间角度,作为合成攻角与攻角α进行区分的标识;
机体坐标轴系ox1y1z1包括:
原点o:质心;x1轴正向:沿纵轴指向来流;y1轴正向:在对称面内,向上;z1轴正向:按照右手法则确定z1轴正向;
Cm(γ)表示指当前滚转角γ下,俯仰力矩系数Cm的值;
γ表示滚转角,即来流速度在横截面的投影与指定子午面的夹角,所述横截面为垂直于纵轴的切面;
γ(t)表示t时刻的滚转角;
ω2表示绕自身纵轴旋转的圆频率,要求ω2=nω1,n为大于2的正整数;
表示滚转角的初始相位;
表示当前滚转角γ下,俯仰力矩系数的值;
表示将Cm对αΦ、q及αΦ、q的变化率进行泰勒展开后的第一项,即除了γ之外的变量均取零值时的Cm;
q表示俯仰角速度。
优选地,俯仰阻尼导数计算公式及小振幅俯仰振荡运动规律包括:
小振幅俯仰振荡运动规律为:
其中
αΦ(t)表示t时刻的瞬时合成攻角αΦ的值;
表示俯仰振荡的初始合成攻角,下标0表示变量合成攻角αΦ的初值,即t=0时的αΦ(t)值;
θ0表示俯仰振荡运动合成攻角的幅值;
周期平均俯仰阻尼导数计算公式:
其中:
表示周期平均俯仰阻尼导数;
表示俯仰阻尼导数,指Cm对的偏导数;
表示无量纲的俯仰角速度;
表示的周期平均值。
优选地,周期平均滚转阻尼组合导数计算公式及小振幅滚转振荡运动规律包括:
小振幅滚转振荡运动规律为:
Φ表示瞬时滚转角,即速度来流在横截面的分量与对称面的夹角,后视图从对称面出发逆时针为正;
Φ0表示滚转振荡的滚转角幅值;
周期平均组合导数的计算公式:
其中
表示的周期平均值,即表示周期平均组合导数;
表示滚转阻尼导数,即Cl对的偏导数;
表示周向洗流时差导数,即Cl对的偏导数;
Cl表示绕纵轴的滚转力矩系数;
表示绕纵轴的滚转力矩;
表示无量纲的滚转角速度;
p表示绕纵轴的滚转角速度;
表示无量纲的滚转角变化率;
表示滚转角变化率;
C1(γ)表示当前滚转角γ下,滚转力矩系数Cl的值;
表示将C1对αΦ、q及αΦ、q的变化率进行泰勒展开后的第一项,即除了γ之外的变量均取零值时的Cl。
根据本发明提供的一种旋转条件下的动导数风洞试验方法,根据上述任一项所述方法获得的动导数数据,进行旋转条件下的动导数风洞试验。
根据本发明提供的一种旋转条件下的动导数风洞试验方法,根据上述任一项所述的旋转条件下的动导数计算检测方法,进行旋转条件下的动导数风洞试验。
优选地,包括:
步骤A:根据需要计算的动导数计算公式,选择相应的运动规律;
步骤B:根据所述相应的运动规律,对旋转复合振荡的耦合运动进行风洞试验。
与现有技术相比,本发明具有如下的有益效果:
1、本发明首次设计了考虑旋转角速度的旋转复合振荡的耦合运动,巧妙利用Fourier展开和三角函数系的正交性,建立了一种旋转条件下的动导数计算方法,本发明全面考虑旋转非定常效应,提高了数学模型精度,为解决旋转飞行器关键飞行问题提供了技术基础。
2、本发明全面考虑旋转流场对绕自身纵轴旋转飞行器动态特性的影响,为解决旋转飞行器动态特性评估提供了手段,提高了非定常气动力地面建模精度,解决关键点性能可能颠覆的问题
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1为根据本发明的优选例中旋转条件下的动导数计算检测方法中坐标系和相关角度定义的后视图示意图
图2为根据本发明的优选例中旋转条件下的动导数计算检测方法中坐标系和相关角度定义的三维示意图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变化和改进。这些都属于本发明的保护范围。
根据本发明提供的一种旋转条件下的动导数计算检测方法,包括:
流场解算步骤:根据需要计算的动导数计算公式,选择相应的确定运动规律,根据所述运动规律,对旋转复合振荡的耦合运动进行流场解算,获得计算结果;
动导数数据计算步骤:根据所述计算结果,利用动导数计算公式及相应的运动规律,计算获得周期平均动导数数据。
进一步地,所述旋转条件为绕自身纵轴旋转。
具体地,所述计算结果包括以下任一个或任多个结果:密度、压强、速度。
具体地,所述动导数计算公式及相应的运动规律包括:
-周期平均纵向洗流时差导数计算公式及沿z轴的小振幅振荡运动规律;
-周期平均俯仰阻尼导数计算公式及小振幅俯仰振荡运动规律;或者
-周期平均滚转阻尼组合导数计算公式及小振幅滚转振荡运动规律。
具体地,周期平均纵向洗流时差导数计算公式及沿z轴的小振幅振荡运动规律包括:
沿z轴的小振幅振荡运动规律为:
其中
z表示气流坐标轴系oxyz的z坐标;气流坐标轴系oxyz为:原点o:质心;x轴正向:沿飞行速度方向;z轴正向:在飞行器对称面内垂直于ox轴指向下方;y轴正向:按照右手法则确定y轴正向;进一步地,气流坐标轴系是自于《GB/T 16638.2-2008空气动力学概念、量和符号》第2部分中的概念。
z(t)表示t时刻的z坐标值;
t表示时间;
z0表示振荡运动的振幅;
ω1表示振荡运动的圆频率;
表示振荡运动的初始相位;
周期平均纵向洗流时差导数计算公式:
其中
表示的周期平均值,即周期平均纵向洗流时差导数;
表示纵向洗流时差导数,指Cm对的偏导数;
Cm表示绕气流坐标轴系oxyz的oy轴的力矩系数;
My表示绕气流坐标轴系oxyz的oy轴的力矩,即俯仰力矩;
ρ表示空气密度;
V表示飞行器速度;
Sr表示参考面积;
Lr表示参考长度;
表示无量纲的合成攻角变化率,;
表示合成攻角变化率,即合成攻角对时间的导数;
表示减缩频率;
π表示圆周率;
cos αΦ表示αΦ的余弦值;
αΦ表示合成攻角,即纵轴与来流速度的夹角,仅取正;
下标Φ表示一个空间角度,作为合成攻角与攻角α进行区分的标识;
机体坐标轴系ox1y1z1包括:
原点o:质心;x1轴正向:沿纵轴指向来流;y1轴正向:在对称面内,向上;z1轴正向:按照右手法则确定z1轴正向;进一步地,所述纵轴平行于机身轴线。进一步地,机体坐标轴系是自于《GB/T 16638.2-2008空气动力学概念、量和符号》第2部分中的概念。
Cm(γ)表示指当前滚转角γ下,俯仰力矩系数Cm的值;
γ表示滚转角,即来流速度在横截面的投影与指定子午面的夹角,所述横截面为垂直于纵轴的切面;
γ(t)表示t时刻的滚转角;
ω2表示绕自身纵轴旋转的圆频率,要求ω2=nω1,n为大于2的正整数;
表示滚转角的初始相位;
表示当前滚转角γ下,俯仰力矩系数的值;
表示将Cm对αΦ、q及αΦ、q的变化率进行泰勒展开后的第一项,即除了γ之外的变量均取零值时的Cm;
q表示俯仰角速度。
具体地,俯仰阻尼导数计算公式及小振幅俯仰振荡运动规律包括:
小振幅俯仰振荡运动规律为:
其中
αΦ(t)表示t时刻的瞬时合成攻角αΦ的值;
表示俯仰振荡的初始合成攻角,下标0表示变量合成攻角αΦ的初值,即t=0时的αΦ(t)值;
θ0表示俯仰振荡运动合成攻角的幅值;
周期平均俯仰阻尼导数计算公式:
其中:
表示周期平均俯仰阻尼导数;
表示俯仰阻尼导数,指Cm对的偏导数;
表示无量纲的俯仰角速度;进一步地,俯仰角速度是自于《GB/T 16638.2-2008空气动力学概念、量和符号》的名词,本领域技术人员参照现有技术可以得出俯仰角速度,在此不予赘述。
表示的周期平均值。
具体地,周期平均滚转阻尼组合导数计算公式及小振幅滚转振荡运动规律包括:
小振幅滚转振荡运动规律为:
Φ表示瞬时滚转角,即速度来流在横截面的分量与对称面的夹角,后视图从对称面出发逆时针为正;
Φ0表示滚转振荡的滚转角幅值;
周期平均组合导数的计算公式:
其中
表示的周期平均值,即表示周期平均组合导数;
表示滚转阻尼导数,即Cl对的偏导数;
表示周向洗流时差导数,即Cl对的偏导数;
Cl表示绕纵轴的滚转力矩系数;
表示绕纵轴的滚转力矩;
表示无量纲的滚转角速度;
p表示绕纵轴的滚转角速度;
表示无量纲的滚转角变化率;
表示滚转角变化率;
Cl(γ)表示当前滚转角γ下,滚转力矩系数Cl的值;
表示将Cl对αΦ、q及αΦ、q的变化率进行泰勒展开后的第一项,即除了γ之外的变量均取零值时的Cl。
根据本发明提供的一种旋转条件下的动导数风洞试验方法,根据上述任一项所述方法获得的动导数数据,进行旋转条件下的动导数风洞试验。
根据本发明提供的一种旋转条件下的动导数风洞试验方法,根据上述任一项所述的旋转条件下的动导数计算检测方法,进行旋转条件下的动导数风洞试验。
具体地,包括:
步骤A:根据需要计算的动导数计算公式,选择相应的运动规律;
步骤B:根据所述相应的运动规律,对旋转复合振荡的耦合运动进行风洞试验。
下面通过优选例,对本发明进行更为具体的说明。
实施例1:
本发明一种旋转条件下的动导数计算方法包括以下步骤:
(1)基于Fourier展开和旋转流场特性,得到旋转条件下的动导数计算公式;
(2)根据步骤(1)中得到的旋转条件下的动导数计算公式,选择合适的振荡频率,采用步骤(1)中的运动规律,对旋转+振荡的耦合运动进行流场解算,获得计算结果;
(3)根据步骤(2)中得到的输入,按照步骤S1中旋转条件下的动导数计算公式求得动导数数据。
如图1、图2所示,是旋转条件下的动导数计算公式方法的坐标系和相关角度的定义,包括气流坐标轴系oxyz、机体坐标轴系ox1y1z1、弹体10、鸭舵20和尾翼30,其中有2个鸭舵和4个尾翼。
本实施例中涉及的定义如下:
坐标系:
气流坐标轴系oxyz——原点:质心,x轴正向:沿飞行速度方向,z轴正向:在飞行器对称面内垂直于ox轴指向下方,y轴正向:按照右手法则确定;
机体坐标轴系ox1y1z1——原点:质心,x1轴正向:沿纵轴指向来流,y1轴正向:在对称面内,向上,z1轴正向:按照右手法则确定;
符号:
z——气流坐标轴系oxyz的z坐标;
t——时间;
z0——振荡运动的振幅;
ω1——振荡运动的圆频率;
——振荡运动的初始相位;
的周期平均值;
——纵向洗流时差导数,Cm对的偏导数;
Cm——绕气流坐标轴系oxyz的oy轴的力矩系数,
My——绕气流坐标轴系oxyz的oy轴的力矩,即俯仰力矩;
——将Cm对αΦ、q及其变化率进行Taylor展开后的第一项,即除了γ之外的变量均取零值时的Cm;
q——俯仰角速度;
ρ——空气密度;
V——飞行器速度;
Sr——参考面积;
Lr——参考长度;
——减缩频率,
——无量纲的合成攻角变化率,
——合成攻角变化率;
αΦ——合成攻角,纵轴与来流速度的夹角,仅取正;
π——圆周率;
γ——滚转角,即来流速度在横截面的投影与指定子午面的夹角,所述横截面为垂直于纵轴的切面,对于绕自身纵轴旋转的飞行器,可表示为如下公式:
ω2——绕自身纵轴旋转的圆频率,要求ω2=nω1,n为大于2的正整数;
——滚转角的初始相位;
αΦ(t)——瞬时合成攻角;
——俯仰振荡的初始合成攻角;
θ0——俯仰振荡运动合成攻角的幅值;
的周期平均值;
——俯仰阻尼导数,Cm对的偏导数;
一一无量纲的俯仰角速度;
Φ——瞬时滚转角;
Φ0——滚转振荡的滚转角幅值;
的周期平均值;
——滚转阻尼导数,Cl对的偏导数;
——周向洗流时差导数,Cl对的偏导数;
Cl——绕纵轴的滚转力矩系数,
——绕纵轴的滚转力矩;
——无量纲的滚转角速度,
p——绕纵轴的滚转角速度;
——无量纲的滚转角变化率,
——滚转角变化率。
本实施例的一种旋转条件下的动导数计算方法包括以下步骤:
(1)基于Fourier展开和旋转流场特性,得到旋转条件下的动导数计算公式:
a)纵向洗流时差导数
采用沿z轴的小振幅振荡运动规律:
周期平均纵向洗流时差导数的计算公式:
b)俯仰阻尼导数
采用如下小振幅俯仰振荡运动规律:
周期平均组合导数的计算公式:
周期平均俯仰阻尼导数的计算公式:
c)滚转阻尼组合导数
采用如下小振幅滚转振荡运动规律:
周期平均组合导数的计算公式:
(2)根据步骤(1)中得到的旋转条件下的动导数计算公式,考虑飞行器实际的旋转角速度,选择合适的振荡频率,特别要注意ω2应是ω1三倍频以上的关系;然后,从需要计算的动导数类型出发,选择步骤(1)中所描述的运动规律,用合适的计算方法对旋转+振荡的耦合运动进行流场解算,获得计算结果。
(3)根据步骤(2)中得到的计算结果,按照步骤(1)中旋转条件下的动导数计算公式,选取一个周期的数值进行积分,即可求得动导数数据。特别要说明的是,通过CFD求解获得的流场结果,可用类似的方法获得其他类型的动导数,只是由于在工程中量值较小,一般予以忽略,故本发明中未予以一一列举。
根据需求动导数的类型,选择相应的运动规律,实施必需的计算或风洞试验,依据结果积分获得动导数数据。本发明可应用于绕自身纵轴旋转的各类飞行器。
在本申请的描述中,需要理解的是,术语“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本申请和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本申请的限制。
本领域技术人员知道,除了以纯计算机可读程序代码方式实现本发明提供的系统、装置及其各个模块以外,完全可以通过将方法步骤进行逻辑编程来使得本发明提供的系统、装置及其各个模块以逻辑门、开关、专用集成电路、可编程逻辑控制器以及嵌入式微控制器等的形式来实现相同程序。所以,本发明提供的系统、装置及其各个模块可以被认为是一种硬件部件,而对其内包括的用于实现各种程序的模块也可以视为硬件部件内的结构;也可以将用于实现各种功能的模块视为既可以是实现方法的软件程序又可以是硬件部件内的结构。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变化或修改,这并不影响本发明的实质内容。在不冲突的情况下,本申请的实施例和实施例中的特征可以任意相互组合。
Claims (9)
1.一种旋转条件下的动导数计算检测方法,其特征在于,包括:
流场解算步骤:根据需要计算的动导数计算公式,选择相应的确定运动规律,根据所述运动规律,对旋转复合振荡的耦合运动进行流场解算,获得计算结果;
动导数数据计算步骤:根据所述计算结果,利用动导数计算公式及相应的运动规律,计算获得周期平均动导数数据。
2.根据权利要求1所述的旋转条件下的动导数计算检测方法,其特征在于,所述计算结果包括以下任一个或任多个结果:密度、压强、速度。
3.根据权利要求1所述的旋转条件下的动导数计算检测方法,其特征在于,所述动导数计算公式及相应的运动规律包括:
-周期平均纵向洗流时差导数计算公式及沿z轴的小振幅振荡运动规律;
-周期平均俯仰阻尼导数计算公式及小振幅俯仰振荡运动规律;或者
-周期平均滚转阻尼组合导数计算公式及小振幅滚转振荡运动规律。
4.根据权利要求3所述的旋转条件下的动导数计算检测方法,其特征在于,周期平均纵向洗流时差导数计算公式及沿z轴的小振幅振荡运动规律包括:
沿z轴的小振幅振荡运动规律为:
其中
z表示气流坐标轴系oxyz的z坐标;气流坐标轴系oxyz为:原点o:质心;x轴正向:沿飞行速度方向;z轴正向:在飞行器对称面内垂直于ox轴指向下方;y轴正向:按照右手法则确定y轴正向;
z(t)表示t时刻的z坐标值;
t表示时间;
z0表示振荡运动的振幅;
ω1表示振荡运动的圆频率;
表示振荡运动的初始相位;
周期平均纵向洗流时差导数计算公式:
其中
表示的周期平均值,即周期平均纵向洗流时差导数;
表示纵向洗流时差导数,指Cm对的偏导数;
Cm表示绕气流坐标轴系oxyz的oy轴的力矩系数;
My表示绕气流坐标轴系oxyz的oy轴的力矩,即俯仰力矩;
ρ表示空气密度;
V表示飞行器速度;
Sr表示参考面积;
Lr表示参考长度;
表示无量纲的合成攻角变化率,;
表示合成攻角变化率,即合成攻角对时间的导数;
表示减缩频率;
π表示圆周率;
cos αΦ表示αΦ的余弦值;
αΦ表示合成攻角,即纵轴与来流速度的夹角,仅取正;
下标Φ表示一个空间角度,作为合成攻角与攻角α进行区分的标识;
机体坐标轴系ox1y1z1包括:
原点o:质心;x1轴正向:沿纵轴指向来流;y1轴正向:在对称面内,向上;z1轴正向:按照右手法则确定z1轴正向;
Cm(γ)表示指当前滚转角γ下,俯仰力矩系数Cm的值;
γ表示滚转角,即来流速度在横截面的投影与指定子午面的夹角,所述横截面为垂直于纵轴的切面;
γ(t)表示t时刻的滚转角;
ω2表示绕自身纵轴旋转的圆频率,要求ω2=nω1,n为大于2的正整数;
表示滚转角的初始相位;
表示当前滚转角γ下,俯仰力矩系数的值;
表示将Cm对αΦ、q及αΦ、q的变化率进行泰勒展开后的第一项,即除了γ之外的变量均取零值时的Cm;
q表示俯仰角速度。
5.根据权利要求4所述的旋转条件下的动导数计算检测方法,其特征在于,俯仰阻尼导数计算公式及小振幅俯仰振荡运动规律包括:
小振幅俯仰振荡运动规律为:
其中
αΦ(t)表示t时刻的瞬时合成攻角αΦ的值;
表示俯仰振荡的初始合成攻角,下标0表示变量合成攻角αΦ的初值,即t=0时的αΦ(t)值;
θ0表示俯仰振荡运动合成攻角的幅值;
周期平均俯仰阻尼导数计算公式:
其中:
表示周期平均俯仰阻尼导数;
表示俯仰阻尼导数,指Cm对的偏导数;
表示无量纲的俯仰角速度;
表示的周期平均值。
6.根据权利要求5所述的旋转条件下的动导数计算检测方法,其特征在于,周期平均滚转阻尼组合导数计算公式及小振幅滚转振荡运动规律包括:
小振幅滚转振荡运动规律为:
Φ表示瞬时滚转角,即速度来流在横截面的分量与对称面的夹角,后视图从对称面出发逆时针为正;
Φ0表示滚转振荡的滚转角幅值;
周期平均组合导数的计算公式:
其中
表示的周期平均值,即表示周期平均组合导数;
表示滚转阻尼导数,即Cl对的偏导数;
表示周向洗流时差导数,即Cl对的偏导数;
Cl表示绕纵轴的滚转力矩系数;
表示绕纵轴的滚转力矩;
表示无量纲的滚转角速度;
p表示绕纵轴的滚转角速度;
表示无量纲的滚转角变化率;
表示滚转角变化率;
Cl(γ)表示当前滚转角γ下,滚转力矩系数Cl的值;
表示将Cl对αΦ、q及αΦ、q的变化率进行泰勒展开后的第一项,即除了γ之外的变量均取零值时的Cl。
7.一种旋转条件下的动导数风洞试验方法,其特征在于,根据权利要求1至6中任一项所述方法获得的动导数数据,进行旋转条件下的动导数风洞试验。
8.一种旋转条件下的动导数风洞试验方法,其特征在于,根据权利要求1至6中任一项所述的旋转条件下的动导数计算检测方法,进行旋转条件下的动导数风洞试验。
9.根据权利要求8所述的旋转条件下的动导数风洞试验方法,其特征在于,包括:
步骤A:根据需要计算的动导数计算公式,选择相应的运动规律;
步骤B:根据所述相应的运动规律,对旋转复合振荡的耦合运动进行风洞试验。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811160514.8A CN109063391B (zh) | 2018-09-30 | 2018-09-30 | 旋转条件下的动导数计算检测方法及动导数风洞试验方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811160514.8A CN109063391B (zh) | 2018-09-30 | 2018-09-30 | 旋转条件下的动导数计算检测方法及动导数风洞试验方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109063391A true CN109063391A (zh) | 2018-12-21 |
CN109063391B CN109063391B (zh) | 2020-10-09 |
Family
ID=64767393
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811160514.8A Active CN109063391B (zh) | 2018-09-30 | 2018-09-30 | 旋转条件下的动导数计算检测方法及动导数风洞试验方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109063391B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111473949A (zh) * | 2020-04-14 | 2020-07-31 | 日照坤仑智能科技有限公司 | 一种测量飞行器掉翼、滚摇和洗流时差的装置及方法 |
CN111651833A (zh) * | 2020-05-11 | 2020-09-11 | 上海机电工程研究所 | 一种旋转类飞行器流场分析方法及系统 |
CN114088330A (zh) * | 2021-11-19 | 2022-02-25 | 中国航空工业集团公司哈尔滨空气动力研究所 | 一种纵向动导数低速风洞连续测量试验方法 |
CN114608786A (zh) * | 2022-05-11 | 2022-06-10 | 中国空气动力研究与发展中心设备设计与测试技术研究所 | 一种飞行器动导数试验数据处理方法 |
CN115343012A (zh) * | 2022-07-07 | 2022-11-15 | 中国航空工业集团公司哈尔滨空气动力研究所 | 一种非稳态大幅振荡试验方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2189962A1 (en) * | 2007-08-16 | 2010-05-26 | Indra Sistemas, S.A. | Method for real-time simulation of a helicopter rotor |
CN103400035A (zh) * | 2013-07-29 | 2013-11-20 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种高可信度快速预测飞行器滚转动导数的方法 |
CN105182989A (zh) * | 2015-09-25 | 2015-12-23 | 北京航空航天大学 | 一种风场影响下的飞机姿态控制方法 |
CN106227971A (zh) * | 2016-08-03 | 2016-12-14 | 中国人民解放军63821部队 | 基于谐波平衡法的动导数快速预测技术 |
US20170046968A1 (en) * | 2015-08-11 | 2017-02-16 | The Boeing Company | Flight simulation modeling of aircraft dynamic stall aerodynamics |
CN107357976A (zh) * | 2017-06-27 | 2017-11-17 | 四川腾盾科技有限公司 | 一种飞行器的动导数的计算方法 |
CN108256264A (zh) * | 2018-02-08 | 2018-07-06 | 北京航空航天大学 | 一种基于地面频响试验的气动伺服弹性稳定性预测方法 |
-
2018
- 2018-09-30 CN CN201811160514.8A patent/CN109063391B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2189962A1 (en) * | 2007-08-16 | 2010-05-26 | Indra Sistemas, S.A. | Method for real-time simulation of a helicopter rotor |
CN103400035A (zh) * | 2013-07-29 | 2013-11-20 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种高可信度快速预测飞行器滚转动导数的方法 |
US20170046968A1 (en) * | 2015-08-11 | 2017-02-16 | The Boeing Company | Flight simulation modeling of aircraft dynamic stall aerodynamics |
CN105182989A (zh) * | 2015-09-25 | 2015-12-23 | 北京航空航天大学 | 一种风场影响下的飞机姿态控制方法 |
CN106227971A (zh) * | 2016-08-03 | 2016-12-14 | 中国人民解放军63821部队 | 基于谐波平衡法的动导数快速预测技术 |
CN107357976A (zh) * | 2017-06-27 | 2017-11-17 | 四川腾盾科技有限公司 | 一种飞行器的动导数的计算方法 |
CN108256264A (zh) * | 2018-02-08 | 2018-07-06 | 北京航空航天大学 | 一种基于地面频响试验的气动伺服弹性稳定性预测方法 |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111473949A (zh) * | 2020-04-14 | 2020-07-31 | 日照坤仑智能科技有限公司 | 一种测量飞行器掉翼、滚摇和洗流时差的装置及方法 |
CN111651833A (zh) * | 2020-05-11 | 2020-09-11 | 上海机电工程研究所 | 一种旋转类飞行器流场分析方法及系统 |
CN114088330A (zh) * | 2021-11-19 | 2022-02-25 | 中国航空工业集团公司哈尔滨空气动力研究所 | 一种纵向动导数低速风洞连续测量试验方法 |
CN114088330B (zh) * | 2021-11-19 | 2023-03-17 | 中国航空工业集团公司哈尔滨空气动力研究所 | 一种纵向动导数低速风洞连续测量试验方法 |
CN114608786A (zh) * | 2022-05-11 | 2022-06-10 | 中国空气动力研究与发展中心设备设计与测试技术研究所 | 一种飞行器动导数试验数据处理方法 |
CN115343012A (zh) * | 2022-07-07 | 2022-11-15 | 中国航空工业集团公司哈尔滨空气动力研究所 | 一种非稳态大幅振荡试验方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109063391B (zh) | 2020-10-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109063391A (zh) | 旋转条件下的动导数计算检测方法及动导数风洞试验方法 | |
CN106249745B (zh) | 四轴无人机的控制方法 | |
CN105760587B (zh) | 一种双轴非定常气动力建模方法及装置 | |
CN106933104A (zh) | 一种基于dic‑pid的四旋翼飞行器姿态与位置的混合控制方法 | |
CN103940442B (zh) | 一种采用加速收敛算法的定位方法及装置 | |
CN109614633A (zh) | 一种复合式旋翼飞行器非线性建模方法及配平方法 | |
CN103307938B (zh) | 一种旋转弹气动参数获取方法 | |
CN110162933B (zh) | 一种共轴多旋翼仿真方法及系统 | |
CN109612676A (zh) | 基于飞行试验数据的气动参数反算方法 | |
CN107992063B (zh) | 基于变参数章动阻尼的变速倾侧动量轮进动控制方法 | |
CN104281150A (zh) | 一种姿态机动的轨迹规划方法 | |
CN106527122A (zh) | 平流层飞艇定高飞行非线性pid控制方法 | |
CN106248082B (zh) | 一种飞行器自主导航系统及导航方法 | |
CN105043414A (zh) | 一种三轴惯性稳定平台系统的台体控制参数计算方法 | |
CN107357976A (zh) | 一种飞行器的动导数的计算方法 | |
Song et al. | New intermediate quaternion based control of spacecraft: part I—almost global attitude tracking | |
CN102506862B (zh) | 基于2阶不可交换误差补偿模型的圆锥补偿方法 | |
CN106774385A (zh) | 一种采用自适应变结构的飞艇定点悬停控制方法 | |
Millane et al. | Range-inertial estimation for airborne wind energy | |
Li et al. | Simulation method for wind tunnel based virtual flight testing | |
CN109062242A (zh) | 一种新的旋翼无人机控制方法 | |
CN108508905A (zh) | 一种基于最短空间轴的姿态机动控制和导引律规划方法 | |
Gardecki et al. | An adequate mathematical model of four-rotor flying robot in the context of control simulations | |
CN105905296B (zh) | 基于自旋稳定的单翼回旋飞行器升力优化设计方法 | |
Wang et al. | Influence of yawing force frequency on angular motion and ballistic characteristics of a dual-spin projectile |
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 |