CN104391364B - 基于光缆敷设的共轴双旋翼无人直升机参数辨识试验*** - Google Patents
基于光缆敷设的共轴双旋翼无人直升机参数辨识试验*** Download PDFInfo
- Publication number
- CN104391364B CN104391364B CN201410693197.1A CN201410693197A CN104391364B CN 104391364 B CN104391364 B CN 104391364B CN 201410693197 A CN201410693197 A CN 201410693197A CN 104391364 B CN104391364 B CN 104391364B
- Authority
- CN
- China
- Prior art keywords
- optical cable
- theta
- cable
- cos
- platform
- 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
- 239000000835 fiber Substances 0.000 title claims abstract description 8
- 230000005540 biological transmission Effects 0.000 claims abstract description 11
- 230000001133 acceleration Effects 0.000 claims abstract description 10
- 238000005259 measurement Methods 0.000 claims abstract description 9
- 230000003287 optical effect Effects 0.000 claims description 138
- 238000000034 method Methods 0.000 claims description 37
- 238000013095 identification testing Methods 0.000 claims description 25
- 230000007246 mechanism Effects 0.000 claims description 14
- 239000013598 vector Substances 0.000 claims description 14
- 238000012937 correction Methods 0.000 claims description 12
- 230000005284 excitation Effects 0.000 claims description 10
- 230000000737 periodic effect Effects 0.000 claims description 10
- 238000004458 analytical method Methods 0.000 claims description 7
- 230000004044 response Effects 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 5
- 238000009795 derivation Methods 0.000 claims description 4
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 230000036544 posture Effects 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000006073 displacement reaction Methods 0.000 claims description 2
- 230000009977 dual effect Effects 0.000 claims description 2
- 230000005484 gravity Effects 0.000 claims description 2
- 238000004804 winding Methods 0.000 claims description 2
- 238000009434 installation Methods 0.000 abstract 2
- 238000013461 design Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 6
- 238000012360 testing method Methods 0.000 description 6
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 description 4
- 238000007781 pre-processing Methods 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 125000004122 cyclic group Chemical group 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000005096 rolling process Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000000087 stabilizing effect Effects 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000004148 unit process Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B6/00—Light guides; Structural details of arrangements comprising light guides and other optical elements, e.g. couplings
- G02B6/46—Processes or apparatus adapted for installing or repairing optical fibres or optical cables
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C23/00—Combined instruments indicating more than one navigational value, e.g. for aircraft; Combined measuring devices for measuring two or more variables of movement, e.g. distance, speed or acceleration
- G01C23/005—Flight directors
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Optics & Photonics (AREA)
- Laying Of Electric Cables Or Lines Outside (AREA)
- Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)
Abstract
基于光缆敷设的共轴双旋翼无人直升机参数辨识试验***,包括飞行器平台和地面设备两个部分,所述飞行器平台上分别加装惯性测量单元、全球定位***、磁航向计、激光高度表和数据记录仪和控制器,所述惯性测量单元与所述全球定位***、所述磁航向计、所述激光高度表、所述数据记录仪均连接至所述控制器,并由控制器将信号处理,并输出至控制单元,所述地面设备包括无线电遥控发射机、数传终端和地面控制计算机***,该地面控制计算机***与所述RC遥控发射机连接,该地面控制计算机***还与数传终端连接,所述飞行器平台上集成有传感器,该传感器将测量飞行器平台的三轴角速度、线加速度、航向角、三轴位置、三轴线速度。
Description
技术领域
本发明涉及一种用于共轴双旋翼无人直升机的敷设柔性光缆的参数辨识试验***,属于航空技术领域。
背景技术
无人直升机由于其具备起降便捷、定点悬停、机动性强的特点,因此作战用途广泛,具有非常可观的作战潜力和发展前景。国外开展无人直升机研制比较早的国家有美国、英国、加拿大、德国等。其中美国上世纪50年代就开始研究无人直升机,60年代QH-50A反潜遥控直升机试飞成功,在越南战争中美国陆军就使用QH-50D无人直升机执行侦察、战场监视、炮兵目标观测等任务。由于无人直升机飞行控制技术存在诸多难点,所以无人直升机相对于固定翼无人机而言,技术和应用的发展一直较为迟缓。
在无人直升机中,共轴双旋翼是一种重要的直升机布局形式,共轴双旋翼无人直升机和常规无人机布局有很大差别,共轴双旋翼无人直升机具有绕同一轴线正反旋转的上下两幅旋翼,由于两幅旋翼旋转方向相反,所以其产生的扭矩在航向恒定的情况下相互抵消,因而取消了尾桨;同时通过上下旋翼总距的差动操作来实现无人直升机的航向操纵,两幅旋翼在直升机的飞行过程中既是升力面又是纵横向和航向的操纵面。共轴双旋翼无人直升机由于独特的布局形式,非常适合光缆一类柔性线缆的快速敷设,但由于光缆快速敷设涉及柔性电缆及其放线机构,所以***的动力学建模必须综合考虑拖拽式的柔性光缆对飞行器动力学特性的影响,因此其建模方法比较复杂。
发明内容
本发明的目的是针对用于光缆快速敷设的无人直升机飞行器平台、光缆释放机构和光缆所组成的综合体进行统一建模,从而获得无人直升机在进行快速光缆敷设时的动力学模型,为该无人直升机的飞行控制***研制提供支持。
本发明包括以下两部分内容
本发明的第一方面提供将恒张力释放状态下的柔性光缆转换为多刚体动力学建模的具体方法。
本发明第二方面提供一种可以面向飞行控制***设计实际应用的光缆敷设无人直升机动力学模型,并设计一个安全的参数辨识试验***。
本发明的第一方面提供的恒张力释放状态下的柔性光缆动力学建模方法如下:
1恒张力释放状态下的柔性光缆动力学建模方法
由于光缆的一端固定于直升机平台之上,另一端为自由端,光缆在恒张力释放装置控制下在保证光缆张力恒定的情况下释放。因此,我们在进行光缆受力分析时,可以认为光缆的一端为固定端,该固定端静止或匀速运动。同时为了使用多刚体***动力学方法描述光缆运动,设定:
a)整条待布设的光缆由N个光缆段构成;
b)各段逐次从释放机构中释放,正在释放的光缆段的速度与刚释放的光缆段的速度一致;
c)各光缆段均视作刚性杆,其质量集中于各光缆段的一个端点;
d)各个刚性杆之间通过绞进行连接。
单个光缆段的受力计算方法
考虑到直升机用于光缆敷设的工程需求,第j段绳子的阻力系数CDj和升力系数CLj由工程经验公式(3.1)和(3.2)进行计算,其中αj为第j段绳子的攻角,攻角也称迎角,为翼弦与来流速度之间的夹角,翼弦抬头为正,翼弦低头为负,常用符号α表示:
CDj≈0.022+1.1sin3αj (3.1)
CLj≈1.1sin2αj cosαj (3.2)
光缆段相对于风的速度为Vwind,第J个光缆段质心的速度为VJ,则根据力的平衡和力的分解原理,第J个光缆段相对于风的速度的大小为:
其中,---第J个光缆段相对于风的速度;
VJ---第J个光缆段质心的速度;
VJ-1---第J-1个光缆段质心的速度;
Vwind---风速;
---第J个光缆段相对于风的速度在水平方向的分量;
---第J个光缆段相对于风的速度在竖直方向的分量;
i---水平方向的方向向量;
j---竖直方向的方向向量;
设第J段光缆的气动阻力系数为CDJ、气动升力系数为CLJ,每个光缆段上的气动力主要包括气动阻力和气动升力,根据受力分析并结合式(3.1)、(3.2)得,第J段光缆的气动阻力和气动升力分别为:
式(3.4)、(3.5)中
---第J光缆段的气动阻力;
---第J光缆段的气动升力;
ρ---空气密度;
CD---光缆段的气动阻力系数;
CL---光缆段的气动升力系数;
d---光缆段直径;
---第J光缆段相对于风的速度;
eD---第J光缆段的气动阻力方向单位矢量;
eL---第J光缆段的气动升力方向单位矢量。
根据以上的推导可得出第J段光缆质点处的气动力的值为:
其中:---第J光缆段质点处的气动力;
---第J光缆段的气动阻力;
---第J光缆段的气动升力;
---第J+1光缆段的气动阻力;
---第J+1光缆段的气动升力;
均为矢量。
1.2整体光缆动力学模型建立方法
根据前述设定,整条光缆看作是由多个刚性杆通过绞进行连接的结构体,针对无人直升机布设光缆的实际状态,设定如下建模条件:
1.2.1由于光缆的第一段盘绕在恒张力铺设机构上,所以设定该第一段光缆段为固定端且固定于直升机平台之上,光缆最后一段为贯序着陆的自由端;
1.2.2自由端所受到的在x轴、y轴上力分量通过测量获得;
1.2.3直升机在正常布设光缆的过程中,其飞行高度保持固定,即在光缆全部释放之前,直升机平台所拖曳的停留在空中的光缆段的长度固定;
将各个光缆段间的运动约束关系转化为各个光缆段张力的递推关系,最后得出光缆段约束力Tj与***外力的线性关系。
对于结点Mj-1,其质量为mj-1,结点Mj,其质量为mj,设这两个结点所受的气动力分别为和则各个光缆段结点的运动方程为
其中,mj-1为j-1段光缆的质量,a(j-1)x为j-1段光缆x方向的加速度,a(j-1)y为j-1段光缆y方向的加速度,Tj-1为j-1段光缆所受的外力,Tj为j段光缆所受的外力,θj为j段光缆与飞行器平台z轴的夹角;
和分别是j-1所受的空气动力在x和y方向的分量;g为重力常数;式(3.9)描述最后一段的受力关系,其中Flastx、Flasty是自由端所受到的力在x和y方向的分量。
各个光缆段的加速度约束关系为:
其中:ajx为j段光缆在x方向的加速度,ajy为j段光缆在y方向的加速度,lj为j段光缆的长度,θj为j段光缆与飞行器平台z轴的夹角,ωj为j段光缆与飞行器平台z轴夹角的角速度。
将式(3.7)带入式(3.10)得到式(3.12),式(3.8)和(3.9)分别代入(3.11)式中,各个式子中的加速度项被消掉,可以获得各段光缆之间的铰约束力递推关系:
解上述线性方程组(3.12)(3.13)(3.14),经过化简后得到各光缆段结点的张力约束关系如下的线性方程:
Tj+1=Aj+1Tj+Bj+1Tj-1+Cj+1 (j=2,3,...n-1) (3.15)
Tj代表结点Mj所受的张力。
根据(3.12)式,可以列出如下方程组:
其中:Aj+1、Bj+1、Cj+1,均为系数,其取值与光缆段的长度、质量和姿态有关。从上述约束关系的推导过程,可以得出三个系数仅与光缆段的位移和速度有关。这说明在定高飞行情况下,当光缆释放速度、无人机平台速度以及风速恒定时,整个光缆及光缆段所受到的张力也恒定,这也是我们通过控制光缆释放速度来控制张力的原理。这个结果也表明,恒张力控制不仅仅是为了满足光缆强度要求Tj小于光缆规定的张力,同时保持光缆的张力使其具备规定的空中姿态,确保不与旋翼发生卷绕或干涉,此外也是为了降低光缆对飞行器平台的力与力矩的扰动,有利于***的操稳特性。
本发明第二方面提供的面向飞行控制***的光缆敷设无人直升机动力学模型及参数辨识试验***如下:
2.1面向控制的光缆敷设共轴双旋翼无人直升机平台动力学模型
直升机运动可以表示成如下矢量方程:
这里x是直升机的状态矢量,即数学意义上的x的导数,u是控制输入变量,这里x=[u,v,w,φ,θ,ψ,p,q,r],其中各参数代表的含义分别为:滚转角速率p,俯仰角速率q,航向角速率r,滚转角φ,俯仰角θ,航向角ψ,前向速度u,侧向速度v,垂直速度w。u是直升机的控制输入变量,[ulat,ulon,ucol,uyaw],其中ulat横向周期变矩,ulon纵向周期变矩,ucol主桨距总矩,uyaw上下旋翼差动总距。
直升机线运动方程为:
其中m为直升机质量,∑X、∑Y、∑Z为作用在直升机机体轴x,y,z方向上的合力,这里我们将光缆整体的张力在各个坐标轴上的分量对直升机平台运动的影响明确列出,并表示为Tx,Ty,Tz,
式(3.19)中,Ixx是直升机对OX轴的转动惯量,Iyy是直升机对OY轴的转动惯量,Izz是直升机对OZ轴的转动惯量。Ixz是直升机对OX轴的惯性积。∑L、∑M、∑N分别为绕三轴的力矩和,同样的,LTX、MTY、NTZ分别表示光缆组件绕无人直升机平台三轴的力矩。
对于用于光缆敷设的无人机平台来说,其动力学特征主要包括光缆的动力学和共轴双旋翼直升机各组成部件的动力学特性共同决定,附图显示共轴双旋翼无人直升机动力学建模中,光缆以及各个飞行器平台的部件气动特性与飞行器平台总体运动特性之间的关系。
不考虑柔性光缆,直升机平台本身的动力学模型非常复杂,仅采用上述机理建模方式,较为准确反应直升机悬停状态的的线性化模型就需要超过70个参数,至于非线性模型的参数就更加复杂了。一般而言,这些模型可以用于模拟和仿真,直接用于控制***设计就十分困难,以至于不能实现,所以保证模型参数准确性的同时,简化模型的复杂度对控制***的设计至关重要。
我们这里将光缆敷设过程中光缆对于***的扰动和共轴双旋翼无人机一并考虑,提出了一个实用的、面向线性控制***设计的一体化模型。
为了简化,仅反映直升机的水平面姿态和速度特性的动力学特性,是一个2输入、6输出的***,式(3.14)线性化后,将直升机的状态空间方程表达如下:
其中:b-横向旋翼挥舞角;a-纵向旋翼挥舞角;
这里处理光缆影响的方法是将其作为直升机平台的一个外在扰动,通过对平台模型中相应气动导数增加修正量的方式,来获得与现实一致的模型响应。气动导数及其修正量可以通过下面的参数辨识的方法统一获得。
根据上述方法,确定的***模型结构为:
其中,Ba,Ab,La,Lb,Lu,Lv,Ma,Mb,Mu,Mv,Xu,Yv,Xa,Yb均为气动导数,这些也就是我们需要确定的参数,特别地,α/τ,为光缆组件对旋翼挥舞角影响的修正量,g/k为光缆组件对飞行平台姿态角影响的修正量。
2.2基于光缆敷设的共轴双旋翼无人直升机参数辨识试验***
为了高效率实施拖曳柔性线缆的共轴双旋翼无人机上述气动导数的辨识,我们构建一个具体的辨识试验***。
该试验***包括飞行器平台和地面设备两个部分:所述飞行器平台上分别加装惯性测量单元(Inertial Measuring Unit,简称IMU)、全球定位***(Global PositioningSystem,简称GPS)、磁航向计、激光高度表和数据记录仪、控制器,所述惯性测量单元IMU与、所述GPS、所述磁航向计、所述激光高度表、所述数据记录仪均连接至所述控制器,并由控制器将信号处理,并输出至控制单元,从而控制飞行器平台的动作。由于我们需要在飞行器平台多个工作点进行参数辨识试验,当飞行器平台处于大速度(大速度是指大于20m/s的速度)前飞时,无人直升机很快会超出操作手视线,威胁飞行安全,因此本***中设置两种工作模式,一是透明模式,即直接由操作手控制执行机构,二是自动模式,即利用直升机内置稳定算法进行闭环控制共轴双旋翼的姿态。
所述地面设备包括无线电遥控发射机(简称RC遥控发射机)、数传终端和地面控制计算机***,该地面控制计算机***与所述RC遥控发射机连接,该地面控制计算机***还与数传终端连接。
所述飞行器平台上集成有传感器,该传感器将测量飞行器平台的三轴角速度、线加速度、航向角、三轴位置、三轴线速度,以及相对高度及高度变化率,经过卡尔曼滤波之后将形成飞行器平台的位置、三轴角速率和欧拉角数据。这些飞行器平台参数、遥控接收机的控制给定、执行机构的行程都将被控制器打上时标,统一存储于数据记录仪中。同时将该数据下传至地面计算机***。地面计算机采集航模操纵手的输入,并记录数据终端传回的飞行器平台控制及飞行参数。
在起始状态下,航模操纵手直接控制飞行器平台,航模操纵手施加给飞行器平台的激励为横向周期变距、纵向周期变距,该横向周期变距和纵向周期变距直接传递给飞行器平台的控制器。同时,这两个信号首先被地面计算机采集,地面计算机通过数据终端获得了飞行器平台的飞行参数。当飞行器平台的姿态角大于给定的预警值时,此时飞行器处于比较危险的状态,控制器安全切换逻辑单元将断开接收机所获得的航模操纵手的激励信号,直接将飞行器平台处于控制器的闭环控制之下,保证飞行器平台的安全。
2.3基于光缆敷设的共轴双旋翼无人直升机参数辨识试验***的参数辨识步骤
2.3.1在悬停和前飞状态下,通过遥控器在各个通道输入给定激励信号并采集飞行器飞行状态变量:
操纵信号为阶跃信号,分别输入纵向通道和横向通道,由机载记录仪记录直升机的三个欧拉角和三个角速率,即姿态模型中的状态变量,其中,欧拉角的单位为度,角速率的单位为度/秒。机载***还对操纵输入信号进行采集,操纵输入信号表示为舵面的角度,单位也是度。
机载***以50HZ的采集速度对悬停状态下对飞行手的操纵输入和对应的状态数据进行采集,飞行试验共采集2000组飞行数据。
飞行器在自动模式下将前飞速度控制在10m/s,然后选择透明模式,机载***以50Hz的采集速率对前飞状态下对飞行手的操纵输入和对应的状态数据进行采集,飞行试验共采集2000组飞行数据。
2.3.2数据预处理
对上述数据进行预处理,无人直升机的各个通道之间总会存在一定程度的耦合,具体地,横向周期变距是影响横滚角速度的主要因素,但是其他通道的输入对横滚角速度也有影响。因此,进行辨识之前用偏相干分析法对数据进行进一步处理,除去其他通道对横滚通道的影响。
2.3.3预测误差法进行姿态模型参数辨识
采用MATLAB包含的预报误差方法(简称PEM)进行姿态模型辨识,预报误差方法是一种以极小化模型输出和实际输出的二次误差为目标的线性***辨识方法。
更多操作对于本领域普通技术人员而言都是已知的,不再赘述。
本发明所提供的基于光缆敷设的共轴双旋翼无人直升机参数辨识试验***的技术方案包括上述各部分的任意组合,上述各部分组件的简单变化或组合仍为本发明的保护范围。
本发明所提供的光缆敷设的共轴双旋翼无人直升机参数辨识试验***为用于光缆或其他线缆布设的共轴双旋翼无人直升机的控制提供理论方法基础和工程技术指导。
本发明提供的基于光缆敷设的共轴双旋翼无人直升机参数辨识试验***安全高效,极大地提高参数辨识精度,具有积极意义。
共轴双旋翼无人直升机的参数辨识试验***为无人机的使用提供了精确的控制方式,提高柔性线缆敷设精度。
附图说明
图1为按照本发明的共轴双旋翼无人直升机平台动力学模型的一优选实施例的建模框架示意图。
图2为按照本发明第二方面的共轴双旋翼无人直升机的参数辨识试验***的一优选实施例的示意图。
图3为按照本发明第一方面的共轴双旋翼无人直升机拖拽柔性线缆的建模方法的一优选实施例的整体受力的简化模型示意图。
图4为图3所示实施例的柔性线缆释放过程中的受力分析示意图。
图5为图3所示实施例的绳段单元所受的气动力示意图。
图6为图3所示实施例在不同风速下线缆张力计算结果示意图。
图7为图3所示实施例相同横向变距输入下的实飞和仿真滚转响应。
图8为图3所示实施例相同纵向变距输入下的实飞和仿真纵向响应。
图9为图3所示实施例相同航向变距下的实飞和仿真航向响应。
图10为图3所示实施例不同风速下模型准确度分析示意图。
具体实施方式
为了更好地理解本发明,下面结合附图分别详细描述按照本发明的基于光缆敷设的共轴双旋翼无人直升机参数辨识试验***的优选实施例。
实施例1:如图3-5所示,先建立面向控制的光缆敷设共轴双旋翼无人直升机平台动力学模型,过程如下:
直升机运动可以表示成如下矢量方程:
这里x是直升机的状态矢量,即数学意义上的x的导数,u是直升机的控制输入变量,这里x=[u,v,w,φ,θ,ψ,p,q,r],其中各参数代表的含义分别为:滚转角速率p,俯仰角速率q,航向角速率r,滚转角φ,俯仰角θ,航向角ψ,前向速度u,侧向速度v,垂直速度w。u是直升机的控制输入变量,[ulat,ulon,ucol,uyaw],其中符合的意思分别为:ulat横向周期变矩,ulon纵向周期变矩,ucol主桨距总矩,uyaw上下旋翼差动总距。
直升机线运动方程为:
其中m为直升机质量,∑X、∑Y、∑Z为作用在直升机机体轴x,y,z方向上的合力,这里将光缆整体的张力在各个坐标轴上的分量对直升机平台运动的影响明确列出,并表示为Tx,Ty,Tz,。
式(3.19)中,Ixx是直升机对OX轴的转动惯量,Iyy是直升机对OY轴的转动惯量,Izz是直升机对OZ轴的转动惯量。Ixz是直升机对OX轴的惯性积。∑L、∑M、∑N分别为绕三轴的力矩和,同样的,LTX、MTY、NTZ分别表示光缆组件绕无人直升机平台三轴的力矩。
对于用于光缆敷设的无人机平台来说,其动力学特征主要包括光缆的动力学和共轴双旋翼直升机各组成部件的动力学特性共同决定,附图显示共轴双旋翼无人直升机动力学建模中,光缆以及各个飞行器平台的部件气动特性与飞行器平台总体运动特性之间的关系。
不考虑柔性光缆,直升机平台本身的动力学模型非常复杂,仅采用上述机理建模方式,较为准确反应直升机悬停状态的的线性化模型就需要超过70个参数,至于非线性模型的参数就更加复杂了。一般而言,这些模型可以用于模拟和仿真,直接用于控制***设计就十分困难,以至于不能实现,所以保证模型参数准确性的同时,简化模型的复杂度对控制***的设计至关重要。
这里将光缆敷设过程中光缆对于***的扰动和共轴双旋翼无人机一并考虑,提出了一个实用的、面向线性控制***设计的一体化模型。
为了简化,仅反映直升机的水平面姿态和速度特性的动力学特性,是一个2输入、6输出的***,式(3.14)线性化后,将直升机的状态空间方程表达如下:
其中:b-横向旋翼挥舞角;a-纵向旋翼挥舞角;
这里处理光缆影响的方法是将其作为直升机平台的一个外在扰动,通过对平台模型中相应气动导数增加修正量的方式,来获得与现实一致的模型响应。气动导数及其修正量可以通过下面的参数辨识的方法统一获得。
根据上述方法,确定的***模型结构为:
其中,Ba,Ab,La,Lb,Lu,Lv,Ma,Mb,Mu,Mv,Xu,Yv,Xa,Yb均为气动导数-。这些也就是需要确定的参数。特别的,α/τ,为光缆组件对旋翼挥舞角影响的修正量,g/k为光缆组件对飞行平台姿态角影响的修正量。
如图2所示,建立基于光缆敷设的共轴双旋翼无人直升机的参数辨识试验***
为了高效率实施拖曳柔性线缆的共轴双旋翼无人机上述气动导数的辨识,构建一个具体的辨识试验***。该试验***包括飞行器平台和地面设备两个部分。
所述飞行器平台上分别加装惯性测量单元(Inertial Measuring Unit,简称IMU)、全球定位***(Global Positioning System,简称GPS)、磁航向计、激光高度表和数据记录仪、控制器,所述惯性测量单元IMU与、所述GPS、所述磁航向计、所述激光高度计分别连接至滤波器,所示滤波器又连接所述数据记录仪,该数据记录仪连接至所述控制器安全切换逻辑单元,该控制器安全切换逻辑单元将信号处理,并输出至执行机构,从而控制飞行器平台的动作,所述飞行器平台的数传链路与所述控制器安全切换逻辑单元连接,且他们之间的数据互相传递,所述控制器安全切换逻辑单元还连接接收机,所述控制器安全切换逻辑单元能与滤波器之间交换。
由于需要在飞行器平台多个工作点进行参数辨识试验,当飞行器平台处于大速度(大速度是指大于20m/s的速度)前飞时,无人直升机很快会超出操作手视线,威胁飞行安全,因此本***中设置两种工作模式,一是透明模式,即直接由操作手控制执行机构,二是自动模式,即利用直升机内置稳定算法进行闭环控制共轴双旋翼的姿态。
所述地面设备包括无线电遥控发射机(简称RC遥控发射机)、数传终端和地面控制计算机***,该地面控制计算机***与所述RC遥控发射机连接,该地面控制计算机***还与数传终端连接,该数传终端将地面操控信号通过数传链路传递至所述控制器安全切换逻辑单元。所述RC遥控发射机发射信号,并由飞行器平台的接收机接收,然后传递至所述控制器安全切换逻辑单元进而控制飞行器平台的动作。
所述飞行器平台上集成有传感器,该传感器将测量飞行器平台的三轴角速度、线加速度、航向角、三轴位置、三轴线速度,以及相对高度及高度变化率,经过卡尔曼滤波之后将形成飞行器平台的位置、三轴角速率和欧拉角数据。这些飞行器平台参数、接收机的控制给定、执行机构的行程都将被控制器打上时标,统一存储于数据记录仪中。同时将该数据下传至地面计算机***。地面计算机采集航模操纵手的输入,并记录数据终端传回的飞行器平台控制及飞行参数。
在起始状态下,航模操纵手直接控制飞行器平台,航模操纵手施加给飞行器平台的激励为横向周期变距、纵向周期变距,该横向周期变距和纵向周期变距直接传递给飞行器平台的控制器。同时,这两个信号首先被地面计算机采集,地面计算机通过数据终端获得了飞行器平台的飞行参数。当飞行器平台的姿态角大于给定的预警值时,此时飞行器处于比较危险的状态,控制器安全切换逻辑单元将断开接收机所获得的航模操纵手的激励信号,直接将飞行器平台处于控制器的闭环控制之下,保证飞行器平台的安全。
下面给出上述建模的一个实际结果。本模型的结构如3.21所示,我们采用上面所述之方法,建立了一个面向控制的光缆敷设共轴双旋翼无人直升机平台动力学模型:
为了验证建模方法的有效性和模型的准确度,对共轴双旋翼无人直升机动力学模型和实际***施加同样激励,比较了两者在相同激励下(横向周期变距、纵向周期变距以及上下旋翼半差动变距下)的响应,结果如图7-图9所示。
可以看出,两者具有很高的一致性,相似度均大于60%,这对于控制***的设计已经足够了。同时我们分析了风场扰动对***模型准确性的影响,如图10所示结果表明,风速越大,数学模型和实际***的形似度越低,当风速大于10m/s时,模型的纵向特性与实际***的相似度将小于60%。
基于光缆敷设的共轴双旋翼无人直升机参数辨识试验***的参数辨识步骤如下:2.3.1在悬停和前飞状态下,通过遥控器在各个通道输入给定激励信号并采集飞行器飞行状态变量:操纵信号为阶跃信号,分别输入纵向通道和横向通道,由机载记录仪记录直升机的三个欧拉角和三个角速率,即姿态模型中的状态变量,其中,欧拉角的单位为度,角速率的单位为度/秒。机载***还对操纵输入信号进行采集,操纵输入信号表示为舵面的角度,单位也是度。
机载***以50HZ的采集速度对悬停状态下对飞行手的操纵输入和对应的状态数据进行采集,飞行试验共采集2000组飞行数据。
飞行器在自动模式下将前飞速度控制在10m/s,然后选择透明模式,机载***以50Hz的采集速率对前飞状态下对飞行手的操纵输入和对应的状态数据进行采集,飞行试验共采集2000组飞行数据。
2.3.2数据预处理
对上述数据进行预处理,无人直升机的各个通道之间总会存在一定程度的耦合,具体地,横向周期变距是影响横滚角速度的主要因素,但是其他通道的输入对横滚角速度也有影响。因此,进行辨识之前用偏相干分析法对数据进行进一步处理,除去其他通道对横滚通道的影响。2.3.3预测误差法进行姿态模型参数辨识
采用MATLAB包含的预报误差方法(简称PEM)进行姿态模型辨识,预报误差方法是
一种以极小化模型输出和实际输出的二次误差为目标的线性***辨识方法。
更多操作对于本领域普通技术人员而言都是已知的,不再赘述。
本发明所提供的基于光缆敷设的共轴双旋翼无人直升机参数辨识试验***的技术方案包括上述各部分的任意组合,上述各部分组件的简单变化或组合仍为本发明的保护范围。
本发明提供的基于光缆敷设的共轴双旋翼无人直升机参数辨识试验***安全高效,极大地提高参数辨识精度,具有积极意义。
Claims (3)
1.基于光缆敷设的共轴双旋翼无人直升机参数辨识试验***,包括飞行器平台和地面设备两个部分,所述飞行器平台上分别加装惯性测量单元、全球定位***、磁航向计、激光高度表和数据记录仪和控制器,所述惯性测量单元与所述全球定位***、所述磁航向计、所述激光高度表、所述数据记录仪均连接至所述控制器,并由控制器将信号处理,并输出至控制单元,所述地面设备包括无线电遥控发射机、数传终端和地面控制计算机***,该地面控制计算机***与RC遥控发射机连接,该地面控制计算机***还与数传终端连接,所述飞行器平台上集成有传感器,该传感器将测量飞行器平台的三轴角速度、线加速度、航向角、三轴位置、三轴线速度,以及相对高度及高度变化率,经过卡尔曼滤波之后将形成飞行器平台的位置、三轴角速率和欧拉角数据;这些飞行器平台参数、遥控接收机的控制给定、执行机构的行程都将被控制器打上时标,统一存储于数据记录仪中,同时将该数据下传至地面计算机***,地面计算机采集航模操纵手的输入,并记录数据终端传回的飞行器平台控制及飞行参数,其特征在于:所述基于光缆敷设的共轴双旋翼无人直升机参数辨识试验***设置两种工作模式,以便不同场合选择不同的工作模式,一是透明模式,即直接由操作手控制执行机构,二是自动模式,即利用直升机内置稳定算法进行闭环控制共轴双旋翼的姿态;
在起始状态下,航模操纵手直接控制飞行器平台,航模操纵手施加给飞行器平台的激励为横向周期变距、纵向周期变距,该横向周期变距和纵向周期变距直接传递给飞行器平台的控制器;同时,这两个信号首先被地面计算机采集,地面计算机通过数据终端获得了飞行器平台的飞行参数;当飞行器平台的姿态角大于给定的预警值时,此时飞行器处于比较危险的状态,控制器安全切换逻辑单元将断开接收机所获得的航模操纵手的激励信号,直接将飞行器平台处于控制器的闭环控制之下,保证飞行器平台的安全;
所述基于光缆敷设的共轴双旋翼无人直升机参数辨识试验***采用的建模方法包括,单个光缆段的受力计算方法和整体光缆动力学模型建立方法;其中单个光缆段的受力计算方法包括第j段绳子的阻力系数CDj和升力系数CLj由工程经验公式(3.1)和(3.2)进行计算,其中αj为第j段绳子的攻角,翼弦抬头为正,翼弦低头为负,
CDj≈0.022+1.1sin3αj (3.1)
CLj≈1.1sin2αjcosαj (3.2)
光缆段相对于风的速度为Vwind,第J个光缆段质心的速度为VJ,则根据力的平衡和力的分解原理,第J个光缆段相对于风的速度的大小为:
其中,---第J个光缆段相对于风的速度;
VJ---第J个光缆段质心的速度;
VJ-1---第J-1个光缆段质心的速度;
Vwind---风速;
---第J个光缆段相对于风的速度在水平方向的分量;
---第J个光缆段相对于风的速度在竖直方向的分量;
i---水平方向的方向向量;
j---竖直方向的方向向量;
设第J段光缆的气动阻力系数为CDJ、气动升力系数为CLJ,每个光缆段上的气动力主要包括气动阻力和气动升力,根据受力分析并结合式(3.1)、(3.2)得,第J段光缆的气动阻力和气动升力分别为:
式(3.4)、(3.5)中
---第J光缆段的气动阻力;
---第J光缆段的气动升力;
ρ---空气密度;
CD---光缆段的气动阻力系数;
CL---光缆段的气动升力系数
d---光缆段直径;
---第J光缆段相对于风的速度;
eD---第J光缆段的气动阻力方向单位矢量;
eL---第J光缆段的气动升力方向单位矢量;
根据以上的推导可得出第J段光缆质点处的气动力的值为:
---第J光缆段质点处的气动力;
---第J光缆段的气动阻力;
---第J光缆段的气动升力;
---第J+1光缆段的气动阻力;
---第J+1光缆段的气动升力;
VJ、VJ-1、Vwind、均为矢量。
2.如权利要求1所述的基于光缆敷设的共轴双旋翼无人直升机参数辨识试验***,其特征在于:所述整体光缆动力学模型建立方法包括:将整条光缆看作是由多个刚性杆通过绞进行连接的结构体,针对无人直升机布设光缆的实际状态,设定如下建模条件:
1.2.1由于光缆的第一段盘绕在恒张力铺设机构上,所以设定该第一段光缆段为固定端且固定于直升机平台之上,光缆最后一段为贯序着陆的自由端;
1.2.2自由端所受到的在x轴、y轴上力分量通过测量获得;
1.2.3直升机在正常布设光缆的过程中,其飞行高度保持固定,即在光缆全部释放之前,直升机平台所拖曳的停留在空中的光缆段的长度固定;
将各个光缆段间的运动约束关系转化为各个光缆段张力的递推关系,最后得出光缆段约束力Tj与***外力的线性关系;
对于结点Mj-1,其质量为mj-1,结点Mj,其质量为mj,设这两个结点所受的气动力分别为和则各个光缆段结点的运动方程为
其中,mj-1为j-1段光缆的质量,a(j-1)x为j-1段光缆x方向的加速度,a(j-1)y为j-1段光缆y方向的加速度,Tj-1为j-1段光缆所受的外力,Tj为j段光缆所受的外力,θj为j段光缆与飞行器平台z轴的夹角;
和分别是j-1所受的空气动力在x和y方向的分量;g为重力常数;式(3.9)描述最后一段的受力关系,其中Flastx、Flasty是自由端所受到的力在x和y方向的分量;
各个光缆段的加速度约束关系为:
其中:ajx为j段光缆在x方向的加速度,ajy为j段光缆在y方向的加速度,lj为j段光缆的长度,θj为j段光缆与飞行器平台z轴的夹角,ωj为j段光缆与飞行器平台z轴夹角的角速度;
将式(3.7)带入式(3.10)得到式(3.12),式(3.8)和(3.9)分别代入(3.11)式中,各个式子中的加速度项被消掉,获得各段光缆之间的铰约束力递推关系:
解上述线性方程组(3.12)(3.13)(3.14),经过化简后得到各光缆段结点的张力约束关系如下的线性方程:
Tj+1=Aj+1Tj+Bj+1Tj-1+Cj+1(j=2,3,...n-1) (3.15)
Tj代表结点Mj所受的张力;
根据(3.12)式,可以列出如下方程组:
Aj+1、Bj+1、Cj+1,均为系数,其取值与光缆段的长度、质量和姿态有关;从上述约束关系的推导过程,可以得出三个系数仅与光缆段的位移和速度有关;这说明在定高飞行情况下,当光缆释放速度、无人机平台速度以及风速恒定时,整个光缆及光缆段所受到的张力也恒定,这也是我们通过控制光缆释放速度来控制张力的原理;这个结果也表明,恒张力控制不仅仅是为了满足光缆强度要求Tj小于光缆规定的张力,同时保持光缆的张力使其具备规定的空中姿态,确保不与旋翼发生卷绕或干涉,此外也是为了降低光缆对飞行器平台的力与力矩的扰动。
3.如权利要求2所述的基于光缆敷设的共轴双旋翼无人直升机参数辨识试验***,其特征在于:对于用于光缆敷设的无人机平台来说,其动力学特征主要包括光缆的动力学和共轴双旋翼直升机各组成部件的动力学特性共同决定,共轴双旋翼无人直升机动力学建模中,光缆以及各个飞行器平台的部件气动特性与飞行器平台总体运动特性之间的关系;
将式(3.14)线性化后,将直升机的状态空间方程表达如下:
x=Ax+Bu,y=Cx (3.20)
其中:b-横向旋翼挥舞角;a-纵向旋翼挥舞角;
处理光缆影响的方法是将其作为直升机平台的一个外在扰动,通过对平台模型中相应气动导数增加修正量的方式,来获得与现实一致的模型响应;气动导数及其修正量可以通过下面的参数辨识的方法统一获得;
根据上述方法,确定的***模型结构为:
其中,Ba,Ab,La,Lb,Lu,Lv,Ma,Mb,Mu,Mv,Xu,Yv,Xa,Yb均为气动导数;这些也就是我们需要确定的参数,α/t,为光缆组件对旋翼挥舞角影响的修正量,g/k为光缆组件对飞行平台姿态角影响的修正量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410693197.1A CN104391364B (zh) | 2014-11-26 | 2014-11-26 | 基于光缆敷设的共轴双旋翼无人直升机参数辨识试验*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410693197.1A CN104391364B (zh) | 2014-11-26 | 2014-11-26 | 基于光缆敷设的共轴双旋翼无人直升机参数辨识试验*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104391364A CN104391364A (zh) | 2015-03-04 |
CN104391364B true CN104391364B (zh) | 2017-06-09 |
Family
ID=52609288
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410693197.1A Active CN104391364B (zh) | 2014-11-26 | 2014-11-26 | 基于光缆敷设的共轴双旋翼无人直升机参数辨识试验*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104391364B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106774282B (zh) * | 2017-02-09 | 2019-09-06 | 中国航天空气动力技术研究院 | 一种在线参数辨识电路 |
CN107238907B (zh) * | 2017-07-18 | 2024-05-10 | 沈阳亨通光通信有限公司 | 架空光缆敷设***及架空光缆敷设方法 |
CN111504341B (zh) * | 2020-04-30 | 2023-09-19 | 中国直升机设计研究所 | 一种直升机飞行状态识别方法 |
CN114253296B (zh) * | 2021-12-22 | 2024-06-07 | 中国人民解放军国防科技大学 | 高超声速飞行器机载轨迹规划方法、装置、飞行器及介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0966081A1 (en) * | 1998-06-19 | 1999-12-22 | ABB Sae SpA | Device for laying a towing cable for overhead power lines on supporting towers by means of a helicopter and related method |
CN201551845U (zh) * | 2009-09-05 | 2010-08-18 | 魏承赟 | 固定翼航模飞行辅助*** |
CN102566446A (zh) * | 2012-01-16 | 2012-07-11 | 南京航空航天大学 | 基于线性模型组的无人直升机全包线数学模型构建方法 |
CN103197683A (zh) * | 2013-04-18 | 2013-07-10 | 南京航空航天大学 | 四旋翼无人直升机飞行高度与姿态智能控制方法 |
CN103760906A (zh) * | 2014-01-29 | 2014-04-30 | 天津大学 | 神经网络与非线性连续无人直升机姿态控制方法 |
CN103885450A (zh) * | 2014-03-03 | 2014-06-25 | 天津大学 | 无人直升机姿态非线性控制方法及验证平台 |
-
2014
- 2014-11-26 CN CN201410693197.1A patent/CN104391364B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0966081A1 (en) * | 1998-06-19 | 1999-12-22 | ABB Sae SpA | Device for laying a towing cable for overhead power lines on supporting towers by means of a helicopter and related method |
CN201551845U (zh) * | 2009-09-05 | 2010-08-18 | 魏承赟 | 固定翼航模飞行辅助*** |
CN102566446A (zh) * | 2012-01-16 | 2012-07-11 | 南京航空航天大学 | 基于线性模型组的无人直升机全包线数学模型构建方法 |
CN103197683A (zh) * | 2013-04-18 | 2013-07-10 | 南京航空航天大学 | 四旋翼无人直升机飞行高度与姿态智能控制方法 |
CN103760906A (zh) * | 2014-01-29 | 2014-04-30 | 天津大学 | 神经网络与非线性连续无人直升机姿态控制方法 |
CN103885450A (zh) * | 2014-03-03 | 2014-06-25 | 天津大学 | 无人直升机姿态非线性控制方法及验证平台 |
Also Published As
Publication number | Publication date |
---|---|
CN104391364A (zh) | 2015-03-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Quindlen et al. | Flush air data sensing for soaring-capable UAVs | |
CN104503258B (zh) | 基于光缆敷设的共轴双旋翼无人直升机建模方法 | |
Arifianto et al. | Development and modeling of a low-cost unmanned aerial vehicle research platform | |
Owens et al. | Overview of dynamic test techniques for flight dynamics research at NASA LaRC | |
CN104391364B (zh) | 基于光缆敷设的共轴双旋翼无人直升机参数辨识试验*** | |
Tzanetos et al. | Future of Mars rotorcraft-Mars science helicopter | |
Tian et al. | UAV flight test evaluation of fusion algorithms for estimation of angle of attack and sideslip angle | |
Saderla et al. | Parameter estimation of UAV from flight data using neural network | |
Chambers | Longitudinal dynamic modeling and control of powered parachute aircraft | |
Carnduff | System identification of unmanned aerial vehicles | |
Jurisson et al. | Aeroservoelastic flight testing platform development for system identification | |
Nshuti et al. | Modeling, Simulation and Flight Testing to Support Proof of a Stratospheric Dual Aircraft Platform Concept | |
Jager | Test and evaluation of the Piccolo II autopilot system on a one-third scale Yak-54 | |
Mihai et al. | A practical method to estimate the aerodynamic coefficients of a small-scale paramotor | |
Richards et al. | Airworthiness evaluation of a scaled joined-wing aircraft | |
Howard | Automatic flight controls in fixed wing aircraft–The first 100 years | |
Krueger et al. | Flight loads analysis and measurements of external stores on an atmospheric research aircraft | |
Nagy et al. | Unmanned measurement platform for paragliders | |
Liwauddin et al. | Development of a Dodecacopter using Pixhawk 2.4. 8 autopilot flight controller | |
Hamada et al. | Development of an Unmanned Aerial Vehicle (UAV) Research Platform for Flutter Analysis | |
Rose et al. | Development and validation of flight dynamics model of a uav airplane | |
Sim | A correlation between flight-determined derivatives and wind-tunnel data for the X-24B research aircraft | |
Setlak et al. | Control system of the multi-rotor in flight in the presence of strong wind | |
Nshuti | Design of Flight Control Laws for a Novel Stratospheric Dual-Aircraft Platform | |
Schafer | Flight Test Measurement of Ship Airwake Disturbances using Small-Scale Rotorcraft |
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 |