共轴双旋翼动力学建模与验证
Dynamics Modeling and Validation of Coaxial Lifting Rotors
通讯作者: 黄健哲,男,副教授,博士生导师,电话(Tel.):021-34204112;E-mail:huangdeng9@sjtu.edu.cn.
责任编辑: 石易文
收稿日期: 2021-02-5
基金资助: |
|
Received: 2021-02-5
作者简介 About authors
胡金硕(1988-),男,山东省滨州市人,助理研究员,主要从事无人机动力学与控制.
建立共轴双旋翼动力学模型可用于共轴多旋翼飞行器的控制器设计、飞行模拟等研究,需要同时兼顾模型的计算实时性和精度.首先,基于伴随理论推导包含旋翼尾迹的全流场诱导速度计算模型.然后,在单旋翼有限状态入流模型基础上进行扩展,建立考虑尾迹倾斜的共轴双旋翼有限状态动力学模型.最后,推导共轴双旋翼悬停状态上下旋翼推力计算公式,并进行试验验证.研究结果表明:所建立的动力学模型计算复杂度在可接受的范围内,推力计算值与测试值在一定转速范围内基本吻合,并可很好地反映出共轴双旋翼推力损失趋势.
关键词:
The dynamics model for coaxial lifting rotors can be used to study the controller design and flight simulation for coaxial-rotor aerial vehicles. However, both the computational efficiency and the accuracy should be considered. First, the computational model of the induced velocity of lifting rotor everywhere including the wake region is derived based on adjoint theorem. Then, the finite state dynamics model for coaxial lifting rotors with wake skew considered is developed by extending the finite state inflow model for single rotor. Finally, the equations for calculating the thrust of coaxial lifting rotors in the hover condition are given, and the test is conducted. The results show that the computational complexity of the proposed dynamics model for coaxial lifting rotors is acceptable, and the computational thrusts are almost close to the test results when the rotational speed is within a certain range, which can also reflect the trend of the thrust lost for coaxial lifting rotors.
Keywords:
本文引用格式
胡金硕, 黄健哲.
HU Jinshuo, HUANG Jianzhe.
不同于旋翼计算流体力学(CFD)的研究,针对旋翼飞行器控制设计和飞行仿真的动力学模型需要有更好的实时性,同时对计算精度也有一定的要求.1989年文献[3]开发的Peters-He模型是目前最广泛应用的旋翼有限状态模型,可以用于计算旋翼盘上任意位置的诱导速度,进而实时计算单旋翼的推力.该理论隐含了著名的动态入流理论、Pradtl/Galdstein静态入流分布和Theodorsen理论,并且计算结果和试验数据有很好的吻合度.此后,一种被称为Morillo模型[4]的动态入流模型由Peters教授团队提出,它对旋翼上方流场的求解有其优越性.由于共轴双旋翼之间流场的干扰,上下旋翼的拉力性能将下降[5],所以用于计算单旋翼诱导速度和推力的Peters-He模型和针对旋翼上方流场计算的Morillo模型,它们均无法直接用于对共轴双旋翼的动力学建模.基于当前的有限入流建模理论,文献[6]利用关联上下旋翼入流系数和压强场系数的入流影响系数矩阵来计算每个旋翼产生的诱导速度.该入流影响系数矩阵包含了旋翼自诱导速度信息和上下旋翼之间互诱导速度信息,由于当时缺少互诱导速度的封闭计算公式,所以互诱导速度需要通过数值计算得到.其开发出的模型在XV-15倾转旋翼机上得到了试验验证.由于该模型需要进行数值计算,所以其计算的实时性能还有一定的提升空间.2015年,文献[7,8]提出了混合势流模型,建立了Peters-He模型和Morillo模型非线性混合函数,对旋翼盘上和旋翼上方流场的计算精度做了进一步提升,更重要的是,该模型对旋翼尾迹进行了建模,给出了旋翼尾迹流场计算的封闭形式解.文献[9]在旋翼悬停地效建模上应用了该混合势流模型,并与试验数据和RotCFD计算结果进行了对比,验证了模型有效性.文献[10]基于运动嵌套网格方法,建立了一套适合于悬停状态下共轴刚性双旋翼非定常干扰流场分析的计算流体力学方法.该研究运用了计算流体力学方法,可以给出共轴双旋翼更加细致的动态流场信息,提高了计算精度,但是存在很大的计算复杂度.文献[11]对倾转过渡状态共轴倾转旋翼进行了研究,考虑了尾流倾斜和收缩的上下旋翼的相互干扰,通过能量方程推导了诱导速度的计算公式,最后利用迭代方法求解了上下旋翼的载荷.由于整个计算过程需要迭代,所以通过此方法进行共轴多旋翼动力学建模也具有一定的计算复杂度.文献[12]结合叶素理论和Pitt-Peters 动态入流模型对共轴双旋翼的气动干涉和挥舞运动进行了近似建模,其他未准确建模的部分用控制算法进行补偿,并在Simulink搭建了控制器进行联合仿真.该方法利用动态入流模型计算诱导速度,仿真实时性好.但是Pitt-Peters 模型结构简单,是Peters教授团队最初提出的三状态动态入流模型[13],对于一些复杂的飞行状态难以保证较好的计算效果.此外,若仅采用诱导速度干扰系数表示上、下旋翼之间的气动干涉,需要较为充分的理论支撑.现有研究中主要考虑了上下旋翼的流场干扰,并对其进行了相应的数学建模,但还未见到有基于时滞伴随理论的势流模型来建立共轴双旋翼动力学模型的研究.
针对共轴双旋翼上旋翼入流受到下旋翼上游流场影响,以及下旋翼入流受到上旋翼尾流影响的问题,本文提出一种基于势流模型的共轴双旋翼动力学建模方法,并设计了共轴双旋翼的推力测试装置,验证了悬停状态下共轴双旋翼的推力计算准确性.本文研究可为共轴双旋翼飞行器的旋翼初步设计、控制和飞行仿真提供理论基础.
1 共轴双旋翼有限状态模型
根据文献[4],上下旋翼在其各自平面和上方位置 (ν≥0)产生的诱导速度v(u)(ν, η, ψ, t)和v(l)(ν, η, ψ, t) 分别表示为
式中:
式中:x,y,z为计算点在体坐标系中的正则化坐标位置,则在旋翼边缘有A=1.
由于上旋翼的尾迹会对下旋翼产生影响,所以需要给出上旋翼的尾迹速度的计算公式.根据文献[8],上旋翼在 (ν, η, ψ)位置的尾迹速度v(u)(ν, η, ψ, t)(v<0) 表示为
式中: v(u)(0, η0, ψ0, t-τ)为上旋翼在(0, η0, ψ0)位置的诱导速度;ˉω(u)(0, η0, ψ0+π, t-τ)为上旋翼在(0, η0, ψ0+π)位置的伴随速度;τ为时滞量.由于椭圆坐标系与体坐标系的转换比较复杂,η0、ψ0与η、ψ 之间的关系可用如下关系式进行表述:
如果(x, y, 0)→(0, η, ψ),那么(x+
其中:
式中:
如果诱导速度和伴随速度计算式包含 N次谐波,即m=0, 1, …, N, 定义诱导速度因子和伴随速度因子的状态向量为 a(u)、a(l)、Δ(u). 根据文献[7],压强势用于计算旋翼周围的压强,而压强势可用压强因子
其中,当 m+n为偶数时,
式中: Q和R分别为旋翼的叶片数和半径;ρ为空气密度;
通过低速流体的质量守恒和动量守恒定律公式,代入用勒让德函数表示的速度和压强表达式[7],并通过伽辽金变换可得到一组以诱导速度因子和伴随速度因子为状态量的常微分方程:
其中
K(u)=DV(u)L−1(u)M
K(1)=DV(l)L−1(l)M
˜τ(u)=(−1)n+1τm(u)n
式中:M和D分别为表观质量矩阵和阻尼矩阵;
其中
式中: Ua和Va分别为体坐标系下相对气流速度的x和y分量;
由于式(9)中的入流矩阵需要 τ时刻前的伴随速度因子信息,所以需要同时对式(9)进行时间正向数值积分,并对式(10)采取计算时刻t1为初始的时间反向数值积分,同时求得速度因子和伴随速度因子的时间序列信息.对于式(10)的时间反向数值积分,需要在每一个计算时刻进行时长为τ的反向数值迭代,计算量比较大.通过以下计算步骤可以通过更少的计算量,求得在t1时刻伴随速度因子τ时刻前的数值Δ(u)(t1-τ).
步骤1 求解 R-1M-1K(u)的特征值和特征向量分别为λi(i=1, 2, …, G),φ, 使得:
φ-1(R-1M-1K(u))φ=
其中:G为矩阵M的维数.
步骤2 给出 t1时刻前
步骤3 通过下式计算Δ(u)(t1-τ):
式(8)中,求解
式中:
2 计算结果验证
为了验证共轴双旋翼诱导速度动力学模型的正确性,选取悬停状态上、下旋翼的推力作为研究对象,搭建试验验证平台,平台实物图如图1所示.右侧旋翼为固定安装,旋翼轴高度和轴向位置不可调.左侧旋翼可沿轨道向前向后移动,用于调整双旋翼的轴向距离;左侧旋翼的离地高度也可进行调整,当高度调整至右侧旋翼高度相同时,为共轴状态.选取型号为1024的PVC螺旋桨作为测试桨叶,该桨叶直径约为254 mm,距桨叶中心点5 mm半径范围内为平直安装面(矩形翼型、无扭转角).桨叶桨距为119.38 mm,弦长沿桨叶径向位置连续可变,为了简化计算,将桨叶从平直安装面结束段开始,沿径向方向分为5部份,分别测量每一部分边缘位置的弦长,数据如表1所示.沿桨叶径向方向距离桨叶中心为r的截面,其扭角可根据下式计算:
式中:ϑ为桨叶桨距.两侧旋翼分别由两套完全独立的动力系统提供旋转动力,两套动力系统的电机也由完全独立的控制系统来分别调节脉冲宽度调制(PWM)值,进而实现电机转速的独立控制.
图1
表1 测试旋翼叶片径向弦长分布
Tab.1
r/mm | c/mm |
---|---|
5.0 | 13.0 |
25.4 | 20.8 |
50.8 | 29.7 |
71.2 | 32.5 |
91.6 | 31.0 |
117.0 | 24.8 |
调整左、右电机旋转方向和旋翼桨叶安装方式,使得左、右旋翼旋转方向相反,且尾迹均朝左.此时,左旋翼处于右旋翼的尾迹(下方流场)中,可视为共轴双旋翼的下旋翼.首先,只为一侧旋翼电机通电,使其旋转并测量和记录不同转速下产生的推力,用于校验计算模型对单旋翼计算的准确性, 也为双旋翼计算提供准备.由于一侧旋翼旋转时,产生的诱导速度会成为另一侧旋翼的入流,带动另一侧旋翼随之旋转.所以,可连通右侧旋翼, 即使左侧旋翼被动旋转,右侧旋翼处于左侧旋翼的上游流场,产生的影响较小.图2(a)和2(b)分别给出了单轴旋翼旋转时转速和对应推力测试结果的时间历程图,其中:F为推力测试值;ω为旋翼转速.测试采取逐渐增加转速的方式,并在每个转速停留一段时间使流场充分发展.截取每一转速段中间时刻0.8 s长度数据,分别计算平均转速和平均推力,并用上一节提出的方法对不同转速下的单轴旋翼进行计算.这里估计旋翼升力线斜率为6.11 rad-1,且因为大于 5000 r/min后转速波动过于激烈,而 1800 r/min时旋翼推力极小,使得力传感器几乎感受不到旋翼产生的推力,因此只对 3000 r/min和 5000 r/min之前的数据进行计算和对比.计算值与测试值的对比如图3所示.在不同转速下,推力的计算值和测试值的偏差如表2所示.其中:
图2
图3
采用MATLAB中的Curve Fitting Tool工具包,对数据进行拟合,可得到转速与推力的特性关系式,如下式所示:
表2 单旋翼不同转速下的推力
Tab.2
ω/(r·min-1) | F/g | eF/g | |
---|---|---|---|
3000 | 65.7 | 111.80 | 46.10 |
3500 | 128.0 | 152.23 | 24.23 |
4000 | 200.0 | 198.81 | -1.19 |
5000 | 393.0 | 310.59 | -82.41 |
对于测试数据,拟合的误差平方和(SSE)为 0.7526,而对于计算数据,拟合的误差平方和为1.908×10-7.由式(16)可知,推力系数的计算值和测试值仅有9%的偏差.
对于共轴双旋翼的测试,设置旋翼轴向距离为一个旋翼半径长度,即127 mm.同时调节左、右旋翼电机PWM值,使得左、右旋翼转速基本稳定在某一设定值附近,并分别记录下对应的PWM值.采用试验出的PWM值,同时对左、右旋翼进行升转速测试,并测试对应的左、右旋翼推力值,如图4所示.由于实际条件限制,左、右电机不可能同时开机,并且升转速的起始点也无法做到完全同步,从时序图上看,出现不完全同步的现象.因此,可选取左、右电机转速同步时间段的转速和推力数据进行平均,获得不同转速下左、右旋翼平均推力的测试值.利用上一节推导的模型,对相应的测试工况进行计算,测试值与计算值的对比如图5所示.不同转速下左、右旋翼的计算值、测试值和偏差如表3所示.其中:
图4
图5
图5
共轴双旋翼推力计算结果(轴距127 mm)
Fig.5
Computational thrusts results of coaxial lifting rotor (axial distance=127 mm)
表3 双旋翼不同转速下的推力
Tab.3
ω/(r·min-1) | F(l)/g | F(u)/g | ||||
---|---|---|---|---|---|---|
1700 | 39.30 | 33.11 | -6.19 | 12.75 | 20.99 | 8.24 |
3000 | 129.21 | 103.11 | -26.10 | 64.34 | 65.36 | 1.02 |
3700 | 227.02 | 156.87 | -70.15 | 135.82 | 99.44 | -36.38 |
4400 | 307.86 | 221.80 | -86.06 | 216.77 | 140.60 | -76.17 |
4900 | 395.95 | 275.13 | -120.82 | 276.72 | 174.41 | -102.31 |
由图5可知,在 3000~5000 r/min转速范围内,相较于单轴旋翼,对于共轴双旋翼左、右旋翼在相同转速下均出现推力损失,并且由于左旋翼位于右旋翼尾迹之中,造成的推力损失更大.所提出的模型可以较好的反应出这种规律,并且趋势与测试值大致相同.和单轴旋翼类似,随着转速的提高,推力测试值越高于计算值.同样地,通过数据拟合,可以分别得到左、右旋翼的转速与推力特性关系式,如下式所示:
对于测试数据,右旋翼拟合的SSE为 0.02175, 左旋翼拟合的SSE为 0.1522.对于计算数据,右旋翼拟合的SSE为1.354×10-7,左旋翼拟合的SSE为6.538×10-8.
3 结语
本文研究了基于有限状态模型的共轴双旋翼动力学建模方法,利用伴随理论推导了上旋翼尾迹诱导速度表达式,拓展了传统的单轴旋翼有限状态模型.上旋翼动力学模型中的入流矩阵包含了上旋翼的平均自诱导速度信息和下旋翼的平均互诱导速度信息,下旋翼动力学模型中的入流矩阵包含了下旋翼的平均自诱导速度信息和上旋翼的平均互诱导速度信息,其中上旋翼的平均互诱导速度需要利用上旋翼尾迹诱导速度计算式来求取.动力学模型的载荷输入表达式也包含了诱导速度信息,由此可见,该动力学模型是入流和载荷、上旋翼和下旋翼耦合的强非线性动力学模型.跟一般的经验模型不同,此模型需要输入共轴双旋翼桨叶几何信息、气动参数、轴向间距、转速等实际物理信息,因此更具有普适性,能够适应更多的应用环境.通过试验数据验证,此模型具有一定的计算准确度,可用于共轴双旋翼的设计初步计算、控制器设计与实时仿真,有一定的工程应用价值.
参考文献
共轴双旋翼直升机的技术特点及发展
[J].
Technology characteristic and development of coaxial rotor helicopter
[J].
Advanced UAV aerodynamics, flight stability and control
[M].
Finite-state induced-flow model for rotors in hover and forward flight
[J].DOI:10.4050/JAHS.34.5 URL [本文引用: 1]
Velocity field above a rotor disk by a new dynamic inflow model
[J].DOI:10.2514/2.3019 URL [本文引用: 2]
悬停共轴双旋翼干扰流动数值模拟
[J].
Numerical simulation of interaction unsteady flows around co-axial rotors in hover
[J].
Potential-flow inflow model including wake distortion and contraction
[D].
Converged velocity field for rotors by a blended potential flow method
[J].
悬停状态共轴刚性双旋翼非定常流动干扰机理
[J].
Unsteady flow interaction mechanism of coaxial rigid rotors in hover
[J].
共轴倾转旋翼性能计算方法
[J].
Mathematical model for performance of coaxial tilt-rotor
[J].
一种共轴双旋翼飞行器悬停控制联合仿真
[J].
Hover control co-simulation of a coaxial dual-rotor aircraft
[J].
Theoretical prediction of dynamic-inflow derivatives
[J].
/
〈 |
|
〉 |
