农用多旋翼无人机作业能耗白盒模型构建与试验_李继宇.pdf

2 0 2 3 年 8 月农 业 机 械 学 报第 54 卷 第 8 期doi:106041/j issn 1000-1298 2023 08 008农用多旋翼无人机作业能耗白盒模型构建与试验李继宇张植勋赵胤伟沈卓刘知杰吕佳(华南农业大学工程学院，广州 510642)摘要:为了减少因能量有限而导致无人机飞行时间有限的问题，研究农用无人机在飞行时的作业能耗是十分必要的。本研究以农用多旋翼无人机为研究对象，采用白盒建模的方法，依照动力系统各部件原理，构建了农用无人机的飞行总效率模型，再根据农用无人机的飞行特点，建立了一种针对农用无人机的作业能耗模型，并通过对农用无人机的速度、载荷和航程这 3 个动态参数分组试验得到的数据进行了验证。结果表明，模型具有较高的精度，最大平均误差约为 6.582%，最大绝对误差中位数约为 7.654%。最后对模型的各个参数进行分析，引入模型修正系数对模型进行修正，修正后的模型精度相比修正前的最大平均误差减少约 3.092 个百分点，误差中位数减少约 3.612个百分点，修正后效果显著。本研究构建的理论模型可以用于农用无人机作业能耗的计算和预测，在规划无人机任务前就可以得到相应的作业能耗并进行优化，也可以适配不同控制器的其他旋翼机型，具有一定的应用价值。关键词:农业;无人机;能耗;白盒模型中图分类号:S252文献标识码:A文章编号:1000-1298(2023)08-0087-10OSID:收稿日期:2022 09 09修回日期:2022 12 10基金项目:广东省自然科学基金项目(2021A1515010736)和广州市重点研发计划项目(2023B03J1323)作者简介:李继宇(1979)，男，教授，博士生导师，主要从事精准农业航空研究，E-mail:lijiyu scau edu cn通信作者:吕佳(1979)，女，讲师，博士，主要从事精准农业航空研究，E-mail:lvjia scau edu cnConstruction and Experiment of White Box Model for EnergyConsumption of Agricultural Multi-rotor UAV OperationLI JiyuZHANG ZhixunZHAO YinweiSHEN ZhuoLIU ZhijieL Jia(College of Engineering，South China Agricultural University，Guangzhou 510642，China)Abstract:In order to reduce the problem of limited flight time of UAVs due to limited energy，it isnecessary to study the operational energy consumption of agricultural UAVs during flightTakingagricultural multi-rotor UAV as the research object，the white-box modeling method was used to constructthe total flight efficiency model of agricultural UAV according to the principle of each component of thepower system，and then an operational energy consumption model was established for agricultural UAVaccording to the flight characteristics of agricultural UAV，and the data obtained were verified by groupingtests of three dynamic parameters，namely speed，load and range of agricultural UAV The resultsshowed that the model had a high accuracy，with a maximum average error of about 6.582%and amedian maximum absolute error of about 7.654%Finally，each parameter of the model was analyzedand the model correction coefficient was introduced to correct the model，and the accuracy of thecorrected model was reduced by about 3.092 percentage points of the maximum mean error and about3.612 percentage points of the median error compared with the maximum mean error before thecorrection，which was significant after the correction The theoretical model constructed can be used forthe calculation and prediction of agricultural UAV operation energy consumption，and the correspondingoperation energy consumption can be obtained and optimized before planning the UAV mission，and it canalso be adapted to other rotor models with different controllers，which had certain application valueKey words:agriculture;UAV;energy consumption;white-box model0引言当前我国农用无人机在所有机型中超过 80%，农用无人机的优势主要体现在:用途广泛，无人机在农业领域可以进行喷药施肥1、授粉播种2、农情监测3、生物检测4 和放牧喂食5 等任务。使用方便效率高，其作业地形不受限制，还可以通过航线完成监测和撒药等作业，植保用药效率高，可以最大程度减少病虫害;经济效益明显，例如其单位植保费用相比传统方法低 30%50%，效率是人工的20 50 倍。农用无人机根据类型可以分为单旋翼、多旋翼、固定翼，按旋翼类型目前主要有 4 旋翼、6 旋翼、8 旋翼及更多旋翼构型，目前大规模使用的是 4 旋翼和6 旋翼。根据能源类型可以分为电动和油动，其中无人机电动转换效率为 70%，油动转换效率为30%，因此本文针对以锂电池为动力的农用无人机进行研究。目前商品锂电池的能量密度普遍低于200 W h/kg，使得农用无人机配套的电池可供续航时间最多为 20 min，难以满足当前农用无人机对长航程、长航时、大负载的作业需求。虽然增加电池数量和减少载荷可以提高续航能力，但前者增加重量，导致额外的能量消耗，后者降低了效率6。为了降低因电池容量有限而导致无人机飞行时间有限的影响，对农用无人机在飞行时的作业能耗进行研究是十分必要的。目前针对农用无人机作业能耗研究主要围绕黑盒或白盒建模。黑盒建模主要通过数据预处理、回归法、实验法和测试法等针对整体系统的各个模块实际的作业能耗数据进行采集，并拟合成有关函数进行建模。ABEYWICKAMA 等7 通过使用实验法对无人机各种状态下电池使用情况进行研究，得到一个完整的无人机能耗模型，但该方法只针对单款无人机在单种电池使用情况下进行研究，缺乏对无人机系统的原理性分析。PASETIA 等8 通过使用数据收集、数据预处理和回归法提出了一种基于任务的无人机能耗预测模型，然而该方法缺乏对无人机能耗规律的研究。林晋立等9 通过测试法搭建测试平台构建了针对多旋翼农用无人机功率能耗计算的理论模型，不过该方法并未通过实际飞行测试，只使用试验架进行固定测试。ABUK 等10 通过黑盒模型中经验测量法得出每个子系统能耗对整体能耗的影响，但该研究缺少对电机这个消耗最大的子系统进行进一步原理性分析。TAUB11 通过测试法使用不同种类锂电池进行风洞研究以评估电池放电行为，得到一个准确的续航估计值，通过研究电池在无人机飞行过程的能耗消耗情况对农用无人机作业能耗的研究具有一定的借鉴意义。黑盒建模研究虽然可以更加准确得到无人机的模型，但只是针对被测无人机进行建模，且需要大量重复的试验，存在一定局限性。白盒建模主要通过研究整体性系统的各个模块(包含电机、螺旋桨、电池、电调等)或者飞行任务下的原理推导进行建模。DOLING 等12 通过整体建模方法推导得出多旋翼无人机在悬停时的能耗，主要重点是将消耗的功率与无人机重量联系起来，没有 分 析 在 飞 行、起 飞 或 着 陆 期 间 的 能 耗。ACHTELIK 等13 通过研究微型飞行器上 4 个以上螺旋桨和电机的影响，针对几种类型的多旋翼进行建模，以比较它们的动力学、效率和冗余，然而该方法只针对螺旋桨和电机进行建模分析，缺乏对整个动力系统进行建模和分析。SHI 等14 通过研究螺旋桨、电机、电调、电池这 4 个动力系统的主要部件，得到适合多种机型性能综合评价方法，不过该方法仅针对无人机进行静态模型建模，没有对有关动态模型建模的进一步分析。BEZZO 等15 研究飞行任务 3 个阶段的能耗特点，然而该研究最后只专注于匀速段路径，忽略了路径的加速部分和减速部分。白盒建模研究相较于黑盒建模减少了比较多的重复性试验，从原理性分析可以更好地了解农用无人机任务飞行的能耗变化特点，且通过该方法建立的模型可以更好的规划无人机任务以降低能耗。国内外学者已从不同的角度对无人机作业能耗建模进行研究，提供了较高的参考价值。对于目前针对农用无人机飞行任务下不同飞行参数下的作业能耗缺乏有效评估的问题，本文在前人研究的基础上，通过白盒建模方法针对农用无人机动力系统各个能耗部件进行建模，并结合农用无人机任务中的飞行过程得到作业能耗模型，通过建立一个可靠且完整的农用无人机作业能耗模型，来评估飞行任务下的作业能耗。同时上述国内外研究所使用的无人机大多数轴距和尺寸较小，比较适用于航拍、巡检等简单任务，没有针对农用无人机的任务特点进行研究。因此，本研究针对农用无人机飞行任务中的速度、航程和载荷这 3 个动态参数设计试验，并进行试验验证，以验证模型的有效性。1白盒理论模型构建理论模型主要是由电池输入功率和螺旋桨输出功率的总效率模型以及克服重力所做的功、加速减速段所做的功和克服空气阻力所做的功的做功模型组成，其中总效率模型表示电池输入功率到螺旋桨输出功率的转化关系。如图 1 所示，首先根据农用无人机所使用的动力系统的螺旋桨、电机、电调和电池这 4 个组件的原理进行白盒建模16，得到总效率模型，然后根据所得到的无人机总效率模型结合无人机飞行过程的动态飞行参数，计算出无人机匀速直线飞行所需的能量，包括保持无人机高度所需能88农业机械学报2 0 2 3 年量与克服空气阻力所需能量的总和，最终得到农用无人机的作业能耗模型。图 1总效率模型示意图Fig 1Schematic of total efficiency model当前作业能耗模型的输入参数为动力系统参数和飞行参数。其中动力系统参数是来自所使用无人机供应商提供的参数，飞行参数是来自飞行任务时所设定的参数，阻力参数中无人机受空气阻力影响的有效面积近似为无人机最大截面面积17，阻力系数为无人机整体的阻力系数18。1.1总效率模型构建为了得到总效率模型，需要对农用无人机所使用的动力系统螺旋桨、电机、电调和电池这 4 个组件进行建模。在建模时首先需要计算螺旋桨输出功率，再根据螺旋桨输出功率对应的螺旋桨转速和转矩，求出电池输入功率。1.1.1螺旋桨输出功率农用无人机在不同任务下有不同的高度，处于不同环境也存在不同的温度，不同高度和温度下的大气压强和空气密度不同，这两个因素都直接影响无人机螺旋桨所提供的推力，进而对无人机的能耗存在影响。将无人机当前高度 h 和温度 Tt输入，求解出当前大气压强 pa19 和空气密度 为pa(=101 3251 0.006 5h273+T)t5.256 1(1)=273pa101 325(273+Tt)0(2)式中0 标准空气密度农用无人机通常采用定直径和定螺距的螺旋桨，螺旋桨转速 N20 为N=60TD4PCT槡(3)其中T=m0gnp(4)CT=1432BPK0p 0A+K0(5)式中T 单个螺旋桨所提供的拉力DP 螺旋桨直径CT 螺旋桨拉力系数m0 整机质量g 重力加速度np 螺旋桨个数 叶片翼型面积修正系数 位置系数BP 桨叶数量K0 升力曲线斜率 下洗气流校正因子p 螺旋桨桨叶角0 螺旋桨零升力角A 展弦比螺旋桨转矩是指螺旋桨运转时克服空气阻力所需的转矩，则螺旋桨转矩 M21 为M=CQ(N)602D5P(6)其中CQ=18A2Cd2B2P(7)Cd=Cfd+AK20e(p 0)2(A+K0)2(8)式中CQ 螺旋桨转矩系数，是螺旋桨转矩的无纲量系数Cd 螺旋桨叶片截面翼型的阻力系数Cfd 零升阻力系数e 奥斯瓦尔德效率因子螺旋桨桨叶角定义为中螺旋桨的旋转平面(A B)与螺旋桨桨叶弦线(A C)之间的夹角。假设其沿其半径方向保持不变，则桨叶角 p为p=arctanHPDP(9)式中HP 螺旋桨螺距将上面建模得到的螺旋桨拉力特性和转矩特性与常见二叶螺旋桨试验数据进行比较，综合考虑不同螺旋桨的型号和工艺，选取一组适应性较广的平均参数22，则取值为:A=5，=0.85，=0.75，=0.5，e=0.83，Cfd=0.015，0=0，K0=6.11。综合以上公式可计算螺旋桨输出功率 Pp为Pp=30npMN(10)1.1.2电池输入功率在飞行中，螺旋桨和电机可以视为一个整体，所以螺旋桨转矩等价于电机负载转矩、螺旋桨转速等价于电机空载转速，而且农用无人机一般使用无刷直流电机，无刷电机可以等效为一个永磁直流电机模型23，结合电机基本参数可得出电机等效电压Um和电机等效电流 Im24 为Um(=MKT+Im)0m+KEN(11)Im=MKT+Im0(12)式中KT 电机转矩常数98第 8 期李继宇 等:农用多旋翼无人机作业能耗白盒模型构建与试验KE 电机反电势常数Im0 电机空载电流m 电机内阻电机转矩常数 KT和电机反电势常数 KE存在关系式24，为KT=9.55KE(13)其中KE=Um0 Im0mKV0Um0(14)式中Um0 电机空载电压KV0 电机空载 KV 值电调调制下无刷电机的转速范围主要取决于电机负载转矩与电池电压，所以电调调节占空比后的等效直流电压 Ue为Ue=Um+Ime(15)式中e 电调内阻电调输出电压占空比可以等效反映油门大小，所以电调输出电压占空比 和电调输入电流Ie为=UeUb(16)Ie=Im(17)式中Ub 电池电压电池除了提供动力系统的电流，还提供其他机载设备(如飞行控制器、机载计算机等)电流，同时电池电流在传输过程会出现电流损耗，根据测算这部分所消耗的电流 Io一般约为 1 A。所以电池电流Ib为Ib=nrIe+Io(18)综合以上公式得电池输入功率 Pb为Pb=UbIb(19)1.1.3总效率模型最终根据螺旋桨输出功率 Pp和电池输入功率Pb得到总效率 total为total=PpPb=30nrMNUbIb(20)1.2作业能耗模型构建为了构建最终的作业能耗模型，对农用无人机的整个飞行过程的作业能耗参数进行研究。其中作业能耗参数中无人机的载荷主要影响飞行时克服重力所做的功，速度和无人机航程主要影响整段飞行任务的作业能耗。建模时结合总效率模型后，再对飞行任务中速度、航程和载荷的这 3 个动态参数进行分析，得到最终的作业能耗模型。如图 2 所示，无人机飞行时克服重力所做的功为 EG，在带载荷时还会受到额外的重力影响。再根据动能定理，得到加速减速段所做的功为 EK，最后结合阻力所做的功为ED，经过总效率模型转换得到所需作业能耗 E。图 2作业能耗示意图Fig 2Schematic of dynamic energy consumption1.2.1无人机克服重力所做的功 EG根据激励盘理论，螺旋桨可以简化为一个无厚度可穿透的圆盘25，所有螺旋桨圆盘的总面积At为At=nr(DP)22(21)在飞行时，空气流过螺旋桨，向下推动空气，螺旋桨每秒向下推动的空气质量 m为ma=Atva(22)其中va(=2m0gA)t12(23)式中va 流过螺旋桨圆盘的空气速度所有螺旋桨对空气质量 ma施加的力 FP等于无人机的重力，FP为FP=mava=m0g(24)最后可得功率 P026 为P0=Atv3a2(=2A)t12(m0g)32(25)图 3加速减速段所做功示意图Fig 3Schematic of work done by accelerationand deceleration sections农用无人机沿着直线飞行，如图 3 所示。在 t0时刻，从 v0开始加速直到在 t1达到速度 v1。保持速度直到 t2，在 t3减速到 v0，假设速度 v1和 v2相同，速度都等于 v，无人机克服重力所做的功 EG为EG(=dv+v)aP0(26)式中d 飞行距离v 飞行速度a 加速度1.2.2加速减速段所做的功 EK根据动能定理，如图3 所示，假设速度 v1和 v2相同，且等于设定速度 v，则加速减速段所做的功相09农业机械学报2 0 2 3 年同，EK为EK=m0v2(27)1.2.3克服空气阻力所作的功 ED在飞行过程中，无人机受到空气阻力影响，无人机所受阻力 FD27 为FD=12CDAev2(28)式中Ae 无人机受空气阻力影响的有效面积CD 阻力系数农用无人机沿着直线飞行，会持续受到空气阻力，则克服空气阻力所做的功 ED为ED=12dCDAev2(29)1.2.4农用无人机作业能耗根据叠加原理，综合克服重力所作的功为 EG、加速减速段所做的功为 EK、阻力所做的功为 ED，再经过总效率模型转换得到所需要提供的动态能量，得到最终的作业能耗模型 E 为E=30nrMNUbI(bdv+v)aP0+m0v2+12dCDAev2(30)2模型验证设计不同的试验来验证所构建作业能耗模型的有效性。对于农用无人机飞行任务下的应用场景，农用无人机应用在不同任务下需要携带不同的设备即挂载质量不同的载荷28，如测绘需要挂载相机，遥感需要挂载高光谱成像仪，施药需要挂载药箱等，对于不同任务下需要设定不同的速度，不同地块也需要设计不同长度的航程。因此针对农用无人机飞行任务中速度、航程和载荷的这 3个动态参数，设计试验对比不同速度、航程和载荷下机载采集设备采集得到的作业能耗数据与相同条件下理论模型计算的作业能耗数据，以此来验证模型的有效性。2.1试验场地试验在华南农业大学华山运动场和华南农业大学教学科研基地进行。试验期间使用 Kestrel 的NK5500 型气象站记录风速和风向，如图 4 所示。试验时风速基本集中在 3.3 m/s 以下，根据风力等级划分，风速 3.3 m/s 以下为轻风，气象环境较为稳定，符合试验要求29。2.2材料和设备为了验证所得到的作业能耗模型，所有试验主要使用自行组装的专用型无人机 A(型号为飞越图 4风速及风向Fig 4Wind speed and directionX6，华南农业大学无人机创新实验室)，最大挂载载荷为 2.5 kg，可挂载摄像机、高光谱相机或播种装置等农用设备。目前该机型设计最大前飞速度为12 m/s，单程最长飞行距离为 6 km。可满足试验的需求。本次试验还使用了通用型无人机 B(型号为大疆 S1000，大疆创新科技)，为 2014 年发布的通用消费级航拍无人机，而该无人机为当年的通用级载机，因其稳定安全的飞行性能，广泛应用于各种无人机应用领域，如图 5 所示。两架无人机的动力系统详细的配置参数是厂商提供动力系统各个部件的理论参数，如表 1 所示。图 5试验无人机Fig 5Tested UAV其中机载采集设备采用机载计算机(型号为NVIDIA Jetson TX2，英伟达)，用来运行机器人操作系统(OS)，使用该系统外部获取功率计采集的电池电压和电流信息以及 TK 定位系统的定位数据。其中功率计是一款库仑计(型号为 95V/200A，战锂智能科技)，可以检测电池输出的电流、电压和负载功率。电压测量范围为 0 95 V，数据采样频率5 Hz。功率计所配的分流器为 200 A/75 mV。其中TK 定位系统采用全系全频点高精度 GNSS 板卡(型号为 K705，司南导航)，支持北斗全球信号的接收，配合地面基站支持 TK 定位。2.3试验方法针对农用无人机飞行任务中速度、航程和载荷这3 个动态参数来设计试验。在动态参数选择上，19第 8 期李继宇 等:农用多旋翼无人机作业能耗白盒模型构建与试验表 1试验农用无人机基本参数Tab1Basic parameters of experimental agricultural UAV无人机类型参数数值A螺旋桨参数电机参数电调参数电池参数阻力参数其他参数DP/in15HP/in5.5Bp2np6Im0/A0.7Um0/V25KV0/kV450m/0.091e/0.024Ub/V25CD0.317 1Ae/m20.142 993m0/kg6.2h/m50/(kg m3)1.293Tt/15B螺旋桨参数电机参数电调参数电池参数阻力参数其他参数DP/in15HP/in5.5Bp2np8Im0/A0.7Um0/V25KV0/kV400m/0.086e/0.016Ub/V25CD0.362 2Ae/m20.129 21m0/kg7.4h/m50/(kg m3)1.293Tt/25结合当前农用无人机的作业特点及调查了解，对于速度动态参数选择上，目前如大疆和极飞等农用无人机在携带载荷的情况下速度最快在 10 m/s 左右，本试验所用专用型无人机 A 最大前飞速度为12 m/s，所以在设计试验时将最快速度设定为 12 m/s，从4 m/s 开始，以 2 m/s 为梯度进行递增。对于载荷动态参数选择上，农用无人机在飞行时的能耗与产生的推力成正比13，在一定范围内，农用无人机所携带的载荷越大，能耗越高，同时考虑到飞行时的安全因素和易实现性，所以在设计试验时选用 0 1.5 kg进行试验，以 0.5 kg 为梯度进行递增。对于航程动态参数选择上，无人机在飞行任务时主要是对每段航程进行规划，所有的航程组成飞行任务，所以在设计试验时选用 150 m 和 300 m 进行试验。如图 6 所示，针对速度参数，在保证航程和载荷不变的情况下，使用无人机 A 和 B 分别进行速度组的 8 组试验;针对航程参数，在保证载荷不变的情况下，设计使用无人机 A 进行航程组的 10 组试验;针对载荷参数，在保证航程不变的情况下，设计使用无人机 A 进行载荷组的 12 组试验，得到对应参数的作业能耗数据。试验无人机的飞行参数如表 1 所示。在飞行前在地面站设置好参数，其中所有试验组的加速度均设置为1 m/s2，在飞行试验时的无人机会根据目标设定的加速度加速和减速飞行，根据对应的速度进行设置，飞行的航线根据对应航程进行设置，并上传到飞行控制器中，开始飞行任务。图 6试验方法Fig 6Test method3试验结果对比分析为了验证模型的有效性，在试验后将速度、航程和载荷组分别由机载采集设备采集得到的试验作业能耗数据进行处理后，结合表 1 试验农用无人机的基本参数，在相同条件下与使用工程计算软件Maple 计算理论模型得到的作业能耗进行对比，得到两者的误差。进一步对模型的各个参数进行分析，引入模型修正系数对模型的输入参数进行修正，得到实际试验的速度和加速度数据，再代入模型计算出修正后的理论模型作业能耗，并与修正前进行对比，得到修正后的结果。3.1理论计算作业能耗与试验作业能耗对比首先通过模型计算速度、航程和载荷组的作业能耗，根据表 1 使用工程计算软件 Maple 计算理论模型对应的作业能耗。为了求出速度、航程和载荷组的作业能耗，需要对机载采集设备采集的数据进29农业机械学报2 0 2 3 年行处理。机载采集设备运行的 OS 系统在试验时会采集北斗 TK 和对应能耗数据，并将其自动处理及保存，根据 OS 系统采集得到的数据计算得出每个组的作业能耗，如表 2 所示。表 2不同设定速度组在不同航程下的理论计算和试验作业能耗Tab2Theoretical calculation and experimental dynamic energy consumption of different set speed groupsunder different voyagesJ速度/(m s1)速度组总能耗航程组总能耗载荷组总能耗100 m理论(A)100 m实际(A)100 m理论(B)100 m实际(B)150 m理论150 m实际300 m理论300 m实际0.5 kg理论0.5 kg实际1.0 kg理论1.0 kg实际1.5 kg理论1.5 kg实际426 29628 76228 90131 61137 58139 93371 43574 00729 46230 76232 76434 22236 19938 834620 82122 80222 88424 20428 38929 13951 09348 34123 30924 80225 90425 95728 60428 861819 41521 23421 31021 71524 93025 64442 15439 24621 68323 43424 04925 57626 51027 2751019 14720 28521 04721 34423 82723 80737 86535 98521 38422 98523 71624 67526 14028 2851223 95522 89335 97633 994将理论计算的作业能耗和对应的试验作业能耗数据进行汇总，如图 7 所示，可以看出理论计算的作业能耗和实际试验得到的作业能耗的变化情况。将 3 组的理论计算和试验作业能耗进行对比，代入误差计算公式e=|x r|x100%(31)式中x 理论计算的作业能耗r 实际试验得到的作业能耗可以得到理论模型计算的作业能耗数据与实际试验得到的作业能耗数据的误差，如图 8 所示。将所得到的误差数据绘制为箱线图，如图9 所示，由 均 值 点 可 得，速 度 组 的 平 均 误 差 约 为6.582%，航程组的平均误差约为 4.264%，载荷组的平均误差约为 5.057%。由中位线可得，速度组的误差中位数约为 7.654%，航程组的误差中位数约为 4.7%。载荷组的误差中位数约为 5.399%，可以看出在当前输入的动态参数下，该模型计算的结果不管是在平均误差还是误差中位数，都与实际数据结果相近，证明了模型的有效性。再将该模型与基于电池性能研究的黑盒模型的误差4.3%进行对比7，虽然不及黑盒建模准确，但也表现出较高的准确度。图 7不同速度、航程、载荷时理论计算和试验作业能耗变化曲线Fig 7Theoretical calculations and experimentaloperational energy consumption for three groups图 8理论计算和试验作业能耗的误差Fig 8Errors of theoretical calculations and experimental operational energy consumption for three groups39第 8 期李继宇 等:农用多旋翼无人机作业能耗白盒模型构建与试验图 9理论计算和试验作业能耗误差箱线图Fig 9Error box plots of theoretical calculations andexperimental operational energy consumption for three groups3.2模型误差分析及修正方法3.2.1模型误差分析模型的误差分析主要是对模型的计算结果和主要参数进行进一步分析，公式(30)中总效率主要是从电池所提供电能的输入功率到螺旋桨所消耗机械能的输出功率转化关系，影响整个模型，总效率模型是由厂商提供的动力系统各个部件的理论参数计算的结果。公式(30)中 P0主要影响无人机克服重力所做的功 EG，也是由厂商提供的参数计算而来。这两个参数在模型中为固定值。其次是对公式(30)中的阻力系数和无人机受空气阻力影响的有效面积进行分析，主要受环境风速和飞行速度影响，其中无人机受空气阻力影响的有效面积是由其俯视面的面积来估算，阻力系数为对应面积的阻力系数，以图 10 中无人机 B的阻力能耗占作业能耗的百分比为例，这两个参数影响的作业能耗只占很小一部分，且试验时的环境风速如图 4 所示，大多数是轻风等级以下风速，轻风对无人机所受空气阻力影响很小，主要还是受到飞行速度影响，随着速度的增大，无人机的阻力所做功占作业能耗的比例逐渐增大。图 10无人机 B 的阻力所做功与作业能耗对比Fig 10Comparison of resistance work and dynamicenergy consumption of UAV最后是对模型的速度和加速度参数进行分析，其中速度参数影响整个作业能耗，加速度影响无人机克服重力所做的功，这两个参数在理想模型中为定值。为了研究这两个参数对模型的影响，首先对环境风进行分析，在试验时环境风基本处于 3.2 m/s，而在环境风速 5 m/s 下，在设定速度最高 10 m/s 的飞行速度下，环境风对无人机迎风倾角的影响可以忽略不计，而无人机倾角影响其飞行速度，所以在当前试验环境下，对加速度和速度的影响也可以忽略不计29。其次对机载采集设备采集的数据做进一步分析，在当前飞行控制器控制下，无人机在试验时由于本身存在控制误差，所以不能达到所需的速度。对于飞行控制器而言，不同控制策略会导致无人机不同的能量消耗30，进而影响无人机飞行的作业能耗。对于飞行控制器而言，影响最大就是速度和加速度。原因有:当前飞行控制器下，无人机在飞行时并不能达到预定的速度，而且设定速度越大时，其与实际的速度的误差越大，以无人机 B 的试验数据为例进行分析，如图 11 所示，匀速段的设定平均速度为 4、6、8 m/s，但实际所达到的平均速度只有3.840 8、5.788 5、7.699 9 m/s，而对于设定平均速度为 10 m/s 的情况，实际的试验速度最高只达到约9.175 5 m/s。在当前飞行控制器下，无人机在加速和减速航段中，也不能达到设定的加速度，同样以无人机 B 的试验数据为例进行分析，得到设定平均速度为 10 m/s 的 4 组试验加速减速段的平均加速度为0.8679、0.8534、0.8711、0.8363 m/s2，低于设定的平均加速段 1 m/s2。图 11无人机 B 的试验速度变化曲线Fig 11Test speed of UAV综合以上对模型各个参数的分析，模型的总效率和功率由厂商提供的参数计算得出，阻力系数和无人机受空气阻力影响的有效面积对模型结果影响较小，且主要是受速度参数影响，模型的速度和加速度参数对模型的结果影响较大，将主要对这两个参数进行修正。3.2.2模型修正的方法及修正后结果对比综合 3.2.1 节对模型各个参数进行分析的结49农业机械学报2 0 2 3 年论，对模型的速度和加速度参数进行修正，模型修正方法主要是引入修正系数来对模型的误差进行修正。然后针对速度和加速度输入参数引入速度和加速度修正系数来修正模型，对于飞行控制器，需要试验修正系数来修正该模型的输入参数，使得模型的输入参数更加接近当前飞行控制器的实际控制效果。E=30nrMNUbI(bdv+v)aPG+m0(v)2+12dCDAe(v)2(32)式中 速度修正系数 加速度修正系数这两个修正系数需要对当前使用的飞行控制器输出进行分析，得到实际的速度和加速度，进而得到对应的修正系数。通过修正系数修正后的模型就可以得出更加准确的模型计算值和实际值之间的误差，以此来验证模型。这两个修正系数的获取方法为:对于当前飞行控制器，根据采集得到北斗 TK 数据的经纬度数据和时间数据，进而计算得出航程和时间，最后通过速度和加速度公式求出每一个匀速段的实际平均速度和实际平均加速度，得到两个输入参数的修正系数，如表 3 所示。表 3所有组在不同航程下的实际平均速度、加速度及对应修正系数Tab3Actual average speed，acceleration and corresponding correction factor of all groups under different voyages理想速度/(m s1)速度组航程组加速度组速度/(m s1)修正系数加速度/(m s2)修正系数速度/(m s1)修正系数加速度/(m s2)修正系数速度/(m s1)修正系数加速度/(m s2)修正系数43.8570.964 250.8830.8833.9010.975 250.8520.8523.8580.964 500.8810.88165.8280.971 330.8820.8825.8500.975 000.8910.8915.7770.962 630.8770.87787.7640.970 500.8900.8907.8150.976 880.8890.8897.7360.967 000.8660.866100.8560.8569.7460.974 600.8940.8940.8740.8741211.7170.976 420.9010.901使用修正系数对模型的输入参数进行修正后，得到最终的作业能耗模型误差，如图 12 所示。从均值点可以看出，修正后速度组、航程组和载荷组相比修正前的平均误差减小约 2.701、1.118、3.092 个百分点。从误差中位线可以看出，误差中位数减少约 3.267、1.617、3.612 个百分点。修正模型输入参数后的模型也避免了因为飞行控制器带来输入参数误差的问题，保证模型的输入量更加接近于图12所有组修正前后的模型计算和实际的作业能耗误差Fig 12Corrected model calculations and actual dynamicenergy consumption errors for all groups真实输入量，并且修正模型输入参数后的理论计算数据相比修正前更加接近于试验得到的数据。4结论(1)建立了一种针对农用无人机的作业能耗模型，并通过农用无人机的分组试验得到的数据进行了验证，结果表明，验证试验中的最大平均误差约为6.582%，最大误差中位数约为 7.654%，模型具有较高的精度。(2)通过对模型参数的误差分析和修正，分析了不同参数对当前模型的影响，使用模型修正系数来对速度和加速度的误差进行修正，修正后的模型精度相比修正前的最大平均误差减少约 3.092 个百分点，误差中位数减少约 3.612 个百分点。(3)使用专用型无人机 A 和通用无人机 B，针对速度、载荷和航程这 3 个动态参数进行试验，表明该模型可扩展到不同任务，如授粉播种、农情监测等任务。对飞行控制器的输入参数进行修正的方法，可以使模型适用于其他飞行控制器，进而扩展到其他旋翼机型。参考文献 1 高志政，彭孝东，林耿纯，等 无人机撒播