无人潜航器无动力下潜轨迹仿真分析
- 格式:pdf
- 大小:1.46 MB
- 文档页数:3
国外水下无人潜航器及其通信技术发展综述一、本文概述随着科技的快速发展,无人潜航器(Unmanned Underwater Vehicle,UUV)作为海洋探索与利用的重要工具,已经吸引了全球科研人员和工程师的广泛关注。
作为无人潜航器的重要组成部分,水下通信技术对于实现潜航器的远程控制、数据传输、多潜航器协同作业等功能具有关键作用。
本文旨在综述国外水下无人潜航器及其通信技术的发展现状与趋势,分析当前主流通信技术的优缺点,并探讨未来可能的研究方向和应用前景。
通过对国外相关文献的梳理和分析,本文旨在为国内外从事水下无人潜航器及通信技术研究的学者和工程师提供有益的参考和启示。
二、国外AUV的发展现状近年来,随着科技的飞速发展,国外在自主水下航行器(AUV)领域取得了显著的进步。
AUV作为水下无人潜航器的一种,其自主导航、环境感知、任务执行等能力不断增强,为海洋科学研究、海底资源勘探、水下搜救等领域提供了有力支持。
在硬件设计方面,国外的AUV技术日趋成熟。
许多先进的AUV已经实现了小型化、模块化、高度集成化,以适应不同复杂度的水下环境。
例如,某些AUV采用了先进的推进系统,包括矢量喷水推进器、机械式螺旋桨等,以提高其机动性和稳定性。
同时,为了应对深海高压、低温等极端环境,AUV的耐压壳体和材料技术也在不断更新,确保了AUV的安全性和可靠性。
在软件技术方面,国外的AUV已经实现了高度智能化和自主化。
通过集成先进的算法和人工智能技术,AUV可以自主完成路径规划、环境感知、目标识别等任务。
随着深度学习技术的发展,AUV在图像识别、声呐信号处理等方面也取得了显著突破,进一步提升了其在水下复杂环境中的作业能力。
在通信技术方面,国外的AUV同样取得了长足的进步。
为了实现在水下环境中的数据传输和远程控制,研究人员开发了一系列高效、稳定的水下通信技术。
例如,某些AUV采用了高速水声通信技术,实现了与水面基站或卫星的实时数据传输;还有研究团队在探索利用电磁波或光学通信技术在水下环境中实现数据传输的可能性。
海洋机器人与无人潜水器的比较分析在探索海洋的进程中,海洋机器人和无人潜水器逐渐成为了重要的工具。
它们各自具有独特的特点和优势,为我们了解和利用海洋资源提供了有力的支持。
接下来,让我们对海洋机器人和无人潜水器进行一个详细的比较分析。
首先,从定义和功能上来看,海洋机器人是一种能够在海洋环境中自主或半自主执行任务的智能设备。
它的功能十分广泛,可以进行海洋观测、数据采集、海底地形测绘、海洋资源勘探等工作。
而无人潜水器则通常是指能够在水下自主航行和作业的潜水设备,其重点在于对特定区域或目标进行深入的探测和研究。
在结构设计方面,海洋机器人的形态多样。
有的像小型船只,有的则类似水下飞行器。
其结构设计需要考虑到海洋环境的复杂性,如水流、水压、腐蚀等因素。
同时,为了实现各种功能,还需要配备先进的传感器、通信设备和动力系统。
相比之下,无人潜水器的结构相对较为紧凑,通常以流线型为主,以减少水阻,提高航行效率。
其外壳材料要求更高的耐压性能,以适应不同深度的水下环境。
动力系统是两者的一个关键区别。
海洋机器人的动力来源多样,包括电能、燃油、太阳能等。
一些海洋机器人依靠电池供电,适合短时间、近距离的任务;而一些大型海洋机器人则可能采用燃油发动机,具备更长的续航能力和更强的动力。
无人潜水器的动力系统则通常以电池为主,因为其在水下作业时间相对较短,对动力的需求相对较小。
但对于需要在深海长时间工作的无人潜水器,也会采用一些特殊的能源供应方式,如热液能源等。
在操控方式上,海洋机器人通常可以通过卫星通信或岸基基站进行远程控制,操作人员可以在陆地上对其进行实时监控和指令下达。
部分先进的海洋机器人还具备一定的自主决策能力,能够根据预设的程序和环境变化自主调整行动方案。
无人潜水器的操控则相对较为复杂,由于其在水下工作,通信受到很大限制,往往需要提前规划好任务路径和动作,在水下依靠自身的控制系统完成任务。
再来看它们的应用领域。
海洋机器人广泛应用于海洋科学研究、海洋环境监测、海上石油和天然气勘探等领域。
第26卷第2期 水下无人系统学报 Vol.26No.22018年4月JOURNAL OF UNMANNED UNDERSEA SYSTEMS Apr. 2018收稿日期: 2017-08-05; 修回日期: 2018-01-22.基金项目: 发改委海洋工程装备研发及产业化项目(20151409).作者简介: 黄苏和(1987-), 男, 工程师, 主要从事水中兵器系统研究.[引用格式] 黄苏和, 王凯帅, 刘星. 潜用AUV 自航发射弹道建模与仿真[J]. 水下无人系统学报, 2018, 26(2): 129-132.潜用AUV 自航发射弹道建模与仿真黄苏和, 王凯帅, 刘 星(中国船舶工业系统工程研究院, 北京, 100094)摘 要: 采用鱼雷发射管自航发射自主式水下航行器(AUV)是潜艇发射AUV 的最理想方式。
为研究潜用AUV 自航发射技术, 文章基于后端有效补水的发射装置, 建立了AUV 自航发射弹道模型, 包括AUV 发射管内弹道模型和管口弹道模型, 对AUV 自航发射出管的可行性进行了仿真验证。
仿真结果表明, 自航发射时艇速、AUV 出管速度以及AUV 自身正浮力或负浮力是影响AUV 自航出管可行性的主要因素。
该结果可为采用鱼雷发射管自航发射AUV 技术提供参考。
关键词: 自主式水下航行器(AUV); 潜艇; 自航发射; 弹道模型中图分类号: U674.941; TJ635; TP391.9 文献标识码: A 文章编号: 2096-3920(2018)02-0129-04DOI: 10.11993/j.issn.2096-3920.2018.02.005Modeling and Simulation of Trajectory for AUV Swim-outLaunch on SubmarineHUANG Su-he , WANG Kai-shuai , LIU Xing(Systems Engineering Research Institute, Beijing 100094, China)Abstract: Using torpedo launch tube for swim-out launch is the best method for submarine to launch an autonomous undersea vehicle(AUV). This study establishes the mathematical trajectory models of the AUV swim-out launch based on the launcher with backend water replenishment. The models include a submodel of internal trajectory of launch tube and a submodel of trajectory at the tube orifice. The feasibility of the AUV swim-out launch out of tube is verified by simulation, and the results show that the speed of a submarine, the launching speed, and the positive or negative buoy-ancy of the AUV are the main factors influencing the AUV swim-out launch. This study may provide a reference for swim-out launch of AUV through torpedo launch tube.Keywords: autonomous undersea vehicle(AUV); submarine; swim-out launch; trajectory model0 引言自主式水下航行器(autonomous undersea ve-hicle, AUV), 由于具有适应性强、远航程、大潜深、零伤亡以及隐蔽性好等一系列特点, 被誉为水下作战优势力量的倍增器[1-2]。
无人水面艇仿真系统设计与实现胡辛明;张鑫;钟雨轩;彭艳青;杨毅;姚骏峰【摘要】无人水面艇(unmanned surface vehicle,USV)的海上调试成本高、难度大、效率低.为了减少海上调试工作量,设计了仿真系统,规划了仿真流程,为无人水面艇控制系统的测试和分析提供了平台.该仿真系统包括5个模块:显控模块用于全局路径规划和综合信息监控;数据仿真模块模拟障碍检测传感器与位置和姿态传感器数据;障碍处理模块对检测数据进行预处理、栅格化、聚类和拟合,并对拟合后的动态障碍进行跟踪;导航和避障模块嵌入视线制导(line of sight,LOS)导航和椭圆聚类-碰撞锥推演的动态避障算法;运动控制模块采用广义预测控制-比例积分(generalized predictive control-proportion-integral-derivative,GPC-PID)串级控制方法.各个仿真模块具有独立性,可根据不同需求替换.仿真实验结果表明了该系统设计的正确性和合理性.【期刊名称】《上海大学学报(自然科学版)》【年(卷),期】2017(023)001【总页数】12页(P56-67)【关键词】无人水面艇;仿真系统;全局路径规划;导航和避障;运动控制【作者】胡辛明;张鑫;钟雨轩;彭艳青;杨毅;姚骏峰【作者单位】上海大学机电工程与自动化学院,上海200072;上海大学机电工程与自动化学院,上海200072;上海大学机电工程与自动化学院,上海200072;中国人民解放军理工大学理学院,南京210007;上海大学机电工程与自动化学院,上海200072;上海大学机电工程与自动化学院,上海200072【正文语种】中文【中图分类】TP242.3无人水面艇(unmanned surface vehicle,USV)平台具有体积小、速度快、机动灵活等特点,根据不同的搭载平台,无人水面艇具有不同的功能,如海事管理、海上搜救和水文环境监测等.然而,无人水面艇控制系统比较复杂,海上调试存在成本高、难度大和效率低等问题.为了减少海上调试时间,需要设计仿真系统,对无人水面艇控制系统进行测试,为海上调试做好充分的准备工作.目前,国内外学者已经对无人水面艇进行了比较深入的仿真研究.挪威科技大学的Fossen和Breivik等研究了无人水面艇的航迹跟随[1-3]、动态定位[4-5]、目标跟踪[6]等问题,阐述了基于视线制导(line of sight,LOS)导航方法;并利用CS2(cybership 2)船模型搭建了半实物仿真平台[7],对航线和艏向的数据进行分析,证明了导航算法的有效性.Larrazabal等[8]利用Nomoto模型建立了艇体的辨识模型,设计了基于遗传算法的增益比例积分控制最优化(gain scheduling proportional-integral-derivate control optimized by genetic algorithms,GS-PIDGAs)控制器,结合模糊逻辑控制器来解决动力学的不确定问题,仿真实验结果表明,航线和航向角偏差较小,满足要求.吉林大学的杨树仁[9]基于电子海图设计了单机非实时船舶运动控制仿真系统,该仿真系统由5个部分组成:航迹控制器、舵机执行机构、舵机反馈机构、船舶运动模拟器和电子海图,并对A,B,C类船进行了3级海浪干扰仿真实验,结果可以看出航线和航向偏差都满足航迹控制要求.哈尔滨工程大学的卢艳爽[10]设计的仿真系统包括全局路径规划、航迹跟踪和动态避障等,并规划了仿真流程.本工作基于A*启发式搜索算法(A*算法)实现全局路径规划,采用“跟随目标”的方法实现无人水面艇的航线无障碍跟踪,提出“远域协商”动态避障方法,并进行了动、静态障碍物避障仿真实验.已有的无人水面艇仿真主要针对单一特定功能,系统较为简单,而本工作从全局路径规划、传感器数据模拟、障碍检测和跟踪、运动控制、导航和避障等方面系统性地提出了相对完备的仿真方案,为无人水面艇控制系统提供测试和分析平台,可减少30%~40%的海上调试工作量,对缩短无人水面艇研发周期,提高海上调试效率有重要意义.本工作主要贡献如下:(1)提出了无人水面艇仿真系统设计方案;(2)基于生物激励神经网络的完全遍历路径规划方法,结合电子海图提供的障碍信息,实现了全局路径规划;(3)仿真罗经和GPS数据,根据Nomoto模型和广义预测控制-比例积分(generalized predictive control-proportion-integral-derivative,GPC-PID)串级航迹控制方法,提出了无人水面艇的位置和姿态更新方法;(4)仿真激光数据,根据连通域标记法、K-means算法思想实现障碍的检测,基于Kalman滤波算法实现动态障碍跟踪;(5)在无人水面艇仿真系统中加入了海浪干扰并进行了仿真实验,实验结果证明了本系统设计的合理性.通常,无人水面艇系统包括无人艇艇体、显控台、GPS、罗经、激光、4G雷达、前视声纳、运动控制计算机等[11],其中显控台与无人水面艇之间通过无线进行信号传递,显控台根据任务类型进行路径规划,并监控无人水面艇的综合状态信息;GPS和罗经采集位置和姿态信息;激光和4G雷达用于海面障碍检测;前视声纳用于海下障碍检测;运动控制计算机内部包含障碍处理算法、导航和避障算法、运动控制算法,结合航线、位置和姿态、障碍信息,计算并控制喷泵的输出和舵机的转动,从而控制无人水面艇的运动.由以上无人水面艇系统信息传递过程规划出仿真系统的流程(见图1).设置全局任务路径、障碍物和无人水面艇的初始位置和姿态,当激光检测到障碍物,经过障碍处理,将障碍信息传递给导航和避障模块;随后导航和避障模块根据路径、障碍信息以及位置和姿态信息计算处理,得到无人水面艇的航向和速度,而将航向和速度变化作为运动控制的输入,输出舵角和油门,舵角和油门用于位置和姿态的更新;位置更新后,判断无人水面艇是否到达路径终点,若没有,则在新的位置继续进行障碍检测,并将新的位置和姿态信息、障碍信息作为导航和避障的输入,更新数据,直到无人水面艇到达路径终点.由无人水面艇仿真流程可确定仿真系统设计方案.仿真系统软件分为5个模块:显控模块、数据仿真模块、障碍处理模块、导航和避障模块、运动控制模块,其中显控模块实现全局路径规划和监控,数据仿真模块仿真激光、罗经和GPS数据,导航和避障模块用于控制无人水面艇航向和速度,而运动控制模块则是根据航向和速度变化控制舵角和油门.2.1 显控模块显控模块是基于电子海图实现的.电子海图具有地理信息获取和显示、船舶速度航向监控、报警等功能.显控模块主要进行全局路径规划和无人水面艇的综合信息显示.全局路径规划分为任务区域覆盖扫测和多点规划,其中任务区域覆盖扫测是在给定任务扫测区域规划出绕开障碍物的最优扫测路径;多点规划是用于无人水面艇布放和回收时,在海图上人工选取路径点,自动生成绕开障碍物的最短航行路线.电子海图提供的静止障碍物包括海岸线、码头、浮标和岛屿等.任务区域的覆盖扫测是基于生物激励神经网络的完全遍历路径规划方法实现[12],其中生物激励神经网络用于移动机器人工作环境建模,而完全遍历路径规划包括模板模型法和A*算法,能够实现复杂岛礁区域周围环境的遍历.图2(a)为在电子海图上实现的存在岛礁的任务区域覆盖扫测,只要给定扫测范围、起点、方向和距离就可以自动规划出最优扫测路径;图2(b)为在海图上实现的多点规划(多点规划是基于A*算法实现的,能够实现代价最小的路径规划),只要给定路径点就可以规划出一条绕开障碍物的最短航行路径.2.2 数据仿真模块数据仿真模块包括仿真障碍检测、位置和姿态数据,其中障碍检测数据指的是激光扫到物体后得到的距离和角度,位置和姿态数据是由GPS和罗经采集得到的位置和姿态信息.2.2.1 激光数据仿真激光测距[13]的方法有多种,但目的都是为了获取目标的角度和距离.常用的方法是脉冲法,即在工作时向目标射出3维很细的激光(由光电元件接收目标反射的激光束),计时器测定激光束从发射到接收的时间差,从而得到距离.式中:D为激光传感器与被测障碍物之间的距离;c为光在真空中的传播速度(3×108m/s);∆T为激光两路脉冲信号之间的时间间隔;ς是激光传输介质的平均反射率,在空气中取1.通常,激光的扫描距离D、扫描范围ϕ、角度分辨率∆θ都为已知参数.将激光安装在无人水面艇上时,激光中垂面与艇的中垂面重合.在电子海图里以激光位置为原点,艏向为x轴,建立平面直角坐标系.设置圆形的障碍物位置和大小,Pn(x,y)为激光束打在障碍物上的点(见图3),则可以得到每一个激光束与x轴夹角θn,即获得了激光角度信息.联立激光束直线方程和圆的方程,可以得到激光束与圆的交点Pn(x,y)的坐标值,求得距离D,取其中较小值;若无交点,则距离D=0,即可得到激光的距离信息.依此类推,可以得到所有激光束数据,激光数据仿真即可完成.2.2.2 位置和姿态数据仿真无人水面艇的位置坐标和速度由GPS采集,姿态信息由罗经采集.姿态信息包括航向角、艏向角、俯仰角、横滚角、航向角速度、艏向角速度、俯仰角速度和横滚角速度等.由于与导航和避障相关的数据主要包括位置坐标、航向角、航向角速度和速度,故仿真系统针对这4个参数进行仿真.航向角速度由舵机运动决定,并受到艇体本身结构、惯性和风浪的影响.常用的舵角与航向角速度的辨识模型为经典的Nomoto模型[14],本工作采用简化的一阶Nomoto模型,即舵角变化δ和航向角速度r的关系为式中,S为传递函数,K为静态航向角速度增益,T为时间常数.航向角变化量∆ψ由航向角速度r积分得到,由此可以得到航向角ψt:式中,f为更新频率,ψt−1为上一帧航向值.无人水面艇的速度大小由油门控制,油门与速度可以用简单的线性关系表示,即式中,k为油门系数.无人水面艇的位置变化量∆P由速度V和航向ψ决定,由此可以得到位置坐标:式中,ψt−1为上一帧航向值.2.3 障碍处理模块激光数据处理模块包括障碍检测和跟踪.障碍检测是将激光采集的数据处理成椭圆障碍,其过程包括预处理、栅格化、聚类和拟合等.预处理是去除激光里面的特殊点,对于多线激光,还需将激光点投影到同一个平面;栅格化的方法比较多,本工作采用最大最小高度图法[15-16]对激光数据进行栅格化;采用通常用在二值图像处理的连通域标记算法[17-18]对栅格进行聚类,得到障碍物块;海上障碍多为船只,船只长和宽的比例较大,可以近似看成椭圆,基于K-means算法思想[19],对聚类的障碍物块的形状进行近似拟合,形成椭圆的障碍信息.激光数据处理过程如图4所示.图5为对实际场景中采集的激光数据进行处理的结果,其中(a)为多线激光检测实际场景,(b)为多线激光检测障碍物得到的立体点云,(c)为经多线激光数据处理形成的椭圆障碍信息.现实中的障碍物多为动态,对障碍物位置和速度进行预测必不可少.障碍目标跟踪有两个作用:一是减少或消除传感器的检测误差,二是为动态避障提供预测位置和速度信息.Kalman滤波器是一种线性递归滤波器[20],是用于状态最优估计的常用算法,通过动态的状态和观测方程描述系统,并采用递归滤波方法预测下一个状态的最优估计.由于仿真激光数据采样时间间隔短,处理得到的椭圆目标运动状态变化小,故将目标运动定为匀速直线运动.将运动椭圆的中心轨迹作为输入,利用最近邻域法将测量值与预测值进行关联,并采用Kalman滤波进行最优状态估计,从而得到下一时刻的椭圆障碍物位置和速度信息.图6为基于Kalman滤波算法的障碍目标跟踪过程.2.4 导航与避障模块在电子海图上利用全局路径规划方法生成任务路径,其中导航算法是控制无人水面艇在复杂的海洋环境干扰下,按照既定路径进行航行,故导航算法是无人水面艇执行任务的前提,导航算法的优劣直接影响到任务的执行结果.LOS是一种直线循迹的方法,该算法首先在航行路段上虚拟出一个追踪目标,然后引导船舶沿着船舶到追踪目标的视线方向航行.目前,常用的虚拟追踪目标确定方法主要包括封闭圈法和超前置位法.本仿真系统采用封闭圈法,即以船舶为圆心建立一个半径固定不变的封闭的圆圈,将封闭圈与计划航线的交点设为虚拟目标点.随着船舶不断接近计划航线,船舶航向偏差也随之减小,这样就达到同时减小位置和航向偏差的目的.无人水面艇在海上执行任务时会遇到其他船舶、岛礁、灯塔等障碍物,因此避障功能是其安全航行的重要保障.海上航行的船舶需要遵循《国际海上避碰规则》,已有很多基于《国际海上避碰规则》的导航避障算法,如神经网络算法、进化算法、模糊逻辑算法、2D栅格地图等.本工作中的避障算法遵循《国际海上避碰规则》,基于椭圆聚类-碰撞锥推演的无人水面艇海洋动态避障控制方法实现无人水面艇避障.该方法是基于Chakravarthy等[21]和Fiorini[22]的点与点、点与圆的碰撞锥方法,其本质是在未知轨迹的移动目标之间进行碰撞检测,根据已计算出的碰撞锥区间采取相应的避碰措施.无人水面艇采取避障措施后会偏离原先设定的路径,因此利用LOS导航算法可实现其航迹回归.2.5 运动控制模块导航和和避障算法的输出是期望速度和航向(速度和航向的变化是由油门和舵机决定的),运动控制模块主要是接收航迹跟踪和避障算法输出的期望速度和航向,计算处理后产生油门和舵角值.油门值可由式(4)反推得到,即而舵角值的确定则需利用GPC-PID串级控制方法[23],将广义预测理论与PID控制器相结合,来克服无人水面艇系统中存在的时滞,故该系统具有较强的鲁棒性和抗干扰能力.GPC-PID运动控制原理如图8所示.图8中,ψr为航迹跟踪和自主避障模块输出的期望航向角,ψ为实际航向角,ψr和ψ作为GPC的输入量输出控制量δr,而δr作为PID输入量输出舵角δ,G1为2.2.2节提到的一阶Nomoto模型,其输出为航向角速度γ.2.6 航向干扰无人水面艇在航行过程中会受到风、浪、流的干扰,其中风浪会对无人水面艇的艏摇、横摇和纵摇产生影响,故本工作加入2阶风浪模型对艏摇产生干扰[24]:式中,ω(s)为高斯白噪声,h(s)为2阶海浪传递函数,其中λ为阻尼系数,ω0为海浪频率,σ为海浪强度.在完成仿真系统中的流程和方案设计后,为了验证仿真系统设计的正确性和合理性,以航线跟踪精度和避障效果为指标对仿真系统进行实验.本仿真实验为半实物仿真,硬件平台有显控、仿真和控制计算机等.3.1 算法步骤利用显控计算机进行全局路径规划,由仿真计算机更新位置、仿真障碍信息等,由控制计算机根据任务航线、障碍信息、输出期望油门和舵角给仿真计算机,以便进行位置和姿态更新.仿真实验算法步骤如下:(1)设定全局路径L、障碍物Φ,初始化无人水面艇位置和姿态,如位置P(x,y)、航向ψ、速度V等;(2)激光检测Φ,输出激光数据,经过障碍处理后得到障碍物Φ′;(3)导航和避障根据L,Φ′,P(x,y),ψ,计算输出期望航向ψr和速度Vr.(4)将ψr,Vr作为运动控制的输入量,输出舵角δ和油门E;(5)将δ,E作为位置和姿态更新的输入,输出更新后的位置和姿态信息P(x,y),ψ,V,航向角速度γ;(6)循环执行步骤(2)~(5),直到到达路径结束点.3.2 仿真对象本仿真实验以上海大学“精海3号”无人水面艇的控制系统为仿真对象,艇体和激光参数如表1和2所示.经测试,“精海3号”无人水面艇静态航向角速度增益K=0.3176,时间常数T=0.2417,舵角与航向角速度关系为根据Fossen的理论[25],通常情况下海浪干扰参数Kω=4693.9,λ=0.12,ω0=0.8,则3.3 仿真实验结果和分析仿真实验基于电子海图来实现,图9(a)~(c)依次为无人水面艇任务区域覆盖扫测、动态避障和多障碍避障实验结果.图中黑色同心圆为无人水面艇,红色圆为障碍物,蓝色半圆点云为激光打在障碍物上的点,褐色圆为障碍处理后得到的障碍信息,红色轨迹线为障碍信息的运动轨迹,红色虚线圆为激光检测范围.可以看出,无人水面艇在遇到障碍物后能够按照正确的方向进行避障并回归航线.如图10所示,在导航循迹实验中,取50 s,500帧数据进行分析,其中(b)和(c)为航迹、航向偏差情况.实验得出最大航迹偏差为13.1 m,最大航向偏差为10.7°.结果满足IEC 6206标准[26]航向偏差不超过15°,航迹偏差不超过35 m的航迹控制要求.取“精海3号”在东海2级海况下测绘数据进行分析,得到最大航迹偏差为9.7 m,最大航向偏差7.8°.经分析可知,仿真结果和实际实验结果相近,仿真结果偏差较大的原因是受海浪干扰模型影响较大.本工作为提高海上调试效率,为无人水面艇的控制系统提供了测试和分析平台,搭建了仿真系统;提出了无人水面艇仿真系统设计方案,即将仿真系统分为5个模块,利用电子海图信息和生物激励神经网络的完全遍历路径规划方法实现全局路径规划,并根据无人水面艇在海上航行实际情况,实现了无人水面艇的位置和姿态更新;仿真了激光数据,进行障碍检测处理,并对动态障碍进行运动预测.在进行导航和避障实验中加入了2阶风浪干扰,实验结果和数据分析表明本仿真系统能对无人水面艇控制系统性能进行测试,验证了所设计系统的合理性和正确性.【相关文献】[1]FOssEN T I,BREIVIK M,SKJETNE R.Line-of-sight path following of underactuated marine craft[C]//Proceedings of the 6th IFAC MCMC.2003:244-249.[2]PAVLOV A,NORDAHL H,BREIVIK M.MPC-based optimal path following for underactuated vessels[J].IFAC International Conference on Manoeuvring and Control of Marine Craft,2009, 42(18):340-345.[3]BREIVIK M,HOVsTEIN V E,FOssEN T I.Straight-line target tracking for unmanned surface vehicles[J].Modeling,Identifcation and Control,2008,29(4):131-149.[4]KJERsTAD K,BREIVIK M.Weather optimal positioning control for marine surface vessels[J]. IFAC Conference on Control Applications in Marine Systems,2010,43(20):114-119.[5]BREIVIK M,STRAND J P,FOssEN T I.Guided dynamic positioning for fully actuated marine surface vessels[C]//Proceedings of the 7th IFAC MCMC.2006:1-6.[6]BREIVIK M,FOssEN T I.Guidance laws for planar motion control[J].IEEE Conference on Decision and Control,2008,16(5):570-577.[7]BREIVIK M.Nonlinear maneuvering control of underactuatedships[D].Trondheim:Norwegian University of Science and Technolgy,2003:75-79.[8]LARRAZABAL J M,PE˜NAs M S.Intelligent rudder control of an unmanned surface vessel[J]. Expert Systems with Applications,2016,55:106-117.[9]杨树仁.船舶航迹控制算法的研究与仿真[D].长春:吉林大学,2014.[10]卢艳爽.水面无人艇路径规划算法研究[D].哈尔滨:哈尔滨工程大学,2010.[11]杨学斌.无人艇体系结构研究[D].大连:大连海事大学,2011.[12]朱博,邓三鹏,王英飞,等.基于生物激励神经网络的移动机器人遍历路径规划[J].装备制造技术,2014(12):30-32.[13]YOON C.Human following of a mobile service robot with single laser rangefnder[J].Journal of Measurement Science and Instrumentation,2012,3(1):94-98.[14]Xu F,XIAO T,XING X W.Identifcation of Nomoto models with integral sample structure for identifcation[C]//Chinese Control Conference.2014:6721-6725.[15]谌彤童.三维激光雷达在自主车环境感知中的应用研究[D].长沙:国防科学技术大学,2011.[16]HIMMELsBACH M,M¨uLLER A,L¨uTTEL T,et al.LIDAR-based 3D object perception[C]// Proceedings of 1st International Workshop on Cognition for Technical Systems.2008. [17]APPIAH K,HuNTER A,DICKINsON P,et al.A run-length based connected component algorithm for FPGA implementation[C]//International conference on Field-Programmable Technology. 2008:177-184.[18]JOHNsTON C T,BAILEY D G.FPGA implementation of a single pass connected components algorithm[C]//4th IEEE International Symposium on Electronic Design,Test and Applications. 2008:228-231.[19]LINGRAs P,WEsT C.Interval set clustering of web users with rough K-means[J].Journal of Intelligent Information Systems,2004,23(1):5-16.[20]JABBARIAN M,BIZAKI H K.Target tracking in pulse-doppler MIMO radar by extended Kalman flter using velocity vector[J].20th Iranian Conference on Electrical Engineering.2012:1373-1378.[21]CHAKRAVARTHY A,GHOsE D.Collision cones for quadric surfaces[J].IEEE Transactions on Robotics,2011,27(6):1159-1166.[22]FIORINI P.Motion planning in dynamic environments using velocityobstacles[J].International Journal of Robotics Research,1998,17(7):760-772.[23]彭艳,吴伟清,刘梅,等.无人艇航迹跟踪GPC-PID串级控制[J].控制工程,2014,21(2):245-248.[24]THOMAs S.Guidance system for autonomous surfacevehicles[D].Trondheim:Norwegian University of Science and Technology,2015.[25]FOssEN T I.Handbook of marine craft hydrodynamics and motion control[J].IEEE Control Systems,2016,36(1):78-79.[26]刘廷伟.船舶航迹控制及气象导航方法的研究[D].长春:吉林大学,2015.。
现代电子技术Modern Electronics Technique2024年5月1日第47卷第9期May 2024Vol. 47 No. 90 引 言无人潜航器(Unmanned Underwater Vehicle, UUV )在20世纪50年代开始出现,当时主要应用于民用领域,比如:海洋科学领域的考察、水下打捞和水下电缆的铺设等[1]。
随着科学技术的发展,计算机技术、通信技术、大规模集成电路、材料、导航技术[2⁃3]及能源[4]等技术的跨越式发展,无人潜航器的功能更加丰富、下潜深度更深,应用领域也从民用扩展到了军用。
当前人工智能的飞速发展,无人潜航器也从遥控式发展到了自主式[5],其应用也更加广泛。
无人潜航器具有体积小、重量轻、灵活性高等特点,非常适用于水下目标的探测和识别,也被广泛应用于军事用途,成为世界各国研制的重要水下装备。
目前国外有美国、俄罗斯、英国、德国、法国、印度、韩国、日本等十几个国家二十几个研究机构进行无人潜航器研制,足以看出在未来海洋战争中无人潜航器将发挥重要的作用,研究无人潜航器的军事应用有着重要的意义[6⁃8]。
1 无人潜航器在水下作战的应用1.1 无人潜航器在水下侦察搜索中的应用随着传感器技术、能源技术及人工智能技术的发展,无人潜航器探测能力更强、续航距离更远、自主化程度更高,能够在水下作战中发挥更大的作用。
因无人潜航器自身特有的体积小、灵活性高、隐蔽性强的特点,无人潜航器能够完成水下的侦察与搜索任务,可以对重要的海域或者重要的海上运输线完成搜索和侦察,也可组成无人潜航器集群对重要海域的海洋水文环境资料进行搜集,并完成搜索数据的回传,以便战时能够提前掌握海域环境[9]。
无人潜航器解算目标运动要素的一种方法张 严, 朱伟良, 程 菲(海军潜艇学院, 山东 青岛 266199)摘 要: 随着计算机技术、人工智能、能源技术的发展,无人潜航器也被应用于军事领域。
在分析了无人潜航器水下侦察搜索、反水雷战、水下诱骗战应用的基础上,提出潜艇采用线导方式控制双无人潜航器进行目标运动要素解算的方法,并分别通过图解法和计算法分析了其解算的原理。
舰船科学技术SHIP SCIENCE AND TECHNOLOGYVol. 44, No. 7Apr.,
2022
第44卷第
7期
2022年4月
潜水器动态收放缆索过程的动力学仿真研究郑鹏,王小庆,
陈伟
(中国船舶科学研究中心水动力学国家重点实验室,江苏无锡
214082)
摘 要:调研了解到部分上述高端装备或其他水下平台在水下作业过程中,可能遇到的载荷吊放、
载荷拖
曳、载荷回收等工况下缆索张力动态非线性变化的问题。基于Ablow提岀的缆索偏微分控制方程,通过在采用有限 差分法求解过程中动态调整缆索微元的长度进行缆索收放过程仿真,建立了潜水器-缆索-载荷三者间相互耦合的动
力学模型,并基于建立的仿真方法对潜水器拖曳载荷进行变深度航行、潜水器变深度直航过程中起吊载荷过程中缆 索-载荷系统对潜水器的影响及缆索自身的动态响应特性进行了研究。结果表明:稳定拖曳过程中,由于缆索-载荷 系统的阻力,导致速度降低;拖曳起吊过程中的缆索张力变化明显,且缆索张力对收缆的速度变化敏感;拖曳稳定 航行段潜水器速度与理论值一致,表明了本文方法的正确性。
关键词:潜水器;缆索张力;拖曳回收;多体系统;变长度单元
中图分类号:
U674.941 文献标识码
:
A
文章编号:1672 - 7649(2022)07 -0055 -07 doi: 10.3404/j.issn.l672 - 7649.2022.07.011
Research on coupled dynamics simulation of towing and lifting mission load by
submersible
ZHENG Peng, WANG Xiao-qing, CHEN Wei(China Ship Scientific
Research Center, State Key Laboratory of hydrodynamics, Wuxi 214082, China)
无人船编队控制仿真开发无人船编队控制仿真开发随着科技的不断发展,无人船的应用越来越广泛,如水下勘察、海上运输、环境监测、海洋资源开发等。
无人船在实际应用中常常需要组建编队,且编队控制精度和互动问题一直是无人船技术领域关注的热点之一。
在对无人船编队控制算法进行实验之前,需要进行仿真模拟,以评估算法的性能及调试控制参数。
无人船编队控制仿真是通过计算机软件模拟无人船行驶、通信、控制等方面的场景,进行无人船编队控制算法的验证和测试,为实际应用提供数据参考。
无人船编队控制仿真开发包括仿真软件、模型设计、仿真参数调试等方面。
一、仿真软件仿真软件是实现无人船编队控制仿真的核心,主要用来绘制无人船模型、场景、控制器等。
目前,常用的无人船编队控制仿真软件有MATLAB、Simulink、Virtual Marine Technology等。
其中,MATLAB和Simulink是较为常用的仿真软件,支持多种物理模型和仿真场景。
1.MATLABMATLAB是一种高级计算机编程语言及交互式环境。
它拥有广泛的工具箱和可视化工具,方便用户进行各种科学和工程计算、绘图、数据分析和交互式应用程序开发。
MATLAB支持各种编程语言,包括C++、Python和Java。
在无人船编队控制仿真开发中,MATLAB常用于建立无人船物理模型,并通过Simulink流程建模和仿真,实现控制器的仿真算法。
同时,MATLAB也可以利用自带的图形绘制库进行无人船编队的可视化。
2.SimulinkSimulink是MATLAB的附属工具,是一种用于建模和仿真系统的图形化仿真环境。
它与MATLAB集成,提供了多种系统模型的设计,模拟和分析工具。
同时,Simulink的可视化界面也为用户提供了近似实时的仿真分析空间。
在无人船编队控制仿真开发中,Simulink常用于建立无人船编队控制模型,通过简单易用的图形界面进行仿真模拟及控制器设计。
3.Virtual Marine TechnologyVirtual Marine Technology是一款功能强大的船舶仿真软件。
第45卷第6期2023年12月指挥控制与仿真CommandControl&SimulationVol 45㊀No 6Dec 2023文章编号:1673⁃3819(2023)06⁃0112⁃09无人艇运动控制仿真试验训练系统设计研究刘㊀峰1,许二旭2(1.中国人民解放军92941部队43分队,辽宁葫芦岛㊀125000;2.中国人民解放军91976部队,广东广州㊀510080)摘㊀要:针对靶场无人装备仿真试验手段缺乏问题,在构建无人艇动力学模型和运动模型的基础上,基于CoppeliaSim仿真软件,设计了一套无人艇运动控制仿真试验训练系统,构建了PID㊁滑膜和模糊PID三种控制器,并采用Navmesh寻路算法实现无人艇仿真平台的路径规划㊂该系统弥补了靶场无人装备仿真试验的不足,解决了无人艇模拟训练手段缺乏的难题,可为无人装备数字化试验鉴定提供基础支撑㊂关键词:无人艇;运动控制;仿真系统;控制器;路径规划中图分类号:U664 82㊀㊀㊀㊀文献标志码:A㊀㊀㊀㊀DOI:10.3969/j.issn.1673⁃3819.2023.06.017DesignandresearchofsimulationtesttrainingsystemforunmannedboatmotioncontrolLIUFeng1,XUErxu2(1.Navy92941Unit43ofPLA,Huludao125000;2.Unit91976ofPLA,Guangzhou510080,China)Abstract:Aimingatthelackofsimulationtestmeansforunmannedequipmentatthefiringrange,thispaperdesignsasim⁃ulationtesttrainingsystemforunmannedvehiclemotioncontrolbasedonCoppeliaSimsimulationsoftwareandonthebasisoftheconstructionofunmannedvehicledynamicsmodelandmotionmodel,andconstructsthreecontrollers,PID,synoviumandfuzzyPID.TheNavmeshroutingalgorithmisselectedtorealizethepathplanningoftheunmannedboatsimulationplat⁃form.Thesystemfillsthegapinthesimulationtestofunmannedequipmentatthefiringrange,solvestheproblemofthelackofsimulationtrainingmeansofunmannedboat,andcanprovidebasicsupportforthedigitaltestevaluationofunmannede⁃quipment.Keywords:USV;motioncontrol;simulationsystem;controllers;pathplanning收稿日期:2023⁃06⁃12修回日期:2023⁃07⁃07作者简介:刘㊀峰(1983 ),男,工程师,研究方向为导航㊁两栖陆战装备㊂许二旭(1979 ),男,副教授㊂㊀㊀无人装备模型构建和运动控制仿真技术是近年来国内外军工靶场的研究重点,也是无人装备数字化鉴定及模拟训练的核心关键,其鉴定考核具有标准依据少㊁方法设计难等特点,传统试训方法不能直接观察到控制系统各参数变化对运动稳定性的影响,试训人员对控制参数的调试理解不深,使得训练效果不佳,因此,通过模型构建及仿真技术来验证无人装备运动控制参数和性能,对其融入作战体系及能力生成具有重要意义[1]㊂本文基于CoppliaSim仿真软件设计无人艇运动控制仿真试验训练系统,基于构建的无人艇动力学模型和运动学模型,设计构建稳健比例积分(PID)控制㊁滑膜控制和模糊PID控制三种控制器,仿真系统与实装艇通过通信系统连接,试训人员能够通过仿真系统调节参数,结合实船演示和仿真可视化的方法验证不同控制器下的轨迹变化特性,为无人艇试验训练提供手段方法和基础支撑[2]㊂设计过程包括CoppeliaSim中的仿真运动过程和软件界面设计整合过程㊂CoppeliaSim中的仿真过程包括路径规划㊁控制器设计㊁运动学仿真三大部分,仿真全过程的逻辑框图如图1所示㊂在实现运动仿真后,软件界面对系统功能进行整合,以运动模式㊁行程点和控制参数为输入量,经运动学仿真后,将运动状态输出显示并进行运动渲染,其具体逻辑框图如图2所示㊂1㊀无人艇实体模型及环境搭建为仿真试验提供可视化效果,需要构建相应的无人艇实体模型和水面环境模型[3]㊂1 1㊀无人艇实体模型考虑仿真平台运行时的性能和无人艇的运动特性,本系统采用CoppliaSim仿真软件中的简单立方体代替无人艇实体模型进行仿真㊂螺旋桨也用两个简单的立方体代替,并用力传感器将替代螺旋桨的立方体与无人艇体进行连接,虚拟无人艇组成部分从属关系图如图3所示㊂为了画面效果,依据试验用实物艇建立无人艇的实体建模,将模型放置在显示的图层,而将仿真的实体放在隐藏的图层,如图4所示㊂第6期指挥控制与仿真113㊀图1㊀仿真过程逻辑图Fig 1㊀Logicdiagramofsimulationprocess图2㊀软件逻辑框架Fig 2㊀Softwarelogicframework图3㊀虚拟无人艇组成部分从属关系图Fig 3㊀Dependencydiagramofvirtualunmannedsurfacevehiclecomponents1 2㊀环境模型搭建CoppliaSim软件提供的模型不具备如水浮力㊁水阻力等物理特征,而本系统关注的是无人艇的水面运动图4㊀导入完成的全部模型显示Fig 4㊀Showsofimportedmodel状态,即我们是基于三自由度无人艇模型进行后续设计工作,因此,对模型进行简化,将无人艇实体放在一个支撑大平面上进行仿真㊂水阻力以及风浪等干扰可以通过在调用软件中添加力和力矩的函数进行仿真㊂2㊀动力学模型构建2 1㊀无人艇坐标系无人艇在水面上航行时,可将运动行为看成在六个自由度上的运动㊂运动行为分为线性运动和旋转运动㊂线性运动为沿轴x0㊁y0㊁z0方向的运动,旋转运动则是沿这三轴的往复运动㊂为了更加详细地描述无人艇的运动状态和位置,我们建立了描述无人艇位置的惯性坐标系和描述其运动姿态的机体坐标系,如图5所示㊂惯性坐标系也称为大地坐标系,其原点选为地球表面某个运动参考点,通常情况下为无人艇某个时刻的重心㊂o0x0轴的正方向指向正北方,o0y0轴的正方向指向正东方,o0z0轴的正方向垂直向下指向地心㊂机114㊀刘㊀峰,等:无人艇运动控制仿真试验训练系统设计研究第45卷图5㊀坐标系构建Fig 5㊀Constructionofcoordinatesystem体坐标系也称为运动坐标系,用于描述无人艇的运动信息[4]㊂o为无人艇的重心,该坐标系会跟随艇一起进行移动㊂规定:ox轴正方向指向艇首,oy轴正方向指向艇体的右舷,oz轴正方向垂直向下㊂无人艇状态变量定义为η1=[x,y,z],表示无人艇在惯性坐标系的位置,η2=[φ,θ,ψ]表示无人艇在惯性坐标系中的姿态;v1=[u,v,w]为无人艇在运动坐标系中的线速度,分别对应了纵荡㊁横荡和垂荡速度,v2=[p,q,r]则是代表在同一坐标系中的角速度,分别对应了横摇㊁纵摇和艏摇角速度;τ1=[X,Y,Z]代表了无人艇在运动坐标系下所受到的各方向外力,τ2=[K,M,N]则代表在同一坐标系下的各方向力矩㊂综上,定义η=[η1,η2]为无人艇的位置向量,v=[v1,v2]为速度向量,τ=[τ1,τ2]为其所受外力/力矩[5]㊂2 2㊀无人艇动力学模型无人艇线速度在两坐标系之间转换关系为η㊃1=J1(η2)v1(1)其中,J1(η2)为转换矩阵,表示为J1η2()=cosθcosψsinψsinθcosψ-cosφsinψcosφsinθcosψ+sinφsinψcosθsinψsinψsinθsinψ+cosφcosψcosφsinθsinψ-sinφcosψ-sinθsinφcosθcosφcosθéëêêêùûúúú(2)㊀㊀无人艇角速度在两坐标系之间转换关系为η㊃2=J2(η2)v2(3)其中,J2(η2)为转换矩阵,表示为J2η2()=1sinφtanθcosφtanθ0cosφ-sinθ0sinφsecθcosφsecθéëêêêùûúúú(4)结合式(1)和(3)可得无人艇水面运动模型:η㊃1η㊃2éëêêùûúú=J1η2()03ˑ303ˑ3J2η2()éëêêùûúúv1v2éëêêùûúú(5)式(5)可简化为η㊃=Jη()v(6)将无人艇视为刚体,由牛顿-欧拉方程推导出的刚体力和力矩平衡方程为m[v㊃1+v2v1+v㊃2rg+v㊃2(v㊃2rg)]=τ1I0v㊃2+v2I0v2+mrgv㊃1+v2v1()=τ2{(7)其中,rg=[xg,yg,zg]为重心坐标,I0为刚体转动惯量㊂将式(7)改写为MRBv㊃+CRBv()v=τRB(8)式中,v=u,v,w,p,q,r[],τRB=[X,Y,Z,K,M,N]㊂MRB为刚体惯性矩阵:MRB=m000m000m0mzg-myg-mzg0mxgmyg-mxg00-mzgmygmzg0-mxg-mygmxg0Ix-Ixy-Ixz-IyxIy-Iyz-Izx-IzyIzéëêêêêêêêêêùûúúúúúúúúú(9)CRBv()是刚体科里奥利向心力矩阵,令CRBv()=[C1,C2],则:C1=000000000-m(ygq+zgr)m(ygp+w)m(zgp-v)m(xgq-w)-m(zgr+xgp)m(zgq+u)m(xgr+v)m(ygr-u)-m(xgp+ygq)éëêêêêêêêêùûúúúúúúúú第6期指挥控制与仿真115㊀C2=m(ygq+zgr)-m(xgp-w)-m(xgp+v)-m(ygp+w)m(zgr+xgp)-m(ygp-u)-m(zgp-v)-m(zgp+u)m(xgp+ygq)0-Iyzq-Ixzp+IzzrIyzr+Ixyp-IyyqIyzq+Ixzp+Izzr0-Ixzr-Ixyq+Ixxp-Iyzr-Ixyp+IyyqIxzr+Ixyq-Ixxp0éëêêêêêêêêêùûúúúúúúúúú(10)㊀㊀向量τRB由流体力和力矩向量τH㊁外部环境干扰力和力矩向量τE㊁推进力和力矩向量τ组成㊂下面重点介绍τH的组成㊂由流体力学可知,作用于刚体上的流体动力和力矩是线性叠加的[6]㊂流体动力和力矩由附加质量㊁水动力阻尼和重力与浮力引起的回复力三部分组成,因此,τH可以写成τH=-MAv㊃-CAv()v-Dv()v-g(η)(11)MA是附加质量矩阵㊂MA=-X㊃uX㊃vX㊃wX㊃pX㊃qX㊃rY㊃uY㊃vY㊃wY㊃pY㊃qY㊃rZ㊃uZ㊃vZ㊃wZ㊃pZ㊃qZ㊃rK㊃uK㊃vK㊃wK㊃pK㊃qKrM㊃uM㊃vM㊃wM㊃pM㊃qM㊃rN㊃uN㊃vN㊃wN㊃pN㊃qN㊃réëêêêêêêêêêùûúúúúúúúúú(12)CAv()是流体力学科里奥利向心力矩阵㊂CAv()=0000-a3a2000a30-a1000-a2a10a30-a1b30-b1a30-a1b30-b1-a2a10-b2b0éëêêêêêêêêêùûúúúúúúúúú(13)其中,a1a2a3b1b2b3éëêêêêêêêêêùûúúúúúúúúú=X㊃uX㊃vXwX㊃pX㊃qX㊃tY㊃uY㊃vYwY㊃pY㊃qY㊃㊃rZ㊃uZvZ㊃wZ㊃pZ㊃qZrK㊃uK㊃vKwK㊃pK㊃qK㊃rM㊃uM㊃vM㊃wM㊃pM㊃qM㊃rN㊃uN㊃vN㊃wN㊃pN㊃qN㊃réëêêêêêêêêêùûúúúúúúúúúuvwpqréëêêêêêêêêùûúúúúúúúú(14)Dv()是阻尼矩阵,可写为线性和非线性阻尼矩阵的和㊂Dv()=D+Dn(v)(15)其中,线性阻尼矩阵D为D=XuXvXwXpXqXrYuYvYwYpYqYrZuZvZwZpZqZrKuKvKwKpKqKrMuMvMwMpMqMrNuNvNwNpNqNréëêêêêêêêêêùûúúúúúúúúú(16)g(η)是回复力和力矩矩阵向量㊂gη()=W-B()sinθ-W-B()cosθsinφ-W-B()cosθcosφ-ygW-ybB()cosθsinφ+zgW-zbB()cosθsinφzgW-zbB()sinθ+xgW-xbB()cosθcosφ-xgW-xbB()cosθsinφ-ygW-ybB()sinθéëêêêêêêêêùûúúúúúúúú(17)其中,W为重力,B为浮力㊂W=mg,B=ρgρᶄ(18)式中,m为无人艇质量,g为重力加速度,ρ为水密度,ρᶄ为排水密度㊂2 3㊀三自由度无人艇运动模型在无人艇的实际航行过程中,主要关注运动坐标系纵向运动㊁横向运动以及转艏运动(水平面上的运动),因此,可以将六自由度无人艇运动数学模型简化为三自由度的平面运动数学模型,使运动控制问题得到简化,水平坐标系如图6所示㊂图6㊀水平坐标系Fig 6㊀Horizontalcoordinatesystem假设无人艇的几何中心与重心重合,同时忽略风㊁浪和流等外界干扰因素,则可以得到无人艇水面三自由度运动数学模型:116㊀刘㊀峰,等:无人艇运动控制仿真试验训练系统设计研究第45卷η㊃=J(η)vMv㊃+C(v)v+D(v)v=τ{(19)其中,位置矢量η=[x,y,ψ],速度矢量v=[u,v,r],而仿真对象是双桨无人艇,则其所受推力和力矩矢量τ=[T1+T2,0,(T1-T2)B/2]㊂T1和T2分别是无人艇左右桨提供的推力,B是桨间的水平距离㊂对于螺旋桨提供的推力,用以下数学模型进行描述:T=(1-tP)ρn2D4PKT(JP)(20)式中,ρ表示水密度;JP()为桨推力系数,其中,JP为进速系数;tP为伴流系数;DP表示螺旋桨直径;n是螺旋桨转速㊂转换矩阵J(η)可表示为Jη()=cosψ()-sinψ()0sinψ()cosψ()0001éëêêêùûúúú(21)惯性矩阵M=MRB+MA,科里奥利力矩阵C(v)=CRB(v)+CA(v),阻尼矩阵D(v)=D+Dn(v)㊂其中:MRB=m0-myg0mmxg-mygmxgIzéëêêêêùûúúúú(22)MA=-X㊃u000-Y㊃v-Y㊃r0-N㊃v-N㊃réëêêêêùûúúúú(23)CRBv()=00-mxgr+v()00mumxgr+v()-mu0éëêêêêùûúúúú(24)CAv()=00Y㊃vv+12Y㊃r+N㊃v()r00-X㊃uu-Y㊃vv-12Y㊃r+N㊃v()rX㊃uu0éëêêêêêêùûúúúúúú(25)D=-Xu000-Yv-Yr0-Nv-Nréëêêêêùûúúúú(26)Dnv()=-Xuuu000-Yvvv000-Nrrréëêêêêùûúúúú(27)前文中假设无人艇重心与几何中心重合,因此,(xg,yg)=(0,0)㊂在不考虑运动坐标系原点变化的情况下,假设无人艇结构前后对称时,M㊁D(v)均采取对角线型结构,这表示Y㊃r=0,Yr=0,N㊃v=0,Nv=0㊂结合式(22)和(23)可得M的具体表达形式:M=m-X㊃u000m-Y㊃v000Iz-N㊃réëêêêêùûúúúú=m11000m22000m33éëêêêêùûúúúú(28)结合式(24)和(25)可得C(v)的具体表达形式:Cv()=00-m22v00m11um22v-m11u0éëêêêêùûúúúú(29)阻尼矩阵包含线性阻尼矩阵D以及非线性阻尼矩阵Dn(v),忽略高于二阶的流体动力阻尼项,则非线性阻尼矩阵中各项系数均为零,因此:Dv()=-Xu000-Yv000-Nréëêêêêùûúúúú(30)综上,将式(19)展开为u㊃=m22m11vr-d11m11u+1m11(T1+T2)v㊃=-m11m22ur-d22m22vr㊃=m11-m22m33uv-d33m33r+1m33(T1-T2)B/2x㊃=ucosψ-vsinψy㊃=usinψ+vcosψψ㊃=r(31)无人艇基本参数及无人艇水动力模型参数如表1㊁表2所示㊂表1㊀无人艇基本参数Tab 1㊀Basicparametersofunmannedsurfacevehicle艇长1 25m推进器到纵向轴距15cm艇宽0 61m推进器到船尾距16cm艇高0 21m重心到船首距90cm质量33 25kg最大线速度2 9ʃ0 2m/s桨叶直径9cm最大角速度ʃ2 0ʃ0 1rad/s㊀表2㊀无人艇水动力模型参数Tab 2㊀HydrodynamicmodelparametersofunmannedsurfacevehicleX㊃u-6 86Nr-1 26Y㊃v-17 5Xuu-12 91N㊃r-1 2Yvv-24 2Xu3 02Nrr0Yv-20 5Iz2 70第6期指挥控制与仿真117㊀3㊀控制器构建无人水面艇在海上航行,由于受到风㊁浪㊁流等环境的影响,给定的航向会不可避免地产生偏移㊂精准的航向控制是解决航迹跟踪㊁自主航行㊁自动避碰等问题的前提㊂相对民船而言,无人艇的重量小,速度快,动力有限,风㊁浪㊁流等海洋环境的干扰显得非常突出,其控制要求更高㊂因此,必须研究出更加先进的无人艇运动控制技术来代替传统的控制方法,以满足控制要求[7]㊂无人艇的运动控制器由航向控制器㊁航迹控制器和速度控制器三个控制器组成,为了让试训人员更直观地体验控制器的设计,我们分别采用了PID控制器㊁滑模控制器和模糊PID控制器来进行无人艇的仿真航向控制,并预留可供修改的控制器参数来体验控制器的特性,控制逻辑如图7所示㊂图7㊀无人艇控制逻辑Fig 7㊀Controllogicofunmannedboat㊀㊀我们使用的无人艇是一种双推进器无人艇㊂对于该艇,航迹控制器的输入为无人艇与目标航迹的距离偏差(e),输出为预期航向角(ψd);航向控制器的输入为航向误差(Δψ),输出为推进器转速差值;速度控制器输入为速度误差(ΔV),输出为推进器转速均值㊂航向控制器起到了调整航向的作用,在目标航向已知的情况下向目标航向靠拢;航迹控制器在航向控制器表现较好时才能发挥更好的作用,其会根据航行实际位置与规划航迹综合分析,对无人艇航行航向进行直接决策,扮演指令下达者的角色;速度控制器用来控制无人水面艇航速,在不同水面情况下会有不同航速要求,速度控制器负责达到预期的航速目标㊂设计三种不同的控制器,分别是PID控制器㊁滑模控制器和鲁棒控制器(模糊PID控制器)㊂三种控制器各有优劣㊂其中,PID控制器和模糊PID控制器在存在未知干扰的情况下会有更好的控制效果,具有更高的鲁棒性能㊂而滑模控制器对模型精准度要求更高,但在模型确定的情况下控制精度更好㊂下面是我们用三个控制器进行航向角控制,结果如图8 图10所示㊂由仿真结果可知,三个控制器都可以达到较好的结果,三个控制器对航向角阶跃信号的响应时间和控制精度分别为:PID控制器响应时间:21 65s,收敛后的波动范围图8㊀PID航向角控制结果Fig 8㊀PIDheadinganglecontrolresults图9㊀滑模航向角控制结果Fig 9㊀Slidingmodeheadinganglecontrolresults是(179 0413,180 4245);滑膜控制器响应时间:26 35s,收敛后的波动范围是(169 5212,170 7525);模糊PID控制器响应时间是21 55s,收敛后的波118㊀刘㊀峰,等:无人艇运动控制仿真试验训练系统设计研究第45卷图10㊀模糊PID航向角控制结果Fig 10㊀FuzzyPIDheadinganglecontrolresults动范围是(175 1142,175 6587)㊂三个控制器中,PID控制器最为简单,具有一定的鲁棒性;模糊PID鲁棒性最高,有自适应能力;滑膜控制器对模型的精度要求最高,参数调整较为困难,若相关参数调整不合理,容易出现抖震现象㊂4㊀Navmesh寻路算法Navmesh(NavigationMesh)寻路算法是目前比较主流的一种寻路算法,是一种基于凸多边形网格的寻路方式[8]㊂Navmesh是对基于导航网格寻路体系的统称,实现方式复杂多样㊂但无论具体使用怎样的方法,Navmesh寻路算法的全部流程至少包含两个部分:导航网格构建和寻路算法㊂本平台使用的Navmesh寻路算法原理如图11所示㊂图11㊀Navmesh寻路算法实现原理Fig 11㊀Navmeshpathfindingalgorithmimplementationprinciple4 1㊀导航网格构建导航网格构建(NavigationMeshConstruction)是指根据得到的地图信息中障碍物和陆地的位置,在地图上静态地构建可供寻路算法使用的导航网格的过程㊂导航网格构建的基本流程如下:体素化场景:将提供的源地图数据通过体素化方式建立实体Span,与此同时,对体素化模型进行简单过滤,以去除无人艇可能无法到达的位置;创建Regions:将第一步生成的体素模型描述的可行走区域划分为重叠的2D区域,检测可跨越的Span,并为其生成简单的轮廓数据;创建Polygon:根据轮廓数据对区域进行化简,获得导航多边形,再对得到的导航多边形进行处理,得到凸多边形㊂4 2㊀寻路算法寻路算法(PathFindingAlgorithm)是指根据前面构建的导航网格体,迅速生成一条可供无人艇在自动驾驶模式下无碰撞行驶的合理路径的过程㊂寻路算法的实现主要分为以下三个阶段㊂4 2 1㊀确定优化目标为构建的凸多边形网格体添加一个 成本 值,这样在搜索路径时,寻路算法会尝试找到总成本最低的路径㊂4 2 2㊀确定凸多边形集合以经过导航网格的总 成本 值最小为目标,使用其他寻路算法可以获得无人艇达到终点所需要经过的导航网格集合㊂这里所使用的寻路算法通常为A∗算法[9]㊂A∗算法是一种应用十分广泛的寻路算法,不一定能找到最短的路径,但可以在较短的时间内搜索到一条相对较短的路径,在Dijkstra算法基础上引入启发因子,寻找那些 看起来离终点更近 的节点㊂A∗算法的运行结果如图12所示㊂图12㊀A∗算法的运行结果∗注:绿色网格为起点,红色网格为终点,灰色网格为障碍物,黄色折线为搜索到的最短路径㊂Fig 12㊀OperationresultsoftheA∗algorithm4 3㊀算法仿真应用通过使用Navmesh寻路算法,系统可以生成无人艇到达终点的合理路径㊂本文设计的路径规划方法在仿真平台中实现了以下目标:根据较复杂的场景中的障碍物信息,构建导航网格体,并在短时间内迅速生成一条成本相对较低的路径;将得到的路径转化为保障无人艇不接触障碍的路径点,同时有利于提高控制器的控制效果;根据路径的整体情况以及船的转弯㊁制动㊁加速性能得到速度矩阵,并根据路径点的分布㊁位置㊁距离等参数进行速度矩阵的再次修正,最终传输给控制器进行每一段速度的控制,如图13所示㊂第6期指挥控制与仿真119㊀图13㊀路径的初步规划Fig 13㊀Preliminaryplanningofthepath4 4㊀路径生成优化4 4 1㊀不同速度下的寻路优化在实际运行过程中,随着速度的增大,转弯半径也会增大,进而影响无人艇的路径规划[10]㊂为了能够在仿真平台中体现出这一特点,我们为障碍物的碰撞模型增添了一段安全距离,以确保有一定宽度的无人艇可以安全通过㊂设置安全距离后的效果如图14所示㊂图中,绿色方格是构建的导航网格体,而不规则的灰色环围成的范围则是设置的安全距离㊂由于转弯半径与速度有关,随着速度的增大,其安全距离也应随之增大,在寻路算法运行的时候,会把安全距离以内的区域认定为被阻挡的区域,进而保证船在高速状态下也可选择合适且安全的路径㊂图14㊀障碍物实际碰撞范围影响导航网格体Fig 14㊀Actualcollisionrangeofobstaclesaffectsthenavigationgridbody图15㊁16是在搭建合适的场景后,无人艇行驶速度慢与快对路径生成带来的影响㊂在慢速情况下,无人艇的转弯半径较小,因此,寻路算法可以在障碍物比较密集的地方生成路径;而在快速情况下,无人艇的转弯半径增大,为了行驶的安全以及保持较快的行驶速度,此时路径轨迹的生成则会选择在宽阔的海面上㊂4 4 2㊀减少船的转弯在实际运行过程中,频繁的转弯会影响无人艇的运行速度,同时可能会影响控制器的效果㊂从图16可以看出,使用A∗算法得到的路径效果并不是真正意义上的最短路径,这是由A∗算法所使用的启发式函数决定的,而之后使用的路径平滑算法也无法减少转弯的图15㊀较低速度下的路径规划Fig 15㊀Pathplanningatlowerspeeds图16㊀较高速度下的路径规划Fig 16㊀Pathplanningathigherspeeds次数,如图17所示㊂因此,需要对生成的路径点进行额外的处理㊂下面将对处理的方式进行介绍:图17㊀仿真平台中减少转弯次数的算法示意图Fig 17㊀SchematicdiagramofthealgorithmforreducingthenumberofturnsinthesimulationplatformStep1:对于寻路算法得到的点集,如果点集中的元素个数不少于3个,则转入Step2处理;否则不需要处理;Step2:取出点集中的第一个点与第三个点,记为C1和C3,由C1向C3发射一条射线,如果射线没有命中障碍物,说明无人艇可以通过直线直接到达C3,可以把二者间的点删除,再取出点集中的第三个点进行检测;如果命中障碍物,说明C1和C3间的点需要保留,此时取出C2和C4进行检测;Step3:重复上述过程,直到取出点集中的倒数第二个点为止㊂在减少无人艇转弯次数时,如果删除过多点,可能会导致无人艇在某次转弯过程中转弯幅度过大,会对行驶造成不好的影响㊂因此,我们基于一定规则,对本该删除的点进行必要的取舍,以保证在不降低无人艇120㊀刘㊀峰,等:无人艇运动控制仿真试验训练系统设计研究第45卷转弯平稳性的同时,减少无人艇转弯次数㊂最后达到的路径规划效果如图18所示㊂相对图16,无人艇自动生成的路径中转弯次数出现了较明显减少㊂图18㊀使用优化方法减少路径中的转弯次数Fig 18㊀Usingoptimizationmethodstoreducethenumberofturnsinthepath5㊀结束语本文设计的仿真试验训练系统,弥补了靶场无人艇仿真试验鉴定领域不足,基本解决了无人艇模拟训练手段缺乏的难题㊂基于该系统开展试验训练工作,可加深对无人装备运动控制理论知识的理解掌握,有利于试训人员对各种控制算法及数字化试验鉴定方法的理解㊂同时,基于该系统,能够对无人装备运动控制试验训练方法进行创新设计,有助于提升海军靶场试验训练综合保障能力㊂参考文献:[1]㊀刘星,杨冰,霍清华,等.无人艇工程控制课程实践教学平台开发[J].造船技术,2023,51(1):35⁃39,63.LIUX,YANGB,HUOQH,etal.Developmentofpracticeteachingplatformofunmannedboatengineeringcontrolcourse[J].MarineTechnology,2023,51(1):35⁃39,63.[2]㊀迟刚,胡晓峰,吴琳.异构模型系统协同仿真与联合运行研究[J].系统仿真学报,2014,26(11):2704⁃2708.CHIG,HUXF,WUL.Researchofcoordinatedsimulationandjointactionamongheterogenousmodelsys⁃tems[J].JournalofSystemSimulation,2014,26(11):2704⁃2708.[3]㊀董慧颖,段云波.水面无人艇运动控制系统建模与仿真[J].沈阳理工大学学报,2017,36(1):77⁃84.DONGHY,DUANYB.Modelingandsimulationofmo⁃tioncontrolsystemofunmannedsurfacevessel[J].JournalofShenyangLigongUniversity,2017,36(1):77⁃84.[4]㊀谢钢.GPS原理与接收机设计[M].北京:电子工业出版社,2009:25⁃33.XIEG.PrinciplesofGPSandreceiverdesign[M].Bei⁃jing:PublishingHouseofElectronicsIndustry,2009:25⁃33.[5]㊀丁晓伟.无人水面艇航行控制器的研究与实现[D].大连:大连海事大学,2018:55⁃78.DINGXW.ResearchandimplementationofUSVnaviga⁃tioncontroller[D].Dalian:DalianMaritimeUniversity,2018:55⁃78.[6]㊀SAUNDERSRNCS.HIS.Jane sfightingships[M].Coulsdon:UnitedKingdom,2016:251⁃255.[7]㊀张维明,黄松平,黄金才,等.多域作战及其指挥控制问题探析[J].指挥信息系统与技术,2020,11(1):1⁃6.ZHANGWM,HUANGSP,HUANGJC,etal.Analysisonmulti⁃domainoperationanditscommandandcontrolproblems[J].CommandInformationSystemandTechnology,2020,11(1):1⁃6.[8]㊀李桥.基于NavMesh自动寻路算法在动态寻路问题上的研究[D].武汉:华中科技大学,2019.LIQ.ResearchondynamicroutingproblembasedonNavMeshautomaticroutingalgorithms[D].Wuhan:Hua⁃zhongUniversityofScienceandTechnology,2019.[9]㊀武善平,黄炎焱,陈天德.改进A∗算法的水面舰艇静态航路规划[J].计算机工程与应用,2022,58(23):307⁃315.WUSP,HUANGYY,CHENTD.StaticrouteplanningofsurfaceshipsbasedonimprovedA∗algorithm[J].ComputerEngineeringandApplications,2022,58(23):307⁃315.[10]杜威,丁世飞.多智能体强化学习综述[J].计算机科学,2019,46(8):1⁃8.DUW,DINGSF.Overviewonmulti⁃agentreinforcementlearning[J].ComputerScience,2019,46(8):1⁃8.(责任编辑:张培培)。
总第299期收稿日期:2018年11月14日,修回日期:2018年12月26日作者简介:雷阳,男,研究方向:UUV总体技术。王占勇,男,研究方向:仿真计算。∗
1引言随着近年对深远海领域的可靠研究加大,各类型的深海水下无人潜航器发展迅速。但深海潜航器的总体技术难度巨大,不同于浅水的潜器的研究开发成熟。其需要的保障及试验花费巨大,同时面临的海洋环境更加恶劣。在深海潜航器总体设计之处,运用软件对即将加工的潜航器进行各参数的仿真具有非常重要的意义。仿真结果不仅能够验证设计方案的合理性,同时也能指导总体方案的优化设计,是潜航体总体设计之中不可或缺的一环。本文针对一款深海潜航器进行了下潜轨迹的仿真,并结合半实物模型在湖上摸底试验采集的数据进行对比,对于深海潜航器的设计具有非常重要的指导意义[1~2]。2仿真分析针对给定潜航体进行了自由下沉轨迹仿真,潜航体运动采用基于fortran语言的六自由度运动方程程序模拟,水动力通过基于SIMPLE算法的CFD工具获取六分量插值样本,针对不同压载的下潜轨迹进行了仿真,结果表明,在20kg和30kg压载下均能得到较好的下潜轨迹[3~8]。1)计算模型及坐标系定义
潜航体几何坐标系OXYZ:坐标原点O为潜航体底部圆心,OX轴位于纵向对称平面内,与体轴线重合,指向前方为正;OZ轴位于纵向对称平面内,垂直于OX轴指向下方为正;OY轴由右手法则确定。表1潜航体几何参数序号123456参数项目X向质心坐标Y向质心坐标Z向质心坐标X向浮心坐标Y向浮心坐标Z向浮心坐标符号XcgYcgZcgXcfYcfZcf单位mmmmmm量值等效质心坐标0
0.0253.63400
无人潜航器无动力下潜轨迹仿真分析∗雷阳王占勇(中船重工第七一○研究所宜昌443000)
摘要论文主要研究无人潜航器无动力下潜轨迹仿真分析,采用旋转导数的理论对无动力下潜进行了六自由度的下潜轨迹分析。论文以两种下潜压载为力,分别分析了不同的下潜负浮力工况下,无动力下潜的轨迹及动力参数。关键词轨迹仿真;六自由度;水动力中图分类号N945.13DOI:10.3969/j.issn.1672-9730.2019.05.033
SimulationandAnalysisoftheTrajectoryofUnmannedUnderwaterVehicle
LEIYangWANGZhanyong(710thInstitute,ChinaShipbuildingIndustryGroupCo.,Ltd.,Yichang443000)
AbstractThispapermainlystudiesthesimulationanalysisoftheunpowereddivingtrajectoryofanunmannedunderwaterve⁃hicle.Thetheoryofrotationalderivativeisusedtoanalyzethedivetrajectoryofsixdegreesoffreedomforunpowereddive.Basedontwokindsofdiveballasts,thepaperanalyzesthetrajectoryanddynamicparametersofunpowereddiveunderdifferentdiveandneg⁃ativebuoyancyconditions.KeyWordstrajectorysimulation,sixdegreesoffreedom,hydrodynamicClassNumberN945.13
舰船电子工程ShipElectronicEngineering总第299期
2019年第5期
Vol.39No.5
1382019年第5期舰船电子工程表1中的Xcg、Ycg、Zcg为潜航体质心在几何坐标系中的坐标,Xcf、Ycf、Zcf为浮心在几何坐标系中的坐标。2)六自由度运动方程[6~7]惯性系内质心平动的运动学方程:x=u ; y=v ; z=w惯性系内质心平动的动力学方程:u=Fx/m ; v=Fy/m ; w=Fz/m弹体系内绕质心转动的运动学方程:ìíîïïϕ=ωx+(ωysinϕ +ωzcosϕ)tanθ θ=ωycosϕ-ωzsinϕ ψ=(ωysinϕ +ωzcosϕ)/cosθ 弹体系内绕质心转动的动力学方程:ìíîïïïïïïïïïïïïJxxωx-Jxyωy-Jxzωz+ωy(-Jzxωx-Jzyωy+Jzzωz)-ωz(-Jyxωx+Jyyωy-Jyzωz)=Mx-Jyxωx+Jyyωy-Jyzωz+ωz(Jxxωx-Jxyωy-Jxzωz)-ωx(-Jzxωx-Jzyωy+Jzzωz)=My-Jzxωx-Jzyωy+Jzzωz+ωx(-Jyxωx+Jyyωy-Jyzωz)-ωy(Jxxωx-Jxyωy-Jxzωz)=Mz表2弧板几何参数序号123参数项目X向质心坐标Y向质心坐标Z向质心坐标符号XcgYcgZcg单位mmm量值5.5900.546采用四步Runge-Kutta或者Adams方法求解以上12个方程,得到潜航体各时刻的质心位置和姿态角。CFD计算得到水动力Fbx、Fby、Fbz为弹体系下的结果,需要通过以下关系式转换到惯性系下:éëêêùûúúFxFyFz=éëêêùûúúcosψcosθ-sinψcosϕ+cosψsinθsinϕcosψsinθcosϕ+sinψsinϕsinψcosθcosψcosϕ+sinψsinθsinϕsinψsinθcosϕ-cosψsinϕ-sinθcosθsinϕcosθcosϕéëêêêêùûúúúúFbxFbyFbz为了获得弹体系上相对于质心的浮力矩Mbx、Mby、Mbz,需要通过以下关系式将惯性系下的浮力转换到弹体系下,然后再计算弹体系下的浮力矩:éëêêêêùûúúúúfbxfbyfbz=éëêêùûúúcosψcosθsinψcosθ-sinθ-sinψcosϕ+cosψsinθsinϕcosψcosϕ+sinψsinθsinϕcosθsinϕcosψsinθcosϕ+sinψsinϕsinψsinθcosϕ-cosψsinϕcosθcosϕéëêêùûúúfxfyfz3)压载为20kg的轨迹仿真[9~13](1)潜航体X向等效质心坐标计算Xcg=(3350kg×3.634m+20kg×5.59m)/(3350kg+20kg)=3.64561m(2)转动惯量计算假设潜航体为均匀质量分布的长L为7m、半径r为0.4m的圆柱体,根据Ixx=mr2/2,Iyy=Izz=mL2/12计算得到:Ixx=Mr2/2=3370´0.42/2=269.6kgm2Iyy=Izz=ML2/12=3370´72/12=13760.8kgm2这里需要说明的是,流体旋转导数的取值规则根据仿真工况的不同有所差别,具体如下。仿真采用如下旋转导数:Cϖyz=-0.6693ϖy+0.0075mϖzz=-0.2485ϖz+0.0094Cϖzy=0.2688ϖz-0.0085mϖyy=-0.4185ϖy+0.0124(3)轨迹仿真压载20kg,下潜6000m,历时16311s,约4.53h,螺旋直径17m,平衡状态下的弹道倾角46.4°,攻角40.6°,侧滑角-9°。如图1~2所示为压载20kg的下潜仿真示意图[8~10]。4)压载为30kg的轨迹仿真(1)潜航体X向等效质心坐标计算Xcg=(3350kg×3.634m+30kg×5.59m)/(3350kg+30kg)=3.65136m(2)转动惯量计算假设潜航体为均匀质量分布的长L为7m、半径r为0.4m的圆柱体,根据Ixx=mr2/2,Iyy=Izz=mL2/12计算得到:Ixx=Mr2/2=3380´0.42/2=270.4kgm2Iyy=Izz=ML2/12=3380´72/12=13801.6kgm2(3)轨迹仿真压载30kg,下潜6000m,历时12671s,约3.52h,螺旋直径17m,平衡状态下的弹道倾角49.5°,攻角42.3°,侧滑角-11°。如图3~4为30kg下潜压载下潜轨迹仿真图。139总第299期xyz图1轨迹图050001000015000t(s)806040200χαβχαβ图2倾角、攻角、侧滑角时间历程xyz图3倾角、攻角、侧滑角时间历程0500010000t(s)806040200χαβχαβ图4倾角、攻角、侧滑角时间历程3结语采用六自由度运动方程对不同压载的下潜轨迹进行了仿真,水动力通过基于Simple算法的CFD数值计算得到,潜航体旋转带来的阻尼力矩假设为某一常数,仿真结果表明:压载20kg,下潜6000m,历时4.53h,螺旋直径17m,平衡状态下的弹道倾角46.4°,攻角40.6°;压载30kg,下潜6000m,历时3.52h,螺旋直径17m,平衡状态下的弹道倾角49.5°,攻角42.3°。该结果与用理论值计算出的水动力参数相差不大。从运动规律、运动参数等方面验证了实际总体设计方案的合理性。参考文献[1]许诚,王作祥,刘宗杰.基于Simulink的潜空导弹运载器水下弹道仿真[J].海军航空工程学院学报,2012(4):435-438.[2]立栋,张宇文,王伟.潜射导弹无动力运载器规避弹道仿真[J].火力与指挥控制,2012,37(7):189-191.[3]田宝国,耿霆.波浪力作用下无动力运载器水弹道仿真[J].计算机与数字工程,2011,39(1):24-27.[4]蔡群,张宇文,张斌,等.大深度无动力运载器弹道散布影响因素分析[J].火力与指挥控制,2013,38(6):72-76.[5]任磊,贾跃,姜选凯.基于Simulink的无动力自导深弹导引弹道仿真研究[J].指挥控制与仿真,2010,32(3):71-75.[6]张志良,翟景春,张雄.基于Simulink的无动力运载器水弹道仿真[J].硅谷,2011(22):20-21.[7]顾闯,逄洪照,张永.海流对潜布水雷水中弹道的影响研究[J].舰船电子工程,2010,30(2):168-171.[8]朱雄勇,张勇,苏广和.基于SIMULINK的某型潜射水雷水下弹道仿真研究[J].舰船电子工程,2011,32(4):53-55.[9]程载斌,刘玉标,刘兆.导弹水下潜射过程的流体——固体耦合仿真[J].兵工学报,2008,29(2):178-183.[10]崔乃刚,曹春泉,韦常柱.潜射导弹水下运动过程仿真分析[J].弹道学报,2009,21(2):95-99.[11]严卫生.鱼雷航行动力学[M].西安:西北工业大学出版社,2005:10-185.[12]黄景全,张宇文.鱼雷流体力学[M].西安:西北工业大学出版社,1989:85-197.[13]AshrafMA.FuzzyTrack-to_trackAssociationandTrackFusionApproachinDistributedMultisensory-mul⁃titargetMultiple-attributeEnvironment[J].Signalpro⁃cessing,2007,87(6):1472-1492.雷阳等:无人潜航器无动力下潜轨迹仿真分析140