梯度风场滑翔仿生飞行器设计与仿真

摘  要：      将信天翁从梯度风场中获取能量滑翔数千英里的飞行策略应用在飞行器上是突破飞行器续航瓶颈的重要方法之一。 本文结合信天翁身体特性和飞行策略， 设计了一种具有信天翁仿生外形、 气动布局和控制方式的仿信天翁飞行器， 巡航时升阻比为梯形翼飞行器的1.13倍， 最大可达梯形翼飞行器升阻比的1.31倍， 通过翼尖变形阻力最高可减小12%， 翼根载荷减小35%。 根据信天翁在滑翔时不时扑动翅膀的行为提出了间歇动力梯度风滑翔方法， 对仿信天翁飞行器进行了风梯度滑翔操控性能仿真。 结果表明实际飞行路径与风梯度滑翔路径的误差小于5%， 能量消耗为匀速飞行的51.5%， 验证了间歇动力梯度风滑翔策略的合理性。

关键词：     梯度风场； 飞行策略； 仿生无人机； 变形机翼； 操控性能

中图分类号：       TJ760； V221

文献标识码：    A

文章编号：     1673-5048（2024）05-0103-07

DOI： 10.12132/ISSN.1673-5048.2023.0248

0  引  言

在自然界中， 许多大型海鸟可以利用海面上垂直分布的水平梯度风场进行远距离滑翔飞行， 如信天翁可以飞行数千英里并且几乎不拍动翅膀。 所以， 将信天翁从梯度风场中获取能量的飞行策略应用到飞行器上， 可以改善飞行器的巡航时间甚至做到无动力飞行。

目前， 国内的无人机风梯度滑翔策略研究较少， 且主要为滑翔机理、 路径优化等理论分析， 对操控方面的研究很有必要。 2001年， Kicenuik提出滑翔机在水平层流和垂直层流中的最佳操控方式［1-2］。 2011年， Lawrance等总结了信天翁远距离滑翔飞行的策略和方法， 并分析其应用于飞行器动态滑翔的可行性［3-6］。 2012年，  Sachs等研究了信天翁的飞行轨迹， 公布了一些关于信天翁风梯度滑翔的实验数据［7］。 2014～2017年， 高显忠等对梯度风获取能量的运动特征进行了分析， 建立了自主梯度风动态滑翔飞行的简化方法［8-10］。 2018年， Stempeck等根据信天翁等多种鸟类翼翅设计了平面机翼并与标准几何翼形进行对比［11］。 2021年， Mathew等模拟了信天翁和猎鹰的空气动力学特性， 可以作为进一步设计具有最佳空气动力学特性和低能耗的仿生无人机的基础［12］。 同年，  An等建立了信天翁仿生机翼， 并分析了该机翼的气动性能［13］。 2022年， Wang等提出一种动态滑翔轨迹， 为深入研究动态滑翔机理提供飞行试验基础［14］。 以上研究表明， 信天翁的翼翅结构和滑翔策略可以极大地延长飞行器的巡航时间， 但目前的研究所用的飞行器属于传统飞行器， 没有适用于风梯度滑翔的特殊飞行器。 其次， 目前的研究主要集中在运动学分析和对航迹进行分析优化上， 虽然这表明风梯度滑翔可以延长飞行器巡航时间， 但缺少操控性方面的研究。

本文从信天翁的外形中提取了仿生机翼和机身的几何参数， 建立了仿生飞行器的三维模型， 并对飞行器进行气动计算和分析； 根据信天翁姿态控制方式设计变形机翼来控制飞行器的俯仰、 滚转和偏航运动并且实现减阻减载的效果； 建立虚拟海上梯度风场， 对仿信天翁飞

行器从海面上梯度风场中获取能量进行滑翔飞行进行操控仿真分析， 验证了仿生飞行器间歇动力风梯度滑翔策略。

1  飞行器总体设计

1.1  机翼仿生外形设计

信天翁由于性别、 种群、 年龄的差异， 身体参数也有

所不同。 本文选取一种中等体型的信天翁——黑背信天

收稿日期： 2023-12-26

作者简介： 刘建河（1976-）， 男， 河南叶县人， 副研究员， 硕士生导师。

*通信作者： 李加东（1980-）， 男， 山东临沂人， 博士， 研究员， 博士生导师。

翁为仿生对象， 据密歇根大学动物博物馆介绍， 该信天翁身长为80 cm， 体重为2 210～2 800 g， 是一种有极长翼羽的大鸟， 翼展为195～215 cm。 本文以黑背信天翁直线滑翔状态为拟合对象， 建立翼展为2 m、 平均弦长为218 mm、 根弦长为260 mm的仿生机翼模型（总体参数指标如表1所示）， 并对其进行优化。

航空兵器  2024年第31卷第5期

刘建河， 等： 梯度风场滑翔仿生飞行器设计与仿真

机翼制作通常需要知道展长、 弦长、 后掠角和二面角， 为了方便建模， 本文用当前截面翼型的前缘点相对翼根翼型的前缘点的前后偏移表达后掠角， 用上下偏移表达二面角。 为了收集信天翁翼翅的参数， 通过大量信天翁水平滑翔时的腹部仰视图和尾部视图进行对比测量［10-11］。 首先， 根据信天翁滑翔时翼翅的后掠角、 二面角对其进行分段， 如图1所示， 翼一具有较大的上反角和后掠角， 翼二上反角逐渐减小为零且机翼前掠， 翼三机翼下反后掠。 然后， 以展长为比例因子， 计算每个翼段断面的弦长、 偏移量（上下偏移和前后偏移）， 结果记录在表2中。 表中前后偏移和上下偏移指当前截面翼型的前缘点距翼根翼型前缘点的竖直和水平距离。

1.2  翼型确定及优化

信天翁翼翅的尺寸参数确定后， 需要选择合适的翼型建立机翼三维模型。 翼型设计软件中的GOE 174是基于信天翁翼翅截面创建出的翼型， 具有高升阻比的特点， 然而信天翁在飞行中几乎看不到尾巴， 在实际应用中仿照鸟类的无尾布局要求机翼自身具有较强的俯仰稳定性， 所以需要对GOE 174翼型进行俯仰力矩系数方面的优化。

本节使用遗传算法结合有限元分析软件对翼型进行优化计算。 使用Hicks-Henne型函数对翼型上下表面各关键点形状的影响平滑均衡， 根据优化软件提供的邻域培植遗传算法NCGA对翼型进行优化［15］。

首先使用Hicks-Henne型函数对翼型上下曲线进行调整， 然后对翼型进行网格划分， 结果如图2所示。 翼型弦长为1 m， 网格数目为86 366， 最大网格质量为0.427。 解算器设置： 采用标准k-ε模型， 经验常数取默认值， 二阶迎风格式控制方程， 边界条件为速度入口， 迎角为4°， 来流速度为15 m/s， 求出4°迎角时的升阻比和俯仰力矩系数。 多目标优化器优化方法参数为每代个体6， 共8代， 其余参数采用默认值。 设计要求如下：

（1） 4°迎角， 标准SST k-omega模型进行有限元计算。

（2） 约束条件为： 优化后的翼型最大厚度为8.6%～9.8%， 俯仰力矩系数Cm范围为-0.3～0， 4°迎角时的升阻比k范围为65～85。

（3） 要求减小俯仰力矩系数， 同时尽可能保持较大的升阻比。

优化前后翼型对比如图3所示。 可以看出优化后的翼型在俯仰力矩系数方面的表现很好， 保持在0附近， 并且在大迎角情况下也可以稳定在0左右。

1.3  飞行器总体布局

根据优化后的翼型和表2参数， 建立仿生机翼模型， 再根据信天翁身体建立机身模型， 信天翁滑翔时为飞翼布局， 所以本文所设计的飞行器为尾推式飞翼布局， 如图4所示。

2  空气动力学分析

2.1  飞行器性能分析

使用仿真软件对二维翼型分析或者对整个飞机进行气动分析， 基于升力线理论、 涡格法、 三维面元法对机翼进行求解。 对整机进行定速分析， 机翼参数为： 投影面积为0.4 m2， 展弦比为10， 平均气动弦长为218 mm， 翼展为2 m。 速度为15 m/s， 机翼安装角为4°， 尾翼翼型为NACA0012， 由信天翁翼翅在身体上的分布情况确定尾翼前缘与机翼前缘的距离为290 mm， 平尾面积为197 cm2， 平尾力臂为262.5 mm， 重心位置到机翼前缘距离为72.5 mm， 位于机翼三分之一处且在气动中心之前。 仿生飞行器网格划分如图5所示。

飞行器定速分析结果如图6所示。 可以看出， 当俯仰力矩系数为0时， 仿生飞行器迎角为0.4°， 表明飞行器稳定飞行时的迎角为0.4°， 机翼自身迎角等于安装角+飞行器迎角=4.4°， 此时飞行器升力系数0.416， 根据升力公式， 此时机翼提供的升力为22 N， 阻力系数为0.017 2。 仿生飞行器在正迎角下升阻比大于矩形翼飞行器， 0°迎角巡航时升阻比为梯形翼飞行器的1.13倍， 迎角越大仿生飞行器的升阻比相对梯形翼飞行器的倍数越大， 5°时为梯形翼飞行器升阻比的1.31倍， 所以仿生飞行器在爬升时具有较大优势。

稳定性分析结果如图7和表3所示。 从表3数据中可以看出， 纵向模态频率ω12=8.427 92+11.681 62=14.4， ω34=0.55 ； 阻尼比ξ12=0.59， ξ34=0.018； 周期T12=0.44， T34=11.4。 前四个模态是稳定的， 而螺旋模态为不稳定状态， 这是由于飞行器没有方向舵， 无法控制飞行器的偏航姿态。 针对这种情况， 模仿信天翁飞行时翼翅上下折叠控制自身姿态的行为设计了转折变形机翼， 此模态可以通过控制翼尖的折叠运动进行控制。

2.2  翼尖设计

信天翁使用翼翅就可以控制滚转、 俯仰和偏航运动。 图8（a）为机翼初始位置， 在飞行器控制方面以信天翁为蓝本设计了转折翼尖，  如图8（b）所示，  俯仰和滚转运动

均由翼尖进行控制， 增大了舵面积； 另一方面信天翁在滑翔时需要面临不断变化的阵风冲击， 这对机翼的强度和抗疲劳性能具有一定的要求。 为了更好地应对阵风冲击［16］并延长机翼的使用寿命， 本文设计了一种具有转折翼尖机构的变形机翼， 如图8（c）所示。

转折翼尖的展向长度直接影响飞行器的偏航运动， 对不同翼尖展长进行偏航力矩分析， 翼尖转动角度均为20°， 来流速度15 m/s， 得出的结果如图9所示。 当翼尖展长为480 mm时产生的偏航力矩最大为0.43 N·m， 所以翼尖展长选为480 mm。

对两侧翼尖做对称折叠运动进行仿真。 翼尖初始位置如图8（a）所示， 水平向下8°即-8°， 设置翼尖角度每0.5 s增加4°， 结果如图10～11所示。 在翼尖从下反到上反的运动过程中， 阻力最大减少12%， 翼根载荷减少35%。 当翼尖位于初始位置时阻力最大， 翼尖水平时翼根载荷达到最大。

3  梯度风场飞行仿真

信天翁滑翔过程中大部分时间不用拍动翅膀， 只有

在自身动能不足时会扑动翅膀。 基于这种飞行方式， 本文提出了间歇动力梯度风滑翔策略， 即在每个滑翔周期初始阶段进行逆风爬升时， 飞行器继承自上一周期的能量不足以使飞行器爬升到指定高度的情况下需要启动螺旋桨对飞行器进行能量补充， 当飞行器速度达到周期滑翔的标准时再关闭螺旋桨， 一般会在逆风爬升阶段提供动力， 到顶点转弯阶段关闭动力。 下面通过操控仿真对该方法进行验证。

信天翁进行梯度风滑翔的必要条件是该区域存在梯度风场， 为了方便计算， 假设梯度风场为沿水平方向垂直分布的风速快速变化的风场， 风速大小仅与高度相关。 文中梯度风场使用美军标准MIL-F-8785C［17］， 风速与高度的关系可表示为

W=W0ln（z/z0）ln（z1/z0） （1）