近年来,如何解决城市道路交通拥堵,提高包括救援、巡逻在内的各种社会保障工作的效率,以及提升交通生活品质等问题得到了大众的广泛关注。本文就此提出了一种单人全电涵道旋翼飞行器,根据使用要求,该飞行器具有垂直起降、体积较小、低速前飞、安全性高、可靠性高以及操纵简单和控制便捷等优点。
本文主要针对此新型涵道飞行器的总体/飞控一体化设计进行相关研究。首先,针对飞行器的使用要求,确定主要设计要求,提出3种构型方案。利用实验室样机的设计经验,建立总体设计模型和控制模型。然后,在两种模型的基础上,采用层次分析法与灰色模糊综合评估的方法对各方案进行评估,确定方案。
1 总体布局 1.1 设计要求依据飞行器的使用环境和使用方式,对该新型飞行器提出以下设计要求:(1) 可垂直起降;(2) 操控简易可靠;(3) 能安全搭载1名乘员(质量不低于65 kg);(4) 外形小巧,可车载,宽度不大于2.0 m,长度不大于3.5 m[1];(5) 全电设计,提高维修性。
1.2 布局设计为乘员安全飞行,本文采用涵道旋翼作为升力机构。布局设计主要是指升力机构的数目和布置方式。
根据设计要求,可采用以下基本布局:(1) 采用一个涵道共轴双旋翼作为升力装置,将驾驶舱置于其下,通过变距的操纵控制方式实现其反扭矩平衡。(2) 采用2个或3个涵道旋翼系统前后布置,驾驶舱置于中部,如图 1所示。(3) 采用4个涵道旋翼系统,对角线布置,如图 2所示。(4) 采用4个以上涵道旋翼系统,如图 3所示。从控制角度考虑,仅第3、4类布局操控简易可靠,第1、2类布局操纵通道耦合严重,甚至可能会因操纵功效小而无法实现功能,故舍弃第1、2类布局。
上述的方案思路中,相较于变距的操控方式,变转速方式的操纵机构简单,可靠性增强。在涵道旋翼系统中,相同桨盘面积的涵道共轴双旋翼系统产生的升力大于单旋翼涵道系统,但是,其气动效率低于涵道单旋翼系统,且共轴涵道系统高度更长,阻力更大,不利于涵道装置小角度前飞。过多的升力机构会复杂操纵设计、增加结构质量以及降低可靠性。经过初步比较,本文采用4个涵道单旋翼系统作为升力装置,采用4个电机,电池被安装在机体中央下方。其优点是:平衡反扭矩易行;可以只通过调节旋翼转速差动来控制飞行器的姿态,控制方式难度降低,简化了机械部件,减轻了飞行器的质量,提高了操纵效率;动力装置相互独立,没有复杂的传动机构,且维护成本和时间都大大降低。
涵道旋翼基本布局可以是“十”字形、“X”形或“H”形。图 4列出了“十”字形和“H”形。本文相邻的两个涵道设计成“8”字形涵道,使其在受限尺寸内拥有大桨盘面积,获得最大的升力。
2 总体模型 2.1 特性参数分析
本文总体模型综合了气动、质量、性能及尺寸等特性模型,是用来计算分析飞行器的气动布局、重心、飞行性能及总体布置效能的综合模型,该模型受以下参数影响。
(1) 质量分析。总质量是飞行器的主要参数,决定了飞行器的尺寸和大小。本文飞行器样机的质量组成如图 5所示。新型涵道四旋翼飞行器的结构和装置都比较简洁,统计数据的参考价值相对较低,本文主要参考实验室类似涵道旋翼样机的质量组成数据。设计要求有效载荷不低于65 kg,所以,原理样机总质量初步设定为125 kg。
(2) 旋翼参数分析。本文采用经验公式与试验数据相结合的方式进行参数选择。经验公式参考文献[2],反映了桨盘载荷DL、旋翼实度σ等参数和设计起飞重量W0的关系。试验[3]和理论计算[4]发现,涵道单旋翼系统中,涵道升力占总升力的20%以上,涵道内的旋翼拉力要略低于同等条件下的孤立单旋翼的拉力值。据此,本文在叶素动量理论计算旋翼拉力的基础上估算系统总拉力。最终得到旋翼总距、半径R、弦长c和片数N等参数。
(3) 涵道参数分析。单涵道的结构示意图如图 6所示,主要设计参数如下:涵道内径dd、出口直径为Dd、涵道高度L、涵道唇口半径rd、涵道最大壁厚bd、涵道锥角β、桨尖与涵道间隙δtip以及涵道入口与桨盘之间的距离为hd。表 1列出获得较佳气动效率时的涵道各参数[5-7]的取值状况[8]。
两相邻涵道也可以设计成如图 7所示的“8”字形涵道,内部部分重叠的双旋翼模型如图 8所示,该模型主要考虑轴向间距h和纵向间距d。
通过实验室的相关计算结果可以得出:h=0.2R,d大于1.8R时,升力值与孤立状态时的升力值相差不大;d=1.5R,h从大减小到0.2R时,升力变化非常小。与文献[1, 9]的结论基本相符。
上述各参数中,L/dd,bd和h的值对质量影响较大。L/dd的取值决定了涵道体的基本形状,并对涵道体的气动效率、迎风阻力以及动量阻力有一定影响,当其值处于表 1的取值范围时,可以获得较高的气动效率和较低的阻力。适当增大bd可加强涵道体的增升效应,但是bd值越大,涵道本体的质量也越大。h值越大,上下旋翼的气动干扰越小,但是h值的增大一方面会增加旋翼轴和支撑机构的质量,另一方面会增加涵道体的高度,从而增加涵道本体质量。
(4) 功率分析。功率是计算飞行性能的基础,本文以经试验验证过的基于数值积分的叶素动量理论模型计算各状态功率,高精度地描述出飞行器飞行性能。电机的功率及传动损失也是经试验测试获得的。本文涵道单旋翼升力系统在设定起飞质量试验中测得的悬停状态功率载荷为4.58 kg/kW。
2.2 总体设计方案及综合模型在基本布局中,以下列3个方案进行设计分析。方案A:4个涵道旋翼“十”字形布置,前后两个大涵道,左右两个小涵道;方案B:“X”形布置4个升力系统,对角线呈90°;方案C:“H”形布置4个升力系统。各方案的主要参数如表 2所示,部分对计算和评估影响甚微的参数不列出。其中,质量分配部分,针对质量计算与控制要求,在质量分配上留有5%的余量,所以,总质量为120 kg左右,设计以125 kg考虑。
本文总体综合模型如下
$\begin{array}{l} F\left( x \right) = {a_1}{f_1}\left( x \right) + {a_2}{f_2}\left( x \right) + {a_3}{f_3}\left( x \right) + \\ \;\;\;\;\;\;\;\;\;\;\;{a_4}{f_4}\left( x \right) \end{array}$ | (1) |
式中:fi(x)分别对应质量特性、任务能力、空间尺寸和气动特性模型;ai分别为各特性在综合模型中的权重。质量特性模型中主要包含质量效率、转动惯量和最大起飞质量;任务能力模型包含最大速度、有/无地效悬停升限、自转性能、功率载荷等;气动特性模型包含悬停效率及桨盘载荷;空间尺寸模型包含长宽高尺寸。各特性模型均是以特性参数代入专业计算模型获得。
总体模型将与控制模型一起放入评估模型中完成方案选择。
3 控制模型本文飞行器强调全电设计和操控简易可靠,故既要选择易控构型或布局,又要设计好飞控系统。动力学建模是研究涵道旋翼飞行器控制模型的基础,在此基础上进行配平和控制律的设计等工作,并以配平中的各状态配平值和控制效果为控制模型的主要评估因素。
3.1 飞行器动力学模型(1) 坐标系
地轴系ODXDYDZD,即地面坐标系,其定义如下:OD:地面预选一点;XD:水平面上的初始航向,与正北的偏角为航向角;ZD:铅垂线;YD:横轴,由右手法则确定。
体轴系OXYZ的定义如下:O:飞行器重心;X:纵轴,与飞行器的构造轴线平行,指向飞行器的机头方向;Z:在飞行器垂向对称面内,垂直纵轴向下;Y:横轴,由右手法则确定。
(2) 数学模型
设飞行器的质量为m,作用于质心的合外力为F。飞行速度V在体轴系中的xyz三个方向分量为u,v,w,角速度分量为p,q,r。Ix,Iy,Iz为飞行器对体轴系3个方向轴的惯性矩,Izy,Iyz,Ixz为惯性积。9个运动方程如式(2~4) 所示[10]
$\left\{ \begin{array}{l} \dot u = \frac{{{F_x}}}{m} + vr-wq-g\sin \theta \\ \dot v = \frac{{{F_y}}}{m} + wp-ur + g\cos \theta \sin \phi \\ \dot w = \frac{{{F_z}}}{m} + uq - vp + g\cos \theta \cos \phi \end{array} \right.$ | (2) |
$\left\{ \begin{array}{l} \dot p = \frac{{{I_z}\left[{{M_x} + \left( {{I_y}-{I_z}} \right)qr + {I_{xz}}pq} \right] + {I_{xz}}\left[{{M_z} + \left( {{I_x}-{I_y}} \right)pq-{I_{xz}}qr} \right]}}{{\left( {{I_x}{I_z} - I_{xz}^2} \right)}}\\ \dot q = \frac{1}{{{I_y}}}\left[{{M_y} + \left( {{I_z}-{I_x}} \right)pr-{I_{xz}}\left( {{p^2}-{r^2}} \right)} \right]\\ \dot r = \frac{{{I_{xz}}\left[{{M_x} + \left( {{I_y}-{I_z}} \right)qr + {I_{xz}}pq} \right] + {I_x}\left[{{M_z} + \left( {{I_x}-{I_y}} \right)pq-{I_{xz}}qr} \right]}}{{\left( {{I_x}{I_z} -I_{xz}^2} \right)}} \end{array} \right.$ | (3) |
$\left\{ \begin{array}{l} \dot \phi = p + \tan \theta \left( {q\sin \phi + rcos\phi } \right)\\ \dot \theta = q\cos \phi-r\sin \phi \\ \dot \psi = \left( {q\sin \phi + r\cos \phi } \right)/\cos \theta \end{array} \right.$ | (4) |
通常将式(2~4) 的方程组表示为
$\boldsymbol{\dot x = f}\left( {\boldsymbol{x}, \boldsymbol{u}, t} \right)$ | (5) |
式中:
一般非线性模型不利于模型分析,基于小扰动假设[11],可将其线性化,得到
$\left\{ \begin{array}{l} \boldsymbol{\dot x} = \boldsymbol{Ax} + \boldsymbol{Bu}\\ \boldsymbol{y} = \boldsymbol{Cx} + \boldsymbol{Du} \end{array} \right.$ | (6) |
式中:A和B为重点需要获取的矩阵;C为输出矩阵;D一般取零矩阵。
3.2 配平分析和计算配平的任务是根据给定的状态条件,确定涵道旋翼飞行器在该种状态下稳定飞行时的操纵输入量和飞行姿态角。在MATLAB/Simulink中搭建飞行器模型,如图 9所示。
用MATLAB提供的配平函数trim()配平。在不同条件下配平时,配平值是评估控制模型的因素之一,在上述方案中,特别需要观察4个旋翼的转速随飞行器飞行状态的变化情况。
3.3 控制器设计(1) 模型线性化
以配平结果中的状态量和操纵量为基准状态,利用函数linmod()得到线性方程的状态空间形式。
(2) 控制律设计
由于设计处于概念设计阶段,所以本文的控制模型采用PID控制方法。这里采用最常用的串级反馈回路设计控制系统。在Simulink中搭建出增稳回路和姿态回路,如图 10所示,其中,内环增稳控制回路主要是角速度控制回路,增加系统的稳定性。外环姿态回路主要是滚转角、俯仰角和偏航角的控制回路,是整个控制系统的核心部分。
以方案C为例,在滚转角、俯仰角通道和滚转通道加入方波信号,完成参数的整定工作。响应曲线如图 11~13所示。
对控制效果的评估就是从衰减比、静态偏差、超调量以及调节时间等参数出发,对其进行评估。这些指标有些是相互矛盾的,比如超调量和过渡时间。所以,在进行PID整定时,要根据要求对上述指标有所侧重,然后进行控制模型评估。
4 飞行器设计方案评估本节依据上述总体和控制模型进行综合评估,从而选出最佳的设计方案进行下一步研究。
评估体系分为3层:顶层为最终评估结果,中间层为总体模型和控制模型的性能,底层为各模型需要考虑的指标。文中将质量特性、任务能力、空间尺寸和气动特性等参数作为总体模型的底层指标。将配平情况和控制响应情况作为控制模型的底层指标。
总体模型与控制模型存在一定的耦合关系,两者中的部分参数互相影响。一方面总体模型中涵道旋翼的参数设计影响模型的气动分析和质量分析,而控制模型中的动力学建模和配平等计算一般依赖于总体模型传递出来的气动分析和质量分配。另一方面,控制模型中的操纵工作又会对气动力和力矩造成影响。图 9中利用s函数构造的4个升力机构气动模块里的参数主要是总体模型中涵道旋翼的相关参数,运动方程中的转动惯量和惯性积又主要受总体模型中的质量分配影响。在控制模型中,操纵值影响旋翼转速等参数,进而影响气动力和气动力矩。本文通过线性方程体现其耦合关系。
首先,采用层次分析法确定各层指标的权重系数,具体流程如图 14所示。
模型中的指标可分为定量指标和定性指标,为了消除指标量纲对评估的影响,需对其进行无量纲化处理。定量指标无量纲化时,将当前指标值与所有方案中的最优值相比。定性指标的度量采用“很差、较差、差、一般、较好、好、很好”七级打分制,并用梯形模糊数[12]将每个等级表述出来,无量纲化值为其整体期望值。以“一般”和“较好”两个等级为例,无量纲化过程如表 3所示。
在求出权重和无量纲化后的指标值后,采用灰色模糊综合评估的方法进行评估,如图 15所示。
对前3个设计方案进行评估,3个方案的优属度分别为:0.451 7,0.649 2,0.792 1。方案B和方案C在评估中较优,方案C最好。方案C飞行器总体布置如图 16所示。
根据本文评估结果选定方案C。在此基础上,利用本文提出的方法对方案C进行进一步设计。考虑旋翼轴向间距、纵向间距和旋翼半径的值,将其组合成若干方案再次代入总体模型。在控制模型中加入操稳特性的分析,在新的总体/飞控一体化模型上完成方案C的设计、评估和选择。
5 结束语本文提出一种新型涵道四旋翼飞行器方案设计方法。总体模型中对特性参数进行分析,进而对不同构型的涵道四旋翼飞行器进行参数设计。结合总体模型的参数,在MATLAB/SIMULINK平台上完成控制模型的搭建,其子模块的参数易修改,仿真结果易得到。在控制模型中比较了各状态的配平结果和飞行器不同构型的控制效果。利用层次分析法和灰色模糊综合评估的方法对3种构型方案进行了评估。
[1] | AVERA M, KANG H, SINGH R. Performance and controllability assessment of an overlapping Quad-Rotor concept[C]//American Helicopter Society 72nd Annual Forum. West Palm Beach, Florida:[s.n.], 2016. |
[2] |
彭名华. 直升机总体多学科设计优化研究[D]. 南京: 南京南京航空航天大学, 2009.
PENG Minghua. Research on helicopter conceptual/preliminary multidisciplinary design optimization[D]. Nanjing: Nanjing University of Aeronautics & Astronautics, 2009.http://cdmd.cnki.com.cn/Article/CDMD-10287-1011253994.htm |
[3] |
李彬. 高速型无人涵道旋翼飞行器总体设计技术研究[D]. 南京: 南京航空航天大学, 2012.
LI Bin. Research on technologies of high-speed ducted-rotor UAV preliminary design[D]. Nanjing :Nanjing University of Aeronautics & Astronautics, 2012.http://cdmd.cnki.com.cn/Article/CDMD-10287-1012041217.htm |
[4] |
刘沛清.
空气螺旋桨理论及其应用[M]. 北京: 北京航空航天大学出版社, 2006.
LIU Peiqing. Theory and application of the airscrew[M]. Beijing: Beihang University Press, 2006. |
[8] | MARTIN P, TUNG C. Performance and flowfield measurements on a 10-inch ducted rotor VTOL UAV[J]. Journal of the American Helicopter Society, 2004. |
[9] | OTSUKA H, NAGATANI K. Thrust loss saving design of overlapping rotor arrangement on small multirotor unmanned aerial vehicles[C]//IEEE International Conference on Robotics and Automation. Stockholm, Sweden:[s.n.], 2016. |
[10] |
牛金宇. 涵道共轴旋翼飞行器建模与验证[D]. 南京: 南京航空航天大学, 2015.
NIU Jinyu. Modeling and verification of a ducted coaxial rotor aircraft[D].Nanjing: Nanjing University of Aeronautics & Astronautics, 2015.http://cdmd.cnki.com.cn/Article/CDMD-10287-1015952048.htm |
[11] |
高正, 陈仁良.
直升机飞行动力学[M]. 北京: 科学出版社, 2003.
GAO Zheng, CHEN Renliang. Helicopter flight dynamics[M]. Beijing: Science Press, 2003. |
[12] |
侯家麟, 唐慧娟, 于九如.
模糊数量化理论的应用与探索[J]. 系统工程理论与实践, 2003, 6(6): 107–112.
HOU Jialin, TANG Huijuan, YU Jiuru. Exploration and application for the fuzzy quantitative theory[J]. Systems Engineering Theory and Practice, 2003, 6(6): 107–112. |