CN105629725B - 一种后缘舵滑翔飞行器的弹性运动建模方法 - Google Patents
一种后缘舵滑翔飞行器的弹性运动建模方法 Download PDFInfo
- Publication number
- CN105629725B CN105629725B CN201410602382.5A CN201410602382A CN105629725B CN 105629725 B CN105629725 B CN 105629725B CN 201410602382 A CN201410602382 A CN 201410602382A CN 105629725 B CN105629725 B CN 105629725B
- Authority
- CN
- China
- Prior art keywords
- rudder
- force
- disturbance
- angle
- vibration
- 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.)
- Active
Links
Landscapes
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明属于飞行器弹性运动建模技术领域,具体涉及一种后缘舵滑翔飞行器的弹性运动建模方法。该方法包括如下步骤:(1)利用模态正交性可把弹体横向自由振动转化为各个相互独立的主振动的叠加;(2)根据线性小扰动假设,在外力作用下,弹体的横向振动仍可近似用模态叠加来描述,振型函数由弹体结构特性(刚度和质量分布)和弹体的边界条件确定;(3)根据达朗伯原理等方法建立振动微分方程;(4)分析影响滑翔飞行器弹性振动的外力:包括气动力,舵面控制力,姿态喷管控制力以及舵面摆动的惯性力等。本发明对传统弹性运动建模方法进行修正和完善。
Description
技术领域
本发明属于飞行器弹性运动建模技术领域,具体涉及一种后缘舵滑翔飞行器的弹性运动建模方法。
背景技术
随着飞行速度和机动性的不断提高,飞行器弹性结构与气动、控制系统耦合的动力学问题受到日益重视。作为控制对象的弹性飞行器,姿态控制系统的敏感元件除了感应刚体运动外,还感受弹性结构的振动。用于描述飞行器弹性振动的弹性运动方程,其模型的正确性关乎姿控系统设计的成败。而传统的简化弹性运动建模方法(参考文献:《导弹与航天丛书——控制系统(上)》,2.3.4节),对于后缘舵控制的新一代高超声速面对称飞行器,已不再适用。
发明内容
本发明的目的在于提供一种后缘舵滑翔飞行器的弹性运动建模方法,对传统弹性运动建模方法进行修正和完善。
为达到上述目的,本发明所采取的技术方案为:
一种后缘舵滑翔飞行器的弹性运动建模方法,包括如下步骤:(1)利用模态正交性把弹体横向自由振动转化为各个相互独立的主振动的叠加;(2)根据线性小扰动假设,在外力作用下,弹体的横向振动用模态叠加来描述,振型函数由弹体刚度和质量分布、弹体的边界条件确定;(3)根据达朗伯原理建立振动微分方程;(4)分析影响滑翔飞行器弹性振动的外力:包括气动力,舵面控制力,姿态喷管控制力以及舵面摆动的惯性力。
所述的步骤(1)具体为:利用模态正交性把弹体横向自由振动转化为各个相互独立的主振动的叠加,即设
Wi(X)即弹体弹性的第i次固有振型函数,Φ=[W1(X),W2(X),…,Wn(X)],qi(t)为第i次固有振型对应的广义坐标。
所述的步骤(3)具体为:根据达朗伯原理建立振动微分方程为:
把式(1)代入式(2),并前乘ΦT得
式中,Μp=ΦTMΦ,Cp=ΦTCΦ,Kp=ΦTKΦ分别为广义模态质量矩阵、模态阻尼矩阵和模态刚度矩阵,都为对角矩阵,Fp=ΦTf为广义力阵,
式(3)写成
其中ωi——第i次振型的固有频率;
ξi——第i次振型的阻尼系数;
Fi——对应第i次振型的广义力;
Mi——对应第i次振型的广义质量;
且有
式中是沿弹体纵轴受到的外力在弹体坐标系轴Y1上的投影;m(X)是沿弹体纵轴的弹体质量分布;lK为飞行器纵轴总长度;ωi、ξi和Wi(X)是系统的固有属性,由数值方法求得并通过模态试验校验。
所述的步骤(4)具体为:对于俯仰通道,在弹体变形情况下,气动力沿纵轴各点局部攻角不同,气动力须看作沿纵轴的分布力,除弹体气动力外,其它外力都作用在弹体的某一位置上,是集中力;
影响弹体分布气动力的飞行姿态扰动包括攻角扰动Δα、俯仰角速度扰动舵偏角扰动角速度扰动及角加速度扰动将引起舵面控制力、摆动惯性力变化;姿控动力发动机的启用与关闭也产生控制力扰动;
用D1i和D2i分别表征与俯仰角速度扰动和攻角小扰动Δα成比例的广义气动力对第i次振型振动的影响;根据广义力定义和气动力计算方法推导出D1i和D2i的表达式(7)和(8):
式中q——飞行动压;V——飞行速度;
SM——飞行器参考面积;
lK——飞行器纵轴总长度;
——沿弹体纵轴分布的局部法向升力系数(单位:/°);
xZ——飞行器质心距弹头顶点的距离;
Wi(xn)——xn处对应的第i次振型。
后缘升降舵局部坐标系ORXRYRZR定义为:OR为舵轴中心线与控制舵交点处,ORXR平行弹体纵轴指向前,ORYR平行弹体坐标系Y1轴指向上,ORZR沿舵轴中心线,方向由右手法则确定。假设该升降舵偏角为对单个舵受力进行分析:
①舵面气动力和力矩
对小扰动进行线性化,假设扰动前的舵偏角为即
则舵偏小扰动产生的气动力在Y1轴方向的投影为
相应的气动力矩变化量为
其中——单个升降舵法向力系数导数(/°);
——单个升降舵绕ZR轴的气动力矩系数导数(/°),所有气动系数均沿局部坐标系各轴正方向定义;
②舵面惯性力和力矩
舵面质心处的惯性力与弹道轴向过载nx、法向过载ny以及舵偏角速度角加速度有关,将惯性力分解到ORYR轴和ORXR轴两个方向上;
忽略二阶小量,得舵偏小扰动产生的惯性力在Y1轴方向的投影为
相应的惯性力矩变化量
式中mR——单个升降舵质量;
lR——升降舵质心到舵轴的距离;
JR——单个升降舵相对于舵轴的转动惯量;
③伺服机构与舵系统的支反力
当飞行器升降舵采用直线式机电伺服机构驱动时,升降舵轴处受到的扭矩是通过摇臂与伺服机构作用力产生的力矩来平衡,此时舵面受到的外力矩对控制舱后舱体产生的作用力Fh在Y1轴方向的投影为
其中,lRh——升降舵支架到伺服机构后支点的距离;
则舵偏小扰动产生的后支点作用力在Y1轴方向的投影为
升降舵支架处舱体受到的作用力在Y1方向的投影与fhY大小相等,方向相反;
④与舵偏有关的方程式系数
式(10)、(12)和(15)给出了舵偏小扰动引起的作用于弹上的外力在Y1轴上的投影,这些外力均是集中力,忽略差动舵偏角影响时,代入公式(5)得广义力为
式中——升降舵的个数;
若考虑初始差动舵偏角影响,即右舵和左舵的舵偏角δ1和δ2为
考虑初始差动舵偏影响,推导得相应的方程式系数为
假设俯仰通道姿控发动机推力为fRCS,即沿Y1轴方向的投影,则相应的有
当俯仰通道姿控喷管产生的控制力在飞行器体坐标系内为正时,ΔδZK取+1,反之ΔδZKΔδRCS取-1,俯仰喷管不工作时ΔδZKΔδRCS取0;
综合上述与小扰动有关的广义力项,则俯仰通道的弹性振动方程为
偏航通道飞行器的控制是通过方向舵和姿控喷管,其偏航通道扰动运动方程形式同公式(22),其中的飞行姿态扰动换成了偏航角速度扰动侧滑角扰动Δβ、方向舵偏角扰动Δδψ和角加速度扰动
对于滚转通道弹性运动建模方法,其飞行姿态扰动换成了差动舵偏角扰动Δδγ和角加速度扰动舵面差动小扰动产生的绕弹体系X1轴的气动滚转力矩为
设左舵和右舵舵面惯性力的法向分量变化量为f2IY和f1IY,正方向均定义为沿Y1轴正方向,则产生滚转力矩为:
伺服机构后支点处产生的滚转力矩为
前支点处产生的滚转力矩大小相等,方向相反;
则滚转通道弹性振动方程中与舵面操纵有关的广义力为
与姿控喷管有关的方程式系数为
其中MRCS——滚转通道的姿控喷管控制力矩;
Mγ——扭转振型对应的广义质量;
Qγ(xRCS)——姿控喷管xRCS处的扭转振型;
最终推导得到飞行器滚转通道弹性振动方程为
其中与舵面操纵机构有关的方程式系数:
式中nγ——滚转通道控制舵面个数,这里指飞行器升降舵的个数;
——右舵、左舵相对弹体轴X1的滚转力矩系数对差动舵偏的导数;
ZR——升降舵质心距弹体中心对称面的距离;
δγ0——扰动前的等效俯仰舵偏角、差动舵偏角(rad);
Qγ(xR)——舵轴位置xR处的扭转振型;
Zh——伺服机构后支点距弹体中心对称面的距离;
Qγ(xh)——后支点xh处的扭转振型;
lRh——升降舵支架到伺服机构后支点的距离。
本发明所取得的有益效果为:
本发明针对某滑翔飞行器的特点,全面探讨影响该飞行器弹性振动的各种广义力,随后根据小扰动线性化假设、达朗伯原理和模态叠加等方法,建立滑翔飞行器的弹性振动方程,并推导出方程中各项系数的表达式,最终形成了一套适用范围更广、可应用于后缘控制舵滑翔飞行器的弹性运动建模方法。
本发明应用滑翔飞行器模态试验结果,对推导出的弹性运动方程各项系数开展灵敏度分析,总结规律,确定弹性振动方程中各项系数的主要影响因素,为飞行器弹性运动方程的建立和姿控系统设计提供支撑。
本发明有效地保证了型号任务进度,为滑翔飞行器的弹性运动建模和姿控系统设计提供了有力支撑,为后续型号设计奠定了基础。
本发明通过开展灵敏度分析,提高了对高超声速飞行器弹性运动建模的认识,确定了弹性振动方程中各项系数的主要影响因素;分析采用改进方法计算的弹性方程式系数与传统方法计算结果的差异来源,明确了后缘舵的气动铰链力矩项对弹性方程的影响。
附图说明
图1为某后缘舵滑翔飞行器弹体坐标系O1X1Y1Z1定义示意图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
对于战略导弹和运载火箭,在研究和设计其姿态控制系统时,必须考虑弹体的弹性变形影响。本发明所述后缘舵滑翔飞行器的弹性运动建模方法为:
(1)利用模态正交性可把弹体横向自由振动转化为各个相互独立的主振动的叠加;(2)根据线性小扰动假设,在外力作用下,弹体的横向振动仍可近似用模态叠加来描述,振型函数由弹体结构特性(刚度和质量分布)和弹体的边界条件确定;(3)根据达朗伯原理等方法建立振动微分方程;(4)分析影响滑翔飞行器弹性振动的外力:包括气动力,舵面控制力,姿态喷管控制力以及舵面摆动的惯性力等。
在研究和设计运载火箭等飞行器的姿态控制系统时,必须考虑弹体的弹性变形影响。利用模态正交性可把弹体横向自由振动转化为各个相互独立的主振动的叠加。即设
Wi(X)即弹体弹性的第i次固有振型函数,Φ=[W1(X),W2(X),…,Wn(X)]。qi(t)为第i次固有振型对应的广义坐标。
根据达朗伯原理等方法可建立振动微分方程为:
把式(31)代入式(32),并前乘ΦT得
式中,Μp=ΦTMΦ,Cp=ΦTCΦ,Kp=ΦTKΦ分别为广义模态质量矩阵、模态阻尼矩阵和模态刚度矩阵,它们都为对角矩阵。Fp=ΦTf为广义力阵。
式(33)可写成
其中ωi——第i次振型的固有频率;
ξi——第i次振型的阻尼系数;
Fi——对应第i次振型的广义力;
Mi——对应第i次振型的广义质量。
且有
式中是沿弹体纵轴受到的外力在弹体坐标系轴Y1上的投影;m(X)是沿弹体纵轴的弹体质量分布;lK为飞行器纵轴总长度。ωi、ξi和Wi(X)是系统的固有属性,可由数值方法求得并通过模态试验校验。
下面以俯仰通道为例,讨论式(35)中的形式。首先,在弹体变形情况下,气动力沿纵轴各点局部攻角不同,气动力须看作沿纵轴的分布力。除弹体气动力外,其它外力都作用在弹体的某一位置上,是集中力。
影响弹体分布气动力的飞行姿态扰动包括攻角扰动Δα、俯仰角速度扰动另外,舵偏角扰动角速度扰动及角加速度扰动将引起舵面控制力、摆动惯性力等变化。姿控动力发动机的启用与关闭也将产生控制力扰动。
通常用D1i和D2i分别表征与俯仰角速度扰动和攻角小扰动Δα成比例的广义气动力对第i次振型振动的影响。根据广义力定义和气动力计算方法可推导出D1i和D2i的表达式,列于式(7)和(8)。
式中q——飞行动压;V——飞行速度;
SM——飞行器参考面积;
lK——飞行器纵轴总长度;
——沿弹体纵轴分布的局部法向升力系数(单位:/°);
xZ——飞行器质心距弹头顶点的距离;
Wi(xn)——xn处对应的第i次振型。
下面重点探讨与舵面操纵机构有关的广义力。后缘升降舵局部坐标系ORXRYRZR定义为:OR为舵轴中心线与控制舵交点处,ORXR平行弹体纵轴指向前,ORYR平行弹体坐标系Y1轴指向上,ORZR沿舵轴中心线,方向由右手法则确定。假设该升降舵偏角为对单个舵受力进行分析。
①舵面气动力和力矩
对小扰动进行线性化,假设扰动前的舵偏角为即
则舵偏小扰动产生的气动力在Y1轴方向的投影为
相应的气动力矩变化量为
其中——单个升降舵法向力系数导数(/°);
——单个升降舵绕ZR轴的气动力矩系数导数(/°),所有气动系数均沿局部坐标系各轴正方向定义。
②舵面惯性力和力矩
舵面质心处的惯性力与弹道轴向过载nx、法向过载ny以及舵偏角速度角加速度有关,将惯性力分解到ORYR轴和ORXR轴两个方向上。
忽略二阶小量,得舵偏小扰动产生的惯性力在Y1轴方向的投影为
相应的惯性力矩变化量
式中mR——单个升降舵质量;
lR——升降舵质心到舵轴的距离;
JR——单个升降舵相对于舵轴的转动惯量。
③伺服机构与舵系统的支反力
值得注意的是,当XX飞行器升降舵采用直线式机电伺服机构驱动时,升降舵轴处受到的扭矩是通过摇臂与伺服机构作用力产生的力矩来平衡。此时舵面受到的外力矩对控制舱后舱体产生的作用力Fh在Y1轴方向的投影为
其中,lRh——升降舵支架到伺服机构后支点的距离。
则舵偏小扰动产生的后支点作用力在Y1轴方向的投影为
升降舵支架处舱体受到的作用力在Y1方向的投影与fhY大小相等,方向相反。
④与舵偏有关的方程式系数
式(40)、(42)和(45)给出了舵偏小扰动引起的作用于弹上的外力在Y1轴上的投影,这些外力均是集中力,忽略差动舵偏角影响时,代入公式(35)得广义力为
式中——升降舵的个数。
若考虑初始差动舵偏角影响,即右舵和左舵的舵偏角δ1和δ2为
推导得相应的方程式系数(考虑初始差动舵偏影响)为
对于XX飞行器,俯仰通道的控制除了升降舵,通常还有姿控动力发动机。假设俯仰通道姿控发动机推力为fRCS(沿Y1轴方向的投影),则式(22)中相应的有
当俯仰通道姿控喷管产生的控制力在飞行器体坐标系内为正时,ΔδZK取+1,反之ΔδZK取-1,俯仰喷管不工作时ΔδZK取0。
综合上述与小扰动有关的广义力项,则俯仰通道的弹性振动方程为
偏航通道XX飞行器的控制通常是通过方向舵和姿控喷管,因此其偏航通道扰动运动方程与俯仰类似,形式同公式(52),只是其中的飞行姿态扰动换成了偏航角速度扰动侧滑角扰动Δβ、方向舵偏角扰动Δδψ和角加速度扰动
滚转通道弹性运动建模方法与俯仰通道类似,只是其飞行姿态扰动换成了差动舵偏角扰动Δδγ和角加速度扰动舵面差动小扰动产生的绕弹体系X1轴的气动滚转力矩为
设左舵和右舵舵面惯性力的法向分量变化量为f2IY和f1IY(正方向均定义为沿Y1轴正方向),则产生滚转力矩为
伺服机构后支点处产生的滚转力矩为
前支点处产生的滚转力矩大小相等,方向相反。
则滚转通道弹性振动方程中与舵面操纵有关的广义力为
与姿控喷管有关的方程式系数为
其中MRCS——滚转通道的姿控喷管控制力矩;
Mγ——扭转振型对应的广义质量;
Qγ(xRCS)——姿控喷管xRCS处的扭转振型。
最终推导得到飞行器滚转通道弹性振动方程为
其中与舵面操纵机构有关的方程式系数:
式中nγ——滚转通道控制舵面个数,这里指飞行器升降舵的个数;
——右舵、左舵相对弹体轴X1的滚转力矩系数对差动舵偏的导数;
ZR——升降舵质心距弹体中心对称面的距离;
——扰动前的等效俯仰舵偏角、差动舵偏角(rad);
Qγ(xR)——舵轴位置xR处的扭转振型;
Zh——伺服机构后支点距弹体中心对称面的距离;
Qγ(xh)——后支点xh处的扭转振型;
lRh——升降舵支架到伺服机构后支点的距离。
建模方法总结:
下面以俯仰通道为例概要给出对后缘控制舵滑翔飞行器弹性运动建模方法的研究成果。推导得到的弹体俯仰通道弹性运动方程见式(31):
推导得到的弹体俯仰通道弹性运动方程见式(31):
式中qi——第i次振型的广义坐标(m);
ξi——第i次振型的阻尼系数;
ωi——第i次振型的固有频率(rad/s);
D1i——飞行器俯仰(偏航)方向角速度对俯仰(偏航)第i阶振型的影响系数(m/s);
——飞行器(等效刚体纵轴)俯仰角速度的小扰动量(rad/s);
D2i——飞行器攻角(侧滑角)对俯仰(偏航)方向第i阶振型的影响系数(m/s2);
Δα——飞行器(等效刚体纵轴)攻角的小扰动量(rad);
D3i——升降舵(方向舵)对称舵偏角对俯仰(偏航)第i阶振型的影响系数(m/s2);
——升降舵俯仰通道对称舵偏角的小扰动量(rad);
——升降舵(方向舵)舵偏角加速度对俯仰(偏航)方向第i阶振型的影响系数(m);
D3i_ZK——姿控喷管操纵力对俯仰(偏航)方向第i阶振型的影响系数(m/s2);
ΔδZK——当俯仰(偏航)通道姿控喷管产生的控制力在飞行器体坐标系内为正时,ΔδZK取+1;反之取-1;俯仰(偏航)喷管不工作时取为0。
经推导,本文探讨的飞行器其俯仰方向运动方程式系数计算公式见式(32)至(36):
式中q——速度头(单位Pa),ρ:当地空气密度,V:飞行器空速(m/s);
SM——飞行器参考面积(m2);
Wi(x)——飞行器俯仰(偏航)方向第i阶x处的振型;
Mi——飞行器第i次振型的广义质量(kg),
——飞行器部段法向力系数对攻角的导数(/°);
xZ——飞行器质心距弹头顶点的距离(m);
——升降舵的个数,本文中讨论的飞行器
——单个升降舵(左舵、右舵)法向力系数对舵偏角的导数(/°);
——单个升降舵(左舵、右舵)绕舵轴的气动力矩系数对舵偏角的导数(/°);
xR——升降舵舵轴到弹头顶点的距离(m);
lK——飞行器参考长度(m);
mR——单个升降舵质量(kg);
lR——升降舵质心到舵轴的距离(m);
g——当地重力加速度(m/s2);
nx——飞行器轴向过载系数;
ny——飞行器法向过载系数;
——扰动前的升降舵等效俯仰舵偏角、差动舵偏角(rad);
xh——伺服机构后支点到弹头顶点的距离(m);
JR——单个升降舵绕舵轴的转动惯量(kg·m2);
xZK——姿控喷管到弹头顶点的距离(m);
fZK——俯仰(偏航)通道姿控喷管操纵力(N)。
Claims (4)
1.一种后缘舵滑翔飞行器的弹性运动建模方法,其特征在于:包括如下步骤:(1)利用模态正交性把弹体横向自由振动转化为各个相互独立的主振动的叠加;(2)根据线性小扰动假设,在外力作用下,弹体的横向振动用模态叠加来描述,振型函数由弹体刚度和质量分布、弹体的边界条件确定;(3)根据达朗伯原理建立振动微分方程;(4)分析影响滑翔飞行器弹性振动的外力:包括气动力,舵面控制力,姿态喷管控制力以及舵面摆动的惯性力;
所述的步骤(1)具体为:利用模态正交性把弹体横向自由振动转化为各个相互独立的主振动的叠加,即设
Wi(X)即弹体弹性的第i次固有振型函数,Φ=[W1(X),W2(X),…,Wn(X)],qi(t)为第i次固有振型对应的广义坐标;
所述的步骤(3)具体为:根据达朗伯原理建立振动微分方程为:
把式(1)代入式(2),并前乘ΦT得
式中,Mp=ΦTMΦ,Cp=ΦTCΦ,Kp=ΦTKΦ分别为广义模态质量矩阵、模态阻尼矩阵和模态刚度矩阵,都为对角矩阵,Fp=ΦTf为广义力阵,
式(3)写成
其中ωi——第i次振型的固有频率;
ξi——第i次振型的阻尼系数;
Fi——对应第i次振型的广义力;
Mi——对应第i次振型的广义质量;
且有
式中是沿弹体纵轴受到的外力在弹体坐标系轴Y1上的投影;m(X)是沿弹体纵轴的弹体质量分布;lK为飞行器纵轴总长度;ωi、ξi和Wi(X)是系统的固有属性,由数值方法求得并通过模态试验校验;
所述的步骤(4)具体为:对于俯仰通道,在弹体变形情况下,气动力沿纵轴各点局部攻角不同,气动力须看作沿纵轴的分布力,除弹体气动力外,其它外力都作用在弹体的某一位置上,是集中力;
影响弹体分布气动力的飞行姿态扰动包括攻角扰动Δα、俯仰角速度扰动舵偏角扰动角速度扰动及角加速度扰动将引起舵面控制力、摆动惯性力变化;姿控动力发动机的启用与关闭也产生控制力扰动;
用D1i和D2i分别表征与俯仰角速度扰动和攻角小扰动Δα成比例的广义气动力对第i次振型振动的影响;根据广义力定义和气动力计算方法推导出D1i和D2i的表达式(7)和(8):
式中q——飞行动压;V——飞行速度;
SM——飞行器参考面积;
lK——飞行器纵轴总长度;
——沿弹体纵轴分布的局部法向升力系数(单位:/°);
xZ——飞行器质心距弹头顶点的距离;
Wi(xn)——xn处对应的第i次振型。
2.根据权利要求1所述的后缘舵滑翔飞行器的弹性运动建模方法,其特征在于:后缘升降舵局部坐标系ORXRYRZR定义为:OR为舵轴中心线与控制舵交点处,ORXR平行弹体纵轴指向前,ORYR平行弹体坐标系Y1轴指向上,ORZR沿舵轴中心线,方向由右手法则确定;假设该升降舵偏角为对单个舵受力进行分析:
①舵面气动力和力矩
对小扰动进行线性化,假设扰动前的舵偏角为即
则舵偏小扰动产生的气动力在Y1轴方向的投影为
相应的气动力矩变化量为
其中——单个升降舵法向力系数导数(/°);
——单个升降舵绕ZR轴的气动力矩系数导数(/°),所有气动系数均沿局部坐标系各轴正方向定义;
②舵面惯性力和力矩
舵面质心处的惯性力与弹道轴向过载nx、法向过载ny以及舵偏角速度角加速度有关,将惯性力分解到ORYR轴和ORXR轴两个方向上;
忽略二阶小量,得舵偏小扰动产生的惯性力在Y1轴方向的投影为
相应的惯性力矩变化量
式中mR——单个升降舵质量;
lR——升降舵质心到舵轴的距离;
JR——单个升降舵相对于舵轴的转动惯量;
③伺服机构与舵系统的支反力
当飞行器升降舵采用直线式机电伺服机构驱动时,升降舵轴处受到的扭矩是通过摇臂与伺服机构作用力产生的力矩来平衡,此时舵面受到的外力矩对控制舱后舱体产生的作用力Fh在Y1轴方向的投影为
其中,lRh——升降舵支架到伺服机构后支点的距离;
则舵偏小扰动产生的后支点作用力在Y1轴方向的投影为
升降舵支架处舱体受到的作用力在Y1方向的投影与fhY大小相等,方向相反;④与舵偏有关的方程式系数
式(10)、(12)和(15)给出了舵偏小扰动引起的作用于弹体上的外力在Y1轴上的投影,这些外力均是集中力,忽略差动舵偏角影响时,代入公式(5)得广义力为
式中——升降舵的个数;
若考虑初始差动舵偏角影响,即右舵和左舵的舵偏角δ1和δ2为
考虑初始差动舵偏影响,推导得相应的方程式系数为
假设俯仰通道姿控发动机推力为fRCS,即沿Y1轴方向的投影,则相应的有
当俯仰通道姿控喷管产生的控制力在飞行器体坐标系内为正时,ΔδZK取+1,反之ΔδZKΔδRCS取-1,俯仰喷管不工作时ΔδZKΔδRCS取0;
综合上述与小扰动有关的广义力项,则俯仰通道的弹性振动方程为
3.根据权利要求2所述的后缘舵滑翔飞行器的弹性运动建模方法,其特征在于:偏航通道飞行器的控制是通过方向舵和姿控喷管,其偏航通道扰动运动方程形式同公式(22),其中的飞行姿态扰动换成了偏航角速度扰动侧滑角扰动Δβ、方向舵偏角扰动Δδψ和角加速度扰动
4.根据权利要求2所述的后缘舵滑翔飞行器的弹性运动建模方法,其特征在于:对于滚转通道弹性运动建模方法,其飞行姿态扰动换成了差动舵偏角扰动Δδγ和角加速度扰动舵面差动小扰动产生的绕弹体系X1轴的气动滚转力矩为
设左舵和右舵舵面惯性力的法向分量变化量为f2IY和f1IY,正方向均定义为沿Y1轴正方向,则产生滚转力矩为:
伺服机构后支点处产生的滚转力矩为
前支点处产生的滚转力矩大小相等,方向相反;
则滚转通道弹性振动方程中与舵面操纵有关的广义力为
与姿控喷管有关的方程式系数为
其中MRCS——滚转通道的姿控喷管控制力矩;
Mγ——扭转振型对应的广义质量;
Qγ(xRCS)——姿控喷管xRCS处的扭转振型;
最终推导得到飞行器滚转通道弹性振动方程为
其中与舵面操纵机构有关的方程式系数:
式中nγ——滚转通道控制舵面个数,这里指飞行器升降舵的个数;
——右舵、左舵相对弹体轴X1的滚转力矩系数对差动舵偏的导数;
ZR——升降舵质心距弹体中心对称面的距离;
——扰动前的等效俯仰舵偏角、差动舵偏角(rad);
Qγ(xR)——舵轴位置xR处的扭转振型;
Zh——伺服机构后支点距弹体中心对称面的距离;
Qγ(xh)——后支点xh处的扭转振型;
lRh——升降舵支架到伺服机构后支点的距离。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410602382.5A CN105629725B (zh) | 2014-10-31 | 2014-10-31 | 一种后缘舵滑翔飞行器的弹性运动建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410602382.5A CN105629725B (zh) | 2014-10-31 | 2014-10-31 | 一种后缘舵滑翔飞行器的弹性运动建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105629725A CN105629725A (zh) | 2016-06-01 |
CN105629725B true CN105629725B (zh) | 2019-01-08 |
Family
ID=56044795
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410602382.5A Active CN105629725B (zh) | 2014-10-31 | 2014-10-31 | 一种后缘舵滑翔飞行器的弹性运动建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105629725B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106228014B (zh) * | 2016-07-27 | 2019-01-29 | 江西洪都航空工业集团有限责任公司 | 一种导弹气动系数的获取方法 |
CN106950982B (zh) * | 2017-02-16 | 2020-02-14 | 北京临近空间飞行器系统工程研究所 | 再入飞行器姿控动力系统高空力矩特性辨识方法 |
CN107609307B (zh) * | 2017-10-10 | 2019-09-10 | 北京理工大学 | 一种考虑气弹和地球影响的远程火箭弹弹道分析方法 |
CN108256264B (zh) * | 2018-02-08 | 2020-03-31 | 北京航空航天大学 | 一种基于地面频响试验的气动伺服弹性稳定性预测方法 |
CN108958038B (zh) * | 2018-08-16 | 2021-04-23 | 上海航天控制技术研究所 | 一种适应飞行器推力离散特性的控制参数调整方法 |
CN109085847B (zh) * | 2018-09-03 | 2021-12-07 | 中国运载火箭技术研究院 | 一种大静不稳定度飞行器弹性稳定控制系统及控制方法 |
CN109656149B (zh) * | 2018-12-10 | 2021-07-30 | 上海卫星装备研究所 | 星箭耦合多体系统动力学计算试验方法及系统 |
CN110287505B (zh) * | 2019-03-20 | 2020-12-25 | 北京机电工程研究所 | 飞行器稳定性分析方法 |
CN112046761B (zh) * | 2020-08-04 | 2021-10-15 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种基于统计检验和滤波的飞机结冰在线探测方法 |
CN113505434B (zh) * | 2021-06-24 | 2022-10-28 | 上海机电工程研究所 | 基于气动力数学模型的飞行器设计制造方法及其飞行器 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101236573A (zh) * | 2007-05-25 | 2008-08-06 | 清华大学 | 一种柔性翼微型飞行器的流固耦合数值仿真方法 |
US8442705B2 (en) * | 2009-02-27 | 2013-05-14 | Airbus Operations Gmbh | Method and device for determining aerodynamic characteristics of an aircraft |
CN103970957A (zh) * | 2014-05-20 | 2014-08-06 | 北京航空航天大学 | 一种弹性乘波体高超声速飞行器仿真方法 |
-
2014
- 2014-10-31 CN CN201410602382.5A patent/CN105629725B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101236573A (zh) * | 2007-05-25 | 2008-08-06 | 清华大学 | 一种柔性翼微型飞行器的流固耦合数值仿真方法 |
US8442705B2 (en) * | 2009-02-27 | 2013-05-14 | Airbus Operations Gmbh | Method and device for determining aerodynamic characteristics of an aircraft |
CN103970957A (zh) * | 2014-05-20 | 2014-08-06 | 北京航空航天大学 | 一种弹性乘波体高超声速飞行器仿真方法 |
Non-Patent Citations (4)
Title |
---|
Elastic Multibody Models of Transport Aircraft High-Lift Mechanisms;Janos Zierath 等;《JOURNAL OF AIRCRAFT》;20091031;第40卷(第5期);全文 |
乘波布局高超声速巡航飞行器轨迹控制问题研究;郝创洲;《中国优秀硕士学位论文全文数据库工程科技II辑》;20140615(第06期);全文 |
弹性高超声速飞行器建模及精细姿态控制;孟中杰 等;《宇航学报》;20110830;第32卷(第8期);第1684-1685页 |
高超声速巡航飞行器振动建模及精细姿态控制;孟中杰 等;《西北工业大学学报》;20110630;第29卷(第3期);全文 |
Also Published As
Publication number | Publication date |
---|---|
CN105629725A (zh) | 2016-06-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105629725B (zh) | 一种后缘舵滑翔飞行器的弹性运动建模方法 | |
Li et al. | Dynamics modeling and simulation of flexible airships | |
CN102073755A (zh) | 近空间高超声速飞行器运动控制仿真方法 | |
US11104420B2 (en) | System and method for minimising buffeting | |
CN109703768B (zh) | 一种基于姿态/轨迹复合控制的软式空中加油对接方法 | |
CN109703769B (zh) | 一种基于预瞄策略的空中加油对接控制方法 | |
Zhang et al. | Analysis of dynamic characteristics of the rigid body/elastic body coupling of air-breathing hypersonic vehicles | |
Oh et al. | Autonomous helicopter landing on a moving platform using a tether | |
CN113919081B (zh) | 一种考虑惯性耦合的柔性飞行动力学建模与分析方法 | |
Blunt et al. | Trailing edge tabs on folding wingtips (fwts) for aircraft roll control | |
Grauer et al. | Aeroelastic modeling of X-56A stiff-wing configuration flight test data | |
Liu et al. | Continuous dynamic simulation for morphing wing aeroelasticity | |
Tong et al. | Dynamic response analysis under atmospheric disturbances for helicopters based on elastic blades | |
Banazadeh et al. | Control effectiveness investigation of a ducted-fan aerial vehicle using model predictive controller | |
Noll et al. | Active flexible wing program | |
Tuzcu et al. | Aeroelastic modeling and adaptive control of GTM | |
Tang et al. | Aeroelastic scaling laws for gust load alleviation control system | |
Ricci et al. | Active control of three-surface aeroelastic model | |
Li et al. | A comparison of the fixed-axes and the mean-axes modeling methods for flexible aircraft simulation | |
Qiu et al. | Improving aileron effectiveness based on changing the position of aileron connectors | |
Xiaofeng et al. | Overload control system design of a deflectable nose missile | |
Peyada et al. | Mathematical modelling, simulation, and estimation of aircraft parameters using five degree-of-freedom dynamic test rig | |
Du et al. | Flight Control Law Design and Analysis for Mars Airplane in Crosswind | |
Sharma | Development and Application of a Comprehensive Simulation for Modeling Helicopter Ship Landing | |
Kuzmina et al. | Analytical-experimental study on using different control surfaces to alleviate dynamic loads |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |